一种应用于弹性波正演的三阶Higdon阻尼吸收边界方法

专利申请号：CN202310789261.5

公 开 号：CN116822297A

发 明 人：秦志亮 解闯 申恒广 于盛齐 马本俊 任强 刘雪芹 兰文剑 马林伟 孟怡君 

代 理 人：袁崇建

代理机构：长沙准星专利代理事务所(普通合伙)

专利类型：发明专利

申 请 日：20230929

公 开 日：20230630

专利主分类号：G06F30/23

关 键 词：正演模拟 中心波 弹性波方程 人工边界 阻尼吸收 场区域 弹性波 入射波 二阶 三阶 模型实验结果 吸收边界条件 边界反射 边界条件 差分格式 地震勘探 高稳定性 交错网格 模拟区域 边界处 地震波 高阶数 偶数阶 扰动 截断 一阶 吸收 地下 传播 应用 

摘      要：本发明涉及一种应用于弹性波正演的三阶Higdon阻尼吸收边界方法，属于地震勘探正演模拟领域，具体包括以下步骤：在采用有限差分进行一阶应力‑速度弹性波方程正演模拟时，将模拟区域划分为中心波场区域和人工边界区域，首先在中心波场区域中，采用时间二阶、空间任意偶数阶精度的交错网格有限差分格式对弹性波方程进行离散计算，然后在人工边界区域中，采用高阶数、高稳定性的三阶Higdon阻尼吸收边界条件对入射波进行吸收。模型实验结果表明，相比常规二阶Higdon吸收边界条件，本发明方法能够在截断边界处更好地吸收各个角度的入射波，极大地减少了边界反射波对中心波场的扰动影响，其使弹性波正演模拟的结果更接近于地震波在地下传播的真实情况。