一种地震数据的采集脚印的处理方法及装置

技术领域

本发明涉及地震勘探技术领域，具体涉及一种地震数据的采集脚印的处理方法及装置。

背景技术

当前，油气勘探正向广度和深度两个方面发展。前者的目标是要发现新的远景勘探领域油气田；后者要求寻找那些埋藏深度大、复杂程度高的油气藏，同时要求解决精细构造和储层参数精细描述等问题。随着岩性油气藏和隐蔽性油气藏勘探的不断深入，地质目标也从原来的构造找油气转向了岩性找油气，要很好解决油气藏以及储层参数精细描述问题，就要涉及到成像数据信噪比问题。

采集脚印噪声又称采集痕迹噪声，是一种人为因素产生的地震噪声，是地震资料采集与处理过程中留下的人为痕迹，表现方式为地震剖面深度或者时间切片上出现规则性的振幅变化假象，部分采集脚印噪声能量很强，经常会严重影响油气藏的预测精度和效率。对于这种采集脚印的压制，通常叠前一般可以通过组合和插值的办法对采集脚印进行部分衰减；叠后的采集脚印压制方法比较多，包括倾角滤波、F-Kx-Ky滤波方法、小波变换以及自适应类方法等。DMO或叠前偏移过程中采用加权函数也能够部分衰减采集脚印。但是不同采集脚印产生机理不同，不同采集脚印压制方法都具备自身适用性以及方法本身的局限性，因此开展有针对性的采集脚印压制方法研究具有重要的意义。

发明内容

本发明实施例提供一种地震数据的采集脚印的处理方法及装置，针对随机分布或者局部相关分布的采集脚印，解决如何提高地震数据信噪比和分辨率的问题。

一方面，本发明提出一种地震数据的采集脚印的处理方法，包括：

对地震数据进行叠前预处理，获得叠后地震数据；

对所述叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；

对所述补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；

对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；

对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；

对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据。

另一方面，本发明提供一种地震数据的采集脚印的处理装置，包括：

预处理模块，用于对地震数据进行叠前预处理，获得叠后地震数据；

补缺处理模块，用于对所述叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；

旋转模块，用于对所述补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；

压制模块，用于对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；

反旋转模块，用于对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；

恢复处理模块，用于对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据。

再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述地震数据的采集脚印的处理方法的步骤。

又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述地震数据的采集脚印的处理方法的步骤。

本发明实施例提供的地震数据的采集脚印的处理方法及装置，对地震数据进行叠前预处理，获得叠后地震数据，对叠后地震数据进行补缺处理，获得补缺后的叠后地震数据，对补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据，对旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据，对压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据，对反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据，能够保持原始地震数据的波形特征，提高了地震数据的信噪比和分辨率。

附图说明

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

图1是本发明一实施例提供的地震数据的采集脚印的处理方法的流程示意图。

图2是本发明另一实施例提供的地震数据的采集脚印的处理方法的流程示意图。

图3是本发明又一实施例提供的地震数据的采集脚印的处理方法的流程示意图。

图4是本发明再一实施例提供的地震数据的采集脚印的处理方法的流程示意图。

图5是本发明一实施例提供的第1个时间切片数据的部分数据的示意图。

图6是本发明一实施例提供的第1个时间切片数据的部分数据的InLine方向上的积分结果的示意图。

图7是本发明一实施例提供的第1个时间切片数据的部分数据对应的积分时间切片数据的示意图。

图8是本发明一实施例提供的新的第1个时间切片数据的部分数据的示意图。

图9是本发明一实施例提供的采集脚印压制前的叠后地震数据的示意图。

图10是本发明一实施例提供的去噪后的叠后地震数据的示意图。

图11是本发明一实施例提供的采集脚印压制前的叠后地震数据的频率振幅谱的示意图。

图12是本发明一实施例提供的去噪后的叠后地震数据的频率振幅谱的示意图。

图13是本发明一实施例提供的采集脚印压制前的叠后地震数据的切片显示示意图；

图14是本发明一实施例提供的去噪后的叠后地震数据的切片显示示意图。

