焦 貞，張衍福，劉 凱，于 瀟，徐征和，龐桂斌，徐 晶

（1.濟(jì)南大學(xué)，濟(jì)南250022；2.山東省水利綜合事業(yè)服務(wù)中心，濟(jì)南250013）

我國北方水資源短缺，越來越嚴(yán)重地制約著農(nóng)業(yè)生產(chǎn)。山東省人均水資源量低于500 m3，農(nóng)業(yè)用水占總用水量七成以上，使得水資源供需矛盾更加尖銳。利用有限的灌溉水量提高灌水利用效率、優(yōu)化灌溉制度是實(shí)現(xiàn)節(jié)水灌溉的有效途徑[1,2]。

為了提高農(nóng)業(yè)水資源利用效率，眾多學(xué)者對(duì)灌溉制度優(yōu)化策略進(jìn)行了大量研究。馬健琴[3]采用改進(jìn)層次分析的模糊綜合評(píng)價(jià)法優(yōu)選了玉米的灌水方案。張志宇[4]建立冬小麥和夏玉米的水分生產(chǎn)函數(shù)模型，并通過粒子群優(yōu)化算法計(jì)算模型中的敏感指數(shù)，在此基礎(chǔ)上實(shí)現(xiàn)灌溉制度的多目標(biāo)優(yōu)化模型建立。楊習(xí)清[5]通過多作物、多約束的非線性灌溉優(yōu)化模型，實(shí)現(xiàn)了非充分灌溉條件下灌溉水量的優(yōu)化配置。SWAP模型在指導(dǎo)農(nóng)田作物灌溉方面的應(yīng)用更為廣泛。繳錫云[6]建立了田間水層深度模型，劉虎[7]模擬了春小麥全生育期和生育階段的作物騰發(fā)量變化規(guī)律與土壤水分狀況，均發(fā)現(xiàn)模擬值與實(shí)測(cè)值的吻合度較高。袁成福[8-10]利用SWAP 模型在干旱區(qū)進(jìn)行了小麥和玉米的模擬灌溉，并預(yù)測(cè)作物指標(biāo)的長期變化，擬定灌溉方案。Jiang[11]、馮紹元[12]等利用SWAP 模型模擬土壤含水率、水鹽運(yùn)移和相對(duì)產(chǎn)量，提出最適宜的灌溉定額和礦化度區(qū)間。這都為農(nóng)業(yè)灌溉的發(fā)展提供了理論依據(jù)。

于成[2]對(duì)山東省主要作物的水足跡研究發(fā)現(xiàn)冬小麥和夏玉米的生產(chǎn)水足跡較大，說明2種作物水資源利用效率還有較大潛力?；概_(tái)縣人均水資源量僅為354 m3，且時(shí)空分布不均[13]。因此，本研究將利用SWAP模型，致力于冬小麥和夏玉米的灌溉制度優(yōu)化，這對(duì)克服我國農(nóng)業(yè)用水效率低下的缺陷，達(dá)到節(jié)水、增產(chǎn)、高效的目的具有十分現(xiàn)實(shí)的意義。

1 研究區(qū)概況

桓臺(tái)縣隸屬于山東省淄博市，位于山東中北部，全縣面積509 km2，屬暖溫帶大陸性季風(fēng)氣候，雨熱同期，年平均氣溫12.5 ℃，降雨多集中于夏季七八月份。土壤類型包括褐土、潮土和砂姜黑土，其中褐土為主要土壤類型。試驗(yàn)站位于桓臺(tái)縣新城鎮(zhèn)逯家村南（117°59′36.5″N，36°58′30.6″E），占地面積2.0 hm2，站內(nèi)土壤類型以潮土為主，主要作物種類為冬小麥和夏玉米，一年兩作，復(fù)種指數(shù)約1.9。

2 試驗(yàn)設(shè)計(jì)與研究方法

2.1 試驗(yàn)設(shè)計(jì)

