基于ESDA-GIS山東省人口與經(jīng)濟(jì)空間集聚差異特征研究

劉廣亮,吳宗杰,董會(huì)忠,萬(wàn)里洋

(山東理工大學(xué) 商學(xué)院，山東 淄博 255012)

基于ESDA-GIS山東省人口與經(jīng)濟(jì)空間集聚差異特征研究

劉廣亮,吳宗杰,董會(huì)忠,萬(wàn)里洋

(山東理工大學(xué) 商學(xué)院，山東 淄博 255012)

運(yùn)用ESDA和自相關(guān)指數(shù)方法，結(jié)合山東省十七城市2004-2013年人口和GDP數(shù)據(jù)，探索山東人口和經(jīng)濟(jì)空間分布特征動(dòng)態(tài)演化規(guī)律，分析人口、經(jīng)濟(jì)空間集聚差異的原因.結(jié)果表明：山東省十七地市人口和經(jīng)濟(jì)空間集聚特征顯著，人口空間分布和經(jīng)濟(jì)發(fā)展水平在空間上呈現(xiàn)差異性，說(shuō)明山東各區(qū)域整體經(jīng)濟(jì)發(fā)展水平差異加劇，各區(qū)域之間的經(jīng)濟(jì)關(guān)聯(lián)度減弱，山東沿海區(qū)域與內(nèi)陸區(qū)域經(jīng)濟(jì)發(fā)展不均衡，內(nèi)陸區(qū)域人口和經(jīng)濟(jì)發(fā)展不同步.山東省應(yīng)該采取適當(dāng)措施，縮小經(jīng)濟(jì)空間分布差距，實(shí)現(xiàn)區(qū)域經(jīng)濟(jì)均衡發(fā)展.

人口與經(jīng)濟(jì)空間分布；空間集聚；空間自相關(guān)

經(jīng)過(guò)30年粗放式發(fā)展，我國(guó)經(jīng)濟(jì)逐步進(jìn)入產(chǎn)業(yè)結(jié)構(gòu)調(diào)整階段和區(qū)域發(fā)展的轉(zhuǎn)型階段[1].人口和經(jīng)濟(jì)是區(qū)域發(fā)展差異性的重要體現(xiàn)因素[2].人口與經(jīng)濟(jì)相互影響，區(qū)域經(jīng)濟(jì)發(fā)展不平衡引起人口流動(dòng)，人口流動(dòng)又促使區(qū)域經(jīng)濟(jì)發(fā)展不平衡[3].了解和掌握區(qū)域人口和經(jīng)濟(jì)空間分布規(guī)律，對(duì)于實(shí)現(xiàn)區(qū)域經(jīng)濟(jì)空間規(guī)劃和區(qū)域發(fā)展有指導(dǎo)意義.一些學(xué)者從基尼系數(shù)(Gini)、泰爾指數(shù)(Theil)等角度研究[4-5]涉及人口變化[6]、區(qū)域經(jīng)濟(jì)[7]等.這些研究方法強(qiáng)調(diào)整個(gè)研究區(qū)域的特征值分布差異，忽視了區(qū)域內(nèi)部的特征值空間集聚[8-9].探索性空間數(shù)據(jù)分析法ESDA(Exploratory Spatial Data Analysis) 既分析整個(gè)研究區(qū)域特征值的分布特征，又強(qiáng)調(diào)區(qū)域內(nèi)部特征值的空間分布特征.

20世紀(jì)中后期，Moran、Geary、Anselin和Getis等學(xué)者對(duì)ESDA的發(fā)展起到了重要推動(dòng)作用，使ESDA理論和方法逐步完善[10-13].隨著空間統(tǒng)計(jì)、分析軟件的更新，運(yùn)用ESDA方法研究區(qū)域經(jīng)濟(jì)空間分布的成果越來(lái)越多.如運(yùn)用ESDA研究區(qū)域經(jīng)濟(jì)特征和空間環(huán)境的關(guān)系，揭示區(qū)域經(jīng)濟(jì)、人口分布和演化規(guī)律[14].以增量分析、不均衡指數(shù)等統(tǒng)計(jì)分析方法為基礎(chǔ)，從地市級(jí)尺度角度討論區(qū)域人口與經(jīng)濟(jì)的空間集聚特征與均衡過(guò)程[15].本文運(yùn)用ArcGIS10.2及 GeoDa 軟件研究山東省17地市區(qū)域經(jīng)濟(jì)和人口分布特征，研究區(qū)域演化規(guī)律，為山東省均衡發(fā)展和可持續(xù)發(fā)展政策的制定提供理論依據(jù).

1 研究方法

ESDA方法是展示區(qū)域?qū)傩詳?shù)據(jù)的空間分布狀況，描述空間關(guān)聯(lián)度、空間集聚模式，發(fā)現(xiàn)空間差異的一系列技術(shù)方法[16].空間自相關(guān)經(jīng)常使用的測(cè)度指標(biāo)是Moran ′s I.本文運(yùn)用全局Global Moran ′s I指數(shù)測(cè)度山東省17地市人口和GDP的空間關(guān)聯(lián)結(jié)構(gòu)，利用Moran散點(diǎn)圖與空間關(guān)聯(lián)局部性指標(biāo)(LISA)檢測(cè)局部空間的聚集性及分析局部空間的演化規(guī)律.

