Acoustic wave equation modeling with rectangle grid finite difference operator and its linear time space domain solution
Liang Wenquan1, Wang Yanfei2, Yang Changchun2
1. School of Resource Engineering, Longyan University, Longyan, Fujian 364000, China;
2. Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
Abstract:Finite difference methods are widely used in wave equation modeling,reverse time migration,and full waveform inversion.Suppressing the grid dispersion is one of the key points for finite difference approaches.Many current methods determining the finite difference operators assume that the intervals are identical for the spatial partial derivative in different directions.However,to use finite difference operators with different spatial intervals in different directions is more suitable for wave field extrapolation.Therefore,we use a liner method in this paper to determine the finite difference coefficients with rectangle grid.Through dispersion analysis and numerical simulations,the proposed method is proved to be more accurate and efficient.Therefore,the rectangle grid finite difference operator and its linear solution might be applied in wave equation extrapolation instead of conventional methods.
梁文全, 王彦飞, 杨长春. 声波方程数值模拟矩形网格有限差分系数确定法[J]. 石油地球物理勘探, 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.
Dablain M A.The application of high-order differencing to the scalar wave equation.Geophysics,1986,51(1):54-66.
[2]
De Basabe J D and Sen M K.Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations.Geophysics,2007,72(6):T81-T95.
[3]
De Basabe J D and Sen M K.A comparison of finite-difference and spectral-element methods for elastic wave propagation in media with a fluid-solid interface.Geophysical Journal International,2015,200(1):278-298.
[4]
Liu H W,Li B,Liu H et al.The issues of prestack reverse time migration and solutions with graphic processing unit implementation.Geophysical Prospecting,2012,60(5):906-918.
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.
Liang Wenquan,Yang Changchun,Wang Yanfei et al.Acoustic wave equation modeling with new time-space domain finite difference operators.Chinese Journal of Geophysics,2013,56(10):3497-3506.
Zhang Zhiyu,Tan Xianbo,Huang Luyao et al.Anti-dispersion finite difference simulation and reverse-time migration for wave equations.OGP,2014,49(6):1115-1121.
Tang X,Liu H,Zhang H et al.An adaptable 17-point scheme for high-accuracy frequency-domain acoustic wave modeling in 2D constant density media.Geophysics,2015,80(6):T211-T221.
[10]
Chu C and Stoffa P L.Implicit finite-difference simulations of seismic wave propagation.Geophysics,2012,77(2):T57-T67.
[11]
Chen H,Zhou H and Li Y.Application of unsplit convolutional perfectly matched layer for scalar arbitrarily wide-angle wave equation.Geophysics,2014,79(6):T313-T321.
[12]
Liu Y,Sen M K.Acoustic VTI modeling with a time-space domain dispersion-relation-based finite-difference scheme.Geophysics,2010,75(3):A11-A17.
[13]
Yang D,Lu M,Wu R et al.An optimal nearly analytic discrete method for 2D acoustic and elastic wave equations.Bulletin of the Seismological Society of America,2004,94(5):1982-1992.
[14]
Zhang J H,Yao Z X.Optimized finite-difference operator for broadband seismic wave modeling.Geophy-sics,2012,78(1):A13-A18.
[15]
Liu H F,Dai H X,Niu F L et al.An explicit time evolution method for acoustic wave propagation.Geophysics,2014,79(3):T117-T124.
[16]
Wang Y F,Liang W Q,Nashed Z et al.Seismic modeling by optimizing regularized staggered-grid finite-difference operators using a time-space-domain dispersion-relationship-preserving method.Geophysics,2014,79(5):T277-T285.
[17]
Liu Y,Sen M K.A new time-space domain high-order finite-difference method for the acoustic wave equation.Journal of Computational Physics,2009,228(23):8779-8806.
[18]
王彦飞,斯捷潘诺娃I E,提塔连科V N等.地球物理数值反演问题.北京:高等教育出版社,2011.
[19]
Liu Y,Sen M K.A hybrid scheme for absorbing edge reflections in numerical modeling of wave propagation.Geophysics,2010,75(2):A1-A6.