微纳米孔隙网络中天然气充注的三维可视化物理模拟

乔俊程，曾溅辉，夏宇轩，蔡建超，陈冬霞，蒋恕，韩国猛，曹喆，冯枭，冯森

（1. 中国石油大学（北京）油气资源与探测国家重点实验室，北京 102249；2. 中国石油大学（北京）地球科学学院，北京102249；3. 中国地质大学（武汉）地球物理与空间信息学院，武汉 430074；4. 中国地质大学（武汉）构造与油气资源教育部重点实验室，武汉 430074；5. 中国石油大港油田公司，天津 300280；6. 中国石油化工股份有限公司石油勘探开发研究院，北京 102206；7. 中国石油集团工程技术研究院有限公司，北京 102206）

0 引言

随着全球能源需求加剧，非常规油气资源成为油气勘探开发的重要接替领域[1-3]。其中，低渗（致密）砂岩气是当前全球非常规油气勘探与开发的重要对象[4-5]，也是中国油气增储上产的重要战略目标[6-7]。低渗（致密）砂岩气的成藏特征及机理与常规油气藏有明显区别[8-10]，源储紧密交互的空间配置和微纳米级微小孔喉系统使得天然气从烃源岩初次运移进入储集层的充注过程成为低渗（致密）油气成藏的核心环节[11-14]。充注过程的流体流动机理决定了低渗（致密）砂岩中油气的充注、赋存和聚集程度，是充注机理研究的关键[11,15-17]。

针对低渗（致密）砂岩油气储集层流体流动机理，国内外学者从物理模拟、数值模拟和理论推导 3个方面开展了大量研究[15,18-22]。其中，物理模拟主要包括了沙箱模拟、岩心渗流模拟和孔隙尺度渗流物理模拟实验[23-29]。通过沙箱模型，前人在实验室从宏观砂体尺度研究了油气充注过程的流体流动，但其与实际储集层的物性差异大，研究结果难以真实反映充注过程的流体流动机制。因此，国内外学者开展了大量的真实岩心渗流实验[16-17,29-30]，指出低渗（致密）砂岩中油气的流动表现出明显的非达西流，存在明显的启动压力梯度[15,26-27]；Qiao等、Zeng等和曾溅辉等研究发现低渗（致密）砂岩油气成藏过程中的油气充注和运移也具有低速非线性流的特征[14,27,31]。许多学者通过理论推导，得到了描述低速非达西流的渗流方程，但考虑到流体渗流影响因素差别，所得到的渗流方程也具有差异[32-34]。近年来，一些学者通过微米X-CT与油气渗流实验装置相结合，深入开展了微观可视化油气多相渗流模拟；然而，这些研究以油气开发过程中的水驱油、气或渗吸过程为主[35-40]，且研究对象以常规砂岩或碳酸盐岩为主；尽管 Lin等在研究中关注到了非均质砂岩中气驱水的毛管驱替过程[38]，但很少有人关注油气充注成藏过程中油、气驱水多相流体流动。在此基础上，很多学者通过 Lattice Boltzmann、Monte Carlo等多场流体渗流数值模拟方法对低渗（致密）储集层中流体的动态流动特征进行了观测和研究[41]。

作为流动媒介，多孔介质的孔隙几何学和拓扑学特征控制并影响了流体流动和分布特征[23-24,42-44]，低渗（致密）砂岩发育的微纳米级孔喉系统是其流体流动和分布规律复杂的根本原因[22]。随着孔喉分析技术的发展，特别是扫描电镜、核磁共振和微纳米 X-CT等技术的应用，一些学者开始探索油气储集层中流体流动和分布规律及其主控因素[14,29,45-46]。Qiao等通过联合微米 X-CT与岩心渗流实验指出孔喉配置关系控制了致密气充注的流体流动特征[14]。公言杰等指出受微观孔隙结构类型控制，低渗（致密）砂岩中气体的非达西流动具备复合型、下凹型、上凹型和直线型 4模式[29]。然而，这些研究都是基于岩心流动实验与孔隙结构表征的机械耦合；一方面，跨尺度研究会造成系统误差；另一方面，实验结果仍是流体流动和孔隙结构的平均统计，所得出的规律难以直接反映流体在微纳米级孔喉中流动和分布状态及方式[31]，致使目前对于低渗（致密）储集层中流体流动和分布特征、规律及其主控因素的相关认识并不准确。因此，阐明低渗（致密）砂岩油气充注机理的关键在于明确充注过程中流体在微纳米孔喉中的流动与分布特征、规律及主控因素，这必须借助于高精度微纳米孔隙流体流动物理模拟方法[23-25]。

