供需關(guān)系下長(zhǎng)三角城市群 生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)

陳燕，譚堰琴，邱曉敏，宋新山，周鎮(zhèn)宇，萬(wàn)榮榮

(1.東華大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 201620；2.中國(guó)科學(xué)院 流域地理學(xué)重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 210008；3.密蘇里州立大學(xué) 地理地質(zhì)規(guī)劃系，伊利諾伊 斯普林菲爾德 65897)

隨著人口劇增與快速城市化，城市群生態(tài)系統(tǒng)的結(jié)構(gòu)和功能均受到人類(lèi)活動(dòng)的干擾[1]，區(qū)域生態(tài)系統(tǒng)服務(wù)供需失衡，威脅著生態(tài)安全與可持續(xù)發(fā)展.生態(tài)系統(tǒng)服務(wù)供需平衡是區(qū)域生態(tài)環(huán)境健康程度的重要衡量指標(biāo)之一[2]，而生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)是實(shí)現(xiàn)生態(tài)服務(wù)供給和人類(lèi)需求耦合的重要橋梁[3]，對(duì)生態(tài)系統(tǒng)服務(wù)的傳輸、轉(zhuǎn)化和消費(fèi)具有重要的作用[4].因此，基于供給和需求關(guān)系的生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)研究能夠明確生態(tài)系統(tǒng)服務(wù)的產(chǎn)生源頭、耗散過(guò)程及受益區(qū)域，有助于優(yōu)化生態(tài)系統(tǒng)服務(wù)管理，改善區(qū)域生態(tài)環(huán)境狀況.

生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)是指重要生態(tài)服務(wù)功能通過(guò)一定的流通介質(zhì)流轉(zhuǎn)到區(qū)域外的地區(qū)，并產(chǎn)生效應(yīng)[5-6]，是生態(tài)系統(tǒng)服務(wù)供給與需求關(guān)聯(lián)研究的深化[7]，關(guān)于這方面的研究尚處于起步階段，并向空間制圖與定量分析的方向發(fā)展[8-9].目前，生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)研究主要是對(duì)生態(tài)系統(tǒng)服務(wù)從產(chǎn)生到耗散過(guò)程中受益區(qū)接受的生態(tài)系統(tǒng)服務(wù)價(jià)值量進(jìn)行評(píng)估[10-11]，且主要集中于單一區(qū)域供給服務(wù)的流轉(zhuǎn)研究[12]，關(guān)于跨區(qū)域空間流轉(zhuǎn)及綜合考慮供需關(guān)系的研究仍然較少.

長(zhǎng)江三角洲(長(zhǎng)三角)城市群是我國(guó)人口最密集、經(jīng)濟(jì)最發(fā)達(dá)、城市化進(jìn)程最快的地區(qū)之一.2019年，國(guó)家政府工作報(bào)告提出，將長(zhǎng)三角區(qū)域一體化發(fā)展上升為國(guó)家戰(zhàn)略，建設(shè)生態(tài)型城市群[13].近年來(lái)，長(zhǎng)三角地區(qū)在經(jīng)濟(jì)快速發(fā)展的同時(shí)，出現(xiàn)了生物多樣性銳減、河流湖泊水質(zhì)惡化、城市空氣質(zhì)量下降等一系列生態(tài)環(huán)境問(wèn)題，城市生態(tài)空間遭受不同程度的破壞，其實(shí)質(zhì)是區(qū)域生態(tài)系統(tǒng)服務(wù)供需關(guān)系的空間差異或失衡[14].基于此，本文對(duì)供需關(guān)系下長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)進(jìn)行研究.

1 研究區(qū)域及數(shù)據(jù)來(lái)源

1.1 研究區(qū)域

長(zhǎng)三角城市群是由上海市及江蘇省、浙江省、安徽省的25個(gè)地級(jí)市組成(圖1)，面積為21.17萬(wàn)km2，地形以平原和丘陵為主，氣候類(lèi)型主要為亞熱帶季風(fēng)性濕潤(rùn)氣候.

圖1 長(zhǎng)三角城市群地理位置

