陳偉 洪靜蘭 李照令 王亞平 張?zhí)祆? 耿涌

關(guān)鍵詞：生命周期評(píng)價(jià)；大豆；碳排放效率；影響因素；超效率SBM；Tobit模型

全球變暖危及地球生態(tài)安全和人類可持續(xù)發(fā)展，已經(jīng)成為全人類共同面臨的嚴(yán)重威脅［1］。農(nóng)業(yè)生產(chǎn)活動(dòng)是碳排放的重要來源，世界糧食體系的碳排放約占全球人為碳排放的三分之一［2-3］。作為全球農(nóng)業(yè)大國，中國農(nóng)業(yè)碳排放總量由1961年的2.49億t增長(zhǎng)至2018年的8.7億t［4］。若不采取有效的碳減排措施，2060年中國農(nóng)業(yè)碳排放相對(duì)于2018年將增加60%以上［5］。在“雙碳”目標(biāo)背景下，農(nóng)業(yè)面臨更加嚴(yán)峻的碳減排壓力。大豆是中國主要糧食作物之一，同時(shí)也是重要的油料作物以及優(yōu)質(zhì)蛋白的主要來源。隨著人民生活水平的日益提高和養(yǎng)殖業(yè)的快速發(fā)展，對(duì)油脂和飼用蛋白的需求不斷增加［6］。然而，中國大豆對(duì)外依存度高達(dá)85%，提高大豆自給率事關(guān)國家糧食安全。2021年農(nóng)業(yè)農(nóng)村部印發(fā)的《“十四五”全國種植業(yè)發(fā)展規(guī)劃》的通知中提出“到2025年，力爭(zhēng)大豆播種面積達(dá)到1.6億畝左右，產(chǎn)量達(dá)到2300萬t左右”［7］（1畝≈667m2）。然而，2011—2018年中國大豆生產(chǎn)單位產(chǎn)量碳排放和單位產(chǎn)值碳排放分別增長(zhǎng)了18.33%和52.42%［8］；而且種植業(yè)生產(chǎn)過程中農(nóng)業(yè)機(jī)械、農(nóng)藥、化肥等替代勞動(dòng)力使用的增加，將進(jìn)一步提高種植業(yè)碳排放強(qiáng)度。由此可見，加快推動(dòng)大豆低碳生產(chǎn)行動(dòng)責(zé)無旁貸。大豆碳減排要求在既定投入要素下保證大豆產(chǎn)值并減少碳排放。然而，中國地域發(fā)展不平衡，不同地區(qū)大豆種植的自然條件、生產(chǎn)方式、技術(shù)水平等具有較強(qiáng)的異質(zhì)性，不同地區(qū)大豆生產(chǎn)碳排放及碳排放效率存在差異。因此，有必要基于不同地區(qū)要素稟賦結(jié)構(gòu)的異質(zhì)性，考察大豆生產(chǎn)碳排放及碳排放效率。本研究在客觀核算中國大豆生產(chǎn)全生命周期碳排放量的基礎(chǔ)上測(cè)度碳排放效率，探索宏觀層面的影響因素并探尋優(yōu)化路徑，為推進(jìn)大豆乃至農(nóng)業(yè)低碳發(fā)展提供有益支持。

1文獻(xiàn)綜述

從現(xiàn)有文獻(xiàn)來看，農(nóng)業(yè)碳排放的研究成果較為豐富，涵蓋農(nóng)業(yè)碳排放測(cè)度、農(nóng)業(yè)碳排放效率測(cè)算和農(nóng)業(yè)碳排放效率影響因素探究等方面。①學(xué)者們從不同方面測(cè)算了農(nóng)業(yè)碳排放。研究尺度上，金書秦等［4］和田云等［9］針對(duì)國家、吳賢榮等［10］針對(duì)省域、尚杰等［11］針對(duì)市域和李遠(yuǎn)玲等［12］針對(duì)縣域的農(nóng)業(yè)碳排放問題開展了大量研究。研究?jī)?nèi)容上，既有種植業(yè)［13］、畜牧業(yè)［14］和漁業(yè)［15］等產(chǎn)業(yè)碳排放研究，也有針對(duì)單一作物如水稻［16］、小麥［17］、玉米［17］、油菜籽［18］和大豆［19］等的研究。②農(nóng)業(yè)碳排放效率作為衡量農(nóng)業(yè)產(chǎn)出水平的指標(biāo)受到廣泛關(guān)注。吳賢榮等［10］運(yùn)用DEA模型、陳陽等［20］運(yùn)用超效率SBM模型對(duì)不同尺度農(nóng)業(yè)和不同種類農(nóng)作物［21-22］的碳排放效率進(jìn)行測(cè)度，研究尺度涉及國家［23］、省域［24］和縣域［25］等。③在測(cè)度碳排放效率的基礎(chǔ)上對(duì)其影響因素展開探究。吳賢榮等［10］運(yùn)用Tobit模型、尚莉媛等［26］運(yùn)用空間杜賓模型，探討自然、經(jīng)濟(jì)和社會(huì)等因素對(duì)農(nóng)業(yè)碳排放效率的影響，評(píng)價(jià)指標(biāo)涉及產(chǎn)業(yè)結(jié)構(gòu)［24］、農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平［10］、勞動(dòng)力受教育程度［27］、農(nóng)村用電量［24］和農(nóng)作物受災(zāi)程度［23］等。除上述研究外，田云等［28］圍繞農(nóng)業(yè)碳排放效率與經(jīng)濟(jì)增長(zhǎng)的耦合關(guān)系、田云等［29］圍繞農(nóng)業(yè)碳減排補(bǔ)償機(jī)制和王辰璇等［30］圍繞農(nóng)業(yè)碳排放效率門檻效應(yīng)等方面，對(duì)農(nóng)業(yè)碳排放相關(guān)問題展開廣泛探討。