1.1 全局Moran ′s I

全局Moran ′s I是一種測(cè)度總體各觀測(cè)值之間自相關(guān)水平的綜合指標(biāo)，反映了區(qū)域與周邊區(qū)域的差異程度，計(jì)算如公式

(1)

1.2 Moran散點(diǎn)圖

Moran散點(diǎn)圖描述了變量Z與其空間滯后向量wz之間的相關(guān)關(guān)系[16].散點(diǎn)圖的橫軸為變量Z的觀測(cè)值，縱軸為空間滯后向量的值，形成四個(gè)象限.第一象限代表觀測(cè)區(qū)域與周邊區(qū)域的屬性值都比較高，即HH型(high-high)；第二象限代表觀測(cè)區(qū)域?qū)傩灾递^低，周邊區(qū)域?qū)傩灾递^高，即LH型(low-high)；第三象限代表觀測(cè)區(qū)域與周邊區(qū)域?qū)傩灾刀急容^低，即LL型(low-low)；第四象限代表觀測(cè)區(qū)域?qū)傩灾递^高，周邊區(qū)域?qū)傩灾递^低，即HL型(high-low).

1.3 局部Moran ′s I

Anselin在1995年提出局部Moran ′s I的概念，運(yùn)用局部Moran ′s I散點(diǎn)圖可以將區(qū)域空間分布的差異化呈現(xiàn)出來(lái)[16].其定義為

(2)

可寫(xiě)為

(3)

在給定置信水平下，若Ii和Zi都大于0，則區(qū)域i位于HH象限；若Ii大于0且Zi小于0，則區(qū)域i位于LL象限；若Ii小于0且Zi大于0，則區(qū)域i位于HL象限；若Ii和Zi都小于0，則區(qū)域i位于LH象限.

2 實(shí)證研究

2.1 研究區(qū)域概況和數(shù)據(jù)來(lái)源

山東省位居黃河下游，東靠黃海、渤海，北望朝鮮、日本，西北接連河北，西南毗鄰河南，南交安徽、江蘇.山東半島相對(duì)遼東半島，環(huán)抱渤海灣.獨(dú)特的區(qū)域位置，使山東省成為環(huán)渤海經(jīng)濟(jì)區(qū)與沿黃河經(jīng)濟(jì)帶的交匯點(diǎn)，華東和華北區(qū)域的交匯處，山東經(jīng)濟(jì)在我國(guó)經(jīng)濟(jì)中占有重要地位.山東經(jīng)濟(jì)發(fā)展迅速，2007年以后山東經(jīng)濟(jì)總量位居全國(guó)第3位，農(nóng)業(yè)增加值居全國(guó)首位，工業(yè)總產(chǎn)值及增加值位居全國(guó)前三位.

本文所采用的數(shù)據(jù)為山東省地理空間分布Shapefile格式數(shù)據(jù)，主要包括山東省17地市2004-2013年的經(jīng)濟(jì)發(fā)展數(shù)據(jù)(GDP)和人口數(shù)據(jù)，相關(guān)數(shù)據(jù)來(lái)自《山東省統(tǒng)計(jì)年鑒》和17地市統(tǒng)計(jì)年鑒.

2.2 人口和經(jīng)濟(jì)集聚差異格局分析

根據(jù)Global Moran ′s I的公式，將山東省17地市2004年到2013年的數(shù)據(jù)輸入Geoda軟件進(jìn)行全局自相關(guān)分析，得出山東省17地市人口與經(jīng)濟(jì)空間聚集的整體特征及變化(見(jiàn)表1).在正態(tài)假設(shè)條件下，全局Moran ′s I檢驗(yàn)系數(shù)全部大于0.005顯著水平下的臨界值1.96.表1表明，山東省17地市人口和經(jīng)濟(jì)整體呈現(xiàn)空間分布差異特征.2004-2013年Global Moran ′s I指數(shù)中2004年的人口指數(shù)為正數(shù)，其余指數(shù)均為負(fù)數(shù)，表明山東省人口分布具有顯著差異；人口指數(shù)2004年為0.034，其余年份指數(shù)整體波動(dòng)在(-0.0106至-0.0041)范圍內(nèi)，說(shuō)明山東省各地市人口分布與周邊地市有差異，差異性不明顯；山東省GDP指數(shù)中2004、2006、2007年為正數(shù)，其余年份指數(shù)波動(dòng)在-0.0048至-0.0002范圍內(nèi)，說(shuō)明山東省各地市GDP與周邊地市存在差異，差異性不明顯；此外，數(shù)據(jù)說(shuō)明各市人口和GDP存在顯著空間正相關(guān)，各市人口分布和經(jīng)濟(jì)發(fā)展水平表現(xiàn)出相似的空間差異，即人口密度較高(較低)的區(qū)域周邊區(qū)域人口密度相對(duì)低(高)、經(jīng)濟(jì)發(fā)展水平較高(低)的區(qū)域周邊經(jīng)濟(jì)發(fā)展水平較低(高).

