坡面細溝侵蝕斷面形態發育影響因素分析及動力特性試驗

馬小玲，張寬地,，楊 帆，楊明義，董 旭

?

坡面細溝侵蝕斷面形態發育影響因素分析及動力特性試驗

馬小玲1，張寬地1,2※，楊 帆1，楊明義2，董 旭1

（1. 西北農林科技大學水利與建筑工程學院，楊凌712100； 2. 中國科學院水利部水土保持研究所黃土高原土壤侵蝕與旱地農業國家重點實驗室，楊陵 712100）

研究細溝形態發育過程對認識細溝侵蝕具有重要作用，該文采用6種坡度（2°、4°、6°、8°、10°、12°），5種流量（8、16、24、32、40 L/min）下的組合沖刷試驗，系統研究了坡面細溝橫縱斷面形態發育影響機制及動力特性。結果表明：細溝寬深比變化范圍為3.006～4.884，根據水力最佳斷面，細溝水流遠未達到穩定。橫斷面形態系數隨坡度的變化范圍為0.36～0.522，細溝橫斷面形態隨流量、坡度以及沖刷歷時均趨近于梯形水力最佳斷面，即阻力最小的斷面。隨著流程長度的增加，橫斷面形態由寬深逐漸變窄，橫斷面形態系數也隨之減小。細溝縱斷面形態范圍為0.60～11.26，且隨坡度的增大而增大，與流量相關性不大。綜合阻力系數及消能率均與細溝縱斷面形態系數呈良好的冪函數關系。

侵蝕；試驗；土壤；細溝形態；形態指標系數；跌坑；消能率

0 引 言

坡面土壤侵蝕主要分為細溝間侵蝕和細溝侵蝕2種方式。而細溝侵蝕產沙占總輸沙量的50%～70%，細溝侵蝕有別于其他溝道侵蝕的一個顯著特點是伴隨著侵蝕產沙過程，細溝的形態演化迅速。因此細溝形態的演化與水流水動力學特性及土壤侵蝕間呈互饋過程[1-2]。而目前關于細溝形態的研究主要集中在其與侵蝕因子的響應關系上。吳普特等[3]采用稀土元素示蹤技術對細溝的垂直分布特征進行了研究，并采用平均密度和平均深度來表征細溝形態。Bruno等[4]根據細溝長度和橫斷面發展來描述細溝形態，并得出細溝長度與體積之間呈冪函數變化。和繼軍等[5]選用細溝密度和細溝寬深比等指標表征了楊凌塿土和安塞黃綿土的細溝侵蝕特征差異。沈海鷗等[6]選取細溝傾斜度、細溝密度、細溝割裂度和細溝寬深比5個指標，分別分析各指標與細溝形態的相應關系，并提出表征細溝形態的最優指標。除了以上所述形態指標，細溝縱斷面形態的另一個宏觀表現是跌坑的發育，跌坑是坡面細溝開始發育的標志，關于跌坑的形成機理主要有以下幾種觀點：一種認為坡面立地條件差異會影響跌坑的形成，如地形和土壤特性的差異等[7-8]。嚴冬春等[9]通過試驗研究發現，由于土壤顆粒組成成分的不同，紫色土坡面形成的跌坑邊界要比黃土的更為模糊。另一種觀點是從能量的角度分析了跌坑的形成機制，當徑流能量積累到一定程度，坡面侵蝕方式將發生改變，坡面跌坑開始發育[10-13]。坡面上一系列的跌坑是徑流侵蝕最為劇烈的地方，并以此為基準不斷地開始溯源侵蝕。兩相鄰跌坑甚至也會相遇貫通。跌坑的發育是細溝為維持自身穩定而在水流方向做出的形態調整，跌坑能夠集中坡降，增大坡面形狀阻力，并通過水躍來實現高消能率，減少坡面下切侵蝕，從而維持坡面穩定。其消能機理與臺階式溢洪道相似，即通過各級臺階制造跌水，消耗下泄的水流能量，能夠有效地預防河床底面掏蝕，保護下游河床穩定[14]。但目前對跌坑的研究僅限于宏觀描述，沒有定量地研究。

目前關于細溝形態的研究僅僅局限于地學描述以及基本的寬深變化，而對橫縱斷面形態發育過程研究較少。在徑流侵蝕作用下，坡面侵蝕形態發生改變，從而使得徑流的水力學特性也隨之發生變化[15]。定性描述細溝形態，揭示其空間分布規律，對于深刻認識細溝侵蝕規律、促進土壤侵蝕模擬技術發展具有十分重要的意義。

1 材料與方法

1.1 試驗材料

試驗在中國科學院水利部水土保持研究所黃土高原土壤侵蝕與旱地農業國家重點實驗室進行。試驗所用土壤為楊凌當地的塿土，對試驗土壤風干測其機械組成，結果見表1。

表1 土壤顆粒機械組成

1.2 試驗設計

