CN103308943A - 一种海洋地震资料处理中层间多次波衰减的方法及装置 - Google Patents

Abstract

本发明提供了一种海洋地震资料处理中层间多次波衰减的方法及装置，包括：对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)；将动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置；在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)；根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并进行反动校，得到第一波场项P₁(x,z₀;t)；令第三波场项P₃(x,z₀;t)等于第一波场项；根据第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)；在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)；由于该方案通过纯数据驱动进行层间多次波预测和衰减，因而不依赖与速度场的精度，提高了海洋地震资料的处理质量。

Description

一种海洋地震资料处理中层间多次波衰减的方法及装置

技术领域

本发明涉及地震资料处理技术领域，特别涉及一种海洋地震资料处理中层间多次波衰减的方法及装置。

背景技术

随着油气勘探的深入，海洋油气勘探受到了越来越多的重视，海洋地震资料处理成果是海洋油气勘探的重要依据，而多次波是影响海洋地震资料处理的最为关键因素之一。由于海平面和海底都存在强速度差的反射界面，反射系数很大，所以表面相关多次波非常严重，严重的影响了地震资料处理的可靠性，降低了数据的信噪比，干扰了地震资料解释人员对构造的判断。在海上地震资料中，地表相关多次波是其主要形式，当前阶段已经有多种地表相关多次波衰减方法，其中应用最为广泛的就是基于反馈环的地表相关多次波衰减（SRME）方法，该方法已经在工业界得到了广泛的应用。

但是，当地下存在较强反射界面时，就会产生能量不容忽视的层间多次波，如海洋资料中海底与海底下盐丘顶部、陆上资料中表层或浅层火成岩之间，均会有层间多次波与有效波相干涉。与地表面相关多次波相比，层间多次波反射路径更为复杂，规律性更弱，因而在剖面上非常不易识别。

为了消除层间多次波对地震资料的影响，当前工业界主要采用基于周期性的预测反褶积和基于多次波与一次波速度差异拉东域多次波衰减方法，但是当多次波周期性不明显或多次波与一次波速度差异小的情况，这类方法处理效果非常有限。目前现有其它较为先进的层间多次波去除方法有王德利等(2011)提出的CFP技术层间多次波预测及小曲线Curvelet域相减方法，该方法的根源要追溯到Berkhout教授最早提出的共聚焦点(CFP)成像技术。该技术具有一定的抗噪能力，但是由于是基于波场延拓的方法，需要具有较高的速度精度。另一类方法是Weglein等人提出完全数据驱动型的逆散射去层间多次波方法，该方法可从逆散射子序列中一次预测出所有的层间多次波，但运行成本非常高，在工业界至今还没有应用。

层间多次波的衰减方法是当前世界性难题，随着油气勘探开发的不断深入，该项技术的研究有着重要的应用价值，对提高地震资料处理品质有着重要意义。

发明内容

为克服现有技术中存在的问题，本发明提供一种海洋地震资料处理中层间多次波衰减的方法及装置，以实现纯数据驱动进行层间多次波预测和衰减，从而不依赖与速度场的精度，提高海洋地震资料的处理质量。

为达到上述目的，本发明一方面提供了一种海洋地震资料处理中层间多次波衰减的方法，包括：

对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)；

对所述预处理后的地震记录P_PR(x,z₀;t)进行动校正，得到动校正后的地震记录；

将所述动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置；

根据所述界面L的位置，在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)；

根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校，得到第一波场项P₁(x,z₀;t)；

将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)；

将所述第一波场项P₁(x,z₀;t)、第二波场项P₂(x,z₀;t)和第三波场项P₃(x,z₀;t)进行频域变换，得到频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)；

根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)；

将所述层间多次波P_LM(x,z₀;w)进行时域变换，得到时间域的层间多次波P_LM(x,z₀;t)；

在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)。

本发明另一方面提供了一种海洋地震资料处理中层间多次波衰减的装置，包括：

预处理单元，用于对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)；

动校正单元，用于对所述预处理后的地震记录P_PR(x,z₀;t)进行动校正，得到动校正后的地震记录；

