土天時(shí)用總體規(guī)劃指標(biāo)時(shí)空分配研討——以濟(jì)南市長(zhǎng)清區(qū)為例

1、土天時(shí)用總體規(guī)劃指標(biāo)時(shí)空分配研討以濟(jì)南市長(zhǎng)清區(qū)為例1 1引言引言我國(guó)目前已開(kāi)展了兩輪土地利用總體規(guī)劃，目前正在進(jìn)行新一輪土地利用總體規(guī)劃修編工作。自 20 世紀(jì) 80 年代以來(lái)，土地利用規(guī)劃的理論和方法有了較大的發(fā)展1。近幾年，王萬(wàn)茂2、張鳳榮3、吳次芳4等人對(duì)土地利用規(guī)劃理論做了進(jìn)一步提升。隨著優(yōu)化技術(shù)的出現(xiàn)與應(yīng)用，各種線性和非線性規(guī)劃及多目標(biāo)規(guī)劃方法開(kāi)始應(yīng)用于土地利用規(guī)劃過(guò)程中，如層次分析法5、分形模型6、元胞自動(dòng)機(jī)7、理想點(diǎn)法8、生態(tài)足跡方法9,10等，諸多學(xué)者們對(duì)規(guī)劃用地指標(biāo)分解、土地利用類型變化分析、土地利用結(jié)構(gòu)分析和演變預(yù)測(cè)、土地利用規(guī)劃布局、土地利用規(guī)劃方案評(píng)價(jià)等進(jìn)行了探討研究。

2、已有研究雖然探討了各種新方法新技術(shù)在規(guī)劃中的應(yīng)用，但對(duì)規(guī)劃指標(biāo)的時(shí)空分配問(wèn)題涉及較少。而如何增強(qiáng)土地利用總體規(guī)劃時(shí)空分配的合理性和可操作性，正是土地利用總體規(guī)劃修編中必須面對(duì)的一個(gè)核心問(wèn)題和難點(diǎn)問(wèn)題，也是一個(gè)迫切需要突破的問(wèn)題。系統(tǒng)動(dòng)力學(xué)（System Dynamic, SD）模型能從宏觀上反映土地利用系統(tǒng)的結(jié)構(gòu)、功能和行為之間的相互關(guān)系，從而考察系統(tǒng)在不同情景下的變化和趨勢(shì)11,12，多運(yùn)用于區(qū)域時(shí)間尺度可持續(xù)發(fā)展研究中13,14,15,16,17，但是由于缺乏空間概念從而使結(jié)果很難在空間上直觀表達(dá)。而 CLUE-S （Conversionof Land Use and its Effect

3、s at Small regional extent）模型運(yùn)用系統(tǒng)論方法將土地利用系統(tǒng)看作一個(gè)整體，能同時(shí)模擬多種土地利用類型之間的競(jìng)爭(zhēng)關(guān)系并做不同的情景分析18,19，被證明其在模擬區(qū)域土地利用時(shí)空動(dòng)態(tài)變化方面具有明顯優(yōu)點(diǎn)20,21,22,23。但從已有的文獻(xiàn)中，融合上述二者優(yōu)點(diǎn)的研究很少見(jiàn)。那么，能否發(fā)揮 SD 模型對(duì)土地利用總體規(guī)劃指標(biāo)的時(shí)間模擬優(yōu)勢(shì)和 CLUE-S 模型在空間模擬方面的優(yōu)勢(shì)，從而解決土地利用總體規(guī)劃中一直困擾人們的“定量、定時(shí)、定位”問(wèn)題，便成為本文要解決的關(guān)鍵問(wèn)題。鑒于此，我們首先建立 SD 模型進(jìn)行土地利用總體規(guī)劃指標(biāo)的時(shí)間預(yù)測(cè)仿真模擬，然后在眾多影響土地利用的因素

4、中找到主控參量，并借助 GIS 使之空間化，隨后將 SD 預(yù)測(cè)指標(biāo)、主控參量集成到 CLUE-S 模型中進(jìn)行多方案模擬，經(jīng)過(guò)比對(duì)分析，并從中提出建議實(shí)施的規(guī)劃方案。2 2研究思路與技術(shù)路線研究思路與技術(shù)路線土地利用總體規(guī)劃指標(biāo)時(shí)空分配的基本思路是：定量、定時(shí)、定位。定量主要依據(jù)我國(guó)土地利用總體規(guī)劃編制思路，采取自上而下對(duì)總量指標(biāo)進(jìn)行分解；定時(shí)則是將總量指標(biāo)、分解指標(biāo)進(jìn)行不同時(shí)期的落實(shí)；定位則是將分解指標(biāo)具體落實(shí)到地類和地塊上。SD 的特征決定了它既可以解決定量問(wèn)題，也可以解決定時(shí)問(wèn)題，故文中借助 SD 模型解決定量和定時(shí)問(wèn)題。CLUE-S 模型通過(guò)綜合對(duì)土地利用的經(jīng)驗(yàn)分析、空間變異分析以及動(dòng)態(tài)

