張鐵楠, 許為政, 魏 湜, 顧萬榮, 賀 丹, 蘆玉雙, 錢 誠, 陳 琦, 李 晶

(1.東北農業(yè)大學 農學院, 哈爾濱150030; 2.大慶市農業(yè)委員會, 黑龍江 大慶 163000;3.黑龍江省農業(yè)技術推廣站, 哈爾濱150030; 4.黑龍江省農業(yè)科學院綏化分院, 黑龍江 綏化 152000)

?

基于WOFOST模型的東北地區(qū)春小麥水分平衡模擬與驗證

張鐵楠1,2, 許為政3, 魏 湜1, 顧萬榮1, 賀 丹1, 蘆玉雙3, 錢 誠4, 陳 琦4, 李 晶1

(1.東北農業(yè)大學 農學院, 哈爾濱150030; 2.大慶市農業(yè)委員會, 黑龍江 大慶 163000;3.黑龍江省農業(yè)技術推廣站, 哈爾濱150030; 4.黑龍江省農業(yè)科學院綏化分院, 黑龍江 綏化 152000)

對WOFOST模型進行了本地化參數(shù)校正，并對輸入模型的參數(shù)進行了敏感性分析，選取出對結果影響較大的幾個參數(shù)進行調整，建立相應的模型運行數(shù)據(jù)庫，對東北春麥區(qū)具有代表性的哈爾濱、綏化、嫩江、呼瑪和牙克石地區(qū)2012年、2013年春小麥生育期內土壤水分含量情況進行了模擬及精度分析，最終評價WOFOST模型在東北地區(qū)春小麥生長過程中水分平衡模擬的適應性及準確性。結果表明：WOFOST模型對2012年、2013年各地春小麥生育期內土壤水分含量變化模擬值與實測值的絕對平均誤差為2.43%～8.75%，均方根誤差RMSE為6.75%～13.75%，模擬性能指數(shù)EF為-0.844～0.245，殘差聚集指數(shù)CRM為-0.272～0.084，一致性系數(shù)IoA為0.457～0.789，在可信區(qū)間內。從整體看，該模型對土壤水分含量的模擬具備一定的參考價值；同時利用校正后的模型對2012年、2013年各地春小麥作物蒸騰、土壤蒸發(fā)與作物需水量進行了模擬，實現(xiàn)了利用較少參數(shù)達到作物蒸騰與土壤蒸發(fā)分離的目的。

WOFOST模型; 東北地區(qū); 春小麥; 參數(shù)調試; 土壤含水量; 模擬驗證

東北地區(qū)是中國緯度位置最北的區(qū)域，包括黑龍江、吉林、遼寧和內蒙古東部。該地區(qū)雨熱同季，屬溫帶濕潤、半濕潤大陸性季風氣候，地勢平坦、土壤肥沃，以典型黑土、黑鈣土著稱，是中國重要的糧食生產地區(qū)。2012年東北三省(黑龍江、吉林和遼寧)糧食總產量達到11 175萬t，占全國糧食產量的18.95%[1]。該區(qū)為傳統(tǒng)雨養(yǎng)農業(yè)，氣候變化是影響春麥生產及產量的重要因子。長期以來，該區(qū)農業(yè)生產仍處于開放式經營模式，導致生產潛力尚未發(fā)揮其應有水平，缺少數(shù)字模擬技術對生產預測預警的應用研究。

WOFOST模型是在世界糧農組織FAQ的資助下，由荷蘭的瓦赫寧根大學de Wit教授在作物模型理論基礎上不斷完善開發(fā)而成的，是模擬特定的土壤和氣候條件下一年生作物生長的動態(tài)解釋性模型[2]。WOFOST模型的首次在3個非洲國家被應用[3]。秘魯?shù)陌驳谒股矫}地區(qū)應用該模型對灌溉和水分的保持進行了評價[4]，歐洲把它用于主要大田作物的潛在產量評估[5-6]。布基納法索的AGRISK工程中用其進行產量風險研究。在MARS工程[7-8]中，歐洲聯(lián)合研究中心把WOFOST作為產量預測的模型。開發(fā)應用后的30年間，我國多地區(qū)應用WOFOST模型進行模擬驗證試驗[9-10]，寒地作物生長模擬涉及水稻、小麥、玉米、大豆、甜菜、馬鈴薯等[11-13]。王光火，謝文霞[14]等人根據(jù)中國浙江具體情況對WOFOST 模型的部分參數(shù)進行了修改，對該地區(qū)的部分一年生作物進行了產量評估。于強等[15]運用WOFOST模型進行了小麥的水分處理試驗，用來評價該模型在華北平原的適用性。杜春英等[16]運用WOFOST模型對黑龍江省水稻產量動態(tài)預報進行研究，實現(xiàn)了WOFOST 模型在水稻產量動態(tài)預報業(yè)務中的應用。

