范 強(qiáng)，史 月，王戲游

（1.遼寧工程技術(shù)大學(xué) 測(cè)繪與地理科學(xué)學(xué)院，遼寧 阜新 123000；2.廣州南方測(cè)繪科技股份有限公司 導(dǎo)航事業(yè)部，廣東 廣州 510663）

0 引言

城市建成區(qū)用地?cái)U(kuò)張為城市發(fā)展帶來(lái)城市地表溫度變化[1]、PM2.5 污染程度變化[2]和區(qū)域生態(tài)安全性變化[3]等多種城市問(wèn)題.現(xiàn)今國(guó)內(nèi)外學(xué)者針對(duì)城市擴(kuò)張的研究逐漸增多，其中，中國(guó)學(xué)者關(guān)于城市擴(kuò)張?zhí)卣餮芯考胺治鲆灾袊?guó)、珠三角、京津冀、江西省和武漢市等國(guó)家、省和大中發(fā)達(dá)城市為主[4-7]；國(guó)外則以亞利桑那州、美國(guó)和印度等西方南亞國(guó)家為主[8-11].國(guó)內(nèi)外對(duì)于大中城市或城市群以及發(fā)達(dá)地區(qū)的建成區(qū)用地研究逐漸升溫，但長(zhǎng)時(shí)間序列的資源枯竭型城市建成區(qū)用地?cái)U(kuò)張研究較少.資源枯竭型城市作為獨(dú)特的城市類(lèi)型，其建成區(qū)用地?cái)U(kuò)張導(dǎo)致的城市問(wèn)題尤為嚴(yán)重，而且資源枯竭型城市轉(zhuǎn)型是國(guó)家可持續(xù)發(fā)展的重要環(huán)節(jié)，因此，對(duì)資源枯竭型城市建成區(qū)用地?cái)U(kuò)張?zhí)卣鬟M(jìn)行量化分析、合理解讀資源枯竭型城市的擴(kuò)張?zhí)卣?，?duì)資源枯竭型城市規(guī)劃及決策等有著重大意義.

綜上所述，本文利用Landsat TM/OLI 遙感影像，綜合運(yùn)用RS、GIS 和景觀學(xué)相結(jié)合的方式，分析1999—2019 年阜新市建成區(qū)用地?cái)U(kuò)張的時(shí)空特征，對(duì)資源枯竭型城市的發(fā)展及轉(zhuǎn)型有著重要的參考意義.

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

1.1 研究區(qū)概況

由圖1 可知，阜新市位于遼寧省西北部，東經(jīng)121°01′～122°56′，北緯41°41′～42°56′.下設(shè)細(xì)河區(qū)、太平區(qū)、新邱區(qū)、清河門(mén)區(qū)、海州區(qū)、彰武縣和阜新蒙古族自治縣（下文簡(jiǎn)稱(chēng)阜蒙縣）7 個(gè)行政區(qū).礦產(chǎn)資源十分豐富，是中國(guó)重要的礦業(yè)開(kāi)發(fā)基地之一，除礦產(chǎn)資源外，煤層氣及地?zé)岬荣Y源也相繼被開(kāi)發(fā)使用.自1897 年以來(lái)，經(jīng)過(guò)120 多年的開(kāi)采，地下煤礦幾近枯竭，而且造成水污染、大氣污染、土地資源嚴(yán)重破壞等一系列生態(tài)環(huán)境問(wèn)題.2001 年被國(guó)務(wù)院正式認(rèn)定為全國(guó)第一個(gè)資源枯竭型城市.

圖1 研究區(qū)位置Fig.1 location of study area

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

采集1999 年、2004 年、2009 年、2014 年和2019 年5 個(gè)年份的Landsat TM/OLI 遙感影像（來(lái)源于美國(guó)地質(zhì)勘探局網(wǎng)站：http://glovis.usgs.gov/），利用隨機(jī)森林分類(lèi)方法對(duì)5 期遙感影像進(jìn)行解譯.最終提取得到建成區(qū)用地的Kappa 系數(shù)為0.85～0.92，符合研究精度要求.

1.3 用地?cái)U(kuò)張時(shí)空特征

（1）城市擴(kuò)張強(qiáng)度指數(shù)

