一种绕射多次波压制方法及装置

摘要

本发明提供了一种绕射多次波压制方法及装置，涉及地震资料多次波压制的技术领域，包括在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；对所述剩下波场数据进行反动校正处理，得到波场数据；将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。解决了先压制反射波后压制绕射多次波导致的双曲形态被破坏的绕射多次波无法被压制的技术问题，达到有效地使绕射多次波得到压制，提高压制绕射多次波的效率的技术效果。

权利要求

1.一种绕射多次波压制方法，其特征在于，包括以下步骤：在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；

对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；

对所述剩下波场数据进行反动校正处理，得到波场数据；

将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波；

所述在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据，包括：利用与所述原始地震数据对应的目标区域的深度计算所述目标区域中地震波的叠加速度；

根据所述叠加速度将所述原始地震数据进行动校正处理，得到动校正结果；

将所述动校正结果按照预设的偏移距分选出共偏移距道集；

在所述共偏移距道集中根据所述动校正结果选取含有所述绕射多次波的时窗；

利用平面波解构方法在所述时窗内提取所述反射波数据，得到剩余地震数据。

2.根据权利要求1所述的绕射多次波压制方法，其特征在于，所述对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据，包括：利用预设的绕射走时计算公式，计算绕射多次走时；

根据所述绕射多次走时对所述剩余地震数据进行动校正处理，得到校平数据；

在所述校平数据中，利用所述平面波解构方法去除绕射多次波数据，得到所述剩下波场数据。

3.根据权利要求1或2所述的绕射多次波压制方法，其特征在于，所述对所述剩下波场数据进行反动校正处理，得到波场数据，包括：基于所述偏移距、所述叠加速度及预设的时差计算公式，计算时差；

根据所述时差对所述剩下波场数据进行反动校正处理，得到波场数据。

4.一种绕射多次波压制装置，其特征在于，包括：提取模块，用于在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；

动校正处理模块，用于对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；

反动校正处理模块，用于对所述剩下波场数据进行反动校正处理，得到波场数据；

组合模块，用于将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波；

所述提取模块，包括：

第一计算单元，用于利用与所述原始地震数据对应的目标区域的深度计算所述目标区域中地震波的叠加速度；

第一动校正处理单元，用于根据所述叠加速度将所述原始地震数据进行动校正处理，得到动校正结果；

分选单元，用于将所述动校正结果按照预设的偏移距分选出共偏移距道集；

选取单元，用于在所述共偏移距道集中根据所述动校正结果选取含有所述绕射多次波的时窗；

提取单元，用于利用平面波解构方法在所述时窗内提取所述反射波数据，得到剩余地震数据。

5.根据权利要求4所述的绕射多次波压制装置，其特征在于，所述动校正处理模块，包括：第二计算单元，用于利用预设的绕射走时计算公式，计算绕射多次走时；

第二动校正处理单元，用于根据所述绕射多次走时对所述剩余地震数据进行动校正处理，得到校平数据；

去除单元，用于在所述校平数据中，利用所述平面波解构方法去除绕射多次波数据，得到所述剩下波场数据。

6.根据权利要求4或5所述的绕射多次波压制装置，其特征在于，所述反动校正处理模块，包括：第三计算单元，用于基于所述偏移距、所述叠加速度及预设的时差计算公式，计算时差；

反动校正处理单元，用于根据所述时差对所述剩下波场数据进行反动校正处理，得到波场数据。

7.一种电子设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至3任一项所述方法的步骤。

8.一种具有处理器可执行的非易失的程序代码的计算机可读介质，其特征在于，所述程序代码使所述处理器执行所述权利要求1-3任一所述方法。

说明书

一种绕射多次波压制方法及装置

技术领域

[0001] 本发明涉及地震资料多次波压制技术领域，尤其是涉及一种绕射多次波压制方法及装置。

背景技术

[0002] 多次波是地震波研究中会出现的一种干扰波的种类，在探地雷达中也常出现。多次波在海洋地震资料中一般都异常发育，其特点为分布范围广、能量强。根据多次波产生机理的不同可以将不同类型的多次波进行有针对性的压制。

[0003] 目前的多次波压制技术，一般都是先压制反射多次波，然后再对一些剩余的多次波进行处理。这些剩余的多次波有一大部分为崎岖海底或者复杂构造引起的绕射多次波，有一小部分是之前没有去除干净的反射多次波。

