山東省縣域經(jīng)濟空間格局演變及驅動機制研究

蒙莎莎，張曉青，張玉澤，黃彩虹，尹向來,白 晨

（1.山東師范大學 地理與環(huán)境學院，山東 濟南 250358；2.濟南大學 商學院，山東 濟南 250002）

山東省縣域經(jīng)濟空間格局演變及驅動機制研究

蒙莎莎1，張曉青1，張玉澤2，黃彩虹1，尹向來1,白 晨1

（1.山東師范大學 地理與環(huán)境學院，山東 濟南 250358；2.濟南大學 商學院，山東 濟南 250002）

文章以山東省103個縣域單元為研究對象，選取多項指標，利用標準差橢圓、空間變差函數(shù)、探索性空間數(shù)據(jù)分析及動態(tài)面板數(shù)據(jù)模型分析2005-2015年山東省縣域單元經(jīng)濟空間格局演化特征和地域分異驅動機制。結果表明：山東省縣域經(jīng)濟總體格局呈東北—西南走向，重心向東遷移且速度加快，主方向呈先擴張后收縮態(tài)勢；縣域經(jīng)濟空間差異、空間自組織性及結構性分異有所增強，東北—西南向分異顯著，內(nèi)在結構形態(tài)異質性未隨經(jīng)濟發(fā)展而減弱；全局空間關聯(lián)性表明經(jīng)濟空間集聚程度增強，呈顯著核心—外圍二元結構，局部空間關聯(lián)性表明縣域經(jīng)濟增長極與低洼區(qū)兩極分化嚴重，經(jīng)濟發(fā)達區(qū)密布于東部沿海，魯西南與魯西北等經(jīng)濟低洼區(qū)連片分布，空間擴散和梯度轉移是其主要原因；分權化、工業(yè)化水平、市場化和對外開放程度是影響山東省縣域經(jīng)濟格局變遷的重要動力機制。

空間格局；標準差橢圓；空間變差函數(shù)；驅動機制；動態(tài)面板數(shù)據(jù)模型；山東省

一、引 言

區(qū)域經(jīng)濟空間格局是經(jīng)濟活動在空間的表現(xiàn)形式，主要反映經(jīng)濟要素空間分布形態(tài)、經(jīng)濟活動的區(qū)位特點以及在地域空間中的相互作用機理［1］，對實現(xiàn)資源空間高效配置及區(qū)域間協(xié)調發(fā)展具有重要意義。隨著區(qū)域經(jīng)濟全球化發(fā)展，中國經(jīng)濟空間格局在資源要素空間、產(chǎn)業(yè)技術空間和地緣戰(zhàn)略空間等多因素交互作用下呈現(xiàn)地域性和復雜化的發(fā)展趨勢，而縣域單元作為衡量經(jīng)濟發(fā)展水平及各種生產(chǎn)要素有效聚合的區(qū)域經(jīng)濟綜合體，其經(jīng)濟協(xié)調發(fā)展將對區(qū)域資源優(yōu)化配置起決定性作用，其規(guī)模與水平將直接影響我國工業(yè)化與新型城鎮(zhèn)化發(fā)展進程，“國家十三五規(guī)劃綱要”指出：推動新型城鎮(zhèn)化和新農(nóng)村建設，提升縣域經(jīng)濟支撐輻射能力。我國正處于工業(yè)化中后期階段，大中城市要素功能、生產(chǎn)功能逐漸向縣域單元轉化，縣域經(jīng)濟空間依賴性與溢出效應逐漸顯現(xiàn)，局部經(jīng)濟空間集聚特征日益深化，縣域經(jīng)濟空間演化階段性致使空間格局不平衡性問題凸出，因此，對縣域經(jīng)濟空間格局演變過程及影響因素的深入剖析，不僅成為學術界關注的焦點，而且對發(fā)達地區(qū)保持經(jīng)濟競爭力、欠發(fā)達地區(qū)生產(chǎn)要素實現(xiàn)跨區(qū)域流動、區(qū)域經(jīng)濟協(xié)調發(fā)展與空間協(xié)作具有重要戰(zhàn)略和現(xiàn)實意義。

眾多學者已對區(qū)域經(jīng)濟格局變動趨勢展開了卓有成效的研究。配魯增長極理論、核心邊緣理論、極化-涓滴效應、“點—軸”漸進擴散理論等為揭示空間結構演變過程提供強有力理論支撐［2-5］。①研究區(qū)域尺度和類型區(qū)趨于多元化，研究尺度涵蓋國家、東中西三大地帶等宏觀尺度，經(jīng)濟區(qū)、省域等中觀經(jīng)濟單元及縣域等微觀尺度單元［6-11］，地域類型由沿海發(fā)達區(qū)轉向山區(qū)、省際邊緣等貧困經(jīng)濟帶［1，12-14］；②研究方法趨于多元化，主要包括變異系數(shù)、錫爾系數(shù)、基尼系數(shù)等數(shù)理統(tǒng)計方法，近年來，ESDA、空間計量經(jīng)濟學、空間變差函數(shù)、神經(jīng)網(wǎng)絡模型等方法更注重空間因素溢出關聯(lián)效應對空間結構演化的影響［15-18］，研究范式由計量學方法向格局、過程、機理多角度耦合研究發(fā)展［19-20］；③空間格局演變驅動機制由單因素獨立作用發(fā)展為多因素交互影響，主要探討資源稟賦、區(qū)位條件、交通運輸、產(chǎn)業(yè)結構、區(qū)域政策等要素跨區(qū)域流動綜合作用機理［21-22］，研究手段趨于創(chuàng)新，由定性分析向數(shù)學模型定量分析經(jīng)濟格局變遷軌跡轉變［9，23］；④指標體系構建基本以人均GDP單一指標衡量為主，經(jīng)濟發(fā)展多層面構建多指標評價研究較為缺乏［19-20，24］。綜上所述，國內(nèi)外學者對經(jīng)濟空間格局進行了系統(tǒng)化、多元化的理論與實證研究，但仍存在諸多不足之處，雖有部分學者對山東省縣域單元經(jīng)濟空間格局進行探討，但隨著近年來國內(nèi)外經(jīng)濟形勢的變化，山東省縣域經(jīng)濟空間格局呈現(xiàn)出新的特征，空間格局關聯(lián)性與異質性演化特征的研究方法有待進一步革新，且研究過程多采用截面數(shù)據(jù)，割裂了時間過程的連續(xù)性，不能真實地反應區(qū)域空間格局的演化規(guī)律，影響因素的揭示以定性分析為主，借鑒計量經(jīng)濟學動態(tài)面板模型的研究較少。

