加权吸收边界条件的优化设计-石油物探.pdf

第 卷第 期 石 油 物 探 ， ４８ ３ Ｖｏｌ．４８Ｎｏ．３ 　 　 　 年 月 ， ２００９ ５ ＧＥＯＰＨＹＳＩＣＡＬＰＲＯＳＰＥＣＴＩＮＧＦＯＲＰＥＴＲＯＬＥＵＭ Ｍａ ２００９ ２４４ ｙ 文章编号： （ ） １０００ １４４１２００９０３ ０２４４ ０３ 　　 加权吸收边界条件的优化设计 杨 微，陈可洋 　 （大庆石油学院地球科学学院，黑龙江大庆 １６３３１８） 摘要：基于声波波动方程的旁轴近似吸收边界条件可以通过提高泰勒级数展开的旁轴近似阶数获得更好的吸收 效果，但这样不仅增加了差分计算的复杂性，同时还导致计算过程不稳定。针对这一问题，提出了优化加权吸收 边界条件。首先从理论上讨论了优化加权吸收边界条件，提出了加权因子概念，通过调节加权因子来拓展边界 反射完全吸收角度；然后利用脉冲响应和数值模拟对优化加权吸收边界条件的吸收效果进行了分析。结果表 明，选择合适的加权因子可以将边界处的完全吸收角度拓宽到 以上，反射系数误差不超过 。该方法简 ６０° ２．５％ 单易行，计算量小，而吸收效果比常规旁轴近似吸收边界条件更好。 关键词：地震波波动方程；旁轴近似；加权吸收边界条件；加权因子；脉冲响应 中图分类号： 文献标识码： Ｐ６３１．４ Ａ 针对数值模拟的边界效应问题人们进行了大 对加权因子进行定量研究。我们兼顾计算效率、计 　　 量的研究，消除边界效应的方法归纳起来主要有两 算精度和边界吸收效果等因素，在双重吸收边界条 种：一种是完全匹配层吸收边界条件，一种是基于 件基础上，给出了优化的加权吸收边界条件。 声波波动方程的旁轴近似吸收边界条件。 ［］ 完全匹配层吸收边界条件是 １ 于 方法基本理论 Ｂｅｒｅｎｅｒ ｇ １ 　 １９９４年提出的，并将其应用于时域有限差分法 （ ）麦克斯韦电磁场方程的数值模拟，后来被 地震波波动方程的时间域形式为 ＦＤＴＤ 人们扩展到声波和弹性波等地震波场数值模 ２（ ） （） 狌 狏 狌 狌 １ 狋狋＝ 狓狓＋ 狕狕 ［ ］ ２ ６ ～ 式中：为位移；为地震波在地下介质中的传播 拟 中。完全匹配层吸收边界条件