基于PML吸收边界的三维声波数值模拟方法

摘要

本发明公开了一种基于PML吸收边界的三维声波数值模拟方法，其主要包括如下步骤：s1、在内部计算区域应用规则网格，仅在规则网格整数节点上定义波场压力；在PML吸收区域应用交错网格，在交错网格整数节点上定义波场压力，在交错网格半数节点上定义质点的速度分量和介质密度；s2、在内部计算区域采用二阶双曲型偏微分波动方程更新波场压力值；在PML吸收区域采用一阶应力速度‑方程计算质点的速度分量和波场压力值。本发明方法既减少了内存需求量，扩大了可模拟的模型规模，同时又使得PML边界处理更为简洁。本发明方法尤其适用于油气勘探领域中的地震波场数值模拟或逆时偏移。