表1 山東省2004-2013人口經(jīng)濟(jì)Moran’sI指數(shù)

指標(biāo)2004年2005年2006年2007年2008年2009年2010年2011年2012年2013年人口0．0340-0．0041-0．0067-0．0062-0．0069-0．0063-0．0081-0．0094-0．0106-0．0104GDP0．2427-0．00120．01140．0175-0．0038-0．0002-0．0002-0．0048-0．0005-0．0047

2004-2013年山東17地市區(qū)域人口Global Moran ′s I指數(shù)呈現(xiàn)先下降后平穩(wěn)下降的趨勢(shì).2004-2005年菏澤、臨沂、濟(jì)寧人口下降明顯.2005年沿海區(qū)域經(jīng)濟(jì)發(fā)展的速度高于內(nèi)陸區(qū)域發(fā)展速度，沿海區(qū)域相對(duì)于內(nèi)陸區(qū)域提供較多的就業(yè)機(jī)會(huì)，存在著內(nèi)陸人口向沿海區(qū)域遷徙的現(xiàn)象，導(dǎo)致山東省17地市區(qū)域人口分布空間分布存在差異性.2005年以后臨沂、濟(jì)寧經(jīng)濟(jì)較快，人口遷移不明顯，山東省人口分布差異性平穩(wěn).2004-2013年山東省GDPGlobal Moran ′s I指數(shù)整體呈先下降后平穩(wěn)下降的趨勢(shì)，說(shuō)明2005年以后山東省十七地市區(qū)域經(jīng)濟(jì)發(fā)展明顯不平衡，區(qū)域經(jīng)濟(jì)水平分布差異性加劇，各區(qū)域之間的經(jīng)濟(jì)關(guān)聯(lián)度減弱.現(xiàn)實(shí)原因：沿海區(qū)域煙臺(tái)、青島、濰坊區(qū)域經(jīng)濟(jì)發(fā)展速度比日照、威海區(qū)域高；內(nèi)陸區(qū)域濟(jì)南、濟(jì)寧、臨沂經(jīng)濟(jì)發(fā)展速度比周邊區(qū)域高.沿海區(qū)域與內(nèi)陸區(qū)域經(jīng)濟(jì)關(guān)聯(lián)度較低，內(nèi)陸區(qū)域濟(jì)南、臨沂、濟(jì)寧經(jīng)濟(jì)與周邊區(qū)域經(jīng)濟(jì)關(guān)聯(lián)度較低.

2.3 人口空間集聚差異特征分析

利用ArcGis10.2軟件描繪出2004年，2008年和2013年3個(gè)年份為代表的如圖1所示.

從圖1可看出，山東省17地市人口空間集聚特征顯著，人口空間集聚區(qū)域具有明顯差異：山東省地17市區(qū)域人口集聚的高高型(HH)區(qū)域主要集中分布在魯西南、沿海城市.魯西南區(qū)域：菏澤、臨沂、濟(jì)寧、泰安、聊城5地市土地面積接近山東省土地面積1/3，人口超過(guò)全省1/3，魯西南區(qū)域地域遼闊，農(nóng)、牧業(yè)歷史悠久且穩(wěn)步發(fā)展，古代就是人口稠密區(qū).沿海區(qū)域：濰坊、青島、煙臺(tái)地處環(huán)渤海經(jīng)濟(jì)圈，區(qū)位優(yōu)勢(shì)明顯，交通網(wǎng)絡(luò)四通八達(dá).此外，濰坊區(qū)域自然資

(a)2004年 (b)2008年 (c)2013年

圖1 2004、2008、2013年山東省人口LISA集聚圖

源比周邊區(qū)域豐富，環(huán)境優(yōu)美，是我國(guó)蔬菜基地，經(jīng)濟(jì)較發(fā)達(dá).濟(jì)南、淄博、東營(yíng)、萊蕪、棗莊、日照、威海7市，以高低型(HL)和低高型(LH)人口集聚為主.萊蕪地區(qū)多丘陵，鋼鐵企業(yè)對(duì)該區(qū)域經(jīng)濟(jì)發(fā)展起到?jīng)Q定性作用；棗莊是資源型城市，棗莊區(qū)域人口集聚主要集中在建國(guó)后；東營(yíng)，位居山東北，石油資源豐富，多鹽堿地，人口稠密度較低.濱州市人口分布以低低型(LL)為主、濱州地處山東正北，土地以鹽堿地居多，經(jīng)濟(jì)相對(duì)落后、人口密度較低.研究發(fā)現(xiàn)德州2004年人口數(shù)據(jù)位于LISA散點(diǎn)圖的數(shù)軸上，人口集聚模式介于高低型(HL)和高高型(HH)，將其標(biāo)注為過(guò)渡類(lèi)型區(qū)域.2013年德州人口集聚模式從過(guò)渡型發(fā)展成高高型(HH).

