基于夜间灯光数据的中国能源消费碳排放特征及机理_苏泳娴

地理学报

ACTAGEOGRAPHICASINICA

Vol.68,No.11Nov.,2013

基于夜间灯光数据的中国能源消费

碳排放特征及机理

苏泳娴1,3,4,陈修治2,叶玉瑶1,吴旗韬1,张虹鸥1,黄宁生3,匡耀求3(1.广州地理研究所,广州510070;2.中国科学院华南植物园,广州510650;

3.中国科学院广州地球化学研究所,广州510640;

4.中国科学院大学,北京100049)

摘要：本研究基于DMSP/OLS夜间灯光影像实现了1992-2010年以市级为基础单元的我国碳排放估算，弥补了统计数据不全、统计口径不一的缺点。从全国、4个经济区和6大城市群3个层面的碳排放分析结果显示，我国CO2排放总量持续增长，各地区、省市增速各不相同，空间聚集程度越来越明显，基本形成了“东部沿海城市高高集聚，西部欠发达城市低低集聚”的格局。人均碳排放强度基本呈现为“东部>东北部>西部>中部”，单位GDP碳排放强度则呈现为“东北部和西部较高”、“东部和中部较低”。GDP增长是决定CO2排放总量增长的主导因素，而能源结构、能源利用效率、产业结构是影响碳排放强度的主要原因。对于西部和东北部等以能源和重工业为主导产业的城市，其减排策略应着重能源结构优化和能源利用效率的提高。对于东部和中部等以技术、劳动密集型和轻工业为主导产业的城市，其减排策略应侧重于产业结构调整和转型升级。

关键词：碳排放；时空变化；机理；DMSP/OLS夜间灯光影像；遥感DOI:10.11821/dlxb201311007

1引言

碳循环过程是全球气候及环境变化的基本动力源之一。随着陆表系统物质能量交换规模的不断扩大和交换速率的不断加快，碳循环过程的非对称性发展(排放多、吸收少)对全球气候和环境产生越来越明显的负面影响[1]。联合国政府间气候变化专门委员会(IPCC)第四次报告认为，人类活动产生的大规模能源消费所排放的CO2是导致全球气候变化的最主要原因之一[2]。因此，探索经济发展与能源消费等人类活动造成的碳排放成为目前国内外学术界研究的热点问题[3-4]。据统计，2000年以来中国CO2排放量年均增速达10%左右，2008年碳排放总量达83.25亿t，超过美国成为世界第一碳排放大国[5]。巨大的排放总量和较高的增速，使中国面临着巨大的减排压力。由于中国各区域经济社会发展和能源消费结构存在很大差异，因此，科学准确地测算和分析不同区域碳排放的时空变化，揭示其机理问题，是未来明确减排目标、科学制定减排规划的根本前提。

近年来，学者开展了大量的中国碳排放相关研究，主要包含以下几个方面。(1)基于全国能源消费统计数据，测算了国家尺度的能源消费碳排放量、碳排放强度和人均碳排放

收稿日期：2013-04-05;修订日期：2013-10-11

基金项目：十二五”国家科技支撑计划项目(2012BAJ15B02);国家自然科学基金项目(41001385);广东省科技计划项

目(2011B031100003)[Foundation:TheNationalScience&TechnologyPillarProgramduringthe12thFive-yearPlanPeriod,No.2012BAJ15B02;NationalNaturalScienceFoundationofChina,No.41001385;Sci-ence&TechnologyPlanProjectGrant,Guangdongprovince,China,No.2011B031100003]

作者简介：苏泳娴(1985-),女,硕士,助理研究员,主要从事区域生态与区域发展研究。E-mail:suyongxian@163.com通讯作者：张虹鸥(1962-),男,研究员,中国地理学会会员(S110000213M),近年来主要从事城市与区域规划等研究。

E-mail:hozhang@gdas.ac.cn

1513-1526页

1514地理学报68卷

量的变化趋势、空间格局以及不同产业的碳排放差异[1,2,5-11]。(2)结合中国人口、GDP、产业数据、能源消费数据等，分析了影响中国能源消费碳排放量的主要因素[7,12-14]。(3)利用对数平均迪氏分解模型(LogarithmicMeanDivisiaIndex,LMDI)，分解了影响中国能源消费碳排放的主导因素，揭示碳排放影响因素在不同区域的贡献率[15-18]。(4)基于中国碳排放格局、趋势及主要影响因素分析结果，开展了初步的减排策略研究[15,18-20]。另外，部分学者也对省级、区域尺度和部分发达城市的碳排放进行了类似的研究[21-28]。但大部分研究主要依据国家统计局及下属部门发布的统计数据，由于统计数据大多是以国家或省级为单位，市级及更小尺度的统计数据较难收集，给全面了解中国市级及以下尺度碳排放情况带来很大限制。另外，鉴于国家级、省级、市级统计数据的统计口径、计算方法和统计误差等各不相同，三个层面统计数据之间存在较严重的不一致性[28]，也给我国制定准确、系统、差异化的碳减排计划带来较大困难。因此，急需利用新的空间信息获取手段，实现我国碳排放的规范性评估。