本研究對WOFOST模型進行了本地化參數(shù)校正，并對輸入模型的參數(shù)進行了敏感性分析，選取出對結果影響較大的幾個參數(shù)進行調整，建立相應的模型運行數(shù)據(jù)庫，對東北春麥區(qū)具有代表性的哈爾濱、綏化、嫩江、呼瑪和牙克石地區(qū)2012年，2013年春小麥生育期內土壤水分含量情況進行了模擬及精度分析，最終評價WOFOST模型在東北地區(qū)春小麥生長過程中水分平衡模擬的適應性及準確性，以其能為政府決策部門進行農業(yè)生產宏觀調控提供理論依據(jù)，進而保障東北地區(qū)糧食生產的安全與穩(wěn)定。

1　試驗設計與方法

1.1試驗區(qū)自然概況

試驗選擇東北春麥區(qū)5個具有代表性的監(jiān)測地點(表1)，獲取2010—2013年春小麥生育期間田間自然狀態(tài)下作物數(shù)據(jù)及氣象數(shù)據(jù)。

表1　試驗監(jiān)測點概況

1.2試驗材料與設計

試驗品種采用龍麥26，采用戊唑醇種衣劑拌種，基本保苗數(shù)600萬株/hm2，播深3～5 cm，播后鎮(zhèn)壓；應用除草劑封閉除草；施肥同普通大田水平，其他均為常規(guī)大田管理。

1.3數(shù)據(jù)觀測與方法

作物資料：2012年，2013年作物資料來自東北農業(yè)大學小麥栽培研究室及黑龍江其他各站點多年研究記錄數(shù)據(jù)。

氣象資料：2012年，2013年春小麥生育期內逐日氣象數(shù)據(jù)由中國農科院聯(lián)合清華大學、北京農業(yè)信息技術研究中心等機構聯(lián)合開發(fā)的“農業(yè)環(huán)境無線遠程實時監(jiān)控系統(tǒng)”采集記錄。系統(tǒng)配有HAYASHI DENKO公司的Pt100空氣溫度傳感器和土壤溫度傳感器，Automata公司的AQUA-TEL-TDR空氣濕度傳感器和土壤濕度傳感器，還配有太陽輻射測定裝置及風向風速測定裝置，每10 min自動測定試驗地氣候狀況和土壤溫度、濕度變化情況及空氣溫度、濕度變化及太陽輻射情況等上傳至服務器。數(shù)據(jù)缺失部分由中國氣象科學數(shù)據(jù)共享服務網上的數(shù)據(jù)進行補充。試驗站的經度、緯度、海拔高度數(shù)據(jù)通過GPS測量。

土壤資料：本研究所用土壤參數(shù)除部分實測外，其余主要取自《中國土種志》[17]和《黑龍江土壤》[18]；土壤含水量采用“農業(yè)環(huán)境無線遠程實時監(jiān)控系統(tǒng)”土壤濕度傳感器進行測定，每10 min測一次，計算每天土壤濕度平均值。

1.4WOFOST數(shù)據(jù)庫的建立與參數(shù)調試

氣象數(shù)據(jù)庫的建立：常規(guī)觀測數(shù)據(jù)包括2012年，2013年東北地區(qū)各檢測站點設備自動采集，人工整理的逐日最高溫度、最低溫度、日照時數(shù)、降水、平均風速和水汽壓等要素數(shù)據(jù),建立模型可識別的氣象數(shù)據(jù)庫。

表2　各試驗地土壤理化性質

