曹雪松,李和平,鄭和祥,王 軍,馮亞陽
(1.水利部牧區水利科學研究所,呼和浩特 010020;2.中國水利水電科學研究院,北京 100038; 3.內蒙古農業大學,呼和浩特 010018)
水資源不足是中國北方干旱、半干旱地區農業生產的主要限制因素[1,2]。節水灌溉是節約水資源、提高農業用水效率的重要手段[3]。作物需水量是制定流域規劃,地區水利規劃及灌排工程規劃、設計、管理和農田灌排實施的基本依據,在農業生產實踐中占有重要地位[4]。因此,準確地確定一個地區的作物需水量是十分必要的。我國是一個水資源相對缺乏的國家,年用水量約為5 560 億m3,農業用水占72%左右[5],由于工業發展和城市擴大對所需水量大幅度增加,到21世紀中葉,估計農業用水的比例將降至50%左右。然而,我國水資源平均利用率僅為40%,發達國家為80%~90%,單方水生產能力僅為發達國家的42.5%~50%[6],因此對農業用水采用必要的節水措施勢在必行。
近年來,隨著現代節水技術的快速發展,地埋滴灌技術在紫花苜蓿種植上有了較好的應用,但提水的動力都是使用電能。考慮到限制灌溉面積進一步增加的因素,無外乎地區水資源的短缺和灌溉能源動力的短缺,以致使灌溉問題長期得不到解決。隨著太陽能提水技術和現代節水灌溉技術的不斷發展,利用太陽能光伏提水地埋滴灌在西北干旱、半干旱地區具有廣闊的應用前景,對節約能源和水資源有著十分重要的意義。這些地區有著十分豐富的太陽能資源,特別是西北地區,全年平均日照時間3 000 h以上,平均年總輻射量1510 kWh/m2以上[7]。將這些可再生的豐富的太陽能資源與西北地區有限的水資源結合起來開展太陽能提水節水灌溉,使這些地區生態環境的改善、糧食的增產以及畜牧業的發展再次成為可能,這也將會是今后發展節水灌溉的趨勢。
與噴灌和地面灌溉不同,地埋滴灌屬于局部灌溉,因此地埋滴灌條件下作物需水規律有其特殊性。J E Ayars,C J Phene等[8]研究總結了地下滴灌在番茄、甜玉米、紫花苜蓿上的應用,結果顯示地下滴灌不僅能夠提高作物產量,而且能夠提高水分利用率。作物系數是確定作物需水量和制定灌溉制度的重要依據,2006年陳鳳等[9]對陜西楊凌冬小麥和夏玉米的蒸發蒸騰與作物系數進行了研究; 2010年趙娜娜等[10]在北京大興區試驗基地對夏玉米作物系數的計算與耗水量進行了研究;2011年李久生等[11]在內蒙古達拉特旗試驗基地對干旱區玉米滴灌的需水規律進行了研究。目前針對太陽能光伏提水地埋滴灌紫花苜蓿的研究鮮有報道,本文旨在利用太陽能光伏提水地埋滴灌紫花苜蓿條件下的試驗數據,對單作物系數法和雙作物系數法計算的作物需水量進行比較,分析其差異性和原因,以便為準確地計算太陽能光伏提水地埋滴灌紫花苜蓿條件下的作物需水量及制定科學合理的灌溉制度提供依據。
試驗在內蒙古鄂爾多斯市鄂托克前旗昂素鎮哈日根圖嘎查巴圖巴雅爾節水示范戶進行。鄂托克前旗位于內蒙古自治區鄂爾多斯市西南端,海拔1 300~1 400 m,東經106°30′~108°30′,北緯37°38′~38°45′,屬中溫帶半干旱大陸性氣候,年平均氣溫7.9 ℃,年平均降水量261 mm,年平均蒸發量2 498 mm,年平均風速2.6 m/s,年平均日照時數2 958 h;年平均無霜期171 d,最大凍土層深度1.54 m。對試驗區0~40 cm深土壤進行了室內顆粒分析實驗,實驗結果表明土壤容重為1.62 g/cm3,比重為2.71,確定土壤類型為砂土。利用環刀在試驗區取原狀土,進行了田間持水量室內測定,并與田間的測定結果進行對比,確定了試驗區0~40 cm土層的田間持水量為22.86%。采用HOBO地下水位自動測定儀(美國)測定試驗區地下水水位變化,確定了試驗區地下水埋深在1.2~2.0 m。
試驗按灌水定額不同設置3個處理,分別稱之為高水分處理(30 mm)、中水分處理(22.5 mm)和低水分處理(15 mm),每個處理的尺寸為5 m×5 m,每個處理設置3次重復,為了避免每個處理相互之間的影響,每個處理之間設置2 m寬的隔離帶,相同處理不同重復之間相距1 m。試驗處理周圍種植的均為紫花苜蓿。每個處理灌水日期和灌水次數相同,灌水周期均為4~5 d,遇降雨時,灌水日期順延,灌水日期根據中水分處理的適宜含水率下限計算確定,每次灌水的灌水量采用水表計量。各處理的灌水日期與灌水量列入表1。

