深埋地下厂房涌水COMSOL仿真

关键词：COMSOL；三维渗流；地下厂房；涌水量预测；参数反演

一、文章简要介绍

深埋地下水电站修在地下几百米，最大的风险之一就是涌水。这篇文章以一个深埋地下厂房为背景，把地质模型和水流模型耦合起来，用 COMSOL 模拟了开挖前、开挖期和运行期三个阶段的渗流场，并预测了开挖期的涌水量。工程背景很典型：厂房紧邻河道和一条冲沟，场区里分布着多条断层和一条韧性剪切带，构造控水特征明显。我们重点看它的 COMSOL 建模与参数反演部分（DOI：10.3390/w18091000）。

图1 三维模型与边界条件设置（COMSOL）

二、COMSOL 建模与参数反演

模型建在 COMSOL Multiphysics 6.0 里，物理内核是三维达西渗流，求解非均质参数场下的连续性方程和达西定律。边界条件分三类：河道和冲沟按水位插值给定水头，并用地形梯度构造随位置变化的函数；地表和模型上部边界设成质量流边界，用来施加降雨入渗；模型底面和基岩设成不透水边界。网格用 COMSOL 自带的流体动力学模块生成加密四面体网格，最终 68.8 万个单元、11.9 万个节点，兼顾了精度和计算量。

图2 研究区三维地质模型

图3 模型网格划分（68.8 万四面体单元）

参数怎么来？作者先把场区按岩性和构造分成多个水力分区，用压水试验和抽水试验给出渗透系数的初值区间：基岩只有 10⁻⁶ 量级，强卸荷带和几条断层能到 10⁻⁴ 量级，差了近三个数量级。再用 MATLAB 2024b 写了机器学习模型，在区间内生成渗透系数，通过 MATLAB-COMSOL 实时链接反复迭代，用 5 个钻孔的实测水位做约束，最终 RMSE 只有 8.9 米，R² 达到 0.916，相对 2000 多米的总水头误差不到 0.4%。

图4 钻孔位置分布（5 个水位观测孔）

三、关键结果与发现

先看开挖前的天然渗流场。场区地下水整体从西北流向东南，水头向西北方向梯度变大。三种降雨工况下水平流场都比较稳定，但降雨越强，垂向渗流的路子变化越明显；暴雨时地表高水头会把水压向地下，剖面上甚至出现局部的地下水分水岭。

图5 天然状态下不同降雨工况的渗流场云图

开挖期是重点。厂房开挖后地下水位明显下降，形成降落漏斗，最大降幅约 200 米。作者把厂房分成球阀室、主厂房、主变洞和尾水管闸室四段，分别预测涌水量：无雨时储水区最大，约 10583 立方米每天；暴雨时升到 13043.7，是全场最大值；而增幅最大的是主变洞，暴雨下增加 49.68%。设了排水廊道之后，各段涌水量明显回落，暴雨下储水区的增幅被压到 40.9%。

图6 开挖期渗流场云图（平面与剖面）

运行期换成算外水压力。衬砌和灌浆把洞室当成不透水边界，水头逐步恢复。沿输水线路统计，暴雨工况下主厂房到储水区段的外水压力最大，达到 882.78 米水头，比无雨状态增加 43.7%；下游输水管的增速最大，达到 58.9%。这说明管线设计必须按极端降雨留足压力余量。

图7 运行期渗流场云图（平面与剖面）

四、我们提供的服务

如果你在做地下厂房、隧道或边坡的渗流仿真，我们可以提供 COMSOL 三维渗流建模、地质分区与边界条件设置、网格划分与无关性验证、渗透系数反演与模型标定、降雨工况对比分析、涌水量预测，以及外水压力评估等服务。如需相关服务，欢迎私信聊聊。

如有需要，欢迎通过公众号“320科技工作室”联系我们。