[发明专利]地震勘探记录中机器噪声压制方法及系统

说明书

公开了一种地震勘探记录中机器噪声压制方法及系统。该方法可以包括：通过主通道输入含有机器噪音的地震信号；确定滤波器的阶数与初始滤波值，进而获得t时刻对应的滤波值；通过参考通道输入机器噪音，根据机器噪音与t时刻对应的滤波值，获得t时刻对应的滤波器输出值；根据t时刻对应的滤波器输出值与主通道输入的地震信号，获得t时刻对应的消噪输出值，进而获得每个时刻对应的消噪输出值。本发明通过获取主通道与参考通道的地震信号，并进行计算消噪，可以有效抵消地震记录中的机器噪声，大幅提高资料的信噪比，降低有效信号的失真。

技术领域

本发明涉及地球物理勘探领域，更具体地，涉及一种地震勘探记录中机器噪声压制方法及系统。

背景技术

伊尔马滋(地震资料分析[M].北京：石油工业出版社.2006.)指出在f-k(频率-波数)滤波方法是目前压制地震记录中相干噪声的主流技术。f-k滤波的方法主要利用噪声与有效信号之间的视速度(同相轴倾角)差异。当相干噪声与有效信号具有不同的视速度时，在f-k域噪声与信号将表现为不同的倾角，因此，容易进行噪声的压制。而当噪声与有效信号的视速度比较接近时，在f-k域将无法准确区分噪声与信号，因此，该方法将不再适用(李振春、张军华.地震数据处理方法[M].东营：中国石油大学出版社.2004.)。

在工厂、老油区等特殊环境进行地震数据采集时，经常会受到周围机械设备的强噪声影响，致使地震记录噪声干扰严重，资料品质低下，因此,必须设法压制机器噪声。图1a和图1b分别示出了采集的地震记录和采集地震数据之前同一位置记录的背景噪声的示意图，从圈中部分可以看出机器噪声近似规则，是一种相干噪声，并且视速度与有效信号近似。这一特点使得主流的f-k滤波无法用于压制这种相干噪声。因此，有必要开发一种地震勘探记录中机器噪声压制方法及系统。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种地震勘探记录中机器噪声压制方法及系统，其能够通过获取主通道与参考通道的地震信号，并进行计算消噪，可以有效抵消地震记录中的机器噪声，大幅提高资料的信噪比，降低有效信号的失真。

根据本发明的一方面，提出了一种地震勘探记录中机器噪声压制方法。所述方法可以包括：通过主通道输入含有机器噪音的地震信号；确定滤波器的阶数与初始滤波值，进而获得t时刻对应的滤波值；通过参考通道输入机器噪音，根据所述机器噪音与所述t时刻对应的滤波值，获得t时刻对应的滤波器输出值；根据所述t时刻对应的滤波器输出值与所述主通道输入的地震信号，获得t时刻对应的消噪输出值，进而获得每个时刻对应的消噪输出值。

优选地，还包括：根据所述t时刻对应的滤波值、滤波器的收敛步长、所述t时刻对应的消噪输出值与所述机器噪音，获得t+1时刻对应的滤波值，进而获得t+1时刻对应的消噪输出值。

优选地，所述t时刻对应的滤波值为：

ω(t)＝[ω₀(t),ω₁(t),...,ω_L-1(t)]^T (4)

其中，ω(t)表示t时刻对应的滤波值，ω₀(t)、ω₁(t)、ω_L-1(t)均表示滤波器的权重系数，L表示滤波器的阶数。

优选地，所述t时刻对应的滤波器输出值为：

y(t)＝ω^T(t)*X(t) (5)

其中，y(t)表示t时刻对应的滤波器输出值，X(t)表示参考通道输入的机器噪音。

优选地，所述t时刻对应的消噪输出值为：

e(t)＝d(t)-y(t) (6)