本研究通过微纳米孔隙三维可视化在线物理模拟实验开展天然气充注物理模拟，结合铸体薄片与扫描电镜分析，在系统表征低渗（致密）砂岩储集层微纳米孔喉网络结构特征的基础上，应用视渗透率理论、孔隙尺度原位叠算技术和孔隙网络模拟方法（Pore Network Modelling，PNM），对低渗（致密）砂岩气充注过程中微纳米孔隙网络中的气水流动与分布特征及其变化规律进行刻画，并在此基础上探讨气水流动及分布的微观影响因素，以期能够揭示低渗（致密）砂岩气充注机理，为低渗（致密）砂岩气的勘探开发提供借鉴和参考。

1 实验介绍

本文的低渗（致密）岩心样品来源于渤海湾盆地黄骅凹陷中部歧口凹陷歧北斜坡深层古近系沙河街组二段（简称沙二段）典型致密砂岩气储集层。歧口凹陷沙二段中烃源岩与储集层紧密交互[47-49]，实验岩心样品取自歧北斜坡板桥斜坡带BS35井沙二段3 790.49 m井深，位于致密气由烃源岩向储集层充注的源储界面位置，岩性为岩屑中砂岩，岩心长度为6.56 cm，直径为2.51 cm，气测孔隙度为 11.42%，渗透率为 3.94×10-3μm2，实测润湿角为 51.90°，是研究低渗（致密）砂岩气充注的一块代表性样品。

实验方法及其流程如图1所示，其中，关键的微纳米孔隙三维可视化在线天然气充注物理模拟的详细实验方法与流程如下：

图1 实验方法与流程

1.1 实验装置

实验采用中国石油大学（北京）油气资源与探测油气成藏物理模拟实验室研制的非常规油气储集层油气运移与聚集物理模拟实验装置（见图2a）。该装置包含Zeiss Xradia Versa 510 X-μCT扫描仪、CoreLab岩心驱替系统（Core Flooding System，CFS）及Peek岩心夹持系统（见图2a、图2b），微米X-CT最高分辨率为650 nm；CFS岩心驱替系统轴压与围压最高可达70 MPa，温度可达120 ℃；流体计量系统的计量精度为0.001 μL/min；岩心夹持系统的抗压强度大于45 MPa，2 cm的夹持器外径保证模拟实验成像分辨率可达约2 μm（见图2b），该系统实现了对真实地层温压条件下低渗（致密）砂岩储集层微纳米孔隙网络中流体流动的实时监测、成像及分析。

图2 非常规油气储集层油气运移与聚集物理模拟实验装置及原理示意图

1.2 实验流程

本次实验基于稳态法开展天然气充注物理模拟。实验全程在线，实现对样品的原位扫描（见图2c），实验步骤如下。

①对样品进行洗油、洗盐预处理，钻取直径为5 mm、长度为22 mm的柱塞样品，将其装入特制胶套，装填进Peek岩心夹持器中，将夹持器置入微米X-CT载物台，连接管线，接入CFS系统，检查实验系统的整体密封性。

②通过CFS系统向岩心夹持系统中加载合适的围压（2 MPa），对样品进行第1次原始状态X-CT扫描（分辨率为2.00 μm），用以获取样品原始的孔喉结构信息。

③对岩心夹持系统进行抽真空，排驱样品中空气，通过CFS系统向样品中持续注入质量分数为25%碘化钾（KI）溶液72 h以上，KI溶液对X射线的吸收能力强，在灰度图像上往往呈现出明显高于岩石基质和含气孔隙的灰度，可以有效区分岩石基质与孔隙中的气体和水，在样品达到地层水饱和状态后，进行 X-CT扫描。

