張玲玲，頡耀文，2*，穆亞超(.蘭州大學資源環(huán)境學院，甘肅蘭州730000;2.蘭州大學西部環(huán)境教育部重點實驗室，甘肅蘭州730000)

土地利用變化是一個復雜的過程，是自然因素和人為因素共同作用的結果，如何準確和合理地模擬土地利用的變化過程是許多人研究的重點，也是一個難點［1］。與GIS結合的元胞自動機(CA)模型被越來越多地運用于對城市空間擴展的模擬研究中［2-3］。CA模型通過轉換規(guī)則的適當定義，不僅可以很好地模擬土地利用的時空復雜性，也可模擬城市發(fā)展的某些規(guī)劃場景。該文以重慶市沙坪壩區(qū)為例，基于GeoSOS平臺，運用Logistic—CA模型模擬重慶市沙坪壩區(qū)土地利用擴展的過程及預測重慶市沙坪壩地區(qū)未來不同情景下的土地利用發(fā)展格局，從而為區(qū)域可持續(xù)發(fā)展提供決策支持。

1 資料與方法

1.1 研究區(qū)域概況 沙坪壩區(qū)是重慶市主城區(qū)和都市圈的重要組成部分，位于重慶的中心［4］，地處 106°14'36″～106°31'35″E，29°27'13″～ 29°46'36″N，面積 396.20 km2，轄小龍坎、沙坪壩、磁器口等11個街道，覃家崗、井口、歌樂山等12個鎮(zhèn)和聯(lián)芳村經(jīng)濟園區(qū)(圖1)。東與江北區(qū)、渝北區(qū)隔江相望，東南與渝中區(qū)接壤，東北與北碚區(qū)相連，南與九龍坡區(qū)相靠，西與壁山縣相依。

1.2 數(shù)據(jù)來源及處理 該數(shù)據(jù)源包括2004、2006年的土地利用數(shù)據(jù)(圖2、3)、2008年道路現(xiàn)狀、鎮(zhèn)政府點分布以及高程。借助GIS平臺對數(shù)據(jù)進行處理及格式轉換，轉換為二進制文件(.txt)，使之能在GeoSOS中進行相應的模擬操作。首先，將2004、2006年的土地利用數(shù)據(jù)進行格式轉換，土地利用類型有農(nóng)田、森林、城市用地、果園和水域五大類型。其次，從2008年道路現(xiàn)狀圖中提取普通公路(Road)、快速道路(fastRoad)2個二進制文件。坡度影響因子從高程圖中提取并轉換，鎮(zhèn)政府點分布圖直接進行轉換。最終得到離鎮(zhèn)政府點距離、快速道路距離、普通公路距離3個距離變量、坡度變量。注意在進行處理及格式轉換過程中，要進行標準化處理、投影轉換，保持投影及研究區(qū)范圍的一致性。

圖1 重慶市沙坪壩區(qū)鄉(xiāng)鎮(zhèn)界圖

1.3 CA模型 元胞自動機(Cellular Automata，CA)具有強大的空間運算能力，常用于自組織系統(tǒng)演變過程的研究。它是一種時間、空間、狀態(tài)都離散，空間相互作用和時間因果關系都為局部的網(wǎng)格動力學模型，具有模擬復雜系統(tǒng)時空演化過程的能力［5］。

CA的一個主要特征是CA與GIS的耦合。CA和GIS的耦合能使CA模擬出與實際情況更為接近的模擬結果。目前也有多個系列的CA模擬軟件，在模擬城市格局變化時常把GIS作為集成平臺，用GIS軟件提供的二次開發(fā)宏語言來開發(fā)城市CA，如Arc/Info GRID模塊開發(fā)城市CA［6］。

圖2 重慶市沙坪壩區(qū)2004年土地利用類型圖

圖3 重慶市沙坪壩區(qū)2006年土地利用類型圖

2 結果與分析

2.1 基于GeoSOS平臺的CA模型土地利用變化模擬 GeoSOS平臺提供的轉換規(guī)則有多準則判斷、邏輯回歸、主成分分析、神經(jīng)網(wǎng)絡、決策樹5種方法。管超等［7］利用該系統(tǒng)進行了基于耦合地理模擬優(yōu)化系統(tǒng)GeoSOS的旅順口區(qū)城市空間擴展等應用并得到了很好的效果，Logistic回歸方法將土地利用類型中的某一小塊作為一個元胞，只需判斷某個元胞的狀態(tài)在［t，t+1］時刻間是否從0變化到1，就可判斷該元胞所代表的地塊是否城市化，這就是邏輯回歸CA模型的最基本思想，用其在元胞空間的狀態(tài)來確定［8］。因此該文將采用GeoSOS來對重慶市沙坪壩區(qū)土地利用變化進行模擬，并以Logistic—CA作為理論模型，對重慶市沙坪壩區(qū)2006年土地利用情況模擬，定義1代表的是耕地，2代表的是林地，3代表的是城市用地，4代表的是園地，5代表的是水系，并設置好顏色。其中農(nóng)田、果園可轉換為城市用地，森林、水域不可轉換為城市用地，自動計算各空間變量的回歸系數(shù)。

