首页期刊简介编委会编辑部过刊浏览征稿简则期刊信息下载中心欢迎订阅English
引 en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

CO2提高天然气采收率的重力效应与气藏压力影响模拟

  • 刘树阳1,2
  • ， 孙宝江1,2
  • ， 宋永臣3
  • ， 张 毅3

1. 非常规油气开发教育部重点实验室(中国石油大学(华东)),山东青岛 266580； 2. 中国石油大学(华东)石油工程学院, 山东青岛 266580； 3. 大连理工大学能源与动力学院,辽宁大连 116024；

Simulation on gravity effect and reservoir pressure influence analysis in CO2 enhanced gas recovery

  • LIU Shuyang1,2
  • ， SUN Baojiang1,2
  • ， SONG Yongchen3
  • ， ZHANG Yi3

1. Key Laboratory of Unconventional Oil & Gas Development (China University of Petroleum (East China)), Ministry of Education, Qingdao 266580 , China； 2. School of Petroleum Engineering in China University of Petroleum(East China), Qingdao 266580 , China； 3. School of Energy and Power Engineering in Dalian University of Technology, Dalian 116024 , China；

作者简介:

刘树阳(1990-),男,博士,研究方向为二氧化碳强化非常规油气开发。E-mail:shuyang_liu@126.com。

通讯作者:

孙宝江(1963-),男,教授,博士,博士生导师,长江学者特聘教授,国家“973"计划项目首席科学家,研究方向为海洋石油工程、多相流理论及应用、非常规油气开发。E-mail:sunbj1128@vip.126.com。

中图分类号:TQ28

文献标识码:A

文章编号:1673-5005(2020)03-0081-09

DOI:10.3969/j.issn.1673-5005.2020.03.009

