徐家屯，馬海姣，陳 慧，朱 艷，王云霏，蔡煥杰(西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室，陜西 楊凌 712100)

畦灌是大田作物廣泛采用的一種灌水方法，其灌水質(zhì)量評價及合理的灌水技術(shù)參數(shù)優(yōu)化組合一直是人們長期研究的重要課題[1,2]。而在冬小麥-夏玉米輪作周期內(nèi)，由于土地翻耕、灌溉、降雨、農(nóng)作物生長等因素作用，使得土壤入滲特性及綜合糙率系數(shù)存在著明顯的時空變異性，農(nóng)業(yè)灌溉措施在這種變異條件下，其灌溉效果也必然發(fā)生較大的變化，采用相同的灌水技術(shù)參數(shù)則顯著地降低了畦灌灌溉質(zhì)量[3-5]。在實際應(yīng)用中，對地面灌溉灌水質(zhì)量評價及灌水方案設(shè)計時，應(yīng)充分考慮這種時空變異性?；诖?，國內(nèi)外學(xué)者對其進行了大量的研究[6-9]，取得了一定的成果，但大多數(shù)學(xué)者集中在土壤入滲特性及糙率系數(shù)的空間變異性方面的研究，較少研究夏玉米-冬小麥輪作周期內(nèi)土壤的入滲特性及糙率系數(shù)的時間變異性對灌水過程及灌溉性能的影響。因此，本文以涇惠渠灌區(qū)大田試驗為基礎(chǔ)，分析了不同灌季條件下土壤入滲模型參數(shù)值以及田面綜合糙率系數(shù)的時間變異程度，并利用WinSRFR軟件對其影響下的灌水質(zhì)量評價指標(biāo)以及對灌水過程及灌水效果的影響進行了模擬分析，從而為提高農(nóng)業(yè)灌溉灌水質(zhì)量和灌溉水有效利用程度提供依據(jù)，達到節(jié)約用水、科學(xué)用水的目的。

1 材料與方法

1.1 試驗概況

試驗在陜西省涇陽縣橋底鎮(zhèn)涇惠渠灌區(qū)進行，該地區(qū)土壤類型為塿土，土質(zhì)為粉砂質(zhì)黏壤土，試驗區(qū)面積約2 hm2。作物種植模式為夏玉米-冬小麥輪作(夏玉米種植時間為6月-10月，冬小麥種植時間為10月-6月)冬小麥種植品種為西農(nóng)979，夏玉米種植品種為武科2號。畦灌試驗時間為2012年10月到2015年6月，不同灌水時期按照當(dāng)?shù)毓嗨?jīng)驗進行灌溉。具體試驗布設(shè)圖及指標(biāo)測定方法見文獻[10]及文獻[11]。

1.2 WinSRFR軟件對灌水質(zhì)量的評價

WinSRFR模型是美國水保實驗室開發(fā)的一維地面灌溉模擬模型軟件，它集地面灌溉評價、設(shè)計和模擬為一體[12]。該軟件包含零慣量模型(Zero-Ineria)和運動波模型(Kinematic-wave)。當(dāng)?shù)孛嫫露刃∮?.004時采用零慣量模型(Zero-Ineria)，大于0.004時采用運動波模型[13]。由于涇惠渠灌區(qū)畦灌試驗地塊地面坡度較小，且主要研究尾部閉合條件下的灌水質(zhì)量，所以本文選用零慣量模型進行模擬分析[11]。

模型的輸入?yún)?shù)為：畦長L，畦寬W，畦田坡度i，畦尾開閉的狀態(tài)；入滲參數(shù)K及α，田面糙率系數(shù)n；作物需水量M，灌水流量Q，灌水時間t或改水成數(shù)等。模型輸出結(jié)果包括水流推進、消退曲線、灌水效率、儲水效率及灌水均勻度等。

1.3 灌水質(zhì)量評價指標(biāo)

畦灌灌水質(zhì)量評價指標(biāo)通常包括：灌水效率Ea、灌水均勻度Ed及儲水效率Es，各灌水質(zhì)量評價指標(biāo)表達式如下：

(1)

(3)

