孟 成，潘琛玲，徐 寅

（1.華中科技大學(xué) 公共管理學(xué)院，湖北 武漢 430074；2.武漢市國(guó)土資源和規(guī)劃信息中心，湖北 武漢 430014）

用地布局模擬預(yù)測(cè)動(dòng)態(tài)性研究

孟 成1,2，潘琛玲2，徐 寅2

（1.華中科技大學(xué) 公共管理學(xué)院，湖北 武漢 430074；2.武漢市國(guó)土資源和規(guī)劃信息中心，湖北 武漢 430014）

分析了用地模擬預(yù)測(cè)中的空間因子動(dòng)態(tài)性，并通過(guò)主成分分析法將各空間因子變化量融合成新空間因子，再帶入預(yù)測(cè)模型對(duì)模擬結(jié)果進(jìn)行修正。以武漢市為例進(jìn)行了實(shí)例論證，通過(guò)該模型修正的預(yù)測(cè)結(jié)果和傳統(tǒng)靜態(tài)模型的預(yù)測(cè)結(jié)果的擬合度在80%以上，說(shuō)明該動(dòng)態(tài)修正模型具有較高的可靠性，能有效提高模擬預(yù)測(cè)工作效率、模擬精度以及預(yù)測(cè)成果的延用性。

用地規(guī)模；用地布局；模擬預(yù)測(cè)；動(dòng)態(tài)模型

改革開(kāi)放以來(lái)，我國(guó)進(jìn)入快速城鎮(zhèn)化時(shí)期，特別在2000年以后，我國(guó)城鎮(zhèn)面積的年均增長(zhǎng)率達(dá)到7%，隨著城市的快速發(fā)展，有限的土地資源與人類(lèi)對(duì)土地的無(wú)限需求、社會(huì)經(jīng)濟(jì)的快速發(fā)展以及人口不斷增加間的矛盾日益明顯?？茖W(xué)合理地模擬預(yù)測(cè)用地規(guī)模和用地布局有利于人民生活水平的穩(wěn)步提高和國(guó)民經(jīng)濟(jì)的持續(xù)發(fā)展，也能為合理地解決土地資源與人口、經(jīng)濟(jì)、社會(huì)等之間的矛盾提供技術(shù)方法[1]。

現(xiàn)行的用地布局模擬預(yù)測(cè)方法多屬于靜態(tài)預(yù)測(cè)，是在假設(shè)各作用力不變的情況下進(jìn)行的外推計(jì)算。當(dāng)有新數(shù)據(jù)加入樣本時(shí)，靜態(tài)預(yù)測(cè)模擬方法無(wú)法充分挖掘和利用原有預(yù)測(cè)成果與新樣本的聯(lián)系，需重新設(shè)定樣本進(jìn)行預(yù)測(cè)模擬工作，預(yù)測(cè)效率較低、結(jié)果時(shí)效性差，且中遠(yuǎn)期預(yù)測(cè)結(jié)果誤差大。目前，針對(duì)預(yù)測(cè)模型動(dòng)態(tài)修正性的研究主要是將統(tǒng)計(jì)學(xué)理論與不確定性分析、進(jìn)化算法等相結(jié)合，例如陳曉懷[2]等提出基于貝葉斯模型研究的動(dòng)態(tài)誤差實(shí)時(shí)修正方法，陳小祥[3]進(jìn)行的基于貝葉斯概率城市土地利用動(dòng)態(tài)變化模擬研究，以及基于傅里葉變換的質(zhì)量預(yù)測(cè)動(dòng)態(tài)殘差修正模型等[4-6]。這些研究主要集中在統(tǒng)計(jì)學(xué)模型方面，針對(duì)用地空間布局模擬預(yù)測(cè)動(dòng)態(tài)修正方法的研究較少[7-9]。基于此，本文設(shè)計(jì)了一種動(dòng)態(tài)修正的用地布局模擬預(yù)測(cè)方法，通過(guò)挖掘和利用新加入的樣本數(shù)據(jù)和同期預(yù)測(cè)值間的聯(lián)系對(duì)原預(yù)測(cè)結(jié)果進(jìn)行修正，從而減少了預(yù)測(cè)工作量，提高了預(yù)測(cè)效率和預(yù)測(cè)成果的延用性。

1 研究思路及方法

1.1 用地模擬預(yù)測(cè)動(dòng)態(tài)性分析

空間因子指對(duì)城市用地發(fā)展和變化有較大影響，具有空間特征的影響因子，主要包括地形地貌因子、區(qū)域因子、交通因子等。地形地貌因子主要包括城市坡度、坡向，山體等；區(qū)域因子主要包括城市中心點(diǎn)、城市行政區(qū)中心點(diǎn)等；交通因子主要包括城市環(huán)線、軌道交通線、城市主干路網(wǎng)、鐵路線、機(jī)場(chǎng)、交通換乘點(diǎn)等。

