[发明专利]基于多样品整体虚拟测量的多毛管力曲线平均化方法

权利要求书

1.基于多样品整体虚拟测量的多毛管力曲线平均化方法，依次包括以下步骤：

第一步：取得m种岩石样品，得到每种样品的毛管力-饱和度曲线p_c(S_w)-S_w、外观体积V_b和孔隙度φ，将样品j的毛管力-饱和度曲线记为外观体积记为孔隙度记为φ^j；

第二步：检查每种样品的毛管力-饱和度曲线质量并对曲线端点进行预处理，过程如下：

(1)确保毛管力-饱和度曲线中毛管力随饱和度S_w的变化满足单调函数特征；

(2)确保毛管力-饱和度曲线的最大毛管力值大于平均化毛管力-饱和度曲线所要求的最大毛管力值p_cmax；

(3)将毛管力-饱和度曲线中平缓段外推延长至湿相饱和度S_w＝1处，所对应的毛管力即为排驱压力，取所有样品中最小排驱压力为p_cmin，并将外推延长得到的曲线段代替相同饱和度区间原来的曲线段；

(4)将经过上述处理后的毛管力-饱和度曲线命名为

第三步：在实数区间[p_cmin,p_cmax]按等间距或不等间距取出一系列足够多的数值，记为{p_ci}，0≤i≤N，其中下标i表示第i个数据点，且p_c0＝p_cmin，p_cN＝p_cmax，计算多个样品在整体虚拟测量下，不同毛管力值p_ci所对应的平均湿相饱和度过程如下：

(1)当毛管力-饱和度曲线采用解析函数表示：

由于毛管力随饱和度的变化满足单调性，毛管力-饱和度函数有反函数：

采用下式计算多样品的平均湿相饱和度：

(2)当毛管力-饱和度曲线采用图形曲线表示，采用图形查找方式，将序列{p_ci|p_ci∈[p_cmin,p_cmax],p_c0＝p_cmin,p_cN＝p_cmax,i∈{0,1,..,N}}中每一毛管力值p_ci在样品j(j∈{1,2,...,N})的毛管力-饱和度曲线上找到对应的湿相饱和度根据下式计算多样品的平均湿相饱和度

(3)当毛管力-饱和度曲线采用数据表描述，根据数据表所确定的函数关系表示样品j毛管力数据表中数据点的总数，采用内部插值方式确定p_ci在样品j的毛管力-饱和度曲线上的值根据下式计算多样品的平均湿相饱和度

第四步：以平均湿相饱和度为横坐标，以毛管力p_ci为纵坐标，将数据点标在图形上，光滑地连接所有的数据点，即得到平均化的毛管力曲线。

2.如权利要求1所述的基于多样品整体虚拟测量的多毛管力曲线平均化方法，其特征在于，所述第二步中，确保所有毛管力-饱和度曲线的最大值大于平均化毛管力-饱和度曲线所要求的最大毛管力值p_cmax，是指将不满足要求的毛管力-饱和度曲线的端点曲线采用外推法延长，直至毛管力-饱和度曲线的最大毛管力值超过p_cmax。

3.如权利要求1所述的基于多样品整体虚拟测量的多毛管力曲线平均化方法，其特征在于，所述第三步中，当毛管力-饱和度曲线采用图形曲线表示或当毛管力-饱和度曲线采用数据表描述，如果p_ci小于样品j的毛管力-饱和度曲线上的最小毛管力值，取

4.如权利要求1所述的基于多样品整体虚拟测量的多毛管力曲线平均化方法，其特征在于，所述第三步中，当毛管力-饱和度曲线采用数据表描述，根据数据点作出毛管力-饱和度曲线图形，将其转换为第(2)情况处理。

5.如权利要求1所述的基于多样品整体虚拟测量的多毛管力曲线平均化方法，其特征在于，所述第三步中，当毛管力-饱和度曲线采用函数、图形和数据表混合方式表示，按本步骤(1)～(3)种情况分别得到平均湿相饱和度