趙 藝，鄒雨伽，張玉芳，3

（1.中國氣象局成都高原氣象研究所/高原與盆地暴雨旱澇災(zāi)害四川省重點實驗室，成都 610072；2.四川省農(nóng)業(yè)氣象中心，成都 610072；3.南方丘區(qū)節(jié)水農(nóng)業(yè)研究四川省重點實驗室，成都 610066）

中國是油料消費大國，油料供給不足情況一直未得到有效解決，近年來油料的消費增速高于供給增速，導(dǎo)致供給日趨緊張［1］。油料是中國對國際市場需求量最大的農(nóng)產(chǎn)品，過度的進口對國家食用油供給安全造成威脅［1］。目前，中國生產(chǎn)和消費的油料主要以草本油料作物為主，主要包括油菜（Brassica napus L.）、花生、大豆、芝麻、向日葵等，其中油菜是中國第一大油料作物，種植面積及產(chǎn)量均位居世界首位［1-3］。四川省作為全國油菜生產(chǎn)大省，近年來產(chǎn)業(yè)穩(wěn)定發(fā)展并逐步提升，播種面積及產(chǎn)量均位列全國第一。四川省油菜以冬油菜為主，其中四川盆區(qū)的種植面積和總產(chǎn)占主要地位。因此，針對四川盆區(qū)開展油菜產(chǎn)量預(yù)報對保障國家糧食安全具有重要意義。

近年來，國內(nèi)外學(xué)者對作物產(chǎn)量動態(tài)預(yù)報已有較多研究，主要包括利用作物生長模擬模型預(yù)測作物產(chǎn)量的方法［4］，基于遙感技術(shù)的產(chǎn)量動態(tài)預(yù)報方法［5］，利用積分回歸法預(yù)報產(chǎn)量的動態(tài)預(yù)報方法［6］，構(gòu)建溫度、降水、日照及綜合氣候適宜度指數(shù)的產(chǎn)量動態(tài)預(yù)報方法［7-10］，篩選作物不同時段關(guān)鍵氣象影響因子的產(chǎn)量動態(tài)預(yù)報方法［11，12］，以及根據(jù)作物產(chǎn)量氣象影響指數(shù)判斷相似年型的歷史豐歉氣象影響指數(shù)方法［13，14］等。四川省作物產(chǎn)量預(yù)報的傳統(tǒng)方法主要是根據(jù)實地調(diào)研及相關(guān)統(tǒng)計分析進行預(yù)報，雖然準確率較高，但不適用于產(chǎn)量動態(tài)預(yù)報。作物產(chǎn)量預(yù)報動態(tài)化可以提高預(yù)報時效性，及時反映極端氣候事件對作物產(chǎn)量的影響。本研究采用歷史豐歉氣象影響指數(shù)的方法，對四川盆區(qū)油菜產(chǎn)量進行動態(tài)預(yù)報模型構(gòu)建，以期為四川省油菜產(chǎn)量預(yù)報業(yè)務(wù)服務(wù)提供參考。

1 材料與方法

1.1 資料來源

產(chǎn)量資料來源于四川省農(nóng)業(yè)統(tǒng)計年鑒，包括各代表站點1983—2012 年油菜平均單產(chǎn)（kg/hm2）、種植面積（hm2）、總產(chǎn)量（t）數(shù)據(jù)。生育期資料來源于四川盆區(qū)13 個農(nóng)業(yè)氣象觀測站的油菜觀測資料。氣象資料包括各代表站點1982—2012 年逐日平均氣溫（℃）、降水量（mm）和日照時數(shù)（h）。

1.2 代表站點選取

