伴随着生产建设活动引起的生态环境问题以及水土保持工作的迫切需求, 工程堆积体边坡侵蚀已成为水土保持研究的热点之一^[1]。人类通过挖方、埋压以及爆破等不同的扰动方式对生产建设项目区的地貌进行了重塑, 产生大量含粒径不一砾石的工程堆积体地貌单元, 其边坡具有坡度大、结构疏松及土壤有机质缺乏等特点^[2]。在暴雨作用下, 工程堆积体水土流失形式常以“点状”或“线型”、单一或综合出现^[3], 造成局部水土流失剧增^[4], 使其成为新增水土流失的重要的物源地之一^[5]。

砾石的存在改变了下垫面的土壤理化性质及地表形态特征, 影响边坡的水文特性和侵蚀过程^[6]。对于传统含砾石坡面, 众多的国内外学者对含砾石土壤的入渗过程^[7-8]、水动力学特性^[9-10]以及产流产沙特征^[11]等方面已开展了大量的试验研究。然而, 对于人为扰动形成的含砾石工程堆积体边坡, 其径流产沙过程及侵蚀动力机制则更为特殊和复杂。史东梅等^[12]比较含砾石紫色土和黄沙壤堆积体边坡径流产沙过程, 表明偏土质(土石比4:1)边坡累计产流量和产沙量均高于土石质混合边坡(土石比3:2), 且产沙率呈连续性多峰多谷变化；李建明等^[13]通过室内模拟降雨试验发现, 坡度在25°左右时砾石对塿土堆积体边坡产沙的促进作用最大；丁亚东等^[14]研究黄绵土散乱锥状堆积体表明了砾石对边坡的土壤侵蚀具有抑制作用；康宏亮等^[15]对北方风沙土堆积体边坡的研究表明, 不同雨强下砾石对边坡侵蚀产沙具有完全相反的作用, 1.0 mm/min雨强时, 砾石促进降雨侵蚀, ＞1.0 mm/min雨强时, 砾石具有显著的减沙效益；以往在工程堆积体水土流失预测模型(经验模型)方面^[16]上取得了一些成果；对于高速公路弃土堆积体^[1]和碾压型堆积体^[17]边坡土壤侵蚀动力机制也建立了控制方程, 均认为径流功率是描述工程堆积体侵蚀动力机制最合理的参数。然而, 这些研究多集中在黄土高原和西南地区, 而对于曾有“红色沙漠”之称的南方红壤区的工程堆积体边坡径流产沙鲜见报导, 且以往针对不同砾石质量分数的土壤侵蚀动力机制的研究较少, 对于不同砾石质量分数条件下的径流水动力学参数特征研究也需进一步加深。

基于众多学者研究的基础上, 本试验采取室内人工模拟降雨的方法, 对不同降雨强度和不同砾石质量分数条件下边坡的侵蚀过程进行观测, 分析不同降雨强度对含砾石红壤边坡径流特性及侵蚀过程的影响, 揭示不同砾石质量分数条件下红壤边坡的径流水动力参数特性及侵蚀动力机制, 以期为红壤区工程堆积体水土流失预测模型提供理论基础。

1 材料与方法 1.1 试验设计

通过对东北黑土区、北方石质山区、西北黄土高原区、黄淮海平原区、贵州土石山区和江西红壤区的368座不同类型工程堆积体的坡长、坡度以及下垫面物质组成等进行实地调查, 经统计分析后表明：2—8 m坡长的堆积体坡面占调查总数的78.4%；坡度多集中在20°—40°；砾石粒径2—15 mm占总样本数的32%, 15—25 mm占总样本数的46%, 粒径＞25 mm占总样本数的22%, 各粒径的质量比大约为3:5:2；砾石质量分数＜40%的工程堆积体占调查样本数的90%以上。在野外调查数据的分析基础上, 结合模拟试验的实际条件以及调查中边坡砾石可被侵蚀搬运的粒径范围对工程堆积体进行室内概化：将本试验砾石筛分为小(2—15 mm)、中(15—25 mm)、大(25—50 mm) 3个粒径等级, 分别过2、15、25、50 mm筛, 然后按质量比3:5:2均匀混合后备用；设计砾石质量分数0(土质边坡)、10%、20%和30%, 砾石与红壤按质量比0:10、1:9、2:8和3:7用搅拌机进行均匀混合, 随后用于填装试验槽；根据李宏伟等^[18]明确定义的我国工程堆积体水蚀测算模型标准小区坡度为25°, 以及已有研究表明坡度影响侵蚀产沙突变的临界坡度大约在25°^[19-20], 再结合野外调查, 最终设计试验坡度为25°, 试验坡长为5 m。

