產(chǎn)業(yè)空間集聚中心引力指數(shù)算法的設(shè)計(jì)及應(yīng)用
——基于長(zhǎng)三角一體化視角的企業(yè)微觀數(shù)據(jù)驗(yàn)證

嵇正龍，宋 宇

(西北大學(xué) 經(jīng)濟(jì)管理學(xué)院，陜西 西安 710127)

一、問(wèn)題提出

促進(jìn)經(jīng)濟(jì)增長(zhǎng)和發(fā)展的政策措施效果在很大程度上取決于經(jīng)濟(jì)活動(dòng)在地理上的空間集聚。長(zhǎng)三角一體化、國(guó)家中心城市和戶(hù)籍制度改革等國(guó)家戰(zhàn)略規(guī)劃決策，意圖進(jìn)一步促進(jìn)資源的集聚，形成經(jīng)濟(jì)增長(zhǎng)的新動(dòng)能。根據(jù)2018年各省市國(guó)民經(jīng)濟(jì)與社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)，長(zhǎng)江三角洲地區(qū)(滬、蘇、浙、皖)在約35.91萬(wàn)平方公里的國(guó)土面積上，實(shí)現(xiàn)地區(qū)生產(chǎn)總值 21.15萬(wàn)億元，工業(yè)增加值 8.83萬(wàn)億元，常駐人口2.25億人，分別占全國(guó)的3.74%、23.49%、28.95%、16.15%。從國(guó)土面積占比與三項(xiàng)經(jīng)濟(jì)社會(huì)發(fā)展的核心指標(biāo)的全國(guó)占比形成的反差，可以看出長(zhǎng)三角地區(qū)經(jīng)濟(jì)活動(dòng)的空間集聚程度較高?？傮w來(lái)說(shuō)，長(zhǎng)三角地區(qū)經(jīng)濟(jì)空間集聚慢于人口空間集聚，創(chuàng)新的空間集聚和溢出效應(yīng)表現(xiàn)出明顯的區(qū)域內(nèi)部差異[1]，生產(chǎn)者服務(wù)業(yè)發(fā)展和密集的區(qū)域高鐵網(wǎng)絡(luò)等對(duì)制造業(yè)的空間集聚產(chǎn)生了顯著的正向促進(jìn)作用[2]。

國(guó)內(nèi)關(guān)于地區(qū)經(jīng)濟(jì)空間集聚的測(cè)度大多是對(duì)傳統(tǒng)方法的應(yīng)用，也有研究文獻(xiàn)是利用微觀企業(yè)層面的數(shù)據(jù)應(yīng)用較為前沿的DO指數(shù)測(cè)度空間集聚特征[3]，而對(duì)空間集聚方法本身缺乏更深入的探討。已有空間集聚衡量方法的發(fā)展脈絡(luò)是怎么樣的？有哪些優(yōu)點(diǎn)和不足？能否設(shè)計(jì)一個(gè)新的衡量方法對(duì)已有方法進(jìn)行有益的補(bǔ)充？為回答這些問(wèn)題，本文聚焦于經(jīng)濟(jì)活動(dòng)空間集聚衡量新方法的設(shè)計(jì)。

二、文獻(xiàn)回顧

自產(chǎn)業(yè)集聚理論發(fā)展之初，經(jīng)濟(jì)集聚程度的測(cè)度方法的研究與應(yīng)用就成為了經(jīng)濟(jì)集聚理論發(fā)展的核心組成部分之一。回顧經(jīng)濟(jì)集聚測(cè)度方法的相關(guān)研究成果，主要可以分為基于產(chǎn)業(yè)集群的離散空間測(cè)度法和基于點(diǎn)對(duì)距離的連續(xù)空間測(cè)度法兩類(lèi)，并有學(xué)者對(duì)二者的融合做出了探索?？偟膩?lái)說(shuō)，兩類(lèi)空間集聚程度測(cè)度方法共存于現(xiàn)有研究文獻(xiàn)中。

(一)基于產(chǎn)業(yè)集群的離散空間集聚指數(shù)

該類(lèi)指數(shù)，基本沿著馬歇爾提出的外部性理論和規(guī)模收益遞增解釋集聚，至關(guān)重要的經(jīng)濟(jì)集聚的空間特征被簡(jiǎn)化處理。假設(shè)空間由可識(shí)別的離散地理單元構(gòu)成，基于產(chǎn)業(yè)集群視角，應(yīng)用產(chǎn)值或者人口等數(shù)據(jù)構(gòu)造集聚度指標(biāo)，試圖將空間維度作為經(jīng)濟(jì)變量引入主流分析框架。

早期的研究，由于公司的確切位置數(shù)據(jù)的獲取限制，所有關(guān)于經(jīng)濟(jì)活動(dòng)地理集中的研究都使用行政區(qū)域的匯總數(shù)據(jù)來(lái)測(cè)量空間集聚。這方面的研究有大量的文獻(xiàn)，主要涉及區(qū)域比較，但公司或研究機(jī)構(gòu)等經(jīng)濟(jì)行為主體的實(shí)際空間位置無(wú)從體現(xiàn)。這類(lèi)指數(shù)以熵指數(shù)、赫芬達(dá)爾指數(shù)、空間基尼系數(shù)、泰勒指數(shù)、EG指數(shù)等為代表，將經(jīng)濟(jì)活動(dòng)空間按照國(guó)家或者行政管理區(qū)劃等分為若干區(qū)域，評(píng)估經(jīng)濟(jì)活動(dòng)的空間集中程度[4]。其中，EG指數(shù)與其他指標(biāo)不同之處在于，通過(guò)構(gòu)造產(chǎn)業(yè)集聚的三因素模型，考慮企業(yè)規(guī)?？傮w分布，對(duì)基礎(chǔ)數(shù)據(jù)進(jìn)行處理，然后計(jì)算出產(chǎn)業(yè)集聚新指標(biāo)[5]，對(duì)集聚判斷更為準(zhǔn)確。該類(lèi)指標(biāo)的優(yōu)點(diǎn)表現(xiàn)為，基于經(jīng)濟(jì)集聚的理論分析，將容易獲取的宏觀數(shù)據(jù)或區(qū)域經(jīng)濟(jì)總量指標(biāo)，應(yīng)用簡(jiǎn)單易懂的公式估計(jì)經(jīng)濟(jì)集聚程度，尤其是計(jì)算量較小。因此，此類(lèi)指標(biāo)的應(yīng)用，不僅在早期的實(shí)證研究中較為常見(jiàn)，當(dāng)前的研究中依然被廣泛使用。

