Search

Over the recent decades, a number of fast approximate solutions of Lippmann–Schwinger equation, which are more accurate than classic Born and Rytov approximations, were proposed in the field of electromagnetic modeling. Those developments could be naturally extended to acoustic and elastic fields; however, until recently, they were almost unknown in seismology. This paper presents several solutions of this kind applied to acoustic modeling for both lossy and lossless media. We evaluated the numerical merits of those methods and provide an estimation of their numerical complexity. In our numerical realization we use the matrix-free implementation of the corresponding integral operator. We study the accuracy of those approximate solutions and demonstrate, that the quasi-analytical approximation is more accurate, than the Born approximation. Further, we apply the quasi-analytical approximation to the solution of the inverse problem. It is demonstrated that, this approach improves the estimation of the data gradient, comparing to the Born approximation. The developed inversion algorithm is based on the conjugate-gradient type optimization. Numerical model study demonstrates that the quasi-analytical solution significantly reduces computation time of the seismic full-waveform inversion. We also show how the quasi-analytical approximation can be extended to the case of elastic wavefield. [ABSTRACT FROM AUTHOR]

We solve the 2D transverse-electrically polarized domain-integral equation in a layered background medium by applying a Gabor frame as a projection method. This algorithm employs both a spatial and a spectral discretization of the electric field and the contrast current in the direction of the layer extent. In the spectral domain we use a representation on the complex plane that avoids the poles and branchcuts found in the Green function. Because of the special choice of the complex-plane path in the spectral domain and because of the choice to use a Gabor frame to represent functions on this path, fast algorithms based on FFTs are available to transform to and from the spectral domain, yielding an O ( N log ⁡ N ) scaling in computation time. [ABSTRACT FROM AUTHOR]

In this paper we deal with the inverse scattering problem for an impedance obstacle. The aim is to recover both the impedance function and the scatterer simultaneously. Based on boundary integral equations, our method splits the inverse problem into a well-posed direct problem followed by a smaller ill-posed problem which has advantages both in understanding the inverse problem and in the numerical reconstructions. [ABSTRACT FROM AUTHOR]