山東省縣域經(jīng)濟空間格局演變及趨勢分析

摘要：文章運用空間自相關(guān)、標(biāo)準(zhǔn)差橢圓及灰色預(yù)測模型等工具對2013-2022年山東省縣域經(jīng)濟空間格局演變及趨勢進行研究分析。研究結(jié)果表明，研究期內(nèi)山東省縣域經(jīng)濟發(fā)展空間聚集分布顯著，但是聚集程度總體上不斷減弱，縣域經(jīng)濟空間差異不斷縮小，且縮小趨勢由加速轉(zhuǎn)向平緩；山東省縣域經(jīng)濟局部空間關(guān)聯(lián)性強，極化效應(yīng)比較顯著，形成熱點區(qū)、冷點區(qū)連片集中聚集分布態(tài)勢，形成半島及藍黃經(jīng)濟區(qū)“兩翼一廊”經(jīng)濟隆起帶；縣域經(jīng)濟時空演變特征上，縣域經(jīng)濟重心向東北方向遷移，但遷移幅度有限，重心位置變化不大。空間分布范圍略微擴張，空間分布方向呈東北－西南走向，沿空間分布方向的縣域經(jīng)濟溢出效應(yīng)略為增強；未來山東省縣域經(jīng)濟重心向西北方向移動，空間分布主軸仍然保持在東北－西南方向，但是沿長軸方向的空間溢出效應(yīng)減弱，沿短軸方向空間溢出效應(yīng)增強，縣域經(jīng)濟差異性將進一步降低。

關(guān)鍵詞：縣域經(jīng)濟；空間格局；演變；預(yù)測

縣域經(jīng)濟是區(qū)域經(jīng)濟的一種特定形態(tài)和基本經(jīng)濟空間單元。縣域是區(qū)域經(jīng)濟空間格局研究的重要層面之一，較之省域、市域更加微觀化，其微觀視角有利于拓展區(qū)域經(jīng)濟研究。縣域經(jīng)濟空間格局是某一行政或者經(jīng)濟地理區(qū)域內(nèi)縣域單元經(jīng)濟活動在空間上的表現(xiàn)形式，主要反映域內(nèi)縣域經(jīng)濟要素空間分布形態(tài)、特征及域內(nèi)縣域經(jīng)濟要素相互作用的機理。一些研究者以山東省縣域經(jīng)濟空間格局為研究對象，使用探索性空間數(shù)據(jù)分析方法（ESDA）等空間計量模型對其演變特征進行分析，主要研究內(nèi)容為空間分布、聚類分布和時空格局變遷。劉玉等（2012）認(rèn)為山東縣域經(jīng)濟差異變化幅度較小，且縣域經(jīng)濟差異先擴大再縮小；經(jīng)濟重心向東北方向遷移；縣域發(fā)展空間自相關(guān)較為顯著，高高集聚區(qū)、低低集聚區(qū)分別集中分布于山東半島和魯西南、魯中南。白全民等（2013）對1987年、2000－2010年山東省縣域經(jīng)濟空間格局演變及其動態(tài)影響機制進行了研究，認(rèn)為山東省縣域經(jīng)濟空間自相關(guān)指數(shù)呈先降后升趨勢，縣域關(guān)聯(lián)性較強，且形成黃河南岸及“青煙威”半島區(qū)域兩大經(jīng)濟增長極。王慧等（2014）認(rèn)為1995-2012年山東省縣域經(jīng)濟絕對差異顯著擴大，相對差異緩慢上升，極化現(xiàn)象顯著；縣域經(jīng)濟空間自相關(guān)性顯著，空間自相關(guān)集聚模式較為穩(wěn)定。杜霞等（2015）認(rèn)為2000-2010年山東省縣域經(jīng)濟差異逐漸加大，高高型區(qū)域主要集中分布于山東東部沿海地區(qū)，低低型區(qū)域主要集中分布在魯西和魯南地區(qū)，兩類區(qū)域分別逐漸向魯北、魯中地區(qū)擴散；縣域經(jīng)濟發(fā)展水平呈高低交錯空間格局即北高南低、東高西低，極化效應(yīng)明顯；熱點區(qū)主要集中分布于以濟南、青島為核心的魯中和沿海地帶，以及北部黃河三角洲地帶，冷點區(qū)則主要集中分布于魯西、魯西北及魯中南地區(qū)。蒙莎莎等（2017）研究了2005－2015年山東省縣域經(jīng)濟空間格局演化特征和地域分異驅(qū)動機制，認(rèn)為山東省縣域經(jīng)濟總體格局呈東北—西南走向、重心向東遷移且速度加快，縣域經(jīng)濟空間差異有所增強，經(jīng)濟空間集聚程度增強、呈顯著核心—外圍二元結(jié)構(gòu)，縣域經(jīng)濟增長極與低洼區(qū)兩極分化嚴(yán)重。這些研究揭示了過去一段時期山東省縣域經(jīng)濟空間格局演變的共性特征：縣域經(jīng)濟呈空間聚集分布，且聚集程度增強；高高聚集區(qū)、低低聚集區(qū)及熱點區(qū)、冷點區(qū)分化；縣域經(jīng)濟空間格局時空演變特征明顯。但是，縣域經(jīng)濟聚集程度趨勢變化，聚類分布演變趨勢及時空格局經(jīng)濟重心遷移、空間方向等演變趨勢也得出了不同結(jié)論。

