坡面植被覆盖度对泥沙输移的影响特性

孙 一，刘晓燕，田 勇，伊晓燕，王协康

(1.黄河水利科学研究院，河南 郑州 450003；2.黄河水利委员会，河南 郑州 450000；3.四川大学 水力学与山区河流开发保护国家重点实验室，四川 成都 610065)

近年来，黄河泥沙大幅减少，减幅接近90%[1]，主要原因在于泥沙主要来源区黄土高原的产沙环境发生了显著变化。黄河主要产沙区的林草梯田有效覆盖率已由1978年的不足18%提高至2018年的超过60%[2]。由此可见，研究植被对坡面产沙的作用具有重要意义。

植被减水减沙作用相关研究已有大量成果[3-6]，如：植被冠叶消减雨滴的动能，降低雨滴的溅蚀效果，减小了对坡面土壤的侵蚀作用[7]；植被的截流作用减小产流，从而减少产沙量[8]；地表植被减缓水流速度，减小了土壤的水蚀作用，直接降低水流带走泥沙的能力[9]；植被地下根系增大土壤渗透性及生物固土效果，从而增大土壤抗蚀性[10]；枯落物的覆盖效果增大水流阻力，减小坡面产沙[11]；达到一定覆盖度的植被能够对泥沙起动流速产生一定影响并改变局部床面地形，从而减少输沙量[12]。

植被覆盖度是表示坡面上植被多少的一个重要参数。通过分析植被覆盖度与流域产沙量的关系可知：当植被覆盖度小于30%时，流域产沙量随植被覆盖度的减少而迅速加剧[13]；当植被覆盖度大于70%时，流域产沙量趋于稳定[14]。更进一步地，当林草有效覆盖率小于20%时，植被的减沙作用不太稳定；当林草有效覆盖率达到55%～65%时，植被可以有效遏制流域产沙，且该覆盖率在不同区域的取值不尽相同[15]，类似成果也常有报道[16-17]。可见，植被覆盖度是影响流域产沙量的一个重要指标。为此，学者们专门针对植被覆盖度这一指标进行了坡面水沙运动方面的研究[18-22]。Zhang等[23]给出了植被覆盖度与坡面流阻力系数关系式，植被覆盖度以指数形式存在于公式中。Shang等[24]采用植被密度表示植被覆盖度，认为坡面粗糙高度与植被密度具有线性关系。孙一等[25-26]通过水槽和小区试验，拟合了考虑植被覆盖度的指数形式的坡面流速公式和输沙率公式。Cai等[27]通过水槽试验数据，采用多元非线性回归拟合出植被覆盖下的阻力系数公式，该公式中，植被覆盖度、雷诺数、流量以乘积的形式出现。

然而，现有成果多是由资料分析中得出的偏经验性的结论，所得成果以拟合关系式居多。为此，本文从泥沙运动的角度出发，基于Einstein泥沙运动理论，分析了植被覆盖度对坡面泥沙输移的影响程度，研究成果可从理论层面阐明黄土高原植被覆盖度的有效减沙作用阈值。

1 理论方法

1.1 Einstein输沙率公式

式中： Φ为输沙强度；gb为 输沙率，kg·s-1；g为重力加速度，m·s-2；γs、γ分别为泥沙、水的容重，N·m-3；D为粒径，mm。

根据Einstein泥沙运动理论，床面泥沙运动的全过程可描述为运动—静止—再运动，输沙率实际上取决于泥沙颗粒在床面上停留时间的长短。对于一个特定的沙粒，其进入运动状态或沉积下来的概率在床面各处都相同。任意沙粒在两次连续沉积之间的平均运动距离λD（λ≈100）取决于沙粒的大小、形状，与水流条件无关。因此，无植被时的输沙率公式可写成：

式中：A*为系数，与泥沙形状、运动有关；P为泥沙运动的概率。

1.2 植被作用下的公式建立

图1 局部坡面示意图Fig. 1 Diagram of local slope surface

1.3 相对输沙率特征

图2和3分别为固定泥沙粒径（D=0.05 mm）条件、不同株径Dv下的植被覆盖度-相对输沙率变化曲线及固定株径（Dv=3 mm）条件、不同泥沙粒径D下的植被覆盖度-相对输沙率变化曲线。由图2和3可知：对于某一特定的泥沙粒径，在植被覆盖度相同的条件下，株径越细小，植被的减沙作用越显著；株径越大，植被发挥相同减沙作用所需达到的覆盖度越大。对于某一特定的株径，其对粗泥沙的减沙效果大于细泥沙。

图2 不同株径下的相对输沙率随植被覆盖度的变化曲线Fig. 2 Variation curves of relative sediment transport rate with vegetation coverage under different plant diameters

