一种地震采集脚印噪音压制方法及系统

摘要

本申请实施例公开了一种地震采集脚印噪音压制方法及系统，所述方法包括：获取原始叠后地震数据；调整原始叠后地震数据的排列方式，得到标准原始叠后地震数据；获取标准原始叠后地震数据中第一时间切片的地震数据，对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据；将拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子；将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱；利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。可以提高自动识别并压制采集脚印噪音以及保持有效地质信号的能力。

权利要求

1.一种地震采集脚印噪音压制方法，其特征在于，包括：

获取原始叠后地震数据；

调整所述原始叠后地震数据的排列方式，得到标准原始叠后地震数据；

获取所述标准原始叠后地震数据中第一时间切片的地震数据，对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据；

将所述拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱，根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子；

将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱；

利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。

2.根据权利要求1所述一种地震采集脚印噪音压制方法，其特征在于，在利用所述第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制之后，所述方法还包括：将采集脚印噪音压制后的第一时间切片的地震数据波数谱从波数域转换到空间域，得到中间叠后地震数据；

调整所述中间叠后地震数据的排列方式，得到目标叠后地震数据。

3.根据权利要求2所述一种地震采集脚印噪音压制方法，其特征在于，通过矩阵转置调整所述中间叠后地震数据的排列方式。

4.根据权利要求1所述一种地震采集脚印噪音压制方法，其特征在于，所述根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子，包括：根据所述拉普拉斯变换地震数据波数谱，获取所述拉普拉斯变换地震数据波数谱中的最大波数值和最小波数值；

根据预设规则以及所述最大波数值和所述最小波数值，调整所述拉普拉斯变换地震数据波数谱中的波数值；

将所述调整后的拉普拉斯变换地震数据波数谱作为第一时间切片的滤波因子。

5.根据权利要求4所述一种地震采集脚印噪音压制方法，其特征在于，所述根据预设规则以及所述最大波数值和所述最小波数值，调整所述拉普拉斯变换地震数据波数谱中的波数值，包括：采用下述公式实现：

公式中，Ai表示所述调整后的拉普拉斯变换地震数据波数谱中第i个波数值；A0max表示所述调整前的拉普拉斯变换地震数据波数谱中的最大波数值；A0min表示所述调整前的拉普拉斯变换地震数据波数谱中的最小波数值；A0i表示所述调整前的拉普拉斯变换地震数据波数谱中第i个波数值。

6.根据权利要求1所述一种地震采集脚印噪音压制方法，其特征在于，所述利用所述第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制，包括：将第一时间切片的滤波因子与第一时间切片的地震数据波数谱在波数域相乘。

7.根据权利要求1所述一种地震采集脚印噪音压制方法，其特征在于，所述原始叠后地震数据包括：时间、振幅值、主测线号和联络测线号。

8.根据权利要求7所述一种地震采集脚印噪音压制方法，其特征在于，所述调整原始叠后地震数据的排列方式，包括：调整所述原始叠后地震数据的至少两个维度的排列方式；

具体地，将所述原始叠后地震数据的至少两个维度中的时间维度设置为最慢维度；所述原始叠后地震数据的至少两个维度包括：所述时间对应的时间维度、所述主测线号对应的主测线维度和所述联络测线号对应的联络测线维度。

9.根据权利要求1所述一种地震采集脚印噪音压制方法，其特征在于，通过矩阵转置调整所述原始叠后地震数据的排列方式。

10.一种地震采集脚印噪音压制系统，其特征在于，包括：原始叠后地震数据获取单元、排列方式调整单元、拉普拉斯变换单元、滤波因子获取单元、傅里叶变换单元和噪音压制单元；其中，所述原始叠后地震数据获取单元，用于获取原始叠后地震数据；

所述排列方式调整单元，用于调整所述原始叠后地震数据的排列方式，得到标准原始叠后地震数据；

