不同復墾模式下鐵尾礦坡面產流產沙與水力特性

呂春娟，郭巖松，畢如田，梁建財，郭星星，陳 丹，王 煜，許彩彩

不同復墾模式下鐵尾礦坡面產流產沙與水力特性

呂春娟，郭巖松，畢如田，梁建財，郭星星，陳 丹，王 煜，許彩彩

（1. 山西農業大學資源環境學院，太谷 030801；2. 山西農業大學農業資源與環境國家級實驗教學示范中心，太谷 030801）

為了揭示降雨條件下鐵尾礦復墾坡面的產流產沙特征與水力特性，該文采用室內人工模擬降雨方法研究了不同降雨強度（60、90、120 mm/h）和不同復墾模式（T，鐵尾礦；TSH，鐵尾礦與高比例生土；TM，鐵尾礦與菌糠；TSM，鐵尾礦、菌糠與生土；TSL，鐵尾礦與低比例生土；植被類型均為胡枝子）25°坡度條件下的侵蝕過程。結果表明：1）產流速率在降雨初期迅速增長，中后期趨于穩定。添加菌糠的復墾模式產流速率高于未添加菌糠的復墾模式；2）產沙速率隨降雨歷時變化呈現兩種模式，添加菌糠的復墾模式為減少型，未添加菌糠的復墾模式為增加型，TSL的減沙效果最為顯著，T對降雨強度的響應最為敏感且高雨強下產沙量最大；3）所有復墾坡面平均徑流深度范圍在0.23～0.93 mm，水流流態均為層流，大部分情況下水流流型為緩流，但TSH在中高雨強下以及T在高雨強下為急流。添加菌糠模式的曼寧糙率系數和水流阻力系數高于未添加菌糠的復墾模式，流速低于未添加菌糠的復墾模式，添加菌糠的基質表面粗糙度較大；4）各復墾坡面的侵蝕形式以濺蝕和片蝕為主，高雨強下T出現細溝侵蝕，其余復墾模式未觀察到明顯細溝侵蝕；5）產流速率與水流功率、雷諾數（2＞0.998）呈極顯著的線性關系，產沙速率與水流功率（2＞0.733）、雷諾數（2＞0.744）呈極顯著的冪函數、指數函數關系。該研究可為鐵尾礦坡面的復墾模式選擇提供理論依據。

復墾；徑流；侵蝕；鐵尾礦；水力特性

0 引 言

鐵礦開采是關系國計民生的重要活動。根據2018年自然資源部發布的《中國礦產資源》記載，鐵礦石產量為12.3億t，2017年鐵礦石消費量為15.7億t，并且每年都呈增長趨勢[1]。由于中國的鐵礦資源多屬貧礦，綜合利用率不高，開采過程中形成大量露天堆棄的鐵尾礦固體廢棄物[2]，其鐵尾礦顆粒松散、容重較大、透氣不良、干板濕黏、養分貧瘠，已成為限制鐵礦復墾地區復墾效果的重要影響因素[3-4]。鐵尾礦產生的重金屬物質隨徑流滲入鄰近農地及河道并污染水系，其尾礦堆積體坡度大且極易引發侵蝕和滑坡等水土流失現象，氣候變暖和極端降雨事件增加又進一步加劇風險，嚴重影響區域生態系統的安全與服務功能[5]。因此，探究鐵尾礦的水土流失特征，對于理解和控制鐵尾礦復墾區水土保持和固體廢棄物污染具有重要意義。

土壤重構和植被重建是礦山復墾措施的重要手段，也是確保礦山復墾方向和實現水土保持效果的關鍵[6]。土壤重構可以在短時間內提供植被生長的土壤條件[7]；植被重建可以通過減緩流速，減少雨滴動能，改善土壤結構的穩定性和滲透性，來達到控制水土流失的目標[8]，還可以通過根莖來吸收和固定礦區土壤中的重金屬從而有效控制土壤污染[9]。因此，礦山復墾措施很大程度上影響了土壤侵蝕過程。例如，Taylor等[10-11]研究了復墾措施對美國肯塔基州東部的露天礦場坡面產流和侵蝕的影響，發現采用松散堆積覆土的方式可以提高滲透率和降低產流產沙的峰值，陡坡條件下也能有效的提高入滲。Luna等[12]通過對復墾土壤基質施加城市污泥、堆肥等改良劑和礫石、木屑等地表覆蓋物來提高基質的有機質含量和增加地表粗糙度，發現施加改良劑能夠提高土壤孔隙率，從而提高滲透和減少水蝕。牛耀彬等[13]通過對不同植被重建的整地措施進行水力學性質的研究，表明防護措施能夠有效減少土壤可蝕性且水平階的防護效果優于魚鱗坑。郭星星等[14]研究了在不同雨強和坡度下鐵尾礦砂裸坡復墾坡面的產流產沙特征，結果發現雨強對產流影響較大，而雨強和坡度對產沙的影響顯著，并且鐵尾礦砂與菌糠混合物的基質減沙效果相比其他類型更好。綜上所述，將土壤重構和植被恢復結合的復墾模式已經成為指導生產實踐的有效舉措。

水蝕是由降雨和徑流等侵蝕力引起土壤物質分離、運輸和沉積的過程[15]。在斜坡上這些過程高度依賴土壤基質的水力特性。一些學者從水力學角度來研究土壤侵蝕過程并進行土壤侵蝕預測[16]。Nearing等[17]認為水流功率是邁阿密粉質壤土或者塞西爾沙壤土的坡面侵蝕的最佳預測指標。康宏亮等[18]發現單位徑流功率是描述風沙區土質和10%礫石質量分數工程堆積體侵蝕產沙的最佳指標，徑流功率是描述20%、30%礫石質量分數工程堆積體土壤侵蝕參數的最優指標。Guo等[19]研究表明在預測黃土坡面土壤的侵蝕量和溶質運移速率時，曼寧糙率系數與平均水流深度的比值是比水流功率或雷諾數更好的預測指標。有效模擬土壤侵蝕過程，需要考慮坡面的具體情況和研究目標，選擇合理的水動力參數。