設置轉換總量為1 000 000個柵格單元，迭代次數(shù)為2次，模擬結果精度為95.2%，得到的結果利用Arcgis分析，把原始的2006年的土地利用類型圖像(圖3)與基于邏輯回歸的CA模型模擬的2006年結果圖像(圖4)進行疊加分析，發(fā)現(xiàn)模擬結果是極其接近的(表1)。

圖4 重慶市沙坪壩區(qū)2006年土地利用類型模擬

2.2 模型的評價 城市用地的數(shù)量在不斷增多，而在無較強約束和管理條件下，對區(qū)域生態(tài)安全起重要作用的森林和水域也不可避免地受到城市建設用地的侵占。從土地利用類型分布上看，分布基本上接近，尤其在沙坪壩東南部的井口鎮(zhèn)、井口街道、沙坪壩、歌樂山鎮(zhèn)的模擬精度極高。此外，值得注意的是在曾家鎮(zhèn)及虎溪鎮(zhèn)交界地帶及其周邊有一塊集中的城市區(qū)，在模擬圖像上是沒有的。說明模擬的結果仍然與實際的土地利用/覆蓋有一定的誤差。從城市用地的發(fā)展走向來看，除曾家鎮(zhèn)及虎溪鎮(zhèn)的交界地帶，重慶市沙坪壩區(qū)的發(fā)展大體還是集中在該區(qū)域的東南部，即沙坪壩、井口鎮(zhèn)、井口街道、歌樂山鎮(zhèn)，其他地區(qū)主要依托先前的地區(qū)周邊逐步發(fā)展起來。

表1 2004、2006年研究區(qū)土地利用各類型面積比較 km2

從模擬結果可推測，未來沙坪壩區(qū)的土地利用變化仍然會按照2004～2006年的趨勢，依托現(xiàn)有城市逐步實現(xiàn)鄉(xiāng)村、果園等向城市用地的轉換。曾家鎮(zhèn)及虎溪鎮(zhèn)交界地帶也依然會依托自身優(yōu)勢發(fā)展周邊，逐步完成城鄉(xiāng)轉換，統(tǒng)籌發(fā)展。并且按照2004～2006年的發(fā)展趨勢，未來沙坪壩區(qū)建設用地面積將大規(guī)模增加，其他用地類型面積均有不同程度減少，其中耕地流失量最大，林地、果園、農(nóng)田等建設用地也有一定程度的減少，水域面積也有一定縮減。經(jīng)分析發(fā)現(xiàn)，由城市化進程導致的城市建設用地擴張在未來土地利用演變過程中仍占據(jù)主導地位。其中，農(nóng)田(耕地)、果園及其他非建設用地是城市建設用地的主要來源。從城市擴展格局上看，未來城市擴展將集中在現(xiàn)有城鎮(zhèn)邊緣，沿道路發(fā)展的趨勢也較為明顯。此外，虎溪大學城區(qū)域也將得到較高程度的發(fā)展。

2.3 土地利用功能性發(fā)展分析 考慮不同的城市發(fā)展政策，調整模型參數(shù)，實現(xiàn)不同管理情景下土地利用變化預測。除上述目前趨勢發(fā)展情景，還設計了城市規(guī)劃與區(qū)域發(fā)展情景和生態(tài)環(huán)境保護管理情景。

