基于主成分分析的西藏自治區(qū)城鎮(zhèn)化空間分布類型

余 波， 劉 松， 熊 鷹， 曹永旺， 周春山*

(1. 中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院，廣東省城市化與地理環(huán)境空間模擬重點(diǎn)實(shí)驗(yàn)室，廣州 510275；2. 國(guó)家開發(fā)銀行西藏分行，拉薩 850000)

基于主成分分析的西藏自治區(qū)城鎮(zhèn)化空間分布類型

余 波1,2， 劉 松1， 熊 鷹1， 曹永旺1， 周春山1*

(1. 中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院，廣東省城市化與地理環(huán)境空間模擬重點(diǎn)實(shí)驗(yàn)室，廣州 510275；2. 國(guó)家開發(fā)銀行西藏分行，拉薩 850000)

利用主成分分析法，選取表征西藏自治區(qū)城鎮(zhèn)化的人口、產(chǎn)業(yè)、用地、社會(huì)和自然環(huán)境等方面的變量，提取了西藏自治區(qū)城鎮(zhèn)化特征的主成分，并通過(guò)聚類分析劃分了西藏自治區(qū)城鎮(zhèn)化類型. 結(jié)論如下：(1)西藏自治區(qū)城鎮(zhèn)化內(nèi)容可歸納為7大主因子：“援藏因子”、“消費(fèi)因子”、“農(nóng)業(yè)因子”、“自然因子”、“城建因子”、“二三產(chǎn)業(yè)因子”和“教育因子”，體現(xiàn)了西藏自治區(qū)發(fā)展的特點(diǎn)；(2)2015年西藏自治區(qū)城鎮(zhèn)化可劃分為“中心地服務(wù)型”、“農(nóng)牧業(yè)生產(chǎn)型”、“二三產(chǎn)業(yè)發(fā)展型”、“援助依賴型”、“幫扶援建型”、“自然環(huán)境限制型”和“自然環(huán)境依托型”等7種類型，不同類型的城鎮(zhèn)化特征明顯；(3)西藏自治區(qū)城鎮(zhèn)化模式不同于全國(guó)大部分地區(qū)，城鎮(zhèn)化建設(shè)不能照搬其他地區(qū)經(jīng)驗(yàn).

城鎮(zhèn)化類型； 主導(dǎo)因子； 空間分布； 主成分分析； 西藏自治區(qū)

從模式角度研究城鎮(zhèn)化問題已經(jīng)成為城鎮(zhèn)化研究領(lǐng)域中的重要內(nèi)容[1]. 國(guó)內(nèi)城鎮(zhèn)化模式研究始于20世紀(jì)80年代后期，在90年代趨于豐富. 學(xué)者們從不同的視角研究城鎮(zhèn)化模式問題，如超前型和滯后型城鎮(zhèn)化模式[2-3]、集中型和分散型城鎮(zhèn)化模式[4-5]、大中小城鎮(zhèn)化模式、自上而下型和自下而上型城鎮(zhèn)化模式[6-7]、內(nèi)生型和外生型城鎮(zhèn)化模式[8]、城鄉(xiāng)分割型和城鄉(xiāng)統(tǒng)籌型城鎮(zhèn)化模式[9]，以及攀枝花、珠三角、蘇南、溫州等模式[10].

關(guān)于城鎮(zhèn)化模式概念，盛廣耀[11]認(rèn)同城鎮(zhèn)化模式是國(guó)家、地區(qū)城鎮(zhèn)化的實(shí)現(xiàn)途徑和方式，是對(duì)城鎮(zhèn)化的演進(jìn)過(guò)程、表現(xiàn)特征、動(dòng)力機(jī)制、戰(zhàn)略選擇的概括與總結(jié)；邊雪等[12]認(rèn)為已有的城鎮(zhèn)化模式研究是在總結(jié)成功經(jīng)驗(yàn)、識(shí)別典型問題、趨近高質(zhì)高效路徑，并利用人口城鎮(zhèn)化水平、產(chǎn)業(yè)非農(nóng)化水平和用地城鎮(zhèn)化水平等3個(gè)維度之間的協(xié)調(diào)關(guān)系來(lái)劃分城鎮(zhèn)化關(guān)系模式，為城鎮(zhèn)化模式研究打開了思路. 學(xué)者們也采用主成分和系統(tǒng)聚類方法，研究各地區(qū)城鎮(zhèn)化模式[13-15]，這些研究既豐富了城鎮(zhèn)化模式研究?jī)?nèi)容，也為研究西藏自治區(qū)城鎮(zhèn)化模式提供經(jīng)驗(yàn). 青海省經(jīng)濟(jì)研究院課題組對(duì)青藏高原地區(qū)城鎮(zhèn)化模式做出了一次有益嘗試，其研究認(rèn)為有別于內(nèi)地城市經(jīng)濟(jì)帶動(dòng)型模式，青藏高原地區(qū)城鎮(zhèn)化是依托宗教、民族文化、政權(quán)建設(shè)而建立，可以歸納為3種模式：宗教與政治經(jīng)濟(jì)中心融合型城鎮(zhèn)、政府主導(dǎo)的邊塞型城鎮(zhèn)、資源開發(fā)型城鎮(zhèn)[10].

由于城鎮(zhèn)化機(jī)制的復(fù)雜性，不同地區(qū)實(shí)現(xiàn)城鎮(zhèn)化的過(guò)程不同，從而形成相互差別的發(fā)展模式[11]. 2014年，西藏自治區(qū)第一、二、三產(chǎn)業(yè)比重為9.4∶36.7∶53.9，城鎮(zhèn)化率達(dá)到25.8%. 在這種情況下，西藏自治區(qū)城鎮(zhèn)化的地域分布特征值得研究. 已有西藏自治區(qū)城鎮(zhèn)化的研究較多地集中于城鎮(zhèn)化動(dòng)力機(jī)制[16-17]、特色產(chǎn)業(yè)發(fā)展[18]和產(chǎn)業(yè)結(jié)構(gòu)演進(jìn)[19]等方面，較少涉及模式. 鑒于此，本文根據(jù)城鎮(zhèn)化過(guò)程的特征，選取表征城鎮(zhèn)化的人口、產(chǎn)業(yè)、用地、社會(huì)和自然環(huán)境等方面的變量，分析西藏自治區(qū)城鎮(zhèn)化空間分布類型.

