深埋地下厂房涌水COMSOL仿真
关键词：COMSOL；三维渗流；地下厂房；涌水量预测；参数反演
一、文章简要介绍
深埋地下水电站修在地下几百米，最大的风险之一就是涌水。这篇文章以一个深埋地下厂房为背景，把地质模型和水流模型耦合起来，用 COMSOL 模拟了开挖前、开挖期和运行期三个阶段的渗流场，并预测了开挖期的涌水量。工程背景很典型：厂房紧邻河道和一条冲沟，场区里分布着多条断层和一条韧性剪切带，构造控水特征明显。我们重点看它的 COMSOL 建模与参数反演部分（DOI：10.3390/w18091000）。
二、COMSOL 建模与参数反演
模型建在 COMSOL Multiphysics 6.0 里，物理内核是三维达西渗流，求解非均质参数场下的连续性方程和达西定律。边界条件分三类：河道和冲沟按水位插值给定水头，并用地形梯度构造随位置变化的函数；地表和模型上部边界设成质量流边界，用来施加降雨入渗；模型底面和基岩设成不透水边界。网格用 COMSOL 自带的流体动力学模块生成加密四面体网格，最终 68.8 万个单元、11.9 万个节点，兼顾了精度和计算量。
参数怎么来？作者先把场区按岩性和构造分成多个水力分区，用压水试验和抽水试验给出渗透系数的初值区间：基岩只有 10⁻⁶ 量级，强卸荷带和几条断层能到 10⁻⁴ 量级，差了近三个数量级。再用 MATLAB 2024b 写了机器学习模型，在区间内生成渗透系数，通过 MATLAB-COMSOL 实时链接反复迭代，用 5 个钻孔的实测水位做约束，最终 RMSE 只有 8.9 米，R² 达到 0.916，相对 2000 多米的总水头误差不到 0.4%。
三、关键结果与发现
先看开挖前的天然渗流场。场区地下水整体从西北流向东南，水头向西北方向梯度变大。三种降雨工况下水平流场都比较稳定，但降雨越强，垂向渗流的路子变化越明显；暴雨时地表高水头会把水压向地下，剖面上甚至出现局部的地下水分水岭。
开挖期是重点。厂房开挖后地下水位明显下降，形成降落漏斗，最大降幅约 200 米。作者把厂房分成球阀室、主厂房、主变洞和尾水管闸室四段，分别预测涌水量：无雨时储水区最大，约 10583 立方米每天；暴雨时升到 13043.7，是全场最大值；而增幅最大的是主变洞，暴雨下增加 49.68%。设了排水廊道之后，各段涌水量明显回落，暴雨下储水区的增幅被压到 40.9%。
运行期换成算外水压力。衬砌和灌浆把洞室当成不透水边界，水头逐步恢复。沿输水线路统计，暴雨工况下主厂房到储水区段的外水压力最大，达到 882.78 米水头，比无雨状态增加 43.7%；下游输水管的增速最大，达到 58.9%。这说明管线设计必须按极端降雨留足压力余量。
