水动力学基本微分方程课件

水动力学基本微分方程水动力学基本微分方程2-1

地下水动力学基本微分方程Date2Date2一、含水层的弹性理论1.含水层的弹性释水（以承压水为例）

从承压含水层中抽出的水，由两部分组成含水层所贮存水的弹性释放侧向补给（来自远方）现取一处于平衡状态的承压含水层柱体。设含水层上覆岩层对含水层中的固体颗粒和地下水产生的应力为σ，骨架上的反压强为，水的顶托力为P。当水处于平衡状态时，当从含水层中抽水、水头下降△h时，图Date3一、含水层的弹性理论1.含水层的弹性释水（以承压水为例）图D（1）H下降△h，水体积膨胀，从而释放出一定体积的水；（2）σ保持不变，骨架所受压力增加，因为固体颗粒接近于刚性体不可压缩，所以压力增加引起含水层压缩，使含水层空隙中的部分地下水被挤出。这两点就是弹性释水的机理。

Date4（1）H下降△h，水体积膨胀，从而释放出一定体积的水；Dat上述分析表明：H降低，承压含水层释放部分地下水；H增大，承压含水层贮存部分地下水，这部分水量称为弹性贮存量。

弹性贮水量的大小与含水层的岩性和结构有关，为了表征含水层弹性释水（储水）的能力，下面将给出弹性贮水率和贮水系数的概念。Date5上述分析表明：H降低，承压含水层释放部分地下水；H2.含水层的贮水率和贮水系数１.贮水率（Specificstorativity）用表示

当含水层中的水头降低（或升高）一个单位时，单位时间内在单位体积含水层中，由于水的弹性膨胀（或压缩）及含水层的弹性压缩（或膨胀）释放（或贮存）的水量，称为贮水率,也称单位贮存量，量纲为［L-1］。Date62.含水层的贮水率和贮水系数１.贮水率（Specifics

２.贮水系数（storagecoefficient）

m为含水层厚度，用于二维流计算。无量纲，大部分含水层介于10-5～10-3之间物理意义：在单位面积、厚度为m的含水层柱体中，当水头降低（或升高）一个单位时，单位时间内从含水层中释放（或贮存）的水量。Date7２.贮水系数（storagecoefficient）

3.给水度

对潜水含水层而言，当水头下降时，引起两部分排水：①含水层下部饱水部分的弹性释水，其释水能力用表示；②上部潜水面下降部分引起的重力疏干排水，这部分给水能力用给水度

（Specificyield）表示；给水度的物理意义：当含水层中水头下降一个单位时，在单位体积含水层中，由重力疏干所排出的水量。Date83.给水度Date84.贮水率与给水度的区别①弹性释水由减压引起，为压力变化所给出的水量，为重力疏干排出的水量；②贮水率与整个含水层厚度上的岩性、液体性质有关，给水度仅与水位波动带的岩性、液体性质有关；③弹性释放瞬时完成；重力疏干具明显的滞后效应；④数量级：约10-5～10-3；约0.1～0.3；∵＞＞∴在潜水含水层中，通常只考虑重力疏干，忽略弹性释水。Date94.贮水率与给水度的区别Date9二、承压含水层中渗流基本微分方程（一）方程的建立设：①地下水、含水层均为弹性体；②弹性释水瞬时完成；③为常数（constant）；在渗流场中取一无限小的平行六面体，作为均衡单元，如图示，六面体边长分别为dx、dy、dz，下面分析dt时段内，均衡单元中的质量守恒问题。

均衡区均衡期Date10二、承压含水层中渗流基本微分方程均衡区均衡期Date10

依据质量守恒和能量守恒定律，建立承压含水层中渗流基本微分方程。

Ｑx为单位时间内通过abcd断面流入的水量。在dt内，沿x方向通过abcd断面流入均衡单元的水量,通过a＇b＇c＇d＇断面从均衡单元流出的水量为Date11依据质量守恒和能量守恒定律，建立承压含水层中渗流基dt时段内分别沿oy、oz进入单元的水量为而流出的水量为

在dt内，均衡单元贮存量的变化量为:

据水均衡原理得：

Date12dt时段内分别沿oy、在dt内，均衡单元贮存量的变化化简为：根据达西定律：Date13化简为：根据达西定律：Date13

上式就是非均质各向异性承压含水层中地下水三维非稳定运动的基本微分方程。对各向异性介质，取坐标轴方向与主渗透方向一致。（二）方程的化简和讨论

