楊璟旎，方 然

（武漢大學資源與環(huán)境科學學院，湖北武漢 430079）

我國進入快速城市化發(fā)展階段以來，城市用地盲目擴張，不僅造成土地資源浪費，而且造成城市及其周邊地區(qū)生態(tài)環(huán)境的破壞，如何科學合理地確定城鎮(zhèn)建設和保護范圍，成為土地學科研究的重點。我國在“十一五”規(guī)劃中提出了主體功能區(qū)劃，將區(qū)域空間劃分為重點開發(fā)區(qū)、優(yōu)化開發(fā)區(qū)、限制開發(fā)區(qū)和禁止開發(fā)區(qū);在城鎮(zhèn)區(qū)域總體規(guī)劃中，將土地空間劃分為已建區(qū)、適建區(qū)、限建區(qū)和禁建區(qū)［1］。而遙感技術宏觀、快速、準確獲取數(shù)據的優(yōu)勢，以及越來越高的分辨率和越來越短的重復觀測周期，正好適應了城市變化快速、結構復雜的特點［2］。因此，利用遙感與GIS技術研究城鎮(zhèn)發(fā)展適宜性時空演變特征，對提高土地利用綜合效益具有重要實用價值，其結果能為土地規(guī)劃、城市規(guī)劃等提供依據，尤其是在建設用地的布局與開發(fā)方面具有重要的現(xiàn)實意義。

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

1.1 研究區(qū)概況 武漢是湖北省省會，地處我國腹地中心，江漢平原東部(113°41'～115°05'E、29°58'～31°22'N)。世界第三大河長江及其最大支流漢水橫貫市境中央，將武漢一分為三，形成了武昌、漢口、漢陽3鎮(zhèn)隔江鼎立的格局。武漢市幅員面積8 494.41 km2，下轄7個中心城區(qū)和6個遠城區(qū)，2013年常住人口1 022萬。武漢是中部6省唯一的副省級城市，也是國務院批準的中部地區(qū)中心城市，按照城市常住人口、GDP總量、建成區(qū)面積等指標，武漢均為中部第一、全國前十的特大城市，也是世界各大城市中人均擁有水量最多的城市。2014年武漢市實現(xiàn)地區(qū)生產總值10 069.48億元，位列全國所有城市的第八位［3］。

1.2 數(shù)據來源 該研究數(shù)據主要來源于遙感影像數(shù)據、地形數(shù)據以及武漢市歷年土地利用現(xiàn)狀數(shù)據。其中遙感影像數(shù)據分別選取2001、2007與2013年3期的30 m×30 m Landsat7 ETM數(shù)據，地形數(shù)據選取武漢市DEM數(shù)據與坡度數(shù)據，同時以武漢市歷年土地利用現(xiàn)狀數(shù)據以及各級路網數(shù)據、城市軌道交通網絡數(shù)據作為參考。3期Landsat7 ETM影像數(shù)據采用5、4、3三波段進行假彩色合成，均轉換至西安80投影坐標系。在Envi4.8中對3期遙感影像進行空間校正、拼接、裁剪，采用統(tǒng)一的解譯標志(Region of Interest)進行分類，建立2001、2007、2013年3個時期的土地利用空間、屬性數(shù)據庫。

1.3 評價指標體系的構建

1.3.1 評價因子的選取。根據武漢市的自然環(huán)境特點以及當前的社會經濟發(fā)展狀況，結合近年來國內外學者在城鎮(zhèn)發(fā)展適宜性評價領域的主要研究成果，將影響武漢市城鎮(zhèn)發(fā)展的因素分為自然環(huán)境因素、區(qū)位因素以及生態(tài)環(huán)境因素3大類，并在每一大類中細分若干因子，共同構建武漢市城鎮(zhèn)發(fā)展適宜性評價指標體系。

(1)自然環(huán)境因素。主要選取了坡度和高程2大因素。坡度和高程是土地資源固有的2個重要環(huán)境因子，對土地利用和土地生產力有直接作用。坡度和高程基本上決定了土地利用方向和利用方式，從而影響到經濟效益和發(fā)展?jié)摿ΑＦ露鹊拇笮こ探ㄔO的投入、施工難度和城市供水等產生重要影響。因此，應該盡量避免在較大坡度區(qū)域進行城市建設。