本研究于2018年6月-2019年9月在山東省桓臺(tái)縣試驗(yàn)站進(jìn)行。試驗(yàn)在田間共設(shè)2 排測(cè)坑，每排18 個(gè)，共36 個(gè)。測(cè)坑?xùn)|部設(shè)自動(dòng)氣象站，觀測(cè)降水、氣溫、地溫、相對(duì)濕度和日照時(shí)數(shù)、風(fēng)速和大氣壓強(qiáng)等氣象資料，西部放置防雨棚。小麥供試品種為“濟(jì)麥17”，玉米供試品種為“鄭黃糯2 號(hào)”。冬小麥和夏玉米選擇灌3 水的灌溉處理。冬小麥種植過程中，測(cè)坑上方使用防雨棚消除自然降雨影響，設(shè)置4種不同灌溉定額的處理（G3-1、G3-2、G3-3、G3-4），均在苗期、抽穗和灌漿期進(jìn)行灌溉。夏玉米生長過程中對(duì)降雨的依賴較大，因此未進(jìn)行擋雨，設(shè)未灌溉的對(duì)照處理（CK）和3 種不同灌溉定額的處理（M3-1、M3-2、M3-3），均在拔節(jié)、抽雄和灌漿期進(jìn)行灌溉。作物的具體灌水方案見表1和表2。

表1 2018-2019年冬小麥試驗(yàn)方案Tab.1 Experimental scheme of winter wheat from 2018 to 2019

表2 2019年夏玉米試驗(yàn)方案Tab.2 Experimental scheme for summer maize in 2019

在作物不同生育期隨機(jī)選取采樣點(diǎn)，利用土鉆采集土樣，采樣深度1 m，每20 cm 采集一次土樣，烘干法測(cè)定土壤含水率。選取有代表性的測(cè)點(diǎn)，挖掘土壤剖面，環(huán)刀取土，測(cè)定土壤干容重；利用環(huán)刀原狀土吸水實(shí)驗(yàn)測(cè)定土壤田間持水量。土壤顆粒組成采用比重法進(jìn)行測(cè)量，土壤物理性質(zhì)見表3。

表3 土壤物理性質(zhì)Tab.3 Physical and chemical properties of soil

人工觀測(cè)收集不同生育期內(nèi)小麥和玉米的株高、根長和葉面積指數(shù)。采用LI-6200葉面積儀測(cè)定冬小麥和夏玉米的單株葉面積，每次選擇3株，取均值，并計(jì)算該時(shí)期的葉面積指數(shù)。利用卷尺測(cè)定不同生長期株高。定期（苗期－拔節(jié)、拔節(jié)－抽雄、抽雄－灌漿）采用直徑10 cm 的根鉆取冬小麥和夏玉米根系，取樣后立即沖洗、挑根、測(cè)定根長。作物收獲后統(tǒng)計(jì)穗數(shù)、穗粒數(shù)、千粒重等產(chǎn)量構(gòu)成因子。

2.2 研究方法

SWAP是一個(gè)基于物理的、詳細(xì)的農(nóng)業(yè)水文模型，可模擬土壤-大氣-植物環(huán)境系統(tǒng)中整個(gè)生長季節(jié)的垂直土壤水流、溶質(zhì)和熱傳遞過程[14]。該模型的上下邊界分別位于植物冠層之上和地下水飽和層之上，綜合考慮了自然大氣環(huán)境和地下水系統(tǒng)的動(dòng)態(tài)變化，在國內(nèi)外干旱半干旱地區(qū)得到廣泛應(yīng)用。

SWAP 模型簡單模型下作物生長的函數(shù)需要Doorenbos 和Kassam(1979年)提出的作物各生長階段的產(chǎn)量響應(yīng)因子作為發(fā)育階段的函數(shù)，以獲得實(shí)際產(chǎn)量。計(jì)算如下：

式中：Ky,k為生育期k的產(chǎn)量響應(yīng)因子；Tp,k為潛在蒸騰量，cm；Ta,k為k期實(shí)際蒸騰量，cm；指數(shù)k為生育期；n為確定的作物生育期數(shù)；Ea,k和Ep,k分別為各生育階段的實(shí)際產(chǎn)量和潛在產(chǎn)量。