地形地貌因子是較為穩(wěn)定的一類(lèi)空間因子，城市坡度、坡向，山體等屬于城市的自然屬性，一般不會(huì)輕易發(fā)生變化。城市邊界、城市行政區(qū)劃等區(qū)域因子，在一定時(shí)期內(nèi)變化較小，也屬較為穩(wěn)定的一類(lèi)空間因子。與其他兩種因子相比，交通因子具有非穩(wěn)定性特點(diǎn)，隨著城市經(jīng)濟(jì)發(fā)展，人口不斷聚集，機(jī)動(dòng)車(chē)保有量不斷上升，城市的交通壓力逐步增大，城市交通也逐步從路面交通向地下、路面、地上3層的立體式交通轉(zhuǎn)變，地下通道、高架立交不斷增加，輕軌地鐵也進(jìn)入快速建設(shè)時(shí)期。因此，空間因子動(dòng)態(tài)性主要體現(xiàn)在交通因子的動(dòng)態(tài)性變化上。

1.2 用地布局模擬預(yù)測(cè)動(dòng)態(tài)修正模型

本文將城市中心點(diǎn)、行政區(qū)中心點(diǎn)、都市發(fā)展區(qū)中心點(diǎn)、環(huán)線、主干路網(wǎng)、交通路網(wǎng)、軌道交通線、鐵路線、高程9個(gè)因子作為空間因子進(jìn)行研究。首先，通過(guò)空間主成分分析對(duì)以上空間因子進(jìn)行融合。根據(jù)§1.1中分析可知，空間因子的動(dòng)態(tài)性主要集中在環(huán)線道路、主干交通路網(wǎng)和軌道交通等交通因子方面，因此，當(dāng)交通因子發(fā)生變化時(shí)，將各交通因子中新增的部分單獨(dú)提出并融合為一類(lèi)空間因子，記為PC_traffic，將PC_traffic作為新主成分帶入CA模型進(jìn)行運(yùn)算。具體研究流程如圖1所示。

設(shè)原交通因子為P1、P2、…、Pn；新交通因子為變化的交通因子為p1、p2、…、pn；原有主成分空間因子為PC_DEM、PC1、PC2、…、PCn。

圖1 空間因子動(dòng)態(tài)性修正流程圖

空間因子動(dòng)態(tài)修正模型的驗(yàn)證主要包括兩方面：

1）計(jì)算新融合的交通因子變化量PC_traffic與原主成分空間因子間的相關(guān)性，記為R。|R|介于0～1之間，通常|R|≥0.8表示高度相關(guān)；0.8＞|R|≥0.5表示中度相關(guān)；|R|≤0.5表示低度相關(guān)。本文中，若相關(guān)性較低，則說(shuō)明PC_traffic可作為新的主成分空間因子帶入模型計(jì)算。

2）將新交通因子與城市區(qū)域因子和城市地形地貌因子進(jìn)行主成分分析，融合為新主成分空間因子[10]。

3）將融合后得到的主成分空間因子帶入預(yù)測(cè)模型的模擬結(jié)果與將新融合的交通因子變化量PC_traffic與原有主成分空間因子作為變量帶入CA預(yù)測(cè)模型得到的結(jié)果進(jìn)行比較，計(jì)算兩種結(jié)果的相似度，以驗(yàn)證修正模型的可靠性。CA模型公式為：

4）評(píng)價(jià)模擬結(jié)果精度，通常采用Kappa系數(shù)表示，其計(jì)算公式為：

設(shè)柵格總像元數(shù)為n；真實(shí)柵格單元為1的像元數(shù)為a1，為0的像元數(shù)記為a0；模擬柵格單元記為1的像元數(shù)為b1，為0的像元數(shù)為b0；兩柵格對(duì)應(yīng)的相等像元數(shù)為s；則

Kappa計(jì)算結(jié)果通常是落在 0～1之間，可分5組表示不同級(jí)別的一致性：0～0.20時(shí)為極低一致性，0.21～0.40時(shí)為一般一致性，0.41～0.60時(shí)為中等一致性，0.61～0.80時(shí)為高度一致性，0.81～1時(shí)則幾乎為完全一致。

2 模型檢驗(yàn)與結(jié)果分析

2.1 研究區(qū)概況

