Krugman模型在城市生長(zhǎng)模擬研究中的應(yīng)用

杜艷琴，武澤江，張 梁

（1.陜西華地勘察設(shè)計(jì)咨詢有限公司，陜西 西安 710020）

Krugman模型在城市生長(zhǎng)模擬研究中的應(yīng)用

杜艷琴1，武澤江1，張 梁1

（1.陜西華地勘察設(shè)計(jì)咨詢有限公司，陜西 西安 710020）

基于Krugman理論模型，利用GIS技術(shù)，從地理和經(jīng)濟(jì)兩個(gè)視角，探索城市體系的演變過(guò)程，并分析城市數(shù)量和分布的理論規(guī)律。結(jié)果顯示，在地理和經(jīng)濟(jì)的非均質(zhì)空間下，城市數(shù)量與制造業(yè)勞動(dòng)力份額、需求價(jià)格彈性和可達(dá)性都有很大的相關(guān)性，且城市易于在人口比較密集的地區(qū)形成。

Krugman模型；GIS；非均質(zhì)空間；城市數(shù)量和分布

城市生長(zhǎng)是一個(gè)比較復(fù)雜的時(shí)空動(dòng)態(tài)演化過(guò)程，其主要表現(xiàn)在城市空間的演化上。Krugman模型作為當(dāng)代研究城市生長(zhǎng)的一種理論與方法，用于探索城市的根源及城市分布狀況。盡管該模型立意很高，但由于模型假設(shè)條件和模擬手段等原因，一直沒(méi)有被真正用于實(shí)證研究中。隨著Stelder引入“地理空間”對(duì)該模型進(jìn)行了一系列的改造[1]，該模型才從理論走向?qū)嵶C應(yīng)用。由于Stelder是經(jīng)濟(jì)學(xué)方面的專家，因此他在GIS環(huán)境下演繹的歐洲城市空間格局分布顯得比較淺顯，將GIS與城市生長(zhǎng)模型進(jìn)行集成有待進(jìn)一步研究。

1 研究?jī)?nèi)容與技術(shù)路線

本文以Krugman理論模型為基礎(chǔ)，采用GIS技術(shù)，構(gòu)造城市體系演繹的實(shí)驗(yàn)?zāi)Ｐ?，通過(guò)改變地理和經(jīng)濟(jì)參數(shù)，實(shí)現(xiàn)城市體系的動(dòng)態(tài)模擬和分析，獲得城市數(shù)量和分布的規(guī)律。技術(shù)路線如圖1所示。

2 構(gòu)建模型分析框架

2.1 空間數(shù)據(jù)集獲取

在進(jìn)行城市生長(zhǎng)模擬前，首先要解決虛擬城市區(qū)位的構(gòu)建問(wèn)題。本文沿用Stelder的方法，在GIS環(huán)境下，采用空間劃分技術(shù)，根據(jù)研究區(qū)域的起始位置及步長(zhǎng)，生成一組長(zhǎng)、寬均相等的格網(wǎng)點(diǎn)。每個(gè)格網(wǎng)點(diǎn)代表一個(gè)潛在的城市，它們均有可能演變成為未來(lái)的城市。構(gòu)造虛擬城市點(diǎn)位及空間網(wǎng)絡(luò)的流程如圖2所示。

圖2 生成虛擬城市點(diǎn)位流程圖

2.2 非均質(zhì)空間構(gòu)造

Krugman模型缺乏實(shí)證性在一定程度上是由它的均質(zhì)空間假設(shè)造成的，即沒(méi)有將空間的非均質(zhì)性納入到理論模型中加以分析。因此，為了使Krugman模型走向現(xiàn)實(shí)，需要假設(shè)非均質(zhì)空間來(lái)實(shí)現(xiàn)城市生長(zhǎng)模擬，進(jìn)而探索城市體系的演變過(guò)程。本文主要從地理參數(shù)與經(jīng)濟(jì)參數(shù)兩個(gè)方面對(duì)地理空間和經(jīng)濟(jì)空間進(jìn)行改造，生成的非均質(zhì)空間作為Krugman城市生長(zhǎng)模型的初始空間假設(shè)。

