地下水动力学—地下水水文学(持续更新)
部分内容转载自南方科技大学梁修雨教授地下水水文学讲课PPT
表面张力:水等液体会产生使表面尽可能缩小的力。
水文地质学英文及其含义
英文 | 中文 |
---|---|
stream | 流动,流出,小河,溪流 |
flow | 流动 |
streamflow | 河道汇流部分,单位m^3/s |
runoff | 径流深,单位:mm |
discharge | 流量 |
baseflow | 基流 |
evaporation | 蒸发 |
recharge | 补给 |
rainfall | 降雨 |
catchment | 流域、集水处、聚水盆地 |
throughfall | 穿透雨 |
leakage | 渗漏 |
infiltration | 入渗 |
freshwater | 淡水 |
continental | 大陆的,大洲的 |
vadose zone | 包气带、渗流带 |
overflow | 溢流 |
overland flow | 地表径流 |
diffuse | 四散的,散开的,弥散的 |
watershed | 流域、分水岭、集水区 |
water table | 潜水面 |
subsurface | 次表面的,地下的 |
the lumped parameter | 集总参数 |
drainage basin/catchment | 流域 |
watershed<br / divide |
分水岭 |
creek | 小溪,小河 |
outflow constant | |
recession constant | 退水常数 |
specific yield | Sy 给水度 |
specific storage | Ss |
tributary | 支流 |
fluvial | 河的、河流的、河中的 |
S–specific capitiy | 贮水系数 |
Ss–specific storage | 贮水率 |
T–transmissivity | 导水系数 |
Sy–specific yeild | 给水度 |
D- hydraulic diffusivity | 水力扩散系数 |
K–Permeability | 渗透系数 |
傅里叶变化的基础 | |
$$ | |
\omega=2\pi f | |
$$ | |
ω \omega ω是角频率,f是频率 |
土壤水分特征曲线
土壤水分特征曲线是土壤含水量 θ \theta θ与土壤水吸力 ϕ \phi ϕ的关系曲线该曲线能够间接反映土壤孔隙大小的分布,分析不同质地土壤的持水性等。目前负压计法是测量土壤水吸力最简单最直观的方法,而时域反射法是测量土壤体积含水率的最常用最便捷的方法。
基模势是土壤含水量的函数,干燥土壤,吸力最大,随着土壤含水量的增加,吸力逐渐减小,当土壤含水量达到田间持水量时,吸力变为0,通常基模势的负值为吸力 ϕ \phi ϕ,同样的吸力,在脱水和吸水两个过程中,土壤吸持的水分数量是不同的,脱水过程过程吸持的水分要大于吸水过程吸持的水分,或者说,当土壤含水量一定时,脱水过程中相应的土壤吸力要大于吸水过程中相应的土壤吸力。
将土壤含水量和土壤吸力绘制成散点图,然后用经验公式拟合
下渗率:单位时间通过单位面积的土壤层面渗入到土壤的水量称为下渗率,常用单位为mm/min,mm/h。
下渗容量:
下渗容量随时间的变化曲线称为下渗曲线。
饱和土壤的达西定量
v = K I v=KI v=KI
渗透系数也称水力传导系数,渗透系数具有速度的量纲,常用cm/s或m/d
控制非饱和土壤水流运动的因素与控制饱和土壤水流运动的的因素有所不同,一是两者的总势组成不同,在饱和土壤中,总势由重力势和静水压力势组成,而在非饱和土壤中,总势则由重力势和基模势组成;二者的水力传导度不同,饱和水力传导度是一个常数,非饱和水力传导度是土壤含水量的函数,干燥土壤的水力传导度最小,随着土壤含水量的增加,水力传导度也增加,当土壤含水量达到饱和时,水力传导度也达到最大,此时即为饱和水力传导度。
线性和非线性偏微分方程
如果偏微分方程中,未知函数及它的所有偏导数都是线性的,且方程中的系数都依赖于自变量(或者是常数),那么这样的偏微分方程就称为线性偏微分方程,特别的,如果方程中的系数都是常数,则称为常系数偏微分方程,显然,如果方程中的系数是自变量的函数,则称为变系数偏微分方程,方程中如果出现未知函数及偏导数不是线性的,则称为非线性偏微分方程。
各阶微分项有次数高于1的微分方程即为非线性偏微分方程。
扩散率为常数时的下渗方程:
{ ∂ θ ∂ t = D ∂ 2 θ ∂ Z 2 θ ( z , 0 ) = θ 0 θ ( 0 , t ) = θ n lim z → ∞ θ ( z , t ) = θ 0 \left\{ \begin{aligned} \frac{\partial{\theta}}{\partial{t}}=D\frac{\partial^2{\theta}}{\partial{Z^2}} \\ \theta(z,0)=\theta_0 \\ \theta(0,t)=\theta_n\\ \lim_{z\to\infty}{\theta(z,t)}=\theta_0 \end{aligned} \right. ⎩
⎨
⎧∂t∂θ=D∂Z2∂2θθ(z,0)=θ0θ(0,t)=θnz→∞limθ(z,t)=θ0
上式可以使用laplace变换进行求解。
土壤中三相关系
V t V_t Vt为土块的总体积
V s V_s Vs为固相体积
V w V_w Vw为液相体积
V a V_a Va为气相体积
V f V_f Vf土壤中孔隙的容积
则有:
V f = V w + V a V_f=V_w+V_a Vf=Vw+Va
V t = V s + V w + V a V_t=V_s+V_w+V_a Vt=Vs+Vw+Va
M t M_t Mt为土块的总质量
M s M_s Ms为固相质量
M w M_w Mw为液相质量
M a M_a Ma为气相质量,其中 M a = 0 M_a=0 Ma=0
与土壤三相有关的土壤物理量
固体密度:土壤中固体颗粒的质量与固相容积的比值
干容重:土壤中固体颗粒的质量与土块总体积的比值
孔隙度:土块中孔隙容积与土块总容积之比
孔隙比:土块中孔隙容积与固体颗粒容积之比
土壤含水量
土壤含水量描述土壤中含有水分的情况,又称土壤湿度,他有几种表示方法,
- 质量含水率
- 容积含水量
- 饱和度
土壤水的存在形态
分子力**
土壤水分常数
最大吸湿量(吸湿系数):在水汽达到饱和的空气中,干燥土壤的吸湿水达到最大数量时的土壤水量
最大分子持水量:
凋萎系数
毛管断裂含水量
田间持水量
饱和含水量
土壤水分剖面
土壤含水率沿深度方向的变化曲线称为土壤水分剖面
它描述了土壤含水率在深度方向的分布情况,故又称土壤含水率垂向分布,根据土壤水分剖面,可以计算出土壤中任一土层,以水深计的含水量。
土壤水分剖面图
下渗曲线
下渗现象的定量表示是下渗率,单位时间通过单位面积的土壤层面渗入到土壤中的水量称为下渗率。单位mm/min或mm/h
如果供水强度充分大,则下渗速率将达到同初始土壤含水量和同土壤质地,结构条件下的最大值,称为下渗能力或下渗容量,
下渗容量随时间的变化曲线称为下渗曲线。
下渗曲线是一条递减曲线,分为
- 渗润阶段:土壤含水量较小,下渗容量较大,下渗容量随时间迅速递减
- 渗漏阶段:土壤含水量不断增加,下渗容量明显减小,下渗容量随时间递减变得缓慢
- 渗透阶段:土壤含水量达到田间持水量以上,下渗容量变得稳定,达到下渗容量最小值,为稳定下渗率。
土壤水分剖面的四分带
最上层为饱和带,一般厚度不大,不到1.5cm,饱和带以下为水分传递带,这是一个土壤含水量沿深度分布比较均匀,、厚度较大的非饱和土层,厚度随供水是时间的增长而不断增加,土壤含水量介于田间持水量和饱和含水量之间,水分传递带以下为湿润带,他是连接水分传递带和湿润峰的水分带。在这一带中,土壤含水量随深度迅速减小,湿润带与下渗水未涉及到的土壤的交界面称为湿润峰。在湿润峰处,土壤含水量梯度很大。
控制非饱和土壤水流运动的因素与控制饱和土壤水流运动的因素有所不同。一是二者的总势组成不同,在饱和土壤中,总势由重力势和静水压力势组成,而在非饱和土壤中,总势由重力势和基模势组成;二是二者的水力传导度不同,饱和水力传导度是一个常数,非饱和水力传导度是土壤含数量的函数,干燥土壤的水力传导度最小,随着土壤含水量的增加,水力传导度也增加,当土壤含水量达到饱和时,水力传导度也达到最大,此即为饱和水力传导度。
对于理查兹方程,只要已知土壤水分特征曲线 Ψ − θ \Psi-\theta Ψ−θ,和水力传导度与土壤含水量的关系曲线 K ( θ ) − θ K(\theta)-\theta K(θ)−θ
并给定适当的初始条件和边界条件,就可由Richards方程求解各种非饱和水流运动问题。
拉普拉斯变换
F ( s ) = L { f ( t ) } = ∫ 0 + ∞ f ( t ) e − s t d t F(s)=\mathcal{L}\left\{
{f(t)} \right\} =\int_{0}^{+\infty}{f(t)e^{-st}dt} F(s)=L{
f(t)}=∫0+∞f(t)e−stdt
使用拉普拉斯变换求解常微分方程
x ′ ′ ′ + 3 x ′ ′ + 3 x ′ + x = 6 e − 1 x^{'''}+3x^{''}+3x^{'}+x=6e^{-1} x′′′+3x′′+3x′+x=6e−1
玻尔兹曼变换解微分方程
水循环(hydrologic cycle)
降水落到地表后进入到不同的水循环途径,一部分暂时以冰或雪的形式留在地表或汇入洼地,后者称洼地蓄水(depression storage),部分河流中的水来自地表径流或者渗入河床的地下水,河流中来自地下水的部分称为基流(baseflow)
导水系数:含水层的渗透系数与其厚度的乘积
含水层响应时间:
the characteristic aquifer response time
是指由降雨到含水层接受补给的时间,在干旱区特征响应时间可能会特别长。
系统在单位冲激函数激励下引起的零状态响应被称之为该系统的冲激响应,他与系统的传递函数互为傅里叶变换的关系。
时间尺度
水分特征曲线斜率的倒数,表示单位基质势的变化所引起的土壤含水率的变化。称为容水度。
C = d θ d ϕ m C=\frac{d\theta}{d\phi_m} C=dϕmdθ
田间持水量:在长时间的重力排水后仍然保留在土中的水量
土壤含水率 θ \theta θ与非饱和导水率 K K K、垂向位置z,土壤基质势 Φ \Phi Φ,土壤重力势 h h h,时间t的相互关系
- 土壤含水率是位置和时间的函数——
- 非饱和导水率是位置和时间的函数——
- 土壤基质势是——
- 土壤重力势是———
含水层的释水能力为S,(storativity),分为specific storage 和 specific yield
给水度:单位面积的含水层,当潜水面下降一个单位长度时,在重力作用下所能释放出的水量
比蓄水量: 是指蓄水层的一部分在保持完全饱和的情况下,每单位质量或体积的蓄水层,每单位水头变化释放的水量。
承压含水层基本微分方程
throughfall净降水量
baseflow 激流
streamflow 径流,河流流量,河川流量
runoff flow径流
1. Darcy’s law
hydraulic head
hydraulic gradient
seepage velocity 渗透速率
Q = − K A I = − K A d h d l Q=-KAI=-KA\frac{dh}{dl} Q=−KAI=−KAdldh
v = K I v=KI v=KI
gravity 重力势
pressure gradient 压力梯度
kinetic energy 动能
gravitational potential energy
fluid pressural energy
理想流体,管道流,伯努利方程
Φ = g z + P ρ + v 2 2 \Phi=gz+\frac{P}{\rho}+\frac{v^2}{2} Φ=gz+ρP+2v2
把机械能除以g,单位重量的机械能
Φ = z + P ρ g + v 2 2 g \Phi=z+\frac{P}{\rho g}+\frac{v^2}{2g} Φ=z+ρgP+2gv2
elevation head
pressure head
velocity head