基于面板Logistic增長模型中國城鎮(zhèn)化演進特征與趨勢分析

劉洪濤，尚 進，蒲學吉

（西北工業(yè)大學管理學院，西安 710072）

一、引言

城鎮(zhèn)化伴隨工業(yè)化與現(xiàn)代化進程發(fā)展，是社會發(fā)展的長遠方向和必然趨勢。而當今世界上主要的工業(yè)大國已經(jīng)基本上完成了城鎮(zhèn)化進程，少數(shù)國家甚或已出現(xiàn)逆城鎮(zhèn)化的情形。自改革開放以來，中國城鎮(zhèn)化發(fā)展迅猛。1978年以來，中國城鎮(zhèn)化快速發(fā)展，但改革開放初期發(fā)展較為粗放，這一階段國家開始逐漸減少對農(nóng)村勞動力轉移至城鎮(zhèn)的限制；1996年起，中國城鎮(zhèn)化逐漸進入加速發(fā)展階段，至2004年以后，城鎮(zhèn)中基礎設施等投資增大，城鎮(zhèn)化發(fā)展水平亦顯著提升?？梢?，伴隨城鎮(zhèn)化發(fā)展，城鎮(zhèn)建設不斷擴張，投資力度亦持續(xù)增強。城鎮(zhèn)化作為支撐國家經(jīng)濟增長和社會發(fā)展的重要引擎受到廣泛重視，黨的十八大報告亦指出要加快推進中國城鎮(zhèn)化進程，深入開展以人為本的新型城鎮(zhèn)化。明確中國當前城鎮(zhèn)化狀況并據(jù)此預測城鎮(zhèn)化未來發(fā)展方向，將有利于認清當前發(fā)展形勢并深刻理解黨中央相關政策，進而促進中國城鎮(zhèn)化健康持續(xù)發(fā)展。

二、文獻綜述

Northam（1979）研究表明城鎮(zhèn)化的增長過程一般呈現(xiàn)一條“S”型曲線的增長規(guī)律，且包含三個階段：起始期（Initial stage），城鎮(zhèn)化率低且增長速度較為緩慢；加速期（Acceleration stage），城鎮(zhèn)化率增長速度較快；末尾期（Terminal stage），城鎮(zhèn)化率高但增長速度逐漸降低且城鎮(zhèn)化率增長趨于平緩，該研究為Logistic增長模型的發(fā)展和城鎮(zhèn)化率演進規(guī)律的預測奠定了基礎[1]。王遠飛和張超（1997）探討了用Logistic模型來描述城市化過程的合理性，預測了中國城市化率的發(fā)展趨勢[2]。喬松珊和孫成金（2016）依據(jù)城鄉(xiāng)人口增長率的變化特點建立優(yōu)化的Logistic模型并進行實證分析[3]。此外，亦有相關學者將改進的Logistic增長模型運用于省域或縣域數(shù)據(jù)以分析城鎮(zhèn)化發(fā)展狀況[4][5]。上述文獻對Logistic增長模型進行優(yōu)化，均站在整體角度對研究對象城鎮(zhèn)化演進狀況進行分析。

陳彥光和周一星（2005），馮年華（2005）修正了諾瑟姆曲線，得到一個四階段的劃分結果并進行地理空間解釋，進而對中國城鎮(zhèn)化進程進行預測[6][7]。王建軍（2009）運用高等數(shù)學方法，推出城鎮(zhèn)化Logistic模型關于時間變量t的方程表達式，并推導出曲線的三個特征點及其數(shù)學表達式[8]。陳素平（2015）等采用回歸分析方法對城鎮(zhèn)化發(fā)展各階段特征點進行探索[9]。張樂勤和張勇（2015）通過Logistic模型分析了中國城鎮(zhèn)化快速演進期的時長及增長速度[10]。上述文獻將Logistic增長模型劃分為若干階段，描述了各階段的城鎮(zhèn)化研究特征，有利于了解中國城鎮(zhèn)化發(fā)展進程各階段的發(fā)展動向以及當前發(fā)展態(tài)勢，為中國城鎮(zhèn)化下一階段發(fā)展方向、發(fā)展速度提供參考。

