摘要
大沽河是青岛市主要水源地之一,是青岛市生活和工农业用水的重要水源。大沽河洪水大部分发生在7-8月,历时一般为3-5d。干流洪水对地下水的渗漏补给直接影响着行洪流量和地下水补给量的大小。研究大沽河河道地下水的补排规律,对建立流量(特别是洪峰流量)与地下水渗漏量的时空关系,为大沽河流域洪水库河联合调度系统的建立提供科学的依据,也为其它季节河流域调度的研究提供借鉴。
本研究对大沽河沿岸进行了气象、水文、地质等方面的调查工作,在了解和掌握国内外有关河道渗漏研究成果的基础上,采用试验分析和数值模拟相结合的方法,探讨在不同条件下河道的入渗规律及其影响因素。本文选取了大沽河沿岸上海坝和移风店的砂样作为研究对象,测定了砂样非饱和渗流特性;通过室内土柱、土槽试验和野外原位试验,模拟了河道水分的非饱和入渗;建立了不同条件下的包气带水分数值模型,对室内试验的结果进行模拟验证;利用VS2DT软件分析了水分非饱和运移时含水率和入渗通量随时间变化的特点和规律。通过上述一系列的研究,得出了一些认识和结论,主要成果如下:
(1)选取了大沽河沿岸上海坝和移风店的砂样作为研究对象,通过试验获得了砂样的水分特征曲线和水力传导度,利用Van Genuchten模型进行了参数的拟合研究,并对参数拟合结果进行了分析;
(2)通过土柱、土槽模拟试验,研究了在定水头入渗的条件下,不同质地、不同初始含水率的含水介质中,水分饱和-非饱和运移的规律和入渗通量随时间的变化,对孔隙中水分的迁移规律进行了探讨。通过试验数据的整理和比较,可知砂样的垂向和侧向入渗,通量随时间的变化可分为2个阶段,最初是一个较大的值,随时间迅速减小至一个较为稳定的值;然后随时间的变化趋于稳定。粗砂样在初始含水率为0.44%和9.30%时垂向稳定入渗通量约为0.003ml/(s*cm2)和0.002ml/(s*cm2)不同,细砂样在初始含水率为2.50%和9.65%时垂向稳定入渗通量约为0.002ml/(s*cm2)和0.0018ml/(s*cm2);粗砂样在初始含水率为2.60%、7.35%和13.98%时侧向稳定入渗通量约为分别约为0.002ml/(s*cm~2)、0.0011 ml/(s*cm2)、0.0008ml/(s*cm2),细砂样在初始含水率为3.20%、7.61%和16.20%时侧向稳定入渗通量约为分别约为0.0026ml/(s*cm2)、0.002ml/(s*cm2)、0.0016ml/(s*cm2)。
(3)在河道上中下游选取代表性地段进行非饱和入渗试验,通过野外原位试验,测定河道上中游不同位置的入渗通量规律。三次原位试验测得于家小里稳定通量为0.028ml/(s*cm2);王璧稳定通量为0.010ml/(s*cm2);南沙梁村稳定通量为0.023ml/(s*cm2)。河道上游的砂土入渗通量比中游的大;同为上游,距离河道越近,包气带渗透性越好,入渗通量越大。
(4)通过试验确定的参数、上下边界条件及初始条件等,在遵循大沽河河道特征的条件下,利用VS2DT软件建立了包气带水分运移数值模型,模拟了包气带水分的运移,验证了模型所采用的参数可以用来计算大沽河河道非饱和入渗;
(5)利用VS2DT软件,对包气带水分运移的理论模型进行了初步应用研究,并模拟了初始含水率和入渗水头对垂向/侧向入渗通量大小随时间变化的影响。
Daguhe River is one of the main source of water and an important sources of industrial and agricultural in Qingdao. The floods of Daguhe River usually occurred from July to August, which lasted for 3 to 5 day.The seepage of riverway to groundwater directly impacts the flux of floods and the supply of groundwater. Take a test to the seepage from Daguhe River can provide scientific basis to build a system of flow(particularly peak flow) and groundwater seepage, as well as offer a reference to other research of scheduling season river.
In this study, we carried out the meteorology, hydrology, geology and other aspects of the investigation along the Daguhe River. Take both experimental study and numerical simulation to explore the rule and effect of seepage in different riverway. We made Shangshai dam and Yifengdian dam as objects, mensurated the unsaturated characteristics of samples; simulated the unsaturated infiltration with laboratory studies and field studies; established the numerical model under different conditions of vadose zone to validate the results of the laboratory tests; analyzed the characteristic and rule in unsaturated infiltration with the software VS2DT. As a result of these studies, we obtained some new understanding and conclusions as follows:
(1) We selected the Shangshai dam and Yifengdian dam as studying objects,tested the curves of soil water characteristic and hydraulic conductivities of the samples, fitted the parameters and taked analyses by Van Genuchten model.
(2) We studied the saturated-unsaturated infiltration and changes of flux with time under the fixed head, through the soil column and trough in the laboratory studies, discussed the rule of moisture transfer. Through the data’s collation and comparison,we know that the moisture infiltrations not only had sameness but also had their own characteristics, which were related with the way of infiltration, the river rock, the moisture content of vadose zone, the time of infiltration.The vertical flux of coarse sand were 0.003ml/(s*cm2) and 0.002ml/(s*cm2) when the intial moisture contain was 0.44% and 9.30%,and the vertical flux of fine sand were 0.002ml/(s*cm2) and 0.0018ml/(s*cm2) when the intial moisture contain was 2.50% and 9.65%. The side flux of coarse sand were 0.002ml/(s*cm2),0.0011ml/(s*cm2) and 0.0008 ml/(s*cm2), when the intial moisture contain was 2.60%,7.35% and 13.98%,and the side flux of fine sand were 0.0026ml/(s*cm2),0.002ml/(s*cm2) and 0.0016ml/(s*cm2) when the intial moisture contain was 3.20%,7.61% and 16.20%.
(3) Through field studies,we determined the change of flux of seepage. According to the laboratory test results.We tested that the steady flux of Yujiaxiaoli was 0.028ml/(s*cm2), the flux of Wangbi was 0.010ml/(s*cm2) and the flux of Nashaliang was 0.023ml/(s*cm2).
(4) By using the parameters, upper and lower boundary conditions, and initial conditions, which were determined by studies, we used the software VS2DT to establish a vadose zone numerical model to simulate the moisture infiltrations, got the validation that the parameters could be used to calculate the unsaturated Seepage from the Daguhe River.
(5) We used the software VS2DT to take the vadose zone moisture infiltrations theoretical model in a preliminary applied research. We established a generalized model and analyzed the effects of intial moisture contain and the hand to the vertical and flux.
引文
[1] Darcy.H, Les fontoines publiques .de la Ville de Dijon.Dalmont, Paris.1856
[2] Richards, L.A., Capillary conduction of liquids through porous mediums, Physics,1931,318-333
[3]雷志栋,胡和平,杨诗秀.土壤水研究进展与评述.水科学进展,1999,NO.3
[4] Green,W.H. and Ampt, G..A., studies on soil physics: I .Flow of air and water through soils, Agr. Sci. 1911, 4.1~24
[5] Philip, J.R., The theory of infiltration:4., Sorptivity and algebraic infiltration equations, Soil Sci., 1957c, 84,257~264
[6] Aboukhled A,Alfaro,M Smith,Lysimeters Food and Agriculture organization of the united Nations FAO Irrigation and Drainage Paper 1981,39
[7] F.J.Molz,Model of water transport in the soil-plant system:A review, Water Resour Res,17,1981
[8] Dyer,J.A and A.R.Mack,The versatile soil moisture budget,Tech,Bull,87, Agromeorology section,Research Branch,1984
[9] Allison,G.B,A review of some of the physical,chemical and isotopic techniques available for estimation groundwater recharge,Reidel,Boston,Massachusetts,USA, 1988
[10] Van der lee, J.Q Gehrels, J.C,Modelling Aquifer Recharge, Introduction to the lumped parameter Model EARTH,Inst,Earth Sci,1990
[11] Shigeru ogawa Yosake Kishi, Studies on the infiltration-discharge of rain water and translation phenomena in soil, Journal of Hydrology,1992
[12] Jinquan Wu,R Zhang,J,Yang,Analysis of rainfall-recharge relationships,Journal of Hydrology,1996,177:143-160
[13] Jackson,T.J,J.Schmugge,Remote sensing Application to Hydrology-Soil Moisture,Hydrol sci.J.41, 1996
[14] James.BONTA and Mike Moller,Evaluation of the Glugla method for estimation evapotranspiration and ground water recharge,Hydrol sci J.44,1999
[15]庄季屏.四十年来的中国土壤水分研究.土壤学报,1989,26(3):241~247
[16]雷志栋,杨诗秀,等.土壤水动力学.北京:清华大学出版社,1988.276P
[17]张蔚榛等.地下水与土壤水动力学.北京:中国水利水电出版社,1996.287P
[18]康绍忠,刘晓明,等.土壤-植物-大气连续体水分传输理论及其应用.北京:水利电力出版社,1994
[19]李韵珠,李保国.土壤溶质运移.北京:科学出版社,1998.312P
[20]杨邦杰,隋红建.土壤水热运动模型及其应用.北京:中国科学技术出版社,1997.224P
[21]荆恩春等.土壤水分通量法实验研究.北京:地震出版社,1994.183P
[22]萧树铁.面入渗水分在包气带的运动.水文地质工程地质.1981,NO.5
[23]雷志栋,杨诗秀.非饱和土壤水一维流动的数值计算.土壤学报,1982,NO.2
[24]杨诗秀,雷志栋.均质土壤降雨喷洒入渗模型的数值计算.水利学报,1983,NO.5
[25]张蔚榛、张瑜芳.垂向一维均质土壤水分运动的数值模拟在降雨入诊和蒸发条件下的运用.工程勘探,1984,NO.5
[26]吴金全.降雨入渗条件下土壤水分运动及降雨入渗补给的研究.武汉水利电力大学研究生论文,1986
[27]解光城.雨入渗条件下层状非均质介质中的土壤水分运动及降雨对浅层地下水补给的初步研究.武汉水利电力大学研究生论文,1987
[28]任理.野外条件下土壤水从浅层地下水动态的模拟和预报.武汉水利电力大学研究生论文,1988
[29]雷志栋等.田间土壤水量平衡与定位通量法的应用.水利学报,1988,NO.5
[30]吴金全.地表水、土壤水及地下水相互转化关系的研究.武汉水利电力入学博士论文,1989
[31]邱景唐.非饱和土壤水零通量面的研究.水利学报,1992,NO.5
[32]张德桢,徐世民.大气水-地下水一土壤水-地下水相互转化关系的试验研究.水文地质工程地质,1993,NO.5
[33]高斯科,等.入渗条件下非饱和土壤水分运动的数值分析.西北水资源与水工程,1995.NO.4
[34]龚元石,李保间.应用农田水量平衡模型估算土壤水渗漏量.水科学进展,1995,NO.1
[35]刘群昌,谢森传.华北地区夏玉米田间分转化规律研究.水利学报,1998,NO.1
[36]杨建锋,地下水-土壤水-大气水界面水分转化研究综述,水科学进展,1998,NO.2
[37] Philip, J.R., 1957. The theory of infiltration: 2. The profile at infinity. Soil Science, 83:435-448
[38]杨诗秀,雷志栋,谢森传.匀质土壤一维非饱和流动通用程序.土壤学报,1985,(1):24~34
[39]周维博.降雨入渗和蒸发条件下野外层状土壤水分运动的数值模拟.水利学报,1991,(9):32~36
[40]彭建萍,邵爱军,段生贵.田间非饱和土壤水分运动的数值模拟.河北地质学院学报,Vol.18 No.6 1995
[41]韩雄,张强.低定额灌溉条件下农田土壤水蒸散模型的研究.安华北水利水电学院学报Vol.21,No.1,2000年3月
[42]陈喜.引黄灌区盐碱化预测模型.河海大学学报,1997年11月Vol.25,No.6,Nov.1997
[43]黄兴法,曾德超.土壤水热盐运动模型的建立与初步验证. Transactions of the CSAE 1997年9月
[44]骆祖江,付延玲,等.非饱和带水气二相渗流动力学模型.煤田地质与勘探,1999年10月第27卷第5期
[45]木拉提·胡塞因,虎胆·吐马尔拜.土壤水分运动数学模型的建立及应用.新疆农业大学学报,2002,25(1):60-62
[46]中华人民共和国国家技术监督局. GB3100~3102.中华人民共和国国家标准-量与单位.北京:中国标准出版社,1994-11-01
[47]邵爱军,沈荣开等.一种确定原状土样非饱和水力传导度的方法.土壤学报,2000年8月第三期,419-423
[48]郑西来,程善福等.滨海地下水库利用与保护.地质出版社:2007.1
[49]王铁军,郑西来,崔俊芳.莱西地区施肥对地下水硝酸盐污染的过程.中国海洋大学学报,2006年3月,第36卷第2期
[50]王铁军,郑西来,赵淑梅.土壤非饱和带确定性数值模型研究进展,海洋湖沼通报,2005年第2期
[51]刘贯群,朱新军等.孪井灌区主要作物节水灌溉模式的研究.中国海洋大学学报,2006年9月,第36卷第5期
[52]朱新军.内蒙古孪井灌区土壤水分运移及节水灌溉模式的研究:[硕士学位论文].中国海洋大学,2004
[53]张红梅.饱和-非饱和土壤中氟运移规律动态试验及数值模拟研究:[博士学位论文].河海大学,2005
[54]乔冈.天山北麓平原区包气带水分运移机理与数值分析:[硕士学位论文].长安大学,2006
[55]杨庆.银川平原包气带试验及数值模拟研究:[硕士学位论文].中国地质大学,2006
[56]张建丰.黄土区层状土入渗特性及其指流的实验研究:[博士学位论文].西北农林科技大学,2004
[57]张颀.灌区农田土壤水分运移模型研究与预测:[硕士学位论文].新疆农业大学,2005
[58] Yeh G T.FEMEATER:A finite-element model of water flow through saturated-unsaturated porous media-first revision. Environment Science Division Publication,No.2943.ORNL-5567/R1,USA,1987
[59] Brooks R H,Corey A T.Hydraulic properties of porous media. Colorada:Colorada State Univ,Fort Colins,1964
[60] Mualem Y.A new model for predicting the hydraulic conductivity of unsaturated porous media.Water Resour Res,1976.12:593-622
[61] Van Genuchten M Th.A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J,1980,44:892-898
[62] Kosugi K.Lognormal distribution model for unsaturated soil hydraulic properties.Water Resour Res,1996,32(9),2697-2703
[63] Stormont J C and zhou Shenxiong. Improving Pavement Sub-surface Drainage Systems by Considering Unsaturated Water Flow. Department of Civil Engineering, University of New Mexico, Aibuquerque, New Mexico, 2001: 1-74
[64]刘建立,徐绍辉,刘慧.估计土壤水分特征曲线的间接方法研究进展.水利学报,2004年2月第2期:68-78
[65]李少龙,杨金忠,蔡树英.基于Van Genuchten-Mualem模型的饱和非饱和介质流动随机数值分析.水利学报,2006年1月37卷第1期:33-39
[66]黄强,刘玉芸,李生秀,宋郁东.塔克拉玛干沙漠土壤水分运动参数的计算.干旱区地理,25卷第1期,2002年3月:75-78
[67]卢小慧,靳孟贵,汪丙国.栾城农业生态系统试验站土壤水分特征曲线分析.中国农村水利水电,2006年第l2期:30-32
[68]张喜英,张橹,刘昌明.太行山前平原土壤水分特征曲线拟合参数的确定.华北农学报2001年第16(2)卷:75-82
[69]张富仓,张-平.用-步出流法测定土壤的导水参数.西北农业大学学报,1993,21(增):129-133
[70]刘建立,徐绍辉,刘慧.估计土壤水分特征曲线的间接方法研究进展.水利学报,2004年2月第2期:68-78
[71]赵鸿铎,马网兔,吴青峰.公路中央分隔带渗水规律研究.公路交通科技,2005,22(3):47-50
[72]张华,陈善雄,陈守义.非饱和土入渗的数值模拟,岩土力学,2003年l0月第24卷第5期,715-718
[73]王红闪,黄明斌.四种方法推求土壤导水参数的差别与准确性研究,干旱地区农业研究,2004年6月第22卷第2期,76-80
[74]侯宪东,汪志荣.非饱和土壤水分运动数值模拟研究综述,水资源与水工程学报,2006年8月第17卷第4期,41-49
[75]张耀峰,张德生,武新乾.一维非饱和土壤水分运动的数值模拟.纺织高校基础科学学报,2004,17(2),123-127
[76]张宏仁.有限差分法的改进.水文地质工程地质,1994,2:27-30
[77]张思聪,惠示博,雷志栋.水平非饱和土壤水二维运动的准解析解.水利学报,1986,17(3):54-60
[78]寰镒吾.水平非饱和土壤水二维运动的近似解析解.四川工业学院学报,1990,9(12)
[79]齐春英,刘克岩.沿程渗漏河道的洪水流量演算模型.水文,1997,6
[80]孙淑珍.平原河道渗漏实验模型研究.河北工程技术高等专科学校学报,2003,6(2)