参考文献 1
CONINCK H C,CARBO M C,MIKUNDA T,et al.Car-bon capture and storage in industrial applications:tech-nology synthesis report[R/OL].(2010-12-09)[2019-03-27].https://publications.tno.nl/publication/34632076/JT069c/o10043.pdf. 
查找原文
参考文献 2
孙宝江,王金堂,孙文超,等.非常规天然气储层超临界 CO2 压裂技术基础研究进展[J].中国石油大学学报(自然科学版),2019,43(5):82-91.SUN Baojiang,WANG Jintang,SUN Wenchao,et al.Advances in fundamental research of supercritical CO2 fracturing technology for unconventional natural gas reser-voirs[J].Journal of China University of Petroleum(Edi-tion of Natural Science),2019,43(5):82-91. 
查找原文
参考文献 3
任韶然,张莉,张亮.CO2 地质埋存:国外示范工程及其对中国的启示[J].中国石油大学学报(自然科学版),2010,34(1):93-98.REN Shaoran,ZHANG Li,ZHANG Liang.Geological storage of CO2:overseas demonstration projects and its implications to China[J].Journal of China University of Petroleum(Edition of Natural Science),2010,34(1):93-98. 
查找原文
参考文献 4
van der BURGT M J,CANTLE J,BOUTKAN V K.Car-bon dioxide disposal from coal-based IGCC＇s in depleted gas fields [J].Energy Conversion and Management,1992,33(5/6/7/8):603-610. 
查找原文
参考文献 5
OLDENBURG C M,BENSON S M.CO2 injection for en-hanced gas production and carbon sequestration [ R].SPE 74367-MS,2002. 
查找原文
参考文献 6
OLDENBURG C M,PRUESS K,BENSON S M.Process modeling of CO2 injection into natural gas reservoirs for carbon sequestration and enhanced gas recovery [J].En-ergy & Fuels,2001,15(2):293-298. 
查找原文
参考文献 7
OLDENBURG C M,STEVENS S H,BENSON S M.E-conomic feasibility of carbon sequestration with enhanced gas recovery(CSEGR)[J].Energy,2004,29(9/10):1413-1422. 
查找原文
参考文献 8
GANZER L,REITENBACH V,PUDLO D,et al.Exper-imental and numerical investigations on CO2 injection and enhanced gas recovery effects in Altmark gas field(Cen-tral Germany)[J].Acta Geotechnica,2013,9:39-47. 
查找原文
参考文献 9
HUSSEN C,AMIN R,MADDEN G,et al.Reservoir simulation for enhanced gas recovery:an economic evalu-ation [J].Journal of Natural Gas Science and Engineer-ing,2012,5:42-50. 
查找原文
参考文献 10
JIKICH S A,SMITH D H,SAMS W N,et al.En-hanced gas recovery(EGR)with carbon dioxide se-questration:a simulation study of effects of injection strategy and operational parameters [R].SPE 84813,2003. 
查找原文
参考文献 11
REBSCHER D,OLDENBURG C M.Sequestration of carbon dioxide with enhanced gas recovery-case study altmark,North German Basin [ R ].LBNL-59033(2005-12-10)[2019-03-27].https://www.escholar-ship.org/uc/item/7sj298mc. 
查找原文
参考文献 12
AL-HASAMI A,REN S,TOHIDI B.CO2 injection for enhanced gas recovery and geo-storage:reservoir simula-tion and economics[R].SPE 94129,2005. 
查找原文
参考文献 13
HOU Z,GOU Y,TARON J,et al.Thermo-hydro-me-chanical modeling of carbon dioxide injection for en-hanced gas-recovery(CO2-EGR):a benchmarking study for code comparison [J].Environmental Earth Sciences,2012,67(2):549-561. 
查找原文
参考文献 14
LUO F,XU R N,JIANG P X.Numerical investigation of the influence of vertical permeability heterogeneity in stratified formation and of injection/production well per-foration placement on CO2 geological storage with en-hanced CH4 recovery [J].Applied Energy,2013,102:1314-1323. 
查找原文
参考文献 15
REGAN M.A numerical investigation into the potential to enhance natural gas recovery in water-drive gas reser-voirs through the injection of CO2 [D].Adelaide:The University of Adelaide,2010. 
查找原文
参考文献 16
CLASS H,EBIGBO A,HELMIG R,et al.A bench-mark study on problems related to CO2 storage in geolog-ic formations [J].Computational Geosciences,2009,13(4):409-434. 
查找原文
参考文献 17
SEO J G,MAMORA D D.Experimental and simulation studies of sequestration of supercritical carbon dioxide in depleted gas reservoirs [J].Journal of Energy Re-sources Technology,2005,127:1-6. 
查找原文
参考文献 18
ROBINSON B A,VISWANATHAN H S,VALOCCHI A J.Efficient numerical techniques for modeling multicom-ponent ground-water transport based upon simultaneous solution of strongly coupled subsets of chemical compo-nents [J].Advances in Water Resources,2000,23(4):307-324. 
查找原文
参考文献 19
JIANG Y.Techniques for modeling complex reservoirs and advanced wells [D].Palo Alto:Stanford Universi-ty,2007. 
查找原文
参考文献 20
PRUESS K.The TOUGH codes—a family of simulation tools for multiphase flow and transport processes in per-meable media [J].Vadose Zone Journal,2004,3:738-746. 
查找原文
参考文献 21
PRUESS K,OLDENBURG C M,MORIDIS G J.TOUGH2 user＇s guide version 2.0[Z/OL].Rep LBNL-43134(1999-11-01)[2019-03-27].https://escholarship.org/uc/item/4df6700h. 
查找原文
参考文献 22
PATEL M J,MAY E F,JOHNS M L.High-fidelity res-ervoir simulations of enhanced gas recovery with super-critical CO2 [J].Energy,2016,111:548-559. 
查找原文
参考文献 23
BIAGI J,AGARWAL R,ZHANG Z.Simulation and optimization of enhanced gas recovery utilizing CO2 [J].Energy,2016,94:78-86. 
查找原文
参考文献 24
SPAN R,WAGNER W.A new equation of state for car-bon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa [J].Journal of Physical and Chemical Reference Data,1996,25(6):1509-1596.
查找原文
目录contents

    摘要

    利用 TOGUH2 开展均质天然气藏模型中实施 CO2 提高天然气采收率的模拟研究,分析重力效应作用机制和气藏压力对 CO2 运移和 CH4 采收率变化规律的影响,为 CO2 提高天然气采收率的工程应用提供指导。 模拟结果表明: 重力效应抑制 CO2 向上运移,减缓了 CO2 与天然气纵向方向的混相,间接促进水平方向运移,导致气藏内 CO2 水平运移速率显著大于纵向运移速率,但是 CH4 采收率曲线不仅受重力效应影响,也受 CO2 注入流速和气藏压力的综合影响;气藏内的 CH4 生产时间及最终采收率随气藏压力降低而增加,在枯竭状态的低气藏压力条件下实施 CO2 提高天然气采收率能够获得更好的 CH4 采收率,但生产时间相对较长;较高气藏压力条件下,CO2 为超临界状态,重力效应作用明显,CO2 快速沉降并优先驱替气藏底层 CH4 ,但是容易从生产井底部自下而上突破至生产井,提高天然气采收率效果不佳。

    Abstract

    In this study, TOGUH2 simulator was applied to simulate CO2 enhanced natural gas recovery (CO2 -EGR) with a homogeneous natural gas reservoir model, in which the gravity effect and influence of reservoir pressure on CO2 migration and CH4 recovery were studied in order to provide theoretical support for the application of the CO2 -EGR technique. The simula- tion results show that the gravity effect can inhibit upward migration of CO2 and slow down the vertical mixing of CO2 and CH4 , which can indirectly promote the horizontal migration and make the horizontal displacing rate significantly greater than the vertical one. However, the produced curve of CH4 is not only determined by gravity effect, but also affected by CO2 in- jection rate and reservoir pressure. The CH4 production time and final recovery efficiency increase with the decrease of reser-voir pressure, and better CH4 recovery efficiency can be achieved at the depleted reservoir pressure than other operation con- ditions, while the production duration can be relatively longer. When CO2 is in a supercritical state at relatively high reservoir pressure, the gravity effect can have relatively significant effect on the displacement process. CO2 can rapidly sink and pref- erentially displace CH4 in the bottom layer of the gas reservoir, but it can break through vertically from the bottom of the pro- duction well, resulting in a poor natural gas recovery.

  • CO2 提高天然气采收率(CO2 enhanced gas re-covery, CO2-EGR)技术既能提高天然气产量,又可以地质封存CO2,是一项潜力巨大的二氧化碳捕集、利用与封存(CCUS) 技术[1-3]。自20世纪90 年代CO2-EGR技术概念[4]提出至今,研究人员首先通过大量模拟研究,证明了CO2 -EGR技术及其经济可行性[5-9]。研究表明过早地注入CO2 有增大产气速率的效果[10],而储层非均质性会导致快速流道的形成,造成CO2 过早突破[4],地层水的存在,可以溶解CO2,缓解非均质性引起的CO2 过早突破不利影响[11-13]。井网布局对采收率的提高也发挥重要作用,增大CO2 注入井与CH4 产气井之间的距离,可以增加天然气产量[13],注入井和生产井射孔均布置在最低渗透率层中,有利于提高CO2 埋存量[14]。 此外,水驱型气藏中黏性重力比也是影响CO2 驱替天然气的重要因素[15]。上述研究探究了诸多因素的影响,但是CO2 驱替天然气过程中重力效应的作用机制尚未研究透彻;不同气藏压力条件下实施CO2-EGR技术的研究也较为缺乏。针对以上问题,笔者通过开展均质气藏中CO2-EGR模拟,剖析重力效应、CO2 注入时机对气藏内CO2 运移规律及CH4 采收率的影响,探究最佳注入压力,为CO2-EGR的实施提供理论参考与数据支撑。

  • 1 天然气藏模型与模拟方案

  • 1.1 气藏模型

  • 本研究选择文献[16]中的天然气藏模型作为模拟对象,模拟区域为方形的五点井网布置(201.19 m× 201.19 m×45.72 m)的1/4。 CO2 注入井位于五点井网气藏模型区域的中心,天然气(假定全部为CH4 )生产井位于4 个角,模型的边界均绝热,且垂直边界方向无流动,沿边界方向为对称流动,如图1 所示。 根据前人研究经验[5],CO2 注入气藏底部有利于CO2 -EGR,因此本研究将CO2 注入井设置于模型底层, CH4 生产井位于模型顶层,网格划分为44×44×10,共19360 个边长为4.572 m的正方体单元。模拟中,注入井以恒定质量流量将CO2 注入气藏底部,CH4 生产井维持压力不变进行开采。

  • 图1 气藏模型

  • Fig.1 Model of gas reservoir

  • 气藏模型为均质储层模型,岩性数据取自美国德克萨斯州北部的碳酸盐岩储层[17],假设不含水, 孔隙度为0.23,水平渗透率为5.0×10 -14 m 2,垂向渗透率为5.0×10-15 m 2 ,温度为66.7℃,CO2 与天然气之间扩散系数设定为6.0×10 -7 m 2·s-1。 在具体模拟研究中,生产井产气中含有CO2 质量分数超过20%时,认为天然气被CO2 大规模污染,不再具有开采价值,为生产井关井时刻。 为探究CO2 注入对CH4 采收率的影响,计算定义关井时刻的采收率 η 与产出率 θ 分别为

  • η=mCH4,prodmCH4,orig
    (1)

  • θ=mCH4,prodmCO2,inj
    (2)

  • 式中, mCH4,prod为关井时CH4 累积产量; mCH4,orig为CO2 注入前气藏内CH4 的量;mCO2,inj为关井时CO2 累积注入量。

  • 1.2 模拟软件选择与改进

  • 针对CO2 地质利用与封存,研究人员提出了许多模拟器, 如COORES、 DUMUX、 FEHM、 GPRS及TOUGH2 等[16,18-20],其中TOUGH2 被多个专业科研机构应用于CO2 地质封存数值模拟研究,如澳大利亚联邦科学与工业研究组织(CSIRO) 的CO2CRC、 法国地质调查局的BRGM及德国亚琛大学的E.ON ERC等项目。 因此本研究应用TOUGH2 开展CO2 提高天然气采收率过程的模拟。

  • TOUGH2 由美国劳伦斯伯克利国家实验室开发[20],是一款适用于研究地层多孔介质内多组分多相流动过程的通用型模拟器。 TOUGH2 中,采用积分有限差分法对空间进行离散,而时间则使用隐式的一阶有限差分法离散。 在TOUGH2 软件架构中, 针对不同的CO2 地质封存过程使用不同的流体属性模块进行模拟。 选择EOS7C [22]模块模拟器,该模块具有专门计算CO2—CH4 流体属性的热力学性质的部分,适用于CO2 提高天然气采收率模拟研究。

  • 1.3 控制方程

  • 1.3.1 质量守恒方程

  • 地下流体多组分—多相流体质量守恒方程[21] 为

  • ddt∫VnMκdVn=∫VnFκndΓn+∫VnqκdVn
    (3)

  • 式中,Vn 为任意的控制体积; Γn为控制体积Vn 的闭合表面;n为指向体积Vn 的表面元素dΓn的法向矢量; Mκ 为单位体积中的质量积累; Fκ 为控制体积的净质量通量; qκ 为控制体积内的质量源/汇。

  • 式(3)中任意组分的质量积累项[23]表示为

  • Mκ=φ∑βSβρβXκ,β
    (4)

  • 式中, φ 为多孔介质的孔隙度; Sβ和ρβ分别为 β 相的饱和度和密度; Xκ,β为 β 相中 κ 组分的质量分数。

  • 质量通量的广义形式也表示为各相的总和,

  • Fκ=∑βXκ,βFβ
    (5)

  • 其中

  • Fβ=ρβvβ=-kkrβρβμβ(∇pβ-ρβg),pβ=p+pcapβ

  • 式中, Fβ 为 β 相的质量通量; vβ、μβ 和 ρβ 分别为 β 相的达西流速、黏度和压力;k为绝对渗透率;krβ为对于 β 相的相对渗透率;g为自由落体加速度;p和pcap茁分别为参考相压力和毛细管压力。

  • 1.3.2 能量守恒方程

  • 地下流体多组分—多相流体能量守恒方程[21]为

  • ddt∫VnEκdVn=∫ΓnDκndΓn+∫VnrκdVn
    (6)

  • 式中, Eκ 为单位体积中的能量积累;Dκ  为控制体积的能量通量;rκ  为控制体积内的能量源/汇。

  • 式(6)中任意组分的能量积累项[23]表示为

  • Eκ=(1-φ)ρRCRT+φ∑βSβρβuβ
    (7)

  • 式中,ρR 和CR 分别为岩石颗粒密度和比热容;T为温度;uβ 为 β 相的比内能。

  • 热通量 Dκ 由导热和对流热流组成:

  • Dκ=-λ∇T+∑βhβDβ
    (8)

  • 式中,λ 为热导率; ∇T为温度梯度;hβ 和 Dκ 分别为 β 相的比焓和热通量。

  • 1.4 模拟可靠性验证

  • 文献[16]中的天然气藏储层模型被广泛作为基准模型用于CO2-EGR过程的模拟研究, CO2CRC [16]、Luo等[14] 和Patel等[22] 都应用基准模型进行过CO2-EGR模拟研究,本研究也采用该基准模型进行模拟,并与以上研究结果进行对比,验证本研究中CO2-EGR模拟的可靠性。

  • 图2 为CO2CRC [16]、Luo等[14]、Patel等[22] 和本研究对基准模型模拟的产气中CO2 与CH4 质量流

  • 图2 基准模拟的产气中CO2 与CH4 质量流量变化

  • Fig.2 CO2 and CH4 mass flow rates at production well in benchmark model

  • 量变化。 CO2CRC与本研究均应用TOUGH2 软件开展模拟研究, Luo等[14] 和Patel等[22] 分别应用FLUENT和COMSOL进行模拟研究。 结果比较表明,基准模型模拟中Luo等[14] 的CO2 突破最早而Patel等[22]的CO2 突破最晚。 本研究的模拟结果与CO2CRC基本相似,产气中CH4 的质量流量在初期略有不同,可能由于基准模型储层的迂曲度并未明确规定,导致结果稍有差异。

  • 1.5 模拟方案

  • 重力效应对CO2-EGR过程中驱替运移规律、 提高天然气采收率发挥着重要作用。 超临界CO2 与CH4 之间存在较大的密度、黏度差异,受重力作用,CO2-EGR中超临界CO2 会向气藏储层下部沉降,即在较轻的天然气下形成“垫气冶 [4],在垂直方向上减缓混合、延迟突破,进而增加CH4 的采收率, 但具体过程仍需深入探究。 通过是否考虑重力的对比模拟,探究重力效应对驱替运移的作用,并且分析不同气藏压力条件下注入CO2 对天然气采收率的影响,为CO2-EGR实施提供指导。

  • 开展CO2 提高天然气采收率模拟的详细参数如表1 所示。 其中根据Biagi等[23]模拟结果,3.55 MPa枯竭状态下均质气藏模型的最优注入流速为4.67×10 -3 m 3·s -1 ,本研究中除基准模拟及其对应的不考虑重力模拟,其他模拟均采用此体积注入速率。 体积注入速率和质量注入流速换算时用到的CO2 密度数据通过NIST中专门用于CO2 热力学性质预测的Span和Wagner状态方程[24] 计算得到。

  • 表1 气藏模型中CO2-EGR模拟参数

  • Table 1 Simulation cases of CO2-EGR in homogeneous gas reservoir model

  • 2 重力效应

  • 2.1 产气流量对比

  • 图3 为不同条件下是否考虑重力模拟的产气中CO2 与CH4 流量曲线对比。 由图3 可知:不考虑重力,模拟7 #参数组合条件下产气中CO2 出现较晚,并且CH4 流量在产气后期下降幅度也较小,此时的CH4 采收率达到63.3%,明显高于普通模拟1 #参数组合条件下的59.6%;相较于普通模拟2 #,不考虑重力模拟的8 #产气中CO2 出现更早,且上升幅度更大,同时,CH4 流量在产气后期下降也更早,此时8 #的CH4 采收率(61.3%)明显低于2 # (62.9%);压力为9.0 MPa、注入流速为4.67×10 -3 m 3·s -1条件下,CO2 为超临界态,相比于普通模拟4 #,不考虑重力模拟的9 #参数组合条件下的产气中CO2 出现较晚,产气后期CH4 流量曲线下降趋势也相应的推迟了,造成9 #的CH4 采收率(63.6%)大于4 #的采收率(59.0%)。

  • 图3 不同条件下重力对模拟产气中CO2 与CH4 质量流量的影响

  • Fig.3 Effect of CO2 and CH4 mass fluxes between normal case and no-gravity case

  • 3组对比模拟结果表明,气藏中是否考虑重力的CO2-EGR模拟产气中的CO2 和CH4 流量曲线之间并没有一致性变化规律,其不仅受重力效应作用, 同时也受CO2 注入速率和气藏压力的综合影响。CH4 采收率的变化规律则与产气中CO2 和CH4 流量的变化密切相关。

  • 此外,对比普通模拟1 #和2 #发现,相同气藏压力条件下,CO2 注入流速增大,CH4 流量在产气前期增长速度变慢。 对比普通模拟2 #和4 #,相同的体积注入流速条件下,高压(CO2 超临界) 时,CH4 流量在产气前期增长速度加快,产气速率提高,主要原因是超临界CO2 与CH4 之间的密度、黏度差异远大于气态CO2 与CH4 差异,CO2 沉降明显,形成“垫气冶, CH4 在浮力驱动下向上流动加快,因而产气前期CH4 质量流量增长速度加快。

  • 2.2 气藏内气体组分分布对比

  • 图4 为气藏压力3.55 MPa、不同CO2 注入速率条件下,考虑重力普通模拟(1 #和2 #)和不考虑重力模拟(7 #和8 #)的气藏内CO2 和CH4 三维分布。 其中,wCO2和wCH4分别为CO2 和CH4 的质量分数。

  • 图4 重力及注入流速对模拟气藏内CO2 与CH4 分布的影响

  • Fig.4 Effect of gravity and injection rate on CO2 and CH4 distribution in gas reservoir

  • 由图4(a)、(b)可以看出,注入流速1.59×10 -3 m 3·s -1条件下不考虑重力模拟7 #中CO2 在纵向扩散速度比普通模拟1 #更快,500 d时CO2 已经运移至注入井上部的顶层区域,而此时1 #的该区域仍为纯CH4,7 #中CO2 水平运移明显落后于1 #,且随着时间的推移,这种现象越发明显。 2000 d时,7 #内CO2 从气藏顶层流动至生产井实现突破,但生产井下部区域仍然由未产出的CH4 占据,且CO2 驱替过渡带前缘甚至未到达模拟区域的底层侧边界;1 #中CO2 从生产井底部向上运移至生产井实现突破,模拟区域内底层CH4 全部被CO2 驱替产出,但是顶层大量CH4 未能产出。 1 #和2 #模拟中CO2 驱替未波及到的顶层区域或生产井下部底层区域的体积影响最终采收率,前者的体积大于后者,所以模拟1 #的CH4 采收率小于7 #的。

  • 与图4(a)、(b)中1 #与7 #相似,图4(c)主要区别在于2 #顶层区域内被驱替产出的CH4 量更多,2 #的CH4 采收率大于8 #和1 #。 这说明较大的注入速率增大CO2 在气藏顶层的驱替波及区域面积,提高天然气采收率。

  • 图5和图6分别为气藏压力大于CO2 临界压力条件下是否考虑重力模拟4 #与9 #中气藏内CO2 与CH4 分布的三维视图和注入井所在的二维x-z纵向视图。与图4(c)、(d)中气藏压力3.55 MPa条件下模拟2 #与不考虑重力模拟8 #相比,图5 中普通模拟4 #中气藏底层的CH4 被超临界CO2 驱替产出的更加彻底,这主要是由于超临界CO2 密度、黏度与CH4 差异更大,受重力影响CO2 沉降至气藏底部的速度更快,抑制CO2 垂直向上运移的作用更大,促进CO2 优先驱替气藏底部CH4。图5中9 # 与图4(d)中8 #的表现极为类似,即在不考虑重力效应时, 不同压力条件下模拟的CO2 运移规律相似。图6二维x-z纵向视图中,重力分异的影响表现得更为明显。 800 d时,不考虑重力模拟9 #气藏底层CH4仍未完全被CO2 驱替,驱替过程中的混相过渡带呈近似球形弧面,与底层水平面呈近似直角,CO2 运移的纵向分量与水平分量相近。普通模拟4 #中,受重力效应作用,纵向运移受到抑制而远小于水平运移。400 d时,气藏底层CH4 已经完全被CO2 驱替,混相过渡带与底层水平面呈小于30°倾角。 此外,不考虑重力模拟9 # 的混相过渡带厚度均匀一致,且较厚;考虑重力模拟4 #的混相过渡带在左上顶部位置较薄,而在右下位置较厚,且明显薄于不考虑重力模拟9 #,说明重力效应起着抑制CO2 扩散运移、减小混相的作用。

  • 图59.0 MPa条件下是否考虑重力模拟的气藏内CO2 与CH4 三维分布

  • Fig.5 Three-dimensional CO2 and CH4 distribution in gas reservoir of normal case and no-gravity case at 9.0 MPa

  • 图69.0 MPa条件下是否考虑重力模拟的气藏内CO2 与CH4 分布的二维x-z纵向视图

  • Fig.6 Two-dimensional x-z portrait view of CO2 and CH4 distribution in gas reservoir of normal case and no-gravity case at 9.0 MPa

  • 综上所述,不考虑重力模拟中,CO2 向上运移速率增大,最终水平方向和纵向方向的运移速率相似, 驱替混相过渡带呈近似球形弧面;考虑重力的普通模拟中,CO2 与CH4 之间较大的密度和黏度差异造成了重力分异,在重力效应作用下CO2 向气藏底部沉降,抑制向上运移,减缓CO2 与天然气纵向混相, 气藏内CO2 纵向运移速率远小于水平运移速率;此外,重力效应间接促进CO2 在水平方向的运移,尤其在超临界条件下优先完成底层水平驱替,最终将驱替过程演变为类似垂向驱替。

  • 3 气藏压力影响

  • 探讨气田不同开采时期( 即不同气藏储层压力)注入CO2 对天然气采收率的影响,开展天然气生产阶段(气藏压力为18、15、9 和6 MPa)及气田枯竭阶段(气藏压力为3.55 MPa) CO2 提高天然气采收率模拟研究,对应表1 中的2 #~6 #模拟。

  • 3.1 产气曲线与CH4 采收率变化规律

  • 图7 为不同气藏压力、相同CO2 体积注入流速条件下CO2-EGR模拟产气中CO2 与CH4 流量变化曲线。 随气藏压力升高,产气初期CH4 流量增长速度加快,流量达到稳定的时间减少,产气后期CH4 流量下降前移,与此同时CO2 突破时间也相应地提前。 最大体积流量随着气藏压力增大而增大,增大的幅度在3.55~9.0 MPa情况下不明显,但是15.0 和18.0 MPa情况下的最大体积流量要远大于低压情况。

  • 图7 不同气藏压力模拟的产气中CO2 与CH4 流量变化

  • Fig.7 CO2 and CH4 volume fluxes for cases at different reservoir pressures

  • 图8 为生产井关井时间、采收率及产出率随气藏压力的变化。可以发现,虽然CO2 体积注入速率相同,但是随着气藏压力增大,生产时间变短,关井时刻提前;CH4 采收率和产出率与生产井关井时间变化规律一致,随气藏压力增大而降低。采收效果最优的气藏压力为枯竭阶段的3.55 MPa,对应的CH4 采收率最高为62.9%。

  • 图8 关井时间、采收率及产出率随注入压力变化

  • Fig.8 Well shut-in time, recovery efficiency and production efficiency at different reservoir pressures

  • 3.2 气藏内CO2 和CH4 分布

  • 图9 和图10 分别为不同气藏压力模拟中的气藏内CO2 与CH4 分布的三维视图和二维纵向x-z前视图。

  • 图9 不同气藏压力条件下的气藏内CO2 与CH4 三维分布

  • Fig.9 Three-dimensional CO2 and CH4 distribution in gas reservoir of cases with different reservoir pressure

  • 图10 不同气藏压力条件下的气藏内CO2 与CH4 分布的二维x-z纵向视图

  • Fig.10 Two-dimensional x-z portrait view of CO2 and CH4 distribution in gas reservoir of cases with different reservoir pressures

  • 对于气藏压力为15 和18 MPa的模拟5 #和6 #, 生产井关井时间均小于600 d,所以本研究只分析其在注入CO2 后600 d的气藏内CO2 分布;其他模拟的生产井关井时间均为600~800 d,因此也对800 d时气藏内CO2 分布进行描述。 从图9 可以看出,体积注入速率相同的模拟2 #~6 #中,随着气藏压力下降,CO2 驱替CH4 到达生产井的速度变慢,并且模拟2 # ~6 #中的CO2 均从底部自下向上突破至生产井,顶层中部区域内的CH4 均未能被有效驱替产出。 在模拟4 #~6 #中,注入的CO2 为超临界状态,4 # 中气藏底层的CH4 在800 d时几乎全部被CO2 驱替产出,5 #和6 #中藏底层的CH4 在600 d内均被驱替产出;而气藏压力为3.55 和6.0 MPa的2 #和3 #中, CO2 为气态,直至到达关井时刻生产井下部区域仍然存在一些未被驱替的CH4。 如图10 所示,CO2 的水平运移距离显著大于纵向,气藏压力越高,该现象越明显。 生产井关井前,低气藏压力条件下CO2 与CH4 之间的驱替过渡带为倾斜面,而高气藏压力条件下则为近似水平面。 这表明CO2 超临界条件下, 优先完成气藏底层CH4 的水平驱替,然后形成近似垂直向上驱替过程,能够更加彻底地将气藏底层中的天然气驱替产出。 这主要是由于超临界CO2 的密度及黏度远大于气藏内天然气,重力效应明显,超临界CO2 向气藏底层沉降的速度更快,水平运移速率显著大于纵向运移速率,优先将气藏底层的CH4 驱替至生产井产出,最终形成形似自下向上的垂向驱替过程。 但是气藏压力较高时,CO2 沿气藏底层驱替CH4 从而更早地运移至生产井下部区域,受压力驱动自下向上突破至生产井的时间更早,使靠近注入井上部区域内CH4 未能被有效驱替,造成高气藏压力模拟的CH4 采收率小于低气藏压力模拟。

  • 综上所述,生产井关井时间、CH4 采收率和产出率均随气藏压力降低而增加,在气藏枯竭状态下实施CO2-EGR能够获得更好的CH4 采收率,但生产时间相对较长。 较高气藏压力下的CO2-EGR模拟中,注入气藏的CO2 为超临界状态,重力效应作用明显,CO2 快速向气藏底层沉降,优先完成气藏底层CH4 的水平驱替,最终形成形似自下向上的垂向驱替过程,但是针对本研究中的气藏区块,易从生产井底部自下而上突破至生产井,进而导致高压情况下CO2-EGR的天然气采收效率效果相对低压情况下略差。

  • 4 结论

  • (1)是否考虑重力CO2 -EGR模拟的产气中CO2 和CH4 流量曲线之间并没有一致性的变化规律,受CO2 注入速率及气藏压力的综合影响,而CH4 采收率变化与产气中CO2 和CH4 流量的变化规律密切相关。

  • (2)不考虑重力模拟中,CO2 纵向与水平运移速率相近,驱替混相过渡带呈近似球形弧面;考虑重力模拟中,重力效应抑制CO2 向上运移,减缓与天然气纵向混相,间接促进水平运移,CO2 水平运移速率显著大于纵向运移速率。

  • (3)生产井关井时间、CH4 采收率和产出率随气藏压力降低而增加,在气藏枯竭状态下实施CO2-EGR能够获得更好的CH4 采收率,但生产时间相对较长。

  • (4)较高气藏压力条件下实施CO2-EGR,CO2 为超临界状态,重力效应作用明显,CO2 快速沉降并优先水平驱替气藏底层CH4,但是容易从生产井底部自下而上突破至生产井,造成提高天然气采收率效果不佳。

  • 参考文献

    • [1] CONINCK H C,CARBO M C,MIKUNDA T,et al.Car-bon capture and storage in industrial applications:tech-nology synthesis report[R/OL].(2010-12-09)[2019-03-27].https://publications.tno.nl/publication/34632076/JT069c/o10043.pdf. 

    • [2] 孙宝江,王金堂,孙文超,等.非常规天然气储层超临界 CO2 压裂技术基础研究进展[J].中国石油大学学报(自然科学版),2019,43(5):82-91.SUN Baojiang,WANG Jintang,SUN Wenchao,et al.Advances in fundamental research of supercritical CO2 fracturing technology for unconventional natural gas reser-voirs[J].Journal of China University of Petroleum(Edi-tion of Natural Science),2019,43(5):82-91. 

    • [3] 任韶然,张莉,张亮.CO2 地质埋存:国外示范工程及其对中国的启示[J].中国石油大学学报(自然科学版),2010,34(1):93-98.REN Shaoran,ZHANG Li,ZHANG Liang.Geological storage of CO2:overseas demonstration projects and its implications to China[J].Journal of China University of Petroleum(Edition of Natural Science),2010,34(1):93-98. 

    • [4] van der BURGT M J,CANTLE J,BOUTKAN V K.Car-bon dioxide disposal from coal-based IGCC＇s in depleted gas fields [J].Energy Conversion and Management,1992,33(5/6/7/8):603-610. 

    • [5] OLDENBURG C M,BENSON S M.CO2 injection for en-hanced gas production and carbon sequestration [ R].SPE 74367-MS,2002. 

    • [6] OLDENBURG C M,PRUESS K,BENSON S M.Process modeling of CO2 injection into natural gas reservoirs for carbon sequestration and enhanced gas recovery [J].En-ergy & Fuels,2001,15(2):293-298. 

    • [7] OLDENBURG C M,STEVENS S H,BENSON S M.E-conomic feasibility of carbon sequestration with enhanced gas recovery(CSEGR)[J].Energy,2004,29(9/10):1413-1422. 

    • [8] GANZER L,REITENBACH V,PUDLO D,et al.Exper-imental and numerical investigations on CO2 injection and enhanced gas recovery effects in Altmark gas field(Cen-tral Germany)[J].Acta Geotechnica,2013,9:39-47. 

    • [9] HUSSEN C,AMIN R,MADDEN G,et al.Reservoir simulation for enhanced gas recovery:an economic evalu-ation [J].Journal of Natural Gas Science and Engineer-ing,2012,5:42-50. 

    • [10] JIKICH S A,SMITH D H,SAMS W N,et al.En-hanced gas recovery(EGR)with carbon dioxide se-questration:a simulation study of effects of injection strategy and operational parameters [R].SPE 84813,2003. 

    • [11] REBSCHER D,OLDENBURG C M.Sequestration of carbon dioxide with enhanced gas recovery-case study altmark,North German Basin [ R ].LBNL-59033(2005-12-10)[2019-03-27].https://www.escholar-ship.org/uc/item/7sj298mc. 

    • [12] AL-HASAMI A,REN S,TOHIDI B.CO2 injection for enhanced gas recovery and geo-storage:reservoir simula-tion and economics[R].SPE 94129,2005. 

    • [13] HOU Z,GOU Y,TARON J,et al.Thermo-hydro-me-chanical modeling of carbon dioxide injection for en-hanced gas-recovery(CO2-EGR):a benchmarking study for code comparison [J].Environmental Earth Sciences,2012,67(2):549-561. 

    • [14] LUO F,XU R N,JIANG P X.Numerical investigation of the influence of vertical permeability heterogeneity in stratified formation and of injection/production well per-foration placement on CO2 geological storage with en-hanced CH4 recovery [J].Applied Energy,2013,102:1314-1323. 

    • [15] REGAN M.A numerical investigation into the potential to enhance natural gas recovery in water-drive gas reser-voirs through the injection of CO2 [D].Adelaide:The University of Adelaide,2010. 

    • [16] CLASS H,EBIGBO A,HELMIG R,et al.A bench-mark study on problems related to CO2 storage in geolog-ic formations [J].Computational Geosciences,2009,13(4):409-434. 

    • [17] SEO J G,MAMORA D D.Experimental and simulation studies of sequestration of supercritical carbon dioxide in depleted gas reservoirs [J].Journal of Energy Re-sources Technology,2005,127:1-6. 

    • [18] ROBINSON B A,VISWANATHAN H S,VALOCCHI A J.Efficient numerical techniques for modeling multicom-ponent ground-water transport based upon simultaneous solution of strongly coupled subsets of chemical compo-nents [J].Advances in Water Resources,2000,23(4):307-324. 

    • [19] JIANG Y.Techniques for modeling complex reservoirs and advanced wells [D].Palo Alto:Stanford Universi-ty,2007. 

    • [20] PRUESS K.The TOUGH codes—a family of simulation tools for multiphase flow and transport processes in per-meable media [J].Vadose Zone Journal,2004,3:738-746. 

    • [21] PRUESS K,OLDENBURG C M,MORIDIS G J.TOUGH2 user＇s guide version 2.0[Z/OL].Rep LBNL-43134(1999-11-01)[2019-03-27].https://escholarship.org/uc/item/4df6700h. 

    • [22] PATEL M J,MAY E F,JOHNS M L.High-fidelity res-ervoir simulations of enhanced gas recovery with super-critical CO2 [J].Energy,2016,111:548-559. 

    • [23] BIAGI J,AGARWAL R,ZHANG Z.Simulation and optimization of enhanced gas recovery utilizing CO2 [J].Energy,2016,94:78-86. 

    • [24] SPAN R,WAGNER W.A new equation of state for car-bon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa [J].Journal of Physical and Chemical Reference Data,1996,25(6):1509-1596.

  • 基本信息

    中图分类号: TQ28
    文献标识码: A
    DOI: 10.3969/j.issn.1673-5005.2020.03.009
    文章编号: 1673-5005(2020)03-0081-09

    基金信息

    国家自然科学基金项目(U1762216)；
    教育部“长江学者冶创新团队项目(IRT_14R58)；

    引用信息

    刘树阳,孙宝江,宋永臣,等. CO2 提高天然气采收率的重力效应与气藏压力影响模拟[J]. 中国石油大学学报(自然科学版),2020,44(3):81-89.

    LIU Shuyang, SUN Baojiang, SONG Yongchen, et al. Simulation on gravity effect and reservoir pressure influence analysis in CO2 enhanced gas recovery[J]. Journal of China University of Petroleum(Edition of Natural Science),2020,44(3):81-89.

    稿件历史

    收稿日期: 2019-06-10

  • 参考文献

    • [1] CONINCK H C,CARBO M C,MIKUNDA T,et al.Car-bon capture and storage in industrial applications:tech-nology synthesis report[R/OL].(2010-12-09)[2019-03-27].https://publications.tno.nl/publication/34632076/JT069c/o10043.pdf. 

    • [2] 孙宝江,王金堂,孙文超,等.非常规天然气储层超临界 CO2 压裂技术基础研究进展[J].中国石油大学学报(自然科学版),2019,43(5):82-91.SUN Baojiang,WANG Jintang,SUN Wenchao,et al.Advances in fundamental research of supercritical CO2 fracturing technology for unconventional natural gas reser-voirs[J].Journal of China University of Petroleum(Edi-tion of Natural Science),2019,43(5):82-91. 

    • [3] 任韶然,张莉,张亮.CO2 地质埋存:国外示范工程及其对中国的启示[J].中国石油大学学报(自然科学版),2010,34(1):93-98.REN Shaoran,ZHANG Li,ZHANG Liang.Geological storage of CO2:overseas demonstration projects and its implications to China[J].Journal of China University of Petroleum(Edition of Natural Science),2010,34(1):93-98. 

    • [4] van der BURGT M J,CANTLE J,BOUTKAN V K.Car-bon dioxide disposal from coal-based IGCC＇s in depleted gas fields [J].Energy Conversion and Management,1992,33(5/6/7/8):603-610. 

    • [5] OLDENBURG C M,BENSON S M.CO2 injection for en-hanced gas production and carbon sequestration [ R].SPE 74367-MS,2002. 

    • [6] OLDENBURG C M,PRUESS K,BENSON S M.Process modeling of CO2 injection into natural gas reservoirs for carbon sequestration and enhanced gas recovery [J].En-ergy & Fuels,2001,15(2):293-298. 

    • [7] OLDENBURG C M,STEVENS S H,BENSON S M.E-conomic feasibility of carbon sequestration with enhanced gas recovery(CSEGR)[J].Energy,2004,29(9/10):1413-1422. 

    • [8] GANZER L,REITENBACH V,PUDLO D,et al.Exper-imental and numerical investigations on CO2 injection and enhanced gas recovery effects in Altmark gas field(Cen-tral Germany)[J].Acta Geotechnica,2013,9:39-47. 

    • [9] HUSSEN C,AMIN R,MADDEN G,et al.Reservoir simulation for enhanced gas recovery:an economic evalu-ation [J].Journal of Natural Gas Science and Engineer-ing,2012,5:42-50. 

    • [10] JIKICH S A,SMITH D H,SAMS W N,et al.En-hanced gas recovery(EGR)with carbon dioxide se-questration:a simulation study of effects of injection strategy and operational parameters [R].SPE 84813,2003. 

    • [11] REBSCHER D,OLDENBURG C M.Sequestration of carbon dioxide with enhanced gas recovery-case study altmark,North German Basin [ R ].LBNL-59033(2005-12-10)[2019-03-27].https://www.escholar-ship.org/uc/item/7sj298mc. 

    • [12] AL-HASAMI A,REN S,TOHIDI B.CO2 injection for enhanced gas recovery and geo-storage:reservoir simula-tion and economics[R].SPE 94129,2005. 

    • [13] HOU Z,GOU Y,TARON J,et al.Thermo-hydro-me-chanical modeling of carbon dioxide injection for en-hanced gas-recovery(CO2-EGR):a benchmarking study for code comparison [J].Environmental Earth Sciences,2012,67(2):549-561. 

    • [14] LUO F,XU R N,JIANG P X.Numerical investigation of the influence of vertical permeability heterogeneity in stratified formation and of injection/production well per-foration placement on CO2 geological storage with en-hanced CH4 recovery [J].Applied Energy,2013,102:1314-1323. 

    • [15] REGAN M.A numerical investigation into the potential to enhance natural gas recovery in water-drive gas reser-voirs through the injection of CO2 [D].Adelaide:The University of Adelaide,2010. 

    • [16] CLASS H,EBIGBO A,HELMIG R,et al.A bench-mark study on problems related to CO2 storage in geolog-ic formations [J].Computational Geosciences,2009,13(4):409-434. 

    • [17] SEO J G,MAMORA D D.Experimental and simulation studies of sequestration of supercritical carbon dioxide in depleted gas reservoirs [J].Journal of Energy Re-sources Technology,2005,127:1-6. 

    • [18] ROBINSON B A,VISWANATHAN H S,VALOCCHI A J.Efficient numerical techniques for modeling multicom-ponent ground-water transport based upon simultaneous solution of strongly coupled subsets of chemical compo-nents [J].Advances in Water Resources,2000,23(4):307-324. 

    • [19] JIANG Y.Techniques for modeling complex reservoirs and advanced wells [D].Palo Alto:Stanford Universi-ty,2007. 

    • [20] PRUESS K.The TOUGH codes—a family of simulation tools for multiphase flow and transport processes in per-meable media [J].Vadose Zone Journal,2004,3:738-746. 

    • [21] PRUESS K,OLDENBURG C M,MORIDIS G J.TOUGH2 user＇s guide version 2.0[Z/OL].Rep LBNL-43134(1999-11-01)[2019-03-27].https://escholarship.org/uc/item/4df6700h. 

    • [22] PATEL M J,MAY E F,JOHNS M L.High-fidelity res-ervoir simulations of enhanced gas recovery with super-critical CO2 [J].Energy,2016,111:548-559. 

    • [23] BIAGI J,AGARWAL R,ZHANG Z.Simulation and optimization of enhanced gas recovery utilizing CO2 [J].Energy,2016,94:78-86. 

    • [24] SPAN R,WAGNER W.A new equation of state for car-bon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa [J].Journal of Physical and Chemical Reference Data,1996,25(6):1509-1596.

    • 版权所有 中国石油大学学报（自然科学版）编辑部 Copyright©2008 All Rights Reserved
    主管单位：中华人民共和国教育部 主办单位：中国石油大学（华东）
    地址： 青岛市黄岛区长江西路66号中国石油大学期刊社 邮编：266580 电话：0532-86983553 E-mail： journal@upc.edu.cn
    本系统由：北京勤云科技发展有限公司设计