黃 彭，郝 妙，杜永華，陳一輝

（1宜賓學院生命科學與食品工程學院，四川宜賓644000；2宜賓學院香料植物資源開發(fā)與利用四川省高校重點實驗室，四川宜賓644000；3中共宜賓市委黨校，四川宜賓644000；4屏山縣農(nóng)業(yè)和林業(yè)局，四川宜賓644000）

基于GM(1,1)模型的四川糧食產(chǎn)量預(yù)測研究

黃 彭1,2，郝 妙3，杜永華2，陳一輝4

（1宜賓學院生命科學與食品工程學院，四川宜賓644000；2宜賓學院香料植物資源開發(fā)與利用四川省高校重點實驗室，四川宜賓644000；3中共宜賓市委黨校，四川宜賓644000；4屏山縣農(nóng)業(yè)和林業(yè)局，四川宜賓644000）

中國人口眾多，糧食安全關(guān)系到國計民生，加強糧食產(chǎn)量預(yù)測有利于確保糧食安全。根據(jù)2001—2015年四川糧食產(chǎn)量的歷史數(shù)據(jù)，運用灰色系統(tǒng)理論，建立基于弱化緩沖算子的GM(1,1)預(yù)測模型，通過殘差、級比偏差、關(guān)聯(lián)度、后驗差檢測、模擬數(shù)據(jù)檢查對模型的合理性和精度進行誤差檢驗，并應(yīng)用模型預(yù)測未來3年的糧食產(chǎn)量。研究結(jié)果表明，灰色系統(tǒng)理論GM(1,1)適用于糧食產(chǎn)量預(yù)測且具有較高的精度。預(yù)測了2016、2017、2018年的糧食產(chǎn)量同比增長分別為-2.11%、-0.39%和1.21%，由此得出未來糧食產(chǎn)量將在波動中增長。

灰色理論系統(tǒng)；糧食產(chǎn)量；預(yù)測；四川

0 引言

國以民為本，民以食為天，糧食是安天下穩(wěn)民心的戰(zhàn)略物質(zhì)，是經(jīng)濟社會發(fā)展的重要物質(zhì)基礎(chǔ)[1]。糧食生產(chǎn)是關(guān)系到國計民生的頭等大事，準確預(yù)測糧食產(chǎn)量有利于優(yōu)化配置農(nóng)業(yè)資源，合理調(diào)整農(nóng)業(yè)結(jié)構(gòu)，促進農(nóng)業(yè)經(jīng)濟健康發(fā)展[2]。除無法確定的因素外，糧食產(chǎn)量預(yù)測還受自然環(huán)境、氣候條件、社會經(jīng)濟3類因素的影響，是復(fù)雜的農(nóng)學和統(tǒng)計學相互交叉的問題[3]。

目前應(yīng)用于農(nóng)業(yè)經(jīng)濟管理領(lǐng)域的預(yù)測方法有：神經(jīng)網(wǎng)絡(luò)預(yù)測法、遙感預(yù)測模型、氣象生產(chǎn)力模型預(yù)測法、灰色預(yù)測模型、多元回歸模型、時間序列模型等[4-5]。與其他預(yù)測方法相比，灰色系統(tǒng)理論無論在理論上，還是在方法以及實際應(yīng)用上，均具有一定的優(yōu)勢[6]?；谌趸彌_算子的灰色系統(tǒng)理論GM(1,1)模型通過弱化數(shù)據(jù)序列的隨機性，充分挖掘系統(tǒng)演化規(guī)律，使系統(tǒng)所需原始信息量減少，簡化計算過程，從而提高預(yù)測精度[7-8]。目前，該模型廣泛地運用于價格和產(chǎn)量預(yù)測中。程偉、黃學文等[9-10]學者證實了灰色系統(tǒng)理論GM(1,1)模型用于預(yù)測農(nóng)林產(chǎn)品的產(chǎn)量的可行性。應(yīng)用GM(1,1)模型預(yù)測黃河三角洲地區(qū)糧食產(chǎn)量，關(guān)聯(lián)度、均方差比值及小誤差概率均判定該模型的精度為一級，誤差較小，預(yù)測結(jié)果可靠[11]。四川作為糧食生產(chǎn)大省，但以往文獻對其糧食產(chǎn)量進行預(yù)測的研究不多。本研究根據(jù)四川糧食產(chǎn)量歷史數(shù)據(jù)，采用灰色系統(tǒng)理論GM(1,1)的理論和方法建立預(yù)測模型，并對擬合模型的合理性和可靠行進行檢驗分析，進而為糧食產(chǎn)量預(yù)測提供科學的理論依據(jù)。