縣域經(jīng)濟發(fā)展的空間差異性、不均衡性，不斷推動縣域經(jīng)濟空間格局演變。縣域經(jīng)濟空間格局演變研究揭示了空間上縣域經(jīng)濟發(fā)展動態(tài)特征，是制定和調(diào)整區(qū)域發(fā)展政策的重要依據(jù)之一。“十三五”以來，山東省實施重點區(qū)域帶動戰(zhàn)略，不斷優(yōu)化迭代“兩區(qū)一圈一帶”“一群兩心三圈”區(qū)域發(fā)展總體布局，加快區(qū)域經(jīng)濟一體化，影響著縣域經(jīng)濟空間格局演變。為此，本文運用空間自相關(guān)、標(biāo)準(zhǔn)差橢圓等工具，研究2013－2022年山東省縣域經(jīng)濟空間分布模式、聚類模式和時空格局變遷，探索新時期山東省縣域經(jīng)濟空間格局新的演變特征。引入灰色預(yù)測模型，對未來一段時期山東省縣域經(jīng)濟發(fā)展格局趨勢進行了測度和分析。

一、研究方法及數(shù)據(jù)來源

（一）空間自相關(guān)

空間自相關(guān)分析是探索地理對象之間空間關(guān)聯(lián)性的重要工具，本文用以揭示山東省縣域經(jīng)濟空間分布格局及其差異特征。

1. 全局空間自相關(guān)。全局空間自相關(guān)（Global Moran’I）用以觀測整個研究區(qū)域內(nèi)觀測值的空間自相關(guān)性，描述某個地理區(qū)域內(nèi)觀測值是否彼此相關(guān)及這種相關(guān)性的強度和顯著性。如果全局 Moran’I顯著不為零，表明整個研究區(qū)域內(nèi)的觀測點在空間上不是隨機分布的。全局 Moran’I大于零，表示空間正相關(guān)性，即呈現(xiàn)空間聚集分布，其值越大，空間自相關(guān)性越強。

2. 局部空間自相關(guān)。局部空間自相關(guān)用以識別特定位置及其領(lǐng)域的空間自相關(guān)性，識別熱點、冷點或者高值聚集區(qū)、低值聚集區(qū)等局部模式，能夠反映特定位置與領(lǐng)域之間的空間差異程度和顯著性。熱點分析（Getis’s-Ord Gi*）是空間自相關(guān)分析常用工具之一，揭示局部要素空間分布相對于總體的模式特征，用以發(fā)掘數(shù)據(jù)特征在空間模式上是否存在任何統(tǒng)計顯著性聚類。與局部Moran’I定量分析每個空間單元與鄰域內(nèi)所有空間單元之間的相互聯(lián)系不同，Getis’s-Ord Gi*則定量分析這些關(guān)聯(lián)的具體程度。