④完成饱和水状态的X-CT扫描后，开始以较低的驱替压力（0.1 MPa）进行恒压驱替，根据稳态法的实验原理，待出口端稳定出气，气流速度稳定且与气体注入速度相等后，再次进行X-CT扫描。

⑤继续提高驱替压力确保每个驱替压力下出口的气体流速达到稳定后，在相应的驱替压力节点分别依次进行 X-CT扫描，直至提高驱替压力出口端不再出水，气体流速稳定与注入速度相等后停止实验。

⑥详细记录气体充注物理模拟实验过程中的流体流动数据，对各个驱替压力节点X-CT扫描结果进行数据分析。

1.3 数据分析

为准确获取不同充注阶段气水两相的流体流动和赋存特征，采用FEI Avizo 9.0.1图形处理软件对各个充注压力节点获得的 X-CT扫描图像进行孔隙尺度原位叠算和孔隙网络模拟，具体流程如下。

①使用 Zeiss重构软件对不同充注压力节点的X-CT原始扫描数据进行重构，准确校正图像中心位移，消除射线硬化，获得灰度断层图像；采用基于自相似图像结构像素灰度加权平均的三维非局部中值滤波（3D Non Local Median Filter）方法对灰度断层图像进行滤波，降低图像噪音，重复进行多次降噪，直至图像中颗粒表面光滑，颗粒、孔隙、气和水多相边缘和界面清晰，保证灰度分割的准确性（见图3a—图3c）。

图3 微米X-CT图像滤波、空间原位校正与原位叠算

②利用 FEI Avzio 9.0.1图像配准（Image Registration）模块，依据相邻充注压力节点灰度断层图像中特征矿物，通过校准特征矿物的微米级空间位移，使相邻充注压力节点的灰度断层图像在空间上处于原位。据此将所有充注压力节点灰度图像两两配准，之后进行图像的原位叠加计算（见图3a—图3c）。

③对原始干燥节点扫描获得灰度断层图像进行灰度阈值分割（Interactive Threshold Segmentation），以样品的气测孔隙度为标准，提取岩石中的孔喉空间，使 Avizo9.0.1中 Material Statistics模块计算得到的X-CT孔隙度与气测孔隙度基本一致，进而获得二值化孔喉网络模型图像；通过孔隙连通性分割，获取孔喉网络的三维连通模型；采用孔隙网络模拟技术，得到孔喉网络三维骨架模型，并计算得到相应的原始孔隙结构参数。

④在灰度图像空间位移校正的基础上，进行原位数字岩心叠算，以饱和 KI溶液的扫描图像为基准节点，利用 Avizo 9.0.1的 Arithmetic算法对充注压力0.10 MPa和 0.15 MPa节点的扫描图像进行差值计算（见图3d、图3e），将充注压力 0.10 MPa状态下与0.15 MPa状态下的KI溶液的分布图像（见图像中高亮部分）相减，可以得到充注压力0.15 MPa节点岩心孔隙中被气驱出的水相二值化图像，即该阶段充注进入岩心中气体的二值化图像（见图3f），将其与充注压力0.10 MPa节点的气体二值化图像叠加即可准确获得充注压力0.15 MPa节点孔喉中的气体分布；通过阈值分割充注压力0.15 MPa节点的高亮度KI溶液图像可以获得该节点的孔隙水分布二值化图像；将获得的空间气水分布叠合，可以得到该充注压力下孔喉中的空间气水分布状态（见图3g）。以此类推，可以分别获得各个充注压力下的气水分布状态与特征。基于气水图像的阈值分割，利用体积分数（Volume Fraction）方法可以分别计算岩心在饱和KI溶液状态下的含水体积分数和不同充注压力节点的含气体积分数，通过（1）式可以计算不同充注压力节点的含气饱和度。同时，对不同充注压力节点的气相分布图像进行孔隙网络模拟分析，可获得相应节点的含气孔喉的孔喉结构参数。

2 实验结果

2.1 微纳米孔喉网络孔隙结构特征