作物數(shù)據(jù)庫的建立：WOFOST模型作物模塊共計18個參數(shù)，其中7個默認參數(shù)：TSUMEN，TSUM1，TSUM2，TDWI，SLATB，SPAN，AMAXTB，11個需調參CLV，CVS，CVO 、CVR，RGRLAI(最大葉面積增長量)、KDIFTB(散射光消光系數(shù))、EFFTB(單葉光利用效率)、RML，RMO，RMR，RMS，將參數(shù)輸入模型，建立作物數(shù)據(jù)庫。

土壤數(shù)據(jù)庫的建立：11個參數(shù)，8個需調參：SMTAB(土壤含水量)、SMW(凋萎系數(shù))、SMFCF(田間持水量)、SMO(土壤飽和含水量)、CONTAB(土壤導水率)、KO(飽和土壤導水率)、RDI(初始扎根深度)、RDM(最大扎根深度)，3個默認參數(shù)：CRAIRC(土壤通氣臨界空氣含量)、SOPE(根層最大入滲速率)、KSUB(根層以下土壤最大入滲速率)，參數(shù)輸入模型，建立土壤數(shù)據(jù)庫。

準確的參數(shù)是模擬結果準確和穩(wěn)定的前提，需根據(jù)歷史資料，需要根據(jù)實際情況對作物參數(shù)進行調整。

1.5檢驗模型有效性的方法和指標

為評價各項模擬值與實測值之間的差距和模擬效果，選擇以下參數(shù)進行評價:線性回歸方程R2，平均相對誤差(AE)、相對平均誤差(MAE)、均方根誤差(RMSE)、模擬性能指數(shù)(EF)、殘差聚集系數(shù)(CRM)、一致性系數(shù)(IoA)。計算公式如下:[19-21]

R2越接近1，說明模擬值與實測值的相關性越好。AE，MAE越小，說明模型的模擬性越好。RMSE取值范圍≥0，越靠近0說明模擬效果越好。EF取值范圍≤1，取值越接近1表明模擬效果越好。CRM可為正值或負值，正值表示相對實測值偏高，負值表示相對實測值偏低。IoA越大，表示誤差越小，模型的預測性越好。

2　結果與分析

運用調試后的WOFOST模型對哈爾濱、綏化、嫩江、呼瑪、牙克石等地2012年，2013年的春小麥生育期內土壤水分含量變化進行模擬，結合當時實測的試驗數(shù)據(jù)，對比分析模型對水分平衡模擬結果的有效性，同時對小麥生育期內作物需水量、作物蒸騰量、土壤蒸發(fā)量進行模擬，最終評價WOFOST模型對東北地區(qū)春小麥水分平衡模擬的適用性。

2.1土壤含水量模擬結果及精度分析

運用WOFOST模型對2012年、2013年哈爾濱、綏化、嫩江、呼瑪、牙克石等地的春小麥生育期內土壤水分含量進行模擬，與小麥苗情檢測系統(tǒng)的土壤濕度的實測值進行對比分析，如圖1所示。

從圖1中對土壤水分含量的模擬值與實測值的對比可以看出，WOFOST模型基本能反映出小麥生育期內土壤含水量整體變化趨勢，模擬值與實測值之間的誤差在允許范圍之內，土壤濕度的模擬值整體低于實測值。

圖1　2012-2013年各地春麥生育期內土壤體積含水量實測值與模擬值

表3　WOFOST模型對土壤水分含水量模擬結果精度分析

從表3中可以看出，WOFOST模型對2012年，2013年各地春小麥生育期內土壤水分含量變化模擬值與實測值的絕對平均誤差在2.43%～8.75%，RMSE在6.75%～13.75%，EF在-0.844～0.245，CRM在-0.272～0.084，IoA在0.457～0.789，在可信區(qū)間內。從整體看，該模型對土壤水分含量的模擬具備一定的參考價值。

2.2生育期內作物需水量、作物蒸騰量、土壤蒸發(fā)量的模擬

