(完整版)第二章土的渗透性和渗流问题要点

第二章 土的渗透性和渗流问题

第一节 概 述

土是多孔介质，其孔隙在空间互相连通。当饱和土体中两点之间存在能量差时，水就通过土体的孔隙从能量高的位置向能量低的位置流动。水在土体孔隙中流动的现象称为渗流；土具有被水等液体透过的性质称为土的渗透性。

土的渗透性是土的重要力学性质之一。在水利工程中，许多问题都与土的渗透性有关。渗透问题的研究主要包括以下几个方面：

1．渗流量问题。 例如对土坝坝身、坝基及渠道的渗漏水量的估算（图2－la 、b ），基坑开挖时的渗水量及排水量计算（图2－1C ），以及水井的供水量估算（图2－1d ）等。渗流量的大小将直接关系到这些工程的经济效益。

2．渗透变形（或称渗透破坏）问题。 流经土体的水流会对土颗粒和土体施加作用力，这一作用力称为渗透力。当渗透力过大时就会引起土颗粒或土体的移动，从而造成土工建筑物及地基产生渗透变形。渗透变形问题直接关系到建筑物的安全，它是水工建筑物和地基发生破坏的重要原因之一。由于渗透破坏而导致土石坝失事的数量占总失事工程数量的25%~30%。

3．渗流控制问题。 当渗流量和渗透变形不满足设计要求时，要采用工程措施加以控制，这一工作称为渗流控制。

渗流会造成水量损失而降低工程效益；会引起土体渗透变形，从而直接影响土工建筑物和地基的稳定与安全。因此，研究土的渗透规律、对渗流进行有效的控制和利用，是水利工程及土木工程有关领域中的一个非常重要的课题。

第二节 土的渗透性

一、土的渗透定律—达西定律

（一）渗流中的总水头与水力坡降

液体流动除了要满足连续原理外，还必须要满足液流的能量方程，即伯努里方程。在饱和土体渗透水流的研究中，常采用水头的概念来定义水体流动中的位能和动能。水头是指单位重量水体所具有的能量。按照伯努里方程，液流中一点的总水头h ，可用位置水头Z 、压力水头w u

γ和流速水头g v 22

