不同地面處理下的降雨入滲規(guī)律研究

祝陳夢媛，付玉娟，姜國輝

(沈陽農(nóng)業(yè)大學(xué)水利學(xué)院，沈陽 110866)

土壤水是指從地面向下至地下水含水層以上土壤層中的水分，廣泛分布于陸地表層，是地表水與地下水相互轉(zhuǎn)化的紐帶，又是植物生長的必要水源[1]。土壤水循環(huán)的過程是從入滲開始的，國內(nèi)土壤水入滲的研究始于土壤入滲與其外界因素的影響研究，其后應(yīng)用了大量的入滲方程以及模擬方法，力求將土壤水入滲研究從定性描述轉(zhuǎn)化為定量描述。例如，徐克輝和盛平(1992年)[2]等人應(yīng)用Mein-Larson入滲模型來模擬上海郊區(qū)低洼旱區(qū)的入滲規(guī)律，為排水模型提供了入滲參數(shù)；陳麗華(1995年)[3]使用philip半解迭代式遞推的方法求解出了黃土林地區(qū)的土壤入滲率理論值，并通過理論值與實測值比較，發(fā)現(xiàn)二者吻合較好；趙西寧(2002年)[4]以黃土溝壑地區(qū)為實驗區(qū)，結(jié)合達西公式和kostiakov方程，計算了不同耕作類型的入滲率，表明等高耕作能大幅提高入滲率；吳欽孝和韓冰(2004年)[5]等人采用“環(huán)刀法”對森林小流域和草灌小流域土壤水分入滲特征進行了研究，結(jié)果表明森林土壤的入滲速率明顯高于草地；孫元元(2006年)[6]為了研究不同土壤在不同恒定降雨入滲條件下的土壤非飽和特性，對恒定降雨有限含水層深度的土壤在垂直和水平入滲條件下進行分析，研究發(fā)現(xiàn)：其非飽和特性值隨含水層的透水性和雨強而變化，透水性越小，雨強越大，則非飽和特征曲線凹形越甚，反之亦然；宇宙，王勇(2015年)[7]為了探索赤峰地區(qū)玉米的最佳滴灌灌溉定額和節(jié)水效果，進行的田間試驗研究了覆膜滴灌下土壤含水率變化情況，研究表明:膜下滴灌灌溉制度定在灌水5次為最佳，土壤含水率在玉米拔節(jié)期極速升高，在灌漿階段和成熟期土壤含水率維持在較低水平；趙靖丹，李瑞平等(2016年)[8]通過大田試驗研究了有膜滴灌與無膜滴灌對玉米土壤水分和溫度的變化，結(jié)果表明:地膜覆蓋條件下的玉米生育前期應(yīng)增加土壤淺層溫度，提高土壤表層含水率，為玉米的生長發(fā)育提供良好條件。除了建立入滲模型，利用入滲方程研究外，也有不少學(xué)者分析了不同土壤質(zhì)地、耕作植物、耕作方式對土壤水入滲規(guī)律的影響。隨著膜下滴灌在半干旱區(qū)的大面積推廣，覆膜起壟改變了下墊面狀態(tài)及其田間微地形，勢必會影響到土壤水分入滲過程，為了分析膜下滴灌的覆膜起壟對降雨入滲過程的改變以及對降雨利用程度的影響，試驗利用HOBO采集器及土壤含水量探針動態(tài)監(jiān)測不同下墊面、雨強及初始土壤含水量下不同剖面土壤含水量的變化情況，進而對其入滲規(guī)律進行分析。

1 材料與方法

1.1 試驗區(qū)概況

試驗在遼寧省朝陽市建平縣灌溉試驗站完成。建平屬于半干旱區(qū)，雨熱同季，全年平均氣溫7.6 ℃，最高氣溫37 ℃，最低氣溫-36.9 ℃，多年平均降水量614.7 mm，多集中在6-8月，多年平均蒸發(fā)量為1 850～2 200 mm，多年平均徑流深50～81 mm。

試驗在該試驗站測坑中進行，每個測坑2 m×2.5 m，深2 m，采用C25鋼筋混凝土澆筑，坑之間用8 mm厚鋼板隔開，鋼板高出地面0.15 m，測坑中壟溝溝尾設(shè)置排水孔，一旦有徑流產(chǎn)生，立即排出。測坑上方有遮雨棚可以避免天然降雨對試驗的影響，測坑內(nèi)的土壤物理特性見表1。

