公开(公告)号：CN104765948A

公开(公告)日：2015-07-08

申请号：CN201510096974.9

申请日：2015-03-05

Applicant: 山东科技大学

Inventor： 林年添 ,  王守进 ,  丁仁伟 ,  张建彬 ,  赵传伟 ,  魏乾乾 ,  杨修超 ,  刘向超 ,  文博 ,  张栋

IPC: G06F19/00

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

公开(公告)号：CN104765948B

公开(公告)日：2019-01-22

申请号：CN201510096974.9

申请日：2015-03-05

Applicant: 山东科技大学

Inventor： 林年添 ,  王守进 ,  丁仁伟 ,  张建彬 ,  赵传伟 ,  魏乾乾 ,  杨修超 ,  刘向超 ,  文博 ,  张栋

IPC: G06F17/10

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