之和表示，即 1)-(2 22

g v u

z h w ++=γ 式(2—1)中各项的物理意义均代表单位重量液体所具有的各种机械能，其量纲为长度。 对于流经土体中A 、B 二点渗流（图2－2），按照式（2－1），A 、B 两点的总水头可分别表示为：

g v u z h g

v u z h B w B B A w A A 22222

1++=++=γγ

且 h h h ?+=21

式中，A Z 、B Z — 分别为A 点和B 点相对于任意选定的基准面的高度。代表单位重量液体所具有的位能，故称z 为位置水头。

u A 、、u B — 分别为A 和B 两点的水压力即孔隙水压力（kN/m 2），代表单位重量液体所具有的压力势能。将它们除以水的容重w γ（kN/m 3）后，w A

u γ和w B u γ就分别代表A 、B 两点孔隙水压力的水柱高度，因此称w u γ为压力水头。

A v 、

B v — 分别为A 点和B 点处的渗流流速（m/s ），g 为重力加速度（m/s 2）。g v 22

即代表单位重量液体所具有的动能，故称g v 22为流速水头。

1h 、2h — 分别代表A 点和B 点单位重量液体所具有的总机械能，故称之为总水头。 h ?— A 、B 二点间的总水头差，代表单位重量液体从A 点向B 点流动时，为克服阻力而损失的能量。 另外，我们常将位置水头与压力水头之和w u

z γ+称为测压管水头。如果将二根测压管

分别安装在点A 和点B 处时，测压管中的水面将会分别上升至A Z +

w A u γ和B Z +w B u γ的标高

处。所以，测压管水头代表的是单位重量液体所具有的总势能。 饱和土体中两点间是否会出现渗流是由总水头差h ?（=1h -2h ）决定的。只有当两点间的总水头差h ?>0时，才会发生水从总水头高的点向总水头低的点流动。由于土中渗流阻力大，流速v 在一般情况下都很小，因此流速水头也很小，为简便起见可以忽略。这样，渗流中任一点的总水头就可用测压管水头来代替，式（2—1）可简化为

5)-(2 w u

z h γ+=

将土体中A 、B 两点的测压管水头连接起来（图2－2），就得到测压管水头线（又称水力坡降线）。由于渗流过程中存在能量损失，测压管水头线沿渗流方向下降。A 、B 两点间的水头损失，可用无量纲的形式来表示，即

6)-(2 L

h i ?= 这里，i 称为水力坡降，L 为A 、B 两点间渗流路径的长度。水力坡降i 代表单位渗流长度上的水头损失。

（二）渗透试验与达西定律

达西利用图2－5所示的试验装置对均匀砂进行了渗流试验，得出了层流条件下，土中水的渗流速度与能量（水头）损失之间的渗流规律，即达西定律。

达西试验装置的主要部分是一个上端开口的直立圆筒，下部放碎石，碎石上放一块多孔滤板，滤板上面放置颗粒均匀的土样，其断面积为A ，长度为L 。筒的侧壁装有两支测压管，分别设置在土样两端的两个过水断面处。水由上端进水管注入圆筒，并以溢水管保持简内为恒定水位。透过土样的水从装有控制阀门的弯管流入容器中。

当筒的上部水面保持恒定以后，通过砂土的渗流是恒定流，测压管中的水面将恒定不变。取图2—5中的0—0界面为基准面，1h 、2h 分别为上下断面处的测压管水头；h ?即为渗流流经L 长度砂样后的水头损失。

达西根据对不同类型及长度的土样所进行的试验发现，渗出水流量Q [L 3 T -1]与圆筒断面积A [L 2]和水力坡降i [L L -1]成正比，且与土的透水性质有关，即

8)-(2 kAi Q =

或 9)-(2 ki A

Q v == 式（2一8）或（2一9）称为达西定律。

式中，v —断面平均渗透速度[L T -1]，单位mm/s 或m/day

k —土的渗透系数，它反映了土的透水性能的大小。渗透系数相当于水力坡降i =1时的渗透速度[L T -1]，故其量纲与流速相同，mm/s 或m/day 。

达西定律说明，在层流状态的渗流中，渗透速度v 与水力坡降i 的一次方成正比，并与土的性质有关。

渗透流速v 并不是土孔隙中水的实际平均流速。在公式推导中采用的是土样的整个断面积，其中包括了土粒骨架所占的部分面积在内。土粒本身是不能透水的，实际的过水面积A v 应小于A ，从而实际平均流速v s 应大于v 。一般称v 为达西流速。v 与v s 的关系可通过水流连续原理建立。

按照水流连续原理，

Q=vA=v s A v （2－10）

若均质砂土的孔隙率为n （砂土孔隙率为0.28~0.35；粘性土0.6~0.7），则A v =nA ，

11)-(2 /n v nA

vA v s ==

（三）达西定律的适用范围

达西定律是描述层流状态下渗透流速与水头损失之间关系的规律，即渗流速度v 与水力坡降i 成线性关系只适用于层流范围。在水利工程中，绝大多数渗流，无论是发生于砂土中或一般的粘性土中，均属于层流范围，故达西定律均可适用。

但须注意的是，在纯砾等粒径很粗的土中例如堆石体中的渗流，当水力坡降较大时，流态已不再是层流而是紊流，达西定律不再适用，此时渗流速度v 与水力坡降i 之间的关系不再保持直线而变为曲线关系（图2－6a ）。层流进入紊流的界限就为达西定律适用的上限。一般可用临界流速cr v =0.3~0.5cm/s 来划分这一界限。当v ＞v cr 后达西定律可修改为：

12)-(2 1 m ki v m <=

在粘性很强的致密土体中，渗透特征也偏离达西定律。此时v ～i 关系（图2－6b ）也

呈曲线规律，且不通过原点。使用时，可将曲线简化为如图虚线所示的直线关系。截距i 0称为起始坡降。这时，达西定律可修改为：

13)-(2 )(0i i k v -=

式（2－13）说明，当坡降很小即i< i 0时，没有渗流发生。因为密实粘土颗粒的外围具有较厚的结合水膜，它占据了土体内部的过水通道（图2－7），因此只有在较大的水力坡降作用下，挤开结合水膜的堵塞后才能发生渗流。起始水力坡降i 0是用以克服结合水膜阻力所消耗的能量。i= i 0就是达西定律适用的下限。

二、渗透系数的测定和影响因素

渗透系数k 是一个代表土的渗透性强弱的定量指标。不同种类的土，k 值差别很大。

（一）渗透系数的测定方法

渗透系数的测定方法主要分室内测定和野外现场测定两大类。

1．实验室测定法

目前在实验室中测定渗透系数k 的仪器种类和试验方法很多，但从试验原理上大体可分为常水头法和变水头法两种。

常水头试验法就是在整个试验过程中保持水头为一常数，从而水头差也为常数。试验时（图2－9a ），在透明塑料筒中装填截面为A 、长度为L 的饱和试样，打开阀门，使水自上而下流经试样，并自出水口处排出。待水头差?h 和渗出流量Q 稳定后，量测经过一定时间t 内流经试样的水量V ，则vAt Qt V ==

根据达西定律ki v =，则At L

h k

V ?= 从而得出 14)-(2 ht A VL k ?= 常水头试验适用于测定透水性大的砂性土的渗透系数。

变水头试验法就是试验过程中水头差一直在随时间而变化，其装置示意图见图2－9b 。水流从一根直立的带有刻度的玻璃管和U 形管自下而上流经土样。试验时，将玻璃管充水至需要的高度后，开动秒表，测记起始水头差?h 1，经过时间t 后，再测记终了水头差?h 2。通过建立瞬时达西定律，即可推出渗透系数k 的表达式。

设试验过程中任意时刻t 时，作用于试样两端的水头差为?h ；经过dt 时段后，管中水位下降dh ，则dt 时间内流入试样的水量为dh a dV e ?-=，式中a 为玻璃管断面积；右端的负号表示水量随?h 的减少而增加。

根据达西定律，dt 时间内流出试样的渗流量为：

dt A L

h k dt A i k dV o ????=???= 式中，A —试样断面积；L —试样长度。

根据水流连续原理，应有 o e dV dV =，即

h dh kA aL dt dt A L h k dh a ??-=????

=?- 等式两边各自积分?????-=210h h t

h

dh kA aL dt 得： 2

1ln h h kA aL t ??= 从而得到土的渗透系数：

15)-(2 ln 2

1h h At aL k ??= 改用常用对数表示，则上式可写为

16)-(2 lg 3.22

1h h At aL k ??= 通过选定几组不同的1h ?、2h ?值，分别测出它们所需的时间t ，利用式（2－15）或式（2-16）计算它们的渗透系数k ，然后取平均值，作为该土样的渗透系数。

实验室内测定渗透系数k 的优点是设备简单，费用较省。但是，由于土的渗透性与土的结构有很大的关系，地层中水平方向和垂直方向的渗透性往往不一样；再加之取样时的扰动，不易取得具有代表性的原状土样，特别是砂土。因此，室内试验测出的k 值常常不能够很好地反映现场中土的实际渗透性质。为了量测地基土层的实际渗透系数，可直接在现场进行k 值的原位测定。

2．现场测定法

在现场进行渗透系数k 值的测定时，常用现场井孔抽水试验或井孔注水试验的方法。 现场井孔抽水试验示意图见图2—10。在现场打一口试验井，贯穿要测定k 值的砂土层，并在距井中心不同距离处设置两个观测孔。然后自井中以稳定流量连续抽水。抽水时，井周围地下水位形成一个以井孔为中心的降落漏斗。测压管水头差形成的水力坡降使水流流向井内。假定水流是水平流向时，则流向水井的渗流过水断面应是一系列的同心圆柱面。待出水量和井中的动水位稳定后，测定抽水量Q 以及观测孔内的水位高度h 1和h 2。如果设观测孔距井轴线的距离分别为r 1和r 2，通过达西定律即可求出土层的平均k 值。

围绕井轴取一过水断面，该断面距井中心距离为r ，水面高度为h ，则过水断面积A 为rh A π2=。

假设该过水断面上各处水力坡降为常数，且等于地下水位线在该处的坡度时，则dr

dh i = 根据达西定律，单位时间自井内抽出的水量为 dr dh k

rh Aki Q ?==π2 整理得 khdh r

dr Q π2=

等式两边进行积分 ??=212

12h h r r hdh k r

dr Q π 得 )(ln 21221

2h h k r r Q -=π 所以 17)-(2 )()/ln(212212h h r r Q k -=

π

用常用对数表示时则为 18)-(2 )()