2018年統(tǒng)計(jì)年鑒數(shù)據(jù)表明，長(zhǎng)三角城市群城鎮(zhèn)化率為68%；2018年土地利用分類(lèi)數(shù)據(jù)中，長(zhǎng)三角城市群建設(shè)用地比例約為20%，城市用地類(lèi)型及空間結(jié)構(gòu)變化顯著.快速的城市化在促進(jìn)社會(huì)經(jīng)濟(jì)高速發(fā)展的同時(shí)，也削弱了生態(tài)系統(tǒng)為城市提供的重要生態(tài)服務(wù)功能[15]，城市群空間的高強(qiáng)度開(kāi)發(fā)導(dǎo)致生態(tài)系統(tǒng)服務(wù)供需失衡等一系列環(huán)境問(wèn)題愈演愈烈，協(xié)調(diào)生態(tài)環(huán)境保護(hù)與區(qū)域經(jīng)濟(jì)發(fā)展成為長(zhǎng)三角城市群綠色一體化發(fā)展的重要任務(wù).

1.2 數(shù)據(jù)來(lái)源

使用的4類(lèi)數(shù)據(jù)來(lái)源如下.

1)土地利用數(shù)據(jù).該數(shù)據(jù)包括長(zhǎng)三角城市群2010年、2018年的土地利用數(shù)據(jù)(30 m×30 m)，來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http:∥www.resdc.cn).2010年的土地利用數(shù)據(jù)以當(dāng)年Landsat TM遙感影像為主信息源，通過(guò)人工目視解譯生成；2018年的土地利用數(shù)據(jù)是在2015年土地利用遙感監(jiān)測(cè)數(shù)據(jù)的基礎(chǔ)上，基于Landsat 8遙感影像，通過(guò)人工目視解譯生成.兩期遙感影像的獲取時(shí)間均為云覆蓋較少(云量小于 10%)的5—9月，解譯精度達(dá)90%以上.

2)社會(huì)經(jīng)濟(jì)數(shù)據(jù).該數(shù)據(jù)包括2010年、2018年全國(guó)與長(zhǎng)三角城市群城鄉(xiāng)居民人口、人均可支配收入、地區(qū)生產(chǎn)總值、糧食產(chǎn)值與產(chǎn)量等，來(lái)源于國(guó)家統(tǒng)計(jì)局、長(zhǎng)三角城市群各市統(tǒng)計(jì)局.

3)植被凈初級(jí)生產(chǎn)力(NPP)、年降水量數(shù)據(jù).該數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心.

4)長(zhǎng)三角城市群主要河流、主要道路、行政邊界等矢量數(shù)據(jù).該數(shù)據(jù)來(lái)源于長(zhǎng)江三角洲科學(xué)數(shù)據(jù)中心(http:∥nnu.geodata.cn:8008/#).

2 研究方法

2.1 生態(tài)系統(tǒng)服務(wù)供需核算

2.1.1 生態(tài)系統(tǒng)服務(wù)供給價(jià)值核算 以長(zhǎng)三角城市群的區(qū)/縣為供給計(jì)算單元，主要依據(jù)文獻(xiàn)[16-17]的研究成果進(jìn)行生態(tài)系統(tǒng)服務(wù)供給價(jià)值(Pes)核算，有

(1)

式(1)中：EVj為生態(tài)系統(tǒng)第j項(xiàng)服務(wù)功能的單位面積總價(jià)值，元·km-2；Pno為無(wú)人力投入的自然生態(tài)系統(tǒng)的單位面積農(nóng)田糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值，元·km-2，通常認(rèn)為無(wú)人力投入的自然生態(tài)系統(tǒng)的單位面積農(nóng)田糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值約為單位面積農(nóng)田提供糧食作物生產(chǎn)服務(wù)經(jīng)濟(jì)價(jià)值的1/7；ei，j為第i類(lèi)土地利用類(lèi)型生態(tài)系統(tǒng)的第j項(xiàng)服務(wù)功能相對(duì)于農(nóng)田生態(tài)系統(tǒng)提供生態(tài)服務(wù)價(jià)值的當(dāng)量系數(shù)；Pg為研究區(qū)每年的糧食總價(jià)值，元；Ag為糧食播種面積，km2.

由于采用全國(guó)參數(shù)表征長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)供給價(jià)值具有一定的誤差，因此，對(duì)研究區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量系數(shù)進(jìn)行調(diào)整：采用社會(huì)經(jīng)濟(jì)調(diào)節(jié)因子修正糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值[18]，采用NPP時(shí)空調(diào)節(jié)因子和降水時(shí)空調(diào)節(jié)因子修正價(jià)值當(dāng)量系數(shù)[17].

