Finite-difference contrast source inversion based on hybrid fast conjugate gradient method
(2.State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum (Beijing), Beijing, China 102249)
【Abstract】Finite-difference contrast source inversion (FDCSI) is a solution to inverse scattering. The background model remains unchanged in the inversion, and forward modeling is performed for only one time; this reduces the workload of computation. FDCSI translates the problem of inverse scattering into that of optimization. The cost function is optimized with the conjugate gradient method, which suffers from a low convergence rate and low efficiency. After the study of FDCSI based on acoustic wave equation in frequency domain, we develop a FDCSI algorithm based on a hybrid fast conjugate gradient method to improve the efficiency of inversion. The hybrid fast conjugate gradient method is modified from the fast iterative shrinkage-thresholding algorithm and is feasible for FDCSI. The cost function could converge quickly without more computation in a single iteration; this guarantees fast and robust convergence of FDCSI.
【Keywords】 inverse scattering; contrast source inversion; acoustic wave equation in frequency domain; fast iterative shrinkage-thresholding algorithm; hybrid fast conjugate gradient method;
 Tarantola A. A strategy for nonlinear elastic inversion of seismic reflection data [J]. Geophysics, 1986, 51 (5): 1893–1903.
 Pratt R G, Worthington M H. Inverse theory applied to multisource cross-hole tomography: Part 1, Acoustic wave-equation method [J]. Geophysical Prospecting, 1990, 38 (3): 287–310.
 LI Haishan, YANG Wuyang, YONG Xueshan. Three-dimensional full waveform inversion based on the first-order velocity-stress acoustic wave equation [J]. Oil Geophysical Prospecting, 2018, 53 (4): 730–736 (in Chinese).
 Choi Y, Shin C, Min D J, et al. Efficient calculation of the steepest descent direction for source-independent seismic waveform inversion: an amplitude approach [J]. Journal of Computational Physics, 2005, 208: 455–468.
 Jang U, Min D J, Shin C. Comparison of scaling methods for waveform inversion [J]. Geophysical Prospecting, 2009, 57 (1): 49–59.
 Chi B, Dong L, Liu Y. Full waveform inversion method using envelope objective function without low frequency data [J]. Journal of Applied Geophysics, 2014, 109: 36–46.
 Pratt R G. Seismic waveform inversion in the frequency domain, Part 1: Theory and verification in a physical scale model [J]. Geophysics, 1999, 64 (3): 888–901.
 Kamei R, Pratt R G, Tsuji T. Misfit functions in Laplace–Fourier domain waveform inversion, with application to wide-angle ocean bottom seismograph data [J]. Geophysical Prospecting, 2014, 62 (5): 1054–1074.
 Pratt R G, Shin C, Hick G J. Gauss–Newton and full Newton methods in frequency space seismic waveform inversion [J]. Geophysical Journal International, 1998, 133 (2): 341–362.
 Brossier R, Operto S, Virieux J. Seismic imaging of complex onshore structures by 2D elastic frequency domain full waveform inversion [J]. Geophysics, 2009, 74 (6): WCC105–WCC118.
 Métivier L, Bretaudeau F, Brossier R, et al. Full waveform inversion and the truncated Newton method: quantitative imaging of complex subsurface structures [J]. Geophysical Prospecting, 2014, 62 (6): 1353–1375.
 Kim W K, Min D J. A new parameterization for frequency-domain elastic full waveform inversion for VTI media [J]. Journal of Applied Geophysics, 2014, 109: 88–110.
 van den Berg P M, Ralph E K. A contrast source inversion method [J]. Inverse Problems, 1997, 13 (6): 1607–1620.
 Tarek M H, Michael L O. Simultaneous nonlinear reconstruction of two-dimensional permittivity and conductivity [J]. Radio Science, 1994, 29 (4): 1101–1118.
 van Dongen, Koen W A, William M D. A full vectorial contrast source inversion scheme for 3D acoustic imaging of both compressibility and density profiles [J]. Journal of the Acoustical Society of America, 2007, 121: 1538–1549.
 Wang S D, Wu R S, Liu Y F. The contrast source inversion for reflection seismic data [C]. Extended Abstracts of 78th EAGE Conference and Exhibition, 2016, 01540.
 Abubakar, Hu W, van den Berg P M, et al. A finite-difference contrast source inversion method [J]. Inverse Problems, 2008, 24 (6): 065004–065020.
 Abubakar, Hu W, Tarek M H, et al. Application of the finite-difference contrast-source inversion algorithm to seismic full-waveform data [J]. Geophysics, 2009, 74 (6): WCC47–WCC58.
 Abubakar, Pan G, Li M, et al. Three-dimensional seismic full-waveform inversion using the finite-difference contrast source inversion method [J]. Geophysical Prospecting, 2011, 59 (5): 874–888.
 Han B, He Q, Chen Y, et al. Seismic waveform inversion using the finite-difference contrast source inversion method [J]. Journal of Applied Mathematics, 2014, 109: 1–11.
 He Q, Han B, Chen Y, et al. Application of the finite-difference contrast source inversion method to multi-parameter reconstruction using seismic full-waveform data [J]. Journal of Applied Geophysics, 2016, 124: 4–16.
 DUAN Xiaoliang, ZHAI Hongyu, WANG Yibo, et al. Effect of stratigraphic attenuation on inverse-scattering seismic velocity inversion [J]. Chinese Journal of Geophysics, 2016, 59 (10): 3788–3797 (in Chinese).
 Beck A, Teboulle M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems [J]. SIAM Journal on Imaging Sciences, 2009, 2 (1): 183–202.
 Chambolle A, Dossal C. On the convergence of the iterates of the “fast iterative shrinkage/thresholding algorithm” [J]. Journal of Optimization Theory and Applications, 2015, 166 (3): 968–982.
 CHEN Shaoli, YANG Min. An improved fast iterative shrinkage-thresholding algorithm with variable step size [J]. Computer Technology and Development, 2017, 27 (10): 69–73 (in Chinese).
 ZHANG Guangzhi, SUN Changlu, PAN Xinpeng, et al. Fast conjugate gradient method for frequency domain acoustic wave full waveform inversion [J]. Oil Geophysical Prospecting, 2016, 51 (4): 730–737 (in Chinese).
 Jo C H, Shin C S, Suh J H. An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator [J]. Geophysics, 1996, 61 (2): 529–537.
 DONG Shiqi, HAN Liguo, HU Yong, et al. Acoustic equation high-accuracy modeling in the frequency domain based on MCPML absorbing boundary [J]. Oil Geophysical Prospecting, 2018, 53 (1): 47–54 (in Chinese).
 CHENG Jingwang, LYU Xiaochun, GU Hanming, et al. Full waveform inversion with Cauchy distribution in the frequency domain [J]. Oil Geophysical Prospecting, 2014, 49 (5): 940–945 (in Chinese).
 Polyak B T. The conjugate gradient method in extremal problems [J]. USSR Computational Mathematics & Mathematical Physics, 1969, 9 (4): 94–112.
 Sirgue L, Pratt R G. Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies [J]. Geophysics, 2004, 69 (1): 231–248.