王 璇，史佳璐，慈福義

(山東師范大學(xué)a.地理與環(huán)境學(xué)院,b.經(jīng)濟(jì)學(xué)院，中國(guó) 濟(jì)南 250358)

現(xiàn)階段我國(guó)經(jīng)濟(jì)由高速增長(zhǎng)轉(zhuǎn)向高質(zhì)量發(fā)展，旅游經(jīng)濟(jì)仍保持較快增長(zhǎng)[1]。2013年習(xí)近平總書(shū)記在俄羅斯“中國(guó)旅游年”開(kāi)幕式的致辭中強(qiáng)調(diào)：“旅游是綜合性產(chǎn)業(yè)，是拉動(dòng)經(jīng)濟(jì)發(fā)展的重要?jiǎng)恿Α?。根?jù)2019年文化和旅游發(fā)展統(tǒng)計(jì)公報(bào)，我國(guó)2019年全年實(shí)現(xiàn)旅游總收入66 300億元，同比增長(zhǎng)11.1%，旅游發(fā)展再上新臺(tái)階。與此同時(shí)，作為旅游業(yè)總體規(guī)模居于全國(guó)第4位的山東省，因其擁有獨(dú)特的地理位置和豐富的旅游資源，旅游經(jīng)濟(jì)也再創(chuàng)新高，2019年山東省旅游總收入達(dá)11 087.3億元，比上年增長(zhǎng)12.1%。山東省作為旅游資源較為豐富的省份，東有青煙威的濱海旅游；南有沂蒙紅色文化、“三孔”；北部有黃河口文化；西有大運(yùn)河、水滸文化；中有“天下泉城”、泰山等。然而，2018年印發(fā)的《山東省全域旅游發(fā)展總體規(guī)劃(2018—2025年)》卻指出，山東省旅游業(yè)的規(guī)模和發(fā)展水平還沒(méi)有達(dá)到應(yīng)有的優(yōu)化程度，在省內(nèi)仍存在明顯的東西差異。因此，探究區(qū)域旅游經(jīng)濟(jì)發(fā)展水平空間格局演變特征對(duì)于認(rèn)識(shí)區(qū)域間旅游經(jīng)濟(jì)發(fā)展差異，優(yōu)化旅游空間結(jié)構(gòu)，實(shí)現(xiàn)山東省旅游經(jīng)濟(jì)高質(zhì)量發(fā)展具有重要的現(xiàn)實(shí)意義。

隨著旅游產(chǎn)業(yè)地位的不斷提升，旅游經(jīng)濟(jì)時(shí)空差異的研究也與日俱增，已成為國(guó)內(nèi)外旅游研究的熱點(diǎn)問(wèn)題之一。國(guó)外有關(guān)旅游經(jīng)濟(jì)的研究起步較早，研究較為完善[2-5]，相關(guān)研究側(cè)重于經(jīng)濟(jì)與旅游間的內(nèi)在關(guān)聯(lián)等方面[6-8]，研究方法包括主成分分析[9]、SBM-DE模型[10]等。國(guó)內(nèi)相關(guān)研究各異：從指標(biāo)的選擇來(lái)看，有旅游總收入、國(guó)內(nèi)旅游收入等單因子指標(biāo)[11-14]，也有多因子復(fù)合指標(biāo)[15-17]；在研究尺度上，仍以省域[18-21]或市域[22,23]為主；在研究方法上，有ESDA[1,24,25]、空間變差函數(shù)[24,26]、空間計(jì)量經(jīng)濟(jì)學(xué)模型[27]、網(wǎng)絡(luò)分析方法[28-30]等。

現(xiàn)階段有關(guān)旅游經(jīng)濟(jì)的研究不斷細(xì)化，研究區(qū)以中國(guó)省域?yàn)橹?，有關(guān)山東省旅游經(jīng)濟(jì)研究，特別是早期的旅游經(jīng)濟(jì)研究極少，而較長(zhǎng)時(shí)間序列的分析更能凸顯研究對(duì)象的變化特征。因此，考慮到數(shù)據(jù)的可獲取性和便捷性，本文選取1995—2018年山東省旅游經(jīng)濟(jì)發(fā)展數(shù)據(jù)，從長(zhǎng)時(shí)間序列中探究旅游經(jīng)濟(jì)發(fā)展?fàn)顩r。已有研究還表明，從時(shí)空交互視角下進(jìn)行旅游經(jīng)濟(jì)發(fā)展的研究極為稀缺，如今普遍使用的ESDA方法沒(méi)有納入時(shí)間因素，時(shí)間維度上的分析存在嚴(yán)重不足，易造成時(shí)間分析和空間分析二者間的割裂。因此，本文探索性地引入Rey提出的探索性時(shí)空數(shù)據(jù)分析ESTDA(exploratory time-space data analysis)法和時(shí)空掃描統(tǒng)計(jì)分析法，探討局部旅游經(jīng)濟(jì)的時(shí)空結(jié)構(gòu)、時(shí)空變遷和時(shí)空關(guān)系，實(shí)現(xiàn)時(shí)間和空間的共同測(cè)度?；诖?，本文以山東省16地市為研究對(duì)象，運(yùn)用標(biāo)準(zhǔn)差橢圓、ESTDA和時(shí)空掃描統(tǒng)計(jì)分析等方法，從時(shí)空交互視角探究山東省旅游經(jīng)濟(jì)發(fā)展水平的時(shí)空格局演變特征，以期為今后各地實(shí)現(xiàn)旅游經(jīng)濟(jì)高質(zhì)量發(fā)展及因地制宜地制定區(qū)域旅游發(fā)展戰(zhàn)略提供參考。

