宮倩楠，王瑩瑩，童玉芬

(1. 首都經(jīng)濟(jì)貿(mào)易大學(xué) 人口經(jīng)濟(jì)研究所，北京 100070；2.廣東財經(jīng)大學(xué) 經(jīng)濟(jì)學(xué)院，廣東 廣州 510000)

一、問題提出

京津冀地區(qū)是中國北方經(jīng)濟(jì)發(fā)展最快、人口規(guī)模最大、城市群和產(chǎn)業(yè)群最為密集的區(qū)域。以北京、天津和河北的石家莊、保定等13個城市組成的京津冀地區(qū)面積達(dá)到216 000平方公里，2018年末常住人口1.13億，實現(xiàn)生產(chǎn)總值8.51萬億元。該區(qū)域以全國2.25%的地域面積承載了全國8.08%的人口，提供了全國9.46%的經(jīng)濟(jì)總量。近年來，隨著京津冀經(jīng)濟(jì)發(fā)展水平的提高和人口的快速集聚，交通日益擁堵，環(huán)境持續(xù)惡化，資源加速枯竭?！毒┙蚣絽f(xié)同發(fā)展規(guī)劃綱要》明確提出：到2030年，京津冀區(qū)域一體化格局基本形成，區(qū)域經(jīng)濟(jì)結(jié)構(gòu)更加合理，生態(tài)環(huán)境質(zhì)量總體良好，公共服務(wù)水平趨于均衡，成為具有較強(qiáng)國際競爭力、影響力的重要區(qū)域。因此，緩解京津冀地區(qū)較為緊張的人口資源環(huán)境壓力、發(fā)展低碳經(jīng)濟(jì)成為當(dāng)前實現(xiàn)京津冀協(xié)同發(fā)展的重要任務(wù)。

由于自然稟賦、經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、人口要素、消費水平以及能源利用效率等存在較大的空間差異，因此區(qū)域內(nèi)部不同城市地區(qū)的碳排放水平也存在較大差異。人口規(guī)模的擴(kuò)張以及由于生活水平提高導(dǎo)致的居民能源消費模式的改變，使得京津冀地區(qū)實現(xiàn)碳減排目標(biāo)面臨巨大壓力，測算人口規(guī)模變動對碳排放的壓力并分析其空間分布格局的特征具有重要的理論和實踐意義。開展京津冀地區(qū)人口碳排放壓力空間分布格局的研究可以深入了解京津冀地區(qū)的人均碳排放狀況，揭示當(dāng)前人口給京津冀地區(qū)節(jié)能減排帶來的壓力，從而為差別化的低碳城市建設(shè)、人口規(guī)模調(diào)控、碳排放峰值控制提供數(shù)據(jù)支撐，并提高政策制定的針對性和可操作性，為探索不同類型的低碳發(fā)展道路和模式提供指導(dǎo)。

二、文獻(xiàn)綜述

溫室氣體排放帶來全球氣候變暖等一系列的氣候異常問題，在全球范圍內(nèi)給生態(tài)安全、糧食安全和水資源安全等均造成巨大的壓力，嚴(yán)重影響人類社會經(jīng)濟(jì)發(fā)展，節(jié)能減排的低碳經(jīng)濟(jì)成為實現(xiàn)全球可持續(xù)發(fā)展的必然選擇。國內(nèi)外學(xué)術(shù)界紛紛對低碳經(jīng)濟(jì)開展研究，碳排放成為研究重點，其中對碳排放的影響因素的研究較多。學(xué)者們主要關(guān)注經(jīng)濟(jì)增長、人口因素、產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)和效率、國際貿(mào)易、外商直接投資等要素的影響，研究方法多是基于IPAT模型、STIRPAT模型等進(jìn)行相關(guān)性分析，測算各個要素對碳排放規(guī)模的影響程度。其中，經(jīng)濟(jì)要素對碳排放的影響研究，尤其是庫茲涅茲曲線研究較多，如張慶宇等(2019)運用環(huán)境庫茲涅茨模型和STIRPAT模型分析了中國人均GDP和人均碳排放的關(guān)系以及影響碳排放的因素，指出中國人均GDP和人均碳排放的關(guān)系符合環(huán)境庫茲涅茨曲線“倒U型”的規(guī)律，經(jīng)濟(jì)增長是影響中國碳排放的最大要素，其后依次是化石燃料使用、城鎮(zhèn)人口、可替代能源與核使用、人口總數(shù)等[1]。人口要素對碳排放的影響同樣引發(fā)關(guān)注[2-4]，有學(xué)者指出過去20年全球人口變化與二氧化碳排放增長之間存在著較大關(guān)聯(lián)，而且人口變化對發(fā)展中國家碳排放的影響比對發(fā)達(dá)國家更為明顯[5]。朱勤和彭希哲(2010)研究發(fā)現(xiàn)居民消費水平、人口城市化率和人口規(guī)模會顯著影響中國的碳排放總量[6]。余明成(2018)等研究發(fā)現(xiàn)人口對山西的CO2排放影響最大，其次是城鎮(zhèn)化率和第二產(chǎn)業(yè)占GDP比重[7]?？偟膩碚f，當(dāng)前研究主要側(cè)重各要素與碳排放的關(guān)系，其中人口也是影響碳排放的重要因素之一，但關(guān)于人口對碳排放壓力的測定研究相對較少。

