奉莉軍，鄭 琪，周家俊，郭 浩，梁勤歐

（浙江師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院，浙江 金華 321004）

研究未來(lái)土地利用變化對(duì)城市規(guī)劃和政策制定有重要意義[1]，其中，注重城市生態(tài)環(huán)境保護(hù)的土地利用格局尤為重要。模型模擬是研究土地利用變化的重要工具之一[2-3]，構(gòu)建土地利用模型有助于理解土地利用變化過(guò)程，解釋政策因素、人為因素以及自然因素等對(duì)未來(lái)土地利用的影響機(jī)制[4-5]。戴爾阜和馬良將土地利用模型分為6類(lèi)：機(jī)器學(xué)習(xí)和統(tǒng)計(jì)模型、經(jīng)濟(jì)學(xué)模型、元胞自動(dòng)機(jī)模型、多主體模型、部門(mén)經(jīng)濟(jì)學(xué)模型以及混合方法模型[7]。元胞自動(dòng)機(jī)（CA）模型能夠?qū)⒖臻g作用和時(shí)間關(guān)系相結(jié)合模擬復(fù)雜的時(shí)空演變過(guò)程，利用全局影響機(jī)制和局部轉(zhuǎn)移規(guī)則自上而下進(jìn)行模擬建模，從而能夠提取復(fù)雜的轉(zhuǎn)移規(guī)則進(jìn)行土地利用模擬[8-10]，被國(guó)內(nèi)外學(xué)者廣泛采用。然而，CA模型在土地利用模擬過(guò)程中側(cè)重于挖掘自然環(huán)境因素對(duì)土地變化的影響[11-12]，無(wú)法考慮土地利用變化過(guò)程中經(jīng)濟(jì)、人口等人為因素的影響[13-14]。基于神經(jīng)網(wǎng)絡(luò)方法的元胞自動(dòng)機(jī)（ANN-CA）模型利用神經(jīng)網(wǎng)絡(luò)算法提取土地利用變化的歷史規(guī)則[15-16]，在模擬過(guò)程中既能有效考慮自然因素又能考慮人為因素，能有效控制數(shù)據(jù)挖掘過(guò)程中所面臨的多因子干擾，提高模型模擬精度[17-18]。

在未來(lái)的中國(guó)城市群研究中，研究城市土地使用政策下的城市群空間變化尤其重要[19]。城市建設(shè)用地的擴(kuò)張速度不僅受到區(qū)域位置等自然因素的影響，還受到政策的調(diào)控。浙中城市群作為浙江省的第三大城市群，是浙江省參與長(zhǎng)江三角洲地區(qū)競(jìng)爭(zhēng)與合作的載體之一，城市化水平達(dá)到70%以上，其建設(shè)用地?cái)U(kuò)張速度受到浙中城市群規(guī)劃（2008—2020年）的土地利用政策的調(diào)控[20]。因此，只考慮歷史規(guī)則模擬未來(lái)土地利用情景不能科學(xué)解釋未來(lái)浙中城市群不同城市規(guī)劃下的土地空間利用，需要設(shè)置多種情景來(lái)解決這一問(wèn)題。本研究在分析2005—2019年浙中城市群建設(shè)用地?cái)U(kuò)張?zhí)卣鞯幕A(chǔ)上，以將浙中城市群建設(shè)為宜居城市為目標(biāo)，模擬浙中城市群建設(shè)用地在不同情景下的擴(kuò)張情況，以期為未來(lái)浙中城市群的城市規(guī)劃提供參考。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

研究區(qū)范圍包括浙江省金華市的整個(gè)行政區(qū)和衢州市的龍游、麗水市的縉云2個(gè)縣，總面積約為1.36萬(wàn)km2，如圖1。浙中城市群南部和北部地勢(shì)較高，中部地勢(shì)平坦且河流眾多，適合城市發(fā)展。作為浙江省第三大城市群，浙中城市群對(duì)浙江省中部地區(qū)城市建設(shè)高質(zhì)量發(fā)展具有引導(dǎo)作用。浙中城市群建設(shè)用地主要分布在中部地區(qū)，包括婺城區(qū)、金東區(qū)、義烏市、永康市以及東陽(yáng)市等地。

圖1 研究區(qū)位置、行政區(qū)劃和2019年浙中城市群建設(shè)用地Figure 1 Location of the study area，administrative divisions and the construction land of the mid-Zhejiang urban agglomeration in 2019

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

