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

摘要

边界条件问题长期以来一直是困扰地震波数值模拟研究者的一个难题,许多人都提出了有效的方法,但是该问题仍然需要研究.本文针对伪谱法地震波场数值模拟的特殊要求,即所有网格点相互耦合,且傅立叶变换的周期性使得方程在网格边界上很难得到一个吸收的替代方程,研究衰减边界条件中衰减系数的变化对边界反射波衰减效果的影响.数值实验结果说明,随着衰减系数的增大,振幅衰减加快,用很少的过渡带网格就可以使边界上的反射波能量很小.但是如果衰减率过大,就会在传播区和过渡带产生干扰波场.因此,在衰减系数的选择上应采取折中办法,即在取某一较小的衰减系数的前提下,设置较少的过渡带网格数使边界的能量反射最少.