求解弹性波方程辛RKN格式.pdfVIP

求解弹性波方程的辛 RKN 格式 1 1 2 3 刘少林 , 李小凡 ,汪文帅 , 朱童  1 中国科学院地质与地球物理研究所，中国科学院地球与行星物理重点实验室，北京 100029 2 宁夏大学数学与计算机学院，银川 750021 3 中国石油化工股份有限公司石油物探技术研究院, 南京 210014 摘要 将弹性波方程变换至Hamilton 体系, 构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nyström (RKN) 格式, 运用根数理论得到此格式的阶条件方程组. 通过给定系数的限定条件, 得到方程的对称解. 为了使时间离 散误差达到极小，提出数值频率与真实频率比较，通过Taylor 展开，得到关于辛系数的限定方程，求解方程组得 到最小频散辛RKN 格式. 对比分析时间演进方程的稳定性，得到使库朗数达到极大值的限定方程，求解方程组得 到最稳定辛RKN 格式. 发现此两种格式为同一格式. 新得到的辛RKN 格式不依赖于空间离散方法，为了对比的 需要，选取有限差分法进行空间离散. 在频散、稳定性分析中，与常见辛格式对比，从理论上分析了本文提出的 格式在数值频散压制、稳定性提升等方面的优势, 数值实验进一步证实了理论分析的正确性. 关键词 弹性波方程, 辛RKN 格式, 稳定性条件, 频散关系, 有限差分法 A symplectic RKN scheme for solving elastic wave equations 1 1 2 3 Liu Shao-Lin , Li Xiao-Fan , Wang Wen-Shuai , Zhu Tong 1 Key Laboratory of Earth and Planetary Physics, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China 2 School of Mathematics and Computer Science, Ningxia University, Yinchuan 750021, China 3 SINOPEC Geophysical Research Institute, Nanjing 210014, China Abstract Following the transform of the elastic seismic wave equation into Hamiltonian system, a explicit second-order symplectic Runge-Kutta-Nyström (RKN) scheme is proposed for high-efficient elastic wave modeling, and the order conditions are obtained by the rooted trees theory. To obtain symmetric solutions, we assign a constraint condition to the symplectic coefficients. In order to minimize the error in temporal space discretization, we propose a method that compares numerical angular frequency with real frequency, and obtain a constraint equation based on Taylor expansion. A symplectic RKN with minimized numerical dispersion is developed by solving the order con