前人在研究過程中大都通過直接測算整體城鎮(zhèn)化水平來考量研究對象的城鎮(zhèn)化進程，或者簡單依照自然地理狀況，將中國各省市分成東部、中部、西部三個地區(qū)。直接采用整體城鎮(zhèn)化水平籠統(tǒng)地分析中國城鎮(zhèn)化的發(fā)展進程會將高杠桿值（省份）包含于其中，產(chǎn)生以偏概全的狀況；簡單分類后再分析所得到的結論亦顯粗放，這些均不利于對全國總體形勢的認識。本文試圖將全國30個省、直轄市、自治區(qū)分組，通過面板數(shù)據(jù)聚類分析，根據(jù)城鎮(zhèn)化率特征重新構造距離函數(shù)并以此衡量各省市的相似性，進而建立面板數(shù)據(jù)Logistic增長模型，分別對各組的城鎮(zhèn)化發(fā)展狀況進行分析，排除高杠桿值（省份），減小了其對全國城鎮(zhèn)化率可能產(chǎn)生的影響，兼顧了各地區(qū)城鎮(zhèn)化進程的特殊性與一般性，使對城鎮(zhèn)化率的測度方法更系統(tǒng)化、客觀化。并據(jù)此探索中國城鎮(zhèn)化發(fā)展的特征點，預測城鎮(zhèn)化發(fā)展方向，以達到全面了解中國城鎮(zhèn)化發(fā)展態(tài)勢的目的，為把握當前城鎮(zhèn)化發(fā)展形勢，明確中國城鎮(zhèn)化所處階段，制定未來城鎮(zhèn)化戰(zhàn)略性指導政策提供科學理論依據(jù)。

三、研究方法與數(shù)據(jù)來源

（一）面板數(shù)據(jù)的聚類分析

Bonzo（2002）定義聚類標準時以概率連接函數(shù)代替?zhèn)鹘y(tǒng)的距離函數(shù)，將聚類過程作為優(yōu)化問題處理[11]。Ren（2009）基于Fisher次序集群理論，通過Frobenius準則重建了Ward函數(shù)提出多變量面板數(shù)據(jù)序聚類方法[12]。李因果和何曉群（2010）基于上述研究，提出在對面板數(shù)據(jù)進行聚類分析時需要考慮：（1）各樣本指標間的靜態(tài)距離；（2）時間序列的動態(tài)發(fā)展特征距離[13]。其綜合考量了面板數(shù)據(jù)“絕對指標”，“增量指標”及“時序波動”等特征，重新構建面板數(shù)據(jù)相似性測度的距離函數(shù)，并通過Ward聚類算法計算類與類之間的距離，以此為基礎提出面板數(shù)據(jù)聚類分析方法。本文沿用其思想，仍采用歐式距離計算面板數(shù)據(jù)相似性指標的距離測度。鑒于本文采用了單指標面板數(shù)據(jù)，無法考慮“時序波動”因素，又考慮到城鎮(zhèn)化率的演變過程中應受到政策等推動因素的影響，故本文加入了加速度歐氏距離指標，以更好的描繪城鎮(zhèn)化發(fā)展過程中推動力的作用效果。在聚類分析中綜合考慮了城鎮(zhèn)化發(fā)展的歷史與現(xiàn)狀、經(jīng)濟和政治等因素，使可以更完善地描繪城鎮(zhèn)化率的演變過程，以更好地實現(xiàn)聚類，則根據(jù)指標特征重新定義距離測度公式如下。

對面板數(shù)據(jù)中各數(shù)值用雙下標變量表示，即xi,t,i=1,2,…,N,t=1,2,…,T，表示第i個個體中t時間的數(shù)據(jù)，其中N表示面板數(shù)據(jù)包含的個體數(shù)；T表示各橫截面時間序列最大值。

定義1 個體間“直接歐氏距離”（Direct Eu?clidean Distance），可記為di,j(DED)

di,j(DED)描繪了整個時間測量范圍內個體間數(shù)值上距離的遠近，在本文中表示對i省份與j省份城鎮(zhèn)化率的差異（歷史和發(fā)展現(xiàn)狀因素差異）測度。

定義2 個體間“速度歐氏距離”（Speed Eu?clidean Distance），可記為di,j(SED)

其中 ?xi,t=xi,t+1-xi,t，?xj,t=xj,t+1-xj,t。?xi,t，?xj,t表示相鄰時期個體的絕對值差異。di,j( )SED描繪了整個時間測量范圍內隨時間變化個體間速度的差異，在本文中表示對i省份與j省份城鎮(zhèn)化率增長速度的差異（經(jīng)濟因素差異）測度。

