劉珂瑤，吳普特，朱德蘭，蔡耀輝，宋小林，王建平

(1. 西北農林科技大學水利與建筑工程學院，陜西 楊凌 712100； 2. 西北農林科技大學中國旱區(qū)節(jié)水農業(yè)研究院，陜西 楊凌 712100； 3. 西北農林科技大學水土保持研究所，陜西 楊凌 712100； 4. 中國科學院水利部水土保持研究所，陜西 楊凌 712100； 5. 西北農林科技大學園藝學院，陜西 楊凌 712100； 6. 延安市寶塔區(qū)果業(yè)技術推廣和營銷服務中心，陜西 延安 716000)

涌泉根灌是一種針對高耗水經濟林果進行局部、大流量灌溉的地下灌溉技術[1-3].該技術通過微管-灌水器將水肥直接輸送至樹木根區(qū)，形成適宜其生長的水、肥、氣、熱條件[4-7].灌水器一般采用迷宮式流道，流道較長、過水斷面較大且抗堵塞能力較強[8-11]，外部設置有PVC或多孔材料套管，可有效解決負壓吸泥和根系入侵導致的堵塞問題[4].

歷經10余年，涌泉根灌技術在黃土高原應用后取得了良好的節(jié)水、增產效果，提高了作物的水分利用效率，減少了表層土壤水分蒸發(fā)[12-15].大田試驗中，多以蘋果和紅棗為研究對象.ZHONG等[16]發(fā)現對蘋果從開花至結果期進行中度水分虧缺處理可提高果實品質、產量和水分利用效率，與漫灌相比，涌泉根灌的節(jié)水效果更明顯.曾健[17]對果樹進行調虧灌溉，分析了蘋果各生育期的耗水規(guī)律以及土壤含水率對果樹生長和果實產量、品質的影響，建立了陜北山地涌泉根灌蘋果樹灌溉制度評價體系.強敏敏等[18]發(fā)現利用涌泉根灌技術對棗樹進行輕度和中度水分虧缺處理可有效抑制棗吊和新梢生長，促進果實生長.車銀偉等[19]通過試驗發(fā)現每株灌水80 L為應用涌泉根灌技術時棗樹的最優(yōu)灌水定額.但在實際應用中，由于灌水技術參數設置不合理導致的表層蒸發(fā)損失、深層滲漏問題依然存在.因此，通過明晰涌泉根灌入滲特性，確定其應用參數，將為該技術的大面積應用提供理論和技術支撐.

近年來，相關學者將室內試驗與數值模擬相結合對涌泉根灌土壤水分入滲特性進行了深入研究.牛文全等[20]研究發(fā)現初始含水率對涌泉根灌濕潤體大小有顯著影響，初始含水率越大，濕潤體體積越大.代智光[21]建立了涌泉根灌土壤水分入滲模型，并利用Hydrus-2D軟件進行求解，發(fā)現模擬值與實測值具有較好的一致性.李耀剛等[22]運用Hydrus-3D軟件對涌泉根灌土壤水分入滲進行模擬，發(fā)現灌水量相同時，隨灌水器流量增加，濕潤體內高含水區(qū)域的面積明顯增大.現有研究大多僅揭示了涌泉根灌技術在單因素影響下的入滲特性，缺乏對其在套管和灌水器參數交互影響下入滲特性的研究.

在實際應用中，涌泉根灌配套套管多采用PVC打孔管，人工打孔費時費力，質量控制難度較大.研究中也有學者將PP棉濾芯作為配套套管，其是以聚丙烯樹脂纖維纏繞粘結而成的多孔滲水管，具有可工業(yè)化大批量生產、價格低廉和質量穩(wěn)定的優(yōu)點.陳俊英等[4]應用PVC和PP棉2種材料制成涌泉根灌套管，發(fā)現2種套管在土壤表面形成的水平濕潤半徑與時間的自然對數呈線性關系，并建立了由表面水平濕潤半徑反演最大水平濕潤半徑和濕潤深度的預測模型.但涌泉根灌入滲特性不僅與套管參數直接相關，更受灌水器流量和套管參數交互作用的影響.

