呂春娟，張 徐，畢如田，牛耀彬，夏 露，郭星星，陳 丹

（山西農業(yè)大學資源環(huán)境學院，農業(yè)資源與環(huán)境國家級實驗教學示范中心，太谷030801）

0 引 言

2019 年9 月18 日習總書記在黃河流域生態(tài)保護和高質量發(fā)展座談會上明確將黃河流域生態(tài)環(huán)境和高質量發(fā)展上升為重大國家戰(zhàn)略，強調“治理黃河，重在保護，要在治理”。山西省作為黃河流域中游的流經省份，河谷縱橫，地貌類型復雜多樣，夏季忽晴忽雨、晴雨夾雜的天氣經常出現，且多暴雨，一年中土壤侵蝕量往往由幾次較大的暴雨造成[1]；同時山西省也是礦產資源大省，鐵礦儲量居全國第四位，礦產資源的開發(fā)利用嚴重擾動了原地表土壤、破壞植被，產生大量的廢棄土石混合堆積體，組成成分復雜，結構松散、黏聚力差[2]，極易成為黃河流域泥沙來源，加速侵蝕，因此減少水土流失，加快礦區(qū)生態(tài)保護和修復已成為山西省及黃河流經市縣生態(tài)保護和高質量發(fā)展的重點任務[3]。

工程堆積體是由生產建設活動排棄的廢棄物重塑而形成的松散堆積體，外在的地貌形態(tài)和內在物質組成等與原生地表都有很大的區(qū)別，其土壤侵蝕速率和產沙量遠超原地貌[2,4]。基于野外各類工程堆積體的調查[5]和室內模擬沖刷試驗[4]發(fā)現工程堆積體的侵蝕形式因物質組成不同有所差異，依據堆積體的物質組成和顆粒大小，可以分為土質堆積體、土多石少混合堆積體、土少石多混合堆積體石質堆積體、尾礦堆積體等。其中工程建設活動產生的土石混合堆積體的侵蝕關注度最高，主要圍繞石礫含量與堆積體產流產沙和水動力學展開研究。土壤中的礫石改變了堆積體的流速、坡面粗糙度、弗汝德數流態(tài)等而影響侵蝕過程。不同研究得出的礫石對坡面堆積體侵蝕作用不盡相同[2,6-8]，但均認為可能存在礫石含量、粒徑大小的“臨界值”，小于這個臨界值，會促進坡面侵蝕，高于臨界值，會減小坡面侵蝕[9-10]。也有研究表明礫石的侵蝕影響與降雨強度有關系[11]，或者與礫石在堆積體中分布的位置，在堆積體的表面還是處在堆積體內部有關[12]。張樂濤等[13-14]等研究了以土質為主體的高速公路堆積體的水力學特性及治理措施對土壤剝蝕率的影響。這些研究中堆積體坡面均是在系列放水流量或者系列雨強單一作用下的侵蝕過程觀測，野外坡面往往要經歷多次間歇性非均勻降雨或者非均勻連續(xù)降雨的侵蝕，因此也有學者關注持續(xù)降雨條件或者雨型改變對坡面侵蝕的影響，Parsons 等[15-16]根據雨強變化均設計了均勻型、上升型、下降型、上升下降型和下降上升型5 種暴雨模式；呂佼容等[10]研究了連續(xù)3 次恒雨強條件下礫石含量對堆積體剝蝕率和坡面侵蝕的影響；張樂濤等[17]在總來水量一致條件下，對比了連續(xù)均勻性與3 種連續(xù)性非均勻上方來水對土質堆積體坡面侵蝕過程的影響；霍云云等[18]對土質堆積體坡面在7 次恒雨強降雨下的細溝動態(tài)發(fā)育過程進行了詳細的監(jiān)測。持續(xù)降雨更接近真實氣象條件，易于反映坡面的侵蝕發(fā)展演化；研究發(fā)現，與均勻降雨相比，非均勻降雨條件增加了坡面侵蝕的時空隨機性，對總徑流量影響不大，但是顯著影響坡面產沙量。對于堆積體侵蝕的水動力學研究，大多文獻從坡面流流速、阻力系數、徑流剪切力、水流功率、過水斷面單位能量、弗汝德數、雷諾數等多個參數去分析，可能因研究對象和試驗條件的改變，不同堆積體坡面侵蝕描述的水動力參數不同，但大多數研究結果都認同的是水流功率能很好地刻畫侵蝕產沙量，對產流量的水動力學研究較少[7-8,11,13-14,17]。

