地下水数值模拟模型简介
有限差分法(Finite Difference Method)
有限单元法(Finite Element Method)
Numerical Methods
数值法是一种近似解法。
它以地下水运动的微分方程的定解问题为基础，将表示水位随时
间和空间连续变化的函数离散化，求得函数在有限节点（或结点） 上的近似值。
系统的参数和输入来表示地下水系统的输出，
也称为建立地下水模拟模型。
建模目的
识别 预报 探测
Solving GW Flow Equation
求解方法：
解析法：裘布依公式、泰斯公式（条件和局 限性） 2 Q μr s w ( ) 4 T 4 Tt 电网络方法(已被淘汰)
数值法
3D Groundwater Flow Equation
h h h h (hs h) Ky Kz S Kx t x x y y z z 2 2 2 h h h h h K Kz Kz p p Ky x t x z z y h(x, y, z,t) t0 h0 h 0 n 1 h Kn 2 q(x, y,t) n x, y, z , t 0 x, y, z 0, t 0 x, y, z , t 0 x, y, z 1, t 0 x, y, z 2, t 0
Groundwater Model_Math Model
变量：
自变量（输入变量、决策变量） 因变量（输出变量、状态变量） K、T、K’等 结构（几何参数） 初始条件 边界条件：边界类型、边界值
参数
定解条件
Solving GW Flow Equation
求解：即建立地下水系统状态方程，用地下水
Problems about GW Development
Groundwater System
Structure: Aquifer, Aquitard (distribution, Parameters) Input: Recharge, Discharge Output: Water table, Chemical Concentration,…
Problems about GW Development
Functions of GW
An Important Water Supply for Human An Important Environmental Factor
Water Shortage Ecologic problems Saline soil Sea Water Intrusion Land Subsidence Groundwater Contamination ……
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
Groundwater Mo GW Development
How Can Control GW System? Groundwater Model
To Calibrate the System To Predict the Behave To Control or Manage the System
地下水数值模拟模型简介
For Group of Groundwater Resources
By Shao Jingli Institute of Water Resources and Environment China University of Geosciences, Beijing
Contents
Groundwater Model_Math Model
2D, 2Layers Groundwater Flow Equation
h1 h1 h1 Kz S ( h b ) K ( h b ) K (h2 h1) Q 1 1 1 1 1 1 x, y 1, t 0 x y y L t x h2 h2 h2 Kz S ( MK ) ( MK ) (h2 h1) Q2 x、y 2, t 0 2 2 2 t x x y y L 0 h1(x, y, t) t0 h1 (x, y) x, y 1, t 0 0 x, y 2, t 0 h2 (x, y, t) t0 h2 (x, y) h1 (h1 b)K1 x, y , t 0 q 1( x, y, t) n h2 x, y , t 0 MK q2 (x, y, t) 2 n x, y 1, t 0, hs h1 smax ε ε0( (hs h1)
What is Groundwater Model? Problems about GW Development Groundwater Model_Math Model Solving GW Flow Equation Numerical Methods Calibration of the Model
《地下水数值模拟》课件
CHAPTER 04
地下水数值模拟的案例分析
案例一：某地区地下水污染模拟
总结词
该案例展示了如何运用地下水数值模拟技术 预测和评估某地区地下水污染情况。
详细描述
该案例首先介绍了该地区的地下水分布和流 向，然后通过建立数值模型，模拟了不同污 染源对地下水的影响，并预测了污染扩散的 范围和程度。最后，根据模拟结果，提出了 相应的污染防治措施。
VS
有限体积法适用于不规则的网格系统 和复杂的边界条件，能够得到相对准 确的结果，计算量适中，适用于较大 的模型规模。
CHAPTER 03
地下水数值模拟的步骤
建立数学模型
01
确定研究区域和边界条件
02
描述地下水流动和物质传输过程
03
建立数学方程，包括连续性方程、动量方程、源汇 项等
模型离散化
1
地下水数值模拟的应用
地下水数值模拟广泛应用于水资源管理、环境保护、地质 灾害防治等领域。
通过模拟地下水动态变化，可以预测未来地下水资源量、 评估地下水污染风险、研究地下水与地质灾害的关系等， 为相关决策提供科学依据。
CHAPTER 02
地下水数值模拟的基本方法
有限差分法
有限差分法是一种将偏微分方程离散 化为差分方程的方法，通过在时间和 空间上将偏微分方程近似为差分方程 ，从而将连续的物理量离散化为离散 的数值。
随着数值计算技术的发展，地下水数值模型将越来越复杂，能够 模拟更多的物理过程和化学反应。
参数优化和数据同化
通过人工智能和机器学习技术，对模型参数进行自动优化和数据同 化，提高模拟精度和可靠性。
多尺度模拟
从微观到宏观的多尺度模拟将成为一个重要方向，能够更好地揭示 地下水系统的复杂性和规律性。
地下水数值模拟
一、模型概化-内部结构
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、研究区边界
研究区应尽可能以自然边界为计算边界，最好是以 完整的水文地质单元作为计算区。
地下水数值模拟基础
地下水数值模拟基础
地下水数值模拟是指通过建立数学模型，模拟地下水流动和污染传输
的过程，以预测和评估地下水资源的利用和保护。
在地下水资源管理
和环境保护中，地下水数值模拟是一种重要的工具。
地下水数值模拟的基础是地下水流动方程和质量传输方程。
地下水流
动方程描述了地下水的流动过程，包括水流速度、水位变化等；质量
传输方程描述了污染物在地下水中的传输过程，包括污染物浓度变化等。
通过建立这些方程的数学模型，可以模拟地下水的流动和污染传
输过程。
地下水数值模拟的过程包括模型建立、参数确定、模拟计算和结果分
析等步骤。
模型建立是指根据实际情况，建立地下水流动和污染传输
的数学模型。
参数确定是指确定模型中的各项参数，包括地下水的渗
透系数、孔隙度、水头等参数，以及污染物的扩散系数、降解速率等
参数。
模拟计算是指利用计算机对模型进行数值计算，得到地下水流
动和污染传输的结果。
结果分析是指对模拟结果进行分析和评估，以
确定地下水资源的利用和保护策略。
地下水数值模拟在地下水资源管理和环境保护中具有重要的应用价值。
它可以预测地下水资源的开发和利用情况，评估地下水资源的可持续
利用性，指导地下水资源的合理开发和利用。
同时，它也可以预测地下水污染的扩散和影响范围，评估污染物的风险和危害程度，指导地下水环境的保护和治理。
总之，地下水数值模拟是一种重要的工具，它可以为地下水资源管理和环境保护提供科学依据和技术支持。
在未来的发展中，地下水数值模拟将继续发挥重要作用，为地下水资源的可持续利用和环境保护做出更大的贡献。
地下水溶质运移数值模型
地下水溶质运移数值模型(资料性附录)水是溶质运移的载体，地下水溶质运移数值模拟宜在地下水流场模拟基础上，因此地下水溶质运移数值模型包括水流模型和溶质运移模型两部分。
DJ 地下水水流模型非均质、各向异性、空间三维结构、非稳定地下水流系统:1)控制方程σ∂h ∂hy 3(“∂h}∂(∂h ∖S,—=—K v —+—K Y —+—K ——+/∂t 3xI ∂x)为('∂y JAzI ~∂z)式中：SS 一一给水度［I/］；h --- 水位［1］；Kχf Ky,Kz ——分别为X,y,Z 方向上的渗透系数［EΓ∣］;T 一一时间［T ］；Qs 一一源汇项m注：方括号［］中的符号为量纲，以下同。
2)初始条件h(x y y 9z y t)=Zz 0(x,y,z)(x,y,z)∈Ω,/=O 式中：4*，y ，z)——已知水位分布：Q ——模型模拟区。
3)边界条件：第一类边界： 〃(x,y,z√)∣「=Λ(x,y,z√)(x,y,z)∈Γ1,r≥O式中：r '一一类边界； h(x,y,z,t)一一类边界上的己知水位函数。
第二类边界：式中：「2 --- 二类边界;∂nq(x,y,Z) (x,y,z)∈Γ2κ——三维空间上的渗透系数张量;nn——边界r2的外法线方向；q(x,y t z)——二类边界上已知流量函数。
第三类边界：r(k(h-z)-+ah)=q(x,y,z)加r3式中：0一一系数；「3一—二类边界；k一一三维空间上的渗透系数张量；n——边界G的外法线方向;q(x,y f z)——三类边界上已知流量函数。
D.2地下水水质模型1)控制方程R啜喘［吗(他C)Fe—/〜元式中：R——迟滞系数，无量纲Pb SC~Θ~∂Cph——介质密度IM1-3］；θ——介质孔隙度，无量纲；C——组分的浓度［M1，］；亍一一介质骨架吸附的溶质浓度［M1，］；t——时间［T］；X,y,Z一—空间位置坐标［1］；Dij——水动力弥散系数张量［1?T」］；Vi——地下水渗流速度张量［EΓ∣］;q s——源和汇［T∣］;CJ一一源或汇水流中组分的浓度［M1"；4一一溶解相一级反应速率［T」］；4一一吸附相一级反应速率［Tj］。
地下水系统模拟与数值模拟方法
地下水系统模拟与数值模拟方法地下水系统是指地下水的流动、贮存和分布所构成的地下水环境系统。
地下水资源是人类生存所必需的一种重要自然资源，对生态环境和社会经济发展有着重要意义。
在地下水资源的管理和保护过程中，需要对地下水系统进行模拟和数值模拟，以便更准确地预测和评估地下水系统的水文地质特征，判断地下水资源的开发潜力和合理利用方案，为实现地下水资源的可持续利用奠定科学依据。
地下水系统模拟方法地下水系统模拟是指通过对地下水系统进行数学模型的建立和仿真，以预测和分析地下水流动、污染传递等水文地质过程的方法。
常用的地下水系统模拟方法有分析模型和数值模型两种。
分析模型是基于对地下水流动或污染物扩散方程的解析求解，分析模型简单、易懂，计算速度快，但只适用于简单的地下水流动或污染扩散问题。
其主要方法包括平衡分析法、线性和非线性回归分析法、空间分析法等。
数值模型则是通过计算机技术，以数值方法求解数学模型的过程，将地下水系统划分成有限的离散单元，通过离散化的方法，将求解大型、复杂的地下水流动或污染扩散方程转化为大量小规模的计算，从而得到地下水流动或污染传输的详细状况。
常用的数值模型包括有限差分法、有限元法、边元法等。
数值模拟方法的分析优势相比分析模型，数值模型在复杂的地下水流动和污染扩散问题中表现出更强的分析优势。
一是用途广泛。
数值模型可以应用于各种类型的地下水问题，如地下水资源、污染物传输、地下水入渗、河流与地下水交互作用、地下水流场演变分析等。
二是精度高。
数值模型可以准确地反映地下水系统的水文地质特征，得到非常细致的地下水流动和污染扩散情况，为分析和预测地下水资源的分布和变化趋势提供了更多的信息。
三是可视化强。
数值模型的结果可以通过数据处理和可视化技术轻松呈现，利用图表、三维可视化等手段，可以帮助决策者更加直观地了解地下水系统、掌握地下水资源的动态变化。
四是模型灵活性高。
数值模型可以对不同地区、不同时间段的地下水系统进行模拟和分析，从而更好地应对不同地区、不同年份的地下水管理、调控和保护问题。
地下水数值模拟模型简介
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
地下水数值模拟的研究与应用进展
地下水数值模拟的研究与应用进展1. 引言1.1 地下水数值模拟简介地下水数值模拟是指利用数学模型和计算机技术对地下水系统进行模拟和预测的方法。
通过模拟地下水系统的水文地质特征、水文动力过程和水文化学过程，可以更好地理解地下水运动规律，预测地下水资源的变化趋势，指导地下水资源的合理开发和利用。
地下水数值模拟的基本原理包括建立地下水数学模型、确定模型参数、选择数值计算方法、进行模拟计算和模拟结果分析。
地下水数值模拟常用的模型包括地下水流模型、地下水热盐模型、地下水污染迁移模型等，可以根据实际问题的不同选择合适的模型进行建模。
地下水数值模拟在水资源管理、环境保护、地质灾害防治等领域有着重要的应用价值。
通过地下水数值模拟，可以预测地下水位变化、地下水资源补给和排泄规律，为科学合理地开发利用地下水资源提供参考依据。
地下水数值模拟还可以用于评估地下水污染风险、指导地下水污染防治，保护地下水资源环境。
地下水数值模拟是一种强大的工具，为研究人员提供了深入理解地下水系统运行机制和分析地下水问题的方法。
通过不断地研究和应用，地下水数值模拟将在未来发展中发挥更加重要的作用。
1.2 地下水数值模拟的重要性地下水作为重要的水资源之一，对人类生存和发展具有重要意义。
地下水数值模拟是研究地下水流动规律和预测地下水变化的重要手段。
其重要性主要体现在以下几个方面：1.优化地下水资源管理：地下水数值模拟可通过对地下水流动模式的研究和模拟，优化地下水资源的开发和利用。
通过模拟可以更好地预测地下水位变化、水质变化等情况，有助于科学合理地规划地下水资源的开发和利用方案。
2.保护地下水环境：地下水数值模拟可以帮助研究人员识别地下水受到威胁和污染的情况，从而采取合适的措施进行保护和修复。
通过模拟可以及时发现地下水受到污染的源头和扩散路径，指导环境保护工作的开展。
3.灾害预警和防范：地下水数值模拟可以用于预测地下水位变化、地下水涌出、地下水泛滥等情况，为灾害预警和防范提供科学依据。
地下水数值模拟的研究与应用进展
地下水数值模拟的研究与应用进展地下水数值模拟是指利用计算机和数学模型对地下水系统进行模拟和预测的一种方法。
通过数值模拟，可以预测地下水的水位、水量、水质以及地下水与地表水和土壤水之间的相互作用等情况。
本文将探讨地下水数值模拟的研究和应用进展。
地下水数值模拟的研究主要集中在以下几个方面：第一，模型建立。
地下水数值模拟的第一步是建立数学模型。
常用的地下水数值模型有有限元法、有限差分法和边界元法等。
这些方法可以将地下水系统离散化，并通过计算机模拟地下水系统的运动规律。
第二，参数估计。
地下水数值模拟需要大量的参数来描述地下水系统的特性，如孔隙度、渗透率、水头等。
这些参数往往难以直接测量，需要通过试验或其他方法进行估计。
目前研究者们提出了一系列参数估计的方法，如反问题求解、遗传算法等。
数值算法。
由于地下水系统的非线性和复杂性，数值模拟需要高效、稳定的算法来求解方程。
近年来，随着计算机技术的进步，地下水数值模拟中出现了一些新的数值算法，如并行计算、多核计算等。
除了研究方面，地下水数值模拟也广泛应用于实际工程和科学研究中。
以下是一些地下水数值模拟的应用进展：第一，地下水资源管理。
地下水是重要的水资源，然而地下水资源的开发和利用存在一定的风险。
通过地下水数值模拟，可以模拟地下水系统的响应，帮助决策者科学地管理地下水资源，避免资源的过度开发和污染。
第二，地下水污染控制。
地下水污染是当前面临的重要环境问题之一。
通过地下水数值模拟，可以对地下水污染的来源、传输和演化进行模拟和预测，为地下水污染控制提供科学参考。
地下水排水和灌溉。
地下水数值模拟可以帮助工程师科学地设计地下排水和灌溉系统，提高系统的效率和可靠性。
通过模拟地下水的水动力行为，可以优化排水和灌溉的方案，减少水资源的浪费。
第四，地下水地热利用。
地下水中的热量可以被用于供暖和制冷，被广泛应用于地热能利用。
地下水数值模拟可以模拟地热系统的热量传递过程，优化地下水热交换器的设计，提高地热能利用效率。
地下水数值模拟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 (
地下水数值模拟任务步骤及常用软件
地下水数值模拟任务步骤及常用软件地下水数值模拟是指通过建立数学模型和运用计算机方法，利用计算机模拟地下水的水文过程，预测地下水的动态变化，并定量分析地下水资源的开发利用。
地下水数值模拟在地下水资源管理、环境保护、地下水污染防治等领域具有广泛的应用。
1.建立地下水数学模型：根据地下水的特征和要研究的问题，建立合适的数学方程和边界条件，描述地下水系统的基本运动规律。
2.选择合适的计算方法：根据模型的特征和要求，选择合适的数值计算方法，如有限差分法、有限元法、边界元法等。
3.模型参数的确定：对于地下水数学模型中的一些参数，如渗透率、初始压力等，需要通过现场实测或实验室测试获得，并进行合理的插值和外推处理。
4.数值模拟的实施和验证：利用计算机软件进行数值计算，模拟地下水系统的动态变化，并通过对模拟结果的与实测数据的比较，验证模型的可靠性和准确性。
5.模型的应用和优化：在模型建立和验证的基础上，利用模型进行不同方案的对比研究，优化地下水资源的管理和利用方式。
1.MODFLOW：是美国地质调查局开发的地下水流动模型，是目前最常用的三维地下水数值模拟软件之一、具有强大的建模和计算功能，可以模拟各种地下水问题。
2. FEFLOW：是德国DHIGmbH公司开发的强大的地下水和污染物运移模拟软件，可模拟多孔介质中的多个相（水、气和污染物）的运动和相互作用，广泛应用于地下水资源管理和环境保护领域。
3.MODPATH：是美国地质调查局开发的地下水路径分析软件，可以模拟地下水流动路径，并用于评估污染物传输路径和确定水源保护区等。
4.SEAWAT：是美国地质调查局开发的海岸带地下水模拟软件，结合了MODFLOW和MT3DMS，可以模拟地下水和盐水的运动、混合和溶解反应等。
5. GMS（Groundwater Modeling System）：是美国Aquaveo公司开发的集成地下水模型软件平台，集成了多个地下水模型的功能和算法，提供了友好的图形界面和强大的后处理功能。