1）裝土前在試驗土槽底部裝填20 cm厚的細沙，細沙上鋪設土工布2層，土槽底部采用梅花形布置打孔，以模擬天然土壤透水狀況，填土采用分層裝填的方法，每層10 cm，共30 cm，根據實際野外土壤狀況，土壤平均干體積質量控制在1.09 g/cm3。

2）試驗小區采用可調坡度鋼槽，小區結構尺度長×寬×深＝6 m×0.4 m×0.6 m，試驗段有效長度為5.8 m。

3）坡度大于5°～6°的坡耕地會發生明顯的土壤侵蝕現象，而小于5°～6°的坡耕地亦有輕微的土壤侵蝕。考慮到陡坡耕地實行退耕還林政策，試驗主要為緩坡耕地土壤侵蝕提供依據。試驗取土地點耕地坡度較緩，故本次試驗選取2°、4°、6°、8°、10°、12°，共6個坡度進行研究，以期為緩坡耕地土壤流失治理提供依據。

4）放水流量沖刷試驗應與野外實際相近似，由于室內試驗土槽尺寸的限制，因此需要考慮上方匯水因素。放水流量由雨強、上方匯流面積、徑流系數等確定。根據野外實際降雨情況及匯流面積、雨強特征，設計試驗沖刷流量為8、16、24、32、40 L/min。

5）為了消除試驗誤差，共設置5個觀測斷面，以槽頂為零點，向下取長度0+0.5、0+1.5、0+2.5、0+3.5、0+4.5 m，并在槽側作出標記，每段測流長度為1.0 m，采用KMnO4染色示蹤法觀測斷面表面優勢流速，采用精度為mm的鋼尺測量各斷面水深和水面寬度（重復測量3次），最終平均化處理；

6）試驗過程中，間隔20 min停止放水，測量跌坑深度（測得跌坑最深處高度和上游結皮段高度，二者的高度差即為跌坑深度）及跌坑長度（鋪設于土槽邊沿精度為mm的米尺）。用測針排架在各斷面處測量細溝橫斷面形態，坡面侵蝕劇烈的部位進行加測（如圖1b所示）。

7）試驗歷時80 min，土壤前期含水量控制在8.32%～12.54%，每場試驗2次重復，并取其均值。

a. 細溝發育a.Rill developmentb. 橫斷面形態測量B .Cross-sectional shape measurement

2 橫斷面形態研究

2.1 細溝寬深比變化規律

最小耗能原理是自然界一切運動的普遍原理，基于此理論，可以斷定細溝水流在流動過程中亦向能耗最小的的趨勢發展。徑流能量受自身流動，剝離泥沙，輸移泥沙等因素共同調配。從細溝侵蝕的發育過程來看，都是有效能量減少、系統趨于穩定的過程。所謂穩定狀態即床面泥沙處于不沖不淤狀態，服從廣義熱力學熵增原理[16]。細溝流道趨于阻力最小的形態發展。在水力學中，水力最佳斷面也是阻力最小，過流量最大的過水斷面。根據明渠均勻流公式

式中為流量，L/min；為過水斷面面積，m2；為謝才系數，m1/2/s；為水力半徑，m；為底坡；為綜合糙率；為濕周，m。

當、、一定時，濕周越小，過水斷面面積越小。而當過水斷面面積一定時，濕周越小，則阻力越小[17]。即水力最佳斷面亦阻力最小的斷面，這一點與徑流對細溝流道的比選不謀而合。沙玉清[18]經過研究得出河道穩定曲線斷面接近邊坡系數=1的梯形。暫且以邊坡系數1:1的梯形斷面為研究對象，由《水力學》[17]知，梯形水力最佳斷面寬深比僅與邊坡系數有關

式中為梯形斷面底寬，m;為梯形斷面深，m;為邊坡系數；β梯形斷面寬深比。

當邊坡系數=1時，寬深比β=0.828。可見，當渠道下墊面相對穩定時，最優水力斷面對應的寬深比為0.828，為探究細溝水流斷面寬深比的演化規律，列出了本試驗條件下不同流量和坡度下的細溝橫斷面寬深比（表2）。

表2 不同流量和坡度下細溝寬深比

由表2可以看出，細溝寬深比均比較大，在3.006～4.884之間變化。與水力最佳斷面寬深比β=0.828相差甚遠，這說明細溝水流系統遠未達到穩定狀態。細溝侵蝕可分為：下切侵蝕、溯源侵蝕、溝壁坍塌，這3種侵蝕方式反映在細溝形態上分別為細溝加深、拓長、加寬。寬深比隨流量和坡度的增加而減小，說明細溝在發育過程中主要以下切侵蝕和溯源侵蝕為主，溝壁坍塌侵蝕作用相對較小。細溝侵蝕產沙過程伴隨著溝道形態的迅速演化。沙玉清[18]認為橫斷面的變形，主要是決定于斷面內流速的分布。當斷面濕周上的各泥沙質點所受作用力和該處泥沙本身的穩定性無法平衡時，沖刷仍將繼續。