綜上，目前堆積體的侵蝕研究：1）多集中在公路鐵路等工程建設活動產生的棄土棄渣；2）多以產流產沙研究為主，對于坡面形態(tài)變化的觀測較少；3）單一降雨條件較多，有少量連續(xù)降雨條件但土壤多于飽和狀態(tài)或者高含水率狀態(tài)，土壤在自然含水率的非飽和狀態(tài)下的間歇性降雨研究較少；4）基質類型多以純土或者土石混合體為主，類似尾礦的基質研究很少。研究區(qū)鐵尾礦堆積體，顆粒細，與土質堆積體或者巖土堆積體組成有很大差異，且坡度陡，入滲率速度慢，有機質缺乏，坡面溝壑縱橫，嚴重制約生態(tài)恢復，初步研究表明尾礦中添加土壤和菌糠可以顯著提高其滯留貯水量，暴雨時快速吸納積水[19]。根據研究區(qū)的降雨特征，模擬雨季自然氣象條件下非均勻間歇性降雨，基于少土復墾思路，選擇添加少量土壤和菌糠作為尾礦改良模式，研究不同重構坡面的產流產沙、坡面微形態(tài)變化及水動力學特征，以期為尾礦坡面的生態(tài)恢復提供理論支撐。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)垣曲縣位于黃土高原西南部，山西省南段 （35°14′14″～35°14′54″N，111°40′24″～111°41′09″E），境內地勢西北高、東南低，海拔范圍176.7～2 358 m，高差大，主要河流有亳清河、允西河、板澗河、西陽河、五福澗河等，均由北而南注入黃河；屬于溫帶大陸性季風氣候，年均降雨量631 mm，季節(jié)分布不均，主要集中在6—8 月，且短歷時暴雨較為頻發(fā)，存在水土流失的潛在客觀條件。亳清河流域是垣曲縣面積最大、人口最多、耕地最多的流域，也是小浪底庫區(qū)的主要匯水區(qū)之一；上游零散分布有很多鐵礦區(qū)，試驗基地為位于亳清河流域中游的泉子溝干排尾礦庫，尾礦砂自然堆積坡度大多在30°～40°之間，表層容重1.70 g/cm3左右，水分入滲率遠低于農田土壤入滲率[19]，坡面侵蝕嚴重，因此尾礦砂的水土流失治理和生態(tài)修復對于小浪底庫區(qū)的水質和生態(tài)環(huán)境具有非常重要的意義。

1.2 試驗材料與試驗設計

試驗材料包括鐵尾礦砂、土壤和菌糠。鐵尾礦砂采集于泉子溝干排尾礦庫，土壤采自尾礦周圍的農田，菌糠來源于當地蘑菇種植基地，采集的材料均置于陰涼處，自然風干后，尾礦砂和土壤過2 mm 篩，碾碎后的菌糠過5 mm 篩，以備裝填土槽所用。

工程堆積體坡度一般在30°～40°之間[13,20]，結合對鐵尾礦區(qū)的調查，試驗設計25°、30°、35°3 個坡度；根據黃土高原點雨量站記錄的暴雨強度[21]和汛期多暴雨的降雨特點，同時參考土高原侵蝕相關文獻中的雨強變化[20,22-23]，本研究設置3 場間歇性非均勻降雨，降雨強度依次為60、90、120 mm/h（定義為第1 次，第2次，第3 次），在同一基質坡面上依次間斷進行，每2場降雨之間的時間間隔不固定，主要依據坡面水分風干的快慢，保證每次降雨前表層的含水率都在8%左右（接近尾礦砂夏季自然含水率），一般間隔約1 周；中國尾礦區(qū)大多是土石山區(qū)，復墾土源缺乏，因此選擇添加少量土壤或者添加菌糠的無土和少土基質改良模式，共6 種鐵尾礦砂重構坡面：松散鐵尾礦砂（LT），模擬新排棄的尾礦坡面；結殼鐵尾礦砂（CT），模擬長期堆積壓實的尾礦坡面；鐵尾礦砂與菌糠混合物（質量比10∶1，TM），鐵尾礦砂 、菌糠和土壤混合物（質量比10∶1∶4，TSM），鐵尾礦砂和土壤混合物（質量比2.6∶1，TS），純土（S）模擬鐵尾礦添加外源物后的坡重構面；每個處理重復3次。根據國際土壤質地分類制，LT 和CT 為砂質壤土，TM、TSM 和TS 均為砂質黏壤土，S 為黏壤土。

所有試驗材料按照試驗設計比例混合均勻后，以5 cm 一層分層裝填的方式裝填至設計坡度，下部為尾礦砂，表層20 cm 為6 種重構基質。參照野外復墾小區(qū)容重實測值[19]，設計LT、CT、TM、TSM、TS 和S 裝填容重分別為1.65、1.65、1.35、1.4、1.5、1.3 g/cm3；為保證垂直剖面均勻性，層間采用鋼絲刷打毛以消除分層效應。

