基于COMSOL软件的相场方法模拟多孔介质两相驱替及拓展研究：附算例与毛管数、饱和度计算方法

提供基于comsol中相场方法模拟多孔介质两相驱替（水气、油水等等）的算例（也可以定做水平集驱替的算例），可在此基础上学会利用comsol软件进行两相流驱替的模拟，拓展研究，具体参考算例附后。 附赠基于相场方法模拟驱替时的毛管数计算方法和饱和度计算方法

咱们今天唠唠怎么用COMSOL的相场法整两相驱替仿真。搞石油的兄弟都知道，水驱油、气驱水这些场景里，界面运动跟毛细血管似的，传统方法处理这种弯弯绕绕的界面能累死个人。相场法这玩意儿妙就妙在不用实时追踪界面，直接靠浓度场过渡，特别适合新手操作。

先看模型骨架——在COMSOL里新建多物理场，勾选层流、相场和达西定律这三个模块。材料库调个水的黏度参数，多孔介质渗透率建议先设1e-12 m²试水。这里有个关键点：相场参数里的界面厚度得手动调，太薄了容易发散，太厚了界面糊成一片。老司机经验是初始网格用常规尺寸，在界面附近手动加密三层。

看这段相场控制的PDE代码：

phi_t + u·grad(phi) = gamma*(epsilon*laplacian(phi) - (phi*(1-phi)*(1-2*phi))/epsilon + 2*lambda*epsilon*phi*(1-phi)*grad(c))

这堆符号看着唬人，其实核心就是控制相场变量phi在0到1之间平滑过渡。gamma是迁移率，控制相界面移动速度，新手建议先设1e-3。epsilon是界面厚度系数，通常取网格最小尺寸的1/5。重点注意最后那个grad(c)项，这是耦合毛管力的关键，漏了这项毛细现象直接摆烂。

提供基于comsol中相场方法模拟多孔介质两相驱替（水气、油水等等）的算例（也可以定做水平集驱替的算例），可在此基础上学会利用comsol软件进行两相流驱替的模拟，拓展研究，具体参考算例附后。 附赠基于相场方法模拟驱替时的毛管数计算方法和饱和度计算方法

驱替过程的后处理有讲究。饱和度计算直接在派生变量里写：

S_w = (phi > 0.5) ? 1 : 0; //二值化处理 integrate(S_w, pore_volume)/total_volume //孔隙度修正

实际跑仿真时会发现界面处有过渡区，建议用时间平均法消抖。毛管数Ca的计算要抓取界面曲率，在COMSOL里直接用内置的pfm.curvature变量，配合表面张力系数sigma：

Ca = (mu * velocity_magnitude) / (sigma * abs(pfm.curvature))

注意曲率正负号代表界面凹凸方向，做统计时要取绝对值。有个坑点：COMSOL的曲率计算在三维和二维中差个系数，二维模型记得手动乘2。

调试时最容易崩的是相场与流动的耦合顺序。建议先冻结流速场跑稳相场，再开启双向耦合。碰到发散别急着调步长，把相场方程的gamma调小一个量级试试。有次我仿水驱油，界面跑着跑着突然分叉，最后发现是入口流速设太大导致毛细数超临界值，把注入速度砍半就正常了。

想进阶的话可以玩点花的：在达西定律里加Forchheimer修正项模拟高速流动，或者在相场方程里掺入随机扰动模拟介质非均质性。最近帮个课题组做了个裂缝性储层的版本，关键是在几何里随机生成裂缝网络，然后给裂缝和基质的相场参数设不同权重系数，出来的指进现象那叫一个带感。

水平集方法其实也能搞，但相场更适合多孔介质这种复杂边界。两者切换时注意重新标定表面张力参数，水平集里的δ函数宽度得和相场的epsilon对应起来。有个取巧办法：用同一组几何模型，把相场结果导出作为水平集的初始条件，对比两种方法的指进形态差异。