样品中存在原生孔隙和次生孔隙两类孔隙类型[50]，其中原生孔隙主要为原生残余粒间孔，孔隙形状规则，大小在几十至几百微米（见图4a）。次生孔隙包括了溶蚀粒间孔、溶蚀粒内孔、晶间孔和微裂缝，溶蚀粒间孔为颗粒边缘溶蚀形成的粒间孔隙，形成于长石和岩屑等可溶蚀性矿物的边缘，形状不规则，大小在几十至几百微米（见图4a、图4b）；溶蚀粒内孔形成于长石、岩屑等溶蚀性矿物内部，形状不规则，多在十几至几十微米（见图4b）；晶间孔为黏土矿物晶体之间形成的微小空间，大小介于几百纳米至几微米（见图4c）；微裂缝为切穿多个岩石颗粒的狭长孔隙，延伸较远，能够连接多个互不连通的孔隙（见图4a）。样品表现出与致密砂岩储集层相似的孔隙成因构成[50]，粒间孔隙发育，占比高达 80%左右，其中原生粒间孔和粒间溶蚀孔占比相近，在岩屑和长石颗粒中存在小规模的粒内溶蚀孔隙，占比约为 15%左右，发育少量的微裂缝和晶间孔（见图4a）。

图4 BS35井沙二段3 790.49 m井深低渗（致密）砂岩储集层孔隙成因类型

微米X-CT分析表明，空间上，形状规则的连通大孔隙构成了低渗（致密）砂岩的微纳米孔喉网络的主体，形状不规则的小孔隙一般靠近孔喉网络主体分布，连通程度低，而球状的微小孔隙既可分布在孔喉网络主体中，也可分布在远离网络主体的区域，其中位于孔网络主体中的微小孔隙往往能够桥接大孔隙和小孔隙，而远离网络主体的微小孔隙多呈孤立分布（见图5a—图5d）。基于Dong和Blunt开发的孔隙网络模拟方法[51]计算表明，样品的孔喉半径分布范围较广（见图5d），孔隙半径为 10～120 μm，集中分布在 20～60 μm，平均值为12.34 μm（见图5e），喉道半径为2～120 μm，集中分布在 10～30 μm，平均值为9.67 μm（见图5f），平均孔喉比为 1.28，平均配位数为 1.31，孔喉连通性相对较好（见图5d）。

图5 低渗（致密）储集层孔隙形态学与几何学特征

2.2 低渗（致密）砂岩气充注过程的气体流动特征

低渗（致密）砂岩气充注过程的天然气流动表现为非线性流，与低渗（致密）砂岩中的气体单相非达西流相似，存在明显的“启动压力平方梯度”，即天然气充注的临界压力平方梯度，为 0.01 MPa2/cm。流速随压力平方梯度增加而增加，流速增幅先快后稳，流动曲线可分为上凹段和线性段两段（见图6a）。含气饱和度表现出快速至稳定的二段式增长模式，初始阶段，含气饱和度随压力梯度增长快速增加，压力梯度超过0.29 MPa/cm，含气饱和度逐渐稳定（见图6b）。

图6 低渗（致密）砂岩气充注过程的气体渗流与含气饱和度变化特征

气相视渗透率多被广泛地应用于表征两相流体渗流中气体单相的流动特征变化[14,27]，气相视渗透率，如（2）式分析的气相视渗透率表明，充注过程中的气相视渗透率表现出随压力平方梯度增加先增加后稳定的两段式增长特征，与充注过程气体流动和含气饱和度的变化规律相符（见图6c）。

2.3 充注过程孔隙尺度气水分布与气相充注通道变化