經(jīng)計(jì)算，社會(huì)經(jīng)濟(jì)調(diào)節(jié)因子、NPP時(shí)空調(diào)節(jié)因子和降水時(shí)空調(diào)節(jié)因子分別為1.32，1.73，1.57.修正的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量系數(shù)(er)，如表1所示.通過(guò)修正的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量系數(shù)，可得長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)供給價(jià)值.

表1 修正的單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量系數(shù)

2.1.2 生態(tài)系統(tǒng)服務(wù)需求核算 在文獻(xiàn)[19]的基礎(chǔ)上，將生態(tài)系統(tǒng)服務(wù)需求理解為人類(lèi)消耗或偏好獲取的生態(tài)系統(tǒng)服務(wù)數(shù)量，考慮到長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)變化的驅(qū)動(dòng)因素主要為土地開(kāi)發(fā)利用的強(qiáng)度與人類(lèi)活動(dòng)的干擾程度[20]，選取生態(tài)系統(tǒng)服務(wù)需求指數(shù)表征研究區(qū)域的生態(tài)系統(tǒng)服務(wù)需求.

研究區(qū)域的生態(tài)系統(tǒng)服務(wù)需求指數(shù)(Ied)的計(jì)算公式為

Ied=δ×lgρ×lgP.

(2)

式(2)中：δ為建設(shè)用地比例，%；ρ為人口密度，人·km-2；P為地均國(guó)內(nèi)生產(chǎn)總值(GDP)，萬(wàn)元·km-2.

2.2 生態(tài)系統(tǒng)服務(wù)供需空間格局

采用生態(tài)系統(tǒng)服務(wù)供需比率(RESD)量化分析供需平衡度，揭示研究區(qū)域生態(tài)系統(tǒng)服務(wù)供給和需求在數(shù)量上的盈余或赤字關(guān)系[21].當(dāng)RESD為正時(shí)，表示生態(tài)系統(tǒng)服務(wù)供過(guò)于求；當(dāng)RESD為0時(shí)，表示生態(tài)系統(tǒng)服務(wù)供需平衡；當(dāng)RESD為負(fù)時(shí)，表示生態(tài)系統(tǒng)服務(wù)供不應(yīng)求.

RESD的計(jì)算公式為

(3)

式(3)中：采用歸一化方法對(duì)生態(tài)系統(tǒng)服務(wù)的供給數(shù)據(jù)和需求數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化計(jì)算；S，D分別為研究區(qū)域各區(qū)/縣生態(tài)系統(tǒng)服務(wù)的供給值與需求值；Smax為生態(tài)系統(tǒng)服務(wù)供給最大值；Dmax為生態(tài)系統(tǒng)服務(wù)需求最大值.

2.3 基于供需關(guān)系的生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)

2.3.1 生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)路徑的確定 生態(tài)源地具有高供給價(jià)值，是生態(tài)系統(tǒng)服務(wù)流動(dòng)和傳遞的源頭[1]；生態(tài)需求地具有高需求指數(shù)，是生態(tài)系統(tǒng)服務(wù)的主要受益區(qū).生態(tài)廊道是由阻力較小的景觀類(lèi)型構(gòu)成的條帶狀通道，是區(qū)域內(nèi)物質(zhì)交換及能量流動(dòng)的載體[24].

基于生態(tài)系統(tǒng)服務(wù)供需空間格局確定生態(tài)源地與生態(tài)需求地，采用最小累積阻力(MCR)模型提取生態(tài)廊道[25].將熱點(diǎn)分析結(jié)果識(shí)別出的供需盈余區(qū)判定為生態(tài)源地，供需赤字區(qū)判定為生態(tài)需求地，利用 ArcGIS平臺(tái)的成本路徑工具，以所有生態(tài)源地質(zhì)心為源，以生態(tài)需求地質(zhì)心為目標(biāo)，生成阻力最小的成本路徑(最小成本路徑)，作為兩地之間的潛在生態(tài)廊道；通過(guò)重力模型計(jì)算生態(tài)源地與生態(tài)需求地之間的相互作用矩陣，提取最小成本路徑作為生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)的最優(yōu)路徑.