山東省作為縣域經(jīng)濟競爭力和發(fā)展?jié)摿^大的東部沿海強省，是國家十三五時期“支持東部地區(qū)率先發(fā)展”政策實施與響應的重要省份，研究其縣域經(jīng)濟空間格局演化規(guī)律不僅有利于深化山東省經(jīng)濟發(fā)展水平，且對其他省份區(qū)域經(jīng)濟發(fā)展具有一定的啟示作用。隨著山東省“兩區(qū)一圈一帶”空間發(fā)展戰(zhàn)略的提出，自然稟賦、區(qū)位條件、經(jīng)濟基礎存在較大差異的縣域單元空間格局及經(jīng)濟協(xié)調發(fā)展問題亟待研究。因此，本文以山東省103個縣為研究對象，運用2005-2015年面板數(shù)據(jù)，基于標準差橢圓、空間變差函數(shù)、空間自相關分析等多元空間分析方法及動態(tài)面板數(shù)據(jù)模型揭示山東省縣域經(jīng)濟時空格局演變規(guī)律及影響機制，以期全面衡量縣域經(jīng)濟發(fā)展差異，為山東省經(jīng)濟發(fā)展政策提供理論支持。

二、指標體系與研究方法

（一）指標選取與數(shù)據(jù)來源

1.指標選取

為揭示山東省縣域單元經(jīng)濟發(fā)展真實水平，指標體系構建需全面充分考慮山東省發(fā)展背景和區(qū)域經(jīng)濟結構演進特點?？紤]到各縣域單元資源稟賦，區(qū)位條件等方面的差異，根據(jù)各區(qū)域經(jīng)濟結構內(nèi)在聯(lián)系，從指標選取科學性、客觀性等方面出發(fā)，以人均GDP（gdp）代表經(jīng)濟發(fā)展程度，二、三產(chǎn)業(yè)結構（ti）表征產(chǎn)業(yè)結構水平，工業(yè)發(fā)展程度以人均工業(yè)總產(chǎn)值（ind）來衡量，分權化（dec）主要通過人均地方財政支出表征，對外開放程度用人均出口額（exp）表示，以人均固定資產(chǎn)投資（fix）表征區(qū)域發(fā)展戰(zhàn)略，以人均社會消費品零售額（mar）表示市場化發(fā)展程度。

2.數(shù)據(jù)來源

本文以2005-2015年面板數(shù)據(jù)為主要分析樣本，根據(jù)最新的山東省行政區(qū)劃方案，山東省共有137個縣級行政區(qū)，其中包括51個市轄區(qū)、28個縣級市和58個縣，由于相鄰市轄區(qū)經(jīng)濟發(fā)展和生產(chǎn)要素流動相似性較高，因此將市轄區(qū)歸并為一個單元，同時考慮到數(shù)據(jù)單元一致性，將2015年以前的部分縣級市（兗州市、沾化縣、陵縣、文登市）由獨立縣域單元劃入相應市轄區(qū)，最后得到103個縣域單元。數(shù)據(jù)主要來源于《中國城市統(tǒng)計年鑒》、《中國區(qū)域經(jīng)濟統(tǒng)計年鑒》、《中國縣（市）社會經(jīng)濟統(tǒng)計年鑒》及各地市統(tǒng)計年鑒，為克服價格指數(shù)影響，以2000年價格為依據(jù)，其他年份按當年山東省GDP增長指數(shù)計算縣域單元不變價格的人均GDP。

（二）研究方法

1.標準差橢圓

本文利用標準差橢圓揭示山東省縣域單元經(jīng)濟格局的總體特征，標準差橢圓是揭示地理要素空間分布特征的有效方法。SDE主要由轉角θ、長軸標準差及短軸標準差等3要素組成，以多視角反映地理要素的空間分布形態(tài)、中心性和方向性等特點。

2.空間變差函數(shù)模型

本文引入空間變差函數(shù)分析山東省縣域經(jīng)濟發(fā)展空間變異規(guī)律?？臻g變差函數(shù)又稱為半變異函數(shù)，它是描述區(qū)域化變量隨機性和結構性的基本工具。設區(qū)域化變量Z（xi）滿足二階平穩(wěn)假設的隨機化變量，Z（xi）和Z（xi+h）分別是區(qū)域化變量Z（x）在空間位置xi和xi+h的實際值［i=1，2，3，…，N（h）］，空間變差函數(shù)的公式為：

