不同流量級組合下細溝侵蝕產(chǎn)沙機理研究

許強強，張寬地，，任 涓

（1.西北農(nóng)林科技大學黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室，陜西 楊凌 712100；2.西北農(nóng)林科技大學水利與建筑工程學院，陜西 楊凌 712100）

0 引 言

細溝侵蝕是黃土高原地區(qū)主要的土壤侵蝕方式之一，系指在坡面形成跌坑和細小溝道的情況下，徑流對其分散、沖刷和搬運的過程；其侵蝕量占坡面侵蝕量的70%，占總侵蝕量的45.3%，其發(fā)生取決于坡面流水動力學特性及下墊面特征[1-3]。

不同雨型時坡面流水動力學特性存在明顯差異，進而對坡面土壤侵蝕產(chǎn)生顯著影響[4-6]。秦鳳等[4]研究了紫色土地表延遲型和減弱型兩種雨型下的侵蝕量，得出減弱型降雨侵蝕量比延遲型低26.38%。安娟等[5]基于不同雨型下的模擬降雨試驗，對褐色土橫壟坡面的產(chǎn)沙情況進行了定量研究，發(fā)現(xiàn)侵蝕量最小的雨型是增加型，減弱型、增加-減弱型和減弱-增加型產(chǎn)沙量較增加型分別增加3.03%、43.03%和10.30%，并且在細溝發(fā)育階段增加幅度最大。同一雨強處于降雨過程的不同階段，其對土壤侵蝕的貢獻率具有差異性[6]。鄔鈴莉等[6]通過對北方土石山區(qū)棕壤坡面2013-2015年的105 場自然降雨進行水文調(diào)查，發(fā)現(xiàn)隨著最大雨強出現(xiàn)的時間不同，侵蝕量由小到大的雨型分別是均值型、遞減型、峰值型和遞增型。可見雨型對坡面土壤侵蝕有重要的影響，其中在細溝侵蝕階段的影響最大。不同雨型下，現(xiàn)有的研究已經(jīng)對多種不同土質(zhì)（紫土、褐土和棕壤土）的坡面土壤侵蝕過程進行了探討，已取得豐碩的研究成果[7]。然而，坡面土壤侵蝕較為嚴重的黃土高原黃綿土地區(qū)，缺乏不同雨型條件下的研究成果來指導生產(chǎn)實踐，防止暴雨對坡耕地品質(zhì)產(chǎn)生影響。根據(jù)黃土高原地區(qū)的暴雨特性及其分布規(guī)律，暴雨的基本雨型可以歸納為三類，即猛降型、增加型和間歇型，尤其以增加型和間歇型雨型危害最大[8]。本研究將闡明黃土高原地區(qū)黃綿土土質(zhì)坡面在不同雨型（遞增型和間歇型）對應流量類型下的細溝侵蝕產(chǎn)沙過程。

土壤侵蝕主要包括土壤分離、泥沙輸移和泥沙沉積3個過程。土壤剝蝕過程主要包括邊壁掏蝕、下切侵蝕和溯源侵蝕等，是將土壤顆粒(包括團聚體)與土體/土壤母體分離變?yōu)槟嗌车倪^程，該過程與水流強度密切相關。土壤剝蝕率是表征土壤剝蝕程度的一個重要參數(shù)，主要受到含沙量、細溝長度、徑流強度、坡度等的影響[9,10]。研究土壤剝蝕率與坡面水流水動力學參數(shù)的關系可以將水動力特性和坡面土壤特性聯(lián)系起來，更好地認識坡面土壤侵蝕過程。由流量、坡度、流速及水深等水動力學參數(shù)集成的水流功率、單位水流功率和水流剪切力等水流強度參數(shù)，具有較強的綜合性[11]。恒定流工況下，相關學者對土壤剝蝕率和水流強度參數(shù)（水流剪切力、水流功率、單位水流功率、單寬能耗等）的關系研究取得了豐富的成果[12]。美國WEPP 模型、歐洲EUROSEM 和LISEM 模型以及澳大利亞GUEST 模型所采用的水蝕因子分別為水流剪切力、單位水流功率以及水流功率。肖海等[13]基于室內(nèi)放水沖刷試驗，發(fā)現(xiàn)黃綿土土壤剝蝕率與水流剪切力、水流功率和單位水流功率均能較好地建立對應的擬合關系；同時得出對土壤剝蝕率影響程度最大的侵蝕因子是水流功率。而Li 等[14]認為與紫色土土壤剝蝕率顯著相關的水蝕因子是水流剪切力和水流功率，其較單位水流功率和單寬能耗更能準確預測紫色土土壤剝蝕率。Niu 等[15]認為粉砂壤土土壤剝蝕率與水流功率顯著相關，而與水流剪切力、單位能量和徑流動能等水蝕因子相關性較差。雖然目前對土壤剝蝕率與水流強度參數(shù)之間的耦合關系研究較為豐富，但主要集中于恒定流工況，且沒有得出統(tǒng)一的結(jié)論[13-15]?？梢姡捎谕寥捞卣?、徑流條件等方面的差異，土壤剝蝕率與水流強度參數(shù)之間的耦合關系還需進一步深入研究。

