[发明专利]一种确定成藏期致密砂岩古孔隙度方法

说明书

技术领域

本发明涉及地质勘探技术领域，具体而言，本发明涉及一种确定成藏期致密砂岩古孔隙度方法。

背景技术

致密砂岩气因其具有储量巨大，分布广泛的特点而成为非常规天然气中最重要的一种类型，由于气藏的成藏动力学机制和成藏过程存在根本差异，因此可将致密砂岩气藏的特征分为两种，一种是流体在剖面上表现为不服从重力分异原理的气、水倒置关系，气藏的形态、大小和边界不受构造等高线的控制，即Masters所提出的深盆气藏；另一种是圈闭类型主要为背斜型气藏或地层岩性气藏，气藏的形态和气水分布通常与构造密切相关的致密砂岩。目前有“先成型”和“后成型”两种致密砂岩气藏的概念来解释上述两种特征的差异。先成型是指储层先致密化后才发生大规模的天然气聚集，成藏的根本原因是在致密储层中气藏底部缺少自由流动的地层水，致使储层中缺少对于气的浮力作用，气体在低渗砂岩中上浮受毛细管力的阻碍而导致对气的封闭，这就使得在勘探部署时需寻找储层的甜点；后成型是指天然气大规模运移发生在储层致密化之前，而后期构造作用和成岩作用才导致储层致密，气藏的成藏动力来自早期储层物性较好时的浮力作用，这就使得在勘探部署时需同时兼顾古、今构造特征。因此，在对某一气藏进行勘探部署时需要先判断其成因类型，这就需要研究成藏期致密砂岩气的古孔隙度是否已达到致密化储层的标准，到目前为止国内外学者提出了多种孔隙演化的预测方法和模型。Athy（1930）最早提出根据储层埋藏深度来预测孔隙度的关系式，Maxwell等根据实验数据推导了温度和时间对孔隙度演化的影响，Scherer考虑了储层的埋藏深度、石英颗粒的体积分数、颗粒分选性、地温梯度和沉积物年龄。Athy、Maxwell和Scherer等提出的孔隙度预测方法均只考虑了压实作用的影响，对于溶解作用、胶结作用和自生矿物的生长引起的孔隙度的变化在定量与储层的孔隙度变化模型中均没有加以讨论。以上研究均未给出基于成岩演化分析的以地质历史为背景的孔隙度演化史科学的定量化模型，很难从地质历史时期上将致密砂岩气的成藏期与储层孔隙度的演化相匹配起来，因此无法准确判断成藏期致密砂岩的古孔隙度。

发明内容

为了解决现有技术中无法准确判断成藏期致密砂岩气的古孔隙度的问题，提出了一种确定成藏期致密砂岩古孔隙度方法。

本发明实施例提供了一种确定成藏期致密砂岩古孔隙度方法，包括：

获取储层岩性、分选系数和成岩作用特征；

根据所述成岩序列中不同成岩作用特征的组合，利用初始孔隙度、孔隙度损失模型、孔隙度增加模型和裂缝孔隙度模型生成储层的孔隙度综合演化模型；

根据储层成岩作用阶段，结合储层初始孔隙度和储层孔隙度综合演化模型建立储层孔隙度演化史模型；

将待测的致密砂岩气藏的成藏时期与储层孔隙度演化模型进行匹配进而确定成藏期致密砂岩古孔隙度。

根据本发明实施例所述的一种确定成藏期致密砂岩古孔隙度方法的一个进一步的方面，所述成岩作用特征包括，压实和压溶作用、胶结作用和溶蚀作用、压实压溶产生微裂缝作用。

根据本发明实施例所述的一种确定成藏期致密砂岩古孔隙度方法的再一个进一步的方面，所述孔隙度损失模型中包括压实损失和胶结损失，其中压实损失通过以下公式得到：

Φ_C=Ф₀*e^（-C*Z）               （2）

其中，C为压实因子，Z为埋藏深度，Ф₀为初始孔隙度，Ф₀=20.91+（22.9/S₀），其中，S₀是Trask分选系数；

胶结损失是由于胶结作用所导致的孔隙度损失，其孔隙度损失量Φ_Ce大体与胶结物的含量相当。

根据本发明实施例所述的一种确定成藏期致密砂岩古孔隙度方法的另一个进一步的方面，所述孔隙度增加模型中，溶蚀增加孔隙度Φ_d是指总储集空间中溶蚀孔所占据的那部分储集空间：

溶蚀增加孔隙度Φ_d=溶蚀孔面孔率/总面孔率×实测孔隙度           （3）。

根据本发明实施例所述的一种确定成藏期致密砂岩古孔隙度方法的另一个进一步的方面，所述裂缝孔隙度是指在单位体积的岩石中裂缝体积所占的比例：

Φ_f=W_f*A_f*Q_f/V                     (4)