基于多DEA-Gini的紡織業(yè)碳排放效率時(shí)空分異研究

馬玉超，張永林

基于多DEA-Gini的紡織業(yè)碳排放效率時(shí)空分異研究

馬玉超，張永林

（武漢紡織大學(xué) 管理學(xué)院，湖北 武漢 430200）

基于2010年到2019年我國30個(gè)?。ㄊ?，自治區(qū)）紡織業(yè)的面板數(shù)據(jù)，采用多種DEA模型并結(jié)合Gini準(zhǔn)則對(duì)其碳排放效率進(jìn)行測度，并利用全局趨勢分析法、探索性空間數(shù)據(jù)分析法和標(biāo)準(zhǔn)差橢圓技術(shù)揭示其時(shí)空分異特征。研究發(fā)現(xiàn)：觀察期內(nèi)，我國各省份紡織業(yè)碳排放效率差距明顯，存在集聚趨勢；紡織業(yè)碳排放效率在東西和南北方向均呈現(xiàn)倒“U”型趨勢，且東西（南北）方向上的差距有所縮?。〝U(kuò)大）；紡織業(yè)碳排放效率自2017年開始表現(xiàn)出顯著的正向空間自相關(guān)性；紡織業(yè)碳排放效率區(qū)位分布由偏東北—偏西南轉(zhuǎn)變?yōu)槠珫|南—偏西北，分布重心向南移動(dòng)。最后，根據(jù)研究發(fā)現(xiàn)提出若干政策建議。

紡織業(yè)；Gini準(zhǔn)則；碳排放效率；時(shí)空分異

我國在“十三五”期間已經(jīng)基本實(shí)現(xiàn)了從紡織大國到紡織強(qiáng)國的轉(zhuǎn)變，我國紡織業(yè)在全球價(jià)值鏈中的位置穩(wěn)步提升，整體競爭力進(jìn)一步增強(qiáng)。有關(guān)資料顯示，我國在2020年的紡織纖維加工總量占世界纖維加工總量的一半以上，化纖產(chǎn)量更是達(dá)到了70%以上的占比，我國紡織品服裝出口額穩(wěn)居世界首位[1]。隨著全球氣候治理形勢的日益緊迫，我國提出了“爭取二氧化碳排放量在2030年前達(dá)到峰值，在2060年前實(shí)現(xiàn)碳中和”的目標(biāo)愿景，這對(duì)高能耗、高污染的紡織業(yè)形成了剛性要求。因此，在“十四五”期間，紡織業(yè)作為國民經(jīng)濟(jì)和社會(huì)發(fā)展的支柱產(chǎn)業(yè)，迫切需要實(shí)現(xiàn)高質(zhì)量發(fā)展。紡織業(yè)碳排放效率的測算是實(shí)現(xiàn)我國紡織業(yè)高質(zhì)量發(fā)展的重要決策依據(jù)，中國地大物博，由于各地區(qū)資源稟賦不同，各地區(qū)紡織業(yè)發(fā)展水平也不盡相同，紡織業(yè)碳排放效率也可能存在時(shí)空分異現(xiàn)象。于是對(duì)我國紡織業(yè)碳排放效率進(jìn)行科學(xué)測算，探索其時(shí)空分異，對(duì)于促進(jìn)我國紡織業(yè)高質(zhì)量發(fā)展與我國“雙碳”目標(biāo)的實(shí)現(xiàn)有著重大意義。

1 文獻(xiàn)回顧

