摘要：起伏地表下地震波场模拟有助于解释主动源和被动源地震探测中穿过山脉和盆地的测线所获得的资料．然而传统的有限差分法处理起伏的自由边界比较困难，为了克服这一困难，我们将笛卡尔坐标系的各向异性介质弹性波方程和自由边界条件变换到曲线坐标系中，采用一种稳定的、显式的二阶精度有限差分方法离散（曲线坐标系）VTI介质中的弹性波方程；对地表自由边界条件处理时采用了一种修饰的差分算子来计算弹性波方程中的混合导数项在自由边界上的法向导数．兰姆问题的解析解与本文的数值解对比结果表明该方法可以有效地处理自由地表边界条件．模拟实例表明：起伏地表对地震波场有重要影响，各向异性导致弹性波波前形状复杂且具有明显的方向性．

关键词：起伏地表 各向异性 曲线网格 波场模拟 有限差分 

单位：中国科学院地质与地球物理研究所 岩石圈演化国家重点实验室 北京100029 中国科学院研究生院 北京100049