This review is devoted to some inverse problems arising in the context of linear elasticity, namely the identification of distributions of elastic moduli, model parameters or buried objects such as cracks. These inverse problems are considered mainly for three-dimensional elastic media under equilibrium or dynamical conditions, and also for thin elastic plates. The main goal is to overview some recent results, in an effort to bridge the gap between studies of a mathematical nature and problems defined from engineering practice. Accordingly, emphasis is given to formulations and solution techniques which are well suited to general-purpose numerical methods for solving elasticity problems on complex configurations, in particular the finite element method and the boundary element method. An underlying thread of the discussion is the fact that useful tools for the formulation, analysis and solution of inverse problems arising in linear elasticity, namely the reciprocity gap and the error in constitutive equation, stem from variational and virtual work principles, i.e., fundamental principles governing the mechanics of deformable solid continua. In addition, the virtual work principle is shown to be instrumental for establishing computationally efficient formulae for parameter or geometrical sensitivity, based on the adjoint solution method. Sensitivity formulae are presented for various situations, especially in connection with contact mechanics, cavity and crack shape perturbations, thus enriching the already extensive known repertoire of such results. Finally, the concept of topological derivative and its implementation for the identification of cavities or inclusions are expounded.

/pdf/inverse-problems-in-elasticity-1dpbui74d3.pdf

Inverse problems in elasticity

This paper presents an overview and a detailed description of the key logic steps and mathematical-physics framework behind the development of practical algorithms for seismic exploration derived from the inverse scattering series. There are both significant symmetries and critical subtle differences between the forward scattering series construction and the inverse scattering series processing of seismic events. These similarities and differences help explain the efficiency and effectiveness of different inversion objectives. The inverse series performs all of the tasks associated with inversion using the entire wavefield recorded on the measurement surface as input. However, certain terms in the series act as though only one specific task, and no other task, existed. When isolated, these terms constitute a task-specific subseries. We present both the rationale for seeking and methods of identifying uncoupled task-specific subseries that accomplish: (1) free-surface multiple removal; (2) internal multiple attenuation; (3) imaging primaries at depth; and (4) inverting for earth material properties. A combination of forward series analogues and physical intuition is employed to locate those subseries. We show that the sum of the four task-specific subseries does not correspond to the original inverse series since terms with coupled tasks are never considered or computed. Isolated tasks are accomplished sequentially and, after each is achieved, the problem is restarted as though that isolated task had never existed. This strategy avoids choosing portions of the series, at any stage, that correspond to a combination of tasks, i.e., no terms corresponding to coupled tasks are ever computed. This inversion in stages provides a tremendous practical advantage. The achievement of a task is a form of useful information exploited in the redefined and restarted problem; and the latter represents a critically important step in the logic and overall strategy. The individual subseries are analysed and their strengths, limitations and prerequisites exemplified with analytic, numerical and field data examples.

https://iopscience.iop.org/article/10.1088/0266-5611/19/6/R01/pdf

Inverse scattering series and seismic exploration

Optical tomography has been widely investigated for biomedical imaging applications. In recent years optical tomography has been combined with digital holography and has been employed to produce high-quality images of phase objects such as cells. In this paper we describe a method for imaging 3D phase objects in a tomographic configuration implemented by training an artificial neural network to reproduce the complex amplitude of the experimentally measured scattered light. The network is designed such that the voxel values of the refractive index of the 3D object are the variables that are adapted during the training process. We demonstrate the method experimentally by forming images of the 3D refractive index distribution of Hela cells.

/pdf/learning-approach-to-optical-tomography-571t3nwj7t.pdf

Learning approach to optical tomography

Seismic data regularization, which spatially transforms irregularly sampled acquired data to regularly sampled data, is a long-standing problem in seismic data processing. Data regularization can be implemented using Fourier theory by using a method that estimates the spatial frequency content on an irregularly sampled grid. The data can then be reconstructed on any desired grid. Difficulties arise from the nonorthogonality of the global Fourier basis functions on an irregular grid, which results in the problem of “spectral leakage”: energy from one Fourier coefficient leaks onto others. We investigate the nonorthogonality of the Fourier basis on an irregularly sampled grid and propose a technique called “antileakage Fourier transform” to overcome the spectral leakage. In the antileakage Fourier transform, we first solve for the most energetic Fourier coefficient, assuming that it causes the most severe leakage. To attenuate all aliases and the leakage of this component onto other Fourier coefficients, the data component corresponding to this most energetic Fourier coefficient is subtracted from the original input on the irregular grid. We then use this new input to solve for the next Fourier coefficient, repeating the procedure until all Fourier coefficients are estimated. This procedure is equivalent to “reorthogonalizing” the global Fourier basis on an irregularly sampled grid. We demonstrate the robustness and effectiveness of this technique with successful applications to both synthetic and real data examples.

Antileakage Fourier transform for seismic data regularization

/pdf/marchenko-imagingmarchenko-imaging-2mlkjgjqc1.pdf

Marchenko imagingMarchenko imaging

An Inverse-Scattering Sub-Series for Predicting the Spatial Location of Reflectors without the Precise Reference Medium and Wave Velocity

Flow in a reservoir is controlled predominantly by connectivity of permeability extremes, such as those associated with clear sand channels and shale layers. These elements usually feature complex spatial patterns which are difficult to describe with two-point statistics. Furthermore, specific relationships between the facies are often an important factor in reservoir geology, requiring the use of simulation methods capable of reproducing these associations in order to generate reliable reservoir models. In this work, we were able to bestow physical realism to the geostatistical realizations of a reservoir composed by submarine fans. Multiple-point geostatistics (MPS) relies on training images to model the spatial structure of variables. The MPS simulations were conditioned to seismic and geology data and yielded realistic maps of distributary channels within sand lobes interleaved with shale layers. We concluded that MPS enhances data conditioning and uncertainty assessment with reproduction of specific geometry and facies relationships, making it suitable for geometry-sensitive applications like flow simulations.

Paulo M. Carvalho

Papers

An Inverse-Scattering Sub-Series for Predicting the Spatial Location of Reflectors without the Precise Reference Medium and Wave Velocity

MPS Facies Modelling of a Submarine Fan Reservoir in Southeast Brazil with SNESIM

Automatic variogram model fitting of a variogram map based on the Fourier integral method