本研究采用歐洲航天局（https://earth.esa.int/）土地利用數(shù)據(jù)，經(jīng)過(guò)剪裁、重分類(lèi)等預(yù)處理得到浙中城市群2005、2010、2015和2019年的土地利用分類(lèi)數(shù)據(jù)，歸并為農(nóng)田、林地、草地、建設(shè)用地和水域5類(lèi)，空間分辨率為300 m×300 m。

影響土地利用的驅(qū)動(dòng)因素包括自然因素和社會(huì)因素，本研究選取包括自然因素和社會(huì)因素的9個(gè)驅(qū)動(dòng)因子來(lái)解釋土地利用的變化。地形因子由地理空間數(shù)據(jù)云（www.gscloud.cn）獲得，包括高程數(shù)據(jù)（DEM）以及由高程數(shù)據(jù)計(jì)算得到的坡度與坡向；社會(huì)因子直接由柵格數(shù)據(jù)形成，包括人口（https://www.worldpop.org/）和國(guó)民生產(chǎn)總值（https://www.ngdc.noaa.gov/）；區(qū)位因子為到城鎮(zhèn)中心的距離，由資源環(huán)境科學(xué)與數(shù)據(jù)中心（http://www.resdc.cn/）獲得；交通因子為距公路和鐵路的距離，與距水系距離一樣均通過(guò)對(duì)距離矢量數(shù)據(jù)進(jìn)行歐式距離運(yùn)算得到。對(duì)以上9個(gè)驅(qū)動(dòng)因子進(jìn)行數(shù)據(jù)預(yù)處理：先將數(shù)據(jù)重采樣到300 m，尺度與土地利用數(shù)據(jù)一致；再利用ArcGIS采用模糊隸屬度的方法進(jìn)行數(shù)據(jù)歸一化處理，消除量綱影響。

1.3 研究方法

1.3.1 GeoSOS系統(tǒng)與模擬評(píng)價(jià)

由中山大學(xué)黎夏教授等提出的地理模擬與優(yōu)化系統(tǒng)（GeoSOS）對(duì)土地利用變化歷史規(guī)則提取有較好的解釋?zhuān)荞詈系乩砟M和空間優(yōu)化對(duì)土地利用變化進(jìn)行提取[16]。GeoSOS 系統(tǒng)提供多種復(fù)合CA 模型，包括ANN-CA模型、DT-CA模型和Logistic-CA模型，本研究選取ANN-CA模型進(jìn)行浙中城市群的土地利用空間格局變化預(yù)測(cè)，利用神經(jīng)網(wǎng)絡(luò)的方法進(jìn)行歷史規(guī)則提取，按照規(guī)律計(jì)算各單元的轉(zhuǎn)移概率，模擬土地利用空間格局，預(yù)測(cè)未來(lái)的土地利用情景[16-18]。

首先按照百分比對(duì)研究區(qū)進(jìn)行采樣，設(shè)置鄰域參數(shù)，訓(xùn)練神經(jīng)網(wǎng)絡(luò)后提取網(wǎng)絡(luò)權(quán)重值；然后設(shè)置土地利用類(lèi)型轉(zhuǎn)移矩陣，判定土地類(lèi)型轉(zhuǎn)移的可轉(zhuǎn)移性；最后規(guī)定自適應(yīng)系數(shù)、迭代次數(shù)等參數(shù)計(jì)算轉(zhuǎn)移結(jié)果，當(dāng)達(dá)到某類(lèi)用地的實(shí)際需求后停止模擬。ANN-CA模型轉(zhuǎn)移規(guī)則計(jì)算公式[18]為

其中，P(k,t,l)為轉(zhuǎn)移概率，(1+(-lnγ)α)為不可控因素，Pann(k,t,l)為某地類(lèi)轉(zhuǎn)移概率，為城市建設(shè)用地密度，con()為兩種土地利用類(lèi)型之間轉(zhuǎn)移可能性。P(k,t,l)為1 說(shuō)明可轉(zhuǎn)移，P(k,t,l)為0 說(shuō)明不能發(fā)生轉(zhuǎn)移。

Kappa系數(shù)是衡量模型模擬精度的常用指標(biāo)。本研究利用Kappa系數(shù)進(jìn)行預(yù)測(cè)結(jié)果與實(shí)際分類(lèi)的一致性檢驗(yàn)，一般認(rèn)為數(shù)值大于0.75時(shí)即通過(guò)檢驗(yàn)，說(shuō)明模擬效果較好[21]。