遥感技术在时空一致性方面独具优势，逐渐成为地理学领域监测空间信息时空变化的重要手段。美国军事气象卫星DefenseMeteorologicalSatelliteProgram(DMSP)搭载的OperationalLinescanSystem(OLS)传感器，能够有效地探测到城市夜间灯光甚至小规模居民地、车流等产生的低强度夜间灯光，是监测人类活动强度的良好数据源[29]。国内外许多学者已成功地将DMSP/OLS夜间灯光影像应用于城镇扩展监测[30]、社会经济因子估算[31-32]以及能源消耗估算[33]等领域。由于人类活动是碳排放的主要来源，而DMSP/OLS夜间灯光影像恰能有效地反映人类活动强度。因此，DMSP/OLS夜间灯光影像可用于碳排放量的估算，这一结论已被国外学者证明[34-35]。但目前利用DMSP/OLS夜间灯光影像估算碳排放量的研究较少，仍处于起步阶段，研究尺度大多为全球或国家层面，省市级尺度较少，且研究内容仅停留于碳排放总量的探讨。

本文基于1992-2010年DMSP/OLS夜间灯光影像，以TM遥感影像分类结果为验证，精确提取研究区建成区范围，结合30个省或直辖市以及66个地级市的统计数据，建立DMSP/OLS灯光数值与碳排放统计量之间的关系方程。在此基础上，以市级为基础单元，模拟碳排放量，并分别从中国大陆、4大经济区和6大城市群3个层面，分析1992-2010年我国碳排放总量、碳排放密度的时空变化特征，揭示影响碳排放的内在机理，为我国制定差异化的减排策略提供重要参考依据。

2研究区域

选取我国大陆(香港、澳门特别行政区和台湾地区除外)为研究区域(图1)，从市级尺度分析我国CO2排放的时空分布特征。选取东部、中部、西部和东北部4大经济区，分析各区域之间碳排放的差异。另外，近年来逐步成型和迅速发展的6大城市群，集聚了我国大部分的人口和经济活动，是目前碳减排重点关注的地区，因此，本文进一步选取珠三角、辽中南、京津唐、山东半岛、川渝地区、沪宁杭等6大城市群做更深入的研究。

3数据来源和处理

数据来源包括DMSP/OLS夜间灯光影像、LandsatTM影像、全国30个省或直辖市以及66个地级市的统计数据。

(1)DMSP/OLS夜间灯光影像。采用的1992-2010年DMSP/OLS数据，来自美国国家海洋和大气管理局NOAA(NationalOceanicandAtmosphericAdministration)下属的国家地球物理数据中心NGDC(NationalGeophysicalDataCenter)，该数据消除了云及火光等偶然

11期苏泳娴等：基于夜间灯光数据的中国能源消费碳排放特征及机理1515

图1研究区域

Fig.1Locationsofthestudyareas

噪声影响，数据灰度值范围1~63，空间分辨率为0.008333度。由于同一时期不同传感器

获取的夜间灯光数据之间有轻微差异，且存在一些光噪声，本文采用Liu等[10]建立的方法对夜间灯光遥感影像进行融合去噪、切割、相对辐射标定以及地理坐标转换等预处理。

(2)LandsatTM影像。选取覆盖珠三角、长三角、北京、上海、重庆等多个城市的1995、2000、2005、2009、2010年的多幅质量较好的TM遥感影像进行地表类型分类，检验和校正DMSP/OLS夜间灯光影像提取的建成区精度。

(3)统计数据。考虑到采集统计数据的可能性和准确性，选取1995、2000、2005、2009、2010年30个省或直辖市(除澳门、香港、台湾和西藏)及66个地级市(哈尔滨、齐齐哈尔、鸡西、鹤岗、双鸭山、大庆、伊春、佳木斯、七台河、牡丹江、黑河、绥化、太原、大同、阳泉、长治、晋城、朔州、运城、忻州、临汾、武汉、黄石、十堰、宜昌、鄂州、荆门、孝感、咸宁、随州、仙桃、潜江、天门、神农架、南昌、景德镇、萍乡、九江、新余、鹰潭、赣州、吉安、宜春、上饶、厦门、南平、宁德、广州、深圳、珠海、汕头、佛山、韶关、河源、梅州、惠州、东莞、中山、江门、湛江、茂名、肇庆、清远、潮州、揭阳、云浮)的统计年鉴数据。其中用于计算CO2排放量的能源消费数据来自能源消费平衡表中的一次能源终端消费量。

(4)其他辅助数据。主要包括中国国家基础地理信息中心发布的1：400万比例尺的国家、省级、市级的矢量行政界线。

4研究方法

(1)基于DMSP/OLS夜间灯光影像提取建设用地。DMSP/OLS夜间灯光影像与DEM高程数据相似，存在着明显的“峰”和“陷”，在此分别称为“亮源”和“暗源”。亮源与暗源之间存在一个明显的突变区域，即建设用地与非建设用地的分界线。本文采用ArcGIS的NeighborhoodStatistics表面分析工具，借鉴地形起伏度分析方法，以5个栅格为单元，识别建设用地与非建设用地的分界线，提取建设用地范围。利用TM遥感影像分类

1516地理学报68卷

图2DMSP/OLS与TM影像建设用地提取结果对比图

Fig.2ComparisonofurbanareasextractedfromDMSP/OLSdatawiththosefromTMimageries

结果进行验证精度达95.27%，kappa系数0.91。此方法提取的建设用地不会漏掉灯光值比较暗的小城镇建成区，且能够有效提取较宽的道路主干道以及城郊之间分布于道路两侧的居民点。受篇幅所限，仅展示2010年珠三角局部地区的建设用地提取结果(图2)。

(2)基于统计数据计算CO2排放量。采用IPCC公布的《2006温室气体排放清单》来计算能源产生的CO2排放量，公式如下：

9

44CO2=×KiEi(1)12∑i=1

式中：Ei为能源i消费量，按标准煤计，104t；Ki为能源i碳排放系数，(104t碳)/(104t标准煤)；i为能源种类，Ki值来源于IPCC碳排放计算指南缺省值，原始数据单位为J，为与统计数据单位一致，将其转化成标准煤，转化系数为1×104t标准煤等于2.93×105GJ。各种能源的标准煤折算系数和碳排放系数见表1。

