杜久升，馮云超，王 羽，趙貝貝

(河南理工大學 測繪與國土信息工程學院，河南 焦作 454003)

鄉(xiāng)村旅游是以鄉(xiāng)土民情、農耕文明、鄉(xiāng)野特色和田園風光等農村資源為依托，以村莊、山野為空間形成的集休閑觀光、度假學習、農居體驗、農事操作、傳統(tǒng)手工藝制作等活動于一體的新興旅游形式.2015—2023年中央一號文件連續(xù)提及鄉(xiāng)村旅游，為鄉(xiāng)村旅游發(fā)展提供了方向和動力.近年來，有關部門先后制定鄉(xiāng)村旅游名錄，標志著鄉(xiāng)村旅游建設逐漸步入體系化、制度化的新階段.河南省鄉(xiāng)村旅游發(fā)展于20世紀90年代中期[1]，包括洛陽市欒川縣重渡溝村、新鄉(xiāng)市輝縣市郭亮村等均屬于較早一批的鄉(xiāng)村旅游地.經過近30年的建設與發(fā)展，河南省鄉(xiāng)村旅游地已經覆蓋全省18個地市，其形式也由最初的農家餐飲服務逐步發(fā)展為以體驗、康養(yǎng)、休閑為主的旅游服務.2020年5月，河南省人民政府辦公廳印發(fā)《關于加快鄉(xiāng)村旅游發(fā)展的意見》，提出到2025年，打造1 000個鄉(xiāng)村旅游特色村，實現(xiàn)全省鄉(xiāng)村旅游年接待游客4億人次，年經營總收入達到3 000億元的目標.然而，不同村落的資源稟賦和區(qū)位條件迥異，導致鄉(xiāng)村旅游特色村的空間分布形態(tài)復雜多樣，其空間異質性為鄉(xiāng)村旅游的建設和發(fā)展帶來一定的困擾.因此，立足實際，探究鄉(xiāng)村旅游特色村的空間布局和結構類型，厘清影響鄉(xiāng)村旅游特色村空間分布的影響因素十分必要.

國外關于鄉(xiāng)村旅游的研究起步較早，目前已建立了較為完備的研究體系，主要圍繞鄉(xiāng)村旅游的定義[2]、空間格局[3]和驅動機制[4]等方面進行探討.國內關于鄉(xiāng)村旅游的研究起源于20世紀90年代，以發(fā)展模式[5]、發(fā)展問題與對策[6]、發(fā)展驅動機制[7]等方面為主.近年來研究重點逐漸轉向鄉(xiāng)村旅游空間分布及影響因素的研究.研究方法涉及平均最近鄰指數[8]、核密度估計法[9]、網絡維數分析法[10]、地理聯(lián)系率[11]、相關性分析[12]等方法.總的來看，學者從不同視角對鄉(xiāng)村旅游進行了廣泛探討，并取得了豐富的理論研究成果.然而已有研究對鄉(xiāng)村旅游空間分布特點的研究較少，對鄉(xiāng)村旅游空間分布影響因素的量化分析、影響因素之間的相互關系以及不同因素的疊加作用分析較為欠缺.

基于此，以河南省723個鄉(xiāng)村旅游特色村為研究對象，運用空間分析方法和地理探測器模型，分析鄉(xiāng)村旅游特色村空間布局及影響因素，以期為河南省鄉(xiāng)村旅游的開發(fā)與建設以及高質量發(fā)展提供參考.

1 數據來源與研究方法

1.1 數據來源