現(xiàn)有農(nóng)業(yè)碳排放相關(guān)文獻(xiàn)和實(shí)踐研究成果較為豐富，對(duì)于理解大豆生產(chǎn)碳排放效率具有重要啟發(fā)和借鑒意義。文獻(xiàn)梳理表明：①不同地區(qū)的農(nóng)業(yè)在自然條件、生產(chǎn)方式和技術(shù)水平等方面具有較強(qiáng)的異質(zhì)性，且不同農(nóng)作物在投入產(chǎn)出方面有較大差異。因此，需要立足具體的農(nóng)作物構(gòu)建區(qū)域異質(zhì)性的作物生命周期清單，測(cè)度其碳排放量、碳排放效率及其影響因素。然而，目前全國尺度大豆生產(chǎn)碳排放研究相對(duì)較少，且鮮有關(guān)于大豆生產(chǎn)碳排放效率的系統(tǒng)研究。②農(nóng)業(yè)碳排放效率研究中碳排放測(cè)算多基于碳排放系數(shù)法，所用碳源的碳排放系數(shù)來源多樣，而碳源的碳排放系數(shù)受生產(chǎn)技術(shù)、管理水平等影響，具有明顯的區(qū)域異質(zhì)性。③目前尚沒有關(guān)于大豆生產(chǎn)碳排放效率影響因素的研究。因此，本研究采用生命周期評(píng)價(jià)法客觀測(cè)度2011—2020年全國及各主產(chǎn)區(qū)大豆生產(chǎn)全生命周期碳排放，繼而采用兼顧期望產(chǎn)出和非期望產(chǎn)出的超效率SBM模型測(cè)度大豆生產(chǎn)碳排放效率，并探究其時(shí)空演化態(tài)勢(shì)，而后利用Tobit模型探尋影響大豆碳排放效率變動(dòng)的主要因素，并提出促進(jìn)大豆生產(chǎn)低碳發(fā)展的對(duì)策建議。

2研究方法、變量選取與數(shù)據(jù)說明

2.1研究方法

2.1.1大豆生產(chǎn)碳排放測(cè)算方法

（1）研究功能單位與系統(tǒng)邊界。國際標(biāo)準(zhǔn)化組織（ISO）于2013年發(fā)布ISO14067，將生命周期評(píng)價(jià)確定為產(chǎn)品碳足跡核算的技術(shù)方法［31-32］，確定功能單位和系統(tǒng)邊界是開展生命周期評(píng)價(jià)的基礎(chǔ)。農(nóng)業(yè)生命周期評(píng)價(jià)研究通常以單位農(nóng)產(chǎn)品產(chǎn)量為功能單位，多采取“搖籃到大門”的系統(tǒng)邊界［33］。因此，本研究的功能單位為1t大豆產(chǎn)品，系統(tǒng)邊界為“搖籃到大門”，即測(cè)算大豆從播種到收獲全過程各項(xiàng)投入造成的直接和間接溫室氣體排放，具體包括：①生產(chǎn)及運(yùn)輸各項(xiàng)農(nóng)資（化肥、種子和農(nóng)藥等）產(chǎn)生的間接碳排放。②柴油消耗過程的直接碳排放。

（2）清單分析。生命周期清單分析是對(duì)產(chǎn)品、工藝或活動(dòng)整個(gè)生命周期階段的能源等資源消耗和向環(huán)境的排放進(jìn)行數(shù)據(jù)量化分析［34］。大豆種植數(shù)據(jù)收集分為兩個(gè)部分：一是實(shí)景流程（即表層清單）數(shù)據(jù)收集，包括播種、灌溉、施肥等過程的資源投入，以及向環(huán)境中的排放（例如柴油消耗導(dǎo)致的直接碳排放）；二是背景流程（即背景清單）數(shù)據(jù)收集，即投入資源（電力、種子和化肥等）的生產(chǎn)數(shù)據(jù)。大豆生產(chǎn)碳排放相關(guān)生命周期表層清單見表1。

（3）碳排放計(jì)算。SimaPro軟件由荷蘭Leiden大學(xué)環(huán)境科學(xué)中心開發(fā)，是現(xiàn)階段常用的生命周期影響評(píng)價(jià)軟件，該軟件集成了ReCiPe等多種生命周期評(píng)價(jià)方法。本研究選用SimaPro軟件ReCiPe2016模型計(jì)算大豆生產(chǎn)全生命周期溫室氣體排放（以CO2為當(dāng)量進(jìn)行表征，本研究測(cè)算結(jié)果統(tǒng)一用碳排放表示）。

2.1.2基于非期望產(chǎn)出的超效率SBM模型

DEA是評(píng)價(jià)決策單元有效性的常用模型。傳統(tǒng)DEA模型無法測(cè)度非期望產(chǎn)出對(duì)效率的影響，且忽視了對(duì)多個(gè)有效決策單元的對(duì)比［35］。因此，Tone對(duì)DEA模型進(jìn)行改進(jìn)，提出了基于非期望產(chǎn)出的非徑向、非角度超效率DEA方法，即基于非期望產(chǎn)出的超效率SBM（Slacks?BasedMeasure，SBM）模型［36］。超效率SBM模型既能將松弛變量納入考慮，又能克服傳統(tǒng)SBM模型難以區(qū)分效率值大于1的問題，在碳排放效率測(cè)算方面得到廣泛應(yīng)用［1］。

2.1.3Tobit模型

在對(duì)大豆種植碳排放效率進(jìn)行測(cè)算后，為進(jìn)一步探究影響大豆碳排放效率變動(dòng)的深層原因，有必要開展大豆種植碳排放效率影響因素分析。但由于超效率SBM模型測(cè)算出的農(nóng)業(yè)碳排放效率值是受限因變量（以0為下限），而Tobit模型能夠較好地規(guī)避最小二乘法導(dǎo)致的估計(jì)偏差［37］，因此該方法在碳排放效率影響因素探究中得到廣泛使用。參考張展等［38］的研究，選擇Tobit模型來解決這類受限因變量的問題。該模型已經(jīng)被學(xué)界熟知，在農(nóng)業(yè)碳排放效率影響因素分析中得到了較廣泛應(yīng)用［10，24，38］，為此對(duì)其公式不進(jìn)行詳細(xì)介紹。

