土的渗透性和渗流问题

官方地理老师·34 页·深入(追求细节与边界)·0 次浏览·2 天前
渗透系数渗流破坏达西定律工程安全

土的渗透性和渗流问题

看懂渗流如何决定工程安全,掌握从达西定律到渗透破坏的完整链条

按 空格/→ 演示下一步

1 / 34 页

全部页面点击任意一页,跳回舞台从这页播放

渗透系数渗流破坏达西定律工程安全

土的渗透性和渗流问题

看懂渗流如何决定工程安全,掌握从达西定律到渗透破坏的完整链条

1第 1 页 · 土的渗透性和渗流问题

什么是渗流

下雨后地面不会瞬间干透——水慢慢往下渗进土里。河边挖一个浅坑,坑底很快会渗出水来。这种「水在土体内部缓慢流动」的现象,就是渗流。

渗流的定义
水在土体内部连通孔隙通道中的缓慢流动,是水与土骨架耦合的过程
发生条件
两个条件缺一不可:存在水头差(势能差),且孔隙相互连通
流动特征
流速极小(mm/s 量级)、连续、稳定,绝大多数情况下为层流
咖啡滴滤对应 →土中渗流

咖啡颗粒间有成千上万条细小水通道,水靠重力(水头差)从上向下慢慢渗过

2第 2 页 · 什么是渗流

渗流问题的工程背景

渗流本身只是水的流动,但工程里「水的流动」往往意味着麻烦:大坝会漏、基坑会涌水、边坡会滑动——三类真实场景背后是同一个核心矛盾。

水利工程
土石坝蓄水后水从坝体渗出,可能带走细颗粒甚至溃坝
基坑工程
开挖到地下水位以下,地下水涌入基坑,需要降水或止水
边坡工程
雨水或地下水位上升使土体饱和,抗剪强度降低,诱发滑坡
核心矛盾
流量问题(水多水少)与稳定问题(渗透破坏)必须同时解决
船底裂缝进水对应 →土体渗透破坏

水从薄弱处渗入并带走材料;土中水渗出会冲走细颗粒、掏空结构

3第 3 页 · 渗流问题的工程背景

渗流引起的工程病害

流土、管涌、地面沉降的典型实例

渗流引起的工程病害
流土、管涌、地面沉降的典型实例
4第 4 页 · 渗流引起的工程病害

为什么要研究渗流

上一页我们看到了管涌、流土、滑坡、基坑失稳……这些灾害提醒我们:渗流绝非小事。面对它,工程界不能只做事后补救,必须建立系统能力——预测、控制、防治,正是研究渗流的根本出发点。

预测——识别风险
通过计算求出孔压、水力坡降和渗流量,在事故发生前圈出危险区
控制——干预渗流
用防渗墙截渗、井点排水、反滤层保护出口,让渗流按设计意图发展
防治——源头规避
选址避让、参数留余量、施工严管控,从设计阶段消除灾害条件
完整应对链
预测定位风险,控制处理当下,防治避免重现,三者缺一不可
医生看病的流程对应 →渗流研究的逻辑

诊断对应预测,手术用药对应控制,保健防复发对应防治

5第 5 页 · 为什么要研究渗流

达西定律的物理意义

上一节我们知道了土里确实有水在流。那水在土里流得有多快?什么决定快慢?达西用一个实验回答了:流速与水头差成正比,土越密实流得越慢。

线性正比关系
渗流速度 v 与水力梯度 i 成正比:水头差越大,渗流越快
水头是驱动力
i = Δh/L 是单位距离的水头损失,即推动水流的'势能坡'
土骨架是阻力
渗透系数 k 综合反映土颗粒对水流的阻碍,越密实的土 k 越小
v = k·i
流速 = 渗透系数 × 水力梯度,驱动力与阻力共同决定
水过滤网对应 →水过土骨架

水压是推力,滤网越细水流越慢;土越密 k 越小水越慢

v=ki=kΔhLv = k \cdot i = k \cdot \dfrac{\Delta h}{L}
6第 6 页 · 达西定律的物理意义

达西实验的原理

沿箭头看水流路径:上方水箱稳住水头,土样筒是研究对象,量筒记录渗出水量。

图解渲染中…
a2水箱水位始终保持h不变,是实验的关键控制点a3土样两端用透水石夹紧,长度L、截面积A已知a4多进水就溢走,保证水头h恒定不变a8k=QL/(Aht),由测得的Q、L、A、h、t反算
7第 7 页 · 达西实验的原理