河南省鄉(xiāng)村旅游特色村是指旅游特色鮮明，文化旅游資源富集，生態(tài)環(huán)境和諧，組織管理有序，具備一定基礎設施和接待服務能力，形成一定旅游規(guī)模、效益和市場吸引力的村落[13].文中所使用的鄉(xiāng)村旅游特色村數據來源于2018—2022年河南省文化和旅游廳公布的《河南省鄉(xiāng)村旅游特色村名單》，整合上述數據，得到河南省鄉(xiāng)村旅游特色村共723個.村落坐標及其他點坐標通過“百度地圖API坐標拾取器”提?。缓幽鲜⌒姓^(qū)劃數據、路網數據等矢量數據來源于國家基礎地理信息數據庫和OpenStreetMap數據集；數字高程模型數據來源于地理空間數據云；所涉及的統(tǒng)計數據來源于《中國縣域統(tǒng)計年鑒》及《河南省統(tǒng)計年鑒》.

1.2 研究方法

1.2.1 平均最近鄰指數 最近鄰距離能夠衡量鄉(xiāng)村旅游特色村在空間上的鄰近程度，平均最近鄰指數R能夠反映鄉(xiāng)村旅游特色村的空間集散特征，即“平均觀測最近鄰距離”與“理論平均最近鄰距離”的比值.計算公式為

(1)

1.2.2 核密度分析 核密度分析用于識別鄉(xiāng)村旅游特色村的空間集聚區(qū)域.核密度估計法認為，事件的發(fā)生可以在地理空間的任意位置，但是在不同位置發(fā)生的概率是不同的[15].核密度值的高低反映的是研究對象在空間分布上的集聚程度.計算公式為

(2)

其中，f(x)為鄉(xiāng)村旅游特色村的核密度值；n為搜索半徑內鄉(xiāng)村旅游特色村的數量；h為搜索半徑；k為權重值；(x-xi)為鄉(xiāng)村旅游特色村x到樣本鄉(xiāng)村旅游特色村xi的距離[16].

1.2.3 地理探測器 地理探測器是探測空間分異性，揭示其背后驅動力的統(tǒng)計學方法.其核心思想是，如果某個自變量對因變量有重要影響，那么自變量和因變量的空間分布應該具有相似性.目前已被廣泛運用于旅游[17]和環(huán)境[18]等多個領域.地理探測器包括4個部分：因子探測、交互作用探測、風險區(qū)探測和生態(tài)探測.其中因子探測和交互作用探測使用q值來度量因子X對變量Y的解釋力，值域為[0,1]，q值越大，說明因子X對變量Y的解釋力越強，反之則越弱.計算公式為

(3)

2 空間分布

2.1 空間分布類型

通過平均最近鄰分析，得到河南省鄉(xiāng)村旅游特色村理論最近鄰距離為6 402 m，實際最近鄰距離為9 795 m，平均最近鄰指數為0.654，通過了1%的顯著性檢驗，說明河南省鄉(xiāng)村旅游特色村具有集聚分布的特征.這種空間集聚特征有利于村與村之間通過內部合作，共同利用外部資源獲得集聚效益[20].

2.2 空間分布格局

為揭示鄉(xiāng)村旅游特色村的空間集聚區(qū)域，利用核密度分析工具，可視化河南省鄉(xiāng)村旅游特色村空間分布格局，如圖1所示.河南省鄉(xiāng)村旅游特色村呈現(xiàn)以“多核心-鏈狀”集聚分布格局.高密度區(qū)主要集中在豫北地區(qū)和豫中地區(qū).其中豫北地區(qū)和豫南地區(qū)分別呈現(xiàn)以安陽-鶴壁-新鄉(xiāng)、焦作、濟源、信陽為核心的“多核心”集聚分布；豫中地區(qū)和豫西地區(qū)呈現(xiàn)以鄭州-許昌-平頂山、三門峽-洛陽為“鏈狀”的集聚分布.總體來看，河南省鄉(xiāng)村旅游特色村空間分布不均衡，呈現(xiàn)“北密南疏，西密東疏”的空間分布格局，各地市差異顯著.

圖1 鄉(xiāng)村旅游特色村核密度

3 影響因素分析

3.1 指標構建與處理

