摘要：在地震波场的波动方程数值模拟中，由于计算量的限制，必须加入人为的边界，使模拟计算可以在一定的空间范围内进行，由于边界节点上的波场值不能像模拟区域内部的节点一样使用中心差分来计算，使其计算精度大大降低，从而产生边界反射．为了消除边界反射，本文提出了延迟边界方法，根据弹性波在传播方向上等距离质点的等相位延迟性质和振幅衰减特性，由内部波场的时空分布，推算出边界波场的相位延迟的大小和振幅衰减系数，从而提高边界节点上的波场值计算精度，消除边界反射的产生。

关键词：地震波场 波动方程 数值模拟 延迟边界 相位延迟 

单位：中国科学院地质与地球物理研究所，北京100029; 中国地震局地球物理研究所，北京100081