陳曉華，王千辰

(1.安徽建筑大學 建筑與規(guī)劃學院，安徽 合肥 230601；2.安徽省城鎮(zhèn)化發(fā)展研究中心，安徽 合肥 230601)

一、引 言

在工業(yè)化、城鎮(zhèn)化快速推進過程中，鄉(xiāng)村大量勞動力、土地、公共資源、資本等要素向城市集聚[1]，促使城鄉(xiāng)發(fā)展二元分化，導致城鄉(xiāng)差距拉大、空間生產結構失衡[2-3]、鄉(xiāng)村空心化[4]、老齡化[5]等一系列問題。為加快破除城鄉(xiāng)分割的體制機制，重塑新型城鄉(xiāng)關系，促進城鄉(xiāng)融合一體化發(fā)展，黨的十九屆五中全會提出，全面實施鄉(xiāng)村振興戰(zhàn)略，強化以工補農、以城帶鄉(xiāng)，推動形成工農互促、城鄉(xiāng)互補、協(xié)調發(fā)展、共同繁榮的新型工農城鄉(xiāng)關系。城鄉(xiāng)融合發(fā)展是社會經濟結構轉變、城鄉(xiāng)空間重構的過程，通過城鄉(xiāng)之間各要素合理流動和有效配置，實現城鄉(xiāng)要素趨同和城鄉(xiāng)之間協(xié)同發(fā)展，促進城鄉(xiāng)之間多維度的重構，實現城鄉(xiāng)地域平等共生的新型城鄉(xiāng)關系[6]。近年，相關學者們從不同的角度和尺度對城鄉(xiāng)融合發(fā)展進行廣泛的研究：一是以城鄉(xiāng)關系為演變背景闡述城鄉(xiāng)融合發(fā)展的內涵[7]；二是探析城鄉(xiāng)融合背景下城鄉(xiāng)發(fā)展存在的問題和對策[8]；三是通過探索城鄉(xiāng)融合發(fā)展空間格局和影響因素分析[9-10]，揭示城鄉(xiāng)空間融合過程的動力機制[11]。既有研究多關注宏觀區(qū)域尺度城鄉(xiāng)融合發(fā)展水平測度，較少關注縣域城鄉(xiāng)融合發(fā)展的時空格局及其影響因素。本文以皖北地區(qū)35個縣級單元為例，構建城鄉(xiāng)融合發(fā)展水平評價指標體系，測度皖北地區(qū)2010-2019年城鄉(xiāng)融合發(fā)展水平，運用空間計量分析等多種方法揭示皖北地區(qū)城鄉(xiāng)融合發(fā)展的時空特征與影響因素。

二、研究區(qū)域與研究設計

1.研究區(qū)域與數據來源

(1)研究區(qū)域概況。皖北地區(qū)位于安徽省北部，地處皖蘇魯豫四省交匯區(qū)，地勢平坦，以平原地形為主，行政區(qū)劃范圍包括安徽省淮河沿岸及以北的6個地級市，包括宿州市、淮北市、亳州市、阜陽市、蚌埠市、淮南市。皖北地區(qū)總面積約為4.278萬平方公里，耕地面積約2.711萬平方公里，耕地占比63.371%；2019年末城鎮(zhèn)常住人口1 417.724萬人，鄉(xiāng)村常住人口1 421.725萬人，城鎮(zhèn)化率為49.930%；第一產業(yè)產值1 313.371億元，第二、三產業(yè)產值9 550.656億元。皖北地區(qū)是鄉(xiāng)村人口占比大、城鄉(xiāng)經濟差距較大的傳統(tǒng)農業(yè)發(fā)展區(qū)，對研究城鄉(xiāng)融合發(fā)展的問題具有典型意義(圖1)。

圖1 研究區(qū)域①

(2)數據來源與處理。數據來源于《安徽省統(tǒng)計年鑒》《安徽省建設統(tǒng)計年鑒》以及各市統(tǒng)計年鑒、統(tǒng)計公報和郵政業(yè)發(fā)展公報等。對于部分缺失的數據，采用均值填補法對缺失項進行填補。

2.研究設計

(1)研究框架。要素流動理論是在吸收早期西方經濟學的資源稟賦理論、比較優(yōu)勢理論以及勞動價值論的認知基礎之上形成的，主要指生產要素在地域空間的位移現象，表現為生產要素在地域空間上優(yōu)化配置[12]。通過城鄉(xiāng)要素雙向流動，打破因城鄉(xiāng)要素單向流動導致的城鄉(xiāng)二元發(fā)展格局，實現城市與鄉(xiāng)村之間要素雙向流動與優(yōu)勢互補，促使各要素在城市與鄉(xiāng)村之間合理配置，優(yōu)化城市與鄉(xiāng)村地域的空間、人口、經濟、社會、生態(tài)結構，推動新型城鎮(zhèn)化質量提升與城鄉(xiāng)融合發(fā)展(圖2)。

