高压压汞在致密气藏孔喉分布表征和早期产能评价中的应用

朱晴， 乔向阳， 张磊

(陕西延长石油(集团)有限责任公司研究院， 陕西 西安 710065)

准确评价储层岩石的孔喉(包括孔隙和喉道)结构特征，是渗流机理研究和产能预测的基础。目前分析岩石孔喉结构特征的主要方法包括：铸体薄片、扫描电镜、恒速压汞、氮气吸附、计算机断层扫描(CT)三维重建、核磁共振和高压压汞等，但不同方法皆存在一定的应用局限性。铸体薄片和扫描电镜仅能实现在一定尺度下对岩心剖面的观察和描述[1-3]；恒速压汞法不能测量半径小于0.1μm的孔喉[4]；氮气吸附法主要测定半径小于0.1μm的孔喉[5-7]；CT三维重建测试成本较高[8-9]；核磁共振无法单独实现定量表征，需联合其他实验[10]；高压压汞法不能测量半径小于0.003μm的孔喉。影响致密气藏开发的有效孔喉半径下限通常大于0.05μm[11]，CT三维重建、核磁共振和高压压汞等方法均能实现对有效孔喉半径的完整表征，但CT三维重建和核磁共振的测试成本较高，而高压压汞法的测试成本较低[12-13]，故被广泛采用。

高压压汞法的研究起源于20世纪40年代，最早由Purcell[14]在1949年提出了采用高压压汞测量毛管压力，随后国内外学者开展了大量关于高压压汞法的研究。例如，Wardlaw等[15-16]、Pittman[17]先后开展了毛管压力曲线的变化特征和通过毛管压力曲线评价孔喉尺寸相关性等方面的研究；唐仁骐等[18]、洪世铎等[19]开展了关于退汞效率的研究；Rashid等[20]、Salimifard等[21]、廖明光等[22]分别研究了高压压汞与渗透率的关系；Schmitt等[23]、Clarkson等[24]开展了联合高压压汞法和氮气吸附法表征孔喉结构的研究；Gane等[25]、Marschall等[26]、公言杰等[27]、宁传祥等[28]和李爱芬等[29]联合高压压汞和核磁共振对储层定量评价进行了研究；欧阳思琪等[30]、严强等[31]对高压压汞联合扫描电镜表征致密砂岩储层孔喉结构进行了研究。在油气储层的孔喉结构表征和物性特征研究等方面，高压压汞发挥了重要作用。

鄂尔多斯油气藏资源丰富[32-41]，东南部山西组致密气藏开发潜力巨大[42-46]，前人已开展了相关研究。如胡林楠[47]对东南部山西组致密砂岩的成藏规律进行了研究；秦波等[48]通过扫描电镜等手段对东南部山西组致密砂岩的有效储层识别和定量评价进行了分析；朱晴等[4]采用恒速压汞评价了东南部山西组致密砂岩的喉道特征；魏虎[49]采用高压压汞对东南部山西组致密砂岩的孔喉分布形态和孔隙结果分类进行了深入研究。目前高压压汞法在东南部山西组致密砂岩储层中的应用很少，通过高压压汞分析东南部山西组致密砂岩孔隙和喉道分布特征亦报道不多，基于准确表征孔喉分布的重要性和高压压汞评价孔喉分布的实用性，亟需开展高压压汞实验在鄂尔多斯盆地东南部致密砂岩储层中的应用研究。

本文选取了鄂尔多斯山西组致密砂岩储层岩心进行高压压汞实验，根据毛管压力曲线和孔喉半径分布曲线研究孔喉发育情况，分析孔喉特征参数和物性参数的主要控制因素，基于高压压汞实验结果，实现了完整孔喉结构的表征，进而提出了致密气藏早期产能评价的方法。

1 实验部分

1.1 实验样品