武漢市地處江漢平原東部，長(zhǎng)江中下游地帶與漢水交匯處，總面積為8 569.36 km2，截至2012年，武漢市共有常住人口1 012萬(wàn)，國(guó)民生產(chǎn)總值8 003.82億 元（數(shù)據(jù)來(lái)源于《武漢市統(tǒng)計(jì)年鑒2013》）。本文依據(jù)2007年版《土地利用現(xiàn)狀分類(lèi)》對(duì)土地進(jìn)行分類(lèi)，根據(jù)研究需要將耕地、園地、林地和草地歸并為一類(lèi)，交通運(yùn)輸用地除街巷用地外作為一類(lèi)，水域及水利設(shè)施用地作為一類(lèi)，其他土地中除空閑地以外的部分作為一類(lèi)，一級(jí)地類(lèi)中的公共管理與公共服務(wù)用地、住宅用地、商服用地、工礦倉(cāng)儲(chǔ)用地、特殊用地及二級(jí)地類(lèi)中的街巷用地、空閑地歸并為一類(lèi)稱(chēng)為城鎮(zhèn)村及工礦用地，同時(shí)對(duì)該類(lèi)地進(jìn)行二級(jí)劃分，分為城市用地、建制鎮(zhèn)用地、村莊用地、采礦和風(fēng)景名勝用地及特殊用地5類(lèi)，如表1所示。

表1 武漢市2002～2012年各類(lèi)土地?cái)?shù)據(jù)表/km2

2.2 空間因子動(dòng)態(tài)修正模型論證

作為國(guó)家中心城市之一，近年來(lái)，武漢市正在大規(guī)模進(jìn)行城市建設(shè)。城市基礎(chǔ)道路、快速道路、地鐵等軌道交通設(shè)施處于快速建設(shè)階段，按照武漢軌道交通建設(shè)的近遠(yuǎn)期規(guī)劃，截至2020年主城區(qū)線網(wǎng)規(guī)模將達(dá)到333 km，建成7條過(guò)江通道。位于漢口地區(qū)的武漢軌道交通1號(hào)線的一期線路于2004-09-28投入運(yùn)營(yíng)，二期工程也于2010-07-29投入運(yùn)營(yíng)；武漢地鐵2 號(hào)線一期工程于2012-12-28投入運(yùn)營(yíng)；地鐵4號(hào)線一期工程于2013年12月28日投入運(yùn)營(yíng)；目前已形成了約70 km的“工”字形軌道路網(wǎng)骨架，武漢三鎮(zhèn)除漢陽(yáng)區(qū)外的7個(gè)主城區(qū)以及東西湖區(qū)實(shí)現(xiàn)軌道路網(wǎng)的連通。由于軌道交通1號(hào)線一期于2004年通車(chē)，二期于2009 年鋪軌，2010年正式通車(chē)；地鐵2號(hào)線于2006年底開(kāi)建，2007年7月全線開(kāi)工，2012年12月正式運(yùn)行。因此研究中以2008年、2009年、2010年、2011年和2012年的土地利用現(xiàn)狀數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)。

空間因子動(dòng)態(tài)修正模型論證思路如下：首先將軌道交通1號(hào)線，地鐵2號(hào)線連同其他空間因子進(jìn)行主成分分析，以2008年土地利用數(shù)據(jù)為基期數(shù)據(jù)，2009 年土地利用數(shù)據(jù)為末期數(shù)據(jù)，預(yù)測(cè)2012年各類(lèi)土地規(guī)模與分布，記為方案1；再將軌道交通1號(hào)線與其他空間因子進(jìn)行主成分分析，以地鐵2號(hào)線作為空間因子變量PC_traffic，以2008年土地利用數(shù)據(jù)為基期數(shù)據(jù)，2009年土地利用數(shù)據(jù)為末期數(shù)據(jù)，預(yù)測(cè)2012年各類(lèi)土地規(guī)模與分布，記為方案2；最后比較兩種預(yù)測(cè)結(jié)果，驗(yàn)證空間因子動(dòng)態(tài)修正模型的可靠性。方案1主成分空間因子和方案2新增主成分空間因子PC_traffic如圖2所示。

圖2 主成分空間因子圖集

分別計(jì)算PC_traffic與PC1、PC2、PC3、PC4、PC5之間的相關(guān)系數(shù)，以確定PC_traffic與原各主成分空間因子間的相關(guān)性。PC_traffic與PC1之間相關(guān)系數(shù)為0.59，與PC2之間相關(guān)系數(shù)為-0.43，與PC3之間相關(guān)系數(shù)為-0.16，與PC4之間相關(guān)系數(shù)為-0.47，與PC5之間相關(guān)系數(shù)為0.54?？芍?，PC_traffic與PC1和PC5正相關(guān)，與PC2、PC3、PC4負(fù)相關(guān)，且相關(guān)系數(shù)均小于0.6，因此PC_traffic與原主成分空間因子之間相關(guān)性較低。

方案1預(yù)測(cè)的2012年的各類(lèi)土地規(guī)模布局如圖3。

圖3 空間因子動(dòng)態(tài)性修正結(jié)果對(duì)比模擬圖