分區(qū)域來看，近年來，隨著北京非首都功能疏解、京津冀一體化和生態(tài)文明建設(shè)的推進(jìn)，關(guān)于京津冀地區(qū)碳排放研究的內(nèi)容逐漸增多，主要集中在京津冀地區(qū)碳排放的區(qū)域差異和時空特征、碳排放足跡和強(qiáng)度分析、碳排放的分解及影響因素等方面。具體來看，京津冀地區(qū)碳排放的區(qū)域差異和時空特征方面，汪浩等(2013)基于縣域尺度運用比較統(tǒng)計地圖和空間自相關(guān)分析方法對1990—2009年京津冀都市圈的碳排放空間格局和空間依賴性進(jìn)行研究，發(fā)現(xiàn)京津冀都市圈CO2排放量逐年增加，CO2排放格局呈現(xiàn)較明顯的分層聚集現(xiàn)象，高值區(qū)域空間分布由“雙中心”發(fā)展至“四中心，大部分區(qū)域基本保持不變或減弱，唐山等地CO2排放增強(qiáng)[8]。孫雷剛等(2016)采用碳排放系數(shù)法對2000—2013年環(huán)京津區(qū)域5個城市的碳排放效應(yīng)及時空格局進(jìn)行測算分析，結(jié)果顯示環(huán)京津區(qū)域的城市凈碳排放總量逐年遞增，唐山的凈碳排放總量仍然最高，其次為石家莊、張家口、保定和廊坊[9]。趙玉煥等(2018)也分析了京津冀地區(qū)碳排放的時空特性，發(fā)現(xiàn)河北對京津冀CO2排放的貢獻(xiàn)最大[10]。也有學(xué)者對京津冀地區(qū)部分重點行業(yè)部門的碳排放空間格局進(jìn)行分析，如張克勇(2019)等對2000—2015年京津冀交通部門CO2排放的時空特征進(jìn)行分析，發(fā)現(xiàn)京津冀交通部門CO2增長較快，其中北京對京津冀CO2排放貢獻(xiàn)最大，河北和天津貢獻(xiàn)較小，京津冀各地區(qū)CO2排放差異不斷擴(kuò)大[11]。

碳排放足跡和強(qiáng)度分析方面，閆云鳳(2016)構(gòu)建多區(qū)域投入產(chǎn)出-結(jié)構(gòu)分解分析(MRIO-SDA)模型測度京津冀碳足跡的演變趨勢及影響因素，發(fā)現(xiàn)京津冀碳足跡整體呈現(xiàn)增加趨勢，人口增長、人均消費量增加和消費結(jié)構(gòu)變化是主要的驅(qū)動因素[12]。孫麗文等(2019)利用IPCC法測算了京津冀工業(yè)部門六大高耗能行業(yè)的碳足跡值，并從經(jīng)濟(jì)規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)、能源強(qiáng)度和能源結(jié)構(gòu)四方面對其整體及各行業(yè)進(jìn)行時空差異分析，發(fā)現(xiàn)河北近年碳足跡促增趨勢上升，成為降低碳足跡的重點控制地區(qū)，其中經(jīng)濟(jì)規(guī)模效應(yīng)是首要促增因素[13]。

碳排放的分解及影響因素方面，學(xué)者主要從經(jīng)濟(jì)增長、人口增長、能源消費、產(chǎn)業(yè)轉(zhuǎn)移等不同角度分解碳排放的驅(qū)動因素。如宗剛等(2016)[14]、宋健和趙怡芳(2018)[15]運用對數(shù)平均迪氏指數(shù)法(LMDI)將京津冀三地區(qū)的能源消耗碳排放分解為人口、經(jīng)濟(jì)、產(chǎn)業(yè)結(jié)構(gòu)、能源強(qiáng)度、能耗結(jié)構(gòu)五方面原因，其中人口和經(jīng)濟(jì)對京津冀地區(qū)碳排放均具有拉動作用。閆慶友和尹潔婷(2017)則利用廣義迪氏指數(shù)分解法(GDIM)對京津冀地區(qū)的相關(guān)碳排放數(shù)據(jù)進(jìn)行分析，研究結(jié)果表明經(jīng)濟(jì)發(fā)展對京津冀地區(qū)碳排放的促進(jìn)作用最大，人口增長對碳排放具有拉動作用，但由于人口城鎮(zhèn)化比率和消費習(xí)慣不同，人均碳排放在北京表現(xiàn)為抑制作用，在天津和河北表現(xiàn)為拉動作用[16]。也有學(xué)者基于STIRPAT模型采用回歸等方法對人口要素對碳排放的影響進(jìn)行研究，如趙濤等(2017)研究了1997—2012年京津冀經(jīng)濟(jì)圈人口規(guī)模、人口結(jié)構(gòu)、經(jīng)濟(jì)水平及技術(shù)水平對碳排放的影響，發(fā)現(xiàn)人口規(guī)模、人口年齡結(jié)構(gòu)和人口就業(yè)結(jié)構(gòu)均會促進(jìn)碳排放增加[17]。于洋(2017)則引入人口結(jié)構(gòu)要素，對京津冀城鎮(zhèn)化、人口老齡化與碳排放關(guān)系進(jìn)行實證研究，發(fā)現(xiàn)人口規(guī)模與碳排放呈正相關(guān)關(guān)系，老齡化與碳排放呈顯著的倒U型關(guān)系[18]。

總體來看，關(guān)于人口與碳排放的關(guān)系，現(xiàn)有研究主要將人口要素作為京津冀地區(qū)碳排放影響因素的一方面進(jìn)行分析，雖然也有學(xué)者研究京津冀地區(qū)人口與資源環(huán)境壓力的關(guān)系，但并未將碳排放納入考慮[19]，目前對于京津冀地區(qū)人口對碳排放壓力的直接測定仍然較少。此外，當(dāng)前空間視角下京津冀地區(qū)的碳排放研究仍停留在碳排放規(guī)模的時空變動層面，人口對碳排放壓力的空間格局變動分析較為缺乏。因此，本文將在既有研究的基礎(chǔ)上，計算京津冀地區(qū)各城市區(qū)域的碳排放規(guī)模、測定京津冀地區(qū)面臨的人口碳排放壓力，并基于空間視角分析京津冀地區(qū)人口碳排放壓力的空間格局，進(jìn)一步豐富對京津冀地區(qū)碳排放的相關(guān)研究。

