CN101630019B - 一种校正预测的海底多次波大时差的方法 - Google Patents

Abstract

本发明是地震资料处理中校正预测的海底多次波大时差的方法。用多次波预测算子生成预测的海底反射波和多次波模型道集，自动拾取海底反射波的旅行时进行中值滤波，计算原始炮集和预测的海底反射波的旅行时差，对多次波模型道集做时差校正，应用时差校正后的多次波模型道集与原始数据进行数据匹配，从原始数据中减去匹配后的多次波道集，输出压制多次波后的地震剖面。本发明是在常规数据匹配之前先对大的时差进行校正的方法，因而能够适应存在拖缆漂移的采集数据，提高自适应减去法的效率，改善压制多次波的效果。

Description

一种校正预测的海底多次波大时差的方法

技术领域

本发明涉及地球物理勘探技术，属于地震资料处理中数据匹配的技术范畴，是一种校正预测的海底多次波大时差的方法。

背景技术

在海上地震勘探中，由于作为自由界面的海水面和硬的海底面都是很强的反射界面，其中海水面的反射系数可达0.9，因此在地震记录中不可避免的会记录到海底和海水面产生的各种各样的多次波。多次波的存在会干扰人们对有效波的识别，影响后续的速度分析、偏移成像、解释的真实性和可靠性，因此必须加以压制。

压制多次波的方法一般分为两大类：一类是基于有效波和多次波之间差异的滤波方法，简称为滤波方法；另一类是基于波动方程的预测减去法，即首先利用原始地震数据，根据波动原理，预测出多次波，然后对预测的多次波与原始数据中的多次波进行匹配，从原始数据中减去匹配后的多次波就可以得到一次有效波。

在基于波动方程的预测减去方法中，根据波动原理利用原始地震数据预测的多次波，与原始数据中的多次波之间会存在振幅、相位和旅行时差异，因此需要对预测的多次波以原始数据中的多次波为目标进行匹配，然后用原始炮集减去匹配后的预测的多次波就达到消除或压制多次波的目的。

数据匹配一般要求匹配滤波算子的长度较短，因为长的匹配滤波算子要求目标数据，包括预测的多次波和原始地震数据的长度也要长，而当目标数据较长时，就有可能包含不期望参与数据匹配的信号和干扰波，进而影响匹配效果。然而在实际海上拖缆地震数据采集中，受海流的影响，拖缆漂移是不可避免的，因此一方面由拖缆漂移的地震数据得到的海底高程会存在误差，另一方面当把拖缆漂移的地震道投影到测线时，也不能真实地反映地震道的位置，这样就有可能导致预测的海底多次波与实际资料中的海底多次波存在较大的时差，这种时差客观上需要很长的匹配滤波算子来进行处理。因此，为了减少匹配算子的长度，提高匹配的效率与精度，在数据匹配之前，必须消除预测的海底多次波大的时差。

发明内容

本发明目的在于提供一种能够适应存在拖缆漂移的采集数据，提高自适应减去法的效率，改善压制多次波的效果的校正预测的海底多次波大时差的方法。

本发明的具体实施步骤如下：

1)采集炮集地震记录；

2)用多次波预测算子对输入原始炮集地震记录进行运算，生成预测的海底反射波和多次波模型道集；

步骤2)中所述的输入的原始炮集记录中要保留直达波；

3)利用原始炮集和预测的海底反射波和多次波模型道集自动拾取海底反射波的旅行时；

步骤3)中拾取海底反射波的旅行时是拾取海底反射波同相轴的初至对应的时间；

4)对拾取的海底反射波旅行时进行中值滤波：

步骤4)中所述的中值滤波是指，将某点为中点的窗口中的所有点的值由小到大进行排序，中间值作为该点的输出值。

步骤4)中所述的中值滤波窗口的长度取7-11个采样点。

5)计算原始炮集和预测的海底反射波的旅行时差，对多次波模型道集做时差校正；

步骤5)所述的计算原始炮集和预测的海底反射波的旅行时差是利用原始炮集中确定的海底反射波旅行时减去预测的海底反射波的旅行时。

步骤5)所述的时差校正逐道进行，按照多次波模型道与原始数据的海底反射时差整道上移或下移，向上移动时地震道尾部所缺部分充零，向下移动时地震道头部所缺部分充零。

6)应用时差校正后的多次波模型道集与原始数据进行数据匹配；

7)从原始数据中减去匹配后的多次波道集；

8)输出压制多次波后的地震剖面。

