摘 要：利用田間實測數(shù)據對ＡｑｕａＣｒｏｐ 模型進行校驗，獲取了模型的關鍵參數(shù)，并量化不同施氮量水平與肥脅迫相關的壓力系數(shù)；以此為基礎，提出了半定量方法的參數(shù)化方案，實現(xiàn)了華北冬小麥施氮量的優(yōu)化管理。結果表明：①模型校驗過程中冠層覆蓋度的決定系數(shù)Ｒ２大于０．８７６、ＮＲＭＳＥ 小于１５．５５４％；土壤含水率模擬值與實測值的Ｒ２ 大于０．８４０、ＮＲＭＳＥ 小于９．９６０％；所有處理下產量模擬值與實測值的Ｒ２大于０．９６４。說明經校驗的模型能較好地模擬不同施氮處理下作物生長及產量產出過程，可利用該模型進行施氮量優(yōu)化管理研究。②利用半定量方法實現(xiàn)了氮肥的多情景模擬，根據模擬結果建立了施氮量與產量之間的關系，當施氮量達到３５５．３０ ｋｇ／ ｈｍ２后，繼續(xù)增加施氮量，冬小麥的產量出現(xiàn)了降低的趨勢。

關鍵詞：冬小麥；施氮量；氮肥利用效率；ＡｑｕａＣｒｏｐ 模型；半定量方法

中圖分類號：Ｓ５１３； Ｓ５１２．１ 文獻標志碼： Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２３．０４．０１９

引用格式：牟彥文，彭致功，仵峰，等．基于半定量方法的華北冬小麥施氮量優(yōu)化管理研究［Ｊ］．人民黃河，２０２３，４５（４）：１０９－１１４，１２５．

華北平原是中國重要的糧食生產基地，其中冬小麥播種面積和產量均居中國首位，在國家糧食安全中具有重要作用［１］ 。氮肥是作物生長所需氮素的重要來源，對葉片光合能力和產量形成過程具有顯著的調控作用［２］ 。近年來，我國過量施用氮肥的現(xiàn)象越來越普遍，嚴重影響了氮肥利用效率，并導致作物出現(xiàn)不同程度的減產，同時也引發(fā)嚴重的環(huán)境安全問題［３－５］ 。因此，探究不同施氮量對作物生長的影響，確定適宜的施氮量，對提高氮肥利用效率和農業(yè)的可持續(xù)發(fā)展具有重要作用。

由于田間試驗耗時費力，且受外界因素的影響較大，因此尋求模擬精度高、適用范圍廣的作物模型成為當下研究者們關注的重點［６－７］ 。目前，較為常見的作物模型有ＡＰＳＩＭ［８］ 、ＷＯＦＯＳＴ［９］ 、ＤＳＳＡＴ － ＣＥＲＥＳ［１０］等，然而這些模型的研究機理和輸入參數(shù)較復雜，造成了模型在實際應用中受限。與其他模型相比，ＡｑｕａＣｒｏｐ 模型具有輸入參數(shù)少、界面直觀清晰、模擬精度高等優(yōu)點［１１］ 。國外的一些研究表明，ＡｑｕａＣｒｏｐ模型能夠較好地模擬大麥［１２］ 、小麥［１３］ 、玉米［１４］ 等作物產量對缺氮脅迫的響應，但受當前試驗處理的限制，難以尋求精度較高的施氮制度［１５］ 。國內對該模型與田間施氮量相關的研究較少［１６］ ，對于指導作物田間施氮量管理研究的借鑒意義不大。

本文以華北地區(qū)冬小麥為研究對象，結合田間試驗的分析結果，利用ＡｑｕａＣｒｏｐ 模型的半定量方法進行多情景模擬，以尋求華北平原保障冬小麥高產高效的優(yōu)化施氮制度，為提高氮素利用效率和實現(xiàn)模型的氮肥優(yōu)化管理提供理論依據。

１ 研究區(qū)概況和試驗設計

１．１ 研究區(qū)概況

