碳酸盐岩溶洞纵向尺度计算方法及系统

技术领域

本发明涉及地震资料储层预测领域，更具体地，涉及一种碳酸盐岩溶洞纵向尺度计算方法及系统。

背景技术

碳酸盐岩缝洞型储层在地震上表现为明显的串珠状反射特征，由于该类储层埋藏很深，地震资料的应用存在一定的限制，加大了储层描述的难度。目前利用地震技术，能对溶洞进行定性的描述，但是由于溶洞的尺度相对较小，精确描述仍然存在很大的困难，这给此类储层的勘探开发造成了较大的影响。因此，有必要开发一种碳酸盐岩溶洞纵向尺度计算方法及系统。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种碳酸盐岩溶洞纵向尺度计算方法及系统，其通过将叠前分频振幅属性与理论实践振幅解释图版结合，实现溶洞纵向尺度的精确计算，提高储层的预测精度，降低勘探风险。

根据本发明的一方面，提出了一种碳酸盐岩溶洞纵向尺度计算方法。所述方法可以包括：根据分频数据体，获得分频振幅数组；通过正演模拟，获得不同主频的厚度-振幅解释曲线；根据厚度-振幅解释曲线，确定不同主频对应的理论振幅值，进而获得理论振幅数组；针对所述理论振幅数组与所述分频振幅数组进行相关计算，确定所述溶洞的纵向尺度。

优选地，根据分频数据体，获得分频振幅数组包括：根据所示分频数据体，计算溶洞的均方根振幅数据，进而获得所示分频振幅数组。

优选地，通过公式(1)计算所述均方根振幅数据：

其中，a为地震数据，i为采样点序号，N为地震数据总采样点数。

优选地，通过公式(2)进行相关计算：

其中，X为代表溶洞特征的分频振幅数组，Y

优选地，以相关性最大处对应的地层厚度为所述溶洞的纵向尺度。

根据本发明的另一方面，提出了一种碳酸盐岩溶洞纵向尺度计算系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：根据分频数据体，获得分频振幅数组；通过正演模拟，获得不同主频的厚度-振幅解释曲线；根据厚度-振幅解释曲线，确定不同主频对应的理论振幅值，进而获得理论振幅数组；针对所述理论振幅数组与所述分频振幅数组进行相关计算，确定所述溶洞的纵向尺度。

优选地，根据分频数据体，获得分频振幅数组包括：根据所示分频数据体，计算溶洞的均方根振幅数据，进而获得所示分频振幅数组。

优选地，通过公式(1)计算所述均方根振幅数据：

其中，a为地震数据，i为采样点序号，N为地震数据总采样点数。

优选地，通过公式(2)进行相关计算：

其中，X为代表溶洞特征的分频振幅数组，Y

优选地，以相关性最大处对应的地层厚度为所述溶洞的纵向尺度。

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

附图说明

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

图1示出了根据本发明的碳酸盐岩溶洞纵向尺度计算方法的步骤的流程图。

图2示出了根据本发明的一个实施例的厚度-振幅解释曲线的示意图。

图3示出了根据本发明的一个实施例的分频振幅数组与理论振幅数组的振幅属性叠合的示意图。

图4示出了根据本发明的一个实施例的理论振幅数组与分频振幅数组的相关性曲线的示意图。

具体实施方式

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

图1示出了根据本发明的碳酸盐岩溶洞纵向尺度计算方法的步骤的流程图。

在该实施例中，根据本发明的碳酸盐岩溶洞纵向尺度计算方法可以包括：步骤101，根据分频数据体，获得分频振幅数组；步骤102，通过正演模拟，获得不同主频的厚度-振幅解释曲线；步骤103，根据厚度-振幅解释曲线，确定不同主频对应的理论振幅值，进而获得理论振幅数组；步骤104，针对理论振幅数组与分频振幅数组进行相关计算，确定溶洞的纵向尺度。

在一个示例中，根据分频数据体，获得分频振幅数组包括：根据所示分频数据体，计算溶洞的均方根振幅数据，进而获得所示分频振幅数组。

在一个示例中，通过公式(1)计算均方根振幅数据：

其中，a为地震数据，i为采样点序号，N为地震数据总采样点数。

在一个示例中，通过公式(2)进行相关计算：

其中，X为代表溶洞特征的分频振幅数组，Y

在一个示例中，以相关性最大处对应的地层厚度为溶洞的纵向尺度。

具体地，根据本发明的碳酸盐岩溶洞纵向尺度计算方法可以包括：

根据分频数据体，通过公式(1)计算溶洞的均方根振幅数据，进而获得分频振幅数组；通过正演模拟，获得不同主频的厚度-振幅解释曲线；根据厚度-振幅解释曲线，确定地层厚度为H

本方法通过将叠前分频振幅属性与理论实践振幅解释图版结合，实现溶洞纵向尺度的精确计算，提高储层的预测精度，降低勘探风险。

为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。

在西北某探区，针对碳酸盐岩溶洞开展了针对性目标处理。通过叠前处理，得到4套分频数据，频带宽度分别为8-20Hz、21-30Hz、31-41Hz、42-100Hz。提取T74层以下80ms均方根振幅属性，选取典型井A，并统计井点处均方根振幅值，如表1所示。

表1

图2示出了根据本发明的一个实施例的厚度-振幅解释曲线的示意图。

图3示出了根据本发明的一个实施例的分频振幅数组与理论振幅数组的振幅属性叠合的示意图。

图4示出了根据本发明的一个实施例的理论振幅数组与分频振幅数组的相关性曲线的示意图。

根据分频数据体，通过公式(1)计算溶洞的均方根振幅数据，进而获得分频振幅数组；通过正演模拟，获得主频分别为15Hz、25Hz、35Hz和45Hz的厚度-振幅解释曲线，如图2所示；根据厚度-振幅解释曲线，确定地层厚度为H

综上所述，本发明通过将叠前分频振幅属性与理论实践振幅解释图版结合，实现溶洞纵向尺度的精确计算，提高储层的预测精度，降低勘探风险。

本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。

根据本发明的实施例，提供了一种碳酸盐岩溶洞纵向尺度计算系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：根据分频数据体，获得分频振幅数组；通过正演模拟，获得不同主频的厚度-振幅解释曲线；根据厚度-振幅解释曲线，确定不同主频对应的理论振幅值，进而获得理论振幅数组；针对理论振幅数组与分频振幅数组进行相关计算，确定溶洞的纵向尺度。

在一个示例中，根据分频数据体，获得分频振幅数组包括：根据所示分频数据体，计算溶洞的均方根振幅数据，进而获得所示分频振幅数组。

在一个示例中，通过公式(1)计算均方根振幅数据：

其中，a为地震数据，i为采样点序号，N为地震数据总采样点数。

在一个示例中，通过公式(2)进行相关计算：

其中，X为代表溶洞特征的分频振幅数组，Y

在一个示例中，以相关性最大处对应的地层厚度为溶洞的纵向尺度。

本系统通过将叠前分频振幅属性与理论实践振幅解释图版结合，实现溶洞纵向尺度的精确计算，提高储层的预测精度，降低勘探风险。

本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。

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