限制入滲域下滴灌土壤水分時空變異特征研究

寧莎莎，徐利崗，張 林，王懷博

（1.西北農林科技大學水利與建筑工程學院，陜西咸陽712100；2.寧夏水利科學研究院，銀川750021；3.西北農林科技大學水土保持研究所，陜西咸陽712100；4.寧夏回族自治區水利科學研究院，銀川750021）

0 引 言

根域限制栽培是根據作物根系的分布范圍開挖相應的栽植溝，在溝渠邊壁及底部鋪墊隔膜材料，并在底部預留透氣排水的空隙或鋪設波紋管，栽植溝內填充秸稈、農家肥及土壤混合物后種植果樹的一種栽培方式。根域限制栽培方式能抑制作物根系生長、抑制地上部分營養生長和矮化植株，并且采用滴灌點源局部水肥供給，以提高水肥利用效率。方金豹等人研究了不同根域限制方式對桃樹生長的影響，得出桃樹適宜的根域限制方式和根域體積[1]；王世平等人綜合了在根域限制條件下葡萄、蘋果、桃、櫻桃的研究進展和存在問題，肯定了果樹根域限制的控長促花以及提高果實品質的效果[2]；謝兆森等人從葡萄果實發育動態、果實膨壓及動力學參數來探討根域限制對葡萄果實發育及糖分積累影響的機制,研究了根域限制對葡萄果實發育動態、膨壓及動力學參數的影響[3]。但是限根技術在滴灌枸杞種植中應用的研究較少。土壤濕潤體的形狀、大小及含水率分布等是確定灌溉制度的重要依據[4]，許多專家和學者對滴灌的入滲特性等各方面進行了大量試驗并取得了一定成果。侯立柱等人通過室內模擬試驗得出滴灌對作物根系蓄水保墑效果良好[5]；李道西等人研究了不同滴頭流量對不同質地的土壤水分運動規律的影響[6]。然而當前研究對根域限制條件下滴灌的土壤水分遷移規律鮮有報道。本研究以土箱為工具模擬根域限制立地條件，開展室內試驗，探究在限定土體邊壁條件下，不同滴頭流量點源入滲的土壤水分時空變異性特征，為根域限制栽培作物提供一定的數據參考。

1 試驗材料與方法

1.1 試驗裝置與材料

1.1.1 試驗裝置

試驗在西北農林科技大學中國旱區節水農業研究院灌溉水力試驗大廳進行。試驗裝置由水源及加壓系統、滴灌帶、土箱及土柱、自動土壤水分監測系統組成(圖1)，加壓裝置為可提供穩定壓力的清水泵，其吸程為9 m，功率為0.75 kW；試驗土箱為有機玻璃制作，尺寸為85 cm×60 cm×60 cm，采用一行兩帶的方式鋪設方式布置滴灌帶，在土箱底部及一側邊壁用棚膜鋪蓋(棚膜厚度10～12絲)，另一側不鋪膜作為對照。

1.1.2 試驗材料及準備

試驗土壤取自陜西渭河三級階地，按試驗所需土樣將其風干、碾壓、混合均勻后過2 mm 篩網后盡可能曬干土壤備用(實際含水率以烘干法測定為準)。土壤顆粒組成采用激光粒度分析儀(MS2000 型，馬爾文，英國)測定，按國際制土壤質地分類標準，試供土壤土質為黏壤土，采用環刀法測得容重，烘干法測得田間持水率及飽和含水率，結果如表1所示。按照設計容重1.327 g/cm3進行填土，每層5 cm，共12 層，將兩土層之間進行打毛，使土壤顆粒充分接觸。兩條滴灌帶平行放置在土壤表面(間距20 cm)。

試驗設置不同流量的滴頭。毛管安裝在Φ25 PE 管上，上級管道連接至水泵，水泵出水口安裝壓力表及回流裝置，管網系統尾部設排氣閥(試驗開始前進行排氣處理)，根據所需流量調節水壓。重復試驗依次進行。

表1 試驗土樣基本性質Tab.1 Physical properties of tested soil

1.2 試驗方法及測定內容

1.2.1 試驗設計

為對比分析限制入滲域條件及滴灌流量對點源入滲的土壤水分時間、空間變異性特征的影響，本研究采用初始含水率θ0=2.45%(占干土重)的黏壤土，開展地表滴灌土箱試驗。試驗設置了限制入滲域條件與常規條件，限制入滲域條件的土箱底部鋪設棚膜，常規條件的土箱涂抹凡士林以防止壁面流。試驗采用不同的滴頭流量，分別為3 L/h、4 L/h，壓力水頭為10 m。