图15是本发明一实施例提供的地震数据的采集脚印的处理装置的结构示意图。

图16是本发明另一实施例提供的地震数据的采集脚印的处理装置的结构示意图。

图17是本发明又一实施例提供的地震数据的采集脚印的处理装置的结构示意图。

图18是本发明一实施例提供的电子设备的实体结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。

图1是本发明一实施例提供的地震数据的采集脚印的处理方法的流程示意图，如图1所示，本发明实施例提供的地震数据的采集脚印的处理方法，包括：

S101、对地震数据进行叠前预处理，获得叠后地震数据；

具体地，在野外进行地震数据的采集，获得地震数据。对所述地震数据进行叠前预处理，获得叠后地震数据。其中，所述叠前预处理包括但不限于对地震数据进行置观测系统、静校正、速度分析、能量补偿，噪声衰减、切除等方式，根据实际需要进行设置，本发明实施例不做限定。在叠前预处理的过程中，应按照正常的处理流程，避免由于叠前预处理的不合理引入新的人为痕迹及脚印。所述地震数据可以是共中心点道集(Common MiddlePoint，简称CMP)。本发明实施例提供的地震数据的采集脚印的处理方法的执行主体包括但不限于计算机。

S102、对叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；

具体地，在获得所述叠后地震数据之后，对所述叠后地震数据进行补缺处理，以补齐所述叠后地震数据在主测线(InLine)和联络测线(CrossLine)方向上缺失的地震道数据，获得补缺后的叠后地震数据，补缺后的叠后地震数据在InLine和CrossLine方向形成一个规则的矩阵。其中，缺失的地震道数据可以用零值补齐。

例如，由于观测系统分布不均匀，InLine方向线号与CrossLine方向线号同CMP号不一致，那么需要对叠后地震数据进行补缺处理，使补缺后的叠后地震数据在整个工区形成一个规则矩阵，便于后续进行数据处理。在补齐地震道数据的同时缺失地震道数据的InLine和CrossLine被记录，以方便后续进行数据的恢复。

S103、对补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；

具体地，在获得补缺后的叠后地震数据之后，可以对补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据，旋转后的叠后地震数据包括多个切片数据，所述切片数据可以是时间切片数据或者深度切片数据。旋转后的叠后地震数据可以存储到一个临时文件中，方便后续进行数据读取。

例如，补缺后的叠后地震数据表示为(InLine，CrossLine，Z)，Z可以为时间或者深度，补缺后的叠后地震数据旋转后得到的旋转后的叠后地震数据表示为(Z，InLine，CrossLine)，旋转后的叠后地震数据包括多个切片数据，每个切片数据对应一个Z值，以切片数据的形式进行数据读取，适合对海量数据的处理，能够提高数据的读取效率。

S104、对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；

具体地，在获得所述旋转后的叠后地震数据之后，可以对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据。其中，在本发明实施例中，所述切片数据为时间切片数据或者深度切片数据。

S105、对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；

具体地，在获得所述压制后的叠后地震数据之后，对所述压制后的叠后地震数据进行反旋转，数据旋转的方向与步骤S103旋转的方向相反，是步骤S103的逆过程，可以获得反旋转的叠后地震数据。

例如，所述压制后的叠后地震数据是切片的形式(Z，InLine，CrossLine)，反旋转后得到的反旋转的叠后地震数据表示为(InLine，CrossLine，Z)，成为排列分布的数据。

S106、对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据；

具体地，在获得所述反旋转的叠后地震数据之后，对所述反旋转的叠后地震数据进行恢复处理，从所述反旋转的叠后地震数据中去除步骤S102补齐的地震道数据对应的数据，获得恢复后的叠后地震数据。

本发明实施例提供的地震数据的采集脚印的处理方法，对地震数据进行叠前预处理，获得叠后地震数据，对叠后地震数据进行补缺处理，获得补缺后的叠后地震数据，对补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据，对旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据，对压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据，对反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据，能够保持原始地震数据的波形特征，并压制原始地震数据中的采集脚印，提高了地震数据的信噪比和分辨率。