定義3 個體間“加速度歐氏距離”（Accelerat?ed Speed Euclidean Distance），可記為di,j(ASED)

定義4 個體間“綜合加權歐式距離”（Com?prehensive Weighted Euclidean Distance），可記為di,j(CWED)

其中，α+β+γ=1，α、β、γ代表相應距離的權重系數(shù)。綜合加權歐式距離中權重系數(shù)可以通過實際問題進行主觀判定或者客觀判定，常用方法如專家評分法、層次分析法、模糊判別法等。通過觀察上式可知，若β=γ=0，則綜合加權歐式距離等同于直接歐式距離，即傳統(tǒng)歐式距離測度。

本文采用系統(tǒng)聚類法中的Ward法，其基本思想為方差分析，若分類得當則同類樣本間的離差平方和最小，不同類樣本間離差平方和相對較大。

（二）面板數(shù)據(jù)的Logistic增長模型

韓本毅（2011）在研究中表明中國城鎮(zhèn)化進程從1978年之后逐步體現(xiàn)Logistic曲線的特征，此后眾多學者亦采取類似方法對不同范圍的城鎮(zhèn)化演進過程進行探索[14]。對此，本文亦采用Logistic模型考察中國城鎮(zhèn)化進程的發(fā)展規(guī)律，但又依據(jù)聚類分析結果將全國30個省、直轄市、自治區(qū)分成M類，選取每一類數(shù)據(jù)的重心值形成一個橫截面，建立M個橫截面T時期的面板Logistic模型。

傳統(tǒng)的Logistic模型表達式為：

其中：U代表城鎮(zhèn)化水平，UMax代表城鎮(zhèn)化水平的上限，a代表積分常數(shù)，b代表增長系數(shù)，t代表時間，t=1,2,…,T。為使便捷地確定回歸模型的參數(shù)，令a=lnUMax·c，b=-lnd，即：

其中：c代表模型回歸常數(shù)，d代表回歸系數(shù)。

本文基于上式，將其擴展為面板模型如下：

其中：Ui代表第i個橫截面的城鎮(zhèn)化水平，ci代表第i個橫截面的模型回歸常數(shù)，di代表第i個橫截面的回歸系數(shù)，UMaxi代表第i個橫截面的城鎮(zhèn)化水平上限。為方便進行擬合，可將上述模型線性化后再進行回歸。

（三）Logistic增長模型的階段劃分以及分界點測量

圖1 Logistic增長模型示意圖

Logistic增長模型的總體趨勢可直觀地表示如圖1。

圖2 Logistic增長模型一階導數(shù)示意圖

同樣地，亦可依據(jù)圖1作出Logistic增長模型導數(shù)及二階導數(shù)圖像如圖2、圖3。

圖3 Logistic增長模型二階導數(shù)示意圖

本文依據(jù)圖2、圖3可在圖1中找到三點(TM,UM)，(TU,UU)，(TL,UL)。并可依據(jù)三點將城鎮(zhèn)化發(fā)展的整體進程劃分為三個階段，找到城鎮(zhèn)化發(fā)展速度最快的時期。其中，(TM,UM)是使一階導數(shù)達到最大值的點，也即城鎮(zhèn)化增長速度最快的時間點以及其對應的城鎮(zhèn)化水平；(TU,UU)是使二階導數(shù)達到最大值的點，也即城鎮(zhèn)化增長加速度達到最大值的時間點以及其對應的城鎮(zhèn)化水平；(TL,UL)是使二階導數(shù)達到最小值的點，也即城鎮(zhèn)化增長加速度達到最小值的時間點以及其對應的城鎮(zhèn)化水平。

第一階段（(TU,UU)之前）：該階段城鎮(zhèn)化水平相對較低，增長速度較緩慢，加速度不斷增加直至最大值；第二階段（(TU,UU)，(TL,UL)之間）：該階段城鎮(zhèn)化水平快速提高，增長速度不斷增大至(TM,UM)達到最大值之后開始逐漸減小，同時該階段加速度不斷減小，(TM,UM)之前加速度為正值，(TM,UM)之后加速度為負值并在(TL,UL)達到最小值；第三階段（(TL,UL)之后）：該階段城鎮(zhèn)化率處于較高水平且仍繼續(xù)緩慢增長，增長速度逐漸減小并趨近于零，加速度由最小值逐漸提高且向零逼近。

