基于多年降水的豫東地區夏玉米灌溉制度優化

孫晉鍇，馮躍華，張子敬，高子樂(河南省豫東水利工程管理局惠北水利科學試驗站，河南 開封 475002)

水，是一種基礎性自然資源，也是是一種公共性社會資源，在今天更是一種重要的戰略資源。一直以來，我國對于水資源的利用一直處于高消耗、低效率、低效益的狀態，再加上我國降水時空分布不均，水資源對農業產生的制約效應越來越大。由于可用于灌溉的水資源是有限的，在有限的灌溉水量條件下，通過對自然降水和玉米各生育期的需水量進行耦合，在非充分灌溉理論研究指導下，提高灌水利用效率，使夏玉米灌溉制度進一步優化，從而達到用更少的水生產更多的糧食的目的。

節水灌溉的途徑有很多，其中最重要的是灌溉制度的優化。灌溉制度的優化的早期方案主要是動態規劃法[1]，但動態規劃法有一個缺點，就是隨著變量的增加、劃分段數增多，其解法愈加復雜并且不能保證所得解為最優。付強等采用改進的RAGA算法與多維動態規劃相結合，以灌水量作為優化變量提出了RAGA-DP模型，解決了反復試算、局部最優解與全局最優解不同、未涉及具體灌溉日期等情況[2]。郄志紅等通過對NSGA-Ⅱ算法進行改進，采用遺傳編碼設計并以灌溉日期和灌水量作為決策變量，對相對產量及灌水總量進行優化，建立了同時對灌水日期與灌水量進行的多目標優化模型[3]。上述研究多為靜態研究，將降雨作為固定值考慮，即針對特定降水年的降水情況進行優化，對于其余年份一般并不適用，這在某種程度上限制了其在現實生活中的應用與推廣。

節水農作制度的研究從20世紀80年代以來，國家對其重視程度越來越高，得出了很多不同的理論，通過理論聯系實際，在實踐中又獲取了很多的經驗。不同作物的耗水量與有效降水量時空分布的耦合程度存在顯著差異[4-7]，這是節水型種植結構調整的依據，根據當地水資源條件選取合適的農作物，增加對雨水的利用，可有效節水。這種情況我國學術界一般認為，不同的農業種植結構有不同的需水量要求，應充分考慮當地的水資源與自然情況，結合實際，選用作物生育期與自然降水吻合程度好的、需水少的、對水分利用率高的品種。

本文在惠北水利科學試驗站多年降雨觀測數據支持下，通過對降水進行綜合分析并對灌溉水量、灌溉日期優化，進一步和當地玉米生育周期及各生育期需水量進行耦合，取得最優解，以期得到適應性較強的灌溉制度。

1 試驗方案與方法

1.1 試驗方案及地區背景

本次試驗在河南省開封市祥符區惠北水利科學試驗站節水試驗基地的測坑及大田進行，6-9月是河南省豫東地區的雨季，也是夏玉米生長發育的季節。豫東地區此時由于海洋季風控制，各地普遍多有，此階段的降雨量占到全年降雨量的50%，夏玉米生育期總降雨量可達300～800 mm，最大年降水量為1 180 mm，最小年降水量為179.2 mm，多年平均降水量為659 mm，變率為21%。豫東地區夏玉米生育期太陽輻射總量為1 535.4～2 435.1 MJ/m2，日照時數為650～750 h，全生育期積溫2 375～2 741 ℃。降雨指標基本上可以達到玉米生育期所需水量，積溫也達到中早熟品種的指標。降雨和高溫同時發生，雖然在一定程度上提高了玉米對于水分、熱量和光能輻射的綜合利用效率，但是，由于豫東地區的降水時空分布與玉米的生長周期不協調，往往導致玉米整個生育期需水期與自然降雨的不協調，需水的時候旱、該旱的時候澇。

影響豫東地區夏玉米產量的主要氣象災害有：5月下旬至6月中上旬，正是夏玉米播種的季節，此時豫東地區干旱的頻率很高，從而導致出苗水不足，如果只靠自然降水補水，將會造成晚播減產；在夏玉米孕穗、抽雄和開花吐絲期，正是玉米生長發育需水最多的時期，此時一般為7月下旬到8月上旬，尤其最近幾年豫東地區此階段往往大面積干旱，如果不及時灌水，將造成玉米大量缺粒、禿尖，阻礙玉米的灌漿從而導致產量下降。