图2是本发明另一实施例提供的地震数据的采集脚印的处理方法的流程示意图，如图2所示，在上述各实施例的基础上，进一步地，所述对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据包括：

S1041、沿主测线方向对所述切片数据的振幅进行积分，获得主测线方向上的积分切片数据；

具体地，对于所述旋转叠后地震数据的每个切片数据，可以沿主测线方向对所述切片数据的振幅进行积分，积分的起点为最小InLine数，可以获得InLine方向上的积分切片数据。

S1042、沿联络测线方向对所述主测线方向上的积分切片数据的振幅进行积分，获得积分切片数据；

具体地，在获得InLine方向上的积分切片数据之后，沿联络测线方向对所述主测线方向上的积分切片数据的振幅进行积分，获得积分切片数据。其中，上述积分切片数据可以为积分时间偏移数据或者积分深度偏移数据。

S1043、根据设定主测线填充长度在主测线方向上对所述积分切片数据填充切片样点值，并根据设定联络测线填充长度在联络测线方向上对所述积分切片数据填充切片样点值，获得填充后的积分切片数据；

具体地，在获得所述积分切片数据之后，可以根据设定主测线填充长度在主测线方向上对所述积分切片数据填充切片样点值，即沿所述积分切片数据的每个InLine最小CMP方向填充主测线填充长度个切片样点值，填充长度为主测线填充长度，同样沿所述积分切片数据的每个InLine最大CMP方向填充主测线填充长度个切片样点值，填充长度同样为主测线填充长度，并根据设定联络测线填充长度在联络测线方向上对所述积分切片数据填充切片样点值，即沿所述积分切片数据的每个CrossLine最小line方向填充联络测线填充长度个切片样点值，填充长度为联络测线填充长度，同样沿切片的每个CrossLine最大line方向填充联络测线填充长度个样点值，填充长度为同样为联络测线填充长度，获得填充后的积分切片数据。其中，所述主测线填充长度和所述联络测线填充长度是预先设定的，根据实际需要进行设置，本发明实施例不做限定。

例如，对于所述积分切片数据，设定InLine方向和CrossLine方向两个空间变化步长，分别为InLine方向步长a和CrossLine方向步长b，a和b可以相等，也可以不相等，通常采集脚印能量分布越强的方向，选择的步长应该更长，采集脚印能量强弱分布可以通过积分切片数据观察获得。设定主测线填充长度为c，c为InLine方向步长的一半，即c＝a/2，设定联络测线填充长度为d，d为CrossLine方向步长b的一半，即d＝b/2。在对所述积分切片数据进行填充时，沿所述积分切片数据的每个InLine最小CMP方向填充c个切片样点值，填充长度为c，同样沿切片的每个InLine最大CMP方向填充c个切片样点值，填充长度同样为c，并沿切片的每个CrossLine最小line方向填充d个切片样点值，填充长度为d，同样沿切片的每个CrossLine最大line方向填充d个切片样点值，填充长度为同样为d，获得填充后的积分切片数据。其中，填充的样点值可以采用零值。

S1044、根据所述填充后的积分切片数据对所述积分切片数据中的每个样点数据进行对角振幅处理，获得积分样点振幅差值切片数据。

具体地，在获得所述填充后的积分切片数据之后，可以根据所述填充后的积分切片数据对所述积分切片数据中的每个样点数据进行对角振幅处理，即对于所述积分切片数据中的每个样点数据都可以从所述填充后的积分切片数据中获得对应的四个矩形角点数据，对四个矩形角点数据的振幅进行数据处理，可以获得积分样点振幅差值切片数据。对所述旋转叠后地震数据的每个切片数据，都进行步骤S1041、S1042、S1043和S1044的数据处理过程，即可获得压制后的叠后地震数据。

图3是本发明又一实施例提供的地震数据的采集脚印的处理方法的流程示意图，如图3所示，在上述各实施例的基础上，进一步地，所述根据所述填充后的积分切片数据对所述积分切片数据中的每个样点数据进行对角振幅处理，获得积分样点振幅差值切片数据包括：

