李云濤 ，王海英 ，陶 犁 ※，白廷斌 ，3

（1.首都師范大學(xué)資源環(huán)境與旅游學(xué)院，北京 100048；2.首都師范大學(xué)首都旅游研究院，北京 100048；3.云南財經(jīng)大學(xué)經(jīng)濟(jì)研究院，昆明 650221）

0 引言

鄉(xiāng)村旅游是保護(hù)鄉(xiāng)村生態(tài)環(huán)境，改善鄉(xiāng)村非農(nóng)就業(yè)結(jié)構(gòu)，有效推進(jìn)鄉(xiāng)村一二三產(chǎn)業(yè)融合發(fā)展的重要途徑之一，是推進(jìn)城鄉(xiāng)一體化建設(shè)、鄉(xiāng)村振興戰(zhàn)略的重要抓手之一。鄉(xiāng)村旅游地作為鄉(xiāng)村旅游的重要組成部分和鄉(xiāng)村旅游活動的依托載體，其空間分布情況可以反映一個地區(qū)的鄉(xiāng)村旅游發(fā)展水平。同時，鄉(xiāng)村旅游發(fā)展過程中諸多問題的解決，需建立在正確分析鄉(xiāng)村旅游地空間分布及其影響因素的基礎(chǔ)上。

鄉(xiāng)村旅游地的空間分布演化及其影響因素研究受到國內(nèi)外學(xué)者關(guān)注。從研究內(nèi)容上看，國外學(xué)者從鄉(xiāng)村旅游地的空間結(jié)構(gòu)特征[1]、空間演化規(guī)律[2]、影響因素[3,4]、驅(qū)動機(jī)制[5]、開發(fā)經(jīng)營模式[6]等方面進(jìn)行了研究。國內(nèi)相關(guān)研究相較于國外起步較晚，主要圍繞鄉(xiāng)村旅游地的景點類型劃分[7]、空間分布特征[8]、空間聚集特征[9,10]、影響因素[11]、空間分布與時空演變機(jī)制研究[12,13]、發(fā)展模式[14,15]、綜合評價研究[16,17]等方面。從研究方法上看，國內(nèi)外學(xué)者都逐漸由定性研究轉(zhuǎn)向注重定量與定性相結(jié)合的研究方法，定量方法主要使用了非期望產(chǎn)出模型[18]、空間點格局分析法[19]、地理集中指數(shù)[20]、空間基尼系數(shù)[21]、地理聯(lián)系率[22]、緩沖區(qū)分析[23]、全局空間自相關(guān)方法[24]、格網(wǎng)分析法[25]等數(shù)理模型和空間分析方法。從研究的空間尺度上涵蓋了全國[26,27]、區(qū)域[28,29]、省域[30,31]、市級[32]，也有從微觀視角出發(fā)[33,34]，通過選取具有典型性和代表性的鄉(xiāng)村旅游案例地進(jìn)行剖析。我國鄉(xiāng)村旅游研究起步較晚，加之2013年以前旅游主管部門對于鄉(xiāng)村旅游的統(tǒng)計數(shù)據(jù)較少，這使得多數(shù)宏觀中觀尺度的研究在數(shù)據(jù)獲取方面就遭遇瓶頸，近年來雖有增多的趨勢，但對影響因素的研究以單一定性分析為主。此外，目前對于鄉(xiāng)村旅游地空間分布的研究多為靜態(tài)分析，動態(tài)演化分析的相對較少。

北京鄉(xiāng)村旅游興起于20世紀(jì)90年代后期，經(jīng)過30余年的積極探索，形成了在區(qū)域特色、行業(yè)管理、業(yè)態(tài)創(chuàng)新等方面極具特色的鄉(xiāng)村旅游“北京模式”，為中國鄉(xiāng)村旅游的發(fā)展探索出了一條獨特的發(fā)展道路。截止2020年8月，北京市共有32個全國重點鄉(xiāng)村旅游村、34個中國美麗休閑鄉(xiāng)村、21個國家休閑農(nóng)業(yè)與鄉(xiāng)村旅游示范點、274個星級民俗村、6 042個星級民俗接待戶[35]。北京作為我國鄉(xiāng)村旅游發(fā)展較早且較成熟的地區(qū)，使用定量分析方法對北京郊區(qū)鄉(xiāng)村旅游地空間分布演化及其影響因素進(jìn)行的研究尚不多見。且近年來，北京不斷完善鄉(xiāng)村旅游評定標(biāo)準(zhǔn)，加之全市經(jīng)濟(jì)社會的快速發(fā)展，都將對京郊鄉(xiāng)村旅游地的空間分布格局產(chǎn)生重要影響。在此背景下，分析京郊鄉(xiāng)村旅游地的空間分布演化特征及其影響因素，對促進(jìn)北京鄉(xiāng)村旅游發(fā)展具有重要的現(xiàn)實意義。

鑒于此，文章以2010—2020年的京郊鄉(xiāng)村旅游特色業(yè)態(tài)和鄉(xiāng)村民俗旅游村為研究對象，綜合運用計量地理學(xué)模型和ArcGIS 10.2空間分析工具，分析京郊鄉(xiāng)村旅游地的空間分布特征及演化過程，定量剖析其影響因素，旨在探索京郊鄉(xiāng)村旅游地空間分布演化特征，明晰其影響因素，為京郊鄉(xiāng)村旅游地空間結(jié)構(gòu)的優(yōu)化，促進(jìn)京郊鄉(xiāng)村旅游的健康發(fā)展提供科學(xué)依據(jù)。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

北京郊區(qū)地處華北平原西北邊緣，西部為西山屬太行山脈，北部和東北部為軍都山屬燕山山脈，四面與河北省接壤，東南面與天津市接壤，具體包含了門頭溝區(qū)、房山區(qū)、通州區(qū)等10個北京市下轄的區(qū)級行政單元（圖1）。京郊地區(qū)擁有豐富的自然資源和獨特的歷史人文景觀，不僅是拱衛(wèi)北京中心城區(qū)的綠色生態(tài)屏障，也是北京市民休閑度假、親近自然、體驗民俗文化的重要旅游目的地。鄉(xiāng)村旅游是京郊旅游的重要組成部分，京郊鄉(xiāng)村旅游憑借鄉(xiāng)村特色產(chǎn)業(yè)、優(yōu)美的自然風(fēng)光和獨特的民俗文化，逐漸形成了包括民俗旅游村、休閑度假村、觀光農(nóng)業(yè)示范園等眾多極具特色的京郊鄉(xiāng)村旅游目的地。

