地下水利用讲义.ppt

LOGO地下水利用水利与环境学院：市政教研室LOGO本课程主要讲解内容本课程主要讲解内容第一章第一章 管管 井井 出出 水水 量量 计计 算算1第二章第二章 地下水资源计算与评价地下水资源计算与评价2第三章第三章 地下水资源计算的数值法地下水资源计算的数值法3第四章第四章 地下水开采建筑物设计地下水开采建筑物设计45 556第六章第六章 地下水资源管理地下水资源管理第五章第五章 井井 灌灌 区工程区工程 规规 划划LOGO 第一章第一章 管井出水量管井出水量计算算 v1.开采地下水井的分类和地下水种类开采地下水井的分类和地下水种类v2.单井出水量的稳定流计算（裘布依模型）单井出水量的稳定流计算（裘布依模型）v3.单井出水量的非稳定流计算（泰斯模型）单井出水量的非稳定流计算（泰斯模型）LOGO 第一节第一节.井的分类和地下水种类井的分类和地下水种类v 根据井径大小和构造的不同分根据井径大小和构造的不同分为管井、筒井与大口井、管井、筒井与大口井、辐射井射井。v1 1、管井：直径、管井：直径较小（小（250250450450）mm mm，深度，深度较大大(50(50800800)m,)m,井井壁采用壁采用钢管，管，铸铁管，管，钢筋混凝土管或塑料管加固的井型。筋混凝土管或塑料管加固的井型。v 一般采用一般采用钻机施工，水机施工，水泵机机组抽水，又称作机井。城抽水，又称作机井。城镇与工与工业区用水机井多开采深区用水机井多开采深层承承压水；水；农业机井一般开采浅机井一般开采浅层的潜水的潜水含水含水层。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,井壁采用井壁采用钢筋筋混凝土管或混凝土管或砖加固的井型。加固的井型。一般开采埋深一般开采埋深较小且含水小且含水层厚度小厚度小的潜水含水的潜水含水层，安装水，安装水泵机机组抽水。抽水。LOGO 第一节第一节.井的分类和地下水种类井的分类和地下水种类v 3 3、辐射井射井v 由垂直的集水井和沿集水井四周水平由垂直的集水井和沿集水井四周水平辐射状的射状的进水管水管组成。成。集水井直径（集水井直径（1.51.51010）m m，较浅浅(5(550)m50)m。适用于开采埋深小，。适用于开采埋深小，含水含水层薄，富水性差的地区，如黄土区，薄，富水性差的地区，如黄土区，砾卵石含水卵石含水层，河漫，河漫滩地。地。v 根据开采地下含水根据开采地下含水层的不同分的不同分为：潜水完整井，潜水非完整：潜水完整井，潜水非完整井，承井，承压完整井，承完整井，承压非完整井。非完整井。v（1 1）潜水含水）潜水含水层：自由浸：自由浸润面到不透水地板之面到不透水地板之间的的饱含重力含重力v 水的地水的地质构造构造层。v（2 2）承）承压含水含水层：上、下不透水地板之：上、下不透水地板之间的的饱含承含承压水的水的v 地地质构造构造层。v（3 3）完整井与非完整井：穿）完整井与非完整井：穿过整个含水整个含水层直达隔水直达隔水层的井的井v 为完整井；否完整井；否则为非完整井。非完整井。v LOGO 第二第二节 单井出水量的井出水量的稳定流定流计算算v 一、潜水完整井出水量一、潜水完整井出水量计算裘布依公式算裘布依公式v 1 1、裘布依模型、裘布依模型v（1 1）潜水含水）潜水含水层特征：特征：圆岛柱体，均柱体，均质同性，等厚，同性，等厚，v 地板水平地板水平v（2 2）边界条件界条件:圆岛四周定水四周定水头供水供水边界界v（3 3）渗流特征：符合达西定律和）渗流特征：符合达西定律和连续方程；方程；轴对称称v 均匀无阻力均匀无阻力汇入井中入井中v（4 4）稳定井流：抽水形成降落漏斗，当降落漏斗定井流：抽水形成降落漏斗，当降落漏斗扩v 展到供水展到供水边界界时，抽水流量与，抽水流量与边界界v 供水流量相等，降落漏斗也供水流量相等，降落漏斗也稳定。定。v LOGO 第二第二节 单井出水量的井出水量的稳定流定流计算算v 2 2、公式推、公式推导v由达西定律由达西定律 Qr Qr2rhKdh/dr 2rhKdh/dr v 流量流量连续方程方程QrQrQ Qv 积分得分得 h hQLnr/KQLnr/KC C ，当，当r rR R时,h,hH Hv C CHHQLnR/KQLnR/Kv 不同断面不同断面处 Q QK(HK(Hh)/Ln(R/r)h)/Ln(R/r)v 井出流量井出流量 Q QK(HK(Hh h0 0)/Ln(R/r)/Ln(R/r0 0)v Q Q1.364K(H1.364K(Hh 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.364K(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、实测降深较计算值大、实测降深较计算值大 计算值忽略下列水力损失：（计算值忽略下列水力损失：（1 1）井壁泥浆堵塞）井壁泥浆堵塞 （2 2）花管水力损失）花管水力损失 （3 3）向上转弯水头损失）向上转弯水头损失 （4 4）沿吸水管至泵进口的沿程水头损失，）沿吸水管至泵进口的沿程水头损失，因此，承压完整井出水量与降深并非呈直线关系。因此，承压完整井出水量与降深并非呈直线关系。4 4、水跃、水跃 ：水流流入井壁时存在水位比井中水位高。：水流流入井壁时存在水位比井中水位高。井壁附近水流流线为曲线，存在水跃水才能流入井中。井壁附近水流流线为曲线，存在水跃水才能流入井中。LOGO 第二第二节 单井出水量的井出水量的稳定流定流计算算v 5 5、影响半径、影响半径R R0 0v 实际含水层，抽水一定时间后，由于抽水时间有限，实际含水层，抽水一定时间后，由于抽水时间有限，存在补水条件，当抽水量与补给量达到平衡时，地下渗存在补水条件，当抽水量与补给量达到平衡时，地下渗流为稳定流；另外，抽水影响范围有限，从实际观测不流为稳定流；另外，抽水影响范围有限，从实际观测不到地下水位下降处到井中心的水平距离确定为稳定抽水到地下水位下降处到井中心的水平距离确定为稳定抽水影响半径影响半径R R0 0，R R0 0处流量为零，此处不存在水力坡降。处流量为零，此处不存在水力坡降。v四、群井干扰抽水出水量稳定流计算四、群井干扰抽水出水量稳定流计算v 已知两口相同已知两口相同承压完整井承压完整井，间距为，间距为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/rQLg(R/r0 0)/2.73KM )/2.73KM t tQLg(R/2b)/2.73KM QLg(R/2b)/2.73KM 求得实际降深求得实际降深S SS S0 0t tQLg(R/rQLg(R/r0 02b)/2.73KM 2b)/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/2b)/2.73KM QLg(R/2b)/2.73KM 有效降深有效降深S SS S0 0t tQLg(2b/rQLg(2b/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-tn)tn)v 为1 1号井保持出水量号井保持出水量Q1Q1不不变，实际降深。降深。为1 1号井保持降号井保持降深深S S0101不不变，有效降深。，有效降深。t2t3-tn t2t3-tn 分分别为其它井抽水在其它井抽水在1 1号号井引起的降深。井引起的降深。v 2 2号井：号井：S2S2S S0202(t1(t1t3-t3-t)t)v（4 4）N N眼潜水井眼潜水井干干扰抽水（井距近）抽水（井距近）v1 1号井号井单独抽水独抽水时的流量和降深的流量和降深为Q1Q1，S01S01v Q1 Q11.364K1(H11.364K1(H1h01)/LgR1/r01 h01)/LgR1/r01 v Q1 Q11.364K12Hp1S01/LgR1/r01 1.364K12Hp1S01/LgR1/r01 v S01 S01Q1LgQ1Lg（R1/r01R1/r01）/2.73K1Hp1 /2.73K1Hp1 v在在2 2号井引起降深号井引起降深t2t2Q2Lg(R2/2B2)/2.73K2Hp2 Q2Lg(R2/2B2)/2.73K2Hp2 v S1 S1S01(t2+t3-+t)S01(t2+t3-+t)v 为1 1号井保持出水量号井保持出水量Q1Q1不不变，实际降深。降深。为1 1号井保持降号井保持降深深S01S01不不变，有效降深。，有效降深。v 2 2号井号井S2S2S02(t1S02(t1t3-t3-t)t)LOGO 第二第二节 单井出水量的井出水量的稳定流定流计算算v五五 单井出水量计算的抽水试验经验公式法单井出水量计算的抽水试验经验公式法 v 由于实际含水层与裘布依模型条件差别大，计算由于实际含水层与裘布依模型条件差别大，计算结果误差大，可采用经验公式法，根据抽水试验资料，结果误差大，可采用经验公式法，根据抽水试验资料，比较反应复杂的实际情况，建立流量与降深的关系曲比较反应复杂的实际情况，建立流量与降深的关系曲线。常用的曲线有：线。常用的曲线有：v11直线型直线型Q=qSQ=qS，回归系数，回归系数q q为水井单位出水量为水井单位出水量v m m(h.m)(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 第三第三节 单井出水量的非井出水量的非稳定流定流计算算一一 潜水完整井出水量非潜水完整井出水量非稳定流微分方程式（泰斯模型）定流微分方程式（泰斯模型）(1)(1)潜水含水潜水含水层特征：均特征：均质、同性、等厚、地板水平且无限延伸。、同性、等厚、地板水平且无限延伸。(2)(2)渗流特征：符合达西定律和渗流特征：符合达西定律和连续方程；方程；轴对称均匀无阻力称均匀无阻力汇 入井中。入井中。(3)(3)垂直方向无垂直方向无补给。求解模型：求解模型：取微元体：半径取微元体：半径为r r及及r+dr,r+dr,高度高度为H H及及H+dHH+dH的的圆筒体筒体1 1、由、由圆筒内壁流入井中的流量筒内壁流入井中的流量为Q,Q,由由圆筒外壁流入的流量筒外壁流入的流量为 Q Q Q/Q/r r，则流出与流入的流量差：流出与流入的流量差：Q Q(Q(Q Q/Q/rdr)rdr)Q/Q/rdrrdr2 2、重力、重力释水量水量 2rdr2rdr H/H/t u,ut u,u潜水潜水层给水度水度LOGO 第三第三节 单井出水量的非井出水量的非稳定流定流计算算 3 3、水量平衡：、水量平衡：Q/Q/rdrrdr2rdr2rdr H/H/tu tu 4 4、由达西公式流出内筒面的渗流速度、由达西公式流出内筒面的渗流速度V VK K H/H/r r Q QAVAV2rHV2rHV H H(Hr+H(r+dr)/2(Hr+H(r+dr)/2 5 5、Q/Q/r r2kH2kH(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/trutru 令令 T TKHKH，H/H/r rr r H/H/rr ru/T ru/T H/H/t t LOGO 第三第三节 单井出水量的非井出水量的非稳定流定流计算算 令令a aT/u T/u，H/H/r rr r H/H/rr H/H/t/at/a 在抽水在抽水过程中若有降水入渗程中若有降水入渗补给或蒸或蒸发排泄排泄时，在水量平衡，在水量平衡计算算时计入，入，设在在单位位时间在在单位面位面积上的上的补给率或蒸率或蒸发率率为，上式改写，上式改写为：a a（H/H/r rr r H/H/rr）+/u u H/H/t tu-u-潜水含水潜水含水层给水度：水水度：水头下降一个下降一个单位位时，从，从单位含水位含水层面面积高度等于含水高度等于含水层厚度的柱体中由于重力作用厚度的柱体中由于重力作用释放出的水量。放出的水量。u*-u*-承承压含水含水层弹性性释水系数：当承水系数：当承压含水含水层的水的水头降低一个降低一个单位位时，从，从单位含水位含水层面面积高度等于含水高度等于含水层厚度的柱体中厚度的柱体中释放出的水放出的水量。量。T=KMT=KM为承承压含水含水层的的导水系数水系数a-a-在潜水含水在潜水含水层称称为水位水位传导系数，在承系数，在承压含水含水层称称为压力力传导，表示含水表示含水层水位或水水位或水头变化的化的传导速度。速度。二二 承压完整井出水量非稳定流计算泰斯公式承压完整井出水量非稳定流计算泰斯公式 泰斯物理模型泰斯物理模型:(1)(1)潜水含水层特征：均质、同性、等厚、地板水平潜水含水层特征：均质、同性、等厚、地板水平且无限延伸且无限延伸 (2)(2)渗流特征：符合达西定律和连续方程；轴对称均渗流特征：符合达西定律和连续方程；轴对称均匀无阻力汇入井中匀无阻力汇入井中 (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 第三第三节 单井出水量的非井出水量的非稳定流定流计算算 解解:偏微分方程偏微分方程组 v1 1、二元、二元变量量转化化为一元一元变量（偏微分方程量（偏微分方程变为常微分方程）常微分方程）v构造一个二元复合构造一个二元复合变量量 U=ru*/4Tt U=ru*/4Tt v U/U/r rru*/2Tt ru*/2Tt v U/U/t tru*/4Ttru*/4Ttv H/H/r r U/U/rdH/du=ru*/2Tt dH/durdH/du=ru*/2Tt dH/duv H/H/r r (H/H/r)/r)/r rv u*/2TtdH/duu*/2TtdH/duru*/4TtdH/du ru*/4TtdH/du v H/H/t tru*/4Tt*dH/duru*/4Tt*dH/duv代入模型方程得到：代入模型方程得到：v u dH/du+(1+u)dH/du u dH/du+(1+u)dH/du0 0 LOGO 第三第三节 单井出水量的非井出水量的非稳定流定流计算算 2 2、二、二阶常微分方程先降常微分方程先降阶v 令令dH/dudH/duG G 则 dG/du dG/du(1(1u)/u*Gu)/u*G0 0v 分离分离变量量积分：分：dG/GdG/G(1(11/u)du 1/u)du v lnG lnGu ulnulnuC1 C1 v lnGu lnGuu+C1 u+C1，v v G=*c/uG=*c/uv v 分离分离变量量积分：当分：当 r r时 ，u,HH0u,HH0v v dH dHC /uduC /udu H=H H=H0 0-C /udu-C /udu 3 3、d/du(H0d/du(H0C /udu)C /udu)d/du(-C /udu)d/du(-C /udu)Cd/duCd/du（/udu)/udu)C /u C /u 2C /r2C /r r r 2C2C (r )(r )(2C )(2C )2C2CQ/2T CQ/2T CQ/4T Q/4T H HH0H0Q/4T /udu Q/4T /udu S SH0H0H HQ/4T /udu Q/4T /udu 令井函数令井函数W(u)W(u)/udu /udu ，S SQW(u)/QW(u)/（4T4T）(泰斯公式）泰斯公式）LOGO 第三第三节 单井出水量的非井出水量的非稳定流定流计算算 三三 泰斯公式泰斯公式讨论v 1 1、井函数、井函数W(U)W(U)可查表，可查表，W(U)W(U)随随U U值增大而减小；当承压井定流量抽水时，值增大而减小；当承压井定流量抽水时，降深降深S S随距离随距离r r的增大而减小，随抽水时间延长而增大。的增大而减小，随抽水时间延长而增大。v井函数为负指数积分函数，展开后为一级收敛级数，当井函数为负指数积分函数，展开后为一级收敛级数，当u u值很小时忽略展值很小时忽略展开式第二项以后部分，化简后得：开式第二项以后部分，化简后得：v W(u)W(u)2.3lg2.25at/r2.3lg2.25at/r vS S2.3Q/4Tlg2.25at/r2.3Q/4Tlg2.25at/r u=0.01 u=0.010.10.1v2 2、降深、降深S S下降速度下降速度 Q /4Tt Q /4Tt v v 在一定距离范围内，当抽水时间足够长时，在一定距离范围内，当抽水时间足够长时，Q/4T/tQ/4T/t。v 说明在一定距离范围内，当抽水时间足够长时，水头下降速度与距说明在一定距离范围内，当抽水时间足够长时，水头下降速度与距离无关，等幅下降。离无关，等幅下降。v v3 3、流量变化规律、流量变化规律 Q/2TrQ/2Trv 由达西定律由达西定律:r:r处断面流量处断面流量QrQr2rKM =Q Q2rKM =Q Qv 不同过水段面流量不等。不同过水段面流量不等。v 4 4、井径无限小、井径无限小 (r )(r )Q/2T Q/2T 的目的使的目的使 趋近于趋近于1 1 U=ru*/4Tt U=ru*/4Tt 0.990.991 1，当，当U=0.01U=0.01时，由于井径无限小引起时，由于井径无限小引起的相对误差不超过的相对误差不超过1 1 。四四 潜水完整井非稳定流计算仿泰斯公式潜水完整井非稳定流计算仿泰斯公式1 1、当、当S=0.1HS=0.1H0 0时，时，T TKMKM改为改为T TKHKH0 0 代入泰斯公式。代入泰斯公式。2 2、当、当0.1H0.1H0 0S0.3HS0.3HS0.3H0 0时不能用泰斯公式。时不能用泰斯公式。五五 群井干扰抽水非稳定流计算群井干扰抽水非稳定流计算1 1、承压完整井群干扰抽水时抽水降深计算、承压完整井群干扰抽水时抽水降深计算 S1 S1S01(t2S01(t2t3t3tn),tn),井灌区任一点井灌区任一点P P水位消减值：水位消减值：SP SPS1S1S2S2S3S3SNSN Q1W(u1)/4TQ1W(u1)/4TQ2W(u2)/4TQ2W(u2)/4TQNW(un)/4TQNW(un)/4T当当u=0.01u=0.01时时 可用简化公式计算可用简化公式计算2 2、潜水完整井群干扰抽水时抽水降深计算、潜水完整井群干扰抽水时抽水降深计算 当当S=0.1H0S=0.1H0时，可按上式计算时，可按上式计算,其中其中T TKH0KH0 当当0.1H0S0.3H00.1H0S0.3H0时，采用仿泰斯公式按上式计算时，采用仿泰斯公式按上式计算LOGO 第二章第二章 地下水地下水资源源计算与算与评价价v第一第一节 基本概念基本概念v1 1、地下水、地下水资源源评价：价：对一个地区地下水一个地区地下水资源的源的质量，数量，量，数量，时v 空分布特征和开空分布特征和开发利用技利用技术要求作出定量分析，并要求作出定量分析，并评价其价其 v 开采价开采价值。它是地下水。它是地下水资源合理开源合理开发和科学管理的依据。和科学管理的依据。v2 2、评价任价任务：v（1 1）水）水质评价：水价：水质分析分析评价其可用性，水价其可用性，水质监测与防与防护措施。措施。v（2 2）水量）水量评价：价：计算各种水算各种水资源量，确定允源量，确定允许开采量及用水开采量及用水v 保保证率。率。v（3 3）开采技）开采技术条件条件评价：分析价：分析论证在在长期开采条件下是否会期开采条件下是否会v 引起不良的工程地引起不良的工程地质问题，并提出相，并提出相应的的预防措施。防措施。LOGO 第二章第二章 地下水地下水资源源计算与算与评价价v3 3、以水均衡、以水均衡为基基础的分的分类法（三分法）法（三分法）v 某一地下水含水某一地下水含水层单元，在某均衡元，在某均衡时段内，地下水的循段内，地下水的循环总是表是表现为补给-排泄排泄-储存量存量变化三种形式，它化三种形式，它们之之间在数量在数量上的均衡关系上的均衡关系为V V补VV排排=VV，该方法方法较大大储量分量分类法合理。法合理。v（1 1）补给量：量：v 单位位时间内内进入某一入某一单元含水元含水层的重力水体的重力水体积，分，分为天然天然v 补给量，人工量，人工补给量，开采量，开采补给量。量。v 天然天然补给量：天然状量：天然状态下下进入含水入含水层的水量的水量v （降雨入渗（降雨入渗补给，地表水渗漏，地表水渗漏，邻区地下水区地下水测渗）。渗）。v 人工人工补给量：人工引水入渗量：人工引水入渗补给地下水的水量。地下水的水量。v 开采开采补给量：开采条件下，除天然量：开采条件下，除天然补给量之外，量之外，额外外获得得的的补给量。开采形成降落漏斗使地表水体向地下含水量。开采形成降落漏斗使地表水体向地下含水层补水，水，或定托渗流或定托渗流变为自由渗流加大自由渗流加大补给量，或使量，或使邻区地下水加大区地下水加大侧渗量。渗量。v需要需要说明：明：计算算补给量是地下水量是地下水资源源评价的核心内容；价的核心内容；一般先一般先计算算现状条件下的状条件下的补给量，再估算开采后可能量，再估算开采后可能获得的得的额外外补给量。量。LOGO 第二章第二章 地下水地下水资源源计算与算与评价价v（2 2）排泄量）排泄量v 单位位时间内从某一内从某一单元含水元含水层排泄出的重力水体体排泄出的重力水体体积，分，分为天然排泄量和人工开采量。天然排泄量和人工开采量。v 天然排泄量指潜水蒸天然排泄量指潜水蒸发，补给地表水体，地表水体，测向流入向流入临区。区。v 人工开采量是取水构筑物从含水人工开采量是取水构筑物从含水层取出来的地下水量。取出来的地下水量。v允允许开采量：开采量：指通指通过技技术经济合理的取水建筑物，在整个开采合理的取水建筑物，在整个开采期内水量和水位不超期内水量和水位不超过设计要求，水要求，水质水温水温变化在允化在允许范范围内，内，不影响已建水源地正常生不影响已建水源地正常生产，不，不发生危害性工程地生危害性工程地质现象的前象的前提下，提下，单位位时间内从水文地内从水文地质单元中能元中能够取得水量。取得水量。v允允许开采量的大小取决于开采量的大小取决于补给量，一般比量，一般比补给量小。量小。v（3 3）储存量存量v 储存在含水存在含水层的重力水体的重力水体积。容。容积储存量存量G GuHFuHF：潜水含水：潜水含水层中的重力水体中的重力水体积。弹性性储存量：承存量：承压含水含水层水水头下降下降顶板以板以上某一位置上某一位置时，由于含水，由于含水层弹性性压缩和水体和水体弹性膨性膨胀所所释放出放出的水量。的水量。G GuuSF SF S-S-承承压水位降低水位降低值m mLOGO 第二章第二章 地下水地下水资源源计算与算与评价价v4 4、评价区价区v 若某区域的水文气象条件、地若某区域的水文气象条件、地质构造条件、地貌条件、水构造条件、地貌条件、水文地文地质条件、岩性条件等比条件、岩性条件等比较接近，整个区域作接近，整个区域作为一个一个计算算评价区，否价区，否则划分划分计算算评价价单元。元。v水利部水文局水利部水文局编制的地下水制的地下水资源源调查评价工作价工作细则指出；指出；v 按地形地貌特征和地下水按地形地貌特征和地下水类型，把型，把评价区划分价区划分为平原区、平原区、山丘区、沙漠区和内山丘区、沙漠区和内陆闭合盆地平原区，称合盆地平原区，称为一一级区。区。v 一一级区按水文地区按水文地质条件又划分条件又划分为若干个水文地若干个水文地质区，区，为二二级区。区。v 二二级区按地下水埋深、包气区按地下水埋深、包气带岩性的不同，再分岩性的不同，再分为若干若干个个均衡均衡计算区，算区，为三三级区。均衡区。均衡计算区是算区是计算水算水资源的最小源的最小单元。元。LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价一一、水均衡法水均衡法（一）适用条件（一）适用条件 地下水埋藏地下水埋藏较浅，地下水的浅，地下水的补给和排泄条件容易和排泄条件容易查清的地区。清的地区。对于干旱或半干旱山前洪于干旱或半干旱山前洪积平原和喀斯特地区，某些河谷地平原和喀斯特地区，某些河谷地区以及封区以及封闭的自流盆地，使用效果的自流盆地，使用效果较好。好。对深深层承承压含水含水层或山或山区基岩裂隙含水区基岩裂隙含水层不宜使用。不宜使用。（二）基本原理（二）基本原理 水量平衡方程式：水量平衡方程式：Q Q补Q Q排排uFh/tuFh/t Q Q补地下水地下水总补给量（量（m);m);Q Q排排地下水地下水总排泄量（排泄量（m)m)；uFh/t uFh/t单位位时间内含水内含水层中中储存量的存量的变化量（化量（m/d)m/d)t t均衡期（年包括干旱年，平水年，丰水年）均衡期（年包括干旱年，平水年，丰水年）h h均衡期均衡期评价区地下水位平均价区地下水位平均变幅（幅（m m）F F均衡区面均衡区面积（m m)LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v（三）开采条件下的均衡方程（三）开采条件下的均衡方程v 天然状天然状态下，地表水，地下水，大气水依靠自然条件相互下，地表水，地下水，大气水依靠自然条件相互转化，在一定周期内，多年化，在一定周期内，多年调节，天然，天然补给量与天然消耗量近量与天然消耗量近似相等，似相等，维持持动态平衡。开采条件下，打破天然状平衡。开采条件下，打破天然状态下的下的动态平衡，在天然渗流平衡，在天然渗流场的基的基础上叠加了开采渗流上叠加了开采渗流场即降落漏斗。即降落漏斗。由于人工开采，使得天然由于人工开采，使得天然补给量增加，天然消耗量减少。量增加，天然消耗量减少。v在一个均衡期内，合理的开采会建立新的在一个均衡期内，合理的开采会建立新的动态平衡。否平衡。否则地下地下水不断消耗，水位持水不断消耗，水位持续下降，不存在平衡。下降，不存在平衡。v水量均衡方程式水量均衡方程式v （Q Q入入Q Q出出）（）（W WQ Q开开）uFh/tuFh/tv（Q Q入入Q Q出出）侧向向补给量与排泄量之差量与排泄量之差（m/d)m/d)（W WQ Q开开）垂直）垂直补给量与排泄量之差量与排泄量之差（m/d)m/d)v u u含水含水层平均平均给水度水度 ，其它参数同上式，其它参数同上式v W WQ Q雨雨Q Q河河Q Q越越.Q Q蒸蒸LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v（四）（四）计算步算步骤v1 1 划分均衡区，确定均衡期划分均衡区，确定均衡期v 均衡期一般均衡期一般为一个水文年，包括丰水年，平水年和枯水年一个水文年，包括丰水年，平水年和枯水年在内的一个多年均衡期。在内的一个多年均衡期。v2 2 确定均衡要素，建立均衡方程确定均衡要素，建立均衡方程v3 3 地下水地下水资源源评价价v Q Q开开（Q Q入入Q Q出出）W WuFh/tuFh/tv 给出均衡期地下水位允出均衡期地下水位允许变幅的条件下，幅的条件下，计算地下水开采量，算地下水开采量，分析分析评价地下水价地下水资源源对用水的保用水的保证程度。程度。LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v二、二、开采开采试验法法v 适用于完全没有水文地适用于完全没有水文地质资料，水文地料，水文地质资料料难以以查清的清的地区，当急需确定地下水允地区，当急需确定地下水允许开采量开采量时。可采用打井或利用。可采用打井或利用现有井按需要的开采流量有井按需要的开采流量进行抽水行抽水试验，依据，依据试验结果果计算地下算地下水允水允许开采量。开采量。对潜水和承潜水和承压水，新水源和旧水源地均适用，水，新水源和旧水源地均适用，由于群井抽水由于群井抽水费用很高，只适用于中小型水源地。一般在枯水用很高，只适用于中小型水源地。一般在枯水季季节按开采条件做抽水按开采条件做抽水试验，延，延续几个月几个月。v(一）抽水出一）抽水出现稳定状定状态时允允许开采量开采量Q QP P与影响半径与影响半径R R0 0确定确定v 按按设计开采量开采量长时间抽水，若水位达到抽水，若水位达到设计降深降深S S0 0后，一后，一直能保持直能保持稳定，且停抽后水位又能定，且停抽后水位又能较快恢复到原水位，快恢复到原水位，说明抽明抽水量小于开采条件下的水量小于开采条件下的补给量，按量，按这样的流量开采是有保的流量开采是有保证的。的。v 延延长三个三个观测孔降深曲孔降深曲线（稳定定时）与静水位交点，至抽）与静水位交点，至抽水井中心的距离水井中心的距离为影响半径影响半径R R0 0。LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v(二二)抽水出现非稳定状态时允许开采量抽水出现非稳定状态时允许开采量Q QP P与影响半径与影响半径R R0 0确定确定v 按按设计开采量开采量长时间抽水，若水位达到抽水，若水位达到设计降深降深S S0 0后，并后，并不不稳定，停抽后水位有所恢复，但不恢复到原水位，定，停抽后水位有所恢复，但不恢复到原水位，说明抽水明抽水量大于开采条件下的量大于开采条件下的补给量，按量，按这样的流量开采是没有保的流量开采是没有保证的。的。v在水位持在水位持续下降的下降的过程中，当大部分漏斗出程中，当大部分漏斗出现等幅下降后，任等幅下降后，任意意时段的水量平衡段的水量平衡为：（：（Q Q抽抽-Q-Q补）t=uFSt=uFSv uF-uF-水位下降水位下降1 1米米时储存量的减少量存量的减少量m m v S-t S-t时段的水位降低段的水位降低m mv Q Q抽抽-t-t时段的平均抽水流量段的平均抽水流量m m/d/dv Q Q补-开采条件下的开采条件下的补给流量流量m m/d/dv水量平衡方程水量平衡方程 Q Q抽抽Q Q补uFS/tuFS/tv认为Q Q补和和uFuF变化不大，当作常数，化不大，当作常数，则将抽水流量比将抽水流量比较稳定，水定，水位下降比位下降比较均匀的若干均匀的若干时段的段的资料代入上式解方程料代入上式解方程组求常数。求常数。LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v用水位下降用水位下降时的的资料求允料求允许开采量开采量QpQpv Q Q抽抽1 1Q Q补uFS1/t1uFS1/t1v Q Q抽抽2 2Q Q补uFS2/t2 -uFS2/t2 -v Q Q抽抽n nQ Q补uFSn/tnuFSn/tnv 求求Q Q补和和uFuF的平均的平均值。v再用水位恢复再用水位恢复时的的资料料检验校核校核Q Q补的可靠性。的可靠性。v若在抽水若在抽水过程中减少抽水量使程中减少抽水量使Q Q抽抽QQ补,则地下水位会等幅回升，地下水位会等幅回升，v Q Q补Q Q抽抽uFS/tuFS/tv uF uF前面求出的前面求出的单位位储量的平均量的平均值m m v S/t S/t水位等幅回升速度水位等幅回升速度m/dm/dv当停止抽水当停止抽水时，Q Q抽抽0 0，则Q Q补uFS/tuFS/t。v允允许开采量开采量QpQpQ Q补LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v例例题一一某水源地位于基岩裂隙水的富水段，在某水源地位于基岩裂隙水的富水段，在0.2km0.2km 面面积内打内打了了1212个个钻孔，最大孔距不超孔，最大孔距不超过300300米。在其中米。在其中3 3孔中孔中进行了行了4 4个多个多月的抽水月的抽水试验，观测数据数据见下表。下表。试验表明，在水位迅速下降表明，在水位迅速下降阶段段结束后，开始等幅持束后，开始等幅持续下降，停抽或抽水量减少下降，停抽或抽水量减少时，水位，水位都有等幅回升都有等幅回升现象，象，这表明正常抽水已表明正常抽水已经大于大于实际补给量。量。v求解：求解：选用用5 5月月1 1日至日至6 6月月3030日各日各时段的段的资料分料分别代入公式代入公式v 3169 3169Q Q补0.47uF 27730.47uF 2773Q Q补0.09uF0.09uFv 3262 3262Q Q补0.94uF 30710.94uF 3071Q Q补0.54uF0.54uFv 2804 2804Q Q补0.19uF0.19uF时间（月日）时间（月日）5.15.255.266.26.76.106.116.196.206.30平均抽水流量m/d31692773326230712804水位平均下降速度m/d0.470.090.940.540.19LOGO 第二节第二节 地下水资源数量计算与评价地下水资源数量计算与评价v Q Q补和和uFuF计算表算表v 水位恢复水位恢复计算表算表 v 联立方号联立方号（1 1）和）和（2 2）（3 3）和）和（