/lg(3.2212212h h r r Q k -=π

现场测定法可获得较为可靠的渗透系数值，但试验所需费用较多。

注水试验的原理与抽水试验类似，我们不再进一步讨论。

（二）影响渗透系数的因素

渗透系数k 综合反映了水在土体中运动的难易程度。土的性质和水的性质均对其有一定的影响。

1．土的性质对k 值的影响

土的粒径大小与级配、孔隙比、矿物成分、结构以及饱和度等性质均对k 值有很大的影响，其中粒径大小和孔隙比对k 的影响最大。粒径越大，k 值越大；孔隙比越大，k 值越大；凝聚结构比分散结构具有更大的透水性；土层水平方向的透水性大于垂直方向的透水性。渗透系数随饱和度的增加而增大。

2．渗透水的性质对k 值的影响

水的性质对渗透系数k 值的影响主要是由于粘滞度不同所引起。温度高时，水的粘滞性降低，k 值变大；反之k 值变小 。表2—1列出了各类土渗透系数的变化范围，可供参考。

三、层状地基的等效渗透系数

大多数天然沉积土层是由渗透系数不同的几层土所组成。在计算渗流量时，为简单起见，常常把几个土层等效为厚度等于各土层之和、渗透系数为等效渗透系数的单一土层。但等效渗透系数的大小与水流的方向有关。

（－）水平渗流情况

设一座建造在多层透水地基上的水闸（图2－13a ），已知地基内各层土的渗透系数分别为k 1、k 2、k 3、…，厚度相应为H 1、H 2、H 3…，总土层厚度即等效土层厚度为H 。当渗透水流自断面1一1流至断面2－2时，渗透距离为L ，水头损失为?h 。这种平行于各层面的水平渗流具有如下的特点：

（1）各层土中的水力坡降i （L

h ?=）与等效土层的平均水力坡降i 相同。 （2）垂直X －Z 面取单位宽度，通过等效土层H 的总的单宽渗透流量等于各土层渗流流量之和，即

21)-(2 1321∑==+++=n

j jx x x x x q q q q q K

将达西定律代入式（2－21）可得

∑∑====n j n

j j j j j x H k i iH k iH k 11

消去i 后，即可得出沿水平方向的等效渗透系数k x 。

∑==n j j j x H k

H k 1 22)-(2 1

即平行于层面的等效渗透系数k x 值是各土层渗透系数按厚度的加权平均值。

(二)竖直渗流情况

对于渗流垂直于土层的情况图（2-13b ），设承压水流流经土层H 厚度的总水头损失为?h ，流经每一层土的水头损失为1h ?、2h ?、3h ?……。这种垂直于各层面的渗流具有如下特点：

（1）根据水流连续原理，流经各土层的流速与流经等效土层的流速相同，即

23)-(2 321v v v v ====K

（2）流经等效土层H 的总水头损失面?h 等于各层土的水头损失之和，即

∑=?==?=?=?=?n

j j h h h h h 132124)-(2 K

将达西定律代入式（2－23），则

v H h k H h k H h k H h k j

j j =?==?=?=?K 333222111 从而可解出

25)-(2 222111j

j j k vH h k vH h k vH h =?=?=?Λ、、 设竖直等效渗透系数为k z ，对等效土层有H h k v z

?=，从而可得 26)-(2 z

k vH h =? 将式（2－26）和（2－25）代入式（2－24）得

∑==n j j j z k vH k vH 1 消去v ，即可得出垂直于土层方向的等效渗透系数k z

27)-(2 1∑==

n j j

j z k H H k

上式表明，垂直于层面的等效渗透系数k z 主要由渗透系数小的土层所控制。

第三节 二维渗流与流网

当遇到一维渗流问题时，可直接利用达西定律进行渗流计算。但工程上遇到的渗流问题常常属于边界条件较为复杂的二维或三维渗流问题，例如闸坝下透水地基的渗流、土坝坝身的渗流等（图2－14），此时达西定律需要用微分形式来表达。为了求解和评价渗流在地基或坝体中是否会造成有害影响，需要知道整个渗流场中各处的测压管水头、渗透坡降和渗流速度。当闸坝很长且断面轮廓一致时，可按二维平面渗流问题处理。

一、平面渗流的基本方程

设混凝土坝上下游水位差h ?保持恒定（图2－14a ），坝基下的渗流即为稳定渗流。这时，渗流场中的测压管水头h 以及流速v 等渗流要素仅是位置的函数而与时间无关，即),(),(z x f v z x f h ==，。

从稳定渗流场中任意点A 取一微元土体，其面积为dx ·dz 、厚度为dy=1，在x 和z 方向的流速为v x 和v z （图2－15）。

单位时间内流入这个微元体的水量为dq e ，则11?+?=dx v dz v dq z x e

单位时间内流出这个微元体的水量为dq o ，则

11???? ????++???? ?

???+=dx dz z v v dz dx x v v dq z z x x o 假定水体不可压缩，则根据水流连续原理，单位时间内流入和流出微元体的水量应相等，即dq e =dq o

从而得出 28)-(2 0=??+??z

v x v z x 式（2－28）即为二维渗流连续方程。

根据达西定律，对于各向异性土，

30)-(2 29)-(2 z h k i k v x h k i k v z z z z x

x x x ??==??== 式中，k x 、k z —分别为x 和z 方向的渗透系数；h —测压管水头。

将式（2－29）和（2－30）代入式（2－28）可得出：

31)-(2 02222=??+??z

h k x h k z x 对于各向同性的均质土，k x =k z ，则式（2－31）可表达为：

32)-(2 02222=??+??z

h x h 式（2－32）即为著名的拉普拉斯方程，它是平面稳定渗流的基本方程式。该方程描述

了渗流场内部测压管水头h 的分布。通过求解一定边界条件下的拉普拉斯方程，即可求得该条件下的渗流场。

拉普拉斯方程式的求解方法大致可分为四种类型：

1．数学解析法

根据具体边界条件，以解析法求式（2－32）的解。一般通解较易得到。满足拉普拉斯微分方程的解是两个共轭调和函数，即势函数Φ(x,z)和流函数ψ(x,z)，该函数描绘出两族相互正交的曲线即等势线和流线。该法的缺点是当边界条件复杂时定解较难求得。

2．数值解法

数值解法是一种近似方法，常用的数值解法有有限差分法和有限元法。随着计算机的发展，数值解法的精度愈来愈高，因而数值解法的应用也愈来愈广。

3．实验法

实验法即采用一定比例的模型来摸拟真实的渗流场，用实验手段测定渗流场中的渗流要素。例如，电比拟法就是利用渗流与电流现象存在着的比拟关系，来实测渗流等势线族的一种实验方法。其它的实验法还有电网络法等。

4．图解法

图解法即用绘制流网的方法求解拉普拉斯方程的近似解。该法具有简便、迅速的优点，并能用于建筑物边界轮廓较复杂的情况。只要满足绘制流网的基本要求，精度就可以得到 保证，因而该法在工程上得到广泛应用。

二、流网的绘制及应用

在透水地基上混凝土坝下的渗流流网图见图2-16，图中标有号码①②③……的线表示流线。在稳定渗流场中，流线表示水质点的运动路线。图中标有号码1、2、3、…的线代表等势线，等势线是渗流场中势能或测压管水头的等值线。如果在同一等势线上的不同点处，例如图中的a 、b 两点处，安放测压管，则管中水位将升至相同的高度。由流线和等势线所组成的曲线正交网格称为流网。

（一）绘制流网的基本要求

绘制流网时必满足下列几个条件：

（1）流线与等势线必须正交。

（2）流线与等势线构成的各个网格的长宽比应为常数，即?l /?s =C 。当取s l ?=?时，网格应呈曲线正方形，这是绘制流网时最方便和最常见的一种流网图形。

（3）必须满足流场的边界条件，以保证解的唯一性。

（二）流网的绘制方法

(1) 首先根据渗流场的边界条件，确定边界流线和边界等势线。在图2-16中，渗流是有压渗流，因而坝基轮廓线A —B —C —D 是第一条流线；不透水层面0－0也是一条边界流线。上下游透水地基表面1一A 和D —11则是两条边界等势线。

（2）根据绘制流网的另外两个要求，即流线与等势线必须正交以及流线与等势线构成的各个网格的长宽比应为常数，初步绘制流网。按边界趋势先大致画出几条流线如②、③、④，彼此不能相交，且每条流线都要和上下游透水地基表面（等势线）正交。然后再自中央向两边画等势线，图2－16中先绘中线6，再绘5和7，如是向两侧推进。每根等势线要与流线正交，并弯曲成曲线正方形。

（3）对初绘的流网进行反复修改，直至大部分网格满足曲线正方形为止。因为边界形状不规则，在边界突变处很难画成正方形，而可能是三角形或五边形。这主要是由于流网图

中流线和等势线的根数有限所造成的。只要网格的平均长度和宽度大致相等，就不会影响整个流网的精度。

（三）流网的应用

流网绘出后，即可求得渗流场中各点的测压管水头、水力坡降、渗透流速和渗流量。现仍以图2－16所示的流网为例。

1．测压管水头

根据流网特征可知，任意两相邻等势线间的势能差相等，即水头损失相等，从而可算出相邻两条等势线之间的水头损失?h ，即

33)-(2 1)-n (N 1

=-?=?=?n H N H h 式中，?H —上、下游水位差，也就是水从上游渗到下游的总水头损失；

N —等势线的间隔数；

n —等势线的个数。

在图2-16中，n =11，N=10，?H=5.0m ，故每一个等势线间隔所消耗的水头?h=5/10=0.5m 。有了?h 就可求出任意点的测压管水头。例如求a 点的测压管水头h a ：以0－0为基准面，h a =h ua ＋z a ；z a 为a 点的位置高度，为已知值。由于a 点位于第2条等势线上，所以测压管水位应比上游水位降低一个?h ，故其测压管水位应在上游地表面以上的6.0-0.5=5.5 m 处。压力水头h ua 的高度可自图中按比例直接量出。

2．孔隙水压力

渗流场中各点的孔隙水压力等于该点以上测压管中的水柱高度h u 乘以水的容重γw ，故a 点的孔隙水压力为

34)-(2 w ua a h u γ?=

应当注意，图中所示a 、b 两点位于同一根等势线上，其测压管水头虽然相同，即 h a =h b ，但压力水头h ua ≠h ub ，所以其孔隙水压力不相等，即u a ≠u b 。

3．水力坡降 流网中任意网格的平均水力坡降l

h i ??=，?l 为该网格处流线的平均长度，可自图中量出。根据水力坡降的定义可知，流网中网格越密处，其水力坡降越大。因此在图2－16中，下游坝趾水流渗出地面处（图中CD 段）的水力坡降最大。该处的水力坡降称为逸出坡降，它常是地基渗透稳定的控制坡降。

4．渗透流速

各点的水力坡降已知后，渗透流速的大小可根据达西定律求出，即v =ki ，其方向为流线的切线方向。

5．渗透流量 因为s l

h k

s ki A v q ???=???=?=?1 当取?l =?s 时， 35)-(2 h k q ??=?

由于?h 是常数，故?q 也是常数。因此，流网中任意两相邻流线间的单宽流量?q 是相等的。

通过坝下渗流区的总单宽流量：

36)-(2 h Mk q M q q ??=??=?=∑

式中M 为流网中的流槽数，数值上等于流线数减1。在图2－16中M=4。

通过坝底的总渗透流量

37)-(2 qL Q =

式中L 为坝基的长度。

第四节 渗透力和渗透变形

一、渗透力和临界水力坡降

（一）渗透力

设定水头渗透破坏试验装置的土样长度为L （图2－21）、面积为A ＝1，士样两端各安装一测压管，其测压管水头相对0－0基准面分别为h 1与h 2。当h 1=h 2时，土中水处于静止状态，无渗流发生。若将左侧的贮水器向上提升，使h 1＞h 2，则由于存在水头差，土中将产生向上的渗流。水头差?h 是渗流穿过L 长土体时所损失的能量。能量损失说明土粒对水流施加了阻力，因此渗透水流必然对每个土颗粒施加推动和拖拉作用力（图2－22）。定义每单位土体内土颗粒所受的渗流作用力为渗透力，用j 表示。为了求出渗透力的大小，对图2－21中的土体进行受力分析。

1．土体的受力分析

首先，取土和水整体作为隔离体进行受力分析。则作用在土柱上的力（图 2－23a ）有：

（1）土和水总重量W=γsat L ；

（2）土柱两端的边界水压力：w w h γ和1h w γ；

（3）土柱下部滤网的支承反力R 。

此时，土粒与水之间的作用力为内力，在土柱的受力分析中不出现。因此，为了求解渗透力，须把土骨架和水分开来取隔离体（图2－23 b 、c ）。显然，这两种隔离体取法的总效果是一样的，即图2－23中(a)= (b)+(c)。

对于土骨架隔离体（图2－23b ），由于土骨架浸没于水中，土颗粒受浮力作用，其值等于排开同体积的水重，因此在计算土粒重量时应用浮容重γ’。另外，由于将土骨架与水体分开，则土颗粒上受到的水流作用力—渗透力即成为外力。因此，作用在土柱内土骨架上的作用力有：

（1）土粒有效重量W’=γ’L ；

（2）总渗透力J=jL ，方向竖直向上；

（3）下部支承反力R 。

对于土柱中的孔隙水隔离体（图2－23C ），作用在其上的力有：

（1）孔隙水重量与土粒浮力的反力之和。土粒浮力的反力应等于与土颗粒同体积的水重。因此有，

w w w s w v w L V V V W γγγγ==+=

即W w 就为L 长度的水柱重量。

（2）水柱上下两端面的边界水压力，w w h γ和1h w γ；

（3）土柱内土粒对水流的阻力，其大小应和渗透力相等，方向相反。设单位土体内土粒给水流的阻力为j ’，则总阻力J’=j’L ，方向竖直向下。

2．渗透力的计算

现考虑水体隔离体（图2－23c ）的平衡条件，可得1'h J W h w w w w γγ=++，则

i L h

L h h L L h h L L h h j h L j L h w w w w w w w w w w w γγγγγγγγ=?=-=+-=--==++)

