侯賀平 劉艷芳 李紀偉 孔雪松

(1.武漢大學資源與環(huán)境科學學院，湖北武漢430079;2.武漢大學教育部地理信息系統(tǒng)重點實驗室，湖北武漢430079)

基于改進輻射模型的鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡研究

侯賀平1，2劉艷芳1，2李紀偉1，2孔雪松1，2

(1.武漢大學資源與環(huán)境科學學院，湖北武漢430079;2.武漢大學教育部地理信息系統(tǒng)重點實驗室，湖北武漢430079)

鄉(xiāng)鎮(zhèn)尺度的人口流動聯(lián)系和空間結(jié)構(gòu)研究對統(tǒng)籌區(qū)域城鎮(zhèn)體系規(guī)劃和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整具有重要意義。本文以復雜網(wǎng)絡理論為研究視角，將空間可達性引入傳統(tǒng)輻射模型進行參數(shù)修正，構(gòu)建鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡;從節(jié)點和社團結(jié)構(gòu)以及無標度和小世界特征等方面，探討在地域差異影響下人口流動的空間分布格局和復雜網(wǎng)絡特征，并以荊門市55個鄉(xiāng)鎮(zhèn)為例進行實證研究。結(jié)果表明:①荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的節(jié)點度分布以高發(fā)展水平鄉(xiāng)鎮(zhèn)為核心，呈現(xiàn)“一橫二縱”的空間分異格局，與荊門市城鎮(zhèn)體系和產(chǎn)業(yè)布局存在極強的耦合關(guān)系;②全網(wǎng)包含四個社團，各社團內(nèi)部聯(lián)系緊密且空間結(jié)構(gòu)特征各異，社團間聯(lián)系相對松散。③受空間距離和地理環(huán)境影響，網(wǎng)絡的無標度特征不夠明顯;④網(wǎng)絡具有一定的小世界特征，整體聯(lián)通性和局部集聚性較好。本文通過對傳統(tǒng)模型進行方法改進，實現(xiàn)了異質(zhì)空間下鄉(xiāng)鎮(zhèn)人口流動聯(lián)系及分布的可視化表達，實例研究結(jié)果可為荊門市城鎮(zhèn)規(guī)劃和產(chǎn)業(yè)布局提供決策思路。

鄉(xiāng)鎮(zhèn)人口流動;復雜網(wǎng)絡;輻射模型;可達性;荊門市

人口流動是市場經(jīng)濟資源優(yōu)化配置的內(nèi)在要求，是人口現(xiàn)代化、城鎮(zhèn)一體化和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整的必然結(jié)果。人口流動的規(guī)模和空間格局直接影響著區(qū)域的社會經(jīng)濟結(jié)構(gòu)和資源的組合分配，關(guān)系到區(qū)域的城鎮(zhèn)體系規(guī)劃和產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級。人口流動研究一直是學術(shù)界廣泛關(guān)注的熱點，相關(guān)研究主要集中在驅(qū)動因素［1］、效應分析［2］、分布特征［3］和政策措施［4］等方面。已有的人口流動研究在內(nèi)容上多把出發(fā)地和目的地孤立分析，缺少對不同地區(qū)間的人口流動聯(lián)系及其空間特征研究;研究方法多采用人口統(tǒng)計指標、回歸分析、空間自相關(guān)等，利用復雜數(shù)學模型探討不同地區(qū)之間人口流動聯(lián)系的研究并不多見。

新興的復雜網(wǎng)絡理論為地區(qū)間的人口流動研究提供了一個全新的視角。復雜網(wǎng)絡以復雜系統(tǒng)實體以及實體間的相互作用或關(guān)聯(lián)構(gòu)建網(wǎng)絡，節(jié)點之間不同的聯(lián)結(jié)方式和聯(lián)結(jié)程度將形成具有不同結(jié)構(gòu)和功能的復雜網(wǎng)絡，利用統(tǒng)計物理學分析網(wǎng)絡結(jié)構(gòu)及其動力學特征。復雜網(wǎng)絡模型主要包括小世界網(wǎng)絡和無標度網(wǎng)絡。1998年Strogatz和wattsk提出了融合一般網(wǎng)絡優(yōu)點的小世界網(wǎng)絡模型，其既具有規(guī)則網(wǎng)絡的高聚集性，又有類似隨機網(wǎng)絡的較小平均路徑距離，能夠更好地反映真實網(wǎng)絡系統(tǒng)［5］;1999年Barabási和Albert發(fā)現(xiàn)world wide web網(wǎng)絡的點和度分布符合冪律特征，將此稱之為無標度網(wǎng)絡［6］，其與小世界網(wǎng)絡共同表現(xiàn)出新增節(jié)點連接的優(yōu)先選擇特征。為深入探討復雜網(wǎng)絡的結(jié)構(gòu)及其演化規(guī)律，眾多學者對基礎模型進行了不同程度的改進［7－8］，并廣泛應用于神經(jīng)組織［9］、知識網(wǎng)絡［10］、人際關(guān)系［11］、交通網(wǎng)絡［12］等領(lǐng)域。

已有復雜網(wǎng)絡模型著重于網(wǎng)絡的拓撲結(jié)構(gòu)，在體現(xiàn)人口流動網(wǎng)絡的空間地域特征方面具有一定局限性。人口流動網(wǎng)絡的空間特征表現(xiàn)在人口流動的空間關(guān)聯(lián)上，受到節(jié)點的地理位置和其他因素的綜合影響。區(qū)域人口流動網(wǎng)絡屬于整體網(wǎng)絡，人口流動聯(lián)系的數(shù)據(jù)收集一般不能采用抽樣調(diào)查法［13］，在缺乏完整“O － D”(Origin-Destination)人口流動統(tǒng)計數(shù)據(jù)的情況下，空間數(shù)理模型模擬方法得到了最為廣泛的應用。將物理定律引入到人口流動研究中的重力模型是具有代表性的模型之一［14］;介入機會模型假定兩地間的人口流動強度正比于目的地提供的機會，反比于目的地距離半徑內(nèi)的其他機會［15］;隨機效用模型通過比較不同地區(qū)產(chǎn)生的綜合效用大小決定目的地選擇［16］。上述模型的一些假設條件缺少量化基礎、參數(shù)設置需要大量的經(jīng)驗數(shù)據(jù)進行標定，計算過程也較為復雜，有待進一步改進。輻射模型能夠較好地解決上述問題［17］，它通過出發(fā)地、目的地、影響范圍(以出發(fā)地為圓心、兩地距離為半徑的圓，出發(fā)地和目的地除外)的人口分布確定兩地之間的人口流動強度，需要指出的是輻射模型影響范圍的確定是基于均質(zhì)的理想空間，忽視了不同空間地域格局對人口流動的影響，存在一定的應用局限性。

基于此，本文首先嘗試以復雜網(wǎng)絡理論為指導，借助基于空間可達性的輻射模型測度鄉(xiāng)鎮(zhèn)間的人口流動強度，將鄉(xiāng)鎮(zhèn)抽象為節(jié)點、人口流動強度設為邊構(gòu)建鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡，其次，按照節(jié)點、局部網(wǎng)絡、整體網(wǎng)絡的順序?qū)W(wǎng)絡的節(jié)點結(jié)構(gòu)、社團結(jié)構(gòu)、無標度和小世界等整體特征進行分析，深入探討在不同的地域格局影響下，地區(qū)間的人口流動聯(lián)系及其空間結(jié)構(gòu)特征。以荊門市55個鄉(xiāng)鎮(zhèn)為例進行實證研究，為城鎮(zhèn)體系規(guī)劃和產(chǎn)業(yè)結(jié)構(gòu)布局提供科學決策。

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