城市規(guī)劃與區(qū)域發(fā)展情景，即考慮重慶市城市整體發(fā)展規(guī)劃、沙坪壩區(qū)發(fā)展規(guī)劃及其他建設區(qū)域開發(fā)政策對未來土地利用的可能影響。生態(tài)環(huán)境保護管理情景，限制城市用地侵占連續(xù)成片的林地、基本農(nóng)田、水體和東部水源涵養(yǎng)地等對區(qū)域生態(tài)安全起重要作用的生態(tài)系統(tǒng)，保護生態(tài)環(huán)境，尋求城市區(qū)域的可持續(xù)發(fā)展。城市規(guī)劃與區(qū)域發(fā)展情景和生態(tài)環(huán)境保護情景下，土地利用類型面積變化趨勢與目前趨勢發(fā)展情景下的發(fā)展變化相似，即城市建設用地面積增加，其他用地類型面積減少。相比而言，在生態(tài)環(huán)境保護情景下，應鼓勵城鎮(zhèn)沿現(xiàn)有建設用地向外邊緣式擴展，城市建設用地面積增加較慢，但是耕地、林地、水域和未利用土地類型損失面積比例較低，說明生態(tài)環(huán)境保護政策對重要的生態(tài)系統(tǒng)和自然資源起到良好的保護作用。城市規(guī)劃與區(qū)域發(fā)展情景下，城市用地擴展受到較強的空間限制，邊緣增長類型受到約束，城市增長面積低于目前趨勢發(fā)展情景，但是耕地的損失比例較高。從城市用地擴展空間格局上來看，城市規(guī)劃與區(qū)域發(fā)展情景和生態(tài)環(huán)境保護情景與目前趨勢發(fā)展情景存在明顯的空間差異。城市規(guī)劃與區(qū)域發(fā)展情景下，城市用地擴展集中政策開發(fā)區(qū)內(nèi)，新城鎮(zhèn)不斷產(chǎn)生，且城鎮(zhèn)沿道路擴展的趨勢較為明顯。相比而言，生態(tài)環(huán)境保護情景下，城市用地格局非常緊湊?？偟膩碚f，2種管理情景下城市建設用地對具有重要生態(tài)功能的林地威脅較小，但對非基本農(nóng)田之外的耕地或農(nóng)田占用較多。

3 結論與討論

(1)研究表明，重慶市沙坪壩區(qū)城市的空間擴張以外延式為主，集中在現(xiàn)有城鎮(zhèn)邊緣，沿道路發(fā)展的趨勢也較為明顯。對重慶市沙坪壩區(qū)2004～2006年的土地利用情況模擬結果表明，運用Logistic—CA模型能夠較好地模擬土地利用變化，總體精度達到95.2%。在城市規(guī)劃與區(qū)域發(fā)展情景和生態(tài)環(huán)境保護管理情景2種管理情景下，未來的重慶市沙坪壩區(qū)城市建設用地將持續(xù)增加，大量的耕地資源被侵占。該文基于Logistic—CA模型對重慶市沙坪壩區(qū)土地利用模擬主要采用歷史數(shù)據(jù)來優(yōu)化參數(shù)，是一種靜態(tài)的轉換規(guī)則，這種方法對城市規(guī)劃調整、重大事件、重大建設項目和行政區(qū)劃調整影響無法明確體現(xiàn)，尤其在考慮大學城區(qū)域的政策規(guī)劃下模擬城市的擴展，這是導致沙坪壩區(qū)西南部曾家鎮(zhèn)及虎溪鎮(zhèn)交界處土地利用模擬誤差的主要原因。

(2)由于存在數(shù)據(jù)源、CA模型不確定性、未考慮到的政策規(guī)劃［9］以及空間尺度［10］利用等不同程度的誤差，今后可以通過優(yōu)化CA模型參數(shù)，定量化社會、經(jīng)濟、政策因素等提高模型精度，也可以發(fā)揮GIS軟件的功能，使之與GeoSOS集成，這樣將會給用戶更大的使用空間。

［1］朱守輝.基于CA模型的土地利用變化研究——以新疆鄯善縣為例［J］.安徽農(nóng)業(yè)科學，2009，37(15):6887-6891.

［2］鄭占.基于CA模型的城市用地擴張模擬研究［D］.武漢:華中農(nóng)業(yè)大學，2010.

［3］徐昔保.基于GIS與元胞自動機的城市土地利用動態(tài)演化模擬與優(yōu)化研究——以蘭州市為例［D］.蘭州:蘭州大學，2007.

［4］江華麗，李華民.重慶市沙坪壩區(qū)土地集約利用評價研究［J］.安徽農(nóng)業(yè)科學，2008，36(24):10593-10595.

［5］周成虎，孫戰(zhàn)利，謝一春，等.地理元胞自動機研究［M］.北京:科學出版社，2001.

［6］李愛民，呂安民，隋春玲，等.集成GIS的元胞自動機在城市擴展模擬中的應用［J］.測繪科學技術學報，2009，26(3):165-169.

［7］管超，哈斯.基于GeoSOS的城市擴展模擬——以旅順口區(qū)為例［J］.測繪與空間地理信息，2014，37(3):49-54.

［8］岳彩榮，盧鵬，趙憲文，等.基于CA模型的香格里拉縣土地利用變化模擬［J］.安徽農(nóng)業(yè)科學，2011，39(4):2380-2384.

［9］黎夏，葉嘉安，劉小平，等.地理模擬系統(tǒng):元胞自動機與多智能體［M］.北京:科學出版社，2007:57-61.

［10］黃慶旭，何春陽，史培軍，等.城市擴展多尺度驅動機制分析——以北京為例［J］.經(jīng)濟地理，2009，29(5):46-49.