表1 試驗各處理的灌水日期與灌水量 mm
太陽能光伏提水設備型號為GFJ-3P;太陽能板每塊50 W,共66塊;水泵型號為QS15-26型潛水泵,額定電壓380 V,額定電流6 A,額定功率3.0 kW。試驗均采用以色列生產的貼片式滴灌帶:滴灌帶的壁厚0.4 mm,流量1.2 L/h,滴頭間距0.3 m;每條滴灌帶控制2行紫花苜蓿,滴灌帶間距60 cm,埋設深度為20 cm。
太陽能電池板接收太陽輻射將光能轉換為電能,并帶動水泵工作將水從水源提出,輸送至地埋滴灌系統進行灌溉。試驗裝置見圖1。

圖1 太陽能提水試驗裝置示意圖
試驗區設有HOBOU30型農田氣象站,在作物生育期內對氣溫、降雨量、風速、相對濕度、氣壓、風向等要素進行了觀測。紫花苜蓿全生育期內的降雨日期和有效降雨量列入表2。

表2 紫花苜蓿全生育期內的降雨日期和有效降雨量 mm
(1)紫花苜蓿耗水量采用下式計算:
ETa=P+I-ΔW-Q
(1)
式中:ETa為各時段內耗水量,mm;P為相應時段內的有效降水量,mm;I為相應時段內的灌水量,mm;ΔW為相應時段內的土壤貯水變化量,mm;Q為相應時段內的下邊界水分通量,mm。
(2)水分生產率(Wa)。水分生產率指作物消耗單位水量的產出,其值等于作物產量與作物凈耗水量之比值。作物水分生產率采用下式計算:
(2)
式中:Wa為水分生產率,kg/m3;Y為作物產量,kg/hm2;其他符號同上。
2.2.1 參考作物蒸發蒸騰量計算
參考作物蒸發蒸騰量采用FAO-56推薦的 Penman-Monteith法計算[12]:
(3)
式中:ET0為參照騰發速率, mm/d;Rn為作物表面上的凈輻射,MJ/(m2·d);G為土壤熱通量,MJ/(m2·d);T為2 m高處日平均氣溫,℃;u2為2 m高處的風速,m/s;es為飽合水汽壓,kPa;ea為實際水汽壓,kPa;Δ為飽和水汽壓曲線的傾率;γ為濕度計常數,kPa/℃。
2.2.2 單作物系數法計算
FAO-56中可查到某種作物在標準條件下的作物系數。所謂的標準條件是指無病蟲害,具有最優水土條件,施肥量適宜,在給定氣候條件下能獲全收成的大面積作物系數。FAO-56建議把作物的生育周期劃分為4個階段:生長初期、發育期、生長中期和生長后期。從FAO-56中查得紫花苜蓿的不同生育期作物系數分別為:Kcini= 0.40,Kcmid= 1.20,Kcend= 1.15。
由于試驗區氣候、水土等條件與標準條件下不同,需要對推薦的作物系數進行修正。在生長初期,ET0=3.09 mm/d,從FAO中可得修正后的作物系數為Kcini= 0.42, 推薦的Kcmid和Kcend需要按下式進行修正:
Kcmid=Kcmid(推薦)+
(4)
Kcend=Kcend(推薦)+
(5)
式中:u2為作物生長中、后期2 m高處的日平均風速,m/s;RHmin為作物生長中、后期的日最小相對濕度,%;h為作物生長中、后期的平均株高,m。
2.2.3 雙作物系數法計算
雙作物系數法是把Kc分成兩個系數來表示,一個表征作物蒸騰,即基礎作物系數(Kcb),另一個是表征土壤蒸發的系數(Ke),因此有:
ETc=(Kcb+Ke)ET0
(6)
FAO-56推薦的紫花苜蓿各階段的基礎作物系數分別為Kcbini= 0.30,Kcbmid= 1.15,Kcbend= 1.10。當中、后期的最小相對濕度的平均值RHmin≠ 45%,2 m高處的日平均風速u2≠ 2.0 m/s,并且Kcbend≥0.45時,推薦的Kcbmid和Kcbend需要按下式進行修正:
(7)
式中:u2為作物生長中、后期2 m高處的日平均風速,m/s;RHmin為作物生長中、后期的日最小相對濕度,%;h為作物生長中、后期的平均株高,m。
土壤蒸發系數(Ke)則是用來描述ETc中的土壤蒸發部分,當土壤表面由于降雨或者灌溉較濕潤時,Ke達到最大;當土壤表面干燥時,由于土壤表面沒有可用于蒸發的水分,Ke值很小甚至為零。Ke值用下式計算:
Ke=Kr(Kcmax-Kcb)≤fewKcmax
(8)
Kcmax=max({1.2+[0.04(u2-2)-
0.004(RHmin-45)](h/3)0.3},(Kcb+0.05))
(9)
式中:Kcmax是任何種有作物地表的騰發的上限;Kr是土壤蒸發減小系數;few是裸露與濕潤土壤表面的比值。
其中few可用下式計算:
few=min(1-fc,fw)
(10)

