施建剛， 段鍇豐， 吳光東

(1. 同濟大學 經濟與管理學院， 上海 200092; 2. 重慶大學 公共管理學院, 重慶 400044)

一、 引 言

破除城鄉(xiāng)二元結構壁壘，解決區(qū)域發(fā)展不平衡不充分問題，是緩解“人民日益增長的美好生活需要和不平衡不充分的發(fā)展之間的矛盾”的關鍵所在，也是實現(xiàn)共同富裕的必然要求。城鄉(xiāng)關系是最基本的社會經濟關系[1]，在區(qū)域發(fā)展過程中，必須將城市和鄉(xiāng)村作為一個有機整體，縮小城鄉(xiāng)發(fā)展差距，促進城鄉(xiāng)良性互動、協(xié)調共生，進而實現(xiàn)由城鄉(xiāng)分離走向城鄉(xiāng)融合之轉變。為此，中國政府相繼于2002年十六大報告中提出“城鄉(xiāng)統(tǒng)籌發(fā)展”、2012年十八大報告中提出“城鄉(xiāng)發(fā)展一體化”以及2017年十九大報告中提出“城鄉(xiāng)融合發(fā)展”等重大戰(zhàn)略性論述。從“城鄉(xiāng)統(tǒng)籌發(fā)展”到“城鄉(xiāng)發(fā)展一體化”，再到“城鄉(xiāng)融合發(fā)展”，這一過程本質上是“城鄉(xiāng)融合”思想在不同歷史階段的繼承與升華(1)“城鄉(xiāng)統(tǒng)籌發(fā)展”“城鄉(xiāng)發(fā)展一體化”和“城鄉(xiāng)融合發(fā)展”三者之間既有區(qū)別又有聯(lián)系。其區(qū)別在于，城鄉(xiāng)統(tǒng)籌發(fā)展強調城鄉(xiāng)雙方的個體地位，重視政府在統(tǒng)籌城鄉(xiāng)協(xié)調發(fā)展過程中的主體作用，側重點在于推動城鄉(xiāng)關系演進的方式及手段；城鄉(xiāng)發(fā)展一體化弱化了城鄉(xiāng)雙方的個體地位，強調縮小城鄉(xiāng)差距，實現(xiàn)城鄉(xiāng)一體化發(fā)展，側重點在于城鄉(xiāng)關系演進的最終目標；城鄉(xiāng)融合發(fā)展則在承認城鄉(xiāng)雙方個體地位的基礎上，強調城鄉(xiāng)之間的雙向互動融合和城鄉(xiāng)發(fā)展的平等性、可持續(xù)性，側重點在于城鄉(xiāng)關系演進的狀態(tài)和過程。其聯(lián)系在于，三者的理念精髓和基本目標是一致的，均強調城鄉(xiāng)發(fā)展的整體性和不可分割性，均以實現(xiàn)城鄉(xiāng)良性互動、協(xié)調發(fā)展和深度融合為目標。(魏后凱：《深刻把握城鄉(xiāng)融合發(fā)展的本質內涵》，《中國農村經濟》，2020年第6期，第5-8頁)，這種政策導向的演進反映了我國對形塑新型城鄉(xiāng)關系的不懈努力和逐步深化，順應了不同歷史階段城鄉(xiāng)關系演進的時代特征及發(fā)展趨勢[2-3]。

在中國城鄉(xiāng)關系由二元結構走向融合發(fā)展的演進過程中，如何科學定量地衡量城鄉(xiāng)關系的演進狀態(tài)成為學術界需要重點關注的課題，而城鄉(xiāng)融合發(fā)展水平的科學測度是城鄉(xiāng)關系問題由定性分析轉向定量評估的基礎。現(xiàn)有研究主要通過構建綜合評價指標體系對城鄉(xiāng)關系演進狀態(tài)進行定量測度。例如，周江燕等[4]、周佳寧等[5]、趙德起等[6]對全國省域層面的城鄉(xiāng)一體化、融合狀態(tài)進行了定量評價，從全國整體視角探索了城鄉(xiāng)關系的演進狀態(tài)。由于我國幅員遼闊，不同地區(qū)自然條件和人文地理環(huán)境差異明顯，故城鄉(xiāng)關系研究亦有必要落實到區(qū)域層面，因地制宜地探索區(qū)域城鄉(xiāng)融合發(fā)展的現(xiàn)實狀態(tài)，方能有效制定促進不同區(qū)域城鄉(xiāng)融合的方針政策。為此，在區(qū)域層面，漆莉莉[7]、劉融融等[8]、張新林等[9]、張海朋等[10]分別對我國中部地區(qū)、西北地區(qū)、淮海經濟區(qū)和環(huán)首都地區(qū)的城鄉(xiāng)融合、協(xié)調發(fā)展程度進行了定量測度。