所述拉普拉斯变换单元，用于获取所述标准原始叠后地震数据中第一时间切片的地震数据，对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据；

所述滤波因子获取单元，用于将所述拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱，根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子；

所述傅里叶变换单元，用于将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱；

所述噪音压制单元，用于利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。

11.根据权利要求10所述一种地震采集脚印噪音压制系统，其特征在于，所述系统还包括：中间叠后地震数据获取单元和目标叠后地震数据获取单元；其中，所述中间叠后地震数据获取单元，用于将采集脚印噪音压制后的第一时间切片的地震数据波数谱从波数域转换到空间域，得到中间叠后地震数据；

所述目标叠后地震数据获取单元，用于调整所述中间叠后地震数据的排列方式，得到目标叠后地震数据。

说明书

一种地震采集脚印噪音压制方法及系统

技术领域

[0001] 本申请涉及石油地震勘探中的地震数据处理技术领域，特别涉及一种地震采集脚印噪音压制方法及系统。

背景技术

[0002] 在地震数据的采集过程中所使用的观测系统会对采集到的地震数据产生干扰，可以形象地将这种干扰被喻为在地震数据中留下了脚印，即采集脚印。这种脚印表现为在叠后地震数据的时间切片上出现有规律的振幅变化，是调制在地层反射信号上的一种系统性噪声，它严重地影响了地震属性反演的精确度。

[0003] 通常有三种途径衰减采集脚印噪音，第一种途径是通过采集系统参数设置使不同偏移距的道数变化最小；第二种途径是通过叠前处理使要叠加的道集间差异最小；第三种途径是通过叠后处理压制采集脚印。地震数据处理实例表明，通过优化采集参数以及叠前处理措施后，采集脚印依然发育。因此，需要对叠后地震数据进行噪音压制，以提高地震数据信噪比和突出有效反射信号。

[0004] 目前基于叠后地震勘探数据压制采集脚印噪音的方法主要是频率波数(F-Kx-Ky)滤波方法，这种方法是通过二维傅里叶变换将原始地震数据由空间域变换至波数域，得到原始地震数据波数谱，根据采集脚印周期性特点设计滤波因子，利用所述滤波因子对原始地震数据波谱数进行采集脚印噪音压制，得到目标地震数据波谱数，通过二维傅里叶反变换将所述目标地震数据波谱数由波数域变换至空间域，完成采集脚印噪音压制。

[0005] 发明人发现现有技术中至少存在如下问题：

[0006] 部分采集脚印，例如海上采集脚印噪音，具有随时间和空间变化的特征，在整个地震数据中没有明显的周期性分布特点。由于现有技术中频率波数滤波方法是根据采集脚印周期性特点设计滤波因子来完成采集脚印噪音压制的，所以，该方法很难有效识别和压制所述随时间和空间变化的采集脚印，对地质信号保持能力较差。

发明内容

[0007] 本申请实施例的目的是提供一种地震采集脚印噪音压制方法及系统，以提高自动识别并压制采集脚印噪音以及保持有效地质信号的能力。

[0008] 为解决上述技术问题，本申请实施例提供一种地震采集脚印噪音压制方法及系统是这样实现的：

[0009] 一种地震采集脚印噪音压制方法，包括：

[0010] 获取原始叠后地震数据；

[0011] 调整所述原始叠后地震数据的排列方式，得到标准原始叠后地震数据；

[0012] 获取所述标准原始叠后地震数据中第一时间切片的地震数据，对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据；

[0013] 将所述拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱，根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子；

[0014] 将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱；

[0015] 利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。

[0016] 优选方案中，在利用所述第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制之后，所述方法还包括：

[0017] 将采集脚印噪音压制后的第一时间切片的地震数据波数谱从波数域转换到空间域，得到中间叠后地震数据；

[0018] 调整所述中间叠后地震数据的排列方式，得到目标叠后地震数据。

[0019] 优选方案中，通过矩阵转置调整所述中间叠后地震数据的排列方式。

