第三讲 地下水系统分布参数数学模型
合集下载
第三讲 地下水数值模拟原理及建模方法和步骤_xiugai
➢ 方法二．在渗流区剖分的基础上，直 接由达西定律和水均衡原理，建立各 个均衡区的水均衡方程，从而得到差 分方程。适用于矩形网格、三角形网 格。
矩形网格 多边形网格
方法一：差商代替微商
1、网格划分的基本类型
（1）先划格线，格点位
均
于网格中心
衡
网
格
（2）先规定格点位置，
再垂直平分两相邻结点的连
节
线作格线，形成的网格即为 水均衡区
H (x0
x, t0 )
2H (x0 , t0 ) (x)2
H (x0
x, t0 )
一维控制方程差分格式
方法一
控制方程
T
2h( x, t ) x2
h( x, t ) t
网格剖分nx个
显式差分格式
H (x0 x, t0 ) 2H (x0 , t0 ) H (x0 x, t0 ) H (x0 , t0 t) H (x0 ,t0 )
(x) 2
t
Hn i 1
2
H
n i
Hn i 1
(x)2
H
n1 i
H
n i
t
i 2,3,, nx 1
隐式差分格式
H (x0 x, t0 t) 2H (x0 , t0 t) H (x0 x, t0 t) H (x0 , t0 t) H (x0 , t0 )
(x) 2
t
H n1 iБайду номын сангаас1
数值方法很多，但是最简单实用的是有限差分法：
✓ 有限差分法
✓ 有限单元法
✓ 积分有限差分法
✓ 半解析半数值法
✓ 边界元法
✓ 有限体积法
只讲有限差分法
一、有限差分法的基本原理
矩形网格 多边形网格
方法一：差商代替微商
1、网格划分的基本类型
（1）先划格线，格点位
均
于网格中心
衡
网
格
（2）先规定格点位置，
再垂直平分两相邻结点的连
节
线作格线，形成的网格即为 水均衡区
H (x0
x, t0 )
2H (x0 , t0 ) (x)2
H (x0
x, t0 )
一维控制方程差分格式
方法一
控制方程
T
2h( x, t ) x2
h( x, t ) t
网格剖分nx个
显式差分格式
H (x0 x, t0 ) 2H (x0 , t0 ) H (x0 x, t0 ) H (x0 , t0 t) H (x0 ,t0 )
(x) 2
t
Hn i 1
2
H
n i
Hn i 1
(x)2
H
n1 i
H
n i
t
i 2,3,, nx 1
隐式差分格式
H (x0 x, t0 t) 2H (x0 , t0 t) H (x0 x, t0 t) H (x0 , t0 t) H (x0 , t0 )
(x) 2
t
H n1 iБайду номын сангаас1
数值方法很多，但是最简单实用的是有限差分法：
✓ 有限差分法
✓ 有限单元法
✓ 积分有限差分法
✓ 半解析半数值法
✓ 边界元法
✓ 有限体积法
只讲有限差分法
一、有限差分法的基本原理
《地下水数值模拟》课件
CHAPTER 04
地下水数值模拟的案例分析
案例一：某地区地下水污染模拟
总结词
该案例展示了如何运用地下水数值模拟技术 预测和评估某地区地下水污染情况。
详细描述
该案例首先介绍了该地区的地下水分布和流 向，然后通过建立数值模型，模拟了不同污 染源对地下水的影响，并预测了污染扩散的 范围和程度。最后，根据模拟结果，提出了 相应的污染防治措施。
VS
有限体积法适用于不规则的网格系统 和复杂的边界条件，能够得到相对准 确的结果，计算量适中，适用于较大 的模型规模。
CHAPTER 03
地下水数值模拟的步骤
建立数学模型
01
确定研究区域和边界条件
02
描述地下水流动和物质传输过程
03
建立数学方程，包括连续性方程、动量方程、源汇 项等
模型离散化
1
地下水数值模拟的应用
地下水数值模拟广泛应用于水资源管理、环境保护、地质 灾害防治等领域。
通过模拟地下水动态变化，可以预测未来地下水资源量、 评估地下水污染风险、研究地下水与地质灾害的关系等， 为相关决策提供科学依据。
CHAPTER 02
地下水数值模拟的基本方法
有限差分法
有限差分法是一种将偏微分方程离散 化为差分方程的方法，通过在时间和 空间上将偏微分方程近似为差分方程 ，从而将连续的物理量离散化为离散 的数值。
随着数值计算技术的发展，地下水数值模型将越来越复杂，能够 模拟更多的物理过程和化学反应。
参数优化和数据同化
通过人工智能和机器学习技术，对模型参数进行自动优化和数据同 化，提高模拟精度和可靠性。
多尺度模拟
从微观到宏观的多尺度模拟将成为一个重要方向，能够更好地揭示 地下水系统的复杂性和规律性。
地下水数值模型设计步骤及对资料的要求
一、数值模型研究一般步骤（续）
建立数值模型
网格剖分：根据确定的数值方法和软件，对研究区进行剖分。 对于平面二维流问题，一般将研究区剖分成矩形或三角形网格； 对于剖面二维流问题，一般也是将剖面区域剖分成矩形或三角形； 对于三维流问题，一般先在垂向上分成若干层，而在每层剖分成矩 形或三角形。
边界条件： 初始条件： 含水层参数：渗透系数、储水系数、给水度、孔隙度 源汇项：降雨入渗、河流补给、蒸发排泄、地表水体、沟渠渗漏、
上述数学模型是一个偏微分方程定解问题，通常只能用数值方法 求解，常用的数值方法有：有限差分法和有限单元法。
目前有一些软件可以直接用于求解地下水流动问题，如果不是自 己编写程序，可以选择合适的软件，建立数值模型。
如果用软件，则需对软件功能作简要介绍，论述软件的适用性。
3
第七章 数值模型一般步骤及对 勘查资料的要求
模型设计者应出具有较高理论水平和丰富经验的水文地质工作者担 任。设计者应精细地分析有关资料, 以获得较符合实际条件的分区 图。
18
三、数值模型设计中一些特殊问题
（一）抽水试验设计 （二）抽水试验数值模拟设计 （一）含水层剖分注意事项
（一）抽水试验设计
19
三、数值模型设计中一些特殊问题
（一）含水层剖分注意事项
12
二、数值模型设计及对资料和水文地质勘探的要求
（四）边界条件的确定
边界条件是与计算区的范围同时确定的
考虑计算区的范围时, 必须同时确定边界条件的性质； 反之, 边界的位置一旦确定, 计算区的范围自然也就确定下来了。
合水层的边界分为自然边界与人为边界两类
当研究的合水层系统(包括弱透水的含水层)与非含水层相接触时,其界 面(线)称为自然边界；
建立数值模型
网格剖分：根据确定的数值方法和软件，对研究区进行剖分。 对于平面二维流问题，一般将研究区剖分成矩形或三角形网格； 对于剖面二维流问题，一般也是将剖面区域剖分成矩形或三角形； 对于三维流问题，一般先在垂向上分成若干层，而在每层剖分成矩 形或三角形。
边界条件： 初始条件： 含水层参数：渗透系数、储水系数、给水度、孔隙度 源汇项：降雨入渗、河流补给、蒸发排泄、地表水体、沟渠渗漏、
上述数学模型是一个偏微分方程定解问题，通常只能用数值方法 求解，常用的数值方法有：有限差分法和有限单元法。
目前有一些软件可以直接用于求解地下水流动问题，如果不是自 己编写程序，可以选择合适的软件，建立数值模型。
如果用软件，则需对软件功能作简要介绍，论述软件的适用性。
3
第七章 数值模型一般步骤及对 勘查资料的要求
模型设计者应出具有较高理论水平和丰富经验的水文地质工作者担 任。设计者应精细地分析有关资料, 以获得较符合实际条件的分区 图。
18
三、数值模型设计中一些特殊问题
（一）抽水试验设计 （二）抽水试验数值模拟设计 （一）含水层剖分注意事项
（一）抽水试验设计
19
三、数值模型设计中一些特殊问题
（一）含水层剖分注意事项
12
二、数值模型设计及对资料和水文地质勘探的要求
（四）边界条件的确定
边界条件是与计算区的范围同时确定的
考虑计算区的范围时, 必须同时确定边界条件的性质； 反之, 边界的位置一旦确定, 计算区的范围自然也就确定下来了。
合水层的边界分为自然边界与人为边界两类
当研究的合水层系统(包括弱透水的含水层)与非含水层相接触时,其界 面(线)称为自然边界；
地下水数值模拟03
差分公式对比
名称
前差
公式
df f ( x x) f ( x) dx x
截断误差
O(x) O(x)
一 阶 导 数
后差 中心差
df f ( x) f ( x x) dx x
df f ( x x) f ( x x) dx 2x
O(x )
2
二阶导数
该问题属于一维承压非稳定流的定解问题。求解该问题，需要对空间 和时间进行离散，形成的差分网格称为时空网格。
网格剖分
• • •
以等距剖分为例 将研究区域[0，L]用直线等分为n份，把时间段[0,T]用直线等分成m份 以 H ik 表示结点（i,k）处的水头
2H * H T 2 x t H ( x, t ) H 0 ( x ) t 0 H ( x, t ) (t ) 0 x 0 H ( x, t ) x L L (t )
开始
输入初始参数 t=0 对i=0,1,…,n循环执行
hi0 H0 ix
t=t+Δt 对i=1,…,n-1循环执行
hik 1 hik1 (1 2 )hik hik1
• 用差分方程的解作为微分方程的近似解。
一、有限差分法的基本思想
• 2、求解步骤
选取网格； 对微分方程及定解条件选择差分近似，列出
差分格式；
求解差分格式； 讨论差分格式解对于微分方程解的收敛性及 误差估计。
差分网格
结点 格点
y
x
Δx，Δy——空间步长 Δt——时间步长
一、有限差分法的基本思想
显式差分方程的求解
差分方程与离散后的定界条件构成了数学模型 的显式差分方程问题。其求解步骤如下： 由初始条件给出t0时刻各结点的水头值h00， h10，…hl0；再根据差分方程，在k=0时， 分别取i=1，i=2,…i=n-1,便可求得t1时刻各 内结点的水头值h11，…hn-10。 由以上计算的h11， h21 …hn-11值及由边界 条件计算的h01和hl1,再次利用差分公式 （取k=1， i=1，i=2,…i=n-1，便可计算得 t2时刻各结点的水头值。如此重复，便可 计算出t3 ，t4 ，…各时刻的水头分布值。
地下水数值模拟
一、模型概化-内部结构
1、含水介质 2、含水层空间分布
确定含水层类型，查明含水层在空间的分布形状。对承压 水，可用顶底板等值线图或含水层等厚度图来表示；对潜水， 则可用底板标高等值线图来表示；
一、模型概化-内部结构
1、含水介质 2、含水层空间分布 4、地下水运动状态
（1）层流、紊流 一般情况下，在松散含水层及发育较均匀的裂隙、岩溶含水层中的地下
三维有限差分模型与MODFLOW
**使用VMODFLOW建立模型 1. 准备数据资料（划分含水层，顶底板高度，
渗透参数，存贮参数，初始地下水位） 2. 划分平面网格，输入顶、底板高度 3. 输入模拟层参数 4. 输入边界属性和模块参数 5. 是稳定流还是非稳定流 6. 非稳定流的初始水头分布 7. 运行模型，查看结果
一模型概化内部结构11含水介质含水介质2含水层空间分布含水层空间分布3地下水运动状态地下水运动状态4水文地质参数水文地质参数1时间概化时间概化水文地质参数是慢时变的在一定时期和外部条件下可以近水文地质参数是慢时变的在一定时期和外部条件下可以近似地看作恒定不变建立概念模型时将参数概化为随时间不变似地看作恒定不变建立概念模型时将参数概化为随时间不变一模型概化内部结构11含水介质含水介质22含水层空间分布含水层空间分布33地下水运动状态地下水运动状态44水文地质参数水文地质参数11时间概化时间概化2空间概化空间概化查明含水层的导水性储水性及主渗透方向的变化规律
模拟步骤
建立概念模型 选择数学模型 将数学模型进行
数值化 模型校正(识别) 模型验证 预测
一、模型概化-边界条件
1、研究区边界
研究区应尽可能以自然边界为计算边界，最好是以 完整的水文地质单元作为计算区。
第六章(3) 地下水数值模拟模型的应用实例xiugai
图5-51 距垃圾场中心约150m（cell Id：4227） 溶解氧浓度变化曲线
3、硫酸根在地下水中的浓度分布规律
图5-60 垃圾场中心（cell Id：4867）硫酸根浓度变化曲线
图5-63 距垃圾场中心约200m污染羽边缘（cell Id：4147） 硫酸根浓度变化曲线
在所有点中硫酸根都随季节有轻微的 波动，这可能是由于雨季降雨对地下 水中随酸根的稀释作用。 在填埋场中心，由于硫酸根是补给源， 硫酸根的浓度是持续上升的，并在 6000d左右时基本稳定。 随着距填埋场的距离增加，硫酸根上 升的速率逐渐降低。 在污染羽边缘（约200m）处，硫酸根 先是降低然后增加，这是因为一开始 该区域固有的硫酸根先与BTEX反应， 而垃圾渗滤液中的硫酸根还没有补给 到这一区域，所以硫酸根离子先降低， 随后硫酸根逐渐补给，故其浓度又开 始逐渐增加。 在污染羽范围之外，硫酸根离子则是 逐渐降低。
a：概念模型
场区为三面环山地形，与单斜地层基本构成一独立的水文地质单元。 填埋场地下方分布有新生界第四系及中生界侏罗系、古生界二叠系地层。 本区第四系地层分布不甚规律，缺少全新统。上更新统下段洪积层（莲花组 PLQ3）以黄色粘土、亚粘土组成，缺少砂砾层位。平均厚度约3米左右；中更 新统洪积层（之江组PLQ2）以棕红色粘土、亚粘土组成，结构紧密，层厚一 般为7米左右。 第四系下覆有埋藏型裂隙岩溶水，含水层由古生界二叠系下统栖霞组 （P3q）、石炭系上统船山组（C3c）、中统黄龙组（Czh）灰岩、白云质灰岩 组成。
（2）厌氧状态下在反硝化细菌作用下，NO3-作为电子接受体，接 受电子氧化BTEX。 （3）厌氧状态下在铁还原菌的作用下，Fe(Ⅲ)还原为Fe(Ⅱ)接受 电子，氧化BTEX。 （4）厌氧状态下在硫酸根还原菌的作用下，S6+还原为S2-接受电 子，氧化BTEX。
地下水数值模拟模型简介
Problems about GW Development
Functions of GW
An Important Water Supply for Human An Important Environmer Shortage Ecologic problems Saline soil Sea Water Intrusion Land Subsidence Groundwater Contamination ……
基 本 思 路
运用地下水模型计算各组参数情况下地下水位（hsij)
根据目标函数评价各组参数值的优劣
输出结
Y
两次迭代的差足够小？
N
考虑到变量上下限等约束，根据一定的准则， 生成新的N组参数
确定决策变量、目 标函数和约束条件
研究工作区 水文地质条件
产生初始种群N 建立水文地质 概念模型
编
码
求解水均衡约束： W(X,)=0 得到地下水系统状 态
Calibration of the Model
目标：使模拟期内所有观测点的水位误差的平方和最小。 变量： Sk、Kk (k=1,…,L) 约束：
各参数分区参数的上下限约束 连续性方程约束
m n
2 obj . z ( hs hr ) ij ij j 1 i 1
s .t.
Groundwater Model_Math Model
Mass Balance
Q Q Q input output
Darcy’s Law
dh q K dl
Boundary Conditions
First Type Second Type Mixed Type
地下水渗流基本方程及数学模型总结
V Qt vt m V vt
§5 描述地下水运动的数学模型及解算方法
X方向流入流出差
( v x ) |( x , y , z ,t ) yzt ( v x ) |( x x , y , z ,t ) yzt
y方向流入流出差
( v y ) |( x , y , z ,t ) xzt ( v y ) |( x , y y , z ,t ) xzt
流入－ 流出＝V
研究对象可以是大区域的，也可以是微分单元体 大区域的水均衡计算经常用于区域的水资源评价 地下水渗流的连续性方程就是质量守恒方程，也称为水
均衡方程。现在基于微分单元体做水均衡分析，推导渗 流连续性方程。
（一）渗流连续性方程的推导
在渗流场中任意取一点P(x, y, z)，以P为 中心沿直角坐标轴取一微小的六面体，体 积为 ，称为特征单元体，设单元 体无限小，但保证单元体中地下水穿过介 质骨架和空隙。 假设：水是可压缩的，固体颗粒不能被压 缩，多孔介质骨架在垂直方向（Z）可压 缩，x、y为常量。因此，只有水的密 度、孔隙度n和单元体高度z三个量随压 力而变化。
( v y ) xyzt y
y 0
z方向流入流出差 单元体内地下 水质量变化量 地下水运动的 连续性方程：
( v z ) xyzt z
z 0
( v x ) ( v y ) ( v z ) ( nz ) x y z xy y K xx ) x x x
H ( K yy ) y y ( v z ) H ( K zz ) z z z
( v y ) y
§5 描述地下水运动的数学模型及解算方法
( nz ) 第二步：化简方程右端项： t
地下水数值模拟02_地下水运动的数学模型
2
H 0
n 2
——隔水边界
第三类边界条件 H aH b n
例：弱透水边界
K H Hn H 0 n m1 / K1
溶质运移问题的边界条件
第一类边界条件
c(x,
y, z,t) 1
c1(x,
y, z,t)
——给定浓度边界
第二类边界条件 c
Di, j x j ni 2 f2 (xi , t)
u(x, y, z,t) t0 0(x, y, z)
• 2、边界条件
第一类边界条件 u(x, y, z,t) 1 1(x, y, z,t)
第二类边界条件
u n
2
1(x, y, z,t)
第三类边界条件
u
u n
3
3x,
y, z,t
水流问题的边界条件
Reynolds数小于1~10
• 有些情况下，用液体压强表示更为方便
– 例如：油水两相流动
vx
K
H x
vy
K
H y
vz
K
H z
K g k
H z p
g
k p
vx
x
v y
k
p y
vz
k
K ( d
)
dhc
C
t
x
K( )
x
y
K
(
)
y
z
K (
地下水运动的数学模型
第四章 地下水运动的数值模型解析解虽然具有精确可靠的特点，但采用解析解反映自然状态和复杂人类活动干扰下的地下水运动是相当困难的。
因此，当含水层的条件严重偏离现有解析模型的简化假设时，人们通过数值模型来获得近似的地下水流场及演变趋势。
第一节 地下水流数值方法概述地下水流的数学模型采用偏微分方程描述地下水流的时间和空间连续状态，而数值模型则是采用离散（非连续）时空模型中水头的分布与演变对数学模型进行近似描述。
从精确数学模型到近似数值模型的转化，虽然会损失一些精度，但使复杂地下水流问题的分析得以通过机械计算实现，而且误差也是可控的。
把偏微分方程求解的数值方法引入到地下水流问题的求解始于20世纪70年代，主要方法包括有限差分法、有限元法和边界元法，此后又发展了有限分析法、多重网格法和无网格法等。
这些方法的共同特点是将模型空间及边界离散为由一系列的节点以及联系这些节点的单元（无网格法除外），含水层的水头在这些节点上定义，从而实现了水头分布空间连续函数向离散变量的转化，表示为2121211122111221202()02()02()002(0)k k k k k k k k k k k k k k k k k k k k k k k kkk f f f f a b c e x L x x t t t t f x f f f f a b c e x L x x t t f f f f a b c x L e x xd f dfe ef a b f c x L dx dx t t f x u---------∂∂-++=<<∂∂∆∆=-∂∂-++=<<∂∂∆∆∂∂=+++<<∂∂+-++=<<∆∆==，，，，{}(,,);1,2,3,,p H x y z H p M ⇒=⋅⋅⋅ （4.1.1）式中；H 为含水层的水头；x 、y 、z 为空间坐标；p 为数值模型的节点；M 为节点的数目。
地下水流数值模型设计与应用
地下水流数值模型设计与应用地下水流数值模型是指利用数值方法对地下水流动过程进行模拟和预测的模型。
在地下水资源管理、地下水开采与补给、地下水污染治理等领域，地下水流数值模型具有重要的应用价值。
本文将介绍地下水流数值模型的设计原理、建模步骤和应用案例。
1. 地下水流数值模型的设计原理地下水流数值模型是基于地下水流动的物理原理和方程建立的数学模型。
其设计原理主要包括以下几个方面：1.1 地下水流动方程地下水流动方程是描述地下水在地下各层介质中运动规律的基本方程。
常用的地下水流动方程有Richards方程、Darcy方程等。
1.2 数值方法地下水流数值模型的建立离不开数值方法的支持。
常用的数值方法包括有限差分法、有限元法、边界元法等。
1.3 参数估计与校正地下水流数值模型需要估计和校正地下水系统的参数，包括渗透系数、压力头等，以保证模型的准确性和可靠性。
2. 地下水流数值模型的建模步骤地下水流数值模型的建模过程一般包括以下几个步骤：2.1 地下水系统的规划与数据收集首先需要对地下水系统进行规划和确定研究范围，然后收集相关的地下水系统数据，包括地质条件、水位、水质等信息。
2.2 模型网格的建立根据地下水系统的特性和数据建立模型的网格结构，划分各个单元，并确定边界条件。
2.3 模型参数的估计与校正根据实测数据和地质资料，对模型中的各项参数进行估计和校正，以提高模型的准确性。
2.4 模型求解与验证利用数值方法对建立的地下水流数值模型进行求解，并与实际观测数据进行对比验证，调整模型参数以提高模型的拟合程度。
3. 地下水流数值模型的应用案例3.1 地下水资源管理通过建立地下水流数值模型，可以预测地下水资源的分布和变化规律，为地下水资源的合理开发利用提供科学依据。
3.2 地下水开采与补给地下水流数值模型可以模拟地下水的开采和补给过程，帮助调整地下水的开采量和位置，保护地下水资源的可持续利用。
3.3 地下水污染治理地下水流数值模型可以模拟地下水中污染物的传输规律，对地下水污染的源头和扩散路径进行分析，为地下水污染治理提供技术支持。
地下水资源管理模型及模型
决非取之不尽，用之不竭。
（2）地下水既是一种地质营力，又是一种重要的环境
要素。
（3）人类通过开发地下水资源参与社会经济生活，不 仅仅是满足生产和生活需要，从社会全局利益出发，还存在 着一个经济效益的问题。也就是说，对于有限的地下水资源，
全社会应如何分配，如何合理利用，才能取得最佳的经济效
益。
5.1.3系统分析的特点及其在地下水资源管理中的应用
近几十年来，世界各国学者都十分重视地下水资源管理模型的建立，
并已建立了许多地下水资源管理模型，许涓铭教授将其归纳为以下十 种：
①区域国民经济的产业结构规划与水资源管理
②地下水政策评估与分配 ③地下水—地表水系统优化调度与开发 ④含水层—河流系统优化开发与运行 ⑤城市发展的地下水规划与管理 ⑥供水水源地优化设计与管理 ⑦灌溉地下水和含盐量控制的规划与管理 ⑧ 排水的最优化设计与管理
（1）系统性
地下水系统的概念，可归纳为：
①是由若干个具有一定独立性而又相互联系、互有影响的
不同等级的子系统或次子系统所组成；
②是水文系统的一个组成部分，与降水—径流子系统和地
表水子系统存在密切联系，互相转化；地下水系统的演变，
很大程度上受这些子系统的输入与输出子系统的控制；
③每个地下水系统都具有各自的特征和演变规律，包括各自
五、为专门目的服务的地下水管理内容
1.为城市生活和工业用水服务的地下水管理
城市供水存在的问题：
① 集中供水造成水源地的地下水位下降过快。 ② 水文地球化学环境的改变造成地下水矿化度、硬度的 升高。
③ 过量开采引起地面沉降。
城市地下水管理的主要内容有：
（2）地下水既是一种地质营力，又是一种重要的环境
要素。
（3）人类通过开发地下水资源参与社会经济生活，不 仅仅是满足生产和生活需要，从社会全局利益出发，还存在 着一个经济效益的问题。也就是说，对于有限的地下水资源，
全社会应如何分配，如何合理利用，才能取得最佳的经济效
益。
5.1.3系统分析的特点及其在地下水资源管理中的应用
近几十年来，世界各国学者都十分重视地下水资源管理模型的建立，
并已建立了许多地下水资源管理模型，许涓铭教授将其归纳为以下十 种：
①区域国民经济的产业结构规划与水资源管理
②地下水政策评估与分配 ③地下水—地表水系统优化调度与开发 ④含水层—河流系统优化开发与运行 ⑤城市发展的地下水规划与管理 ⑥供水水源地优化设计与管理 ⑦灌溉地下水和含盐量控制的规划与管理 ⑧ 排水的最优化设计与管理
（1）系统性
地下水系统的概念，可归纳为：
①是由若干个具有一定独立性而又相互联系、互有影响的
不同等级的子系统或次子系统所组成；
②是水文系统的一个组成部分，与降水—径流子系统和地
表水子系统存在密切联系，互相转化；地下水系统的演变，
很大程度上受这些子系统的输入与输出子系统的控制；
③每个地下水系统都具有各自的特征和演变规律，包括各自
五、为专门目的服务的地下水管理内容
1.为城市生活和工业用水服务的地下水管理
城市供水存在的问题：
① 集中供水造成水源地的地下水位下降过快。 ② 水文地球化学环境的改变造成地下水矿化度、硬度的 升高。
③ 过量开采引起地面沉降。
城市地下水管理的主要内容有：
地下水数值模拟 第三讲 地下水数值模拟原理及建模方法和步骤_xiugai
N Hn i, j i, j1
(Wi, j
Ei, j
Si, j
Ni, j
)
H
n i,
j
i,
j
Ai, j
Ci,
j
H
n i,
j
]/
Ci, j
不等距矩形网格有限差分 y
(i,j+1)
法：以二维承压水流为例（续3）
(i-1,j)
(i,j) (i+1,j)
节点（i,j）的均衡区
如果 H 取tn+1时刻的值
Hn i, j1
H
n i,
j
y j y j1
xi1
2
xi1
(KM
)i, j1/ 2
Hn i, j1
H
n i,
j
y j1 y j
xi1
2
xi1
i,
j
Ai, j
ei, j
H n1 i, j
H
n i,
j
tn1 tn
Ai, j
Aij:节点(i,j)均衡区的面积
e 表示：当水头下降一个单位时，从单位底面积饱水空隙介质中 释放的水量（体积）
xi1
2
xi1
i,
j
Ai,
j
ei, j
H n1 i, j
H
n i,
j
tn1 tn
Ai, j
H 可取tn+1或 tn时 刻的值
Aij:节点(i,j)均衡区的面积
e 表示：当水头下降一个单位时，从单位底面积饱水空隙介质中 释放的水量（体积）
不等距矩形网格有限差分 y
(i,j+1)
法：以二维承压水流为例（续1）
地下水位数学模型构建及其在水资源管理中的应用
地下水位数学模型构建及其在水资源管理中的应用地下水是地球上最重要的水资源之一。
然而，由于人类活动和自然灾害的影响，地下水资源面临着严重的威胁。
因此，建立数学模型是水资源管理中不可忽视的重要工具。
本文将探讨地下水位数学模型的构建及其在水资源管理中的应用。
一、地下水位数学模型的构建地下水位数学模型是通过将实际的地下水位变化与可能影响地下水位变化的各种因素，例如降雨和蒸发量等进行比较和分析，来预测未来地下水变化的数学模型。
利用该模型，我们可以准确地估计地下水资源的可持续性，以及接下来的输送和使用过程中可能会遇到的问题。
（一）地下水位数学模型的类型地下水位数学模型通常根据其研究所要建立的领域和具体目的而划分。
常见的类型如下：1.分布式模型该模型被用于描述某一特定区域的地下水位变化。
由于其强大的数据分析功能，它可以准确地研究地下水位的变化规律和可能出现的问题，从而使我们能够有针对性地规划地下水资源的开发。
2.系统动力学模型该模型适用于考虑多种不同因素对地下水位变化的影响，例如人类活动和自然灾害的作用。
该模型利用现代计算机技术，模拟各种可能的因素，并预测它们在特定区域内的作用与影响。
3.时间序列模型该模型适用于以时间为基础的研究。
它将地下水位的变化作为时间序列进行建模，使我们可以推断出未来的变化趋势或可能出现的问题。
（二）地下水位数学模型的参数确定地下水位数学模型的准确性高度依赖于其参数的选择和确定。
这些参数包括地质条件、地下水位历史数据、降雨量、蒸发量等。
为了确保模型的精度，我们需要采取以下步骤：1.确定地下水位的基本特征，如地下水位的年际变化和季节变化等。
2.分析地下水位与可能影响它的各种影响因素之间的关系，可以是人类活动或自然界。
3.建立数学模型并验证准确性。
验证模型的准确性方法可以采用模型的误差分析方法、交叉验证等方法评估预测结果，使模型更加可靠。
二、地下水位数学模型在水资源管理中的应用地下水位数学模型在水资源管理中具有重要的应用价值。
地下水资源管理模型及模型
1.地下水位的控制与优化 地下水位的变化主要受地下水的补给量和开采量控制，所以 控制地下水位最有效的办法是合理分配开采量。
首先应该做的一点就是查明该地区的水文地质条件，应用合
理的方法评价地下水的天然补给量，这样在分配开采量时心中有
数。
高：盐渍化，增大开采量
地下水位控制
低：区域低，减采 局部低，合理分配开采量
一般为：
目标函数
约束条件
g i x bi
i 1,2, m
min
max F X
5.2.2地下水资源管理模型的分类 在进行地下水资源规划与管理时，首先要根据研究区的条件
及规划管理的目标要求，选定所需建立的管理模型的类型。
地下水资源管理模型分类如下：
1）按系统管理目标划分 ⑪多目标系统管理模型—主要用于地下水系统综合管理，模
1955年哈佛大学开始修订水资源大纲，其目的是为规划 和设计复杂的现代水资源系统工程研究一种方法。
50—60年代的水资源管理，大多是将地下水系统作为集
中参数系统来处理。
1962年哈佛大学发表 “水资源系统规划”一书。
20世纪70年代，将描述地下水状态的数值模拟模型与
优化技术结合起来，形成分布参数地下水管理模型。
的水动力系统、水化学系统等；
④含水层系统与地下水系统代表两种不同的概念，前者具有
固定边界，而后者的边界是自由可变的；
⑤地下水系统的时空分布与演变规律，既受天然条件的控制，
又受社会环境、特别是人类活动的影响而发生变化。
由此可知：地下水系统是一个错综复杂，包括各种天然因素、
人为因素所控制的，具有不同等级的互有联系又互相影响， 在时空分布上具有四维性质和特征，不断运动演化的若干独 立单元的统一体。
第三讲 地下水数值模拟原理及建模方法和步骤
的迁移机理及数学模型和求解方法
地下水数值模拟
绪论
地下水数值方法在水文地质学中的位置
地下水动力学主要内容
连续性原理、达西定律、水均衡原理 地下水流基本方程 几类特殊水文地质问题数学模型及解析解
地下水向沟渠河中的流动 园岛模型 泰斯模型 有越流的不稳定井流（Hantush and Jacob） 无越流的潜水含水层不稳定井流（ Neuman ）
为f(x)在x0处的二阶中心差商，O(x)2
为截断误差。
方法一
(2)有限差分方程建立(续)
对于偏导数（偏微商），类似可以得到相应的差商：
H (x0,t0 ) H (x0,t0 t) H (x0,t0 )
t
t
H (x0 ,t0 ) H (x0 x,t0 ) H (x0 ,t0 )
x
x
2H (x0 , t0 ) x 2
➢ 方法二．在渗流区剖分的基础上，直 接由达西定律和水均衡原理，建立各 个均衡区的水均衡方程，从而得到差 分方程。适用于矩形网格、三角形网 格。
矩形网格 多边形网格
方法一：差商代替微商
1、网格划分的基本类型
（1）先划格线，格点位
均
于网格中心
衡
网
格
（2）先规定格点位置，
再垂直平分两相邻结点的连
节
线作格线，形成的网格即为 水均衡区
绪论
数值解与解析解
数值方法是地下水动力学的完善和补充或延续 数值解的特点：
只是求出研究区某些空间点和某些时间点处的水头值 适用于所有的问题 具备水文地质基础和线性代数知识 已有数值模拟专门软件（或自己编程） 需要有高性能计算机 对实际问题的刻画比较精确
因此，其应用非常广泛
绪论
地下水数值模拟
绪论
地下水数值方法在水文地质学中的位置
地下水动力学主要内容
连续性原理、达西定律、水均衡原理 地下水流基本方程 几类特殊水文地质问题数学模型及解析解
地下水向沟渠河中的流动 园岛模型 泰斯模型 有越流的不稳定井流（Hantush and Jacob） 无越流的潜水含水层不稳定井流（ Neuman ）
为f(x)在x0处的二阶中心差商，O(x)2
为截断误差。
方法一
(2)有限差分方程建立(续)
对于偏导数（偏微商），类似可以得到相应的差商：
H (x0,t0 ) H (x0,t0 t) H (x0,t0 )
t
t
H (x0 ,t0 ) H (x0 x,t0 ) H (x0 ,t0 )
x
x
2H (x0 , t0 ) x 2
➢ 方法二．在渗流区剖分的基础上，直 接由达西定律和水均衡原理，建立各 个均衡区的水均衡方程，从而得到差 分方程。适用于矩形网格、三角形网 格。
矩形网格 多边形网格
方法一：差商代替微商
1、网格划分的基本类型
（1）先划格线，格点位
均
于网格中心
衡
网
格
（2）先规定格点位置，
再垂直平分两相邻结点的连
节
线作格线，形成的网格即为 水均衡区
绪论
数值解与解析解
数值方法是地下水动力学的完善和补充或延续 数值解的特点：
只是求出研究区某些空间点和某些时间点处的水头值 适用于所有的问题 具备水文地质基础和线性代数知识 已有数值模拟专门软件（或自己编程） 需要有高性能计算机 对实际问题的刻画比较精确
因此，其应用非常广泛
绪论
1-7描述地下水运动的数学模型及其解法
2020/4/15
12
三、建立数学模型的实例
例１ 潜水非稳定流方程
(Kh H ) (Kh H ) W H
x
x y
y
t
H c1 H(t)
H c2 z(t)(位置水头)
H c3 z(t（) 渗出面）
H c4 hw (t（) 定水头边界）
H n
c5 （ 0 隔水边界）
H (x, y,t) t0 (x, y)
➢区域的抽水井、注水井或疏干巷道也可作为 给定水头边界处理；
➢无限边界 H(x, y,t)
x2 y2
H
亦为第一类边界；
0
➢潜水面任一点的水位已知时，抽水井井壁水
位为一类边界。
2020/4/15
9
（2）第二类边界条件
当已知渗流区某部分边界上的流量分布时，称这 部分边界为第二类边界或给定流量边界。相应的边界 条件表示为：
（1）第一类边界条件
若在渗流区的某部分边界上各点在每一时刻的水头是已知
的，则称这部分边界为第一类边界或给定水头边界，常表
示为：
H(x, y, z,t) S1
1(x, y, z,t),
H (x, y,t) 1 2 (x, y,t),
（x, y, z) S1
（x, y) 1
分别表示在三维和 二维条件下边界上 2的020点/4/1在5 t时刻的水头
H
K n
s2
q (x, y, z,t), 1
（x, y, z) S2
或
T H n
2
q (x, 2
y,t)， （x, y) 2
式中：n为边界 S2 或 2的外法线方向； q1和q2为已知函数，分别表示S2 上单位面积和 2上单位
地下水数值模拟
j
i
~ fij
j
j1
Gi
d
,
gij
j1 Gi d , j n
g~ij
j1 Gi d , j n
Ei
D
W T
Gi dxdy
iHi
n j 1
1 j1 j
j1gij
g~ij
Hj
g~ij j gij
H
j 1
j1
fij
~ fij
H n
j
~ fij
j
fij
iHi
n j 1
M j M j1
H
Gi n
Gi
H n
ds
W DT
Gi dxdy
n j 1
j1 j
j1 j1
j
H
j j j1 j
H
Gi
j1
n
j1 j1
j
H n
j j j1 j
H n
j 1
Gi
d
D
W T
Gi dxdy
令：fij
G d , j1
地下水数值模拟
一、基本原理
• 基本思想
——将微分方程得基本解化为边界积分方程, 将边界剖分为有限个单元,在离散得区域边 界上将边界积分方程化为代数方程求解。
• 边界元 ——区域内满足控制方程,边界上近似满足边界条件
• 有限元、有限差 ——区内近似满足控制方程,边界上满足边界条件
一、基本原理
• 特点
u x
v x
u y
v y
dxdy
v
2u x 2
2u y 2
dxdy
v
u n
ds
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
b 有限单元法：基本思想是用有限个单元的集合代替渗透区，用比较 简单的函数代替每个单元上的实际水头函数，然后通过某种变换处理， 把偏微分方程的定解问题转化为求解线性方程组的问题。
c 边界元法：边界元的剖分仅仅在边界上进行，然后利用边界积分方 程求出边界上的未知量，并进而计算区域内的未知量。
4.2 地下水流动数学模型求解的有限单元法
目前在实际工作中，研究途径主要有： （1）北方地区，由于岩溶发育程度相对较弱，多数情况，地下水的运动基本上 服从达西定律，所以常用孔隙水的数学模型或双重介质模型近似地进行地下水研究。 （2）南方地区，可考虑同时存在裂隙流和管道流两种水流情况，裂隙流服从达 西定律，管道流则不服从，裂隙流区与管道流区的分界线上存在着水量交换，这样 可推导出相应的地下水流模型。
地下水系统是一个受多因素控制的系统，这些因素中有些是随机的、 不确定的。严格地说，地下水流动数学模型应该用随机的微分方程（或 积分方程）来表示。
随着地下水赋存介质的特征不同，地下水的流动特点也就不相同，相 应的地下水流数学模型的建立和研究途径也就存在着差异。
4.1 地下水流动数学模型及其求解方法
1 孔隙水 在孔隙介质中，地下水的运动服从达西定律。
第一节 地下水系统的控制 第二节 地下水系统的管理 第三节 目标函数、约束条件、决策变量、状态变量
第五章 地下水系统的控制与管理
第一节 地下水系统的控制
地下水系统由自然系统转为自然-人工复合系统，其行为、功能及活 动方式受到人为活动的影响，使其具有目的性并具有适应目标价值多样 化的特点。为使地下水的开采利用不至于造成系统的失稳，又使其活动 和状态朝着我们所希望的方向发展，就需对施加给地下水系统的（人为) 作用进行选择。这种选择过程就是控制。
一 基本知识 1 泛函 2 变分原理 （1）描述地下水运动的偏微分方程的定解问题等价于求某一泛函的 极小值问题；（2）求二次泛函的极小值问题与求解某一线性代数方 程的解是等价的。 3 剖分与插值
4.2 地下水流动数学模型求解的有限单元法
二 稳定流问题的解法 三 非稳定流问题的解法
第五章 地下水系统的控制与管理
要沿裂隙进行，而地下水的储存主要在孔隙中，孔隙和裂隙中各有不同 的水头，它们之间存在着水量交换，由此建立的数学模型称为双重介质 模型。
③把裂隙岩石看作是一种按一定几何规律分布的裂隙系统所割裂的 介质，此时，研究裂隙的参数和建立裂隙介质的水力学（层流或紊流模 型）就应该考虑单个裂隙的形状、规模及排列。
4.1 地下水流动数学模型及其求解方法
二、地下水流数学模型的内容及其求解方法 1 一般包含两方面的内容： （1）有一个或一组能描述地下水运动规律的偏微分方程。一般，可把它 理解为是水量（或质量）平衡方程的一种微分表达形式；【注意对水文地质 条件的分析】 （2）有相应的定解条件。定解条件包括初始条件和边界条件。当模拟稳 定流时，只需要边界条件；当模拟非稳定流时，则要同时给出初始条件和边 界条件。 初始条件是指某一选定的初始时刻渗流区内地下水的水头分布状况；边 界条件是指渗流区的边界位置及边界上地下水应满足的条件，一般包括有已 知水头边界、已知流量边界和混合边界等三类。
第四章 地下水系统分布参数数学模型
4.1 地下水流动数学模型及其求解方法 4.2 地 地下水流动数学模型及其求解方法
一 地下水流数学模型 分布参数地下水流数学模型可分为不同的类型。按赋存介质的不同，
可分为孔隙水、裂隙水和岩溶水等模型；按流动要素的特点，可分为稳 定流、非稳定流、承压流、非承压流、饱和流、非饱和流、单相流、多 相流、一维流、二维流和三维流等模型。
4.1 地下水流动数学模型及其求解方法
2 地下水流数学模型的求解方法 （1）解析法是利用经过严格的推理和数学推演获得的各种计算公
式进行计算的方法，裘布依公式和泰斯公式即属此类。其特点是：在 一定的假设条件下所进行的数学推导是非常严密的，所获得的解是精 确的。但水文地质条件非常复杂，用此法求解，误差非常大。
2 裂隙水 （1）裂隙介质的特点 ①其渗透性具有较强的非均质性和各向异性，在一个计算区内往往可以划出许
多非均值区，每个区之间的导水系数或渗透系数有时可以相差数百倍。 ②渗透性往往随深度发生变化，在不同地点或不同的岩层，渗透系数随深度的
变化规律也不相同。因此，在厚层裂隙介质中确定含水层的厚度比较困难。 ③储水能力比孔隙介质小。 ④裂隙介质的储水性主要取决于微裂隙和岩层的孔隙（如粒间孔隙、气孔等），
有限单元法可分为变分有限元法和迦辽金有限元法。变分有限元法的 基本思想是：利用变分原理，把地下水流偏微分方程的定解问题转化 为某一泛函I[H]的极值问题，再通过剖分和插值，可将该泛函I[H]转 化为一个普通的二次函数。这个二次函数的极值问题等价于一个线性 代数方程的求解问题，求解该方程组，即可得到地下水流问题的解。
（2）数值法是一种只能求出计算区内有限个点在某些时刻近似解 的计算方法。也就是说，数值法不能求出解的解析表达式，只能求出 解（未知量）在若干个离散点的近似值。
数值法在处理渗透区的非均质性和不规则性等方面具有很大的灵 活性。
4.1 地下水流动数学模型及其求解方法
a 有限差分法：原理就是在空间域时间离散的基础上用差商代替微商， 最后把地下水运动的偏微分方程的求解问题转化为差分方程的求解问 题；
4.1 地下水流动数学模型及其求解方法
3 岩溶水
岩溶水与裂隙水在许多方面有共同之处，因而裂隙水的某些特点对岩溶水也同 样使用。但是，岩溶介质比裂隙介质具有更强的非均质性，使建立岩溶水流数学模 型更困难。
北方与南方岩溶水在分布、径流、补给、动态等方面存在着差异。南方岩溶地区， 地下水流动以管道流为主；北方地区，岩溶发育相对较弱，多数裂隙-岩溶水，含水 层中的地下水往往具有统一的地下水面。
而导水性则主要取决于一些较大的裂隙。 ⑤由于压力对裂隙张开的程度有明显的影响，故压力变化与岩层的透水性之间
存在着密切的关系。
4.1 地下水流动数学模型及其求解方法
（2）裂隙介质渗透性研究的途径 ①把裂隙介质视为连续介质，认为裂隙水的运动基本上服从达西定
律，可以用与孔隙水相同的数学模型来计算裂隙水流。 ②仍然认为裂隙水的运动符合达西定律，但地下水的传导和运动主
c 边界元法：边界元的剖分仅仅在边界上进行，然后利用边界积分方 程求出边界上的未知量，并进而计算区域内的未知量。
4.2 地下水流动数学模型求解的有限单元法
目前在实际工作中，研究途径主要有： （1）北方地区，由于岩溶发育程度相对较弱，多数情况，地下水的运动基本上 服从达西定律，所以常用孔隙水的数学模型或双重介质模型近似地进行地下水研究。 （2）南方地区，可考虑同时存在裂隙流和管道流两种水流情况，裂隙流服从达 西定律，管道流则不服从，裂隙流区与管道流区的分界线上存在着水量交换，这样 可推导出相应的地下水流模型。
地下水系统是一个受多因素控制的系统，这些因素中有些是随机的、 不确定的。严格地说，地下水流动数学模型应该用随机的微分方程（或 积分方程）来表示。
随着地下水赋存介质的特征不同，地下水的流动特点也就不相同，相 应的地下水流数学模型的建立和研究途径也就存在着差异。
4.1 地下水流动数学模型及其求解方法
1 孔隙水 在孔隙介质中，地下水的运动服从达西定律。
第一节 地下水系统的控制 第二节 地下水系统的管理 第三节 目标函数、约束条件、决策变量、状态变量
第五章 地下水系统的控制与管理
第一节 地下水系统的控制
地下水系统由自然系统转为自然-人工复合系统，其行为、功能及活 动方式受到人为活动的影响，使其具有目的性并具有适应目标价值多样 化的特点。为使地下水的开采利用不至于造成系统的失稳，又使其活动 和状态朝着我们所希望的方向发展，就需对施加给地下水系统的（人为) 作用进行选择。这种选择过程就是控制。
一 基本知识 1 泛函 2 变分原理 （1）描述地下水运动的偏微分方程的定解问题等价于求某一泛函的 极小值问题；（2）求二次泛函的极小值问题与求解某一线性代数方 程的解是等价的。 3 剖分与插值
4.2 地下水流动数学模型求解的有限单元法
二 稳定流问题的解法 三 非稳定流问题的解法
第五章 地下水系统的控制与管理
要沿裂隙进行，而地下水的储存主要在孔隙中，孔隙和裂隙中各有不同 的水头，它们之间存在着水量交换，由此建立的数学模型称为双重介质 模型。
③把裂隙岩石看作是一种按一定几何规律分布的裂隙系统所割裂的 介质，此时，研究裂隙的参数和建立裂隙介质的水力学（层流或紊流模 型）就应该考虑单个裂隙的形状、规模及排列。
4.1 地下水流动数学模型及其求解方法
二、地下水流数学模型的内容及其求解方法 1 一般包含两方面的内容： （1）有一个或一组能描述地下水运动规律的偏微分方程。一般，可把它 理解为是水量（或质量）平衡方程的一种微分表达形式；【注意对水文地质 条件的分析】 （2）有相应的定解条件。定解条件包括初始条件和边界条件。当模拟稳 定流时，只需要边界条件；当模拟非稳定流时，则要同时给出初始条件和边 界条件。 初始条件是指某一选定的初始时刻渗流区内地下水的水头分布状况；边 界条件是指渗流区的边界位置及边界上地下水应满足的条件，一般包括有已 知水头边界、已知流量边界和混合边界等三类。
第四章 地下水系统分布参数数学模型
4.1 地下水流动数学模型及其求解方法 4.2 地 地下水流动数学模型及其求解方法
一 地下水流数学模型 分布参数地下水流数学模型可分为不同的类型。按赋存介质的不同，
可分为孔隙水、裂隙水和岩溶水等模型；按流动要素的特点，可分为稳 定流、非稳定流、承压流、非承压流、饱和流、非饱和流、单相流、多 相流、一维流、二维流和三维流等模型。
4.1 地下水流动数学模型及其求解方法
2 地下水流数学模型的求解方法 （1）解析法是利用经过严格的推理和数学推演获得的各种计算公
式进行计算的方法，裘布依公式和泰斯公式即属此类。其特点是：在 一定的假设条件下所进行的数学推导是非常严密的，所获得的解是精 确的。但水文地质条件非常复杂，用此法求解，误差非常大。
2 裂隙水 （1）裂隙介质的特点 ①其渗透性具有较强的非均质性和各向异性，在一个计算区内往往可以划出许
多非均值区，每个区之间的导水系数或渗透系数有时可以相差数百倍。 ②渗透性往往随深度发生变化，在不同地点或不同的岩层，渗透系数随深度的
变化规律也不相同。因此，在厚层裂隙介质中确定含水层的厚度比较困难。 ③储水能力比孔隙介质小。 ④裂隙介质的储水性主要取决于微裂隙和岩层的孔隙（如粒间孔隙、气孔等），
有限单元法可分为变分有限元法和迦辽金有限元法。变分有限元法的 基本思想是：利用变分原理，把地下水流偏微分方程的定解问题转化 为某一泛函I[H]的极值问题，再通过剖分和插值，可将该泛函I[H]转 化为一个普通的二次函数。这个二次函数的极值问题等价于一个线性 代数方程的求解问题，求解该方程组，即可得到地下水流问题的解。
（2）数值法是一种只能求出计算区内有限个点在某些时刻近似解 的计算方法。也就是说，数值法不能求出解的解析表达式，只能求出 解（未知量）在若干个离散点的近似值。
数值法在处理渗透区的非均质性和不规则性等方面具有很大的灵 活性。
4.1 地下水流动数学模型及其求解方法
a 有限差分法：原理就是在空间域时间离散的基础上用差商代替微商， 最后把地下水运动的偏微分方程的求解问题转化为差分方程的求解问 题；
4.1 地下水流动数学模型及其求解方法
3 岩溶水
岩溶水与裂隙水在许多方面有共同之处，因而裂隙水的某些特点对岩溶水也同 样使用。但是，岩溶介质比裂隙介质具有更强的非均质性，使建立岩溶水流数学模 型更困难。
北方与南方岩溶水在分布、径流、补给、动态等方面存在着差异。南方岩溶地区， 地下水流动以管道流为主；北方地区，岩溶发育相对较弱，多数裂隙-岩溶水，含水 层中的地下水往往具有统一的地下水面。
而导水性则主要取决于一些较大的裂隙。 ⑤由于压力对裂隙张开的程度有明显的影响，故压力变化与岩层的透水性之间
存在着密切的关系。
4.1 地下水流动数学模型及其求解方法
（2）裂隙介质渗透性研究的途径 ①把裂隙介质视为连续介质，认为裂隙水的运动基本上服从达西定
律，可以用与孔隙水相同的数学模型来计算裂隙水流。 ②仍然认为裂隙水的运动符合达西定律，但地下水的传导和运动主