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

权利要求书

1.一种确定成藏期致密砂岩古孔隙度方法，其特征在于包括：

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

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

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

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

2.根据权利要求1所述的一种确定成藏期致密砂岩古孔隙度方法，其特征在于，所述成岩作用特征包括，压实和压溶作用、胶结作用和溶蚀作用、压实压溶产生微裂缝作用。

3.根据权利要求1所述的一种确定成藏期致密砂岩古孔隙度方法，其特征在于，所述孔隙度损失模型中包括压实损失和胶结损失，其中压实损失通过以下公式得到：

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

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

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

4.根据权利要求1所述的一种确定成藏期致密砂岩古孔隙度方法，其特征在于，所述孔隙度增加模型中，溶蚀增加孔隙度Φ_d是指总储集空间中溶蚀孔所占据的那部分储集空间：

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

5.根据权利要求1所述的一种确定成藏期致密砂岩古孔隙度方法，其特征在于，所述裂缝孔隙度是指在单位体积的岩石中裂缝体积所占的比例：

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

式中，Φ_f为裂缝孔隙度；W_f为平均裂缝宽度；A_f为平均裂缝面面积；Q_f为岩石内裂缝的数量；V为岩石体积。

6.根据权利要求5所述的一种确定成藏期致密砂岩古孔隙度方法，其特征在于，所述储层孔隙度综合演化模型如下：

式中，Ф为某一时期孔隙度；Ф₀为初始孔隙度；C为压实因子；Z为埋藏深度；Z_A为成岩演化阶段初始埋藏深度；Z_B为成岩演化阶段最大埋藏深度；Φ_l为岩石中绿泥石含量导致的储层的孔隙度增加；Φ_Ce为胶结作用损失的孔隙度；Φ_d为溶蚀增加孔隙度；Φ_f为裂缝孔隙度。

7.根据权利要求6所述的一种确定成藏期致密砂岩古孔隙度方法，其特征在于，在根据储层成岩作用阶段，结合储层初始孔隙度和储层孔隙度综合演化模型建立储层孔隙度演化史模型中具体包括：

针对早成岩作用的A期的储层孔隙度演化史模型为：

Ф=Ф₀*e^（-C*Z）+Z/1500*Ф_l            (6)

其中，Z：0-1500m；Z_B=1500m；

针对早成岩作用的B期的储层孔隙度演化史模型为：

Ф=Ф₀*e^（-C*Z）+Ф_l-（Z-1500）/（2500-1500）*Ф_Q+（Z-1500）/（2500-1500）*Ф_d       (8)

其中，Z：1500-2500m；Z_B=2500m；

针对晚成岩作用的A期的储层孔隙度演化史模型为：

Ф=Ф₀*e^（-C*Z）+Ф_l-（Z-2500）/（4800-2500）*Ф_Q+（Z-2500）/（4800-2500）*Ф_d-（Z-2500）/（4800-2500）*Ф_A           (9)

其中，Z：2500-4800m；Z_B=4800m；

针对晚成岩作用的B期的储层孔隙度演化史模型为：

Ф=Ф₀*e^（-C*Z）+Ф_l-Ф_Q+Ф_d-Ф_A+Ф_f            (10)

其中，Z=Z_B>4800m；

Ф为某一埋藏深度下的岩石的初始孔隙度；Ф₀为初始孔隙度；Z为埋藏深度；Z_B为成岩作用阶段的最大埋藏深度；C为压实因子；Ф_l为岩石中绿泥石含量导致的储层的孔隙度增加；Ф_Q为岩石中石英次生加大和自生石英的含量导致的孔隙度的减低量；Ф_d为溶蚀增加孔隙度；Ф_f为裂缝孔隙度；Ф_A为岩石中晚期碳酸盐胶结物的含量。