CN110824561B - 度量时移地震数据可重复性的方法 - Google Patents

Abstract

本申请公开了一种度量时移地震数据可重复性的方法，属于油田勘探技术领域。该方法包括：先将开始数据和监控数据作差得到差数据，再利用凯泽窗希尔伯特算子分别得到开始数据、监控数据和差数据对应的凯泽窗希尔伯特道数据，进而得到地震数据的包络数据，利用包络数据计算包络和，并得到包络特征可重复性因子，最终确定截断门槛值，利用该截断门槛值将包络特征可重复性因子转化为包络特征可重复性指示因子，这样，能够根据取值为0或1的包络特征可重复性指示因子，直观地度量地震数据是否具有可重复性。

Description

度量时移地震数据可重复性的方法

技术领域

本申请涉及油田勘探技术领域，特别涉及一种四维地震数据可重复性度量的方法。

背景技术

时移地震技术是一种地震勘探技术，即在不同时间对同一区域激发人工地震，并利用监测仪器，按相同的观测方式进行重复性观测以获取地震数据，最后利用获取到的地震数据研究油气藏的动态变化。

从理论上讲，地震数据包括与油气藏的动态性质相关的数据和与油气藏的静态性质相关的数据，动态性质指可能随时间发生变化的性质，在不同时刻下的值通常是不一样的，静态性质指基本不随时间发生变化的性质，在不同时刻下的值通常是一样的。例如：动态性质可以包括流体的饱和度、压力和温度，静态性质可以包括地层构造和岩性。对在不同时间获取到的地震数据进行一定处理，使与油气藏的静态性质相关的数据具有可重复性，保留与油气藏的动态性质相关的数据的差异性，这样，将之后得到的地震数据与最开始得到的地震数据作差，得到的差数据能反应油气藏随时间的变化情况。

但在实际中，时移地震数据的采集是间隔性的，很难保证在不同时间下的采集条件能完全一致，而采集条件的变化会对地震数据产生很多干扰，导致在不同时间获取到的地震数据的可重复性差，这样就难以准确地得到油气藏随时间的变化情况，因此需要采取一定的数据质量控制措施，来改善地震数据的可重复性。

相关技术中，通常采用计算两期地震数据的相干性估计值的方法，来度量地震数据的可重复性。相干性估计值越大，说明两期地震数据的相似程度越高，地震数据的可重复性越好。

在实现本申请的过程中，发明人发现相关技术中至少存在以下问题：

上述方法只能定性度量地震数据的可重复性，对地震数据可重复性好和差的区分界限比较模糊，难以直观地度量地震数据是否具有可重复性。

发明内容

有鉴于此，本申请实施例提供了一种度量时移地震数据可重复性的方法，能够直观地度量时移地震数据的可重复性。所述技术方案如下：

本申请实施例提供了一种度量时移地震数据可重复性的方法，所述方法包括：

步骤1，根据开始数据B_i，j和监控数据M_i，j，计算差数据D_i，j，i是地震道的顺序号，i＝1，2，…，I，I是地震道的总数，j是时间样点的顺序号，j＝1，2，…，J，J是时间样点的总数，所述开始数据B_i，j为第一次观测时获取到的地震数据，所述监控数据M_i，j为后续观测时获取到的地震数据；

步骤2，利用凯泽窗希尔伯特算子，对所述开始数据B_i，j、所述监控数据M_i，j和所述差数据D_i，j进行变换，分别得到所述开始数据的凯泽窗希尔伯特道数据U_i，j、所述监控数据的凯泽窗希尔伯特道数据V_i，j和所述差数据的凯泽窗希尔伯特道数据W_i，j；

步骤3，分别根据所述开始数据B_i，j和所述开始数据的凯泽窗希尔伯特道数据U_i，j、所述监控数据M_i，j和所述监控数据的凯泽窗希尔伯特道数据V_i，j以及所述差数据D_i，j和所述差数据的凯泽窗希尔伯特道数据W_i，j，得到所述开始数据的包络数据EB_i，j、所述监控数据的包络数据EM_i，j和所述差数据的包络数据ED_i，j；

