当CO₂变成‘数字流体’：相场法解锁压裂仿真新维度

【引言】

在页岩油开采的竞技场上，水力压裂技术曾被誉为“解锁地下黑金的钥匙”，但随之而来的水资源消耗、化学污染和低效裂缝预测等问题，正让这把钥匙逐渐生锈。当全球能源行业将目光投向更清洁的超临界CO₂压裂技术时，一个更棘手的难题浮出水面：如何驯服这种介于气液之间的“暴躁流体”，精准预判它在千米地层下的裂岩轨迹？

传统实验手段如同“盲人摸象”——物理模拟成本高昂，现场试错风险巨大，而经典数值模型又难以刻画CO₂与页岩间复杂的相变交互。直到相场法（Phase-Field Method）的出现，这场博弈迎来了转机。这项起源于材料科学的数学工具，正将超临界CO₂转化为可被方程描述的“数字流体”，在虚拟空间中重构裂缝生长的每一个细节：从CO₂分子穿透岩石孔隙的微观动力学，到宏观裂缝网络的混沌分形演化，原本不可见的流体暴力被解构成万亿次计算的优雅舞蹈。

当算法的精度突破物理实验的边界，页岩压裂正在从“经验驱动”迈向“预测驱动”的新纪元。这场由相场法引领的仿真革命，或将重新定义非常规油气开采的底层逻辑——用数字孪生代替盲目试错，用计算预见性取代经验不确定性。而我们，正站在这场技术范式转移的临界点上。

【超临界CO2压裂模拟存在难点】

难点一：CO₂的“流体变身记”——超临界态下的物性捉迷藏

超临界CO₂既非气体也非液体，密度接近液体，黏度类似气体，扩散性极强。这种“双重人格”导致它在页岩孔隙中的流动规律难以捉摸。传统本构方程（如Forchheimer方程）在水力压裂中表现稳定，但面对超临界CO₂时，其渗透率-压力关系会因微小的温压波动（±5℃或±0.5MPa）而发生剧变，就像在仿真中埋下无数个“数值地雷”。更棘手的是，CO₂与页岩有机质接触时可能引发萃取效应（如沥青质溶解），进一步改变岩石力学参数——流体在改造岩石，岩石也在反向驯化流体。

难点二：裂缝演化的“混沌预言”——多物理场耦合陷阱

CO₂压裂并非单纯的流体驱动裂缝问题，而是热-流-固-化（THMC）四场耦合的混沌系统：

*热场：CO₂注入时的焦耳-汤姆逊效应会引发近井地带温度骤降（可能达30℃），导致岩石脆性增强；

*化学场：CO₂-地层水-页岩的化学反应（如碳酸盐化）会动态改变孔隙结构和裂缝表面能；

*流固耦合：超临界CO₂的低黏度特性使其更易渗入纳米级孔隙，诱发“裂缝分支”（多级次生裂缝无序增殖）。

难点三：相场法的“精度围城”——从方程到现实的最后一公里 尽管相场法通过序参量（Phase-Field Parameter）优雅地统一了裂缝形核、扩展与交汇的数学描述，但其在CO₂压裂中的落地仍面临三重围剿：

*参数：CO₂-页岩界面能、裂缝表面耗散率等关键参数在高温高压下极难通过实验获取，许多模型被迫使用“借来的数据”（如借用甲烷或水的数据）；

*尺度：纳米孔隙（~10nm）中的CO₂分子吸附效应与宏观裂缝（~10m）的扩展需跨10个数量级建模，常规网格剖分技术如同“用渔网捞病毒”；

*验证困局：地下压裂过程无法直接观测，即便仿真结果与微地震监测数据吻合，也可能只是“错误的胜利”（不同机理导致相似表象）。

【采用COMSOL模拟实现】

采用COMSOL实现相场法模拟裂缝扩展过程，考虑以下几个因素：

（1）超临界CO2物性随压力变化，包括流体粘度和密度；

（2）超临界CO2流动能力强，很容易进入孔缝之中，改变岩体有效应力，减弱岩体的剪切强度，因此，需要考虑拉伸和剪切强度变化，建立相场驱动方程，从而实现多分支裂缝扩展形态；

（3）进一步，可以分析各类因素的综合影响，包括：天然裂缝、地应力、纹层以及多簇裂缝应力干扰作用等。

采用滑溜水压裂，求解获得简单裂缝形态

采用超临界CO2压裂，求解获得多分支裂缝形态

https://mp.weixin.qq.com/s/i42LhN_9amAAOlMHBc7lTA