试验土壤采自江西省南昌市新建县撂荒地非耕层土壤。南昌市位于江西省中北部, 属于典型的红壤分布区, 土壤类型为第四纪粘土发育的土壤, 采集的试验土壤为研究区内典型红壤, 具有较好的代表性和适用性。将所取土壤过6 mm筛网剔除土中的树根、砾石等杂物后平铺至室外阴凉处自然风干, 测得土样的土壤颗粒粒径(mm)＜0.002, 0.002—0.02, 0.02—0.05, 0.05—0.1, 0.1—2.0的机械组成分别为17.87%, 32.84%, 21.24%, 9.71%, 18.34%。试验所用砾石采自陕西省商洛市山阳县高速公路旁山体自然滑坡体块石(砾石表面几乎不下渗水分), 据试验粒径要求机械粉碎筛选获取。

试验在中国科学院水利部水土保持研究所模拟降雨大厅进行, 降雨系统为下喷式, 降雨高度18 m, 能满足所有雨滴达到终点速度, 且降雨均匀度能达到80%以上, 可模拟的降雨强度范围在30—350 mm/h之间, 有效降雨面积为27 m×18 m^[21]。试验采用规格为5 m(长)×1 m(宽)×0.5 m(高)的移动式液压可变坡度钢制土槽, 试验装置布设简图如图 1所示。为模拟降雨入渗以及重力水的自由流出, 土槽底层铺厚度为10 cm的石英砂并在土槽底部开有均匀分布的小孔。试验装土深度为40 cm, 装土分为2层, 下层和表层装土深度均为20 cm, 下层进行夯实, 容重控制在1.40—1.55 g/cm³之间, 上层土壤不夯实但对边坡表面进行打磨平整, 对边壁处的边坡尽可能压实以减小边界的影响, 容重控制在1.15—1.30 g/cm³之间, 土层之间进行打毛处理以防止滑塌, 土石混合体经24 h自然沉降, 试验前使用环刀在边坡(上、中、下)分别取样, 控制每场试验土壤质量含水量保持在10%左右。每次试验的装土人员与试验人员固定, 以保证装土工艺一致。

基于对红壤区多年的降雨气象资料的统计分析, 选择暴雨频率较高的30 min最大降雨强度作为本模拟降雨试验降雨强度的设计来源, 取整后设计为1.0、1.5、2.0、2.5 mm/min共4个降雨强度级别, 降雨历时45 min。试验场次共计16场。

1.2 试验步骤

试验开始前, 用遮雨布盖住整个试验土槽, 在土槽的四角分别放置一个雨量筒, 打开降雨器开关进行雨强率定, 率定结果与试验设计雨强误差不超过5%, 雨强率定后快速掀开遮雨布, 直至边坡出现线状水流流出集流槽出口, 开始产流, 并从0开始计时直至45 min整场试验结束。边坡的测流断面选定在距小区顶部1.0—2.5 m和2.5—4.0 m的位置, 产流开始后, 前5 min每1 min在集流槽出口用1000 mL量筒接取径流泥沙样, 5 min后每5 min接取一个径流泥沙样, 并记录接样时间、读取其体积和用电子天平称取质量, 接样同时在每个观测断面分别测定同时刻的流速和流宽。接取的径流泥沙样静置后去除上层清液, 放入105 ℃烘箱内烘干至恒重后称取泥沙质量。

1.3 试验数据处理

文中各参数使用Excel 2013进行数据统计, 使用SPSS 20.0进行数据分析, 使用Origin2018b进行图表绘作。

1) 流速(V)和流宽(L)：流速测定采用KMnO₄染色示踪法, 秒表精度为0.01秒, 测距为1.5 m, 流速乘以修正系数K^[22]得到边坡径流的平均流速；流宽采用自制的尺片测量, 精度为0.01 m。

2) 径流强度(R_r)：单位时间通过单位过水断面的的水流流深, 计算公式为：

(1)

式中, R_r为径流强度, mm/min；Q为径流率, m³/min；b为过水断面宽, m；L为坡长, m。