1.2.2 檢測指標及方法

(1)試驗過程中監測在土箱正面中心處滴頭在單平面上的濕潤鋒演進過程，根據濕潤鋒的演進過程隨時標記濕潤鋒位置并同時記錄灌水時間，以滴頭為中心在土箱表面將用記號筆描摹濕潤線并記錄濕潤峰的坐標。

(2)采用美國Decagon 公司生產的EM50 監測系統實時監測土壤含水率，數據步長1 min，測定深度為10、20、40、60 cm 共5 層，監測水平位置為滴頭正下方、沿滴灌帶方向距離滴頭10、20 cm 處，垂直滴灌帶方向距離滴頭10、20 cm(貼近限制入滲域邊壁)處。

1.3 數據處理

該試驗各處理數據采用Excel2016 進行數據處理，Origin2017、surfer8.0繪圖。

2 結果與分析

2.1 限制入滲域條件對濕潤體形狀及范圍的影響

濕潤體的形狀大小與土壤質地、土壤容積密度、土壤初始含水率、灌水時間、灌水量和滴頭流量[7]等眾多因素有關[8]，本試驗各處理土壤的質地、容重、初始含水率均相同。按照設計分別采用滴頭流量為3 L/h和4 L/h的灌水器進行滴灌，實時監測土壤濕潤鋒的演進過程，并分別繪制不同灌水時間下限制入滲域條件與常規條件的濕潤體變化情況(圖2)。圖中原點為滴頭位置，濕潤峰水平運移距離為R，濕潤峰垂向運移距離為H。

從圖2左右對比可以看出，同樣流量、同樣入滲時長，限制入滲域條件與常規條件下，前期的濕潤體的形狀及大小差別不大，濕潤峰在同一水平面垂直于滴灌帶方向和平行于滴灌帶方向及其他方向的運移距離基本相同，而入滲后期，限制入滲域條件下濕潤體呈不規則柱形，前后形變較大。分析濕潤峰數據推斷水平濕潤峰運移距離與垂直濕潤峰運移距離存在線性關系(表2)，且判定系數(R2)均大于0.96。各處理的回歸直線的斜率均小于1，這表明在試驗開始120 min 內，各個處理濕潤峰水平運移距離均大于同時間點的垂向濕潤峰運移距離。綜合濕潤峰水平及垂向運移特性，濕潤體呈現出以滴頭為中心的半橢球體。

2.2 限制入滲域條件對濕潤峰運移的影響

流量3 L/h、4 L/h 的限制入滲域條件與常規條件下邊壁濕潤峰運移距離與時間的關系見圖3。

當土壤水運移至邊壁位置后(入滲時長＞120 min)，限制入滲域條件與常規條件對水分入滲的影響非常明顯。如圖3所示，在限制入滲域條件下鋪膜的邊壁水分入滲明顯慢于常規條件下的入滲，入滲速率下降幅度在2.30%～40.38%，說明在底部鋪膜后對于灌水后期水分的邊壁入滲及其再分布影響很大，明顯降低了水分下滲速率，進而在灌水停止后水分再分布過程中水分深層滲漏減少。對不同處理邊壁土壤水運移與時間數據進行分析得出二者符合冪函數關系，故濕潤峰運移距離與時間關系的方程可設為：

表2 水平濕潤峰運移距離與垂直濕潤峰運移距離的關系Tab.2 Relationship between level wetting front and vertical downward wetting front

式中：Y為邊壁濕潤峰的運移距離；t為灌水時間；a、b、y0、t0為擬合參數。

根據試驗實測數據對邊壁濕潤峰運移距離與時間關系公式(1)進行擬合，擬合結果如表3所示。

表3 不同限制條件的邊壁濕潤峰運移與時間關系的擬合參數Tab.3 Fitting parameters of boundary wetting front migration distance with time under different emitter discharges

根據表3擬合得到的濕潤峰運移的冪函數參數，可預測各時刻濕潤峰運移情況，對比結果如表4所示。根據表4可知，由式(1)得到的擬合值與邊壁濕潤峰運移距離的實測值的誤差均在2.09%～5.95%，說明該模型能較好的預測不同時間邊壁濕潤峰運移距離變化。

表4 濕潤峰運移距離實測值與預測值對比Tab.4 Comparison between predicted and measured wetting front migration distance

2.3 限制入滲域條件對土壤水空間變化的影響