达西定律的公式与适用条件

上一页达西实验告诉我们:流量与水头差成正比、与渗径成反比。现在把这些观察写成方程,问清楚它在什么条件下成立、什么时候会失效。

达西定律公式
v = ki,渗流速度 = 渗透系数 × 水力梯度;也写为 Q = kiA 求流量
三个变量的含义
v 是单位面积流速,k 是土本身的透水能力,i = Δh/L 是驱动势
由实验到方程
Q ∝ Δh,Q ∝ 1/L,比例常数除以面积 A 即得 v = ki
适用范围
层流状态,主要适用于砂、粉土、黏性土等细粒土,Reynolds 数约小于 1~10
失效边界
粗粒土(砾石、堆石)流速过大转捩为紊流,v–i 关系偏离线性,达西定律不再适用
挤过拥挤走廊对应 →层流与紊流的分界

人多时只能随大流匀速走(v 与推力 i 成正比);人少能冲刺时,速度不再随推力线性增长

v=kiQ=kiAv = ki \qquad Q = kiA
8第 8 页 · 达西定律的公式与适用条件

渗透系数的定义

上一节我们写出达西定律 v=ki,k是个比例系数——但它物理上到底是什么?为什么单位恰好是速度?这页我们拆开看。

比例系数
k是达西定律里的比例系数,反映土体允许水通过的难易程度:k越大水越易渗过
物理意义
令水力坡降 i=1,则 v=k——k在数值上等于「单位水力坡降下的渗透流速」
量纲分析
[k]=[v]/[i]=(L/T)/1=L/T,与速度同量纲,常用 cm/s 或 m/s
依赖因素
k不仅取决于土骨架(孔隙比、粒径),还取决于水的黏滞系数——同种土水温升高,k也变大
电路中的电导率 σ对应 →渗透系数 k

欧姆定律 J=σE 与 达西定律 v=ki 形式完全对应,都是「通量=传导系数×驱动力」

k=vi,[k]=LT1k = \dfrac{v}{i}, \quad [k] = L \cdot T^{-1}
9第 9 页 · 渗透系数的定义

达西定律适用性的判断

雷诺数是否低于临界值,是达西定律线性关系成立与否的分水岭。

层流:达西定律适用
  • 流体质点沿固定路径有序流动
  • 雷诺数 Re 低于临界值(约 1~10)
  • 水头损失与流速呈线性关系
  • 达西定律严格成立
紊流:达西定律失效
  • 流体质点相互混掺、运动紊乱
  • 雷诺数 Re 超过临界值
  • 水头损失与流速呈非线性关系
  • 需改用 Forchheimer 等非线性公式
用雷诺数裁决:Re 低于临界值为层流,达西定律成立;超过后进入紊流,线性关系失效。
10第 10 页 · 达西定律适用性的判断

渗透系数的测定方法

测渗透系数没有万能仪器,要先看土的颗粒粗细,再选对应试验方法。

1
判土选法
粗粒土流量大易测用常水头;细粒土流量极微改用变水头,把小流量放大为水位变化
2
常水头法
维持水头差H恒定,量筒测时段t内通过土样的累计流量Q
3
变水头法
细管初始水头h₁,记录水头下降至h₂对应的时间t
4
公式求k
常水头k=QL/(AHt);变水头k=(aL)/(At)·ln(h₁/h₂)
11第 11 页 · 渗透系数的测定方法

影响渗透系数的因素

土性、颗粒级配、孔隙比、温度的综合影响

影响渗透系数的因素
土性、颗粒级配、孔隙比、温度的综合影响
12第 12 页 · 影响渗透系数的因素

各类土的渗透性差异

砂土、黏土、裂隙岩体的数量级对比

各类土的渗透性差异
砂土、黏土、裂隙岩体的数量级对比
13第 13 页 · 各类土的渗透性差异

渗流场的基本概念

达西定律告诉我们流量和水力坡降成正比,这个「水力坡降」本质就是两点之间的水头差。

势能 ↔ 水头
单位体积水的势能除以 γw,量纲由能量变为长度,得到的就是水头
位置水头 z
该点到选定基准面的高差,对应位能的大小
压力水头 u/γw
孔隙水压力换算成的水柱高度,对应压能的大小
水头差 Δh
上、下游两点的总水头之差,是渗流的根本驱动力
插入土层中的测压管对应 →总水头的物理图像