以h為橫坐標、γ（h）為縱坐標，得到半變異函數(shù)如圖1所示。

圖1 理論方差

其中，C0為塊金值，表示區(qū)域化變量非連續(xù)性變異；C+C0為基臺值，表示隨距離增加半變異函數(shù)達到相對穩(wěn)定時的常數(shù)值；a為變程，表征變異函數(shù)γ（h）達到穩(wěn)定時的距離；C0/C+C0為塊金系數(shù)，反映隨機部分引起空間異質性的程度。

分維數(shù)是表征變異函數(shù)特性的第四個變量，γ（h）和間隔距離h關系為：

分維數(shù)D為雙對數(shù)直線回歸方程斜率，表示變異函數(shù)曲線曲率，可作為隨機變異的量度。其值越接近于2，表明空間呈隨機均衡分布。

3.空間自相關分析

空間自相關分析是衡量兩個相鄰空間單元經(jīng)濟發(fā)展水平變化趨勢的空間統(tǒng)計分析方法，本文主要選取2005年，2010年和2015年三個時間節(jié)點分析山東省縣域單元空間關聯(lián)性格局的變動趨勢?？臻g自相關包括全局空間自相關和局部空間自相關。

（1）全局空間自相關。Global Moran’s I是衡量研究對象在空間上是否存在顯著關聯(lián)的方法［25］。其公式為：

（2）局部空間自相關。Local Moran’s I值是從局部探索經(jīng)濟發(fā)展的高值與低值，對于縣域單元i。局部自相關公式為：

Ii可以將縣域單元劃分為高高、低低、低高和高低四種空間關聯(lián)格局。由于山東省存在海島縣，各縣域單元面積相差懸殊，如果采用queen法則確定權重，容易產(chǎn)生較大誤差，在計算各縣平均鄰域數(shù)目后，6鄰法（K6）較為合理。

三、山東省縣域經(jīng)濟發(fā)展空間格局演變特征

（一）總體分布特征：重心向東遷移，經(jīng)濟總體呈東北—西南走向，強集聚性

倘若不考慮自然要素和經(jīng)濟社會要素的作用，社會經(jīng)濟、人口等的地理空間分布應該是均衡、隨機的［26］。各年份標準差橢圓中心均為該年份的經(jīng)濟重心，從其遷移軌跡看，2005-2015年山東省縣域經(jīng)濟重心向東遷移，10年間總體遷移了19.06km，總體變動幅度較小。標準差橢圓東至龍口、棲霞、海陽、即墨，西連東平、東阿、汶上和曲阜，南接青島市轄、五蓮、莒縣、沂南，北到商河、惠民、濱州市轄、墾利，涵蓋省會城市都市圈大部分地區(qū)且經(jīng)濟水平相對較高。重心均位于青州市境內(nèi)，2010年縣域單元經(jīng)濟重心以1.59km/Y的速度向西南方向遷移，2015年則以2.78km/Y的速度遷移至東部，遷移速度明顯加快，說明經(jīng)濟發(fā)展向東部傾斜。重心遷移方向及遷移速度表明山東省經(jīng)濟發(fā)展呈現(xiàn)明顯階段性特征，魯東地區(qū)經(jīng)濟平穩(wěn)增長到魯東魯西地區(qū)協(xié)調發(fā)展，再至魯東地區(qū)經(jīng)濟快速發(fā)展，這與山東省不同時期區(qū)域發(fā)展戰(zhàn)略調整密切相關。

轉角變化范圍在63.5°～65.5°之間，整體上，山東省縣域經(jīng)濟空間格局呈現(xiàn)東北-西南走向。2005-2015年，轉角 θ由 65.34°縮小至 63.75°，再增加至65.11°，東北—西南方向呈現(xiàn)先強化后弱化趨勢，至2015年，東西向經(jīng)濟發(fā)展趨勢進一步增強。主軸標準差2005-2015年由235.43km減少至233.94km，說明此階段縣域經(jīng)濟基本保持東北—西南走向，且呈先擴張后收縮的態(tài)勢。其中2005-2010年呈現(xiàn)緩慢擴張趨勢，年均增長0.56km，2015年快速收縮至233.94km，年均下降1.07km，說明縣域經(jīng)濟在主方向出現(xiàn)極化現(xiàn)象。輔軸標準差則由2005年的120.88km下降至2015年118.63km，下降速度為0.23km/Y，說明縣域經(jīng)濟在西北—東南向呈相對較弱極化分布。具體如圖2所示。

圖2 2005-2015年山東省經(jīng)濟空間格局標準差橢圓

（二）縣域經(jīng)濟空間變異性格局：空間結構自組織性與內(nèi)在結構異質性增強