[0020] 优选方案中，所述根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子，包括：

[0021] 根据所述拉普拉斯变换地震数据波数谱，获取所述拉普拉斯变换地震数据波数谱中的最大波数值和最小波数值；

[0022] 根据预设规则以及所述最大波数值和所述最小波数值，调整所述拉普拉斯变换地震数据波数谱中的波数值；

[0023] 将所述调整后的拉普拉斯变换地震数据波数谱作为第一时间切片的滤波因子。

[0024] 优选方案中，所述根据预设规则以及所述最大波数值和所述最小波数值，调整所述拉普拉斯变换地震数据波数谱中的波数值，包括：

[0025] 采用下述公式实现：

[0026]

[0027] 公式中，Ai表示所述调整后的拉普拉斯变换地震数据波数谱中第i个波数值；A0max表示所述调整前的拉普拉斯变换地震数据波数谱中的最大波数值；A0min表示所述调整前的拉普拉斯变换地震数据波数谱中的最小波数值；A0i表示所述调整前的拉普拉斯变换地震数据波数谱中第i个波数值。

[0028] 优选方案中，所述利用所述第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制，包括：将第一时间切片的滤波因子与第一时间切片的地震数据波数谱在波数域相乘。

[0029] 优选方案中，所述地震数据包括：时间、振幅值、主测线号和联络测线号。

[0030] 优选方案中，所述调整原始叠后地震数据的排列方式，包括：调整所述原始叠后地震数据的至少两个维度的排列方式；

[0031] 具体地，将所述原始叠后地震数据的至少两个维度中的时间维度设置为最慢维度；所述原始叠后地震数据的至少两个维度包括：所述时间对应的时间维度、所述主测线号对应的主测线维度和所述联络测线号对应的联络测线维度。

[0032] 优选方案中，通过矩阵转置调整所述原始叠后地震数据的排列方式。

[0033] 一种地震采集脚印噪音压制系统，包括：原始叠后地震数据获取单元、排列方式调整单元、拉普拉斯变换单元、滤波因子获取单元、傅里叶变换单元和噪音压制单元；其中，[0034] 所述原始叠后地震数据获取单元，用于获取原始叠后地震数据；

[0035] 所述排列方式调整单元，用于调整所述原始叠后地震数据的排列方式，得到标准原始叠后地震数据；

[0036] 所述拉普拉斯变换单元，用于获取所述标准原始叠后地震数据中第一时间切片的地震数据，对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据；

[0037] 所述滤波因子获取单元，用于将所述拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱，根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子；

[0038] 所述傅里叶变换单元，用于将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱；

[0039] 所述噪音压制单元，用于利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。

[0040] 优选方案中，所述系统还包括：中间叠后地震数据获取单元和目标叠后地震数据获取单元；其中，

[0041] 所述中间叠后地震数据获取单元，用于将采集脚印噪音压制后的第一时间切片的地震数据波数谱从波数域转换到空间域，得到中间叠后地震数据；

[0042] 所述目标叠后地震数据获取单元，用于调整所述中间叠后地震数据的排列方式，得到目标叠后地震数据。

[0043] 本申请提供了一种地震采集脚印噪音压制方法及系统，基于时间切片地震数据的拉普拉斯变换和波数滤波结合的方式，所述方法对叠后地震数据进行采集脚印压制，能够根据海上拖缆地震数据采集脚印随空间和时间变化的特点，自动识别并压制采集脚印噪音，并且能够有效保持地质信号。

附图说明

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

[0045] 图1是本申请一种地震采集脚印噪音压制方法实施例的流程图；

[0046] 图2是本申请第1000毫秒处的时间切片的地震数据及其波数谱的示意图；

[0047] 图3是本申请第1000毫秒处的时间切片的地震数据拉普拉斯变换后的拉普拉斯变换地震数据和拉普拉斯变换地震数据波数谱的示意图；

