基于人工模擬降雨的黃潮土坡面侵蝕規(guī)律研究

許 滸，王 健，索梅芹，王振龍，胡永勝，周 超

(1.安徽省(水利部淮委)水利科學研究院水利水資源安徽省重點實驗室，安徽 蚌埠 233000；2.河北工程大學，河北 邯鄲 056021)

0 引 言

淮北平原多坡地，平原坡地占比為68%，地勢西北高而東南低。其中，黃潮土占淮北平原總面積的33%，主要分布在蕭縣、碭山、亳州等平原區(qū)和低山殘丘坡地區(qū)，最大坡度在10°左右。該地區(qū)自然降水資源的70%以上集中在夏季，多以暴雨形式出現(xiàn)，地表徑流量大，易造成土壤侵蝕。杜明成[1]等通過人工模擬降雨實驗分析了黃潮土在多雨強變坡度產(chǎn)流產(chǎn)沙規(guī)律，認為雨強對產(chǎn)流產(chǎn)沙的影響略大于坡度。雨強作為與土壤侵蝕相關(guān)的重要的條件，已有許多研究成果，An[2]、Wu[3]等通過對不同雨強和不同性質(zhì)土壤的降雨侵蝕實驗認為雨強是與土壤侵蝕關(guān)系最密切的降雨因子。張夢[4]等通過人工降雨試驗對黃土坡面進行侵蝕研究表明，在相同坡度下，坡面產(chǎn)流量和產(chǎn)沙量與降雨強度呈正相關(guān)，降雨過程中侵蝕強度上下波動，降雨強度越大，波動值越大。姜芃等利用稀土示蹤法在人工降雨條件下在南方坡耕地水土流失研究中得出坡面侵蝕量隨著降雨強度的增大而增加；不同坡段的侵蝕量有所不同，上坡段占總侵蝕量的10%左右，而中、下坡段占總侵蝕量的90%左右[5]。坡面土壤侵蝕地形條件包括坡度、坡長等，地形條件與侵蝕量關(guān)聯(lián)性很大[6-10]。在不同坡度下紫色土土壤侵蝕的研究顯示，產(chǎn)沙率隨坡度的增大速率先上升而后趨于穩(wěn)定，坡度對產(chǎn)沙量的影響最大[11]。張少博等基于前人對臨界坡度的研究，分析了坡度對伊犁河流域土壤侵蝕的影響，發(fā)現(xiàn)該流域土壤侵蝕區(qū)域主要集中在10°～15°的雪線以下至洪積平原由陡變緩的區(qū)域，此區(qū)域土壤團粒結(jié)構(gòu)易被破壞，進而形成土壤侵蝕[12]。陳思旭[13]等通過對南方丘陵區(qū)土壤進行降雨侵蝕研究得出海拔200～500 m坡度在15°～25°范圍的地區(qū)的侵蝕最為嚴重。本文采用人工模擬降雨方法，通過對不同雨強和坡度下黃潮土產(chǎn)流產(chǎn)沙及入滲分析，得出降雨侵蝕規(guī)律，對淮北平原水土保持和環(huán)境保護具有參考意義。

1 材料與方法

1.1 實驗設計

實驗時間為2020年6月～7月，在五道溝水文實驗站人工降雨實驗場進行。人工模擬降雨裝置為西安清遠QYJY—503型移動式人工模擬降雨系統(tǒng)，徑流小區(qū)為固定坡長可調(diào)坡度的變坡鋼槽，尺寸為8 m×4 m×2 m，回填黃潮土取自蕭縣楊樓鎮(zhèn)農(nóng)田，取土后分層回填并自然密實，底部設置30 cm過濾層。實驗區(qū)近30年實測最大雨強為78 mm/h，本次實驗選定40、60、80 mm/h 3種雨強和3°、6°、9°、12°等4種坡度。每組實驗重復2次，共進行了36組人工模擬降雨實驗，每次降雨歷時60 min。

1.2 實驗過程

