描述地下水运动的数学模型及其解法课件

1、1.7 描述地下水运动的数学模型及其解法7/27/20221一、数学模型的有关概念 数学模型：描述某一研究区地下水流运动的数学方程与其定解条件共同构成的表示某一实际问题的数学结构。亦即从物理模型出发，用简洁的数学语言，即一组数学关系式来刻画它的数量关系和空间形式，从而反映所研究地质体的地质、水文地质条件和地下水运动的基本特征，达到复制或再现一个实际水流系统基本状态的目的的一种数学结构。 其中微分方程表示地下水的流动规律，定解条件表明研究对象所处的特定环境条件，即所研究的地下水流的真实状态。 定解问题是给定了方程（或方程组）和相应定解条件的数 学物理问题。 建立模型是指建立数学模型的过程。7/2

2、7/20222二、描述渗流问题的数学模型 前面介绍的承压含水层、越流系统和潜水含水层中渗流基本微分方程，称为渗流的支配方程（Governing equation）：反映了含水层系统中各变量之间的关系。 一个含水层系统中的渗流状况不仅取决于各变量之间的关系，而且还取决于渗流区域的边界条件和初始条件，边界条件和初始条件统称为含水层系统渗流的定解条件。7/27/20223 支配方程数学模型 初始条件 定解条件 边界条件 确定型(学习的重点)常见数学模型 随机型7/27/20224三、定解条件.初始条件（Initial conditions） 初始条件：是指给定的某选定时刻（通常t=）渗流区内各点的水

3、头H的分布情况。 H0，H0为空间和平面区域上的已知函数。如研究平面问题，某一时刻所测的等水头线即可作为初值。7/27/20225注意： a、初始条件对计算结果的影响，随计算时间的延长而减弱； b、初始条件并非地下水的原始状态或未开发以前的状态； c、初始条件可根据需要任意选取。7/27/20226.边界条件（Boundary conditions） 指渗流区域几何边界上的水力性质（所处的条件），用来表示水头或渗流量在渗流区边界上所应满足的条件，也就是渗流区水流与其周围环境相互制约的关系。 （1）第一类边界条件 若在渗流区的某部分边界上各点在每一时刻的水头是已知的，则称这部分边界为第一类边界或

4、给定水头边界，常表示为：分别表示在三维和二维条件下边界上的点在t时刻的水头分别表示三维和二维条件下边界上的水位已知函数7/27/20227常见的第一类边界有： 河流或湖泊切割含水层，二者有密切的水力联系，此时，河湖的水位是已知的，水头 或 是由河湖水位的统计资料得到的关于t的函数；但要注意，某些河、湖底部及两侧沉积有一些粉砂、亚粘土和粘土，使地下水和地表水的直接水力联系受阻，就不能作为第一类边界条件来处理。 在自然界，这种情况很少见。就是附近有河流、湖泊，也不一定能处理为定水头边界，还要视河流、湖泊与地下水水力联系的情况，以及这些地表水体本身的径流特征而定。在没有充分依据的情况下，不要随意把某

5、段边界确定为定水头边界，以免造成很大误差。7/27/20228泉水溢出带：其标高即为水位资料，但必须保证溢出带不消失；区域的抽水井、注水井或疏干巷道也可作为给定水头边界处理；无限边界 亦为第一类边界；潜水面任一点的水位已知时，抽水井井壁水位为一类边界。7/27/20229 （2）第二类边界条件 当已知渗流区某部分边界上的流量分布时，称这部分边界为第二类边界或给定流量边界。相应的边界条件表示为： 式中：n为边界 或 的外法线方向； q1和q2为已知函数，分别表示 上单位面积和 上单位宽度的侧向补给量。 7/27/202210常见的二类边界条件有（给定流量边界）: 隔水边界为零通量边界： 地下分水

6、岭：无水流通过，也为零通量边界； 流线为零通量边界； 研究区内的抽水井或注水井，需要时作为内边界，取井壁为 第二类边界： 由达西定律可知： 式中：r为径向距离； Q为抽水井流量（Q0为注水井流量） 由于此时外法线方向n趋近于井中心，所以上式写为：7/27/202211（3）第三类边界条件（混合边界）若某边界 或 上的和 的线性组合为已知，即 其中： 为已知函数，这种类型的边界条件称为第三类边界条件或混合边界条件。注意： 1求解非稳定渗流问题时，数学模型应包括支配方程、初始条件和边界条件。 2对于稳定流问题，数学模型仅有支配方程和边界条件。7/27/202212三、建立数学模型的实例例 潜水非稳定流方程 W7/27/202213例2 设想如图所示的水文地质条件 设w（x , y ,t）为单位时间、单位面积上的补给量，即补给强度，P (x , y, t)为开采区单位面积上的抽水量，试写出其数学模型。7/27/202214Hz7/27/202215四、数学模型的解法A.解析法：适用于边界条件简单的地区，一般径流 影响范围为200300m，其解为精确解，常用来验 证数值解是否正确。B.数值法：适用于较复杂的边界条件和数学模型。有限单元法、有限差分法、边界元法其解为近似解，所以模型需