圖1 北京郊區(qū)區(qū)位

1.2 研究方法

1.2.1 標(biāo)準(zhǔn)差橢圓

標(biāo)準(zhǔn)差橢圓法是分析空間分布方向性特征的經(jīng)典方法之一[36]，可從中心性、方向性、空間形態(tài)等方面，直觀地刻畫出研究對象在區(qū)域空間分布上的整體性特征。該文借助標(biāo)準(zhǔn)差橢圓，以揭示京郊鄉(xiāng)村旅游地空間分布的時空演變特征。標(biāo)準(zhǔn)差橢圓法主要包括中心點、長半軸、短半軸和方位旋轉(zhuǎn)角等基本參數(shù)。其中，方位旋轉(zhuǎn)角表示京郊鄉(xiāng)村旅游地分布的主趨勢方向；長半軸表示京郊鄉(xiāng)村旅游地分布在主趨勢方向的離散程度；短半軸表示其在次趨勢方向的離散程度。具體測算過程參考文獻(xiàn)[37]的研究。

1.2.2 重心坐標(biāo)

重心坐標(biāo)表示某一區(qū)域范圍內(nèi)，某種屬性值在空間內(nèi)各個方向上能夠保持均衡的點的坐標(biāo)[38]。通過觀察重心的空間變遷，可以了解該種屬性值在區(qū)域范圍內(nèi)的發(fā)展方向和平衡問題。假設(shè)某一較大區(qū)域由n個小區(qū)域構(gòu)成，第i個小區(qū)域的中心坐標(biāo)為Ni(xi，yi)，ui為該小區(qū)域的某種屬性值，則該種屬性值意義下的區(qū)域重心坐標(biāo)為)，計算公式為：

1.2.3 核密度估算法

核密度估算法是一種非參數(shù)估計方法，即借助一個移動的單元格或窗口，對給定區(qū)域范圍內(nèi)的點或線格局的密度進(jìn)行估計[39]。通過核密度估算法測算所得的核密度值及其可視化后的形態(tài)特征，能清晰反映出京郊鄉(xiāng)村旅游地在空間上的集聚或離散特征及時空過程。核密度估算公式為：

式（2）中，Kn(x)為京郊鄉(xiāng)村旅游地空間分布的核密度估算值；n為京郊范圍內(nèi)的鄉(xiāng)村旅游地數(shù)量；h（h＞0）為指定帶寬的平滑參數(shù)；為核函數(shù)；k為權(quán)函數(shù)，表示鄉(xiāng)村旅游地x的權(quán)重，若x與中心點的距離越大，其權(quán)重越??；x-xi表示鄉(xiāng)村旅游地x到鄉(xiāng)村旅游地xi的距離。

1.2.4 地理集中指數(shù)

地理集中指數(shù)是某一節(jié)點要素在研究區(qū)域范圍內(nèi)各子區(qū)域的空間集中分布程度，是衡量研究對象區(qū)域分布集中程度的重要指標(biāo)[40]，其計算公式為：

式（3）中，G為京郊鄉(xiāng)村旅游地的地理集中指數(shù)；xi為京郊第i個區(qū)縣的鄉(xiāng)村旅游地的數(shù)量；t為京郊鄉(xiāng)村旅游地的總數(shù)；n為北京郊區(qū)區(qū)縣的數(shù)量。G的取值范圍為[0, 100]，G值越大，表示京郊鄉(xiāng)村旅游地的分布越集中。

1.2.5 基尼系數(shù)

基尼系數(shù)是經(jīng)濟(jì)學(xué)中用于衡量一個國家或地區(qū)居民收入均衡程度的重要指標(biāo)，地理學(xué)中運用該方法刻畫某一地理要素的空間分布離散程度[41]，其計算公式為：

式（4）（5）中，G為京郊鄉(xiāng)村旅游地的基尼系數(shù)；pi為京郊第i個區(qū)縣的鄉(xiāng)村旅游地數(shù)量占京郊鄉(xiāng)村旅游地總數(shù)的比重；N為北京郊區(qū)區(qū)縣的數(shù)量。G的取值范圍為[0, 1]，G值越大，表示京郊鄉(xiāng)村旅游地的分布越集中。C為均勻分布度，G值越大，則C值越小，表示均勻度越低。

1.2.6 地理探測器

地理探測器是探測空間分異性，揭示其背后驅(qū)動因子的一種計量分析方法，因其具有無線性假設(shè)、共線性免疫、使用靈活等優(yōu)點，被廣泛用于探測分析地理對象空間異質(zhì)性的影響因素和作用機(jī)理[42]。具體做法是比較該影響因子在不同分區(qū)內(nèi)的區(qū)內(nèi)方差和全區(qū)的總方差，其計算模型為：

式（6）（7）中，q為某個影響因子對某一地理現(xiàn)象空間分異的影響作用能力的探測值；h=1，2，…；L為自變量X和因變量Y的分層；Nh和N分別為層h和全區(qū)的單元數(shù)量；σ2h和σ2分別為層H和全區(qū)的因變量Y的方差。SSW和SST分別為層內(nèi)方差之和與全區(qū)總方差。q的取值范圍為[0, 1]，q值越大，表示影響因素X對京郊鄉(xiāng)村旅游地時空分布影響作用力就越大。

1.3 數(shù)據(jù)來源與數(shù)據(jù)預(yù)處理