S10441、从所述积分切片数据中获取一个样点数据，并从所述填充后的积分切片数据中获取所述样点数据对应的矩形角点；其中，所述样点数据对应的矩形角点有四个，构成矩形，所述样点数据位于矩形的中心；

具体地，从所述积分切片数据中获取一个样点数据，在所述填充后的积分切片数据中可以获得所述样点数据对应的矩形角点，所述样点数据对应的矩形角点有四个，四个矩形角点相连可以构成矩形，所述样点数据位于矩形的中心。

例如，在所述积分切片数据中获取任意一个样点数据，以这个样点数据为中心，可以在所述填充后的积分切片数据中上找到一个边长分别为a和b的一个矩形，并且这个样点数据恰好是矩形的中心。

S10442、计算第一矩形角点和相对的第三矩形角点对应的积分样点振幅值之和与第二矩形角点和相对的第四矩形角点对应的积分样点振幅值之和的差值，获得中间值；其中，所述样点数据对应的矩形角点包括所述第一矩形角点、所述第二矩形角点、所述第三矩形角点和所述第四矩形角点；

具体地，所述样点数据对应的矩形角点包括第一矩形角点、第二矩形角点、第三矩形角点和第四矩形角点，所述第一矩形角点和所述第三矩形角点相对，为矩形对角线的两个端点，所述第二矩形角点和所述第四矩形角点相对，为矩形另一个对角线的两个端点。计算所述第一矩形角点的积分样点振幅值与所述第三矩形角点的积分样点振幅值之和p，并计算第二矩形角点的积分样点振幅值与所述第四矩形角点的积分样点振幅值之和q，再计算q减去q的差值，将上述差值作为中间值。其中，所述第一矩形角点是积分InLine方向和CrossLine方向的相对起点，是所述样点数据对应的矩形角点中积分次数最少的那个点。

S10443、根据所述中间值、所述主测线填充长度、所述联络测线填充长度以及预设公式，计算获得所述样点数据对应的新样点数据。

具体地，在获得所述中间值之后，将所述中间值、所述主测线填充长度、所述联络测线填充长度带入到预设公式中，可以计算获得所述样点数据对应的新样点数据。对于所述积分切片数据中的每个样点数据，都进行步骤S10441、S10442和S10443的处理过程，获得每个样点数据的新样点数据，各个样点数据的新样点数据构成压制后的叠后地震数据。

例如，所述预设公式为E＝M/[(2c-1)(2d-1)]，其中，E为样点数据对应的新样点数据，M为中间值，c为主测线填充长度，d为联络测线填充长度。

在上述各实施例的基础上，进一步地，本发明实施例提供的地震数据的采集脚印的处理方法还包括：

根据所述叠后地震数据和恢复后的叠后地震数据进行混波组合处理，获得去噪后的叠后地震数据。

具体地，在获得所述恢复后的叠后地震数据之后，根据所述叠后地震数据和恢复后的叠后地震数据进行混波组合处理，可以获得去噪后的叠后地震数据。

图4是本发明再一实施例提供的地震数据的采集脚印的处理方法的流程示意图，如图4所示，在上述各实施例的基础上，进一步地，所述根据所述叠后地震数据和恢复后的叠后地震数据进行混波组合处理，获得去噪后的叠后地震数据包括：

S401、获取每个样点位置对应的混波百分比；

具体地，由于地震数据的采集脚印在浅层比较明显，随着地层的深度增加，采集脚印会变弱，因此可以为不同深度的地震数据设置不同的混波百分比，深度越深混波百分比越大。所述叠后地震数据通常是以深度或者时间表示的，相同的深度或者时间对应相同的混波百分比，每个样点位置对应深度或者时间，可以获取每个样点位置对应的混波百分比。

例如，设置所述叠后地震数据中混波百分比为1和0的样点位置，然后通过插值的方式获取混波百分比为1和百分比0的样点位置之间的各个样点位置对应的混波百分比。

S402、根据每个样点位置对应的混波百分比计算所述叠后地震数据和恢复后的叠后地震数据的样点数据的混波组合结果，作为去噪后的叠后地震数据的样点数据。