試驗地點位于中國水利水電科學研究院大興節(jié)水灌溉試驗研究基地（北緯３９°３７．２５′、東經１１６°２５．５１′，海拔３０ ｍ），屬半濕潤大陸性季風氣候區(qū)，豐富的光熱條件適宜冬小麥、玉米等作物生長。土壤類型為沙壤土，試驗期間研究區(qū)年均降水量約５４０ ｍｍ，多年平均氣溫１２．１ ℃，年均水面蒸發(fā)量在１ ８００ ｍｍ 以上，全年日照時數(shù)約為２ ６００ ｈ。

１．２ 試驗設計

２０１９—２０２０ 年冬小麥的播種和收獲日期分別為２０１９ 年１０ 月１１ 日和２０２０ 年６ 月１３ 日，供試品種為京生麥１ 號，生育期內降水量為１５７ ｍｍ；２０２０—２０２１ 年冬小麥的播種和收獲日期分別為２０２０ 年１０ 月２３ 日和２０２１ 年６ 月１５ 日，供試品種為京冬１７ 號，生育期內降水量為９２ ｍｍ。本試驗設計５ 種施氮處理，施氮總量分別為０ ｋｇ／ ｈｍ２ （Ｎ１ 處理）、１５０ ｋｇ／ ｈｍ２ （Ｎ２ 處理）、２２５ ｋｇ／ ｈｍ２（Ｎ３ 處理）、３００ ｋｇ／ ｈｍ２（Ｎ４ 處理）、４２０ ｋｇ／ ｈｍ２（Ｎ５ 處理）；在作物生育期施肥分２ 次完成，每次施肥量相同，其中播前施復合肥（含Ｎ 量１５％、Ｐ２Ｏ５量１５％、Ｋ２Ｏ量１５％）作為底肥施入，在拔節(jié)期施入尿素（含Ｎ 量５２％）。各處理灌溉制度保持一致，考慮播前、播種—越冬、返青—拔節(jié)、拔節(jié)—抽穗４ 個灌水階段，每次灌水定額均為６０ ｍｍ。各處理試驗區(qū)面積均為５６ ｍ２，行間距為２０ ｃｍ，為避免各試驗區(qū)間土壤水分側滲，設有寬１．０ ｍ的保護行，其他管理措施如播種、耕作及鋤草等與當?shù)剞r民習慣保持一致。

１．３ 數(shù)據獲取

氣象數(shù)據來自節(jié)水灌溉研究基地的氣象站，各氣象資料見圖１。采用ＴＲＩＭＥＲ－Ｔ３／ ＩＰＨ 型土壤水分測定儀測定０～１００ ｃｍ 深度土層的土壤含水率，每２０ ｃｍ深度測定１ 個數(shù)據，觀測間隔時間為５～７ ｄ，灌水前后和降雨后分別加測１ 次。冬小麥返青后，每１０ ｄ 測１次生物量，各處理小區(qū)選取具有代表性的３ 株，從莖部剪下、烘干后稱重，并根據種植密度和面積折算成單位面積地面上部生物量。在小麥成熟時，各試驗小區(qū)選?。?個有代表性的１ ｍ２樣方進行考種測產，經自然風干后，測定籽粒質量（下文簡稱產量）。采用凱氏定氮儀測定各葉片以及植株含氮量［１７］ 。各處理小區(qū)選?。?株長勢均勻的植株，利用直尺測量各葉片的長和寬，利用折算系數(shù)法推算冬小麥葉面積指數(shù)。

２ 結果與分析

２．１ 施氮量對冬小麥產量和氮素吸收與利用效率的影響

