地下水资源评价

1、地下水水量评价： 是对地下水源地或某一地区、某个含水层的补给量、储存量，允许开采量进行计算（包括地表水、地下水的地的基础上，对所用计算方法的适宜性、水文地质参数的可靠性、资源计算结果精度、开采资源保证程度所做出的全面评价。水资源调查评价工作，就是要回答一个地区或流域有多少水量区分布、时间变化、质量标准、可靠程度） 。同时还要研究社会经济发展需要多少水量（各种用水的现状，近期和远景预测） ，以及供需平衡存在的问题。地下水资源评价方法 用于确定地下水资源数量的方法很多， 这里主要介绍一下 4 种评价方法： 开 采试验法、 补给疏干法、水文分析法、开采强度法。1、开采试验法在地下水的非补给期 （ 或

2、枯水期 ）按接近取水工程设计的开采条件进行较长时间的抽水试验，然后根据抽水量、水位降深动态或开采条件下的水量均衡方程求解出水源地枯季补给量，并以此量作为水源地的允许开采量。1、1适用条件在水文地质条件复杂地区，如果一时很难查清补给条件而又急需做出评价是，则可打勘探开采孔，并按开采条件（开采降深和开采量）进行抽水试验，根据试验结果可以直接评价开采量，这种评价方法，对潜水或承压水，对新水源地或旧水源地扩建都能适用。对于含水性不均匀的岩溶地区最为常用。主要适用于中小型水源地。该方法的缺点是不能做区域性的水资源评价。1 、2 计算方法完全按开采条件抽水，最好从旱季开始，延续一至数月，从抽水到恢复水位进

3、行全面贯彻，结果可能出现两种情形：1）稳定状态：在长期抽水过程中，如果水位达到设计降深并趋于稳定状态，抽水量大于或等于需水量；抽停后，水位又能较快恢复到原始水位。则说明抽水量小于开采条件下的补给量，按需水量开采是有补给保证的，这时，实际的抽水量就是要求的开采量。2）非稳定状态 : 如果水位达到设计降深并不稳定，继续下降；停抽后，虽然水位有所恢复，但始终达不到原始水位，测说明抽水量已经超过开采条件下的补给量，按需水量开采是没有保证的，这时，可按下列方法评价开采量：在水位持续下降过程中，只有大部分漏斗开始等幅下降，降速大小同抽水量成比例，则任意时段的水量均衡应满足下式：卩FJ= （Q抽-Q补）t卩

4、单位储存量，m3S t时段的水位降，m3Q抽一平均抽水量m ?d3Q补一开采条件下的补给量 m ?d 由此得岀：:应取负号，故,其中抽水量有两部分组成：一是开采条件下的补给量；二是含水层中消耗的储存量。在抽水过程中，如果抽水量小于补给量，则水位应发生等幅回升，这时SQ 补=Q 抽 + jiF其中卩F取已求的平均值;S为等幅回升速度。停抽时，Q抽=0,由此得根据以上所求的Q补，结合水文地质条件和需水量即可评价开采量，但由此求得的Q补评价是偏保守的，因为，旱季抽水只能确定一年中最小的补给量，所以Q补用年平均补给量或多年平均补给量进行评价。1、3实例某水源地位于基岩裂隙水的富水地段，在面积内打了12

5、个钻孔，最大孔距不超过300叶 在其中的三个孔中进行了四个多月的开采抽水试验，观测数据见表1 1。1中的资料可列岀五个方程式:表1 1时段(月日)3平均抽水(m ?d)31692773307130712804平均降速 (m ?d)这些数据表明，在水位急速下降阶段结束后，开始等幅持续下降，停抽或暂时中断抽水以及抽水量减少时，都发现水位有等幅回升现象。这说明抽水量大于补给量。利用表 3169=QyF 2773=Q+F 3262=Q+F 3071=Q+F 2804=Q补+卩F为了全面考虑，把五个方程搭配联解，求岀Q补和卩值，结果见表1 2。联立方程和和和和平均值Q补26792813268826592

6、7101042473611763723从计算结果看，由不同时段组合所求岀的补给量相差不大，但卩F直变化较大，可能是由于裂隙发育不均，降落漏斗扩展速度不匀所致。再利用水位恢复资料进行复核Q补，数据及计算结果见表1 3。表1 3时段（月 日）水位恢复值（m；(m/d)平均抽水量（ m/d）公式补给量（m/d）SQ# =卩匚SQ补=Q抽+平均值2657从以上计算结果看，该水源地旱季的补给量在m/d之间，以此作为开采量是完全有保证的。若不能满足需水量的要求，还可以利用年内暂时储存量，适当增大允许开采量。但还应考虑总的降深大小及评价开采后对环境的影响。2、补给疏干法根据水均衡的原理和以丰补欠的原则，把丰