式中:Ws為灌后儲存于土壤計劃濕潤層中的水量，mm；Wf為灌入田間的總水量，mm；Zmin為田間入滲水量最少部分田塊的平均入滲水深，m；Zav為整個畦(溝)長的平均入滲深度，m；Wn為計劃濕潤層理論需水量，mm。

2 結(jié)果與分析

2.1 土壤入滲模型參數(shù)值及田面糙率系數(shù)的時間變異性

蔡煥杰等[10]根據(jù)試驗區(qū)實測的水流推進及消退數(shù)據(jù)資料，采用WinSRFR4.1對涇惠渠灌區(qū)考斯加科夫土壤入滲模型兩參數(shù)值及田面綜合糙率系數(shù)值進行了模擬優(yōu)化求解，得出2013-2015年夏玉米-冬小麥輪作體系下不同灌水時期的優(yōu)化參數(shù)值，為了探究夏玉米-冬小麥輪作周期內(nèi)的變異特征，選取2014-2015年進行分析，如表1所示。

引入變異系數(shù)CV評價土壤入滲模型參數(shù)值及綜合糙率系數(shù)的時間變異性。當(dāng)CV≤0.1時，為弱變異性；0.1

表1 涇惠渠灌區(qū)不同灌水時期入滲參數(shù)及綜合糙率系數(shù)優(yōu)化值Tab.1 Soil infiltration parameters and surface roughness coefficients in different irrigation date in Jinghui Canal irrigation area

2.2 時間變異性因子對灌水質(zhì)量的影響

為了探究土壤入滲參數(shù)值及田面綜合糙率系數(shù)等單個時間變異性對灌水質(zhì)量的影響，本文利用WinSRFR4.1軟件對其進行模擬研究。模擬畦田長度及寬度取當(dāng)?shù)剌^為常見的組合，即120 m、4 m；計劃灌水定額取100 mm，改水成數(shù)為當(dāng)?shù)毓嗨畛Ｓ玫?成改水法。采用WinSRFR4.1軟件模擬在不同時間變異性因子作用下，灌水效率Ea、儲水效率Es，灌水均勻度Ed的變化規(guī)律。模擬時控制其他灌水技術(shù)參數(shù)不變，探求單個變異性因子對灌水質(zhì)量指標(biāo)的影響，模擬結(jié)果如圖1所示。

圖1 不同時間變異性因子對畦灌灌水質(zhì)量的影響Fig.1 Effects of different temporal variation factors on border irrigation quality

(1)入滲系數(shù)K。由圖1(a)可見，灌水效率Ea、灌水均勻度Ed和儲水效率Es均隨著入滲系數(shù)值的增大呈現(xiàn)先增大后減小的趨勢。其原因為入滲系數(shù)值較小時，土壤入滲速度小，水流推進速度較快，從而導(dǎo)致儲存在計劃濕潤層中的水量遠小于計劃灌水量，并且在畦尾部形成大量積水，深層滲漏率較大，使得畦尾入滲水量遠大于其他部分，故灌水效率Ea、灌水均勻度Ed和儲水效率Es均較低；隨著入滲系數(shù)值的增大，土壤入滲速度增大，儲存在計劃濕潤層中的水量增大，畦尾部積水逐漸減少，田面受水時間較充分且相對均勻，使得灌水效率Ea、儲水效率Es及灌水均勻度Ed增大；當(dāng)入滲系數(shù)增大到一定程度時，水流推進速度變慢，使畦田上游段入滲水量過大并導(dǎo)致大量深層滲漏但畦尾部灌水不足，使得灌水效率Ea、灌水均勻度Ed、儲水效率Es再次降低。

(2)入滲指數(shù)α。由圖1(b)可知，灌水效率Ea隨著入滲指數(shù)的增大總體呈現(xiàn)增大的趨勢，灌水均勻度Ed呈現(xiàn)減小的趨勢，儲水效率Es呈現(xiàn)先不變后減小的趨勢。其原因是入滲指數(shù)較小時，水流推進慢，再加上畦田受水時間長，土壤入滲量大，使得儲存在計劃濕潤層中的水量超過計劃灌水定額，深層滲漏量較大，故儲水效率Es達到100%，灌水效率Ea較低；隨著入滲指數(shù)的增大，水流運動加快，使得深層滲漏量減少，灌水效率Ea增大，儲存在計劃濕潤層中的水分繼續(xù)減少，使得儲水效率Es不斷減小；8成改水使畦尾部一直存在積水，當(dāng)計劃濕潤層中的水分不斷降低時，相應(yīng)的灌水均勻度Ed也呈現(xiàn)出減小的趨勢。