但是該類(lèi)指標(biāo)忽略了空間因素的關(guān)鍵缺陷也是顯而易見(jiàn)的。第一，忽略了空間因素。作為衡量經(jīng)濟(jì)活動(dòng)空間集聚的指標(biāo)卻沒(méi)有考慮地理距離的影響，而地理距離正是知識(shí)溢出的空間局限性來(lái)源。盡管技術(shù)進(jìn)步和社會(huì)發(fā)展了，但是知識(shí)溢出仍然受到空間距離的限制。也正是因?yàn)楹?jiǎn)單化的處理，這些指數(shù)反映的是某個(gè)單元的產(chǎn)業(yè)集聚，而非真正意義上的空間集聚，同時(shí)缺乏判斷一個(gè)地區(qū)的專(zhuān)門(mén)化程度是否存在集群的閾值，被認(rèn)為更多的是反映專(zhuān)業(yè)化的趨勢(shì)，不能確定集群的空間維度，只能確定專(zhuān)業(yè)化程度較高的區(qū)域。第二，高估了集聚的作用。由于該類(lèi)指標(biāo)反映的是產(chǎn)值在較大的地理尺度上的相對(duì)集中，而產(chǎn)值與經(jīng)濟(jì)增長(zhǎng)等發(fā)展指標(biāo)具有顯著的同向性，即具有較強(qiáng)的內(nèi)生性，難以通過(guò)控制變量解決，因此會(huì)正向高估集聚程度及其對(duì)經(jīng)濟(jì)增長(zhǎng)的影響。第三，集聚程度受限于區(qū)域劃分。這類(lèi)指標(biāo)所確定的產(chǎn)業(yè)集聚程度取決于區(qū)域邊界的主觀劃分，即分區(qū)邊界的位置對(duì)觀測(cè)水平有影響[6]，產(chǎn)生所謂的區(qū)域單位可修改問(wèn)題，簡(jiǎn)稱(chēng)MAUP問(wèn)題。理想情況下，應(yīng)該在連續(xù)的空間中進(jìn)行所有的數(shù)據(jù)分析，以避免主觀的空間單元?jiǎng)澐侄a(chǎn)生的邊界偏差[7]。有些研究假設(shè)離散測(cè)度的空間集聚程度不受區(qū)域排列的影響，嘗試將離散測(cè)度與自相關(guān)測(cè)度相結(jié)合限制MAUP的作用；有些研究試圖通過(guò)整合空間集聚程度指數(shù)的自相關(guān)程度，在一定程度上糾正空間集聚程度結(jié)果。前者是試圖尋找限制MAUP效應(yīng)的方法，后者是遵循MAUP的要求，不做任何的空間分區(qū)，但都不是很成功[8]。因此，關(guān)注地理指標(biāo)因素，發(fā)展新的空間集聚程度指數(shù)是非常有必要的。

(二) 基于點(diǎn)對(duì)距離的連續(xù)空間集聚指數(shù)

基于點(diǎn)對(duì)距離的連續(xù)空間集聚測(cè)度方法，假設(shè)空間是連續(xù)的非均勻狀態(tài)，不可分割，基于地理空間距離，描述經(jīng)濟(jì)活動(dòng)在空間上不同單元的聯(lián)合行為所體現(xiàn)的整體分布特征，并借助于經(jīng)典假設(shè)對(duì)該分布特征進(jìn)行統(tǒng)計(jì)檢驗(yàn)。

隨著數(shù)據(jù)的收集和處理技術(shù)的進(jìn)步，空間經(jīng)濟(jì)學(xué)的實(shí)證研究需要高效利用越來(lái)越容易獲取的詳實(shí)的微觀數(shù)據(jù)的新工具。同時(shí)，為解決經(jīng)濟(jì)學(xué)家傳統(tǒng)上使用的空間集聚衡量方法和指標(biāo)存在的缺陷，地理位置信息成為構(gòu)造衡量集聚程度新指標(biāo)的重要數(shù)據(jù)來(lái)源?；诰嚯x的方法，利用微觀經(jīng)濟(jì)數(shù)據(jù)，將每個(gè)企業(yè)視為地圖上的一個(gè)點(diǎn)，并借用所謂的點(diǎn)模式分析方法研究它們的空間分布，其中最具有開(kāi)創(chuàng)性意義的是使用點(diǎn)對(duì)數(shù)據(jù)的K函數(shù)測(cè)度聚類(lèi)[9]，成為后來(lái)大部分基于距離衡量經(jīng)濟(jì)活動(dòng)空間集聚指標(biāo)設(shè)計(jì)的起點(diǎn)。但是，直到DO指數(shù)方法的出現(xiàn)[10]，基于距離的經(jīng)濟(jì)活動(dòng)空間集聚程度的衡量方法才得到迅速推廣。DO指數(shù)開(kāi)創(chuàng)性的提出了一種基于公司間距離分布的集聚度量方法，不依賴(lài)分區(qū)，提供關(guān)于所有尺度的集聚分布的信息。DO指數(shù)很好的發(fā)展了K函數(shù)方法，將企業(yè)的地理位置信息納入空間集聚的衡量研究中，應(yīng)用核密度函數(shù)對(duì)經(jīng)濟(jì)活動(dòng)的集聚程度做出估計(jì)和檢驗(yàn)，成為后來(lái)同類(lèi)研究的典范?；贒O指數(shù)，相繼發(fā)展出具有基準(zhǔn)的累積指數(shù)版本的BT指數(shù)；基于K函數(shù)的累積密度函數(shù)估計(jì)的MP指數(shù)等代表性方法[11-12]。

上述指數(shù)研究利用微觀數(shù)據(jù)，把企業(yè)看作是經(jīng)濟(jì)空間中無(wú)量綱的點(diǎn)。但在經(jīng)濟(jì)空間中觀察到的點(diǎn)并非無(wú)量綱，企業(yè)的員工數(shù)量、產(chǎn)品、資本等方面具有不同的量綱特征?；诖?，利用點(diǎn)過(guò)程理論的工具，采用基于模型的方法，使用基于K函數(shù)加權(quán)版本的指數(shù)考慮企業(yè)特征指標(biāo)構(gòu)建了EGA指數(shù)[13]。DO指數(shù)以及MP指數(shù)和EGA指數(shù)都是受到了標(biāo)記點(diǎn)過(guò)程理論的啟發(fā)，但是只有EGA指數(shù)明確地將產(chǎn)業(yè)集聚測(cè)量與一個(gè)確定的統(tǒng)計(jì)參數(shù)聯(lián)系起來(lái)，可以不使用蒙特卡羅方法來(lái)虛擬檢驗(yàn)[14]。此外，受到波特的產(chǎn)業(yè)集群理論啟示，使用公司級(jí)聚類(lèi)指數(shù)檢測(cè)空間聚類(lèi)進(jìn)行區(qū)域研究，在計(jì)算和解釋上都不同于現(xiàn)有的基于距離的度量，是識(shí)別沒(méi)有預(yù)先確定邊界的集群的一種方法，被稱(chēng)為Cluster指數(shù)法[15]。

針對(duì)上述方法都是計(jì)算歐幾里德直線(xiàn)距離，無(wú)法跨越自然屏障的缺陷，部分研究成果采用了規(guī)劃距離對(duì)上述的一些方法進(jìn)行了優(yōu)化。關(guān)于距離的測(cè)算：有定義并研究了Network K函數(shù)，分析沿著街道網(wǎng)絡(luò)測(cè)量的距離；也有考慮現(xiàn)實(shí)世界的地理位置，需要道路聯(lián)結(jié)的街道和橋梁穿過(guò)的河流形成的路徑[16]。點(diǎn)過(guò)程的數(shù)學(xué)性質(zhì)通常依賴(lài)于歐幾里德距離，因此總的來(lái)說(shuō)還是歐幾里德距離的計(jì)算方法較多。