然而，使用上述水力學參數在降雨條件下對鐵尾礦復墾坡面的研究較為缺乏。為此，本文采用人工模擬降雨的方法，研究不同降雨強度下不同復墾模式鐵尾礦坡面的土壤侵蝕過程，探索鐵尾礦復墾坡面的產流產沙規律和水力特性，并據此建立土壤侵蝕預報模型，以期為提升鐵尾礦復墾效果和坡面水土流失治理提供理論依據。

1 材料與方法

1.1 試驗材料

鐵尾礦研究區位于山西省運城市垣曲縣，介于35°14′14″～35°14′54″N，111°40′24″～111°41′09″E，屬暖溫帶半干旱、半濕潤大陸性季風氣候，如圖1所示。年平均氣溫為13.5 ℃，年平均蒸散量為2090 mm。當地植被除因采礦活動造成植被破壞嚴重外，其余地區植被覆蓋良好，植被以灌木叢和草本植物占優勢。該礦區位于黃河支流——亳清河流域范圍內，小浪底壩濕地自然保護區的上游，年均降雨量為631 mm，降水一般集中在5－9月，短歷時暴雨頻發，在鐵尾礦復墾地區引發水土流失和滑坡事件的風險較大。

圖1 研究區示意圖

供試材料是鐵尾礦，土壤和菌糠。鐵尾礦采集于山西省運城市垣曲縣國泰礦業有限公司泉子溝干排尾礦庫。鐵尾礦帶回實驗室后，放置在陰涼地方自然風干，待風干后過2 mm的篩。菌糠來源于蘑菇種植基地的菌糠棒，菌糠棒主要材料由玉米芯、木屑、麩皮等組成，菌糠棒經過粉碎處理，過5mm篩。供試土壤類型屬于褐土，與采礦區原始生境的土壤類型相似，待風干后過5 mm的篩，去除植被、石塊和雜物后準備作為裝填土槽的土壤材料。參照美國制土壤質地分類標準（黏粒＜0.002 mm，粉粒0.002～<0.05 mm，砂粒0.05～2 mm），粒徑組成為：鐵尾礦的砂粒、粉粒和黏粒質量分數分別為29.7%、57.8%、12.5%，生土的砂粒、粉粒和黏粒質量分數分別為21.56%、57.2%、21.24%，質地均屬于粉壤土。

1.2 試驗裝置

在山西農業大學資源與環境學院的模擬降雨廳對不同復墾模式的坡面進行了模擬降雨試驗。試驗裝置由2個獨立部分組成：模擬降雨設備和土槽。模擬降雨設備（QYJY-501T）由西安清遠西安測控技術有限公司定制，系統通過改變噴嘴尺寸和水壓來調節降雨強度。混凝土土槽規格為：長2 m × 寬1 m × 高（0.6～2.3 m），土壤坡面裝填坡度可變范圍在0°～40°之間，土槽前端設計有集流口，使用徑流桶于集流口處收集徑流泥沙樣。Wang等[20]基于同類規格土槽開展巖石覆蓋對粗骨土土壤流失和水文響應影響的研究。本研究共有3種降雨強度：60、90、120 mm/h，下文均簡稱為低雨強、中雨強和高雨強。這是位于中國山西省南部以及黃土高原東部山地的中條山地區典型的普通暴雨、中型暴雨、大型暴雨等降雨等級[21]，經常被用于黃土質坡地的降雨試驗研究[22-23]，對于鐵尾礦復墾地區來說也是同樣受限制的降雨強度。模擬降雨器經過多次校準后，測得的降雨強度分別為61.9～64.6，88.9～94.5和119.8～122.8 mm/h，模擬精度為±5 mm/h，在允許誤差范圍內，有效降雨范圍內的均勻度超過80%。

1.3 試驗設計與樣品采集

試驗設計5個土槽復墾坡面，分別對應野外試驗基地5種鐵尾礦的土壤和植被恢復組合。5個復墾模式的基質類型如下：鐵尾礦（iron tailings, T），鐵尾礦與高比例土壤摻混物（iron tailings and high proportion of soil, TSH），鐵尾礦與菌糠摻混物（iron tailings and mushroom, TM），鐵尾礦、菌糠與土壤摻混物（iron tailings and mushroom and soil, TMS），鐵尾礦與低比例土壤摻混物（iron tailings and low proportion of soil, TSL）。土槽底部40 cm均為鐵尾礦，表層20 cm為不同復墾基質類型。復墾基質的裝填采取分層法[24]，保持垂直斜面以5 cm深度為間隔，裝填與壓實同步進行，各層間通過抓毛措施以確保容重均勻分布且防止出現分層現象。所有復墾坡面的平均坡度均為25°，為野外鐵尾礦復墾坡面的代表性坡度。在模擬降雨開始前，在土槽斜坡表面不同位置進行取樣以測定基質的理化性質，包括土壤顆粒組成[25]、容重、飽和持水量、田間持水量和土壤有機碳等[26]。具體參數和結果詳見表1。

表1 復墾模式基本概況

注：T指鐵尾礦；TSH指鐵尾礦與高比例土壤摻混物；TM指鐵尾礦與菌糠摻混物；TMS指鐵尾礦、菌糠與土壤摻混物；TSL指鐵尾礦與低比例土壤摻混物，下同。

Note: T means substrate of iron tailings; TSH means blend substrate of iron tailings and high proportion of soil; TM means blend substrate of iron tailings and mushroom; TMS means blend substrate of iron tailings、mushroom and soil; TSL means blend substrate of iron tailings and low proportion of soil, the same as below.

