限水灌溉下冬小麥最佳灌溉施肥制度研究

杜娟娟，李粉嬋，王仰仁

(1.山西省水利水電科學研究院, 太原 030002; 2. 天津農學院水利工程系， 天津 300384)

0 引 言

冬小麥在我國具有悠久的種植歷史，也是目前我國的主要糧食作物之一，在現行糧食結構中占有較大的比例[1]。光、熱、水、氣以及養分是冬小麥生長發育過程中的重要因子，其中水分、養分是冬小麥生產過程中的重要資源，其獲取小部分依靠自然資源，大部分依靠相應的田間管理的供給，如灌水和施肥等[2]。但是，長期以來我國在灌溉和施肥方面仍存在較大的問題，不合理的灌溉施肥方式不僅會造成資源的浪費，也會導致水分及肥料利用率的下降，同時有可能造成生態環境的惡化。因此，在我國目前灌溉供水量普遍不能滿足作物需水要求的現狀下，建立科學的灌溉施肥制度，改變傳統的灌溉施肥方式是現代農業發展的必然趨勢[3]。

如何對農田管理過程中的水分養分進行綜合調控，全面提升水肥利用效率成為當前研究的重點課題。目前，國內已進行了大量研究，如水肥耦合增產效應及其機理研究、水肥生產函數的研究等，在此基礎上提出了不同作物、不同土壤肥力、不同種植環境條件(設施和大田)、以及不同氣候條件下的優化灌溉施肥方案[4-8]，這些研究成果屬于充分供水灌溉情況，對農業生產發揮了重要指導作用。但對于普遍存在的有限供水灌溉，尚未見到相關報道。據此，本文在田間試驗和理論研究的基礎上，分析確定限水供水條件下作物優化灌溉施肥方案，實現農田灌溉施肥的動態的、優化的管理，對豐富和發展現代灌溉施肥理論，推動精準灌溉施肥技術的發展，具有重要的科學意義。

1 灌溉施肥制度優化模型

1.1 優化模型

以PS-123模型(作物生長模型)[9]作為冬小麥生長過程模擬的基礎，在此基礎上構建限水灌溉制度優化模型，該模型的目標函數為單位面積純收益，主要變量包括產量、灌水次數、施肥量，約束條件主要受灌水時間的限制，由此求得不同水文年限水灌溉的優化灌溉施肥制度，其中本次研究為了簡化計算，不考慮冬小麥生長期內灌水定額的變化，取固定值75 mm，施肥量主要針對追肥，追肥時間統一為第二次灌水之前，追肥后立即灌水(追肥為尿素)。本優化灌溉施肥制度，主要是對有限供水條件下不同可供水量時相應的灌水次數、灌水時間以及追肥數量所做的優化。其數學模型如下：

maxB=Pcy-PwmJ/η/1.5-Pff-C0

(1)

式中：B為單位面積的純收益，元/hm2；產量y用作物水肥生產函數計算，kg/hm2；J為返青到收獲期的灌水次數；m為灌水定額，取m=75 mm；η為灌溉水利用系數，這里取η=0.5；f為追肥量，kg/hm2；Pc、Pw、Pf分別為冬小麥產品價格、灌溉水價格、尿素價格；C0為除灌溉水外的其他農業投入，元/hm2， 取C0=3 150 元/hm2。

1.2 模型求解

采用序列無約束最小化技術中的內點法求解上述灌溉施肥制度優化模型，針對某一水文年，假定可供灌溉水量為75 mm，通過優化計算，可確定最優灌水時間、施肥量。依次假定可供灌溉水量為150、225、300 mm…，分別計算不同可供水量下的優化灌水時間、施肥量、產量和效益，直到效益開始減小為止。由此可確定該水文年不同可供水量條件下的優化灌溉施肥制度(對應效益最大的灌水次數、灌水時間及施肥量)以及在優化灌溉基礎上的產量、施肥量與供水量的關系。

2 計算實例

本研究利用山西省臨汾市灌溉試驗站開展的2011-2012、2012-2013、2013-2014三個年度的冬小麥田間試驗資料進行了作物生長模型參數的率定及檢驗，構建適合該地區的土壤水分、養分特征參數。同時以該站1954-2014年共計61 a的氣象資料為依據，分析研究限水灌溉下冬小麥的優化灌水時間及施肥量。

