京津冀都市圈城鎮(zhèn)擴展時空過程及其未來情景預(yù)測

徐新良 通拉嘎 鄭凱迪 張學(xué)霞

(1.中國科學(xué)院地理科學(xué)與資源研究所資源與環(huán)境信息系統(tǒng)國家重點實驗室，北京100101;2.北京林業(yè)大學(xué)水土保持與荒漠化防治教育部重點試驗室，北京100083)

近年來，京津冀都市圈城市化發(fā)展迅速［1－2］，城市建成區(qū)范圍飛速擴張，雖然快速城市化帶來了一定的經(jīng)濟效益，但伴隨其中的是大量優(yōu)質(zhì)耕地資源被侵占，人地矛盾日益尖銳［3－4］，同時帶來了環(huán)境污染、交通擁擠、居住環(huán)境惡化等一系列問題，從而影響到區(qū)域和國家的可持續(xù)發(fā)展［5－8］。因此分析都市圈城市空間擴展時空特征，模擬、預(yù)測其動態(tài)發(fā)展演變，深入揭示其內(nèi)在機制和成因，從而有效地避免所帶來的矛盾和問題已經(jīng)成為當(dāng)前城市化研究最需要解決的問題之一。城市擴展模擬及多情景預(yù)測，無疑是解決上述問題的最佳途徑之一。近年來，國內(nèi)外學(xué)者通過對城市擴展的模擬研究，發(fā)展了大量的模型，如元胞自動機 CA 模型［9－10］、邏輯回歸模型［11］、CLUE － S 模型［12］、LTM 模 型［13］、SLEUTH 模 型［14－18］等。 目 前，SLEUTH模型已在全球范圍得到了廣泛應(yīng)用，如Reza Rafiee運用SLEUTH模型模擬預(yù)測了伊朗馬什哈德市不同情景下的城市擴展過程［19］，李明杰等運用SLEUTH模型開展了廣州市海珠區(qū)高密度城區(qū)擴展的模擬預(yù)測［20］，劉勇在SLEUTH模型支持下開展了杭州市城市擴展的模擬與評估［14］，吳曉青等利用SLEUTH模型開展了沈陽市城市擴展與土地利用變化的多情景模擬預(yù)測［15］，朱飛鴿等運用SLEUTH模型分析了不同情景下的上海市城市增長方式［16］。

本文以京津冀都市圈為研究區(qū)，利用遙感和GIS技術(shù)手段，分析了1990年以來京津冀都市圈城鎮(zhèn)擴展的時空演變過程，并在此基礎(chǔ)上，運用SLEUTH模型，模擬了不同發(fā)展場景下京津冀都市圈城鎮(zhèn)擴展的未來情景，旨在預(yù)測分析城市擴展的時空演變特征，以及對耕地資源的影響，以期為京津冀都市圈未來城市發(fā)展空間布局規(guī)劃、土地資源有效管理以及城市可持續(xù)發(fā)展戰(zhàn)略和決策的制定，提供科學(xué)依據(jù)和決策參考。

1 研究資料和方法

1.1 研究區(qū)概況

京津冀都市圈包括北京、天津兩個直轄市以及河北省的石家莊、唐山、保定、廊坊、秦皇島、滄州、張家口和承德8個地級市［21］，總面積18.26萬平方公里，占全國國土面積的1.9%。2003年京津冀都市圈總?cè)丝?310萬，占全國總?cè)丝诘?.65%;GRP(地區(qū)生產(chǎn)總值)為11808.19億元，占全國總GRP的9.73%;人均GRP為16153元，是全國平均值的1.97倍［4］。京津冀都市圈不僅是我國最重要的政治、文化中心，也是我國北方最大、發(fā)育最好、現(xiàn)代化程度最高的城鎮(zhèn)和工業(yè)密集區(qū)以及經(jīng)濟核心區(qū)，同時還是我國參與國際經(jīng)濟交流與合作的重要樞紐與門戶［22］。京津冀都市圈的穩(wěn)定與進步對全國的發(fā)展十分重要。因此，研究京津冀都市圈城鎮(zhèn)擴展的時空過程以及未來發(fā)展趨勢具有重要的意義。