關(guān)于碳排放效率，國內(nèi)外相關(guān)研究已取得不少成果。碳排放效率可以分為單要素碳排放效率和全要素碳排放效率。早期研究主要以單要素碳排放效率為主，例如以碳指數(shù)[1]、碳排放強(qiáng)度[3]、能源強(qiáng)度[4]等指標(biāo)來衡量碳排放效率，僅考慮了碳排放量和GDP、人口等單個(gè)變量的關(guān)系。而近年來學(xué)者們開始關(guān)注于全要素碳排放效率，將經(jīng)濟(jì)、能源和環(huán)境納入測算指標(biāo)，使得碳排放效率的測算指標(biāo)更為全面。如今關(guān)于碳排放效率測度的對(duì)象，已有學(xué)者將國家或者地區(qū)作為對(duì)象測算其碳排放效率，同時(shí)也有眾多學(xué)者針對(duì)不同產(chǎn)業(yè)開展了豐富的研究，例如工業(yè)[5]、交通運(yùn)輸業(yè)[6, 7]、制造業(yè)[8]、物流業(yè)[9]、農(nóng)業(yè)[10]等。然而，目前學(xué)界對(duì)紡織業(yè)的關(guān)注度相對(duì)較低，呂明元等測算了我國紡織業(yè)的生態(tài)效率并對(duì)其結(jié)構(gòu)性因素進(jìn)行分析[11]。Zhao和Lin測算了我國紡織業(yè)在碳排放約束下的能源效率[12]。Lin和Bai測算了我國各省的紡織業(yè)能源績效，發(fā)現(xiàn)在不同的時(shí)空尺度上，中國三大地理區(qū)域的動(dòng)態(tài)能源效率呈上升趨勢[13]。Gai等認(rèn)為紡織業(yè)子行業(yè)對(duì)環(huán)境造成的影響存在差異，于是將重點(diǎn)放在棉化纖紡織行業(yè)、非棉紡織行業(yè)和印染行業(yè)，對(duì)三個(gè)子行業(yè)的生態(tài)效率進(jìn)行了測度，結(jié)果發(fā)現(xiàn)三個(gè)子行業(yè)之間的差異較大[14]。從研究尺度來看，已有關(guān)于碳排放效率測算的研究主要基于國家、區(qū)域、省域和縣域?qū)用?。從測算方法來看，如今的主流方法主要有三類。第一類是通過構(gòu)建綠色技術(shù)創(chuàng)新評(píng)價(jià)指標(biāo)體系，然后采用熵值法等客觀評(píng)價(jià)方法測算碳排放效率[15, 16]，這類方法具有簡單易懂且能夠避免人為因素導(dǎo)致結(jié)果出現(xiàn)偏差等優(yōu)點(diǎn)，但是忽略了指標(biāo)本身的重要程度。第二類方法是采用隨機(jī)前沿分析法（SFA）測算碳排放效率[17, 18]，SFA是參數(shù)方法的典型代表，其優(yōu)點(diǎn)是能夠考慮隨機(jī)因素對(duì)產(chǎn)出的影響，其缺點(diǎn)是需要事先對(duì)生產(chǎn)前沿的具體形式進(jìn)行確定，且模型中只能有一個(gè)產(chǎn)出指標(biāo)。第三類方法則是數(shù)據(jù)包絡(luò)分析法（DEA），其為非參數(shù)方法的代表，這類方法可以有效解決SFA模型的局限之處，因此得到廣泛應(yīng)用，且發(fā)展出了各種形式，例如BCC[19]、超效率DEA[20]、SBM[21]、EBM[22]、MinDS[23]、DDF等模型[24]。

空間分異最初來源于地理學(xué)，是指某一屬性在不同地理區(qū)域之間存在的差異。本文增加時(shí)間維度的考量，將紡織業(yè)碳排放效率時(shí)空分異界定為“紡織業(yè)碳排放效率在地理區(qū)域之間存在的差異，以及差異在時(shí)間維度上的變化規(guī)律”。關(guān)于碳排放效率的空間分異，已有學(xué)者從不同角度開展了研究，但是目前以紡織業(yè)為對(duì)象的研究卻鳳毛麟角。王興民等基于198個(gè)地級(jí)市數(shù)據(jù)，發(fā)現(xiàn)中國城市層面的CO2排放存在明顯的空間分異[25]。牛亞文等基于縣域數(shù)據(jù)測算了長株潭地區(qū)的土地利用碳排放并分析其空間分異現(xiàn)象[26]。高新才和韓雪研究了黃河流域碳排放的空間分異，發(fā)現(xiàn)黃河流域碳排放效率呈現(xiàn)出較強(qiáng)的空間相關(guān)性且有加強(qiáng)趨勢，碳排放分布存在較強(qiáng)的轉(zhuǎn)移惰性[27]。韓傳峰等則采用多種空間數(shù)據(jù)分析方法對(duì)長三角地區(qū)城市層面的碳排放空間分異進(jìn)行分析[28]。

綜上所述，可以看出現(xiàn)有關(guān)于碳排放效率測算及其空間分異的研究已經(jīng)取得了豐碩的研究成果，但仍然存在一些不足之處：首先，現(xiàn)有關(guān)于碳排放效率測算的研究對(duì)紡織業(yè)的關(guān)注較少，而紡織業(yè)是我國經(jīng)濟(jì)與社會(huì)發(fā)展的支柱產(chǎn)業(yè)，同時(shí)也是我國開展國際合作和融合發(fā)展的優(yōu)勢產(chǎn)業(yè)，其對(duì)于民生和美化生活具有重要作用；其次，在碳排放效率測算方法上，現(xiàn)有絕大多數(shù)研究多采用單一DEA模型，由于每個(gè)測算模型都具有從自身角度出發(fā)的有價(jià)值的信息，采用單一模型可能會(huì)導(dǎo)致結(jié)果出現(xiàn)偏誤或者不可比的情況；最后，對(duì)于紡織業(yè)碳排放效率的時(shí)空分異特征，目前還幾乎沒有學(xué)者進(jìn)行探討。鑒于此，本文為了豐富紡織業(yè)碳排放效率的相關(guān)研究，從以下方面對(duì)現(xiàn)有研究進(jìn)行拓展：一是參考相關(guān)研究，采用多DEA-Gini方法對(duì)紡織業(yè)碳排放效率進(jìn)行測算，避免使用單一模型存在的風(fēng)險(xiǎn)；二是基于ArcGIS平臺(tái)，采用多種空間分析方法實(shí)證考察我國紡織業(yè)碳排放效率的時(shí)空分異特征。從碳排放效率視角對(duì)紡織業(yè)進(jìn)行研究，可以為政府部門制定有針對(duì)性的紡織業(yè)低碳發(fā)展戰(zhàn)略提供理論和實(shí)證依據(jù)，進(jìn)而有利于紡織業(yè)的轉(zhuǎn)型升級(jí)和高質(zhì)量發(fā)展。

