标量波方程广义有限差分正演模拟

常规有限差分法地震波场正演模拟受限于固定的网格步长,不可避免地会出现网格剖分与实际速度界面不一致的情形,进而带来起伏界面处的阶梯状绕射以及反射波旅行时不准确等问题.广义有限差分法是一种无网格方法,它基于泰勒函数展开和加权最小二乘拟合,将微分方程中未知参数的偏导数表示为相邻节点函数值的线性组合,可根据不同地质体模型建立适用的场节点分布形式,克服了传统有限差分法对网格的依赖性.文中采用沿层布置节点的非均匀剖分方法,使生成的节点与起伏的地表或边界贴合.针对不同模型的测试结果表明,广义有限差分法能有效解决虚假反射、反射波旅行时不准确等问题,且具有稳定性.