(3)建立夜间灯光数值与CO2排放量统计值之间的关系方程。利用建设用地边界切割DMSP/OLS夜间灯光影像，获取各省市建设用地范围内的夜间灯光数据，并统计各省市的夜间灯光总值，然后与相应的CO2排放量统计值进行拟合分析，结果显示(图3和公式2)，夜间灯光数据总值(SDN)与CO2排放统计值具有较好的线性相关关系，在0.01水平上显著相关，R2为0.818。

SDN=0.041×CO2(2)

(4)人均碳排放强度、单位GDP碳排放强度计算方法。为全面了解碳排放格局和机理，本文进一步结合各市人口及经济统计数据，研究了全国人均碳排放强度和单位GDP碳排放强度的时空分布特征及影响机理。

人均碳排放强度=CO2/p(3)单位GDP碳排放强度=CO2/GDP(4)

表1能源碳排放系数表[36-37]

能源种类换算成标准煤

0.71430.97141.42861.47141.47141.45711.4286

(t标准煤/t)碳排放系数

0.75590.8550.58570.55380.57140.59210.6185

(104t碳/104t标准煤)

*注：热力换算成标准煤单位为kg标煤/百万kJ，电力的换算系数单位为kg/kWh。

Tab.1Carbonemissionfactorfordifferenttypesoffuels原煤焦炭原油汽油煤油柴油燃料油

天然气1.330.4483

热力34.12*0.67

电力0.3450.272*

11期苏泳娴等：基于夜间灯光数据的中国能源消费碳排放特征及机理1517

１０００００式中：CO2为二氧化碳排９００００放总量，万t；p采用统计年８００００７００００鉴中的年末常住人口数据，６００００万人。GDP为地区生产总５００００４００００yｙ=0.041x值，换算为以1992年为基期　＝　０．０４１ｘ３００００Ｒ？＝　０．８１８R=0.818２００００的可比价格，万元。

１００００(5)利用Moran'sI指数分０󰀓󰀘󰀓󰀓󰀓󰀓󰀓󰀔󰀓󰀓󰀓󰀓󰀓󰀓󰀔󰀘󰀓󰀓󰀓󰀓󰀓󰀕󰀓󰀓󰀓󰀓󰀓󰀓󰀕󰀘󰀓󰀓󰀓󰀓󰀓析中国CO2排放的空间聚集夜间灯光影像数据总值ＳＤＮ特征。

Moran'sI指数方法是用图3夜间灯光数据总值(SDN)与CO2排放量统计值的拟合关系于统计分析空间分布特征和Fig.3RelationshipbetweentotallightvaluesofDMSP/OLSdataandCO

emissionscalculatedusingstatisticaldata揭示空间分布差异的有效研

究手段[38]。Moran'sI指数包括全局Moran'sI指数和局部Moran'sI指数(LISA)。全局Moran'sI指数能有效反应区域整体集聚情况，值介于-1~1，大于0为正相关，且越接近1，正相关性越强，即邻接空间单元之间具有很强的相似性；小于0为负相关，且越接近-1，负相关性越强，即邻接空间单元之间具有很强的差异性；接近0则表示邻接空间单元不相关。LISA指数能识别局部区域集聚特性。

(6)利用倾向值法(SLOPE)分析中国CO2排放的变化趋势。利用一元线性回归模型计算1992-2010年19年间的变化斜率，即倾向值(SLOPE)，来分析每个城市CO2排放的线性倾向[39]，计算公式为：

能源CO2排放量(104t)22

2(5)æö

n×∑xi2-ç∑xi÷i=1èi=1ø

式中：n是总年份数，等于19；xi是第i年(1992为第一年)；Ci代表第i年对应的CO2

排放量。CO2排放增长类型等级划分标准见表2。

n

n

SLOPE=

n×∑xiCi-∑xi∑Ci

i=1

i=1

i=1

nnn

5结果与分析

5.1CO2排放量模拟值的精度检验

为确保CO2排放量模拟值的可靠性，将其与基于统计数据计算的30个省及直辖市以及66个地级市的CO2排放量进行比照，结果发现(图4)，模拟的CO2排放量与统计数据计算值的均方根误差(RMSE)为943.79万t，相对误差(RE)为7.65%。可见，基于夜间灯光数据模拟所得CO2排放量的精度良好，可用来开展全国CO2排放量的时空特征及机理研究。5.2全国碳排放时空格局

1992-2010年，我国CO2排放总量持续增长，由1992年的14.92亿t增长到2010年的101.23亿吨，年均增速11.22%。总体来看，主要经历4个阶段(图5)：1992-1999年增长较为平缓，年均增长率为7.59%；1999-2002年增速有所加快，年均增长率为12.82%；2002-2007年增长速度再度放缓，年均增长率为4.47%；2007-2010年CO2排放总量迅猛增

表2CO2倾向值(SLOPE)等级划分标准

Tab.2ClassificationcriteriaofCO2emissionSLOPE

增长类型SLOPE值

缓慢增长型

ˉ-0.5s较慢增长型

ˉ-0.5s~Cˉ+0.5sC中速增长型

ˉ+0.5s~Cˉ+1.5sC较快增长型

ˉ+1.5s~Cˉ+2.5sC迅猛增长型

ˉ+2.5s>C

ˉ是指某城市19年CO2排放量平均值；s是指某城市19年CO2排放量标准差。注：C

1518地理学报68卷

长，仅4年间就增加了42.27亿t，相当于过去16年的增加量，年均增长率达19.72%。