1 指標(biāo)體系與研究方法

1.1 度量西藏自治區(qū)城鎮(zhèn)化的指標(biāo)體系

基于城鎮(zhèn)化內(nèi)涵[20]和西藏自治區(qū)的獨(dú)特地理環(huán)境，考慮到指標(biāo)的代表性、差異性和易得性，本文從人口、產(chǎn)業(yè)、土地、社會(huì)城鎮(zhèn)化以及自然環(huán)境等5個(gè)方面衡量西藏自治區(qū)城鎮(zhèn)化發(fā)展水平，并構(gòu)建了指標(biāo)體系(表1),共包括13個(gè)一級(jí)指標(biāo)、31個(gè)二級(jí)指標(biāo). 研究數(shù)據(jù)主要來(lái)源于西藏自治區(qū)統(tǒng)計(jì)局和國(guó)家統(tǒng)計(jì)局西藏調(diào)查總隊(duì)對(duì)西藏自治區(qū)74個(gè)縣2015年經(jīng)濟(jì)社會(huì)發(fā)展基本情況的統(tǒng)計(jì)，少量缺失數(shù)據(jù)由西藏自治區(qū)5市2區(qū)統(tǒng)計(jì)年鑒補(bǔ)充，采用對(duì)數(shù)化處理以減少不同量綱對(duì)計(jì)算造成的不利影響.

表1 2015年西藏自治區(qū)城鎮(zhèn)化發(fā)展評(píng)價(jià)指標(biāo)體系

Table 1 The evaluation index of urbanization of Tibet Autonomous Region in 2015

指標(biāo)類別一級(jí)指標(biāo)二級(jí)指標(biāo)A：人口城鎮(zhèn)化A1：人口A11：戶籍人口密度/(人·km-2)A2：勞動(dòng)力A21：每萬(wàn)人城鎮(zhèn)在崗職工人數(shù)/人B：產(chǎn)業(yè)城鎮(zhèn)化B1：經(jīng)濟(jì)水平B11：人均GDP/元B12：工業(yè)總產(chǎn)值/萬(wàn)元B2：產(chǎn)業(yè)結(jié)構(gòu)B21：糧食總產(chǎn)量/噸B22：人均糧食產(chǎn)量/kgB23：油菜籽產(chǎn)量/噸B24：肉奶總產(chǎn)量/噸B25：農(nóng)林牧漁總產(chǎn)值/萬(wàn)元B26：第二產(chǎn)業(yè)產(chǎn)值比重B27：第三產(chǎn)業(yè)產(chǎn)值比重C：用地城鎮(zhèn)化C1：農(nóng)用地C11：耕地覆蓋率/‰C12：森林覆蓋率/‰C2：建設(shè)用地C21：城區(qū)人均用地/m2D：社會(huì)城鎮(zhèn)化D1：交通D11：公路密度/(km·km-2)D2：教育D21：萬(wàn)人普通小學(xué)在校學(xué)生數(shù)/人D22：萬(wàn)人普通中學(xué)在校學(xué)生數(shù)/人D3：醫(yī)療D31：萬(wàn)人擁有醫(yī)療床位數(shù)/張D4：社保D41：萬(wàn)人擁有社會(huì)福利收養(yǎng)性單位床位數(shù)/張D5：居民生活水平D51：農(nóng)村居民人均可支配收入/元D52：萬(wàn)人移動(dòng)電話年末用戶/戶D53：人均居民儲(chǔ)蓄存款余額/元D54：批發(fā)和零售業(yè)貿(mào)易總額/萬(wàn)元D55：住宿和餐飲業(yè)貿(mào)易總額/萬(wàn)元D6：政府行為D61：人均公共預(yù)算財(cái)政支出/元D62：全社會(huì)固定資產(chǎn)投資總額/萬(wàn)元D63：人均財(cái)政轉(zhuǎn)移支付量/元E：自然環(huán)境E1：自然條件E11：城鎮(zhèn)海拔高度/mE12：年均降水量/mLE13：年均日照時(shí)數(shù)/hE14：年平均氣溫/℃

1.2 分類模型選擇和主要步驟

由于表征城鎮(zhèn)化的變量多，且變量之間具有一定聯(lián)系，盡管能夠反映西藏自治區(qū)城鎮(zhèn)化的空間分布差異，但也增加了分類難度和復(fù)雜度. 為了從眾多變量中抓住西藏自治區(qū)城鎮(zhèn)化的整體特征，科學(xué)地劃分城鎮(zhèn)化類型，揭示西藏自治區(qū)城鎮(zhèn)化水平的區(qū)域差異，本文采用適合多變量分析的主成分分析法[21-22]，主要步驟如下：

首先，運(yùn)用SPSS軟件標(biāo)準(zhǔn)化31個(gè)二級(jí)指標(biāo)原始數(shù)值，形成31個(gè)彼此之間不相關(guān)的新的綜合變量，通過(guò)選擇其中影響成分大并具有代表意義的部分新變量(也就是下文所說(shuō)的因子)，對(duì)西藏自治區(qū)城鎮(zhèn)化發(fā)展進(jìn)行分析評(píng)價(jià).

