一种地震数据中规则干扰的压制方法及装置

摘要

本发明提供一种地震数据中规则干扰的压制方法及装置。所述方法包括：将地震道按照地理位置顺序关系排列形成待处理地震数据炮集记录；对存在规则干扰覆盖的地震数据区域进行规则噪声频谱分析获取规则干扰的估算主频；分离出规则干扰形成规则干扰数据记录，并转换得到FK域地震数据记录；能量扫描确定地震数据中能量最强的前K个倾角；在规则干扰数据记录中，若在所述倾角中连续统计到N个正的振幅值，则判断该倾角上存在规则干扰，然后对该倾角方向进行滤波形成新的干扰数据记录；用待处理地震数据炮集记录减去新的干扰数据记录。利用本发明实施例，可以有效减小地震数据规则干扰去噪时的有效信号损伤，提高噪声压制处理后的有效信号的保真度。

权利要求

1.一种地震数据中规则干扰的压制方法，其特征在于，所述方法包括：将采集的地震数据中的地震道按照地理位置顺序关系进行排列，形成待处理地震数据炮集记录；

对所述待处理地震数据炮集记录中判断为存在规则干扰覆盖的地震数据区域进行规则干扰噪声频谱分析，获取所述规则干扰的估算主频；

通过设计的滤波器从所述待处理地震数据炮集记录中分离出频率为所述估算主频的规则干扰，形成规则干扰数据记录；对所述规则干扰数据记录进行FK变换，得到FK域地震数据记录；

在预先设定的规则干扰倾角范围内对所述FK域地震数据记录进行能量扫描，确定地震数据记录中能量最强的前K个倾角，K≥1；

在所述规则干扰数据记录中分别对所述K个倾角进行地震数据振幅统计；若在所述倾角中连续统计到N个正的振幅值，则判断所述倾角对应的地震数据记录存在规则干扰，N≥

2；

在所述待处理地震数据炮集记录中，沿判断存在规则干扰的倾角方向进行滤波，提取所述存在规则干扰的倾角方向上的规则干扰波，形成新的干扰数据记录；

用所述待处理地震数据炮集记录减去所述新的干扰数据记录，形成规则干扰压制后的地震数据。

2.如权利要求1所述的一种地震数据中规则干扰的压制方法，其特征在于，所述方法还包括：利用预设频段的带通滤波器对所述新的干扰数据记录进行滤波，获取滤波后干扰数据记录；所述预设频段包括获取的用户设置的认为存在有效信号的频段范围；

相应的，所述用所述待处理地震数据炮集记录减去所述新的干扰数据记录包括：用所述待处理地震数据炮集记录减去所述滤波后干扰数据记录。

3.如权利要求1或2所述的一种地震数据中规则干扰的压制方法，其特征在于，所述N的取值范围为：5≤N≤32。

4.如权利要求3所述的一种地震数据中规则干扰的压制方法，其特征在于，所述N的取值为：16。

5.如权利要求1或2所述的一种地震数据中规则干扰的压制方法，其特征在于，所述K的取值范围为：1≤K≤100。

6.如权利要求5所述的一种地震数据中规则干扰的压制方法，其特征在于，所述K的取值为：20。

7.一种地震数据中规则干扰的压制装置，其特征在于，所述装置包括：排序处理模块，用于将采集的地震数据中的地震道按照地理位置顺序关系进行排列，形成待处理地震数据炮集记录；

干扰主频估算模块，用于对所述待处理地震数据炮集记录中判断为存在规则干扰覆盖的地震数据区域进行规则干扰噪声频谱分析，获取所述规则干扰的估算主频；

干扰分离及变换模块，用于通过设计的滤波器从所述待处理地震数据炮集记录中分离出频率为所述估算主频的规则干扰，形成规则干扰数据记录，以及对所述规则干扰数据记录进行FK变换，得到FK域地震数据记录；

扫描模块，用于获取设置的K值，以及，在预先设定的规则干扰倾角范围内对所述FK域地震数据记录进行能量扫描，确定地震数据记录中能量最强的前K个倾角，K≥1；

振幅统计模块，用于获取设置的N值，以及，在所述规则干扰数据记录中分别对所述K个倾角进行地震数据振幅统计；若在所述倾角中连续统计到N个正的振幅值，则判断所述倾角对应的地震数据记录存在规则干扰，N≥2；