該研究以京郊10區(qū)范圍內(nèi)的鄉(xiāng)村旅游特色業(yè)態(tài)、鄉(xiāng)村民俗旅游村為研究對象，主要來源于2010—2020年北京市旅游主管部門公布的北京市鄉(xiāng)村旅游特色業(yè)態(tài)名錄和鄉(xiāng)村民俗旅游村名錄，共1 063家。經(jīng)濟(jì)、社會發(fā)展等數(shù)據(jù)主要來源于京郊10區(qū)政府網(wǎng)站發(fā)布的2010—2020年的國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報、政府工作報告，部分補(bǔ)充數(shù)據(jù)來源于2010—2020年的《北京區(qū)域統(tǒng)計年鑒》。北京郊區(qū)矢量數(shù)據(jù)、河流數(shù)據(jù)從全國地理信息資源目錄服務(wù)系統(tǒng)獲取（https://www.webmap.cn）；數(shù)字高程模型（DEM）數(shù)據(jù)來源于地理空間數(shù)據(jù)云（http://www.gscloud.cn/）；氣溫、日照時數(shù)、風(fēng)速、空氣相對濕度等氣象觀測月均值數(shù)據(jù)來源于國家氣象科學(xué)數(shù)據(jù)中心（http://data.cma.cn/）北京市郊區(qū)的14個基準(zhǔn)臺站，時間范圍為2010—2020年；空氣質(zhì)量等級數(shù)據(jù)主要來源于中國環(huán)境監(jiān)測總站的歷史發(fā)布值（http://www.cnemc.cn/），時間范圍為2010—2020年。

數(shù)據(jù)預(yù)處理方面，借助“百度地圖坐標(biāo)拾取系統(tǒng)”，通過網(wǎng)絡(luò)爬蟲獲取京郊鄉(xiāng)村旅游特色業(yè)態(tài)、鄉(xiāng)村民俗旅游村的經(jīng)緯度坐標(biāo)，共1 063處，使用ArcGIS10.2軟件的坐標(biāo)轉(zhuǎn)點工具，將其轉(zhuǎn)化為點要素矢量圖層。對氣象觀測數(shù)據(jù)的處理，使用ArcGIS10.2軟件的地統(tǒng)計分析模塊選用樣條函數(shù)插值法對各項氣象數(shù)據(jù)進(jìn)行空間插值。

2 京郊鄉(xiāng)村旅游地空間分布演化

2.1 京郊鄉(xiāng)村旅游地空間分布范圍特征

借助ArcGIS10.2軟件繪制2010年、2013年、2016年、2019年和2020年的京郊鄉(xiāng)村旅游地空間分布重心、標(biāo)準(zhǔn)差橢圓，直觀地反映出京郊鄉(xiāng)村旅游地的主要空間分布范圍演化趨勢，如圖2、表1。

圖2 2010年、2013年、2016年、2019年和2020年京郊鄉(xiāng)村旅游地重心及標(biāo)準(zhǔn)差橢圓分布

表1 2010年、2013年、2016年、2019年和2020年京郊鄉(xiāng)村旅游地標(biāo)準(zhǔn)差橢圓參數(shù)

從重心軌跡的計算結(jié)果可以看出，2010年、2013年、2016年、2019年和2020年的京郊鄉(xiāng)村旅游地重心分布所在區(qū)域較為穩(wěn)定，均位于昌平區(qū)內(nèi)，這表明京郊鄉(xiāng)村旅游地密度較大區(qū)域在昌平區(qū)。重心移動軌跡為：昌平區(qū)內(nèi)，2010—2013年重心向東南方向遷移，2013—2016年向東北方向遷移后，重心總體保持在東北方位。

（1）從標(biāo)準(zhǔn)差橢圓計算結(jié)果來看，標(biāo)準(zhǔn)差橢圓的形態(tài)變化幅度不大，方位與北京郊區(qū)的形狀、方向大致相似。2010年的標(biāo)準(zhǔn)差橢圓位置差別相對較為顯著，偏向于西南方向，其余年份呈小幅變化態(tài)勢，并未出現(xiàn)橢圓方向、面積較大變化差異的情況。

（2）標(biāo)準(zhǔn)差橢圓面積逐漸減小，空間分布越來越集中。從橢圓的面積變化來看，2010年的標(biāo)準(zhǔn)差橢圓面積最大，為14 465.48km2，此時京郊鄉(xiāng)村旅游地空間分布最為分散。2013年橢圓面積開始逐年下降，2020年橢圓面積略微增大。整體上看，橢圓面積主要呈減小趨勢，表示京郊鄉(xiāng)村旅游地主要分布范圍逐漸由分散趨向于集中。

（3）標(biāo)準(zhǔn)差橢圓長軸長度呈“減小—增大”的過程特征，由2010年的87.73km逐漸減少為2016年的78.38km，2019年和2020年小幅增大，可看出京郊鄉(xiāng)村旅游地在南北分布主方向上呈現(xiàn)“分散—集中”的變化過程。標(biāo)準(zhǔn)差橢圓短軸長度呈小幅波動，但總體變化不明顯，2010年的短軸長度最小，為52.49km，2016年后維持小幅度波動變化，這表明京郊鄉(xiāng)村旅游地在東西分布次方向上的空間變化并不明顯。

2.2 京郊鄉(xiāng)村旅游地分布集中程度

借助地理集中指數(shù)來衡量京郊鄉(xiāng)村旅游地在京郊各區(qū)的空間分布集中程度。根據(jù)地理集中指數(shù)模型，計算2010—2020年京郊鄉(xiāng)村旅游地的實際分布的地理集中指數(shù)和均勻分布的地理集中指數(shù)。假設(shè)第i年京郊鄉(xiāng)村旅游地數(shù)量為n個，且平均分布于京郊的10個區(qū)內(nèi)，則均勻分布的地理集中指數(shù)計算公式為Zi=n/10，計算結(jié)果見表2。

根據(jù)表2的計算結(jié)果可以看出：2010—2012年實際分布的地理集中指數(shù)G略大于均勻分布的地理集中指數(shù)Z，表明京郊鄉(xiāng)村旅游地在各區(qū)分布較為集中，且集中程度變化較小。2013—2020年實際分布的地理集中指數(shù)G則均小于均勻分布的地理集中指數(shù)Z，且G值總體呈逐年減小趨勢，這表明該時期內(nèi)，京郊區(qū)級層面的鄉(xiāng)村旅游地分布集中程度呈連續(xù)小幅下降。總體而言，2010—2020年京郊鄉(xiāng)村旅游地在區(qū)級層面的分布是不均衡的，且呈逐年小幅增加態(tài)勢。

表2 2010—2020年京郊鄉(xiāng)村旅游地地理集中指數(shù)

2.3 京郊鄉(xiāng)村旅游地分布均衡程度