（二）標(biāo)準(zhǔn)差橢圓

標(biāo)準(zhǔn)差橢圓（Standard Deviational Ellipse，SDE）是一種定量描述地理要素某屬性的中心趨勢、離散趨勢和方向趨勢等空間分布特征，探索地理要素空間分布狀態(tài)演變的方法。通過確定研究區(qū)域標(biāo)準(zhǔn)差橢圓的空間分布范圍和平均中心、長軸、短軸、方位角等基本參數(shù)，能夠定量描述相關(guān)經(jīng)濟屬性空間分布特征。其中，標(biāo)準(zhǔn)差橢圓的分布范圍即相關(guān)經(jīng)濟屬性空間分布的主要范圍。標(biāo)準(zhǔn)差橢圓的平均中心即相關(guān)經(jīng)濟屬性空間分布的重心，長軸方向代表相關(guān)經(jīng)濟屬性分布較多的方向，短軸方向代表相關(guān)經(jīng)濟屬性分布較少的方向，長短軸之間的差值越大（扁率越大），說明相關(guān)經(jīng)濟屬性的方向性就越強。

（三）灰色預(yù)測模型

灰度預(yù)測模型能夠?qū)r間序列短、統(tǒng)計數(shù)據(jù)少、信息不完全的系統(tǒng)進行預(yù)測，這是灰色預(yù)測模型區(qū)別于一般動態(tài)預(yù)測方法的獨特優(yōu)勢。本文采用灰色預(yù)測模型的GM（1，1）模型對于標(biāo)準(zhǔn)差橢圓平均中心坐標(biāo)的經(jīng)度和緯度、長半軸、短半軸和旋轉(zhuǎn)角進行預(yù)測，從而揭示未來山東省縣域經(jīng)濟發(fā)展重心變遷和空間分布狀況。

（四）數(shù)據(jù)來源

人均GDP是經(jīng)濟研究廣泛使用的、體現(xiàn)經(jīng)濟發(fā)展水平的主要指標(biāo)。本文2013-2022年人均GDP數(shù)據(jù)，由中經(jīng)網(wǎng)縣域經(jīng)濟年度庫的縣市區(qū)地區(qū)生產(chǎn)總值和年末戶籍人口數(shù)量計算得到，個別缺失的縣域歷史統(tǒng)計數(shù)據(jù)由有關(guān)地市統(tǒng)計年鑒及縣市區(qū)國民經(jīng)濟和社會發(fā)展情況公報補充。本文使用的空間計量分析工具為ArcGis10.8軟件，建立山東省縣域空間統(tǒng)計數(shù)據(jù)庫，然后運用ArcGis10.8相關(guān)模塊進行實證分析。

二、山東省縣域經(jīng)濟空間格局演變分析

（一）縣域經(jīng)濟發(fā)展的空間集聚演變分析

1. 空間分布模式分析。研究期內(nèi)山東省縣域經(jīng)濟發(fā)展空間聚集分布顯著。2013-2022年山東省縣域經(jīng)濟發(fā)展水平全局Moran’I值均為正值，且通過1%置信水平的顯著性水平檢驗，這表明山東省縣域經(jīng)濟發(fā)展顯著空間正相關(guān)，整體上呈現(xiàn)空間聚集分布特征，即經(jīng)濟發(fā)展水平相近的縣域在空間上聚集分布。如表1所示，全局Moran’I值總體呈下降走勢，且具有明顯的階段性特征：除個別年份外，2013-2019年全局Moran’I值下行加速，2020-2022年則下行趨緩。這表明，研究期內(nèi)山東省縣域經(jīng)濟空間聚集程度總體上不斷減弱，縣域經(jīng)濟空間差異總體上不斷縮小，且縮小趨勢由加速轉(zhuǎn)向平緩。白全民等（2013）、蒙莎莎等（2017）則認(rèn)為在所研究時期內(nèi)山東省縣域經(jīng)濟空間聚集效應(yīng)呈增強態(tài)勢。這說明，“十二五”以來，隨著山東省實施重點區(qū)域帶動戰(zhàn)略，山東黃河高效生態(tài)經(jīng)濟區(qū)和山東半島藍色經(jīng)濟區(qū)（“藍黃兩區(qū)”）上升為國家戰(zhàn)略，不斷優(yōu)化“兩區(qū)一圈一帶”“一群兩心三圈”區(qū)域發(fā)展布局，支持特色經(jīng)濟區(qū)加快發(fā)展，加快區(qū)域一體化進程，縣域經(jīng)濟發(fā)展更加協(xié)調(diào)均衡，從而呈現(xiàn)出新的空間分布趨勢和特征。