四川盆區(qū)油菜產(chǎn)區(qū)幅員遼闊，地形復(fù)雜，不同氣候類型對油菜生產(chǎn)及產(chǎn)量的影響存在差異。因此，根據(jù)地形及氣候類型差異，在保持縣級行政區(qū)完整的基礎(chǔ)上，將四川盆區(qū)油菜種植區(qū)劃分為5 個區(qū)域，分別為盆周山地產(chǎn)區(qū)、盆東北丘陵產(chǎn)區(qū)、盆中淺丘產(chǎn)區(qū)、盆西平原產(chǎn)區(qū)、盆南丘陵產(chǎn)區(qū)。每個產(chǎn)區(qū)選取5～7 個代表站點，共選取30 個代表站點。盆周山地產(chǎn)區(qū)代表站點分別為蒼溪縣、廣元市中區(qū)、雨城區(qū)、天全縣、洪雅縣和通江縣，盆東北丘陵產(chǎn)區(qū)代表站點分別為巴州區(qū)、儀隴縣、廣安區(qū)、渠縣和達縣，盆中淺丘產(chǎn)區(qū)代表站點分別為梓潼縣、三臺縣、樂至縣、安岳縣、威遠縣和南部縣，盆西平原產(chǎn)區(qū)代表站點分別為江油市、什邡縣、夾江縣、溫江區(qū)、綿竹縣、名山縣和仁壽縣，盆南丘陵產(chǎn)區(qū)代表站點分別為宜賓縣、合江縣、古藺縣、納溪區(qū)、榮縣和富順縣。

1.3 資料處理

1.3.1 產(chǎn)量資料處理 對油菜平均單產(chǎn)進行如下處理：

式（1）中，i 代表1984—2012 年的任一年，i-1為第i 年的前一年，Yi及Yi-1分別為第i 年及第i-1年的油菜平均單產(chǎn)，ΔYi為相鄰兩年油菜平均單產(chǎn)的增減百分率（即產(chǎn)量豐歉氣象影響指數(shù)）。

1.3.2 生育期資料處理 生育期資料采用2010—2012 年3 年的油菜平均發(fā)育期資料。

1.3.3 氣象資料處理 氣象資料的處理主要包括盆區(qū)油菜各產(chǎn)區(qū)的積溫、有效溫度累積、累積降水量、累積日照時數(shù)、分段累積降水量、分段累積日照時數(shù)、標準化降水量、標準化日照時數(shù)、分段標準化降水量和分段標準化日照時數(shù)，具體計算方法參照文獻［15］。

由于引起相鄰兩年作物平均單產(chǎn)變化的主要因素是相鄰兩年氣象條件的差異，因此將相鄰兩年氣象因子的變化表示為：

式（2）中，i 代表1983—2012 年的任一年，i-1為第i 年的前一年，ΔXi為第i 年與第i-1 年的氣象要素差值，Xi與Xi-1分別為第i 年和第i-1 年的積溫、有效溫度累積、累積降水量、累積日照時數(shù)、分段累積降水量、分段累積日照時數(shù)、標準化降水量、標準化日照時數(shù)、分段標準化降水量和分段標準化日照時數(shù)。

1.4 模型構(gòu)建的基本原理及預(yù)報方法

1.4.1 基本原理 影響油菜產(chǎn)量形成的因素主要有氣象條件、土壤條件、管理措施及作物品種類型等。在一定區(qū)域范圍內(nèi)，相鄰兩年的油菜品種及土壤肥力相對較為穩(wěn)定，氣象條件的差異是引起油菜單產(chǎn)變化的主要原因，兩者之間的關(guān)系可表示為：

式（3）中，ΔY 為相鄰兩年油菜單產(chǎn)的增減百分率（即產(chǎn)量豐歉氣象影響指數(shù)），ΔX 為相鄰兩年氣象條件的變化。

歷史豐歉氣象影響指數(shù)模型的構(gòu)建主要取決于歷史相似年型的尋找，本研究將利用氣象綜合診斷指標進行相似年型判定，具體方法如下：

歐式距離：

式（4）至式（6）中，i 為歷史上任意一年，j 為氣象要素序號，k 為預(yù)報年，ΔXij為第i 年同一時段第j個氣象要素差異，N 為樣本長度，ΔXkj為第k 年作物播種至發(fā)布預(yù)報時第j 個氣象要素差異，-Xi與-Xk分別為第i 序列與第k 序列N 個樣本的平均值，dik為第k 年與第i 年的歐式距離，rik為第k 年與第i 年的相關(guān)系數(shù)，Cik為第k 年與第i 年的綜合聚類指標。