實驗過程中產(chǎn)流量和泥沙濃度數(shù)據(jù)由位于鋼槽底部傳感器間隔1 min自動采集，侵蝕量由測量的泥沙濃度與產(chǎn)流量的乘積得出。不同埋深的土壤含水率由預先埋設的TDR315H土壤水分傳感器測出。鋼槽內(nèi)沿坡面方向上、中、下3個位置，分別在垂向距地表0.1、0.2、0.3、0.5、0.7、0.9、1.5 m深處埋設了TDR315H土壤水分傳感器共21個，采集頻率為1次/10 min。實驗在夜晚或凌晨無風的環(huán)境下進行，人工降雨裝置實驗前進行雨強率定和均勻度檢驗，實驗前24 h預先采用25 mm/h雨強進行預降雨，坡面即將產(chǎn)流時停止降雨，使每次實驗前表層土壤含水率相差不超過2%。

2 結(jié)果與分析

2.1 不同坡度、雨強對坡面產(chǎn)流過程的影響

不同雨強、坡度實測產(chǎn)流過程如圖1所示。由圖1可知，各雨強產(chǎn)流量隨坡度的變化均有較好的規(guī)律性：在降雨初期，經(jīng)過短暫的初損時間后坡面流會迅速增大直至穩(wěn)定產(chǎn)流，這與杜明成[1]等實驗結(jié)果相似，產(chǎn)生這一現(xiàn)象是因為隨著土壤入滲使含水率增大，土壤入滲率迅速減小達到穩(wěn)滲，而雨強保持不變，因此會形成穩(wěn)定產(chǎn)流。3種雨強穩(wěn)定產(chǎn)流時間會隨著坡度的增大而減少，其中坡度為3°的穩(wěn)定產(chǎn)流時間近20 min，而6°、9°和12°條件下的穩(wěn)定產(chǎn)流時間則大幅縮短，尤其在大雨強下迅速達到穩(wěn)定產(chǎn)流。3種雨強3°和6°的穩(wěn)定徑流量都十分接近，在實驗過程中可以發(fā)現(xiàn)實驗進行20 min以后3°坡面底部1 m長的坡段內(nèi)會產(chǎn)生積水。這是由于在雨強相同情況下3°的坡度較小導致坡面流流速低，降雨產(chǎn)生的徑流不能迅速從集水口排除。雖然6°時鋼槽底部未產(chǎn)生積水，但流速和3°時相比幾乎沒有變化，導致3°和6°情況下坡面徑流量出現(xiàn)接近的情況。而隨著坡度的增加達到9°時坡面流流速迅速增加，徑流流經(jīng)坡面的時間縮短，也減少了入滲水量，這兩方面原因使徑流量在9°時出現(xiàn)了較大幅度的增長，3種雨強穩(wěn)定坡面徑流量分別達到了30.02、32.22 L/min和36.51 L/min，較9°時徑流量增加了41.6%、41.3%和39.7%。而隨著坡度繼續(xù)增大至12°，3種雨強坡面徑流量都出現(xiàn)了較9°時減小的情況，原因在于變坡鋼槽在坡度增加的過程中有效承雨面積會相應減少。根據(jù)三角函數(shù)余弦公式可知，在角度較小時坡度的增大坡面的承雨面積變化量輕微，隨著角度的逐漸增大，坡面承雨面積減小的幅度會迅速增加，坡度增大導致的承接雨量減少不能忽略。

圖1 相同雨強不同坡度坡面徑流過程

2.2 不同坡度、雨強對坡面侵蝕過程的影響

相同雨強不同坡度的坡面侵蝕過程規(guī)律較產(chǎn)流過程復雜，侵蝕量和徑流量、泥沙濃度大小關(guān)系密切，坡面侵蝕過程如圖2所示。在40 mm/h雨強下，坡面侵蝕率與坡度呈正相關(guān)，在3°坡和6°坡侵蝕產(chǎn)沙量由0增至穩(wěn)定值需要20 min，上升過程較為平緩，而隨著坡度繼續(xù)增大至9°和12°，坡面出現(xiàn)徑流并迅速加大，產(chǎn)沙量也隨著迅速增大至穩(wěn)定。這一過程的泥沙濃度見表1。在6°時由13.11 g/L增至16.37 g/L，增加了24.86%，在9°和12°時泥沙濃度有所下降。12°時前期產(chǎn)沙量有較大幅度的波動，這是因為隨著坡度的增大，土壤顆粒受重力方向的分離也隨之增大，降雨初期土坡面表層的土壤結(jié)構(gòu)比較松散，抗沖能力差，更容易被運輸至坡底，坡度的增大也會形成更快的流速，水流對泥沙的運輸能力增強，從而導致了前期坡面侵蝕產(chǎn)沙量數(shù)值變化較大，經(jīng)過20 min表面松散土壤顆粒被降雨徑流沖刷殆盡后侵蝕量波動減小，逐漸趨于穩(wěn)定。