5個復墾模式的植被類型均為胡枝子（），胡枝子的枝葉茂盛和根系發達，可有效地保持水土和減少地表徑流。播種前進行坡面整地，沿等高線方向以15 cm為間隔開挖約2 cm寬的水平階，如圖2所示，其作用主要是提高土壤水分入滲，攔蓄泥沙，以及降低植被生長期間坡面表層土壤侵蝕的程度[27]。降雨廳為半開放降雨廳，能夠滿足胡枝子所需的光照條件，生長期內通過人工澆水保持充足的供水條件。胡枝子出苗1個月，開始定苗，確保每個土槽的胡枝子密度相同，以品字形均勻分布在坡面上。在2018年5月開始裝填土槽播種，2018年6月對坡面植被進行定苗，模擬降雨試驗于2018年9月進行，在模擬降雨試驗后，對每個土槽植株的覆蓋率[28]和生物量進行測定（表1）。

圖2 坡面水平階與植物分布圖

對于同一個土槽，按照降雨強度從小到大的順序逐次降雨。降雨開始前，使用土壤水分測定儀對土槽基質含水率進行測定，通過預降雨和晾曬等方式將前期土壤含水率控制在8%±2%；每次降雨前后使用溫度計測量水溫，以計算運動粘度。降雨開始后，以坡面開始產流的時間為計時起點，以3 min為間隔在土槽前端V形排水口使用水桶收集泥沙樣本，單場降雨產流后歷時30 min，共收集10個泥沙樣本；坡面流速采用高錳酸鉀染料示蹤法進行測量[29]。高錳酸鉀染料示蹤法測量的是最大表面速度，乘以層流流態的校正因子0.67得到平均流速[30]。降雨結束后，使用土壤水分測定儀測定土槽基質的表層0～5 cm土壤含水率；使用量筒測量每個泥沙樣品的體積，待泥沙樣品靜置沉淀澄清后，倒掉上部清液，將泥沙轉移到燒杯中，將泥沙和燒杯放入通風烘箱在105 ℃下持續12 h至恒質量并稱質量[31]。

1.4 數據分析處理

平均徑流深度、運動黏度的計算公式[32]為

式中為平均徑流深度，m；R為產流速率，L/(m2·min)；為土槽的寬度，m；為平均流速，m/s；為實測水溫，℃。為重力加速度，m/s2，通常取9.8 m/s2。

式中為水力半徑，m，可近似用平均水流深度代替；為水力坡度，可用坡度的正弦值近似替代；為水的密度，g/cm3。

使用Excel 2010、Origin 9.1與R 3.4.4軟件進行數據處理與繪圖。

2 結果與分析

2.1 復墾坡面產流特征

圖3為不同復墾模式的產流速率隨降雨強度的變化。隨降雨時間的延長，產流速率的變化趨勢基本相似。在低雨強下，不同處理隨時間變化的變異程度較低，基本呈穩定趨勢。在中雨強和高雨強下，坡面起流的初始階段（0~9 min）其產流速率迅速增長，中后期增長變緩或趨于穩定。表2結果表明，在低降雨強度下，TMS的產流速率最大，TM和TSH次之，T和TSL最小。中雨強和高雨強下，添加菌糠復墾模式產流速率與未添加菌糠的復墾模式差異顯著，TM和TMS的產流速率均大于其余的復墾措施的產流速率。表明在低雨強下，添加較高比例的生土和只添加菌糠的這兩種復墾模式對產流速率變化的影響其差異不顯著，但在中雨強度和高降雨強度下，添加菌糠的復墾模式產流速率更高。比較不同降雨強度間產流速率的增長率變化，未添加菌糠的復墾模式其平均產流速率為0.34 L/(m2·min)，在中降雨強度和高降雨強度下其平均增長率為187.32%和46.7%。添加菌糠的復墾模式其平均產流速率為0.65 L/(m2·min)，在中降雨強度和高降雨強度下其平均增長率為109.26%和25.9%。上述結果表明未添加菌糠的復墾模式在不同降雨強度下產流速率的增幅大于添加菌糠的基質，但是其產流速率總體較低于添加菌糠的復墾模式。

圖3 不同降雨強度不同復墾措施產流速率變化

表2 不同降雨強度不同復墾措施侵蝕和水流特征

注：不同字母代表各處理間的差異顯著（<0.05），下同。

Notes: Different letters indicate significant differences at the 0.05 level, and the same as below.

土壤基質的入滲能力對坡面產流速率的大小會有很大影響。表2結果表明，隨著降雨強度的增加，其降雨后表層含水率也相應遞增。TMS在低雨強條件下降雨后含水率較高，而在中高雨強下降雨后含水率與TM接近。再者，與對應復墾基質的持水能力相比，在單次模擬降雨結束之后，T、TSH和TSL的表層含水率在高降雨強度下均能達到田間持水量與飽和持水量之間，而TM和TMS在不同降雨強度下均未能達到田間持水量的水平（表1和表2）。上述結果表明，未添加菌糠的基質其在降雨過程中入滲能力更強，因而其產流速率相比添加菌糠的基質更低，添加菌糠的基質更容易在坡面表面形成致密的不透水層，從而降低入滲能力，提高坡面的產流速率。

2.2 復墾坡面產沙特征

由圖4可知，不同復墾模式的產沙速率隨時間和降雨強度的變化其特征也具有明顯的差異。隨降雨時間的延長，產沙速率的變化呈現兩種模式，即增加型和減少型。在低雨強下，產沙速率隨時間變化其變異程度較低，穩定性較高，除TSH外其余復墾模式的變化趨勢整體趨穩。在中雨強和高雨強下，T、TSH和TSL隨時間延長其產沙速率增加，TM和TMS隨時間延長其產沙速率降低。不同降雨強度下，TSH與TSL的產沙速率具有明顯差異且TSH大于TSL，TM和TMS之間的產沙速率差異不大。雖然隨降雨強度的增加復墾坡面的產沙速率增加，但不同復墾模式對其的響應程度具有差異特征。T在低降雨強度其平均產沙速率為0.56 g/(m2·min)，在中雨強和高雨強下其平均增長率為2690.71%和164.88%，隨降雨強度的變化其增長率變化幅度最大。TM和TMS的產沙速率在低中降雨強度的平均增長率為907.41%且增幅較高，在中高雨強的平均增長率為26.17%且增幅較小。TSH和TSL在低中雨強的平均增長率為332.05%，與其他復墾模式相比較低但增幅較高，在中高雨強的平均增長率為112.2%且增幅趨緩。上述結果表明，純尾礦坡面雖然在低降雨強度下其侵蝕程度相較其他復墾模式更低，但是對降雨強度的響應更為靈敏，尤其是在中高降雨強度下需要注意；添加菌糠的復墾模式坡面產沙規律類似，與其余復墾模式相比無論從時間響應特征還是侵蝕程度均有明顯差異；添加較低比例的砂土混合模式TSL其減沙效果最為明顯，隨降雨強度的增加其變化程度最低，在高雨強下其產沙速率也能維持在較低水平；另外，在TSH和TSL的植被覆蓋率和生物量相比其他復墾模式而言更為茂盛，但在不同降雨強度下TSL的產沙速率也比TSH要低，說明復墾基質的作用還是較為明顯。

