基于nmr测井和毛管压力资料的储层渗透率研究.docxVIP

基于nmr测井和毛管压力资料的储层渗透率研究 目前，许多方法和储层渗透法已经开发出。然而, 这些模型都是一些经验的关系式, 不能进行大规模的推广应用。自从核磁共振测井技术诞生以来, 许多学者相继提出了利用核磁共振测井资料计算储层渗透率的经典SDR模型和Tim-Coates模型。通过对这两种模型的分析不难发现, 要利用核磁共振测井资料准确计算储层渗透率必须准确地确定储层孔隙度, 孔隙可动流体体积和束缚流体体积。在实际储层评价中, 孔隙可动流体体积和束缚流体体积的确定主要依赖于核磁共振测井T2截止值的求取, 然而到目前为止, 还没有一种方法能准确计算出T2截止值。往往将T2截止值取单一值, 在砂泥岩中T2截止值取值为33 ms, 在碳酸盐岩中取值为92 ms。油田实际分析表明,T2截止值取单一值在中国大多数油田是不成立的, 而且在这两种模型中需要确定的参数较多, 不适合在储层评价中普遍应用。为了解决现有模型存在的问题, 笔者提出了结合核磁共振测井和毛管压力资料准确计算储层渗透率的新方法, 并将其应用于低孔隙度、低渗透率砂岩中以连续地计算储层渗透率。 1 毛管压力曲线 一般情况下, 毛管压力曲线作图采用半对数坐标系, 进汞饱和度为线性横坐标, 进汞压力为对数纵坐标。如果将毛管压力曲线作于双对数坐标系下时, 其形态近似于双曲线, 毛管压力曲线形态特征如图1所示。 该毛管压力双曲线函数的数学表达式为 lg(pcpd)×lg(SΗgSΗg∞)=-C(1)lg(pcpd)×lg(SHgSHg∞)=?C(1) 式中:pc为进汞压力, MPa;pd为驱替压力, MPa;SHg为进汞饱和度, %;SHg∞为无限大毛管压力下挤入的非润湿相汞饱和度, %;C是孔隙几何因子。 通过pd和SHg∞作两条分别垂直于纵坐标和横坐标的直线即为该双曲线的两条渐近线。Swanson等人通过对大量的高孔隙度、高渗透率的岩心数据分析发现, 通过不同岩心样品实验观察分析得到的有效控制流体流动的主孔隙系统时的汞饱和度值与双对数坐标系下的毛管压力曲线拐点是相对应的。在拐点出现前, 非润湿相占据了有效的相互连通的孔隙空间;过了拐点, 非润湿相开始进入更细小的孔隙空间或占据不规则孔隙的一些边角部分, 非润湿相在其中的流动能力明显下降。因此, 拐点处的汞饱和度指对流体流动有效贡献的那部分孔隙空间的体积, 而对应的毛管压力是指连通整个有效孔隙空间的最小喉道大小。毛管压力曲线的拐点在双对数坐标系下就是双曲线的顶点 (见图1中A点) 。在A点处, 进汞饱和度SHg与毛管压力pc的比值SΗgpcSHgpcA(Swanson参数) 比毛管压力曲线上其他部分都要高, 这表示在该点处单位压力内进汞量最多。 2 sanson参数与储层参数的关系 基于上述分析, 通过对实际岩心压汞毛管压力资料进行处理, 笔者发现Swanson参数与储层孔隙度、渗透率等参数之间具有良好的相关性。根据此相关性, 可以建立相应的模型以利用毛管压力资料求取岩石的渗透率。为了建立广泛适用的模型, 选取渗透率范围介于 (0.002～1 150) ×10-3μm2岩心样品, 通过对其形态进行分析, 提取Swanson参数, 并将其与孔隙度、渗透率进行交会, 建立了Swanson参数与岩石孔隙度、渗透率等参数之间的相关关系, 其对应关系如下: √Κ?=-8×10-8×(SΗgpc)2A+0.0011×(SΗgpc)A+0.1394(2)K???√=?8×10?8×(SHgpc)2A+0.0011×(SHgpc)A+0.1394(2) 式中:√Κ?K???√为岩石综合物性指数;SΗgpcA为Swanson参数。 从式 (2) 可以看到, Swanson参数与岩石综合物性指数之间具有很好的二次函数关系, 二者相关系数达到0.997 1, 证明利用该模型可以用Swanson参数来准确地计算岩石综合物性指数, 进而可以计算出渗透率, 并将其推广到实际储层条件下, 用以准确地计算出储层的渗透率。 3 核磁共振tm与sanson参数的相关性 由式 (2) 可知, 要利用毛管压力资料来计算储层渗透率, 须利用压汞毛管压力资料来准确地求取Swanson参数。由于受到压汞实验费用昂贵以及实验周期长等因素的影响, 实际油田的压汞实验分析资料较少, 只会在可能的目的层段进行有限的取心压汞实验分析。为了实现利用压汞毛管压力资料在整个井段连续地计算储层渗透率, 笔者提出了利用核磁共振测井横向弛豫时间T2几何平均值 (T2lm) 来获取Swanson参数的方法。 选取36块岩心样品进行核磁共振测井实验和常规物性分析。通过对测量的核磁共振回波数据进行解谱、计算和整理, 得到了核磁共振T2lm、束缚水饱和度、岩石孔隙度和渗透率等信息。核磁共振T2lm是核磁共