1.3.2 情景構(gòu)建

根據(jù)浙中城市群的歷史發(fā)展規(guī)律，建立3種規(guī)劃情景，利用ANN-CA 模型設(shè)置規(guī)劃參數(shù)，對(duì)3種情景下土地利用格局進(jìn)行預(yù)測(cè)。

（1）慣性發(fā)展下的土地利用情景 基于浙中城市群土地利用歷史發(fā)展規(guī)律，本研究根據(jù)2005—2015年浙中城市群土地利用面積轉(zhuǎn)移矩陣設(shè)置土地利用類(lèi)型轉(zhuǎn)移矩陣，進(jìn)行慣性發(fā)展下的土地利用預(yù)測(cè)。表1是根據(jù)2005—2015 年浙中城市群土地利用面積轉(zhuǎn)移地類(lèi)變化情況建立的慣性情景下土地利用類(lèi)型轉(zhuǎn)移矩陣[22]。

表1 土地利用類(lèi)型轉(zhuǎn)移矩陣Table 1 Land use type transfer matrix

（2）經(jīng)濟(jì)優(yōu)先發(fā)展下的土地利用情景 此情景下的土地利用類(lèi)型轉(zhuǎn)移矩陣設(shè)置方法與慣性發(fā)展情景相似。在此情景下，制定城市用地規(guī)劃時(shí)，優(yōu)先考慮經(jīng)濟(jì)發(fā)展和城市建設(shè)，保證城鎮(zhèn)化發(fā)展速度最快，建設(shè)用地?cái)U(kuò)張面積快速增加。

（3）生態(tài)優(yōu)先發(fā)展下的土地利用情景 浙中城市群的建設(shè)目標(biāo)是打造重要生態(tài)屏障，促進(jìn)城鄉(xiāng)之間協(xié)調(diào)發(fā)展，爭(zhēng)取做生態(tài)城市示范區(qū)[23]。因此，該情景在慣性發(fā)展的前提下，著重保護(hù)耕地、林地和水域，提高森林覆蓋率，限制建設(shè)用地?cái)U(kuò)張速度，注重城市綠色發(fā)展。

2 結(jié)果與分析

2.1 浙中城市群城市擴(kuò)張變化特征

2.1.1 浙中城市群土地利用格局

研究區(qū)的土地類(lèi)型主要包括林地、耕地和建設(shè)用地，林地面積最大，耕地次之，建設(shè)用地處于快速擴(kuò)張狀態(tài)。如圖2，林地主要呈東北—西南走向，以條帶狀分布在浙中城市群地區(qū)，與耕地交錯(cuò)分布；耕地呈條帶狀分布于浙中城市群中部；草地零星分布于研究區(qū)內(nèi)；水域主要分布在耕地內(nèi)部，與耕地相互交錯(cuò)；建設(shè)用地主要分布于中部地區(qū)的婺城區(qū)、金東區(qū)、義烏市、永康市和東陽(yáng)市。

圖2 2005和2015 年浙中城市群土地利用格局Figure 2 Distribution of land use in the mid-Zhejiang urban agglomeration in 2005 and 2015

2005—2015年浙中城市群土地利用變化情況見(jiàn)表2，建設(shè)用地流入量最大，主要由耕地、草地和林地流入，其中耕地流入量高達(dá)389.67 km2；除建設(shè)用地和水域外，其他土地利用類(lèi)型面積均減少，2015年耕地面積比2005年減少了342.50 km2。

表2 2005—2015年浙中城市群土地利用面積轉(zhuǎn)移矩陣Table 2 Transfer matrix of the mid-Zhejiang urban agglomeration from 2005 to 2015km2

2.1.2 浙中城市群建設(shè)用地?cái)U(kuò)張?zhí)卣?/p>

由2005—2019年浙中城市群土地利用數(shù)據(jù)單獨(dú)提取建設(shè)用地?cái)?shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，得到浙中城市群建設(shè)用地2005—2019年面積變化以及增長(zhǎng)速度，結(jié)果如表3所示。2005—2019年，浙中城市群建設(shè)用地總面積持續(xù)增加，2019年建設(shè)用地總面積是2005年的2.78倍。浙中城市群建設(shè)用地增長(zhǎng)率不斷減小，從40.20%減小至14.13%。近年來(lái)，浙中城市群以建設(shè)宜居城市為目標(biāo)，注重生態(tài)環(huán)境保護(hù)和生態(tài)屏障建設(shè)，建設(shè)用地?cái)U(kuò)張速度逐漸下降。