具体地，在获得每个样点位置的混波百分比之后，根据每个样点位置对应的混波百分比计算所述叠后地震数据和恢复后的叠后地震数据的样点数据的混波组合结果，即计算所述叠后地震数据和所述恢复后的叠后地震数据中相同样点位置的样点数据加权平均值的结果，作为去噪后的叠后地震数据的相同样点位置的样点数据。

例如，样点位置X对应的混波百分比为g，所述叠后地震数据在样点位置X的样点数据为m，所述恢复后的叠后地震数据在样点位置X的样点数据为n，那么去噪后的叠后地震数据在样点位置X的样点数据为m(1-g)+ng。

下面以一个具体的实施例来说明本发明实施例提供的地震数据的采集脚印的处理方法的实现过程。

对采集后地震数据进行叠前预处理。对地震数据进行置观测系统、静校正、速度分析、能量补偿和噪声衰减处理，获得叠后地震数据。叠前预处理的目的是要避免产生新的采集脚印，也可以减弱甚至压制原有数据采集脚印，最终为叠后采集脚印压制创造条件。

对叠后地震数据进行数据补缺处理，补齐叠后地震数据在InLine和CrossLine方向上缺失的地震道数据，获得补缺后的叠后地震数据，补缺后的叠后地震数据在InLine和CrossLine方向上形成一个规则矩阵，以便于后续进行采集脚印的压制。

对补缺后的叠后地震数据进行旋转，把(InLine，CrossLine，T)排列数据转换成以(T，InLine，CrossLine)排列分布的数据，T表示时间，这样处理后的主要目的是方便海量数据切片处理，防止数据量太大，程序内存不足的问题。将旋转后的叠后地震数据以时间切片的方式写入一个临时文件，

读入以切片方式保存的临时文件中旋转后的叠后地震数据的第1个时间切片数据，第1个时间切片数据的部分数据如图5所示，在第1个时间切片数据上首先沿每个InLine方向分别对第1个时间切片数据上的振幅进行积分，积分起点为最小InLine数，形成InLine方向上的积分时间切片数据。图5所示的第1个时间切片数据的部分数据的InLine方向上的积分结果如图6所示，边缘上的零值是填充零值后的积分值，图中只显示一行或者一列。然后沿CrossLine方向对InLine方向上的积分时间切片数据上的振幅再次进行积分，积分起点为最小CrossLine数，获得第1个时间切片数据对应的积分时间切片数据。图6所示的InLine方向上的积分时间切片数据，沿CrossLine方向的积分结果如图7所示。

设定InLine方向和CrossLine方向两个空间变化步长，分别为InLine方向步长6和CrossLine方向步长6。两个空间步长分别除以2，获得主测线填充长度为3，联络测线填充长度为3。沿第1个时间切片数据对应的积分时间切片数据的每个InLine最小CMP方向填充3个切片样点值，填充长度为3，同样沿第1个时间切片数据对应的积分时间切片数据的每个InLine最大CMP方向填充3个切片样点值，填充长度同样为3，边界扩展填充是本发明实施例考虑采集脚印滤波边界效应对采集脚印压制带来的影响，都用零填充。同样沿第1个时间切片数据对应的积分时间切片数据的每个CrossLine最小line方向填充3个切片样点值，填充长度为3，同样沿第1个时间切片数据对应的积分时间切片数据的每个CrossLine最大line方向填充3个切片样点值，填充长度为同样为3。经过上述数据填充处理，可以获得第1个时间切片数据对应的填充后的积分切片数据。

