乘匹配算法将两者进行子波匹配，得到电缆鬼波压制算子I;
[0022]B)将得到的电缆鬼波压制算子I作为滤波器对未压制鬼波的地震数据进行滤波，得到仅压制了电缆鬼波的地震数据；
[0023]C)将不含鬼波的地震子波作为期望输出，将仅含震源鬼波的地震子波作为最小二乘匹配算法的输入，通过最小二乘匹配算法对两者进行子波匹配，得到震源鬼波压制算子2 ；
[0024]D)将得到的震源鬼波压制算子2作为滤波器对步骤B)中得到的仅压制了电缆鬼波的地震数据进行滤波，得到压制了震源鬼波的地震数据，即得到同时压制了电缆鬼波和震源鬼波的地震数据；
[0025]如图2所示，方法二的具体过程为:
[0026]a)将地震资料信息中不含鬼波的地震子波作为期望输出，含电缆鬼波和震源鬼波的地震子波作为最小二乘匹配算法的输入，采用最小二乘法对两者进行子波匹配，得到电缆鬼波和震源鬼波压制算子；
[0027]b)将得到的电缆鬼波和震源鬼波压制算子作为滤波器对未压制鬼波的地震数据进行滤波，得到同时压制了电缆鬼波和震源鬼波的地震数据；
[0028]3)输出同时压制了电缆鬼波和震源鬼波的地震数据，用于地震资料处理的后续环
-K-T。
[0029]下面通过具体实施例对本发明基于子波处理的鬼波压制方法的效果进行详细说明:
[0030]实施例1:
[0031]如图3所示，首先，模拟反射系数模型，其结果如(a)所示。在未考虑子波衰减的情况下，将带震源鬼波和电缆鬼波的子波与反射系数模型合成地震记录，其结果如(b)所示，其结果显示合成地震记录中出现多轴多相位现象，不能够真实表示反射系数模型，因此需要对鬼波进行压制，使用本发明对合成地震记录(b)进行鬼波压制，结果如(d)所示，其结果显示鬼波压制后同相轴相位变单一，能够清楚表示反射系数。
[0032]地震波在地下传播过程中高频信号会逐渐衰减，在考虑子波衰减情况下，将带震源鬼波和电缆鬼波的子波与反射系数模型合成地震记录，其结果如(C)所示，对合成的地震记录采用本发明进行鬼波压制，结果如(e)所示，其结果显示，在考虑了子波衰减这种接近地震波真实传播规律的情况下，本发明依然能够获得较好的鬼波压制效果，压制鬼波后的地震记录仍能反映地下反射系数，表明了本发明对于鬼波压制的有效性。
[0033]实施例2:
[0034]本发明在实际地震资料中的应用效果如图4和图5所示。
[0035]如图4所示，(a)中箭头指出未进行鬼波压制的海底反射波，海底反射波表现为白-黑-白三个相位的反射特征，其中第一个白色相位是一次反射波，黑色相位为震源鬼波和电缆鬼波的叠置，第二个白色相位是电缆鬼波。使用本发明对(a)进行压制鬼波处理，得到如(b)所示的仅压制电缆鬼波的结果，以及如(c)所示的压制电缆鬼波和震源鬼波的结果。从图4可知，随着鬼波被压制，箭头所指的海底反射波由多相位逐渐变为单一相位。
[0036]如图5所示，使用本发明压制实际地震资料中的鬼波，并放大深层的反射同相轴，从(a)、(b)、(c)中可以看出，随着鬼波被压制，箭头指示的反射层位由多相位变为单一相位，剖面波组特征变清楚，鬼波压制后的剖面更接近地下的反射系数。如图6所示，随着电缆鬼波和震源鬼波的压制，频谱的低频端和高频段都得到了拓展，频带变宽，这对后期的地震资料解释、沉积相分析、储层预测、综合地质评价等研宄意义重大。
[0037]现阶段的鬼波压制方法有基于频率域L2优化的波场分离技术，该技术纯数据驱动，主要用于压制海上拖缆数据中的电缆鬼波。其压制电缆鬼波的剖面(如图7所示)相比于图5中(a)显示的未压制电缆鬼波的剖面，箭头指示的反射层的相位变少，说明电缆鬼波得到了压制，相比于图5中(b)可以发现，基于频率域L2优化的波场分离技术达到了本发明对电缆鬼波压制的效果，但是现有的基于频率域L2优化的波场分离技术仅能压制电缆鬼波，对震源鬼波却难以压制。
[0038]上述各实施例仅用于说明本发明，其中各个步骤等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
【主权项】
1.一种基于子波处理的鬼波压制方法，包括以下步骤: 1)输入海上地震资料信息；其中，地震资料信息中包含未压制鬼波的地震数据和三种地震子波，三种地震子波包括含电缆鬼波和震源鬼波的地震子波、仅含震源鬼波的地震子波和不含鬼波的地震子波； 2)通过地震子波对未压制鬼波的地震数据进行鬼波压制； 3)输出同时压制了电缆鬼波和震源鬼波的地震数据。
2.如权利要求1所述的一种基于子波处理的鬼波压制方法，其特征在于:所述步骤2)通过地震子波对未压制鬼波的地震数据进行鬼波压制，包括以下步骤: A)将地震资料信息中仅含震源鬼波的地震子波作为期望输出，含电缆鬼波和震源鬼波的地震子波作为最小二乘匹配算法的输入，通过最小二乘匹配算法将两者进行子波匹配，得到电缆鬼波压制算子； B)将得到的电缆鬼波压制算子作为滤波器对未压制鬼波的地震数据进行滤波，得到仅压制了电缆鬼波的地震数据； C)将不含鬼波的地震子波作为期望输出，将仅含震源鬼波的地震子波作为最小二乘匹配算法的输入，通过最小二乘匹配算法对两者进行子波匹配，得到震源鬼波压制算子； D)将得到的震源鬼波压制算子作为滤波器对步骤B)中得到的仅压制了电缆鬼波的地震数据进行滤波，得到压制了震源鬼波的地震数据，即得到同时压制了电缆鬼波和震源鬼波的地震数据。
3.如权利要求1所述的一种基于子波处理的鬼波压制方法，其特征在于:所述步骤2)通过地震子波对未压制鬼波的地震数据进行鬼波压制，包括以下步骤: a)将地震资料信息中不含鬼波的地震子波作为期望输出，含电缆鬼波和震源鬼波的地震子波作为最小二乘匹配算法的输入，采用最小二乘法对两者进行子波匹配，得到电缆鬼波和震源鬼波压制算子； b)将得到的电缆鬼波和震源鬼波压制算子作为滤波器对未压制鬼波的地震数据进行滤波，得到同时压制了电缆鬼波和震源鬼波的地震数据。
【专利摘要】本发明涉及一种基于子波处理的鬼波压制方法，包括以下步骤：1)输入海上地震资料信息；其中，地震资料信息中包含未压制鬼波的地震数据和三种地震子波，三种地震子波包括含电缆鬼波和震源鬼波的地震子波、仅含震源鬼波的地震子波和不含鬼波的地震子波；2)通过地震子波对未压制鬼波的地震数据进行鬼波压制；3)输出同时压制了电缆鬼波和震源鬼波的地震数据。本发明采用了子波匹配的算法压制鬼波，从地震资料处理入手压制鬼波，而不是从采集入手压制鬼波，因此处理过程易于实现，节约成本。本发明能够压制海上地震资料中的震源鬼波和电缆鬼波，可以广泛应用于海洋地震勘探过程中。
【IPC分类】G01V1-28, G01V1-37
【公开号】CN104536045
【申请号】CN201510023338
【发明人】陈宝书, 李松康, 汪小将, 仝中飞, 张连群, 陶杰, 杨小椿, 薛冬
【申请人】中国海洋石油总公司, 中海油研究总院
【公开日】2015年4月22日
【申请日】2015年1月16日