步骤4，分别根据所述开始数据的包络数据EB_i，j、所述监控数据的包络数据EM_i，j和所述差数据的包络数据ED_i，j，得到所述开始数据的包络和SEB_i、所述监控数据的包络和SEM_i以及所述差数据的包络和SED_i；

步骤5，根据所述开始数据的包络和SEB_i、所述监控数据的包络和SEM_i以及所述差数据的包络和SED_i，得到包络特征可重复性因子NES_i；

步骤6，根据所述包络特征可重复性因子NES_i，确定截断门槛值NSECut；

步骤7，如果第i地震道的所述包络特征可重复性因子NES_i小于或等于所述截断门槛值NSECut，则判定所述第i地震道的地震数据具有可重复性。

可选的，所述步骤1包括：

根据开始数据B_i，j和监控数据M_i，j，按照公式D_i,j＝B_i,j-M_i,j计算差数据D_i，j。

可选的，所述步骤2包括：

步骤21，按照以下公式计算所述凯泽窗希尔伯特算子H_k，

其中，K指示所述凯泽窗希尔伯特算子的样点数，(2K+1)是所述凯泽窗希尔伯特算子的样点总数，

和X₀(α)是第一类零阶修正贝塞尔函数，α是凯泽窗函数的形状参数；

步骤22，根据所述凯泽窗希尔伯特算子H_k，按照以下公式，计算所述开始数据的凯泽窗希尔伯特道数据U_i，j；

按照以下公式，计算所述监控数据的凯泽窗希尔伯特道数据V_i，j；

按照以下公式，计算所述差数据的凯泽窗希尔伯特道数据W_i，j；

其中，当j-k≤0时，B_i，j-k、M_i，j-k和D_i，j-k的取值均为0。

可选的，所述步骤3包括：

根据所述开始数据B_i，j和所述开始数据的凯泽窗希尔伯特道数据U_i，j，按照以下公式，得到所述开始数据的复数道E_i，j；

E_i,j＝B_i,j+LU_i,j

根据所述监控数据M_i，j和所述监控数据的凯泽窗希尔伯特道数据V_i，j，按照以下公式，得到所述监控数据的复数道F_i，j；

F_i,j＝M_i,j+LV_i,j

根据所述差数据D_i，j和所述差数据的凯泽窗希尔伯特道数据W_i，j，按照以下公式，得到所述差数据的复数道G_i，j；

G_i,j＝D_i,j+LW_i,j

其中，L是虚数单位，L²＝-1。

可选的，所述步骤3还包括：

按照公式

对所述开始数据的复数道E_i，j取模，得到所述开始数据的包络数据EB_i，j；

按照公式

对所述开始数据的复数道F_i，j取模，得到所述开始数据的包络数据EM_i，j；

按照公式

对所述差数据的复数道G_i，j取模，得到所述开始数据的包络数据ED_i，j。

可选的，所述步骤4包括：

根据所述开始数据的包络数据EB_i，j，按照公式

计算所述开始数据的包络和SEB_i；

根据所述监控数据的包络数据EM_i，j，按照公式

计算所述监控数据的包络和SEM_i；

根据所述差数据的包络数据ED_i，j，按照公式

计算所述差数据的包络和SED_i。

可选的，在所述步骤4之后，所述方法还包括：

根据所述差数据的包络和SED_i以及所述开始数据的包络和SEB_i，按照公式

计算所述差数据与所述开始数据的包络特征差异化因子RDB_i；

根据所述差数据的包络和SED_i以及所述监控数据的包络和SEM_i，按照公式

计算所述差数据与所述监控数据的包络特征差异化因子RDM_i；

根据所述监控数据的包络和SEM_i以及所述开始数据的包络和SEB_i，按照公式

计算所述监控数据与所述开始数据的包络特征差异化因子RMB_i。

可选的，所述步骤5包括：

根据所述开始数据的包络和SEB_i、所述监控数据的包络和SEM_i以及所述差数据的包络和SED_i，按照公式

得到包络特征可重复性因子NES_i。

可选的，所述步骤6包括：

根据第一地震道的包络特征可重复性因子NES₁，按照公式NSECut＝nNES₁ 0<n<1，确定截断门槛值NSECut。

可选的，所述步骤7包括：