目前，選用PP棉濾芯作為涌泉根灌套管，探究其入滲特性對流量和套管長度反饋機制的研究尚未見報道.因此，文中研究選用PP棉濾芯作為涌泉根灌配套套管進行土壤入滲試驗，設置不同的套管長度和灌水器流量，揭示兩者交互作用下土壤水分入滲特性，明確灌水器的出流機制，提出技術應用參數(流量和套管長度)，為該技術的實際應用提供技術參考.

1 材料與方法

1.1 試驗裝置

涌泉根灌入滲試驗在西北農林科技大學旱區(qū)節(jié)水農業(yè)研究院灌溉水力學試驗大廳進行，試驗裝置以及灌水器和PP棉套管細部結構如圖1所示.試驗裝置由供水裝置、灌水裝置、土箱和土壤水分傳感器4部分組成.供水裝置由水箱、水泵和管道等組成.灌水裝置由涌泉根灌灌水器和孔徑為5 μm的PP棉套管組成.灌水器埋置于土箱一角，以模擬實際大田中1/4的土壤水分運移情況，進水口距土壤表層10 cm.矮化蘋果樹的根系主要集中在0～80 cm土層深度，上層根系分布密集，下層稀疏[23].為確定濕潤體區(qū)域與果樹根系分布范圍的匹配效果，選擇尺寸為45 cm×45 cm×75 cm(長×寬×高)的有機玻璃土箱.

圖1 試驗裝置以及灌水器和PP棉套管細部結構示意圖

1.2 試驗土壤

試驗土壤取自陜西省延安市寶塔區(qū)，土壤類型為黃綿土.試驗前,土壤經自然風干后碾壓、粉碎過2 mm篩后混合均勻制成試驗土樣備用.土壤顆粒組成由馬爾文激光粒度分析儀MS2000測定，其中黏粒、粉粒和砂粒體積分數分別為16.15%，28.36%和55.49%，土壤質地為砂壤土.將試驗土樣按1.35 g/cm3的容重裝入土箱，每層10 cm，層間打毛.試驗過程中用塑料薄膜覆蓋土壤表面，以防止土壤水分蒸發(fā).

1.3 試驗方法與觀測內容

試驗過程中，各處理的總灌水量相同，計算式為

m=γzp(θmax-θmin)，

(1)

式中:m為灌水定額，mm；γ為土壤容重，取1.35 g/cm3；z為計劃濕潤層深度，取0.75 m；p為濕潤比，取0.35；θmax，θmin分別為土壤體積含水率上、下限(占土體積百分比)，分別取41.08%，1.89%.

經計算，總灌水量為7 L.試驗設置了3個不同的套管長度I(20，40，60 cm)和2個灌水器流量Q[18,19,24](4，12 L/h)，共6個處理.試驗處理見表1.

表1 試驗處理表

由于試驗時將套管和灌水器置于土箱一角，模擬大田試驗1/4濕潤體情況，故取大田實際使用流量的1/4，分別為1和3 L/h進行土箱試驗.

土壤經自然沉降24 h后開始試驗.入滲開始后,采用秒表計時，定時觀測所形成濕潤體的水平濕潤半徑和垂直向上入滲距離.按照先密后疏的原則，在有機玻璃外壁描出不同時刻所對應的濕潤體形狀，并用鋼尺測量濕潤鋒距離.隨后應用美國Decagon公司的EM50土壤水分監(jiān)測系統(tǒng)監(jiān)測土壤含水率的動態(tài)變化.

2 試驗結果與分析

2.1 對濕潤體形狀的影響

圖2為不同處理水平(L)和垂直向上(R)濕潤鋒隨灌溉時間(t)的變化情況：變化趨勢基本相同，先快速增加，而后速率逐漸減慢.入滲初期，基質吸力較大，灌水器流量為4和12 L/h時，分別在0～50和0～20 min時間段內濕潤鋒的變化速率最快.隨入滲時間增加，灌水器周圍一定范圍內的土壤體積含水率接近體積飽和含水率，毛管懸著水達到最大，基質吸力減小，使得濕潤鋒變化速率明顯下降.

