不同長度小麥秸稈覆蓋下黃土耕地坡面流水動力學特性

李朝棟，李占斌，馬建業，馬 波

?

不同長度小麥秸稈覆蓋下黃土耕地坡面流水動力學特性

李朝棟1,2，李占斌1，馬建業1，馬 波1※

（1. 西北農林科技大學水土保持研究所黃土高原土壤侵蝕與旱地農業國家重點實驗室，楊凌 712100；2. 西北農林科技大學資源環境學院，楊凌 712100）

為了系統研究秸稈長度對坡面流水動力學特性的影響，利用室內人工模擬沖刷定床阻力試驗，在4個坡度、5個流量條件下分別對3～5、8～10和13～15 cm長度秸稈覆蓋坡面流水動力學參數變化特征進行分析，闡明不同長度秸稈覆蓋坡面水流阻力與雷諾數的關系。結果表明，秸稈長度對坡面流水動力學參數影響顯著。不同長度秸稈覆蓋條件下，坡面流雷諾數變化范圍為166～558，當流量≤7.5 L/min時，坡面流態為層緩流；流量為9.0 L/min時，坡面流態為過渡緩流。秸稈覆蓋條件下，坡面流具有較小的流速和較大的水深、Darcy-Weisbach阻力系數及曼寧糙率系數。隨著秸稈長度增加，坡面流流速隨之增加，而水深、Darcy-Weisbach阻力系數及曼寧糙率系數隨之降低。隨著秸稈長度和水力坡度增加，流態指數值總體呈現降低趨勢。當秸稈長度由3～5 cm增加至13～15 cm時，流態指數平均由0.716降至0.501。研究可為秸稈覆蓋條件下水土流失阻控機理研究提供科學依據。

流型；水動力學；秸稈；坡面流；定床試驗；黃土耕地

0 引 言

在活動頻繁的坡耕地上，由坡面徑流造成的土壤侵蝕已成為一個嚴重的問題。中國坡耕地水土流失量占該區域總流失量的比例大部分均高于60%，坡耕地土壤侵蝕已成為中國生態環境建設和可持續發展的首要問題之一[1]。秸稈覆蓋是中國重要的保護性耕作措施之一，對減少坡面土壤侵蝕具有積極的作用。目前針對秸稈覆蓋的研究主要集中于秸稈還田的土壤水肥效應和防蝕效益等方面，對其侵蝕動力機制研究較少[2-8]。秸稈覆蓋不僅有效攔截消耗降雨動能、降低擊濺侵蝕、減少土壤顆粒分離與擴散，還可有效增加地表糙率和水分入滲，從而對減少徑流、降低流速、降低對土壤的沖刷作用產生積極作用；與此同時，秸稈在降雨過程中形成低洼蓄水區促使泥沙沉積，可有效減少土壤侵蝕[9]。秸稈覆蓋作為保護性耕作的一種典型措施，當其地表秸稈覆蓋為30%時便可達到減少農田土壤風蝕或水蝕等目的[10]。秸稈覆蓋能顯著減少地表徑流73.9%～86.2%，但增加壤中流15.4%～156.4%，使徑流總量降低32.5%～66.6%，并極顯著降低土壤侵蝕總量96.4%～98.1%[11]。楊青森等[2]對黑土區的秸稈覆蓋研究認為，秸稈覆蓋可有效減少坡面土壤侵蝕，秸稈覆蓋量為2 kg/m2時，坡面徑流量減少87%，侵蝕量相應減少99.86%；而秸稈覆蓋量為4 kg/m2時，坡面幾乎不產生侵蝕。但秸稈覆蓋對土壤侵蝕的作用受秸稈量影響較大，研究認為，當秸稈覆蓋率>40%時能有效地控制水土流失，但當秸稈覆蓋率<40%時作用不明顯[9]。

秸稈覆蓋不僅對降雨有截留持水作用，而且對降雨侵蝕能量的削減和在地表的分布產生復雜的影響，進而有可能對坡面侵蝕及細溝發育產生影響。裸露坡面細溝流侵蝕力與其水力學參數關系密切，不同的徑流系數可以用來衡量細溝流侵蝕力；但是當坡面存在作物秸稈時，在裸露地表上得出的各種物理參數關系均無法得到較好應用[10,12]。坡面流是由降雨形成并沿坡面流動的淺層水流，相比較明水流，坡面流底坡陡，水深淺，受坡度、雨強和地表狀況等多種因素影響，這使得坡面流的研究較為困難[13]。目前，對坡面流水動力學特征的研究主要集中于裸露坡面上，針對植被、枯落物和秸稈覆蓋等地表覆蓋措施的坡面水動力學特征的研究較為少見。不同的覆蓋，在降低坡面土壤流失量的同時，還可顯著改變地表徑流水力學特性，例如，植被殘留物、秸稈和碎石等覆蓋對降低流速均有顯著效果[10]。由于秸稈覆蓋和土壤緊密接觸，尤其是淺耕半固定混合式覆蓋，強烈地改變了土壤表面粗糙度，因此對地表水力學參數影響很大。而在細溝侵蝕方面作物秸稈對細溝水流水力學特性和侵蝕力影響程度的計算方法與當前對細溝流水力學特性及其和細溝形態間關系的研究結果并不相符[10]。Foster等[14]曾指出，在土中混入秸稈可有效防止細溝向下深切，且表土中秸稈與土壤結合越緊密對細溝的影響更甚，當混入土中的秸稈由于徑流沖刷暴露在外，可降低土壤的剪應力。由此可見，秸稈覆蓋條件下侵蝕動力及水動力學特性與裸坡條件下存在諸多差異，難點和疑問也較多。本文通過室內秸稈覆蓋人工模擬徑流沖刷定床阻力試驗，對不同長度秸稈覆蓋下的坡面流水動力學特性進行系統研究，以期闡明秸稈覆蓋條件下坡面流水動力學參數的變化特征，為秸稈覆蓋條件下水土流失阻控機理研究提供科學依據。