2.3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量的構(gòu)建 生態(tài)系統(tǒng)服務(wù)具有一定的方向性和區(qū)域性，不同位置發(fā)生的效用不同，發(fā)生流轉(zhuǎn)的價(jià)值量基本遵循距離衰減規(guī)律[12].根據(jù)生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)的特點(diǎn)，借鑒康維斯在引力模型基礎(chǔ)上提出的斷裂點(diǎn)公式，模擬生態(tài)系統(tǒng)服務(wù)從生態(tài)源地向生態(tài)需求地流動(dòng)的邊界，再根據(jù)場(chǎng)強(qiáng)理論計(jì)算單位面積內(nèi)生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)強(qiáng)度[6]，可得生態(tài)源地到生態(tài)需求地的生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量(Es，d).

Es，d的計(jì)算公式為

(4)

式(4)中：Is，d為生態(tài)系統(tǒng)服務(wù)從生態(tài)源地到生態(tài)需求地的流轉(zhuǎn)強(qiáng)度；A為生態(tài)系統(tǒng)服務(wù)從生態(tài)源地流轉(zhuǎn)到生態(tài)需求地的輻射面積，km2；ks，d為影響生態(tài)源地到生態(tài)需求地的生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)影響因子，取值為[0～1]，受水流、風(fēng)速及風(fēng)向等流動(dòng)因子的影響，ks，d一般取0.6[26]；Ds為生態(tài)源地質(zhì)心與斷裂點(diǎn)的距離，km；Ds，d為生態(tài)源地質(zhì)心與生態(tài)需求地質(zhì)心的距離，km；Es，Ed分別為生態(tài)源地與生態(tài)需求地的生態(tài)系統(tǒng)服務(wù)價(jià)值，元.

3 研究結(jié)果與分析

3.1 長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)供需空間特征

3.1.1 生態(tài)系統(tǒng)服務(wù)供給空間特征 2010年和2018年生態(tài)系統(tǒng)服務(wù)供給價(jià)值，如圖2所示.由圖2可知：2010年和2018年長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)供給價(jià)值分別為65 694，67 312元·km-2，鹽城、紹興、臺(tái)州等個(gè)別區(qū)/縣供給價(jià)值增加，常州和上海的部分區(qū)/縣供給價(jià)值減少，研究年份內(nèi)，生態(tài)系統(tǒng)服務(wù)供給價(jià)值變化微??；長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)供給呈南高北低的分布形態(tài)，供給高值區(qū)主要分布于杭州市淳安縣、安慶市宿松縣及蘇州市吳中區(qū)，該區(qū)域土地利用類(lèi)型主要為林地、草地和水域，提供生態(tài)系統(tǒng)服務(wù)的能力較強(qiáng)，蘇州市吳中區(qū)擁有超過(guò)180 km的親水岸，占太湖水域的3/5，水資源供給和水源涵養(yǎng)服務(wù)價(jià)值極高，供給價(jià)值達(dá)1 326元·km-2，是供給值最高的區(qū)/縣單元；供給低值區(qū)主要分布于南京市秦淮區(qū)、杭州市下城區(qū)及上海市區(qū)，供給最低值區(qū)為上海市靜安區(qū)(6.5 元·km-2)，土地利用類(lèi)型主要是生態(tài)服務(wù)功能極少的建設(shè)用地.

(a)2010年 (b)2018年

