不同质地重塑土坡面水沙定量关系研究*

倪世民 张德谦 冯舒悦 王军光 蔡崇法

（华中农业大学水土保持研究中心，农业农村部长江中下游耕地保育重点实验室，武汉 430070）

水力侵蚀是引起土地退化和水土流失的主要驱动因子，也是水体面源污染的主要来源之一[1]。水力侵蚀包括由降雨溅蚀和地表径流引起的土壤分离、泥沙输移、泥沙沉积3 个过程[2-3]，研究各过程间相互转化、相互影响的机理是建立土壤侵蚀物理模型的前提条件[4]。土壤分离是指径流将土壤颗粒从土壤基质表面剥离的过程[5]，土壤分离速率（soil detachment rate）被定义为单位时间单位面积土壤表面被径流剥蚀掉的土壤质量[6]。

在已有的研究中，一些学者通过定床水槽试验的方法，对室内直径或宽度为10 cm 的小尺寸样品进行径流冲刷来获取土壤分离速率[7-8]。在这种方法中，土壤分离被假设为一个单独的土壤侵蚀系统，考虑在限定性细沟内集中水流对土壤表面的侵蚀作用，并没有涉及土壤侵蚀过程中细沟形态、水力学特性及床面粗糙度等因素的动态变化。在自然条件下，坡面侵蚀是从层状面蚀开始逐步向细沟状面蚀转变，侵蚀形成的细沟是一个粗糙、复杂、不规则的形态，细沟由下切侵蚀、侧向坍塌等引起的形态变化对细沟中的泥沙输移与沉积过程、水力学特性具有显著的影响[9]。研究表明，在土壤侵蚀过程中径流对土壤的剥蚀能力随着径流含沙量的增加呈线性降低，当径流泥沙载荷等于输沙能力时，土壤分离速率达到最小值[10]。此外一部分学者采用动床试验的方法获取土壤分离速率，考虑了土壤侵蚀过程中各个过程的相互作用与动态变化：坡面的径流水力学特性、侵蚀产沙及侵蚀形态三者相互影响、相互制约，是一个包含土壤分离、泥沙输移及泥沙沉积三个过程的动态侵蚀系统[11-12]。在GUEST 模型中，土壤分离与泥沙沉积被认为是一个同时发生的连续过程，坡面的侵蚀产沙是土壤分离、泥沙沉积、泥沙输移相互作用的结果[13]。土壤分离能力是细沟侵蚀中下切侵蚀、溯源侵蚀、侧蚀等各个部位分量的总和。对于一个完整的坡面侵蚀系统而言，动床试验考虑了土壤分离、土壤沉积、土壤输移3 个过程，而不是注重水力侵蚀过程中单独的土壤分离过程。可见，定床试验并不能充分地反映自然土壤坡面上的土壤分离过程，相关研究结果的应用受到限制与质疑[4，11-12]；而动床试验更贴合自然条件下土壤侵蚀过程，两者之间的主要差异体现在水流阻力、水流形态、侵蚀泥沙来源等方面。动床试验条件下，土壤分离速率的概念被考虑为整个侵蚀面的土壤侵蚀速率（soil erosion rate），定义为单位时间单位侵蚀面积坡面土壤在水流侵蚀动力的作用下被剥蚀的土壤颗粒质量[10]。土壤分离速率与土壤侵蚀速率两者的概念相同，但是土壤侵蚀速率（或土壤剥蚀率）更适合描述动床条件下的坡面土壤侵蚀过程。

回顾以往研究，诸多学者使用了不同类型的土壤进行关于土壤分离的研究，如黄土[9，11-12]、培肥土壤[8]、红壤[7]、工程堆积体[15]、农田[16]等。这些研究中，土壤分离速率通常通过水力学参数进行表达，如在WEPP 模型中使用水流剪切力对其进行描述，在EUROSEM模型中使用单位水流功率进行描述，在GUEST 模型中使用水流功率进行描述等[13，17]，但由于试验材料与条件的差异，究竟哪个是最适合的参数尚无定论。同时，已有的研究大多集中在单一质地或者粒径范围较窄的土壤上，而关于土壤侵蚀速率随土壤质地变化规律的相关研究较为有限；且由于试验土壤的差异，不同数据集之间所得的结果仅仅适用于特定的土壤，建立的模型在不同的土壤中应用时需要进行校准[18]。因此，在进行不同质地土壤间土壤侵蚀速率的系统试验研究具有重要的理论意义和现实价值。鉴于上述研究背景与存在问题，本研究采用了红壤（黏土）和砂土混合配制的不同质地的土壤为研究对象，通过动床试验的方法，探讨了不同质地土壤坡面的土壤侵蚀速率，分析了土壤侵蚀速率与坡面的细沟形态特征之间的响应关系，确定了水力学参数与土壤侵蚀速率间的定量关系，以期为坡面土壤侵蚀预报模型研究提供理论基础。