1 材料與方法

1.1 試驗方法

本研究在西北農林科技大學水土保持工程實驗室進行(34°14′～34°20′N、107°59′～108°08′E)，試驗采用人工模擬徑流沖刷法。供試秸稈為當年自然曬干的小麥秸稈，于2015年5月底人工收割而得，再根據當地秸稈還田典型長度并適當分級，利用鍘刀裁剪成3～5、8～10和13～15 cm共3種試驗所需的秸稈長度，坡面秸稈覆蓋量為3.5 t/hm2。本試驗采用定床阻力試驗，如此既可以模擬天然地表糙度，又可以消除床面形態變化、底沙交換對水流紊動的影響[15]。人工模擬徑流沖刷定床阻力試驗坡度依據黃土區坡耕地典型坡度設為3°、5°、10°和15°；試驗設計放水流量依據黃土高原地區侵蝕性降雨范圍并考慮秸稈抵抗流水沖刷產生漂移的承受能力，最終將放水流量設計為3.0、4.5、6.0、7.5和9.0 L/min。3種長度秸稈覆蓋坡面及裸坡在不同坡度和流量下分別進行試驗，重復測量3次，共計240個處理。

試驗裝置（圖1）由供水水泵、流量計、穩流槽和沖刷槽等構成，其中沖刷槽采用變坡式鋼板水槽，結構尺寸長2 m、寬0.3 m、深0.3 m，其中填土深度為0.25 m，坡度可在0°～15°范圍內任意調節。由于本試驗采用定床阻力試驗，故需將沖刷水槽進行填充壓實并固化床面，采用當地的土壤（旱耕土墊人為土），將土壤自然風干后過5 mm篩，填土設計容重約為1.3 g/cm3，可以代表黃土高原坡耕地土壤自然狀態。填土過程采取分層填土、分層壓實的方法，每層填土深度為5 cm，共計填土5層。將試驗土槽調至水平。為了使土槽內土地自然沉降，便于下一步的固化床面。在20 mm/h降雨強度對土槽進行前期降雨約1 h，使水分充分入滲，將其土壤自然沉降壓實，也可消除人為裝填土壤的松散狀況，此后再將沉降的部分填至水平。將處理好的土槽進行自然風干，然后利用噴漆槍將清漆均勻地噴灑并完全覆蓋在土槽表面以固化床面。將直徑為2～5 mm的土壤顆粒均勻地灑在噴灑過清漆的土槽表面，再次噴灑清漆，使其包裹住表面的土壤顆粒。將土壤顆粒粘于坡面土壤表面，可在床面形成穩定的糙率，且在試驗中水流無法改變坡面糙率；這樣模擬天然土壤地表糙率，不僅可消除下墊面糙率變化對水流阻力等的影響[16]，又可防止床面過于光滑導致覆蓋在上面的秸稈沿水流方向移動。

試驗開始前調整坡度和放水流量，流量大小由流量計和閥門控制，用稱質量法對所需流量進行率定。將剪裁好的秸稈在水中浸泡12 h，使其充分吸水，以防止其隨水流漂浮移動，同時還防止其吸水影響水槽出口流量。浸泡好的秸稈將其瀝干。將沖刷槽調整至水平，在距離床面30 cm高度處，雙手將秸稈揚起使其自然墜落于床面，覆蓋層內保持秸稈交錯縱橫的隨機位置。鋪撒完畢后人工調整床面不同位置的秸稈，使其能夠均勻分布于沖刷槽表面，并形成均勻厚度的秸稈覆蓋。3～5 cm長的秸稈坡面厚度約為1.5 cm，13～15 cm長的秸稈由于其秸稈較長，秸稈層空隙較大，其厚度約為2.5 cm。試驗按水流方向設3個斷面，斷面間距為0.6 m，距沖刷槽頂端0.1 m處為第1觀測斷面，依次為第2和第3觀測斷面，每個觀測斷面按橫向方向分左、中、右采用常規高錳酸鉀染色法測得流速。試驗的同時將水銀溫度計（精度0.1 ℃）放置于穩流槽內測定水流溫度，以計算雷諾數。

1.2 水動力學參數計算方法

利用沖刷槽3個斷面多次觀測的表面流速，求取平均值后乘以修正系數0.75，得到水流斷面平均流速[17]。文中涉及的水深（）、雷諾數（）、佛汝德數（）、運動黏滯力系數（）、Darcy-Weisbach阻力系數（）、曼寧糙率系數（）的計算公式如下

=/（1）

=/（2）

=0.017 75/(1+0.033 7+0.000 2212) （3）

=/()0.5（4）

=8/2（5）

=-12/31/2（6）

式中為水深，m；為單寬流量，L/(m·s)；為斷面流速，m/s；為水力半徑，m，其取值可近似等于水深，m；為運動黏滯系數，m2/s；為水溫，℃；為Darcy-Weisbach阻力系數；為重力加速度，取9.81 m/s2；為水力坡度，其值為坡度的正弦值，即=sin，為坡度，(°)；為曼寧糙率系數；為水流容重，10 kN/m3。

2 結果與分析

2.1 秸稈覆蓋坡面流流態及流型

