基于ESDA的山東省四維城鎮(zhèn)化協(xié)調(diào)發(fā)展空間差異

吳云青 王多多 郭佳

摘要：在構(gòu)建評價(jià)指標(biāo)體系的基礎(chǔ)上，運(yùn)用多指標(biāo)綜合評價(jià)法、協(xié)調(diào)發(fā)展度模型、探索性空間分析等方法探討山東省人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間差異。結(jié)果表明，大多數(shù)地市的城鎮(zhèn)化都處于中度及以上協(xié)調(diào)發(fā)展?fàn)顟B(tài);人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度存在顯著的正空間自相關(guān)特點(diǎn)，城鎮(zhèn)化協(xié)調(diào)發(fā)展的均質(zhì)性區(qū)域明顯多于異質(zhì)性區(qū)域，城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間依賴性和差異性同時(shí)存在;H-H聚集類型區(qū)主要集中在東部沿海和魯中地區(qū)，H-L聚集類型區(qū)內(nèi)只有濟(jì)南市;L-L聚集類型區(qū)、L-H聚集類型區(qū)主要分布在魯西和魯南地區(qū)。

關(guān)鍵詞：城鎮(zhèn)化;協(xié)調(diào)發(fā)展;空間差異;探索性空間分析

中圖分類號： F291.3 ?文獻(xiàn)標(biāo)志碼： A ?文章編號：1002-1302（2019）13-0337-05

長期以來，我國的城鎮(zhèn)化采用傳統(tǒng)模式，質(zhì)量不高，存在諸多限制性因素，尤其表現(xiàn)為城鎮(zhèn)化內(nèi)部各要素不協(xié)調(diào)、不匹配[1]，經(jīng)濟(jì)城鎮(zhèn)化不穩(wěn)，土地城鎮(zhèn)化冒進(jìn)，人口城鎮(zhèn)化滯后，社會城鎮(zhèn)化脆弱[2]。黨的十八屆三中全會提出要堅(jiān)持走中國特色新型城鎮(zhèn)化道路，而新型城鎮(zhèn)化是內(nèi)涵式發(fā)展的城鎮(zhèn)化，離不開其內(nèi)部要素相互協(xié)調(diào)、良性互動(dòng)，城鎮(zhèn)化內(nèi)部各要素能否協(xié)調(diào)發(fā)展及其協(xié)調(diào)發(fā)展程度決定了新型城鎮(zhèn)化進(jìn)程的可持續(xù)性及其效益。國內(nèi)相關(guān)研究大多將城鎮(zhèn)化劃分為人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化4個(gè)要素;國外并沒有我國這種“要素式”城鎮(zhèn)化內(nèi)涵界定與分類，更沒有城鎮(zhèn)化內(nèi)部各要素協(xié)調(diào)發(fā)展的研究，但他們對人口遷移[3]、城市空間擴(kuò)張[4]、城市人口增長[5]的研究成果可為本研究提供有益借鑒。近年來，著眼于“健康、有序、綠色、可持續(xù)”城鎮(zhèn)化，田時(shí)中等對人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展?fàn)顩r開展了系列研究，但這些研究多停留在協(xié)調(diào)度測量層面，聚焦于要素協(xié)調(diào)發(fā)展空間差異的研究并不多見[6-8]。另外，已有為數(shù)不多的城鎮(zhèn)化協(xié)調(diào)發(fā)展空間差異的研究也大多采用傳統(tǒng)分析方法，忽視了空間因素的作用，缺乏對空間信息的有效體現(xiàn)，對城鎮(zhèn)化協(xié)調(diào)發(fā)展在地理位置上的關(guān)聯(lián)度、聚集性的研究較少，難以揭示人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間特征。因此，本研究以山東省城鎮(zhèn)化發(fā)展的實(shí)際狀況為例，構(gòu)建人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化各要素評價(jià)指標(biāo)體系，利用協(xié)調(diào)發(fā)展度模型、探索性空間分析（ESDA）研究山東省人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間差異，這對于有針對性地制定創(chuàng)新相關(guān)政策，促進(jìn)新型城鎮(zhèn)化健康、可持續(xù)發(fā)展，提升新型城鎮(zhèn)化質(zhì)量具有重要的現(xiàn)實(shí)意義。