1 模型原理

1982年，鄧聚龍[12]首次提出灰色系統(tǒng)理論（Grey Model，簡稱GM），1985年成立灰色系統(tǒng)研究會。灰色系統(tǒng)是指系統(tǒng)中既含有已知信息，又含有不確定信息，以“部分信息已知，部分信息未知”的“小樣本、貧信息”不確定性系統(tǒng)為研究對象，對部分已知信息的生成開發(fā)，提取有價值的信息，實現(xiàn)對系統(tǒng)運行行為、演化規(guī)律的正確描述和有效控制，具有所需原始信息量少、計算簡單、通用性強及預(yù)測精度高等優(yōu)點[13-14]?；疑碚撓到y(tǒng)(GM)對試驗觀察數(shù)據(jù)沒有特殊的要求和限制，在建模過程中將具有隨機性的原始數(shù)據(jù)經(jīng)過科學的處理，弱化雜亂無章的原始數(shù)據(jù)，轉(zhuǎn)化為具有較強規(guī)律性的數(shù)據(jù)，從而弱化不確定性和強化規(guī)律性。經(jīng)過超過30年的發(fā)展，在自然科學、經(jīng)濟社會領(lǐng)域里，解決生產(chǎn)、生活和科研等實際問題上GM模型分析預(yù)測成果顯著[15]。

2 模型的建立

近年來，四川糧食產(chǎn)量基本維持在2000萬t以上，糧食預(yù)測是政府制定和實施農(nóng)業(yè)經(jīng)濟政策的重要依據(jù)。研究選取了2001—2012年的糧食單產(chǎn)量作為輸出因子，2013—2015年的糧食產(chǎn)量數(shù)據(jù)作為模型驗證數(shù)據(jù)。從圖1可知，2005年和2006年分別為近年來糧食產(chǎn)量最高和最低的兩個年份。2004年底，國家、省市連續(xù)出臺了一系列優(yōu)惠政策，積極推進糧食品種改良、提高植保水平，調(diào)整農(nóng)業(yè)結(jié)構(gòu)，大力促進糧食的恢復(fù)與生產(chǎn)，使糧食供求偏緊的狀況很快得到緩解，2005年四川省糧食豐收并創(chuàng)造歷史最高為3409.2萬t[16]。然而，2006年川渝地區(qū)受到特大旱災(zāi)，導(dǎo)致大春減產(chǎn)，年產(chǎn)量達到歷年來最低產(chǎn)量為2859.8萬t，同比下降16.11%。2006年后糧食產(chǎn)量不斷增加，但同比增長率不斷減少，波動范圍越來越小。

由圖1得到糧食產(chǎn)量的原始序列：X0=(3056.54,3275.20,3183.30,3326.53,3409.19,2859.85,3026.92,3140.31,3194.68,3223.52,3292.25,3315.70)，根 據(jù) 可 行性檢驗，因為覆蓋范圍不全在(e-2/(12+1),e2/(12+1))即(0.9260,1.1663)區(qū)間內(nèi)，對原始序列做構(gòu)造弱化緩沖算子處理：Xd={Xd(1),Xd(2),…,Xd(n)}，其中，經(jīng)過弱化緩沖算子作用，生成新序列為：Xd={3192.00,3204.31,3197.22,3198.77,3182.80,3150.46,3198.89,3233.29,3256.54,3277.15,3303.97,3315.70}。

圖1 2001—2015年四川糧食產(chǎn)量波動走勢

對Xd作緊鄰均值生成：Zt={Zt(2),Zt(3),…Zt(n)}；其中，生成緊鄰 均 值Zt＝(4794.15,7994.92,11192.92,14383.70,17550.34,20725.01,23941.10,27186.02,30452.86,33743.43,37053.26)，計算GM(1,1)參數(shù)a＝-0.0041，b＝3143.5168。確定GM(1,1)的擬合模型為：Xd(k)-0.0041Zt(k)＝3143.5168。

因此，GM(1,1)的白化響應(yīng)式為[17]：可 求 得Xd(k+ 1) =769903.41e0.0041k-766711.41，Xd(k+1)=X?d(k+1)-X?d(k),(k=1,2,3,…,n)。

3 結(jié)果與分析

3.1 誤差檢驗

誤差檢驗結(jié)果見表1[18]。

3.1.3 關(guān)聯(lián)度檢驗 設(shè) min=min{|δk|}，max=max{|δk|}；則min=1.3551，max=64.9306；設(shè)則η(k)=(0.5696,0.4071,0.7894,0.1788,0.3672,1.3276,1.1144,1.0000,0.9046,0.7958,0.7290)；關(guān)聯(lián)度；根據(jù)經(jīng)驗判斷G(1,1)模型較為滿意[19]。