[0048] 图4是本申请采集脚印噪音压制后的第1000毫秒处的时间切片的地震数据及其波数谱的示意图；

[0049] 图5是本申请采集脚印噪音压制前后的第1000毫秒处的时间切片的地震数据差及其波数谱的示意图；

[0050] 图6是本申请一种地震采集脚印噪音压制系统实施例的组成结构图。

具体实施方式

[0051] 本申请实施例提供一种地震采集脚印噪音压制方法及系统。

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

[0053] 图1是本申请一种地震采集脚印噪音压制方法实施例的流程图。如图1所示，所述地震采集脚印噪音压制方法，包括以下步骤。

[0054] 步骤S101：获取原始叠后地震数据。

[0055] 具体地，可以将原始叠后地震数据读入计算机内存。其中，将地表接受到的来自地下同一反射点的地震数据经过动校正后进行叠加，可以有效提高地震数据的信噪比，叠加后的地震数据，称为叠后地震数据。所述原始叠后地震数据可以包括：时间、振幅值、主测线(Inline)号和联络测线(Xline)号。所述原始叠后地震数据为三维地震数据，可以包括：时间方向、Inline方向和Xline方向。其中，所述原始叠后地震数据中Inline方向和Xline方向相互垂直，时间方向与Inline方向和Xline方向均垂直。

[0056] 例如，所述原始叠后地震数据中的振幅值可以为所述原始叠后地震数据中Inline方向上的Inline号5574和Xline方向上的Xline号3221在时间方向上的时间1000毫秒处对应的时间样点的振幅值。

[0057] 步骤S102：调整所述原始叠后地震数据的排列方式，得到标准原始叠后地震数据。

[0058] 具体地，所述原始叠后地震数据的维度可以包括：所述时间对应的时间维度、所述主测线号对应的主测线维度和所述联络测线号对应的联络测线维度。所述原始叠后地震数据是按照Inline维度-Xline维度-时间维度的序列进行排列的。由于采集脚印沿时间切片具有条带状等规律性特征，可以通过矩阵转置将所述原始叠后地震数据的排列方式调整为时间维度-Inline维度-Xline维度的序列。即可以将所述时间维度设置为最慢维度。可以将所述排列方式调整后的原始叠后地震数据作为标准原始叠后地震数据。

[0059] 例如，所述矩阵转置的过程可以采用下述方式实现：

[0060] 定义矩阵A转置为n×m阶矩阵B，满足B＝a(j，i)，即b(i，j)＝a(j，i)(B的第i行第j列元素是A的第j行第i列元素)，记AT＝B。二维矩阵 转置 同理三维数组的转置，保持一维不变，另外两维数组按照二维矩阵转置方式实现转置。

[0061] 步骤S103：获取所述标准原始叠后地震数据中第一时间切片的地震数据，对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据。

[0062] 具体地，可以将所述标准原始叠后地震数据由时间域转换到时间切片域，可以得到第一时间处对应的第一时间切片的地震数据。可以对获得的第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据，可以使所述标准原始叠后地震数据中的采集脚印噪音在拉普拉斯变换地震数据中突显出来。其中，第一时间切片的地震数据可以是所述标准原始叠后地震数据中任意时间处对应的任意时间切片的地震数据。

[0063] 例如，如图2所示，图2中左图为第1000毫秒处的时间切片的地震数据。如图3所示，图3中左图为图2中第1000毫秒处的时间切片的地震数据拉普拉斯变换后的拉普拉斯变换地震数据。从图3中左图中可以看出，第1000毫秒处的时间切片的地震数据经过拉普拉斯变换以后，条带状的采集脚印噪音在所述拉普拉斯变化地震数据中突显出来，从而使采集脚印噪音得到有效表征。

[0064] 步骤S104：将所述拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱，根据所述拉普拉斯变换地震数据波数谱，获取第一时间切片的滤波因子。