2.2指標(biāo)選取

2.2.1碳排放效率指標(biāo)選取

本研究選取勞動(dòng)、土地和資本三種投入要素，各指標(biāo)的數(shù)據(jù)特征見表2。將大豆生產(chǎn)用工天數(shù)（農(nóng)戶家庭成員勞動(dòng)時(shí)長(zhǎng)和雇傭工人勞動(dòng)時(shí)長(zhǎng)之和）作為勞動(dòng)投入指標(biāo)；以大豆播種面積為土地投入指標(biāo)；以大豆生產(chǎn)所需物質(zhì)與服務(wù)費(fèi)用為資本投入指標(biāo)，包括與大豆種植相關(guān)的種子費(fèi)、化肥費(fèi)、農(nóng)藥費(fèi)、租賃作業(yè)費(fèi)、工具材料費(fèi)和修理維護(hù)費(fèi)等直接費(fèi)用，以及保險(xiǎn)費(fèi)和管理費(fèi)等間接費(fèi)用，并將所需物質(zhì)與服務(wù)費(fèi)用調(diào)整為以2011年為基期的不變價(jià)。

本研究將產(chǎn)出要素指標(biāo)分為期望產(chǎn)出和非期望產(chǎn)出兩種指標(biāo)。以大豆產(chǎn)值為期望產(chǎn)出，并將其調(diào)整為2011年的不變價(jià)格；將大豆種植碳排放作為非期望產(chǎn)出，并采用生命周期評(píng)價(jià)法評(píng)估大豆生產(chǎn)碳排放。

2.2.2影響因素指標(biāo)選取

厘清碳排放效率影響因素是把握優(yōu)化方向的重要前提。根據(jù)相關(guān)文獻(xiàn)中對(duì)農(nóng)業(yè)碳排放效率影響的研究經(jīng)驗(yàn)啟示，研究中重點(diǎn)關(guān)注種植結(jié)構(gòu)、農(nóng)戶文化水平、機(jī)械化水平等與種植業(yè)生產(chǎn)相關(guān)的因素對(duì)碳排放效率的影響，選取以下變量探討大豆種植碳排放效率的優(yōu)化路徑。

（1）作物種植結(jié)構(gòu)。種植結(jié)構(gòu)反映地區(qū)某類作物種植的比例，種植結(jié)構(gòu)的變化會(huì)導(dǎo)致農(nóng)業(yè)投入要素結(jié)構(gòu)的變化，進(jìn)而影響作物生產(chǎn)的碳排放效率［39］。參考張利國等［40］的研究，選用大豆播種面積占糧食總播種面積的比例表征作物種植結(jié)構(gòu)，預(yù)計(jì)該指標(biāo)與大豆生產(chǎn)碳排放效率呈正相關(guān)。

（2）農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)。農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)能夠體現(xiàn)第一產(chǎn)業(yè)各生產(chǎn)部門在整個(gè)農(nóng)業(yè)系統(tǒng)中所占比例。參考吳賢榮等［10］的研究，選擇農(nóng)業(yè)產(chǎn)值占農(nóng)林牧漁總產(chǎn)值的比重來表示農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，預(yù)計(jì)該指標(biāo)與大豆生產(chǎn)碳排放效率呈負(fù)相關(guān)。

（3）農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平。農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平是反映地區(qū)農(nóng)業(yè)發(fā)展水平和農(nóng)民生活水平的重要指標(biāo)［24］。農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平高的地區(qū)農(nóng)用物資投入量相對(duì)較高，但同時(shí)農(nóng)業(yè)現(xiàn)代化進(jìn)程較快，有利于農(nóng)業(yè)碳排放效率提升［24］。借鑒侯孟陽等［39］的研究，選用農(nóng)村人均可支配收入作為農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平的表征，預(yù)計(jì)該指標(biāo)與大豆生產(chǎn)碳排放效率呈正相關(guān)。

（4）財(cái)政支農(nóng)水平。財(cái)政支農(nóng)水平反映地區(qū)財(cái)政支農(nóng)的力度。財(cái)政支農(nóng)增加農(nóng)民收入，幫助提高農(nóng)業(yè)生產(chǎn)過程中的物資投入，改進(jìn)生產(chǎn)設(shè)備和生產(chǎn)方式，有利于大豆生產(chǎn)碳排放效率的提高。借鑒胡青江等［41］的研究，選擇農(nóng)林水支出占財(cái)政總支出的比重作為描述財(cái)政支農(nóng)水平的指標(biāo)，預(yù)計(jì)該指標(biāo)與大豆生產(chǎn)碳排放效率呈正相關(guān)。

（5）農(nóng)業(yè)機(jī)械化水平。一方面，農(nóng)業(yè)機(jī)械利用會(huì)消耗柴油等物資，造成碳排放；另一方面，農(nóng)業(yè)機(jī)械能夠通過提升農(nóng)業(yè)產(chǎn)出效率以減少農(nóng)產(chǎn)品的碳排放，進(jìn)而提高碳排放效率［42］。因難以預(yù)測(cè)農(nóng)業(yè)機(jī)械化水平對(duì)大豆生產(chǎn)碳排放效率的作用方向，所以需進(jìn)行實(shí)證檢驗(yàn)。參考吳昊玥等［43］的研究，以地均農(nóng)業(yè)機(jī)械動(dòng)力對(duì)農(nóng)業(yè)機(jī)械化水平進(jìn)行表征。