1.2 試驗設(shè)計

試驗設(shè)置3個影響因素，分別為下墊面處理、降雨強度及初始土壤含水量，每個因素設(shè)置3個水平，其中下墊面3個處理分別為覆膜起壟、起壟、平整裸地，降雨強度3個水平分別為10、20、30 mm/h，初始土壤含水量3個水平分別為田間持水量的50%、60%及70%。試驗為3因素3水平，采用正交方案進行實施，每組試驗的影響因素及水平見表2。試驗從2015年6月開始，同年9月中旬結(jié)束。

表1 試驗地土壤物理性質(zhì)

表2 試驗組數(shù)的具體參數(shù)

注：土壤初始含水率控制為田間持水量的百分比。

降雨通過人工降雨模擬器模擬不同強度的降雨過程，利用HOBO U30 NRC數(shù)據(jù)采集器動態(tài)監(jiān)測不同深度土壤含水率的變化過程。在下墊面為平整裸地時布設(shè)1個監(jiān)測剖面；在覆膜起壟及起壟兩種下墊面情況下，因為入滲過程是二維入滲，所以布設(shè)了3個監(jiān)測剖面，用來監(jiān)測不同剖面位置的土壤含水量的變化情況。測坑內(nèi)不同下墊面處理下的土壤含水率探針具體布設(shè)位置見圖1，平整裸地的傳感器布置與覆膜起壟相同。各剖面土壤水分傳感器埋設(shè)深度分別為10、30、50、70、90以及120 cm，監(jiān)測時間間隔為10 min，從降雨開始到降雨結(jié)束后48 h連續(xù)監(jiān)測。

圖1 剖面布置圖(單位：cm)

2 各因素對降雨入滲過程影響程度分析

為了分析降雨強度、土壤初始含水率、農(nóng)田下墊面處理對降雨入滲量的影響程度，根據(jù)監(jiān)測到的各剖面的土壤含水量的數(shù)據(jù)，分別選取降雨開始后3、6及12 h的土壤含水量與降雨前的土壤含水量的差值得到各剖面的降雨入滲量。

根據(jù)HOBO數(shù)據(jù)采集器中的含水量數(shù)據(jù)，利用正交試驗數(shù)據(jù)分析方法，得出降雨強度、農(nóng)田下墊面處理以及土壤初始含水率對降雨入滲量的影響程度。正交試驗數(shù)據(jù)分析結(jié)果見表3。

表3 各因素對降雨入滲量的正交分析結(jié)果(降雨3 h后)

圖2 各因素對入滲量的影響程度

通過圖3(a)可以看出，在降雨3 h時，3個因素對降雨入滲量的影響程度分別是：降雨強度>農(nóng)田下墊面處理>田間持水率，降雨強度對土壤水入滲過程有顯著影響。在農(nóng)田下墊面處理方式中，平整裸地條件下降雨入滲量最多，起壟條件次之，覆膜起壟條件下最少。覆膜起壟條件下入滲量最小的原因有兩個方面，首先由于覆膜的不透水作用，透水面大幅減少，導(dǎo)致降雨的入滲通道減少；另外覆膜區(qū)域的降雨會匯集到不覆膜區(qū)域，間接增加了降雨強度，導(dǎo)致地面產(chǎn)流提前發(fā)生，兩方面原因疊加導(dǎo)致覆膜起壟區(qū)域降雨入滲量最小。起壟條件下，由于壟溝微地形改變的原因，降雨開始后，溝內(nèi)雨水迅速匯集，也會提前產(chǎn)生徑流，導(dǎo)致土壤入滲量減少。土壤初始含水率對降雨入滲量的影響相對最小，土壤初始含水率越低，土壤的儲水能力越強，土壤中能夠保存的雨量也越多。

在降雨開始后6及12 h的分析結(jié)果與3 h的情況相似，差異是降雨強度對土壤入滲量的影響略有不同。12 h時，在30 mm/h雨強的條件下，土壤入滲量增幅有所減緩，推測其原因是由于在相同的降雨歷時內(nèi)，30 mm/h的雨強條件使表層土壤更快達到飽和，導(dǎo)致產(chǎn)流提前發(fā)生，土壤入滲速率下降，所以趨勢線有所改變。

