多孔介质两相驱替模拟：相场法原理与COMSOL实践

1. 多孔介质两相驱替模拟的工程挑战

在油气开采、地下水污染治理等领域，多孔介质中的两相驱替过程模拟一直是工程仿真中的难点。当非混相流体在多孔介质中发生驱替时，界面不稳定性会导致复杂的指进现象——就像把水倒入干燥的沙堆时，液体不会均匀渗透，而是形成树枝状的分叉通道。这种粘性指进现象会显著影响驱替效率，在三次采油中可能导致高达30%的采收率差异。

传统模拟方法（如VOF）在处理复杂界面拓扑变化时面临巨大挑战。相场法通过引入扩散界面概念，将尖锐界面转化为厚度可控的过渡区域，配合Cahn-Hilliard方程描述相场演化，成为处理这类问题的有力工具。COMSOL Multiphysics提供的相场模块，让我们能够相对便捷地实现这一物理过程的仿真。

2. 相场法建模核心架构

2.1 Cahn-Hilliard方程解析

相场法的核心控制方程包含三个关键部分：

matlab复制// 双阱势函数（Double-well potential）
theta = epsilon^2*phi*(1-phi)*(1-2*phi)  

// 化学势定义
mu = theta - epsilon^2*laplacian(phi)    

// 相场演化方程
d(phi,t) + div(phi*u) = mobility*laplacian(mu)  

其中界面厚度参数epsilon的选取需要满足：

code复制epsilon ≥ 2.5*(网格尺寸)

在实际操作中，我们通常取网格尺寸的1/5~1/10。例如当使用0.1mm的基础网格时，epsilon建议设置在0.02~0.05mm范围。过大的epsilon会导致界面过度模糊，而过小则会引起数值不稳定。

迁移率mobility的物理意义是相场演化的速率系数。根据我们的实测经验：

• 对于水驱油系统：1e-6 ~ 1e-5 [m³·s/kg]
• 对于气驱水系统：1e-7 ~ 1e-6 [m³·s/kg]

2.2 多孔介质特性建模

多孔介质的非均质性是指进现象的主要诱因。我们采用随机场方法生成孔隙率分布：

matlab复制// 孔隙率随机场生成（保证可重复性）
random_seed = 12345;  // 固定随机种子
phi_porous = 0.6 + 0.2*r