灌水結束時，相同灌水器流量下，套管長度為40 cm所形成的水平和垂直向上濕潤鋒均最大.套管長度相同時，灌水器流量為4 L/h所形成的水平濕潤鋒大于12 L/h的；垂直向上濕潤鋒則相反.在套管長度和灌水器流量的交互影響下，水平和垂直向上方向最有利于土壤水分運移的配套組合分別是處理2和5.

圖2 水平濕潤鋒和垂直向上濕潤鋒運移距離與時間關系圖

研究表明[4]L和R與灌溉時間t之間的關系可用冪函數表示，即

(2)

擬合參數a,b,c,d見表2，相關系數R2均不低于0.97，各參數擬合結果較好，能較好說明涌泉根灌土壤入滲濕潤鋒與入滲時間的關系.

表2 不同套管長度和灌水器流量下濕潤鋒運移距離與入滲時間關系的擬合結果

灌水結束后土壤水分仍在繼續(xù)運移.與灌水過程中的運移速率相比，此時的運移速率大大降低.表3為再分布過程中水平和垂直向上方向上濕潤鋒運移距離.由表可知，從灌水結束時至灌水結束后24 h，灌水器流量為4 L/h的水平和垂直向上方向濕潤鋒的變化量均大于灌水器流量為12 L/h的.再分布過程中，處理2的濕潤鋒運移距離變化最明顯，水平和豎直向上方向濕潤鋒的變化量分別為4.6和4.3 cm.

表3 再分布過程中水平和垂直向上方向上濕潤鋒運移距離

2.2 對濕潤體內土壤水分分布的影響

2.2.1 灌水結束時土壤水分分布規(guī)律

圖3為灌水結束時不同處理所形成濕潤體剖面的土壤水分分布圖，圖中θ為土壤體積含水率；x,y分別為水平和垂直位置.相同套管長度、不同灌水器流量所形成濕潤體的形狀以及土壤水分分布規(guī)律相似，土壤體積含水率均從最大值處向四周遞減，分布曲線由疏到密變化，濕潤體高含水區(qū)出現在灌水器下方偏右位置.灌水器流量相同時，套管長度為40 cm所形成濕潤體的體積最大且水分分布更加均勻.灌水器流量為12 L/h的濕潤體體積含水率最大值均大于灌水器流量為4 L/h的.在流量和套管長度的交互作用下，處理1所形成濕潤體高含水區(qū)域范圍以及含水率的最大值均小于其他處理.處理6由于土箱尺寸的限制，濕潤體內土壤水分積聚現象最明顯，含水率的最大值明顯大于其他處理.處理5所形成濕潤體的體積最大，水分分布更加均勻，最有利于土壤水分運移.

圖3 不同處理灌水結束時土壤剖面水分分布圖

2.2.2 灌水結束24 h后土壤水分再分布規(guī)律

圖4為灌水結束24 h后濕潤體剖面的土壤水分分布圖.

圖4 不同處理灌水結束24 h后土壤剖面水分分布圖

再分布過程中，濕潤體高含水區(qū)的含水率減小，土壤水分分布趨于均勻，濕潤體體積明顯增大.與大流量相比，小流量灌水歷時長，基質勢對其水分擴散的影響大，灌水過程中已開始進行水量平衡，因此再分布過程中濕潤體內體積含水率變化程度較小.在流量和套管長度交互作用下經過24 h再分布，處理6土壤水分消退現象最明顯.處理5濕潤體體積變化最大，水分分布最均勻.

3.3 對套管內積水深度和入滲速率的影響

灌水過程中，灌水器流量為4 L/h的處理，其入滲速率始終小于或等于土壤入滲速率，套管內并未形成積水.因此，此處僅討論灌水器流量為12 L/h的情況.

圖5為灌水器流量為12 L/h時套管內積水深度h和入滲速率v隨時間的變化圖.0～30 min套管內的積水深度快速增加，30～60 min增速減緩，趨于穩(wěn)定，80 min后基本為定值(圖5a).套管長度為40 cm時，套管內積水深度最大，但僅較其他2個處理高2 cm.由圖5b可以看出,20～30 min入滲速率迅速降低，30～60 min入滲速率緩慢降低，60 min后基本穩(wěn)定.穩(wěn)定時，套管長度40 cm的入滲速率最小為0.225 cm/min.