基于距離的集聚程度衡量方法的優(yōu)點(diǎn)是具有統(tǒng)計(jì)理論支撐，能夠在全域空間上檢測(cè)地區(qū)所發(fā)生顯著地理集中或分散的尺度。以DO指數(shù)為代表，基于空間距離的指數(shù)研究都構(gòu)造了統(tǒng)計(jì)檢驗(yàn)，比傳統(tǒng)的空間集聚程度測(cè)量方法更為穩(wěn)健，成為了當(dāng)前基于距離的空間集聚衡量的評(píng)估方法的主流和前沿。與基于集群離散空間的衡量方法不同，基于距離的方法不以特定的方式對(duì)問(wèn)題區(qū)域進(jìn)行分區(qū)，而是考慮連續(xù)空間。與只在單一尺度上描述經(jīng)濟(jì)活動(dòng)位置的測(cè)量方法不同，基于距離的方法可以在任意尺度上檢測(cè)空間結(jié)構(gòu)。

基于點(diǎn)對(duì)距離的連續(xù)空間集聚指數(shù)法，存在企業(yè)異質(zhì)性考察不足、置信區(qū)間設(shè)置較為隨意、集聚程度的具體數(shù)值缺乏應(yīng)用和具體集聚位置不明等缺陷。第一，基于點(diǎn)對(duì)的地理距離構(gòu)建指數(shù)，理論上可以考慮企業(yè)質(zhì)量因素做加權(quán)處理，但是實(shí)際應(yīng)用中基本都是以未加權(quán)的地理距離進(jìn)行評(píng)估。因此，未能體現(xiàn)企業(yè)異質(zhì)性帶來(lái)的空間距離影響差異。隨著經(jīng)濟(jì)發(fā)展，交通和通信水平都在改善，經(jīng)濟(jì)聯(lián)系越來(lái)越緊密，單純的靜態(tài)地理距離衡量產(chǎn)業(yè)集聚程度卻是變化不明顯的?？臻g集聚的異質(zhì)性還體現(xiàn)在從許多中小企業(yè)的高度集聚到幾個(gè)大企業(yè)集聚，或者是一個(gè)大企業(yè)和許多中小企業(yè)聚集在一起，把這種現(xiàn)象稱(chēng)為企業(yè)集群和經(jīng)濟(jì)活動(dòng)集群[13]。第二，置信區(qū)間的隨機(jī)選擇位置缺乏科學(xué)依據(jù)，因?yàn)槊總€(gè)行業(yè)的選址具有不同的真實(shí)要求，不同行業(yè)位置并不是可以相互替代的，比如機(jī)械制造選址要求與化工產(chǎn)業(yè)就明顯不能替代。第三，除了EGA指數(shù)和Cluster指數(shù)方法，其他指標(biāo)由于其設(shè)定的評(píng)估方法非常理論，加權(quán)的計(jì)算比較復(fù)雜，并依賴(lài)函數(shù)檢驗(yàn)和圖形判斷集聚情況，而對(duì)于集聚程度的具體數(shù)值卻不重視，因此限制了將集聚程度作為一個(gè)重要的經(jīng)濟(jì)變量應(yīng)用于更為廣泛和深入的空間經(jīng)濟(jì)分析，特別是空間計(jì)量的可能性。第四，由于考察的結(jié)論是在某個(gè)尺度上是否集聚的情況，而對(duì)真實(shí)的地理空間的集聚位置并未指明，顯得模糊不清[15，17]。

此外，該類(lèi)函數(shù)的思想很簡(jiǎn)單，但主要是受限于計(jì)算量的約束?；诰嚯x的方法基本都是數(shù)據(jù)密集型和計(jì)算機(jī)依賴(lài)型。數(shù)據(jù)雖然容易得到了，但是算法過(guò)于細(xì)致導(dǎo)致速度很慢，系統(tǒng)資源占用很高，與基于集群的傳統(tǒng)指數(shù)相比，正是數(shù)據(jù)密集型和計(jì)算機(jī)依賴(lài)型限制了基于距離的方法的應(yīng)用[15]。所有這些函數(shù)的計(jì)算以及它們各自的置信區(qū)間估計(jì)非常耗時(shí)。計(jì)算時(shí)間與點(diǎn)的平方數(shù)(計(jì)算所有對(duì)點(diǎn)之間的距離)和零假設(shè)的虛擬次數(shù)成正比，內(nèi)存需求與點(diǎn)數(shù)的平方成正比，計(jì)算中涉及的點(diǎn)的數(shù)量取決于函數(shù)[18]。最近的發(fā)展已經(jīng)提出大數(shù)據(jù)集的處理的挑戰(zhàn)以及如何處理數(shù)據(jù)的重要計(jì)算需求[15]。

(三)產(chǎn)業(yè)集聚和點(diǎn)對(duì)距離融合的探索

單純根據(jù)地理位置距離設(shè)計(jì)的集聚程度指標(biāo)更多反映的是空間集聚的地理尺度，而我們感興趣的問(wèn)題是空間經(jīng)濟(jì)活動(dòng)的集聚程度及其地理位置。實(shí)際上，每個(gè)公司都有許多經(jīng)濟(jì)質(zhì)量特征(比如員工數(shù)量或資本)。利用企業(yè)選址概率模型對(duì)EG指數(shù)修正的研究是最早關(guān)于綜合集聚和位置的嘗試[19]。真正意義上綜合權(quán)衡集聚和距離的研究是建立一個(gè)集聚概率模型，混合使用基于聚類(lèi)的方法和基于距離的方法，識(shí)別產(chǎn)業(yè)集聚的空間尺度、空間位置和類(lèi)型的帶狀集聚理論[17]。上述關(guān)于產(chǎn)業(yè)集聚與地理位置融合的各類(lèi)模型擴(kuò)展具有極大的創(chuàng)新意義，但是囿于其模型設(shè)計(jì)的復(fù)雜，計(jì)算量的龐大，限制了其應(yīng)用推廣，基本沒(méi)有這些方法推廣應(yīng)用的研究文獻(xiàn)。

(四)總體述評(píng)

通過(guò)對(duì)上述兩條研究路線(xiàn)的梳理和比較發(fā)現(xiàn)，兩條研究路線(xiàn)存在較大的差異，各有優(yōu)劣。第一條研究路線(xiàn)，在按照主觀劃分地理單元的空間上，基于產(chǎn)值或者就業(yè)人數(shù)的評(píng)估，主要是提供一個(gè)數(shù)值。該數(shù)值更多體現(xiàn)的是專(zhuān)業(yè)化水平，而非集聚程度，且與產(chǎn)出等指標(biāo)存在明顯的共線(xiàn)性。由于缺乏閾值進(jìn)行標(biāo)準(zhǔn)的統(tǒng)計(jì)檢驗(yàn)，該數(shù)值難以做出是否集聚的判斷。第二條研究路線(xiàn)，基于微觀的地理數(shù)據(jù)，試圖對(duì)全域做出判斷，并對(duì)集聚做出嚴(yán)苛的統(tǒng)計(jì)檢驗(yàn)，反映的是空間集聚的地理尺度而非具體集聚位置。理論上可以綜合考慮就業(yè)人員或者產(chǎn)出的加權(quán)，但是實(shí)際應(yīng)用中很少做加權(quán)處理，原因在于計(jì)算的復(fù)雜度呈現(xiàn)幾何級(jí)別的增加。因此，大多是從整體分布角度描述集聚程度的分布，同時(shí)對(duì)于集中程度指數(shù)的不重視也限制了計(jì)量分析的應(yīng)用。因此，綜合考慮集群和距離等多種因素，融合企業(yè)多維信息的企業(yè)活動(dòng)空間集聚衡量的新方法，符合中國(guó)區(qū)域發(fā)展的差異性和競(jìng)爭(zhēng)性的實(shí)際，也是經(jīng)濟(jì)活動(dòng)空間集聚指標(biāo)研究的一個(gè)有價(jià)值的探索方向。