水流流態是坡面薄層徑流水力學特征的基本表征參數之一，判別方法一般采用雷諾數；雷諾數是衡量水流紊動程度的重要指標，其值為水流慣性力與黏滯力之比[18-19]。根據明渠均勻流的基本理論，當雷諾數小于500時，水流屬于層流；當雷諾數大于2 000時，水流屬于紊流；介于二者之間為過渡流。當水流呈層流時，黏滯力的作用大于慣性力，黏滯力對液流質點的運動起抑制作用；反之，水流呈紊流時，慣性力作用大于黏滯力，慣性力對液流質點的運動起推動作用[19]。一般而言，在流量恒定的條件下，雷諾數越大，坡面流的紊動性越大，水流剝蝕以及挾沙能力也越強[18]。此外，水流流態也可用佛汝德數進行表征，它反映了流速和徑流深的對比關系。在定流量條件下，越大，表明徑流挾沙能力越強、坡面徑流剪切力越大。依據明渠流的判別標準，佛汝德數≥1，坡面流態為急流，<1則為緩流[17]。本研究特定試驗條件下，通過雷諾數和佛汝德數來判別水流流態與流型。不同長度秸稈覆蓋條件下坡面流流態流型相關參數如表1所示。

表1 不同流量和秸稈長度條件下水流流態流型相關水力參數

注：均值為4個坡度下的平均值。采用鄧肯檢驗，不同字母表示處理間差異顯著(<0.05)。

Note: Values are mean of 4 slope gradients. Different letter means the significant difference among treatments by Duncan’s test (<0.05).

由表1可知，在試驗設定的流量和坡度范圍內，秸稈覆蓋坡面流雷諾數變化介于166~558之間，且在試驗流量和坡度范圍內變化可分為2個階段。不同長度秸稈覆蓋及裸坡坡面流在流量≤7.5 L/min時，雷諾數均小于500，其流態均為層流；當流量為9.0 L/min時，其流態則為過渡流。秸稈覆蓋對雷諾數的影響較小，進而對流態的影響也相應較小。與裸坡相比，當流量為3.0 L/min時，秸稈覆蓋坡面流雷諾數較裸坡有較大提升，3～5、8～10和13～15 cm長度條件下，雷諾數分別較裸坡增加了27.41%、24.16%和32.11%（<0.000 1）。而當流量≥4.5 L/min時，同一流量條件下，秸稈覆蓋對雷諾數主要表現為降低效應，但是降低幅度較小；最大降幅為6.61%（<0.05，流量為7.5 L/min，秸稈長度為8～10 cm）。可以看到，當流量為3.0 L/min時，秸稈覆蓋條件下坡面流雷諾數較裸坡增加了24.16%～32.11%（<0.000 1）。在小流量情況下（3.0 L/min）使水深提升后基于最小能耗原理穿行于秸稈之間，增加了水流與粗糙單元接觸和碰撞的概率；加之水深無法沒過秸稈層，在水流表面與秸稈碰撞已形成壅浪和尾渦，進而增加了水流的紊動性[20]，使之在相同流量條件下較裸坡具有較高的雷諾數。但隨著流量增大（≥4.5 L/min），水深進一步增加，使大部分秸稈淹沒于水流之下，減少了表面水流接觸秸稈后增加紊動性的可能性；水體內部流速垂線分布也會趨于均勻，有利于降低雷諾數。相同秸稈長度下，雷諾數隨流量增大呈現極顯著增大（<0.000 1）。

由流型的判定標準可知，在本研究試驗范圍內，裸露條件下，坡面流流型均為急流；而秸稈覆蓋條件下，坡面流流型均表現為緩流。由于秸稈覆蓋對水流形成較大的阻礙，流速減慢且水深加深，呈現出層緩流的狀態（其中流量為9.0 L/min時為過渡緩流）。

2.2 秸稈覆蓋坡面流水深及平均流速

水深是坡面徑流水動力學最基本的參數之一，其變化常受到微地形、覆被狀況、降雨和其他外界擾動因素的影響，而本研究采用固定床面的方法，使坡面糙率保持穩定，水流中不含泥沙，不存在坡面微地形變化等問題，因而可準確測定坡面流水深。在本研究特定試驗條件下，水深與單寬流量和坡度的關系如圖2所示。

秸稈覆蓋條件下，坡面流水深較裸坡均有所增加（圖2）。在當前試驗坡度和流量范圍內，秸稈覆蓋較裸坡水深增加了107.69%～178.64%，平均為145.64%（<0.000 1）。秸稈長度對水深的影響較為顯著，隨秸稈長度增加，對水深的影響逐漸減弱，水深隨之減小。相同坡度下，水深與流量呈正相關關系；而相同流量下，隨坡度增加，平均水深呈減小趨勢。同一秸稈長度覆蓋條件下，隨著流量增大，水深呈增加趨勢；隨坡度增大，水深呈減小趨勢（圖2）。秸稈覆蓋條件下，隨流量增加，不同坡度間水深差異逐漸擴大；隨坡度增加，不同流量間水深差異較小；這種變化趨勢在3~5 cm秸稈長度覆蓋坡面表現最為顯著，即隨著秸稈長度增加，不同流量下的水深差異逐漸降低。

不同長度秸稈覆蓋條件下斷面流速與單寬流量和坡度的關系如圖3所示。由圖可知，無論覆蓋秸稈與否，流速隨流量和坡度的增加均呈現增加趨勢。秸稈覆蓋對流速具有顯著的減緩作用，與裸坡相比，流速平均降低55.54%，其中3~5、8~10和13~15 cm 3種長度秸稈覆蓋坡面流流速分別較裸坡降低了60.64%、57.67%和48.30%（<0.05）。裸坡和秸稈覆蓋坡面間流速差異隨流量和坡度增大逐漸增大。

江忠善等[21]對國內外諸多學者的坡面流試驗資料進行總結分析，認為坡面流速是單寬流量與水力坡度的冪函數形式，即=··J，式中、、均為系數，經回歸分析求得為0.5，為0.35。本研究中裸坡和不同長度秸稈覆蓋坡面平均流速與單寬流量和水力坡度的冪函數關系如下所示：

=0.0220.363·-0.3042=0.72（裸坡） （7）

=0.165sm0.209·0.348·0.2952=0.86（秸稈覆蓋）（8）

