[发明专利]预测砂泥岩薄互层砂体累积厚度的方法

说明书

本发明提供一种预测砂泥岩薄互层砂体累积厚度的方法，该预测砂泥岩薄互层砂体累积厚度的方法包括：步骤1，加载地震数据；步骤2，确定提取地震属性时窗范围；步骤3，提取均方根振幅属性和质心频率属性；步骤4，建立均方根振幅、质心频率与薄互层砂体相对累积厚度的关系式；步骤5，计算砂泥岩薄互层砂体相对累积厚度并求取校正系数；步骤6，计算砂泥岩薄互层砂体累积厚度及成图。该预测砂泥岩薄互层砂体累积厚度的方法直接利用地震数据本身的振幅和频率信息，直接预测薄互层储层砂体累积厚度，避免了提频、时频分析等方法介入带来的误差，具有很高的可靠性。

技术领域

本发明涉及油气地球物理勘探技术领域，特别是涉及到一种预测砂泥岩薄互层砂体累积厚度的方法。

背景技术

陆相含油气盆地砂泥岩薄互层广泛分布。薄互层是指不同岩性的薄层交互出现的一种地质现象，它的一个重要特征是每相邻两层的岩性不同，相同岩性的薄层隔层重复出现。砂泥岩薄互层储层纵向上单砂体厚度薄，横向上砂体通透性差。这种储层地震剖面上的一个或几个同相轴，往往是一个薄互层的反射。其中砂泥岩薄互层的砂体累积厚度是储层预测的重要指标。如今众多的薄层定量描述方法通常基于单砂体楔状模型，利用薄层地震属性可以识别厚度小于调谐厚度的孤立薄层。然而，利用上述适用于孤立薄层的定量描述方法不能有效预测薄互层砂体累积厚度。薄互层砂体厚度预测已成为制约薄互层油气勘探的核心问题。

基于时频分析的谱分解技术是一种基于单砂体楔状模型的预测薄层厚度的方法，用一系列离散的频率(尤其是高频)成分对地震资料进行解析成像，通过频率扫描从整体上反映储集层的形态、厚度，揭示地层纵向变化规律及沉积相带的演变，进而进行薄砂体的厚度预测。但是该方法对砂泥岩薄互层进行累积厚度预测是会产生较大的误差。

地震属性切片技术广泛应用研究薄层结构及其沉积演化过程，由于地震反射的干涉效应，某一时刻的地震切片并不能反映该深度地层的横向变化和展布情况，并且地震属性切片解释存在多解性。

地震属性分析技术是利用地震波的反射特征对薄层结构进行间接的分析和预测。地震属性是指地震数据经过数学变化后导出的关于地震波的几何特征、运动学特征、动力学特征和统计学特征的特殊量度。其中振幅和频率是最为稳定的基本属性。由于子波和薄层结构决定了地震波的反射特征及其地震属性。因此,地震属性中的确包含着薄层的结构特征信息。但是，薄互层条件下的地震属性(振幅和频率)与薄互层砂体累积厚度之间的预测关系尚缺乏相应研究。为此我们发明了一种新的预测砂泥岩薄互层砂体累积厚度的方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种直接利用地震属性准确预测砂泥岩薄互层砂体累积厚度的方法。

本发明的目的可通过如下技术措施来实现：预测砂泥岩薄互层砂体累积厚度的方法，该预测砂泥岩薄互层砂体累积厚度的方法包括：步骤1，加载地震数据；步骤2，确定提取地震属性时窗范围；步骤3，提取均方根振幅属性和质心频率属性；步骤4，建立均方根振幅、质心频率与薄互层砂体相对累积厚度的关系式；步骤5，计算砂泥岩薄互层砂体相对累积厚度并求取校正系数；步骤6，计算砂泥岩薄互层砂体累积厚度及成图。

本发明的目的还可通过如下技术措施来实现：

在步骤1中，加载地震数据，将待分析的地震数据加载到系统中。

在步骤2中，以目标砂体所在的地震层位为时窗中心，以地震子波的周期为时窗长度，确定地震属性提取时窗。

在步骤2中，时窗选择原则是采用包含砂泥岩薄互层在内的时窗，提取振幅和频率属性对薄互层砂体累积厚度进行预测。

在步骤2中，在提取相应属性时，薄互层的位置与参考标准层在一个沉积旋回内，保证提取的属性体现相同的沉积环境，利于分析已知井同一旋回砂层变化与属性变化的最佳匹配。