[发明专利]菲涅尔带约束的浅水多次波衰减方法

说明书

一种菲涅尔带约束的浅水多次波衰减方法，包括以下步骤：一，对原始的地震数据进行预处理；二，重构近偏移距信息；三，对格林函数进行构建；四，地震数据与格林函数褶积；其采用浅水多次波预测的公式为：五，菲涅尔带约束叠加；六，使用自适应匹配进行相减，其公式为：p(t)＝d(t)‑f(t)*m(t)；七，得到浅水多次波衰减后的地震数据。本发明不仅能够有效地进行浅水多次波衰减，消除地震资料中的陷频问题，拓展地震资料频带；而且，还为后续构造解释和储层研究提供重要的资料基础，对海上油气田勘探开发具有重要意义。

技术领域

本发明属于油气勘探地震资料处理领域，尤其涉及菲涅尔带约束的浅水多次波衰减方法。

背景技术

在进行地震勘探时，地震波在不同地层之间、薄层之间、裂缝内部与裂缝之间会发生多次散射或折射后被检波器接收，其被统称为多次波或多次散射。如果多次波处理不当，很容易混淆、干扰一次有效反射波的能量，进而在很大程度上，降低地震资料的偏移成像质量，甚至导致对地下构造的错误认识，直接影响后续的井位部署和勘探成功率。且由于在浅水情况下，在水层中完整传播后，被检波器接收到的浅水多次波的旅行时，与一次反射波的旅行时的时差较小，波形相互混叠，严重影响了互层结构下，对储层发育模式的识别。因此，如何有效地解决浅水多次波的衰减问题是地震资料处理中的重点、难点和热点。

当前，有关深水情况下的多次波预测与衰减，已有相对成熟的理论和方法，且得到了大规模的应用，而浅水情况下的多次波预测和衰减技术还不完善，仍然面临着巨大的挑战。

目前，浅水多次波的衰减方法主要分为滤波类方法和波动方程类方法两种。其中，波动方程类方法又包括：逆散射级数法、反馈迭代法、波场延拓法和格林函数法等方法；逆散射级数法和反馈迭代方法是基于数据驱动，但由于在浅水条件下，其无法区分直达波和海底反射波，因此，不仅导致近偏移距处反射波估算困难，很难获得规则采样的全波场数据；而且，浅水多次波的多阶性使得计算过程需要多次迭代，容易引入振幅误差，存在伤害有效波的风险。而基于波场延拓的确定性浅水多次波衰减方法，其是通过对炮记录进行完整的上、下行波场延拓来模拟出多次波；但它是基于瑞利积分对波场延拓的假设，因此，对于弯曲反射界面的波场延拓，其结果不是完全正确的，而且，在上行延拓和下行延拓中，都没有很好地考虑海底反射系数的影响。

发明内容

本发明目的在于提供一种菲涅尔带约束的浅水多次波衰减方法,以解决在浅水条件下，自由表面多次波衰减的要求，为后续油气田勘探开发提供重要的资料基础的技术问题。

为实现上述目的，本发明的菲涅尔带约束的浅水多次波衰减方法的具体技术方案如下：

一种菲涅尔带约束的浅水多次波衰减方法，包括以下步骤：

第一步，对原始的地震数据进行预处理；

第二步，重构近偏移距信息；

第三步，对格林函数进行构建；

第四步，地震数据与格林函数褶积

将第二步得到的地震数据和第三步得到的时间域格林函数变换到频率域；根据相同的偏移距信息，在频率域将每一个地震道与格林函数进行对应相乘，从而得到每一个地震道经过一次完整的水层传播后的结果；且采用浅水多次波预测如下公式：

第五步，菲涅尔带约束叠加

将选取的地震数据进行叠加得到预测的浅水多次波，以提升目标地震道浅水多次波预测精度和信噪比；

第六步，自适应相减

将原始地震数据与第五步得到的浅水多次波预测结果进行自适应匹配相减，从而，最终实现浅水多次波的衰减；使用自适应匹配相减的方法公式为：p(t)＝d(t)-f(t)*m(t) (2)；

第七步，得到浅水多次波衰减后的地震数据。