式中sm為平均秸稈長度，cm，取不同秸稈長度范圍之平均值分別為4、9和14 cm。上式中流量和坡度的指數均略小于江忠善的0.5和0.35[21]，同時也略小于張光輝等[16-17]的結果。

由于秸稈的存在，加之水深較淺，坡面流各點的流速并非一致，而是基于最小能耗原理繞過秸稈并流動于秸稈之間，因此在局部區域紊動性可能較強。Horton假定坡面流是層流和紊流流態相互交錯的混合流，即在完全的湍流中間點綴有層流[22]。基于此，張寬地等[20]基于流體力學理論推導出平均流速與單寬流量和水力坡度的關系=·1-m·J，式中、和均為系數。其中斷面平均流速與單寬流量的（1-）次方呈正比，與水力坡度的次方呈正比；張寬地等將定義為流態指數，表征了單寬流量對流速的影響程度；認為當值較大時，流體中能量的消耗主要以阻力做功為主，當值較小時，流體主要通過增加流體流速進行能量轉化，阻力做功較弱[20]。不同長度秸稈覆蓋條件下，各水力坡度下擬合冪函數關系中流態指數值如表2所示。

表2 不同水力坡度下流態指數值

由表2可知，秸稈覆蓋坡面流態指數值均大于同等條件下的裸坡。隨著秸稈長度增加到15 cm和水力坡度增加到0.259，值總體呈現降低趨勢。當秸稈長度由3～5 cm增加至13～15 cm時，流態指數平均由0.716降至0.501。當秸稈長度較短時，單位面積上的秸稈個體小而更多，內部交叉縱橫結構更為復雜，空隙更小，對坡面覆蓋的效果也更好，從而使坡面阻力增強。濕周增大，水流與秸稈個體接觸的面積和頻率增加，水流各流層之間的水體交換愈加頻繁，流速在垂線上的分布更加趨于均勻[20]。在這種情況下，阻力做功耗能處于上峰，而流速增加轉化能量處于次要因素，水流紊動性減弱。

2.3 秸稈覆蓋坡面流阻力

徑流沿程運動必然會受到床面起伏變化和覆被物的阻滯作用，水流所受阻力強弱可影響徑流的流速，進而影響徑流的侵蝕力[23]。由于床面為穩定糙率設計，因此在水流阻力分析中可忽略床面沙粒、徑流泥沙以及床面微地形改變等因素的影響，但仍需要考慮秸稈覆蓋對水流的阻礙作用。本研究采用Darcy-Weisbach阻力系數和曼寧糙率系數來分析秸稈覆蓋對坡面流阻力的影響。秸稈覆蓋對Darcy-Weisbach阻力系數的影響見圖4。

由圖4可知，在試驗設定的流量和坡度范圍內，裸坡阻力系數的變化范圍為0.14～1.47，而在3～5、8～10和13～15 cm 3種長度秸稈覆蓋坡面上的變化范圍分別為2.67～8.00、2.39～6.36和1.31～4.28，平均分別較裸坡增加了20.89、15.69和8.17倍（<0.000 1），由此說明秸稈覆蓋對坡面流具有較強的阻礙作用。不同長度的秸稈，對水流的阻力作用差異較大，其中3～5 cm長度秸稈對水流的阻力最強，13～15 cm則最弱。

丁文峰等[23-24]的研究認為，Darcy-Weisbach阻力系數與雷諾數之間也存在=1·-p2的冪函數關系，即阻力系數與雷諾數呈反比例關系。與其研究相似，秸稈覆蓋條件下，阻力系數與雷諾數的關系均可用冪函數來描述。其中裸坡和13～15 cm長度秸稈覆蓋條件下，冪函數指數為負值，即阻力系數與雷諾數成反比例關系；3～5和8～10 cm長度秸稈覆蓋條件下，指數為正值，即阻力系數與雷諾數呈正比例關系。說明在秸稈長度較小的情況下，阻力系數隨雷諾數增大而呈增加趨勢。

秸稈覆蓋使坡面流阻力系數表現出較大的不確定性，這主要是由于徑流沖刷使秸稈在坡面上重新分布和堆疊，造成局部分布不均，從而增加了阻力系數的不確定性。但是總體而言，秸稈覆蓋對增加坡面流阻力系數是有積極作用的。由方差分析可知，秸稈長度、流量和坡度對阻力系數的影響均達到了極顯著水平（<0.000 1）。

曼寧糙率系數表征地表粗糙程度對水流流動過程中的阻滯作用。由圖5所示，在試驗設定的流量和坡度范圍內，裸坡曼寧糙率系數的變化范圍為0.013～0.044，而在3～5、8～10和13～15 cm 3種長度秸稈覆蓋坡面上的變化范圍分別為0.071～0.132、0.067～0.118和0.049～0.083，平均分別較裸坡增加了4.38、3.64和2.33倍（<0.000 1）。

秸稈長度不同，對曼寧糙率系數的影響也有所不同。隨著秸稈長度增加，對水流的阻滯作用逐漸減弱，曼寧糙率系數隨之降低；其中3～5 cm長度處理下曼寧糙率系數是8～10 cm長度處理的1.15倍（<0.000 1），是13～15 cm長度處理的1.61倍（<0.000 1），這說明當秸稈較短時，對坡面水流的阻滯作用更加強烈。曼寧糙率系數與雷諾數的關系與阻力系數較為相似。裸坡坡面上，糙率系數隨雷諾數增加而減小。秸稈覆蓋坡面上，曼寧糙率系數隨雷諾數增加而增加。由于短秸稈相對于長秸稈具有更多的秸稈個體數量，個體之間空隙較小、密度較大，在減少過水斷面面積的同時又增加了水流的濕周，導致阻力變化較大；而長秸稈由于單位面積上個體數量較少，個體間空隙較大，對水流過水面的減少不及短秸稈，表現出較低的阻滯作用，這與王俊杰等在不同覆蓋度下得出的結論較為相似[25]。由方差分析可知，秸稈長度、流量和坡度對曼寧糙率系數的影響均達到了極顯著水平（<0.000 1）。

