CN115034489A

CN115034489A - 一种考虑溶解作用的气藏co2埋存潜力预测方法

Info

Publication number: CN115034489A
Application number: CN202210699347.4A
Authority: CN
Inventors: 秦佳正; 刘梦云; 汪勇; 何佑伟; 汤勇
Original assignee: Southwest Petroleum University
Current assignee: Southwest Petroleum University
Priority date: 2022-06-20
Filing date: 2022-06-20
Publication date: 2022-09-09
Anticipated expiration: 2042-06-20
Also published as: CN115034489B

Abstract

本发明公开了一种考虑溶解作用的气藏CO₂埋存潜力预测方法，该方法基于新的物质的量平衡方法，考虑CO₂在气藏本身的地层水以及外部侵入的水中的溶解作用，建立气藏CO₂埋存量预测模型。使用软件编程运算，结合具体气藏实例，预测气藏CO₂埋存量潜力。对比计算结果显示，新的物质的量平衡方法考虑溶解的气藏CO₂埋存量比传统物质平衡法不考虑溶解的气藏CO₂埋存量多1.0969倍。本发明解决了传统物质平衡法中未考虑溶解的问题，使气藏CO₂埋存潜力预测结果更精确。

Description

一种考虑溶解作用的气藏CO2埋存潜力预测方法

技术领域

本发明涉及CO₂埋存量技术领域，尤其是一种考虑溶解作用的气藏CO₂埋存潜力预测方法。

背景技术

随着世界范围内化石燃料的使用，CO₂排放量在不断增加，过量排放CO₂将会带来严重的环境影响。CO₂气体具有吸热和隔热的功能。过量排放CO₂的结果是在地球表面形成一种无形的玻璃罩，使太阳辐射到地球上的热量无法向外层空间发散，其结果是地球表面变热。这种变化也被称为温室效应，CO₂也被称为温室气体。温室效应会导致海平面上升、气候反常，海洋风暴增多、土地干旱，沙漠化面积增大和病虫害增加等。为了解决和应对温室效应，国内外的研究学者提出了很多种解决方法，其中CO₂捕集、利用和埋存技术(CCUS)被认为是国家和全球减排方案的核心组成部分。有学者认为CCUS将会成为对CO₂减排贡献最大的单个技术，并预测其贡献将从2020年占总减排量的3％提高到2030年的10％，再到2050年的20％左右。这一观点得到国内外的普遍认同。CO₂埋存的经济效益同样不可忽略。在常规煤化工、钢铁、水泥等高CO₂排放行业，处理排放出来的CO₂需要花费的经济远超于埋存的成本。同样，对于温室效应带来的严重环境影响上，解决这些问题所消耗的社会成本和经济成本远大于CO₂埋存技术的成本。因此CCUS技术的发展将给世界带来巨大的环境效益和经济效益。

CO₂埋存地点包括：枯竭油藏、枯竭气藏、盐水层、煤层等多种储层。针对不同的埋存地点，埋存机理也不尽相同，埋存机理主要有构造埋存、溶解埋存、残余气埋存、矿化埋存、水动力埋存等。天然气比液态的石油具有更大的压缩性，且气藏采收率一般比油藏采收率高，所以与枯竭油藏相比，枯竭气藏具有更多的CO₂埋存空间，中国天然气储量排名世界第六，探明储量达12.468万亿立方米，气藏资源丰富，因此气藏为CO₂埋存的有利场所。

气藏CO₂埋存量计算方法包括容积法、压缩系数法、类比法、概率预测法、数值模拟法、传统物质平衡法等。容积法主要考虑储层属性如岩性、孔隙度、渗透率、面积、厚度等，结合CO₂埋存效率即注入CO₂所占的体积占据孔隙体积的比例；类比法通过对比不同气藏的属性特征，简单估算相似气藏的埋存量的大致范围，得到理论埋存量，因此计算具有很大的误差；概率预测法通过建立概率模型来求解埋存量，其中考虑了输入参数的不确定性，但其所需储层数据量过大，部分参数在矿场应用中很难获取；数值模拟法通过建立数值模型模拟埋存过程，建模过程复杂、模拟时间长，且参数获取难度较大。因此，上述方法在矿场应用中适用性较差。

