我國牧草產(chǎn)品進(jìn)口價格及進(jìn)口量預(yù)測

董寧+趙偉

[提要] 目前，國內(nèi)牧草種植的產(chǎn)量和質(zhì)量均不能滿足國內(nèi)畜牧業(yè)需求，進(jìn)口優(yōu)質(zhì)牧草成為解決需求的一個重要途徑。而進(jìn)口價格和進(jìn)口量變動會影響國內(nèi)進(jìn)口企業(yè)以及使用企業(yè)相關(guān)決策，影響國內(nèi)牧草產(chǎn)業(yè)體系的建立。本文運用2008年1月至2016年12月苜蓿草進(jìn)口價格以及進(jìn)口量數(shù)據(jù)，使用ARIMA模型預(yù)測苜蓿草進(jìn)口價格及進(jìn)口量。

關(guān)鍵詞：苜蓿草進(jìn)口；ARIMA模型；預(yù)測

基金項目：山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系牧草產(chǎn)業(yè)經(jīng)濟(jì)崗位專家項目（SDAIT-23-06）；“雙一流”獎補(bǔ)資金（SYL2017XTTD16）；通信作者：趙偉

中圖分類號：F326.3 文獻(xiàn)標(biāo)識碼：A

收錄日期：2017年12月19日

一、文獻(xiàn)回顧

當(dāng)前，國內(nèi)人們的食物來源結(jié)構(gòu)正在改變，對于草食型產(chǎn)物如牛羊肉、奶蛋等產(chǎn)品的需求日益增加，這就增加了人們對于草食畜產(chǎn)品食品安全的關(guān)注度，由于眾多食品安全事件人們逐漸認(rèn)識到優(yōu)質(zhì)安全的牧草才能生產(chǎn)出優(yōu)質(zhì)安全的草食畜產(chǎn)品。國內(nèi)已經(jīng)逐步建立了使用苜蓿等優(yōu)質(zhì)牧草的意識，許多大型養(yǎng)殖企業(yè)為提高畜產(chǎn)品數(shù)量和質(zhì)量開始注重苜蓿等優(yōu)質(zhì)牧草的使用。但是我國種植苜蓿草的產(chǎn)量和質(zhì)量均不能滿足國內(nèi)畜牧業(yè)巨大的需求，進(jìn)口優(yōu)質(zhì)苜蓿草成為解決需求的一個重要途徑。如目前進(jìn)口苜蓿草產(chǎn)品使用量前十的企業(yè)有現(xiàn)代牧業(yè)、上海奶牛集團(tuán)、東營澳亞、蒙牛富源牧業(yè)、福建長富乳業(yè)、上海光明荷斯坦牧業(yè)、南京衛(wèi)崗乳業(yè)、新希望乳業(yè)以及天津嘉立牧業(yè)等，2014年這十家企業(yè)的使用份額約占苜蓿草進(jìn)口量的48.8%。

國內(nèi)種植苜蓿的經(jīng)濟(jì)效益自2014年有所下降，雖然較其他牧草產(chǎn)品純收益有優(yōu)勢但是相較糧食作物優(yōu)勢不明顯，對國內(nèi)草產(chǎn)品的穩(wěn)定性和積極性產(chǎn)生負(fù)面影響。根據(jù)乳品的收購價格，企業(yè)會將國產(chǎn)苜蓿草的價格控制在1，800～2，000元/噸（如山東、陜西等省份），部分省份低至600元/噸左右（如河南省）。加上國內(nèi)運輸費用使得國產(chǎn)苜蓿價格高于進(jìn)口價格，致使國內(nèi)的苜蓿草使用企業(yè)對進(jìn)口苜蓿草依賴度過高，如現(xiàn)代牧業(yè)2014年使用進(jìn)口苜蓿草14萬噸，但使用國產(chǎn)苜蓿草只有3萬噸。因此，進(jìn)口價格的變動將直接關(guān)系著國內(nèi)養(yǎng)殖戶的收入，關(guān)系著國內(nèi)牧草種植業(yè)、畜產(chǎn)品養(yǎng)殖業(yè)、畜產(chǎn)品加工業(yè)等行業(yè)的穩(wěn)定，并最終關(guān)系著我國國民生計。

二、我國苜蓿產(chǎn)品貿(mào)易現(xiàn)狀