根据所述截断门槛值NSECut与所述包络特征可重复性因子NES_i的值，按照以下公式，确定包络特征可重复性指示因子INSE_i；

如果INSE_i＝1，则判定所述第i地震道的地震数据具有可重复性，如果INSE_i＝0，则判定所述第i地震道的地震数据不具有可重复性。

本申请实施例提供的技术方案至少可以包括以下有益效果：

在本申请实施例中，使用上述的度量时移地震数据可重复性的方法时，先将开始数据和监控数据作差得到差数据，再利用凯泽窗希尔伯特算子分别得到开始数据、监控数据和差数据对应的凯泽窗希尔伯特道数据，进而得到地震数据的包络数据，利用包络数据计算包络和，并得到包络特征可重复性因子，最终确定截断门槛值，利用该截断门槛值将包络特征可重复性因子转化为包络特征可重复性指示因子，这样，能够根据取值为0或1的包络特征可重复性指示因子，直观地度量地震数据是否具有可重复性。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。

附图说明

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

图1是本申请实施例提供的一种度量时移地震数据可重复性的方法的流程图；

图2是本申请实施例提供的一种开始数据的示意图；

图3是本申请实施例提供的一种监控数据的示意图；

图4是本申请实施例提供的一种开始数据的包络数据的示意图；

图5是本申请实施例提供的一种监控数据的包络数据的示意图；

图6是本申请实施例提供的一种差数据的包络数据的示意图；

图7是本申请实施例提供的一种开始数据的包络和的示意图；

图8是本申请实施例提供的一种监控数据的包络和的示意图；

图9是本申请实施例提供的一种差数据的包络和的示意图；

图10是本申请实施例提供的一种差数据与开始数据的包络特征差异化因子的示意图；

图11是本申请实施例提供的一种差数据与监控数据的包络特征差异化因子的示意图；

图12是本申请实施例提供的一种监控数据与开始数据的包络特征差异化因子的示意图；

图13是本申请实施例提供的一种包络特征可重复性因子的示意图；

图14是本申请实施例提供的一种截断门槛值的示意图；

图15是本申请实施例提供的另一种截断门槛值的示意图；

图16是本申请实施例提供的一种包络特征可重复性指示因子的示意图；

图17是本申请实施例提供的另一种包络特征可重复性指示因子的示意图。

具体实施方式

为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。

在油田勘探技术领域，技术人员常利用时移地震勘探技术对油气藏进行勘探。在不同时间下对目标地区进行观测并获取地震数据，对两期地震数据进行处理，使得其中能够反应油气藏不随时间变化的性质的部分具有可重复性，再将两期地震数据作差，就能够分析油气藏随时间变化的情况。但在实际中，很难保证在不同时间下的采集条件能完全一致，这就有可能导致地震数据的可重复性差，因此需要对地震数据的可重复性进行度量。相关技术中，通常采用计算两期地震数据的相干性估计值的方法，来度量地震数据的可重复性。相干性估计值越大，地震数据的可重复性越好。但上述方法只能定性度量地震数据的可重复性，对地震数据可重复性好和差的区分界限比较模糊，难以直观地度量地震数据是否具有可重复性。且这种方法无法具体分析地震数据中的振幅和相位等因素的情况，具有一定的局限性。本申请提供了一种度量时移地震数据可重复性的方法，能够解决上述问题。

图1是本申请实施例提供的一种度量时移地震数据可重复性的方法的流程图，该方法的执行主体可以是用于处理地震数据的终端，参见图1，该度量时移地震数据可重复性的方法包括以下步骤：

在步骤1中，终端根据开始数据B_i，j和监控数据M_i，j，计算差数据D_i，j。

其中，i是地震道的顺序号，一个地震道代表一个观测点，i＝1，2，…，I，I是地震道的总数。j是时间样点的顺序号，在一次地震数据的观测中，地震会持续一段时间，在这段时间内都可以观测到地震数据，确定数据采样的时间间隔，就会得到若干个时间样点，j＝1，2，…，J，J是时间样点的总数，时间样点数越多，采集精度越高。开始数据B_i，j为第一次观测时获取到的地震数据，可以参考图2，横坐标表示地震道，纵坐标表示时间样点；监控数据M_i，j为后续观测时获取到的地震数据，可以参考图3，横坐标表示地震道，纵坐标表示时间样点。