城市擴(kuò)張強(qiáng)度指數(shù)（UEI）表示同一城市在不同時(shí)間段的擴(kuò)張強(qiáng)度，對(duì)于研究城市擴(kuò)張強(qiáng)度及擴(kuò)張?zhí)卣骶哂辛己玫谋碚饕饬x，公式[12]為

式中，UEIn為第n年城市擴(kuò)張強(qiáng)度指數(shù)；分別為t2和t1時(shí)間段的城市建成區(qū)面積，m2；t為t2和t1的間隔年份.

單一時(shí)間段的城市擴(kuò)張強(qiáng)度指數(shù)可以較好地反映出研究區(qū)域各縣區(qū)建成區(qū)用地的擴(kuò)張強(qiáng)度，但是缺乏縣區(qū)間的橫向?qū)Ρ龋蕵?gòu)建城市擴(kuò)張強(qiáng)度均指數(shù)（AUEI）為

式中，N為年份時(shí)間間隔個(gè)數(shù).

（2）Voronoi 圖變異系數(shù)

離散的點(diǎn)集在空間上一般分為均勻、隨機(jī)和聚集3 種空間布局分布方式.離散點(diǎn)集的Voronoi 圖可以較為直觀地描述區(qū)域內(nèi)建成區(qū)用地的分布情況.本文引用Voronoi 圖的變異系數(shù)來(lái)描述城市建成區(qū)用地的空間分布情況.DUYCKAERTS[13]等認(rèn)為Voronoi 圖的變異系數(shù)小于33%時(shí)，點(diǎn)集呈均勻分布；變異系數(shù)在33%～64%時(shí)，點(diǎn)集呈隨機(jī)分布；變異系數(shù)大于64%時(shí)，點(diǎn)集呈聚集分布.本文采用的離散點(diǎn)集是由建成區(qū)用地的質(zhì)心集合而成.

式中，CV為Voronoi 圖的變異系數(shù)；Ai為第i個(gè)Voronoi 多邊形的面積，m2；n為Voronoi 多邊形的總數(shù).

（3）景觀格局指數(shù)

隨著時(shí)間的演變，城市建成區(qū)用地逐漸擴(kuò)張，從而導(dǎo)致景觀格局隨之變化.利用景觀學(xué)中指數(shù)分析方式進(jìn)行城市建成區(qū)情況分析，可以較好地反映建成區(qū)用地的演變規(guī)律.本文選用最大斑塊指數(shù)、斑塊類(lèi)型比例、景觀形狀指數(shù)、結(jié)合度指數(shù)和分維數(shù)共5 個(gè)景觀格局指標(biāo).

2 城市建成區(qū)用地?cái)U(kuò)張時(shí)空格局

2.1 城市擴(kuò)張強(qiáng)度時(shí)空分析

將1999—2004 年、2004—2009 年、2009—2014年和2014—2019 年4 期阜新市各縣區(qū)建成區(qū)用地?cái)U(kuò)張強(qiáng)度進(jìn)行量化分析，見(jiàn)圖2.

圖2 阜新市各縣區(qū)用地?cái)U(kuò)張強(qiáng)度Fig.2 UEI of Fuxin each county area

1999—2019 年，阜新市各縣區(qū)的建成區(qū)面積不斷增加，不同時(shí)間段內(nèi)各縣區(qū)擴(kuò)張強(qiáng)度存在較大差異.1999—2019 年，彰武縣的UEI指數(shù)在阜新市各縣區(qū)內(nèi)維持在最高數(shù)值.1999—2004 年，細(xì)河區(qū)的UEI指數(shù)位于彰武縣之后，阜蒙縣次之，海州區(qū)、清河門(mén)區(qū)、太平區(qū)和新邱區(qū)UEI指數(shù)均處于0.04%以下.2004—2009 年，阜蒙縣UEI指數(shù)超過(guò)細(xì)河區(qū)，細(xì)河區(qū)為0.0472%，其余海州區(qū)、清河門(mén)區(qū)、太平區(qū)和新邱區(qū)四區(qū)依然維持在0.04%以下.2009—2014 年，阜蒙縣和細(xì)河區(qū)UEI指數(shù)為0.092%和0.043 6%，清河門(mén)區(qū)首次達(dá)到0.04%，海州區(qū)、太平區(qū)和新邱區(qū)依然處于0.04%以下.2014—2019 年，海州區(qū)、新邱區(qū)和太平區(qū)的擴(kuò)張強(qiáng)度維持在0.04%以下，其余各縣市的擴(kuò)張強(qiáng)度處于0.04%之上.另外通過(guò)AUEI指數(shù)橫向?qū)Ρ雀沸率懈骺h區(qū)建成區(qū)用地?cái)U(kuò)張強(qiáng)度，可以得知，20 年間，阜新市各縣區(qū)的建成區(qū)用地?cái)U(kuò)張強(qiáng)度均指數(shù)以彰武縣為首；阜蒙縣和細(xì)河區(qū)在其之后；海州區(qū)、清河門(mén)區(qū)、太平區(qū)和新邱區(qū)的AUEI指數(shù)較小.