从充注模拟实验中选取了岩心在原始状态下、饱和KI溶液状态下以及在充注压力分别为0.10，0.15，0.20，0.50 MPa的4种状态下共计6个特征节点进行分析（见图7），获得了不同充注压力节点下气水分布的二维和三维X-CT原始灰度图像。KI溶液使得饱和水的孔隙在灰度图像中显示为亮白色，岩石基质呈灰色，而气体充注进入孔隙后会使得存在气体的孔隙部分呈现黑色，与Armstrong等和Khishvand等学者所开展的孔隙尺度多相流模拟实验中饱和节点的成像特征基本一致[23,36]。基于岩石基质、气体和水溶液的灰度差异以及不同节点气水分布变化导致的图像灰度分布变化，通过孔隙尺度原位叠算技术和阈值分割技术，获得不同充注压力节点的微纳米孔隙网络中的气水空间分布图像（见图7g—图7l）；进一步地，对气相占据的孔喉空间进行分离提取（见图8a—图8d），利用孔隙网络模拟技术构建了不同充注节点气体充注通道的孔喉网络骨架模型（见图8e—图8h），并计算相应节点气体充注通道的孔隙结构参数。

图7 低渗（致密）砂岩气充注过程中二维与三维气水分布变化

图8 低渗（致密）砂岩气充注过程气相充注通道变化

从二维和三维原始灰度图像与三维气水分布图像可以看出，从饱和状态至充注压力0.10 MPa，天然气优先充注进入较大粒间孔隙的中央（见图7c、图9a、图9d），气体主要赋存在大孔隙中央，水溶液在大孔隙中呈薄膜状吸附在孔隙边缘；这与 Lin等在毛管驱替X-CT成像实验中观察到氮气优先驱替非均质砂岩中大孔隙的现象相似[38-39]，也同 Alizadeh等在含碳酸水实时 X-CT驱替成像实验中发现 CO2析出后优先占据孔喉中央的现象相似[35]；但本实验阶段仍有大量的水溶液在小孔隙中呈现饱和式充填的赋存形式。在空间上，气相和水相都呈网簇状分布（见图7i），计算含气饱和度为 40%。在该阶段，被气体充注的粒间大孔隙由较粗的喉道相连接，孔喉系统中形成了气体充注通道的初始基础格架（见图8a、图8e），基础格架中气体的连续流动和孔喉系统中含水饱和度降低使得气相视渗透率突变式增加（见图6c）。

图9 低渗（致密）砂岩气充注过程中不同孔径孔喉中气水流动与分布变化

随充注压力增至0.15 MPa，较大粒间孔隙边缘的KI溶液被驱出（见图7d、图9b、图9e），水膜厚度不断减小，大孔隙基本被气体完全饱和；同时，孔径较小的粒间孔隙中央的水溶液也被驱出，较小孔隙开始参与流动。空间上，网簇状气相的规模不断扩大，水相逐渐变为以薄膜式赋存为主（见图7j、图7k），孔喉网络的含水程度不断降低，含气饱和度增加至78%（见图6d），气相充注通道的规模不断扩张（见图8b、图8c）；相应地，骨架模型中主体格架的球体和圆柱体半径不断增加；同时，更多较小的球体和圆柱体与主体格架相连通（见图8f、图8g）。在气体充注通道扩张过程中，气体流速增大，含水饱和度降低，致使气相视渗透率快速增加。

充注压力增加至0.50 MPa，连通孔喉中的可动水基本被完全驱出（见图7f），即使充注压力再增加，气体也难以驱替连通孔隙边缘和微小粒内和晶间孔隙中的地层水（见图9c、图9f）孔喉网络中形成束缚水，水相在孔喉网络中整体表现出薄膜状零散分布，而气相则表现出集中式网簇状分布形式（见图7l），含气饱和度到达95%（见图6d），气相充注通道扩张至极限规模，天然气沿稳定的通道进行运移。该阶段孔喉系统中气体多表现为稳定的线性流动，气相视渗透率增幅逐渐降低并趋于稳定。稳态充注过程中，天然气充注通道呈现稳定连续的变化，这与 Spurn等在 CO2注入过程的快速同步辐射 X-CT扫描实验中发现的瞬态流动过程中间断性流体通道变化规律具有差异，而与其在稳态流动过程中发现的流体通道变化和含气饱和度变化规律十分相似[52]。