在实施中，可以利用检波器对地震进行观测，获得地震数据。地震数据包括一系列信号，这些信号中含有振幅、频率和相位等参数的信息，对地震数据进行处理，能够将信号中含有的与特定参数有关的信息提取出来进行分析。

可选的，可以按照以下公式计算差数据。

D_i,j＝B_i,j-M_i,j

在步骤2中，终端利用凯泽窗希尔伯特算子，对开始数据B_i，j、监控数据M_i，j和差数据D_i，j进行变换，分别得到开始数据的凯泽窗希尔伯特道数据U_i，j、监控数据的凯泽窗希尔伯特道数据V_i，j和差数据的凯泽窗希尔伯特道数据W_i，j。

其中，利用凯泽窗希尔伯特算子对开始数据、监控数据和差数据进行变换，能够使地震数据中的信号的相位移相90度，以方便后续对地震数据进行处理。

在实施中，终端利用希尔伯特算子对地震数据进行变换时，需要进行无穷次计算，而这是无法实现的，因此引入凯泽窗函数，将希尔伯特算子限定为凯泽窗希尔伯特算子，这样，数据变换只需进行有限次的计算。

可选的，步骤2具体可以包括：

步骤21，终端按照以下公式计算凯泽窗希尔伯特算子H_k。

其中，K指示凯泽窗希尔伯特算子的样点数，对计算精度要求越高，样点数越大，K的取值也就越大，这样，能够达到更好的计算效果。(2K+1)是凯泽窗希尔伯特算子的样点总数，

和X₀(α)是第一类零阶修正贝塞尔函数，α是凯泽窗函数的形状参数。

步骤22，终端根据凯泽窗希尔伯特算子H_k，按照以下公式，计算开始数据的凯泽窗希尔伯特道数据U_i，j。

按照以下公式，计算监控数据的凯泽窗希尔伯特道数据V_i，j。

按照以下公式，计算差数据的凯泽窗希尔伯特道数据W_i，j。

其中，当j-k≤0时，B_i，j-k、M_i，j-k和D_i，j-k的取值均为0。经计算得到的U_i，j、V_i，j和W_i，j的相位分别比B_i，j，M_i，j和D_i，j的相位偏移了90度。

在步骤3中，终端分别根据开始数据B_i，j和开始数据的凯泽窗希尔伯特道数据U_i，j、监控数据M_i，j和监控数据的凯泽窗希尔伯特道数据V_i，j以及差数据D_i，j和差数据的凯泽窗希尔伯特道数据W_i，j，得到开始数据的包络数据EB_i，j、监控数据的包络数据EM_i，j和差数据的包络数据ED_i，j。

其中，包络数据可以代表地震数据中包含的振幅信息，开始数据的包络数据可以参考图4，监控数据的包络数据可以参考图5，差数据的包络数据可以参考图6，横坐标表示地震道，纵坐标表示时间样点。

可选的，步骤3具体可以包括：

终端根据开始数据B_i，j和开始数据的凯泽窗希尔伯特道数据U_i，j，按照以下公式，得到开始数据的复数道E_i，j。

E_i,j＝B_i,j+LU_i,j

终端根据监控数据M_i，j和监控数据的凯泽窗希尔伯特道数据V_i，j，按照以下公式，得到监控数据的复数道F_i，j。

F_i,j＝M_i,j+LV_i,j

终端根据差数据D_i，j和差数据的凯泽窗希尔伯特道数据W_i，j，按照以下公式，得到差数据的复数道G_i，j。

G_i,j＝D_i,j+LW_i,j

其中，L是虚数单位，L²＝-1。将地震数据的原始数据作为实部，将地震数据通过变换得到的凯泽窗希尔伯特道数据作为虚部，构建地震数据的复数道，方便了后续对包络数据的计算。

可选的，步骤3还可以包括：

按照公式

对开始数据的复数道E_i，j取模，得到开始数据的包络数据EB_i，j；

按照公式

对开始数据的复数道F_i，j取模，得到开始数据的包络数据EM_i，j；

按照公式