1 材料与方法

1.1 供试材料

试验材料为人工混合配制的重塑土。配制土样所用的黏土采自鄂南丘陵区湖北省咸宁市咸安区贺胜桥镇（114°42′12″E，29°99′35″N），采样深度为0～15 cm，土壤类型为第四纪红黏土发育的红壤，土地利用方式为林地；试验所用的沙土为普通工程沙土。土样被自然风干并过5 mm 筛，以剔除石块、根系等杂质。不同质量分数的黏土与沙土在干燥条件下按设置的含沙量梯度被均匀混合，配置不同质地类型的重塑土。根据已有研究，试验设置了5 个含沙量梯度：沙土质量百分数依次为0、30%、50%、70%、100%[19]。本研究中，5 种质地的重塑土（含沙量由低至高）分别被定义为试验组S1、S2、S3、S4、S5，根据美国土壤质地分类制，土壤质地依次为粉质黏土、黏土、砂质黏壤土、砂质壤土、砂土，涵盖了从黏土至砂土的质地范围。土样的抗剪强度（τ0）通过ZJ 型应变式控制直接剪切仪测定，环刀含水率控制在30%，垂直法向压力分别为 50、100、150 和 200 kPa，剪切速率为0.8 mm·min-1[19]。根据库仑-摩尔公式：τ0=C+λtanφ，计算土壤黏聚力C（kPa）和内摩擦角φ（°），其中λ为法向压力[19]。土样的基本理化性质采用常规方法测定（表1）[19]。

1.2 研究方法

试验通过室内模拟径流冲刷动床试验的方法进行。试验装置主要由试验土槽、供水装置和接样装置组成（图1）。土槽为长3.00 m、宽1.00 m、高0.35 m的变坡式大型钢制土槽，分为相同的两个平行土槽，且可在0～30°的范围调节坡度。土槽后端安装了内部有不锈钢消能板的稳流箱，以保证水流平稳流入土槽。供水设备由蓄水池、恒压水泵、阀门组和水管组成，在试验过程中持续向蓄水池注水，放水流量由阀门控制。

土槽的填土容重控制在1.35 g·cm-3，接近田间的自然状况。填土采取分层填土的方式，边填边压实，且在填装上层土壤前，将下层土壤的表层进行抓毛处理，以保证填土尽量均匀。同时，坡顶15 cm处被覆盖纱布以削弱稳流箱出口跌水能量与边缘效应的影响。在每次试验开始前，使用降雨器以30 mm· h-1的雨强对土槽进行预湿润处理，放置后当土壤含水率降至30%时开始进行试验。试验设置了2 个坡度（5°与15°）和4 个流量（2、4、6、8 L·min-1），每组试验设置两个平行，试验时间设置为15 min。在前期预试验中，由于在15° 坡度、8 L·min-1流量条件下坡面侵蚀强烈，造成短时间内坡面下段发生大量沉积，因此未进行该条件下的试验。坡面产流后，沿坡面每隔20 cm（共15 个断面）使用直尺获取径流水宽；径流表面流速采用高锰酸钾染色法分段测定，水温通过稳流槽内的温度计读取。试验过程中在出水口收集泥沙和径流样品，以获取含沙量、产沙量和径流量，接样容器为接样瓶（间隔30 s 接1 次）和接样桶（间隔1 min 接1 次）。在试验结束后，通过直尺法沿坡面获取细沟（15 个断面）的形态参数（沟深、沟宽、宽深比）。

表1 试验土样的基本理化性质 Table1 Basic physicochemical properties of the soil samples used in the experiment

图1 试验装置 Fig.1 Sketch of the experiment apparatus

1.3 参数获取方法

综合有关研究，本研究选取下列相关参数。土壤侵蚀速率（土壤剥蚀率）为单位时间单位侵蚀面积坡面土壤在水流侵蚀动力作用下被剥蚀的土壤颗粒质量[11]，表征径流对坡面土壤的分离能力，其求解方法如下：

式中，Dr为土壤侵蚀速率，kg·m-2·min-1；mt为累计产沙量，kg；B为水宽，cm；L为沟长，m；T为时间，min。

平均流速（V）通过试验过程中染色法测得的表面流速根据相应的流态乘以换算系数得到：