其次，檢驗(yàn)主成分分析可行性. 通過(guò)檢驗(yàn)KMO系數(shù)和Bartlett系數(shù)進(jìn)行判斷. 經(jīng)SPSS軟件計(jì)算，本文構(gòu)建的西藏自治區(qū)城鎮(zhèn)化31個(gè)二級(jí)指標(biāo)的相關(guān)矩陣KMO檢驗(yàn)系數(shù)為0.731>0.700，根據(jù)常用的分析標(biāo)準(zhǔn)[23]，在可接受進(jìn)行主成分分析的范圍之內(nèi)；Bartlett檢驗(yàn)系數(shù)(Sig)<0.01，拒絕球形檢驗(yàn)的零假設(shè)，說(shuō)明適宜主成分分析(表2).

表2 2015年西藏自治區(qū)城鎮(zhèn)化主成分分析的KMO和Bartlett檢驗(yàn)

Table 2 KMO and Bartlett test of the principal components analysis of urbanization of Tibet Autonomous Region in 2015

采樣充足度的KMO檢驗(yàn)Bartlett球形度檢驗(yàn)近似卡方自由度(df)顯著性(Sig)0．7312472．2044650．000

注：數(shù)據(jù)由SPSS軟件計(jì)算得出.

第三，對(duì)本文構(gòu)建的西藏自治區(qū)城鎮(zhèn)化31個(gè)二級(jí)指標(biāo)的相關(guān)矩陣，運(yùn)用SPSS軟件，得到各因子的初始特征值，包括特征根、方差貢獻(xiàn)率和累計(jì)方差貢獻(xiàn)率(表3). 并根據(jù)確定主成分的分析標(biāo)準(zhǔn)(特征根>1)，確定前7個(gè)因子為分析西藏自治區(qū)城鎮(zhèn)化發(fā)展類型的7個(gè)主因子. 該7個(gè)主因子的累計(jì)旋轉(zhuǎn)方差貢獻(xiàn)率達(dá)到80.155%，概括了原有變量的大部分信息.

第四，對(duì)選取的7個(gè)主因子的初始特征值采用方差最大化正交旋轉(zhuǎn)，使共性因子載荷系數(shù)向更大或更小的方向變化，并保證原共性因子之間的正交性和共性方差總和不變. 經(jīng)過(guò)旋轉(zhuǎn)，得到正交旋轉(zhuǎn)之后的各成分特征值(表4)，同時(shí)得到了正交旋轉(zhuǎn)之后的2015年西藏自治區(qū)城鎮(zhèn)化主因子旋轉(zhuǎn)載荷矩陣(表5).

2 西藏自治區(qū)城鎮(zhèn)化發(fā)展類型分析

2.1 主因子命名與解釋

第1因子與“人均公共預(yù)算財(cái)政支出”、“人均財(cái)政轉(zhuǎn)移支付量”變量呈顯著正相關(guān)(相關(guān)系數(shù)分別為0.881、0.878)，與“農(nóng)林牧漁業(yè)總產(chǎn)值”、“肉奶總產(chǎn)量”變量呈顯著負(fù)相關(guān)(相關(guān)系數(shù)分別為-0.862、-0.822)，該因子與外部資金援助量密切相關(guān)，可命名為“援藏因子”. 第1主因子的表現(xiàn)是：當(dāng)因子值越大時(shí)，西藏自治區(qū)的外部資金援助量越大，農(nóng)業(yè)產(chǎn)值越小，反之亦然. 該因子得分較高的地區(qū)主要分布在西藏邊疆地區(qū)(圖1A).

表3 2015年西藏自治區(qū)城鎮(zhèn)化因子初始特征值

Table 3 The initial eigenvalues of the factors of urbanization of Tibet Autonomous Region in 2015

因子特征根l方差貢獻(xiàn)率/%累計(jì)貢獻(xiàn)率/%18．68128．00328．00325．42017．48545．48833．50611．31056．79842．8349．14365．94151．7175．53771．47961．5945．14176．62071．0963．53580．15580．8172．63582．79090．7922．55685．346100．5541．78987．135110．5111．65088．784120．4851．56390．347130．4011．29491．641140．3611．16692．807150．3131．01193．818160．2710．87594．694170．2500．80695．500180．2270．73396．232190．2100．67896．910200．1830．58997．499210．1530．49497．993220．1340．43298．425230．1230．39698．821240．0980．31799．138250．0710．23099．368260．0700．22799．595270．0410．13499．729280．0390．12699．855290．0350．11499．969300．0070．02299．991310．0030．009100．000

表4 2015年西藏自治區(qū)城鎮(zhèn)化主因子正交旋轉(zhuǎn)后特征值

Table 4 The eigenvalues of the principal factors of urbanization after orthogonal rotation of Tibet Autonomous Region in 2015

主因子特征根l方差貢獻(xiàn)率/%累計(jì)貢獻(xiàn)率/%15．47517．66217．66225．10216．45834．12034．42614．27748．39743．11110．03558．43152．4777．99166．42362．4467．89074．31371．8115．84280．155

表5 2015年西藏自治區(qū)城鎮(zhèn)化主因子旋轉(zhuǎn)載荷矩陣Table 5 The principal factor rotated component matrix of urbanization of Tibet Autonomous Region in 2015

第2因子與“人均GDP”、“農(nóng)村居民人均可支配收入”、“住宿和餐飲業(yè)貿(mào)易總額”、“批發(fā)和零售業(yè)貿(mào)易總額”、“全社會(huì)固定資產(chǎn)投資總額”、“萬(wàn)人移動(dòng)電話年末用戶”、“人均居民儲(chǔ)蓄存款余額”呈正相關(guān)(相關(guān)系數(shù)分別為0.909、0.767、0.705、0.691、0.677、0.656、0.633)，其中與“人均GDP”的相關(guān)性最強(qiáng)，第2主因子與消費(fèi)和經(jīng)濟(jì)內(nèi)容關(guān)系密切，可命名為“消費(fèi)因子”. 該因子的表現(xiàn)是：當(dāng)因子值越大時(shí)，人均GDP等7個(gè)變量越高，居民消費(fèi)能力越強(qiáng)，地區(qū)發(fā)展水平越好，反之亦然. 該因子得分較高的地區(qū)主要分布在西藏中部地區(qū)，包括卡若區(qū)、巴宜區(qū)、城關(guān)區(qū)、乃東區(qū)、桑珠孜區(qū)、那曲縣和噶爾縣(圖1B).