綜上所述，本研究旨在通過變流量下的徑流沖刷試驗，分析黃土高原地區(qū)黃綿土土質(zhì)坡面細溝侵蝕產(chǎn)沙機理：①探索不同流量類型下的細溝侵蝕產(chǎn)沙量；②探究各流量類型在不同沖刷歷時的侵蝕產(chǎn)沙特征；③量化各流量類型條件下的土壤剝蝕率與水流強度參數(shù)的關系；④明確各流量類型條件下的土壤剝蝕率與坡面流能耗的關系；⑤計算各流量類型條件下的細溝流挾沙系數(shù)和判別細溝流挾沙的飽和程度，可望為黃土高原黃綿土侵蝕地區(qū)侵蝕預報模型的建立和侵蝕防治提供理論支持。

1 材料與方法

1.1 試驗土壤

本試驗土壤為黃綿土，取自陜西安塞（109°19′23″E，36°51′30″N），土壤質(zhì)地為粉質(zhì)壤土。表1為試驗土粒徑分布特征[16]。

表1 試驗土粒徑分布特征Tab.1 Particle size distribution characteristics of test soil

1.2 試驗設計

（1）本次試驗通過在土槽內(nèi)設置徑流小區(qū)來進行沖刷試驗。土槽長、寬和高分別為6.0 m、0.4 m 和0.5 m，坡度可調(diào)范圍為0°～15°（圖1）。

圖1 峰值型流量類型沖刷下坡面發(fā)育過程Fig.1 Development process of downslope surface under peak flow type

（2）為模擬自然坡面水文條件，土槽底部布置有梅花形孔；在土槽底部裝填20 cm 厚的細砂，細砂之上鋪設兩層土工布，繼續(xù)分3 層裝填30 cm 厚的試驗土。試驗土裝填前需烘干并過10 mm 的篩以剔除碎石等雜物；試驗土分層裝填時平整壓實，并用環(huán)刀法隨時測量其堆積密度[17]，控制其為1.20 g/cm3[18]。試驗土裝填完成后，對其進行充分灑水養(yǎng)護，保證其達到飽和含水率。

（3）為避免野外降雨試驗流量等難以控制和測量，提高試驗精度，本次試驗采用室內(nèi)徑流沖刷。根據(jù)黃土高原地區(qū)侵蝕性降雨的強度標準（I15≥0.852 mm/min，I10≥1.055 mm/min，I5≥1.520 mm/min）[19]，設計3 個放水流量（2.5 L/min，5.0 L/min，7.5 L/min）。流量由蠕動泵提供（YZ35 型號，ZG600x 泵頭）。通過調(diào)研黃土高原地區(qū)天然降雨類型和地表徑流狀況，增加型和間歇型降雨類型對該地區(qū)的危害較大，而間歇型可近似分為凹陷型和峰值型[6]，所以增加型、凹陷型、峰值型和均勻型（對照組）等4種流量組合類型在該地區(qū)具有很好的代表性。值得注意的是，天然降雨雨強是漸變的，長歷時內(nèi)波動較大，而短歷時內(nèi)在較小范圍內(nèi)波動，因此本文借鑒五點概化過程線法原理，利用三階段放水沖刷過程來模擬徑流流量和徑流總量。故本次試驗設計了增加型（2.5-5.0-7.5 L/min）、凹陷型（5.0-2.5-7.5 L/min）、峰值型（2.5-7.5-5.0 L/min）和均勻型（5.0-5.0-5.0 L/min）4 種流量類型。不同流量類型各沖刷階段（起始、中間和結(jié)束階段）的流量布置如表2所示。根據(jù)中國科學院黃土高原綜合科學考察隊對黃土高原地區(qū)坡耕地統(tǒng)計結(jié)果，坡度小于15°的耕地占耕地總面積的72.16%，因此本文選取4°、8°、12°等3 個較緩坡度進行研究[20]。

