一种针对非常规储层样品的相对渗透率曲线计算方法与流程

本发明属于石油开采领域，特别是涉及一种基于岩心成像的获取非常规储层样品的相对渗透率曲线的计算方法。

背景技术

非常规储层作为一种有潜力的勘探目的层，具有广阔的开发前景。以页岩、致密砂岩为代表的非常规储层孔隙结构复杂，孔隙种类众多，不同尺度特征不一，非均质性较强，流体在非常规储层中运移规律复杂。特别是在开发过程中，由于钻井液、压裂液、地层水和非常规储层中的油气资源相互影响，非常规储层中的多相流一直是国内外学者和工程师关注的重点和难点。现有的各种技术方案对计算量巨大，需要建立一种快速计算非常规储层样品相对渗透率的新方法，以探寻非常规储层两相流的渗流规律。

技术实现要素：

本发明的目的是提供一种基于岩心成像的获取非常规储层样品的相对渗透率曲线的计算方法。

特别地，本发明提供一种针对非常规储层样品的相对渗透率曲线计算方法，包括如下步骤：

步骤100，使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样的矿物颗粒类型划分为组分特征，将矿物颗粒中的粒内孔，矿物颗粒间的粒间孔信息划分为孔隙特征，将岩样中的其它特征划分为辅助特征；

步骤200，将组分特征根据矿物颗粒的类型不同而细分为a1、a2……an个对应的子组分特征；将孔隙特征对应所述子组分特征的不同而细分为b1、b2……bn个对应的子孔隙特征；将辅助特征根据与不同所述子组分特征的关系而细分为c1、c2……cn个对应的子辅助特征；利用混合分形方法分别拟合每种子组分特征的分形参数；

步骤300，建立立方体构成的表征单元体，将表征单元体整体作为子组分特征a1，再将a2、a3按照分形重构的块体大小和个数顺序放进表征单元体，a2、a3的块体所占的几何空间的组分信息替换掉a1的组分信息，再将子孔隙特征b1随机放进表征单元体，统计每个孔的壁面的最大组成成分，作为该孔隙的组成成分；

步骤400，设定模型为油驱水，油水界面张力30mn/m，油粘度1.39mpa.s，水粘度1.05mpa.s，岩心表观渗透率1000md，孔隙度22％，接触角0°，建立连通概率矩阵，通过计算即可得到该模型的相渗曲线图。

在本发明的一个实施方式中，对于所述组分特征，如果子组分特征a1、a2……an的块体大小分布符合分形特征，则先提取a1的块体大小分布频谱，使用混合分形方法对a1的块体大小分布频谱进行拟合，得到相应a1的混合分形参数af1，使用混合分形参数af1重构出a1的块体大小分布频谱得到等同于与a1的块体大小分布频谱；依次对a2、a3做同样处理。

在本发明的一个实施方式中，对于所述组分特征，如果子组分特征a1、a2……an的块体大小分布不符合分形特征，则直接统计a1、a2、……an的面积或体积占比，即ai的占比为

其中为ai的面积或体积占比，vai为ai的面积或体积，j为整数。

在本发明的一个实施方式中，对于所述孔隙特征，提取子孔隙特征b1的孔隙等效直径分布频谱，使用混合分形方法对b1的孔隙等效直径分布频谱进行拟合，得到b1的混合分形参数bf1，使用混合分形参数bf1重构出b1的孔隙等效直径分布频谱，得到等同于b1的孔隙等效直径分布频谱；依次对b2、b3做同样处理。

在本发明的一个实施方式中，所述的分形参数包括迭代次数i，分形标度f，第i次迭代产生的孔隙/块体边长λi，最大孔隙/块体边长λmax，第i次迭代中产生的孔隙/块体个数ni，第i次迭代中不参与迭代的孔隙/块体个数ni-soild。

在本发明的一个实施方式中，所述a2、a3采用随机的方式放进表征单元体空间内，且其中任意块体占用空间互相不重叠；

所述b1采用随机的方式放进表征单元体空间内，且其中所述b1所占的空间可以和任意块体占用的空间相互重合；

统计所述b1所占的空间中各个子组成特征a1、a2……an的占比，以赋予所述b1最大占比的组成成分所对应的矿物性质。

在本发明的一个实施方式中，所述b1采用随机的方式放进表征单元体空间内，统计所述b1所占的空间中各个子组成特征a1、a2……an的面积和体积所占比例，等比例赋予所述b1壁面矿物性质；所述壁面性质包括壁面润湿性和壁面粗糙度。

在本发明的一个实施方式中，将所述b1和所述c1综合，按照等效直径从大到小排序，等效直径相同时，按对应毛管力从大到小排序，构造连通概率矩阵，所述连通概率矩阵第一列为等效直径，第二列为等效直径所对应的个数，第三列连通概率，第四列为该行驱扫概率，第五列为该行导通概率，第六列为该行壁面性质。

在本发明的一个实施方式中，计算所述相渗曲线的公式为：