1 評價(jià)指標(biāo)體系構(gòu)建及研究方法

1.1 指標(biāo)體系構(gòu)建

基于研究目的，遵循評價(jià)指標(biāo)構(gòu)建的系統(tǒng)性、科學(xué)性、代表性、可比性等原則，利用理論分析法、頻度統(tǒng)計(jì)法、專家咨詢法甄選人口、土地、經(jīng)濟(jì)、社會城鎮(zhèn)化評價(jià)指標(biāo)體系。為了能夠獲得客觀的原始數(shù)據(jù)，盡可能選擇統(tǒng)計(jì)年鑒、統(tǒng)計(jì)公報(bào)中有統(tǒng)計(jì)數(shù)據(jù)的指標(biāo)或依據(jù)統(tǒng)計(jì)年鑒、統(tǒng)計(jì)公報(bào)經(jīng)數(shù)學(xué)計(jì)算能夠得到數(shù)據(jù)的指標(biāo)。為便于比較不同區(qū)域城鎮(zhèn)化發(fā)展?fàn)顩r，指標(biāo)盡量少采用總量指標(biāo)，較多采用人均和相對指標(biāo)，構(gòu)建評價(jià)指標(biāo)體系（表1）。

1.2 研究方法

1.2.1 人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化發(fā)展評價(jià)方法

1.2.1.1 基于組合賦權(quán)法的評價(jià)指標(biāo)權(quán)重計(jì)算 本研究將主觀賦權(quán)的三角模糊法與客觀賦權(quán)的熵值法結(jié)合計(jì)算指標(biāo)組合權(quán)重。按照模糊數(shù)學(xué)原理，三角模糊法計(jì)算的標(biāo)權(quán)重取決于專家的模糊判斷。首先，由專家給出最悲觀估計(jì)、最可能估計(jì)和最樂觀估計(jì)，形成三角模糊合成矩陣;其次，運(yùn)用算子 M（，） 進(jìn)行模糊合成;最后，根據(jù)各指標(biāo)的模糊得分計(jì)算出最終權(quán)重，所得權(quán)重向量可表示為：

1.2.1.2 指標(biāo)數(shù)據(jù)的標(biāo)準(zhǔn)化處理 采用極差標(biāo)準(zhǔn)化方法對各個(gè)評價(jià)指標(biāo)數(shù)據(jù)進(jìn)行規(guī)范化處理，從而得到標(biāo)準(zhǔn)化值[10]。

1.2.1.3 人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化指數(shù)計(jì)算 根據(jù)各指標(biāo)標(biāo)準(zhǔn)化值及其對應(yīng)的組合權(quán)重，運(yùn)用多指標(biāo)綜合評價(jià)法分別計(jì)算人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化指數(shù)值。

1.2.2 人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度計(jì)算 引入?yún)f(xié)調(diào)發(fā)展度模型對山東省人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展情況進(jìn)行定量評價(jià)[11]。

式中：P表示綜合城鎮(zhèn)化指數(shù);D表示城鎮(zhèn)化耦合系數(shù);E表示城鎮(zhèn)化協(xié)調(diào)發(fā)展度;F1、F2、F3、F4分別表示人口、土地、經(jīng)濟(jì)、社會城鎮(zhèn)化指數(shù);a、b、c、d分別表示人口、土地、經(jīng)濟(jì)、社會城鎮(zhèn)化對城鎮(zhèn)化總系統(tǒng)的貢獻(xiàn)度，考慮到各子系統(tǒng)對城鎮(zhèn)化發(fā)展同等重要，故取a=b=c=d=1/4。

1.2.3 探索性空間數(shù)據(jù)分析（ESDA）

1.2.3.1 全局空間自相關(guān)分析 全局空間自相關(guān)體現(xiàn)了觀測變量在研究區(qū)域內(nèi)空間相關(guān)性的總體特征，可以定量表征某種地理現(xiàn)象在空間上的總體關(guān)聯(lián)與差異程度。本研究運(yùn)用Moran提出的Global Morans I指數(shù)衡量全局自相關(guān)情況[12]，計(jì)算公式為：