（4）本試驗于2021年4月8日至6月1日在西北農(nóng)林科技大學黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室降雨大廳完成。本試驗為3 個坡度和4 個不同流量類型的組合試驗，3次重復，共36場。每場試驗內(nèi)由2臺蠕動泵持續(xù)提供3個流量級，即第1 臺蠕動泵提供起始階段和結(jié)束階段流量，第2 臺蠕動泵提供中間階段流量，一臺蠕動泵關閉的同時開啟另外一臺蠕動泵（表2）。沿土槽下部設置觀測斷面5個，分別距離坡頂0.5 m、1.5 m、2.5 m、3.5 m 和4.5 m，并在各斷面處用精度為1 mm 的鋼尺測量水深和水寬。同時采用高錳酸鉀顏色示蹤法測量表面流速。流速由流速修正系數(shù)乘以表面流速得到，本文流速修正系數(shù)取0.7[21]。徑流沖刷時，每2 min在徑流小區(qū)出口用量筒收集含沙水流，記錄時間和體積。含沙水流樣本烘干后稱重可得產(chǎn)沙量[22]。本次主要研究對象為單條細溝侵蝕產(chǎn)沙過程，故為引導細溝發(fā)育，在試驗土表面布置細溝雛形（圖1）。細溝雛形長6.0 m，沿觀測斷面剖開為等腰梯形；該梯形頂寬為10 cm，底寬為8 cm，深度為1 cm。在土槽頭部設置有與土槽等寬的駝峰堰，水流經(jīng)過駝峰堰穩(wěn)定流入土槽床面，經(jīng)過床面匯流區(qū)域后（圖1），均勻流入細溝雛形[23]。一場試驗結(jié)束將試驗土挖除，重新填充，進行下一場次試驗。

表2 流量類型Tab.2 Discharge flow types

1.3 計算公式

徑流剪切力[24]為徑流分離土壤顆粒和輸移泥沙的沖刷動力，公式為：

式中：τ為徑流剪切力，Pa；γ為水流容重，N/m3；g為重力加速度，m/s2；R為徑流水力半徑，m；J為水力坡度。

水流功率[25]為單位面積水體勢能隨時間的變化率，公式為：

式中：W為水流功率，N/(m·s)；q為單寬流量，m3/(m·s)；h為水流深度，m；U為斷面平均流速，m/s。

有效水流功率[25]為水流功率與臨界水流功率（水流開始剝蝕床面的水流功率）之差，公式為：

式中：φ為有效水流功率，N/(m·s)；Wc為臨界水流功率，N/(m·s)；τr為起動切應力，Pa；ur為泥沙起動流速，m/s；γs為土壤顆粒容重，N/m3；θ為坡度；γb為床面泥沙容重，N/m3；γb0為床面泥沙壓實容重，N/m3；A為常數(shù)，本文取2.9×10-3g/cm3；μ為機械阻力系數(shù)；d為沙粒粒徑，mm，本文取d=d50；ε為坡面土壤孔隙率，黃綿土為52%。

消能率[26]為徑流總能量的變化率，公式為：式中：ΔE為消能率；E和E0分別為水流進出口斷面總能量，m；v1和v2分別為進出口斷面處平均流速，m/s；h1和h2分別為進出口斷面處平均水深，m；H為進出口斷面處的高度差，m；a1和a2分別為兩斷面的動能修正系數(shù)，本次近似取為1.0；g為重力加速度，m/s2。

土壤剝蝕率[25]為單位時間、單位面積上被剝蝕的土壤質(zhì)量，公式為：

式中：Dr為土壤剝蝕率，g/(m2·s)；M為泥沙樣本烘干后質(zhì)量，g；T為泥沙樣本接取時間，s；L為坡長，m；b為徑流寬度，m。

挾沙能力[27]為一定水力和邊界條件下，水流能夠輸移的泥沙量，公式為：

式中：Tc為挾沙能力，kg/m3；ω為泥沙顆粒沉速，m/s；K為挾沙系數(shù)，表示水流的挾沙狀態(tài)，與床面的供沙條件有關；n為指數(shù)（J＞0.002 5，n＝3）；v為水流的運動黏滯性系數(shù)，m2/s；k1本文取值2.0；k2本文取值3.0。