Cik越大，則預(yù)報年的氣象條件變化與歷史上某一年越相似。但由于農(nóng)業(yè)技術(shù)水平的不斷提高、油菜種植品種的變化以及預(yù)報發(fā)布之后氣象條件的不斷變化，使得預(yù)報年的油菜產(chǎn)量變化與氣象條件變化最相似年份的油菜產(chǎn)量變化并不一定最接近。因此，為確定預(yù)報年ΔY 的計算方法，在歷史相似年型的選取中，分別選擇溫度、降水和日照3 類氣象因子中Cik最大的3 個歷史相似年型，對3 個歷史相似年型的油菜產(chǎn)量變化與預(yù)報年的油菜實際產(chǎn)量變化進行分析，以確定預(yù)報方法，最終建立油菜產(chǎn)量動態(tài)預(yù)報模型。

1.4.2 預(yù)報方法 分別選取溫度（積溫、有效溫度累積）、降水（累積降水量、分段累積降水量、標準化降水量、分段標準化降水量）和日照（累積日照時數(shù)、分段累積日照時數(shù)、標準化日照時數(shù)、分段標準化日照時數(shù)）3 類氣象要素中的1 個氣象因子，組成32 種組合，通過計算得到每種組合中3 個因子的Cik，選取每個因子歷史上3 個最大的Cik，每種組合對應(yīng)得到9 個ΔY，采用大概率分析法和加權(quán)平均分析法［15］對各產(chǎn)區(qū)預(yù)報結(jié)果與實際產(chǎn)量變化進行對比分析，得到不同時段盆區(qū)油菜各產(chǎn)區(qū)的產(chǎn)量預(yù)報因子及方法，最后根據(jù)各產(chǎn)區(qū)預(yù)報結(jié)果，以面積為權(quán)重進行加權(quán)集成，得到四川盆區(qū)油菜產(chǎn)量的最終預(yù)報結(jié)果。

2 結(jié)果與分析

2.1 四川盆區(qū)油菜產(chǎn)量動態(tài)預(yù)報模型建立

分別對2001—2010 年四川盆區(qū)5 個油菜產(chǎn)區(qū)進行2—4 月的逐月動態(tài)預(yù)報試驗，通過分析預(yù)報結(jié)果，確定了各產(chǎn)區(qū)2 月28 日、3 月31 日及4 月30日油菜產(chǎn)量逐月預(yù)報因子及方法，最后將各產(chǎn)區(qū)預(yù)報結(jié)果按種植面積的權(quán)重進行加權(quán)集成，得到四川盆區(qū)油菜產(chǎn)量動態(tài)預(yù)報方法，結(jié)果見表1。其中，盆周山地產(chǎn)區(qū)的權(quán)重系數(shù)為14.84%，盆東北丘陵產(chǎn)區(qū)的權(quán)重系數(shù)為20.39%，盆中淺丘產(chǎn)區(qū)的權(quán)重系數(shù)為29.77%，盆西平原產(chǎn)區(qū)的權(quán)重系數(shù)為26.72%，盆南丘陵產(chǎn)區(qū)的權(quán)重系數(shù)為8.28%。

表1 油菜產(chǎn)量動態(tài)預(yù)報因子及方法

2.2 動態(tài)預(yù)報模型回代檢驗

利用表1 所選的預(yù)報因子及預(yù)報方法，對2001—2010 年四川盆區(qū)5 個產(chǎn)區(qū)油菜單產(chǎn)豐歉趨勢進行了逐月動態(tài)預(yù)報回代檢驗，預(yù)報結(jié)果的平均正確率見表2。由表2 可以看出，各產(chǎn)區(qū)油菜產(chǎn)量趨勢預(yù)報的平均正確率均在80% 及以上，各產(chǎn)區(qū)集成后的預(yù)報平均正確率均在70% 及以上。