本文利用趋势分析法计算了1992-2010年我国各城市的CO2排放总量倾向值，并根据分级标准(表2)，将各城市CO2排放总量增长趋势划分为缓慢增长型、较慢增长型、中速增长型、较快增长型、迅猛增长型等5种类型。结果显示(图6)，全国有12个城市属于迅猛增长型、32个城市属于较快增长型，分别占了国土总面积的2.11%和10.94%，主要集中

分布于东部沿海和经济较发达的图4能源消费CO2排放量统计值与DMSP/OLS影像模拟值散点图

Fig.4Scatterdiagrambetweenenergy-relatedCOemissionscalculated

内陆地区；另外有126个城市属

fromstatisticaldataandCOemissionssimulated

于较慢增长型、168个城市属于fromDMSP/OLSnighttimelightdata缓慢增长型，分别占了国土总面

积的28.33%和48%，主要集中分布于西部地区以及东北部、中部的欠发达地区。

由于各地区CO2排放量的差异化增长，至今已形成新的CO2排放量空间分布格局。Moran'sI空间相关性分析结果显示，1995、2000、2005、2010年全国CO2排放量全局Moran'sI指数分别为0.3848、0.3978、0.4303和0.4781，表明全国CO2排放量分布一直存在正的空间相关性，且空间相关性逐步增强，可见，1992-2010年期间我国CO2排放量在整体上空间聚集现象越来越明显。进一步利用LISA指数对1995、2000、2005、2010年全国CO2排放量进行局部分析发现(图7)。1995年在中部和西部形成了西藏—青海—甘肃—四川—云南—贵阳—广西—湖南—江西—安徽连片的低低集聚区，在海南省形成了另一个低低集聚区。而在张家口—北京—唐山—天津—廊坊—保定、晋中—长治、南通—苏州—泰州—上海、广州—佛山—中山—珠海—东莞—深圳等形成了4个高高集聚区，但此时的高高聚集规模较小，空间分散。兰州、成都、重庆与其周边城市表现出高低邻近的负相关性，承德市、湖州市、清远市则分别与其周边城市表现出低高邻近的负相关性。至2000年，中部和西部连片的低低集聚区逐渐破碎化、规模逐渐减小，而东部和中部的4个高高集聚区规模不断增大，且在潍坊市开始出现新的高高集聚区。至2010年，最终形成了以京津唐城市群为核心的11个城市高高集聚区、以长三角城市群为核心的13个城市高高集聚区、以山东半岛城市群为核心的7个城市高高集聚区、以珠三角城市群为核心的8个城市高高集聚区等；并形成了西藏—甘肃—青海、湖北—湖南—云南—广西、海南岛等3大低低集聚区；兰州、成都、重庆、长沙、武汉分别与其周边城市表现出高低邻近的负相关性，而承德、衡水、日照、宜城、厦门则分别与其周边城市表现出低高邻近的负相关性。碳排放的高高集聚主要发生在区域经济水平较高的东部地区，低低集聚则主要分布在经济发展水平落后的西部和中部及海南岛，高低、低高集聚现象则主要出现在经济发展水平差异较大的相邻城市。可见，区域经济发展水平差异是碳排放空间格局的决定因素。5.3四大经济区碳排放变化趋势

由于不同时期各地区的经济发展水平、政策导向、产业结构的不同，我国东部、中部、西部、东北部四大经济区碳排放量占全国总排放量的比例不断变化(图8)。(1)东部是中国CO2排放的主要地区，1992-1993年，东部地区碳排放在全国的比重大体保持在

模拟的能源CO2排放量(104t)2

2

11期苏泳娴等：基于夜间灯光数据的中国能源消费碳排放特征及机理1519

0100200400600800100012001400160018002000220025003000350040005000600070008000104t

图51992-2010年全国CO2排放总量时空分布图

Fig.5SpatialandtemporaldistributionsofnationalCO2emissionsinChinafrom1992to2010

1520地理学报68卷

图61992-2010年全国CO2排放总量增长等级图

Fig.6ClassificationofCO2emissionincrementspeedsinChinafrom1992to2010

图71995、2000、2005、2010年全国CO2排放总量集聚分布图

Fig.7ClusterconditionsoftheCO2emissionsinChinain1995,2000,2005and2010

50%左右。继1992年邓小平南巡讲话和1994年党的十四大召开，东部地区经济建设进入新的阶段，大量生产要素涌入，使之成为我国人口和经济的聚集地，1994年CO2排放比重也随之猛升到54.28%。1995-2008年，我国经济发展战略逐步转向区域经济的协调发展，东部的碳排放比重也呈稳步略有下降趋势。2008年之后东部地区经济发展重心逐步转为经济结构优化升级和提升经济发展质量，此期间东部的碳排放比重下降至48%。(2)中部

11期苏泳娴等：基于夜间灯光数据的中国能源消费碳排放特征及机理1521

年地区的碳排放比重始终比较稳

２０１０定，基本处于17%~19%之间，间２００９２００８接反映中部地区经济发展速度和２００７００６政策导向处于比较稳定的局面。２２００５(3)西部地区的碳排放比重逐年２００４２００３上升。20世纪90年代初，西部碳２００２００１排放比重仅约为15%，2000年国２２０００家开始实施西部大开发战略，西１９９９１９９８部碳排放比重瞬间增长至１９９７１９９６18.71%，但由于初期资金投入不东部１９９５足和投资效益不高，西部经济并１９９４中部１９９３没得到很快的增长，2000-2005１９９２年期间碳排放比重增长缓慢(仅０％２０％４０％６０％８０％１００％各经济区ＣＯ排放总量占全国的比例２为1.03%)。而“十一五”时期，图81992-2010年四大经济区CO2排放量占全国的比例西部投资进入高速增长阶段，

