水动力学基本微分方程(共39张PPT).pptx

水动力学基本微分方程第一页，共三十九页。2-1 地下水动力学基本微分方程第二页，共三十九页。一、含水层的弹性理论1.含水层的弹性释水（以承压水为例） 从承压含水层中抽出的水，由两部分组成 含水层所贮存水的弹性释放 侧向补给（来自远方） 现取一处于平衡状态的承压含水层柱体。设含水层上覆岩层对含水层中的固体颗粒和地下水产生的应力为σ，骨架上的反压强为 ，水的顶托力为P。 当水处于平衡状态时， 当从含水层中抽水、水头下降△h时，图第三页，共三十九页。（1）H下降△h，水体积膨胀，从而释放出一定体积的水；（2）σ保持不变，骨架所受压力增加，因为固体颗粒接近于刚性体不可压缩，所以压力增加引起含水层压缩，使含水层空隙中的部分地下水被挤出。 这两点就是弹性释水的机理。 第四页，共三十九页。 上述分析表明：H降低，承压含水层释放部分地下水；H增大，承压含水层贮存部分地下水，这部分水量称为弹性贮存量。 弹性贮水量的大小与含水层的岩性和结构有关，为了表征含水层弹性释水（储水）的能力，下面将给出弹性贮水率和贮水系数的概念。第五页，共三十九页。１.贮水率（Specific storativity）用 表示 当含水层中的水头降低（或升高）一个单位时，单位时间内在单位体积含水层中，由于水的弹性膨胀（或压缩）及含水层的弹性压缩（或膨胀）释放（或贮存）的水量，称为贮水率,也称单位贮存量，量纲为［L-1］。第六页，共三十九页。　２.贮水系数（ storage coefficient ） m为含水层厚度，用于二维流计算。 无量纲，大部分含水层 介于10-5～10-3之间 物理意义：在单位面积、厚度为m的含水层柱体中，当水头降低（或升高）一个单位时，单位时间内从含水层中释放（或贮存）的水量。第七页，共三十九页。 3.给水度 对潜水含水层而言，当水头下降时，引起两部分排水： ①含水层下部饱水部分的弹性释水，其释水能力用 表示； ②上部潜水面下降部分引起的重力疏干排水，这部分给水能力用给水度 （Specific yield）表示； 给水度的物理意义：当含水层中水头下降一个单位时，在单位体积含水层中，由重力疏干所排出的水量。第八页，共三十九页。4.贮水率与给水度的区别① 弹性释水由减压引起， 为压力变化所给出的水量， 为重力疏干排出的水量；② 贮水率与整个含水层厚度上的岩性、液体性质有关，给水度仅与水位波动带的岩性、液体性质有关；③ 弹性释放瞬时完成；重力疏干具明显的滞后效应；④ 数量级： 约10-5～10-3； 约～；∵ ＞＞∴在潜水含水层中，通常只考虑重力疏干，忽略弹性释水。第九页，共三十九页。二、承压含水层中渗流基本微分方程（一）方程的建立设： ① 地下水、含水层均为弹性体； ② 弹性释水瞬时完成； ③ 为常数（constant）； 在渗流场中取一无限小的平行六面体，作为均衡单元，如图示，六面体边长分别为dx、dy、dz，下面分析dt时段内，均衡单元中的质量守恒问题。 均衡区均衡期第十页，共三十九页。 Ｑx为单位时间内通过abcd断面流入的水量。在dt内，沿x方向通过abcd断面流入均衡单元的水量,通过a＇b＇c＇d＇断面从均衡单元流出的水量为 依据质量守恒和能量守恒定律，建立承压含水层中渗流基本微分方程。 第十一页，共三十九页。 dt时段内分别沿oy、 oz进入单元的水量为 而流出的水量为 在dt内，均衡单元贮存量的变化量为: 据水均衡原理得： 第十二页，共三十九页。化简为：根据达西定律：第十三页，共三十九页。 上式就是非均质各向异性承压含水层中地下水三维非稳定运动的基本微分方程。对各向异性介质，取坐标轴方向与主渗透方向一致。（二）方程的化简和讨论 １.对于均质各向同性含水层，Ｋ为常数，这时 简化为：第十四页，共三十九页。２.对于二维的情况，常用 和T表示( 各项均乘以m） ３.当含水层有垂直水量交换时，其量常用Ｗ表示，称为源（汇）项，含水层的源（汇）项可是t和位置的函数 W=W(x,y,z,t)。 当从含水层中抽水或从垂直方向有水流出含水层时，Ｗ为负，称为汇； 当给含水层中注水或从垂直方向有水流入含水层时，Ｗ为正，称为源；第十五页，共三十九页。三维：W为单位时间从单位体积含水层中流入或流出的水量 二维：W为单位时间从垂直方向单位面积含水层中流入或流出的水量(补给强度或蒸发强度) 第十六页，共三十九页。第十七页，共三十九页。 可得到三维、二维相应的稳定状态下的渗流基本微分方程。 对于均质各向同性的三维流来说：　5.若化为柱坐标（三维各向同性介质）第十八页，共三十九页。三、越流含水层中渗流基本微分方程（一）什么是越流？ 1.半承压含水层：一个主含水层的上层和（或）下层为弱透水层，主含水层通过弱透水层与相邻