Mixed staggered grid finite difference numerical simulation and reverse time migration of velocity-stress acoustic equation
HU Ziduo1,2, LIU Wei1,2, SONG Jiawen3, ZENG Qingcai4, TIAN Yancan1,2, HAN Linghe1,2
1. PetroChina Research Institute of Petroleum Exploration & Development-Northwest, Lanzhou, Gansu 730020, China; 2. PetroChina Key Laboratory of Reservoir Description, Lanzhou, Gansu 730020, China; 3. GRI, BGP Inc., CNPC, Zhuozhou, Hebei 072750, China; 4. PetroChina Research Institute of Petroleum Exploration and Development, Beijing 100083, China
Abstract:The conventional staggered grid finite difference method (C-SFD) is widely used to simulate the velocity-stress acoustic equation. Although the spatial difference operator can reach the order of 2M difference accuracy (M represents the number of sets of axial grid points equidistant from the center point of the spatial difference operator), the discretized difference acoustic equation only has the second-order difference accuracy, resulting in low simulation accuracy and poor stability. In this paper, the axial grid points and off-axial grid points are used to construct the spatial difference operator to approximate the first-order spatial partial derivative, and a mixed staggered grid finite difference method (M-SFD) suitable for the velocity-stress acoustic equation simulation is constructed. Based on the time-space dispersion relationship and Taylor series expansion, the difference coefficient solution equations are established, and the analytical solution of the difference coefficient is derived. The discretized difference acoustic equation given by M-SFD can reach 4th-order, 6th-order, 8th-order, or even any even order difference accuracy. The dispersion analysis shows that under the condition of specific Courant condition number value (such as Courant Condition number r=0.3), C-SFD cannot control the numerical dispersion error within 1% by adjusting the value of M. By adjusting the values of M and N(N represents the number of sets of off-axial grid points equidistant from the center point of the spatial difference operator), M-SFD can basically control the dispersion error within 1‰ or even within 0.1‰. Stability analysis shows that when M values are the same, the stability of M-SFD is stronger than that of C-SFD. Numerical simulation examples show that when the computational efficiency is almost the same, M-SFD can more effectively suppress numerical dispersion than C-SFD, resulting in higher simulation accuracy. M-SFD can also use a larger time sampling interval than C-SFD to achieve higher computational efficiency while holding higher simulation accuracy. This paper further extends M-SFD to inverse time migration. As a wave field propagation operator in inverse time migration, M-SFD can effectively eliminate imaging artifacts caused by numerical dispersion, thereby improving the imaging accuracy and resolution of deep structures.
CARCIONE J M.Wave Fields in Real Media:Wave Propagation in Anisotropic, An Elastic, Porous and Electromagnetic Media[M].3rd ed.Elsevier, Amsterdam, 2015.
[2]
吴建鲁, 吴国忱, 王伟, 等.基于等效交错网格的流固耦合介质地震波模拟[J].石油地球物理勘探, 2018, 53(2):272-279.WU Jianlu, WU Guochen, WANG Wei, et al.Seismic forward modeling in fluid-solid media based on equivalent staggered grid scheme[J].Oil Geophysical Prospecting, 2018, 53(2):272-279.
[3]
VIRIEUX J, CALANDRA H, PLESSIX R É.A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging[J].Geophysical Prospecting, 2011, 59(5):794-813.
[4]
BERKHOUT A J.Review paper:an outlook on the future of seismic imaging, part I:forward and reverse modelling[J].Geophysical Prospecting, 2014, 62(5):911-930.
[5]
GERHARD P, SHIN C, HICKS.Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J].Geophysical Journal International, 1998, 133(2):341-362.
[6]
VIRIEUX J, OPERTO S.An overview of full-waveform inversion in exploration geophysics[J].Geophysics, 2009, 74(6):WCC1-WCC26.
[7]
RESHEF M, KOSLOFF D, EDWARDS M, et al.Three-dimensional acoustic modeling by the Fourier method[J].Geophysics, 1988, 53(9):1175-1183.
[8]
MITTET R.On the pseudospectral method and spectral accuracy[J].Geophysics, 2021, 86(3):T127-T142.
[9]
MARFURT K J.Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations[J].Geophysics, 1984, 49(5):533-549.
[10]
韩德超, 刘卫华, 司文朋.三角网格间断有限元法弹性波模拟精度分析[J].石油地球物理勘探, 2021, 56(4):758-770.HAN Dechao, LIU Weihua, SI Wenpeng.Analysis of elastic wave simulation accuracy with discontinuous Galerkin finite element method based on triangular meshes[J].Oil Geophysical Prospecting, 2021, 56(4):758-770.
[11]
王辉, 何兵寿, 邵祥奇.一阶速度-胀缩-旋转弹性波方程交错网格数值模拟[J].石油地球物理勘探, 2022, 57(6):1352-1361.WANG Hui, HE Bingshou, SHAO Xiangqi.Numerical simulation of first-order velocity-dilatation-rotation elastic wave equation with staggered grid[J].Oil Geophysical Prospecting, 2022, 57(6):1352-1361.
[12]
王绍文, 宋鹏, 谭军, 等.弹性波数值模拟中的高斯型混合吸收边界条件及其GPU并行[J].石油地球物理勘探, 2021, 56(3):485-495.WANG Shaowen, SONG Peng, TAN Jun, et al.Gaussian-type weighted hybrid absorbing boundary for elastic wave simulation and its acceleration on GPU[J].Oil Geophysical Prospecting, 2021, 56(3):485-495.
[13]
REN Z, LIU Y.Acoustic and elastic modeling by optimal time-space-domain staggered-grid finite-difference schemes[J].Geophysics, 2015, 80(1):T17-T40.
[14]
彭更新, 刘威, 郭念民, 等.基于时空域交错网格有限差分法的应力速度声波方程数值模拟[J].石油物探, 2022, 61(1):156-165.PENG Gengxin, LIU Wei, GUO Nianmin, et al.A time-space domain dispersion-relationship-based staggered-grid finite-difference scheme for modeling the stress-velocity acoustic wave equation[J].Geophysical Prospecting for Petroleum 2022, 61(1):156-165.
[15]
胡建林, 宋维琪, 张建坤, 等.交错网格有限差分正演模拟的联合吸收边界[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.
[16]
DABLAIN M A.The application of high-order differencing to the scalar wave equation[J].Geophysics, 1986, 51(1):54-66.
[17]
LIU Y, SEN M K.A new time-space domain high-order finite-difference method for the acoustic wave equation[J].Journal of Computational Physics, 2009, 228(23):8779-8806.
[18]
LIU Y, SEN M K.Finite-difference modeling with adaptive variable-length spatial operators[J].Geophysics, 2011, 76(4):T79-T89.
[19]
任志明, 戴雪, 包乾宗, 等.有限差分的时间和空间频散及其对逆时偏移和全波形反演的影响[J].地球物理学报, 2021, 64(11):4166-4180.REN Zhiming, DAI Xue, BAO Qianzong, et al.Time and space dispersion in finite difference and its influence on reverse time migration and full-waveform inversion[J].Chinese Journal of Geophysics, 2021, 64(11):4166-4180.
[20]
FORNBERG B.Classroom note:calculation of weights in finite difference formulas[J].SIAM Review, 1998, 40(3):685-691.
[21]
LIU Y, SEN M K.Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes[J].Bulletin of the Seismological Society of America, 2011, 101(1):141-159.
[22]
GELLER R J, TAKEUCHI N.Optimally accurate second-order time-domain finite difference scheme for the elastic equation of motion:one-dimensional cases[J].Geophysical Journal International, 1998, 135(1):48-62.
[23]
CHU C, STOFFA P L.Determination of finite-difference weights using scaled binomial windows[J].Geophysics, 2012, 77(3):W17-W26.
[24]
ZHANG J, YAO Z.Optimized finite-difference operator for broadband seismic wave modeling[J].Geophysics, 2013, 78(1):A13-A18.
[25]
LIU Y.Globally optimal finite-difference schemes based on least squares[J].Geophysics, 2013, 78(4):T113-T132.
[26]
LIU Y.Optimal staggered-grid finite-difference schemes based on least-squares for wave equation modelling[J].Geophysical Journal International, 2014, 197(2):1033-1047.
[27]
LIU Y, SEN M K.Time-space domain dispersion-relation-based finite-difference method with arbitrary even-order accuracy for the 2D acoustic wave equation[J].Journal of Computational Physics, 2013, 232(1):327-345.
[28]
WANG E, LIU Y, SEN M K.Effective finite-difference modelling methods with 2-D acoustic wave equation using a combination of cross and rhombus stencils[J].Geophysical Journal International, 2016, 206(3):1933-1958.
[29]
胡自多, 贺振华, 刘威, 等.旋转网格和常规网格混合的时空域声波有限差分正演[J].地球物理学报, 2016, 59(10):3829-3846.HU Ziduo, HE Zhenhua, LIU Wei, et al.Scalar wave equation modeling using the mixed-grid finite-difference method in the time-space domain[J].Chinese Journal of Geophysics, 2016, 59(10):3829-3846.
[30]
JO C H, SHIN C, SUH J H.An optimal 9-point, finite-difference, frequency-space, 2-D scalar wave extrapolator[J].Geophysics, 1996, 61(2):529-537.
[31]
SHIN C, SOHN H.A frequency-space 2-D scalar wave extrapolator using extended 25-point finite-difference operator[J].Geophysics, 1998, 63(1):289-296.
[32]
胡自多, 刘威, 雍学善, 等.三维波动方程时空域混合网格有限差分数值模拟方法[J].地球物理学报, 2021, 64(8):2809-2828.HU Ziduo, LIU Wei, YONG Xueshan, et al.Mixed-grid finite-difference method for numerical simulation of 3D wave equation in the time-space domain[J].Chinese Journal of Geophysics, 2021, 64(8):2809-2828.
[33]
TAN S, HUANG L.An efficient finite-difference method with high-order accuracy in both time and space domains for modelling scalar-wave propagation[J].Geophysical Journal International, 2014, 197(2):1250-1267.
[34]
REN Z, LI Z, LIU Y, et al.Modeling of the acoustic wave equation by staggered-grid finite-difference schemes with high-order temporal and spatial accuracy[J].Bulletin of the Seismological Society of America, 2017, 107(5):2160-2182.
[35]
REN Z, LI Z.Temporal high-order staggered-grid finite-difference schemes for elastic wave propagation[J].Geophysics, 2017, 82(5):T207-T224.
[36]
ZHOU H, LIU Y, WANG J.Time-space domain scalar wave modeling by a novel hybrid staggered-grid finite-difference method with high temporal and spatial accuracies[J].Journal of Computational Physics, 2022, 455:111004.
[37]
ZHOU H, LIU Y, WANG J.Elastic wave modeling with high-order temporal and spatial accuracies by a selectively modified and linearly optimized staggered-grid finite-difference scheme[J].IEEE Transactions on Geoscience and Remote Sensing, 2022, 60:1-22.