表1 誤差檢驗表

3.1.5 模型實際檢測 通過模擬方程Xd(k)-0.0041Zt(k)＝3143.5168和白化響應(yīng)方差Xd(k+1)=769903.41e0.0041k，Xd(k+1)=X?d(k+1)-X?d(k) (k=1,2,3,…,n)，計算出2013、2014和2015年四川糧食產(chǎn)量的模擬數(shù)據(jù)，與實際數(shù)據(jù)對比分析，誤差均小于5%，說明模擬出來的數(shù)據(jù)正確性達到95%以上，對實際數(shù)據(jù)與模擬數(shù)據(jù)進行F檢驗，，說明2013—2015年四川糧食產(chǎn)量模擬數(shù)據(jù)與實際數(shù)據(jù)方差差異不顯著，即模擬數(shù)據(jù)可靠，檢驗結(jié)果見表2[21-22]。

從上面2.1各項檢測可知，擬合模型GM(1,1)的級比偏差檢驗、殘差檢驗、關(guān)聯(lián)度檢驗、后驗差檢驗和模擬實際數(shù)據(jù)檢測都能通過，說明模型有足夠的精度滿足要求，可作為預(yù)測模型對四川糧食產(chǎn)量進行預(yù)測。

3.2 預(yù)測應(yīng)用

基于擬合模型GM(1,1)式Xd(k)-0.0041Zt(k)＝3143.5168和白化響應(yīng)方差，白化響應(yīng)式Xd(k+1)=769903.41e0.0041k-766711.41，Xd(k+1)=X?d(k+1)-X?d(k)(k=1,2,3,…,n)和緊鄰生成值得到2016—2018年四川糧食產(chǎn)量預(yù)測值，分別為3370.28，3357.16，3397.92（見表3），即預(yù)測2016年糧食產(chǎn)量為3370.28萬t，2017年糧食產(chǎn)量為3357.16萬t，2018年糧食產(chǎn)量為3397.92萬t。由此可見，在未來3年內(nèi)糧食產(chǎn)量仍存在一定的波動，同比增長分別為-2.11%、-0.39%和1.21%，未來四川糧食產(chǎn)量將在波動中增長。

表3 2016—2018年糧食產(chǎn)量預(yù)測

4 結(jié)論

通過對2001—2015年四川糧食產(chǎn)量歷史數(shù)據(jù)的統(tǒng)計分析和基于弱化緩沖算子的GM(1,1)建模預(yù)測未來3年的糧食產(chǎn)量，結(jié)果如下。

（1）通過對歷史數(shù)據(jù)緩存弱化處理，得到GM(1,1)擬合模型Xd(k)-0.0041Zt(k)＝3143.5168和白化響應(yīng)式 ：Xd(k+ 1) =769903.41e0.0041k- 766711.41，Xd(k+1)=X?d(k+1)-X?d(k)(k=1,2,3,…,n)。

（2）誤差檢驗糧食產(chǎn)量灰色系統(tǒng)理論GM(1,1)，P0=99.30%＞90%、ρ(avg)=0.8783%＜10%、r=0.7439＞0.6和P=0.8182，此方法具有較高的預(yù)測精度，不需要殘差修正；模型實際檢測，F(xiàn)=3.064＜F0.05=19.00，與實際數(shù)據(jù)對比分析可知，誤差均小于5%，模擬出來的數(shù)據(jù)正確性達到95%以上。

（3）預(yù)測2016、2017、2018年的糧食產(chǎn)量分別為3370.28、3357.16、3397.92 萬 t。同比增長分別為-2.11%、-0.39%和1.21%。

5 討論

灰色理論GM(1,1)模型是通過對原始數(shù)據(jù)整理變換來尋求變化趨勢。雖然GM(1,N)和GM(r,h)等預(yù)測結(jié)果準確度更高，但本研究從理論和模型檢測已論證用GM(1,1)模型預(yù)測糧食產(chǎn)量是可行的，不需要進行更復(fù)雜的多個變量數(shù)。2005年四川糧食豐產(chǎn)受到2004年底國家和地方優(yōu)惠政策，農(nóng)業(yè)結(jié)構(gòu)調(diào)整，農(nóng)民種糧積極性提高，產(chǎn)量達到歷史最高值；然而2006年四川省糧食產(chǎn)量在近15年來最低，是由于川渝地區(qū)受自然災(zāi)難（特大伏旱）的影響所致。這說明糧食產(chǎn)量受政策扶持、糧農(nóng)主觀能動性、自然災(zāi)難、金融危機、人口流動、科技投入水平、耕地面積和社會經(jīng)濟發(fā)展水平等因素的影響，則使預(yù)測結(jié)果的準確性將會產(chǎn)生較大偏差[23]。

