何 靖, 白 丹, 郭 霖, 王新端

(西安理工大學 水利水電學院, 西安 710048)

?

豎管地下灌溉粉質(zhì)壤土入滲濕潤體的試驗研究

何 靖, 白 丹, 郭 霖, 王新端

(西安理工大學 水利水電學院, 西安 710048)

豎管灌水器是一種在低壓(0.6～2.0 m)狀況下對作物進行根部地下灌溉的新型灌水技術的核心部件。為了研究該灌水器在不同影響因素不同組合情況下濕潤體特征參數(shù)值的變化，為豎管灌水器的進一步研究提供借鑒。試驗采用正交試驗設計安排了9組試驗對壓力水頭、土壤初始含水率、豎管灌水器的豎管直徑和土壤容重(每個因素取3個水平)4個因素進行了室內(nèi)試驗研究。入滲試驗7 h之后獲得的結(jié)果表明：濕潤體的濕潤半徑大小排序為Y負≥X≥Y正，濕潤體的含水率分布在水平方向呈現(xiàn)為圓形擴散，而豎直方向呈現(xiàn)為橢圓形擴散，表明重力對水分入滲有一定的影響。豎管地下灌溉濕潤體平均含水率變化范圍為7.5%～33.3%，能夠更好地滿足不同作物的需水要求。

豎管灌水器; 濕潤體; 地下灌溉; 含水率; 正交試驗

豎管灌水器是豎管地下灌溉技術的核心部件，在灌水系統(tǒng)工作時，豎管灌水器埋于地下，下端敞口處與土壤接觸形成一個水土接觸面[1]，水分經(jīng)此接觸面進入土壤，作用到作物的根系活動土層，為作物提供生長所需的水分；該灌水器技術的灌水壓力為0.6～2 m，有效的減少了地下灌溉的能源消耗；豎管灌水器的豎管有一定的長度且其直徑較大(4，8，12 mm)，在灌水結(jié)束時可以允許少量的土壤淤積，能夠很好防止灌水器的堵塞現(xiàn)象的發(fā)生。豎管灌水器入滲研究來源于國家自然科學基金資助項目《豎管地下灌溉毛管水力特性及土壤水分布對土壤空間變異的動態(tài)響應》。地下灌溉濕潤體的特征值的變化規(guī)律和定量化解法，是特定的灌水技術進行灌溉系統(tǒng)設計和對作物需水量的供給進行管理的基礎[2]。近年來，國內(nèi)外大量的學者對滴灌[3-8]、膜孔灌[9-11]、無壓灌溉[12-13]、負壓灌溉[14-15]、微潤灌溉[16-17]條件下，土壤濕潤體的特征進行了研究。趙穎娜等[18]研究了不同流量和不同歷時入滲情況下，土壤濕潤體的特征變化，研究結(jié)果表明入滲時間與濕潤體的擴散距離為顯著的冪函數(shù)關系，濕潤體的體積和含水率的分布受灌水器流量控制。王超等[19]研究了地下滴灌條件下，不同滴頭流量和灌水量對濕潤體特性的影響，流量和灌水量都選取了3個水平進行了研究，結(jié)果表明濕潤峰的運移與流量的時間之間為正相關關系。王平等[20]研究了不同流量和濕潤比條件下濕潤體特性的變化，采用2種流量和3種濕潤比進行了研究，研究結(jié)果表明濕潤體的體積變化與時間成線性變化關系，大流量條件下的濕潤體的濕潤比大，灌水均勻度高?，F(xiàn)有的對單點源入滲的濕潤體特性研究大多集中在流量、灌水歷時、灌水量等因素方面，豎管灌水器地下灌溉作為一種新型的地下灌溉技術，本試驗在壓力水頭、豎管灌水器的豎管直徑、土壤初始含水率和土壤容重變化的情況下，初步研究該灌水器的濕潤體特性。

1 試驗方案與裝置

1.1 試驗方案

參照地下滴灌的入滲試驗，本試驗選定壓力水頭、豎管灌水器的豎管直徑、土壤初始含水率和土壤容重4個因素為初步試驗因素，每個因素選取3個水平進行室內(nèi)試驗。試驗中，水土接觸面埋深均為20 cm。在本次試驗中，每組入滲試驗的入滲時間均為7 h，試驗安排見表1。

1.2 試驗土壤