2. 空間聚類分析。研究期內(nèi)山東省縣域經(jīng)濟空間演進特征明顯，熱點區(qū)分布擴散，冷點區(qū)趨向集中，局部空間關(guān)聯(lián)性增強，極化效應(yīng)增強。參照白全民等（2013）研究方法，使用Quantile分位法對2013-2022年山東省縣域經(jīng)濟發(fā)展Getis’s-Ord Gi*統(tǒng)計值從低到高進行排列，依次劃分為冷點區(qū)、次冷點區(qū)、次熱點區(qū)和熱點區(qū)等四類聚類分布類型。研究結(jié)果顯示，2013-2022年，山東省縣域經(jīng)濟熱點區(qū)和冷點區(qū)發(fā)生較為明顯的變化。從空間聚集演變看，熱點區(qū)初始區(qū)域主要位于黃河流域、膠東半島地區(qū)，呈現(xiàn)左右兩翼分布態(tài)勢。隨著濰坊沿渤海縣域單元逐步成為熱點區(qū)，逐漸產(chǎn)生連接兩翼熱點區(qū)的“走廊”地帶，最終形成山東省藍黃經(jīng)濟帶“兩翼一廊”連片熱點區(qū)，該連片熱點區(qū)與山東省藍黃兩區(qū)或者省會經(jīng)濟圈和膠東經(jīng)濟圈地理范圍基本一致。冷點區(qū)初始區(qū)域主要位于魯西南和魯南地區(qū)，此后魯西北冷點區(qū)擴大，魯南冷點區(qū)逐步縮小，最終形成魯西南和魯西北兩大連片冷點區(qū)。從冷熱點區(qū)分布演變趨勢看，熱點區(qū)逐步向山東半島方向擴大，而冷點區(qū)逐步向西魯西南、魯西北方向集中。白全民等（2013）認(rèn)為熱點區(qū)位于黃河南岸和半島“青煙威”地區(qū)，并呈“一”字形分布，本文研究表明新時期山東省熱點區(qū)進一步擴散，藍黃兩區(qū)聯(lián)動效應(yīng)明顯增強，縣域經(jīng)濟發(fā)展聚集呈魯西北和魯西南低洼帶、藍黃兩區(qū)經(jīng)濟隆起帶的二元聚集特征。

（二）縣域經(jīng)濟發(fā)展的時空格局變遷分析

1. 空間分布重心。標(biāo)準(zhǔn)差橢圓平均中心即經(jīng)濟重心。如表2所示，從遷移軌跡看，研究期內(nèi)山東省縣域經(jīng)濟重心總體上向東北方向遷移。2013年山東省縣域經(jīng)濟重心在（118.262°E，36.533°N）處，2022年則遷移到（118.387°E，36.566°N）處，縣域經(jīng)濟重心以平均1.31千米/年的速度向東北方向遷移，總移動距離為11.75km，遷移速度、幅度均較小。研究期內(nèi)，重心均位于青州市境內(nèi)。2016年經(jīng)濟重心主要向東遷移，總遷移距離為1.80km，其中向東遷移1.79km，向南遷移0.22km。2019年經(jīng)濟重心向主要向東遷移，總遷移距離為8.91km，其中向東遷移8.85km，向北遷移1.11km。2022年經(jīng)濟重心主要向北遷移，總遷移距離為2.83km，其中向東遷移0.84km，向北遷移2.78km。從遷移過程看，經(jīng)濟重心向東遷移距離大于向北遷移距離，但向東遷移趨緩、向北遷移加速。這表明山東東部經(jīng)濟經(jīng)歷了快速發(fā)展到加速發(fā)展再到平穩(wěn)發(fā)展的過程，山東北部經(jīng)濟處于較小幅度加速發(fā)展過程中，東西差距擴大趨勢得到遏制，南北差距相對小幅擴大，區(qū)域經(jīng)濟整體協(xié)調(diào)性增強，這與山東省“兩區(qū)一圈一帶”“一群兩心三圈”等區(qū)域發(fā)展戰(zhàn)略密切相關(guān)。