1.3 間歇模擬降雨試驗方法

間歇模擬降雨試驗在山西農業(yè)大學資源環(huán)境學院試驗站進行。降雨模擬系統為QYJY-501 全自動下噴式模擬降雨器，由西安清遠測控技術有限公司定制，通過調控噴頭，降雨強度可在15～200 mm/h 范圍內變化；試驗土槽長為2 m，寬為1 m，前端高為0.6 m、后端高為2.3 m，通過前后端試驗材料裝填深度，坡度可在0～40°之間調節(jié)（圖1）。每場降雨過程中均以坡面產流開始計時，以3 min 為間隔，在土槽出水口連續(xù)收集30 min 徑流、泥沙樣，樣品充分沉降后，用量筒量測徑流體積，計算徑流率；泥沙樣在105 ℃烘干稱質量，以計算產沙率，流速用高錳酸鉀示蹤法測定，降雨開始前后分別測定試驗水溫，計算水流的運動黏滯系數[24]。

圖1 人工模擬降雨示意圖 Fig.1 Schematic diagram of artificial simulated rainfall

1.4 數據處理分析

1.4.1 坡面形態(tài)變化

坡面形態(tài)微變化采用測釬法、網格坐標法和照相法相結合。第1 次降雨前在坡面上以10 cm×10 cm 網格布設測釬，保證測釬頂部和坡面剛好齊平，每次降雨結束后采用精度為1 mm 的鋼尺測量各測釬三維坐標，如果有細溝在網格中間出現，則采用坐標加密方式對細溝的坐標進行測量；利用AutoCAD 2010 軟件對照片數據矢量化，再根據測量數據進行校核，形成坡面及細溝的閉合圖，用來計算面積及長度。選擇地形起伏度RA（mm）、溝壑密度GD（km/km2）和溝壑面積GA（%）3 個指標來反映降雨過程中的坡面微地形變化[25]。

式中 iH 為測釬出露地表的高度，mm；m 為測釬數量；A 代表細溝面積，m2；S 為坡面面積，m2；L 為坡面細溝總長度，m；1 000 為單位換算系數。

1.4.2 坡面水動力學參數計算

坡面流是坡面侵蝕產沙的主要動力，水層一般很薄，并且受下墊面影響，水深、流速沿程不斷變化；水流底面坡度一般也比較陡峭，在一定條件下便會失穩(wěn)，產生滾波，并且產流過程中受降雨雨滴的擾動，水流在時間和空間上都是變化的，這些復雜因素交織在一起，使得對坡面水流的研究有很大難度[26-27]。目前大多數學者的做法是將坡面流過程進行簡化處理，應用明渠流理論進行研究[13-14,27]。本研究中用雷諾數Re和弗如德數Fr來判斷坡面水流流態(tài)，用Darcy-weisbach 阻力系數f 和Manning 粗糙系數n 作為坡面流阻力指標，坡面流流速ν 、徑流剪切力τ、水流功率ω 作為侵蝕動力指標，來分析不同基質在間歇性降雨條件下的侵蝕水動力特征差異。

式中h 為坡面平均徑流深，m；Q 為坡面產流量，L；S為土槽面積，m2；Re雷諾數，無量綱；V 為坡面平均流速，m/s；R 為水力半徑，m，在薄層水流下用水深h 代替R；υ 為流體運動黏滯系數，cm2/s，t 為水溫，℃；rF為坡面流弗汝德數，無量綱；f 為Darcy-Weisbach 阻力系數；n 為Manning 粗糙系數；ω 為水流功率，N/(m·s)；g為重力加速度9.8 m/s2；τ 為徑流剪切力，N/m2；γ 為水的重度，N/m3；ρ 為水的容重，kg/m3；J 為水力坡度，取坡度的正弦值。

2 結果與分析

2.1 坡面產流變化

3 種坡度重構坡面產流率變化趨勢相似，以30°重構坡面為例闡述產流率隨降雨場次的變化，見圖2。3 次降雨過程中，6 種基質的產流率均隨著產流歷時波動增加，且隨著降雨次數的增加，產流率增大；但是3 次降雨中產流率的增幅（產流率均值之差）不同，前2 次降雨各基質的增幅為0.55～2.52 L/min；后2 次降雨雖然雨強增大，但是產流率并沒有按比例增加，增幅反而減小，均小于1.77 L/min，主要是因為第1 次降雨時，基質表面整體比較疏松，基質入滲能力大，坡面主要是濺蝕和面蝕，第1 次降雨結束后，表面被擊實，形成結皮，基質孔隙也部分被堵塞，降雨入滲減少，而且降雨強度增大，結合觀察到的坡面形態(tài)變化可知，第2 次降雨過程很多坡面出現大量細溝，因此綜合作用下，產流率大幅增加；第3 次降雨，坡面流主要集中于細溝，由于細溝的溝壁坍塌，溝道下部也出現很多泥沙沉積物，溝道內坡度反而變緩，因此溝道徑流有了更多入滲機會，產流率增幅也相應減小。