2009年的碳排放比重明显增长，Fig.8COemissionproportionsofthefourmajoreconomicregions

from1992to2010

较2008年增加了1.65%。(4)由于

90年代以来东北老工业基地逐步衰落，东北地区经济与沿海发达地区差距不断扩大，碳排放量比重不断下降。2003年党中央、国务院正式实施振兴东北老工业基地重大战略决策，2008年其生产总值占全国的比重升至8.62%，增长率为13.4%，碳排放比重同时也呈现转折性的改变。综上可见，GDP是决定地区碳排放水平的主导因素。

1992-2010年期间，四大经济区人均碳排放强度总体呈上升趋势(图9)，基本呈现“东部>东北部>西部>中部”的格局。由于东部地区产业规模大，经济发展程度高，能源消费总量大，人均碳排放强度最大；东北部由于主要发展重工业产业类型，能源利用效率低，消费量大，人均碳排放强度也相对较高。从增长速度来看，四大经济区的人均碳排放强度增长速度越来越快，特别是在2007/2008年之后，出现了明显的增长拐点，其中，东北部地区的人均碳排放强度增长趋势最大，年均增长率达30.05%，其次为西部地区，年均增长率达29.62%，中部地区和东部地区相对较低。这种增长特征与四大经济区的主导

2

６．００　１２．００　１０．００　８．００　单位GDP碳排放强度(t/万元)５．００　４．００　３．００　２．００　１．００　０．００　１３５７９１１１３１５１７１９６．００　４．００　２．００　０．００　东部单位ＧＤＰ碳排放强度东北部单位ＧＤＰ碳排放强度西部人均碳排放强度中部单位ＧＤＰ碳排放强度东部人均碳排放强度东北部人均碳排放强度西部单位ＧＤＰ碳排放强度中部人均碳排放强度图91992-2010年四大经济区人均碳排放强度与单位GDP碳排放强度

Fig.9CO2emissionintensitypercapitaandperGDPofthefourmajoreconomicregionsfrom1992to2010

人均碳排放强度(t)1522地理学报68卷

产业类型有密切关系，由于东北和西部以资源和重工业为主导，相对于以轻工业和制造业为主的东部和中部地区，东北和西部的经济增长带来更快的碳排放增速，而其人口增长速率较东部和中部较缓，导致呈现相对较高的人均碳排放强度。

与人均碳排放强度不同，单位GDP碳排放强度呈现“东北部和西部较高”、“东部和中部较低”的格局(图10)。主要由于东部和中部以轻工业和高端制造业为主，能源效率高，而西部和东北部分别以能源和重工业为主，能源效率相对较低，可见，东部和中部的产业结构较西部和东北部地区更为“低碳”。从时间上看，四大经济区单位GDP碳排放强度总体均逐年减少，特别是2002-2008年期间，各区域单位GDP碳排放强度降幅明显加快，这与中国“十五”期间大力倡导节能减排，积极引入高新技术提高劳动效率与物质能源使用效率有着密切的关系。而2008年之后，东北和西部地区单位GDP碳排放强度突然反弹，这恰好与西部和东北部地区经济快速发展的拐点时间吻合，也间接说明以能源产业和重工业为主导产业的我国西部和东部地区，能源利用效率较低，其经济增长导致单位GDP碳排放强度的增加。

进一步在4大经济区选取典型省市构建了能源消耗总量与碳排放总量的双累计曲线(图10)。以北京和天津为代表的京津唐省市，能源结构一直处于稳定状态，而以广东和上

碳排放累计值(104t)碳排放累计值(104t)碳排放累计值(104t)能源消耗量累计值(104t标准煤)

能源消耗量累计值(104t标准煤)能源消耗量累计值(104t标准煤)

碳排放累计值(104t)碳排放累计值(104t)能源消耗量累计值(104t标准煤)能源消耗量累计值(104t标准煤)

碳排放累计值(104t)能源消耗量累计值(104t标准煤)

碳排放累计值(104t)碳排放累计值(104t)能源消耗量累计值(104t标准煤)能源消耗量累计值(104t标准煤)

碳排放累计值(104t)能源消耗量累计值(104t标准煤)

碳排放累计值(104t)图10四大经济区部分省市能源消耗量与碳排放总量

双累计曲线图

Fig.10Doublemasscurvesbetweenenergyconsumptions

andCO2emissionsofsomeprovincesandcities

inthefourmajoreconomicregions

能源消耗量累计值(104t标准煤)

11期苏泳娴等：基于夜间灯光数据的中国能源消费碳排放特征及机理1523