1.1 基于改進輻射模型的網(wǎng)絡構(gòu)建方法

1.1.1 輻射模型及其改進

輻射模型(radiation model)是2012年美國學者Simini借鑒固體物理學中物質(zhì)運動的發(fā)散和吸收過程提出的，彌補了重力模型應用于人口流動強度模擬存在的多項不足，并利用人口流、物流、信息流等進行了實例驗證。與重力模型相比，輻射模型與實測數(shù)據(jù)的擬合程度較高，能夠更趨近現(xiàn)實地模擬兩地之間的流動強度。輻射模型將人口流動看做一個受聯(lián)合概率支配的隨機過程，取決于出發(fā)地、目的地和影響范圍的人口分布。經(jīng)過嚴格的理論邏輯推導，得到輻射模型的方程表達式(式1、2)，如下所示。模型在應用過程中不需要進行復雜的參數(shù)設置，在缺少歷史統(tǒng)計資料的情況能夠?qū)Φ貐^(qū)間的人口流動強度進行有效估算。

式中:Tij為i，j兩地之間的預期人口流動強度，mi、nj分別為出發(fā)地i和目的地j的總?cè)丝?，Sij為i到j之間人口流動的影響范圍總?cè)丝?以出發(fā)地為圓心，兩地距離為半徑的圓，ij兩地除外)，Ti為i鎮(zhèn)的流動人口數(shù)，V為i鎮(zhèn)的流動人口比重。

在Simini等學者的研究中，影響范圍的確定基礎是均質(zhì)的理想空間(圖1a)，而實際上，人口流動行為以真實的地理空間為載體，受到交通條件、地形地貌、土地利用狀況等多種因素及其分布特征的綜合影響(圖1b)。當出行的時間成本一定時，流動人口在出發(fā)地的各個方向上的出行距離各不相同(圖1c)，輻射模型用均質(zhì)的圓形作為影響范圍忽視了真實地理空間的異質(zhì)性，有待進一步改進。本文嘗試基于空間地域差異，從空間可達性的角度利用出發(fā)地的等時圈確定人口流動的異質(zhì)性影響范圍。根據(jù)輻射模型中對影響范圍的定義，若i鎮(zhèn)到j鎮(zhèn)的時間成本為tij，則i鎮(zhèn)的tij等時圈為i，j兩鎮(zhèn)之間人口流動的影響范圍，稱為異質(zhì)影響范圍。范圍內(nèi)的總?cè)丝赟’ij(式3)為等時圈內(nèi)各個鄉(xiāng)鎮(zhèn)的人口總和(i，j兩鎮(zhèn)的人口除外)，i鎮(zhèn)到異質(zhì)性影響范圍內(nèi)的任一鄉(xiāng)鎮(zhèn)的時間成本小于鎮(zhèn)i到鎮(zhèn)j的時間成本。

式中:S’ij為i鎮(zhèn)到j鎮(zhèn)的異質(zhì)影響范圍總?cè)丝冢琌k為異質(zhì)影響范圍內(nèi)k鄉(xiāng)鎮(zhèn)的人口，tik為i鎮(zhèn)到其他鄉(xiāng)鎮(zhèn)k的時間成本，tij為i鎮(zhèn)到j鎮(zhèn)的時間成本。

1.1.2 鄉(xiāng)鎮(zhèn)空間可達性的測度

目前對于空間可達性的測度方法主要有四種思路:①道路網(wǎng)分析［18］:通過構(gòu)建道路網(wǎng)矢量數(shù)據(jù)庫中點狀、線狀等要素的拓撲關(guān)系，形成O－D矩陣確定每個節(jié)點的可達性，不足是對路網(wǎng)密度要求高，且沒有考慮路網(wǎng)的非節(jié)點;②基于柵格的道路網(wǎng)分析［19］:在道路網(wǎng)分析的基礎之上把矢量數(shù)據(jù)柵格化，通過空間插值把一維空間(線)可達性計算推廣到二維空間(面)，雖然考慮了非節(jié)點但是其精度有待提高;③柵格分析［20］:借助ArcGIS軟件的空間分析功能計算最短成本加權(quán)距離，可以比較好的模擬現(xiàn)實，工作量較小，但是無法描述封閉道路(高速公路、鐵路)的出入口問題;④集成法［21］:道路網(wǎng)分析和柵格分析相結(jié)合，提高計算精度的同時解決了封閉道路的特殊性，但是計算步驟復雜。考慮到鄉(xiāng)鎮(zhèn)區(qū)域的路網(wǎng)密度較低，高速公路、鐵路的出入口一般在縣級及以上地區(qū)且沿線均有國道或者省道通過等因素，這里不予考慮，故在此采用柵格分析方法測度鄉(xiāng)鎮(zhèn)的可達性。

圖1 改進輻射模型示意圖Fig.1 Sketch map of the improved method for radiation model

由于人們在不同土地利用類型上的通勤方式和出行速度不同，可達性測度的實現(xiàn)需要設定不同土地利用類型的時間成本。根據(jù)《土地利用現(xiàn)狀分類(GB/T21010－2007)》將土地利用類型分為交通運輸用地(鐵路、機場用地、港口碼頭用地、管道運輸用地除外)、水域及水利設施用地和其他用地(指交通設施用地、水域及水利設施用地之外連續(xù)的其他土地利用類型部分)，交通運輸用地根據(jù)道路等級分為國道、省道、縣道和農(nóng)村道路。參考《中華人民共和國公路工程技術(shù)標準(JTGB01－2003)》，結(jié)合路網(wǎng)密度和路網(wǎng)質(zhì)量，實際運行速度分別設為國道75 km/h、省道50 km/h、縣道30 km/h、農(nóng)村道路15 km/h;水域及水利設施用地對可達性影響較為顯著，設為不能通行的隔離層;根據(jù)其他用地對通行的限制程度將其分為Ⅰ類、Ⅱ類、Ⅲ類①注:Ⅰ類:商服用地、工礦倉儲用地、住宅用地、公共管理與公共服務用地、特殊用地、城鎮(zhèn)村及工礦用地、交通運輸用地中的機場用地、港口碼頭用地;Ⅱ類:耕地、園地、其他土地中的空閑地、設施農(nóng)用地、田坎;Ⅲ類:林地、草地、其他農(nóng)用地中的鹽堿地、沼澤地、沙地、裸地。，并賦予相應的運行速度:Ⅰ類 15 km/h、Ⅱ類4 km/h、Ⅲ類3 km/h。綜上所述，得到不同空間對象每100 m的時間成本:國道 0.08 min、省道 0.12 min、縣道0.2 min、農(nóng)村道路0.4 min、Ⅰ類 0.4 min、Ⅱ類1.5 min、Ⅲ類2 min。