盡管現(xiàn)有研究對城鄉(xiāng)關系演進狀態(tài)開展了較為豐富的定量探索，但其研究尺度多集中在國家或省域層面，對城市尺度的分析較為缺乏，且缺少對區(qū)域內部城鄉(xiāng)融合發(fā)展水平時空特征的考察。此外，鮮有研究對長江三角洲(以下簡稱“長三角”)地區(qū)的城鄉(xiāng)融合發(fā)展水平進行定量測度。長三角地區(qū)是我國綜合實力最強的區(qū)域之一，在我國社會主義現(xiàn)代化建設全局中具有重要的戰(zhàn)略地位和突出的帶動作用。該地區(qū)以上海，江蘇省的南京、蘇州、無錫、常州、鹽城、泰州、揚州、鎮(zhèn)江、南通，浙江省的杭州、溫州、寧波、湖州、金華、嘉興、紹興、臺州、舟山，安徽省的合肥、馬鞍山、蕪湖、銅陵、滁州、安慶、池州、宣城27座城市為中心區(qū)，輻射帶動整個地區(qū)高質量發(fā)展。城鄉(xiāng)融合發(fā)展作為解決長三角地區(qū)城鄉(xiāng)矛盾的重要手段，對促進該地區(qū)區(qū)域一體化發(fā)展具有重要意義。鑒于此，本文以長三角地區(qū)為研究對象，在科學把握城鄉(xiāng)融合發(fā)展時代內涵的基礎上，構建了一套城鄉(xiāng)融合發(fā)展水平測度指標體系，并基于2003—2018年該地區(qū)27座中心區(qū)城市的面板數據對城鄉(xiāng)融合發(fā)展水平進行評價，進而對其城鄉(xiāng)融合發(fā)展水平的時空特征進行了分析。本文的研究結果不僅可以豐富現(xiàn)有城鄉(xiāng)融合理論體系，亦可為制定促進長三角地區(qū)城鄉(xiāng)融合發(fā)展的方針政策提供科學依據。

二、 城鄉(xiāng)融合發(fā)展的內涵闡釋

當前，城鄉(xiāng)融合發(fā)展已經明確成為主導中國未來城鄉(xiāng)關系演進的核心方向[11]。伴隨著新時代工業(yè)化、信息化、城鎮(zhèn)化和農業(yè)現(xiàn)代化的快速推進，城鄉(xiāng)經濟社會聯(lián)系日益緊密，空間關聯(lián)不斷增強，城鄉(xiāng)融合發(fā)展被不斷賦予新的時代內涵。不同學科領域通常從不同角度對城鄉(xiāng)融合發(fā)展的本質進行解讀。以劉易斯二元經濟論為代表的新古典經濟學理論從城鄉(xiāng)資源配置的角度探討城鄉(xiāng)融合發(fā)展的內涵，認為城鄉(xiāng)融合應追求城鄉(xiāng)資源要素的優(yōu)化組合以實現(xiàn)生產效益最優(yōu)[12-13]。社會學理論強調城鄉(xiāng)區(qū)域以及城鄉(xiāng)參與主體之間的公平性，認為城鄉(xiāng)融合發(fā)展意味著城鄉(xiāng)居民享受平等的社會福利、均等的公共服務以及同質的生活條件[14]。生態(tài)學理論認為，城鄉(xiāng)融合發(fā)展是指城鄉(xiāng)區(qū)域生態(tài)群落合理分布，城鄉(xiāng)經濟社會發(fā)展與生態(tài)環(huán)境和諧共生，共同實現(xiàn)可持續(xù)發(fā)展[15-16]。經濟地理學在傳統(tǒng)經濟學的基礎上繼續(xù)拓展，強調城鄉(xiāng)融合發(fā)展必須注重城鄉(xiāng)區(qū)域的空間關聯(lián)[17-18]，認為城鄉(xiāng)融合是以一體化基礎設施體系和交通體系為支撐，城鄉(xiāng)空間實現(xiàn)動態(tài)均衡的過程[19-20]。系統(tǒng)論則將城鄉(xiāng)融合發(fā)展視為城鄉(xiāng)相互作用形成的多層次、多要素的復雜關聯(lián)系統(tǒng)，主張城鄉(xiāng)融合應是經濟、人口(人類)、社會、區(qū)域(空間)、環(huán)境(自然)等多維系統(tǒng)的協(xié)調有序發(fā)展[1][21]。