文化生態(tài)學理論認為，文化事物的產生、發(fā)展受到自然、社會等環(huán)境的影響[21].通過對河南省鄉(xiāng)村旅游特色村推薦條件的解讀，發(fā)現(xiàn)河南省鄉(xiāng)村旅游特色村的推薦條件包括產品特色突出、基礎設施完善、服務功能齊全、生態(tài)環(huán)境良好、旅游效益突出和組織管理有序.借鑒如表1所示的已有研究，參考河南省鄉(xiāng)村旅游特色村的推薦條件，兼顧數據的可獲取性，運用數理統(tǒng)計、疊加分析、緩沖區(qū)分析等研究方法模擬探測可能的驅動因子.

表1 鄉(xiāng)村旅游空間分布影響因素文獻總結

1)海拔高度

海拔高度對旅游資源的空間分布具有較強的制約作用，一定程度上影響著鄉(xiāng)村旅游資源的形成和分布.提取鄉(xiāng)村旅游特色村的高程信息，得到河南省鄉(xiāng)村旅游特色村分布在海拔31～1 314 m范圍內.具體來看，河南省鄉(xiāng)村旅游特色村大多分布于地形起伏多變的山區(qū)，包括豫北地區(qū)太行山，豫西地區(qū)小秦嶺、王屋山、崤山、伏牛山，豫南地區(qū)桐柏山、大別山和豫中地區(qū)的嵩山，如圖2所示.以200 m為間隔將海拔高度劃分為6個等級，統(tǒng)計不同海拔范圍內鄉(xiāng)村旅游特色村的數量，結果如圖3所示.可以看出，隨著海拔高度的上升，鄉(xiāng)村旅游特色村數量逐漸減少，其中海拔高度在32～200 m之間的村落數量最多，占總數的56.43%；海拔高度高于1 000 m的村莊數量最少，僅占總數的2.22%，說明河南省鄉(xiāng)村旅游特色村具有“依山”的空間分布特點，也表明海拔高度一定程度上影響著鄉(xiāng)村旅游特色村的分布狀況.

圖2 鄉(xiāng)村旅游特色村與山脈、河流

圖3 鄉(xiāng)村旅游特色村海拔高度分級統(tǒng)計圖

2)河流水系

河流水系不僅能夠為周圍村莊的農業(yè)生產和生態(tài)建設提供水源，還可以開發(fā)劃船、垂釣等鄉(xiāng)村旅游項目.以河南省等級河流為基礎，在河流兩側分別建立以2 km為間隔的多級緩沖區(qū)，再與鄉(xiāng)村旅游特色村進行疊加分析，統(tǒng)計緩沖區(qū)范圍內的村落個數，如圖4所示.可以看出，2 km范圍內鄉(xiāng)村旅游特色村的數量為312個，占總數的43.15%；4 km范圍內鄉(xiāng)村旅游特色村數量為484個，占總數的66.94%.距離河流越遠，村落數量越少，表明河南省鄉(xiāng)村旅游特色村具有“傍水”的空間分布特點，河流水系與鄉(xiāng)村旅游特色村的分布密切相關.

圖4 鄉(xiāng)村旅游特色村沿河流分布統(tǒng)計

3)交通干線

良好的交通條件能夠促進城鄉(xiāng)要素的流通，大大壓縮游客的出行時間.以2022年河南省主要公路為基礎，在公路兩側分別建立以2 km為間隔的多級緩沖區(qū)，并與鄉(xiāng)村旅游特色村進行疊加分析，統(tǒng)計緩沖區(qū)范圍內的村落個數，如圖5所示.可以看出，2 km范圍內鄉(xiāng)村旅游特色村的數量為575個，占總數的79.52%，4 km范圍內鄉(xiāng)村旅游特色村的數量為680個，占總數的94.05%，表明河南省鄉(xiāng)村旅游特色村具有“沿路”的空間分布特點，鄉(xiāng)村旅游特色村沿主要公路分布的特征明顯.

圖5 鄉(xiāng)村旅游特色村沿公路分布統(tǒng)計

