付云鵬，馬樹(shù)才

?

基于空間計(jì)量的碳排放量影響因素分析

付云鵬，馬樹(shù)才

(遼寧大學(xué)經(jīng)濟(jì)學(xué)院，遼寧沈陽(yáng)，110036)

在對(duì)2000?2014年中國(guó)30個(gè)地區(qū)的碳排放量進(jìn)行估算的基礎(chǔ)上，利用Moran’s I檢驗(yàn)法對(duì)區(qū)域碳排放量進(jìn)行了空間相關(guān)性檢驗(yàn)。在此基礎(chǔ)上以擴(kuò)展的STIRPAT模型為理論基礎(chǔ)，借助空間回歸模型研究了人口結(jié)構(gòu)、產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)、技術(shù)水平等因素對(duì)碳排放量的影響效應(yīng)。結(jié)果顯示：產(chǎn)業(yè)結(jié)構(gòu)、人口規(guī)模和技術(shù)水平是碳排放量的主要影響因素，產(chǎn)業(yè)結(jié)構(gòu)和人口規(guī)模對(duì)碳排放量的影響是正向的，技術(shù)水平對(duì)碳排放量的影響是負(fù)向的。

人口規(guī)模；產(chǎn)業(yè)結(jié)構(gòu)；技術(shù)水平；空間自相關(guān)；空間回歸模型

一、引言及相關(guān)文獻(xiàn)綜述

人口與資源環(huán)境的關(guān)系非常密切、復(fù)雜，人口規(guī)模、結(jié)構(gòu)、質(zhì)量等因素都與資源環(huán)境密切相關(guān)，因此，不少學(xué)者開(kāi)始關(guān)注人口與碳排放之間的關(guān)系問(wèn)題，而隨著空間計(jì)量經(jīng)濟(jì)學(xué)的發(fā)展，也使得對(duì)碳排放問(wèn)題的研究考慮了空間因素的影響。中國(guó)是一個(gè)地廣人多的國(guó)家，經(jīng)濟(jì)發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)、人口分布、消費(fèi)結(jié)構(gòu)等都存在著一定的地域差異。因此在研究碳排放量的影響因素時(shí)考慮空間因素的影響有著非常重要的理論和現(xiàn)實(shí)意義。Ehrlich和Holdren[1]為定量地研究人口對(duì)環(huán)境的影響提出了著名的IPAT方程。Dietz和Rose等[2]將IPAT方程改寫(xiě)成隨機(jī)形式并稱之為STIRPAT模型，該模型克服了IPAT方程中解釋變量同比例變動(dòng)的不足，是一個(gè)含有多個(gè)解釋變量的指數(shù)形式的模型，模型系數(shù)可以理解為環(huán)境變量對(duì)各影響因素的彈性，STIRPAT模型在應(yīng)用過(guò)程中比較靈活，同時(shí)還可以根據(jù)研究需要加入其他影響環(huán)境變量的解釋變量。該模型為研究人口因素、社會(huì)財(cái)富和技術(shù)水平等因素對(duì)環(huán)境變量的影響提供了新的依據(jù)，為研究人口與環(huán)境間的關(guān)系奠定了新的理論基礎(chǔ)。Mackellar等[3]將IPAT中的“P”用家庭戶規(guī)模(households)代替，將IPAT改寫(xiě)為IHAT模型，并用這兩種模型分析了人口因素對(duì)全國(guó)能源消費(fèi)量的影響。York等[4]利用STIRPAT模型研究了人口、經(jīng)濟(jì)生產(chǎn)、城市化和地理因素對(duì)環(huán)境的影響，并解釋了跨國(guó)排放對(duì)環(huán)境的影響。肖宏偉等[5]基于擴(kuò)展 STIRPAT模型和空間杜賓模型研究了人口、產(chǎn)業(yè)結(jié)構(gòu)、能源價(jià)格等因素對(duì)碳排放量和碳排放強(qiáng)度的影響。宋德勇和劉習(xí)平[6]基于1978－2010年中國(guó)各地區(qū)碳排放數(shù)據(jù)，對(duì)碳基尼系數(shù)和碳洛倫茲曲線進(jìn)行了重新定義，并對(duì)中國(guó)各地區(qū)2020年的碳排放量進(jìn)行了空間分配。Wang Can等[7]利用中國(guó)1957－2000年的碳排放數(shù)據(jù)，采用LMDI方法研究了碳排放量的影響因素，結(jié)果顯示能源強(qiáng)度、燃料轉(zhuǎn)換和可再生能源的利用是減少碳排放量的積極因素。Wu等[8]利用LMDI分解方法將碳排放強(qiáng)度的影響因素分解為產(chǎn)業(yè)結(jié)構(gòu)、能源強(qiáng)度、能源結(jié)構(gòu)和發(fā)射系數(shù)四種。Ren等[9]利用Divisia指數(shù)法探討了制造業(yè)碳排放量的影響因素。范丹[10]利用擴(kuò)展Johan 恒等式對(duì)中國(guó)能源消費(fèi)人均碳排放的影響因素進(jìn)行了分析，利用LMDI分解法對(duì)產(chǎn)業(yè)能源消耗碳排放進(jìn)行分解，根據(jù)公式測(cè)算出1995－2010年我國(guó)能源結(jié)構(gòu)、潛在能源強(qiáng)度、能源績(jī)效、能源技術(shù)進(jìn)步等因素對(duì)碳排放的影響效應(yīng)。許海平[11]采用2000－2008年我國(guó)29個(gè)省域的空間樣本數(shù)據(jù)研究了城市化水平、貿(mào)易開(kāi)放、非農(nóng)化就業(yè)人員的比重和技術(shù)進(jìn)步對(duì)碳排放量的空間影響。Stern D. I. 和Common M. S.[12]應(yīng)用73個(gè)國(guó)家1960－1990年的樣本數(shù)據(jù)借助環(huán)境EKC研究了人均SO2排放量與人均財(cái)富的關(guān)系，研究結(jié)果表明：減排與時(shí)間有關(guān)系而不是與收入有關(guān)系。Marzio Galeottia和Alessandro Lanzab[13]利用1960?1995年108個(gè)國(guó)家的CO2作為環(huán)境污染的代理變量，研究了人均CO2與人均GDP之間的關(guān)系，證明了環(huán)境EKC的“倒U”型關(guān)系存在。Amy K. Richmond和Robert K. Kaufmann[14]研究了收入和能源使用、收入和碳排放之間的關(guān)系，結(jié)果顯示：對(duì)于OECD國(guó)家，收入與人均能源使用或者碳排放之間的拐點(diǎn)存在；對(duì)于非OECD國(guó)家，收入與人均能源使用或者碳排放之間的拐點(diǎn)不存在。Song T等[15]利用中國(guó)1985－2005年的省級(jí)面板數(shù)據(jù)，基于環(huán)境庫(kù)茲涅茨曲線假說(shuō)，以廢氣、廢水、固體廢物作為環(huán)境指標(biāo)，以GDP作為經(jīng)濟(jì)指標(biāo)，對(duì)廢氣、廢水、固體廢物與GDP之間的關(guān)系進(jìn)行了實(shí)證分析。數(shù)據(jù)的協(xié)整檢驗(yàn)顯示：廢氣、廢水、固體廢物的人均排放與人均GDP之間均存在著長(zhǎng)期的協(xié)整關(guān)系，三種污染物與GDP之間均呈“倒U”型關(guān)系。王建軍和周曉唯[16]以空間計(jì)量模型為基礎(chǔ)研究了人口、能源消費(fèi)和碳排放量對(duì)我國(guó)經(jīng)濟(jì)的空間溢出效應(yīng)和空間關(guān)聯(lián)性。杜慧濱等[17]利用Moran 指數(shù)和空間面板數(shù)據(jù)模型研究了1997－2009年我國(guó)29個(gè)地區(qū)碳排放績(jī)效的空間分布特征與影響因素，結(jié)果表明：我國(guó)區(qū)域碳排放績(jī)效存在著較強(qiáng)的區(qū)域集群現(xiàn)象。陳德湖和張津[18]利用2000－2009年中國(guó)30個(gè)省的數(shù)據(jù)，采用空間面板計(jì)量經(jīng)濟(jì)模型研究了經(jīng)濟(jì)增長(zhǎng)與碳排放之間的關(guān)系，結(jié)果表明：經(jīng)濟(jì)發(fā)展與碳排放之間存在著顯著的空間效應(yīng)。本文將在前人研究成果的基礎(chǔ)上，利用Moran’s I檢驗(yàn)法檢驗(yàn)中國(guó)區(qū)域碳排放量的空間相關(guān)性，在空間相關(guān)性檢驗(yàn)的基礎(chǔ)上將空間因素考慮到碳排放量影響因素模型的構(gòu)建中，研究人口規(guī)模、人口年齡結(jié)構(gòu)、城市化率、經(jīng)濟(jì)產(chǎn)出、產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)水平和能源消費(fèi)結(jié)構(gòu)等因素對(duì)中國(guó)區(qū)域碳排放量的影響效應(yīng)。