3) 径流深(h)：测量时段内边坡的平均径流深, 由于水流断面的径流深度较小, 直接人工测量误差较大, 因此采用公式进行计算。计算公式为：

(2)

式中, h为径流深, m；q为测量接样时间T(s)内的径流量, m³。

4) 雷诺数(Re)：判断边坡径流为层流紊流的定量标准, 计算公式为：

(3)

式中, R为水力半径, m；υ为水流黏滞性系数, m²/s。

5) 弗汝德数(Fr)：判别缓急流的定量标准, 表征水流惯性力与重力的比值, 计算公式为：

(4)

式中, g为重力加速度, m/s²。

6) 径流剪切力(τ)：反映径流在流动过程中对表层土壤颗粒分离的能力, 计算公式为：

(5)

式中, τ为径流剪切力, N/m²或Pa；γ_m浑水密度, kg/m³；J为水力坡降, m/m, 用坡度的正弦值近似替换。

7) 单位径流功率(U)：作用于泥沙床面的单位重量水体所消耗的功率, 计算公式为：

(6)

式中, U为单位径流功率, m/s。

8) 径流功率(ω)：表示作用于单位面积水流所需要消耗的功率, 计算公式为：

(7)

式中, ω为径流功率, N m^-1 s^-1。

9) 过水断面单位能(E)：以过水断面最低点作基准面的单位水重的动能及势能之和, 计算公式为：

(8)

式中, E为过水断面单位能, m；a为校正系数, 取值为1。

10) 侵蚀速率(Er)：边坡径流在单位时间单位面积上所搬运的泥沙质量, 计算公式为：

(9)

式中, Er为侵蚀速率, g m^-2 s^-1；M为测量时段T(s)内的产沙量, g。

2 结果与分析 2.1 含砾石红壤工程堆积体边坡径流特征 2.1.1 径流强度变化特征

不同试验条件下边坡径流强度随产流历时的变化过程见图 2。各雨强条件下, 径流强度可描述为先迅速增大后保持稳定的变化趋势, 1.0、1.5、2.0、2.5 mm/min雨强的径流强度分别在产流时间20、10、5、5 min后保持相对稳定, 整体上, 降雨强度越大, 径流强度越快达到稳定状态。在保持稳定期间径流强度的变化差异较小, 其变异系数变化范围在0.01—0.06之间, 属微弱变异程度。

对稳定径流强度值进行了计算和分析, 不同砾石质量分数边坡的次降雨稳定径流强度值如图 3所示。稳定径流强度随雨强的增大而增大, 经分析二者呈极显著线性关系(P＜0.01)。1.0 mm/min雨强时, 10%砾石质量分数稳定径流强度较土质边坡稳定径流强度增大8.4%, 砾石质量分数＞10%时, 稳定径流强度随砾石质量分数的增大而减小, 减小幅度在4.5%—11.8%, 且砾石质量分数为10%时取极大值；雨强＞1.0 mm/min时, 10%砾石质量分数稳定径流强度较土质边坡稳定径流强度减小7.3%—14.2%, 砾石质量分数＞10%时, 稳定径流强度随砾石质量分数的增大而增大, 增大幅度在1.3%—17.0%, 且砾石质量分数为10%时取极小值, 说明在砾石质量分数为10%时, 砾石对稳定径流强度的影响最大且存在一个阈值(10%左右)。

2.1.2 径流水动力学参数变化特征

不同试验条件下边坡径流水动力学参数见表 1。

1) 流速。各砾石质量分数条件下, 流速随雨强的增大而增大。1.0 mm/min雨强条件下, 土质边坡流速最小(0.049 m/s), 含砾石边坡较土质边坡的流速增大23.3%—68.6%；2.0、2.5 mm/min雨强条件下, 土质边坡流速最大(0.108、0.127 m/s), 含砾石边坡较土质边坡的流速分别减小1.3%—12.1%、0.5%—20.9%；而在1.5 mm/min雨强条件下, 砾石质量分数＜10%时, 流速随砾石含量的增大而减小, 砾石质量分数＞10%时, 流速随砾石含量的增大而增大。由此可知, 不同雨强和砾石质量分数均会对流速产生了重要影响, 雨强＜1.5 mm/min时, 砾石的存在促进边坡径流的流动, 流速随砾石含量的增大而增大, 雨强＞1.5 mm/min时, 砾石的存在抑制边坡径流的流动, 流速随砾石含量的增大而减小。

