公开(公告)号：CN116822297A

公开(公告)日：2023-09-29

申请号：CN202310789261.5

申请日：2023-06-30

Applicant: 哈尔滨工程大学

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

IPC: G06F30/23 ,  G06F17/10 ,  G06F119/14

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

公开(公告)号：CN116660979A

公开(公告)日：2023-08-29

申请号：CN202310649157.6

申请日：2023-06-02

Applicant: 哈尔滨工程大学

Inventor： 秦志亮 ,  解闯 ,  申恒广 ,  任强 ,  马本俊 ,  于盛齐 ,  刘雪芹 ,  兰文剑 ,  马林伟

IPC: G01V1/28 ,  G01V1/30

Abstract: 本发明涉及一种基于Kaiser时窗积分的OBN资料全波形反演方法，属于地震波形反演领域，所述方法为首先采用Kaiser时窗对震源子波和OBN观测地震记录分别进行积分，然后采用声波方程有限差分正演模拟技术得到正传波场和模拟地震记录，将该模拟地震记录与积分后的观测地震记录的残差记录平方和作为反演的目标函数，并采用声波方程有限差分逆时延拓技术得到反传波场，通过正传波场与反传波场的互相关得到目标函数的梯度，最后对速度模型进行多次迭代更新，从而实现OBN资料全波形反演。本发明方法可有效压制积分处理中的高频旁瓣；且在缺失2Hz以下低频信息的情况下，可为后续的常规全波形反演提供可靠的初始模型，有效避免反演算法陷入局部极值。

公开(公告)号：CN119105076A

公开(公告)日：2024-12-10

申请号：CN202411253836.2

申请日：2024-09-09

Applicant: 哈尔滨工程大学

Inventor： 解闯 ,  孟怡君 ,  蔡文伯 ,  秦志亮 ,  马本俊 ,  申恒广 ,  毕成浩 ,  于盛齐 ,  刘雪芹 ,  缪可言

IPC: G01V1/28 ,  G01V1/38 ,  G01V1/30 ,  G06F30/23 ,  G06F111/10

Abstract: 本发明涉及一种应用于低频舰船地震波数值模拟的联合吸收边界方法，属于舰船地震波数值模拟领域，所述方法在数值模拟过程中，首先将模拟区域划分为中心波场区域和人工边界区域；在中心波场区域内，将雷克子波作为舰船噪声声源，基于雷克子波表达式模拟出声源扰动，并采用交错网格有限差分法进行舰船地震波数值模拟；舰船地震波波动方程由二维一阶速度‑应力弹性波方程来描述；在中心波场区域内，对弹性波方程进行交错网格有限差分离散，采用时间2阶、空间2N阶精度的展开形式；在人工边界区域内，添加L层完全匹配层作为吸收边界，在PML层内采用高斯型衰减函数；在外边界处联合使用二阶Higdon吸收边界。本发明方法可在不增加存储消耗的情况下，有效提高整体边界吸收效果。

公开(公告)号：CN116660979B

公开(公告)日：2024-01-30

申请号：CN202310649157.6

申请日：2023-06-02

Applicant: 哈尔滨工程大学

Inventor： 秦志亮 ,  解闯 ,  申恒广 ,  任强 ,  马本俊 ,  于盛齐 ,  刘雪芹 ,  兰文剑 ,  马林伟

IPC: G01V1/28 ,  G01V1/30

Abstract: 本发明涉及一种基于Kaiser时窗积分的OBN资料全波形反演方法，属于地震波形反演领域，所述方法为首先采用Kaiser时窗对震源子波和OBN观测地震记录分别进行积分，然后采用声波方程有限差分正演模拟技术得到正传波场和模拟地震记录，将该模拟地震记录与积分后的观测地震记录的残差记录平方和作为反演的目标函数，并采用声波方程有限差分逆时延拓技术得到反传波场，通过正传波场与反传波场的互相关得到目标函数的梯度，最后对速度模型进行多次迭代更新，从而实现OBN资料全波形反演。本发明方法可有效压制积分处理中的高频旁瓣；且在缺失2Hz以下低频信息的情况下，可为后续的常规全波形反演提供可靠的初始模型，有效避免反演算法陷入局部极值。

公开(公告)号：CN116822297B

公开(公告)日：2024-01-16

申请号：CN202310789261.5

申请日：2023-06-30

Applicant: 哈尔滨工程大学

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

IPC: G06F30/23 ,  G06F17/10 ,  G06F119/14

Abstract: 本发明涉及一种应用于弹性波正演的三阶Higdon阻尼吸收边界方法，属于地震勘探正演模拟领域，具体包括以下步骤：在采用有限差分进行一阶应力‑速度弹性波方程正演模拟时，将模拟区域划分为中心波场区域和人工边界区域，首先在中心波场区域中，采用时间二阶、空间任意偶数阶精度的交错网格有限差分格式对弹性波方程进行离散计算，然后在人工边界区域中，采用高阶数、高稳定性的三阶Higdon阻尼吸收边界条件对入射波进行吸收。模型实验结果表明，相比常规二阶Higdon吸收边界条件，本发明方法能够在截断边界处更好地吸收各个角度的入射波，极大地减少了边界反射波对中心波场的扰动影响，其使弹性波正演模拟的结果更接近于地震波(56)对比文件杨慧珠,张友生,陶果.井眼条件下弹性波传播问题的三维有限差分数值模拟.地球物理学进展.2003,(第02期),全文.