




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第二章
水动力学基本微分方程12/21/202212-1
地下水动力学基本微分方程12/21/20222一、含水层的弹性理论1.含水层的弹性释水(以承压水为例)
从承压含水层中抽出的水,由两部分组成含水层所贮存水的弹性释放侧向补给(来自远方)现取一处于平衡状态的承压含水层柱体。设含水层上覆岩层对含水层中的固体颗粒和地下水产生的应力为σ,骨架上的反压强为,水的顶托力为P。当水处于平衡状态时,当从含水层中抽水、水头下降△h时,图12/21/20223上述分析表明:H降低,承压含水层释放部分地下水;H增大,承压含水层贮存部分地下水,这部分水量称为弹性贮存量。
弹性贮水量的大小与含水层的岩性和结构有关,为了表征含水层弹性释水(储水)的能力,下面将给出弹性贮水率和贮水系数的概念。12/21/202252.含水层的贮水率和贮水系数1.贮水率(Specificstorativity)用表示
当含水层中的水头降低(或升高)一个单位时,单位时间内在单位体积含水层中,由于水的弹性膨胀(或压缩)及含水层的弹性压缩(或膨胀)释放(或贮存)的水量,称为贮水率,也称单位贮存量,量纲为[L-1]。12/21/20226
2.贮水系数(storage
coefficient)
m为含水层厚度,用于二维流计算。无量纲,大部分含水层介于10-5~10-3之间物理意义:在单位面积、厚度为m的含水层柱体中,当水头降低(或升高)一个单位时,单位时间内从含水层中释放(或贮存)的水量。12/21/20227
3.给水度
对潜水含水层而言,当水头下降时,引起两部分排水:①含水层下部饱水部分的弹性释水,其释水能力用表示;②上部潜水面下降部分引起的重力疏干排水,这部分给水能力用给水度
(Specificyield)表示;给水度的物理意义:当含水层中水头下降一个单位时,在单位体积含水层中,由重力疏干所排出的水量。12/21/20228二、承压含水层中渗流基本微分方程(一)方程的建立设:①地下水、含水层均为弹性体;②弹性释水瞬时完成;③为常数(constant);在渗流场中取一无限小的平行六面体,作为均衡单元,如图示,六面体边长分别为dx、dy、dz,下面分析dt时段内,均衡单元中的质量守恒问题。
均衡区均衡期12/21/202210dt时段内分别沿oy、oz进入单元的水量为而流出的水量为
在dt内,均衡单元贮存量的变化量为:
据水均衡原理得:
12/21/202212
上式就是非均质各向异性承压含水层中地下水三维非稳定运动的基本微分方程。对各向异性介质,取坐标轴方向与主渗透方向一致。(二)方程的化简和讨论
1.对于均质各向同性含水层,K为常数,这时简化为:12/21/2022142.对于二维的情况,常用和T表示(各项均乘以m)
3.当含水层有垂直水量交换时,其量常用W表示,称为源(汇)项,含水层的源(汇)项可是t和位置的函数W=W(x,y,z,t)。
当从含水层中抽水或从垂直方向有水流出含水层时,W为负,称为汇;当给含水层中注水或从垂直方向有水流入含水层时,W为正,称为源;12/21/202215
二维:W为单位时间从单位体积含水层中流入或流出的水量
三维:W为单位时间从垂直方向单位面积含水层中流入或流出的水量(补给强度或蒸发强度)
12/21/20221612/21/202217
可得到三维、二维相应的稳定状态下的渗流基本微分方程。对于均质各向同性的三维流来说:
5.若化为柱坐标(三维各向同性介质)12/21/20221812/21/202220(二)越流含水层中渗流基本微分方程
1.假定a.忽略弱透水层的弹性释水;b.水流在弱透水层中是垂向运动,而在主含水层中折射为水平运动;
2.方程的建立
在主含水层中取一微分柱体(其长宽分别为dx、dy,高为含水层厚度m)作为均衡单元。下面分析在dt时段内,微分柱体的水均衡问题。12/21/202221沿x方向流入单元体的水量:
流出:
沿y方向流入单元体的水量:
流出:
沿z方向流入单元体的水量:
流出:
12/21/202223流入量-流出量=:单元体内贮存量的变化为:根据水均衡原理得:12/21/202224(三)越流因素和越流系数1.越流因素(Leakagefactor)
上式中,若介质为均质各向同性介质,T=Constant式中:T为主含水层的导水系数Ki,mi分别为弱透水层的渗透系数和厚度12/21/202226B1、B2分别称为上、下两弱透水层的越流因素
越流因素B是反映弱透水层隔水性能的参数。B越大,越流量Q越小,对于隔水层B无穷大,越流量为零。
12/21/2022272.越流系数(Coefficientofleakage)
因为物理意义:当主含水层和相邻含水层间的水头差等于一个长度单位时,通过单位面积含水层上的越流量。越流系数反映越流量的大小,越大,相同水头下的越流量也越大。
12/21/202228
相当于忽略了渗透速度的垂直分量,代替,在铅垂面上各点的水头都是相等的;或者说,水头不随深度而变化,同一铅直面上各点的水力坡度和渗透速度都相等,渗透速度可表示为:
在此条件下,通过宽度为B的铅直平面沿x方向的流量Qx为:
h-含水层厚度,当隔水底板水平时,h=H12/21/202230(二)潜水含水层中渗流基本微分方程1.方程的建立
在Dupuit假定下,考虑一维问题,取平行于xoz平面的单位宽度进行研究。
首先取一微分柱体(其长度为∆x,宽为1,高为整个含水层厚度)作为均衡单元,下面分析在dt时段内,微分柱体的水均衡问题。12/21/202231从上游断面流入:从下游断面流出:
在∆t时间内,垂直方向的补给量为:
由于潜水面的上升而引起的均衡区内的水的增量为:
12/21/202232根据连续性原理,上面两个增量应相等,即代入上式得到
上式为非均质各向异性潜水含水层中一维地下水非稳定运动的基本微分方程,也称为J.Boussinesq方程。12/21/202233对于均质含水层,K=Constant
式中::水位下降时称为给水度,水位上升时称为饱和差;W:降雨入渗强度(+)或蒸发强度(—);h:含水层厚度;H:含水层的水位(平均值);K:含水层渗透系数;注意:
a.H为整个含水层厚度上的平均值;b.H、h均为未知,所以该方程为二阶非线性偏微分方程;c.该方程不适于水力梯度较大地段;不能计算任一点的H。12/21/202234
2.方程讨论
a.二维情况:b.均质含水层,K=Constantc.当隔水层底板水平时,H=h;12/21/202235d.若需计算渗流区内任一点水位,如土坝渗流等,则需采用下式计算:
其中潜水面要用边界条件考虑,不包括在方程之内。e.无蒸发、无入渗时,潜水位H变化很小,则
为稳定流。12/21/20223612/21/202237
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 四年级信息技术上册 回收站(一)教学设计 教科版
- 全国电子工业版初中信息技术第二册第2单元2.3活动1《认识网络新媒体》教学设计
- 北师大版数学二年级下册《测量》单元整体教学设计
- 2025至2030年中国低边拼密铝扣板数据监测研究报告
- 内江2025年四川内江市东兴区招募特聘动物防疫专员12人笔试历年参考题库附带答案详解
- 八上第二单元大单元教学设计
- Module5 Unit2 They were young(教学设计)2023-2024学年外研版(三起)英语四年级下册
- 佛山2025年广东佛山市技师学院招聘坐班制聘用人员(第五批)笔试历年参考题库附带答案详解
- 电子栅栏施工方案
- 人教版八年级下册历史与社会教案
- 2022年中小学体育课堂教学规范
- 新人教版八年级下册英语全册教案(教学设计)
- 2022年河南省郑州市中考二模语文试卷
- 东莞市卫生与健康十三五规划
- 土壤分析技术规范(第二版)
- 3力浮力答案第1讲难题型密度计
- 地下车库交通标志标线及地坪漆工程施工组织设计
- 专题一电磁感应与电路ppt课件
- GDFJ005修改个人信息申请表
- JJF 1363-2019硫化氢气体检测仪型式评价大纲(高清版)
- 氟喹诺酮类抗菌药物的不良反应和临床应用概要
评论
0/150
提交评论