二、碳排放量的STIRPAT模型的構(gòu)建

(一)變量選取及數(shù)據(jù)來(lái)源

目前大多數(shù)國(guó)內(nèi)外學(xué)者對(duì)碳排放問(wèn)題的研究已由理論分析轉(zhuǎn)向?qū)嵶C分析，主要是利用IPAT方程[19?20]及其擴(kuò)展模型STIRPAT[21?23]、LMDI[24?25]分解法，分析碳排放量的影響因素及影響效應(yīng)。碳排放量的影響因素包括經(jīng)濟(jì)規(guī)模、人口規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)和能源效率等。在借鑒相關(guān)研究成果的基礎(chǔ)上，選取2000－2014年中國(guó)30個(gè)地區(qū)的省際樣本數(shù)據(jù)，由于西藏地區(qū)大量能源消費(fèi)數(shù)據(jù)缺失，因此不包含西藏地區(qū)。選取人均 GDP、人口規(guī)模(年末人口總數(shù))、人口年齡結(jié)構(gòu)(勞動(dòng)年齡人口占總?cè)丝诘谋戎?、產(chǎn)業(yè)結(jié)構(gòu)(第二產(chǎn)業(yè)增加值占地區(qū)生產(chǎn)總值的比重)、城市化率(城鎮(zhèn)人口占總?cè)丝诘谋戎?、技術(shù)水平(專利申請(qǐng)授權(quán)數(shù))和能源消費(fèi)結(jié)構(gòu)(煤炭消費(fèi)占能源消費(fèi)總量的比重)7個(gè)指標(biāo)作為碳排放量的影響因素來(lái)研究各因素對(duì)碳排放量的影響效應(yīng)。研究過(guò)程中考慮空間因素對(duì)碳排放量的影響，通過(guò)構(gòu)建空間回歸模型進(jìn)行實(shí)證分析。人均GDP、產(chǎn)業(yè)結(jié)構(gòu)、人口年齡結(jié)構(gòu)、城鎮(zhèn)化率、技術(shù)水平數(shù)據(jù)來(lái)源于《中國(guó)統(tǒng)計(jì)年鑒》歷年數(shù)據(jù)，能源消費(fèi)結(jié)構(gòu)數(shù)據(jù)由《中國(guó)能源統(tǒng)計(jì)年鑒》歷年數(shù)據(jù)經(jīng)簡(jiǎn)單計(jì)算獲得，由于碳排放量沒(méi)有統(tǒng)計(jì)數(shù)據(jù)，在借鑒李丹丹等[26]、程葉青等[27]研究成果的基礎(chǔ)上，利用焦炭、煤炭、煤油、柴油、汽油、原油、燃料油和天然氣8類能源的消費(fèi)量來(lái)對(duì)碳排放量進(jìn)行估算，公式如下：

式中：2表示碳排放量，單位為萬(wàn)噸；C表示各類能源的碳排放系數(shù)；E表示第種能源的標(biāo)準(zhǔn)煤消費(fèi)量。《中國(guó)能源統(tǒng)計(jì)年鑒》中提供的是各類能源消費(fèi)的實(shí)物消費(fèi)量，在利用公式(1)進(jìn)行碳排放量的估算之前需將實(shí)物消費(fèi)量都轉(zhuǎn)換為標(biāo)準(zhǔn)煤消費(fèi)量，各類能源的碳排放系數(shù)和標(biāo)準(zhǔn)煤折算系數(shù)數(shù)據(jù)具體見(jiàn)表1。