其中，每级孔对应的等效孔径为λ1,λ2,…,λi,…,λm，每级孔对应的个数分别为n1,n2,…,ni,…,nm；在排序后的连通概率矩阵中使用连通孔隙概率算法得到第i级连通概率p1,i；模型为非润湿相驱替润湿相，在第i级的ni各孔中，驱扫概率为p2,i，使每级孔中驱扫成功的孔隙个数为p2,i×ni，驱扫概率p2,i通过真实岩心压汞曲线和孔隙成像结果综合分析给出；第i级孔导通概率为p3,i，krw(j)为润湿相相对渗透率，kro(j)为非润湿相相对渗透率，srw(j)为润湿相饱和度，ci为第i级孔隙形状因子；其中1≤j≤m且j为整数。

在本发明的一个实施方式中，当驱替动力大于不等径孔隙连接处阻力时，孔隙导通(p3,i＝1)，否则孔隙不能导通(p3,i＝0)，影响p3,i的参数包括毛管力、润湿性、驱替压力、孔隙等效直径和界面张力参数。

本发明是基于真实岩石的成像技术，利用岩石孔隙空间的分形特性，使用混合分形单元法提取分形特征，将孔隙壁面赋予不同润湿参数，使用连通概率，驱扫概率，导通概率共同构建基于三维混合分形的两相流模型。该模型能够有效表征不同尺度下的非常规岩石样品孔隙结构特征，能够快速计算非常规岩石样样品相渗曲线。与传统方法相比，该方法涵盖信息多(孔隙分布、润湿性)、计算速度快，便于工程应用，为获得非常规储层相对渗透率开辟了一种新方法。

附图说明

图1是本发明一个实施方式评价方法流程图；

图2是本发明一个实施方式得到的相渗曲线图。

具体实施方式

如图1所示，本发明一个实施方式公开一种针对非常规储层样品的相对渗透率曲线计算方法，包括如下步骤：

步骤100，使用岩石成像技术对真实岩样进行扫描和识别后得到相应图像，将图像中岩样的矿物颗粒类型划分为组分特征，将矿物颗粒中的粒内孔，矿物颗粒间的粒间孔信息划分为孔隙特征，将岩样中的其它特征划分为辅助特征；

这里的岩石成像技术包括但不限于sem(扫描电子显微镜)，fib-sem(聚焦离子束)，micro-ct(微计算机断层扫描技术),nano-ct、背散射扫描电镜等技术。

在上述各成像技术中，既可以只采用一种方式成像后进行分析，也可以同时采用两种以上的成像技术进行结合分析，因为有的成像技术观测范围广，有的成像技术精度高，可以综合起来得到块体大小分布和孔隙分布。

岩石可以是页岩、致密砂岩或者其它非常规储层岩样。

步骤200，将组分特征根据矿物颗粒的类型不同而细分为a1、a2……an个对应的子组分特征；将孔隙特征对应所述子组分特征的不同而细分为b1、b2……bn个对应的子孔隙特征；将辅助特征根据与不同所述子组分特征的关系而细分为c1、c2……cn个对应的子辅助特征；利用混合分形方法分别拟合每种子组分特征的分形参数；

在划分时，将岩样图像中的矿物，如，有机质块体、黄铁矿块体、石英，长石，蒙脱石等一类矿石类型统一划分为组分特征，并可对该组分特征中的每个具体矿物按子组分特征a1、子组分特征a2、子组分特征a3……an的编号进行细分类；将图像中的粒间孔、粒内孔等孔隙全部划分为孔隙特征中，其中根据孔隙的分布位置如有机质块体中的有机孔，黄铁矿块体中的孔，基质中的无机孔等，按子孔隙特征b1、子孔隙特征b2、子孔隙特征b3……bn进行细分类；辅助特征是除前述两种特征外的其它特征，如裂纹特征、密度特征，本实施方式是以裂缝作为辅助特征，裂缝根据类型可分为自然裂缝和人工裂缝，其根据裂缝密度、缝长、缝宽等信息将辅助特征分为子辅助特征c1、子辅助特征c2、子辅助特征c3、……cn。

步骤300，建立立方体构成的表征单元体，将表征单元体整体作为子组分特征a1，再将a2、a3按照分形重构的块体大小和个数顺序放进表征单元体，a2、a3的块体所占的几何空间的组分信息替换掉a1的组分信息，再将子孔隙特征b1随机放进表征单元体，统计每个孔的壁面的最大组成成分，作为该孔隙的组成成分；

立方体的边长可以为0.4mm，对于组分特征，如果子组分特征a1、a2……an的块体大小分布符合分形特征，则先提取a1的块体大小分布频谱，使用混合分形方法对a1的块体大小分布频谱进行拟合，得到相应a1的混合分形参数af1，使用混合分形参数af1重构出a1的块体大小分布频谱得到等同于与a1的块体大小分布频谱；依次对a2、a3做同样处理。

对于组分特征，如果子组分特征a1、a2……an的块体大小分布不符合分形特征，则直接统计a1、a2、……an的面积或体积占比，即ai的占比为