可以求出(TM,UM)，(TU,UU)，(TL,UL)三點的具體值如下。

首先對Logistic回歸方程求二階導數(shù)，(TM,UM)點對應橫坐標值應為使二階導數(shù)等于零的自變量值。

Logistic回歸方程為：

一階導數(shù)表達式為：

二階導數(shù)表達式為：

則 當 1-UMax·c·elnd·t=0 時 ，也 即 當 TM=-時，此時城鎮(zhèn)化率增長速度達到最大值，此時的城鎮(zhèn)化率為：

又，(TU,UU)，(TL,UL)的橫坐標值應是使三階導數(shù)等于零的自變量值，對（10）式求導，

三階導數(shù)表達式為：

令T=UMax·c·elnd·t，則當 T2-4T+1=0時可求得，也即，

將TU，TL代入（12）式得到此時對應的城鎮(zhèn)化率為：

則：

（四）數(shù)據(jù)來源

鑒于韓本毅在研究中表明，中國城鎮(zhèn)化進程從1978年之后逐步體現(xiàn)Logistic曲線的特征[14]，本文采用1979～2014年全國各省、直轄市、自治區(qū)（將重慶市相關數(shù)據(jù)并入四川省以方便計算）年末城鎮(zhèn)人口與總人口比值計算各省市城鎮(zhèn)化率，所有數(shù)據(jù)均來自國家統(tǒng)計局發(fā)布的《新中國60年統(tǒng)計資料匯編》，缺失的數(shù)據(jù)由各省統(tǒng)計年鑒數(shù)據(jù)予以補足。

四、實證分析

（一）中國各省、直轄市、自治區(qū)城鎮(zhèn)化率的面板數(shù)據(jù)聚類分析

本文對中國30個省、直轄市、自治區(qū)的36年間城鎮(zhèn)化率相關數(shù)據(jù)采取前述方法進行聚類分析。通過專家評分法，本文將綜合加權歐式距離的測算中所需權重系數(shù)分別賦予α=0.71，β=0.24，γ=0.05，利用Java編寫程序進行聚類，并生成譜系聚類樹狀圖如圖4。

圖4 基于綜合加權歐式距離與Ward法的中國各省城鎮(zhèn)化率聚類分析譜系樹狀圖

本文最終將30個省、直轄市、自治區(qū)分為6類。表1給出了基于綜合加權歐式距離聚類的分類結果：第一類：北京、上海；第二類：西藏；第三類：內蒙古、遼寧、天津、黑龍江、吉林、新疆、青海、寧夏、山西；第四類：云南、貴州；第五類：河北、海南、福建、江蘇、湖北、山東、廣西、湖南、河南、四川、陜西、江西、甘肅、安徽；第六類：浙江、廣東。

表1 基于綜合加權歐式距離與Ward法的中國各省城鎮(zhèn)化率聚類結果

（二）中國城鎮(zhèn)化率的Logistic增長模型回歸分析、階段劃分及預測

1.中國城鎮(zhèn)化率的Logistic增長模型回歸分析及階段劃分

關于UMax的確定方法較多，國內外相關學者均對該值進行了相關推導與預測，國務院發(fā)展研究中心在《中國城鎮(zhèn)化——前景戰(zhàn)略與政策（2010）》中預測了中國城鎮(zhèn)化水平的峰值，約為80%。本文立足于中國國情，亦參考這一權威數(shù)據(jù)。在本文中，假定各類省份的城鎮(zhèn)化率上限均相等且等于全國城鎮(zhèn)化率上限，又考慮到第一類省份（北京、上海）的平均城鎮(zhèn)化率在2014年就已經(jīng)達到87%以上，故本文將城鎮(zhèn)化率上限設定在90%。

利用第三章第一節(jié)中聚類分析分類結果（表1），本文采用每類城鎮(zhèn)化率重心代表該類的城鎮(zhèn)化整體水平，以時間為自變量，以各類城鎮(zhèn)化率水平重心為因變量，以類別劃分橫截面建立面板Logistic增長模型，擬合圖像如圖5所示。

可得到回歸方程為：

圖5 1979～2014年各類省份面板數(shù)據(jù)Logistic增長模型回歸曲線擬合

表2 各類省份城鎮(zhèn)化特征及現(xiàn)階段城鎮(zhèn)化水平排序