1.2 數(shù)據(jù)來源及處理

1.2.1 土地利用數(shù)據(jù)

土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心土地利用現(xiàn)狀數(shù)據(jù)庫。該數(shù)據(jù)庫是在國家科技支撐計劃、中國科學(xué)院知識創(chuàng)新工程重要方向項目等多項重大科技項目的支持下經(jīng)過多年的積累而建立的覆蓋全國陸地區(qū)域的多時相1:10萬比例尺土地利用現(xiàn)狀數(shù)據(jù)集［7－8，23－26］。該數(shù)據(jù)集是以 Landsat TM/ETM 遙感影像為主要數(shù)據(jù)源，通過人工目視解譯生成［27］。土地利用類型包括耕地、林地、草地、水域、居民地和未利用土地6個一級類型以及25個二級類型。本研究利用GIS技術(shù)，從1990、1995、2000、2005和 2008 年五期土地利用數(shù)據(jù)中分別提取居民地中的城鎮(zhèn)用地作為京津冀都市圈城鎮(zhèn)擴展時空過程研究的數(shù)據(jù)基礎(chǔ)。

1.2.2 SLEUTH 模型輸入數(shù)據(jù)

本研究運用SLEUTH模型來開展京津冀都市圈城鎮(zhèn)擴展未來情景模擬和預(yù)測。SLEUTH模型的輸入數(shù)據(jù)包括土地利用、坡度、dem、交通道路以及水域、濕地、自然保護區(qū)等。交通道路數(shù)據(jù)是在1990年Landsat TM遙感影像的基礎(chǔ)上通過遙感解譯生成的交通道路圖，然后在此基礎(chǔ)上，通過1990年和2000年、2000年和2005年遙感影像的對比分析，通過人工目視解譯生成2000年和2005年的交通道路圖。水域、濕地數(shù)據(jù)、自然保護區(qū)是作為SLEUTH模型輸入數(shù)據(jù)的排除層，用來設(shè)置不同場景城市增長的限制因素，用來減緩或改變城市化速率。該層對模型的運行非常重要，在很大程度上控制著模型的運行。本研究利用GIS軟件ArcGIS和遙感圖像處理軟件Envi，通過遙感解譯、分類和格式轉(zhuǎn)化等功能，建立模型運行所要求的gif格式的京津冀都市圈城市擴展時空數(shù)據(jù)集。

1.3 研究方法

1.3.1 城鎮(zhèn)擴展過程分析

基于以往土地利用變化研究相關(guān)成果，在GIS技術(shù)支持下，本文運用單一土地利用類型動態(tài)度、土地利用類型轉(zhuǎn)移矩陣等模型以及空間統(tǒng)計分析方法，開展1990—2008年京津冀都市圈城鎮(zhèn)擴展時空過程分析。

單一土地利用類型動態(tài)度，即某一土地利用類型面積的年變化速率，本文運用城鎮(zhèn)擴展速率、耕地減少速率，來反映研究區(qū)一定時間范圍內(nèi)城鎮(zhèn)面積和耕地面積的變化情況，其表達式為［28－29］:

R為某一土地利用類型動態(tài)度，Ua、Ub分別為研究期初及研究期末該土地利用類型的面積，T為研究時段長，當(dāng)T的時段設(shè)定為年時，R的值就是該研究區(qū)某一土地利用類型面積的年變化速率。

1.3.2 城鎮(zhèn)擴展未來情景預(yù)測