3.1.2 生態(tài)系統(tǒng)服務(wù)需求空間特征 2010年和2018年生態(tài)系統(tǒng)服務(wù)需求指數(shù)，如圖3所示.由圖3可知：2010年和2018年長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)需求指數(shù)分別為520.48，798.50，呈明顯的增長(zhǎng)趨勢(shì)，南京、紹興、金華、寧波等地區(qū)大部分區(qū)/縣需求指數(shù)顯著增加；生態(tài)系統(tǒng)服務(wù)需求與區(qū)域經(jīng)濟(jì)發(fā)展水平密切相關(guān)，各區(qū)/縣差異顯著，整體呈東高西低的分布形態(tài)；需求高值區(qū)主要分布于南京、杭州、合肥、上海等經(jīng)濟(jì)發(fā)達(dá)城市的中心城區(qū)，需求最高值區(qū)為上海市黃浦區(qū)(地均GDP高達(dá)1 107 630.17萬(wàn)元·km-2，人口密度為31 955人·km-2，建設(shè)用地比例為91.70%)；需求低值區(qū)主要分布于安慶、池州、宣城及杭州西南部的區(qū)/縣，需求最低值區(qū)為安慶市岳西縣(地均GDP為171.98萬(wàn)元·km-2，人口密度為136人·km-2，建設(shè)用地比例僅為0.2%)，土地的開(kāi)發(fā)程度明顯不足，各項(xiàng)指標(biāo)均遠(yuǎn)低于平均值.

(a)2010年 (b)2018年

3.2 生態(tài)系統(tǒng)服務(wù)供需平衡分析

2010年和2018年生態(tài)系統(tǒng)服務(wù)熱點(diǎn)分析結(jié)果，如圖4所示.由圖4可知：2010年和2018年分別有9個(gè)和11個(gè)大斑塊為供需失衡的關(guān)鍵區(qū)域；2010年有4個(gè)斑塊為供需盈余區(qū)，主要分布于池州、安慶、銅陵3市交界處的區(qū)/縣，總面積為13 023.56 km2，2010年有5個(gè)斑塊為供需赤字區(qū)，分布于上海、南京、杭州等經(jīng)濟(jì)發(fā)達(dá)城市的主城區(qū)，總面積為20 082.03 km2；2018年有5個(gè)斑塊為供需盈余區(qū)，總面積為19 921.69 km2，比2010年增加6 898.13 km2，增加的斑塊主要位于杭州市桐廬縣、杭州市建德市和寧波市寧海縣，桐廬縣2018年森林覆蓋率達(dá)75.31%，林地面積顯著增加，總供給價(jià)值增加48.15元·km2，建德市2010年至2018年城鎮(zhèn)經(jīng)濟(jì)發(fā)展分散，需求增長(zhǎng)不顯著，生態(tài)服務(wù)供給富余，寧?？h近8年間生態(tài)治理成效顯著，生態(tài)系統(tǒng)服務(wù)供給充足，2018年有6個(gè)斑塊為供需赤字區(qū)，總面積為24 079.65 km2，比2010年增加3 997.62 km2，增加的斑塊主要位于南京市江寧區(qū)、鎮(zhèn)江市句容市、蘇州市張家港市、無(wú)錫市江陰市，自2010年后，江寧區(qū)進(jìn)入工業(yè)化后期階段，區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展，生態(tài)需求顯著提升，句容市大力推進(jìn)城鎮(zhèn)化建設(shè)，研究年份內(nèi)建設(shè)用地比例增加12.5%，江陰市和張家港市位列2018年中國(guó)縣域經(jīng)濟(jì)百?gòu)?qiáng)排行榜前三，產(chǎn)業(yè)結(jié)構(gòu)以第二產(chǎn)業(yè)為主導(dǎo)，工業(yè)發(fā)達(dá)，部分生態(tài)用地被占用，區(qū)域生態(tài)系統(tǒng)服務(wù)供需失衡.

(a)2010年 (b)2018年

3.3 生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)路徑與價(jià)值流轉(zhuǎn)量