三、研究方法與數(shù)據(jù)處理

(一)京津冀地區(qū)人口碳排放壓力指標(biāo)的構(gòu)建

人口變化對碳排放產(chǎn)生壓力已經(jīng)成為不爭的事實。然而如何衡量人口對碳排放的壓力，學(xué)術(shù)界沒有統(tǒng)一的方法，目前這方面的研究方法主要有以下幾種：

一是因素分解法。采用基于投入產(chǎn)出表的結(jié)構(gòu)分解方法SDA(structural decomposition analysis)和基于擴(kuò)展和改進(jìn)的Kaya恒等式的指數(shù)分解方法IDA(index decomposition analysis)對碳排放進(jìn)行分解，測度人口要素包括人口規(guī)模和結(jié)構(gòu)等對碳排放的貢獻(xiàn)度，其中對數(shù)平均迪氏指數(shù)法(LMDI)使用較多。該類方法可以用來測定人口要素對碳排放值變化的具體貢獻(xiàn)度，但人口被作為社會經(jīng)濟(jì)要素之一進(jìn)行分解研究，并非是對人口碳排放壓力的單獨測算。

二是基于STIRPAT模型，采用面板數(shù)據(jù)回歸估計法、最小二乘法和狀態(tài)空間模型方法等進(jìn)行回歸分析研究人口對碳排放的影響。該類方法雖然可以采用時間序列動態(tài)考察要素變動的影響，仍停留在研究人口要素對碳排放的影響程度上，不能反映對碳排放量變化的具體貢獻(xiàn)，也忽視了不同碳排放單元的空間聯(lián)系和相關(guān)性[20]。

三是采用情景分析方法對不同人口規(guī)模、經(jīng)濟(jì)發(fā)展程度、產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)水平和不同減排政策下的碳排放量進(jìn)行組合，該方法可以在一般均衡條件下模擬各社會經(jīng)濟(jì)要素與碳排放之間的關(guān)系，雖然人口變動也會對不同的碳減排目標(biāo)產(chǎn)生影響，但該方法以模擬預(yù)測為主，參數(shù)設(shè)定具有較大的動態(tài)變化性，并不能實際測量人口對碳排放的壓力。

本文認(rèn)為，構(gòu)建人口碳排放壓力指數(shù)應(yīng)該遵循如下幾個基本原則：(1)簡潔；(2)數(shù)據(jù)易得而準(zhǔn)確；(3)直接反映人口與資源環(huán)境尤其是碳排放的對比關(guān)系；(4)不僅能夠反映人口與資源環(huán)境因素的對比關(guān)系，而且能夠衡量壓力是否存在及其相對大小。鑒于此，本文在已有研究的基礎(chǔ)上，考慮將某地實際的人均碳排放與碳排放的某一個人均標(biāo)準(zhǔn)值進(jìn)行比較，來反映人口對碳排放的壓力。具體測算公式為：

人口碳排放壓力指標(biāo)=(碳排放實際占有量/實際人口規(guī)模)/碳排放的標(biāo)準(zhǔn)人均值

(1)

這里涉及的標(biāo)準(zhǔn)人均值，可以是可利用的人均標(biāo)準(zhǔn)或公認(rèn)的國際或國內(nèi)標(biāo)準(zhǔn)，也可以是全國或地區(qū)的平均值。本文根據(jù)碳排放的自身特點，選取全國平均的人均碳排放量作為標(biāo)準(zhǔn)，京津冀各城市區(qū)域的人均碳排放量越大，其人口碳排放壓力也就越大。

(二)數(shù)據(jù)來源及處理

本文數(shù)據(jù)均來自2005—2015年北京、天津、河北11個城市的城市統(tǒng)計年鑒。

關(guān)于碳排放的計算，主要采用IPCC排放系數(shù)計算，其中河北11個地級市由于缺少分品種能源消耗數(shù)據(jù)，無法采用一般的碳排放計算方法，因此本文通過對2005—2015年河北的能源消耗總量和碳排放總量關(guān)系擬合，估算出二者之間的擬合方程，從而利用河北省11個地級市2005—2015年的能源消耗量估算出其碳排放量。而北京和天津則根據(jù)煤品、油品以及天然氣的消耗量乘以IPCC排放系數(shù)加總得到。

圖1 河北碳排放兩種計算方法對比

從圖1可以看出，根據(jù)IPCC排放系數(shù)計算的河北碳排放數(shù)量與根據(jù)能源消費估算出的碳排放數(shù)據(jù)極為接近，而河北能源消費總量等于各個地級市的能源消費量之和，河北碳排放總量等于各個地級市的碳排放量的加總，因此通過擬合得到河北總能源消費與碳排放之間的數(shù)量關(guān)系，可以推導(dǎo)出各個地級市能源消費量與碳排放量之間的數(shù)量關(guān)系，即通過上述擬合得到河北總的能源消費與碳排放之間的數(shù)量關(guān)系為：

C=0.729 5E+85.497

(2)

(3)

由此得到每個城市的能源消費與碳排放之間的數(shù)量關(guān)系式ci=0.729 5ei+85.497,從而估算出的各地級市各年份碳排放數(shù)據(jù)，由于篇幅有限，僅列出部分年份數(shù)據(jù)，如表1所示。

表1 2005—2015年京津冀各城市區(qū)域碳排放規(guī)模 單位：萬噸標(biāo)準(zhǔn)煤

注：資料來自《北京統(tǒng)計年鑒2016》《天津統(tǒng)計年鑒2016》《河北經(jīng)濟(jì)統(tǒng)計年鑒2016》;其中北京市和天津市數(shù)據(jù)根據(jù)IPCC系數(shù)計算。

四、京津冀地區(qū)人口規(guī)模對碳排放壓力的空間格局分析