2) 流型与流态。试验条件下, 雷诺数随降雨强度的增大而增大, 边坡径流以层流为流动形态。弗汝德数随随雨强的增大而增大。1.0 mm/min雨强条件下, 弗汝德数变化范围为0.526—0.988, 边坡径流为缓流, 此时, 砾石的存在增大了径流的弗汝德数, 相比土质边坡增大了24.5%—87.8%；2.0、2.5 mm/min雨强下, 各砾石质量分数边坡径流以急流为主, 土质边坡的弗汝德数达最大, 边坡的砾石降低了径流的弗汝德数, 降低幅度为4.2%—13.0%。

3) 径流功率。径流功率随雨强的增大而增大, 经分析二者呈极显著幂函数关系(P＜0.01)。1.0、1.5 mm/min雨强条件下, 各边坡径流功率由大到小对应的砾石质量分数依次为10%、20%、0、30%和30%、20%、0、10%；2.0、2.5 mm/min雨强条件下, 土质边坡径流功率为最大分别为0.539、0.723 N m^-1 s^-1, 含砾石边坡的径流功率较土质边坡的降低幅度分别为6.6%—13.6%和10.2%—17.9%。

2.2 含砾石红壤工程堆积体边坡侵蚀特征 2.2.1 侵蚀速率随产流历时变化过程

不同试验条件下边坡侵蚀速率随产流历时的变化过程见图 4。1)1.0、1.5 mm/min雨强下, 侵蚀速率先上升(1—3 min)后下降(4—35 min), 最后趋于稳定(40—45 min), 各边坡侵蚀速率分别在3—5、1—3 min达到最大, 且含砾石边坡最大侵蚀速率是土质边坡的1.12—2.59倍。在下降和稳定阶段, 1.0 mm/min雨强时, 土质边坡相对于含砾石边坡的侵蚀速率低, 土质边坡侵蚀速率稳定在0.19 g m^-2 s^-1左右, 含砾石边坡稳定在0.22 g m^-2 s^-1左右；而在1.5 mm/min雨强时, 土质边坡侵蚀速率整体上较含砾石边坡高, 土质边坡侵蚀速率的稳定值在1.06 g m^-2 s^-1左右, 含砾石边坡侵蚀速率稳定在0.59 g m^-2 s^-1左右。2)2.0、2.5 mm/min雨强条件下, 侵蚀速率先迅速下降(1—5 min)后波动变化(10—45 min)。在波动阶段, 2.0 mm/min雨强时, 土质和10%砾石含量边坡侵蚀速率波动增大, 而20%、30%砾石含量边坡持续下降, 且各边坡侵蚀速率由大到小对应的砾石质量分数为0、10%和20%(30%)；2.5 mm/min雨强时, 30%砾石含量边坡侵蚀速率逐渐降低, 土质、10%和20%砾石含量边坡侵蚀速率波动剧烈, 变异系数分别为0.25—0.27, 但其变化趋势并不明显。试验条件下, 雨强越大, 砾石质量分数越低的边坡侵蚀速率变化多呈多峰多谷。

2.2.2 次降雨侵蚀速率变化特征

次降雨侵蚀速率随砾石质量分数和雨强的变化见图 5。1.0 mm/min雨强条件下, 各边坡的侵蚀速率由小到大按砾石质量分数表现为0＜30%＜20%＜10%, 砾石质量分数为10%、20%和30%边坡的侵蚀速率分别是土质边坡的1.51、1.50、1.29倍；1.5—2.5 mm/min雨强条件下, 土质边坡的侵蚀速率最大, 在0.95—4.32 g m^-2 s^-1之间, 3种雨强条件下10%、20%、30%砾石质量分数边坡的侵蚀速率较土质边坡的分别降低21.2%、13.1%、5.0%和13.1%、44.8%、44.9%和21.1%、40.2%、64.4%。由此可知, 1.0 mm/min雨强时, 砾石的存在加剧土壤的侵蚀, 1.5—2.5 mm/min雨强时, 砾石的存在降低土壤的侵蚀速率。