2.2 細溝橫斷面形態指標變化

2.2.1 橫斷面形態指標隨沖刷歷時變化

徑流在坡面流動過程中，不斷地對細溝溝壁、溝底及溝頭的土壤進行沖刷和運移，即徑流的側蝕、下切侵蝕及溯源侵蝕作用。不同的侵蝕方式塑造了形態各異的細溝斷面形態，同時，細溝的斷面形態又對細溝侵蝕起制約作用，即細溝斷面形態和侵蝕產沙是相互影響、相互制約的結果。由于細溝發育隨機性較大，流道左右擺動，細溝斷面形態在水流方向上差異較大。為了能客觀評價細溝形態的變化特征，張科利等[19]提出了細溝橫斷面形態指標，即細溝斷面面積與等深同寬的矩形面積之比。圖2為12°條件下細溝橫斷面形態隨沖刷歷時變化規律。

由圖2可見，沖刷初期（0～40 min），溝床主要以下切侵蝕和細溝邊壁坍塌為主，斷面不斷加寬拓深。是因為沖刷初期，水流水層較薄，侵蝕能力有限，隨著細顆粒土壤的啟動，粗顆粒泥沙逐漸暴露，床面開始粗化，水流逐漸集中，產生束水攻沙現象，溝壁臨空面逐漸增加，邊坡穩定性降低。當沖刷歷時為40 min后，細溝寬度基本不再變化，床面以下切侵蝕為主。經過水流的不斷沖刷，細溝由上而下呈階梯狀，水流在各級階梯處跌落，重力作用使得階梯底部不斷加深。沖刷初期（0～40 min），細溝橫斷面形態指標增加了19.3%，而當沖刷歷時為40 min后，橫斷面指標只增加了10.4%。在流量較小時，流量對細溝的溝深有較大影響，流量為8、16L/min時，最大溝深分別為14和19 cm。而當流量大于16 L/min后，隨著流量的增加，最大溝深均為18 cm左右。

圖3為24 L/min流量下，不同坡度下的細溝橫斷面輪廓圖。由圖3可以發現，隨著坡度的增大，細溝橫斷面也隨之加寬拓深，斷面形態也愈加不規則，細溝橫斷面指標隨坡度增大，變化范圍為0.36～0.522。原因是當坡度增大時，徑流沿流動方向的重力分力增加，促使徑流流速增加，最終導致橫斷面不規則變化。同時根據圖2和圖3，橫斷面形態隨歷時和坡度的變化都是由倒三角形向倒梯形轉變，故前面以梯形的水力最佳斷面來分析細溝斷面是可行的。

細溝橫斷面形態變化多樣，如對其進行概化，大致可以分為矩形、倒梯形、倒三角形等形態。當斷面形態近似為矩形時，斷面面積最大，倒三角斷面面積最小，而倒梯形形體斷面面積居中。即值越大，代表細溝斷面形態越接近矩形。值越小，表示細溝斷面形態越接近于倒三角形[19]。經計算，斷面形態為邊坡系數1∶1的梯形水力最佳斷面時，對應的形態指標η為0.646。表3為圖2中各細溝橫斷面所對應的形態指標系數隨時間的變化過程，變化范圍為0.361～0.574。將表3的計算結果與圖2進行對比發現，橫斷面形態指標能夠很形象地描述橫斷面變化過程，指標隨流量和時間的增大大體呈增大趨勢。即值隨流量、坡度以及沖刷歷時均趨近于梯形水力最佳斷面對應的橫斷面指標η。

表3 不同流量下細溝橫斷面形態指標η隨時間的變化

2.2.2 橫斷面形態指標隨坡長變化

在試驗過程中觀察到，細溝侵蝕速率隨著流程長度的增加而逐漸減小。沿細溝縱斷面方向上，按流程距坡頂距離將坡面分為坡上（0～0+1.5 m）、坡中（0+2.0～0+3.5 m）、坡下（0+4.0～0+5.5 m）3個坡段。本次選取侵蝕最為劇烈的坡度12°進行分析。圖4為12°、不同流量條件下坡上、坡中、坡下3個坡段的橫斷面形態。

表4為圖4所對應的各坡段處細溝橫斷面形態指標。

表4 不同流量下橫斷面形態指標η隨坡段的變化

