关于CO2地下盐水层封存的Comsol复刻之旅

【模型】:CO2地下盐水层封存comsol二氧化碳封存 【软件】:COMSOL 【类型】:复刻某sciQ2区-关于低渗透盐水含水层中与孔隙度和渗透率演变相关的超临界二氧化碳封存性能的数值研究 【物理场景】:研究了二氧化碳在低渗透和高渗透盐水含水层中渗透率模量的影响考虑温度孔隙度、渗透率、注入速度、渗透率各向异性、残余水饱和度和残余气体饱和度的影响 【复现效果】:基本复现效果如下对比图 【延伸方向】:后续可考虑CO2在盐水中的低溶解性以及不同气体驱替效果 【难点】:复现过程中出现体积分数的奇异值已解决最近在研究CO2地下盐水层封存相关课题用Comsol软件复刻了一篇发表在sci Q2区的论文主题是“低渗透盐水含水层中与孔隙度和渗透率演变相关的超临界二氧化碳封存性能的数值研究”和大家分享下过程与心得。一、物理场景剖析这篇论文聚焦在研究二氧化碳在低渗透和高渗透盐水含水层中渗透率模量的影响。在这个过程中需要考虑众多因素如温度、孔隙度、渗透率、注入速度、渗透率各向异性、残余水饱和度和残余气体饱和度等。每个因素都对CO2封存效果有着微妙且关键的作用。比如温度它会影响CO2在盐水中的溶解度以及流体的粘性等物理性质。在Comsol中我们可以通过定义热传递模块来考虑温度场的变化。假设我们有一个简单的二维模型定义温度场的代码片段可能如下// 定义温度场 model.geom.create(geom1, 2); model.physics.create(ht, HeatTransferInSolids, geom1); model.physics(ht).par(T0, 300[K]); // 初始温度 model.physics(ht).bc(bc1).set(T, T0);这里通过HeatTransferInSolids模块来模拟固体中的热传递设定了初始温度T0并在边界条件bc1中应用。温度的变化会连锁反应般影响到孔隙度和渗透率等其他参数比如温度升高可能导致岩石膨胀进而改变孔隙结构影响渗透率。【模型】:CO2地下盐水层封存comsol二氧化碳封存 【软件】:COMSOL 【类型】:复刻某sciQ2区-关于低渗透盐水含水层中与孔隙度和渗透率演变相关的超临界二氧化碳封存性能的数值研究 【物理场景】:研究了二氧化碳在低渗透和高渗透盐水含水层中渗透率模量的影响考虑温度孔隙度、渗透率、注入速度、渗透率各向异性、残余水饱和度和残余气体饱和度的影响 【复现效果】:基本复现效果如下对比图 【延伸方向】:后续可考虑CO2在盐水中的低溶解性以及不同气体驱替效果 【难点】:复现过程中出现体积分数的奇异值已解决孔隙度和渗透率更是核心参数。孔隙度决定了岩石能容纳多少流体而渗透率则决定了流体在岩石中的流动难易程度。在Comsol里定义渗透率张量时可以考虑各向异性的情况代码示例model.physics.create(darcy, DarcyFlow, geom1); model.physics(darcy).prop(k, [1e - 12, 0; 0, 1e - 13]); // 各向异性渗透率张量上述代码创建了DarcyFlow模块来模拟渗流并定义了一个简单的各向异性渗透率张量x方向渗透率为1e - 12$m^2$y方向为1e - 13$m^2$。注入速度则直接影响着CO2在含水层中的运移速度和分布范围调整注入速度参数就能观察到不同的CO2扩散形态。二、复现效果及对比经过一番努力基本复现了论文中的效果。通过对比图可以清晰看到模型模拟出的CO2在含水层中的分布与论文中的结果趋势相符。这意味着我们在参数设置、物理场耦合等方面的理解和操作是正确的。这种对比验证不仅是对自己工作的肯定也是进一步深入研究的基础。三、攻克难点体积分数奇异值问题在复现过程中遇到了体积分数出现奇异值的问题。这可把我折腾了好一阵子。经过一番排查发现是在多相流模型中某些参数的边界条件设置不合理导致的。比如说在定义残余水饱和度和残余气体饱和度时如果设置的值与实际物理情况偏差较大就可能引发数值计算的不稳定出现奇异值。解决办法是重新审视这些参数的物理意义和取值范围结合实际地质数据进行调整。同时在Comsol的多相流模块中对体积分数的计算方法和限制条件进行了仔细检查和优化。例如在设置多相流的相间传质时确保传质系数等参数的取值在合理区间内避免因传质过快或过慢导致体积分数计算异常。model.physics.create(mp, MultiphaseFlow, geom1); model.physics(mp).set(residual_saturation_w, 0.2); // 合理设置残余水饱和度 model.physics(mp).set(residual_saturation_g, 0.1); // 合理设置残余气体饱和度通过这样调整残余饱和度等关键参数成功解决了体积分数奇异值的问题让模拟能够顺利进行。四、延伸方向展望后续研究可以考虑CO2在盐水中的低溶解性这一特性。CO2在盐水中并非无限溶解了解其溶解极限以及溶解过程对封存效果的长期影响至关重要。此外研究不同气体驱替效果也很有意义。比如使用氮气等其他气体进行驱替对比与CO2驱替效果的差异这对于优化地下封存方案有着重要的实际应用价值。或许在未来的研究中能找到更高效、更安全的CO2地下封存方式为应对全球气候变化贡献一份力量。希望这次的分享能给同样在这个领域探索的小伙伴们一些启发一起加油