限制入滲域的處理(兩側及底部鋪膜)影響了土壤水入滲特征，也導致了各土層土壤含水率空間分布的異質性。繪制濕潤鋒到達土箱邊壁(入滲時長120 min)前后以滴頭為原點的縱剖面土壤水空間分布圖(圖4)。從圖4（a）～圖4（c）可以看出灌溉水到達邊壁之前含水率等值線以滴頭位置為圓心，呈1/4 橢圓形，并隨入滲時間增加向外側擴散，土壤含水率變幅0.08～0.094m3/m3；灌水超過120 min 后的土壤含水率變化見圖4（d）～圖4（f），從圖4可以看出限制入滲域邊壁對水分的滯留現象顯著，且對含水率等值線形狀即土層水分空間分布影響較大。土壤含水率等值線外側明顯向里凹陷，說明邊壁鋪膜能有效滯留水分，減緩下滲且阻止水分的水平擴散及深層滲漏，將有限的灌溉水限制在有限的空間范圍，提高作物對水分的利用。

2.4 滴頭流量對濕潤峰運移規律的影響

2.4.1 流量對濕潤峰水平運移和垂直運移規律的影響

以滴頭位置為原點，分別作出各處理灌溉水到達邊壁之前的濕潤峰水平運移距離與時間的關系圖(圖5)。從圖5可以看出，不同流量濕潤峰入滲初期(灌水40 min 以內)在水平方向和垂直方向的運移距離的增長趨勢相同，前期增長較快，后期逐漸趨于穩定狀態。初期變幅分別為0.15～0.89 cm/min 和0.15～1.1 cm/min；灌水40 min 以后，流量為4 L/h 的濕潤峰水平運移距離明顯大于3 L/h的，流量為3 L/h的濕潤峰水垂向運移距離大于流量為4 L/h的。

由圖5可知，不同流量下的濕潤峰水平與垂直運移距離與時間為冪函數關系，擬合參數見表5，擬合結果的判定系數(R2)均大于0.98，說明該擬合結果能夠較好的反映濕潤峰運移距離與時間的關系。流量為4 L/h 的水平濕潤峰運移指數均大于流量為3 L/h 的水平濕潤峰運移指數，最大為0.273 7，對于濕潤峰運移系數相差不大，在5.6～5.9 之間，最大值出現在流量為4 L/h的處理中為5.977 5；且流量為4 L/h的垂直濕潤峰運移距離的指數均大于流量為3 L/h 的，最大為0.460 2；而流量為3 L/h 的垂直濕潤峰運移距離的系數要大于流量為4 L/h 的，最大為2.755 6。

2.4.2 流量對邊壁濕潤峰垂向運移規律的影響

不同流量下在限制入滲域邊壁的濕潤峰運移變化情況見圖6，從圖6可以看出灌水120 min 時流量為4 L/h 的濕潤峰最先移動到限制入滲域的邊壁，流量為3 L/h 的濕潤峰在隨后的10～20 min 內到達限制入滲域邊壁。根據圖6可知各流量下濕潤峰變化趨勢大致相同，濕潤峰初到邊壁時，其移動速率最大，且隨時間推移該移動速率逐漸減小，后期在相對穩定的最小值上下波動。流量4 L/h 時濕潤峰運移速率約為0.138 4 cm/min，流量為3 L/h 時濕潤峰運移速率約為0.119 6 cm/min，顯然流量為4 L/h的濕潤峰運移速率大于流量為3 L/h的，相差約15.70%，說明較大流量條件下水分下滲更快，由于灌水量是一定的，所以流量為4 L/h 時，在灌水期間保證水分不會深層滲漏的情況下，相比流量3 L/h 的情況可以更快地下滲到作物可吸收的土層范圍，可以更快達到后期水分再分布穩定狀態。

3 討論

3.1 限制入滲域條件對濕潤體形狀及土壤含水量分布的影響

在灌水初期，限制入滲域條件對濕潤體形狀及大小影響非常小，濕潤體形狀近似為半橢球體[9-11]，且大小相近，這是因為土壤水分在各個水平方向入滲較為均勻[12]，而垂直入滲速率主要取決于兩個方面：灌水速率和土壤的入滲能力[13]，明顯灌水速率大于滲水能力，所以垂直入滲速率等于土壤的滲水能力，而表面形成積水，促進了水分水平入滲，所以土壤水分的水平入滲大于垂直入滲，在剖面圖中呈現出半橢圓形，這與張振華等人研究的地表滴灌入滲深度[14]得出的結論明顯不一致，這是因為兩者在土質、土壤入滲速率、試驗條件不同，導致土壤入滲能力不同，因而影響到試驗結果不同。在灌水后期，由于底部鋪棚膜制造出限制區域，水分入滲在限制區域內入滲，使濕潤體形成不規則的柱狀。