1.2.3.2 局部空間自相關(guān)分析 局部空間自相關(guān)可以分析觀測變量的局部空間變化，能具體測量每個(gè)區(qū)域與周邊區(qū)域的局部空間關(guān)聯(lián)和空間差異程度，用Moran散點(diǎn)圖和Local Morans I（LISA）指數(shù)來衡量。Moran散點(diǎn)圖的橫軸是變量在不同位置的觀測值向量，縱軸是該向量空間滯后向量的所有取值，每個(gè)空間滯后是該區(qū)域周邊觀測值的加權(quán)平均值，Moran散點(diǎn)圖分為4個(gè)象限。然而，Moran散點(diǎn)圖無法定量揭示各研究區(qū)域與其相鄰區(qū)域的自相關(guān)程度及其顯著性，Local Morans I（LISA）指數(shù)恰恰可彌補(bǔ)這一缺陷，可用于定量分析研究區(qū)域i與j之間的空間自相關(guān)性，計(jì)算公式為：

式中：m0=∑ni=1（xi-x）2/n，xi、xj分別表示觀測變量x在第i、j研究區(qū)域的觀測值;wij表示研究區(qū)域之間的空間權(quán)重矩陣，如果研究區(qū)域i與j相鄰，則wij=1，否則wij=0;x 表示觀測變量x的均值。統(tǒng)計(jì)檢驗(yàn)也采用標(biāo)準(zhǔn)化Z檢驗(yàn)。如果Morans Ii>0且顯著，表示某研究區(qū)域與周圍區(qū)域之間差異小，存在高-高值或低-低值空間聚集;如果Morans Ii<0且顯著，表示研究區(qū)域與周圍區(qū)域之間差異大，存在高（H）-低（L）值或低-高值空間聚集。將這些通過顯著性檢驗(yàn)的區(qū)域可視化到地圖上即可形成LISA集聚圖，可測量局部空間的異質(zhì)性，識別局部空間聚集的“熱點(diǎn)”與“冷點(diǎn)”。

2 實(shí)證結(jié)果與分析

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

作為全國的農(nóng)業(yè)大省、工業(yè)大省和經(jīng)濟(jì)大省，山東省近年來城鎮(zhèn)化進(jìn)程不斷推進(jìn)。但山東省城鎮(zhèn)化也具有“速度快、協(xié)調(diào)性不夠、區(qū)域差異明顯”的特點(diǎn)，其城鎮(zhèn)化進(jìn)程具有較強(qiáng)的現(xiàn)實(shí)背景和代表性，所得研究結(jié)論和啟示對于其他地區(qū)的城鎮(zhèn)化協(xié)調(diào)發(fā)展具有重要的參考意義。

本試驗(yàn)是基于2015年山東省地級市尺度的研究，共有17個(gè)研究單元，數(shù)據(jù)主要來源于2016年《中國城市統(tǒng)計(jì)年鑒》、2016年《山東省統(tǒng)計(jì)年鑒》、山東省2015年國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計(jì)公報(bào)等。

2.2 城鎮(zhèn)化協(xié)調(diào)發(fā)展度計(jì)算

2.2.1 評價(jià)指標(biāo)權(quán)重的計(jì)算 首先，分別運(yùn)用三角模糊法、熵值法計(jì)算指標(biāo)權(quán)重，得到權(quán)重向量u1j、u2j;其次，利用公式（4）至公式（6）可以得到人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化評價(jià)中三角模糊法和熵值法的組合系數(shù)θk分別為0.364、0.636，0.198、0.802，0.244、0.756，0.194、0.806;最后，依據(jù)公式（3）可以計(jì)算得到各個(gè)指標(biāo)的組合權(quán)重（表1）。

2.2.2 城鎮(zhèn)化指數(shù)計(jì)算 在對各評價(jià)指標(biāo)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理的基礎(chǔ)上，可計(jì)算出2015年山東省人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化指數(shù)值（表2）。

2.2.3 城鎮(zhèn)化協(xié)調(diào)發(fā)展度計(jì)算 利用公式（8）至公式（10）可計(jì)算出山東省2015年綜合城鎮(zhèn)化指數(shù)、城鎮(zhèn)化耦合系數(shù)和城鎮(zhèn)化協(xié)調(diào)發(fā)展度（表2）。本研究將山東省人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度（E）劃分為4個(gè)等級：0

2.3 城鎮(zhèn)化協(xié)調(diào)發(fā)展度計(jì)算結(jié)果分析