[0004] 剩余的绕射多次波往往就是一些随机剩余能量和频率等。常规的剩余绕射多次波去除方法基本上是在常规压制反射多次波(压制过程使用自由表面多次波衰减(Surface-Related Multiple Elimination，SRME)与高精度Radon变换)之后，然后再压制残余的绕射多次波。

[0005] 由于常规的反射多次波压制方法根据反射多次波的形态、传播规律等对具有反射多次波形态的多次波进行了压制，同时对绕射多次波的基本形态会有一定的破坏作用，很难保留绕射多次波的一些基本特征。因此，在压制反射多次波的同时，会破坏绕射多次波在空间和时间上的形态规律，使绕射多次波的双曲形态不完整，导致双曲形态被破坏的绕射多次波无法被压制，降低压制绕射多次波的效率。

发明内容

[0006] 有鉴于此，本发明的目的在于提供一种绕射多次波压制方法及装置，以缓解现有技术中存在的先压制反射多次波再压制绕射多次波，破坏了绕射多次波的双曲形态，使绕射多次波的双曲形态不完整，导致双曲形态被破坏的绕射多次波无法被压制，降低压制绕射多次波的效率的技术问题。

[0007] 第一方面，本发明实施例提供了一种绕射多次波压制方法，包括以下步骤：

[0008] 在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；

[0009] 对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；

[0010] 对所述剩下波场数据进行反动校正处理，得到波场数据；

[0011] 将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。

[0012] 结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，所述在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据，包括：

[0013] 利用与所述原始地震数据对应的目标区域的深度计算所述目标区域中地震波的叠加速度；

[0014] 根据所述叠加速度将所述原始地震数据进行动校正处理，得到动校正结果；

[0015] 将所述动校正结果按照预设的偏移距分选出共偏移距道集；

[0016] 在所述共偏移距道集中根据所述动校正结果选取含有所述绕射多次波的时窗；

[0017] 利用平面波解构方法在所述时窗内提取所述反射波数据，得到剩余地震数据。

[0018] 结合第一方面，本发明实施例提供了第一方面的第二种可能的实施方式，其中，所述对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据，包括：

[0019] 利用预设的绕射走时计算公式，计算绕射多次走时；

[0020] 根据所述绕射多次走时对所述剩余地震数据进行动校正处理，得到校平数据；

[0021] 在所述校平数据中，利用所述平面波解构方法去除绕射多次波数据，得到所述剩下波场数据。

[0022] 结合第一方面，本发明实施例提供了第一方面的第三种可能的实施方式，其中，对所述剩下波场数据进行反动校正处理，得到波场数据，包括：

[0023] 基于所述偏移距、所述叠加速度及预设的时差计算公式，计算时差；

[0024] 根据所述时差对所述剩下波场数据进行反动校正处理，得到波场数据。

[0025] 第二方面，本发明实施例还提供一种绕射多次波压制装置，包括：

[0026] 提取模块，用于在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；

[0027] 动校正处理模块，用于对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；

[0028] 反动校正处理模块，用于对所述剩下波场数据进行反动校正处理，得到波场数据；

[0029] 组合模块，用于将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。

[0030] 结合第二方面，本发明实施例提供了第二方面的第一种可能的实施方式，其中，所述提取模块，包括：

[0031] 第一计算单元，用于利用与所述原始地震数据对应的目标区域的深度计算所述目标区域中地震波的叠加速度；

[0032] 第一动校正处理单元，用于根据所述叠加速度将所述原始地震数据进行动校正处理，得到动校正结果；

[0033] 分选单元，用于将所述动校正结果按照预设的偏移距分选出共偏移距道集；

[0034] 选取单元，用于在所述共偏移距道集中根据所述动校正结果选取含有所述绕射多次波的时窗；

[0035] 提取单元，用于利用平面波解构方法在所述时窗内提取所述反射波数据，得到剩余地震数据。

[0036] 结合第二方面，本发明实施例提供了第二方面的第二种可能的实施方式，其中，所述动校正处理模块，包括：

[0037] 第二计算单元，用于利用预设的绕射走时计算公式，计算绕射多次走时；

[0038] 第二动校正处理单元，用于根据所述绕射多次走时对所述剩余地震数据进行动校正处理，得到校平数据；