2 結(jié)果與分析

2.1 各流量類型工況下的土壤剝蝕率

不同流量類型沖刷下，不同坡度坡面土壤剝蝕率變化情況如圖2所示。由圖2可以看出，同一坡度下，土壤剝蝕率最大的流量類型為峰值型，均勻型流量沖刷下的剝蝕率次之（除坡度4°外）；峰值型、凹陷型和增加型的土壤剝蝕率分別為均勻型的1.36倍、0.91倍、0.95倍。因此，用恒定流（均勻型流量類型）來研究黃土高原地區(qū)黃綿土土質(zhì)坡面細溝侵蝕實際產(chǎn)沙量時，需在恒定流產(chǎn)沙量的基礎上乘以1.36 的安全系數(shù)，即能得到較為準確的土壤剝蝕率。而溫磊磊等[28]認為黑土土質(zhì)土壤剝蝕率最小的雨型是峰值型，其次是均勻型。本文與溫磊磊等[28]得出完全相反的結(jié)論，可能是本文試驗土為黃綿土，并且為變坡沖刷試驗，兩者在試驗設計和布置上存在差異。因此，除流量類型外，土壤的物理性質(zhì)也是造成細溝侵蝕差異的因素[29]。

圖2 不同坡度下各流量類型產(chǎn)沙特征Fig.2 Sediment production characteristics of different flow types under different slopes

隨著坡度由4°增加到12°，峰值型、增加型、均勻型和凹陷型流量類型下的土壤剝蝕率均有不同程度的增加，分別增加131.02%、28.37%、280.34%和122.51%；其中增幅最大的流量類型為均勻型，其土壤剝蝕率大小順序由第四增加至第二，增幅最小的流量類型為增加型，其土壤剝蝕率大小順序由第二降為第四。另外，對各流量類型工況下的土壤剝蝕率與坡度進行皮爾遜相關性分析（圖2），發(fā)現(xiàn)二者顯著相關（P＜0.05）。

通過觀察細溝發(fā)育過程也可以得出上述結(jié)論。較低坡度時，跌坑發(fā)育遲緩，跌坑數(shù)量較少（圖1）。坡度逐漸增大到12°時，跌坑的發(fā)育愈加明顯，數(shù)量增多，深度加深（圖1）。水流跌入跌坑的過程中，會產(chǎn)生漩渦，對跌坑和細溝的側(cè)壁進行剝蝕[30]。同時，坡度較大時，坡面臨界切應力顯著減低，坡面土壤易被水流剝離[25]。加之坡度越大，徑流所含的能量越大，紊動性越強，大量的泥沙便會被裹挾成為含沙水流的一部分，運送至徑流小區(qū)的出口[31]。由此可見，坡度對四種流量類型下的產(chǎn)沙過程均有顯著影響[32]。

2.2 各流量類型在不同沖刷歷時的產(chǎn)沙量

流量類型不同，其產(chǎn)沙量顯著不同，因此其在各沖刷階段（起始、中間和結(jié)束階段）的產(chǎn)沙情況也將存在差異。繪制各流量類型不同沖刷階段產(chǎn)沙量及其貢獻率如表3所示。由表3可以看出，凹陷型、峰值型和增加型均在沖刷結(jié)束階段（沖刷歷時為40～60 min 時）產(chǎn)沙量最多，貢獻率范圍為41.26%～62.36%。究其原因，主要是隨著徑流沖刷的持續(xù)，跌坑的出現(xiàn)、發(fā)育和貫穿顯著影響了產(chǎn)沙進程[33]。隨著徑流沖刷的進行，流路上逐漸出現(xiàn)一些跌坑，這些跌坑是侵蝕產(chǎn)沙的最活躍點；由于跌坑上游方向有垂直的陡壁，該區(qū)域被剝離的方式不同于平整坡面，同時受到水流剪切應力和自身重力的雙重影響。此時跌坑處的土粒呈塊狀滾落進入水流，并與水流一起運動，侵蝕產(chǎn)沙較為強烈；逐漸跌坑會演化為細溝的溝頭，持續(xù)發(fā)生溯源侵蝕。同時點狀細溝不斷加長，多處溝頭之間相互貫通，最終在坡面形成完整細溝。坡面細溝產(chǎn)生及發(fā)育的速度隨著坡度、流量類型的變化而明顯不同。在緩坡或凹陷型、增加型流量類型時，跌坑形成所需要的時間較長（主要集中在起始、中間階段），溯源侵蝕速度較慢；而在陡坡或峰值型流量類型時，坡面上很快便會形成一系列的跌坑（主要集中在起始階段），溝頭迅速發(fā)育貫通，形成完整的輸沙通道－細溝。形成細溝后，坡面流進一步集中，導致在沖刷的結(jié)束階段（沖刷歷時為40～60 min時）產(chǎn)沙量最多。