1 數(shù)據(jù)來(lái)源與研究方法

1.1 指標(biāo)選取和數(shù)據(jù)來(lái)源

根據(jù)2019年國(guó)務(wù)院對(duì)撤銷萊蕪市劃歸濟(jì)南市管轄的批復(fù)，本文以山東省16個(gè)地級(jí)單元作為研究樣本。鑒于數(shù)據(jù)的可獲得性及連續(xù)性，采用旅游總收入來(lái)反映區(qū)域旅游經(jīng)濟(jì)發(fā)展水平，其中本文所研究的濟(jì)南市旅游總收入由濟(jì)南市旅游總收入和萊蕪市旅游總收入求和得到。研究期缺失旅游總收入數(shù)據(jù)的部分年份或地區(qū)根據(jù)國(guó)內(nèi)旅游收入與旅游外匯收入相加得到[25,31]，旅游外匯收入根據(jù)當(dāng)年美元兌換人民幣匯率進(jìn)行換算。本研究所需數(shù)據(jù)均來(lái)源于《山東旅游統(tǒng)計(jì)便覽(2001—2018)》、《山東統(tǒng)計(jì)年鑒(1995—2019)》、《中國(guó)區(qū)域經(jīng)濟(jì)統(tǒng)計(jì)年鑒》、《山東旅游年鑒》、《中國(guó)旅游年鑒》、各地級(jí)市年鑒、各地級(jí)市統(tǒng)計(jì)年鑒、各地級(jí)市國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)等，部分缺失數(shù)據(jù)根據(jù)相鄰年份補(bǔ)齊。

1.2 研究方法

1.2.1 標(biāo)準(zhǔn)差橢圓及重心 標(biāo)準(zhǔn)差橢圓是定量描述地理要素空間分布特征的一種空間統(tǒng)計(jì)方法，其長(zhǎng)軸與短軸的方向分別反映了地理要素在二維空間上分布的主、次趨勢(shì)方向[32]。該方法主要以重心、轉(zhuǎn)角θ和X軸標(biāo)準(zhǔn)差及Y軸標(biāo)準(zhǔn)差等參數(shù)定量描述旅游經(jīng)濟(jì)要素的空間分布特征。具體計(jì)算公式詳見(jiàn)參考文獻(xiàn)[33,34]，計(jì)算均依托ArcGIS 10.2實(shí)現(xiàn)。

1.2.2 探索性時(shí)空數(shù)據(jù)分析(ESTDA) 本文運(yùn)用探索性時(shí)空數(shù)據(jù)分析(ESTDA)[35]描述市域單元旅游經(jīng)濟(jì)發(fā)展水平的時(shí)空結(jié)構(gòu)、時(shí)空變遷和時(shí)空關(guān)系特征。相比于探索性數(shù)據(jù)分析(ESDA)，ESTDA可以有效避免以往測(cè)量中出現(xiàn)的時(shí)空割裂，時(shí)間維度分析不足的問(wèn)題，實(shí)現(xiàn)時(shí)間和空間的共同測(cè)度。ESTDA方法主要包括LISA時(shí)間路徑、LISA時(shí)空躍遷和時(shí)空網(wǎng)絡(luò)3部分[36]。

1)LISA時(shí)間路徑

通過(guò)記錄Moran散點(diǎn)圖中由旅游經(jīng)濟(jì)變量及其空間滯后量組成的LISA坐標(biāo)，運(yùn)用LISA時(shí)間路徑測(cè)量坐標(biāo)的穩(wěn)定性水平，解釋旅游經(jīng)濟(jì)在局部空間內(nèi)的時(shí)空協(xié)同變化和動(dòng)態(tài)性特征[37,38]，局域空間動(dòng)態(tài)軌跡的幾何特征主要包括相對(duì)長(zhǎng)度和彎曲度兩方面，其計(jì)算公式[39]為

(1)

(2)

式中：Di為相對(duì)長(zhǎng)度；Fi為彎曲度；N為地級(jí)市數(shù)量；T為年度間隔；Li,t為第t年地級(jí)市i在Moran散點(diǎn)圖中的位置坐標(biāo)；d(Li,t,Li,t+1)為地級(jí)市i從第t年到t+1年的移動(dòng)距離；d(Li,1,Li,T)為地級(jí)市i從首年1995年到末年2018年的移動(dòng)距離。LISA時(shí)間路徑長(zhǎng)度用以體現(xiàn)區(qū)域局部空間結(jié)構(gòu)的動(dòng)態(tài)性特征[37]，Di值越高表明地級(jí)市i的旅游經(jīng)濟(jì)在散點(diǎn)圖上移動(dòng)的時(shí)間路徑越長(zhǎng)，Di值越低表明局部空間結(jié)構(gòu)較為穩(wěn)定；LISA時(shí)間路徑彎曲度用以體現(xiàn)區(qū)域局部空間結(jié)構(gòu)的波動(dòng)性特征[26]，F(xiàn)i值越高表明地級(jí)市i在散點(diǎn)圖上的移動(dòng)路徑越彎曲，F(xiàn)i值越低表明局部空間結(jié)構(gòu)波動(dòng)較小。