3 討 論

秸稈覆蓋是重要的水土保持措施，其不僅能攔截消耗降雨動能、降低濺蝕，還能有效較少坡面產流產沙。本研究通過對不同秸稈長度條件下坡面流水動力學特性的研究，通過水動力學參數的變化特征來闡述和解釋秸稈覆蓋對土壤侵蝕阻控動力機制問題。秸稈覆蓋是重要的水土保持措施，不僅能攔截降雨動能、降低濺蝕，又可有效減少坡面產流產沙。作為地表覆被的一種類型，其水動力學特性與植被覆蓋措施具有諸多相似之處。當地表覆蓋有秸稈時，水流內部原有結構產生較大變化，增加了水流阻力，不僅能有效攔截阻滯徑流和泥沙，還影響了坡面的產匯流過程，這與植被覆蓋條件下是相似的[20]。

目前，有關坡面流水動力學特性的研究主要是在裸露坡面上開展的，針對地表覆蓋坡面水動力學特性的研究還相對較少。由于地表覆蓋差異較大，開展此類研究亦相對較為困難，目前主要集中于植被覆蓋坡面的水動力學特性研究，且分歧較大。其中就流態而言，一般認為植被覆蓋條件下，坡面流流態為層緩流，隨流量增加逐漸轉變為過渡緩流，并不排除隨流量進一步增大轉變為紊流的可能[20,26-30]。本研究雖以秸稈覆蓋為主，但是不同長度秸稈覆蓋條件下，坡面流流態變化規律與植被覆蓋坡面較為一致，均為隨流量增加，坡面流由層緩流逐漸轉變為過渡緩流。張寬地等[20]研究認為，植被覆蓋條件下，當雷諾數<500時，水流流態不同于傳統意義上的層流。由于植被阻擋，水流基于最小能耗原理繞行于植被之間，并在迎水面行程壅浪，在背水面產生尾渦，故將植被覆蓋條件下的層流定義為“虛擬層流”。而在本研究中，由于秸稈相對植被覆蓋度更高，秸稈內部又交錯縱橫穿插，對水流的穿梭繞行影響更大。因此，在試驗過程中不僅可以觀測到水流穿行于秸稈個體之間，還能觀測到表面水流越過秸稈上方流動，同樣產生了壅浪和尾渦。這與張寬地等在植被覆蓋條件下所觀察到的現象相似而又有所不同。因此，秸稈覆蓋條件下，不論秸稈長度如何變化，其層流也可認為是區別于傳統層流意義上的“虛擬層流”。

坡面流是坡面土壤侵蝕發生發展的重要原因，而坡面流的侵蝕力主要由流速決定。作為水土保持重要措施之一的秸稈覆蓋，能夠有效減緩坡面流流速，這對防止坡面水土流失具有重要的意義。而秸稈覆蓋坡面流流速減緩主要體現在坡面阻力的增加上。一般認為，坡面流阻力由顆粒阻力、波阻力、形態阻力和降雨阻力4部分構成。由于本研究采用定床阻力試驗方法，故可忽略降雨阻力的影響。由于波阻力是坡面流受地表障礙物阻擋而產生的阻力，因此在秸稈覆蓋條件下，坡面流阻力主要受波阻力影響最大，但由于坡面水深較淺，還需要考慮顆粒阻力產生的影響。秸稈覆蓋條件下，坡面流阻力顯著降低，這與肖培青等[31]、王俊杰等[25]和張冠華[32]在植被覆蓋條件下的研究結果一致。其中肖培青等[31]研究發現草地和灌木坡面阻力系數分別是裸坡的5.58～7.45倍和5.61～6.26倍；而本研究中3種長度秸稈覆蓋坡面阻力系數分別較裸坡增加了20.89、15.69和8.17倍，充分顯示了秸稈覆蓋的“增阻”效應。

植物殘體覆蓋土表是當前普遍采用的一種保土措施，在中國的坡耕地上使用也較為廣泛；中國干旱地區，播種后往往會迎來當年的雨季，此時地表植被覆蓋稀疏，秸稈覆蓋將對水土保持起到重要的作用[33]，而對秸稈覆蓋條件下坡面水動力學特性的研究，對在坡耕地上合理布設水土保持措施具有重要的意義。在實際生產建設活動中，由于人為破壞了原生地質地貌、生態水文、土壤植被等，為土壤侵蝕的發生創造了條件[34]；在人工堆積體土石混合介質表面布設秸稈覆蓋措施可有效防止土壤流失的產生。秸稈是否與表土緊密接觸是秸稈發揮防蝕作用的重要條件，加之工程堆積體土體較為松散、坡度相對較陡，因此秸稈越短越有利于在表土層形成致密的秸稈層，對坡面徑流的阻力更大，更有利于防止水土流失。秸稈覆蓋影響坡面土壤侵蝕的研究目前較多，但多側重于秸稈覆蓋量和秸稈覆蓋度等因素，而針對秸稈長度的研究還較為少見。由于秸稈長度不同，決定了單位質量秸稈個體數量的多寡；秸稈越小，單位個體越多；反之則越少。與此同時，短秸稈更利于在地表行程交叉縱橫、相互疊加的空間結構，減緩流速，增加了坡面流阻力；但缺點是當流量較大時更容易隨水流發生空間位置的改變和沿水流方向產生短距離遷移。例如橫向位置的秸稈在水流作用下產生轉向，從而對水流產生導流作用，進而使局部流速加快。而長秸稈條件下，水流阻力相對較弱，水流繞流而行，對秸稈的拖拽力也相對較小，反而不易發生空間位置的改變。由此也說明秸稈長度變化對坡面流水動力學參數的影響較為深刻，進而也會對坡面土壤侵蝕過程產生一定的影響，需要引起重視。此外，每次人工模擬沖刷試驗所用秸稈并不是固定在沖刷槽床面上的，而是每次試驗時更換新的秸稈。因此，每次試驗的坡面秸稈空間結構均會有所不同，可能會對試驗結果產生一定的影響。后續的研究還需考慮將秸稈均勻鋪撒并黏貼于沖刷槽床面，以便更精確地分析秸稈覆蓋下坡面流水動力學參數的變化特征。