表1 能源的碳排放系數(shù)和標(biāo)準(zhǔn)煤折算系數(shù)表

注：標(biāo)準(zhǔn)煤折算系數(shù)來(lái)源于《中國(guó)能源統(tǒng)計(jì)年鑒2013》，碳排放系數(shù)來(lái)源于IPCC碳排放計(jì)算指南；*表示天然氣的標(biāo)準(zhǔn)煤折算系數(shù)單位為千克標(biāo)準(zhǔn)煤/立方米

(二)模型構(gòu)建

IPAT 方程是用來(lái)研究人類活動(dòng)對(duì)環(huán)境壓力的影響的模型，目前該模型已被廣泛地應(yīng)用于研究環(huán)境問(wèn)題，Dietz 和Rosa[2]將 IPAT方程表示成隨機(jī)形式，即 STIRPAT模型：

式中：表示環(huán)境變量；表示人口規(guī)模；表示人均財(cái)富；表示技術(shù)水平；表示隨機(jī)誤差；表示模型系數(shù)；，，表示各因素的參數(shù)。將模型(2)兩邊同時(shí)取對(duì)數(shù)，得到如下線性模型：

ln=ln+ln+ln+ln+ln(3)

在人口對(duì)環(huán)境的影響中，人口規(guī)模是最主要的影響因素。首先，人口作為生產(chǎn)者通過(guò)生產(chǎn)活動(dòng)對(duì)環(huán)境產(chǎn)生影響，人口規(guī)模越大，在就業(yè)率不變的情況下，就業(yè)人口越多，從而將會(huì)有更多的生產(chǎn)者從事生產(chǎn)活動(dòng)，生產(chǎn)規(guī)模的擴(kuò)大，勢(shì)必帶動(dòng)更多的能源消費(fèi)、更多的碳排放和更多的其他環(huán)境污染物的排放；其次，人口同時(shí)也是消費(fèi)者，人口消費(fèi)的各種物質(zhì)生活資料在生產(chǎn)過(guò)程中產(chǎn)生間接的碳排放量。用年末人口數(shù)()代表人口規(guī)模。

除人口規(guī)模外，與碳排放量密切相關(guān)的人口因素還包括人口結(jié)構(gòu)，其中最主要的是人口年齡結(jié)構(gòu)和人口城鄉(xiāng)結(jié)構(gòu)。近年來(lái)中國(guó)的老齡化現(xiàn)象凸顯，人口年齡結(jié)構(gòu)通過(guò)就業(yè)、儲(chǔ)蓄率、消費(fèi)結(jié)構(gòu)等對(duì)碳排放量產(chǎn)生影響，勞動(dòng)年齡人口比重越大，從事生產(chǎn)活動(dòng)的人口越多，從而帶來(lái)更多的碳排放量，同時(shí)，年輕人的消費(fèi)理念和消費(fèi)方式與老年人不同，年輕人更傾向于汽車等高排放的生活消費(fèi)品的消費(fèi)，因此，人口年齡結(jié)構(gòu)對(duì)碳排放量也存在影響。用勞動(dòng)年齡人口占總?cè)丝诘谋戎卮砣丝谀挲g結(jié)構(gòu)；人口城鄉(xiāng)結(jié)構(gòu)()對(duì)碳排放量的影響體現(xiàn)在城鄉(xiāng)生活方式的不同；用城市化率()代表人口城鄉(xiāng)結(jié)構(gòu)。

人均財(cái)富對(duì)碳排放量的影響不言而喻，經(jīng)濟(jì)產(chǎn)出越多，人均財(cái)富越多。但產(chǎn)出多的同時(shí)，能源消費(fèi)越多，從而能源消費(fèi)帶來(lái)的碳排放量越多，用人均地區(qū)生產(chǎn)總值代表人均財(cái)富()。

第二產(chǎn)業(yè)相對(duì)于農(nóng)業(yè)和第三產(chǎn)業(yè)需要更多的能源消費(fèi)，因此第二產(chǎn)業(yè)比重越大，碳排放量越大；用第二產(chǎn)業(yè)增加值占地區(qū)生產(chǎn)總值的比重代表產(chǎn)業(yè)結(jié)構(gòu)()。

技術(shù)水平對(duì)碳排放量的影響主要通過(guò)提高能源利用效率、改變生產(chǎn)方式等途徑來(lái)實(shí)現(xiàn)，減少碳排放量的主要手段只能依賴技術(shù)進(jìn)步。目前大部分學(xué)者采用研發(fā)投入和專利數(shù)量作為技術(shù)水平的代理變量，本文采用專利申請(qǐng)授權(quán)數(shù)代表技術(shù)水平()。

能源產(chǎn)品中煤炭類能源的碳排放系數(shù)最大，因此煤炭類能源消費(fèi)比重越大，碳排放總量越大。本文用能源消費(fèi)總量中煤炭類能源的比重代表能源消費(fèi)結(jié)構(gòu)()。

根據(jù)上述分析，以STIRPAT模型為理論基礎(chǔ)，以區(qū)域碳排放量(2)作為環(huán)境變量，以人口規(guī)模()、人口年齡結(jié)構(gòu)()、城市化率()作為影響環(huán)境的人口因素，以人均GDP()表示人均財(cái)富，同時(shí)考慮產(chǎn)業(yè)結(jié)構(gòu)()、技術(shù)水平()和能源消費(fèi)結(jié)構(gòu)()對(duì)碳排放量的影響，構(gòu)建如下模型：

ln2i=0+1lnpop+2lnpp+3lncpp+4lngdp+

5lnind+6lnpa+7lnes+(4)

式中：表示中國(guó)30個(gè)不同的地區(qū)。

三、區(qū)域碳排放量的空間回歸分析

(一)空間相關(guān)性檢驗(yàn)

1. 全局空間相關(guān)性檢驗(yàn)

全局空間相關(guān)性用以分析空間數(shù)據(jù)在整個(gè)系統(tǒng)內(nèi)所表現(xiàn)出的分布特征，目前國(guó)內(nèi)外學(xué)者大多 采用Moran’s I來(lái)度量變量的全局空間相關(guān)性，全局Moran’s I的定義為：