從表4可以觀察到，由坡上段到坡下段，橫斷面指標逐漸減小。即在徑流侵蝕的作用下，隨著流程長度的增加，細溝斷面形態由寬深狀逐漸變窄，同時細溝侵蝕深度也減小。出現這種現象的原因有2個：徑流在沖刷過程中，除了自身流動所需的動能外，還需要剝離和搬運泥沙，徑流從坡頂流下時，首先對坡頂土壤進行剝離，徑流中攜帶了大量的泥沙，徑流因輸運泥沙而消耗了大量能量，以致其流經坡中和坡下部位時，用于剝蝕土壤的能量相對較小，故而對細溝的侵蝕能力減弱[20]；整個沖刷過程中，水流能量的分配是不斷調整變化的，同時被水流輸移的泥沙顆粒和溝床泥沙不斷地進行交替變化。隨著流程長度的增加，徑流對坡面土壤顆粒由侵蝕作用逐漸過渡為侵蝕和沉積交替進行[21]。本文對不同流量和坡度條件下沖刷結束后細溝的各斷面形態進行分析，發現細溝橫斷面形態指標變化范圍在0.29～0.69之間，其變化范圍服從正態分布，并得出各個坡段橫斷面形態指標頻率分布圖（圖5）。由圖5可以看出，相較于坡下部，坡上部和坡中部細溝橫斷面形態指標范圍更大，即坡上部細溝橫斷面形態呈寬深式，隨著距坡頂距離的增加，細溝橫斷面形態逐漸變窄。這與沈海鷗等的研究結果有所不同，沈海鷗試驗侵蝕強度變化為：坡中部>坡下部>坡上部[6]。產生這種差異的原因在于其試驗是在降雨條件下開展的，坡上部匯水面積相對較小，斷面形態為窄淺式。而其試驗中，坡中部位（530～600 cm）處于侵蝕增強帶，坡下（840～860）處于侵蝕減弱帶。

3 縱斷面形態研究

3.1 細溝縱斷面形態指標

坡面侵蝕過程中，隨著徑流的持續沖刷作用，細溝縱斷面方向逐漸出現一系列跌坑，水流在跌坑處發生跌水，并不斷地進行溯源侵蝕，故跌坑的發育是細溝侵蝕的重要環節[22]。為了描述細溝縱斷面方向上跌坑的發育程度，現引出無量綱參數100/，特定義其為細溝縱斷面形態系數。其中為跌坑平均深度，cm；為各組次相鄰兩跌坑最深處沿坡面方向距離的平均值，cm；表5為不同流量和坡度下的細溝縱斷面形態系數。

表5 不同試驗組次縱斷面形態系數δ

本試驗條件下縱斷面形態系數范圍為0.60<<11.26。由表5，隨坡度增加了17.8倍，而隨流量最大增加了1.45倍。即相較于流量，隨坡度的增加更為顯著。是因為縱斷面形態的發育代表了床面侵蝕的劇烈程度，而徑流的侵蝕能力受流量影響不大，主要與坡度有關。試驗過程中觀察到，隨著坡度的增大，床面跌坑平均間距減小，跌坑深度增大，故也隨之增大。從能量角度來分析，徑流在坡面流動過程中，能量不斷地調整變化，當坡度愈大，徑流能量也隨之加大，需要更多的跌坑來消耗水流能量，故而縱斷面形態系數增大。

3.2 縱斷面形態系數與糙率系數之間的響應關系

坡面細溝在水流的沖刷作用下，跌坑與結皮逐漸發育，跌坑和結皮是細溝土壤對水流沖刷起抑制作用的一種手段[23]，徐江等[24]通過水槽試驗來模仿河流階梯-深潭系統，發現河床由粗化層、沙波發展形成階梯-深潭，即通過增加床面阻力來消耗徑流能量。而在階梯-深潭河段，顆粒阻力和形狀阻力平均占總阻力的4%[25]，因此水流阻力主要與跌坑發育程度有關。本次選用糙率系數來描述阻力與縱斷面發育系數的耦合關系。

由圖6a可以看出，綜合阻力系數隨縱斷面形態系數的增加而增大。細溝徑流是細溝侵蝕發育的主要作用力，而坡面形態的發育又對徑流侵蝕有一定的控制作用，故細溝水力特性和細溝形態之間呈相互聯系、相互制約的作用與反作用系統[26]。對于水槽試驗，徑流能量主要消耗在剝離和輸移泥沙顆粒以及坡面結構2方面，且這2部分呈此消彼長的關系[27]。在坡面徑流剛開始沖刷時，坡面較為平整，細溝還未完全發育，徑流開始大量剝蝕土壤顆粒，隨著越來越多的泥沙被輸移出坡面，坡面跌坑開始發育，跌坑數量和深度不斷增加，即增大，坡面形狀阻力也隨之增大。經分析，糙率系數與呈良好的冪函數關系：=0.0240.356(2=0.869 4)。

3.3 縱斷面形態系數與消能率Δ之間的響應關系

根據能量守恒方程，以跌坑上下段能量差與跌坑上游總能量的比值作為消能率。計算公式為

式中Δ為總水頭差，cm; Δ為跌坑深度，cm;為流經跌坑的流量，L/min;1為跌坑上游水深，cm;2為跌坑下游水深，cm; Δ為消能率，%;為重力加速度，m/s2；為細溝溝寬，cm。

徑流在流動過程中，勢能不斷地轉化為動能，水流流入跌坑內時，動能通過水流的強烈紊動轉化為紊動能，能量因此大量耗散。本試驗與消能率Δ的關系如圖6b所示。