(3)田面糙率系數(shù)n。由圖1(c)可知，灌水效率Ea隨著田面綜合糙率系數(shù)n的增大呈現(xiàn)減小趨勢，而灌水均勻度Ed呈現(xiàn)先增大后減小的趨勢，儲水效率Es呈現(xiàn)先增大后不變的趨勢。田面糙率系數(shù)較小時，水流推進阻力小，速度快，無深層滲漏率，故具有較高的灌水效率Ea，隨著糙率系數(shù)的增大，水流阻力變大，灌水時間增長，深層滲漏率逐漸變大，故灌水效率Ea呈減小趨勢；在較小的田面糙率情況下，容易導(dǎo)致畦尾部灌水不足，灌水均勻度Ed較低，當(dāng)田面糙率較大時，灌水時間增長，容易導(dǎo)致畦尾產(chǎn)生積水，灌水均勻度Ed減小；同樣的，田面糙率系數(shù)小，灌水時間短，田間灌水不足，儲水效率Es低，當(dāng)田面糙率系數(shù)較大時，灌水時間增大，使得儲存在計劃濕潤層中的水量遠遠大于計劃灌水定額，儲水效率Es達到100%。

2.3 土壤入滲模型參數(shù)值及糙率系數(shù)綜合時間變異性對畦灌灌水的影響

由上述分析得出，土壤入滲模型參數(shù)值及田面糙率系數(shù)等單個時間變異性因子對各灌水質(zhì)量指標(biāo)具有明顯的影響。在生產(chǎn)實踐中，土壤入滲模型參數(shù)值及田面糙率系數(shù)等因子往往是同時變化的，為了研究這種變異性對灌水質(zhì)量的影響作用，有必要更深入地在入滲參數(shù)及田面糙率系數(shù)的綜合變異作用下，研究其對灌水過程及灌水質(zhì)量的影響。由于田間實際灌水過程中，灌水流量很難保持一致，因此利用WinSRFR4.1對涇惠渠灌區(qū)不同灌水時期田塊畦灌過程進行模擬，模擬地塊長度為120 m，寬為4 m，地面坡度為0.003，計劃灌水定額為100 mm，采用8成改水法。圖2為2014-2015年土壤入滲模型參數(shù)值及其田面綜合糙率系數(shù)的時間變異性對畦灌水流推進、水流消退過程及土壤入滲受水時間的影響。

圖2 土壤入滲參數(shù)及田面糙率系數(shù)時間變異性對畦灌過程的影響Fig.2 Effects of Combination temporal variation factors of infiltration parameters and surface roughness coefficient on irrigation process

(1)土壤入滲參數(shù)及田面糙率系數(shù)時間變異性對畦灌過程的影響。圖2(a)～(c)反映了2014-2015年夏玉米-冬小麥輪作周期內(nèi)，土壤入滲模型參數(shù)值及田間糙率系數(shù)的綜合時間變異性對水流推進時間、水流消退時間及土壤受水時間的影響。受土壤入滲參數(shù)及田面糙率系數(shù)的綜合影響，灌溉水流推進到畦尾部的時間差異較明顯，在輪作周期內(nèi)不同灌水時期的推進時間值在15.6～34.3 min之間變化；水流消退時間變化在前期差異性較小，受畦尾部不同積水程度的影響，灌溉水流消退至100 m之后差異性較明顯，直至畦尾部積水消失時消退時間值在9～219.7 min之間劇烈變化；土壤受水時間即灌溉水流消退時間與推進時間的差值，受灌溉水流推進及消退過程的共同作用，輪作周期內(nèi)的差異較明顯。隨著土壤入滲系數(shù)K值的增加，水流在向畦田下游推進的過程中，入滲速率變大，水流推進速度減慢，土壤受水時間減小，入滲系數(shù)K值較小時，入滲速率小，水流推進較快，同時大量灌溉水積累在畦尾，使得水流短時間內(nèi)難以消退；入滲指數(shù)α增大，水流推進時，入滲速率變小，水流推進速度加快，土壤受水時間減少，同時水分入滲量減少，水流消退速度加快；隨著田面糙率系數(shù)n的增加，田面水流推進過程阻力加大，水流推進速度減慢，土壤受水時間增加，水流消退時間加大。