回歸方程調整R2為0.978，除第一個面板的回歸常數(shù)外其余各參數(shù)均在0.01顯著水平下顯著，證明回歸模型擬合較好。分析圖5可得出如下信息：第一類（北京、上海）的城鎮(zhèn)化水平在全國各類省份中處于最高水平；第二類（西藏）、第四類（云南、貴州）的城鎮(zhèn)化水平在全國各類省份中處于最低水平；其余各類均處在中間水平。

進一步分析圖5可知，第一類（北京、上海）城鎮(zhèn)化水平從本文測度范圍的起始點處就遠高于其余各類省份，且該類省份在上世紀90年代之前一直保持較高的增長速度，從本世紀初開始，該類省份已逐漸達到發(fā)達國家的城鎮(zhèn)化水平，城鎮(zhèn)化率繼續(xù)增長，增長速度逐漸放緩；第三類（內蒙古、遼寧等）城鎮(zhèn)化水平在本文測量范圍的起始點處處于中高水平，且在整個測度范圍內城鎮(zhèn)化率增長速度保持相對穩(wěn)定；第五類（河北、海南等）、第六類（浙江、廣東）城鎮(zhèn)化水平從本文測量始點起處于較低水平，在整個測度范圍內兩類省份城鎮(zhèn)化率增長速度均不斷增加，而第五類省份城鎮(zhèn)化率增速相對第六類省份變化較緩，截至2014年，第五類省份逐漸逼近第三類省份的城鎮(zhèn)化水平，接近世界平均城鎮(zhèn)化水平，而第六類省份在測度時期內城鎮(zhèn)化率始終保持較高速度增長，至2014年已接近一般發(fā)達國家的城鎮(zhèn)化水平；第二類（西藏）、第四類（云南、貴州）城鎮(zhèn)化水平自測量始點起就處于全國最低水平，兩類省份在整個測量時期內均保持較緩慢的增長，且增長速度均有所增加。各類省份城鎮(zhèn)化特征描述及現(xiàn)階段城鎮(zhèn)化水平排序列于表2。

第一類、第二類以及第四類省份的城鎮(zhèn)化率值與其余各省份偏離較遠，在擬合回歸中如同高杠桿值，可能導致回歸曲線拉向高杠桿值方向從而造成回歸不準確?？紤]到第一類省份受國家政策導向，城鎮(zhèn)化建設均早于全國其他省份；第二類、第四類省份地處全國最西南，其地勢多為山區(qū)，城鎮(zhèn)化建設較難且起步較晚，故以上三類省份均不適宜代表中國城鎮(zhèn)化發(fā)展狀況。為了較準確地測度中國整體城鎮(zhèn)化水平，本文利用上述聚類分析結果，本文剔除城鎮(zhèn)化率遠高于大多數(shù)省份的第一類（北京、上海）省份和較低于其他省份的第二類（西藏）和第四類（云南、貴州）省份，選取其余各省份的城鎮(zhèn)化率重心作為中國整體城鎮(zhèn)化率水平，以此方法計算，兼顧了各省市城鎮(zhèn)化水平的獨特性和全國整體水平的一般性。根據(jù)以上分析，回歸后可得到中國城鎮(zhèn)化水平的Logistic回歸曲線如圖6。

得到回歸方程為：

回歸方程調整R2為0.969，模型擬合較好。

計算由本文所述方法測定的中國城鎮(zhèn)化水平與國家統(tǒng)計局依照傳統(tǒng)方法測定數(shù)據(jù)的偏差如表3所示。

可見，通過本文所述測量方法得到的城鎮(zhèn)化率整體趨勢和國家統(tǒng)計局數(shù)據(jù)較為接近，具有一定科學性，同時又排除了高杠桿值影響，對國家統(tǒng)計局數(shù)據(jù)進行了適度修正。

根據(jù)前面的結論又可計算城鎮(zhèn)化率增長曲線的三個分界點的坐標：

由計算結果可知，約在2007年前后，中國城鎮(zhèn)化率增長加速度為零，此時增長速度達到頂峰，全國城鎮(zhèn)化率為45%。中國城鎮(zhèn)化快速發(fā)展階段為1982年至2034年共約53年，這一階段始終保持較高的城鎮(zhèn)化率增長速度，城鎮(zhèn)化率加速度不斷減小，至2034年，城鎮(zhèn)化率加速度達到最小。1982年之前，中國城鎮(zhèn)化處于起始期，城鎮(zhèn)化率加速度不斷提高。2034年之后，中國城鎮(zhèn)化開始逐漸進入末尾期，這一階段中國城鎮(zhèn)化率處于較高水平，城鎮(zhèn)化率增速仍不斷減小。綜上所述，當前中國正處于城鎮(zhèn)化快速發(fā)展階段后期，這一階段城鎮(zhèn)化率仍保持相對較高的增長速度，但增長速度逐漸下降，城鎮(zhèn)化率的加速度已為負值且不斷減小。

