鄭騰輝，邢媛媛，何凱旋，滕　飛，李光錄,

(1 西北農(nóng)林科技大學(xué) a水土保持研究所,b資源環(huán)境學(xué)院，陜西 楊凌 712100；2 中國科學(xué)院 新疆生態(tài)與地理研究所，新疆 烏魯木齊 830011)

?

雨滴擊濺與薄層水流混合侵蝕的輸沙機理

鄭騰輝1a，邢媛媛1b,2，何凱旋1b，滕飛1b，李光錄1a,1b

(1 西北農(nóng)林科技大學(xué) a水土保持研究所,b資源環(huán)境學(xué)院，陜西 楊凌 712100；2 中國科學(xué)院 新疆生態(tài)與地理研究所，新疆 烏魯木齊 830011)

[摘要]【目的】 針對雨滴擊濺與薄層水流混合侵蝕過程輸沙機理研究中存在的問題,研究雨滴擊濺對坡面徑流輸沙的影響?！痉椒ā?采用室內(nèi)人工模擬降雨試驗,以土婁土、黃綿土和黑壚土為供試土壤，在不同雨滴直徑(直徑2.22，2.68和3.04 mm)、不同坡度(0°，2°和4°)和薄層水流厚度(0，2，4和6 mm)條件下，分析單雨滴擊濺時，泥沙濺蝕量與水層厚度的關(guān)系；多雨滴擊濺時，泥沙濺蝕量與薄層水流厚度的關(guān)系，分別對二者關(guān)系進行擬合；同時，對于多雨滴擊濺時，不同薄層水流厚度下黃綿土泥沙濺蝕量和雨滴動能的關(guān)系進行了分析?！窘Y(jié)果】 單雨滴擊濺時，泥沙濺蝕量隨水層厚度的增加而減少，且雨滴直徑越大其擾動水層厚度越大，當(dāng)雨滴直徑由2.22 mm增加到3.04 mm時，平均擾動水層厚度由10 mm增加到14.67 mm。多雨滴擊濺時，在相同的坡度(0°～4°)，不同土壤和雨滴直徑下，隨薄層水流厚度的增加，泥沙濺蝕量呈先升后降趨勢；當(dāng)薄層水流厚度相同時，隨雨滴直徑的增加，泥沙濺蝕量呈增大趨勢。當(dāng)雨滴動能不變時，泥沙濺蝕量隨著薄層水流厚度的增加而增大；當(dāng)薄層水流厚度一定時，泥沙濺蝕量隨雨滴動能的增加而增大?！窘Y(jié)論】 建立了單雨滴擊濺和多雨滴擊濺條件下擾動水層厚度關(guān)與泥沙濺蝕量的關(guān)系式。

[關(guān)鍵詞]模擬降雨；薄層水流；混合侵蝕；雨滴動能；雨滴擊濺；輸沙機理