传统物质平衡法与数值模拟法和概率预测法等动态方法相比，计算过程更简单，效率更高；且计算所需参数均可根据矿场实际参数来取值，解决了数值模拟法和概率预测法中参数不易获取的问题。此外，传统物质平衡法相较于容积法、压缩系数法、类比法等静态法而言，考虑影响埋存量的因素更全面，计算结果更贴近实际，且具有更高的精度。但传统物质平衡方法存在一定问题：未考虑到CH₄、CO₂的溶解作用对埋存量的影响。当温度为37℃、压力为0.1MPa时，CO₂在地层水中的溶解度能够达到50kg/m³，而地层温度和压力下更有利于CO₂的溶解。因此，在进行气藏CO₂埋存时，溶解捕集CO₂的埋存机制不能忽略。

在矿场实际应用中，溶解作用对CO₂埋存潜力的影响不能忽略。因此，亟需研究基于物质的量平衡原理考虑CO₂溶解作用预测气藏CO₂埋存潜力的有效方法。

发明内容

本发明的目的是提供一种基于物质的量平衡原理考虑CO₂溶解作用预测气藏CO₂埋存潜力的方法。该方法采用以下三个公式之一计算CO₂埋存量G_CO2或者CO₂埋存量与累计采出井流物比值R_C-N：

其中，

式中，

是剩余混合气时的地层压力，MPa；c_e是有效压缩系数；p_i是原始地层压力，MPa；G_i是地质储量，m³；Z_fgi是原始自由气偏差因子；W_e是水侵量，m³；W_p是累计产水量，m³；B_w是混合状态下地层水体积系数，m³/m³；Z_sc是标况下偏差因子；p_sc是标况下压力，MPa；c_w是地层水压缩系数；M是水体倍数；S_wc是束缚水饱和度；

是混合状态下溶解气水比，m³/m³；S_wi是原始含水饱和度；R_swi是原始状态下溶解气水比，m³/m³；B_wi是原始状态下地层水体积系数，m³/m³；G_p是累计井流物产量，m³；

是剩余混合气偏差因子；

是剩余混合气温度，K；T_sc是标况下温度，K；T_i是原始地层温度，K；

是CO₂埋存量，m³；R_C-N是CO₂埋存量与累计采出井流物比值；R_N-I是累计采出井流物与地质储量的比值。

根据上述的CO₂埋存量计算公式，采用软件编程计算。

针对具体气藏，对CO₂埋存量计算公式中的参数进行求解和赋值，计算得到CO₂埋量。

上述计算公式中，水侵量W_e的计算公式如下：

W_e＝G_pB_g+W_pB_w-N_g(B_g-B_gi)

式中，N_g是气藏的原始地质储量，m³；G_p是累积产气量，m³：B_gi是气藏天然气的原始体积系数，m³/m³；B_g是气藏天然气的目前的体积系数，m³/m³；W_p是累计产水量，m³；B_w是混合状态下地层水体积系数，m³/m³。

进一步的，气藏天然气的原始体积系数B_gi和目前的体积系数B_g的计算公式分别为：

式中，p_sc是标况下压力，MPa；Z_i是原始天然气偏差系数；T_i是气藏初始温度，K；P_i是原始地层压力，MPa；T_sc是标况下温度，K；Z是目前天然气偏差系数；T是目前的气藏温度，K；P是目前地层压力，MPa。

有效压缩系数c_e由以下公式计算得到：

式中，c_f是岩石压缩系数；c_w是地层水压缩系数；S_wc是束缚水饱和度；S_wi是原始含水饱和度。

与现有技术相比，本发明的有益之处在于：

目前缺乏基于物质的量平衡原理考虑CO₂溶解预测气藏CO₂埋存潜力的有效方法，导致评估结果出现偏差，在矿场实际应用中，应重点考虑溶解对CO₂埋存量计算产生的影响。本方法开展的研究使人们更好地理解了CO₂埋存量计算过程，这对优化现场CO₂埋存量计算具有重要意义。本发明提出一种考虑CO₂溶解作用的新的物质的量平衡的方法。与传统物质平衡法相比，重点考虑了CO₂在地层水中的溶解效应对气藏CO₂埋存量的影响。新的物质的量平衡方法中对溶解的考虑是动态全面的，包括流体组成、温度、压力等参数对气体溶解于地层水中的影响，同样考虑偏差系数、体积系数随着压力和温度的变化。最终的实例计算结果表明，考虑溶解效应的新的物质的量平衡方法比传统不考虑溶解的物质平衡法多1.0969倍的CO₂埋存量，解决了传统方法预测CO₂埋存量偏低的问题。在传统物质平衡法的基础上，考虑CO₂溶解埋存机制后，埋存量更大，CO₂埋存量预测结果更精确。