2 研究設(shè)計(jì)

2.1 多DEA-Gini準(zhǔn)則模型

本文采用了多種DEA模型對(duì)2010年到2019年中國30個(gè)省、直轄市和自治區(qū)（不含西藏和港澳臺(tái)）的紡織業(yè)碳排放效率進(jìn)行了測度，然后采用數(shù)據(jù)挖掘分類回歸決策樹理論中的Gini準(zhǔn)則來對(duì)各DEA模型賦予權(quán)重，最終計(jì)算出結(jié)合多種DEA模型所得結(jié)果的綜合值[29]，以避免重要信息被忽略和效率結(jié)果不可比的情況。本文首先基于規(guī)模報(bào)酬不變和可變假設(shè)，然后考慮到DEA結(jié)果可能會(huì)出現(xiàn)很多“1”從而導(dǎo)致無法對(duì)決策單元進(jìn)行進(jìn)一步的比較，于是最終選擇了規(guī)模報(bào)酬不變假設(shè)下的SBM模型和超效率SBM模型、規(guī)模報(bào)酬可變假設(shè)下的SBM模型和超效率SBM模型。根據(jù)這四種不同的DEA模型測算出的不同結(jié)果，通過Gini準(zhǔn)則，得到綜合效率值。

2.1.1 標(biāo)準(zhǔn)效率和超效率SBM模型

本文采用的標(biāo)準(zhǔn)效率和超效率SBM模型的具體形式如下所示。

2.1.2 Gini準(zhǔn)則

2.1.3 多DEA-Gini準(zhǔn)則方法步驟

參考薛暉等[29]的研究，本文采用的多DEA- Gini準(zhǔn)則方法步驟如下。

2.2 全局趨勢分析法

2.3 探索性空間數(shù)據(jù)分析法

探索性空間數(shù)據(jù)分析（ESDA）是一種具有識(shí)別功能的空間數(shù)據(jù)分析方法，主要用于探測某些變量的空間相關(guān)性和集聚性。目前，該方法已經(jīng)在多個(gè)領(lǐng)域得到了廣泛應(yīng)用，最為常見的有全局莫蘭指數(shù)和局部莫蘭指數(shù)（LISA），在本文中可以檢驗(yàn)紡織業(yè)碳排放效率在整體上的空間自相關(guān)性和局部上的集聚性。全局莫蘭指數(shù)大于0則表示空間正相關(guān)，值越接近1表示具有相似屬性的地區(qū)越傾向于聚集在一起（即高值與高值鄰接，低值與低值鄰接）；全局莫蘭指數(shù)小于0表示空間負(fù)相關(guān)，值越接近-1表示具有相異屬性的地區(qū)越傾向于聚集在一起（即高值與低值鄰接，低值與高值鄰接）；全局莫蘭指數(shù)接近0時(shí)，則表示該屬性是隨機(jī)分布的，不存在空間相關(guān)性。與全局莫蘭指數(shù)相似，局部莫蘭指數(shù)（LISA）為正，則表示一個(gè)高值被高值所包圍（高—高），或者是一個(gè)低值被低值所包圍（低—低）；若指數(shù)為負(fù)，則表示一個(gè)低值被高值所包圍（低—高），或者是一個(gè)高值被低值所包圍（高—低）[31]。兩種指數(shù)的計(jì)算方法見式5和6。

2.4 標(biāo)準(zhǔn)差橢圓

標(biāo)準(zhǔn)差橢圓是一種衡量地理要素空間分布特征的統(tǒng)計(jì)技術(shù)，但該方法也在諸多其他領(lǐng)域得到廣泛應(yīng)用。在本文中，采用標(biāo)準(zhǔn)差橢圓來考察我國紡織業(yè)碳排放效率的空間集聚和區(qū)位特征。標(biāo)準(zhǔn)差橢圓分析中的重要參數(shù)包括分布重心、橢圓長短半軸和方位角等，具體計(jì)算公式參考劉華軍等[32]。

2.5 指標(biāo)設(shè)定與數(shù)據(jù)來源

根據(jù)柯布-道格拉斯生產(chǎn)函數(shù)的思想，并結(jié)合本文的特點(diǎn)，在文中主要有資本、勞動(dòng)力和能源三種投入變量。其中資本投入由各省份紡織業(yè)固定資產(chǎn)投資額表示，勞動(dòng)力投入由各省份紡織業(yè)就業(yè)人數(shù)表示，由于能源消耗量和碳排放量具有線性關(guān)系，如今國家提倡節(jié)能減排，碳排放量越小越好，符合作為DEA投入指標(biāo)的條件，因此在本文中用各省份紡織業(yè)二氧化碳排放量代替能源投入。產(chǎn)出變量采用各地區(qū)紡織業(yè)主營業(yè)務(wù)收入代表。

