CN105425300B - 一种剩余静校正方法 - Google Patents

Abstract

本申请实施例公开了一种剩余静校正方法。所述方法包括：获取地震数据，并对所述地震数据进行动校正处理；对动校正处理后的地震数据进行叠加，生成叠加剖面，通过所述叠加剖面获取反射层位；获取地震道的时间窗，将所述时间窗内反射层位的地震子波进行叠加，并将叠加后形成的地震道作为超级道；获取每个炮点所对应地震道的超级道，并基于该炮点的超级道计算该炮点的剩余静校正量；获取每个检波点所对应地震道的超级道，并基于该检波点的超级道计算该检波点的剩余静校正量。本申请实施例的剩余静校正方法可以减少剩余静校正计算所需的时间，同时可以增强剩余静校正的抗噪声性能。

Description

一种剩余静校正方法

技术领域

本申请涉及地球物理勘探技术领域，特别涉及一种剩余静校正方法。

背景技术

剩余静校正是对地震资料所作的校正，用于补偿由高程、风化层厚度以及风化层速度产生的影响，把资料校到一个指定的基准面上。它是一种地震勘探中消除近地表影响从而提高地震反射波成像质量的主要数据处理方法。

反射波法是一种常用的剩余静校正方法。反射波法一般是利用地震数据的反射波信息，通过互相关法或者能量最大法等准则，利用统计的方法来计算剩余静校正量。目前的反射波剩余静校正方法主要包括地表一致性的剩余静校正方法和自动剩余静校正方法。例如，2005年9月李全胜等人在“吐哈油气”第10卷第3期中公开了一种基于模型道的地表一致性剩余静校正方法，2003年2月井西利在“石油地球物理勘探”第38卷第1期中公开了一种大剩余静校正量求解的两步法。

在实现本申请过程中，发明人发现现有技术中至少存在如下问题：

地表一致性剩余静校正方法受噪声影响较大，抗噪性较差，尤其是在低信噪比地区的校正结果不尽如人意。而自动剩余静校正方法的计算过程比较复杂，计算耗时较长，不适用于大数据量的计算。尤其是在覆盖次数较高、测线较长的情况下，自动剩余静校正方法一般需要三至四天的时间才能够得到校正结果。

发明内容

本申请实施例的目的是提供一种剩余静校正方法。该方法可以减少剩余静校正计算所需的时间，同时可以增强剩余静校正的抗噪声性能。

为解决上述技术问题，本申请实施例提供的一种剩余静校正方法是这样实现的：

一种剩余静校正方法，包括：

获取地震数据，并对所述地震数据进行动校正处理；

对动校正处理后的地震数据进行叠加，生成叠加剖面，通过所述叠加剖面获取反射层位；

获取时间窗以及地震道，将所述时间窗内反射层位的地震子波进行叠加，并将叠加后形成的地震道作为超级道；

获取每个炮点所对应地震道的超级道，并基于该炮点的超级道计算该炮点的剩余静校正量；

获取每个检波点所对应地震道的超级道，并基于该检波点的超级道计算该检波点的剩余静校正量。

由以上本申请实施例提供的技术方案可见，本申请实施例通过超级道数据计算炮点和检波点的剩余静校正量，在减少剩余静校正计算所需时间的同时，可以提升地震资料的信噪比，从而增强剩余静校正的抗噪声性能，使剩余静校正的计算结果更加准确。

附图说明

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

图1为本申请实施例剩余静校正方法的流程图；

图2为工区剩余静校正前的叠加剖面；

图3为工区应用常规剩余静校正方法校正后的叠加剖面；

图4为工区应用本申请实施例的剩余静校正方法校正后的叠加剖面。

具体实施方式

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

本申请实施例提供的剩余静校正方法，如图1所示，具体包括以下步骤：

S101：获取地震数据，并对所述地震数据进行动校正处理。

具体地，可以获取工区内的地震数据，并对所述地震数据进行动校正处理。