雨滴濺蝕是以雨滴打擊為主要動力，打擊地表引起土壤分散、躍移的一種侵蝕方式[1]。薄層徑流則是通過對土壤分散、懸浮和搬運而形成侵蝕[2]。雨滴擊濺與薄層徑流混合作用于某一坡面時,其侵蝕過程更為復(fù)雜,且侵蝕強度明顯增加。雨滴侵蝕量一般占總侵蝕量的70%以上，薄層水流侵蝕所產(chǎn)生的泥沙絕大部分是由雨滴擊濺所致[3]。雨滴從空中垂直下落，對坡面上的泥沙產(chǎn)生2種作用，一種在垂直坡面方向，雨滴對泥沙顆粒產(chǎn)生向下的打擊作用，使坡面泥沙更加密實，更難于起動;另一種沿坡面水流方向，雨滴對坡面泥沙顆粒產(chǎn)生分散作用，使坡面泥沙松動，易于起動[4]。黃土坡面土壤侵蝕發(fā)生發(fā)展過程包括細溝間侵蝕過程、細溝侵蝕過程和淺溝侵蝕過程,三者的發(fā)生發(fā)展及其在坡面侵蝕產(chǎn)沙中的作用是目前坡面侵蝕機制研究的核心內(nèi)容,也是坡面侵蝕預(yù)報模型所考慮的關(guān)鍵所在[5]。野外和室內(nèi)試驗也證實,雨滴擊濺對侵蝕和泥沙輸移具有重要的影響[6-11]。坡面薄層水流并不是簡單的遵循雷諾規(guī)律，它應(yīng)該既不是層流，也不是充分紊流，而是一種介于它們之間的特殊流動[12]。雨滴打擊能夠使坡面薄層水流由層流變?yōu)槲闪骰驍_動層流,從而增加其侵蝕能力[6-8]。 Kinnel[13-17]對非黏性沙土的雨滴擊濺試驗結(jié)果顯示, 隨薄層水流深度的增加, 泥沙剝蝕量減??；Gabet等[18]研究認為,坡面薄層水流與濺蝕關(guān)系明顯,對侵蝕產(chǎn)沙有較大制約作用。坡面是構(gòu)成山地丘陵和破碎高原的重要景觀單元,而坡面薄層水流產(chǎn)生的拖拽力和雨滴擊濺力是導(dǎo)致坡面土樣侵蝕的主要營力,雨滴擊濺擾動作用下坡面薄層水流的泥沙輸移機理十分復(fù)雜,一直是人們研究的難點和熱點問題[19-20]。為此，本研究采用人工模擬降雨試驗方法，選用3種不同類型的土壤，測定了單雨滴擊濺時不同雨滴直徑與水層厚度條件下不同土壤類型的泥沙濺蝕量，以及多雨滴擊濺時不同雨滴直徑與薄層水流厚度條件下不同土壤類型的泥沙濺蝕量和不同雨滴動能與薄層水流厚度下黃綿土的泥沙濺蝕量，以期為揭示不同條件下雨滴擊濺與薄層水流混合侵蝕的過程及其機理提供參考。

1材料與方法

1.1試驗裝置

1.1.1雨滴發(fā)生裝置雨滴發(fā)生裝置如圖1-a所示，其整體為長方形的槽形裝置，上部開口，長70 cm，寬20 cm，高25 cm。其左邊是一個進水裝置，試驗時運用虹吸原理向長槽里注水，降雨時可通過調(diào)節(jié)進水管高度來控制水槽中的水頭高度；在長槽底部按1 cm等間距設(shè)置直徑1 cm的圓孔，圓孔用橡膠塞封閉，每個橡膠塞可插一個針頭。通過不同針號以產(chǎn)生不同直徑的單雨滴和多雨滴。為防止在試驗過程中，水中雜質(zhì)使針頭堵塞，故而選用純凈水。裝置設(shè)計完成后，通過測定降雨器各個部位的雨強來檢驗雨滴發(fā)生裝置的穩(wěn)定性。

1.1.2單雨滴擊濺裝置單雨滴擊濺試驗裝置如圖1-b所示，采用直徑為10 cm、高1.5 cm的蒸發(fā)皿。

1.1.3多雨滴徑流收集裝置多雨滴徑流收集裝置如圖1-c所示，該裝置是一個長90 cm、寬15 cm、高15 cm、上開口的長方體，水槽前方進水口用量程較小的壓力表調(diào)節(jié)進水流量，以便控制水槽前段水箱中的水頭高度；緩流多孔板主要起到保持水流穩(wěn)定性的作用；在水槽前后兩處各裝1個探針, 用于測定薄層厚度；在土槽底部裝有直徑1 cm的土樣飽和供水管,以便使供試土樣能夠緩慢飽和。整個裝置采用單循環(huán)系統(tǒng), 由水泵、閥門、壓力表、進水箱、探針和飽和管組成。

1.2土樣的采集及裝填

試驗使用土壤為陜西楊凌的土婁土、延安的黃綿土和志丹縣的黑壚土，各土樣理化性狀見表1。供試土樣采集后進行風(fēng)干，去除植物枯落物、根系，過篩。單雨滴試驗裝填3 mm厚過篩(1 mm篩孔)后的土樣，向蒸發(fā)皿中緩慢加水待土樣飽和；多雨滴試驗前先在試驗槽底均勻鋪設(shè)5 cm厚的河沙(過5 mm篩孔), 再裝填厚度為5 cm、體積為4 500 cm3過篩(2.5 mm篩孔)后的土樣。裝填結(jié)束后, 用特制平尺將表面刮平,打開土樣飽和供水管，緩慢供水，使土樣飽和6 h,保證試驗過程中無下滲水流出現(xiàn)。