1992～2008年我國苜蓿草產(chǎn)品出口大于進(jìn)口，我國苜蓿產(chǎn)業(yè)在這期間發(fā)展速度較快，主要得益于國際市場對苜蓿產(chǎn)品的強(qiáng)勁需求。當(dāng)時大型苜蓿產(chǎn)品加工企業(yè)，其目標(biāo)市場主要瞄準(zhǔn)日本、韓國和東南亞國家。這樣，苜蓿的價格高低主要決定于這些國家經(jīng)濟(jì)狀況和對苜蓿產(chǎn)品需求量，同樣也決定于這些國家的貿(mào)易政策。由于這些國家經(jīng)濟(jì)波動較為頻繁，且對我國農(nóng)產(chǎn)品出口經(jīng)常設(shè)置各種貿(mào)易壁壘，造成國內(nèi)苜蓿產(chǎn)品價格的劇烈波動，這就挫傷了農(nóng)戶種植苜蓿的積極性。由于2004年后糧食直補(bǔ)政策，農(nóng)戶種植糧食面積擠壓了苜蓿種植面積，造成了我國目前苜蓿基本依賴進(jìn)口的局面。特別是2009年以來，由于國內(nèi)苜蓿需求量不斷增加，苜蓿產(chǎn)品進(jìn)口量大幅增加。尤其在1993～2000年，出口明顯大于進(jìn)口，苜蓿草產(chǎn)品出口與進(jìn)口之差年平均為24.60萬噸；2009年和2010年苜蓿草產(chǎn)品進(jìn)口與出口之差分別為5.34萬噸和21.2萬噸。

此前，我國許多地方的苜蓿生產(chǎn)以出口日本、韓國和中國臺灣地區(qū)為主要目標(biāo)，外向型苜蓿生產(chǎn)為主要方向，實踐證明出口貿(mào)易存在許多困難，這些地區(qū)對所進(jìn)口的苜蓿草產(chǎn)品要求質(zhì)量高、批量大、貨源穩(wěn)定，而我國苜蓿加工企業(yè)所生產(chǎn)的產(chǎn)品基本上達(dá)不到要求。我國苜蓿草產(chǎn)品貿(mào)易格局逐步由出口轉(zhuǎn)向進(jìn)口為主。1997～2007年美國出口到中國的牧草量非常小，相比其牧草總出口量可以忽略不計，但2008年后，美國出口到中國的牧草呈井噴式上升趨勢。從2007年到2012年的5年間苜蓿干草捆進(jìn)口量增長了100倍，引起了美國專家、草商極大的重視和興趣，中國苜蓿市場和產(chǎn)業(yè)的興起與美國的苜蓿產(chǎn)品進(jìn)口中國同步進(jìn)行，是一個雙贏的局面。但是目前，國內(nèi)進(jìn)口苜蓿已經(jīng)占到了全國苜蓿總銷售量的40%以上，對國外尤其是美國苜蓿的依賴性較高。

三、數(shù)據(jù)來源及研究方法

（一）數(shù)據(jù)來源。考慮數(shù)據(jù)的可獲得性，本文以2008年1月至2016年12月的苜蓿草進(jìn)口價格以及進(jìn)口量月度數(shù)據(jù)作為研究對象。數(shù)據(jù)來源于荷斯坦網(wǎng)站（http：//www.hesitan.com/），數(shù)據(jù)序列相對完整。（表1）

（二）研究方法。ARIMA模型，即自回歸積分滑動平均模型，就是指將非平穩(wěn)時間序列轉(zhuǎn)化為平穩(wěn)時間序列，然后將因變量僅對它的滯后值以及隨機(jī)誤差項的現(xiàn)值和滯后值進(jìn)行回歸所建立的模型。其中，ARIMA（p，d，q）稱為差分自回歸移動平均模型，AR是自回歸，p為自回歸項；MA為移動平均，q為移動平均項數(shù)，d為時間序列成為平穩(wěn)時所做的差分次數(shù)。ARIMA模型根據(jù)原序列是否平穩(wěn)以及回歸中所含部分的不同，包括移動平均過程（MA）、自回歸過程（AR）、自回歸移動平均過程（ARMA）以及ARIMA過程。

