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

    呼包鄂榆資源型城市群“水-能-碳”耦合變化及影響因素

    2023-11-29 01:41:44張園園
    中國環(huán)境科學(xué) 2023年11期
    關(guān)鍵詞:足跡城市群耦合

    楊 屹,張園園

    呼包鄂榆資源型城市群“水-能-碳”耦合變化及影響因素

    楊 屹*,張園園

    (西安理工大學(xué)經(jīng)濟與管理學(xué)院,陜西 西安 710054)

    在描述呼包鄂榆城市群“水-能-碳”時空演變特征基礎(chǔ)上分析了耦合協(xié)調(diào)發(fā)展水平,并建立灰色關(guān)聯(lián)模型識別了耦合協(xié)調(diào)變化的關(guān)鍵影響因素.結(jié)果顯示:2005~2019年,城市群水資源消耗、能源消費和碳排放呈波動增長態(tài)勢,“水-能-碳”之間耦合協(xié)調(diào)度較低,城鎮(zhèn)化率、人口規(guī)模、經(jīng)濟發(fā)展水平是影響耦合協(xié)調(diào)變化的關(guān)鍵因素,且各城市“水-能-碳”耦合協(xié)調(diào)度的關(guān)鍵影響因素呈現(xiàn)明顯的空間異質(zhì)性.研究將為揭示黃河流域資源型城市生態(tài)系統(tǒng)動態(tài)變化過程與驅(qū)動機制、制定差異化環(huán)境治理政策提供科學(xué)依據(jù).

    黃河流域;資源型城市群;“水-能-碳”;耦合協(xié)調(diào);影響因素

    水資源短缺、能源風險和碳排放已被列入本世紀全球面臨的十大風險[1],與之對應(yīng)的清潔水(SDG6)、清潔能源(SDG7)和氣候(SDG13)成為聯(lián)合國可持續(xù)發(fā)展目標的主要關(guān)注點[2].流域作為一個包含水、能和碳等多個關(guān)聯(lián)子系統(tǒng)的復(fù)雜生態(tài)系統(tǒng),水資源和能源的投入推動了社會經(jīng)濟發(fā)展,但同時造成水資源短缺、能源危機和高碳排放等問題.黃河流域作為我國主要的資源能源聚集區(qū)域和生態(tài)功能保障區(qū)域[3-4],《黃河流域生態(tài)保護和高質(zhì)量發(fā)展規(guī)劃綱要》強調(diào)了流域水資源生態(tài)保護的重要性,而“雙碳”目標的提出給流域低碳轉(zhuǎn)型與水資源生態(tài)保護帶來了嚴峻挑戰(zhàn).因為采用碳捕獲和碳封存等綠色低碳技術(shù)或提升生態(tài)系統(tǒng)固碳能力等碳減排行動會增加用水量,采用節(jié)水技術(shù)也將導(dǎo)致額外的能源消耗和碳排放[5].水、能和碳作為主要資源環(huán)境要素,單獨地尋求某一資源要素的優(yōu)化可能會誘發(fā)更深層次的環(huán)境沖突[6].呼包鄂榆城市群作為黃河流域資源型城市群,依托豐富的礦產(chǎn)資源發(fā)展能源重化工產(chǎn)業(yè)獲得了經(jīng)濟的迅速增長,但也導(dǎo)致“水-能-碳”之間的約束性特征日益加劇,嚴重制約了流域?qū)崿F(xiàn)高質(zhì)量發(fā)展.耦合協(xié)調(diào)關(guān)系能夠有效分析不同子系統(tǒng)之間的相互權(quán)衡和協(xié)同作用,被廣泛應(yīng)用于社會經(jīng)濟和生態(tài)環(huán)境保護協(xié)調(diào)發(fā)展的政策制定中[7].因此,通過測算水資源消耗、能源消費和碳排放之間的耦合協(xié)調(diào)度來評估黃河流域資源型城市群資源環(huán)境狀況,并進一步識別關(guān)鍵影響因素是探索新時期黃河流域資源型城市群高質(zhì)量發(fā)展路徑的關(guān)鍵.

    為了應(yīng)對水資源短缺、能源安全和氣候變化問題,許多學(xué)者開始關(guān)注水、能和碳等單一環(huán)境要素下的區(qū)域資源環(huán)境狀況[8-10],但單一要素的研究難以解決多要素之間的約束性問題[11].因此,在環(huán)境治理中需要協(xié)調(diào)不同環(huán)境要素之間的關(guān)系[12].“水-能”之間的關(guān)系最早被關(guān)注[13],能源開發(fā)、生產(chǎn)和運輸?shù)入x不來水資源消耗,而取水、用水、排水和污水處理都需要能源,考慮水和能源的系統(tǒng)關(guān)系能夠以較低的成本促進城市資源可持續(xù)利用[14-15].“能-碳”和“水-碳”關(guān)系也越來越多的被關(guān)注,能源消費是影響碳排放增長的主要原因,且能源消費結(jié)構(gòu)的變化能夠遏制碳排放[16].水資源和碳排放通過與生產(chǎn)活動相關(guān)的社會經(jīng)濟網(wǎng)絡(luò)交織在一起[17],在環(huán)境治理過程中要考慮水資源和碳減排的管理政策在社會經(jīng)濟系統(tǒng)中的耦合機制[18].隨著研究的深入,逐漸拓展到對“水-能-碳”關(guān)系的探索.如,Venkatesh等[19]提出水務(wù)管理部門要考慮能源強度和溫室氣體排放;Zhang等[5]指出中國電力低碳發(fā)展和水資源消耗之間存在“空間錯配”的問題.“水-能-碳”為核心的資源環(huán)境問題也逐漸引起了學(xué)者們對影響資源環(huán)境變化的關(guān)鍵驅(qū)動因素的關(guān)注[20],旨在通過識別關(guān)鍵的影響因素來有效改善資源環(huán)境狀況.但是現(xiàn)有研究主要集中在對水、能和碳等單一環(huán)境要素影響因素的探索,研究發(fā)現(xiàn)人口規(guī)模、城鎮(zhèn)化水平是影響水污染排放的主要因素[20],經(jīng)濟增長是驅(qū)動碳排放增長的關(guān)鍵因素[21],產(chǎn)業(yè)結(jié)構(gòu)、人口規(guī)模、對外貿(mào)易等對能源效率呈現(xiàn)積極推動作用[22].

    總體而言,在從水、能和碳等單要素分析演進到“水-能”,“水-碳”和“能-碳”兩要素關(guān)系研究的過程中,以“水-能-碳”為核心的多要素關(guān)系研究才剛剛起步,仍缺乏對驅(qū)動資源環(huán)境要素耦合協(xié)調(diào)度變化因素的深入討論.此外,資源型城市的持續(xù)性較差,資源環(huán)境問題更為突出[23],促進其擺脫資源依賴困境,實現(xiàn)城市轉(zhuǎn)型是高質(zhì)量發(fā)展的關(guān)鍵[24].鑒于此,本研究創(chuàng)新之處在于,不僅結(jié)合流域資源型城市群的特點,以呼包鄂榆城市群為研究對象,從“水-能-碳”三維環(huán)境要素探討水資源利用、能源消耗和碳排放的空間分布、空間差異及動態(tài)演變特征.而且通過建立耦合協(xié)調(diào)度模型分析水、能和碳資源環(huán)境要素之間的協(xié)調(diào)關(guān)系,并進一步挖掘了耦合協(xié)調(diào)度的關(guān)鍵驅(qū)動因素,旨在為流域資源型城市群制定差異性城市環(huán)境治理政策提供依據(jù).

    1 研究方法與數(shù)據(jù)來源

    1.1 研究區(qū)域

    呼包鄂榆城市群處于黃河流域上游(圖1),包括內(nèi)蒙古呼和浩特市、包頭市、鄂爾多斯市和陜西省榆林市,在推進西部大開發(fā)和新型城鎮(zhèn)化中具有重要地位[25],是國內(nèi)唯一的國家級能源重化工基地[26]. 2020年,城市群生產(chǎn)總值13211.38億元,呼包鄂“金三角”GDP為9121.72億元,占內(nèi)蒙古的52.54%;榆林市GDP為4089.66億元,占陜西省的15.6%.雖然城市群經(jīng)濟增長較好,但是城市群地處西部黃土高原,一方面作為典型的缺水地區(qū),水土流失嚴重,生態(tài)環(huán)境脆弱;另一方面作為重要的能源和資源供應(yīng)地區(qū)[27],煤炭資源占全國的1/6,為全國提供40%的能源供給[28],導(dǎo)致碳排放量持續(xù)上升.如何緩解城市群社會經(jīng)濟發(fā)展和資源環(huán)境保護之間的矛盾,實現(xiàn)水資源、能源安全和碳減排的協(xié)同目標,是推動呼包鄂榆城市群可持續(xù)發(fā)展的關(guān)鍵.

    圖1 研究區(qū)域

    1.2 研究方法

    1.2.1 水足跡、能源消費和碳足跡核算 (1)水足跡

    水足跡是在一定的物質(zhì)條件生活標準下生產(chǎn)出一定居民消費的產(chǎn)品和服務(wù)所必需的水資源量[29].根據(jù)水足跡網(wǎng)絡(luò)制定的水足跡評價國際標準[30],參照Cai等[31]的研究成果,水足跡的計算公式如下:

    式中:為總水足跡;WF為農(nóng)業(yè)水足跡;WF為工業(yè)水足跡;WF為居民日常生活水足跡;WF為生態(tài)環(huán)境水足跡;WFWF分別為進口和出口水足跡;WF為灰水足跡.采用各地區(qū)的總水足跡與地區(qū)生產(chǎn)總值的比值代表水足跡強度.

    將水行星邊界和水足跡結(jié)合可以建立淡水資源占用是否可持續(xù)的基準[32].基于公平配額原則采用人均平等分配法[33]計算縮小后的水行星邊界,計算公式如下:

    式中:WPB為縮小后城市的水行星邊界;POP為城市的人口數(shù),POP為世界總?cè)丝跀?shù),WPB為全球水行星邊界.設(shè)置為40000億m[34].當WF大于WPB時,區(qū)域取水量超過安全范圍,產(chǎn)生水赤字.反之,產(chǎn)生盈余,表明水資源可持續(xù)利用水平越高[35].

    (2)能源消費

    能源消費直接采用各市統(tǒng)計年鑒里面公布的能源消費數(shù)據(jù).能源足跡采用碳吸收法核算吸收能源消耗產(chǎn)生的碳排放所需要的生產(chǎn)性土地面積[36],雖然也能反映能源消費,但同本研究采用的碳足跡計算中的能源消耗賬戶有所重合,可能會造成結(jié)果分析誤差.

    (3)碳足跡

    聯(lián)合國政府間氣候變化專門委員會(IPCC)制定了溫室氣體排放的核算方法[37]為計算碳足跡提供了科學(xué)參考.因此采取IPCC排放清單法計算碳足跡,計算公式如下:

    式中:CFE是能源消耗賬戶,主要測算了包括原煤、焦炭、原油、汽油等固定能源消耗產(chǎn)生的碳足跡,以及私人汽車、摩托車、載客和載貨汽車等移動能源消費產(chǎn)生的碳足跡.CFP是工業(yè)生產(chǎn)過程賬戶,主要計算硅酸鹽水泥熟料生產(chǎn)過程中產(chǎn)生的碳足跡. CFW是污染排放賬戶,主要計算了固體廢物、廢水等污染處理產(chǎn)生的碳足跡.CFL是畜禽賬戶,核算了牛、馬、豬等畜禽產(chǎn)生的碳足跡.此外,采用區(qū)域內(nèi)林地、草地等各種綠色植被的固碳能力計算碳承載.當碳足跡大于碳承載力時,區(qū)域為碳赤字,反之,為碳盈余.采用地區(qū)的碳足跡與國內(nèi)生產(chǎn)總值的比值代表碳足跡強度.

    1.2.2 耦合協(xié)調(diào)度模型 耦合協(xié)調(diào)度模型被廣泛用于探討資源、能源和環(huán)境之間的相互作用分析[38].因此,采用耦合協(xié)調(diào)度模型計算水資源消耗、能源消費和碳排放系統(tǒng)之間的耦合協(xié)調(diào)水平.計算步驟如下:

    第一,對指標進行標準化處理.計算公式如下:

    第二,測算綜合指數(shù).水足跡能有效衡量一定時期內(nèi)區(qū)域消耗的水量[39].碳足跡已經(jīng)被廣泛用于地區(qū)碳排放總量測度[40].考慮到能源足跡會與碳足跡部分賬戶存在重疊,選取能源消費總量表征能源消耗水平.因此,對水足跡、能源消費和碳足跡進行標準化的指標即為各個系統(tǒng)的綜合評價指數(shù).

    第三,計算耦合度.耦合度能夠判斷水、能和碳等子系統(tǒng)之間的相互影響程度[41].計算公式如下:

    第四,計算多系統(tǒng)耦合協(xié)調(diào)度.耦合協(xié)調(diào)度模型能更好反映兩系統(tǒng)和三系統(tǒng)耦合協(xié)調(diào)關(guān)系[42],有利于進一步分析城市群“水-能”、“水-碳”、“能-碳”和“水-能-碳”的協(xié)調(diào)發(fā)展水平.計算公式如下:

    第五,劃分協(xié)調(diào)等級.參照Li等[43]的研究成果,將耦合協(xié)調(diào)度劃分為5個等級(表1).

    1.2.3 灰色關(guān)聯(lián)模型 (1)影響因素選取

    根據(jù)現(xiàn)有研究[44-45],選取可能影響“水-能-碳”耦合協(xié)調(diào)程度的經(jīng)濟和社會發(fā)展指標作為驅(qū)動因素.其中,用年末常住人口數(shù)表征常住人口規(guī)模,城鎮(zhèn)人口比例以及第二三產(chǎn)業(yè)就業(yè)人口比例表征城鎮(zhèn)化水平,地區(qū)人均GDP表征經(jīng)濟發(fā)展水平,單位面積非農(nóng)業(yè)產(chǎn)出表征經(jīng)濟集聚程度,全社會固定資產(chǎn)投資額表征地區(qū)資本投資規(guī)模,建成區(qū)面積占城市土地總面積比例以及人均道路鋪設(shè)面積表征城市擴張.

    表1 耦合協(xié)調(diào)等級的劃分

    (2)灰色關(guān)聯(lián)模型的構(gòu)建

    灰色關(guān)聯(lián)度是根據(jù)參考數(shù)列與比較數(shù)列曲線幾何形狀接近程度來判斷不同因素的差異性對現(xiàn)象變化影響程度的相對大小[46],對樣本數(shù)量無限制,且不要求樣本服從任何概率分布[47].因此,參考劉宇峰等[44]的研究采用灰色關(guān)聯(lián)度模型分析水、能和碳耦合協(xié)調(diào)變化的關(guān)鍵影響因素.計算步驟如下:

    第二,采用公式(4)對各指標進行標準化處理.

    第三,計算指標灰色關(guān)聯(lián)度:

    1.3 數(shù)據(jù)來源

    水足跡、能源消耗、碳足跡、城市發(fā)展指數(shù)以及城鎮(zhèn)化各指標計算相關(guān)數(shù)據(jù)主要來自于2006~ 2020年《內(nèi)蒙古自治區(qū)統(tǒng)計年鑒》《陜西省統(tǒng)計年鑒》、呼包鄂榆城市群各市的統(tǒng)計年鑒、國民經(jīng)濟與社會發(fā)展統(tǒng)計公報以及政府網(wǎng)站、中國城市統(tǒng)計年鑒和中國能源統(tǒng)計年鑒.水足跡中的工業(yè)、居民生活和生態(tài)用水的數(shù)據(jù)來源于陜西省和內(nèi)蒙古自治區(qū)的水資源公報.

    2 結(jié)果分析

    2.1 “水-能-碳”時空演變特征分析

    2.1.1 水足跡時空演變特征分析 2005~2019年,呼包鄂榆城市群水足跡總量呈“先上升后下降再上升”的變化態(tài)勢,由2005年的162.47億m3增至2012年的222.85億m3,接著降至2016年的196.23億m3后又增至2019年的223.74億m3,總體增長了37.71%(圖2).2005~2012年,水足跡的“上升”主要是由于人口和經(jīng)濟規(guī)模的擴大,該時間段常住人口增長了105.98萬人,增長率為10.6%,經(jīng)濟快速發(fā)展和人口規(guī)模迅速擴大會導(dǎo)致居民生活、農(nóng)業(yè)和工業(yè)生產(chǎn)用水的增加[48],居民生活、農(nóng)業(yè)、工業(yè)、灰水足跡分別增加了37%、50%、33%和12%.2012~2016年,城市群水足跡的“下降”主要體現(xiàn)在呼和浩特市和包頭市水足跡的快速下降,農(nóng)業(yè)、工業(yè)、灰水足跡分別下降了20.8%、22.5%、80.6%和22.1%、12.0%、94.98%,在節(jié)水、治水方面取得了良好成效.2016~ 2019年,呼包鄂榆城市群水足跡的上升主要體現(xiàn)在榆林市和鄂爾多斯市水足跡的上升.尤其是榆林市2018~2019年水足跡總量超過呼和浩特市,成為城市群水赤字最為嚴重的城市.主要原因是榆林不僅是能源化工強市,并且在“十三五”期間積極推進特色農(nóng)業(yè)建設(shè),導(dǎo)致該時段內(nèi)農(nóng)業(yè)和工業(yè)用水量的快速增加,分別增加了68.8%和49.7%.值得注意的是, 2005~2019年,在呼包鄂榆城市群的水足跡強度由700.84m3/萬元降至168.89m3/萬元,下降了75.9%(圖2),說明水資源經(jīng)濟效率在不斷提升.其中,呼和浩特市和包頭市的水足跡強度雖然均高于城市群,但一直呈下降態(tài)勢.而鄂爾多斯市和榆林市的水足跡強度均呈“先下降后上升”的態(tài)勢,說明兩市水資源消費經(jīng)濟效率近幾年有所降低.提高水資源利用效率是當前呼包鄂榆城市群急需解決的問題.從各賬戶來看(圖3a),呼包鄂榆城市群農(nóng)業(yè)水足跡占比最大,但逐漸降低,其次是工業(yè)水足跡,且占比逐年升高.生態(tài)水足跡和灰水足跡一增一減,表明生態(tài)建設(shè)和環(huán)境保護得到了越來越多的重視.城市群一直呈虛擬水凈出口加劇了資源型城市水資源壓力,但占比呈降低態(tài)勢.此外,水行星邊界計算結(jié)果顯示(圖3b),呼包鄂榆城市群2005~2019年水行星邊界變化不大,且均小于水足跡總量,從而產(chǎn)生了嚴重的水赤字,表明流域資源型城市水資源利用在絕對可持續(xù)評估中也遠遠超過了水的安全使用范圍,阻礙SDG6目標的實現(xiàn).

    圖2 2005~2019年呼包鄂榆城市群水足跡和水足跡強度時空演變

    圖3 2005~2019年呼包鄂榆城市群水足跡各賬戶(a)和水足跡、水行星邊界與水赤字動態(tài)變化(b)

    2.1.2 能源消耗時空演變特征分析 圖4顯示,雖然2005~2019年呼包鄂榆城市群及各市能源消耗總量呈波動增長的態(tài)勢,但是增長速度變化明顯. 2005~2012年,城市群能源消耗高速增長,由5982.49萬tce增至11777.08萬tce,年平均增長率為12.11%, 2012年之后能源消耗增速降低,由11777.08萬tce增至2019年的15356.25萬tce,年平均增長率為4.34%.2005~2012年呼包鄂榆城市群依靠其豐富的煤炭資源實現(xiàn)了經(jīng)濟的快速增長,但是造成了資源能源的過度消耗.2012年之后受煤炭價格回落影響,經(jīng)濟增速放緩,使得資源能源消費增速下降.包頭市是城市群中年平均能源消耗量最多的城市,但增速較慢,2005~2019年由2410.86萬tce增至4658.86萬tce,增長了93.24%,其次是鄂爾多斯市.榆林市的能源消費量增長幅度最大,從2005年的634.98萬tce增至2019年的4327.41萬tce,增長了581.50%,榆林市是國家級能源化工基地,近些年憑借其富集的能源礦產(chǎn)資源實現(xiàn)了經(jīng)濟快速發(fā)展.呼和浩特市是城市群中年平均能源消費最低的城市.2005~2019年,呼包鄂榆城市群能源消費強度整體呈下降趨勢,由2005年的2.35tce/萬元降至2019年的0.86tce /萬元.其中,呼和浩特市能源消費強度最低,平均能源消費強度為1.13tce/萬元,這是因為該市2020年三次產(chǎn)業(yè)比例為4.5:29.1:66.4,主要以旅游、金融、商貿(mào)等服務(wù)業(yè)帶動經(jīng)濟增長.包頭市、鄂爾多斯市和榆林市經(jīng)濟發(fā)展以資源能源相關(guān)產(chǎn)業(yè)為主,能源消費強度相對較高,其中,包頭市能源消耗強度最高,但是其近十五年的降幅也最大.總體來看,2005~2012年間能源消費強度在逐年下降,但是整體尚未達到全國單位GDP能耗水平,能源利用率需進一步提升.

    圖4 2005~2019年呼包鄂榆城市群能源消費和能源消費強度時空演變

    2.1.3 碳足跡時空演變特征分析 2005~2019年,呼包鄂榆城市群碳足跡呈“先降后升”的變化態(tài)勢(圖5).其中,2005~2007年由12.17億tCO2eq降至6.42億tCO2eq,下降了47.25%,主要是晉陜蒙寧交界區(qū)域的環(huán)境整治行動一定程度上減少了碳排放.2008年之后,個別年份出現(xiàn)小幅的下降,但整體呈波動上升態(tài)勢,2019年增至16.51億tCO2eq,增長了157.17%,主要與能源消費波動有關(guān).呼包鄂榆城市群碳足跡賬戶中,能源消耗賬戶占比最大,除2007和2008年外,其他年份該賬戶碳足跡均占碳足跡總量的90%以上,其變化態(tài)勢與城市群碳足跡總量變化趨勢一致,說明城市群碳排放主要來源于能源消耗.城市群碳承載力雖然在2005~2019略微增長,但始終遠遠低于碳足跡呈現(xiàn)碳赤字.其中,鄂爾多斯市碳足跡最高,是城市群最大的碳排放源區(qū), 2005~2019年平均碳足跡為5.44億tCO2eq,城市群平均碳足跡為10.5億tCO2eq,鄂爾多斯市就導(dǎo)致了城市群一半的排放量,鄂爾多斯市的發(fā)展高度依賴稀土高嶺土資源,以及豐富的煤炭天然氣,以工業(yè)為主的發(fā)展模式導(dǎo)致碳排放較高.2005~2019年,榆林的平均碳足跡為2.13億tCO2eq,占城市群平均碳足跡的20.25%,是碳足跡增長最快的城市,從2005年的0.43億tCO2eq增至2019年的4.38億tCO2eq,增長了918.6%,榆林的產(chǎn)業(yè)結(jié)構(gòu)存在嚴重失衡問題,2020年三產(chǎn)業(yè)占比分別為6.7:62.5:30.8,是典型的煤炭資源型城市.呼和浩特市和包頭市2005~2019年的平均碳足跡相對較小,而且呼和浩特市是城市群唯一碳足跡降低的城市,下降了51.7%.2005~2019年呼包鄂榆城市群碳足跡強度呈“先下降后上升”態(tài)勢,由2005年的52.51tCO2eq/萬元降至2016年的9.87tCO2eq/萬元,但2017年之后又出現(xiàn)增長,增至2019年的12.46tCO2eq/萬元.

    圖5 2005~2019年呼包鄂榆城市群碳足跡和碳足跡強度時空演變

    2.2 “水-能-碳”耦合協(xié)調(diào)關(guān)系分析

    為了進一步分析“水-能-碳”為核心的資源環(huán)境問題,本研究分別構(gòu)建“水-能”、“水-碳”、“能-碳”和“水-能-碳”的多維耦合協(xié)調(diào)模型.圖6結(jié)果顯示,從城市群整體來看,2005~2019年呼包鄂榆城市群各系統(tǒng)之間平均耦合協(xié)調(diào)度均沒有達到高度協(xié)調(diào),且耦合協(xié)調(diào)水平的變化較大.“水-能-碳”三元耦合協(xié)調(diào)水平呈“減-增-減”的弱倒“N”型變化態(tài)勢,從2005年0.9降至2012年的0.39,再增至2016年的0.59,后又降至2019年的0.01,整體波動下降,表現(xiàn)為“高度協(xié)調(diào)-中度失調(diào)-中度協(xié)調(diào)-嚴重失調(diào)”.變化主要分以下三個階段:第一個階段是2005~ 2012年,正處于煤炭經(jīng)濟的黃金時期,經(jīng)濟發(fā)展嚴重依賴于煤炭資源的開采和加工,導(dǎo)致水資源消耗、能源消費和碳排放的大量增加,“水-能-碳”耦合協(xié)調(diào)水平大幅下降.第二階段是2013~2016年,在此期間煤炭需求下降導(dǎo)致全國煤炭市場供大于求,全國煤炭產(chǎn)量也在快速減少.此外,《全國資源型城市可持續(xù)發(fā)展規(guī)劃(2013~2020年)》的發(fā)布要求加快資源型城市轉(zhuǎn)變經(jīng)濟發(fā)展方式.在此背景下,呼包鄂榆城市群水資源消耗降低,能源消費增速減慢,碳排放也在2013年和2014年有一定程度上的改善,“水-能-碳”的耦合協(xié)調(diào)水平有所提升,說明政策的實施一定程度上促進了資源環(huán)境耦合協(xié)調(diào)發(fā)展.第三階段是2017~2019年,2017年全國原煤產(chǎn)量自2014年以來首次出現(xiàn)恢復(fù)性增長,全國整體煤炭產(chǎn)量波動增加,呼包鄂榆城市群作為中國典型的資源型城市群,能源和水資源消耗也快速提升,并由此導(dǎo)致碳排放量增加,資源環(huán)境約束加劇,“水-能-碳”耦合協(xié)調(diào)度下降.

    進一步研究結(jié)果發(fā)現(xiàn),“水-能”、“水-碳”和“能-碳”耦合協(xié)調(diào)度變化趨勢存在明顯的差異.“水-能”和“水-碳”的耦合協(xié)調(diào)度均呈“減-增-減”變化態(tài)勢,主要是受水資源消耗量變動的影響,2013~2016年水資源消耗降低,促使耦合協(xié)調(diào)度在該時間段內(nèi)有所提升,但是2005~2019年整體呈倒“N”型的波動降低變化態(tài)勢,這與王風初等[15]、李波等[49]的研究結(jié)果存在差異.主要原因是資源型城市群的資源依賴發(fā)展模式對以第二產(chǎn)業(yè)為主的產(chǎn)業(yè)結(jié)構(gòu)具有鎖定效應(yīng),產(chǎn)業(yè)轉(zhuǎn)型困難,導(dǎo)致資源環(huán)境的耦合協(xié)調(diào)度依然較低.“能-碳”耦合協(xié)調(diào)水平呈波動下降態(tài)勢,從2005年的1.00降至2019年的0.01,表現(xiàn)為“高度協(xié)調(diào)-嚴重失調(diào)”,主要是因為能源和碳排放之間關(guān)系較為密切,能源消費的增加加劇了碳排放.與趙金輝等[50]的研究結(jié)果一致,認為城市群粗放的發(fā)展模式導(dǎo)致能源消耗和以碳排放為主的環(huán)境污染嚴重,進而使兩者之間的耦合協(xié)調(diào)水平下降.從呼包鄂榆城市群各市情況來看,呼和浩特市的資源環(huán)境要素之間的耦合協(xié)調(diào)度雖然波動較大,但2012年開始耦合協(xié)調(diào)水平不斷提升,尤其是2016年之后各要素之間均呈“高度協(xié)調(diào)”狀態(tài).這主要是因為2012年呼和浩特市開始推進大青山和敕勒川草原生態(tài)綜合治理,環(huán)境治理效果明顯.包頭市資源環(huán)境要素之間的耦合協(xié)調(diào)度變化態(tài)勢與城市群整體基本一致,除了“能-碳”耦合協(xié)調(diào)度持續(xù)下降外,其余均呈“減-增-減”變化態(tài)勢.主要是由于包頭市2014年開始積極推行大氣污染防治、城市供水改造和生態(tài)環(huán)境整治行動,從而有效的改善了資源環(huán)境問題.鄂爾多斯市和榆林市水、能和碳之間的耦合協(xié)調(diào)度均呈波動下降態(tài)勢,說明這兩市的資源環(huán)境在不斷惡化,嚴重阻礙了流域資源型城市群實現(xiàn)“雙碳”目標和流域高質(zhì)量發(fā)展.

    2.3 影響因素分析

    城市群整體分析結(jié)果顯示(表2),本研究選取的社會和經(jīng)濟發(fā)展指標與水、能和碳耦合協(xié)調(diào)度的關(guān)聯(lián)度均大于0.6,說明城市群社會經(jīng)濟發(fā)展與資源環(huán)境消耗的關(guān)聯(lián)度較強.其中,城鎮(zhèn)化率是水、能和碳之間耦合協(xié)調(diào)度的關(guān)鍵影響因素,第二三產(chǎn)業(yè)就業(yè)人口占比和城鎮(zhèn)人口比例與耦合協(xié)調(diào)度的關(guān)聯(lián)度分別為0.82和0.70,說明呼包鄂榆城市群城市人口擴張和城鎮(zhèn)化水平的提升加劇了資源能源利用和污染排放.年末常住人口、地區(qū)人均GDP和單位面積非農(nóng)業(yè)產(chǎn)出與耦合協(xié)調(diào)度的關(guān)聯(lián)度均為0.69,表明常住人口規(guī)模、經(jīng)濟發(fā)展水平和經(jīng)濟集聚程度對耦合協(xié)調(diào)的影響程度次之.人口規(guī)模的擴大增加了水資源、能源消耗和污染排放,不利于水、能和碳的耦合協(xié)調(diào),而經(jīng)濟集聚過程中產(chǎn)生的規(guī)模經(jīng)濟、技術(shù)和知識溢出效應(yīng)會有利于水、能和碳的耦合協(xié)調(diào)的提升.全社會固定投資、人均道路鋪設(shè)面積和建成區(qū)面積占土地總面積比例與水、能和碳的耦合協(xié)調(diào)度的灰色關(guān)聯(lián)度分別為0.67、0.65和0.64,表明投資規(guī)模和城市擴張對水、能和碳的耦合協(xié)調(diào)度產(chǎn)生了影響.此外,研究結(jié)果顯示,“水-能”、“水-碳”、“能-碳”和“水-能-碳”耦合協(xié)調(diào)度影響最大的均是城鎮(zhèn)化水平,表明快速城鎮(zhèn)化過程中城市人口的迅速增長以及產(chǎn)業(yè)結(jié)構(gòu)的轉(zhuǎn)變帶來的資源能源消耗以及居民生活和工業(yè)污染排放是影響水、能和碳耦合協(xié)調(diào)的重要原因.

    城市群各城市分析結(jié)果顯示(表3~表6),水、能和碳耦合協(xié)調(diào)度的影響因素呈現(xiàn)明顯的空間異質(zhì)性.具體而言,呼和浩特市年末常住人口(0.75)和城鎮(zhèn)人口比例(0.74)與水、能和碳耦合協(xié)調(diào)度的關(guān)聯(lián)度最高,包頭市地區(qū)人均GDP(0.73)、單位面積非農(nóng)業(yè)產(chǎn)出(0.73)和人均道路鋪設(shè)面積(0.73)與水、能和碳耦合協(xié)調(diào)度的關(guān)聯(lián)度最強,鄂爾多斯市建成區(qū)面積占土地總面積比例(0.84)和城鎮(zhèn)人口比例(0.74)與水、能和碳耦合協(xié)調(diào)度的關(guān)聯(lián)度最高,榆林市第二三產(chǎn)業(yè)就業(yè)人口占比(0.75)和建成區(qū)面積占土地總面積比例(0.73)與水、能和碳耦合協(xié)調(diào)度的關(guān)聯(lián)度最高.影響因素的差異主要與各市的資源稟賦和產(chǎn)業(yè)結(jié)構(gòu)的差異相關(guān).呼和浩特市作為內(nèi)蒙古自治區(qū)的省會城市,是城市群的公共服務(wù)和現(xiàn)代服務(wù)中心,常住人口規(guī)模和城鎮(zhèn)化水平的提升,伴隨著收入水平的提高和生活方式改變,促進了資源能源的集約利用,降低了CO2的排放,促進了水、能和碳耦合協(xié)調(diào)度的提升.包頭市作為老工業(yè)基地,目前正在大力發(fā)展科技創(chuàng)新,向現(xiàn)代工業(yè)城市轉(zhuǎn)型,經(jīng)濟發(fā)展水平的提升往往需要消耗大量的資源能源并產(chǎn)生碳排放,對水、能和碳的耦合協(xié)調(diào)度產(chǎn)生不利的影響,而經(jīng)濟集聚程度的提升可以通過促進規(guī)模經(jīng)濟來節(jié)約資源能源消耗和降低污染排放,促進水、能和碳的耦合協(xié)調(diào)度的提升.鄂爾多斯市在近十五年間進行了快速的城市擴張,建成區(qū)面積占比增長了236.64%,城市人口增長和城市建設(shè)消耗了大量的資源和能源,并造成大量碳排放,因此城市擴張和城鎮(zhèn)化是影響水、能和碳耦合協(xié)調(diào)關(guān)鍵因素.榆林市是我國重要的煤炭和石油后備基地,社會經(jīng)濟發(fā)展嚴重依賴高碳排放產(chǎn)業(yè),第二三產(chǎn)業(yè)就業(yè)人口增加,第二三產(chǎn)業(yè)的產(chǎn)值規(guī)模相較于第一產(chǎn)業(yè)會排放更多CO2.

    表2 呼包鄂榆城市群“水-能-碳”耦合協(xié)調(diào)度的影響因素

    表3 呼和浩特市“水-能-碳”耦合協(xié)調(diào)的影響因素

    表4 包頭市“水-能-碳”耦合協(xié)調(diào)的影響因素

    表5 鄂爾多斯市“水-能-碳”耦合協(xié)調(diào)的影響因素

    表6 榆林市“水-能-碳”耦合協(xié)調(diào)的影響因素

    3 結(jié)論

    3.1 呼包鄂榆城市群以水、能和碳為核心的資源環(huán)境問題突出.2005~2019年,呼包鄂榆城市群水足跡、能源消耗和碳足跡整體上呈現(xiàn)增長態(tài)勢,且變化態(tài)勢呈現(xiàn)明顯的差異性和波動性.呼和浩特市水資源壓力最大,而包頭市和鄂爾多斯市能源消耗和碳排放最大,榆林市是城市群中“水-能-碳”資源環(huán)境壓力增長最快的城市.整體看,呼包鄂榆城市群水資源消耗和能源消費需求較大,碳排放問題嚴重,資源環(huán)境壓力主要受能源產(chǎn)業(yè)發(fā)展的影響.

    3.2 呼包鄂榆城市群水、能和碳的耦合協(xié)調(diào)度較低且呈動態(tài)復(fù)雜的變化特征.2005~2019年,城市群“水-能”、“水-碳”、“能-碳”和“水-能-碳”之間的平均耦合協(xié)調(diào)度較低,“水-能”、“水-碳”和“水-能-碳”耦合協(xié)調(diào)度均呈“減-增-減”的倒“N”型變化態(tài)勢,“能-碳”耦合協(xié)調(diào)度呈波動降低態(tài)勢.此外,城市群水資源、能源和碳排放系統(tǒng)之間的耦合協(xié)調(diào)度存在顯著的空間差異.呼和浩特市資源環(huán)境之間的耦合協(xié)調(diào)度在不斷提升,并實現(xiàn)了高度協(xié)調(diào).包頭市的耦合協(xié)調(diào)度雖然在一定時期內(nèi)有所提升,但是整體依然呈下降態(tài)勢.鄂爾多斯市和榆林市的水、能和碳之間的耦合協(xié)調(diào)度最低且不斷下降,目前均呈現(xiàn)“嚴重失調(diào)”.說明城市群水資源消耗、能源消費和碳排放之間的約束性不斷加劇.

    3.3 社會和經(jīng)濟發(fā)展對呼包鄂榆城市群各市水、能和碳耦合協(xié)調(diào)度存在顯著的影響,且影響因素存在顯著的差異性.城鎮(zhèn)化是影響水、能和碳耦合協(xié)調(diào)度的關(guān)鍵影響因素,常住人口規(guī)模、經(jīng)濟發(fā)展水平和經(jīng)濟集聚程度對耦合協(xié)調(diào)的影響程度次之,同時投資規(guī)模和資本擴張也對水、能和碳耦合協(xié)調(diào)度產(chǎn)生顯著的影響.城市群內(nèi)不同城市水、能和碳耦合協(xié)調(diào)度之間的關(guān)鍵影響因素存在顯著的差異.呼和浩特市、包頭市、鄂爾多斯市和榆林市的關(guān)鍵影響因素分別是常住人口規(guī)模、經(jīng)濟發(fā)展水平、城市擴張和城鎮(zhèn)化水平.

    4 政策建議

    流域資源型城市群往往陷入“資源詛咒”而難以擺脫“剛性”的資源依賴性[51-52],其系統(tǒng)變革與轉(zhuǎn)型發(fā)展不但是國內(nèi)的難點,在世界上也是難題.這些城市群面臨著“雙碳”目標同水資源保護與利用的雙重困境,揭示新時期流域資源型城市群“水-能-碳”關(guān)系耦合變化過程及影響機理,探討提升綠色低碳現(xiàn)代城市服務(wù)功能的路徑是體現(xiàn)新發(fā)展理念的關(guān)鍵.

    4.1 城市群應(yīng)加強對水資源、能源利用以及碳排放的管理,一方面,政府可以通過提供適當?shù)募夹g(shù)創(chuàng)新獎勵促進企業(yè)淘汰掉老舊且高耗能的設(shè)備,積極引進先進技術(shù)以及自主研發(fā)來提升技術(shù)水平和改進生產(chǎn)工藝,從而提升資源利用效率,另一方面,通過政府干預(yù)來減少污染排放,從而保證水、能和碳耦合協(xié)調(diào)發(fā)展.

    4.2 城市群應(yīng)構(gòu)建資源綜合管理平臺促進不同環(huán)境要素管理部門的協(xié)同管理,在“雙碳”目標下建立水資源、能源集約利用的約束機制,提升生活生產(chǎn)過程中的水能利用效率,制定水能系統(tǒng)碳減排路徑,加強不同部門和產(chǎn)業(yè)的協(xié)同減排,以降低用水過程的能源消費和碳排放,并降低碳減排過程中的水資源消耗.

    4.3 城市群整體在快速城鎮(zhèn)化過程中要充分考慮城市人口增加對資源環(huán)境的脅迫作用,處理好城市發(fā)展和生態(tài)保護之間的關(guān)系,推動生態(tài)宜居的新型城鎮(zhèn)化建設(shè),促進產(chǎn)業(yè)之間基于跨產(chǎn)業(yè)鏈合作.此外,城市群內(nèi)各市要根據(jù)自身獨特的資源稟賦、發(fā)展模式和功能定位針對性的制定措施.如,呼和浩特市通過優(yōu)化能源結(jié)構(gòu)、改善住房條件等來緩解城市人口壓力,包頭市應(yīng)該繼續(xù)推動科技創(chuàng)新發(fā)展和現(xiàn)代工業(yè)城市轉(zhuǎn)型,提升資源環(huán)境效率,鄂爾多斯市在城市土地空間的擴張時應(yīng)充分考慮生態(tài)環(huán)境效益,限制擴張范圍,實現(xiàn)土地資源的最優(yōu)化利用.榆林市第二三產(chǎn)業(yè)就業(yè)人口占比成為影響其資源環(huán)境耦合協(xié)調(diào)的關(guān)鍵驅(qū)動因素,應(yīng)該積極推動二三產(chǎn)業(yè)供給側(cè)改革,通過二三產(chǎn)業(yè)的技術(shù)創(chuàng)新來促進企業(yè)節(jié)能減排.

    本研究也存在一定局限性.第一,由于呼包鄂榆城市群僅包含四座城市,不滿足回歸分析大樣本的要求.因此,采用對樣本量和樣本分布情況沒有限制的灰色關(guān)聯(lián)度模型進行影響因素分析,沒有反映具體的因果關(guān)系,我們在后續(xù)研究中將擴大資源型城市的范圍、選擇黃河流域資源型城市進行研究,提升研究成果的普適性;第二,由于能源足跡的大多采用吸收能源消耗產(chǎn)生CO2所需土地面積進行計算,與碳足跡的能源消耗賬戶計算重合,因此本文直接采用能源消費總量進行測度,資源利用量測度指標有待于進一步研究.

    [1] Liang M S, Huang G H, Chen J P, et al. Energy-water-carbon nexus system planning: A case study of Yangtze River Delta urban agglomeration, China [J]. Applied Energy, 2022,308,118144.

    [2] The United Nations. 2015. Sustainable Development Goals. 2022- 03-10, http://www.un.org/sustainabledevelopment/zh/sustainable- development-goals/.

    [3] 孫久文,崔雅琪,張 皓.黃河流域城市群生態(tài)保護與經(jīng)濟發(fā)展耦合的時空格局與機制分析[J]. 自然資源學(xué)報, 2022,37(7):1673-1690. Sun J W, Cui Y Q, Zhang H. Spatio-temporal pattern and mechanism analysis of coupling between ecological protection and economic development of urban agglomerations in the Yellow River Basin [J]. Journal of Natural Resources, 2022,37(7):1673-1690.

    [4] 趙武生,石培基.基于生態(tài)系統(tǒng)服務(wù)的蘭西城市群復(fù)合生態(tài)系統(tǒng)耦合協(xié)調(diào)關(guān)系研究[J]. 中國環(huán)境科學(xué), 2022,42(11):5358-5368. Zhao W S, Shi P J. Research on the the coupling and coordination of the complex ecosystem in the Lanzhou-Xining urban agglomeration based on ecosystem services [J]. China Environmental Science, 2022, 42(11):5358-5368.

    [5] Zhang C, He G, Johnston J, et al. Long-term transition of China’s power sector under carbon neutrality target and water withdrawal constraint [J]. Journal of Cleaner Production, 2021,329,129765.

    [6] Wang X, Kleme? J, Ouyang X et al. Regional embodied Water- Energy-Carbon efficiency of China [J]. Energy, 2021,224,120159.

    [7] 饒清華,林秀珠,李家兵,等.流域社會經(jīng)濟與水環(huán)境質(zhì)量耦合協(xié)調(diào)度分析[J]. 中國環(huán)境科學(xué), 2019,39(4):1784-1792. Rao Q H, Lin X Z, Li J B, et al. Analysis of coupling coordination between social economy and water environment quality in river basin [J]. China Environmental Science, 2019,39(4):1784-1792.

    [8] 任玉芬,蘇小婉,賀玉曉,等.中國生態(tài)地理區(qū)城市水資源利用效率及影響因素[J]. 生態(tài)學(xué)報, 2020,40(18):6459-6471. Ren Y F, Su X W, He Y X, et al. Urban water resource utilization efficiency and its influencing factors in eco-geographic regions of China [J]. Acta Ecologica Sinica, 2020,40(18):6459-6471.

    [9] 杜昱東,劉 戀,王佳斌,等.基于碳排放約束的北京市社區(qū)生態(tài)效率評價[J]. 生態(tài)學(xué)報, 2023,(2):1-13. Du Y D, Liu L, Wang J B, et al. Evaluation of community’ s ecological efficiency in Beijing based on carbon emission constraint [J]. Acta Ecologica Sinica, 2023,(2):1-13.

    [10] 王雅晴,譚德明,張佳田,等.我國城市發(fā)展與能源碳排放關(guān)系的面板數(shù)據(jù)分析[J]. 生態(tài)學(xué)報, 2020,40(21):7897-7907. Wang Y Q, Tan D M, Zhang J T, et al. The impact of urbanization on carbon emissions: Analysis of panel data from 158cities in China [J]. Acta Ecologica Sinica, 2020,40(21):7897-7907.

    [11] 孫才志,魏亞瓊,趙良仕.干旱區(qū)水—能源—糧食紐帶系統(tǒng)協(xié)同演化——以中國西北地區(qū)為例[J]. 自然資源學(xué)報, 2022,37(2):320-333. Sun C Z, Wei Y Q, Zhao L S. Co-evolution of water-energy-food nexus in arid areas: Take Northwest China as an example [J]. Journal of Natural Resources, 2022,37(2):320-333.

    [12] Ahmad S, Jia H, Chen Z, et al. Water-energy nexus and energy efficiency: A systematic analysis of urban water systems [J]. Renewable and Sustainable Energy Reviews, 2020,134,110381.

    [13] Wilkinson R. Methodology for analysis of the energy intensity of California’s water systems and an assessment of multiple potential benefits through integrated water-energy efficiency measures [D]. University of California Santa Barbara. 2000.

    [14] de Oliveira G C, Bertone E, Stewart R A. Challenges, opportunities, and strategies for undertaking integrated precinct-scale energy–water system planning [J]. Renewable and Sustainable Energy Reviews, 2022,161,112297.

    [15] 王風初,曹建軍,王 寧,等.近20年我國虛擬水、能消耗及耦合和需求預(yù)測[J]. 中國環(huán)境科學(xué), 2022,42(10):4919-4930. Wang F C, Cao J J, Wang N, et al. Consumption and coupling of virtual water and virtual energy in the past 20years and predicting their demand for the next decade in China [J]. China Environmental Science, 2022,42(10):4919-4930.

    [16] 王長建,汪 菲,張虹鷗.新疆能源消費碳排放過程及其影響因素——基于擴展的Kaya恒等式[J]. 生態(tài)學(xué)報, 2016,36(8):2151-2163. Wang C J, Wang F, Zhang H O. The process of energy-related carbon emissions and influencing mechanism research in Xinjiang [J]. Acta Ecologica Sinica, 2016,36(8):2151-2163.

    [17] Wang P, Huang G, Li Y. A factorial stepwise-clustering input-output model for unveiling water-carbon nexus from multi-policy perspectives [J]. Science of The Total Environment, 2023,866:161315.

    [18] Fang D, Chen F. Linkage analysis for water-carbon nexus in China [J]. Applied Energy, 2018,225:682-695.

    [19] Venkatesh G, Chan A, Bratteb? H. Understanding the water-energy- carbon nexus in urban water utilities: Comparison of four city case studies and the relevant influencing factors [J]. Energy, 2014,75:153–166.

    [20] 周 侃,伍健雄,樊 杰,等.長江經(jīng)濟帶環(huán)境污染脅迫的驅(qū)動因素及空間效應(yīng)[J]. 環(huán)境科學(xué)學(xué)報, 2021,41(5):1996-2004. Zhou K, Wu J X, Fan J, et al. 2021. Driving factors of environmental pollution stress and their spatial effects in the Yangtze River Economic Belt [J]. Acta Scientiae Circumstantiae, 2021,41(5):1996-2004.

    [21] 王 正,樊 杰.能源消費碳排放的影響因素特征及研究展望[J]. 地理研究, 2022,41(10):2587-2599. Wang Z, Fan J. The characteristics and prospect of influencing factors of energy-related carbon emissions: Based on literature review [J]. Geographical Research, 2022,41(10):2587-2599.

    [22] 陳菁泉,連欣燕,馬曉君,等.中國全要素能源效率測算及其驅(qū)動因素[J]. 中國環(huán)境科學(xué), 2022,42(5):2453-2463. Chen J Q, Lian X Y, Ma X J, et al. Total factor energy efficiency measurement and drivers in China [J]. China Environmental Science, 2022,42(5):2453-2463.

    [23] 邊志強,鐘順昌.資源型城市經(jīng)濟增長目標對碳排放效率的影響[J]. 中國環(huán)境科學(xué), 2023,43(8):4395-4408 Bian Z Q, Zhong S C. The impact of economic growth target on carbon emission efficiency in resource-based cities [J]. China Environmental Science, 2023,43(8):4395-4408.

    [24] 吳 康,張文忠,張平宇,等.中國資源型城市的高質(zhì)量發(fā)展:困境與突破[J]. 自然資源學(xué)報, 2023,38(1):1-21. Wu K, Zhang W Z, Zhang P Y, et al. High-quality development of resource-based cities in China: Dilemmas and breakthroughs [J]. Journal of Natural Resources, 2023,38(1):1-21.

    [25] 國家發(fā)展改革委關(guān)于印發(fā)呼包鄂榆城市群發(fā)展規(guī)劃的通知 [N]. 中央人民政府網(wǎng)站,http://www.gov.cn/xinwen/2018-03/07/content_ 5271788.htm, 2018-03-07. Notice of the National Development and Reform Commission on Printing and Distributing the Development Plan of Hohhot-Baotou- Ordos-Yulin Urban Agglomeration [N]. The website of the Central People's Government, http://www.gov.cn/xinwen/2018-03/07/content_ 5271788.htm, 2018-03-07.

    [26] 崔遠政,姜 磊,張為師,等.中國西部人為源NO2污染的時空分布及影響因素研究[J]. 環(huán)境科學(xué)學(xué)報, 2019,39(12):4152-4161. Cui Y Z, Jiang L, Zhang W S, et al. Spatio-temporal variations and influential factors of anthropogenic NO2pollution over Western China [J]. Acta Scientiae Circumstantiae, 2019,39(12):4152-4161.

    [27] 康 哲,李 巍,劉 偉.黃河流域城市群工業(yè)減污降碳影響因素及協(xié)同推進策略研究[J]. 中國環(huán)境科學(xué), 2023,43(4):1946-1956. Tang Z, Li w, Liu W. Influencing factors and synergistic promotion strategy of industrial pollution and carbon reduction at urban agglomerations of the Yellow River Basin [J]. China Environmental Science, 2023,43(4):1946-1956.

    [28] 鄒亞鋒,段 偉,郭娜娜.呼包鄂榆城市群經(jīng)濟和環(huán)境協(xié)調(diào)發(fā)展研究[J]. 內(nèi)蒙古大學(xué)學(xué)報(哲學(xué)社會科學(xué)版), 2020,52(5):42-49. Zou Y F, Duan W, Guo N N. Research on Coordinated Development of Economy and Environment in Hohhot-Baotou-Erdos-Yulin Economic Zone [J]. Journal of Inner Mongolia University (Philosophy and Social Sciences), 2020,52(5):42-49.

    [29] Hoekstra A Y. Virtual water an introduction. Virtual water trade: Proceedings of the International expert meeting on virtual water trade. Value of water research report series No. 12 [R]. Delft. The Netherlands, IHE, 2003:12-13.

    [30] Water footprint network. Water footprint assessment manual - Setting the global standard. 2019. https://www.waterfootprint.org/en/water- footprint/global-water-footprint-standard/

    [31] Cai J, He Y, Xie R, et al. A footprint-based water security assessment: An analysis of Hunan province in China [J]. Journal of Cleaner Production, 2020,245,118485.

    [32] Bunsen J, Berger M, Finkbeiner M. Planetary boundaries for water - A review [J]. Ecological Indicators, 2021,121,107022.

    [33] Lucas P L, Wilting H C, Hof A F, et al. Allocating planetary boundaries to large economies: Distributional consequences of alternative perspectives on distributive fairness [J]. Global Environmental Change, 2020,60,102017.

    [34] Steffen W, Richardson K, Rockstrm J, et al. Planetary boundaries: Guiding human development on a changing planet [J]. Science, 2015,347,1259855.

    [35] 方 愷.基于足跡家族和行星邊界的主要國家環(huán)境可持續(xù)性多維評價[J]. 生態(tài)環(huán)境學(xué)報, 2014,11:1868-1875. Fang K. Multidimensional assessment of national environmental sustainability based on footprint family and planetary boundaries [J]. Ecology and Environmental Sciences, 2014,23(11):1868-1875.

    [36] 孫艷芝,沈 鐳.關(guān)于我國四大足跡理論研究變化的文獻計量分析[J]. 自然資源學(xué)報, 2016,31(9):1463-1473. Sun Y Z, Shen L. Bibliometric analysis on research progress of four footprint methodologies in China [J]. Journal of Natural Resources, 2016,31(9):1463-1473.

    [37] Eggleston H S, Buendia L, Miwa K, et al. IPCC guidelines for national greenhouse gas inventories [J]. Hayama: Institute for Global Environmental Strategies, 2006,2:48-56.

    [38] 張 寧,楊 肖,陳 彤.中國西部地區(qū)水-能源-糧食系統(tǒng)耦合協(xié)調(diào)度的時空特征[J]. 中國環(huán)境科學(xué), 2022,42(9):4444-4456. Zhang N, Yang X, Chen T. Research on the coupling coordination of water-energy-food system and its temporal and spatial characteristics [J]. China Environmental Science, 2022,42(9):4444- 4456.

    [39] Xian C, Fan Y, Zhang J, et al. Assessing sustainable water utilization from a holistic view: A case study of Guangdong, China [J]. Sustainable Cities and Society, 2022,76,103428.

    [40] Tomás M, López LA, Monsalve F. Carbon footprint, municipality size and rurality in Spain: Inequality and carbon taxation [J]. Journal of Cleaner Production, 2020,266,121798.

    [41] 趙武生,石培基.基于InVEST模型的復(fù)合生態(tài)系統(tǒng)耦合協(xié)調(diào)關(guān)系研究——以蘭西城市群為例[J]. 中國環(huán)境科學(xué), 2023,43(4):1883-1894. Zhao W S, Shi P J. The coupling and coordination of complex ecosystem based on the InVEST model- A case study in the Lanzhou- Xining urban agglomeration [J]. China Environmental Science, 2023, 43(4):1883-1894.

    [42] 徐 輝,王億文,張宗艷,等.黃河流域水-能源-糧食耦合機理及協(xié)調(diào)發(fā)展時空演變[J]. 資源科學(xué), 2021,43(12):2526-2537.Xu H, Wang Y W, Zhang Z Y, et al. Coupling mechanism of water- energy-food and spatiotemporal evolution of coordinated development in the Yellow River Basin [J]. Resources Science, 2021,43(12):2526-2537.

    [43] Li W, Wang Y, Xie S, et al. Coupling coordination analysis and spatiotemporal heterogeneity between urbanization and ecosystem health in Chongqing municipality, China [J]. Science of The Total Environment, 2021,791,148311.

    [44] 劉宇峰,原志華,郭玲霞,等.陜西省城市綠色增長水平時空演變特征及影響因素解析[J]. 自然資源學(xué)報, 2022,37(1):200-220. Liu Y F, Yuan Z H, Guo L X, et al. Spatio-temporal characteristics of urban green growth level and its influencing factors in Shaanxi province [J]. Journal of Natural Resources, 2022,37(1):200-220.

    [45] 郭向陽,穆學(xué)青,丁正山,等.長三角多維城市化對PM2.5濃度的非線性影響及驅(qū)動機制[J]. 地理學(xué)報, 2021,76(5):1274-1293. Guo X Y, Mu X Q, Ding Z S, et al. Nonlinear effects and driving mechanism of multidimensional urbanization on PM2.5concentrations in the Yangtze River Delta [J]. Acta Geographica Sinica, 2021,76(5): 1274-1293.

    [46] 梁建飛,陳松林.環(huán)境約束下的福建省城市建設(shè)用地利用效率及驅(qū)動因素[J]. 自然資源學(xué)報, 2020,35(12):2862-2874. Liang J F, Chen S L. Research on the land-use efficiency and driving factors of urban construction in Fujian province under environmental constraints [J]. Journal of Natural Resources, 2020,35(12):2862-2874.

    [47] 樊鵬飛,馮淑怡,蘇 敏,等.基于非期望產(chǎn)出的不同職能城市土地利用效率分異及驅(qū)動因素探究[J]. 資源科學(xué), 2018,40(5):946-957. Fan P F, Feng S Y, Su M, et al. Differential characteristics and driving factors of land use efficiency in different functional cities based on undesirable outputs [J]. Resources Science, 2018,40(5):946-957.

    [48] 劉 寧,李華姣,邊志強,等.基于空間面板STIRPAT模型的水足跡影響因素分析——以山東省為例[J]. 生態(tài)學(xué)報, 2022,42(22):9335- 9347. Liu N, Li H J, Bian Z Q. Influence factors of water footprint based on spatial panel STIRPAT Model: A case study of Shandong Province [J]. Acta Ecologica Sinica, 2022,42(22):9335-9347.

    [49] 李 波,嚴建飛.黃河流域水-能源-環(huán)境系統(tǒng)動態(tài)耦合協(xié)調(diào)發(fā)展研究[J]. 人民黃河, 2022,44(7):59-63. Li B, Yan J F. Research on the dynamic coupling coordinative degree among water-energy-environment system in the Yellow River Basin [J]. Yellow River, 2022,44(7):59-63.

    [50] 趙金輝,田 林,李思源,等.黃河流域能源與環(huán)境—經(jīng)濟—生態(tài)耦合協(xié)調(diào)發(fā)展研究[J]. 人民黃河, 2022,44(11):13-19.Zhao J H, Tian L, Li S Y, et al. Research on the coupling and coordination development of energy and environment-economy-ecology in the Yellow River Basin [J]. Yellow River, 2022,44(11):13-19.

    [51] 鄭賀允,葛力銘.資源型城市可持續(xù)發(fā)展對碳排放的影響研究——基于資源依賴的視角[J]. 中國環(huán)境科學(xué), 2022,42(6):2955-2964. Zheng H Y, Ge L M. Impact of sustainable development in resource-based cities on carbon emissions: From the perspective of resource dependence [J]. China Environmental Science, 2022,42(6): 2955-2964.

    [52] 秦炳濤,彭 涔,葛力銘,等.資源依賴、政府廉政水平與綠色技術(shù)創(chuàng)新——來自中國資源型城市的經(jīng)驗證據(jù)[J]. 中國環(huán)境科學(xué), 2023,43(7):3835-3847. Qin B T, Peng C, Ge L M, et al. Resource dependence, government integrity and green technology innovation——Evidence from resource-based cities in China [J]. China Environmental Science, 2023, 43(7):3835-3847.

    Analysis of the coupling process and the influencing factors of the “water-energy-carbon” relationship in the Hohhot-Baotou-Ordos-Yulin resource-based urban agglomeration.

    YANG Yi*, ZHANG Yuan-yuan

    (School of Economics and Management, Xi’an University of Technology, Xi’an 710054, China)., 2023,43(11):6212~6224

    To describe the spatiotemporal evolution of “water-energy-carbon” in the Hohhot-Baotou-Ordos-Yulin urban agglomeration, the level of coupled and coordinated development was analyzed, and a gray correlation model was established to identify the key driving factors affecting the change in the coupled and coordinated levels. The results showed that the water consumption, energy consumption and carbon emissions of urban agglomerations had fluctuating growth from 2005 to 2019, and the coupling coordination of “water-energy-carbon” was low, with urbanization rate, population size and economic development level being the key driving factors affecting coupling coordination. The key drivers of coupling coordination of “water-energy-carbon” in each city had obvious spatial heterogeneity. This study provides a scientific basis for revealing the dynamic change process and driving mechanism of urban ecosystems and formulating differentiated urban environmental management policies in resource-based urban agglomerations in the Yellow River Basin.

    Yellow River Basin;resource-based urban agglomeration;“water-energy-carbon”;coupling coordination;influencing factor

    X321,X24

    A

    1000-6923(2023)11-6212-13

    楊 屹(1969-),男,貴州鳳岡人,教授,博士,主要從事生態(tài)環(huán)境管理與可持續(xù)發(fā)展研究.發(fā)表論文87篇.yangyi_ nwpu@xaut.edu.cn.

    楊 屹,張園園.呼包鄂榆資源型城市群“水-能-碳”耦合變化及影響因素 [J]. 中國環(huán)境科學(xué), 2023,43(11):6212-6224.

    Yang Y, Zhang Y Y.Analysis of the coupling process and the influencing factors of the “water-energy-carbon” relationship in the Hohhot-Baotou-Ordos-Yulin resource-based urban agglomeration [J]. China Environmental Science, 2023,43(11):6212-6224.

    2023-03-13

    教育部人文社會科學(xué)研究規(guī)劃基金項目(21YJA790073)

    * 責任作者, 教授, yangyi_nwpu@xaut.edu.cn

    猜你喜歡
    足跡城市群耦合
    非Lipschitz條件下超前帶跳倒向耦合隨機微分方程的Wong-Zakai逼近
    成長足跡
    長三角城市群今年將有很多大動作
    我國第7個城市群建立
    建筑科技(2018年1期)2018-02-16 04:05:36
    足跡
    心聲歌刊(2018年6期)2018-01-24 00:56:12
    把省會城市群打造成強增長極
    足跡
    社會與公益(2016年2期)2016-04-13 02:49:06
    從國外經(jīng)驗看我國城市群一體化組織與管理
    基于“殼-固”耦合方法模擬焊接裝配
    大型鑄鍛件(2015年5期)2015-12-16 11:43:20
    春的足跡
    小主人報(2015年5期)2015-02-28 20:43:29
    一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩人妻高清精品专区| 免费无遮挡裸体视频| 三级经典国产精品| 精品国产露脸久久av麻豆 | 亚洲性久久影院| 91精品一卡2卡3卡4卡| 少妇的逼水好多| 久久久精品大字幕| 国产色爽女视频免费观看| 国产欧美日韩精品一区二区| 午夜免费男女啪啪视频观看| 久久精品国产亚洲av涩爱| 国产探花极品一区二区| 老司机影院毛片| 天堂中文最新版在线下载 | 天堂√8在线中文| 亚洲精品色激情综合| 久久久久久大精品| 成人av在线播放网站| 国产色爽女视频免费观看| 夫妻性生交免费视频一级片| 亚洲综合色惰| 国产精华一区二区三区| 观看免费一级毛片| 欧美日韩综合久久久久久| 亚洲欧洲国产日韩| 日韩欧美精品免费久久| 十八禁国产超污无遮挡网站| 成年女人看的毛片在线观看| АⅤ资源中文在线天堂| 亚洲国产色片| 国产一区有黄有色的免费视频 | 欧美区成人在线视频| 日韩中字成人| 国产免费视频播放在线视频 | 99久久精品一区二区三区| 国产免费福利视频在线观看| 亚洲熟妇中文字幕五十中出| 春色校园在线视频观看| 观看免费一级毛片| 一级黄色大片毛片| 99久久精品热视频| 综合色av麻豆| 99久久精品国产国产毛片| 国产精品一二三区在线看| 十八禁国产超污无遮挡网站| 免费看光身美女| 欧美激情国产日韩精品一区| 亚洲最大成人中文| 18禁裸乳无遮挡免费网站照片| 亚洲av熟女| 毛片一级片免费看久久久久| 国产淫片久久久久久久久| 久热久热在线精品观看| 免费看av在线观看网站| 国产综合懂色| 国产成人一区二区在线| 高清午夜精品一区二区三区| 成人毛片a级毛片在线播放| 97人妻精品一区二区三区麻豆| 搡女人真爽免费视频火全软件| 免费电影在线观看免费观看| 国产一区二区在线av高清观看| 18禁在线播放成人免费| 禁无遮挡网站| 最近的中文字幕免费完整| 99久久成人亚洲精品观看| 中国美白少妇内射xxxbb| 久久精品91蜜桃| 日韩精品青青久久久久久| 亚洲精品日韩在线中文字幕| 久久婷婷人人爽人人干人人爱| 亚洲18禁久久av| av在线播放精品| 亚洲国产色片| 特大巨黑吊av在线直播| 日本欧美国产在线视频| 成人无遮挡网站| 国内精品宾馆在线| 亚洲四区av| 国产在线一区二区三区精 | 国产91av在线免费观看| 嘟嘟电影网在线观看| 久久草成人影院| 久久久国产成人精品二区| 成人亚洲精品av一区二区| 国产乱人偷精品视频| 少妇的逼水好多| 欧美区成人在线视频| 自拍偷自拍亚洲精品老妇| 菩萨蛮人人尽说江南好唐韦庄 | 97热精品久久久久久| 免费在线观看成人毛片| 高清毛片免费看| 国产精品不卡视频一区二区| 国产免费福利视频在线观看| 久久久久精品久久久久真实原创| 欧美激情在线99| 亚洲不卡免费看| 精品不卡国产一区二区三区| 国产一区有黄有色的免费视频 | 亚洲人成网站在线播| 一级爰片在线观看| 久久久久久久久久成人| 亚洲欧美精品自产自拍| 婷婷六月久久综合丁香| 欧美一区二区国产精品久久精品| 99国产精品一区二区蜜桃av| 精品欧美国产一区二区三| 国产v大片淫在线免费观看| 久久久精品大字幕| 国产白丝娇喘喷水9色精品| 日韩av在线免费看完整版不卡| 男人舔女人下体高潮全视频| 国产一区二区在线av高清观看| 日韩av在线大香蕉| 色视频www国产| 国产精品乱码一区二三区的特点| 国产精品av视频在线免费观看| 丝袜美腿在线中文| 免费观看精品视频网站| 久久精品人妻少妇| 91在线精品国自产拍蜜月| 成年av动漫网址| 精品不卡国产一区二区三区| 国产一区二区亚洲精品在线观看| 51国产日韩欧美| 日韩在线高清观看一区二区三区| 黄片wwwwww| 亚洲高清免费不卡视频| 秋霞伦理黄片| 久久精品夜夜夜夜夜久久蜜豆| 卡戴珊不雅视频在线播放| 中文字幕人妻熟人妻熟丝袜美| 在线观看66精品国产| 亚洲av福利一区| 中文字幕精品亚洲无线码一区| 亚洲成av人片在线播放无| av国产免费在线观看| 97热精品久久久久久| 国产色婷婷99| 少妇人妻一区二区三区视频| av国产久精品久网站免费入址| 直男gayav资源| 精品熟女少妇av免费看| 午夜亚洲福利在线播放| 亚洲av电影不卡..在线观看| 亚洲欧美精品自产自拍| 狠狠狠狠99中文字幕| 久久久久久久久久久丰满| 永久网站在线| 国产精品精品国产色婷婷| 国产单亲对白刺激| av福利片在线观看| 成人漫画全彩无遮挡| 国产91av在线免费观看| 中文字幕av成人在线电影| 青春草亚洲视频在线观看| 亚洲av电影不卡..在线观看| 三级国产精品欧美在线观看| 日日干狠狠操夜夜爽| 中文字幕免费在线视频6| 九九在线视频观看精品| 舔av片在线| 午夜精品国产一区二区电影 | 久久99热6这里只有精品| 内射极品少妇av片p| 色综合亚洲欧美另类图片| 69av精品久久久久久| 别揉我奶头 嗯啊视频| 精品一区二区三区人妻视频| 偷拍熟女少妇极品色| 免费黄网站久久成人精品| 亚洲欧美精品自产自拍| 69av精品久久久久久| 观看免费一级毛片| 只有这里有精品99| 大香蕉97超碰在线| 少妇高潮的动态图| 长腿黑丝高跟| 欧美成人精品欧美一级黄| 亚洲国产高清在线一区二区三| 一边亲一边摸免费视频| 欧美丝袜亚洲另类| 欧美日韩一区二区视频在线观看视频在线 | 国产精华一区二区三区| 亚洲精华国产精华液的使用体验| 久久久久精品久久久久真实原创| 插阴视频在线观看视频| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久成人| 床上黄色一级片| 桃色一区二区三区在线观看| 美女脱内裤让男人舔精品视频| 一边亲一边摸免费视频| 在线播放国产精品三级| 久久99精品国语久久久| 美女cb高潮喷水在线观看| 日韩在线高清观看一区二区三区| 男女下面进入的视频免费午夜| 国模一区二区三区四区视频| 亚洲国产日韩欧美精品在线观看| 草草在线视频免费看| 国产精品熟女久久久久浪| 色5月婷婷丁香| 国产成人91sexporn| 亚洲欧美一区二区三区国产| 亚洲经典国产精华液单| 午夜亚洲福利在线播放| 国产探花在线观看一区二区| 亚洲av成人精品一二三区| 国产黄色视频一区二区在线观看 | 97热精品久久久久久| 午夜免费男女啪啪视频观看| 亚洲国产欧洲综合997久久,| 国产91av在线免费观看| 午夜日本视频在线| 国产片特级美女逼逼视频| 国产精品乱码一区二三区的特点| 国产黄a三级三级三级人| 久久精品久久久久久久性| 亚洲中文字幕日韩| 日日摸夜夜添夜夜爱| 大又大粗又爽又黄少妇毛片口| 国产午夜精品久久久久久一区二区三区| 天天躁日日操中文字幕| 亚洲国产日韩欧美精品在线观看| 黑人高潮一二区| 久久久国产成人精品二区| 国产av码专区亚洲av| 中文字幕熟女人妻在线| 尤物成人国产欧美一区二区三区| 免费看a级黄色片| 国语对白做爰xxxⅹ性视频网站| 99久久精品一区二区三区| 99久久精品热视频| 亚洲伊人久久精品综合 | 国产精品一区二区三区四区久久| 国产黄片视频在线免费观看| 成人鲁丝片一二三区免费| 国产免费男女视频| 欧美日韩在线观看h| 69av精品久久久久久| 成人亚洲精品av一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 汤姆久久久久久久影院中文字幕 | 久久精品夜夜夜夜夜久久蜜豆| 久久6这里有精品| 国产精品久久久久久久久免| 国产黄片美女视频| 男女下面进入的视频免费午夜| 亚洲aⅴ乱码一区二区在线播放| 成人亚洲精品av一区二区| www.av在线官网国产| 亚洲美女视频黄频| 中文字幕精品亚洲无线码一区| 久久久久久久久久成人| or卡值多少钱| 国产乱人偷精品视频| 中文欧美无线码| 在线免费十八禁| 国产精品一区二区三区四区免费观看| 国产高清国产精品国产三级 | 久久久成人免费电影| 日本黄色片子视频| 亚洲av中文av极速乱| 毛片女人毛片| 亚洲经典国产精华液单| 女人十人毛片免费观看3o分钟| 免费人成在线观看视频色| 色综合亚洲欧美另类图片| 亚洲欧美一区二区三区国产| av在线蜜桃| 欧美成人a在线观看| kizo精华| 成人特级av手机在线观看| 国产色爽女视频免费观看| 色噜噜av男人的天堂激情| 51国产日韩欧美| av.在线天堂| 日日摸夜夜添夜夜爱| 亚洲国产日韩欧美精品在线观看| 久久人人爽人人爽人人片va| 国产高清三级在线| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av日韩在线播放| 亚洲精品影视一区二区三区av| 18+在线观看网站| 爱豆传媒免费全集在线观看| 国产乱人偷精品视频| 精华霜和精华液先用哪个| 亚洲国产精品成人久久小说| 国产毛片a区久久久久| 亚洲国产欧美在线一区| 午夜老司机福利剧场| 一个人免费在线观看电影| 午夜免费激情av| 亚洲国产精品久久男人天堂| 欧美激情久久久久久爽电影| 村上凉子中文字幕在线| 一个人观看的视频www高清免费观看| 毛片一级片免费看久久久久| 熟女人妻精品中文字幕| 成人性生交大片免费视频hd| 又爽又黄无遮挡网站| 日本黄色片子视频| 九九爱精品视频在线观看| 免费观看精品视频网站| 国产伦精品一区二区三区四那| 日韩视频在线欧美| 小蜜桃在线观看免费完整版高清| 欧美变态另类bdsm刘玥| 欧美高清成人免费视频www| 日韩中字成人| 亚州av有码| 最近的中文字幕免费完整| 国产成人freesex在线| av播播在线观看一区| 一区二区三区乱码不卡18| 看片在线看免费视频| 天天躁夜夜躁狠狠久久av| 一级二级三级毛片免费看| 亚洲欧美精品综合久久99| 白带黄色成豆腐渣| 日韩欧美精品v在线| 国产黄色小视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲精品一区蜜桃| 精品久久久噜噜| 中国美白少妇内射xxxbb| 精品国产一区二区三区久久久樱花 | 精品午夜福利在线看| 一区二区三区免费毛片| 亚洲av中文字字幕乱码综合| 久久精品国产自在天天线| 国产精品嫩草影院av在线观看| 中文在线观看免费www的网站| 男的添女的下面高潮视频| 国产伦精品一区二区三区视频9| 丰满乱子伦码专区| 午夜精品国产一区二区电影 | 国产真实伦视频高清在线观看| 国产精品久久久久久精品电影| 99久久精品国产国产毛片| 91精品一卡2卡3卡4卡| 成人国产麻豆网| 在线免费观看不下载黄p国产| 国产免费男女视频| 国产精品久久久久久精品电影| 在线免费十八禁| 亚洲国产精品合色在线| 亚洲伊人久久精品综合 | 老师上课跳d突然被开到最大视频| 男的添女的下面高潮视频| 欧美激情在线99| 国产精品电影一区二区三区| 国产精品无大码| 国产亚洲精品久久久com| 欧美激情国产日韩精品一区| 久久久久久久久久黄片| 精品无人区乱码1区二区| 女人久久www免费人成看片 | av在线老鸭窝| 小说图片视频综合网站| 校园人妻丝袜中文字幕| 少妇高潮的动态图| 丰满少妇做爰视频| 日本黄色视频三级网站网址| 91在线精品国自产拍蜜月| 两性午夜刺激爽爽歪歪视频在线观看| 日本免费a在线| 国产精品不卡视频一区二区| 我要搜黄色片| 日韩国内少妇激情av| 青春草亚洲视频在线观看| 91精品一卡2卡3卡4卡| 舔av片在线| 国产伦理片在线播放av一区| 亚洲av一区综合| 国产伦在线观看视频一区| 久久久精品94久久精品| 欧美日韩国产亚洲二区| 日本av手机在线免费观看| 亚洲最大成人中文| 日本爱情动作片www.在线观看| 欧美日本视频| 99九九线精品视频在线观看视频| 亚洲无线观看免费| 国产黄色小视频在线观看| 国产伦在线观看视频一区| 最近视频中文字幕2019在线8| 麻豆一二三区av精品| 国产成人福利小说| 国产精品不卡视频一区二区| 久久99热这里只频精品6学生 | 九色成人免费人妻av| 午夜精品一区二区三区免费看| av免费在线看不卡| 国产日韩欧美在线精品| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看| 日本免费a在线| 伦精品一区二区三区| 青春草亚洲视频在线观看| 国产伦在线观看视频一区| 一级二级三级毛片免费看| 内射极品少妇av片p| 欧美激情在线99| .国产精品久久| 日韩欧美精品免费久久| 国产免费一级a男人的天堂| 亚洲国产精品成人久久小说| 日本五十路高清| 国产69精品久久久久777片| 黄片无遮挡物在线观看| 老司机福利观看| 国产伦精品一区二区三区视频9| 九九爱精品视频在线观看| 国产精品久久久久久精品电影| 99久久人妻综合| 七月丁香在线播放| 美女被艹到高潮喷水动态| 男女国产视频网站| 天天一区二区日本电影三级| 亚洲国产精品成人综合色| 成年版毛片免费区| 2021少妇久久久久久久久久久| 永久网站在线| 韩国av在线不卡| 久久国产乱子免费精品| 免费黄网站久久成人精品| 99久久中文字幕三级久久日本| 久久人妻av系列| 久久精品影院6| 可以在线观看毛片的网站| 国产亚洲91精品色在线| 国产亚洲精品av在线| 偷拍熟女少妇极品色| 亚洲成人中文字幕在线播放| 亚洲av电影不卡..在线观看| 久久婷婷人人爽人人干人人爱| 国产乱人偷精品视频| 五月玫瑰六月丁香| 国产亚洲精品av在线| 视频中文字幕在线观看| 在线a可以看的网站| 亚洲欧美日韩高清专用| 国产三级在线视频| 极品教师在线视频| 精品一区二区免费观看| 亚洲丝袜综合中文字幕| 欧美一区二区亚洲| 在线观看66精品国产| 精品欧美国产一区二区三| 中文天堂在线官网| 久久草成人影院| 午夜a级毛片| 久热久热在线精品观看| 日韩大片免费观看网站 | 97人妻精品一区二区三区麻豆| 国产精品蜜桃在线观看| 免费电影在线观看免费观看| 啦啦啦韩国在线观看视频| 色哟哟·www| 特级一级黄色大片| 18禁裸乳无遮挡免费网站照片| 亚洲国产高清在线一区二区三| 好男人在线观看高清免费视频| 在线观看av片永久免费下载| 蜜桃亚洲精品一区二区三区| 丰满人妻一区二区三区视频av| 国产成人a区在线观看| 国国产精品蜜臀av免费| 少妇熟女aⅴ在线视频| 日韩av在线大香蕉| 日本免费在线观看一区| 久久精品国产亚洲网站| 亚洲美女视频黄频| av专区在线播放| 精品免费久久久久久久清纯| 伦理电影大哥的女人| 大香蕉97超碰在线| av播播在线观看一区| 久久久久久久久久久丰满| 国产黄a三级三级三级人| 亚洲欧美成人精品一区二区| 男女视频在线观看网站免费| 国产高潮美女av| 国产一级毛片在线| av免费在线看不卡| 亚洲欧美成人精品一区二区| 亚洲精品影视一区二区三区av| 久久精品91蜜桃| 亚洲人与动物交配视频| 日本午夜av视频| 老师上课跳d突然被开到最大视频| 又粗又硬又长又爽又黄的视频| 国产在线男女| 日韩av在线大香蕉| 国产一级毛片七仙女欲春2| 国产视频内射| 26uuu在线亚洲综合色| 欧美一级a爱片免费观看看| 偷拍熟女少妇极品色| 小说图片视频综合网站| 99热全是精品| 久久99热这里只有精品18| 黄色欧美视频在线观看| 综合色丁香网| 日本免费在线观看一区| 久久99热6这里只有精品| 国产成人91sexporn| 看黄色毛片网站| 两个人视频免费观看高清| 九九在线视频观看精品| a级毛色黄片| 草草在线视频免费看| 天堂网av新在线| 日本免费a在线| 在线a可以看的网站| 免费黄色在线免费观看| 亚洲怡红院男人天堂| a级一级毛片免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 你懂的网址亚洲精品在线观看 | 免费观看a级毛片全部| 老女人水多毛片| 一级二级三级毛片免费看| 精品国产露脸久久av麻豆 | 国产精品一二三区在线看| 亚洲综合色惰| 国产伦理片在线播放av一区| 男插女下体视频免费在线播放| 久久久久久九九精品二区国产| 婷婷色av中文字幕| 久久久久久久久久黄片| 国产精品不卡视频一区二区| 亚洲va在线va天堂va国产| av在线观看视频网站免费| 色尼玛亚洲综合影院| 久久久久久大精品| 看十八女毛片水多多多| 热99re8久久精品国产| 菩萨蛮人人尽说江南好唐韦庄 | 日本wwww免费看| 中文字幕av在线有码专区| 午夜福利在线观看吧| 亚洲精品久久久久久婷婷小说 | 亚洲自拍偷在线| 天堂影院成人在线观看| 91久久精品国产一区二区三区| 伊人久久精品亚洲午夜| av女优亚洲男人天堂| 久久精品久久久久久噜噜老黄 | 青青草视频在线视频观看| 精品熟女少妇av免费看| 欧美成人精品欧美一级黄| 九色成人免费人妻av| 免费无遮挡裸体视频| 亚洲欧美中文字幕日韩二区| 搡女人真爽免费视频火全软件| 成人一区二区视频在线观看| 91在线精品国自产拍蜜月| 欧美一区二区精品小视频在线| 男人狂女人下面高潮的视频| 久久人人爽人人爽人人片va| 伦理电影大哥的女人| 少妇人妻精品综合一区二区| 高清日韩中文字幕在线| 亚洲性久久影院| 国产精品熟女久久久久浪| 日本色播在线视频| 精品久久久噜噜| 91久久精品国产一区二区成人| 人妻夜夜爽99麻豆av| 精品久久久久久成人av| 亚洲欧洲日产国产| 亚洲成人av在线免费| 日本黄色片子视频| 亚洲av成人av| 亚洲国产精品合色在线| 亚洲国产精品成人久久小说| 草草在线视频免费看| 欧美变态另类bdsm刘玥| 亚洲va在线va天堂va国产| 亚洲精品日韩av片在线观看| 联通29元200g的流量卡| 亚洲精品影视一区二区三区av| 国产乱人视频| 最近最新中文字幕大全电影3| 一级二级三级毛片免费看| 精品一区二区三区人妻视频| 免费黄网站久久成人精品| 我要看日韩黄色一级片| 国产午夜福利久久久久久| 麻豆av噜噜一区二区三区| 中文精品一卡2卡3卡4更新| 白带黄色成豆腐渣| 亚洲人与动物交配视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99久国产av精品国产电影| 26uuu在线亚洲综合色| 亚洲精华国产精华液的使用体验| 欧美成人一区二区免费高清观看| 久久久午夜欧美精品| 国产成人精品一,二区| 一二三四中文在线观看免费高清| ponron亚洲| 别揉我奶头 嗯啊视频| 综合色丁香网| 久久久亚洲精品成人影院| 国产色爽女视频免费观看|