再次读入第1个时间切片数据，在第1个时间切片数据上获取任意一个样点数据，以这个样点数据为中心，可以在第1个时间切片数据对应的填充后的积分切片数据上找到一个边长分别为3和3的矩形，如图5中第5列第5行的样点数据为3.1。图5中的第5列第5行的样点在第1个时间切片数据对应的填充后的积分切片数据上对应第3行第3列、第8行第8列、第8行第3列和第3行第8列的中心，如图7所示，第5.5行第5.5列位置的黑色虚线方框，因为第1个时间切片数据对应的填充后的积分切片数据的样点值的位置与第1个时间切片数据的样点值的位置相差0.5列和0.5行，以这个黑色虚线方框为矩形中心，形成以长和宽为6的矩形，用矩形的第一矩形角点上对应的第3行第3列的样点振幅值10.5加上其对角上的第8行第8列的样点振幅值126.8，再减去另外两个对角上第3行第8列的样点振幅值39.2与第8行第3列的样点振幅值36.9的和，获得以长和宽为6矩形的对角积分样点振幅差值，即样点数据3.1对应的中间值，中间值大小为126.8+10.5-39.2-36.9＝61.2。用61.2/(6-1)/(6-1)＝2.4，用2.4替换第1个时间切片数据上第5行第5列上的样点数据3.1，就获得样点数据3.1对应的新样点数据2.4。第1个时间切片数据的所有样点数据都用上述方法进行替换就形成一个新的第1个时间切片数据，新的第1个时间切片数据的部分数据如图8所示。

用一个新的临时文件保存好这个新的第1个时间切片数据。然后采用对第1个时间切片数据同样的处理方法，依次处理第2个时间切片数据、第3个时间切片数据、第4个时间切片数据，直到旋转后的叠后地震数据中的每个时间切片数据都处理完成，从而获得压制后的叠后地震数据。

从新的临时文件中读出压制后的叠后地震数据，进行反旋转，把切片方式(T，InLine，CrossLine)数据旋转成以(InLine，CrossLine，T)排列分布的数据，获得反旋转的叠后地震数据。从反旋转的叠后地震数据中剔除补缺处理时补齐的地震数据所对应的数据，获得恢复后的叠后地震数据。

最后根据采集脚印存在时间段上的严重程度预测100ms处混波百分比为10％，1000ms处为60％，5000ms为98％，并对这些不同时间点的混波百分比进行插值，保证每个地震道样点上都有一个对应混波百分比。读出叠后地震数据和恢复后的叠后地震数据进行混波组合处理，获得去噪后的叠后地震数据，去噪后的叠后地震数据就是去除采集脚印后的地震数据。

图5和图8中的数据去噪前后采集脚印效果图可以对比图9至图14。对比图9和图10，可以发现采集脚印压制后同相轴更光滑，信噪比分辨率得到显著提高。对比图11和图12可以发现采集脚印压制后频谱上主要能量基本保持一致，只是40Hz-80Hz的频谱能量变化较大，说明采集脚印能量主要分布在40Hz-80Hz范围内，同时也说明本发明对采集脚印能量在InLinehe和CrossLine方向分布趋势以及在频率域分布范围判断准确，采集脚印压制后保振幅能力较强。对比图13和图14可以发现经过本发明实施例提供的地震数据的采集脚印的处理方法对采集脚印压制后，采集脚印切片更光滑，在图14中已经很难用肉眼识别到采集脚印，说明本发明实施例提供的地震数据的采集脚印的处理方法提高了地震数据时间切片上信噪比和分辨率，消除了与油气藏预测无关的采集脚印干扰，为下一步叠后解释，油气藏成功预测创造有利条件。

图15是本发明一实施例提供的地震数据的采集脚印的处理装置的结构示意图，如图15所示，本发明实施例提供的地震数据的采集脚印的处理装置，包括预处理模块1510、补缺处理模块1520、旋转模块1530、压制模块1540、反旋转模块1550和恢复处理模块1560，其中：

预处理模块1510用于对地震数据进行叠前预处理，获得叠后地震数据；补缺处理模块1520用于对所述叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；旋转模块1530用于对所述补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；压制模块1540用于对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；反旋转模块1550用于对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；恢复处理模块1560用于对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据。

具体地，在野外进行地震数据的采集，获得地震数据。预处理模块1510对所述地震数据进行叠前预处理，获得叠后地震数据。其中，所述叠前预处理包括但不限于对地震数据进行置观测系统、静校正、速度分析、能量补偿，噪声衰减、切除等方式，根据实际需要进行设置，本发明实施例不做限定。在叠前预处理的过程中，应按照正常的处理流程，避免由于叠前预处理的不合理引入新的人为痕迹及脚印。