不同降雨場次中，6 種重構坡面間產流率的相對大小也是變化的（圖2）。第1 次降雨中，25°、30°和35°3個坡度下重構坡面的平均產流率變幅為 0.49 ～0.65 L/min，CT、TM 和TS 的產流率最高；第2 次降雨中，3 個坡度下6 種基質間的平均產流率變幅增大為1.12～1.89 L/min；第3 次降雨中，平均產流率變幅分別為1.47～1.73 L/min；第2、3 次降雨中，LT、TM 和TS表現出較高的產流率，S 表現出較低的產流率。根據3 個坡度下多次降雨累計產流量發(fā)現，TM 的累計產流量最高，S 累計產流量最低；但方差分析表明，任何坡度下，重構坡面對產流率影響都不顯著。

2.2 坡面產沙變化

產沙率變化與產流率變化完全不同。在3 個坡度下，隨著降雨場次增加，產沙率并沒有持續(xù)增大，以30°重構坡面產沙率為例（圖3），不同基質表現不同：LT、CT 和TS 產沙率整體表現出持續(xù)快速增大的趨勢，TM、TSM 和S 在大多數試驗條件下變化緩慢甚至有減少的趨勢。說明在同樣間歇降雨條件下，不同重構坡面基質抵抗徑流剝蝕的能力不同，坡面侵蝕基本上經歷擊濺侵蝕-面蝕-跌坎-細溝侵蝕4 個階段，試驗中觀測到 LT、CT和TS 坡面在第2 次降雨中已進入細溝的發(fā)展階段，第3次降雨細溝繼續(xù)向長向寬發(fā)展，產沙率快速增大，而TM、TSM 和S 基本上還處在面蝕-跌坎的過渡階段或者細溝的初始階段，而且這個階段一直持續(xù)到第3 次降雨，發(fā)展特別緩慢，因此產沙率增加比較慢甚至減小。

圖2 30°重構坡面產流率隨降雨場次的變化 Fig.2 Runoff rate variations of 30° reconstructed slopes with rainfall events

在每場降雨中，隨著產流歷時增加，各重構坡面的產沙率波動變化，變異系數最大為70.7%，最小為7.4%，與大多數土質或者偏土質工程堆積體類似[4,8]，由于細溝中的重力侵蝕頻發(fā)，導致泥沙突然增加或減少，出現波動，但是本研究中不同的是6 種重構坡面產沙率的波動變化趨勢不同，有3 種類型：單峰型、增加型和減少型。S 在25°坡度下第2 次降雨，30°和35°坡度下的第2、3次降雨整體表現為單峰型；TM 在第2 次和第3 次降雨中多表現為減小型；其余情況下各坡面產流率隨產流歷時變化趨勢均為增加型。重構坡面產沙率也表現出相似的分異現象（圖3），方差分析表明重構坡面類型對產沙率影響顯著（P<0.05）。為了更科學地探索產沙的分異規(guī)律，依據產沙率對重構坡面進行系統聚類分析，6 種重構坡面可分為2 類（圖4）：1）LT、CT 和TS，表現為高產沙型；2）TM、TSM 和S，表現為低產沙型。6 種重構坡面中TM 的產沙量最低，LT 最高，但CT 幾乎在所有坡度和所有降雨場次中，產沙率都要低于LT，累計產沙量約為LT 的67%，進一步說明下墊面對產沙有很大的影響；煤礦區(qū)的棄土棄渣（35°和40°）也有近似的結果，不同下墊面的輸沙速率由高到低為棄土體、沙多石少棄渣體、沙少石多棄渣體、煤矸石，棄土體的輸沙速率最高，是煤矸石的100～1 000 倍[28]。

圖3 30°重構坡面產沙率隨降雨場次變化 Fig.3 Variations of sediment yield rate of 30° reconstructed slopes with rainfall events

大多數重構坡面在坡度25°下的產沙量最大，坡度增大產沙量反而減小。這可能與很多文獻中水力侵蝕的臨界坡度有關，王瑄等[29]在3°～30°間設置10 個塿土坡度梯度的放水沖刷試驗，發(fā)現21°左右是徑流剝蝕土壤的臨界坡度，小于這個坡度產沙率隨坡度增大，大于這個坡度產沙率隨坡度減小。因此推測本研究中的基質在25°和30°之間可能存在土壤流失的臨界坡度，有待于更多的坡度梯度試驗來進一步驗證。