以山東17個(gè)市作為橫軸標(biāo)簽，以表2中人口城鎮(zhèn)化指數(shù)、土地城鎮(zhèn)化指數(shù)、經(jīng)濟(jì)城鎮(zhèn)化指數(shù)、社會城鎮(zhèn)化指數(shù)、綜合城鎮(zhèn)化指數(shù)為主縱坐標(biāo)軸，以城鎮(zhèn)化耦合系數(shù)、城鎮(zhèn)化協(xié)調(diào)發(fā)展度為次縱坐標(biāo)軸，可得到2015年山東省17個(gè)市城鎮(zhèn)化協(xié)調(diào)發(fā)展度變量的二維坐標(biāo)圖（圖1）。

由圖1可知，山東17個(gè)市之間城鎮(zhèn)化耦合度、協(xié)調(diào)發(fā)展度差異明顯，但耦合系數(shù)和協(xié)調(diào)發(fā)展度2條折線的折拐方向基本一致，說明山東省城鎮(zhèn)化耦合度、協(xié)調(diào)發(fā)展度的空間布局特點(diǎn)較相似。濟(jì)南市城鎮(zhèn)化協(xié)調(diào)發(fā)展水平最高，其次依次是青島市、東營市和威海市，其協(xié)調(diào)發(fā)展度均大于0.8，屬于優(yōu)度協(xié)調(diào)發(fā)展。這4個(gè)市不僅各系統(tǒng)的城鎮(zhèn)化水平在全省位居前列，且發(fā)展較均衡，無明顯短板。淄博、煙臺、濰坊、萊蕪等市城鎮(zhèn)化協(xié)調(diào)發(fā)展度均大于0.6但小于0.8，屬于良好協(xié)調(diào)發(fā)展。其中，淄博市土地城鎮(zhèn)化水平較高，位居全省第2位，人口、經(jīng)濟(jì)和社會城鎮(zhèn)化水平均居全省上游;煙臺市土地、經(jīng)濟(jì)、社會城鎮(zhèn)化水平居全省中上游且發(fā)展較均衡，但人口城鎮(zhèn)化發(fā)展相對滯后，列全省第11位。濰坊市、萊蕪市均屬于在城鎮(zhèn)化某單一系統(tǒng)表現(xiàn)突出，但在其他方面發(fā)展較差。濱州、日照、泰安、棗莊、濟(jì)寧等市處于中度協(xié)調(diào)發(fā)展?fàn)顟B(tài)，其各系統(tǒng)

城鎮(zhèn)化水平大多位列全省中下游水平，但發(fā)展較均衡。臨沂、德州、聊城、菏澤等市處于低度協(xié)調(diào)發(fā)展?fàn)顟B(tài)，其各維度城鎮(zhèn)化水平都處于中下游位置，其中菏澤市更是各方面都居末位。

2.4 城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間差異分析

2.4.1 全局空間自相關(guān)分析 基于上述協(xié)調(diào)發(fā)展度計(jì)算結(jié)果，依據(jù)公式（10）并借助OpenGeoDa軟件可以計(jì)算出2015年山東省人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度的Global Morans I指數(shù)為0.191，運(yùn)用公式（11）的Z值法對空間自相關(guān)進(jìn)行顯著性檢驗(yàn)，結(jié)果表明Global Morans I指數(shù)的正態(tài)統(tǒng)計(jì)量Z值為2.00，大于0.05置信水平的臨界值1.96，顯著性檢驗(yàn)通過。說明山東省人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度存在顯著的正空間自相關(guān)特點(diǎn)。