圖4 不同降雨強度不同復墾模式產沙速率

2.3 復墾坡面水力特性

通過水流流速、水流深度、雷諾數和弗汝德數來反映坡面流的特征。由圖5可知，在低雨強下隨降雨時間的延長除TSH外的不同復墾模式的坡面流速均表現平緩的趨勢且變異程度較低。在中雨強下，TSH的坡面流速隨降雨時間的延長其變化趨勢為先增加后穩定且變異程度較高，但在高降雨強度下TSH的變化趨勢整體平穩，其時間序列結果與產沙速率的相關性很強，低雨強流速的波動也導致了產沙速率的波動。值得注意的是，T的坡面流速在高雨強下隨降雨時間的延長其中后期（15～30 min）會發生較大的突變增長且增幅為93.9%，這說明在降雨的后期坡面發生較大程度的崩解和侵蝕從而造成流速的突變增長。由表2可知，隨著降雨強度增加其水流流速增加，TSH均與其他復墾模式表現出明顯的差異且坡面流速的值最高，T的坡面流速在高降雨強度下會與TM、TMS和TMS表現出明顯的差異且增幅最大。在水流深度方面，隨著降雨強度的增加其水流深度增加，范圍在0.23～0.93 mm之間，水流深度較淺，說明復墾坡面植被恢復整地措施諸如水平階的作用增加了水流彎曲度，從而降低了水流深度。根據坡面淺層均勻流[34]的標準，Re和的臨界值分別為900和1。由表3可知，在不同降雨強度不同復墾模式的Re值在4.23～27.36之間，當Re＜900時則表示水流流態為層流，因此本研究的水流流態均為層流。姚文藝[34]認為，盡管雨滴已對水流有所擾動，水質點有局部摻混現象，但由于雷諾數和阻力系數的關系仍符合在通常明渠層流內的變化規律，表明此時水流的黏滯力的作用仍大于慣性力的作用，整體水流仍處于層流狀態，據此提出“偽層流”的概念。弗勞德數是用來判斷坡面水流是急流還是緩流的參數，＞1，說明重力小于慣性力作用，水流流型為急流；反之，＜1則水流流型為緩流。在本研究中大部分情況下不同降雨強度不同復墾模式的＜1，但TSH在中雨強和高雨強下以及T在高雨強下的＞1，則表現為急流的特征。

一般情況下，曼寧糙率系數隨降雨強度的增加而減少，在復墾模式T、TSH和TSL均觀察到此類現象，但是TMS隨降雨強度的增加而增加，TM為先增后減，表3結果也表明TM和TMS的曼寧糙率系數高于其他復墾模式。阻力系數反映了坡面流在流動過程中所受到的阻力大小，阻力系數越大，說明水流克服坡面阻力所消耗的能量就越大，則用于坡面侵蝕和泥沙輸移的能量就越小，坡面產沙就越少。在低雨強下，TSL的阻力系數較高且變異程度較大，其余復墾模式的阻力系數之間不具有明顯差異。在中雨強下，不同復墾模式的阻力系數不具有明顯的差異且變異程度較大，但在高雨強下阻力系數的排序為：TMS＞TM&TSL＞T&TSH。在低雨強下，不同復墾模式水流剪切力的排序為：TMS＞TM&TSL＞T&TSH，不同復墾模式水流功率的排序為：TMS＞TM&TSH＞T&TSL；在中雨強和高雨強下，TMS和TM的水流剪切力和水流功率的值與其余復墾模式的差異較大且排序靠前，這說明在不同降雨強度下添加菌糠的復墾基質均與未添加菌糠的復墾基質的水流剪切力和水流功率具有明顯的差異且排序靠前。

圖5 不同降雨強度不同復墾模式水流流速

表3 不同降雨強度不同復墾模式的水力參數

2.4 復墾坡面侵蝕特征與水力特性的擬合

表4為不同降雨強度下侵蝕特征與水力參數之間的相關關系。產流方面，不同降雨條件下產流速率與水力參數的相關性程度有差異，總體來看，水流功率、雷諾數和水流剪切力與產流速率的相關性較高，因此將這三個參數作為產流速率的擬合指標。產沙方面，在低中雨強下，雷諾數、水流功率與產沙速率呈極顯著相關性，但在高雨強條件下，水流流速、弗汝德數與產沙速率呈極顯著相關性，雷諾數、水流功率和產沙速率的相關性不顯著，水流剪切力在中雨強下與產沙速率的相關性顯著。因此初步將水流功率、雷諾數和水流剪切力作為產沙速率的擬合指標并進行探索。

在確定擬合指標的基礎上，進一步分析每個擬合指標在不同復墾模式下的產流產沙擬合效果。表5為侵蝕特征與水力參數的擬合關系。擬合結果表明，不同復墾模式下水流功率和雷諾數對產流速率的線性擬合均呈極顯著（＜0.01），而水流剪切力和產流速率除TSH外均呈極顯著（＜0.01）的二次多項式關系，TSH的擬合效果較差，其水流剪切力和產流速率呈顯著（＜0.05）的二次多項式關系。不同復墾模式下產沙速率與水流功率、雷諾數、水流剪切力呈冪函數和指數函數關系，產沙速率和雷諾數、水流功率之間的指數函數、冪函數關系均呈極顯著（＜0.01），而水流剪切力與產沙速率的擬合效果較差，其中T和TSH呈顯著（＜0.05）的冪函數關系，其余復墾模式則均呈極顯著（＜0.01）的冪函數關系。綜上所述，水流功率和雷諾數均可以作為產流速率和產沙速率的水力指標，而水流剪切力的效果相對較差。