管中水柱升到某高度后稳定,这个高度恰好就是该点的总水头

h=z+uγwh = z + \dfrac{u}{\gamma_w}
14第 14 页 · 渗流场的基本概念

连续性方程的推导

达西定律告诉我们水流有多快,但渗流场里水从哪儿来、到哪儿去?必须满足质量守恒——这就是连续性方程要做的事。

代表性微元体
在渗流场中取边长 dx·dy·dz 的小立方体作为研究对象
渗流速度分量
沿 x、y、z 三个方向各有渗流速度 vx、vy、vz
三个方向的流量差
每个方向流入量减流出量,等于该方向上水质量的变化率
连续性方程
三项差累加,对稳态不可压缩水流即得连续性方程
家庭账本收支平衡对应 →微元体的质量守恒

微元体就是一本账,每方向入流是收入、出流是支出,收支相抵就是连续性方程

vxx+vyy+vzz=0\frac{\partial v_x}{\partial x} + \frac{\partial v_y}{\partial y} + \frac{\partial v_z}{\partial z} = 0
15第 15 页 · 连续性方程的推导

达西定律与连续性方程的结合

把达西定律代入连续性方程,就能推导出饱和土渗流的基本控制方程。

1
写出达西定律
用流速v与水力坡降i的关系描述渗流
2
取微元体
在饱和土中取一个微小平行六面体单元
3
建立连续性方程
根据质量守恒,进流量等于出流量加上储存量变化
4
代入并化简
把达西定律代入连续性方程,整理得到渗流基本方程
5
得到控制方程
稳态时为Laplace方程,非稳态时为扩散方程
16第 16 页 · 达西定律与连续性方程的结合

拉普拉斯方程的物理意义

上一页把连续性方程和达西定律合到一起,得到了稳定渗流的控制方程——拉普拉斯方程。它为什么能描述整个渗流场?我们来拆开看它的物理含义。

方程形式
∇²h = 0,即水头 h 对各空间坐标的二阶偏导数之和恒等于零
质量守恒
本质是稳态下每个微元体「流入 = 流出」,水量不随时间积累
无源无汇
渗流场内部不存在水的产生或消失,水量只能从边界进出
调和函数
水头 h 是调和函数,内部无极值,最值只能出现在边界上
边值定解
方程本身有无穷多解,必须配合边界条件才有唯一渗流场
稳态温度场对应 →稳定渗流场

温度↔水头、热流密度↔达西流速,稳态下满足同一方程

2h=2hx2+2hy2+2hz2=0\nabla^2 h = \frac{\partial^2 h}{\partial x^2} + \frac{\partial^2 h}{\partial y^2} + \frac{\partial^2 h}{\partial z^2} = 0
17第 17 页 · 拉普拉斯方程的物理意义

什么是流网

上一页讲到拉普拉斯方程描述了渗流场中水头的分布,但要画出整个水头场太麻烦。工程师用两组曲线把渗流场"画"出来,这就是流网。

流线
水粒子实际走过的路径,每点切线方向即流速方向
等势线
水头(势能)相等的点连成的曲线
流网
流线与等势线交织形成的网格
正交性
流线与等势线在每个交点处都垂直
曲边正方形
各向同性土中每个网格单元近似对边中距相等
地形图等高线与山谷溪流对应 →等势线与流线

等高线相当于势能相等的位置,溪流路径就是水流真实走向;二者正交,等高线越密代表坡度越陡

18第 18 页 · 什么是流网

流网的特性与绘制规则

上页讲到流网由流线与等势线交织,但「随手画两条」不算流网——它有严格骨架:正交、曲边正方形、还有边界条件。这一页拆开看。

正交性
流线与等势线处处正交,夹角恒为 90°
曲边正方形
网格由四条曲边围成,相邻两边近似相等,像被捏过的方格
等流量等水头
相邻流线间渗流量相等;相邻等势线间水头差相等
边界条件
不透水边界为流线;进、出水面为等势线
等高线地形图对应 →流网

等高线对应等势线,最陡下降路径对应流线,故二者必然垂直

19第 19 页 · 流网的特性与绘制规则

流网的绘制步骤

从定边界到画网格再到算渗流量,流网绘制是一套需要反复修正的手艺。