2004-2013年山東省17地市人口集聚特征整體變化不明顯.2004年山東省17地市地區(qū)人口集聚高高型、高低型、低高型、低低型和過(guò)渡型的區(qū)域分別為8、1、6、1和1個(gè)，到2013年該區(qū)域人口高高型、高低型、低高型、低低型和過(guò)渡型區(qū)域分別為9、1、6、1和0個(gè)，在這10年間德州從過(guò)渡型城市轉(zhuǎn)變成人口集聚高高型區(qū)域，其他區(qū)域沒(méi)有發(fā)生變化.

2.4 經(jīng)濟(jì)空間集聚差異特征分析

以2004年，2008年和2013年3個(gè)年份為代表分析山東省GDP空間集聚特征如圖2所示.

(a)2004年 (b)2008年 (c)2013年

圖2 2004、2008、2013年山東省GDP LISA集聚圖

從圖2看出，2013年山東省經(jīng)濟(jì)區(qū)域分布圖與2004年比較有明顯變化.沿海區(qū)域經(jīng)濟(jì)整體以HH型集聚發(fā)展；濟(jì)南、臨沂、濟(jì)寧區(qū)域經(jīng)濟(jì)是HL型集聚；部分內(nèi)陸、沿海區(qū)域經(jīng)濟(jì)是LL型集聚.沿海區(qū)域：青島、濰坊、煙臺(tái)、東營(yíng)經(jīng)濟(jì)以明顯高于周邊區(qū)域的速度發(fā)展(從2004年到2013年GDP分別增長(zhǎng)了2.53、2.69、2.44、2.64倍)，同為沿海區(qū)域的威海經(jīng)濟(jì)發(fā)展速度比其他沿海區(qū)域慢(威海增長(zhǎng)了1.66倍)，沿海區(qū)域的濱州、日照經(jīng)濟(jì)起步較晚，經(jīng)濟(jì)發(fā)展水平遠(yuǎn)低于其他沿海區(qū)域，沿海區(qū)域經(jīng)濟(jì)發(fā)展水平在空間分布上形成差異；內(nèi)陸地區(qū)：濟(jì)南的省會(huì)地位是濟(jì)南經(jīng)濟(jì)高速發(fā)展的區(qū)域優(yōu)勢(shì)，臨沂的物流基地促使經(jīng)濟(jì)飛速崛起，濟(jì)寧經(jīng)濟(jì)憑著交通、運(yùn)河、文化旅游等優(yōu)勢(shì)逐步提升，這三個(gè)區(qū)域與周邊內(nèi)陸區(qū)域經(jīng)濟(jì)產(chǎn)生差距，內(nèi)陸區(qū)域經(jīng)濟(jì)空間分布差異呈擴(kuò)大趨勢(shì).研究表明，各區(qū)域經(jīng)濟(jì)發(fā)展水平差異加劇，各區(qū)域之間的經(jīng)濟(jì)關(guān)聯(lián)度呈減弱趨勢(shì)，區(qū)域經(jīng)濟(jì)發(fā)展不平衡.沿海區(qū)域經(jīng)濟(jì)發(fā)展水平整體高于內(nèi)陸經(jīng)濟(jì)發(fā)展水平；沿海區(qū)域中青島、煙臺(tái)、濰坊區(qū)域經(jīng)濟(jì)發(fā)展水平高于周邊區(qū)域；內(nèi)陸區(qū)域中濟(jì)南、濟(jì)寧、臨沂區(qū)域經(jīng)濟(jì)發(fā)展水平高于周邊區(qū)域.此外，研究發(fā)現(xiàn)2004年淄博位于散點(diǎn)圖的數(shù)軸上.主要原因：位于淄博東北的濰坊、東營(yíng)經(jīng)濟(jì)較發(fā)達(dá)；位于其西北處的濱州、西南處的萊蕪，臨沂經(jīng)濟(jì)較落后.到2013年淄博周邊的區(qū)域中：濟(jì)南、臨沂、東營(yíng)、濰坊經(jīng)濟(jì)發(fā)展水平較高，濱州、萊蕪相對(duì)落后，淄博經(jīng)濟(jì)發(fā)展水平由介于HH型和HL型集聚區(qū)域發(fā)展成HH型集聚區(qū)域.

2004-2013年山東省17地市經(jīng)濟(jì)發(fā)展水平在空間分布上呈現(xiàn)差異性，經(jīng)濟(jì)發(fā)展存在明顯差異.2004年HH型6個(gè)，HL型和LH型為3個(gè)，LL型7個(gè)，過(guò)渡型1個(gè).2013年HH型6個(gè)，HL型和LH型為11個(gè)，LL型0個(gè)，過(guò)渡型0個(gè).沿海地區(qū)經(jīng)濟(jì)整體以較快的速度發(fā)展，濱州、日照和威海地區(qū)經(jīng)濟(jì)發(fā)展相對(duì)較慢；山東省內(nèi)陸區(qū)域經(jīng)濟(jì)整體發(fā)展較慢，濟(jì)南、濟(jì)寧和臨沂地區(qū)經(jīng)濟(jì)發(fā)展相對(duì)較快.