基于上述理論背景，多數學者認為，城鄉(xiāng)要素的雙向自由流動是城鄉(xiāng)融合發(fā)展的關鍵，城鄉(xiāng)融合發(fā)展本質上是在城鄉(xiāng)發(fā)展要素自由流動、公平和共享的基礎上實現(xiàn)城鄉(xiāng)協(xié)調、一體化發(fā)展[1][6][11][20][22-23]。例如，高帆[23]從四個方面描繪了城鄉(xiāng)融合發(fā)展的時代內涵，其認為城鄉(xiāng)融合發(fā)展意味著：(1)城鄉(xiāng)要素流動性與再配置功能增強；(2)城鄉(xiāng)產業(yè)多樣化特征增強；(3)城鄉(xiāng)空間交錯化程度攀升；(4)城鄉(xiāng)居民社會福利均等化程度提高。根據流動要素的類型，部分學者認為城鄉(xiāng)融合發(fā)展的核心在于實現(xiàn)“人”“地”以及“資本”的融合[8][24]。此外，學者們普遍認為，城鄉(xiāng)融合發(fā)展應強調從經濟一維向社會(生活)、經濟、空間、生態(tài)環(huán)境等多維度拓展，并重視不同維度之間的平衡發(fā)展[9][16]。例如，何仁偉[20]應用城鄉(xiāng)空間均衡模型和等值化曲線刻畫了城鄉(xiāng)融合發(fā)展的內涵與機理，指出城鄉(xiāng)融合發(fā)展是基于空間結構優(yōu)化及制度供給創(chuàng)新的社會、經濟、環(huán)境全面融合的發(fā)展，其最終目標是為了實現(xiàn)城鄉(xiāng)多維均衡發(fā)展、居民生活品質相當。楊志恒[25]認為，城鄉(xiāng)融合發(fā)展是以城鄉(xiāng)資源要素流動為主線，通過協(xié)調制度、產業(yè)、設施和生態(tài)環(huán)境等方面以實現(xiàn)經濟協(xié)作聯(lián)合、空間平等有序、生態(tài)本底約束以及社會相互認同的發(fā)展。周佳寧等[5][26-27]將“流空間”和“城鄉(xiāng)等值化”的理念置于中國城鄉(xiāng)關系發(fā)展背景下，提出城鄉(xiāng)融合發(fā)展應該體現(xiàn)為經濟、社會、人口、空間和生態(tài)環(huán)境的多維空間重構。

基于上述對城鄉(xiāng)融合發(fā)展的內涵闡釋，本文認為，城鄉(xiāng)融合發(fā)展的本質內涵在于，以尊重城鄉(xiāng)功能特征及分工差異為前提，以城鄉(xiāng)要素自由流動和資源公平共享為基礎，通過經濟自發(fā)力量與政策制度力量相結合，破除城鄉(xiāng)二者在經濟、社會、空間等維度的二元對立關系，實現(xiàn)城鄉(xiāng)經濟、人口、社會、空間及生態(tài)多維度全面協(xié)調共生、良性互動、深度融合和共同繁榮。

三、 評價指標體系構建和數據獲取

(一) 評價指標體系構建

城鄉(xiāng)融合發(fā)展既是城鄉(xiāng)關系演進的目標，也是一種狀態(tài)和過程。因此，城鄉(xiāng)融合發(fā)展水平測度指標體系中應同時包含反映“城鄉(xiāng)差距縮小”的對比類指標以及反映“城鄉(xiāng)共同繁榮”的狀態(tài)類指標。為了科學測度長三角地區(qū)各中心區(qū)城市2003—2018年的城鄉(xiāng)融合發(fā)展水平，基于上述城鄉(xiāng)融合發(fā)展的本質內涵，參考當前涉及城鄉(xiāng)融合發(fā)展水平評價的現(xiàn)有文獻[5][17]，同時遵循系統(tǒng)性、科學性、全面性和可獲得性等原則，本文從“經濟融合”“人口融合”“社會融合”“空間融合”和“生態(tài)融合”五個維度設計構建了城鄉(xiāng)融合發(fā)展水平測度指標體系，如表1所示。

表1 城鄉(xiāng)融合發(fā)展水平測度指標體系

(二) 數據獲取

城鄉(xiāng)融合發(fā)展水平測度所使用的相關原始數據主要來源于長三角地區(qū)各中心區(qū)城市歷年統(tǒng)計年鑒、《中國城市統(tǒng)計年鑒》、政府網站以及統(tǒng)計公報等；對于部分缺失的數據采用線性插值法的方式進行補全。

四、 研究方法

(一) 加速遺傳算法投影尋蹤模型(RAGA-PPC)

考慮到城鄉(xiāng)融合發(fā)展水平評價指標體系的多維特征，本文采用加速遺傳算法投影尋蹤模型(RAGA-PPC)對長三角地區(qū)城鄉(xiāng)融合發(fā)展水平進行科學評價。該方法將投影尋蹤模型(Projection Pursuit Classification, PPC)和基于實數編碼的加速遺傳算法(Real Coded Accelerating Genetic Algorithm, RAGA)相結合，綜合了PPC和RAGA方法各自的優(yōu)勢[28]。同傳統(tǒng)的主、客觀賦值方法相比，該方法可有效避免主觀因素的影響，同時得出的評價結果也相對客觀，因而更具適用性、穩(wěn)健性和準確性。應用RAGA-PPC方法測度長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的主要步驟如下：

步驟1評價樣本指標集的歸一化處理。以長三角地區(qū)各中心區(qū)城市作為研究對象，構建各指標的初始樣本集合：{X*(i，j)|i=1，2，3，…，n；j=1，2，3，…，p}。其中X*(i，j)為第i個樣本城市的第j個評價指標值，n，p分別為樣本城市的個數和評價指標的數目。為消除量綱對整體評價結果的影響，對采集的數據進行歸一化處理：