本发明不是采用增加匹配滤波算子长度的方法来校正匹配数据间的大的时差，而是在常规数据匹配之前先对大的时差进行校正的方法，因而能够适应存在拖缆漂移的采集数据，提高自适应减去法的效率，改善压制多次波的效果。

附图说明

图1原始炮集(左)及预测的海底反射波和海底多次波道集(右)。

图2原始炮集(左)及预测的海底反射波和海底多次波道集(右)中自动拾取的海底反射波的旅行时曲线。

图3原始炮集中的海底反射波及预测的海底反射波的旅行时之差。

图4对预测的海底反射波和海底多次波道集进行大的时差校正后的结果。

图5对预测的海底多次波进行大的时差校正和数据匹配后的结果。

图6原始炮集。

图7是原始炮集及大的时差校正后多次波自适应减去后的结果。

具体实施方案

在预测海底或自由表面多次波时，人们往往要把直达波切除，因为多次波预测算子会把直达波转换成海底反射波，这样在多次波自适应减去时，海底反射波也会被消除。

本发明的不切除原始炮集中的直达波，这样多次波预测算子就会把直达波转换成海底反射波，即预测的海底反射波；自动拾取原始炮集中和预测的海底反射波的旅行时，计算两者之间的时差；设原始炮集中和预测的海底多次波之间的时差与原始炮集中和预测的海底反射波之间的时差是相同的，这样就可以利用原始炮集中和预测的海底反射波旅行时之差，对预测的海底多次波进行大的时差校正。

本发明的具体实施方式为：

以下是本发明具体实例：

1)取得原始炮集地震记录。

2)按照步骤(2)对原始地震数据行处理得到预测的海底反射波和多次波模型道集

图1原始炮集及预测的海底反射波和海底多次波道集。可以看出原始炮集中的海底反射波和预测的海底反射波有较大的旅行时误差。

3)按照步骤3)-步骤4)分别在原始炮集、预测的海底反射波和海底多次波道集上拾取海底反射波的旅行时。

图2是分别在原始炮集及预测的海底反射波和海底多次波道集中自动拾取的海底反射波的旅行时曲线。

4)按照步骤5)计算原始炮集中的海底反射波和预测的海底反射波的旅行时误差，然后对预测的海底反射波和海底多次波道集进行大的时差校正。

图3是原始炮集的海底反射波及预测的海底反射波的旅行时之差，可以看出误差比较大，一般在-100.0ms-100.0ms之间。

图4是对预测的海底反射波和海底多次波道集进行大的时差校正后的结果。可以看出，经过校正，预测的多次波与原始数据中的多次波已经不存在大的旅行时误差。

5)按照步骤6)-步骤8)对时差校正后的预测多次波进行匹配，然后从原始数据中减去匹配后的多次波就可以得到有效波。

图5是对预测的海底多次波进行大的时差校正和数据匹配后的结果。

图6是原始炮集。

图7是原始炮集及大的时差校正后多次波自适应减去后的结果。可以看出，多次波的能量得到了很好的压制。

Claims (3)

1.一种校正预测的海底多次波大时差的方法，其特征在于实施步骤如下：

1)采集炮集地震记录；

2)用多次波预测算子对输入原始炮集地震记录进行通常的运算，生成预测的海底反射波和多次波模型道集，所述的输入原始炮集地震记录中要保留直达波；

3)利用原始炮集和预测的海底反射波和多次波模型道集自动拾取海底反射波的旅行时；其中所述拾取海底反射波的旅行时是拾取海底反射波同相轴的初至对应的时间；

4)对拾取的海底反射波旅行时进行中值滤波：

5)计算原始炮集和预测的海底反射波的旅行时差，对多次波模型道集做时差校正；

所述的计算原始炮集和预测的海底反射波的旅行时差是利用原始炮集中确定的海底反射波旅行时减去预测的海底反射波的旅行时；

所述的时差校正逐道进行，按照多次波模型道集与原始数据的海底反射时差整道上移或下移，向上移动时地震道尾部所缺部分充零，向下移动时地震道头部所缺部分充零；

6)应用时差校正后的多次波模型道集与原始数据进行数据匹配；

7)从原始数据中减去匹配后的多次波模型道集；

8)输出压制多次波后的地震剖面。

2.根据权利要求1所述的校正预测的海底多次波大时差的方法，其特征在于步骤4)中所述的中值滤波是指，将某点为中点的窗口中的所有点的值由小到大进行排序，中间值作为该点的输出值。

3.根据权利要求1所述的校正预测的海底多次波大时差的方法，其特征在于步骤4)中所述的中值滤波窗口的长度取7-11个采样点。

Images