“雙碳”目標背景下我國畜牧業(yè)碳排放效率實證分析

王維新 李焱求

（1.山東畜牧獸醫(yī)職業(yè)學院，山東 濰坊 261061；2.培材大學，韓國 大田廣域 340934）

0 引言

隨著全球溫室效應的加劇，二氧化碳排放問題逐漸得到重視。全球氣候變化危及人類安全和健康，影響所有人賴以生存的自然生態(tài)系統(tǒng)和資源。近年來，碳排放問題愈發(fā)嚴峻，除了新型冠狀病毒肺炎疫情導致2020年碳排放量下降外，自2015年巴黎締約方大會召開以來，全球碳排放量每年都在上升。我國作為碳排放的主要國家，2020年承諾實現(xiàn)2030年“碳達峰”與2060年“碳中和”目標，要求全社會科學有序降碳，促進經(jīng)濟社會發(fā)展全面綠色轉型。黨的十八大召開后，我國的發(fā)展理念發(fā)生了根本性變化，環(huán)境保護不再被視為經(jīng)濟發(fā)展的負擔，而是被視為高質(zhì)量發(fā)展的驅動力。畜牧業(yè)是公認的全球氣候變暖的貢獻者之一，并非傳統(tǒng)認知中的“異味產(chǎn)業(yè)”，其碳排放量約占全球總排放量的1/7，倡導畜牧業(yè)節(jié)能減排已達成共識。我國是世界上最大的畜牧生產(chǎn)國，隨著城鄉(xiāng)居民膳食結構和營養(yǎng)水平的不斷升級、農(nóng)業(yè)農(nóng)村經(jīng)濟發(fā)展和農(nóng)民增收的迫切需求，我國畜牧業(yè)的碳排放量正不斷向高位趨近，亟待依托產(chǎn)業(yè)結構調(diào)整和效率變革，加快畜牧經(jīng)濟增長向低碳化轉型。因此，開展畜牧業(yè)碳排放效率研究對“雙碳”目標背景下促進經(jīng)濟社會發(fā)展具有一定緊迫性和必要性。

1 研究現(xiàn)狀

目前，國內(nèi)學界對農(nóng)業(yè)碳排放、碳排放效率及其影響因素的研究較多，但針對全國畜牧業(yè)碳排放效率實證分析的相關研究報道非常少，且存在一定的不足。例如，實證研究基于我國省域、某一特定區(qū)域或地市級視角，全面性不足；研究方法以傳統(tǒng)的DEA模型為主，未將非期望產(chǎn)出和環(huán)境變量納入研究范疇，嚴謹性不足；采用的數(shù)據(jù)相對較早，科學性不足等。鑒于此，筆者運用非期望產(chǎn)出的三階段超效率SBMDEA模型和Malmquist指數(shù)，分別對2002—2020年國家層面、四大經(jīng)濟帶（東部、西部、中部和東北）和31個?。ㄗ灾螀^(qū)、直轄市）的畜牧業(yè)碳排放效率進行測度分析，以期為“雙碳”目標下畜牧業(yè)減排和畜牧業(yè)高質(zhì)量發(fā)展提供實證依據(jù)。

2 研究方法

2.1 三階段超效率SBM-DEA模型

為了解決傳統(tǒng)數(shù)據(jù)包絡分析（Data Envelopment Analysis，DEA）的DUM效率值大于1而無法有效評價的情況，筆者采用包含非意愿產(chǎn)出基于松弛變量的三階段超效率DEA模型（Slacks-Based Measure，SBM）。三階段超效率SBM模型如公式（1）所示：

2.2 隨機前沿（SFA）模型

為剔除隨機因素和環(huán)境因素對效率值的影響，構建SFA模型，如公式（2）所示：

3 指標選取與數(shù)據(jù)來源

該研究的實證數(shù)據(jù)來源為2002—2020年《中國畜牧業(yè)統(tǒng)計年鑒》、國家統(tǒng)計局網(wǎng)站、《山東省畜禽養(yǎng)殖場用地參考指標》和市場調(diào)研。筆者選取畜牧業(yè)從業(yè)人員數(shù)量、畜牧業(yè)固定資產(chǎn)投資額、畜禽飼養(yǎng)土地面積作為投入指標；選取人均畜牧業(yè)產(chǎn)值作為意愿產(chǎn)出指標，選取畜牧業(yè)碳排放量作為非意愿產(chǎn)出指標；選取產(chǎn)業(yè)結構、經(jīng)濟發(fā)展水平和區(qū)域創(chuàng)新能力作為環(huán)境指標（見表1）。

