《.次降雨侵蚀量的计算》由会员分享,可在线阅读,更多相关《.次降雨侵蚀量的计算(40页珍藏版)》请在三一办公上搜索。
1、欧洲自由行 摘要:用最大30分钟雨强(I30)、径流量(Q)或者坡度(S)建立侵蚀量(Qs)的单因子或多因子方程。用内蒙古自治区伊克昭盟五分地沟、五不进沟及河北省张家口市的坡度小区观测资料进行计算,比较其效果。结果表明:Qs=kQmSn和Qs=kQm(坡度一定时)用来计算次降雨侵蚀量较好;用I30代替以上方程中Q的结果不理想。 关键词:雨强 径流量 侵蚀量 坡度 影响侵蚀量的因素很多,如降雨情况、地形(坡度、坡长、坡形)、地面状况(植被、土壤性质)等。在建立侵蚀量的方程时,常用的变量是降雨强度、坡度、坡长、植被覆盖度、径流量等。有用单因素的1,2,有用双因子的以至多因素的3,4。多引进变量一般
2、能提高预测精度,但资料的收集也更为困难。以内蒙古伊克昭盟五分地沟、五不进沟及河北张家口坡度小区的观测资料为基础,本文旨在选择一种较好的计算次降雨侵蚀量的方法。1 试验区基本情况张家口试验小区位于张家口市郊沈家屯镇马家沟流域郭家梁试验场西南坡耕地上,东经11450,北纬4047,海拔822m,土壤为黄土。张家口属温带大陆性季风气候,多年平均降雨量400mm,其中80%90%集中在79月份。试验场内共设7个试验小区,由坡度9的耕地通过简单的填挖方,改为坡度试验小区,坡度分别为0、5、10、15、20、25、30。小区面积均为10m2。小区边界用混凝土板围成,下部有集水池。小区常年休耕,耕层土质为粉
3、质沙壤土,各小区的土壤粒度组成以及有机质含量见表1。表1 张家口试验小区土壤性质分析Soil properties of runoff plots at Zhangjiakou gully 粒度组成(%) 有机质含量小区坡度砂粒粉粒粘粒(%)2537.550.012.50.5331044.845.211.00.5741541.549.09.50.6652038.052.59.50.5562532.555.711.80.4373044.549.06.50.38五分地沟和五不进沟位于内蒙古自治区准格尔旗境内。准格尔旗属于中温带大陆性气候,年平均气温6.28.7,从西北向东南逐步升高。全年降水少而集
4、中,多集中在79月。降雨年际变化大,最低年份仅143.5mm,最高年份为636.5mm,平均400mm。分地沟小区坡度分别为6、9、12、15,小区面积50m2。土壤为黄土,地表无植被覆盖。五不进沟小区坡度为9、12、15、17、20,小区面积50m2。地表物质为砒砂岩,无植被覆盖。天然降雨观测是每次降雨后,测集水池中水位,取水沙样。然后将水沙样过滤,烘干,称重。最后计算含沙量、径流量、侵蚀量。2 观测结果2.1 雨强在地表状况变化不大时,降雨因素就成为侵蚀量的决定因子。在黄土高原,次降雨降雨量与侵蚀量相关性差,与降雨强度则存在很好的相关关系。陈永宗、王万忠、蔡强国等发现,黄土高原坡面次降雨流
5、失量与最大30分钟降雨强度的相关性最好57。因此,在有雨强资料的五不进沟、五分地沟试验地,用每个径流小区的次降雨最大30min雨强与侵蚀量作回归运算,方程形式为Qs=kIm30(1)Qs=kI30+c(2)式中s是侵蚀量(kg),I30是最大30分钟雨强(mm/30min),k,m,c是常数。方程(1)是取对数后运算的,即lgQs=lgk+mlgI30以下幂函数方程都是这样计算的,方程(1)的结果见表2、表3。表2 五不进沟侵蚀量与最大30分钟雨强回归方程Equations for soil loss based on rainfall intensity (I30)in Wubujingou
6、 gully 方程坡度主程R2次数显著性Qs=kIm309Qs=0.0018I3.21300.6713非常显著 12Qs=0.0097I2.70300.5112非常显著 15Qs=0.0108I2.77300.2712不显著 17Qs=0.0009I3.69300.7713非常显著 20Qs=0.005I3.97300.6313非常显著Qs=kI30+c9Qs=1.27I30-6.660.7313非常显著 12Qs=2.31I30-11.590.4112不显著 15Qs=3.77I30-19.370.6912非常显著 17Qs=4.89I30-28.810.2913不显著 20Qs=19.86
7、I30-94.450.7813非常显著注:显著、不显著指在0.05水平检验的结果,非常显著指在0.01水平下显著,下同。方程(1)中m值与坡度的关系不明显,k值则随坡度增大。五不进沟17小区方程Qs=kIm30中k值异常,经计算,年侵蚀量17小区略小于15小区,可能是因为17小区土壤性质与其它小区有差别。表3 五分地沟侵蚀量与最大30分钟雨强回归方程Equations for soil loss based on rainfall intensity (I30) in Wubujingou gully 方程形式坡度方程R2次数检验结果Qs=kIm306Qs=0.16I1.34300.619显著
8、 9Qs=0.22I1.61300.929非常显著 12Qs=0.86I0.96300.599显著 15Qs=0.95I1.29300.729非常显著Qs=kI30+c6Qs=0.86I30-4.520.649显著 9Qs=1.75I30-7.640.929非常显著 12Qs=1.00I30-1.920.719非常显著 15Qs=3.41I30-13.950.819非常显著1 2 3 4下一页 浅谈城市排水工程施工 3 讨论3.1 方程的选择方程选择的标准主要是回归方程的R值。3.1.1 单因子方程指坡度或坡长等因素一定时,所建立的侵蚀量与某一变量的方程,本文指就每个径流小区建立的方程。将各回
9、归方程的R2按各地计算均值及均方差,结果见表6。表6 R2均值及均方差Average R2 and square roots of their variances 地点方程R2平均值R2均方差 Qs=kQm0.770.13五不进沟Qs=kIm0.570.19 Qs=kI+c0.580.22 Qs=kQm0.700.07五分地沟Qs=kIm0.710.15 Qs=kI+c0.770.12张家口Qs=kQm0.690.11从表6及表2、表3可知,方程Qs=kIm30和Qs=kI30+c结果则从不显著到非常显著都有,R2值变化很大,前者R2值小一些,但均方差也小。而方程Qs=kQm各小区回归结果都非
10、常显著(表4),R2值变化较小,平均约0.7(表6)。因此,用径流量计算侵蚀量最好。3.1.2 双因子以及三因子方程多因子回归方程的结果(表5)都非常显著,以R2为标准,方程(5)的结果差一些。方程(6)比方程(4)多引进一个变量,R2值略有增加,但计算量增大,且Q与I30本身相关。因此,方程(4)应优于(6)。将三地侵蚀量的对数的实测值与计算值(据方程(4)计算)作图于图1、图2、图3,实测值与计算值符合较好。从图1可以看到,五不进沟17小区多数点计算值大于实测值。前已述及,年侵蚀量17小区略小于15小区。将17小区作为异常值去掉,用其它小区资料得出 Qs=6760.83Q1.82S2.23
11、 R2=0.84 n=50图1 五不进沟lgQs实测及计算值Observed and calculated values of soil loss in Wubujingou gully 图2 五分地沟lgQs实测及计算值图3 张家口lgQs实测及计算值Observed and calculated values of soil loss at WufendigouObserved and calculated values of soil loss at Zhangjiakou图2中五分地沟12小区也有类似现象,用其它三个小区资料得Qs=863.38Q0.77S1.95 R2=0.81 n=
12、27作上述处理后,R2值明显增大。3.2 雨强与径流量的关系一些研究者用坡长(L)及雨强(I)两者组合起来,代替径流量。用LI代替Q时,其中含有假设:坡面土层不可渗透、或达到饱和、或全坡面均一入渗,从而径流量与LI呈线性关系。本文中五分地沟及五不进沟小区坡长L一定(均为10m),因此用I30代替Q。尽管小区当初是经过平整的、土壤性质较为一致的直形坡,但方程Qs=kIm30Sn与Qs=kQmSn、Qs=kIm30与Qs=kQm相比,后者结果好得多,说明用LI代替Q的前提条件不易满足。Qs=kIm30中m值明显大于Qs=kQm中m值,且前者具有随坡度增大的趋势,后者与坡度之间关系不明显,也表明I3
13、0不能代替Q。3.3 系数的大小将Qs=kQmSnIp30看作总方程,以上其它方程就是该方程的某种简化形式。由于方程之间有联系,上述方程系数之间有一些规律,如方程(1)、(2)k值随坡度增大。各地方程Qs=kQmSn与Qs=kQm中m值大小比较接近。Mathier发现方程s=kQm中m随坡度增大13,本文中m值与坡度没有明显的关系。由于多数研究都用方程(4),下面主要分析系数方程(4)中系数k、m、n大小及其变化。k表示地面粗糙度、颗粒粒度等因素的影响。由于各自所用单位不同,大小比较没有意义。表7 m、n值表Values of m and n 研究者来源雨强(mm/h)mnWainwright
14、野外模拟试验1001.140.24Bissonais等室内模拟试验300.981.28 Mathier等野外模拟试验2151.50.9Band野外观测62.070.84Averaert室内模拟试验601.103.101.772.96本文五不进沟17.721.771.58 五分地沟20.920.721.72 张家口-0.800.57注:五不进沟、五分地沟雨强用两倍最大30分钟雨强。张家口无雨强资料。*用方程(3)计算的结果。一般认为,m值表示水流的输送能力,而n值反映坡度在水流侵蚀中的作用。Julien和Simons总结出m取值范围1.42.4,n值范围1.21.94。不少研究者将该结果当作标准
15、,将所得结果与之比较(表7)。Wainwright解释其实验所得m、n较小的原因在于土壤的石质及坡面的不规则11。Band认为:在土壤粒径随坡度变化的地方,n值较小(多见于干旱气候下无植被覆盖的坡面);在土质相同时n值应该较大(植被较好发育,地表水流不足以造成地表物质沿坡面的分化常见)14。Everaert发现:m、n值随土壤粒径增加而增加15。与Julien和Simons的标准相比,本文中五不进沟m、n值在范围内,五分地沟m值、张家口m、n值偏小。由于观测都是在土壤质地相对均一、坡面较平整的小区进行的,以上因素(粒度、覆盖、坡形)对m、n值影响不大。对张家口数据分成每年计算,得出m、n值见表
16、8。表8 张家口小区各年计算m、nCalculated m and n for data collected during every year at plots of Zhangjiakou 年份kmnR2次数19922248068069078901993321881270810892419941290313102707059由表8可见,1993年、1994年m值接近Julien和Simons的标准,1992年m值则较小。Bissonnains认为Q=kQm中m受土壤湿度的影响10。Bryan通过实验发现,水流输送能力是坡度的多项式方程。方程系数变化很大,受土壤性质及试验小区尺寸的影响13。
17、可能因为这种复杂性,m值变化较大,张家口范围为0.681.31。因此,五分地沟和张家口计算所的m值也属正常。张家口n值各年变化较大,没有明显的趋势(表8),与Julien和Simons的标准相比较小。前面提到,方程Qs=kIm30中m值及Qs=kI30+c中k值随坡度增加而增大(个别小区k值异常)。因此,坡度间侵蚀量的差别随雨强的增大而增加,表现在方程Qs=kQmsn中为大雨强降雨次数多时n值增大。从表7可知,除Wainwright的研究是在不规则坡面进行的外,其它研究都表明,大雨强试验(观测)取得较大的n值。本文中五不进沟及五分地最大30分钟雨强平均值分别为8.86mm/30min、10.4
18、6mm/min。张家口没有雨强资料,但张家口观测次数多,如果大雨强降雨次数少时n值也会较小。张家口三年所得结果不同可能原因在于此。从以上分析来看,n值受到雨强的影响。4 结论综上所述,可得出以下结论:1.方程Qs=kQmSn和Qs=kQm用来计算侵蚀量较为可靠。2.坡度小区(坡长一定)用I30代替Q结果差一些,用I30代替Q时必须注意满足所需条件。3.方程Qs=kQmSn中m值影响因素复杂,其值变化较大,n值受雨强的影响。 1 2 3 4下一页 浅谈城市排水工程施工 参考文献1 吴普特,周佩华。地表坡度与薄层水流侵蚀关系的研究。水土保持通报,1993,13(3),P1-5.2 Liu,B.Y.
19、,Nearing,M.A.and Risse,L.M.Slope gradient effects on soil loss for steep slopes.Transaction of American Society of Agricultrue Engineers,1994,37(6),P1835-1840.3 江忠善,王志强,刘志。黄土丘陵区小流域土壤侵蚀空间变化定量研究。土壤侵蚀与水土保持学报,1996,2(1),P1-9.4 Julien,P.Y.,Simons,D.B.Soil loss capacity of overland flow.Trans.Am。Soc.Agric.
20、Eng.28,P755-766,1985.5 陈永宗,景可,蔡强国。黄土高原现代侵蚀与治理。科学出版社,1988.6 王万忠,焦菊英。黄土高原降雨侵蚀产沙与黄河输沙。科学出版社,1996.7 蔡强国,吴淑安等。坡耕地表土结皮对降雨径流和侵蚀产沙过程的影响。见:山西省水土保持科学研究所等编,晋西黄土高原土壤侵蚀规律实验研究文集,水利电力出版社,1990.8 蔡强国,马绍嘉等。黄土坡耕地上坡长对径流侵蚀产沙过程的影响。见:蔡强国等编,水土流失规律与坡地改良利用,环境科学出版社,1995.9 蔡强国,马绍嘉等。白洋淀平原地区不同地类径流与侵蚀产沙模拟降雨试验的初步研究。见:章申等著,白洋淀区域水污染
21、控制研究,科学出版社,1995.10 Le Bissionnais,Y.,Renaux,B.,Delouche,H.Interactions between soil properties and moisture content in crust formation,runoff and interrill erosion from tilled loess soils.Catena 25,P33-46,1995.11 Wainwright,J.Infiltration,runoff and erosion characteristics of agricultural land in ex
22、treme storm events,SE France.Catena 26,P27-47,1996.12 Carson,M.A.& Kirkby M.J. Hillslope form and process.Cambridge University Press,1972.13 Mathier,L.,Roy,A.G., & Pare,J.P.The effect of slope gradient and length on the parameters of a sediment transport equation for sheetwash.Catena 16,P545-558,1989.14 Band,L.Field parameterization of an empirical sheetwash transport euqation.Catena 12,281-290,1985.15 Everaert,W.Empirical relations for the sediment transport capacity of interrill flow.Earth surface process and landforms 16,P513-532,1991. 1 2 3 4下一页 浅谈城市排水工程施工