圖 1　雨滴發(fā)生器試驗裝置示意圖(單位：mm)

土類Soil母質(zhì)Matter有機質(zhì)/(g·kg-1)Organiccarbon全氮/(g·kg-1)Totalnitrogen全磷/(g·kg-1)Totalphosphorus顆粒組成/%Particlecomposing<2μm≥2～<10μm≥10～<50μm≥50～<100μm≥100μm黃綿土Huangmiansoil黃土Loess27.01±13.891.14±0.61.21±0.665.949.124.712.28.1土婁土Lousoil黃土Loess14.61.040.623.128.751.17.99.2黑壚土Heilusoil沙黃土SandyLoess7.5±2.50.065±0.0350.16±0.0111.150.721.310.46.5

1.3試驗方法

1.3.1單雨滴擊濺試驗單雨滴擊濺試驗時，需將圖1-a裝置放置一定高度h(h可根據(jù)實際降雨時雨滴終點速度算出)，不同的針頭號根據(jù)試驗需要安插好，左邊進水管高度根據(jù)試驗需要調(diào)節(jié)。因為是單雨滴試驗所以只插入1個針頭。

針頭雨滴直徑采用單雨滴法測定，因為一滴雨滴從高處落下極易受到風(fēng)力的影響，因此要盡量減少風(fēng)對試驗的影響，試驗時用PVC管是為了使風(fēng)力對雨滴的影響降到最低。取試驗針頭，將雨滴通過豎直的PVC管從6.8 m的高處降落，使雨滴盡量滴在蒸發(fā)皿的圓心處，記錄滴數(shù)后,測算單個雨滴質(zhì)量,根據(jù)水溫查出相應(yīng)水的密度,計算單個雨滴體積, 根據(jù)球體體積,計算單個雨滴直徑。每個針頭試驗設(shè)3個不同水頭高度，每一高度重復(fù)測定5次，結(jié)果取平均值，分別用5、6、7、9、12、16號針頭重復(fù)試驗，共測90組數(shù)據(jù)，選出雨滴直徑相差基本為0.4 mm的5#、7#、9#針頭數(shù)據(jù)，結(jié)果如表2所示。

將飽和后土樣保持一定厚度的水層，根據(jù)選定的針頭進行試驗，試驗水層厚度分別為3，5，7，10 cm。然后換針頭換土樣重復(fù)上述試驗，測定泥沙濺蝕量，每組試驗重復(fù)3次，試驗共108組。

1.3.2多雨滴試驗多雨滴人工降雨的降雨強度采用直接法測定。對每一針號在水槽底部插264個針頭,連續(xù)降雨10 min,測定降水體積,根據(jù)降雨面積和降雨時間計算降雨強度。每一針號重復(fù)測定5次，取平均值即為該針號降雨強度。試驗共15組，最終選取0.057，0.228，0.724 mm/min雨強進行試驗。

雨滴下降過程中在豎直方向受重力和空氣阻力的作用，水平方向受到氣流的作用。為此,通過雨滴降落高度的反復(fù)試驗,并根據(jù)沙玉清修正的牛頓終點速度公式[21]，當(dāng)雨滴直徑d>1.9 mm時，其終點速度為：

(1)

式中：v為終點速度，m/s；d為雨滴直徑，mm。

根據(jù)李光錄[19]降雨能計算公式，單位面積產(chǎn)生的雨滴動能為：

(2)

式中：A為降雨面積,m2；ρ為水的密度 ,1×103kg/m3；I為一場降雨的雨強,m/s；t為降雨歷時,s；θ為地表坡度, °。

人工降雨雨滴的主要參數(shù)如表2所示。

表 2　人工模擬降雨雨滴的主要參數(shù)

1.3.3薄層水流的模擬如圖1-c所示的多雨滴徑流收集裝置，在其進水口上側(cè)放置1個壓力表，在壓力表2側(cè)各設(shè)置1個調(diào)節(jié)閥，以調(diào)節(jié)進水量大小，進而調(diào)節(jié)儲水量高度以控制薄層水流厚度。在試驗前，利用水槽后部飽和管進水口緩慢進水，先浸透水槽內(nèi)沙土及試驗土樣，讓坡面土樣入滲飽和，以便試驗中薄層水流厚度均勻光滑。水槽通過下部支架調(diào)節(jié)到試驗所需要的坡度(試驗坡度分別為0°，2°，4°)。