5、模擬來(lái)實(shí)現(xiàn)土地利用需求的空間分配，而這種分配總是和土地需求以及該地區(qū)的自然環(huán)境和社會(huì)經(jīng)濟(jì)狀況處在動(dòng)態(tài)平衡中。故文中將 SD 模型模擬結(jié)果作為土地需求，從導(dǎo)致土地利用需求變化的驅(qū)動(dòng)力中篩選自然和社會(huì)經(jīng)濟(jì)因素，并根據(jù)土地利用格局和備選驅(qū)動(dòng)因子，運(yùn)用邏輯斯蒂逐步回歸對(duì)每一個(gè)土地利用單元可能出現(xiàn)的各土地利用類型的概率進(jìn)行診斷。通過(guò) SD+CLUE-S 的數(shù)據(jù)融合實(shí)現(xiàn)土地利用總體規(guī)劃指標(biāo)的時(shí)空配置。使用 ROC（Relative Operating Characteristics）方法27對(duì)配置結(jié)果進(jìn)行檢驗(yàn)。隨 ROC 值的增加，回歸方程對(duì)土地利用空間分布格局的解釋能力逐漸增強(qiáng)。當(dāng) ROC0.7 時(shí)，可

6、認(rèn)為驅(qū)動(dòng)因子具有良好的解釋能力。3 3土地利用格局與驅(qū)動(dòng)因子分析土地利用格局與驅(qū)動(dòng)因子分析3.1 研究區(qū)概況及數(shù)據(jù)預(yù)處理濟(jì)南市長(zhǎng)清區(qū)處于北緯 3601至3732，東經(jīng) 11611至 11744之間，濟(jì)南市西南部，距濟(jì)南市區(qū)約22km，總面積 1178km2。2001 年 6 月撤縣設(shè)區(qū)。所用原始數(shù)據(jù)來(lái)源于長(zhǎng)清區(qū) 2006 年土地利用現(xiàn)狀變更數(shù)據(jù)（土地利用類型為過(guò)渡時(shí)期分類標(biāo)準(zhǔn)）。我們首先將土地利用類型在原有的二級(jí)十類的基礎(chǔ)上歸并為六類，以便進(jìn)行空間分析和下一步的時(shí)空分配模擬研究。這六類分別為耕地、園林地、牧草地及其他農(nóng)用地、建設(shè)用地、未利用土地和其他土地。然后分別從自然、社會(huì)、經(jīng)濟(jì)三個(gè)方面對(duì)長(zhǎng)

7、清區(qū)土地利用變化驅(qū)動(dòng)因子進(jìn)行分析。根據(jù)已有的相關(guān)數(shù)據(jù)和資料，共選取了 6 個(gè)驅(qū)動(dòng)因子城市化水平、高程、距城鎮(zhèn)距離、距主干道距離、坡度和人口密度。其中，城市化水平和人口密度數(shù)據(jù)以鄉(xiāng)鎮(zhèn)為單元進(jìn)行區(qū)域賦值，距城鎮(zhèn)距離和距主干道距離采用緩沖區(qū)進(jìn)行賦值，高程采用數(shù)字高程模型進(jìn)行區(qū)域化，坡度采用 GIS 軟件中的坡度生成工具進(jìn)行區(qū)域化。數(shù)據(jù)處理在 ArcView3.3 平臺(tái)上進(jìn)行，得到長(zhǎng)清區(qū) 2006 年土地利用現(xiàn)狀及各驅(qū)動(dòng)因子?xùn)鸥駡D（柵格大小為 250m250m）。3.2 回歸分析我們將 Grid 文件在 ArcView 下通過(guò) Export Data Source命令轉(zhuǎn)出為*.asc 文件，然后在 C