在回顧已有空間集聚衡量研究方法的基礎(chǔ)上，本文提出了基于引力模型的產(chǎn)業(yè)空間集聚衡量方法，定義為中心引力指數(shù)法，即基于企業(yè)微觀層面數(shù)據(jù)，綜合考慮了企業(yè)的地理位置和質(zhì)量指標(biāo)，構(gòu)造相對(duì)置信區(qū)間做統(tǒng)計(jì)檢驗(yàn)和設(shè)定絕對(duì)判斷閾值，并追蹤產(chǎn)業(yè)空間集聚的地理位置動(dòng)態(tài)。本文的主要?jiǎng)?chuàng)新之處在于：將企業(yè)的分布看作是多維的立體空間，而不是二維的平面地圖分布，考慮了多維分布的聚集問(wèn)題；試圖兼顧基于產(chǎn)值集聚程度衡量的簡(jiǎn)單易用和基于距離集聚程度衡量方法的統(tǒng)計(jì)檢驗(yàn)要求；通過(guò)產(chǎn)業(yè)視角，基于企業(yè)的微觀地理數(shù)據(jù)衡量集聚程度，能夠?qū)Ξa(chǎn)業(yè)的空間集聚程度、空間分布特征以及空間轉(zhuǎn)移過(guò)程做出合理的判斷和解釋?zhuān)_定產(chǎn)業(yè)集聚的空間地理位置。

三、中心引力指數(shù)算法設(shè)計(jì)

本文借鑒已有集聚指數(shù)的計(jì)算方法，包括基于產(chǎn)業(yè)和基于距離兩種研究路線(xiàn)，結(jié)合引力模型理論，提出衡量經(jīng)濟(jì)活動(dòng)空間集聚的中心引力指數(shù)算法。

(一)虛擬中心企業(yè)設(shè)定

中心企業(yè)的設(shè)定是中心引力指數(shù)核算的關(guān)鍵，也是確定產(chǎn)業(yè)空間集聚地理位置的關(guān)鍵。傳統(tǒng)的空間集聚程度的測(cè)度方法是使用產(chǎn)值或者人口規(guī)?？疾?，缺乏對(duì)空間因素尤其是距離的考察。為了修正該缺陷，同時(shí)得益于數(shù)據(jù)獲取分析技術(shù)的進(jìn)步，基于距離的方法異軍突起。關(guān)注地理空間考察，理論上可以兼顧產(chǎn)值和人口等質(zhì)量因素，但是現(xiàn)實(shí)應(yīng)用卻極其罕見(jiàn)，并且對(duì)集聚的真實(shí)地理空間位置語(yǔ)焉不詳。為彌補(bǔ)上述兩類(lèi)方法的缺陷，融合產(chǎn)業(yè)集聚和地理距離的研究做出了有益探索，但是受限于公式復(fù)雜和計(jì)算量的龐大，實(shí)際應(yīng)用較少。因此，綜合已有研究成果，基于引力模型，同時(shí)考察多個(gè)因素，特別是解決集聚的地理位置問(wèn)題，本文提出了中心引力指數(shù)方法。設(shè)定中心企業(yè)，為空間集聚的地理位置定位提供基礎(chǔ)，彌補(bǔ)了DO指數(shù)等無(wú)法定位的缺陷。具體設(shè)定方法如下：

第一步，劃定經(jīng)濟(jì)集聚核心企業(yè)集合?；诰嚯x的經(jīng)濟(jì)聚集測(cè)度的核心是經(jīng)濟(jì)活動(dòng)的集聚，其次才是空間分布問(wèn)題。由于全域企業(yè)的分布規(guī)模大小不一，位置較為分散，因此，經(jīng)濟(jì)集聚核心企業(yè)集合的選擇，更多的側(cè)重于經(jīng)濟(jì)活動(dòng)的集聚，在考察地理距離因素時(shí)需避免離群點(diǎn)極端值的影響。假設(shè)Z行業(yè)中有N家企業(yè)，企業(yè)的考察指標(biāo)涵蓋多個(gè)質(zhì)量指標(biāo)，統(tǒng)稱(chēng)為指標(biāo)集合Xi。按照經(jīng)濟(jì)質(zhì)量Xi指標(biāo)集合對(duì)行業(yè)所有企業(yè)分別排序，分別取每個(gè)指標(biāo)降序排列中位數(shù)以上的企業(yè)，共同組成一個(gè)企業(yè)集合，并剔除重復(fù)樣本，最后形成一個(gè)包含n家企業(yè)的經(jīng)濟(jì)集聚核心企業(yè)集合z。第二步，設(shè)定中心企業(yè)。此處的中心企業(yè)，不僅僅考慮地理位置的中心，還需要考慮經(jīng)濟(jì)活動(dòng)的中心。因此，取經(jīng)濟(jì)集聚核心企業(yè)集合z的企業(yè)各項(xiàng)指標(biāo)的中位數(shù)，作為一個(gè)虛擬中心企業(yè)的各項(xiàng)對(duì)應(yīng)指標(biāo)值，同時(shí)對(duì)各個(gè)企業(yè)的經(jīng)緯度坐標(biāo)取中位數(shù)，作為中心企業(yè)的地理位置坐標(biāo)，記為O點(diǎn)，用CO標(biāo)記虛擬中心企業(yè)。

根據(jù)上述設(shè)定規(guī)則，中心企業(yè)的各項(xiàng)指標(biāo)均值計(jì)算公式如下：

XO=median(Xi)

(1)

(二)引力模型介紹

1.兩點(diǎn)之間簡(jiǎn)單引力模型

設(shè)有兩個(gè)點(diǎn)，分別為點(diǎn)1和點(diǎn)2，各自的質(zhì)量(經(jīng)濟(jì)分析中可以用人口、產(chǎn)出和GDP等表示)分別為m1和m2，兩點(diǎn)之間的距離為d1，2(經(jīng)濟(jì)分析中可以是地理距離，也可以是其他具有空間性質(zhì)的變量)，則兩點(diǎn)之間的引力為g1，2，公式表示如下：

(2)

其中φ代表引力系數(shù)(1)φ在牛頓引力公式中是常數(shù)，但是在做經(jīng)濟(jì)分析時(shí)，在不同的分析情況下是變化的。。

2.兩點(diǎn)之間單因素引力模型一般形式

設(shè)有空間點(diǎn)的集合P，共有N個(gè)點(diǎn)，任意兩點(diǎn)i點(diǎn)和j點(diǎn)之間的引力為：

(3)

其中：α，β和γ是大于0的常數(shù)。

3.兩點(diǎn)之間多因素引力模型

在經(jīng)濟(jì)分析中，往往是由多點(diǎn)組成的空間，也同時(shí)考察多個(gè)因素共同作用，因此，兩點(diǎn)多因素引力模型變得更為復(fù)雜。

Gi，j=f(Xi，Xj，Si，j)

(4)

其中：X表示考察的因素的集合；S表示點(diǎn)之間的具有空間性質(zhì)的向量。假設(shè)考察k個(gè)因素，那么兩點(diǎn)之間引力合計(jì)為：

(5)

(三)中心引力指數(shù)估計(jì)及檢驗(yàn)