倾角滤波模块，用于在所述待处理地震数据炮集记录中，沿判断存在规则干扰的倾角方向进行滤波，提取所述存在规则干扰的倾角方向上的规则干扰波，形成新的干扰数据记录；

干扰处理模块，用于用所述待处理地震数据炮集记录减去所述新的干扰数据记录，形成规则干扰压制后的地震数据。

8.如权利要求7所述的一种地震数据中规则干扰的压制装置，其特征在于，所述装置还包括：规则干扰数据滤波模块，用于利用预设频段的带通滤波器对所述新的干扰数据记录进行滤波，获取滤波后干扰数据记录；所述预设频段包括获取的用户设置的认为存在有效信号的频段范围；

相应的，所述干扰处理模块中所述用所述待处理地震数据炮集记录减去所述新的干扰数据记录包括：用所述待处理地震数据炮集记录减去所述滤波后干扰数据记录。

9.如权利要求7或8所述一种地震数据中规则干扰的压制装置，其特征在于，所述N的取值范围为：5≤N≤32。

10.如权利要求9所述的一种地震数据中规则干扰的压制装置，其特征在于，所述N的取值为：16。

11.如权利要求7或8所述的一种地震数据中规则干扰的压制装置，其特征在于，所述K的取值范围为：1≤K≤100。

12.如权利要求11所述的一种地震数据中规则干扰的压制装置，其特征在于，所述K的取值为：20。

说明书

一种地震数据中规则干扰的压制方法及装置

技术领域

[0001] 本发明涉及地球物理勘探数据处理技术，是一种地震数据处理中通过压制规则干扰来达到提高地震数据信噪比和分辨率的方法。

背景技术

[0002] 地震数据勘探采集时，野外布置的检波器接收到的地震信号通常含有多种扰动波，其中包括可用于解决地质任务的有效波,也包括影响地震数据成像的干扰波。

[0003] 一般的，根据干扰波的分布特征可以将干扰波分成规则干扰波和随机干扰波。规则干扰波又可以被称为相干噪声，这种波在地震数据时间和空间分布上通常具有一定的规律性。而随机干扰波通常由随机噪声引起，在地震数据上时间和空间上的分布没有任何规律。这两类干扰波通常会在地震数据处理中影响地震数据成像质量,为了提高地震数据信噪比和分辨率，经常需要对这两类干扰波进行压制，而压制效果也将直接影响地震数据成像质量。因此对于低信噪比地震数据处理噪声压制，尤其是对相干噪声的压制已经是地震数据处理中的一个重要环节。在地震勘探领域，常用的相干去噪声压制方法包括比如FK域规则噪声压制技术、KL变换规则噪声压制技术、t-p域规则噪声压技术、T-X域和F-X域规则噪声压制技术等，而这些规则干扰压制技术会在一定程度上削弱有效信号能量，尤其是会对低频有效信号造成较大的损伤，甚至有时会破坏有效信号，导致给后期处理的地震数据不够准确，甚至带来无法估量的损失，不适应当前正在发展的低频勘探技术。

[0004] 现有技术中的规则噪声的压制技术通常会对地震数据中的有效信号尤其是对低频有效信息造成较大的损伤，导致给后期处理的地震数据不够准确，甚至带来无法估量的损失，不适应当前正在发展的低频勘探技术。

发明内容

[0005] 本发明目的在于提供一种地震数据中规则干扰的压制方法及装置，可以有效减小地震数据规则干扰去噪时的有效信号损伤，在压制规则干扰的同时尽可能的保留数据中的有效信号，提高噪声压制处理后的有效信号的保真度。

[0006] 本发明提供的一种地震数据中规则干扰的压制方法及装置是这样实现的：

[0007] 一种地震数据中规则干扰的压制方法，所述方法包括：

[0008] 将采集的地震数据中的地震道按照地理位置顺序关系进行排列，形成待处理地震数据炮集记录；

[0009] 对所述待处理地震数据炮集记录中判断为存在规则干扰覆盖的地震数据区域进行规则噪声频谱分析，获取所述规则干扰的估算主频；

[0010] 通过设计的滤波器从所述待处理地震数据炮集记录中分离出频率为所述估算主频的规则干扰，形成规则干扰数据记录；对所述规则干扰数据记录进行FK变换，得到FK域地震数据记录；