3 討 論

垂直向上方向受基質吸力與反方向的重力勢作用，隨套管內積水深度增加，垂直向上濕潤鋒與積水上邊界的水分梯度減小，將增大該方向上濕潤鋒的運移速率.因此，灌水結束時灌水器流量為12 L/h時的垂直向上濕潤鋒要大于4 L/h的，這與費良軍等[25]研究結果一致.

濕潤體的位置主要受套管長度影響.套管長度越長，濕潤體越接近底層土壤，且在重力勢和基質勢的共同作用下，高含水率區(qū)域也會發(fā)生下移和右移現象，這與地下滴灌濕潤體的特征有一定區(qū)別[26].文中研究由于受土箱尺寸的限制，套管長度為60 cm所形成的濕潤體主要位于土壤底層，土箱底板與分層土壤中的不透水層類似[27]，可有效阻隔土壤水向底層運移，減少深層滲漏風險，提高灌溉水利用效率.因此，在土壤質地較輕、蒸發(fā)強烈的地區(qū)可將較長的套管和隔離層聯合使用，使有限的水分匯聚至作物根系集中的土層.GUO等[28]發(fā)現通過土壤壓實可減小土壤孔隙度，降低水分入滲速率，有效防止灌溉水向下層滲漏.因此，實際應用中可先將與套管底部接觸的土壤區(qū)域進行壓實處理再放置套管進行填埋，以降低深層滲漏，提高保水能力.

對于同一種土壤，灌水器流量和土壤特性的交叉作用會對其穩(wěn)定入滲速率產生影響，從而使積水深度發(fā)生變化.灌水器流量相同時，由于套管長度不同會導致濕潤體體積和濕潤體內高含水區(qū)域范圍不同，則在土壤水勢作用下會對積水深度產生一定影響.套管長度為40 cm時所形成濕潤體的體積和高含水區(qū)域范圍均最大，使得土壤水勢對套管內積水入滲的抑制作用最明顯.灌水過程中套管內腔并未被水完全充滿，仍存留有少量空氣，對套管頂部和底部進行密封處理后，空氣壓力也會對灌溉水出流產生影響.其中60 cm時壓縮空氣對出流的促進作用最明顯.因此，在土壤水勢和壓縮空氣的共同作用下，套管長度為40 cm時積水深度最大.

套管內積水深度發(fā)生變化時，積水上沿的灌溉水會迅速入滲，使得積水深度降低，基本上維持平衡，使得穩(wěn)定入滲后入滲速率維持相同.但相對于地下滴灌點源入滲，灌水器出口處形成正壓[29]，使得灌水器流量降低，涌泉根灌套管內部的積水顯著提高了土壤入滲的水勢梯度，形成有壓壓差水頭入滲，提高了土壤水分的運移速率，促進土壤水分擴散.當土壤水分入滲達到穩(wěn)定后，套管內積水深度不再變化，表現出典型的柱狀出流面源入滲.這與李耀剛等[30]的初步研究結論一致.

4 結 論

1) 在套管長度和灌水器流量的交互影響下，最有利于土壤水分在水平和垂直向上方向上運移的配套組合：灌水器流量為12 L/h、套管長度為40 cm.

2) 灌水器流量主要影響濕潤體的特征，而套管長度主要影響濕潤體的位置.灌水器流量為12 L/h、套管長度為40 cm時所形成濕潤體的體積更大，水分分布更加均勻，最有利于土壤水分的運移.

3) 涌泉根灌入滲時，套管內部積水會使入滲面增大，提高土壤入滲水勢梯度，形成有壓入滲，從而提高土壤水分的運移速率.流量和土壤特性對積水深度影響較大，套管長度對積水深度影響較小.

4) 根據不同果樹主根系的分布范圍合理確定套管長度和埋深，選擇適宜的灌水器流量，并進行土壤壓實處理，可降低深層滲漏和表層蒸發(fā)損失，提高水分利用效率和保水能力.