本次試驗采用的試驗用土取自西安周邊，用MS2000型激光分析粒度儀測定其顆粒級配，為粉質(zhì)壤土。該土壤顆粒粒徑小于0.001的含量為3.817%，小于0.005的含量為19.857%，小于0.01的含量為34.065%，小于0.1的含量占到總量的99.868%，小于0.2的含量為99.997%。

1.3 試驗裝置

試驗裝置主要由馬氏瓶(7×6×60 cm長×寬×高)、試驗土箱(四周為10 mm鋼化玻璃，底部為5 mm鋼板，尺寸：長41 cm×寬31 cm×高55 cm)、鐵架及核心部件豎管灌水器(豎管直徑分別為4，8，12 mm)。

表1 試驗方案

1.4 數(shù)據(jù)記錄

土壤濕潤體含水率分布測量：以水土接觸面為參考平面，建立直角坐標系，豎直方向為Y軸，水平方向為X軸。試驗結(jié)束后，將濕潤體取出，以Y軸為對稱軸從上向下將濕潤體剖開，用直尺測量豎直和水平方向的濕潤半徑，其中Y軸測量水土接觸面上下兩邊的濕潤半徑；由于水平方向濕潤面為圓，以直徑為軸兩邊對稱，故只測一邊濕潤半徑。在坐標軸上每隔2 cm取一個點(不足2 cm的按照所余距離取最后一點)，取一定量的土壤用烘干法測量其土壤含水率。

2 結(jié)果與分析

2.1 影響因素顯著性分析

2.1.1 直觀分析 數(shù)據(jù)的直觀分析，應用Excel 2010軟件進行，獲得的直觀分析結(jié)果見表2。

從表2中可以看出，壓力水頭的極差值最大，表明壓力水頭對豎管灌水器累計入滲量的影響程度最大；土壤容重的極差值最小，表明其對累計入滲量的影響程度最??；豎管直徑和土壤初始含水率的極差值介于上述二者之間，表明對累計入滲量有一定的影響程度，且豎管直徑的影響程度大于土壤初始含水率的影響程度。

2.1.2 方差分析 根據(jù)正交試驗的方差分析，將極差值最小的因素作為誤差項進行處理，獲得方差分析結(jié)果(置信度為95%)見表3。

由表3中的結(jié)果可以看出：在置信度為95%的情況下，只有壓力水頭呈顯著性，其余兩個因素呈非顯著性，豎管灌水器的豎管直徑的影響作用大于土壤初始含水率(F比值大小不同)，表明壓力水頭對豎管灌水器入滲影響起到?jīng)Q定性作用。

表2 直觀分析值

注：表2中，Ⅰj，Ⅱj，Ⅲj表示表示各因素在每一水平下的平均累計入滲量。

表3 方差分析值

注：F0.05(2,2)=19。

2.2 濕潤體含水率分析

根據(jù)表1安排的試驗方案，進行室內(nèi)試驗，獲得因素不同組合情況下濕潤體X軸和Y軸的濕潤半徑，見表4。

表4 濕潤體半徑 cm

試驗結(jié)束后，測得不同試驗組中濕潤體不同位置的含水率。表5中給出了部分試驗方案的濕潤體含水率。

由表4和表5中的數(shù)據(jù)可以看出：距水土接觸面的濕潤半徑：X軸濕潤半徑大于Y正濕潤半徑、小于Y負濕潤半徑；土壤濕潤體含水率的變化：隨Y軸和X軸半徑土水結(jié)合面濕潤半徑越遠，濕潤體土壤含水率越??；水分受重力影響，Y負方向的入滲驅(qū)動力略大于正向，使?jié)駶欝wY負方向略凸出，濕潤半徑增大。

表5 濕潤體不同點含水率

注：表中距離為距水土接觸面的距離,水土接觸面處探頭計入X軸。

由表5可以看出，水土接觸面處的含水率最大，并在水土接觸面周圍形成一個土壤含水率趨于飽和[21-22]的飽和帶，在該過程中灌溉水受基質(zhì)勢和壓力勢的共同作用，基質(zhì)勢的作用逐漸減小，壓力勢的作用逐漸增大；當土壤達到飽和含水率時，基質(zhì)勢為零，此時灌溉水入滲主要受壓力勢的作用。隨著入滲時間的不斷增大，土壤含水率飽和圈逐漸增大，水分向周圍土壤的擴散速度逐漸減小，最終趨于一個穩(wěn)定值。此時，濕潤體與試驗用土壤產(chǎn)生的接觸面處，土壤含水率產(chǎn)生跳躍式減小，形成一層干濕交接面。