（6）農(nóng)村勞動(dòng)力受教育程度。農(nóng)村勞動(dòng)力的受教育水平一方面影響勞動(dòng)力對(duì)低碳農(nóng)業(yè)的認(rèn)識(shí)，另一方面影響勞動(dòng)力對(duì)新型農(nóng)業(yè)生產(chǎn)技術(shù)的運(yùn)用，進(jìn)而影響農(nóng)業(yè)碳排放效率。但馬大來［44］的研究結(jié)果表明，當(dāng)勞動(dòng)力以初、高中文化素質(zhì)為主時(shí)，將對(duì)農(nóng)業(yè)碳排放效率產(chǎn)生負(fù)向影響。該變量可能的作用方向難以確定，因此需進(jìn)行實(shí)證檢驗(yàn)。參考田云等［27］的研究，選用初中及以上文化程度勞動(dòng)者占比對(duì)農(nóng)村勞動(dòng)力受教育程度進(jìn)行表征。

（7）農(nóng)業(yè)技術(shù)水平。在農(nóng)業(yè)生產(chǎn)活動(dòng)中，電力能夠提供動(dòng)力，替代部分柴油等能源。農(nóng)業(yè)產(chǎn)業(yè)化經(jīng)營模式的改變會(huì)導(dǎo)致農(nóng)村用電量的增加［45］。農(nóng)村用電量的增加意味著農(nóng)村農(nóng)用機(jī)械設(shè)備使用的增加和技術(shù)水平的進(jìn)步［46］。參考吳妍［47］的研究，選擇農(nóng)村人均用電量這一指標(biāo)探究農(nóng)業(yè)技術(shù)水平對(duì)大豆碳排放效率的影響。農(nóng)業(yè)現(xiàn)代化程度越高，農(nóng)業(yè)產(chǎn)出水平也會(huì)越高，預(yù)計(jì)該指標(biāo)與大豆生產(chǎn)碳排放效率呈正相關(guān)。

2.3研究區(qū)域與數(shù)據(jù)來源

鑒于數(shù)據(jù)可得性和連續(xù)性，本研究選取2011—2020年為研究年限，基于國家尺度以及中國10個(gè)大豆主產(chǎn)區(qū)（河北、山西、內(nèi)蒙古、遼寧、吉林、黑龍江、安徽、山東、河南和陜西）開展大豆生產(chǎn)碳排放效率相關(guān)研究。2020年，10個(gè)主產(chǎn)區(qū)大豆總播種面積和產(chǎn)量分別占全國大豆播種面積和產(chǎn)量的80.82%和79.16%，具有良好的代表性。

研究所需的大豆生產(chǎn)各類碳排放源、大豆生產(chǎn)勞動(dòng)投入、資本投入等數(shù)據(jù)來自《全國農(nóng)產(chǎn)品成本收益資料匯編》［48］；碳排放測(cè)算所需的碳排放系數(shù)等數(shù)據(jù)來自基于過程的中國生命周期清單數(shù)據(jù)庫［8］；大豆產(chǎn)量、大豆產(chǎn)值、農(nóng)村人均可支配收入、大豆播種面積、糧食總播種面積、耕地面積、農(nóng)業(yè)機(jī)械總動(dòng)力、農(nóng)業(yè)產(chǎn)值、農(nóng)林牧漁業(yè)總產(chǎn)值和農(nóng)村用電量數(shù)據(jù)來自《中國農(nóng)村統(tǒng)計(jì)年鑒》［49］；農(nóng)林水支出和財(cái)政總支出來自各省統(tǒng)計(jì)年鑒；農(nóng)村居民消費(fèi)價(jià)格指數(shù)來自國家統(tǒng)計(jì)局［50］；農(nóng)村初中及以上文化程度勞動(dòng)者數(shù)據(jù)來自《中國人口和就業(yè)統(tǒng)計(jì)年鑒》［51］。

3研究結(jié)果分析

3.1大豆生產(chǎn)碳排放時(shí)空演變分析

3.1.1中國大豆生產(chǎn)歷年碳排放分析

本研究基于生命周期評(píng)價(jià)法，測(cè)度中國大豆生產(chǎn)全生命周期碳排放（表3）。結(jié)果表明，大豆生產(chǎn)碳排放總量由2011年的4.07×107t增至2020年的6.02×107t，年均增長(zhǎng)率為4.43%；2011—2012年碳排放量稍有增加，年均增長(zhǎng)率為0.20%；2013—2015年碳排放量逐年降低，年均遞減率為4.49%；2016—2020年碳排放量呈上升態(tài)勢(shì)，由4.60×107t增加到6.02×107t，年均增長(zhǎng)率為6.96%。

生產(chǎn)1t大豆的碳排放由2011年的2.81t增至2020年的3.07t，提高了9.25%。2011—2013年中國大豆生產(chǎn)單位產(chǎn)量碳排放呈上升趨勢(shì)；隨后生產(chǎn)1t大豆產(chǎn)生的碳排放在2.99～3.55t間波動(dòng)，其中2016年每噸大豆生產(chǎn)產(chǎn)生的碳排放達(dá)到最大值（3.55t）。羅燕等［52］研究指出每噸大豆生產(chǎn)碳排放為1.53～2.17t，Castanheira等［53］研究指出每噸大豆種植碳排放在0.1～17.8t間變化。數(shù)據(jù)比較表明，本研究結(jié)果在已有研究結(jié)果范圍內(nèi)。通過表3可知，每萬元產(chǎn)值大豆碳排放由2011年的5.64t增至2020年的7.48t，年均增長(zhǎng)率為3.20%。2011—2016年中國大豆生產(chǎn)單位產(chǎn)值碳排放總體呈上升趨勢(shì)；2016—2020年呈波動(dòng)下降趨勢(shì)。其中，2017—2018年出現(xiàn)短暫上升，并于2018年達(dá)到最大值（10.17t/萬元）。大豆生產(chǎn)的碳排放強(qiáng)度（單位面積碳排放量）由2011年的5.16t/hm2增至2020年的6.09t/hm2，年均增長(zhǎng)率為1.85%。2011—2013年中國大豆生產(chǎn)碳排放強(qiáng)度呈上升趨勢(shì)；隨后在5.34～6.38t/hm2間波動(dòng)。數(shù)據(jù)分析表明，2011—2020年中國大豆生產(chǎn)碳排放總體呈上升趨勢(shì)，但中國大豆生產(chǎn)尚未實(shí)現(xiàn)碳達(dá)峰。