其中为ai的面积或体积占比，vai为ai的面积或体积，j为整数。

子组分特征a2、a3采用随机的方式放进表征单元体空间内，且其中任意块体占用空间互相不重叠。

对于孔隙特征，提取子孔隙特征b1的孔隙等效直径分布频谱，使用混合分形方法对b1的孔隙等效直径分布频谱进行拟合，得到b1的混合分形参数bf1，使用混合分形参数bf1重构出b1的孔隙等效直径分布频谱，得到等同于b1的孔隙等效直径分布频谱；依次对b2、b3做同样处理。

孔隙特征b1同样采用随机的方式放进表征单元体空间内，且其中子孔隙特征b1所占的空间可以和任意块体占用的空间相互重合；统计子孔隙特征b1所占的空间中各个子组成特征a1、a2……an的占比，以赋予子孔隙特征b1最大占比的组成成分所对应的矿物性质。

另一对b1赋予壁面性质的方法是：对子孔隙特征b1赋予壁面性质，过程如下：将子孔隙特征b1采用随机的方式放进表征单元体空间内，统计所述b1所占的空间中各个子组成特征a1、a2……an的面积和体积所占比例，等比例赋予所述b1壁面性质；壁面性质包括壁面润湿性和壁面粗糙度。

其中的分形参数包括迭代次数i，分形标度f，第i次迭代产生的孔隙/块体边长λi，最大孔隙/块体边长λmax，第i次迭代中产生的孔隙/块体个数ni，第i次迭代中不参与迭代的孔隙/块体个数ni-soild。

在构造连通概率矩阵时，将子孔隙特征b1和子辅助特征c1综合，按照等效直径从大到小排序，等效直径相同时，按对应毛管力从大到小排序，构造完成的连通概率矩阵中，其第一列为等效直径，第二列为等效直径所对应的个数，第三列连通概率，第四列为该行驱扫概率，第五列为该行导通概率，第六列为该行壁面性质。

即每级孔对应的等效孔径为λ1,λ2,…,λi,…,λm，每级孔对应的个数分别为n1,n2,…,ni,…,nm；在排序后的连通概率矩阵中使用连通孔隙概率算法得到第i级连通概率p1,i；模型为非润湿相驱替润湿相，在第i级的ni各孔中，驱扫概率为p2,i，使每级孔中驱扫成功的孔隙个数为p2,i×ni，驱扫概率p2,i通过真实岩心压汞曲线和孔隙成像结果综合分析给出；第i级孔导通概率为p3,i，当驱替动力大于不等径孔隙连接处阻力时，孔隙导通(p3,i＝1)，否则孔隙不能导通(p3,i＝0)，影响p3,i的参数有毛管力，润湿性，驱替压力、孔隙等效直径，界面张力等参数。

步骤400，设定模型为油驱水，油水界面张力30mn/m，油粘度1.39mpa.s，水粘度1.05mpa.s，岩心表观渗透率1000md，孔隙度22％，接触角0°，建立连通概率矩阵，通过计算即可得到该模型的相渗曲线图。

计算相渗曲线的公式为：

krw(j)为润湿相相对渗透率，kro(j)为非润湿相相对渗透率，srw(j)为润湿相饱和度，ci为第i级孔隙形状因子；其中1≤j≤m且j为整数，其它的参数参见前述限定。

如图2所示，图中左侧则低至高的曲线为水相相渗曲线，而由左侧向右侧月逐渐降低的曲线为油相相渗曲线，其中的点为实验值，可以验证本方法计算相渗曲线结果有效。

本实施方式是基于真实岩石的成像技术，利用岩石孔隙空间的分形特性，使用混合分形单元法提取分形特征，将孔隙壁面赋予不同润湿参数，使用连通概率，驱扫概率，导通概率共同构建基于三维混合分形的两相流模型。该模型能够有效表征不同尺度下的非常规岩石样品孔隙结构特征，能够快速计算非常规岩石样样品相渗曲线。与传统方法相比，该方法涵盖信息多(孔隙分布、润湿性)、计算速度快，便于工程应用，为获得非常规储层相对渗透率开辟了一种新方法。

下表为相渗方法下一种优选的连通矩阵示意：

其中块体λb(a1)-1代表λb孔隙赋予子组分特征a1性质之后的第一个孔隙的等效直径，λb(a1)-2代表λb孔隙赋予子组分特征a1性质之后的第2个孔隙的等效直径，块体λb(a2)-1代表λb孔隙赋予子组分特征a2性质之后的第一个孔隙的的等效直径。

孔隙等效直径满足第一列从大到小。直径相等时使用毛管力比大小，非润湿相驱替润湿相时毛管力从大到小排，润湿相驱替非润湿相时毛管力从小到大排。

至此，本领域技术人员应认识到，虽然本文已详尽示出和描述了本发明的多个示例性实施例，但是，在不脱离本发明精神和范围的情况下，仍可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其他变型或修改。因此，本发明的范围应被理解和认定为覆盖了所有这些其他变型或修改。