圖2(d)顯示了土壤入滲參數(shù)及田面糙率系數(shù)時間變異性分別對水流推進、消退過程及土壤入滲受水時間沿畦長方向各點的變異系數(shù)值的變化規(guī)律，由圖2(d)可以看出，土壤入滲參數(shù)及田面糙率系數(shù)的時間變異性對水流推進、消退過程及土壤受水時間影響的變異程度均達到了中等變異水平以上。其中，對水流消退過程及土壤入滲受水時間的變異程度較大且在畦尾部均達到了極高的變異程度。其原因主要是，在相同的單寬流量及改水成數(shù)下，不同的入滲參數(shù)及田面糙率系數(shù)的組合導(dǎo)致不同灌水時期土壤水分入滲不均勻性的加劇，土壤入滲參數(shù)大且田面糙率系數(shù)較大時，深層滲漏量較大，并且在畦尾部產(chǎn)生較多的積水，土壤入滲參數(shù)小且田面糙率系數(shù)也較小時，深層滲漏量較少甚至無深層滲漏量，畦尾部產(chǎn)生較少積水，從而使得這種變異程度加劇。

(2)土壤入滲參數(shù)及田面糙率系數(shù)時間變異性對畦灌質(zhì)量的影響。土壤入滲參數(shù)及田面糙率系數(shù)的時間變異性對畦田水流推進、消退及土壤受水時間的影響直接導(dǎo)致土壤入滲水量的分布，從而直接影響到畦灌灌水質(zhì)量。圖3是在相同單寬流量及改水成數(shù)下，對夏玉米-冬小麥輪作周期內(nèi)不同灌水時期模擬的入滲水分在土壤中沿畦長方向的分布圖，可以較為直觀地顯示出灌溉水進入田間后沿畦長方向不同位置的水分分布狀況。在設(shè)定灌水量為100 mm的條件下，2014年11月3日以及2015年1月14日兩次灌水沿畦長方向前半段水分分布明顯不足，而在后半段又造成了較嚴重的深層滲漏；其他3次灌水則在前半段造成了嚴重的深層滲漏，后半段又明顯灌水不足。由于入滲參數(shù)值及田面綜合糙率系數(shù)的時間變異性影響，造成了畦田嚴重的灌水水量分布不均。

圖3 模擬不同灌溉時期土壤入滲水量分布圖Fig.3 Simulated distribution of infiltration water in different irrigation date

表2為2014-2015年模擬的不同灌水時期，相同單寬流量及改水成數(shù)下灌水質(zhì)量指標(biāo)值及統(tǒng)計特征值分布表。由表2可以看出，在其他條件保持一致的條件下，土壤入滲參數(shù)及田面綜合糙率系數(shù)的時間變異性對灌水質(zhì)量指標(biāo)值的影響較大，均呈現(xiàn)中等變異程度。在兩個輪作周期內(nèi)，對儲水效率的影響程度均最大，極差值及變異系數(shù)值分別為46，0.24和51,0.32；對灌水均勻度及灌水效率的影響程度則相對較小，極差值及變異系數(shù)值最大分別為37和0.19。

表2 2014-2015年不同灌水時期灌水質(zhì)量指標(biāo)值及其統(tǒng)計特征值表Tab.2 Simulated values of water quality indexes and its statistical characteristic values in different irrigation date during 2014-2015

3 結(jié) 語

本文基于大田灌水試驗數(shù)據(jù)，分析了不同灌季條件下土壤入滲模型參數(shù)值以及田面綜合糙率系數(shù)的時間變異程度，并利用WinSRFR軟件對其影響下的灌水質(zhì)量評價指標(biāo)以及對灌水過程及灌水效果的影響進行了模擬分析，主要得到以下結(jié)論。

(1)涇惠渠灌區(qū)冬小麥-夏玉米輪作下，不同灌水時期土壤水分入滲模型系數(shù)K、入滲指數(shù)α及田面綜合糙率系數(shù)值n差異性較為顯著，變異系數(shù)值分別為0.27,0.18,0.51.均達到了中等變異程度。

