声波方程数值模拟矩形网格有限差分系数确定法-石油地球物理勘探.PDF

年 月 第 卷 第 期 ２０１７ ２ ５２ １ 　 　 　 ·地震模拟 · 文章编号： （ ） １０００７２１０２０１７０１００５６０７ 声波方程数值模拟矩形网格 有限差分系数确定法 ① ② ② 梁文全 王彦飞 杨长春 　 　 （ 龙岩学院资源工程学院，福建龙岩 ３６４０００； 中国科学院地质与地球物理研究所， ① ② 中国科学院油气资源研究重点实验室，北京 １０００２９） 梁文全，王彦飞，杨长春 声波方程数值模拟矩形网格有限差分系数确定法 石油地球物理勘探， ， （）： ． ． ２０１７５２１ ５６６２． 摘要 压制数值频散是有限差分方法的关键问题之一。目前压制数值频散的方法大多假设不同方向空间偏导 　 数的空间步长相同，导致算法精度低，计算效率低。为此，提出使用线性方法压制声波方程矩形网格有限差分 算子的数值频散，并进行了稳定性分析、频散分析和数值模拟。通过频散分析和数值模拟，验证了本文方法能 够有效压制矩形网格有限差分数值频散，相较于泰勒展开方法和最小二乘方法，线性方法计算有限差分系数的 效率更高，可以替代传统的正方形有限差分网格和相应的系数用于声波方程数值延拓。 关键词 声波模拟 时间—空间域 有限差分格式 矩形网格 数值频散 　 　 　 　 　 中图分类号： 文献标识码： ： ／ Ｐ６３１ Ａ 犱狅犻１０．１３８１０ ．ｃｎｋｉ．ｉｓｓｎ．１０００７２１０．２０１７．０１．００９ 　　 　　 ｊ 了 点矩形网格有限差分格式在频率—空间域声 １７ ［ ］ １０ １ 引言 波方程中的应用。Ｃｈｕ等 讨论了矩形网格（立方 　 体）隐格式有限差分地震波数值模拟。但是鲜见对 地震波数值模拟延拓是地震学和勘探地震学的 时间域矩形网格（立方体）有限差分算子系数如何确 ［］ １ 定的讨论，并且一般采用泰勒展开法确定有限差分 重要基础，在理论和应用中都发挥了重要的作用 。