表4 侵蝕特征與水力參數的相關分析

注：為水流流速；為雷諾數；為弗汝德數；為阻力系數；為水流剪切力；為水流功率；為曼寧糙率系數；為曼寧糙率系數除以水流深度。*表示顯著相關，**表示極顯著相關。

Note:is the flow velocity;is the Reynolds number;is the Froude Number;is the Darcy-Weisbach friction ;is the hydraulic shear stress;is the stream power;is the Manning roughness coefficient;is the Manning roughness coefficient divided by water flow depth. * is significant correlation at 0.05 level, ** is significant correlation at 0.01 level.

表5 侵蝕特征與水力參數的關系

3 討 論

本文比較了鐵尾礦5種不同復墾模式的產流產沙特征和水力特性，其結果表明復墾模式間均表現明顯的差異特征。在對比產流速率、產沙速率隨時間變化的結果表明，添加和未添加菌糠的復墾模式之間的差異最大。產流速率產生差異的原因在于，未添加菌糠復墾模式基質其水分滲透性更強，其表土層的抗侵蝕能力和團聚體的穩定性能夠有效防止土壤顆粒分離。添加菌糠復墾模式基質其水分滲透性較差，其產沙量也高于未添加菌糠的復墾基質，這可能和侵蝕過程產生的泥沙顆粒組成和復墾基質的質地顆粒粒徑有關。Thomaz等[35]認為，降雨分離的土壤顆粒，其中細顆粒堵塞土壤表面的孔隙，導致徑流的增加和泥沙的運移，降低了土壤的滲透性，表土的粘結性對坡面徑流具有重要影響。產沙速率隨時間序列的延長有兩種模式：增加型和減少型。對于增加型[36]的模式，是在降雨早期，侵蝕主要由雨滴濺蝕引起，土壤顆粒脫落從而增加侵蝕，隨著降雨持續，雨滴沖擊使得土壤水分增加和懸浮顆粒沉積在土壤表面，逐漸導致在粘壤質土壤表面形成限制性密封層，密封層具有提高土壤抗剪強度的作用[30]，從而降低土壤可蝕性，進而導致產沙速率隨累積降雨量的增加而減少，徑流的穩定進一步穩定產沙速率。對于減少型模式，其產沙速率隨時間變化的曲線從現象來說與Zhao等[36]的研究類似，Zhao等的研究表明草或植被根系的存在對沙質土壤的土壤侵蝕過程有影響，降雨開始時本身就有大量原始的松散物質沿斜坡水流運移，從而迅速增加產沙量。然而，本研究中添加菌糠的復墾模式其植被長勢較差，遠低于砂土混合的復墾模式，砂土混合的復墾模式其侵蝕特征為逐漸增長模式，這說明植被措施對土壤侵蝕過程的影響較小。坡面的水力特性結果表明，含有菌糠復墾模式曼寧糙率系數和水流阻力系數高于未添加菌糠復墾模式，流速較低于未添加菌糠的復墾模式，Ding等[37]的研究表明，粗糙的土壤表面與光滑的土壤表面相比在降雨初期其產沙量迅速增加，在20 min左右到達峰值后降低為恒定水平，表面粗糙度與植被措施改良坡面土壤的相比其結論更為貼切，因此本研究中逐漸降低的模式是降雨初期較高的表面粗糙度更易受到降雨和坡面流的沖刷而導致產沙量的增加，在降雨后期侵蝕細顆粒堵塞土壤孔隙，對含菌糠基質的坡面其產流和產沙均有重要影響。綜上所述，斜坡表面粗糙度是影響復墾坡面土壤侵蝕特征差異的重要因素。

坡面侵蝕均可廣泛觀察到濺蝕、片蝕和細溝侵蝕等侵蝕形式，這些侵蝕形式會單獨或同時地發生[38]。Zhang等[39]認為，陡坡的侵蝕形式在上坡（0～80 cm）主要是濺蝕，雨滴濺蝕[40]在土壤表面水分接近飽和時達到最大值，之后隨積水深度增加而減少；中坡（80～160 cm）為濺蝕與片蝕的過渡帶，下坡（160～200 cm）由片蝕主導，而劃分的依據主要是坡面水流深度和流速在斜坡位置的分布。Parsons等[41]認為，水深小于1 mm的條件下，無法有效降低或消除雨滴濺蝕的影響。植被覆蓋率也對濺蝕過程產生影響。本研究中在不同雨強下所有坡面的平均水流深度范圍在0.23～0.93 mm，水深較淺，極易引發雨滴濺蝕。再者，由于本研究的植被恢復措施中采取水平階整地措施，通過增加水流彎曲度來對水流流速產生阻延作用。在試驗過程中，高雨強下純尾礦基質出現細溝侵蝕但程度較輕，形成一條1.4 m長，2 cm寬的細溝，連通性較好，其余復墾模式未觀察到明顯細溝侵蝕。在這些條件的基礎上，本研究推斷坡面的侵蝕形式受到植被分布的影響較大，沒有植被覆蓋的地方，上坡由濺蝕主導，中下坡受到濺蝕和片蝕的共同作用；有植被覆蓋的地方，主要由片蝕主導。值得注意的是，純尾礦基質在中高雨強條件下降雨后期會發生較大程度的崩解和細溝侵蝕，如果這種情況發生在野外，隨著降雨時間的持續發生嚴重土壤侵蝕乃至滑坡的風險會陡然增加。