1.中心引力指數(shù)估計(jì)

對(duì)于有N個(gè)企業(yè)的Z行業(yè)，企業(yè)與虛擬中心企業(yè)形成N個(gè)引力點(diǎn)對(duì)關(guān)系。一個(gè)企業(yè)與虛擬中心企業(yè)之間形成多因素引力點(diǎn)對(duì)關(guān)系，中心點(diǎn)對(duì)距離記為dO，i。所有因素的中心引力加總，得到該企業(yè)的中心引力總值GO，i。

(6)

行業(yè)中所有企業(yè)的中心引力加總為GZ。

(7)

設(shè)定中心企業(yè)所在位置為圓心O點(diǎn)，距離r為半徑形成一個(gè)圓形區(qū)域，r∈[mindO，i，maxdO，i]，其中半徑r均勻的從圓心O點(diǎn)逐漸延長(zhǎng)到最遠(yuǎn)的企業(yè)距離maxdO，i。該圓形區(qū)域內(nèi)的企業(yè)組成經(jīng)濟(jì)集聚核心區(qū)企業(yè)集合，記為ZO，所有企業(yè)的中心引力值，記為GO。

(8)

行業(yè)Z的產(chǎn)業(yè)集聚程度，記為CRZ，用經(jīng)濟(jì)集聚核心區(qū)ZO中的企業(yè)中心引力總值GO與行業(yè)Z中所有企業(yè)的中心引力總值GZ的比值表示。

(9)

2.經(jīng)濟(jì)集聚檢驗(yàn)

產(chǎn)業(yè)是否集聚以及集聚的空間分布特征等，并不能簡(jiǎn)單的從前述行業(yè)空間集聚程度指標(biāo)的絕對(duì)值做出可靠的判斷。借鑒已有的空間集聚測(cè)度的檢驗(yàn)思路，需要設(shè)定一個(gè)科學(xué)的置信區(qū)間做統(tǒng)計(jì)檢驗(yàn)以及集聚判斷閾值，并作圖分析。基于對(duì)產(chǎn)業(yè)集聚和分散審慎的判斷，本文的檢驗(yàn)分兩步進(jìn)行：第一步，根據(jù)前述經(jīng)濟(jì)核心區(qū)面積和行業(yè)覆蓋區(qū)域面積指標(biāo)估計(jì)相對(duì)閾值指標(biāo)，構(gòu)造置信區(qū)間；第二步，提出空間集聚程度的絕對(duì)值標(biāo)準(zhǔn)作為判斷閾值。檢驗(yàn)的原則：只有依次符合上述兩步判斷，才能確認(rèn)集聚或者分散。

第一步，關(guān)于相對(duì)值判斷的置信區(qū)間設(shè)定。假設(shè)，核心區(qū)域面積記為AO；行業(yè)Z覆蓋的圓形區(qū)域面積，記為AZ；二者的比值，記為μO(píng)。此處的μO(píng)為集聚與分散的基準(zhǔn)值。

(10)

(11)

第二步，關(guān)于絕對(duì)值判斷的閾值設(shè)定。集聚與否，不等同于集聚程度，因此需要對(duì)不同的集聚程度做出更為細(xì)致的判斷，便于考察產(chǎn)業(yè)的集聚結(jié)構(gòu)特征。本文采用一定的比例將集聚程度分為低度集聚、中度集聚和高度集聚三個(gè)等級(jí)。本文選定的標(biāo)準(zhǔn)是將0.25和0.75兩個(gè)分位點(diǎn)作為三個(gè)等級(jí)劃分的臨界值。

(12)

四、中心引力指數(shù)算法應(yīng)用

(一)數(shù)據(jù)整理

中心引力指數(shù)算法應(yīng)用例證，使用的空間范圍為長(zhǎng)江三角洲地區(qū)三省一市，即江蘇省、浙江省、安徽省和上海市；研究的對(duì)象為制造業(yè)(兩位標(biāo)準(zhǔn)產(chǎn)業(yè)代碼SIC13-SIC43)，數(shù)據(jù)來(lái)源《中國(guó)工業(yè)企業(yè)數(shù)據(jù)庫(kù)(1998—2013)》中制造業(yè)企業(yè)微觀數(shù)據(jù)，其中2010年企業(yè)質(zhì)量指標(biāo)基本缺失，故剔除。根據(jù)長(zhǎng)三角地區(qū)的微觀企業(yè)數(shù)據(jù)的實(shí)際情況，選擇資產(chǎn)總計(jì)大于0、固定資產(chǎn)1 000萬(wàn)元以上、年銷(xiāo)售額500萬(wàn)元以上、就業(yè)人數(shù)30人以上，且企業(yè)地址相對(duì)完整的制造業(yè)企業(yè)為研究樣本。

數(shù)據(jù)的統(tǒng)計(jì)口徑和兩位數(shù)產(chǎn)業(yè)代碼在2002年和2012年發(fā)生兩次調(diào)整，因此整體數(shù)據(jù)存在匹配問(wèn)題。需要進(jìn)行行業(yè)大類(lèi)名稱(chēng)匹配，而不是簡(jiǎn)單的代碼匹配。數(shù)據(jù)匹配調(diào)整有三種基本法方案：一是前向調(diào)整，以1998年統(tǒng)計(jì)口徑為基準(zhǔn)，2003年之后數(shù)據(jù)做前向調(diào)整。二是后向調(diào)整，以2012年的統(tǒng)計(jì)口徑對(duì)之前的數(shù)據(jù)做相應(yīng)調(diào)整。三是以2002年的統(tǒng)計(jì)口徑為基準(zhǔn)做前后向靈活調(diào)整。本文以合并數(shù)據(jù)為首選調(diào)整原則，兼顧指標(biāo)的均衡處理，選擇2002年的統(tǒng)計(jì)口徑為基準(zhǔn)做數(shù)據(jù)匹配處理(2)具體數(shù)據(jù)匹配調(diào)整處理方法以及后續(xù)的數(shù)據(jù)驗(yàn)證、圖形分析，有需要的讀者，聯(lián)系作者索取。。煙草制品行業(yè)(SIC16)和武器彈藥(SIC39)屬于特殊行業(yè)，故排除在考察范圍之外。

研究指標(biāo)根據(jù)前述引力指數(shù)算法模型，選取企業(yè)地址、工業(yè)銷(xiāo)售產(chǎn)值、固定資產(chǎn)和資產(chǎn)總計(jì)四個(gè)指標(biāo)因素。企業(yè)的精準(zhǔn)經(jīng)緯度坐標(biāo)值是將企業(yè)詳細(xì)地址通過(guò)R語(yǔ)言的RCurl包對(duì)接高德地圖的API接口獲??；而虛擬中心企業(yè)的地理位置根據(jù)經(jīng)緯度坐標(biāo)，逆向處理。不同于現(xiàn)有大部分研究使用歐幾里得距離，本文應(yīng)用航海上廣泛運(yùn)用的蘭勃特公式計(jì)算大圓距離。本文的數(shù)據(jù)清洗、挖掘和計(jì)算是基于R語(yǔ)言3.5.3版本完成。

(二)計(jì)算步驟

1.設(shè)定中心企業(yè)

