热-流耦合与热-流-固耦合作用下的水气交替及间歇注入对咸水层CO2溶解封存的影响

张杰城, 杜鑫芳, 赫文豪, 张来斌

张杰城, 杜鑫芳, 赫文豪, 张来斌. 热-流耦合与热-流-固耦合作用下的水气交替及间歇注入对咸水层CO2溶解封存的影响[J]. 华南师范大学学报(自然科学版), 2024, 56(1): 9-17. DOI: 10.6054/j.jscnun.2024002
引用本文: 张杰城, 杜鑫芳, 赫文豪, 张来斌. 热-流耦合与热-流-固耦合作用下的水气交替及间歇注入对咸水层CO2溶解封存的影响[J]. 华南师范大学学报(自然科学版), 2024, 56(1): 9-17. DOI: 10.6054/j.jscnun.2024002
ZHANG Jiecheng, DU Xinfang, HE Wenhao, ZHANG Laibin. Impact of Water Alternating Gas and Intermittent Injection on CO2 Dissolution During CO2 Sequestration in Saline Aquifers Considering Thermal-Hydraulic Coupling and Thermal-Hydraulic-Mechanical Coupling Effects[J]. Journal of South China Normal University (Natural Science Edition), 2024, 56(1): 9-17. DOI: 10.6054/j.jscnun.2024002
Citation: ZHANG Jiecheng, DU Xinfang, HE Wenhao, ZHANG Laibin. Impact of Water Alternating Gas and Intermittent Injection on CO2 Dissolution During CO2 Sequestration in Saline Aquifers Considering Thermal-Hydraulic Coupling and Thermal-Hydraulic-Mechanical Coupling Effects[J]. Journal of South China Normal University (Natural Science Edition), 2024, 56(1): 9-17. DOI: 10.6054/j.jscnun.2024002

热-流耦合与热-流-固耦合作用下的水气交替及间歇注入对咸水层CO2溶解封存的影响

基金项目: 

国家重点研发计划项目 2022YFC2806401

国家自然科学基金项目 52301354

中国石油大学(北京)科研基金项目 2462023SZBH004

中国博士后科学基金项目 2023M743875

详细信息
    通讯作者:

    赫文豪,Email:hwh@cup.edu.cn

  • 中图分类号: P641