8、LUE-S 模型中進(jìn)行數(shù)據(jù)轉(zhuǎn)換。繼而在系統(tǒng)分析軟件 SPSS 下對(duì)每一種土地利用類型和備選驅(qū)動(dòng)因子，運(yùn)用 Logistic 回歸方程進(jìn)行分析和診斷，并對(duì)所得結(jié)果進(jìn)行 ROC 檢驗(yàn)?；貧w方程主要參數(shù)及 ROC 檢驗(yàn)值。4 4土地利用總體規(guī)劃指標(biāo)時(shí)空分配土地利用總體規(guī)劃指標(biāo)時(shí)空分配4.1 時(shí)間分配4.1.1 構(gòu)建 SD 時(shí)間分配模型我們?cè)O(shè)置土地利用總體規(guī)劃指標(biāo)時(shí)間分配模型的時(shí)間界限為 2006-2020 年，以研究區(qū) 1996-2006 年歷史數(shù)據(jù)作為設(shè)定各參數(shù)的主要依據(jù)，在 SD 模擬軟件 Powersim 下進(jìn)行仿真模擬。該模型系統(tǒng)是一個(gè)具有層次性的復(fù)雜開(kāi)放系統(tǒng)，包括人口、經(jīng)濟(jì)和土地利用三個(gè)子

9、系統(tǒng)。其中，人口子系統(tǒng)主要包括研究區(qū)總?cè)丝?、暫住人口、常住人口、城市化水平、城?zhèn)人口等變量；經(jīng)濟(jì)子系統(tǒng)主要包括國(guó)內(nèi)生產(chǎn)總值、建設(shè)用地對(duì)經(jīng)濟(jì)的拉動(dòng)系數(shù)、固定資產(chǎn)投資及固定資產(chǎn)土地增長(zhǎng)系數(shù)等變量；土地利用子系統(tǒng)主要包括各類用地面積構(gòu)成及其變化量、變化速率等變量。4.1.2 模型參數(shù)確定本文采用的系統(tǒng)參數(shù)確定方法主要有：（1）利用系統(tǒng)歷史統(tǒng)計(jì)資料作算數(shù)平均數(shù)，主要用于計(jì)算每年差異不大的用地類型的年均變化率，如牧草地增長(zhǎng)率、水利設(shè)施用地增長(zhǎng)率、未利用地減少率；（2）發(fā)展趨勢(shì)推算法，主要借助 SPSS 軟件對(duì)歷年數(shù)據(jù)進(jìn)行回歸分析，以建立合適的數(shù)學(xué)模型進(jìn)行預(yù)測(cè)，包括線性回歸預(yù)測(cè)（可用于計(jì)算固定資產(chǎn)投資土地

10、增長(zhǎng)系數(shù)、除以上三類用地之外的其他各類用地的年均變化率）和非線性回歸預(yù)測(cè)（可用于計(jì)算城市化水平、固定資產(chǎn)投資占 GDP 比重、城鎮(zhèn)人口）。4.1.3 模型有效性檢驗(yàn)我們以長(zhǎng)清區(qū) 1996 年的土地利用類型數(shù)據(jù)及人口社會(huì)經(jīng)濟(jì)數(shù)據(jù)為基數(shù)，通過(guò)設(shè)定變量參數(shù)，模擬 1997-2006 年的研究區(qū)各類土地的數(shù)量，發(fā)現(xiàn)農(nóng)用地、建設(shè)用地和未利用地的模擬結(jié)果與實(shí)際值之間的相對(duì)誤差均小于 1%，而農(nóng)用地中的耕地、園地、林地、其他農(nóng)用地的相對(duì)誤差最大為 2.67%；建設(shè)用地中的居民點(diǎn)及獨(dú)立工礦用地、交通用地、水利設(shè)施用地的相對(duì)誤差最大為 3.28%；常住人口、總?cè)丝?、?guó)民生產(chǎn)總值的相對(duì)誤差均小于1%。4.1.4

11、時(shí)間分配結(jié)果以 2006 年為基期年，2020 年為模擬終止年，步長(zhǎng)為5 年，得到 2006-2020 年規(guī)劃期間長(zhǎng)清區(qū)調(diào)整后六類用地的三個(gè)時(shí)間分配方案。4.2 空間分配4.1.1 分配原理空間分配動(dòng)態(tài)模擬是在綜合分析土地利用的空間分布概率適宜圖、土地利用變化規(guī)則和研究初期土地利用現(xiàn)狀圖的基礎(chǔ)上，根據(jù)相對(duì)適宜度（式）對(duì)土地利用需求進(jìn)行空間分配的過(guò)程。i u i u u u TPROP = P + ELAS + ITER , , 式中， i u TPROP, 柵格 i 適合土地利用類型 u 的相對(duì)適宜度， i u P, 由式求得， u ITER 土地利用類型 u 的迭代變量， u ELAS 根據(jù)