4)區(qū)位條件

鄉(xiāng)村旅游特色村的區(qū)位條件決定著游客出行的時間成本和消費成本， 進而影響鄉(xiāng)村旅游的發(fā)展.

以區(qū)(縣)政府所在地代表城區(qū)，建立多級緩沖區(qū)，計算不同距離范圍內鄉(xiāng)村旅游特色村的數量，如圖6所示.可以看出，隨著與城區(qū)之間距離的增加，鄉(xiāng)村旅游特色村的數量呈現(xiàn)先增加后減少的變化特征.在距離城區(qū)30 km的范圍內涵蓋了665個鄉(xiāng)村旅游特色村，占總數的91.97%，表明鄉(xiāng)村旅游特色村主要集聚在城區(qū)周邊.其中10～15 km區(qū)間范圍內鄉(xiāng)村旅游特色村數量最多，占總數的25.59%，該區(qū)域多為城區(qū)與鄉(xiāng)村的過渡地帶，擁有豐富的自然景觀和便利的出行條件，能夠滿足城市游客逃離城市快節(jié)奏生活，親近大自然的旅游的需求.同時表明河南省鄉(xiāng)村旅游特色村具有“環(huán)城”的空間分布特點.

圖6 鄉(xiāng)村旅游特色村距城區(qū)距離分級統(tǒng)計

5)A級景區(qū)

A級景區(qū)擁有完善的基礎設施、豐富的旅游資源和良好的品牌效應，能夠形成旅游地增長級，通過乘數效應推動其他鄉(xiāng)村旅游地的發(fā)展.運用緩沖區(qū)分析，計算不同距離緩沖區(qū)內鄉(xiāng)村旅游特色村的數量， 結果如圖7所示， 隨著距離A級景區(qū)的距離不斷增加，鄉(xiāng)村旅游特色村的數量呈現(xiàn)減少趨勢.其中0～5 km范圍內特色村數量最多，達到276個，占總數的38.17%.大于30 km的特色村數量最少，僅為4個，僅占總數的0.89%，表明河南省鄉(xiāng)村旅游特色村具有“近景”的空間分布特點.A級景區(qū)與鄉(xiāng)村旅游特色村的分布密切相關，周圍的村莊可以依托A級景區(qū)實現(xiàn)聯(lián)動發(fā)展.

圖7 鄉(xiāng)村旅游特色村距A級景區(qū)距離分級統(tǒng)計

6)經濟發(fā)展水平

經濟基礎是進行旅游活動的前提條件.鄉(xiāng)村旅游作為特定區(qū)域的旅游活動，其空間分布受到經濟規(guī)模的影響.選取地區(qū)生產總值、第三產業(yè)占比、旅游收入等指標計算其與鄉(xiāng)村旅游特色村數量的地理聯(lián)系率.結果顯示，地區(qū)生產總值、第三產業(yè)占比、旅游收入與鄉(xiāng)村旅游特色村數量的地理聯(lián)系率分別為99.73,99.84,99.62，非常接近100，表明經濟發(fā)展水平與鄉(xiāng)村旅游特色村具有高度的正向關聯(lián).

通過對可能的驅動因子進行模擬探測，結合已有研究，最終確定自然環(huán)境、區(qū)位交通、社會經濟3個方面共12項指標.自然環(huán)境方面包括海拔高度、距河流水系距離、年均AQI、林地面積、耕地面積；區(qū)位交通方面包括距縣級行政中心距離、距主要交通線距離、距高級別景區(qū)距離；社會經濟方面包括地區(qū)生產總值、旅游收入、第三產業(yè)占比、一般財政預算支出.