对不同充注压力节点气相充注通道孔隙结构参数的分析表明，充注通道的最大孔隙半径和喉道半径基本不随充注动力增加而变化，仅在充注压力0.15 MPa节点最大孔隙半径略有增加，最大孔隙半径基本稳定在116 μm左右，最大喉道半径基本稳定在78 μm左右，表明大孔喉在充注的起始阶段就基本完全参与流动，对气相充注通道的贡献到达极值（见图9g）。然而，平均孔隙半径、喉道半径、配位数和喉道长度随充注动力增加而降低（见图9h），表明随气相充注通道扩张，较小的孔喉对天然气充注的贡献逐渐增加。平均孔隙结构参数在充注压力为0.10～0.15 MPa阶段下降最明显，表明该阶段气相充注通道的扩张程度最大，也说明较小的孔喉主导了充注中间阶段气相充注通道的扩张幅度。可以推断，在低渗（致密）砂岩气充注过程中，粒间大孔喉主导了天然气充注的起始阶段，较小的粒间和连通粒内孔喉主导了充注中后期气相充注通道的扩张过程，决定了气相充注通道的极限规模（见图7—图9）。

3 讨论

3.1 微纳米孔喉结构对充注过程的影响

通过（3）式计算不同充注压力节点下不同孔径范围孔喉中的含气体积分数，可以得到充注过程中不同孔径孔喉中的含气性变化曲线，如图10a、图10b所示：

图10 不同充注压力下孔隙（a）和喉道（b）的含气体积分数曲线及其对不同充注压力下气相充注通道扩张的贡献率直方图（c、d）

通过对比可以发现，在充注的起始阶段（充注压力0.10 MPa），半径大于20 μm孔隙的含气体积占孔隙总体含气体积的90.35%，半径大于20 μm的喉道含气体积占总喉道含气体积的76.90%；其中，半径大于80 μm的孔隙含气体积分数在起始阶段就到达最大值，是充注初期的天然气流动方式、气水赋存和含气饱和度增长特征的重要微观影响因素。当充注压力增至0.15 MPa，含气体积分数的增长主要集中在半径为20～70 μm的孔隙和喉道中，该部分孔隙和喉道的含气体积分数增长分别占各自总含气体积增量的80.21%和80.31%（见图10c、图10d）；充注压力增至0.20 MPa，半径大于50 μm孔喉中基本没有含气体积分数增长，半径20～50 μm的孔喉主导了含气体积的增长，该部分孔隙和喉道含气体积分数增长分别占该阶段各自总含气体积分数增加量的85.42%和93.29%（见图10c、图10d），说明随充注动力增加，半径较小的孔隙逐渐成为气水分布变化和含气饱和度增长的主控因素；充注压力增至0.50 MPa，含气体积分数的增加主要来源于半径为3～50 μm的孔喉，其中孔隙增量占该阶段孔隙含气体积总增量的 84.43%，喉道增量占其含气体积总增量的85.86%（见图10c、图10d）。充注压力由0.15 MPa增至0.50 MPa，半径小于20 μm的孔隙和喉道的含气量持续增加，孔隙的含气体积分数增量从0.224%增加至0.428%，喉道的含气体积分数增量从 0.54%增加至0.65%；可见，半径小于20 μm的孔喉控制了充注后期的气水分布与含气饱和度变化（见图10）。

不同充注压力节点气相充注通道孔喉结构参数与气相视渗透率的相关关系分析表明，视渗透率变化与平均孔隙半径、喉道半径、喉道长度和配位数呈明显的负相关（见图11），表明充注过程气体流动状态变化受连通孔喉网络的整体孔喉结构特征的影响。随着充注程度提高，半径较小的连通孔喉不断参与气体流动，成为气相渗透率增长的重要影响因素。

图11 低渗（致密）砂岩气充注过程中视渗透率与气相充注通道孔喉结构参数变化的相关关系

总体而言，在低渗（致密）砂岩气充注过程中，半径大于20 μm的孔喉中，特别是半径大于80 μm的孔喉，先形成天然气充注的基础格架，是气体充注临界条件存在的重要因素，对充注初期的天然气流动方式、气水分布和含气饱和度增长特征具有重要影响；随充注动力的增加，半径为20～50 μm和小于20 μm的孔喉依次主导了气相充注通道的扩张过程，是主要充注期气水流动方式、气水分布变化、含气饱和度增长特征的主控因素之一；半径小于20 μm的孔喉是充注后期阶段稳定气水分布状态形成和最大含气饱和度的重要影响因素（见图10、图11）。

