張甲森，劉旭飛，張春華，張 林

（1.西北農林科技大學水利與建筑工程學院，陜西 楊凌 712100；2.西北農林科技大學水土保持研究所，陜西 楊凌 712100）

0 引 言

滲灌是將具有無數微小孔隙的地下多孔滲水管埋入表層土壤以下一定深度處并將灌溉水緩慢滲入到作物根區(qū)土壤的地下節(jié)水灌溉技術[1]，通過將多孔滲水管埋于作物根系土壤，直接向作物根系連續(xù)低流量灌溉，可為根系層土壤營造適宜穩(wěn)定的土壤水分環(huán)境[2,3]，并基于多孔滲水管的滲灌系統可在水流自身重力和土壤毛細管作用力影響，將水分擴散到表層土壤以下，可有效降低水分蒸發(fā)損失，抑制雜草生長，提高水分利用效率[4]。滲灌作為先進的地下節(jié)水灌溉技術，對于現代農業(yè)生產發(fā)展具有十分重要的意義。

塑料滲灌產品大多由不可降解的化學原料制備，在田間應用過程中若回收處理不當，會造成微塑料污染，嚴重時會造成土壤中有毒物質積累，并最終危害國家糧食和人民生命健康安全。農田土壤中的微塑料會顯著破壞土壤細菌群落的多樣性和生態(tài)結構，對土壤理化性質、土壤動植物生長等均產生不利影響，損害土壤健康[5-7]。因此，微孔陶瓷滲灌作為節(jié)水、節(jié)能的新型地下灌溉方式，通過將微孔陶瓷灌水器或管道埋置于地下，進行微壓或零壓灌溉，可持續(xù)性地為農田灌溉提供幫助，緩解環(huán)境負擔[8]。

近年來，大量學者利用微孔陶瓷滲灌開展了田間研究，陶瓷灌水器的主要原料為天然黏土或石英砂，不僅對土壤環(huán)境更為友好，而且會極大簡化當灌溉系統超過使用年限時的后處理步驟，在生菜、番茄、菠菜、枸杞、蘋果等作物上表現出了顯著節(jié)水增產效果[2,9-12]。但目前利用人工拉坯法和機械擠壓法制備的多孔陶瓷滲灌管，存在著坯體壁厚難以掌握，坯體變形程度大等缺點[13,14]。因此，將綠色環(huán)保、適用性強、成型技術易掌握的注漿成型法引入到多孔陶瓷滲灌管道研發(fā)中，可能成為提高陶瓷滲灌管材料和水力性能，保證滲灌系統安全穩(wěn)定高效運行的重要新思路，在生態(tài)環(huán)境和農業(yè)灌溉領域具有廣闊發(fā)展前景。

目前，研究土壤水分入滲過程是評估滲灌管道水力性能的重要方式，土壤水分入滲以管道內外水勢差為主要驅動力[15]，基質勢和壓力勢是影響微孔陶瓷滲灌管道水分入滲的主要因素?，F有的滲灌試驗研究大多從濕潤鋒運移情況、土壤含水量變化、累計入滲量和灌水器出流量角度分析土壤水分運移規(guī)律，Wei等[16]利用土箱試驗深入研究了低壓塑料滲灌管的土壤水分運動規(guī)律并以最優(yōu)根水匹配度優(yōu)化了滲灌管的工作水頭和埋深；陳高聽[17]通過HYDRUS-2D 模型對土壤水分入滲規(guī)律進行動態(tài)模擬，研究了不同工作壓力下的土壤水分運動規(guī)律，但研究對象均為塑料滲灌管。Liu 等[9]試驗表明利用微孔陶瓷灌水器調節(jié)土壤含水率，通過負壓自動控制的方法滿足作物生長溫濕度需求，對加快作物生長和提高作物產量有顯著影響，但其研究也尚未對管道周圍土壤水分時空變化進行系統性研究。Vaghei 等[18]在干旱半干旱地區(qū)采用多孔陶瓷材料管道，研究了10～100 kPa 靜水壓力下管道流量和土壤水分分布情況，僅得出土壤濕潤體形狀呈垂直橢球形，并未對滲灌管道濕潤鋒運移進行深入研究。由此可見，目前大多研究并未采用土箱試驗深入研究微孔陶瓷滲灌管道在土壤水分入滲過程水分空間分布特征，其土壤水分入滲特性尚未明晰。

