相场法在水力压裂模拟中越来越火，尤其是COMSOL这种多物理场耦合神器。今天咱们拆解几个典型工况，手把手看裂缝怎么在代码里“长“出来。先拿最简单的单裂缝开刀——

COMSOL 相场法与水力压裂 案例一：单一裂缝延伸； 案例二：两簇压裂； 案例三：三簇压裂-对称； 案例四：三簇压裂-完全； 案例五：水力裂缝与垂直天然裂缝相交； 案例六：水力裂缝与倾斜裂缝相交。

案例一的核心代码段长这样：

PhaseField.epsilon = 0.05; //界面厚度参数 PhaseField.Gc = 2.7e3; //断裂能 physics.set('d', 1, 'DispersionType', 'Anisotropic');

epsilon控制相场界面过渡区的宽度，这个值太小会导致计算爆炸，太大又会让裂缝变"胖"。Gc是断裂能，相当于材料的抗裂能力，这里设2700J/m²对应页岩参数。各向异性分散设置让裂缝沿着最大主应力方向劈开。

进阶到案例二两簇压裂时，射孔位置的初始化很关键：

initCrack1 = (x>=2.5).*(x<=3).*(y<=0.1); initCrack2 = (x>=7).*(x<=7.5).*(y<=0.1); model.phase = 1 - sqrt(initCrack1.^2 + initCrack2.^2);

用两个矩形区域定义初始裂缝簇，注意平方开根号的写法比直接相加更稳定。当泵注压力上来后，两簇裂缝会出现竞争扩展，这时候看应力阴影效应特别明显——先起裂的裂缝会把周围应力场抬高，压制另一簇发展。COMSOL的自动自适应网格这时会疯狂在裂缝尖端加密，记得在求解器配置里打开误差估计器。

COMSOL 相场法与水力压裂 案例一：单一裂缝延伸； 案例二：两簇压裂； 案例三：三簇压裂-对称； 案例四：三簇压裂-完全； 案例五：水力裂缝与垂直天然裂缝相交； 案例六：水力裂缝与倾斜裂缝相交。

案例四的三簇非对称压裂更有意思，需要处理裂缝转向问题。下面这段流固耦合设置是精髓：

physics.create('FSI', 'LaminarFlow', 'SolidMechanics'); physics.set('p', 'outlet', 'PumpPressure', 'sin(t/10)*1e6'); interface.set('DisplacementContinuity', 'on');

时变泵压用正弦函数模拟现实中的段塞注入，位移连续性条件确保流体压力能传递到固体变形。当中间裂缝比两侧早0.5秒起裂时，会产生不对称的应力场，导致两侧裂缝向外弯曲成"八"字形。这时候相场变量的演化方程里，驱动力的各向异性系数要调高到1.2以上才能捕捉这种转向。

遇到天然裂缝（案例五、六）时，接触条件处理是难关。倾斜裂缝相交的代码魔改：

ContactPair1.set('PhaseFieldThreshold', 0.7); SolidMechanics.add('Friction', 'mu', 0.6);

设置相场阈值0.7作为接触判据，相当于当损伤变量超过70%时启动摩擦接触。摩擦系数mu设为0.6对应砂岩表面粗糙度。有意思的是，当水力裂缝与倾斜天然缝以55度角相交时，会出现三种模式：直接穿透、沿天然缝滑移、或者产生分岔裂缝，这取决于注入速率和水平应力差。在结果后处理中，用等值线叠加粒子流线图，能清晰看到流体优先选择低应力通道的特征。

跑这些案例时有个坑：相场法对网格尺寸敏感得一匹！建议在裂缝路径上预定义扫掠网格，边界层至少分三层。时间步长用自动调整，别头铁用固定步长，否则分分钟发散。最后说个骚操作——把相场变量导出到Paraview，用GLSL着色器实时渲染裂缝生长过程，比自带的后处理器酷炫十倍。