裴青寶 劉偉佳 張建豐 王海偉

(1.西安理工大學(xué)水利水電學(xué)院， 西安 710048； 2.南昌工程學(xué)院江西省水工程安全與資源高效利用工程研究中心， 南昌 330099；3.江西省水利規(guī)劃設(shè)計院， 南昌 330020)

引言

南方紅壤丘陵地區(qū)多種植柑橘、臍橙等經(jīng)濟作物，近年來在丘陵地區(qū)這些作物的灌溉多采用滴灌水肥一體化灌溉方式，本文根據(jù)對柑橘、臍橙等作物根系的調(diào)查和紅壤剖面土壤容重的測定，通過試驗研究不同滴頭流量、間距以及土壤容重下交匯入滲規(guī)律以及土壤水分氮素分布特性，以期為紅壤柑橘、臍橙滴灌設(shè)計參數(shù)選擇提供依據(jù)。

1 材料與方法

1.1 供試土壤

試驗土壤于2016年3月采自江西省尋烏縣澄江鎮(zhèn)汶口村臍橙果園內(nèi)，取地表100 cm的種植土，土壤容重在1.36～1.41 g/cm3之間，將試驗土樣取回，于試驗室風(fēng)干、過2 mm篩、通風(fēng)保存,土壤顆粒分析用馬爾文激光粒度儀測定，粒徑0～0.002 mm、 0.002～0.02 mm、 0.02～2 mm之間的顆粒含量分別為44.25%、38.39%、17.36%，采用國際制土壤質(zhì)地分級標(biāo)準(zhǔn)，最終確定土壤質(zhì)地為壤質(zhì)粘土。試驗在南昌工程學(xué)院灌溉排水實驗中心展開，時間為2016年6—12月，室內(nèi)溫度保持為18～23℃。

1.2 試驗裝置

試驗裝置采用西安理工大學(xué)水資源所設(shè)計生產(chǎn)的土壤水分運動試驗系統(tǒng)，系統(tǒng)由馬氏瓶、有機玻璃土箱、支架等組成；土箱設(shè)計為矩形結(jié)構(gòu)，長×寬×高為40 cm×40 cm×60 cm。土壤含水率的測量采用埋設(shè)在土體內(nèi)不同深度處的TDR探頭測定，TDR探頭2組共13個，分別在1號滴頭水平距離為5 cm,垂直距離為5、10、15、20、25、30 cm處和在濕潤鋒交匯面水平距離為20 cm，垂直距離為0、5、10、15、20、25、30 cm處埋設(shè)TDR探頭，記錄不同時間點的土壤含水率，其他位置通過烘干法測定。Hydrus-3D建模過程中與物理模型相應(yīng)的位置插入觀測點，采用馬氏瓶供水，醫(yī)用針頭模擬滴灌出水，用流量計控制流量。裝置示意圖如圖1所示。

圖1 試驗裝置示意圖Fig.1 Sketch of experimental setup1、8.支架 2、9.出水口 3、10.進氣口 4、11.馬氏瓶 5.1號滴頭 6.2號滴頭 7.土箱

1.3 試驗內(nèi)容

1.4 試驗過程

1.5 數(shù)學(xué)模型

1.5.1水分運動基本方程

滴灌土壤水分運動為三維流動問題。假定土壤各向同性,地表沒有蒸發(fā),土壤初始含水率相同,并不考慮滯后現(xiàn)象,則三維軸對稱點源水分入滲Richard方程可表示為[14,17]

(1)

式中θ——土壤體積含水率，cm3/cm3

h——土壤負壓水頭，cm

r、y、z——坐標(biāo)(z坐標(biāo)向下為正)，cm

t——時間，min

Kh——土壤非飽和導(dǎo)水率，cm/min

Hydrus-3D模擬時需要的非飽和土壤水分特征曲線θh、土壤非飽和導(dǎo)水率Kh采用van Genuchten模型表示[14,18]，并不考慮滯后效應(yīng)

(2)

θh=θs(h≥0)

(3)

(4)

其中

(5)

式中θs、θr——土壤飽和含水率和殘余含水率，cm3/cm3

Ks——土壤飽和導(dǎo)水率，cm/min

l——孔隙連通性參數(shù)，對大多數(shù)土壤來說可取0.5

a、n、m——擬合經(jīng)驗參數(shù)

Se——有效含水率(飽和度)

(6)

其中

(7)

(8)

(9)

qr——縱向上的土壤水分通量

qz——橫向上的土壤水分通量

