楊 晶，段 萌，張寶忠，段喜明

（1.山西農(nóng)業(yè)大學林學院，山西晉中 030801；2.中國水利水電科學研究院流域水循環(huán)模擬與調(diào)控國家重點實驗室，北京 100038）

0 引 言

優(yōu)化灌溉制度是根據(jù)作物的需水特征和缺水敏感因子、以及當?shù)氐臍夂?、土壤、農(nóng)業(yè)和灌溉技術等方面在灌溉水量有限時，制定的最佳灌水方案[1]。由于可用于灌溉的水資源日益稀缺，在作物生長階段實行非充分灌溉，可在一定程度上提高水分生產(chǎn)率[2]。

優(yōu)化灌溉制度是我國北方農(nóng)業(yè)節(jié)水的研究重點。田間試驗得到的結果較為可靠，但因周期長，投入比較高，且受各種不確定因素的影響，實用性較差；模型研究可降低試驗成本，縮短試驗周期，減少各種因素的干擾。與其他作物生長模型相比，Aquacrop 模型具有輸入?yún)?shù)少、操作簡便、模擬精度尚可的特點[3]。為保證模型模擬結果的可靠性，需利用田間試驗結果進行模型精度的驗證與評價，常用做優(yōu)化灌溉制度的手段。

以往研究多集中在冬小麥和夏玉米單作模式下優(yōu)化灌溉制度，較少在冬小麥-夏玉米輪作（簡稱“麥玉輪作”）模式下優(yōu)化灌溉制度得到耗水規(guī)律和產(chǎn)量分析[4]，并提出合理的灌水措施，并且較少考慮不同灌溉制度下的土壤水分變化[5]。為此，本文以Aquacrop 模型為技術支持，通過田間觀測結果來校驗作物生長模型，以確保模型精度。在輪作模式下對冬小麥和夏玉米的灌溉制度進行優(yōu)化和分析，從而為華北地區(qū)的農(nóng)業(yè)高效灌溉提供科學依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

本試驗在中國水科院大興節(jié)水灌溉試驗研究基地（39°37′25″N,116°25′51″E）進行[6]，多年平均氣溫12.1 ℃，全年日照時數(shù)約2 600 h，屬半干旱大陸性季風氣候，地下水深10 m。

1.2 觀測項目

模型的輸入模塊主要包括作物、氣象、土壤、地下水、田間管理和初始條件等。根據(jù)試驗站的大氣平均溫度、平均相對濕度、飽和水氣壓差、太陽凈輻射，平均氣壓，2 m 高度處風速，采用FAO 推薦的Penman-Monteith 公式計算參照作物騰發(fā)量ET0[7]。土壤體積含水率采用TRIMER-T3/IPH 型土壤水分測定儀測定，每隔3～5d 在平均間隔為10 cm 且1.1 m 的深度的Trime 管進行測量，降雨和灌溉后加測一次。冬小麥和夏玉米出苗后，每10 d 測定一次地上生物量，在各小區(qū)（面積為5.5 m×5.5 m）選取有代表性3株植株，從莖基部剪下，獲得完整的冠部，烘干后稱量，并根據(jù)同期分蘗情況，推算單位面積的生物量[8]。收獲時考種測產(chǎn)，測定各小區(qū)經(jīng)自然風干后籽粒質(zhì)量，并折算成每公頃的產(chǎn)量[9]。葉面積指數(shù)(LAI)平均每5 d觀測一次，選取各小區(qū)內(nèi)生長均勻、具有代表性的3株植物，對其株高、葉面積進行標記，并進行定點觀察，用鋼尺測量了3株有效葉片的長度和寬度，并計算了每株有效葉面積的平均面積，將其乘以折算系數(shù)K，除以株距與行距。冬小麥和夏玉米的冠層覆蓋度(C)通過葉面積指數(shù)獲得，計算公式如下[10]：

土壤數(shù)據(jù)來源于2020年大興試驗站的土壤粒徑分析，試驗區(qū)土質(zhì)以砂壤土為主，土壤理化特性見表1[11]。

表1 大興試驗站土壤特性Tab.1 Soil characteristics of Daxing experiment station

1.3 AquaCrop模型

利用實測的作物生長數(shù)據(jù)、氣象數(shù)據(jù)、土壤數(shù)據(jù)、灌溉制度、田間管理措施等資料建立模型基礎數(shù)據(jù)庫，分別選擇土壤含水率、冠層覆蓋度、生物量與產(chǎn)量等指標對模型參數(shù)進行校驗。其中選擇2019年10月-2020年10月麥玉生育期間的實測數(shù)據(jù)對模型參數(shù)進行率定；利用2020年10月-2021年10月的實測數(shù)據(jù)對模型進行驗證。

1.4 情景方案設置