地理探測器要求如果自變量X為數值量，則需要對數據進行離散化處理.傳統(tǒng)地理探測器方法通常采用先驗知識，但這種方法的主觀性會對結果產生干擾.文中運用參數最優(yōu)地理探測器(Optimal parameters-based geographical detector, OPGD)，借助R語言中的GD包，選取q值最高的離散化方法和區(qū)間數量進行離散化處理，離散化的方法包括相等間隔分類法(Equal breaks)、自然間斷點分類法(Natural breaks)、分位數分類法(Quantile breaks)、幾何間隔分類法(Geometric breaks).使用ArcGIS構建3×3 km的漁網作為研究單元，共獲取19 073個采樣點.核密度值最適合揭示鄉(xiāng)村旅游特色村的空間分異，因此將核密度值和各自變量與漁網進行空間賦值，以平均值代表各漁網的因子屬性.

3.2 單因子影響力分析

單因子探測結果顯示如表2所示，所有因素P值均為0，表明各因素對鄉(xiāng)村旅游特色村的空間分布均有顯著的影響.各因素對鄉(xiāng)村旅游特色村空間分布的影響力從強到弱依次為：一般財政預算支出(X12)>旅游收入(X10)>地區(qū)生產總值(X9)>第三產業(yè)占比(X11)>年均AQI(X3)>距高級別景區(qū)距離(X8)>林地面積(X4)>海拔高度(X1)>耕地面積(X5)>距縣級行政中心距離(X6)>距河流水系距離(X2)>距主要交通線距離(X7).總的來看，社會經濟對鄉(xiāng)村旅游特色村空間分布影響最強，自然環(huán)境次之，區(qū)位交通較弱.地方經濟狀況奠定了鄉(xiāng)村旅游發(fā)展的基礎，能夠為鄉(xiāng)村旅游的發(fā)展提供較好的支撐條件，在鄉(xiāng)村旅游特色村的組織形態(tài)和演化格局中發(fā)揮著主導作用.具體來看，自然環(huán)境中影響力最強的是空氣質量，空氣質量是區(qū)域生態(tài)環(huán)境好壞的最直接表征，良好的生態(tài)環(huán)境是河南省鄉(xiāng)村旅游特色村推薦條件之一；區(qū)位交通中影響力最強的是距高級別景區(qū)距離，A級景區(qū)的開發(fā)與完善能夠為周邊的村莊帶來大量的客源，周邊鄉(xiāng)村可依托A級景區(qū)帶來的市場優(yōu)勢開發(fā)多樣的鄉(xiāng)村旅游產品，從而帶動鄉(xiāng)村旅游的發(fā)展；社會經濟中影響力最強的是一般財政預算支出，政府支持力度的加強能夠對鄉(xiāng)村旅游的空間布局起到積極的導向作用和調控作用.

表2 影響因素及因子探測結果

3.3 雙因子影響力分析

各因素對鄉(xiāng)村旅游特色村空間分布的影響力大小線性等間隔反映如圖8所示.可以看出，不同影響因子交互作用的影響力均大于因子單獨作用時的影響力，不存在相互獨立的因子.各因子間具有較強的雙因子增強作用，即兩個影響因子在交互作用時產生的影響力大于各自單獨作用時的影響力之和，甚至部分因子間表現(xiàn)出非線性增強作用，即兩個影響因子在交互作用時產生的影響力大于它們單獨作用時的影響力.說明12個影響因子中任意兩個因子的交互作用對鄉(xiāng)村旅游特色村的空間格局的影響力更為顯著，也就是說鄉(xiāng)村旅游特色村的空間分布受到不同維度不同因子的共同影響.

具體來看，各因子兩兩交互呈現(xiàn)出解釋力最強的模式為“海拔高度(X1)∩旅游收入(X10)”，q值為0.593，顯著高于海拔高度(X1)的q值(0.122)與旅游收入(X10)的q值(0.342)之和.在q值大于0.5的交互作用中與海拔高度(X1)的交互最為頻繁.在鄉(xiāng)村旅游特色村建設初期，鄉(xiāng)村旅游的發(fā)展更多依賴于自然資源稟賦，“地”在一定程度上決定了村莊的規(guī)模，而海拔較低，坡度較緩的山區(qū)地勢起伏多變，景觀多樣，能夠為鄉(xiāng)村旅游特色村提供良好的條件，體現(xiàn)了自然環(huán)境中的海拔高度對于鄉(xiāng)村旅游特色村空間分布的重要作用.