為此，本研究采用注漿成型法制備微孔陶瓷滲灌管道，通過土箱試驗和數值模擬相結合的方法，研究不同工作水頭條件下微孔陶瓷滲灌管道土壤水分運移規(guī)律。在此基礎上，結合3種典型作物根系分布規(guī)律，以最優(yōu)根水匹配度為優(yōu)化目標，對微孔陶瓷滲灌管道工作水頭進行優(yōu)化，研究結果可為微孔陶瓷滲灌管道的推廣應用提供重要理論依據。

1 材料與方法

1.1 試驗材料

本研究主要以微孔陶瓷滲灌管道為研究對象，該滲灌管采用注漿成型法制備，主要原料為黏土、爐渣、水和硅酸鈉。黏土為陶瓷滲灌管的骨架材料，爐渣為造孔劑、水為黏結劑、硅酸鈉為解凝劑。取陜西楊凌耕作土壤（0.05～0.45 m），將取得土壤經粉碎機粉碎10 min后，過50目（孔徑0.355 mm）篩，得到粒徑范圍為0.18～0.355 mm 的黏土原料。爐渣取自當地陜西渭河某發(fā)電廠，其主要的化學成分是SiO2、Al2O3等氧化物，為粒徑5～20 mm 的灰黑色固體顆粒。硅酸鈉購自陜西咸陽某化工企業(yè)，在陶瓷坯體制備過程中摻入少量（質量比約0.1%）即可改善漿料流動性，保持黏度適中避免出現縮釉等現象，防止?jié){料絮凝。將上述黏土和爐渣按照4∶1 的質量比充分混合，再置于攪拌機中低速攪拌均勻，攪拌30 min 后，按照固體材料：水=66∶34 的質量比，將水放入攪拌機中攪拌30 min，再加入0.1%的解凝劑，得到漿料。而后，將漿料注入石膏模具中，使其形成圓管狀濕坯，再將濕坯置于105 ℃的烘干箱內烘干90 min，再放入高溫電窯以1 085 ℃高溫燒結2 h，最終得到滲灌管道成品。注漿成型法制備而得的微孔陶瓷滲灌管道材料性能優(yōu)良，收縮率、抗彎強度、抗拉強度和抗剪強度分別為4.67%、16.9 MPa、13.4 MPa和8.7 MPa，且此類微孔陶瓷能較好地滿足田間應用要求[19]。

馬氏瓶整體高度1.25 m，安置于0.5 m 的試驗臺上提供恒定水頭；毛管采用管徑2 cm 的PE 管，連接馬氏瓶與微孔陶瓷滲灌管道；有機玻璃土箱主要材質為聚甲基丙烯酸甲酯，規(guī)格為100 cm×100 cm×90 cm（長×寬×高），土箱四周透明且能夠較為明顯地觀察到土壤水分運動情況；本試驗采用的黏壤土來自于楊凌區(qū)及附近的試驗大田，于土壤表層0～15 cm 處選取土樣，土樣經自然晾干后，過篩去除根系等雜質分層填入土箱，每填入5 cm 試驗土樣后壓實，保證每層試驗土樣密實度均勻，填土高度為80 cm，試驗土樣物理特性如表1 所示。微孔陶瓷滲灌管道總長度為100 cm，管道埋深為30 cm，距土箱底部50 cm，于土箱中水平鋪設，管道中心距土箱左右兩側50 cm，管道布置見圖1。

圖1 微孔陶瓷滲灌管道土壤入滲試驗示意圖（單位：cm）Fig.1 Microporous ceramic infiltration pipe soil infiltration test

表1 土壤物理特性Tab.1 Physical properties of the soil

1.2 試驗設計

為研究微孔陶瓷滲灌管道土壤入滲過程，本研究在陜西省楊凌區(qū)現代農業(yè)創(chuàng)新園（108°02′E，34°18′N）的內進行土箱入滲試驗。試驗裝置如圖1 所示，主要由馬氏瓶、止水閥、毛管、有機玻璃土箱、供試土壤、微孔陶瓷滲灌管道和土壤水分采集系統組成（圖1）。微孔陶瓷滲灌管道單寬流量及累計入滲量通過記錄相應時間內馬氏瓶內液面高度差推算[20]。通過觀察記錄有機玻璃土箱外壁濕潤區(qū)域以測量水平、垂直向上和垂直向下濕潤鋒[21]。土壤含水量利用EM50 系列數據采集器（ECH2O 土壤含水量監(jiān)測系統核心部件）進行監(jiān)測，在距土壤表面30 cm處安裝3個ECH2O土壤水分傳感器探針，各傳感器相距6 cm，監(jiān)測時間間隔為5 min，傳感器埋設位置為30 cm 土層深度，待土壤濕潤區(qū)域不再變化時結束試驗，試驗時間48 h。此外，土箱入滲試驗共設置3 組，分別為0.25 m、0.5 m、0.75 m，每組入滲時間均為48 h，該試驗每組同時進行3組重復。

1.3 土壤水分入滲模型

1.3.1 滲透系數

利用馬氏瓶提供的不同工作壓力調節(jié)工作水頭，制得的陶瓷滲灌管道流量額定流量為0.068 2 L/(h·m)（0.2 m 工作水頭），并且流量隨工作水頭的增加而增大（圖2）。采用達西定律可描述不同工作壓力水頭下微孔陶瓷滲灌管道流量：

圖2 不同工作水頭下微孔陶瓷滲灌管道流量Fig.2 The flow rate of microporous ceramic infiltration pipe under different working heads

式中：Q為管道單位寬度流量，cm3/(h·cm)；Ks為水力傳導系數，cm/h；H為提供的工作壓力水頭，m；δ為微孔陶瓷滲灌管道的厚度，cm，δ=1 cm；d為微孔陶瓷滲灌管道的直徑，cm，d=2 cm。工作壓力水頭和管道單寬流量之間的線性函數關系如圖2 所示，管道流量與工作壓力水頭函數關系式見式（2）。

在空氣中進行微孔陶瓷滲灌管道出流試驗，微孔陶瓷滲灌管道流量和工作壓力水頭呈正相關，因此結合公式（1）、（2） 可得出微孔陶瓷滲灌管道水力傳導系數Ks=5.5 × 10-3cm/h。

1.3.2 土壤水分入滲規(guī)律數學模型

本試驗利用HYDRUS-2D模擬軟件，模擬了水從微孔陶瓷滲灌管道向供試土壤滲透過程。三組試驗均將微孔陶瓷滲灌管道看作半圓柱形恒定水頭水源，將管道內水分入滲過程假設為軸對稱，由于供試土壤密實度均勻，進而整個土壤水分入滲過程可以用Richards方程[25]來描述：

式中：θ為土壤體積含水量，cm3/cm3；z和r分別為徑向和垂直空間坐標；K(h)為非飽和土壤導水率，cm/min；t為時間，min。

本試驗模擬的是恒定水頭下微孔陶瓷滲灌管道的土壤水分入滲規(guī)律，土壤體積含水量和非飽和土壤導水率可以用van Genuchten模型[22]的常規(guī)表達式描述：

式中：θs為土壤飽和含水量，cm3/cm3；θr為土壤殘余含水量，cm3/cm3；α為經驗參數，1/cm；h為土壤負壓水勢，cm；Ks為滲透系數，cm/min；Se為土壤有效含水量，cm3/cm3；m和n分別是擬合參數，m= 1 - 1/n，n> 1。土壤水力特性參數見表2。

表2 土壤水力特性Tab.2 Hydraulic properties of soil

1.3.3 系統邊界與初始條件

大多數間作系統中作物根系在豎直方向上主要分布范圍為0～80 cm 深度[23]，為研究灌溉過程中的深層滲漏，因此，土壤水分入滲模擬中使用深度為80 cm 的土壤剖面，模擬區(qū)域為50 cm × 80 cm 的矩形，其中最上層為大氣邊界，寬度為50 cm，最底層為深層滲漏邊界，微孔陶瓷滲灌管道管徑2 cm，厚度1 cm，提供恒定水頭作為恒定水頭邊界，由于左右兩側對稱，本研究只模擬右側土壤剖面土壤水分運移。HYDRUS-2D 模擬中網狀幾何和邊界條件如圖3 所示。本試驗土箱內的供試土樣初始含水量均為0.12 cm3/cm3，運用式（4）將含水量轉換為負壓水頭，公式中參數值如表2所示，計算得到微孔陶瓷滲灌管道恒定水頭邊界處負壓水頭為-8 638.83 cm。

圖3 網狀幾何和邊界條件（單位：cm）Fig.3 Mesh geometry and boundary conditions

1.3.4 模型評價

本模型利用平均絕對值誤差(MAE)、均方根誤差(RMSE)和決定系數R2評估HYDRUS-2D 模型預測水平與實測數據的一致性[24]。

2 結果

2.1 微孔陶瓷滲灌管道出流特性

圖4為微孔陶瓷滲灌管道實測和模擬濕潤鋒、累計入滲量和土壤含水量關系。從圖4 中可以看出，數據點接近于1∶1線，MAE和RMSE分別小于0.8 和0.9，R2大于0.75，這說明所構建的HYDURS-2D模型能以較高精度反映微孔陶瓷滲灌管的濕潤鋒、累計入滲量和土壤含水量。

圖4 實測和模擬濕潤鋒、累計入滲量、土壤含水量間數學關系Fig.4 The mathematical relationships between measured and simulated wetting front, cumulative infiltration and soil water content

圖5（a）為不同工作水頭下微孔陶瓷滲灌管道流量隨時間的變化。由于微孔陶瓷滲灌管道周圍土壤快速接近飽和，滲灌管道周圍土壤含水量對滲灌管道流量具有反饋調節(jié)作用，導致流量下降[26]，從圖5中可以看出，微孔陶瓷滲灌管道流量在1 h 內迅速下降，于48 h 后達到穩(wěn)定出流狀態(tài)，并以相對穩(wěn)定的流量持續(xù)出流（此處稱為穩(wěn)定出流量）。例如，當工作壓力為0.50 m 時，微孔陶瓷滲灌管道流量先在1 h 內從初始的2.164 L/(h·m)降低為初始值的10.76%，然后在1-36 h內流量小幅度下降，最終于48 h后微孔陶瓷滲灌管道保持0.212 L/(h·m)的穩(wěn)定出流量持續(xù)穩(wěn)定出流，這說明微孔陶瓷滲灌管道可以在作物生育期內以極低的流量進行連續(xù)灌溉。此外，微孔陶瓷滲灌管道穩(wěn)定出流量隨工作水頭的增大而增大，當工作水頭從0.25 m 增加至0.75 m 時，微孔陶瓷滲灌管道穩(wěn)定出流量由0.157 L/(h·m)增加至0.270 L/(h·m)，增加了71.97%。

圖5 不同工作水頭下微孔陶瓷滲灌管道流量與時間的關系、不同工作水頭與滲灌管道穩(wěn)定流量的關系Fig.5 The relationship between flow rate and time of microporous ceramic infiltration pipe under different working heads、Relationship between different working heads and stable outflow rate of microporous ceramic infiltration pipe

通過HYDRUS-2D數值模擬的方法，進一步地繪制了如圖5（b）所示的不同工作水頭下微孔陶瓷滲灌管道在穩(wěn)定出流狀態(tài)的流量變化情況。由圖5（b）可知，微孔陶瓷滲灌管道在0.05 ～0.75 m 工作水頭的微壓條件下，由于滲灌管道流量變化是滲灌管道外基質勢和壓力勢作用所致，且滲灌管道在穩(wěn)定出流階段其周圍的土壤水勢為定值，則微孔陶瓷滲灌管道穩(wěn)定出流量與工作水頭呈顯著的線性正相關關系[19,27]。當工作水頭從0.05 m 增加至0.75 m 時，微孔陶瓷穩(wěn)定出流量由0.157 L/(h·m)增加至0.270 L/(h·m)。這說明在實際應用條件下，微孔陶瓷滲灌管道可通過改變工作壓力調節(jié)穩(wěn)定出流量，以適應不同生育期作物的灌溉需水量要求，實現主動灌溉。

2.2 土壤水分的時空變化

由于微孔陶瓷滲灌管道在土箱中入滲形成的濕潤體剖面為近似于左右對稱的圓形，故在此僅展示1/2 濕潤鋒說明運移情況。圖6為不同工作水頭下入滲48 h內微孔陶瓷滲灌管道土壤水分分布圖。從圖6中可以看出，微孔陶瓷滲灌管道附近土壤含水量較高，接近達到土壤飽和含水量，并沿濕潤鋒運移方向土壤含水量逐漸降低。當入滲時間為48 h 時，工作水頭為0.25 m 的微孔陶瓷滲灌管道水平、垂直向上和垂直向下的濕潤鋒運移距離分別為0.253 7 m、0.250 9 m 和0.258 9 m，濕潤面積為0.203 4 m2，并且濕潤面積隨入滲時間的延續(xù)小幅度增大。此外，入滲初期，土壤水分入滲的主要驅動力為土壤吸力，水勢差受工作壓力影響較小[28]，當入滲48 h，工作水頭由0.25 m 增加至0.75 m 時，穩(wěn)定出流狀態(tài)的濕潤面積從0.203 4 m2增加至0.345 6 m2，增幅為69.91%。如圖6 所示，當入滲時間為48 h 時，濕潤體內平均土壤含水量也隨工作水頭的增加而增大，平均土壤含水量從0.25 m的0.243 cm3/cm3增加至0.75 m 的0.406 cm3/cm3，增幅為67.07%。這說明通過調節(jié)工作水頭，可以控制微孔陶瓷滲灌管道土壤水分入滲過程，以實現最優(yōu)水分供給。

圖6 模擬3種工作水頭下48 h后的微孔陶瓷滲灌管道濕潤鋒圖Fig.6 The wetting front of microporous ceramic infiltration pipe after 48 h under three working heads was simulated

為精確描述不同工作水頭下微孔陶瓷滲灌管道土壤水分入滲過程，繪制了如圖7所示的垂直向上、垂直向下和水平濕潤鋒運移距離以及濕潤體內平均土壤含水量隨時間的變化規(guī)律。從圖7（a）～圖7（c）可以看出，在入滲初期，隨著土壤入滲時間增加，垂直向上、垂直向下和水平濕潤鋒運移距離逐漸增加，但隨著時間延長，濕潤鋒在入滲后期運移距離差值逐漸減小，濕潤鋒各方向運移距離變化率與時間呈負相關，濕潤面積變化越來越小，而后基本維持不變[8,29]。此外，垂直向上、垂直向下和水平濕潤鋒運移距離隨工作水頭的增加而增加。以入滲時間48 h 為例，當工作水頭由0.25 m 增加至0.75 m 時，垂直向上濕潤鋒運移距離從25.09 cm 增加至28.6 cm，垂直向下濕潤鋒運移距離從25.89 cm 增加至29.1 cm，水平濕潤鋒運移距離從25.37 cm 增加至28.67 cm。以0.75 m 工作水頭為例，在入滲后期即第48 h 時，垂直向上、垂直向下和水平方向的濕潤鋒分別以0.05 cm/h、0.08 cm/h 和0.04 cm/h 的速率擴散，濕潤面積基本維持穩(wěn)定。由于土壤水分受重力作用，垂直向下濕潤鋒運移距離僅比垂直向上濕潤鋒運移距離和水平濕潤鋒運移距離多0.5 cm和0.43 cm。

圖7 微孔陶瓷滲灌管道濕潤鋒與濕潤體平均土壤含水量隨時間變化規(guī)律Fig.7 The microporous ceramic infiltration pipe wetting front and average soil water content of the wetting body of the changes with time

由圖7（d）可知，濕潤體內平均土壤含水量在1 h內迅速上升，入滲48 h 后，濕潤體內平均土壤含水量保持在穩(wěn)定狀態(tài)。例如，當工作壓力為0.25 m 時，濕潤體內平均土壤含水量先在1 h 內由初始土壤含水量0.120 cm3/cm3迅速增長至0.238 cm3/cm3，然后隨著入滲時間增加，當達到穩(wěn)定出流狀態(tài)，濕潤體內平均土壤含水量保持在0.243 cm3/cm3。此外，濕潤體內平均土壤含水量隨工作水頭增加而增加，入滲48 h后， 0.25 m、0.50 m和0.75 m工作水頭下濕潤體內平均含水量分別為0.243 cm3/cm3、0.348 cm3/cm3和0.406 cm3/cm3。

2.3 微孔陶瓷滲灌管道工作水頭確定

根據HYDRUS-2D 模型模擬結果，表3 給出了入滲48 h 后微孔陶瓷滲灌管道穩(wěn)定出流量、水平、垂直向上、垂直向下濕潤鋒運移距離以及濕潤體內平均土壤含水量與工作水頭間的數學關系。為確定微孔陶瓷滲灌管道最優(yōu)工作水頭，以根水匹配度最高為優(yōu)化指標，建立目標函數：

表3 入滲48 h后土壤水分入滲參數與工作水頭間的數學關系Tab.3 The mathematical relationship between soil water infiltration parameters and working heads after infiltration for 48 h

式中：q為微孔陶瓷滲灌管道穩(wěn)定出流量，L/(h·m)；I為灌水量，mm；S為種植面積，hm2；N為株數，棵；D為灌溉天數；X為水平濕潤鋒運移距離，cm；Z+為垂直向上濕潤鋒運移距離，cm；Z-為垂直向下濕潤鋒運移距離，cm；Lz為根系深度，cm；Lx為根系深度，cm。

約束條件：

通過上述目標函數，結合不同作物需水量及根系空間分布情況，進而根據土壤水分入滲參數，獲得微孔陶瓷滲灌管道最優(yōu)工作水頭。表4給出了在與試驗土樣物理特性相似的試驗條件下常見的冬小麥、玉米和番茄3種根系分布特征密植作物的微孔陶瓷滲灌管道適宜工作水頭。表4中不同密植作物的適宜灌水量、生育期時間、根系深度和寬度可見參考文獻。

表4 典型密植型作物的工作水頭優(yōu)化Tab.4 Working head optimization of typical close planting crops

3 討 論

由于微孔陶瓷滲灌管道依靠其內部相互連通的微米級孔隙出流，其入滲過程可理解為多孔介質出流，因此影響滲灌管道流量變化的主要因素是工作水頭與外部土壤水勢[33]。蔡耀輝等[34]通過分析微孔陶瓷灌水器出流過程中的力能關系，發(fā)現當土壤干燥時灌水器同時受到工作水壓和土壤水吸力的共同作用，并隨土壤含水量的增加，土壤水吸力逐漸降低，最終表現為灌水器出流量隨入滲時間的延續(xù)逐漸降低。本研究中微孔陶瓷滲灌管道出流量也呈現相同的變化趨勢，這說明陶瓷灌水器結構形式并不改變其內部出流機理。然而，Liu 等[9]利用微孔陶瓷灌水器進行溫室番茄灌溉，發(fā)現作物生育期內微孔陶瓷灌水器出流量連續(xù)隨機波動，而本試驗中微孔陶瓷滲灌管道出流量基本維持穩(wěn)定，這可能是土箱試驗條件下土壤水分環(huán)境相對穩(wěn)定所導致的。近年來，Cai 等[35]運用HYDRUS-2D 數值模擬軟件預測了微孔陶瓷灌水器在不同運行管理參數條件下流量變化規(guī)律，并利用灌水器穩(wěn)定出流量對微孔陶瓷灌水器工作水頭進行了優(yōu)化，這說明本研究中以最優(yōu)根水匹配進行滲灌管道工作水頭優(yōu)化是可行的。

王辰元[36]與歐陽淼[37]在地下滴灌土壤水分入滲試驗中發(fā)現了相似結果，然而地下滴灌灌水器所形成的土壤濕潤體近似于橢圓，且垂向濕潤距離顯著大于水平濕潤距離，這主要是工作水頭較高、滴頭流量較大所導致的。微孔陶瓷滲灌與地下滴灌的土壤水分運移規(guī)律和機理不同，在灌溉過程中，微孔陶瓷滲灌在低壓條件下依靠土壤基質勢動態(tài)調節(jié)，濕潤體內土壤含水量變化較小，提供了恒定的水分環(huán)境，地下滴灌雖然可以采取干濕交替的被動恒壓灌溉方式[8]，但土壤含水量不易控制，且操作復雜。低壓微孔陶瓷滲灌管道入滲過程主要受管道內外水勢差影響，入滲后期由于滲灌管道周圍土壤接近飽和[8]，當微孔陶瓷滲灌管道入滲48 h后，其土壤入滲基本維持穩(wěn)定，此時土壤基質勢較低，微孔陶瓷滲灌管道穩(wěn)定出流量受工作壓力影響較為明顯，這一現象與牛文全研究結果吻合[27]，因此可根據作物的需水量，通過調節(jié)工作水頭實現實時補水的主動灌溉。而微壓灌溉條件下的濕潤鋒尚未接觸到深層滲漏邊界，垂直向下運移的水分深度未能達到深層滲漏條件，這說明滲灌管道在控制土壤濕潤區(qū)域和降低土壤水分深層滲漏方面可能存在明顯優(yōu)勢，Liu 等[9]也發(fā)現了相似現象。與傳統微孔陶瓷灌水器相比，由于微孔陶瓷滲灌管道需考慮邊界問題，所形成的濕潤體為平行于土壤表面的圓柱體，故將其作為柱狀面源入滲進行分析[8,38]，因此在實際應用中，為充分利用微孔陶瓷滲灌的入滲特點與提供穩(wěn)定土壤水分環(huán)境的優(yōu)點，可將研究所提出的滲灌管道應用于密植型作物。

進一步地，本研究以最優(yōu)根水匹配為優(yōu)化目標，對冬小麥、玉米和番茄3種典型密植型作物的微孔陶瓷滲灌管道工作水頭進行了優(yōu)化，發(fā)現番茄適宜工作水頭最高，而冬小麥作物工作水頭最低（表4），工作水頭優(yōu)化結果是可以預見的，這是因為與玉米和冬小麥相比，番茄作為一種典型密植型經濟作物，由于根系發(fā)達且吸水能力強，對土壤水分敏感，因此需要更高的灌溉水量以保證產量[32]，而根據過往研究發(fā)現，張子卓[39]認為利用微潤帶灌溉番茄適宜的工作水頭為120～200 cm，蔡耀輝[19]研究發(fā)現微孔陶瓷灌水器無壓或微壓灌溉番茄的工作水頭合理取值則為0～20 cm，因此分別在0～100 cm的工作水頭中設置3組梯度進行研究。然而，本研究僅對黏壤土下微孔陶瓷滲灌管道的土壤水分入滲過程進行了研究，未來尚需對不同土質下的滲灌管道入滲過程及田間應用效果進行研究，以實現該技術的實際推廣應用。

4 結 論

本研究采用注漿成型法制備了一種新型材料性能優(yōu)良的綠色微孔陶瓷滲灌管道，并以此為研究對象，采用試驗研究和數值模擬的方法，利用HYDRUS-2D模擬微孔陶瓷滲灌管道土壤水分入滲過程，為實現最優(yōu)根水匹配，優(yōu)化了微孔陶瓷滲灌管道工作水頭。本研究的主要結論如下：

（1）HYDRUS-2D 模擬軟件預測的濕潤鋒運移距離、累計入滲量和土壤含水量與試驗結果吻合度較高，驗證了模型的適用性與有效性，可較好地預測微孔陶瓷滲灌管道在不同運行管理條件下土壤水分入滲參數變化規(guī)律，為后續(xù)優(yōu)化工作水頭提供了可行性。

（2）微孔陶瓷滲灌管道流量在1 h 內迅速下降，待穩(wěn)定后以穩(wěn)定出流量持續(xù)出流，且滲灌管道穩(wěn)定出流量隨工作水頭的增大而增大，呈顯著的線性正相關關系。此外，滲灌管道附近土壤含水量最高，并沿濕潤鋒運移方向逐漸降低，濕潤面積隨入滲時間的延續(xù)逐漸增大，待入滲穩(wěn)定后變化越來越小，濕潤體內平均土壤含水量和濕潤面積基本維持不變，僅隨工作水頭增加而增大。

（3）本研究對冬小麥、夏玉米和番茄3種典型密植型作物的微孔陶瓷滲灌管道工作水頭進行了優(yōu)化，結合密植作物的根系分布和作物生長需求情況，冬小麥推薦工作水頭為0.423 m，夏玉米推薦工作水頭為0.529 m，番茄推薦工作水頭為0.642 m。