圖6 1979～2014年中國城鎮(zhèn)化水平的Logistic回歸曲線擬合

2.中國城鎮(zhèn)化率趨勢預測

根據(jù)前文中得出的結論，對中國城鎮(zhèn)化未來近40年的發(fā)展態(tài)勢進行預測如圖7。

根據(jù)圖7，結合第三章第二節(jié)第一部分中相關分析可以看出，中國城鎮(zhèn)化水平在未來40年內仍會持續(xù)增長，增長速度不斷放緩，預計2034年達到中國城鎮(zhèn)化進程中的(TL,UL)點，即城鎮(zhèn)化率加速度最低點。2015年至2034年的時間段內，中國城鎮(zhèn)化率加速度始終為負值且不斷降低。2034年之后，中國城鎮(zhèn)化率加速度仍為負值但逐漸增大。按本文所建立的模型預測，直至2050年中國總體城鎮(zhèn)化率可達到80.48%，逐漸接近中等發(fā)達國家的城鎮(zhèn)化水平。

表3 本文方法測量數(shù)據(jù)對國家統(tǒng)計局數(shù)據(jù)的修正比例

圖7 中國城鎮(zhèn)化率預測的Logistic回歸曲線擬合

五、結論與啟示

（一）結論

借鑒了前人的研究方法與思路，考慮到按地理因素對省份進行分類的不足之處，本文對中國30個省份36年的城鎮(zhèn)化率數(shù)據(jù)進行聚類分析。依據(jù)綜合加權歐氏距離將各省份分為六類，即：第一類（京、滬），第二類（藏），第三類（蒙、遼、津、黑、吉、新、青、寧、晉），第四類（云、貴），第五類（冀、瓊、閩、蘇、鄂、魯、桂、湘、豫、川、山、贛、甘、皖）和第六類（浙、粵）。對各類省份的城鎮(zhèn)化率進行面板數(shù)據(jù)Logistic增長模型回歸后，分析了各省份的城鎮(zhèn)化發(fā)展特點，其中第一類省份由于“起點高且增速快”位于各類省份中第一梯隊；第三、五、六類省份由于“起點高或增速快”位于各類省份中第二梯隊；而第二、四類省份則由于“起點低且增速慢”位于各類省份中第三梯隊。由于第一、二、四類與其余各類省份城鎮(zhèn)化水平差距較大，可視為高杠桿值，可能對最終回歸結果產(chǎn)生影響。在排除高杠桿值問題的基礎上，依據(jù)模型求出中國城鎮(zhèn)化發(fā)展過程中的階段分界點，說明了中國正處于城鎮(zhèn)化發(fā)展的快速發(fā)展階段后期。最后，本文對未來40年中國城鎮(zhèn)化的發(fā)展進程進行預測。數(shù)據(jù)表明，在未來40年內，中國城鎮(zhèn)化水平仍會不斷提升，但提升速度逐漸放緩，至2050年中國城鎮(zhèn)化率將逐漸接近中等發(fā)達國家的城鎮(zhèn)化水平。

本文對面板聚類分析中距離測度方式進行改進，加入了加速度歐氏距離指標，在聚類分析中綜合考慮了城鎮(zhèn)化發(fā)展的歷史與現(xiàn)狀、經(jīng)濟和政治等因素，使可以更完善地描繪城鎮(zhèn)化率的演變過程，以更好地實現(xiàn)聚類。在此基礎上又對中國城鎮(zhèn)化率的測度方法進行了優(yōu)化，減小了高杠桿值對全國城鎮(zhèn)化率的影響，兼顧了各地區(qū)城鎮(zhèn)化進程的特殊性與一般性，使對城鎮(zhèn)化率的測度方法更系統(tǒng)化、客觀化。本文所述方法將有助于把握當前城鎮(zhèn)化發(fā)展形勢，明確中國城鎮(zhèn)化所處階段，為制定未來城鎮(zhèn)化戰(zhàn)略性指導政策提供科學理論依據(jù)。本文在對中國30個省份的城鎮(zhèn)化率數(shù)據(jù)進行面板聚類分析時僅考慮到各省份城鎮(zhèn)化率的關聯(lián)而沒有考慮到地域聯(lián)系，下一步研究可以在本文的分類方法基礎上進一步分析區(qū)域城鎮(zhèn)化發(fā)展，為各區(qū)域可持續(xù)發(fā)展和城鎮(zhèn)化提供更為詳細的政策制定參考。