空間變差函數(shù)能對區(qū)域化變量空間結構和空間變異規(guī)律的隨機性和結構性進行有效分析［27］，數(shù)據(jù)符合正態(tài)分布是采用空間變差函數(shù)模型分析各經(jīng)濟單元空間演化格局的前提，數(shù)據(jù)非正態(tài)分布會增加變異函數(shù)估計誤差［8］。經(jīng)Shapiro-Wilk法檢驗，3個時間節(jié)點數(shù)據(jù)均為非正態(tài)分布，經(jīng)過對數(shù)變換后均通過檢驗，將其作為空間變量賦值給各縣市區(qū)幾何中心，并對各時間節(jié)點數(shù)據(jù)采用擬合度最高的模型進行模擬，并計算各年份不同方向分維數(shù)，進行Kriging插值，由基臺值、塊金值和塊金系數(shù)指標變化看，山東省經(jīng)濟發(fā)展空間差異不斷擴大，塊金系數(shù)單調遞增，變程在3個年份的既定步長下變化范圍較大，在80～2 133km之間波動，表明在經(jīng)濟空間差異不斷增大過程中，空間關聯(lián)產(chǎn)生的結構性分異隨不同時期數(shù)據(jù)變異成分變化而有所增強。利用OLS擬合空間變差模型3個年份的模型分別為指數(shù)模型、球面模型和高斯模型，可知各縣域單元在不同經(jīng)濟發(fā)展時期空間結構異質性較強，且模型擬合程度越來越高（R2越接近1，擬合程度越高），說明隨著時間推移，空間經(jīng)濟結構差異趨于明顯，經(jīng)濟發(fā)展空間自組織性逐年增強。具體見表1所列。

表1 山東省空間變差函數(shù)擬合參數(shù)

由變差函數(shù)的分維數(shù)來看，全方位分維數(shù)先減后增，擬合程度逐年增強，但變動幅度較小，在各方向分維數(shù)中，東北—西南分維數(shù)最低，擬合效果相對較好，表明東北-西南方向縣域單元空間格局脫離均質分布理想值，空間差異最大。東南—西北分維數(shù)逐年下滑，擬合度不斷增高，可見魯西北地區(qū)與魯東南地區(qū)與空間均質分布理想值偏離度增大，經(jīng)濟空間分異性增強。2015年各方向分維數(shù)有一定程度提升，表明“兩區(qū)一圈一帶”等政策實施對協(xié)調縣域經(jīng)濟有顯著作用。具體見表2所列。

表2 山東省空間格局變差函數(shù)分維數(shù)

克里金插值3D擬合圖演化過程展現(xiàn)了山東省縣域單元經(jīng)濟發(fā)展空間結構變異程度和分布格局。2005年，縣域單元經(jīng)濟格局高低起伏，落差較大，山東半島青煙威地區(qū)隆起為經(jīng)濟至高點的極化中心，周邊地區(qū)凹陷顯著，核心邊緣形狀凸顯，經(jīng)濟空間集聚態(tài)勢趨于增強。2010年縣域單元經(jīng)濟高低變異性與2005年類似，經(jīng)濟發(fā)展極核中心向西部縣域單元轉移，魯東沿海地帶仍為經(jīng)濟高地密布區(qū)，而經(jīng)濟低谷區(qū)由魯西南菏澤向魯西北德州等地區(qū)擴散。2015年經(jīng)濟空間格局得到一定程度延續(xù)，曲面起伏增強，表明區(qū)域經(jīng)濟內(nèi)在結構形態(tài)的異質性與波動性增強。具體如圖3所示。

圖3 山東省經(jīng)濟空間格局演變

（三）縣域經(jīng)濟空間關聯(lián)性分析

1.全局空間關聯(lián)性分析

山東省縣域單元經(jīng)濟發(fā)展空間集聚趨勢明顯。基于Geoda軟件測算2005-2015年103個縣域單元經(jīng)濟發(fā)展程度全局Moran’s I值，所有值均為正，標準化檢驗統(tǒng)計量Z值均大于5%置信水平的閾值1.96，說明山東省縣域經(jīng)濟發(fā)展在三個時間節(jié)點表現(xiàn)為顯著正相關，空間集聚態(tài)勢明顯。由表3可知，2005-2015年Moran’s I值呈逐年增長趨勢，且空間關系無明顯波動，隨著各縣域單元城鎮(zhèn)化、工業(yè)化、信息化發(fā)展及生產(chǎn)力水平的大幅提升，區(qū)域經(jīng)濟發(fā)展較強的核心地區(qū)地位凸顯，區(qū)域經(jīng)濟空間結構呈現(xiàn)核心—外圍二元結構形態(tài)，空間集聚程度與空間依賴性逐步增強。

表3 2005-2015年山東省縣域經(jīng)濟發(fā)展全局Moran’s I值

2.局部空間關聯(lián)性分析

縣域單元局部空間關聯(lián)性增強，空間集聚和極化效應顯著。由表4可見，山東省縣域空間自相關類型數(shù)量逐年上升，低-低類型比例在13.5%～21.4%之間，高-高類型區(qū)占比5.8%～9.8%之間，其余兩類型之和為3.8%～5.8%之間，縣域單元空間分布格局空間集聚明顯。同質單元由2005年的19.42%增加至2015年的30.1%，異質性單元由2005年的4.85%增加至2015年的5.82%，同質性單元比重及增長速度遠高于異質性縣域單元，表明經(jīng)濟空間格局分布不均衡且具有明顯極化效應。低-低類型區(qū)主要集中在魯西平原地帶的菏澤、德州、聊城南部地區(qū)，且連片分布。此類地區(qū)以第一產(chǎn)業(yè)為主，經(jīng)濟基礎薄弱，交通通達性低，成為縣域單元經(jīng)濟滯后區(qū)。高-高類型區(qū)主要分布在煙臺北部龍口和蓬萊，威海市東部榮成和乳山，呈集簇狀，集聚特征顯著，由于地處沿海，出口貿(mào)易發(fā)達，人均GDP處于縣域單元經(jīng)濟金字塔頂端。高-低與低-高類型區(qū)則分布較為分散，主要有魯西地區(qū)德州、聊城及臨沂市轄區(qū)。