本研究運用SLEUTH模型來模擬和預(yù)測京津冀都市圈城鎮(zhèn)擴展未來情景。在模型模擬過程中，首先初始化繁衍系數(shù)(Breed_coefficient)、擴散系數(shù)(Spread_coefficient)、坡度阻抗系數(shù)(Slope_coefficient)、傳播系數(shù)(Dispersion_coefficient)和道路引力系數(shù)(Road_coefficient)的值，然后通過粗校準(zhǔn)、精校準(zhǔn)、終校準(zhǔn)，逐步細化各參數(shù)范圍，每個校準(zhǔn)階段包含多次蒙特卡羅迭代，并用Lee－Sallee指數(shù)來評價模擬結(jié)果與實際情況的空間匹配程度，從而確定下一步的參數(shù)范圍，最后得出一組最佳參數(shù)，用于模型模擬。本研究將京津冀都市圈2005和2008年兩期土地利用類型空間分布、城鎮(zhèn)用地空間分布，以及道路、坡度等數(shù)據(jù)輸入SLEUTH模型，采用激活LCDM模塊，對城市擴展進行模擬，經(jīng)過模型校準(zhǔn)，最后得到的模擬參數(shù)分別為擴散系數(shù)2、繁殖系數(shù)41、擴展系數(shù)3、坡度系數(shù)15、道路引力系數(shù)30。通過Kappa系數(shù)檢驗，SLEUTH模型模擬的2008年京津冀都市圈土地利用類型分布與遙感監(jiān)測結(jié)果的Kappa系數(shù)為0.742。可見，SLEUTH模型能較好的模擬預(yù)測京津冀都市圈城鎮(zhèn)擴展的未來情景。

為了對京津冀都市圈城鎮(zhèn)擴展未來情景進行預(yù)測，本研究在SLEUTH模型中，通過設(shè)置排除層，制定了3種城鎮(zhèn)擴展未來情景:1)現(xiàn)有發(fā)展趨勢情景(MSQ)，此情景只有水體(河流、湖泊、水庫、坑塘)被設(shè)置為排除層(100%的概率不被城市化)，其他因素對城市開發(fā)不構(gòu)成限制。2)生態(tài)環(huán)境保護管理情景(ES)，此情景考慮生態(tài)環(huán)境保護，將城鎮(zhèn)周邊的水體、濕地、自然保護區(qū)設(shè)置為100%不被城市擴展侵占，將基本農(nóng)田以及連續(xù)成片分布的林地設(shè)置為80%不被城市擴展侵占。3.緊湊城市發(fā)展情景(CC)，此情景設(shè)置距城區(qū)1－2km之間的緩沖區(qū)內(nèi)有50%的概率不被開發(fā)，2km以上的緩沖區(qū)有70%的概率不被開發(fā)。本研究基于上述SLEUTH模型校準(zhǔn)獲取的模型參數(shù)，分別開展了3種情景下京津冀都市圈未來2020年、2030年城鎮(zhèn)擴展的模擬預(yù)測。

2 結(jié)果與分析

2.1 1990～2008年京津冀都市圈城鎮(zhèn)擴展時空特征

1990－2008 年期間，京津冀都市圈城鎮(zhèn)擴展明顯(表1)，18年間，城鎮(zhèn)用地面積由1990年的1.17萬km2增長到2008 年的1.57 萬 km2，凈增加 4022.37km2，城鎮(zhèn)擴展速率平均為每年1.91%。從城鎮(zhèn)擴展的時間過程看，1990年以來，京津冀都市圈城鎮(zhèn)擴展速率呈現(xiàn)逐步減緩的趨勢，其中1990－2000年，城鎮(zhèn)擴展速率為2.07%，2000－2005年城鎮(zhèn)擴展速率下降為1.74%，而2005－2008年城鎮(zhèn)擴展速率僅為0.79%。從北京、天津和河北城鎮(zhèn)擴展的地區(qū)差異看(見表1)，北京城鎮(zhèn)擴展速率最高，從1990年到2008年城鎮(zhèn)面積擴展了1328.68km2，18年間，城鎮(zhèn)年平均擴展速率為5.05%，其次為天津，城鎮(zhèn)面積擴展504.62km2，年平均擴展速率為1.6%，而河北8地市城鎮(zhèn)擴展速率最低，城鎮(zhèn)面積擴展了2189.07km2，年平均擴展速率僅為1.43%。在時間過程上，北京、天津和河北三地區(qū)城鎮(zhèn)擴展速率均呈減緩趨勢，其中北京、河北均是在1990－2000年城鎮(zhèn)擴展速率最高，分別為5.32%和1.85%，而天津在2000－2005年城鎮(zhèn)擴展速率最高，為4.28%。