2.2 城市建成區(qū)用地分布模式

根據(jù)阜新市各縣區(qū)建成區(qū)用地質(zhì)心構(gòu)成的點(diǎn)集生成Voronoi 圖，計(jì)算變異系數(shù)Cv值見(jiàn)圖3.

圖3 各縣區(qū)變異系數(shù)Fig.3 variation index of each county district

Voronoi 圖的變異系數(shù)Cv可以較為直觀地表示城市建成區(qū)用地分布情況，由圖3 可以看出阜蒙縣在1999—2019 年，呈現(xiàn)均勻分布模式，變異系數(shù)值均小于33%；海州區(qū)的變異系數(shù)在各年趨近于64%，總體為隨機(jī)分布模式；清河門(mén)區(qū)在1999—2004 年從聚集模式轉(zhuǎn)為隨機(jī)模式后，在2004—2019 年保持聚集模式分布；太平區(qū)的變異系數(shù)在1999—2019 年逐年遞減，但總體分布模式為聚集模式；細(xì)河區(qū)的變異系數(shù)維持在33%以下，屬均勻分布模式；新邱區(qū)1999—2019 年的變異系數(shù)均大于64%，保持在聚集模式；而彰武縣的整體分布模式是由均勻分布轉(zhuǎn)為隨機(jī)分布，變異系數(shù)數(shù)值呈“V”型分布.

2.3 景觀生態(tài)效應(yīng)

由圖4 可以看出，1999—2019 年阜新市各縣區(qū)的景觀格局變化趨勢(shì).海州區(qū)、太平區(qū)、新邱區(qū)和清河門(mén)區(qū)的LPI指數(shù)在1999—2019 年保持穩(wěn)定態(tài)勢(shì)，阜蒙縣和彰武縣呈現(xiàn)下降趨勢(shì)，細(xì)河區(qū)則呈現(xiàn)上升趨勢(shì)，表示海州區(qū)、太平區(qū)、新邱區(qū)和清河門(mén)區(qū)的建成區(qū)用地優(yōu)勢(shì)斑塊地位幾乎沒(méi)有改變；阜蒙縣和彰武縣的優(yōu)勢(shì)斑塊地位有所下降，而細(xì)河區(qū)則出現(xiàn)上升態(tài)勢(shì).所有縣區(qū)的PLAND指數(shù)呈現(xiàn)逐年增長(zhǎng)趨勢(shì)，其中細(xì)河區(qū)和清河門(mén)區(qū)增長(zhǎng)的線條更為陡峭，說(shuō)明20 年間，城市建成區(qū)的面積相對(duì)自身來(lái)說(shuō)，增長(zhǎng)的幅度較大.海州區(qū)、細(xì)河區(qū)、新邱區(qū)、太平區(qū)和清河門(mén)區(qū)的LSI指數(shù)變化較為平緩，1999—2019年，形狀幾乎沒(méi)有太大的改變；而彰武縣和阜蒙縣呈現(xiàn)的波動(dòng)性十分巨大，二者的數(shù)值都有增大后穩(wěn)定的趨勢(shì)，彰武縣變化幅度較大.海州區(qū)、太平區(qū)和新邱區(qū)在COHESION指數(shù)上表現(xiàn)為變化較為平緩，但是數(shù)據(jù)值較大，說(shuō)明三區(qū)內(nèi)部建成區(qū)用地的連通性較好；細(xì)河區(qū)建成區(qū)用地的內(nèi)部連通性也在不斷地增大；隨著時(shí)間的變化，阜蒙縣、彰武縣和清河門(mén)區(qū)的建成區(qū)用地內(nèi)部連通性較差，其中以彰武縣變化最為明顯.PAFRAC指數(shù)作為一個(gè)衡量建成區(qū)用地受人為干擾程度的指數(shù)，在1999—2019 年表示為清河門(mén)區(qū)和海州區(qū)受人為干擾程度振蕩遞減；其余5 縣區(qū)，受人為干擾程度總體呈上升趨勢(shì)，其中彰武縣在2009 年以后受人為干擾程度最大.