在砾石质量分数一定的条件下, 侵蚀速率随雨强的增大而增大, 二者呈显著幂函数关系(P＜0.05), 但各砾石质量分数边坡的侵蚀速率随雨强增大的增长幅度存在差异, 相比于1.0 mm/min雨强条件下, 其中土质边坡侵蚀速率可增大4.35—19.78倍, 砾石质量分数10%、20%和30%的边坡侵蚀速率分别可增大2.28—10.37倍、2.53—7.89倍和3.22—5.47倍, 相比于土质边坡, 整体上含砾石边坡的侵蚀速率增大倍数较小。通过对侵蚀速率与砾石质量分数和雨强进行回归分析, 结果表明侵蚀速率与砾石质量分数和雨强呈极显著二元线性函数关系(式10, P＜0.01, R²=0.84, N=16)：

(10)

式中, I为降雨强度, mm/min；G为砾石质量分数, %。

2.3 含砾石红壤工程堆积体边坡侵蚀动力机制

表 2为次降雨侵蚀速率与各水动力学参数之间的相关性分析。结果显示边坡侵蚀速率(Er)与各水动力学参数之间均有良好的相关关系, 相关系数具体表现为：r(ω)＞r(Re)＞r(E)＞r(V)＞r(Fr)＞r(τ)。

图 6描述了次降雨土壤侵蚀速率与侵蚀动力参数的关系。20%和30%砾石质量分数边坡的侵蚀速率(Er)与剪切力(τ)可用幂函数表示；0和10%砾石质量分数边坡的侵蚀速率(Er)与单位径流功率(U)存在线性关系；各砾石质量分数边坡侵蚀速率(Er)与径流功率(ω)、过水断面单位能(E)均呈显著性线性关系(P＜0.05)。由各土壤侵蚀动力机制控制方程分析可知, 不同砾石质量分数条件下, 侵蚀速率与径流水动力学参数关系式中的横坐标截距(发生侵蚀临界值)与斜率(可蚀性参数)具有差异性同时也具有规律性, Er与ω、E之间的关系中发生侵蚀的临界值及可蚀性参数表现出一致性, 其大小按砾石质量分数排列均为10%＞0＞20%＞30%, 砾石质量分数为0、10%、20%、30%的边坡发生侵蚀的临界径流功率分别为0.213、0.238、0.212、0.084 N m^-1 s^-1, 对应的可蚀性参数分别为8.13、8.60、6.26、2.90 s²/m²；临界过水断面单位能分别为1.0×10^-3、1.1×10^-3、0.9×10^-3、0.6×10^-3 m, 对应的可蚀性参数分别为3726、4481、3017、1224 g m^-3 s^-1。10%的砾石质量分数边坡发生侵蚀的临界功率及临界过水断面单位能分别为0、20%、30%的砾石质量分数边坡的1.14、1.14、3倍和1.10、1.22、1.83倍, 可蚀性参数分别增大5.8%—196.6%和20.3%—266.0%。

3 讨论 3.1 砾石对边坡径流特性的影响

