一种基于岩心图像的非常规储层渗透率评价方法技术

本发明专利技术提供一种基于岩心图像的非常规储层渗透率评价方法，包括如下步骤：使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样划分为相应的特征；分别提取每种特征的频谱，然后以混合分形方法对分析结果进行拟合，得到各相应特征的分形参数；将各拟合后的分形参数进行组合，得到基于混合分形理论的岩石孔隙模型；通过孔隙连通性算法对岩石孔隙模型中不同种类孔隙之间的连通概率按连通概率矩阵排列后进行计算，即可得到当前非常规储层的表观渗透率的评估结果。本发明专利技术从真实岩样数据出发，结合孔隙连通性算法，得到了不同种类孔隙之间的连通概率矩阵，最终实现快速准确地评估非常规储层表观渗透率。

【技术实现步骤摘要】
一种基于岩心图像的非常规储层渗透率评价方法
本专利技术属于石油开采领域，特别是涉及一种基于岩心图像的能够快速准确地评估非常规储层表观渗透率的评价方法。
技术介绍
页岩气是从页岩储层中开采出的天然气，成分以甲烷为主，是一种重要的新兴能源。然而成功开发页岩气需要解决很多难题，其中的基础问题之一是快速而有效地评估页岩的表观渗透率。影响页岩气表观渗透率的因素有很多，主要有以下几点：第一，页岩非常致密，孔隙度和渗透率低，渗流空间小，主要孔隙直径在纳米量级；第二，渗流空间种类繁多，连通性较差，尺度不一，包括纳米级有机孔，微米级裂缝无机孔；第三，气体渗流和赋存机理复杂，包括游离态(以自由气形式存在于岩石孔隙与裂隙中)，吸附态(吸附于岩石孔隙壁面)和溶解态(溶解于干酪根、沥青质和液态原油中)等。现有的各种技术方案对非常规储层样品的渗透率计算时计算量巨大，需要建立一种快速评估非常规储层样品渗透率的新方法。目前表征岩石孔隙空间的方法有毛管束法和孔隙空间重构法，其中毛管束法是表征岩石孔隙空间的经典方法之一。该方法将岩石孔隙简化为一系列平行的毛管束，以研究流体在高孔隙度岩石渗流过程中的渗流过程。该方法模型简单，计算简便。在此基础上发展了添加迂曲度、更改毛管束直径以及增加滑移项等一系列方程，以模拟不同种类岩石孔隙空间和流体流动过程。但毛管束法的缺点在于未考虑岩石中复杂的孔隙网络结构，也未考虑孔喉连通性，模型的建立基于毛管束理论模型，这类模型的孔隙网络组成与真实岩石有一定差距。而孔隙空间重构法则基于统计来表征岩石孔隙空间，该方法首先统计岩石中的孔隙分布频谱以及配位数。进一步建立一个网状联通管道组成的网络，管道模拟喉道，两个管道的交点模拟孔隙，将统计得到的孔径按比例分散进网络的交点上，然后按统计的配位数随机封堵管道，最终完成孔隙空间的重构。该方法的缺点在于将复杂的孔隙空间简化为标准的方形网状结构模型后，打散后随机分散进网络，取消了孔隙和喉道在空间上分布的独立性。此外，喉道的迂曲度、以及喉道本身的孔径分布也未加以考虑。
技术实现思路
本专利技术的目的是提供及一种基于岩心图像的能够快速准确地评估非常规储层表观渗透率的评价方法。特别地，本专利技术提供一种基于岩心图像的非常规储层渗透率评价方法，包括如下步骤：步骤100，使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样的有机质块体和矿石颗粒类型划分为组分特征，将有机质块体，矿石颗粒和岩石基质中的孔隙信息划分为孔隙特征，将岩样中的其它特征划分为辅助特征；步骤200，分别提取每种特征的频谱，然后以混合分形方法对分析结果进行拟合，得到各相应特征的分形参数；步骤300，将各拟合后的分形参数进行组合，得到基于混合分形理论的岩石孔隙模型；步骤400，通过孔隙连通性算法对岩石孔隙模型中不同种类孔隙之间的连通概率按连通概率矩阵排列后进行计算，即可得到当前非常规储层的表观渗透率的评估结果。在本专利技术的一个实施方式中，所述步骤100中的所述孔隙特征是指有机质块体、矿石颗粒中的有机质孔隙和矿石中的粒间孔或粒内孔；所述辅助特征是指所述岩样中的天然裂缝或人工裂缝特征。在本专利技术的一个实施方式中，所述的岩石成像技术包括SEM，FIB-SEM，Micro-CT,nano-CT。在本专利技术的一个实施方式中，所述组分特征根据矿石颗粒的类型不同而细分为多个对应的子组分特征；所述孔隙特征对应所述子组分特征的不同而细分为多个对应的子孔隙特征；所述辅助特征根据与不同所述子组分特征的关系而细分为多个对应的子辅助特征。在本专利技术的一个实施方式中，所述组分特征的频谱是建立在其中各子组分特征中各自的子组分的等效大小分布频谱上，所述孔隙特征的频谱是建立在各子孔隙特征的等效孔径的分布频谱上，所述辅助特征的频谱是建立在子辅助特征所表示的裂缝密度、缝长和缝宽上。在本专利技术的一个实施方式中，每个所述子组分特征的分形参数得到的每级子组分特征的等效大小和其对应个数共同构成的大小分布频谱，等同于该子组分特征的大小分布频谱；每个所述子孔隙特征的分形参数得到的每级孔隙的等效直径和其对应个数共同构成的孔径大小分布频谱，等同于该子孔隙特征的孔径大小分布频谱；每个所述子辅助特征的分形参数得到的每级子辅助特征的等效大小和其对应个数共同构成的裂缝频谱，与等同于子辅助特征的裂缝频谱。在本专利技术的一个实施方式中，使用连通孔隙概率算法建立连通概率矩阵的方式如下：将三种特征中的各子特征按照等效边长或等效直径的方式从大到小排序，其中，第一列为等效边长或等效直径大小，第二列为等效边长或等效直径所对应的个数，第三列连通概率，第四列为该行组分类型。在本专利技术的一个实施方式中，在排序前需要进行聚并操作。在本专利技术的一个实施方式中，计算渗透率的步骤如下：其中，迭代次数i，分形标度F，第i次迭代产生的孔隙/块体边长λi，最大孔隙/块体边长λmax，第i次迭代中产生的孔隙/块体个数Ni，第i次迭代中不参与迭代的孔隙/块体个数Ni-soild，A为岩心表征单元面的截面积，Ci为第i行块体/孔隙形状因子，ki(pavg)为在平均圧力pavg下第i行块体/孔隙的渗透率。在本专利技术的一个实施方式中，ki(pavg)的计算公式需要根据矩阵中各行的类型进行判断；①设，组分特征为A，其子组分特征分别以A1、A2……An的顺序依次编号，孔隙特征为B，其子孔隙特征以B1、B2……Bn的顺序依次编号，辅助特征为C，其子辅助特征以C1、C2……Cn的顺序依次编号；组分特征A中，每一个组成子组分特征的块体都为多孔介质；当连通概率矩阵中第i行第四列记录该行为组分特征A时，ki(pavg)表达式为：其中，knD(j,pavg)为平均圧力pavg下，第j级子孔隙特征的渗透率值；计算公式如下：当j＝1时当j＝2时当j＝3时以此类推，当j＝n时其中，λ1,λ2,…,λj,…,λn为被聚并操作后的被包含在第i行所代表的多孔介质中的孔隙的等效直径；N1,N2,…,Nj,…,Nn为被聚并操作后的被包含在第i行所代表的多孔介质中的孔隙的个数；其中λ1≥λ2≥…≥λj≥…≥λn；②当连通概率矩阵中第i行第四列记录该行为未被聚并的孔时，渗透率项ki(pavg)为③当连通概率矩阵中第i行第四列记录该行为子辅助特征C1时，渗透率项ki(pavg)为④所述公式(1)中基质渗透率项kmatrix为以上各式中，μ为气体粘度，M为气体摩尔数，R为理想气体常数，T为温度，ρavg为平均密度，pavg为平均圧力，α为切向动量调节系数，为n种基质矿物百分比，k1,k2,…,kn为n种基质矿物的表观渗透率。本专利技术从真实岩样数据出发，提取主要组成成分，然后以混合分形理论为基础，针对每一种组分，分别使用分形参数拟合，最后将各组分进行组合，得到基于混合分形理论的岩石孔隙模型。进一步结合孔隙连通性算法，得到了不同种类孔隙之间的连通概率矩阵，最终构建一种非常规储层渗透率评价方法，实现快速准确地评估非常规储层表观渗透率。附图说明图1是本专利技术一个实施方式评价方法流程图。具体实施方式如图1所示，本专利技术一个实施方式公开的非常规储层渗透率评价方法一般性包括如下步骤：步骤100，使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样的有机质块体和矿石颗粒类型划分为本文档来自技高网...

