[发明专利]一种碳酸盐岩储层渗透率识别方法、系统及装置

权利要求书

1.一种碳酸盐岩储层渗透率识别方法，其特征在于，该识别方法包括以下步骤：

1)利用碳酸盐岩储层岩心各点的孔隙结构参数和岩心各点的孔隙类型建立各种孔隙类型和孔隙结构参数取值范围的相互关系，其中，岩心各点的孔隙结构参数是根据岩心各点的横波测井数据、测井密度数据、孔隙度和岩石骨架剪切模量计算得到，所述孔隙类型包括粒间孔隙、粒内孔隙、混合孔和溶蚀孔洞裂缝；

2)根据岩心各点的孔隙类型建立每种孔隙类型对应的渗透率模型，所述渗透率模型用于表征渗透率与孔隙度的相关关系；

3)获取待识别碳酸盐岩储层的孔隙度，并确定待识别碳酸盐岩储层的孔隙结构参数，其中，待识别碳酸盐岩储层的孔隙结构参数是根据待识别碳酸盐岩储层的横波测井数据、测井密度数据、孔隙度和岩石骨架剪切模量计算得到，采用的计算公式为：

其中γ_μ为孔隙结构参数，为孔隙度，f_μ为岩石骨架柔韧性参数，无量纲，V_s为横波速度，ρ为体积密度，μ_s为岩石骨架剪切模量；

4)根据孔隙类型和孔隙结构参数取值范围的相互关系，利用待识别碳酸盐岩储层的孔隙结构参数确定待识别碳酸盐岩储层的孔隙类型；

5)根据待识别碳酸盐岩储层的孔隙类型选取相应的渗透率模型，将待识别碳酸盐岩储层的孔隙度代入所选取的渗透率模型，确定待识别碳酸盐岩储层的渗透率。

2.根据权利要求1所述的碳酸盐岩储层渗透率识别方法，其特征在于，当待识别碳酸盐岩储层缺少横波测井数据时，所述识别方法还包括：

利用中子测井数据、纵波测井数据和深浅侧向电阻率测试数据中的至少一种确定待识别碳酸盐岩储层的孔隙类型；

根据孔隙类型选取相应的横波正演模型，并将待识别的纵波测井数据代入所选取的横波正演模型，获得待识别碳酸盐岩储层的横波正演数据，将该横波正演数据作为横波测井数据来计算待识别碳酸盐岩储层的孔隙结构参数，所述横波正演模型是根据岩心各点的横波测井数据、岩心各点的纵波测井数据和岩心各点的孔隙类型建立的。

3.根据权利要求1所述的碳酸盐岩储层渗透率识别方法，其特征在于，当待识别碳酸盐岩储层为井间储层时，所述识别方法还包括：

利用有测井数据的井和地震数据，采用确定性反演法得到有测井数据和地震数据各点的波阻抗、弹性波阻抗和拉梅系数；

分层和分孔隙类型分别统计有测井数据和地震数据各点的波阻抗、弹性波阻抗和拉梅系数的变差函数及概率密度函数，得到各个层段的各种孔隙类型的波阻抗、弹性波阻抗和拉梅系数的变差函数及概率密度函数；

根据各个层段的各种孔隙类型的波阻抗、弹性波阻抗和拉梅系数的变差函数和概率密度函数，通过随机模拟产生井间各点的波阻抗、弹性波阻抗和拉梅系数；

根据井间各点的波阻抗、弹性波阻抗和拉梅系数获得井间各点的孔隙度和孔隙结构参数；

根据井间各点的孔隙结构参数识别井间各点的孔隙类型；

根据井间各点的孔隙类型选取相应的渗透率模型，将井间各点的孔隙度代入选取的渗透率模型，获得井间各点的渗透率。

4.根据权利要求1所述的碳酸盐岩储层渗透率识别方法，其特征在于，所述步骤1)中岩心的孔隙类型是对表征碳酸盐岩储层微观孔隙特征的影响因子进行聚类划分得到，所述的影响因子包括排驱压力、平均喉道半径、相对分选系数、偏度系数和最大进汞饱和度。