在获得所述叠后地震数据之后，补缺处理模块1520对所述叠后地震数据进行补缺处理，以补齐所述叠后地震数据在主测线(InLine)和联络测线(CrossLine)方向上缺失的地震道数据，获得补缺后的叠后地震数据，补缺后的叠后地震数据在InLine和CrossLine方向形成一个规则的矩阵。其中，缺失的地震道数据可以用零值补齐。

在获得补缺后的叠后地震数据之后，旋转模块1530可以对补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据，旋转后的叠后地震数据包括多个切片数据，所述切片数据可以是时间切片数据或者深度切片数据。旋转后的叠后地震数据可以存储到一个临时文件中，方便后续进行数据读取。

在获得所述旋转后的叠后地震数据之后，压制模块1540可以对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据。其中，在本发明实施例中，所述切片数据为时间切片数据或者深度切片数据。

在获得所述压制后的叠后地震数据之后，反旋转模块1550对所述压制后的叠后地震数据进行反旋转，可以获得反旋转的叠后地震数据。

在获得所述反旋转的叠后地震数据之后，恢复处理模块1560对所述反旋转的叠后地震数据进行恢复处理，从所述反旋转的叠后地震数据中去除补齐的地震道数据对应的数据，获得恢复后的叠后地震数据。

本发明实施例提供的地震数据的采集脚印的处理装置，对地震数据进行叠前预处理，获得叠后地震数据，对叠后地震数据进行补缺处理，获得补缺后的叠后地震数据，对补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据，对旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据，对压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据，对反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据，能够保持原始地震数据的波形特征，并压制原始地震数据中的采集脚印，提高了地震数据的信噪比和分辨率。

图16是本发明另一实施例提供的地震数据的采集脚印的处理装置的结构示意图，如图16所示，在上述各实施例的基础上，进一步地，压制模块1540包括第一积分单元1541、第二积分单元1542、填充单元1543和对角振幅处理单元1544，其中：

第一积分单元1541用于沿主测线方向对所述切片数据的振幅进行积分，获得主测线方向上的积分切片数据；第二积分单元1542用于沿联络测线方向对所述主测线方向上的积分切片数据的振幅进行积分，获得积分切片数据；填充单元1543用于根据设定主测线填充长度在主测线方向上对所述积分切片数据填充切片样点值，并根据设定联络测线填充长度在联络测线方向上对所述积分切片数据填充切片样点值，获得填充后的积分切片数据；对角振幅处理单元1544用于根据所述填充后的积分切片数据对所述积分切片数据中的每个样点数据进行对角振幅处理，获得压制后的叠后地震数据。

具体地，对于所述旋转叠后地震数据的每个切片数据，第一积分单元1541可以沿主测线方向对所述切片数据的振幅进行积分，积分的起点为最小InLine数，可以获得InLine方向上的积分切片数据。

在获得InLine方向上的积分切片数据之后，第二积分单元1542沿联络测线方向对所述主测线方向上的积分切片数据的振幅进行积分，获得积分切片数据。其中，上述积分切片数据可以为积分时间偏移数据或者积分深度偏移数据。

在获得所述积分切片数据之后，填充单元1543可以根据设定主测线填充长度在主测线方向上对所述积分切片数据填充切片样点值，即沿所述积分切片数据的每个InLine最小CMP方向填充主测线填充长度个切片样点值，填充长度为主测线填充长度，同样沿所述积分切片数据的每个InLine最大CMP方向填充主测线填充长度个切片样点值，填充长度同样为主测线填充长度，并根据设定联络测线填充长度在联络测线方向上对所述积分切片数据填充切片样点值，即沿所述积分切片数据的每个CrossLine最小line方向填充联络测线填充长度个切片样点值，填充长度为联络测线填充长度，同样沿切片的每个CrossLine最大line方向填充联络测线填充长度个样点值，填充长度为同样为联络测线填充长度，获得填充后的积分切片数据。其中，所述主测线填充长度和所述联络测线填充长度是预先设定的，根据实际需要进行设置，本发明实施例不做限定。