2. 空間分布范圍。如表2所示，研究期內(nèi)山東省縣域經(jīng)濟發(fā)展標(biāo)準(zhǔn)差橢圓分布范圍總體上呈擴張態(tài)勢，由8.36×104km2增長到8.55×104km2，這表明標(biāo)準(zhǔn)差橢圓之外的縣域經(jīng)濟發(fā)展增速更快，縣域經(jīng)濟發(fā)展空間分布范圍更加分散。山東省推進半島藍色經(jīng)濟區(qū)、黃河高效生態(tài)經(jīng)濟區(qū)發(fā)展及推進省會經(jīng)濟圈、膠東半島經(jīng)濟圈、魯南經(jīng)濟圈區(qū)域一體化進程，城市群和都市圈高效集約發(fā)展，對區(qū)內(nèi)及周邊縣域產(chǎn)生帶動輻射效應(yīng)。空間分布范圍基本上東至青島市區(qū)、萊陽市、棲霞市、龍口市，西至陽谷縣、鄆城縣、梁山縣、嘉祥縣，南起平邑縣、費縣、沂南縣、莒縣，北抵惠民縣、濱州市區(qū)、利津縣、東營市區(qū)，涵蓋膠東經(jīng)濟圈、省會經(jīng)濟圈大部分地區(qū)以及魯南經(jīng)濟圈一部分地區(qū)。

3. 空間分布方向。研究期內(nèi)山東省縣域經(jīng)濟發(fā)展標(biāo)準(zhǔn)差橢圓旋轉(zhuǎn)角變化范圍在72.143°～70.496°，縣域經(jīng)濟空間分布方向為東北—西南方向。標(biāo)準(zhǔn)差橢圓旋轉(zhuǎn)角由72.143°逆時針旋轉(zhuǎn)為70.496°，表明研究期內(nèi)西北—東南方向（短軸方向）縣域經(jīng)濟發(fā)展在推動山東省縣域經(jīng)濟時空演變上起到了主要拉力作用。研究期內(nèi)橢圓長半軸由221.97km增長到227.11km，短半軸由119.95km增長到120.14km，山東省縣域經(jīng)濟表現(xiàn)在東北-西南方向（長軸方向）、西北－東南方向（短軸方向）均擴張，但東北－西南方向（長軸方向）擴張更為明顯，山東省沿空間分布方向的縣域經(jīng)濟溢出效應(yīng)更為顯著。

（三）縣域經(jīng)濟發(fā)展的空間分布格局預(yù)測

對山東省縣域經(jīng)濟發(fā)展標(biāo)準(zhǔn)差橢圓各參數(shù)構(gòu)建灰色預(yù)測模型GM（1，1）時間序列模型，級比檢驗值均處于標(biāo)準(zhǔn)范圍區(qū)間，各參數(shù)適合進行GM（1，1）模型構(gòu)建。采用后驗差檢驗法對預(yù)測結(jié)果進行檢驗，后驗差比均小于0.35，小誤差概率均大于0.9，說明GM（1，1）模型預(yù)測結(jié)果可信度較高。從GM（1，1）模型預(yù)測結(jié)果中，選取2030、2035年山東省縣域經(jīng)濟發(fā)展標(biāo)準(zhǔn)差橢圓相關(guān)預(yù)測參數(shù)，見表3。