運用WOFOST模型對各地春小麥生育期內作物需水量、作物蒸騰量、土壤蒸發(fā)量進行模擬，作物蒸散及區(qū)分作物蒸騰與土壤蒸發(fā)量的準確模擬，是預測生物量的關鍵。在生物量模擬達到一定準確度的前提下，作物蒸騰與土壤蒸發(fā)的模擬結果具有一定的可信度。2012年和2013年黑龍江各監(jiān)測站點的生育期內土壤蒸發(fā)、作物蒸騰、作物需水量模擬結果如圖2所示。

從 圖2 中可以看出在作物生長發(fā)育前期，土壤中的水分消耗主要以土壤水分蒸發(fā)為主，原因是由于生長發(fā)育前期植株剛出苗，對地面的覆蓋度較小，土壤直接接受太陽直照及風力侵蝕，土壤水分蒸發(fā)較快，哈爾濱、綏化、嫩江、呼瑪、牙克石各站點日平均土壤蒸發(fā)量分別為為1.8，1.81，1.99，1.96，2.62 mm/d，出苗后，隨著植株逐漸長大，對地面的覆蓋率逐漸增大，減少太陽直接照射和阻擋風力，植株對地面起到冠層覆蓋作用，使土壤水分蒸發(fā)強度逐漸下降，出苗到成熟期日平均土壤蒸發(fā)量分別為0.98，0.73，0.94，0.88，0.69 mm/d。

作物蒸騰量在整個生育期的變化曲線均呈脈沖狀，主要是不同天氣和氣象條件影響的結果。每次降雨后的晴天，日蒸騰量和蒸騰速率有明顯的上升趨勢，而后依次減少。從出苗到成熟的生長季節(jié)中，小麥的蒸騰量隨時間的推移發(fā)生明顯的變化。一般趨勢是苗期較低,拔節(jié)后逐漸上升，明顯高于苗期，抽穗到灌漿期蒸騰量繼續(xù)上升，達到高峰階段，日平均蒸騰量為4.32 mm/d，灌漿以后逐漸下降。經過模擬所得，各個站點生育期內日平均作物蒸騰量分別為2.75，3.24，3.2，3.08，2.72 mm/d。

因降雨、溫度等因素的綜合影響，作物需水量曲線在生育期內波動幅度較為頻繁，生育期內作物日平均需水量的模擬結果分別為3.74，3.98，4.16，3.96，3.41 mm/d。

從總體來看，在春小麥整個生育期內，隨著春小麥冠層覆蓋度的日益增大，土壤蒸發(fā)量下降幅度增大，在此后一段時間內，土壤蒸發(fā)量變化波動較小，但遇到強降雨之后的土壤蒸發(fā)量依然變化強烈，生育末期由于植株葉片衰老，地面覆蓋度減小，土壤水分蒸發(fā)較為活躍。受葉片衰老與溫度的影響，作物需水量、作物蒸騰和土壤蒸發(fā)三者在成熟期末整體均呈下降趨勢。通過模型的模擬，哈爾濱、綏化、嫩江、呼瑪、牙克石各站點生育期內每日的土壤蒸發(fā)量變化區(qū)間為0.85～2.95，0.15～3.9，0.15～3.8，0.1～4.1，0.15～3.15 mm/d；每日的作物蒸騰量變化區(qū)間為0.25～5.7，0.5～6.05，0.15～6.95，0.25～5.9，0.15～6.15 mm/d；每日的作物需水量變化區(qū)間為0.95～6.85，1.15～6.85，0.95～7.1，1.1～6.7，0.55～6.4 mm/d。

3　結論與討論

3.1結 論

本研究采用模型模擬和田間試驗相結合的方法，利用整理好的各地氣象數(shù)據(jù)、土壤數(shù)據(jù)和生育期內各種作物觀測資料，對WOFOST模型進行參數(shù)校準，利用本地化調參后的WOFOST模型對黑龍江哈爾濱、綏化、嫩江、呼瑪和內蒙古牙克石地區(qū)五個地點春小麥生育期內土壤含水量變化過程進行模擬，研究結果表明，模擬結果與實測結果較為接近，各參數(shù)指標誤差均在可信區(qū)間，說明WOFOST模型可用于區(qū)域范圍內土壤含水量變化的研究，對農業(yè)水分利用和預測方面具有一定的參考價值；同時通過對各地春小麥作物蒸騰、土壤蒸發(fā)與作物需水量進行了模擬，實現(xiàn)了利用較少參數(shù)達到作物蒸騰與土壤蒸發(fā)分離的目的。