3.1.2大豆生產(chǎn)主產(chǎn)區(qū)碳排放時(shí)空演變分析

大豆生產(chǎn)主產(chǎn)區(qū)碳排放測(cè)度結(jié)果表明大豆生產(chǎn)碳排放時(shí)空差異明顯，各主產(chǎn)區(qū)大豆生產(chǎn)碳排放總量（表4）、單位產(chǎn)值碳排放（表5）和碳排放強(qiáng)度（表6）均存在明顯的差異。就碳排放總量來看，排在首位的黑龍江2020年大豆生產(chǎn)碳排放高達(dá)34.84百萬t，占全國大豆生產(chǎn)碳排放的57.91%；位居第2的內(nèi)蒙古2020年大豆生產(chǎn)碳排放為7.45百萬t，占全國大豆生產(chǎn)碳排放的12.38%。黑龍江、內(nèi)蒙古是中國大豆生產(chǎn)最主要的產(chǎn)區(qū)，大豆種植面積和產(chǎn)量高于其他地區(qū)，因此，碳排放量高。研究期內(nèi)，除吉林和河北2020年大豆生產(chǎn)碳排放總量較2011年分別下降20.86%和5.37%外，其余各主產(chǎn)區(qū)大豆生產(chǎn)碳排放總量均呈增長(zhǎng)趨勢(shì)。這可能是由于吉林、河北糧食作物結(jié)構(gòu)發(fā)生變化，大豆種植面積和產(chǎn)量減少［54-55］，導(dǎo)致其碳排放量減少；而山東、黑龍江和內(nèi)蒙古等地在單產(chǎn)保持增長(zhǎng)的同時(shí)種植面積也在擴(kuò)大［56］，因此碳排放總量呈上升趨勢(shì)。其中，山東2020年大豆生產(chǎn)碳排放總量較2011年增長(zhǎng)了119.81%，黑龍江（94.97%）、內(nèi)蒙古（52.69%）分列第二、三位。

碳排放總量與生產(chǎn)規(guī)模密切相關(guān)，而單位產(chǎn)值碳排放和碳排放強(qiáng)度更能客觀反映區(qū)域碳排放水平。從表5可以看出，2020年黑龍江大豆每萬元產(chǎn)值碳排放高達(dá)9.07t，內(nèi)蒙古（8.99t）、陜西（8.63t）分列第二、三位。研究期內(nèi)，單位產(chǎn)值大豆碳排放最高達(dá)到27.34t/萬元（2016年，內(nèi)蒙古），最低為1.77t/萬元（2012年，安徽；2015年，河南）。吉林單位產(chǎn)值大豆碳排放2011年6.03t/萬元，2020年5.41t/萬元，下降幅度為10.27%，其余主產(chǎn)區(qū)單位產(chǎn)值大豆碳排放均呈上升趨勢(shì)。其中，大豆產(chǎn)值最高的黑龍江，單位產(chǎn)值大豆碳排放由2011年的6.10t/萬元增至2020年的9.07t/萬元，增幅為48.58%；山西單位產(chǎn)值大豆碳排放由2011年的2.86t/萬元增至2020年的8.25t/萬元，增幅為187.84%，是增幅最大的地區(qū)；山東單位產(chǎn)值大豆碳排放由2011年的2.86t/萬元增至2020年的4.90t/萬元，增幅為71.30%；河南、陜西、遼寧、安徽、河北和內(nèi)蒙古單位產(chǎn)值大豆碳排放增幅分別為93.84%、65.18%、56.56%、32.76%、15.68%和5.87%。

就碳排放強(qiáng)度而言，2020年河北大豆生產(chǎn)碳排放強(qiáng)度高達(dá)9.64t/hm2，遼寧（7.60t/hm2）、黑龍江（7.21t/hm2）分列第二、三位（表6）。研究期內(nèi)，碳排放強(qiáng)度最高達(dá)到9.64t/hm2（2018年，2020年，河北），最低為1.49t/hm2（2015年，河南）。吉林大豆生產(chǎn)碳排放強(qiáng)度2011年6.68t/hm2，2020年5.02t/hm2，內(nèi)蒙古大豆生產(chǎn)碳排放強(qiáng)度2011年7.09t/hm2，2020年6.20t/hm2，降幅分別為24.86%和12.63%，其余主產(chǎn)區(qū)大豆生產(chǎn)碳排放強(qiáng)度均呈上升趨勢(shì)。產(chǎn)量最高的黑龍江，大豆生產(chǎn)碳排放強(qiáng)度2011年5.58t/hm2，2020年7.21t/hm2，提高了29.19%；山西大豆生產(chǎn)碳排放強(qiáng)度2011年2.96t/hm2，2020年7.17t/hm2，增幅為142.36%；山東大豆生產(chǎn)碳排放強(qiáng)度2011年3.42t/hm2，2020年6.22t/hm2，增幅為81.95%；河南大豆生產(chǎn)碳排放強(qiáng)度2011年2.03t/hm2，2020年3.61t/hm2，增幅為77.84%；安徽、河北和陜西大豆生產(chǎn)碳排放強(qiáng)度也均有增加，增幅分別為47.94%、43.90%和40.56%。