表1 1990－2008年京津冀都市圈城鎮(zhèn)擴展面積和速率統(tǒng)計表

1990－2008 年京津冀都市圈城鎮(zhèn)擴展對耕地資源的占用特別明顯(見表2)，18年間，耕地面積減少了3507.62 km2，耕地年減少速率為0.23%。從時間過程看，1990－2000年城鎮(zhèn)擴展占用耕地資源尤其明顯，10年間耕地減少2258.60km2，年減少速率達到了為0.26%。北京、天津和河北三地區(qū)城鎮(zhèn)擴展對耕地資源的影響有所差別，其中北京城鎮(zhèn)擴展占用耕地資源最為突出，18年間耕地面積下降1419.47km2，年下降速率為1.34%，而天津城鎮(zhèn)擴展對耕地的影響較小，1990－2008年耕地面積減少了459.13km2，年減少速率為0.35%。在時間過程上，北京、河北均是在1990－2000年耕地減少速率最高，分別為1.6%和0.15%，而天津在2000－2005年耕地減少速率最高，為0.57%。

表2 1990－2008年京津冀都市圈耕地減少面積和速率統(tǒng)計表

從城鎮(zhèn)擴展的空間格局看，京津冀都市圈依托北部燕山山脈和西部太行山脈的生態(tài)屏障作用，城鎮(zhèn)擴展主要集中在東南部的平原城市群，城鎮(zhèn)擴展呈現(xiàn)顯著的“雙核、一帶和兩軸、多中心”的空間格局特征。北京和天津是京津冀地區(qū)發(fā)育和發(fā)展的重要動力源，城鎮(zhèn)擴展最為突出。近20年來，北京城市空間以舊城為核心，呈單中心圈層狀不斷向外擴大，四環(huán)、五環(huán)、六環(huán)的修建更是起到了推波助瀾的作用，從而使得整個市區(qū)呈現(xiàn)出一種“攤大餅”的外延發(fā)展的模式。天津城市擴展也是以向周邊團塊蔓延為主，其中以西北—東南走向沿北運河—海河軸向擴展一直是城市擴展的主軸，近年來濱海軸擴展速度也在加快。受北京和天津兩核“圈層輻射力”和20世紀(jì)90年代以來的“外向型城市化力”的共同影響，形成了城鎮(zhèn)擴展的“一帶”，即“北京—廊坊—天津—塘沽走廊帶”。城鎮(zhèn)擴展的“兩軸”則包括傳統(tǒng)的山前城鎮(zhèn)密集軸(即北京—保定—石家莊—邢臺—邯鄲等太行山脈山前地區(qū)和北京—薊縣—遵化—秦皇島的燕山山脈山前城鎮(zhèn)密集軸)、沿海新興城鎮(zhèn)密集軸(以秦皇島—京唐港—曹妃甸港—天津港—黃驊港為主要節(jié)點的城鎮(zhèn)密集軸)。在兩核、一帶和兩軸的共同作用了，京津冀都市圈呈現(xiàn)出以北京、天津、唐山、保定、石家莊等多中心城市擴展為主體，沿京石、京秦和京津等交通線軸向延伸和扇形擴張的總體發(fā)展趨勢(圖1)。

圖1 1990－2008年京津冀都市圈城鎮(zhèn)擴展空間分布圖