7、水期多余的地下水补给量（即大于开采量的那一部分补给量）平均分配到枯水期进行开采的资源评价方法。2、1适用条件补偿疏干法适用于蓄水范围不大，仅有季节性补给，且有一定储存量，能够其调节作用的季节性的调节水源地。在半干旱地区，降雨季节性分布极不均匀，雨季时间短、降雨集中，地下水开采在旱季以来于消耗含水层的储存量而在雨季以回填被疏干的地下库容的形式进行补给。开采量多少取决于允许降深范围如何最大限度地利用储存量的调节库容。采用这种评价方法时， 它要求具备以下两个条件： 一是可借用的储存量必须满足旱季的连续稳定开采;二是雨季补给必须在平衡当时开采的同时，保证能全部补偿借用的储存量而非部分补偿。2、2计算方

8、法用补偿疏干法评价，要进行抽水试验，要求有两点：抽水量大小，必须造成动水位等幅下降，以便观测代表整个漏斗的下降值；抽水时间，应包括观测到整个漏斗的等幅上升值。在旱季漏斗斗幅下降过程中，任意时间段内储存量的变化值，应该等于该时段抽岀的水体积，即:卩FS= QtS时段t内漏斗的等幅下降值;Q为抽水量rm/d卩F单位储存量;卩一给水度;F漏斗面积;当漏斗扩展全区时，卩值接近常量，则:Qit Qi (ti - to) !1F=Ss1 -SoQ1 旱季的定量抽水量;So 水位急速下降结束时刻to的水位降;据求岀的卩值，分两步对开采量进行评价。计算开米量，旱季可能借用的储存量，必须保证整个旱季连续开采，所

9、以旱季末期形成的最大水位降深不得超过设计的允许降深。设允许降深为Smax，S = Smax - S);旱季开采时间设为t开，则t= t开-to。由此可以得岀开采量：Q开=”Smax -S 0Smax -S 0iF tH t开-t ot开因为t开to，略去t0更安全些。用上式求岀的 Q开，可保证旱季连续开采，不会中断，但不一定有补给保证。(2)计算补给量和评价，等幅回升时的单位补偿量和水位下降时的单位储存量相等。设雨季抽水过程中测得水位回升值为S，经过时间为t，则单位时间内补偿的水体积为卩fF。如用t补表示雨季的总补给时间，则雨季补给的水体积为(hF+ Q2 ) t补。把这个体积分配为全年开采时

10、：即得年平均补给量:s°补=365(Hf+Q2)q2 雨季开采量，为了供水安全，考虑到可能岀现旱年系列时，应从多年气象周期岀发，采用安全系 数r=。这时t补=rT补，T补为勘察年的时间补给时间。10km 2。灰岩分布区有间歇2、3实例某新建水源地，据勘探查明：含水层为厚层灰岩，呈条带状分布，面积约水试验，试验资料归纳如图性河，故岩溶水的补给来源主要是季节性河水渗漏和降水渗入。为了评价开采量，在整个旱季做了长期抽1所示，勘察年的旱季时间t开=253天，两季补给时间为 T补=112天，允许降深规定为Smax = 23m。解：按旱季抽水资料求岀(iF直,3?d14.53 - 5Q1(t1-

11、t0)= 1761.7 丄150 - 10) = 25880 m M 地下径流模数，m ?skm2? S1 - S0把允许降深作为旱季末期的最大降深，令t开=253天，则Smax - S0Q开=i=25880 23 - F含水层面积，km2 由此可知，地下径流模数是评价区域地下水资源的重要指标，它受区域地下水的补给、径流、排泄条 件所控制。因此结合不同的水文地质特征采用不同的方法进行评价: = 1841.2 m3?d t 开253d取安全系数r= t补=rT补=0.7 x 112 = 88.6天，得岀Q 补=邑(ifF+ Q2 ) = 886 (25880 x 11+ 1900 ) = 196