[0011] 在预先设定的规则干扰倾角范围内对所述FK域地震数据进行能量扫描，确定地震数据中能量最强的前K个倾角，K≥1；

[0012] 在所述规则干扰数据记录中分别对所述K个倾角进行地震数据振幅统计；若在所述倾角中连续统计到N个正的振幅值，则判断所述倾角对应的地震数据存在规则干扰，N≥2；

[0013] 在所述待处理地震数据炮集记录中，沿所述判断存在规则干扰的倾角方向进行滤波，提取所述存在规则干扰的倾角方向上的规则干扰波，形成新的干扰数据记录；

[0014] 用所述待处理地震数据炮集记录减去所述新的干扰数据记录，形成规则干扰压制后的地震数据。

[0015] 优选的实施中，所述方法还包括：利用预设频段的带通滤波器对所述新的干扰数据记录进行滤波，获取滤波后干扰数据记录；所述预设频段包括获取的用户设置的认为存在有效信号的频段范围；

[0016] 相应的，所述用所述待处理地震数据炮集记录减去所述新的干扰数据记录包括：用所述待处理地震数据炮集记录减去所述滤波后干扰数据记录。

[0017] 优选的实施中，所述N的取值范围为：5≤N≤32。

[0018] 优选的实施中，所述N的取值为：16。

[0019] 优选的实施中，述K的取值范围为：1≤K≤100。

[0020] 优选的实施中，所述K的取值为：20。

[0021] 一种地震数据中规则干扰的压制装置，所述装置包括：

[0022] 排序处理模块，用于将采集的地震数据中的地震道按照地理位置顺序关系进行排列，形成待处理地震数据炮集记录；

[0023] 干扰主频估算模块，用于对所述待处理地震数据炮集记录中判断为存在规则干扰覆盖的地震数据区域进行规则噪声频谱分析，获取所述规则干扰的估算主频；

[0024] 干扰分离及变换模块，用于通过设计的滤波器从所述待处理地震数据炮集记录中分离出频率为所述估算主频的规则干扰，形成规则干扰数据记录，以及对所述规则干扰数据记录进行FK变换，得到FK域地震数据记录；

[0025] 扫描模块，用于获取设置的K值，以及，在预先设定的规则干扰倾角范围内对所述FK域地震数据进行能量扫描，确定地震数据中能量最强的前K个倾角，K≥1；

[0026] 振幅统计模块，用于获取设置的N值，以及，在所述规则干扰数据记录中分别对所述K个倾角进行地震数据振幅统计；若在所述倾角中连续统计到N个正的振幅值，则判断所述倾角对应的地震数据存在规则干扰，N≥2；

[0027] 倾角滤波模块，用于在所述待处理地震数据炮集记录中，沿所述判断存在规则干扰的倾角方向进行滤波，提取所述存在规则干扰的倾角方向上的规则干扰波，形成新的干扰数据记录；

[0028] 干扰处理模块，用于用所述待处理地震数据炮集记录减去所述新的干扰数据记录，形成规则干扰压制后的地震数据。

[0029] 优选的实施中，所述装置还包括：规则干扰数据滤波模块，用于利用预设频段的带通滤波器对所述新的干扰数据记录进行滤波，获取滤波后干扰数据记录；所述预设频段包括获取的用户设置的认为存在有效信号的频段范围；

[0030] 相应的，所述干扰处理模块中所述用所述待处理地震数据炮集记录减去所述新的干扰数据记录包括：用所述待处理地震数据炮集记录减去所述滤波后干扰数据记录。

[0031] 优选的实施中，所述N的取值范围为：5≤N≤32。

[0032] 优选的实施中，所述N的取值为：16。

[0033] 优选的实施中，所述K的取值范围为：1≤K≤100。

[0034] 优选的实施中，所述K的取值为：20。

[0035] 本发明提供的一种地震数据中规则干扰的压制方法及装置，可以将地震数据中的地震道按照采集数据时地震道的地理位置关系重新排序，这样可以保持规则干扰的线性分布规律，增加规则干扰的相干性，使得可以有效、准确的识别出干扰信息。然后，本发明可以根据规则干扰噪声与有效信号的倾角和能量等方面的一些差异，在FK域判断出炮记录上规则干扰的倾角。进一步的可以根据FK域内判断出的噪声倾角，然后采用规则干扰同相轴在时间域内波峰连续分布的特点，在时间域重新识别和最终确定规则干扰波，并通过滤波获得相对准确可靠的规则干扰数据记录。然后可以采用从原始炮记录减去规则干扰数据记录的方法达到压制噪声的目的。本发明方法可以提高地震数据信噪比和分辨率，并且与现有技术中单纯的FK域线性噪音去噪方法相比，去噪后有效信号保真度更高，能为后期地震数据保真处理奠定坚实的基础。