4 結 論

本文利用室內人工模擬徑流沖刷手段，采用定床阻力試驗方法，系統研究不同長度的秸稈覆蓋條件下坡面流水動力學參數隨流量和坡度的變化規律。結果表明，秸稈覆蓋條件下，坡面流以層緩流為主，并隨流量增加逐漸轉變為過渡緩流。不同長度的秸稈覆蓋坡面流雷諾數較為接近。秸稈覆蓋坡面流水深較裸坡增加145.64%，并隨秸稈長度增加而降低；流速較裸坡降低了55.54%，隨秸稈長度增加呈增加趨勢。秸稈覆蓋條件下，坡面流流態指數值隨秸稈長度和坡度的增加總體呈降低趨勢；當秸稈長度由3～5 cm增加至13～15 cm時，流態指數平均由0.716降至0.501。坡面流Darcy-Weisbach阻力系數隨秸稈長度增加而降低，并隨秸稈長度增加，阻力系數與雷諾數的冪函數關系由正比例關系轉變為反比例關系。曼寧糙率系數隨秸稈長度增加而顯著降低，且3種秸稈長度條件下，曼寧糙率系數隨雷諾數呈增加趨勢。

[1] 唐克麗. 中國土壤侵蝕與水土保持學的特點及展望[J]. 水土保持研究，1999，6(2)：3－8.

Tang Keli. Characteristics and perspectives on scientific discipline of soil erosion and soil and water conservation in China[J]. Research of Soil and Water Conservation, 1999, 6(2): 3－8. (in Chinese with English abstract)

[2] 楊青森，鄭粉莉，溫磊磊，等. 秸稈覆蓋對東北黑土區土壤侵蝕及養分流失的影響[J]. 水土保持通報，2011，31(2)：1－5.

Yang Qingsen, Zheng Fenli, Wen Leilei, et al. Effect of mulch cover on hillslope soil erosion and nutrient loss in black soil region of Northeast China[J]. Bulletin of Soil Water Conservation, 2011, 31(2): 1－5. (in Chinese with English abstract)

[3] Bhatt R, Khera K L. Effect of tillage and mode of straw mulch application on soil erosion in the submontaneous tract of Punjab, India[J]. Soil and Tillage Research, 2006, 88 (1/2): 107－115.

[4] Jordán A, Zavala L M, Gil J. Effects of mulching on soil physical properties and runoff under semi-arid conditions in southern Spain[J]. Catena, 2010, 81(1): 77－85.

[5] Edwards L, Burney J R, Richter G, et al. Evaluation of compost and straw mulching on soil-loss characteristics in erosion plots of potatoes in Prince Edward Island, Canada[J]. Agriculture, Ecosystems & Environment, 2000, 81(3): 217－222.

[6] 嚴坤，王玉寬，徐佩，等. 秸稈覆蓋對三峽庫區坡面侵蝕的影響[J]. 水土保持通報，2016，36(1)：6－10.

Yan Kun, Wang Yukuan, Xu Pei, et al. Effects of mulch cover on slope erosion in Three Gorges Reservoir area[J]. Bulletin of Soil and Water Conservation, 2016, 36(1): 6－10. (in Chinese with English abstract)

[7] 徐錫蒙，鄭粉莉，吳紅艷，等. 玉米秸稈覆蓋緩沖帶對細溝侵蝕及其水動力學特征的影響[J]. 農業工程學報，2015，31(24)：111－119.