3 不同下墊面下土壤含水率變化特征分析

覆膜和起壟都會將降雨過程由一維入滲轉(zhuǎn)換為二維入滲，因此對本試驗中不同下墊面下的土壤含水量變化特征進行分析，由于土壤初始含水率的差異對降雨入滲過程影響較小，后面分析不同下墊面對土壤含水量的變化主要分析下墊面處理及雨強對入滲過程的影響。

3.1 平整裸地條件下的垂向含水率變化

為了進行對比分析，首先分析平整裸地下墊面條件下各組試驗的土壤含水率變化情況。選取雨前和降雨開始后3、6、12以及24 h的土壤含水率監(jiān)測數(shù)據(jù)進行分析；利用3次試驗過程中5個時間點的土壤含水率數(shù)據(jù)繪制散點圖，如圖3所示。

觀察上述3個圖，可以看出在降雨過程中土壤含水率有明顯的增加，對于試驗A1B1C1和A1B2C2，10 cm處土壤含水率在降雨開始后6 h達到基本飽和，后期雖有增加，但增幅不大，后面的降雨主要進行垂向運動，從圖3可以清晰看出入滲濕潤鋒的運移過程。A1B3C3試驗的雨強為30 mm/h，從降雨開始，各層土壤含水量由上到下快速增加，降雨開始后6 h之后的土壤含水率變化不大。表1中10 cm處的土壤飽和含水率為0.427，對比圖3(c)可知10 cm處的土壤已經(jīng)達到飽和含水率，因此含水率不再變化。由降雨開始后24 h的土壤含水量變化情況可以看出，30 cm以上土壤含水率都有明顯減少，這是因為降雨開始后24 h，降雨已經(jīng)停止一段時間，一部分水量向下繼續(xù)運移使得表層的土壤含水量逐漸減少。另外表層土壤含水量存在一定量的蒸發(fā)，這也是表層土壤含水量減少的另一個原因。

圖3 A1B1C1、A1B2C2、A1B3C3含水率變化曲線

圖3(a)中90 cm以下的土壤含水率一直沒有變化，主要原因是由于雨強小，雨量少，土壤入滲沒有達到深層；圖3(b)在降雨開始后6 h，70 cm以下土壤含水率變化不明顯，直到降雨開始后12 h，70 cm以下土壤含水率才開始逐漸增加，說明只有在雨量充分的情況下，土壤水才能下滲至深層；在降雨開始后24 h，此時降雨已經(jīng)結(jié)束12 h，各層土壤含水率都有所減小，說明土壤水仍繼續(xù)向120 cm以下入滲，甚至有深層滲漏產(chǎn)生。

圖3(b)與圖3(c)對應(yīng)的兩個試驗的降雨歷時相同，但降雨強度分別為20和30 mm/h，但圖3(c)在90 cm以下的土壤含水率增加不如圖3(b)的明顯，推測其原因可能是圖6的試驗雨強大，土壤表層在短時間內(nèi)迅速達到飽和，提早產(chǎn)流，此時降雨強度>土壤下滲速率，使得土壤中的入滲量反而減少，這也表明只有在一定范圍內(nèi)，雨強與土壤入滲量呈正相關(guān)。

3.2 不覆膜起壟條件下的垂向含水率變化

不覆膜起壟試驗進行了3組試驗，由于壟溝的集雨作用，除了降雨的垂直入滲之外，還存在水量的側(cè)向運移，因此布置了3個監(jiān)測剖面，分別是臺中、臺邊和溝中。分別選取A2B3C1、A2B1C2、A2B2C33組試驗的3個監(jiān)測剖面在雨前、降雨開始后3、6、12及24 h的土壤含水率數(shù)據(jù)，利用這些土壤含水率數(shù)據(jù)繪制散點圖，見圖4～圖6。

圖4 A2B3C1不同含水率變化曲線圖

圖5 A2B1C2不同含水率變化曲線

圖6 A2B2C3不同剖面含水率變化曲線

首先分析A2B3C13個剖面的土壤含水量隨著降雨入滲的變化過程。對比上述3個圖的差異，在3 h時，溝中的含水量增幅最大，臺中次之，臺邊最少。分析其主要原因，臺邊由于坡度的影響，雨量會有不同程度向溝中匯集的過程，因此臺邊最小，溝中最大，臺中是正常的垂直入滲過程，不受壟溝邊皮坡度影響。其后其他時間的土壤含水量變化趨勢基本與3 h一致，都有不同程度的增加，但都顯示其增加幅度為溝中>臺中>臺邊。由此可以看出壟溝有明顯的匯集降雨的作用。