資本投入數(shù)據(jù)來源于2011年到2020年的《中國固定資產(chǎn)統(tǒng)計(jì)年鑒》和《中國紡織工業(yè)發(fā)展報(bào)告》，勞動(dòng)力投入數(shù)據(jù)來源于2011年到2020年的《中國勞動(dòng)統(tǒng)計(jì)年鑒》，碳排放量數(shù)據(jù)來源于中國碳核算數(shù)據(jù)庫（CEADs）（https://www.ceads.net.cn/），與IPCC公布的碳排放數(shù)據(jù)相比，中國碳核算數(shù)據(jù)庫基于表觀能耗數(shù)據(jù)和更新之后的中國的碳排放因子，從而提供的數(shù)據(jù)更加準(zhǔn)確[33]。產(chǎn)出數(shù)據(jù)來源于2011年到2020年的《中國工業(yè)統(tǒng)計(jì)年鑒》。

3 中國紡織業(yè)碳排放效率評(píng)價(jià)結(jié)果

多DEA-Gini準(zhǔn)則方法計(jì)算出的各個(gè)模型的權(quán)重值如表1所示。在四種DEA模型中，基于規(guī)模報(bào)酬可變假設(shè)的超效率SBM模型的信息純度最大，為0.0049，因此應(yīng)該給予該模型較大的權(quán)重值；基于規(guī)模報(bào)酬可變假設(shè)的標(biāo)準(zhǔn)效率SBM模型的信息純度最小，為0.0037，因此應(yīng)該給予該模型較小的權(quán)重值。最終，基于規(guī)模報(bào)酬不變和可變假設(shè)的標(biāo)準(zhǔn)效率SBM模型、基于規(guī)模報(bào)酬不變和可變假設(shè)的超效率SBM模型的權(quán)重值分別為0.228、0.223、0.254、0.295。將四種DEA模型的結(jié)果根據(jù)權(quán)重值加總所得綜合效率結(jié)果如表2所示。受限于篇幅，本文只選取2010、2013、2016、2019年的效率結(jié)果進(jìn)行匯報(bào)。

多DEA-Gini方法測算結(jié)果顯示，在觀察期內(nèi)我國紡織業(yè)碳排放效率均值有所上升，但上升幅度較小，我國絕大多數(shù)省份的紡織業(yè)碳排放效率小于1，具有較大的提升空間。為了將效率測度結(jié)果進(jìn)行直觀展示，本文借助ArcGIS10.5軟件對(duì)我國30個(gè)省份的紡織業(yè)碳排放效率進(jìn)行可視化展示。具體方法為，將我國30個(gè)省份的紡織業(yè)碳排放效率數(shù)據(jù)導(dǎo)入ArcGIS10.5軟件，采用自然斷點(diǎn)法（Natural Breaks）將各省份按照紡織業(yè)碳排放效率高低分為5個(gè)等級(jí)，分別為高、較高、中等、較低和低效率，由此得到我國紡織業(yè)碳排放效率的空間分布地圖（圖1）。從圖1中可以發(fā)現(xiàn)，高效率水平地區(qū)有青海、寧夏、海南，效率均值較高的地區(qū)有北京、內(nèi)蒙古、江蘇、安徽、河南、湖北、江西、貴州等地；低效率水平地區(qū)有黑龍江、吉林、山西、廣西和云南，效率均值較低的地區(qū)有遼寧、天津、廣東、重慶、四川、甘肅、新疆等地；其余省份效率水平居中。從空間分布格局的演變來看，隨著時(shí)間的推移，紡織業(yè)碳排放效率呈現(xiàn)出一定的聚集趨勢，具體表現(xiàn)為效率較高地區(qū)發(fā)生聚集（主要分布在東南沿海和中部地區(qū)），效率較低的地區(qū)發(fā)生聚集（主要分布在東北地區(qū)，部分北方地區(qū)和西南地區(qū)）。

表1 各模型權(quán)重值

表2 紡織業(yè)碳排放效率結(jié)果

續(xù)表2

圖1 中國紡織業(yè)碳排放效率空間分布

注：該圖基于國家自然資源部標(biāo)準(zhǔn)地圖（審圖號(hào)：GS(2019)1822號(hào)）繪制，底圖無修改。

4 中國紡織業(yè)碳排放效率時(shí)空分異分析

4.1 全局趨勢分析