鄂尔多斯盆地东南部为典型的致密砂岩气藏，主力开采层位为山西组，根据山西组致密砂岩储层的物性特征，选取6口典型气井进行取心，取心深度为2158.21～2916.13m，取心层位为山西组1段(简称“山1段”)和山西组2段(简称“山2段”)，钻取柱塞样岩心，进行孔隙度和渗透率测试。本次实验岩心测试的基本参数如表1所示，孔隙度分布范围为5.84%～10.64%，渗透率分布范围为0.02～4.72mD。

表1 实验岩心基本参数

Table 1 Basic parameters of experiment cores

样品编号取心深度(m)层位孔隙度(%)渗透率(mD)岩心12158.21山1段9.880.02岩心22825.71山1段7.070.05岩心32672.64山2段6.070.25岩心42916.13山1段10.640.08岩心52708.70山2段6.070.82岩心62871.33山2段5.844.72

1.2 测试方法

对岩心样品1～6进行高压压汞实验，获取岩心的毛管压力曲线和孔喉半径分布。本次实验在西南石油大学油气田开发工程实验室进行，主要实验仪器为AutoPore Ⅳ 9505全自动压汞仪，设备最高实验压力为228MPa。实验流程按照国家标准《岩石毛管压力曲线的测定》(GB/T 29171—2012)进行，具体测试条件为：测试温度25.4℃，测试湿度38%RH，表面张力0.48N/m，润湿接触角140°。

2 结果与讨论

2.1 毛管压力曲线特征

图1a为实验岩心的进汞曲线。进汞曲线呈现出明显的三段式特征，进汞曲线反映了孔喉的整体发育及连通情况[5,29]：初始段的进汞饱和度反映出岩心的大尺度孔隙发育相对很少；中间段进汞曲线说明随着进汞饱和度增大，孔喉之间的连通性经历了变好再变差的过程；末端段的上翘特征表明不同岩心细小尺度孔喉的连通性差异显著。

图1 毛管压力曲线特征及不同压力区间的饱和度分布Fig.1 Capillary pressure curve characteristics and saturation distribution in different pressure intervals

图1b为实验岩心的退汞曲线。Wardlaw等[15]和严强等[31]提出了退汞曲线代表喉道的退汞过程，根据这一认知，残余汞饱和度代表孔隙体积占比，退汞饱和度代表喉道体积占比。实验结果显示：残余汞饱和度分布范围介于43.39%～70.43%之间，平均值56.18%；退汞饱和度介于9.43%～42.10%之间，平均值31.22%；残余汞饱和度和退汞饱和度差异均较大，表明不同岩心的孔隙和喉道发育情况差异较大。

借鉴十进制孔隙分类标准的思路[5,7]，采用十进制将毛管压力分成四个区间：100～1000MPa、10～100MPa、1～10MPa和0.1～1MPa。

图1c为进汞过程中不同毛管压力区间的进汞饱和度分布。结果表明：整体孔喉发育的毛管压力区间依次为0.1～1MPa、1～10MPa、10～100MPa和100～1000MPa，孔喉发育数量随着毛管压力增大呈现减少的趋势。毛管压力在0.1～1MPa区间孔喉整体最为发育，毛管压力在100～1000MPa区间孔喉整体最不发育。由此可见，0.1～1MPa和1～10MPa为鄂尔多斯盆地东南部山西组致密砂岩储层整体孔喉发育的主要毛管压力区间。

图1d为退汞过程中不同毛管压力区间的退汞饱和度分布。结果显示：单独喉道发育的毛管压力区间依次为0.1～1MPa、1～10MPa、10～100MPa和100～1000MPa，喉道发育数量随着毛管压力增大呈现减少的趋势。毛管压力在0.1～1MPa区间喉道整体最为发育，毛管压力在100～1000MPa区间喉道整体最不发育。由此可见，0.1～1MPa和1～10MPa为鄂尔多斯盆地东南部山西组致密砂岩储层单独喉道发育的主要毛管压力区间。

2.2 孔喉半径分布特征