Drr、Dzz、Drz——水動力彌散系數(shù)張量的分量

q——土壤水通量的絕對值

DL、DT——溶質(zhì)的縱向和橫向彌散度，L

Dw——自由水中的分子擴散系數(shù)

τ——溶質(zhì)的彎曲系數(shù)，通常表示為土壤體積含水率的函數(shù)

1.6 模型區(qū)域以及邊界條件和初始條件

模型求解區(qū)域以實際土箱尺寸為基礎(chǔ)建立幾何模型，建模區(qū)域如圖2所示，Hydrus-3D 不能描述移動的水分邊界,但可以模擬邊界條件隨時間變化的過程，因此在試驗過程中測定滴頭下方積水的區(qū)域及水頭高度。試驗開始后15 min時達到飽和區(qū)域的58%、積水深度為0.18 cm(用游標(biāo)卡尺通過測量滴頭下方積水深度得到);灌水后25 min達到飽和區(qū)域的80%、積水深度為0.50 cm;38 min 后，飽和區(qū)穩(wěn)定, 積水深度為0.85 cm。根據(jù)試驗測試結(jié)果，飽和區(qū)域在很短的時間內(nèi)達到穩(wěn)定，所以在模擬過程中不考慮其范圍的變化，只輸入最終的范圍值，設(shè)定飽和區(qū)半徑為定值Rs(滴頭地表最小濕潤半徑)如圖2所示。試驗過程中無表面蒸發(fā)，則水分運動的上邊界條件可以表示為

h=0.85 cm (0≤r≤R,z=Z,t≥0)

(10)

(11)

式中R、Z——模擬區(qū)域邊界(裝置物理邊界)在徑向和垂直方向的坐標(biāo)

圖2 模型求解區(qū)域示意圖Fig.2 Schematic diagram of solving model

本次模擬輸入邊界條件中水頭隨時間的變化個數(shù)為4，按照變化的時間，在各段時間上輸入相應(yīng)的水頭值。由于滴頭周圍有積水產(chǎn)生,因此溶質(zhì)運移的上邊界條件采用一類邊界條件

C(r,z)=Ca(0≤r≤Rs,z=Z,t>0)

(12)

側(cè)面為不透水邊界

(13)

(14)

下邊界為自由排水邊界

(15)

(16)

θ(r,z)=θ0(0≤r≤R,0≤z≤Z,t=0)

(17)

C(r,z)=C0(0≤r≤R,0≤z≤Z,t=0)

(18)

式中θ0——土壤初始含水率，cm3/cm3

2 模型參數(shù)

2.1 土壤水力學(xué)參數(shù)

依據(jù)1.5.1節(jié)所述土壤水分特性的 VG 模型參數(shù)，結(jié)合不同容重的水分特征曲線以及飽和導(dǎo)水率試驗結(jié)果，并根據(jù)顆粒分析的成果輸入Hydrus-3D里計算結(jié)果并進行對比分析，最終確定合適的參數(shù)如表1所示。

表1 不同紅壤容重van Genuchten模型參數(shù)Tab.1 Van Genuchten model parameters of red soil with different bulk densities

圖3 滴灌交匯入滲濕潤鋒運移變化Fig.3 Changes of drip irrigation, intersection infiltration and wetting front transport

2.2 溶質(zhì)運移參數(shù)

3 結(jié)果與分析

3.1 容重變化影響下濕潤鋒模擬值與試驗值的對比分析

試驗開始后在滴頭下方形成一個積水區(qū)域，隨著時間的變化，積水區(qū)域趨于穩(wěn)定。水氮開始在土壤中運動，將不同時刻濕潤鋒的變化在有機玻璃板上用記號筆繪出。圖3為滴灌交匯入滲試驗濕潤鋒變化過程實測及模擬照片。

圖4 不同容重條件下濕潤鋒變化模擬值與實測值Fig.4 Simulated and measured values of wetting front under different bulk densities