(2)區(qū)位因素。區(qū)位適宜度反映新布局建設用地的便利程度，新增建設用地具有向交通線路以及現(xiàn)狀城鎮(zhèn)等輻射效應性強的地段集聚的態(tài)勢。城鎮(zhèn)的擴張與現(xiàn)狀城鎮(zhèn)以及交通網絡的輻射效應直接相關，因此社會經濟因素選取與現(xiàn)狀城鎮(zhèn)的距離、與城市軌道交通線路的距離、與國道的距離、與省道的距離以及與市內主干道距離5大因素。

(3)生態(tài)環(huán)境因素。生態(tài)適宜度反映生態(tài)用地對城鎮(zhèn)建設用地的外部性環(huán)境影響程度，建設用地靠近生態(tài)用地能獲得更多的生態(tài)效益。評價中選取了與湖泊、水庫、林地、江河的距離等4個因子。但是，值得注意的是，生態(tài)因子對建設用地的布局也會起到限制性作用，布局新增建設用地不能損害生態(tài)效益，所以在水體保護區(qū)和自然與人文景觀保護區(qū)所在的區(qū)域，在分級指標量化時將這些區(qū)域賦空值(NoData)，即作為城鎮(zhèn)發(fā)展用地適宜性的禁區(qū)而存在，以避免分析過程中出現(xiàn)因子疊加后這些區(qū)域適宜或部分適宜的現(xiàn)象，保證評價結果的客觀現(xiàn)實性［4］。

1.3.2 權重的確定。城鎮(zhèn)發(fā)展適宜性評價各因素因子的權重通過層次分析法(AHP)結合專家打分法(Delphi)確立。層次分析法是一種多目標決策方法，基于系統(tǒng)的層次性原理，將復雜的問題分成若干層次，逐步分解比較，把人的主觀判斷用數(shù)量的形式來表達和處理［5-7］。專家打分法是指彼此相互獨立的多名專家通過匿名的方式對同一目標給出意見，經過多輪的意見反饋與修正后綜合成基本一致的意見與決策，并以此作為最終的評價結果［8-11］。

根據各因素因子對于城鎮(zhèn)發(fā)展適宜性的影響機理，對各因子在不同距離范圍的輻射作用進行評分，各因子層適宜性分值介于0～100，分值越高，單因子在該輻射范圍內的城鎮(zhèn)發(fā)展適宜度越高。評價指標體系與各因子在不同范圍內的得分及權重如表1所示。

1.4 模型構建與評價 在ArcGIS9.3中通過武漢市數(shù)字高程模型(DEM)獲取高程數(shù)據，通過空間分析模塊中的坡度提取功能獲得坡度數(shù)據;將區(qū)位因素與生態(tài)環(huán)境因素中的各因子按指標體系中的距離進行緩沖區(qū)分析，通過重分類功能將每一因子層的不同緩沖區(qū)域賦值，分別得到11個指標的單因子適宜度評價圖層。

該研究使用多因子加權疊加的方式得到最終的評價結果。以各單項因子的分析為基礎，結合權重值，計算各單項因子的適宜性得分，多因子加權疊加匯總確定城鎮(zhèn)發(fā)展用地適宜性級別，其綜合適宜性評價公式為［12-14］:

式中，Sij為第ij個格網的綜合適宜性分值;gij(m)表示第m個評價因子在第ij個格網的適宜性分值;wm為表示第m個評價因子的權重。

2 結果與分析

在ArcGIS9.3中分別對11個因子層進行緩沖區(qū)分析，根據指標體系中各評價因子所處區(qū)域范圍的不同分別賦予0～100的分值，生成11個單因子評價層，利用柵格計算器(raster calculator)加權疊加各單因子適宜性分級圖層，得到綜合評價結果，并將其按照適宜性分值劃分為4個等級:高適宜區(qū)、中適宜區(qū)、低適宜區(qū)、不適宜區(qū)。最后將生態(tài)與水體保護區(qū)嵌入綜合評價圖層作為禁止建設的區(qū)域，得到2001、2007、2013年3期武漢市城鎮(zhèn)發(fā)展用地適宜性評價的最終結果(圖1)以及不同級別城鎮(zhèn)發(fā)展適宜區(qū)面積統(tǒng)計表(表2)。