表3 各流量類型不同沖刷階段的產(chǎn)沙量及其貢獻率Tab.3 Sediment yield and its contribution rate in different scouring stages of different flow types

同時表3也反映出，凹陷型和增加型流量類型不同沖刷階段的侵蝕量與其流量變化趨勢保持一致，表現(xiàn)為“水大沙大”[34]。從能量的角度分析，流量越大則徑流的能量越大；而土壤侵蝕的過程需要消耗徑流的能量，可用于侵蝕的能量越多，產(chǎn)沙量隨之增加，土壤侵蝕愈發(fā)嚴重[35]。

與凹陷型流量類型不同，峰值型流量類型在整個沖刷過程產(chǎn)沙量均在增大。主要原因包括：①起始階段的放水流量較小，容易在細溝表面形成土壤結(jié)皮，增加了土壤的抗蝕性，此時產(chǎn)沙量較小[28]。這種土壤結(jié)皮在大流量下很容易被破壞，土壤抗蝕性將顯著降低。因此，中間階段的7.5 L/min 最大流量首先會破壞坡面表面形成的結(jié)皮，加速對細溝溝槽、細溝內(nèi)壁的掏刷，同時擁有較強的水流挾沙力，將泥沙運送至徑流小區(qū)出口，侵蝕量增大。②結(jié)束階段流量有所減小但仍大于起始階段，隨著中間階段7.5 L/min 最大流量的沖刷，細溝斷面由寬淺型過渡到窄深型[26,36]，過水斷面顯著減小，細溝水流進一步集中。由水流連續(xù)性方程可得，徑流過水斷面減少的同時其流速必然增加。過水斷面減小造成的流速增加量補充了流量降低造成的流速減少量，并且有盈余，所以結(jié)束階段流速仍大于中間階段流速，結(jié)束階段產(chǎn)沙量大于中間階段產(chǎn)沙量。此外，坡度的增大也不同程度的加劇了各階段的侵蝕作用。

2.3 各流量類型工況下的土壤剝蝕率與水流強度參數(shù)的關系

由前述已知，不同流量類型的土壤剝蝕率存在顯著差異，各沖刷階段（起始、中間和結(jié)束階段）產(chǎn)沙量并不都與流量保持一致。因此，土壤剝蝕率與各水流強度參數(shù)之間的變化關系在不同流量類型工況下也存在差異。基于此，本文點繪了土壤剝蝕率與各水流強度指標間的分布關系如圖3所示。由圖3可以看出，土壤剝蝕率整體上隨水流強度參數(shù)的增加而增加。然而水流強度參數(shù)在不同范圍時的土壤剝蝕率增長幅度不同；當水流強度參數(shù)較小時，數(shù)據(jù)點相對密集，土壤剝蝕率隨著水流強度的增加而迅速增加，水流剪切力表現(xiàn)為小于3 Pa，水流功率表現(xiàn)為小于1.0 N/(m·s)，有效水流功率表現(xiàn)為小于0.2 N/(m·s)；而當水流強度參數(shù)超過上述范圍時，土壤剝蝕率將隨著水流強度參數(shù)的增加緩慢增大。這與肖海等[13]得出的結(jié)論相一致，說明上述水流強度參數(shù)仍然對不同流量類型工況下的細溝侵蝕產(chǎn)沙有很好的預測效果。

圖3 土壤剝蝕率與水流剪切力、有效水流功率和水流功率的關系Fig.3 Relationship between soil erosion rate and water flow shear force,effective water flow power and water flow power