容重變化影響下濕潤鋒模擬值與試驗值如圖4所示(圖中γ表示容重)。由圖4可見，通過試驗測得的濕潤鋒過程和Hydrus-3D模擬的值在觀測時間點上有同步性，其相對偏差在9%以內(nèi)。5個不同容重的濕潤鋒模擬值和實測值經(jīng)過顯著性檢驗后差異均不顯著，這說明數(shù)值模擬可以反映不同容重條件下濕潤鋒的運移規(guī)律，這也與李久生等[14]、張林等[18]的研究結(jié)果相同。模擬的濕潤鋒以滴頭為中心擴散，濕潤半徑更為圓滑均勻，從圖3b的模擬圖也可以看出這個現(xiàn)象，而實測的濕潤鋒某一時刻濕潤鋒面值則出現(xiàn)凹凸現(xiàn)象；分析認為Hydrus在設(shè)定不同的容重后軟件認為土壤內(nèi)部是均一的各向同性，且不受邊壁滯后效應(yīng)的影響。所以Hydrus-3D模擬的濕潤鋒面推進較快且均勻，而實測值由于土壤為擾動土，裝土的質(zhì)量和土層的結(jié)合等情況影響著濕潤鋒的運移；所以實際工作中如果需要采用模型來預(yù)測濕潤鋒的運移時需要考慮到滯后效應(yīng)。容重從1.2 g/cm3到1.4 g/cm3，雙滴頭滴灌下出現(xiàn)了濕潤鋒的交匯現(xiàn)象，而交匯的時間則與容重有關(guān)，容重為1.2 g/cm3的濕潤鋒在灌水后273 min時發(fā)生交匯，而容重越大交匯過程越緩慢，容重1.4 g/cm3在461 min時才發(fā)生交匯。試驗土為紅壤，容重較小時土壤中的孔隙率大，導(dǎo)水速率快，所以濕潤鋒推進得比較快。當(dāng)容重增加時，孔隙率減少，加上紅壤遇水后顆粒變得粘稠，堵塞了土壤中的孔隙，造成濕潤鋒不能推進。試驗取土的江西省贛州市尋烏縣臍橙果園內(nèi)，現(xiàn)場調(diào)查其距地表100 cm以內(nèi)的土壤容重在1.36～1.41 g/cm3之間，根據(jù)模擬和實測結(jié)果可知，高容重情況下土壤濕潤鋒難以推進，交匯發(fā)生較緩慢，影響到臍橙的灌溉。所以在灌溉之前需要對土壤結(jié)構(gòu)進行前期處理，疏松土層，便于水肥能夠入滲到根系周圍。

3.2 濕潤體內(nèi)土壤含水率的分布模擬

結(jié)合圖5、6土壤含水率模擬和實測值，兩者的平均相對誤差為9.5%以內(nèi)。其中，90%的時間節(jié)點上的模擬值和實測值的差值在3.5%之內(nèi)，而且各觀測點誤差情況類似。差異較大主要在初始時刻和試驗結(jié)束時段，初始觀測點上實測的含水率值均低于模擬值；實測值通過埋設(shè)在觀測點上的TDR探頭測定，探頭在使用之前設(shè)置在不同的土壤容重中，在灌水后通過與烘干法進行含水率測定率定，并對探頭進行修正。所以在試驗過程中探頭的值可以反映出實際情況，而模擬值出現(xiàn)的偏差主要是在建模過程中模型為理想狀態(tài)，輸入顆分數(shù)據(jù)和土壤吸力后神經(jīng)網(wǎng)絡(luò)預(yù)測出一個含水率值，從而造成與實際土壤初始含水率之間的誤差。而在灌水結(jié)束時出現(xiàn)的部分實測值大于模擬值的原因是，模型在該時刻模擬結(jié)束，含水率定格在這個狀態(tài)，不會變化，模型忽視了土壤水分再分布這一過程，所以實測值大于模擬值。除此之外，模型模擬的精度較高，可以反映出滴灌交匯入滲過程中濕潤體內(nèi)土壤含水率的分布情況，這與李久生等[14-15]的研究結(jié)果相同。

圖5 不同容重條件下距離滴頭5 cm處含水率變化Fig.5 Changes of water content under condition of different densities with 5 cm away from emitter

圖6 不同容重條件下交匯處含水率模擬值與實測值Fig.6 Confluence of moisture content of simulated and measured values under condition of different densities

圖5中，總體而言，隨著入滲時間變化，觀測點的土壤含水率呈現(xiàn)增加的現(xiàn)象，土壤深度越深，含水率增加值越小。容重1.2 g/cm3垂直距離5 cm的觀測點土壤含水率在入滲開始后就快速增加，200 min后接近飽和，并維持穩(wěn)定狀態(tài)；其他4個不同容重該位置處的觀測點模擬和實測值均顯示該現(xiàn)象。其余觀測點隨著深度的增加，距離滴頭間距越大含水率越小，并受到容重的影響，同一位置處的觀測點，容重越大，含水率相應(yīng)地減小，這表明，高容重壤土阻礙了水分的入滲。