海为代表的东部沿海省市，１４００００　其能源结构在2005-2006年１２００００　(a)CO2排放总量得到有效改善，碳排放系数明显减少。相反，西部(陕１０００００　西和云南)、东北部(黑龙江和辽宁)、中部(河南和湖８００００　北)在2009-2010年碳排放系６００００　数明显增高，说明其能源结构水平有朝高能源消耗水平４００００　发展的趋势，西部比较明显，东北部次之，中部趋势２００００　最不明显，主要因为西部的０　能源密集型产业开始蓬勃发１９９２１９９４１９９６１９９８２０００２００２２００４２００６２００８２０１０展，东北部的重工业基地开珠江三角洲辽中南京津唐山东半岛川渝地区沪宁杭２０．００　始快速崛起，中部则因为承接着东部地区的产业转移，１８．００　(b)人均碳排放强度引入了大批较高能源消耗水１６．００　平的产业，导致其碳排放系１４．００　数有所增加。１２．００　5.4六大城市群碳排放变化１０．００　趋势８．００　珠三角、辽中南、京津６．００　唐、山东半岛、川渝地区和４．００　沪宁杭等6大城市群的CO2２．００　排放总量、人均碳排放强度０．００　总体呈现增长态势，单位１９９２１９９４１９９６１９９８２０００２００２２００４２００６２００８２０１０GDP碳排放强度则呈现减少１２．００　珠江三角洲辽中南京津唐山东半岛川渝地区沪宁杭趋势(图11)。珠三角城市群(c)单位GDP碳排放强度作为改革开放的先锋地区，１０．００　20世纪90年代碳排放具有８．００　“三高”特点，即高CO2排放量、高人均碳排放强度、高６．００　单位GDP碳排放强度。进入21世纪，随着国家经济发展４．００　战略转移，珠三角经济发展遇到瓶颈，开始寻求高新技２．００　术发展模式，大力开展结构调整和产业升级，提高能源０．００　１９９２１９９４１９９６１９９８２０００２００２２００４２００６２００８２０１０利用效率，此时珠三角CO2排放量增长减缓，单位GDP珠江三角洲辽中南京津唐山东半岛川渝地区沪宁杭碳排放强度明显降低，到图111992-2010年6大城市群CO2排放趋势。2010年其单位GDP碳排放强Fig.11COemissionconditionsofthesixurbanagglomerations度1.49t/万元，是6大城市群

的最低。沪宁杭城市群CO2排放总量一直保持较高的增长率，自1997年起超过珠三角城市

CO2排放强度(t/万元)人均CO2排放量(t)CO2排放总量(104t)2

1524地理学报68卷

群，成为我国CO2排放总量最高的城市群。由于沪宁杭城市群始终以发展低碳的“信息化和高新技术产业”为主，其单位GDP碳排放强度一直处于较低水平。京津唐城市群与依靠“工商业发展、先行对外开放”的珠三角城市群和以“信息化、技术化”为主的沪宁杭城市群发展路径不同，其受全国资源向都城的集中程度影响较大，特别是2008年奥运会的举办，给京津唐城市群带来快速的经济增长，其碳排放总量和人均碳排放强度增长趋势凸显。山东半岛城市群、辽中南城市群和川渝城市群起步较晚，无论是CO2排放总量、人均碳排放强度，还是单位GDP碳排放强度都均处于相对较低的水平。

6结论与讨论

利用DMSP/OLS夜间灯光影像能有效估算我国市级CO2碳排放总量，具有独特的时空连续性，弥补了传统统计数据不全、统计口径不一的缺点，是一种可靠的新型碳排放监测手段。分析结果显示，我国CO2排放总量持续增长，各地区增速各不相同，东部、中部、西部、东北部四大经济区碳排放总量变化与经济发展趋势(GDP)、国家经济政策调整密切相关，六大城市群作为东部、中部、西部、东北部四大经济区域最有代表性的地区，也具有类似的碳排放变化趋势。全国碳排放总量空间聚集程度越来越明显，基本形成了“东部沿海城市高高集聚，西部欠发达城市低低集聚”的格局，人均碳排放强度基本呈现“东部>东北部>西部>中部”，单位GDP碳排放强度则呈现“东北部和西部较高”、“东部和中部较低”的格局，这主要归因于地区能源结构、能源利用效率、产业结构的差异。因此，GDP增长是决定CO2排放总量增长的主导因素，而能源结构、能源利用效率、产业结构是影响碳排放强度的主要原因。

通过调整碳排放强度来实现减排目标是减排的重要手段。由于各地的经济发展水平、产业结构、能源利用水平各不相同，所面临的碳减排瓶颈也不一样，应因地制宜地实施差异化的减排措施。(1)对于西部和东北地区等以能源和重工业为主导产业的城市，能源结构和能源利用效率是影响地区碳排放强度的主导因素，减排策略应该主要放在能源结构优化和能源高效利用新技术的引入。逐步减少煤炭等高碳排放系数能源的使用，大力发展太阳能、风能和地热能等清洁能源，稳步推进新能源基础设施的建设。(2)对于东部和中部地区以技术、劳动密集型和轻工业产业为主导的城市，其减排策略应侧重于产业结构调整，加快淘汰钢铁、有色、化工、建材、电力等高耗能行业的落后生产能力、工艺装备，加快传统高碳产业的升级改造，大力发展高技术为代表的低能耗、低污染、技术密集型产业，增强自主创新能力，不断降低单位产值能耗；着力培育以低碳为特征的产业体系。参考文献(References)

[1]ZhangLei.AchangingpatternofregionalCO2emissionsinChina.GeographicalResearch,2006,25(1):1-9.[张雷.中

国一次能源消费的碳排放区域格局变化.地理研究,2006,25(1):1-9.]

[2]ZhaoRongqin,HuangXianjin,ZhongTaiyang.Researchoncarbonemissionintensityandcarbonfootprintofdifferent

industrialspacesinChina.ActaGeographicaSinica,2010,65(9):1048-1056.[赵荣钦,黄贤金,钟太洋.中国不同产业空间的碳排放强度与碳足迹分析.地理学报,2010,65(9):1048-1056.]

[3]SoytasaU,SariR,EwingBT.Energyconsumption,income,andcarbonemissionsintheUnitedStates.Ecological

Economics,2007,62(3/4):482-489.

[4]ZhuYongbin,WangZheng,PangLietal.SimulationonChina'seconomyandpredictiononenergyconsumptionand

