确定地震勘探中采集的地震记录的类型的方法

摘要

本发明提供一种确定地震勘探中采集的地震记录的类型的方法，所述方法包括：(A)获取地震记录，其中，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅；(B)根据所述地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间；(C)判断接收到的地震波的初至时间是否满足预定条件；(D)当接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。根据本发明的示例性实施例的方法能够准确地识别地震记录的类型。

说明书

技术领域

本发明属于地震勘探领域，更具体地说，涉及一种确定地震勘探中采集的地震记录的类型的方法。

背景技术

在地震勘探生产过程中，每进行几次单炮记录后就要进行一次背景记录(即，背景噪声的记录)。这样，获取的地震记录中就夹杂着不同类型的地震记录，一种是用来确定地下地质构造的单炮记录，另一种是用来分析背景噪声的背景记录。然而，随着地震数据采集的高效化，迫切需要对野外现场采集的地震记录进行分类以快速提取所需的地震记录。目前，现有的确定地震记录的类型的方法无法准确地识别地震记录的类型。

发明内容

本发明的示例性实施例的目的在于提供一种确定地震勘探中采集的地震记录的类型的方法，以准确地识别地震记录的类型。

本发明的示例性实施例提供一种确定地震勘探中采集的地震记录的类型的方法，所述方法包括：(A)获取地震记录，其中，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅；(B)根据所述地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间；(C)判断接收到的地震波的初至时间是否满足预定条件；(D)当接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。

可选地，所述方法还包括：(E)当接收到的地震波的初至时间不满足预定条件时，所述地震记录为背景记录。

可选地，所述预定条件为所述设置在一个排列上的全部检波器中有连续的预定数量的检波器接收到的地震波的初至时间递增或递减。

可选地，确定一个检波器接收到地震波的初至时间的步骤包括：选取所述一个检波器接收到的所有地震振幅中对应的时间位于预定时窗内的多个地震振幅；通过将预定长度的滑动时窗在所述多个地震振幅上进行滑动，将所述多个地震振幅划分为多个单元，其中，在所述滑动时窗的各个位置上的滑动时窗内的地震振幅构成一个单元，滑动步长为预定时间间隔，所述预定长度为所述预定采样间隔的奇数倍；将所述多个单元按单元的起始时间从小到大的顺序排列，分别计算排列后的多个单元中相邻的两个单元之间的能量比值，其中，能量比值为后一单元的能量值比上前一单元的能量值；确定计算的所述能量比值中最大的能量比值对应的单元，所述单元中位于中间位置的地震振幅对应的时间为所述一个检波器接收到地震波的初至时间。

可选地，所述预定时间间隔为所述预定采样间隔。

可选地，通过如下公式计算所述能量比值：

$>A_{i+1}={\left\{\right[\underset{t=T_{i+\mathrm{1,1}}}{\overset{T_{i+\mathrm{1,2}}}{\Sigma }}{x_{i+1}}^2\left(t\right)]/[\underset{t=T_{i,1}}{\overset{T_{i,2}}{\Sigma }}{x_i}^2\left(t\right)\left]\right\}}^{1/2}$>

其中，A_i+1为与第i+1单元对应的能量比值，x_i+1(t)为第i+1单元内的时间为t时的地震振幅，T_i+1,1为第i+1单元的起始时间，T_i+1,2为第i+1单元的终止时间，x_i(t)为第i单元内的时间为t时的地震振幅，T_i,1为第i单元的起始时间，T_i,2为第i单元的终止时间，i为大于0的整数，其中，第i单元与第i+1单元为相邻的两个单元，第i单元为前一单元，第i+1单元为后一单元。

可选地，在地震勘探为三维地震勘探时，步骤(B)中所述一个排列为所述所有检波器设置到的多个排列中距离所述一个炮点最近的一个检波器所在的一个排列。

根据本发明的示例性实施例的一种确定地震勘探中采集的地震记录的类型的方法，能够准确地识别地震记录的类型。

将在接下来的描述中部分阐述本发明另外的方面和/或优点，还有一部分通过描述将是清楚的，或者可以经过本发明的实施而得知。

附图说明

通过下面结合附图进行的对实施例的描述，本发明的上述和/或其它目的和优点将会变得更加清楚，其中：

图1是示出根据本发明示例性实施例的确定地震勘探中采集的地震记录的类型的方法的流程图。

具体实施方式

现将详细描述本发明的示例性实施例，所述实施例的示例在附图中示出，其中，相同的标号指示相同的部分。以下将通过参照附图来说明所述实施例，以便解释本发明。

图1是示出根据本发明示例性实施例的确定地震勘探中采集的地震记录的类型的方法的流程图。

如图1所示，在步骤S100，获取地震记录。这里，地震记录为一个炮点对应的所有检波器在预定时间内按预定采样间隔接收到的地震数据，地震数据包括地震振幅。换言之，地震记录就是在预定时间内使用对应一个炮点的所有检波器分别对地震波进行采样而得到的地震数据，采样间隔为所述预定采样间隔。所述地震波可以为人工产生的地震波(即在所述所有检波器对应的炮点产生的地震波)，也可以为由干扰源产生的地震波(即背景噪声)。

