陳紅梅 劉洪良

摘要 利用合川區(qū)2004～2013年枇杷產(chǎn)量和同期氣象資料，用線性回歸方法對(duì)枇杷產(chǎn)量進(jìn)行趨勢(shì)產(chǎn)量和氣象產(chǎn)量的分離，氣象產(chǎn)量與各生育期43個(gè)氣象要素進(jìn)行相關(guān)分析，利用逐步回歸方法，找出影響枇杷生產(chǎn)的關(guān)鍵氣象因子為3月最低氣溫、10月降水量、1月上旬平均氣溫、1月中旬平均氣溫，建立合川區(qū)枇杷單產(chǎn)的預(yù)測(cè)模型，通過(guò)對(duì)歷史回代、2014年預(yù)測(cè)試驗(yàn)，預(yù)測(cè)效果較好。

關(guān)鍵詞 枇杷；產(chǎn)量；農(nóng)業(yè)氣象預(yù)報(bào)；合川區(qū)

中圖分類(lèi)號(hào) S161 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2014）26-09076-02

枇杷為薔薇科枇杷屬植物，果實(shí)外形美觀，質(zhì)地柔軟多汁，味道鮮美，富含人體所需的多種營(yíng)養(yǎng)元素，其花、葉、樹(shù)皮、根均具有極高的藥用價(jià)值。合川地處亞熱帶季風(fēng)氣候區(qū)，冬暖夏熱，氣候溫和、雨量充沛、日照充足、無(wú)霜期長(zhǎng)、土地酸堿度適宜，是枇杷生長(zhǎng)的最適宜區(qū)之一。近年來(lái)枇杷產(chǎn)業(yè)在合川區(qū)發(fā)展較快，已成為效益較好的農(nóng)業(yè)支柱產(chǎn)業(yè)之一，每年經(jīng)濟(jì)效益超過(guò)2億元。但枇杷的產(chǎn)量直接受氣象條件的影響較大，因此開(kāi)展枇杷產(chǎn)量的農(nóng)業(yè)氣象預(yù)報(bào)研究，能進(jìn)一步提高枇杷生產(chǎn)的氣象服務(wù)水平，為枇杷的生產(chǎn)、儲(chǔ)運(yùn)、銷(xiāo)售、加工和外貿(mào)提供及時(shí)、準(zhǔn)確的服務(wù)信息，對(duì)進(jìn)一步促進(jìn)合川枇杷產(chǎn)業(yè)持續(xù)健康發(fā)展具有重要意義。筆者利用合川區(qū)2004～2013年枇杷產(chǎn)量和同期氣象資料，對(duì)氣象產(chǎn)量與各生育期氣象要素進(jìn)行相關(guān)分析，采用逐步回歸方法，建立合川區(qū)枇杷氣象產(chǎn)量預(yù)報(bào)模型，并通過(guò)歷史回代和預(yù)測(cè)試驗(yàn)進(jìn)行檢驗(yàn)。

1 資料與方法

1.1 資料來(lái)源

2004～2014年氣象資料及枇杷生育期觀測(cè)資料由重慶市合川區(qū)氣象局提供。2004～2014年合川區(qū)枇杷每年的實(shí)際種植面積、實(shí)際總產(chǎn)量等資料由重慶市合川區(qū)統(tǒng)計(jì)局、合川區(qū)農(nóng)委提供。

1.2 氣象產(chǎn)量的分離方法

枇杷產(chǎn)量受品種、土壤、肥力、病蟲(chóng)害、管理措施和氣象條件等因素的綜合制約，隨著社會(huì)的發(fā)展、農(nóng)業(yè)生產(chǎn)水平的提高，從總的趨勢(shì)來(lái)說(shuō)枇杷產(chǎn)量隨年代增加呈現(xiàn)上升趨勢(shì)，但會(huì)因溫度、降水、光照等氣象條件的影響出現(xiàn)較大幅度波動(dòng)。故把枇杷的實(shí)際產(chǎn)量（Y）分解為趨勢(shì)產(chǎn)量（Yt）、氣象產(chǎn)量（Yw）和隨機(jī)噪聲（Ye）3部分，用式子表示為：

Y=Yt+Yw+Ye（1）

其中，趨勢(shì)產(chǎn)量Yt是反映歷史時(shí)期生產(chǎn)力發(fā)展水平的長(zhǎng)周期產(chǎn)量分量，也被稱(chēng)為技術(shù)產(chǎn)量；氣象產(chǎn)量Yw是受氣候要素為主的短周期變化因子（農(nóng)業(yè)氣象災(zāi)害為主）影響的波動(dòng)產(chǎn)量分量；隨機(jī)噪聲Ye則對(duì)產(chǎn)量影響不大，為小量，一般可以忽略不計(jì)^[1-2]。