表4 2005-2015年山東省縣域經(jīng)濟局部空間集聚狀態(tài)

由圖4可見，山東省縣域單元經(jīng)濟格局呈現(xiàn)典型核心-邊緣二元結構，經(jīng)濟兩極分化嚴重，經(jīng)濟增長極與經(jīng)濟低洼區(qū)差距較大，馬太效應嚴重。

縣域單元空間格局演變是經(jīng)濟擴散與經(jīng)濟空間梯度轉移共同作用的結果?？臻g集聚類型演化反映縣域經(jīng)濟發(fā)展速度差異，根本動因是經(jīng)濟基礎和區(qū)域發(fā)展政策分異。2010年高-高值類型區(qū)在保持2005年基礎上增加了墾利縣，而2015年萊州、膠州、平度、威海市等轉移為經(jīng)濟增長極。相對于2005年，2010年低-低類型區(qū)由魯西南地區(qū)向魯西北和魯南地區(qū)轉移且數(shù)量增多，主要分布在德州，聊城，臨沂和濟寧等魯西地區(qū)，2015年低-低類型區(qū)再次向魯西南地區(qū)轉移，說明山東省在西部經(jīng)濟隆起帶和魯南經(jīng)濟帶政策影響下，交通基礎設施軸線帶動作用明顯，加之東部地區(qū)經(jīng)濟擴散效應延伸，經(jīng)濟增效顯著。高低和低高類型空間單元數(shù)量逐步減少，說明山東省縣域單元經(jīng)濟極化現(xiàn)象主導作用有增強的趨勢。

圖4 2005-2015年山東省經(jīng)濟格局LISA圖

綜上所述，空間關聯(lián)類型發(fā)展歷程和規(guī)模數(shù)量演進表明魯西、魯西北和魯西南地區(qū)是山東省經(jīng)濟發(fā)展的動力來源，山東半島藍色經(jīng)濟區(qū)則依靠區(qū)位優(yōu)勢，技術先進和對外貿(mào)易通達度高等優(yōu)勢獲得優(yōu)先發(fā)展，成為經(jīng)濟熱點區(qū)。魯南臨沂山地區(qū)及魯西北德州、濱州等部分縣域單元，由于地形、氣候等自然條件影響，加之周邊區(qū)域競爭壓力大等問題，成為經(jīng)濟發(fā)展的冷點地區(qū)。經(jīng)濟極化和擴散效應是影響區(qū)域經(jīng)濟空間結構演進的主要原因，隨著經(jīng)濟技術進步與發(fā)展，中心和外圍地區(qū)將不斷轉移。

四、縣域經(jīng)濟空間結構演變的驅動機制

（一）模型構建

以上研究表明，山東省縣域經(jīng)濟發(fā)展空間格局存在顯著的地理空間分異現(xiàn)象，為進一步探討產(chǎn)生這種空間分異格局作用機理，此部分以人均GDP為被解釋變量，其余變量為解釋變量，構建的動態(tài)面板數(shù)據(jù)模型闡述縣域單元發(fā)展產(chǎn)生明顯地域分異的原因，并利用差分GMM和系統(tǒng)GMM兩種估計方法進行刻畫，由于經(jīng)濟增長具有滯后性，故將被解釋變量的滯后項作為解釋變量引入模型，同時為消除指數(shù)趨勢和異方差的影響，對除二、三產(chǎn)業(yè)結構以外的變量做對數(shù)轉換。構建動態(tài)面板模型形式如下：

（二）結果分析