從大量研究數(shù)據(jù)分析可知，每當糧食產(chǎn)量連續(xù)增長幾年后，由于受內(nèi)外界各因素的影響，增長率將降到一個相對較低的水平，并在此水平上持續(xù)徘徊幾年，從而形成了糧食產(chǎn)量具有波動性和周期性的規(guī)律特征，這與本研究的模型預(yù)測相符合[24]。未來中央和地方政策、生產(chǎn)環(huán)境的變化對四川地區(qū)糧食產(chǎn)量的影響，還需根據(jù)實際情況即時對以上GM(1,1)模型進一步完善修正，使擬合數(shù)據(jù)與實際值趨向基本一致[25]。

6 討論

習近平總書記提出，只要糧食不出大問題，中國的事就穩(wěn)得住。糧食安全既是經(jīng)濟問題，也是政治問題，是國家發(fā)展穩(wěn)定的基本前提。尤其是應(yīng)對當前風云變幻的國際形勢，糧食安全更不能麻痹松懈。從預(yù)測結(jié)果來看，未來四川省的糧食產(chǎn)量會在波動中增長。因此，為確保糧食生產(chǎn)安全和產(chǎn)量增長，四川省應(yīng)切實增強新形勢下的糧食安全意識；從堅守耕地紅線、建設(shè)高標準農(nóng)田、提高糧食生產(chǎn)科技水平等方面著手，鞏固和提高糧食生產(chǎn)能力；落實糧食扶持生產(chǎn)政策、抓好糧食收購、提高種糧比較收益等工作，保護糧農(nóng)的種糧積極性，提高糧農(nóng)種植收入。

[1]張元紅,劉長全,國魯來.中國糧食安全狀況評價與戰(zhàn)略思考[J].中國農(nóng)村觀察,2015(1):2-14.

[2]楊克磊,張振宇,和美.應(yīng)用灰色GM(1,1)模型的糧食產(chǎn)量預(yù)測研究[J].重慶理工大學學報:自然科學版,2015(4):124-127.

[3]蘇博,劉魯,楊方廷.GM(1,N)灰色系統(tǒng)與BP神經(jīng)網(wǎng)絡(luò)方法的糧食產(chǎn)量預(yù)測比較研究[J].中國農(nóng)業(yè)大學學報,2006(4):99-104.

[4]龔波,肖國安,張四梅.基于灰色系統(tǒng)理論的湖南糧食產(chǎn)量預(yù)測研究[J].湖南科技大學學報:社會科學版,2012,15(5):62-65,79.

[5]史衛(wèi)亞,何明揚.基于智能信息的糧食產(chǎn)量趨勢預(yù)測模型[J].信息通信,2016(09):151-152.

[6]王成智.灰色系統(tǒng)模型對水產(chǎn)品產(chǎn)量的實證分析[D].大連:大連海洋大學.

[7]亓雪龍,孫洪雁,陳進展.基于灰色系統(tǒng)理論的山東省水果產(chǎn)量預(yù)測模型[J].山東農(nóng)業(yè)科學,2010(3):16-19.

[8]姚宏軍.基于灰色系統(tǒng)理論的黃/紅麻產(chǎn)量預(yù)測研究[D].福州:福建農(nóng)林大學,2014.

[9]黃文學,劉凌,閆爭亮,等.基于灰色系統(tǒng)理論的我國木質(zhì)產(chǎn)品產(chǎn)量預(yù)測模型的建立[J].東北林業(yè)大學學報,2012(10):156-158.

[10]Yuan L I,Bo J.Study on the Yield Prediction Model of Processing Tomato Based on the Grey System Theory[J].Agricultural Science&Technology,2011,33(5):632-633,642.