利用空間分析軟件GeoDa對(duì)中國(guó)30個(gè)地區(qū)2000－2014年的區(qū)域碳排放量數(shù)據(jù)進(jìn)行空間自相關(guān)分析，歷年全局空間自相關(guān)的Moran’s I值見(jiàn)表2。

表2 2000－2014年中國(guó)區(qū)域碳排放量的Moran’s I值

從表2中可以看出2000－2014年中國(guó)區(qū)域碳排放量間存在明顯的空間自相關(guān)性，空間相關(guān)性指數(shù)Moran’s I值總體上呈現(xiàn)先下降，后上升，再下降的趨勢(shì)。Moran’s I值最大的是2001年，為0.2907，此后Moran’s I值呈下降趨勢(shì)，下降至2003年的0.2517，2004年又開(kāi)始上升，2006年之后呈現(xiàn)逐漸下降的趨勢(shì)。這說(shuō)明各地區(qū)碳排放量的空間集聚程度有所減弱，中國(guó)區(qū)域碳排量的總體空間差異有所擴(kuò)大，但仍然呈現(xiàn)比較明顯的空間自相關(guān)性。

2. 局部空間相關(guān)性檢驗(yàn)

全局空間自相關(guān)指數(shù) Moran’s I揭示了碳排放量的整體空間依賴程度，尚未指出哪些局部地區(qū)出現(xiàn)了高值或低值的空間聚集情況。為進(jìn)一步檢驗(yàn)不同年份各地區(qū)碳排放量是否具有高值或低值在空間上的集聚現(xiàn)象，還需要進(jìn)行局部空間自相關(guān)檢驗(yàn)。局部Moran’sI(或稱LISA)是檢驗(yàn)局部地區(qū)之間是否存在相似或相異的觀察值聚集情況的一種指數(shù)。局部Moran’s I的定義為：

若I大于0，表示特征值相似的地區(qū)是集聚的(“高-高”或“低-低”)；若I小于0，表示特征值不相似的地區(qū)是集聚的(“高-低”或“低-高”)。根據(jù)不同地區(qū)的時(shí)空分布特點(diǎn)，選取了具有代表性的年份(2001、2003和2012年) 用GeoDa軟件進(jìn)行分析，可以得到這3年中國(guó)區(qū)域碳排放量的Moran’s I散點(diǎn)圖，如圖1所示。

圖1 2001、2003和2012年中國(guó)區(qū)域碳排放量散點(diǎn)圖

圖1中，2001年的散點(diǎn)圖中遼寧省、山西省、山東省、河北省、江蘇省、浙江省、黑龍江省、內(nèi)蒙古自治區(qū)、河南省、上海市和安徽省共有11個(gè)區(qū)域位于散點(diǎn)圖的第一象限，碳排放量高的地區(qū)與碳排放量高的地區(qū)聚集在一起，呈現(xiàn)“高-高”的正向空間相關(guān)關(guān)系；吉林省、陜西省、福建省、寧夏回族自治區(qū)、北京市和天津市共有6個(gè)區(qū)域位于第二象限，碳排放量低的地區(qū)被周圍碳排放量高的地區(qū)包圍，呈現(xiàn)“低-高”聚集的空間關(guān)系；江西省、湖南省、廣西壯族自治區(qū)、海南省、重慶市、四川省、貴州省、云南省、甘肅省、青海省和新疆維吾爾自治區(qū)共有11個(gè)區(qū)域位于第三象限，碳排放量低的地區(qū)聚集在一起，呈現(xiàn)“低-低”聚集的空間關(guān)系；湖北省和廣東省共有2個(gè)區(qū)域位于第四象限，碳排放量高的地區(qū)被周圍碳排放量低的地區(qū)包圍，呈現(xiàn)“高-低”聚集的空間關(guān)系。2001年的Moran’s I散點(diǎn)圖中“高-高”和“低-低”聚集的類型居于主導(dǎo)地位。2003年的散點(diǎn)圖中，遼寧省、山西省、山東省、河北省、江蘇省、浙江省、黑龍江省、內(nèi)蒙古自治區(qū)、河南省和上海市共10個(gè)區(qū)域依舊位于第一象限；吉林省、陜西省、福建省、寧夏回族自治區(qū)、北京市和天津市共有6個(gè)區(qū)域依舊位于第二象限；江西省、湖南省、廣西壯族自治區(qū)、海南省、重慶市、貴州省、云南省、甘肅省、青海省和新疆維吾爾自治區(qū)共有10個(gè)區(qū)域位于第三象限；四川省、湖北省和廣東省共有3個(gè)區(qū)域位于第四象限。從2001年到2003年各地區(qū)的空間聚集現(xiàn)象差別不大，只有安徽省和四川省2個(gè)地區(qū)發(fā)生了變化。2012年的散點(diǎn)圖中，山西省、山東省、河北省、江蘇省、浙江省、內(nèi)蒙古自治區(qū)、河南省和上海市共8個(gè)區(qū)域位于第一象限；吉林省、上海市、北京市、天津市、安徽省、黑龍江省和寧夏回族自治區(qū)共7個(gè)區(qū)域位于第二象限；福建省、江西省、湖南省、廣西壯族自治區(qū)、海南省、重慶市、四川省、貴州省、云南省、甘肅省、青海省和新疆維吾爾自治區(qū)共有12個(gè)區(qū)域位于第三象限；遼寧省和廣東省位于第四象限。上述3個(gè)時(shí)點(diǎn)的Moran’s I散點(diǎn)圖總體上呈現(xiàn)出東部和西部地區(qū)表現(xiàn)出相似的空間關(guān)聯(lián)性(“高-高”或“低-低”)，而中部地區(qū)則表現(xiàn)出非相似的空間關(guān)聯(lián)性(“低-高”或“高-低”)，這說(shuō)明中國(guó)區(qū)域碳排放量的空間依賴性和差異性是并存的，但以空間依賴性為主。

(二)區(qū)域碳排放量的空間回歸分析

