李 燕，賈春林，趙 偉，孫雅飛，姜文娟

（1. 山東農(nóng)業(yè)大學(xué)經(jīng)濟管理學(xué)院, 山東 泰安 271018；2. 山東省農(nóng)科院休閑農(nóng)業(yè)研究所, 山東 濟南 250100；3. 山東省地礦工程勘察院, 山東 濟南 250014；4. 山東省地質(zhì)博物館, 山東 濟南 250014）

在共同富裕和消費升級的背景下，發(fā)展牧草產(chǎn)業(yè)成為解決食物安全問題的重要途徑[1]。我國牧草需求旺盛，但國內(nèi)牧草產(chǎn)業(yè)的生產(chǎn)在技術(shù)與政策支持等方面存在限制，供需缺口近年來一直存在，需要以進口彌補國內(nèi)供給缺口[2-4]。牧草生產(chǎn)中苜蓿(Medicago sativa)相比其他種類牧草具有低成本高效益的特點[5]，在牧草進口中占據(jù)主導(dǎo)地位，同時苜蓿草進口價格存在約36個月的周期性變化特征[6]。

2020年新冠肺炎疫情發(fā)生后，我國農(nóng)業(yè)包括畜牧養(yǎng)殖業(yè)的生產(chǎn)與供應(yīng)受到一系列影響，一度出現(xiàn)飼料短缺、畜產(chǎn)品運輸受阻、依靠外調(diào)原料生產(chǎn)的養(yǎng)殖業(yè)生產(chǎn)中斷、進口產(chǎn)品人工及運輸成本明顯上升、國際市場農(nóng)產(chǎn)品價格波動預(yù)期不平穩(wěn)等一系列的問題[7-9]。苜蓿草生產(chǎn)遭受著本國疫情及進口因素的沖擊，產(chǎn)生新的價格變動，為此，研究疫情影響下苜蓿草價格波動，以探明疫情對牧草產(chǎn)業(yè)的實際影響，有利于幫助牧草產(chǎn)業(yè)在疫情中穩(wěn)定生產(chǎn)，尋找新的發(fā)展契機，并為日后牧草產(chǎn)業(yè)應(yīng)對突發(fā)事件提供經(jīng)驗建議。

1 材料與方法

1.1 數(shù)據(jù)來源

本研究采用2018年1月至2021年6月的月度數(shù)據(jù)，從疫情影響、苜蓿草進口情況和下游產(chǎn)品需求方面考慮設(shè)定變量，變量分別是苜蓿草價格(AP)、苜蓿草進口量(AI)、苜蓿草進口價格增長率(AG)、國內(nèi)生鮮乳價格(RP)、國內(nèi)月增長新冠疫情確診人數(shù)(ND)，數(shù)據(jù)序列完整。其中苜蓿草進口量與進口價格增長率數(shù)據(jù)來自中國海關(guān)并進行了整理，苜蓿草價格來自內(nèi)蒙古牧草產(chǎn)業(yè)發(fā)展協(xié)會、草都飼草料研究院，月增長疫情確診人數(shù)來自中華人民共和國國家衛(wèi)生健康委員會，國內(nèi)生鮮乳價格來自農(nóng)業(yè)農(nóng)村部監(jiān)測數(shù)據(jù)。

2020年1月至2021年6月我國苜蓿草價格頻繁波動，2020年上半年價格先漲后跌，2020年1月至3月價格大幅上漲，平均增長幅度在30%以上，苜蓿草價格達到最高為2 792 元·t-1，2020年后最低價為2 221 元·t-1，波動幅度達到26% (圖1)。 2020年7月至2021年6月苜蓿草價格持續(xù)波動，但總體存在上升趨勢。2020年以來生鮮乳價格與苜蓿草價格大致保持一致的變化趨勢，總的來說生鮮乳價格先降后升，2020年5月至2021年2月價格增長幅度為20%。

圖1 新冠肺炎疫情下苜蓿草價格與生鮮乳價格Figure 1 Line chart of the prices of alfalfa and raw milk for the COVID-19 pandemic

1.2 研究方法