由圖6b可知，隨著坡度的增大，消能率增加了15.3%。消能率Δ隨的增大而增大，二者呈冪函數關系增長（2=0.805）。越大，表示跌坑發育愈成熟，即跌坑對于消能有明顯的影響。水流進入跌坑后，水流強烈紊動并摻氣，此階段紊動的主要來源是水流跌入跌坑后所形成的尾流和漩滾作用，在沒有跌坑發育的床面上，水流主要通過床面膚面摩擦來產生紊動。但相對于膚面摩擦，由跌坑產生的紊動消能效果更為顯著[28]。跌坑的這種消能原理與階梯-深潭系統有類似之處，即消能率與階梯高度呈正比，但兩者之間也存在一些差別。階梯-深潭是由山區河流中較大的卵石形成的獨特構造，而一些學者利用水槽試驗模擬階梯-深潭系統時，也會選用粒徑較粗的泥沙顆粒，土壤級配不均勻，如徐江[24]選用泥沙粒徑范圍為1～3 cm。兩者的破壞形式也不同，階梯-深潭是由粒徑較大的卵石所構成的互鎖結構，其穩定性主要靠階梯上最關鍵石塊的受力情況決定，一旦階梯失穩，整個系統都將崩塌破壞[29]。而跌坑在徑流的持續掏蝕作用下，表現為溯源發展，跌坑間距逐漸減小，通過增加床面阻力來抑制徑流繼續侵蝕。

4 討 論

由于黃土坡面土壤流失情況日益加劇，使得水土流失治理成為當務之急，現行的研究多基于坡面細溝水流水動力學特性之上。所以本文嘗試從細溝縱橫斷面形態指標入手，研究細溝形態指標的變化過程，探究影響細溝形態指標的影響機制，以期為構建黃土坡面細溝侵蝕預報模型提供一定的理論支撐。結果表明，本試驗研究結果與野外試驗結果基本相同，但由于試驗條件的不同也存在一些差異。

已有野外試驗結果表明，描述細溝形態的基本參數細溝寬和細溝深均隨流量和坡度的增加而增加，而寬深比卻呈減小趨勢[30]。這一結果與本文結論一致。這表明坡面細溝在流量和坡度較小時以溝壁坍塌侵蝕為主，而在大流量和大坡度條件下以下切侵蝕為主。其原因是隨著流量的增大，徑流動能增大，隨著坡度的增大，流速增大，徑流剪切應力增大；在野外試驗中，細溝橫斷面隨流程長度的增加呈現出先增后減的變化趨勢[31]。是因為細溝侵蝕沿坡面呈現出強弱交替的波動趨勢。而降雨試驗研究表明，細溝橫斷面隨流程長度的增加而增加[6]。其原因主要與匯水面積的增加有關。這與本文研究結果存在差異，室內沖刷試驗條件下，徑流能量既要用于剝蝕床面，還要輸移泥沙顆粒，隨著流程長度的增加，徑流輸移的泥沙顆粒逐漸增多，對床面的剝蝕能力則減小。

由于目前試驗條件和方法所限，細溝發育過程中瞬時細溝形態變化還無法準確獲得，而野外動態監測細溝侵蝕過程的方法和技術還不成熟，缺乏對細溝形態特征參數量化效果的評價。因此今后關于細溝斷面形態的研究需要不斷改進試驗方法，來準確測量瞬時細溝斷面形態等指標，為進一步認識細溝侵蝕機理奠定基礎。

5 結 論

細溝縱橫斷面形態是水流沖刷的結果，水流沖刷后所呈現的細溝縱橫斷面形態也較復雜。徑流侵蝕造就了不同的斷面形態，而細溝斷面形態宏觀上對徑流侵蝕有一定的抑制作用。通過對縱橫斷面形態參數的分析來研究細溝形態變化規律，以及坡面形態與侵蝕的內在耦合機制。

1）橫斷面方向上，細溝寬深比隨坡度的增加而減小，且與梯形水力最佳斷面所對應的寬深比相差較遠，即床面還未達到穩定狀態。細溝橫斷面形態指標（細溝斷面面積與等深同寬的矩形面積之比）隨沖刷歷時和坡度的增加而增大，并逐漸趨近于梯形水力最佳斷面對應的橫斷面形態指標η。即斷面形態隨著沖刷歷時和坡度的增加向寬深式發育。

2）跌坑平均深度與縱向間距比值范圍為：0.60<<11.26。隨坡度的增大而增大，而與流量相關性不大。隨著坡度的增大，床面跌坑平均間距減小，跌坑深度增大，也隨之增大。

3）糙率系數和消能率均與呈良好的冪函數關系。細溝水力特性和細溝形態之間呈相互聯系、相互制約的作用與反作用系統。越大，即床面形狀阻力增大，跌坑也隨之發育愈成熟。水流在跌坑內強烈紊動并摻氣，對于消能有顯著的影響。