可達性測算的具體實現(xiàn)步驟為:借助ArcGIS10.1軟件，在土地利用現(xiàn)狀數(shù)據(jù)中提取所有參與柵格分析的圖層，進行重分類。對各圖層賦予相應的時間成本后做柵格化處理，本文確定的柵格大小為100 m×100 m?？臻g疊加各圖層取最小值得到各空間地物的時間成本柵格。以高程、坡向數(shù)據(jù)作為修正因素，在 ArcGIS中運行 Path Distance命令生成各鄉(xiāng)鎮(zhèn)可達性空間分布圖。在ArcGIS二次開發(fā)程序的支持下，實現(xiàn)對兩兩鄉(xiāng)鎮(zhèn)之間的可達性值提取，構(gòu)成可達性時間成本矩陣，以此確定兩兩鄉(xiāng)鎮(zhèn)對應的等時圈和異質(zhì)影響范圍，根據(jù)各個異質(zhì)性影響范圍與其他鄉(xiāng)鎮(zhèn)的空間關(guān)系確定其范圍內(nèi)的人口總數(shù)。

1.1.3 鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的構(gòu)建程序

根據(jù)復雜網(wǎng)絡分析的思路，將各個鄉(xiāng)鎮(zhèn)抽象為節(jié)點、鄉(xiāng)鎮(zhèn)之間的人口流動強度設為聯(lián)接邊構(gòu)建鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡。具體構(gòu)建步驟:①網(wǎng)絡節(jié)點定位:將節(jié)點設置在各鄉(xiāng)鎮(zhèn)鎮(zhèn)政府所在地;②鄉(xiāng)鎮(zhèn)人口流動強度矩陣:借助改進的輻射模型和各鄉(xiāng)鎮(zhèn)人口統(tǒng)計數(shù)據(jù)，得到鄉(xiāng)鎮(zhèn)人口流動強度矩陣;③矩陣標準化:為消除不同人口規(guī)模帶來的影響，將鄉(xiāng)鎮(zhèn)人口流動強度矩陣進行標準化處理;④二值矩陣轉(zhuǎn)化:選擇正確的切分點是二值矩陣轉(zhuǎn)換的關(guān)鍵，合適的切分值不至于使網(wǎng)絡聯(lián)接邊數(shù)目過多或過少，對網(wǎng)絡結(jié)構(gòu)的穩(wěn)定性影響較小，有利于正確描述網(wǎng)絡結(jié)構(gòu)特征。轉(zhuǎn)換規(guī)則是:當i鎮(zhèn)和j鎮(zhèn)之間的人口流動強度高于切分值，則賦值矩陣單元(i，j)的數(shù)值轉(zhuǎn)換成1，反之則為0;⑤對稱化處理:雖然鄉(xiāng)鎮(zhèn)間的人口流動過程具有方向性，即兩地間的流動人口從出發(fā)地流向目的地。本文旨在討論鄉(xiāng)鎮(zhèn)之間的人口流動聯(lián)系及其空間網(wǎng)絡結(jié)構(gòu)，兩地之間一旦發(fā)生人口流動，無論方向如何，即建立了人口流動聯(lián)系，對兩地同時產(chǎn)生影響。因此在文本中忽略兩地間人口流動方向的差異性，采取對稱化處理的方式，將二值矩陣轉(zhuǎn)化為對稱矩陣。

1.2 網(wǎng)絡復雜性分析方法

網(wǎng)絡復雜性的分析指標主要包括無標度(scalefree)、小世界(small－world)等復雜網(wǎng)絡特征的判斷，中心性(Centrality)、結(jié)構(gòu)洞(Structural holes)、網(wǎng)絡的層次(Hierarchy)結(jié)構(gòu)、社團結(jié)構(gòu)(Community)等網(wǎng)絡靜態(tài)結(jié)構(gòu)分析，以及可靠性(Reliability)或魯棒性(Robustness)、傳播性(Epidemic)、同步性(Synchronization)等網(wǎng)絡動力學分析。本文重點考察人口流動網(wǎng)絡節(jié)點結(jié)構(gòu)和社團結(jié)構(gòu)的空間分異規(guī)律以及無標度、小世界等復雜網(wǎng)絡特征，其中節(jié)點結(jié)構(gòu)空間分異規(guī)律由節(jié)點度標度，社團結(jié)構(gòu)特征由社團化指數(shù)標度，無標度性由度分布形態(tài)標度，小世界性由平均路徑距離和簇系數(shù)標度。

1.2.1 無標度特征

度和度分布:度表示與節(jié)點i連接的邊數(shù)目，記為ki，節(jié)點度的分布情況用概率分布函數(shù)p(k)來描述，當網(wǎng)絡規(guī)模較小時常用P(k)表示p(k)的度分布的累計概率分布(Cumulative Degree Distribution)函數(shù)(式4)。度是描述網(wǎng)絡節(jié)點結(jié)構(gòu)特性的基本參數(shù)，度分布則反映網(wǎng)絡系統(tǒng)的宏觀統(tǒng)計特征，當P(k)或者p(k)為冪函數(shù)時，節(jié)點度分布不均衡，少數(shù)節(jié)點的度較高，而大部分節(jié)點則較低，網(wǎng)絡具有無標度特征。近年的實證研究表明多數(shù)現(xiàn)實網(wǎng)絡屬于冪律特征為P(k)∝k－r(2≤r≤3)的無標度網(wǎng)絡［22］。

1.2.2 小世界特征

平均路徑距離:距離dij為網(wǎng)絡中連接i，j兩點之間的最短路徑的邊數(shù)，網(wǎng)絡的平均路徑距離L為任意兩節(jié)點之間距離的平均值(式5)。平均路徑距離可衡量整個網(wǎng)絡

簇系數(shù):一個節(jié)點i的簇系數(shù)Ci為ki個鄰節(jié)點之間聯(lián)接邊的數(shù)目占最大可能連邊數(shù)目的比值(式6)，Ci∈［0，1］，反映網(wǎng)絡節(jié)點集聚情況的局部屬性指標，節(jié)點的簇系數(shù)越大，表示相鄰節(jié)點內(nèi)部聯(lián)系越緊密。整體網(wǎng)絡的簇系數(shù)定義為所有節(jié)點簇系數(shù)的平均值(式7)，C值越大，表示整個網(wǎng)絡中各節(jié)點之間形成短距離聯(lián)系的程度越大。的傳輸性能與效率，L值越小，表示網(wǎng)絡中任意節(jié)點之間的拓撲距離越小，網(wǎng)絡的整體可達性越好。

平均路徑距離和簇系數(shù)是衡量小世界特征的測度指標，同時具有短路徑和高聚類特征的網(wǎng)絡屬于小世界網(wǎng)絡，在節(jié)點數(shù)目大而各點的度有限(一般小于節(jié)點數(shù)目)的情況下，任意兩個節(jié)點之間仍可通過相對較短的路徑相聯(lián)系，網(wǎng)絡的傳輸效率較高，交通網(wǎng)絡、社會關(guān)系等現(xiàn)實網(wǎng)絡的小世界特征較為明顯。

社團是指網(wǎng)絡中的節(jié)點之間相對較強、直接、緊密、經(jīng)常的或者積極的關(guān)系集合，屬于網(wǎng)絡的子結(jié)構(gòu)。社團結(jié)構(gòu)分析的目的主要是為了從宏觀上揭示網(wǎng)絡的子結(jié)構(gòu)內(nèi)部以及子結(jié)構(gòu)之間的關(guān)系，實現(xiàn)網(wǎng)絡整體結(jié)構(gòu)的簡要表達。社團的合理劃分是社團結(jié)構(gòu)分析的基礎，根據(jù)劃分規(guī)則的不同，可以分為凝聚算法、分裂算法、搜索算法以及其他。劃分效果的好壞利用社團化指數(shù)Q(式8)判斷，Q的取值范圍為 0.3 －0.7 為佳［23］。

