弹性波传播的低数值频散波场模拟方法
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
压制数值频散以提高计算精度是检验地震波数值模拟方法的一个重要标准。基于弹性波传播方程,建立了低数值频散波场模拟的八阶FNRK方法。该方法以Runge-Kutta方法对时间导数进行三阶离散,以近似解析离散算子替代差分算子对空间偏导数进行八阶离散,结合通量校正传输技术消除离散后的数值频散。弹性波场模拟结果表明,与高阶有限差分方法相比,该方法能在压制数值频散方面具有明显的优势,计算精度提高,且适应于地震波在大规模复杂介质中传播的波场模拟
Suppressing the numerical dispersion to enhance the computational accuracy is an important standard to check a numerical simulation method of seismic wave.Based on the elastic wave equations,the eighth-order FNRK method is established using the low numerical dispersion.The new method employs a third-order Runge-Kutta method to discretize the temporal derivatives,uses the eighth-order nearly-analytic discrete operator instead of the finite-difference operator to discretize the high-order spatial derivatives,and uses the flux-corrected transport technique to eliminate the discretized numerical dispersion.The results of the elastic wave-field simulation show that this method enhances the computational accuracy and has an obvious superiority in eliminating the numerical dispersion,compared with the traditional high-order finite-difference method.At the same time,this new method has a strong adaptability for simulating the seismic wave propagation in the large-scale complex media.
引文
[1]Chapman C H,Shearer P M.Ray tracing in azimuthally anisotropic media-Ⅱ:Quasi shear wave coupling[J].Geophys J,1989,96:65-83.
    [2]杨顶辉.双相各向异性介质中弹性波方程的有限元解法及波场模拟[J].地球物理学报,2002,45(4):575-583.Yang D H.Finite element method of the elastic wave equation and wave field simulation in two-phase anisotropic media[J].Chinese Journal of Geophysics,2002,45(4):575-583.(In Chinese)
    [3]Komatitsch D,Vilotte J P.The spectral element method:An efficient tool to simulate the seismic responses of 2Dand 3Dgeological structures[J].Bull Seism Soc Am,1998,88:368-392.
    [4]Chen X F.A systematic and efficient method of computing normal modes for multi-layered half space[J].Geophys J Int,1993,115:391-409.
    [5]Kelly K R,Wave R W,Treitel S.Synthetic seismograms:A finite-difference approach[J].Geophysics,1976,41:2-27.
    [6]Dablain M A.The application of high-order differencing to scalar wave equation[J].Geophysics,1986,51:54-66.
    [7]董良国,马在田,曹景忠,等.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,2000,43(3):411-419.Dong L G,Ma Z T,Cao J Z,et al.A staggered-grid high-order difference method of one-order elastic wave equation[J].Chinese Journal of Geophysics,2000,43(3):411-419.(In Chinese)
    [8]Book D L,Boris J P,Hain K.Flux-corrected transport,II:generalization of the method[J].J Comput Phys,1975,18(3):248-283.
    [9]Fei T,Larner K.Elimination of numerical dispersion in finite-difference modeling and migration by fluxcorrected transport[J].Geophysics,1995,60(6):1830-1842.
    [10]杨顶辉,滕吉文.各向异性介质中三分量地震记录的FCT有限差分模拟[J].石油地球物理勘探,1997,32(2):181-190.Yang D H,Teng J W.FCT finite difference modeling of three-component seismic records in anisotropic medium[J].OGP,1997,32(2):181-190.(In Chinese)
    [11]Yang D H,Liu E,Zhang Z J,et al.Finitedifference modeling in two-dimensional anisotropic media using a flux-corrected transport technique[J].Geophys J Int,2002,148:320-328.
    [12]郑海山,张中杰.横向各向同性(VTI)介质中非线性地震波场模拟[J].地球物理学报,2005,48(3):660Zheng H S,Zhang Z J.Synthetic seismograms of nonlinear seismic waves in isotropic(VTI)media[J].Chinese Journal of Geophysics,2005,48(3):660-671.(In Chinese)
    [13]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[J].地球物理学进展,2005,20(1):58-65.Wu G C,Wang H Z.Analysis of numerical dispersion in wave-field simulation[J].Progress in Geophysics,2005,20(1):58-65.(In Chinese)
    [14]Zheng H S,Zhang Z J,Liu E.Non-linear seismic wave propagation in anisotropic media using fluxcorrected transport technique[J].Geophys J Int,2006,165:943-956.
    [15]王珺,杨长春,冯英杰.用优化通量校正传输技术压制数值模拟的频散[J].勘探地球物理进展,2007,30(4):252-256.Wang J,Yang C C,Feng Y J.Elimination of numerical dispersion in finite-difference modeling by optimized flux-corrected transport[J].PEG,2007,30(4):252-256.(In Chinese)
    [16]陈可洋,杨微,刘洪林,等.二阶弹性波动方程高精度交错网格波场分离数值模拟[J].物探与化探,2009,33(6):700-704.Chen K Y,Yang W,Liu H L,et al.Wave field separation numerical modeling of second order elastic wave equation by high-precision staggered-grid finite difference scheme[J].Geophysical&Geochemical exploration,2009,33(6):700-704.(In Chinese)
    [17]Kondoh Y,Hosaka Y,Ishii K.Kernel optimum nearly-analytical discretization algorithm applied to parabolic and hyperbolic equations[J].Computers Math Appl,1994,27(3):59-90.
    [18]Yang D H,Teng J W,Zhang Z J,et al.A nearlyanalytic discrete method for acoustic and elastic wave equations in anisotropic media[J].Bull Seism Soc Am,2003,93(2):882-890.
    [19]Yang D H,Peng J M,Lu M,et al.Optimal nearlyanalytic discrete approximation to the scalar wave equation[J].Bull Seism Soc Am,2006,96(3):1114-1130.
    [20]Yang D H,Song G J,Chen S,et al.An improved nearly analytical discrete method:an efficient tool to simulate the seismic response of 2-D porous structures[J].J Geophys Eng,2007,4:40-52.
    [21]Yang D H,Wang N,Chen S,et al.An explicit method based on the implicit Runge-Kutta algorithm for solving the wave equations[J].Bull Seism Soc Am,2009,99(6):3340-3354.
    [22]王磊,杨顶辉,邓小英.非均匀介质中地震波应力场的WNAD方法及其数值模拟[J].地球物理学报,2009,52(6):1526-1535.Wang L,Yang D H,Deng X Y.A WNAD method for seismic stress-field modeling in heterogeneous media[J].Chinese J Geophys,2009,52(6):1526-1535.(In Chinese)
    [23]Tong P,Yang D H,Hua B L,et al.A high-order stereo-modeling method for solving wave equations[J].Bulletin of the Seismological Society of America,2013,103(2):811-833.
    [24]Zhang C Y,Li X,Ma X,et al.A Runge-Kutta method with using eighth-order nearly-analytic spatial discretization operator for solving a 2Dacoustic wave equation[J].Journal of Seismic Exploration,2014,23(3):279-302.
    [25]Zhang C Y,Ma X,Yang L,et al.Symplectic partitioned Runge-Kutta method based on the eighthorder nearly analytic discrete operator and its wavefield simulations[J].Applied Geophysics,2014,11(1):89-106.
    [26]Qiu J X,Li T G,Khoo B C.Simulations of compressible two-medium flow by Runge-Kutta discontinuous Galerkin methods with the ghost fluid method[J].Commun Computs Phys,2008,3:479-504.

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