圖2 相同雨強不同坡度下坡面侵蝕過程

表1 不同坡度雨強下穩(wěn)定產(chǎn)流泥沙濃度

在60 、80 mm/h兩種雨強下坡度對坡面侵蝕的影響具有極高的相似性，坡度從3°增至6°后產(chǎn)沙量出現(xiàn)了大幅度增加，從表中可以看出6°時泥沙濃度達到最大值，分別較3°時增加了122.68%和205.76%，泥沙濃度的陡增使得侵蝕量也大幅增加，在3°～12°坡度間達到了最大值。隨著坡度的繼續(xù)增大到9°，雖然坡面徑流量較6°時有所增加，但此時泥沙濃度與3°時相比較為接近；因此，坡面產(chǎn)沙量出現(xiàn)了大幅回落。當坡度達到12°時，侵蝕率增幅不明顯，坡面產(chǎn)流使初期土壤表面的不穩(wěn)定的顆粒被迅速沖刷至坡底，坡面侵蝕產(chǎn)沙量出現(xiàn)短時間大幅度變化。隨著坡面徑流量的減小，產(chǎn)流穩(wěn)定后的侵蝕量進一步減小。60、80 mm/h兩種雨強下6°坡度時降雨侵蝕量達到峰值，而40 mm/h雨強下雖然侵蝕量未達到最大值，但此時的泥沙濃度最大，說明6°為黃潮土的臨界坡度。

2.3 不同坡度累計產(chǎn)流量與累計產(chǎn)沙量的關(guān)系

由降雨引起的坡面土壤侵蝕量與坡面徑流量有著十分密切的關(guān)系，不同雨強和坡度實測1 h累計產(chǎn)流量和侵蝕量數(shù)據(jù)見表2。

表2 不同雨強和坡度累計產(chǎn)流量和產(chǎn)沙量

累計產(chǎn)流量與累計侵蝕量擬合曲線如圖3所示，點子的分布可以按坡度分為兩部分分別擬合線性曲線。其中，坡度6°的產(chǎn)流量與侵蝕量單獨擬合，其余3種坡度的產(chǎn)流量與產(chǎn)沙量一起擬合效果較好，擬合精度分別達到了0.999 7和0.894 6，擬合公式分別為y=0.113x-111.9和y=0.019x-6.216。式中，x為總產(chǎn)流量，L；y為總侵蝕量，kg。

圖3 累計產(chǎn)流量與累計侵蝕量擬合曲線

2.4 不同坡度、雨強坡面入滲特征分析

降雨經(jīng)過植物截留、雨間蒸發(fā)、填洼及入滲后，剩余水量形成地表徑流，在裸土的降雨徑流過程中，入滲量為主要的徑流損失。產(chǎn)流形式分為超滲產(chǎn)流和蓄滿產(chǎn)流，不論哪種形式的產(chǎn)流，都要滿足降雨量超過入滲量。在降雨初期，土壤較為干燥，降雨受分子力、毛管力和重力共同作迅速入滲，隨著降雨的持續(xù)，毛管孔隙逐漸充滿，土壤含水率增大，入滲速率逐漸降低，當土壤孔隙充水達到飽和時，入滲速率主要受重力影響。降雨強度對入滲的影響目前仍然沒有統(tǒng)一的結(jié)論，一部分學者通過實驗得出土壤入滲率和穩(wěn)定入滲率會隨著降雨強度的增大而增大，而另一部分學者認為降雨強度與降雨動能密切相關(guān)，而較大的降雨動能會加速土壤結(jié)皮，從而影響降雨入滲。實驗選取60 mm/h雨強時3°、6°、9°、12° 4個坡度和40、60 mm/h和80 mm/h 3個雨強9°坡度時的降雨入滲過程進行分析(見圖4)。

圖4 不同坡度和雨強下降雨入滲過程