3.3.1 生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)路徑 將生態(tài)系統(tǒng)服務(wù)供需盈余和供需赤字的斑塊分別判定為生態(tài)源地和生態(tài)需求地，考慮到部分斑塊區(qū)/縣數(shù)不止一個(gè)，且少數(shù)相鄰區(qū)/縣集中連片，文中生態(tài)源地和生態(tài)需求地的名稱用其中面積最大且聚集程度最高的區(qū)/縣名稱表示.2010年，有4個(gè)主要生態(tài)源地(安慶市太湖縣、池州市貴池區(qū)、銅陵市樅陽(yáng)縣和常州市溧陽(yáng)市)，5個(gè)主要生態(tài)需求地(南京市六合區(qū)、湖州市德清縣、嘉興市桐鄉(xiāng)市、杭州市蕭山區(qū)和上海全區(qū)/縣).2018年，有5個(gè)主要生態(tài)源地(安慶市宿松縣、池州市貴池區(qū)、銅陵市樅陽(yáng)縣、杭州市建德市和寧波市寧?？h)，6個(gè)主要生態(tài)需求地(杭州市蕭山區(qū)、嘉興市海寧市、嘉興市桐鄉(xiāng)市、無(wú)錫市江陰市、南京市江寧區(qū)及上海全區(qū)/縣).由ArcGIS平臺(tái)的成本路徑工具生成生態(tài)源地和生態(tài)需求地之間的多條潛在生態(tài)廊道.2010年和2018年生態(tài)廊道、生態(tài)源地和生態(tài)需求地的空間分布，如圖5所示.

(a)2010年 (b)2018年

根據(jù)重力矩陣可得2010年和2018年生態(tài)斑塊的相互作用矩陣，如表2，3所示.表2，3中：F為各斑塊之間的作用力，無(wú)量綱.由表2可知：2010年的5個(gè)生態(tài)需求地均與貴池區(qū)的相互作用力最強(qiáng)，表明貴池區(qū)與各生態(tài)需求地進(jìn)行生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)的可能性最大，進(jìn)而提取出5條重要生態(tài)廊道.由表3可知：2018年的6個(gè)生態(tài)需求地均與建德市的相互作用力最強(qiáng)，由此提取出6條重要生態(tài)廊道作為生態(tài)系統(tǒng)服務(wù)從生態(tài)源地到生態(tài)需求地的流轉(zhuǎn)路徑.重要生態(tài)廊道將生態(tài)需求地與生態(tài)源地進(jìn)行了有效連通，以此作為生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)的最佳路徑，建立起自然生態(tài)系統(tǒng)與城鎮(zhèn)需求空間的聯(lián)系紐帶[1].

表2 2010年生態(tài)斑塊的相互作用矩陣

表3 2018年生態(tài)斑塊的相互作用矩陣

3.3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量 根據(jù)斷裂點(diǎn)公式和場(chǎng)強(qiáng)模型，可得2010年和2018年生態(tài)源地向生態(tài)需求地流轉(zhuǎn)的生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量及相關(guān)參數(shù)，如表4，5所示.

表4 2010年生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量及相關(guān)參數(shù)

由表4可知：由于貴池區(qū)到桐鄉(xiāng)市的距離最近，流轉(zhuǎn)強(qiáng)度最大，兩地之間的生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量最多(73.82億元)；貴池區(qū)到六合區(qū)輻射面積遠(yuǎn)小于其他需求地，因而生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量最少(6.74億元)；貴池區(qū)到德清縣、蕭山區(qū)、上海全區(qū)/縣的流轉(zhuǎn)量分別為23.15，10.26，13.13億元；2010年生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)總量為127.1億元.

由表5可知:建德市到桐鄉(xiāng)市具有較大的輻射面積及流轉(zhuǎn)強(qiáng)度，兩地之間的生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量最多(92.15億元)；建德市到上海全區(qū)/縣的兩斑塊質(zhì)心之間的距離最長(zhǎng)，因而流轉(zhuǎn)量最少(15.11億元)，建德縣到蕭山區(qū)、海寧市、江陰市和江寧區(qū)的流轉(zhuǎn)量分別為74.60，40.45，28.20，22.09億元；2018年生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)總量為272.6億元，2018年的流轉(zhuǎn)量大約是2010年的兩倍，說(shuō)明隨著經(jīng)濟(jì)社會(huì)的發(fā)展，城市生態(tài)需求日益增多，不同區(qū)域間生態(tài)系統(tǒng)服務(wù)供需匹配的差異增大.

表5 2018年生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量及相關(guān)參數(shù)

4 討論與結(jié)論