以正東和正北方向作為X軸和Y軸，以紡織業(yè)碳排放效率作為 Z軸，以此直觀顯示2010到2019年我國紡織業(yè)碳排放效率的空間分布趨勢（圖2）。從趨勢線圖2中可以發(fā)現(xiàn)，我國紡織業(yè)碳排放效率在東西方向和南北方向上均呈現(xiàn)出倒“U”型趨勢，這表明我國紡織業(yè)碳排放效率存在明顯的空間非均衡狀態(tài)，總體呈現(xiàn)出中心高，周圍低的空間分布趨勢。隨著時(shí)間的推移，效率空間分布趨勢也發(fā)生了變化，圖中趨勢曲線的弧度越大，表明紡織業(yè)碳排放效率間的差異越大。從趨勢曲線的弧度變化可以發(fā)現(xiàn)，在觀察期內(nèi)，東西方向的趨勢曲線弧度有所減小，南北方向的趨勢曲線弧度有所增大，兩種方向上仍然保持著倒“U”型趨勢，這表明我國紡織業(yè)碳排放效率在東西方向上的差距有減少的趨勢，南北方向上的差距有擴(kuò)大的趨勢，且我國紡織業(yè)碳排放效率仍然顯現(xiàn)中心高，四周低的空間分布趨勢。

圖2 中國紡織業(yè)碳排放效率空間趨勢線

4.2 空間自相關(guān)性分析

全局空間自相關(guān)性檢驗(yàn)結(jié)果見表3?？梢园l(fā)現(xiàn)，我國紡織業(yè)碳排放效率的全局莫蘭指數(shù)在2010年到2016年期間為負(fù)，且未通過顯著性水平檢驗(yàn)，表明我國紡織業(yè)碳排放效率在此期間不存在顯著的空間自相關(guān)性。然而，自2017年起，莫蘭指數(shù)由負(fù)轉(zhuǎn)正，且通過了顯著性水平檢驗(yàn)，表明我國紡織業(yè)碳排放效率開始表現(xiàn)出顯著的正向空間自相關(guān)性，即高效率地區(qū)與高效率地區(qū)聚集，低效率地區(qū)與低效率地區(qū)聚集，這驗(yàn)證了前文通過直觀觀察所得出的結(jié)論。

表3 全局莫蘭檢驗(yàn)結(jié)果

全局莫蘭指數(shù)檢驗(yàn)無法反映局部地區(qū)的集聚情況，本文進(jìn)一步對(duì)局部空間自相關(guān)性進(jìn)行檢驗(yàn)，結(jié)果采用LISA集聚地圖形式展示，經(jīng)過分別對(duì)2010年、2013年、2016年和2019年的紡織業(yè)碳排放效率局部空間自相關(guān)檢驗(yàn)發(fā)現(xiàn)，2010年和2013年不存在顯著的局部空間自相關(guān)性，2016年和2019年存在顯著的局部空間自相關(guān)性。因此，限于篇幅，僅展示2016和2019年的LISA集聚地圖，如圖3所示。在2016年，只有安徽省位于高—高值集聚區(qū)，同時(shí)在5%的水平上顯著；處于高—低集聚區(qū)的只有貴州省，且通過了5%的顯著性水平檢驗(yàn)。在2019年，江西省成為唯一位于高—高值集聚區(qū)的省份，在1%的水平上顯著；內(nèi)蒙古、黑龍江、吉林、遼寧等地屬于低—低值集聚區(qū)，分別在5%、5%、1%、5%的水平上顯著。

圖3 中國紡織業(yè)碳排放效率LISA集聚圖

注：該圖基于國家自然資源部標(biāo)準(zhǔn)地圖（審圖號(hào)：GS(2019)1822號(hào)）繪制，底圖無修改。

4.3 標(biāo)準(zhǔn)差橢圓分析

我國紡織業(yè)碳排放效率的標(biāo)準(zhǔn)差橢圓及重心轉(zhuǎn)移情況如圖4所示，相關(guān)參數(shù)見表4。從橢圓形狀來看，2010年、2013年、2016年和2019年的平均形狀指數(shù)分別為0.88、0.84、0.91、0.95，呈現(xiàn)出增加態(tài)勢，形狀接近于正圓，這說明向心力較弱，意味著集聚區(qū)內(nèi)各省份的紡織業(yè)碳排放效率更加均衡。從方位角來看，紡織業(yè)碳排放效率的標(biāo)準(zhǔn)差橢圓在觀察期內(nèi)從偏東北—偏西南向偏東南—偏西北方向旋轉(zhuǎn)了112.314度，表明紡織業(yè)碳排放效率的空間分布方向具有較大的偏移趨勢。具體從長短半軸來看，長半軸在觀察期內(nèi)有所縮短，短半軸在觀察期內(nèi)有所延長，結(jié)合方位角可知起初我國紡織業(yè)碳排放效率在偏東北—偏西南方向上的空間集聚特征有所增加，在偏東南—偏西北方向上的空間集聚特征有所下降，隨后在偏東南—偏西北方向上的空間集聚特征有所增加，在偏東北—偏西南方向的空間集聚特征有所下降。從標(biāo)準(zhǔn)差橢圓的面積來看，2019年與2010年紡織業(yè)碳排放效率的標(biāo)準(zhǔn)差橢圓面積之比為1.014，變化程度較小，這說明紡織業(yè)碳排放效率的空間聚集特征未發(fā)生明顯變化。從橢圓重心位置來看，紡織業(yè)碳排放效率的分布重心均位于河南省境內(nèi)，但重心位置存在轉(zhuǎn)移趨勢，2010年位于欒川縣和嵩縣邊界處，2013年和2016年位于內(nèi)鄉(xiāng)縣境內(nèi)，2019年移動(dòng)至鄧州市境內(nèi)。整體來看，重心東西移動(dòng)范圍小于南北移動(dòng)范圍，意味著我國紡織業(yè)碳排放效率空間分布重心更傾向于在南北方向發(fā)生轉(zhuǎn)移。