2.2 京津冀都市圈城鎮(zhèn)擴展未來情景

本研究利用SLEUTH模型對三種情景(現(xiàn)有發(fā)展趨勢情景(MSQ)、生態(tài)環(huán)境保護管理情景(ES)、緊湊城市發(fā)展情景(CC))下的京津冀都市圈2008－2030年城鎮(zhèn)擴展進行了模擬預(yù)測。從預(yù)測結(jié)果看(見表3)，三種情景下，2008－2020年京津冀都市圈城鎮(zhèn)擴展速率較歷史時期均有一定減緩，其中在現(xiàn)有發(fā)展趨勢情景(MSQ)下，京津冀都市圈城鎮(zhèn)面積將進一步擴大，2020年城鎮(zhèn)面積將達到17766km2，未來12年城鎮(zhèn)面積將擴大2035km2，年擴展速率為1.08%，而在生態(tài)環(huán)境保護管理情景(ES)和緊湊城市發(fā)展情景(CC)下，城鎮(zhèn)擴展速率將進一步減緩，分別下降為0.09%和0.15%，這主要是由于在生態(tài)環(huán)境保護管理情景(ES)下，城鎮(zhèn)擴展過程中特別注重保護水體、濕地、自然保護區(qū)，以及城市周邊的基本農(nóng)田和連續(xù)成片分布的林地等，因此城鎮(zhèn)擴展被抑制，城鎮(zhèn)擴展速率明顯下降;而在緊湊城市發(fā)展情景(CC)下，城鎮(zhèn)發(fā)展過程中集中于城市內(nèi)部填充和調(diào)整，因此城鎮(zhèn)擴展速率也明顯放緩。對比不同的未來城鎮(zhèn)擴展情景，生態(tài)環(huán)境保護管理情景(ES)、緊湊城市發(fā)展情景(CC)將極大的抑制城鎮(zhèn)邊緣擴張，因此，在未來京津冀都市圈城鎮(zhèn)擴展過程中，通過有效的生態(tài)環(huán)境保護管理和制訂相應(yīng)的城市發(fā)展策略，能有效的控制都市圈城鎮(zhèn)擴展速率，避免城鎮(zhèn)邊緣盲目擴張所帶來的一系列資源、環(huán)境問題和矛盾。

從北京、天津、河北三地市2008－2020年城鎮(zhèn)擴展的發(fā)展趨勢看(見表3)，在現(xiàn)有發(fā)展趨勢情景(MSQ)下河北的城鎮(zhèn)擴展速率最高，年擴展速率為1.1%，而天津的城鎮(zhèn)擴展速率相對較低，年擴展速率為0.99%，在生態(tài)環(huán)境保護管理情景(ES)下，北京、河北的城鎮(zhèn)擴展速率均為0.1%。而天津的城鎮(zhèn)擴展速率相對較低，為0.06%，而在緊湊城市發(fā)展情景(CC)下，天津的城鎮(zhèn)擴展速率最高，為0.17%，北京的城鎮(zhèn)擴展速率最最低，年擴展速率為0.13%。

從城鎮(zhèn)擴展對耕地資源的影響看(見表4)，2008－2020年京津冀都市圈城鎮(zhèn)擴展仍然將侵占一定的耕地資源，其中在現(xiàn)有發(fā)展趨勢情景(MSQ)下，耕地面積將減少1792km2，年減少速率為0.18%，而在生態(tài)環(huán)境保護管理情景(ES)和緊湊城市發(fā)展情景(CC)下耕地減少的態(tài)勢有所緩解，耕地面積將分別減少65 km2和270km2，年減少速率分別下降為0.01%和0.03%。從北京、天津、河北未來城鎮(zhèn)擴展對耕地資源的影響看，北京耕地面積減少的速率最高，在MSQ、ES和CC三種情景下，耕地年減少速率分別為0.54%、0.03% 和0.08%，而河北城鎮(zhèn)擴展對耕地的影響最弱，在三種情景下，耕地年減少速率僅為0.15%、0.004% 和 0.02%。