借助基尼系數(shù)可以衡量京郊鄉(xiāng)村旅游地在京郊各區(qū)內(nèi)部的分布均勻程度。根據(jù)基尼系數(shù)和均勻分布度計算模型，計算2010—2020年京郊鄉(xiāng)村旅游地的基尼系數(shù)和均勻分布度，計算結(jié)果見圖3。

根據(jù)圖3可以看出，2010—2020年基尼系數(shù)均保持較高水平，維持在0.92以上，除2011年有小幅下降，總體呈逐年上升趨勢；均勻分布度總體保持較低水平，除2011年小幅上升，其余年份呈逐年下降趨勢。這表明2010—2020年京郊鄉(xiāng)村旅游地在京郊各區(qū)內(nèi)部的空間分布均維持極不均衡狀態(tài)。

圖3 2010—2020年京郊鄉(xiāng)村旅游地基尼系數(shù)和均勻分布度

2.4 京郊鄉(xiāng)村旅游地密度分布特征

借助ArcGIS 10.2的核密度分析工具，計算2010年、2013年、2016年、2019年和2020年京郊鄉(xiāng)村旅游地的空間分布核密度，對其進(jìn)行可視化呈現(xiàn)（圖4）。由圖4可以看出，京郊鄉(xiāng)村旅游地的空間分布核密度變化較為明顯，呈現(xiàn)出一定的階段性特征，具體如下。

圖4 2010年、2013年、2016年、2019年和2020年京郊鄉(xiāng)村旅游地空間分布核密度

（1）2010—2013年高密度區(qū)數(shù)量逐漸增加，呈現(xiàn)以點狀輻射分布為主的特征。2010年京郊鄉(xiāng)村旅游地已初步形成了位于昌平和密云的兩個高密度核心區(qū)。2011—2013年京郊鄉(xiāng)村旅游地相比2010年時數(shù)量明顯增加，形成多個高密度區(qū)，其中，昌平的高密度核向北延伸拓展，密云的高密度核轉(zhuǎn)化為高密度區(qū)。這些高密度區(qū)奠定了之后京郊鄉(xiāng)村旅游地密度變化的主要區(qū)域，具體主要集中在昌平、密云、懷柔、延慶以及大興。

（2）2014—2016年高密度區(qū)呈連片拓展特征，形成昌平、懷柔雙高密度核心區(qū)。2014—2016年京郊鄉(xiāng)村旅游地數(shù)量規(guī)模迅速擴(kuò)大，高密度區(qū)向周邊擴(kuò)展延伸，在延慶、密云、大興形成多個高密度連片區(qū)域，懷柔形成高密度核心區(qū)并向東輻射拓展。

（3）2017—2020年形成以昌平、懷柔為雙高密度核心區(qū)，密云、延慶、大興多個高密度連片區(qū)并存的格局。2017年以來，昌平、懷柔的高密度核心區(qū)核持續(xù)擴(kuò)張，昌平的高密度核心區(qū)向南向擴(kuò)展，懷柔高密度核持續(xù)向東擴(kuò)展。密云、延慶、大興三區(qū)的高密度區(qū)形成較大范圍的連片區(qū)域，門頭溝、房山、通州、順義連片區(qū)域較少。

3 京郊鄉(xiāng)村旅游地空間分布的影響因素研究

3.1 影響因素的選取及數(shù)據(jù)處理

3.1.1 影響因素的選取

很多學(xué)者對不同地區(qū)鄉(xiāng)村旅游地空間分布演化的影響因素進(jìn)行了探討，見表3，可以看出，鄉(xiāng)村旅游地的形成與發(fā)展不僅需借勢于鄉(xiāng)村自身的自然生態(tài)環(huán)境條件，同時經(jīng)濟(jì)、文化、政治、交通條件等多個因素也驅(qū)動了鄉(xiāng)村旅游地的空間分布演化。

表3 鄉(xiāng)村旅游地空間分布的影響因素文獻(xiàn)總結(jié)

參考相關(guān)已有研究成果，結(jié)合京郊鄉(xiāng)村旅游發(fā)展和鄉(xiāng)村旅游地空間分布的實際情況，該文認(rèn)為，京郊鄉(xiāng)村旅游地的空間分布演化是多種因素共同作用，總體上可將影響因素綜合為自然環(huán)境因素、產(chǎn)業(yè)經(jīng)濟(jì)因素和社會政治因素三大類。

自然環(huán)境因素是鄉(xiāng)村旅游開發(fā)和建設(shè)的基礎(chǔ)，海拔低、地勢起伏緩和、臨近河流、生態(tài)環(huán)境優(yōu)美的鄉(xiāng)村，往往易于形成鄉(xiāng)村旅游地。產(chǎn)業(yè)經(jīng)濟(jì)因素包括了國民經(jīng)濟(jì)發(fā)展?fàn)顩r、旅游產(chǎn)業(yè)發(fā)展?fàn)顩r、農(nóng)業(yè)發(fā)展?fàn)顩r等，產(chǎn)業(yè)經(jīng)濟(jì)發(fā)展?fàn)顩r好的地區(qū)，經(jīng)濟(jì)發(fā)展水高，產(chǎn)業(yè)結(jié)構(gòu)合理，旅游產(chǎn)業(yè)發(fā)達(dá)，推動著鄉(xiāng)村旅游的形成與發(fā)展，從而影響鄉(xiāng)村旅游地的形成與分布。社會政治因素主要包括旅游政策導(dǎo)向、交通基礎(chǔ)設(shè)施、旅游地常住人口等，旅游政策指導(dǎo)鄉(xiāng)村旅游地的發(fā)展方向，較好的交通基礎(chǔ)設(shè)施提升鄉(xiāng)村旅游地的可達(dá)性，旅游地常住人口為鄉(xiāng)村旅游地提供了客源，支撐著鄉(xiāng)村旅游地的形成與發(fā)展?；诖?，該文構(gòu)建京郊鄉(xiāng)村旅游地空間分布影響指標(biāo)（表4）。

表4 京郊鄉(xiāng)村旅游地空間分布影響因素