2.4.2 局部空間自相關(guān)分析 上述全局空間自相關(guān)分析無法有效識別空間聚類的內(nèi)部狀態(tài)，也可能會隱藏偏離整體布局特征的空間關(guān)系。本研究利用Moran散點(diǎn)圖和Local Morans I（LISA）指數(shù)進(jìn)一步分析山東省人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展的局部空間特征。以城鎮(zhèn)化協(xié)調(diào)發(fā)展度的標(biāo)準(zhǔn)化向量ZE為橫坐標(biāo)，以城鎮(zhèn)化協(xié)調(diào)發(fā)展度滯后向量 W_ZE 為縱坐標(biāo)，以散點(diǎn)橫縱坐標(biāo)的均值為中心坐標(biāo)，可以得到包括4個(gè)基本象限的Moran散點(diǎn)圖（圖2），每個(gè)象限對應(yīng)不同的城鎮(zhèn)化協(xié)調(diào)發(fā)展度局部空間聚集類型。由圖2可知，山東省人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度散點(diǎn)主要分布在H-H聚集類型區(qū)、L-L聚集類型區(qū)，約占全部散點(diǎn)的 70.59%，而H-L聚集類型區(qū)、L-H聚集類型區(qū)的散點(diǎn)較少，約占全部散點(diǎn)的29.41%，山東省城鎮(zhèn)化協(xié)調(diào)發(fā)展的均質(zhì)性區(qū)域明顯多于異質(zhì)性區(qū)域，這也驗(yàn)證了上述全局空間自相關(guān)的分析結(jié)果。

利用OpenGeoDa軟件計(jì)算出2015年山東省城鎮(zhèn)化協(xié)調(diào)發(fā)展度的LISA值，在Z檢驗(yàn)的基礎(chǔ)上（P<0.05）繪制出LISA聚集圖（圖3）。由圖3可知，山東省大部分地市的城鎮(zhèn)化協(xié)調(diào)發(fā)展度與相鄰市之間的關(guān)聯(lián)性并不顯著，只有青島市、煙臺市、濟(jì)南市和濟(jì)寧市的LISA值通過5%顯著性水平的Z統(tǒng)計(jì)檢驗(yàn)，相對于其他地市，這4個(gè)市對山東省城鎮(zhèn)化協(xié)調(diào)發(fā)展度空間自相關(guān)程度的貢獻(xiàn)更大。

（1）H-H聚集類型。屬于該類型的地市共有7個(gè)，分別是青島市、東營市、威海市、淄博市、煙臺市、濰坊市、萊蕪市，以面狀形式連片分布在山東半島區(qū)域和魯中地區(qū)。總體來說，該區(qū)域內(nèi)各市自身的城鎮(zhèn)化協(xié)調(diào)發(fā)展水平較高，周圍相鄰區(qū)域的城鎮(zhèn)化協(xié)調(diào)發(fā)展度也較高，成為高值集簇區(qū)。但值得

注意的是，盡管從全省層面看，該區(qū)域內(nèi)各市城鎮(zhèn)化協(xié)調(diào)發(fā)展度較高，但這7個(gè)市也存在較大差異，青島、東營、威海等市各維度城鎮(zhèn)化指數(shù)均位居全省前列，協(xié)調(diào)發(fā)展度均列全省前4位，但濰坊市和萊蕪市城鎮(zhèn)化各系統(tǒng)發(fā)展卻差異較大，優(yōu)點(diǎn)和缺點(diǎn)均較明顯，這2個(gè)市城鎮(zhèn)化協(xié)調(diào)度比較靠近全省協(xié)調(diào)度均值線，隨時(shí)都有可能落入L-H聚集類型區(qū)。LISA指數(shù)分析結(jié)果表明，該聚集類型區(qū)的7個(gè)市中，只有青島市和煙臺市的LISA值通過了顯著性檢驗(yàn)。一方面，說明青島市和煙臺市本身城鎮(zhèn)化協(xié)調(diào)發(fā)展能力較強(qiáng)，向東、西2個(gè)方向輻射到了威海、東營、濰坊、淄博等市;另一方面，也說明青島市和煙臺市與周邊各市同處山東半島沿海和膠濟(jì)鐵路沿線，空間聯(lián)系密切、文化風(fēng)俗相似，在城鎮(zhèn)化進(jìn)程中相互協(xié)作，良性競爭，共同推動(dòng)青島市和煙臺市成為城鎮(zhèn)化協(xié)調(diào)發(fā)展的“熱點(diǎn)”區(qū)。濰坊、東營、淄博、萊蕪等市的LISA值沒有通過顯著性檢驗(yàn)，說明這4個(gè)地市城鎮(zhèn)化協(xié)調(diào)發(fā)展能力還較弱，還不能有效輻射帶動(dòng)周邊的濱州、日照、臨沂、泰安等市區(qū)域城鎮(zhèn)化的協(xié)調(diào)發(fā)展。

