地下水动力学地下水流基本微分方程及定解条件1.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 近似处理成均质的，然后建立每个单元地下水流动的关系 式。把形状不规则的、非均质的问题转化为形状规则的均质 问题。根据研究需要，确定单元划分数量，对于非稳定流还 要对时段进行划分。最后，把局部整合起来，加上定解条 件。 数值解(Numerical solution)是 用数值方法求得的数值解， 是一种近似解。 (3)模拟法：利用物理现象与水流的相似性，在实验室内采用 模拟的方法求解。 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 对于等厚承压含水层，且属于平面二维流 Txx和Tyy为主方向的含水层导水系数(L2/T)；M为承压含水层厚度(L)； ?*为承压含水层的储水系数或弹性给水度。 其物理意义是：单位水平面积承压含水层柱体，当水头下降一个单位时所释放的水量(无量纲)。 极坐标下均质、等厚、各向同性承压含水层轴对称流(径向流) 稳定流条件 当存在源汇项时 ? 和W分别为三维流和平面二维流的源汇。分别定义为单位体积含水层和单位水平面积含水层柱体中，单位时间内产生（为正值）或消耗（为负值）的水量。 总结各种形式，当存在源汇项时左端加上? 原形 均质 二度各向异性 轴对称问题 各向同性介质 稳定流条件 原形 均质 二度各向异性 轴对称问题 各向同性介质 稳定流条件 2.4 潜水流动的布西涅斯克微分方程 渗流的垂直分流速度vz远远小于水平分流速度vx和vy,可忽略vz,既假定等水头面是铅垂面。 一、裘布依假定(Dupuit Assumption): 裘布依微分方程 对于剖面二维流，其中H仅仅随x而变，与z无关，即H=H(x,t)。对于单宽流量qx,有 裘布依微分方程应用 裘布依假定的应用: 使剖面二维流问题(x,z)降阶为水平一维问题近似处理 使三维问题(x,y,z)降阶为水平二维(x,z)问题处理 使潜水面边界处理的简单化,直接近似地在微分方程中处理 裘布依假定的局限性: 当潜水面存在垂向补给、排泄或潜水呈不稳定流时，即使潜水面坡度很小，能否引出裘布依假定而将等水头面视为铅垂面则要视条件具体分析 二、布西涅斯克(Boussinesq)微分方程 假定: 研究的潜水流满足裘布依假定 水和骨架不可压缩（含义：无弹性储存）?=const, ?s=0 潜水层隔水底板水平 潜水面上存在水量的交换W W0， 入渗 W0， 蒸发 W定义为潜水面处单位水平面积、单位时间的入渗量 布西涅斯克方程推导 研究剖面二维流（x-z) 均衡单元体：长度为?x，宽度为1个单位的含水层柱体 均衡 ?V=?流入－?流出 均衡时段 ?t x方向流入流出差 z方向流入 上边界－－潜水面边界 下边界－－隔水边界 V=0 x,z方向流入流出差 单元体内水体积的增量表现为水位在该时段内的上升或下降 水均衡方程: 布西涅斯克方程推导 注意: ?d为重力给水度，是在重力作用下单位水平面积的潜水含水层柱体在其潜水面下降或上升单位水头时释放或储存的水量。 上述推导运用了裘布依假定，忽略垂向分流速，因此可近似用?x代替?s. 布西涅斯克方程将刻画潜水面边界问题简单化为用W直接表示。 布西涅斯克微分方程 线性化布西涅斯克方程 线性化布西涅斯克方程 线性化布西涅斯克方程 2.5 地下水流动定解条件及数学模型 地下水流动控制微分方程 潜水二维不稳定流动控制方程 承压水二维不稳定流动控制方程 定解条件 边界条件 初始条件 2.5 地下水流动定解条件及数学模型 1 数学模型的有关概念 同一形式的偏微分方程代表了整个一大类的地下水流的运动规律，而对于不同边界性质、不同边界形状的含水层，水头的分布是不同的。而且对于偏微分方程而言，方程本身并不包含反映特定渗流区条件的全部信息，方程可能存在无数个解，如需要从大量的可能解中求得与特定区域条件相对应的唯一特解，就必须提供反映特定区域特征的信息。 这些信息包括： （1）微分方程中的有关参数， 当这些参数确定后，微分方程才能被确定下来。 （2）渗流区范围和形状，当微分方程所对应的区域被确定之后才能对方程求解。 （3）边界条件(boundary condi