在一些实施方式中，在获取到工区内的地震数据后，可以对获取的地震数据进行能量均衡、地表一致性震幅补偿、地表一致性反褶积、滤波、震幅谱增益分析和速度分析等预处理，并对预处理后的地震数据进行动校正处理。

S102：对动校正处理后的地震数据进行叠加，生成叠加剖面，通过所述叠加剖面获取反射层位。

所述叠加剖面可以是地震反射剖面。具体地，可以对地震反射剖面进行时间-地层单元划分，从而得到反射层位。一般地，反射波经过地下地层的反射可以形成震幅层。由于震幅层与地下地层具有对应关系，因此可以根据反射波的旅行时间和该反射波所形成的震幅层的震幅，从所述叠加剖面上获取反射层位。例如，在所述叠加剖面上，可以将旅行时间和震幅大致相同的层位作为一个反射层位。

在一些实施方式中，可以从叠加剖面上手工拾取旅行时间和震幅大致相同的层位，从而得到反射层位。

在另一些实施方式中，可以采用扫描计算的方法从叠加剖面上拾取旅行时间和震幅大致相同的层位，从而得到反射层位。

S103：获取地震道的时间窗，将所述时间窗内反射层位的地震子波进行叠加，并将叠加后形成的地震道作为超级道。

具体地，可以获取待计算的目标地层，然后获取该目标地层所对应的反射波旅行时间范围，并划定时间窗。一般地，时间窗内地震子波是不变的，因此，可以把地震子波叠加起来，从而形成超级道。

超级道可以大幅度提高地震数据的抗噪性和计算的高效性。需要说明的是，超级道一般不同于超长道。超长道通常在模拟退火方法中使用。具体地，将时间窗内的地震数据经过扫描，得到多个强同相轴子波，然后将所述强同相轴子波叠加可以得到超级道。将炮集内或者检波点内道集的所有地震记录串联可以得到超长道。地震记录串联的过程通常是将地震记录首尾相连接的过程。

S104：获取每个炮点所对应地震道的超级道，并基于该炮点的超级道计算该炮点的剩余静校正量。

具体地，可以获取工区内的炮点。对于每个炮点，可以进行如下的处理，以计算该炮点的剩余静校正量，进而得到工区内每个炮点的剩余静校正量。所述处理具体包括：

1)获取该炮点所对应的超级道，将该炮点所对应的超级道数据进行叠加，然后将叠加后的超级道数据作为该炮点的超级炮数据。

2)获取该炮点所对应的地震道集，将该炮点所对应的各个地震道数据进行叠加，然后将叠加后的地震道数据作为该炮点的超级模型道数据。

3)将该炮点的超级炮数据和超级模型道数据进行互相关，从而得到该炮点的剩余静校正量。

其中，步骤3)中所述的互相关也称为“互协方差”，一般用来表示两个信号之间相似性的一个度量。

在步骤3)中，通过将该炮点的超级炮数据和超级模型道数据进行互相关，可以直接得到该炮点的剩余静校正量。而常规的剩余静校正方法一般是将地震子波叠加形成模型道，然后将炮点的CMP道集(Common Middle Point，共中心点道集)和该CMP道集对应的模型道进行互相关，得到每道地震数据的时差，最后对每道地震数据集的时差进行分解，得到该炮点的剩余静校正量。因此，与常规的剩余静校正方法相比，在步骤3)中，通过将炮点的超级炮数据和超级模型道数据进行互相关，可以无需对每道地震数据集的时差进行分解，从而可以提高剩余静校正量的计算效率。

S105：获取每个检波点所对应地震道的超级道，并基于该检波点的超级道计算该检波点的剩余静校正量。

具体地，可以获取工区内的检波点。对于每个检波点，可以进行如下的处理，以计算该检波点的剩余静校正量，进而得到工区内每个检波点的剩余静校正量。所述处理具体包括：

4)获取该检波点所对应的超级道，将该检波点所对应的超级道数据进行叠加，然后将叠加后的超级道数据作为该检波点的超级检波点数据。