表2 油菜單產(chǎn)豐歉趨勢動態(tài)預(yù)報平均正確率（單位：%）

根據(jù)油菜單產(chǎn)趨勢預(yù)報結(jié)果，對油菜單產(chǎn)進行了定量預(yù)報，預(yù)報平均準確率見表3。由表3 可以看出，除盆南丘陵產(chǎn)區(qū)4 月30 日的預(yù)報平均準確率不足90%，其余均在92.48% 及以上，各產(chǎn)區(qū)集成后的預(yù)報平均準確率均在96.16% 及以上。

表3 油菜單產(chǎn)動態(tài)預(yù)報平均準確率 （單位：%）

2.3 動態(tài)預(yù)報模型預(yù)報檢驗

參照表1 所選預(yù)報因子及預(yù)報方法，對2011年、2012 年四川盆地5 個產(chǎn)區(qū)油菜單產(chǎn)豐歉氣象影響指數(shù)進行了逐月動態(tài)預(yù)報檢驗，預(yù)報結(jié)果見表4。

表4 2011 年、2012 年油菜單產(chǎn)豐歉氣象影響指數(shù)（單位：%）

根據(jù)油菜單產(chǎn)趨勢預(yù)報結(jié)果，對油菜單產(chǎn)進行了定量預(yù)報，預(yù)報準確率見表5。由表5 可以看出，盆區(qū)各產(chǎn)區(qū)油菜單產(chǎn)逐月動態(tài)定量預(yù)報準確率在89.11%～99.39%，其中準確率在98.00% 以上的概率為40.00%，在95.00% 以上的概率為73.33%，在90.00%以上的概率為96.67%。盆區(qū)集成后的預(yù)報準確率在95.84%～100.00%，其中準確率在98.00% 以上的概率為50.00%，在95.00% 以上的概率為100.00%。

表5 2011 年、2012 年油菜單產(chǎn)動態(tài)預(yù)報準確率（單位：%）

3 小結(jié)與討論

本研究方法為四川省油菜產(chǎn)量預(yù)報的客觀性、準確性及動態(tài)性提供了支撐，且整體的預(yù)報準確率較高，有助于推動本地產(chǎn)量預(yù)報業(yè)務(wù)的發(fā)展。但該方法還存在一定的局限性，即未考慮到氣象因子的致災(zāi)量值變化對油菜產(chǎn)量的影響，使預(yù)報結(jié)果的穩(wěn)定性受到影響。每一種產(chǎn)量預(yù)報方法都存在一定的局限性。因此，作物的產(chǎn)量預(yù)報研究可以從作物生長模擬、遙感技術(shù)以及其他數(shù)理統(tǒng)計等諸多方法入手，從多種角度進行預(yù)報，充分發(fā)揮各種模型方法的優(yōu)勢，以進一步提高油菜產(chǎn)量預(yù)報的客觀性、時效性及準確性，為四川省貧區(qū)產(chǎn)量預(yù)報業(yè)務(wù)服務(wù)提供參考。

本研究基于歷史豐歉氣象影響指數(shù)，確定了四川省5 個油菜主產(chǎn)區(qū)及整個盆區(qū)油菜產(chǎn)量的動態(tài)預(yù)報方法，最終建立了油菜產(chǎn)量動態(tài)預(yù)報模型。該模型的回代檢驗結(jié)果表明，盆區(qū)各產(chǎn)區(qū)的動態(tài)趨勢預(yù)報平均正確率在80% 及以上，各產(chǎn)區(qū)集成的平均正確率在70% 及以上；盆地各產(chǎn)區(qū)的動態(tài)定量預(yù)報平均準確率在89.37% 及以上，盆區(qū)集成的平均準確率在96.16% 及以上。預(yù)報檢驗結(jié)果表明，盆地各產(chǎn)區(qū)的動態(tài)定量預(yù)報準確率在89.11% 及以上，盆區(qū)集成的準確率在95.84% 及以上。預(yù)報結(jié)果的準確率較高，可應(yīng)用于實際業(yè)務(wù)服務(wù)。