carbonemissionunderoptimalgrowthpath.ActaGeographicaSinica,2009,64(8):935-944.[朱永彬,王铮,庞丽等.基于经济模拟的中国能源消费与碳排放高峰预测.地理学报,2009,64(8):935-944.]

[5]LiuZhancheng,WangAnjian,YuWenjiaetal.ResearchonregionalcarbonemissionsinChina.ActaGeoscientica

Sinica,2010,31(5):727-732.[刘占成,王安建,于汶加等.中国区域碳排放研究.地球学报,2010,31(5):727-732.][6]ZhangLF.AnalysisofcarbonemissioncharacteristicsofChina.PhysicsProcedia,2012,33:1304-1308.

[7]HuCZ,HuangXJ.CharacteristicsofcarbonemissioninChinaandanalysisonitscause.ChinaPopulation,Resources

11期苏泳娴等：基于夜间灯光数据的中国能源消费碳排放特征及机理1525

andEnvironment,2008,18(3):38-42.

[8]YouJ.China'senergyconsumptionandsustainabledevelopment:ComparativeevidencefromGDPandgenuine

savings.RenewableandSustainableEnergyReviews,2011,15:2984-2989.

[9]HuangB,MengLN.ConvergenceofpercapitacarbondioxideemissionsinurbanChina:Aspatio-temporal

perspective.AppliedGeography,2013,40:21-29.

[10]ZhangXP,ChengXM.Energyconsumption,carbonemissions,andeconomicgrowthinChina.Ecological

Economics,2009,68:2706-2712.

[11]GengYH,TianMZ,ZhuQetal.Quantificationofprovincial-levelcarbonemissionsfromenergyconsumptionin

China.RenewableandSustainableEnergyReviews,2011,15:3658-3668.

[12]LiuDU,ChuW,CaiSH.EconomicdevelopmentandcarbondioxideemissionsinChina:Provincialpaneldata

analysis.ChinaEconomicReview,2012,23:371-384.

[13]ZhouXY,ZhangJ,LiJP.IndustrialstructuraltransformationandcarbondioxideemissionsinChina.EnergyPolicy,

2012,http://dx.doi.org/10.1016/j.enpol.2012.07.017.

[14]ZhangYJ.Theimpactoffinancialdevelopmentoncarbonemissions:AnempiricalanalysisinChina.EnergyPolicy,

2011,39:2197-2203.

[15]ZhangY,ZhangJY,YangZFetal.RegionaldifferencesinthefactorsthatinfluenceChina'senergy-relatedcarbon

emissions,andpotentialmitigationstrategies.EnergyPolicy,2011,39:7712-7718.

[16]TianX,ChangM,TanikawaHetal.StructuraldecompositionanalysisofthecarbonizationprocessinBeijing:A

regionalexplanationofrapidincreasingcarbondioxideemissioninChina.EnergyPolicy,2013,53:279-286.

[17]ZhuQ,PengXZ.TheimpactsofpopulationchangeoncarbonemissionsinChinaduring1978-2008.Environmental

ImpactAssessmentReview,2012,36:1-8.

[18]DongJun,ZhangXu.Decompositionofcarbonemissionsandlowcarbonstrategiesforindustrialsectorenergy

consumptioninChina.ResourcesScience,2010,32(10):1856-1862.[董军,张旭.中国工业部门能耗碳排放分解与低碳策略研究.资源科学,2010,32(10):1856-1862.]

[19]ChenWY,WuZX,HeJKetal.CarbonemissioncontrolstrategiesforChina:Acomparativestudywithpartialand

generalequilibriumversionsoftheChinaMARKALmodel.Energy,2007,32:59-72.

[20]WangT,WatsonJ.ScenarioanalysisofChina'semissionspathwaysinthe21stcenturyforlowcarbontransition.

EnergyPolicy,2010,38:3537-3546.

[21]WangYS,MaWC,TuW.AstudyoncarbonemissionsinShanghai2000-2008,China.EnvironmentalScienceand

Policy,2013,27:151-161.

[22]ZhaoGuanwei,ChenJianfei,CuiHaishanetal.CarbonEmissioninenergyconsumptioninGuangZhouduring1992

to2007.ResourcesandIndustries,2012,12(6):179-184.[赵冠伟,陈健飞,崔海山等.1992-2007年广州市能源消费碳排放研究.资源与产业,2012,12(6):179-184.]

[23]LiuJ,FengTT,YangX.TheenergyrequirementsandcarbondioxideemissionsoftourismindustryofWestern

China:AcaseofChengducity.RenewableandSustainableEnergyReviews,2011,15:2887-2894.

[24]NiuXG,NiuJG,MenD.Energyconsumptionandcarbonemissions:Analysisandprediction-thecaseofHebei

provinceinChina.EnergyProcedia,2011,5:2271-2277.

[25]ZhangLeqin,LiRongfu,ChenSupingetal.Trendpredictionandanalysisofdrivingfactorsofcarbonemissionsfrom

energyconsumptionduringtheperiod1995-2009inAnhuiProvincebasedontheSTIRPATmodel.ResourcesScience,2012,34(2):316-327.[张乐勤,李荣富,陈素平等.安徽省1995年-2009年能源消费碳排放驱动因子分析及趋势预测:基于STIRPAT模型.资源科学,2012,34(2):316-327.]

[26]WuBiai,GaoJianhua,XuChong.ThedecompositionofcarbonemissionsinHenanprovincebasedontheindustrial

structureandenergystructure.EconomicGeography,2010,30(11):1902-1907.[吴彼爱,高建华,徐冲.基于产业结构和能源结构的河南省碳排放分解分析.经济地理,2010,30(11):1902-1907.]

[27]ZhangJiekun.FactordecompositionofcarbonemissionsfromenergyconsumptionofShandongProvincebasedon