圖4 重構坡面產沙量聚類樹狀圖 Fig.4 Clustering tree of sediment yield of reconstructed slopes

2.3 坡面微地形侵蝕特征

觀察圖5a～圖5c 重構坡面溝壑密度和圖5d～圖5f 溝壑面積的變化，發(fā)現3 種坡度下，第1 次降雨后的溝壑密度和溝壑面積都很小，第2 次降雨后均急劇增加，但第3次降雨后，大多數重構坡面溝壑密度表現為減小，溝壑面積表現為增加。這與細溝的形態(tài)特征有很大的關系，LT、CT 和TS 的細溝寬深比一般比較小，細溝狹長，上下易貫通，S 的細溝寬深比較大，細溝一般中部最寬，上下不易貫通，TM 和TSM 坡面多為小坑穴，偶然局部有少量細溝，但都比較淺。大多坡面是在第2 次降雨后產生細溝并發(fā)育，細溝在向長、向寬、向深的三向發(fā)展過程中，會有相鄰細溝的合并，上下細溝的貫通和溝道的橫向擴展，因此第3 次降雨后，溝壑面積均增加，但溝壑密度不一定增加。

圖5 重構坡面微地形特征隨降雨場次的變化 Fig.5 Microtopography variations of reconstructed slopes with rainfall events

隨著降雨場次的增加，各坡面的地形起伏度均增大，但增大速率不同，第1 次降雨中，3 個坡度下的地表起伏度差別不大，但是隨著降雨場次增加，LT、CT 和TS 坡面的起伏度增加速度遠大于TM、TSM 和S，特別是在30°和35°下，2 組重構坡面間的分異更明顯，LT、CT 和TS 起伏度變化范圍為10.2～38.9 mm，TM、TSM 和S 起伏度變化僅在3.2～10.6 mm 范圍內。主要是因為鐵尾礦砂抗蝕性差，陡坡下的LT、CT 和TS 極易產生細溝，造成坡面的高低起伏，坡面徑流匯集細溝，導致細溝的縱向發(fā)展，進一步加劇坡面的起伏度。S 坡面主要是因為強大的入滲能力，徑流流量和流速均減小，坡面細溝發(fā)展緩慢，因此隨降雨場次的增加，起伏度也沒有太大的變化。TM、TSM 坡面第1 次降雨后，表層尾礦顆粒流失之后，坡面殘留的菌糠類似不透水的小片石鑲嵌在尾礦中形成均勻分布的菌糠結殼和小坑穴結構，坡面粗糙度增加，坡面徑流被約束形成“網狀流”[30]沿著菌糠間隙曲折流動，路徑延長，流速減小，動能降低，雖然整個坡面上都有水流在流動，但由于無法匯集成大的股流，因此坡面上的細溝很少，類似鱗片狀面蝕，起伏度變化也很?。挥绕銽M在30°和35°坡度下3 場次降雨中，均沒有明顯細溝產生。

重構坡面對侵蝕形態(tài)影響很大。根據溝壑面積、溝壑密度和地表起伏度微地形指標，對6 種重構坡面進行系統聚類，可以分為2 類（圖6）：1）LT、CT 和TS，侵蝕形態(tài)嚴重；2）TM、TSM 和S，侵蝕形態(tài)輕微。坡面基質類型對侵蝕形態(tài)的影響在堆積體研究中也有發(fā)現，礫石堆積體在連續(xù)降雨條件下，隨著表土剝離礫石覆蓋率增大，坡面沒有細溝產生[10]；而黃壤堆積體一般形成單一主溝，紫色土堆積體則形成較多細溝[4]。

比較3 個坡度之間微地形侵蝕特征，沒有明顯變化趨勢，也并非坡度越小，溝壑面積和溝壑密度和起伏度越小，甚至有些重構坡面在25°坡度下反而比30°和35°下侵蝕特征嚴重，與前面的產沙量變化一致，可能也是受臨界坡度的影響。

圖6 重構坡面微地形聚類分析樹狀圖 Fig.6 Clustering tree of microtopography of reconstructed slopes

2.4 坡面侵蝕的水力學特征分析

2.4.1 坡面流流速