圖4 中國紡織業(yè)碳排放效率標(biāo)準(zhǔn)差橢圓

注：該圖基于國家自然資源部標(biāo)準(zhǔn)地圖（審圖號(hào)：GS(2019)1822號(hào)）繪制，底圖無修改。

表4 中國紡織業(yè)碳排放效率標(biāo)準(zhǔn)差橢圓相關(guān)參數(shù)

5 結(jié)論與啟示

5.1 研究結(jié)論

本文基于省級(jí)層面的紡織業(yè)數(shù)據(jù)，采用多種DEA模型并結(jié)合Gini準(zhǔn)則對(duì)我國30個(gè)?。ㄊ?，自治區(qū)）的紡織業(yè)碳排放效率進(jìn)行了測算，并采用多種空間數(shù)據(jù)分析方法檢驗(yàn)了紡織業(yè)碳排放效率的時(shí)空分異特征。本文主要有以下發(fā)現(xiàn)：（1）多DEA-Gini方法測算結(jié)果顯示，2010年到2019年期間我國整體紡織業(yè)碳排放效率有小幅上升，絕大多數(shù)省份的效率值小于1，具有較大的提升空間。高水平地區(qū)有青海、寧夏、海南等地，低水平地區(qū)有黑龍江、吉林、山西、廣西和云南等地，且隨著時(shí)間推移，紡織業(yè)碳排放效率存在聚集趨勢。（2）全局趨勢分析結(jié)果顯示，我國紡織業(yè)碳排放效率在東西和南北方向上均呈現(xiàn)出倒“U”型趨勢，表明我國紡織業(yè)碳排放效率存在明顯的空間非均衡狀態(tài)，觀察期內(nèi)東西（南北）方向上的差距有縮?。〝U(kuò)大）的趨勢。（3）空間自相關(guān)分析發(fā)現(xiàn)，我國紡織業(yè)碳排放效率空間自相關(guān)性在2010年到2016年期間不顯著，但從2017年開始顯著，局部空間自相關(guān)檢驗(yàn)發(fā)現(xiàn)，2016年只有安徽省位于高—高值集聚區(qū)，貴州省處于高—低值集聚區(qū)，2019年江西省成為高—高值集聚區(qū)，內(nèi)蒙古和東北三省為低—低值集聚區(qū)。（4）標(biāo)準(zhǔn)差橢圓分析發(fā)現(xiàn)，在觀察期內(nèi)，集聚區(qū)內(nèi)各省份紡織業(yè)碳排放效率更加均衡，空間分布方向具有較大的偏移趨勢，紡織業(yè)碳排放效率空間分布方向由偏東北—偏西南方向發(fā)散分布轉(zhuǎn)向偏東南—偏西北方向發(fā)散分布，偏東南—偏西北方向的集聚程度有所增加，偏東北—偏西南方向的集聚程度有所下降，分布重心在觀察期內(nèi)表現(xiàn)出向南移動(dòng)的趨勢。

5.2 政策建議

本文基于研究發(fā)現(xiàn)和現(xiàn)實(shí)背景，提出以下政策建議：（1）加強(qiáng)低碳技術(shù)創(chuàng)新，推進(jìn)節(jié)能低碳發(fā)展。紡織業(yè)高質(zhì)量發(fā)展需要持續(xù)加強(qiáng)基礎(chǔ)研究投入，實(shí)現(xiàn)關(guān)鍵技術(shù)突破，同時(shí)還要優(yōu)化能源結(jié)構(gòu)，鼓勵(lì)企業(yè)使用清潔能源，提高二次能源消費(fèi)比重，以達(dá)到碳排放量減少的目的。（2）提升區(qū)域協(xié)調(diào)發(fā)展水平。我國各省份紡織業(yè)發(fā)展相對(duì)孤立，高效率水平地區(qū)輻射帶動(dòng)能力不足，因此要加強(qiáng)地區(qū)之間的聯(lián)系與合作，打破區(qū)域保護(hù)壁壘和要素流動(dòng)藩籬，高效率水平地區(qū)和低效率水平地區(qū)要做到合理分工，實(shí)現(xiàn)優(yōu)勢互補(bǔ)。（3）因地制宜地制定發(fā)展模式。我國各地區(qū)之間紡織業(yè)碳排放效率具有明顯差距，各省份應(yīng)該充分考慮自身比較優(yōu)勢，探索不同的紡織業(yè)發(fā)展路徑。比如，東部地區(qū)要借助自身良好的經(jīng)濟(jì)和技術(shù)基礎(chǔ)，大力提高產(chǎn)業(yè)低碳水平，為其他地區(qū)樹立產(chǎn)業(yè)標(biāo)準(zhǔn)；東北地區(qū)可以借助自身的工業(yè)基礎(chǔ)，優(yōu)化紡織業(yè)空間布局，挖掘地區(qū)低碳發(fā)展?jié)撃?。?）推進(jìn)高水平產(chǎn)業(yè)集聚發(fā)展。在京津冀城市群、長三角地區(qū)、粵港澳大灣區(qū)和長江經(jīng)濟(jì)帶等戰(zhàn)略區(qū)域落實(shí)行業(yè)發(fā)展重點(diǎn)，優(yōu)先建設(shè)具有世界領(lǐng)先的創(chuàng)新能力、制造能力和可持續(xù)發(fā)展能力的紡織產(chǎn)業(yè)集群，發(fā)揮龍頭企業(yè)跨區(qū)域布局的作用以帶動(dòng)中小紡織企業(yè)發(fā)展。