2)LISA時(shí)空躍遷

Rey等將Moran’s Ⅰ散點(diǎn)圖中各空間單元在特定時(shí)間間隔內(nèi)的移動(dòng)距離等屬性嵌入傳統(tǒng)馬爾可夫鏈中，提出轉(zhuǎn)移概率矩陣和時(shí)空躍遷理論[35]。本文運(yùn)用時(shí)空躍遷方法測(cè)度各地級(jí)市旅游經(jīng)濟(jì)的空間關(guān)聯(lián)類型的動(dòng)態(tài)變遷過(guò)程，根據(jù)各區(qū)域自身與鄰域單元之間的旅游經(jīng)濟(jì)發(fā)展水平轉(zhuǎn)移狀態(tài)將時(shí)空躍遷劃分為4種類型[40]：Ⅰ型：自身和鄰域均穩(wěn)定(HHt→HHt+1，HLt→HLt+1，LLt→LLt+1和LHt→LHt+1)；Ⅱ型：自身躍遷，鄰域穩(wěn)定(HHt→LHt+1，LHt→HHt+1，HLt→LLt+1和LLt→HLt+1)；Ⅲ型：自身穩(wěn)定，鄰域躍遷(HHt→HLt+1，LHt→LLt+1，HLt→HHt+1和LLt→LHt+1)；Ⅳ型：自身和鄰域均躍遷(HHt→LLt+1，LLt→HHt+1，LHt→HLt+1和HLt→LHt+1)，若自身與鄰域的躍遷方向一致，則為Ⅳ(1)型(HHt→LLt+1和LLt→HHt+1)，若自身與鄰域的躍遷方向相反，則為Ⅳ(2)型(LHt→HLt+1和HLt→LHt+1)。各地級(jí)市旅游經(jīng)濟(jì)的Moran’s Ⅰ散點(diǎn)圖空間穩(wěn)定性、時(shí)空變遷、時(shí)空凝聚和相對(duì)移動(dòng)率依據(jù)可能存在的躍遷數(shù)和各躍遷類型計(jì)算，可能存在的躍遷總數(shù)n=(2018-1995)×16=368，具體計(jì)算公式詳見(jiàn)參考文獻(xiàn)[41,42]。

3)時(shí)空交互可視化

相對(duì)于傳統(tǒng)的ESDA，ESTDA可以通過(guò)圖論來(lái)顯示旅游經(jīng)濟(jì)發(fā)展的時(shí)空交互格局，進(jìn)而揭示出傳統(tǒng)時(shí)空格局耦合分析中不顯著的時(shí)空關(guān)聯(lián)特征[40]。地級(jí)市間的旅游經(jīng)濟(jì)時(shí)空交互的網(wǎng)絡(luò)流通過(guò)計(jì)算臨接地級(jí)市LISA時(shí)間移動(dòng)軌跡的協(xié)方差相關(guān)系數(shù)得出[43]，并根據(jù)各地級(jí)市之間的連線進(jìn)行可視化表達(dá)。關(guān)聯(lián)強(qiáng)度大于0表明區(qū)域間為協(xié)同發(fā)展態(tài)勢(shì)；關(guān)聯(lián)強(qiáng)度小于0表明區(qū)域間為競(jìng)爭(zhēng)態(tài)勢(shì)[42]。

1.2.3 時(shí)空重排掃描統(tǒng)計(jì)分析 時(shí)空掃描統(tǒng)計(jì)量聚類方法是探測(cè)時(shí)空范圍內(nèi)聚集性是否顯著增加的聚類方法，其為美國(guó)學(xué)者Kulldorff等[44]于2005年提出，主要應(yīng)用于流行病學(xué)及社會(huì)經(jīng)濟(jì)學(xué)等領(lǐng)域，近年來(lái)開(kāi)始在地理學(xué)領(lǐng)域中應(yīng)用。文獻(xiàn)中廣泛使用的空間自相關(guān)分析方法由于沒(méi)有納入時(shí)間因素，不能判斷地理空間上聚集區(qū)域隨時(shí)間變化的趨勢(shì)，而時(shí)空掃描分析同時(shí)考慮了時(shí)間和空間兩方面，可以揭示地理空間上要素聚集區(qū)域隨時(shí)間變化的規(guī)律。因此，本文探索性地將這一方法在旅游經(jīng)濟(jì)時(shí)空分布問(wèn)題中加以應(yīng)用，以探測(cè)旅游經(jīng)濟(jì)的聚集性及其精確的地理位置。

2 旅游經(jīng)濟(jì)的時(shí)空發(fā)展?fàn)顩r

2.1 旅游經(jīng)濟(jì)時(shí)間特征