[1] 張樂濤，李占斌，王杉杉. 坡面徑流類型對侵蝕產沙及水沙傳遞關系的影響[J]. 農業工程學報，2015，31(20)：124－131.

Zhang Letao, Li Zhanbin, Wang Shanshan.Impect of runoff regimes on sediment yield and sediment flow behavior at slope scale[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(20): 124－131. (in Chinese with English abstract)

[2] Dong Y Q，Zhuang X H, Lei T W,et al.A method for measuring erosive flow velocity with simulated rill[J]. Geoderma, 2014, 232(11): 556－562.

[3] 吳普特，周佩華，武春龍，等. 坡面細溝侵蝕垂直分布特征研究[J]. 水土保持研究，1997，4(2)：47－56.

Wu Pute, Zhou Peihua, Wu Chunlong, et al. Research on the spatial distribution characters of slope rill erosion[J]. Research of Soil and Water Conservation, 1997, 4(2): 47－56. (in Chinese with English abstract)

[4] Bruno C, Stefano C D, Ferro V. Field investigation on rilling in the experimental Sparacia area, South Italy[J]. Earth Surface Processes and Landforms, 2008, 33(2): 263－279.

[5] 和繼軍，宮輝力，李小娟，等. 細溝形成對坡面產流產沙過程的影響[J]. 水科學進展，2014，25(1)：90－97.

He Jijun, Gong Huili, Li Xiaojuan, et al. Effects of rill development on runoff and yielding processes[J]. Advances in Water Science, 2014, 25(1): 90－97. (in Chinese with English abstract)

[6] 沈海鷗，鄭粉莉，溫磊磊，等. 降雨強度和坡度對細溝形態特征的綜合影響[J]. 農業機械學報，2015，46(7)：162－170.

Shen Haiou, Zheng Fenli, Wen Leilei, et al.The intensity of rainfall and slope comprehensive effects of rill morphological characteristics[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(7): 162－170. (in Chinese with English abstract)

[7] 唐克麗，鄭世清. 杏子河流域坡耕地的水土流失及其防治[J]. 水土保持通報，1984，4：5－8.

[8] 蔡強國. 坡面細溝發生臨界條件研究[J]. 泥沙研究，1998(1)：52－59

Cai Qiangguo. A study on critical condition about rill taken place on slope[J]. Journal of Sediment Research, 1998(1): 52－59. (in Chinese with English abstract)

[9] 嚴冬春，文安邦，史忠林，等. 紫色土坡面水流跌坑形態特征及其成因[J]. 山地學報, 2012, 30(1):30－35.

Yan Dongchun, Wen Anbang, Shi Zhonglin,et al. Purple soil water flow down pit morphological characteristics and its causes[J]. Journal of Mountain Research, 2001, 30(1): 30－35. (in Chinese with English abstract)

[10] 雷阿林，唐克麗. 坡溝系統土壤侵蝕研究回顧與展望[J].水土保持通報，1997，17(3)：37－43.

Lei Alin, Tang Keli. Retrospect and prospect of soil erosion research slope gully system[J]. Journal of soil and water conservation bulletin, 1997, 17(3): 37－43. (in Chinese with English abstract)

[11] 李占斌. 黃土地區坡溝系統暴雨侵蝕試驗及小流域產沙模型研究[D]. 西安：西安理工大學，1991.

Li Zhanbin. Heavy rain Erosion Slope Gully Loess Region System Test and Small Watershed Sediment Yield Model Research[D]. Xi’an: Xi’an University of Science and Technology, 1991. (in Chinese with English abstract)

[12] Church M, Zimmermann A. Form and stability of step‐pool channels: Research progress[J]. Water Resources Research, 2007, 43(43):10－1029.

[13] 徐江，王兆印. 階梯-深潭的形成及作用機理[J]. 水利學報，2004，35(10)：48－56.

Xu Jiang, Wang Zhaoyin. Ladder-deep pools formation and mechanism of action[J]. Journalofhydraulicengineering, 2004, 35(10): 48－56. (in Chinese with English abstract)

[14] 田嘉寧. 臺階式泄水建筑物水力特性試驗研究[D]. 西安：西安理工大學，2005.

Tian Jianing. Experimental Study Terraced Building Drainage Hydraulic Characteristics[D]. Xi’an: Xi’an University of Science and Technology, 2005. (in Chinese with English abstract)

[15] 肖培青，鄭粉莉，姚文藝. 坡溝系統坡面徑流流態及水力學參數特征研究[J]. 水科學進展，2009，20(2)：236－240.

Xiao Peiqing, Zheng Fenli, Yao Wenyi, et al. Slope gully slope surface runoff flow and hydraulics parameters characteristics research[J]. Advances in Water Science, 2009, 20(2): 236－240. (in Chinese with English abstract)

[16] 張攀，姚文藝，唐洪武，等. 模擬降雨條件下坡面細溝形態演變與量化方法[J]. 水科學進展，2015，26(1)：51－58.

Zhang Pan, Yao Wenyi, Tang Hongwu, et al. Evolution and quantization methods of rill morphology on the slope under rainfall simulation[J]. Advances in Water Science, 2015, 26(1): 51－58. (in Chinese with English abstract)

[17] 吳持恭. 水力學(上冊)[M]. 北京：高等教育出版社，1982：249－251.

[18] 沙玉清. 沖淤平衡穩定渠道設計方法[J]. 水利學報，1959(4)：25－44，85.

[19] 張科利，唐克麗. 黃土坡面細溝侵蝕能力的水動力學試驗研究[J]. 土壤學報，2000，37(1)：9－15.

Zhang Keli, Tang Keli. A study on hydraulic characteristic of flow for prediction of rill detachment capacity on loess slope[J]. Acta Pedologica Sinica, 2000, 37(1): 9－15. (in Chinese with English abstract)

[20] 張樂濤，高照良，李永紅，等. 模擬徑流條件下工程堆積體陡坡土壤侵蝕過程[J]. 農業工程學報，2013，29(8)：145－153.

Zhang Letao, Gao Zhaoliang, Li yonghong, et al. Soil erosion prosess of engineering accumulation in steep slope under simulated runoff conditions[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(8): 145－153. (in Chinese with English abstract)

[21] 丁文峰，李勉，張平倉，等. 坡溝系統侵蝕產沙特征模擬試驗研究[J]. 農業工程學報，2006，22(3)：10－14.

Ding Wenfeng, Li Mian, Zhang Pingcang, et al. Experimental study on the sediment yield characteristics in slope-gully system[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2006, 22(3): 10－14. (in Chinese with English abstract)

[22] 白清俊，馬樹升. 細溝侵蝕過程中水流跌坑的發生機理探討[J]. 水土保持學報，2001，15(6)：62－65.

Bai Qingjun, Ma Shusheng. Flow down to the pit in the process of rill erosion mechanism[J]. Journal of soil and water conservation. 2001, 15(6): 62－65. (in Chinese with English abstract)

[23] 沙際德，白清俊. 粘性土坡面細溝流的水力特性試驗研究[J].泥沙研究，2001(6)：39－43.

Sha Jide, Bai Qingjun. Stickness soil slope surface of rill flow hydraulic characteristics experimental study[J]. Journal of Sediment Research, 2001(6): 39－43. (in Chinese with English abstract)

[24] 徐江，王兆印. 山區河流階梯-深潭的發育及其穩定河床的作用[J]. 泥沙研究，2003(5)：21－26.

Xu Jiang, Wang Zhaoyin. Ladder-deep pools of mountainous rivers development and its effect on the stability of riverbed[J]. Journal of Sediment Research, 2003(5): 21－26. (in Chinese with English abstract)

[25] Curran Janet H, Wohl Ellen E. Large woody debris and flow resistance in step 2 pool channels[J]. Cascade Range, Washington, Geomorphology, 2003, 51(1/2/3): 141－157.

[26] 張霞，李鵬，李占斌，等. 黃土高原丘陵溝壑區臨界地貌侵蝕產沙特征[J]. 農業工程學報，2015，31(4)：129－136.

Zhang Xia, Li Peng, Li Zhanbin, et al. Characteristics of erosion and sediment yield under critical landform in hill-gully area of loess Plateau[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(4): 129－136. (in Chinese with English abstract)

[27] 王瑄，李占斌，李雯，等. 土壤剝蝕率與水流功率關系室內模擬實驗[J]. 農業工程學報，2006，22(2)：185－187.

Wang Xuan, Li Zhanbin, Li Wen, et al. Indoor simulation experiment of the relationshio between soil detachment rate and stream power[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2006, 22(2): 185－187. (in Chinese with English abstract)

[28] 李文哲，王兆印. 階梯-深潭系統消能機理的實驗研究[J].水利學報，2014，45(5)：537－546.

Li Wenzhe, Wang Zhaoyin, et al. Step-deep pools system energy dissipation mechanism of experimental research[J]. JournalofHydraulicEngineering, 2014, 45(5): 537－546. (in Chinese with English abstract)

[29] 李志威，王兆印，張晨笛，等. 人工階梯-深潭破壞案例與穩定性分析[J]. 水科學進展，2015，26（6）：820－828.

Li Zhiwei, Wang Zhaoyin, Zhang Chendi, et al. Artificial ladder-deep pools damage cases and stability analysis[J]. Advances in Water Science, 2015, 26(6): 820－828. (in Chinese with English abstract)

[30] 吳普特，周佩華，武春龍，等. 坡面細溝侵蝕垂直分布特 征研究[J]. 水土保持研究，1997，4(2)：47－56.

Wu Pute, Zhou Peihua, Wu Chunlong, et al. Research on the spatial distribution characteristics of slope rill erosion[J]. Research of Soil and Water Conservation, 1997, 4(2): 47－56. (in Chinese with English abstract)

[31] Di Stefano C, Ferro V, Pampalone V, et al. Field investigation of rill and ephemeral gully erosion in the Sparacia experimental area, South Italy [J]. Catena, 2013, 101(2): 226－234.

Influencing factor analysis of rill erosion section morphology development on slope and its dynamic characteristic experiment

Ma Xiaoling1, Zhang Kuandi1,2※, Yang Fan1, Yang Mingyi2, Dong Xu1

(1.,,712100,; 2.,,712100,)

Rill erosion is the process of soil particles separating in the groove and handling by rill flow, and it is one type of surface corrosion, and the rill formed can be completely filled by normal farming activities. In the loess plateau region, rill erosion of slope accounts for more than 70% of the total erosion, and rill erosion is a major erosion in the loess plateau hilly land, and is also an important part of mechanism of slope erosion prediction. Rills are small channels which function as both sediment sources and sediment transport vehicles on hill slopes. It is a quantum leap during the process of soil erosion on slopes, and is also the beginning of qualitative changes in the process of soil erosion. The rill morphology development plays an important role in the understanding of rill erosion.Rill morphology changes have many influence factors in the process of rill erosion, so systematically studying the relationship between rill shape coefficient and influencing factors can make us better understand the process of rill erosion and erosion mechanism. For the study of rill morphological development process, the combination scouring tests of 6 slopes (2°, 4°, 6°, 8°, 10°, and 12°) and 5 discharges (8, 16, 24, 32, and 40 L/min) were carried out, and the effect mechanism of rill transverse profile morphological development on the loess slope surface was investigated, as well as the dynamic characteristics of the impact. The results showed that rill width-depth ratio ranged from 3.006 to 4.841, and according to the hydraulic optimum section, the rill flow was far from steady state. In rill erosion, downward erosion and lateral erosion gradually transited to downward erosion. Because of incised and lateral erosion, rill cross-sectional shape was gradually widened and deepened with the increase of slope and time.Along with the change of slope, cross-sectional shape coefficient varied within the range of 0.36-0.522.Research showed that: The rill cross-sectional shape was approximating to the trapezoidal hydraulic best section with the increase of the flow, slope and scour lasted time, namely the section with the least resistance.With the increase of water flow length, cross section shape changed gradually from wide and deep to narrow, and the cross-sectional shape coefficient was reduced. Through calculation, the groove profile form ranged within 0.6-11.26, and increased with the increase of slope, but had no significant correlation with discharge.With the increase of slope, the average spacing of drop pit on bed face decreased, the depth of drop pit increased, and the groove profile shape coefficient increased. Comprehensive resistance coefficient and energy dissipation showed a good power function relationship with groove profile shape coefficient, and with the increase of profile shape coefficient, the slope shape resistance increased. Rill hydraulic characteristics and groove shape showed the role of mutual connection and restriction and reaction system.Groove profile shape coefficient was larger, the bed surface shape resistance increased, and drop pit also would be developed more maturely.The water in drop pit had strong turbulence and aeration, which had significant effect on energy dissipation.

erosion; experiment; soils; rill morphology; form index coefficient; drop pit; energy dissipation efficiency

10.11975/j.issn.1002-6819.2017.04.029

S279.2

A

1002-6819(2017)-04-0209-08

2016-06-22

2017-01-03

國家自然科學基金資助項目“坡面含沙水流波流耦合特性及泥沙輸移機理研究”（51579214），中央高校基本科研業務費專項資金“黃土坡面含沙水流能量耗散規律研究”（Z109021564）；黃土高原土壤侵蝕與旱地農業國家重點實驗室主任基金“基于超聲測量技術研究坡面流流動特性”（A314021402-1619）；流域水循環模擬與調控國家重點實驗室開放基金“坡面滾波演化及挾沙能力研究”（IWHR-SKL-2014）

馬小玲，女，甘肅白銀人，博士生，主要從事水土保持及土壤侵蝕方面的研究。楊凌 西北農林科技大學水利與建筑工程學院，712100。Email：15002981284@163.com

張寬地，男，寧夏隆德人，副教授，工學博士，主要從事水工水力學及坡面水流研究。楊凌 西北農林科技大學水利與建筑工程學院，712100。Email：zhangkuandi428@126.com

馬小玲，張寬地，楊 帆，楊明義，董 旭.坡面細溝侵蝕斷面形態發育影響因素分析及動力特性試驗[J].農業工程學報，2017，33(4)：209－216. doi：10.11975/j.issn.1002-6819.2017.04.029 http://www.tcsae.org

Ma Xiaoling, Zhang Kuandi, Yang Fan, Yang Mingyi, Dong Xu. Influencing factor analysis of rill erosion section morphology development on slope and its dynamic characteristic experiment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(4): 209－216. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.04.029 http://www.tcsae.org