[0065] 具体地，将图3中左图中拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波谱数。根据图3中右图中所述拉普拉斯变换地震数据波数谱，可以获取所述拉普拉斯变换地震数据波数谱中的最大波数值和最小波数值。根据预设规则以及最大波数值和最小波数值，可以调整所述拉普拉斯变换地震数据波数谱中的波数值。所述预设规则可以采用下述公式实现：

[0066]

[0067] 公式中，Ai表示所述调整后的拉普拉斯变换地震数据波数谱中第i个波数值；A0max表示所述调整前的拉普拉斯变换地震数据波数谱中的最大波数值；A0min表示所述调整前的拉普拉斯变换地震数据波数谱中的最小波数值；A0i表示所述调整前的拉普拉斯变换地震数据波数谱中第i个波数值。按照上述公式对所述拉普拉斯变化地震数据波数谱的波数值进行调整，将调整后的拉普拉斯变换地震数据波数谱作为第一时间切片的滤波因子。所述得到的第一时间切片的滤波因子是基于第一时间切片的地震数据得到的。因此，每一个时间切片的滤波因子可以随着拉普拉斯变换后识别的每一个时间切片的地震数据中的采集脚印噪音的变化而变化，可以对每一个时间切片的地震数据中的采集脚印噪音进行有效压制。

[0068] 例如：如图3所示，图3中右图为图3中左图中所述拉普拉斯变换地震数据经过二维傅里叶变换后，得到的拉普拉斯变换地震数据波谱数。按照所述公式对图3中右图中的拉普拉斯变换地震数据波谱数进行波数值调整，将调整后的拉普拉斯变换地震数据波数谱作为第1000毫秒处的时间切片的滤波因子。

[0069] 步骤S105：将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱。

[0070] 具体地，将所述标准原始叠后地震数据中第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱，以便利用上述步骤中的滤波因子在波数域对第一时间切片的地震数据进行采集脚印压制。

[0071] 例如，将图2中左图中第1000毫秒处的时间切片的地震数据进行二维傅里叶变换，得到第1000毫秒处的时间切片的地震数据波数谱。图2中右图为经过二维傅里叶变换后的第1000毫秒处的时间切片的地震数据的波谱数。

[0072] 需要说明的是，步骤S105可以在步骤S103和/或步骤S104之前或之后，本申请对此并不作出限定。

[0073] 步骤S106：利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。

[0074] 具体地，可以将第一时间切片的滤波因子与第一时间切片的地震数据波数谱在波数域相乘，对第一时间切片的地震数据波数谱进行采集脚印压制，得到所述采集脚印压制后的第一时间切片的地震数据波数谱。

[0075] 例如，图4是本申请采集脚印噪音压制后的第1000毫秒处的时间切片的地震数据及其波数谱的示意图。如图4所示，图4中左图为采集脚印噪音压制后的第1000毫秒处的时间切片的地震数据，图4中右图为经过二维傅里叶变换后的采集脚印噪音压制后的第1000毫秒处的时间切片的地震数据波数谱。对比图2中左图和图4中左图，可以看出，条带状的采集脚印噪音可以得到有效压制。从图4中右图中可以看出，低、高波数类采集脚印噪音均可以得到压制。

[0076] 图5是本申请采集脚印噪音压制前后的第1000毫秒处的时间切片的地震数据差及其波数谱的示意图。如图5所示，图5中左图为采集脚印噪音压制前后的第1000毫秒处的时间切片的地震数据差，图5中右图为经过二维傅里叶变换后的采集脚印噪音压制前后的第1000毫秒处的时间切片的地震数据差波数谱。从图5可以看出，采集脚印压制前后地质信息得到有效保持。

[0077] 在另一个实施方式中，还可以将采集脚印噪音压制后的第一时间切片的地震数据波数谱从波数域转换到空间域，得到中间叠后地震数据。可以调整所述中间叠后地震数据的排列方式，得到目标叠后地震数据。