預(yù)測結(jié)果顯示，未來山東省縣域經(jīng)濟重心遷移趨勢為西北方向。其中，2035年山東省縣域經(jīng)濟發(fā)展重心由（118.387°E，36.566°N）遷移到（118.292°E，36.874°N），呈西北向移動，分別向西移動8.48km，向北移動34.25km，表明山東北部環(huán)渤海地區(qū)和西部地區(qū)將影響山東省縣域經(jīng)濟發(fā)展空間分布的整體格局。預(yù)測結(jié)果還顯示，未來山東省縣域經(jīng)濟標(biāo)準(zhǔn)差橢圓長半軸標(biāo)準(zhǔn)差、短半軸標(biāo)準(zhǔn)差、旋轉(zhuǎn)角分別呈減小、增加及減小態(tài)勢。其中，2035年旋轉(zhuǎn)角度由2022年的70.50°縮小為64.83°，這表明山東省縣域經(jīng)濟發(fā)展空間分布東北—西南格局進一步強化。標(biāo)準(zhǔn)差橢圓面積由2022年的8.55×104km2縮小到7.94×104km2，縮小7.13%，長半軸標(biāo)準(zhǔn)差由2022年的221.97km下降到2035年的202.88km，下降8.6%，短半軸標(biāo)準(zhǔn)差由2022年的120.14km增加到2035年的124.51km，增加3.6%，且扁平率從2022年的0.53升至2035年的0.61，表明未來山東省縣域經(jīng)濟發(fā)展沿長軸方向的空間溢出效應(yīng)弱化，而沿短軸方向空間溢出效應(yīng)增強，縣域經(jīng)濟差異性會整體降低。

三、結(jié)論與討論

第一，研究期內(nèi)山東省縣域經(jīng)濟發(fā)展水平整體上呈顯著空間正相關(guān)，聚集分布特征顯著，縣域經(jīng)濟空間聚集效應(yīng)總體呈平緩減弱趨勢。山東省縣域經(jīng)濟局部空間關(guān)聯(lián)性強，極化效應(yīng)比較顯著，形成熱點區(qū)、冷點區(qū)連片集中聚集分布態(tài)勢。

第二，研究期內(nèi)山東省縣域經(jīng)濟時空演變特征顯著。縣域經(jīng)濟重心總體上向東北方向遷移，向東遷移趨緩、向北遷移加速，但遷移幅度有限，重心位置變化不大。空間分布范圍略呈擴張態(tài)勢，擴張幅度較小。空間分布方向呈東北－西南走向，沿空間分布方向的縣域經(jīng)濟溢出效應(yīng)較短軸方向的溢出效應(yīng)略強。這些演變特征表明，新時期山東省縣域經(jīng)濟空間差異性擴大不顯著，且變化趨向穩(wěn)定。

第三，從未來變化看，山東省縣域經(jīng)濟重心向西北方向移動，空間分布主軸仍然保持在東北－西南方向，但是沿長軸方向的空間溢出效應(yīng)減弱，沿短軸方向的空間溢出效應(yīng)增強，縣域經(jīng)濟差異性將進一步降低。

參考文獻：

[1]劉玉，潘瑜春，陳秧分.山東省縣域經(jīng)濟發(fā)展的時空動態(tài)研究[J].地理經(jīng)濟，2012，32（05）：43-48.

[2]蒙莎莎，張曉青，張玉澤，黃彩虹，尹向來，白晨.山東省縣域經(jīng)濟空間格局演變及驅(qū)動機制研究[J].華東經(jīng)濟管理，2017，31（12）：27-34.

[3]杜霞，錢宏勝，吳殿廷.山東省縣域經(jīng)濟的空間分異及其成因[J].城市問題，2015（08）：97-103.

[4]王慧，孫晉芳，余靜.山東省縣域經(jīng)濟時空分異演化特征研究[J].地域研究與開發(fā)，2014，33（05）：16-20+33.

（作者單位：中共臨朐縣委黨校）