[1] 中國紡織工業(yè)聯(lián)合會(huì). 紡織工業(yè)“十四五”發(fā)展綱要[J]. 紡織科學(xué)研究, 2021, (07): 40-49.

[2] Mielnik O, Goldemberg J. The evolution of the “carbonization index” in developing countries[J]. Energy Policy, 1999, 27(5): 307-308.

[3] Sun J W. The decrease of CO2emission intensity is decarbonization at national and global levels[J]. Energy Policy, 2005, 33(8): 975-978.

[4] Ang B W. Is the energy intensity a less useful indicator than the carbon factor in the study of climate change?[J]. Energy Policy, 1999, 27(15): 943-946.

[5] Gao P, Yue S J, Chen H T. Carbon emission efficiency of China’s industry sectors: From the perspective of embodied carbon emissions[J]. Journal of Cleaner Production, 2021, 283: 124655.

[6] Zhao P J, Zeng L G, Li P L, et al. China’s transportation sector carbon dioxide emissions efficiency and its influencing factors based on the EBM DEA model with undesirable outputs and spatial Durbin model[J]. Energy, 2022, 238: 121934.

[7] Shao H Q, Wang Z F. Spatial network structure of transportation carbon emission efficiency in China and its influencing factors[J]. Chinese Journal of Population, Resources and Environment, 2021, 19(04): 295-303.

[8] Li J X, Cheng Z. Study on total-factor carbon emission efficiency of China’s manufacturing industry when considering technology heterogeneity[J]. Journal of Cleaner Production, 2020, 160: 121021.

[9] 范建平，肖慧，樊曉宏. 考慮非期望產(chǎn)出的改進(jìn)EBM-DEA三階段模型——基于中國省際物流業(yè)效率的實(shí)證分析[J]. 中國管理科學(xué), 2017, 25(08): 166-174.

[10] 銀西陽，賈小娟，李冬梅. 農(nóng)業(yè)產(chǎn)業(yè)集聚對(duì)農(nóng)業(yè)綠色全要素生產(chǎn)率的影響——基于空間溢出效應(yīng)視角[J/OL]. 中國農(nóng)業(yè)資源與區(qū)劃, https://kns.cnki.net/kcms/ detail/11.3513.S.20220129.1448.016.html.

[11] 呂明元，孫獻(xiàn)貞，陳維宣. 供給側(cè)結(jié)構(gòu)性因素對(duì)產(chǎn)業(yè)生態(tài)效率影響的實(shí)證研究——基于中國紡織業(yè)的DEA-Tobit分析[J]. 科技管理研究, 2018, 38(06): 240-245.

[12] Zhao H L, Lin B Q. Assessing the energy productivity of China’s textile industry under carbon emission constraints[J]. Journal of Cleaner Production, 2019, 228: 197-207.

[13] Lin B Q, Bai R. Dynamic energy performance evaluation of Chinese textile industry[J]. Energy, 2020, 199: 117388.

[14] Gai Y X, Qiao Y B, Deng H J, et al. Investigating the eco-efficiency of China’s textile industry based on a firm-level analysis[J]. Science of the Total Environment, 2022, 833: 155075.

[15] 范莉莉，江玉國. 基于熵值法的鋼鐵企業(yè)低碳競爭力評(píng)價(jià)[J]. 軟科學(xué), 2016, 30(08): 42-46.

[16] 曹陽，胡繼亮. 調(diào)整工業(yè)結(jié)構(gòu)發(fā)展三峽低碳經(jīng)濟(jì)——以湖北省宜昌市為例[J]. 長江流域資源與環(huán)境, 2012, 21(01): 17-23.

[17] 苗成林，孫麗艷，楊力. 能源消耗與碳排量約束下區(qū)域技術(shù)效率研究[J]. 科研管理, 2016, 37(02): 1-8.

[18] 朱磊，張建清，孫元元，等. 環(huán)境治理約束與中國經(jīng)濟(jì)增長——以控制碳排放為例的實(shí)證分析[J]. 中國軟科學(xué), 2018, (06): 163-171.