將方案2和方案1的預(yù)測(cè)結(jié)果進(jìn)行比較得到，城市用地相似度為93.62%，建制鎮(zhèn)用地相似度為98.63%，村莊用地相似度為93.50%，采礦用地相似度為96.93%，風(fēng)景名勝及特殊用地相似度為95.72%，交通運(yùn)輸用地相似度為81.81%，水域及水利設(shè)施用地相似度為98.74%，農(nóng)林地相似度為99.57%，其他用地面相似度為99.48%，兩方案的整體相似度為89.04%。綜上所述，經(jīng)空間因子動(dòng)態(tài)修正模型修正的模擬結(jié)果與采用原方法進(jìn)行模擬的結(jié)果之間具有較高的擬合度，證明該修正模型具有較高可靠性。

3 結(jié) 語(yǔ)

本文分析了用地布局模擬中的空間動(dòng)態(tài)性因素，提出了用地布局模擬預(yù)測(cè)動(dòng)態(tài)修正模型。經(jīng)過(guò)實(shí)例論證，該模型具有較高的可靠性，經(jīng)過(guò)修正的結(jié)果具有較高的精度，能替代基于新數(shù)據(jù)重新預(yù)測(cè)的結(jié)果。本文對(duì)用地布局模擬預(yù)測(cè)動(dòng)態(tài)修正方面進(jìn)行了初探。下一步可在此基礎(chǔ)上對(duì)國(guó)家、地方層面的相關(guān)政策、法規(guī)等控制因子的動(dòng)態(tài)性進(jìn)行研究；還可將本文的研究思路和方法應(yīng)用到人口、經(jīng)濟(jì)、交通等其他領(lǐng)域的模擬預(yù)測(cè)研究中。

[1] 李宗華,汪如民,童秋英.基于優(yōu)化組合模型的建設(shè)用地增長(zhǎng)趨勢(shì)預(yù)測(cè)[C]//中國(guó)地理信息系統(tǒng)協(xié)會(huì).中國(guó)地理信息產(chǎn)業(yè)論壇暨第二屆教育論壇就業(yè)洽談會(huì)論文集.北京:中國(guó)地理信息系統(tǒng)協(xié)會(huì),2009:7

[2] 陳曉懷,程真英,費(fèi)業(yè)泰.基于貝葉斯模型的動(dòng)態(tài)誤差實(shí)時(shí)修正方法研究[J].計(jì)量學(xué)報(bào),2008(增刊1):31-34

[3] 陳小祥.基于貝葉斯概率城市土地利用動(dòng)態(tài)變化模擬研究[D].長(zhǎng)沙:中南大學(xué),2007

[4] 甘麗珍,張銘鑫,劉明周,等.基于傅里葉變換的質(zhì)量預(yù)測(cè)動(dòng)態(tài)殘差修正模型[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,32(6):806-810

[5] 孟成,彭明軍,鄂棟臣.一種可動(dòng)態(tài)修正的建設(shè)用地預(yù)測(cè)方法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2014,39(1):95-99

[6] 王良健,師迎春,林目軒.BP神經(jīng)網(wǎng)絡(luò)結(jié)合小波處理在城市擴(kuò)張預(yù)測(cè)中的應(yīng)用:以長(zhǎng)沙市區(qū)為例[J].中國(guó)土地科學(xué),2008,22(1):39-47

[7] 郜允兵,高秉博,李曉嵐,等.長(zhǎng)時(shí)間序列土地利用數(shù)據(jù)時(shí)空索引技術(shù)研究[J].中國(guó)土地科學(xué),2015,29(9):34-41

[8] 馬世發(fā),蔡玉梅,念沛豪,等.土地利用規(guī)劃模型研究綜述[J].中國(guó)土地科學(xué),2014,28(3):34-38

[9] 李建春,張軍連,周楊,等.基于GIS的土地利用規(guī)劃空間布局合理性評(píng)價(jià)[J].中國(guó)土地科學(xué),2013,27(12):32-38

[10] 孟成,周鵬,李青.基于主成分分析的空間因子融合方法研究[C]//中國(guó)城市規(guī)劃協(xié)會(huì),武漢市國(guó)土規(guī)劃資源和規(guī)劃局組織編寫(xiě).城市規(guī)劃信息化.武漢:武漢出版社,2015:78-80

P208

B

1672-4623（2016）07-0006-03

10.3969/j.issn.1672-4623.2016.07.002

孟成，博士，主要研究方向?yàn)橥恋刭Y源管理和地理信息。

2016-01-22。

項(xiàng)目來(lái)源：教育部人文社會(huì)科學(xué)研究項(xiàng)目青年基金資助項(xiàng)目（15YJC630093）；中國(guó)博士后科學(xué)基金資助項(xiàng)目（2015M582229）。