由表4可知，土壤剝蝕率與水流強度參數(shù)之間的關系均能很好地用冪函數(shù)描述。峰值型流量類型下水流強度參數(shù)擬合公式系數(shù)項（4.81～29.80）和決定系數(shù)R2（0.77～0.94）均最大，而均勻型條件下決定系數(shù)R2（0.37～0.59）最小，系數(shù)項較?。?.06～21.96）。這表明峰值型流量類型工況下的徑流單位能量的產(chǎn)沙量最高，對土壤有較強的剝蝕、分散和搬運能力，且二者之間相關性顯著。值得注意的是，這并不能表明均勻型流量類型工況下的徑流對土壤的剝蝕、分散和搬運能力較弱。目前學者對均勻型（恒定流）流量類型侵蝕產(chǎn)沙過程的研究較為明晰[15]。徑流沖刷開始時，含沙量很快增加到最大值，在維持高位值一定時間后，快速降為較低的含沙量，此時徑流沖刷達到了動態(tài)平衡。上述產(chǎn)沙過程約占沖刷全歷時的70%，剩下的沖刷歷時內(nèi)，幾乎不再產(chǎn)沙。本文均勻型流量類型侵蝕產(chǎn)沙過程與前人的研究成果完全相同，但是本文在數(shù)據(jù)處理時，對均勻型流量類型整個放水過程進行了研究，導致產(chǎn)生較大誤差。因此，沖刷歷時對侵蝕產(chǎn)沙預測模型的建立有重要影響。

表4 土壤剝蝕率與水流強度參數(shù)的擬合方程Tab.4 Fitting equation of soil erosion rate and water flow intensity parameters

均勻型流量類型徑流在沖刷結(jié)束階段（沖刷歷時為40～60 min時）產(chǎn)沙最少，主要表現(xiàn)在以下兩個方面：一是水流強度參數(shù)的降低，二是土壤抗侵蝕能力的增強。均勻型流量類型徑流在侵蝕的階段，促使細溝斷面發(fā)育，此時徑流的過水斷面顯著增加[20]。由水流連續(xù)性方程可知，當流量一定時，過水斷面面積增加，必然會導致水流流速的降低。此時集中流發(fā)散，水流強度指標降低。同時，細溝深度越深，其容重越大，進而抗侵蝕能力增強[25]。

2.4 各流量類型工況下的土壤剝蝕率與消能率的關系

土壤侵蝕的3個子過程中，土壤分離、泥沙輸移和泥沙沉積均需消耗能量[13]，而細溝侵蝕產(chǎn)沙及微地貌產(chǎn)生是其外在表現(xiàn)。因此，需從能耗角度探究其對土壤剝蝕率的影響；繪制了土壤剝蝕率和能耗關系曲線如圖4所示。由圖4可知，各流量類型的土壤剝蝕率與消能率均呈正相關，且峰值型的土壤剝蝕率與消能率間的相關性最強，其次是增加型和均勻型，凹陷型最弱。

不同流量類型之間土壤剝蝕率與消能率相關關系存在差異，可能是由能耗內(nèi)部構(gòu)成差異造成的。單位耗能率或徑流切應力是從床面阻力或系統(tǒng)耗能的角度出發(fā)來研究水流的侵蝕產(chǎn)沙。白清俊[36]經(jīng)過物理公式推導認為徑流切應力對單位水體所做的功就等于徑流的單位耗能率。單位水體水流在沿流道向下游流動時，其能耗主要用來克服阻力做功，能耗的結(jié)果主要表現(xiàn)在兩個方面，一是部分能耗轉(zhuǎn)化為熱能由水體表面散發(fā)出去；二是通過邊界阻力的擾動作用增加了水流的紊動性，紊動的結(jié)果使水體內(nèi)局部瞬時流速增大而具有較大的能量，但這些能量強點很快向四周散發(fā)自身能量而調(diào)節(jié)整體水流速度。由表4和圖4可知，土壤剝蝕率與消能率（水流剪切力）的擬合總是存在一定的誤差。事實上，水流的侵蝕挾沙與轉(zhuǎn)化為熱能散發(fā)出去的能耗沒有關系，只有轉(zhuǎn)化為紊動的那部分能耗才與水流的侵蝕挾沙密切相關。凹陷型與均勻型流量類型工況下的土壤剝蝕率與消能率（水流剪切力）之間的擬合決定系數(shù)僅為0.47（0.61）和0.55（0.59），而增加型與峰值型流量類型工況下的土壤剝蝕率與消能率（水流剪切力）之間的擬合決定系數(shù)可達到0.76（0.70）和0.84（0.94）。土壤剝蝕率與消能率（水流剪切力）之間的擬合決定系數(shù)越大，表明能耗中用于侵蝕挾沙的能量越多。這就說明，水體表面散發(fā)出去的能耗和增加水流紊動性的能耗大小不是相對固定的，存在此消彼長的過程。因此，研究增加型與峰值型流量類型工況下的土壤剝蝕率與消能率之間的關系是有重要意義的，而在凹陷型與均勻型流量類型采用單位耗能率來分析細溝流的侵蝕挾沙有較大誤差產(chǎn)生，因此凹陷型與均勻型流量類型工況下的土壤剝蝕率與能耗之間的關系還需后續(xù)進一步的研究。