式中，V为平均流速，m·s-1；Vs为表面流速，m·s-1；k为相关的换算系数，根据对应的流态（层流、过渡流、紊流）分别取值0.67、0.70、0.80[20]。

坡面径流的流态通过水流内部的紊乱指标来描述，通常使用的判别参数为雷诺数与弗汝德数[21]：

式中，R为水力学半径，m，本试验条件下可近似用水深h代替；v为黏滞系数，cm2·s-1；g为重力加速度，g=9.8 m·s-2；h为径流水深，cm。

水力学参数选用水流剪切力、水流功率、单位水流功率，分别通过下式计算[22]：

式中，τ为水流剪切力，Pa；ρ为水流容重，kg·m-3；J为水力坡度，m；ω为水流功率，N·m-1s-1；Pr为单位水流功率，m·s-1。

细沟断面宽深比是无量纲参数，表示细沟的断面形态在水平方向和竖直方向上尺寸的相对大小，其数学表达式为：

式中，μ为细沟断面宽深比；dr为细沟沟宽，cm；hr为细沟沟深，cm。宽深比μ的值越大，断面形状越趋近“宽浅式”，宽深比μ的值越小，断面形状越趋近“窄深式”。

本文采用方程决定系数和均方根误差来评估拟合函数的表现，筛选试验结果，其求解形式为：

式中，R2为方程决定系数，RMSE 为均方根误差，Mi为实测值，为实测值的平均值，Pi为预测值，为预测值的平均值，n为样本数。

1.4 数据处理方法

数据处理与分析在Excel 2016、Origin 2017 与SPSS 19.0 软件中完成，相关分析采用Spearman 双侧显著性检验。

2 结果与讨论

2.1 坡面侵蚀产沙特征

各试验条件下重塑土坡面的土壤侵蚀速率如图2所示。土壤侵蚀速率随着坡度、流量与土壤质地的变化呈现明显的规律性：土壤侵蚀速率随着流量的增加而增加，随着坡度的增加呈现剧烈增加的趋势，表明坡度对土壤侵蚀速率的影响占主要地位。这与张光辉等[22]的研究结果不同，原因与在动床条件下坡度的增加使坡面土体的稳定性大幅降低有关。对于不同质地的土壤而言，随着土壤质地的偏砂化土壤侵蚀速率呈现“峰状”分布，即先增加后减小的趋势，并在S3（50%含沙量）坡面处取极值。坡度和流量是坡面水流侵蚀力变化的源动力，大量研究表明，坡面产沙与坡度、流量呈正相关性，在一定范围内，坡度和流量越大产沙量越大[8，12，15]。

土壤性质（尤其是土壤质地）对细沟形态特征、土壤侵蚀速率具有重要的影响[23-24]。相比黏重的土壤，砂质的土壤由于土壤黏结性差而易被径流所剥离。同时，土壤侵蚀速率取决于径流输移能力和径流剥蚀速率中的较小者[25]。在相同的径流侵蚀力条件下，土壤侵蚀速率呈现“峰状”分布，表明质地黏重的土壤的产沙过程主要受径流剥蚀能力的限制，质地偏砂的土壤的产沙过程主要受径流搬运能力的限制。与S4（70%含沙量）、S5（100%含沙量）坡面相比，虽然S1（含沙量为0）坡面产沙略高，但显而易见，S1 坡面的产沙过程受径流剥蚀能力的限制，这表明红壤较高的黏结性与较好的团粒结构可以在一定程度上增加坡面的抗蚀性[7]；而S4、S5坡面则表现为径流剥离的土壤大于径流的输沙能力，从而在坡面下段发生较为明显的沉积。

2.2 坡面断面形态对侵蚀产沙的响应

图2 不同试验条件下土壤侵蚀速率 Fig.2 Soil erosion rate relative to treatment

坡面的侵蚀形态受控于坡面的径流侵蚀力与土 壤抗蚀性的共同作用[26]，明确坡面侵蚀形态与坡面侵蚀产沙的关系具有重要意义。图3所示为各试验条件下的细沟断面宽深比。细沟断面宽深比随坡度的增大呈现减小的趋势，断面形态由“宽浅式”趋向于“窄深式”。在缓坡条件下，流量对细沟断面形态没有显著的影响，在陡坡条件下，随着流量的增加细沟断面宽深比随之减少。土壤质地的差异使不同坡面的细沟断面形态呈现明显的规律性，质地越砂化（坡面含沙量越高）断面形态越趋近于“宽浅式”（图4）。土壤的颗粒组成是影响土壤抗蚀性的重要因素，颗粒组成越细的土壤，黏结力越强，土壤抗蚀性和形成的细沟的稳定性越强[27]。细沟侵蚀为径流输沙提供了输移通道，一方面黏结力较差的土壤坡面细沟沟壁稳定性差，沟壁坍塌引起的侧向侵蚀频繁，使细沟断面形态趋向“宽浅式”；另一方面宽浅的细沟引起径流输沙能力减小，坡面上方侵蚀的泥沙因径流搬运能力的限制而在下方发生沉积，致使断面形态进一步趋向“宽浅式”。