表3 2005—2019年浙中城市群建設(shè)用地統(tǒng)計(jì)數(shù)據(jù)Table 3 Statistics on the construction land of the mid-Zhejiang urban agglomeration from 2005 to 2019

2.2 不同情景模型模擬精度檢驗(yàn)

本研究選取2005—2015年為時(shí)間步長(zhǎng)進(jìn)行ANN-CA 模型參數(shù)設(shè)置，提取土地利用數(shù)據(jù)總柵格數(shù)的5%為訓(xùn)練樣本，鄰域選擇為5時(shí)模擬結(jié)果較好，模擬精度較高[23]。對(duì)不同土地利用情景下模型進(jìn)行設(shè)置，模擬土地利用空間布局。模擬結(jié)果以2019年慣性發(fā)展情景為例展示（圖3），2019年慣性發(fā)展情景下的各類(lèi)土地利用空間分布格局與2019年真實(shí)土地利用空間分布格局基本一致，模擬效果較好。

圖3 2019年浙中城市群真實(shí)和慣性發(fā)展情景下土地利用對(duì)比圖Figure 3 Comparison of land use under real and inertial development scenari?os of the mid-Zhejiang urban agglomerations in 2019

3種情景下的模擬結(jié)果Kappa系數(shù)、FoM精度值和總體精度同時(shí)進(jìn)行檢驗(yàn)，結(jié)果如表4所示。慣性發(fā)展情景、經(jīng)濟(jì)優(yōu)先發(fā)展情景以及生態(tài)優(yōu)先發(fā)展情景下的Kappa系數(shù)均高于0.8，F(xiàn)oM精度值在0.04左右，總體精度達(dá)90%以上，均通過(guò)精度檢驗(yàn)。

表4 不同情景土地利用模擬精度比較Table 4 Comparison of simulation accuracy of land use in different scenarios

因本研究主要探究浙中城市群建設(shè)用地的擴(kuò)張情況，故對(duì)建設(shè)用地精度進(jìn)行單獨(dú)統(tǒng)計(jì)驗(yàn)證。3種情景下建設(shè)用地Kappa 系數(shù)見(jiàn)表5，均不低于0.75。其中，生態(tài)優(yōu)先情景下Kappa 系數(shù)最高，為0.78，該情景是最適合研究區(qū)建設(shè)用地模擬的情景?？傮w而言，精度檢驗(yàn)結(jié)果顯示ANN-CA模型適合浙中城市群建設(shè)用地?cái)U(kuò)張模擬，能夠滿(mǎn)足研究精度需求。

表5 3種情景下建設(shè)用地Kappa系數(shù)比較Table 5 Comparison of Kappa coefficients of construction land under three scenarios

2.3 2030年浙中城市群土地利用空間格局預(yù)測(cè)

3種情景下2030年浙中城市群的土地利用空間格局模擬結(jié)果如圖4所示。對(duì)比3種情景下2030年浙中城市群的土地利用格局發(fā)現(xiàn)，各情景下總體土地利用格局相似，生態(tài)優(yōu)先情景下建設(shè)用地面積低于其他2種情景。不同情景下2030 年浙中城市群各類(lèi)土地利用面積如表6 所示，生態(tài)優(yōu)先情景下的建設(shè)用地面積為904.70 km2，明顯低于其他2種情景；耕地面積為3 738.27 km2，明顯大于其他2種情景，達(dá)到了生態(tài)優(yōu)先發(fā)展的規(guī)劃要求，較好滿(mǎn)足了浙中城市群的建設(shè)目標(biāo)需要。因此，本研究選取生態(tài)優(yōu)先情景為最優(yōu)情景預(yù)測(cè)結(jié)果。

圖4 2030年浙中城市群3種情景土地利用格局模擬結(jié)果Figure 4 Simulated land use pattern of three scenarios in the mid-Zhejiang urban agglomeration in 2030

表6 不同情景下2030年浙中城市群土地利用面積Table 6 Land use area statistics of different scenarios in the mid-Zhejiang urban agglomeration in 2030 km2

