[发明专利]一种致密砂岩储层孔隙度的预测方法

说明书

本发明公开了一种致密砂岩储层孔隙度的预测方法，包括：确定砂岩的纵波阻抗三维数据体；确定划分薄储层和厚储层的纵波阻抗门槛值；获得厚储层的孔隙度的预测结果；获得薄储层的孔隙度的预测结果；将厚储层的孔隙度和薄储层的孔隙度的预测结果合并成一个数据体，得到全储层的孔隙度的预测结果。本发明公开的致密砂岩储层孔隙度的预测方法，能够有效预测厚储层的孔隙度分布，兼顾薄储层的孔隙度分布，经实钻井验证具有较高的符合率，实现了致密砂岩储层孔隙度的高精度预测。通过分类型预测的方式，使储层孔隙度的刻画更加准确，为致密砂岩储层高效开发提供有效成果。

技术领域

本发明属于油气勘探开发领域，尤其涉及一种致密砂岩储层孔隙度的预测方法。

背景技术

鄂尔多斯盆地东缘致密砂岩气勘探前景广阔，属于非常规气藏。天然气主要储集在致密砂岩的相对高孔储层(孔隙度一般为7-15％)中。储层孔隙度的高精度预测对于在致密砂岩中寻找相对高孔的甜点尤为重要。孔隙度与纵波阻抗有一定的相关性，常规储层(孔隙度一般为12-25％)可以利用地震数据反演计算纵波阻抗，建立纵波阻抗与孔隙度线性拟合关系，进而对砂岩储层孔隙度进行预测。

然而影响纵波阻抗的因素众多，例如，致密砂岩储层孔隙度较低，因而纵波阻抗变化范围小；气层纵波阻抗通常低于差气层和干层，因而影响纵波阻抗的因素还包括含气性等；同时，由于以I类气层为主的I类储层与以II类气层为主的II类储层斜率不同，因而不同类型储层孔隙度与纵波阻抗拟合关系差异较大……以上多种原因，共同导致了应用纵波阻抗与孔隙度关系进行线性拟合预测孔隙度误差较大。此外，致密砂岩储层往往分布于陆表海缓坡沉积环境的三角洲与水下分流河道砂体中，多层砂体复杂叠置，单层砂体有效厚度仅为3-10m，储层厚度薄。反演频带与地震有效频带密切相关，通常为地震有效频带的最高频以下，例如0-80Hz，薄储层纵波阻抗在反演频带内与干层纵波阻抗存在叠置，直接应用纵波阻抗线性拟合不利于薄储层孔隙度的准确预测。

因此，根据纵波阻抗与孔隙度线性拟合关系对砂岩储层孔隙度进行预测的方法不适用于致密砂岩储层孔隙度预测。

发明内容

本发明的目的在于提供一种致密砂岩储层孔隙度的预测方法，用以解决目前致密砂岩储层孔隙度不能有效预测的问题。

本发明提供一种致密砂岩储层孔隙度的预测方法，包括：

确定砂岩的纵波阻抗三维数据体；

确定划分薄储层和厚储层的纵波阻抗门槛值b；

根据低于所述划分薄储层和厚储层的纵波阻抗门槛值b的厚储层所对应的纵波阻抗三维数据体，采用分斜率线性拟合的方式，获得厚储层的孔隙度的预测结果；

根据高于所述划分薄储层和厚储层的纵波阻抗门槛值b的薄储层所对应的纵波阻抗三维数据体，采用贝叶斯判别分析和协模拟的方式，获得薄储层的孔隙度的预测结果；

将所述厚储层的孔隙度和所述薄储层的孔隙度的预测结果合并成一个数据体，得到全储层的孔隙度的预测结果。

具体地，所述确定砂岩的纵波阻抗三维数据体具体包括以下步骤：

步骤S11：根据目标区的测井曲线绘制纵横波速度比与纵波阻抗的交会图，以岩相曲线为色标，确定砂岩对应的纵横波速度比的门槛值a，低于所述砂岩对应的纵横波速度比的门槛值a为砂岩，高于所述砂岩对应的纵横波速度比的门槛值a为泥岩；

步骤S12：通过叠前约束稀疏脉冲反演技术计算纵波阻抗三维数据体和纵横波速度比三维数据体；

步骤S13：将所述纵波阻抗三维数据体和所述纵横波速度比三维数据体进行交会，并根据交会结果确定纵横波速度比低于所述砂岩对应的纵横波速度比的门槛值a时所对应的纵波阻抗的三维空间分布范围，并将其输出的数据确定为砂岩的纵波阻抗三维数据体。

具体地，所述确定划分薄储层和厚储层的纵波阻抗门槛值b具体包括以下步骤：