Joint absorbing boundary in the staggered-grid finite difference forward modeling simulation
Hu Jianlin1, Song weiqi1, Zhang Jiankun2, Xing Wenjun2, Xu Wenhui2
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266555, China; 2. Research Institute of Petroleum Exploration and Development, Jidong Oilfield Company, PetroChina, Tangshan, Hebei 063004, China
Abstract:The boundary condition is one of the key issues discussed in the simulation of the staggered-grid finite difference of the 3D acoustic equations.The perfectly matched layer absorbing boundary (PML) is stable and effective,but it needs a certain boundary thickness to make perfectly absorption which makes the model bigger and requires a huge computation.Higdon absorbing boundary condition can stably eliminate the incident wave at any angle on the boundary,but it becomes very complex after discretization when the higher order condition used.Based on the plane wave propagation in a PML medium,we obtain a new Higdon absorbing operator in the PML,and then get joint absorbing boundary conditions.In the test of the joint absorbing conditions,it is found that the joint absorbing conditions are effective with a thinner PML layer and smaller computation.
黄金强,李振春,黄建平等.TTI介质Lebedev网格高阶有限差分正演模拟及波型分离.石油地球物理勘探, 2017,52(5):915-927.Huang Jinqiang,Li Zhenchun,Huang Jianping et al.Lebedev grid high-order finite-difference modeling and elastic wave-mode separation for TTI media.OGP,2017,52(5):915-927.
[5]
姚振岸,孙成禹,谢俊法等.黏弹TTI介质旋转交错网格微地震波场模拟.石油地球物理勘探,2017,52(2):253-263.Yao Zhenan,Sun Chengyu,Xie Junfa et al.Micro-seismic forward modeling in visco-elastic TTI media based on rotated staggered grid finite-difference me-thod.OGP,2017,52(2):253-263.
[6]
李庆洋,黄建平,李振春等.起伏地表贴体全交错网格仿真型有限差分正演模拟.石油地球物理勘探,2015,50(4):633-642.Li Qingyang,Huang Jianping,Li Zhenchun et al.Undulating surface body-fitted grid seismic modeling based on fully staggered-grid mimetic finite diffe-rence.OGP,2015,50(4):633-642.
[7]
殷学鑫,刘洋.二维随机介质模型正演模拟及其波场分析.石油地球物理勘探,2011,46(6):862-872.Yin Xuexin,Liu Yang.Random medium 2-D modeling and its wavefield analysis.OGP,2011,46(6):862-872.
[8]
毛小平,何大伟,辛广柱等.厚层火山岩地震响应研究.石油地球物理勘探,2005,40(6):670-676.Mao Xiaoping,He Dawei,Xin Guangzhu et al.Study of seismic response for thick-layered volcanic rock.OGP,2005,40(6):670-676.
[9]
Smith W D.A nonreflecting plane boundary for wave propagation problems.Journal of Computational Physics,1974,15(4):492-503.
[10]
Cerjan C,Dan K,Kosloff R et al.A nonreflecting boundary condition for discrete acoustic and elastic wave equations.Geophysics,1985,50(4):705-708.
[11]
Burns D R.Acoustic and elastic scattering from seamounts in three dimensions:A numerical modeling study.Journal of the Acoustical Society of America,1992, 92(5):2784-2791.
[12]
Kosloff R,Kosloff D.Absorbing boundaries for wave propagation problems.Journal of Computational Physics,1986,63(2):363-376.
[13]
Clayton R,Engquist B.Absorbing boundary condi-tions for acoustic and elastic wave equations.Bulletin of the Seismological Society of America,1977,67(6):1529-1540.
[14]
Liao Z P,Wong H L,Yang B P et al.A transmitting boundary for transient wave analyses.Science in China Series A,1984,27(10):1063-1076.
[15]
Randall C J.Absorbing boundary condition for the elastic wave equation.Geophysics,1988,53(5):611-624.
[16]
Keys R G.Absorbing boundary conditions for acoustic media.Geophysics,1985,50(6):892-902.
[17]
Higdon R L.Absorbing boundary conditions for diffe-rence approximations to the multidimensional wave equation.Mathematics of Computation,1986,47(176):437-459.
[18]
Bécache E,Givoli D,Hagstrom T.High-order absorbing boundary conditions for anisotropic and convective wave equations.Journal of Computational Physics,2010,229:1099-1129
[19]
Berenger J P.A perfectly matched layer for the ab-sorption of electromagnetic waves.Academic Press Professional,1994.
[20]
Chew W C,Liu Q H.Perfectly matched layers forelsas to dynamics:a new absorbing boundary condition.Journal of Computational Acoustics,1996,4(4):72-79.
[21]
Roden J A,Gedney S D.Convolutional PML (CPML):an efficient FDTD implementation of the CFS-PML for arbitrary media.Microwave and Optical Technology Letters,2000,27(5):334-339.
[22]
Gedney S D,Zhao B.An auxiliary differential equation formulation for the complex-frequency shifted PML.Transactions on Antennas and Propagation,2010,58(3):838-847.
[23]
Liu Y and Sen M K.A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation.Geophysics,2010,75(2):A1-A6.
[24]
Liu Y and Sen M K.A hybrid absorbing boundary condition for elastic staggered-grid modeling.Geophysical Prospecting,2012,60(6):1114-1132.
[25]
Gao Y,Song H,Zhang J et al.Comparison of artificial absorbing boundaries for acoustic wave equation mode-lling.Exploration Geophysics,2017,48(1):76-93.
[26]
Huang Mingxi,Xue Tingxiao,Wang Youxue.Seismic wave field simulation on the symplectic and FCT algorithm.Progress in Geophysics,2015,30(2):565-570.
[27]
Zhang X B,Song P,Li J S et al.Analysis of the impact of the theoretical reflection coefficient on the perfectly matched layer absorbing boundary.Periodical of Ocean University of China,2016,46(6):84-89.
[28]
Higdon R L.Numerical absorbing boundary condi-tions for the wave equation.Mathematics of Computation,1987,49(179):65-90.