全生育期相對(duì)產(chǎn)量按以下公式計(jì)算：式中：Ya為全生育期累積實(shí)際產(chǎn)量，kg/hm2；Yp為全生育期累積潛在產(chǎn)量，kg/hm2。

SWAP模型需要輸入的資料包括氣象資料、灌溉資料、作物生長數(shù)據(jù)、土壤參數(shù)資料、初始和邊界條件等。氣象資料（降水、氣溫、地溫、相對(duì)濕度、日照時(shí)數(shù)、風(fēng)速和大氣壓強(qiáng)等）由2018-2019年試驗(yàn)站自動(dòng)氣象站記錄所得；灌溉資料和土壤參數(shù)資料均采用田間實(shí)測(cè)數(shù)據(jù)；作物生長數(shù)據(jù)的線性函數(shù)用插值法構(gòu)造，出苗、開花、成熟時(shí)分別為0、1、2 表示，其余生長階段進(jìn)行線性插值，主要輸入數(shù)據(jù)有株高、根長、葉面積指數(shù)等，與作物生長階段（DVS）構(gòu)成相對(duì)應(yīng)的函數(shù)來模擬作物生長[15]，平均插值結(jié)果見表4；SWAP 模型的上邊界條件為大氣邊界，下邊界條件為地下水埋深；試驗(yàn)前的土壤含水率作為初始條件輸入模型。

表4 作物生長階段（DVS）與作物生理指標(biāo)的關(guān)系Tab.4 Relationship between crop growth stage(DVS)and crop physiological indexes

利用率定后的模型模擬不同水平年下不同灌溉定額的作物產(chǎn)量，結(jié)合作物耗水量計(jì)算水分利用效率。水分利用效率（WUE）是農(nóng)田蒸散消耗單位重量水所制造的干物質(zhì)量，可以用來表征作物的節(jié)水灌溉水平：

式中：Yact為產(chǎn)量；Tact為作物耗水量，用蒸騰量ET來表示。

2.3 模型評(píng)價(jià)標(biāo)準(zhǔn)

模擬值和實(shí)測(cè)值的吻合度采用均方誤差（RMSE）和平均相對(duì)誤差（MRE）來評(píng)價(jià)：

式中：Pi為第i個(gè)模擬值；Qi為第i個(gè)實(shí)測(cè)值；N為觀測(cè)值個(gè)數(shù)。RMSE與MRE值越小表示模型模擬效果越好。

3 結(jié)果與討論

3.1 模型率定與驗(yàn)證

利用冬小麥的G3-1 處理和夏玉米的M3-1 處理對(duì)土壤含水率和實(shí)際產(chǎn)量進(jìn)行率定，首先使用初始土壤參數(shù)和給定的溶質(zhì)擴(kuò)散系數(shù)、水分作物阻力系數(shù)等參數(shù)進(jìn)行模擬，輸入數(shù)據(jù)，將模擬出的土壤含水率和產(chǎn)量與實(shí)測(cè)值對(duì)比，適當(dāng)調(diào)整參數(shù)，使模擬值與實(shí)測(cè)值的差距最小。

3.1.1 土壤含水率

冬小麥（G3-1 處理）播種后第17、154、166、186、220和235 d，土壤不同深度處含水率（0～20、20～40、40～60、60～80、80～100 cm）隨時(shí)間的變化如圖1（a）所示。由圖1（a）可知，冬小麥由于接受了苗期、抽穗和灌漿期的灌溉，各層土壤含水率均有明顯上升，0～20、20～40、40～60、60～80和80～100 cm 深度土壤含水率率定過程的RMSE分別為0.837、0.912、0.782、0.771 和0.634；MRE分別為3.480%、3.750%、3.719%、2.839%和2.687%。夏玉米（M3-1 處理）播種后第7、35、51、66、72、79、86 和96 d，土壤不同深度處含水率（0～20、20～40、40～60、60～80、80～100 cm）隨時(shí)間的變化如圖1（b）所示。由圖1（b）可知，夏玉米由于生長期內(nèi)降雨充足，土壤含水率波動(dòng)程度更大，0～20、20～40、40～60、60～80 和80～100 cm 深度土壤含水率率定過程的RMSE分別為0.012、0.010、0.018、0.016 和0.016；MRE分別為5.188%、3.573%、6.197%、4.560%和4.430%。結(jié)果表明2 種作物土壤含水率實(shí)測(cè)值和模擬值的變化趨勢(shì)一致，吻合度較高。

