三相流仿真：CO₂地质封存中毛细管捕获的数值透视

三相流仿真在地质封存工程评估中处于一个无法被回避的位置——当超临界CO₂注入含咸水砂岩地层后，CO₂-咸水-岩石三相系统的毛细管捕获效率直接决定了封存容量的估算上限。我们参与过一个深部咸水层CO₂封存可行性评估项目，核心任务是通过三相流仿真预测十年注入期内CO₂的运移范围和捕获机制分配比例，整个工程论证的底层逻辑就压在这套仿真的可信度上。

项目在哪里碰壁

第一轮三相流仿真采用TOUGH2软件，网格分辨率50米，注入井模型给出的CO₂气相饱和度分布在三年后出现了非物理振荡——注入井附近的网格饱和度跳变超过0.4，明显是网格分辨率不足以捕捉毛细管压力前沿。项目在这里卡了三周，因为饱和度振荡会导致相对渗透率插值失真，进而影响整个运移预测。更棘手的是，粗网格下毛细管捕获份额被严重低估，结构捕获占比仅为8%，而文献中同等条件下的实测值通常在15-25%之间。

方法的性格：为什么选择OpenFOAM而非TOUGH2

认定OpenFOAM配合VOF方法更适合这个场景的高分辨率需求，理由并不复杂。TOUGH2的优势在于热力学一致性，但它的积分有限差分格式在处理三相界面时的数值弥散难以控制。OpenFOAM的VOF方法通过几何重构追踪界面，在网格足够密时能够清晰分辨CO₂-咸水界面，毛细管压力的梯度计算精度显著优于有限差分方案。

代价是计算规模——OpenFOAM模型采用非结构化网格，近井区域加密至0.5米分辨率，远场区域渐变至20米，总网格数达到380万。三相流仿真的时间步长受CFL条件约束，注入初期需要控制在0.01天以内，十年模拟累计超过36万个时间步，在192核集群上运行了11天。关于VOF方法的理论框架，[OpenFOAM官方文档](https://www.openfoam.com/documentation/guides/latest/doc/guide-applications-solvers-multiphase.html)提供了多相流求解器的详细说明，而[CFD Online](https://www.cfd-online.com/Wiki/Volumetric_particle_tracking)上关于界面重构方法的讨论也有参考价值。

相对渗透率曲线的赋值是另一个决定性的参数。项目组采用Corey模型，气相指数n_g=2.5，液相指数n_l=3.0，残余气饱和度S_gr=0.15，端点相对渗透率k_rg0=0.35。这些参数来自目标地层的实际岩芯驱替实验，而非文献默认值——因为不同砂岩的孔喉分布差异极大，直接套用文献曲线的三相流仿真结果几乎不可信。

数据的态度：毛管数揭示捕获机制切换

OpenFOAM模型给出的十年预测结果中，毛细管捕获占比升至18%，溶解捕获占比9%，结构捕获占比3%，与文献区间吻合。但真正有价值的发现藏在毛管数（Nc = μv/σ）的敏感性分析里——当注入速率从0.1 kg/s提升到0.5 kg/s时，毛管数从10⁻⁷升至10⁻⁶，毛细管捕获份额从18%下降到11%。差距不会说谎：注入速率越高，CO₂越容易突破毛细管束缚，运移到更远区域被结构捕获或溶解。

这个结论对工程设计的意义直接而明确：降低单井注入速率、增加注入井数量，虽然在投资上不经济，但毛细管捕获效率的提升可以显著延长封存安全边际。项目组据此向甲方建议了五井分布注入方案替代三井集中注入方案。

回过头看的判断

项目最终采用OpenFOAM高分辨率模型支撑了封存容量评估报告，十年预测的CO₂运移范围与后续监测井的压力响应数据基本一致，算是被部分验证了的结论。值得警醒的是，如果项目组满足于TOUGH2粗网格的结果，毛细管捕获被低估40%的结论将直接误导工程设计决策。三相流仿真的价值不在于跑出一个漂亮的饱和度云图，而在于它能否在网格分辨率和物理参数上做到足够诚实——这两点决定了结论的工程可信度。

对于需要解析多相流体在多孔介质中运移行为的工程评估项目，三相流仿真提供了从机理理解到工程优化的完整数值链路。只要在网格分辨率、界面追踪方法和本构关系参数上保持审慎，它的结论足以支撑重大工程决策。