2019, Vol. 31
2. 油气藏地质及开发工程国家重点实验室·西南石油大学, 成都 610500;
3. 延长油田股份有限公司, 陕西 延安 716000;
4. 中国石油集团测井有限公司 测井应用研究院, 西安 710077;
5. 中国石油集团渤海钻探工程有限公司 第一录井公司, 天津 300280
2. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China;
3. Yanchang Oil Field Co., Ltd., Yan'an 716000, Shaanxi, China;
4. Logging Application Research Institute, CNPC Logging, Xi'an 710077, China;
5. No.1 Mudlogging Company, Bohai Drilling Engineering Company, CNPC, Tianjing 300280, China
多孔介质中流体的微观渗流（尤其是两相渗流）对地下油气资源开发、河流堤坝防渗、环境保护、化学工程等领域均具有重要意义[1-5]。在油气资源开发领域，微观两相渗流对流体渗流机理、提高采收率等研究具有重要影响[6-8]。国内外学者针对岩石微观两相渗流已开展了大量的研究，所采用的研究手段分为物理实验观测方法[9-15]和数值模拟方法[16-21]两大类。
为了便于观察，目前的岩石微观两相渗流物理实验观测方法都是基于二维孔隙结构的，主要包括宏观驱替CT扫描实验法[9-12]、岩心切片微观驱替实验法[2]、玻璃刻蚀微观驱替实验法[3, 22-32]等。宏观驱替CT扫描实验法是利用CT扫描技术对柱塞岩心的流体驱替过程进行观察，可以获得流体在柱塞岩心中的分布及驱替优势通道，但驱替过程中所采用的CT扫描技术分辨率往往较低，无法清晰地观察到孔隙尺度的微观驱替过程及渗流特征，难以满足油气开发精细研究的需要。随着岩石两相渗流机理研究的不断深入，有学者针对宏观驱替实验的不足，提出了岩心切片微观驱替实验法和玻璃刻蚀微观驱替实验法等微观驱替实验方法[2-3, 22-32]。这2种微观实验方法采用的驱替和观测手段相似，主要区别在于驱替和观测的对象（物理模型）不同。岩心切片微观驱替实验法直接将岩心切成薄片再进行驱替和观测，能尽量保留岩石的孔隙结构和润湿性特征，且制作技术难度较小，成本低，但模型透明程度差，无法清晰地观测流体的微观渗流特征[2]。玻璃刻蚀微观驱替实验法则是基于岩心铸体薄片制作透明的微观孔隙网络模型再进行驱替和观测，具有较强的可视化能力，便于直接观测微观孔隙中流体的动态渗流特征，因此在水驱油的机理研究中运用非常广泛[31-32]。近年来，众多学者在微观驱替实验的理论和技术上取得了不少成果[2, 27-29, 31]，针对常规岩心驱替实验和玻璃刻蚀微观驱替实验，不仅实现了显微镜下直接观察两相流体渗流动态特征，还进行了一些两相流体渗流参数定量分析的有益尝试。
近些年，随着微观渗流理论的进步及计算机性能的提升，岩石微观渗流数值模拟方法也迅速发展。目前使用较多的岩石微观渗流数值模拟方法主要包括格子玻尔兹曼方法[16-17]和基于Navier-Stokes方程的数值模拟方法[34-36]。这2种方法均以岩石微观孔隙网络为基础，对流体的运动进行模拟。格子玻尔兹曼方法起源于20世纪80年代中期，是一种以格子气自动机为基础发展起来的流场模拟方法，该类方法采用离散介质模型来模拟流体的运动。基于Navier-Stokes方程的数值模拟方法则将流体视为不间断的整体，用微积分方法来计算流体运动的参数，具体的计算方法则包括有限差分法（FDM）[18-19]和有限元法（FEM）[34-37]。一般来说，有限差分法是以微分方程为基础，通过离散化处理求解区域后，用差分和差商近似代替微分和微分商[18]，因此求解微分方程和边界条件，可以归结为求解一个线性方程组以得到数值解。有限差分法具有直观、网格划分简单等优点[18]，但在某些情况下对于复杂的边界是不适用的，由于有限差分法的缺点，它几乎被网格划分更加灵活、边界适应性更强的有限元法所取代，因此有限元数值模拟在流体渗流模拟中也逐渐得到广泛的应用[18]。
现有的岩石微观两相渗流研究主要针对油水两相流体，相对而言，针对气水两相流体的岩石微观渗流研究还比较少，并且到底气水两相流体的岩石微观渗流研究更适合物理实验观测方法还是数值模拟方法尚无定论。笔者利用相同的孔隙结构模型，分别采用玻璃刻蚀微观驱替实验法和基于Navier-Stokes方程的有限元数值模拟方法对气水两相流体的微观渗流特征进行研究，并对比分析物理实验观测方法和数值模拟方法的优缺点，以期为气水两相流体的岩石微观渗流机理研究和气藏的精细高效开发提供参考。
1 玻璃刻蚀驱替实验方法 1.1 玻璃刻蚀模型制备在油气田开发实验中所采用的微观模型种类很多，微观玻璃刻蚀模型是微观模型中常用的模型之一。本文基于岩心铸体薄片照片构建玻璃刻蚀模型，如图 1所示：第一，将岩石铸体薄片[图 1（a）]在ImageJ软件中转换为灰度图[图 1（b）]；第二，将灰度图在Avizo软件中作分割处理，划分出孔隙和岩石颗粒[图 1（c）]；第三，基于分割图在Coreldraw软件中绘制刻蚀模型图像[图 1（d）]；第四，将刻蚀模型图像刻蚀到玻璃上，制作玻璃刻蚀模型[图 1（e）]。微观玻璃刻蚀模型的制备工艺（尤其是光刻技术）直接决定了玻璃刻蚀模型的精度，进而影响渗流实验的结果。由于目前光刻技术主要应用于光电技术领域，在油气资源开发中应用较少，一些学者[38-40]在微观玻璃刻蚀模型制作理论与技术上做了深入研究，为微观渗流研究提供了理论指导。本文所采用的微观玻璃刻蚀模型制备工艺的具体操作步骤分3步：①清洗匀胶。打开真空泵，吸住清洗后的玻璃片；先在玻璃片上滴打底胶作为增粘剂；烘烤10 min后滴上光刻胶，匀胶30 s，得到一定厚度的光刻胶膜。此处采用底胶和刻蚀胶配合使用的办法，能有效防止刻痕浅和孔隙“V”形槽变形。②曝光显影。在烤胶板上进行前烘10 min，然后进行曝光；利用光刻机进行曝光，曝光时间为80~120 s；将曝光后的玻璃片，放入显影液中，显影时间为5 min，后烘时间为10 min。③刻蚀去胶。在烤胶板上后烘并封蜡，在酸蚀液中腐蚀5~15 min；溶解经光照后的光刻胶，获得刻蚀模型。
|
下载eps/tif图 图 1 微观玻璃刻蚀模型制作流程图 Fig. 1 Flow chart for fabrication of micro glass etching model |
目前，基于微观玻璃刻蚀模型的驱替实验主要针对油水两相流开展研究，油和水的颜色差异较大，可以比较清晰地对实验过程进行观测[24, 27, 41-42]。对于气水两相流研究，气体和水的颜色相近，因此气水两相流体的界面不易观察。为了更好地观测驱替实验过程，采用亚甲基蓝染料对水进行染色，以区分水和气体（氮气）。准备好驱替模型和驱替流体后，即可开始驱替实验。驱替实验装置如图 2所示，具体实验步骤如下：①将微观玻璃刻蚀模型安装固定好，给玻璃刻蚀模型施加合适的围压，并对模型抽真空，再从入口端注入水（亚甲基蓝染色水）直至模型完全饱和水。②在24 ℃的恒温条件下，在入口端启动微量恒速泵，将初始驱替速度设置为0.05 mL/min，将气体（氮气）持续（20 min）注入玻璃刻蚀模型孔隙中。③当气水分布达到稳定状态后，逐渐加大入口端流量，分别以0.1 mL/min流量驱替20 min，0.2 mL/min流量驱替10 min，0.5 mL/min流量驱替5 min等，直至出口基本无水流出，驱替结束。整个实验过程中利用微摄像机拍摄实验过程的图像。
|
下载eps/tif图 图 2 气水两相流玻璃刻蚀驱替实验装置示意图 1.空压机；2.压力缓冲罐；3-5.中间容器；6.六通阀；7.光源；8.微观模型及加持容器；9.摄像头；10.采集电脑；11.录像系统；12.集液容器 Fig. 2 Schematic diagram of glass etching displacement experiment device for gas-water two-phase flow |
不同驱替时刻的驱替图像可以由图像处理系统得到，为了使两相流界面更易区分及分析渗流特征，需对图片作进一步处理（图 3）：首先将实验得到的驱替图片导入ImageJ软件，转换成灰度图；再将灰度图导入Avizo软件，进行分割处理，在气体和水所占孔隙空间填充不同颜色加以区分。
|
下载eps/tif图 图 3 驱替图像处理流程示意图 Fig. 3 Flow diagram of displacement image processing |
二维有限元模型构建（图 4）与玻璃刻蚀模型一样，也是基于岩心铸体薄片照片[图 4（a）]，先经过灰度图转化[图 4（b）]和分割处理[图 4（c）]这2个步骤。有限元模型制作的第三步，是利用AutoCAD软件，基于分割后的孔隙网络图按照真实比例构建孔隙结构模型[图 4（d）]；第四步是将孔隙结构模型导入COMSOL Mutilphysics软件来定义求解域[图 4（e）]，模型左端为入口，右端为出口，闭合曲线所围区域为岩石颗粒，其余连通区域为孔隙和喉道，连通区域设置为求解域；第五步是采用自由三角形对求解域进行网格划分[图 4（f）]，小喉道网格划分较为密集，大孔隙中间网格划分相对稀疏，越靠近孔隙壁网格划分越密集。
|
下载eps/tif图 图 4 有限元模型构建流程示意图 Fig. 4 Flow chart of finite element model construction |
二维单相流体渗流数值模拟实验采用缓慢的蠕动流接口求解Stokes方程。由于模型尺寸较小，所以在数值模拟过程中认为孔隙内流体是恒温等密度流。假设流体不可压缩，且不考虑流体的重力影响。模型入口在左端，出口在右端。在恒定压力条件下，流体从模型入口流入，单相流体渗流Stokes方程和连续性方程可表示为
| $ {\nabla \left\{ { - P\mathit{\boldsymbol{I}} + \mu \left[ {\nabla \mathit{\boldsymbol{u}} + {{(\nabla \mathit{\boldsymbol{u}})}^T}} \right]} \right\} + \mathit{\boldsymbol{F}}=0} $ | (1) |
| $ {\rho \nabla (\mathit{\boldsymbol{u}})=0} $ | (2) |
式中：P为压力，Pa；I为单位矩阵；μ为流体动态黏度，Pa·s；u为流体速度向量，m/s；F为体积力向量，N/m3，ρ为流体密度，kg/m3；T为绝对温度（293.15 K）。
边界条件：模型左端各孔隙入口压力相等且恒定不变，右端各孔隙出口压力也相等且恒定不变；除出入口外的其他边界均无流体的流入和流出；孔隙壁为中性润湿，流体流速在孔隙壁方向的法向分量为零。
2.3 两相流模拟方法气水两相流数值模拟实验采用Navier-Stokes方程并结合Level set方法构建微观两相流数学模型，并利用有限元方法求解方程。模拟实验采用层流两相流接口，在模型中，指定材料1（润湿相）为水，饱和孔隙和喉道，材料2（非润湿相）为气体，气体从模型左端入口注入，驱替孔隙内的水，水由出口流出。用水平集函数φ定义气水两相流各自的体积分数，为水时，φ=0，为气体时，φ=100%，并且利用重新初始化的水平集函数φ值追踪两相流界面，研究气水两相流界面动态特征。气水两相流界面的动态方程可表示为
| $ \frac{{\partial \varphi }}{{\partial t}} + \mathit{\boldsymbol{u}}\varphi = \gamma \nabla [\varepsilon \nabla \varphi - \varphi (1 - \varphi )]\frac{{\nabla \varphi }}{{|\nabla \varphi |}} $ | (3) |
式中：φ为水平集函数；t为两相流驱替时间，s；γ是初始化参数（默认设置为1 m/s）；ε为控制界面厚度的参数（设置为hmax /2，hmax为两相流界面流经区域的最大网格单元的大小）。
除了能定义流体界面，水平集函数φ也被用于描述流体的特性（密度和黏度跃变），即
| $ {\rho = {\rho _{{\rm{ air }}}} + \left( {{\rho _{{\rm{ water }}}} - {\rho _{{\rm{ air }}}}} \right)\varphi } $ | (4) |
| $ {\mu = {\mu _{{\rm{ air }}}} + \left( {{\mu _{{\rm{ water }}}} - {\mu _{{\rm{ air }}}}} \right)\varphi } $ | (5) |
式中：ρ为密度，kg/m3；ρwater和ρair分别为水和气体的密度，kg/m3；μwater和μair分别是水和气体的动态黏度，Pa·s。
流体的流速与声音的速度相比特别小，可以忽略不计，所以可以假设水和气体均为不可压缩的。Navier-Stokes方程可以用来描述不可压缩两相流体的质量以及动量的转输，并且考虑了界面张力，即气水两相流动的Navier-Stokes方程和连续性方程为
| $ {\rho \frac{{\partial \mathit{\boldsymbol{u}}}}{{\partial t}} = \nabla \left\{ { - p\mathit{\boldsymbol{I}} + \mathit{\boldsymbol{\mu }}\left[ {\nabla \mathit{\boldsymbol{u}} + {{(\nabla \mathit{\boldsymbol{u}})}^T}} \right]} \right\} + \rho \mathit{\boldsymbol{g}} + {\mathit{\boldsymbol{F}}_{{\rm{st}}}} + \mathit{\boldsymbol{F}}} $ | (6) |
| $ {\nabla \mathit{\boldsymbol{u}} = 0} $ | (7) |
式中：g为重力，N/m；Fst为气水界面张力，N/m。
两相流界面上的表面张力可表示为
| $ {{\mathit{\boldsymbol{F}}_{{\rm{st}}}} = \nabla T} $ | (8) |
| $ {T = \sigma \left( {\mathit{\boldsymbol{I}} - \mathit{\boldsymbol{n}}{\mathit{\boldsymbol{n}}^T}} \right)\delta } $ | (9) |
式中：σ为界面张力系数，N/m；n为垂直于界面的单位向量；δ为界面上的狄拉克函数，1/m。
垂直于界面的单位向量为
| $ \mathit{\boldsymbol{n}} = \frac{{\nabla \varphi }}{{|\nabla \varphi |}} $ | (10) |
狄拉克函数δ可以根据如下光滑函数近似计算得到
|
(11) |
川西地区蓬莱镇组致密砂岩储层按照物性及孔隙结构类型可分为3类（图 5）：Ⅰ类储层，孔隙度＞11%，渗透率＞0.8 mD，核磁共振呈明显双峰分布，可动峰幅度较高，压汞排驱压力低于0.3 MPa，孔喉类型为中孔-细喉道[图 5（a）]。Ⅱ类储层，孔隙度为9%~11%，渗透率为0.1~0.8 mD，核磁共振呈明显双峰分布，可动峰幅度中等，压汞排驱压力为0.3~1.0 MPa，孔喉类型为中孔-大喉道[图 5（b）]。Ⅲ类储层，孔隙度为6%~9%，渗透率为0.03~0.10 mD，核磁共振双峰分布不明显，可动峰幅度较低，压汞排驱压力为1~2 MPa，孔喉类型为小孔-大喉道[图 5（c）]。针对Ⅰ类储层、Ⅱ类储层和Ⅲ类储层各选取了一张典型铸体薄片照片，同时开展玻璃刻蚀驱替实验、单相流渗流模拟和气水两相流渗流模拟。
|
下载eps/tif图 图 5 川西地区蓬莱镇组致密砂岩储层孔隙结构特征 Fig. 5 Pore structure characteristics of tight sandstone reservoir of Penglaizhen Formation in western Sichuan |
微观可视化玻璃刻蚀驱替实验选用3张具有不同孔隙结构类型的铸体薄片制作玻璃刻蚀模型进行气驱水实验，研究气水两相流渗流特征及驱替后的气、水分布规律。玻璃刻蚀气驱水实验初始阶段，水饱和孔隙空间，3个模型均是左端为入口，右端为出口；将气体由左端入口注入模型进行驱替实验，随着实验的进行，被驱替的水由右端出口流出，气体与水的基本参数如表 1所列。
|
下载CSV 表 1 气体和水的参数 Table 1 Parameters of gas and water |
在气驱水实验中，由于孔隙和喉道的半径存在微观非均质性，不同孔喉半径具有不同的毛管压力，气体驱替时所受毛管阻力也不同，驱替前缘速度不同，存在明显的微观指进现象。即将气体由左端入口注入模型，当驱替压力大于喉道的毛管阻力时，气体率先通过喉道；当驱替压力小于毛管阻力时，气体先停止向前移动一段时间，然后发生运移跳跃；驱替前缘沿着毛管阻力较小的连通路径向前推进并率先到达出口，形成优势通道。当气体到达模型出口后，注入的气体多半沿着优势通道流动。从图 6可看出，Ⅰ类储层驱替前缘向前推进较为均匀，指进现象较弱，驱替效果较好[图 6（a）]；Ⅱ类储层模型出口端见气时，含水饱和度为56.87%，最终含水饱和度为52.29%[图 6（b）]；Ⅲ类储层模型出口端见气时，含水饱和度为57.80%，最终含水饱和度为46.25%，含水饱和度降低幅度较小，驱替效果差，指进现象较为明显[图 6（c）]。实验结束后，剩余水主要分布在小喉道、小喉道所控制的大孔隙及孔隙盲端。
|
下载eps/tif图 图 6 玻璃刻蚀驱替实验过程中不同驱替状态的两相分布图 Fig. 6 Two-phase flow distribution of different displacement states in glass etching displacement experiment |
构建的Ⅰ类、Ⅱ类、Ⅲ类储层模型的尺寸分别为11.70 mm×10.98 mm，11.76 mm×10.97 mm，12.21 mm×10.93 mm。单相流数值模拟时，流体在恒压条件下从模型左端入口被注入到孔隙空间，气体首先占据大喉道和大喉道控制的孔隙，接着气体逐渐填充小喉道和小喉道控制的孔隙，随着时间的推移，模拟结果逐渐接近于稳态。如图 7（a）—（c）所示，模型中流体的压力从左（入口）向右（出口）逐渐降低，由于3类储层模型孔隙结构不同，毛管阻力不同，压力等值线分布也就不同。Ⅰ类储层模型中压力大于0.67 kPa的区域较大，Ⅱ类储层模型中压力为0.57~0.67 kPa的区域较大，Ⅲ类储层模型中压力为0.26~0.57 kPa的区域较大。Ⅰ类储层模型比Ⅱ类储层模型和Ⅲ类储层模型具有更大的孔隙半径，因此，它的压力沿着驱替的方向降低的最慢。如图 7（d）—（f）所示，在流速图中，蓝色表示低流速，红色表示高流速。通道中心处的流速比入口处的更大，越靠近孔隙壁流速越小。Ⅰ类储层模型的孔喉半径最大，所以毛管阻力更小，流动速度更大。
|
下载eps/tif图 图 7 3类储层模型单相流数值模拟结果 Fig. 7 Numerical simulation results of one-phase flow of three types of reservoir models |
为了便于实验结果对比分析，气水两相流数值模拟实验采用了和微观玻璃刻蚀气驱水实验相同的孔隙结构模型，数值模拟的模型出入口条件、气体和水的相关参数也都与玻璃刻蚀驱替实验保持一致，气体和水的相关参数如表 1所列。
气驱水数值模拟实验过程中存在明显的微观指进现象，即气体沿着一条或多条阻力较小、由大喉道和连通孔隙组成的连续路径向前运移并率先到达出口（图 8—图 10）。模拟实验初始阶段，气体由左端注入，当驱替前缘到达模型左端时，模型中Ⅰ类、Ⅱ类、Ⅲ类储层的含水饱和度分别为34.69%[图 8（c）]、56.87%[图 9（e）]和57.80%[图 10（e）]。Ⅰ类储层模型比Ⅱ类和Ⅲ类储层模型具有更大的孔喉半径，驱替时间更快，当两相流界面前缘到达模型出口后，驱替范围继续变大（图 8）。由于Ⅱ类储层和Ⅲ类储层模型孔喉半径较小，毛管阻力较大，模型右端出口见气后，流动通道中的驱替阻力减小，被注入的气体基本沿着优势通道向前运移，波及范围不再扩大，驱水效率基本不再增加，指进现象更为明显（图 9，图 10）。
|
下载eps/tif图 图 8 Ⅰ类储层模型气水两相流数值模拟实验驱替过程 Fig. 8 Displacement process of gas-water two-phase flow numerical simulation experiment of type Ⅰ reservoir model |
|
下载eps/tif图 图 9 Ⅱ类储层模型气水两相流数值模拟实验驱替过程 Fig. 9 Displacement process of gas-water two-phase flow numerical simulation experiment of type Ⅱ reservoir model |
|
下载eps/tif图 图 10 Ⅲ类储层模型气水两相流数值模拟实验驱替过程 Fig. 10 Displacement process of gas-water two-phase flow numerical simulation experiment of type Ⅲ reservoir model |
为了对比玻璃刻蚀驱替实验结果和有限元数值模拟结果的异同点，对玻璃刻蚀实验和数值模拟实验的驱替时间进行归一化处理。玻璃刻蚀驱替实验和有限元数值模拟结果对比表明：
（1）归一化后的数值模拟结果与玻璃刻蚀驱替实验结果基本一致（图 11）。相同含水饱和度条件下，数值模拟与玻璃刻蚀驱替实验的两相流体界面前缘基本相同[图 6（h），图 10（e）]，数值模拟和玻璃刻蚀驱替实验都存在明显的指进现象。当驱替前缘到达模型入口端时，Ⅰ类、Ⅱ类和Ⅲ类储层模型的含水饱和度分别为34.69%，56.87%和57.80%，这与图 11中曲线的拐点基本一致。数值模拟和玻璃刻蚀驱替实验中，Ⅰ类储层模型比Ⅱ类和Ⅲ类储层模型的孔喉半径更大，驱替时间更快。当两相流界面前缘到达模型出口后，Ⅰ类储层模型驱水效率继续增加，驱替范围继续变大，Ⅱ类和Ⅲ类储层模型驱水效率基本不再增加，波及范围不再扩大。
|
下载eps/tif图 图 11 玻璃刻蚀驱替实验与有限元数值模拟的含水饱和度变化规律对比 Fig. 11 Comparison of water saturation variation between glass etching displacement experiment and finite element numerical simulation |
（2）数值模拟和玻璃刻蚀驱替实验均为微观可视化实验。玻璃刻蚀驱替实验可直接观测整个连续的驱替过程，比如气体在通过狭窄喉道时的运移跳跃现象。数值模拟则能较好地观测不同时刻的动态驱替特征。在驱替过程的初始阶段，气体率先驱替分布在大喉道及与大喉道相连的孔隙内的水，在驱替结束后，残余水主要分布在小喉道、小喉道控制的大孔隙和孔隙盲端。
（3）微观玻璃刻蚀驱替实验方法能较好地观测气体在通过狭窄喉道时的运移跳跃现象，但玻璃刻蚀驱替实验受仪器和刻蚀工艺精度不够等不利因素的影响，容易导致驱替压力不稳定、刻蚀成型的孔隙尺寸与设计尺寸有一定误差等现象。整体而言，玻璃刻蚀模型制备的制作成本高，模型精度较低，孔喉容易变形、成功率太低，其工艺技术还需进一步提高。
（4）数值模拟不受仪器和刻蚀工艺精度不够等不利因素的影响，可完美地重现岩石的微观孔隙结构特征，但数值模拟只考虑了有限的边界条件，因此无法准确模拟一些特殊的微观渗流现象。整体而言，有限元数值模拟方法则具有操作简单、可重复实验、运算能力强、成本低的特点，为微观可视化渗流研究提供了一种新方法。
（5）目前，玻璃刻蚀驱替实验和有限元数值模拟各有优缺点，结合二者开展气水两相渗流研究的效果更佳。
5 结论（1）微观玻璃刻蚀驱替实验方法与有限元数值模拟方法均可用于岩石孔隙中气水两相流体的渗流研究，在采用相同孔隙结构模型的前提下，这2种方法所得到的流体驱替前缘形态和渗流规律大致相似。
（2）气水两相流玻璃刻蚀驱替实验过程中存在明显的指进现象；驱替结束后，分布在大孔喉和与大孔喉相连的孔隙中的水被驱替，未被驱替的水主要分布在小孔喉、小孔喉控制的大孔隙和孔隙盲端。
（3）在单相渗流模拟实验达到稳态时，流动通道中心的流速高于入口端，孔隙半径越大，流动通道中心的速度越高，压力下降越慢。两相流数值模拟实验过程中也存在明显的指进现象，气体沿着由大喉道及大喉道控制的孔隙组成的连续通道向前驱替，达到驱替稳态后，残余水主要分布在小孔喉、与小孔喉连接的大孔隙及孔隙盲端。
（4）微观玻璃刻蚀驱替实验方法能较好地观测微观指进现象和气体在通过狭窄喉道时的运移跳跃现象，但玻璃刻蚀模型制备时，孔喉容易变形、成功率太低。有限元数值模拟方法则具有操作简单、可重复实验、运算能力强、成本低的特点，为微观可视化渗流研究提供了一种新的选择。
| [1] |
赵习森, 党海龙, 庞振宇, 等. 特低渗储层不同孔隙组合类型的微观孔隙结构及渗流特征:以甘谷驿油田唐157井区长6储层为例. 岩性油气藏, 2017, 29(6): 11-17. ZHAO X S, DANG H L, PANG Z Y, et al. Microscopic pore structure and seepage characteristics of different pore assemblage types in ultra low permeability reservoir:a case of Chang 6 reservoir in Tang 157 well area, Ganguyi Oilfield. Lithologic Reservoirs, 2017, 29(6): 11-17. |
| [2] |
KELLER AA, BLUNT M J, ROBERTS A P V. Micromodel observation of the role of oil layers in three-phase flow. Transport in Porous Media, 1997, 26(3): 277-297. DOI:10.1023/A:1006589611884 |
| [3] |
于明旭, 朱伟耀, 宋洪庆. 低渗透储层可视化微观渗流模型研制. 辽宁工程技术大学学报(自然科学版), 2013, 32(12): 1646-1650. YU M X, ZHU W Y, SONG H Q. Development of microscopic visualization flow model of low-permeability reservoir. Journal of Liaoning Technical University(Natural Science), 2013, 32(12): 1646-1650. DOI:10.3969/j.issn.1008-0562.2013.12.014 |
| [4] |
何文祥, 杨亿前, 马超亚. 特低渗透率储层水驱油规律实验研究. 岩性油气藏, 2010, 22(4): 109-115. HE W X, YANG Y Q, MA C Y. Experimental study on waterflooding in ultra-low permeability reservoirs. Lithologic Reservoirs, 2010, 22(4): 109-115. DOI:10.3969/j.issn.1673-8926.2010.04.020 |
| [5] |
李继庆, 刘曰武, 黄灿, 等. 页岩气水平井试井模型及井间干扰特征. 岩性油气藏, 2018, 30(6): 138-144. LI J Q, LIU Y W, HUANG C, et al. Multi-stage fracturing horizontal well interference test model and its application. Lithologic Reservoirs, 2018, 30(6): 138-144. |
| [6] |
李友全, 韩秀虹, 阎燕, 等. 低渗透油藏CO2吞吐压力响应曲线分析. 岩性油气藏, 2017, 29(6): 119-126. LI Y Q, HAN X H, YAN Y, et al. Pressure transient analysis on CO2 huff and puff in low permeability reservoir. Lithologic Reservoirs, 2017, 29(6): 119-126. DOI:10.3969/j.issn.1673-8926.2017.06.015 |
| [7] |
邓学峰. 致密低渗油藏压裂水平井合理生产压差优化设计. 岩性油气藏, 2017, 29(1): 135-139. DENG X F. Optimization of reasonable production pressure difference of fractured horizontal well in low permeability tight reservoirs. Lithologic Reservoirs, 2017, 29(1): 135-139. DOI:10.3969/j.issn.1673-8926.2017.01.018 |
| [8] |
严向阳, 王腾飞, 何双喜, 等. 过量顶替液作业下压裂水平气井的产能模拟. 岩性油气藏, 2017, 29(1): 140-146. YAN X Y, WANG T F, HE S X, et al. Productivity simulation for fracturing horizontal gas wells considering overdisplacing operation. Lithologic Reservoirs, 2017, 29(1): 140-146. DOI:10.3969/j.issn.1673-8926.2017.01.019 |
| [9] |
吕伟峰, 冷振鹏, 张祖波, 等. 应用CT扫描技术研究低渗透岩心水驱油机理. 油气地质与采收率, 2013, 20(2): 87-90. LYU W F, LENG Z P, ZHANG Z B, et al. Study on waterflooding mechanism in low-permeability cores using CT scan technology. Petroleum Geology and Recovery Efficiency, 2013, 20(2): 87-90. DOI:10.3969/j.issn.1009-9603.2013.02.022 |
| [10] |
曹永娜. CT扫描技术在微观驱替实验及剩余油分析中的应用. CT理论与应用研究, 2015, 24(1): 47-56. CAO Y N. Application CT scanning technology analysis microflooding experiments and the residual oil. Computerized Tomography Theory and Applications, 2015, 24(1): 47-56. |
| [11] |
张虹, 董国臣, 王松洋, 等. 浅析显微CT技术在油气储层研究中的应用. 资源与产业, 2013, 15(6): 98-102. ZHANG H, DONG G C, WANG S Y, et al. Application of microscopic CT in oil-gas reservoir study. Resources & Industries, 2013, 15(6): 98-102. DOI:10.3969/j.issn.1673-2464.2013.06.015 |
| [12] |
吕菲, 刘建立, 张佳宝, 等. 利用随机网络模型和CT数字图像预测近饱和土壤水分特征曲线. 灌溉排水学报, 2009(6): 18-21. LYU F, LIU J L, ZHANG J B, et al. Prediction of near saturated soil water retention curve using CT images and random network model. Journal of Irrigation and Drainage, 2009(6): 18-21. |
| [13] |
潘少伟, 梁鸿军, 李良, 等. 微观剩余油仿真研究进展. 岩性油气藏, 2013, 25(1): 16-20. PAN S W, LIANG H J, LI L, et al. Research progress on simulation of microscopic remaining oil. Lithologic Reservoirs, 2013, 25(1): 16-20. DOI:10.3969/j.issn.1673-8926.2013.01.004 |
| [14] |
尚庆华, 王玉霞, 黄春霞, 等. 致密砂岩油藏超临界与非超临界CO2驱油特征. 岩性油气藏, 2018, 30(3): 153-158. SHANG Q H, WANG Y X, HUANG C X, et al. Supercritical and non-supercritical CO2 flooding characteristics in tight sandstone reservoir. Lithologic Reservoirs, 2018, 30(3): 153-158. |
| [15] |
车洪昌, 任耀宇, 刘汉平, 等. 龙虎泡油田活性水驱油室内实验研究. 岩性油气藏, 2011, 23(2): 128-132. CHE H C, REN Y Y, LIU H P, et al. Laboratory study on oil displacement with active water in Longhupao Oilfield. Lithologic Reservoirs, 2011, 23(2): 128-132. DOI:10.3969/j.issn.1673-8926.2011.02.025 |
| [16] |
LI R, YANG Y S, PAN J, et al. Lattice Boltzmann modeling of permeability in porous materials with partially percolating voxels. Physical Review E, 2014, 90(3): 1-10. |
| [17] |
HUO F, HONG Z, HUO F, et al. MRT-LBM-based numerical simulation of seepage flow through fractal fracture networks. Science China-Technological Sciences, 2013, 56(12): 3115-3122. DOI:10.1007/s11431-013-5402-3 |
| [18] |
凌跃胜, 张雄, 陈子痛, 等. 有限差分法的一种边界处理方法. 电机与控制学报, 1999(1): 18-21. LING Y S, ZHANG X, CHEN Z T, et al. A new method of boundary treatment in finite difference method. Electric Machines and Control, 1999(1): 18-21. |
| [19] |
黄朝琴, 高博, 王月英, 等. 基于模拟有限差分法的离散裂缝模型两相流动模拟. 中国石油大学学报(自然科学版), 2014, 38(6): 97-105. HUANG Z Q, GAO B, WANG Y Y, et al. Two-phase flow simulation of discrete fracture model using a novel mimetic finite difference method. Journal of China University of Petroleum (Edition of Natural Science), 2014, 38(6): 97-105. DOI:10.3969/j.issn.1673-5005.2014.06.015 |
| [20] |
ZHANG Q G, JU Y, GONG W B, et al. Numerical simulations of seepage flow in rough single rock fractures. Petroleum, 2015, 1(3): 200-205. DOI:10.1016/j.petlm.2015.09.003 |
| [21] |
AFIF M, AMAZIANE B. Numerical simulation of two-phase flow through heterogeneous porous media. Numerical Algorithms, 2003, 34(2/4): 117-125. |
| [22] |
高辉, 孙卫, 路勇, 等. 特低渗透砂岩储层油水微观渗流通道与驱替特征实验研究:以鄂尔多斯盆地延长组为例. 油气地质与采收率, 2011, 18(1): 58-62. GAO H, SUN W, LU Y, et al. Experimental study on microflow channel and displacement characteristic of ultra-low permeability sandstone reservoir:Taking Yanchang formation, Ordos Basin as example. Petroleum Geology and Recovery Efficiency, 2011, 18(1): 58-62. DOI:10.3969/j.issn.1009-9603.2011.01.017 |
| [23] |
刘太勋, 徐怀民. 扇三角洲储层微观剩余油分布模拟试验. 中国石油大学学报(自然科学版), 2011, 35(4): 20-26. LIU T X, XU H M. Micro-remaining oil distribution simulation test of fan delta reservoir. Journal of China University of Petroleum(Edition of Natural Science), 2011, 35(4): 20-26. DOI:10.3969/j.issn.1673-5005.2011.04.004 |
| [24] |
SIRIVITHAYAPAKORN S, KELLER A. Transport of colloids in saturated porous media:a pore-scale observation of the size exclusion effect and colloid acceleration. Water Resources Research, 2003, 39(4): 1-11. |
| [25] |
LI G, REN W, MENG Y, et al. Micro-flow kinetics research on water invasion in tight sandstone reservoirs. Journal of Natural Gas Science & Engineering, 2014, 20(2): 184-191. |
| [26] |
付晓燕, 孙卫. 低渗透储集层微观水驱油机理研究:以西峰油田庄19井区长82储集层为例. 新疆石油地质, 2005, 26(6): 681-683. FU X Y, SUN W. Study of microscopic oil-water displacement mechanism of low-permeability reservoir:an example of Chang 82 reservoir in well block Zhuang-19 of Xifeng Oilfield. Xinjiang Petroleum Geology, 2005, 26(6): 681-683. DOI:10.3969/j.issn.1001-3873.2005.06.023 |
| [27] |
朱华银, 周娟, 万玉金, 等. 多孔介质中气水渗流的微观机理研究. 石油实验地质, 2004, 26(6): 571-573. ZHU H Y, ZHOU J, WAN Y J, et al. Microscopic mechanism study of gas-water flow in porous media. Petroleum Geology & Experiment, 2004, 26(6): 571-573. DOI:10.3969/j.issn.1001-6112.2004.06.012 |
| [28] |
李登伟, 张烈辉, 周克明, 等. 可视化微观孔隙模型中气水两相渗流机理. 中国石油大学学报(自然科学版), 2008, 32(3): 80-83. LI D W, ZHANG L H, ZHOU K M, et al. Gas-water two-phase flow mechanism in visual microscopic pore mode. Journal of China University of Petroleum(Edition of Natural Science), 2008, 32(3): 80-83. DOI:10.3321/j.issn:1673-5005.2008.03.017 |
| [29] |
解伟, 赵蕾, 孙卫, 等. 利用微观水驱油模型实验对储层进行流动单元的划分. 吉林大学学报, 2008, 38(5): 745-748. XIE W, ZHAO L, SUN W, et al. Flow units division by applying the sandstone micro-model for water-oil displacement experiment. Journal of Jilin University, 2008, 38(5): 745-748. |
| [30] |
滕起, 杨正明, 刘学伟, 等. 特低渗透油藏平面渗流规律研究. 科学技术与工程, 2012, 12(36): 9820-9823. TENG Q, YANG Z M, LIU X W, et al. 2 D seepage flow rule of the ultra-low permeability reservoir. Science and Technology and Engineering, 2012, 12(36): 9820-9823. DOI:10.3969/j.issn.1671-1815.2012.36.005 |
| [31] |
周克明, 李宁, 张清秀, 等. 气水两相渗流及封闭气的形成机理实验研究. 天然气工业, 2002, 22(增刊1): 122-125. ZHOU K M, LI N, ZHANG Q X, et al. Experimental research on gas-water two-phase flow and confined gas formation mechanism. Natural Gas Industry, 2002, 22(Suppl 1): 122-125. |
| [32] |
吴聃, 鞠斌山, 陈常红, 等. 基于微观驱替实验的剩余油表征方法研究. 中国科技论文, 2015, 10(23): 2707-2710. WU D, JU B S, CHEN C H, et al. Research on characterization method for residual oil based on microscopic displacement experiments. China Science Paper, 2015, 10(23): 2707-2710. DOI:10.3969/j.issn.2095-2783.2015.23.003 |
| [33] |
秦梓钧, 刘保君, 张雪, 等. COMSOL Multiphysics有限元软件数值模拟气液两相流的可行性研究. 当代化工, 2016, 45(5): 916-919. QIN Z J, LIU B J, ZHANG X, et al. Study on the feasibility of using COMSOL multiphysics software in numerical simulation of gas-liquid two-phase flow. Contemporary Chemical Industry, 2016, 45(5): 916-919. DOI:10.3969/j.issn.1671-0460.2016.05.015 |
| [34] |
BATTAGLIA L, STORTI M A, D'ELIA J. Simulation of freesurface flows by a finite element interface capturing technique. International Journal of Computational Fluid Dynamics, 2010, 24(3/4): 121-133. |
| [35] |
MIAO X, NARSILIO G A, WU A, et al. A 3D dual pore-system leaching model. Part 1:Study on fluid flow. Hydrometallurgy, 2017, 167: 173-182. DOI:10.1016/j.hydromet.2016.11.015 |
| [36] |
TOURÉ M K, FAHSI A, SOULAÏMANI A. Stabilised finiteelement methods for solving the level set equation with mass conservation. International Journal of Computational Fluid Dynamics, 2016, 30(1): 18. |
| [37] |
TAN H, PILLAI K M. Finite element implementation of stressjump and stress-continuity conditions at porous-medium, clearfluid interface. Compute and Fluids, 2009, 38(6): 1118-1131. DOI:10.1016/j.compfluid.2008.11.006 |
| [38] |
刘文, 赵士友, 杜军军.微观驱油玻璃模型的制备方法: 中国, 2015108687064. 2018-05-22. LIU W, ZHAO S Y, DU J J. Preparation method of microscopic oil displacement glass model: China, 2015108687064. 2018-0522. |
| [39] |
宋考平, 张继成, 张涛, 等.微观驱油用三维玻璃多孔介质模型及制造方法: 中国, 2007100983281. 2010-12-22. SONG K P, ZHANG J C, ZHANG T, et al. Three-dimensional glass porous medium model for micro-oil displacement and its manufacturing method: China, 2007100983281. 2010-12-22. |
| [40] |
施雷庭, 王瑶, 舒政, 等.一种微观可视化刻蚀低渗透模型制备方法: 中国, 2016109536297. 2017-01-18. SHI L T, WANG Y, SHU Z, et al. A preparing method for microscopic visual etching low permeability model: China, 2016109536297. 2017-01-18. |
| [41] |
皮彦夫. 特殊人造岩心的类型及应用. 科学技术与工程, 2010, 10(31): 7724-7727. PI Y F. Type and application of special artificial cores. Science Technology and Engineering, 2010, 10(31): 7724-7727. DOI:10.3969/j.issn.1671-1815.2010.31.024 |
| [42] |
李芳芳, 杨胜来, 高旺来, 等. 大尺寸石英砂环氧树脂胶结人造岩心制备技术研究及应用. 科学技术与工程, 2013, 13(3): 685-689. LI F F, YANG S L, GAO W L, et al. The manufacturing methods of large size heterogeneity artificial cores. Science and Technology and Engineering, 2013, 13(3): 685-689. DOI:10.3969/j.issn.1671-1815.2013.03.029 |