图2a为孔喉半径分布特征曲线。结果显示全部岩心都呈现出典型的双峰特征[33]。左峰峰值半径分布范围介于0.09～2.87μm之间，平均值1.12μm；左峰峰值占比分布范围介于11.26%～33.18%之间，平均值21.40%。分析表明：不同岩心孔喉峰值分布两极分化现象明显，峰值半径越小的岩心其峰值占比越低，峰值半径越大的岩心其峰值占比越高，整体呈现出递增的趋势。

右峰与左峰完全分离，不同岩心的右峰跨度和右峰峰值半径均比较接近，右峰峰值占比分布范围介于0.96%～4.89%之间，平均值2.44%，左峰峰值比例远远超过右峰。对比表明鄂尔多斯盆地东南部山西组致密砂岩储层的小尺度孔喉发育差异很大，大尺度孔喉发育比较相似，表现出独特的孔喉发育特征。

基于全部岩心均发育双峰且两峰几乎完全分离，将左峰定义为小孔喉区，将右峰定义为大孔喉区，如图2b所示。采用积分计算左峰和右峰的面积，结果显示：小孔喉相对体积占比分布范围介于84.51%～94.88%之间，平均值91.59%；大孔喉相对体积占比分布范围介于5.12%～15.49%之间，平均值8.41%。表明鄂尔多斯盆地东南部山西组致密砂岩储层以小孔喉为主，小孔喉所占空间远远超过大孔喉。

2.3 孔喉分布对参数的影响

图3a和图3b为平均孔喉体积比的影响因素。平均孔喉体积比为退汞饱和度与残余汞饱和度的比值，代表了孔隙体积和喉道体积的比值。分析平均孔喉体积比的影响因素，结果显示：平均孔喉体积比与最大进汞饱和度和残余汞饱和度之间的相关性较差，与分选系数和退出汞饱和度之间的相关性很好。对比表明喉道决定了岩心的孔喉配比情况。

图3c和图3d为退汞效率的影响因素。退汞效率为退出汞饱和度与最大进汞饱和度的比值，在一定程度上反映了喉道体积占孔喉总体积的比例。分析退汞效率的影响因素，结果显示：退汞效率与最大进汞饱和度、残余汞饱和度和分选系数之间的相关性一般，与退出汞饱和度和平均孔喉体积比之间的相关性很好。表明喉道体积决定了岩石的退汞效率，随着孔喉发育情况越来越均匀，孔喉配比越接近1，退汞效率越高。

图3e和图3f为渗透率的影响因素。分析表明：渗透率与小孔喉体积占比和大孔喉体积占比之间均无明显的相关性，与小孔喉峰值半径和孔喉半径均值之间的相关性很好。由于小孔喉体积远超过大孔喉，小喉道决定了孔喉半径均值，渗透率代表了岩石的渗流能力，揭示储层渗流能力主要受到小孔喉发育程度的控制。

图3g和图3h为孔隙度的影响因素。分析表明：孔隙度与小孔喉区峰值占比和小孔喉区体积占比之间均无明显的相关性，与大孔喉区峰值占比和大孔喉区体积占比之间的相关性很好。与渗透率主要受到小孔喉区控制相反，孔隙度主要受到大孔喉区峰值占比和大孔区喉体积占比的影响，与小孔喉区参数之间的相关性较弱，孔隙度代表了岩石的储容性，表明储容性主要受到大孔喉区发育程度的控制。

图3 孔喉分布对参数的影响Fig.3 Effect of pore-throat distribution on parameters

3 实验结果的应用

3.1 完整孔喉分布表征

本次实验最大进汞压力为200.33MPa，对应的孔喉半径为3.67nm，表明小于3.67nm的喉道及其所控制的孔喉空间通过本次压汞实验无法测量，这亦是所有岩心最大进汞饱和度均小于100%的根本原因。

平均孔喉体积比为相对体积占比，结合最大进汞饱和度和孔隙度，分别计算得到小孔喉区和大孔喉区的绝对体积占比，拼接出岩心的高压压汞完整孔喉分布[7]，如图4所示。结果显示：岩心样品5的未进汞孔喉绝对体积占比最高，占岩心体积的2.86%；岩心样品2的大孔喉绝对体积占比最高，占岩心体积的0.60%；岩心样品4的小孔喉绝对体积占比最高，占岩心体积的9.89%。