根據(jù)模型參數(shù)選擇的原則，參數(shù)d表示為使時間序列平穩(wěn)的差分階數(shù)。確定p和q值具體來說主要有兩種方法：一是通過數(shù)據(jù)的統(tǒng)計特征，自相關(guān)和偏相關(guān)函數(shù)來判斷，如表2所示；二是根據(jù)信息準(zhǔn)則來判斷。通常采用AIC和BIC準(zhǔn)則，AIC準(zhǔn)則又稱為最小信息準(zhǔn)則，它是由日本學(xué)者赤池提出應(yīng)用于AR模型的分析定階中，BIC準(zhǔn)則比AIC準(zhǔn)則增加了更嚴(yán)格的懲罰條件。理論上，AIC和BIC要盡可能的小，當(dāng)模型擬合優(yōu)度上升時，AIC和BIC的值會趨近于負(fù)無窮，即AIC和BIC值模型1小于模型2，則稱模型1優(yōu)于模型2。（表2）

四、實證分析

（一）進(jìn)口量預(yù)測。繪制出原始數(shù)據(jù)苜蓿草進(jìn)口量的時間序列，可以初步判定苜蓿草進(jìn)口量序列是非平穩(wěn)的。進(jìn)一步對其進(jìn)行ADF單位根檢驗，表明苜蓿草進(jìn)口量是非平穩(wěn)的。為了對序列進(jìn)行分析，需要使其平穩(wěn)化，進(jìn)行差分處理結(jié)果表明一階差分后苜蓿草進(jìn)口量序列是平穩(wěn)的。（表3）

分析一階差分后的苜蓿草進(jìn)口量的自相關(guān)偏相關(guān)函數(shù)，如圖1所示，結(jié)合表3和AIC準(zhǔn)則根據(jù)AIC和SC取值最小準(zhǔn)則，多次試驗對比AIC和SC的值選出最適合的模型，最終建立模型ARIMA（3，1，3）。（表4）

利用Eviews做出ARIMA（3，1，3）殘差圖，如圖1為殘差、實際值和擬合值的時序圖，大體可以看出殘差圍繞0值上下波動，進(jìn)一步對殘差做平穩(wěn)性檢驗，結(jié)果表明殘差通過了適應(yīng)性檢驗，說明模型很好的提取了信息。（圖1、表5）

根據(jù)ARIMA（3，1，3）模型對苜蓿草進(jìn)口量進(jìn)行動態(tài)預(yù)測，預(yù)測期間為12個月，預(yù)測結(jié)果如表5。進(jìn)口量在經(jīng)過2016年底的波谷后開始上升，2017年進(jìn)口量持續(xù)上升，總之我國對進(jìn)口苜蓿草的依賴在加重。（表6）

（二）進(jìn)口價格預(yù)測。繪制苜蓿草進(jìn)口量的時間序列，從表7可以初步判定苜蓿草進(jìn)口價格序列是非平穩(wěn)的。進(jìn)一步對其進(jìn)行ADF單位根檢驗，表明苜蓿草進(jìn)口價格是非平穩(wěn)的。為了對序列進(jìn)行分析，需要使其平穩(wěn)化，進(jìn)行差分處理結(jié)果表明一階差分后苜蓿草進(jìn)口價格序列是平穩(wěn)的。（表7）

分析一階差分后的苜蓿草進(jìn)口價格的自相關(guān)偏相關(guān)函數(shù)，結(jié)合表3和AIC準(zhǔn)則根據(jù)AIC和SC取值最小準(zhǔn)則，多次試驗對比AIC和SC的值選出最適合的模型，最終建立模型ARIMA（2，1，2），如表8所示。（表8）

利用Eviews做出ARIMA（2，1，2）殘差圖，如圖2為殘差、實際值和擬合值的時序圖，大體可以看出殘差圍繞0值上下波動，進(jìn)一步對殘差做平穩(wěn)性檢驗，結(jié)果表明殘差通過了適應(yīng)性檢驗，說明模型很好地提取了信息。（圖2、表9）

根據(jù)ARIMA（2，1，2）模型對苜蓿草進(jìn)口價格進(jìn)行動態(tài)預(yù)測，預(yù)測期間為12個月，預(yù)測結(jié)果如表10所示。（表10）

五、結(jié)論