【技术保护点】
1.一种基于岩心图像的非常规储层渗透率评价方法，其特征在于，包括如下步骤：步骤100，使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样的有机质块体和矿石颗粒类型划分为组分特征，将有机质块体，矿石颗粒和岩石基质中的孔隙信息划分为孔隙特征，将岩样中的其它特征划分为辅助特征；步骤200，分别提取每种特征的频谱，然后以混合分形方法对分析结果进行拟合，得到各相应特征的分形参数；步骤300，将各拟合后的分形参数进行组合，得到基于混合分形理论的岩石孔隙模型；步骤400，通过孔隙连通性算法对岩石孔隙模型中不同种类孔隙之间的连通概率按连通概率矩阵排列后进行计算，即可得到当前非常规储层的表观渗透率的评估结果。

【技术特征摘要】
1.一种基于岩心图像的非常规储层渗透率评价方法，其特征在于，包括如下步骤：步骤100，使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样的有机质块体和矿石颗粒类型划分为组分特征，将有机质块体，矿石颗粒和岩石基质中的孔隙信息划分为孔隙特征，将岩样中的其它特征划分为辅助特征；步骤200，分别提取每种特征的频谱，然后以混合分形方法对分析结果进行拟合，得到各相应特征的分形参数；步骤300，将各拟合后的分形参数进行组合，得到基于混合分形理论的岩石孔隙模型；步骤400，通过孔隙连通性算法对岩石孔隙模型中不同种类孔隙之间的连通概率按连通概率矩阵排列后进行计算，即可得到当前非常规储层的表观渗透率的评估结果。2.根据权利要求1所述的非常规储层渗透率评价方法，其特征在于，所述步骤100中的所述孔隙特征是指有机质块体、矿石颗粒中的有机质孔隙和矿石中的粒间孔或粒内孔；所述辅助特征是指所述岩样中的天然裂缝或人工裂缝特征。3.根据权利要求1所述的非常规储层渗透率评价方法，其特征在于，所述的岩石成像技术包括SEM，FIB-SEM，Micro-CT,nano-CT。4.根据权利要求3所述的非常规储层渗透率评价方法，其特征在于，所述组分特征根据矿石颗粒的类型不同而细分为多个对应的子组分特征；所述孔隙特征对应所述子组分特征的不同而细分为多个对应的子孔隙特征；所述辅助特征根据与不同所述子组分特征的关系而细分为多个对应的子辅助特征。5.根据权利要求4所述的非常规储层渗透率评价方法，其特征在于，所述组分特征的频谱是建立在其中各子组分特征中各自的子组分的等效大小分布频谱上，所述孔隙特征的频谱是建立在各子孔隙特征的等效孔径的分布频谱上，所述辅助特征的频谱是建立在子辅助特征所表示的裂缝密度、缝长和缝宽上。6.根据权利要求5所述的非常规储层渗透率评价方法，其特征在于，每个所述子组分特征的分形参数得到的每级子组分特征的等效大小和其对应个数共同构成的大小分布频谱，等同于该子组分特征的大小分布频谱；每个所述子孔隙特征的分形参数得到的每级孔隙的等效直径和其对应个数共同构成的孔径大小分布频谱，等同于该子孔隙特征的孔径大小分布频谱；每个所述子辅助特征的分形参数得到的每级子辅助特征的等效大小和其对应个数共同构成的裂缝频谱，与等同于子辅助特征的裂缝频谱。7.根据权利要求6所述的非常规储层渗...

【专利技术属性】
技术研发人员：李曹雄，林缅，江文滨，姬莉莉，曹高辉，
申请(专利权)人：中国科学院力学研究所，
类型：发明
国别省市：北京,11