結(jié)合2011年和2020年大豆種植碳排放量、單位產(chǎn)值碳排放和碳排放強(qiáng)度差異，借助ArcGIS軟件將各主產(chǎn)區(qū)碳排放情況劃分為5種類型（表7），分別為“高高高”型，即大豆種植碳排放高、單位產(chǎn)值碳排放高、碳排放強(qiáng)度高的地區(qū)；“低高高”型，即大豆種植碳排放低、單位產(chǎn)值碳排放高、碳排放強(qiáng)度高的地區(qū)；“低低高”型，即大豆種植碳排放低、單位產(chǎn)值碳排放低、碳排放強(qiáng)度高的地區(qū)；“低高低”型，即大豆種植碳排放低、單位產(chǎn)值碳排放高、碳排放強(qiáng)度低的地區(qū)；“低低低”型，即大豆種植碳排放低、單位產(chǎn)值碳排放低、碳排放強(qiáng)度低的地區(qū)。

由表7可知，黑龍江屬于“高高高”型地區(qū)，大豆種植碳排放總量高，單位產(chǎn)值碳排放和碳排放強(qiáng)度同樣較高；內(nèi)蒙古、遼寧、河北屬于“低高高”型地區(qū)，大豆種植碳排放總量較低，但單位產(chǎn)值碳排放和碳排放強(qiáng)度較高；由于單位產(chǎn)值碳排放和碳排放強(qiáng)度較高，上述兩類地區(qū)的大豆種植碳排放壓力較大。安徽和河南屬于“低低低”型地區(qū)，大豆種植碳排放總量低，單位產(chǎn)值碳排放和碳排放強(qiáng)度同樣較低，這類地區(qū)大豆種植碳排放壓力相對(duì)較小。吉林由2011年的“低高高”型地區(qū)轉(zhuǎn)變?yōu)?020年的“低低低”型，單位產(chǎn)值碳排放和碳排放強(qiáng)度有所降低，大豆種植碳排放壓力呈下降趨勢(shì)。陜西由2011年的“低高低”型轉(zhuǎn)變?yōu)?020年的“低高高”型；山東由2011年的“低低低”型轉(zhuǎn)變?yōu)?020年的“低低高”型；山西由2011年的“低低低”型轉(zhuǎn)變?yōu)?020年的“低高高”型；上述地區(qū)大豆種植碳排放壓力呈上升趨勢(shì)。

3.2大豆生產(chǎn)碳排放效率分析

3.2.1中國大豆生產(chǎn)碳排放效率分析

運(yùn)用基于非期望產(chǎn)出的超效率SBM模型測(cè)算中國大豆生產(chǎn)碳排放效率，結(jié)果如圖1所示。總體來看，2011—2020年中國大豆生產(chǎn)碳排放效率呈先降后升的演變趨勢(shì)，平均值為0.78。其中，2020年大豆生產(chǎn)碳排放效率為1.00，較最小值0.55（2016年）提升了0.46，較最大值1.05（2012年）降低了0.05。2011—2012年，大豆生產(chǎn)碳排放效率由0.90上升至1.05。該階段單位產(chǎn)值大豆的勞動(dòng)投入下降了13.98%，同時(shí)土地投入和資本投入也分別下降11.44%和1.40%。2012—2016年，碳排放效率大幅下降。該階段單位產(chǎn)值大豆生產(chǎn)所投入的勞動(dòng)、土地和資本分別上升39.35%、62.93%和48.32%，大豆生產(chǎn)碳排放總量上升12.64%，而大豆產(chǎn)值呈顯著的下降趨勢(shì)。2016—2017年，大豆生產(chǎn)碳排放效率有所提高。該階段大豆產(chǎn)值增幅為29.71%，高于大豆生產(chǎn)碳排放增幅（0.15%）。2017—2018年，大豆生產(chǎn)碳排放效率明顯降低。該階段生產(chǎn)單位產(chǎn)值大豆所投入的勞動(dòng)、土地和資本分別提高5.69%、15.80%和14.74%，大豆生產(chǎn)碳排放量增加709.38萬t，大豆生產(chǎn)的物資利用效率再度下降。2018—2020年，碳排放效率顯著提高。該階段大豆生產(chǎn)的勞動(dòng)、土地和資本投入分別降低了21.31%、23.64%和27.43%，而大豆生產(chǎn)碳排放量?jī)H增加了13.22%。

3.2.2大豆生產(chǎn)主產(chǎn)區(qū)的碳排放效率

大豆主產(chǎn)區(qū)碳排放效率測(cè)度結(jié)果表明，大豆生產(chǎn)碳排放效率空間差異明顯（表8）。從多年均值來看，內(nèi)蒙古大豆生產(chǎn)碳排放效率值最高（1.15），山西最低（0.50）；大豆產(chǎn)量最高的黑龍江碳排放效率為0.93。除吉林、山西碳排放效率年均變化率為負(fù)值外，其余地區(qū)大豆生產(chǎn)碳排放效率年均變化率均為正。其中，內(nèi)蒙古、河北初始值較低但增速快，年均增長(zhǎng)率為7.40%和3.68%，屬于低位高速增長(zhǎng)類型；河南和安徽碳排放效率初始值較高且增速較快，年均增長(zhǎng)率分別為15.14%和8.12%，屬于高位高速增長(zhǎng)類型。此外，山西、陜西大豆生產(chǎn)歷年碳排放效率與最優(yōu)前沿面有較大差距，屬于碳排放效率低值地區(qū)，其碳排放效率低的原因在于其“有效產(chǎn)出”能力不強(qiáng)，從而使其碳排放效率與前沿面距離較遠(yuǎn)，山西和陜西在大豆低碳種植方面具有較大的提升空間。

3.3大豆生產(chǎn)碳排放效率的影響因素