若該指標屬性為“正”：

(1)

若該指標屬性為“負”：

(2)

其中，Xmin(j)和Xmax(j)分別為第j個評價指標的最小值和最大值，X(i，j)為評價指標值的歸一化處理結果。

步驟2投影指標函數的構建。PPC可以將p維數據{X(i，j)|i=1，2，3，…，n；j=1，2，3，…，p}綜合成以a={a(1)，a(2)，…，a(p)}為最佳投影方向的一維投影值Z(i)。

(3)

其中，a為單位長度向量。投影指標函數Q(a)的表達式為：

Q(a)=Sz·Dz

(4)

其中，Sz和Dz分別為投影值Z(i)的標準差和局部密度，計算公式如下：

(5)

(6)

其中，E(z)為投影值Z(i)的平均值，R為局部密度的窗口半徑。R的取值要確保窗口內投影點的個數不能過少，避免滑動平均偏差太大，可根據多次試驗來確定。r(i,j)代表樣本之間的距離，r(i,j)=|Z(i)-Z(j)|。u(R-r(i,j))為單位躍階函數。當R≤r(i,j)時，該函數值為0；當R≥r(i,j)時，該函數值為1。

步驟3投影指標函數的優(yōu)化。對于確定的樣本指標集，投影方向a的變化會導致投影指標函數Q(a)發(fā)生變化，投影指標函數取最大值時，對應的向量反映了最佳投影方向，即最大概率暴露多維度數據結構特征的投影方向。最佳投影方向可以通過應用RAGA方法求解高維全局非線性最優(yōu)問題進行估計。在本文中，選定父代初始種群規(guī)模N= 400，交叉概率ρ= 0.80，變異概率ρm=0.80，優(yōu)秀個體數目選定為n=15個，a=0.05，加速次數為30。經運算得出最佳投影值a*，然后將其代入式(4)可以得到最大化目標函數，即：

MaxQ(a)=Sz·Dz

(7)

(8)

步驟4綜合評價分析。將步驟3求得的最佳投影向量a*代入式(3)，可得到各樣本城市的投影值Z*(i)，即長三角地區(qū)27座中心區(qū)城市的城鄉(xiāng)融合發(fā)展水平值。

(二) 核密度估計(KDE)

為探究長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的動態(tài)變化趨勢，本文應用核密度估計法(Kernel Density Estimation, KDE)對其時序動態(tài)演進特征進行分析。核密度估計法可以通過探究函數曲線收斂范圍以及收斂度的變化，對城鄉(xiāng)融合發(fā)展水平總體差異在不同階段的動態(tài)演進狀況進行可視化展示。其優(yōu)勢在于：(1)可通過連續(xù)的密度曲線對隨機變量的數據分布形態(tài)以及時序動態(tài)演進狀況進行直觀呈現(xiàn)，而不需要進行任何形式的模型參數假設。(2)可直接對隨機變量本身的數據特征進行分析，避免了參數模型估計形式的主觀影響。結合分布形式以及核密度函數圖的位置、形狀和峰度特征，可以對長三角地區(qū)城鄉(xiāng)融合發(fā)展水平在不同時期的變化情況進行判斷，從而探究其時序動態(tài)演進特征。

(三) 探索性空間分析(ESDA)

探索性空間分析(Exploratory Spatial Data Analysis, ESDA)是一種揭示變量之間的空間異質性、空間依賴性以及空間自相關性的方法[29]。地理學第一定律認為，事物之間存在一定程度上的相互關聯(lián)。隨著區(qū)域城鄉(xiāng)之間聯(lián)系日益緊密，要素流動不斷加強，同一區(qū)域內的城鄉(xiāng)發(fā)展必定存在緊密關聯(lián)[5]，故此本文采用ESDA方法中主流的Moran’s I指數來度量長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的空間相關性。Moran’s I指數可分為全局Moran’s I指數和局域Moran’s I指數，分別可用來探究樣本整體和局域的空間集聚或離散狀況。二者的計算公式分別為：

(9)

(10)

(四)趨勢面分析

趨勢面分析是一種利用數學曲面模擬呈現(xiàn)地理系統(tǒng)要素在空間上的分布規(guī)律和變化趨勢的方法。其本質上是通過運用回歸分析的原理，以最小二乘法擬合一個二維非線性函數，進而直觀地呈現(xiàn)地理要素在空間上的分布規(guī)律和變化趨勢。趨勢面分析方法可以幫助揭示長三角地區(qū)各中心區(qū)城市城鄉(xiāng)融合發(fā)展水平的空間分異規(guī)律。

五、 結果分析

(一) 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的總體變化特征分析