2.1 試驗設計概況

試驗站位于山西省臨汾市堯都區劉村鎮東宜村西(111°18′38″E，36°06′00″N)，海拔高程449 m。全年平均氣溫為12.1 ℃，降雨量為486.7 mm左右。試驗田土壤為中壤土，田間持水率為25.7%，地下水埋深在3 m以內，0～100 cm土壤剖面平均干體積質量為1.42 g/cm3。

2011-2012、2012-2013、2013-2014三個年度冬小麥各個處理都進行了土壤含水量、土壤硝態氮、以及葉面積、莖、葉和籽粒干物質量等的測試。試驗灌水方式為地面畦灌，單次灌水定額設定為75 mm，底肥為復合肥，追肥為尿素，追肥隨第二次灌水進行。

按照水分和肥料兩因素進行區組設計，灌水時間、中等施肥水平以及施肥時間均參照當地農戶的管理措施進行。灌水設置高水、中水以及不灌水3個水平，對應的灌溉水量為225、150和0 mm，施肥設置4個水平分別為高肥、中肥、低肥以及不追肥，底肥數量對應為900、600、300、600 kg/hm2，追肥量分別為337.5、225、112.5、0 kg/hm2。采用田間小區對比的方法進行小區試驗，小區規格為20 m×3.3 m，每個處理設3個重復小區，重復小區相鄰布置。

2.2 作物生長模型參數率定及檢驗

2.2.1 作物生長期土壤水分養分參數

本研究采用2011-2012年度臨汾市冬小麥田間試驗資料進行參數的率定，以2012-2013、2013-2014兩個年度試驗資料進行參數的檢驗。首先以模擬含水率與實測含水率誤差平方和最小為目標函數，經過多次迭代計算，求得作物系數和土壤水分運動參數最優值；在此基礎上，再以硝態氮模擬值與實測值誤差平方和最小為目標函數，求得土壤養分(氮素)運移和轉化參數最優值，土壤水分和氮素的模擬值與實測值比較見圖1、2。

圖1 土壤含水率模擬值與實測值比較圖

圖2 土壤硝態氮模擬值與實測值比較圖

2.2.2 作物生長過程參數

利用冬小麥零水不追肥處理、高水高肥處理、高水低肥處理的冬小麥莖干重、葉干重和籽粒干重及其地上部總干物重等試驗資料，依據參數反演法求得其水分虧缺敏感指數σ=3.0，養分虧缺敏感指數λ=3.0，光合產物轉化效率Yg=0.5。圖3給出了2012-2013年度高水低肥處理冬小麥其地上部總干物質重、籽粒干重、莖干重及葉干重模擬計算值和實測值隨時間的變化過程及其比較的散點圖。

由圖1～圖3可知，土壤水分、氮素、干物質重、籽粒干重、莖干重及葉干重模擬計算值與實測值都比較一致，表明所確定的參數值合理可靠，可用于模擬水分和養分虧缺對冬小麥生長過程和產量的影響。

圖3 2012-2013高水低肥處理模擬值與實測值比較圖

2.3 優化灌溉施肥制度求解結果

本研究以山西省臨汾市灌溉試驗站作物生長模型為依據，以山西省臨汾市1951-2014年氣象資料為背景，將冬小麥全生育期(10月1日到第二年6月10日，共計253天)的降水量做頻率分析，確定了P=95%、75%、50%、25%和5%五個水文年，以此為依據求得了5種水文年不同灌溉供水量條件下的優化灌水時間、經濟用水灌溉制度(定量灌水定額下對應效益最大的灌水次數及灌水時間)以及追肥量，見表2，以及在優化灌溉基礎上的產量和效益與供水量和追肥量的關系，見圖4、5。表2灌水時間為從播種日算起的天數，降水量、可供灌水量為從灌水日算起到作物收獲期間的合計值，施肥量包括了底肥和追肥，底肥為定值600 kg/hm2，其余為追肥量。