12、上述土地利用轉(zhuǎn)變規(guī)則設(shè)置的參數(shù)。4.1.2 模型參數(shù)確定4.1.2.1 土地利用類型轉(zhuǎn)換彈性系數(shù) CLUE-S 模型下轉(zhuǎn)換參數(shù) ELAS 的值介于 01 之間，參數(shù)值越大，土地利用就越穩(wěn)定，土地利用類型發(fā)生轉(zhuǎn)變的概率就越小22,28。經(jīng)過(guò)調(diào)試，本研究設(shè)定長(zhǎng)清區(qū)土地利用類型數(shù)量為 6，選取驅(qū)動(dòng)因子個(gè)數(shù)為 6，各地類回歸方程中入選驅(qū)動(dòng)因子最大個(gè)數(shù)為 5，柵格單元大小為 250m250m。各土地類型編碼為：0耕地，1園林地，2牧草地及其他農(nóng)用地，3建設(shè)用地，4未利用土地，5其他土地。它們的轉(zhuǎn)換彈性系數(shù)分別為 0.2，0.3，0.4，0.7，0.1，0.2。4.1.2.2 土地利用類型轉(zhuǎn)化矩陣土地類型

13、轉(zhuǎn)化矩陣行列數(shù)相等，均等于土地類型的數(shù)量。行表示現(xiàn)在的土地利用類型，列表示將來(lái)可能的土地類型，1表示現(xiàn)在和將來(lái)之間可以轉(zhuǎn)化，0 表示不能轉(zhuǎn)化。通過(guò)咨詢相關(guān)專家確定長(zhǎng)清區(qū)土地類型轉(zhuǎn)化矩陣。4.1.3 空間分配結(jié)果在 CLUE-S 模型下，根據(jù) Logistic 回歸方程診斷出的長(zhǎng)清區(qū) 2006 年各土地利用類型分布情況與其影響因素之間的相互關(guān)系，以及 2006 年土地利用現(xiàn)狀，結(jié)合該區(qū)未來(lái) 15 年的各土地利用類型數(shù)量需求及時(shí)間分配方案，根據(jù)上述迭代過(guò)程對(duì)其變化進(jìn)行空間定位分配，即可得到長(zhǎng)清區(qū)2006-2020 年規(guī)劃期間三個(gè)方案下的土地利用總體規(guī)劃指標(biāo)空間分配結(jié)果。5 5分析與討論分析與討論5

14、.1 指標(biāo)時(shí)間變化分析從 SD 模型時(shí)間模擬三個(gè)方案預(yù)測(cè)結(jié)果來(lái)看，長(zhǎng)清區(qū)未來(lái) 15 年內(nèi)，土地利用結(jié)構(gòu)中農(nóng)用地、建設(shè)用地、未利用地三大用地類型總體變化趨勢(shì)與 1996-2006 年期間的變化趨勢(shì)大體一致。伴隨著農(nóng)用地和未利用地面積的減少，建設(shè)用地面積增加?？紤]到今后城市的發(fā)展和人口的增長(zhǎng)，耕地后備資源會(huì)逐年減少。5.2 指標(biāo)空間變化分析將 2006 年土地利用現(xiàn)狀圖分別與三個(gè)方案下 2020年規(guī)劃指標(biāo)空間分配模擬圖進(jìn)行比較，可以發(fā)現(xiàn)六類土地中耕地、園林地、建設(shè)用地、未利用土地的空間變化比較大。總的來(lái)說(shuō)，耕地轉(zhuǎn)出主要集中在城鎮(zhèn)建設(shè)用地周邊地區(qū)及南部坡度較大和地勢(shì)較高的地區(qū)。轉(zhuǎn)出的耕地一部分為今后