自1995年以來(lái)，山東省旅游經(jīng)濟(jì)發(fā)展迅速(見(jiàn)圖1)，其旅游總收入從1995年的170.45億元增至2018年的10 461.3億元。除2003年“非典”影響以外，每年均比往年有大幅增長(zhǎng)，旅游經(jīng)濟(jì)發(fā)展持續(xù)攀升。山東省作為較早發(fā)展旅游業(yè)的省份，旅游經(jīng)濟(jì)發(fā)育日趨成熟，其中1999—2011年旅游經(jīng)濟(jì)收入增幅較大，增長(zhǎng)率大于20%，2012年后增幅放緩。由歷年地級(jí)市旅游總收入計(jì)算的變異系數(shù)可知：變異系數(shù)整體下降，在1995—2009年期間的波動(dòng)下降，其后降幅歸于平穩(wěn)。表明研究初期山東省各地級(jí)市旅游經(jīng)濟(jì)呈現(xiàn)顯著的空間差異性，其后逐步減弱，這是由于各地級(jí)市依托自身的旅游資源的優(yōu)勢(shì)，大力發(fā)展旅游業(yè)所致，同時(shí)各市不斷提高對(duì)旅游業(yè)的重視程度，地區(qū)間旅游經(jīng)濟(jì)發(fā)展更趨均衡。

圖1 1995—2018年山東省旅游經(jīng)濟(jì)發(fā)展?fàn)顩rFig. 1 Development of tourism economy in Shandong province from 1995 to 2018

2.2 旅游經(jīng)濟(jì)空間特征

本文探究了山東省各地級(jí)市1995—2018年旅游經(jīng)濟(jì)時(shí)空分布狀況，發(fā)現(xiàn)在各地區(qū)飛速發(fā)展的同時(shí)，各年份旅游經(jīng)濟(jì)發(fā)展水平的空間格局較為穩(wěn)定，高值區(qū)和低值區(qū)的基本格局變化較小，因此本文選取24年各地區(qū)旅游總收入均值進(jìn)行可視化(見(jiàn)圖2)。由圖2可知，山東省旅游經(jīng)濟(jì)發(fā)水平高值區(qū)為魯東地區(qū)的青島市、煙臺(tái)市和魯中地區(qū)的濟(jì)南市，空間上形成“雙核心”的格局。青島市作為山東省旅游經(jīng)濟(jì)發(fā)展的龍頭，其旅游經(jīng)濟(jì)總收入為位居第二的濟(jì)南市的1～2倍(見(jiàn)圖1)，近10年來(lái)與濟(jì)南市差距逐漸拉開(kāi)，在省內(nèi)形成逐步凸顯的極化格局。此外與青島市毗鄰的煙臺(tái)市旅游經(jīng)濟(jì)發(fā)展水平僅次于青島市和濟(jì)南市，青島市和煙臺(tái)市旅游經(jīng)濟(jì)總收入占全省旅游總收入的28%以上，二者成為省內(nèi)旅游經(jīng)濟(jì)的一大核心。在萊蕪市并入濟(jì)南市之前，濟(jì)南市旅游經(jīng)濟(jì)發(fā)展水平與煙臺(tái)市基本持平，并入萊蕪后濟(jì)南市成為僅次于青島市的省內(nèi)旅游經(jīng)濟(jì)發(fā)展的另一核心。旅游經(jīng)濟(jì)發(fā)展水平第二梯隊(duì)為濰坊市、泰安市、臨沂市、威海市、濟(jì)寧市和淄博市，在空間上連接成“軸”。泰安市、濟(jì)寧市等傳統(tǒng)旅游城市雖不及兩大核心，但也保持快速發(fā)展。同時(shí)，臨沂市、棗莊市、淄博市等新興旅游城市近些年來(lái)增長(zhǎng)速度也較快，其依托當(dāng)?shù)鬲?dú)特的自然及人文景觀大力發(fā)展旅游業(yè)，打造5A景區(qū)，提升景區(qū)規(guī)范化，不斷吸引了大量游客。同時(shí)德州市、濱州市、東營(yíng)市、日照市、聊城市、棗莊市和菏澤市，因其經(jīng)濟(jì)發(fā)展水平、旅游資源稟賦、交通區(qū)位條件和區(qū)域文化等方面的原因，旅游經(jīng)濟(jì)發(fā)展較緩慢。“兩核一軸”格局引領(lǐng)，周邊發(fā)展較弱，共同塑造出山東省旅游經(jīng)濟(jì)空間發(fā)展格局。

圖2 1995—2018年山東省各地級(jí)市旅游總收入均值Fig. 2 Average value of total tourism income of cities in Shandong province from 1995 to 2018

2.3 旅游經(jīng)濟(jì)空間格局演變的標(biāo)準(zhǔn)差橢圓分析

為了宏觀地了解研究時(shí)段山東省旅游經(jīng)濟(jì)空間格局特征，本文采用山東省旅游經(jīng)濟(jì)發(fā)展水平重心軌跡和標(biāo)準(zhǔn)差橢圓進(jìn)行分析。選取研究期內(nèi)的6個(gè)特征時(shí)點(diǎn)繪制山東省旅游經(jīng)濟(jì)空間分布格局演化趨勢(shì)，計(jì)算結(jié)果如圖3和表1所示。

圖3 山東省旅游經(jīng)濟(jì)空間分布格局演化趨勢(shì)Fig. 3 Evolution trend of spatial distribution pattern of tourism economy in Shandong province

表1 山東省旅游經(jīng)濟(jì)發(fā)展空間分布格局的標(biāo)準(zhǔn)差橢圓參數(shù)Tab. 1 Standard deviation ellipse parameter of spatial distribution pattern of tourism economic development in Shandong province

