三维非均匀介质中弹性波传播的数值模拟

提出了一种三维非均匀介质中弹性波传播数值模拟的方法 ,文中称为三维格子法 .该算法是二维格子法 (一种二维非均匀介质中P SV波传播的数值模拟算法 )向三维非均匀介质情况的推广 .在空间离散上该文方法与有限元方法类似 ,容许根据连续体的形状和介质分界面任意剖分网格 ,且自然满足自由表面边界条件 .不同于常规有限差分法在各个节点上满足动力学微分方程 ,该算法通过满足各节点周围格子的整体平衡 (积分平衡方程 )来对问题进行求解 .三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当 .算例表明 ,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射 .