１.对于均质各向同性含水层，Ｋ为常数，这时简化为：Date14上式就是非均质各向异性承压含水层中地下水２.对于二维的情况，常用和T表示(各项均乘以m）

３.当含水层有垂直水量交换时，其量常用Ｗ表示，称为源（汇）项，含水层的源（汇）项可是t和位置的函数W=W(x,y,z,t)。

当从含水层中抽水或从垂直方向有水流出含水层时，Ｗ为负，称为汇；当给含水层中注水或从垂直方向有水流入含水层时，Ｗ为正，称为源；Date15２.对于二维的情况，常用和T表示(各项均乘以m）３

二维：W为单位时间从单位体积含水层中流入或流出的水量

三维：W为单位时间从垂直方向单位面积含水层中流入或流出的水量(补给强度或蒸发强度)

Date16二维：W为单位时间从单位体积含水层中流入或流出的水量三维Date17Date17

可得到三维、二维相应的稳定状态下的渗流基本微分方程。对于均质各向同性的三维流来说：

5.若化为柱坐标（三维各向同性介质）Date18可得到三维、二维相应的稳定状态下的渗流基本三、越流含水层中渗流基本微分方程（一）什么是越流？

1.半承压含水层：一个主含水层的上层和（或）下层为弱透水层，主含水层通过弱透水层与相邻含水层发生水力联系，但它本身具有承压性，主含水层称为半承压含水层。2.越流：当半承压含水层与相邻含水层间存在水头差时，地下水便通过弱透水层从高水头含水层向低水头含水层产生垂向流动，这种现象称为越流。简言之，相邻含水层在水头差的作用下，通过弱透水层与主含水层发生水力联系的现象称为越流。

越流的方向：由两相邻含水层的水位决定。Date19三、越流含水层中渗流基本微分方程（一）什么是越流？Date1Date20Date20（二）越流含水层中渗流基本微分方程

1.假定a.忽略弱透水层的弹性释水；b.水流在弱透水层中是垂向运动，而在主含水层中折射为水平运动；

2.方程的建立

在主含水层中取一微分柱体（其长宽分别为dx、dy，高为含水层厚度m）作为均衡单元。下面分析在dt时段内，微分柱体的水均衡问题。Date21（二）越流含水层中渗流基本微分方程Date21P(x，y）Date22P(x，y）Date22沿x方向流入单元体的水量：

流出：

沿y方向流入单元体的水量：

流出：

沿z方向流入单元体的水量：

流出：

Date23沿x方向流入单元体的水量：沿y方向流入单元体的水量：沿z方向流入量－流出量＝：单元体内贮存量的变化为：根据水均衡原理得：Date24流入量－流出量＝：单元体内贮存量的变化为：Date24代入上式得：这就是不考虑弱透水层弹性释水条件下，非均质各向异性越流含水层系统地下水非稳定运动的基本微分方程。Date25代入上式得：这就是不考虑弱透水层弹性释水条件下，非均质（三）越流因素和越流系数１.越流因素（Leakagefactor）

上式中，若介质为均质各向同性介质，T＝Constant式中：Ｔ为主含水层的导水系数Ｋi，mi分别为弱透水层的渗透系数和厚度Date26（三）越流因素和越流系数式中：Ｔ为主含水层的导水系数DatＢ1、Ｂ2分别称为上、下两弱透水层的越流因素

越流因素Ｂ是反映弱透水层隔水性能的参数。Ｂ越大，越流量Ｑ越小，对于隔水层Ｂ无穷大，越流量为零。

Date27Ｂ1、Ｂ2分别称为上、下两弱透水层的越流因素越流因素Ｂ２.越流系数（Coefficientofleakage）

因为物理意义：当主含水层和相邻含水层间的水头差等于一个长度单位时，通过单位面积含水层上的越流量。越流系数反映越流量的大小，越大，相同水头下的越流量也越大。

Date28２.越流系数（Coefficientofleakage）四、潜水含水层中渗流基本微分方程（一）Dupuit假定

潜水面通常不是水平面，潜水含水层中存在着流速的垂直分量，潜水面本身又是渗流区边界，随时间而变化。为了建立潜水含水层中渗流基本微分方程，引出了Dupuit假定：假设潜水面比较平缓，潜水面上任意一点P有：Date29四、潜水含水层中渗流基本微分方程（一）Dupuit假定假设潜

相当于忽略了渗透速度的垂直分量，代替，在铅垂面上各点的水头都是相等的；或者说，水头不随深度而变化，同一铅直面上各点的水力坡度和渗透速度都相等，渗透速度可表示为：

在此条件下，通过宽度为B的铅直平面沿x方向的流量Ｑx为：

h-含水层厚度，当隔水底板水平时，h=HDate30相当于忽略了渗透速度的垂直分量（二）潜水含水层中渗流基本微分方程1.方程的建立

在Dupuit假定下，考虑一维问题，取平行于xoz平面的单位宽度进行研究。

首先取一微分柱体（其长度为∆x，宽为１，高为整个含水层厚度）作为均衡单元，下面分析在dt时段内，微分柱体的水均衡问题。Date31（二）潜水含水层中渗流基本微分方程首先取一微分柱从上游断面流入：从下游断面流出：

在∆t时间内，垂直方向的补给量为：

由于潜水面的上升而引起的均衡区内的水的增量为：

Date32从上游断面流入：在∆t时间内，垂直方向的由于潜根据连续性原理，上面两个增量应相等，即代入上式得到

上式为非均质各向异性潜水含水层中一维地下水非稳定运动的基本微分方程，也称为Ｊ.Boussinesq方程。Date33根据连续性原理，上面两个增量应相等，即代入上式得到上对于均质含水层，Ｋ＝Constant

式中：：水位下降时称为给水度，水位上升时称为饱和差；W：降雨入渗强度（+）或蒸发强度（—）；h：含水层厚度；H：含水层的水位（平均值）；K：含水层渗透系数；注意：

a.Ｈ为整个含水层厚度上的平均值；b.Ｈ、h均为未知，所以该方程为二阶非线性偏微分方程；c.该方程不适于水力梯度较大地段；不能计算任一点的Ｈ。Date34对于均质含水层，Ｋ＝Constant式中：：水位下降时称

2.方程讨论

a.二维情况：b.均质含水层，Ｋ＝Constantc.当隔水层底板水平时，H＝h；Date352.方程讨论a.