其中，Q為模塊性指數(shù)，NM為模塊數(shù)量，L為城鎮(zhèn)網(wǎng)絡中有效聯(lián)接邊數(shù)，ls為社團s內(nèi)部城鎮(zhèn)間有效聯(lián)接邊數(shù)，ds為社團s內(nèi)部所有城鎮(zhèn)的度值總和。

1.3 研究區(qū)域與數(shù)據(jù)來源

本文以荊門市為研究區(qū)域，該市位于湖北省中部，漢江從中穿過。東北、西北多低山丘陵，中、南部為河流、湖泊淤積平原?，F(xiàn)轄5個街道辦事處、50個鎮(zhèn)、2個鄉(xiāng)。本文將城鎮(zhèn)發(fā)展程度較高的荊門市市轄區(qū)中東寶區(qū)泉口街辦、龍泉街辦和掇刀區(qū)白廟街辦合并為東寶區(qū)，其余為掇刀區(qū)，從而得到55個研究單元。人口數(shù)據(jù)來源于荊門市統(tǒng)計局提供的2009年人口統(tǒng)計報表;空間行政邊界、土地利用現(xiàn)狀矢量數(shù)據(jù)來自2009年荊門市1∶50 000土地利用現(xiàn)狀數(shù)據(jù)庫;高程和坡度數(shù)據(jù)來自于空間分辨率為30 m×30 m標準分幅DEM數(shù)據(jù)。城鎮(zhèn)體系和產(chǎn)業(yè)布局數(shù)據(jù)來自荊門市城鄉(xiāng)規(guī)劃局提供的《荊門市城市總體規(guī)劃(2005－2020)》及相關(guān)資料。

2 結(jié)果分析

2.1 網(wǎng)絡節(jié)點結(jié)構(gòu)特征分析

2.1.1 網(wǎng)絡空間結(jié)構(gòu)特征

將荊門市各鄉(xiāng)鎮(zhèn)的鎮(zhèn)政府所在地抽象為節(jié)點，鄉(xiāng)鎮(zhèn)間人口流動強度矩陣設為聯(lián)接邊構(gòu)建荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡。為分析網(wǎng)絡節(jié)點的度分布格局和復雜性特征，首先，對賦值矩陣進行極值標準化處理;然后，經(jīng)過多次敏感性測試，選擇“0.018 8”為切分點進行二值矩陣轉(zhuǎn)化，共得到424條邊;最后，通過UCINET軟件對二值矩陣進行對稱化處理，以ArcGIS10.1為平臺實現(xiàn)網(wǎng)絡結(jié)構(gòu)圖的可視化表達(見圖2)。采用Natural Breaks(Jenks)方法，結(jié)合網(wǎng)絡節(jié)點度的分布特征，將節(jié)點度分為高(16－31)、中(8－14)、低三級(2－7)，如圖2所示。為分析節(jié)點度的地區(qū)差異，將荊門市以漢江為界分為東、西兩區(qū)，網(wǎng)絡空間分布特征如下:

網(wǎng)絡以高節(jié)點度鄉(xiāng)鎮(zhèn)為核心形成“一橫二縱”的空間格局。①高節(jié)點度鄉(xiāng)鎮(zhèn)東寶區(qū)(31)、郢中街辦(22)、新市鎮(zhèn)(16)處于網(wǎng)絡的核心位置，在全市中部形成“一橫”向人口流動網(wǎng)絡密集區(qū)。根據(jù)節(jié)點度的大小和相聯(lián)系鄉(xiāng)鎮(zhèn)的分布，三鄉(xiāng)鎮(zhèn)的輻射區(qū)域具有等級差異。東寶區(qū)節(jié)點度最大，相聯(lián)系鄉(xiāng)鎮(zhèn)在全區(qū)分布較為均勻(東、西區(qū)比例42∶58)，是全市網(wǎng)絡的核心;郢中街辦受漢江影響，多與東部鄉(xiāng)鎮(zhèn)有人口流動(東、西區(qū)比例91∶7)，為東區(qū)網(wǎng)絡的核心;新市鎮(zhèn)地處偏遠且發(fā)展水平有限，輻射范圍限于本縣，與其聯(lián)系鄉(xiāng)鎮(zhèn)中本縣占總數(shù)的75%，屬于縣級網(wǎng)絡核心。②中節(jié)點度鄉(xiāng)鎮(zhèn)主要分布于中部平原地區(qū)，以東寶區(qū)、郢中街辦為核心形成“二縱”格局?！暗谝豢v”，胡集鎮(zhèn)(13)、子陵鋪鎮(zhèn)(8)、掇刀區(qū)(11)、團林鋪鎮(zhèn)(11)、五里鋪鎮(zhèn)(8)、十里鋪鎮(zhèn)(8)以東寶區(qū)為中心沿207國道向南北延伸，形成貫穿全市的人口流動密集帶，涉及鄉(xiāng)鎮(zhèn)最多，影響范圍最大;“第二縱”:豐樂鎮(zhèn)(8)、洋梓鎮(zhèn)(9)、柴湖鎮(zhèn)(14)、舊口鎮(zhèn)(13)以郢中街辦為中心沿漢江向南北延伸。③低節(jié)點度鄉(xiāng)鎮(zhèn)多地處東部以及西部、南部邊緣，東部鄉(xiāng)鎮(zhèn)多低山丘陵，與京山縣縣政府所在地新市鎮(zhèn)形成放射狀網(wǎng)絡格局。西、南部邊緣湖泊、丘陵相間分布，與外界聯(lián)系不便，導致這些地區(qū)的鄉(xiāng)鎮(zhèn)與其他鄉(xiāng)鎮(zhèn)之間聯(lián)系較少，節(jié)點度較低。

2.1.2 網(wǎng)絡空間結(jié)構(gòu)與城鎮(zhèn)體系和產(chǎn)業(yè)布局的相互關(guān)系

圖2 荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡結(jié)構(gòu)圖Fig.2 Network structure of population mobility among towns in Jingmen city

將荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的空間結(jié)構(gòu)圖與《荊門市城市總體規(guī)劃(2005－2020)》中的城鎮(zhèn)體系結(jié)構(gòu)和產(chǎn)業(yè)發(fā)展布局圖進行空間疊加(見圖3)，可以看到兩者的擬合程度較高，在人口流動網(wǎng)絡密集的地方往往培育出了全市的中心、副中心城鎮(zhèn)以及聯(lián)系密切的產(chǎn)業(yè)布局體系。例如:網(wǎng)絡高、中節(jié)點度鄉(xiāng)鎮(zhèn)中東寶區(qū)、郢中街辦、新市鎮(zhèn)、胡集鎮(zhèn)、沙洋鎮(zhèn)的節(jié)點度分布格局與《荊門市城市總體規(guī)劃(2005－2020)》“一主四副”的市域城鎮(zhèn)體系結(jié)構(gòu)相符合;人口流動聯(lián)系“一橫二縱”的分布格局與“一橫三縱”的產(chǎn)業(yè)布局在空間上相匹配:“一橫”對應沿長荊鐵路形成的新興重點產(chǎn)業(yè)發(fā)展帶，“第一縱”、“第二縱”分別對應沿襄荊高速形成的傳統(tǒng)產(chǎn)業(yè)發(fā)展帶、沿漢江形成的生態(tài)產(chǎn)業(yè)發(fā)展帶。需要說明的是沿隨岳高速形成的綜合產(chǎn)業(yè)發(fā)展帶與以新市鎮(zhèn)為核心的局部網(wǎng)絡結(jié)構(gòu)沒有有效疊合。主要原因是隨岳高速公路旁邊沒有國道或者省道通過，道路網(wǎng)絡體系不完整，限制了當?shù)氐娜丝诹鲃有袨椋灾恍纬闪艘孕率墟?zhèn)為核心的放射狀網(wǎng)絡，而沒有出現(xiàn)聯(lián)系密切的帶狀結(jié)構(gòu)。