表3 2008－2020年京津冀都市圈城鎮(zhèn)擴展面積和速率統(tǒng)計表

表4 2008－2020年京津冀都市圈耕地減少面積和速率統(tǒng)計表

2020年以后京津冀都市圈城鎮(zhèn)擴展的速率將進一步減緩(表5)，其中在現(xiàn)有發(fā)展趨勢情景(MSQ)下2030年京津冀都市圈城鎮(zhèn)面積將達到18201km2，10年間城鎮(zhèn)面積將擴大435km2，年擴展速率為0.24%，而在生態(tài)環(huán)境保護管理情景(ES)和緊湊城市發(fā)展情景(CC)下，城鎮(zhèn)擴展速率將進一步減緩，分別下降為0.08%和0.1%。從北京、天津、河北三地市2020－2030年城鎮(zhèn)擴展的發(fā)展趨勢看(表5)，在現(xiàn)有發(fā)展趨勢情景(MSQ)下北京、河北的城鎮(zhèn)擴展速率均為0.25%，而天津的城鎮(zhèn)擴展速率相對較低，為0.23%，在生態(tài)環(huán)境保護管理情景(ES)下，河北的城鎮(zhèn)擴展速率相對較高，為0.09%。而天津的城鎮(zhèn)擴展速率相對較低，為0.05%，而在緊湊城市發(fā)展情景(CC)下，天津的城鎮(zhèn)擴展速率最高，為0.13%，河北的城鎮(zhèn)擴展速率最最低，為0.09%。2020年以后京津冀都市圈城鎮(zhèn)擴展對耕地資源的影響已非常微弱(表6)，其中在現(xiàn)有發(fā)展趨勢情景(MSQ)下，京津冀都市圈耕地面積將減少379km2，耕地年減少速率為0.05%，而在生態(tài)環(huán)境保護管理情景(ES)和緊湊城市發(fā)展情景(CC)年減少速率僅為0.01% 和0.02%。從北京、天津、河北三地市耕地減少速率差異看，北京耕地減少的速率仍然最高，在MSQ、ES和CC三種情景下，耕地年減少速率為分別為0.16%、0.03%和 0.06%。注:MSQ現(xiàn)有發(fā)展趨勢情景、ES生態(tài)環(huán)境保護管理情景、CC緊湊城市發(fā)展情景

表5 2020－2030年京津冀都市圈城鎮(zhèn)擴展面積和速率統(tǒng)計表

表6 2020－2030年京津冀都市圈耕地減少面積和速率統(tǒng)計表

3 結(jié)論與討論

從1990－2008年城鎮(zhèn)擴展時空過程的分析看，歷史時期京津冀都市圈城鎮(zhèn)擴展時空過程具有以下顯著特征:

(1)京津冀都市圈城鎮(zhèn)擴展明顯，1990－2008年城鎮(zhèn)面積凈增加4022.37km2，年平均擴展速率為1.91%。在時間過程上，城鎮(zhèn)擴展呈現(xiàn)逐步減緩的趨勢，1990－2000年，年擴展速率為2.07%，2000－2005年下降為1.74%，而2005－2008年僅為0.79%。

(2)在城鎮(zhèn)擴展的空間格局上，京津冀都市圈城鎮(zhèn)擴展呈現(xiàn)顯著的“雙核、一帶和兩軸、多中心”的空間格局特征。雙核為北京和天津，一帶為“北京—廊坊—天津—塘沽走廊帶”，兩軸為傳統(tǒng)的山前城鎮(zhèn)密集軸和沿海新興城鎮(zhèn)密集軸，多中心則包括北京、天津、唐山、保定、石家莊等多個核心城市。