使用VAR模型探究不同變量對苜蓿草價格波動的影響，VAR模型即向量自回歸模型(vector autoregressive, VAR)，VAR(p)模型的一般形式為：

式中：t= 1, 2, …,n；Yt為d維內(nèi)生變量的列向量，Yt-p為Yt的p階滯后變量，Xt為k維外生變量的列向量，n為樣本個數(shù)，p為模型滯后階數(shù)，B1…BP與H都是待估計的系數(shù)矩陣，εt是d維擾動列向量。本研究中Yt= (AP，AI，AG，ND，RP)T，對于不含外生變量的非限制向量自回歸模型可以表示為：

式中：t= 1, 2, …,n。公式(2)將系統(tǒng)中每一個內(nèi)生變量作為系統(tǒng)中所有內(nèi)生變量的滯后值的函數(shù)來構(gòu)造模型，VAR模型在經(jīng)濟學(xué)中適用于經(jīng)濟系統(tǒng)動態(tài)分析，能解釋各種沖擊對經(jīng)濟變量的影響。在試驗中通過格蘭杰因果關(guān)系檢驗分析影響因素的傳導(dǎo)過程，再利用脈沖響應(yīng)函數(shù)觀察不同期各影響因素對苜蓿草價格的影響形式，最后通過方差分解得到動態(tài)情況下不同因素對苜蓿草價格影響的貢獻率[10]。結(jié)果應(yīng)用Eviews10分析。

市場上的供給與需求變動始終是農(nóng)產(chǎn)品價格波動的主導(dǎo)因素，供給與需求的差異方向與差異程度影響價格波動方向與波動程度[11]。我國苜蓿草供給80%以上依賴進口，在供給上影響最為顯著，本研究從進口供給沖擊以及國內(nèi)需求拉動兩條路徑考察疫情因素下苜蓿草價格的波動原因及傳導(dǎo)過程，以苜蓿草進口量及苜蓿草進口價格增長率作為進口供給影響因素，苜蓿草生產(chǎn)最終支持畜產(chǎn)品生產(chǎn)，在產(chǎn)業(yè)鏈的作用下下游產(chǎn)品價格通過需求拉動上游生產(chǎn)資料價格，或者上游產(chǎn)品通過成本控制下游產(chǎn)品價格[12]，研究以生鮮乳價格體現(xiàn)畜牧業(yè)生產(chǎn)中對苜蓿草的需求。

2 實證分析過程

2.1 平穩(wěn)性檢驗

為防止出現(xiàn)偽回歸問題，使用ADF單位根檢驗法檢驗變量的平穩(wěn)性(表1)，結(jié)果表明，在5%的顯著性水平下各變量均符合同階平穩(wěn)要求，可以使用VAR模型進行研究。

表1 變量平穩(wěn)性Table 1 Stationarity test of variables

2.2 格蘭杰因果關(guān)系檢驗

用變量苜蓿草價格(AP)、苜蓿草進口量(AI)、苜蓿草進口價格增長率(AG)、國內(nèi)生鮮乳價格(RP)、國內(nèi)月增長新冠疫情確診數(shù)(ND)構(gòu)建五維向量自回歸模型，對變量進行格蘭杰因果檢驗來判斷變量間影響的傳導(dǎo)機制，使用AIC 與SC最小信息準(zhǔn)則判斷滯后階數(shù)，最終選擇滯后階數(shù)為5進行檢驗(表2)，由結(jié)果得出的變量間影響傳導(dǎo)過程如圖2所示。

圖2 變量間影響的傳導(dǎo)過程Figure 2 Transmission mechanism of influence between variables

表2 格蘭杰因果檢驗結(jié)果Table 2 Granger causality test results

2.3 脈沖響應(yīng)函數(shù)

VAR模型是穩(wěn)定的才可以保證結(jié)果的準(zhǔn)確性，使用AR方法進行模型的穩(wěn)定性檢驗，經(jīng)檢驗所有特征根都在單位圓范圍內(nèi)，說明模型穩(wěn)定，可以利用脈沖響應(yīng)函數(shù)進行分析。因為使用的是月度數(shù)據(jù)，響應(yīng)期便設(shè)置為35期。結(jié)果表明，通過建立生鮮乳價格、苜蓿草進口價格增長率和苜蓿草價格分別沖擊剩余變量的脈沖響應(yīng)函數(shù)，進口價格增長率、生鮮乳價格和疫情沖擊均會對苜蓿草價格產(chǎn)生短期顯著影響。

