• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    京津冀人口對碳排放的壓力:空間格局及變動分析

    2020-06-10 01:24:46宮倩楠王瑩瑩童玉芬
    關(guān)鍵詞:京津冀人口區(qū)域

    宮倩楠,王瑩瑩,童玉芬

    (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ū)域人口碳排放壓力的平衡。

    猜你喜歡
    京津冀人口區(qū)域
    《世界人口日》
    人口轉(zhuǎn)型為何在加速 精讀
    英語文摘(2022年4期)2022-06-05 07:45:12
    人口最少的國家
    1723 萬人,我國人口數(shù)據(jù)下滑引關(guān)注
    關(guān)于四色猜想
    分區(qū)域
    京津冀大聯(lián)合向縱深突破
    基于嚴(yán)重區(qū)域的多PCC點暫降頻次估計
    電測與儀表(2015年5期)2015-04-09 11:30:52
    京津冀一化
    養(yǎng)老“京津冀一體化”謹(jǐn)慎樂觀看
    免费看a级黄色片| 在线观看免费日韩欧美大片| 国产又爽黄色视频| 老司机靠b影院| 搡老熟女国产l中国老女人| 国产单亲对白刺激| 狠狠狠狠99中文字幕| 亚洲三区欧美一区| 老司机午夜福利在线观看视频| 久久九九热精品免费| 亚洲熟妇中文字幕五十中出 | 亚洲 欧美 日韩 在线 免费| 国产精品国产高清国产av| 欧美精品啪啪一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 在线观看免费视频日本深夜| 黄色成人免费大全| 老司机靠b影院| 久久欧美精品欧美久久欧美| 成人国产一区最新在线观看| 亚洲自拍偷在线| 亚洲精品久久成人aⅴ小说| 亚洲av美国av| 俄罗斯特黄特色一大片| 欧美在线一区亚洲| 免费不卡黄色视频| 国产伦一二天堂av在线观看| 亚洲国产精品999在线| 国产精品爽爽va在线观看网站 | 在线免费观看的www视频| 国产成人精品久久二区二区91| 最新在线观看一区二区三区| 男男h啪啪无遮挡| 色哟哟哟哟哟哟| a级片在线免费高清观看视频| 99国产精品一区二区蜜桃av| 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 丁香欧美五月| 欧美在线一区亚洲| 在线播放国产精品三级| 亚洲精品成人av观看孕妇| 男女做爰动态图高潮gif福利片 | 亚洲人成电影免费在线| 国产1区2区3区精品| 精品一区二区三区视频在线观看免费 | 欧美一区二区精品小视频在线| 免费在线观看完整版高清| 久久久久久久久久久久大奶| 最新美女视频免费是黄的| 九色亚洲精品在线播放| 欧美老熟妇乱子伦牲交| 成年人黄色毛片网站| 中文字幕最新亚洲高清| 亚洲精品国产区一区二| 国产精品九九99| 一区在线观看完整版| xxx96com| 国产欧美日韩一区二区三区在线| 无人区码免费观看不卡| 亚洲精品av麻豆狂野| 777久久人妻少妇嫩草av网站| 免费在线观看视频国产中文字幕亚洲| 宅男免费午夜| 日韩精品免费视频一区二区三区| 亚洲男人天堂网一区| 国产成人精品在线电影| 亚洲国产欧美一区二区综合| 午夜老司机福利片| 黄色a级毛片大全视频| 热99国产精品久久久久久7| 女人被躁到高潮嗷嗷叫费观| 久久国产乱子伦精品免费另类| 国产免费av片在线观看野外av| 久久中文字幕人妻熟女| a级片在线免费高清观看视频| 日韩精品青青久久久久久| 国产精品国产高清国产av| 精品国内亚洲2022精品成人| 夜夜夜夜夜久久久久| 中国美女看黄片| 女性生殖器流出的白浆| 国产不卡一卡二| 人人妻,人人澡人人爽秒播| 免费在线观看完整版高清| 国产av精品麻豆| 免费高清在线观看日韩| 天堂中文最新版在线下载| 又大又爽又粗| 高清毛片免费观看视频网站 | 欧美丝袜亚洲另类 | 啦啦啦免费观看视频1| 国产精品久久久人人做人人爽| 老汉色av国产亚洲站长工具| 18禁美女被吸乳视频| 青草久久国产| 亚洲第一欧美日韩一区二区三区| 夜夜躁狠狠躁天天躁| 欧美国产精品va在线观看不卡| 欧美日韩国产mv在线观看视频| 久久久久亚洲av毛片大全| 亚洲国产精品999在线| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久毛片微露脸| a级毛片在线看网站| 国产精品秋霞免费鲁丝片| 级片在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 一级毛片高清免费大全| 一区二区三区激情视频| 亚洲国产毛片av蜜桃av| 琪琪午夜伦伦电影理论片6080| a在线观看视频网站| 久久草成人影院| 免费观看精品视频网站| 日日摸夜夜添夜夜添小说| 老司机亚洲免费影院| 午夜两性在线视频| 亚洲全国av大片| 男女做爰动态图高潮gif福利片 | 亚洲五月天丁香| 又大又爽又粗| 香蕉国产在线看| 久久久国产成人精品二区 | 国产免费av片在线观看野外av| 久久久精品欧美日韩精品| 在线观看免费日韩欧美大片| 亚洲第一欧美日韩一区二区三区| 亚洲精品av麻豆狂野| 国产视频一区二区在线看| 啦啦啦免费观看视频1| 成人免费观看视频高清| 一边摸一边抽搐一进一小说| www国产在线视频色| 一区二区三区国产精品乱码| 精品国产超薄肉色丝袜足j| 久久精品影院6| 精品国产一区二区久久| 一边摸一边抽搐一进一出视频| 精品人妻在线不人妻| 一级,二级,三级黄色视频| av天堂久久9| 久久婷婷成人综合色麻豆| 琪琪午夜伦伦电影理论片6080| 在线av久久热| 中文字幕精品免费在线观看视频| 少妇被粗大的猛进出69影院| 青草久久国产| 9191精品国产免费久久| 亚洲黑人精品在线| 久久久久久久久免费视频了| 色在线成人网| 香蕉国产在线看| 久久久久久久久中文| 国产精品久久久av美女十八| 男人的好看免费观看在线视频 | 久久中文字幕一级| 制服人妻中文乱码| 亚洲成人久久性| 99国产综合亚洲精品| 精品午夜福利视频在线观看一区| 超碰97精品在线观看| 久久久久久大精品| 男女高潮啪啪啪动态图| av中文乱码字幕在线| 日韩av在线大香蕉| 亚洲 欧美一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产精品 国内视频| 午夜免费鲁丝| 亚洲国产毛片av蜜桃av| 十八禁人妻一区二区| 日日爽夜夜爽网站| 国产不卡一卡二| 欧美精品一区二区免费开放| 在线十欧美十亚洲十日本专区| 精品国产美女av久久久久小说| 午夜精品久久久久久毛片777| 黑人操中国人逼视频| 日韩 欧美 亚洲 中文字幕| 在线观看免费日韩欧美大片| 咕卡用的链子| 99热只有精品国产| 少妇的丰满在线观看| 又紧又爽又黄一区二区| 琪琪午夜伦伦电影理论片6080| 精品欧美一区二区三区在线| 欧美日韩亚洲高清精品| 人成视频在线观看免费观看| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 欧美日韩福利视频一区二区| 97人妻天天添夜夜摸| 天堂动漫精品| 欧美一级毛片孕妇| 久久婷婷成人综合色麻豆| 91av网站免费观看| 欧美久久黑人一区二区| 亚洲第一av免费看| 美女 人体艺术 gogo| 中文字幕精品免费在线观看视频| 国产精品九九99| 老司机福利观看| 色哟哟哟哟哟哟| 国产伦人伦偷精品视频| 精品久久久久久久久久免费视频 | 国产欧美日韩一区二区精品| 国产精品98久久久久久宅男小说| 日本黄色日本黄色录像| 久久久久精品国产欧美久久久| 91字幕亚洲| 亚洲七黄色美女视频| 两个人看的免费小视频| 丝袜在线中文字幕| 在线天堂中文资源库| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区精品视频观看| 香蕉丝袜av| 美女 人体艺术 gogo| 国产亚洲精品综合一区在线观看 | 日韩欧美三级三区| 在线观看免费视频日本深夜| 新久久久久国产一级毛片| 国产高清激情床上av| 欧美黄色片欧美黄色片| 久久人妻福利社区极品人妻图片| 亚洲人成电影免费在线| 精品第一国产精品| 久久伊人香网站| 一区在线观看完整版| 午夜精品久久久久久毛片777| 日本wwww免费看| 黑人欧美特级aaaaaa片| 麻豆成人av在线观看| 少妇的丰满在线观看| 久久精品亚洲av国产电影网| 国产成人欧美在线观看| netflix在线观看网站| 久久婷婷成人综合色麻豆| 久久香蕉国产精品| 99久久国产精品久久久| 中文字幕高清在线视频| 波多野结衣高清无吗| 99国产综合亚洲精品| 十分钟在线观看高清视频www| 欧美日韩国产mv在线观看视频| 最好的美女福利视频网| 欧美日韩乱码在线| 美女扒开内裤让男人捅视频| 亚洲精品一卡2卡三卡4卡5卡| 黄片小视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 窝窝影院91人妻| 韩国精品一区二区三区| 黑丝袜美女国产一区| 狂野欧美激情性xxxx| 欧美日韩精品网址| 国产精品免费一区二区三区在线| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 欧美在线黄色| 中文字幕人妻丝袜制服| 精品无人区乱码1区二区| 丝袜美足系列| 热re99久久精品国产66热6| 亚洲精品一二三| 一级,二级,三级黄色视频| 欧美在线一区亚洲| 一级片'在线观看视频| 婷婷丁香在线五月| 在线观看免费午夜福利视频| 久久国产精品人妻蜜桃| 欧美性长视频在线观看| 精品无人区乱码1区二区| 成人18禁在线播放| 在线观看免费视频日本深夜| 国产精品1区2区在线观看.| 99国产极品粉嫩在线观看| 精品乱码久久久久久99久播| 超色免费av| 国产熟女午夜一区二区三区| 国产激情欧美一区二区| 91成人精品电影| 日日干狠狠操夜夜爽| 亚洲男人的天堂狠狠| 午夜激情av网站| 久久精品国产清高在天天线| 母亲3免费完整高清在线观看| 日本一区二区免费在线视频| 亚洲精品av麻豆狂野| 最好的美女福利视频网| 午夜亚洲福利在线播放| 欧美丝袜亚洲另类 | 在线观看免费视频网站a站| 桃红色精品国产亚洲av| 淫妇啪啪啪对白视频| 女性被躁到高潮视频| 精品欧美一区二区三区在线| 欧美乱码精品一区二区三区| 啦啦啦免费观看视频1| 香蕉丝袜av| 欧美另类亚洲清纯唯美| 日本免费一区二区三区高清不卡 | 久久精品国产亚洲av高清一级| 精品第一国产精品| 三上悠亚av全集在线观看| 久久久久久久午夜电影 | 国产深夜福利视频在线观看| av在线天堂中文字幕 | 亚洲熟女毛片儿| 国产亚洲欧美在线一区二区| 操美女的视频在线观看| 欧美日韩国产mv在线观看视频| 国产亚洲精品第一综合不卡| 99热只有精品国产| 国产精品秋霞免费鲁丝片| 老司机深夜福利视频在线观看| 两个人免费观看高清视频| 亚洲av成人av| 日韩有码中文字幕| 精品卡一卡二卡四卡免费| 日本欧美视频一区| 69av精品久久久久久| 欧美中文综合在线视频| 亚洲五月色婷婷综合| 美女高潮到喷水免费观看| 麻豆国产av国片精品| 精品高清国产在线一区| 人人妻,人人澡人人爽秒播| 大型av网站在线播放| 国产精品自产拍在线观看55亚洲| 日韩av在线大香蕉| 亚洲精品国产色婷婷电影| 午夜福利影视在线免费观看| 男女午夜视频在线观看| 国产有黄有色有爽视频| 日本撒尿小便嘘嘘汇集6| 国产成人精品久久二区二区免费| av天堂在线播放| 成人亚洲精品一区在线观看| 桃色一区二区三区在线观看| 欧美成人午夜精品| 免费高清在线观看日韩| 欧美成人免费av一区二区三区| 亚洲欧美日韩另类电影网站| 99热只有精品国产| 在线观看免费午夜福利视频| 亚洲色图综合在线观看| 久久精品aⅴ一区二区三区四区| 久久香蕉激情| 丁香六月欧美| 国产又色又爽无遮挡免费看| 日韩高清综合在线| 亚洲av片天天在线观看| 亚洲成人精品中文字幕电影 | 美女福利国产在线| 国产xxxxx性猛交| 18禁裸乳无遮挡免费网站照片 | 岛国视频午夜一区免费看| 久久久精品国产亚洲av高清涩受| 50天的宝宝边吃奶边哭怎么回事| 久久久久亚洲av毛片大全| 欧美日韩乱码在线| 久久国产精品人妻蜜桃| 亚洲一区二区三区色噜噜 | 多毛熟女@视频| 丝袜在线中文字幕| 日本a在线网址| 成年女人毛片免费观看观看9| 中文欧美无线码| 日本三级黄在线观看| 亚洲一区中文字幕在线| 国产成人免费无遮挡视频| 自线自在国产av| 日本a在线网址| 国产亚洲欧美精品永久| 水蜜桃什么品种好| 国产精品日韩av在线免费观看 | 日日干狠狠操夜夜爽| 国产激情欧美一区二区| 午夜免费观看网址| 性少妇av在线| 久久久久久免费高清国产稀缺| 少妇裸体淫交视频免费看高清 | 国产高清videossex| av国产精品久久久久影院| 国产蜜桃级精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 激情视频va一区二区三区| 日本a在线网址| 咕卡用的链子| 久久久水蜜桃国产精品网| 亚洲国产中文字幕在线视频| 亚洲欧美激情综合另类| 香蕉久久夜色| 午夜a级毛片| 亚洲一区二区三区欧美精品| 一进一出抽搐gif免费好疼 | 青草久久国产| 90打野战视频偷拍视频| 女警被强在线播放| 亚洲一区中文字幕在线| 亚洲av日韩精品久久久久久密| 级片在线观看| 亚洲自拍偷在线| 男人的好看免费观看在线视频 | 国产精品野战在线观看 | 亚洲欧美精品综合久久99| 国产成人欧美| 90打野战视频偷拍视频| 最近最新中文字幕大全电影3 | 欧美激情高清一区二区三区| 多毛熟女@视频| 人人妻,人人澡人人爽秒播| 高清黄色对白视频在线免费看| 亚洲av片天天在线观看| 另类亚洲欧美激情| 无人区码免费观看不卡| 又紧又爽又黄一区二区| 日本免费a在线| 91老司机精品| 在线视频色国产色| 母亲3免费完整高清在线观看| 在线观看舔阴道视频| 国产免费av片在线观看野外av| 又大又爽又粗| 啦啦啦在线免费观看视频4| 亚洲精品在线观看二区| 香蕉国产在线看| 久久人妻av系列| 久久久久久人人人人人| 黑丝袜美女国产一区| 法律面前人人平等表现在哪些方面| 在线观看免费视频日本深夜| 久久精品国产亚洲av香蕉五月| e午夜精品久久久久久久| 美女国产高潮福利片在线看| 男女下面进入的视频免费午夜 | 国产一区二区激情短视频| 91大片在线观看| 精品久久久久久,| 日韩大码丰满熟妇| 日日爽夜夜爽网站| 神马国产精品三级电影在线观看 | 丝袜在线中文字幕| 五月开心婷婷网| 精品欧美一区二区三区在线| 999久久久精品免费观看国产| 亚洲熟妇中文字幕五十中出 | 极品教师在线免费播放| 亚洲狠狠婷婷综合久久图片| 99久久综合精品五月天人人| 亚洲一区中文字幕在线| 午夜福利,免费看| 国产区一区二久久| 在线观看www视频免费| 老汉色av国产亚洲站长工具| 亚洲成国产人片在线观看| 老司机在亚洲福利影院| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 亚洲精品中文字幕在线视频| 日本a在线网址| 怎么达到女性高潮| 三上悠亚av全集在线观看| 大陆偷拍与自拍| 国产av一区在线观看免费| 亚洲精品在线观看二区| 一级,二级,三级黄色视频| 国产精品野战在线观看 | 无人区码免费观看不卡| 激情视频va一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产在线精品亚洲第一网站| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 一本综合久久免费| 日韩欧美免费精品| netflix在线观看网站| 日本精品一区二区三区蜜桃| 精品久久久久久成人av| 男人舔女人的私密视频| 欧美乱妇无乱码| 精品国内亚洲2022精品成人| 人妻丰满熟妇av一区二区三区| 亚洲欧美日韩无卡精品| 一进一出抽搐动态| 激情在线观看视频在线高清| 国产成人av激情在线播放| 国产av在哪里看| 中亚洲国语对白在线视频| 一本综合久久免费| 久久久水蜜桃国产精品网| 中文字幕最新亚洲高清| 午夜免费观看网址| 夫妻午夜视频| 在线观看一区二区三区| 欧美一级毛片孕妇| 日韩成人在线观看一区二区三区| 热99re8久久精品国产| 成人av一区二区三区在线看| 18禁裸乳无遮挡免费网站照片 | 少妇被粗大的猛进出69影院| 最新美女视频免费是黄的| 男女下面进入的视频免费午夜 | 成人三级黄色视频| 在线观看免费午夜福利视频| 亚洲精品中文字幕在线视频| 久久天躁狠狠躁夜夜2o2o| 久久天躁狠狠躁夜夜2o2o| 91精品国产国语对白视频| 一进一出好大好爽视频| 正在播放国产对白刺激| 欧美黑人精品巨大| 国产亚洲欧美精品永久| avwww免费| 欧美黑人精品巨大| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久亚洲av毛片大全| 国产精华一区二区三区| 色老头精品视频在线观看| 日韩欧美免费精品| 午夜日韩欧美国产| 亚洲精品国产精品久久久不卡| 怎么达到女性高潮| 18禁国产床啪视频网站| 色尼玛亚洲综合影院| 亚洲激情在线av| 久久午夜综合久久蜜桃| 久久久久久免费高清国产稀缺| 好男人电影高清在线观看| 国产精品偷伦视频观看了| www日本在线高清视频| 免费女性裸体啪啪无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品99久久99久久久不卡| 中国美女看黄片| 黑人猛操日本美女一级片| 窝窝影院91人妻| 一二三四社区在线视频社区8| 久久午夜综合久久蜜桃| 国产精品香港三级国产av潘金莲| 亚洲成国产人片在线观看| av天堂在线播放| 亚洲成a人片在线一区二区| 美女大奶头视频| 国产成人欧美在线观看| 国产又爽黄色视频| 老司机福利观看| 久久天堂一区二区三区四区| 一级毛片女人18水好多| 亚洲精品粉嫩美女一区| 国产1区2区3区精品| 19禁男女啪啪无遮挡网站| 99国产精品99久久久久| 成人亚洲精品一区在线观看| 日韩免费高清中文字幕av| 久久99一区二区三区| 亚洲一区中文字幕在线| 国产精品 国内视频| 欧美成人午夜精品| 久9热在线精品视频| av天堂在线播放| 99久久综合精品五月天人人| av在线播放免费不卡| 久久久久国产一级毛片高清牌| 午夜a级毛片| 国内久久婷婷六月综合欲色啪| 免费av毛片视频| xxxhd国产人妻xxx| 黑丝袜美女国产一区| 国产免费男女视频| av超薄肉色丝袜交足视频| 欧美日韩瑟瑟在线播放| 51午夜福利影视在线观看| 亚洲一区中文字幕在线| www.www免费av| 999久久久精品免费观看国产| 欧美乱妇无乱码| 国产深夜福利视频在线观看| 性欧美人与动物交配| 久久狼人影院| 欧洲精品卡2卡3卡4卡5卡区| 两个人免费观看高清视频| 中文字幕最新亚洲高清| 一二三四社区在线视频社区8| 在线观看免费高清a一片| 亚洲九九香蕉| 老司机深夜福利视频在线观看| 高清黄色对白视频在线免费看| 老熟妇乱子伦视频在线观看| a在线观看视频网站| 12—13女人毛片做爰片一| 女人精品久久久久毛片| 久久亚洲精品不卡| 最近最新中文字幕大全电影3 | 亚洲五月天丁香| 亚洲精品在线美女| 精品国产一区二区久久| 国产xxxxx性猛交| 欧美日韩亚洲国产一区二区在线观看| 丝袜在线中文字幕| 国产精品成人在线| 久久久久久人人人人人| 99热国产这里只有精品6| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清国产精品国产三级| 欧美日韩中文字幕国产精品一区二区三区 | av免费在线观看网站| 丝袜美腿诱惑在线| 又黄又爽又免费观看的视频| 如日韩欧美国产精品一区二区三区| 大型黄色视频在线免费观看| 中出人妻视频一区二区| 亚洲第一av免费看|