第一，展現(xiàn)了較強的企業(yè)資源屬性與價值性、稀缺性、難以被復制性等。構(gòu)建了有形和無形的資源，整合了多樣化的生產(chǎn)技能，打造了企業(yè)核心競爭力。第二，人力資源，主要是在知識、技能、能力、稀缺的人力資源本身上，構(gòu)建的核心資源。因此，組織中的稀缺人力資源，屬于核心資源。第三，資源創(chuàng)造核心競爭優(yōu)勢；展現(xiàn)了有形或者無形、人力資源能帶來核心競爭優(yōu)勢［3］。

由此看出，城鎮(zhèn)體系和產(chǎn)業(yè)布局與人口流動網(wǎng)絡的空間結(jié)構(gòu)存在著極強的耦合關(guān)系。鄉(xiāng)鎮(zhèn)對人口的聚集和輻射作用反映出鄉(xiāng)鎮(zhèn)的城市化發(fā)展進程，以及在城鎮(zhèn)體系和產(chǎn)業(yè)布局中的功能和地位，城鎮(zhèn)體系的發(fā)展和產(chǎn)業(yè)布局的形成為改善人們的生產(chǎn)條件和生活水平提供條件，是人口空間流動的內(nèi)在驅(qū)動力;勞動力資源在空間上的重新配置為生產(chǎn)發(fā)展提供人力和技術(shù)資本，是區(qū)域城鎮(zhèn)發(fā)展和產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級的重要支撐系統(tǒng)。

2.2 網(wǎng)絡社團結(jié)構(gòu)特征分析

2.2.1 社團內(nèi)部結(jié)構(gòu)分析

圖3 荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡與城鎮(zhèn)規(guī)劃體系和產(chǎn)業(yè)帶的空間耦合圖Fig.3 Spatial unite of population mobility network among towns and the urban system planning，industrial belt

本研究基于鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡內(nèi)的有效聯(lián)系，采用Guimera 等［24－25］提出的基于模擬退火 算法(Simulated Annealing-Genetic Algorithm)的復雜網(wǎng)絡聚類方法，借助netcarto軟件實現(xiàn)對網(wǎng)絡社團的有效劃分。經(jīng)過多次選代計算，社團指數(shù)最高可達到0.49，此時社團劃分結(jié)果較好，整體網(wǎng)絡劃分為四大社團。各社團內(nèi)部以及社團之間的關(guān)系如圖4所示。

各社團內(nèi)部呈現(xiàn)不同的空間結(jié)構(gòu)特征。受自然、社會、經(jīng)濟等多種因素的綜合影響，各個社團內(nèi)部的聯(lián)系特征差異顯著，根據(jù)影響網(wǎng)絡格局的主要因素可將社團分為不同的指向性集合。①社團1位于荊門市西北部，包含10個鄉(xiāng)鎮(zhèn)，38條內(nèi)部有效聯(lián)接邊，密度0.38。社團1屬于資源指向性網(wǎng)絡，胡集鎮(zhèn)、石橋驛鎮(zhèn)、子陵鋪鎮(zhèn)等礦產(chǎn)資源豐富，工業(yè)企業(yè)發(fā)達，成為社團內(nèi)人口流動的主要集散地。②社團2位于荊門市西南部，包含18個鄉(xiāng)鎮(zhèn)，136條內(nèi)部有效聯(lián)接邊數(shù)，密度0.42。社團2屬于道路指向性網(wǎng)絡，社團內(nèi)部地形平坦、交通發(fā)達，位于交通干道上的鄉(xiāng)鎮(zhèn)聯(lián)系較為密切。③社團3位于中部漢江附近，包含14個鄉(xiāng)鎮(zhèn)，72條內(nèi)部有效聯(lián)接邊，密度0.37。社團3屬于地形指向性網(wǎng)絡，隨著與漢江距離的增加，即從平原到丘陵網(wǎng)絡聯(lián)系逐漸稀疏。④社團4位于荊門市東部，包含13個鄉(xiāng)鎮(zhèn)，78條內(nèi)部有效聯(lián)接邊，密度0.46。社團4屬于政治經(jīng)濟指向性網(wǎng)絡，社團內(nèi)多山地丘陵，人口流動聯(lián)系不便，形成以京山縣縣政府所在地新市鎮(zhèn)為一級中心、經(jīng)濟發(fā)展水平相對較高的羅店鎮(zhèn)、水漋鎮(zhèn)為次核心的縱向網(wǎng)絡結(jié)構(gòu)。

2.2.2 社團外部結(jié)構(gòu)分析

社團內(nèi)部聯(lián)系密切、外部聯(lián)系松散，且聯(lián)系強度具有地區(qū)差異。第一、通過兩兩社團間實際聯(lián)接邊數(shù)與理想聯(lián)接邊數(shù)之比得到各社團間的密度，形成密度矩陣。社團內(nèi)部平均密度0.407遠大于社團間平均密度0.049，說明社團的內(nèi)部聯(lián)系強度明顯高于社團間的聯(lián)系強度，網(wǎng)絡社團結(jié)構(gòu)呈現(xiàn)內(nèi)部聯(lián)系密切、外部對松散的特征。第二、為分析社團間聯(lián)系強度的差異，借助Natural Breaks(Jenks)方法，將密度矩陣值大于0.049定義為強聯(lián)系，反之為弱聯(lián)系。從圖4可以看出:強聯(lián)系主要存在于相鄰社團之間，說明距離較近的鄉(xiāng)鎮(zhèn)之間人口流動強度較大。社團之間的強聯(lián)系體現(xiàn)了網(wǎng)絡內(nèi)部的開放性，有利于不同指向性社團的交流與合作，整體推動地區(qū)發(fā)展;弱聯(lián)系則分別聯(lián)接社團1與社團3、4，橫跨漢江(需要說明的是社團1—社團4之間也是跨江聯(lián)系，但位于郢中街辦的鐘祥漢江大橋加強了大橋附近的漢江兩岸鄉(xiāng)鎮(zhèn)之間的聯(lián)系)，空間距離較遠。遠距離鄉(xiāng)鎮(zhèn)之間相互聯(lián)系，有利于提高網(wǎng)絡的聯(lián)接廣度和整體聯(lián)通效率;社團2與社團3之間不存在聯(lián)系，主要受到空間距離的影響以及河流阻隔，且兩社團的經(jīng)濟發(fā)展水平相差較小，整體聯(lián)接能力有限。

