基于优化方法的时间-空间域隐格式有限差分算子确定方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
声波方程数值模拟构成了地震逆时偏移成像技术和全波形反演的基础。对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应。如何有效地压制数值频散是有限差分方法研究的关键所在。为了进一步抑制数值频散,利用隐式有限差分比显式有限差分更能压制数值频散的特点,采用前人提出的新的有限差分模板(在保持相同精度的情况下增大了时间步长),应用信赖域优化方法在时间-空间域确定隐格式有限差分系数。频散分析和数值模拟试算的结果表明,这种新模板隐格式有限差分优化方法既提高了声波数值模拟精度又提高了计算效率。
Numerical simulation of acoustic wave equation is the basis of the reverse time migration and full waveform inversion.With some stability conditions,grid dispersion often exists because of the discretization of the time and the spatial derivatives in the wave equation.How to suppress the grid dispersion is therefore a key problem for finite difference approaches.Finite-difference scheme of partial differential equations are either explicit or implicit.Implicit finite difference scheme can suppress the grid dispersion further.Meanwhile,a new finite difference stencil was proposed to increase the time step while preserving the accuracy by previous researchers.Therefore,we propose to determine the optimal implicit finite difference operators in time-space domain with a new finite difference stencil.Dispersion analysis and seismic numerical simulation demonstrate that the accuracy of acoustic wave numerical simulation and the computation efficiency can be improved by the method.
引文
[1]Dablain M A.The application of high-order differencing to the scalar wave equation[J].Geophysics,1986,51(1):54-66
    [2]Liu Y.Globally optimal finite-difference schemes based on least squares[J].Geophysics,2013,78(4):T113-T132
    [3]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
    [4]Liu H,Dai N,Niu F,et al.An explicit time evolution method for acoustic wave propagation[J].Geophysics,2014,79(3):T117-T124
    [5]Liu Y,Sen M K.A practical implicit finite difference method:examples from seismic modeling[J].Journal of Geophysics Engineering,2009,6(3):231-249
    [6]赵建国,史瑞其,陈竞一,等.辅助微分方程完全匹配层在声波方程数值模拟中的应用[J].吉林大学学报(地球科学版),2014,44(2):675-682Zhao J G,Shi R Q,Chen J Y,et al.Application of auxiliary differential equation perfectly matched layer for numerical modeling of acoustic wave equations[J].Journal of Jilin University(Earth Science Edition),2014,44(2):675-682
    [7]Du X,Bancroft J C,Lines L R.Anisotropic reversetime migration for tilted TI media[J].Geophysical Prospecting,2007,55(6):853-869
    [8]Liu H W,Li B,Liu H,et al.The issues of prestack reverse time migration and solutions with Graphic Processing Unit implementation[J].Geophysical Prospecting,2012,60:906-918
    [9]Xie W,Yang D,Liu F,et al.Reverse-time migration in acoustic VTI media using a high-order stereo operator[J].Geophysics,2014,79(3):WA3-WA11
    [10]Yan H,Liu Y.Acoustic VTI modeling and pre-stack reverse-time migration based on the time-space domain staggered-grid finite-difference method[J].Journal of Applied Geophysics,2013,90(1):41-52
    [11]Yan H,Liu Y.Visco-acoustic prestack reverse-time migration based on the time-space domain adaptive high-order finite-difference method[J].Geophysical Prospecting,2013,61(5):941-954
    [12]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
    [13]梁文全,杨长春,王彦飞,等.用于声波方程数值模拟的时间-空间域有限差分系数确定新方法[J].地球物理学报,2013,56(10):3497-3506Liang W Q,Yang C C,Wang Y F,et al.Acoustic wave equation modeling with new time-space domain finite difference operators[J].Chinese Journal of Geophysics,2013,56(6):840-850
    [14]Liu Y,Sen M K.Time-space domain dispersion-relation-based finite-difference method with arbitrary even-order accuracy for the 2Dacoustic wave equation[J].Journal of Computational Physics,2013,232(1):327-345
    [15]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
    [16]Tan S,Huang L.A staggered-grid finite-difference scheme optimized in the time-space domain for modeling scalar-wave propagation in geophysical problems[J].Journal of Computational Physics,2014,276:613-634
    [17]Wang Y F,Yuan Y X.Convergence and regularity of trust region methods for nonlinear ill-posed inverse problems[J].Inverse Problems,2005,21(3):821-838
    [18]王彦飞.反演问题的计算方法及其应用[M].北京:高等教育出版社,2007:1-370Wang Y F.Computational methods for Inverse problems and applications[M].Beijing:Higher Education Press,2007:1-370

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心