1）地理非均質(zhì)空間。在城市形成理論中，通過(guò)構(gòu)建經(jīng)濟(jì)學(xué)模型來(lái)揭示規(guī)模經(jīng)濟(jì)和需求外部性的重要性。模型中假設(shè)城市區(qū)位是同質(zhì)的，城市的形成由規(guī)模經(jīng)濟(jì)和交通成本決定，地理因子在城市演變中具有一定的影響[2]。目前，用分形的方法來(lái)研究地表面起伏是近年來(lái)最為流行的應(yīng)用領(lǐng)域之一。李后強(qiáng)、劉繼生、陳彥光等在分形方法的理論探討和區(qū)域驗(yàn)證中都做了深入探討[3-5]。構(gòu)造分形地形的方法很多，最為常用的是基于格網(wǎng)的隨機(jī)中點(diǎn)位移法和基于三角形的隨機(jī)中點(diǎn)位移法。前者適用于具有規(guī)則格網(wǎng)的地形數(shù)據(jù)，由于本文中數(shù)據(jù)是規(guī)則的，因此采用后者建立分形地形表面。

2）經(jīng)濟(jì)非均質(zhì)空間。勞動(dòng)力在城市演變中起著舉足輕重的作用,在一定的地域環(huán)境下，勞動(dòng)力呈現(xiàn)一定的空間分布形態(tài)?？臻g自相關(guān)是空間場(chǎng)中數(shù)值聚集程度的一種度量，是認(rèn)識(shí)空間分布特征的一種最常用的方法[6]。通過(guò)對(duì)城市勞動(dòng)力分布的空間相關(guān)性進(jìn)行分析，能揭示出城市勞動(dòng)力的空間分布形態(tài)。因此，本文構(gòu)造出符合空間自相關(guān)的數(shù)據(jù)作為勞動(dòng)力的分布。將勞動(dòng)力空間分布作為城市點(diǎn)位的屬性字段Z值，構(gòu)造人口密度圖。鑒于本文城市點(diǎn)位為均勻分布，且勞動(dòng)力存在空間相關(guān)性。因此，選用Kriging插值方法來(lái)完成柵格面的構(gòu)建。

2.3 Krugman模型

Krugman模型包含短期均衡求解和長(zhǎng)期均衡求解兩個(gè)求解過(guò)程。短期均衡是指在較短時(shí)期和區(qū)域內(nèi)的制造業(yè)勞動(dòng)力不可移動(dòng)的前提下，每個(gè)區(qū)域制造業(yè)的供求平衡，生產(chǎn)廠商沒(méi)有額外的利潤(rùn)可以賺取?！爸圃鞓I(yè)勞動(dòng)力不可移動(dòng)”就意味著制造業(yè)勞動(dòng)力和生產(chǎn)廠商的分布是不能通過(guò)模型計(jì)算而得到，必須事先給定一系列制造業(yè)勞動(dòng)力分布的初始值。

在給定的規(guī)模經(jīng)濟(jì)、不完全競(jìng)爭(zhēng)、位置、外部經(jīng)濟(jì)以及勞動(dòng)力不能流動(dòng)的前提下，短期均衡公式為：

其中，每個(gè)變量的含義如表1所示。

表1 變量含義

2.4 城市生長(zhǎng)模擬實(shí)現(xiàn)

本文采用完全集成原理，即在GIS軟件中，用Python腳本語(yǔ)言實(shí)現(xiàn)整個(gè)城市生長(zhǎng)模擬研究。采用該方法可以直接利用GIS軟件的所有功能，無(wú)需額外編寫此類程序，省時(shí)省力；缺點(diǎn)是運(yùn)行效率較低，無(wú)法完成更為復(fù)雜的擴(kuò)展。