（2）L-H聚集類型。屬于該類型的市有4個(gè)，分別是濱州市、泰安市、德州市和日照市，以點(diǎn)狀形態(tài)零散分布在魯西和魯南地區(qū)。從歷史基礎(chǔ)條件來看，這4個(gè)區(qū)域一直都是山東省的農(nóng)業(yè)大市，工業(yè)化、城鎮(zhèn)化基礎(chǔ)較差，本身的城鎮(zhèn)化協(xié)調(diào)發(fā)展度較低;而其周邊的東營、淄博、濟(jì)南、萊蕪、濰坊等市盡管城鎮(zhèn)化協(xié)調(diào)發(fā)展度較高，但其本身的輻射力還不夠強(qiáng)大，對于周邊區(qū)域的帶動(dòng)能力有限，因此這4個(gè)市成為城鎮(zhèn)化協(xié)調(diào)發(fā)展的低高空心區(qū)。

（3）L-L聚集類型。屬于該類型的市有5個(gè)，分別是臨沂市、棗莊市、濟(jì)寧市、聊城市、菏澤市，以面狀形態(tài)分布在魯西和魯南地區(qū)。該區(qū)域內(nèi)各市自身的城鎮(zhèn)化協(xié)調(diào)發(fā)展水平較低，而相鄰區(qū)域的水平也較低，相對落后的聚集特征明顯，成為低值蕭索區(qū)。一直以來，這些市都是山東省的欠發(fā)達(dá)地區(qū)，工業(yè)化基礎(chǔ)薄弱，城鎮(zhèn)化起步晚，產(chǎn)業(yè)支撐不足;這些市在地理位置上均是三面環(huán)繞山東省的欠發(fā)達(dá)地區(qū)，而南部又與河北省、河南省、江蘇省的欠發(fā)達(dá)市相鄰，缺乏發(fā)達(dá)城市體的輻射帶動(dòng)。LISA指數(shù)分析結(jié)果表明，該類型區(qū)的5個(gè)市中，只有濟(jì)寧市的LISA值通過了顯著性檢驗(yàn)，成為“冷點(diǎn)”區(qū)域，主要是因?yàn)橄鄬τ谄渌幱贚-L聚集類型區(qū)的地市，濟(jì)寧市的周邊為泰安市、聊城市、菏澤市、棗莊市、臨沂市，而這些地市均為城鎮(zhèn)化協(xié)調(diào)度的低值區(qū)域，低水平聚集的特征更加明顯。

（4）H-L聚集類型。屬于該類型的市只有濟(jì)南市，其LISA值也通過了顯著性檢驗(yàn)，說明濟(jì)南市自身城鎮(zhèn)化協(xié)調(diào)發(fā)展水平較高，位列全省第1位，但其周邊的德州、聊城、泰安、濱州等市協(xié)調(diào)發(fā)展度卻較低，成為高值孤立區(qū)。這種空間相關(guān)格局表明濟(jì)南市城鎮(zhèn)化發(fā)展的極化效應(yīng)顯著，長期以來濟(jì)南市都是山東省的政治、文化中心，但其城鎮(zhèn)化協(xié)調(diào)發(fā)展能力并沒有強(qiáng)大到可以有效帶動(dòng)周邊區(qū)域，與相鄰市沒有形成“近域擴(kuò)散”效應(yīng);相反，由于自身的優(yōu)勢地位，濟(jì)南市吸引了周邊區(qū)域的人才、資本、技術(shù)向其聚集，“馬太效應(yīng)”明顯，影響了相鄰區(qū)域的發(fā)展。

3 結(jié)論與建議

本研究基于2015年山東省的相關(guān)數(shù)據(jù)，運(yùn)用ESDA分析探討人口、土地、經(jīng)濟(jì)與社會城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間差異，結(jié)果表明，山東省大多數(shù)地市的城鎮(zhèn)化都處于中度及以上協(xié)調(diào)發(fā)展?fàn)顟B(tài);山東省人口、土地、經(jīng)濟(jì)和社會城鎮(zhèn)化協(xié)調(diào)發(fā)展度存在顯著的正空間自相關(guān)特點(diǎn);山東省城鎮(zhèn)化協(xié)調(diào)發(fā)展的均質(zhì)性區(qū)域明顯多于異質(zhì)性區(qū)域，城鎮(zhèn)化協(xié)調(diào)發(fā)展的空間依賴性和差異性同時(shí)存在，東部沿海和西部內(nèi)陸不均衡現(xiàn)象突出。H-H聚集類型區(qū)內(nèi)有7個(gè)市，主要集中在山東東部沿海和魯中地區(qū)，是城鎮(zhèn)化協(xié)調(diào)發(fā)展的“熱點(diǎn)”地區(qū);H-L聚集類型區(qū)域內(nèi)只有濟(jì)南市;L-L聚集類型區(qū)內(nèi)有5個(gè)市，主要分布在魯西和魯南地區(qū)，是城鎮(zhèn)化協(xié)調(diào)發(fā)展的“冷點(diǎn)”地區(qū);L-H聚集類型區(qū)內(nèi)只有4個(gè)市，且其顯著性不明顯，零星散布在魯西和魯南地區(qū)。