由于中國(guó)區(qū)域碳排放量間空間相關(guān)性的存在，使得線性回歸模型中空間均質(zhì)的假設(shè)遭到破壞，因此有必要構(gòu)建空間回歸模型對(duì)中國(guó)區(qū)域碳排放量的影響因素及影響程度進(jìn)行分析，在模型(4)的基礎(chǔ)上考慮空間因素對(duì)區(qū)域碳排放量的影響，構(gòu)建空間回歸模型。常見(jiàn)的空間回歸模型有空間誤差模型(SEM)和空間滯后模型(SLM)兩種。

1. 空間誤差模型(Spatial Error Model，SEM)

空間誤差模型用不同地區(qū)間的空間協(xié)方差反映隨機(jī)誤差過(guò)程，考慮空間因素的碳排放量的SEM為：

ln2i=0+1lnpop+2lnpp+3lncpp+4lngdp+

5lnind+6lnpa+7lnes++(7)

式中：0為待估計(jì)常數(shù)項(xiàng)；參數(shù)(=1, 2, …, 7)為各解釋變量對(duì)碳排放量的影響系數(shù)；參數(shù)為回歸殘差之間的空間相關(guān)性強(qiáng)度，用來(lái)衡量樣本觀測(cè)值中的空間依賴作用；為回歸殘差；是正態(tài)分布的隨機(jī)誤差。

2. 空間滯后模型(Spatial Lag Model, SLM)

空間滯后模型考慮的是某一空間對(duì)象的因變量不僅與該對(duì)象的自變量有關(guān)，同時(shí)也與相鄰對(duì)象的因變量有關(guān)。考慮空間因素的碳排放量的SLM為：

ln2i=ln2i+0+1lnpop+2lnpp+

3lncpp+4lngdp+5lnind+6lnpa+

7lnes+(8)

式中：為滯后變量系數(shù)，用來(lái)反映相鄰的空間對(duì)象之間的空間溢出效應(yīng)；0為待估計(jì)常數(shù)項(xiàng)；(=1, 2, …, 7)為各解釋變量對(duì)碳排放量的影響系數(shù)；w,j為空間權(quán)重，反映了空間距離對(duì)區(qū)域行為的作用；為隨機(jī)誤差。

3. 空間回歸模型的適用性檢驗(yàn)

在利用空間回歸模型進(jìn)行分析前，需要對(duì)模型的適用性進(jìn)行檢驗(yàn)?？臻g回歸模型的適用性檢驗(yàn)標(biāo)準(zhǔn)為：若空間滯后模型的拉格朗日統(tǒng)計(jì)量LM-lag比空間誤差模型的拉格朗日統(tǒng)計(jì)量LM-error顯著，則采用空間滯后模型進(jìn)行回歸分析；反之，則采用空間誤差模型進(jìn)行回歸分析。

4. 空間回歸分析

空間滯后模型和空間誤差模型都是用于分析截面數(shù)據(jù)的模型，但是單純應(yīng)用截面數(shù)據(jù)進(jìn)行分析容易忽略變量之間的滯后效應(yīng)和忽略某些數(shù)據(jù)信息，因此借鑒呂健等[28]學(xué)者的觀點(diǎn)采用2000－2014年數(shù)據(jù)均值進(jìn)行空間回歸分析。將樣本數(shù)據(jù)帶入模型(4)，用GeoDa軟件先對(duì)模型進(jìn)行普通最小二乘估計(jì)(OLS)，估計(jì)結(jié)果見(jiàn)表3。

表3 中國(guó)區(qū)域碳排放量的OLS和SEM估計(jì)結(jié)果表

變量7個(gè)變量3個(gè)變量 OLSSEMOLSSEM C?20.669**(0.037)?27.572**(0.037)?0.340(0.867)0.045(0.791) lnpop1.046***(0.000)1.091***(0.000)0.486***(0.001)0.580***(0.000) lnind1.170**(0.024)0.973**(0.020)1.641***(0.002)1.096***(0.000) lnpa?0.258(0.117)?0.362***(0.001)?0.113***(0.001)?0.134*(0.060) lnpp2.123(0.483)2.799(0.273) lncpp0.274(0.791)?0.714(0.415) lngdp0.773(0.133)1.557***(0.000) lnes0.221(0.452)0.546**(0.015) λ0.780***(0.002)0.571***(0.000) R20.9840.9900.9870.988 F294.311570.533 AIC35.48430.54838.48828.599 SC47.69542.76044.59334.705

注：系數(shù)下面的括號(hào)里為變量顯著性檢驗(yàn)的P值，***、**和*分別表示1% 、5%和10%的顯著性水平

表3中可以看出，普通最小二乘估計(jì)的模型2為0.984，模型的擬合優(yōu)度很好。模型的F統(tǒng)計(jì)量為294.311，在1%的顯著水平下通過(guò)模型的線性關(guān)系檢驗(yàn)。解釋變量的顯著性t檢驗(yàn)顯示：變量ln通過(guò)1%顯著性水平的變量顯著性檢驗(yàn)，變量ln通過(guò)5%顯著性水平的變量顯著性檢驗(yàn)，變量ln通過(guò)10%顯著性水平的顯著性檢驗(yàn)，變量ln、ln、ln和ln沒(méi)有通過(guò)檢驗(yàn)，說(shuō)明這4個(gè)變量對(duì)被解釋變量的解釋能力不顯著，可能存在多重共線性。進(jìn)一步進(jìn)行多重共線性檢驗(yàn)發(fā)現(xiàn)：ln的方差膨脹因子為14.672，ln的方差膨脹因子為14.672，ln的方差膨脹因子為11.239，都大于臨界值10，自變量間存在較高的多重共線性，不能同時(shí)作為ln2的解釋變量進(jìn)入模型。采用逐步回歸的方法來(lái)消除解釋變量間的多重共線性，得到只含有l(wèi)n、ln和ln這3個(gè)變量的模型。模型的OLS估計(jì)結(jié)果見(jiàn)表3。含有3個(gè)變量的模型的共線性檢驗(yàn)顯示：3個(gè)變量均通過(guò)了多重共線性檢驗(yàn)，并且3個(gè)變量均通過(guò)了5%的變量顯著性檢驗(yàn)，同時(shí)模型的統(tǒng)計(jì)量和2也都優(yōu)于包含7個(gè)解釋變量的模型，說(shuō)明含有3個(gè)解釋變量的模型更加適用于分析碳排放量的影響因素。