運用加速遺傳算法投影尋蹤模型，得出長三角地區(qū)27座中心區(qū)城市城鄉(xiāng)融合發(fā)展水平均值，如圖1所示?？梢园l(fā)現(xiàn)，2003—2018年16年間長三角地區(qū)城鄉(xiāng)融合發(fā)展水平總體上呈現(xiàn)波動上升趨勢。其中，2003—2008年為城鄉(xiāng)融合發(fā)展水平快速上升階段，而2008—2018年為城鄉(xiāng)融合發(fā)展水平波動階段。總體來看，長三角地區(qū)城鄉(xiāng)融合發(fā)展水平上升幅度不大，年增長率僅為0.49%。

圖1 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平變化(2003—2018)

運用核密度估計法，本文分別對2003年、2008年、2013年和2018年長三角地區(qū)城鄉(xiāng)融合發(fā)展水平進行了核密度估計，以此從全域視角刻畫長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的時序演進特征(見圖2)。

圖2 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的

首先，從核密度年度曲線的重心位置來看，2003—2008年重心位置向右遷移，2008—2013年向左遷移，2013—2018年則基本保持不變。這表明在研究期內長三角地區(qū)城鄉(xiāng)融合發(fā)展水平整體呈現(xiàn)出先上升后下降繼而水平波動的演進特征。

其次，從核密度曲線主峰的波峰高度來看，2003—2008年有所上升，2008—2013年略微下降，2013—2018年則出現(xiàn)較大幅度的回升。這說明長三角地區(qū)各城市間城鄉(xiāng)融合發(fā)展水平差異呈現(xiàn)先縮小再擴大繼而再縮小的變化趨勢。

再次，從核密度曲線的波峰數量來看，2008年和2018年出現(xiàn)主峰和次峰雙峰共存的情況，說明這兩年長三角地區(qū)各城市的城鄉(xiāng)融合發(fā)展水平呈現(xiàn)一定程度的兩極分化格局。2003年和2013年的波峰個數僅為1個，說明未出現(xiàn)兩極分化。

最后，從核密度曲線的左右拖尾來看，各年對應曲線的左側拖尾大于右側拖尾，呈現(xiàn)左偏分布的特征，且2003—2018年左側拖尾呈加長趨勢。這表明左側城市的數量較右側更多，即研究期內城鄉(xiāng)融合發(fā)展水平值小于平均值的城市個數更多，且低值區(qū)城市占比有所增加，城鄉(xiāng)融合發(fā)展水平低值集聚分布明顯。

由核密度分析可以發(fā)現(xiàn)，總體來看，長三角地區(qū)城鄉(xiāng)融合發(fā)展在整體發(fā)展水平、城市間差異、極化程度等方面均體現(xiàn)了隨時間波動變化的動態(tài)演進特征。

(二) 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的城市層面分析

根據長三角地區(qū)27座中心區(qū)城市的城鄉(xiāng)融合發(fā)展水平測度結果，繪制各城市2003—2018年城鄉(xiāng)融合發(fā)展水平箱形圖，如圖3所示。

圖3 長三角地區(qū)各中心區(qū)城市城鄉(xiāng)融合發(fā)展水平箱形圖

箱形圖的優(yōu)勢在于它不僅可以反映數據的離散分布情況，還可以進行多組數據分布特征的比較。圖3中每座城市對應的箱形圖最上部和最下部兩條水平線分別代表箱形圖的上邊緣和下邊緣；箱形的上下兩邊分別代表上四分位點和下四分位點；箱形內部的水平橫線代表的是城鄉(xiāng)融合發(fā)展水平的中位數，可以用于比較各城市城鄉(xiāng)融合發(fā)展程度的平均水平；箱形圖附近的圓點則代表異常值。

由圖3可以看出，安慶、滁州、合肥、上海、溫州、舟山等城市對應的箱形較長，說明這些城市的城鄉(xiāng)融合發(fā)展水平值在2003—2018年期間波動較大；相比之下，常州、金華、南京、蘇州、揚州等城市的城鄉(xiāng)融合發(fā)展水平16年內的變化幅度不大，保持了穩(wěn)健的發(fā)展態(tài)勢。對比各城市的城鄉(xiāng)融合發(fā)展水平中位數可以發(fā)現(xiàn)：蘇州、無錫、上海三座城市的城鄉(xiāng)融合發(fā)展水平明顯高于其他城市；其余城鄉(xiāng)融合發(fā)展水平相對領先的城市有常州、杭州、湖州、嘉興、南京、寧波、揚州、鎮(zhèn)江、舟山等；而安慶、池州和滁州則為城鄉(xiāng)融合發(fā)展水平排名最后三位的城市，有必要因地制宜進行政策引導，促進其城鄉(xiāng)融合發(fā)展。從省級層面來看，江蘇和浙江兩省的城鄉(xiāng)融合發(fā)展水平高于安徽省，這與區(qū)域城鄉(xiāng)經濟社會發(fā)展及城鎮(zhèn)化水平有一定關系。

(三) 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的時空格局演變

1. 全局空間格局演變特征