15、城鎮(zhèn)建設(shè)的發(fā)展擴(kuò)張所占用，另一部分因?yàn)樽匀坏乩項(xiàng)l件不適宜農(nóng)業(yè)耕作，需退耕轉(zhuǎn)變?yōu)閳@林地或者未利用土地。而轉(zhuǎn)入的耕地大部分集中在北部坡度小、地勢(shì)低的地區(qū)，且以園林地和未利用土地為主要轉(zhuǎn)入地類，主要是因?yàn)橄啾容^南部丘陵山地而言，這些地區(qū)自然地理?xiàng)l件更適宜作為耕地。轉(zhuǎn)出的園林地主要分布在現(xiàn)有建設(shè)用地的周邊地區(qū)，主要有兩個(gè)去向：一是城鎮(zhèn)、交通用地等建設(shè)用地周?chē)霓D(zhuǎn)出園林地今后將轉(zhuǎn)化為建設(shè)用地；二是距離城鎮(zhèn)等建設(shè)用地較遠(yuǎn)但是自然條件適宜耕作的那部分將轉(zhuǎn)化為耕地。轉(zhuǎn)入的園林地主要集中在南部丘陵地區(qū)。園林地的空間分布變化基本符合土地利用結(jié)構(gòu)調(diào)整的原則。建設(shè)用地主要由耕地、園林地及未利用土地等其他用地類型轉(zhuǎn)入，

16、且轉(zhuǎn)入的建設(shè)用地多分布于城市或建制鎮(zhèn)周邊及原交通主干道附近。而未利用土地主要被開(kāi)發(fā)為耕地，被建設(shè)用地?cái)D占，或轉(zhuǎn)化為園林用地。從模擬的結(jié)果看，長(zhǎng)清區(qū)牧草地及其他農(nóng)用地和其他土地的空間分布變化情況不明顯?？傮w來(lái)說(shuō)，比較明顯的變化是建設(shè)用地的擴(kuò)張和耕地的減少，并且建設(shè)用地的擴(kuò)張多是在原有城鎮(zhèn)建設(shè)用地的基礎(chǔ)上向四周擴(kuò)展，侵占的多是其周邊的耕地、園林地及未利用土地。而三個(gè)方案空間分配模擬結(jié)果的差異主要體現(xiàn)在用地?cái)U(kuò)張或轉(zhuǎn)出的幅度上，圖 10 為三個(gè)方案下建設(shè)用地?cái)U(kuò)展對(duì)比圖。最后經(jīng)過(guò)比對(duì)，將方案二作為建議實(shí)施方案，已應(yīng)用于濟(jì)南市土地利用總體規(guī)劃修編中。6 6結(jié)論結(jié)論長(zhǎng)期以來(lái)，制約我國(guó)土地利用總體規(guī)劃編制的一

17、個(gè)重要因素，就是從技術(shù)上沒(méi)有真正解決“定量、定時(shí)、定位”問(wèn)題。SD 模型可以全面考慮驅(qū)動(dòng)因子和影響因素，全面反映土地利用系統(tǒng)結(jié)構(gòu)、功能和行為之間的相互關(guān)系，但其缺乏對(duì)空間因素的處理能力。而 CLUE-S 模型可以實(shí)現(xiàn)土地利用變化過(guò)程和結(jié)果在空間上的表達(dá)，具有更高的可信度和更強(qiáng)的解釋能力。本研究將 SD 和 CLUE-S模型相結(jié)合，發(fā)揮其各自優(yōu)勢(shì)，突破了一直困擾我國(guó)土地利用總體規(guī)劃編制的一個(gè)技術(shù)“瓶頸”。并通過(guò)濟(jì)南市長(zhǎng)清區(qū)的實(shí)驗(yàn)取得成功，為縣級(jí)土地利用總體規(guī)劃修編提供了一個(gè)范例。通過(guò)本文案例研究發(fā)現(xiàn)，將 SD 和 CLUE-S 模型進(jìn)行融合應(yīng)用，不僅可以解決土地利用總體規(guī)劃指標(biāo)的時(shí)空分配問(wèn)題，而

18、且還可以得到不同土地利用類型的具體變化情況，如各種土地利用類型的轉(zhuǎn)入-轉(zhuǎn)出分析等，通過(guò)不同時(shí)期的斷面分析，把握該地區(qū)未來(lái)土地利用格局的變化情況。同時(shí)，借助本文的研究思路，可以動(dòng)態(tài)地進(jìn)行多方案的模擬，使每個(gè)方案的數(shù)量在空間上可視化，提高了規(guī)劃方案的優(yōu)選和決策輔助。研究也發(fā)現(xiàn)，由于 CLUE-S 模型是基于柵格數(shù)據(jù)結(jié)構(gòu)進(jìn)行分析的，其在進(jìn)行后期的規(guī)劃管理時(shí)不夠靈活，制作的規(guī)劃圖也不夠美觀。十年論文機(jī)構(gòu)京都名師論文中心，正規(guī)全面的論文刊物為您提供職稱論文，畢業(yè)論文，碩士論文，醫(yī)學(xué)論文，教育論文等各類論文發(fā)表服務(wù)。參考文獻(xiàn)參考文獻(xiàn)1 孟曉晨，趙星爍.中國(guó)土地利用總體規(guī)劃實(shí)施中主要問(wèn)題及成因分析J. 中國(guó)