夏玉米需水量及灌水制度試驗在測坑進行，灌水處理設計方案如表1和表2所示。

表1 玉米需水非關鍵期測坑田間持水量控制指標表Tab.1 Corn water requirement of non-critical test pit fieldwater holding capacity control indicator table

表2 玉米需水關鍵期測坑田間持水量控制指標表Tab.2 Corn water requirement of key indicators pit fieldwater holding capacity control table

注：田間持水量控制指標為土壤最大持水量的百分比。

1.2 主要觀測內容

(1)降水、蒸發、風速、日照、氣溫、水汽壓等氣象要素。

(2)土壤含水率觀測。

(3)夏玉米灌水量觀測。

(4)葉面積指數、株高、頸粗等生態指標調查。

(5)夏玉米經濟產量、生物產量、百粒重等考種內容。

2 農田水分狀態模擬與相對產量模型

(1)基于水量平衡的農作物耗水量公式。農作物的耗水量是參考作物蒸散量來確定的，作物耗水量的計算公式為：

W=R+Q+ΔW

(1)

ΔW=1.408H×1×(W1-W2)

(2)

式中：W為玉米生育期內的耗水量，m3/hm2；R為生育期內的有效降雨量，m3/hm2；Q為生育期內的灌水量，m3/hm2； △W為土壤水變化量，m3/hm2；H為計算深度，m；W1、W2為土壤含水率觀測初、末含水率，%。

降雨前會將測坑上方的擋雨棚關閉，所以上式中R為0；另外，測坑為有底測坑，不考慮地下水利用量。

(2)基于作物耗水量的作物產量模型。目前，在國內作物耗水量與作物產量的研究中多采用二次拋物線模型：

Y=aET2c+bET2c+c

(3)

式中：Y為作物產量，kg/hm2；ETc為作物耗水量，mm；a、b、c為回歸系數。

(3)基于水量平衡的農田水分狀況的動態模擬。

Si,j=Si-1,j-Ei,j-Ti,j+Pi,j+Ii,j-Di,j-Qi,j

(4)

式中：Si-1,j、Si,j分別為i時段初、末土層j土壤含水量；Ei,j為i時段土層j的蒸發量，mm；Ti,j為i時段土層j的蒸騰量，mm；Pi,j為i時段進入土層j的降雨量，mm；Ii,j為i時段進入土層j的灌溉水量，mm；Di,j為i時段土層j的排水量，mm；Qi,j為根系層底部的水分交換量,以向下滲漏為正，mm。

(4)相對產量及優化模型。本文采用應用較普遍的Jensen乘法模型：

(5)

式中:y為夏玉米的實際產量，kg；ym為夏玉米理論最高產量，kg；λi為夏玉米在i生育期的缺水敏感指數Ei、EMi分別為夏玉米在i生育階段的實際蒸發量和潛在蒸發量，mm；Ti、TMi分別為夏玉米i生育階段的實際蒸騰量和潛在蒸騰量，mm；N為夏玉米生育階段數。

(5)基于改進的GNSGA-Ⅱ與多年降雨資料的灌溉制度優化模型。

(6)

約束條件：

灌溉日期約束：

dmin≤di≤dmax

(7)

土壤計劃濕潤層深度的約束：

Hmin≤Hi≤Hmax

(8)

式中：f1為多年夏玉米相對產量均值；f2為多年夏玉米全生育期的灌溉水量，mm；f3為多年夏玉米相對產量方差；xi為每次灌水的灌水量，mm；dmin、dmax分別為生育期的最小及最大天數；di為從生育期開始到第i次灌水時的天數；Hmin、Hmax為第i生長階段內允許的土壤最小及最大濕潤深度。

3 優化結果與分析

設玉米全生育階段最多可以進行N次灌水，再假設每次灌水量xi以及灌水間隔時間ti作為優化參數，則待優化的決策變量Y=(xi，ti)。采用基于改進的GNSGA-Ⅱ的作物灌溉制度優化的多目標優化算法對模型進行求解。

3.1 基于惠北水利科學試驗站36 a降雨直接優化的結果

根據惠北水利科學試驗站36 a降雨資料， 按優化公式(6)～(8)進行優化， 其結果如圖1和圖2所示。由圖1可知灌水總量達到100 mm時玉米的相對產量達到97%，其后隨著灌水量增加、玉米相對產量上升趨于平緩。與此同時根據圖2的線條趨勢可以看出玉米相對產量的均方差在減小，說明數據的離散性、波動性在降低。當灌水量達到160 mm時，相對產量均方差數量級達到10-5，說明相對產量為零的情況基本不存在。