土壤蒸發減小系數Kr是隨著蒸發水量增加而下降,分為2個階段:能量限制階段和蒸發遞減階段。在第一階段,Kr為1;第二階段,Kr按下式計算:
(11)
TEW=1 000(θFC-0.5θWP)Ze
(12)
式中:TEW為表層土壤在完全濕潤時可供蒸發水量的最大深度, mm,取為10 mm;REW為第一階段末的蒸發累計深度,mm,取為6 mm;De,i-1為計算前一天末土壤表層蒸發量的累積深度, mm;θFC為田間持水率,m3/m3;θWP為凋萎含水率,m3/m3;Ze為蒸發深度,一般取0.10~0.15 m,本文取0.10 m。
作物覆蓋的土壤面積比fc隨作物葉面積指數增加而上升,在沒有實測資料時可用下式進行計算:
(13)
式中:Kcmin為沒有地表覆蓋的干燥土壤的最小作物系數,一般取0.15~0.20,本文取0.15。
表3是按式(1)計算出的不同水分處理下的耗水量。表中數據表明,高水分處理的耗水量最大,為477.2 mm,其次是中水分處理,為454.1 mm,低水分處理的耗水量最小,為418.4 mm。隨著灌水量的增加,耗水量也隨之增加。從低水分處理到高水分處理,灌水量增加了一倍,耗水量增加了14.05%,由此可知一次性灌水不宜過多,紫花苜蓿適宜采用中水分處理(22.5 mm)的灌水定額。