[11]李靜,朱農(nóng),李鳳桂,等.近10年黃河三角洲地區(qū)糧食產(chǎn)量及灰色預(yù)測[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(5):15-19.

[12]李智鵬,李江宏.基于灰色系統(tǒng)理論的新疆水果產(chǎn)量預(yù)測[J].中國科技投資,2013(9):92-126.

[13]程偉,劉國璧.灰色系統(tǒng)理論在糧食產(chǎn)量預(yù)測中的應(yīng)用[J].北京電子科技學院學報,2008,16(2):62-64.

[14]Sifeng L,Julong D.GM(1,1)Coding for Exponential Series[J].The Joural of Grey,1999,2:142.

[15]唐子喬.湖南省生豬生產(chǎn)預(yù)測預(yù)警研究與實證分析[D].長沙:湖南師范大學,2014.

[16]魯曉東.2005年糧食產(chǎn)銷形勢分析[J].中國農(nóng)村小康科技,2005(1):45-46.

[17]彭道民,王飛.浙江省海洋捕撈產(chǎn)量及漁船數(shù)量動態(tài)預(yù)測與分析——基于灰色系統(tǒng)理論[J].安徽農(nóng)業(yè)科學,2015(18):259-261.

[18]亓雪龍,孫洪雁,陳進展.基于灰色系統(tǒng)理論的山東省水果產(chǎn)量預(yù)測模型[J].山東農(nóng)業(yè)科學,2010(3):16-19.

[19]郝妙,傅新紅,陳蓉.灰色系統(tǒng)理論在生豬價格預(yù)測中的應(yīng)用[J].中國農(nóng)學通報,2014(14):310-314.

[20]徐興梅,劉茹,常政,等.基于灰色系統(tǒng)與時序分析方法的吉林省糧食產(chǎn)量預(yù)測[J].農(nóng)業(yè)機械,2015(23):93-96.

[21]史亞楠.我國糧食產(chǎn)量影響因素的實證研究[D].荊州:長江大學,2014.

[22]宰松梅,溫季,仵峰,等.基于灰色關(guān)聯(lián)分析的遼寧省糧食產(chǎn)量預(yù)測模型[J].節(jié)水灌溉,2011(05):64-66.

[23]呂書文,楊婧,李榮.灰色系統(tǒng)理論在瀾滄縣松脂產(chǎn)量預(yù)測中的應(yīng)用[J].綠色科技,2015(2):186-187.

[24]苗森.黑龍江省四種主要糧食作物單產(chǎn)的分析與預(yù)測[D].哈爾濱:東北林業(yè)大學,2013.

[25]馬天祥.黑龍江省糧食生產(chǎn)影響因素研究[D].哈爾濱:東北農(nóng)業(yè)大學,2014.

Prediction of Grain Yield of Sichuan Based on GM(1,1)Model

Huang Peng1,2,Hao Miao3,Du Yonghua2,Chen Yihui4
(1College of Life Science and Food Engineering,Yibin University,Yibin 644000,Sichuan,China;2Key Lab of Aromatic Plant Resources Exploitation and Utilization in Sichuan Higher Education,Yibin University,Yibin 644000,Sichuan,China;3Party School of Yibin CPC,Yibin 644000,Sichuan,China;4Agricultural and Forestry Bureau of Pingshan County,Yibin 644000,Sichuan,China)

China has a large population,and food security is a major event relating to people’s livelihood.Strengthening grain yield prediction can ensure food security.Based on historical data of grain yield in Sichuan from 2001 to 2015,we adopted grey system theory and set up GM(1,1)prediction model on the basis of the weakening buffer operator,then used four error checking methods,including residual error,grade ratio deviation,correlation degree and posteriori error to test the reasonability and the prediction accuracy.Then we applied GM(1,1)prediction model to predicting the grain yield in future three years.Test results indicate that the grey system theory is appropriate when it is applied to grain yield prediction and has high prediction accuracy.According to GM(1,1)prediction model,the grain yield of the year 2016,2017,2018 will increase by-2.11%,-0.39%and 1.21%respectively year on year.It can be concluded that the grain yield of Sichuan will increase in fluctuation in the near future.

Grey System Theory;Grain Yield;Prediction;Sichuan

S114

A論文編號：cjas16120010

四川省農(nóng)村發(fā)展研究中心項目(CR1622)；固態(tài)發(fā)酵資源利用四川省重點實驗室項目(2016GTJ008)；宜賓學院科研項目(2016PY02)。

黃彭，男，1989年出生，四川儀隴人，助教，碩士，研究方向為食品安全與技術(shù)。通信地址：644000四川省宜賓市酒圣路8號宜賓學院生命科學與食品工程學院，Tel：0831-3545069，E-mail：pehuang@sina.com。

郝妙，女，1988年出生，四川漢源人，講師，碩士，研究方向：農(nóng)林經(jīng)濟理論與政策。通信地址：644000四川省宜賓市翠屏區(qū)蜀南大道中段69號中共宜賓市委黨校，Tel：0831-2379878，E-mail：937402414@qq.com。

2016-12-06，

2017-08-19。