基于能量包络的地震波初至拾取方法及计算机可读存储介质与流程

本发明涉及地球物理勘探技术领域，更具体地，涉及一种基于能量包络的地震波初至拾取方法，其可应用于油气地震勘探及非常规地震勘探的资料处理中。

背景技术：

地震勘探资料中存在大量复杂地表的勘探资料，在对地震勘探资料进行处理时，其静校正问题十分复杂，层析静校正等方法的应用效果都依赖于初至拾取的质量。当前初至拾取的方法一般是采用软件自动拾取，然后通过人工予以修正。但由于地表条件复杂，初至波受多种因素影响，自动拾取的质量往往不是很高，很大程度上需要人工拾取修正，工作量巨大，处理周期延长。为此，期待一种能够提高地震波初至拾取的质量和效率的地震波初至拾取方法。

技术实现要素：

本发明的目的是提出一种基于能量包络的地震波初至拾取方法，其可克服现有的初至拾取方法工作量巨大、处理周期长的缺点。

本发明提出了一种基于能量包络的地震波初至拾取方法，包括：

对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；

去除所述首次拾取的地震波初至中的异常值，缩小所述首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；

基于所述原始地震资料计算振幅能量包络的地震道数据；

在所述二次初至拾取时窗内，对所述振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。

优选地，所述基于能量包络的地震波初至拾取方法还包括：

输入所述原始地震资料和空间属性定义，对所述原始地震资料进行滤波和增益处理。

优选地，所述基于能量包络的地震波初至拾取方法还包括：

对所述最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。

优选地，所述基于原始地震资料计算振幅能量包络的地震道数据包括：

基于所述原始地震资料获得地震道数据，对所述地震道数据的数值做平方运算，得到当前地震道的能量；

将所述地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；

将所述当前地震道的能量与所述相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。

优选地，通过炮相关方法对所述原始地震资料进行首次初至拾取。

本发明还提出一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：

对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；

去除所述首次拾取的地震波初至中的异常值，缩小所述首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；

基于所述原始地震资料计算振幅能量包络的地震道数据；

在所述二次初至拾取时窗内，对所述振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。

优选地，所述程序被处理器执行时还实现以下步骤：

输入所述原始地震资料和空间属性定义，对所述原始地震资料进行滤波和增益处理。

优选地，所述程序被处理器执行时还实现以下步骤：

对所述最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。

优选地，所述基于原始地震资料计算振幅能量包络的地震道数据包括：

基于所述原始地震资料获得地震道数据，对所述地震道数据的数值做平方运算，得到当前地震道的能量；

将所述地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；

将所述当前地震道的能量与所述相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。

优选地，通过炮相关方法对所述原始地震资料进行首次初至拾取。

本发明的有益效果在于：振幅能量包络的地震道数据的抗干扰能力增强，易于进行初至识别，因此根据本发明的基于能量包络的地震波初至拾取方法能够提高地震波初至拾取的效率和准确率，为提高后续的静校正精度奠定了良好的基础。

本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。

附图说明

通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的附图标记通常代表相同部件。

图1示出了根据本发明实施例的基于能量包络的地震波初至拾取方法的流程图。

图2示出了根据本发明实施例的地震波初至拾取方法中的地震道数据。

图3出了根据本发明实施例的地震波初至拾取方法中的地震道的能量。

图4示出了根据本发明实施例的地震波初至拾取方法中将地震道数据的相位转换90°之后的地震道数据。

图5出了根据本发明实施例的地震波初至拾取方法中的相位转换后的地震道能量。

图6示出了根据本发明实施例的地震波初至拾取方法中的振幅能量包络的地震道数据。

具体实施方式

下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。

图1示出了根据本发明实施例的地震波初至拾取方法的流程图，如图1所示，根据本发明实施例的基于能量包络的地震波初至拾取方法包括以下步骤：

步骤1：对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至。

可使用炮相关等方法对原始地震资料进行快速的初至拾取，得到首次拾取的地震波初至，首次拾取的地震波初至可为后续的二次初至拾取提供拾取时窗依据；