（1）自然環(huán)境因素：海拔高度一定程度上影響著鄉(xiāng)村及鄉(xiāng)村旅游資源的形成與分布，從而影響鄉(xiāng)村旅游地的形成。河流為村莊提供水源保障，河流沿岸土壤肥沃易于形成農(nóng)業(yè)種植生態(tài)景觀，同時河流水系也可作為鄉(xiāng)村旅游觀賞景點。氣候舒適的地區(qū)有利于農(nóng)業(yè)生產(chǎn)活動的開展，而游客對氣候體感舒適度也影響著游客的旅游意愿。良好的空氣質(zhì)量可提升游客在鄉(xiāng)村旅游過程中的旅游體驗，空氣質(zhì)量較差，尤其霧霾天氣的出現(xiàn)，會極大影響游客體驗，損害鄉(xiāng)村旅游地形象。林木綠化率是衡量地區(qū)自然生態(tài)環(huán)境質(zhì)量的重要指標(biāo)，良好生態(tài)環(huán)境的鄉(xiāng)村往往會受到游客的青睞，較易發(fā)展成為鄉(xiāng)村旅游地。

（2）經(jīng)濟(jì)產(chǎn)業(yè)因素：地區(qū)生產(chǎn)總值反映了地區(qū)整體經(jīng)濟(jì)發(fā)展水平，經(jīng)濟(jì)發(fā)展水平高的地區(qū)，對鄉(xiāng)村旅游的發(fā)展起到重要的推動作用，進(jìn)而影響鄉(xiāng)村旅游地的形成與空間分布。第三產(chǎn)業(yè)占比反映地區(qū)的第三產(chǎn)業(yè)發(fā)展水平及產(chǎn)業(yè)結(jié)構(gòu)狀況，鄉(xiāng)村旅游作為第三產(chǎn)業(yè)的新興行業(yè)，第三產(chǎn)業(yè)的發(fā)展?fàn)顩r勢必會影響到鄉(xiāng)村旅游地的形成與發(fā)展。旅游收入可以體現(xiàn)地區(qū)旅游業(yè)發(fā)展水平，鄉(xiāng)村旅游的發(fā)展離不開當(dāng)?shù)乜傮w旅游業(yè)的支撐，當(dāng)?shù)芈糜萎a(chǎn)業(yè)越發(fā)達(dá)，對鄉(xiāng)村旅游發(fā)展的支撐力度就越強(qiáng)。A級旅游景區(qū)擁有較好的旅游品牌和豐富的旅游資源，鄉(xiāng)村旅游地可依附于景區(qū)實現(xiàn)聯(lián)動發(fā)展。農(nóng)林牧漁業(yè)總產(chǎn)值反映了地區(qū)農(nóng)業(yè)發(fā)展水平，農(nóng)業(yè)發(fā)達(dá)的地區(qū)，往往鄉(xiāng)村發(fā)展程度較高，可拓展多種發(fā)展業(yè)態(tài)，易于形成鄉(xiāng)村旅游地。

（3）社會政治因素：交通是鄉(xiāng)村旅游地形成和發(fā)展的基礎(chǔ)，發(fā)揮著連接旅游客源地和旅游目的地橋梁的作用，是鄉(xiāng)村旅游系統(tǒng)中非常重要的組成部分。本地城鎮(zhèn)居民的可支配收入水平影響著他們的旅游出行意愿和旅游消費能力，從而影響到鄉(xiāng)村旅游地的發(fā)展。本地的常住城鎮(zhèn)居民是鄉(xiāng)村旅游的主要客源市場，常住城鎮(zhèn)人口較多，潛在的鄉(xiāng)村旅游客源規(guī)模便越大。政府政策為鄉(xiāng)村旅游的發(fā)展提供了支持與引導(dǎo)調(diào)控，決定著鄉(xiāng)村旅游地的發(fā)展方向、規(guī)模和空間布局。農(nóng)林水事務(wù)支出反映了政府財政用于支持農(nóng)業(yè)、農(nóng)村發(fā)展的資金投入情況，積極的財政扶持可加快鄉(xiāng)村旅游地的發(fā)展，但財政的地區(qū)傾斜會影響鄉(xiāng)村旅游地的形成與分布。

3.1.2 數(shù)據(jù)處理

數(shù)據(jù)的定量化處理。如表5所示，對部分不能直接獲取量化數(shù)據(jù)的影響因素進(jìn)行量化處理。

表5 京郊鄉(xiāng)村旅游地空間分布二級影響因素數(shù)據(jù)等級類型劃分

（1）對自然環(huán)境因素的定量化處理。海拔高度（X1）使用ArcGIS 10.2軟件從京郊的DEM數(shù)據(jù)中提取。河流水系（X2）借助ArcGIS 10.2軟件對京郊河流水系構(gòu)建0～2km范圍以內(nèi)的4級河流緩沖區(qū)。氣候舒適度（X3）參考文獻(xiàn)[50]的研究，使用溫濕指數(shù)（THI）和風(fēng)寒指數(shù)（K）的加權(quán)求和進(jìn)行測算，權(quán)重分別為0.270 5、0.485 6。溫濕指數(shù)多用于評價較高氣溫時濕度對人體舒適度的影響度，風(fēng)寒指數(shù)多用于評價較低氣溫時氣流對人體舒適度的影響度。溫濕指數(shù)（THI）和風(fēng)寒指數(shù)（K）的計算公式為：

式（8）至（10）中，THI為溫濕指數(shù)和K為風(fēng)寒指數(shù)；T為華氏溫度（°F）；t為攝氏溫度（°C），f為空氣相對濕度（%）；v為風(fēng)速（m/s）；s為日照時數(shù)（h/d）?？諝赓|(zhì)量（X4）使用空氣質(zhì)量優(yōu)良率進(jìn)行刻畫，為京郊各區(qū)全年空氣質(zhì)量等級為優(yōu)良以上的天數(shù)占比。

（2）對社會政治因素的定量化處理。交通通達(dá)度（X11）使用京郊各區(qū)年路網(wǎng)密度進(jìn)行表征，為京郊各區(qū)年道路通車?yán)锍膛c京郊各區(qū)的區(qū)域面積之比。相關(guān)政策（X14）選取京郊各區(qū)年政府工作報告中提及“鄉(xiāng)村旅游”一詞的頻數(shù)。