表1 碳排放效率測算指標

依據(jù)聯(lián)合國政府間氣候變化專門委員會（Intergovernmental Panel on Climate Change，IPCC）提出的碳排放因子法，測算《中國畜牧業(yè)統(tǒng)計年鑒》中包含的牛、馬、家禽、兔、騾、驢、豬和駱駝等畜禽的CO排放當量。測算公式（3）如下：

畜禽具體溫室氣體排放系數(shù)見表2。

表2 畜牧業(yè)溫室氣體排放系數(shù) kg/（頭·a）

4 實證分析

運用MaxDEA軟件對數(shù)據(jù)進行第一階段的效率測算后，在第二階段納入產(chǎn)業(yè)結構、經(jīng)濟發(fā)展水平和科技水平3個環(huán)境變量，用Frontier4.1軟件對2002—2020年國家、31個?。ㄗ灾螀^(qū)、直轄市）和4個經(jīng)濟帶的投入松弛變量進行SFA回歸，利用Excel軟件剔除環(huán)境因素和管理無效率對松弛變量的影響，并重新調(diào)整投入變量的值，將各個DUM置于相同的水平下。第三階段同樣利用超效率SBM模型對第二階段調(diào)整后的投入值運行效率重新進行測度，結果如表3所示。

表3 第三階段碳排放效率

由表3可知，非期望產(chǎn)出的三階段超效率SBMDEA模型測算結果顯示的畜牧業(yè)碳排放效率與該地區(qū)畜牧業(yè)飼養(yǎng)規(guī)模并不完全對應，運算結果代表的是一種投入產(chǎn)出比，并不能體現(xiàn)各個區(qū)域在產(chǎn)業(yè)發(fā)展方面是否具有競爭性。另外，對于測算碳排放效率占優(yōu)的區(qū)域，并不代表其沒有改進空間。

為更加清晰地展示2002—2020年國家層面和四大經(jīng)濟帶畜牧業(yè)碳排放效率，繪制年平均值為單元的折線圖（見圖1）。

圖1 2002—2020年我國年均碳排放效率

從圖1整體來看，全國碳排放效率呈現(xiàn)出先降后升趨勢，2006—2011年增長趨緩，可能與國家頒布的《節(jié)能目標責任和評價考核實施方案》取得一定效果有關；2012—2018年呈平滑下降狀態(tài)，可能與國家頒布《畜禽規(guī)模養(yǎng)殖污染防治條例》和助力畜牧業(yè)發(fā)展的相關政策導致小規(guī)模散養(yǎng)戶淘汰，大規(guī)模養(yǎng)殖場蓬勃發(fā)展相關；2018年到波谷的0.044后開始回暖，可能與2018年起中美貿(mào)易戰(zhàn)出口受挫，禽畜疾病頻發(fā)，畜牧業(yè)一片蕭條，產(chǎn)業(yè)結構升級和新舊動能轉換開始發(fā)力有關。

從各經(jīng)濟帶來看，四大經(jīng)濟帶的變化趨勢與國家層面變化基本一致，東部、中部和西部的碳排放效率都高于國家平均值，東北地區(qū)的碳排放效率遠低于東部、中部、西部地區(qū)和國家平均水平。東部地區(qū)研究期內(nèi)初期碳排放效率領先于其他地區(qū)，但2009年后被中部地區(qū)追平，甚至在有的年份被超越，在2018年以后其效率值又被西部地區(qū)趕超。這可能與東部地區(qū)部分產(chǎn)業(yè)轉移、中西部地區(qū)受益沿海產(chǎn)業(yè)轉移的技術擴散，產(chǎn)業(yè)結構調(diào)整相關；而東北地區(qū)相對土地資源豐富、畜牧業(yè)固定資產(chǎn)投資不足和經(jīng)濟發(fā)展落后，在一定程度上制約了當?shù)匦竽翗I(yè)碳排放效率的提高。