本文依據(jù)上文計算得出的京津冀13個城市2005—2015年的碳排放量、常住人口規(guī)模以及中國人均碳排放量的標(biāo)準(zhǔn)值，采用式(1)進(jìn)一步測算出京津冀地區(qū)的人均碳排放壓力狀況如表2所示。圖2和圖3分別表示2005—2015年京津冀地區(qū)人口碳排放壓力和2015年京津冀地區(qū)人碳排放壓力排序。

表2 2005—2015年京津冀地區(qū)人口碳排放壓力值

圖2 2005—2015年京津冀地區(qū)人口碳排放壓力

圖3 2015年京津冀地區(qū)人口碳排放壓力排序

(一)京津冀地區(qū)內(nèi)部人口碳排放壓力的時空演變分析

首先，總體變化狀況顯示，2005—2015年京津冀地區(qū)長期面臨人口碳排放壓力，但壓力程度近年來呈下降趨勢。2005—2015年京津冀地區(qū)的人口碳排放壓力值始終大于1，長期高于全國人均碳排放水平，說明京津冀地區(qū)仍面臨著一定程度的人口碳排放壓力。但從變動上看，從2005—2015年，京津冀整體的人口碳排放壓力值從1.670下降至1.336，下降幅度達(dá)到25%。內(nèi)部各區(qū)域人口碳排放壓力的最大極值也從2005年的5.111下降至3.839。此外，與2005年相比，2015年京津冀人口碳排放壓力值低于1(即低于全國人均碳排放水平)的城市數(shù)量由4個增加至6個。雖然大多數(shù)區(qū)域仍面臨著人口碳排放壓力，但京津冀地區(qū)人口碳排放壓力相對來說已經(jīng)得到緩解。

與全國水平相比，2005年、2010年和2015年京津冀人均碳排放量分別為2.23噸標(biāo)準(zhǔn)煤、2.56噸標(biāo)準(zhǔn)煤和2.63噸標(biāo)準(zhǔn)煤，2005—2015年人均碳排放量的增幅為18%，而全國2005年和2015年的人均碳排放規(guī)模分別為1.34噸標(biāo)準(zhǔn)煤和1.97噸標(biāo)準(zhǔn)煤，增幅達(dá)到47%，京津冀地區(qū)人均碳排放規(guī)模的增長幅度遠(yuǎn)低于全國增長幅度，說明近年來京津冀地區(qū)的節(jié)能減排工作取得了一定的進(jìn)展。京津冀地區(qū)作為中國北方最大的城市群，始終面臨著嚴(yán)重的環(huán)境問題，為了有效推進(jìn)生態(tài)環(huán)境建設(shè)，京津冀地區(qū)推動區(qū)域產(chǎn)業(yè)結(jié)構(gòu)和能源結(jié)構(gòu)的優(yōu)化調(diào)整，調(diào)控人口規(guī)模、控制溫室氣體排放，協(xié)同推進(jìn)區(qū)域大氣污染治理，從多角度出發(fā)實現(xiàn)京津冀地區(qū)的生態(tài)環(huán)境保護(hù)，京津冀地區(qū)面臨的人口碳排放壓力逐步得到緩解。

其次，從京津冀內(nèi)部看，2005—2015年京津冀地區(qū)的人口碳排放壓力存在較大的空間差異性。圖3顯示，當(dāng)前2015年京津冀地區(qū)人口碳排放壓力值超過1的城市地區(qū)共計7個，分別為唐山、邯鄲、天津、承德、秦皇島、張家口和石家莊，其中唐山的人口碳排放壓力值最大，達(dá)到3.839，遠(yuǎn)高于排名第二和第三的邯鄲與天津。唐山2016年常住總?cè)丝谝?guī)模為784.36萬人，人口密度達(dá)到582人/平方公里；地區(qū)生產(chǎn)總值居全省首位，為6 354.87億元；雖然單位GDP能耗從2010年的1.638噸標(biāo)準(zhǔn)煤/萬元下降至2016年的1.252 9萬噸標(biāo)準(zhǔn)煤/萬元，但全社會能源消費量從2010年的7 320.93萬噸標(biāo)準(zhǔn)煤上升至8 163.4萬噸標(biāo)準(zhǔn)煤，較大的人口規(guī)模和較高的能源消費量使得唐山面臨著巨大的人口碳排放壓力。其他存在人口碳排放壓力的城市壓力值相差不大，基本分布在[1，2]的區(qū)間內(nèi)。整體來看，當(dāng)前面臨人口碳排放壓力的城市區(qū)域除邯鄲外，其他主要集中在京津冀地區(qū)的東北部。而目前壓力值小于1的城市地區(qū)共計6個，分別為滄州、邢臺、廊坊、北京、衡水和保定，主要集中在京津冀地區(qū)的中南部。其中保定的人口碳排放壓力最小，雖然人口規(guī)模較大，但保定較早成為中國低碳城市發(fā)展的試點城市，新能源產(chǎn)業(yè)發(fā)展迅速，以“中國電谷”和“太陽能之城”享譽(yù)海外，逐漸走上了生態(tài)文明建設(shè)下的綠色發(fā)展道路。

最后，從京津冀內(nèi)部人口碳排放壓力變化來看，也呈現(xiàn)出明顯的差異特征。2005—2015年，承德、廊坊以及滄州人口碳排放壓力上升，但其他多數(shù)城市的人口碳排放壓力呈現(xiàn)持續(xù)下降的趨勢。2005—2015年以來，京津冀地區(qū)人口碳排放壓力的變化趨勢主要存在兩種類型：