研究運(yùn)用Tobit模型對(duì)大豆生產(chǎn)碳排放效率影響因素進(jìn)行回歸分析。在模型估計(jì)前，首先對(duì)變量進(jìn)行多重共線性檢驗(yàn)，結(jié)果顯示，最大的方差膨脹因子（VIF）值為2.72，小于10，故不存在共線性問題。為避免變量非平穩(wěn)可能導(dǎo)致的偽回歸，對(duì)面板數(shù)據(jù)進(jìn)行單位根檢驗(yàn)，結(jié)果顯示各變量均通過檢驗(yàn)。由于固定效應(yīng)的Tobit模型找不到個(gè)體異質(zhì)性的充分統(tǒng)計(jì)量，從而無法進(jìn)行條件最大似然估計(jì)［57］，因此本研究考慮隨機(jī)效應(yīng)Tobit模型和混合Tobit模型。在隨機(jī)效應(yīng)Tobit中，由于LR檢驗(yàn)無法拒絕原假設(shè)，最終選擇混合Tobit模型，得出結(jié)果見表9。

作物種植結(jié)構(gòu)影響。大豆播種面積占比對(duì)大豆生產(chǎn)碳排放效率具有顯著正向影響，但相關(guān)程度較弱?？赡艿脑蚴牵啾扔衩?、小麥和水稻等糧食作物，大豆價(jià)格較高，致使大豆單位產(chǎn)值所產(chǎn)生的碳排放較其他糧食作物的單位產(chǎn)值碳排放低。因此，基于目前的種植技術(shù)和生產(chǎn)條件，提高大豆播種面積占糧食總播種面積的比重能夠使大豆碳排放效率實(shí)現(xiàn)小幅提高。

農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)影響。農(nóng)業(yè)產(chǎn)值占農(nóng)林牧漁總產(chǎn)值比重對(duì)大豆生產(chǎn)碳排放效率具有顯著的負(fù)向影響，即農(nóng)業(yè)產(chǎn)值占農(nóng)林牧漁總產(chǎn)值比重越高，大豆生產(chǎn)碳排放效率越低。這可能是因?yàn)橄噍^于第一產(chǎn)業(yè)的其他部門，農(nóng)業(yè)生產(chǎn)活動(dòng)中需要投入更多的農(nóng)用物資，從而導(dǎo)致碳排放量較高，碳排放效率相應(yīng)較低。

農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平影響。農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平對(duì)大豆生產(chǎn)碳排放效率呈負(fù)向影響，但未通過顯著性檢驗(yàn)?？赡艿脑蚴牵r(nóng)民收入的提升促使更多資本和能源投入到農(nóng)業(yè)生產(chǎn)中，從而導(dǎo)致碳排放量增加，碳排放效率降低。但由于中國大豆的集約化生產(chǎn)方式還未成型，因此規(guī)模效應(yīng)難以體現(xiàn)，所以該指標(biāo)對(duì)大豆生產(chǎn)碳排放效率呈現(xiàn)不顯著的負(fù)向影響。

財(cái)政支農(nóng)水平影響。農(nóng)林水支出占財(cái)政總支出的比重對(duì)大豆生產(chǎn)碳排放效率具有顯著的正向影響?？赡艿脑蚴牵ㄟ^財(cái)政支農(nóng)政策加強(qiáng)對(duì)農(nóng)業(yè)生產(chǎn)的支持，在一定程度上提升農(nóng)業(yè)管理水平，從而提高大豆碳排放效率。但李勇［58］指出，由于惠農(nóng)補(bǔ)貼存在對(duì)象錯(cuò)配、范圍和額度有限等問題，導(dǎo)致財(cái)政支農(nóng)對(duì)農(nóng)業(yè)生產(chǎn)者的支持力度不夠。因此，應(yīng)加大財(cái)政補(bǔ)貼力度，并提升精準(zhǔn)管理的保障能力，促進(jìn)大豆低碳生產(chǎn)。

農(nóng)業(yè)機(jī)械化水平影響。農(nóng)業(yè)機(jī)械化水平對(duì)大豆生產(chǎn)碳排放效率具有顯著正向影響，即農(nóng)業(yè)機(jī)械發(fā)展水平越高，大豆生產(chǎn)碳排放效率越高。農(nóng)用機(jī)械的使用雖然導(dǎo)致溫室氣體排放，但農(nóng)業(yè)機(jī)械化生產(chǎn)能提高物資的利用效率和作物的產(chǎn)出效率，從而達(dá)到提高大豆生產(chǎn)碳排放效率的目的。

農(nóng)村勞動(dòng)力受教育程度影響。農(nóng)村勞動(dòng)力受教育程度對(duì)大豆生產(chǎn)碳排放效率具有顯著的負(fù)向影響，該結(jié)果表明，初中及以上文化水平的農(nóng)村勞動(dòng)力占比越高，大豆生產(chǎn)碳排放效率反而越低。其成因可能是先進(jìn)的農(nóng)業(yè)生產(chǎn)技術(shù)一般需要大專及以上學(xué)歷的高素質(zhì)勞動(dòng)力［27］。然而目前從事農(nóng)業(yè)生產(chǎn)的勞動(dòng)力以初、高中學(xué)歷為主，相比較而言，該部分勞動(dòng)力更傾向于選擇“高投入、高產(chǎn)出”的農(nóng)業(yè)生產(chǎn)模式［27］。因此，農(nóng)村勞動(dòng)力受教育程度對(duì)大豆生產(chǎn)碳排放效率未能產(chǎn)生正向影響。

農(nóng)業(yè)技術(shù)水平。農(nóng)村人均用電量對(duì)大豆生產(chǎn)碳排放效率具有顯著的正向影響。這可能是由于電力農(nóng)機(jī)代替部分柴油農(nóng)機(jī)能夠減少碳排放，同時(shí)在農(nóng)業(yè)生產(chǎn)中，各類機(jī)械的使用能夠提高生產(chǎn)效率，在一定程度上促進(jìn)大豆生產(chǎn)碳排放效率的提高。

4結(jié)論與討論

4.1結(jié)論

