油润湿致密砂岩核磁共振弛豫机理与流体识别方法

2020-06-23张海涛方驭洋李高仁谭茂金

石油物探 2020年3期

张海涛，方驭洋，李高仁，谭茂金，王谦

（1.中国石油天然气股份有限公司长庆油田分公司勘探开发研究院，陕西西安730001，2.中国地质大学（北京）地球物理与信息技术学院，北京100083）

鄂尔多斯盆地致密砂岩储层物性差，孔隙结构和油水关系复杂，流体性质判识难，测井解释难度大，符合率低[1-2]。核磁共振（NMR）测井在储层孔隙度、渗透率、含油性及孔隙结构评价中具有明显优势[3-6]。然而，常用的核磁共振测井流体识别方法大多是针对亲水地层，在低孔低渗、油润湿的储层中，流体识别效果差，精度不高[7-9]。

近年来，国内外专家学者在低孔低渗砂岩储层核磁共振流体识别方法研究方面进行了大量探索和尝试。肖立志[10]从理论上剖析了不同湿润性、不同黏度油的岩石核磁共振响应特征，并给出了相应的示意图。谢然红等[11]通过实验研究了原油的核磁共振弛豫特性随原油粘度、NMR 仪器测量频率以及测量回波间隔（TE）等的不同而发生的变化，对利用核磁共振测井资料准确识别和评价油水层具有重要的指导意义，以实验为基础提出了核磁共振测井流体识别中需要关注的影响因素。目前，对油润湿条件下低孔低渗砂岩储层的核磁共振弛豫机理的研究较少，导致多数方法在油润湿储层中应用效果较差，流体识别精度不高。谭茂金等[12]在双等待时间（TW）的观测模式下，利用可动流体体积的差构建了流体指示参数，在低对比度油层判识中取得了较好的应用效果，此外，其研究表明利用遗传算法全局搜索油气的纵向弛豫时间（T1）和横向弛豫时间（T2）值，再用阻尼最小二乘法（LSQR）反演标准T2分布和差谱，可以有效解决一维核磁共振双TW观测数据存在的非线性反演简化、信噪比低的问题，实现油气水的准确识别和定量评价[13]。吕婕[14]利用核磁共振测井资料，结合交会图版、伪毛管压力曲线等方法进行孔隙结构分析和流体识别，提高了解释符合率。该方法的关键在于结合储层特点和核磁共振弛豫机理寻找对油水敏感的测井数据以及相关储层参数，通过优选和构建能够判识流体类型的识别因子来提高流体识别的精度。胡法龙等[15]提出用核磁共振测井数据构建水谱法进行流体识别，该方法利用长等待时间、短回波间隔模式下测量的T2分布构建完全含水状态下长等待时间、长回波间隔模式下的T2分布，通过对比测量T2分布和构建完全含水T2分布之间的差异确定储集层的流体类型。

图1 岩石孔隙系统中的润湿性[20]a水润湿（水与岩石颗粒直接接触）；b油润湿（油与岩石颗粒直接接触）；c混合润湿（岩石颗粒部分与水接触，部分与油接触）

本文基于油湿、水湿以及混合润湿条件下的核磁共振驰豫机理，利用数值模拟方法分析了不同润湿条件下的核磁共振响应特征。选取对流体敏感的核磁共振差谱几何均值与有效孔隙度差构建交会图版，依据流体识别图版建立流体类型测井判别方法，并在鄂尔多斯延长组长8储层中进行了检验。

1 储层润湿性特征

出现在自由弛豫位置，位于T2分布后部，T2分布呈双峰特征。这说明，在水湿条件下，水受表面弛豫和自由弛豫影响，弛豫速率变快，油只受自由弛豫影响，弛豫速率慢。

在油湿条件下，分别设置含水饱和度Sw＝95%，80%，40%，10%，当含水饱和度为95%时指示为水层。模拟结果如图3b所示，油峰出现在T2谱前部，水峰则出现在自由弛豫位置，位于T2谱后部，故T2分布呈双峰特征。这说明，在油湿条件下，油受表面弛豫影响，弛豫速率变快，水只受自由弛豫的影响，弛豫速率慢。

在混合润湿条件下，分别设置含水饱和度Sw＝95%，80%，40%，10%，当含水饱和度为95%时指示为水层。同时，设置水湿比、油湿比均为0.5。模拟结果如图3c所示，T2呈单峰分布。这说明，在混合润湿条件下，油受自由弛豫和表面弛豫共同影响，但自由弛豫对其弛豫时间的影响要远小于表面弛豫，油峰向T2分布的前部移动，而水基本上仅受自由弛豫影响，仍在其自由弛豫位置，因此，整个T2分布趋于单峰特征。

图3 不同润湿条件下砂岩模型核磁共振数值模拟与T 2 分布特征a水湿；b 油湿；c混合润湿