模型的適用性檢驗(yàn)結(jié)果表明：空間誤差模型的拉格朗日統(tǒng)計(jì)量LM-error比空間滯后模型的拉格朗日統(tǒng)計(jì)量LM-lag更加顯著，應(yīng)采用空間誤差模型進(jìn)行空間回歸分析，將樣本數(shù)據(jù)帶入模型(7)，用GeoDa軟件進(jìn)行分析，得到表3中的空間誤差模型的估計(jì)結(jié)果，為了便于比較同時(shí)給出了含有7個(gè)變量的空間誤差模型的估計(jì)結(jié)果，但分析過(guò)程以含有3個(gè)變量的模型為主。

從表3中空間誤差模型的回歸分析結(jié)果可以看出，人口規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)和技術(shù)水平是影響中國(guó)區(qū)域碳排放量的主要因素。人口規(guī)模和產(chǎn)業(yè)結(jié)構(gòu)對(duì)碳排放量的影響都是正的，說(shuō)明人口規(guī)模和第二產(chǎn)業(yè)比重的增加都將帶來(lái)碳排放量的增加，產(chǎn)業(yè)結(jié)構(gòu)的影響系數(shù)為1.096，大于人口規(guī)模的影響系數(shù)0.580，說(shuō)明產(chǎn)業(yè)結(jié)構(gòu)對(duì)碳排放量的影響要大于人口規(guī)模對(duì)碳排放量的影響。技術(shù)水平對(duì)碳排放量的影響是負(fù)的，說(shuō)明技術(shù)水平的提高將帶來(lái)碳排放量的減少，技術(shù)水平對(duì)碳排放量的影響系數(shù)為?0.134，說(shuō)明技術(shù)水平每提高1%，碳排放量將下降0.134%。此外，碳排放量空間誤差模型的空間相關(guān)程度系數(shù)為0.571，并且通過(guò)1%顯著性檢驗(yàn)，說(shuō)明碳排放量的空間溢出效應(yīng)為正，鄰近地區(qū)的碳排放量對(duì)本地區(qū)的碳排放量存在正向的空間溢出效應(yīng)，本地區(qū)的碳排放量除了受到本地區(qū)人口、經(jīng)濟(jì)、產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)水平等影響外，還受到鄰近地區(qū)空間溢出效應(yīng)的影響，鄰近地區(qū)的碳排放量增加1%，將會(huì)對(duì)本地區(qū)的碳排放量帶來(lái)0.571%的空間溢出效應(yīng)。這是傳統(tǒng)線性模型無(wú)法反映的。

四、結(jié)論與建議

在對(duì)2000?2014年中國(guó)區(qū)域碳排放量進(jìn)行估算的基礎(chǔ)上，采用Moran’s I檢驗(yàn)法對(duì)區(qū)域碳排放量的空間相關(guān)性進(jìn)行了檢驗(yàn)，在此基礎(chǔ)上利用空間誤差模型研究了碳排放量的空間依賴性以及人口規(guī)模、人均GDP、城市化率、人口年齡結(jié)構(gòu)、產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)水平、能源消費(fèi)結(jié)構(gòu)等因素對(duì)碳排放量的影響效應(yīng)，得到以下結(jié)論，并提出相應(yīng)的對(duì)策建議。

1. 中國(guó)30個(gè)地區(qū)2000?2014年的碳排放量間存在著明顯的空間自相關(guān)性

全局空間自相關(guān)的Moran’s I值最大的是2001年，此后空間自相關(guān)的Moran’s I值總體上呈現(xiàn)下降趨勢(shì)。東部地區(qū)碳排放量“高-高”聚集的現(xiàn)象比較明顯，碳排放量高的地區(qū)其鄰近省份的碳排放量也較高。西部地區(qū)碳排放量“低-低”聚集的現(xiàn)象比較明顯。

在碳排放量的Moran’s I散點(diǎn)圖中，對(duì)于位于第一象限呈現(xiàn)“高-高”聚集現(xiàn)象的區(qū)域，這些地區(qū)自身的碳排放量較高，同時(shí)周圍地區(qū)對(duì)其存在空間溢出效應(yīng)，因此在減排政策制定的時(shí)候應(yīng)先從這些地區(qū)入手，將這些地區(qū)視為減排的重點(diǎn)區(qū)域；對(duì)于位于第二象限呈現(xiàn)“低-高”聚集現(xiàn)象的區(qū)域，這些地區(qū)中心區(qū)域碳排放量較低，鄰近地區(qū)碳排放量較高，在減排政策制定的時(shí)候應(yīng)采取保守的減排態(tài)度，由于空間溢出效應(yīng)的存在鄰近地區(qū)的碳排放可能會(huì)對(duì)中心地區(qū)造成負(fù)面的影響，因此這些地區(qū)在制定減排政策的時(shí)候應(yīng)密切關(guān)注周圍地區(qū)的減排政策的制定和實(shí)施情況；對(duì)于位于第三象限呈現(xiàn)“低-低”聚集現(xiàn)象的區(qū)域，碳排放量較低的地區(qū)聚集在一起，可以視為減排的緩沖區(qū)，這些地區(qū)碳排放量相對(duì)較低，減排潛力弱，這些地區(qū)自身和鄰近地區(qū)碳排放量都比較低，高排放地區(qū)可以借鑒這些地區(qū)的減排經(jīng)驗(yàn)，實(shí)現(xiàn)減排的目標(biāo)；對(duì)于位于第四象限呈現(xiàn)“高-低”聚集現(xiàn)象的區(qū)域，可視為重點(diǎn)的減排區(qū)域，減排潛力大，應(yīng)制定積極的減排政策，以結(jié)構(gòu)調(diào)整促進(jìn)節(jié)能減排，對(duì)高能耗行業(yè)的新建和擴(kuò)建進(jìn)行適度的限制和置換，發(fā)展清潔能源和再生能源，同時(shí)大力推廣和發(fā)展節(jié)能環(huán)保行業(yè)。

2. 產(chǎn)業(yè)結(jié)構(gòu)是最重要的影響因素