圖4 各流量類型工況下土壤剝蝕率與消能率的關系Fig.4 Relationship between soil erosion rate and energy dissipation rate under different flow types

2.5 各流量類型工況下的細溝水流挾沙系數(shù)

沙玉清[27]在計算水流挾沙能力時，采用了有效流速類指標，見公式（9）。公式（9）里的K定義為挾沙系數(shù)。K是一個綜合系數(shù)，在一定程度上反映的是含沙水流挾沙的飽和狀態(tài)以及細溝侵蝕過程中床面底部的供沙特性。高飽和時可取K值為400，對應不淤保證率為15.9%；中飽和時可取K值為200，對應不淤保證率為50%；低飽和時，K值可取為91，相應不淤保證率為84.1%。黃土高原地區(qū)挾沙系數(shù)相對較高，如渭惠渠平均粒徑為0.02 mm，最大極限挾沙系數(shù)為3 390，渭河懸移質(zhì)平均粒徑為0.026 mm，平均挾沙系數(shù)為408[27]。

各流量類型工況下的試驗實測細溝水流含沙量見表5。將公式（6）、（10）代入公式（9）中，利用實測含沙量來反推試驗中細溝侵蝕的挾沙系數(shù)K，即可判別坡面細溝挾沙水流的飽和程度，計算結(jié)果見表5。從表5可以看出，各流量類型工況下，挾沙系數(shù)隨著坡度的增加而增加。各流量類型的挾沙系數(shù)在低坡度時相差較大，而在高坡度時二者差異較小。各流量類型工況下挾沙系數(shù)K變化范圍為：20.84～83.57，取其平均值52.21，相應的不淤保證率為95%左右，即細溝流處于低飽和狀態(tài)。

表5 各流量類型工況下的細溝流挾沙系數(shù)Tab.5 Sediment-carrying coefficient of rill flow under different flow types

3 結(jié) 論

通過室內(nèi)放水沖刷試驗，研究了不同流量類型工況下黃土高原地區(qū)黃綿土土質(zhì)坡面侵蝕產(chǎn)沙過程，主要結(jié)論為：

（1）相同坡度時土壤剝蝕率最大的流量類型為峰值型；均勻型流量類型沖刷下的土壤剝蝕率次之（除坡度為4°外）。土壤剝蝕率與坡度顯著相關（P＜0.05）；峰值型、增加型、均勻型和凹陷型流量類型工況下的土壤剝蝕率均隨坡度的增加而增加，其中均勻型流量類型工況下的土壤剝蝕率增幅最大。

（2）4 種流量類型產(chǎn)沙過程中產(chǎn)沙量最多的階段為沖刷結(jié)束階段，產(chǎn)沙貢獻率范圍為41.26%～62.36%。凹陷型和增加型流量類型不同沖刷階段的侵蝕量與其流量變化趨勢保持一致，而峰值型流量類型在整個沖刷過程的產(chǎn)沙量均在增大。

（3）土壤剝蝕率隨著水流強度的增加首先迅速增加，超過一定范圍時，增速放緩。土壤剝蝕率與水流強度參數(shù)之間的關系均能很好地用冪函數(shù)描述（R2=0.37～0.94）。峰值型流量類型下水流強度參數(shù)擬合公式系數(shù)項（4.81～29.80）和決定系數(shù)R2（0.77～0.94）均最大，單位能量的產(chǎn)沙量最高，二者間的相關性最強。沖刷歷時對均勻型流量類型細溝侵蝕產(chǎn)沙預測模型的建立有重要影響。

（4）各流量類型工況下的土壤剝蝕率與消能率均呈正相關（冪函數(shù)關系），且峰值型的土壤剝蝕率與消能率間的相關性最強，其次是增加型和均勻型，凹陷型最弱。

（5）細溝流挾沙系數(shù)K變化范圍為：20.84～83.57，平均值為52.21，相對應的不淤保證率為95%，細溝流處于低飽和狀態(tài)。