一種是持續(xù)上升型，如承德、廊坊和滄州的人口碳排放壓力均在2015年有明顯增加，增長幅度分別約為27.41%、24.12%和51.65%。其中承德的人口碳排放壓力值從2015年的1.222上升至1.557，人均碳排放量達(dá)到全國人均水平的1.5倍，承德近年來人口規(guī)模相對變動不大，這主要與當(dāng)前承德產(chǎn)業(yè)結(jié)構(gòu)比重不合理、能源消費結(jié)構(gòu)仍以高污染的煤炭為主、高耗能行業(yè)較多有關(guān)，如果繼續(xù)提高的話，承德將面臨較大的碳排放壓力，因此提高能源資源利用效率、降低碳強(qiáng)度成為發(fā)展方向。而廊坊和滄州兩個城市的人口碳排放壓力盡管也在持續(xù)上升，但目前仍低于1，處于相對無壓區(qū)。

另一種是持續(xù)下降型。2005—2015年北京、天津、石家莊、張家口、秦皇島、唐山、保定、衡水、邢臺和邯鄲呈下降態(tài)勢，其中北京、衡水和石家莊年均下降率較高，分別為8.5%、4.76%和4.31%。尤其北京成為重要的人口碳排放壓力下降拉動區(qū)，人口碳排放壓力值已持續(xù)下降到1以下，這主要與北京開展人口規(guī)模調(diào)控、非首都功能疏解和生態(tài)文明城市建設(shè)等有關(guān)。此外，唐山、天津和邯鄲的人口碳排放壓力也得到較為明顯的緩解，人口碳排放壓力值均在原來基礎(chǔ)上下降1～2個等級，如唐山由2005年的5.21下降到2015年的3.84，天津和邯鄲壓力值也由大于2下降到2以下，這表明隨著京津冀地區(qū)社會經(jīng)濟(jì)發(fā)展水平的提高、技術(shù)的改進(jìn)、產(chǎn)業(yè)結(jié)構(gòu)和能耗結(jié)構(gòu)日趨合理化以及對節(jié)能減排的重視，碳排放得到一定程度的抑制。

本文借鑒王韶華(2018)[21]等對于降耗核心-外圍區(qū)的分析思路，進(jìn)一步將京津冀各城市區(qū)域劃分為降壓拉動區(qū)、降壓中間區(qū)和降壓滯后區(qū)。其中降壓拉動區(qū)為人口碳排放壓力明顯低于京津冀地區(qū)整體人口碳排放壓力、對京津冀地區(qū)人口碳排放壓力的下降起幫助作用的地區(qū)；降壓滯后區(qū)為人口碳排放壓力大于京津冀地區(qū)整體人口碳排放壓力、對京津冀地區(qū)人口碳排放壓力起明顯抑制作用的地區(qū)；處于兩者之間的地區(qū)稱為降壓中間區(qū)。降壓拉動區(qū)、中間區(qū)和滯后區(qū)的判斷用人口碳排放壓力的標(biāo)準(zhǔn)化值為標(biāo)準(zhǔn)，計算公式為：

(4)

依照式(4)可分別得到2005年、2010年和2015年京津冀13個城市區(qū)域人口碳排放壓力的標(biāo)準(zhǔn)化值，按照降壓拉動區(qū)、降壓中間區(qū)和降壓滯后區(qū)的標(biāo)準(zhǔn)進(jìn)行聚類分析，如表3所示。2005—2015年，降壓拉動區(qū)的城市數(shù)量不斷減少，只有保定、衡水、北京和廊坊的人口碳排放壓力標(biāo)準(zhǔn)化值始終穩(wěn)定在0.5以上，成為京津冀地區(qū)的降壓拉動區(qū)。降壓滯后區(qū)則較為穩(wěn)定，唐山、天津和邯鄲長期處于降壓滯后區(qū)，尤其是唐山和天津人口碳排放壓力標(biāo)準(zhǔn)值均在-2以上，可見對緩解京津冀人口碳排放壓力產(chǎn)生較大的負(fù)向作用的區(qū)域主要位于京津冀東北部地區(qū)。降壓中間區(qū)則由2005年的3個城市增加到2015年的6個城市，包括邢臺、滄州、石家莊、張家口、秦皇島、承德，其中滄州、邢臺和廊坊從降壓拉動區(qū)落至降壓中間區(qū)，各區(qū)域人口碳排放壓力的差距整體呈縮小態(tài)勢。

表3 京津冀地區(qū)人口碳排放降壓聚類分析

(二)京津冀地區(qū)人口碳排放壓力的空間相關(guān)性分析

地理學(xué)第一定律認(rèn)為任何空間事物及其屬性都相關(guān)，在分析京津冀地區(qū)人口碳排放壓力空間格局演化的過程中，如果忽略空間自相關(guān)性，將碳排放單元視為相互獨立的均質(zhì)個體，勢必會對研究結(jié)果造成影響或偏差。許多學(xué)者認(rèn)為空間差異和地理鄰近是決定區(qū)域碳排放的重要因素，碳排放存在空間溢出性，人口碳排放壓力的變動也應(yīng)考慮人口規(guī)模變動和碳減排活動對其的直接影響以及對鄰近區(qū)域的空間溢出影響。因此，本文根據(jù)2005—2015年京津冀各城市區(qū)域的統(tǒng)計年鑒，利用GeoDA空間計量軟件使用一階Rook’s準(zhǔn)則構(gòu)建空間權(quán)重矩陣，繪制出莫蘭(Moran)散點圖和LISA集聚圖，對京津冀地區(qū)人口碳排放壓力進(jìn)行空間相關(guān)性分析，并將空間集聚格局可視化，進(jìn)一步檢測京津冀地區(qū)人口碳排放壓力的“熱點”區(qū)和“冷點”區(qū)，從而結(jié)合當(dāng)前京津冀協(xié)同發(fā)展問題，為指導(dǎo)京津冀地區(qū)發(fā)展低碳經(jīng)濟(jì)、科學(xué)制定節(jié)能減排政策、實現(xiàn)各城市區(qū)域人口分布與功能布局的優(yōu)化、推動京津冀一體化提供依據(jù)。