第3因子與“人均糧食產(chǎn)量”、“糧食總產(chǎn)量”、“油菜籽產(chǎn)量”、“年平均氣溫”、“耕地覆蓋率”呈顯著的正相關(guān)(相關(guān)系數(shù)分別為0.866、0.802、0.762、0.721、0.670)，其中與“人均糧食產(chǎn)量”的相關(guān)性最強(qiáng)，第3主因子與農(nóng)業(yè)發(fā)展內(nèi)容相關(guān)，可命名為“農(nóng)業(yè)因子”. 該因子的表現(xiàn)是：當(dāng)因子值越大，地區(qū)耕地面積越大，糧食和油菜產(chǎn)量越高，年均氣溫更高，反之亦然. 該因子得分較高的地區(qū)主要分布在“一江兩河”地區(qū)，范圍從日喀則市到山南市(圖1C).

圖1 2015年西藏自治區(qū)城鎮(zhèn)化因子得分空間分布特征

第4因子與“年均日照時(shí)數(shù)”、“城鎮(zhèn)海拔高度”、“森林覆蓋率”、“年均降水量”4個(gè)變量具有較強(qiáng)的相關(guān)性，其中與“年均日照時(shí)數(shù)”、“城鎮(zhèn)海拔高度”呈強(qiáng)烈的負(fù)相關(guān)(相關(guān)系數(shù)分別為-0.836、-0.800)，與“森林覆蓋率”、“年均降水量”呈強(qiáng)烈的正相關(guān)(相關(guān)系數(shù)分別為0.777、0.701)，第4主因子與自然環(huán)境內(nèi)容關(guān)系密切，可命名為“自然因子”. 該因子的表現(xiàn)是：當(dāng)因子值越大，地區(qū)海拔越低、日照時(shí)數(shù)越短，而森林覆蓋率越密、降水越豐富，反之亦然. 該因子得分較高的地區(qū)主要分布在藏東、藏北地區(qū)，包括昌都市、林芝市和那曲地區(qū)，因子得分低的地區(qū)主要集中在藏西北地區(qū)(圖1D).

第5因子與“公路密度”、“中心城區(qū)人均建設(shè)用地面積”、“人口密度”相關(guān)性較強(qiáng)，其中與“公路密度”、“人口密度”呈正相關(guān)(相關(guān)系數(shù)分別為0.812、0.625)，與“中心城區(qū)人均建設(shè)用地”呈負(fù)相關(guān)(相關(guān)系數(shù)為-0.691)，第5因子與城市建設(shè)內(nèi)容關(guān)系密切，可命名為“城建因子”. 該因子的表現(xiàn)是：當(dāng)因子值越大，地區(qū)公路里程越長(zhǎng)，交通越發(fā)達(dá)，人口密度越大，城區(qū)人均用地越少，反之亦然. 該因子得分較高的地區(qū)較分散，城關(guān)區(qū)、卡若區(qū)、乃東區(qū)和噶爾縣因子得分較高(圖1E).

第6因子與“第三產(chǎn)業(yè)產(chǎn)值比重”、“第二產(chǎn)業(yè)產(chǎn)值比重”2個(gè)因子呈顯著相關(guān)，其中與“第二產(chǎn)業(yè)產(chǎn)值比重”呈正相關(guān)(相關(guān)系數(shù)為0.772)、與“第三產(chǎn)業(yè)產(chǎn)值比重”呈強(qiáng)烈的負(fù)相關(guān)(相關(guān)系數(shù)為-0.901)，第6因子反映第二三產(chǎn)業(yè)產(chǎn)值占比狀況，可命名為“二三產(chǎn)業(yè)因子”. 該因子的表現(xiàn)是：當(dāng)因子值越大，地區(qū)第二產(chǎn)業(yè)產(chǎn)值比重越大，第三產(chǎn)業(yè)產(chǎn)值比重越低，產(chǎn)業(yè)結(jié)構(gòu)以二產(chǎn)為主，反之產(chǎn)業(yè)結(jié)構(gòu)以三產(chǎn)為主. 該因子得分較高的地區(qū)集中在拉薩市、山南市周邊縣，如堆龍德慶區(qū)、曲水縣、墨竹工卡縣和桑日縣因子得分較高(圖1F).

第7因子與“萬(wàn)人普通中學(xué)在校學(xué)生數(shù)”、“萬(wàn)人普通小學(xué)在校學(xué)生數(shù)”呈顯著的正相關(guān)關(guān)系(相關(guān)系數(shù)分別為0.874、0.813)，第7因子反映了地區(qū)的基礎(chǔ)教育狀況，可命名為“教育因子”. 該因子的表現(xiàn)是：當(dāng)因子值越大，接受基礎(chǔ)教育的人數(shù)越多，反之亦然. 該因子得分較高的地區(qū)空間分散，措美縣、聶拉木縣、吉隆縣、薩嘎縣、措勤縣、仲巴縣和改則縣等地區(qū)因子得分較高(圖1G).

2.2 城鎮(zhèn)化發(fā)展類型劃分

根據(jù)主因子得分，采用組間聯(lián)結(jié)、組內(nèi)聯(lián)結(jié)、最近鄰元素和Ward方法對(duì)主因子作聚類分析并比較聚類結(jié)果，發(fā)現(xiàn)Ward方法的結(jié)果最佳. 根據(jù)聚類結(jié)果，2015年西藏自治區(qū)城鎮(zhèn)化劃分為7種類型，其特征見表6、空間分布見圖2.