數(shù)據(jù)的區(qū)屬處理。海拔高度（X1）、河流水系（X2）按照具體數(shù)據(jù)的區(qū)間類型進(jìn)行分級，其中，海拔高度（X1）的區(qū)間設(shè)置參考不同地形的海拔高度進(jìn)行劃分；河流水系（X2）緩沖區(qū)半徑設(shè)置依據(jù)漁網(wǎng)格點距離河流的距離，由于62.45%的格點位于2km范圍內(nèi)的緩沖區(qū)，構(gòu)建0～2km的4級河流緩沖區(qū)，據(jù)此進(jìn)行等級類型劃分。影響因素X3～X15均為區(qū)域統(tǒng)計數(shù)據(jù)，按照行政區(qū)劃進(jìn)行歸屬，借助ArcGIS 10.2軟件使用自然間斷點分級法對每年的具體數(shù)值進(jìn)行空間離散化處理，分為5個等級，見表6，將數(shù)值量轉(zhuǎn)化為類型量。使用ArcGIS 10.2軟件對京郊10區(qū)進(jìn)行格點化處理，設(shè)置漁網(wǎng)為2km×2km，生成3 750個格點，而后匹配每個格點所在區(qū)域的鄉(xiāng)村旅游地核密度值和15個影響因素的等級類型，分別作為地理探測器的因變量Y和自變量X數(shù)據(jù)進(jìn)行提取。

表6 2010—2020年京郊鄉(xiāng)村旅游地空間分布格局二級影響因素因子探測

3.2 影響因素的地理探測結(jié)果分析

借助地理探測器對2010—2020年京郊鄉(xiāng)村旅游地空間分布演化的二級影響因素進(jìn)行風(fēng)險探測、因子探測、生態(tài)探測和交互探測分析。

（1）風(fēng)險探測：可識別不同影響因素對京郊鄉(xiāng)村旅游地空間分布的影響是否具有差異性。探測結(jié)果表明在0.05的顯著性水平檢驗下，對京郊鄉(xiāng)村旅游地空間分布具有顯著影響的因素主要有交通通達(dá)度（X11）、海拔高度（X1）、城鎮(zhèn)居民可支配收入（X12）、常住城鎮(zhèn)人口數(shù)（X13）、地區(qū)生產(chǎn)總值（X6）和A級旅游景區(qū)數(shù)量（X9）。

（2）生態(tài)探測：用于分析比較兩個影響因素對京郊鄉(xiāng)村旅游地空間分布的影響作用能力是否具有顯著差異。探測結(jié)果表明海拔高度（X1）、空氣質(zhì)量（X4）、A級旅游景區(qū)數(shù)量（X9）、常住城鎮(zhèn)人口數(shù)（X13）的影響作用能力要強(qiáng)于其他影響因素。其他影響因素之間的影響作用能力差異在統(tǒng)計上不顯著。

（3）因子探測：因子探測的q值大小可表示該影響因素對京郊鄉(xiāng)村旅游地空間分布的影響作用能力強(qiáng)弱。根據(jù)表6的因子探測結(jié)果可以看出，2010—2020年海拔高度對京郊鄉(xiāng)村旅游地空間分布的影響作用能力最強(qiáng)，從2015年開始，海拔高度的影響作用能力均保持在首位。北京郊區(qū)地形多以山地為主，山區(qū)占北京市總面積的62%，京郊鄉(xiāng)村旅游地大多分布于海拔較低、坡度緩和的自然溝域地帶。京郊實施溝域經(jīng)濟(jì)開發(fā)戰(zhàn)略帶動了溝域鄉(xiāng)村交通基礎(chǔ)設(shè)施建設(shè)和鄉(xiāng)村旅游資源開發(fā)，易于形成鄉(xiāng)村旅游地。

通過對2010—2020年各影響因素的q值求和，發(fā)現(xiàn)旅游收入、A級旅游景區(qū)數(shù)量、地區(qū)生產(chǎn)總值和交通通達(dá)度也具有較強(qiáng)的影響作用能力。旅游收入反映了區(qū)域整體旅游產(chǎn)業(yè)的發(fā)展水平，鄉(xiāng)村旅游作為區(qū)域旅游產(chǎn)業(yè)的重要組成部分，旅游發(fā)展水平較高的地區(qū)，具有完善的旅游服務(wù)體系和較強(qiáng)的旅游吸引力，能夠為鄉(xiāng)村旅游的發(fā)展提供較好的支撐條件。地區(qū)生產(chǎn)總值高的區(qū)域經(jīng)濟(jì)發(fā)展水平高，公共基礎(chǔ)設(shè)施完備，居民收入和消費水平高，投資環(huán)境好，能夠促進(jìn)鄉(xiāng)村旅游地的建設(shè)與發(fā)展。A級旅游景區(qū)往往具有較好的品牌效應(yīng)和市場規(guī)模效應(yīng)，旅游基礎(chǔ)服務(wù)設(shè)置較為完善，可以對景區(qū)周邊的鄉(xiāng)村產(chǎn)生較強(qiáng)的輻射帶動作用，有利于周邊的鄉(xiāng)村迅速開展旅游活動。交通通達(dá)度決定著游客前往鄉(xiāng)村旅游地的便利程度，尤其是在鄉(xiāng)村旅游地形成發(fā)展初期，會直接影響游客對鄉(xiāng)村旅游地的選擇，因此，2012年的交通通達(dá)度的q值為第一位。

（4）交互探測：用于識別兩個影響因素共同作用于京郊鄉(xiāng)村旅游地空間分布時，其影響作用能力會否增強(qiáng)或者削弱。結(jié)果表明，2010—2020年15個影響因素中任意兩個影響因素交互作用后其影響作用能力均呈非線性加強(qiáng)和雙因子加強(qiáng)狀態(tài)，其中非線性加強(qiáng)狀態(tài)最多。如2020年海拔高度與其他影響因素交互作用后呈非線性加強(qiáng)狀態(tài)，交互產(chǎn)生的影響作用能力均大于兩個因素影響作用能力之和。這也進(jìn)一步說明京郊鄉(xiāng)村旅游地空間分布差異是多影響因素綜合作用的結(jié)果。

3.3 影響因素的影響作用能力演化分析