本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。

附图说明

图1、物质平衡方法解释图。

图2、实施例中采出程度、相对压力和水侵系数的关系曲线图。

图3、CO₂埋存量随地层压力的变化图。

图4、CO₂埋存量与累计采出井流物的比值与R_C-N关系曲线图。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

如图1所示，给出了物质平衡方法的具体解释图，其中从左到右的矩形分别代表气藏从生产期到枯竭，再到埋存的时间顺序，各个字母代表含义：(a)：剩余混合气；(b)：剩余水溶混合气；(c)：原始气；(d)：原始水溶烃；(e)：采出井流物；(f)：注入CO₂量；(g)：生产至枯竭时地层残余气。

物质平衡方程及其推导过程如下：

摩尔守恒基本关系为：

剩余混合气摩尔量+剩余水溶混合气摩尔量＝原始气相摩尔量+原始水溶烃气摩尔量-采出井流物摩尔量+注入CO₂摩尔量，即表达式如下：

式中：

剩余混合气摩尔量，kmol；

剩余水溶混合气摩尔量，kmol；

n_fgi：原始气相摩尔量，kmol；

n_sgi：原始水溶气摩尔量，kmol；

n_p：采出的井流物摩尔量，kmol；

注入CO₂摩尔量，kmol。

剩余混合气摩尔量

计算公式如下：

式中：

剩余混合气摩尔量，kmol；

剩余混合气时的地层压力，MPa；

剩余混合气偏差因子；

当前混合气孔隙体积，m³；

剩余混合气时的地层温度，K；

R：通用气体常数，MPa·m³/(kmol·K)；

V_HCi：原始烃类孔隙体积，m³；

W_p：累计产水量，m³；

p_i：原始地层压力，MPa；

c_e：有效压缩系数；

W_e：水侵量；

B_w：混合状态下地层水体积系数；

c_f：岩石压缩系数；

c_w：地层水压缩系数；

S_wc：束缚水饱和度；

S_wi：原始含水饱和度。

剩余水溶气摩尔量

计算公式如下：

式中，

剩余水溶气摩尔量，kmol；

p_sc：标况下压力，MPa；

M：水体倍数，定义为水体体积与气藏原始孔隙体积之比；

Z_sc：标况下偏差因子；

R：通用气体常数，MPa·m³/(kmol·K)；

T_sc：标况下温度，K；

混合状态下溶解气水比，m³/m³；

B_w：混合状态下地层水体积系数，m³/m³；

p_i：原始地层压力，MPa；

W_e：水侵量。

原始气相摩尔量n_fgi根据公式(6)计算：

式中，

n_fgi：原始气相摩尔量，kmol；

p_i：原始地层压力，MPa；

V_HCi：原始烃类孔隙体积，m³；

Z_fgi：原始自由气偏差因子；

R：通用气体常数，MPa·m³/(kmol·K)；

T_i：原始地层温度，K。

原始水溶气摩尔量n_sgi按照公式(7)计算：

式中，

n_sgi：原始水溶气摩尔量，kmol；

V_AQ：水体体积，m³；

V_wc：束缚水体积，m³；

R_swi：原始状态下溶解气水比，m³/m³；

B_wi：原始状态下地层水体积系数，m³/m³；

M：水体倍数，定义为水体体积与气藏原始孔隙体积之比。

采出的井流物摩尔量n_p按照公式(9)得到：

n_p：采出的井流物摩尔量，kmol；

p_sc：标况下压力，MPa；

G_p：累计井流物产量，m³；

Z_sc：标况下偏差因子；

R：通用气体常数，MPa·m³/(kmol·K)；

T_sc：标况下温度，K。

注入CO₂摩尔量

根据公式(10)得到：

注入CO₂体积，m³；

注入CO₂摩尔量；

p_sc：标况下压力，MPa；

Z_sc：标况下偏差因子；

R：通用气体常数，MPa·m³/(kmol·K)；

T_sc：标况下温度，K。

将(2)、(3)、(5)、(6)、(7)、(9)、(10)代入(1)式中，得到公式(11)：

由(12)、(13)可以得到：

V_HCi＝G_i·B_gi (14)

将(14)代入(11)式中：

将(16)代入(15)式中：

消去相同项得：

消去

项和

得：

可以得到

为：

或合并G_i：

其中：