利用全局Moran’s I指數計算公式，借助GeoDa軟件對2003—2018年長三角地區(qū)各年城鄉(xiāng)融合發(fā)展水平的全局Moran’s I指數進行測算，測度結果如圖4所示。

圖4 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平全局Moran’s I指數演變特征

結果顯示：2003—2018年長三角地區(qū)各年城鄉(xiāng)融合發(fā)展水平的全局Moran’s I指數均為正值，在0.209—0.618之間，且各年的P-Value均小于0.05，即各年均通過了顯著性檢驗。這一結果表明，研究期內長三角地區(qū)城鄉(xiāng)融合發(fā)展水平空間分布存在顯著的正向全局空間自相關性，即本地城市的城鄉(xiāng)融合發(fā)展水平不僅會對相鄰城市產生影響，也會受到相鄰城市的影響?？傮w來看，2003—2018年全局Moran’s I指數呈現(xiàn)出先上升后下降再上升的變化趨勢，說明在整個研究期內長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的全局空間自相關性呈現(xiàn)出“弱相關→強相關→弱相關→強相關”的螺旋演進態(tài)勢。

2. 局域空間格局演變特征

全局Moran’s I指數雖然可以從宏觀上驗證長三角地區(qū)城鄉(xiāng)融合發(fā)展水平存在正向全局空間自相關性，但是無法對該地區(qū)的局域空間格局特征進行具象化展示。故此，有必要利用局域Moran’s I指數來刻畫長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的局域空間格局特征。為了考察該地區(qū)城鄉(xiāng)融合發(fā)展水平的局域空間分異變化，本文繪制了2003年、2008年、2013年和2018年城鄉(xiāng)融合發(fā)展水平的局域Moran’s I散點圖，如圖5所示。

圖5通過四個象限將局域空間格局形態(tài)劃分為四類：(1)第一類為高—高集聚型(HH，位于第一象限)。該類型代表本區(qū)域城鄉(xiāng)融合發(fā)展水平較高，其周圍相鄰城市的城鄉(xiāng)融合發(fā)展水平也相對較高，表現(xiàn)為“中心高，四周高”的高水平空間均衡關聯(lián)集聚狀態(tài)。(2)第二類為低—高集聚型(LH，位于第二象限)。該類型代表本區(qū)域城鄉(xiāng)融合發(fā)展水平較低，但其周圍相鄰城市的城鄉(xiāng)融合發(fā)展水平相對較高，表現(xiàn)為“中心低，四周高”的空間非均衡關聯(lián)集聚狀態(tài)。(3)第三類為低—低集聚型(LL，位于第三象限)。該類型代表本區(qū)域城鄉(xiāng)融合發(fā)展水平較低，其周圍相鄰城市的城鄉(xiāng)融合發(fā)展水平也相對較低，表現(xiàn)為“中心低，四周低”的低水平空間均衡關聯(lián)集聚狀態(tài)。(4)第四類為高—低集聚型(HL，位于第四象限)。該類型代表本區(qū)域城鄉(xiāng)融合發(fā)展水平相對較高，但是其周圍相鄰城市的城鄉(xiāng)融合發(fā)展水平相對較低，表現(xiàn)為“中心高，四周低”的空間非均衡關聯(lián)集聚狀態(tài)。參考已有類似研究成果[31]，將HH型地區(qū)定義為擴散互溢區(qū)，LH型地區(qū)定義為極化效應區(qū)，LL型地區(qū)定義為低速增長區(qū)，HL型地區(qū)定義為落后過渡區(qū)。由圖5可以看出：長三角地區(qū)城鄉(xiāng)融合發(fā)展在各階段均存在一定程度上的局部空間集聚現(xiàn)象，擴散互溢區(qū)和低速增長區(qū)兩種空間分異特征占據該地區(qū)城鄉(xiāng)融合發(fā)展演變的主導；城鄉(xiāng)融合發(fā)展的整體水平形成了以上海、蘇州、無錫、嘉興等城市為代表的高值簇(HH)和以池州、安慶、蕪湖、馬鞍山等城市為代表的低值簇(LL)為主的局部空間集聚特征，在空間上大致呈現(xiàn)為以上海和安慶為端點，由東向西逐級遞減的階梯狀分布格局。

從研究期內各年四種類型城市數量的變化情況來看(圖6)，2003—2018年長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的局域空間格局總體上呈現(xiàn)“HH>LL>HL>LH”的數量特征。由局域空間格局演變特征不難發(fā)現(xiàn)，長三角地區(qū)城鄉(xiāng)融合發(fā)展水平存在較大程度上的空間鄰近同伴效應(HH和LL類型城市數量明顯高于HL和LH類型城市數量)，即當本地區(qū)周圍相鄰城市的城鄉(xiāng)融合發(fā)展水平較高(低)時，本地區(qū)也更容易成為較高(低)城鄉(xiāng)融合發(fā)展水平的城市，由此則形成了“弱者恒弱、強者恒強”的局域空間集聚格局。該空間格局形成的主要原因在于：長三角地區(qū)相鄰城市在經濟、社會、文化、資源稟賦等方面相似度較高，空間距離較近，交通成本較低，信息獲取方便快捷，城市間相互聯(lián)系頻繁，因而在城鄉(xiāng)融合發(fā)展相關政策措施的執(zhí)行方面存在較大程度上的關聯(lián)性，故而形成了地理上的空間溢出效應?？傮w來看，HH類型城市數量在波動中略微上升；LL類型城市數量則呈較明顯的波動下降趨勢；HL類型和LH類型城市基本維持在低數量波動態(tài)勢中，該結果反映了研究期內長三角地區(qū)城鄉(xiāng)融合發(fā)展空間格局中低水平集聚狀態(tài)的城市正在朝高水平集聚狀態(tài)的方向緩慢發(fā)展。