由表３ 可知，在試驗中Ｎ４ 處理的產量達到最大值，２０１９—２０２０ 年度、２０２０—２０２１ 年度分別為６ １５２． ３３、５ ３９０．１０ ｋｇ／ ｈｍ２，所以在０ ～３００ ｋｇ／ ｈｍ２ 的施氮量范圍內，增施氮肥能夠顯著提高冬小麥的產量。但是施氮量超過３００ ｋｇ／ ｈｍ２后，對增產反而起反作用，與Ｎ４ 處理相比，２０１９—２０２０ 年度、２０２０—２０２１ 年度在Ｎ５ 處理下的產量分別降低了２．２８％、６．９２％。氮素積累總量都隨施氮量的增加而增加，分３ 個層次達到顯著水平，但Ｎ４ 和Ｎ５ 處理之間差異不顯著。各處理的氮素干物質生產效率隨施氮量的增加呈遞減趨勢，不施氮或者施氮量低的處理與較高施氮量處理之間差異顯著（顯著性水平ｐ＜０．０５）。Ｎ３、Ｎ４ 和Ｎ５ 處理的氮素生產效率差異不顯著，與處理Ｎ１ 的差異顯著（ｐ＜０．０５）。各處理氮素農藝效率隨施氮量的增加呈現(xiàn)先上升再下降的趨勢，其中Ｎ４ 處理最高、Ｎ２ 處理最低。４ 種施氮處理的氮肥偏生產力都隨施氮量的增加呈現(xiàn)遞減趨勢， 各處理間差異顯著（ｐ＜０．０５）。

２．２ 模型校驗及應用

２．２．１ 冬小麥冠層覆蓋度與土壤含水率的校驗

模型率定情況（見圖２ 和表４）中所有處理冠層覆蓋度模擬值與實測值的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０．９１８、９．９３５％、０．８０６、０．９４０，各處理評價指標的變化范圍為０．９４３ ＜Ｒ２ ＜ ０．９７０、６．７８３％ ＜ＮＲＭＳＥ ＜ １１．３８９％、０．７１２＜ＥＮＳ ＜０．８７９、０．９００＜ｄ＜０．９６０；所有處理土壤含水率模擬值與實測值的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０．７８８、１４．７９１％、０．５４６、０．８５９，各處理評價指標的變化范圍為０．８４３＜Ｒ２ ＜０．９６５、７．４６８％＜ＮＲＭＳＥ＜９．１９１％、０．６７３＜ＥＮＳ ＜０．８７９、０．９３０＜ｄ＜０．９６９。模型率定情況表明模型模擬冠層覆蓋度和土壤含水率的誤差較小。模型驗證情況（見圖３ 和表４）中所有處理冠層覆蓋度模擬值與實測值的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０．８７４、８．１４５％、０．７１２、０．９４４，且各評價指標的范圍分別為０． ８７６ ＜ Ｒ２ ＜０．９５２、９． ６８３％ ＜ ＮＲＭＳＥ ＜ １５． ５５４％、０． ４３５ ＜ ＥＮＳ ＜０．６４５、０．８３９＜ｄ＜０．８７９；所有處理土壤含水率模擬值與實測值的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別０．８４８、８．７５４％、０．６２５、０．９３１，各處理評價指標的變化范圍為０．８４０＜ Ｒ２ ＜０．９４２、７． ９０１％ ＜ＮＲＭＳＥ ＜ ９． ９６０％、０． ５４２ ＜ＥＮＳ ＜０．８１７、０．９１７＜ｄ＜０．９５８?？梢姡ｒ灪蟮哪Ｐ湍茌^好地模擬生育期內冠層覆蓋度與土壤含水率的動態(tài)變化過程。

２．２．２ 冬小麥最終生物量與產量的校驗

模型率定期冬小麥最終生物量和產量的校驗結果見圖４，所有處理下最終生物量的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０．９９０、４．０９７％、０．９０１、０．９７９，產量的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０． ９７４、４． ９５６％、０． ９２７、０．９７８；模型驗證期冬小麥最終生物量和產量的校驗結果見圖５，所有處理下最終生物量的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０．９６４、４．６６０％、０．８７７、０．９７５，產量的Ｒ２、ＮＲＭＳＥ、ＥＮＳ、ｄ 分別為０． ９６８、７． ７７３％、０． ７９４、０．９４８。模型校驗結果表明ＡｑｕａＣｒｏｐ 模型能夠較好地模擬冬小麥最終生物量與產量。各處理的相對誤差見表５，最終生物量相對誤差的最小值和最大值分別為０．８０２％、９．１６７％，所有處理平均相對誤差為４．２２７％；產量相對誤差的最小值與最大值分別為０．８５０％、１１． ９２３％， 所有處理平均相對誤差為５．９７２％。驗證期相對誤差偏大，主要原因是模型校驗過程中的肥脅迫參數(shù)直接由各處理施氮量確定，沒有兼顧土壤背景值差異，造成了模擬誤差的偏差，但模擬誤差總體在可接受范圍內。