生態(tài)優(yōu)先情景下，2030年浙中城市群建設(shè)用地的擴(kuò)張主要集中在婺城區(qū)、金東區(qū)、義烏市、永康市、東陽(yáng)市（表7）。上述地區(qū)是浙中城市群經(jīng)濟(jì)發(fā)展的主要區(qū)域，建設(shè)用地空間需求較大。2005年浙中城市群各市（縣、區(qū)）整體建設(shè)用地面積為321.76 km2，2015 年為721.34 km2。根據(jù)預(yù)測(cè)，2030 年在生態(tài)優(yōu)先情景下的浙中城市群各市（縣、區(qū)）整體建設(shè)用地面積將增加到904.70 km2，建設(shè)用地需求不斷增加。

表7 浙中城市群各市（縣、區(qū)）建設(shè)用地面積Table 7 Construction land area requirements of the counties in the mid-Zhejiang urban agglomerationkm2

如表8 所示，在生態(tài)優(yōu)先情景下，2015—2030 年浙中城市群建設(shè)用地?cái)U(kuò)張總量為183.86 km2，與2005—2015年擴(kuò)張總量399.58km2相比下降較多，擴(kuò)張速度逐漸減小。由于浙中城市群秉持生態(tài)優(yōu)先思想，建設(shè)生態(tài)廊道，在保護(hù)原有生態(tài)屏障的同時(shí)確保總體土地利用格局不發(fā)生變化。同時(shí)，有研究發(fā)現(xiàn)，城鎮(zhèn)化建設(shè)達(dá)到一定水平后城市建設(shè)速度會(huì)逐漸減緩，建設(shè)用地?cái)U(kuò)張速度下降會(huì)相對(duì)明顯[24]，浙中城市群建設(shè)用地?cái)U(kuò)張速度變化符合上述規(guī)律?？傮w來(lái)看，與2005—2015年相比，2015—2030年浙中城市群建設(shè)用地?cái)U(kuò)張面積明顯減少，與浙中城市群建設(shè)用地?cái)U(kuò)張?zhí)卣鳎ū?）一致。

表8 浙中城市群各市（縣、區(qū)）建設(shè)用地?cái)U(kuò)張面積Table 8 Expansion area of built-up land in various counties in the mid-Zhejiang urban agglomerationkm2

3 結(jié)論

（1）基于浙中城市群2005、2010、2015和2019年四期土地利用遙感影像和自然因子、社會(huì)因子數(shù)據(jù)，建立了基于慣性發(fā)展、經(jīng)濟(jì)優(yōu)先以及生態(tài)優(yōu)先3種不同發(fā)展情景；3種情景下2019年浙中城市群土地利用空間格局的模擬精度Kappa系數(shù)均達(dá)到0.8以上，模擬效果較好，ANN-CA模型適用于浙中城市群的土地利用空間模擬。

（2）經(jīng)模擬，2030年研究區(qū)在生態(tài)優(yōu)先情景下的建設(shè)用地面積需求均少于其他2種情景，土地利用空間格局最優(yōu)。因此，該情景下的土地利用格局最適合研究區(qū)的城市規(guī)劃發(fā)展。

（3）生態(tài)優(yōu)先情景下，2030年浙中城市群建設(shè)用地?cái)U(kuò)張以婺城區(qū)、金東區(qū)、義烏市、永康市和東陽(yáng)市為主導(dǎo)地區(qū)，其他市（縣、區(qū)）緊隨其后且均處于增長(zhǎng)狀態(tài)。2005—2030 年建設(shè)用地?cái)U(kuò)張?zhí)幱诔掷m(xù)增長(zhǎng)狀態(tài)，與2005—2015 年相比，2015—2030年建設(shè)用地?cái)U(kuò)張速度逐漸變緩，與現(xiàn)實(shí)中浙中城市群建設(shè)用地?cái)U(kuò)張?zhí)卣骰疽恢隆?/p>

本研究采用模型模擬方法進(jìn)行土地利用情景設(shè)置時(shí)，由于不同預(yù)測(cè)者經(jīng)驗(yàn)不同，考慮的政策等因素不同，預(yù)測(cè)的結(jié)果帶有一定的主觀性。本研究選取土地利用數(shù)據(jù)的空間分辨率為300 m×300 m，可能與高精度影像數(shù)據(jù)模擬結(jié)果存在偏差，未來(lái)可提高數(shù)據(jù)精度以進(jìn)一步提高模擬結(jié)果的可靠性[25]。