界面L拾取单元，用于将所述动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置；

第二波场项计算单元，用于根据所述界面L的位置，在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)；

第一波场项计算单元，用于根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校，得到第一波场项P₁(x,z₀;t)；

第三波场项计算单元，用于将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)；

频域变换单元，用于将所述第一波场项P₁(x,z₀;t)、第二波场项P₂(x,z₀;t)和第三波场项P₃(x,z₀;t)进行频域变换，得到频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)；

层间多次波预测单元，用于根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)；

时域变换单元，用于将所述层间多次波P_LM(x,z₀;w)进行时域变换，得到时间域的层间多次波P_LM(x,z₀;t)；

多次波衰减单元，用于在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)。

本发明能够达到的有益效果：通过纯数据驱动型的层间多次波预测和衰减，能够不依赖与速度场的精度，提高了海洋地震资料的处理质量。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1是本发明的原理示意图；

图2是具体实施方式的简单模型速度场；

图3中，a是包含层间多次波的预处理后的地震记录；b是产生层间多次波高速层的一次反射波，也就是第二波场项；c是产生层间多次波高速层以下层界面的反射波场，是第一波场项，也近似为第三波场项；

图4中，a是包含层间多次波的原始单炮记录；b是匹配衰减完层间多次波后的记录；c是预测出的层间多次波；

图5是本发明实施例一种海洋地震资料处理中层间多次波衰减的方法的流程图；

图6是本发明实施例一种海洋地震资料处理中层间多次波衰减的装置的结构图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。

图1是本发明的简单原理图，地表有实际激发震源S，接收记录R，地下存在强反射界面L。层间多次波传播路径是SO₁OO₂R，S₁与R₁分别是虚震源和虚接收点。将波场分为三个波场项来构建该层间多次波传播路径，第一波场项传播路径是SO₁R₁，第二波场项S₁OR₁和第三波场项S₁O₂R，这三个波场项都是可以在地表接收到的，层间多次波SO₁OO₂R可以通过第一波场项SO₁R₁与第三波场项S₁O₂R的褶积再与第二波场项S₁OR₁的相关来得到。

基于强反射界面间会产生可识别的层间多次波这一认识，建立如图2所示的理论模型，该模型第二层为一高速层，速度为4000m/s，层间多次波主要由该层顶界面与其下反射层产生，采用二阶声波有限差分模拟合成记录，子波主频20Hz，每炮101道接收，道距20米，中间放炮，炮距20米，共模拟单炮101炮，如图3中，(a)为2000米处激发的单炮合成记录，其中,M1，M2和M3均为与该强反射层相关的层间多次波。为突出层间多次波，模拟时在表面加吸收边界，该记录中不存在自由表面多次波。按照本发明的方法原理，利用(1)式预测高速层顶界面产生的层间多次波，需要从去掉干扰波后的数据中提取产生层间多次波界面的一次反射波也就是第二波场项，以及该界面以下的所有反射波即第一波场项。为防止直接切除损伤部分有效波，先将炮集记录抽成CMP道集，动校正后选择矩形窗切除，然后反动校抽成炮集，即可获得所需要的两类波场，图3中，(b)是层间多次波相关界面的一次反射波P₂(x,z₀;t)，将该界面以上所有反射切除后的波场为第二波场项P₁(x,z₀;t)，即图3中的(c)，同时也作为第三波场项P₃(x,z₀;t)。

图4中的(a)是根据式(2)相减后得到的衰减层间多次波后的记录P_ME(x,z₀;t)，图4中的(b)是根据式(1)预测出的与4000m/s速度层相关的层间多次波，从图中可以看出，层间多次波M1，M2和M3均得到了有效压制。

具体通过以下实施例进行说明：

实施例一

如图5所示，为本实施例一种海洋地震资料处理中层间多次波衰减的方法的流程图，包括以下步骤：

步骤501，对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)。

对原始地震记录进行预处理，包括面波去除、随机噪音衰减、野值去除和地表相关多次波衰减，得到预处理后的地震记录P_PR(x,z₀;t)。