表3 不同水分處理下紫花苜蓿的耗水量計算結果 mm
圖2為不同水分處理下第二茬紫花苜蓿株高在生育期內的變化情況。從圖2中可以看出,低水分處理的株高明顯低于其他處理。如7月30日第二茬收獲時高、中、低水分處理下的株高分別為52、46和36 cm,高水分處理的株高分別比中、低水分處理的高出13.04%和44.44%。因此灌水量從中水分處理(22.5 mm)提高到高水分處理(30 mm)對株高影響不太明顯,但當灌水量從中水分處理(22.5 mm)降到低水分處理(15 mm)時,紫花苜蓿的植株高度顯著降低,隨之產量也大幅下降。第一、三茬紫花苜蓿株高在生育期內也有相同的變化趨勢。
表4為不同水分處理下的紫花苜蓿的產量和水分生產率,對表中數據進行分析后可以看出由低水分處理到中水分處理,再到高水分處理,紫花苜蓿總產量分別增加33.75%和37.58%,可知影響紫花苜蓿產量的主要因素是灌水。低、中、高水分處理下紫花苜蓿水分生產率分別為1.73、2.14和2.11 kg/m3,可知中水分處理下紫花苜蓿水分生產率最大。

圖2 不同水分處理下第二茬紫花苜蓿株高在生育期內的變化情況

處理耗水量/mm產量/(kg·hm-2)第一茬第二茬第三茬總產量水分生產率/(kg·m-3)高水分477.23304.53862.5289810065.02.11中水分454.13055.53864.027819700.52.14低水分418.42281.52757.022147252.51.73
圖3是用單、雙作物系數法計算出的地埋滴灌紫花苜蓿作物需水量與田間實測值在全生育期的逐日變化過程(圖3中P+I表示的是降雨和灌溉)。從圖中可以看出,紫花苜蓿全生育期的作物需水量很不均衡,第二茬紫花苜蓿需水量最大,第一茬和第三茬需水量相對較小,但兩者相差不大。經過圖中計算出的ETC與實測值進行比較表明,全生育期大多數情況下用雙作物系數法計算出的紫花苜蓿蒸發蒸騰量ETC值和實測值比較接近,變化趨勢亦相同,而用單作物系數法計算出的紫花苜蓿蒸發蒸騰量ETC值和實測值存在一定的偏差。由于實測值會受偶然因素的干擾,圖3中個別實測值與計算值相比波動較大,計算值則比較穩定。在生長初期,用單、雙作物系數法計算的ETC值比實測值稍低;在其他三個生長階段,計算的ETC值和實測值很接近,特別是在發育期和生長中期,用雙作物系數法計算的ETC值與實測值吻合得更好。

圖3 紫花苜蓿全生育期單、雙作物系數計算的ETC與實測ETC
表5分別是用單、雙作物系數法計算出來的地埋滴灌紫花苜蓿不同生育階段作物需水量和實測值的絕對偏差以及相對偏差。從表中可以看出,用雙作物系數法比單作物系數法得出的作物需水量更加接近實測值,第一茬、第二茬和第三茬的絕對偏差分別是6.71、-3.05和3.73 mm,相對偏差分別為4.72%、-1.82%和2.51%,而用單作物系數法計算出來的作物需水量與實測值相比偏差相對較大,從全生育期來講,用雙作物系數法計算出來的ETC值更加接近田間實測值,說明雙作物系數法可以更好的描述地埋滴灌紫花苜蓿條件下降雨或者灌溉后ETC的變化。根據FAO-56 推薦的Penman-Monteith法以及單、雙作物系數法對地埋滴灌紫花苜蓿需水量的計算結果,確定了研究區地埋滴灌紫花苜蓿全生育期的作物需水量為460 mm左右,其中第二茬作物需水量最大,第一茬和第三茬作物需水量較小,但兩者相差不大。