由結(jié)果可知，研究時(shí)段內(nèi)旅游經(jīng)濟(jì)發(fā)展重心均在濰坊市境內(nèi)，隨時(shí)間的推移，由東北向西南方向移動(dòng)，且移動(dòng)的角度和距離逐漸縮小?？傮w上，旅游經(jīng)濟(jì)發(fā)展重心逐年向山東省幾何重心靠攏。說(shuō)明西部地區(qū)旅游經(jīng)濟(jì)發(fā)展水平增長(zhǎng)速度高于東部地區(qū)，而在此期間東部地區(qū)也并非止步不前，因此2010年后重心移動(dòng)速率變慢。山東省旅游經(jīng)濟(jì)重心位于幾何重心的東北方向，究其原因，主要是東部沿海地區(qū)，特別是青島市和煙臺(tái)市，擁有優(yōu)越的自然環(huán)境條件與豐富的海洋旅游資源，經(jīng)濟(jì)發(fā)展水平明顯優(yōu)于西部地區(qū)，對(duì)游客有更強(qiáng)的吸引力。

研究期內(nèi)的標(biāo)準(zhǔn)差橢圓均呈現(xiàn)出長(zhǎng)軸遠(yuǎn)大于短軸，橢圓整體東西延伸成“梭形”趨勢(shì)。一方面，與山東省行政區(qū)劃東西向的空間形態(tài)有關(guān)；另一方面，也與山東省內(nèi)旅游經(jīng)濟(jì)“兩核一軸”的空間布局有關(guān)。此外，同旅游經(jīng)濟(jì)發(fā)展重心一樣，標(biāo)準(zhǔn)差橢圓隨時(shí)間推移由東北向西南方向移動(dòng)，且移動(dòng)幅度減小，橢圓逐步從“青—煙核”向“濟(jì)南核”移動(dòng)，山東省旅游大市青島市一直位于其內(nèi)部。研究期橢圓長(zhǎng)軸縮短，短軸增大，橢圓面積逐漸擴(kuò)大，橢圓扁率逐漸減小。這表明，近些年山東省旅游經(jīng)濟(jì)區(qū)域差異逐步減小，這一差異的縮小不但存在于東部與西部之間，也存在于南部和北部之間。在同一時(shí)期內(nèi)西部地區(qū)稍高的增長(zhǎng)率致使標(biāo)準(zhǔn)差橢圓逐漸西移。尤其是2000年以來(lái)，各地區(qū)對(duì)旅游的重視程度加強(qiáng)，西部地區(qū)的旅游業(yè)相對(duì)于以往落后的狀況有明顯改觀，因而在發(fā)展速度上表現(xiàn)為相對(duì)差距在縮小。

3 旅游經(jīng)濟(jì)的時(shí)空交互特征

3.1 LISA時(shí)間路徑幾何特征分析

通過(guò)記錄1995—2018年山東省16地市旅游經(jīng)濟(jì)發(fā)展水平Moran散點(diǎn)圖的坐標(biāo)，計(jì)算各市域單元LISA時(shí)間路徑的長(zhǎng)度和彎曲度，運(yùn)用自然斷裂點(diǎn)法將山東省旅游經(jīng)濟(jì)發(fā)展水平的長(zhǎng)度和彎曲度分別分為4個(gè)等級(jí)，可視化結(jié)果如圖4所示。

圖4 LISA時(shí)間路徑幾何特征空間分布Fig. 4 Spatial distribution of geometric characteristics of Lisa time path

研究期山東省市域旅游經(jīng)濟(jì)發(fā)展水平LISA時(shí)間路徑表明：相對(duì)長(zhǎng)度小于1的地級(jí)市共9個(gè)，占全部地級(jí)市的56.25%，表明山東省旅游經(jīng)濟(jì)發(fā)展水平的局部空間結(jié)構(gòu)具有一定的穩(wěn)定性，其中相對(duì)長(zhǎng)度值最高的為青島市(1.62)，其次為濰坊市(1.46)和煙臺(tái)市(1.44)，最低的地級(jí)市為棗莊市(0.55)，路徑長(zhǎng)度極差小，地區(qū)空間結(jié)構(gòu)動(dòng)態(tài)性差別小。從整體來(lái)看，LISA時(shí)間路徑大致表現(xiàn)出由東部沿海向西部?jī)?nèi)陸遞減，由中部向周圍遞減的趨勢(shì)，相對(duì)長(zhǎng)度與旅游經(jīng)濟(jì)發(fā)展水平的空間格局具有較強(qiáng)的相似度，尤其是低值的分布區(qū)與旅游經(jīng)濟(jì)低水平地區(qū)高度吻合。這表明1995—2018年魯東和魯中地區(qū)的旅游經(jīng)濟(jì)空間演化更具有動(dòng)態(tài)性，而魯北和魯西空間格局變化相對(duì)較小，相對(duì)于旅游經(jīng)濟(jì)較發(fā)達(dá)的地區(qū)，旅游經(jīng)濟(jì)較弱的地區(qū)空間結(jié)構(gòu)具有更強(qiáng)的穩(wěn)定性。這是由于研究期內(nèi)該地區(qū)旅游經(jīng)濟(jì)一直處于相對(duì)較低的水平，構(gòu)成相對(duì)穩(wěn)定的局部空間結(jié)構(gòu)。魯東和魯中地區(qū)長(zhǎng)期有效的區(qū)域間協(xié)作使旅游經(jīng)濟(jì)增速快，局部結(jié)構(gòu)不穩(wěn)定。