19、土地科學(xué)，2007，21(3):19-25.2 王萬(wàn)茂. 土地利用規(guī)劃學(xué)M. 北京：科學(xué)出版社，2006.3 張鳳榮. 持續(xù)土地利用管理的理論與實(shí)踐M. 北京：北京大學(xué)出版社，1996.4 吳次芳. 土地利用規(guī)劃M. 北京：地質(zhì)出版社，2000.5 堯德明，陳玉福，張福剛，等. 層次分析法在土地利用總體規(guī)劃用地指標(biāo)分解中的應(yīng)用以海南省為例J. 安徽農(nóng)業(yè)科學(xué)，2007，35(34):11175-11176,11178.6 沈非，查良松，李軍利，等. 基于分?jǐn)?shù)維和穩(wěn)定性指數(shù)的蘇州市土地利用類型變化分析J. 土壤，2007，39(6): 964-967.7 楊小雄，劉耀林，王曉紅，等. 基于約束條件

20、的元胞自動(dòng)機(jī)土地利用規(guī)劃布局模型J. 武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2007，32(12): 1164-1167.8 湯江龍，趙小敏，師學(xué)義. 理想點(diǎn)法在土地利用規(guī)劃方案評(píng)價(jià)中的應(yīng)用J. 農(nóng)業(yè)工程學(xué)報(bào)，2005，21(2):56-59.9 賴力，黃賢金. 全國(guó)土地利用總體規(guī)劃目標(biāo)的生態(tài)足跡評(píng)價(jià)研究J.農(nóng)業(yè)工程學(xué)報(bào)，2005，21(2): 66-71.10 符海月，李滿春，毛亮，等. 基于生態(tài)足跡的土地利用規(guī)劃生態(tài)成效定量分析以河北省廊坊市為例J. 自然資源學(xué)報(bào)，2007，22(2): 225-235.11 Mohapatra P, Mandal P, and Bora M. Introducti

21、on to systemdynamics modeling. Hyderabad, India: OrientLongman Ltd, 1994.12 王愛(ài)萍，鄭新奇. 基于 SD 的土地利用總體規(guī)劃J. 山東師大學(xué)報(bào)(自然科學(xué)版)，2001，16(3): 438-443.13 王其藩. 系統(tǒng)動(dòng)力學(xué)M. 北京：清華大學(xué)出版社，1998.14 李景剛，何春陽(yáng)，史培軍，等. 近 20 年中國(guó)北方 13 省耕地變化與驅(qū)動(dòng)力J. 地理學(xué)報(bào)，2004，59(2):274-282.15 Saysel A K, Barlas Y, Yemgun O. Long term sustainability ina

22、n agricultural development project: a systemdynamics approachJ.Journal of Environmental Management, 2002, 64: 247-260.16 孫新新，沈冰，于俊麗，等. 寶雞市水資源承載力系統(tǒng)動(dòng)力學(xué)仿真模型研究J. 西安建筑科技大學(xué)學(xué)報(bào)（自然科學(xué)版），2007，39(1): 72-77.17 程葉青，李同升，張平宇. SD 模型在區(qū)域可持續(xù)發(fā)展規(guī)劃中的應(yīng)用J. 系統(tǒng)工程理論與實(shí)踐，2004，12: 13-18.18 彭建，蔡運(yùn)龍，Verburg P H. 喀斯特山區(qū)土地利用/覆被變化情景模擬J.

23、 農(nóng)業(yè)工程學(xué)報(bào)，2007，23(7):64-73.19 盛晟，劉茂松，徐馳，等. CLUE-S 模型在南京市土地利用變化研究中的應(yīng)用J. 生態(tài)學(xué)雜志，2008, 27(2): 235-239.20 Veldkamp A and Verburg P H. Methodology for simulating thespatial dynamics of land use change in forestfringes. 2002.21張永民，趙士洞，Verburg P H. 科爾沁沙地及其周?chē)貐^(qū)土地利用變化的情景分析J. 自然資源學(xué)報(bào)，2004，19(1): 29-38.22 蔡玉梅，劉彥隨，宇振榮，等. 土地利用變化空間模擬的進(jìn)展CLUE-S 模型及