本文通過優(yōu)化適線法擬合Pearson-Ⅲ型曲線[12]，分析北京站1981-2021年作物生育期內(nèi)的累積降雨量，確定氣象典型年豐水年、平水年和枯水年分別為1997年、2009年和1987年，降雨量分別為690 mm、527 mm 和432 mm。根據(jù)土壤水分占田間持水量（FC）的不同進行灌溉制度情景設定，以耗水量、產(chǎn)量和水分生產(chǎn)率為目標函數(shù)，調(diào)整灌水量、灌水時間和灌水次數(shù)對灌溉制度進行優(yōu)化決策。設計土壤水分下限達到70%FC、60%FC、50%FC時，根據(jù)查閱文獻和當?shù)毓喔冉?jīng)驗設計進行灌溉的灌水定額分別為60 mm、80 mm、100 mm。麥玉輪作的灌水時間階段可分為：冬灌A，返青-拔節(jié)B，拔節(jié)-抽穗C，抽穗-灌漿D，灌漿-乳熟E，玉米苗期F。設置38 種灌溉情景（T1～T38）如表2所示。

表2 麥玉輪作灌溉制度情景方案Tab.2 Scenarios of winter wheat-summer maize rotation irrigation system

2 結果與分析

2.1 模型率定與驗證

本文主要采用均方根誤差(RMSE)、決定系數(shù)（R2）、Nash效率系數(shù)(EF)指標進行模型精度評價。運用率定后的AquaCrop模型用來進行灌溉制度優(yōu)化。麥玉生育期內(nèi)土壤墑情、冠層覆蓋度的動態(tài)變化過程及其生物量和產(chǎn)量見圖1、圖2和表3；模型校驗后的參數(shù)見表4。

表3 模型生產(chǎn)潛力率定驗證評價指標Tab.3 Model production potential calibration verification evaluation index

表4 麥玉AquaCrop模型參數(shù)列表Tab.4 Parameter list of AquaCrop model for winter wheat and summer maize

圖1 麥玉土壤含水率模擬值與實測值比較Fig.1 Comparison between simulated and measured soil moisture content of wheat and maize

圖2 麥玉冠層覆蓋度模擬值與實測值比較Fig.2 Comparison between simulated and measured canopy coverage of wheat and maize

2.2 不同灌溉制度下麥玉輪作的特征指標變化規(guī)律

2.2.1 總耗水量

不同灌溉制度下麥玉輪作總耗水量見圖3。在平水年條件下，灌溉定額為0（雨養(yǎng)條件下）、180～240 mm、300～320 mm、360～400 mm 時的總耗水量分別為555.9 mm、630.2～707.7 mm、738.7～772.4 mm、770.3～796.8 mm。

圖3 不同灌溉制度下麥玉輪作的總耗水量Fig.3 The total water consumption of wheat-maize rotation under different irrigation systems

灌溉定額為180～240 mm、300～320 mm、360～370 mm 時，對應的最高和最低總耗水量分別為T30（707.7 mm）和T17（630.2 mm）、T33 （772.4 mm） 和T31 （738.7 mm）、T38（800.2 mm）和T1（773.3 mm）。

2.2.2 總產(chǎn)量

不同灌溉制度下麥玉輪作的總產(chǎn)量見圖4。在平水年條件下，灌溉定額為0（雨養(yǎng)條件下）、180～240 mm、300～320 mm、360～400 mm 時對應的麥玉總產(chǎn)量分別為8 599 kg/hm2、10 241～14 681 kg/hm2、15 417～16 817 kg/hm2、16 053～16 119 kg/hm2?？梢娍偖a(chǎn)量隨灌溉量的增加而增大。

圖4 不同灌溉制度下麥玉輪作的總產(chǎn)量Fig.4 The total yield of wheat-maize rotation under different irrigation systems

灌溉定額為180～240 mm、300～320 mm、360～370 mm 時，對應的最高和最低總產(chǎn)量分別為T9（14 681 kg/hm2）和T17（10 241 kg/hm2）、T6（16 817 kg/hm2）和T35（15 417 kg/hm2）、T38（16 898 kg/hm2）和T36（16 807 kg/hm2）。

2.2.3 總水分生產(chǎn)率

不同灌溉制度下麥玉輪作的總水分生產(chǎn)率見圖5。在平水年條件下，灌溉定額為0（雨養(yǎng)條件下）、180～240 mm、300～320 mm、360～370 mm 時，對應的總水分生產(chǎn)率分別為1.54 kg/m3、1.62～2.14 kg/m3、2.00～2.22 kg/m3、2.08～2.11 kg/m3。

圖5 不同灌溉制度下麥玉輪作的總水分生產(chǎn)率Fig.5 Total water productivity of wheat-maize rotation under different irrigation systems

2.2.4 總土壤水分變化

