[发明专利]可控震源独立同步激发地震数据邻炮干扰压制方法和装置

说明书

技术领域

本发明涉及地球物理勘探技术领域，特别涉及一种可控震源独立同步激发（Independent Simultaneous Sweeping，ISS）地震数据邻炮干扰压制方法和装置。

背景技术

可控震源ISS地震数据在采集的过程中是通过多组可控震源独立作业、同步激发，并通过连续观测的方式将多组可控震源同步激发产生的地震波场记录到一个母记录中的方式获得的。这种获取可控震源ISS地震数据的方法极大地压缩了地震数据的数据量、提高了获取地震数据的效率。

然而，可控震源ISS地震数据在采集中可控震源同步激发产生的地震波场之间相互会产生干扰、混叠，从而导致地震数据的准确性不高。

发明内容

本发明提供了一种可控震源ISS地震数据邻炮干扰压制方法和装置，以达到提高地震数据准确性的目的。

本发明实施例提供了一种可控震源ISS地震数据邻炮干扰压制方法，包括：在可控震源ISS地震数据中分选出非共炮点道集地震数据；对分选出的非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据；从邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据。

在一个实施例中，对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：在时域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，在频域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，在时域与频域对分选出的非共炮点道集地震数据进行矢量中值滤波。

在一个实施例中，在时域对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：沿时间轴按预定长度将非共炮点道集地震数据分成段；将各段非共炮点道集地震数据作为矢量，按以下公式按地震数据道空间分布对各段非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据：

Σi=1n||Xm,k-Xi,k||p≤Σi=1n||Xj,k-Xi,k||p]]>

其中，X_m,k∈{X_i,k|i=1,...,n}，j=1,...,n，X_m,k为输出的结果，k为对非共炮点道集地震数据分段后各段的顺序号，X_i,k为分段后各段的非共炮点道集地震数据，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波的道数，为矢量Lp范数。

在一个实施例中，所述预定长度包括：一个子波或者一个波形对应时间轴上的长度。

在一个实施例中，在频域对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：按照以下公式对分选出的非共炮点道集地震数据进行傅里叶变换，得到频域的非共炮点道集地震数据：

X(ω)_i=FFT[X(t)_i]

其中，X(t)_i为时域的非共炮点道集地震数据，FFT[]为一维的傅里叶变换，ω为角频率，X(ω)_i为傅里叶变换后得到的频域的非共炮点道集地震数据；