2.5 人口與經(jīng)濟(jì)空間集聚差異分析

2.5.1 人口與經(jīng)濟(jì)空間集聚差異類(lèi)型

結(jié)合山東省2004年、2008年和2013年3個(gè)截面人口和GDP LISA集聚圖，對(duì)山東省各地區(qū)人口和經(jīng)濟(jì)發(fā)展?fàn)顩r進(jìn)行定性分類(lèi).依據(jù)山東省17地市人口Moran散點(diǎn)圖，將位于第一、第四象限的集聚模式的區(qū)域劃分人口密度較高的區(qū)域，將位于第二、第三象限的集聚模式的區(qū)域劃分為人口密度較低的區(qū)域；以此類(lèi)推，根據(jù)山東省17地市GDPMoran 散點(diǎn)圖，將位于第一、第四象限的集聚模式的區(qū)域劃分為經(jīng)濟(jì)發(fā)展水平較高區(qū)域，將位于第二、第三象限的集聚模式的區(qū)域劃分為經(jīng)濟(jì)發(fā)展水平較低區(qū)域.依據(jù)高低重組，形成新的四個(gè)象限，分別代表四類(lèi)地區(qū).劃分原則：以協(xié)調(diào)發(fā)展為主要參考因素，以區(qū)域轉(zhuǎn)型的可能性為次要因素(該區(qū)域10年人口分布變化不明顯，經(jīng)濟(jì)變化較大，所以經(jīng)濟(jì)變動(dòng)作為次要因素).因此，第一象限代表一類(lèi)地區(qū)；第二象限代表四類(lèi)地區(qū)；第三象限代表二類(lèi)地區(qū)；第四象限代表三類(lèi)地區(qū).將介于兩象限之間(即在數(shù)軸上或緊鄰數(shù)軸的)劃為過(guò)渡區(qū).具體特征解釋見(jiàn)表2和圖3).

從山東省人口和經(jīng)濟(jì)空間集聚類(lèi)型分析得出山東省人口分布和經(jīng)濟(jì)發(fā)展的區(qū)域特征：一類(lèi)地區(qū)基本分布在山東半島藍(lán)色經(jīng)濟(jì)區(qū)的中部區(qū)域，包括煙臺(tái)、青島、濰坊三個(gè)區(qū)域.該區(qū)域具有優(yōu)越的地理區(qū)位條件、具有廣闊的海域空間、相對(duì)豐富的海洋資源和深厚的海洋經(jīng)濟(jì)與科技，為工業(yè)生產(chǎn)提供了的基礎(chǔ)；此外，該區(qū)域地形以平原為主，土地肥沃、人口密集、交通便利，為蔬菜產(chǎn)業(yè)提供了自然基礎(chǔ)，人口與經(jīng)濟(jì)同步發(fā)展，經(jīng)濟(jì)發(fā)展帶動(dòng)人口增加，人口增加推動(dòng)經(jīng)濟(jì)發(fā)展.二類(lèi)地區(qū)包括濱州、萊蕪、棗莊、日照、威海區(qū)域.濱州、日照、威海為沿海地區(qū)，濱州土地為鹽堿地、日照、威海多山嶺；萊蕪為內(nèi)陸地區(qū)，該區(qū)域多丘陵，山陵，經(jīng)濟(jì)依托鋼鐵企業(yè)；棗莊為資源型城市，依托當(dāng)?shù)孛禾堪l(fā)展經(jīng)濟(jì)，人口相對(duì)稀疏，隨著資源枯竭，經(jīng)濟(jì)發(fā)展動(dòng)力不足.三類(lèi)地區(qū)主要集中在西部經(jīng)濟(jì)隆起帶區(qū)域.魯西南區(qū)域人口稠密，經(jīng)濟(jì)發(fā)展相對(duì)于沿海區(qū)域較慢.四類(lèi)地區(qū)包括東營(yíng)、淄博.淄博早期工業(yè)經(jīng)濟(jì)起步較早，東營(yíng)經(jīng)濟(jì)依托豐富的石油資源，人口發(fā)展滯后于經(jīng)濟(jì)增長(zhǎng).

表2 山東省人口與經(jīng)濟(jì)集聚類(lèi)型劃分