綜上，提出如下政策建議：第一，踐行科學(xué)的城鎮(zhèn)化理念，包括協(xié)調(diào)發(fā)展理念、共享理念、順應(yīng)市場規(guī)律理念。第二，實(shí)施差異化城鎮(zhèn)化協(xié)調(diào)發(fā)展戰(zhàn)略，提升山東半島城市群各市的城鎮(zhèn)化協(xié)調(diào)發(fā)展能力;加強(qiáng)宏觀調(diào)控，扶持魯西、魯南各市城鎮(zhèn)化協(xié)調(diào)發(fā)展。第三，創(chuàng)新保障城鎮(zhèn)化協(xié)調(diào)發(fā)展的體制機(jī)制，如戶籍制度、城鄉(xiāng)一體化的社會保障體系、進(jìn)城農(nóng)民就業(yè)培訓(xùn)體系、集體建設(shè)用地入市制度等。

參考文獻(xiàn)：

[1]林愛文，樊 星. 湖北省人口城鎮(zhèn)化與土地城鎮(zhèn)化協(xié)調(diào)發(fā)展分析[J]. 地域研究與開發(fā)，2015，34（6）：14-18.

[2]姚士謀，陸大道，陳振光，等. 順應(yīng)我國國情條件的城鎮(zhèn)化問題的嚴(yán)峻思考[J]. 經(jīng)濟(jì)地理，2012，32（5）：1-6.

[3]Diebold W，Despard L E.The evolution of the international economic order[J]. International Economic Policies and Their Theoretical Foundations，1992，9（4）：3-25.

[4]The price effects of international airline alliances[J]. Journal of Law and Economics，2000，43（2）：503-546.

[5]Henderson V.The urbanization process and economic growth：the so-what question[J]. Journal of Economic Growth，2003，8（1）：47-71.

[6]田時(shí)中，方眉玉，李光龍. 人口-土地-經(jīng)濟(jì)-社會城鎮(zhèn)化水平測度及耦合協(xié)調(diào)分析——基于31個(gè)省份面板數(shù)據(jù)的實(shí)證[J]. 河北地質(zhì)大學(xué)學(xué)報(bào)，2017，40（4）：49-57.

[7]多淑杰. 城鎮(zhèn)化四維分析模型的構(gòu)建與協(xié)調(diào)發(fā)展評價(jià)[J]. 企業(yè)經(jīng)濟(jì)，2014（8）：126-130.

[8]劉 勇，田 蕾，金 浩，等. 中國四維城鎮(zhèn)化協(xié)調(diào)性空間格局及演化研究[J]. 天津大學(xué)學(xué)報(bào)（社會科學(xué)版），2015，17（6）：513-517.

[9]張 目，黃春燕，李 巖. 基于相對熵和可變模糊集理論的戰(zhàn)略性新興產(chǎn)業(yè)企業(yè)信用評價(jià)[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識，2014，44（13）：18-26.

[10]遲國泰，齊 菲，張 楠. 基于最優(yōu)組合賦權(quán)的城市生態(tài)評價(jià)模型及應(yīng)用[J]. 運(yùn)籌與管理，2012，21（2）：183-191.

[11]劉耀彬，李仁東，宋學(xué)鋒. 中國城市化與生態(tài)環(huán)境耦合度分析[J]. 自然資源學(xué)報(bào)，2005，20（1）：105-112.

[12]Moran P A P.Notes on continuous stochastic phenomena[J]. Biometrika，1950，37（1/2）：17-23.