為確保計量模型估計的有效性，并避免偽回歸的存在，本部分主要利用LLC、ADF-fisher和HT檢驗3種方法對數(shù)據(jù)平穩(wěn)性進行檢驗，檢驗在1%顯著水平上均拒絕存在單位根的原假設，且序列均為零階單整序列I（0），數(shù)據(jù)平穩(wěn)性較好?；谡`差修正模型的協(xié)整分析，Gt、Pt、Pa在5%顯著水平上，Ga在10%顯著水平上拒絕不存在協(xié)整關系的原假設，因此可以進行動態(tài)面板模型的構建。

根據(jù)表5模型整體分析結果，差分GMM和系統(tǒng)GMM都存在一階序列自相關及二階序列不相關，且Sargen檢驗表明模型工具變量選取較為合理，不存在過度識別問題，且從模型顯著水平來看，系統(tǒng)GMM由于考慮水平和差分方程，對多數(shù)變量的擬合效果明顯優(yōu)于差分GMM。lngdp的滯后項對經(jīng)濟水平增長有較強的正效應，反映了縣域層次經(jīng)濟發(fā)展具有較強的自我強化機制，經(jīng)濟發(fā)展程度高的縣域單元經(jīng)濟發(fā)展動力越強。分權化和工業(yè)化水平對山東省經(jīng)濟空間格局影響效力最大，市場化因素與對外開放程度次之，產(chǎn)業(yè)結構對空間格局異質性影響最弱，而區(qū)域發(fā)展戰(zhàn)略則間接影響縣域單元空間演變，作用效果還不十分顯著。

表5 差分GMM與系統(tǒng)差分模型分析結果

產(chǎn)業(yè)結構在差分GMM和系統(tǒng)GMM模型中影響系數(shù)為正，但差分GMM不顯著，系統(tǒng)GMM1%水平顯著，說明產(chǎn)業(yè)結構調整和變遷是影響縣域單元資源配置的重要因素，魯西菏澤、德州、聊城等縣域一產(chǎn)職能強度較高在0.8以上，第三產(chǎn)業(yè)對經(jīng)濟增長拉動效應不顯著，經(jīng)濟發(fā)展后勁不足。而魯東丘陵地區(qū)在“一核、兩極、三帶、三組團”政策聯(lián)動效應下，縣域產(chǎn)業(yè)結構協(xié)調度、產(chǎn)城融合和資源高效配置程度相對領先，產(chǎn)業(yè)職能分工優(yōu)勢較大。工業(yè)發(fā)展水平對經(jīng)濟增長的影響較大，工業(yè)化每增長1%將帶來人均GDP 0.128%的提升，在地域分布上，魯中地區(qū)淄博、東營等地區(qū)以工業(yè)為主的縣域產(chǎn)業(yè)發(fā)展水平較高，如東營廣饒縣工礦業(yè)和橡膠制造業(yè)對經(jīng)濟增長的貢獻達69%以上，工業(yè)為縣域經(jīng)濟發(fā)展提供了強有力支撐。

分權化因素對縣域經(jīng)濟增長的彈性系數(shù)為正且推動作用最大，人均財政支出每增加1%將推動縣域GDP增長0.15%，表明經(jīng)濟新常態(tài)政策下，經(jīng)濟發(fā)展對政府財政投入產(chǎn)生較大路徑依賴，地方財政支出是縣域單元經(jīng)濟增長主要推動力，由于山東省財政支出傾向于市轄區(qū)，而對縣域單元資金投入?yún)T乏，致使市轄區(qū)與縣域單元尺度上分權作用顯著。市場化因素與經(jīng)濟發(fā)展的彈性系數(shù)在差分和系統(tǒng)GMM模型中均為0.09，且在10%水平下顯著，市場化是經(jīng)濟平穩(wěn)運行的重要保障，自1990年以來，民營經(jīng)濟成為影響山東省縣域經(jīng)濟發(fā)展格局重要形式，市場化在作用于全要素生產(chǎn)和資源配置過程中影響經(jīng)濟增長，因而促使縣域經(jīng)濟單元投入要素流向資本配置效率更高的民營企業(yè)，將有利于縣域經(jīng)濟發(fā)展水平的提升。

對外開放程度與經(jīng)濟增長在差分和系統(tǒng)GMM中呈正相關，且通過5%和10%顯著性檢驗。三駕馬車之一的“出口”對縣域單元經(jīng)濟發(fā)展水平提升具有重要作用，魯東丘陵地區(qū)依靠優(yōu)越區(qū)位條件和“山東半島藍色經(jīng)濟區(qū)”等政策扶持發(fā)展外向型經(jīng)濟所必需的勞動力、技術等生產(chǎn)要素，并通過影響生產(chǎn)要素流動方向引起縣域單元區(qū)域增長分異，如以青島為中心的“東翼”已成為山東省經(jīng)濟集聚增長的熱點區(qū)。而魯西北德州、聊城和魯西南菏澤則處于外部生產(chǎn)要素隔離狀態(tài)，雖先后有“西部經(jīng)濟隆起帶”和“西突破菏澤”等政策扶持，但其經(jīng)濟發(fā)展水平和規(guī)模與東部沿海地區(qū)相差甚遠。

區(qū)域經(jīng)濟發(fā)展戰(zhàn)略與政策是區(qū)域經(jīng)濟格局演化的重要驅動力，人均固定資產(chǎn)投資在兩模型中對經(jīng)濟發(fā)展水平彈性系數(shù)影響較小且不顯著，其原因很可能因為人均固定資產(chǎn)投資并不是直接作用于經(jīng)濟增長，而是通過推進產(chǎn)業(yè)結構調整促進經(jīng)濟增長。政府對基礎設施等固定資產(chǎn)投資將更有利于縣域單元整體經(jīng)濟發(fā)展水平的提升。山東省“西翼”聊城、菏澤等內(nèi)陸地區(qū)在經(jīng)濟格局中處于滯后地位，加之獲得政府投資相對較少，經(jīng)濟發(fā)展緩慢；區(qū)位條件優(yōu)越的“東翼”青煙威縣域單元依托較高物質資本戰(zhàn)略性輸入作用于經(jīng)濟發(fā)展要素而成為經(jīng)濟發(fā)展高地。因此，山東省政府投資戰(zhàn)略傾向性及資本在區(qū)域間不均衡流動和配置是導致區(qū)域經(jīng)濟增長及空間分異的潛在因素。

五、結論及建議

本文采用標準差橢圓、空間變差函數(shù)和探索性空間分析技術分析2005-2015年山東縣域單元經(jīng)濟空間格局演化過程，并利用動態(tài)面板模型揭示其驅動機制。研究主要得出以下結論：

（1）2005-2015年縣域經(jīng)濟重心向東遷移，遷移速度逐漸加快；轉角變化范圍小，主方向始終保持東北—西南方向，縣域經(jīng)濟呈現(xiàn)先擴張后收縮態(tài)勢；

（2）空間變差函數(shù)表明山東省經(jīng)濟空間差異、空間自組織性及空間關聯(lián)產(chǎn)生的結構性分異有所增強，東北—西南向經(jīng)濟分異顯著，空間形態(tài)中經(jīng)濟至高點極核中心密布于青煙威地區(qū)且與周邊經(jīng)濟凹陷帶落差較大，經(jīng)濟低谷區(qū)分布于魯西南菏澤及魯西北德州、聊城等縣域，經(jīng)濟曲面起伏增大表明區(qū)域經(jīng)濟內(nèi)在結構形態(tài)異質性增強。

（3）從全局空間關聯(lián)性看，2005-2015年縣域單元經(jīng)濟空間集聚程度逐年增強，區(qū)域經(jīng)濟發(fā)展較強的核心區(qū)地位凸顯，區(qū)域經(jīng)濟空間結構呈現(xiàn)核心—外圍二元結構形態(tài)。從局部空間關聯(lián)性看，縣域經(jīng)濟兩極分化嚴重，高-高類型縣域單元主要分布在煙臺北部龍口和蓬萊、威海東部榮成及乳山，集聚特征顯著，低-低類型區(qū)主要集中在魯西平原地帶的菏澤、德州、聊城南部地區(qū)，且連片分布，亦呈典型核心邊緣結構?？臻g關聯(lián)類型南北空間轉移格局是經(jīng)濟擴散與經(jīng)濟空間梯度轉移共同作用的結果。

（4）從影響因素看，分權化因素和工業(yè)水平是主導山東省縣域單元空間格局演變的兩類最重要因素，市場化因素影響力趨于增強，縣域民營企業(yè)投資數(shù)量和方向不均衡性加劇經(jīng)濟空間異質性分布；區(qū)域發(fā)展戰(zhàn)略是資源稟賦空間、技術空間流轉的主要驅動力，但其通過作用于產(chǎn)業(yè)結構間接影響經(jīng)濟增長速度。對外開放程度是影響經(jīng)濟空間地域分布差異的重要原因，產(chǎn)業(yè)結構地域分異顯著，二產(chǎn)比重過高，第三產(chǎn)業(yè)發(fā)展仍需提質增效。

因此，未來經(jīng)濟發(fā)展中山東省應從縣域經(jīng)濟空間結構優(yōu)化、升級、重組等方面采取不同措施促進區(qū)域協(xié)調發(fā)展。魯西地區(qū)應加快對外開放力度，構筑內(nèi)陸地區(qū)經(jīng)濟高地，促進生產(chǎn)要素跨區(qū)域流動和資源優(yōu)化配置，延伸產(chǎn)業(yè)鏈條，強化縣域單元整合與集群式產(chǎn)業(yè)發(fā)展戰(zhàn)略；魯東地區(qū)應立足供給側改革，避免產(chǎn)能過剩和資源邊際效益下降，打破東西經(jīng)濟交流壁壘，形成區(qū)域協(xié)作發(fā)展新局面。同時，山東省應強化區(qū)域經(jīng)濟增長核心，依托“三縱三橫”及膠濟鐵路構筑綜合經(jīng)濟走廊，加強縣域經(jīng)濟發(fā)展聯(lián)動效應，優(yōu)化產(chǎn)業(yè)結構，改變政府投資拉動經(jīng)濟增長現(xiàn)狀，探索多元化、市場化經(jīng)濟發(fā)展模式。

［1］米楠，卜曉燕，米文寶.寧夏六盤山區(qū)縣域經(jīng)濟空間結構演化［J］.經(jīng)濟地理，2015（4）：40-46.

［2］Perroux F.Economic space：theory and application［J］.The Quarterly Journal of Economics，1950，64（1）：89-104.

［3］Friedman J R.Regional Development Policy-A Case Study of Venezuela［M］.Cambridge：The MIT Press，1966.

［4］馬浩，王慶金.山東經(jīng)濟非均衡協(xié)調發(fā)展研究［M］.北京：中國社會科學出版社，2015：10.

［5］陸大道.關于“點-軸”空間結構系統(tǒng)的形成機理分析［J］.地理科學，2002（1）：1-6.

［6］李小建，喬家君.20世紀90年代中國縣際經(jīng)濟差異的空間分析［J］.地理學報，2001（2）：136-145.

［7］徐建華，魯鳳，蘇方林，等.中國區(qū)域經(jīng)濟差異的時空尺度分析［J］.地理研究，2005（1）：57-68.

［8］周揚，李寧，吳文祥，等.1982-2010年中國縣域經(jīng)濟發(fā)展時空格局演變［J］.地理科學進展，2014（1）：102-113.

［9］馮興華，鐘業(yè)喜，陳琳，等.長江經(jīng)濟帶縣域經(jīng)濟空間格局演變分析［J］.經(jīng)濟地理，2016（6）：18-25.

［10］潘少奇，李亞婷，高建華.中原經(jīng)濟區(qū)經(jīng)濟聯(lián)系網(wǎng)絡空間格局［J］.地理科學進展，2014（1）：92-101.

［11］王鑫磊，陳斐，王海萍.中部地區(qū)經(jīng)濟格局動態(tài)演變的空間分析［J］.華東經(jīng)濟管理，2012（7）：39-43.

［12］王美霞，李民，周國華，等.武陵山片區(qū)縣域經(jīng)濟空間格局演化與優(yōu)化對策［J］.經(jīng)濟地理，2015（11）：45-53.

［13］張學波，楊成鳳，宋金平，等.中國省際邊緣縣域經(jīng)濟差異空間格局演變［J］.經(jīng)濟地理，2015（7）：30-38.

［14］章青果，鄭林.江西省區(qū)域經(jīng)濟差異及其空間格局研究［J］.華東經(jīng)濟管理，2006（3）：9-13.

［15］靳誠，陸玉麒.基于空間變差函數(shù)的長江三角洲經(jīng)濟發(fā)展差異演變研究［J］.地理科學，2011（11）：1329-1334.

［16］陳利，朱喜鋼，楊陽，等.基于空間計量的云南省縣域經(jīng)濟空間格局演變［J］.經(jīng)濟地理，2017（1）：40-49.

［17］趙磊，方成，黃武龍.浙江省縣域經(jīng)濟發(fā)展差異時空演變分析［J］.華東經(jīng)濟管理，2014（3）：6-11.

［18］張金萍，秦耀辰，張麗君，等.黃河下游沿岸縣域經(jīng)濟發(fā)展的空間分異［J］.經(jīng)濟地理，2012（3）：16-21.

［19］蔣天穎，華明浩，張一青.縣域經(jīng)濟差異總體特征與空間格局演化研究——以浙江為實證［J］.經(jīng)濟地理，2014（1）：35-41.

［20］王芳，高曉路.內(nèi)蒙古縣域經(jīng)濟空間格局演化研究［J］.地理科學，2014（7）：818-824.

［21］崔長彬，姜石良，張正河.河北縣域經(jīng)濟影響因素的空間差異分析——基于貝葉斯地理加權回歸方法［J］.經(jīng)濟地理，2012（2）：39-45.

［22］方葉林，黃震方，涂瑋，等.基于地統(tǒng)計分析的安徽縣域經(jīng)濟空間差異研究［J］.經(jīng)濟地理，2013（2）：33-38.

［23］紀小美，付業(yè)勤，陶卓民，等.福建省縣域經(jīng)濟差異的時空動態(tài)與變遷機制［J］.經(jīng)濟地理，2016（2）：36-44.

［24］李紅，丁嵩，劉光柱.邊緣省區(qū)縣域經(jīng)濟差異的空間格局演化分析——以廣西為例［J］.經(jīng)濟地理，2012（7）：30-36.

［25］程鈺，任建蘭，陳延斌，等.中國環(huán)境規(guī)制效率空間格局動態(tài)演變及其驅動機制［J］.地理研究，2016，35（1）：123-136.

［26］趙璐，趙作權.基于特征橢圓的中國經(jīng)濟空間分異研究［J］.地理科學，2014，34（8）：979-986.

［27］靳誠，陸玉麒.基于縣域單元的江蘇省經(jīng)濟空間格局演化［J］.地理學報，2009（6）：713-724.

A Study on Spatial Pattern Evolution and Driving Mechanism of County-level Economy in Shandong Province

MENG Sha-sha1,ZHANG Xiao-qing1,ZHANG Yu-ze2,HUANG Cai-hong1,YIN Xiang-lai1,BAI Chen1
（1.School of Geography and Environment,Shandong Normal University,Jinan 250358,China;2.Business School,University of Jinan,Jinan 250002,China）

Based on the standard deviation ellipse(SDE),spatial variogram,the exploratory spatial data analysis(ESDA)and dynamic pan?el data model,this paper analyzes the evolution characteristics of economic spatial pattern and the driving mechanism of region differentia?tion of the 103 county-level units in Shandong province from 2005 to 2015.The results show that:The overall pattern of county-level economy in Shandong province is northeast-southwest direction.The gravity of county-level economy center shifts to the east,and the speed is accelerating.The main direction presents the trend of expansion first and contraction subsequent;The spatial differences,spatial self-organization and structural differentiation of county-level economy have been enhanced,and the northeast-southwest direction has an obvious divergence.The heterogeneity of internal structure is not weakened with the economic development;The global spatial correla?tion indicates that the degree of spatial economic agglomeration is generally enhanced,which is a significant core-peripheral binary struc?ture.The local spatial correlation shows that county-level economic growth pole and low-lying areas are polarized seriously.The economi?cally developed areas are densely distributed in the eastern coastal areas,and the economic low-lying areas of the southwest and the north?west in Shandong province are contiguous.The spatial diffusion and the gradient transfer are the main reasons;Decentralization,industrial?ization level,marketization and the degree of opening to the outside world are an important dynamic mechanism,which affects the change of county-level economy pattern in Shandong province.

spatial pattern;standard deviation ellipse;spatial variogram;driving mechanism;dynamic panel data model;Shandong prov?ince

F061.5

A

1007-5097（2017）12-0027-08

10.3969/j.issn.1007-5097.2017.12.004

2017-04-24

山東省自然科學基金項目(ZR2014DM001)

蒙莎莎（1992-），女，山東德州人，碩士研究生，研究方向：區(qū)域可持續(xù)發(fā)展；

張曉青（1971-），女，山東泰安人，教授，博士生導師，通訊作者，研究方向：區(qū)域與城市經(jīng)濟；

張玉澤（1989-），男，山東濟南人，講師，博士，研究方向：區(qū)域生態(tài)環(huán)境與可持續(xù)發(fā)展；

黃彩虹（1992-），女，山東萊蕪人，博士研究生，研究方向：區(qū)域經(jīng)濟與生態(tài)環(huán)境協(xié)調發(fā)展;

尹向來（1994-），男，山東煙臺人，碩士研究生，研究方向：區(qū)域可持續(xù)發(fā)展；

白 晨（1993-），女，山東青島人，碩士研究生，研究方向：區(qū)域可持續(xù)發(fā)展。

［責任編輯：余志虎］