一种基于子波处理的鬼波压制方法
【技术领域】
[0001]本发明涉及海洋地震勘探领域，特别是关于一种基于子波处理的鬼波压制方法。
【背景技术】
[0002]鬼波是海洋地震勘探中始终存在的难题，长期没有得到根本的解决，严重制约了地震勘探技术的进步。按照鬼波形成过程的不同，将鬼波分为震源鬼波和电缆鬼波。将气枪震源沉放于海水中的一定深度激发，有一部分能量直接向上传播到海面，然后再向下反射向海底深部传播，遇到地下界面后再反射回来被电缆检波器所接收，此波恰好尾随在正常一次反射波之后构成干扰，称作震源鬼波。而与此对应，震源激发后能量向下反射，遇到地下界面后再反射回来至海面，再由海面反射回来被电缆检波器所接收，此波也同样在正常一次反射波之后构成干扰，称作电缆鬼波。鬼波的存在造成地震资料中出现陷频，压缩地震资料的频带，同时导致地震资料的多轴、多相位现象，从而无法准确地获得反射系数序列，因而严重影响构造解释、储层反演、油气检测和开发地震的应用效果。
[0003]目前压制鬼波的主流思路是从采集入手，例如:海底电缆采集(OBC)、上下缆采集、上下源采集、双传感器采集、斜缆采集等。这些采集技术配合其相对应的高精度处理技术都能在一定程度上压制鬼波，如针对OBC资料和双传感器资料的水陆检合并技术；针对上下缆上下源资料的上下缆、上下源合并技术；针对斜缆资料的镜像偏移技术、联合反褶积技术等。然而，通过采集手段压制鬼波不仅费用昂贵，而且配套的高精度处理过程复杂、难度较大。与此同时这些技术大部分都只能压制电缆鬼波，对于震源鬼波并没有很好的压制效果。成本高、处理技术难度大、难以压制震源鬼波是此类采集技术压制鬼波的弊端。因此，研宄一种同时压制海上地震资料中震源鬼波和电缆鬼波的处理方法十分必要。
【发明内容】
[0004]针对上述问题，本发明的目的是提供一种能够同时压制海上地震资料中的震源鬼波和电缆鬼波的基于子波处理的鬼波压制方法。
[0005]为实现上述目的，本发明采取以下技术方案:一种基于子波处理的鬼波压制方法，包括以下步骤:1)输入海上地震资料信息；其中，地震资料信息中包含未压制鬼波的地震数据和三种地震子波，三种地震子波包括含电缆鬼波和震源鬼波的地震子波、仅含震源鬼波的地震子波和不含鬼波的地震子波；2)通过地震子波对未压制鬼波的地震数据进行鬼波压制；3)输出同时压制了电缆鬼波和震源鬼波的地震数据。
[0006]所述步骤2)通过地震子波对未压制鬼波的地震数据进行鬼波压制，包括以下步骤:A)将地震资料信息中仅含震源鬼波的地震子波作为期望输出，含电缆鬼波和震源鬼波的地震子波作为最小二乘匹配算法的输入，通过最小二乘匹配算法将两者进行子波匹配，得到电缆鬼波压制算子；B)将得到的电缆鬼波压制算子作为滤波器对未压制鬼波的地震数据进行滤波，得到仅压制了电缆鬼波的地震数据；C)将不含鬼波的地震子波作为期望输出，将仅含震源鬼波的地震子波作为最小二乘匹配算法的输入，通过最小二乘匹配算法对两者进行子波匹配，得到震源鬼波压制算子；D)将得到的震源鬼波压制算子作为滤波器对步骤B)中得到的仅压制了电缆鬼波的地震数据进行滤波，得到压制了震源鬼波的地震数据，即得到同时压制了电缆鬼波和震源鬼波的地震数据。
[0007]所述步骤2)通过地震子波对未压制鬼波的地震数据进行鬼波压制，包括以下步骤:a)将地震资料信息中不含鬼波的地震子波作为期望输出，含电缆鬼波和震源鬼波的地震子波作为最小二乘匹配算法的输入，采用最小二乘法对两者进行子波匹配，得到电缆鬼波和震源鬼波压制算子；b)将得到的电缆鬼波和震源鬼波压制算子作为滤波器对未压制鬼波的地震数据进行滤波，得到同时压制了电缆鬼波和震源鬼波的地震数据。
[0008]本发明由于采取以上技术方案，其具有以下优点:1、本发明由于从地震资料处理入手压制鬼波，而不是从采集入手压制鬼波，因此节约成本。2、本发明由于采用了含震源鬼波和电缆鬼波的地震子波、仅含震源鬼波的地震子波、不含鬼波的地震子波之间的匹配处理对未压制鬼波的地震数据进行鬼波压制，因此既可以压制电缆鬼波，又可以压制震源鬼波，进而恢复地下真实反射信号，获得准确的反射系数序列。3、本发明由于可以先利用含震源鬼波和电缆鬼波的地震子波与仅含震源鬼波的地震子波之间的匹配压制电缆鬼波，再利用仅含震源鬼波和不含鬼波的地震子波之间的匹配压制震源鬼波，以实现所有鬼波的压制；也可以直接利用含震源鬼波和电缆鬼波的地震子波与不含鬼波的地震子波之间的匹配实现所有鬼波的压制，因此本发明实现方式灵活。4、本发明由于采用了子波匹配的算法压制鬼波，因此处理过程易于实现。5、本发明计算效率高，相比于现有技术中采用最优化算法的方法，本发明计算速度快60倍以上，因此本发明更加适用于对实际地震资料的处理中。本发明可以广泛应用于海洋地震勘探过程中。
【附图说明】
[0009]图1是本发明方法一的流程示意图；
[0010]图2是本发明方法二的流程示意图；
[0011]图3是本发明实施例1模拟的合成地震记录中的压制鬼波的效果示意图；其中，(a)是模拟的反射系数模型示意图，(b)是未考虑子波衰减的带震源鬼波和电缆鬼波的子波与反射系数模型合成的地震记录示意图，(C)是子波衰减情况下带震源鬼波和电缆鬼波的子波与反射系数模型合成的地震记录示意图，(d)是对图(b)中合成的地震记录使用本发明压制鬼波的结果示意图，(e)是对图(C)中合成的地震记录使用本发明压制鬼波的结果示意图；
[0012]图4是本发明实施例2实际地震资料鬼波压制的海底反射局部放大示意图；其中，(a)是未进行鬼波压制的叠加剖面示意图，(b)是仅压制电缆鬼波后的叠加剖面示意图，(C)是压制电缆鬼波和震源鬼波后的叠加剖面示意图；
[0013]图5是本发明实施例2实际地震资料鬼波压制的深层反射局部放大示意图；其中，(a)是未进行鬼波压制的叠加剖面示意图，(b)是仅压制电缆鬼波后的叠加剖面示意图，(C)是压制电缆鬼波和震源鬼波后的叠加剖面示意图；
[0014]图6是本发明实施例2鬼波压制前后的频谱对比示意图；其中，黑色虚线表示未压制鬼波的频谱，灰色实线表示仅压制电缆鬼波后的频谱，黑色实线表示压制震源鬼波和电缆鬼波后的频谱；
[0015]图7是现有的基于频率域L2优化的波场分离技术实际地震资料鬼波压制的深层反射局部放大示意图。
【具体实施方式】
[0016]下面结合附图和实施例对本发明进行详细的描述。
[0017]如图1、图2所示，本发明的基于子波处理的鬼波压制方法，包括以下步骤:
[0018]I)输入海上地震资料信息；其中，地震资料信息中包含未压制鬼波的地震数据和三种地震子波，三种地震子波包括含电缆鬼波和震源鬼波的地震子波、仅含震源鬼波的地震子波和不含鬼波的地震子波；
[0019]2)通过地震子波对未压制鬼波的地震数据进行鬼波压制，具体的压制方法包括以下两种:
[0020]如图1所示，方法一的具体过程为:
[0021]A)将地震资料信息中仅含震源鬼波的地震子波作为期望输出，含电缆鬼波和震源鬼波的地震子波作为现有的最小二乘匹配算法的输入，通过最小二