摘要：用最大30分钟雨强(I30)、径流量(Q)或者坡度(S)建立侵蚀量(Qs)的单因子或多因子方程。
用内蒙古自治区伊克昭盟五分地沟、五不进沟及河北省张家口市的坡度小区观测资料进行计算，比较其效果。
结果表明：Qs=kQmSn和Qs=kQm(坡度一定时)用来计算次降雨侵蚀量较好；用I30代替以上方程中Q的结果不理想。
关键词：雨强径流量侵蚀量坡度影响侵蚀量的因素很多，如降雨情况、地形(坡度、坡长、坡形)、地面状况(植被、土壤性质)等。
在建立侵蚀量的方程时，常用的变量是降雨强度、坡度、坡长、植被覆盖度、径流量等。
有用单因素的[1，2］，有用双因子的以至多因素的[3，4]。
多引进变量一般能提高预测精度，但资料的收集也更为困难。
以内蒙古伊克昭盟五分地沟、五不进沟及河北张家口坡度小区的观测资料为基础，本文旨在选择一种较好的计算次降雨侵蚀量的方法。
1 试验区基本情况张家口试验小区位于张家口市郊沈家屯镇马家沟流域郭家梁试验场西南坡耕地上，东经114°50′，北纬40°47′，海拔822m，土壤为黄土。
张家口属温带大陆性季风气候，多年平均降雨量400mm，其中80%～90%集中在7～9月份。
试验场内共设7个试验小区，由坡度9°的耕地通过简单的填挖方，改为坡度试验小区，坡度分别为0°、5°、10°、15°、20°、25°、30°。
小区面积均为10m2。
小区边界用混凝土板围成，下部有集水池。
小区常年休耕，耕层土质为粉质沙壤土，各小区的土壤粒度组成以及有机质含量见表1。
表1 张家口试验小区土壤性质分析Soil properties of runoff plots at Zhangjiakou gully粒度组成(%) 有机质含量小区坡度砂粒粉粒粘粒(%)2 5°37.550.0 12.5 0.53 3 10°44.845.2 11.0 0.57 4 15°41.5 49.0 9.5 0.66 5 20°38.0 52.5 9.5 0.5525°32.5 55.7 11.8 0.43 7 30°44.5 49.0 6.50.38五分地沟和五不进沟位于内蒙古自治区准格尔旗境内。
准格尔旗属于中温带大陆性气候，年平均气温6.2～8.7℃，从西北向东南逐步升高。
全年降水少而集中，多集中在7～9月。
降雨年际变化大，最低年份仅143.5mm，最高年份为636.5mm，平均400mm。
分地沟小区坡度分别为6°、9°、12°、15°，小区面积50m2。
土壤为黄土，地表无植被覆盖。
五不进沟小区坡度为9°、12°、15°、17°、20°，小区面积50m2。
地表物质为砒砂岩，无植被覆盖。
天然降雨观测是每次降雨后，测集水池中水位，取水沙样。
然后将水沙样过滤，烘干，称重。
最后计算含沙量、径流量、侵蚀量。
2 观测结果2.1 雨强在地表状况变化不大时，降雨因素就成为侵蚀量的决定因子。
在黄土高原，次降雨降雨量与侵蚀量相关性差，与降雨强度则存在很好的相关关系。
陈永宗、王万忠、蔡强国等发现，黄土高原坡面次降雨流失量与最大30分钟降雨强度的相关性最好[5～7]。
因此，在有雨强资料的五不进沟、五分地沟试验地，用每个径流小区的次降雨最大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 gully方程坡度主程R2次数显著性Qs=kIm309°Qs=0.0018I3.2130 0.6713非常显著12°Qs=0.0097I2.7030 0.5112非常显著15°Qs=0.0108I2.7730 0.2712不显著17°Qs=0.0009I3.69300.7713非常显著20°Qs=0.005I3.9730 0.6313非常显著Qs=kI30+c9°Qs=1.27I30-6.66 0.7313非常显著12°Qs=2.31I30-11.59 0.4112不显著15°Qs=3.77I30-19.37 0.6912非常显著17°Qs=4.89I30-28.81 0.2913不显著20°Qs=19.86I30-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=kIm306°Qs=0.16I1.3430 0.619显著9°Qs=0.22I1.6130 0.929非常显著12°Qs=0.86I0.9630 0.59显著15°Qs=0.95I1.2930 0.729非常显著Qs=kI30+c6°Qs=0.86I30-4.520.649显著9°Qs=1.75I30-7.64 0.929非常显著12°Qs=1.00I30-1.92 0.719非常显著15°Qs=3.41I30-13.95 0.819非常显著1 2 3 4下一页浅谈城市排水工程施工3 讨论3.1 方程的选择方程选择的标准主要是回归方程的R２值。
3.1.1 单因子方程指坡度或坡长等因素一定时，所建立的侵蚀量与某一变量的方程，本文指就每个径流小区建立的方程。
将各回归方程的R2按各地计算均值及均方差，结果见表6。
表6 R2均值及均方差Average R2 and square roots of their variances地点方程R2平均值R2均方差Qs=kQm 0.77 0.13五不进沟Qs=kIm 0.57 0.19Qs=kI+c 0.580.22Qs=kQm 0.70 0.07五分地沟Qs=kIm 0.71 0.15Qs=kI+c0.770.12张家口Qs=kQm0.690.11从表6及表2、表3可知，方程Qs=kIm30和Qs=kI30+c结果则从不显著到非常显著都有，R2值变化很大，前者R2值小一些，但均方差也小。
而方程Qs=kQm各小区回归结果都非常显著(表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 R2=0.84 n=50图1 五不进沟lgQ s实测及计算值Observed and calculated values of soil loss in Wubujingou gully图2 五分地沟lgQ s实测及计算值图3 张家口lgQ s实测及计算值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=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值，且前者具有随坡度增大的趋势，后者与坡度之间关系不明显，也表明I30不能代替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 野外模拟试验1001.140.24 Bissonais等室内模拟试验300.98～1.28Mathier等野外模拟试验2～151.50.9Band野外观测62.070.84 Averaert室内模拟试验601.10～3.101.77～2.96 本文五不进沟17.721.771.58五分地沟20.920.721.72张家口-0.800.57注：五不进沟、五分地沟雨强用两倍最大30分钟雨强。