圖6中交匯處的土壤含水率，在發(fā)生交匯后迅速增加，并與深度和容重有關(guān)，深度較深容重較大，含水率增加的值較小。交匯面處含水率低于同一位置處滴頭下方的觀測點的值，由此可知在濕潤體內(nèi)同一深度處，含水率的分布是不均勻的，滴頭下方處土壤含水率高于交匯處的，容重越大這一現(xiàn)象愈加明顯。對于臍橙等根系發(fā)達的經(jīng)濟作物，在滴頭間距相同時，兩滴頭交匯處的土壤含水率較低，會影響這一區(qū)域內(nèi)根系吸水，要加強交匯區(qū)域的灌水量。

圖7 不同容重條件下距離滴頭5 cm處質(zhì)量濃度模擬值與實測值Fig.7 N-N mass concentration simulated and measured values with 5 cm away from emitter under different bulk densities

圖8 不同容重條件下交匯處質(zhì)量濃度模擬值與實測值Fig.8 N-N simulated and measured values under different bulk densities intersection

4 結(jié)論

(1)Hydrus-3D可以較好地模擬出濕潤鋒的推移過程，模擬值與實測值的精度在9%以內(nèi)，軟件的模擬值更加理想化，在運用到實際中模擬濕潤體范圍時需要考慮紅壤的空間差異性。

(2)含水率的模擬值和實測值具有很高的一致性，模型可以用于交匯情況下土壤濕潤體內(nèi)含水率分布的模擬。容重對含水率的分布產(chǎn)生較大的影響，交匯面處的含水率低于同一深度滴頭下方的值，模擬和實測結(jié)果均表明，高紅壤容重阻礙濕潤鋒的推進和含水率的增加。

1 CHARLES M B. Rapid field evaluation of drip and microspray distribution uniformity[J]. Irrigation and Drainage Systems, 2004, 18(4): 275-297.

2 ZHOU Jianbin，XI Jingen，LI Shengxiu. Leaching and transformation of nitrogen fertilizers in soil after application of N with irrigation: a soil column method[J].Pedosphere,2006,16(2):245-252.

3 YANG Lijuan, ZHANG Yulong, XU Hui, et al. Effect of irrigation methods on soil water loss and water use efficiency of tomato in greenhouse[J]. Journal of Irrigation and Drainage, 2004, 23(3): 49-51.

4 李就好，譚穎，張志斌，等.滴灌條件下磚紅壤水分運動試驗研究[J]. 農(nóng)業(yè)工程學(xué)報，2005，21(6)：36-39.

LI Jiuhao, TAN Ying, ZHANG Zhibin, et al. Experimental study on water movement of latosol under drip irrigation[J].Transactions of the CSAE, 2005, 21(6): 36-39.(in Chinese)

5 李憑峰,譚煌,王嘉航,等.滴灌水肥條件對櫻桃產(chǎn)量、品質(zhì)和土壤理化性質(zhì)的影響[J/OL]．農(nóng)業(yè)機械學(xué)報，2017,48(7):236-246．http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20170730&flag=1&journal_id=jcsam.DOI:10.6041/j.issn.1000 1298.2017.07.030.

LI Pingfeng, TAN Huang, WANG Jiahang, et al. Effect of water and fertilizer conditions under drip irrigation on yield, quality of Cherry and physic-chemical properties of soil [J/OL].Transactions of the Chinese Society for Agricultural Machinery，2017,48(7): 236-246． (in Chinese)

6 孫海燕, 李明思, 王振華,等.滴灌點源入滲濕潤鋒影響因子的研究[J]. 灌溉排水學(xué)報，2004,23(3):14-18.

SUN Haiyan，LI Mingsi，WANG Zhenhua, et al. Influence factors on soil wetting front under point drip irrigation [J].Journal of Irrigation and Drainage,2004,23(3):14-18. (in Chinese)

7 張振華,蔡煥杰,郭永昌,等.滴灌土壤濕潤體影響因素的實驗研究[J].農(nóng)業(yè)工程學(xué)報,2002,18(2):17-20.

ZHANG Zhenhua, CAI Huanjie, GUO Yongchang, et al. Experimental study on factors effecting soil wetted volume of clay loam under drip irrigation[J]. Transactions of the CSAE,2002, 18(2): 17-20. (in Chinese)

8 關(guān)紅杰,李久生, 栗巖峰．干旱區(qū)滴灌均勻系數(shù)對土壤水氮分布影響模擬[J/OL]．農(nóng)業(yè)機械學(xué)報，2014,45(3):107-113．http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1&file_no=20140319&journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.03.019.