附图说明

[0036] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

[0037] 图1是本发明一种地震数据中规则干扰的压制方法一种实施例的流程示意图；

[0038] 图2是本发明所述一种地震数据中规则干扰的压制方法另一种实施例的流程示意图；

[0039] 图3是本发明一种应用场景中规则干扰压制前原始单炮地震记录图；

[0040] 图4是本发明一种应用场景中利用小波分频从图1中滤出的规则干扰波炮记录；

[0041] 图5是通过本发明方法对图4所示数据判断出的滤波后干扰数据记录；

[0042] 图6是图3的原始地单炮地震记录炮记录减去图5的滤波后干扰数据记录形成规则干扰压制后的炮记录结果；

[0043] 图7是使用常规FK滤波对图1数据规则干扰压制结果；

[0044] 图8是所述常规FK滤波规则干扰压制输出的噪声记录；

[0045] 图9是图1、图6及图7数据记录的频谱对比图；

[0046] 图10是本发明所述地震数据中规则干扰的压制装置一种实施例的模块结构示意图；

[0047] 图11是本发明所述地震数据中规则干扰的压制装置另一种实施例的模块结构示意图。

具体实施方式

[0048] 为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

[0049] 下面结合附图对本发明所述的地震数据中规则干扰的压制方法进行详细的说明。图1是本发明所述一种地震数据中规则干扰的压制方法一种实施例的方法流程图。虽然本发明提供了如下述实施例或附图所示的方法操作步骤，但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中，这些步骤的执行顺序不限于本发明实施例提供的执行顺序。所述的方法的在实际中的数据迁移过程中或者装置执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。

[0050] 具体的本发明提供的一种地震数据中规则干扰的压制方法的一种实施例如图1所示，所述方法可以包括：

[0051] S1：将采集的地震数据中的地震道按照地理位置顺序关系进行排列，形成待处理地震数据炮集记录。

[0052] 野外采集的地震数据通常可以包括按照道序或者时序等一定方式排序的数据。本发明中获取采集地震数据后，可以按照地震数据中地震道的地理位置顺序关系进行重新排列，形成待处理地震数据炮集记录。具体的，例如可以在测线中按照地震道炮检距从小到大的顺序进行排列等。当然，本领域人员也可以采取其他的方式按照地震道地理位置顺序进行排列，以实现增强规则干扰噪声的相干性，保持规则干扰的线性分布规则，便于规则干扰噪声的识别和预测。

[0053] S2：对所述待处理地震数据炮集记录中判断为存在规则干扰覆盖的地震数据区域进行规则噪声频谱分析，获取所述规则干扰的估算主频。

[0054] 本发明中可以在所述待处理地震数据中判断查找到存在规则干扰覆盖的地震数据区域。所述的判断方法可以包括采用的预定的算法或数据处理方法对所述待处理地震数据运算处理后确认得出，例如基于对地震数据估算的视速度或地震数据的振幅或相干性计算等。本发明的一种实施方式中，也可以通过作业人员根据经验在所述待处理地震数据炮集记录上识别出的规则干扰和分布进而确定出的存在规则干扰。

[0055] 然后，可以在所述规则干扰覆盖的地震数据区域进行规则噪声频谱分析，获取所述规则干扰的估算主频。具体的例如可以将估算出的规则干扰的主频作为所述规则干扰的估算主频。

[0056] S3：通过设计的滤波器从所述待处理地震数据炮集记录中分离出频率为所述估算主频的规则干扰，形成规则干扰数据记录，对所述规则干扰数据记录进行FK变换，得到FK域地震数据记录。

[0057] 本实施例中可以根据估算出的噪声的估算主频设计一个小波变换滤波器，滤波的主频可以与所述估算主频相同，或者在以所述估算主频为中心的预定阈值范围内。本实施例可以通过小波分频对地震数据进行分频滤波处理，通过分频把规则干扰从炮记录中分离出来，形成一个以规则干扰为主的规则干扰数据记录。然后可以对这个规则干扰数据记录做FK变换，形成FK域地震数据记录。