參考文獻(xiàn)[51]的計算方法，將各一級影響因素下所對應(yīng)的5個二級影響因素的因子探測q值進(jìn)行相加，將其作為一級影響因素的q值（圖5），分析各級影響因素的q值演化特征。

圖5 2010—2020年一級影響因素q值探測結(jié)果及變化

自然環(huán)境因素的影響作用能力波動上升，總體保持第三位。自然環(huán)境因素的因子探測q值及變化趨勢如圖6所示：自然環(huán)境因素的影響作用能力在2011年略下降后，2012年快速上升，2013—2020年總體保持著波動上升的發(fā)展趨勢。京郊由于山區(qū)廣布，鄉(xiāng)村旅游地的形成與后續(xù)發(fā)展在考慮自然生態(tài)環(huán)境、河流水系分布的同時，更注重對地形地貌的考量。海拔低、地勢相對平坦區(qū)域，利于對外交通路線的建設(shè)和農(nóng)用地的連片集中使用，更有利于鄉(xiāng)村旅游發(fā)展。因此，2010—2020年自然環(huán)境因素的影響作用能力總體呈上升趨勢，其中，海拔高度、河流水系等二級影響因素的影響作用能力逐年上升，2015年后海拔高度的q值居所有二級影響因子的第一位。

圖6 2010—2020年自然環(huán)境因素q值探測結(jié)果及變化

經(jīng)濟(jì)產(chǎn)業(yè)因素的影響作用能力雖呈下降，但總體居第一位。經(jīng)濟(jì)產(chǎn)業(yè)因素的因子探測q值及變化趨勢如圖7所示，可以看出：2010—2013年經(jīng)濟(jì)產(chǎn)業(yè)因素的影響作用能力呈穩(wěn)步上升趨勢，2014年開始下降后直至2019年才有明顯回升。從整體來看，經(jīng)濟(jì)產(chǎn)業(yè)因素影響作用能力最高的一級影響因素。旅游收入反映了區(qū)域旅游業(yè)發(fā)展的整體水平，2010—2020年旅游收入的q值在2016年和2017年處于15個二級影響因子影響作用力的第2位、第3位。地區(qū)生產(chǎn)總值和第三產(chǎn)業(yè)占比為區(qū)域旅游業(yè)的發(fā)展提供基礎(chǔ)，地區(qū)生產(chǎn)總值越高的地區(qū)經(jīng)濟(jì)水平越大發(fā)，第三產(chǎn)業(yè)占比越高區(qū)域旅游服務(wù)水平越高，有利于鄉(xiāng)村旅游地的發(fā)展和管理，所以這兩個二級影響因素的影響作用力總體分別處于第4位、第6位。農(nóng)林牧漁業(yè)總產(chǎn)值是鄉(xiāng)村發(fā)展的基礎(chǔ)，農(nóng)林牧漁業(yè)總產(chǎn)值越高，鄉(xiāng)村經(jīng)濟(jì)發(fā)展水平越高，有利于鄉(xiāng)村旅游的開展和鄉(xiāng)村旅游地的形成，其影響作用力總體分別處于第7位。

圖7 2010—2020年經(jīng)濟(jì)產(chǎn)業(yè)因素q值探測結(jié)果及變化

社會政治因素的影響作用能力呈波動下降，總體居第二位。社會政治因素的因子探測q值及變化趨勢，由圖8可以看出：社會政治因素的影響作用能力2010—2012年為上升趨勢，2013—2014年有所下降，2015—2016年短暫回升后，2017—2020年保持緩慢下降趨勢。旅游交通決定了到達(dá)京郊旅游地的便利程度，隨著交通路網(wǎng)的建設(shè)與完善，區(qū)域交通通達(dá)上升，其對鄉(xiāng)村旅游地空間分布的影響作用能力呈下降趨勢，q值由2012年的0.059 3下降為2020年的0.014 6。旅游政策和鄉(xiāng)村建設(shè)投入引導(dǎo)鄉(xiāng)村旅游地的發(fā)展方向，但隨著鄉(xiāng)村旅游的繁榮，鄉(xiāng)村旅游地的建設(shè)由注重數(shù)量的增減，逐漸轉(zhuǎn)向為注重數(shù)量的增加和兼顧質(zhì)量的提升，其q數(shù)總體呈波動上升的趨勢。京郊鄉(xiāng)村旅游主要以吸引城鎮(zhèn)游客群體到鄉(xiāng)村游玩，但隨著北京鄉(xiāng)村旅游市場的繁榮，鄉(xiāng)村旅游的游客群體主要以北京市城6區(qū)的游客為主，京郊地區(qū)的城鎮(zhèn)居民對鄉(xiāng)村旅游地分布的影響逐漸減弱。

圖8 2010—2020年社會政治因素q值探測結(jié)果及變化

4 結(jié)論與建議

4.1 結(jié)論

該文以京郊鄉(xiāng)村旅游特色業(yè)態(tài)和鄉(xiāng)村民俗旅游村為研究對象，借助ArcGIS 10.2軟件，運用重心標(biāo)準(zhǔn)差橢圓、地理集中指數(shù)、基尼系數(shù)、核密度分析等空間分析和計量地理學(xué)模型，對2010—2020年京郊鄉(xiāng)村旅游地空間分布演化進(jìn)行了分析。并在此基礎(chǔ)上，借助地理探測器的定量分析，探討了影響京郊鄉(xiāng)村旅游地空間分布演化的影響因素，得出以下主要結(jié)論。

（1）2010—2020年京郊鄉(xiāng)村旅游地發(fā)展的重心變化較為穩(wěn)定，均位于昌平區(qū)范圍內(nèi)，重心遷移軌跡總體呈現(xiàn)“先東南方向，后東北方向”的遷移變化特征。2010—2020年京郊鄉(xiāng)村旅游地的標(biāo)準(zhǔn)差橢圓的形態(tài)變化幅度不大，橢圓方位基本呈現(xiàn)“東北—西南”方向分布，長軸呈“減小—增大”的過程，短軸變化平穩(wěn)。

（2）2010—2020年京郊鄉(xiāng)村旅游地在京郊各區(qū)內(nèi)部的空間分布呈現(xiàn)較為顯著的集聚分布、不均衡分布特征，且呈逐年上升的趨勢。