[0039] 去除单元，用于在所述校平数据中，利用所述平面波解构方法去除绕射多次波数据，得到所述剩下波场数据。

[0040] 结合第二方面，本发明实施例提供了第二方面的第三种可能的实施方式，其中，反动校正处理模块，包括：

[0041] 第三计算单元，用于基于所述偏移距、所述叠加速度及预设的时差计算公式，计算时差；

[0042] 反动校正处理单元，用于根据所述时差对所述剩下波场数据进行反动校正处理，得到波场数据。

[0043] 第三方面，本发明实施例还提供一种电子设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，其中，所述处理器执行所述计算机程序时实现第一方面所述的方法。

[0044] 第四方面，本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，其中，所述程序代码使所述处理器执行第一方面所述方法。

[0045] 本发明实施例带来了以下有益效果：本发明实施例通过首先在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；然后对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；再对所述剩下波场数据进行反动校正处理，得到波场数据；最后可以将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。

[0046] 本发明实施例能够通过先从原始地震数据中提取出反射波数据，得到不包含反射波数据的剩余地震数据(剩余地震数据中包含绕射波数据)，通过对剩余地震数据的动校正处理及反动校正处理，得到仅包含反射多次波数据的波场数据(相当于压制了绕射多次波)，也即能够在压制反射多次波之前先压制绕射多次波，保证绕射多次波双曲形态完整，避免先压制反射多次波过程中绕射多次波的双曲形态被破坏的情况，所以，本发明实施例可以有效地使绕射多次波得到压制，提高压制绕射多次波的效率。

[0047] 本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。

[0048] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。

附图说明

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

[0050] 图1为本发明实施例提供的一种绕射多次波压制方法的流程图；

[0051] 图2为图1中步骤S101的流程图；

[0052] 图3为图1中步骤S102的流程图；

[0053] 图4为本发明实施例提供的一种绕射多次波压制装置的结构图。

[0054] 图标：11-提取模块；12-动校正处理模块；13-反动校正处理模块；14-组合模块。

具体实施方式

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

[0056] 发明人发现，目前多次波压制技术，一般都是先压制反射多次波，然后再对一些剩余的多次波进行处理，压制反射多次波的过程使绕射多次波的双曲形态不完整，导致双曲形态被破坏的绕射多次波无法被压制，降低压制绕射多次波的效率。基于此，本发明实施例提供的一种绕射多次波压制方法及装置，可以通过先从原始地震数据中提取出反射波数据，得到不包含反射波数据的剩余地震数据(剩余地震数据中包含绕射波数据)，通过对剩余地震数据的动校正处理及反动校正处理，得到仅包含反射多次波数据的波场数据(相当于压制了绕射多次波)，也即能够在压制反射多次波之前先压制绕射多次波，保证绕射多次波双曲形态完整，避免先压制反射多次波过程中绕射多次波的双曲形态被破坏的情况，所以，本发明实施例可以有效地使绕射多次波得到压制，提高压制绕射多次波的效率。

[0057] 为便于对本实施例进行理解，首先对本发明实施例所公开的一种绕射多次波压制方法进行详细介绍，如图1所示，所述绕射多次波压制方法可以包括以下步骤：

[0058] 步骤S101，在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；

[0059] 在本发明实施例中，由于反射波具有曲率小和平滑度较高的特性，所以可以利用平面波结构方法提取反射波数据，最终在提取出反射波数据之后得到反射波数据及不包含反射波数据的剩余地震数据。在本发明实施例中，可以利用平面波解构滤波器提取反射波数据，平面波解构滤波器的局部平面波微分方程可以定义为：

[0060]

[0061] 式中，P(t,x)表示原始地震数据；σ为局部倾角，可随时间和空间变化。

[0062] 若倾角σ与时间无关，可将公式(1)转换到频率域，形式如下：

[0063]

[0064] 频率域局部平面波微分方程(2)，解的形式可表示为：

[0065]

[0066] 其中， 为P(t,x)的傅里叶表达式，指数中的复数项表示时移。若采用F-X域预测误差滤波器，则上式(3)可表示为：

[0067]

[0068] 其中，a0＝1,a1＝-eiωσ。

[0069] 对于倾角随时间变化的情况，可由上面的推导得出一定的结论。平面波传播特性是，传播过程中总能量恒定。该特性可由频率域公式(4)中复数指数振幅谱值为1得到进一步证实。利用Z变换，全通滤波形式可表示为：

