资源描述
地下水资源评价 Document number:WTWYT-WYWY-BTGTT-YTTYU-2018GT 地下水资源评价 地下水水量评价:是对地下水源地或某一地区、某个含水层的补给量、储存量,允许开采量进行计算的基础上,对所用计算方法的适宜性、水文地质参数的可靠性、资源计算结果精度、开采资源保证程度所做出的全面评价。水资源调查评价工作,就是要回答一个地区或流域有多少水量(包括地表水、地下水的地区分布、时间变化、质量标准、可靠程度)。同时还要研究社会经济发展需要多少水量(各种用水的现状,近期和远景预测),以及供需平衡存在的问题。地下水资源评价方法:用于确定地下水资源数量的方法很多,这里主要介绍一下 4种评价方法:开采试验法、补给疏干法、水文分析法、开采强度法。1、开采试验法 在地下水的非补给期(或枯水期)按接近取水工程设计的开采条件进行较长时间的抽水试验,然后根据抽水量、水位降深动态或开采条件下的水量均衡方程求解出水源地枯季补给量,并以此量作为水源地的允许开采量。1、1 适用条件 在水文地质条件复杂地区,如果一时很难查清补给条件而又急需做出评价是,则可打勘探开采孔,并按开采条件(开采降深和开采量)进行抽水试验,根据试验结果可以直接评价开采量,这种评价方法,对潜水或承压水,对新水源地或旧水源地扩建都能适用。对于含水性不均匀的岩溶地区最为常用。主要适用于中小型水源地。该方法的缺点是不能做区域性的水资源评价。1、2 计算方法 完全按开采条件抽水,最好从旱季开始,延续一至数月,从抽水到恢复水位进行全面贯彻,结果可能出现两种情形:(1)稳定状态:在长期抽水过程中,如果水位达到设计降深并趋于稳定状态,抽水量大于或等于需水量;抽停后,水位又能较快恢复到原始水位。则说明抽水量小于开采条件下的补给量,按需水量开采是有补给保证的,这时,实际的抽水量就是要求的开采量。(2)非稳定状态:如果水位达到设计降深并不稳定,继续下降;停抽后,虽然水位有所恢复,但始终达不到原始水位,测说明抽水量已经超过开采条件下的补给量,按需水量开采是没有保证的,这时,可按下列方法评价开采量:在水位持续下降过程中,只有大部分漏斗开始等幅下降,降速大小同抽水量成比例,则任意时段的水量均衡应满足下式:F?S=(Q抽 Q补)t F单位储存量,m3 St时段的水位降,m Q抽平均抽水量m3d Q补开采条件下的补给量m3d 由此得出:Q抽=Q补+FSt 其中抽水量有两部分组成:一是开采条件下的补给量;二是含水层中消耗的储存量。在抽水过程中,如果抽水量小于补给量,则水位应发生等幅回升,这时St应取负号,故,Q补=Q抽+FSt 其中F取已求的平均值;St为等幅回升速度。停抽时,Q抽=0,由此得 Q补=FSt 根据以上所求的Q补,结合水文地质条件和需水量即可评价开采量,但由此求得的Q补评价是偏保守的,因为,旱季抽水只能确定一年中最小的补给量,所以Q补用年平均补给量或多年平均补给量进行评价。1、3 实例 某水源地位于基岩裂隙水的富水地段,在面积内打了12 个钻孔,最大孔距不超过 300m。在其中的三个孔中进行了四个多月的开采抽水试验,观测数据见表11。表 11 时段(月日)平均抽水(m3d)3169 2773 3071 3071 2804 平均降速 (md)这些数据表明,在水位急速下降阶段结束后,开始等幅持续下降,停抽或暂时中断抽水以及抽水量减少时,都发现水位有等幅回升现象。这说明抽水量大于补给量。利用表 1 中的资料可列出五个方程式:3169=Q补+F 2773=Q补+F 3262=Q补+F 3071=Q补+F 2804=Q补+F 为了全面考虑,把五个方程搭配联解,求出Q补和F值,结果见表 12。表 12 联立方程 和 和 和 和 平均值 Q补 2679 2813 2688 2659 2710 F 1042 473 611 763 723 从计算结果看,由不同时段组合所求出的补给量相差不大,但F值变化较大,可能是由于裂隙发育不均,降落漏斗扩展速度不匀所致。再利用水位恢复资料进行复核Q补,数据及计算结果见表 13。表 13 时段(月日)水位恢复值(m)st(m/d)平均抽水量(m3/d)公式 补给量(m3/d)Q补=FSt Q补=Q抽+FSt 平均值 2657 从以上计算结果看,该水源地旱季的补给量在m3/d 之间,以此作为开采量是完全有保证的。若不能满足需水量的要求,还可以利用年内暂时储存量,适当增大允许开采量。但还应考虑总的降深大小及评价开采后对环境的影响。2、补给疏干法 根据水均衡的原理和以丰补欠的原则,把丰水期多余的地下水补给量(即大于开采量的那一部分补给量)平均分配到枯水期进行开采的资源评价方法。2、1 适用条件 补偿疏干法适用于蓄水范围不大,仅有季节性补给,且有一定储存量,能够其调节作用的季节性的调节水源地。在半干旱地区,降雨季节性分布极不均匀,雨季时间短、降雨集中,地下水开采在旱季以来于消耗含水层的储存量而在雨季以回填被疏干的地下库容的形式进行补给。开采量多少取决于允许降深范围如何最大限度地利用储存量的调节库容。采用这种评价方法时,它要求具备以下两个条件:一是可借用的储存量必须满足旱季的连续稳定开采;二是雨季补给必须在平衡当时开采的同时,保证能全部补偿借用的储存量而非部分补偿。2、2 计算方法 用补偿疏干法评价,要进行抽水试验,要求有两点:抽水量大小,必须造成动水位等幅下降,以便观测代表整个漏斗的下降值;抽水时间,应包括观测到整个漏斗的等幅上升值。在旱季漏斗斗幅下降过程中,任意时间段内储存量的变化值,应该等于该时段抽出的水体积,即:FS=Qt S时段t内漏斗的等幅下降值;Q为抽水量 m3/d F 单位储存量;给水度;F漏斗面积;当漏斗扩展全区时,F 值接近常量,则:F=Q1ts=Q1(t1 t0)s1 s0 Q1旱季的定量抽水量;s0水位急速下降结束时刻t0的水位降;s1旱季末时刻的水位降;(见下图抽水试验过程图)根据求出的F值,分两步对开采量进行评价。(1)计算开采量,旱季可能借用的储存量,必须保证整个旱季连续开采,所以旱季末期形成的最大水位降深不得超过设计的允许降深。设允许降深为Smax,s=Smax S0;旱季开采时间设为t开,则t=t开 t0。由此可以得出开采量:Q开=FSmaxS0t开t0 FSmaxS0t开 因为t开t0,略去t0更安全些。用上式求出的Q开,可保证旱季连续开采,不会中断,但不一定有补给保证。(2)计算补给量和评价,等幅回升时的单位补偿量和水位下降时的单位储存量相等。设雨季抽水过程中测得水位回升值为S,经过时间为t,则单位时间内补偿的水体积为Fst。如用t补表示雨季的总补给时间,则雨季补给的水体积为(Fst+Q2)t补。把这个体积分配为全年开采时:即得年平均补给量:Q补=t补365(Fst+Q2)Q2雨季开采量,为了供水安全,考虑到可能出现旱年系列时,应从多年气象周期出发,采用安全系数 r=。这时t补=rT补,T补为勘察年的时间补给时间。2、3 实例 某新建水源地,据勘探查明:含水层为厚层灰岩,呈条带状分布,面积约10km2。灰岩分布区有间歇性河,故岩溶水的补给来源主要是季节性河水渗漏和降水渗入。为了评价开采量,在整个旱季做了长期抽水试验,试验资料归纳如图 1 所示,勘察年的旱季时间t开=253 天,两季补给时间为T补=112 天,允许降深规定为Smax=23m。解:按旱季抽水资料求出F 值,F=Q1(t1 t0)s1 s0=1761.7(150 10)14.53 5=25880m3d 把允许降深作为旱季末期的最大降深,令t开=253 天,则 Q开=FSmax S0t开=2588023 5253=1841.2m3d 取安全系数 r=,t补=rT补=0.7 112=88.6 天,得出 Q补=t补365(Fst+Q2)=88.6365(2588011.549+1900)=1963.32m3d 由此可得,Q补 Q开,故Q开=1841m3d,是有补给保证又能取出来的开采量。3、水文分析法 在查明水文地质条件的基础上,充分利用水文测流资料和测流控制区的含水层面积,直接求出地下径流模数,即单位时间点位面积含水层的补给量或地下径流量。3、1 适用条件 在水文地质勘察的基础上,需查明地下水的天然补给量,作为有保证的区域地下水资源,评价区域地下水资源的方法较多,但目前国内采用研究地表径流的水文分析发比较成功。尤其在水文地质条件复杂、研究程度又相对较低的岩溶水或裂隙水分布区,用这种方法评价比较简单有效。3、2 计算方法 根据地下径流模数,可以间接推算区域地下水的天然补给量或地下径流量:Q=M?F Q地下径流量,m3s M地下径流模数,m3s?km2 F含水层面积,km2 由此可知,地下径流模数是评价区域地下水资源的重要指标,它受区域地下水的补给、径流、排泄条件所控制。因此结合不同的水文地质特征采用不同的方法进行评价:1、地下河系发育的岩溶区 根据这种水文地质特征,可选择有控制性的暗河出口或泉群,测定其枯水期流量,同时圈定对应的地下流域面积,取流量和地下流域面积之比,就是要求的地下径流模数。2、地表河系发育的非岩溶区 对于裂隙水或岩溶裂隙水和积极交替带的孔隙水,补给量形成地下径流后,直接排入河谷变成河水流量的组成部分,故可充分利用水文站现成的河流水文图来确定地下径流模数。河水通常是由大气降水和地下水补给,在枯水期,河水流量几乎全由地下水维持,而洪水期河水流量的大部分为降水补给,地下水补给量相对减少,甚至河水倒流补给地下水,因此,利用河流水文图时,必须从实际水文地质条件出发,将地下径流量分割出来。目前,分割界限常由经验确定。对岩性单一,集水面积较小的水文站,在流量过程图上涨部分的起涨点至退水部分的退水转折点之间连线,把该线以下部分作为基流量;对岩性非均一,集水面积大的水文站,以枯水期平均流量代表基流量;在没有水文站时,也可沿河流上下游断面布置简易测流法,由上下游断面的流量差可求的控制区的地下径流量和相应的地下径流模数。当一个含水层和另一个模数已知的含水层一起被河流排泄时,可按下式计算未知含水层的模数,M2=Q M1F1F2 M2未知含水层的径流模数,m3s?km2;F2对应M2的含水层面积,km2;Q含水综合体排泄地段上的基流量,m3s;M1和F1已知的含水层面积和径流模数;3、3 实例 我国广西水文地质队,在地苏、大化等岩溶地区采用水文分析评价地下水资源,同时用实测流量进行了检验。结果,平均准确度达 86%。具体见表 31。表 31 编号 露头类型 地下径流模数 m3/skm2 补给面积 km2 计算流量 m3/s 实测流量 m3/s 准确度%43 63 4 10 11 40 平均 地下河天窗 155 60 65 18 14 44 91 83 77 90 92 88 地下河出口 4、开采强度法:在大范围的平原开采区,可将井位分布较均匀、水井流量相差不大的区域概化成一个或几个规则形状的开采区,将分散井群的总流量概化为开采强度。然后按非稳定流的面积井公式去推算设计水位降深条件下的开采量或给定开采量条件下某一时刻开采区中心的水位降深。这种方法即为开采强度法。4、1 适用条件 在井数很多,井位分散、开采面积很大的地区(这是农业供水的特点),采用开采强度法计算开采量比较方便。4、2 计算方法 以无界承压含水层中的矩形开采区为例,在矩形开采区内,以()点为中心,取一微面积 dF=dd,并把它看成开采量为 dQ 的一个井点,在此点井作用下,开采区内外将形成水位降深的非稳定场,对任一点 A(x、y)引起的水位降 ds,用点函数表示:ds=dQ4Ter24at0d T导水系数;A导压系数;t时间;r点井到 A点的距离;A点的总水位降:S(x,y,t)=4a(e(x)24alxlxde(y)24alxlyd)t0d 开采强度公式:S(x,y,t)=t4S(1,1)+Sx(1,2)+Sx(2,1)+Sx(2,2)1=lxx2at,2=lx+x2at,1=lyy2at,2=ly+y2at Sx(,)=()()10d,(z)=2x ez2dzz0几分概率 S(,)的数值查表。在资源评价中,人们最关心的地方时开采区的中心降深最大的部位,这里最易超过允许降深引起掉泵停产,故令 x=y=0,=Sx(,),则 S(t)=tSx(,)其中=lx2at ,=ly2at,如果浅水层厚度 H过大,而水位将 S 相对较小,即SH时,则可以直接近似用于无界含水层,计算结果不会过分歪曲实际。如果SH时,要用12hc(H2 h2)代替 S,用给水度代替,结果得:H2 h2=t2hcSx(1,1)+Sx(1,3)+Sx(2,1)+Sx(2,2)H2 h02=t2hcSx(,)其中的hc=12(H h),表示开采漏斗内浅水层的平均厚度;h表示任一点的动水位;h0代表开采区中心的动水位。4、3 实例 河北省冀县、枣强、衡水地区,位于河北平原中部,有巨厚的第四纪沉积层,形成良好的储水条件。其中有两个承压水含水组,是目前工农业供水的主要开采层。上部含水组在地表下 150250m,下部含水组在 250350m 之间。二含水组均为中细砂组成。随着工农业的发展,开采量逐年扩大,已经形成以衡水为中心的巨大开采漏斗。实践证明,由于距补给区很远,主要消耗弹性储存量,所以形成非稳定开采动态:历年水位下降大于水位回升,每年平均下降,开采量已经失去补给保证。同时,下部含水组的水位下降快而回升慢,水位高于上部含水组,两组的开采漏斗也不重合。所以,两个含水组之间的水力联系并不明显,而有一定的独立性。为了满足农田水利化 2030%的规划要求,应对两个含水组中的地下水资源作出评价。为了简单起见,本例仅摘录上部含水组的计算结果,说明计算和评价方法。上部含水组的历年开采资料统计在表 41 中。表 41 年 月 日 天数 统计开采量 开采面积 开采强度 漏斗中心水位降深 矩形开采区概化图 222 2076 44 283 16506 44 385 29636 44 304 1316 508 1316 222/1316 在边界条件没有完全查清以前,现有开采面积虽已超过 1000km2以上,但同河北平原面积相比还是很小的一部分,而且离补给区很远,含水层可视为无限大。所以,仍属局部开采区,采用开采强度法计算比较合适。(1)确定水文地质参数。把表 41 中第一和第二两行中的数据带入公式 得:6.7=0.000212222Sx(55002a222,20002222)7.05=0.000212 505Sx(55002 a 505,20002 505)+(0.000132 0.000212)283Sx55002 a 283,20002 283 可得上部含水组的参数:=7.5 104m2d,=(2),计算 19861973年的开采量,验证所求参数的可靠性。分两种情况计算:开采区有同一开采强度的 19861970年。开采面积为 44km2。见表 41中的图示(1)。1968 年,t=222 日,a=104m2d,lx=5500m,ly=2000m。求得公式中的 =lx2t=55002 7.5 222 104=0.372 =ly2at=0.2362 查的Sx(,)=0.272由此可得 6.7=68 2220.00258 0.272 68=0.0002 该年总开采量为68 t F=0.0002 222 44=196.33 万m2y。当年的统计开采量为a。二者相比,计算的比实际的偏小%。同理,可求得 1969年和 1970 年的开采量,列入表 42中。开采区有不同开采强度的 19711973 年。开采面积为 1316km2,见表 41中(2)所示。这时开采强度不但历年不同,不同地段也不一样。所以对 1971年来说,虚线地段的中心水位降,按迭加原理为:S71=68t1194Sx(lx2at1194,ly2at1194)+(6968)(t1194t222)Sx(lx2at972,ly2at972)+(7069)(t1194t505)Sx(lx2at689,ly2at689)+(7170)(t1194t890)2Sx(lx2a(t1194t890),lyy2a(t1194t890)+Sx(lx2a(t1194t890),ly+y2a(t1194t890)由此式可求出1。同理也可求出2和3。换成年总开采量后,结果也列入表 42 中。表 42 年份 项目 1968 1969 1970 1971 1972 1973 统计资料 漏斗中心水位降深(m)统 计 开 采 量(万m3d)/开采天数(d)222 283 385 304 508 222 计算结果 44 km2 开采强度 开采量(万m3d)1316 km2 开采强度 开采量(万m3d)误差%+/表中数字比较证明,计算结果和统计结果资料很相近,最大误差均在 10%以内。可见,所求参数和采用的公式基本上符合本区的实际情况。(3),按规划的需水量预测漏斗中心水位降深,根据规划要求,水利化程度为20%,灌溉标准为 200-300m3/y亩时,需水量和预测的水位降深,列入表 43 中。表 43 开采规划 预计降深 水 利 化程度%灌溉标准m3/y亩 需水量m3/ykm2 开 采 强度 m/d 75年 76年 77年 78年 79年 80年 20 200 40000 41348 49432 58768 63193 68980 73500 300 60000 48000 60642 72838 81474 90880 98010 目前采用的取水工具,主要是吸程60m的深井泵,去掉平均埋深后,允许降深以50m左右为宜。从表中数字可见,从1977年开始以后的水位降深均以超过允许降深,部分井将发生抽空吊泵,除非更换设备,否则不可能保持正常开采。(4),按控制降深 50m计算开采量和回灌量。为了保持正常开采,要用人工补给法控制水位降深。按规划要求,在1980年前,漏斗中心水位应当控制在50m以内。因此,每年平均允许下降约3m。计算结果列入表 44 中。表 44 年份 1974年 1975年 1976年 1977年 1978年 1980年 控制水位降深(m)计 算 开 采 量(m2y?km2)27900 25600 22600 22600 22400 21600 从历年规划的需水量中减去上表中的开采量,即得历年缺少的水量,这就是应当进行的回灌量。计算结果列入表 45 中。表 45 规划需水量 必需回灌量 水利化程度%灌溉标准(m3/y?亩)需水量(m2y?km2)1974年 1975年 1976年 1977年 1978年 1980年 20 200 40000 12100 14400 17400 17400 17600 18400 300 60000 32100 34400 37400 37400 37600 38400 结果表明,1980 年前必需的回灌量占需水量的 3164%。按这个比例进行回灌,才能保证规划的需水量。否则,就要中断开采。但是,按这个比例回灌能否成功,还有待实践研究。
展开阅读全文