由表1可見，P=95%、75%、50%、25%和5%典型年的最優灌水次數分別為4次、4次、3次、3次、2次；在同一水文年條件下，隨著灌溉供水量的增加，產量在增加，但效益卻先增大后減小，存在最優值；在同一灌溉供水量條件水平下，干旱程度越大(頻率越大)，總的趨勢是產量越小，但是，也有例外，如灌溉供水量大于150 mm時，25%和50%年份的產量小于75%年份產量，主要原因是25%、50%年份的日照時數偏小，僅有1 356、1 583 h，較75%的日照時數小300～400 h，干旱脅迫的影響要小于光照不足脅迫的影響。

表1 五種水文年不同灌溉供水量條件下的優化灌溉制度

續表1 五種水文年不同灌溉供水量條件下的優化灌溉制度

對不同灌溉供水條件下的追肥量進行研究(見表2)，表明隨著灌溉供水量的增加，追肥量也在增加，且在同一灌溉供水量水平時，隨著干旱程度的增大，追肥量基本上呈現了先增加后減少的趨勢，說明水分不足直接影響了小麥對養分的吸收；同一可供灌溉水量條件下，不同典型年的追肥數量是不同的，如可供灌溉水量為150 mm時，追肥的數量的變化范圍在180～360 kg/hm2，平均值為4 050 kg/hm2，不同典型年不同可供灌溉水量的追肥量平均值為248 kg/hm2；不同典型年某一可供灌溉水量條件下，隨著追肥量的增加，產量在逐漸增大，但增加幅度逐漸減少，效益先增大后減小，存在最優值(見圖4～圖5，以P=25%年份可供灌溉水量為75 mm為例)。

表2 五種水文年不同灌溉供水量條件下的追肥量

3 結 論

本文以山西省臨汾市灌溉試驗站3年冬小麥全生育期實測田間試驗資料為依據，以PS-123模型作為冬小麥生長過程模擬的基礎，分析確定限水供水條件下作物優化灌溉施肥方案，主要結論如下：

(1)充分利用試驗測試的田間土壤含水率、土壤養分資料，綜合灌水和施肥的耦合增產效應，率定和檢驗了PS-123模型的參數。研究結果表明，土壤水分和土壤氮素模擬結果都比較滿意，相關系數R2分別達到0.389 3、0.589 8。

(2)利用PS-123模型，以水肥耦合增產理論為依據，通過優化技術將水肥生產函數、土壤水分養分動態變化規律，與田間土壤水分養分監測技術有機地結合在一起，形成動態、優化的灌溉施肥決策方法。

(3)在同一水文年條件下，隨著灌溉供水量的增加，追肥量在增加，產量也在增加，但效益卻先增大后減小，存在最優值；在同一灌溉供水量條件水平下，干旱程度越大(頻率越大)，追肥量基本上呈現了先增加后減少的趨勢，說明水分不足直接影響了小麥對養分的吸收，且產量也基本上呈現減少的趨勢。

(4)不同典型年某一可供灌溉水量條件下，隨著追肥量的增加，產量在逐漸增大，但增加幅度逐漸減少，效益先增大后減小，存在最優值。

[1] 魏成熙. 喀斯特山區冬小麥施肥量與產量的關系[J]. 干旱地區農業研究, 2006,24(3): 60-63.

[2] 劉增進, 李寶萍, 李遠華, 等. 冬小麥水分利用效率與最優灌溉制度的研究[J]. 農業工程學報, 2004,20(4):58-63.

[3] 王仰仁, 杜秀文, 張紹強. 限量供水條件下灌溉預報增產效益分析[J].中國農村水利水電，2016,(12)：1-7.

[4] 汪德水. 旱地農田肥水協同效應與耦合模式[M].北京: 氣象出版社, 1999.

[5] 徐振劍, 華 珞, 蔡典雄,等. 農田水肥關系研究現狀[J].首都師范大學學報(自然科學版), 2007,28(1):83-88.

[6] 肖自添, 蔣衛杰, 余宏軍. 作物水肥耦合效應研究進展[J].作物雜志, 2007,(6):18-22.

[7] 謝 偉, 黃 璜, 沈建凱.植物水肥耦合研究進展[J].作物研究, 2007,21(5):541-546.

[8] 翟丙年,李生秀.冬小麥產量的水肥耦合模型[J].中國工程科學, 2005,4(9):69-74.

[9] 王仰仁,孫小平. 山西農業節水理論與作物高效用水模式[M].北京: 中國科學技術出版社,2003.