[0070]

[0071] 其中， 为相应地震道的Z变换形式，B(Zt)/B(1/Zt)为全通数字滤波器，需要逼近时移项eiωσ。

[0072] 通过基于泰勒级数展开的推导,可以得出任意点数的中心滤波器，以三点中心滤波器B3(Zt)为例，形式如下：

[0073]

[0074] 对于低频信号，拟合滤波器随着高阶项的加入，精度会提高，对于二维情况，方程(6)，可以表示为：

[0075]

[0076] 为了避免多项式相除的情况，(7)式可以改进为如下形式：

[0077] C(Zt,Zx)＝A(Zt,Zx)B(1/Zt)＝B(1/Zt)-ZxB(Zt)  (8)

[0078] 将反射波数据C(Zt，Zx)转化为C(σ)，当满足条件：C(σ)≈0时，即可完成提取反射波数据，得到不含反射波的剩余地震数据。

[0079] 所述原始地震数据是由地震仪采集且未经数字处理技术处理的地震数据；

[0080] 所述剩余地震数据，包括绕射多次波和剩下波场数据两部分；

[0081] 所述反射波数据涉及到与反射波相关的各种数据，包含速度，位移，波长，振幅等物理量，但不仅限于上述给出的信息。

[0082] 步骤S102，对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；

[0083] 所述动校正处理是为了将绕射多次波数据拉平，使绕射多次波为近似水平的形态，进而可以分离绕射多次波和剩余波场，便于用平面波解构方法进行去除绕射多次波；

[0084] 示例性的，预设一串剩余地震数据有10、20、40、41、40、40、41、42、41、300和400等，动校正量预设为40，经动校正处理后的数据分别对应着-30、-20、0、1、0、0、1、2、1、260和360。显然，在经动校正处理后的数据中可以截取到“0、1、0、0、1、2和1”，这段数据显示的是绕射多次波近似水平的形态。在这里示出和描述的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为限制，因此，示例性实施例的其他示例可以具有不同的值。

[0085] 步骤S103，对所述剩下波场数据进行反动校正处理，得到波场数据；

[0086] 所述反动校正处理是将已经做过动校正的剩下波场数据还原成动校正前的剩下波场数据。

[0087] 所述步骤S103可以包括以下步骤：

[0088] 基于所述偏移距、所述叠加速度及预设的时差计算公式，计算时差；所述时差为t的计算公式为：

[0089]

[0090] 式中，x为检波点到成像中心点的距离，h为偏移距，v为叠加速度。

[0091] 根据所述时差对所述剩下波场数据进行反动校正处理，得到波场数据。示例性的，预设某地震数据为100、120、140和160，动校正量预设为40，经动校正处理后的剩下波长数据为60、80、100和120，反动校正根据计算得到的时差将数据还原为100、120、140和160。在这里示出和描述的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为限制，因此，示例性实施例的其他示例可以具有不同的值。

[0092] 步骤S104，将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。

[0093] 本发明实施例通过首先在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；然后对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；再对所述剩下波场数据进行反动校正处理，得到波场数据；最后可以将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。

[0094] 本发明实施例能够通过先从原始地震数据中提取出反射波数据，得到不包含反射波数据的剩余地震数据(剩余地震数据中包含绕射波数据)，通过对剩余地震数据的动校正处理及反动校正处理，得到仅包含反射多次波数据的波场数据(相当于压制了绕射多次波)，也即能够在压制反射多次波之前先压制绕射多次波，保证绕射多次波双曲形态完整，避免先压制反射多次波过程中绕射多次波的双曲形态被破坏的情况，所以，本发明实施例可以有效地使绕射多次波得到压制，提高压制绕射多次波的效率。

[0095] 在本发明的又一实施例中，如图2所示，所述步骤S101可以包括以下步骤：

[0096] 步骤S201，利用与所述原始地震数据对应的目标区域的深度计算所述目标区域中地震波的叠加速度；

[0097] 叠加速度是利用原始地震数据对应的目标区域的深度通过计算速度谱求得的速度。

[0098] 步骤S202，根据所述叠加速度将所述原始地震数据进行动校正处理，得到动校正结果；

[0099] 所述原始地震数据在得到叠加速度之后，经过动校正处理使反射波曲线拉成平行x轴的直线，即反射波数据为零，其他非零数据为动校正数据。

