下载文档原格式
渗流模型知识点总结图渗流模型是描述地下水流动和传输的数学模型，它可以帮助我们理解和预测水在地下的流动情况。
渗流模型可以应用于地下水资源管理、地下水污染治理、水文地质等领域，具有重要的实用价值。
下面是关于渗流模型的一些重要知识点总结。
1. 渗流方程渗流模型的数学描述基于渗流方程，它描述了地下水在多孔介质中的流动规律。
渗流方程通常采用达西定律和杜安-卡丁方程进行描述，它们可以用来描述地下水的渗流速度、渗透率、孔隙度等参数之间的关系。
2. 边界条件在渗流模型中，边界条件是描述模型边界上的地下水流动情况的重要参数。
常见的边界条件包括：Dirichlet边界条件、Neumann边界条件和混合边界条件。
这些边界条件可以帮助我们对地下水流动的边界条件进行准确描述，是渗流模型计算的基础。
3. 初始条件渗流模型中的初始条件是指模型开始计算时的地下水流动情况。
初始条件通常是指地下水位和地下水流动速度的初始数值，它们是模型计算的起点。
在模型计算中，初始条件的准确性对计算结果具有重要影响。
4. 离散化方法为了解决渗流方程，通常需要将其离散化。
常见的离散化方法包括有限差分法、有限元法和边界元法等。
这些方法可以将连续的渗流方程转化为离散的问题，通过计算机进行数值计算，得到地下水流动的数值解。
5. 模型验证渗流模型的验证是指利用现场观测数据来验证模型的准确性和可靠性。
验证通常包括比对模型计算结果和现场观测数据，评估模型的拟合程度，以及对模型参数的敏感性分析等。
模型验证可以帮助我们了解模型的适用范围和局限性，提高模型的预测准确性。
6. 模型应用渗流模型在地下水资源管理、地下水污染治理、水文地质和地下水开采等领域有着广泛的应用。
通过渗流模型，我们可以模拟地下水流动过程，预测地下水位和地下水流向，并为地下水资源的合理开发和保护提供科学依据。
此外，渗流模型也可以帮助我们理解地下水污染的传播规律，优化地下水治理方案。
总的来说，渗流模型是描述地下水流动和传输的重要工具，它可以帮助我们理解地下水资源的分布和变化规律，为地下水资源管理和保护提供科学依据。
第四章 地下水运动的数值模型解析解虽然具有精确可靠的特点，但采用解析解反映自然状态和复杂人类活动干扰下的地下水运动是相当困难的。
因此，当含水层的条件严重偏离现有解析模型的简化假设时，人们通过数值模型来获得近似的地下水流场及演变趋势。
第一节 地下水流数值方法概述地下水流的数学模型采用偏微分方程描述地下水流的时间和空间连续状态，而数值模型则是采用离散（非连续）时空模型中水头的分布与演变对数学模型进行近似描述。
从精确数学模型到近似数值模型的转化，虽然会损失一些精度，但使复杂地下水流问题的分析得以通过机械计算实现，而且误差也是可控的。
把偏微分方程求解的数值方法引入到地下水流问题的求解始于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、承压含水层中的完整井流一、定流量抽水时的Theis(泰斯公式)假定条件：（1）含水层均质各向同性、等厚，侧向无限延伸，产状水平； （2）抽水前天然状态下水力坡度为零； （3）完整井定流量抽水，井径无限小； （4）含水层中水流服从Darcy 定律；（5）水头下降引起的地下水从贮存量中的释放是瞬时完成的。
教学模型的建立和求解：抽水形成以井轴为对称轴的下降漏斗，将坐标原点放在含水层底板抽水井的井轴处，井轴为z 轴，建立坐标系。
则以降深表示的微分方程为：教学模型为： 其解为：其中：s —抽水影响范围内，任一点任一时刻的水位降深；Q —抽水井的流量； T —导水系数；t —自抽水开始导计算时刻的时间； r —计算点到抽水井的距离； u —含水层的贮水系数。
利用上述W(u)和u 的关系制定P 93表，W(u)可查表得。
首先由 计算出u 値，然后查表得相应的W(u)，再求r 处的s 値。
tsT urs r rs ∂∂=∂∂+∂∂*221⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧-=∂∂>=∂∂=∞∞<<=∞<<>∂∂=∂∂+∂∂→∞→T r s r t r s t s r r s r t t s T u rs r r s r r π2Qlim 00,0),(00)0,(0,010*22()()⎰∞-===uy dy y e u W Tt u r u u W T Q s 44*2πTtu r u 4*2=二、流量变化时的计算公式流量随时间的变化，可分为阶梯变化和连续渐变，相应的流量过程线为台阶状和连续光滑的曲线。
连续光滑的曲线应概括成阶梯状折线。
概化的原则：矩形面积等于曲线与横坐标所围成的面积。
每一个阶梯视为定流量，用Theis 公式计算降深,然后将各降深叠加起来，得流量变化的总降深値。
如图，连续抽水，概括为4个阶梯，若求t 时刻的水位降深，则可分解为四个亚问题，第一个亚问题以Q 1流量抽水从t 0→t ；第二个亚问题以Q 2－Q 1流量抽水，从t 1→t ；第三个亚问题以Q 3－Q 2流量抽水从t 2→t ；第四个亚问题以Q 4－Q 3流量抽水从t 3→t 。
第四章 地下水运动的数值模型解析解虽然具有精确可靠的特点，但采用解析解反映自然状态和复杂人类活动干扰下的地下水运动是相当困难的。
因此，当含水层的条件严重偏离现有解析模型的简化假设时，人们通过数值模型来获得近似的地下水流场及演变趋势。
第一节 地下水流数值方法概述地下水流的数学模型采用偏微分方程描述地下水流的时间和空间连续状态，而数值模型则是采用离散（非连续）时空模型中水头的分布与演变对数学模型进行近似描述。
从精确数学模型到近似数值模型的转化，虽然会损失一些精度，但使复杂地下水流问题的分析得以通过机械计算实现，而且误差也是可控的。
把偏微分方程求解的数值方法引入到地下水流问题的求解始于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,2,3,,;;;1,2,3,,;;k k k k k p p p M t t p M t t t H H +=⋅⋅⋅=→=⋅⋅⋅=+∆ （4.1.2） 式中：k 为时刻；t 为时间；k t ∆为时间步长。
数值模型的核心部分，是在地下水流数学模型的基础上确定节点在流场中的相互关系，并形成如下形式的方程组11111212,1,2,k k k k kp M p p p pp pM p M C C C C H H H H F +++++++⋅⋅⋅+==⋅⋅⋅⋅⋅⋅+（4.1.3） 式中：系数pL C 反映了节点P 与节点L 在流场中的关系，可称为关联系数；k p F 通常是由前一个时刻决定的常量。
式（4.1.3）除包括地下水流控制方程的离散近似以外，还包括了边界条件的离散化处理。
当k=0时，k p F 中包含了初始条件。
如果模拟稳定流，则k p F 与时间无关，上标k 可以略去。
使用矩阵符号，式（4.1.3）可缩写为[]{}{}M M C H F ⨯= （4.1.4）实际在处理模型时，一个节点通常只与相隔最近的若干节点有直接的关系，在式（4.1.3）中大多数关联系数为零，而且一般有对称性，即PL LP C C =。
因此，矩阵[]C 通常是对称的稀疏矩阵，降低了方程组求解的困难。
尽管如此，区域地下水流数值模型的节点数往往很大，而且在某些非线性条件下，[]C 本身是待求水头的函数，方程组的求解需要使用迭代法而耗费相当多的时间。
不同类型的地下水流数值模型，最核心的差异在于形成式（4.1.3）的方式，关联系数具有不同的计算公式。
原则上，对于完全相同的地下水流数学模型，以及完全相同的离散节点序列和离散时间序列，不同的数值模型应给出近乎同样的结果，只是计算精度略有差异。
已有研究表明，某种类型的有限差分法与某种类型的有限元法存在一定的等价性，但更多的对比分析尚待进一步探讨。
实际上，研究者在检验数值模型的可靠性时，总是与解析解的模型对比。
当不同方法的数值模型用相同的解析解作验证时，他们之间的等价性就可以得到一定程度的证明。
然而，人们对数值模型的偏好并不取决于数值模型之间的等价性，而往往取决于某种数值方法的熟悉程度，或者基于某种数值方法专业软件的方便快捷。
专业软件在地下水流数值模型的推广应用中发挥了重要作用，如有限差分模拟程序MODFLOW （McDonald and Harbaugh ，1988）和HST3D （Kipp ，1987）、有限元模拟程序FEMWATER （Lin et al ．，1997）和FEFLOW （Diersch and Kolditz ，1998）等，国内有陈崇希等开发的有限差分模拟程序PGMS （陈崇希等，2007）。
其中，以美国地质调查局发布的有限差分模拟程序MODFLOW 最为著名，本章最后将介绍其特点。
第二节 地下水流有限差分模型有限差分法是地下水流数值模拟最早使用的一种方法。
这种方法古典而简洁、物理意义明了、编程容易，因而被广泛流传。
地下水流的有限差分模型已经从最初的一维均质模型发展为目前的三维非均质模型，并有专业软件加以实现。
一、有限差分法的基本原理有限差分法是求解偏微分方程边值问题和初值问题的一种数值方法，其实质是把连续的模型空间离散化为规则或不规则的网格点，利用导数的差分近似形式代替偏微分方程形成差分方程组，通过求解方程组得到离散点的待求变量作为连续场的一种近似结果。
模型空间的离散化从形态上主要分为规则网格与不规则网格。
在二维平面空间，规则网格由一系列与坐标轴平行的直线组成，直线交叉形成的单元（cell ）为矩形；在三维空间，规则网格由一系列与坐标轴正交的平面组成，平面相互切割形成的块体（block ）为长方体。
不规则网格在平面上由一系列互不重叠的多边形填充模型空间而形成，一般通过分层建立三维空间网格。
一维问题空间离散化很简单，就是在坐标轴上生成一系列具有一定间隔的点，或称节点（node ），相邻两个节点之间为线段，待求变量在这些节点上的分布代表了在坐标轴上连续分布的近似状态。
有限差分法中的“有限”是指网格中的节点、单元或块体数目是有限的，每个线段、单元或块体的尺度也是有限的。
首先用一维问题来阐述导数的差分格式。
设待求变量f 沿x 轴的分布为二阶连续函数()f x ，现将x 轴离散化为节点012(,,,,,,)i N x x x x x ⋅⋅⋅⋅⋅⋅，它们的待求变量为012(,,,,,,)i N f f f f f ⋅⋅⋅⋅⋅⋅，则()f x 在节点i x 处的一阶导数表示为[]11i i i i ii x x f f fx x x x ++=-∂=+∞∆∂- （4.2.1） 式中：1i i i x x x +∆=-，为节点的间距。
去掉式（4.2.1）中的Taylor 展开式余项，就得到一阶导数前向差分格式。
一阶导数还可以表示成后向差分格式为11i i i i i x x f f fx x x --=-∂≈∂- （4.2.2） 式中：11i i i x x x --∆=-。
显然，第一个节点0x 处的一阶导数只能表示为前向差分格式，而最后一个节点N x 处的一阶导数只能表示为后向差分格式。
有了一阶导数的差分格式之后，再次使用中心差分格式，就可以建立二阶导数的差分格式为211211()()i i i i i x x f f fx x x x x +-+-=∂∂-∂∂∂∂-≈（4.2.3）其中：一阶导数在节点1i x +采用后向差分格式，在节点1i x -采用前向差分格式，则可以得到二阶导数的差分格式为211111121111()()()()()()i i i i i i i i i i i i i i i i x x x x f x x f x x f fx x x x x x x -++-+-+-+-=---+-∂∂---≈（4.2.4） 如果相邻节点的间距都相等，即1i i x x x -∆=∆=∆，则上式可简化为2112222i i i i x x f f f fx x +-=-+∂∂∆≈（4.2.5） 建立了一阶和二阶导数的差分格式，就可以用它们来改写所研究问题的偏微分方程。
下面用一个简单的例子来说明有限差分模型的建立方法。
设有一边值问题的数学模型为220,02d f df a b c x L dx dx++=<< （4.2.6） (0)f x u == （4.2.7）2x Ldfv dx == （4.2.8）在建立有限差分模型时，把模型空间[0，2L]离散化为三个节点013(0,,2)x x L x L ===。
对第二个节点1x 建立控制方程（4.2.6）的差分格式为210212202f f f f f a b c L L-+-++= （4.2.9） 第一个节点0x 为边界节点，直接利用边界条件式（4.2.7）得0f u = （4.2.10）第三个节点2x 也是边界节点，把边界条件式（4.2.8）表示成差分格式为21f f v L-= （4.2.11） 这样，由式（4.2.6）至式（4.2.8）描述的数学模型就转化为了由方程组[式（4.2.9）至式（4.2.11）]构成的差分模型。
求解这个差分方程组，得22220122222()(2)bL cL bL cL f u f L v u f L v u a a a a==---=---，， （4.2.12） 待求变量的差分模型结果012()f f f ，，就反映了数学模型精确解()f x 得特征。
如果我们遇到的是一个初边值问题，例如，把上述关于()f x 的边值问题改写为如下关于()f x t ，的数学模型2202f f f a b c e x L x x t∂∂∂++=<<∂∂∂， （4.2.13） (0,0)f x t u =>= （4.2.14）2(0)x Lf x t v x =∂>=∂， （4.2.15） (0)00f x t x L ==≤≤，， （4.2.16）其中，式（4.2.16）是初始条件。
对于这种初边值问题，我们需要先解决关于时间的偏导数。
把时间坐标轴离散化为时步012M ()k t t t t t ⋅⋅⋅⋅⋅⋅，，，，，，，则第k 时刻的控制方程可以写为21202k k k k kf f f f a b c e x L x x t -∂∂-++=<<∂∂∆， （4.2.17） 式中：1k k k t t t -∆=-为时间步长。
时间偏导数取后向差分格式，方程左侧偏导数的自变量全部取第k 时刻的数值，这种时间上的处理称为隐式差分格式。