類(lèi)別基本特征人口集聚模式經(jīng)濟(jì)集聚模式2004年2008年2013年一類(lèi)地區(qū)區(qū)域GDP較高，人口密度較高，人口經(jīng)濟(jì)同步發(fā)展HHHLHHHL366二類(lèi)地區(qū)區(qū)域GDP較低，人口密度較低，人口經(jīng)濟(jì)同步發(fā)展LLLHLLLH255三類(lèi)地區(qū)區(qū)域GDP較低，人口密度較高，人口經(jīng)濟(jì)發(fā)展不同步HHHLLLLH644四類(lèi)地區(qū)區(qū)域GDP較高，人口密度較低，人口經(jīng)濟(jì)發(fā)展不同步LLLHHHHL522過(guò)渡區(qū)介于一類(lèi)與三類(lèi)，一類(lèi)與四類(lèi)，二類(lèi)與三類(lèi)、二類(lèi)與四類(lèi)的區(qū)域L-HL-H100

從2004到2013年，一類(lèi)地區(qū)在原有區(qū)域的基礎(chǔ)上新增了濟(jì)南、濟(jì)寧、臨沂三個(gè)區(qū)域.濟(jì)南、濟(jì)寧、臨沂自2004年依托各自?xún)?yōu)勢(shì)經(jīng)濟(jì)發(fā)展速度相對(duì)周邊區(qū)域較快，由三類(lèi)地區(qū)轉(zhuǎn)為一類(lèi)地區(qū).過(guò)渡型區(qū)域德州轉(zhuǎn)變成三類(lèi)地區(qū)，該區(qū)域人口稠密，經(jīng)濟(jì)發(fā)展較慢.山東省人口、經(jīng)濟(jì)分布格局是：藍(lán)色半島經(jīng)濟(jì)區(qū)中心和內(nèi)陸特色區(qū)域?yàn)楦吒呒鄣囊活?lèi)區(qū)域；以土地貧瘠的區(qū)域?yàn)榈偷图鄣亩?lèi)區(qū)域；以位于魯西、魯南區(qū)域?yàn)楦叩图鄣娜?lèi)區(qū)域；以經(jīng)濟(jì)較發(fā)達(dá)人口較少的東營(yíng)、淄博為低高集聚的四類(lèi)區(qū)域.

(a)2004年 (b)2008年 (c)2013年

圖3 山東省人口與經(jīng)濟(jì)空間集聚類(lèi)型劃分

2.5.2 人口和經(jīng)濟(jì)空間集聚差異改善

人口與經(jīng)濟(jì)協(xié)調(diào)發(fā)展是經(jīng)濟(jì)全球化的趨勢(shì)，是區(qū)域轉(zhuǎn)型的內(nèi)在要求.提升區(qū)域人口和區(qū)域經(jīng)濟(jì)的空間關(guān)聯(lián)度，減緩分布差異，為區(qū)域經(jīng)濟(jì)持久協(xié)調(diào)發(fā)展提供內(nèi)在基礎(chǔ).山東省沿海區(qū)域經(jīng)濟(jì)發(fā)達(dá)、人口密集，沿海區(qū)域人口與經(jīng)濟(jì)的協(xié)調(diào)發(fā)展對(duì)山東省發(fā)展有較強(qiáng)的推動(dòng)作用.實(shí)現(xiàn)沿海區(qū)域與內(nèi)陸區(qū)域協(xié)調(diào)同步發(fā)展，人口和經(jīng)濟(jì)空間集聚差異需要改善.

一類(lèi)地區(qū)：在保持區(qū)域人口和經(jīng)濟(jì)的空間HH集聚的發(fā)展趨勢(shì)的基礎(chǔ)上，充分利用經(jīng)濟(jì)綜合實(shí)力，加強(qiáng)該區(qū)域基礎(chǔ)建設(shè)，完善城市保障體系，以該類(lèi)區(qū)域?yàn)橹行?，帶?dòng)周邊區(qū)域經(jīng)濟(jì)發(fā)展.山東17地市人口與經(jīng)濟(jì)高高集聚區(qū)域近1/3，區(qū)域數(shù)量以及區(qū)域面積相對(duì)穩(wěn)定.以沿海區(qū)域?yàn)槲沂〗?jīng)濟(jì)發(fā)展的推動(dòng)力，以?xún)?nèi)陸區(qū)域的一類(lèi)區(qū)域?yàn)槲沂〗?jīng)濟(jì)發(fā)展的帶動(dòng)力，共同推動(dòng)周邊非均衡區(qū)域的協(xié)調(diào)發(fā)展.

二類(lèi)地區(qū)：二類(lèi)區(qū)域較多，分布相對(duì)分散，人口密度與經(jīng)濟(jì)的發(fā)展主要受到該區(qū)域的自然條件限制.該區(qū)農(nóng)業(yè)落后，受到山嶺、鹽堿地等區(qū)域問(wèn)題約束；人口密低，產(chǎn)業(yè)結(jié)構(gòu)單一，經(jīng)濟(jì)落后等是該區(qū)域的基本特征.該區(qū)域是我省經(jīng)濟(jì)全面快速協(xié)調(diào)發(fā)展的限制區(qū)域，是我省人口經(jīng)濟(jì)協(xié)調(diào)發(fā)展的重要區(qū)域.該區(qū)域需要規(guī)劃與實(shí)現(xiàn)人口經(jīng)濟(jì)分布格局，依據(jù)地區(qū)特色，發(fā)展區(qū)域特色經(jīng)濟(jì)(例濱州、日照可利用港口特色，發(fā)展運(yùn)輸經(jīng)濟(jì))；依托區(qū)域資源，發(fā)展資源經(jīng)濟(jì)(威?？梢劳泻０讹L(fēng)景，發(fā)展旅游業(yè))；依照區(qū)域土地資源，發(fā)展特色農(nóng)業(yè)(棗莊、萊蕪可利用其山嶺，發(fā)展特色農(nóng)業(yè)).