1.4試驗步驟

試驗在有雨滴和無雨滴打擊下連續(xù)進行。首先調(diào)整水槽坡度,啟動水泵,調(diào)節(jié)閥門,待水流穩(wěn)定后, 測定無雨滴打擊時前后測針(測針精度為0.02 mm)水流深度各3次, 取平均值即可得流量和坡度水深，然后每間隔1 min取該水流下泥沙，經(jīng)沉淀、分離、烘干并稱質(zhì)量，獲得水沙量，重復(fù)5次。 每個測次均需更換表層供試土樣，試驗采用坡度分別為0°，2°和4°，雨強分別為0，0.055 7，0.228和0.724 mm/min, 水深為0，2，4，6 mm。一個土樣試驗完成后換另一種土樣重復(fù)上述試驗，3個土樣測完后再換另一號針頭重復(fù)試驗。每組重復(fù)3次，共進行144組試驗。泥沙收集按《水土保持試驗規(guī)范》[22]進行，試驗過程中均無細溝形成。

2結(jié)果與分析

雨滴擊濺坡面薄層水流,將靜止坡面土粒濺起進入到水流中, 泥沙的輸移過程便從此開始。進入水流的土樣顆粒上升到一定高度后, 便會下落到坡面, 在土粒上升和下落的交替過程中, 流動水流對土樣顆粒進行搬運。一旦濺起土粒下落到下游坡面, 再一次的啟動將決定于下一個雨滴的擊濺擾動作用。

2.1單雨滴擊濺時泥沙濺蝕量與水層厚度的關(guān)系

圖2顯示了單雨滴擊濺時，黑壚土、黃綿土、土婁土在雨滴直徑分別為2.22，2.68，3.04 mm時，隨著水層厚度的增加泥沙濺蝕量的變化。由圖2-a可知，在雨滴直徑為2.22 mm時，隨水層厚度的增加，黑壚土泥沙濺蝕量由0.018 g/(cm2·min)減小為0 g/(cm2·min)，黃綿土由0.013 g/(cm2·min)減小為 0 g/(cm2·min)，土婁土由0.006 g/(cm2·min)減小為0 g/(cm2·min)；在水層厚度為7 mm時，黑壚土、黃綿土、土婁土泥沙濺蝕量分別為0.05，0.04和0.012 g/(cm2·min)，說明雨滴直徑為2.22 mm時擾動最大水層厚度是介于7 和10 mm之間的。由圖2-b可知，在雨滴直徑為2.68 mm時，黑壚土泥沙濺蝕量由0.032 g/(cm2·min)減小為0.005 g/(cm2·min)，黃綿土由0.021 g/(cm2·min)減小為0.005 g/(cm2·min)，土婁土由0.009 g/(cm2·min)減小為0.002 g/(cm2·min)，說明擾動最大水層厚度大于10 mm。由圖2-c可知，在雨滴直徑為3.04 mm時，黑壚土泥沙濺蝕量由0.032 g/(cm2·min)減小為0.006 g/(cm2·min)，黃綿土由0.040 g/(cm2·min)減小為0.016 g/(cm2·min)，土婁土由0.009 g/(cm2·min)減小為0.002 g/(cm2·min)，說明擾動最大水層厚度大于10 mm。當(dāng)雨滴直徑由2.22 mm增加到3.04 mm時，平均擾動水層厚度由10 mm增加到14.67 mm。從趨勢上可以看出，隨著水層厚度的增加，泥沙濺蝕量呈減少趨勢。運用SPSS 21.0軟件對單雨滴泥沙濺蝕量和水層厚度進行擬合，發(fā)現(xiàn)二者呈現(xiàn)對數(shù)函數(shù)關(guān)系：

q=alnh+b。

(3)

式中：q為單雨滴泥沙濺蝕量，g/(cm2·min)；h為水層厚度，mm；a，b為經(jīng)驗系數(shù)。

圖 2　單雨滴擊濺時泥沙濺蝕量與水層厚度的關(guān)系

單雨滴擊濺時泥沙濺蝕量與水層厚度的擬合關(guān)系系數(shù)見表3。

表 3　單雨滴擊濺時泥沙濺蝕量與水層厚度的擬合系數(shù)