由前所述，侵蝕特征與水力參數的擬合結果表明，水流功率和雷諾數可以作為產流速率與產沙速率的預測指標。水流剪切力與水流功率的差異主要來源于水流流速的參與[31]，因此在本研究中，流速導致了水流功率成為比水流剪切力更好的預測指標。沙玉清[42]認為，無黏性沙隨沙坡起伏前進，受到剪切力影響較大，而黏性土壤顆粒在起動時會在坡面形成緩和波形，波頂土壤顆粒隨水流而動。本研究復墾基質的質地均為粉壤土，坡面土壤顆粒起動受到的剪切力的影響小，因此水力剪切力在土壤坡面侵蝕過程中影響并不顯著，水流功率和雷諾數等基于水流流速的指標更能反映粉壤質坡面的侵蝕規律。研究結果與高晨燁等[43]、She等[30]和Nearing等[17]等的結論一致。由于本中所采取的植被恢復措施時間較短，坡面土壤性質受到植被措施改良效果的影響較小，本研究可以為鐵尾礦復墾初期不同復墾模式坡面土壤侵蝕特征提供理論依據。

4 結 論

本文研究了不同雨強和復墾模式下鐵尾礦坡地產流產沙特征和水力特性，結果表明：

1）產流規律在降雨初期迅速增長，中后期趨于穩定。添加菌糠的復墾模式其產流速率高于未添加菌糠的復墾模式。產沙規律呈現兩種模式，添加菌糠的復墾模式為減少型，未添加菌糠的復墾模式為增加型，鐵尾礦與低比例土壤摻混物的減沙效果最為顯著，純鐵尾礦對降雨強度的響應最為敏感且高雨強下產沙量最大。

2）各復墾模式下，水流流態均為層流，大部分情況下水流流型為緩流，但鐵尾礦與高比例土壤摻混物在中高雨強下以及純鐵尾礦在高雨強下為急流。不同雨強下所有復墾坡面的平均水流深度范圍在0.23～0.93 mm，水深較淺。添加菌糠復墾模式曼寧糙率系數和水流阻力系數高于未添加菌糠的復墾模式，流速低于未添加菌糠的復墾模式。

3）各復墾坡面的侵蝕形式以濺蝕和片蝕為主。侵蝕形式受到植被分布的影響較大，沒有植被覆蓋的地方，上坡由濺蝕主導，中下坡受到濺蝕和片蝕的共同作用；有植被覆蓋的地方，主要由片蝕主導。高雨強下純鐵尾礦出現細溝侵蝕但程度較輕，其余復墾模式未觀察到明顯細溝侵蝕。

4）水流功率和雷諾數均可以作為產流速率和產沙速率的水力指標。不同復墾模式下產流速率與水流功率、雷諾數呈極顯著的線性關系，不同復墾模式下產沙速率與水流功率、雷諾數呈極顯著的冪函數、指數函數關系。

[1]Sun Sizhong, Anwar Sajid. R&D activities and FDI in China’s iron ore mining industry[J]. Economic Analysis and Policy, 2019, 62: 47－56.

[2]鄧文，江登榜，楊波，等. 我國鐵尾礦綜合利用現狀和存在的問題[J]. 現代礦業，2012，27(9)：1－3.

Deng Wen, Jiang Dengbang, Yang Bo, et.al. Comprehensive utilization status and existing problems of iron tailings in China[J]. Morden Mining, 2012, 27(9): 1－3. (in Chinese with English abstract)

[3]陳穎，李晨光，徐學華，等. 鐵尾礦廢棄地火炬樹林水土保持功能[J]. 水土保持學報，2015，29(3)：106－111.

Chen Ying, Li Chenguang, Xu Xuehua, et al. Soil and water conservation function of. forest in iron tailings[J]. Journal of Soil and Water Conservation, 2015, 29(3): 106－111. (in Chinese with English abstract)

[4]呂春娟，畢如田，陳衛國，等. 土壤結構調理劑PAM對復墾鐵尾礦砂物理性狀的影響[J]. 農業工程學報，2017，33(6)：240－245.