2018年供需失衡區(qū)面積的增加導(dǎo)致生態(tài)服務(wù)價(jià)值流轉(zhuǎn)量的增多，以及生態(tài)源地和生態(tài)需求地空間分布格局的改變，其關(guān)鍵原因在于區(qū)域經(jīng)濟(jì)水平的提升使生態(tài)需求指數(shù)顯著增長(zhǎng)，第二、三產(chǎn)業(yè)是區(qū)/縣經(jīng)濟(jì)發(fā)展的主要?jiǎng)恿?，為滿足人們對(duì)更高生態(tài)質(zhì)量的需求提供各種原材料和消費(fèi)品，生態(tài)系統(tǒng)服務(wù)供不應(yīng)求.此外，城鎮(zhèn)化進(jìn)程改變?cè)型恋乩媒Y(jié)構(gòu)，建設(shè)用地?cái)U(kuò)張占用部分生態(tài)用地，自然生態(tài)系統(tǒng)受到人為活動(dòng)干擾程度加深，導(dǎo)致生態(tài)系統(tǒng)服務(wù)功能下降，供需空間異質(zhì)性增強(qiáng).結(jié)合長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)供需關(guān)系，通過(guò)構(gòu)建生態(tài)源地和生態(tài)需求地之間的重要生態(tài)廊道確定生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)路徑，并計(jì)算生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量，可為相關(guān)部門(mén)制定生態(tài)補(bǔ)償?shù)瓤鐓^(qū)域環(huán)境協(xié)調(diào)治理政策提供一定的科學(xué)依據(jù)；根據(jù)生態(tài)源地流轉(zhuǎn)到生態(tài)需求地的生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量比例，地方政府可以進(jìn)行生態(tài)服務(wù)價(jià)值補(bǔ)償出資比例的安排，從而保障地區(qū)公平的發(fā)展權(quán).

基于價(jià)值量、熱點(diǎn)分析、MCR模型和距離衰減等多種模型和方法，量化分析長(zhǎng)三角城市群2010年和2018年的生態(tài)系統(tǒng)服務(wù)供給和需求，揭示其供需空間格局特征，并確定生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)路徑和生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)量，可得以下3個(gè)結(jié)論.

1)長(zhǎng)三角城市群2010年和2018年生態(tài)系統(tǒng)服務(wù)供給變化微小，而生態(tài)系統(tǒng)服務(wù)需求顯著增長(zhǎng)；在空間分布上，生態(tài)系統(tǒng)服務(wù)供給呈南高北低的分布形態(tài)；需求大致呈現(xiàn)東高西低的分布形態(tài)，沿海、沿江明顯高于內(nèi)陸地區(qū).

2)供需盈余區(qū)主要分布于皖南山地丘陵，供需赤字區(qū)主要分布于上海、南京、杭州等經(jīng)濟(jì)發(fā)達(dá)城市的主城區(qū)，研究年份內(nèi)，供需平衡區(qū)面積下降，供需失衡區(qū)面積增加，經(jīng)濟(jì)發(fā)展水平與土地利用方式的變化導(dǎo)致供需空間異質(zhì)性增強(qiáng).

3)將生態(tài)源地和生態(tài)需求地之間的重要生態(tài)廊道作為生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)路徑，可得2010年和2018年長(zhǎng)三角城市群生態(tài)系統(tǒng)服務(wù)價(jià)值流轉(zhuǎn)總量為127.1，272.6億元.

生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)是一個(gè)動(dòng)態(tài)變化的過(guò)程，不同廊道間傳遞的服務(wù)流量各有差異.文中構(gòu)建的生態(tài)系統(tǒng)服務(wù)空間流轉(zhuǎn)模型有一定的不足：首先，評(píng)估模型中用于反映自然影響因素的修正因子的選取多為經(jīng)驗(yàn)值，隨著對(duì)生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)規(guī)律認(rèn)識(shí)的加深，需進(jìn)一步提高其準(zhǔn)確性；其次，在生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)過(guò)程中僅考慮了自然因素，忽略了人為活動(dòng)的影響.今后的研究還需借助多種技術(shù)手段與數(shù)據(jù)信息，對(duì)生態(tài)系統(tǒng)服務(wù)流轉(zhuǎn)的路徑方向、流量大小、衰減特征等進(jìn)行具體、準(zhǔn)確的描述.