图3 不同粒径下的相对输沙率随植被覆盖度的变化曲线Fig. 3 Variation curves of relative sediment transport rate with vegetation coverage under different grain sizes

2 径流冲刷试验

2.1 试验设置

采用野外径流冲刷试验对式（6）进行验证。试验地点在黄河水土保持西峰治理监督局南小河沟试验场。选取一个天然坡面建造径流试验区，坡度为15°，如图4所示。试验区由多个宽1 m、长20 m的坡道组成，每个坡道上按一定的覆盖度人工种植苜蓿。覆盖度范围为0～100%，取0、20%、30%、40%、50%、60%、80%、100%共8个等级。图5为部分典型植被覆盖度的照片。坡道上游为进水口，设置6个不同的流量级，范围在0.14～1.40 L·s-1，分别为0.14、0.28、0.56、0.83、1.11、1.40 L·s-1；坡道下游出口处每3 min取样一次含沙量及试验总水量。试验共计进行48组坡面径流冲刷组次。

图4 试验小区示意图及照片Fig. 4 Schematic diagram and photos of the test area

图5 部分典型植被覆盖度的照片Fig. 5 Photographs of some typical vegetation coverage

采用电磁流量计对流量进行测量，待流量稳定后，读取坡道进口处电磁流量计的读数；在坡道出口处收集一定时间内的水量，用Q=Vw/t进行复核（Q为径流流量，L·s-1；Vw为水的体积，L；t为接水时间，s）。采用比重瓶法测定含沙量。试验中，自水流流出坡道时开始，每3 min采集1次浑水样本，每次试验共采集6次浑水样本。

2.2 试验结果

图6为典型流量Q=1.40 L·s-1条件下，不同植被覆盖度的实测含沙量变化过程。由图6可知，坡面流含沙量总体上随着植被覆盖度的增大而显著减小，最后趋于稳定。造成含沙量这种变化趋势的原因是：试验初期，坡面上已经受到侵蚀的土壤最容易被水流带走，导致试验前期的含沙量较大；随着水流的持续冲刷，这部分已侵蚀的土壤逐渐被水流带走，新的土壤出露，但由于粗化现象及植被等多种作用，其相对不易被水流侵蚀，因此含沙量逐渐减小并趋于稳定。此外，当植被覆盖度小于60%时，含沙量随试验时间变化较为明显；当植被覆盖度大于60%时，含沙量很小，且几乎没有明显变化。

图7为野外径流冲刷试验不同流量下各组次平均含沙量随植被覆盖度的变化。由图7可知：随着植被覆盖度的增大，含沙量呈现出明显的减小趋势；当植被覆盖度达到60%时，含沙量就非常小了。

图7 不同流量下平均含沙量随植被覆盖度的变化Fig. 7 Variation of average sediment concentration with vegetation coverage under different flow rates

3 公式适用分析

对于恒定流试验，可将输沙率gb转换成平均含沙量C进行分析，即：后，与试验资料进行对比。

图8为本次径流冲刷试验最大流量Q=1.40 L·s-1、中等流量Q=0.83 L·s-1及最小流量Q=0.14 L·s-1这3个典型流量下的含沙量实测值与计算值的对比。由图8可知：在大流量条件下，计算值与实测值吻合较好；在小流量条件下，计算值与实测值有一定偏差，尤其在植被覆盖度较小的情况，偏差较大。同时，计算值较实测值在整体上也偏大。其原因在于，由于式（6）仅由植被覆盖度Vc、株径Dv、泥沙粒径D决定，与水流条件无关，参考艾里定律，河流中推移质的重量与水流速度的6次方成正比，流量Q的减小及植被的增阻作用势必导致流速降低，进而导致水流强度减弱，实际输沙率也将会进一步减小。

图8 不同流量条件下含沙量计算值与实测值对比Fig. 8 Comparison between calculated and measured values under the conditions of different flow discharge

4 结 论

2）由本文的坡面输沙率公式可得出：在相同的坡面地形及植被覆盖度条件下，株径较小的植被较株径较大的植被更有利于减沙；在粒径较粗的粗泥沙来源区，坡面植被的减沙效果更为明显，同一种植被在粗泥沙来源区的坡面上只需达到一个较小的覆盖度，便可以开始发挥有效的减沙效益；在植被覆盖度较低时，发挥有效减沙效果的覆盖度阈值由坡面泥沙粒径和植被株径决定。

3）本文的坡面输沙率公式的计算值较试验实测值在大流量条件下的吻合较好，在小流量条件下有一定偏差。其原因在于，公式中没有考虑植被覆盖度对水流强度的影响，仍有待于进一步深入研究。