由表3可知，擬合相關(guān)系數(shù)均大于0.899，所以相關(guān)性很高。運用擬合公式，計算直徑2.68和 3.04 mm雨滴擊濺下，泥沙濺蝕量為0時的最大擾動水層厚度，結(jié)果如表4所示。

表 4　單雨滴擊濺時泥沙濺蝕量為0時的水層厚度

由圖2可知，直徑2.22 mm的雨滴最大擾動水層厚度為7～10 mm，當(dāng)雨滴直徑為2.68和3.04 mm時，其最大擾動水層厚度均大于10 mm。根據(jù)擬合的單雨滴泥沙濺蝕量與水層厚度函數(shù)方程關(guān)系計算結(jié)果(表4)可知，雨滴直徑為2.68 mm時，黑壚土、黃綿土、土婁土的最大擾動水層厚度分別為 12.749，13.053和14.343 mm；雨滴直徑為3.04 mm時，黑壚土、黃綿土、土婁土的最大擾動水層厚度分別為14.585，14.647和15.094 mm。

2.2多雨滴擊濺時泥沙濺蝕量與薄層水流厚度的關(guān)系

圖3顯示，坡度為0°條件下,土婁土、黃綿土、黑壚土在雨滴直徑為2.22，2.68，3.04 mm時泥沙濺蝕量隨薄層水流厚度的增加而變化。土婁土在0°坡面上,當(dāng)雨滴直徑為3.04 mm，薄層水流厚度從0增加到6 mm，降雨侵蝕擾動所產(chǎn)生的泥沙濺蝕量由0.29 g/(cm2·min)增加至21.65 g/(cm2·min)又減小至7.03 g/(cm2·min)；當(dāng)雨滴直徑為2.68 mm，薄層水流厚度從0增加到6 mm，泥沙濺蝕量由0.36 g/(cm2·min)增加至13.04 g/(cm2·min)又減小至 3.60 g/(cm2·min)；當(dāng)雨滴直徑為2.22 mm時, 薄層水流厚度從0增加到6 mm，泥沙濺蝕量由0.37 g/(cm2·min)增加至7.36 g/(cm2·min)又減小至1.98 g/(cm2·min)。在不同雨滴直徑下，黃綿土、黑壚土的濺蝕量變化趨勢與土婁土相同。說明在無薄層水流時，泥沙主要由擊濺產(chǎn)生，且產(chǎn)生泥沙量較小，隨薄層水流厚度的增大，雨滴擊濺對土樣的侵蝕作用逐漸減小，泥沙主要由薄層水流的沖刷作用產(chǎn)生。

圖 3　0°坡面上不同雨滴直徑下泥沙濺蝕量與薄層水流厚度的關(guān)系

根據(jù)上述分析，采用回歸模型對0°坡面上不同雨滴直徑下的泥沙濺蝕量與薄層水流厚度的關(guān)系進行擬合，得到多元回歸方程為：

QS=ah′2+bh′+c。

(4)

式中：QS為多雨滴泥沙濺蝕量，g/(cm2·min)；h′為薄層水流厚度，mm；a、b為經(jīng)驗系數(shù)；c為經(jīng)驗常數(shù)。

0°坡面上泥沙濺蝕量與薄層水流厚度的擬合關(guān)系系數(shù)見表5。

表 5　0°坡面上泥沙濺蝕量與薄層水流厚度的擬合系數(shù)

圖4顯示，2°坡面上，土婁土、黃綿土、黑壚土在雨滴直徑分別為2.22，2.68，3.04 mm時泥沙濺蝕量隨薄層水流厚度的增加其變化趨勢與坡度為0°時相同，說明坡度因子對泥沙的產(chǎn)生有重要影響，其泥沙濺蝕量也可以用方程(2)表示。鑒于4°坡面上，土婁土、黃綿土、黑壚土在雨滴直徑分別為2.22，2.68，3.04 mm時泥沙濺蝕量隨薄層水流厚度的增加的變化趨勢與坡度為0°和2°時相同，故在此未列出試驗結(jié)果。

2.3多雨滴擊濺時泥沙濺蝕量與雨滴動能的關(guān)系

鑒于多雨滴擊濺時黃綿土、土婁土、黑壚土泥沙濺蝕量與雨滴動能關(guān)系的變化一致，故本研究僅對黃綿土在不同薄層水流厚度下泥沙濺蝕量與雨滴動能關(guān)系進行分析，結(jié)果見圖5。

