2. 渗流理论基础与尾矿特性

尾矿库渗流分析的第一步不是打开软件建模，而是把"水在尾矿里怎么流"这件事说清楚。渗流计算的所有结果——浸润线位置、单宽渗流量、出逸坡降——都建立在两个基础之上：一是描述水与土骨架之间宏观运动关系的渗流定律，二是描述尾矿本身透水能力的颗粒级配与渗透系数。本节围绕达西定律的适用边界、非达西渗流的判别方法、级配与渗透系数的定量关系，以及渗流分析中关键参数的选取原则展开。

一、达西定律及其适用条件

达西定律的表达式为 v = k·i，即渗流表观流速与水力坡降成正比，比例系数为渗透系数 k。这里必须明确参照约定：v 是单位土体截面积上的表观流速（又称达西流速），并非孔隙中的真实流速；真实孔隙流速 v_p = v/n，n 为孔隙率。水力坡降 i = Δh/L，Δh 为上下游水头差，L 为渗流路径长度，水头以位置水头加压力水头计，向下渗流时高程降低而水头损失，符号约定全文统一。

达西定律成立的前提是层流占主导。工程上常用雷诺数判别：Re = v·d/ν，其中 d 为颗粒代表粒径，ν 为水的运动黏性系数（20 ℃时约 1.0×10⁻⁶ m²/s）。一般认为 Re 小于约 1～10 的范围内达西定律有效。尾矿砂粒径细、流速低，多数坝体内部渗流处于达西区；但在排渗设施附近、粗粒尾矿堆存区、堆石排渗体内部以及水力坡降较大的局部区域，流速显著增大，可能出现偏离线性的非达西渗流。

二、非达西渗流的判别与处理

非达西渗流主要表现为两种形式：低坡降下的起始坡降现象（细粒尾矿、黏性尾矿中水需克服结合水膜阻力才开始流动），以及高流速下的非线性偏离（惯性力不可忽略，v 与 i 不再成正比）。工程判别可按以下思路进行：先按达西定律反算表观流速与雷诺数，若 Re 明显超出层流上限，或现场注水/抽水试验得到的 v—i 曲线明显弯曲、不通过原点，则应改用非线性模型。

常用的非线性描述为 Forchheimer 型关系：i = a·v + b·v²，其中 a 项对应黏性阻力、b 项对应惯性阻力。当 b 项可忽略时退化为达西定律，a 即 1/k。对存在起始坡降 i₀ 的细粒尾矿，可采用 i = i₀ + v/k 的形式。实际评价中，若排渗体附近按达西计算得到的流速偏高、浸润线偏低，往往提示应考虑非线性阻力，此时渗流量会小于达西估算值，但局部水头损失更集中，需关注排渗体与尾矿接触面的坡降。

三、尾矿颗粒级配与渗透系数的关系

尾矿是选矿后排出的细粒物料，其级配受矿石类型、磨矿细度和排放方式影响很大。级配通过三个途径控制渗透性：细粒含量决定孔隙通道的连通性，不均匀系数 Cu = d₆₀/d₁₀ 反映粒径分布的宽窄，特征粒径 d₁₀、d₆₀ 则直接进入经验公式。工程上常用的估算式包括：

• 太沙基式：k = C·d₁₀²·e³/(1+e)，k 单位 cm/s，d₁₀ 单位 cm，C 取 1～1.5，e 为孔隙比；
• 哈增式：k = c·d₁₀²，适用于砂性尾矿的快速估算，c 为经验系数；
• 当细粒（d<0.075 mm）含量超过约 30%～35% 时，粗颗粒骨架被细粒充填，渗透系数会陡降，此时经验公式误差明显增大，应以室内渗透试验为准。

需要强调，尾矿的渗透系数不是常数。沉积滩面在排放口附近粗、远离排放口细，形成由粗到细的分带；坝前放矿与周边放矿的沉积规律不同；尾矿在堆存过程中还会持续固结，孔隙比减小使 k 下降。因此渗流分析中应按分区分层给定 k，而不是全库取一个值。

四、渗流分析关键参数的选取

参数选取直接决定计算结果的可靠性，应遵循"以试验为主、经验校核、分区赋值、考虑不利工况"的原则。下表给出某尾矿库示例（假设参数，仅用于说明取值方法）的关键参数选取思路：

以示例参数做一次简单校核：若坝前细粒区渗径 L = 20 m，水头差 Δh = 8 m，则 i = 0.4，按 k = 2×10⁻⁶ cm/s = 2×10⁻⁸ m/s 计算，表观流速 v = k·i = 8×10⁻⁹ m/s，单宽渗流量 q = v·H（H 为过水断面高度）。这一量级说明细粒尾矿本身透水性极弱，渗流控制的关键往往不在坝体内部，而在排渗设施能否有效导出渗水、降低浸润线。反过来，若排渗体渗透系数取值偏大而忽略淤堵折减，计算会高估排渗能力、低估浸润线，这是评价中常见的安全隐患。

综上，达西定律给出了渗流计算的基本框架，非达西判别决定了该框架在何处需要修正，级配与渗透系数的关系提供了参数取值的物理依据，而关键参数的合理选取与反演校核，则是把理论转化为可靠评价结论的最后一环。三者环环相扣，缺一不可。