試驗A2B1C2與A2B3C1相比，降雨強度由10 mm/h增加到20 mm/h，分析3個土壤含水量變化圖5可知，隨著雨強的增加，土壤水分運移的過程及濕潤鋒運移更加明顯。3個剖面進行對比可以看出它們的土壤含水量的增加趨勢基本一致，但土壤含水量的增幅明顯不同。在3 h時，土壤含水量增幅變化由大到小分別為溝中、臺中和臺邊。到6 h及其以后，土壤含水量增幅由大到小變?yōu)闇现?、臺邊、臺中，分析原因主要是由于本次降雨雨強較大，壟溝的集雨作用不僅使溝中剖面的土壤含水量增加最大，也使得臺邊剖面由于水分的側(cè)向運移得到了一定的補給，因此導(dǎo)致后期臺邊含水量增幅大于臺中。

對于A2B2C3試驗，降雨強度增加到30 mm/h，對比A2B2C3試驗3個剖面的土壤含水量的變化過程，不論是土壤含水量的增加過程還是增加幅度，都顯示極強的一致性，由此可以看出，在雨強及雨量足夠大以后，壟溝的集雨作用不再顯著。

通過上述分析可以看出，在起壟下墊面處理下，雨強對壟溝的集雨作用有明顯的影響，在雨強很小時，即使在臺邊也沒有微徑流的產(chǎn)生，壟溝的集雨作用基本沒有；雨強小幅增加以后，臺邊會有少量的微徑流，壟溝的集雨作用開始顯現(xiàn)，溝中的入滲量明顯大于其他剖面，但由于雨強還不夠大，溝中雨量增幅不足以進行明顯的側(cè)向運移，因此也不能通過溝內(nèi)的雨水補給其他剖面；如果雨強繼續(xù)增大，壟溝集雨量大幅增加，不僅可以增加溝內(nèi)的垂直入滲，也可以側(cè)向補給臺邊剖面下方的土壤；當(dāng)雨量繼續(xù)增加，整個田間都有積水時，壟溝的集雨作用會得到削弱，各個剖面的降雨入滲量差異不明顯。

3.3 覆膜起壟下墊面處理的垂向含水率變化

覆膜起壟下墊面對應(yīng)的3組試驗分別是A3B2C1、A3B3C2和A3B1C3，3個監(jiān)測剖面分別是溝中、膜邊和膜中。同樣利用不同剖面監(jiān)測的雨前及雨后土壤含水量繪制散點圖，見圖7～圖9，來分析土壤水的運移情況。

A3B2C1試驗的降雨強度為10 mm/h，分析該組實驗3個剖面的土壤含水量變化過程。在降雨開始3 h時，膜邊及溝中的土壤含水量在30 cm以上的范圍都有增加，其中膜邊增幅最大；降雨開始6 h后，膜中剖面的表層土壤含水量才有小幅增加，其他2個剖面土壤含水量有顯著增加；24 h時，3個剖面土壤含水量分布趨勢的大小基本一致，即3個剖面的土壤含水量分布基本達到一致化。由此可以看出該組試驗在初級階段由于膜上水分的匯集，膜邊的土壤含水量增加最快，隨著時間的推移，膜邊的土壤含水量不僅向下運移，也向膜中進行側(cè)向的水平運移，到24 h時各剖面達到一致。

圖7 A3B2C1不同剖面含水率變化曲線

圖8 A3B3C2不同剖面含水率變化曲線圖

圖9 A3B1C3不同剖面含水率變化曲線圖

A3B3C2試驗的降雨強度為20 mm/h，由于雨強的增加，各剖面的降雨入滲過程有一定的不同。3 h時，溝中及膜邊土壤含水量已經(jīng)開始增加，溝中增幅最大，膜中土壤含水量沒有變化；6 h時也是溝中的土壤含水量增幅最大，膜邊次之，膜中也有一定程度的增加，其后此趨勢一直延續(xù)；到24 h時，各剖面土壤水分經(jīng)過垂向運移和水平運移基本達到一致。