圖1 率定土壤含水率結(jié)果Fig.1 Results of soil water content calibration

利用率定后的參數(shù)對(duì)冬小麥（G3-2）和夏玉米（M3-2）的土壤含水率進(jìn)行驗(yàn)證，結(jié)果如圖2所示。驗(yàn)證結(jié)果表明，冬小麥0～20、20～40、40～60、60～80和80～100 cm深度土壤含水率驗(yàn)證過程的RMSE分別為0.897、1.041、0.449、0.774和0.832；MRE分別為3.504%、4.068%、3.224%、3.019%和3.299%。夏玉米0～20、20～40、40～60、60～80和80～100 cm深度土壤含水率驗(yàn)證過程的RMSE分別為0.009、0.009、0.015、0.010和0.013；MRE分別為3.832%、3.471%、5.536%、3.010%和3.557%。土壤含水率的率定與驗(yàn)證過程中RMSE較小，均不超過1.0，MRE基本上在6%以下。不同灌水方案下土壤含水率模擬值與實(shí)測(cè)值吻合較好，模擬值基本上反映了實(shí)測(cè)值的變化趨勢(shì)。

圖2 土壤含水率模擬結(jié)果Fig.2 Simulation results of soil water content

3.1.2 產(chǎn) 量

利用冬小麥（G3-1處理）和夏玉米（M3-1處理）的產(chǎn)量進(jìn)行率定。將作物產(chǎn)量的模擬值和實(shí)測(cè)值進(jìn)行比較，并對(duì)參數(shù)進(jìn)行適當(dāng)調(diào)整，使得模擬值更加接近實(shí)測(cè)值。本次模型中采用的是簡單作物模塊模擬產(chǎn)量，無法直接得出產(chǎn)量的模擬值，而是給出相對(duì)產(chǎn)量，即實(shí)際產(chǎn)量與潛在產(chǎn)量的比值（Ya/Yp）。假定2019年G3-1 處理得到冬小麥產(chǎn)量（12 763.27 kg/hm2）和M3-1 處理得到的夏玉米產(chǎn)量（7 489.21 kg/hm2）為最大實(shí)際產(chǎn)量，模擬產(chǎn)量根據(jù)模擬的相對(duì)產(chǎn)量與最大實(shí)際產(chǎn)量的換算公式進(jìn)行計(jì)算。冬小麥和夏玉米率定和驗(yàn)證結(jié)果如圖3所示。

圖3 作物產(chǎn)量率定和驗(yàn)證結(jié)果Fig.3 Results of crop yield calibration and validation

由率定和驗(yàn)證結(jié)果可知，作物產(chǎn)量的實(shí)測(cè)值和模擬值較為接近。冬小麥率定過程的RMSE和MRE分別為439.430 kg/hm2、3.443%，驗(yàn)證過程的RMSE和MRE分別為583.696 kg/hm2、6.215%；夏玉米率定過程的RMSE和MRE分別為172.919 kg/hm2、2.309%，驗(yàn)證過程的RMSE和MRE分別為238.931 kg/hm2、4.091%。作物產(chǎn)量率定與驗(yàn)證過程中的RMSE值均不超過最大實(shí)際產(chǎn)量的10%，且MRE均在7%以下，符合誤差精度要求。這說明不同灌水方案下作物產(chǎn)量的模擬值與實(shí)測(cè)值吻合較好，模擬值能夠反映實(shí)測(cè)值的變化趨勢(shì)，模型可靠。