（二）啟示

改革開放以來，中國城鎮(zhèn)化發(fā)展取得了可喜的成果，經(jīng)濟發(fā)展態(tài)勢不斷向好，總結改革開放前以及當今的經(jīng)驗教訓，符合中國經(jīng)濟發(fā)展需求的新型城鎮(zhèn)化發(fā)展道路是推動城鎮(zhèn)化健康發(fā)展的關鍵。根據(jù)本文相關結論，對中國城鎮(zhèn)化發(fā)展提出以下建議：

第一，著力開展集約型城鎮(zhèn)化建設。當前中國處于城鎮(zhèn)化快速發(fā)展階段后期，城鎮(zhèn)化發(fā)展速度較快，但增速逐年減緩。面對城鎮(zhèn)化增速下降的必然趨勢，粗放擴張式城鎮(zhèn)化路線將越來越不可取。改革開放以來，我國采用此方式開展城鎮(zhèn)化建設，并取得優(yōu)異的建設成果。然而，此方式帶給我們的不僅僅是城鎮(zhèn)化水平和國內生產(chǎn)總值的快速增長，亦透支了環(huán)境福利和土地福利等。隨著城鎮(zhèn)化建設進一步深入，在已經(jīng)到來的城鎮(zhèn)化發(fā)展新時期和城鎮(zhèn)化增速減緩新趨勢下，粗放式發(fā)展只能帶來城鎮(zhèn)化建設效益遞減以及遠高于產(chǎn)出的生態(tài)環(huán)境破壞和治理成本。為此，應著力開展集約型城鎮(zhèn)化建設，依據(jù)城市地理區(qū)位、資源稟賦、市場條件和交通運輸?shù)纫蛩?，合理進行城鎮(zhèn)化建設規(guī)劃與定位，以提高城市城鎮(zhèn)化建設投入產(chǎn)出比和城鎮(zhèn)化建設效率。同時，應正確處理好人口、資源與環(huán)境等諸多方面的制約關系，做到早發(fā)現(xiàn)早治理，開展可持續(xù)城鎮(zhèn)化建設。

第二，加快并擴大城市群建設。從依據(jù)中國各省市城鎮(zhèn)化發(fā)展狀況的分類結果中可以看出，除第二類（西藏）外，其余各類均包含不少于兩個省份。每類省份的城鎮(zhèn)化發(fā)展模式相似，其社會、經(jīng)濟、地理、人文等方面必然存在相似之處，同類型的省份共同開展城鎮(zhèn)化建設，相互依托相互扶持，可以使城鎮(zhèn)化建設由以省為單位擴展為以區(qū)域為單位，加大城市中心輻射范圍，增強城市規(guī)模效應，更有利促進城鎮(zhèn)化建設與區(qū)域經(jīng)濟發(fā)展。結合因地制宜的城鎮(zhèn)化發(fā)展政策，在每類省份中求同存異，兼顧了每個省份的特殊性與每一類省份的一般性，從而使全國城鎮(zhèn)化水平整體性提升。

第三，城鎮(zhèn)化建設重心應逐漸向內陸轉移。第一類（京、滬）和第二類（浙、粵）省份城鎮(zhèn)化水平在全國范圍內領跑。分析中國各省份城鎮(zhèn)化指標也可發(fā)現(xiàn)，沿海省份城鎮(zhèn)化水平較內陸地區(qū)普遍偏高。這在一定程度上是由于改革開放、區(qū)位交通因素等諸多政策及地理條件而引起的。隨著交通、物流、運輸業(yè)等的發(fā)展，區(qū)位因素的重要性逐漸下降；而伴隨著“振興東北老工業(yè)基地”、“西部大開發(fā)”等政策的頒布與實施，也逐漸平衡了內陸與沿海地區(qū)的政策支持。整體提升全國城鎮(zhèn)化水平僅依靠沿海地區(qū)是不夠的，更需要依靠廣闊的內陸地區(qū)。在城鎮(zhèn)化建設的新時期，應逐漸縮小沿海省份與內陸省份城鎮(zhèn)化差距，沿海省份帶動內陸省份，實現(xiàn)沿海、內陸城鎮(zhèn)化協(xié)調增長。