表6 2015年西藏自治區(qū)城鎮(zhèn)化發(fā)展類型特征Table 6 Characteristics of urbanization development pattern of Tibet Autonomous Region in 2015

圖2 2015年西藏自治區(qū)城鎮(zhèn)化發(fā)展類型空間分布

Figure 2 Spatial distribution of urbanization development pattern of Tibet Autonomus Region in 2015

第Ⅰ類型城鎮(zhèn)化的因子特點(diǎn)是：消費(fèi)、農(nóng)業(yè)、城建因子為正值，援藏、自然、二三產(chǎn)業(yè)、教育因子為負(fù)值，且消費(fèi)因子正值最大、城建因子次之，自然因子負(fù)值最大、援藏因子次之、二三產(chǎn)業(yè)和教育因子負(fù)值最小. 說(shuō)明該類型城鎮(zhèn)化的特點(diǎn)是經(jīng)濟(jì)發(fā)達(dá)、人口密集、產(chǎn)業(yè)以第三產(chǎn)業(yè)為主，并且此類型作為地區(qū)首府(拉薩市城關(guān)區(qū)、山南市乃東區(qū)和日喀則市桑珠孜區(qū))承擔(dān)為地區(qū)服務(wù)作用，為“中心地服務(wù)型”城鎮(zhèn)化. 這些地區(qū)的GDP總量、人均GDP和人口總數(shù)在西藏自治區(qū)74個(gè)縣區(qū)里排在前列，也是西藏自治區(qū)傳統(tǒng)的政治、經(jīng)濟(jì)和文化中心.

第Ⅱ類型城鎮(zhèn)化的因子特點(diǎn)是：農(nóng)業(yè)、教育因子的平均值為正值，城建、二三產(chǎn)業(yè)、消費(fèi)、援藏、自然因子的平均值為負(fù)值. 農(nóng)業(yè)因子的平均值與平方和均值最大，自然、援藏、消費(fèi)、城建、二三產(chǎn)業(yè)因子的平方和均值較大而平均值為負(fù)，說(shuō)明該類型城鎮(zhèn)化的特點(diǎn)是以農(nóng)業(yè)生產(chǎn)活動(dòng)為主，人均GDP、地區(qū)發(fā)展水平較低，為“農(nóng)牧業(yè)生產(chǎn)型”城鎮(zhèn)化. 這一類型城鎮(zhèn)化包括：拉薩市林周縣、當(dāng)雄縣、尼木縣，日喀則市南木林縣、白朗縣、江孜縣、康馬縣、薩迦縣、拉孜縣、定日縣、昂仁縣、仲巴縣、謝通門縣，山南市扎囊縣、貢嘎縣和浪卡子縣，即“一江兩河地區(qū)”，是西藏自治區(qū)傳統(tǒng)的農(nóng)牧業(yè)發(fā)達(dá)地區(qū).

第Ⅲ類型城鎮(zhèn)化的因子特點(diǎn)是：援藏、消費(fèi)、農(nóng)業(yè)、城建、二三產(chǎn)業(yè)因子的平均值為正值，自然、教育因子的平均值為負(fù)值. 二三產(chǎn)業(yè)因子平均值與平方和均值最大，消費(fèi)、農(nóng)業(yè)、援藏、城建因子平方和均值較大且平均值為正，自然、教育因子平方和均值較大但平均值為負(fù). 說(shuō)明該類型城鎮(zhèn)化以二三產(chǎn)業(yè)發(fā)展為主要特征，尤其是第二產(chǎn)業(yè)的發(fā)展明顯優(yōu)于第三產(chǎn)業(yè)，同時(shí)人均GDP、地區(qū)發(fā)展水平較高，為“二三產(chǎn)業(yè)發(fā)展型”城鎮(zhèn)化. 這一類型城鎮(zhèn)化包括:拉薩市曲水縣、堆龍德慶區(qū)、達(dá)孜縣、墨竹工卡縣，山南市桑日縣、加查縣和隆子縣. 曲水縣、堆龍德慶區(qū)和達(dá)孜縣作為拉薩近郊縣區(qū)，承接了從拉薩外溢出來(lái)的加工制造業(yè)和物流倉(cāng)儲(chǔ)服務(wù)；墨竹工卡縣擁有國(guó)內(nèi)儲(chǔ)量和加工能力排名前列的大型銅礦；桑日縣和加查縣是當(dāng)前西藏自治區(qū)水電資源重點(diǎn)開發(fā)地區(qū)；隆子縣鉛鋅礦資源豐富，已形成規(guī)?；a(chǎn). 這些區(qū)域也是自治區(qū)重點(diǎn)開發(fā)區(qū)域.

第Ⅳ類型城鎮(zhèn)化的因子特點(diǎn)是：援藏因子的平均值與平方和均值都大且為正值，農(nóng)業(yè)、城建因子的平方和均值最大但平均值為負(fù)值，消費(fèi)、自然、二三產(chǎn)業(yè)、教育因子平均值為負(fù)值. 說(shuō)明該類型城鎮(zhèn)化對(duì)外界援助的依賴很強(qiáng)，自身經(jīng)濟(jì)和城鎮(zhèn)化發(fā)展有限，為“援助依賴型”城鎮(zhèn)化. 這一類型城鎮(zhèn)化包括:阿里地區(qū)普蘭縣、札達(dá)縣、日土縣、革吉縣、改則縣、措勤縣，那曲地區(qū)安多縣、申扎縣、尼瑪縣、雙湖縣和山南市錯(cuò)那縣. 這些區(qū)域是西藏自治區(qū)地勢(shì)最高、人口最稀少的地區(qū)，生態(tài)環(huán)境相對(duì)脆弱，且大部分區(qū)域?yàn)橄拗崎_發(fā)區(qū)域，城鎮(zhèn)的發(fā)展對(duì)外界援助的依賴程度很大.