２．２．３ 不同肥脅迫水平情景模擬

依據模型校驗后各處理壓力系數(shù)與施氮量的關系，并結合華北地區(qū)施氮量實際情況，設置了１５ 種施肥情景方案，施氮量范圍為０～４２０ ｋｇ／ ｈｍ２，梯度為３０ ｋｇ／ｈｍ２，對應的序號依次為Ｆ１、Ｆ２、…、Ｆ１５。在保證各處理不受水分脅迫下，１５ 種情景方案下模擬獲得的小麥產量與施氮量呈拋物線關系（見圖６），冬小麥產量隨施氮量的增加而增大，但當施氮量超過某一范圍時，冬小麥產量的增長趨勢減緩甚至出現(xiàn)了降低，Ｆ１３、Ｆ１４、Ｆ１５ 處理相較于Ｆ１２ 處理分別下降了０． ６７１％、４．４７６％、７．５９３％。根據ｙ ＝－０．０２１ ５ｘ２ ＋１５．２７８ｘ＋３ ３７６．４（Ｒ２ ＝０．９８０ ２）計算得出，當施氮量為３５５．３０ ｋｇ／ ｈｍ２時，冬小麥的產量達最大值，為６ ０９０．５５ ｋｇ／ ｈｍ２。

２．３ 討 論

作物模型已成為實現(xiàn)農業(yè)精準施氮管理的有效途徑［１８－１９］ 。目前已有學者將ＡｑｕａＣｒｏｐ 模型與其他模型進行了對比，發(fā)現(xiàn)該模型沒有考慮作物的生理及養(yǎng)分循環(huán)與平衡，而是通過半定量方法來描述作物受到的肥脅迫程度，模擬誤差仍在可接受的范圍內［２０］ 。本研究采用田間試驗與模型相結合的方法，通過ＡｑｕａＣｒｏｐ模型進行了不同施氮水平的多情景模擬，基于模擬產量與施氮量的相關曲線得到華北平原冬小麥適宜施氮量，對于指導區(qū)域施氮管理提供了理論依據。結果表明，冬小麥的產量隨施氮量的增加表現(xiàn)為先增加后減少的趨勢，當施氮量達到３５５．３０ ｋｇ／ ｈｍ２時，冬小麥的產量達到了最大值，其氮肥吸收與利用效率均在合理的范圍內，與李瑞奇等［３］ 、郭曾輝等［２１］ 、張定一等［２２］ 的研究結果基本一致。

３ 結 論

（１）ＡｑｕａＣｒｏｐ 模型能夠較好地模擬華北冬小麥生育期內冠層覆蓋度、土壤含水率、地上部生物量與產量，且模型校驗評價指標均在合理的范圍內，表明可以利用模型的半定量方法進行施氮量管理研究。

（２）基于施氮量與肥脅迫參數(shù)之間的關系，來量化任一施氮量水平與肥脅迫相關的參數(shù)值，破解了采用半定量方法進行施氮量優(yōu)化管理時模型參數(shù)值的確定難題。

（３）基于模型的半定量方法進行了多情景模擬，建立了施氮量與產量之間的二次拋物線函數(shù)關系，得到施氮量在０～３５５．３０ ｋｇ／ ｈｍ２ 時冬小麥的產量隨施氮量的增加而增大，超過這一區(qū)間，冬小麥的產量出現(xiàn)了下降的趨勢，此結果與田間試驗結果一致，過量施用氮肥會造成氮肥利用效率降低。

參考文獻：

［１］ 孫爽，楊曉光，張鎮(zhèn)濤，等．華北平原不同等級干旱對冬小麥產量的影響［Ｊ］．農業(yè)工程學報，２０２１，３７（１４）：６９－７８．

［２］ 劉婷婷，張興燕，褚鵬飛，等．微噴灌條件下減量施氮對冬小麥光合特性和產量的影響［Ｊ］．山東農業(yè)科學，２０２２，５４（７）：４６－５２．