[19] 馬飛，胡江艷，孫啟鵬，等. 中國省際物流業(yè)碳排放績效測度及驅(qū)動(dòng)因素研究[J]. 生態(tài)經(jīng)濟(jì), 2021, 37(09): 27-33+39.

[20] 王恩旭，武春友. 基于超效率DEA模型的中國省際生態(tài)效率時(shí)空差異研究[J]. 管理學(xué)報(bào), 2011, 8(03): 443-450.

[21] 曹炳汝，鄧?yán)蚓? 長江經(jīng)濟(jì)帶物流業(yè)效率增長影響因素[J]. 經(jīng)濟(jì)地理, 2019, 39(07): 148-157.

[22] 劉華軍，郭立祥，喬列成，等. 中國物流業(yè)效率的時(shí)空格局及動(dòng)態(tài)演進(jìn)[J]. 數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究, 2021, 38(05): 57-74.

[23] 董梅，李存芳. 中國低碳技術(shù)效率、技術(shù)差距與低碳化進(jìn)程研究——基于MinDS-Luenberger方法的實(shí)證分析[J]. 華東經(jīng)濟(jì)管理, 2020, 34(11): 81-89.

[24] 邵帥，范美婷，楊莉莉. 經(jīng)濟(jì)結(jié)構(gòu)調(diào)整、綠色技術(shù)進(jìn)步與中國低碳轉(zhuǎn)型發(fā)展——基于總體技術(shù)前沿和空間溢出效應(yīng)視角的經(jīng)驗(yàn)考察[J]. 管理世界, 2022, 38(02): 46-69+4-10.

[25] 王興民，吳靜，白冰，等. 中國CO2排放的空間分異與驅(qū)動(dòng)因素——基于198個(gè)地級(jí)及以上城市數(shù)據(jù)的分析[J]. 經(jīng)濟(jì)地理, 2020, 40(11): 29-38.

[26] 牛亞文，趙先超，胡藝覺. 基于NPP-VIIRS夜間燈光的長株潭地區(qū)縣域土地利用碳排放空間分異研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2021, 41(09): 3847-3856.

[27] 高新才，韓雪. 黃河流域碳排放的空間分異及影響因素研究[J]. 經(jīng)濟(jì)經(jīng)緯, 2022, 39(01): 13-23.

[28] 韓傳峰，宋府霖，滕敏敏. 長三角地區(qū)碳排放時(shí)空特征、空間聚類與治理策略[J]. 華東經(jīng)濟(jì)管理, 2022, 36(05): 24-33.

[29] 薛暉，鄭中華，謝啟偉. 基于多種DEA模型和Gini準(zhǔn)則的效率評(píng)價(jià)方法——兼對(duì)我國高校運(yùn)營績效的評(píng)價(jià)[J]. 中國管理科學(xué), 2014, 22(4): 98-104.

[30] 李少惠，韓慧. 我國省域公共數(shù)字文化服務(wù)供給效率空間分異及驅(qū)動(dòng)因素研究[J]. 圖書情報(bào)工作, 2021, 65(7): 67-77.

[31] Anselin L. Local indicators of spatial association——LISA[J]. Geographical Analysis, 1995, 27(2): 93-115.

[32] 劉華軍，王耀輝，雷名雨. 中國戰(zhàn)略性新興產(chǎn)業(yè)的空間集聚及其演變[J]. 數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究, 2019, 36(07): 99-116.

[33] 劉華軍，郭立祥，喬列成，等. 中國物流業(yè)效率的時(shí)空格局及動(dòng)態(tài)演進(jìn)[J]. 數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究, 2021, 38(05): 57-74.

Research on the Spatial and Temporal Variation of Carbon Emission Efficiency in Textile Industry based on Multi-DEA-Gini

MA Yu-chao, ZHANG Yong-lin

(School of Management, Wuhan Textile University, Wuhan Hubei 430200, China)

Based on the panel data of textile industry in 30 provinces (municipalities, autonomous regions) in China from 2010 to 2019, the carbon emission efficiency of textile industry was measured by using multiple DEA models combined with Gini criterion, and its spatial-temporal differentiation characteristics were revealed by using global trend analysis, exploratory spatial data analysis and standard deviation ellipse technology. The results show that during the observation period, the carbon emission efficiency gap of textile industry in all provinces of China is obvious, and there is an agglomeration trend; The carbon emission efficiency of textile industry showed an inverted "U" trend in both east-west and north-south directions, and the gap narrowed (widened) in east-west (north-south) direction. The carbon emission efficiency of textile industry has shown a significant positive spatial autocorrelation since 2017. The location distribution of carbon emission efficiency of textile industry changed from northeast-southwest to southeast-northwest, and the distribution center moved southward. Finally, some policy suggestions are put forward based on the research findings.

textile industry; gini criterion; carbon emission efficiency; spatial and temporal variation

馬玉超（1970-），女，教授，碩士，研究方向：綠色供應(yīng)鏈管理.

F425

A

2095－414X(2022)06－0015－09