The eighth-order frequency-domain NAD method and full-waveform inversion
HAN Rubing1,2, LANG Chao1
1. Key Laboratory of Deep-Earth Dynamics of Ministry of Natural Resources, Institute of Geology, Chinese Academy of Geological Sciences, Beijing 100037, China; 2. Institute of Geophysics and Geomatics, China University of Geosciences(Wuhan), Wuhan, Hubei 430074, China
Abstract:The full waveform inversion (FWI) is a high precise imaging method which uses wave equation and optimal algorithm to obtain physical parameters of underground media.The forward modeling is the basis of inversion.In order to further improve the efficiency of forward modelings,this paper proposes an eighth-order nearly-analytic discretization (NAD) method to discretize a 2D acoustic wave equation in the frequency domain.The construction of high-order NAD method is deduced in detail and a class of inexact rotated block triangular preconditioned Krylov subspace method is used to solve large-scale sparse linear algebraic equations obtained after the discretization.The wavefield simulation and numerical dispersion ana-lysis show the advantages of the proposed method in suppressing the numerical dispersion and improving the computational efficiency.Specifically the wavefield can be accurately recovered when the sampling points per wave length is less than two.This exceeds the limit of sampling rate.Finally,inversion researches are carried out with the constructed forward-modeling algorithm.The frequency domain full waveform inversion is performed in two typical layered models and the Marmousi mo-del,and the high resolution,high fidelity results are obtained.These test results and inversion error curves illustrate the effectiveness and applicability of the proposed method.
韩如冰, 郎超. 频率域八阶NAD有限差分模拟及全波形反演[J]. 石油地球物理勘探, 2019, 54(6): 1254-1266.
HAN Rubing, LANG Chao. The eighth-order frequency-domain NAD method and full-waveform inversion. Oil Geophysical Prospecting, 2019, 54(6): 1254-1266.
Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49(8):1259-1266.
[2]
Mora P.Nonlinear two-dimensional elastic inversion of multi offset seismic data[J].Geophysics,1987,52(9):1211-1228.
[3]
Pratt R G.Frequency-domain elastic wave modeling by finite differences:A tool for crosshole seismic imaging[J].Geophysics,1990,55(5):626-632.
[4]
Pratt R G.Inverse theory applied to multi-source cross-hole tomography[J].Geophysical Prospecting,1990,38(3):311-329.
[5]
Liu S L,Li X F,Wang W S,et al.A modified symplectic scheme for seismic wave modeling[J].Journal of Applied Geophysics,2017,99:28-36.
[6]
Liu S L,Li X F,Wang W S,et al.A new kind of optimal second-order symplectic scheme for seismic wave simulations[J].Science China:Earth Sciences,2014,57(4):751-758.
[7]
廖建平,刘和秀,戴世鑫,等.二维时间空间域和频率空间域声波全波形速度反演方法的对比研究[J].地球物理学进展,2017,32(5):2029-2034.LIAO Jianping,LIU Hexiu,DAI Shixin,et al.Research on comparisons of 2D acoustic wave full waveform velocity inversion in time-space domain and frequency-space domain[J].Progress in Geophysics,2017,32(5):2029-2034.
[8]
Lang C,Yang D H.A nearly analytic discrete method for solving the acoustic-wave equations in the frequency domain[J].Geophysics,2016,82(1):T43-T57.
[9]
张广智,孙昌路,潘新朋,等.快速共轭梯度法频率域声波全波形反演[J].石油地球物理勘探,2016,51(4):730-737.ZHANG Guangzhi,SUN Changlu,PAN Xinpeng,et al.Acoustic full waveform inversion in the frequency domain based on fast conjugate gradient method[J].Oil Geophysical Prospecting,2016,51(4):730-737.
[10]
王毓玮,董良国,黄超,等.降低弹性波全波形反演强烈非线性的分步反演策略[J].石油地球物理勘探,2016,51(2):288-294.WANG Yuwei,DONG Liangguo,HUANG Chao,et al.A multi-step strategy for mitigating severe nonli-nearity in elastic full-waveform inversion[J].Oil Geophysical Prospecting,2016,51(2):288-294.
[11]
Liu S L,Yang D H,Ma J.A modified symplectic PRK scheme for seismic wave modeling[J].Computers and Geosciences,2017,99:28-36.
[12]
张文生,庄源.频率域声波方程全波形反演[J].数值计算与计算机应用,2017,38(3):167-196.ZHANG Wensheng,ZHUANG Yuan.Full-waveform inversion based on the acoustic wave equation in the frequency[J].Journal on Numerical Methods and Computer Applications,2017,38(3):167-196.
[13]
Alterman Z S.Finite difference solutions to geophysical problems[J].Journal of Physics of the Earth,1968,16(Special):113-128.
[14]
姚振岸,孙成禹,唐杰,等.基于不同震源机制的黏弹各向异性微地震波场模拟[J].石油地球物理勘探,2017,52(1):63-70.YAO Zhen'an,SUN Chengyu,TANG Jie,et al.Micro-seismic forward modeling in viscoelastic anisotropic media based on different focal mechanisms[J].Oil Geo-physical Prospecting,2017,52(1):63-70.
[15]
Lysmer J,Drake L A.A finite element method for seismology[J].Methods of Computational Physics,1972,11:181-216.
[16]
Liu S L,Li X F,Wang W S,et al.A mixed-grid finite element method with PML absorbing boundary conditions for seismic wave modeling[J].Journal of Geophysics and Engineering,2014,11(5):055009.
[17]
Dan D K,Baysal E.Forward modeling by a Fourier method[J].Geophysics,1982,47(10):1402-1412.
[18]
Komatitsch D,Tsuboi S,Tromp J.The Spectral-element method in seismology//Seismic Earth:Array Analysis of Broadband Seismograms[M].American Geophysical Union,2005,205-227.
[19]
Yang D H,Wang L,Deng X Y.An explicit split-step algorithm of the implicit Adams method for solving 2D acoustic and elastic wave equations[J].Geophysical Journal International,2010,180(1):291-310.
[20]
Yang D H,Teng J W,Zhang Z J,et al.A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media[J].Bulletin of the Seismological Society of America,2003,93(2):882-890.
[21]
Yang D H,Peng J M,Lu M,et al.Optimal nearly analytic discrete approximation to the scalar wave equation[J].Bulletin of the Seismological Society of America,2006,96(3):1114-1130.
[22]
胡建林,宋维琪,张建坤,等.交错网格有限差分正演模拟的联合吸收边界[J].石油地球物理勘探,2018,53(5):914-920.HU Jianlin,SONG Weiqi,ZHANG Jiankun,et al.Joint absorbing boundary in the staggered-grid finite difference forward modeling simulation[J].Oil Geophysical Prospecting,2018,53(5):914-920.
[23]
Lang C,Ren Z R.Inexact rotated block triangular preconditioners for a class of block two-by-two matrices[J].Journal of Engineering Mathematics,2015,93(1):1-12.
[24]
Liu S L,Suardi I,Yang D H,et al.Tele seismic tra-veltime tomography of northern Sumatra[J].Geophy-sical Research Letters,2018,45(24):13231-13239.
Liu Q H.The PSTD algorithm:A time-domain me-thod requiring only two cells per wavelength[J].Microwave and Optical Technology Letters,1997,15(3):158-165.
[27]
Jerri A J.The Shannon sampling theorem-Its various extensions and applications:A tutorial review[J].Proceedings of the IEEE,2005,65(11):1565-1596.
[28]
Versteeg R.The Marmousi experience:Velocity model determination on a synthetic complex data set[J].The Leading Edge,1994,13(9):927-936.
[29]
Sirgue L,Pratt R G.Efficient waveform inversion and imaging:A strategy for selecting temporal frequencies[J].Geophysics,2004,69(24):231-248.