（3）2010—2020年京郊鄉(xiāng)村旅游地的密度分布高密度區(qū)數(shù)量逐年增多并呈連片拓展，逐漸形成以昌平、懷柔為雙高密度核心區(qū)，密云、延慶、大興多個高密度連片區(qū)并存的空間格局。

（4）京郊鄉(xiāng)村旅游地空間分布演化受自然環(huán)境因素、經(jīng)濟(jì)產(chǎn)業(yè)因素和社會政策因素的共同影響。自然環(huán)境因素次一級的二級影響因素包括：海拔高度、河流水系、氣候舒適度、空氣質(zhì)量和林木綠化率。經(jīng)濟(jì)產(chǎn)業(yè)因素次一級的二級影響因素包括；地區(qū)生產(chǎn)總值、第三產(chǎn)業(yè)占比、旅游收入、A級旅游景區(qū)數(shù)量和農(nóng)林牧漁業(yè)總產(chǎn)值。社會政治因素次一級的二級影響因素包括：交通通達(dá)度、城鎮(zhèn)居民可支配收入、常住城鎮(zhèn)人口數(shù)、相關(guān)政策和農(nóng)林水事務(wù)支出。

（5）海拔高度、旅游收入、A級旅游景區(qū)數(shù)量、地區(qū)生產(chǎn)總值和交通通達(dá)度是影響京郊鄉(xiāng)村旅游地空間分布的主要影響因素；任意兩個影響因素交互作用后對京郊鄉(xiāng)村旅游地空間分布的影響作用能力均強(qiáng)于單個因素的作用?？傮w來看，自然環(huán)境因素的影響作用能力逐漸上升，居第三位；社會政治因素的影響能力呈波動下降，居第二位；經(jīng)濟(jì)產(chǎn)業(yè)因素的影響能力逐漸下降，仍居第一位。

4.2 建議

以往鄉(xiāng)村旅游地空間分布影響因素的研究，大多以單一定性分析為主，籠統(tǒng)地將經(jīng)濟(jì)發(fā)展水平、資源稟賦、地方政策等影響因素并置，缺少主次關(guān)系的判別，且研究時間多為單一靜態(tài)年份。該文借助計量地理學(xué)模型和GIS空間分析方法，嘗試對2010—2020年京郊鄉(xiāng)村旅游地的空間分布演化及其影響因素進(jìn)行分析梳理，根據(jù)影響因素作用能力的演化，綜合判別影響因素的前后關(guān)系，為京郊鄉(xiāng)村旅游地空間布局的優(yōu)化提供決策依據(jù)。同時，將氣候舒適度、空氣質(zhì)量納入到影響因素中進(jìn)行定量化分析，豐富了原有討論要素的框架范疇。但該文在影響因素的選取方面，暫未選取北京6個主城區(qū)的外力作用影響，后續(xù)研究可對其進(jìn)行深入分析。此外，該文僅對于京郊鄉(xiāng)村旅游地的空間分布演化現(xiàn)狀進(jìn)行了分析，對于京郊鄉(xiāng)村旅游地未來可能的空間分布發(fā)展趨勢未做出預(yù)測，后續(xù)研究可持續(xù)關(guān)注。對京郊鄉(xiāng)村旅游地的發(fā)展與空間布局優(yōu)化，提出以下幾點建議。

（1）優(yōu)化空間布局，加強(qiáng)區(qū)域聚集發(fā)展。依據(jù)產(chǎn)業(yè)集聚和區(qū)域平衡的發(fā)展布局原則，在京郊鄉(xiāng)村旅游地分布稀疏、集中程度不高的地區(qū)，如順義、通州，扶持鄉(xiāng)村農(nóng)業(yè)發(fā)展與村莊基礎(chǔ)設(shè)施建設(shè)，加強(qiáng)鄉(xiāng)村旅游資源的挖掘，逐步提升鄉(xiāng)村旅游發(fā)展實力。在高密度核心區(qū)的昌平和懷柔，高密度連片區(qū)域的密云、延慶、大興，強(qiáng)化其輻射擴(kuò)散效應(yīng)的發(fā)揮，整合周邊鄉(xiāng)村旅游資源優(yōu)勢，優(yōu)化各類資源要素的配置，帶動周邊鄉(xiāng)村的發(fā)展，以形成集聚效應(yīng)，實現(xiàn)鄉(xiāng)村旅游地的“雁陣”式發(fā)展格局。

（2）構(gòu)建溝域綠道系統(tǒng)，增強(qiáng)鄉(xiāng)村旅游地發(fā)展活力。北京郊區(qū)部分區(qū)域，如延慶、密云，鄉(xiāng)村旅游地聚集分于低海拔的溝域地帶，可根據(jù)溝域地勢特點構(gòu)建綠道系統(tǒng)，將區(qū)域范圍內(nèi)的各個村莊景點進(jìn)行串聯(lián)，實現(xiàn)村落景觀化、農(nóng)戶景觀化，助推鄉(xiāng)村旅游地的健康發(fā)展。同時，借助綠道建設(shè)，增強(qiáng)更多村莊、景點的宣傳機(jī)會和可進(jìn)入性，推動更多鄉(xiāng)村旅游地的形成與發(fā)展。

（3）以景區(qū)為依托，帶動景區(qū)周邊鄉(xiāng)村旅游地的“景村融合”發(fā)展。北京郊區(qū)擁有眾多高品質(zhì)的國家A級旅游景區(qū)，如八達(dá)嶺長城風(fēng)景名勝區(qū)、京東石林峽風(fēng)景區(qū)等，可就近將鄉(xiāng)村旅游地和國家A級旅游景區(qū)結(jié)合在一起，構(gòu)建“景村融合”協(xié)調(diào)發(fā)展機(jī)制。景區(qū)帶動周邊鄉(xiāng)村圍繞景區(qū)拓展業(yè)態(tài)，實現(xiàn)景區(qū)與毗鄰鄉(xiāng)村旅游地的優(yōu)勢互補(bǔ)、共同發(fā)展，同時，景區(qū)與周邊村莊聯(lián)動打包宣傳，共享旅游服務(wù)配套設(shè)施，共享客源市場，逐步形成景區(qū)村莊一體化融合的發(fā)展格局。