工程堆积体作为一种特殊的地貌单元, 其径流动力特征具有特殊性和复杂性。本研究表明, 砾石质量分数对稳定径流强度的影响存在一个10%左右的阈值, 雨强为1.0 mm/min时, 10%的砾石质量分数稳定径流强度最大, 当砾石质量分数增至10%时, 表层砾石的增加使径流弯曲度对于稳定径流强度的增大作用占据了优势, 延长了水分下渗路径, 导致稳定径流强度增大, 而在砾石质量分数大于10%时, 砾石的存在显著增多了土壤大孔隙通道, 促进水分入渗, 导致稳定径流强度减小^[23]；而雨强大于1.0 mm/min时, 砾石质量分数为10%时径流强度最小, 这与边坡形成的细沟有关, 当砾石质量分数增至10%时, 边坡表层土壤逐渐被剥蚀, 凸出表层土壤的砾石增多, 坡面水流沿着凸出表层土壤的砾石下渗, 导致稳定径流强度减小, 当砾石质量分数继续增大至20%、30%时, 此时, 细沟中的砾石较多, 在砾石的汇流和入渗面相对减少的共同作用下, 使稳定径流强度增大^[24]。土壤入渗特性与径流强度的关系密不可分。朱元骏等^[25]通过模拟降雨研究表明, 砾石质量分数为10%时土壤入渗率处于临界状态, 当砾石质量分数低于10%时, 入渗率与砾石质量分数呈正相关, 当碎石含量超过10%时, 入渗率反而降低；时忠杰等^[8]研究指出, 在0—40 cm土层, 砾石体积分数与稳渗速率的临界值大致处于15%—20%, 这些研究均与本研究结果具有相似性, 但得出的临界值大小具有一定的差异, 这可能与土壤类型有关。边坡径流的流动速度受降雨和砾石的影响。1.0 mm/min雨强时, 降雨动能较小, 边坡上只有小部分的砾石完全凸出, 而平滑的小砾石镶嵌在表层土壤促进了径流的流动^[26], 这与在1.0 mm/min雨强条件下对含砾石紫色土^[24]、塿土^[27]的研究结果一致；雨强为1.5 mm/min时, 砾石对径流流速的影响并不明显, 原因是砾石对径流的阻碍和促进作用相互抵消；雨强＞1.5 mm/min时, 降雨动能较大, 众多的砾石很快地从表层土壤凸出, 增大了表层土壤的粗糙度, 消耗了径流的动能, 从而延缓径流的流动^[28], 这与Ban等^[10]在放水条件下对冻土的研究结果类似。综上, 砾石对边坡径流特性的影响主要有以下2个方面：1)土壤表层及以下的砾石通过改变表层土壤粗糙度和水分的入渗来影响边坡径流^[29]；2)土壤表层的砾石或部分嵌入表层土壤的砾石影响降水拦截、岩石流产生以及表层径流的分散与汇合, 从而决定边坡产流情况^[30]。

3.2 砾石对边坡产沙的影响

不同雨强及砾石质量分数条件下边坡的侵蚀过程具有差异性。试验条件下, 1.0、1.5 mm/min雨强时, 边坡前期侵蚀方式主要为层状面蚀, 而后期主要受表层结皮的影响, 从而导致侵蚀速率随产流历时表现为先迅速暂上升而后缓慢下降至趋于稳定的变化趋势；雨强增大后(2.0 mm/min、2.5 mm/min), 随降雨侵蚀的进行, 在产流10 min左右后, 边坡径流侵蚀方式由面蚀逐渐向断续细沟转变, 最后伴随着跌坑和侵蚀穴的出现而形成连续细沟^[31], 连续细沟的出现直接导致了土壤侵蚀速率成倍速度的增加^[32], 而细沟的发育过程中同时会伴随着重力侵蚀的发生, 细沟中边壁土体的崩塌是导致侵蚀速率迅速增大和产生波动的重要原因^[33]。但雨强为2.0 mm/min时, 边坡主要形成短、窄、浅型细沟(长＜60 cm, 宽＜7 cm, 深＜3 cm), 细沟较浅, 发生重力崩塌的几率较低, 侵蚀速率发生波动次数少；雨强增大至2.5 mm/min时, 径流拥有更大的动能和冲刷力, 边坡形成长、宽、深型细沟(长100—150 cm, 宽10—15 cm, 深3—5 cm), 此时细沟较深, 易发生重力崩塌, 侵蚀速率发生多次波动, 但砾石质量分数为30%时, 细沟内砾石较多, 可被输移的土壤颗粒少, 从而造成侵蚀速率保持相对稳定, 这也是较其他砾石质量分数条件下侵蚀速率小的重要原因。本研究的侵蚀速率变化过程与塿土堆积体边坡^[34]和黄绵土堆积体边坡^[14]侵蚀速率变化过程具有明显差异, 由此可见, 不同的雨强、土壤类型、边坡形态均会影响含砾石边坡的侵蚀过程。对于含砾石工程堆积体边坡的侵蚀产沙而言, 砾石的存在既可以增加土壤侵蚀量^[35], 也能与土壤侵蚀量呈负相关关系^[14], 在不同雨强条件下, 砾石对侵蚀产沙的影响也有所不同^[26]。本研究表明, 1.0 mm/min雨强时, 砾石的存在加快了径流的流动, 边坡径流具有更大的动能, 且边坡上凸出的砾石较少、砾石覆盖度较低, 此时砾石对边坡的保护作用较弱并且径流对土壤产生“涡旋侵蚀”^[36], 进而加剧土壤侵蚀, 与含砾石风沙土边坡研究结果一致^[26], 而与含砾石紫色土边坡的研究结果存在差异^[24], 试验结果的差异可能与坡度有关；1.5—2.5 mm/min雨强时, 砾石对边坡侵蚀的抑制作用主要有以下两个方面：1)大量砾石从表层土壤凸出, 增大地表粗糙度, 削弱了径流冲刷能力；2)径流的集中冲刷, 边坡侵蚀方式以细沟侵蚀为主, 造成砾石的大量裸露, 使边坡可被侵蚀输移的土壤量减少, 从而导致侵蚀量与砾石的关系转换为负相关关系^[15]。细沟侵蚀在工程堆积体边坡的侵蚀中产生了重要作用, 随着表层土壤逐渐被剥离, 细沟不断发育, 砾石在边坡上的覆盖率也随之增大, 在一定程度上, 砾石质量分数与砾石覆盖度对边坡土壤侵蚀的作用具有相似性^[37]。

