Finite-difference contrast source inversion based on hybrid fast conjugate gradient method
WANG Doudou1,2, WANG Shoudong2, ZOU Shaofeng1, GAO Yanxia1, LIU Han1
1. Sinopec Geophysical Research Institute, Nanjing, Jiangsu 211103, China; 2. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum(Beijing), Beijing 102249, China
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 reduced the workload of computation.FDCSI translates the problem of inverse scattering into a problem of optimization.The cost function is optimized using the conjugate gradient method,which suffers from small convergence rate and low efficiency.After the study of FDCSI using 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 shrin-kage 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.
Tarantola A.A strategy for nonlinear elastic inversion of seismic reflection data[J].Geophysics,1986,51(5):1893-1903.
[2]
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.
[3]
李海山,杨午阳,雍学善.三维一阶速度-应力声波方程全波形反演[J].石油地球物理勘探,2018,53(4):730-736.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.
[4]
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.
[5]
Jang U,Min D J,Shin C.Comparison of scaling me-thods for waveform inversion[J].Geophysical Prospecting,2009,57(1):49-59.
[6]
Chi B,Dong L,Liu Y.Full waveform inversion me-thod using envelope objective function without low frequency data[J].Journal of Applied Geophysics,2014,109:36-46.
[7]
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.
[8]
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.
[9]
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.
[10]
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.
[11]
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.
[12]
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.
[13]
van den Berg P M,Ralph E K.A contrast source inversion method[J].Inverse Problems,1997,13(6):1607-1620.
[14]
Tarek M H,Michael L O.Simultaneous nonli-near reconstruction of two-dimensional permittivity and conductivity[J].Radio Science,1994,29(4):1101-1118.
[15]
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.
[16]
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.
[17]
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.
[18]
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.
[19]
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.
[20]
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.
[21]
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.
[22]
段晓亮,翟鸿宇,王一博,等.地层衰减对地震波速度逆散射反演的影响研究[J].地球物理学报,2016,59(10):3788-3797.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.
[23]
Beck A,Teboulle M.A fast iterative shrinkage-thre-sholding algorithm for linear inverse problems[J].SIAM Journal on Imaging Sciences,2009,2(1):183-202.
[24]
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.
[25]
陈少利,杨敏.改进变步长快速迭代收缩阈值算法[J].计算机技术与发展,2017,27(10):69-73.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.
[26]
张广智,孙昌路,潘新朋,等.快速共轭梯度法频率域声波全波形反演[J].石油地球物理勘探,2016,51(4):730-737.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.
[27]
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.
[28]
董士骑,韩立国,胡勇,等.基于MCPML边界条件的频率域声波方程高精度模拟[J].石油地球物理勘探,2018,53(1):47-54.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.
[29]
成景旺,吕晓春,顾汉明,等.基于柯西分布的频率域全波形反演[J].石油地球物理勘探,2014,49(5):940-945.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.
[30]
Polyak B T.The conjugate gradient method in extremal problems[J].USSR Computational Mathematics & Mathematical Physics,1969,9(4):94-112.
[31]
Sirgue L,Pratt R G.Efficient waveform inversion and imaging:A strategy for selecting temporal frequencies[J].Geophysics,2004,69(1):231-248.