地下水运动基本定律、基本微分方程与数学模型.ppt

第二章 地下水运动基本定律、 基本微分方程和数学模型;介质：地下水赋存于岩土的空隙中，并在其中运动。我们将赋存地下水的岩土称为介质。 渗透：地下水在岩土空隙中的运动称为渗透。;第一节 地下水运动的基本定律:达西定律;第一节 地下水运动的基本定律:达西定律;假设水流必须符合下列条件：;第一节 达西定律;二、达西定律;达西定律，又称线性渗透定律;用微分来表示，即：;达西定律的实质是水流在流动过程中消耗的能量与流速和渗流长度成正比，与含水层的渗透系数成反比。;达西定律的适用范围;渗透系数K 水力坡度为1时的渗透流速。 ;意义：潜水位上升(下降)一个单位时,从单位面积含水层增加(减少)的水量。;第二节 数学模型;第二节 数学模型;第二节 渗流的基本微分方程和数学模型;抽水后，水头降低ΔH，此时上覆岩层总压力不变，为保持平衡，水头降低减少的压力与骨架增加的压力相等， ; 对于含水层骨架有：;当水头降低总共可得到的弹性水量为：;意义：水平面为一个单位面积，高为含水层全厚度M 的含水层柱体中，当水头降低一个单位时弹性释放出来的水量。;意义：潜水位上升(下降)一个单位时,从单位面积含水层增加(减少)的水量。;导水系数T 当水力坡度为1时，通过整个含水层上的单位宽度流量。即： ;水的状态方程 对于给定质量的水体积，增加一个压力dPw，水体积产生一定的压缩，根据质量守恒定律： ρ?Vw=常数 取全微分有： ρdVw+Vwdρ=0 由于dPw=?dH ;水的状态方程;颗粒骨架的状态方程 颗粒骨架在压力的作用下主要表现为垂直方向的空隙变形，而颗粒本身体积基本不变，因此，颗粒的体积为一常量： Vs-nVs=常量 取全微分 dVs-ndVs-Vsdn=0 ; 假定含水层骨架仅在垂直方向变形有：;（2）渗流连续性方程 根据对渗流的假说，渗流场全部空间都被连续水流充满。在渗流场中取任一微小单元体，其坐标为（X、Y、Z），边长分别为Δx、Δy、Δz，地下水的密度为ρ，在X、Y、Z方向上的渗流速度为u、v、w。 在Δt时间???，沿X方向，流入单元的水量为： ρQxΔt=ρuΔyΔzΔt;第二节 数学模型; 而沿X方向流出的水量为：;沿Y方向增加的水量; 因此，在Δt时间内，单元体增加的水量为：; 根据达西定律，水流流速在X、Y、Z方向有;地下水渗流连续性方程 表示：在渗流场中的任何局部，都必须满足质量守恒和能量守恒。; 对于稳定渗流，且假定n、ρ不变，则为地下水稳定流的连续性方程：;（3）承压含水层的基本微分方程 在承压含水层中，含水层产生变形时，主要是在垂直方向（Δz）上，而Δx、Δy近似为不变，因此，连续性方程为;第二节 数学模型;承压含水层三维非稳定渗流的微分方程; 若承压含水层水平等厚，渗透水流作水平二维流，则有：; 对均质各向同性的承压水作二维非稳定流时;对于有注水或抽水时，ε表示单位时间对单位面积含水层抽出或注入的水量，则;（4）潜水含水层的基本微分方程 潜水含水层的顶部是潜水面，在非稳定渗流的过程中潜水面的位置在不断变化中，因此要精确处理这类问题比较困难，但实际问题中，地下水流接近水平流动，潜水面也比较平缓.因此，可作出一定假定条件来近似求解，即裘布衣假定。 ;潜水面比较平缓，为缓变流； 渗流速度的水平分量u、v沿z高度没有变化，仅为x、y坐标和时间t的函数，即垂直流速可忽略不计或水头不随深度变化； 过水断面近似为一个垂直平面。;;dt时间内，从上游断面流入的水量：;Δx土体水量增量是：;水量的变化会引起潜水面的上升或下降，在dt时间内潜水面变化：;水量平衡：;潜水二维流的微分方程：;当水头变化很小时，即ΔH0.1h时，对均质各向同性的潜水有;若隔水底板水平，以隔水底板为基准面，此时有，H=h，当无其他补给和排泄时，均质各向同性的潜水二维稳定流方程为：;（5）定解条件和数学模型 初始条件 开始时刻水头函数在渗流场的分布规律，如果开始时刻渗流场内任意一点的水头已知为H0，则初始条件的数学表达式为：; 注：对于稳定流来说，定解条件中没有初始条件，因为地下水作稳定流时其运动要素是不随时间而变化的。 边界条件 是指在各个计算时刻，边界上某些函数的变化规律是已知的。边界条件主要有两种：; ①第一类边界条件（已知水头边界，Dirichlet条件）