1
定边界条件
画出不透水边界和自由水面,标出进出水口位置
2
试绘流线
凭经验画几条流线,初步划分出若干流道
3
画等势线
在流线之间画等水头线,每条都垂直于流线
4
检查正方形
观察每个网格是否近似正方形,不行就调整流线
5
迭代修正
反复调整流线和等势线,直到网格都接近正方形
6
数Nf与Nd
数出流槽数Nf和等势带数Nd,据此算渗流量
20第 20 页 · 流网的绘制步骤

流网绘制的数值方法

有限差分法求解拉普拉斯方程

流网绘制的数值方法
有限差分法求解拉普拉斯方程
21第 21 页 · 流网绘制的数值方法

流网的应用实例

从堤坝剖面出发,沿流程看流网如何一步步产出水头信息并用于工程决策。

图解渲染中…
G下游坡面逸出处的水力坡降,最危险位置Hi允由土的抗渗强度决定,详见流土/管涌判别
22第 22 页 · 流网的应用实例

从流网求渗流量

画好流网后,工程师真正想知道的只有一件事:水每秒能渗出多少?这页给出一个简洁却严谨的算法。

单宽流量公式
q=k·H·Nf/Nd,单位时间单位宽度渗过的水量
流道等流量
每个完整流道携带的流量Δq相等,与流道宽窄无关
等势线等水头差
相邻等势线间水头差Δh=H/Nd,沿整个流场均匀
四步计算流程
画流网→数Nf、Nd→取k、H→代入公式求q
餐厅并行的传菜传送带对应 →流网中的多个完整流道

每条送菜量相等对应每个流道Δq相等;菜从厨房到餐桌的高差对应总水头H

q=kHNfNdq = k \cdot H \cdot \frac{N_f}{N_d}
23第 23 页 · 从流网求渗流量

渗透力的概念

上一节我们用流网算出了渗流量——知道了水「流多少」。但工程上更要紧的另一问:水在土里流过时对土颗粒施加了多大的力?这个力一旦大到能托起土粒,灾难就发生了——这就是渗透力。

渗透力的定义
渗流施加在土骨架上的拖曳力,由水流通过摩擦与冲击传递给土颗粒
体积力性质
单位体积力 j,方向沿流线,与渗流方向一致
渗透力公式
j = i·γw,i 为水力坡降,γw 为水的重度
临界坡降
渗透力等于有效自重时发生流土:icr = (Gs−1)/(1+e)
对有效应力的影响
向上的渗透力相当于加大浮力,使有效应力下降,土骨架变弱
河水冲刷河床沙粒对应 →渗流对土颗粒的拖曳

流速越快(坡降越大),沙粒被冲走的力越大;与 j = i·γw 同源

j=iγwj = i\gamma_w
24第 24 页 · 渗透力的概念

渗透力的计算

上一页我们说渗透力是水拖拽土颗粒的力。这个力到底多大?方向又朝哪?这一页用一个最简单的力学模型把它算出来。

力学模型的建立
取饱和土柱单元,水头差 Δh 驱动渗流,水对骨架产生拖拽
驱动力的推导
水柱两端水压力差 = γw·Δh·A,沿渗流方向作用于水
渗透力公式
单位体积渗透力 j = γw·i,其中 i = Δh/L 为水力坡降
方向判断
渗透力方向与渗流方向一致;向上渗流时易引发流土破坏
物理含义
j 是体积力(kN/m³),可与土的有效重度比较判稳
河水沿坡推移砂砾对应 →渗流沿水头差拖拽土骨架

坡降对应 i;水重度对应 γw;推移力对应渗透力 j

j=γwij = \gamma_w \cdot i
25第 25 页 · 渗透力的计算

流土与管涌的对比

流土和管涌都因渗流引起,但破坏方式截然不同——一个瞬间整体上浮,一个渐进掏空通道。混淆会导致防治策略错位。

流土
  • 破坏机理:渗透力超过土体有效自重,整体悬浮隆起
  • 发生条件:水力坡降 i ≥ 临界坡降 ic(ic≈1)
  • 适用土类:均匀细砂、粉砂等松散无黏性土
  • 工程表现:突发性整体破坏,伴随翻砂冒水
管涌
  • 破坏机理:细颗粒被渗流逐步带出,形成渗流通道
  • 发生条件:土体级配不连续,细粒能在粗粒孔隙中流动
  • 适用土类:缺乏中间粒径的不均匀级配砂砾料
  • 工程表现:渐进式发展,最终贯通形成集中渗流管道