5)获取该检波点所对应的地震道集，将该检波点所对应的各个地震道数据进行叠加，然后将叠加后的地震道数据作为该检波点的超级模型道数据。

6)将该检波点的超级检波点数据和超级模型道数据进行互相关，从而得到该检波点的剩余静校正量。

本申请实施例的剩余静校正方法，通过超级道数据计算炮点和检波点的剩余静校正量，在减少剩余静校正计算所需时间的同时，可以提升地震资料的信噪比，从而增强剩余静校正的抗噪声性能，使剩余静校正的计算结果更加准确。

进一步地，本申请实施例的剩余静校正方法，综合了地表一致性剩余静校正方法和自动剩余静校正方法的优点，从而使得该剩余静校正方法具有自动剩余静校正的高抗噪性和地表一致性剩余静校正的高效性。

图2为工区剩余静校正前的叠加剖面，图3为该工区内应用常规剩余静校正方法校正后的叠加剖面，图4为该工区内应用本申请实施例的剩余静校正方法校正后的叠加剖面。通过图2、图3和图4的对比可以发现，本申请实施例的剩余静校正方法具有更好的校正效果。

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

Claims (6)

1.一种剩余静校正方法，其特征在于，包括：

获取地震数据，并对所述地震数据进行动校正处理；

对动校正处理后的地震数据进行叠加，生成叠加剖面，通过所述叠加剖面获取反射层位；

获取地震道的时间窗，将所述时间窗内反射层位的地震子波进行叠加，并将叠加后形成的地震道作为超级道；

获取每个炮点所对应地震道的超级道，并基于该炮点的超级道计算该炮点的剩余静校正量；

获取每个检波点所对应地震道的超级道，并基于该检波点的超级道计算该检波点的剩余静校正量。

2.如权利要求1所述的方法，其特征在于，所述获取每个炮点所对应地震道的超级道，并基于该炮点的超级道计算该炮点的剩余静校正量，具体包括：

获取每个炮点所对应的超级道，将该炮点所对应的超级道数据进行叠加，并将叠加后的超级道数据作为该炮点的超级炮数据；以及，

获取该炮点所对应的地震道集，将该炮点所对应的各个地震道数据进行叠加，并将叠加后的地震道数据作为该炮点的超级模型道数据；以及，

将该炮点的超级炮数据和超级模型道数据进行互相关，得到该炮点的剩余静校正量。

3.如权利要求1所述的方法，其特征在于，所述获取每个检波点所对应地震道的超级道，并基于该检波点的超级道计算该检波点的剩余静校正量，具体包括：

获取每个检波点所对应的超级道，将该检波点所对应的超级道数据进行叠加，并将叠加后的超级道数据作为该检波点的超级检波点数据；以及，

获取该检波点所对应的地震道集，将该检波点所对应的各个地震道数据进行叠加，并将叠加后的地震道数据作为该检波点的超级模型道数据；以及，

将该检波点的超级炮数据和超级模型道数据进行互相关，得到该检波点的剩余静校正量。

4.如权利要求1所述的方法，其特征在于，所述获取地震道的时间窗，具体包括：

获取待计算的目标地层，以及该目标地层所对应的反射波的旅行时间范围，根据该旅行时间范围在所述叠加剖面上划定时间窗。

5.如权利要求1所述的方法，其特征在于，所述通过所述叠加剖面获取反射层位，具体包括：

根据反射波的旅行时间和该反射波所形成的震幅层的震幅，从所述叠加剖面上获取反射层位。

6.如权利要求1所述的方法，其特征在于，在获取地震数据后，所述方法还包括：

对所述地震数据进行预处理，所述预处理包括能量均衡处理、地表一致性震幅补偿处理、地表一致性反褶积处理、滤波处理、震幅谱增益分析和速度分析；

相应地，所述对所述地震数据进行动校正处理，具体包括：

对预处理后的地震数据进行动校正处理。