求解声波方程的辛可分Runge-Kutta方法  

Symplectic partitioned Runge-Kutta method for solving the acoustic wave equation

作　　者：马啸[1] 杨顶辉[1] 张锦华[2] 

机构地区：[1]清华大学数学科学系,北京100084 [2]昆明工业职业技术学院,昆明650302

出　　处：《地球物理学报》2010年第8期1993-2003,共11页Chinese Journal of Geophysics

基　　金：国家杰出青年科学基金(40725012)资助

摘　　要：本文基于声波方程的哈密尔顿系统,构造了一种新的保辛数值格式,简称NSFRK方法.该方法在时间上采用二阶辛可分Runge-Kutta方法,空间上采用近似解析离散算子进行离散逼近.针对本文发展的新方法,我们给出了NSPRK方法在一维和二维情况下的稳定性条件、一维数值频散关系以及二维数值误差,并在计算效率方面与传统辛格式和四阶LWC方法进行了比较.最后,我们将本文方法应用于声波在三层各向同性介质和异常体模型中的波传播数值模拟.数值结果表明,本文发展的NSPRK方法能有效压制粗网格或具有强间断情况下数值方法所存在的数值频散,从而极大地提高了计算效率,节省了计算机内存.In this paper, we develop a new symplectic numerical scheme based on Hamihonian system of the acoustic wave equation, which is called the NSPRK method in brief. The NSPRK method uses the nearly-analytic discrete operators to approximate the high-order differential operators, and employs the second-order symplectic partitioned Runge-Kutta method to numerically solve the Hamiltonian system. For the proposed NSPRK method in this paper, we obtain the stability conditions for 1D and 2D cases, the numerical dispersion relation for the 1D case and 2D numerical errors. Meanwhile, we compare the NSPRK against the conventional sympleetic method and the fourth-order LWC method in computational efficiency. Finally, we apply the NSPRK method to model acoustic wave propagating in a three-layer isotropic medium and an abnormal body model. The promising numerical results illustrate that the NSPRK method can effectively suppress the numerical dispersion caused by discretizing the acoustic-wave equation when coarse grids are used or models have large velocity contrasts between adjacent layers.Therefore, the NSPRK method can greatly increase the computational efficiency and save computer memory.

关 键 词：哈密尔顿系统 辛算法 近似解析离散 波场模拟 数值频散 

分 类 号：P631[天文地球—地质矿产勘探] O29[天文地球—地质学]