圖2　2012年、2013年各地春麥生育期內土壤蒸發(fā)、作物蒸騰、作物需水量模擬

3.2討 論

WOFOST模型中，與作物生長發(fā)育相關的參數(shù)有數(shù)十個，加之模型參數(shù)之間的相關性，對全部的參數(shù)一一進行校正，其工作量巨大，且是不現(xiàn)實的[22]。因此，需要針對其中較為重要的一部分參數(shù)進行調整。首先需要根據(jù)研究目的，對各主要參數(shù)進行敏感性分析，選取出對結果影響較大的幾個參數(shù)，經檢驗WOFOST模型最終的模擬結果以及過程最大值對AMAXTB2，CVO，CVS，KDIFTB，Q10，SPAN，TSUM1和TSUM2這幾個參數(shù)的敏感性都較強，均超過了1%，最高的SPAN甚至達到了223.35%，因此需要根據(jù)資料，對這些參數(shù)進行調整。杜春英等[19]運用WOFOST模型在黑龍江對水稻產量預報進行研究，結果表明水稻單產的擬合準確率最大值為100%，最小值為84%，平均值為94.3%，產量預測的精度與本研究運用WOFOST模擬春小麥產量的精度基本一致，模擬結果較為可靠，但還有提升的空間，需要對此模型進行進一步的參數(shù)調整及驗證工作，增強模型在本地模擬的穩(wěn)定性及準確性。

總體來說，經過校正的WOFOST模型在東北地區(qū)春小麥的生產中水分平衡的模擬具有較強的準確性及適應性，可以用于指導東北地區(qū)春小麥的生產，具有一定的參考價值。

[1]中國統(tǒng)計年鑒. http:∥www. stats. gov. cn/tjsj/ndsj/2013/indexch. Htm.2013.

[2]Boogaarcl H L, Diepen C A V, Rotter R P, et al. WOFOST 7.1: user's guide for the WOFOST 7.1 crop growth simulation model and WOFOST Control Center 1.5[M]. Wageningen: DLO Winand Staring Centre, 1998.

[3]謝文霞,王光火,張奇春.WOFOST模型的發(fā)展及應用[J].土壤通報,2006,37(1):154-158.

[4]Berkhout J A A, Huygen J, Azzali S, et al. MARS definition study: results of the preparatory phase [R]. Wgeningen,1988.

[5]Van Lanen H A J, Van Diepen C A, Reinds G J, et al. Physical land evaluation methods and GIS to explore the crop growth potential and its effects within the European Communities[J]. Agricultural Systems, 1992,39(3):307-328.

[6]De Koning G H J, Janssen H, Van Keulen H. Input and output coefficients of various cropping and livestock systems in the European Communities[M]. The Hague: Netherlands Scientific Council for Government Policy, 1992.

[7]Diepen C A, Van Der Wal T. Crop growth monitoring and yield forecasting at regional and national scale[M]. Publication EUR, 1994.

[8]Vossen P, Rijks D, Hough M N. Early crop yield assessment of the EU countries: the system implemented by the Joint Research Centre[M]. Office for Official Publications of the European Communities, 1995.

[9]張雪芬,佘衛(wèi)東,王春乙,等.WOFOST模型在冬小麥晚霜凍害評估中的應用[J].自然災害學報,2006,15(6):337-341.

[10]馬玉平,王石立,張黎,等.針對華北小麥越冬的WOFOST模型改進[J].中國農業(yè)氣象,2005,26(3):145-149.

[11]高永剛,那濟海,顧紅,等.黑龍江省馬鈴薯氣候生產力特征及區(qū)劃[J].中國農業(yè)氣象,2007,28(3):275-280.

[12]李秀芬,王育光,季生太,等.作物長勢監(jiān)測系統(tǒng)(CGMS)在黑龍江省的應用[J].中國農業(yè)氣象,2005,26(3):155-157.

[13]高永剛,南瑞,顧紅,等.黑龍江省甜菜氣候生產力模擬和種植氣候區(qū)劃[J].生態(tài)學雜志,2009,28(1):27-31.

