[发明专利]地震数据反演无拉伸动校正方法、装置、电子设备及介质

说明书

公开了一种地震数据反演无拉伸动校正方法、装置、电子设备及介质。该方法可以包括：参考道贝叶斯反演定位反射系数位置；以反射系数位置为约束，确定稀疏反射系数反演的目标函数；针对稀疏反射系数反演的目标函数通过最小二乘法求解，计算动校正后的反射系数；通过动校正后的反射系数计算动校正后的地震数据。本发明获得无拉伸畸变的动校正结果，达到提高角道集资料质量，更易于后续处理解释工作的目的。

技术领域

本发明涉及石油、天然气能源的地球物理勘探处理领域，更具体地，涉及一种地震数据反演无拉伸动校正方法、装置、电子设备及介质。

背景技术

在地震资料处理中，动校正的目的是消除炮间距对反射波旅行时的影响，校平共深度点反射波时距曲线的轨迹，增强利用叠加技术压制干扰的能力，减小叠加过程引起的反射波同相轴畸变。常规数字动校正在远偏移距处存在波形拉伸的现象，其主频降低，影响高分辨率叠加和高精度AVO分析，严重时需要切除，而切除会使得浅层远道信息缺失，覆盖次数降低。动校正子波拉伸会使得主频向低频移动，子波拉伸是AVO分析中的不利因素。

由于动校正子波拉伸现象会对AVO分析产生影响，国内外的地球物理工作者研究了很多消除子波拉伸及其影响的方法，例如：

通过对数据块整体移动进行动校正，避免了动校正量不同而导致的拉伸畸变，但在炮检距较大时易出现过校正和欠校正。先自动检测反射波垂直反射时间和速度、后计算待校正炮检距的反射时间并将该时间处波形整体移动的方法，较好地消除了拉伸畸变，但该方法对速度分析的精度要求非常高。在分析动校正拉伸对属性影响的基础上对截距和梯度进行校正，但这种方法有一定的局限性。在叠加过程中应用抛物Radon变换消除动校正拉伸现象，但需在Radon域找到合适路径，不利于工业应用。无拉伸动校正方法需找准每个同相轴的初至时间，很难选取合适的时窗和时间因子。频谱代换无拉伸动校正方法理论上可以消除时差，但实际中不同同相轴在进行相位替换后会相互影响。多次迭代的方法，每次迭代时进行部分动校正并利用反褶积来修正拉伸的影响，减小拉伸畸变，但每次迭代过程中采用的是常规动校正方法，动校正结果仍然会引入拉伸畸变。这些方法不能从根本上消除动校正对子波拉伸的影响，而且实际处理中难以应用，稳定性欠佳。

因此，有必要开发一种基于横向几何特征约束的地震数据反演无拉伸动校正方法、装置、电子设备及介质。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种地震数据反演无拉伸动校正方法、装置、电子设备及介质，其能够获得无拉伸畸变的动校正结果，达到提高角道集资料质量，更易于后续处理解释工作的目的。

第一方面，本公开实施例提供了一种地震数据反演无拉伸动校正方法，包括：

参考道贝叶斯反演定位反射系数位置；

以所述反射系数位置为约束，确定稀疏反射系数反演的目标函数；

针对所述稀疏反射系数反演的目标函数通过最小二乘法求解，计算动校正后的反射系数；

通过所述动校正后的反射系数计算动校正后的地震数据。

优选地，参考道贝叶斯反演定位反射系数位置包括：

建立基于褶积模型的地震数据矩阵；

计算反射系数的后验概率密度分布；

计算超参数，根据所述超参数与所述反射系数的后验概率密度分布，确定反射系数位置。

优选地，反射系数的后验概率密度分布为：