王維新 李焱求
(1.山東畜牧獸醫(yī)職業(yè)學(xué)院,山東 濰坊 261061;2.培材大學(xué),韓國(guó) 大田廣域 340934)
隨著全球溫室效應(yīng)的加劇,二氧化碳排放問(wèn)題逐漸得到重視。全球氣候變化危及人類安全和健康,影響所有人賴以生存的自然生態(tài)系統(tǒng)和資源。近年來(lái),碳排放問(wèn)題愈發(fā)嚴(yán)峻,除了新型冠狀病毒肺炎疫情導(dǎo)致2020年碳排放量下降外,自2015年巴黎締約方大會(huì)召開(kāi)以來(lái),全球碳排放量每年都在上升。我國(guó)作為碳排放的主要國(guó)家,2020年承諾實(shí)現(xiàn)2030年“碳達(dá)峰”與2060年“碳中和”目標(biāo),要求全社會(huì)科學(xué)有序降碳,促進(jìn)經(jīng)濟(jì)社會(huì)發(fā)展全面綠色轉(zhuǎn)型。黨的十八大召開(kāi)后,我國(guó)的發(fā)展理念發(fā)生了根本性變化,環(huán)境保護(hù)不再被視為經(jīng)濟(jì)發(fā)展的負(fù)擔(dān),而是被視為高質(zhì)量發(fā)展的驅(qū)動(dòng)力。畜牧業(yè)是公認(rèn)的全球氣候變暖的貢獻(xiàn)者之一,并非傳統(tǒng)認(rèn)知中的“異味產(chǎn)業(yè)”,其碳排放量約占全球總排放量的1/7,倡導(dǎo)畜牧業(yè)節(jié)能減排已達(dá)成共識(shí)。我國(guó)是世界上最大的畜牧生產(chǎn)國(guó),隨著城鄉(xiāng)居民膳食結(jié)構(gòu)和營(yíng)養(yǎng)水平的不斷升級(jí)、農(nóng)業(yè)農(nóng)村經(jīng)濟(jì)發(fā)展和農(nóng)民增收的迫切需求,我國(guó)畜牧業(yè)的碳排放量正不斷向高位趨近,亟待依托產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和效率變革,加快畜牧經(jīng)濟(jì)增長(zhǎng)向低碳化轉(zhuǎn)型。因此,開(kāi)展畜牧業(yè)碳排放效率研究對(duì)“雙碳”目標(biāo)背景下促進(jìn)經(jīng)濟(jì)社會(huì)發(fā)展具有一定緊迫性和必要性。
目前,國(guó)內(nèi)學(xué)界對(duì)農(nóng)業(yè)碳排放、碳排放效率及其影響因素的研究較多,但針對(duì)全國(guó)畜牧業(yè)碳排放效率實(shí)證分析的相關(guān)研究報(bào)道非常少,且存在一定的不足。例如,實(shí)證研究基于我國(guó)省域、某一特定區(qū)域或地市級(jí)視角,全面性不足;研究方法以傳統(tǒng)的DEA模型為主,未將非期望產(chǎn)出和環(huán)境變量納入研究范疇,嚴(yán)謹(jǐn)性不足;采用的數(shù)據(jù)相對(duì)較早,科學(xué)性不足等。鑒于此,筆者運(yùn)用非期望產(chǎn)出的三階段超效率SBMDEA模型和Malmquist指數(shù),分別對(duì)2002—2020年國(guó)家層面、四大經(jīng)濟(jì)帶(東部、西部、中部和東北)和31個(gè)?。ㄗ灾螀^(qū)、直轄市)的畜牧業(yè)碳排放效率進(jìn)行測(cè)度分析,以期為“雙碳”目標(biāo)下畜牧業(yè)減排和畜牧業(yè)高質(zhì)量發(fā)展提供實(shí)證依據(jù)。
為了解決傳統(tǒng)數(shù)據(jù)包絡(luò)分析(Data Envelopment Analysis,DEA)的DUM效率值大于1而無(wú)法有效評(píng)價(jià)的情況,筆者采用包含非意愿產(chǎn)出基于松弛變量的三階段超效率DEA模型(Slacks-Based Measure,SBM)。三階段超效率SBM模型如公式(1)所示:
為剔除隨機(jī)因素和環(huán)境因素對(duì)效率值的影響,構(gòu)建SFA模型,如公式(2)所示:
該研究的實(shí)證數(shù)據(jù)來(lái)源為2002—2020年《中國(guó)畜牧業(yè)統(tǒng)計(jì)年鑒》、國(guó)家統(tǒng)計(jì)局網(wǎng)站、《山東省畜禽養(yǎng)殖場(chǎng)用地參考指標(biāo)》和市場(chǎng)調(diào)研。筆者選取畜牧業(yè)從業(yè)人員數(shù)量、畜牧業(yè)固定資產(chǎn)投資額、畜禽飼養(yǎng)土地面積作為投入指標(biāo);選取人均畜牧業(yè)產(chǎn)值作為意愿產(chǎn)出指標(biāo),選取畜牧業(yè)碳排放量作為非意愿產(chǎn)出指標(biāo);選取產(chǎn)業(yè)結(jié)構(gòu)、經(jīng)濟(jì)發(fā)展水平和區(qū)域創(chuàng)新能力作為環(huán)境指標(biāo)(見(jiàn)表1)。
表1 碳排放效率測(cè)算指標(biāo)
依據(jù)聯(lián)合國(guó)政府間氣候變化專門(mén)委員會(huì)(Intergovernmental Panel on Climate Change,IPCC)提出的碳排放因子法,測(cè)算《中國(guó)畜牧業(yè)統(tǒng)計(jì)年鑒》中包含的牛、馬、家禽、兔、騾、驢、豬和駱駝等畜禽的CO排放當(dāng)量。測(cè)算公式(3)如下:
畜禽具體溫室氣體排放系數(shù)見(jiàn)表2。
表2 畜牧業(yè)溫室氣體排放系數(shù) kg/(頭·a)
運(yùn)用MaxDEA軟件對(duì)數(shù)據(jù)進(jìn)行第一階段的效率測(cè)算后,在第二階段納入產(chǎn)業(yè)結(jié)構(gòu)、經(jīng)濟(jì)發(fā)展水平和科技水平3個(gè)環(huán)境變量,用Frontier4.1軟件對(duì)2002—2020年國(guó)家、31個(gè)省(自治區(qū)、直轄市)和4個(gè)經(jīng)濟(jì)帶的投入松弛變量進(jìn)行SFA回歸,利用Excel軟件剔除環(huán)境因素和管理無(wú)效率對(duì)松弛變量的影響,并重新調(diào)整投入變量的值,將各個(gè)DUM置于相同的水平下。第三階段同樣利用超效率SBM模型對(duì)第二階段調(diào)整后的投入值運(yùn)行效率重新進(jìn)行測(cè)度,結(jié)果如表3所示。
表3 第三階段碳排放效率
由表3可知,非期望產(chǎn)出的三階段超效率SBMDEA模型測(cè)算結(jié)果顯示的畜牧業(yè)碳排放效率與該地區(qū)畜牧業(yè)飼養(yǎng)規(guī)模并不完全對(duì)應(yīng),運(yùn)算結(jié)果代表的是一種投入產(chǎn)出比,并不能體現(xiàn)各個(gè)區(qū)域在產(chǎn)業(yè)發(fā)展方面是否具有競(jìng)爭(zhēng)性。另外,對(duì)于測(cè)算碳排放效率占優(yōu)的區(qū)域,并不代表其沒(méi)有改進(jìn)空間。
為更加清晰地展示2002—2020年國(guó)家層面和四大經(jīng)濟(jì)帶畜牧業(yè)碳排放效率,繪制年平均值為單元的折線圖(見(jiàn)圖1)。
圖1 2002—2020年我國(guó)年均碳排放效率
從圖1整體來(lái)看,全國(guó)碳排放效率呈現(xiàn)出先降后升趨勢(shì),2006—2011年增長(zhǎng)趨緩,可能與國(guó)家頒布的《節(jié)能目標(biāo)責(zé)任和評(píng)價(jià)考核實(shí)施方案》取得一定效果有關(guān);2012—2018年呈平滑下降狀態(tài),可能與國(guó)家頒布《畜禽規(guī)模養(yǎng)殖污染防治條例》和助力畜牧業(yè)發(fā)展的相關(guān)政策導(dǎo)致小規(guī)模散養(yǎng)戶淘汰,大規(guī)模養(yǎng)殖場(chǎng)蓬勃發(fā)展相關(guān);2018年到波谷的0.044后開(kāi)始回暖,可能與2018年起中美貿(mào)易戰(zhàn)出口受挫,禽畜疾病頻發(fā),畜牧業(yè)一片蕭條,產(chǎn)業(yè)結(jié)構(gòu)升級(jí)和新舊動(dòng)能轉(zhuǎn)換開(kāi)始發(fā)力有關(guān)。
從各經(jīng)濟(jì)帶來(lái)看,四大經(jīng)濟(jì)帶的變化趨勢(shì)與國(guó)家層面變化基本一致,東部、中部和西部的碳排放效率都高于國(guó)家平均值,東北地區(qū)的碳排放效率遠(yuǎn)低于東部、中部、西部地區(qū)和國(guó)家平均水平。東部地區(qū)研究期內(nèi)初期碳排放效率領(lǐng)先于其他地區(qū),但2009年后被中部地區(qū)追平,甚至在有的年份被超越,在2018年以后其效率值又被西部地區(qū)趕超。這可能與東部地區(qū)部分產(chǎn)業(yè)轉(zhuǎn)移、中西部地區(qū)受益沿海產(chǎn)業(yè)轉(zhuǎn)移的技術(shù)擴(kuò)散,產(chǎn)業(yè)結(jié)構(gòu)調(diào)整相關(guān);而東北地區(qū)相對(duì)土地資源豐富、畜牧業(yè)固定資產(chǎn)投資不足和經(jīng)濟(jì)發(fā)展落后,在一定程度上制約了當(dāng)?shù)匦竽翗I(yè)碳排放效率的提高。
從省域來(lái)看,我國(guó)31個(gè)?。ㄗ灾螀^(qū)、直轄市)在研究期間的碳排放效率表現(xiàn)出顯著的區(qū)域差異性,最高浙江省的效率均值為3.221,最低黑龍江省的效率均值為0.027。從面板數(shù)據(jù)來(lái)看,浙江省和山東省多年的效率值都大于1,是其他?。ㄗ灾螀^(qū)、直轄市)的追趕目標(biāo)。效率值低于0.1的遼寧省、寧夏回族自治區(qū)、甘肅省、內(nèi)蒙古自治區(qū)、新疆維吾爾自治區(qū)、海南省和黑龍江省,受區(qū)域創(chuàng)新能力、產(chǎn)業(yè)結(jié)構(gòu)和經(jīng)濟(jì)發(fā)展水平影響,畜牧業(yè)發(fā)展方式以傳統(tǒng)的粗放型為主,發(fā)展模式分散且小,產(chǎn)業(yè)化、規(guī)?;?、標(biāo)準(zhǔn)化水平低,污染源點(diǎn)多且面廣,且受地理和自然環(huán)境等要素稟賦限制,該地區(qū)畜牧業(yè)碳排放效率一直處于較低水平。但該區(qū)域畜牧業(yè)碳排放效率提升空間大,可以通過(guò)采用飼料新配方、引進(jìn)新技術(shù)和發(fā)展養(yǎng)殖新模式等提高畜牧業(yè)碳排放效率。
筆者基于三階段超效率SBM-DEA模型測(cè)算得到國(guó)家層面、四大經(jīng)濟(jì)帶和31個(gè)?。ㄗ灾螀^(qū)、直轄市)的畜牧業(yè)二氧化碳排放效率值。根據(jù)結(jié)果分析,得出以下結(jié)論:①我國(guó)畜牧業(yè)的二氧化碳排放效率值出現(xiàn)高—低—高的趨勢(shì),區(qū)域間碳排放效率差異較顯著,剔除新型冠狀病毒肺炎疫情因素,近年來(lái)東部地區(qū)、中部地區(qū)和西部地區(qū)的效率排放基本持平,說(shuō)明畜牧業(yè)的碳排放問(wèn)題開(kāi)始出現(xiàn)好轉(zhuǎn);②2002—2020年,我國(guó)畜牧業(yè)全國(guó)層面碳排放效率偏低,均值為0.134,節(jié)能減排潛力較大,而各?。ㄗ灾螀^(qū)、直轄市)的碳排放效率存在顯著差異,其中黑龍江省和海南省效率最低,平均排放效率僅為0.030左右,需要重點(diǎn)關(guān)注; ③我國(guó)幅員遼闊,各地資源稟賦不一,使得畜牧業(yè)產(chǎn)業(yè)結(jié)構(gòu)、從業(yè)人員規(guī)模、地區(qū)經(jīng)濟(jì)發(fā)展水平和城鎮(zhèn)化率等成為碳排放率不高的主要影響因素。從上述分析結(jié)果可以看出,我國(guó)畜牧業(yè)碳排放實(shí)現(xiàn)2030年達(dá)峰這一目標(biāo)仍任重而道遠(yuǎn)。
為有效控制畜牧業(yè)碳排放量,進(jìn)一步提高排放效率,爭(zhēng)取早日實(shí)現(xiàn)“雙碳”目標(biāo),根據(jù)上述結(jié)論,筆者提出以下建議。一是各級(jí)政府應(yīng)以《關(guān)于完整準(zhǔn)確全面貫徹新發(fā)展理念做好碳達(dá)峰碳中和工作的意見(jiàn)》和《2030年前碳達(dá)峰行動(dòng)方案》為綱領(lǐng),以《“十四五”全國(guó)畜牧獸醫(yī)行業(yè)發(fā)展規(guī)劃》為抓手,結(jié)合實(shí)際制訂和落實(shí)區(qū)域發(fā)展規(guī)劃。二是各區(qū)域畜牧業(yè)減排工作的開(kāi)展要因地而異、因類而異和因地制宜,并以綠色和循環(huán)發(fā)展理念指導(dǎo)發(fā)展畜牧業(yè),實(shí)現(xiàn)環(huán)境和經(jīng)濟(jì)協(xié)調(diào)發(fā)展。三是一方面加大畜牧業(yè)財(cái)政支出力度,建立支持畜牧業(yè)低碳綠色循環(huán)發(fā)展的長(zhǎng)效機(jī)制,實(shí)施積極的技術(shù)補(bǔ)貼措施,對(duì)引進(jìn)、采用和研發(fā)新技術(shù)的規(guī)?;髽I(yè)給予財(cái)政補(bǔ)貼和稅收減免;另一方面加大畜牧業(yè)固定資產(chǎn)投資力度和畜牧業(yè)從業(yè)人員培訓(xùn)投資力度,為畜牧業(yè)高質(zhì)量發(fā)展 提供支持。