计算各岩心的绝对残余汞体积和绝对退汞体积，绝对残余汞体积和绝对退汞体积分别代表岩心的孔隙绝对体积和喉道绝对体积，分析大小孔喉区绝对体积占比与孔喉绝对体积之间的相关性，如图5所示。结果显示：大孔喉区绝对体积占比与孔喉绝对体积之间无明显相关性，小孔喉区绝对体积占比与孔隙(图5a)和喉道(图5b)绝对体积之间具有较好的相关性，表明小孔喉区发育程度决定了储层孔隙和喉道的发育程度。

3.2 早期产能评价

统计高压压汞实验研究结果、可动流体饱和度和对应层位的压裂规模。表2结果显示可动流体饱和度随着大孔喉区绝对体积占比增加而增大，表明大孔喉区发育程度决定了储层的可动流体饱和度。根据渗流力学原理，气井投产早期产能主要来自压裂裂缝和大孔喉中可动流体的贡献，在压裂规模相当的情况下，早期产能主要由大孔喉区发育程度决定。实验岩心对应气井在相应层位的压裂液量相当，推测大孔喉区绝对体积占比越高的气井其早期产能越高，气井的试气效果证实了这一现象。气井钻完井未压裂试气之前，通过分析储层岩心大孔喉区的发育程度，可实现对气井早期产能的初步评价。

图4 岩心完整孔喉分布

Fig.4 Complete pore-throat distribution of cores

表2显示山2段储层的大孔喉区绝对体积占比普遍超过山1段储层，山2段储层的可动流体饱和度整体好于山1段储层，在压裂规模相当的情况下，揭示了鄂尔多斯盆地东南部山西组致密砂岩气藏山2段储层的整体试气效果好于山1段储层，这与目前开发实践相吻合。

图5 小孔喉区绝对体积占比对孔喉绝对体积的影响Fig.5 Effect of the absolute volume ratio in the small pore-throat area on the absolute volume of pore-throat

表2 实验结果和评价参数

Table 2 Experimental results and evaluation parameters

样品编号层位小孔喉区绝对体积占比(%)大孔喉区绝对体积占比(%)可动流体饱和度(%)压裂规模(m3)岩心1山1段89.854.8526.77365.2岩心2山1段79.748.4757.75356.5岩心3山2段88.166.1248.62367.5岩心4山1段92.955.0643.07360.6岩心5山2段44.648.1862.93399.7岩心6山2段88.108.2981.44366.1

4 结论

本文以鄂尔多斯东南部山西组致密砂岩储层岩心为例，开展了高压压汞实验在孔喉结构表征和早期产能评价中的应用研究。结果表明：不同岩心小尺度孔喉和喉道的发育特征差异较大，总体孔喉和单独喉道发育的主要毛管压力区间均为0.1～1MPa和1～10MPa，山西组致密砂岩储层以小孔喉为主，喉道决定了岩心的孔喉配比情况和退汞效率，小孔喉区发育程度控制了储层渗流能力，大孔喉区发育程度控制了储层储容性，孔隙和喉道的发育程度均由小孔喉区发育程度决定，在压裂规模相当的情况下早期产能主要由大孔喉区发育程度决定。

高压压汞实验为广泛采用的储层孔喉评价方法，采用退汞曲线代表喉道退汞过程的分析思路，拓展了利用高压压汞实验研究储层孔喉特征的评价范围。本文根据孔喉半径分布曲线划分出小孔喉区和大孔喉区，进而实现了完整孔喉结构的表征，丰富了对鄂尔多斯东南部山西组致密砂岩孔喉分布特征的认知。结合可动流体饱和度和压裂规模等参数，进一步提出了通过大孔喉区评价气井早期产能的方法，为致密气藏科学开发提供了新的思路。