［３］ 李瑞奇，李雁鳴，何建興，等．施氮量對冬小麥氮素利用和產量的影響［Ｊ］．麥類作物學報，２０１１，３１（２）：２７０－２７５．

［４］ 趙亞南，宿敏敏，呂陽，等．減量施肥下小麥產量、肥料利用率和土壤養(yǎng)分平衡［Ｊ］．植物營養(yǎng)與肥料學報，２０１７，２３（４）：８６４－８７３．

［５］ 蔡媛媛，王瑞琪，王麗麗，等．華北平原不同施氮量與施肥模式對作物產量與氮肥利用率的影響［Ｊ］．農業(yè)資源與環(huán)境學報，２０２０，３７（４）：５０３－５１０．

［６］ 楊靖民，楊靖一，姜旭，等．作物模型研究進展［Ｊ］．吉林農業(yè)大學學報，２０１２，３４（５）：５５３－５６１．

［７］ 朱秀芳，李宜展，潘耀忠，等．ＡｑｕａＣｒｏｐ 作物模型研究和應用進展［Ｊ］．中國農學通報，２０１４，３０（ ８）：２７０－２７８．

［８］ 尹嘉德，侯慧芝，張緒成，等．基于ＡＰＳＩＭ 模型分析不同降水年型下施氮深度對旱地小麥的產量效應［Ｊ］．水土保持學報，２０２２，３６（１）：２４７－２５４．

［９］ 秦雅倩，孫嘉玥，王巖，等．基于ＷＯＦＯＳＴ 模型的河南省冬小麥產量差的時空特征分析［Ｊ］．江蘇農業(yè)科學，２０２２，５０（５）：１９１－１９８．

［１０］ 李毅，張思遠，劉慶祝，等．基于ＤＳＳＡＴ－ＣＥＲＥＳ－Ｗｈｅａｔ的黃土高原西部春小麥干旱影響研究［Ｊ］．農業(yè)機械學報，２０２２，５３（６）：３３８－３４８．

［１１］ 孫仕軍，張琳琳，陳志君，等．ＡｑｕａＣｒｏｐ 作物模型應用研究進展［Ｊ］．中國農業(yè)科學，２０１７，５０（１７）：３２８６－３２９９．

［１２］ ＡＢＩ ＳＡＡＢ Ｍ Ｔ，ＴＯＤＯＲＯＶＩＣ Ｍ，ＡＬＢＲＩＺＩＯ Ｒ．ＣｏｍｐａｒｉｎｇＡｑｕａＣｒｏｐ ａｎｄ ＣｒｏｐＳｙｓｔ Ｍｏｄｅｌｓ ｉｎ Ｓｉｍｕｌａｔｉｎｇ Ｂａｒｌｅｙ Ｇｒｏｗｔｈａｎｄ Ｙｉｅｌｄ Ｕｎｄｅｒ Ｄｉｆｆｅｒｅｎｔ Ｗａｔｅｒ ａｎｄ Ｎｉｔｒｏｇｅｎ Ｒｅｇｉｍｅｓ，Ｄｏｅｓ Ｃａｌｉｂｒａｔｉｏｎ Ｙｅａｒ Ｉｎｆｌｕｅｎｃｅ ｔｈｅ Ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ＣｒｏｐＧｒｏｗｔｈ Ｍｏｄｅｌｓ？ ［Ｊ］．Ａｇｒｉｃｕｌｔｕｒａｌ Ｗａｔｅｒ Ｍａｎａｇｅｍｅｎｔ，２０１５，１４７：２１－３３．

［１３］ ＡＢＥＤＩＮＰＯＵＲ Ｍ．Ｔｈｅ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ＤＳＳＡＴ?ＣＥＲＥＳ ａｎｄＡｑｕａＣｒｏｐ Ｍｏｄｅｌｓ ｆｏｒ Ｗｈｅａｔ Ｕｎｄｅｒ?Ｎｉｔｒｏｇｅｎ Ｉｎｔｅｒａｃｔｉｏｎｓ［Ｊ］．Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ ｉｎ Ｓｏｉｌ Ｓｃｉｅｎｃｅ ａｎｄ Ｐｌａｎｔ Ａｎａｌｙｓｉｓ，２０２１，５２（１７）：２００２－２０１７．