圖4 荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡社團結(jié)構(gòu)圖Fig.4 The Community structure of population mobility among towns in Jingmen City

2.3 復雜網(wǎng)絡特征分析

2.3.1 無標度特征分析

對度累計概率分布與度數(shù)進行函數(shù)擬合(見圖5)發(fā)現(xiàn):冪分布函數(shù)為:P(k)=7.80k－1.65，R2=0.89，說明網(wǎng)絡的度累計概率分布與度數(shù)呈現(xiàn)一定的冪律關(guān)系，但是指數(shù)1.65＜2，無標度特征不夠明顯。網(wǎng)絡的整體趨勢更符合指數(shù)分布(P(k)=1.26e－0.16k，R2=0.92)。在網(wǎng)絡節(jié)點連接中，冪律分布和指數(shù)分布都表示度分布的不均衡特征，分布概率從低節(jié)點度向高節(jié)點度遞減，但是指數(shù)分布的遞減速度較為緩慢，即節(jié)點度分布在高、低節(jié)點度之間有相對較多的中節(jié)點度作為過渡，使度分布的不均衡特征沒有冪律分布明顯。節(jié)點連接的整體擇優(yōu)行為由于一定客觀因素的制約，只能在局部實現(xiàn)相對擇優(yōu)，從而出現(xiàn)一定數(shù)量的節(jié)點度處于中間水平，使節(jié)點度分布偏向于指數(shù)特征［7］。從圖4可以看出，在荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡中，度分布概率從低節(jié)點度(2－7)向高節(jié)點度(16－31)的遞減現(xiàn)象明顯，但是中節(jié)點度(8－14)的度分布相對集中，遞減速度較為緩慢，說明荊門市鄉(xiāng)鎮(zhèn)間的人口流動以東寶區(qū)、郢中街辦等高水平發(fā)展鄉(xiāng)鎮(zhèn)為核心，同時受空間距離、地形地貌等因素的影響呈現(xiàn)局部擇優(yōu)特征，目的地選擇傾向于一定范圍中，即所在社團內(nèi)的核心節(jié)點鄉(xiāng)鎮(zhèn)。荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡距離具有全局擇優(yōu)性質(zhì)的無標度網(wǎng)絡仍有一定距離，在城鎮(zhèn)體系規(guī)劃中需要采取有效措施改善全市范圍內(nèi)的人口流動的基礎設施，同時積極培育全市的“增長極”或“發(fā)展極”鄉(xiāng)鎮(zhèn)，在提高其社會經(jīng)濟發(fā)展水平的基礎上增強其對全市的輻射帶動作用。

圖5 荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡節(jié)點度累計概率分布Fig.5 Cumulative probability distribution of degree in population mobility network

2.3.2 小世界特征分析

(1)平均路徑距離分析。荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的平均路徑距離為2.34，即在網(wǎng)絡中任一節(jié)點到其他節(jié)點需要經(jīng)過2.34條邊。通過最短路徑距離統(tǒng)計分析，荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡中14.3%的鄉(xiāng)鎮(zhèn)之間最短路徑距離為1，主要存在于發(fā)展條件較好、交通便利的鄉(xiāng)鎮(zhèn)與其他鄉(xiāng)鎮(zhèn)，以及鄰近鄉(xiāng)鎮(zhèn)之間，這些鄉(xiāng)鎮(zhèn)一般有道路直接相通;43.8%的鄉(xiāng)鎮(zhèn)之間最短路徑距離為2，此類路徑多位于漢江一邊的鄉(xiāng)鎮(zhèn)之間，少數(shù)跨江路徑主要與漢江附近的發(fā)達鄉(xiāng)鎮(zhèn)相聯(lián)系;36.2%的鄉(xiāng)鎮(zhèn)之間最短路徑距離為3，主要為跨江路徑，且多連接網(wǎng)絡空間邊緣鄉(xiāng)鎮(zhèn);5.6%的鄉(xiāng)鎮(zhèn)之間最短路徑距離為4，其中52%存在于西北部仙居鄉(xiāng)、東北部綠林鎮(zhèn)、南部官垱鎮(zhèn)、李市鎮(zhèn)、毛李鎮(zhèn)四個鄉(xiāng)鎮(zhèn)與其他邊遠鄉(xiāng)鎮(zhèn)之間;0.2%的鄉(xiāng)鎮(zhèn)之間最短路徑距離為5，全部發(fā)生在綠林鎮(zhèn)與官垱鎮(zhèn)、李市鎮(zhèn)、毛李鎮(zhèn)之間。由統(tǒng)計結(jié)果可知，58.1%的鄉(xiāng)鎮(zhèn)之間最多經(jīng)過1個節(jié)點與其他鄉(xiāng)鎮(zhèn)有人口流動聯(lián)系，94.2%的鄉(xiāng)鎮(zhèn)之間則最多需要經(jīng)過2個節(jié)點，較長的最短路徑距離只集中于少數(shù)邊遠鄉(xiāng)鎮(zhèn)之間，荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的傳輸效率較高。網(wǎng)絡傳輸?shù)亩搪窂街饕獨w結(jié)于網(wǎng)絡中相距較遠的節(jié)點之間存在一定的“捷徑”，例如:東寶區(qū)—新市鎮(zhèn)、郢中街辦—水漋鎮(zhèn)等，使網(wǎng)絡整體節(jié)點之間的分離程度較小，有利于鄉(xiāng)鎮(zhèn)之間的人口流動以及以人口流動為載體的信息擴散、技術(shù)傳播等，以此促進鄉(xiāng)鎮(zhèn)之間的交流與合作，從而帶動欠發(fā)達鄉(xiāng)鎮(zhèn)的城市化進程。

(2)簇系數(shù)分析。荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的簇系數(shù)為0.66，表現(xiàn)出一定的局部集聚性，表明網(wǎng)絡中各個節(jié)點之間形成短距離聯(lián)系的可能性比較大。整體網(wǎng)絡的簇系數(shù)與節(jié)點度的統(tǒng)計結(jié)果(見圖6)進行比較分析發(fā)現(xiàn):①高、中節(jié)點度鄉(xiāng)鎮(zhèn)的簇系數(shù)相對較小，平均簇系數(shù)為0.48。其中東寶區(qū)、郢中街辦、新市鎮(zhèn)的簇系數(shù)分別為0.21、0.27、0.25。②低節(jié)點度鄉(xiāng)鎮(zhèn)的簇系數(shù)較高，兩極分化嚴重，平均簇系數(shù)為0.69。其中10個鄉(xiāng)鎮(zhèn)簇系數(shù)為1，平均度數(shù)僅3.4，主要位于山地丘陵和邊遠地區(qū);楊集鎮(zhèn)的簇系數(shù)為0，其節(jié)點度為2，與其相連的客店鎮(zhèn)與新市鎮(zhèn)之間距離較遠，且有山地阻隔，兩者人口流動強度較小，沒有達到切分點水平。以上現(xiàn)象主要歸結(jié)于兩點:第一、高、中節(jié)點度鄉(xiāng)鎮(zhèn)相對于低節(jié)點度鄉(xiāng)鎮(zhèn)一般社會經(jīng)濟發(fā)展水平較高，輻射范圍較大，因為高、中節(jié)點度鄉(xiāng)鎮(zhèn)與低度節(jié)點之間有較大“勢能”(經(jīng)濟因素、人口因素等)存在，在可達性條件較好的情況下易與距離較遠的節(jié)點直接建立聯(lián)系。第二、節(jié)點聯(lián)接由于受到空間距離及地形等條件的約束傾向于局部擇優(yōu)，在較小區(qū)域內(nèi)聯(lián)系緊密，而在較大區(qū)域內(nèi)較為松散。因此，基于高、中節(jié)點度鄉(xiāng)鎮(zhèn)與低節(jié)點度鄉(xiāng)鎮(zhèn)的輻射范圍差異以及節(jié)點聯(lián)接性質(zhì)，網(wǎng)絡簇系數(shù)表現(xiàn)為高、中節(jié)點度鄉(xiāng)鎮(zhèn)低于低節(jié)點度鄉(xiāng)鎮(zhèn)。