三類(lèi)地區(qū)：該區(qū)域基本分布在我省西部經(jīng)濟(jì)隆起帶上，該區(qū)域人口密度高而經(jīng)濟(jì)發(fā)展較慢.該區(qū)域缺少主導(dǎo)產(chǎn)業(yè)，經(jīng)濟(jì)起步相對(duì)沿海區(qū)域較晚.改善該區(qū)域現(xiàn)狀首先要引進(jìn)資金、技術(shù)，發(fā)展該區(qū)域其他產(chǎn)業(yè)；其次，依托科技興農(nóng)，轉(zhuǎn)變農(nóng)業(yè)發(fā)展模式，實(shí)現(xiàn)農(nóng)業(yè)規(guī)?；蛔詈?，推進(jìn)新農(nóng)村和城鎮(zhèn)化建設(shè)，引導(dǎo)農(nóng)村剩余勞動(dòng)力轉(zhuǎn)移；從而推動(dòng)該區(qū)域經(jīng)濟(jì)快速發(fā)展.

四類(lèi)地區(qū)：該區(qū)域包括東營(yíng)、淄博兩個(gè)區(qū)域.東營(yíng)、淄博經(jīng)濟(jì)比較發(fā)達(dá)，東營(yíng)土地鹽堿、沙化比較嚴(yán)重，淄博部分區(qū)域是山區(qū)，兩地人口密度相對(duì)周邊區(qū)域較低.該區(qū)域可以加強(qiáng)基礎(chǔ)建設(shè)，完善城市保障體系，提高城市居住環(huán)境，優(yōu)化城市布局，提供更多的就業(yè)崗位，吸引人才集聚.

3 結(jié)論與建議

運(yùn)用Moran ′s I指數(shù)，得出山東省17地市2004-2013年人口分布和GDP空間自相關(guān)指數(shù).說(shuō)明山東省人口與經(jīng)濟(jì)發(fā)展空間分布差異加大；通過(guò)山東省人口和經(jīng)濟(jì)集聚差異性分析得到：人口以HH型集聚的區(qū)域主要分布在山東半島藍(lán)色經(jīng)濟(jì)區(qū)中部以及魯南部分區(qū)域；人口以HL和LH型集聚的區(qū)域主要分布在魯北區(qū)域；經(jīng)濟(jì)以HH型集聚的區(qū)域主要分布在山東半島藍(lán)色經(jīng)濟(jì)區(qū)中部以及魯南部分區(qū)域.濟(jì)南、臨沂、濟(jì)寧區(qū)域經(jīng)濟(jì)為HL型集聚；其他區(qū)域經(jīng)濟(jì)為L(zhǎng)L型集聚.

少數(shù)區(qū)位優(yōu)勢(shì)明顯、資源豐富的區(qū)域經(jīng)濟(jì)發(fā)展水平較高，應(yīng)該以這些區(qū)域?yàn)榛A(chǔ)點(diǎn)帶動(dòng)整個(gè)區(qū)域的經(jīng)濟(jì)發(fā)展.一類(lèi)地區(qū)主要分布在東部沿海區(qū)域和部分內(nèi)陸區(qū)域，該區(qū)域人口與經(jīng)濟(jì)協(xié)調(diào)發(fā)展，呈高高集聚特征，應(yīng)推動(dòng)城市集群化發(fā)展，培育龍頭企業(yè)，提升城市群整體實(shí)力，帶動(dòng)山東人口與經(jīng)濟(jì)共同發(fā)展；二類(lèi)地區(qū)分布比較廣泛，該區(qū)域發(fā)展需要因地制宜，發(fā)展特色經(jīng)濟(jì)，發(fā)展資源經(jīng)濟(jì)，發(fā)展特色農(nóng)業(yè).吸引資金、技術(shù)發(fā)展符合生態(tài)環(huán)境條件的產(chǎn)業(yè)；三類(lèi)地區(qū)基本分布在魯西、魯南等區(qū)域，需要吸引技術(shù)、資金，發(fā)展區(qū)域經(jīng)濟(jì)；需要實(shí)現(xiàn)農(nóng)業(yè)效益規(guī)?；瑢?shí)現(xiàn)農(nóng)村城鎮(zhèn)化，促進(jìn)區(qū)域格局轉(zhuǎn)變化.四類(lèi)地區(qū)包括東營(yíng)、淄博，需要加強(qiáng)該區(qū)域基礎(chǔ)建設(shè)，提供就業(yè)機(jī)會(huì)，健全城市保障體系，吸引人口集聚，促進(jìn)區(qū)域轉(zhuǎn)型.