LMDI.ResourcesScience,2012,34(1):35-41.[宋杰鲲.基于LMDI的山东省能源消费碳排放因素分解.资源科学,2012,34(1):35-41.]

[28]GuanDB,LiuZ,GengYetal.ThegigatonnegapinChina'scarbondioxideinventories.NatureClimateChange

Letters.2012,2:672-675.

[29]ElvidgeCD,CinzanoP,PetDRetal.Thenightsatmissionconcept.InternationalJournalofRemoteSensing,2007,

28(12):2645-2670.

[30]LiuZF,HeCY,ZhangQFetal.ExtractingthedynamicsofurbanexpansioninChinausingDMSP-OLSnighttime

lightdatafrom1992to2008.LandscapeandUrbanPlanning,2012,106:62-72

[31]ChristopherNH,DollJPM,JeremyGM.Mappingregionaleconomicactivityfromnight-timelightsatellite

imagery.EcologicalEconomics,2006,57:75-92.

[32]ChristopherNH,ShonaliPachauri.Estimatingruralpopulationswithoutaccesstoelectricityindevelopingcountries

throughnight-timelightsatelliteimagery.EnergyPolicy,2010,38:5661-5670.

[33]SilvanaA,GilbertoC,AntonioMVMetal.EstimatingpopulationandenergyconsumptioninBrazilianAmazonia

1526

[34][35][36][37][38]

地理学报68卷

[39]

usingDMSPnight-timesatellitedata.Computers,EnvironmentandUrbanSystems,2005,29:179-195.

ElvidgeCD,ImhoffML,BaughKE,etal.Nighttimelightsoftheworld:1994-1995.ISPRSJournalofPhotogrammetryandRemoteSensing,2001,56(2):81-99.

RaupachM,RaynerP,PagetM.Regionalvariationsinspatialstructureofnightlights,populationdensityandfossilfuelCO2emissions.EnergyPolicy,2009,2:61-65.

IPCC.2006IPCCguidelinesfornationalgreenhousegasinventories:volume2[EB/OL].Japan:theInstituteforGlobalEnvironmentalStrategies,2008[2008-07-20].http://www.ipcc.ch/ipccreports/Methodology-reports.htm.

SongJiekun.EmpiricalresearchonfactorsdrivingeconomicgrowthofAnhuiprovince.TechnologyEconomics,2012,31(01):82-85.[宋杰鲲.山东省能源消费碳排放预测.技术经济,2012,31(1):82-85.]

WangQian,WangXiaobin,ZhouRenkang.Spatialdifferencesanditsdrivingfactorsofemergyindicesoncultivatedlandeco-economicsysteminHebeiProvence.ActaEcologicaSinica,2011,31(1):247-256.[王千,金晓斌,周寅康.河北省耕地生态经济系统能值指标空间分布差异及其动因.生态学报,2011,31(1):247-256.]

LuHu,ShiJiancheng.ReconstructionandanalysisoftemporalandspatialvariationsinsurfacesoilmoistureinChinausingremotesensing.ChineseScienceBulletin,2012,57(16):1412-1422.[卢麾,施建成.基于遥感观测的21世纪初中国区域地表土壤水及其变化趋势分析.科学通报,2012,57(16):1412-1422.]

Thecharacteristicsandmechanismsofcarbonemissions

fromenergyconsumptioninChinausingDMSP/OLSnightlightimageries

SUYongxian1,3,4,CHENXiuzhi2,YEYuyao1,WUQitao1,ZHANGHong'ou1,HUANGNingsheng3,KUANGYaoqiu3

(1.GuangzhouInstituteofGeography,Guangzhou510070,China;2.SouthChinaBotanicalGarden,CAS,Guangzhou510650,China;3.GuangzhouInstituteofGeochemistry,CAS,Guangzhou510640,China;3.UniversityofChineseAcademyofSciences,Beijing100049,China)

Abstract:ItiscriticalforChinatomaketheemissionreductiontargetsanddevelopmentofthescientificemissionreductionplanninginthefuture.OnthebasisoftheDMSP/OLSnightlightimageries,thisresearchestimatestheChina'scity-levelcarbonemissionsfrom1992to2010.Thismakesupthevacanciesofstatisticalcarbonemissiondataandovercomestheinconsistenceofstatisticalcarbonemissionmethods.Analysisresultsfromthreescales(thewholemainlandofChina,4economicregionsand6urbanagglomerations)showthatthenationalCO2emissionsgrewcontinually,butvariedfromplacetoplace.Whatismore,thespatialagglomerationofChina'sCO2emissionshasbecomemoreandmoreobvious,whichhaveledtothecurrentCO2emissionpattern—\"high-highconcentrationineasterncoastalcitiesandlow-lowconcentrationinwesternundevelopedcities\".Thecarbonemissionintensityofpercapitabasicallymaintainsasthe\"Eastern>Northeastern>Western>Central\"pattern.ThecarbonemissionintensityofperGDPshowsthe\"HigherinNortheasternandWesternChina\"and\"LowerinEasternandCentralChina\"pattern.GrowthrateofGDPisthemajorfactoraffectingtheincreasingspeedofcarbonemissions.Energystructures,energyefficiencies,industrialstructuresarethethreemainfactorsinfluencingthecarbonintensities.AsfortheWesternandNortheasternregions,whoseindustriesaremainlyenergy-relatedandheavyones,theirbestmitigationpoliciesshouldbeoptimizingtheenergystructureandincreasingenergyuseefficiently.

Keywords:carbonemission;temporalandspatialvariation;mechanism;DMSP/OLSnightlightimagery;remotesensing