圖 4　2°坡面上不同雨滴直徑下泥沙濺蝕量與薄層水流厚度的關(guān)系

圖 5不同薄層水流厚度下黃綿土

泥沙濺蝕量與雨滴動能的關(guān)系

◆.無薄層水流；■.2 mm薄層水流；▲.4 mm薄層水流

Fig.5Relationship between sediment splash mass and

raindrop kinetic energy with different shallow flow

depths of Huangmian soil

◆.Shallow flow depth of 0 mm；■.Shallow flow

depth of 2 mm；▲.Shallow flow depth of 4 mm

圖5顯示，當(dāng)雨滴動能由68.18×10-3J增加到94.97×10-3J時，0 mm薄層水流下泥沙濺蝕量由8.23 g/(cm2·min)增加到14.73 g/(cm2·min)，2 mm薄層水流下由27.72 g/(cm2·min)增加到42.35 g/(cm2·min),4 mm薄層水流下由84.07 g/(cm2·min)增加到104.68 g/(cm2·min)。運用SPSS 21.0軟件對不同薄層水流條件下的泥沙濺蝕量與雨滴動能的關(guān)系進行擬合，發(fā)現(xiàn)二者呈冪函數(shù)關(guān)系，可用以下方程表示：

Qs=aeEs×b。

(5)

式中：Qs為多雨滴泥沙濺蝕量，g/(cm2·min)；Es為雨滴動能，×10-3J；a、b為經(jīng)驗系數(shù)。

黃綿土在不同薄層水流厚度下泥沙濺蝕量與雨滴動能的擬合關(guān)系系數(shù)見表6。由圖5及表6可知，對黃綿土而言，當(dāng)薄層水流厚度一定時，泥沙濺蝕量隨雨滴動能的增加而增大。雨滴動能由68.16×10-3J增加到94.97×10-3J，在無薄層水流時泥沙濺蝕量由8.23 g/(cm2·min)增加到 17.42 g/(cm2·min)，且薄層水流厚度越高泥沙濺蝕量增加越多。這是由于薄層水流的沖刷是泥沙產(chǎn)生的主要動力，雨滴擊濺擾動薄層水流進而增加泥沙沖刷能力，當(dāng)雨滴動能增加時其擾動能力越強，因此泥沙濺蝕量增加越多。

表 6　黃綿土在不同薄層水流厚度下泥沙

3結(jié)論

1)單雨滴擊濺時，泥沙濺蝕量(q)隨水層厚度(h)的增加而增大，雨滴直徑越大則擾動水層厚度越大，其關(guān)系可用對數(shù)函數(shù)關(guān)系式q=alnh+b表示，當(dāng)雨滴直徑為2.22 mm時擾動最大水層厚度為7～10 mm；當(dāng)雨滴直徑為2.68和3.04 mm時，其擾動最大水層厚度均大于10 mm。根據(jù)擬合函數(shù)計算可得，直徑2.68 mm雨滴擾動的黑壚土、黃綿土、土婁土的最大水層厚度分別為12.749，13.053，14.343 mm；3.04 mm雨滴擾動的黑壚土、黃綿土、土婁土的最大水層厚度分別為14.585，14.647，15.094 mm。

2)坡度相同時，不同土壤(黃綿土、土婁土、黑壚土)及不同雨滴直徑(2.22，2.68，3.04 mm)條件下，當(dāng)薄層水流厚度(h′)從0增加到6 mm時，其與降雨擊濺混合侵蝕下泥沙濺蝕量(Qs)曲線總體呈先上升后下降的趨勢，其關(guān)系可用二次拋物線方程Qs=ah′2+bh′+c表示。

3)薄層水流厚度一定時，泥沙濺蝕量(Qs)隨雨滴動能(Es)的增加而增大，Es由68.16×10-3J增加到94.97×10-3J時，無薄層水流的泥沙濺蝕量由8.23 g/(cm2·min)增加到17.42 g/(cm2·min)，且薄層水流厚度越高泥沙濺蝕量增加越多，其關(guān)系可用冪函數(shù)Qs=aeEs×b表示。這是由于薄層水流的沖刷是泥沙產(chǎn)生的主要動力，雨滴擊濺擾動薄層水流進而增加泥沙沖刷能力，當(dāng)雨滴動能增加時其擾動能力越強，因此泥沙濺蝕量增加越多。