对差数据的复数道G_i，j取模，得到开始数据的包络数据ED_i，j。

这样，能够利用复数道的构建，计算出可以代表地震数据中包含的振幅信息的包络数据。

在一些可能的实施例中，终端还可以利用构建出的地震数据的复数道，计算开始数据、监控数据和差数据的瞬时相位。相应的，计算瞬时相位的方法包括：

终端根据开始数据B_i，j和开始数据的凯泽窗希尔伯特道数据U_i，j，按照以下公式，得到开始数据的瞬时相位PB_i，j。

终端根据监控数据M_i，j和监控数据的凯泽窗希尔伯特道数据V_i，j，按照以下公式，得到监控数据的瞬时相位PM_i，j。

终端根据差数据D_i，j和差数据的凯泽窗希尔伯特道数据W_i，j，按照以下公式，得到差数据的瞬时相位PD_i，j。

这样，能够利用地震数据的复数道，对地震数据中包含的相位信息进行提取和分析。

在步骤4中，终端分别根据开始数据的包络数据EB_i，j、监控数据的包络数据EM_i，j和差数据的包络数据ED_i，j，得到开始数据的包络和SEB_i、监控数据的包络和SEM_i以及差数据的包络和SED_i。

其中，开始数据的包络和可以参考图7，监控数据的包络和可以参考图8，差数据的包络和可以参考图9，横坐标表示地震道，纵坐标表示包络和的数值。

可选的，步骤4具体可以包括：

终端根据开始数据的包络数据EB_i，j，按照公式

计算开始数据的包络和SEB_i；

终端根据监控数据的包络数据EM_i，j，按照公式

计算监控数据的包络和SEM_i；

终端根据差数据的包络数据ED_i，j，按照公式

计算差数据的包络和SED_i。

在步骤5中，终端根据开始数据的包络和SEB_i、监控数据的包络和SEM_i以及差数据的包络和SED_i，得到包络特征可重复性因子NES_i。

可选的，步骤5具体可以包括：

终端根据开始数据的包络和SEB_i、监控数据的包络和SEM_i以及差数据的包络和SED_i，按照公式

得到包络特征可重复性因子NES_i。

其中，包络特征可重复性因子可以参考图13，横坐标表示地震道，纵坐标表示包络特征可重复性因子的取值。包络特征可重复性因子NES_i的数值范围为0～200％。例如，如果第i地震道的开始数据和监控数据完全相同，则该第i地震道的包络特征可重复性因子NES_i取值为0；如果第i地震道的开始数据和监控数据的相位差为180度，或其中一组数据的值全部为0，则该第i地震道的包络特征可重复性因子NES_i取值为200％；如果第i地震道的开始数据的振幅是监控数据的振幅的一半，则该第i地震道的包络特征可重复性因子NES_i取值为66.7％。

在一些可能的实施例中，终端还可以利用地震数据的包络和，计算包络特征差异化因子。具体的可以为：

终端根据差数据的包络和SED_i以及开始数据的包络和SEB_i，按照公式

计算差数据与开始数据的包络特征差异化因子RDB_i。

终端根据差数据的包络和SED_i以及监控数据的包络和SEM_i，按照公式

计算差数据与监控数据的包络特征差异化因子RDM_i。

终端根据监控数据的包络和SEM_i以及开始数据的包络和SEB_i，按照公式

计算监控数据与开始数据的包络特征差异化因子RMB_i。

其中，差数据与开始数据的包络特征差异化因子可以参考图10，差数据与监控数据的包络特征差异化因子可以参考图11，监控数据与开始数据的包络特征差异化因子可以参考图12，横坐标表示地震道，纵坐标表示包络特征差异化因子的数值。包络特征差异化因子可以用于比较第i地震道的开始数据、监控数据和差数据之间的包络和的差异大小，使地震数据之间的差异更为直观。

在步骤6中，终端根据包络特征可重复性因子NES_i，确定截断门槛值NSECut。

可选的，步骤6具体可以包括：

终端根据第一地震道的包络特征可重复性因子NES₁，按照公式NSECut＝nNES₁ 0<n<1，确定截断门槛值NSECut。