［１４］ ＲＡＨＩＭＩＫＨＯＯＢ Ｈ，ＳＯＨＲＡＢＩ Ｔ，ＤＥＬＳＨＡＤ Ｍ．ＳｉｍｕｌａｔｉｎｇＣｒｏｐ Ｒｅｓｐｏｎｓｅ ｔｏ Ｎｉｔｒｏｇｅｎ?Ｄｅｆｉｃｉｅｎｃｙ Ｓｔｒｅｓｓ Ｕｓｉｎｇ ｔｈｅ Ｃｒｉｔ?ｉｃａｌ Ｎｉｔｒｏｇｅｎ Ｃｏｎｃｅｎｔｒａｔｉｏｎ Ｃｏｎｃｅｐｔ ａｎｄ ｔｈｅ ＡｑｕａＣｒｏｐＳｅｍｉ?Ｑｕａｎｔｉｔａｔｉｖｅ Ａｐｐｒｏａｃｈ ［ Ｊ］． Ｓｃｉｅｎｔｉａ Ｈｏｒｔｉｃｕｌｔｕｒａｅ，２０２１，２８５：１１０１９４．

［１５］ ＡＢＥＤＩＮＰＯＵＲ Ｍ，ＳＡＲＡＮＧＩ Ａ，ＲＡＪＰＵＴ Ｔ Ｂ Ｓ，ｅｔ ａｌ．Ｐｅｒ?ｆｏｒｍａｎｃｅ Ｅｖａｌｕａｔｉｏｎ ｏｆ Ａｑｕａ Ｃｒｏｐ Ｍｏｄｅｌ ｆｏｒ Ｍａｉｚｅ Ｃｒｏｐ ｉｎａ Ｓｅｍｉ?Ａｒｉｄ Ｅｎｖｉｒｏｎｍｅｎｔ［Ｊ］．Ａｇｒｉｃｕｌｔｕｒａｌ Ｗａｔｅｒ Ｍａｎａｇｅ?ｍｅｎｔ，２０１２，１１０：５５－６６．

［１６］ 張濤，孫偉，張鋒偉，等．旱地全膜雙壟溝玉米生產的ＡｑｕａＣｒｏｐ 模型模擬及管理措施優(yōu)化［Ｊ］．應用生態(tài)學報，２０１７，２８（３）：９１８－９２６．

［１７］ 溫云杰，李桂花，黃金莉，等．連續(xù)流動分析儀與自動凱氏定氮儀測定小麥秸稈全氮含量之比較［Ｊ］．中國土壤與肥料，２０１５（６）：１４６－１５１．

［１８］ 代娜，史文嬌，石曉麗．華北平原多目標約束下冬小麥最優(yōu)施氮量［Ｊ］．中國生態(tài)農業(yè)學報（中英文），２０２１，２９（９）：１５１２－１５２３．

［１９］ 冷旭，李仙岳，郭宇，等．ＤＮＤＣ 模型和ＴＯＰＳＩＳ 法聯(lián)合確定生物降解地膜覆蓋農田合理施氮量［Ｊ］．水土保持學報，２０２２，３６（４）：３２５－３３２，３４１．

［２０］ 戴嘉璐，李瑞平，李聰聰，等．河套灌區(qū)不同水肥處理對玉米生長影響的ＡｑｕａＣｒｏｐ 模型模擬［Ｊ］．水土保持學報，２０２１，３５（３）：３１２－３１９．

［２１］ 郭曾輝，雒文鶴，劉朋召，等．微噴灌條件下不同水氮處理對冬小麥水氮利用和產量的影響［Ｊ］．麥類作物學報，２０２１，４１（１０）：１２５６－１２６５．

［２２］ 張定一，黨建友，王姣愛，等．施氮量對不同品質類型小麥產量、品質和旗葉光合作用的調節(jié)效應［Ｊ］．植物營養(yǎng)與肥料學報，２００７，１３（４）：５３５－５４２．

【 責任編輯 簡 群】