A matrix‐free variant of the distorted Born iterative method for seismic full‐waveform inversion
Kui Xiang,Morten Jakobsen,Kjersti Solberg Eikrem,Geir Nævdal
DOI: https://doi.org/10.1111/1365-2478.13323
2023-01-27
Geophysical Prospecting
Abstract:The distorted Born iterative method reduces a nonlinear inverse scattering problem to a sequence of (ill‐posed) linear inverse scattering problems that can be solved using a regularized least‐squares formulation. This method was originally applied to two‐dimensional electromagnetic problems, but has been implemented to solve acoustic and electromagnetic problems in three dimensions. It has also been applied to seismic problems, but only for moderately large 2D models. Previous applications of the distorted Born iterative method to seismic inverse scattering were based on a matrix representation of the relevant integral operators. The matrix‐based implementation is simple and transparent, but not very suitable for large‐scale computations, since the memory requirement and computational cost scales like and , where N is the number of grid blocks. In this paper, we introduce a matrix‐free variant of the distorted Born iterative method, which is much more suitable for large‐scale problems, since the memory requirements and computational cost has been reduced to N and N log(N), respectively. Our matrix‐free implementation utilizes a fast‐Fourier‐transform‐accelerated iterative method to solve the linear system that results after discretization. In the computation of Fréchet derivatives, we avoid the direct computation of Green's functions for heterogeneous media, by solving an equivalent direct scattering problem. The adjoint of the Fréchet derivative operator is also computed in an efficient matrix‐free manner, by exploiting the physical interpretation and reciprocity of the Green's function. We illustrate the potential of the matrix‐free variant of the distorted Born iterative method using synthetic waveform data for large 2D and 3D models. The main conclusion is that the distorted Born iterative method has been transformed into a more practical tool for seismic as well as electromagnetic and acoustic nonlinear inverse scattering. This article is protected by copyright. All rights reserved
geochemistry & geophysics