12、3.32 m3?乂补 365 t 236549d3由此可得，Q补 Q开，故Q开=1841 m ?d，是有补给保证又能取岀来的开采量。3、水文分析法在查明水文地质条件的基础上，充分利用水文测流资料和测流控制区的含水层面积，直接求岀地下径 流模数，即单位时间点位面积含水层的补给量或地下径流量。3、1适用条件在水文地质勘察的基础上，需查明地下水的天然补给量，作为有保证的区域地下水资源，评价区域地 下水资源的方法较多，但目前国内采用研究地表径流的水文分析发比较成功。尤其在水文地质条件复杂、 研究程度又相对较低的岩溶水或裂隙水分布区，用这种方法评价比较简单有效。3、2计算方法根据地下径流模数，可以间接推

13、算区域地下水的天然补给量或地下径流量:Q= MF3Q地下径流量，m ?s1、地下河系发育的岩溶区根据这种水文地质特征，可选择有控制性的暗河岀口或泉群，测定其枯水期流量，同时圈定对应的 地下流域面积，取流量和地下流域面积之比，就是要求的地下径流模数。2、地表河系发育的非岩溶区对于裂隙水或岩溶裂隙水和积极交替带的孔隙水,补给量形成地下径流后，直接排入河谷变成河水流量的组成部分，故可充分利用水文站现成的河流水文图来确定地下径流模数。河水通常是由大气降水和地下水补给，在枯水期，河水流量几乎全由地下水维持，而洪水期河水流量的 大部分为降水补给，地下水补给量相对减少，甚至河水倒流补给地下水，因此，利用河流

14、水文图时，必须 从实际水文地质条件岀发，将地下径流量分割岀来。目前，分割界限常由经验确定。对岩性单一，集水面积较小的水文站，在流量过程图上涨部分的起涨点至退水部分的退水转折点之间连线，把该线以下部分作为基流量; 对岩性非均一，集水面积大的水文站，以枯水期平均流量代表基流量; 在没有水文站时，也可沿河流上下游断面布置简易测流法，由上下游断面的流量差可求的控制区的地下径流量和相应的地下径流模数。 当一个含水层和另一个模数已知的含水层一起被河流排泄时，可按下式计算未知含水层的模数,Q- MiFiF2M 23M2 未知含水层的径流模数，m ?skm2 ；F2 对应M2的含水层面积，km 2;Q含水综合

15、体排泄地段上的基流量,M1和Fj 已知的含水层面积和径流模数;3、3实例我国广西水文地质队，在地苏、大化等岩溶地区采用水文分析评价地下水资源，同时用实测流量进行了检验。结果，平均准确度达86%。具体见表3 1。编号露头类型地下径流模数m/s km补给面积km计算流量m7s实测流量m/s准确度%43地下河天窗1559163608346577101140平均地下河出开采强度法:在大范围的平原开采区，可将井位分布较均匀、水井流量相差不大的区域概化成一个或几个规则形状 的开采区，将分散井群的总流量概化为开采强度。然后按非稳定流的面积井公式去推算设计水位降深条件 下的开采量或给