在获得所述填充后的积分切片数据之后，对角振幅处理单元1544可以根据所述填充后的积分切片数据对所述积分切片数据中的每个样点数据进行对角振幅处理，即对于所述积分切片数据中的每个样点数据都可以从所述填充后的积分切片数据中获得对应的四个矩形角点数据，对四个矩形角点数据的振幅进行数据处理，可以获得积分样点振幅差值切片数据。

图17是本发明又一实施例提供的地震数据的采集脚印的处理装置的结构示意图，如图17所示，在上述各实施例的基础上，进一步地，本发明实施例提供的填充单元1543包括获取子单元15431、第一计算子单元15432和第二计算子单元15433，其中：

获取子单元15431用于从所述积分切片数据中获取一个样点数据，并从所述填充后的积分切片数据中获取所述样点数据对应的矩形角点；其中，所述样点数据对应的矩形角点有四个，构成矩形，所述样点数据位于矩形的中心；第一计算子单元15432用于计算第一矩形角点和相对的第三矩形角点对应的积分样点振幅值之和与第二矩形角点和相对的第四矩形角点对应的积分样点振幅值之和的差值，获得中间值；其中，所述样点数据对应的矩形角点包括所述第一矩形角点、所述第二矩形角点、所述第三矩形角点和所述第四矩形角点；第二计算子单元15433用于根据所述中间值、所述主测线填充长度、所述联络测线填充长度以及预设公式，计算获得所述样点数据对应的新样点数据。

具体地，获取子单元15431从所述积分切片数据中获取一个样点数据，在所述填充后的积分切片数据中可以获得所述样点数据对应的矩形角点，所述样点数据对应的矩形角点有四个，四个矩形角点相连可以构成矩形，所述样点数据位于矩形的中心。

所述样点数据对应的矩形角点包括第一矩形角点、第二矩形角点、第三矩形角点和第四矩形角点，所述第一矩形角点和所述第三矩形角点相对，为矩形对角线的两个端点，所述第二矩形角点和所述第四矩形角点相对，为矩形另一个对角线的两个端点。第一计算子单元15432计算所述第一矩形角点的积分样点振幅值与所述第三矩形角点的积分样点振幅值之和p，并计算第二矩形角点的积分样点振幅值与所述第四矩形角点的积分样点振幅值之和q，再计算q减去q的差值，将上述差值作为中间值。其中，所述第一矩形角点是积分InLine方向和CrossLine方向的相对起点，是所述样点数据对应的矩形角点中积分次数最少的那个点。

在获得所述中间值之后，第二计算子单元15433将所述中间值、所述主测线填充长度、所述联络测线填充长度带入到预设公式中，可以计算获得所述样点数据对应的新样点数据。

本发明实施例提供的装置的实施例具体可以用于执行上述各方法实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。

图18是本发明一实施例提供的电子设备的实体结构示意图，如图18所示，该电子设备可以包括：处理器(processor)1801、通信接口(Communications Interface)1802、存储器(memory)1803和通信总线1804，其中，处理器1801，通信接口1802，存储器1803通过通信总线1804完成相互间的通信。处理器1801可以调用存储器1803中的逻辑指令，以执行如下方法：对地震数据进行叠前预处理，获得叠后地震数据；对所述叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；对所述补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据。

此外，上述的存储器1803中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：对地震数据进行叠前预处理，获得叠后地震数据；对所述叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；对所述补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据。

本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：对地震数据进行叠前预处理，获得叠后地震数据；对所述叠后地震数据进行补缺处理，获得补缺后的叠后地震数据；对所述补缺后的叠后地震数据进行旋转，获得旋转后的叠后地震数据；对所述旋转后的叠后地震数据的每个切片数据进行采集脚印压制，获得压制后的叠后地震数据；对所述压制后的叠后地震数据进行反旋转，获得反旋转的叠后地震数据；对所述反旋转的叠后地震数据进行恢复处理，获得恢复后的叠后地震数据。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。