Source-receiver Compression Approach For 3D Full-waveform Inversion With an Iterative Forward Solver

We apply a source-receiver compression approach to reduce the computational time and memory usage of threedimensional acoustic full-waveform inversion. By detecting and quantifying the extent of redundancy in the data, we assemble a reduced set of simultaneous sources and receivers that are weighted sums of the physical sources and receivers used in the survey. Because the number of these simultaneous sources and receivers can be significantly less than those of the physical sources and receivers, the computational time and memory usage of any inversion method can be tremendously reduced. The scheme is based on decomposing the data into their principal components using a singular value decomposition scheme and the data reduction is done through the elimination of the small eigenvalues. For the inversion method, we employ the multiplicative-regularized Gauss-Newton approach, in which the regularization parameter is automatically and adaptively determined by the inversion process. For the forward solver, we employ a frequency-domain iterative solver using the perfectly matched layer absorbing boundary condition and a multigrid preconditioner. We demonstrate the advantage of our approach by using the 3D SEG/EAGE salt model.

[1]  Curtis C. Ober,et al.  Faster shot-record depth migrations using phase encoding , 1998 .

[2]  A. Abubakar,et al.  Preconditioned non‐linear conjugate gradient method for frequency domain full‐waveform seismic inversion , 2011 .

[3]  Changsoo Shin,et al.  Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion , 2001 .

[4]  R. Plessix Three-dimensional frequency-domain full-waveform inversion with an iterative solver , 2009 .

[5]  Aria Abubakar,et al.  Application of the finite-difference contrast-source inversion algorithm to seismic full-waveform data , 2009 .

[6]  Denes Vigh,et al.  3D prestack plane-wave, full-waveform inversion , 2008 .

[7]  P. M. Berg,et al.  Contrast Source Inversion Method: State of Art , 2001 .

[8]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[9]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[10]  A. Abubakar,et al.  Simultaneous multifrequency inversion of full-waveform seismic data , 2009 .

[11]  Aria Abubakar,et al.  2.5D forward and inverse modeling for interpreting low-frequency electromagnetic measurements , 2008 .

[12]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[13]  René-Édouard Plessix,et al.  A Helmholtz iterative solver for 3D seismic-imaging problems , 2007 .

[14]  Martin H. Gutknecht,et al.  Variants of BICGSTAB for Matrices with Complex Spectrum , 1993, SIAM J. Sci. Comput..

[15]  J. Bérenger,et al.  Application of the CFS PML to the absorption of evanescent waves in waveguides , 2002, IEEE Microwave and Wireless Components Letters.

[16]  Aria Abubakar,et al.  Full‐waveform seismic inversion using the source‐receiver compression approach , 2010 .

[17]  J. Krebs,et al.  Fast full-wavefield seismic inversion using encoded sources , 2009 .

[18]  P. M. Berg,et al.  Iterative forward and inverse algorithms based on domain integral equations for three-dimensional electric and magnetic objects , 2004 .

[19]  Jean Virieux,et al.  Efficient 3D Frequency-domain Full Waveform Inversion (FWI) with Phase Encoding , 2009 .

[20]  P. M. van den Berg,et al.  Seismic applications of acoustic reciprocity , 1993 .

[21]  A. Abubakar,et al.  A General Framework for Constraint Minimization for the Inversion of Electromagnetic Measurements , 2004 .

[22]  Aria Abubakar,et al.  Towards non‐linear inversion for characterization of time‐lapse phenomena through numerical modelling , 2003 .

[23]  Cornelis Vuik,et al.  A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems , 2005, SIAM J. Sci. Comput..

[24]  Hicks,et al.  Gauss–Newton and full Newton methods in frequency–space seismic waveform inversion , 1998 .

[25]  S. Operto,et al.  3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study , 2007 .

[26]  C. Shin,et al.  An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator , 1996 .

[27]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..