44、COMSOL模拟二维裂隙流压裂水平井裂缝性油藏离散裂缝网络模型COMSOL数值模拟案例
压裂水平井的裂缝网络模拟总让人头大，尤其是面对复杂地质条件的时候。最近用COMSOL折腾了个二维离散裂缝模型，发现几个实战技巧特别有意思，尤其是如何用代码生成随机裂隙——这事儿要是纯靠鼠标点界面得点到腱鞘炎吧。
咱们先画个二维正方形当油藏基底。重点来了，裂缝得随机分布对吧？这时候直接上COMSOL的JAVA API搞个脚本：
double density = 0.1; //裂缝密度 int numFrac = (int)(density * width * height); for (int i=0; i<numFrac; i++) { double x = width * Math.random(); double y = height * Math.random(); double angle = 2*Math.PI * Math.random(); double length = 5 + 10*Math.random(); model.geom().create("frac"+i, "Line"); model.geom("frac"+i).set("x", x + length*Math.cos(angle)); model.geom("frac"+i).set("y", y + length*Math.sin(angle)); }这段骚操作直接在模型里生成了随机位置、随机方向、长度5-15米的裂缝。特别要注意的是density参数，实际项目里得根据岩芯数据调整，搞大了计算量爆炸，搞小了模型失真。
物理场设置才是重头戏。裂隙流用达西定律，但主岩和裂缝的渗透率能差三个数量级。这时候千万别手滑把单位搞错——见过有人把裂缝渗透率设成10e-12 m²结果水流得比高速公路还快，直接变消防栓现场。
边界条件设置有个坑：水平井筒别直接画成线！得用点源边界条件，配合井筒半径参数。见过同行用线源导致压力梯度异常，结果模拟出个喷泉效果，油没采到反而模拟出个趵突泉。
网格划分必须上手动控制。裂缝周围用边界层网格，主岩区域用自由三角形。有个取巧办法是用特征尺寸函数：
% 裂缝附近加密 model.mesh("mesh1").feature("size").set("custom", "on"); model.mesh("mesh1").feature("size").set("hcurve", "0.1+0.5*(z<0.2)");这个表达式让裂缝周围网格尺寸骤降到0.1米，其他地方放宽到0.5米。计算时间直接砍半，精度还能保持。不过要注意过渡区设置，突然的网格尺寸变化会让求解器抽风。
后处理阶段别光盯着压力云图。流速矢量图里藏着大秘密——有次发现某条裂缝流速异常，追查发现是渗透率张量方向设反了，活生生把泄流通道搞成拦路虎。流量积分时记得选对边界，有哥们儿把积分区域选成整个模型域，结果累计产量比油田储量还高，这就很魔幻了。
最后说个血泪教训：瞬态模拟时时间步长别迷信默认设置。特别是压裂液刚开始注入的阶段，用自适应步长容易漏掉压力突变点。手动设置前0.1秒用0.01秒步长，后面逐步放开，能捕捉到明显的压力突破信号。
搞完这套模型最大的成就感，是看着压力云图里那些红色高压区像血管一样顺着裂缝网络蔓延——虽然只是仿真，但那种地下流体被驯服的感觉，比玩任何策略游戏都带劲。