圖2 要素流動視角下城鄉(xiāng)融合發(fā)展

(2)研究方法。研究基于要素流動視角下城鄉(xiāng)融合發(fā)展機制構建城鄉(xiāng)融合發(fā)展的指標體系與評價模型，繼而運用空間自相關等方法揭示城鄉(xiāng)融合發(fā)展時空演變特征，最后運用地理加權回歸模型探析城鄉(xiāng)融合發(fā)展主要影響因素的作用強度。

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

①城鄉(xiāng)融合發(fā)展評價指標選取與權重賦值?；谇笆鲆亓鲃优c城鄉(xiāng)融合發(fā)展的互動機制，從城鄉(xiāng)系統(tǒng)的人口、土地、社會服務、經濟及生態(tài)等要素配置視角，選擇人口流動、空間配置、社會服務、經濟發(fā)展、生態(tài)環(huán)境等五個方面作為城鄉(xiāng)融合發(fā)展評價的準則層。繼而基于數據的科學性、可獲取性、系統(tǒng)性等原則，參考既有的相關研究成果[13-16]，運用文獻指標借鑒、頻度統(tǒng)計法等方法，最終選取城鄉(xiāng)融合發(fā)展評價指標，包括人口城鎮(zhèn)化水平、城鄉(xiāng)勞動力流動、城鄉(xiāng)人口遷移率、城鄉(xiāng)人口密度比、土地城鎮(zhèn)化水平、城鄉(xiāng)土地配置、城鄉(xiāng)醫(yī)療水平度、城鄉(xiāng)教育發(fā)展度、城鄉(xiāng)物流業(yè)務量、城鄉(xiāng)居民用電量比、城鄉(xiāng)人均可支配收入比、產業(yè)產值比、城鄉(xiāng)消費品零售總額比、城鄉(xiāng)資金流動、人均GDP、森林覆蓋率、單位GDP能耗等17個指標。為避免量綱、數量級以及正負指標的差異，對城鄉(xiāng)融合發(fā)展指標進行無量綱化處理，運用熵值法計算出各指標權重，最終確定城鄉(xiāng)融合發(fā)展水平測度指標(表1)。

②空間自相關模型。地理學第一定律強調事物間相互關聯(lián)性，研究空間中某個位置的觀察值與其相鄰位置是否存在相關性及相關程度的一種空間數據分析方法[17]21-22。本文運用空間自相關模型揭示城鄉(xiāng)融合發(fā)展的空間特征，以全局Moran’sI指數測度城鄉(xiāng)融合發(fā)展整體空間關聯(lián)[18]，計算公式如下：

(1)

(2)

(3)

③地理加權回歸模型。地理加權回歸(GWR)模型反映出變量間關系的局部特征，考慮同一指標對不同空間的回歸結果的不同影響。運用地理加權回歸模型來探究城鄉(xiāng)融合發(fā)展的影響因素隨地理位置變化而呈現的回歸參數大小和方向的具體變化[20]，構建GWR模型。

(4)

式中，γi為城鄉(xiāng)融合發(fā)展水平測度值，xij為城鄉(xiāng)融合發(fā)展的影響因素，(ui,vi)為第i個樣點的空間坐標，βj(ui,vi)是第i個樣本點上的第k個回歸參數。

三、結果與分析

1.皖北地區(qū)城鄉(xiāng)融合發(fā)展水平整體呈上升態(tài)勢

運用城鄉(xiāng)融合發(fā)展評價模型對2010-2019年皖北地區(qū)6市的城鄉(xiāng)融合發(fā)展水平進行測度(圖3)。結果顯示，2010-2019年皖北地區(qū)城鄉(xiāng)融合發(fā)展水平整體呈上升態(tài)勢。根據城鄉(xiāng)融合發(fā)展水平變化速率不同，可將皖北地區(qū)城鄉(xiāng)融合發(fā)展過程劃分為兩個階段：第一階段(2010-2012年)，皖北地區(qū)城鄉(xiāng)融合發(fā)展水平測度值不足0.350，城鄉(xiāng)融合發(fā)展緩慢，城鄉(xiāng)融合發(fā)展水平較低；第二階段(2013-2019年)，城鄉(xiāng)融合發(fā)展水平測度值由0.358上升至0.624，城鄉(xiāng)融合發(fā)展加快，城鄉(xiāng)融合發(fā)展水平較高。