16、定开采量条件下某一时刻开采区中心的水位降深。这种方法即为开采强度法。4、1适用条件，采用开采强度法计算开采量比较方在井数很多，井位分散、开采面积很大的地区（这是农业供水的特点）便。4、2计算方法dF=dS (x,y,t)£4 p? at l -(xt lx e 4a T J ( J0 lxVnT-l y-(y- n) lx e 4a Td nd T v"nT以无界承压含水层中的矩形开采区为例，在矩形开采区内，以（En）点为中心，取一微面积，并把它看成开采量为 dQ的一个井点，在此点井作用下，开采区内外将形成水位降深的非稳定场,对任一点A （x、y）引起的水位降 ds，用点函

17、数表示:t -r 2dQ “ e斗 ds = J - d T4 nT 0 TT导水系数;A导压系数;t 时间;r点井到A点的距离;A点的总水位降:开采强度公式:(X，y，t)=S?( a,闭( a,伍)+ sx( a 302( a 切lx-Xa= 21?Ix+Xly-yly+ya =,伦=於（a，3）= Jq 0（胡 dT2 z 20 (z)=聂Jez dz 几分概率S?（ a ,3的数值查表。在资源评价中，人们最关心的地方时开采区的中心降深最大的部位，这里最易超过允许降深引起掉泵 停产，故令x=y=0,= Sx（ a,）则其中(X =lx2var 'S（t）= 4Sx（a ,3)如果

18、浅水层厚度SH过大，而水位将 S相对较小，即-V时，则可以直接近似用于无界含水层，计算结果不会HS如果 匚V时，要用H过分歪曲实际。丄（H 2（h - h），表示开采漏斗内浅水层的平均厚度;- h2）代替S,用给水度M弋替，结果得:2h cH2- h2 = ：thcSx( a, 3)+ Sx( a, 33) + Sx( a, 3)+ Sx( a, 32)2 MH2-h0=Sx(山其中的hh表示任一点的动水位；ho代表开采区中心的动水位。4、3实例河北省冀县、150 250m，下部含水枣强、衡水地区，位于河北平原中部，有巨厚的第四纪沉积层，形成良好的储水条件。其中有两个承压水含水组，是目前工农业

19、供水的主要开采层。上部含水组在地表下组在250 350m之间。二含水组均为中细砂组成。随着工农业的发展，开采量逐年扩大，已经形成以衡水为中心的巨大开采漏斗。实践证明，由于距补给区很远，主要消耗弹性储存量， 所以形成非稳定开采动态：历年水位下降大于水位回升，每年平均下降,开采量已经失去补给保证。同时，下部含水组的水位下降快而回升慢，水位高于上部含水组，两组的开采漏斗也不重合。所以，两个含水组之间的水力联系并不明显，而有一定的独立性。为了满足农田水利化 2030%的规划要求，应对两个含水组中的地下水资源作岀评价。为了简单起见,本例仅摘录上部含水组的计算结果，说明计算和评价方法。上部含水组的历年开采

20、资料统计在表4 1 中。年 月曰天数统计开采量开采面积开采强度漏斗中心水位降深一222207. 644一283165. 0644一385296. 364430413165081316222/1316矩形开采区概化图在边界条件没有完全查清以前，现有开采面积虽已超过1000km2以上，但同河北平原面积相比还是很小的一部分，而且离补给区很远，含水层可视为无限大。所以，仍属局部开采区，采用开采强度法计算比较合适。(1)确定水文地质参数。把表4 1中第一和第二两行中的数据带入公式得：6.7 =0.000212 X 222 x , ?Sx (550020002Xax 222,2XV222)7.05=0.0

21、00212X 505 sx (二02 X va X 50520002 X Vx 505)可得上部含水组的参数:(0.000132 - 0.000212) +?a2a= 7.5 X 104m ?d，a?=55002 X Va X 28320002 XVX 283(2)，计算1986 1973年的开采量，验证所求参数的可靠性。分两种情况计算:开采区有同一开采强度的1986 1970 年。开采面积为44km2。见表4 1中的图示(1)。1968 年，t=222 日，a= X104m2d, lx=5500m,ly=2000m。求得公式中的lx5500a = = /2 Vat 2 X V 7.5X 22

22、2 X1040.372B= -l= 0.23622 vat查的Sx（ a 3 ）= 0.272由此可得6.7 =:蔦；2 X 0.2720.00258規=0.0002该年总开采量为£68 X t XF = 0.0002 X 222 X 44 = 196.33m2?y?y o当年的统计开采量为xa。二者相比，计算的比实际的偏小%。同理，可求得 1969年和1970年的开采量，列入表42中。开采区有不同开采强度的1971 1973 年。开米面积为1316km2,见表4 1中（2）所示。这时开采强度不但历年不同,不同地段也不一样。所以对1971年来说,虚线地段的中心水位降,按迭加原理为:S

23、71 =(68 t1194?TlxSx( ''2xvat7194'2xvat1194、(69- £68)(t1194-t 222)心 ( lx)+TS (ly2XVat72'2Xvar972) +(60- £69)( t1194 -t 505)E lxly( 61 - £70)( 11194 -t 890 ) xT?( 2XVat689 , 2 X vat 689 )2 t?( 2XVa(t 1194-t 890 )lyly -y2xva(t 1194 -t 890)xlx ly +y4 2 中。S ( 2XVa(t 1194 -t 890),2X Va(t 1194 -t 890)"由此式可求岀£n1。同理也可求岀 勺2和6 3。换成年总开采量后，结果也列入表表项目196819691970197119721973统计 资 料漏斗中心水位降深（m）3统计开采量（万m ?d）/开采天数（d）222283385304508222计算 结 果44km开采强度开采量3（万 m ?d）1316km开采强度开采量3（万 m ?d）误差%+/表中数字比较证明，计算结果和统计结果资料很相近，最大误差均在用的公式基本上符合本区的实际情况。10%以内。可见，所求参数和采20%,灌溉标准为(3)，按规划的需水