其中，進口價格增長率的正向沖擊會給苜蓿草價格帶來顯著的正面效應(yīng)，給苜蓿草進口量帶來顯著的負面效應(yīng)。響應(yīng)顯著的期限均為滯后一期，表明進口價格增長率對苜蓿草價格和苜蓿草進口量的短期影響大。

生鮮乳價格的正向沖擊會使苜蓿草價格產(chǎn)生顯著的正面響應(yīng)，使苜蓿草進口價格增長率產(chǎn)生顯著的負面響應(yīng)，響應(yīng)均為滯后一期，為短期影響，這主要與產(chǎn)品生產(chǎn)成本有關(guān)。

疫情沖擊給苜蓿草價格、進口價格增長率與進口量均帶來正面影響，但影響并不顯著，疫情沖擊使生鮮乳價格產(chǎn)生顯著的負面響應(yīng)，響應(yīng)的顯著持續(xù)時間為6個月，時間較長。

2.4 方差分解

方差分解能分析每一變量對苜蓿草價格(AP)變化的貢獻度，進一步評價不同變量的重要性，對苜蓿草價格變動在10期內(nèi)受其他4個變量的貢獻影響的分解結(jié)果(表3)表明，苜蓿草價格(AP)在最初期自身貢獻率最大，達到100%，但其自身貢獻率在第10期降低到19.49%，不再處于最大貢獻地位。進口量(AI)對苜蓿草價格(AP)的貢獻在逐步增加，在第10期實現(xiàn)最大貢獻，最大貢獻為29.43%。進口價格增長率(AG)對苜蓿草價格(AP)的貢獻率在第2期最高為18.98%，隨后經(jīng)歷波動過程在第10期的貢獻率為12.86%，說明短期內(nèi)進口價格增長率(AG)對苜蓿草價格(AP)的影響程度大于長期，總體上為苜蓿草價格(AP)帶來的影響較小。疫情(ND)對苜蓿草價格(AP)的貢獻率先增加再減少，在第4期達到25.75%，與苜蓿草價格(AP)的自身貢獻率十分接近，基本實現(xiàn)最大貢獻，這說明疫情初期會為苜蓿草價格(AP)帶來較大影響，在疫情長期存在的環(huán)境下對苜蓿草價格(AP)的影響會逐漸削弱。生鮮乳價格(RP)對苜蓿草價格(AP)的貢獻率在第2期達到最大(26.73%)，滯后期增加使得貢獻率降低，但第10期仍保持在20%以上。

表3 方差分解結(jié)果Table 3 Variance decomposition results

3 結(jié)果與討論

苜蓿草價格在疫情環(huán)境下所受影響的傳導(dǎo)機制為疫情影響生鮮乳價格，生鮮乳價格影響進口苜蓿草價格增長率，進口苜蓿草價格增長率影響苜蓿草價格。

疫情影響下生鮮乳價格下降的影響時間為6個月。疫情初期發(fā)生在人口密集的城鄉(xiāng)區(qū)域，社會生產(chǎn)生活相對集中，生鮮乳作為日用消費品在居家隔離以及社會停工停產(chǎn)的環(huán)境下消費需求大幅降低，到2020年2月初云南奶企的銷量下降超過50%。在生產(chǎn)端，2020年春節(jié)過后氣溫回升，牧場在一般情況下產(chǎn)奶量會相對上漲，疫情發(fā)生后政府、奶企與牧場等多方協(xié)作，致力于實現(xiàn)奶業(yè)生產(chǎn)穩(wěn)定，生鮮乳生產(chǎn)沒有受飼料運輸不暢等原因影響而產(chǎn)量下降，相反還略有上升，2020年上半年生鮮乳產(chǎn)量同比增長7.9%，進口量同比增長5.8%，在總體供求失衡的市場環(huán)境下生鮮乳價格略有降低。相比疫情嚴重地區(qū)來說，我國牧區(qū)地廣人稀且對外交流不頻繁，疫情傳播速度較慢，生產(chǎn)生活受限程度小，疫情影響作用減弱且滯后[13]，因此疫情先沖擊更靠近生產(chǎn)鏈下游的生鮮乳價格，對生鮮乳價格影響更直接。