[14]謝文霞,嚴力蛟,王光火.運用WOFOST模型對浙江水稻潛在生長過程的模擬與驗證[J].中國水稻科學,2006,20(3):319-323.

[15]于強,羅毅,鄔定榮,等.作物生長模型WOFOST在華北平原的適用性研究[J].植物生態(tài)學報,2004,27(5):594-602.

[16]杜春英,李秀芬,王晨軼等.基于WOFOST模型的黑龍江省水稻產量動態(tài)預報研究[J].安徽農業(yè)科學,2011,39(24):15093-15095.

[17]全國土壤普查辦公室.中國土種志[M].北京:中國農業(yè)出版社,1991.

[18]黑龍江省土地管理局,黑龍江省土壤普查辦公室.黑龍江土壤[M].北京:農業(yè)出版社,1992.

[19]張忠學,于貴瑞.不同灌水處理對冬小麥生長及水分利用效率的影響[J].灌溉排水學報,2003,22(2):1-4.

[20]商磊,趙軍,祁廣云,等.黑土農田大豆產量形成過程的模擬驗證[J].中國生態(tài)農業(yè)學報,2008,16(4):869-873.

[21]王宗明,張柏,宋開山,等.CropSyst作物模型在松嫩平原典型黑土區(qū)的校正和驗證[J].農業(yè)工程學報,2005,21(5):47-50.

[22]姜志偉,陳仲新,周清波,等.CERES-Wheat作物模型參數(shù)全局敏感性分析[J].農業(yè)工程學報,2011,27(1):236-242.

Simulation and Verification of Spring Wheat Water Balance in Northeast Region Based on WOFOST Model

ZHANG Tienan1,2, XU Weizheng3, WEI Shi1, GU Wanrong1, HE Dan1,LU Yushuang3, QIAN Cheng4, CHEN Qi4, LI Jing1

(1.College of Agriculture, Northeast Agricultural University, Harbin 150030, China;2.CommissionofAgricultureofDaqingCity,Daqing，Heilongjiang163000,China; 3.HeilongjiangProvinceStationforPopularizingAgriculturalTechniques,Harbin150030,China; 4.BranchofSuihuaAcademyofAgriculturalSciencesofHeilongjiangProvince,Suihua，Heilongjiang152000,China)

This research has carried on the localization parameter calibration in WOFOST model, and carried the sensitivity analysis of the input relevant parameters. Then several parameters were picked out which had significant influence on result and made different adjustment scheme. A corresponding model was setted up to run database to simulate spring wheat growth and soil moisture content conditions in Harbin, Suihua, Nenjiang, Huma and Yakeshi area in 2012, 2013. Then the measured values were compared with the simulation values. Finally, the WOFOST model′s adaptability and accuracy of spring wheat in the process of water balance growth in northeast area was evaluated. The results showed that for 2012 and 2013, absolute error of between the soil moisture content simulated by WOFOST model and the measured values ranged from 2.43% to 8.75% on average, RMSE was between 6.75% to 13.75%, EF was between -0.844 to 0.245, CRM was between -0.272 to 0.084, IoA was between 0.457 to 0.789 that was in the confidence interval. On the whole, the model simulation of soil water content had a certain reference value. The calibrated model was used to simulate the spring wheat soil evaporation and crop transpiration and crop water requirement in 2012, 2013. It can be achieved the purpose of using the fewer parameters for soil evaporation and crop transpiration separation.

WOFOST model; northeast area; spring wheat; parameter adjustment; soil moisture content; simulation and verification

2014-08-21

2014-10-10

公益性行業(yè)(農業(yè))科研專項(200903010);黑龍江省博士后基金資助項目(LBH-Z10253);東北農業(yè)大學博士科研資助項目(2010RCB20)

張鐵楠(1989—),男,黑龍江省大慶市人，碩士,從事作物高產栽培研究。E-mail:15846000034@163.com

李晶(1977—),女,黑龍江省哈爾濱市人，博士,副教授,從事作物高產栽培研究。E-mail:jingli1027@163.com

S512.1+2

A

1005-3409(2015)04-0045-07