3.2 微纳米孔喉系统非均质性对充注及含气性的影响

对不同充注压力节点下岩心纵向切片的面含气饱和度的分析发现（见图12a），在充注起始阶段（充注压力为0.10 MPa），由于粒间大孔隙优先形成气相充注通道，纵向含气饱和度整体较为均匀，面含气饱和度以 40%为轴线上下浮动。在充注的中间过程中，纵向面含气饱和度产生明显差异；充注压力增至0.15 MPa，切片的面含气饱和度整体均增加，但不同切片区域增幅差异明显，表现为切片序号为175—400的区域含气饱和度增加明显高于其他部分，这与该区域连通的较大孔隙发育程度高密切相关（见图12b）；至充注压力0.20 MPa，整体面含气饱和度增加幅度较小，说明该阶段天然气主要充注进入连通小孔隙中聚集；至充注压力0.50 MPa，切片面含气饱和度产生一定增幅，可动水基本被驱替完全，达到最大含气饱和度，纵向上含气饱和度差异较充注过程中间阶段降低，面含气饱和度整体均匀，但切片序号为 50—150的区域面含气饱和度较低（见图12a、图12b）。连通大孔喉发育区域，面含气饱和度先快速增长后稳定，而连通小孔喉区域面含气饱和度增长规律十分复杂（见图12a、图12b）。孔隙结构分布的非均质性显著影响了天然气充注过程的气水赋存和含气饱和度的动态变化过程，对低渗（致密）储集层的含气非均一性也具有重要影响。

图12 低渗（致密）砂岩气充注过程中岩心纵向面含气饱和度变化曲线（a）与原始岩心孔隙纵向分布特征（b）

4 结论

低渗（致密）砂岩气充注过程存在气相充注通道扩张和稳定两个阶段。扩张阶段具有粒间大孔隙先于粒间小孔隙，孔隙中央先于边缘的气驱水流动模式；随充注压力增加，孔隙边缘和小孔隙中央的水被驱替，气相形成集中的网簇状分布，水相形成分离的薄膜状分布，气相充注通道规模扩张，通道的孔喉半径、喉道长度和配位数均降低，是含气饱和度增长的主要阶段；在稳定阶段，气相充注通道扩张至极限，通道的孔喉半径、喉道长度和配位数保持稳定，气体沿已形成的稳定通道运移，形成稳定的气水赋存状态，含气饱和度达到最大值。

充注起始阶段，半径大于20 μm的孔喉构成了气体充注通道的基础格架，是充注初期气体流动方式的重要影响因素；半径为30～50 μm的孔喉是气体的主要赋存空间，决定了该阶段气水分布和赋存特征及含气饱和度增长特征；随充注压力提高，半径为20～50 μm和半径小于20 μm的孔喉相继主导了气相充注通道扩张，控制了主要充注期的气水分布变化和含气饱和度增长模式，其中，半径小于20 μm的连通孔喉决定了气相充注通道扩张的极限，控制了稳定气水分布状态的形成和最大含气饱和度。

微纳米孔喉网络的非均质性对充注过程中的气水分布的动态变化过程具有重要影响，充注过程中，连通较大孔喉发育的部位含气性增加最明显。连通孔喉的非均质性对低渗（致密）储集层含气性不均一性具有重要影响。

符号注释：

A——岩心横截面积，cm2；Ka——气相视渗透率，10-3μm2；L——岩心长度，cm；p0——标准大气压，MPa；Δp——充注压力，MPa；Q——气体流量，mL/s；Sgi——第i个充注压力节点的含气饱和度，%；μ——气体黏度，mPa·s；φ——总孔隙度，%；φg——不同孔径孔喉的含气体积分数，%；φgi——i个充注压力节点的气体总空间体积分数，%；φp——不同孔径孔喉的体积分数，%；φwo——饱和水节点的孔隙水总空间体积分数，%。