[发明专利]一种基于孔隙结构的储层临界物性获取方法

摘要

本发明公开了一种基于孔隙结构的储层临界物性获取方法，包括以下方法步骤：步骤1：以压汞实验样品测得的孔隙度Φ、渗透率K、排驱压力P_d为基础；步骤2：计算每一类储层实验样品的平均毛管压力曲线；步骤3：得到J函数，J＝CP；步骤4：重复步骤2、3，得到II、III、IV三类储层的平均毛管压力曲线；步骤5：分别求取四类储层的累积渗透率贡献步骤6：绘制每一类储层孔喉半径随累积渗透率贡献变化关系；步骤7：求取临界渗透率；步骤8：回归得到每一类储层渗透率回归模型：K＝dΦ^e；步骤9：根据每一类储层渗透率回归公式。该方法考虑了不同孔隙结构储层临界孔隙度的差异，显著提高矿井储层尤其是非均质性较强储层的临界物性获取精度，实用性强。

主权项

1.一种基于孔隙结构的储层临界物性获取方法，其特征在于：包括以下方法步骤：步骤1：以压汞实验样品测得的孔隙度Φ、渗透率K、排驱压力Pd为基础，采用K均值聚类方法对实验样品进行聚类分析，将区域储层划分为若干类型，后续步骤以I、II、III、IV共4类储层为例；步骤2：计算每一类储层实验样品的平均毛管压力曲线，以I类储层为例，首先计算I类储层每一块压汞样品常熟C：其中，σ为汞与空气的界面张力，其数值为480mN/m；θ为汞与固体湿相接触角，数值为140°；Φ为样品孔隙度，％；K为样品渗透率，10‑3×μm2；步骤3：将I类每块样品的常数C乘以该样品各测点的压力(P)得到J函数，J＝CP；采用算术平均方法计算得到平均J值与平均汞饱和度曲线，再将平均J函数曲线上各测点除以所有样品的C平均值，得到反映A类储层特征的平均毛管压力曲线；步骤4：重复步骤2、3，得到II、III、IV三类储层的平均毛管压力曲线；步骤5：分别求取四类储层的累积渗透率贡献其中，为累积渗透率贡献；ΔK_i为区间渗透率；ΔS_i～i+1为区间进汞饱和度；P_ci为i测点进汞压力，MPa；步骤6：绘制每一类储层孔喉半径随累积渗透率贡献变化关系，当累积渗透率贡献达99.9％时，对应的孔喉半径即为储层的最小流动孔喉半径Rmin；步骤7：求取临界渗透率，根据区域毛管压力实验样品参数，回归得到储层渗透率与孔喉半径平均值间的关系：将R_min代入即可获得临界渗透率值步骤8：以区域岩心物性分析数据为基础，根据储层类型划分标准将岩心样品分类，回归得到每一类储层渗透率回归模型：K＝dΦe；步骤9：根据每一类储层渗透率回归公式，将K_临代入公式得到每一类储层的临界孔隙度值，