受到生鮮乳價格下降的沖擊，苜蓿草進口價格增長率短期上升。在實際市場中，生鮮乳價格與進口苜蓿草價格之間存在相互影響的現(xiàn)象，生鮮乳價格正向影響進口苜蓿草價格，這種現(xiàn)象穩(wěn)定并且持續(xù)存在[14]，本研究中在疫情的特殊情況下生鮮乳價格與進口價格增長率呈短期反向變動，在生鮮乳價格下降的情況下進口價格增長率呈現(xiàn)上升趨勢。這種情況可能是受疫情造成的國內(nèi)市場恐慌加上國外疫情尚未擴張的影響，一方面生鮮乳價格受需求影響下跌，產(chǎn)品質(zhì)量受更多關(guān)注，進口苜蓿草是提升生鮮乳質(zhì)量的有力保證；另一方面牧場積極備貨防止斷糧，苜蓿草市場需求旺，我國畜牧業(yè)種養(yǎng)分離嚴重，牧草產(chǎn)銷區(qū)域不一致，疫情下國內(nèi)牧草運輸更加困難，在國際海運暫未受影響時進口需求旺盛[15]。從影響的時間為一個月也可以看出生鮮乳價格下降引起苜蓿草進口價格增長率上升的影響并不是長期均衡的。

苜蓿草進口價格增長率提升會在短期影響下使苜蓿草價格增長。從替代品競爭方面來說，進口苜蓿草是國內(nèi)苜蓿草的最強競爭者，因為高蛋白含量的苜蓿草是提升奶牛產(chǎn)奶量的保證，國產(chǎn)苜蓿草質(zhì)量相對進口苜蓿草來說仍有差距，我國苜蓿草供給中從美國、西班牙為代表的國家進口部分占80%以上，進口價格是影響國內(nèi)苜蓿草需求及其價格的重要因素，并且在一定的替代關(guān)系下雙方價格同向變動。美國和西班牙是我國苜蓿草進口的主要來源國，在2020年疫情中苜蓿種植面積略有擴張，美國沒有中斷出口交易，西班牙則略受影響。疫情對苜蓿草進口的影響主要體現(xiàn)在海運上，國際疫情形勢下海運成本上升，到2020年4月船運公司普遍將海運費率上調(diào)50～100 USD，苜蓿草進口成本上漲。疫情影響下我國生產(chǎn)活動大幅減少，出口近乎停滯，往來交易困難，進口訂單順延后出現(xiàn)進口量先暫時穩(wěn)定后扎堆到港的現(xiàn)象，因此苜蓿草進口量的短期影響并不顯著。

對苜蓿草價格來說除自身貢獻外生鮮乳價格一直對其發(fā)揮著重要的影響作用，我國實施振興奶業(yè)苜蓿發(fā)展行動，以提升奶業(yè)發(fā)展質(zhì)量為目的積極調(diào)整苜蓿草的生產(chǎn)，我國苜蓿草生產(chǎn)存在政策效應(yīng)，同時作為經(jīng)濟類作物要以供給側(cè)結(jié)構(gòu)性改革為背景積極實現(xiàn)以需定產(chǎn)，這已成為苜蓿草獲得持續(xù)穩(wěn)定收益的重要保證，奶業(yè)的發(fā)展需要決定著下一步苜蓿草的生產(chǎn)，折射到價格上就是生鮮乳價格對苜蓿草價格的影響作用。進口量對苜蓿草價格的影響程度隨時間延長而逐漸加深，這是受進口運輸過程以及苜蓿草生產(chǎn)需要較長時間所引起的滯后效應(yīng)所引起的，表明應(yīng)重點關(guān)注進口量的長期變化來調(diào)節(jié)生產(chǎn)。疫情帶來的影響在4個月內(nèi)達到頂峰，隨后影響效果減弱，說明我國企業(yè)應(yīng)對疫情的能力在不斷增強，抓緊時間調(diào)控生產(chǎn)、減輕損失的政策已發(fā)揮了重要作用。