在步骤S200，根据步骤S100获取的地震记录中的地震数据，确定所述所有检波器中设置在一个排列上的全部检波器接收到地震波的初至时间。在二维地震勘探中，对应一个炮点的所有检波器都设置在一个排列中，则所述全部检波器为所述所有检波器。在三维地震勘探中，对应一个炮点的所有检波器设置在多个排列中，此时，所述全部检波器可以为所述所有检波器中的设置在多个排列中任一个排列上的检波器。优选地，在三维地震勘探时，步骤S200中所述一个排列为所述所有检波器设置到的多个排列中距离所述所有检波器对应的炮点最近的一个检波器所在的一个排列。由此，可提高后述步骤S300的执行效率，进而提高识别地震记录类型的效率。

可以理解，这里，可使用各种现有的确定初至时间的方法来确定所述全部检波器接收到地震波的初至时间。作为示例，确定一个检波器接收到地震波的初至时间的步骤包括：

首先，选取所述一个检波器接收到的所有地震振幅中对应的时间位于预定时窗内的多个地震振幅。即，从所述所有地震振幅中选取出采样时间在预定时窗内的多个地震振幅。这里，预定时窗位于所述预定时间内，可以根据经验来设定所述预定时窗的位置和所述预定时窗的长度。

然后，通过将预定长度的滑动时窗在所述多个地震振幅上进行滑动，将所述多个地震振幅划分为多个单元。这里，在所述滑动时窗的各个位置上的所述滑动时窗内的地震振幅构成一个单元，滑动步长为预定时间间隔，所述预定长度为所述预定采样间隔的奇数倍。可以理解，所述预定长度小于所述预定时窗的长度。优选地，所述预定时间间隔为所述预定采样间隔，即，所述滑动时窗每次滑动一个采样点(即，地震振幅)。所述预定时间间隔还可设定为所述预定采样间隔的大于1的整数倍。另外，所述滑动可以从所述多个地震振幅中采样时间最小的地震振幅开始滑动到采样时间最大的地震振幅，也可以与前述相反从所述多个地震振幅中采样时间最大的地震振幅开始滑动到采样时间最小的地震振幅。

其次，将所述多个单元按单元的起始时间从小到大的顺序排列，分别计算排列后的多个单元中相邻的两个单元之间的能量比值。这里，能量比值为后一单元的能量值比上前一单元的能量值。这里，所述单元的起始时间为所述单元所对应的位置上的滑动时窗的起始时间。这里，可通过现有的各种方法计算所述能量比值。作为示例，可通过如下公式计算所述能量比值：

$>A_{i+1}={\left\{\right[\underset{t=T_{i+\mathrm{1,1}}}{\overset{T_{i+\mathrm{1,2}}}{\Sigma }}{x_{i+1}}^2\left(t\right)]/[\underset{t=T_{i,1}}{\overset{T_{i,2}}{\Sigma }}{x_i}^2\left(t\right)\left]\right\}}^{1/2}---\left(1\right)$>

这里，A_i+1为与第i+1单元对应的能量比值，x_i+1(t)为第i+1单元内的时间(即，对应的时间)为t时的地震振幅，T_i+1,1为第i+1单元的起始时间，T_i+1,2为第i+1单元的终止时间，x_i(t)为第i单元内的时间(即，对应的时间)为t时的地震振幅，T_i,1为第i单元的起始时间，T_i,2为第i单元的终止时间，i为大于0的整数。这里，第i单元与第i+1单元为相邻的两个单元，第i单元为前一单元，第i+1单元为后一单元。

最后，确定计算的所述能量比值中最大的能量比值对应的单元，所述单元中位于中间位置的地震振幅对应的时间为所述一个检波器接收到地震波的初至时间。这里，能量比值对应的单元为计算所述能量比值时的后一单元。

在步骤S300，判断步骤S200接收到的地震波的初至时间是否满足预定条件。这里，所述预定条件可以为任何能够确定所述地震记录的类型的条件。优选地，所述预定条件为所述设置在一个排列上的全部检波器中有连续的预定数量的检波器接收到的地震波的初至时间递增或递减。这里，所述预定数量可以根据试验来确定。举例说来，在所述预定数量设定为5时，则所述预定条件为所述一个排列中有5个连续的检波器所接收到地震波的初至时间递增或递减。

具体地，在步骤S300进行所述判断时，要依次比较所述一个排列上的全部检波器接收到的地震波的初至时间，当判断接收到的地震波的初至时间满足预定条件时，无论是否已经对所述全部检波器接收到的地震波的初至时间都进行比较，都可停止比较，并进入下一步。

在步骤S400，当步骤S200接收到的地震波的初至时间满足预定条件时，所述地震记录为单炮记录。

在步骤S500，当步骤S200接收到的地震波的初至时间不满足预定条件时，所述地震记录为背景记录。

根据本发明的示例性实施例的一种确定地震勘探中采集的地震记录的类型的方法，能够准确地识别地震记录的类型。

根据本发明的上述方法可以被实现为计算机可读记录介质中的计算机代码。本领域技术人员可以根据对上述方法的描述来实现所述计算机代码。当所述计算机代码在计算机中被执行时实现本发明的上述方法。

本发明的以上实施例仅仅是示例性的，而本发明并不受限于此。本领域技术人员应该理解：在不脱离本发明的原理和精神的情况下，可对这些实施例进行改变，其中，本发明的范围在权利要求及其等同物中限定。