濕潤體含水率的減小在水平方向即X軸上，呈對成型減小，即水平向上的含水率變化呈現(xiàn)圓環(huán)式減?。回Q直方向即Y軸上，其正向含水率減小略小于負向減小，即相同距離時，正向點上的含水率略大于負向點上的含水率，而負向的濕潤半徑略大于正向的濕潤半徑。

微潤灌溉現(xiàn)已廣泛的應用于實際生產(chǎn)中[23]，濕潤體的平均含水率變化范圍為5%～25.1%[15]；豎管地下灌溉濕潤體的平均含水率分別為33.3%，13.4%，10.5%，7.5%，12.0%，18.6%，26.0%，21.0%和20.0%，可變范圍較之微潤灌溉略大，為7.5%～33.3%，表明豎管地下灌溉入滲累計入滲量的變化范圍大，是可行的。在實際生產(chǎn)中可以根據(jù)作物的不同需水要求進行調(diào)控，盡可能的滿足作物對水分的不同需求。

3 結(jié) 論

(1) 壓力水頭對豎管地下灌溉入滲有顯著性影響，豎管管徑次之，其次為土壤初始含水率、土壤容重；

(2) 豎管地下灌溉粉質(zhì)壤土濕潤體含水率的分布在水平方向呈現(xiàn)為圓環(huán)式濕潤鋒，豎直方向上：正向減小略慢于負向，呈現(xiàn)出類橢圓形濕潤鋒，濕潤半徑大小排序：Y負≥X≥Y正；

(3) 豎管地下灌溉濕潤體的平均含水率變化范圍大于微潤灌溉，為7.5%～33.3%，表明豎管地下灌溉入滲累計入滲量的變化范圍大，能夠?qū)崟r調(diào)控滿足作物對水分的不同需求。

[1] Ben-Gal A, Lazorovitch N, Shani U. Subsurface drip irrigation in gravel-filled cavities [J]. Vadose Zone Journal, 2004,3(4):1407-1413.

[2] 趙偉霞,蔡煥杰,陳新明,等.基于土壤濕潤體特征值的無壓灌溉灌水定額模型研究[J].節(jié)水灌溉,2008(12):5-8.

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

[4] 李明思,康紹忠,孫海燕.點源滴灌滴頭流量與濕潤體關系研究[J].農(nóng)業(yè)工程學報,2006.22(4):32-35.

[5] Lubana P P S, Narda N K. SW—Soil and water: modelling soil water dynamics under rrickle emitters：A review[J]. Journal of Agricultural Engineering Research, 2001,78(3):217-232.

[6] 趙偉霞,張振華,蔡煥杰,等.間接地下滴灌土壤濕潤體特征參數(shù)[J].農(nóng)業(yè)工程學報,2010,26(4):87-92.

[7] Hinnell A C, Lazarovitch N, Furman A, et al. Neuro-Drip:estimation of subsurface wetting patterns for drip irrigation using neural networks[J]. Irrigation Science, 2010,28(6):535-544.

[8] Yao W W, Ma X Y, Li J, et al. Simulation of point source wetting pattern of subsurface drip irrigation[J]. Irrigation Science, 2011,29(4):331-339.

[9] 費良軍,李發(fā)文,吳軍虎.膜孔灌單向交匯入滲濕潤體特性影響因素研究[J].水利學報,2003,5(5):62-68.

[10] 費良軍,董玉云,朱興華.膜孔單點源肥液入滲濕潤體特性試驗研究[J].農(nóng)業(yè)工程學報,2006,22(12):78-81.

[11] 范嚴偉,馬孝義,王波雷,等.膜孔灌土壤濕潤體水分分布與入滲特性數(shù)值模擬[J].農(nóng)業(yè)機械學報,2008,39(11):35-41.

[12] 趙偉霞,蔡煥杰,陳新明,等.無壓灌溉土壤濕潤體含水率分布規(guī)律與模擬模型研究[J].農(nóng)業(yè)工程學報,2007,23(3):7-12.