步骤502，对所述预处理后的地震记录P_PR(x,z₀;t)进行动校正，得到动校正后的地震记录；

将整个炮集动校正的目的是为了更容易的拾取处界面L的一次反射，拾取出来后再反动校回去就相当于从炮集数据中提取出了L层的一次反射。

步骤503，将所述动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置。

步骤504，根据所述界面L的位置，在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)。

步骤505，根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校，得到第一波场项P₁(x,z₀;t)。

步骤506，将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)。

步骤507，将所述第一波场项P₁(x,z₀;t)、第二波场项P₂(x,z₀;t)和第三波场项P₃(x,z₀;t)进行频域变换，得到频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)。

步骤508，根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)。

基于地震干涉测量理论，在波场运动学中，褶积过程可以认为是两个波场传播路径的求和，而相关是两个传播路径的相减。在时间域的褶积等价于在频率域的乘积，在时间域的相关等价于在频率域一个信号与另外一个信号的复共轭的乘积。在频率域，定义界面L以下层的所有反射波场为第一波场项P₁(x,z₀;w)，定义界面L的一次反射波场为第二波场项P₂(x,z₀;w)，实际生产中无法确定界面L以下所有构造的一次反射P₃(x,z₀;w)，所以近似第三波场项P₃(x,z₀;w)等于第一波场项P₁(x,z₀;w)。在频率域，第一波场项与第二波场项复共轭的乘积再与第三波场项的乘积就构建了界面L相关的层间多次波场P_LM(x,z₀;w)。公式表示为：

$P_{\mathrm{LM}}(x,z_0;w)=P_1(x,z_0;w)\·P_2^{*}(x,z_0;w)\·P_3(x,z_0;w);$

其中，x是空间横向位置坐标，z₀代表波场在地表接收，w代表频率域，P_LM(x,z₀;w)是预测出的与界面L相关的层间多次波。

步骤509，将所述层间多次波P_LM(x,z₀;w)进行时域变换，得到时间域的层间多次波P_LM(x,z₀;t)。

步骤510，在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)。

预测的层间多次波P_LM(x,z₀;t)与炮集数据P_PR(x,z₀;t)需要进行能量匹配，在相同的数量级后才能进行匹配相减，基于最小二乘匹配滤波方法，得到匹配因子，然后从炮集数据P_pr(x,z₀;t)中减去匹配滤波后的层间多次波，即得到衰减掉的与L层相关层间多次波记录，数学描述如下：

$P_{\mathrm{ME}}(x,z_0;t)=P_{\mathrm{PR}}(x,z_0;t)-\underset{t=0}{\overset{\mathrm{NT}}{\mathrm{\Σ}}}A(x,\mathrm{\τ})\·P_{\mathrm{LM}}(x,z_0;t-\mathrm{\τ})$

其中，t是记录时间，A(x,τ)是匹配因子，NT是整个记录时间的长度。

实施例二

如图6所示，为本实施例一种海洋地震资料处理中层间多次波衰减的装置的结构图，包括以下单元：

预处理单元601，用于对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)；

对原始地震记录进行预处理，包括面波去除、随机噪音衰减、野值去除和地表相关多次波衰减，得到预处理后的地震记录P_PR(x,z₀;t)。

动校正单元602，用于对所述预处理后的地震记录P_PR(x,z₀;t)进行动校正，得到动校正后的地震记录；

将整个炮集动校正的目的是为了更容易的拾取处界面L的一次反射，拾取出来后再反动校回去就相当于从炮集数据中提取出了L层的一次反射。

界面L拾取单元603，用于将所述动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置；

第二波场项计算单元604，用于根据所述界面L的位置，在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)；

第一波场项计算单元605，用于根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校，得到第一波场项P₁(x,z₀;t)；

第三波场项计算单元606，用于将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)；

频域变换单元607，用于将所述第一波场项P₁(x,z₀;t)、第二波场项P₂(x,z₀;t)和第三波场项P₃(x,z₀;t)进行频域变换，得到频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)；

层间多次波预测单元608，用于根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)；