限制入滲域條件影響了土壤水入滲，進而影響到土壤水分的空間分布。在灌溉水到達邊壁之前，限制入滲域條件對土壤水分布無影響，在灌溉水運移到限制入滲域邊壁后，水分停止水平入滲，被限制在一定范圍之內，之后僅有垂直下滲。限制入滲域條件不僅改變了濕潤體形狀，還滯留了水分。停止灌水后，在土壤水分再分布階段，限制入滲域條件使水分僅在棚膜形成的區域以內擴散，致使水分不至于深層滲漏，在大田耕作中能有效提高灌溉水的利用效率，達到節約灌溉用水的目的。

表5 不同流量的水平與垂直濕潤峰運移與時間關系的擬合參數Tab.5 Fitting parameters of vertical level wetting front and vertical downward wetting front migration distance with time under different flow rates

3.2 限制入滲域條件對濕潤峰運移規律的影響

限制入滲域條件對濕潤峰的運移的影響顯著。雖然在水平方向上，由于各個處理土壤和滴灌各方面條件相同，在灌水初期(＜120 min)限制入滲域條件對水平濕潤峰無影響，但是灌水時間大于120 min 后，土壤水停止水平運移，這是因為限制入滲域條件底部鋪有棚膜阻止了水分水平運移；在垂直方向上，無論較大流量還是較小流量，相比于常規條件，限制入滲域條件能有效減緩灌溉水在邊壁的入滲，入滲速率要明顯小于常規條件的入滲速率，能夠將水分長時間保存在植物可吸收范圍內，減少水分深層滲漏。可能是因為棚膜主要成分時高分子聚合物，其本身具有疏水性，所以棚膜的疏水性對水分下滲產生了阻力。相比常規條件，限制入滲域條件下邊壁土壤水不僅受到土壤基質吸力垂直向下的分力和自身重力，而且垂直向上方向還受到棚膜對它的阻力，導致土壤水下滲減慢，有效減少水分深層滲漏。

3.3 滴頭流量對濕潤峰運移規律的影響

滴頭流量對濕潤峰的影響主要表現在灌水后期。雖然在灌水40 min 之前，由于入滲速率主要與土壤入滲能力有關，各處理入滲速率相近，使各處理水平濕潤峰和垂直濕潤峰運移距離相差較小，但是在灌水40～120 min 時，由于較大流量產生積水大于較小流量，導致較大流量的水平濕潤峰運移距離大于較小流量的，這與李思明、康紹忠[15]等人研究結果一致。對于限制入滲域條件下邊壁的土壤水分垂直入滲，由于此時土壤水下滲受到重力、基質勢的作用，而大流量處理的濕潤峰水平運移速率快，致使邊壁濕潤鋒處土壤水勢梯度較大，土壤水分運移速度較快[16]，因此較大流量處理在灌水時期，邊界水分運移速率大于小流量處理。較大流量的濕潤峰在垂直方向的運移距離逐漸小于較小流量的垂直濕潤峰的運移距離，主要由于灌水量一定時，較小流量的處理灌水時間長，所以當較大流量處理停止灌水后，較小流量的處理還在繼續灌水，較大流量處理的土壤水分入滲速率逐漸減慢，較小流量處理的土壤水分垂直入滲速率不變。綜上所述，較小流量的條件更利于土壤水分垂直入滲[17]。

4 結 論

(1)限制入滲域條件下，滴灌濕潤體的形狀主要取決于限制入滲域條件，且對土壤含水率的分布有顯著影響，總體而言，限制入滲域條件能有效將土壤水控制在限制區域內，保持土壤濕度。

(2)限制入滲域條件對邊壁濕潤峰影響非常顯著，棚膜能夠阻礙土壤水的下滲，減慢水分下滲速率，但是并不能阻止土壤水分入滲。相比滴頭流量為4 L/h的處理，流量為3 L/h時更有利于土壤水深層入滲。

(3)滴頭流量對灌溉水入滲及其水分再分布有較大的影響。灌水結束后，較小流量條件下，土壤水垂直運移距離更大，并且經過水分再分布，較大流量更能保持土壤水在限制入滲域范圍內，減少水分深層滲漏。