步骤2：去除首次拾取的地震波初至中的异常值，并缩小首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗。

针对首次拾取的地震波初至进行质量控制，去除异常值，去除异常值后，地震波初至的范围有一定程度的缩小，缩小后的范围即可作为二次初至拾取时窗，在该范围内，不包含异常值。特别地，针对相邻两个地震道的初至值发生异常变化(例如相邻两个地震道的初至值的差异超过预定值)的不合理情况，需要去除异常值。

步骤3：基于原始地震资料计算振幅能量包络的地震道数据。

可以按照以下步骤基于原始地震资料计算振幅能量包络的地震道数据：

子步骤301：基于原始地震资料获得地震道数据，对地震道数据的数值做平方运算，得到当前地震道的能量；

子步骤302：将地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；

子步骤303：将当前地震道的能量与相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。

振幅能量包络的地震道数据的抗干扰能力增强，易于进行初至识别。

步骤4：在二次初至拾取时窗内，对振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。

在步骤2得到的二次初至拾取时窗内，对振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。二次初至拾取可通过现有软件来实现。

在一个示例中，根据本发明实施例的基于能量包络的地震波初至拾取方法还包括以下步骤：

输入原始地震资料和空间属性定义，对原始地震资料进行滤波和增益处理。对原始地震资料的预先处理时在不影响初至形态的前提下进行的。空间属性定义是指在地震资料处理中地震资料与空间位置的对应关系，将地震资料与空间位置相对应，确定二者之间的对应关系，是地震资料处理中必不可少的环节，是本领域技术人员根据现有技术能够实现的。

在一个示例中，根据本发明实施例的基于能量包络的地震波初至拾取方法还包括以下步骤：

对最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。

本发明实施例还提出一种计算机可读存储介质，其上存储有计算机程序，其中，程序被处理器执行时实现以下步骤：

对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；

去除首次拾取的地震波初至中的异常值，缩小首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；

基于原始地震资料计算振幅能量包络的地震道数据；

在二次初至拾取时窗内，对振幅能量包络的地震道数据进行二次初至拾取，获得最终的地震波初至。

在一个示例中，程序被处理器执行时还实现以下步骤：

输入原始地震资料和空间属性定义，对原始地震资料进行滤波和增益处理。

在一个示例中，程序被处理器执行时还实现以下步骤：

对最终的地震波初至按炮检距关系进行分布显示，剔除异常值，然后进行输出。

在一个示例中，基于原始地震资料计算振幅能量包络的地震道数据包括：

基于原始地震资料获得地震道数据，对地震道数据的数值做平方运算，得到当前地震道的能量；

将地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量；

将当前地震道的能量与相位转换后的地震道能量相加，得到振幅能量包络的地震道数据。

在一个示例中，通过炮相关方法对原始地震资料进行首次初至拾取。

实施例

在一个实施例中，基于能量包络的地震波初至拾取方法包括以下步骤：

步骤1：对原始地震资料进行首次初至拾取，获得首次拾取的地震波初至；

步骤2：去除首次拾取的地震波初至中的异常值，缩小首次拾取的地震波初至的范围，将缩小后的范围作为二次初至拾取时窗；

步骤3：基于原始地震资料计算振幅能量包络的地震道数据，该步骤包括以下三个子步骤：

子步骤301：基于原始地震资料获得地震道数据，对地震道数据的数值做平方运算，得到当前地震道的能量，图2显示了地震道数据，图3显示了地震道的能量；

子步骤302：将地震道数据的相位转换90°，通过对相位转换后的地震道数据的数值做平方运算，得到相位转换后的地震道能量，图4显示了将地震道数据的相位转换90°之后的地震道数据，图5显示了相位转换后的地震道能量；

子步骤303：将当前地震道的能量与相位转换后的地震道能量相加，得到振幅能量包络的地震道数据，图6显示了振幅能量包络的地震道数据。

由于振幅能量包络的地震道数据的抗干扰能力增强，易于进行初至识别，因此根据本发明的基于能量包络的地震波初至拾取方法能够提高地震波初至拾取的效率和准确率，为提高后续的静校正精度奠定了良好的基础。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。