Abstract:The Runge-Kutta method is a common technology for solving the Ordinary Differential Equations (ODE) and is characterized by high precision, strong stability and some other advantages. In this paper, based on a new four-stage fourth order Runge-Kutta method, the authors first combined the four-stage calculation formula with a new two-stage iteration formula, thus achieving the purpose of saving computational memory. And then, the time high-order discrete form of the acoustic wave equation was derived and, in combination with the pseudospectral method, some researches on the high-accuracy and high-definition acoustic wave field simulation technology of the four-stage fourth order Runge-Kutta time high-order pseudospectral method were carried out, with an investigation of the stability and dispersion of the method. Finally, homogeneous media, layered media and lens model were selected for wave field simulation test. The simulation results show that the four-stage fourth order Runge-Kutta time high-order pseudospectral method has strong stability and high wave field definition and can effectively remove the dispersion and adapt itself to large simulation parameter range, thus being a high efficient wave field simulation method with great application potential.
杨怀英, 唐小平, 刘宽厚. 一种四级四阶龙格—库塔时间高阶伪谱法声波方程模拟[J]. 物探与化探, 2014, 38(4): 774-780.
YANG Huai-Ying, TANG Xiao-Ping, LIU Kuan-Hou. A four-stage fourth order Runge-Kutta time high-order pseudospectral method for acoustic equation simulation. Geophysical and Geochemical Exploration, 2014, 38(4): 774-780.
[1] Dan D K,Edip B.Forward modeling by a fourier method[J].Geophysics,1982,47:1402-1412. [2] Bum K,George A.Two dimension elastic pseudo-spectral modeling of wide-Aperture seismic array date with application to the Wichita uplift Anadarko basin region of southwestern Oklahoma[J].Bull. Seism. Soc. Am.,1990,80:1677-1695. [3] Tesmer E T, Kosloff D B.A wave propagation simulation in presence of topography[J].Geophysical J Int,1992,108:621-632. [4] Tesmer E,Kosloff T D.3D elastic modeling with surface topography by a cheby-chev spectral method[J].Geophysics,1994,59(3):464-473. [5] Rosemary R,Jochen F A.Pseudospectral Chebychev method for the 2D wave equation with domain stretching and absorbing boundary conditions[J].J. comput. physics,1996,124:324-336. [6] Erding L,Heinz K O.Pseudospectral vs finite difference methods for initial value problems with discontinuous ceocents[J].SLAM.A.SCi.Compute.,1999,20:148-155. [7] Wang Y B,Takenaka H.A ultidomain approach of the Fourier pseudospectral method using discontinuous grid for elastic wave modeling[J].Earth Planets Space,2001,53:149-158. [8] Wang Y B,Hiroshi T,Takashi F.Effect of vertical velocity gradient on ground motion in a sediment-filled basin due to incident SV wave[J].Planetary and Space Sciences,2000,52:13-24. [9] Wang Y B,Hiroshi T,Takashi F.Modelling seismic wave propagation in a two dimensional cylindrical whole-earth model using the pseudospectral method[J].Geophysical.J.Int.,2001,145:689-708. [10] Zhao Z X,Xu J R,Horiuchi S.Differentiation operation in the wave equation for the pseudo-spectral method with a staggered mesh[J].Earth Planets Space,2001,53:327-332. [11] Chu C L.Acoustic anisotropic wave modeling using mormalized pseudolapla-cian//Denver:SEG 2010 Annual Meeting,2010. [12] Sun W B,Sun Z D,Zhu X H.Amplitude preserved VSP reverse time migration for angle-domain CIGs extraction[J].Applied Geophysics,2011,8 (2):141-149. [13] Gottlieb S,Ketcheson D,Shu C W.Strong stability preserving Runge-Kutta and multistep time discretization[M].World scientific publishing,2011. [14] Shu C W.Total-variation diminishing time discretizations[J].SIAM Journal on Scientific and Statistical Computing,1988,9:1073-1084. [15] Yang D H,Chen S,Li J Z.A Runge-Kutta method using high-order interpolation approximation for sloving 2D acoustic and elastic wave equation[J].Journal of seismic exploration,2007,16:331-2353. [16] 陈山,杨顶辉,邓小英.四阶龙格—库塔方法的一种改进算法及地震波场模拟[J].地球物理学报,2010,53(5):1196-1206. [17] Fornberg B.A practical guide to pseudospectral methods[M].Cambridge:Cambridge University Press,1996. [18] Mead J L,Renauty R A.Optimal Runge-Kutta methods for first order pseudo-spectral operators[J].Journal of Computational Physics.1999,152:404-419. [19] Javidia M,Golbabaib A.Spectral collocation method for parabolic partial differential equations with neumann boundary conditions[J].Applied Mathematical Sciences,2007,1(5):211-218. [20] Javidia M.Numerical solution for nonlinear PDEs[J].Int.J.Contemp. Math. Sciences,2007,2(8):373-381. [21] Ogundare B S.On the pseudo-spectral method of solving linear ordinary differential equations[J].Journal of Mathematics and Statistics,2009,5 (2):136-140.