看土的级配——均匀细砂防流土(控坡降),不均匀级配土防管涌(设反滤)。判别土类是选对防治方案的前提。
26第 26 页 · 流土与管涌的对比

流土破坏的过程

从左到右看流土全过程:渗流产生→有效应力下降→达临界→土粒悬浮→结构失稳。中间分支显示未达临界时土体仍稳定。

图解渲染中…
E渗透力:向上渗流对土粒的拖曳力,属体积力F有效应力:粒间骨架压力,决定土体抗剪强度G临界状态:有效应力归零的瞬时分界点J悬浮:粒间接触丧失,整体从固体行为变为液态
27第 27 页 · 流土破坏的过程

临界水力梯度

流土发生的极限条件与安全系数

临界水力梯度
流土发生的极限条件与安全系数
28第 28 页 · 临界水力梯度

管涌的发生条件

上一页讲了流土的临界水力梯度——它是土体整体被托起的破坏。管涌不一样:它是细颗粒被渗流逐个带走,土骨架慢慢架空。发生管涌,土本身得先具备特定结构,再配合水力条件才能启动。

粗颗粒骨架
大颗粒相互咬合搭成稳定骨架,形成连通渗流通道
细颗粒能穿过
细颗粒粒径小于骨架孔隙,能在渗流中被搬运
不连续级配
颗粒级配曲线有明显缺口,缺乏中间粒径,Cu 大
水力+出口条件
渗透力足够大,且渗流方向有畅通出口让细颗粒离开
河床卵石间被冲走的砂对应 →管涌的发生条件

卵石搭骨架=粗骨架;砂粒被冲走=细颗粒流失;洪水能量=渗透力

$$d_f < D_0 \quad \text{且} \quad i > i_{cr}$$
29第 29 页 · 管涌的发生条件

渗流控制的工程措施

三条路径解决渗流问题:上游阻、下游排、中间截。

图解渲染中…
延长渗径让水走更长的路,水力梯度 i=H/L 自然下降降低孔压把内部高水压主动导出,避免顶起上覆土层截断路径用低渗透墙直接堵断渗流通道,缩短过水断面
30第 30 页 · 渗流控制的工程措施

渗透破坏的判断

点击作答

基坑工程中为防止发生流土破坏,设计时对实际水力梯度 i 的要求是?

31第 31 页 · 渗透破坏的判断

知识框架梳理

  • 渗透系数 k 是连接微观孔隙与宏观流量的唯一枢纽参数
  • 达西定律+连续性方程→拉普拉斯方程,渗流场化为边界值问题
  • 流网是拉普拉斯方程的几何语言,正交性是它的核心约束
  • 临界水力梯度 i_cr 是渗透力与土自重的分界点,判破坏与否
  • 工程闭环:k 定量、流网定性、i_cr 定边界,缺一不可
延伸主题:非饱和土渗流各向异性与分层土的渗流耦合渗流的应力与变形分析
32第 32 页 · 知识框架梳理

工程渗流分析的全流程

工程渗流分析分三大阶段推进:①现场勘察→②渗流场求解→③稳定性评价;不满足时返回第二阶段迭代。

图解渲染中…
A3达西定律+连续性方程 → 拉普拉斯方程A7按破坏类型分流:整体上抬→流土;颗粒带出→管涌A11不满足时设计防渗/排水/反滤,返回流网重算
33第 33 页 · 工程渗流分析的全流程

深入思考

先凭直觉写下你的答案,再点开参考答案对照,不必追求完全一致。

1达西定律要求层流,若实际渗流进入湍流区,v=ki 还成立吗?为什么?

参考答案不成立。达西定律基于层流(粘滞力主导),湍流时惯性力主导,v 与 i 不再线性,k 也不再是常数。可用雷诺数 Re 判断边界。

2基坑降水设计中,如何用流网同时获得水头分布和单宽渗流量?两者如何校验?

参考答案流网等势线给出水头分布,沿任一完整流槽用 v=ki 积分求和得单宽流量。两者互相校验:流量突变常预示等势线划分有误。

3临界水力梯度 ic 假设均匀砂土,对成层土或黏性土,这公式还适用吗?边界在哪?

参考答案均匀砂土直接适用。成层土须考虑层间界面,弱透水层处临界梯度降低。黏性土有凝聚力,实际抗渗强度高于 ic 预测值。

34第 34 页 · 深入思考