[0100] 所述动校正处理，是为了形成矩形窗，方便于绕射波的时窗选取。

[0101] 所述动校正数据数据是动校正处理后，得到的校平的剩余地震数据。

[0102] 步骤S203，将所述动校正结果按照预设的偏移距分选出共偏移距道集；

[0103] 所述共偏移距道集是按照同一个偏移距形成的道集。示例性的，若预设的偏移距为1和n，n为正整数，那么共偏移距道集将形成偏移距为1和n的两个道集。

[0104] 步骤S204，在所述共偏移距道集中根据所述动校正结果选取含有所述绕射多次波的时窗；

[0105] 所述时窗是矩形时窗；

[0106] 步骤S205，利用平面波解构方法在所述时窗内提取所述反射波数据，得到剩余地震数据。

[0107] 在本发明实施例中，可以通过选取含有绕射多次波的时窗，并只在视窗中提取发射多次波数据，得到不包含反射多次波数据的剩余地震数据。

[0108] 在本发明的又一实施例中，如图3所示，所述步骤S102可以包括以下步骤：

[0109] 步骤S301，利用预设的绕射走时计算公式，计算绕射多次走时；

[0110] 以海洋地震资料为例，所述预设的绕射走时计算公式为：

[0111]

[0112] 式中，t为绕射多次走时，x为检波点到成像中心点的距离，h为海底面深度，v为海水速度。

[0113] 步骤S302，根据所述绕射多次走时对所述剩余地震数据进行动校正处理，得到校平数据；

[0114] 步骤S303在所述校平数据中，利用所述平面波解构方法去除绕射多次波数据，得到所述剩下波场数据。

[0115] 通过本发明实施例，能够在剩余地震数据中有效去除绕射多次波数据，得到剩下波场数据。

[0116] 在本发明的又一实施例中，如图4所示，所述装置可以包括以下模块：

[0117] 提取模块11，用于在原始地震数据中提取反射波数据，得到反射波数据及不包含所述反射波数据的剩余地震数据；

[0118] 动校正处理模块12，用于对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平的剩下波场数据；

[0119] 反动校正处理模块13，用于对所述剩下波场数据进行反动校正处理，得到波场数据；

[0120] 组合模块14，用于将所述波场数据与所述反射波数据组合，得到压制所述绕射多次波之后的目标地震数据，以用于压制反射多次波。

[0121] 本发明实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。

[0122] 在本发明的又一实施例中，所述提取模块，包括：

[0123] 第一计算单元，用于利用与所述原始地震数据对应的目标区域的深度计算所述目标区域中地震波的叠加速度；

[0124] 第一动校正处理单元，用于根据所述叠加速度将所述原始地震数据进行动校正处理，得到动校正结果；

[0125] 分选单元，用于将所述动校正结果按照预设的偏移距分选出共偏移距道集；

[0126] 选取单元，用于在所述共偏移距道集中根据所述动校正结果选取含有所述绕射多次波的时窗；

[0127] 提取单元，用于利用平面波解构方法在所述时窗内提取所述反射波数据，得到剩余地震数据。

[0128] 在本发明的又一实施例中，所述动校正模块，包括：

[0129] 第二计算单元，利用预设的绕射走时计算公式，计算绕射多次走时；

[0130] 第二动校正处理单元，根据所述绕射多次走时对所述剩余地震数据进行动校正处理，以去除绕射多次波，得到校平数据；

[0131] 去除单元，在所述校平数据中，利用所述平面波解构方法去除绕射多次波数据，得到所述剩下波场数据。

[0132] 在本发明的又一实施例中，所述反动校正模块，包括：

[0133] 第三计算单元，用于基于所述偏移距、所述叠加速度及预设的时差计算公式，计算时差；

[0134] 反动校正处理单元，用于根据所述时差对所述剩下波场数据进行反动校正处理，得到波场数据。

[0135] 在本发明的又一实施例中，还提供一种电子设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法实施例所述方法的步骤。

[0136] 在本发明的又一实施例中，还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，所述程序代码使所述处理器执行方法实施例所述方法。

[0137] 附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

[0138] 本发明实施例所提供的绕射多次波压制方法的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。

[0139] 所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统和装置的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

[0140] 另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

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

[0142] 在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。

[0143] 最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。