第二產(chǎn)業(yè)比重的增加將大幅促進(jìn)碳排放量的增加，其影響系數(shù)為1.096，可見(jiàn)產(chǎn)業(yè)結(jié)構(gòu)和碳排放量之間幾乎是同比例增長(zhǎng)的。第二產(chǎn)業(yè)中以能源密集型產(chǎn)業(yè)為主，長(zhǎng)期的粗放型增長(zhǎng)帶來(lái)大量的能源浪費(fèi)。因此第二產(chǎn)業(yè)比重的增加勢(shì)必帶來(lái)能源消費(fèi)量的增加，從而產(chǎn)生更多的碳排放量。因此各地區(qū)應(yīng)積極優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，大力扶持和發(fā)展第三產(chǎn)業(yè)，降低第二產(chǎn)業(yè)在整個(gè)經(jīng)濟(jì)產(chǎn)出中的比重，以達(dá)到降低污染，減少碳排放量的目的。

除產(chǎn)業(yè)結(jié)構(gòu)外人口規(guī)模對(duì)碳排放量的影響也是正的，我國(guó)是一個(gè)人口大國(guó)，人口基數(shù)大是我國(guó)的一個(gè)基本國(guó)情，近年來(lái)各地區(qū)的人口自然增長(zhǎng)率均已得到控制，但人口自然增長(zhǎng)率降低引發(fā)的新問(wèn)題是人口老齡化現(xiàn)象嚴(yán)重，盡管人口老齡化對(duì)于碳排放量有一定程度的抑制作用，但依靠增加老年人口比重來(lái)減少碳排放量并不是理想的選擇。

3. 技術(shù)水平對(duì)碳排放量的影響是負(fù)的，技術(shù)水平的提高將帶來(lái)碳排放量的減少

技術(shù)水平對(duì)碳排放量的抑制作用主要是通過(guò)提高能源利用效率來(lái)實(shí)現(xiàn)，特別是通過(guò)新技術(shù)的應(yīng)用和新設(shè)備的研發(fā)來(lái)提高能源的利用效率，達(dá)到減少能源使用和減少碳排放量的目的。另一方面新能源的開(kāi)發(fā)和應(yīng)用也將大幅度降低能源的消費(fèi)量，開(kāi)展能源節(jié)約、資源循環(huán)利用、發(fā)展綠色科技，從而達(dá)到降低碳排放量的目的。此外，技術(shù)進(jìn)步也會(huì)帶來(lái)經(jīng)濟(jì)增長(zhǎng)方式的轉(zhuǎn)變，促進(jìn)產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化，從而間接地減少碳排放量。因此應(yīng)加速技術(shù)進(jìn)步的步伐，鼓勵(lì)新技術(shù)的應(yīng)用和新產(chǎn)品的研發(fā)，積極引進(jìn)先進(jìn)的設(shè)備和技術(shù)是減少碳排放量積極有效的途徑。

4. 城市化率對(duì)碳排放量的影響不顯著

城市化率對(duì)碳排放量的影響具有雙向性。一方面，城市化率的提高使得城市基礎(chǔ)設(shè)施建設(shè)增加，使得建筑材料需求增加，居民日常生活和供暖等都需要更多的能源消費(fèi)，從而增大碳排放量。另一方面，由于城市化率的提高，使得人口和經(jīng)濟(jì)要素的集聚，能源利用效率提高，以及城市化所帶來(lái)的知識(shí)和技術(shù)等溢出效應(yīng)使碳排放量得以減少。人口結(jié)構(gòu)對(duì)碳排放量的影響盡管不是很顯著，但是人口老齡化現(xiàn)象嚴(yán)重是我國(guó)目前面臨的另外一個(gè)嚴(yán)峻的人口問(wèn)題。因此各級(jí)地方政府在制定人口政策時(shí)應(yīng)立足于本地區(qū)的實(shí)際情況，根據(jù)自身區(qū)域經(jīng)濟(jì)和人口發(fā)展的特點(diǎn)采取相應(yīng)的人口政策、生育政策來(lái)促進(jìn)本區(qū)域乃至全國(guó)人口結(jié)構(gòu)的優(yōu)化。

5. 能源消費(fèi)結(jié)構(gòu)對(duì)碳排放量的影響也不顯著

煤炭消費(fèi)占能源消費(fèi)總量比重的增加會(huì)帶來(lái)碳排放量的增加。其原因在于煤炭的碳排放系數(shù)較大，因此煤炭消費(fèi)所產(chǎn)生的碳排放量比其他能源消費(fèi)所產(chǎn)生的碳排放量大。盡管我國(guó)各地區(qū)的能源消費(fèi)結(jié)構(gòu)相差較為懸殊，但是大部分地區(qū)的能源消費(fèi)以煤炭為主，水能、電能、風(fēng)能等清潔能源所占的比重較低。北京、上海等經(jīng)濟(jì)發(fā)達(dá)地區(qū)的煤炭消費(fèi)占能源消費(fèi)比重在40%左右，而山西、內(nèi)蒙古、貴州、重慶等地區(qū)的煤炭消費(fèi)占能源消費(fèi)比重在80%以上。應(yīng)鼓勵(lì)提高清潔能源在能源消費(fèi)總量中的比重。此外還應(yīng)該積極提高地?zé)岷吞?yáng)能等清潔、可再生能源的開(kāi)發(fā)力度，達(dá)到減少碳排放量，提高環(huán)境質(zhì)量的目的。

[1] Ehrlich P R, Holdren J P. Impact of population growth[J]. Science, 1971, 171(3977): 1212?1217.

[2] Dietz T, Rosa E A. Rethinking the environmental impacts of population, affluence and technology[J]. Human Ecology Review, 1994(1): 277?300.

[3] Mackellar L F, Lutz W, Prinz C, et al. Population, households and CO2emissions[J]. Population and Development Review, 1995, 21(4): 849?865.

[4] York R, Rosa E A, Dietz T. Footprints on the earth the environmental consequences of modernity[J]. American Sociological Review, 2003,68(2): 279?300.

[5] 肖宏偉, 易丹輝, 張亞雄. 中國(guó)區(qū)域碳排放空間計(jì)量研究[J].經(jīng)濟(jì)與管理, 2013, 27(12): 53?62.

[6] 宋德勇, 劉習(xí)平. 中國(guó)省際碳排放空間分配研究[J]. 中國(guó)人口·資源與環(huán)境, 2013, 23(5): 7?13.