由圖4可知，不同埋深的土層土壤含水率變化差異較大，埋深0.5m以內(nèi)的土壤含水率對降雨的響應較為敏感，土壤含水率隨著降雨過程的持續(xù)而逐漸增大，且表現(xiàn)為土壤埋深越大，變化程度越小，埋深大于0.5m的土壤含水率幾乎不受降雨影響。土壤含水率按照埋深自上而下逐漸增大，降雨后短時間內(nèi)便形成地表徑流，此時上層土壤含水率雖有一定程度增加，仍小于土壤飽和含水率，表現(xiàn)出超滲產(chǎn)流特征。本文選取埋深0.1m處的土壤含水率數(shù)據(jù)，通過計算土壤含水率變化速率分析坡度和雨強對降雨入滲的影響。土壤含水率變化速率

Δωi=(ωi+1-ωi)/t

(1)

式中，Δωi為土壤含水率變化速率，%/min；ωi為測量時段開始時刻的土壤含水率，%；ωi+1為測量時段結(jié)束時刻的土壤含水率，%；t為測量時段的時間，min。

60 mm/h雨強不同坡度下淺層土壤含水率變化速率見表3。由表3可以看出，60 mm/h雨強不同坡度下淺層土壤含水率變化速率總體上隨著降雨的持續(xù)而逐漸減小，符合入滲規(guī)律。觀察降雨期間不同坡度的淺層土壤含水率平均變化速率可以發(fā)現(xiàn)在坡度為6°時入滲速率最大，當坡度繼續(xù)加大，土壤含水率平均變化速率有逐漸減小的趨勢，6°是降雨入滲的臨界坡度，分析40 mm/h和80 mm/h雨強降雨入滲數(shù)據(jù)與60 mm/h雨強表現(xiàn)一致。

表3 雨強60 mm/h時不同坡度淺層土壤含水率變化速率 %

坡度9°時不同雨強下淺層土壤含水率變化速率見表4。由表4可以看出，坡度9°時淺層土壤含水率變化速率與雨強大小的關(guān)系不明顯，坡度為3°和6°時也未發(fā)現(xiàn)顯著的關(guān)聯(lián)性，這說明黃潮土雨強對坡面入滲的影響較為復雜，沒有明顯的規(guī)律性。結(jié)合之前坡度對坡面產(chǎn)沙過程的影響分析中侵蝕量的臨界坡度也為6°，臨界坡度的出現(xiàn)可以認為坡度影響降雨入滲，當達到6°左右時入滲速率最大，而黃潮土黏粒含量偏低，降雨入滲后孔隙水壓力增大，使其結(jié)構(gòu)穩(wěn)定性變差，土壤顆粒易分散，顆粒間有效應力降低，導致土體抗剪強度降低，因此入滲是影響黃潮土降雨侵蝕的重要因素[14]。

表4 坡度9°時不同雨強淺層土壤含水率變化速率 %

3 結(jié) 論

(1)相同雨強坡面產(chǎn)流量隨著坡度的增加先增大后減小，坡度9°時產(chǎn)流量最大。40 mm/h雨強坡面侵蝕量隨坡度的增加而逐漸增大，60、80 mm/h時侵蝕量隨坡度增加表現(xiàn)為先增大而后逐漸減小，坡度6°時侵蝕量最大。泥沙濃度隨著坡度的增加先增大而后降低，在6°時達到峰值，6°為黃潮土降雨侵蝕的臨界坡度。累計產(chǎn)流量與累計侵蝕量在3°和其他坡度時均呈線性函數(shù)關(guān)系(R2分別為0.999 7和0.894 6)。

(2)黃潮土埋深0.5 m以內(nèi)的土壤含水率受降雨響應較為敏感，坡度影響降雨入滲速率，淺層土壤含水率變化速率隨著坡度的增加先增大后逐漸減小，在坡度6°時達到最大值，黃潮土黏粒含量偏低，降雨入滲后結(jié)構(gòu)穩(wěn)定性變差，入滲是影響黃潮土降雨侵蝕的重要因素。

(3)本次實驗模擬降雨過程中保持雨強恒定，這與天然降雨存在差異，降雨過程中雨強的變化對坡面徑流和土壤侵蝕具有一定影響，在以后的土壤侵蝕規(guī)律研究中還應考慮天然降雨條件的影響。

致謝：本文得到“河北省水資源高效利用工程技術(shù)研究中心”的大力幫助，在此深表感謝!