從省域來看，我國31個?。ㄗ灾螀^(qū)、直轄市）在研究期間的碳排放效率表現(xiàn)出顯著的區(qū)域差異性，最高浙江省的效率均值為3.221，最低黑龍江省的效率均值為0.027。從面板數(shù)據(jù)來看，浙江省和山東省多年的效率值都大于1，是其他?。ㄗ灾螀^(qū)、直轄市）的追趕目標。效率值低于0.1的遼寧省、寧夏回族自治區(qū)、甘肅省、內(nèi)蒙古自治區(qū)、新疆維吾爾自治區(qū)、海南省和黑龍江省，受區(qū)域創(chuàng)新能力、產(chǎn)業(yè)結構和經(jīng)濟發(fā)展水平影響，畜牧業(yè)發(fā)展方式以傳統(tǒng)的粗放型為主，發(fā)展模式分散且小，產(chǎn)業(yè)化、規(guī)模化、標準化水平低，污染源點多且面廣，且受地理和自然環(huán)境等要素稟賦限制，該地區(qū)畜牧業(yè)碳排放效率一直處于較低水平。但該區(qū)域畜牧業(yè)碳排放效率提升空間大，可以通過采用飼料新配方、引進新技術和發(fā)展養(yǎng)殖新模式等提高畜牧業(yè)碳排放效率。

5 結論與建議

筆者基于三階段超效率SBM-DEA模型測算得到國家層面、四大經(jīng)濟帶和31個省（自治區(qū)、直轄市）的畜牧業(yè)二氧化碳排放效率值。根據(jù)結果分析，得出以下結論：①我國畜牧業(yè)的二氧化碳排放效率值出現(xiàn)高—低—高的趨勢，區(qū)域間碳排放效率差異較顯著，剔除新型冠狀病毒肺炎疫情因素，近年來東部地區(qū)、中部地區(qū)和西部地區(qū)的效率排放基本持平，說明畜牧業(yè)的碳排放問題開始出現(xiàn)好轉；②2002—2020年，我國畜牧業(yè)全國層面碳排放效率偏低，均值為0.134，節(jié)能減排潛力較大，而各?。ㄗ灾螀^(qū)、直轄市）的碳排放效率存在顯著差異，其中黑龍江省和海南省效率最低，平均排放效率僅為0.030左右，需要重點關注； ③我國幅員遼闊，各地資源稟賦不一，使得畜牧業(yè)產(chǎn)業(yè)結構、從業(yè)人員規(guī)模、地區(qū)經(jīng)濟發(fā)展水平和城鎮(zhèn)化率等成為碳排放率不高的主要影響因素。從上述分析結果可以看出，我國畜牧業(yè)碳排放實現(xiàn)2030年達峰這一目標仍任重而道遠。

為有效控制畜牧業(yè)碳排放量，進一步提高排放效率，爭取早日實現(xiàn)“雙碳”目標，根據(jù)上述結論，筆者提出以下建議。一是各級政府應以《關于完整準確全面貫徹新發(fā)展理念做好碳達峰碳中和工作的意見》和《2030年前碳達峰行動方案》為綱領，以《“十四五”全國畜牧獸醫(yī)行業(yè)發(fā)展規(guī)劃》為抓手，結合實際制訂和落實區(qū)域發(fā)展規(guī)劃。二是各區(qū)域畜牧業(yè)減排工作的開展要因地而異、因類而異和因地制宜，并以綠色和循環(huán)發(fā)展理念指導發(fā)展畜牧業(yè)，實現(xiàn)環(huán)境和經(jīng)濟協(xié)調(diào)發(fā)展。三是一方面加大畜牧業(yè)財政支出力度，建立支持畜牧業(yè)低碳綠色循環(huán)發(fā)展的長效機制，實施積極的技術補貼措施，對引進、采用和研發(fā)新技術的規(guī)模化企業(yè)給予財政補貼和稅收減免；另一方面加大畜牧業(yè)固定資產(chǎn)投資力度和畜牧業(yè)從業(yè)人員培訓投資力度，為畜牧業(yè)高質(zhì)量發(fā)展 提供支持。