为进一步揭示坡面产沙特征与细沟形态特征之间的关系，对土壤侵蚀速率、累计产沙量与细沟形态参数进行了相关分析，以判断其内在的联系。如表2所示，细沟的平均沟深和土壤侵蚀速率、坡面累计产沙量具有极显著的正相关关系（P＜0.01）。相关研究表明，径流对坡面土壤的侵蚀使坡面沟床形态发生改变，沟床形态的改变又反过来影响径流的水力学特性与坡面的侵蚀产沙[28]。径流的下切侵蚀是坡面产沙的主要来源，在本研究中平均沟深能显著地表达坡面的产沙情况。此外，细沟断面宽深比和土壤侵蚀速率、坡面产沙量具有负显著相关关系（P＜0.01），细沟的断面形态越加“窄深”，土壤侵 蚀速率和坡面产沙量越大。坡面的侵蚀过程是一个径流能量耗散的过程[26]，细沟形态的改变、泥沙的输移及坡面产沙伴随着能量的耗散，在相同的流量下，狭窄细沟内形成的集中水流具有较大的径流能量，从而增强了径流的土壤剥蚀能力和搬运能力[24]。

图3 不同试验条件下细沟宽深比 Fig.3 Rill width-depth ratio relative to treatment

图4 不同质地土壤坡面细沟断面形态 Fig.4 Rill width-depth morphology of slope relative to different soil textures

表2 坡面产沙特征与细沟形态参数的相关关系矩阵（n=35） Table2 Correlation matrix between characteristics of slope sediment yield and rill morphological parameters

2.3 水力学参数与土壤侵蚀速率的关系

土壤侵蚀速率是水力学参数与土壤参数的函数，深刻理解坡面侵蚀产沙对水动力学参数的响应关系是建立土壤侵蚀物理模型的基础。表3所示为土壤侵蚀相关参数与和水力学参数的相关关系。土壤侵蚀速率与坡度、水流剪切力、水流功率及单位水流功率呈极显著的相关关系（P＜0.01），与单宽流量、流速、雷诺数、弗汝德数呈显著的相关关系（P＜0.01），相关系数r从大到小依次为：Pr＞ω＞S＞τ＞V＞Fr＞q，表明在相同试验条件下单位水流功率对土壤侵蚀速率影响最大。坡度对坡面水力侵蚀具有重要影响，与土壤侵蚀速率关系密切（r=0.812，P＜0.01），与前人的研究结果一致[11]。在本试验中，坡面侵蚀属于细沟侵蚀的范畴，坡度增加时，坡面土体的稳定性下降，细沟沟壁坍塌频繁，导致产沙量剧增。雷诺数与弗汝德数是表征径流紊乱程度的参数，两者可在一定程度上反映坡面的土壤侵蚀速率和坡面侵蚀形态特征[24]。“窄深式”细沟内径流湍急、紊乱，流态趋向于紊流和急流的范畴，径流的土壤剥蚀能力和搬运能力较强；而“宽浅式”细沟内径流流态趋向于层流和缓流的范畴，土壤剥蚀能力和搬运能力较低。水流剪切力、水流功率、单位水流功率是坡度与流量的函数，是与土壤侵蚀速率密切相关的水力学参数，不同学者选用描述土壤侵蚀速率的参数存在分歧。在本研究中，单位水流功率是描述坡面土壤侵蚀速率的最佳水动力学参数（r=0.911，P＜0.01）。

2.4 建立水力学参数与土壤侵蚀速率的定量关系

在过去的几十年间，诸多学者对水力学参数与土壤侵蚀速率之间的关系进行了定量表述[8，14，29]。土壤侵蚀速率首先是简单水力变量的函数，研究发现基于坡度和单宽流量的幂函数可以很好地描述土壤侵蚀速率：

表3 土壤侵蚀相关参数与水力学参数的相关关系 Table3 Correlation between the relevant parameters of soil erosion and hydraulic parameters