3.2 灌溉制度優(yōu)化

小麥拔節(jié)至抽穗期、玉米的拔節(jié)至抽雄期，是最重要的生長階段。此階段植株生長量劇增，耗水量大增，若水分不足則大幅影響產(chǎn)量[16]。基于優(yōu)化的目的，設(shè)置了以下灌水方案：小麥灌3 水方案選擇在拔節(jié)、抽穗、灌漿期，4 水選擇在苗期、拔節(jié)、抽穗和灌漿期；玉米灌3水選擇苗期、拔節(jié)和抽雄期，2水選擇在拔節(jié)和抽雄期。各生育期內(nèi)的水量分配情況依據(jù)黃靜[17]得出的作物生育階段敏感系數(shù)來確定。根據(jù)當(dāng)?shù)貧庀筚Y料，平、枯降雨年型下，冬小麥生育期降雨量分別為157、135 mm，玉米生育期內(nèi)降水量分別388、324 mm。冬小麥和夏玉米生育期內(nèi)需水量分別為500 和450 mm。每種降雨年型下，分別設(shè)置ETc（作物需水量）、80%ETc和60%ETc條件下的灌水定額，利用率定好的SWAP模型模擬產(chǎn)量，并計(jì)算水分利用效率WUE。具體灌水方案見表5。

表5 不同降雨年型冬小麥-夏玉米優(yōu)選灌溉制度Tab.5 Optimal irrigation system for winter wheat-summer maize in different rainfall years

由于灌溉制度的不同，土壤水分狀況能在一定程度上體現(xiàn)作物的水分脅迫情況。不同降雨年型下冬小麥和夏玉米生育期內(nèi)0～100 cm 深度平均土壤含水率的模擬結(jié)果如圖4和圖5所示。結(jié)果表明，不同灌溉處理的土壤含水率總體上表現(xiàn)為缺水情況越嚴(yán)重，土壤含水率越低。由于冬小麥4 水和3 水的灌溉處理僅區(qū)別于苗期是否灌水，而苗期根系較淺，水分脅迫并不明顯。夏玉米3 水和2 水的灌溉處理區(qū)別是灌漿期是否灌水，3水灌溉情況下，灌漿期后土壤含水率有所提升。有研究表明，水分脅迫并非完全是負(fù)效應(yīng)，特定階段、有限的水分脅迫對(duì)提高產(chǎn)量是有益的[18]，因此灌溉制度的選取也要考慮其產(chǎn)量和水分利用效率。

圖4 不同降雨年型下冬小麥生育期土壤含水率變化Fig.4 Changes of soil moisture content during winter wheat growth stage under different rainfall years

圖5 不同降雨年型下夏玉米生育期土壤含水率變化Fig.5 Changes of soil moisture content during summer maize growth stage under different rainfall years

由表5可知，與ETc相比，其他情況下相對(duì)產(chǎn)量有所降低，但下降幅度微小。平水年冬小麥在灌4 水的方案下，80%ETc相對(duì)產(chǎn)量僅比ETc降低1.7%；夏玉米在灌2 水的方案下，80%ETc的相對(duì)產(chǎn)量僅比ETc降低1.9%，ETc情況下利用水量增多，產(chǎn)量并未顯著提高，反而降低了灌水的邊際效益[4]。灌溉水量為60%ETc時(shí)，蒸發(fā)量減少，水分利用效率有所提升，但相對(duì)產(chǎn)量均有所降低。作物方面，夏玉米的WUE普遍低于冬小麥，這主要是由于夏玉米生育期內(nèi)蒸散發(fā)強(qiáng)烈，耗水量大，導(dǎo)致WUE計(jì)算值較小。