由于合川枇杷產(chǎn)量資料的年限較短，趨勢(shì)產(chǎn)量處理采取線性回歸方程，然后利用下式（2）從實(shí)際產(chǎn)量中分離出氣象產(chǎn)量：

Yw=Y-Yt（2）

1.3 枇杷產(chǎn)量預(yù)測(cè)模型的建立方法

利用分離的氣象產(chǎn)量與各生育期氣象因子進(jìn)行相關(guān)分析，利用逐步回歸方法，找出影響枇杷生產(chǎn)的關(guān)鍵氣象因子，建立氣象產(chǎn)量預(yù)測(cè)模型。

2 枇杷產(chǎn)量預(yù)測(cè)模型的建立

2.1 趨勢(shì)產(chǎn)量預(yù)測(cè)模型

利用SPSS軟件建立合川區(qū)枇杷的趨勢(shì)產(chǎn)量線性回歸方程：

Yt=37.4T-74 738.151（3）

式中，Yt為趨勢(shì)產(chǎn)量，T為年份。此方程的相關(guān)系數(shù)R=0.994，通過(guò)0.05的顯著檢驗(yàn)。

2.2 氣象產(chǎn)量預(yù)測(cè)方程的建立

2.2.1 預(yù)報(bào)因子的選擇

枇杷的氣象產(chǎn)量不僅與單個(gè)的氣象因子有關(guān)，且是多個(gè)氣象因子共同影響的結(jié)果，根據(jù)枇杷的生物學(xué)特性，通過(guò)分析枇杷生長(zhǎng)發(fā)育對(duì)光、溫、水等主要農(nóng)業(yè)氣象條件的要求，結(jié)合合川區(qū)枇杷各關(guān)鍵生育期的氣候條件，對(duì)氣象產(chǎn)量的豐歉年進(jìn)行了農(nóng)業(yè)氣象條件對(duì)比分析，篩選出與各關(guān)鍵生育期關(guān)系密切的各月或各旬的降水、日平均氣溫、最高氣溫、最低氣溫、日照時(shí)數(shù)、無(wú)霜期等43個(gè)相關(guān)氣象因子，作為預(yù)報(bào)氣象產(chǎn)量的自變量因子（表1）。

2.2.2 預(yù)測(cè)方程的建立。

利用SPSS 17.0數(shù)據(jù)統(tǒng)計(jì)軟件進(jìn)行逐步回歸分析，對(duì)相關(guān)氣象因子進(jìn)行篩選，建立多元線性回歸分析模型。

通過(guò)相關(guān)分析（圖1），氣象產(chǎn)量Yw和3月最低氣溫、10月降水量、1月上旬平均氣溫、1月中旬平均氣溫、年降水量、4月降水量和7月降水量等7個(gè)因子相關(guān)性達(dá)到顯著或極顯著，通過(guò)逐步回歸分析，建立氣象產(chǎn)量預(yù)測(cè)模型：

Yw=78.859+51.229X₁-1.586X₂-18.273X₃-6.878X₄

2.2.4 合川區(qū)2014年枇杷單產(chǎn)預(yù)測(cè)檢驗(yàn)。

將2014年枇杷生育期相關(guān)氣象資料帶入式（3）、式（4）中，得出2014年枇杷單產(chǎn)為7 452.75 kg/hm2，2014年實(shí)際產(chǎn)量為7 200.00 kg/hm2，相對(duì)誤差為4%，效果較好。

3 結(jié)論

利用氣象產(chǎn)量與枇杷生育期間的相關(guān)顯著的氣象因子，通過(guò)逐步回歸，建立合川區(qū)枇杷氣象產(chǎn)量預(yù)報(bào)模型，通過(guò)2004～2013年歷史回代檢驗(yàn)、2014年預(yù)測(cè)試驗(yàn)，表明應(yīng)用該預(yù)報(bào)模型對(duì)合川區(qū)枇杷單產(chǎn)進(jìn)行預(yù)測(cè)是可行。

由于資料年代有限，該預(yù)報(bào)模型可能存在一定缺陷，隨著資料的不斷積累，應(yīng)逐步得到完善。

參考文獻(xiàn)

[1] 楊霏云，王建林.晚稻單產(chǎn)動(dòng)態(tài)預(yù)測(cè)方法研究[J].氣象科技，2005，33（5）：433-436.

[2] 宋迎波，王建林，陳暉，等.中國(guó)油菜產(chǎn)量動(dòng)態(tài)預(yù)報(bào)方法研究[J].氣象，2008（3）：93-99.