基于地震干涉测量理论，在波场运动学中，褶积过程可以认为是两个波场传播路径的求和，而相关是两个传播路径的相减。在时间域的褶积等价于在频率域的乘积，在时间域的相关等价于在频率域一个信号与另外一个信号的复共轭的乘积。在频率域，定义界面L以下层的所有反射波场为第一波场项P₁(x,z₀;w)，定义界面L的一次反射波场为第二波场项P₂(x,z₀;w)，实际生产中无法确定界面L以下所有构造的一次反射P₃(x,z₀;w)，所以近似第三波场项P₃(x,z₀;w)等于第一波场项P₁(x,z₀;w)。在频率域，第一波场项与第二波场项复共轭的乘积再与第三波场项的乘积就构建了界面L相关的层间多次波场P_LM(x,z₀;w)。公式表示为：

$P_{\mathrm{LM}}(x,z_0;w)=P_1(x,z_0;w)\·P_2^{*}(x,z_0;w)\·P_3(x,z_0;w)$

其中，x是空间横向位置坐标，z₀代表波场在地表接收，w代表频率域，P_LM(x,z₀;w)是预测出的与界面L相关的层间多次波。

时域变换单元609，用于将所述层间多次波P_LM(x,z₀;w)进行时域变换，得到时间域的层间多次波P_LM(x,z₀;t)；

多次波衰减单元610，用于在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)。

预测的层间多次波P_LM(x,z₀;t)与炮集数据P_PR(x,z₀;t)需要进行能量匹配，在相同的数量级后才能进行匹配相减，基于最小二乘匹配滤波方法，得到匹配因子，然后从炮集数据P_pr(x,z₀;t)中减去匹配滤波后的层间多次波，即得到衰减掉的与L层相关层间多次波记录，数学描述如下：

$P_{\mathrm{ME}}(x,z_0;t)=P_{\mathrm{PR}}(x,z_0;t)-\underset{t=0}{\overset{\mathrm{NT}}{\mathrm{\Σ}}}A(x,\mathrm{\τ})\·P_{\mathrm{LM}}(x,z_0;t-\mathrm{\τ})$

其中，t是记录时间，A(x,τ)是匹配因子，NT是整个记录时间的长度。

本发明的有益效果是：该方法是纯数据驱动型，不依赖与速度场的精度，是一种具有较高应用价值的层间多次波衰减方法。

上述实施例能够达到的有益效果：通过纯数据驱动型的层间多次波预测和衰减，能够不依赖与速度场的精度，提高了海洋地震资料的处理质量。

本领域一般技术人员在此设计思想之下所做任何不具有创造性的改造，均应视为在本发明的保护范围之内。

Claims (8)

1.一种海洋地震资料处理中层间多次波衰减的方法，其特征在于，包括：

对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)；

对所述预处理后的地震记录P_PR(x,z₀;t)进行动校正，得到动校正后的地震记录；

将所述动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置；

根据所述界面L的位置，在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)；

根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校，得到第一波场项P₁(x,z₀;t)；

将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)；

将所述第一波场项P₁(x,z₀;t)、第二波场项P₂(x,z₀;t)和第三波场项P₃(x,z₀;t)进行频域变换，得到频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)；

根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)；

将所述层间多次波P_LM(x,z₀;w)进行时域变换，得到时间域的层间多次波P_LM(x,z₀;t)；

在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)。

2.如权利要求1所述海洋地震资料处理中层间多次波衰减的方法，其特征在于，所述根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)具体为：

$P_{\mathrm{LM}}(x,z_0;w)=P_1(x,z_0;w)\·P_2^{*}(x,z_0;w)\·P_3(x,z_0;w);$

其中，x是空间横向位置坐标，z₀代表波场在地表接收，w代表频率域。

3.如权利要求1所述海洋地震资料处理中层间多次波衰减的方法，其特征在于，所述将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)的算法为：

P₃(x,z₀;t)=P₁(x,z₀;t)。

4.如权利要求1所述海洋地震资料处理中层间多次波衰减的方法，其特征在于，所述在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)具体为：

根据最小二乘法在所述预处理后的地震记录P_PR(x,z₀;t)中获取能量匹配因子A(x,τ)；

