3. 尾矿库渗流模型构建流程

尾矿库渗流模型不是把勘察报告里的地层线原样搬进软件就能算对的。它本质上是一次从真实地质体到数学计算域的等价映射：把连续、复杂、非均质的三维地质体，压缩成一组可求解的网格、材料参数和边界条件。映射过程中每一步的取舍，都会直接反映在浸润线位置、渗流量和排渗设施效果的判断上。本节按建模的实际推进顺序，讲清地质剖面概化、计算域离散、边界条件设置、初始条件与材料分区这五个环节各自要做什么、怎么做、做到什么程度算合格。

一、地质剖面概化：先决定“算几层、算多大”

概化的第一步是确定建模范围。以一座典型上游法尾矿库为例：初期坝为透水堆石坝，坝高约 30 m，堆积坝最终标高比初期坝顶高约 40 m，库区一侧为弱透水基岩，另一侧为第四系覆盖层。建模范围通常取坝轴线向上游至库尾、向下游至排渗设施影响范围以外，两侧以相对隔水边界或分水岭为界。范围过小会截断渗流场，导致下游溢出量偏小；范围过大则网格量激增、参数不确定性放大。

剖面概化的核心是地层合并与简化。勘察报告往往给出十余个亚层，建模时按渗透性量级归并：渗透系数相差不足一个数量级的相邻土层可合并为一层，取厚度加权平均渗透系数；对连续性好、厚度稳定的关键弱透水层（如坝基淤泥质黏土层）必须单独保留，因为它控制着坝基渗流和浸润线抬升。概化结果应形成一张“概化地质剖面图”，标注各层顶底高程、产状和代表性渗透系数，作为后续材料分区的依据。

二、计算域离散：网格要“贴着结构走”

离散是把概化剖面切成有限元或有限差分网格。尾矿库渗流模型的网格有两个硬要求：一是排渗设施附近必须加密，排渗盲沟、褥垫、竖向排渗井的尺寸通常只有 1～3 m，若网格尺度大于设施尺寸，设施在模型中会被“抹平”，排渗效果无法体现；二是自由面附近需要足够分辨率，因为浸润线是迭代求解出来的，网格过粗会导致自由面振荡或收敛困难。

实际操作中可采用“全局粗、局部细”的非结构化网格：坝体和堆积坝区域单元尺寸取 2～5 m，排渗设施周边加密至 0.3～0.5 m，远场区域放宽至 10～20 m。离散完成后必须检查网格质量：单元畸变过大、出现尖角或负体积单元，都会导致求解失败或结果失真。一个可执行的检查是统计最小单元内角与边长比，并对排渗设施轮廓做一次网格贴合性目视核对。

三、边界条件设置：位置和类型都要有依据

渗流模型的边界条件主要有三类：定水头边界、定流量边界和隔水边界。尾矿库模型中最关键的是上游库水位边界和下游溢出边界。库水位边界施加在堆积坝上游坡面与库底，水头值取设计洪水位或正常水位（示例中可取正常水位 45 m、洪水位 47 m 两种工况分别计算）；下游坡面脚部设置为潜在溢出边界（seepage face），即允许水在坡面逸出、但逸出点水头等于位置高程，由迭代确定。

常见错误是把下游坡面直接设为定水头或全隔水。设为定水头会人为压住浸润线，算出偏安全的结果；设为全隔水则水无处可去，自由面会异常抬升。坝基底部和两侧远场边界一般取隔水边界，但若基岩存在区域性补给，应按勘察给出的水头或流量折算为定水头/定流量边界。边界条件设置完成后，应逐条核对：每条边界的位置、类型、数值来源是否能在勘察或设计文件中找到对应依据。

四、初始条件与材料分区：让迭代“有起点、有归属”

稳态渗流计算对初始条件不敏感，但非稳态（如库水位骤降、排渗设施失效过程）必须给出合理初始孔压场。工程上常用做法是先算一个稳态场作为初始条件：以正常水位运行工况收敛后的孔压分布，作为骤降或故障工况的起点。这样既避免了“全零初始场”带来的长时间瞬态，也使后续结果具有明确的物理起点。

材料分区是把概化剖面上的每一层、每一类结构映射为计算域中的材料编号，并赋予渗透系数。尾矿砂、坝体堆石、排渗褥垫、弱透水层应分别编号。参数取值优先采用勘察试验值；缺乏试验时可按经验范围取值并在报告中注明。示例参数如下（属教学合理假设，非实测）：

分区建模时要注意材料界面的连续性：相邻材料区的公共边必须共享节点，否则会出现“界面漏水”或“界面阻水”的数值假象。赋参后建议做一次量级检查——用达西定律估算单宽渗流量 q≈k·i·A，与同类工程经验值对比，若相差两三个数量级，多半是单位或分区映射出了错。

五、建模完成后的自检闭环

五个环节走完后，不要急于看云图。先做三项自检：网格与设施是否贴合（排渗设施是否被足够细的单元覆盖）、边界是否闭合（计算域四周每条边都有明确边界类型）、材料分区是否全覆盖（不存在未赋参的“空白单元”）。三项通过后，再以正常水位工况试算，观察浸润线是否从库水位平滑过渡到下游溢出点、渗流量是否与估算量级一致。若浸润线在排渗设施处出现不合理的陡降或平台，通常说明该处网格或渗透系数赋值需要复核。至此，一个可用于排渗设施评价的渗流模型才算构建完成。