[0058] S4：在预先设定的规则干扰倾角范围内对所述FK域地震数据进行能量扫描，确定地震数据中能量最强的前K个倾角，K≥1。

[0059] 本发明中可以设定一个规则干扰倾角范围。本申请中所述的规则干扰倾角一半的可以包括确定的规则干扰传播的视速度范围。例如本发明的一种实施例中可以将所述视速度设置为选择500m/s-2000m/s，在500m/s-2000m/s范围内对形成的FK域地震数据进行倾角能量扫描，确定地震数据中能量最强的前K个倾角。本发明的一种实施例中，所述的倾角个数K的取值范围可以包括1-00。具体的可以根据初级处理需求进行设置。本发明的一种实施例中，所述K的取值可以为20，即可以扫描出能量最强的20个倾角。然后可以把这些倾角的大小记录下来。

[0060] S5：在所述规则干扰数据记录中分别对所述K个倾角进行地震数据振幅统计；若在所述倾角中连续统计到N个正的振幅值，则判断所述倾角对应的地震数据存在规则干扰，N≥2。

[0061] 本发明的一种实施例，可以在小波分频滤除的规则干扰数据记录中，沿能量最强的20个倾角进行地震数据振幅统计，对连续统计到K个正的振幅值的倾角进行噪声压制，相反就不做任何处理。

[0062] 本发明中所述的N的取值可以根据规则干扰的具体分布情况进行设置、确定。一般的，本发明中所述的N的取值范围可以包括：5≤N≤32。本实施例中所述的N可以取值为16，即对连续统计到16个正的振幅值的倾角进行噪声压制。

[0063] S6：在所述待处理地震数据炮集记录中，沿所述判断存在规则干扰的倾角方向进行滤波，提取所述存在规则干扰的倾角方向上的规则干扰波，形成新的干扰数据记录。

[0064] 可以对判断存在规则干扰的倾角方向上的数据进行噪声压制。本发明一种实施例中，具体的压制方法可以包括：根据判断出的规则干扰倾角以及估算出的噪声主频，设计一个滤波器，沿在噪声炮记录中判断出的规则干扰倾角方向进行滤波，提取所有存在规则干扰倾角上的规则干扰波，形成一个新的干扰数据记录。

[0065] S7：用所述待处理地震数据炮集记录减去所述新的干扰数据记录，形成规则干扰压制后的地震数据。

[0066] 利用上述所述的处理方式获取所述新的干扰数据记录后，可以用去噪前的所述待处理地震数据炮记录减去本发明实施例获取的新的干扰数据记录，可以完成数据记录的噪声压制。

[0067] 本发明的一种优选实施例中，在获取所述新的干扰数据记录后，还可以进一步的在所述新的干扰数据记录中进行筛选，以过滤可能存在的有效信号，进而减少对有效信号的损害，提高数据处理后有效信号的保真度。图2是本发明所述一种地震数据中规则干扰的压制方法另一种实施例的流程示意图，如图2所示，本发明的另一种实施例中，所述方法还可以包括：

[0068] S701：利用预设频段的带通滤波器对所述新的干扰数据记录进行滤波，获取滤波后干扰数据记录；所述预设频段包括获取的用户设置的认为存在有效信号的频段范围；

[0069] 相应的，所述用所述待处理地震数据炮集记录减去所述新的干扰数据记录包括：用所述待处理地震数据炮集记录减去所述滤波后干扰数据记录。

[0070] 具体的实现方式中，可以对上述中获得的新的规则干扰数据记录，再设计一个带通滤波器，用户可以参考所述规则干扰的估算主频确定所述带通滤波的范围。这个滤波器主要是滤除用户认为可能在所述新的干扰数据记录中还存在的有效信号，例如可以为1Hz，5Hz，20Hz或者30Hz。

[0071] 本发明实施例提供的地震数据中规则干扰的压制方法，可以将地震数据中的地震道按照采集数据时地震道的地理位置关系重新排序，这样可以保持规则干扰的线性分布规律，增加规则干扰的相干性，使得可以有效、准确的识别出干扰信息。本发明方法不仅可以提高地震数据信噪比和分辨率，并且与现有技术中单纯的FK域线性噪音去噪方法相比，去噪后有效信号保真度更高，能为后期地震数据保真处理奠定坚实的基础。