根据所述匹配因子对所述层间多次波P_LM(x,z₀;t)进行能量匹配，并在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉进行能量匹配后的层间多次波，得到衰减掉层间多次波的记录P_ME(x,z₀;t)，具体为：

$P_{\mathrm{ME}}(x,z_0;t)=P_{\mathrm{PR}}(x,z_0;t)-\underset{t=0}{\overset{\mathrm{NT}}{\mathrm{\Σ}}}A(x,\mathrm{\τ})\·P_{\mathrm{LM}}(x,z_0;t-\mathrm{\τ});$

其中，t是记录时间，NT是整个记录时间的长度。

5.一种海洋地震资料处理中层间多次波衰减的装置，其特征在于，包括：

预处理单元，用于对原始地震记录进行预处理，得到预处理后的地震记录P_PR(x,z₀;t)；

动校正单元，用于对所述预处理后的地震记录P_PR(x,z₀;t)进行动校正，得到动校正后的地震记录；

界面L拾取单元，用于将所述动校正后的地震记录进行叠加，在叠加后的剖面上获取产生层间多次波的强反射界面L的位置；

第二波场项计算单元，用于根据所述界面L的位置，在所述动校正后的地震记录中拾取界面L产生的一次反射波，并将该一次反射波进行反动校，得到第二波场项P₂(x,z₀;t)；

第一波场项计算单元，用于根据所述界面L的位置，删除所述动校正后的地震记录中该界面L以上的反射波场，并将删除所述界面L以上的反射波场后的所述动校正后的地震记录进行反动校，得到第一波场项P₁(x,z₀;t)；

第三波场项计算单元，用于将所述删除所述界面L以上的反射波场后的所述动校正后的地震记录中所有界面的一次反射波场进行反动校，得到第三波场项P₃(x,z₀;t)；

频域变换单元，用于将所述第一波场项P₁(x,z₀;t)、第二波场项P₂(x,z₀;t)和第三波场项P₃(x,z₀;t)进行频域变换，得到频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)；

层间多次波预测单元，用于根据所述频率域的第一波场项P₁(x,z₀;w)、第二波场项P₂(x,z₀;w)和第三波场项P₃(x,z₀;w)计算得到与所述界面L相关的层间多次波P_LM(x,z₀;w)；

时域变换单元，用于将所述层间多次波P_LM(x,z₀;w)进行时域变换，得到时间域的层间多次波P_LM(x,z₀;t)；

多次波衰减单元，用于在所述预处理后的地震记录P_PR(x,z₀;t)中衰减掉所述时间域的层间多次波P_LM(x,z₀;t)，得到衰减掉层间多次波的记录P_ME(x,z₀;t)。

6.如权利要求5所述海洋地震资料处理中层间多次波衰减的装置，其特征在于，所述层间多次波预测单元根据以下公式计算层间多次波P_LM(x,z₀;w)：

$P_{\mathrm{LM}}(x,z_0;w)=P_1(x,z_0;w)\·P_2^{*}(x,z_0;w)\·P_3(x,z_0;w);$

其中，x是空间横向位置坐标，z₀代表波场在地表接收，w代表频率域。

7.如权利要求5所述海洋地震资料处理中层间多次波衰减的装置，其特征在于，所述第三波场项计算单元计算第三波场项P₃(x,z₀;t)的算法为：

P₃(x,z₀;t)=P₁(x,z₀;t)。

8.如权利要求5所述海洋地震资料处理中层间多次波衰减的装置，其特征在于，所述多次波衰减单元根据以下公式进行衰减：

$P_{\mathrm{ME}}(x,z_0;t)=P_{\mathrm{PR}}(x,z_0;t)-\underset{t=0}{\overset{\mathrm{NT}}{\mathrm{\Σ}}}A(x,\mathrm{\τ})\·P_{\mathrm{LM}}(x,z_0;t-\mathrm{\τ});$

其中，t是记录时间，NT是整个记录时间的长度，A(x,τ)是根据最小二乘法在所述预处理后的地震记录P_PR(x,z₀;t)中获取的能量匹配因子。

Images