瑞雷波有限差分数值模拟中不同自由表面边界条件-石油地球物理勘探.PDF

年 月 第 卷 第 期 ２０１７ １２ ５２ ６ 　 　 　 ·地震模拟 · 文章编号： （ ） １０００７２１０２０１７０６１１５６１４ 瑞雷波有限差分数值模拟中不同自由表面 边界条件的对比 袁士川① ①② ① ① ① 宋先海 蔡 伟 胡 莹 鲁 鹏 　 　 　 　 　 　 　 （ 中国地质大学地球物理与空间信息学院，湖北武汉 ４３００７４； 地球内部多尺度成像湖北省重点实验室，湖北武汉 ４３００７４） ① ② 袁士川，宋先海，蔡伟，胡莹，鲁鹏 瑞雷波有限差分数值模拟中不同自由表面边界条件的对比 石油地球物理勘 ． ． 探， ， （）： ２０１７５２６ １１５６１１６９． 摘要 自由表面边界条件是决定瑞雷波数值模拟效果的关键因素。本文基于标准交错网格高阶有限差分法， 　 在二维各向同性弹性介质背景下，针对应力镜像法（ ）、改进应力镜像法（ ）、横向各向同性介质替换法 ＳＩＭ ＭＳＩＭ （ ）和声学—弹性边界近似法（ ）等四种最具代表性的自由表面边界条件进行了数值模拟，并在均匀半空 ＭＳ ＡＥＡ 间模型中从波场快照、波形曲线和频散曲线三个角度进行了对比分析。在相同条件下，上述四种方法均能生成 符合勘探地球物理规律的波场快照，各自对应的数值解与解析解的拟合误差都随网格剖分精度的提高而减小， 和 数值模拟的稳定性和精度都明显高于 和 。基于层状介质模型的进一步研究表明：对于 ＳＩＭ ＡＥＡ ＭＳＩＭ ＭＳ 简单模型， 和 都能得到比 和 更高精度的数值模拟结果；对于复杂模型， 的精度最高， ＳＩＭ ＡＥＡ ＭＳＩＭ ＭＳ ＡＥＡ 是最适合瑞雷波数值模拟的自由表面边界条件。 关键词 瑞雷波 有限差分模拟 自由表面边界条件 应力镜像法 声学—弹性边界近似法 　 　 　 　 　 中图分类号： 文献标识码： ： ／ Ｐ６３１ Ａ 犱狅犻１０．１３８１０ ．ｃｎｋｉ．ｉｓｓｎ．１０００７２１０．２０１７．０６．００５ 　　 　　 ｊ 等浅地表结构探测中得到广泛应用，其中最具代表 引言