[13] 趙偉霞,蔡煥杰,陳新明,等.無壓灌溉土壤濕潤體特征值變化規(guī)律與經(jīng)驗解[J].灌溉排水學報,2009,28(2):30-33.

[14] 肖娟,江培福,郭秀峰,等.負水頭條件下水質(zhì)對濕潤體運移及水鹽分布的影響[J].農(nóng)業(yè)機械學報,2013,44(5):101-107.

[15] 冀榮華,王婷婷,祁力鈞,等.基于HYDRUS-2 D的負壓灌溉土壤水分入滲數(shù)值模擬[J].農(nóng)業(yè)機械學報,2015,46(4):113-119.

[16] 牛文全,張俊,張琳琳,等.埋深與壓力對微潤灌濕潤體水分運移的影響[J].農(nóng)業(yè)機械學報,2013,44(12):128-134.

[17] 薛萬來,牛文全,羅春艷,等.微潤灌溉土壤濕潤體運移模型研究[J].水土保持學報,2014,28(4):49-54.

[18] 趙穎娜,汪有科,馬理輝,等.不同流量對滴灌土壤濕潤體特征的影響[J].干旱地區(qū)農(nóng)業(yè)研究,2010,28(4):30-35.

[19] 王超,李援農(nóng).地下滴灌條件下濕潤體特性的試驗研究[J].中國農(nóng)村水利水電,2011(3):38-40.

[20] 王平,胡笑濤,喬林然,等.流量與濕潤比對濕潤體特性影響研究[J].中國農(nóng)村水利水電,2015(11):11-15.

[21] 仵峰,吳普特,范永申,等.地下灌溉條件下土壤水能態(tài)研究[J].農(nóng)業(yè)工程學報,2008,24(12):31-35.

[22] 仵峰,李王成,李金山,等.地下滴灌灌水器水力性能試驗研究[J].農(nóng)業(yè)工程學報,2003,19(2):85-88.

[23] 何玉琴,成自勇,張芮,等.不同微潤灌溉處理對玉米生長和產(chǎn)量的影響[J].華南農(nóng)業(yè)大學學報,2012,33(4):566-569.

Study on Wetted Volume of Silt Loam Under Sub-irrigation with Vertical Tube Emitter

HE Jing, BAI Dan, GUO Lin, WANG Xinduan

(InstituteofWaterResourcesandHydroelectricEngineering,Xi′anUniversityofTechnology,Xi′an710048,China)

Vertical tube emitter as a new sub-irrigation emitter which works under a low pressure head (0.6～2.0 m) for the root system growth zone, studying its infiltration characteristic is a necessary. The orthogonal experimental design was adopted to arrange the influence factors including pressure head， initial soil moisture content， diameter of vertical tube and soil bulk density to Table L9(34) of orthogonal experimental design (all the factors have 3 levels) and to study the characteristic parameters of wetted soil volume under sub-irrigation with vertical tube emitter, all the test followed the same study time (7 h). After the test, the radius of wetted soil volume had been measured by getting it out of the soil box and it was a whole because of the matrix potential. The radius of soil wetted volume was measured by manual inspection, the soil wetted volume was divided into 2 sections along axisY. Then, the interval was marked as 2 cm for a test point to collect little wetted soil, oven drying method was used to measure the soil moisture content. The results pointed out that the order of radius wasYnegative axis≥Xaxis≥Ypositive axis, and moisture content had an irregular distribution as like in horizontal plane as a circle, but in vertical plane as a ellipse, as well as that the variation of the average moisture content in soil wetted volume of vertical tube sub-irrigation greater than micro-moist irrigation can have a good satisfaction for water demand of different crops. This conclusion may have important directive significance for the further study on infiltration performance of a vertical tube emitter in subsurface irrigation and the application of a vertical tube emitter.

vertical tube emitter; wetted soil volume; sub-irrigation; soil moisture content; orthogonal experimental design

2015-04-15

2016-04-29

國家自然科學基金(41571222;51279156)

何靖(1992—),男,陜西咸陽人,碩士研究生,主要從事節(jié)水灌溉理論與技術研究。E-mail:694081783@qq.com

白丹(1960—),男,重慶人,博士,教授,博士生導師,主要從事節(jié)水灌溉理論與技術研究。E-mail:baidan@xaut.edu.cn

S275.4

A

1005-3409(2016)06-0069-04