comsol 实际孔隙 两油扩散 驱替老张蹲在实验室盯着CT扫描的岩心图片两个油罐子歪在墙角——这已经是第三次驱替实验失败了。孔隙结构像被揉皱的报纸手工建模根本还原不了真实孔喉比。我咬着冰棍凑过去要不试试COMSOL的Level Set方法带孔隙结构导入的那种。导入真实孔隙模型得先处理STL文件。COMSOL的CAD模块总喜欢把微小孔隙识别成噪点这时候就需要点几何修复的魔法importGeometry(core_sample.stl); // 修复孔隙连接处的微小缝隙 model.geom(geom1).scale(1e-6); // 单位转换μm→m model.geom(geom1).run(repair);操作台的咖啡渍旁边老张的笔记本记着两种原油参数密度892kg/m³ vs 856kg/m³界面张力0.025N/m。在相场模块里设置双相流别忘了迁移方程里的奇偶稳定性修正Physics.create(pf, PhaseField); Physics.feature(pfinit).set(c0, 0.5*(1-tanh((x-1e-4)/2e-5))); //初始界面锋面 Physics.feature(pfeq1).set(Gamma, 0.032*[1,0.8]); //各向异性参数网格划分最要命。孔隙边缘用边界层网格捕捉毛细管力主体区域上自由四面体。有次偷懒全用四面体结果计算到80%直接发散——油水界面像破抹布一样碎成渣。comsol 实际孔隙 两油扩散 驱替等计算结果时顺手调出瞬态求解器的诊断图。当局部Peclet数超过2马上切换全耦合求解器。记得打开自动时间步进的限制器不然驱替前缘容易跳步Study.create(std1, Transient); Study.feature(time).set(tlist, range(0,5,100)); Study.feature(solver).set(stol, 1e-4);后处理时发现角落孔隙存在滞油区。用流线积分器追踪粒子轨迹果然有漩涡生成。老张突然拍大腿注采压差设低了毛细管数卡在临界值下边了调出参数化扫描结果对比滞油体积从23%直降到7%。窗外知了叫得正欢屏幕上两相界面优雅地推进。老张摸着三天没刮的胡子早该上多物理场耦合的白瞎三罐原油。我默默保存mph文件——鬼知道服务器下次啥时候崩。