COMSOL岩石酸化模型：碳酸钙与氧化钙的随机溶解与布林克曼流动

COMSOL岩石酸化模型 碳酸钙氧化钙遇酸溶解，孔隙度随机，酸化路径随机，布林克曼流动，形成雪花状路径。

实验室的玻璃器皿里正上演一场微观地质剧变——当30%盐酸注入碳酸钙基岩，白色固体表面瞬间腾起细密气泡。这看似普通的化学反应，在孔隙尺度下其实藏着复杂的流动-溶解耦合效应。我们今天用COMSOL复现这个岩酸作用过程，看看如何让代码流淌出雪花状的侵蚀路径。

先给模型打地基：创建三维几何时特意用随机函数生成孔隙网络。这里有个小技巧，用高斯随机场定义初始孔隙度分布，0.3的方差让岩石像瑞士奶酪般布满大小不一的孔洞：

% 孔隙度生成脚本 phi0 = 0.2 + 0.3*randn(100,100); phi0(phi0<0.01) = 0.01; //防止负值 surf(phi0,'EdgeColor','none');

运行后得到起伏的孔隙度地形图，像被陨石撞击过的月球表面。这种非均匀结构将引导酸液优先侵蚀脆弱区域。

流动模块选择布林克曼方程是关键——它兼顾了达西流的孔隙渗透和自由流的惯性效应。当设置动态粘度参数时，记得关联实时变化的孔隙度：

// 布林克曼方程参数设置 model.physics('br').set('eta', '0.001*(1+0.5*sin(t))'); //动态粘度 model.physics('br').set('alpha', '1/phi^2'); //渗透率衰减

这种非线性关系会导致流速在侵蚀通道内突然加快，就像峡谷中的河流遇到狭窄处会湍急。

COMSOL岩石酸化模型 碳酸钙氧化钙遇酸溶解，孔隙度随机，酸化路径随机，布林克曼流动，形成雪花状路径。

化学反应模块配置碳酸钙溶解动力学时，采用表面反应速率模型：

reaction_rate = k * (1 - C/C_sat)**n * surf_area

其中表面面积surf_area随孔隙演化动态更新。当某处孔隙扩大到临界值时，程序自动触发网格重构——这相当于岩石在代码世界里实现了自我撕裂。

求解器跑完三小时后，后处理窗口出现令人惊艳的雪花状侵蚀网络。仔细观察分支角度，发现主通道与次级分支呈72°左右夹角，这与现场酸化压裂的CT扫描结果惊人相似。原来酸液在突破孔隙瓶颈后，会沿最大主应力方向分叉，形成类似冰晶生长的形态。

调试过程中有个意外发现：当孔隙度变化步长设为0.05时，模型会出现蝴蝶效应式的路径分异。某次将随机种子从42改为43，原本向右延伸的主通道突然向左突围——这或许解释了现实中相距仅半米的注酸井为何会产生截然不同的增产效果。