4. 尾矿库典型渗流工况建模

尾矿库渗流分析不是只算一个“稳定状态”就结束。同一座库在不同运行阶段，浸润线位置、出逸点、渗流量和排渗设施的减压效果差别很大。建模时通常至少要设置三类典型工况：正常运行水位工况、洪水工况和排渗设施失效工况。三者的几何模型可以共用，但边界条件、材料参数和排渗单元的激活状态必须分别调整，否则计算结果无法对应真实风险。

本节以某上游法尾矿库为例说明设置方法。该库初期坝为透水堆石坝，坝高约35 m，堆积坝最终标高约+120 m，库内设排渗褥垫与竖向排渗井。以下参数均为教学示例的合理假设值，用于说明建模逻辑，不代表实测数据。

一、共用几何模型与材料分区

三类工况应建立在同一套几何模型上，便于结果对比。建模范围取坝体、堆积尾砂、库内沉积滩和坝基覆盖层。材料分区至少包括：尾细砂、尾粉砂、尾黏土、堆石坝体、褥垫排水层和坝基砂砾层。各分区渗透系数按饱和值输入，非饱和部分用土水特征曲线描述。示例取值如下：尾细砂 k=2×10⁻⁵ m/s，尾粉砂 k=5×10⁻⁶ m/s，尾黏土 k=8×10⁻⁸ m/s，堆石坝体 k=1×10⁻² m/s，褥垫排水层 k=5×10⁻³ m/s。尾砂渗透系数具有明显各向异性，水平与垂直方向比值可取3～5，建模时应在材料属性中分别指定 kx、ky。

网格划分上，排渗褥垫、排渗井和坝趾出逸区是渗流梯度变化剧烈的位置，应局部加密。若采用二维断面模型，需明确断面位置代表库长方向的哪一段；若采用三维模型，则排渗井按实际间距布置。无论哪种方式，浸润线观测孔位置都应在模型中标记为监测点，便于后续与实测对比。

二、正常运行水位工况

正常运行水位工况反映库内水位处于设计正常高水位、排渗设施完好时的稳定渗流状态。上游边界设为定水头边界，水头等于正常水位标高，例如+112 m；下游边界设在初期坝趾以外一定距离，取定水头或自由出逸条件。排渗褥垫和排渗井按设计尺寸激活，赋予对应渗透系数。

该工况的核心判断指标是浸润线埋深和坝坡出逸点位置。计算完成后，沿堆积坝外坡提取浸润线与坡面的交点，若出逸点高于坡脚以上2～3 m，或浸润线在坝顶附近与坡面相交，说明排渗能力不足。示例模型中，正常水位下浸润线在堆积坝外坡的出逸点位于坡脚以上约1.2 m，属于可接受范围。此时还应记录单宽渗流量，作为后续工况的对比基准。

三、洪水工况

洪水工况不是简单把上游水头抬高。尾矿库洪水位的持续时间、库水位的骤升过程以及洪水期间排渗设施是否仍正常工作，都会影响孔压响应。建模时通常做两项调整：一是将上游定水头提高至洪水位，例如+116 m；二是根据洪水历时决定采用稳定渗流还是瞬态渗流分析。若洪水持续时间较短、尾砂渗透性较低，孔压来不及充分消散，应采用瞬态分析并输入水位上升曲线。

瞬态分析需要给出初始孔压场，一般以正常运行水位工况的稳定渗流结果作为初始条件，再施加洪水位时间序列。时间步长应覆盖水位上升段和洪峰持续段，例如按1 d步长计算10 d。洪水工况下重点观察坝坡出逸点是否上移、坝趾渗流量是否突增，以及排渗褥垫是否出现局部承压过高。示例中洪水位下出逸点上移至坡脚以上约2.8 m，单宽渗流量约为正常工况的1.6倍，提示需复核坝坡抗渗稳定。

四、排渗设施失效工况

排渗设施失效工况用于评价褥垫堵塞、排渗井淤堵或排水管断裂等情形下的渗流风险。建模时不是删除排渗单元，而是降低其渗透系数或将其改为低渗透材料。例如将褥垫渗透系数由5×10⁻³ m/s降至5×10⁻⁶ m/s，模拟淤堵；或将排渗井按失效比例分段降低导水能力。若需模拟排水管完全失效，可将对应出口边界改为不透水边界。

该工况通常与正常水位或洪水位组合，形成“正常水位+排渗失效”“洪水位+排渗失效”等组合工况。判断重点是浸润线抬升幅度和出逸点是否越过坡面防护层。示例中正常水位下褥垫淤堵后，浸润线整体抬升约4 m，出逸点上移至坡脚以上约5.5 m，已接近坡面反滤层范围，应判定为需要处理的风险工况。

五、参数调整与结果判读的统一做法

三类工况的建模差异集中在三处：上游水头、分析类型和排渗单元状态。建议在模型中建立参数集，把正常水位、洪水位、褥垫渗透系数、排渗井导水系数作为可切换变量，避免每次重建几何。计算完成后统一提取三项指标：坝坡出逸点高程、浸润线在关键断面的埋深、坝趾单宽渗流量。三项指标与正常工况基准值对比，即可判断洪水工况和失效工况的风险增量。

最后需要强调，模型参数应通过现场浸润线观测孔和渗流量监测数据进行反演校核。若正常水位工况计算浸润线与实测偏差超过1 m，应先调整尾砂各向异性比或褥垫渗透系数，再用于洪水与失效工况预测。只有基准工况可信，其余工况的结论才具有工程参考价值。