對于試驗A3B1C3，該組試驗的雨強提高到30 mm/h，降雨開始3 h時，各剖面的土壤含水量的變化趨勢及增量都基本一致，3個剖面沒有明顯差異，且這種一致性延續(xù)到12 h；12 h后，膜邊及膜中的土壤含水量變化趨勢基本一致，但溝中部分的雨水由于不直接進行側(cè)向運移以補充膜中區(qū)域的含水量，因此溝中的土壤含水量最高。

通過覆膜起壟下墊面的3組試驗分析可以看出，不同的降雨條件下，覆膜及起壟對降雨入滲過程的影響也不同。在雨強很小時，由于膜的不透水作用，膜上雨量會運移到膜邊進行入滲；雨強稍微增加后，膜上水分除了向膜邊匯流外，還會向溝中匯集，膜邊的雨水通過側(cè)向運移到達膜中，因此溝中的降雨入滲量總體最高；隨著雨強的持續(xù)增加，膜的不透水作用繼續(xù)存在，而壟溝的集雨作用則不再顯著，覆膜一方面導(dǎo)致膜上水分向膜邊及溝中匯集，另一方面，膜邊的水分會向膜中下方土壤做水平運移，土壤水分的二維入滲過程更加顯著，該過程一直持續(xù)到3個剖面土壤含水量差異不大。

4 結(jié) 語

為了分析滴灌的覆膜及起壟對降雨入滲過程的影響，本文開展了不同下墊面處理下的人工降雨測坑土壤水分變化監(jiān)測試驗研究，通過試驗數(shù)據(jù)分析可以得到以下結(jié)論。

(1)下墊面處理、降雨強度及初始土壤含水率3個影響因素對降雨入滲量的影響從大到小分別為降雨強度、下墊面處理及初始土壤含水率。在降雨強度為20 mm/h時，入滲量最大，當(dāng)降雨強度繼續(xù)增加至30 mm/h時，淺層土壤含水量迅速達到飽和，地表徑流很快發(fā)生，導(dǎo)致入滲量減小。土壤初始含水率對不同下墊面的降雨入滲特征影響不顯著。

(2)在起壟條件下，降雨強度在一定范圍內(nèi)時，壟溝的集雨作用隨著雨強的增加而增加，降雨的二維入滲過程比較明顯；但雨強達到30 mm/h時，壟溝的集雨作用則不再顯著。

(3)在起壟及覆膜起壟條件下，在一定的雨強范圍內(nèi)，膜上水分由于膜的不透水作用及壟溝的集雨作用，分別是膜邊和溝中的降雨入滲量最大；雨強繼續(xù)增加時，壟溝的集雨作用被削弱，但膜的匯流作用仍然存在，覆膜區(qū)域下方的土壤水分都是通過膜邊的土壤水分經(jīng)過側(cè)向運移達到，土壤水分的二維入滲過程十分顯著。

[1] 徐克輝,盛 平. 應(yīng)用Mein-Larson入滲模型計算低洼地區(qū)徑流系數(shù)[J]. 水文學(xué)報,1992,(5):32-34.

[2] 趙西寧. 黃土坡耕地土壤入滲規(guī)律研究[D]. 陜西楊凌:西北農(nóng)林科技大學(xué),2002.

[3] 陳麗華. 黃土地區(qū)水土保持林地土壤入滲規(guī)律的研究[J]. 北京林業(yè)大學(xué)學(xué)報,1995,17(3):51-55.

[4] 吳欽孝,韓 冰. 黃土丘陵區(qū)小流域土壤水分入滲特征研究[J]. 中國水土保持科學(xué),2004,2(2):1-5.

[5] 孫元元. 降雨入滲條件下土壤非飽和特性研究[D]. 南京：南京水利科學(xué)研究院,2006.

[6] 李建興,何丙輝. 生物埂護坡下不同土地利用方式對土壤滲透性的影響[J]. 水土保持學(xué)報，2012,26(6):243-248,253.

[7] 宇 宙，王 勇，羅迪漢,等. 膜下滴灌玉米需水規(guī)律及優(yōu)化灌溉制度研究[J]. 節(jié)水灌溉，2015,(4):10-13,18.

[8] 趙靖丹,李瑞平,史海濱,等. 滴灌條件下地膜覆蓋對玉米田間土壤水熱效應(yīng)的影響[J].節(jié)水灌溉,2016,(1):6-9,15.