圖6 2003—2018年長三角地區(qū)城鄉(xiāng)融合發(fā)展水平局域空間格局演變特征

3. 長三角地區(qū)各中心區(qū)城市城鄉(xiāng)融合發(fā)展水平的時空格局演變

選取研究期內2003年、2008年、2013年和2018年作為四個關鍵時間節(jié)點年份，運用ArcGIS軟件中自然間斷點法將長三角地區(qū)各中心區(qū)城市城鄉(xiāng)融合發(fā)展水平劃分為四個等級梯度，繪制各年的空間分布示意圖，如圖7所示。

圖斑顏色由深到淺分別從高到低對應城鄉(xiāng)融合發(fā)展水平的四個梯度。由圖7可以看出：長三角地區(qū)城鄉(xiāng)融合發(fā)展呈現(xiàn)出較明顯的空間集聚效應，城鄉(xiāng)融合發(fā)展水平接近的區(qū)域大多呈連片分布；城鄉(xiāng)融合發(fā)展水平的相對高值區(qū)主要分布在長三角地區(qū)的中東部區(qū)域，而其相對低值區(qū)則主要分布在西部區(qū)域以及南部部分城市。結合各中心區(qū)城市城鄉(xiāng)融合發(fā)展水平具體值的變化趨勢來看，2003—2008年城市間城鄉(xiāng)融合發(fā)展差距縮小，除上海和銅陵的城鄉(xiāng)融合發(fā)展水平有所降低外，其余城市城鄉(xiāng)融合發(fā)展水平均出現(xiàn)不同程度的提高。其中，泰州和湖州的城鄉(xiāng)融合發(fā)展等級由第三梯度發(fā)展為第二梯度，常州和嘉興則由第二梯度躋身為第一梯度；值得一提的是，盡管這一階段金華的城鄉(xiāng)融合發(fā)展水平有所提高，但由于其增長幅度有限，該城市的城鄉(xiāng)融合發(fā)展等級卻由第三梯度下降為第四梯度。2008—2013年期間城市間城鄉(xiāng)融合發(fā)展差距略微增加，無錫、溫州、湖州、嘉興、合肥、銅陵、滁州7座城市的城鄉(xiāng)融合發(fā)展水平有所增加，其余20座城市的城鄉(xiāng)融合發(fā)展水平均出現(xiàn)下降情況。合肥、金華、溫州的城鄉(xiāng)融合發(fā)展等級由第四梯度上升為第三梯度(2)這一階段金華的城鄉(xiāng)融合發(fā)展水平略微有所降低，但其下降幅度低于長三角地區(qū)整體水平的下降程度，故而該城市的城鄉(xiāng)融合發(fā)展等級由第四梯度上升為第三梯度。，湖州和杭州則由第二梯度發(fā)展為第一梯度。2013—2018年期間，長三角地區(qū)城市間城鄉(xiāng)融合發(fā)展差距再次縮小，南京、鹽城、鎮(zhèn)江、寧波、紹興、舟山、臺州、合肥、蕪湖、安慶和池州11座城市的城鄉(xiāng)融合發(fā)展水平出現(xiàn)不同程度的上升，其余16座城市的城鄉(xiāng)融合發(fā)展水平則出現(xiàn)下降。其中，銅陵和金華的城鄉(xiāng)融合發(fā)展等級由第三梯度下降為第四梯度，合肥、鹽城、臺州由第三梯度上升為第二梯度，南京和寧波由第二梯度躋身為第一梯度，而上海、常州、嘉興、杭州則由第一梯度下降為第二梯度。

4. 長三角地區(qū)城鄉(xiāng)融合發(fā)展水平的趨勢面分析

利用ArcGIS軟件中的全局趨勢面分析工具，將2003—2018年期間長三角地區(qū)城鄉(xiāng)融合發(fā)展水平作為高度屬性Z值，各中心區(qū)城市的地理坐標作為水平屬性X，Y值，采用二階趨勢面函數進行擬合，從全局層面擬合城鄉(xiāng)融合發(fā)展水平同地理位置的函數關系，擬合結果如圖8所示(圖中X軸代表東西方向，Y軸代表南北方向)。