可以看出，在油湿条件下，随着含水饱和度的减小，核磁共振T2分布逐渐往右移动，当含水饱和度为10%时，核磁共振T2分布位于10～200 ms，指示了油层的核磁弛豫特征。根据核磁共振双TW观测方式的工作原理与弛豫机理，由长短等待时间回波串计算的差谱也在这一范围内，利用差谱的T2几何均值可以指示含油性。根据数值模拟结果，地层含油时，差谱几何均值应为10～200 ms，而且，长等待时间核磁共振孔隙度要大于短等待时间核磁共振孔隙度，两者差值越大说明含油饱和度越大，且储层孔隙度也越好。因此，考虑利用核磁共振双TW观测模式下得到的差谱与有效孔隙度差进行流体识别。

2.2 流体识别图版

根据水湿与油湿条件下核磁数值模拟结果，建立了核磁共振双TW观测模式下差谱几何均值与有效孔隙度差的交会图。首先，对目标区域井的核磁共振双TW测井进行数据处理，获取长短不同等待时间TW的回波和有效孔隙度，构建目标井储层核磁共振双TW观测差谱几何均值-有效孔隙度差识别图版，x轴为差谱几何均值T2LM，y轴为有效孔隙度差Δφ，如图4所示，可以看出，该图版中油层和水层的界限明显，油层的差谱几何均值在80～250 ms，有效孔隙度差大于0.5%；水层的差谱几何均值小于80 ms或大于250 ms，有效孔隙度差小于0.5%。

可以看出，核磁共振测井双TW观测模式下差谱的形态及其几何均值反映流体的响应差异，长短等待时间的有效孔隙度差值大小能够反映流体弛豫的差异，而且采用有效孔隙度差值而不是可动流体的差值避开了T2截止值选取的影响，更能够体现流体性质的差异。

2.3 识别方法

图4 核磁共振双T W 观测差谱几何均值-有效孔隙度差识别图版

在利用图4所示的图版进行核磁共振测井解释时，拟采用构建流体指示曲线来指示。首先，对目标井的核磁共振测井数据进行处理，计算差谱几何均值和有效孔隙度差；然后，根据计算得到的差谱几何均值和有效孔隙度差的数值，与流体识别图版中的流体边界值进行比较，确定储层的流体类型。当差谱几何均值和有效孔隙度差Δφ均在图版的油层区域内时，则用数字“1”表示；差谱几何均值或有效孔隙度差有一个不在油层区域内时，则用数字“0”表示；根据数字“0”或“1”绘制流体指示曲线，综合确定储层的流体类型。当储层流体类型指示曲线为连续形态时，判断目标储层为油层；而当指示曲线为单一“脉冲”或无“脉冲”时，判断目标储层为水层。

3 实例与效果分析

鄂尔多斯盆地是我国第二大沉积盆地，中生界延长组富含油气资源。陇东地区位于鄂尔多斯盆地西南部，区域上横跨伊陕斜坡、天环凹陷和西缘冲断带，地质构造复杂。延长组长8段是该区主要的含油层系，储层主要以长石砂岩为主，孔隙度小，主要集中在4%～12%，平均孔隙度为8.72%，基质渗透率低，主要分布在0.01～0.30 m D（1 mD≈0.987×10－3μm2），平均渗透率约为0.95 m D。长8 段储层孔隙类型多样，非均质性强，油水分异较差，而且受润湿性影响，流体性质识别困难，测井解释符合率低。

研究区长8段储层岩心吸油吸水实验表明，岩心存在吸油现象，且吸油量明显大于吸水量（表1），吸油量约占孔隙度的70%，吸水量约占孔隙度的30%。油主要以薄膜状吸附在岩石表面，水则分布在孔隙中央，形成典型的“油包水”特征。图5为长8储层岩心吸油7天后形成的“油包水”分布，反映了长8储层的弱油湿特征[20]。

表1 岩心润湿性实验分析结果

图5 长8储层岩心吸油7天后形成“油包水”分布[24]

针对油润湿的低孔低渗砂岩储层，利用核磁共振测井双TW观测模式下采集的数据，进行差谱分析，计算差谱的T2几何均值T2LM；通过T2反演，分别得到长等待时间下和短等待时间下的有效孔隙度，计算两者的差值Δφ。利用差谱几何均值和有效孔隙度差构建交会图，形成T2LM-Δφ流体识别图版。

图6 Luo99井核磁共振测井解释与流体识别结果

Luo99井为鄂尔多斯盆地的一口探井，主要目的层为长8段，其核磁共振测井解释与流体识别结果如图6所示。图中第4道为核磁共振T2分布，第5道为核磁共振双TW观测模式下反演的差谱，第7道为根据T2LM-Δφ图版输出的流体指示曲线，当其值大于0时进行填充。在该井2812～2815 m 处，差谱几何均值大于250 ms，有效孔隙度差约为0.5%，指示曲线数值大部分为“0”，个别深度点为“1”，表现为单“脉冲”，解释为水层。试油结果见油花，产水39.3 m3/d，为水层，验证了上述解释方法的准确性。

图7 Li189井核磁共振测井解释与流体识别结果

Li189井为研究区的另一口井，其核磁共振测井解释与流体识别结果如图7 所示。第7 道为利用T2LM-Δφ图版输出的流体指示曲线。该井2395～2 400 m 段，指示曲线数值大部分为“1”，个别为“0”，表现为连续“脉冲”，解释为油层，试油结果为产油15.56 t/d，产水0，为油层，证明上述解释是正确的。