3 種坡度下的流速變化趨勢相似，以30°坡面流速變化圖為例：隨降雨次數增加，6 種基質的坡面流平均流速整體呈現增大趨勢，LT、CT 和TS 增速快；TM、TSM和S 增速稍慢，受下墊面影響個別條件下流速甚至會減?。▓D7）。在每場降雨中，隨著產流歷時延長，流速呈現波動變化；坡面不同，變化趨勢不同，有的是波動增加，有的平穩(wěn)波動，特別是LT、CT 和TS 的流速波動幅度最大，流速也較大；相對而言，整個降雨過程中TM 流速平穩(wěn)，在最后1 次降雨中流速也最小。在降雨情況相同的條件下，一般坡面流的流速小于細溝流的流速，試驗中流速觀測，是在坡面微形態(tài)不斷演化過程中測定的，降雨過程中LT、CT 和TS 的坡面侵蝕產沙量最高，坡面的溝壑密度、溝壑面積和起伏度等侵蝕形態(tài)也最嚴重，因此觀測的流速是坡面流向細溝流的轉化過程，特別是第2 次和第3 次降雨中，坡面細溝進入到快速發(fā)展和擴張階段，坡面測定的流速大多是細溝流的流速，而細溝流流速受細溝的發(fā)展影響，側壁崩塌、波涌等現象，流速出現短暫的波動增大或減小趨勢；而TM 在第1 次降雨條件下，表層“浮砂”被沖刷后，自然風干過程中形成“菌糠結殼”，在第2、3 次降雨中，如前所述則會延長徑流流路，消耗能量，因此流速小而且平穩(wěn)。不僅6種重構坡面的徑流流速變化過程不同，方程分析表明6種重構坡面間平均流速差異也顯著（P<0.05）。

圖7 30°重構坡面流流速隨降雨場次變化 Fig.7 Change of flow velocity on 30° reconstructed slopes with rainfall events

2.4.2 坡面流流態(tài)分析

所有坡度所有場次降雨中，各坡面水流的雷諾數變化范圍在0.95～74.88（表1），均小于500，表明黏滯力對流場的影響大于慣性力，流場中流速的擾動會因黏滯力而衰減，流體流動穩(wěn)定，坡面流流型為層流，但隨著降雨場次的增加，特別是第2 次和第3 次降雨中，坡面細溝出現，部分坡面細溝甚至上下貫通，促使流速加快，因此雷諾數表現出增加的趨勢，也表明慣性力相對于黏滯力對坡面流流場的影響越來越大，增強了水流的局部紊動性，在坡面LT 和CT 細溝中也觀測到“徑流浪”現象，驅動更多的泥沙，流型依然為層流，但降雨場次間雷諾數差異顯著（P<0.05）。而在15°平直黃土坡面或者魚鱗坑坡面在5 場間歇性模擬降雨試驗中，隨著降雨場次增加，坡面流均由層流轉變?yōu)槲闪鱗31-32]。與本研究結果不同，分析原因發(fā)現降雨場次增加后，細溝流流速均增大，但文獻中坡面徑流深增加的程度遠大于本研究基質，導致慣性力相對于黏滯力大幅增加，呈現紊流流態(tài)；同樣是土石混合堆積體，25°～40°紫色土堆積體和黃壤堆積體坡面流大多是紊流[4]，北方風沙土堆積體和贛北典型紅壤堆積體在25°坡面流大多是層流[8,11]。坡面流受很多因素影響，不同的試驗條件，研究結果不同，目前關于坡面流流型尚沒有統一的結論。判斷坡面薄層水流的流態(tài)的弗如德數，在大多數基質中均大于1，為急流；只有S、TM、TMS 坡面和TS 在極個別情況下Fr小于1，為緩流（表1）。隨降雨場次和坡度增加，同一重構坡面的弗如德數沒有明顯的變化規(guī)律，但降雨場次和坡度對弗如德數影響顯著（P<0.05）。重構坡面間弗如德數差異顯著（P<0.05），雷諾數差異不顯著；但觀察表1中3 個坡度下的雷諾數和弗如德數變化，均表現出相似的波動性，說明重構坡面對雷諾數也具有較大的影響。

表1 重構坡面的流態(tài)參數與坡面流阻力特征 Table 1 Flow parameters and resistance characteristics of surface flow of reconstructed slope

2.4.3 坡面流阻力特征分析

坡面流阻力反映外界影響因素對坡面薄層水流阻滯能力的大小，是各種因素共同作用的一個綜合系數，其影響因素通常包括地貌形態(tài)、土壤性狀、降雨、礫石、表面覆蓋物的覆蓋程度以及水流運動邊界的形態(tài)特征等[27]，坡面流阻力最大，消耗于坡面侵蝕的能量則越小。

