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

作者:贾宗锋; 吴国忱; 李青阳; 杨凌云; 吴悠
来源:石油地球物理勘探, 2022, 57(01): 101-8.
DOI:10.13810/j.cnki.issn.1000-7210.2022.01.011

摘要

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