Impact of Water Alternating Gas and Intermittent Injection on CO2 Dissolution During CO2 Sequestration in Saline Aquifers Considering Thermal-Hydraulic Coupling and Thermal-Hydraulic-Mechanical Coupling Effects

  • 摘要:

    基于TOUGH+软件架构,结合更为准确的物理性质计算模型,建立了适用于咸水层CO2封存的热-流耦合模拟方法,并使用固定应力分割迭代耦合模型将其与RGMS软件耦合,改进了热-流-固双向耦合模拟方法,建立了准确性更高的热-流-固迭代耦合模拟方法。基于鄂尔多斯盆地地质特征构建的地质模型,使用热-流与热-流-固迭代耦合方法模拟咸水层CO2溶解封存过程,研究了水气交替及间歇注入方案对CO2溶解量、孔隙压力和地层形变的影响。结果表明:热-流-固迭代耦合模拟能帮助设计更加合理的注入方案;仅用水气交替注入方式可提高CO2溶解量;间歇注入有助于孔隙压力与地层形变恢复。研究结果可为咸水层CO2溶解封存提供理论指导。

    Abstract:

    Based on TOUGH+ and a more accurate physical property calculation model, a thermal-hydraulic coupling simulation method suitable for CO2 sequestration in saline aquifers is established. The method was further coupled with RGMS using the fixed-stress split iteration model, improved from the fixed-stress split two-way model, to obtain a more accurate thermal-hydraulic-mechanical iteratively coupling simulation method. Based on the geological characteristics of the Ordos Basin, a geological model is constructed, and the thermal-hydraulic and thermal-hydraulic-mechanical coupling simulation methods mentioned above are used to simulate the process of CO2 dissolution during CO2 sequestration in saline aquifers. The impacts of water alternating gas injection and intermittent injection schemes on CO2 dissolution, pore pressure, and formation deformation are investigated. The following results are found: the thermal-hydraulic-mechanical iteratively coupling model can better assist in designing more reasonable injection schemes; only water alternating gas injection can enhance CO2 dissolution; intermittent injection contributes to the recovery of pore pressure and formation deformation. The results can provide theoretical guidance for the CO2 sequestration in saline aquifers by CO2 dissolution.

  • 聚羧酸系减水剂作为新一代高性能减水剂,具有较好的分散性、无污染生产以及分子结构可设计等优点,在高性能混凝土中得到普遍应用[1-3],但其易引起坍落度损失的问题一直是困扰工程施工的难题,特别是长时间的运输、超高层泵送等会导致混凝土坍落度损失过大,影响混凝土的工作性[4-6], 因此,保坍型减水剂受到广泛关注。提高聚羧酸盐减水剂保坍能力的研究主要集中在合成含有酰胺基、马来酸酐和酯等可水解官能团的聚合物[7-8]。这类缓释型聚羧酸减水剂(PCE)将在水泥浆体的碱性环境中逐渐水解,产生对水泥有分散作用的羧酸阴离子基团,从而达到保坍效果。为了改善水泥浆的流动保持性能,LIU等[9]以丙烯酸羟乙酯为功能单体合成了几种相对分子质量不同的聚羧酸减水剂,相对分子质量较低的减水剂具有良好的分散保持性能。李安等[10]研究了羧酸羟乙基酯、羧酸羟丙基酯对合成缓释减水剂性能的影响,丙烯酸羟丙基酯作为缓释功能单体且与大单体甲基烯丙基聚氧乙烯醚(HPEG)的物质的量之比为5 ∶ 1时,混凝土的坍落度损失较小。LIU等[11]以烯丙基羧酸二甲酯合成了一种缓释效果良好的减水剂,水合物放热曲线表明,缓释减水剂只是延缓了水化,并不影响混凝土的水化热及和易性。但系统研究酯类单体的类型对水泥浆体系分散性能的影响以及寻求更长缓释性能聚羧酸减水剂的文献鲜有报道。

    该研究对羧基官能团进行保护、修饰,采用丙烯酸羟乙酯(HEA)、丙烯酸羟丙酯(HPA)、衣康酸二甲酯(DEI)、富马酸二甲酯(DMF)为功能单体合成了4种缓释型聚羧酸减水剂。通过对比不同功能单体的PCE对水泥浆体分散性能、吸附行为、表面张力的影响规律,研究不同减水剂对混凝土的坍落度保持率的影响,对实际生产中制备高保坍、高性能聚羧酸减水剂具有很好的参考价值。

    主要试剂:甲基烯丙基聚氧乙烯醚(HPEG-2400)为工业级,购于辽宁奥克化学股份有限公司;丙烯酸羟乙酯(HEA)、丙烯酸羟丙酯(HPA)、衣康酸二甲酯(DEI)、富马酸二甲酯(DMF,阿达玛斯)、丙烯酸(AA)、过氧化氢(H2O2,质量分数30%)、抗坏血酸(VC)、巯基乙酸(TGA)、氢氧化钠(NaOH)均为分析纯,市售;去离子水由实验室自制。水泥(P ·O42.5)由尧柏特种水泥集团有限公司提供,粉煤灰为大唐电厂Ⅱ级粉煤灰,砂为细度模数为2.6的天然砂,石料为粒径5~31 mm的连续级配碎石。

    主要仪器:傅里叶变换红外光谱仪(VERTEX-80, 德国Bruker)、核磁共振波谱仪(ADVANCE Ⅲ 600 MHz, 德国Bruker)、凝胶渗透色谱仪(GPC,Waters GPC 1515,美国Waters)、NJ-160A水泥净浆搅拌机(上海东星建材试验设备)、电位粒度分析仪(Zetasizer NANO-ZS90, 英国Malvern)、总有机碳分析仪(vario TOC, 德国Elementar)、表面张力测定仪(德国Dataphysics)、高速离心机(Kratos TG650-WS,上海安亭科学仪器)、纯水制造机(VE-Ds,宏森)。

    采用水溶液自由基聚合制备4种具有不同缓释单体的PCE。将20.5 g HPEG-2400和去离子水加到配有搅拌器和恒压滴液漏斗的三口圆底烧瓶中,搅拌使其溶解,随后加入0.14 g H2O2,待温度升至45 ℃,同时缓慢加入溶液A和B(溶液A为AA及其功能单体的混合溶液,溶液B为VC和链转移剂TGA的水溶液),并控制溶液A滴加3 h,溶液B持续3.5 h。滴完后,继续保温1~2 h。反应完毕后加入NaOH水溶液(质量分数30%)将PH调节为6~8。所得浅黄色或无色粘稠透明溶液为缓释型多元聚羧酸减水剂PCE(40%,固体质量分数)。样品的命名方式:PCE和缓释功能单体的缩写组成。在前期大量实验中,确定了最佳酸醚比n(酸) ∶n(醚)=3.25 ∶ 1,减水剂PCE各单体投料比(物质的量之比)及聚合物多分散性系数(PDI)见表 1, 4种减水剂的合成路线如图 1所示。

    表  1  PCE中各单体的用量及聚合物多分散性系数(PDI)
    Table  1.  The dosage of each monomer in PCE and the polymer dispersity index(PDI)
    样品 投料比/mol Mn/(g·mol-1) Mw/(g·mol-1) PDI
    PCE-HEA 3.25∶1∶3.35 52 651 81 609 1.55
    PCE-HPA 3.25∶1∶3.43 46 752 76 673 1.64
    PCE-DEI 3.25∶1∶2.67 58 963 105 543 1.79
    PCE-DMF 3.25∶1∶2.34 48 721 74 055 1.52
    注:投料比=c(AA) ∶ c(HPEG) ∶ c(HEA、HPA、DEI或DMF)。
    下载: 导出CSV 
    | 显示表格
    图  1  含不同酯类单体PCE的合成路线
    Figure  1.  The synthetic route of PCE containing different esters

    将4种缓释型减水剂PCE分别置于透析袋里透析一段时间, 采用核磁共振波谱仪(以D2O为溶剂), 测定减水剂的核磁共振氢谱(1H NMR)。将干燥后的样品用KBr压片, 采用傅里叶变换红外光谱仪在400~4 000 cm-1的波数范围测定减水剂的红外光谱(FTIR)。

    依照GBT 8077—2012《混凝土外加剂匀质性实验方法》[12],测定水泥的净浆流动度。水灰质量比固定为0.29,PCE掺入量(按固体质量折算)为水泥质量的0.2%。每隔30 min测试其流动度,测试总时间为180 min。

    参照GB 50080—2016《普通混凝土拌合物性能试验方法标准》[13]检测混凝土的坍落度、扩展度等。混凝土各组分质量配比为∶ m(水泥) ∶ m(粉煤灰) ∶m(砂) ∶ m(石) ∶ m(水) ∶ m(外加剂)= 280 ∶ 80 ∶ 800 ∶1 080 ∶ 160 ∶ 7.2。

    分别取0.2 g水泥及20 mL质量浓度为500 mg/L的减水剂溶液置于锥形瓶中,在25 ℃的恒温水浴中振荡不同时间,随后取上层悬浮液用高速离心机分离,并用微孔滤膜过滤。将上层清液稀释至测量范围后,采用总有机碳(TOC)分析仪测定上层清液中PCE的有机碳质量浓度。

    取多份0.2 g水泥于锥形瓶中,分别加入50 mL一系列不同质量浓度的减水剂溶液。混合后,利用恒温摇床在25 ℃下振荡2 h,静置、离心、取少量上层清液,采用电位粒度分析仪测定水泥颗粒表面的Zeta电位。另外,将不加减水剂的水泥处理后作为对照组,测定吸附PCE前后水泥颗粒表面的Zeta电位。

    使用表面张力仪测试缓释减水剂PCE溶液水解前后的表面张力。将加热冷却后的铂金浸入到减水剂溶液中,记录测量值。同一组溶液连续测量3次,取其平均值。在减水剂溶液中加入饱和的NaOH溶液后,在25 ℃下搅拌1~2 h(水解前后的溶液均经标准盐酸溶液调节pH至中性后进行测定)。

    向pH 12.0的NaOH溶液中,加入一定比例的羧酸酯,并用pH 13.0的NaOH溶液中和滴定,记录每次消耗pH 13.0 NaOH溶液的体积Vt,由Vt计算出锥形瓶中未反应的羧酸酯在t时刻的质量浓度[14]

    含不同功能单体减水剂PCE的红外光谱如图 2所示。由于4种聚合物分子红外光谱(FTIR)的出峰位置基本相同,此处以PCE-HEA为代表进行分析。波数3 452 cm-1归属于合成聚合物分子中长侧链末端羟基的伸缩振动峰,2 877 cm-1归属于聚合物分子中饱和C—H键的伸缩振动峰,1 730 cm-1归属于聚合物分子侧链酯基上C=O的伸缩振动吸收峰[15],说明酯类单体成功聚合。1 620~1 680 cm-1范围内的C=C双键特征吸收峰消失,说明所有单体已聚合完全。1 560 cm-1和1 468 cm-1分别归属于羧酸盐的反对称伸缩振动和对称伸缩振动峰,1 344 cm-1归属于甲基的对称变形振动峰。1 106 cm-1归属于聚合物长侧链重复单元的伸缩振动峰,955 cm-1归属于-CH2CH2O-亚甲基的摇摆振动,844 cm-1归属于-CH2CH2O-的变形振动吸收峰。

    图  2  PCE样品的红外光谱
    Figure  2.  The FTIR spectra of PCE samples

    图 3为4种聚合物的核磁共振氢谱,化学位移δ=4.7属于溶剂峰,减水剂在δ=5.0~6.0范围那些属于不饱和碳上氢质子的信号已消失,表明所有单体基本参与了聚合反应。δ=3.8~4.1范围归属于共聚物中PDEI链段上的亚甲基;δ=3.2~3.8的多重强峰对应于减水剂分子长侧链上重复单元的亚甲基氢;δ=0.7~1.8归属于减水剂分子的甲基氢或亚甲基氢,δ=2.2归属于分子主链上的次甲基氢。由于丙烯酸的羧基易形成氢键,很难被检测出,所以核磁共振氢谱中未发现羧基氢的化学位移。综合红外光谱和核磁共振氢谱,证实目标产物分子结构与设计一致。

    图  3  PCE样品的1H NMR谱
    Figure  3.  The 1H NMR spectra of PCE samples

    掺入不同减水剂PCE后,水泥浆的流动性随时间的变化如图 4所示,掺有4种PCE的水泥浆流动度均出现了随时间延长而逐渐增大的趋势(示意图如图 5所示),其中PCE-HEA含有羟乙基酯基,PCE-HPA含有羟丙基酯基,而PCE-DEI和PCE-DMF含有2个酯基。可以看出,掺有聚羧酸减水剂(PCE-HEA和PCE-HPA)的水泥浆流动度随着时间的延长逐渐增大,在60 min时达到最大值,随后逐渐减小,而且PCE-HEA比PCE-HPA的增幅明显。掺有PCE-DEI和PCE-DMF水泥浆的流动度在120 min时达到最大值,缓释效果最好。这是因为,酯在碱性环境下的水解速率与酯基所连接基团的电子效应和空间位阻有关,其中PCE-HPA比PCE-HEA多1个-CH2-,即PCE-HPA的供电子基团和空间位阻比PCE-HEA的大,因此PCE-HPA的水解速率比PCE-HEA的小。此外,受酯基数量的影响,羧酸双酯的水解速率比羧酸单酯的低,所以掺有PCE-DEI和PCE-DMF减水剂的水泥浆体流动度保持效果最好。

    图  4  含不同PCE减水剂水泥浆的流动度
    Figure  4.  The fluidity of cement paste containing different PCE plasticizers
    图  5  PCE共聚酯基的水解及其在水泥颗粒上的吸附
    Figure  5.  The hydrolysis of ester group of PCE and its adsorption on cement particle

    PCE在水泥颗粒表面的吸附动力学曲线如图 6所示,不同功能单体PCE的整体变化趋势基本相同。随着时间的延长,吸附量增大,在30 min左右趋于平缓,说明PCE在水泥颗粒表面的吸附达到平衡,与Langmuir吸附等温方程相符合[16]。同时还可以看出,不同功能单体PCE对水泥的吸附量不同,在PCE掺入量相同时,吸附量从大到小顺序为PCE-DEI、PCE-DMF、PCE-HEA、PCE-HPA。这是因为随着时间的延长,酯基逐渐水解为COO-并吸附到水化的水泥颗粒上[17],PCE-DEI和PCE-DMF结构中含有双酯基,所以吸附能力强于PCE-HEA和PCE-HPA;PCE-HEA由于相同时间内释放的COO-多于PCE-HPA,故吸附能力较强,吸附量较大。这与流动性的实验结果一致。

    图  6  不同PCE减水剂在水泥浆体中的吸附量
    Figure  6.  The adsorption capacity of different PCE plasticizers in cement paste

    Zeta电位可以用来衡量水泥浆体的分散及稳定性能[18]图 7为减水剂质量浓度对水泥颗粒表面ζ电位的影响,无PCE时水泥的ζ=-0.74 mV。当4种减水剂掺入后,ζ的绝对值迅速变大, 说明水泥颗粒表面均与减水剂产生吸附作用。随着减水剂质量浓度的增加,水泥颗粒表面ζ的绝对值略微减小,这主要是由于COO-取代了水泥颗粒表面部分类似于SO42-的高价阴离子[19]。当减水剂质量浓度达到一定值时,ζ绝对值的变化量较小,此时吸附已接近平衡。饱和时水泥颗粒ζ的绝对值从大到小顺序:PCE-DEI(6.12 mV)、PCE-DMF(5.88 mV)、PCE-HEA(5.59 mV)、PCE-HPA(5.21 mV)。ζ绝对值越大,说明水泥胶粒表面的阴离子电荷强度越高[20],而缓释型减水剂的阴离子源于其分子结构中羧酸酯的水解,PCE-DEI和PCE-DMF含有2个酯基,水解后具有较强的吸附性,两者的ζ绝对值比PCE-HEA和PCE-HPA的大。

    图  7  含不同浓度PCE减水剂水泥颗粒的Zeta电位
    Figure  7.  The Zeta potential values of cement particles containing different concentrations of PCE plasticizers

    缓释型PCE(0.35 g/L)水解前后溶液的表面张力测试结果如图 8所示。实验测得纯水的表面张力为73.45 mN/m,当向水中加入4种不同的减水剂后,表面张力均减小,这是由表面活性剂的特性决定的。4种功能单体的亲水性从大到小顺序:丙烯酸羟乙酯(HEA)、丙烯酸羟丙酯(HPA)、衣康酸二甲酯(DEI)、富马酸二甲酯(DMF)。水解前减水剂PCE溶液的表面张力从大到小顺序:PCE-DMF(60.46 mN/m)、PCE-DEI(59.34 mN/m)、PCE-HPA(55.31 mN/m)、PCE-HEA(51.11 mN/m);水解后表面张力从大到小顺序:PCE-HPA(51.11 mN/m)、PCE-HEA(48.76 mN/m)、PCE-DMF(50.42 mN/m)、PCE-DEI(47.44 mN/m)。表面张力与表面活性剂在气/液界面上的平衡吸附密度成反比,随着亲水基团数量的增加,分子之间的斥力减弱,导致减水剂在气液界面的吸附密度增大,表面活性增强,从而导致表面张力减小[21]。水解前减水剂PCE-DEI和PCE-DMF分子结构中疏水基团数量大,使分子间缔合作用增强,亲水结构分布在聚集的疏水结构内部,表现为水溶性降低、表面张力增大。由于缓释型聚羧酸减水剂分子结构中含有疏水的酯基,所以表面张力较大,但酯基水解为COO-后,亲水性增强,表面张力减小。总之,掺有减水剂PCE溶液表面张力的测定结果与净浆流动度的结果一致。

    图  8  减水剂PCE的表面张力
    Figure  8.  The surface tension of PCE plasticizers

    通过混凝土试验对4种缓释型减水剂PCE性能进行对比,具体方法为通过4种缓释型减水剂与减水型减水剂复配制备混凝土添加剂,通过混凝土坍落度与扩展度测试对所制添加剂性能进行验证。4种PCE对混凝土坍落度经时损失的实验数据见表 3

    表  3  混凝土坍落度损失随时间的变化
    Table  3.  The changes of concrete slump loss over time
    样品 t=0 min t=60 min t=120 min
    土坍落 扩展度 土坍落 扩展度 土坍落 扩展度
    PCE-HEA 220 560 220 530 190 475
    PCE-HPA 230 585 220 540 185 470
    PCE-DEI 215 510 210 515 200 490
    PCE-DMF 210 495 195 500 190 480
    下载: 导出CSV 
    | 显示表格

    在减水剂掺入量和混凝土组分配比都相同的情况下,PCE-HEA与PCE-HPA在1 h内坍落度损失较小,而掺有PCE-DEI的混凝土坍落度在2 h内基本无损失。

    HEA、HPA、DEI、DMF水解率随时间的变化趋势如图 9所示。HEA和HPA在1 h内的水解率分别为74.40%、68.69%,而和DEI水解速率较为平缓,在2 h内的水解率分别为64.54%、60.12%。丙烯酸酯的水解为可逆反应,但其在碱性条件下的水解较为彻底,从水解反应机理式(1)可以看出,丙烯酸酯在碱性条件下水解时首先形成1个四面体中间体负离子。因此,羰基附近的碳上带有吸电子基可以使四面体阴离子稳定而促进反应;同时空间位阻越小也越利于水解的进行。如果丙烯酸酯中羰基碳上取代基越多,空间位阻越大,则反应速率越慢。丙烯酸羟丙酯比丙烯酸羟乙酯的与氧相连的烷基碳上多1个-CH2-,使得空间位阻增大,故相同时间内释放羧基的数量少于丙烯酸羟乙酯。富马酸二甲酯与衣康酸二甲酯由于其两端各有1个-CH3,降低了水解速度,所以掺有二甲酯的水泥净浆流动度可以达到2 h以上的缓释分散效果,富马酸二甲酯的水解率高于衣康酸二甲酯,但在水泥体系中缓释效果并不明显,可能是因为该单体反应活性与减水剂基础体系不太匹配。

    (1)
    图  9  丙烯酸酯水解率随时间的变化
    Figure  9.  The change of acrylate hydrolysis rate with time

    利用水溶液自由基聚合反应制备4种不同类型的缓释型聚羧酸减水剂PCE。红外光谱、核磁共振氢谱分析均表明:已成功合成出目标产物,实验中所有单体自由基聚合基本反应完全。吸附实验、Zeta电位及羧酸酯单体水解率测试结果表明:含有羧酸单酯PCE的水解速率高于羧酸双酯,随着时间的延长,羧酸单酯型减水剂在1 h内基本水解完全,羧酸双酯型减水剂可以在2 h内持续水解;缓释型聚羧酸减水剂的掺入使得表面张力都有不同程度的减小,水解前后表面张力的对比结果说明亲水性强的羧酸酯单体水解较快。水泥净浆流动度和混凝土实验表明:在聚羧酸减水剂中引入酯类单体可以改善减水剂的保坍性能,而且引入不同的酯类对性能的提升程度不一致。在4种减水剂中,采用DEI合成的减水剂PCE-DEI对混凝土保坍性能的改善效果最优,PCE-DMF次之,说明双酯单体的引入可以使PCE达到2 h以上的缓释分散效果。PCE-HEA和PCE-HPA型减水剂由于含有的酯类单体水解较快,只能在1 h内维持水泥浆体的流动性。

  • 图  1   固定应力分割迭代耦合模型流程图

    Figure  1.   Iterative coupling scheme with fixed-stress splitting

    图  2   验证模型与结果对比

    Figure  2.   Validation model and results comparison

    图  3   地质模型与注入位置

    Figure  3.   Geological model and injection location

    图  4   CO2溶解量演化

    Figure  4.   Evolution of dissolved CO2 mass

    图  5   咸水层顶部距井筒0.5 m处的孔隙压力演化

    Figure  5.   Evolution of pore pressure in the top grid cell of the saline aquifer at a distance of 0.5 m from the wellbore

    图  6   不同热-流-固耦合作用下的位移演化

    Figure  6.   Evolution of displacement considering different thermal-hydraulic-mechanical coupling effects

    表  1   验证模型关键参数

    Table  1   Validation model key parameters

    类别 注入方式
    模型尺寸/m 110×110
    网格数量 36×36
    渗透率k/m2 4.935×10-14
    孔隙度ϕ 0.25
    弹性模量E/Pa 4.5×108
    泊松比v 0
    比奥系数α 1.0
    流体压缩系数cf/Pa 4.0×10-10
    流体黏度μ/(Pa·s) 1.0×10-3
    顶部作用压力及范围 T1=2.0×107 Pa, a1=0~4 m;
    T2=1.0×107 Pa, a2=4~110 m
    初始压强Pi/Pa 1.0×107
    固结系数cv/(m2·s-1) 0.021 3
    下载: 导出CSV

    表  2   模型关键参数

    Table  2   The results of all-or-nothing assignment method

    类别 模型和参数 文献
    咸水层顶部压力/MPa 16.00 [21]
    压力分布 静水压力分布 [23]
    咸水层顶部温度/℃ 55.00 [21]
    温度梯度/(℃·m-1) 0.025 [23]
    咸水层内NaCl质量分数 0.06 [21]
    渗透率(咸水层) kr = 5.912 × 10-15 m2;
    kz=0.1kr
    [21]
    渗透率(非咸水层) kr=8.16 × 10-21 m2;
    kz=0.1kr
    [22]
    孔隙度(咸水层)ϕ 0.126 [21]
    孔隙度(非咸水层)ϕ 0.006 [22]
    岩石热传导率kθd/(W·m-1·K-1) 2.51 [22]
    岩石比热/(J·kg-1·K-1) 920 [22]
    岩石密度/(kg·m-3) 2 260 [22]
    热传导模型 ˉkθ=kθd+ϕβ=A,GSβkθβ
    相对渗透率 krA=S{1[1(S)1λ]λ}2;S=SASirASmxASirA;λ=0.45krG=(1ˆS)2(1^S2);ˆS=SASirA1SirASirG [22]
    残余液体饱和度SirA 0.30 [21]
    残余气体饱和度SirG 0.05 [21]
    毛细管压力 \\Pcap =P0[(S)1λ1]1λ;S=SASitA SASiAA;λ=0.45$:P_{\max }=1.0 \times 10^7 \mathrm{~Pa};$P0=19.61×103 Pa$:P_{\text {max }}=1.0 \times 10^8 \mathrm{~Pa};$P0=6250.0×103 Pa$ [22]
    渗透率折减经验系数 5.0
    孔隙体积压缩系数/(Pa-1) 0.404×10-10 [24]
    弹性模量E/Pa 咸水层:1.0×1010;非咸水层:5.0×109 [22]
    泊松比v 咸水层:0.25;非咸水层:0.25 [22]
    比奥系数α 咸水层:1.0;非咸水层:1.0 [22]
    下载: 导出CSV

    表  3   注入方案

    Table  3   Injection scheme

    编号 名称 注入方式
    1 水气交替 CO2(12 h)、咸水(12 h)
    2 间歇注入 CO2(12 h)、无(12 h)
    3 水气交替 CO2(6 h)、咸水(6 h)、CO2(6 h)、咸水(6 h)
    4 间歇注入 CO2(6 h)、无(6 h)、CO2(6 h)、无(6 h)
    下载: 导出CSV
  • [1] 蔡博峰, 李琦, 张贤, 等. 中国二氧化碳捕集利用与封存(CCUS)年度报告(2021)——中国CCUS路径研究[R]. 生态环境部环境规划院, 中国科学院武汉岩土力学研究所, 中国21世纪议程管理中心, 2021.
    [2] 张贤, 李阳, 马乔, 等. 我国碳捕集利用与封存技术发展研究[J]. 中国工程科学, 2021, 23(6): 70-80. https://www.cnki.com.cn/Article/CJFDTOTAL-GCKX202106008.htm

    ZHANG X, LI Y, MA Q, et al. Development of carbon capture, utilization and storage technology in China[J]. Strategic Study of CAE, 2021, 23(6): 70-80. https://www.cnki.com.cn/Article/CJFDTOTAL-GCKX202106008.htm

    [3] 张贤, 李凯, 马乔, 等. 碳中和目标下CCUS技术发展定位与展望[J]. 中国人口·资源与环境, 2021, 31(9): 29-33. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGRZ202109005.htm

    ZHANG X, LI K, MA Q, et al. Orientation and prospect of CCUS development under carbon neutrality target[J]. China Population, Resources and Environment, 2021, 31(9): 29-33. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGRZ202109005.htm

    [4] 李海峰, 王强. CCUS中CO2利用和地质封存研究[J]. 现代化工, 2022, 42(10): 86-90;95. https://www.cnki.com.cn/Article/CJFDTOTAL-XDHG202210017.htm

    LI H F, WANG Q. Study on utilization and geological storage of CO2 in CCUS[J]. Modern Chemical Industry, 2022, 42(10): 86-90;95. https://www.cnki.com.cn/Article/CJFDTOTAL-XDHG202210017.htm

    [5]

    ZHANG D, SONG J. Mechanisms for geological carbon Sequestration[J]. Procedia IUTAM, 2014, 10: 319-327. doi: 10.1016/j.piutam.2014.01.027

    [6]

    VIVEK R, KUMAR G S. Numerical investigation on effect of varying injection scenario and relative permeability hysteresis on CO2 dissolution in saline aquifer[J]. Environmental Earth Sciences, 2016, 75(16): 1192/1-14.

    [7] 王涛, 于海洋, 朱旭晨, 等. 水气交替CO2咸水层地质封存数值模拟研究[J]. 中国海上油气, 2023, 35(4): 198-204. https://www.cnki.com.cn/Article/CJFDTOTAL-ZHSD202304020.htm

    WANG T, YU H, ZHU X, et al. Numerical simulation study on geological storage of CO2 in saline aquifers assisted by water alternating gas[J]. China Offshore Oil and Gas, 2023, 35(4): 198-204. https://www.cnki.com.cn/Article/CJFDTOTAL-ZHSD202304020.htm

    [8] 王锐, 李阳, 吕成远, 等. 鄂尔多斯盆地深部咸水层CO2驱水与埋存潜力评价方法研究[J]. 非常规油气, 2021, 8(5): 50-55;69. https://www.cnki.com.cn/Article/CJFDTOTAL-FCYQ202105010.htm

    WANG R, LI Y, LV C, et al. Study on potential evaluation method on CO2 EWR and storage for deep saline layers in Ordos Basin[J]. Unconventional Oil & Gas, 2021, 8(5): 50-55;69. https://www.cnki.com.cn/Article/CJFDTOTAL-FCYQ202105010.htm

    [9]

    PAVAN T N V, GOVINDARAJAN S K. Numerical investigations on performance of sc-CO2 sequestration associated with the evolution of porosity and permeability in low permeable saline aquifers[J]. Geoenergy Science and Engineering, 2023, 225: 211681/1-12.

    [10]

    GUDALA M, GOVINDARAJAN S K, YAN B, et al. Comparison of supercritical CO2 with water as geofluid in geothermal reservoirs with numerical investigation using fully coupled thermo-hydro-geomechanical model[J]. Journal of Energy Resources Technology, 2023, 145: 061302/1-19.

    [11] 龚耕, 李毅, 唐栋, 等. THM耦合条件下CO2地质封存注入方案优化研究[J]. 工程地质学报, 2023, 31(3): 1084-1096. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202303029.htm

    GONG G, LI Y, TANG D, et al. Research on optimization of CO2 injection scheme under THM couplings in CO2 geological storage[J]. Journal of Engineering Geology, 2023, 31(3): 1084-1096. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202303029.htm

    [12]

    AGHIGHI M A, ROSHAN H. Coupled modeling[M]//DAYA SAGAR B S, CHENG Q, MCKINLEY J, et al. Encyclopedia of Mathematical Geosciences. Cham: Springer International Publishing, 2020: 1-5.

    [13] 张冠儒. 咸水层碳封存中CO2迁移转化过程的数值模拟研究[D]. 咸阳: 西北农林科技大学, 2016.
    [14]

    MORIDIS G, MORIDIS G J, KOWALSKY M B, et al. TOUGH+ Hydrate v1.0 user's manual: a code for the simulation of system behavior in hydrate-bearing geologic media: LBNL-149E[R]. Lawrence Berkeley National Lab. Berkeley, 2008.

    [15]

    WAGNER W, PRUβ A. The IAPWS formulation 1995 for the thermodynamic properties of ordinary water substance for general and scientific use[J]. Journal of Physical and Chemical Reference Data, 2002, 31(2): 387-535. doi: 10.1063/1.1461829

    [16]

    SPAN R, WAGNER W. A new equation of state for carbon 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. doi: 10.1063/1.555991

    [17]

    ZHANG J, MORIDIS G J, BLASINGAME T A. Message passing interface (MPI) parallelization of iteratively coupled fluid flow and geomechanics codes for the simulation of system behavior in Hydrate-Bearing geologic media. Part 1: Methodology and validation[J]. SPE Reservoir Evaluation & Engineering, 2022, 25(3): 600-620.

    [18]

    ZHANG J, MORIDIS G J, BLASINGAME T A. Message passing interface (MPI) parallelization of iteratively coupled fluid flow and geomechanics codes for the simulation of system behavior in Hydrate-Bearing geologic media. Part 2: Parallel performance and application[J]. SPE Reservoir Evaluation & Engineering, 2022, 25(3): 621-640.

    [19]

    MCNAMEE J, GIBSON R E. Displacement functions and linear transforms applied to diffusion through porous elastic media[J]. The Quarterly Journal of Mechanics and Applied Mathematics, 1960, 13(1): 98-111. doi: 10.1093/qjmam/13.1.98

    [20]

    MCNAMEE J, GIBSON R E. Plane strain and axially symmetric problems of the consolidation of a semi-infinite clay stratum[J]. The Quarterly Journal of Mechanics and Applied Mathematics, 1960, 13(2): 210-227. doi: 10.1093/qjmam/13.2.210

    [21] 喻英, 李义连, 杨国栋, 等. 储层物性参数对CO2长期封存能力的影响研究[J]. 安全与环境工程, 2017, 24(5): 75-83;89. https://www.cnki.com.cn/Article/CJFDTOTAL-KTAQ201705013.htm

    YU Y, LI Y, YANG G, et al. Influence of reservoir physical parameters on the long-term CO2 storage capacity[J]. Safety and Environmental Engineering, 2017, 24(5): 75-83;89. https://www.cnki.com.cn/Article/CJFDTOTAL-KTAQ201705013.htm

    [22] 胡佩. 多场耦合模型在CO2咸水层封存分析中的应用[D]. 大连: 大连理工大学, 2022.
    [23] 曹默雷, 陈建平. CO2深部咸水层封存选址的地质评价[J]. 地质学报, 2022, 96(5): 1868-1882. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202205022.htm

    CAO M L, CHEN J P. The site selection geological evaluation of the CO2 storage of the deep saline aquifer[J]. Acta Geologica Sinica, 2022, 96(5): 1868-1882. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE202205022.htm

    [24] 王福刚, 孙兆军, 刘红艳, 等. 中-细粒长石石英砂岩低渗储层在围压循环增减条件下渗透率变化规律的试验研究[J]. 水利学报, 2016, 47(9): 1125-1132. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201609004.htm

    WANG F, SUN Z, LIU H, et al. Experimental study on the variation of permeability of medium-fine feldspar-quartz sandstone low-permeability reservoir under the circulatory increasing or reducing conditions of confining pressure[J]. Journal of Hydraulic Engineering, 2016, 47(9): 1125-1132. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201609004.htm

  • 期刊类型引用(4)

    1. 郭海超,唐新德,晁亚楠,马菲,王鑫砚. 聚合条件对聚羧酸减水剂缓释行为的影响. 山东化工. 2024(02): 1-4+8 . 百度学术
    2. 王舒,吴朱亮,郑建民,张阳,邓益强. 乙烯氧基聚氧乙烯醚类保坍型聚羧酸减水剂的合成与应用. 华南师范大学学报(自然科学版). 2024(03): 34-39 . 百度学术
    3. 刘少兵,程绍娟,陈森章,林海阳,崔珊珊,杨国武. EPEG缓释型聚羧酸减水剂的常温制备与性能研究. 混凝土与水泥制品. 2023(02): 23-27 . 百度学术
    4. 赵潜. 一种HPEG型多元复合聚羧酸减水剂的合成与性能研究. 广东建材. 2023(05): 10-13+33 . 百度学术

    其他类型引用(1)

图(6)  /  表(3)
计量
  • 文章访问数:  137
  • HTML全文浏览量:  58
  • PDF下载量:  36
  • 被引次数: 5
出版历程
  • 收稿日期:  2024-01-22
  • 网络出版日期:  2024-04-29
  • 刊出日期:  2024-02-24

目录

/

返回文章
返回