其中，n的取值不是固定的，可以根据地震数据的处理要求进行适应性调整。例如，参见图14，在图13的基础上新增了表示截断门槛值的直线，截断门槛值为第一地震道的包络特征可重复性因子的0.95倍，即n＝0.95。又例如，参见图15，在图13的基础上新增了表示截断门槛值的直线，截断门槛值为第一地震道的包络特征可重复性因子的0.90倍，即n＝0.90。

在步骤7中，终端根据截断门槛值NSECut和第i地震道的包络特征可重复性因子NES_i，对第i地震道的地震数据的可重复性进行判断，如果第i地震道的包络特征可重复性因子NES_i小于或等于截断门槛值NSECut，则判定第i地震道的地震数据具有可重复性。

可选的，步骤7具体可以包括：

根据截断门槛值NSECut与包络特征可重复性因子NES_i的值，按照以下公式，确定包络特征可重复性指示因子INSE_i；

如果INSE_i＝1，则判定第i地震道的地震数据具有可重复性，如果INSE_i＝0，则判定第i地震道的地震数据不具有可重复性。例如，参见图16和图17，分别在图14和图15的基础上，将包络特征可重复性因子大于截断门槛值的地震道的包络特征可重复性指示因子赋值为1，反之赋值为0，能够直观地看出，包络特征可重复性指示因子为1的地震道的地震数据具有可重复性。如果地震数据具有可重复性，则可以根据地震数据分析油气藏的储层条件，进而能够实现对油田的最优化开发，延长油田的开采周期，尽可能多的开采油气。

在本申请实施例中，使用上述的度量时移地震数据可重复性的方法时，先将开始数据和监控数据作差得到差数据，再利用凯泽窗希尔伯特算子分别得到开始数据、监控数据和差数据对应的凯泽窗希尔伯特道数据，进而得到地震数据的包络数据，利用包络数据计算包络和，并得到包络特征可重复性因子，最终确定截断门槛值，利用该截断门槛值将包络特征可重复性因子转化为包络特征可重复性指示因子，这样，能够根据取值为0或1的包络特征可重复性指示因子，直观地度量地震数据是否具有可重复性。

另外，为了得到地震数据的包络数据，构建了地震数据的复数道，可以利用该复数道用来提取地震数据中的振幅和相位信息，进行进一步的分析。还可以利用地震数据的包络和，计算包络特征差异化因子，使地震数据之间的差异更为直观。

以上所述仅为本申请的可选的实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所做的任何修改，等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (9)

1.一种度量时移地震数据可重复性的方法，其特征在于，所述方法包括：

步骤1，根据开始数据B_i，j和监控数据M_i，j，计算差数据D_i，j，i是地震道的顺序号，i＝1，2，…，I，I是地震道的总数，j是时间样点的顺序号，j＝1，2，…，J，J是时间样点的总数，所述开始数据B_i，j为第一次观测时获取到的地震数据，所述监控数据M_i，j为后续观测时获取到的地震数据；

步骤2，利用凯泽窗希尔伯特算子，对所述开始数据B_i，j、所述监控数据M_i，j和所述差数据D_i，j进行变换，分别得到所述开始数据的凯泽窗希尔伯特道数据U_i，j、所述监控数据的凯泽窗希尔伯特道数据V_i，j和所述差数据的凯泽窗希尔伯特道数据W_i，j；

步骤3，分别根据所述开始数据B_i，j和所述开始数据的凯泽窗希尔伯特道数据U_i，j、所述监控数据M_i，j和所述监控数据的凯泽窗希尔伯特道数据V_i，j以及所述差数据D_i，j和所述差数据的凯泽窗希尔伯特道数据W_i，j，得到所述开始数据的包络数据EB_i，j、所述监控数据的包络数据EM_i，j和所述差数据的包络数据ED_i，j；

步骤4，分别根据所述开始数据的包络数据EB_i，j、所述监控数据的包络数据EM_i，j和所述差数据的包络数据ED_i，j，得到所述开始数据的包络和SEB_i、所述监控数据的包络和SEM_i以及所述差数据的包络和SED_i；

步骤5，根据所述开始数据的包络和SEB_i、所述监控数据的包络和SEM_i以及所述差数据的包络和SED_i，得到包络特征可重复性因子NES_i；