同一坡面，隨著降雨場次增加，阻力系數和粗糙系數基本呈現減少的趨勢（表1，P<0.05），主要是因為隨著降雨場次增加，坡面的起伏程度加劇，整個坡面的粗糙度增加，但是在后期的次降雨中，坡面流已不是坡面漫流，而是在坡面細溝或者菌糠間隙均形成固定流路，流速增大，產生阻力的邊界接觸面減小，因此影響水流的阻力還是會減小。重構坡面間的阻力系數和粗糙系數差異不顯著（P>0.05），但因重構坡面物質組成不同，坡面間的阻力系數還是存在一定程度的差異，TM 和TSM和S 阻力系數和粗糙系數較大，LT、CT 和TS 坡面流阻力系數和粗糙系數相對較低，說明TM、TSM 和S 坡面水流流動過程中需要消耗的能量大，則用于坡面侵蝕和泥沙輸移的能量就小，因此坡面侵蝕產沙就少；反之，LT、CT 和TS 坡面水流克服阻力消耗的能量小，則用于侵蝕的能量大，產沙大。

2.4.4 坡面流的徑流剪切力和水流功率特征分析

當坡面水流處于運動狀態(tài)時，水土接觸面方向存在剪切力，當剪切力大于土壤臨界剪切應力時，土壤顆粒被分散剝蝕，卷入水流并挾帶出坡面，徑流剪切力是分離土壤的主要動力[33]。重構坡面間的徑流剪切力差異顯著（P<0.05），最為顯著的特征就是與其他重構坡面相比，TM 的徑流剪切力最高（表2），但是通過前面分析知道，TM 的累計產流率是最高的，累計產沙量又是最低的，據此可以推測本試驗中徑流剪切力值可能尚未達到TM 的臨界剪切應力，說明TM 具有較高的抗剪切能力。

水流功率是水流剪切力和流速的函數，水流功率包含了坡度、徑流深和流速的影響，比水流剪切力具有更廣的水動力特性[29]。表2 也可以看出，同一重構坡面水流功率隨著降雨場次增加整體表現為增加，不同重構坡面間水流功率差異不顯著，但6 種重構坡面的水流功率在3 個坡度下均表現出的相似的波動變化，說明重構坡面對水流功率還是具有較大的影響。

表2 重構坡面徑流剪切力和水流功率特征變化 Table 2 Variations of runoff shear stress and runoff power of reconstructed slopes

2.5 坡面侵蝕特征的水力學驅動因素分析

坡面侵蝕特征與水力學參數多為非線性單調關系[7,9,13-14]，基于此，本研究采用Spearman 相關系數來分析侵蝕特征與水力學參數的相關性，進一步從坡面徑流的水動力學角度科學揭示產流產沙的驅動因素，但基于前面不同坡面間顯著的產流產沙差異和試驗中觀測到的一些現象，若6 種重構坡面放到一起，可能會掩蓋某些坡面的特殊屬性，因此按重構坡面類型分別對侵蝕特征和水力學參數進行Spearman 相關分析（表3）。表3 可以看到，LT、CT、TS 和S 的產流率和產沙率顯著正相關（P<0.05），表現為水大沙大；TM 和TSM的產流率和產沙率之間雖然也是正相關，但是不具有統計學意義的顯著性，特別是TM 的產流率和產沙率的相關系數僅為0.033，與前面高產流、低產沙的試驗結果一致。6 種坡面的產流率、產沙率與水力學參數的相關性不盡相同。LT和CT 的產流率、產沙率受水力學參數的影響相同，均與流速、雷諾數、弗如德數、水流功率顯著正相關，與阻力系數和粗糙系數顯著負相關；TSM、TS 和S 的產流率與流速、雷諾數和水流功率為顯著正相關；TM 的產流率僅與雷諾數和水流功率顯著正相關。同時發(fā)現重構坡面的產流率與水流功率均存在極顯著正相關，相關系數為1，究其原因發(fā)現二者都與試驗中的徑流量有著直接的線性函數關系，所以相關性極高，說明水流功率能很好地預測產流過程，不受基質類型的影響。TS 和S 的產沙率與流速、雷諾數和水流功率為顯著正相關；TSM 和TM 的產沙率僅與流速顯著正相關，其他水力學參數影響均不顯著。

表3 重構坡面的侵蝕特征與水力學參數間的相關系數 Table 3 Correlation coefficient between erosion characteristics and hydraulic parameters of reconstructed slopes

坡面微地形的侵蝕表征指標溝壑密度、溝壑面積和起伏度，因坡面基質不同，與產流產沙之間的相關性表現出極大的差異（表3）。LT、CT 和TS 坡面的溝壑面積和起伏度顯著正相關于其產流率和產沙率，而溝壑密度與產流率和產沙率的相關性不顯著。TM、TMS 和TS的溝壑面積、起伏度和溝壑密度與產流產沙的相關性大多不顯著，說明LT、CT 和TS 的產流產沙主要來源于坡面地形的微變化，細溝的發(fā)展促進徑流的匯集和產沙量的增加；TM、TMS 和S 的產流產沙多來源于坡面的“均勻”流失，與地形的高低起伏變化相關性相對不強。