彎曲度小于均值的地級(jí)市占總體的93.75%，其中彎曲度最高的為煙臺(tái)市(41.14)，其值遠(yuǎn)超其他市域，表明煙臺(tái)市及其鄰域地級(jí)市在旅游經(jīng)濟(jì)增長(zhǎng)的過(guò)程中具有較大的波動(dòng)性。彎曲度最低的為臨沂市(2.56)，極差較大，不同地區(qū)在空間依賴上的波動(dòng)性差別巨大。從總體分布來(lái)看，低彎曲度地區(qū)分布廣泛，表明這些地級(jí)市對(duì)其臨近地區(qū)的空間依賴方向上具有較小的波動(dòng)性和較為穩(wěn)定的空間變遷過(guò)程，旅游經(jīng)濟(jì)相對(duì)穩(wěn)定；高彎曲度地級(jí)市有煙臺(tái)市、濟(jì)寧市等，這些地區(qū)與其鄰域間有明顯的動(dòng)態(tài)變遷和相互作用。

3.2 LISA時(shí)空躍遷分析

雖然LISA時(shí)間路徑分析可以有效地探究區(qū)域旅游經(jīng)濟(jì)發(fā)展的局部空間結(jié)構(gòu)的動(dòng)態(tài)性和波動(dòng)性特征，卻無(wú)法具體反映空間關(guān)聯(lián)類型中的轉(zhuǎn)移特征。因此，本文進(jìn)一步采用時(shí)空躍遷方法，進(jìn)一步刻畫(huà)市域旅游經(jīng)濟(jì)發(fā)展水平Moran散點(diǎn)圖在不同局部類型間的演化過(guò)程，構(gòu)建時(shí)空躍遷和轉(zhuǎn)移概率矩陣(表2)，并計(jì)算其空間穩(wěn)定性、時(shí)空變遷、時(shí)空凝聚和相對(duì)移動(dòng)率。

表2 局部 Moran’s I轉(zhuǎn)移概率矩陣Tab. 2 Local Moran’s I transition probability matrix

1995—2018年山東省旅游經(jīng)濟(jì)Ⅰ型轉(zhuǎn)移概率高達(dá)94.29%，相對(duì)移動(dòng)率為0.061 1，這反映出山東省旅游經(jīng)濟(jì)分布的路徑依賴或鎖定效應(yīng)極強(qiáng)，局部區(qū)域結(jié)構(gòu)穩(wěn)定。研究期內(nèi)從未發(fā)生時(shí)空躍遷的城市有10個(gè)(濟(jì)南市、東營(yíng)市、日照市、青島市、煙臺(tái)市、棗莊市、聊城市、菏澤市、德州市和濱州市)，占總數(shù)的62.5%，該市域單元多為旅游經(jīng)濟(jì)發(fā)展水平極高或者極低的地級(jí)市，其存在較強(qiáng)的轉(zhuǎn)移惰性，市域單元在Moran散點(diǎn)圖中難以改變自身的相對(duì)位置，局部區(qū)域結(jié)構(gòu)極為穩(wěn)定。研究期出現(xiàn)時(shí)空躍遷的地級(jí)市僅有6個(gè)(濟(jì)寧市、濰坊市、泰安市、威海市、淄博市和臨沂市)。這反映出山東省旅游經(jīng)濟(jì)空間凝聚趨勢(shì)非常強(qiáng)，具有高度的凝固性和較低的流動(dòng)性，這些旅游經(jīng)濟(jì)未發(fā)生躍遷的地級(jí)市，其穩(wěn)定性是限制山東省旅游經(jīng)濟(jì)發(fā)展水平整體躍遷的重要地市。研究期內(nèi)Ⅳ型躍遷最少，僅為2次，轉(zhuǎn)移概率僅為0.54%，因此Ⅱ型和Ⅲ型的時(shí)空躍遷地級(jí)市為驅(qū)動(dòng)山東省旅游經(jīng)濟(jì)時(shí)空躍遷的關(guān)鍵地級(jí)市，其時(shí)空變遷為5.16%。

根據(jù)分析可知，山東省旅游經(jīng)濟(jì)發(fā)展水平的局部空間結(jié)構(gòu)極為穩(wěn)定，市域單元存在極強(qiáng)的轉(zhuǎn)移惰性，僅有部分城市在某些年份存在較大的轉(zhuǎn)移概率。山東省旅游經(jīng)濟(jì)局部空間結(jié)構(gòu)受到自身旅游經(jīng)濟(jì)水平和鄰域單元旅游水平的共同影響，因此提升自身旅游經(jīng)濟(jì)水平并促進(jìn)同周邊區(qū)域間的協(xié)同發(fā)展是今后山東省旅游發(fā)展的主題。

3.3 時(shí)空網(wǎng)絡(luò)分析