式（11）中，坡度因子指数（2.59）大于单宽流量（0.54），表明坡度对土壤侵蚀速率的贡献大于流量；决定系数（R2）为0.644，方程的拟合程度良好，坡度和单宽流量的函数可以在一定程度上较好地描述土壤侵蚀速率。根据表3中水力学参数与土壤侵蚀速率的相关性，对单位水流功率（Pr）与土壤侵蚀速率的定量关系进行分析，发现幂函数可以有效地描述两者的关系：

与式（11）相比，式（12）对于土壤侵蚀速率具有更高的表达能力（R2=0.795，RMSE=0.325）。如图5所示，单位水流功率与土壤侵蚀速率之间存在良好的拟合关系，单位水流功率较小时，式（12）的拟合效果最佳，随着单位水流功率的增大，拟合效果相对变差，表明在水深较浅、流速较小、坡度较缓时，式（12）对土壤侵蚀速率具有较为出色的预测能力，其原因可能与紊流、急流状态下（大坡度、大流量）坡面产沙过程频繁的波动性和突变性有关。单位水流功率是一个基于坡度和流量的综合性参数，是一个受表面粗糙度影响的适用性参数[30]。对于坡面侵蚀而言，尤其是细沟侵蚀，与水深相比，流速更易精确地获取，因此在实际应用中使用单位水流功率描述土壤侵蚀速率优于水流剪切力和水流功率。

图5 单位水流功率与土壤侵蚀速率的定量关系 Fig.5 Quantitative relationship between unit stream power and soil erosion rate

Wirtz 等[31]认为单独的水力学参数不足以用来预测土壤分离速率，应考虑到土壤性质对细沟侵蚀过程的影响。Ali 等[32]在建立泥沙输移能力预测方程时引入了土壤的中值粒径D50以表征土壤性质。Kinnell[33]指出分选良好的试验材料在坡面输移过程中受土壤黏聚力C的影响。由于土壤颗粒间黏聚力的阻力，与非黏性土相比，黏性土不易被径流剥离[7]。相比其他影响因素，对土壤黏聚力影响最大的因子是土壤机械组成[34]。基于此，本研究尝试在水力学参数与土壤侵蚀速率定量关系的基础上，引入表征土壤性质的土壤黏聚力参与建立基于土壤性质与水力学参数的定量关系，用于评估坡面的土壤侵蚀速率，见式（13）、式（14）：

对比分析上述式（11）～式（14）可知，在引入土壤黏聚力后定量方程的拟合效果得到了提高，对土壤侵蚀速率具有更好的解释能力。在方程中，黏聚力的指数呈负数，土壤黏聚力对土壤的侵蚀程度呈消极作用，越黏重的土壤，其抗蚀性越强。与应用单独的水力学参数发展的预测模型相比，本研究考虑到了土壤性质对土壤侵蚀过程的定量影响，在预测土壤侵蚀速率时方程的可靠性与适用性更强[11-12]。本研究旨在探讨更具实际意义的土壤侵蚀速率预测方程，但限于试验工作量等因素，试验中只选取了5 种质地的土壤材料，试验材料覆盖范围不足，进一步研究需要考虑更多质地类型的土壤。

3 结 论

本研究在2 个坡度（5°、15°）和4 个流量（2、4、6、8 L·min-1）条件下，分析研究了5 种质地的重塑土坡面的土壤侵蚀速率，建立了水力学参数、土壤性质与土壤侵蚀速率的定量关系。研究得出以下主要结论：（1）坡度和流量对土壤侵蚀速率有显著的积极影响，相对而言，坡度对土壤侵蚀速率的影响较大；土壤侵蚀速率随着土壤质地的变化（黏土至砂土）呈现单峰状分布，在50%含沙量时最大；（2）细沟断面宽深比随着坡度、流量及土壤质地的变化表现出明显的规律性，与土壤侵蚀速率之间具有良好的相关关系，在应用中可以一定程度地反映坡面的侵蚀程度；（3）单位水流功率、水流功率、坡度、水流剪切力、流速、雷诺数及弗汝德数与土壤侵蚀速率具有显著的相关性，其中，单位水流功率是最适合描述土壤侵蚀速率的水力学参数；（4）基于单位水流功率与土壤侵蚀速率的定量关系可以用幂函数表达，在表达式中引入土壤黏聚力后，可以得到更加精确、可靠的定量关系，考虑到参数获取方式的简便，在实际应用中具有较广的适应范围与现实价值。本研究为土壤侵蚀预测模型的发展提供了参考依据，对工程堆积体、崩岗崩积堆等类土沙混合体的土壤侵蚀防治工作也具有指导意义。