[发明专利]储层不可及孔隙体积确定方法及装置

权利要求书

1.一种储层不可及孔隙体积确定方法，其特征在于，包括：

获得测试储层的所有孔隙孔径数据，将所有孔隙孔径数据按照从小到大顺序进行排列；

获得聚合物的所有分子团粒径数据，将所有分子团粒径数据按照从小到大顺序进行排列；

对排序后的所有孔隙孔径数据和排序后的所有分子团粒径数据进行分析：

当分析得到排序后的所有孔隙孔径数据和排序后的所有分子团粒径数据均是离散的，采用如下方式计算不可及孔隙体积IPV：

计算每一个孔隙孔径数据对应的IPV值概率分布律；

根据每一个孔隙孔径对应的分子团粒径临界值以及所有分子团粒径数据，计算每一个孔隙孔径数据点对应的IPV值；

根据每一个孔隙孔径数据对应的IPV值概率分布律、每一个孔隙孔径数据点对应的IPV值，计算测试储层的不可及孔隙体积IPV；

采用如下公式计算不可及孔隙体积IPV：

其中，N为孔隙孔径数据的个数，i＝1,2,…,N，a＝1,2,…,N；x_i为第i个孔隙孔径数据；x_a为第a个孔隙孔径数据；φ_i为第i个孔隙孔径数据对应的IPV值概率分布律；M为分子团粒径数据的个数；P_i为第i个孔隙孔径数据点对应的IPV值；K_i为第i个孔径数据点对应的粒径临界值在所有粒径数据点中的序号；IPV为测试储层的不可及孔隙体积IPV；

当分析得到孔隙孔径的区间和分布频率以及分子团粒径的区间和分布频率时，采用如下方式计算不可及孔隙体积IPV：

根据每一个孔隙孔径区间的分布频率，计算每一个孔隙孔径区间对应的IPV值概率分布律；

根据每一个孔隙孔径区间对应的分子团粒径区间临界值以及所有分子团粒径区间的分布频率，计算每一个孔隙孔径区间对应的IPV值；

根据每一个孔隙孔径区间对应的IPV值概率分布律、每一个孔隙孔径区间对应的IPV值，计算测试储层的不可及孔隙体积IPV；

采用如下公式计算不可及孔隙体积IPV：

其中，N为孔隙孔径区间的个数，i＝1,2,…,N，a＝1,2,…,N；x_i为第i个孔隙孔径区间的均值数据；x_a为第a个孔隙孔径数据；φ_i为第i个孔隙孔径区间对应的IPV值概率分布律；f(x_i)为第i个孔隙孔径区间的分布频率；M为分子团粒径区间的个数；P_i为第i个孔隙孔径区间对应的IPV值；K_i为第i个孔径区间对应的粒径临界值在所有粒径区间中的序号；y_j为第j个分子团粒径区间的均值数据；g(y_j)为第j个分子团粒径区间的分布频率；IPV为测试储层的不可及孔隙体积IPV；

当分析得到孔隙孔径的分布概率密度函数以及分子团粒径的分布概率密度函数时，采用如下方式计算不可及孔隙体积IPV：

根据孔隙孔径的分布概率密度函数，计算任一孔隙孔径对应的IPV值概率密度函数；

根据任一孔隙孔径对应的分子团粒径临界值以及分子团粒径的分布概率密度函数，计算任一孔隙孔径对应的IPV值；

根据任一孔隙孔径对应的IPV值概率密度函数、任一孔隙孔径对应的IPV值，计算测试储层的不可及孔隙体积IPV；

采用如下公式计算不可及孔隙体积IPV：

其中，x为孔隙孔径数据，为连续型随机变量；φ(x)为任一孔隙孔径IPV值概率密度函数；[x_min，x_max]为孔隙孔径的范围；F(x)为孔隙孔径在范围[x_min，x_max]内的概率密度函数；P(x)为孔径为x的孔隙的IPV值；Ω(x)为孔径为x的孔隙对应的分子团粒径临界值；G(y)为分子团粒径在范围[y_min，y_max]内的概率密度函数；y为分子团粒径数据，为连续型随机变量；IPV为测试储层的不可及孔隙体积IPV。