黄土区陡坡侵蚀过程试验研究

土壤与环境 2002, 11(4): 356359 Soil and Environmental Sciences E-mail: 基金项目： 中国科 学院知识创新重大项目 （ KZCX1-10-04） 作者简介： 郑良勇 （ 1977），男，硕士， 从事坡面侵蚀水动力学等方面的研究 。 E-mail: 收稿日期： 2002-04-10 黄土区陡坡侵蚀过程试验研究 郑良勇，李占斌，李 鹏 中国科学院、教育部、西北农林科技大学水土保持与生态环境研究中心，陕西 杨凌 712100 摘要： 通过室内放水冲刷试验，结果表明，陡坡单宽径流能耗与单宽径流产沙率之间存在以下线性关系： Dr=18.672（ E - 0.751），说明土壤的可蚀性参数为 18.672 g/J，发生坡面细沟侵蚀的临界径流能耗为 0.751 J/(ms)。坡面单宽径流能耗随流量增大而增加，随坡度变化呈抛物线趋势，临界坡度出现在 21和 24之间。坡面上各段面单宽径流能耗均随流量的 增大而增加，随着坡度逐渐增加，坡面上部段面能耗渐增，中部段面渐减，下部较稳定。坡面各段面侵蚀产沙也有类似特征。此研究对于深入了解黄土高原陡坡土壤侵蚀过程和机理具有重要意义。 关键词： 径流能耗；陡坡；产沙率 中图分类号： S157 文献标识码： A 文章编号： 1008-181X（ 2002） 04-0356-04黄土高原坡陡沟深，水土流失严重，尤其近年来陡坡开荒严重，更加剧了坡耕地侵蚀。细沟侵蚀是坡面侵蚀的主要形式，其侵蚀量约占坡面总侵蚀量的 70%1。因此研究陡坡细沟侵蚀对于防治坡耕地 水土流失具有重要意义。细沟是在坡面径流差异性侵蚀条件下，在坡面上产生的一种小沟槽地形，其纵剖面与所在斜坡面一致，并能为当年犁耕所平复 2。细沟侵蚀是降雨形成的坡面漫流在顺坡向下流动的过程中逐渐汇集成坡面股流并对坡面进行下切侵蚀，溯源侵蚀和侧蚀的过程。只有当坡面径流达到一定的水力临界后才可能发生细沟侵蚀，水力临界主要用径流量，佛汝德数，径流水流功率，径流剪切力等来描述。 Rauws 和 Govers 研究发现，细沟侵蚀的临界条件可用有效剪切流速与土壤饱和粘滞力之间的线性关系表示，且试验结果表明，35 cm/s 的切 应流速是细沟发生的水力学临界 3。张科利则用佛汝德数 Fr 1 作为细沟发生的水动力临界。细沟水流对土壤的冲刷作用是几种作用的综合。 Elliot and Laflen 把细沟水流的分散能力分为水流冲刷、沟头冲刷、水流淘刷崩塌和沟壁脱落四部分 4。 Nearing 等认为细沟侵蚀是在水流的切应力大于土壤临界切应力且输沙能力大雨实际输沙率的条件下发生的，并提出以下模型： Dr=Dc(1-Gs/Tc) (1) Dc=k( - c) (2) 其中， Dr 为径流分离率； Dc 为径流分离能力；Gs 为输移率； Tc 为输移能力； k 为与土壤性能有关的土壤可蚀性参数； 为径流切应力； c 为临界切应力 5。 尽管目前对坡面细沟侵蚀的研究较多且取得丰硕成果，但上述研究一方面大多是在缓坡进行，难以直接适用于我国坡陡沟深的黄土高原，另一方面细沟发育是一个动态变化过程，细沟水流的水深、流速、切应力等水力要素时空变化很大，以往研究多用整个坡面的平均值来计算，这将不能确切反映细沟发育的实际过程，因此本文通过室内放水冲刷模拟试验利用径流能耗理论对坡面细沟侵蚀过程作 一新的探讨。 1 材料与方法 试验在中科院水土保持研究所降雨大厅利用可调坡度钢制土槽进行，土槽长 5 m、宽 0.33 m、深 0.5 m。试验流量按黄土高原暴雨发生频率在野外标准径流小区上产生的单宽流量换算到试验土槽上的流量得到，分别为 2.5 L/min、 3.5 L/min、 4.5 L/min、 5.5 L/min、 6.5 L/min 五个流量，坡度采用15、 18、 21、 24、 27、 30六级变化。试验土壤体积质量控制在 1.25 g/cm3 左右，试验前土槽洒水使土壤表层达到饱和。每次实验用颜料示踪法分四个段面每分钟测量一次流速，沿土槽纵向依次确定12 个点，每分钟在这些固定点位置用直尺测量水宽，水深由流量、流速和水宽计算获得。每次实验前后用温度计在稳流槽中测水温。在土槽出水口每隔一分钟接取泥沙样以烘干获得含沙量，接取全部径流泥沙样以计算侵蚀量。试验进行两个重复。试郑良勇等：黄土区陡坡侵蚀过程试验研究 357 验填土为杨凌当地娄土，其粒径组成如表 1 所示。 2 结果与讨论 2.1 陡坡细沟发生的能量分析 土壤侵蚀是降雨及降雨汇集的坡面水流对土壤作功的过程，因此必然遵循能量守恒定律。坡面径流在顺坡向下流动的过程中需要克服坡面摩擦力、水流内摩擦力作功，水流剪切分离土壤和携带输移土粒的过程实质上就是一个作功消耗能量的过程，对于一定的土壤来说，消耗相同的能量就会剥离相同的土粒。因此，我们可以通过分析坡面径流的能量变化运用径流能耗理论来研究坡面土壤侵蚀过程。 设单宽径流在坡面顶端所具有的势能为 E 势 = qgLsin ，动能为 E 动 =(1/2) qV12，其总能量为： E=E 势 +E 动 = qgLsin +(1/2) qV12 (3) 如果在理想状态下，径流到达坡面任何部位其总能量应该不变。但由于水流沿程要克服水土界面摩擦力和水流内摩擦力作功，一部分能量转化为热能散发掉，因此在坡面某一段面处径流实际的总能量与理想状态下有很大差异。设在坡面任一段面 x处的动能为 Ex 动 ，势能为 Ex 势 ，总能量为 Ex，则 Ex=Ex 动 +Ex 势 = qxgLxsin +(1/2) qxVx2 (4) 其中： q 为单宽流量（ m2/s）； 为坡度（ ）；L 为坡长（ m）； Lx 为坡面任一段面 x 到坡脚的距离（ m）； V 为坡面顶端流速（ m/s）； Vx 为段面处流速（ m/s）； g 为重力加速度常数（ 9.8 m/s2）。如设径流从坡面顶端到坡面上的任意段面 x 处消耗的能量为E 耗 ，则 E 耗 =E-Ex (5) 径流所消耗的能量一部分用于剥蚀分离土壤，另一部分消耗于携带输移土粒。因此坡面侵蚀产沙率与坡面径流能耗之间必然存在一定关系。如果我们假设单宽径流 在单位长度坡面上的产沙率为Drg/(ms)，单宽径流在坡面上流过单位长度距离消耗的能量为 EJ/(ms)，点绘试验中 Dr 与 E 的关系得到图 1。从图 1 中可以看出，不同坡度下单宽径流在单位长度坡面上的产沙率与单宽径流在坡面上流过单位长度距离消耗的能量之间存在很好的线性关系。我们用下式表示这种关系： Dr=Kr( E-E0) (6) 其中： Kr 为与土壤性质有关的可蚀性参数，不随坡度和流量变化； E0 为发生细沟侵蚀的临界径流能耗；（ E-E0）代表有效径流能耗， 即实际作用于剥离输移土壤的能耗。该式表明只有当径流的能耗大于临界径流能耗时，坡面才会产沙，且产沙率与有效径流能耗呈正比。 对上图各点趋势进行拟合，可得以下关系式： Dr=18.672( E-0.751) R 2 =0.932 (7) 上式表明在试验范围内，所选用土壤的可蚀性参数为 18.672 g/J，即在试验条件下单宽径流消耗 1 J 的能量在单位长度坡面上产沙为 18.672 g。发生坡面细沟侵蚀的临界径流能耗为 0.751 J/(ms)。利用（ 7）式我们可以在已知坡面径流能耗的情况下推算出坡面 的产沙率，这对我们了解细沟发生的时空变异特征具有重要意义。 从图 2 中可以看出，在相同的坡度下，流量越大，径流能耗越大，产沙率也愈大。这是由于在相同的坡度下，流量越大，径流总能量越大，对坡面的侵蚀加剧，消耗的能量也越多。在相同的流量下，径流能耗随坡度呈先增后减的抛物线趋势，其最大临界值出现在 21和 24之间。当流量相同时，坡度越陡，径流总能量越大，对坡面的侵蚀越剧烈，消耗的能量应该越多。但是另一方面，坡度越大，土壤沿坡面方向的重力分力越大，土粒稳定性越差，消耗较少的能量就可以剥离起动土粒，因此径流能耗反而减 少。由图中可知，当坡度由 15增加到 21时，径流能耗随之增加；当坡度大于 24时，径流能耗则随坡度呈递减趋势，即径流能耗变化的临界坡度在 21和 24之间。此临界坡度随着流量的变化略有差异，这种现象可能与坡度和流量对坡面径流能量的影响以及土壤抗冲性的综合作用有关。 表 1 试验土壤的粒径组成 粒径 /mm 1 0.25 0.25 0.05 0.01 0.005 0.05 0.01 0.005 0.001 0.001 百分比 /% 0.12 2.70 6.88 41.13 12.89 36.28 图 1 单 宽 径 流 产 沙 率 与 单 宽 径 流 能 耗 关 系0204060801001200 1 2 3 4 5 6 7 E / ( J m-1s-1)Dr/(gm-1S-1)151821242730358 土壤与环境 第 11 卷第 4 期（ 2002 年 11 月） 2.2 坡面侵蚀分异布特征分析 径流顺着坡面向坡下方流动的过程，对坡面土壤产生剪切分离作用，剥离输移土壤并消耗能量。在坡面不同位置由于径流对土壤的侵蚀作用各异，因此径流能耗在坡面上不同段面具有空间分异特征。另一方面，随着试验时间的延续，在坡面 上形成不同的侵蚀阶段，土壤侵蚀形态各异，消耗的能量也不同，所以坡面径流能耗也具时间分异特征。 在试验中将整个坡面平均分为 4 个段面，每个段面每一分钟分别量测其流速和水宽，由此我们可以分别计算各个段面的径流能耗。图 3 是不同坡度和流量条件下各段面径流能耗的变化图。 在图中横坐标为流量和坡度的综合指标，即2.5/15 代表坡度为 15流量为 2.5 L/min，依次类推；段面 1 在坡面下部与出水口相连，向坡上方 依次为段面 2、段面 3，坡面最上部为段面 4。从图中可以看出在相同的坡度下，随着流量的增加，各段面径流能耗也增大。这是因为随着流量的增加，在相同的坡度下径流总能量也增大，对坡面的侵蚀力加剧，侵蚀量增大，消耗的能量越多。坡度为 15时，段面 3 的能耗最大，段面 1 段面 2 其次，面 4 能耗最小。随着坡度的增加，段面 4 的能耗逐渐增大，段面的能耗 3 则逐渐减小。当坡度增加到 21时，段面 1 段面 2 和段面 3 的能耗基本相同。坡度大于24时，段面 3 的能耗继续降低。这说明坡度较小时，径流在第四段面初始能量较小，须经过一段距离加速到达段面 3 时，能量才达到最大值，此时对坡面的冲刷力最大，径流能耗也最大。随着坡度的增加，径流能量在段面 4 便达到最大，对坡面上部土壤产生剧烈侵蚀，消耗大量能量，产生大量泥沙；因此当径流到达段面 3 时，主要以输移泥沙为主，对坡面的剥蚀分离作用减弱，消耗的能量减少。坡度越陡，这种情况越明显。在坡面下部的 1、 2 段面，径流对土壤的分离输移作用比较稳定，随坡度变化不大。 根据上面所述的理论，在已知各段面径流能耗的情况下，可以根据拟合的公式来推算各段面产沙率，由此可以了解坡面产沙的时空分布特征、细沟的发生发育过程以及坡面剥 蚀、输移和沉积的相互转化规律。图 4 是根据各段面能耗利用拟合公式求算的坡面各段面产沙随坡度流量的变化图。 由图中可以看出，在不同坡度下，随着流量的增大，各段面的产沙率均增大。这是因为坡度相同时，流量越大，径流总能量越多，对坡面的剪切分力作用越强，产沙率也越大。在坡度为 15时，段面 3 产沙率最大，随着坡度的增加，段面 3 的产沙率逐渐降低，段面 4 的产沙率则渐增。另外，当坡度 15由增加到 18时，段面 1 和段面 2 的产沙率出现负值，说明此处泥沙沉积量大于剥蚀量，出现泥沙净沉积。当坡度为 24流量为 2.5 L/min 时，段面3 也开始出现净沉积。随着坡度的继续增大，段面3 的净沉积逐渐增多。由图中可知，坡面下部净沉积一般发生在每个坡度的较小流量时，并且坡度越大，发生净沉积的坡面越多。究其原因是因为流量较小时，径流能量较小，没有足够的能量将在坡面 图 3 各段面径流能耗随坡度流量变化图 图 2 单宽径流能耗随坡度变化 图 4 各段面侵蚀产沙变化图 051015202530352 . 5 / 1 5 5 . 5 / 1 5 3 . 5 / 1 8 6 . 5 / 1 8 4 . 5 / 2 1 2 . 5 / 2 4 5 . 5 / 2 4 3 . 5 / 2 7 6 . 5 / 2 7 4 . 5 / 3 0流 量 / 坡 度能耗E/(Jm-1s-1)段面 4段面 3段面 2段面 1-200 -100 0 100 200 300 400 500 600 2.5/15 6.5/15 5.5/18 4.5/21 3.5/24 2.5/27 6.5/27 5.5/30 流量 /坡度 产沙率/（gm-1 s-1) 断面 4 断面 3 断面 2 断面 1 0123456715 18 21 24 27 30坡度 / oE/(Jm-1s-1)2 . 5 L / m i n3 . 5 L / m i n4 . 5 L / m i n5 . 5 L / m i n6 . 5 L / m i n郑良勇等：黄土区陡坡侵蚀过程试验研究 359 上部剥蚀的土壤带出出水口，因此在坡面下部出现净沉积。随着坡度的增加，一方面径流在坡面上方具有的能量越多，对段面 4 的侵蚀越大，消耗的能量越多，产生的泥沙越多；另一 方面要把这些泥沙都带走需要更多的能量，这是一个矛盾统一的过程。在坡面上部产生的泥沙越多，消耗的能量就越多，因此径流用于输移泥沙和在下部段面剥蚀土壤的能量就越少，在下部段面出现泥沙净沉积。当泥沙部分沉积后，径流用于输移泥沙消耗的能量减少，从而有更多的能量用于剥离土粒，段面产沙率便相应增大。因此在坡度为 27和 30时，段面 3 的产沙率小于段面 1 和段面 2。 3 结论 （ 1）试验的研究结果表明坡面单宽径流能耗与单宽径流输沙率之间存在良好的线性关系，表达式为 Dr=18.672（ E-0.751），所选用土壤的可蚀性参数为 18.672 g/J，发生坡面细沟侵蚀的临界径流能耗为 0.751 J/(ms)。 （ 2）面单宽径流能耗随流量增大而增加，随坡度的变化抛物线趋势，临界极值坡度出现在 21和 24之间。 （ 3）坡面上各段面单宽径流能耗均随流量的增大而增加，随着坡度逐渐增加，坡面上部段面能耗渐增，中部段面渐减，下部较稳定。由此反映坡面各段面侵蚀产沙也有类似特征。 参考文献： 1 朱显谟 . 黄土高原水蚀的主要类型及其有关因素 J. 水土保持通报 , 1981(4): 13-18. 2 郑粉莉 , 高学田 . 黄土坡面土壤侵蚀过程与模 拟 M. 陕西人民出版社 , 2000: 8-9. 3 RAUWS G, GOVERS G. Hydraulic and soil mechanical aspects of rill generation on agricultural soilsJ. Soil Sciences, 1988, 39: 111-124. 4 ELLIOT W J, LAFLEN J M. A process-based rill erosion modelJ. Trans of the ASAE, 1993, 36: 65-72. 5 NEARING M A, FOSTER G R, LANE L J. A Process-based soil erosion model for USDA-Water Erosion Prediction Project technologyJ. Trans of ASAE, 1989, 32(5): 1 587-1 593. Study on the erosion process in steep slope on Loess slope ZHENG Liang-yong, LI Zhan-bin, LI Peng Research Center of Soil and Water Conservation and Ecological Construction, Chinese Academy of Sciences, Education Ministry； Northwest Sci-tech University of Agriculture & Forestry, Yangling, Shaanxi 712100, China Abstract: The results of runoff scouring experiments show that the relationship between unit width