()]

([)

( 2111'1

'所以

故渗透力 49)-(2 'i j j w γ==

从式（2－49）可知，渗透力是一种体积力，量纲与γw 相同。渗透力的大小和水力坡降成正比，其方向与渗流方向一致。

（二）临界水力坡降

若将图2－21中左端的贮水器不断提高，则?h 逐渐增大，从而作用在土体中的渗透力也逐渐增大。当?h 增大到某一数值时，向上的渗透力克服了向下的重力，土体就要发生浮起或受到破坏，这一现象就称为流土。

分析图2－23b 中土骨架隔离体的平衡状态。当发生流土时，土柱压在滤网上的压力

R=0，故0'==-R J W ，即0'=-jL L γ，所以j ='

γ。因为cr w i j γ=，所以 50)-(2 '

w

cr i γγ= i cr 称为临界水力坡降，它是土体开始发生流土破坏时的水力坡降。已知土的浮容重γ’为：

e

G w s +-=1)1('γγ 将其代入式（2－50）后可得

51)-(2 11e

G i s cr +-= 式中Gs 、e 分别为土粒比重及土的孔隙比。由此可知，流土的临界水力坡降取决于土的物理性质。

二、土的渗透变形（或称渗透破坏）

土工建筑物及地基由于渗流作用而出现的变形或破坏称为渗透变形或渗透破坏。渗透变形是水工建筑物发生破坏的重要原因之一。

（－）渗透变形的类型

单一土层土的渗透变形主要有管涌和流土两种基本型式。

1．流土

在向上的渗透水流作用下，表层局部范围内的土体或颗粒群同时发生悬浮、移动的现象称为流土。不论任何类型的土，只要水力坡降达到一定的大小，都会发生流士破坏。

流土常发生在堤坝下游渗流逸出处无保护的情况下。例如，建筑在双层地基上的堤坝（图2－25），地基表层为渗透系数较小的粘性土层，且较薄；下层为渗透性较大的无粘性土层，且k 1<

若地基为比较均匀的砂层，当水位差较大、渗透途径不够长时，下游渗流逸出处也会出现i ＞i cr 的情况。这时地表将会出现小泉眼，冒气泡，土颗粒群会向上鼓起，发生浮动、跳跃，这一现象称为砂沸。砂沸也是流土破坏的一种形式。

2．管涌

在渗透水流作用下，土中的细颗粒在粗颗粒形成的孔隙中移动，以至流失；随着土的孔隙不断扩大，渗透流速不断增加，较粗的颗粒也相继被水流逐渐带走，最终导致土体内形成贯通的渗流管道，造成土体塌陷，这种现象称为管涌（图2－26）。管涌破坏一般是渐进性的破坏。管涌发生在一定级配的无粘性土中，发生的部位可以在渗流逸出处，也可以在土体内部。

（二）渗透破坏类型的判别

土的渗透变形的发生和发展过程有其内因和外因。内因是土的颗粒组成和结构，即常说的几何条件； 外因是水力条件，即作用于土体渗透力的大小。

1．流土可能性的判别

在自下而上的渗流逸出处，任何土，不论粘性土还是无粘性土，只要满足渗透坡降大于临界水力坡降这一水力条件，均要发生流土。因此，只要知道水力坡降和临界水力坡降，即可判别流土的可能性。

若i ＜i cr ，土体处于稳定状态；

i ＞i cr ，土体发生流土破坏；

i=i cr ，土体处于临界状态。

由于流土将造成地基坡坏、建筑物倒塌等灾难性事故，工程上是绝对不允许发生的，故设计时要保证有一定的安全系数，把逸出坡降限制在允许坡降［i ］以内，即

52)-(2 ][s

cr F i i i =≤ 式中F s 为流土安全系数。按照我国《碾压式土石坝设计规范》规定，取F s =1.5~2.0。

2．管涌可能性的判别

土是否发生管涌，首先决定于土的性质。一般粘性土只会发生流土而不会发生管涌，故属于非管涌土；无粘性土中产生管涌必须具备下列两个条件：

（1）几何条件

土中粗颗粒所构成的孔隙直径必须大于细颗粒的直径，才可能让细颗粒在其中移动，这是管涌产生的必要条件。

对于不均匀系数Cu ＜10的较均匀土，颗粒粗细相差不多，粗颗粒形成的孔隙直径不比细颗粒大，因此细颗粒不能在孔隙中移动，也就不可能发生管涌。

对于Cu＞10的不均匀砂砾石土，这种土既可能发生管涌也可能发生流土，主要取决于土的级配情况和细颗粒含量。

对于缺乏中间粒径、级配不连续的土，其渗透变形的型式主要决定于细颗粒含量。当细颗粒含量在25％以下时，细料填不满粗料所形成的孔隙，渗透变形基本上属管涌型；当细料含量在35％以上时，细料足以填满粗料所形成的孔隙，粗细料形成整体，抗渗能力增强，渗透变形是流土型；当细料含量在25~35％之间时，则是过渡型。

对于级配连续的不均匀土，用土的孔隙平均直径D0与最细部分的颗粒粒径d s相比较，以判别土的渗透变形的类型。土的孔隙平均直径D0可以下述经验公式计算：

D0=0.25d20 （2－53）

式中d20为小于该粒径的土质量占总质量的20％。当土中有5％以上的细颗粒小于土的孔隙平均直径时，即D0＞d5 时，破坏形式为管涌；而如果土中小于D0的细粒含量＜3％，即D0＜d3时，可能流失的土颗粒很少，不会发生管涌，则呈流土破坏。

因此，对于无粘性土是否发生管涌的几何条件可用下列准则判别：

A．C u 10的比较均匀的土非管涌土。

B．C u>10的不均匀土

（a）级配不连续的土

细颗粒含量>35％非管涌土；

细料含量<25％管涌土；

细料含量=25~35％过渡型土。

（b）级配连续的土

D0

D0>d5 管涌土；

D0= d3~d5过渡型土。

（2）水力条件

渗透力能够带动细颗粒在孔隙间滚动或移动是发生管涌的水力条件。管涌土的破坏坡降与允许坡降的范围值见表2－3。

（二）渗透变形的防治措施

防治流土的关键在于控制逸出处的水力坡降。为了保证实际的逸出坡降不超过允许坡降，水利工程上常采取下列工程措施。

（1）上游做垂直防渗帷幕，例如混凝土防渗墙、板桩或灌浆帷幕等。其作用主要是延长渗流途径，降低下游逸出处的坡降。

（2）上游做水平防渗铺盖，以延长渗流途径、降低下游的逸出坡降。

（3）下游挖减压沟或打减压井，贯穿渗透性小的粘性土层，以降低作用在粘性土层底面的渗透压力。

（4）下游加透水盖重，以防止土体被渗透力所悬浮。

防止管涌一般可从两方面采取措施：

（1）改变水力条件，降低土层内部和渗流逸出处的渗透坡降。如上游做防渗铺盖或打板桩等。

（2）改变几何条件。在渗流逸出部位铺设层间关系满足要求的反滤层是防止管涌破坏的有效措施。反滤层用以保护基土不让细颗粒带出，同时也具有较大的透水性，使渗流畅通。