Acoustic full waveform inversion in the frequency domain based on fast conjugate gradient method
Zhang Guangzhi1,2, Sun Changlu3, Pan Xinpeng1, Chen Hongliang4, Jiang Lanjie1, Wen Tiemin5
1. School of Geosciences, China University of Petroleum (East China), Qingdao, Shandong 266580, China;
2. Laboratory for Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong 266580, China;
3. China Tianchen Engineering Corporation (TCC), Tianjin 300400, China;
4. EnerTech-Drilling & Production Co., CNOOC, Tianjin 300452, China;
5. GRI, BGP Inc., CNPC, Zhuozhou, Hebei 072751, China
Abstract:Hessian matrix in the full waveform inversion is huge and the convergence of the gradient method is slow. To solve the problem, we propose in this paper a new algorithm, fast conjugate gradient (FCG) method. The method introduces a new variable to transform the conjugate gradient method, which accelerates the convergence and makes convergence more stable. The method needs little more calculation of multiplication. The proposed method is applied to the acoustic full waveform inversion in the frequency domain and it is also tested in the simple depression model and the thinning complex Marmousi model. The tests show that the proposed method can accelerate convergence while the resolution of deep layer is better.
张广智, 孙昌路, 潘新朋, 陈洪亮, 姜岚杰, 温铁民. 快速共轭梯度法频率域声波全波形反演[J]. 石油地球物理勘探, 2016, 51(4): 730-737.
Zhang Guangzhi, Sun Changlu, Pan Xinpeng, Chen Hongliang, Jiang Lanjie, Wen Tiemin. Acoustic full waveform inversion in the frequency domain based on fast conjugate gradient method. OGP, 2016, 51(4): 730-737.
Tarantola A. Inversion of seismic reflection data in the acoustic approximation.Geophysics,1984,49(8):1259-1266.
[2]
Mora P. Nonlinear two-dimensional elastic inversion of multioffset seismic data. Geophysics,1987,52(9):1211-1228.
[3]
Pratt R G, Shin C, Hick G J. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion. Geophysical Journal International, 1998, 133(2): 341-362.
[4]
Bunks C, Saleck F M, Zaleski S et al. Multiscale seismic waveform inversion. Geophysics,1995,60(5):1457-1473.
[5]
Sirgue L, Pratt R G. Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies. Geophysics, 2004, 69(1): 231-248.
[6]
Nocedal J.Updating quasi-Newton matrices with limited storage. Mathematics of Computation, 1980,35(151): 773-782.
[7]
Liu D and Nocedal J.On the limited memory BFGS method for large scale minimization.Mathematical Programming, 1989,5(1):503-528.
[8]
Wang Y, Rao Y. Reflection seismic waveform tomography. Journal of Geophysical Research: Solid Earth (1978-2012), 2009, 114(B3).
[9]
Beck A,Teboulle M. A fast iterative shrinkage-thresholding algorithm for linear inverse problems. SIAM Journal on Imaging Sciences, 2009, 2(1): 183-202.
[10]
Pratt R G. Inverse theory applied to multi-source cross-hole tomography. Geophysical Prospecting, 1990, 38(3): 311-329.
[11]
Plessix R E. A review of the adjoint-state method for computing the gradient of a functional with geophysical applications. Geophysical Journal International, 2006, 167(2): 495-503.
[12]
Davis T A. Algorithm 832: UMFPACK V4.3—an unsymmetric-pattern multifrontal method. ACM Transactions on Mathematical Software (TOMS), 2004, 30(2): 196-199.
Chen Yongrui, Li Zhenchun, Qin Ning et al. Full waveform inversion with wave equation bi-directional illuminationoptimization. Progress in Geophysics,2013,28(6) : 3015-3021.
Liu Lu,Liu Hong,Zhang Heng et al. Full waveform inversion based on modified quasi-Newton equation.Chinese Journal of Geophysics,2013, 56(7): 2447-2451.
Liu Guofeng,Liu Hong,Meng Xiaohong et al. Frequency-related factors analysis in frequency domain waveform inversion.Chinese Journal of Geophysics,2012, 55(4): 1345-1353.
Zhang Guangzhi,Wang Danyang,Yin Xingyao et al.Study on prestack seismic inversion using Markov Chain Monte Carlo.Chinese Journal of Geophysics,2011, 54(11): 2926-2932.