(3)京津冀都市圈城鎮(zhèn)擴展對耕地資源的占用特別明顯，1990－2008年的 18年間，耕地面積凈減少了3507.62km2，年減少速率為0.23%。北京城鎮(zhèn)擴展占用耕地資源最為突出，18年間耕地面積凈減少了1419.47km2，年減少速率為1.34%。

未來京津冀都市圈城鎮(zhèn)擴展將呈現(xiàn)以下特征:

(1)2008－2020年京津冀都市圈城鎮(zhèn)擴展速率較歷史時期有一定減緩，在MSQ、ES和CC三種情景下，城鎮(zhèn)年擴展速率將分別下降為1.08%、0.09%和0.15%。2020年以后京津冀都市圈城鎮(zhèn)擴展的速率將進一步減緩，在MSQ、ES和 CC三種情景下城鎮(zhèn)年擴展速率為0.24%、0.08%和 0.1%。

(2)2008－2020年北京、天津、河北三地市城鎮(zhèn)擴展將呈現(xiàn)差異發(fā)展的態(tài)勢，在MSQ情景下河北的城鎮(zhèn)擴展速率最高，年擴展速率為1.1%;在ES情景下，北京、河北的城鎮(zhèn)擴展速率均較高，為0.1%;在CC情景下，天津的城鎮(zhèn)擴展速率較高，為0.17%。2020以后北京、天津、河北三地市城鎮(zhèn)擴展的速率也都將進一步減緩。

(3)京津冀都市圈未來城鎮(zhèn)擴展仍然將侵占一定的耕地資源，2008－2020年，在MSQ、ES和CC三種情景下，耕地年減少速率將分別下降為0.18%、0.01%和0.03%。2020年以后京津冀都市圈城鎮(zhèn)擴展對耕地資源的影響將非常微弱。

(編輯:徐天祥)

［1］余靜文，趙大利.城市群落的崛起、經(jīng)濟績效與區(qū)域收入差距—基于京津冀、長三角和珠三角城市圈的分析［J］.中南財經(jīng)政法大學(xué)學(xué)報，2010，(4):15－20.

［2］陳潔，陸鋒.京津冀都市圈城市區(qū)位與交通可達性評價［J］.地理與地理信息科學(xué)，2008，24(2):53－56.

［3］紀(jì)純，張士功.近年來中國耕地數(shù)量變化及其途徑和原因分析［J］.現(xiàn)代化農(nóng)業(yè)，2005，(10):24 － 26.

［4］王磊，劉逢媛，李雙成，等.耕地非農(nóng)化格局的演變及其影響因子分析—以京津冀都市圈為例［J］.中國土地科學(xué)，2008，22(1):32－38.

［5］Norman L M，Guertin D P，F(xiàn)eller M .A Coupled Model Approach to Reduce Nonpoint－Source Pollution Resulting from Predicted Urban Growth:A Case Study in the Watershed［J］.Urban Geography，2008，29(5):496－516.

［6］Dietzel C，Clarke K.The effect of disaggregating land use categories in cellular automata during model calibration and forecasting［J］.Computers，Environment and Urban Systems，2006，30(1):78 －101.

［7］劉紀(jì)遠，張增祥，徐新良，等.21世紀(jì)初中國土地利用變化的空間格局與驅(qū)動力分析［J］.地理學(xué)報，2009，64(12):1411－1420.

［8］劉紀(jì)遠，張增祥，莊大方，等.20世紀(jì)90年代中國土地利用變化時空特征及其成因分析［J］.地理研究，2003，22:1－12.

［9］Long Y，Mao Q Z，Dang A R.Beijing Urban Development Model:Urban Growth Analysis and Simulation［J］.Tsinghua Science and Technology，2009，14(6):782 －794.

［10］Silva E A，Clarke K C.Calibration of the SLEUTH urban growth model for Lisbon and Porto，Portugal［J］.Computers，Environment and Urban Systems，2002，26:525 －552.