圖3 2010-2019年皖北地區(qū)城鄉(xiāng)融合發(fā)展水平變化

從城鄉(xiāng)融合發(fā)展水平增長速度看，蚌埠市、淮北市和阜陽市城鄉(xiāng)融合發(fā)展水平高于皖北地區(qū)整體平均水平，2013年以來蚌埠市和阜陽市城鄉(xiāng)融合發(fā)展水平增長最快，2019年蚌埠市城鄉(xiāng)融合發(fā)展水平居于皖北地區(qū)6市首位，高達0.745。近10年淮南市城鄉(xiāng)融合發(fā)展水平提升最慢，2019年居皖北六市最低位次，城鄉(xiāng)融合發(fā)展水平測度值僅為0.536；亳州市和宿州市城鄉(xiāng)融合發(fā)展水平低于皖北地區(qū)整體水平，亳州市城鄉(xiāng)融合發(fā)展水平由2010年的0.156上升至2019年的0.616，城鄉(xiāng)融合發(fā)展水平上升速度最快。

2.縣域城鄉(xiāng)融合發(fā)展水平呈空間差異性

對2019年皖北地區(qū)35個縣級單元運用城鄉(xiāng)融合發(fā)展水平評價模型進行測度分析，得到各縣級行政單元的城鄉(xiāng)人口流動、城鄉(xiāng)空間配置、城鄉(xiāng)社會服務、城鄉(xiāng)經濟發(fā)展、城鄉(xiāng)生態(tài)環(huán)境和綜合城鄉(xiāng)融合發(fā)展水平測度值(表2)。結果顯示，譙城區(qū)的城鄉(xiāng)融合發(fā)展水平測度值最高，達0.574；相山區(qū)、龍子湖區(qū)、禹會區(qū)、蚌山區(qū)、淮上區(qū)、潁州區(qū)等6個縣級單元的城鄉(xiāng)融合發(fā)展水平都處于較高水平，城鄉(xiāng)融合發(fā)展水平測度值均在0.5以上。壽縣的城鄉(xiāng)融合發(fā)展水平測度值最低，為0.201；泗縣、烈山區(qū)、濉溪縣、八公山區(qū)、潘集區(qū)等5個縣級單元的城鄉(xiāng)融合發(fā)展水平均處于較低水平，城鄉(xiāng)融合發(fā)展水平測度值不足0.350。

城鄉(xiāng)融合發(fā)展水平測度在城鄉(xiāng)人口流動、城鄉(xiāng)空間配置、城鄉(xiāng)社會服務、城鄉(xiāng)經濟發(fā)展、城鄉(xiāng)生態(tài)環(huán)境五個維度都有不同程度的表現，各維度的發(fā)展水平呈空間差異性。其中，城鄉(xiāng)人口流動發(fā)展水平測度均值為0.085，最高值太和縣(0.139)與最低值壽縣(0.016)；城鄉(xiāng)空間配置發(fā)展水平測度均值為0.058，最高值龍子湖區(qū)(0.089)與最低值壽縣(0.005)，表明皖北地區(qū)鄉(xiāng)村大量勞動力流向城市，導致鄉(xiāng)村土地資源低效和閑置問題普遍突出；城鄉(xiāng)社會服務發(fā)展水平測度均值為0.121，最高值譙城區(qū)(0.242)與最低值八公山區(qū)(0.020)，社會服務水平差距較大，表明皖北地區(qū)城鄉(xiāng)之間社會公共服務發(fā)展不均等；城鄉(xiāng)經濟發(fā)展水平測度均值為0.129，其中最高值禹會區(qū)(0.218)、最低值臨泉縣(0.075)，皖北地區(qū)縣域城鄉(xiāng)經濟發(fā)展空間差異化明顯；各縣區(qū)城鄉(xiāng)生態(tài)環(huán)境發(fā)展水平測度值較為接近，表明皖北地區(qū)城鄉(xiāng)生態(tài)環(huán)境發(fā)展呈均等化。

表2 2019年皖北地區(qū)縣級單元城鄉(xiāng)融合發(fā)展水平測度結果

3.區(qū)域城鄉(xiāng)融合發(fā)展空間格局呈集聚性

為揭示皖北地區(qū)城鄉(xiāng)融合發(fā)展空間格局，運用ArcGIS10.2進行全局空間自相關分析，得到2019年皖北地區(qū)35個縣級單元的城鄉(xiāng)融合發(fā)展水平的全局空間自相關Moran’sI值(表3)。結果顯示，2019年皖北地區(qū)縣級單元城鄉(xiāng)融合水平Moran’sI指數為0.212，Z值為2.353，P值為0.019。說明皖北地區(qū)的城鄉(xiāng)融合發(fā)展水平在整體上表現為空間正相關，空間集聚性較強。