[7] Wang C, Chen J, Zou J. Decomposition of energy-related CO2emission in China:1957?2000[J]. Energy, 2005, 30(1): 73?83.

[8] Wu L, Zeng W. Research on the contribution of structure adjustment on carbon dioxide emissions reduction based on LMDI method[J]. Procedia Computer Science, 2013, 17: 744?751.

[9] Ren S, Yin H, Chen X. Using LMDI to analyze the decoupling of carbon dioxide emissions by China's manufacturing industry[J]. Environmental Development, 2014, 9(1): 61?75.

[10] 范丹. 中國(guó)能源消費(fèi)碳排放變化的驅(qū)動(dòng)因素研究——基于 LMDI-PDA分解法[J]. 中國(guó)環(huán)境科學(xué), 2013, 33(9): 1705?1713.

[11] 許海平. 空間依賴、碳排放與人均收入的空間計(jì)量研究[J]. 中國(guó)人口·資源與環(huán)境, 2012, 22(9): 149?157.

[12] Stern D I, Common M S. Is there an environmental Kuznets curve for sulfur?[J]. Journal of Environmental Economics and Management, 2001, 41(2): 162?178.

[13] Galeotti M, Lanza A. Desperately seeking environmental Kuznets[J]. Environmental Modelling & Software, 2005, 20(11): 1379?1388.

[14] Richmond A K, Kaufmann R K. Energy prices and turning points: The relationship between income and energy use/carbon emissions[J]. The Energy Journal, 2006, 27(4): 157?179.

[15] Song T, Zheng T, Tong L. An empirical test of the environmental Kuznets curve in China: A panel cointegration approach[J]. China Economic Review, 2008, 19(3): 381?392.

[16] 王建軍, 周曉唯. 人口、能源消耗、碳排放與經(jīng)濟(jì)增長(zhǎng)關(guān)系實(shí)證分析[J]. 統(tǒng)計(jì)與決策, 2013(23): 114?117.

[17] 杜慧濱, 李娜, 王洋洋, 等. 我國(guó)區(qū)域碳排放績(jī)效差異及其影響因素分析——基于空間經(jīng)濟(jì)學(xué)視角[J]. 天津大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版), 2013, 15(5): 412?416.

[18] 陳德湖, 張津. 中國(guó)碳排放的環(huán)境庫(kù)茲涅茨曲線分析——基于空間面板模型的實(shí)證研究[J]. 統(tǒng)計(jì)與信息論壇, 2012, 27(5): 48?53.

[19] York R, Rosa E A, Dietz T. Bridging environmental science with environmental policy: Plasticity of population, affluence, and technology[J]. Social Science Quarterly, 2002, 83(1): 18?34.

[20] Cole M A, Neumayer E. Examining the impact of demographic factors on air pollution[J]. Population and Environment, 2004, 26(1): 5?21.

[21] Dietz T, Rosa E A. Effects of population and affluence on CO2emissions[J]. Proceedings of the National Academy of Sciences, 1997, 94(1): 175?179.

[22] York R, Rosa E A, Dietz T. Stirpat, ipat and impact: Analytic tools for unpacking the driving forces of environmental impacts[J]. Ecological Economics, 2003, 46(3): 351?365.

[23] 彭希哲, 朱勤. 我國(guó)人口態(tài)勢(shì)與消費(fèi)模式對(duì)碳排放的影響分析[J]. 人口研究, 2010, 34(1): 48?58.

[24] Alvydas B, Tomas B, Dalia S. The energy intensity in Lithuania during 1995－2009: A LMDI approach[J]. Energy Policy, 2011, 39(11): 7322?7334.

[25] 黃芳, 江可申, 盧愿清, 等. 中國(guó)碳強(qiáng)度的影響因素分析——基于LIMDI 分解方法[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí), 2012, 42(6): 40?42.

[26] 李丹丹, 劉銳, 陳動(dòng). 基于空間聚類分析的中國(guó)省域能源消費(fèi)碳排放分布特征研究[J]. 北京師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 49(5): 529?533.

[27] 程葉青, 王哲野, 張守志, 等. 中國(guó)能源消費(fèi)碳排放強(qiáng)度及其影響因素的空間計(jì)量[J]. 地理學(xué)報(bào), 2013, 68(10): 1418?1431.

[28] 呂健. 城市化驅(qū)動(dòng)經(jīng)濟(jì)增長(zhǎng)的空間計(jì)量分析: 2000?2009[J]. 上海經(jīng)濟(jì)研究, 2011(5): 3?15.

[編輯: 譚曉萍]

Analysis of factors affecting carbon emissions based on spatial econometrics

FU Yunpeng, MA Shucai

(School of Economics, Liaoning University, Shenyang 110036, China)

The present study, based on the estimate of the carbon emissions of 30 provinces in China from 2000 to 2014, tests the spatial autocorrelation of carbon emissions by exploiting Moran's I Index. On this basis, the extended STIRPAT Model is taken as the theoretical perspective and spatial regression model is employed in studying the effect of demographic structure, industry structure, energy structure, technological level and other factors on carbon emissions. The result shows that industry structure, demographic structure, and technological level are the main influential factors on carbon emissions, and that the effect of industry structure and demographic size is positive while that of technical level is negative.

demographic size; industry structure; technological level; spatial autocorrelation; spatial regression model

F061.5

A

1672-3104(2017)02?0103?08

2016?11?30；

2017?01?14

國(guó)家社科基金青年項(xiàng)目“基于空間計(jì)量分析的人口規(guī)模、結(jié)構(gòu)對(duì)資源環(huán)境的影響效應(yīng)研究”(13CRK027)；中國(guó)博士后科學(xué)基金第60批面上項(xiàng)目“中國(guó)居民消費(fèi)碳排放量及其影響因素的空間計(jì)量分析”(2016M601327)；遼寧省教育廳科學(xué)研究一般項(xiàng)目“遼寧資源環(huán)境承載力的時(shí)空特征研究”(W2015171)

付云鵬(1978?)，女，滿族，遼寧鐵嶺人，統(tǒng)計(jì)學(xué)博士后，遼寧大學(xué)經(jīng)濟(jì)學(xué)院副教授，主要研究方向：資源經(jīng)濟(jì)學(xué)，計(jì)量經(jīng)濟(jì)模型及其 應(yīng)用