[0072] 下面是本发明所述方法运用于实际地震数据单炮记录进行线性噪音衰减的一个应用场景的实施例。

[0073] 图3是应用场景中规则干扰压制前原始单炮地震记录图。图3中所示的原始单炮地震记录已经按地理位置进行了排序，图3横坐标方向显示的是地震道道号，纵坐标方向显示的是每个地震道的样点时间。估算炮记录上的规则干扰噪声和分布区域，然后在规则干扰噪声覆盖的地震数据区域进行规则噪声频谱分析，图3中大致估算噪声的主频是10Hz。设计一个小波变换滤波器，滤波的主频为10Hz，通过小波分频把规则干扰从炮记录中分离出来，形成一个以规则干扰为主的新的噪声炮记录，如图4所示的利用小波分频滤出的规则干扰波炮记录。从图4中可以看出小波分频形成的规则干扰波炮记录中还存在有效信号。因此，本应用场景中可以对这个噪声炮记录数据做FK变换，形成FK域地震数据。给定一个规则干扰倾角范围，一般也是规则干扰传播的视速度范围，本实施例可以选择500m/s-2000m/s，在500m/s-2000m/s范围内对形成的FK数据进行倾角能量扫描。本实施例扫描出了能量最强的
20个倾角，把这些倾角的大小记录下来，在小波分频滤除的规则干扰炮记录中，沿能量最强的20个倾角进行地震数据振幅统计。对连续统计到16个正的振幅值的倾角进行噪声压制，相反就不做任何处理。具体的噪声压制方法可以包括：根据判断出的规则干扰倾角以及估算出的噪声主频，设计一个滤波器，沿在噪声炮记录中判断出的规则干扰倾角方向进行滤波，提取所有存在规则干扰倾角上的规则干扰波，形成一个新的规则干扰波炮记录。对获得的新的规则干扰波炮记录，再设计一个通带范围为(1Hz，5Hz，20Hz，30Hz)的带通滤波器，用这个带通滤波器对新的规则干扰波炮记录进行滤波处理，得到滤波后干扰数据记录。图5是通过本发明方法对图4所示数据判断出的滤波后干扰数据记录。最后可以用去噪前的原始的地震数据炮记录减去发明得到的比较准确的滤波后干扰数据记录，完成单炮地震数据记录的噪声压制。图6是图3的原始地单炮地震记录炮记录减去图5的滤波后干扰数据记录形成规则干扰压制后的炮集记录结果。

[0074] 本发明提供的地震数据中规则干扰的压制方法，可以有效减小地震数据规则干扰去噪时的有效信号损伤，在压制规则干扰的同时尽可能的保留数据中的有效信号，提高噪声压制处理后的有效信号的保真度。图7是使用常规FK滤波对图1数据规则干扰压制结果，图8是所述常规FK滤波规则干扰压制输出的噪声记录，对比图8和图5可以明显看出图8中还存在部分有效信号能量，而图5中几乎看不到有效信号。同时，图9是图1、图6及图7数据记录的频谱对比图，从图9中(line_1为去噪前图1炮记录频谱，line_2为利用发明规则干扰压制后图6炮记录频谱，line_3为常规FK域规则干扰压制后图7炮记录频谱)的对比可以发现，在10Hz左右，本发明方法比常规FK去噪更能保护有效信号能量，本发明规则干扰压制后在高频部分与原始地震数据的频谱也几乎完全一致，而常规FK滤波却存在一点的差异。这可以说明本发明提供的地震数据中规则干扰的压制方法可以提高规则干扰预测精度，在很大程度上提高了地震数据的信噪比和分辨率。同时由于采用在时间域进行噪声衰减，与常规FK域规则噪声压制相比，更能够保持规则干扰压制后地震数据的波形，避免去噪后产生假频。

[0075] 基于本发明所述的地震数据中规则干扰的压制方法，本发明提供一种地震数据中规则干扰的压制装置。图10是本发明所述地震数据中规则干扰的压制装置一种实施例的模块结构示意图，如图10所示，所述装置可以包括：

[0076] 排序处理模块101，可以用于将采集的地震数据中的地震道按照地理位置顺序关系进行排列，形成待处理地震数据炮集记录；