城鄉(xiāng)人口流動、城鄉(xiāng)空間配置、城鄉(xiāng)社會服務、城鄉(xiāng)經濟發(fā)展、城鄉(xiāng)生態(tài)環(huán)境五個維度的空間格局特征不同。城鄉(xiāng)社會服務、城鄉(xiāng)經濟發(fā)展、城鄉(xiāng)生態(tài)環(huán)境Moran’sI指數均大于0，Z值大于2.58，且全部通過1%的顯著性水平檢測，說明皖北地區(qū)的城鄉(xiāng)社會服務、城鄉(xiāng)經濟發(fā)展、城鄉(xiāng)生態(tài)環(huán)境三個維度表現為空間正相關，空間集聚性強；城鄉(xiāng)人口流動和城鄉(xiāng)空間配置未通過5%的顯著性水平檢驗，說明皖北地區(qū)城鄉(xiāng)人口流動和城鄉(xiāng)空間配置數據分布出現隨機分布特征，空間相關性不顯著。

表3 2019年皖北地區(qū)縣級單元城鄉(xiāng)融合發(fā)展水平Moran’sI值

圖4 皖北地區(qū)城鄉(xiāng)融合發(fā)展冷熱點分析(2019年)

4.皖北地區(qū)城鄉(xiāng)融合發(fā)展的影響因素

(1)影響因素變量篩選。為確定皖北地區(qū)城鄉(xiāng)融合發(fā)展水平的主要影響因素，對皖北地區(qū)城鄉(xiāng)融合發(fā)展的空間格局形成原因進行深層次分析。以皖北地區(qū)35個縣級單元為樣本研究，以2019年城鄉(xiāng)融合發(fā)展水平測度值為因變量，從資源、區(qū)位、人口、社會、經濟等維度選取人均水資源(X1)、距所在市中心距離(X2)、距省會城市距離(X3)、城鄉(xiāng)從業(yè)人口率比(X4)、城鄉(xiāng)居民人均用電量比(X5)、城鄉(xiāng)居民醫(yī)療保險參保率(X6)、城鄉(xiāng)戶籍人口遷移率(X7)、第一產業(yè)產值占比(X8)、財政收支占比(X9)、土地墾殖率(X10)、城鎮(zhèn)登記失業(yè)率(X11)、城鄉(xiāng)郵遞線路密度(X12)、耕地面積占比(X13)等13個指標作為影響因素備選變量。

運用最小二乘法OLS回歸模型測算變量對被解釋變量(城鄉(xiāng)融合發(fā)展水平測度值)的影響顯著性水平。剔除顯著性(Robust-Pr)大于0.05以及方差膨脹因子(VIF)大于7.5的指標，結果顯示，距省會城市距離(X3)、財政收支占比(X9)、城鎮(zhèn)登記失業(yè)率(X11)、城鄉(xiāng)郵遞線路密度(X12)共計4個指標呈現影響顯著且未出現多重共線性回歸問題，符合OLS回歸模型，是皖北地區(qū)城鄉(xiāng)融合發(fā)展的主要影響因素(表4)。

表4 OLS回歸模型分析結果

(2)主要影響因素作用強度空間差異。運用地理加權回歸(GWR)對篩選出的4個自變量的回歸系數進行計算，探討主要影響因素對城鄉(xiāng)融合發(fā)展作用強度的空間差異，并用ArcGIS進行可視化表達(圖5)。

圖5 城鄉(xiāng)融合發(fā)展影響因素回歸系數空間分布(2019年)

結果顯示，距省會城市距離與城鄉(xiāng)融合發(fā)展水平呈正相關，在空間上呈現為由北向南逐漸遞增的空間分布特征，回歸系數最大值為壽縣，最小值為碭山縣，說明距省會城市越近，城鄉(xiāng)融合發(fā)展水平受其影響也越大；財政收支占比與城鄉(xiāng)融合發(fā)展水平呈正相關，在空間上表現為由東南向西北逐漸遞減的空間分布特征，回歸系數最大值為五河縣，最小值為界首市，表明財政收支占比對五河縣、泗縣等地的城鄉(xiāng)融合發(fā)展影響較大，對界首市、臨泉縣等北部地區(qū)的城鄉(xiāng)融合發(fā)展影響較??；城鄉(xiāng)郵遞線路密度與城鄉(xiāng)融合發(fā)展水平呈正相關，在空間上呈現為由西北向東南逐漸遞增的空間分布特征，回歸系數最大值為壽縣，最小值為碭山縣，說明城鄉(xiāng)郵遞線路密度對壽縣、大通區(qū)、五河縣等東南地區(qū)的城鄉(xiāng)融合發(fā)展影響高于其他地區(qū)。城鎮(zhèn)登記失業(yè)率與城鄉(xiāng)融合發(fā)展水平呈負相關，在空間分布上總體呈現由東北向西南遞減的特征，回歸系數最大值泗縣，最小值為臨泉縣，表明城鎮(zhèn)登記失業(yè)率對泗縣、靈璧縣、蕭縣等東北部地區(qū)的影響較小，對臨泉縣、阜南縣等西南地區(qū)的影響較大。