長(zhǎng)三角地區(qū)看作一個(gè)整體的經(jīng)濟(jì)區(qū)域，在兩位數(shù)代碼行業(yè)層面進(jìn)行分析。本文根據(jù)前述關(guān)于中心企業(yè)設(shè)定的指標(biāo)的計(jì)算方法，按照工業(yè)銷(xiāo)售產(chǎn)值(Vi)、資產(chǎn)總計(jì)(Ti)和固定資產(chǎn)(Fi)三個(gè)指標(biāo)對(duì)行業(yè)所有企業(yè)分別排序，分別取每個(gè)指標(biāo)降序排列大于中位數(shù)的企業(yè)，共同組成一個(gè)企業(yè)集合，并剔除重復(fù)樣本，最后形成一個(gè)包含n家企業(yè)的核心區(qū)域經(jīng)濟(jì)集聚企業(yè)集合z。中心企業(yè)的各項(xiàng)指標(biāo)由核心區(qū)域經(jīng)濟(jì)集聚企業(yè)集合z中企業(yè)的各項(xiàng)指標(biāo)的中位數(shù)得到，即根據(jù)式(1)計(jì)算。

2.計(jì)算中心點(diǎn)對(duì)的多因素引力

根據(jù)企業(yè)的經(jīng)緯度，應(yīng)用航海上大圓距離公式計(jì)算兩點(diǎn)之間的球面距離。假設(shè)已知兩點(diǎn)經(jīng)緯度P1(φ1，λ1)和P2(φ2，λ2)，應(yīng)用蘭勃特公式計(jì)算大地線(xiàn)長(zhǎng)。

(13)

其中，cosσ=sinφ1sinφ2+cosφ1cosφ2cosΔλ；a、c為橢球的長(zhǎng)半徑和扁率；σ為由地理經(jīng)緯度求得的球面距離。本文使用R語(yǔ)言Geosphere包中基于蘭勃特公式的Distm函數(shù)估計(jì)點(diǎn)對(duì)距離。

根據(jù)前述中心引力指數(shù)算法的設(shè)計(jì)方法核算行業(yè)的空間集聚程度。引力模型中參數(shù)的設(shè)定，參考白俊紅等關(guān)于引力模型參數(shù)設(shè)定的標(biāo)準(zhǔn)[20]，引力系數(shù)φ、引力因素參數(shù)α和β都設(shè)定為1，γ設(shè)定為2。此外，如果企業(yè)的地址與中心企業(yè)相同，或者非常接近，導(dǎo)致距離極小，將該點(diǎn)對(duì)地理距離設(shè)定為1千米，以避免出現(xiàn)中心引力極端值對(duì)估測(cè)的干擾。

3.空間集聚程度與檢驗(yàn)估計(jì)

根據(jù)前述算法設(shè)計(jì)思路，設(shè)定中心企業(yè)CO，將半徑的起始距離設(shè)定為最近的企業(yè)距離mindO，i，然后逐漸增加距離，實(shí)行逐點(diǎn)統(tǒng)計(jì)，連續(xù)估計(jì)經(jīng)濟(jì)空間集聚程度，并同步估計(jì)檢驗(yàn)標(biāo)準(zhǔn)。為了避免極值的影響和深入結(jié)構(gòu)化分析，本文將距離邊界設(shè)定為全域企業(yè)距離的中位數(shù)，選擇行業(yè)距離的中位數(shù)(mediandO，i)為分界點(diǎn)，將距離小于mediandO，i的企業(yè)組成考察企業(yè)集合，記為ZM，企業(yè)所在區(qū)域?yàn)榭疾靺^(qū)域，半徑記為R。已有的大多數(shù)基于距離的集聚指數(shù)的估計(jì)，對(duì)所有的行業(yè)設(shè)置固定考察距離，忽視了行業(yè)差異導(dǎo)致其分布的地理范圍不同的特征。因此，本文針對(duì)不同的行業(yè)以中位數(shù)為標(biāo)準(zhǔn)，也就是不同的行業(yè)選擇了不同的考察距離，符合產(chǎn)業(yè)自身的特征。集聚程度記為CRZ，與式(9)類(lèi)似，用中心區(qū)域的企業(yè)引力總值(GO)除以局域企業(yè)集合的引力總值(GM)得到。

(14)

檢驗(yàn)標(biāo)準(zhǔn)使用中心區(qū)域的面積(記為AO)與考察面積(記為AM)相比得到。

(15)

4.作圖檢驗(yàn)分析

最直觀的方式是按照前述閾值范圍作圖。按照式(11)和式(12)比較以判斷經(jīng)濟(jì)活動(dòng)集聚還是分散，以及集聚的程度分類(lèi)。繪圖可以直觀的判定產(chǎn)業(yè)集聚的空間尺度和集聚程度。繪圖過(guò)程中，將集聚程度使用原值帶入，而閾值分別用上限和下限帶入構(gòu)成置信區(qū)間。綜上，對(duì)于產(chǎn)業(yè)的集聚與否的判定標(biāo)準(zhǔn)具有絕對(duì)標(biāo)準(zhǔn)和相對(duì)標(biāo)準(zhǔn)兩個(gè)檢驗(yàn)，只有同時(shí)滿(mǎn)足兩個(gè)判定才能對(duì)集聚還是分散做出謹(jǐn)慎的斷定。

(三)典型行業(yè)示例

為了更好地理解中心引力指數(shù)法及其結(jié)果的評(píng)價(jià)和應(yīng)用，選取具有代表性特征的木材加工及竹、藤、棕、草制品業(yè)(SIC20)1998年、2002年、2007年和2013年的數(shù)據(jù)作為示例，從時(shí)間上分析空間分布的變化。之所以選擇該行業(yè)，因?yàn)樵撔袠I(yè)屬于傳統(tǒng)的制造業(yè)，具有勞動(dòng)密集和原料依賴(lài)特征，因此其隨著經(jīng)濟(jì)的發(fā)展會(huì)具有明顯的集聚和產(chǎn)業(yè)轉(zhuǎn)移特征。

1.產(chǎn)業(yè)全域分布與虛擬中心轉(zhuǎn)移

產(chǎn)業(yè)(SIC20)在研究時(shí)間跨度中的空間分布變化，同時(shí)表現(xiàn)出集聚程度逐漸提高和集聚的空間轉(zhuǎn)移兩個(gè)顯著特征(見(jiàn)圖1)。從圖1a到圖1d可以看出，隨著時(shí)間的推移，產(chǎn)業(yè)的集聚程度明顯提高，表現(xiàn)為圖中的黑點(diǎn)所代表的規(guī)模以上企業(yè)越來(lái)越密集。集聚的轉(zhuǎn)移的特征也非常顯著，從黑點(diǎn)所形成的簇可以發(fā)現(xiàn)，1998年的時(shí)候，產(chǎn)業(yè)主要集聚在環(huán)太湖地區(qū)，而到了2013年，產(chǎn)業(yè)主要集聚在江蘇北部的宿遷地區(qū)。

1a 1998年

1b 2002年

1c 2007年

1d 2013年圖1 長(zhǎng)三角產(chǎn)業(yè)(SIC20)全域分布與虛擬中心轉(zhuǎn)移圖