(2)通過模擬比較土壤水分入滲參數(shù)值及田面糙率系數(shù)對灌水質(zhì)量的影響，分析得出入滲系數(shù)、入滲指數(shù)、田面綜合糙率系數(shù)值單個時間變異性因子對灌水效率Ea、灌水均勻度Ed、儲水效率Es均有較大的影響。

(3)不同灌水時期，土壤水分入滲參數(shù)值及田面糙率系數(shù)的綜合變異性作用對畦灌水流推進、消退過程及土壤受水時間均達到了顯著性的影響。

垂向調(diào)控措施的作用區(qū)域主要在0～60 cm土層，其中，40～60 cm土層的土壤含水量及根長密度受調(diào)控影響最顯著。以灌溉水利用效率為衡量指標(biāo)，來評價各處理節(jié)水效果。本研究發(fā)現(xiàn)，過大或過小的計劃濕潤層深度都有可能減小灌溉水利用效率，最優(yōu)調(diào)控方案為：灌水下限為65%田間持水量，苗期計劃濕潤層深度為40 cm，拔節(jié)期計劃濕潤層深度為50 cm，抽穗期至成熟期的計劃濕潤層深度為60 cm。

□

[1] 馮廣龍，劉昌明. 人工控制土壤水分剖面調(diào)控根系分布的研究[J]. 地理學(xué)報，1997，52(5):461-469.

[2] 李 夏. 淺析灌溉對作物根系及產(chǎn)量的影響[J]. 生物災(zāi)害科學(xué)，2014，37(2):188-190.

[3] 馬瑞昆. 供水深度與冬小麥根系發(fā)育的關(guān)系[J]. 干旱地區(qū)農(nóng)業(yè)研究，1991，9(3):1-10.

[4] 馮廣龍，劉昌明，王 立，等. 土壤水分對作物根系生長及分布的調(diào)控作用[J]. 生態(tài)農(nóng)業(yè)研究，1996，4(3):5-9.

[5] 李運生，王 菱，劉士平，等. 土壤-根系界面水分調(diào)控措施對冬小麥根系產(chǎn)量的影響[J]. 生態(tài)學(xué)報， 2002，22(10):1 680-1 687.

[6] Proffit A P B, Berliner P R, Oosterhuis D M. A comparative study of root distribution and water extraction efficiency by wheat grown under high-and low-frequency irrigation[J]. Agron. J. 1985,77:655-662.

[7] 肖俊夫，劉戰(zhàn)東，段愛旺，等. 不同土壤水分條件下冬小麥根系分布規(guī)律及其耗水特性研究[J]. 中國農(nóng)村水利水電，2007，(8):18-21.

[8] 周始威，胡笑濤，王文娥，等. 基于RZWQM模型的石羊河流域春小麥灌溉制度的優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報，2016，32(6):121-129.

[9] 廖榮偉，劉晶淼，安順清，等. 基于微根管技術(shù)的玉米根系生長監(jiān)測[J]. 農(nóng)業(yè)工程學(xué)報，2010，26(10):156-161.

[10] 信乃詮. 計算農(nóng)田蒸發(fā)的水量平衡法[J]. 干旱地區(qū)農(nóng)業(yè)研究，1986，4(2):37-43.

[11] 沈榮開，楊路華，王 康. 關(guān)于以水分生產(chǎn)率作為節(jié)水灌溉指標(biāo)的認識[J]. 中國農(nóng)村水利水電，2001，(5):9-11.

[12] 黨根友，魏亦勤，沈強云，等.灌水次數(shù)對春小麥耗水特性及產(chǎn)量的影響[J].西北農(nóng)業(yè)學(xué)報，2014，23(6)：48-55.

[13] 廖榮偉，劉晶淼，白月明，等. 華北平原冬小麥根系在土壤中的分布研究[J]. 氣象與環(huán)境學(xué)報，2014，30(5):83-88.

[14] 王淑芬，張喜英，裴 冬. 不同供水條件對冬小麥根系分布、產(chǎn)量及水分利用效率的影響[J]. 農(nóng)業(yè)工程學(xué)報，2006，22(2):27-32.

[15] 孔箐鋅，張海林，宋振偉，等. 不同灌溉水平對冬小麥耗水構(gòu)成及利用效率的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究，2009，27(4):37-41.