四、結論與建議

1.結論

(1)從要素流動視角對皖北地區(qū)城鄉(xiāng)融合發(fā)展水平進行測度，揭示了近10年皖北地區(qū)城鄉(xiāng)融合發(fā)展水平時空特征。第一，城鄉(xiāng)融合發(fā)展水平整體呈加速上升態(tài)勢，城鄉(xiāng)融合發(fā)展水平的測度值由2010年的0.282上升至2013年的0.358，2019年上升至0.624。第二，縣域城鄉(xiāng)融合發(fā)展水平空間差異性顯著，其中譙城區(qū)城鄉(xiāng)融合發(fā)展水平測度值最高為0.574，壽縣城鄉(xiāng)融合發(fā)展水平最低為0.201。第三，皖北地區(qū)城鄉(xiāng)融合發(fā)展水平Moran’sI指數為0.212，表明區(qū)域城鄉(xiāng)融合發(fā)展空間格局呈較強的集聚性。

(2)城鄉(xiāng)融合發(fā)展是通過城鄉(xiāng)要素流動，推動城鄉(xiāng)地域結構和城鄉(xiāng)發(fā)展結構不斷優(yōu)化的結果。研究發(fā)現，皖北地區(qū)城鄉(xiāng)融合發(fā)展水平受區(qū)域城鄉(xiāng)社會經濟多重要素的綜合影響，距省會城市距離、財政收支占比、城鄉(xiāng)郵遞線路密度及城鎮(zhèn)登記失業(yè)率等四個因素影響最顯著，其中距省會城市距離、財政收支占比、城鄉(xiāng)郵遞線路密度與皖北地區(qū)城鄉(xiāng)融合發(fā)展水平測度呈現正相關性；城鎮(zhèn)登記失業(yè)率與城鄉(xiāng)融合發(fā)展水平測度呈負相關性。這些因素對皖北地區(qū)城鄉(xiāng)融合發(fā)展的影響均具有明顯的空間差異。

2.建議

基于城鄉(xiāng)融合發(fā)展水平的分析結果，為了進一步提高皖北地區(qū)城鄉(xiāng)融合發(fā)展水平，提出皖北地區(qū)城鄉(xiāng)融合發(fā)展的建議。

(1)注重城鄉(xiāng)融合發(fā)展差異，制定分類融合發(fā)展政策。對于譙城區(qū)、蚌山區(qū)等城鄉(xiāng)融合發(fā)展水平較高但生態(tài)環(huán)境水平較低地區(qū)，應注重改善生態(tài)環(huán)境質量，并充分發(fā)揮中心城市的輻射帶動效應，實現城鄉(xiāng)融合高質量發(fā)展。對于壽縣、泗縣等城鄉(xiāng)融合發(fā)展水平較低地區(qū)，加大城鄉(xiāng)基礎設施、公共服務均等化建設力度，注重城市與鄉(xiāng)村協(xié)同發(fā)展。

(2)建立城鄉(xiāng)要素流動機制，促進城鄉(xiāng)要素雙向流動。積極推動皖北地區(qū)土地要素市場化，開發(fā)利用荒廢土地，盤活鄉(xiāng)村土地要素活力；在加快農村剩余勞動力向城市轉移的同時，建立人才、技術下鄉(xiāng)的激勵機制；建立皖北地區(qū)城鄉(xiāng)要素配置、優(yōu)勢互補機制，促進城鄉(xiāng)融合發(fā)展，提升城鄉(xiāng)要素流動效率。

(3)完善財政支出政策，促進城鄉(xiāng)產業(yè)協(xié)調發(fā)展。加大鄉(xiāng)村建設的資金投入力度，深度挖掘多功能農業(yè)，提升農業(yè)現代化水平；加快城鄉(xiāng)第二產業(yè)和第三產業(yè)發(fā)展，促進產業(yè)結構升級，推動就業(yè)崗位穩(wěn)定增長，優(yōu)化城鄉(xiāng)人力資本結構。