本文采用的空間自相關(guān)分析方法根據(jù)研究范圍的大小，可具體劃分為全局空間自相關(guān)和局部空間自相關(guān)。其中全局指標(biāo)用于反映整個京津冀地區(qū)人口碳排放壓力的空間關(guān)聯(lián)模式；局部指標(biāo)則有助于研究京津冀內(nèi)部局部地區(qū)的非典型性特征。

1.全局空間自相關(guān)

本研究采用莫蘭指數(shù)(Moran’sI)來衡量全局空間相關(guān)性，可用式(5)表示：

(5)

式(5)中，n為研究區(qū)域的空間樣本個數(shù)，xi和xj表示空間樣本單元的相應(yīng)屬性值，Wij是空間權(quán)重系數(shù)矩陣，表示空間單元的鄰近關(guān)系，本研究采用基于共同邊界的一階Rook權(quán)重矩陣。Moran’sI的值集中在[-1,1]區(qū)間內(nèi)，值越接近1，則表明具有相似屬性的空間單元產(chǎn)生集聚；值越接近-1，則表明產(chǎn)生集聚的空間單元具有相異屬性；如果值接近或等于0，則表明空間單元屬性處于隨機(jī)分布狀態(tài)。此外，對全局空間自相關(guān)分析的結(jié)果需要進(jìn)行Z值顯著性檢驗，根據(jù)Z值大小，在設(shè)定顯著性水平下做出接受或拒絕零假設(shè)的判斷[22]。

全局自相關(guān)分析結(jié)果顯示，2005年、2010年和2015年京津冀地區(qū)人口碳排放壓力的全局莫蘭指數(shù)分別為0.089 7,0.163 7,0.200 2，總體來看，有所增加，空間“集中”的趨勢逐漸加強(qiáng)。莫蘭散點圖中的“高-高”(HH)和“低-低”(LL)象限表明人口碳排放壓力的觀測值存在較強(qiáng)的空間正相關(guān)，即存在均質(zhì)性；“低-高”(LH)和“高-低”(HL)關(guān)聯(lián)表示存在較強(qiáng)的空間負(fù)相關(guān)，即空間單元存在異質(zhì)性。結(jié)合人口碳排放壓力的莫蘭散點圖來看，2005年和2010年京津冀仍有部分地區(qū)位于第二和第四象限，呈現(xiàn)高低、低高集聚現(xiàn)象；2015年，京津冀人口碳排放壓力的空間正相關(guān)性明顯增強(qiáng)，絕大多數(shù)區(qū)域均處于第一和第三象限，呈現(xiàn)出顯著的低低和高高集聚態(tài)勢。當(dāng)前京津冀區(qū)域人口碳排放壓力存在明顯的空間正相關(guān)關(guān)系，并且這種空間正相關(guān)關(guān)系隨著年份的推移不斷增強(qiáng)，人口碳排放壓力的空間集聚性愈發(fā)顯著。2005—2015年京津冀區(qū)域人口碳排放壓力的全局莫蘭散點圖如圖4所示。

圖4 2005年、2010年和2015年京津冀地區(qū)人口碳排放壓力的莫蘭散點圖

2.局部空間自相關(guān)

(6)

京津冀地區(qū)人口碳排放壓力的局部自相關(guān)結(jié)果(見表4)表明，2005年京津冀人口碳排放壓力的空間集聚狀況不太明顯，13個城市區(qū)域中僅有1個地區(qū)呈現(xiàn)低高集聚現(xiàn)象。承德2005年人口碳排放壓力值為1.222，位于京津冀人口碳排放壓力的低值區(qū)，周邊的張家口、北京、唐山和秦皇島的人口碳排放壓力均高于承德市。承德位于京津冀北部的水源涵養(yǎng)區(qū)，以山地地形為主，森林覆蓋率較高，人口密度相對較低，產(chǎn)業(yè)占比第一產(chǎn)業(yè)較高，自然稟賦和經(jīng)濟(jì)稟賦均區(qū)別于緊鄰的北京、天津和唐山等大城市地區(qū)，這也反映出2005年京津冀地區(qū)的核心城市如北京、天津等均處于人口碳排放壓力的集聚中心。2010年京津冀人口碳排放壓力的低值集聚區(qū)增加為2個，除承德外，保定也處于人口碳排放壓力的“冷點”區(qū)域，而鄰近的石家則處于人口碳排放壓力的“熱點”區(qū)域，呈現(xiàn)高低集聚態(tài)勢。這說明河北中部城市如保定、衡水和邢臺的人口碳排放壓力在下降，人口碳排放壓力向省會城市集聚。2015年京津冀的人口碳排放壓力空間集聚模式以高高集聚和低低稀疏為主。東北部的唐山和秦皇島兩個區(qū)域呈現(xiàn)高高集聚態(tài)勢，說明人口碳排放壓力正在向東北地區(qū)集聚，唐山、天津、承德和秦皇島面臨的人口碳排放壓力較大；而西南地區(qū)的保定、石家莊呈現(xiàn)出低低集聚態(tài)勢。這說明石家莊、保定、衡水和邢臺等河北中部地區(qū)以及廊坊、北京等面臨的人口碳排放壓力整體有所緩解。

表4 2005年、2010年、2015年京津冀人口碳排放壓力局域空間自相關(guān)顯著性結(jié)果

(三)京津冀地區(qū)人口碳排放壓力的重心演變及空間分布