表1 城鎮(zhèn)發(fā)展適宜性評價因素因子體系

表2 2001、2007、2013年城鎮(zhèn)發(fā)展適宜性分級面積統(tǒng)計

2.1 時間演變特征分析 從表2可以看出，伴隨著近十年來武漢市的城市建設和擴張運動，2001～2013年武漢市城鎮(zhèn)發(fā)展高、中適宜區(qū)土地面積總體呈遞增趨勢，其中高適宜區(qū)的面積由2001年的117 087.84 hm2上升到2013年的204 999.75 hm2，增加了75%。新一輪的城市化以及城市人口的爆炸式增長帶來了一系列用地需求的增長，特別是帶來了交通便利程度的增加，因而低適宜區(qū)以及不適宜區(qū)總體呈下降趨勢，其中低適宜區(qū)由2001年的307 028.79 hm2下降到2013年的214 777.98 hm2，下降了30%。同時，2001～2013年中適宜發(fā)展地區(qū)占總面積的比重并未發(fā)生大的改變，僅從23.79%變化到25.52%;而所占比例變化最為明顯的就是高適宜區(qū)，由15.09%增長到26.83%;低適宜區(qū)所占比例則由39.57%大幅下降到28.11%。2001～2013年武漢市城鎮(zhèn)土地發(fā)展適宜區(qū)各類土地面積的變化側面反映了城市發(fā)展所帶來的自然環(huán)境、生態(tài)環(huán)境特別是某些地區(qū)區(qū)位條件的改善，因而使得城鎮(zhèn)總體適宜發(fā)展的土地面積增加［15］。

2.2 空間演變特征分析 從圖1可以看出，2001年武漢市發(fā)展高適宜區(qū)主要集中在中心城區(qū)以及黃陂和新洲等少數(shù)遠城區(qū)地區(qū)。伴隨著接下來幾年中東湖高新技術開發(fā)區(qū)以及吳家山經濟開發(fā)區(qū)的成立、武漢經濟技術開發(fā)區(qū)的發(fā)展，到了2007和2013年，武漢市發(fā)展高適宜區(qū)擴展到了蔡甸區(qū)、漢南區(qū)、東西湖區(qū)以及江夏區(qū)的部分地區(qū)，總體分布呈擴張趨勢。而由于城市邊緣的生態(tài)環(huán)境和自然環(huán)境的低水平，特別是交通和區(qū)位條件的低水平，武漢市城市發(fā)展不適宜區(qū)主要分布在新洲區(qū)的東端、蔡甸區(qū)的西端、江夏區(qū)的南端以及最黃陂區(qū)靠北段。

3 結論與討論

該研究以武漢市為例，基于2001、2007與2013年3時期遙感數(shù)據，從城鎮(zhèn)發(fā)展適宜性時空演變特征角度進行分析，旨在為城市合理發(fā)展提供科學的依據。研究過程中，選擇了城鎮(zhèn)發(fā)展適宜性評價的多個指標，利用GIS技術進行疊加分析，較客觀地實現(xiàn)了城鎮(zhèn)發(fā)展的適宜性評價。結果表明，2001～2013年武漢市城市化進程發(fā)展十分迅速。由于武漢市特殊的地理位置以及多水域的城市城市特征，城鎮(zhèn)發(fā)展高適宜性地區(qū)主要集中在中心城區(qū)以及大型水域附近。同時伴隨著武漢大舉造城，蔡甸區(qū)等幾個遠城區(qū)的城鎮(zhèn)發(fā)展高適宜性地區(qū)面積有所增加。