圖1 相對產量關系優化結果Fig.1 Optimization results of relative yield for the maize

圖2 均方差關系的優化結果Fig.2 Optimization results of relative yield variancefor the maize

3.2 基于豫東地區多年氣象資料和玉米各生育期耗水量優化灌溉制度

夏玉米在不同生育階段對水分的需求有很大不同，苗期至拔節期占全生育期需水量的17.38%～24.43%；拔節期至抽雄期占29.29%～32.86%；抽雄期至灌漿期27.04%～32.98%；灌漿期至成熟占12.15%～15.17%。夏玉米生育期間需水強度的變化趨勢為：苗期-拔節期、灌漿期-成熟期需水強度較小，拔節期-抽雄期，抽雄期-灌漿期需水強度較大，即先升后降。

夏玉米生育期與當地降雨的耦合程度，是評價玉米在當地生態適應型的重要指標，也對優化玉米的節水灌溉制度，節約用水，提高產量具有重要意義。當玉米與當地降雨耦合程度較好時，即便不用澆水，自然降水就能滿足玉米在各個生長階段的需水量。如果玉米與當地降雨耦合不好時，不僅是會需要大量灌水才能滿足玉米生長發育的需求，甚至有時由于水量過大，還需排水保收。

夏玉米生育期與當地降水時序的吻合既表現在滿足玉米全生育期的用水需求，又表現在對于玉米各生育期的用水滿足。從圖3可以看出，豫東地區在玉米各生育期的降水和玉米植株在各生育期的需水量并不耦合。從基于多年降雨優化結果中選取灌溉制度非劣解與豫東地區玉米生育期多年降水平均值相耦合取得最優解集,并在最優解集中確定相對適合農民灌水控制的解為最終解。其優化后的灌溉制度如表3所示。

圖3 玉米各生育階段需水量和當地降水量Fig.3 Water demand and precipitation of maize in different growing stages

灌水量灌水時期/mm總灌水量/mm多年相對產量均值多年相對產量相對均方差8575拔節期至抽雄期抽雄期至灌漿期1600.9823.13×10-5

4 結 語

(1)根據已知的36 a降雨資料，根據灌水量和灌水日期進行優化，得出相對于不同灌水量的多個非劣解，即不同的灌溉制度。隨灌水總量的變化，相對產量及相對產量的均方差隨之變化，曲線平滑、趨勢性強，表明這種灌溉制度適應性較強。

(2)基于既有多年降雨條件的非劣解與豫東地區的玉米生育期耗水量和豫東地區夏玉米生育期降水情況相結合，從優化灌溉制度的非劣解中選取與豫東地區當地氣象條件和玉米生育情況耦合較好的解作為最優解，即在拔節期至抽雄期、抽雄期至灌漿期進行兩次灌水，灌水定額分別為85和75 mm，可以使相對產量達到98.2%。

[1] 操信春, 楊陳玉, 何 鑫,等. 中國灌溉水資源利用效率的空間差異分析[J]. 中國農村水利水電, 2016,(8):128-132.

[2] 付 強，王立坤，門寶輝，等．推求水稻非充分灌溉下優化灌溉制度的新方法----基于實碼加速遺傳算法的多維動態規劃法[J]．水利學報，2003，34(1):123-128.

[3] 郄志紅，韓李明，吳鑫淼．基于改進 NSGA-Ⅱ的作物灌水量與灌溉日期同步優化[J]．農業機械學報，2011，42(5):106-110.

[4] 趙俊芳, 蒲菲癱, 閆偉兄,等. 基于APSIM模型識別氣象因子對內蒙春小麥潛在產量的影響[J]. 生態學雜志, 2017,36(3):757-765.

[5] 鄭鳳杰, 張 超. 玉米生長階段降水分布影響分析[J]. 北方農業學報, 2016, 44(6):82-85.

[6] 趙福成, 譚禾平, 包 斐,等. 作物生長模擬模型及其在玉米上的應用綜述[J]. 農業科技通訊, 2016,(11):12-14.

[7] 懷寶東, 隋文志. 水稻節水灌溉技術試驗研究[J]. 吉林農業, 2017,(2):74.