本文通過(guò)繪制旅游經(jīng)濟(jì)時(shí)空網(wǎng)絡(luò)(圖5)來(lái)揭示臨接地級(jí)市的時(shí)間關(guān)聯(lián)強(qiáng)度大小。由結(jié)果可知，山東省旅游經(jīng)濟(jì)時(shí)間關(guān)聯(lián)強(qiáng)度以弱正向關(guān)聯(lián)和弱負(fù)向關(guān)聯(lián)為主，鄰域關(guān)聯(lián)共29對(duì)，其中弱正向關(guān)聯(lián)占比48.28%，弱負(fù)向關(guān)聯(lián)占比44.83%，強(qiáng)負(fù)向關(guān)聯(lián)占比6.89%，研究期內(nèi)無(wú)強(qiáng)正向關(guān)聯(lián)。區(qū)域間以弱協(xié)同發(fā)展關(guān)系和弱競(jìng)爭(zhēng)發(fā)展關(guān)系為主，含有少部分強(qiáng)競(jìng)爭(zhēng)發(fā)展關(guān)系，無(wú)強(qiáng)協(xié)同發(fā)展關(guān)系。濱州市和淄博市、菏澤市和濟(jì)寧市之間形成強(qiáng)負(fù)向關(guān)聯(lián)，表明地市之間存在極強(qiáng)的競(jìng)爭(zhēng)發(fā)展關(guān)系，區(qū)域旅游經(jīng)濟(jì)發(fā)展具有不協(xié)調(diào)特征。魯中、魯南和魯東地區(qū)內(nèi)旅游經(jīng)濟(jì)時(shí)間關(guān)聯(lián)強(qiáng)度以弱正向關(guān)聯(lián)為主，呈現(xiàn)較弱的正向協(xié)同的空間動(dòng)態(tài)性，這些地區(qū)大多旅游經(jīng)濟(jì)發(fā)展水平較高，特別是在2011年，山東半島藍(lán)色經(jīng)濟(jì)區(qū)正式上升為國(guó)家戰(zhàn)略，區(qū)域旅游產(chǎn)業(yè)也迎來(lái)了新的發(fā)展機(jī)遇。

圖5 山東省旅游經(jīng)濟(jì)的時(shí)空網(wǎng)絡(luò)Fig. 5 Spatio-temporal network of tourism economy in Shandong province

3.4 時(shí)空聚類分析

根據(jù)山東省地級(jí)市空間分布特征，選用回顧性時(shí)空重排掃描統(tǒng)計(jì)分析?？紤]到數(shù)據(jù)特征及數(shù)據(jù)運(yùn)算效果，本文將1995—2018年山東省各地市旅游總收入擴(kuò)大10倍，應(yīng)用SaTScan 9.5軟件進(jìn)行運(yùn)算。研究需要不同地區(qū)的地理坐標(biāo)數(shù)據(jù)、時(shí)間數(shù)據(jù)和旅游總收入數(shù)據(jù)。本文將最大時(shí)間周期設(shè)為研究時(shí)段的50%，最大聚類空間大小設(shè)置為50%，以1年為單位，蒙特卡羅模擬重復(fù)次數(shù)為999，顯著性檢驗(yàn)水平P<0.001，研究結(jié)果借助ArcGIS 10.2軟件進(jìn)行可視化表達(dá)(見(jiàn)圖6)。

圖6 1995—2018年山東省旅游經(jīng)濟(jì)時(shí)空掃描結(jié)果Fig. 6 Spatio-temporal scanning results of tourism economy in Shandong province from 1995 to 2018

研究期內(nèi)旅游經(jīng)濟(jì)并非隨機(jī)分布，而是呈現(xiàn)出2個(gè)明顯的時(shí)空聚類：聚類1以威海市為中心，范圍包括威海市、煙臺(tái)市和青島市，聚集時(shí)間在1996—2007年；聚類2以臨沂市為中心，范圍包括臨沂市、棗莊市和日照市，以及濟(jì)寧市、泰安市、濰坊市和淄博市的部分區(qū)域，聚集時(shí)間在2010—2018年。2個(gè)聚類P值均非常低，具有統(tǒng)計(jì)學(xué)意義，表明這2個(gè)聚類現(xiàn)象隨機(jī)出現(xiàn)的概率極低。聚類所含地級(jí)市均在“兩核一軸”范圍以內(nèi)，二者高度契合。兩大聚類輻射能力較強(qiáng)，范圍覆蓋山東省面積的一半以上，聚類區(qū)旅游經(jīng)濟(jì)發(fā)展水平良好，多為山東省旅游業(yè)發(fā)展的先行地市，旅游業(yè)發(fā)展程度較高，地區(qū)旅游特色和優(yōu)勢(shì)明顯。

4 結(jié)論與討論

4.1 結(jié)論

本文以市級(jí)層面作為研究尺度，選取旅游總收入作為衡量旅游經(jīng)濟(jì)發(fā)展水平的指標(biāo)，利用標(biāo)準(zhǔn)差橢圓、ESTDA和時(shí)空掃描統(tǒng)計(jì)量分析等方法，從時(shí)空交互視角探究1995—2018年山東省旅游經(jīng)濟(jì)的局部時(shí)空動(dòng)態(tài)特征，結(jié)果表明：

(1)研究期內(nèi)山東省旅游經(jīng)濟(jì)發(fā)展水平持續(xù)提升，在空間上形成“兩核一軸”的布局形態(tài)，其中一核為青島市和煙臺(tái)市，另一核為濟(jì)南市；隨著研究時(shí)間的推移，旅游經(jīng)濟(jì)發(fā)展重心和標(biāo)準(zhǔn)差橢圓均從東北向西南方向移動(dòng)，東西地區(qū)旅游經(jīng)濟(jì)差異逐漸縮小。因此，在未來(lái)發(fā)展中還應(yīng)繼續(xù)縮小地區(qū)差異，并提高旅游經(jīng)濟(jì)質(zhì)量，實(shí)現(xiàn)區(qū)域高質(zhì)量發(fā)展。