[參考文獻]

[1]吳普特,周佩華.地表坡度對雨滴濺蝕的影響 [J].水土保持通報,1991,11(3):8-13.

Wu P T,Zhou P H.The effect of land slope upon raindrop splash erosion [J].Bulletin of Soil and Water Conservation,1991,11(3):8-13.(in Chinese)

[2]張洪江.土壤侵蝕原理 [M].北京:中國林業(yè)出版社,2000.

Zhang H J.Principle of soil erosion [M].Beijing:China Forestry Publishing House,2000.(in Chinese)

[3]吳普特,周佩華.雨滴擊濺在薄層水流侵蝕中的作用 [J].水土保持通報,1992,12(4):19-26,47.

Wu P T,Zhou P H.The effect of rain drop splash on the sheet flaw hydraulic friction factor [J].Journal of Soil and Water Conservation,1992,12(4):19-26,47.(in Chinese)

[4]趙春紅,高建恩,王宏杰,等.降雨對坡面薄層水流泥沙起動的影響 [J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2013(6):1057-1069.

Zhao C H,Gao J E,Wang H J,et al.Effect of rainfall on sediment incipient motion of overland flow [J].Journal of Basic Science and Engineering,2013(6):1057-1069.(in Chinese)

[5]鄭粉莉.黃土區(qū)坡耕地細溝間侵蝕和細溝侵蝕的研究 [J].土壤學(xué)報,1998,35(1)：95-103.

Zheng F L.Study on interrill erosion and rill erosion on slope farmland of Loess Area [J].Acta Pedologica Sinica，1998，35(1)：95-103.(in Chinese)

[6]吳普特,周佩華.雨滴擊濺對坡面薄層水流阻力系數(shù)的影響 [J].水土保持學(xué)報,1994,8(2):40-42.

Wu P T,Zhou P H.The effect of rain drop splash on the sheet flaw hydraulic friction factor [J].Journal of Soil and Water Conservation,1994,8(2):40-42.(in Chinese)

[7]吳普特.動力水蝕實驗研究 [M].西安:科學(xué)技術(shù)出版社,1997:9.

Wu P T.Experimental study on soil erosion dynamic mechanism [M].Xi’an:Science Press,1997:9.(in Chinese)

[8]張光輝.坡面薄層流水動力學(xué)特征的實驗研究 [J].水科學(xué)進展,2002,13(2):159-165.

Zhang G H.Study on hydraulic properties of shallow flow [J].Advance in Water Science,2002,13(2):159-165.(in Chinese)

[9]Jayawardena A W,Rezaur R B.Drop size distributions and kinetic energy load of rain storms in Hong Kong [J].Hydrol Processes,2004(14):1069-1082.

[10]Par sons A J,Gadian A M.Uncertainty in modeling the detachment of soil by rainfall [J].Earth Surf Processes Landforms,2000(25):723-728.

[11]Zhang X C,Nearing M A,Miller W P,et al.Modeling interrill sediment delivery [J].Soil Science Society of America Proceedings,1998,62(2):438-444.

[12]閆美華.雨滴擊濺及薄層水徑流動力對坡面侵蝕影響的研究 [J].中國農(nóng)村水利水電,2011(3):66-69

Yan M H. Study on the effect of raindrops splash and runoff power of shallow flow on slope erosion [J].China Rural Water and Hydropower,2011(3):66-69.(in Chinese)

[13]Kinnel P I A.Sediment consent ration resulting from flow depth/drop size interact ions in shallow overland flow [J].Trans ASAE,1993(36):1099-1103.

[14]Kinnel P I A.Sediment transport by shallow flows impacted by pulsed artificial rainy all [J].Aust J Soil Sci,1993(31):199-207.

[15]Kinnel P I A.Sediment transport by medium to large tropism packing flows at sub terminal velocity [J].Soil Sci Soc Am J,2005(69):902-905.

[16]Kinnel P I A.The influence of flow discharge on sediment concentration s in raindrop induced flow transport [J].Aust J Soil Sci,1988(26):575-582.

[17]Kinnel P I A.T he mechanics of raindrop-reduced flow transport [J].Aust J Soil Sci,1990(28):497-516.

[18]Gabet E J,Dunne T.Sediment detachment by rain power [J].Water Resources Research,2003,39(1):1002-1014.

[19]李光錄,吳發(fā)啟,龐小明，等.泥沙輸移與坡面降雨和徑流能量的關(guān)系 [J].水研究進展,2008,11(6):868-874.

Li G L,Wu F Q,Pang X M,et al.Relationship between sediment transport with surface rainfall and runoff energies on sloping [J].Advances in Water Science,2008,11(6):868-874.(in Chinese)

[20]李光錄,吳發(fā)啟,趙小風(fēng)，等.雨滴擊濺下薄層水流的輸沙機理研究 [J].西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版,2009,37(9):149-154.

Li G L,Wu F Q,Zhao X F,et al.The mechanics of sediment transport under raindrops impacting shallow overland flow [J].Journal of Northwest A&F University：Nat Sci Ed,2009,37(9):149-154.(in Chinese)

[21]劉秉正,吳發(fā)啟.土壤侵蝕 [M].西安:陜西人民出版社,1997:34-35.

Liu B Z,Wu F Q.Soil eros ion [M].Xi’an:Shaanxi People Press,1997:34-35.(in Chinese)

[22]中華人民共和國水利電力部.水土保持試驗規(guī)范 [M].北京：中國水電出版社,1991.

Ministry of Water Resources and Electric Power in P R China.Experimental standard of soil and water conservation [M].Beijing:Chinese Water and Electricity Press,1991.(in Chinese)

Rain splash and sheet flow mixed sediment transport mechanism

ZHENG Teng-hui1a,XING Yuan-yuan1b,2,HE Kai-xuan1b,TENG Fei1b,LI Guang-lu1a,1b

(1 aInstituteofSoilandWaterConservation,bCollegeofNaturalResourcesandEnvironment,NorthwestA&FUniversity,Yangling,Shaanxi712100,China；2XinjiangInstituteofEcologyandGeography，ChineseAcademyofSciences，Urumqi，Xinjiang830011,China)

Abstract:【Objective】 Effect of raindrop splash on slope runoff and sediment transport was studied in this paper.【Method】 In the indoor artificial rainfall simulation, different slopes (0°,2° and 4°),different sheep flow thicknesses (0,2,4 and 6 mm) and different raindrop diameters (2.22,2.68 and 3.04 mm) were tested.The relationships between sediment splash erosion amount and the deep layer of water for single raindrop splashing and between sediment splash erosion amount and sheet flow thickness for multiple raindrop splashing were analyzed and fitted.The relationship of sediment splash erosion amount of yellow soil and raindrop kinetic energy under different sheet flow thicknesses was also analyzed.【Result】 For single raindrop splashing,splash amount decreased with the increase of sheet flow thickness, raindrop diameter increased with the increased of disturbance layer thickness,and the average thickness of water resisting disturbance increased from 10 to 14.67 mm when raindrop diameter increased from 2.22 to 3.04 mm.For multiple raindrop splashing,sediment splash erosion amount augmented firstly then decreased as the increased of sheet flow thickness.With same thickness,sediment splash erosion amount increased with the increase of raindrop diameter.When raindrop kinetic energy remained same,sediment splash erosion amount reduced with the increase of shallow flow thickness.When thickness did not change,sediment splash erosion amount increased with the increase of raindrop kinetic energy.【Conclusion】 This study established the relationships between single raindrop splash amount and sediment layer thickness as well as between sediment splash erosion amount and multiple raindrops splash amount.

Key words:artificial precipitation;shallow flow;mixed erosion;rain kinetic energy;raindrop splash;sediment transport mechanism

DOI：網(wǎng)絡(luò)出版時間:2016-02-0209:3710.13207/j.cnki.jnwafu.2016.03.029

[收稿日期]2014-07-08

[基金項目]國家級大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃項目(1210712061)

[作者簡介]鄭騰輝(1990-)，男，河南商丘人，在讀碩士，主要從事水土保持與荒漠化防治研究。E-mail:1364162509@qq.com[通信作者]李光錄(1964-)，男，甘肅永靖人，副教授，博士，主要從事水土保持與土地利用研究。E-mail: guangluli@nwsuaf.edu.cn

[中圖分類號]S157.1

[文獻標志碼]A

[文章編號]1671-9387(2016)03-0211-08