地下水运动的数学模型

《地下水运动的数学模型》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、第四章地下水运动的数值模型解析解虽然具有精确可靠的特点，但采用解析解反映白然状态和复杂人类活动干扰下的地下水运动是相当困难的。因此，当含水层的条件严重偏离现有解析模型的简化假设吋,人们通过数值模型來获得近似的地下水流场及演变趋势。第一节地下水流数值方法概述地下水流的数学模型采用偏微分方程描述地下水流的时间和空间连续状态，而数值模型则是采用离散（非连续）时空模型川水头的分布与演变对数学模型进行近似描述。从精确数学模型到近似数值模型的转化，虽然会损失一些精度，但使复杂地下水流问题的分析得以通过机械计算实现，而且误差也是可控的。把偏微分方程求解的数值方法引入到地下水流问题的求解始于20世纪

2、70年代，主要方法包括•有限差分法、有限元法和边界元法，此后又发展了有限分析法、多重网格法和无网格法等。这些方法的共同特点是将模型空间及边界离散为由一系列的节点以及联系这些节点的单元（无网格法除外），含水层的水头在这些节点上定义，从而实现了水头分布空间连续函数向离散变量的转化，表示为d2fkNa—^r^b^—^c=edx^oxtkfkw-1—,0

3、)=>{巧;p=1,2,3,…,M}式中；H为含水层的水头；x、y、z为空间坐标；p为数值模型的节点；M为节点的数目。与此同时，时间也被离散化为一系列具有先后顺序的时刻，不同时刻节点的水头分布方式也不同，每个时刻节点的水头分布结果总是由前一个时刻演变过来，即｛上人；pi，2,3,…,T｛£/：「;/?=1,2,3,…=4+△—；｝（4.1.2）式中：k为时刻；t为时间；为时间步长。数值模型的核心部分，是在地下水流数学模型的基础上确定节点在流场屮的相互关系，并形成如下形式的方程组…+C丽H：*F；，pi，2…,M（4.1.3）式中：系数反映了节点p与节点l在流场中的关系，可称为关联系数

4、；F；,通常是由前一个时刻决定的常量。式（4.1.3）除包括地下水流控制方程的离散近似以外，还包括了边界条件的离散化处理。当k=o时，中包含了初始条件。如果模拟稳定流，则F打与时间无关，上标k可以略去。使用矩阵符号，式（4.1.3）可缩写为(4.1.4)［CU{H}={F}实际在处理模型吋，一个节点通常只与相隔最近的若干节点有直接的关系，在式(4.1.3)中大多数关联系数为零，而且一般有对称性，即Cpl=Clp°因此，矩阵［C］通常是对称的稀疏矩阵，降低了方程组求解的困难。尽管如此，区域地下水流数值模型的节点数往往很大,而且在某些非线性条件下，［C］本身是待求水头的函数，方程组的求解

5、需要使用迭代法而耗费相当多的时间。不同类型的地下水流数值模型，最核心的差异在于形成式(4.1.3)的方式，关联系数具有不同的计算公式。原则上，对于完全相同的地下水流数学模型，以及完全相同的离散节点序列和离散时间序列，不同的数值模型应给出近乎同样的结果，只是计算精度略有差异。已有研究表明，某种类型的有限差分法与某种类型的有限元法存在一定的等价性，但更多的对比分析尚待进一步探讨。实际上，研究者在检验数值模型的可靠性时，总是与解析解的模型对比。当不同方法的数值模型用相同的解析解作验证时，他们之间的等价性就对以得到一定程度的证明。然而，人们对数值模型的偏好并不取决于数值模型Z间的等价性，而往

6、往取决于某种数值方法的熟悉程度，或者基于某种数值方法专业软件的方便快捷。专业软件在地下水流数值模型的推广应用中发挥了重要作用，如有限差分模拟程序MODFLOW(McDonaldandHarbaugh,1988)和HST3D(Kipp,1987)、有限元模拟程序FEMWATER(Linetal.,1997)和FEFLOW(DierschandKolditz,1998)等，国内有陈崇希等开发的有限差分模拟程序PGMS(陈崇希等，2007)。其屮，以美国地质调查局发布的有限差分模拟程序MODFLOW最为著名，本章最后将介绍其特点。第二节地下水流有限差分模型有限差分法是地下水流数值模拟最早使

7、用的一种方法。这种方法古典而简洁、物理意义明了、编程容易，因而被广泛流传。地下水流的有限差分模型已经从最初的一维均质模型发展为目前的三维非均质模型，并有专业软件加以实现。一、有限差分法的基本原理有限差分法是求解偏微分方程边值问题和初值问题的一种数值方法，英实质是把连续的模型空间离散化为规则或不规则的网格点，利用导数的差分近似形式代替偏微分方程形成差分方程组，通过求解方程组得到离散点的待求变量作为连续场的一种近似结果。模型空间的离散化从形态上主要分为规则网格