選定時(shí)空數(shù)據(jù)庫(kù)的存儲(chǔ)方式。時(shí)空數(shù)據(jù)庫(kù)的語(yǔ)義定義，是當(dāng)今GIS領(lǐng)域的研究熱點(diǎn)。在所有的時(shí)空數(shù)據(jù)模型中，時(shí)間片快照模型數(shù)據(jù)結(jié)構(gòu)比較簡(jiǎn)單，容易在現(xiàn)有的GIS軟件內(nèi)實(shí)現(xiàn)。

基于Krugman理論模型，綜合短期均衡和長(zhǎng)期均衡兩個(gè)求解過(guò)程，得到城市演繹模擬的計(jì)算機(jī)流程如圖3所示。

圖3 城市生長(zhǎng)模擬流程圖

在進(jìn)行城市演繹模擬前，有幾點(diǎn)需要注意：

1）必須明確所要求解的變量。短期均衡可以確定各區(qū)域的收入Yr、價(jià)格指數(shù)Ir和工資率Wr。也就是說(shuō)，要找出一組數(shù)值，使得式（1）、（2）、（3）同時(shí)成立。而長(zhǎng)期均衡則可以得到各個(gè)區(qū)域內(nèi)制造業(yè)勞動(dòng)力的分布狀況。

2）在求解之前需要給定自變量λ0以及變量L、α、β、γ、δ、φ、ρ、ε、τ，參數(shù)的初始值不同，會(huì)導(dǎo)致最終產(chǎn)生的結(jié)果有差異。

3）求解該模型時(shí)用的是順序迭代法。具體步驟為：①給定假設(shè)區(qū)域內(nèi)工資率的初始值（W1,0，W2,0，…，Wr,0），其中0表示第0次迭代。②將參數(shù)（W1,0，W2,0， …， Wr,0）代入式（1）、（2），計(jì)算收入（Y1,0，Y2,0，…， Yr,0）和價(jià)格指數(shù)（I1,0，I2,0，…， Ir,0）。③將前一步中求出的收入和價(jià)格指數(shù)代入式（2）、（3），再求出工資率（W1,1，W2,1，… ，Wr,1）。④重復(fù)第②、③步，直至，其中，σ=0.001；如果無(wú)法滿足該收斂條件，則循環(huán)200次結(jié)束。⑤根據(jù)公式ωr=WrIr-δ和計(jì)算出參數(shù)（w1，w2，…，wr）和平均工資ω。⑥將參數(shù)（w1，w2，…，wr）和平均工資代入公式dλr=λr+ n（ωr-）中， 可求出制造業(yè)勞動(dòng)力分布（λ1，λ2，…，λr），直至當(dāng)?shù)冖懿胶偷冖薏街袃蓚€(gè)判斷條件均滿足時(shí)，可以得到在均衡狀態(tài)下，各區(qū)域勞動(dòng)力λr的分布，其中勞動(dòng)力的聚集地為城市。

4）循環(huán)終止條件為σ = 0．001，只有當(dāng)上述步驟對(duì)所有區(qū)域均滿足時(shí)，求解過(guò)程才可以停止。

3 非均質(zhì)空間下城市生長(zhǎng)分析

3.1 地理非均質(zhì)空間

在地理非均質(zhì)空間網(wǎng)絡(luò)下，通過(guò)調(diào)整制造業(yè)工人的勞動(dòng)力份額μ、可達(dá)性參數(shù)τ和需求價(jià)格彈性指數(shù)ε，驗(yàn)證Krugman理論模型，探討地理非均質(zhì)空間對(duì)城市演變的影響。

