《地下水利用》讲义.ppt
《《地下水利用》讲义.ppt》由会员分享,可在线阅读,更多相关《《地下水利用》讲义.ppt(122页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、LOGO地下水利用地下水利用夏龙兴15838106986水利与环境学院:市政教研室修改日期:2011.02LOGO本课程主要讲解内容本课程主要讲解内容第一章第一章 管管 井井 出出 水水 量量 计计 算算1第二章第二章 地下水资源计算与评价地下水资源计算与评价2第第三三章章 地下水地下水资源计算的数值法资源计算的数值法3第第四四章章 地下水开采建筑物设计地下水开采建筑物设计45 556第六章第六章 地下水资源管理地下水资源管理第第五五章章 井井 灌灌 区区工程工程 规规 划划LOGO 第一章第一章 管井出水量计算管井出水量计算 v1.开采地下水井的分类和地下水种类开采地下水井的分类和地下水种类
2、v2.单井出水量的稳定流计算单井出水量的稳定流计算(裘布依模型)(裘布依模型)v3.单井出水量的单井出水量的非非稳定流计算稳定流计算(泰斯模型)(泰斯模型)LOGO 第一第一节节.井的分类和地下水种类井的分类和地下水种类v 根据井径大小和构造的不同分为管井根据井径大小和构造的不同分为管井、筒井与大口井筒井与大口井、辐射井。辐射井。v1 1、管井:、管井:直径较小(直径较小(250250450450)mm mm,深度较大,深度较大(5050800)m,800)m,井井壁采用钢管,铸铁管,钢筋混凝土管或塑料管加固的井型。壁采用钢管,铸铁管,钢筋混凝土管或塑料管加固的井型。v 一般采用钻机施工,水泵
3、机组抽水,又称作机井。城镇与工一般采用钻机施工,水泵机组抽水,又称作机井。城镇与工业区用水机井多开采深层承压水;农业机井一般开采浅层的潜水业区用水机井多开采深层承压水;农业机井一般开采浅层的潜水含水层。含水层。v2 2、筒井与大口井、筒井与大口井v筒井:直径较大(筒井:直径较大(1.01.01.51.5)m m,深度较浅,深度较浅(30(3050)m,50)m,井壁采用井壁采用钢筋混凝土管或砖加固的井型钢筋混凝土管或砖加固的井型。v大口井:直径大(大口井:直径大(1.51.51010)m m,深度浅,深度浅(10(1030)m,30)m,井壁采用钢筋井壁采用钢筋混凝土管或砖加固的井型。混凝土管
4、或砖加固的井型。一一般开采埋深较小且含水层厚度小般开采埋深较小且含水层厚度小的潜水含水层的潜水含水层,安装水泵机组抽水,安装水泵机组抽水。LOGO 第一第一节节.井的分类和地下水种类井的分类和地下水种类v 3 3、辐射井、辐射井v 由垂直的集水井和沿集水井四周水平辐射状的进水管组成。由垂直的集水井和沿集水井四周水平辐射状的进水管组成。集水井直径(集水井直径(1.51.51010)m m,较浅,较浅(5(550)m50)m。适用于开采埋深小,。适用于开采埋深小,含水层薄,富水性差的地区,如黄土区,砾卵石含水层,河漫滩含水层薄,富水性差的地区,如黄土区,砾卵石含水层,河漫滩地。地。v 根据开采地下
5、含水层的不同分为:潜水完整井,潜水非完整根据开采地下含水层的不同分为:潜水完整井,潜水非完整井,承压完整井,承压非完整井。井,承压完整井,承压非完整井。v(1 1)潜水含水层:自由浸润面到不透水地板之间的饱含重力)潜水含水层:自由浸润面到不透水地板之间的饱含重力v 水的地质构造层。水的地质构造层。v(2 2)承压含水层:上、下不透水地板之间的饱含承压水的)承压含水层:上、下不透水地板之间的饱含承压水的v 地质构造层。地质构造层。v(3 3)完整井与非完整井:穿过整个含水层直达隔水层的井)完整井与非完整井:穿过整个含水层直达隔水层的井v 为完整井;否则为非完整井为完整井;否则为非完整井。v LO
6、GO 第第二二节节 单井出水量的稳定流计算单井出水量的稳定流计算v 一、潜水完整井出水量一、潜水完整井出水量计计算裘布依公式算裘布依公式v 1 1、裘布依模型、裘布依模型v(1 1)潜水含水)潜水含水层层特征:特征:圆岛圆岛柱体,均柱体,均质质同性,等厚,同性,等厚,v 地板水平地板水平v(2 2)边边界条件界条件:圆岛圆岛四周定水四周定水头头供水供水边边界界v(3 3)渗流特征:符合达西定律和)渗流特征:符合达西定律和连续连续方程;方程;轴对轴对称称v 均匀无阻力均匀无阻力汇汇入井中入井中v(4 4)稳稳定井流:抽水形成降落漏斗,当降落漏斗定井流:抽水形成降落漏斗,当降落漏斗扩扩v 展到供水
7、展到供水边边界界时时,抽水流量与,抽水流量与边边界界v 供水流量相等,降落漏斗也供水流量相等,降落漏斗也稳稳定。定。v LOGO 第第二二节节 单井出水量的稳定流计算单井出水量的稳定流计算v 2 2、公式推导、公式推导v由达西定律由达西定律 Qr Qr2rh2rhK Kdh/dr dh/dr v 流量连续方程流量连续方程QrQrQ Qv 积分得积分得 h h QLnr/KQLnr/KC C ,当,当r rR R时时,h,hH Hv C CH H QLnR/KQLnR/Kv 不同断面处不同断面处 Q QK(HK(H h h)/Ln(R/r)/Ln(R/r)v 井出流量井出流量 Q QK(HK(H
8、 h h0 0)/Ln(R/r)/Ln(R/r0 0)v Q Q1.364K(H1.364K(H h h0 0)/Lg(R/r)/Lg(R/r0 0)v二、承压完整井出水量计算裘布依公式二、承压完整井出水量计算裘布依公式 v 裘布依模型条件与潜水完整井裘布依模型相同。裘布依模型条件与潜水完整井裘布依模型相同。v Q Q2.73KM(H2.73KM(Hh h0 0)/Lg(R/r)/Lg(R/r0 0)v 式中式中 M M含水层厚度含水层厚度(m)(m),H H静水头(静水头(m)m)v 作业:公式推导作业:公式推导 三、裘布依公式讨论三、裘布依公式讨论 1 1、实际应用、实际应用 Q1.364
9、K(Hh0)/Lg(R/r0)已知含水层已知含水层K K,H H或或M M,R R,根据井的设计降深,根据井的设计降深S S0 0(S0=Hh0)可求设计流量;根据设计流量可求设计流量;根据设计流量Q Q推求水位降深推求水位降深S S0 0;根据井的;根据井的 稳定流抽水试验资料求水文地质参数稳定流抽水试验资料求水文地质参数K K。2 2、井出水量与井中水位降深的关系:、井出水量与井中水位降深的关系:承压井承压井Q=q SQ=q S0 0 ,直线关系;潜水井二者则呈二次线关系。,直线关系;潜水井二者则呈二次线关系。3 3、实测降深较计算值大、实测降深较计算值大 计算值忽略下列水力损失:(计算值
10、忽略下列水力损失:(1 1)井壁泥浆堵塞)井壁泥浆堵塞 (2 2)花管水力损失)花管水力损失 (3 3)向上转弯水头损失)向上转弯水头损失 (4 4)沿吸水管至泵进口的沿程水头损失,)沿吸水管至泵进口的沿程水头损失,因此因此,承压完整井出水量与降深并非呈直线关系。承压完整井出水量与降深并非呈直线关系。4 4、水跃、水跃 :水流流入井壁时存在水位比井中水位高。:水流流入井壁时存在水位比井中水位高。井壁附近水流流线为曲线,存在水跃水才能流入井中。井壁附近水流流线为曲线,存在水跃水才能流入井中。LOGO 第第二二节节 单井出水量的稳定流计算单井出水量的稳定流计算v 5 5、影响半径、影响半径R R0
11、 0v 实际含水层,抽水一定时间后,由于抽水时间有限,实际含水层,抽水一定时间后,由于抽水时间有限,存在补水条件,当抽水量与补给量达到平衡时,地下渗存在补水条件,当抽水量与补给量达到平衡时,地下渗流为稳定流;另外,抽水影响范围有限,从实际观测不流为稳定流;另外,抽水影响范围有限,从实际观测不到地下水位下降处到井中心的水平距离确定为稳定抽水到地下水位下降处到井中心的水平距离确定为稳定抽水影响半径影响半径R R0 0,R R0 0处流量为零,此处不存在水力坡降。处流量为零,此处不存在水力坡降。v四、群井四、群井干扰抽水干扰抽水出水量稳定流计算出水量稳定流计算v 已知已知两口相同两口相同承压完整井承
12、压完整井,间距为,间距为2b2b,1 1号井单独号井单独抽水时的流量抽水时的流量Q Q和降深和降深S S0 0,引起,引起2 2号井中水位降深号井中水位降深t t;2 2号号井单独抽水时的流量井单独抽水时的流量Q Q和降深和降深S S0 0,引起,引起1 1号井中水位降深号井中水位降深t t。两口井同时干扰(井距近)抽水。两口井同时干扰(井距近)抽水。v 水位消减法:水位消减法:(1 1)若)若1 1号井仍保持出水量号井仍保持出水量Q Q不变,实际降深不变,实际降深S SS S0 0t t Q Q2.73KMS2.73KMS0 0/Lg(R/r/Lg(R/r0 0)S S0 0QLg(R/rQ
13、Lg(R/r0 0)/2.73KM )/2.73KM t tQLg(R/2QLg(R/2b b)/2.73KM)/2.73KM 求得实际降深求得实际降深S SS S0 0t tQLg(RQLg(R/r/r0 02 2b b)/2.73KM )/2.73KM 两井两井同时干扰抽水总流量等于同时干扰抽水总流量等于2Q 2Q (2 2)若)若1 1号井保持降深号井保持降深S0不变,有效降深不变,有效降深S SS S0 0t t Q Q2.73KMS2.73KMS0 0/Lg(R/r/Lg(R/r0 0)S S0 0QLg(R/rQLg(R/r0 0)/2.73KM )/2.73KM t tQLg(R
14、/2QLg(R/2b b)/2.73KM)/2.73KM 有效降深有效降深S SS S0 0t tQLg(2QLg(2b b/r/r0 0)/2.73KM )/2.73KM 1 1号井实际出水量号井实际出水量 QQ QQ2.73KMS/Lg(R/r2.73KMS/Lg(R/r0 0)两井两井同时干扰抽水总流量等于同时干扰抽水总流量等于2QQ2QQ LOGO 第第二二节节 单井出水量的稳定流计算单井出水量的稳定流计算v(3 3)N N眼眼承压承压井同时干扰抽水(井距近)井同时干扰抽水(井距近)v 1 1号井:号井:S1S1S S0101(t2(t2t3-t3-t tn n)v 为为1 1号井保持
15、出水量号井保持出水量Q1Q1不变,实际降深。为不变,实际降深。为1 1号井保持降号井保持降深深S S0101不变,有效降深。不变,有效降深。t2t2t3t3-t-tn n 分别为其它井抽水在分别为其它井抽水在1 1号号井引起的降深。井引起的降深。v 2 2号井:号井:S2S2S S0202(t1(t1t3-t3-t)t)v(4 4)N N眼潜水井眼潜水井干扰抽水(井距近)干扰抽水(井距近)v1 1号井单独抽水时的流量和降深为号井单独抽水时的流量和降深为Q1Q1,S01S01v Q1Q11.364K1(H11.364K1(H1 h01h01)/LgR1/r01 )/LgR1/r01 v Q1 Q
16、11.364K11.364K12Hp12Hp1S01S01/LgR1/r01 /LgR1/r01 v S01S01Q1LgQ1Lg(R1/r01R1/r01)/2.73K1Hp1 /2.73K1Hp1 v在在2 2号井引起降深号井引起降深t2t2Q2Lg(R2/2B2)/2.73K2Hp2 Q2Lg(R2/2B2)/2.73K2Hp2 v S1S1S01S01(t2+t3-+t)(t2+t3-+t)v 为为1 1号井保持出水量号井保持出水量Q1Q1不变,实际降深。为不变,实际降深。为1 1号井保持降号井保持降深深S01S01不变,有效降深。不变,有效降深。v 2 2号井号井S2S2S02S02
17、(t1(t1t3-t3-t)t)LOGO 第第二二节节 单井出水量的稳定流计算单井出水量的稳定流计算v五五 单井出水量计算的抽水试验经验公式法单井出水量计算的抽水试验经验公式法 v 由于实际含水层与裘布依模型条件差别大,计算由于实际含水层与裘布依模型条件差别大,计算结果误差大,可采用经验公式法,根据抽水试验资料,结果误差大,可采用经验公式法,根据抽水试验资料,比较反应复杂的实际情况,建立流量与降深的关系曲比较反应复杂的实际情况,建立流量与降深的关系曲线。常用的曲线有:线。常用的曲线有:v11直线型直线型Q=qSQ=qS,回归系数,回归系数q q为水井单位出水量为水井单位出水量v m m(h.m
18、)(h.m)。v22抛物线型抛物线型S=aQ+bQS=aQ+bQ 回归系数回归系数a,ba,bv33指数曲线型指数曲线型 S=S=(Q/nQ/n)m m 回归系数回归系数n,mn,mv44对数曲线型对数曲线型 Q=a+blgs Q=a+blgs 回归系数回归系数a,ba,bv 具体公式形式可采用图解法或查分法或最小二乘具体公式形式可采用图解法或查分法或最小二乘法求解。一般试验数据处理采用最小二乘法求解(有法求解。一般试验数据处理采用最小二乘法求解(有软件)。软件)。LOGO 第第三三节节 单井出水量的单井出水量的非非稳定流计算稳定流计算一一 潜水完整井出水量非稳定流微分方程式潜水完整井出水量非
19、稳定流微分方程式(泰斯模型泰斯模型)(1)(1)潜水含水层特征:均质、同性、等厚、地板水平且无限延伸潜水含水层特征:均质、同性、等厚、地板水平且无限延伸。(2)(2)渗流特征:符合达西定律和连续方程;轴对称均匀无阻力汇渗流特征:符合达西定律和连续方程;轴对称均匀无阻力汇 入井中入井中。(3)(3)垂直方向无补给垂直方向无补给。求解模型:求解模型:取微元体:半径为取微元体:半径为r r及及r+dr,r+dr,高度为高度为H H及及H H+d+dH H的圆筒体的圆筒体1 1、由圆筒内壁流入井中的流量为、由圆筒内壁流入井中的流量为Q,Q,由圆筒外壁流入的流量为由圆筒外壁流入的流量为 Q Q Q/Q/
20、r r,则则流出与流入的流量差:流出与流入的流量差:Q Q(Q(Q Q/Q/r rdr)dr)Q/Q/r rdrdr2 2、重力释水量、重力释水量 2rdr2rdr H/H/t t u,uu,u潜水层给水度潜水层给水度LOGO 第第三三节节 单井出水量的单井出水量的非非稳定流计算稳定流计算 3 3、水量平衡:、水量平衡:Q/Q/r rdrdr2rdr2rdr H/H/t tu u 4 4、由达西、由达西公式流出内筒面的公式流出内筒面的渗流速度渗流速度V VK K H/H/r r Q QAVAV2rH2rHV V H H(Hr+H(r+dr)/2(Hr+H(r+dr)/2 5 5、Q/Q/r r
21、2kH2kH(H/H/r)/r)/r r 2kH(2kH(H/H/r rr r H/H/r)r)v 代入水量平衡式化简得代入水量平衡式化简得v KH(KH(H/H/r rr r H/H/r)r)H/H/t tr ru u 令令 T TKHKH,H/H/r rr r H/H/r r r ru/T u/T H/H/t t LOGO 第第三三节节 单井出水量的单井出水量的非非稳定流计算稳定流计算 令令a aT/u T/u,H/H/r rr r H/H/r r H/H/t/at/a 在抽水过程中若有降水入渗补给或蒸发排泄时,在水量平衡在抽水过程中若有降水入渗补给或蒸发排泄时,在水量平衡计算时计入,设在
22、单位时间在单位面积上的补给率或蒸发率为计算时计入,设在单位时间在单位面积上的补给率或蒸发率为,上式改写为:,上式改写为:a a(H/H/r rr r H/H/r r)+/u u H/H/t tu u-潜水含水层给水度:水头下降一个单位时,从单位含水层面积潜水含水层给水度:水头下降一个单位时,从单位含水层面积高度等于含水层厚度的柱体中由于重力作用释放出的水量。高度等于含水层厚度的柱体中由于重力作用释放出的水量。u u*-*-承压含水层弹性释水系数:当承压含水层的水头降低一个单位承压含水层弹性释水系数:当承压含水层的水头降低一个单位时,从单位含水层面积高度等于含水层厚度的柱体中释放出的水时,从单位
23、含水层面积高度等于含水层厚度的柱体中释放出的水量。量。T=KMT=KM为承压含水层的导水系数为承压含水层的导水系数a-a-在潜水含水层称为水位传导系数,在承压含水层称为压力传导,在潜水含水层称为水位传导系数,在承压含水层称为压力传导,表示含水层水位或水头变化的传导速度。表示含水层水位或水头变化的传导速度。二二 承压完整井出水量非稳定流计算泰斯公式承压完整井出水量非稳定流计算泰斯公式 泰斯泰斯物理物理模型模型:(1)(1)潜水含水层特征:均质、同性、等厚、地板水平潜水含水层特征:均质、同性、等厚、地板水平且无限延伸且无限延伸 (2)(2)渗流特征:符合达西定律和连续方程;轴对称均渗流特征:符合达
24、西定律和连续方程;轴对称均匀无阻力汇入井中匀无阻力汇入井中 (3)(3)垂直方向无补给垂直方向无补给(4)(4)定流量抽水,井径无限小定流量抽水,井径无限小 泰斯数学模型泰斯数学模型:H/H/r/r r/r H/H/r r H/H/t/a t/a (r0r0)r0r0)H(r,0)H(r,0)H0 H0 (初始条件初始条件r0r,tr0r,t0)0)H(,t)H(,t)H0 H0 (定水头条件定水头条件r0)r0)lim(r lim(r H/H/r)r)Q/2TQ/2T常数常数 r0 r0LOGO 第第三三节节 单井出水量的单井出水量的非非稳定流计算稳定流计算 解解:偏微分方程组偏微分方程组
25、v1 1、二元变量转化为一元变量(偏微分方程变为常微分方程)、二元变量转化为一元变量(偏微分方程变为常微分方程)v构造一个二元复合变量构造一个二元复合变量 U=r U=r u*/4Tt u*/4Tt v U/U/r rru*/2Tt ru*/2Tt v U/U/t tr r u*/4Ttu*/4Tt v H/H/r r U/U/r rdH/du=ru*/2Tt dH/du=ru*/2Tt dH/dudH/duv H/H/r r (H/H/r)/r)/r rv u*/2Ttu*/2TtdH/dudH/dur r u*u*/4T/4T t td d H/dH/d u u v H H/t tr r
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地下水利用 地下水 利用 讲义
限制150内