一种梯度分解的全波形反演方法及其设备

    公开(公告)号:CN117590462A

    公开(公告)日:2024-02-23

    申请号:CN202311454254.6

    申请日:2023-11-03

    IPC分类号: G01V1/28 G01V1/30

    摘要: 本发明公开了一种梯度分解的全波形反演方法及其设备,涉及油气勘探地球物理技术领域,包括以下步骤:步骤S1:数据采集,将设备安装在指定的位置进行记录;步骤S2:数据预处理,对采集到的数据进行进一步的处理;步骤S3:建立初始模型,选择或构建一个初步的地下模型,并将该模型将作为全波形反演的起点;步骤S4:进行正演模拟,利用步骤S3中所选择的地下模型进行模拟以及计算;步骤S5:计算梯度,根据得到的梯度信息来更新地下模型从而来不断的改进模型。本发明通过输入初始速度模型和观测地震记录,同时设置反演参数,根据所述初始速度模型、观测地震记录和反演参数生成全波形反演梯度,利用Poynting矢量分离方法对全波形反演梯度进行梯度分解。

    一种面向构造控制型花岗岩体高波数阻抗地震成像方法

    公开(公告)号:CN117590463A

    公开(公告)日:2024-02-23

    申请号:CN202311454264.X

    申请日:2023-11-03

    IPC分类号: G01V1/28 G01V1/34

    摘要: 本发明公开了一种面向构造控制型花岗岩体高波数阻抗地震成像方法,涉及地质勘探技术领域,包括以下步骤:步骤S1:获取输入参数,获取的参数包括P波速度参数场、观测炮记录;步骤S1包括:步骤S101:计划制定,根据研究区域和目标制定相关的计划;步骤S102:设备布置,部署相关的地震仪器以及传感器;步骤S103:数据采集,通过步骤S102中部署的设备对相关数据进行采集并记录。本发明能够解决传统互相关逆时偏移成像方法产生的低波数伪影问题,进而产生更加清晰的成像结果,开发面向构造控制型深部高温花岗岩体高波数阻抗地震成像技术,为复杂构造区的地震资料处理和后续地震解释工作提供高精度成像的理论依据和技术支持。

    一种降低低波数段频散误差的高精度正演模拟方法

    公开(公告)号:CN117521448A

    公开(公告)日:2024-02-06

    申请号:CN202311454267.3

    申请日:2023-11-03

    摘要: 本发明涉及勘探地球物理学技术领域,公开了一种降低低波数段频散误差的高精度正演模拟方法。该降低低波数段频散误差的高精度正演模拟方法,包括以下步骤,S101,对二阶空间偏导的差分离散格式带入平面波解得到空间域频散关系,S103,通过空间域频散关系绝对误差建立等波纹误差函数,通过雷米兹交换算法计算波纹误差函数下的优化差分系数,和该差分系数的SMW,S105,通过空间域频散关系构建一范数目标函数。本发明通过构建具有约束条件的一范数拉格朗日对偶问题,对低波数段误差进行约束,抑制了优化差分系数的误差,相较于其他空间域差分优化算法,本方法可以在低波数取得较高的模拟精度。

    一种基于优化差分系数的正演模拟方法及其设备

    公开(公告)号:CN115034103B

    公开(公告)日:2024-08-16

    申请号:CN202210506447.0

    申请日:2022-05-11

    摘要: 本说明书实施例公开了一种基于优化差分系数的正演模拟方法及其设备。方法包括:获取输入参数;采用预设矩阵确定优化差分系数;根据所述输入参数和所述优化差分系数对方程进行求解,生成正演模拟结果。通过最小二乘方法最小化不同CFL数下频散误差所求得的差分算子,相比于传统的时空域差分算子,最小化不同CFL数下的频散关系误差方法天然地解决了时空域差分算子对CFL数敏感的问题,实际应用也证明了最小二乘方法最小化不同CFL数下的频散误差求得的差分算子可以很好的适用于不同的CFL数;模拟过程中不存在CFL数对于数值模拟的影响;在CFL数变化时,模拟过程中波场仍可稳定传播,精度更高。

    一种多参数点扩散函数反褶积弹性波最小二乘偏移方法

    公开(公告)号:CN118409357A

    公开(公告)日:2024-07-30

    申请号:CN202410522549.0

    申请日:2024-04-28

    IPC分类号: G01V1/36 G01V1/28

    摘要: 本发明公开了一种多参数点扩散函数反褶积弹性波最小二乘偏移方法,涉及勘探地球物理技术领域,包括如下步骤:基于弹性VTI波动方程计算点扩散函数;利用点扩散函数分量构建滤波器分量;应用波数域高维反褶积,得到成像结果。本发明基于弹性VTI波动方程Born正演和伴随偏移算子,计算4个弹性点扩散函数分量:PP,PS,SP和SS;利用这些点扩散函数分量构建滤波器,用来近似逆Hessian矩阵;利用高斯单元分割,将传统偏移成像结果划分为与单个点扩散函数大小相同的局部成像结果,在波数域进行局部传统偏移成像结果与多参数点扩散函数的高维反褶积,重建高精度地下构造信息。本发明方法成像精度高,且计算成本低。

    面向海底潜山高陡构造的声弹耦合逆时偏移地震成像方法

    公开(公告)号:CN118409356A

    公开(公告)日:2024-07-30

    申请号:CN202410522545.2

    申请日:2024-04-28

    IPC分类号: G01V1/36 G01V1/28

    摘要: 本发明公开了面向海底潜山高陡构造的声弹耦合逆时偏移地震成像方法,涉及勘探地球物理技术领域,包括:采用基于偏应力表征的声弹耦合波动方程对OBN数据正演模拟;应用亥姆霍兹分解方法进行P/S波分离;针对P、S与压力源波场进行上下左右行波分离;采用相同的分解方法对接收波场进行行P/S波和上下左右行波的分离,并设计与倾角相关的权重因子;分别对上下行波和左右行波应用零滞后互相关成像条件,得到两组成像结果并将左右行波成像结果乘以权重因子与上下行波成像结果进行叠加,获得最终的叠加成像结果。该方法可以在保持高分辨率的基础上,显著提高对高陡构造成像的清晰度和准确性同时有效避免低波数伪影的干扰,确保了高质量的成像结果。

    一种面向地震数据的三维自适应双域逐级去噪方法

    公开(公告)号:CN118409355A

    公开(公告)日:2024-07-30

    申请号:CN202410522542.9

    申请日:2024-04-28

    IPC分类号: G01V1/36 G01V1/28

    摘要: 本发明公开了一种面上地震数据的三维自适应双域逐级去噪方法,涉及勘探地球物理领域。本方法利用拉普拉斯掩码估计噪声方差场,得到局部梯度下的噪声方差近似值,并加入鲁棒性噪声估计算子调节局部方差大小,提高对非稳态信号的适应性;并且本方法将传统二维方法拓展到了三维,充分考虑了三维地震数据的空间相干性,提高了地震数据空间连续性,三维梯度场高斯核的滤波函数相较二维对噪声的估计更加准确。该三维自适应双域渐进去噪方法可以对强吸收地震数据体进行自适应降噪,更加适应于存在非平稳性的低信噪比地震数据中,拥有良好的实际应用前景。

    一种黑三角噪音时频域自适应压制方法

    公开(公告)号:CN117724160A

    公开(公告)日:2024-03-19

    申请号:CN202311483079.3

    申请日:2023-11-08

    IPC分类号: G01V1/30 G06F17/14 G01V1/28

    摘要: 本发明公开了一种黑三角噪音时频域自适应压制方法,涉及地球物理技术领域,包括:获取输入数据;计算共炮点道集的线性正值归一化;对归一化后的共炮点道集振幅升序排列;对道集横向每一道进行二次函数拟合;对拟合后的道集沿着横向对每一道计算曲率;对曲率剖面提取最大值,得到自适应阈值百分比曲线;对计算曲率的数据沿着水平方向进行高斯函数拟合,并沿横向计算标准差;对自适应阈值百分比曲线修正,迭代直到取得最优自适应阈值百分比曲线;对原始数据变换到时频域,得到时频数据体;对时频域数据进行自适应软阈值;反变换到时间域,得到黑三角噪音,最后得到结果。本发明能够显著提升信噪比,在不损害有效信号的情况下最大程度衰减噪音。

    一种基于空间域伪差分算子的VTI介质纯qP波正演模拟方法

    公开(公告)号:CN116663351A

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

    申请号:CN202310505436.5

    申请日:2023-05-06

    IPC分类号: G06F30/23 G06F111/10

    摘要: 本发明公开了一种基于空间域伪差分算子的VTI介质纯qP波正演模拟方法,涉及勘探地球物理技术领域,包括:将纯qP波方程化为含有修正伪差分算子的椭圆各向异性方程;对修正伪差分算子近似,得到卷积核库;将椭圆各向异性方程转换为空间域形式;获取输入数据;在震源位置加载震源函数子波;根据非椭圆各向异性区域的各向异性参数索引卷积核库,将卷积核与应力卷积;使用交错网格有限差分计算修正后应力的横向、垂向二阶导数;波场更新,并根据波场快照输出间隔与采样点间隔对波场进行存储;重复上述步骤到数值模拟完成;输出波场快照与炮记录。本发明可对qP波进行高精度模拟,避免傅里叶变换,且适用于区域分解的并行计算,计算成本低,应用前景好。

    一种深层高温花岗岩体的全波形速度反演成像方法

    公开(公告)号:CN116626759A

    公开(公告)日:2023-08-22

    申请号:CN202310631187.4

    申请日:2023-05-31

    IPC分类号: G01V1/30

    摘要: 本说明书实施例公开了一种深层高温花岗岩体的全波形速度反演成像方法。通过获取输入参数,所述输入参数包括初始纵波速度模型、观测系统信息和观测数据dobs;根据所述输入参数,计算地震子波,根据所述地震子波计算合成数据dsyn,并使合成数据与观测数据相减得到残差数据dres;反传来自接收器位置的残差数据来生成反传波场ures;根据所述反传波场和正传波场的零延迟互相关得到梯度g,并对所述梯度采用如下方式进行优化,生成优化后的梯度根据所述优化后的梯度进行反演以最小化残差数据,生成速度模型,从而实现增加反演的低波数成分来缓解传统全波形反演方法的“周波跳跃”问题,使得模型能更快速且正确的收敛。