根據(jù)虛擬中心企業(yè)的位置轉(zhuǎn)移，考察產(chǎn)業(yè)集聚的空間轉(zhuǎn)移。產(chǎn)業(yè)集聚的虛擬中心企業(yè)作為衡量集聚的關(guān)鍵，其空間位置不僅是估計(jì)產(chǎn)業(yè)集聚程度的基礎(chǔ)，并且是產(chǎn)業(yè)轉(zhuǎn)移地理空間路徑的關(guān)鍵標(biāo)志。隨著時(shí)間變化，產(chǎn)業(yè)發(fā)生轉(zhuǎn)移，其集聚中心，也就是虛擬中心企業(yè)必然也會(huì)發(fā)生轉(zhuǎn)移。而中心企業(yè)位置的轉(zhuǎn)移正是代表了產(chǎn)業(yè)集聚在地理空間上的轉(zhuǎn)移。應(yīng)用前述數(shù)據(jù)處理技術(shù)，將所觀察產(chǎn)業(yè)的四個(gè)階段的虛擬中心企業(yè)地理經(jīng)緯度和具體地理位置求解。從圖1中虛擬中心企業(yè)的地址“?”標(biāo)記的變遷進(jìn)行分析發(fā)現(xiàn)，產(chǎn)業(yè)(SIC20)先后從上海到蘇州，再到無(wú)錫，最后到宿遷的產(chǎn)業(yè)轉(zhuǎn)移集聚路徑，表明該勞動(dòng)密集型兼具資源依賴(lài)型的傳統(tǒng)產(chǎn)業(yè)從長(zhǎng)三角的經(jīng)濟(jì)中心地帶向外圍轉(zhuǎn)移。虛擬中心企業(yè)位置轉(zhuǎn)移分析方法同樣適用于其他產(chǎn)業(yè)的轉(zhuǎn)移路徑的追蹤分析。

產(chǎn)業(yè)(SIC20)轉(zhuǎn)移的過(guò)程經(jīng)歷了單中心集聚，到雙中心集聚，再到單中心集聚的完整過(guò)程(見(jiàn)圖1)。從圖1a和圖1b可以看出，產(chǎn)業(yè)在蘇滬杭地區(qū)產(chǎn)生集聚；從圖1c可以看到，產(chǎn)業(yè)在向江蘇北部轉(zhuǎn)移，并逐漸形成新集聚，而原來(lái)的產(chǎn)業(yè)集聚中心依然存在，此時(shí)形成了雙中心集聚；從圖1d可以發(fā)現(xiàn)，最終該產(chǎn)業(yè)在新的地區(qū)形成集聚，回到了單中心集聚。這個(gè)集聚轉(zhuǎn)移的單中心到雙中心，再到單中心的過(guò)程，在中心引力指數(shù)法的應(yīng)用中也將得到驗(yàn)證(見(jiàn)圖2c)。

2.產(chǎn)業(yè)集聚程度與空間尺度變化

木材加工及竹、藤、棕、草制品業(yè)(SIC20)在四個(gè)年份中集聚的時(shí)間和空間變化基本涵蓋了中心引力指數(shù)法所要研究的內(nèi)容(見(jiàn)圖2)。在圖2中，點(diǎn)線(xiàn)組合的粗線(xiàn)是產(chǎn)業(yè)集聚程度；兩條虛線(xiàn)是根據(jù)面積比設(shè)定的置信區(qū)間，也就是相對(duì)判定標(biāo)準(zhǔn)，上面的是上限，下面的是下限；而0.25和0.75的水平實(shí)線(xiàn)是絕對(duì)判斷指標(biāo)臨界值，用以對(duì)集聚程度做絕對(duì)值區(qū)間劃分。遵從前述的判斷標(biāo)準(zhǔn)，從圖2可以看出：該產(chǎn)業(yè)在1998—2013年之間呈現(xiàn)逐漸集聚的特征，可以分為四個(gè)階段。第一階段，從圖2a可以看出，產(chǎn)業(yè)呈現(xiàn)低度集聚特征，空間分布尺度在20千米～50千米；第二階段，從圖2b可以看出，產(chǎn)業(yè)呈現(xiàn)出中度集聚狀態(tài)，集聚的空間尺度在50千米～70千米范圍，相較于上一階段，集聚的程度在提高的同時(shí)，集聚的空間尺度也在明顯擴(kuò)大；第三階段，從圖2c可以看出，自80千米左右的地方開(kāi)始呈現(xiàn)高度集聚，空間集聚的地理尺度也擴(kuò)大到了140千米左右；第四階段，從圖2d可以看出，集聚的空間范圍呈現(xiàn)收縮特征，空間集聚的地理尺度在20千米～90千米范圍呈現(xiàn)高度集聚。

2a 1998年

2b 2002年

2c 2007年

2d 2013年圖2 長(zhǎng)三角產(chǎn)業(yè)(SIC20)集聚程度與空間尺度圖

綜合上述四個(gè)階段，從圖2a到圖2c，我們可以看到該行業(yè)從低度集聚到高度集聚狀態(tài)是一個(gè)地理分布逐漸擴(kuò)大的過(guò)程；從圖2c到圖2d，可以發(fā)現(xiàn)在轉(zhuǎn)向高度集聚的時(shí)候，地理空間分布的地理尺度又呈現(xiàn)收縮，進(jìn)一步強(qiáng)化了集聚。具體體現(xiàn)在圖2c中的集聚曲線(xiàn)出現(xiàn)了明顯的兩個(gè)拐折，分別在50千米～70千米和80千米～140千米兩個(gè)空間尺度，表明出現(xiàn)了兩個(gè)集聚中心。這一點(diǎn)與前述的圖1c所變現(xiàn)的雙中心集聚現(xiàn)象相吻合，有力地證明了中心引力指數(shù)算法在偵測(cè)產(chǎn)業(yè)的多中心集聚方面是有效的。

(四)長(zhǎng)三角地區(qū)產(chǎn)業(yè)集聚的總體特征

通過(guò)對(duì)長(zhǎng)三角地區(qū)兩位代碼的制造業(yè)逐個(gè)分析，數(shù)據(jù)和圖形研究都表明該地區(qū)制造業(yè)表現(xiàn)出高度集聚的特征，并且絕大部分產(chǎn)業(yè)集聚都位于環(huán)太湖地區(qū)，只有少數(shù)的傳統(tǒng)產(chǎn)業(yè)向外圍轉(zhuǎn)移，見(jiàn)圖3。

3a 1998年

3b 2013年圖3 長(zhǎng)三角地區(qū)制造業(yè)1998年與2013年制造業(yè)全域分布對(duì)比圖

產(chǎn)業(yè)集聚主要向三個(gè)方向轉(zhuǎn)移：比如木材加工及竹、藤、棕、草制品業(yè)(SIC20)為代表的傳統(tǒng)勞動(dòng)密集型兼具資源依賴(lài)的產(chǎn)業(yè)在向長(zhǎng)三角北部地區(qū)轉(zhuǎn)移；食品加工業(yè)(SIC13)從蘇州依次轉(zhuǎn)向常州、鎮(zhèn)江和南京，即沿長(zhǎng)江經(jīng)濟(jì)帶向西轉(zhuǎn)移；紡織業(yè)(SIC17)和皮革制造業(yè)(SIC19)從蘇錫常地區(qū)向湖州和嘉興地區(qū)轉(zhuǎn)移。而其他大部分產(chǎn)業(yè)繼續(xù)在長(zhǎng)三角核心地區(qū)集聚強(qiáng)化。在研究期間內(nèi)，集聚的時(shí)間趨勢(shì)較為明顯，都是隨著時(shí)間的推移，集聚程度呈現(xiàn)上升的趨勢(shì)。即使發(fā)生了部分產(chǎn)業(yè)的轉(zhuǎn)移，以蘇滬杭為中心的長(zhǎng)三角核心地區(qū)依然是制造業(yè)中心。產(chǎn)業(yè)不同，集聚的空間尺度存在差異，總體上集聚的空間分布地理尺度在100千米范圍內(nèi)。