4 結論與討論

4.1 結論

以河南省723個鄉(xiāng)村旅游特色村為研究對象，運用平均最近鄰指數、核密度分析、地理探測器等方法，對鄉(xiāng)村旅游特色村的空間分布及影響因素進行研究，得出以下結論：

1)河南省鄉(xiāng)村旅游特色村平均最近鄰指數為0.654，表明鄉(xiāng)村旅游特色村的空間分布類型為集聚型.豫北地區(qū)和豫南地區(qū)呈現(xiàn)“多核心”集聚分布，豫中和豫南地區(qū)呈現(xiàn)“鏈狀”集聚分布.空間分布上具有“依山、傍水、沿路、環(huán)城、近景”的特點.

2)因子探測結果表明，鄉(xiāng)村旅游特色村的空間分布是多個因素作用的結果，社會經濟對鄉(xiāng)村旅游特色村空間分布影響最強，自然環(huán)境次之，區(qū)位交通較弱.空氣質量、距高級別景區(qū)距離和一般財政預算支出分別為自然環(huán)境、區(qū)位交通和社會經濟中影響力最強的.

3)交互作用探測結果表明，不同因子交互作用下呈現(xiàn)雙因子增強或非線性增強.交互作用最強的為“海拔高度(X1)∩旅游收入(X10)”.在q值大于0.5的交互作用中，與海拔高度(X1)的交互最為頻繁.

4.2 討論

發(fā)展鄉(xiāng)村旅游是構建現(xiàn)代鄉(xiāng)村產業(yè)體系的重要任務，是助力鄉(xiāng)村振興的重要抓手.針對河南省鄉(xiāng)村旅游的空間按布局，學者進行了一系列研究.本研究是在吳昕昱[28]和焦士興[29]等研究的基礎上，針對河南省鄉(xiāng)村旅游影響因素研究偏重定性分析，定量分析薄弱的特點以及不同離散化方法對地理探測器模型帶來的誤差影響，按照“推薦條件-文獻梳理-模擬探測-定量分析”的思路，運用參數最優(yōu)地理探測器在最優(yōu)分類的基礎上對影響因素影響力大小進行分析.從研究結果上看，河南省鄉(xiāng)村旅游特色村呈現(xiàn)集聚分布特征，自然環(huán)境、距高級別景區(qū)距離等因素對鄉(xiāng)村旅游地空間分布影響較大，與前人研究結果相一致.目前河南省鄉(xiāng)村旅游特色村空間分布不均衡，集聚分布特征明顯，豫東地區(qū)缺少核心城市帶動鄉(xiāng)村旅游發(fā)展，未來應當深入挖掘豫東地區(qū)鄉(xiāng)村旅游資源，整合區(qū)域自然資源、高級別景區(qū)資源等優(yōu)勢，發(fā)展鄉(xiāng)村旅游核心城市以彌補發(fā)展短板.

由于數據的可獲得性，在潛在影響因素的選擇上，僅對自然環(huán)境、區(qū)位交通、社會經濟3個維度的部分因素進行探討，未對氣候、政策等因素進行研究，今后的研究中還需對這些影響因素進行完善和更新，更加準確地揭示河南省鄉(xiāng)村旅游空間分布的影響因素.此外，鄉(xiāng)村旅游的發(fā)展是一個局部時空動態(tài)發(fā)展的過程，未來研究應當選擇具有時間序列的研究樣本，從時間和空間兩個視角，綜合定性、定量分析方法對鄉(xiāng)村旅游特色村的時空特征進行分析，更好地指導河南省鄉(xiāng)村旅游高質量發(fā)展.