地震勘探记录中机器噪声压制方法及系统技术方案

公开了一种地震勘探记录中机器噪声压制方法及系统。该方法可以包括：通过主通道输入含有机器噪音的地震信号；确定滤波器的阶数与初始滤波值，进而获得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)＝[ω0(t),ω1(t),...,ωL-1(t)]T(4)其中，ω(t)表示t时刻对应的滤波值，ω0(t)、ω1(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)其中，e(t)表示t时刻对应的消噪输出值，d(t)表示主通道输入的地震信号。根据本专利技术的另一方面，提出了一种地震勘探记录中机器噪声压制系统，可以包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：通过主通道输入含有机器噪音的地震信号；确定滤波器的阶数与初始滤波值，进而获得t时刻对应的滤波值；通过参考通道输入机器噪音，根据所述机器噪音与所述t时刻对应的滤波值，获得t时刻对应的滤波器输出值；根据所述t时刻对应的滤波器输出值与所述主通道输入的地震信号，获得t时刻对应的消噪输出值，进而获得每个时刻对应的消噪输出值。优选地，还包括：根据所述t时刻对应的滤波值、滤波器的收敛步长、所述t时刻对应的消噪输出值与所述机器噪音，获得t+1时刻对应的滤波值，进而获得t+1时刻对应的消噪输出值。优选地，所述t时刻对应的滤波值为：ω(t)＝[ω0(t),ω1(t),...,ωL-1(t)]T(4)其中，ω(t)表示t时刻对应的滤波值，ω0(t)、ω1(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)其中，e(t)表示t时刻对应的消噪输出值，d(t)表示主通道输入的地震信号。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1a和图1b分别示出了采集的地震记录和采集地震数据之前同一位置记录的背景噪声的示意图。图2示出了根据本专利技术的地震勘探记录中机器噪声压制方法的步骤的流程图。图3a、图3b、图3c、图3d和图3e分别示出了有效地震信号、机器噪声、主通道输入含有机器噪音的地震信号、参考通道输入的机器噪音以及根据本专利技术的消噪输出值的示意图。图4示出了根据本专利技术的一个实施例的放炮前记录的背景噪声的示意图。图5示出了根据本专利技术的一个实施例的与图4同一位置采集的地震信号的示意图。图6示出了根据本专利技术的一个实施例的消噪输出值的示意图。图7示出了根据本专利技术的一个实施例的剔除掉的机器噪声的示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图2示出了根据本专利技术的地震勘探记录中机器噪声压制方法的步骤的流程图。在该实施例中，根据本专利技术的地震勘探记录中机器噪声压制方法可以包括：步骤101，通过主通道输入含有机器噪音的地震信号。具体地，通过两个通道输入地震信号，主通道输入为被机器噪声污染的地震信号d(t)，含有有效信号s(t)和机器噪声n1(t)。参考通道输入为在采集地震信号之前，在同一位置对机器噪声做的记录，记为x(t)＝n2(t)。由于机器噪声有一定的规律性，而n1(t)和n2(t)是在不同时刻，同一位置记录的机器噪声，所以，n1(t)和n2(t)具有良好相关性。因此，想要求得的输出信号为公式(1)：e＝s+n1-y(1)其中，e表示输出信号，s表示有效信号，n1表示机器噪声，y表示滤波器输出值，则均方误差为公式(2)：E{e2}＝E{(s+n1-y)2}＝E{s2}+2E(s(n1-y))+E{(n1-y)2}(2)由于s与噪声n1和n2是无关的，所以公式(2)中第二项为0，则公式(2)变形为公式(3)：E{e2}＝E{(s+n1-y)2}＝E{s2}+E{(n1-y)2}＝E{s2}+E{(n1-ωTn2)2}(3)从公式(3)容易看出，若调节滤波器ω，使得E{e2}取值最小，必然有n1＝ωTn2。也就是说，在均方误差最小本文档来自技高网...

【技术保护点】
1.一种地震勘探记录中机器噪声压制方法，包括：通过主通道输入含有机器噪音的地震信号；确定滤波器的阶数与初始滤波值，进而获得t时刻对应的滤波值；通过参考通道输入机器噪音，根据所述机器噪音与所述t时刻对应的滤波值，获得t时刻对应的滤波器输出值；根据所述t时刻对应的滤波器输出值与所述主通道输入的地震信号，获得t时刻对应的消噪输出值，进而获得每个时刻对应的消噪输出值。

【技术特征摘要】
1.一种地震勘探记录中机器噪声压制方法，包括：通过主通道输入含有机器噪音的地震信号；确定滤波器的阶数与初始滤波值，进而获得t时刻对应的滤波值；通过参考通道输入机器噪音，根据所述机器噪音与所述t时刻对应的滤波值，获得t时刻对应的滤波器输出值；根据所述t时刻对应的滤波器输出值与所述主通道输入的地震信号，获得t时刻对应的消噪输出值，进而获得每个时刻对应的消噪输出值。2.根据权利要求1所述的地震勘探记录中机器噪声压制方法，其中，还包括：根据所述t时刻对应的滤波值、滤波器的收敛步长、所述t时刻对应的消噪输出值与所述机器噪音，获得t+1时刻对应的滤波值，进而获得t+1时刻对应的消噪输出值。3.根据权利要求1所述的地震勘探记录中机器噪声压制方法，其中，所述t时刻对应的滤波值为：ω(t)＝[ω0(t),ω1(t),...,ωL-1(t)]T(4)其中，ω(t)表示t时刻对应的滤波值，ω0(t)、ω1(t)、ωL-1(t)均表示滤波器的权重系数，L表示滤波器的阶数。4.根据权利要求3所述的地震勘探记录中机器噪声压制方法，其中，所述t时刻对应的滤波器输出值为：y(t)＝ωT(t)*X(t)(5)其中，y(t)表示t时刻对应的滤波器输出值，X(t)表示参考通道输入的机器噪音。5.根据权利要求4所述的地震勘探记录中机器噪声压制方法，其中，所述t时刻对应的消噪输出值为：e(t)＝d(t)-y(t)(6)其中，e(t)表示t时刻对应的消噪输出值，d(t)表示主通道输入的地震信号。6.一种地震勘探记录中机器噪声压制系统，其特征在于，该系统包...

【专利技术属性】
技术研发人员：魏巍，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京,11