［11］楊勇，任志遠，李開宇.基于GIS的西安市城市擴展與模擬研究［J］.人文地理，2010，25(2):95－98.

［12］ Batisani Nnyaladzi，Yarnal Brent.Urban expansion in Centre County，Pennsylvania:Spatial dynamics and landscape transformations［J］.Applied Geography，2009，29:235 －249.

［13］徐新良，閔稀碧，田展.基于GIS與LTM模型的上海市城市空間擴展模擬與預(yù)測［J］.中國人口資源與環(huán)境，2010，20(5):136 －139.

［14］劉勇.SLEUTH模型支持下的城市擴展情景模擬與評估［J］.經(jīng)濟地理，2009，29(10):1661－1664.

［15］吳曉青，胡遠滿，賀紅士，等.沈陽市城市擴展與土地利用變化多情景模擬［J］.地理研究，2009，28(5):1264－1275.

［16］朱飛鴿，胡瀚文，沈興華，等.基于SLEUTH模型的上海城市增長預(yù)測［J］.生態(tài)學(xué)雜志，2011，30(9):2107－2114.

［17］劉小平，黎搖夏，彭曉娟.“生態(tài)位”元胞自動機在土地可持續(xù)規(guī)劃模型中的應(yīng)用［J］.生態(tài)學(xué)報，2007，27(6):2391 －2402.

［18］Clarke K C，Gaydos L J.Loose-coupling a cellular automaton model and GIS:Long-term urban growth prediction for San Francisco and Washington/Baltimore［J］.International Journal of Geographical Information Science，1998，12:699－714.

［19］ Rafiee R，Mahiny A S，Khorasani X，et al.Simulating urban growth in Mashad City，Iran through the SLEUTH model(UGM)［J］.Cities，2009，26(1):19 －26.

［20］李明杰，錢樂祥，吳志峰，等.廣州市海珠區(qū)高密度城區(qū)擴展SLEUTH 模型模擬［J］.地理學(xué)報，2010，65(10):1163 －1172.

［21］樊杰.京津冀都市圈區(qū)域綜合規(guī)劃研究［M］.北京:科學(xué)出版社，2008.

［22］孫文博，苗澤華，董莉.基于DEA方法的京津冀都市圈生態(tài)效率評價［J］.商業(yè)時代，2011，(4):137－138.

［23］劉紀(jì)遠.中國資源環(huán)境遙感宏觀調(diào)查與動態(tài)研究［M］.北京:中國科學(xué)技術(shù)出版社，1996.171－188.

［24］劉紀(jì)遠，布和敖斯?fàn)?中國土地利用變化現(xiàn)代時空特征的研究:基于衛(wèi)星遙感數(shù)據(jù)［J］.第四紀(jì)研究，2000，20(3):229 －239.

［25］劉紀(jì)遠，劉明亮，莊大方，等.中國近期土地利用變化的空間格局分析［J］.中國科學(xué)(D輯)，2002，32(12):1031－1040.

［26］劉紀(jì)遠，張增祥，莊大方，等.20世紀(jì)90年代中國土地利用變化的遙感時空信息研究［M］.北京:科學(xué)出版社，2005.241－243.

［27］徐新良，劉紀(jì)遠，莊大方.國家尺度土地利用/覆被變化遙感監(jiān)測方法［J］.安徽農(nóng)業(yè)科學(xué)，2012，40(4):2365 －2369

［28］Ma，Y，Xu，R.Remote sensing monitoring and driving force analysis of urban expansion in Guangzhou City，China［J］.Habitat International，2010，34(2)，228 –235.

［29］Xiao J Y，Shen Y J，Ge J F，et al.Evaluating urban expansion and Land use change in Shijiazhuang，China，by using GIS and remote sensing［J］.Landscape and Urban Planning，2006，75(1 － 2)，69–80.

［30］鄔建國.景觀生態(tài)學(xué)——格局、過程、尺度與等級［M］.北京:高等教育出版社，2000.