[发明专利]一种页岩储层孔隙分类方法、系统及终端设备

权利要求书

1.一种页岩储层孔隙分类方法，其特征在于，包括：

获取页岩样品的孔径分布曲线；

基于多重分形理论分析所述孔径分布曲线，确定多重分形奇异谱宽和广义分形维数；

根据所述广义分形维数计算赫斯特指数；

组合所述多重分形奇异谱宽随孔径的变化曲线和所述赫斯特指数随孔径的变化曲线；

基于组合后的变化曲线和阈值对所述页岩样品的孔隙进行分类。

2.如权利要求1所述的页岩储层孔隙分类方法，其特征在于，所述基于多重分形理论分析所述孔径分布曲线，得到多重分形奇异谱宽和广义分形维数，包括：

将所述孔径分布曲线转化为累计孔体积-孔径分布的新孔径分布曲线；

将所述新孔径分布曲线的孔径分布范围L划分为2^j个孔径范围为L₁的第一区间，j为正整数且j≥1；

将每个第一区间划分为2^k个尺寸为ε的第一子区间，k为正整数且k≥0；

计算第i个第一子区间的孔体积概率密度函数为：P_i＝V_i/V_t，P_i为概率密度函数，V_i为第i个第一子区间的孔体积，V_t为总孔体积；

根据所述孔体积概率密度函数计算多重分形奇异谱(α～f)，表达式如下；

式中，N为子区间的总数，q为统计矩阶数，–∞q+∞，q取最小值和最大值时的α的差值为所述多重分形奇异谱宽Δα；

根据所述多重分形奇异谱确定所述多重分形奇异谱宽变化曲线；

根据所述孔体积概率密度函数和所述第一子区间的尺寸计算所述广义分形维数D_q，表达式如下：

所述根据所述广义分形维数计算赫斯特指数，包括：

根据H＝(D_q+1)/2计算得出所述赫斯特指数，H为所述赫斯特指数。

3.如权利要求2所述的页岩储层孔隙分类方法，其特征在于，所述方法还包括：获取多重分形奇异谱宽随孔径的变化曲线；

所述获取多重分形奇异谱宽随孔径的变化曲线，包括：

根据每个第一子区间的所述多重分形奇异谱宽Δα和所述赫斯特指数H，确定第一变化曲线，所述第一变化曲线表征所述孔径分布范围L与孔径分布的关系；

将所述孔径分布曲线的孔径分布范围L'划分为2^j-1个孔径范围为L₁'的第二区间，j为正整数且j≥1；

将每个第二区间划分为2^k'个尺寸为ε'的第二子区间，k'为正整数且k'≥0；

根据每个第二子区间的多重分形奇异谱宽Δα'和赫斯特指数H'，确定第二变化曲线，所述第二变化曲线表征所述孔径分布范围L'与孔径分布的关系；

根据所述第一变化曲线和所述第二变化曲线得到所述多重分形奇异谱宽随孔径的变化曲线。

4.如权利要求1所述的页岩储层孔隙分类方法，其特征在于，所述获取页岩样品的孔径分布曲线，包括：

获取页岩样品的第一孔径分布曲线和第二孔径分布曲线，所述第一孔径分布曲线基于所述页岩样品的第一气体吸附实验确定，所述第二孔径分布曲线基于所述页岩样品的第二气体吸附实验确定，所述第一气体吸附实验与所述第二气体吸附实验不同；

基于所述第一孔径分布曲线和所述第二孔径分布曲线，生成所述孔径分布曲线。

5.如权利要求4所述的页岩储层孔隙分类方法，其特征在于，所述第一孔径分布曲线基于所述页岩样品的二氧化碳吸附实验数据得到，所述第二孔径分布曲线基于所述页岩样品的氮气吸附实验数据得到。

6.如权利要求5所述的页岩储层孔隙分类方法，其特征在于，所述基于所述第一孔径分布曲线和所述第二孔径分布曲线，生成所述孔径分布曲线，包括：

截取所述第一孔径分布曲线小于等于第一孔径的部分曲线，和所述第二孔径分布曲线大于等于第一孔径的部分曲线；其中，所述第一孔径为所述第一孔径分布曲线和所述第二孔径分布曲线相交处的孔径大小；

将两部分曲线合成所述孔径分布曲线。