流网渗透流量模拟器 返回
土壤工程

流网渗透流量模拟器

这是一个用流线网(流网)求解通过大坝或导板墙下的渗透流的工具。改变水位差、透水系数、流路数和等势降数,可以实时计算地盘中的渗透流量、出口动水坡度和管涌安全系数。

参数设置
上下游水位差 H
m
上游和下游水面高度的差异
透水系数 k
m/s
地盘的透水性。砂土约为 10⁻³~10⁻⁵ 级别
流路数 Nf
被流线包围的带(流路)的数量
等势降数 Nd
被等势线分隔的段(降)的数量
构筑物深度
m
垂直于渗透流方向的构筑物长度
出口流场长度 l
m
出口最后一个单元格流动方向的长度
计算结果
单位宽度渗透流量 (m³/s/m)
全部渗透流量 (L/s)
单位降水头损失 Δh (m)
出口动水坡度 i_exit
管涌安全系数
管涌判定
导板下的渗透流和流网 — 动画

显示绕过导板的流线(蓝色)和与其直角的等势线(橙色)。粒子沿流线渗透,在下游出口区域(红色)动水坡度最陡。

全部渗透流量 vs 水位差 H
管涌安全系数 vs 出口流场长度 l
理论·主要公式

$$q=k\,H\,\frac{N_f}{N_d},\qquad i_{exit}=\frac{H/N_d}{l},\qquad FOS=\frac{i_{cr}}{i_{exit}}$$

单位宽度渗透流量 q、出口动水坡度 i_exit、管涌安全系数 FOS。k:透水系数,H:水位差,Nf:流路数,Nd:等势降数,l:出口流场长度。

Nf 是流网中被流线包围的流路数,Nd 是被等势线分隔的段数。当出口动水坡度 i_exit 接近临界动水坡度 i_cr(砂土约为 1)时,管涌危险增加。本工具使用代表值 i_cr = 1.0。

什么是流网

🙋
"流网"就是地盘里水流的网格状图吧?那为什么要画这个呢?
🎓
是的,当地下水从大坝或导板墙下渗过时,会画成网格状的图形。地盘中的渗透流遵循"拉普拉斯方程",这是一个偏微分方程。在没有计算机的年代,需要用手工求解。所以人们想出了一种方法,让流线和等势线相互垂直,形成"正方形"的网格。这就是流网法。只要画出图来,就能读出渗透流量和危险位置。
🙋
流线和等势线分别代表什么?
🎓
流线就是水粒子从上游流向下游的路径。等势线是全水头相等的点的连线。全水头是"水压+位置势能"。这两族线总是成直角交叉,把渗透区域分成"正方形"的曲线方格。左边的图中,蓝色线是流线,橙色线是等势线。用滑块改变 Nf 和 Nd,网格的细密度就会改变。
🙋
明白了。那渗透流量怎么算出来的?公式是 q = k·H·(Nf/Nd)…
🎓
画出流网后,只需数出两个整数:流路数 Nf 和等势降数 Nd。Nf 是被流线包围的带的本数,Nd 是被等势线分隔的段数。然后用 q = k·H·(Nf/Nd) 计算单位宽度的渗透流量。乘以构筑物的深度就得到全部的漏水量。比如默认值的情况,每降下降时水头失失 0.5m,全体漏水量是 0.6 L/s。
🙋
不只是漏水量,还有"安全系数"对吧。这个是什么的安全系数?
🎓
这是防止管涌的安全系数。看下游的出口附近,等势线紧密地集中在一起,动水坡度就变陡了。当出口的动水坡度 i_exit 接近临界动水坡度(砂土约为 1)时,向上的渗透力会超过砂粒的有效重量,砂粒就会浮起来,地盘会像"沸腾"一样破坏。这就是管涌,也叫做沸腾。实际上很多大坝决口和开挖事故都是从这里开始的。所以我们计算 FOS = i_cr/i_exit,通常要确保 FOS ≥ 3。
🙋
当安全系数不足的时候,怎么对策?
🎓
最直接的办法就是"加深导板"。这样水的通路就会变长,等势降数增加,出口的坡度就会变缓。或者在下游铺上导流毯来增加出口的流场长度,或者用井点降低地下水位。看左边"出口流场长度 l"的参数,增大这个值,安全系数会显著上升。流网是用一张图同时告诉我们漏水量和"最危险的地方"的非常强大的工具。

常见问题

流网是绘制流线和等势线直角相交形成"正方形"网格的图。单位宽度的渗透流量用 q = k·H·(Nf/Nd) 计算。k 是透水系数,H 是上下游水位差,Nf 是流路数(被流线包围的带数),Nd 是等势降数(被等势线分隔的段数)。Nf 和 Nd 只需绘制流网后数出即可,乘以构筑物的深度就得到全部渗透流量。
在下游出口附近,等势线密集集中,动水坡度变陡。当出口动水坡度 i_exit 接近临界动水坡度 i_cr(砂土约为 1.0)时,向上的渗透力超过砂粒的有效重量,地层浮起形成"沸腾"般的破坏。这就是管涌或沸腾。许多大坝和开挖事故都源于此处。流网用一张图清楚地显示等势线最密集、最危险的位置。
单位降水头损失为 Δh = H/Nd。设出口流场长度为 l,则出口动水坡度为 i_exit = Δh/l = (H/Nd)/l。管涌安全系数为 FOS = i_cr/i_exit。通常确保 FOS ≥ 3,当安全系数不足时,采取加深导板、铺设导流毯或井点降水等措施来增加出口流场长度。
流网法是用手工求解二维拉普拉斯方程的经典方法,在没有计算机的时代就已广泛应用。如今可用渗流有限元法进行详细分析,但流网法作为概算和验证工具仍然有效。对于边界条件简单、地盘等向均质的情况,手工绘制的流网精度足以满足工程要求。如果有限元结果与流网概算相差很大,应怀疑网格划分或透水系数、边界条件的输入错误。

