1、水文地质学第4章 地下水运动的基本规律改.ppt4.1 4.1 基本概念基本概念4.1.1 渗透渗透(渗流渗流) 地下水受重力作用在岩石空隙中的运动地下水受重力作用在岩石空隙中的运动称为渗透称为渗透(渗流渗流) 。 地下水在岩石空隙中运动非常复杂（如地下水在岩石空隙中运动非常复杂（如下图），为了研究地下水的运动，将实际地下图），为了研究地下水的运动，将实际地下水流进行概化。下水流进行概化。 概化条件：概化条件：不考虑骨架，骨架和空隙不考虑骨架，骨架和空隙全部被水充满；不考虑地下水实际运动途全部被水充满；不考虑地下水实际运动途径的迂回曲折，只考虑运动的总体方向。径的迂回曲折，只考虑运动的总体方向
2、。4.1.2 渗透的运动要素渗透的运动要素1、实际流速与渗透流速：、实际流速与渗透流速：实际流速：实际流速：实际水流通过实际水流通过单位空隙过水断面单位空隙过水断面的流的流量称为实际流速量称为实际流速。 u =Q/w渗透流速：渗透流速：假想水流通过假想水流通过单位过水断面单位过水断面的流量称的流量称为渗透流速为渗透流速。 V=Q/ w松散岩石孔隙中的水示意图松散岩石孔隙中的水示意图渗流必须满足下列条件：渗流必须满足下列条件： 1）通过任一断面的渗透流量等于通过该断）通过任一断面的渗透流量等于通过该断面的实际流量；面的实际流量； 2）作用于任一面积的渗透压力或水头等于）作用于任一面积的渗透压力或
3、水头等于作用于该面积的实际压力或水头；作用于该面积的实际压力或水头； 3）渗透通过任一体积所受的阻力等于实际）渗透通过任一体积所受的阻力等于实际水流通过该体积所受到的阻力。水流通过该体积所受到的阻力。实际流速实际流速(u)与渗透流速与渗透流速(V)的关系：的关系：V = ne u u V V = （w / w ） une-有效孔隙度有效孔隙度 Q= u w Q=V w 2、水头与水力梯度：、水头与水力梯度： 水头：水头：渗流场内任一点的测压水头（渗流场内任一点的测压水头（Hn）是该）是该点测压高度（点测压高度（hn）与此点到基准面距离（）与此点到基准面距离（Z）之和）之和。Hn=Z+hn=Z+
4、 P /；P-A点静水压强，点静水压强，-水的容水的容重。重。 水力梯度：水力梯度：渗流通过渗流通过A点单位渗流途径长度点单位渗流途径长度上的水头损失称为水力坡度上的水头损失称为水力坡度。 I = dH/dL 或或 I =（H1H2）/ L水水力力梯梯度度示示意意图图3、流线、迹线、等水头线和流网：、流线、迹线、等水头线和流网：流线：流线：同一时间内不同液流质点的连线同一时间内不同液流质点的连线，这这根连线上的各液流质点速度矢量都与这根连根连线上的各液流质点速度矢量都与这根连线相切线相切。迹线：迹线：某一液流质点在某一时间段内的运动某一液流质点在某一时间段内的运动轨迹轨迹。等水头线：等水头线：
5、水头值相等点的连线水头值相等点的连线。 a) 潜水等水位线：潜水等水位线： 潜水面上任一点的高程称为该点的潜水位，潜水面上任一点的高程称为该点的潜水位，潜水位相等的各点连线称为潜水等水位线潜水位相等的各点连线称为潜水等水位线。 从潜水等水位线图可获得如下信息：从潜水等水位线图可获得如下信息： 确定潜水流向、确定水力坡度、确定潜水确定潜水流向、确定水力坡度、确定潜水与地表水之间的关系、确定潜水埋深，判断与地表水之间的关系、确定潜水埋深，判断泉水、沼泽出露点。泉水、沼泽出露点。b）承压水等水压线：）承压水等水压线： 承压含水层测压承压含水层测压水头相等的各点的连线为承水头相等的各点的连线为承压水等
6、水压线压水等水压线。 从承压水等水压线图可获得如下信息：从承压水等水压线图可获得如下信息： 确定承压水流向、确定承压水流向、 确定承压水水力坡度。确定承压水水力坡度。 等水位（压）线密，等水位（压）线密，水力坡度大；水力坡度大； 等水位（压）线疏，等水位（压）线疏，水力坡度小。水力坡度小。 流网：流网： 在渗流场中，流线与等水头线组成的网格称在渗流场中，流线与等水头线组成的网格称为为流网流网。 在各向同性介质中为在各向同性介质中为正交网正交网； 在各向异性介质中为在各向异性介质中为斜交网斜交网。流网的类型：流网的类型：a) 均质各向同性介质中的流网均质各向同性介质中的流网 均质：均质：在同一含
7、水层中，各处的渗透系数在同一含水层中，各处的渗透系数相等相等。 各向同性：各向同性：含水层中任一点的渗透系数在含水层中任一点的渗透系数在各个方向无变化各个方向无变化。1）在河渠附近的流网：）在河渠附近的流网： 地表水体的断面看作等地表水体的断面看作等水位面，地表水的湿周是水位面，地表水的湿周是一条等水位线。一条等水位线。2）承压含水层：）承压含水层：隔水边隔水边界无水量通过，流线平行界无水量通过，流线平行隔水边界。隔水边界。3）无入渗补给和蒸发排泄）无入渗补给和蒸发排泄时：时：潜水面是一条流线。潜水面是一条流线。4）有入渗补给时：）有入渗补给时：潜水潜水面既不是流线，也不是等面既不是流线，也不
8、是等水位线。水位线。均质各向同性介质中均质各向同性介质中流网示意图流网示意图b) 层状非均质介质中的流网层状非均质介质中的流网: 层状非均质：层状非均质：介质场内各岩层内部渗透性介质场内各岩层内部渗透性均为各向同性，不同层介质的渗透性不同。均为各向同性，不同层介质的渗透性不同。 两层平行等厚渗透两层平行等厚渗透系数分别为系数分别为K1、K2的的岩层，岩层，K2=3K1：等水等水位线间隔一致，流线位线间隔一致，流线密度密度K2为为K1的三倍的三倍 。 流线通过不同渗透流线通过不同渗透性的两套地层性的两套地层K2=3K1，等水位线密度等水位线密度K1为为K2的三倍，流线相等。的三倍，流线相等。非均
9、质介质流网图非均质介质流网图 含水层中有强渗透性透镜体时：含水层中有强渗透性透镜体时：流线向其汇聚。流线向其汇聚。 含水层中有弱渗透性透镜体时：含水层中有弱渗透性透镜体时：流线将绕流。流线将绕流。非均质介质流网图非均质介质流网图河间地块流网图河间地块流网图4.1.3 渗流的分类渗流的分类1、有压流与无压流：、有压流与无压流：有压流：有压流：渗流场中任一点处的压强都不为大气渗流场中任一点处的压强都不为大气压强（一般大于大气压强）的渗流为有压流压强（一般大于大气压强）的渗流为有压流。无压流：无压流：具有自由表面且表面压强为大气压强具有自由表面且表面压强为大气压强的渗流为无压流的渗流为无压流。2、稳
10、定流与非稳定流：、稳定流与非稳定流： 渗流场中任一点处的运动要素（水位、渗流场中任一点处的运动要素（水位、流速、流向等）不随时间变化的渗流称为流速、流向等）不随时间变化的渗流称为稳稳定流。定流。 例如水头：例如水头：H=f(x,y,z)； 否则为非稳定流。否则为非稳定流。例如水头：例如水头：H=f(x,y,z,t)。3、层流与紊流：、层流与紊流：地下水在岩石空隙中渗流时，地下水在岩石空隙中渗流时，水的质点有秩序、互不混杂有规则的运动称为水的质点有秩序、互不混杂有规则的运动称为层流层流；否则为；否则为紊流紊流。 4、一维流、二维流、三维流、一维流、二维流、三维流：一维流：一维流：在渗流场中，速度
11、向量与任一坐标在渗流场中，速度向量与任一坐标轴相一致的渗流称为一维流轴相一致的渗流称为一维流。二维流：二维流：渗流场中，速度向量与某一坐标平渗流场中，速度向量与某一坐标平面平行的渗流称为二维流面平行的渗流称为二维流。三维流：三维流：渗流场中，速度向量不与某一坐标渗流场中，速度向量不与某一坐标平面或轴线平行的渗流为三维流平面或轴线平行的渗流为三维流。4.2 4.2 重力水运动的基本规律重力水运动的基本规律4.2.1 达西定律达西定律 1、实验原理和过程：、实验原理和过程： 渗流的基本定律是渗流的基本定律是法国水力学家达西（法国水力学家达西（Darcy） 于于1856年经年经过大量实验发现建立过大
12、量实验发现建立了地下水层流运动的了地下水层流运动的基本规律。基本规律。 达西实验装置图达西实验装置图 实验过程：实验过程： 实验结果：实验结果：单位时间内通过筒中砂的流量单位时间内通过筒中砂的流量Q与垂直与垂直水流方向的介质面积水流方向的介质面积F及上下测压管的水头差及上下测压管的水头差H成成正比，与渗透长度正比，与渗透长度L成反比。成反比。 Q = K FH /L 因为因为 ： V=Q / F I =H /L V=K I 4-3式式式中：式中：K-渗透系数（渗透系数（m/d），），V-渗透流速渗透流速（ m/d ），），I-水力坡度水力坡度 这就是著名的线性渗透定律这就是著名的线性渗透定律-
13、达西定律。达西定律。 渗透系数渗透系数(K)： 是表示岩土透水性的指标，一般与岩土和渗透是表示岩土透水性的指标，一般与岩土和渗透液体的物理性质有关的常数。液体的物理性质有关的常数。 渗透系数是含水层重要的水文地质参数之一。渗透系数是含水层重要的水文地质参数之一。 根据公式根据公式V=K I可知，当水力坡度可知，当水力坡度I=1时，渗时，渗透系数在数值上等于渗透流速。渗透系数透系数在数值上等于渗透流速。渗透系数K的单位的单位和和V的单位相同，以的单位相同，以m/s或或m/d表示。表示。松散岩石渗透系数参考值见表松散岩石渗透系数参考值见表4-1。岩土渗透岩土渗透 性分级性分级(GB50287-20
14、06)。2、达西定律的讨论与适用范围、达西定律的讨论与适用范围 1)达西定律的讨论达西定律的讨论 公式：公式：V=K I 反映通过任一断面一维流的反映通过任一断面一维流的渗流速度与其水力坡度之间的关系式，这种关系渗流速度与其水力坡度之间的关系式，这种关系的微分表达式：的微分表达式：V=Kdh/dL，对于二维流和三对于二维流和三维流同样适用。维流同样适用。 达西定律是在稳定运动条件下得到的，当地下达西定律是在稳定运动条件下得到的，当地下水为非稳定运动时，渗流中任一点处瞬时流速与水为非稳定运动时，渗流中任一点处瞬时流速与水力坡度之间的关系仍可用上面关系表征，只是水力坡度之间的关系仍可用上面关系表征
15、，只是渗透流速和水力坡度随时间在变化。渗透流速和水力坡度随时间在变化。（1）渗透流速： 根据水力学流速与流量的关系，上式可转化： Q = W V 与（2）式比较 V = KI （3）称为渗透流速（seepage velocity Darcy velocity specific discharge）上式为单位面积上的流量称比流量。由此看出，达西定律中：渗透流速与水力梯度的一次方成正比，故达西定律又称为线性渗透定律渗透流速（V）与过水断面（F） Q = K F I = F V 过水断面F，假想的断面 实际孔隙断面F n， n为孔隙度 实际水流断面F n ， n为有效孔隙度Q/F =V 比照水力学，
16、实际流速 Q/F= u关系：地下水渗透流速V = u ne渗透流速V：是假设水流通过整个岩层断面（骨架+空隙）时所具有的虚拟的平均流速。意义：研究水量时，只考虑水流通过的总量与平均流速，而不去追踪实际水质点的运移轨迹简化的研究过水断面：砂柱的横断面积，包括骨架和空隙在内的断面断面实际水流面积 ：扣除结合水所占据的范围以外的空隙面积实际断面（F）与过水断面（F）过水断面：砂柱的横断面积，包括骨架和空隙在内的断面过水断面：砂柱的横断面积，包括骨架和空隙在内的断面断面实际水流面积断面实际水流面积 ：扣除结合水所占据的范围以外的空隙面积：扣除结合水所占据的范围以外的空隙面积A 过水断面（水流可以穿越颗
17、粒） B 实际过水断面（水流只沿孔隙运动） 实际断面（ F)与过水断面（F）（2）水力梯度（I）（hydraulic gradient）水力学中水力坡度（J）：单位距离的水头损失沿渗透途径上的水头损失与相应的渗流长度之比。即：LhLHLHHI2121物理涵义上来看I：代表着渗流过程中，机械能的损失率，由水力学中水头的概念加以分析：在地下水渗流研究中任意点的水头表达式（2）水力梯度（I）（hydraulic gradient）在达西实验中：其原因是u2/2g 很小而忽略在地下水渗流研究中常：总水头 测压水头渗流过程中总机械能的损耗原因（与水力学相近）流体的粘滞性引起的内摩擦阻力（分子间）固体颗粒
18、表面对水流的反作用力（2）水力梯度（I）（hydraulic gradient）从达西公式: V = KI 来看：当I 增大时，V 也愈大；即流速V 愈大，单位渗流途径上损失的能量也愈大; 反过来，水力梯度I愈大时，驱动水流运动与速度也愈大注意：水头损失一定要与渗流途径相对应（3）渗透系数K（coefficient of permeability）也称为水力传导率（hydraulic conductivity）l定义：水力梯度为I =1 时的渗透流速（V=KI） 具有速度量纲L/ T（m/d；cm/s）l由公式V = K I 分析 当I一定时，岩层的K 愈大，则V 也愈大， Q 大因此，渗透系
19、数K 是表征岩石透水性的定量指标。（3）渗透系数K（coefficient of permeability）l影响因素： 以松散岩石,等径孔隙为例来分析，依据平行板水流理论可以推出：水的重率； 动力粘滞系数l从公式即得出： K与岩石性质有关 与流体物理性质有关l表4-1列出常见岩石渗透系数的参考值表4-1 松散岩石渗透系数参考值达西定律的适用范围达西定律的适用范围地下水必须保持层流运动才符合线性定律。地下水必须保持层流运动才符合线性定律。天然条件下地下水的渗流速度通常很缓慢，绝大部分为层天然条件下地下水的渗流速度通常很缓慢，绝大部分为层流运动，一般可用流运动，一般可用线性定律线性定律描述其运动
20、规律。描述其运动规律。近年来的实验表明，当地下水为流速较大的层流运动，近年来的实验表明，当地下水为流速较大的层流运动，V-V-I I开始偏离一直线，不符合达西定律。开始偏离一直线，不符合达西定律。雷诺数（雷诺数（ReRe）为）为1-10 1-10 的层流才符合达西定律。的层流才符合达西定律。l达西定律的适用范围达西定律的适用范围雷诺数数介于雷诺数数介于1 110 10 之间之间4) 非线性渗透定律非线性渗透定律 通过大量实验得知，当地下水为紊流运动通过大量实验得知，当地下水为紊流运动时，地下水运动服从时，地下水运动服从哲才定律哲才定律，即渗流速度，即渗流速度与水力坡度的平方根成正比：与水力坡度
21、的平方根成正比：V= KmI 1/2。 当地下水运动呈混合流状态时，符合公式当地下水运动呈混合流状态时，符合公式:V= KcI 1/m。 式中：式中： Kc为混合流时的渗透系数；为混合流时的渗透系数；m介于介于12之间。之间。4.2.2 地下水流向流速的测定地下水流向流速的测定 地地下下水水流流向向测测定定流向测定：流向测定：三点法（如图）；三点法（如图）；流速测定：流速测定：一般用试剂法。一般用试剂法。4.2.3 达西定律应用达西定律应用1、 地下水天然流量计算地下水天然流量计算Lhh2hhK2LhhKq21212122 h1h24-104-11式 完整井：完整井：水井贯穿整个含水层，在全部
22、含水井贯穿整个含水层，在全部含水层厚度上都安装过滤器并能全面进水的井。水层厚度上都安装过滤器并能全面进水的井。 非完整井：非完整井：井孔的进水段（过滤器）未穿井孔的进水段（过滤器）未穿透全部含水层，而只穿切含水层的一部分厚度。透全部含水层，而只穿切含水层的一部分厚度。2 、 地下水向完整井的稳定运动地下水向完整井的稳定运动1)、地下水向完整井的稳定运动、地下水向完整井的稳定运动地下水流向潜水完整井示意图地下水流向潜水完整井示意图 过水断面：过水断面： 水力坡度：水力坡度：rh 2 drdhI (1) 、潜水中的完整井、潜水中的完整井 降落漏斗、影响半径。降落漏斗、影响半径。 当遵循直线渗透定律
23、运动时，穿过该断面当遵循直线渗透定律运动时，穿过该断面的流通量为：的流通量为：rhdrdhKKIQ 2 orRlgohHK.Q223661 整理后得：整理后得：orRohHKQln22 或或又因：又因：所以上式还可表达为：所以上式还可表达为：ohHS orRnSSHKQl)2( orRlgS)SH2(K636. 1Q-潜水完整井裘布依公式潜水完整井裘布依公式4-12 对于承压含水层中的完整井，其出水量计对于承压含水层中的完整井，其出水量计算式为：算式为：orRMSKQlg732. 2 式中：式中：M为承压含水层厚度，为承压含水层厚度，m；其它符；其它符号同前。号同前。(2)、承压水中的完整井、
24、承压水中的完整井orRKMSQln2 - 承压水完整井裘布依公式承压水完整井裘布依公式4-13 根据前面公式，可用抽水试验资料求得含根据前面公式，可用抽水试验资料求得含水层的渗透系数水层的渗透系数K。如：。如：潜水完整井潜水完整井承压水完整井承压水完整井22lg732. 0oohHrRQK MSrRQKolg366. 0 SSHrRQKo)2(lg732. 0 or 当有一个观察孔和两个观察孔时：当有一个观察孔和两个观察孔时：orhhKQ101rln22 潜水完整井：潜水完整井：1212rln22rhhKQ 承压水完整井：承压水完整井：orKMSQ1rln2 12rln2rKMSQ 1)、井壁
25、进水的非完整井、井壁进水的非完整井承压水非完整井承压水非完整井完非QCQ1hLhhLC 241潜水非完整井潜水非完整井完非QCQ2式中式中 Q完完为潜水、承压完整井出水量；为潜水、承压完整井出水量；C1、C2分别为潜水、承压非完整井出水量折减系数。分别为潜水、承压非完整井出水量折减系数。MLMMLC 2423、地下水向非完整井的稳定运动、地下水向非完整井的稳定运动非完整井示意图非完整井示意图 这里讨论的结合水主要为结合水中的弱结合水。结这里讨论的结合水主要为结合水中的弱结合水。结合水运动以层流运动形式为主，但不服从牛顿的内摩合水运动以层流运动形式为主，但不服从牛顿的内摩擦定律，必须有外力克服结
26、合水所具有的抗剪强度后擦定律，必须有外力克服结合水所具有的抗剪强度后才能产生流动。才能产生流动。 其其V-I曲线为通过原点向曲线为通过原点向I轴凸出的曲线，即：轴凸出的曲线，即： V=K(I-I0)，其中，其中I0为起始水力坡度，其意义为克服为起始水力坡度，其意义为克服结合水的抗剪强度，使之发生流动所必须具有的水力坡度。结合水的抗剪强度，使之发生流动所必须具有的水力坡度。 上式中，对某一指定粘性土，上式中，对某一指定粘性土，K并非如重力水那样为常并非如重力水那样为常数，而是随数，而是随I的增大而增大，同样的增大而增大，同样I0也并非定值，也随也并非定值，也随I的增的增大而增大。大而增大。4.3
27、 4.3 结合水运动的基本规律结合水运动的基本规律 目前粘性土渗透流速目前粘性土渗透流速V与水力坡度与水力坡度I主要存在主要存在下列三种关系：下列三种关系： (1) V-I关系为通过原点的直线，服从达西定律；关系为通过原点的直线，服从达西定律； (2) V-I曲线不通过原点，水力梯度小于某一值曲线不通过原点，水力梯度小于某一值I0时无渗透；凸出的曲线，然后转为直线；时无渗透；凸出的曲线，然后转为直线； (3) V-I曲线通过原点，曲线通过原点，I小时曲线向小时曲线向I轴凸出，轴凸出，I大时为直线。大时为直线。4.4 4.4 饱水粘性土中水的运动规律饱水粘性土中水的运动规律 不能认为粘性土的渗透
28、特性及结合水的运动不能认为粘性土的渗透特性及结合水的运动规律已得出定论。规律已得出定论。饱水粘性土渗透试验的各类饱水粘性土渗透试验的各类V-I关系曲线关系曲线隐渗流隐渗流 显渗流显渗流4.5.1 毛细现象的实质毛细现象的实质 1、将细小的玻璃管插入水中，水会在管中上升到一、将细小的玻璃管插入水中，水会在管中上升到一定高度才停止，这便是固、液气三相界面上产生的定高度才停止，这便是固、液气三相界面上产生的毛细现象毛细现象。 毛细现象毛细现象 毛细现象的产生与表面张力有关。毛细现象的产生与表面张力有关。 设想在液面上划一根长度为设想在液面上划一根长度为L的线段，此的线段，此线段两边的液面，以一定的力
29、线段两边的液面，以一定的力f相互吸引，力相互吸引，力的作用方向平行于液面而与此线段垂直，大小的作用方向平行于液面而与此线段垂直，大小与线段长度与线段长度L成正比，即成正比，即fL；称为称为表面张表面张力系数力系数，单位为，单位为dyncm。 (1dyn=l05N) 。 由于表面张力的作用，弯曲的液面对液面由于表面张力的作用，弯曲的液面对液面以内的液体产生附加表面压强。以内的液体产生附加表面压强。 附加表面压强总是指向液体表面的曲率中附加表面压强总是指向液体表面的曲率中心方向：心方向： 凸起的弯液面，对液面内侧的液体，附加凸起的弯液面，对液面内侧的液体，附加一个正的表面压强；一个正的表面压强；
30、凹进的弯液面，对液面内侧的液体，附加凹进的弯液面，对液面内侧的液体，附加一个负的表面压强。一个负的表面压强。 表面张力表面张力f=d；液滴所受重力为；液滴所受重力为G=mg。 液滴开始断落时表面张力等于液滴所受重力。液滴开始断落时表面张力等于液滴所受重力。即：即： d=mg =mg/d式中：式中：d为液滴开始断落时液束最细处的直径。为液滴开始断落时液束最细处的直径。估测表面张力系数估测表面张力系数： 附加表面压强是如何引起的？附加表面压强是如何引起的？为了方便起为了方便起见，设想切取一个半径为见，设想切取一个半径为R的半圆球形液面的半圆球形液面(如如下图下图)。显然，在此液面的圆周状边线上都存
31、在。显然，在此液面的圆周状边线上都存在着指向液层内部的表面张力。着指向液层内部的表面张力。 其合力为其合力为2R，垂，垂直于面积为直于面积为R2的投影的投影圆面。由此，表面张力圆面。由此，表面张力所引起的附加表面压强所引起的附加表面压强Pc为：为： Pc =2R/R2 =2/R 因此，此时弯液面下的液体实际承受到因此，此时弯液面下的液体实际承受到的 表 面 压 强的 表 面 压 强 ( 以 下 简 称以 下 简 称 “ 实 际 表 面 压实 际 表 面 压强强)P=P0+Pc，此中，此中P0为大气压强。为大气压强。 任何形状的弯液面所产生的附加表面压强任何形状的弯液面所产生的附加表面压强Pc都
32、可以用拉普拉斯公式表示：都可以用拉普拉斯公式表示： Pc(1/R1+1/R2) 式中：式中：表面张力系数；表面张力系数； R1、R2液体表面的两个主要曲率半液体表面的两个主要曲率半径。径。拉普拉斯公式的函义：拉普拉斯公式的函义： 弯曲的液面将产生一个指向液面凹弯曲的液面将产生一个指向液面凹侧的附加表面压强，附加表面压强与表侧的附加表面压强，附加表面压强与表面张力系数成正比，与表面的曲率半径面张力系数成正比，与表面的曲率半径成反比。成反比。 附加表面压强存在的示意图附加表面压强存在的示意图 凸液面内压强大于液面外的压强；凹液面凸液面内压强大于液面外的压强；凹液面的压强小于液面外的压强。的压强小于
33、液面外的压强。4.5.2 毛细负压毛细负压 凹形弯液面产生的附加压强凹形弯液面产生的附加压强Pc，是个，是个负压强，称为负压强，称为毛细压强。毛细压强。 若将若将Pc换算为水柱高度换算为水柱高度(以以m为单位为单位)，并以并以hc表之，则：表之，则： hc=Pc/g =4/gD =0.03/D式中：式中：水的密度，等于水的密度，等于lgcm3； g重力加速度，等于重力加速度，等于981s2； 表面张力系数，取表面张力系数，取74dyncm；D毛细管直径，单位为毛细管直径，单位为mm。上式称为上式称为茹林公式茹林公式。 压力水头的测定压力水头的测定 hc为毛细压力水头为毛细压力水头(毛细压头毛细
34、压头)，是一个，是一个负的压力水头。就象在饱水带用测压管测负的压力水头。就象在饱水带用测压管测定压力水头定压力水头(测压高度测压高度)一样，可以用张力计一样，可以用张力计测定包气带的毛细压力水头。测定包气带的毛细压力水头。饱水带和包气带的压力水头测定饱水带和包气带的压力水头测定 饱水带中任一点的水头值饱水带中任一点的水头值H可表示为：可表示为： H=Z+hP 包气带中任一点的水头值包气带中任一点的水头值H则为：则为： H=Z-hc 两式中：两式中：Z由指定基准面算起的位置高度由指定基准面算起的位置高度(位置水头位置水头)； hP测压高度测压高度(压力水头压力水头)； hc毛细负压毛细负压(由毛
35、细力引起的负的压力水由毛细力引起的负的压力水头头)。 若取潜水面为基准，则潜水面处任一点饱若取潜水面为基准，则潜水面处任一点饱水带水头值为：水带水头值为： H=Z+hP=0(Z=0，hP=0) 若包气带支持毛细水的弯液面位于潜若包气带支持毛细水的弯液面位于潜水面处，则该点上支持毛细水的水头值为：水面处，则该点上支持毛细水的水头值为：HZhc=0一一hchc 即比周围潜水面水头低即比周围潜水面水头低hc4.5.3 毛细上升高度与悬挂毛细水毛细上升高度与悬挂毛细水 在此水头差驱动下支持毛细水将上升。当在此水头差驱动下支持毛细水将上升。当支持毛细水弯液面上升到支持毛细水弯液面上升到hc处时，弯液面处
36、的处时，弯液面处的水头为：水头为： HZhc= hc一一hc0 此时支持毛细此时支持毛细水处于平衡状态不水处于平衡状态不再上升。最大毛细再上升。最大毛细上升高度即为上升高度即为hc。 表面张力的竖直方向上的表面张力的竖直方向上的合力为合力为f=2 r cos；液柱；液柱所受重力为所受重力为G=ghr2。当。当液面稳定，不在上升时，液液面稳定，不在上升时，液柱所受表面张力等于它所受柱所受表面张力等于它所受的重力，即：的重力，即： 2 r cos=ghr2 h= (2cos) /gr。最大毛细上升高度的影响因素最大毛细上升高度的影响因素毛细现象毛细现象由式由式h= (2cos) /gr可知：可知：
37、 毛细管中液面上升高度与液体的种类、毛细管中液面上升高度与液体的种类、组成毛细管材料及管的内径有关。组成毛细管材料及管的内径有关。土的最大毛细上升高度土的最大毛细上升高度 土的类型土的类型 最大毛细上升高度最大毛细上升高度() 粗砂粗砂 25中砂中砂1235细砂细砂3570粉砂粉砂70150粘性土粘性土200400悬挂毛细水及其受力状况悬挂毛细水及其受力状况 4.5.4 包气带水水分分布及运动包气带水水分分布及运动均质土包气带水分分布均质土包气带水分分布 在包气带中，毛细负压是含水量的函数：在包气带中，毛细负压是含水量的函数：hchc (W) 原因原因：随着含水量降低，毛细水退缩到：随着含水量
38、降低，毛细水退缩到孔隙更加细小处，弯液面的曲率增大孔隙更加细小处，弯液面的曲率增大(曲率半曲率半径变小径变小)，造成毛细负压的负值更大。，造成毛细负压的负值更大。 渗透系数渗透系数(K)随含水量降低而迅速变小，随含水量降低而迅速变小，K也是含水量的函数：也是含水量的函数： KK (W) 原因：原因： (1) 含水量降低，实际过水断面随之减少；含水量降低，实际过水断面随之减少； (2) 含水量降低，水流实际流动途径的弯曲含水量降低，水流实际流动途径的弯曲程度增加；程度增加； (3) 含水量降低，水流在更窄小的孔角通道含水量降低，水流在更窄小的孔角通道及孔隙中流动，阻力增加。及孔隙中流动，阻力增加
39、。 由于上述原因，渗透系数与含水量呈非线性由于上述原因，渗透系数与含水量呈非线性关系。关系。 包气带水的运动，同样可以用达西定律描述，包气带水的运动，同样可以用达西定律描述，但与饱水带的运动相比，有以下三点不同：但与饱水带的运动相比，有以下三点不同： (1) 饱水带只存在重力势，包气带同时存在重饱水带只存在重力势，包气带同时存在重力势与毛细势；力势与毛细势； (2) 饱水带任一点的压力水头是个定值，包气饱水带任一点的压力水头是个定值，包气带的压力水头则是含水量的函数；带的压力水头则是含水量的函数； (3) 饱水带的渗透系数是个定值，包气带的渗饱水带的渗透系数是个定值，包气带的渗透系数随含水量的
40、降低而变小。透系数随含水量的降低而变小。n流网的有关概念流网的有关概念n流网的特点流网的特点n流网的绘制流网的绘制n流网的应用流网的应用n 渗流场：地下水流动（运动）的空间渗流场：地下水流动（运动）的空间. . 流网是描述渗流场中地下水流动状况的有效工具流网是描述渗流场中地下水流动状况的有效工具. .n 流网：渗流场某一典型剖面或切面上，由一系列流网：渗流场某一典型剖面或切面上，由一系列等水头线等水头线（equipotential linesequipotential lines）与）与流线流线（flow linesflow lines）组成的网格，）组成的网格，称流网（称流网（flow ne
41、tflow net）. .n 流线：某时刻在渗流场中画出的一条空间曲线，该曲线上各个水流线：某时刻在渗流场中画出的一条空间曲线，该曲线上各个水质点的质点的流速方向流速方向都与这条曲线相切（某时刻各点流向的连线）都与这条曲线相切（某时刻各点流向的连线）n 迹线：流体水质点在渗流场中某一时间段内的运动轨迹迹线：流体水质点在渗流场中某一时间段内的运动轨迹trace linetrace line）（稳定流条件下流线与迹线重合）（稳定流条件下流线与迹线重合）n 等水头线：在某时刻，渗流场中水头相等各点的连线（水势场的等水头线：在某时刻，渗流场中水头相等各点的连线（水势场的分布）分布）（1 1）基本概念）
42、基本概念n 常见的二维流网图：常见的二维流网图： 平面流网：潜水等水位线图，承压水等测压水位线图平面流网：潜水等水位线图，承压水等测压水位线图 剖面流网：含水量厚度较大时，常需要刻画剖面的水流剖面流网：含水量厚度较大时，常需要刻画剖面的水流l 常见的边界类型常见的边界类型 第一类边界条件：给定水头边界条件，具有无限补给或排泄地第一类边界条件：给定水头边界条件，具有无限补给或排泄地下水的能力，如与地下水具有水力联系的地表河流、湖泊等；下水的能力，如与地下水具有水力联系的地表河流、湖泊等； 第二类边界条件：给定流量边界条件，典型的有隔水边界、地第二类边界条件：给定流量边界条件，典型的有隔水边界、地
43、下水分水岭。下水分水岭。 在各向同性介质中流线与等水头线正交，在各向同性介质中流线与等水头线正交， 在在各向异性介质中流线与等水头线斜交各向异性介质中流线与等水头线斜交 按一定规则绘制的：等水头线按一定规则绘制的：等水头线相邻两条等水相邻两条等水头线间的势差为常量，流线头线间的势差为常量，流线相邻两条流线间相邻两条流线间的通量为常量的通量为常量 等水头线的疏密代表水力坡度的大小，流线的等水头线的疏密代表水力坡度的大小，流线的疏密反映径流强度的大小。疏密反映径流强度的大小。l 等水头线（面）与流线（面）正交；等水头线（面）与流线（面）正交； l 等水头线（面）与流线（面）不是两个独立问题，知道一
44、方就可据等水头线（面）与流线（面）不是两个独立问题，知道一方就可据正交原则推求另一方。正交原则推求另一方。 l 正交网格中，每两条流线间的流量相等。正交网格中，每两条流线间的流量相等。等水头线等水头线水头降水头降各各向向同同性性介介质质Effect of Topography on Flow PatternsLow relief topographyModerate relief topographyHigh relief topographyEffect of Layers on Flow PatternsModerate K ratioHigh K ratioEffect of Dip o
45、n Flow PatternsDip towards recharge area high reliefDip towards discharge area low reliefEffect of Pinchouts on Flow PatternsHigh K below discharge areaHigh K below recharge area（3 3）流网绘制）流网绘制定量方法定量方法精确流网精确流网 定性方法定性方法信手流网信手流网边界条件：定水头边界、隔水边界、潜水面边界等边界条件：定水头边界、隔水边界、潜水面边界等源、汇：源源、汇：源发散流线，汇发散流线，汇吸收流线吸收流线分
46、水线：虚拟隔水边界分水线：虚拟隔水边界流线流线1、分析边界条件,根据边界性质画出易确定的等水头线和流线2、确定源、汇，控制流线的趋向3、画出渗流场周边流线4、中间内插，画其他流线 等单宽流量控制流线根数5、等水头差确定等水头线间隔6、按正交规则画等水头线l首先分析水文地质条件，搞清补给区、排泄区、或源汇项首先分析水文地质条件，搞清补给区、排泄区、或源汇项分布、边界条件等。分布、边界条件等。 l先绘制肯定的流线和等水头线：先绘制肯定的流线和等水头线： 隔水边界是流线隔水边界是流线 无入渗、无蒸发条件下潜水面是流线无入渗、无蒸发条件下潜水面是流线 湖泊、河流边界可看成等水头线，湖泊、河流边界可看成
47、等水头线， l有两个以上排泄点时应确定分水线、面、点有两个以上排泄点时应确定分水线、面、点。河间地块流网的绘制1）寻找已知边界（湿周，隔水边界，水位线） 2）分水线、源、汇的确定 3）画出渗流场周边流线与条件 4）确定等水头值，中间内插（3 3）流网绘制）流网绘制定性流网的绘制（各向同性介质中） 在许多实际工作中，绘制定性流网分析问题很重要 精确流网受许多条件（资料不足等）制约，很难办到 思考：绘制流网需要考虑渗流场的哪些条件？ l用数学方法求解运动方程：求解空间水头分布绘制等用数学方法求解运动方程：求解空间水头分布绘制等水头线；求解流函数绘制流线。水头线；求解流函数绘制流线。 l物理模型模拟
48、：水电比拟。物理模型模拟：水电比拟。 l现场测定（测定水头分布，绘制等水头线，再据正交现场测定（测定水头分布，绘制等水头线，再据正交原则绘制流线）原则绘制流线） l信手流网：据流场边界性质和介质特性，半定量地绘信手流网：据流场边界性质和介质特性，半定量地绘制流网。制流网。zk144zk145zk143zk142zk18zk32zk26zk16zk170zk160zk34zk08zk20zk12zk13zk105zk169zk23zk05zk07zk09zk17zk155zk26zk18zk32zk159zk34zk160zk28zk30zk153zk02zk169zk01zk07zk171zk
49、152zk20zk06zk139n 解释水文地质现象 n 判断地下水系统内部结构 n 分析地下水的补给、排泄、径流特征 n 计算渗流场任意点的水头、压强、水力坡度、渗透流速等 n 据流网选择垃圾填埋场位置等AB1、确定水头值H2、确定水力梯度I 3、确定渗透流速V4、确定流量Q5、了解水质点的渗流途径及长短（当流线与迹线重合，流线近视为水质点的运移轨迹）？比较图中A、B两点的H、I、V的大小； 在何处打井，井水不受污染？它反映了渗流场中地下水的流动状况，同时也是介质场与势场的综合反映，提供这两方面的信息。地表水体地表水体定水头边界定水头边界: :河流湿周为等水头线河流湿周为等水头线隔水边界隔水
50、边界零流量边界：零流量边界：流线流线潜水面边界潜水面边界 稳定的侧向补给：流线稳定的侧向补给：流线 入渗补给：既不是流线入渗补给：既不是流线也不是等水头线也不是等水头线等水头线、流线与各类边界的关系思考题：河流完全切割含水层至隔水底板，其它条件不 变时流网形态？实验：达西渗流实验（演示）实验：达西渗流实验（演示） 实验目的：1通过稳定流条件下的渗流实验，进一步理解渗流基本定律达西定律。2加深理解渗流速度、水力梯度、渗透系数之间的关系，并熟悉实验室测定渗透系数的方法。实验内容：1了解达西实验装置。2根据达西公式： 测定不同试样的渗透系数 K。式中： Q 渗透流量； 过水断面面积； H 上下游过水