標(biāo)準(zhǔn)差橢圓(standard deviational ellipse，SDE)是空間分析中解釋地理要素空間分布特征的有效方法，該方法可以從重心、展布范圍、密集性、方向和形狀等多角度解釋地理要素的空間分布整體特征及其時空演化過程。具體來看，標(biāo)準(zhǔn)差橢圓中心反映經(jīng)濟(jì)要素空間分布整體在二維空間上的相對位置(重心)；長軸的方向(方位角，即正北方向與順時針旋轉(zhuǎn)的長軸之間的夾角)反映其在二維空間上展布的主趨勢方向；長、短軸的比值可以體現(xiàn)要素空間分布的形態(tài)，單位標(biāo)準(zhǔn)差橢圓上分布的空間要素總量體現(xiàn)其在二維空間上展布的密集程度?；诘乩硇畔⒌臉?biāo)準(zhǔn)差橢圓已經(jīng)成為ArcGIS空間統(tǒng)計模塊的常規(guī)統(tǒng)計工具[23]。由于京津冀地區(qū)經(jīng)濟(jì)發(fā)展不平衡、資源分布差異大，人口規(guī)模變動與流動狀況存在差別，人口對碳排放壓力處于不對稱狀態(tài)。采用標(biāo)準(zhǔn)差橢圓方法的重心、橢圓范圍和長短軸等工具，能夠明確當(dāng)前京津冀地區(qū)人口碳排放壓力的重心分布，并有效說明京津冀地區(qū)人口碳排放壓力集聚重心趨勢的變化，為推動京津冀地區(qū)一體化和低碳化發(fā)展提供幫助。

表5 2005—2015年京津冀人口碳排放壓力的重心坐標(biāo)

年份經(jīng)度X緯度Y2005116.62639.1542010116.62739.1662015116.69339.181

本文采用ArcGIS軟件，基于京津冀13個城市的空間區(qū)位(經(jīng)緯度)，用各城市區(qū)域的人口碳排放壓力表示相應(yīng)的權(quán)重，計算京津冀人口碳排放壓力的標(biāo)準(zhǔn)差橢圓，制得表5和表6。表5顯示，2005—2015年京津冀地區(qū)人口碳排放壓力重心的空間軌跡位移不大，緩慢向東北方向移動。在南-北方向上，2005—2010年，京津冀地區(qū)壓力重心呈現(xiàn)明顯的向北移動趨勢；2010—2015年，壓力重心繼續(xù)向北移動；在東-西方向上，京津冀地區(qū)人口碳排放圧力重心基本保持向東移動的趨勢，且2010—2015年向東位移最大。最終，京津冀地區(qū)人口碳排放壓力的重心從2005年的(39.154N,116.626E)移動至2015年的(39.181N,116.693E)，位置落在河北霸州。總體來看，京津冀地區(qū)人口碳排放的壓力逐漸向東北部方向偏移，東北方向的人口碳排放壓力在逐漸增強(qiáng)。

表6 2005—2015年京津冀地區(qū)人口碳排放壓力空間分布的橢圓特征值

結(jié)合表6數(shù)值來看，2015年，京津冀人口碳排放壓力標(biāo)準(zhǔn)差橢圓的方位角(北方向與順時針旋轉(zhuǎn)的長軸之間的夾角)為51.129度，長軸標(biāo)準(zhǔn)距離遠(yuǎn)大于短軸，表明京津冀人口碳排放壓力的空間分布趨勢整體呈現(xiàn)為東北-西南方向。

從時間尺度上看，2005年以來京津冀人口碳排放壓力的空間分布形狀呈扁化趨勢，人口碳排放壓力逐漸向東北地區(qū)集中。表6人口碳排放壓力空間分布的橢圓在京津冀區(qū)域的空間覆蓋面積標(biāo)準(zhǔn)值顯示， 2005—2015年，人口碳排放壓力的空間覆蓋面積整體上保持縮減趨勢，京津冀人口碳排放壓力的空間分布在橢圓范圍內(nèi)逐漸集中。進(jìn)一步分析長短軸距離，2005—2015年，京津冀人口碳排放壓力空間分布標(biāo)準(zhǔn)差橢圓的長軸距離在增加，短軸距離在縮短，長軸和短軸的比值不斷增大。這表明，京津冀地區(qū)的人口碳排放壓力在東北-西南方向上呈擴(kuò)張狀態(tài)，在西北-東南方向上呈收縮狀態(tài)，且東北-西南向的擴(kuò)張趨勢強(qiáng)于西北-東南方向的收縮。此外，標(biāo)準(zhǔn)差橢圓的方位角呈現(xiàn)逆時針轉(zhuǎn)動，方位角的度數(shù)由2005年的53.385度縮小到2015年的51.129度，表明相對于西南部地地區(qū)，位于東北部地區(qū)的城市區(qū)域?qū)θ丝谔寂欧艍毫Φ睦瓌幼饔迷鰪?qiáng)。

最后，結(jié)合城市區(qū)域看，人口碳排放壓力的重心向東北地區(qū)移動、標(biāo)準(zhǔn)差橢圓呈現(xiàn)東北-西南方向的收縮，唐山和天津發(fā)揮了重要的拉動作用。尤其唐山2005年、2010年和2015年人口碳排放壓力值分別達(dá)到5.111、4.102和3.839，雖然整體呈下降趨勢，但歷年來均是京津冀地區(qū)人口碳排放壓力最大的城市區(qū)域，這主要與唐山產(chǎn)業(yè)結(jié)構(gòu)偏向于重工業(yè)、高污染高耗能企業(yè)多、能源消費結(jié)構(gòu)不合理、科技水平相對落后有關(guān)。天津2015年人口碳排放壓力值為1.835，僅次于唐山，天津作為京津冀地區(qū)僅次于北京的特大城市，經(jīng)濟(jì)發(fā)展主要依靠工業(yè)帶動，能源結(jié)構(gòu)優(yōu)化和產(chǎn)業(yè)結(jié)構(gòu)的轉(zhuǎn)型升級還在進(jìn)行，人口規(guī)模較大，面臨的人口碳排放壓力仍較為嚴(yán)重。此外，北部承德人口碳排放壓力2015年超過張家口，達(dá)到1.5以上，面臨的壓力在逐漸增大。相對應(yīng)地，京津冀西南區(qū)域的人口碳排放壓力在緩解，石家莊的人口碳排放壓力值從1.968下降至1.291，長期處于人口碳排放高壓力區(qū)的邯鄲2015年人口碳排放壓力值也下降至2以下，西南地區(qū)的人口碳排放壓力整體呈現(xiàn)下降趨勢，人口碳排放壓力總體向東北地區(qū)集中。

