地下水流数值模拟基础.pptVIP

地下水流数值模拟基础 中国地质大学环境学院 2009.6 与地下水有关的资源、环境、工程问题 地下水资源评价与优化开发 基坑排水预测 坝体渗漏评价 水库渗漏评价 软土路基排水固结沉降预测 地下热水泵、地下空调注采方案评价 海水入侵预测评价 有害废物渗滤液在含水层中迁移范围预测 地下水水位动态对滑坡体稳定性的影响评价 地面塌陷、地面沉降预测 第二章 地下水流动问题概述 2.1 地下水流动微分方程 含水层非均质性与各向异性特征 达西定律 Q=KA(H1-H2)/L=KAI。 式中Q为渗流量，A为过水断面，(H1-H2)为水头损失，L为渗流路径长度，I为水力坡度，K为渗流系数。 关系式表明，水在单位时间内通过多孔介质的渗流量与渗流路径长度成反比，与过水断面面积和水头损失成正比。 达西（渗透）流速与实际流速 达西（渗透）流速：通过某一断面的流量Q等于流速v与过水断面F的乘积，即Q＝Fv，即V=Q/F=KI。 实际水流流动过程中通过的范围是扣除结合水所占据的范围以外的孔隙面积。 实际流速：u＝Q/Fn， 其中，n－有效孔隙度。 两者的关系为：v=nu 渗透系数,K，又称水力传导系数 在各向同性介质中，它定义为单位水力梯度下的单位时间通过单位断面面积的流量。 导水系数,T 承压含水层中的地下水在单位水力梯度作用下，单位时间内通过单位宽度承压含水介质的流量。 T＝KM， 其中K为渗透系数，M为承压含水层厚度。 弹性储水(释水)系数μe 承压含水层中，当水头上升（下降）一个单位时，单位水平面积承压含水层柱体所储存（释放）的水量。 重力给水度μd， 当潜水位下降一个单位时，被疏干含水层单位体积内在重力作用下所能排出的水量。表征了含水层的释水能力。 比弹性储水系数μs μs＝μ/M，M为承压含水层厚度。表征单位体积含水层的储（释）水能力。 Boussinesq假定 在研究潜水的非稳定运动时，假定水是不可压缩的流体，均质岩层中的潜水是缓变运动。 稳定流与不稳定流 水头、渗透速度等任一渗透要素随时间变化的地下水运动为非稳定流，不随时间变化则为稳定流。 地下水水均衡原理 地下水的水量和盐分在收支方面的数量关系，称为地下水均衡。其中水量均衡为水均衡，盐分均衡的话称为盐均衡。 在均衡期中，均衡区的补给量大于排泄量---正均衡。 在均衡期中，均衡区的补给量小于排泄量---负均衡。 地下水的均衡状况是通过建立地下水均衡方程实现的。其原理就是水量平衡原理，一般: △W=X + W1 + Z1 + Y1 + Z2 + W2 + Y2 其中：△W均衡期内地下水量的变化量，X-大气降水的入渗补给量，W1-地下水流入量，Z1-凝结水补给量 ，Y1-地表水入渗补给量，W2-地下水流出量，Z2-地下水蒸发量， Y2-地下水补给地表水量 2.1.2 三维流基本微分方程 2.1.2 三维流微分方程（续1） 2.1.2 三维流微分方程（续2） 2.1.2 三维流微分方程（续3） 2.1.2 三维流微分方程（续4） 2.1.3 承压含水层平面二维流微分方程 2.1.3 承压含水层平面二维流微分方程（续1） 2.1.4 潜水平面二维地下水流动基本微分方程 2.1.4 潜水平面二维地下水流动微分方程（续1） 2.2 边界条件和初始条件 2.2.1 初始条件 2.2.2 边界条件 2.3 地下水运动的数学模型 2.3.1 地下水运动的数学模型结构 2.3.2 地下水运动的数学模型例子 2.3.2 地下水运动的数学模型例子 本章小结 微分方程 定解条件 边界条件 初始条件 已知t=0时的因变量， H(x,y,z,0)=H0(x,y,z) 已知水头边界(I类边界) H(x,y,z,t)=f(x,y,z,t) (x,y,z)?B1 特例：定水头边界 H(x,y,z,t)=C 已知流量边界 特例：隔水边界 数学模型 例：河间地块承压水流模型 设两条河流平行、完全切割乘压含水层，含水层等厚、均质各向同性,无垂向补给或排泄，对于如图所示的坐标系，已知某时刻的含水层各处的水头为20米，自该时刻后，河水位分别为如图所示的函数。试根据条件作合理简化建立其数学模型。 x H(x,t) * * 第一讲 第一章 绪论 第二章 地下水流动问题微分方程及数学模型 （1）建立基本微分方程及数学模型： ● 按空间维数：一维、二维（平面二维、剖面二维）、三维 ● 按含水层类型：承压水流、潜水流、多层（越流联系）等 （2）求解数学模型的数值方法