1、井渠结合灌区地下水动态预报及适宜渠井用水比分析2006 年 2 月灌溉排水JournalofInigati0nand第 25 卷第 1 期文章编号:1000646X(2006)O1 一 OOO604井渠结合灌区地下水动态预报及适宜渠井用水比分析周维博,曾发琛(长安大学环境科学与工程学院,陕西西安 710054)摘要:以陕西注惠渠井渠结合灌区为例,根据灌区多年降雨量 ,渠灌用水量,井灌用水量资料厦渠井灌溉用水量的比值.利用多元非线性相关分析法建立了灌区地下水动态预报的数学模型.对减少渠灌用水量增加井灌用水量,坩加渠灌用水量减少井灌用水量,渠灌与井灌用水量合理比值 3 种情况下地下水动态进行了预测
2、,提出了适宜的渠井用水比例.为灌区地下水合理开采和灌溉水资源优化配置提供依据.关键词:井渠结合灌区;地下水动态预测;非线性模型; 渠井用水比中圈分类号:$273.4 文献标识码:A1 地下水动态预测模型陕西泾惠渠灌区为我国北方典型井渠结合灌区.灌区地下水主要靠灌溉和降雨人渗补给,灌溉和降雨入渗补给分别占总补给量的 52.6%和 32.4%.灌区地下水位的变化主要受大气降水,渠灌用水,井灌用水以及渠井灌溉用水比的影响.考虑到地下水位动态变化与影响因素之间为非线性关系,为此采用非线性回归分析法建立灌区地下水动态预测的模型.将降雨量,渠灌用水量,井灌用水量和渠井灌溉用水比作为预测模型的回归因子,即,
3、0(0)=boll(1)+6022(2)+bo33(3)+bo4u4(4)(1)式中:6.,bb6.4 为回归系数;()为预测的地下水埋深对应的标准正态变量值;,()为降雨量对应的标准正态变量值;() 为渠灌用水量对应的标准正态变量值;,(,) 为井灌用水量对应的标准正态变量值;()为渠井灌溉用水量比对应的标准正态变量值.将各个观测系列分别以相应的正态变量值-j(j=0,I,2,3,4)由公式 (2)转化可求得它们之间的线性相关系数 y*:_Il一式中:,为_,|j系列中序号为 i 的观测值的标准正态变量值.根据灌区 l6 年(19811996 年)实测数据资料进行四元非线性回归分析计算,得到
4、相关系数矩阵:=610.42480.83370.11110.74910.4248l0.13910.42270.36880.83370.1391l0.36230.49670.11110.42270.362310.54770.74910.36880.49670.54771收稿日期:2005113O作者筒介:周维博(1956 一),男,教授,博士.博士生导师.主要从事水资源与水环境及节水灌溉方面的研究.(3)6ol+“ybo3+“yj4bo4=7(.=1,2,3,4)(4)由相关系数矩阵和公式(4)可以列出线性方程组,并最终求解得回归系数:b 们=一0,2295,b:-0.5898,b=-0.015
5、4,bo4=-0.4169.将回归系数代入(1)式,得回归方程:zo(o)=一 0.2295ul(1)0.5898u2(2)一 0.0154u3(3)一 0.4169u4(4)(5)根据各个变量因子 xj(j=0,1,2,3,4)与标准正态变量(. =0,I,2,3,4) 之间的关系,可分别绘制出与的关系图.利用灌区的实测数据和非线性回归分析计算值,分别绘制出图 1图5,其中图 5 为地下水埋深与标准正态变量的关系图.可根据不同年或时段的降雨量渠灌用水量:,井灌用水量和渠井灌溉用水比的大小,则可以由图 1 一图 4 分别查出对应的标准正态变量的值H.,H:,的值,然后代入方程(5), 即可计算
6、出地下水埋深.对应的标准正态变量 H.值,再由图 5 可查出对应的预测地下水埋深.由复相关系数 R 的计算公式(6)和预报的标准误 S 公式(7):R/7ol6oI+7026o2+3o3b+7o4bo4(6)s,=(7)计算得出 R=0.900,S=0.436.根据分析计算中的年数=16 和变量因子个数_=4,当显着性在 0.05水平时,相关系数 R 的检验值 R=0.714,s=0.7,而前者大大优于后者,所以方程预测的结果是可信的.根据灌区 19811996 年地下水多年平均变量值和 1997 年的实际变量因子值,按照图 1 一图 5 中找出对应的H 值带入方程(5), 预测得到灌区地下水
7、多年平均水位埋深为=6.63m,1997 年的地下水位埋深为 Xo=11.60rtl,与实际灌区多年平均地下水埋深(6.24rt1)和 1997 年实际地下水埋深(11.92m)分别相差+0.39m 和一 0.32m,在误差允许范围之内,方程预测的结果与实测值基本一致./,-lf)/】 Jlllll-/,JI:年雨馈 I/mm一/.r+.,/渠滞 JIJ 水蕈/1i113图 1 降雨量.与标准正态变量值关系图图 2 渠灌用水量与 I 标准正态变量值“I 关系图_/一一 llHf,./f水最,/III3焉-,是 n差“_.“厂一./.I 浆井涪川水避比,图 3 井灌用水量 3 与标准正态变量值关
8、系图图 4.渠井灌用水量比与标准正态变量值关系图2 渠井用水量适宜比讨论井渠结合在地下水调控技术措施上可以实现灌区不产生渍涝和不会形成采补失调,并满足耕地的水资源供需平衡,其根本原因在于确定一个适宜的渠井用水比例.在井渠结合灌区内渠灌用水量和井灌用水量及渠井用水比例大小对灌区地下水动态的变化影响较大,渠井用水的适宜比例对调控灌区地下水位和满足作物用水需求以及灌区的可持续发展是很重要的.7i0LI_JIJl_lljlr:三:1-图 5 地下水埋深与标准正态变量值 Io 关系图渠灌用水量和井灌用水量受渠首引水量,地下水开采量,降雨和蒸发量,地下径流排泄等因素影响.近年来由于灌区内井灌水价低于渠灌水
9、价,渠灌不及井灌及时方便,群众则多用井灌致使地下水位大幅下降,一些地区还因井灌的地下水矿化度较高而使土壤积盐,农作物减产,破坏了地下水的采补平衡.自 1983 年以来由于灌区水源地泾河逐年减少引水量,因而使灌区地下水开采量逐年增加而使地下水位持续下降.根据灌区多年渠灌和井灌用水量和地下水位实际资料分析,灌区多年平均降雨量为 538mm,渠灌引水量为4.0 亿 m.,有一定的地下水蒸发和地下径流排泄的情况下 ,地下水埋深为 34m,渠灌与井灌用水量比作为 1.252 左右比较适当的.在地下水埋深大于 67m 以上,地下水蒸发和地下径流排泄较少,地表水灌溉引水薰在 2 亿 m.左右,渠灌与井灌用水
10、比为 1.5 的情况下,才能保持地下水补排平衡.在降雨量小于多年平均值的年份,应通过增加地表水的引水量解决灌溉用水问题.在灌溉季节地表水源不足的情况下超采地下水时,应尽量在当年非灌溉季节或在地表水充足年份及时引用地表水回补地下水,使地下水达到采补平衡引.3 不同灌溉水平下的预测由上述讨论可以看出灌区内渠灌与井灌用水量的比在 1.252 之间比较适宜,最佳比值为 1.5 左右.根据灌区 19811997 年多年灌溉用水量资料分析,多年平均灌溉用水量为 2 亿m 左右.假设灌区在降雨平水年(P=50%)和降雨偏枯水年(P=75%), 年灌溉用水总量在 2 亿 m 的水平下,对灌区内适宜减少渠灌用水
11、量增加井灌用水量,增加渠灌用水量减少井灌用水量和渠灌用水量与井灌用水量在合理灌溉比值的3 种情况下对地下水位动态进行了预测,并与多年平均地下水变化进行了分析比较,以说明 3 种情况对灌区地下水位动态的影响.1)减少渠灌用水量增加井灌用水量情况下地下水位动态变化表 1 不同水平年(P=50%.P=75%)下减少渠灌用水量增加井灌用水量情况下地下水位动态变化状况渠灌与井灌用水量比值(Q)0.70.80.9(P=50%)预测水位/m11.678.698.46(P=75%)预测水位/ml1.829.618.76(P=50%)与多年平均水位差/m5.432.852.22(P=75%)与多年平均水位差/m
12、5.583.372.52在渠灌用水量与井灌用水量的比值 Q 分别为0.7,0.8,0.9 较小的情况下,利用预测模型公式(5)分别计算得出预测结果如表 l 所示.计算结果表明,在减少渠灌用水量增加井灌用水量情况下,不论是平水年(P=50%)还是偏枯水年(P=75%),地下水位均为下降趋势,并且渠灌与井灌比值越小,即井灌用水量越大,地下水位下降越大.如果过量开采势必会造成地下水位大副下降,可使灌区局域形成地下水大的漏斗,造成机井报废,农作物大幅减产,给灌区水源工程和农民的经济均造成巨大损失.2)增加渠灌用水量和减少井灌用水量情况下地下水位动态变化使渠灌用水量与井灌用水量的比值 Q 分别为表 2
13、不同水平年 (P=5o%,P=75%)下增加渠灌用水量和1.75,2.0,2.25 较大的情况下,利用预测模型公式兰!兰竺!兰!竺查兰兰(5)笙果 2.所示 ._计算结果表明,在降雨频率为 P=50%和 P=P=75%;泵 4.28;75%的情况下,增加渠灌用水量(即减少井灌用水j 葫冱了币量),地下水位均为上升,上升幅度随渠灌与井灌用兰型生里 I:!I:兰!:水量比值的增大而增大.因灌区引用的泾河水源含有一定量的泥沙,加之渠灌用水水价较井灌高,所以渠灌用水量与井灌用水量比值不宜太大.高含沙引水灌溉时,若泥沙处理不当,往往会使渠道严重淤积和沙化农田土壤.在灌区地下水埋深较小的情况下,过多的引用
14、地表水灌溉易形成灌区土壤盐渍化问题,又不利于防涝.3)渠灌与井灌用水量在合理灌溉比值情况下地下水位动态变化使渠灌用水量与井灌用水量的比值分别接近最佳灌溉比值的情况下,利用预测模型公式(5)分别计算得出预测结果如表 3 所示.表 3 不同水平年(P-:o,P=75%)下渠灌用水量与井灌用水量在合理灌溉比值情况下地下水位动态变化状况渠灌与井灌用水量比值(Q)1.21.31.41.51.6(P=50%)预测水位/m6.786.575.114.143.9“/(P=75%)预测水位/m7.527.066.395.384.46(P=50%)与多年平均水位差/mJ0.54J0.33f1.13f2.10f2.
15、27(P=75%)与多年平均水位差/m1.280.820.15f0.86f1.78计算结果表明,渠灌用水量与井灌用水量在适宜比值情况下,灌溉后地下水位基本接近多年平均水位,或稍有回升.在降水平水年(P=50%)的最佳渠井灌溉用水比值在i.35 左右,在偏枯水年(P=75%)渠井灌溉用水最佳比是 1.45 左右,二者最佳比值存在差异.考虑到渠首来水的季节变化和田间灌溉需求,以及灌区引用的泾河水源含有一定量的泥沙,渠灌用水水价较井灌取水成本高等诸多问题,应该提倡:冬春灌以渠灌为主 ,夏灌以井灌为主,渠井结合,以渠养井,以井补渠,丰储枯用,达到地下水采补平衡,使灌区水资源可持续利用和生态环境良性循环
16、.4 结论1)利用灌区多年实测资料采用非线性回归分析方法,建立了井渠结合灌区地下水动态的预测模型,通过验证分析,表明预测模型是可靠的.2)在不同灌溉水平下的预测结果表明:在降水平水年 (P=50%)的最佳渠井灌溉用水比值在 i.35 左右,在偏枯水年(P=75%)渠井灌溉用水最佳比是 1.45 左右.灌区内冬春灌应以渠灌为主,夏灌以井灌为主,全年渠灌用水量与井灌用水量的比值为 1.5 左右为宜,做到渠井结合,以渠养井,以井补渠,丰储枯用,采补平衡,使灌区水资源得到合理开发利用.3)降雨量和灌溉用水量大小对灌区地下水动态变化都有着明显的影响,可根据不同水文年来适宜调控渠灌与井灌用水量的大小及比值
17、.4)如果过度开采地下水,会使地下水位,供水保证率均下降,导致部分机井报废,农作物减产,给灌区社会经济和生态环境造成严重损失.参考文献:1陈葆仁,洪再吉 ,等.地下水动态及其预测M.北京:科学出版社.1988.237238.2周维博.李佩成 .井渠结合灌区节水灌溉的有效途径J.沈阳农业大学,2004,35(5-6):4734753张蔚榛.农业节水问题的几点认识J.中国水利,2001,(8):4042.TheGroundwaterLevelForecastandtheRationalRatioAnalysisforIrrigationDistrictwithbothWelIsandCanalsZ
18、HOUWeibo.ZENGFaehen(EnvironmentalScienceandEngineeringCollege,ChanganUniversity,Xian710054,China)Abstract:Usingthenonlinearrelatedmethod,thispaperadoptsflmathematicmodeofforecastingthelevelofundergroundwaterinirrigationdistrictwhichtakesjinghuiqu.irrigationdistrictofshanxiprovinceasanexampleandisbas
19、edonthedataoftheyearsrainfalls,theconsumingamountofirrigationbycanalandwellandtherationalratioofthemintheirrigationdistrict.ThepaperforecaststhedynamicoftheundergroundwaterinthethreesituationsflSfollows:first,addingtheconsumingamountofwaterbywellwhiledecreasingconsumingamountofwateibycanal;second,de
20、creasingtheconsumingamountofwaterbywellwhileaddingconsumingamountofwaterbycanal;third,theproperratioofcanalwellwater,andproposestherationalratiocanal-wellwater.TheresultsCanbeusedasthetoolforthereasonableexploitationofundergroundwaterinirrigationdistrictsandtheoptimizeddispositionofthewaterresourceforirrigation.Keywords:irrigationdistrictwithbothwellsandcanals;thegroundwaterlevelforecast;nonlinearmodel;theratioofcanalwellwater9