全球氣候變暖背景下，農(nóng)業(yè)部門減排與低碳生產(chǎn)是實(shí)現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展的重要途徑，也是中國實(shí)現(xiàn)碳中和目標(biāo)的必要舉措。中國正處于實(shí)現(xiàn)碳達(dá)峰目標(biāo)的關(guān)鍵時(shí)期，各行業(yè)都應(yīng)加快低碳發(fā)展。本研究基于生命周期評(píng)價(jià)法測(cè)度2011—2020年中國大豆生產(chǎn)全生命周期碳排放，采用含有非期望產(chǎn)出的超效率SBM模型測(cè)算2011—2020年中國大豆生產(chǎn)的碳排放效率，并對(duì)不同時(shí)期、不同地區(qū)的效率呈現(xiàn)結(jié)果進(jìn)行分析，最后運(yùn)用Tobit模型分析影響大豆生產(chǎn)碳排放效率的關(guān)鍵因素，主要結(jié)論如下。

（1）2011—2020年中國大豆生產(chǎn)碳排放總量整體呈上升趨勢(shì)，年均增長(zhǎng)率為4.43%。單位產(chǎn)值碳排放和碳排放強(qiáng)度年均增長(zhǎng)率分別為3.20%和1.85%，亟待推進(jìn)大豆低碳化生產(chǎn)。大豆主產(chǎn)區(qū)單位產(chǎn)值碳排放、碳排放強(qiáng)度以及碳排放總量均存在明顯差異。如吉林單位產(chǎn)值大豆碳排放由2011年的6.03t/萬元降低到2020年的5.41t/萬元；而產(chǎn)值最高的黑龍江，單位產(chǎn)值大豆碳排放則由2011年的6.10t/萬元增至2020年的9.07t/萬元。此外，2017—2020年中國大豆生產(chǎn)碳排放持續(xù)增加，中國大豆生產(chǎn)尚未實(shí)現(xiàn)碳達(dá)峰。

（2）大豆生產(chǎn)碳排放效率測(cè)算結(jié)果表明，2011—2020年中國大豆生產(chǎn)碳排放效率呈先降后升的演變趨勢(shì)，研究期內(nèi)碳排放效率整體提高了0.10。大豆主產(chǎn)區(qū)大豆生產(chǎn)碳排放效率空間差異明顯，山西和陜西大豆生產(chǎn)歷年碳排放效率與最優(yōu)前沿面有較大差距，具有較大的提升空間。

（3）運(yùn)用Tobit模型對(duì)大豆生產(chǎn)碳排放效率影響因素進(jìn)行分析發(fā)現(xiàn)，農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)和農(nóng)村勞動(dòng)力受教育程度對(duì)大豆生產(chǎn)碳排放效率有顯著負(fù)向影響；作物種植結(jié)構(gòu)、財(cái)政支農(nóng)水平、農(nóng)業(yè)機(jī)械化水平和農(nóng)業(yè)技術(shù)水平對(duì)大豆生產(chǎn)碳排放效率有顯著正向作用；農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平對(duì)大豆生產(chǎn)碳排放效率作用不顯著。

4.2 討論與啟示

（1）中國大豆生產(chǎn)的碳排放總量、單位產(chǎn)值碳排放和碳排放強(qiáng)度均呈上升趨勢(shì)，且不同地區(qū)大豆生產(chǎn)碳排放特征存在較大差異。因此，因地制宜制定大豆碳減排政策對(duì)于實(shí)現(xiàn)農(nóng)業(yè)領(lǐng)域碳達(dá)峰具有重要意義。各地區(qū)應(yīng)根據(jù)當(dāng)?shù)氐淖匀毁Y源稟賦、技術(shù)和資本等區(qū)域優(yōu)勢(shì)，采用差異化的碳減排策略，促進(jìn)地區(qū)大豆低碳生產(chǎn)。

（2）2011—2020年中國大豆生產(chǎn)碳排放效率波動(dòng)上升，主產(chǎn)區(qū)大豆生產(chǎn)碳排放效率存在區(qū)域異質(zhì)性。未來應(yīng)加強(qiáng)地區(qū)間的合作與交流，注重發(fā)揮山東等大豆生產(chǎn)高效率地區(qū)對(duì)鄰近主產(chǎn)區(qū)的輻射帶動(dòng)作用，積極引導(dǎo)河北等低效率地區(qū)借鑒高效率地區(qū)的先進(jìn)經(jīng)驗(yàn)。

（3）大豆生產(chǎn)涉及農(nóng)資及能源等使用，不可避免導(dǎo)致溫室氣體的排放。研究結(jié)果表明，農(nóng)業(yè)機(jī)械化水平與大豆生產(chǎn)碳排放效率呈顯著正相關(guān)。這是由于實(shí)施農(nóng)業(yè)機(jī)械化生產(chǎn)能夠有效提高作物生產(chǎn)率和物資利用效率，減少大豆生產(chǎn)的碳排放，從而促進(jìn)碳排放效率的提升。因此，應(yīng)加快農(nóng)業(yè)機(jī)械化進(jìn)程，提高農(nóng)資利用效率，促進(jìn)大豆生產(chǎn)低碳化水平的提升。

（4）財(cái)政支農(nóng)水平對(duì)大豆碳排放效率具有顯著的正向影響，即財(cái)政支農(nóng)的增加能夠提高大豆生產(chǎn)碳排放效率。建議加大對(duì)大豆生產(chǎn)領(lǐng)域的精準(zhǔn)補(bǔ)貼，并由各級(jí)政府部門做好預(yù)算指標(biāo)的使用和工作銜接，確保財(cái)政支農(nóng)能夠真正惠及農(nóng)民，更有力地保障大豆生產(chǎn)安全，提高農(nóng)民抵御風(fēng)險(xiǎn)的能力，從而提高大豆生產(chǎn)的碳排放效率。

本研究分析了中國大豆生產(chǎn)的碳排放效率并探究其影響因素，對(duì)于大豆低碳生產(chǎn)政策制定具有一定參考價(jià)值。不足之處在于未能充分考慮各類碳排放源的碳排放因子在時(shí)間和空間上的差異，以上不足有待于今后的深化研究。