不同灌溉制度下冬小麥-夏玉米輪作全年總土壤水分變化見圖6，與其他有灌溉的處理相比，雨養(yǎng)條件下冬小麥-夏玉米的土壤水分變化消耗較大，其水分變化量在69.4 mm 左右。在同一處理條件下，因典型水文年內(nèi)降水分布的差異，輪作下的土壤水分在枯水年條件下與豐水年和平水年相比變化較大，當在平水年條件下，灌溉定額為180～240 mm 下各處理土壤水分變化在63.3～93.7 mm 范圍內(nèi)；灌溉定額為300～320 mm 下各處理土壤水分變化在-23.8～43.7 mm 范圍內(nèi)；灌溉定額為360～370 mm下各處理土壤水分變化在-32.0～37.2 mm范圍內(nèi)；可知，總體上輪作下土壤水分變化量隨灌溉量的增大而減小。

圖6 不同灌溉制度下麥玉輪作的總土壤水分變化Fig.6 Changes of total soil moisture in wheat-maize rotation under different irrigation systems

2.3 大興區(qū)麥玉輪作時灌溉制度的優(yōu)化

在確保麥玉輪作穩(wěn)產(chǎn)條件下，平水年型不同灌溉制度如表5所示。T23處理（灌溉定額為180～240 mm）的產(chǎn)量和水分生產(chǎn)率偏低；而T38 處理（灌溉定額為320～370 mm）時灌溉量偏大；T5和T6水分生產(chǎn)率相對其他情景方案較高，T5處理（灌溉定額為300 mm）與T6 處理（灌溉定額320 mm）相比，灌溉量減少20 mm，但仍能保證99.5%產(chǎn)量。由此遴選出試驗區(qū)冬小麥-夏玉米輪作時，采用T5處理（灌溉定額為300 mm）的優(yōu)化灌溉制度。

表5 推薦冬小麥-夏玉米輪作灌溉制度Tab.5 Recommended wheat-maize rotation irrigation system

3 討 論

（1）農(nóng)業(yè)水資源的稀缺使得大家日益重視農(nóng)業(yè)用水，灌溉制度的研究成為熱點[13,14]，而前人研究大多集中于作物單作模式下灌溉用水更少或追求高產(chǎn)這一方向，較少研究輪作模式下的灌溉制度研究。本文認為作物輪作模式下優(yōu)化灌溉制度能夠更加節(jié)水，而單作模式下優(yōu)化灌溉制度，產(chǎn)量較高但全年整體來看并不省水。輪作時灌溉定額小于300 mm 的耗水量、產(chǎn)量會隨著灌溉定額的增加而增加，這與其他學者的研究成果相一致[15]。

（2）作物單作時僅在生育期內(nèi)優(yōu)化灌溉制度看似省水，但實際上是消耗掉了土壤原本儲存的水分；而作物輪作時優(yōu)化灌溉制度，卻能使水分生產(chǎn)率和產(chǎn)量提高。以麥玉輪作為例[16]，玉米并未進行播前灌溉，而是利用小麥生育期末的灌水和土壤水分，作為玉米的底墑水，這一優(yōu)化灌溉制度能為該地區(qū)農(nóng)業(yè)用水提供指導。

4 結 論

（1）AquaCrop 模型可較好地模擬不同灌溉制度下，麥玉輪作模式的土壤水分和冠層覆蓋度，以及生物量和產(chǎn)量的動態(tài)變化過程，模型的模擬值與觀測值之間的RMSE均小于20%，R2和EF分別大于0.7和0.5，模擬誤差較小。

（2）利用驗證后的模型模擬麥玉輪作灌溉定額與耗水量、產(chǎn)量和水分生產(chǎn)率之間的響應關系，確定和評估麥玉輪作模式下的生產(chǎn)潛力。結果表明，在輪作條件下且灌水量小于300 mm 時，耗水量、產(chǎn)量和水分生產(chǎn)率隨灌溉量的增大而增加，而土壤水分變化量隨灌溉量的增大而減小。灌溉定額和關鍵生育期的選擇，對于保證總年產(chǎn)量和節(jié)水至關重要。

（3）當前華北地區(qū)嚴重缺水，應堅持量水發(fā)展，大力推廣虧缺灌溉制度，建議采用麥玉輪作種植模式，且灌溉定額為300 mm，保障冬小麥在返青-拔節(jié)、拔節(jié)-抽穗、抽穗-灌漿、灌漿-乳熟（T5：300 mm），夏玉米苗期階段灌水，將灌漿-乳熟階段的補充灌水作為夏玉米底墑水，用于滿足玉米種子發(fā)芽的土壤含水率，此時得到的總產(chǎn)量為16 735 kg/hm2，總耗水量為758.2 mm，水分生產(chǎn)率為2.21 kg/m3，因此輪作模式下的該灌溉制度較為節(jié)水并且產(chǎn)量能得到保障。