本发明属于海洋地震勘探领域，特别是涉及一种能够压制船干扰、增强有效波的自适应船干扰压制方法。
背景技术：
海洋地震勘探中往往受到外源的干扰，对后期地震资料处理以及解释产生巨大影响，其中船干扰是常见的外源干扰之一，有效去除船干扰的影响，可以大大提高后续地震资料的处理精度，进而提高地震解释的可靠性，具有十分重要的意义。传统的船干扰压制技术是将含有船干扰的炮集数据转换到其他排列的道集中，在该道集中，船干扰表现为具有非相干性的噪声，可以采用随机噪声压制的方式进行压制,但是该类方法对同相轴顶点附近的噪声能量压制效果较差，且容易损伤有效信号。另一类常用的方法是通过人机交互或外源扫描等方式来确定船干扰的轨迹，然后利用线性滤波方法进行压制,但是此类方法往往需要人工干预，或者需要扫描较大的空间范围以及相关先验信息来确定干扰船的位置，通常需要耗费大量的时间和人力、物力。
技术实现要素：
针对上述问题，本发明的目的是提供一种自适应船干扰压制方法，在高效去除船干扰的同时，有效地保护有效信号。
为实现上述目的，本发明采取以下技术方案：一种自适应船干扰压制方法，其特征在于包括以下步骤：1)基于含船干扰的地震数据，采用互相关的方式确定干扰船的位置参数；2)计算干扰船到各检波点的距离，并将最小距离对应的检波器为参考道，利用互相关计算其他道与该参考道之间的时间延迟，并将该时间延迟进行校正，得到校正后的各地震道集；3)对校正后的各地震道集做滤波处理，进行船干扰压制，继而进行反校正，得到最终船干扰压制后的地震道集。
所述步骤1)中，采用互相关的方式确定干扰船的位置参数的方法为：首先从含船干扰的地震数据中，随机选取三道信号；然后取其中1道为参考道，通过互相关的方式确定另外两道信号与参考道信号之间的时间延迟，建立方程组；最后对建立的方程组进行求解，得到干扰船的位置参数。
所述干扰船的位置参数的计算公式为：
式中，v为海水的速度，(x,y)为干扰船的位置参数，(x0,y0)为参考道r0对应检波点的位置参数，(xi,yi)和(xj,yj)分别为另外两道信号ri,rj对应检波点的位置参数。
所述步骤2)中，干扰船到第k个检波点的距离为：
式中，k为检波点序号，且k＝1,2,...,n，n为检波器的总个数。
本发明由于采取以上技术方案，其具有以下优点：1、本发明采用互相关的方式确定干扰船的位置，无需平缆假设，且不需要假设干扰船在缆线范围内，具有较好的适应性。2、本发明利用互相关的方式确定各地震道集的时间延迟，进行校正、滤波、反校正，容易理解，且快速有效，可以在去除船干扰的同时，保护有效波。因而，本发明可以广泛应用于海洋地震勘探中的船干扰压制领域。
附图说明
图1是本发明自适应船干扰压制方法流程图；
图2是船干扰的示意图；
图3(a)～图3(d)是本发明的处理效果图，其中图3(a)为不含船干扰的地震记录，图3(b)是含有模拟船干扰的地震记录，图3(c)为压制掉的船干扰，图3(d)为船干扰压制后的地震记录。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示，本发明提供一种自适应船干扰压制方法，包括以下步骤：
(1)在含船干扰的地震数据中，随机选取三道信号r0,ri,rj，并取其中1道为参考道r0，通过互相关的方式确定另外两道信号ri,rj与参考道信号r0之间的时间延迟δti,δtj，建立方程组并求解，得到干扰船的位置。
如图2所示，含船干扰地震记录中，假设参考道r0对应检波点的位置参数为(x0,y0)，另外两道ri,rj对应检波点的位置参数分别为(xi,yi)和(xj,yj)，采用互相关的方式对干扰船的位置参数进行计算的公式为：
式中，v为海水的速度，约为1500m/s，(x,y)为干扰船的位置参数。
(2)计算干扰船到各检波点的距离，并将最小距离对应的检波器为参考道，利用互相关计算其他道与该参考道之间的时间延迟，并将该时间延迟进行校正，即将船干扰对应的曲线拉平。其中，对时间延迟进行校正的方法为常规方法，本发明在此不再赘述。
干扰船到第k个检波点的距离为：
式中，k为检波点序号，且k＝1,2,...,n，n为检波器的总个数。记最小距离d*对应的检波器为参考道r*。
(3)对校正后的道集做滤波处理，进行船干扰压制，继而进行反校正，得到最终船干扰压制后的地震道集。其中，对校正后的道集做滤波处理以及反校正的方法均为常规方法，本发明在此不再赘述。
实施例一
如图3(a)～3(d)所示，本实施例以模拟地震资料为例，对本发明方法进行验证。其中，图3(a)为不含船干扰的地震记录，图3(b)为含有模拟船干扰的地震记录，可以看出，船干扰从浅到深展现为一串平行的曲线，可以利用地震信号互相关方法确定各道与参考道之间的时间延迟，进而进行校正以及船干扰的压制。图3(c)为压制掉的船干扰，可以看出，压制掉的船干扰中无有效波的痕迹。图3(d)为船干扰压制后的地震记录，可以看出船干扰得到了有效压制，有效波得到了有效增强，可以提高后续数据成像的精度，有利于地震资料的合理解释。本发明能够有效地去除船干扰的影响，突出有效波，且可以有效地保护有效信息，可有效提高后续资料处理与解释的精度。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。