步骤6，根据所述包络特征可重复性因子NES_i，确定截断门槛值NSECut；

步骤7，如果第i地震道的所述包络特征可重复性因子NES_i小于或等于所述截断门槛值NSECut，则判定所述第i地震道的地震数据具有可重复性；

其中，所述步骤5包括：根据所述开始数据的包络和SEB_i、所述监控数据的包络和SEM_i以及所述差数据的包络和SED_i，按照公式

得到包络特征可重复性因子NES_i。

2.根据权利要求1所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤1包括：

根据开始数据B_i，j和监控数据M_i，j，按照公式D_i,j＝B_i,j-M_i,j计算差数据D_i，j。

3.根据权利要求1所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤2包括：

步骤21，按照以下公式计算所述凯泽窗希尔伯特算子H_k，

其中，K指示所述凯泽窗希尔伯特算子的样点数，(2K+1)是所述凯泽窗希尔伯特算子的样点总数，

和X₀(α)是第一类零阶修正贝塞尔函数，α是凯泽窗函数的形状参数；

步骤22，根据所述凯泽窗希尔伯特算子H_k，按照以下公式，计算所述开始数据的凯泽窗希尔伯特道数据U_i，j；

按照以下公式，计算所述监控数据的凯泽窗希尔伯特道数据V_i，j；

按照以下公式，计算所述差数据的凯泽窗希尔伯特道数据W_i，j；

其中，当j-k≤0时，B_i，j-k、M_i，j-k和D_i，j-k的取值均为0。

4.根据权利要求1所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤3包括：

根据所述开始数据B_i，j和所述开始数据的凯泽窗希尔伯特道数据U_i，j，按照以下公式，得到所述开始数据的复数道E_i，j；

E_i,j＝B_i,j+LU_i,j

根据所述监控数据M_i，j和所述监控数据的凯泽窗希尔伯特道数据V_i，j，按照以下公式，得到所述监控数据的复数道F_i，j；

F_i,j＝M_i,j+LV_i,j

根据所述差数据D_i，j和所述差数据的凯泽窗希尔伯特道数据W_i，j，按照以下公式，得到所述差数据的复数道G_i，j；

G_i,j＝D_i,j+LW_i,j

其中，L是虚数单位，L²＝-1。

5.根据权利要求4所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤3还包括：

按照公式

对所述开始数据的复数道E_i，j取模，得到所述开始数据的包络数据EB_i，j；

按照公式

对所述监控数据的复数道F_i，j取模，得到所述监控数据的包络数据EM_i，j；

按照公式

对所述差数据的复数道G_i，j取模，得到所述差数据的包络数据ED_i，j。

6.根据权利要求1所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤4包括：

根据所述开始数据的包络数据EB_i，j，按照公式

计算所述开始数据的包络和SEB_i；

根据所述监控数据的包络数据EM_i，j，按照公式

计算所述监控数据的包络和SEM_i；

根据所述差数据的包络数据ED_i，j，按照公式

计算所述差数据的包络和SED_i。

7.根据权利要求6所述的度量时移地震数据可重复性的方法，其特征在于，在所述步骤4之后，所述方法还包括：

根据所述差数据的包络和SED_i以及所述开始数据的包络和SEB_i，按照公式

计算所述差数据与所述开始数据的包络特征差异化因子RDB_i；

根据所述差数据的包络和SED_i以及所述监控数据的包络和SEM_i，按照公式

计算所述差数据与所述监控数据的包络特征差异化因子RDM_i；

根据所述监控数据的包络和SEM_i以及所述开始数据的包络和SEB_i，按照公式

计算所述监控数据与所述开始数据的包络特征差异化因子RMB_i。

8.根据权利要求1所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤6包括：

根据第一地震道的包络特征可重复性因子NES₁，按照公式NSECut＝nNES₁ 0<n<1，确定截断门槛值NSECut。

9.根据权利要求1所述的度量时移地震数据可重复性的方法，其特征在于，所述步骤7包括：

根据所述截断门槛值NSECut与所述包络特征可重复性因子NES_i的值，按照以下公式，确定包络特征可重复性指示因子INSE_i；

Images