圖4 各縣區(qū)景觀格局指數(shù)變化Fig.4 landscape pattern of each county area changes exponentially

3 結(jié)論

阜新市作為中國(guó)第一個(gè)資源枯竭型城市，是資源枯竭型城市的典型代表，基于1999—2019 年城市建成區(qū)用地?cái)?shù)據(jù)，綜合運(yùn)用城市擴(kuò)張強(qiáng)度指數(shù)、擴(kuò)張強(qiáng)度均指數(shù)、Voronoi 圖變異系數(shù)Cv和景觀格局指數(shù)對(duì)城市擴(kuò)張強(qiáng)度、用地分布模式及景觀生態(tài)效應(yīng)進(jìn)行分析，得出結(jié)論：

（1）1999—2019 年，隨著阜新市各縣區(qū)建成區(qū)用地面積逐漸增加，建成區(qū)范圍逐漸擴(kuò)大，城市擴(kuò)張強(qiáng)度也隨之發(fā)生變化.從全市視角來(lái)看，東部城市擴(kuò)張強(qiáng)度高于西部，整體的擴(kuò)張中心由西部向東部移動(dòng).其中，彰武縣的擴(kuò)張強(qiáng)度變化程度最大，呈高速擴(kuò)張趨勢(shì).2009 年以后，由于資源枯竭型城市的轉(zhuǎn)型工作不斷推進(jìn)，致使除細(xì)河區(qū)外的所有縣區(qū)城市擴(kuò)張強(qiáng)度變大；而細(xì)河區(qū)的城市擴(kuò)張強(qiáng)度在20 年間，基本保持相對(duì)穩(wěn)定的趨勢(shì).

（2）由于阜新市各縣區(qū)建成區(qū)建設(shè)情況不均.各縣區(qū)的用地分布模式以聚集、隨機(jī)和均勻3 種分布模式共存，其中阜蒙縣和細(xì)河區(qū)呈均勻分布模式；太平區(qū)和新邱區(qū)呈聚集分布模式；海州區(qū)呈隨機(jī)分布模式；而彰武縣和清河門(mén)區(qū)則由于建成區(qū)用地的擴(kuò)張，用地分布模式發(fā)生了轉(zhuǎn)換.

（3）20 年間，阜新市各縣區(qū)建成區(qū)用地面積不斷增大，人類(lèi)與自然相互作用程度也不斷增大，景觀格局隨之發(fā)生較大的變化.其中彰武縣、阜蒙縣和細(xì)河區(qū)的優(yōu)勢(shì)斑塊地位發(fā)生了明顯改變；所有縣區(qū)的建成區(qū)面積均呈上升趨勢(shì)，導(dǎo)致建成區(qū)用地內(nèi)部連通性演變趨勢(shì)也隨之變化；從建成區(qū)用地形狀上來(lái)說(shuō)，彰武縣和阜蒙縣波動(dòng)情況較為明顯，且彰武縣受人為干擾程度較大.

（4）從縣區(qū)尺度上分析阜新市建成區(qū)用地時(shí)空擴(kuò)張?zhí)卣?、用地分布及景觀格局，對(duì)深刻了解阜新市建成區(qū)用地?cái)U(kuò)張?zhí)卣骶哂兄匾饬x.但是資源枯竭型城市的建成區(qū)用地?cái)U(kuò)張是一個(gè)集經(jīng)濟(jì)、社會(huì)、政策和環(huán)境等多種元素影響的過(guò)程，本文以阜新市為例，對(duì)于資源枯竭型城市的用地?cái)U(kuò)張分析的不夠全面，在后續(xù)的研究中將進(jìn)一步對(duì)不同類(lèi)型的資源枯竭型城市進(jìn)行分析，從而全面地分析資源枯竭型城市建成區(qū)用地演變規(guī)律.