王 璇，史佳璐，慈福義

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

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

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

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

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

1.1 指標選取和數(shù)據(jù)來源

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

1.2 研究方法

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

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

1)LISA時間路徑

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

(1)

(2)

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

2)LISA時空躍遷

Rey等將Moran’s Ⅰ散點圖中各空間單元在特定時間間隔內的移動距離等屬性嵌入傳統(tǒng)馬爾可夫鏈中，提出轉移概率矩陣和時空躍遷理論[35]。本文運用時空躍遷方法測度各地級市旅游經(jīng)濟的空間關聯(lián)類型的動態(tài)變遷過程，根據(jù)各區(qū)域自身與鄰域單元之間的旅游經(jīng)濟發(fā)展水平轉移狀態(tài)將時空躍遷劃分為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īng)濟的Moran’s Ⅰ散點圖空間穩(wěn)定性、時空變遷、時空凝聚和相對移動率依據(jù)可能存在的躍遷數(shù)和各躍遷類型計算，可能存在的躍遷總數(shù)n=(2018-1995)×16=368，具體計算公式詳見參考文獻[41,42]。

3)時空交互可視化

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

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

2 旅游經(jīng)濟的時空發(fā)展狀況

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

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

圖1 1995—2018年山東省旅游經(jīng)濟發(fā)展狀況Fig. 1 Development of tourism economy in Shandong province from 1995 to 2018

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

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

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

2.3 旅游經(jīng)濟空間格局演變的標準差橢圓分析

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

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

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

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

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

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

3.1 LISA時間路徑幾何特征分析

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

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

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

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

3.2 LISA時空躍遷分析

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

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

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

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

3.3 時空網(wǎng)絡分析

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

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

3.4 時空聚類分析

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

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

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

4 結論與討論

4.1 結論

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

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

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

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

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

4.2 討論

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

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