实际应用

土坝·堆石坝设计:用土或岩石筑造的大坝中,堤体和基础地盘中通过的渗透流用流网评估。既要抑制漏水量,也要防止下游坡脚出现管涌,保持动水坡度在临界值以下。坝心(防渗区)的厚度和基础的防渗墙(截水)的深度根据流网中能获得足够的等势降数来确定。

导板截水和开挖工事:在河川和港湾工事中,用导板挡水并在内侧进行开挖时,绕过导板的渗透流可能导致开挖底面发生"沸腾"破坏。设计时需用流网检讨导板的根入深度,确认出口动水坡度的安全系数。当安全系数不足时,需加深导板或用井点强制降低地下水位。

堤·水门·取水构造物:混凝土堤或水门的基础地盘中,从上游到下游的渗透流会产生扬压力(托起构筑物的水压)。从流网的等势线可以读出基础底面的扬压力分布,用于结构物的滑动、倾覆稳定性计算。导流毯和截水墙的配置是控制渗透路径的设计手段。

CAE渗流分析的验证:在用 SEEP/W 或 PLAXIS 等渗流有限元进行详细分析之前或之后,用流网手算来确认 q 和 i_exit 的数量级。如果有限元结果与流网概算相差很大,说明可能存在透水系数单位错误或边界条件(水头边界、不透水边界)设置错误。简单的概算是复杂数值分析的可信性基础。

常见误解和注意点

首先,有个很大的误解是"流网的正方形必须是真正的正方形"。流网的方格实际上是曲线正方形,只要对边的中点长度大致相等、对角线互相直交就足够了。不必为了追求完美的正方形而反复重画。即使流路数 Nf 和等势降数 Nd 有小数点(比如 3.5),也没问题。重要的是比值 Nf/Nd,即使网格粗糙,这个比值也能稳定地求出来,这正是流网法的优势。

其次,"直接把等向均质地盘的公式用于异向地盘"。q = k·H·(Nf/Nd) 的公式是基于水平和竖直透水系数相等的等向地盘。实际的堆积地盘通常水平透水性强于竖直(kh > kv),具有成层异向性。这种情况需要在"变换断面"上绘制流网,即把水平距离缩小 √(kv/kh) 倍,流量计算用等效透水系数 k' = √(kh·kv)。忽视异向性会导致漏水量和动水坡度偏离很大。

最后,"只要出口动水坡度的安全系数足够就能防止管涌"的误解。本工具计算的是出口面的渗透力与土的重量比较的"隆起(沸腾)"型安全系数。但实际的管涌还有"后退侵蚀型",即土粒从下游逆向被冲蚀,形成扩大的管状水渗通道。这种情况可能从局部弱点、粒度不连续或动物洞穴等处开始,即使平均动水坡度安全也会发生。所以必须根据粒度分布设置滤层材料,并持续监测渗水的浑浊度和流量。

使用指南

  1. 输入水头差(水位差):输入大坝上下游水位差或导板两侧的水位差,单位为米。例如:上游EL.100m,下游EL.98m时输入2m
  2. 设置透水系数:根据钻孔试验结果输入地盘的透水系数。粘土地盘k=1×10⁻⁷cm/s,砂质地盘k=1×10⁻²cm/s等,单位为m/s
  3. 指定流网网格数:输入等势线数(Nf)和流线数(Nd)。调整 Nf/Nd 比值到 0.3~0.5 左右,反映渗透路径的复杂程度
  4. 运行模拟:检查各输出值,确认出口动水坡度不超过 Terzaghi 的临界动水坡度 ic≈0.5~1.0

具体计算示例

堤体下粘土基础的渗透流分析:水头差Δh=3m,透水系数k=2×10⁻⁶cm/s,流网设置Nf=5、Nd=15的情况下,单位降水头损失0.6m,单位宽度流量q=1.8×10⁻⁵m³/s/m,出口动水坡度i_exit=0.35。在比重γ_w=9.8kN/m³、有效应力σ'_v=150kPa的条件下,管涌安全系数Fs=0.35/0.5≈1.43(安全)。

实务中的注意点

  1. 透水系数的可靠性直接影响结果:使用现场定水位透水试验或钻孔试验得到的k值。估计值有±1个数量级的不确定性,需要用上下限值进行敏感性分析
  2. 确认出口动水坡度不超过 Terzaghi 的临界值 ic:河堤安全标准要求 Fs=ic/i_exit≥1.3。如发现管涌迹象,应考虑延长防渗膜或引入深井降水
  3. 流网网格精度与计算时间的平衡:Nf/Nd=0.33 左右是标准设置。对复杂地形或复合地层,需要更细的网格划分