預(yù)測結(jié)果可見，2017年進(jìn)口量、進(jìn)口價格將是上升的態(tài)勢，綜合人民幣匯率下跌、美國奶業(yè)不景氣、我國對國外進(jìn)口苜蓿草的依存度和國內(nèi)苜蓿草產(chǎn)品競爭力短時間難以改善的現(xiàn)狀更能證實。ARIMA模型只考慮了時間序列本身的特征，而沒有考慮其他一些因素的影響，即使這些因素在模型中是以隨機(jī)項來反映，但在預(yù)測值中是無法反映出來的。企業(yè)今后在苜蓿草進(jìn)口的預(yù)測和決策時注意結(jié)合產(chǎn)地氣候、產(chǎn)地市場變化、國內(nèi)相關(guān)政策、國內(nèi)牧草產(chǎn)量質(zhì)量價格等因素，提高對變動的敏感性和可預(yù)測性，從而得出更科學(xué)、更合理的進(jìn)口量預(yù)測結(jié)果。

主要參考文獻(xiàn)：

[1]王明利，楊春，胡向東，石自忠.關(guān)于苜蓿產(chǎn)業(yè)發(fā)展對糧食安全影響的政治經(jīng)濟(jì)學(xué)思考[J].草業(yè)科學(xué)，2012.29.12.

[2]石自忠，向碧云.我國牧草與草種進(jìn)出口貿(mào)易分析及預(yù)測——基于HP濾波法和ARMA模型[J].中國畜牧雜志，2013.49.6.

[3]王明利.2013中國牧草產(chǎn)業(yè)經(jīng)濟(jì)[M].北京：中國農(nóng)業(yè)出版社，2014.

[4]國家牧草產(chǎn)業(yè)技術(shù)體系.中國牧草主產(chǎn)區(qū)產(chǎn)業(yè)發(fā)展報告[M].北京：中國農(nóng)業(yè)大學(xué)出版社，2013.

[5]楊春，王明利，劉亞釗.中國的苜蓿草貿(mào)易——歷史變遷、未來趨勢與對策建議[J].草業(yè)科學(xué)，2011.28.9.

[6]江影舟，張潔冰，南志標(biāo)，王麗佳.中國苜蓿國際貿(mào)易競爭力分析[J].草業(yè)科學(xué)，2016.33.2.

[7]周忠麗，石自忠.世界苜蓿貿(mào)易的周期識別、預(yù)測及隨機(jī)沖擊影響[J].世界農(nóng)業(yè)，2013.11.

[8]劉玉鳳，王明利，胡向東，石自忠.美國苜蓿產(chǎn)業(yè)發(fā)展及其對中國的啟示[J].農(nóng)業(yè)展望，2014.10.8.

[9]汪武靜，王明利，呂官旺，劉玉鳳，石自忠.美國苜蓿貿(mào)易——趨勢、經(jīng)驗與啟示[J].草業(yè)科學(xué)，2016.33.3.

[10]韓雯.ARIMA模型在貴州省農(nóng)產(chǎn)品價格預(yù)測中的應(yīng)用——以辣椒為例[J].安徽農(nóng)業(yè)科學(xué)，2011.39.21.

[11]高鐵梅.計量經(jīng)濟(jì)分析方法與建模：EViews應(yīng)用及實例[M].北京：清華大學(xué)出版社，2006.

[12]石自忠，王明利，胡向東，崔姹.我國牧草種植成本收益變化與比較[J].草業(yè)科學(xué)，2017.34.4.

[13]豆明.2015中國牧草統(tǒng)計資料[M].北京：荷斯坦出版社，2015.

[14]胡向東，王明利.美國牧草生產(chǎn)及價格[J].草業(yè)科學(xué)，2013.30.5.

[15]劉亞釗，王明利，修長柏.我國牧草產(chǎn)品國際競爭力分析[J].農(nóng)業(yè)經(jīng)濟(jì)問題，2011.7.

[16]國家牧草產(chǎn)業(yè)技術(shù)體系.中國現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)可持續(xù)發(fā)展戰(zhàn)略研究牧草分冊[J].北京：中國農(nóng)業(yè)出版社，2016.

[17]侯冰凌，樊孝鳳.我國天然橡膠價格波動趨勢及短期預(yù)測——基于ARIMA模型的實證分析[J].價格理論與實踐，2016.