Xu Ximeng, Zheng Fenli, Wu Hongyan, et al. Impacts of cornstalk mulching buffer strip on rill erosion and its hydrodynamic character[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(24): 111－119. (in Chinese with English abstract)

[8] 朱高立，黃炎和，林金石，等. 模擬降雨條件下秸稈覆蓋對崩積體侵蝕產流產沙的影響[J]. 水土保持學報，2015，29(3)：27－31.

Zhu Gaoli, Huang Yanhe, Lin Jinshi, et al. Effect of straw mulch in colluvial soil erosion and yield of runoff and sediment under simulated rainfall[J]. Journal of Soil and Water Conservation, 2015, 29(3): 27－31. (in Chinese with English abstract)

[9] 唐濤，郝明德，單鳳霞. 人工降雨條件下秸稈覆蓋減少水土流失的效應研究[J]. 水土保持研究，2008，15(1)：9－11.

Tang Tao, Hao Mingde, Shan Fengxia. Effects of straw mulch application on water loss and soil erosion under simulated rainfall[J]. Research of Soil and Water Conservation, 2008, 15(1): 9－11. (in Chinese with English abstract)

[10] Giménez R, Govers G. Effects of freshly incorporated straw residue on rill erosion and hydraulics[J]. Catena, 2008, 72(2): 14－223.

[11] 林超文，羅春燕，龐良玉，等. 不同耕作和覆蓋方式對紫色丘陵區坡耕地水土及養分流失的影響[J]. 生態學報，2010，30(22)：6091－6101.

Lin Chaowen, Luo Chunyan, Pang Liangyu, et al. Effects of different cultivation and mulching methods on soil erosion and nutrient losses from a purple soil of sloping land[J]. Acta Ecologica Sinica, 2010, 30(22): 6091－6101. (in Chinese with English abstract)

[12] Giménez R, Govers G. Flow detachment by concentrated flow on smooth and irregular beds[J]. Soil Science Society of America Journal, 2002, 66(5): 1475－1483.

[13] Dunkerley D. Flow threads in surface run-off: Implications for the assessment of flow properties and friction coefficients in soil erosion and hydraulics investigations[J]. Earth Surface Processes and Landforms, 2004, 29(8): 1011－1026.

[14] Foster G R, Johnson C B, Moldenhauer W C. Hydraulics of failure of unanchored cornstalks and wheat straw mulch for erosion control[J]. Transactions of the ASAE, 1982, 25(4): 940－947.

[15] 張寬地，王光謙，孫曉敏，等. 坡面薄層水流水動力學特性試驗[J]. 農業工程學報，2014，30(15)：182－189.

Zhang Kuandi, Wang Guangqian, Sun Xiaomin, et al. Experiment on hydraulic characteristics of shallow open channel flow on slope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(15): 182－189. (in Chinese with English abstract)

[16] 張光輝. 坡面水蝕過程水動力學研究進展[J]. 水科學進展，2001，12(3)：395－402. Zhang Guanghui Advances in study of runoff detachment processes based on hydraulics[J]. Advances in Water Science, 2001, 12(3): 395－402. (in Chinese with English abstract)

[17] 吳淑芳，吳普特，原立峰. 坡面徑流調控薄層水流水力學特性試驗[J]. 農業工程學報，2010，26(3)：14－19.

Wu Shufang, Wu Pute, Yuan Lifeng. Hydraulic characteristics of sheet flow with slope runoff regulation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(3): 14－19. (in Chinese with English abstract)

[18] 吳秋菊，吳發啟，王林華. 土壤結皮坡面流水動力學特征[J]. 農業工程學報，2014，30(1)：73－80.

Wu Qiuju, Wu Faqi, Wang Linhua. Hydrodynamic characteristics of overland flow under soil crusts condition[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(1): 73－80. (in Chinese with English abstract)

[19] 楊錦，呂宏興，上官周平. 薄層水流水力特性試驗研究[J]. 灌溉排水學報，2008，27(4)：58－60.

Yang Jin, Lü Hongxing, Shangguan Zhouping. Experiment study on hydraulic properties of overland flow[J]. Journal of Irrigation and Drainage, 2008, 27(4): 58－60. (in Chinese with English abstract)

[20] 張寬地，王光謙，孫曉敏，等. 模擬植被覆蓋條件下坡面流水動力學特性[J]. 水科學進展，2014，25(6)：825－834.

Zhang Kuandi, Wang Guangqian, Sun Xiaomin, et al. Hydraulic characteristic of overland flow under different vegetation coverage[J]. Advances in Water Science, 2014, 25(6): 825－834. (in Chinese with English abstract)

[21] 江忠善，宋文經. 坡面流速的試驗研究[J]. 中國科學院西北水土保持研究所集刊，1988(7)：46－52.

[22] Horton R E. Erosional development of streams and their drainage basins: Hydrophysical approach to quantitative morphology[J]. Geological Society of America Bulletin, 1945, 56(3): 275－370.

[23] 丁文峰，李占斌，丁登山. 坡面細溝侵蝕過程的水動力學特征試驗研究[J]. 水土保持學報，2002，16(3)：72－75.

Ding Wenfeng, Li Zhanbin, Ding Dengshan. Study on hydrodynamic characters of runoff in rill erosion process on slope[J]. Journal of Soil and Water Conservation, 2002, 16(3): 72－75. (in Chinese with English abstract)

[24] Foster G R, Huggins L F, Meyer L D. A laboratory study of rill hydraulics: I. Velocity relationships[J]. Transactions of the ASABE, 1984, 27(3): 790－796.

[25] 王俊杰，張寬地，龔家國，等. 不同覆蓋度條件坡面水流阻力規律[J]. 水土保持學報，2015，29(5)：1－6.

Wang Junjie, Zhang Kuandi, Gong Jiaguo, et al. Overland flow resistance law under different vegetation coverage[J]. Journal of Soil and Water Conservation, 2015, 29(5): 1－6. (in Chinese with English abstract)

[26] Pan C Z, Shangguan Z P, Lei T W. Influences of grass and moss on runoff and sediment yield on sloped loess surfaces under simulated rainfall[J]. Hydrological Processes, 2006, 20(18): 3815－3824.

[27] Hsieh P, Bolton S. Laminar surface water flow over vegetated ground[J]. Journal of Hydraulic Engineering, 2007, 133(3): 335－341.

[28] 郭雨華，趙廷寧，孫保平，等. 草地坡面水動力學特性及其阻延地表徑流機制研究[J]. 水土保持研究，2006，13(4)：264－267.

Guo Yuhua, Zhao Tingning, Sun Baoping, et al. Study on the dynamic characteristics of overland flow and resistance to overland flow of grass slope[J]. Research of Soil and Water Conservation, 2006, 13(4): 264－267. (in Chinese with English abstract)

[29] 李勉，姚文藝，楊劍鋒，等. 草被覆蓋對坡面流流態影響的人工模擬試驗研究[J]. 應用基礎與工程科學學報，2009，17(4)：513－523. Li Mian, Yao Wenyi, Yang Jianfeng, et al. Experimental study on the effect of grass cover on the overland flow pattern in the hillslope-gully side erosion system[J]. Journal of Basic Science and Engineering, 2009, 17(4): 513－523. (in Chinese with English abstract)

[30] 王玲玲，姚文藝，申震洲，等. 草被覆蓋度對坡面流水力學參數的影響及其減沙效應[J]. 中國水土保持科學，2009，7(1)：80－83.

Wang Lingling, Yao Wenyi, Shen Zhenzhou, et al. Effects of grass coverage on shallow flow hydraulic parameters and sediment reduction[J]. Science of Soil and Water Conservation, 2009, 7(1): 80－83. (in Chinese with English abstract)

[31] 肖培青，姚文藝，申震洲，等. 草被覆蓋下坡面徑流入滲過程及水力學參數特征試驗研究[J]. 水土保持學報，2009，23(4)：50－53.

Xiao Peiqing, Yao Wenyi, Shen Zhenzhou, et al. Study on runoff and infiltration process and hydraulic parameters characteristics with grass coverage[J]. Journal of Soil and Water Conservation, 2009, 23(4): 50－53. (in Chinese with English abstract)

[32] 張冠華. 茵陳蒿群落分布格局對坡面侵蝕及坡面流水動力學特性的影響[D]. 楊凌：西北農林科技大學，2012.

Zhang Guanhua. Influence of Patterned Artemisia Capillaris on Slope Erosion and Overland Flow Hydrodynamic Characteristics[D]. Yangling: Northwest A&F University, 2012. (in Chinese with English abstract)

[33] 史志華，張風寶. 土壤侵蝕動力過程與調控[M]. 北京：科學出版社，2015.

[34] 張樂濤. 工程堆積體陡坡土壤侵蝕動力過程試驗研究[D]. 楊凌：中國科學院研究生院，2013.

Zhang Letao. Experimental Study on Dynamic Process of Soil Erosion of Engineering Accumulation in Steep Slope During the Process of Engineering Construction[D]. Yangling: The University of Chinese Academy of Sciences, 2013. (in Chinese with English abstract)

Hydraulic characteristic of overland flow on loess farmland slope under mulch with different wheat straw lengths

Li Chaodong1,2, Li Zhanbin1, Ma Jianye1, Ma Bo1※

(1.712100; 2.712100)

Straw mulching is one of the important conservation tillage measures, and has a positive effect on reducing soil erosion on slopping farmland. In order to study the impact of straw length on the hydrodynamic characteristics of overland flow, indoor artificial fixed bed erosion resistance experiments were carried in the Soil and Water Conservation Engineering Laboratory in Northwest A&F University. The study area was located in the southern Loess Plateau, Shaanxi Province, China (34°14′–34°20′N, 107°59′–108°08′E). By the scour trough, the variations of hydrodynamic parameters of overland flow such as flow velocity, runoff depth, Reynolds number, Froude number, Darcy-Weisbach resistance coefficient and Manning roughness coefficient were determined at 4 slope gradients (3°, 5°, 10° and 15°), 5 flow rates (3.0, 4.5, 6.0, 7.5 and 9.0 L/min) and 3 wheat straw lengths (3-5, 8-10 and 13-15 cm). The slope without straw mulch was considered as the control (CK). The results showed that under the experimental condition, the wheat straw length greatly affected the hydrodynamic parameters. Under straw mulching with different straw lengths, Reynolds number ranged from 166 to 558, the overland flow was primarily laminar and subcritical flow when the flow rate was less than or equal to 7.5 L/min, and transformed into transition and subcritical flow while the flow rate was 9.0 L/min. Under the conditions of straw mulching, compared with bare slope, the overland flow ran with smaller flow velocity and higher runoff depth, Darcy-Weisbach resistance coefficient and Manning roughness coefficient. Straw length affected the runoff depth. Compared with the bare slope, the runoff depth of overland flow on slope with the straw length was increased by 107.69%-178.64%. Straw mulching could decrease the flow velocity by 48.30%-60.64% than the CK. The runoff depth decreased as the straw length increased, and the flow velocity increased with the increase in straw length. The flow state index of overflow with straw mulching was reduced with increased of straw length and slope gradient. When the straw length increased from 3-5 cm to 13-15 cm, the average flow state index was decreased from 0.716 to 0.501. The Darcy-Weisbach resistance coefficient was 0.14-1.47 on bare slope, 2.67-8.00 on slope with 3-5 cm straw mulch, 2.39-6.36 on slope with 8-10 cm straw mulch and 1.31-4.28 on slope with 13-15 cm straw mulch, respectively. It averagely increased by 20.89, 15.69 and 8.17 times respectively compared with the bare slope, indicating that the straw mulch could strongly resist the flow movement. Darcy-Weisbach resistance coefficient was decreased with the increase in straw length. The Manning roughness coefficient was 0.013-0.044 on the bare slope, 0.071-0.132 on slope with 3-5 cm straw mulch, 0.067-0.118 on slope with 8-10 cm straw mulch and 0.049-0.083 on slope with 13-15 cm straw mulch, respectively. It increased by 4.38, 3.64 and 2.33 times respectively compared with the bare slope. Manning roughness coefficient was reduced with straw length increased. This study could provide reliable information for clarifying the mechanism behind soil erosion with straw mulching on farmland slope.

flow patterns;hydrodynamics; straw; overland flow; fixed bed test; Loess farmland

10.11975/j.issn.1002-6819.2017.04.022

S157.1

A

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

2016-04-15

2016-10-10

國家自然科學基金（41401309、41330858）；西北農林科技大學博士科研啟動基金（2452015345）

李朝棟，男，陜西千陽人，主要從事黃土高原土壤侵蝕試驗研究。楊凌西北農林科技大學水土保持研究所，712100。 Email：soilcrop@163.com

馬 波，男，陜西西鄉人，博士，助理研究員，主要從事土壤侵蝕與水土保持研究。楊凌西北農林科技大學水土保持研究所，712100。 Email：mabo_916@126.com

李朝棟，李占斌，馬建業，馬 波. 不同長度小麥秸稈覆蓋下黃土耕地坡面流水動力學特性[J]. 農業工程學報，2017，33(4)：153－160. doi：10.11975/j.issn.1002-6819.2017.04.022 http://www.tcsae.org

Li Chaodong, Li Zhanbin, Ma Jianye, Ma Bo. Hydraulic characteristic of overland flow on loess farmland slope under mulch with different wheat straw lengths[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(4): 153－160. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.04.022 http://www.tcsae.org