總體來看，研究得出的傳導(dǎo)過程是合理的，并且傳導(dǎo)過程在短期內(nèi)均有顯著效果，長期不明顯；苜蓿草價格自身貢獻率在短期內(nèi)最大，在長期苜蓿草進口量及生鮮乳價格的貢獻率更大，疫情給苜蓿草價格帶來的影響會在4個月后減弱。

4 結(jié)論與建議

4.1 結(jié)論

除自身貢獻外，短期苜蓿草價格受進口價格增長率和生鮮乳價格影響較大，苜蓿草價格更多受需求影響。疫情主要通過影響生鮮乳價格來提升苜蓿草價格，影響顯著的時間為6個月，在疫情環(huán)境下生鮮乳價格的平穩(wěn)性是影響苜蓿草價格波動的重要因素，應(yīng)注意生鮮乳價格在公共安全事件發(fā)生時對上游產(chǎn)業(yè)的調(diào)節(jié)作用。從長期來看苜蓿草價格更多地受到進口量的影響。

4.2 建議

1)加快畜牧產(chǎn)業(yè)發(fā)展，實現(xiàn)草飼一體化。

推進種養(yǎng)結(jié)合，實現(xiàn)草飼一體化有利于降低運輸風(fēng)險，同時減輕下游產(chǎn)品價格對苜蓿草價格的沖擊，幫助企業(yè)實現(xiàn)規(guī)模效益與生態(tài)效益[16-17]。這需要國家在政策層面給予讓步與政策支持，建立完善的金融服務(wù)保障機制，在技術(shù)層面致力于破解生產(chǎn)難題，并對企業(yè)加強宣傳指引。在牧草流通服務(wù)支持上仍需建立完備的應(yīng)急運輸服務(wù)機制，避免意外災(zāi)害對牧草運輸產(chǎn)生持續(xù)影響[18]。

2)完善市場風(fēng)險預(yù)警機制，增強風(fēng)險調(diào)控能力。

苜蓿草價格受短期沖擊影響顯著，為保障牧草產(chǎn)業(yè)平穩(wěn)健康發(fā)展，我國仍需建立更加完善的風(fēng)險研究、預(yù)防、監(jiān)測與應(yīng)急處理機制。擴大風(fēng)險監(jiān)控的時空領(lǐng)域，對國內(nèi)外風(fēng)險分級評估，以降低國際市場價格波動的沖擊，加強風(fēng)險預(yù)測，積極采取有效預(yù)防措施，對已發(fā)生的問題進行靶向處理，明確各主體責(zé)任，及時調(diào)控，提升管理效率[19]。

3)建立數(shù)據(jù)云服務(wù)平臺，實現(xiàn)資源全方位共享。

疫情防控經(jīng)驗表明，我國疫情以省市為單位管控的情況下解決疫情沖擊需要精準(zhǔn)導(dǎo)向，推動建立牧草產(chǎn)業(yè)大數(shù)據(jù)服務(wù)平臺，可以了解存在疫情風(fēng)險地區(qū)的訴求，并及時調(diào)度資源，即注意為服務(wù)對象提供精準(zhǔn)服務(wù)，完善資源服務(wù)內(nèi)容，實現(xiàn)全方位數(shù)字化產(chǎn)業(yè)管理。國家要支持數(shù)據(jù)云服務(wù)技術(shù)發(fā)展，聯(lián)結(jié)市場、政府與企業(yè)進行合作，致力于豐富大數(shù)據(jù)內(nèi)容，及時更新，為信息使用者提供權(quán)威性與透明度共存的產(chǎn)業(yè)數(shù)據(jù)[20]，在此基礎(chǔ)上注意平臺向中小企業(yè)開放，調(diào)節(jié)資源分配不均的現(xiàn)狀。