表5 單、雙作物系數法根據修正值計算的作物需水量與實測值比較
(1)通過將FAO提供的作物系數進行修正,用單、雙作物系數法分別對研究區紫花苜蓿的需水量進行計算,并且與作物實測的騰發量進行比較發現,用雙作物系數法計算的作物需水量更切合實際,建議采用雙作物系數法計算該地區紫花苜蓿的作物需水量。
(2)研究區作物耗水量主要取決于灌水量,耗水量隨著灌水量的增加而增大,當灌水量從260.4 mm增加到506.9 mm(增加94.7%)時,作物耗水量從418.4 mm增加到477.2 mm(增加14.1%)。
(3)紫花苜蓿的株高和產量都隨灌水量的增加而增大。低水分處理與中水分處理相比(灌水量從260.4 mm增加到386.9 mm,增加48.6%),產量從7 252.5 kg/hm2增加到9 700.5 kg/hm2,增加33.8%;中水分處理與高水分處理相比(灌水量從386.9 mm增加到506.9 mm,增加31.0%),產量從9 700.5 kg/hm2增加到10 065 kg/hm2,增加3.8%。從節水并且獲得高產的角度出發,建議采用中水分處理(22.5 mm)的灌溉制度,即灌水定額控制在22.5 mm左右,灌水周期4~5 d,遇到降水時灌水日期順延,整個生育期的需水量為460 mm。
(1)太陽能光伏發電技術對天氣的依賴性較強。在晴天太陽輻射較強時,太陽能電池發電充足,水泵運行時間較長,抽水量較大,能夠滿足滴灌需求。而在多云或者陰天時,太陽能電池發電不足以滿足水泵的工作需要。
(2)地埋滴灌作為一種非充分灌溉方式,其灌水均勻性是影響灌水質量的重要指標,并最終影響水分生產率和紫花苜蓿產量。因此有必要開展太陽能光伏提水地埋滴灌紫花苜蓿灌水均勻性研究。
(3)由于紫花苜蓿屬于多年生牧草,其不同生長階段的根系層分布也不同,因此有必要對滴灌帶不同埋設深度做進一步的研究。
(4)由于地埋滴灌是將滴灌帶埋于地下,土壤顆粒有可能進入滴頭處導致滴灌帶滴頭堵塞,從而影響灌水的均勻性,因此地埋滴灌的堵塞問題有待進一步研究。
(5)為了使紫花苜蓿得到高產,在其生長的不同階段可以使用一定的肥料,將肥料隨滴灌帶直接作用于紫花苜蓿根系部位,實現水肥一體化,有待做進一步的研究。
[1] 孫景生,肖俊夫,張寄陽,等. 夏玉米產量與水分關系及其高效用水灌溉制度[J]. 灌溉排水,1998,17(3):17-21.
[2] Ziaei A N, Sepaskhah A R. Model for simulation of winter wheat yield under dry land and irrigated conditions[J]. Agricultural Water Management. 2003,58(1):1-17.
[3] 康紹忠,蔡煥杰,馮紹元. 現代農業與生態節水的技術創新與未來研究重點[J]. 農業工程學報,2004,20(1):1-6.
[4] 樊引琴,蔡煥杰. 單作物系數法和雙作物系數法計算作物需水量的比較研究[J]. 水利學報,2002,(3):50-54.
[5] 張啟舜.從國際水問題看我國節水灌溉革命[J].科技導報,2000,(8):51-54.
[6] 黃修橋.關于我國節水灌溉宏觀發展的幾點思考[J].灌溉排水,1998,17(4):53-57.
[7] 劉 偉,吳永忠. 太陽能光伏提水節水灌溉技術[J]. 節水灌溉,2004,(5):49-50.
[8] J E Ayars,C J Phene,R B Hutmacher,et al. Subsurface drip irrigation of row crops: a review of 15 years of research at the Water Management Research Laboratory[J]. Agricultural Water Management,1999,(42):1-27.
[9] 陳 鳳,蔡煥杰,王 健,等.楊凌地區冬小麥和夏玉米蒸發蒸騰和作物系數的確定[J].農業工程學報,2006,22(5): 191-193.
[10] 趙娜娜,劉 鈺.夏玉米作物系數計算與耗水量研究[J].水利學報,2010,(8):0 953-0 959.
[11] 李久生,饒敏杰.干旱區玉米滴灌需水規律的田間試驗研究[J].灌溉排水學報,2003,22(1):16-21.
[12] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration.Guidelines or Computing Crop Water Requirements[C]∥ Irrigation and Drainage Paper No.56, FAO, Rome, 1998.