起伏地表弹性波传播有限差分法数值模拟

？有限差分是进行地震波传播数值模拟的最常用方法，但该方法处理起伏的自由边界比较困难。为此，通过对不同地形起伏情况下自由边界的具体分析，将整个二维空间离散点划分为24类，对每一类自由边界处的网格点选择了合理的表现方式，实现了起伏地表自由边界条件的数值化。该方法可以模拟出地表起伏情况下弹性波复杂的传播现象，为进行起伏地表地震波传播规律研究、山地地震勘探野外观测系统设计、山地地震勘探干扰波分析和识别、以及静校正研究提供了正演模拟工具。模拟实例表明，地形起伏引起面波、体波等地震波型之间的相互转化，产生了大量的散射p波、散射s波和散射面波，尤其是由沿地表传播的强能量面波，在地表起伏及近地表物性突变处产生了大量的强能量散射面波，同时也产生了相对较弱的散射p波和散射s波，这是造成山地地震资料信噪比低的主要原因。