变为R_C-N和N_P-Z形式：

同除G_p得：

其中：

式中，

是剩余混合气偏差因子；

是CO₂埋存量，m³；R_C-N是CO₂埋存量与累计采出井流物的比值；R_N-I是累计采出井流物与地质储量的比值。

通过推导可以得出式(20)、(21)、(22)都是不同形式的考虑溶解的气藏埋存量计算公式。公式中涉及的参数，除了水侵量W_e外，其他参数对于某个具体气藏来说，都是已知的或者可以通过气藏资料简单计算得出。即W_e属于未知参数，其他参数都属于已知参数。然后通过软件编程计算。

本发明的气藏埋存量计算方法应用到某个具体气藏实例中，方法如下:

目标求出：

：累计注入CO₂量，m³。

已知的参数取值见表1。

表1、某气藏已知参数的取值

未知需要计算的参数是水侵量W_e(m³)，计算得到的值为128943.439。具体计算方法如下：

水侵量W_e计算公式如下：

W_e＝G_pB_g+W_pB_w-N_g(B_g-B_gi)

带入本次气藏埋存量案例计算中的数据进行气藏水侵量计算。按照上述方法进行气藏水侵量计算，整个计算过程从气藏开始开采到产出0.5159×10⁸m³时对应地层压力为18MPa时侵入地层的水侵量大小。

已知：

原始地层压力：P_i＝37MPa；

原始天然气偏差系数：Z_i＝0.9803；

目前地层压力：P＝14.67MPa；

目前天然气偏差系数：Z＝0.9306；

累积产气量：G_p＝0.5159×10⁸m³；

原始地质储量：N_g＝0.8×10⁸m³；

累计产水量：W_p＝8.43×10⁴m³；

混合状态下地层水体积系数:B_w＝1.0383。

计算气藏的地层相对压力

气藏的目前采出程度R_D：

将计算得到的(0.6449,0.4177)点代入图2所示的曲线图，该点在图中的位置明显位于45°对角线之上的三角形区域，所以可以判断出该气藏存在水侵，通过下式计算出水侵体积系数为：

同时从曲线图上也可以判断出水侵体积系数数值在0.14左右附近，代表计算结果可用，再通过下式可以计算出气藏被水侵占据的有效孔隙体积大小：

最后通过式下式计算出气藏水侵量大小：

W_e＝G_pB_g+W_pB_w-N_g(B_g-B_gi)＝0.5159×10⁸×0.00828659+8.43×10⁴×1.0378-0.8×10⁸(0.00828659-0.00346099)＝128943.439

通过软件编程计算，结果如图3和图4所示。图3是埋存量随着压力的变化图。从图3可以明显得出，当储层压力为具体指时，对应的埋存量。图4是CO₂埋存量与累计采出井流物的比值与R_C-N关系曲线图，通过压力和其他一些参数可以计算出N_P-Z大小，带入曲线得出R_C-N大小，由于产出量已知，所以也可以得出注的量。

通过对比分析两种情况下的气藏埋存量计算，即不考虑溶解和考虑溶解的气藏埋存量计算，其中每种又具体分为有水侵量和无水侵量的情况，同时也考虑使用Bachu的方法计算理论埋存量，对比结果见表2。

表2、气藏CO₂埋存量计算结果对比

上表中，不考虑溶解的的气藏埋存量计算公式参考现有文献：Chi-Chung Tseng,Bieng-Zih Hsieh,Shin-Tai Hu,Zsay-Shing Lin,Analytical approach for estimatingCO₂ storage capacity of produced gas reservoirs with or without a waterdrive,International Journal of Greenhouse Gas Control,Volume 9,2012,Pages254-261,ISSN 1750-5836。公式中各参数的含义在该文献中有详细解释。

通过软件编程计算出上述表格中的具体数据。可以看出，两种方法，即考虑不考虑溶解和考虑溶解的计算基础数据一致，但是与传统物质平衡方法相比，新的物质的量平衡法的计算结果更有利于气藏CO₂埋存。两种方法都充分考虑水侵量的影响，不考虑溶解和考虑溶解的气藏CO₂埋存量结算结果分别是97050吨和203500吨。其中考虑溶解的新的物质的量平衡的计算方法的假设基础是充分考虑CO₂在水中的溶解，计算得出考虑溶解的气藏埋存量比不考虑溶解的气藏埋存量多106450吨，多出1.0969倍的CO₂埋存量。

以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。