摘要
针对具有哈密尔顿系统的二维弹性波方程,使用具有八阶精度的近似解析离散算子离散空间高阶偏导数,利用有着二阶精度的辛分部Runge-Kutta方法离散时间导数,得到了八阶NAD-SPRK方法。利用该方法模拟二维弹性波分别在均匀横向各向同性介质、双层均匀各向同性介质、双层非均匀介质和含起伏间断面的三层介质中的传播。波场模拟结果表明,八阶NAD-SPRK方法能有效地压制粗网格和强间断介质条件下引起的数值频散,PML吸收边界条件同八阶NAD-SPRK方法有机结合能有效地衰减来自人工边界的反射波。
- 单位