第Ⅴ類型城鎮(zhèn)化的因子特點(diǎn)是：援藏、農(nóng)業(yè)、城建、教育因子的平均值為正值，消費(fèi)、自然、二三產(chǎn)業(yè)因子的平均值為負(fù)值. 援藏因子的平均值與平方和均值最大且為正值，雖然地區(qū)發(fā)展水平較低，但農(nóng)業(yè)生產(chǎn)和二三產(chǎn)業(yè)產(chǎn)值較大于第Ⅳ類型(“援助依賴型”)，說(shuō)明該類型城鎮(zhèn)化以外部援助為主要特征，同時(shí)自身具備一定的發(fā)展條件，不完全依賴于外界援助，為“幫扶援建型”城鎮(zhèn)化. 這一類型城鎮(zhèn)化包括:阿里地區(qū)噶爾縣;林芝市朗縣，日喀則市仁布縣、亞東縣、崗巴縣、定結(jié)縣、聶拉木縣、薩嘎縣、吉隆縣,山南市瓊結(jié)縣、曲松縣、措美縣和洛扎縣. 這些區(qū)域大部分沿邊境分布，具有傳統(tǒng)的經(jīng)濟(jì)發(fā)展優(yōu)勢(shì)，部分區(qū)域，如瓊結(jié)縣、曲松縣和措美縣的農(nóng)牧業(yè)發(fā)展條件相對(duì)較好.

第Ⅵ類型城鎮(zhèn)化的因子特點(diǎn)是：自然、城建、二三產(chǎn)業(yè)因子的平均值為正值，援藏、消費(fèi)、農(nóng)業(yè)、教育因子的平均值為負(fù)值. 自然因子的平均值和平方和均值最大且為正值，說(shuō)明該類型城鎮(zhèn)化地區(qū)自然因子影響較強(qiáng)，雖然二三產(chǎn)業(yè)有一定的發(fā)展(產(chǎn)業(yè)結(jié)構(gòu)以二產(chǎn)為主)，但產(chǎn)值低于第Ⅶ類型、地區(qū)發(fā)展水平低、農(nóng)業(yè)生產(chǎn)弱，為“自然環(huán)境限制型”城鎮(zhèn)化. 這一類型城鎮(zhèn)化包括:昌都市卡若區(qū)、江達(dá)縣、貢覺縣、類烏齊縣、丁青縣、察雅縣、八宿縣、左貢縣、芒康縣、洛隆縣、邊壩縣，那曲地區(qū)班戈縣、那曲縣、索縣、比如縣、巴青縣、嘉黎縣和聶榮縣. 這些區(qū)域涵蓋了昌都市的所有縣區(qū)和那曲地區(qū)東部的大部分縣區(qū). 昌都市大部分縣區(qū)位于西藏自治區(qū)三江流域(金沙江、瀾滄江和怒江)，城鎮(zhèn)發(fā)展被河流限制；那曲地區(qū)東部縣區(qū)以高山谷地為主，對(duì)城鎮(zhèn)發(fā)展產(chǎn)生較多不利影響.

第Ⅶ類型城鎮(zhèn)化的因子特點(diǎn)是：援藏、消費(fèi)、農(nóng)業(yè)、自然、教育因子的平均值為正值，城建、二三產(chǎn)業(yè)因子的平均值為負(fù)值. 自然因子的平均值與平方和均值最大且為正，說(shuō)明該類型城鎮(zhèn)化地區(qū)自然因子影響最強(qiáng)烈，同時(shí)二三產(chǎn)業(yè)產(chǎn)值較高于第Ⅵ類型，產(chǎn)業(yè)結(jié)構(gòu)以三產(chǎn)為主，外界援助資金較多，人均GDP、地區(qū)發(fā)展水平高，為“自然環(huán)境依托型”城鎮(zhèn)化. 這一類型城鎮(zhèn)化包括:林芝市巴宜區(qū)、工布江達(dá)縣、米林縣、墨脫縣、波密縣和察隅縣. 這些區(qū)域海拔相對(duì)較低，受沿雅魯藏布江上溯的印度洋暖濕氣流影響，森林植被覆蓋率較高，自然環(huán)境相對(duì)較好，近年來(lái)在外部援助的支持下開發(fā)力度不斷加大.

3 結(jié)論與討論

本文根據(jù)西藏自治區(qū)城鎮(zhèn)化評(píng)價(jià)指標(biāo)劃分城鎮(zhèn)化發(fā)展類型，主要結(jié)論如下：

(1)西藏自治區(qū)城鎮(zhèn)化可以總結(jié)為7大主因子，分別是“援藏因子”、“消費(fèi)因子”、“農(nóng)業(yè)因子”、“自然因子”、“城建因子”、“二三產(chǎn)業(yè)因子”和“自然因子”. 城鎮(zhèn)化發(fā)展類型可以劃分為：“中心地服務(wù)型”、“農(nóng)牧業(yè)生產(chǎn)型”、“二三產(chǎn)業(yè)發(fā)展型”、“援助依賴型”、“幫扶援建型”、“自然環(huán)境限制型”和“自然環(huán)境依托型”.