[0078] 具体地，可以采用上述实施方式对所述标准叠后地震数据中除第1000毫秒处的时间切片以外其他时间处对应的时间切片的地震数据进行采集脚印压制，得到采集脚印噪音压制后的所述标准叠后地震数据中所有时间切片的地震数据波数谱。采用二维傅里叶反变换，可以将所述采集脚印噪音压制后的所述标准叠后地震数据中所有时间切片的地震数据波数谱从波数域转换到时间切片域，然后将其再由时间切片域转换到时间域，得到中间叠后地震数据。可以通过矩阵转置将所述中间叠后地震数据的排列方式调整为Inline维度-Xline维度-时间维度的序列，得到目标叠后地震数据。从而实现了整个原始叠后地震数据的采集脚印噪声压制。其中，这一个步骤中矩阵转置的方法与步骤S102中矩阵转置的方法相同。

[0079] 所述地震采集脚印噪音压制方法实施例，基于时间切片地震数据的拉普拉斯变换和波数滤波结合的方式，对叠后地震数据进行采集脚印压制，能够根据海上拖缆地震数据采集脚印随空间和时间变化的特点，自动识别并压制采集脚印噪音，并且能够有效保持地质信号。

[0080] 图6是本申请一种地震采集脚印噪音压制系统实施例的组成结构图。如图6所示，所述地震采集脚印噪音压制系统可以包括：原始叠后地震数据获取单元100、排列方式调整单元200、拉普拉斯变换单元300、滤波因子获取单元400、傅里叶变换单元500和噪音压制单元600。

[0081] 所述原始叠后地震数据获取单元100，可以用于获取原始叠后地震数据。

[0082] 所述排列方式调整单元200，可以用于调整所述原始叠后地震数据的排列方式，得到标准原始叠后地震数据。

[0083] 所述拉普拉斯变换单元300，可以用于获取所述标准原始叠后地震数据中第一时间切片的地震数据。可以对第一时间切片的地震数据进行拉普拉斯变换，得到拉普拉斯变换地震数据。

[0084] 所述滤波因子获取单元400，可以用于将所述拉普拉斯变换地震数据进行二维傅里叶变换，得到拉普拉斯变换地震数据波数谱。根据所述拉普拉斯变换地震数据波数谱，可以获取第一时间切片的滤波因子。

[0085] 所述第二傅里叶变换单元500，可以用于将第一时间切片的地震数据进行二维傅里叶变换，得到第一时间切片的地震数据波数谱。

[0086] 所述噪音压制单元600，可以用于利用第一时间切片的滤波因子对第一时间切片的地震数据波数谱进行采集脚印噪音压制。

[0087] 在另一个实施方式中，所述地震采集脚印噪音压制系统还可以包括：中间叠后地震数据获取单元700和目标叠后地震数据获取单元800。

[0088] 所述中间叠后地震数据获取单元700，可以用于将采集脚印噪音压制后的第一时间切片的地震数据波数谱从波数域转换到空间域，得到中间叠后地震数据。

[0089] 所述目标叠后地震数据获取单元800，可以用于调整所述中间叠后地震数据的排列方式，得到目标叠后地震数据。

[0090] 所述地震采集脚印噪音压制系统实施例与所述地震采集脚印噪音压制方法实施例相对应，可以实现对叠后地震数据进行采集脚印压制，能够根据海上拖缆地震数据采集脚印随空间和时间变化的特点，自动识别并压制采集脚印噪音，并且能够有效保持地质信号。

[0091] 在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable Gate Array，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片2。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等，目前最普遍使用的是VHDL(Very-High-Speed Integrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

[0092] 控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：ARC 625D、Atmel AT91SAM、Microchip PIC18F26K20以及Silicone Labs C8051F320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。

[0093] 本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

[0094] 上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

[0095] 为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

[0096] 通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。

[0097] 本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

[0098] 本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

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

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