Equivalent Capillary-Based Liquid Phase Invasion Model for Low Permeability Gas Reservoirs
-
摘要:
为了从微观尺度上探讨低渗透气藏液相侵入微观流动机理，采用激光刻蚀技术建立了致密砂岩孔隙网络模型，开展了液相侵入微观可视化流动实验，分析了液相侵入和流体返排过程孔隙网络内水相的动态分布；建立了基于等效毛细管束的低渗透气藏液相侵入微观流动模型，采用致密砂岩水相自吸侵入实验验证了模型的可行性。实验发现，孔隙网络内水相侵入与毛细管力侵入规律类似，初期液相主要沿着较大孔隙流动，液相通过与孔隙连通的喉道逐渐推进；较小喉道中的水相难以返排，阻碍气相的流动。研究结果表明，黏滞阻力对致密砂岩水相侵入起主导作用，但在负压差条件下，液相仍能侵入岩心，且岩石越致密，水相最大侵入深度越大。建立的液相侵入模型为低渗透气藏液相侵入损害及保护机理研究提供了理论参考。
Abstract:The goal of this study was to avoid formation damage by finding a better way to map the flow of fluids through pore networks in tight sandstones. In order to investigate the flow mechanism of liquid phase invasion in low permeability gas reservoirs from the microscopic scale, a pore network model of tight sandstone was established by using laser etching technology. In that way, the microscopic visualization flow experiment of liquid phase invasion was carried out, and the dynamic aqueous phase distribution in the pore network during the process and fluid flow back was analyzed. The liquid phase invasion microscopic flow model for low permeability gas reservoirs was established based on equivalent capillary beam, and the model was verified by aqueous phase self-absorption invasion experiment in tight sandstone. Experimental results showed that the rule of aqueous phase invasion in the pore network is similar to that of capillary force invasion. The liquid phase mainly flows through larger pores initially, and then advances through the throats communicating with the pores gradually; It is difficult for the aqueous phase in the smaller throat to flow back, which can hinder the flow of gas phase. The research suggested that viscous drag plays a dominant role in the invasion of aqueous phase in tight sandstone; the liquid phase can still invade the core under negative pressure difference. Further, the denser the rock, the greater the maximum invasion depth of aqueous phase would be. The established liquid phase invasion model will provide a theoretical reference in studying liquid phase invasion damage and protection mechanism of low permeability gas reservoirs.
-
Keywords:
- low permeability gas reservoirs /
- pore network /
- liquid invasion /
- flowmodel /
- microscopic mechanism
-
致密砂岩气藏具有孔渗特性差、水润湿性强、微裂缝发育和毛细管效应显著等特征，在钻井完井及储层改造过程中易遭受水锁损害，大幅度降低油气井产能[1–2]。针对水锁损害，众多学者开展了一系列实验与理论研究，建立了水锁预测模型[3–5]，探讨了水锁损害机理及防治措施[6–10]。赖南君等人[11]利用可对岩心加温加压的水锁损害评价仪研究了毛细管力与正压差对水锁损害的影响；刘建坤等人[12]将低磁场核磁共振T2谱技术与常规流动实验相结合，提出了水锁伤害核磁共振试验评价方法；唐洪明等人[13]利用核磁共振和T2谱技术开展了毛细管自吸与水相返排可视化实验；丁绍卿等人[14]将核磁共振技术应用于压裂液伤害机理研究，分析了黏土吸水效应及水锁效应对岩心渗透率的伤害程度。
以上方法多采用岩心驱替实验进行宏观规律研究，基于微观可视化的研究较少。为此，笔者建立了致密砂岩气藏孔隙网络模型，以期通过室内微观可视化实验分析液相侵入过程中孔隙网络内水相前缘的动态分布，揭示液相侵入过程中的微观流动机理；并在此基础上，建立了基于等效毛细管的低渗透气藏液相侵入微观流动模型，结合实验验证了模型的可行性，以期为低渗透气藏储层保护提供理论支撑。
1. 孔隙网络微观可视化实验
1.1 微观流动实验装置
该实验中，将HC区块须家河组致密砂岩制成铸体薄片，利用图形扫描软件刻画孔隙网格，然后采用激光刻蚀致密砂岩孔隙网格。微观流动实验装置主要由计量管、微观流动实验装置、体视显微镜和量筒组成（见图1）。所选岩心孔隙度为9.23%，渗透率为0.27 mD。刻蚀材料为光学石英玻璃，尺寸30 mm×30 mm，刻蚀模型尺寸11 mm×8 mm。石英玻璃无涂层，孔隙流道亲水，且石英玻璃透光性好，便于采用光学显微镜观察水相的侵入与返排。实验流体为蒸馏水，由于孔隙网络微观模型尺寸小，实验开始时在刻蚀模型注入端预先注入10 cm高的水柱，由于水相侵入开始后不再补充水，随着实验的进行水柱高度下降。
1.2 微观流动实验结果
不同侵入时间下孔隙水相分布实验结果如图2所示。
从图2可以看出，侵入初期（<8 s），在毛细管力和水柱压力作用下，水相侵入速度较快，迅速占据孔隙空间及其喉道；侵入中期，水相侵入速度减小，侵入深度缓慢增加，直至水侵前缘抵达模型右端出口；侵入后期，由于水柱压力降低且黏滞阻力增加，孔隙内水相流动能力下降，直至水侵呈稳定状态，水相侵入基本停止，侵入水相大多以残余水状态分布于孔隙网络模型中。由30 s时的水相分布可知，水相大部分位于孔喉处。根据实验结果，孔隙内水相侵入主要发生在流道中，孔道连接处的影响几乎可以忽略；同时，水相前缘推进与毛细管流动规律类似，可以为建立水相侵入模型提供实验依据。
水侵实验结束后，采用注射器返排侵入水相，当返排压力大于水相黏性阻力时，水相开始返排。返排后孔隙中的水相分布如图3所示。从图3可以看出，大部分孔隙中的水相得以返排，但仍有部分较小喉道中的水相未能返排。未能返排的水相以残余水状态存在，集中于孔喉处，阻碍气相流动。对于致密砂岩气藏，工作液与储层接触后，在井筒压差和毛细管力作用下沿孔隙侵入基质，由于致密砂岩通常具有亚束缚水饱和度特征，侵入水相可能在部分孔道壁面形成滞留水，即使后期采取负压差返排，孔喉处的滞留水仍不容易排出。
![]()
图 3 水相返排后孔隙空间水相分布Figure 3. Aqueous phase distribution in pore space after the flow back of aqueous phase2. 不同阶段水相侵入评价模型
由于水相侵入主要发生在孔隙流道中，孔喉处可以忽略，因此将致密砂岩孔隙网络简化为等效毛细管，忽略孔道壁面的摩擦系数。水相受到毛细管力、孔道壁面黏滞阻力、重力及惯性力作用[15]，根据主要作用力将水相侵入划分为初始阶段和平稳阶段，其中初始阶段包括惯性力侵入阶段、惯性力–黏滞阻力作用阶段及黏滞阻力作用阶段[16]。
2.1 惯性力侵入阶段
水相在惯性力侵入阶段开始进入毛细管，侵入量和侵入距离极小，主要作用力为毛细管力和惯性力，该阶段的瞬时侵入深度为[17]：
l1=t√2σcosθρr+Fpρπr2 (1) 式中：l1为惯性力作用阶段水相侵入深度，m；t为时间，s；σ为表面张力，N/m；θ为接触角，（°）；ρ为流体密度，kg/m3；r为孔隙喉道半径，m；Fp为压差作用力，N。
等效毛细管半径采用Kozeny-Carman公式计算：
r=√8Kϕ (2) 式中：K为渗透率，mD；ϕ为孔隙度。
惯性力作用阶段持续的时间为[17]：
t1=0.0232ρr2μ (3) 式中：t1为惯性力作用阶段持续时间，s。
2.2 惯性力–黏滞阻力作用阶段
惯性力–黏滞阻力作用阶段的黏滞阻力不能忽略，此时水相侵入深度为[18]：
l2=√2ba[t−1a(1−e−at)] (4) 其中a=8μρr2 (5) b=2σcosθρr (6) 式中：l2为惯性力–黏滞阻力作用阶段的水相侵入深度，m。
该阶段的作用时间为[17]：
t2=2.1151ρr2μ (7) 式中：t2为惯性力–黏滞阻力作用持续时间，s。
2.3 黏滞阻力作用阶段
水相侵入进入黏滞阻力作用阶段时，惯性效应可以忽略，动力平衡条件为[15]：
(2σcosθr−Fpπr2)πr2=8π μl3v3+ρgl3πr2sinα (8) 式中：α为自吸方向与水平方向的夹角，（°）；l3为平稳阶段水相侵入深度，m；
v3 为平稳阶段水相侵入速度，m/s；g为重力加速度，m/s2。由式（8）可得水相侵入速度的计算公式：
v3=rσcosθ4μl3−Fp8πμl3−ρgr2sinα8μ (9) 一维线性流动时，水相侵入过程中的压力分布为：
p=pin−(pin−pou)l3L (10) 式中：L为水侵压力波及长度，m；pin为入口处的流体压力，Pa；pou为出口处的流体压力，Pa。
则压差作用力为：
Fp=(pou−pin)l3πr2L (11) 将式（10）和式（11）代入式（9）求解。由于通常无法用解析方法求解，于是采用数值求解方法来求取近似解，笔者采用四阶Runge-Kutta差分格式求解：
{v3=cl3−dc=rσcosθ4μd=ρgr2sinα8μ+(pou−pin)r28μL (12) 式（12）即为建立的平稳阶段侵入模型，l3的初始值取t2时刻的l2。负压差条件下，水相侵入达到动力学平衡时的侵入深度可以认为是最大侵入深度：
l3max (13) 式中：l3max为负压差下水相最大侵入深度，m。
3. 模型验证
通过微观流动实验装置测量刻蚀模型毛细管束平均半径约为149.6 μm，表面张力取0.072 N/m，接触角取30°，黏度为1.0 mPa∙s，计算毛细管束不同水相侵入阶段水相自吸侵入深度，结果如图4所示。由式（2）和式（4）可知：惯性力和惯性力–黏滞阻力作用阶段持续时间均与毛细管半径成正比。由图4可知，惯性力作用阶段持续时间t1为0.52 ms时，水相自吸侵入深度为0.43 mm；惯性力–黏滞阻力作用阶段持续时间t2为47.0 ms时，水相自吸侵入深度为14.41 mm。之后水相侵入进入黏滞阻力作用阶段，水平方向自吸侵入不受重力的影响，侵入深度随时间增长不断增大；对于垂向自吸侵入，在重力作用下存在最大水相自吸侵入高度。算例条件下垂向水相最大自吸侵入高度为84.96 mm。致密砂岩孔喉半径通常介于0.03～2.00 μm[18]，可知惯性力和惯性力–黏滞阻力作用阶段持续时间短，水相自吸侵入深度主要取决于黏滞阻力作用阶段。
![]()
图 4 不同阶段毛细管水相自吸侵入深度Figure 4. Self-absorption invasion depth of capillary aqueous phase at different stages为了验证侵入模型（式（12））的可行性，采用西南石油大学研制的岩心自吸水测量仪，测试了致密砂岩垂向自吸侵入高度，实验装置及测试流程见文献[19]。岩心取自HC地区须家河组致密砂岩，渗透率0.24 mD，孔隙度9.3%，长度50.0 mm。实验结果表明，初始阶段水相侵入速度较快，之后侵入高度趋于平缓，实验数据与模型计算结果吻合度较高（见图5）。
![]()
图 5 模型预测结果与实验数据对比Figure 5. Comparison of model prediction results with experimental data通过式（13）可以计算负压差条件下最大水相侵入深度，此处采用负压差水相自吸实验数据进行验证。实验岩样取自川西蓬莱镇组和川中须家河组，岩样L54渗透率为0.17 mD，孔隙度为6.2%；岩样HE6渗透率为1.88 mD，孔隙度为11.8%[19]，实验欠压值1 MPa，实验结果与计算结果如图6所示。
从图6可以看出，负压差条件下，实验初期仍有水相侵入岩心，随着实验时间增长，水相侵入深度未明显增加。岩样L54和HE6的最大侵入深度的实验测试结果分别为3.1 和1.8 mm，计算结果分别为2.9和2.3 mm。可见，在相同实验条件下，岩石越致密，渗透率越低，水相最大侵入深度越大。
![]()
图 6 负压差条件下水相最大侵入深度Figure 6. The maximum invasion depth of aqueous phase under negative pressure difference4. 结 论
1）采用激光刻蚀技术雕刻了致密砂岩孔隙网络，开展了液相侵入微观可视化实验，分析了液相侵入过程中及返排后孔隙网络内的水相分布。实验结果表明，孔隙介质内的水相侵入主要发生在流道中，水相推进与毛细管流动规律类似，部分较小喉道中的水相不能返排。
2）孔隙网络内水相侵入过程包括惯性力作用、惯性力–黏滞阻力作用和黏滞阻力作用3个阶段。致密砂岩水相侵入的惯性力和惯性力–黏滞阻力作用阶段持续时间短，黏滞阻力对水相自吸侵入起主导作用。
3）将孔隙网络视为等效毛细管束，建立了致密砂岩液相侵入微观动力学模型，通过致密砂岩垂向自吸和负压差条件水平向侵入实验验证了模型的可行性。分析表明，负压差条件下水相仍能侵入岩心，且岩石越致密最大侵入深度越大。
-
![]()
图 3 水相返排后孔隙空间水相分布
Figure 3. Aqueous phase distribution in pore space after the flow back of aqueous phase
![]()
图 4 不同阶段毛细管水相自吸侵入深度
Figure 4. Self-absorption invasion depth of capillary aqueous phase at different stages
![]()
图 5 模型预测结果与实验数据对比
Figure 5. Comparison of model prediction results with experimental data
-
[1] 李皋, 蔡武强, 孟英峰, 等. 不同钻井方式对致密砂岩储层损害评价实验[J]. 天然气工业, 2017, 37(2): 69–76 http://d.old.wanfangdata.com.cn/Periodical/trqgy201702009 LI Gao, CAI Wuqiang, MENG Yingfeng, et al. Experimental evaluation on the damages of different drilling modes to tight sandstone reservoirs[J]. Natural Gas Industry, 2017, 37(2): 69–76 http://d.old.wanfangdata.com.cn/Periodical/trqgy201702009
[2] 杨旭, 孟英峰, 李皋, 等. 考虑水锁损害的致密砂岩气藏产能分析[J]. 天然气地球科学, 2017, 28(5): 812–818 http://d.old.wanfangdata.com.cn/Periodical/trqdqkx201705016 YANG Xu, MENG Yingfeng, LI Gao, et al. Productivity analysis of tight sandstone gas reservoirs considering water blocking dama-ge[J]. Natural Gas Geoscience, 2017, 28(5): 812–818 http://d.old.wanfangdata.com.cn/Periodical/trqdqkx201705016
[3] 蒋官澄, 王晓军, 关键, 等. 低渗特低渗储层水锁损害定量预测方法[J]. 石油钻探技术, 2012, 40(1): 69–73 doi: 10.3969/j.issn.1001-0890.2012.01.014 JIANG Guancheng, WANG Xiaojun, GUAN Jian, et al. The quantitative prediction method of water blocking damage in low and extra-low permeability reservoir[J]. Petroleum Drilling Techniques, 2012, 40(1): 69–73 doi: 10.3969/j.issn.1001-0890.2012.01.014
[4] 蒋官澄, 张弘, 吴晓波, 等. 致密砂岩气藏润湿性对液相圈闭损害的影响[J]. 石油钻采工艺, 2014, 36(6): 50–54 http://d.old.wanfangdata.com.cn/Periodical/syzcgy201406016 JIANG Guancheng, ZHANG Hong, WU Xiaobo, et al. Effect of tight sandstone gas reservoir wettability on liquid traps damage[J]. Oil Drilling & Production Technology, 2014, 36(6): 50–54 http://d.old.wanfangdata.com.cn/Periodical/syzcgy201406016
[5] 张益, 李军刚, 佟晓华, 等. 基于神经网络信息融合技术预测气藏水锁[J]. 特种油气藏, 2011, 18(2): 102–103, 110 doi: 10.3969/j.issn.1006-6535.2011.02.029 ZHANG Yi, LI Jungang, TONG Xiaohua, et al. Prediction of water lock in gas reservoirs based on neural network information fusion[J]. Special Oil & Gas Reservoirs, 2011, 18(2): 102–103, 110 doi: 10.3969/j.issn.1006-6535.2011.02.029
[6] 韦青, 李治平, 白瑞婷, 等. 微观孔隙结构对致密砂岩渗吸影响的试验研究[J]. 石油钻探技术, 2016, 44(5): 109–116 http://d.old.wanfangdata.com.cn/Periodical/syztjs201605019 WEI Qing, LI Zhiping, BAI Ruiting, et al. An experimental study on the effect of microscopic pore structure on spontaneous imbibition in tight sandstones[J]. Petroleum Drilling Techniques, 2016, 44(5): 109–116 http://d.old.wanfangdata.com.cn/Periodical/syztjs201605019
[7] LI Gao, REN Wenxi, MENG Yingfeng, et al. Micro-flow kinetics research on water invasion in tight sandstone reservoirs[J]. Journal of Natural Gas Science and Engineering, 2014, 20: 184–191 doi: 10.1016/j.jngse.2014.06.024
[8] 赵宏波, 贾进孝, 孟令涛, 等. 一种新型降水锁洗井液NDF–1的性能评价及现场试验[J]. 石油钻探技术, 2015, 43(6): 87–92 http://d.old.wanfangdata.com.cn/Periodical/syztjs201506016 ZHAO Hongbo, JIA Jinxiao, MENG Lingtao, et al. Performance evaluation and field application of a novel water lock reducing flushing fluid DNF-1[J]. Petroleum Drilling Techniques, 2015, 43(6): 87–92 http://d.old.wanfangdata.com.cn/Periodical/syztjs201506016
[9] 唐洪明, 徐诗雨, 王茜, 等. 克拉苏气田超致密砂岩气储层水锁损害[J]. 断块油气田, 2017, 24(4): 541–545 http://d.old.wanfangdata.com.cn/Periodical/dkyqt201704023 TANG Hongming, XU Shiyu, WANG Xi, et al. Water blocking damage of hyper-tight sandstone gas reservoir in Kelasu Gas Field[J]. Fault-Block Oil & Gas Field, 2017, 24(4): 541–545 http://d.old.wanfangdata.com.cn/Periodical/dkyqt201704023
[10] 徐新丽. 东风港油田特低渗透油藏微观孔隙结构及渗流特征试验研究[J]. 石油钻探技术, 2017, 45(2): 96–100 http://d.old.wanfangdata.com.cn/Periodical/syztjs201702016 XU Xinli. Experimental study on micro-pore structure and seepage characteristics of ultra-low permeability reservoirs in the Dongfeng-gang Oilfield[J]. Petroleum Drilling Techniques, 2017, 45(2): 96–100 http://d.old.wanfangdata.com.cn/Periodical/syztjs201702016
[11] 赖南君, 叶仲斌, 刘向君, 等. 低渗透致密砂岩气藏水锁损害室内研究[J]. 天然气工业, 2005, 25(4): 125–127 doi: 10.3321/j.issn:1000-0976.2005.04.040 LAI Nanjun, YE Zhongbin, LIU Xiangjun, et al. In-house study on water locking damage of tight sand gas reservoirs with low permeability[J]. Natural Gas Industry, 2005, 25(4): 125–127 doi: 10.3321/j.issn:1000-0976.2005.04.040
[12] 刘建坤, 郭和坤, 李海波, 等. 低渗透储层水锁伤害机理核磁共振实验研究[J]. 西安石油大学学报(自然科学版), 2010, 25(5): 46–49, 53 doi: 10.3969/j.issn.1673-064X.2010.05.011 LIU Jiankun, GUO Hekun, LI Haibo, et al. Experimental study on the water-blocking damage mechanism of low permeability reservoir by nuclear magnetic resonance[J]. Journal of Xi′an Shiyou University (Natural Science Edition), 2010, 25(5): 46–49, 53 doi: 10.3969/j.issn.1673-064X.2010.05.011
[13] 唐洪明, 朱柏宇, 王茜, 等. 致密砂岩气层水锁机理及控制因素研究[J]. 中国科学: 技术科学, 2018, 48(5): 537–547 http://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201805008.htm TANG Hongming, ZHU Baiyu, WANG Xi, et al. Mechanism and control factors of water blocking in tight sandstone gas reservoir[J]. Scientia Sinica (Technologica), 2018, 48(5): 537–547 http://www.cnki.com.cn/Article/CJFDTOTAL-JEXK201805008.htm
[14] 丁绍卿, 郭和坤. 应用核磁共振技术研究压裂液伤害机理[J]. 钻井液与完井液, 2006, 23(3): 60–62 doi: 10.3969/j.issn.1001-5620.2006.03.018 DING Shaoqing, GUO Hekun. Research on the damage mechanism of fracturing fluids through nuclear magnetic resonance technology[J]. Drilling Fluid & Completion Fluid, 2006, 23(3): 60–62 doi: 10.3969/j.issn.1001-5620.2006.03.018
[15] 汪传磊. 致密砂岩孔隙介质内水相侵入微观流动机理研究[D]. 成都: 西南石油大学, 2012. WANG Chuanlei. The study of microscopic flow mechanism of water invading tight sandstone porous media[D]. Chengdu: Southwest Petroleum University, 2012.
[16] FRIES N, DREYER M. The transition from inertial to viscous flow in capillary rise[J]. Journal of Colloid and Interface Science, 2008, 327(1): 125–128 doi: 10.1016/j.jcis.2008.08.018
[17] QUÉRÉ D. Inertial capillarity[J]. Europhysics Letters, 1997, 39(5): 533–538 doi: 10.1209/epl/i1997-00389-2
[18] NELSON P H. Pore-throat sizes in sandstones, tight sandstones, and shales[J]. AAPG Bulletin, 2009, 93(3): 329–340 doi: 10.1306/10240808059
[19] 谢晓永. 致密砂岩气藏欠平衡钻井储层保护适应性评价体系研究[D]. 成都: 西南石油大学, 2010. XIE Xiaoyong. Study on formation protection suitability evaluation system of underbalanced drilling in tight sandstone gas reservoir[D]. Chengdu: Southwest Petroleum University, 2010.
-
期刊类型引用(3)
1. 袁媛，孟英峰，李皋，苏晓明. 致密砂岩储层毛管自吸微观分布特征. 油气地质与采收率. 2020(05): 71-78 . 
百度学术
2. 于浩雨，于明德，李洲，熊亮. 洛伊凹陷西南部边界大断裂发育特征及其对油气成藏的控制作用. 油气地质与采收率. 2020(05): 13-24 . 百度学术
3. 孟也，李相方，何敏侠，蒋明洁. 气泡卡断过程中的喉道液领形态与聚并模型. 断块油气田. 2019(05): 632-637 . 百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 11399
- HTML全文浏览量: 5783
- PDF下载量: 21
- 被引次数: 4