GUAN Hongjie, LI Jiusheng, LI Yanfeng. Modeling distributions of water and nitrate in soil as affected by drip system uniformity under arid conditions[J/OL]．Transactions of the Chinese Society for Agricultural Machinery，2014,45(3):107-113． (in Chinese)

9 BRANT A, BRESLER E, DINER N, et al.Infiltration from a trickle source：I. mathematical models[J].Soil Sciences Society of America Proceedings, 1971, 35:675-682.

10 SINGH D K, RAIPUT T B S, SINGH D K, et al. Simulation of soil wetting pattern with subsurface drip irrigation from line source[J]. Agricultural Water Management, 2006, 83(1): 130-134.

11 許迪，程先軍．地下滴灌土壤水分運動和溶質(zhì)運移數(shù)學(xué)模型的應(yīng)用[J]．農(nóng)業(yè)工程學(xué)報，2002，18(1)：27-31．

XU Di, CHENG Xianjun. Mathematical model for simulating water flow and solute transport during non-steady diffusion from subsurface trickle source[J].Transactions of the CSAE,2002，18(1)：27-31． (in Chinese)

12 李光永,曾德超,鄭耀泉. 地表點源滴灌土壤水運動的動力學(xué)模型與數(shù)值模擬[ J]. 水利學(xué)報, 1998, 29(11): 21-25.

LI Guangyong, ZENG Dechao, ZHENG Yaoquan. The mathematical model of water movement from a surface trickle source[J]. Journal of Hydraulic Engineering, 1998,29(11): 21-25. (in Chinese)

13 COTEC M,BRISTOW K L, CHARLES W P B, et al. Analysis of soil wetting and solute transport in subsurface trickle irrigation[J]. Irrigation Science, 2003, 22(3-4): 143-156.

14 李久生，張建君，饒敏杰．滴灌施肥灌溉的水氮運移數(shù)學(xué)模擬及試驗驗證[J]．水利學(xué)報，2005，36(8)：932-938．

LI Jiusheng, ZHANG Jianjun, RAO Minjie. Model verification of water and nitrate transport from a surface point source[J].Journal of Hydraulic Engineering, 2005, 36(8): 932-938. (in Chinese)

15 陳若男，王全九，楊艷芬. 新疆礫石地葡萄滴灌帶合理設(shè)計及布設(shè)參數(shù)的數(shù)值分析[J]. 農(nóng)業(yè)工程學(xué)報，2010，26(12)：40-46.

16 黃凱，蔡德所，潘偉，等.廣西赤紅壤甘蔗田間滴灌帶合理布設(shè)參數(shù)確定[J].農(nóng)業(yè)工程學(xué)報，2015，31(11)：136-139.

HUANG Kai, CAI Desuo,PAN Wei, et al. Determination of drip tapes layout parameters for irrigation of sugarcane in latosolic red Autonomous Region[J]. Transactions of the CSAE,2015,31(11)：136-139. (in Chinese)

17 RICHARDS L A. Capillary conduction of liquids through porous mediums[J]. Physics, 1931, 1(5): 318-333.

18 張林，吳普特，范興科. 多點源滴灌條件下土壤水分運動的數(shù)值模擬[J].農(nóng)業(yè)工程學(xué)報，2010,26(9):40-46.

ZHANG Lin, WU Pute, FAN Xingke. Numerical simulation of soil water movement with drip irrigation of multiple point source[J]. Transactions of the CSAE, 2010, 26(9): 40-46. (in Chinese)

19 SIMUNEK J, SE J M, VANGENUCHTEN M T.HYDRUS-2D simulating water flow,heat,and solute transport in two-dimensional variably saturated media[M].Calif: International Ground Water Modeling Center,Riverside,1999.

20 李久生,張群君,薛克宗.滴灌施肥灌溉原理與應(yīng)用[M].北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社,2003.

21 張效先.飽和條件下田間土壤縱向及橫向彌散系數(shù)的試驗和計算[J]．水利學(xué)報，1989,20(1):1-9.

ZHANG Xiaoxian.A method for determining the longitudinal dispersion parameters of field soils and lateral under saturated conditfion[J].Journal of Hydraulic Engineering,1989,20(1):1-9. (in Chinese)

22 馬軍花.農(nóng)田土壤硝態(tài)氮運移動態(tài)的數(shù)學(xué)模擬[D].北京:中國農(nóng)業(yè)大學(xué), 2001.