(2)“中心地服務(wù)型”城鎮(zhèn)化地區(qū)主要是拉薩市、日喀則市、山南市；“農(nóng)牧業(yè)生產(chǎn)型”城鎮(zhèn)化地區(qū)主要分布在西藏“一江兩河”地區(qū)，該類型城鎮(zhèn)化的特點(diǎn)是以農(nóng)業(yè)生產(chǎn)活動(dòng)為主，人均GDP、地區(qū)發(fā)展水平較低；“二三產(chǎn)業(yè)發(fā)展型”城鎮(zhèn)化地區(qū)主要分布在拉薩市、山南市外圍地區(qū)，該類型城鎮(zhèn)化以二三產(chǎn)業(yè)發(fā)展為主要特征，尤其是第二產(chǎn)業(yè)的發(fā)展明顯優(yōu)于第三產(chǎn)業(yè)；“援助依賴型”城鎮(zhèn)化地區(qū)主要分布在阿里、那曲地區(qū)，該類型城鎮(zhèn)自身發(fā)展動(dòng)力不足，對(duì)外界援助的依賴很強(qiáng)；“幫扶援建型”城鎮(zhèn)化地區(qū)主要分布在西藏邊境地區(qū)，該類型城鎮(zhèn)化以外部援助為主要特征，同時(shí)自身具備一定的發(fā)展條件，發(fā)展并不完全依賴于援藏資金；“自然環(huán)境限制型”城鎮(zhèn)化地區(qū)主要分布在昌都市、阿里地區(qū)G317國(guó)道沿線，該類型城鎮(zhèn)化發(fā)展受自然環(huán)境限制較大，地區(qū)發(fā)展水平低；“自然環(huán)境依托型”城鎮(zhèn)化地區(qū)主要分布在藏東南林芝市，該類型城鎮(zhèn)化地區(qū)自然環(huán)境較優(yōu)越，人口密度小、城區(qū)人均建設(shè)用地多，農(nóng)業(yè)和二三產(chǎn)業(yè)具有一定規(guī)模.

本文是利用計(jì)量方法研究西藏自治區(qū)城鎮(zhèn)化空間分布類型的一次實(shí)踐. 首先，采用主成分分析方法和綜合指標(biāo)體系在研究西藏自治區(qū)城鎮(zhèn)化空間分布類型中是有效的，如“中心地服務(wù)型”城鎮(zhèn)化地區(qū)為拉薩市、山南市和日喀則市，在地區(qū)首府周邊形成了二三產(chǎn)業(yè)明顯發(fā)展的高地，在一江兩河流域是農(nóng)牧業(yè)占主導(dǎo)的城鎮(zhèn)化類型，這樣的分類結(jié)果與西藏自治區(qū)實(shí)情相符. 其次，主因子的提取結(jié)果彰顯了區(qū)域獨(dú)特性，如援藏因子、自然因子的提取，體現(xiàn)了西藏自治區(qū)城鎮(zhèn)化發(fā)展的特點(diǎn). 再次，西藏自治區(qū)城鎮(zhèn)化發(fā)展模式與全國(guó)大部分地區(qū)不同，甚至西藏自治區(qū)內(nèi)部不同地區(qū)城鎮(zhèn)化發(fā)展模式也不相同，因此，不能拿西藏自治區(qū)外已有的城鎮(zhèn)化發(fā)展模式硬套在西藏自治區(qū)城鎮(zhèn)化建設(shè)上來(lái)，而是要因地制宜地采取措施.

[1] 盛廣耀. 城市化模式研究綜述[J]. 城市發(fā)展研究,2011,18(7):13-19.

SHENG G Y. A literature review on the research of urbani-zation pattern[J]. Urban Studies,2011,18(7):13-19.

[2] 姜斌,李雪銘. 世界城市化模式及其對(duì)中國(guó)的啟示[J]. 世界地理研究,2007(1):40-45.

JIANG B,LI X M. World urbanization modes and their enlightenment to China[J]. World Regional Studies,2007(1):40-45.

[3] 陳明星,陸大道,查良松. 中國(guó)城市化與經(jīng)濟(jì)發(fā)展水平關(guān)系的國(guó)際比較[J]. 地理研究,2009(2):464-474.

CHEN M X,LU D D,ZHA L S. Urbanization and economic development in China:an international comparison based on quadrant map approach[J]. Geographical Research,2009(2):464-474.

[4] 毛蔣興,薛德升. 世界城市化模式及其對(duì)珠江三角洲的啟示[J]. 規(guī)劃師,2006(5):76-79.

MAO J X,XUE D S. The world’s urbanization pattern and its inspiration on the Pearl River Delta[J]. Planners,2006(5):76-79.

[5] 石憶邵. 德國(guó)均衡城鎮(zhèn)化模式與中國(guó)小城鎮(zhèn)發(fā)展的體制瓶頸[J]. 經(jīng)濟(jì)地理,2015(11):54-60;70.

SHI Y S. The balanced urbanization mode of Germany and the system bottlenecks of development of Chinese small towns[J]. Economic Geography,2015(11):54-60;70.

[6] 寧越敏. 新城市化進(jìn)程——90年代中國(guó)城市化動(dòng)力機(jī)制和特點(diǎn)探討[J]. 地理學(xué)報(bào),1998(5):470-477.

NING Y M. New process of urbanization——dynamics and features of urbanization in China since 1990[J]. Acta Geographica Sinica,1998(5):470-477.

[7] 崔功豪,馬潤(rùn)潮. 中國(guó)自下而上城市化的發(fā)展及其機(jī)制[J]. 地理學(xué)報(bào),1999(2):12-21.

CUI G H,MA LAURENCE J C. Urbanization from below in China:its development and mechanisms[J]. Acta Geographica Sinica,1999(2):12-21.

[8] 叢茂昆,張明斗. 內(nèi)生型城鎮(zhèn)化:新型城鎮(zhèn)化的模式選擇[J]. 南京農(nóng)業(yè)大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版),2016(3):30-36;157.

CONG M K,ZHANG M D. Endogenous urbanization:new urbanization’s pattern choice[J]. Journal of Nanjing Agricultural University(Social Sciences Edition),2016(3):30-36;157.

[9] 趙宇剛. 城鎮(zhèn)化視域下城鄉(xiāng)統(tǒng)籌模式的分析和選擇[J]. 江漢學(xué)術(shù),2013(2):14-21.

[10] 青海省經(jīng)濟(jì)研究院課題組. 青藏高原城市化發(fā)展模式研究[J]. 經(jīng)濟(jì)研究參考,2013(25):35-51.