令ε = 5、μ = 0.7、τ在（0,1）區(qū)間，從小到大調(diào)整τ的值，結(jié)果如圖4。城市的生長(zhǎng)經(jīng)歷了一個(gè)從無(wú)到有的過(guò)程，隨著τ不斷增大，可達(dá)性逐漸降低。于是，制造業(yè)經(jīng)濟(jì)活動(dòng)的分散趨勢(shì)開(kāi)始逐漸增加，城市的模擬數(shù)目也相應(yīng)增加，即城市數(shù)量與可達(dá)性成反比。將ε和μ取不同的值，τ的取值范圍不變，進(jìn)行多次實(shí)驗(yàn)如圖5。當(dāng)ε和μ的值都較小時(shí)，可達(dá)性需減小到一定程度才能形成較多的城市，這表明城市數(shù)量與可達(dá)性為反比的關(guān)系，即可達(dá)性越高，城市數(shù)量越少；反之，城市數(shù)量越多，這明顯反映了交通條件對(duì)城市空間格局演變起到了推動(dòng)作用。

圖4 可達(dá)性與城市空間分布圖

圖5 各參數(shù)與城市數(shù)量關(guān)系圖

為了使結(jié)果更清晰，對(duì)區(qū)域內(nèi)的高差進(jìn)行統(tǒng)計(jì)并分類，將其劃分為平坦、較平坦、中等、較陡、陡峭5個(gè)等級(jí)。對(duì)城市出現(xiàn)的概率與高差分布進(jìn)行相關(guān)性驗(yàn)證，結(jié)果如表2、圖6。

表2 高差與城市出現(xiàn)概率相關(guān)系數(shù)

圖6 高差與城市出現(xiàn)概率趨勢(shì)圖

結(jié)果顯示，城市出現(xiàn)概率與高差呈現(xiàn)負(fù)相關(guān)，即越平坦的地區(qū)越易形成城市，而陡峭的地區(qū)，則不易形成城市。具體講，按照城市發(fā)展理論，地理優(yōu)勢(shì)在城市的形成和發(fā)展過(guò)程中發(fā)揮著相當(dāng)重要的作用。與此同時(shí)，城市的分布隨著時(shí)間的變化而變化，而地形是度量地理優(yōu)勢(shì)的關(guān)鍵所在。因此，地形在城市的形成和發(fā)展中起著重要的作用。城市不可能形成在山體太高，坡度太陡的地方，即越平坦的區(qū)域，越利于城市的發(fā)展。

3.2 經(jīng)濟(jì)非均質(zhì)空間

Krugman曾對(duì)城市的形成作出解釋，他指出人口集聚所帶來(lái)的收益遞增是城市形成的重要原因。因此，本文在經(jīng)濟(jì)非均質(zhì)空間下，通過(guò)調(diào)整制造業(yè)勞動(dòng)力份額，需求價(jià)格彈性指數(shù)和可達(dá)性等參數(shù)，對(duì)城市生長(zhǎng)模擬進(jìn)行一系列實(shí)驗(yàn)，分析經(jīng)濟(jì)參數(shù)對(duì)城市演變的影響，模擬過(guò)程同§3．1。

結(jié)果顯示，在經(jīng)濟(jì)非均質(zhì)空間時(shí)，城市數(shù)量隨著制造業(yè)勞動(dòng)力份額的增加而增加，即城市的數(shù)量與制造業(yè)動(dòng)力份額、需求價(jià)格彈性指數(shù)成正比，與可達(dá)性成反比；由于經(jīng)濟(jì)空間的非均質(zhì)性，城市分布均不再呈現(xiàn)完全對(duì)稱狀態(tài)。

將勞動(dòng)力初始空間分布劃分為5個(gè)區(qū)間，從小到大依次為：稀疏、較稀疏、中等、密集、較密集。對(duì)模擬結(jié)果進(jìn)行統(tǒng)計(jì)如表3，然后對(duì)統(tǒng)計(jì)結(jié)果進(jìn)行相關(guān)性驗(yàn)證，所得相關(guān)系數(shù)接近于1，如表4，表明城市出現(xiàn)概率與勞動(dòng)力分布密度表現(xiàn)出顯著的相關(guān)性。勞動(dòng)力比較密集的地區(qū)，城市出現(xiàn)的概率比較大，而稀疏地區(qū)，城市出現(xiàn)概率比較少，說(shuō)明勞動(dòng)力較密集的地區(qū)，其經(jīng)濟(jì)發(fā)展越快，越容易形成城市。從另一個(gè)側(cè)面說(shuō)明，城市的形成與市場(chǎng)經(jīng)濟(jì)的發(fā)達(dá)程度有關(guān)。