[0077] 干扰主频估算模块102，可以用于对所述待处理地震数据炮集记录中判断为存在规则干扰覆盖的地震数据区域进行规则噪声频谱分析，获取所述规则干扰的估算主频；

[0078] 干扰分离及变换模块103，可以用于通过设计的滤波器从所述待处理地震数据炮集记录中分离出频率为所述估算主频的规则干扰，形成规则干扰数据记录，以及对所述规则干扰数据记录进行FK变换，得到FK域地震数据记录；

[0079] 扫描模块104，可以用于获取设置的K值，以及，在预先设定的规则干扰倾角范围内对所述FK域地震数据进行能量扫描，确定地震数据中能量最强的前K个倾角，K≥1；

[0080] 振幅统计模块105，可以用于获取设置的N值，以及，在所述规则干扰数据记录中分别对所述K个倾角进行地震数据振幅统计；若在所述倾角中连续统计到N个正的振幅值，则判断所述倾角对应的地震数据存在规则干扰，N≥2；

[0081] 倾角滤波模块106，可以用于在所述待处理地震数据炮集记录中，沿所述判断存在规则干扰的倾角方向进行滤波，提取所述存在规则干扰的倾角方向上的规则干扰波，形成新的干扰数据记录；

[0082] 干扰处理模块107，可以用于用所述待处理地震数据炮集记录减去所述新的干扰数据记录，形成规则干扰压制后的地震数据。

[0083] 图11是本发明所述地震数据中规则干扰的压制装置另一种实施例的模块结构示意图，如图11所示，所述装置还可以包括：

[0084] 规则干扰数据滤波模块108，可以用于利用预设频段的带通滤波器对所述新的干扰数据记录进行滤波，获取滤波后干扰数据记录；所述预设频段包括获取的用户设置的认为存在有效信号的频段范围；

[0085] 相应的，所述干扰处理模块107中所述用所述待处理地震数据炮集记录减去所述新的干扰数据记录包括：用所述待处理地震数据炮集记录减去所述滤波后干扰数据记录。

[0086] 本发明所述装置中各个模块涉及的具体描述可以参照上述所述方法的描述，例如，所述N的取值范围可以包括：5≤N≤32，优选的实施例中所述N的取值可以为：16。或者，其他的实施例中，所述K的取值范围可以包括：1≤K≤100，优选的实施例中，所述K的取值可以为：20等，各个实施例之间可以互相参照即可，在此不做赘述。所述的装置可以理解为实现本发明所述方法的一种产品的实现方式。

[0087] 本发明提供一种地震数据中规则干扰的压制方法及装置，可以通过小波分频首先分离噪声和有效信号，并且在一定程度上保持了规则干扰噪声的波形特征，是通过FK变换，把相同倾角的规则噪声能量在FK域内叠加起来，增强了规则干扰的能量，便于在FK域内对规则干扰倾角进行判断和识别，然后在时间域利用这些倾角进行噪声再次判断和识别，通过对噪声进行带通滤波，更好保护了有效信号的能量。与常规FK规则干扰滤波相比，能更好地保持去噪后地震数据的波形特征。

[0088] 尽管本发明内容中提到数据滤波、FK变换、能量扫描、振幅统计等的数据变换、处理的描述，但是，本发明并不局限于必须是完全标准或者所提及的方式的数据变换、处理的情况。本发明中各个实施例所涉及的上述描述仅是本发明中的一些实施例中的应用，在某些标准、方法的基础上略加修改后的处理方法也可以实行上述本发明各实施例的方案。当然，在符合本发明上述各实施例的中所述的处理方法步骤的其他无创造性的变形，仍然可以实现相同的申请，在此不再赘述。

[0089] 虽然本发明提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。

[0090] 上述实施例阐明的装置或模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本发明时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然，也可以将实现某功能的模块由多个子模块或子单元组合实现。

[0091] 本发明中所述的方法、装置或模块可以以计算机可读程序代码的控制器按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

[0092] 本发明所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本发明，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

[0093] 通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本发明可借助软件加必需的硬件的方式来实现。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。

[0094] 本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本发明的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

[0095] 虽然通过实施例描绘了本发明，本领域普通技术人员知道，本发明有许多变形和变化而不脱离本发明的精神，希望所附的权利要求包括这些变形和变化而不脱离本发明的精神。