地下水动力学地下水流基本微分方程及定解条件(1)

地下水弹性储存
弹性储存：当地下水水头（水压）降低(或升高）时， 含水层、弱透水层释放（或储存）地下水的性质
物理意义： ➢ 弹性储存与重力储存不同；
给水机制不同
➢ 弹性储存更宜理解为“变形储存”；
➢ 弹性储存这种性质不仅承压含水层具备，层间 弱透水层也有弹性储存
因Vs=constant，故
只在垂直方向上有压缩，
单元体内地下水 质量变化量
m V nxyz
m [( nz) |(x,y,z,tt) (nz) |(x,y,z,t) ]xy
X方向流入流出差
y方向流入流出差
z方向流入流出差 单元体内地下水 质量变化量 地下水连续性方程
( vx ) |(x, y,z,t) ( vx ) |(xx, y,z,t) xyzt
H x
)]
由于在一般情况下，水的密度变化很小，可视近似不变，故
x
K xx
H x
x
K
xx
H x
(vx
x
)
x
(K xx
H x
)
渗流连续性 方程化简
(
v
x
x
)
( v y
y
)
( v z
z
) xyz
(nz)
t
xy
（二）化简方程左端项
(nz) z ( n ) H
t
t
vx
K xx
H x
vy
Vb Vs Vv
Vv＝nVb；
由于骨架部分体积不变 dVb dVv
Vs＝(1-n)Vb
1 dVs 1 dVv 1 n dVs n dVv
Vb d Vb d
Vs d Vv d
Vv＝nVb；Vs＝(1-n)Vb
式中
——多孔介质固体颗粒压缩系数，表示多
孔介质中固体颗粒本身的压缩性的指标，s<<p；
➢ 必须区分两者之间的不同，潜水含水层还存在滞后疏干现象。
承压含水层抽水时，水的释放是由于压力减少造成的，这 一过程是瞬时完成的。只要水头下降不低到隔水顶板以下，水 头降低只引起含水层的弹性释水，可用贮水系数*表示这种释 水的能力。
导压系数
描述含水层水头变化的传导速度的参数，其数值等于含水层 的导水系数与贮水系数之比或渗透系数与贮水率之比。
由于
Vm
V p
V0
V p
V
d(m)
dV V
m
d( 1 ) d
水的压缩方程
dp 1 d
d
dp
(1 5)
多孔介质的压缩方程
假定多孔介质近似地符合弹性变形，依虎克定律，有
1 dVb d Vb
α为岩土的体积弹性压缩系数。
如果上部荷载不变，则 d dp
dp 1 dVb
Vb
的参数，在地下水动力学计算中具有重要的意义。
➢ 贮水率
表示含水层水头变化一个单
位时，从单位体积含水层中，因水体积膨胀（压缩）以
及骨架的压缩（或伸长）而释放（或储存）的弹性水量。
单位1/L。
➢ 贮水系数又称释水系数或储水系数，为含水层水头变化 一个单位时，从底面积为一个单位，高度等于含水层厚 度的柱体中所释放（或贮存）的水量；指面积为一个单 位、厚度为含水层全厚度M的含水层柱体中，当水头改 变一个单位时弹性释放或贮存的水量，无量纲。既适用 于承压含水层，也适用于潜水含水层。
量纲为L2T-1。
2.2 渗流连续性方程
连续性方程就是质量守恒方程，也称为水均衡方程 水均衡的基本思想：
对某一研究对象，流入－ 流出＝V 研究对象可以是大区域的，也可以是微分单元体
大区域的水均衡计算经常用于区域的水资源评价 本课程基于微分单元体做水均衡，推导渗流连续性方程。
为反映含水层地下水运动的普遍规律，我们选定在各向 异性多孔介质中建立地下三维不稳定流动连续性方程。
第二章 地下水流基本微分方程及定解条件
基本理论：连续性假设+达西定律+水均衡原理
➢ 对各种水流问题建立基本微分方程及数学模型： ●按空间维数：一维、二维（平面二维、剖面二维）、三维 ● 按含水层类型：承压水流、潜水流、多层（越流联系）等
➢ 求解数学模型（利用解析法），得到一些典型解析解： ●裘布依稳定井流模型 ●无越流承压含水层中的完整井流（泰斯模型） ●无越流潜水含水层中的完整井流（博尔顿模型-考虑滞后给水、 纽曼模型-考虑流速垂直分量和弹性储量） ●越流系统中的承压完整井流模型
➢ 应用： ●预测抽水条件下的水头变化； ●利用抽水试验资料求含水层参数。
第二章 地下水流基本微分方程及定解条件
教学目标：
➢ 准确理解渗流连续性概念 ➢ 掌握达西定律和质量守恒原理的应用 ➢ 掌握建立地下水基本微分方程的思想方法 ➢ 几种典型的地下水流方程的推导
●潜水剖面二维流、平面二维流 ●承压水二维流 ● 三维流 ➢ 边界条件概化，初始条件确定方法与原则 ➢ 能够用数学模型描述实际问题
➢ 与水体积膨胀所释放出的水量（dV）之和
上述二者之和所释放出的水量为
或
（1-14）
式中 s ——贮水率[释水率]（specific storativity），量纲
[L-1]，为弹性释水[贮水] ； 式中 M——含水层厚度(m)；
*——贮水系数（storativity）。
*=sM ➢ 贮水系数*和贮水率s都是表示含水层弹性释水能力
即：1-65式变为：
渗流连续性方 程化简
(v
x
x
)
(vy )
y
( v z
z
)
xyz
(nz)
t
xy
（二）化简方程左端项
当渗流满足达西定律，且取坐标与各向异性主轴方向一致，有
vx
K xx
H x
vy
K yy
H y
vz
K zz
H z
(vx )xFra bibliotekx(K xx
H x
)
[
x
K xx
H x
x
(K xx
第二章 地下水流基本微分方程及定解条件
主要内容：
➢ 建立连续性方程 ➢ 分析含水层与岩石、流体压缩性关系 ➢ 建立不同含水层地下水流微分方程 ➢ 讨论边界条件及初始条件 ➢ 用数学模型描述实际问题
2.1 水和多孔介质的压缩性
水的压缩方程（地下水的状态方程)
假定水近似地符合弹性变形，依虎克定律，有
➢ 贮水率是描述地下水三维非稳定流或剖面二维流中的 水文地质参数，既适用于承压水也适用于潜水。对于平 面二维非稳定流地下水运动，当研究整个含水层厚度上 的释水情况时，用贮水系数来体现。
上述两参数之间的不同，还在于潜水含水层存在滞后疏干现 象。
➢ 弹性释水与重力给水: 对于含水层而言，由于受埋藏条件的 限制，抽水时，水的给出存在着不同。
...
由于很小，且p变化不大，故
e ( p0 p) 1 ( p0 p)
V V0
1 ( p0
p)
V V0[1 ( p0 p)]
V0 V0 ( p0 p)
V V0 V0 ( p0 p)
V V0 V0
( p0
p)
V V0
( p0
p)
水的压缩方程
dp 1 dV
V
由于V~V0变化不大，故
——多孔介质中孔隙压缩系数 （Compressibility of the pores of a porous medium），表 示多孔介质中孔隙的压缩性的指标。n——多孔介质的孔隙度。
1-8
因
，故
。
1-9
水的压缩方程 多孔介质的压缩方程
dp 1 dV
V V p
V
1 dVb d Vb
嫁到多孔介质固体骨架上，增大有效应力，压缩多孔介质，结果使含 水层介质厚度变薄和空隙率n变小，同时从孔隙中释放地下水； ➢ p减少多孔介质固体颗粒也会膨胀,而有效应力增大又会影响固体颗 粒的变形。综合起来，这种现象比较复杂。考虑到固体颗粒的压缩性 比多孔介质要小得多，因此通常忽略多孔介质固体颗粒的压缩性。
渗流连续性方程推导
X方向流入流出差
(vx ) |(x,y,z,t) yzt (vx ) |(xx,y,z,t) yzt
y方向流入流出差
(v y ) |(x,y,z,t) xzt ( v y ) |(x,yy,z,t) xzt
z方向流入流出差 ( v z ) |( x, y,z,t) xyt ( v z ) |( x, y,zz,t) xyt
e ( p0 p) V V0
水的压缩方程
按Taylor级数展开
f (x) f (x) f (0) x f (0) x2 f (n) (0) xn
1!
2!
n!
ex 1 x x2 x3 ... 2! 3!
e ( p0 p)
1
( p0
p)
2
2!
(
p0
p)2
3
3!
( p0
p)3
渗流连续性方程推导
图2-1-1多孔介质单元水均衡要素图
X方向流入 X方向流出
假设：水是可压缩的，多孔介质骨 架在垂直方向可压缩，但在水平方 向不可变形。
均衡的含义：在t时段内从x,y,z三 个方向共6个单元界面上流入流出 水的净总质量等于单元体内储存量 的变化。
V Qt vt
m V vt
X方向流入流出差 (vx ) |( x, y,z,t) yzt ( vx ) |( xx, y,z,t) yzt
(v
x
x
)
(vy )
y
(vz
z
)
xyz
(nz)
t
xy
渗流连续性方程化简