表3 勞動(dòng)力分布與城市出現(xiàn)概率關(guān)系

表4 勞動(dòng)力分布與城市出現(xiàn)概率相關(guān)系數(shù)

4 結(jié) 語(yǔ)

在地理非均質(zhì)空間下，城市數(shù)量與制造業(yè)勞動(dòng)力份額、需求價(jià)格彈性成正比，與可達(dá)性成反比，城市易于在地勢(shì)較平坦的地區(qū)形成。在經(jīng)濟(jì)非均質(zhì)空間下，城市數(shù)量依然與制造業(yè)勞動(dòng)力份額、需求價(jià)格彈性成正比，與可達(dá)性成反比，而城市易于在人口比較密集的地區(qū)形成，由于空間的非均質(zhì)性，城市分布的對(duì)稱性被破壞。結(jié)論還表明，城市的演繹不僅受到制造業(yè)勞動(dòng)力份額、需求價(jià)格彈性指數(shù)、可達(dá)性3個(gè)模型參數(shù)的影響，同時(shí)還會(huì)受到地形因子與經(jīng)濟(jì)因子的影響。

由于本文所選課題的先前研究成果并不多，另外受個(gè)人水平及客觀條件的限制，故研究工作仍需進(jìn)一步深入探討和完善，主要表現(xiàn)在以下兩方面：①本文模擬所用到的交通網(wǎng)絡(luò)采用的是虛擬生成的網(wǎng)絡(luò)，因此，采用現(xiàn)實(shí)存在的交通網(wǎng)絡(luò)進(jìn)行可達(dá)性度量有待于進(jìn)一步研究；②探索地理邊界對(duì)城市生長(zhǎng)最終結(jié)果的影響。

[1] Stelder D．The Simultaneous Prediction of Place and Size of Cities∶ New Results of a Geographical Agglomeration Model for Europe[C]．In HWWA Conference∶ New Economic Geography,Closing the Gap between Theory and Empirics,2004

[2] 皮永華,徐少君,潘文安．基于區(qū)位角度的城市化與經(jīng)濟(jì)發(fā)展之間互動(dòng)關(guān)系的實(shí)證研究——以浙江省為例[J]．經(jīng)濟(jì)地理,2005,25(4)∶502-505

[3] 李后強(qiáng)．關(guān)于城市演化的非線征收科學(xué)方法[J]．經(jīng)濟(jì)地理,2000,20(2)∶7-9

[4] 劉繼生．城鎮(zhèn)體系空間結(jié)構(gòu)的分形維數(shù)及其測(cè)算方法[J]．地理研究,1999,18(2)∶171-172

[5] 陳彥光．城市規(guī)模分布的三參數(shù)Zipf模型-Davis二倍規(guī)律的理論推廣及其分形性質(zhì)的實(shí)證分析[J]．華中師范大學(xué)學(xué)報(bào)∶自然科學(xué)版,2000,34(1)∶98-101

[6] 湯國(guó)安,楊昕．ArcGIS地理信息系統(tǒng)空間分析實(shí)驗(yàn)教程[M]．北京∶科學(xué)出版社,2006

[7] 劉繼生．城鎮(zhèn)體系空間結(jié)構(gòu)的分形維數(shù)及其測(cè)算方法[J]．地理研究,1999,18(2)∶171-172

P208

B

1672-4623（2016）01-0049-04

10.3969/j.issn.1672-4623.2016.01.015

杜艷琴，碩士，工程師，研究方向?yàn)橥恋乩眉俺鞘邪l(fā)展。

2014-03-24。