伪谱法弹性波场数值模拟中的边界条件.pdfVIP

第 卷 第 期 地 球 物 理 学 进 展 ２２ ５ Ｖｏｌ．２２ Ｎｏ．５ 　 　 　 　 　 　 　 　 年 月（页码： ） ２００７ １０ １３７５ １３７９ ＰＲＯＧＲＥＳＳ ＩＮ ＧＥＯＰＨＹＳＩＣＳ Ｏｃｔ． ２００７ ～ 　 　 　 伪谱法弹性波场数值模拟中的边界条件 ， １２ ２ ２ 李信富 ， 李小凡 ， 张美根 　 　 （ 中国地质大学（北京），地下信息探测技术与仪器教育部重点实验室，北京 ； １． １０００８３ ２．中国科学院地质与地球物理研究所，岩石圈演化国家重点实验室，北京 １０００２９） 摘 要 边界条件问题长期以来一直是困扰地震波数值模拟研究者的一个难题，许多人都提出了有效的方法，但是 　 　 该问题仍然需要研究 本文针对伪谱法地震波场数值模拟的特殊要求，即所有网格点相互耦合，且傅立叶变换的周期 ． 性使得方程在网格边界上很难得到一个吸收的替代方程，研究衰减边界条件中衰减系数的变化对边界反射波衰减效 果的影响 数值实验结果说明，随着衰减系数的增大，振幅衰减加快，用很少的过渡带网格就可以使边界上的反射波 ． 能量很小 但是如果衰减率过大，就会在传播区和过渡带产生干扰波场 因此，在衰减系数的选择上应采取折中办法， ． ． 即在取某一较小的衰减系数的前提下，设置较少的过渡带网格数使边界的能量反射最少． 关键词 伪谱法，数值模拟，边界条件，衰减系数， 　 中图分类号 文献标识码 文章编号 （ ） Ｐ３１５ Ａ １００４２９０３２００７０５１３７５０５ 　 　　　　　 　 　　　　　 　 犅狅狌狀犱犪狉犮狅狀犱犻狋犻狅狀犳狅狉 狊犲狌犱狅狊犲犮狋狉犪犾狀狌犿犲狉犻犮犪犾 狔 狆 狆 狊犻犿狌犾犪狋犻狅狀狅犳狊犲犻狊犿犻犮狑犪狏犲 狉狅犪犪狋犻狅狀 狆 狆犵 ， １２ ２ ２ ， ， ＬＩＸｉｎｆｕ ＬＩＸｉａｏｆａｎ ＺＨＡＮＧＭｅｉｅｎ 　 　 ｇ （ ， ，