Lü Chunjuan, Bi Rutian, Chen Weiguo, et al. Effect of soil structure conditioner PAM on physical properties of iron tailings in reclaiming[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(6): 240－245. (in Chinese with English abstract)

[5]Feng Yu, Wang Jinman, Bai Zhongke, et al. Effects of surface coal mining and land reclamation on soil properties: A review[J]. Earth-Science Reviews, 2019, 62: 47－56.

[6]方星. 礦山土地復墾理論與方法[M]. 北京：地質出版社，2015：318－334.

[7]Shrestha Raj K, Lal Rattan. Carbon and nitrogen pools in reclaimed land under forest and pasture ecosystems in Ohio, USA[J]. Geoderma, 2010, 157(3/4): 196－205.

[8]Feng Tianjiao, Wei Wei, Chen Liding, et al. Assessment of the impact of different vegetation patterns on soil erosion processes on semiarid loess slopes[J]. Earth Surface Processes & Landforms, 2018, 43 (9): 1860－1870.

[9]Mahar Amanullah, Wang Ping, Ali Amjad, et al. Challenges and opportunities in the phytoremediation of heavy metals contaminated soils: A review[J]. Ecotoxicology and Environmental Safety, 2016, 126: 111－121.

[10]Taylor Timothy J, Agouridis Carmen T, Warner Richard C, et al. Runoff curve numbers for loose-dumped spoil in the Cumberland Plateau of eastern Kentucky[J]. International Journal of Surface Mining, Reclamation and Environment, 2009, 23(2): 103－120.

[11]Hoomehr Siavash, Schwartz John S, Yoder Daniel C, et al. Curve numbers for low-compaction steep-sloped reclaimed mine lands in the Southern Appalachians[J]. Journal of Hydrologic Engineering, 2013, 18(12): 1627－1638.

[12]Luna Lourdes, Vignozzi Nadia, Miralles Isabel, et al. Organic amendments and mulches modify soil porosity and infiltration in semiarid mine soils[J]. Land Degradation & Development, 2017, 29(4): 1019－1030.

[13]牛耀彬，高照良，劉子壯，等. 工程措施條件下堆積體坡面土壤侵蝕水動力學特性[J]. 中國水土保持科學，2015，13(6)：105－111.

Niu Yaobin, Gao Zhaoliang, Liu Zizhuang, et al. Hydrodynamic characteristics of soil erosion on deposit slope under engineering measures[J]. Science of Soil and Water Conservation, 2015, 13(6): 105－111. (in Chinese with English abstract)

[14]郭星星，呂春娟，陳丹，等. 雨強和坡度對鐵尾礦砂坡面復墾前后產流產沙的影響[J]. 水土保持研究，2019，26(1)：8－13.

Guo Xingxing, Lü Chunjuan, Chen Dan, et al. Effects of rainfall intensity and slope gradient on runoff and sediment production on the slope before and after recalmation in iron tailings[J]. Research of Soil and Water Conservation, 2019, 26(1): 8－13. (in Chinese with English abstract)

[15]Shi Zhihua, Yan Fengling, Li Lu, et al. Interrill erosion from disturbed and undisturbed samples in relation to topsoil aggregate stability in red soils from subtropical China[J]. Catena, 2010, 81(3): 240－248.

[16]Zhang Yinghu, Zhang Mingxiang, Niu Jianzhi, et al. Rock fragments and soil hydrological processes: Significance and progress[J]. Catena, 2016, 147: 153－166.

[17]Nearing M A, Norton L D, Bulgakov D A, et al. Hydraulics and erosion in eroding rills[J]. Water Resources Research, 1997, 33(4): 865－876.

[18]康宏亮，王文龍，薛智德，等. 北方風沙區礫石對堆積體坡面徑流及侵蝕特征的影響[J]. 農業工程學報，2016，32(3)：125－134.

Kang Hongliang, Wang Wenlong, Xue Zhide, et al. Effect of gravel on runoff and erosion characteristics on engineering accumulation slope in windy and sandy area, northern China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(3): 125－134. (in Chinese with English abstract)

[19]Guo Tailong, Wang Quanjiu, Li Dingqiang, et al. Flow hydraulic characteristic effect on sediment and solute transport on slope erosion[J]. Catena, 2013, 107: 145－153.

[20]Wang Xiaoyan, Li Zhaoxia, Cai Chongfa, et al. Effects of rock fragment cover on hydrological response and soil loss from Regosols in a semi-humid environment in South-West China[J]. Geomorphology, 2012(151/152): 234－242.

[21]張漢雄. 黃土高原的暴雨特性及其分布規律[J]. 地理學報，1983，38(4)：416－425.

Zhang Hanxiong. The characteristics of hard rain and its distribution over the loess plateau[J]. Acta Geographica Sinica, 1983, 38(4): 416－425. (in Chinese with English abstract)

[22]Han Luyan, Jiao Juying, Jia Yanfeng, et al. Seed removal on loess slopes in relation to runoff and sediment yield[J]. Catena, 2011, 85(1): 12－21.

[23]Shi Z H, Fang N F, Wu F Z, et al. Soil erosion processes and sediment sorting associated with transport mechanisms on steep slopes[J]. Journal of Hydrology, 2012(454/455): 123－130.

[24]倪世民，馮舒悅，王軍光，等. 不同質地重塑土坡面細溝侵蝕形態與水力特性及產沙的關系[J]. 農業工程學報，2018，34(15)：149－156.

Ni Shimin, Feng Shuyue, Wang Junguang, et al. Relationship between rill erosion morphology and hydraulic characteristics and sediment yield on artificial soils slope with different textures[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(15): 149－156. (in Chinese with English abstract)

[25]Bieganowski A, Ryzak M, Witkowska-Walczak B. Determination of soil aggregate disintegration dynamics using laser diffraction[J]. Clay Minerals, 2010, 45(1): 23－34.

[26]郭巖松，呂春娟，郭星星，等. 不同復墾措施下鐵尾礦的土壤顆粒多重分形特征[J]. 土壤通報，2019，50(1)：81－88.

Guo Yansong, Lü Chunjuan, Guo Xingxing, et al. Multifractal characteristics of soil particles in iron tailings under different reclamation measures[J]. Chinese Journal of Soil Science, 2019, 50(1): 81－88. (in Chinese with English abstract)

[27]孫秀殿，李純麗，張鳳霞. 胡枝子的栽培利用[J]. 特種經濟動植物，1999(2)：33.

[28]Guijarro M, Pajares G, Riomoros I, et al. Automatic segmentation of relevant textures in agricultural images[J]. Computers & Electronics in Agriculture, 2011, 75(1): 75－83.

[29]Abrantes Joao R C B, Moruzzi Rodrigo B, Silveira Alexandre, et al. Comparison of thermal, salt and dye tracing to estimate shallow flow velocities: Novel triple-tracer approach[J]. Journal of Hydrology, 2018, 557: 362－377.

[30]She Dongli, Liu Dongdong, Meng Jiajia, et al. Hydrological responses on saline-sodic soil slopes in a coastal reclamation area of China[J]. Catena, 2016, 144: 130－140.

[31]Wei Wei, Chen Liding, Fu Bojie, et al. Responses of water erosion to rainfall extremes and vegetation types in a loess semiarid hilly area, NW China[J]. Hydrological Processes, 2010, 23(12): 1780－1791.

[32]Zhang Letao, Gao Zhaoliang, Yang Shiwei, et al. Dynamic processes of soil erosion by runoff on engineered landforms derived from expressway construction: A case study of typical steep spoil heap[J]. Catena, 2015, 128: 108－121.

[33]李少華，何丙輝，李天陽，等. 紫色土丘陵區生物埂不同植被類型土壤分離水動力學特征[J]. 水土保持學報，2019，33(3)：70－75.

Li Shaohua, He Binghui, Li Tianyang, et al. Hydraulic characteristics of soil detachment under the terrace banks with different vegetation types in the purple hilly area[J]. Journal of Soil and Water Conservation, 2019, 33(3): 70－75. (in Chinese with English abstract)

[34]姚文藝. 坡面流阻力規律試驗研究[J]. 泥沙研究，1996(1)：74－82.

Yao Wenyi. Experiment study on hydraulic resistance laws of overland sheet flow[J]. Journal of Sediment Research, 1996(1): 74－82. (in Chinese with English abstract)

[35]Thomaz Edivaldo L, Pereira Adalberto A. Misrepresentation of hydro-erosional processes in rainfall simulations using disturbed soil samples[J]. Geomorphology, 2017, 286: 27－33.

[36]Zhao Chunhong, Gao Jianen, Huang Yuefei, et al. The contribution of astragalus adsurgens roots and canopy to water erosion control in the water-wind crisscrossed erosion region of the loess plateau, China[J]. Land Degradation & Development, 2016, 28(1): 265－273.

[37]Ding Wenfeng, Huang Chihua. Effects of soil surface roughness on interrill erosion processes and sediment particle size distribution[J]. Geomorphology, 2017, 295: 801－810.

[38]Wang Zhijie, Jiao Juying, Rayburg Scott, et al. Soil erosion resistance of “Grain for Green” vegetation types under extreme rainfall conditions on the Loess Plateau, China[J]. Catena, 2016, 141: 109－116.

[39]Zhang X C, Wang Z L. Interrill soil erosion processes on steep slopes[J]. Journal of Hydrology, 2017, 548: 652－664.

[40]Poesen J. Rainwash experiments on erodibility of loose sediments[J]. Earth Surface Processes & Landforms, 2010, 6(3/4): 285－307.

[41]Parsons Anthony J, Abrahams Athol D, Wainwright Jhon. Rainsplash and erosion rates in an interrill area on semi-arid grassland, Southern Arizona[J]. Catena, 1994, 22(3): 215－226.

[42]沙玉清. 泥沙運動學引論[M]. 北京：中國工業出版社，1965.

[43]高晨燁，張寬地，楊明義. 基于無量綱水流強度指標的坡面流輸沙能力計算方法[J]. 農業工程學報，2018，34(17)：134－142.

Gao Chenye, Zhang Kuandi, Yang Mingyi. Overland flow sediment transport capacity calculation method based on non-dimensional flow intensity index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(17): 134－142. (in Chinese with English abstract)

Effects of different reclamation patterns on surface runoff, sediment yield and hydraulic characteristics of slopes in iron ore tailings

Lü Chunjuan, Guo Yansong, Bi Rutian, Liang Jiancai, Guo Xingxing, Chen Dan, Wang Yu, Xu Caicai

(1,030801,;2,,030801,)

Reclamation is a key measure to ameliorate soil erosion from slopes in mining tailings. However, how rainfall intensity and different reclamation patterns combine to affect soil erosion remains elusive, especially in the loess plateau in China. We experimentally investigated these using artificial slopes in iron ore tailings as an example. The experiment was conducted in laboratory using a rainfall simulator. We compared three rainfall intensities, 60, 90 and120 mm/h, and five remediation methods: iron ore tailings without adding any amendment (T), overlaying the iron ore tailings with a high proportion of soil (TSH), planting mushroom in the iron ore tailings (TM), overlaying soil on the iron ore tailings and planting mushroom in the soil (TSM); overlaying a low proportion of soil on the iron ore tailings (TST). The experiment for each treatment was conducted in a 2 m × 1 m × 2.3 m soil flume, in which a layer of iron ore tailings 40 cm thick was underlaid by a 20 cm of above remediating substrate. The cultivar of the vegetation was. The angle between the soil flume to the horizontal was 25o. The soil water content in all treatments prior to rainfall was the same. During the experiment, we measured the water flow velocity and collected sediment samples after each three minutes to calculate the surface runoff and sediment erosion rate. The hydraulic parameters were also calculated in attempts to determine its impact on sediment erosion. The results showed: 1) The surface runoff rate increased asymptotically with time and it was higher in treatments with mushroom than those without mushroom. 2) The mushroom had a significant impact on sediment erosion. Regardless of rainfall density, the sediment erosion rate decreased with time in treatments with mushroom while increased with time in treatments without mushroom, especially in TSL. Soil erosion in T was most sensitive to rainfall intensity and it increased with rainfall intensity. 3) Water flow was laminar in most cases. Except in TSH under moderate and high rainfall intensity and T under high rainfall intensity where water flow was supercritical, water flow in other treatments was subcritical. The average infiltration depth in all treatments was shallow and ranged from 0.23 to 0.93 mm. The Manning roughness coefficient and the Darcy-Weisbach friction of the remediating substrates planted with mushroom were higher than that without mushroom. Mushroom reduced water velocity but increased the surface roughness compared to those without mushroom. 4) Splash erosion and shear erosion was the main cause of sediment erosion. 5) The runoff rate was linearly related to the stream power and to the Reynolds number with2>0.998. The erosion rate was related to the steam power exponentially with2>0.733 and to the Reynolds number in a power-law with2>0.744. This study improves our understanding of how different mediating methods and rainfall density combine to impact surface runoff and sediment erosion in slopes in iron or tailings.

reclamation; runoff; erosion; iron tailings; hydraulic characteristics

呂春娟，郭巖松，畢如田，梁建財，郭星星，陳 丹，王 煜，許彩彩. 不同復墾模式下鐵尾礦坡面產流產沙與水力特性[J]. 農業工程學報，2020，36(2)：156－165.doi：10.11975/j.issn.1002-6819.2020.02.019 http://www.tcsae.org

Lü Chunjuan, Guo Yansong, Bi Rutian, Liang Jiancai, Guo Xingxing, Chen Dan, Wang Yu, Xu Caicai. Effects of different reclamation patterns on surface runoff, sediment yield and hydraulic characteristics of slopes in iron ore tailings[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(2): 156－165. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.02.019 http://www.tcsae.org

2019-07-26

2019-12-28

國家青年科學基金（41401619）；山西農業大學引進人才科研啟動項目（2014ZZ07）

呂春娟，博士，副教授，主要從事土壤侵蝕與土地整治的研究。Email：lcjcwg@126.com

10.11975/j.issn.1002-6819.2020.02.019

S157.1

A

1002-6819(2020)-02-0156-10