圖6 荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的簇系數(shù)、節(jié)點度Fig.6 Statistics of the clustering coefficient and degree in population mobility network

(3)小世界特征判斷。與具有相同節(jié)點和平均度的隨機網(wǎng)絡比較，實際網(wǎng)絡與隨機網(wǎng)絡的平局路徑距離的比值越接近La(L/Lr)1，且簇系數(shù)比值Ca(C/Cr)越大于1，則該網(wǎng)絡的小世界特征越明顯，其中L、C、Lr、Cr分別表示實際網(wǎng)絡和隨機網(wǎng)絡的平均路徑距離(Lr=lnN/ln＜k＞，N、＜k＞分別表示實際網(wǎng)絡的節(jié)點總數(shù)、平均節(jié)點度)、簇系數(shù)(Cr≈ ＜k＞/N)［26］。根據(jù)指標計算可知:La≈1.19 接近1且Ca≈4.71明顯大于1。綜上所述，荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡具有明顯的小世界網(wǎng)絡特性，即該網(wǎng)絡具有較好的聯(lián)通性和強集聚性［5］，具體表現(xiàn)為:任意兩個鄉(xiāng)鎮(zhèn)之間發(fā)生人口流動聯(lián)系最多只需要中轉(zhuǎn)其他1－2個鄉(xiāng)鎮(zhèn)，網(wǎng)絡整體流動順暢;除此之外，較高的族系數(shù)表現(xiàn)為網(wǎng)絡局部聯(lián)系緊密，人口流動多發(fā)生在小范圍內(nèi)，特別是相鄰鄉(xiāng)鎮(zhèn)間，鄉(xiāng)鎮(zhèn)間發(fā)生人口流動關(guān)聯(lián)關(guān)系無需過多中介。

3結(jié)論

(1)本文將復雜網(wǎng)絡分析的理論和方法引入人口流動空間格局研究，借助基于空間可達性的輻射模型測度鄉(xiāng)鎮(zhèn)間的人口流動強度，以鄉(xiāng)鎮(zhèn)為節(jié)點，鄉(xiāng)鎮(zhèn)間的人口流動強度為聯(lián)接邊構(gòu)建了鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡，實現(xiàn)了在異質(zhì)性地域格局下，各鄉(xiāng)鎮(zhèn)在人口流動網(wǎng)絡中的聯(lián)系及其空間分布的可視化表達，并通過度分布、平均路徑距離、簇系數(shù)等指標探討網(wǎng)絡的無標度、小世界等復雜網(wǎng)絡特性。本研究有助于為區(qū)域城鎮(zhèn)體系規(guī)劃和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整提供決策依據(jù)。

(2)從網(wǎng)絡分析的結(jié)果來看:荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡以東寶區(qū)、郢中街辦、新市鎮(zhèn)為核心，呈“一橫兩縱”的空間分布格局，與荊門市城鎮(zhèn)體系和產(chǎn)業(yè)布局存在極強的耦合關(guān)系;全網(wǎng)分為四大社團，受不同因素的綜合作用，社團內(nèi)部的結(jié)構(gòu)特征各異，社團間聯(lián)系較為松散，且相鄰社團間呈現(xiàn)強聯(lián)系;受到空間距離和地形地貌等條件的影響，網(wǎng)絡的無標度冪指數(shù)較正常值偏小，度分布的指數(shù)特征較為顯著，荊門市鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的節(jié)點聯(lián)接具有局部擇優(yōu)性質(zhì);網(wǎng)絡具有較小的平均路徑距離和較高的簇系數(shù)，小世界網(wǎng)絡特性明顯，呈現(xiàn)較好的整體連通性和較強的局部聚集性。

人口流動網(wǎng)絡影響因素的多樣性增加了網(wǎng)絡空間結(jié)構(gòu)的復雜程度，本文雖然考慮了道路、用地類型、地形地貌等影響因子，但并未融入經(jīng)濟發(fā)展水平、主體決策行為等條件的思考，對人口流動網(wǎng)絡微觀機理的描述需要進一步加強;其次，人口流動網(wǎng)絡的形成和演化是一個連續(xù)的動態(tài)過程，由于數(shù)據(jù)限制本文的分析基礎主要是斷面數(shù)據(jù)，缺乏基于多年份數(shù)據(jù)的網(wǎng)絡時空演化過程和發(fā)展規(guī)律研究;再次，本文研究了網(wǎng)絡的節(jié)點和社團結(jié)構(gòu)以及無標度、小世界等靜態(tài)特征，沒有對網(wǎng)絡的動力學特征展開深入討論，這也是需要繼續(xù)改進的地方。在以后的研究中，需融入地區(qū)的社會經(jīng)濟屬性和流動人口的主體決策等因素，同時加強人口流動網(wǎng)絡變化的縱向?qū)Ρ妊芯浚瑥撵o態(tài)和動態(tài)相結(jié)合的復雜網(wǎng)絡評價視角更好地揭示鄉(xiāng)鎮(zhèn)人口流動網(wǎng)絡的空間格局和演化規(guī)律。

References)

［1］朱傳耿，顧朝林，馬榮華，等.中國流動人口的影響要素與空間分布［J］.地理學報，2001，56(5):548－559.［Zhu Chuangeng，Gu Chaolin，Ma Ronghua，et al.The Influential Factors and Spatial Distribution of Floating Population in China［J］.Acta Geographica Sinica，2001，56(5):548 －559.］

［2］段平忠.中國省際間人口遷移對經(jīng)濟增長動態(tài)收斂的影響［J］.中國人口·資源與環(huán)境，2011，18(12):146－152.［Duan Pingzhong.Influence of China’s Population Flow in the Change of Regional Disparity Since［J］.China Population Resources and Environment，2011，18(12):146 －152.］

［3］劉盛和，鄧羽，胡章.中國流動人口地域類型的劃分方法及空間分布特征［J］.地理學報，2010，65(10):1187－1197.［Liu Shenghe，Deng Yu，Hu Zhang.Deng Yu，Hu Zhang.Research on Classification Methods and Spatial Patterns of the Regional Types of China’s Floating Population［J］.Acta Geographica Sinica，2010，65(10):1187 －1197.］

［4］俞憲忠.人口流動規(guī)律及其政策含義［J］.中國人口·資源與環(huán)境，2005，15(1):120－124.［Yu Xianzhong.Population Flow Law And Its Policy Connotation［J］.China Population Resources and Environment，2005，15(1):120 － 124.］