圖8顯示：(1)整體來看，長三角地區(qū)城鄉(xiāng)融合發(fā)展水平在東西方向和南北方向均大致呈現(xiàn)出“倒U型”的空間分異特征。(2)從具體表征來看，東西方向城鄉(xiāng)融合發(fā)展水平從西部區(qū)域安慶、池州、合肥、銅陵等城市向中、東部區(qū)域常州、無錫、湖州、嘉興、蘇州、上海等城市大幅度遞增，而至寧波、舟山等最東部城市出現(xiàn)小幅度回落。南北方向城鄉(xiāng)融合發(fā)展水平則由南部區(qū)域溫州、臺州、金華、紹興等城市逐步過渡到中部上海、湖州、嘉興、蘇州、無錫、常州等城市，進而到北部區(qū)域泰州、揚州、滁州、鹽城等城市；其中以上海、蘇州、無錫等城市為代表的中部區(qū)域城鄉(xiāng)融合發(fā)展水平最高，分布在“倒U型”曲線的波峰位置，且保持穩(wěn)定分布。(3)從各年空間分異的變化趨勢來看，2003—2018年長三角地區(qū)城鄉(xiāng)融合發(fā)展的空間分異規(guī)律基本保持不變。東西方向上，西部、中部和東部區(qū)域城市間城鄉(xiāng)融合發(fā)展差距呈先擴大再縮小的趨勢，至2018年東部區(qū)域城鄉(xiāng)融合發(fā)展的相對水平明顯提高，出現(xiàn)由西向東遞增演變的特征；南北方向上，南部、中部和北部區(qū)域城市間的城鄉(xiāng)融合發(fā)展差距也呈現(xiàn)出先擴大再縮小的趨勢，至2018年北部區(qū)域城鄉(xiāng)融合發(fā)展的相對水平明顯提高，大致呈現(xiàn)由南向北遞增演變的特征。

六、 結論與啟示

本文在科學闡釋城鄉(xiāng)融合發(fā)展時代內涵的基礎上，應用加速遺傳算法投影尋蹤模型對2003—2018年長三角地區(qū)27座中心區(qū)城市的城鄉(xiāng)融合發(fā)展水平進行了測度，進而運用核密度估計、探索性空間分析、趨勢面分析等方法分析了該地區(qū)城鄉(xiāng)融合發(fā)展的時空特征。

研究的主要結論如下：(1)2003—2018年長三角地區(qū)城鄉(xiāng)融合發(fā)展水平總體上呈波動上升趨勢，且體現(xiàn)出先上升后下降繼而水平波動的演進特征；各中心區(qū)城市間城鄉(xiāng)融合發(fā)展水平差異呈先縮小再擴大繼而再縮小的變化趨勢，且低值集聚分布明顯。(2)蘇州、無錫、上海的城鄉(xiāng)融合發(fā)展水平明顯高于其他城市，而安慶、池州和滁州則為城鄉(xiāng)融合水平排名最后三位的城市。(3)研究期內，該地區(qū)城鄉(xiāng)融合發(fā)展水平空間分布存在顯著的正向全局空間自相關性，且呈現(xiàn)出“弱相關→強相關→弱相關→強相關”的螺旋演進態(tài)勢。該地區(qū)城鄉(xiāng)融合發(fā)展的整體水平形成了以上海、蘇州、無錫、嘉興等城市為代表的高值簇(HH)和以池州、安慶、蕪湖、馬鞍山等城市為代表的低值簇(LL)為主的局部空間集聚特征，在空間上大致呈現(xiàn)為以上海和安慶為端點，由東向西逐級遞減的階梯狀分布格局，且低水平集聚狀態(tài)的城市正在朝高水平集聚狀態(tài)的方向緩慢發(fā)展。(4)長三角地區(qū)城鄉(xiāng)融合發(fā)展水平在東西方向和南北方向均大致呈現(xiàn)出“倒U型”的空間分異特征。

根據上述研究結論，本文得出以下幾點政策啟示：第一，應對長三角地區(qū)城鄉(xiāng)融合發(fā)展的現(xiàn)實狀態(tài)及時進行評估，綜合考慮各地區(qū)城鄉(xiāng)資源稟賦、經濟發(fā)展階段和產業(yè)結構特點等方面的差異性，因城施策制定促進各地區(qū)城鄉(xiāng)融合發(fā)展的方針政策。第二，應適當增加對長三角西部和南部區(qū)域城鄉(xiāng)發(fā)展的政策傾斜和資金支持，通過出臺具體政策引導城鄉(xiāng)融合發(fā)展水平相對落后的城市向東部、北部城市借鑒城鄉(xiāng)發(fā)展經驗，并加大資金投入，改善西部和南部地區(qū)農村公共服務和基礎設施供給。第三，應制定區(qū)域內部城鄉(xiāng)發(fā)展定向合作計劃，加強長三角地區(qū)城際協(xié)同合作，充分發(fā)揮東北部地區(qū)對西南部地區(qū)城鄉(xiāng)融合發(fā)展的溢出效應。