第四，加大對西南部省份的城鎮(zhèn)化發(fā)展建設力度。第二類（藏）、第四類（云、貴）省份整體城鎮(zhèn)化水平在全國各省份中處于最低的水平，此類省份大都處于西南邊境地區(qū)，經(jīng)濟發(fā)展多依靠旅游業(yè)帶動。自改革開放以來，該類省份城鎮(zhèn)化增長速度均較為緩慢，應制定相應政策加大對該地區(qū)的城鎮(zhèn)化建設力度，可通過產(chǎn)業(yè)結構調整，利用其資源優(yōu)勢，著力發(fā)展第三產(chǎn)業(yè)，充分發(fā)揮“結構紅利效應”，開發(fā)及合理利用當?shù)刈匀环A賦以調節(jié)生產(chǎn)要素的供需結構，推動旅游業(yè)建設與發(fā)展。在發(fā)展第三產(chǎn)業(yè)同時，積極把握產(chǎn)業(yè)間的關聯(lián)效應和擴散效應以帶動第一、二產(chǎn)業(yè)同步發(fā)展，誘導新技術、新原料等需求的出現(xiàn)以推動經(jīng)濟發(fā)展，以使其城鎮(zhèn)化率增長速度得到提升。從而縮小這一地區(qū)與其他區(qū)域城鎮(zhèn)化水平差異，帶動全國整體城鎮(zhèn)化水平提升。

[1]Northam，R.，M.Urban Geography（2nd ed）[M].New York：John Wiley&Sons，1979.

[2]王遠飛，張超.Logistic模型參數(shù)估計與我國城市化水平預測[J].經(jīng)濟地理，1997，17（4）：8-13.

[3]喬松珊，孫成金.優(yōu)化的Logistic模型及其在城鎮(zhèn)化水平預測中的應用[J].許昌學院院報，2016，35（5）：10-15.

[4]王帥，滕玉成.基于Logistic模型的縣域城鎮(zhèn)化率及其推進預測——以濟南市平陰縣為例[J].西安建筑科技大學學報（社會科學版），2014，33（6）：32-37.

[5]趙文英，葛禮霞.基于改進Logistic模型的黑龍江省城鎮(zhèn)化水平預測[J].數(shù)學的實踐與認識，2013，43（13）：44-49.

[6]陳彥光，周一星.城鎮(zhèn)化Logistic過程的階段劃分及空間解釋——對Northam曲線的修正與發(fā)展[J].經(jīng)濟地理，2005，25（6）：817-822.

[7]馮年華，宣衛(wèi)紅，戴軍，王書明.基于Logistic模型的我國城鎮(zhèn)化演進軌跡分析[J].生態(tài)經(jīng)濟，2015，31（7）：76-79.

[8]王建軍.城鎮(zhèn)化發(fā)展階段劃分[J].地理學報，2009，64（2）：177-188.

[9]陳素平，張樂勤，許信旺.基于Logistic模型的中國城鎮(zhèn)化演進階段特征及趨勢探析[J].干旱區(qū)地理，2015，38（2）：384-390.

[10]張樂勤，張勇.基于Logistic模型的中國城鎮(zhèn)化快速演進期時長、速度及啟示[J].河南大學學報（自然版），2015，45（5）：550-555.

[11]Bonzo D.C.，Hermosilla A.Y.Clustering Panel Data via Perturbed Adaptive Simulated Annealing and Genetic Algo?rithms[J].Advances in Complex Systems，2002（4）：339-360.

[12]Ren J.，Shi Sh L.Multivariable Panel Data Ordinal Clus?tering and Its Application in Competitive Strategy Identifica?tion of Appliance-wiring Listed Companies[C].International Conference on Management Science&Engineering（16 th），Moscow，Russia，2009：253-258.

[13]李因果，何曉群.面板數(shù)據(jù)聚類方法及應用[J].統(tǒng)計研究，2010，27（9）：73-79.

[14]韓本毅.國城市化發(fā)展進程及展望[J].西安交通大學學報（社會科學版），2011，31（3）：18-22.