灌4水的方案下，冬小麥的相對(duì)產(chǎn)量和水分利用效率較灌3水整體偏低，多一次灌溉未達(dá)到高效增產(chǎn)的目的，且增加了灌水成本。同樣，夏玉米灌2水的相對(duì)產(chǎn)量和水分利用效率普遍大于灌3水的結(jié)果。因此，最終選取產(chǎn)量和水分利用效率相協(xié)調(diào)的80%ETc條件，對(duì)冬小麥拔節(jié)、抽穗和灌漿期進(jìn)行灌溉，對(duì)夏玉米拔節(jié)和抽雄期進(jìn)行灌溉為最優(yōu)方案。

本模擬結(jié)果表明作物產(chǎn)量的實(shí)測(cè)值與模擬值變化規(guī)律基本一致，誤差在可接受范圍內(nèi)，模型基本上能夠反映作物生長過程。但在配水過程中的灌溉需水量僅簡單地以作物需水與降水量的差值計(jì)算，而實(shí)際過程是復(fù)雜而多樣的，未來應(yīng)探索更為詳實(shí)的計(jì)算方法，提高準(zhǔn)確度。同時(shí)，康紹忠[18]提出降雨年際和季節(jié)性分配不均，灌溉制度也要依賴于此。孫晉鍇[1]依據(jù)豫東地區(qū)多年降水資料，得到玉米灌溉的最優(yōu)解為拔節(jié)期至抽雄期、抽雄期至灌漿期進(jìn)行2 次灌水， 灌水定額分別為85 和75 mm，相對(duì)產(chǎn)量達(dá)到0.982，然而水分利用效率并未作為其優(yōu)化灌溉制度的選取指標(biāo)。呂夢(mèng)醒[19]在研究濁漳河干流灌溉制度的優(yōu)化時(shí)，首先保證水分利用效率和產(chǎn)量，綜合考慮下泄徑流量，得出了更為完善的灌溉制度。Jensen[20]提出作物本身的生理指標(biāo)較土壤水分狀況能夠更好地指示作物是否發(fā)生水分虧缺，因此未來考慮水分脅迫可以從作物生理方面入手。

提高農(nóng)田水利效率，需要經(jīng)濟(jì)的灌水定額、合理的田間管理、高效的水利設(shè)施，更需要開展灌溉用水的優(yōu)化配置，從而達(dá)到高效、增產(chǎn)又節(jié)水的目的[21]。

4 結(jié) 論

利用SWAP 模型，對(duì)桓臺(tái)縣試驗(yàn)站2018-2019年不同灌水方案下冬小麥、夏玉米的土壤含水率和產(chǎn)量進(jìn)行模擬，并根據(jù)模擬結(jié)果得到優(yōu)化后的作物灌溉配水方案，主要結(jié)論如下。

（1）利用2019年冬小麥G3-1 處理和夏玉米M3-1 處理的土壤含水率和產(chǎn)量數(shù)據(jù)進(jìn)行率定，其他處理進(jìn)行驗(yàn)證，RMSE和MRE檢驗(yàn)冬小麥和夏玉米率定與驗(yàn)證過程。結(jié)果表明：土壤含水率RMSE較小，均不超過1.0，MRE基本上在6%以下；作物產(chǎn)量RMSE值均不超過最大實(shí)際產(chǎn)量的10%，且MRE均在7%以下。符合誤差精度要求。不同灌水方案下土壤含水率和作物產(chǎn)量的模擬值與實(shí)測(cè)值吻合較好，模擬值基本上反映了實(shí)測(cè)值的變化趨勢(shì)，模型可靠。

（2）利用率定好的SWAP 模型模擬ETc、80%ETc、60%ETc下不同降雨年型的作物相對(duì)產(chǎn)量和水分利用效率。結(jié)果表明，冬小麥灌3 水、夏玉米灌2 水能夠提高產(chǎn)量并且顯著提升水分利用效率。最優(yōu)灌溉制度為80%ETc下，平水年，小麥在拔節(jié)、抽穗和灌漿期單次灌水定額分別為87.3、124.7、62.4 mm，枯水年分別為92.9、132.7、66.4 mm；平水年，玉米在拔節(jié)和抽雄期單次灌水定額分別為19.8、29.8 mm，枯水年分別為40.2、60.6 mm。