[11] 盛廣耀. 關(guān)于城市化模式的理論分析[J]. 江淮論壇,2012(1):24-30.

[12] 邊雪,陳昊宇,曹廣忠. 基于人口、產(chǎn)業(yè)和用地結(jié)構(gòu)關(guān)系的城鎮(zhèn)化模式類型及演進(jìn)特征——以長(zhǎng)三角地區(qū)為例[J]. 地理研究,2013(12):2281-2291.

BIAN X,CHEN H Y,CAO G Z. Patterns of regional urbanization and its implications:an empirical study of the Yangtze River Delta region[J]. Geographical Research,2013(12):2281-2291.

[13] 安祥生. 山西省城市化類型的劃分[J]. 經(jīng)濟(jì)地理,2000(4):50-54.

AN X S. Study on classification of urbanization types of Shanxi Province[J]. Economic Geography,2000(4):50-54.

[14] 謝先全,晏路明. 福建城市化類型的數(shù)值劃分[J]. 地域研究與開發(fā),2004(1):21-24.

XIE X Q,YAN L M. The numerical division of urbanization types in Fujian Province[J]. Areal Research and Development,2004(1):21-24.

[15] 寧越敏,李健. 泛長(zhǎng)三角地區(qū)城鎮(zhèn)化的機(jī)制、模式與戰(zhàn)略[J]. 南京社會(huì)科學(xué),2009(5):8-14.

[16] 傅小鋒. 青藏高原城鎮(zhèn)化及其動(dòng)力機(jī)制分析[J]. 自然資源學(xué)報(bào),2000,15(4):369-374.

FU X F. An analysis of urbanization and its dynamical mechanism in Qinghai-Tibet Plateau[J]. Journal of Natu-ral Resources,2000,15(4):369-374.

[17] 唐偉,鐘祥浩,周偉. 西藏高原城鎮(zhèn)化動(dòng)力機(jī)制的演變與優(yōu)化——以“一江兩河”地區(qū)為例[J]. 山地學(xué)報(bào),2011(3):378-384.

TANG W,ZHONG X H,ZHOU W. The evolvement and optimization of urbanization dynamics in “three-river area” on Tibetan Plateau[J]. Journal of Mountain Science,2011(3):378-384.

[18] 孫勇,周煒,倪邦貴,等. 中國(guó)西藏農(nóng)牧業(yè)特色產(chǎn)業(yè)報(bào)告. 2008[M]. 北京:中國(guó)藏學(xué)出版社,2010.

[19] 劉剛,沈鐳. 1951～2004年西藏產(chǎn)業(yè)結(jié)構(gòu)的演進(jìn)特征與機(jī)理[J]. 地理學(xué)報(bào),2007,62(4):364-376.

LIU G,SHEN L. Characteristics and mechanism of Tibet’s industrial structure evolution from 1951 to 2004[J]. Acta Geographica Sinica,2007,62(4):364-376.

[20] 許學(xué)強(qiáng),周一星,寧越敏. 城市地理學(xué)[M]. 2版. 北京:高等教育出版社,2009:55.

[21] 張利,雷軍,張小雷,等. 烏魯木齊城市社會(huì)區(qū)分析[J]. 地理學(xué)報(bào),2012(6):817-828.

ZHANG L,LEI J,ZHANG X L,et al. Analysis of the urban social areas in Urumqi[J]. Acta Geographica Sinica,2012(6):817-828.

[22] 盧紋岱,朱紅兵. SPSS統(tǒng)計(jì)分析[M]. 北京:電子工業(yè)出版社,2015:499-512.

[23] 數(shù)據(jù)分析. KMO和Bartlett球形度檢驗(yàn)[EB/OL]. (2016-05-08)[2017-04-10]. http://www.cdadata.com/8145.

Study on Spatial Distribution Pattern of Urbanization of Tibet Autonomous Region Based on Principal Component Analysis

YU Bo1,2, LIU Song1, XIONG Ying1, CAO Yongwang1, ZHOU Chunshan1*

(1. Guangdong Key Laboratory for Urbanization and Geo-Simulation, School of Geography and Planning, Sun Yat-Sen University, Guangzhou 510275, China; 2. China Development Bank Tibet Branch， Lasa 850000, China)

Using the principal component analysis method, the variables characterising the urbanization of Tibet Autonomous Region selected from population, industry, land, social aspects, and natural environment are standar-dized to extract the principal components representing the characteristics of the urbanization of the Tibet Autonomous Region. The urbanization of the Tibet Autonomous Region is divided into different patterns by clustering analysis. Conclusions are drawn as follows: (1) Factors of Tibet Autonomous Region’s urbanization are summarized as 7 main factors, namely “outside aiding factor”, “consumption factor”, “agricultural factor”, “natural factor”, “construction factor”,“industrial factor” and “education factor”, and the main factor that embodies the characteristics of Tibet Autonomous Region. (2) The urbanization pattern of Tibet Autonomous Region in 2015 can be divided into “service centered type”, “agricultural environment produced type”, “industrial development type”, “aiding dependent type”, “aiding type”, “natural environment restricted type” and “natural environment supporting type”. The characteristics of different types of urbanization are obvious. (3) The pattern of urbanization of Tibet Autonomous Region is different from that in most areas of China. The experience of urbanization construction of other areas cannot be copied to Tibet Autonomous Region.

urbanization pattern; main factor; spatial distribution; principal component analysis; Tibet Autonomous Region

2017-04-20 《華南師范大學(xué)學(xué)報(bào)(自然科學(xué)版)》網(wǎng)址：http://journal.scnu.edu.cn/n

國(guó)家自然科學(xué)基金項(xiàng)目(40971097,41271182)

*通訊作者:周春山，教授，Email:zhoucs@mail.sysu.edu.cn.

K901

A

1000-5463(2017)06-0092-09

【中文責(zé)編：莊曉瓊 英文審校:肖菁】