［5］Watts D J，Strogatz S H.Collective dynamics of‘small-world’networks［J］.1998，393(6684):440 － 442.

［6］Barabási A，Albert R.Emergence of Scaling in Random Networks［J］.Science，1999，286(5439):509 －512.

［7］Li X，Chen G.A local－ world evolving network model［J］.Physica A:Statistical Mechanics and its Applications，2003，328(1 – 2):274－286.

［8］Zhang Z，Rong L，Guo C.A deterministic small－ world network created by edge iterations［J］.Physica A:Statistical Mechanics and its Applications，2006，363(2):567 －572.

［9］Prettejohn B J，Berryman M J，Mcdonnell M D.Methods for generating complex networks with selected structural properties for simulations:a review and tutorial for neuroscientists［J］.Frontiers in computational neuroscience 2011，5(11):1 －18.

［10］林敏，李南，陳婷婷.基于復雜網(wǎng)絡的知識轉(zhuǎn)移模擬與分析［J］.系統(tǒng)工程，2009，27(3):115 －118.［Lin Mei，Li Nan，Chen Tingting.Simulation and Analysis of Knowledge Transfer Based on Complex Network Theory［J］.Systems Engineering，2009，27(3):115 －118.］

［11］Borgatti S P，Mehra A，Brass D J，et al.Network Analysis in the Social Sciences［J］.Science，2009，323(5916):892 －895.

［12］Villas Boas P R，Rodrigues F A，Da F.Costa L.Modeling worldwide highway networks［J］.Physics Letters A，2009，374(1):22－27.

［13］羅家德.社會網(wǎng)分析講義［M］.北京:社會科學文獻出版社，2010.［Luo Jiade.Social Network Analysis［M］.Beijing:Social Sciences Academic Press，2010.］

［14］朱杰.長江三角洲人口遷移空間格局、模式及啟示［J］.地理科學進展，2009，28(3):353 －361.［Zhu Jie.Spatial Structure，Mode and the Meaning of Population Migration in the Yangtze River［J］.Progress in Geography，2009，28(3):353 －361.］

［15］Cheung C，Black J.A Reappraisal of the Intervening Opportunities Model of Commuter Behaviour［J］.Road and Transport Research，2008，17(2):3 －18.

［16］Bastin F，Cirillo C，Toint P L.Estimating Nonparametric Random Utility Models with an Application to the Value of Time in Heterogeneous Populations［J］.Transportation science，2010，44(4):537－549.

［17］Simini F，Gonzalez M C，Maritan A，et al.A universal model for mobility and migration patterns［J］.Nature，2012，484(7392):96－100.

［18］劉少麗，顧小平，裴友法，等.徐州市避震疏散場所的可達性與公平性［J］.經(jīng)濟地理，2012，32(3):69－74.［Liu Shaoli，Gu Xiaoping，Pei Youfa，et al.Accessibility and Equity of Seismic Shelters for Evacuation in Xuzhou［J］.Economic Geography，2012，32(3):69 －74.］

［19］劉俊，陸玉麒.江蘇省公路交通網(wǎng)絡可達性評價研究［J］.南京師大學報(自然科學版)，2008，31(3):129－134.［Liu Jun，Lu Yuqi.Appraisement of Road Transportation Network Accessibility in Jiangsu Province［J］.Journal of Nanjing Normal University:Natural Science Edition，2008，31(3):129 －134.］

［20］鐘業(yè)喜，陸玉麒.基于空間聯(lián)系的城市腹地范圍劃分——以江蘇省為例［J］.地理科學，2012，29(5):536－543.［Zhong Yexi，Lu Yuqi.Measuring Method of Urban Hinterland Based on Spatial Linkage:A Case of Jiangsu Province［J］. Scientia Geographica Sinica，2012，29(5):536 －543.］

［21］蔣海兵，徐建剛，祁毅.京滬高鐵對區(qū)域中心城市陸路可達性影響［J］.地理學報，2010，65(10):1287－1298.［Jiang Haibin，Xu Jiangang，Qi Yi.The Influence of Beijing－ Shanghai High －speed Railways on Land Accessibility of Regional Center Cities［J］.Acta Geographica Sinica，2010，65(10):1287 －1298.］

［22］Newman M.The Structure and Function of Complex Networks［J］.SIAM Review，2003，45(2):167 －256.

［23］郭世澤，路哲明.復雜網(wǎng)絡基礎理論［M］.Beijing:科學出版社，2012:267－268.［Guo Shize，Lu Zheming.Basic theory of Complex Network［M］.Beijing:Science Press，2012:267 －268.］

［24］Guimera R，Amaral L A.Functional cartography of complex metabolic networks［J］.2005，433(7028):895 －900.

［25］Guimera R，Sales－pardo M，Amaral L A N.Classes of complex networks defined by role－to－role connectivity profiles［J］.2007，3(1):63－69.

［26］Brian U，Jarrett S.Collaboration and Creativity:The Small World Problem［J］.American Journal of Sociology，2005，111(2):447 －504.

Study on Population Mobility Network among Towns Based on Improved Radiation Model

HOU He-ping1，2 LIU Yan-fang1，2 LI Ji-wei1，2 KONG Xue-song1，2
(1.School of Resource and Environment Science，Wuhan University，Wuhan Hubei 430079，China;2.Key Laboratory of Geographic Information System，Ministry of Education，Wuhan University，Wuhan Hubei 430079，China)

It is of great significance for the regional urban system planning and industrial structure adjustment to explore the relationship and spatial structure of population mobility on the scale of town.The paper，from the perspective of complex network theory，by introducing spatial accessibility to traditional radiation and correcting parameter，builds a network of population mobility among towns.From the perspectives of node，community structure，scale free and small world characteristics，the paper also discusses the spatial distribution pattern and complex network characteristics of population mobility under the influence of regional differences.On basis of the empirical research of 55 towns in Jingmen City，the results reveal that:firstly，the node degree distribution of population mobility network of towns in Jingmen City centers around the towns who are relatively developed，presents a pattern of spatial disparities with the characteristics of“one vertical line and two horizontal lines”，and shows an obvious couple relationship with the urban system and industrial layout of Jingmen City.Secondly，the whole network consists of four communities whose interior spatial structures are close and different from each other and the relationships among them are relatively loose.Thirdly，impacted by the spatial distance and geographic environment，scale free charateristics are not obvious.Fourthly，to certain extent，the network bears the characteristic of small world whose overall connectivity and local clustering are obvious.By improving the traditional model，the paper realizes the visual presentation of the relationship and distribution of township population mobility in the heterogeneous space，and the results of this empirical research are of great value in decision making of Jingmen’s town planning and industrial layout.

population mobility;complex network;radiation model;accessibility;Jingmen City

K901.3;C922

A

1002 －2104(2013)08 －0107 －09

10.3969/j.issn.1002 －2104.2013.08.016

2013－05－09

侯賀平，博士生，主要研究方向為城鄉(xiāng)區(qū)域規(guī)劃。

劉艷芳，博士，教授，博導，主要研究方向為區(qū)域規(guī)劃與地理信息應用工程。

國家“十二五”科技支撐計劃資助項目“村鎮(zhèn)區(qū)域空間規(guī)劃與集約發(fā)展關(guān)鍵技術(shù)研究”(編號:2012BAJ22B00)。

(編輯:于 杰)