一种四级四阶龙格库塔时间高阶伪谱法声波方程模拟.pdf

　 第３８卷第４期 物　 探　 与　 化　 探 Ｖｏｌ．３８，Ｎｏ．４　 　 ２０１４年８月 ＧＥＯＰＨＹＳＩＣＡＬ ＆ ＧＥＯＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ Ａｕｇ．，２０１４　 ＤＯＩ：１０．１１７２０／ ｗｔｙｈｔ．２０１４．４．２６ － 杨怀英，唐小平，刘宽厚．一种四级四阶龙格—库塔时间高阶伪谱法声波方程模拟［Ｊ］．物探与化探，２０１４，３８（４）：７７４ ７８０．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０． １１７２０／ ｗｔｙｈｔ．２０１４．４．２６ Ｙａｎｇ ＨＹ，ＴａｎｇＸＰ，ＬｉｕＫＨ．Ａｆｏｕｒ⁃ｓｔａｇｅｆｏｕｒｔｈｏｒｄｅｒＲｕｎｇｅ⁃Ｋｕｔｔａｔｉｍｅｈｉｇｈ⁃ｏｒｄｅｒｐｓｅｕｄｏｓｐｅｃｔｒａｌｍｅｔｈｏｄｆｏｒａｃｏｕｓｔｉｃｅｑｕａｔｉｏｎｓｉｍｕｌａｔｉｏｎ［Ｊ］．Ｇｅｏｐｈｙｓ⁃ － ｉｃａｌ ａｎｄ Ｇｅｏｃｈｅｍｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ，２０１４，３８（４）：７７４ ７８０．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１４．４．２６ 一种四级四阶龙格—库塔时间高阶伪谱法声波方程模拟 杨怀英，唐小平，刘宽厚 （中国地质调查局 西安地质调查中心，陕西 西安　 ７１００５４） 摘 要：龙格—库塔法是常用于求解常微分方程（ＯＤＥ）的一项技术，该技术具有精度高、稳定性强等特点。 笔者以 一种新的四级四阶龙格—库塔法为基础，先将其四级计算公式合并为新的两级计算迭代公式，从而达到节约计算 内存的目的；再以此为基础推导出声波方程的时间高阶离散形式，并与伪谱法技术相结合，研究四级四阶龙格—库 塔时间高阶伪谱法声波高精度、高清晰度的波场模拟技术，进而研究该方法的稳定性与频散特性；最后，分别选取 均匀介质、层状介质和透镜体模型进行波场模拟试验。 模拟结果表明，该方法具有稳定性强、能有效去除频散，波 场清晰度高，能适应更宽的模拟参数范围等特点，是一种高效率的、具有较大应用潜力的波场模拟方法。 关键词：四级四阶龙格—库塔法；两级公式；伪谱法；声波方程；频散特性 － － － 中图分类号：Ｐ６３１．４　 　 　 文献标识码：Ａ　 　 　 文章编号：１０００ ８９１８（２０１４）０４ ０７７４ ０７ ［１１］ 　 　 伪谱法是一种精度高、空间频散少的场模拟方 Ｃｈｕ Ｃ Ｌ 提出正则化伪谱—拉普拉斯法实现了ＴＩ 法，常用于各种地震波的模拟与成像。 其基本算法 介质波动方程的模拟，该方法虽为时间二阶误差，却 ［１２］ 是：采用傅立叶变换以求取空间导数部分的计算，而 可有效地改进模拟精度；孙文博，孙赞东 进行了 ［１］ 应用有限差分求取时间部分的导数 。 常用的伪 伪谱法下的ＶＳＰ逆时偏移，该方法采用了空间四阶 谱法在时间导数的计算中，通常采用向前差分（一 伪谱法提升了成像效果。 以上方法存有两点不足： 阶速度—应力方程）或中心差分（二阶位移