[发明专利]毛管压力曲线分类方法

说明书

本发明公开了一种毛管压力曲线分类方法，首先对毛管压力曲线进行拟合，确定毛管压力曲线所反映的孔隙系统数量N和排驱压力P_d，并根据P_d的大小对曲线进行排序；然后比较任意两条曲线的相似性，确保两曲线具有相同的N值，并逐个压力点计算两条毛管压力曲线所对应的饱和度差值的绝对值，如果所有点的饱和度差值都小于容许误差，则将两条曲线划分为同一类别。本发明确保了曲线整体形态的相似性，并可根据需要调整分类的精度，具有较大的灵活性。

技术领域

本发明涉及一种毛管压力曲线分类方法，用于孔隙介质的毛细管压力曲线定量分析，属于储集岩孔隙结构评价领域。

背景技术

对于储集岩孔隙结构研究而言，毛管压力实验是研究喉道大小的基础，不同的毛管压力曲线形态特征反映了岩石不同的孔隙结构特征。在储层评价中，常常根据毛管压力曲线的形态特征对储层进行分类，从而提高储层评价的准确性。因此，对毛管压力曲线进行准确的分类是储层分类评价的基础。

传统的毛管压力曲线分类方法几乎都是利用P_d（排驱压力）、R35（汞饱和度35%对应的喉道半径）、R50（汞饱和度50%对应的喉道半径）和S_wir（束缚水饱和度）等曲线特征参数中的一个或多个，通过聚类的方法对其进行分类。这些分类方法存在明显的不足之处，常常将形态差异较大的曲线划分为同一类，尤其在非均质性较强的碳酸盐岩中的应用效果较差。如图1（a）中的两条毛管压力曲线明显具有不同的形态特征，但两条曲线的R35和R50值都非常接近，如果仅根据R35或R50对其进行分类，则会将两者划分为同一类；同样在图1（b）中，仅根据P_d和S_wir也会将两条差异较大的曲线分为同一类。即使利用P_d、R35、R50和S_wir等多个参数进行综合分类，也可能会将差异较大的两条曲线分为同一类，如图1（c）所示。

由此可见，现有技术对毛管压力曲线的分类仅依据特征参数进行分类，没有考虑曲线整体形态的相似性，除低了毛管压力曲线分类的准确性、快速性和灵活性。

发明内容

本发明的目的在于克服现有技术存在的上述问题，提供一种毛管压力曲线分类方法。本发明确保了曲线整体形态的相似性，并可根据需要调整分类的精度，具有较大的灵活性。

为实现上述目的，本发明采用的技术方案如下：

一种毛管压力曲线分类方法，其特征在于：首先对毛管压力曲线进行拟合，确定毛管压力曲线所反映的孔隙系统数量N和排驱压力P_d，并根据P_d的大小对曲线进行排序；然后比较任意两条曲线的相似性，确保两曲线具有相同的N值，并逐个压力点计算两条毛管压力曲线所对应的饱和度差值的绝对值，如果所有点的饱和度差值都小于容许误差，则将两条曲线划分为同一类别。

所述方法中，利用（1）式对所有原始毛管压力曲线逐一进行拟合，确定每一条曲线所反映的孔隙系统数量N及排驱压力P_d；

（1）

式中：

P_c—注汞压力（bar或Mpa）；

S_w—含水饱和度（小数）；

N—孔隙系统的数量（无量纲）；

a、b、c、d—需要拟合的系数（无量纲）；

tanh—双曲正切函数。

所述方法中，按排驱压力P_d从高到低的原则，对毛管压力曲线进行排序，依次编号，并将所有曲线的类别属性初始化为0。