在武漢市的城市建設用地規(guī)劃過程中應合理考慮該地區(qū)的生態(tài)適宜性，規(guī)劃布局應順應生態(tài)系統(tǒng)的自然規(guī)律，用地應盡量選擇適宜度高的區(qū)域，這樣才能最大程度地保證城市優(yōu)良的自然生態(tài)環(huán)境。在規(guī)劃城市用地的使用上應將市域內生態(tài)敏感區(qū)和重點保護的土地作為禁建區(qū)，同時應注重保護濕地、植被等自然因子。

該研究選擇多個影響因子并利用GIS技術分析了武漢市土地資源的城鎮(zhèn)發(fā)展適宜性，研究結果具有現(xiàn)實意義。但由于影響城鎮(zhèn)發(fā)展適宜性評價的因子較多而且地區(qū)間差異較大，因此在確定評價因子和權重等方面可能還不夠全面和準確。在以后的研究中應進一步完善，以得到更精確的結果。

［1］李亞奇，席廣亮，秦天.基于GIS技術的城鎮(zhèn)建設用地適宜性評價研究——以嘉興市鳳橋鎮(zhèn)為例［J］.海南科學，2010(28):105 -107.

［2］林恒章.遙感與祖國改革開放同行——中國科學院遙感應用研究所成立與發(fā)展回顧［M］.北京:中國科學技術出版社，2004:231.

［3］武漢市資料.百度百科.［EB/OL］.(2015)(2015-02-25)http://baike.baidu.com/link?url=8x1SlwB4yjM7sgjc - UCP3U0fraqZo Wvem3zLrZb9b6Jy314jARw3NMLYiUn1iwKiB7XZ8YfLKaGDwkdubsVDM Pre79UuCB6YoAH-6yl-xRF44gNgqLlDlxUExLQ7yiHok_B-vwl41lRmw MW9lUtCvjJiMtEHqTJsC59a1ZKf-uclCo6rvE3do_GVc_NlqpA_.

［4］張浩，趙智杰.基于GIS的城市建設用地生態(tài)適宜性評價研究——綜合生態(tài)足跡分析與生態(tài)系統(tǒng)服務［J］.北京大學學報:自然科學版，2011，47(3):531 -538.

［5］AZADNIA A H，SAMAN M Z M，WONG K Y.Sustainable supplier selection and order lot-sizing:an integrated multi-objective decision-making process［J］.International Journal of Production Research，2015，53(2):383-408.

［6］鄭宇，胡業(yè)翠，劉彥隨，等.山東省土地適宜性空間分析及其優(yōu)化配置研究［J］.農業(yè)工程學報，2005(2):60-65.

［7］郭金玉，張忠彬，孫慶云.層次分析法的研究與應用［J］.中國安全科學學報，2008(5):148-153.

［8］徐永輝，陳文瑞.土地適宜性評價關鍵問題探討［J］.莆田高等?？茖W校學報，2000(3):23 -25.

［9］YEH D，CHENG C.Recommendation system for popular tourist attractions in Taiwan using Delphi panel and repertory grid techniques［J］.Tourism Management，2015，46:164 -176.

［10］曹雋雋，周勇，葉青清，等.基于模糊灰色物元與克里格插值的縣級耕地質量分等更新方法研究［J］.經濟地理，2012(11):131-137.

［11］鄭宇，胡業(yè)翠，劉彥隨，等.山東省土地適宜性空間分析及其優(yōu)化配置研究［J］.農業(yè)工程學報，2005，21(2):60 -65.

［12］陳然.基于GIS的農村土地生態(tài)適宜性評價及應用研究［D］.南京:南京農業(yè)大學，2011.

［13］王翠萍，劉寶軍，孫景梅，等.基于GIS的臨澤縣土地生態(tài)適宜性評價［J］.林業(yè)資源管理，2011(3):78-82.

［14］陳燕飛，杜鵬飛，鄭筱津.基于GIS的南寧市建設用地生態(tài)適宜性評價［J］.清華大學學報:自然科學版，2006，46(6):801 -804.

［15］肖捷穎.基于GIS和Landsat TM/ETM+遙感資料分析的城市土地利用/覆蓋變化研究——以石家莊為例［D］.石家莊:河北師范大學，2002.