Abstract:The 9-point scheme for the second-order accuracy of the hybrid grid cannot meet the actual needs of acoustic equation modeling in the frequency domain because of the complex structure of the underground medium.A fourth-order 17-point difference scheme based on the 9-point with fourth-order differential accuracy is derived,and a nonsplit multi-axial convoluntion perfectly matched layer (MCPML) absorbing boundary is applied.Compared with conventional perfectly matched layer (PML) absorbing boundary,MCPML combines the advantages of convolution perfectly matched layer (CPML) absorbing boundary and multi-axial perfectly matched layer (MPML) absorbing boundary to improve the absorption of wave field with large incident angle without increasing computational burden.Meanwhile,the non-split method is easy to program.The proposed method is applied to the Marmousi model and its effectiveness is verified.
董士琦, 韩立国, 胡勇, 罗玉钦. 基于MCPML边界条件的频率域声波方程高精度模拟[J]. 石油地球物理勘探, 2018, 53(1): 47-54.
Dong Shiqi, Han Liguo, Hu Yong, Luo Yuqin. Acoustic equation high-accuracy modeling in the frequency domain based on MCPML absorbing boundary. Oil Geophysical Prospecting, 2018, 53(1): 47-54.
Lysmer J,Drake L A.A finite element method for seismology.Methods in Computational Physics:Advances in Research & Applications,1972,11:181-216.
[2]
Marfurt K J.Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations. Geophysics,1984,49(4):533-549.
[3]
Marfurt K J,Shin C S.The future of iterative modeling in geophysical exploration//Supercomputers in Seismic Exploration(Handbook of Geophysical Exploration:Seismic Exploration,Volume 21).Pergamon Press,New York,1989,203-228.
[4]
Pratt R G.Frequency-domain elastic wave modeling by finite difference:A tool for crosshole seismic imaging.Geophysics,1990,55(5):626-632.
[5]
Jo C H,Shin C S,Suh J H.An optimal 9-point,finite difference,frequency-space,2-D scalar wave extrapolator.Geophysics,1996,61(2):529-537.
[6]
Stekl I,Pratt R G.Accurate viscoelastic modeling by frequency-domain finite differences using rotated operators.Geophysics,1998,63(5):1779-1794.
OpertoS,VirieuxJ,RibodettiAetal.Finite-difference frequency-domain modeling of viscoacoustic wave propagation in 2D tilted transversely isotropic(TTI) media.Geophysics,2009,74(5):T75-T95.
[10]
Dupuy B,De Barros L,Garambois S et al.Wave propa-gation in heterogeneous porous media formulated in the frequency-space domain using a discontinuous Galerkin method.Geophysics,2011,76(4):N13-N28.
[11]
Hustedt B,Operto S,Virieux J.Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modelling.Geophysical Journal International,2004,157(3):1269-1296.
[12]
Min D J,Shin C,Kwon B D et al.Improved frequency-domain elastic wave modeling using weighted-averaging difference operators.Geophysics,2000,65(3):884-895.
[13]
马召贵,王尚旭,宋建勇.频率域波动方程正演中的多网格迭代算法.石油地球物理勘探,2010,45(1):1-5.Ma Zhaogui,Wang Shangxu,Song Jianyong.Multi-grid iterative algorithm in frequency domain wave equation modeling.OGP,2010,45(1):1-5.
[14]
韩利,韩立国,李翔等.二阶声波方程频域PML边界条件及频域变网格步长并行计算.吉林大学学报(地球科学版),2011,41(4):1226-1232.Han Li,Han Liguo,Li Xiang et al.PML boundary conditions for second-order acoustic wave equations and variable grid parallel computation in frequency-domain modeling.Journal of Jilin University(Earth Science Edition),2011,41(4):1226-1232.
[15]
高正辉,孙建国,孙章庆等.基于完全匹配层构建新方法的2.5维声波近似方程数值模拟.石油地球物理勘探,2016,51(6):1128-1133.Gao Zhenghui,Sun Jianguo,Sun Zhangqing et al.2.5D acoustic approximate equation numerical simulation with a new construction method of the perfectly matched layer.OGP,2016,51(6):1128-1133.
[16]
蔡志成,顾汉明,成景旺等.大规模变网格三维地震正演MPI并行策略与实现.石油地球物理勘探,2017,52(3):468-476.Cai Zhicheng,Gu Hanming,Cheng Jingwang et al.Parallel large-scale variable grid 3D seismic wave modeling.OGP,2017,52(3):468-476.
[17]
殷文,印兴耀,吴国忱等.高精度频率域弹性波方程有限差分方法及波场模拟.地球物理学报,2006,49(2):561-568.Yin Wen,Yin Xingyao,Wu Guochen et al.The me-thod of finite difference of high precision elastic wave equations in the frequency domain and wave field simulation.Chinese Journal of Geophysics,2006,49(2):561-568.
[18]
刘璐,刘洪,刘红伟.优化15点频率-空间域有限差分正演模拟.地球物理学报,2013,56(2):644-652.Liu Lu,Liu Hong,Liu Hongwei.Optimal 15-point finite difference forward modeling in frequency-space domain.Chinese Journal of Geophysics,2013,56(2):644-652.
[19]
梁文全,王彦飞,杨长春.声波方程数值模拟矩形网格有限差分系数确定法.石油地球物理勘探,2017,52(1):56-62.Liang Wenquan,Wang Yanfei,Yang Changchun.Acoustic wave equation modeling with rectangle grid finite difference operator and its linear time space domain solution.OGP,2017,52(1):56-62.
[20]
周聪,刘江平,罗银河等.二维频率域全波场有限差分数值模拟方法.石油地球物理勘探,2014,49(2):278-287.Zhou Cong,Liu Jiangping,Luo Yinhe et al.2D full-wavefield modeling in frequency domain using finite-difference.OGP,2014,49(2):278-287.
[21]
马晓娜,李志远,谷丙洛等.2D声波频率域数值模拟中几种有限差分方法的对比分析.地球物理学进展,2015,30(2):878-888.Ma Xiaona,Li Zhiyuan,Gu Bingluo et al.Comparisons and analysis of several optimization finite-differencing schemes in 2D acoustic frequency-domain numerical modeling.Progress in Geophysics,2015,30(2):878-888.
[22]
曹红书,陈景波.声波方程频率域高精度正演的17点格式及数值实现.地球物理学报,2012,55(10):3440-3449.Cao Hongshu,Chen Jingbo.A 17-point scheme and its numerical implementation for high-accuracy modeling of frequency-domain acoustic equation.Chinese Journal of Geophysics,2012,55(10):3440-3449.
[23]
Bérenger J P.A perfectly matched layer for the absorp-tion of electromagnetic waves.Journal of Computational Physics,1994,114(9):185-200.
[24]
Chew W C,Weedon W H.A 3-D perfectly matched medium from modified Maxwell's equations with stretched coordinates.Microwave & Optical Techno-logy Letters,1994,7(13):599-604.
[25]
Engquist B,Majda A.Absorbing boundary conditions for numerical simulation of waves.Mathematics of Computation, 1977,31(139):629-651.
[26]
Shin C.Sponge boundary condition for frequency-domain modeling.Geophysics,1995,60(6):1870-1874.
[27]
Higdon R L.Absorbing boundary conditions for elastic wave.Geophysics,1991,56(2):231-241.
[28]
Liao Z P,Wang H L,Yang B P et al.A transmitting boundary for transient wave analysis.Science China(Series A),1984,27(10):1063-1076.
[29]
Kuzuoglu M,Mittra R.Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers.IEEE Microwave & Guided Wave Letters,2002,6(12):447-449.
[30]
Meza-Fajardo K C,Papageorgiou A S.A non-convolutional split-field,perfectly matched layer for wave propagation in isotropic and anisotropic elastic media:stability analysis.Bulletin of the Seismological Society of America,2008,98(4):1811-1836.
[31]
Hu W,Abubakar A,Habashy T.Application of the nearly perfectly matched layer in acoustic wave mo-deling.Geophysics,2006,72(5):SM169-SM175.
[32]
Gao H,Zhang J.Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modeling.Geophysical Journal International,2008,174(3):1029-1036.
[33]
NilssonS,PeterssonN,SjögreenBetal.Stablediffe-rence approximations for the elastic wave equation in second order formulation.Siam Journal on Numerical Analysis,2007,45(5):1902-1936.
[34]
Martin R,Komatitsch D,Bruthiaux E et al.Convolution and nonconvolution perfectly matched layer techniques optimized at grazing incidence for high-order wave propagation modelling.Geophysical Research Abstracts,2010,12:EGU2010-4268-2.
[35]
Meza-Fajardo K C,Papageorgiou A S.On the stability of a non-convolutional perfectly matched layer for isotropic elastic media.Soil Dynamics & Earthquake Engineering,2010,30(3):68-81.
[36]
Zeng C,Xia J H,Miller R D et al.Application of the multiaxial perfectly matched layer (MPML) to near-surface seismic modeling with Rayleigh waves.Geophysics,2011,76(3):T43-T52.
张显文,韩立国,黄玲等.基于递归积分的复频移PML弹性波方程交错网格高阶差分.地球物理学报,2009,52(7):1800-1807.Zhang Xianwen,Han Liguo,Huang Ling et al.A staggered-grid high-order difference method of complex frequency-shifted PML based on recursive intergration for elastic wave equation.Chinese Journal of Geophysics,2009,52(7):1800-1807.
[39]
刘洋.波动方程时空域有限差分数值解及吸收边界条件研究进展.石油地球物理勘探,2014,49(1):35-46.Liu Yang.Research progress on time-space domain finite difference numerical solution and absorption boundary conditions of wave equation.OGP,2014,49(1):35-46.
Drossaert F H,Giannopoulos A.A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves.Geophysics,2007,72(2):T9-T17.
[44]
Zhang W,Shen Y.Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling.Geophysics,2010,75(4):T141-T154.