(2)山東省旅游經(jīng)濟(jì)發(fā)展水平的局部空間結(jié)構(gòu)具有一定的穩(wěn)定性，地區(qū)空間結(jié)構(gòu)動(dòng)態(tài)性差別小，魯東和魯中地區(qū)的旅游經(jīng)濟(jì)空間演化更具有動(dòng)態(tài)性，旅游經(jīng)濟(jì)較弱的地區(qū)空間結(jié)構(gòu)具有更強(qiáng)的穩(wěn)定性；不同地區(qū)在空間依賴上的波動(dòng)性差別巨大，波動(dòng)性較小的地級(jí)市廣泛分布。

(3)山東省旅游經(jīng)濟(jì)發(fā)展水平的局部空間結(jié)構(gòu)整體穩(wěn)定，市域單元存在極強(qiáng)的轉(zhuǎn)移惰性，相對(duì)位置變動(dòng)較困難，僅有部分城市在某些年份存在較大的轉(zhuǎn)移概率。各地級(jí)市應(yīng)進(jìn)一步提高當(dāng)?shù)芈糜谓?jīng)濟(jì)發(fā)展的積極性，尋求自身發(fā)展優(yōu)勢(shì)。在旅游經(jīng)濟(jì)發(fā)展過(guò)程中，區(qū)域間以弱協(xié)同發(fā)展關(guān)系和弱競(jìng)爭(zhēng)發(fā)展關(guān)系為主，魯中、魯南和魯東地區(qū)內(nèi)主要是弱正向關(guān)聯(lián)，在未來(lái)發(fā)展過(guò)程中應(yīng)注重區(qū)域間的協(xié)同發(fā)展。

(4)研究期內(nèi)旅游經(jīng)濟(jì)時(shí)空上呈現(xiàn)出2個(gè)明顯的時(shí)空聚類，覆蓋范圍廣泛，且各地區(qū)旅游特征明顯。

4.2 討論

2018年山東省委、省政府印發(fā)《大力推進(jìn)全域旅游高質(zhì)量發(fā)展實(shí)施方案》，表明在今后，旅游業(yè)越發(fā)成為省內(nèi)戰(zhàn)略性支柱產(chǎn)業(yè)。雖然在研究中發(fā)現(xiàn)，山東省旅游經(jīng)濟(jì)發(fā)展水平持續(xù)提升的同時(shí)地區(qū)之間差異縮小，但整體來(lái)看，地區(qū)差異仍不容樂(lè)觀。自1995年以來(lái)，省內(nèi)旅游經(jīng)濟(jì)高值區(qū)穩(wěn)定于“兩核一軸”地帶，如何打破“兩核一軸”格局，實(shí)現(xiàn)旅游經(jīng)濟(jì)高質(zhì)量發(fā)展；如何發(fā)揮核心地級(jí)市旅游經(jīng)濟(jì)對(duì)周邊的極核輻射作用，以軸線串聯(lián)起各地級(jí)市來(lái)實(shí)現(xiàn)區(qū)域間協(xié)同發(fā)展等問(wèn)題，將是今后山東省旅游經(jīng)濟(jì)發(fā)展面臨的新挑戰(zhàn)?；谏衔慕Y(jié)論及分析，筆者對(duì)山東省旅游經(jīng)濟(jì)發(fā)展提出如下建議：①各地級(jí)市應(yīng)合理規(guī)劃旅游產(chǎn)業(yè)規(guī)模和結(jié)構(gòu)，各旅游部門(mén)加強(qiáng)協(xié)同聯(lián)動(dòng)，因地制宜，發(fā)揮自身旅游資源優(yōu)勢(shì)，杜絕旅游資源千篇一律現(xiàn)象；②提升旅游資源轉(zhuǎn)換效率，提高單位旅游資源投入獲得的經(jīng)濟(jì)和生態(tài)效益，實(shí)現(xiàn)地區(qū)旅游經(jīng)濟(jì)高質(zhì)量發(fā)展；③各地級(jí)市應(yīng)提高旅游經(jīng)濟(jì)增長(zhǎng)質(zhì)量、旅游產(chǎn)品質(zhì)量和旅游環(huán)境質(zhì)量，堅(jiān)持“綠水青山就是金山銀山”的生態(tài)文明建設(shè)思想，以“兩山”理論為指導(dǎo)，實(shí)現(xiàn)旅游資源的永續(xù)利用和區(qū)域可持續(xù)發(fā)展；④優(yōu)化旅游資源配置，增加旅游資源科技含量，實(shí)現(xiàn)創(chuàng)新驅(qū)動(dòng)旅游經(jīng)濟(jì)高質(zhì)量發(fā)展。

此外，本研究仍有諸多不足，由于研究時(shí)段較長(zhǎng)，所獲取數(shù)據(jù)非常有限，研究中還缺少旅游經(jīng)濟(jì)增長(zhǎng)和旅游經(jīng)濟(jì)局域空間格局的穩(wěn)定與動(dòng)態(tài)性影響因素的探討，因此在今后的研究中還應(yīng)選取相應(yīng)的影響因素，進(jìn)一步探究影響旅游經(jīng)濟(jì)發(fā)展的驅(qū)動(dòng)機(jī)制。