五、結(jié)論及建議

本文在構(gòu)建人口碳排放壓力指標(biāo)的基礎(chǔ)上，進(jìn)一步采用聚類分析、空間相關(guān)性分析和標(biāo)準(zhǔn)差橢圓等方法，研究了京津冀地區(qū)人口對碳排放壓力的空間格局及變動，在此基礎(chǔ)上提出以下結(jié)論及相應(yīng)的政策建議：

第一，當(dāng)前京津冀地區(qū)內(nèi)部人口碳排放壓力存在較大的空間差異性，唐山、天津和邯鄲面臨著較為嚴(yán)重的人口碳排放壓力，北京、衡水和保定的人口碳排放壓力值較低，處于無壓區(qū)。總體來看，人口碳排放壓力的整體空間分布呈東北-西南主體方向，壓力重心進(jìn)一步向東北方向集聚。

針對當(dāng)前人口碳排放壓力水平較高的地區(qū)如唐山、天津和邯鄲等，應(yīng)當(dāng)調(diào)整能源消費結(jié)構(gòu)，增加清潔能源和新能源的使用，減少對傳統(tǒng)能源的依賴程度；積極轉(zhuǎn)變經(jīng)濟(jì)發(fā)展方式、優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，降低高能耗產(chǎn)業(yè)的比重；加大科技創(chuàng)新投入，提高能源資源的利用效率，促進(jìn)城市發(fā)展向綠色低碳方向轉(zhuǎn)型。尤其天津作為京津冀城市群的雙核之一，應(yīng)當(dāng)利用制造業(yè)優(yōu)勢與北京形成研發(fā)和人力資本對接，向上延伸制造業(yè)產(chǎn)業(yè)鏈，提高制造業(yè)的附加價值，推動制造業(yè)的精細(xì)化、高端化和創(chuàng)新化，并大力發(fā)展第三產(chǎn)業(yè)，促進(jìn)產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型。而北京、衡水和保定等雖然處于無壓區(qū)，但距離低碳城市仍有較大發(fā)展空間，這些地區(qū)應(yīng)當(dāng)進(jìn)一步加大節(jié)能減排力度，繼續(xù)淘汰過剩產(chǎn)能，提高全民節(jié)能意識。尤其北京要加強(qiáng)與京津冀其他區(qū)域的協(xié)作，充分利用自身優(yōu)勢，包括充足的資金投入、較強(qiáng)的技術(shù)創(chuàng)新能力和豐富的科研單位、高校和企業(yè)資源等，推進(jìn)節(jié)能和減排技術(shù)的研發(fā)推廣。

第二，京津冀目前人口碳排放壓力普遍減弱，但從時間尺度上看內(nèi)部變動趨勢存在差異，一些生態(tài)比較脆弱或者碳排放比較嚴(yán)重的城市如承德、廊坊和滄州等人口碳排放壓力不降反升，應(yīng)該引起高度重視。具體來看，張家口、承德、秦皇島等作為京津冀水源涵養(yǎng)功能區(qū)，應(yīng)當(dāng)依托自身地理位置和豐富的資源優(yōu)勢，推進(jìn)新能源尤其是可再生能源的開發(fā)和利用，構(gòu)建安全清潔高效的能源供應(yīng)和消費體系。打造綠色低碳循環(huán)型的經(jīng)濟(jì)發(fā)展方式，實現(xiàn)經(jīng)濟(jì)集聚區(qū)內(nèi)產(chǎn)業(yè)間的生態(tài)鏈接，并結(jié)合當(dāng)前實現(xiàn)京津冀協(xié)同發(fā)展規(guī)劃的政策優(yōu)勢，建設(shè)國家生態(tài)文明先行示范區(qū)，拓展生態(tài)內(nèi)涵，延展生態(tài)效益，走綠色低碳的發(fā)展道路。此外，河北其他城市也應(yīng)進(jìn)一步減少煤炭消費，提高能源利用效率，嚴(yán)格開展節(jié)能減排工作，淘汰落后產(chǎn)能，改善當(dāng)前資源型傳統(tǒng)產(chǎn)業(yè)占據(jù)主導(dǎo)地位，高耗能行業(yè)主導(dǎo)工業(yè)能源消費的格局。

第三，京津冀地區(qū)人口碳排放壓力的空間相關(guān)性逐漸增強(qiáng)，當(dāng)前呈現(xiàn)出顯著的高高、低低集聚態(tài)勢。碳排放存在空間溢出效應(yīng)，人口也具有集聚效應(yīng)，當(dāng)前人口碳排放壓力的空間集聚狀況要求加強(qiáng)京津冀地區(qū)尤其是人口碳排放壓力較大地區(qū)的內(nèi)部協(xié)作，建立區(qū)域減排責(zé)任分?jǐn)倷C(jī)制和減排補(bǔ)償制度，并結(jié)合當(dāng)前京津冀亟待解決的大城市病包括人口規(guī)模分布不合理等問題 ,通過區(qū)域合作與城市功能布局優(yōu)化和人口空間格局優(yōu)化,實現(xiàn)京津冀不同等級規(guī)模城市的分工合作與資源優(yōu)化組合,并積極應(yīng)用低碳技術(shù),對京津冀內(nèi)部人口碳排放壓力的溢出熱點進(jìn)行控制，強(qiáng)化碳排放的發(fā)散效應(yīng),促進(jìn)京津冀一體化和低碳化發(fā)展,實現(xiàn)區(qū)域人口碳排放壓力的平衡。