五、結(jié)論與展望

(一) 總體結(jié)論

1.中心引力指數(shù)算法是對(duì)已有產(chǎn)業(yè)空間集聚測(cè)度方法的綜合權(quán)衡和補(bǔ)充。現(xiàn)有空間集聚測(cè)度方法大體可以從基于產(chǎn)值和基于距離的兩個(gè)研究路線(xiàn)分為兩類(lèi)，前者忽略距離，后者所謂的產(chǎn)值加權(quán)缺乏實(shí)際應(yīng)用。究其原因，兩條路線(xiàn)的研究都是受制于技術(shù)，前者受限于數(shù)據(jù)采集技術(shù)，后者受限于計(jì)算能力限制。隨著技術(shù)的進(jìn)步，出現(xiàn)了一些基于實(shí)際距離的考察，但依然處于初步階段，未能突破運(yùn)算能力的限制，正是在這樣的背景下，本文提出了中心引力指數(shù)法，有效綜合了現(xiàn)有研究方法的優(yōu)點(diǎn)，規(guī)避了限制。通過(guò)融入成熟的引力模型，綜合了產(chǎn)值和距離，然后采用中心化方法，設(shè)定虛擬中心企業(yè)，有效擺脫了運(yùn)算能力的限制。因此，中心引力指數(shù)算法對(duì)已有的測(cè)度方法是一個(gè)可靠的拓展，形成綜合權(quán)衡和補(bǔ)充。

2.中心引力指數(shù)法能夠高效測(cè)度空間集聚。由于采用較為成熟的引力模型，融合產(chǎn)值和距離，并借鑒了已有測(cè)度方法的許多數(shù)學(xué)處理方法，通過(guò)設(shè)定虛擬中心企業(yè)的做法，在具備堅(jiān)實(shí)的理論基礎(chǔ)的同時(shí)，極大地簡(jiǎn)化了數(shù)學(xué)計(jì)算。特別是通過(guò)相對(duì)指標(biāo)構(gòu)造置信區(qū)間和絕對(duì)指標(biāo)判斷閾值，識(shí)別集聚程度并分類(lèi)，能夠有效地對(duì)經(jīng)濟(jì)活動(dòng)的空間集聚程度和地理空間分布特征做科學(xué)的結(jié)構(gòu)分析。方法的有效性還體現(xiàn)在使用長(zhǎng)三角地區(qū)的微觀企業(yè)數(shù)據(jù)進(jìn)行驗(yàn)證，表明了中心引力指數(shù)法能夠有效地測(cè)度集聚程度和空間分布特征，準(zhǔn)確偵測(cè)產(chǎn)業(yè)多中心集聚，并能夠根據(jù)虛擬中心企業(yè)的地理位置轉(zhuǎn)移，準(zhǔn)確追蹤產(chǎn)業(yè)集聚的轉(zhuǎn)移。

3.長(zhǎng)三角地區(qū)制造業(yè)集聚和空間分布特征分析。樣本數(shù)據(jù)分析表明，長(zhǎng)三角地區(qū)的制造業(yè)同時(shí)呈現(xiàn)集聚和轉(zhuǎn)移兩個(gè)特征。部分傳統(tǒng)產(chǎn)業(yè)，特別是勞動(dòng)密集型和資源依賴(lài)型產(chǎn)業(yè)在向長(zhǎng)三角的外圍轉(zhuǎn)移，虛擬中心企業(yè)的地理位置清楚表明了主要是向江蘇北部、長(zhǎng)江中上游和浙江中部三個(gè)方向轉(zhuǎn)移。特別強(qiáng)調(diào)的是，產(chǎn)業(yè)轉(zhuǎn)移而不是發(fā)散，轉(zhuǎn)移的過(guò)程中集聚程度在逐漸提升，最終在外圍的某個(gè)地區(qū)形成了高水平的集聚。比如江蘇北部的木材、南京的食品加工和嘉興的皮革制造等。蘇滬杭核心區(qū)的制造業(yè)的集聚程度在進(jìn)一步提高，也就是沒(méi)有發(fā)生轉(zhuǎn)移的產(chǎn)業(yè)的集聚程度進(jìn)一步強(qiáng)化。制造業(yè)的集聚基本都經(jīng)歷了集聚程度與地理空間分布同向變化，然后集聚程度提高，地理空間分布收縮的分階段過(guò)程。最終制造業(yè)都呈現(xiàn)高度集聚，同時(shí)在100千米范圍內(nèi)的空間分布特征。

(二)研究展望

中心引力指數(shù)算法能夠兼容現(xiàn)有方法的優(yōu)點(diǎn)，并規(guī)避其限制，應(yīng)用示例效果較好。但是，作為現(xiàn)有研究方法的權(quán)衡和融合，是對(duì)現(xiàn)有方法的有效補(bǔ)充，而非替代。中心引力指數(shù)算法尚存在一些不足，未來(lái)的研究可以從以下幾個(gè)方面拓展。

1.關(guān)于引力模型本身的拓展。本文所引入的引力模型是較為基本的模型，而且考察的企業(yè)的質(zhì)量指標(biāo)都是產(chǎn)值相關(guān)，而對(duì)于就業(yè)規(guī)模、技術(shù)水平和財(cái)稅貢獻(xiàn)等指標(biāo)并未考察。而這些指標(biāo)對(duì)于產(chǎn)業(yè)的集聚也非常關(guān)鍵，因此對(duì)于能夠融合更多因素，全面反映產(chǎn)業(yè)集聚影響因素的引力模型的拓展，將對(duì)中心引力指數(shù)算法的完善產(chǎn)生根本性的影響。

2.關(guān)于企業(yè)之間引力考察不足。中心引力指數(shù)算法，顧名思義，是圍繞著現(xiàn)實(shí)企業(yè)和虛擬中心企業(yè)之間的引力做出的產(chǎn)業(yè)集聚測(cè)度。而產(chǎn)業(yè)集聚是全體企業(yè)之間產(chǎn)生的綜合的關(guān)系，而不是單一的一對(duì)企業(yè)的關(guān)系，也就是所有點(diǎn)對(duì)之間的引力。這正是中心引力指數(shù)算法為了簡(jiǎn)化計(jì)算而產(chǎn)生的效率損失。關(guān)于點(diǎn)對(duì)企業(yè)的引力計(jì)算量將是極其龐大的，但是其衡量集聚程度更為精確。

3.關(guān)于多中心集聚的識(shí)別問(wèn)題。中心引力算法具有一個(gè)先天的缺陷，那就是假定了產(chǎn)業(yè)是單中心集聚。而事實(shí)上有不少產(chǎn)業(yè)是多中心集聚，尤其是產(chǎn)業(yè)在轉(zhuǎn)移的過(guò)程中，必然有一個(gè)多中心集聚的過(guò)程。如果能夠應(yīng)用諸如核密度或者泊松分布等統(tǒng)計(jì)檢驗(yàn)工具對(duì)多中心集聚進(jìn)行識(shí)別和偵測(cè)，結(jié)果會(huì)更為精準(zhǔn)。這應(yīng)該是完善中心引力指數(shù)算法的一個(gè)拓展方向。