3.3 砾石对边坡侵蚀动力机制的影响

通过研究含砾石红壤边坡的侵蚀速率与各水动力学参数的关系发现：“能量”参数径流功率(ω)、过水断面单位能(E)均可用于描述含砾石红壤边坡侵蚀动力机制, 但在最优参数选择方面(R²结合相关性分析), 径流功率(ω)是描述红壤区含砾石工程堆积体侵蚀动力机制更为合理的参数, 这与诸多学者在放水条件或模拟降雨条件下对工程堆积体的研究一致^{[1, 26, 38]}。根据土壤侵蚀速率和径流功率的侵蚀动力机制控制方程得到的红壤区工程堆积体可蚀性参数介于2.90—8.13 s²/m²之间, 是牛耀彬等^[38]对有工程措施堆积体得到的可蚀性参数的1.0—5.8倍, 这说明了工程措施对防治堆积体水土流失的重要性, 整体上比张乐涛等^[1]对陡坡堆积体得到的可蚀性参数小(8.0 s²/m²), 这与土壤在较陡边坡上对侵蚀更为敏感有关。同时, 在不同条件下运用不同的控制方程得出的土壤可蚀性参数大小之间并无必然联系, 在运用通用土壤流失方程等预测模型时, 对土壤可蚀性参数的选择需要根据实际具体条件而定。本研究中得到的侵蚀临界值和可蚀性参数的大小按砾石质量分数排列均为10%＞0＞20%＞30%, 这说明了10%的砾石质量分数边坡在降雨下难以启动侵蚀, 但是又对侵蚀的敏感性高, 这可能与入渗特征有关, 整体上10%砾石质量分数条件下边坡入渗最大, 导致坡面径流强度低, 坡面径流要达到启动侵蚀的能力不易, 但是一旦坡面侵蚀被启动, 最大的入渗作用使边坡含水量增大, 导致抵抗侵蚀的能力减弱^[39]。

4 结论

(1) 砾石质量分数对稳定径流强度的影响存在阈值, 稳定径流强度在雨强＞1.0 mm/min时随砾石质量分数的增加先减小后增大, 在1.0 mm/min雨强时则表现为先增大后减小, 且均在10%砾石质量分数下达到极值。流速随雨强的增大而增大, 边坡径流以层流为流动形态；1.0 mm/min雨强条件下, 边坡径流为缓流, 砾石促进了径流的流动, 使弗汝德数增大24.5%—87.8%；2.0、2.5 mm/min雨强条件下, 边坡径流为急流, 流速随砾石含量的增大而减小, 弗汝德数降低4.2%—13.0%。(2)1.0—1.5 mm/min雨强时, 侵蚀速率在产流前期达到最大, 随雨强增大到2.5 mm/min且砾石质量分数越低时, 边坡越易发生细沟侵蚀并伴随重力崩塌现象, 侵蚀速率呈多峰多谷变化趋势。1.0 mm/min雨强时, 砾石存在加剧了土壤侵蚀, 产沙增幅达28.7%—50.5%；雨强＞1.0 mm/min时, 砾石减沙效益为5.0%—64.4%；侵蚀速率与雨强、砾石质量分数三者之间可用二元线性函数表示。(3)土壤侵蚀速率与径流功率及过水断面单位能之间均存在线性函数关系, 而径流功率是描述红壤区含砾石工程堆积体侵蚀动力机制最为合理的参数, 其可蚀性参数及发生侵蚀临界径流功率从大到小对应的砾石质量分数均为10%、0、20%和30%。