[1]劉慧, 樊杰, 王傳勝. 歐盟空間規(guī)劃研究進(jìn)展及啟示[J]. 地理研究, 2008, 27(6):1 381-1 389.

[2]王磊, 段學(xué)軍, 田方等.長(zhǎng)江三角洲人口與經(jīng)濟(jì)的空間分布關(guān)系研究[J].經(jīng)濟(jì)地理,2009, 29(10):1 619-1 623.

[3] 樊杰, 陶岸君, 呂晨. 中國(guó)經(jīng)濟(jì)與人口重心的耦合態(tài)勢(shì)及其對(duì)區(qū)域發(fā)展的影響[J]. 地理科學(xué)進(jìn)展，2010, 29(1):87-95.

[4]韓增林, 許旭.中國(guó)海洋經(jīng)濟(jì)地域差異及演化過(guò)程分析[J]. 地理研究, 2008, 27(3):613-622.

[5]蘆惠, 歐向軍, 李想，等.中國(guó)區(qū)域經(jīng)濟(jì)差異與極化的時(shí)空分析[J]. 經(jīng)濟(jì)地理, 2013, 33(6):15-21.

[6]ALAHMADIM,ATKINSONP,MARTIND.EstimatingthespatialdistributionofthepopulationofRiyadh,SaudiArabiausingremotelysensedbuiltlandcoverandheightdata[J].Computers,EnvironmentandUrbanSystems, 2013(41):167-176.

[7]BASILER,DURBANM,MINGUEZCR.Modelingregionaleconomicdynamics:Spatialdependence,spatialheterogeneityandnonlinearities[J].JournalofEconomicDynamicsandControl, 2014(48):229-245.

[8]UNWINA,UNWIND.Exploratoryspatialdataanalysiswithlocalstatistics[J].TheStatistician, 1998, 47(3):415-421.

[9]GIUSEPPEA.Theroleofspatialeffectsintheempiricalanalysisofregionalconcentration[J].Journalofgeographicalsystems, 2001, 3(3): 271-281.

[10]MORANPAP.Notesoncontinuousstochasticphenomena[J].Biometrika, 1950, 37(1-2): 17-23.

[11]GEARYRC.TheContiguityratioandstatisticalmapping[J].TheIncorporatedStatistician, 1954, 5(3): 115-145.

[12]ANSELINL.Localindicatorsofspatialassociation-LISA[J].GeographicalAnalysis, 1995, 27(2): 93-115.

[13]GETISA,ORDJK.Theanalysisofspatialassociationbyuseofdistancestatistics[J].GeographicalAnalysis, 1992, 24(3): 189-206.

[14]劉艷華,徐勇,劉毅.基于ESDA方法的黃土高原地區(qū)經(jīng)濟(jì)發(fā)展格局及其演化特征分析[J].地理科學(xué)進(jìn)展. 2011, 30(5):627-634.

[15]蔣子龍,樊杰,陳東.2001-2010年中國(guó)人口與經(jīng)濟(jì)的空間集聚與均衡特征分析.經(jīng)濟(jì)地理[J], 2014，34(5)：9-12.

[16]呂韜,曹有揮,“時(shí)空接近”空間自相關(guān)模型建立及其應(yīng)用——以長(zhǎng)三角區(qū)域經(jīng)濟(jì)差異分析為例[J]. 地理研究, 2010, 29(2):351-360.

(編輯：劉寶江)

Analysis of spatial agglomeration of population and economic in Shandong Province based on ESDA-GIS method

LIU Guang-liang，WU Zong-jie，DONG Hui-zhong，WAN Li-yang

(School of Business, Shandong University of Technology, Zibo 255012, China)

Based on the population and GDP data of 17 cities in Shandong Province during 2004 - 2013, we use the methods of ESDA gravity center and the index of geographic concentration, and regional spatial agglomeration differences，which explored the regional spatial distribution gather characteristics and regular pattern of dynamic evolution in Shandong Province, and analyze the reason of differences of economic and population spatial distribution. The results show that the spatial distributions of population and economy in Shandong Province display significant spatial agglomeration. The spatial distribution of population and economic presents the difference. This suggests that the level of economic development in Shandong regional differences is expanding, and the connections between every region′s economic become weakened. Studies show that the Shandong regional economic development is imbalance, that inland regional population and economic development are not coordinate. Shandong Province should take appropriate measures to narrow the gap between the spatial distributions of economy, and realize the balanced development of regional economy.

the spatial distribution of population and economic; spatial agglomeration; spatial auto-correlation

2016-05-31

國(guó)家自然科學(xué)基金項(xiàng)目(71371112)

劉廣亮，男，1115332867@ qq.com ； 通信作者： 吳宗杰，男，wzj@sdut.edu.cn

1672-6197(2017)03-0030-06

C924.2;F127

A