大多數坡面的坡面微地形侵蝕表征指標溝壑密度、溝壑面積和起伏度大多數情況下與水力學參數平均流速、雷諾數、弗如德數、徑流剪切力、水流功率表現為正相關，與阻力系數和粗糙系數表現為負相關；TM 和TMS 的坡面微地形侵蝕表征指標中個別指標與水力學參數間表現出與其他基質相反的相關關系，且沒有明顯的規(guī)律。說明坡面微地形的變化可能受多種水力學因素的綜合作用，單一因素很難改變坡面的變化，且這些相關關系大多達不到顯著水平，因此對于坡面微地形的侵蝕表征指標與水力學參數之間的關系沒有進一步定量預測的必要。

據上述產流產沙與水力學參數間顯著相關分析結果，利用SPSS 中的非線性回歸法篩選各類基質侵蝕特征的最佳預測模型（表4）。因LT 和CT 的擬合方程比較接近，所以將二者合并到一起，表4 中6 種坡面的產流率和產沙率與水力學參數間的預測模型均通過回歸方程的擬合優(yōu)度檢驗、顯著性檢驗和回歸系數的顯著性檢驗，結果表明水力學參數能通過對數、線性和冪函數較好地預測產流產沙，同時發(fā)現基質類型相近的有近似的函數關系，如TM 和TSM 的產沙與流速的冪函數關系，LT和CT 具有相同的水力影響因素，說明重構坡面類型對侵蝕的水動力學參數有一定的影響。關于坡面侵蝕的水動力學影響因素及函數描述有很多學者關注[9,17,28-29,34]，土壤剝蝕率隨水力學相關參數增大主要以線性函數或冪函數關系增加，即使坡度接近，因基質類型不同，水動力學的影響因素和函數關系也不完全相同；可見單獨的水力學參數無法考慮坡面類型差異對剝蝕率的影響[35]，目前已有學者嘗試將土壤黏聚力作為表征土壤性質差異的指標引入到水力學參數預測土壤剝蝕率的定量關系中[36]，基于水動力學參數融合能代表土壤性質的指標預測土壤剝蝕率將進一步增加預測的科學性，也是本研究進一步要考慮的問題。

3 結 論

在鐵尾礦砂中添加不同比例的土壤、菌糠，共形成6種重構坡面：松散鐵尾礦砂（LT），結殼鐵尾礦砂（CT），鐵尾礦砂與菌糠混合物（TM），鐵尾礦砂、菌糠和土壤混合物（TSM），鐵尾礦砂和土壤混合物（TS），純土（S），在3 次非均勻間歇模擬降雨條件下，研究重構坡面水動力特性，主要結論如下：

1）鐵尾礦砂重構坡面的產流率隨降雨場次均持續(xù)增大；重構坡面間產流差異不顯著。而重構坡面間的產沙率、坡面微地形差異顯著，因重構坡面類型不同表現出明顯的分異規(guī)律：純尾礦砂、鐵尾礦砂和土壤混合物坡面細溝發(fā)展迅速，坡體上下細溝易貫通，細溝寬深比較小，侵蝕形態(tài)嚴重，為高產沙型坡面，產沙率隨降雨場次持續(xù)增大。覆土坡面和在鐵尾礦砂中添加菌糠坡面可不同程度抑制坡面細溝發(fā)展速度，但覆土坡面細溝僅發(fā)生在坡體中下部局部部位，寬深比較大；添加菌糠的坡面細溝很少，多以均勻分布的菌糠小坑穴出現；這些坡面侵蝕形態(tài)輕微，為低產沙型，產沙率隨降雨場次整體表現為先增加后減少。

2）鐵尾礦砂6 種重構坡面產流率均隨產流歷時持續(xù)波動增加，而產沙率隨產流歷時變化表現出3 種類型：單峰型、增加型和減少型。

3）鐵尾礦砂6 種重構坡面的坡面流流型均為層流，流態(tài)大多為急流；隨降雨場次增加，流速和水流功率整體表現為增大的趨勢，阻力系數、曼寧粗糙系數表現出減小的趨勢，降雨場次對這些水動力因素影響顯著（P<0.05）。6 種坡面因重構基質類型不同，與產流率、產沙率顯著相關的水力學參數和最佳預測函數關系不盡相同，主要為對數、線性或冪函數。

4）鐵尾礦砂坡面覆土減流效果最明顯，添加菌糠的減沙效果最好，但不能有效減流，因此對于土源缺乏的鐵尾礦砂復墾區(qū)，可利用農業(yè)有機廢棄物如菌糠和土壤的無土或者少土模式進行基質改良，有效減緩侵蝕，促進植被修復。