白翠玲，雷欣，楊麗花2,，董志良

（1. 河北地質(zhì)大學(xué) 土地科學(xué)與空間規(guī)劃學(xué)院，河北 石家莊 050031；2. 河北省科協(xié)智庫(kù)地學(xué)旅游與鄉(xiāng)村振興研究基地，河北 石家莊 050031；3. 河北地質(zhì)大學(xué) 管理學(xué)院，河北 石家莊 050031；4. 河北地質(zhì)大學(xué) 科技處，河北 石家莊 050031）

隨著時(shí)代的發(fā)展，旅游逐漸成為促進(jìn)社會(huì)發(fā)展、改善生態(tài)環(huán)境的重要戰(zhàn)略性產(chǎn)業(yè)。然而，以環(huán)境為支撐、資源為條件的特性也決定了旅游開發(fā)與生態(tài)環(huán)境保護(hù)之間存在著一定的矛盾性[1]；我國(guó)旅游產(chǎn)業(yè)尚處于發(fā)展階段，過于注重產(chǎn)業(yè)經(jīng)濟(jì)效益可能會(huì)導(dǎo)致社會(huì)效益的損害[2]。因此，如何協(xié)調(diào)旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益之間的關(guān)系，提升區(qū)域發(fā)展質(zhì)量與效率，成為社會(huì)關(guān)注的熱點(diǎn)[3-4]。

大量研究表明旅游產(chǎn)業(yè)與生態(tài)環(huán)境、旅游產(chǎn)業(yè)與社會(huì)效益、社會(huì)效益與生態(tài)環(huán)境之間存在一定的互動(dòng)關(guān)系；旅游產(chǎn)業(yè)以生態(tài)環(huán)境為載體，對(duì)生態(tài)環(huán)境的影響具有二重性，而生態(tài)環(huán)境則為旅游業(yè)發(fā)展提供空間與客觀物質(zhì)條件，兩者互融共生[2,5-7]；旅游產(chǎn)業(yè)與社會(huì)效益之間具有長(zhǎng)期穩(wěn)定的均衡關(guān)系，兩者耦合互動(dòng)效應(yīng)明顯，旅游業(yè)既具有產(chǎn)業(yè)性質(zhì),又具有實(shí)現(xiàn)人民幸福生活目標(biāo)的事業(yè)性質(zhì)[4,8-9]；社會(huì)效益和生態(tài)環(huán)境的關(guān)系隨著時(shí)間與空間的演變呈現(xiàn)差異性，總體為倒 “U” 型結(jié)構(gòu)，過載的經(jīng)濟(jì)開發(fā)必然會(huì)遭到生態(tài)系統(tǒng)的限制[10-12]。在研究尺度上，有學(xué)者基于國(guó)家層面對(duì)全球35 個(gè)貧困國(guó)家旅游與經(jīng)濟(jì)的耦合協(xié)調(diào)水平及其成因進(jìn)行了分析；部分研究從新疆[13]、湖南[9]、甘肅[7]等省域或市域?qū)用娣治雎糜萎a(chǎn)業(yè)和經(jīng)濟(jì)水平、生態(tài)環(huán)境與旅游產(chǎn)業(yè)等之間的發(fā)展關(guān)系；也有從特定的城市群如京津冀[14]、中國(guó)十大城市群[15]、西北五省區(qū)[16]等探討旅游產(chǎn)業(yè)、人居環(huán)境和區(qū)域經(jīng)濟(jì)之間協(xié)調(diào)關(guān)系與時(shí)空演變特征。研究方法包括定性研究和定量評(píng)價(jià)兩大類。學(xué)者們基于數(shù)據(jù)包絡(luò)法[5,17-18]、VAR 模型[19]、VEC 模型[20-21]探索不同系統(tǒng)之間的關(guān)系；基于灰色關(guān)聯(lián)模型[22]預(yù)測(cè)系統(tǒng)之間未來的協(xié)調(diào)發(fā)展程度；基于耦合協(xié)調(diào)模型[23-25]、DPSIR 模型[26-27]等定量分析系統(tǒng)之間耦合協(xié)調(diào)的發(fā)展水平與時(shí)序演變；基于探索性空間數(shù)據(jù)分析、泰爾指數(shù)、核密度分析法等研究不同系統(tǒng)之間耦合協(xié)調(diào)水平的時(shí)空演變。

綜上，目前國(guó)內(nèi)外學(xué)者從多個(gè)角度探討了旅游產(chǎn)業(yè)、生態(tài)環(huán)境與社會(huì)效益之間的關(guān)系，但大多以兩兩分析為主，將三者結(jié)合起來的研究甚少；大多數(shù)研究基于市域、省域、國(guó)家等層面，而從縣域視角出發(fā)的研究較少。本文以河北省長(zhǎng)城沿線59 個(gè)縣域?yàn)檠芯繉?duì)象，將生態(tài)環(huán)境、旅游產(chǎn)業(yè)與社會(huì)效益三者綜合考慮，構(gòu)建了旅游產(chǎn)業(yè)—生態(tài)環(huán)境—社會(huì)效益系統(tǒng)（以下統(tǒng)稱TES 系統(tǒng)），對(duì)相關(guān)理論是一種有益的補(bǔ)充。從研究方法來看，目前缺乏從空間效應(yīng)角度，運(yùn)用空間計(jì)量經(jīng)濟(jì)模型探討耦合協(xié)調(diào)水平的重要影響因素，在方法上具有一定的創(chuàng)新性。本文將為長(zhǎng)城沿線旅游產(chǎn)業(yè)、生態(tài)環(huán)境與社會(huì)效益的協(xié)調(diào)可持續(xù)發(fā)展提供理論依據(jù)，以期推動(dòng)各縣建立長(zhǎng)城保護(hù)與傳承利用的協(xié)同機(jī)制，加強(qiáng)縣域間合理的合作與競(jìng)爭(zhēng)，從而促進(jìn)產(chǎn)業(yè)、生態(tài)和效益的協(xié)調(diào)發(fā)展。

1 研究設(shè)計(jì)

1.1 研究區(qū)域概況

本文選擇河北省為研究對(duì)象。作為長(zhǎng)城國(guó)家文化公園重點(diǎn)建設(shè)區(qū)，河北省長(zhǎng)城資源比例占全國(guó)11.79%，資源認(rèn)定總量位居全國(guó)第二位，是中國(guó)長(zhǎng)城保存最完整、最具代表性的區(qū)段之一，肩負(fù)著引領(lǐng)和示范全國(guó)其他?。ㄊ校﹪?guó)家文化公園建設(shè)的重?fù)?dān)。河北省現(xiàn)存戰(zhàn)國(guó)、漢、北魏、北齊、唐、金、明代等不同時(shí)期的長(zhǎng)城，各個(gè)時(shí)代長(zhǎng)城分布于秦皇島、唐山、承德、張家口、保定、廊坊、石家莊、邢臺(tái)、邯鄲及雄安新區(qū)的59 個(gè)縣域，基本涉及河北省全境，總長(zhǎng)度2 498.54 千米。長(zhǎng)城途經(jīng)縣域海拔為500 ～2 000 米，年降水量為400 ～600 毫米，屬于大陸性季風(fēng)氣候，四季分明，溫差較大，孕育了獨(dú)特的生態(tài)系統(tǒng)。長(zhǎng)城沿線的大部分縣域地處山區(qū)，經(jīng)濟(jì)和社會(huì)發(fā)展情況滯后，生態(tài)環(huán)境脆弱，而當(dāng)?shù)鼐用窈凸芾碚邔?duì)長(zhǎng)城的保護(hù)意識(shí)不足，并未對(duì)其進(jìn)行有效的管理與保護(hù)，對(duì)短期利益的盲目追求使得長(zhǎng)城遭到了嚴(yán)重破壞。因此，選擇此區(qū)域旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益間內(nèi)部互動(dòng)關(guān)系的研究具有典型性，有利于促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)與生態(tài)環(huán)境良好協(xié)調(diào)發(fā)展，同時(shí)對(duì)長(zhǎng)城的保護(hù)與傳承、長(zhǎng)城國(guó)家文化公園的建設(shè)有重要意義[28]。

1.2 指標(biāo)體系構(gòu)建

基于耦合協(xié)調(diào)作用機(jī)理，借鑒相關(guān)文獻(xiàn)資料[4,13-14,29]，結(jié)合河北省長(zhǎng)城沿線縣域的發(fā)展現(xiàn)狀及數(shù)據(jù)可獲取性，從旅游產(chǎn)業(yè)子系統(tǒng)、生態(tài)環(huán)境子系統(tǒng)和社會(huì)效益子系統(tǒng)三個(gè)維度出發(fā)，最終構(gòu)建了包含17 個(gè)二級(jí)指標(biāo)的耦合協(xié)調(diào)度評(píng)價(jià)指標(biāo)體系（表1）。

表1 長(zhǎng)城TES系統(tǒng)評(píng)價(jià)指標(biāo)體系

旅游產(chǎn)業(yè)子系統(tǒng)選取7 項(xiàng)指標(biāo)[30]，包括旅游總收入、人均旅游收入、景區(qū)發(fā)展?fàn)顟B(tài)集聚熵、長(zhǎng)城文化資源集聚熵、星級(jí)飯店數(shù)量、A 級(jí)景區(qū)數(shù)量和公路網(wǎng)密度。其中，長(zhǎng)城文化資源集聚熵以各縣域現(xiàn)存的長(zhǎng)城各類遺址遺跡數(shù)量，如烽火臺(tái)、馬面、敵臺(tái)、關(guān)堡和墻體數(shù)目來代表；景區(qū)發(fā)展?fàn)顟B(tài)集聚熵是以旅游總收入、旅游資源及長(zhǎng)城文化資源三個(gè)方面測(cè)算旅游產(chǎn)業(yè)區(qū)位熵指數(shù)，并取其平均值作為區(qū)位熵指數(shù)來衡量旅游產(chǎn)業(yè)集聚程度。這些指標(biāo)主要反映以長(zhǎng)城為代表的旅游產(chǎn)業(yè)的市場(chǎng)規(guī)模、經(jīng)濟(jì)效益和要素結(jié)構(gòu)。

生態(tài)環(huán)境子系統(tǒng)選取4 項(xiàng)指標(biāo)[31]，包括森林覆蓋率、環(huán)境空氣質(zhì)量綜合指數(shù)、全社會(huì)用電量和水環(huán)境質(zhì)量。其中，森林覆蓋率在生態(tài)子系統(tǒng)中的占比較大，反映了環(huán)境效益對(duì)旅游產(chǎn)業(yè)的狀態(tài)和響應(yīng)。河北省自2015 年以來，森林覆蓋率由31%提高到35%，成為全國(guó)創(chuàng)建森林城市數(shù)量最多的省份。

社會(huì)效益子系統(tǒng)選取6 項(xiàng)指標(biāo)[32]，包括城鎮(zhèn)化率、人均收入、年末參加基本醫(yī)療保險(xiǎn)人數(shù)、城鎮(zhèn)登記失業(yè)率、醫(yī)療衛(wèi)生機(jī)構(gòu)床位數(shù)和普通中學(xué)專任教師。

1.3 研究方法

1.3.1 綜合評(píng)價(jià)指數(shù)

本文采用綜合評(píng)價(jià)模型[33]來測(cè)算旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益三個(gè)子系統(tǒng)的綜合評(píng)價(jià)指數(shù)，具體公式如下：

式中：F1表示旅游產(chǎn)業(yè)子系統(tǒng)的綜合評(píng)價(jià)指數(shù)；F2表示生態(tài)環(huán)境子系統(tǒng)的綜合評(píng)價(jià)指數(shù)；F3表示社會(huì)效益子系統(tǒng)的綜合評(píng)價(jià)指數(shù)：xi、yj、zk分別表示旅游產(chǎn)業(yè)子系統(tǒng)、生態(tài)環(huán)境子系統(tǒng)和社會(huì)效益子系統(tǒng)各要素標(biāo)準(zhǔn)化后的值：αi、bj、ck分別表示旅游產(chǎn)業(yè)子系統(tǒng)、生態(tài)環(huán)境子系統(tǒng)和社會(huì)效益子系統(tǒng)中各指標(biāo)的權(quán)重；n、m、p分別表示旅游產(chǎn)業(yè)子系統(tǒng)、生態(tài)環(huán)境子系統(tǒng)和社會(huì)效益子系統(tǒng)的要素個(gè)數(shù)。

1.3.2 耦合度模型及協(xié)調(diào)度模型

本文借助物理學(xué)中耦合協(xié)調(diào)模型，得到了TES 系統(tǒng)中三個(gè)子系統(tǒng)之間的耦合度與協(xié)調(diào)度[29]：

式中：C表示TES 系統(tǒng)的耦合度，取值范圍為[0, 1]。借鑒中值分段法，將耦合度進(jìn)行類型劃分[34]，標(biāo)準(zhǔn)如下：(0.00, 0.30]表示低水平耦合階段，(0.30, 0.50]表示拮抗階段，(0.50, 0.80]表示磨合階段，(0.80, 1.00]表示高水平耦合階段。D表示TES 系統(tǒng)的耦合協(xié)調(diào)度，T為TES 系統(tǒng)的綜合發(fā)展指數(shù)，αi為待定系數(shù)，其和為1?？紤]到本文主要研究對(duì)象為基于長(zhǎng)城旅游資源及其他旅游資源的TES 系統(tǒng)，因此取α1=0.4、α2=0.3、α3=0.3。參考相關(guān)研究成果[35]，將TES 系統(tǒng)耦合協(xié)調(diào)度進(jìn)行等級(jí)劃分（表2）。

表2 TES系統(tǒng)耦合協(xié)調(diào)度等級(jí)劃分

1.3.3 探索性空間數(shù)據(jù)分析

探索性空間數(shù)據(jù)分析是用來探索社會(huì)現(xiàn)象之間空間關(guān)聯(lián)性的方法，并以空間自相關(guān)為測(cè)度[30,36]。本文采用全局和局部莫蘭指數(shù)研究空間關(guān)聯(lián)效應(yīng)，全局莫蘭指數(shù)公式如下：

局部莫蘭指數(shù)分析研究單元空間的不穩(wěn)定性，公式如下：

1.3.4 空間效應(yīng)模型選取

河北省長(zhǎng)城沿線縣域TES 系統(tǒng)耦合協(xié)調(diào)度在相近縣域之間產(chǎn)生一定的空間交互作用。考慮空間效應(yīng)，構(gòu)建如下三種空間面板模型。

（1）空間滯后模型（SAR）。

式中：空間滯后自回歸系數(shù)為ρ，空間權(quán)重矩陣經(jīng)過行標(biāo)準(zhǔn)化處理后為W，解釋變量的空間滯后項(xiàng)為Wyjt。

（2）空間誤差模型（SEM）。

式中：相鄰縣域觀測(cè)值的誤差項(xiàng)對(duì)本縣域影響程度用ρ表示，表示空間滯后誤差變量。

（3）空間杜賓模型（SDM）。

空間杜賓模型是上述兩種空間面板模型的綜合，認(rèn)為研究單元的耦合協(xié)調(diào)度還與鄰近縣域之間存在空間關(guān)聯(lián)。

1.4 數(shù)據(jù)來源與說明

本文選取河北省長(zhǎng)城沿線2015—2019 年59 個(gè)縣域的面板數(shù)據(jù)。其中，公路網(wǎng)密度、森林覆蓋率、環(huán)境空氣質(zhì)量綜合指數(shù)、全社會(huì)用電量、水環(huán)境質(zhì)量、城鎮(zhèn)化率、人均收入、年末參加基本醫(yī)療保險(xiǎn)人數(shù)、城鎮(zhèn)登記失業(yè)率、醫(yī)療衛(wèi)生機(jī)構(gòu)床位數(shù)、普通中學(xué)專任教師、年末總?cè)丝凇⒌貐^(qū)面積、第三產(chǎn)業(yè)生產(chǎn)總值、出口貿(mào)易總額、技術(shù)人員數(shù)、第二產(chǎn)業(yè)總產(chǎn)值等數(shù)據(jù)主要來源于2016—2020 年的《河北經(jīng)濟(jì)年鑒》、各縣國(guó)民經(jīng)濟(jì)發(fā)展統(tǒng)計(jì)公報(bào)、市級(jí)人民政府官網(wǎng)；旅游總收入、人均旅游收入、長(zhǎng)城各類遺址遺跡數(shù)量、長(zhǎng)城文化資源、長(zhǎng)城旅游資源、星級(jí)飯店數(shù)量和A 級(jí)景區(qū)數(shù)量等數(shù)據(jù)來源于市級(jí)文化和旅游部門官網(wǎng)、縣級(jí)人民政府官網(wǎng)、縣級(jí)文化和旅游部門網(wǎng)站等。對(duì)少數(shù)數(shù)據(jù)缺失情況則采用插值法進(jìn)行補(bǔ)充，并對(duì)所有數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理以消除數(shù)據(jù)的量綱和量級(jí)。

2 實(shí)證分析

2.1 TES系統(tǒng)綜合發(fā)展水平

根據(jù)熵值法得出河北省長(zhǎng)城沿線縣域TES 系統(tǒng)與三大子系統(tǒng)的綜合評(píng)價(jià)指數(shù)（圖1），整體上看，TES 系統(tǒng)的綜合發(fā)展水平較高，且呈現(xiàn)上升趨勢(shì)，取得了良好的成效。同時(shí)，三大子系統(tǒng)的發(fā)展水平也呈現(xiàn)上升趨勢(shì)，且社會(huì)效益子系統(tǒng)＞旅游產(chǎn)業(yè)子系統(tǒng)＞生態(tài)環(huán)境子系統(tǒng)；社會(huì)效益子系統(tǒng)穩(wěn)中有升；旅游產(chǎn)業(yè)子系統(tǒng)雖呈上升趨勢(shì)，但并不太明顯；生態(tài)環(huán)境子系統(tǒng)數(shù)值相對(duì)較低，且有下降趨勢(shì)，說明對(duì)環(huán)境的重視程度仍有待提高，對(duì)長(zhǎng)城沿線旅游資源的開發(fā)不能忽視保護(hù)環(huán)境。

圖1 2015—2019年河北省長(zhǎng)城沿線縣域TES系統(tǒng)、子系統(tǒng)綜合發(fā)展水平

2.2 TES系統(tǒng)耦合協(xié)調(diào)度特征分析

2.2.1 耦合協(xié)調(diào)度分析

經(jīng)測(cè)算，85%的縣域TES 系統(tǒng)間的耦合度在0.8 以上，表明大部分長(zhǎng)城沿線縣域三大子系統(tǒng)間實(shí)現(xiàn)了高度耦合，但仍有少數(shù)縣域系統(tǒng)間的耦合程度較低，三大子系統(tǒng)間的耦合存在較大進(jìn)步空間。

進(jìn)一步運(yùn)用耦合協(xié)調(diào)度模型測(cè)算三大系統(tǒng)間協(xié)調(diào)狀況，得到長(zhǎng)城沿線縣域2015—2019 年的TES 系統(tǒng)耦合協(xié)調(diào)度指數(shù)，均值為0.519 4，整體上三大子系統(tǒng)處于勉強(qiáng)協(xié)調(diào)狀態(tài)，但數(shù)值呈現(xiàn)增長(zhǎng)趨勢(shì)。耦合協(xié)調(diào)度的最大值為0.881，為遷安市，因其具有較高的經(jīng)濟(jì)發(fā)展水平，大力發(fā)展旅游產(chǎn)業(yè)，持續(xù)改善生態(tài)環(huán)境，社會(huì)效益也逐年提升，三個(gè)子系統(tǒng)處于良好協(xié)調(diào)狀態(tài)；最小值為0.281，為安新縣，處于中度失調(diào)狀態(tài)，因其生態(tài)環(huán)境較好，但旅游產(chǎn)業(yè)發(fā)展尚不成熟，對(duì)社會(huì)效益貢獻(xiàn)不足。

2.2.2 耦合協(xié)調(diào)度時(shí)空分異特征

為進(jìn)一步研究長(zhǎng)城沿線縣域TES 系統(tǒng)耦合協(xié)調(diào)水平的空間分異特征，本文選取三個(gè)子系統(tǒng)2015 年、2017年和2019 年的耦合協(xié)調(diào)度進(jìn)行分析。發(fā)現(xiàn)59 個(gè)縣域的耦合協(xié)調(diào)度呈現(xiàn)出中度失調(diào)、輕度失調(diào)、瀕臨失調(diào)、勉強(qiáng)協(xié)調(diào)、初級(jí)協(xié)調(diào)、中級(jí)協(xié)調(diào)和良好協(xié)調(diào)7 種類型，前三類為失調(diào)狀態(tài)，后四類為協(xié)調(diào)狀態(tài)。從整體來看，中度失調(diào)和良好協(xié)調(diào)比重最少，輕度失調(diào)和初級(jí)協(xié)調(diào)占比較少，瀕臨失調(diào)和勉強(qiáng)協(xié)調(diào)比重最高，占到66%以上，呈現(xiàn)出兩頭少中間多的橄欖型結(jié)構(gòu)特征。

從時(shí)間上看（表3），整體上呈現(xiàn)由失調(diào)狀態(tài)向協(xié)調(diào)狀態(tài)逐步演化的趨勢(shì)。2015 年處于失調(diào)狀態(tài)的縣域有26個(gè)，總體占比為44.1%；中度失調(diào)的縣域有3 個(gè)，輕度失調(diào)的縣域有8 個(gè)，瀕臨失調(diào)的縣域有15 個(gè)，其中中度失調(diào)的縣域?yàn)樯辛x縣、雙灤區(qū)和安新縣；處于協(xié)調(diào)狀態(tài)的縣域有33 個(gè)，總體占比為55.9%；勉強(qiáng)協(xié)調(diào)的縣域有24 個(gè)，初級(jí)協(xié)調(diào)的縣域有5 個(gè)，中級(jí)協(xié)調(diào)的縣域有2 個(gè)：遷西縣和武安市，良好協(xié)調(diào)的縣域有2 個(gè)：遵化市和遷安市。與2015 年相比，2017 年失調(diào)狀態(tài)的縣域數(shù)量與協(xié)調(diào)狀態(tài)的縣域數(shù)量相同。不同的是，中度失調(diào)的縣域數(shù)減少1 個(gè)，雙灤區(qū)和尚義縣轉(zhuǎn)化為輕度失調(diào)狀態(tài)，雄縣降為中度失調(diào)；勉強(qiáng)協(xié)調(diào)的縣域數(shù)量減少3 個(gè)，沙河市、寬城滿族自治縣和圍場(chǎng)滿族蒙古族自治縣均演變?yōu)槌跫?jí)協(xié)調(diào)。與2017 年相比，2019 年失調(diào)狀態(tài)的縣域數(shù)量減少4 個(gè)，肥鄉(xiāng)區(qū)、唐縣、雙橋區(qū)和靈壽縣均轉(zhuǎn)化為協(xié)調(diào)狀態(tài)。從數(shù)量上可見，各縣域積極開發(fā)旅游產(chǎn)業(yè)，促進(jìn)經(jīng)濟(jì)提升與生態(tài)環(huán)境協(xié)調(diào)發(fā)展，從而減小了三大系統(tǒng)之間的差異，長(zhǎng)城沿線縣域TES 系統(tǒng)的耦合協(xié)調(diào)性不斷增強(qiáng)，失調(diào)狀態(tài)轉(zhuǎn)向協(xié)調(diào)狀態(tài)的縣域不斷增多。

表3 2015年、2017年、2019年TES系統(tǒng)耦合協(xié)調(diào)發(fā)展的時(shí)間趨勢(shì)

從空間分布上看（表3），東北部和南部的縣域耦合協(xié)調(diào)度均值分別為0.597 3 和0.541 6，處于協(xié)調(diào)狀態(tài)的縣域占比分別達(dá)到77%和65%，TES 系統(tǒng)發(fā)展良好，說明這些縣域旅游產(chǎn)業(yè)、社會(huì)環(huán)境與社會(huì)效益三者協(xié)調(diào)發(fā)展；中部和西北部的縣域耦合協(xié)調(diào)度均值分別為0.472 6和0.425 9，且處于協(xié)調(diào)階段的縣域僅占到了21%和29%，整體上處于瀕臨失調(diào)狀態(tài)，大部分縣域失調(diào)的原因是較好的生態(tài)環(huán)境與低水平的旅游產(chǎn)業(yè)、社會(huì)效益之間形成了較大的差距，三大子系統(tǒng)互相束縛與制約，因此TES 系統(tǒng)的耦合協(xié)調(diào)度有待改善。耦合協(xié)調(diào)度的地域分布格局呈現(xiàn)東北部＞南部＞中部＞西北部的態(tài)勢(shì)。

綜合來看，協(xié)調(diào)狀態(tài)由低層次向高層次轉(zhuǎn)化的縣域主要位于河北省中部和東北部。在2015 年、2017 年和2019 年始終處于中度失調(diào)的為安新縣，始終處于輕度失調(diào)的為容城縣、下花園區(qū)和陽原縣，始終處于初級(jí)協(xié)調(diào)的有懷來縣和橋西區(qū)，這五個(gè)縣域位于河北省的中部和西北部地區(qū)；始終處于良好協(xié)調(diào)的為東北部的遵化市和遷安市。協(xié)調(diào)狀態(tài)由低層次向高層次的轉(zhuǎn)變的有雄縣和尚義縣，由中度失調(diào)轉(zhuǎn)變?yōu)檩p度失調(diào)；位于河北中部的贊皇縣、阜平縣和順平縣等由輕度失調(diào)轉(zhuǎn)為瀕臨失調(diào)；位于河北東北部的圍場(chǎng)滿族蒙古族自治縣、隆化縣、撫寧區(qū)和隆化縣等由勉強(qiáng)協(xié)調(diào)轉(zhuǎn)為初級(jí)協(xié)調(diào)。周邊縣域較高的協(xié)調(diào)水平對(duì)本地良性轉(zhuǎn)變產(chǎn)生了一定的影響，這也反映出旅游產(chǎn)業(yè)發(fā)展和生態(tài)環(huán)境的有效治理產(chǎn)生了良性協(xié)同效應(yīng)。

3 耦合協(xié)調(diào)度空間回歸分析

3.1 全局空間自相關(guān)分析

研究區(qū)域耦合協(xié)調(diào)度的全局Moran’s I（表4），均通過了1%的顯著性水平檢驗(yàn)，且數(shù)值在0.3 ～0.5 的區(qū)間內(nèi)變化，說明耦合協(xié)調(diào)度具有空間正相關(guān)關(guān)系。指數(shù)呈現(xiàn)反復(fù)波動(dòng)趨勢(shì)表明相鄰地區(qū)TES 系統(tǒng)耦合協(xié)調(diào)度空間自相關(guān)性不穩(wěn)定，總體下降的趨勢(shì)表明相鄰地區(qū)耦合協(xié)調(diào)水平逐漸向均衡趨勢(shì)發(fā)展。

表4 耦合協(xié)調(diào)度的全局Moran’s I

3.2 局部空間自相關(guān)分析

全局Moran’s I 顯示整體呈現(xiàn)空間自相關(guān)性，而局部Moran’s I 則反映各縣域之間的空間關(guān)聯(lián)性。首先選取2015 年、2017 年和2019 年三個(gè)時(shí)間斷面的局部Moran’s I，數(shù)值分別為0.379、0.347、0.295，呈下降趨勢(shì)且均大于0.25，說明各縣域間耦合協(xié)調(diào)度存在空間關(guān)聯(lián)性。其次輸出河北省長(zhǎng)城沿線縣域耦合協(xié)調(diào)水平的LISA 集聚結(jié)果（表5）以顯示各縣域的相關(guān)類型及其聚集區(qū)相關(guān)的顯著性，其中四個(gè)象限的縣域均通過了5%的顯著性水平檢驗(yàn)。

表5 LISA集聚結(jié)果

從集聚效果來看，HH 象限的縣域主要為興隆縣、青龍滿族自治縣、遷西縣和盧龍縣，涉縣只在2015 年分布在此象限，寬城滿族自治縣只在2017 年分布在此象限，表明這些縣域的耦合協(xié)調(diào)水平較高，且其周邊縣域的耦合協(xié)調(diào)水平也較高。HL 象限的縣域主要為張北縣，表明其自身的耦合協(xié)調(diào)水平較高，而周邊縣域的發(fā)展水平較低。LL 象限主要為容城縣、雄縣和安新縣，表明這些縣域不僅自身發(fā)展較弱，且周邊縣域發(fā)展也較落后，是需要重點(diǎn)發(fā)展扶持的區(qū)域。徐水區(qū)的變動(dòng)最大，從2015 年的HL 象限到LL 象限，表明徐水區(qū)自身耦合協(xié)調(diào)水平逐漸減弱。LH 象限的縣域顯著性不強(qiáng)，主要為2015 年的肥鄉(xiāng)區(qū)，縣域自身耦合協(xié)調(diào)度發(fā)展水平低，而周邊縣域耦合協(xié)調(diào)水平高，通過周邊高水平縣域的帶動(dòng)作用，逐漸有所發(fā)展。

3.3 空間效應(yīng)分析

耦合協(xié)調(diào)度并不能夠深入分析這三大子系統(tǒng)是如何作用于TES 系統(tǒng)以及作用的程度，為分析河北省長(zhǎng)城沿線縣域TES 系統(tǒng)耦合水平差異的影響因素，在借鑒前人研究成果的基礎(chǔ)上，結(jié)合河北省長(zhǎng)城沿線縣域的具體發(fā)展情況，本文將考察如下解釋變量對(duì)耦合協(xié)調(diào)度的影響程度：人口密度（POP）[37]，以年末總?cè)丝谂c地區(qū)面積的比值進(jìn)行衡量；經(jīng)濟(jì)發(fā)展水平（ECO）[38]，以第三產(chǎn)業(yè)生產(chǎn)總值來衡量；對(duì)外開放程度（OPEN）[39]，以出口貿(mào)易總額占GDP 的比重來衡量；技術(shù)創(chuàng)新（TEC）[40]，以技術(shù)人員數(shù)來衡量；產(chǎn)業(yè)結(jié)構(gòu)（STR）[41-42]，第二產(chǎn)業(yè)總產(chǎn)值占GDP 的比重，對(duì)數(shù)據(jù)進(jìn)行預(yù)處理以減少異方差性對(duì)結(jié)果的影響。

對(duì)空間計(jì)量模型進(jìn)行選擇（表6），豪斯曼檢驗(yàn)在1%的顯著水平下拒絕了原假設(shè)，說明在河北省長(zhǎng)城沿線縣域TES 系統(tǒng)耦合協(xié)調(diào)度空間效應(yīng)研究中固定效應(yīng)優(yōu)于隨機(jī)效應(yīng)；采用拉格朗日乘數(shù)檢驗(yàn)（LM），空間滯后模型與空間誤差模型的LM 檢驗(yàn)與Robust LM 檢驗(yàn)均通過了10%的顯著性檢驗(yàn)，拒絕了原假設(shè)，說明河北省長(zhǎng)城沿線縣域TES 系統(tǒng)耦合協(xié)調(diào)度具有空間相關(guān)性；通過Wald檢驗(yàn)與LR 檢驗(yàn)發(fā)現(xiàn)在10%的顯著水平下拒絕了空間杜賓模型簡(jiǎn)化為空間滯后模型的假設(shè)，Wald 檢驗(yàn)在1%的顯著水平下拒絕了簡(jiǎn)化為空間誤差模型的假設(shè)，因此，基于固定效應(yīng)的空間杜賓模型為本文研究的最佳選擇。

表6 空間面板計(jì)量模型的檢驗(yàn)結(jié)果

本文將不同模型的固定效應(yīng)進(jìn)行了比對(duì)，檢驗(yàn)了模型的穩(wěn)健性。根據(jù)空間誤差模型、空間滯后模型和空間杜賓模型的估計(jì)數(shù)據(jù)可以看出（表7），影響河北省長(zhǎng)城沿線縣域TES 系統(tǒng)耦合協(xié)調(diào)水平的因子存在空間關(guān)系，本縣域TES 系統(tǒng)耦合協(xié)調(diào)度的影響因子也會(huì)影響鄰近縣域TES 系統(tǒng)耦合協(xié)調(diào)水平。

表7 空間面板模型的估計(jì)結(jié)果

加權(quán)人口密度和產(chǎn)業(yè)結(jié)構(gòu)的回歸系數(shù)分別為-0.634和-0.303，且在空間杜賓模型中顯著為負(fù)，在河北省長(zhǎng)城沿線縣域發(fā)展過程中，人口密度越大，可能帶來的資源供應(yīng)不足和環(huán)境破壞等問題越明顯，不利于耦合協(xié)調(diào)水平的提高；同時(shí)，第二產(chǎn)業(yè)的比重越高，給環(huán)境帶來一定負(fù)面影響，也會(huì)對(duì)耦合協(xié)調(diào)度產(chǎn)生負(fù)向作用；經(jīng)濟(jì)發(fā)展水平對(duì)耦合協(xié)調(diào)度的回歸系數(shù)高達(dá)0.06，通過了1%的顯著性檢驗(yàn)，表明經(jīng)濟(jì)發(fā)展水平是增加河北省長(zhǎng)城沿線縣域TES 系統(tǒng)耦合協(xié)調(diào)度的重要因素，同時(shí)也說明本縣域經(jīng)濟(jì)發(fā)展水平的提升在一定程度上也會(huì)促進(jìn)鄰近縣域耦合協(xié)調(diào)水平的提升；對(duì)外開放程度、技術(shù)創(chuàng)新的回歸系數(shù)并不顯著，表明這些因素對(duì)耦合協(xié)調(diào)水平提升的影響較弱，這是由于長(zhǎng)城沿線的縣域發(fā)展水平較低，對(duì)外開放水平較低，對(duì)技術(shù)創(chuàng)新的重視程度也較低，對(duì)耦合協(xié)調(diào)水平的影響較不明顯。

為了進(jìn)一步厘清空間影響機(jī)制，將空間杜賓模型中的空間效應(yīng)進(jìn)一步分解為直接效應(yīng)、間接效應(yīng)和總效應(yīng)（表8）。人口密度的間接效應(yīng)和總效應(yīng)均顯著為負(fù)，說明本縣域的人口密度不僅會(huì)負(fù)向影響當(dāng)?shù)氐腡ES 系統(tǒng)耦合協(xié)調(diào)水平，同時(shí)也會(huì)對(duì)鄰近縣域產(chǎn)生較大的負(fù)面作用；經(jīng)濟(jì)發(fā)展水平在直接效應(yīng)和總效應(yīng)中均顯著為正，經(jīng)濟(jì)發(fā)展水平越高，民眾的環(huán)境保護(hù)意識(shí)越強(qiáng)，對(duì)生活質(zhì)量的要求也會(huì)越高，旅游相關(guān)產(chǎn)業(yè)也會(huì)根據(jù)市場(chǎng)的需求進(jìn)行改進(jìn)與優(yōu)化，從而促進(jìn)TES 系統(tǒng)耦合協(xié)調(diào)水平的提高；技術(shù)創(chuàng)新在總效應(yīng)中顯著為正，說明在提高TES 系統(tǒng)耦合協(xié)調(diào)水平的過程中，技術(shù)的進(jìn)步會(huì)通過勞動(dòng)生產(chǎn)率的提升、企業(yè)結(jié)構(gòu)的優(yōu)化等方面自然地促進(jìn)經(jīng)濟(jì)增長(zhǎng)，減少對(duì)環(huán)境的影響，促進(jìn)TES 系統(tǒng)耦合協(xié)調(diào)水平的提升；產(chǎn)業(yè)結(jié)構(gòu)在間接效應(yīng)和總效應(yīng)中顯著為負(fù)，第二產(chǎn)業(yè)占比較大不僅會(huì)對(duì)本縣域環(huán)境質(zhì)量造成影響，也會(huì)對(duì)鄰近縣域產(chǎn)生一定的影響。

4 結(jié)論與建議

在建立旅游產(chǎn)業(yè)、生態(tài)環(huán)境、社會(huì)效益系統(tǒng)評(píng)價(jià)指標(biāo)體系的基礎(chǔ)上，對(duì)河北省長(zhǎng)城沿線縣域2015—2019 年旅游產(chǎn)業(yè)、生態(tài)環(huán)境與社會(huì)效益三個(gè)子系統(tǒng)之間協(xié)調(diào)發(fā)展的時(shí)空演化特征與空間效應(yīng)進(jìn)行分析，得出以下結(jié)論。

（1）2015—2019 年河北省長(zhǎng)城沿線縣域TES 系統(tǒng)總體發(fā)展的綜合得分呈現(xiàn)上升趨勢(shì)，發(fā)展?fàn)顟B(tài)穩(wěn)步向好。長(zhǎng)城沿線縣域社會(huì)效益子系統(tǒng)和旅游產(chǎn)業(yè)子系統(tǒng)發(fā)展水平總體上呈現(xiàn)出緩慢上升趨勢(shì)，生態(tài)環(huán)境子系統(tǒng)發(fā)展水平較低，增長(zhǎng)趨勢(shì)不明顯，仍有很大的改善空間。

（2）河北省長(zhǎng)城沿線大部分縣域TES 系統(tǒng)處于高耦合低協(xié)調(diào)狀態(tài)。近年來轉(zhuǎn)向協(xié)調(diào)狀態(tài)的縣域數(shù)量有所上升，隨著時(shí)間的推移，總體上三大子系統(tǒng)間的耦合協(xié)調(diào)關(guān)系得到一定程度改善；不同縣域之間耦合協(xié)調(diào)水平呈現(xiàn)出空間性差異，地域分布格局呈現(xiàn)東北部＞南部＞中部＞西北部的態(tài)勢(shì)。

（3）人口密度、經(jīng)濟(jì)發(fā)展水平、對(duì)外開放程度、技術(shù)創(chuàng)新、產(chǎn)業(yè)結(jié)構(gòu)在不同程度上，以不同的方式影響著TES 系統(tǒng)耦合協(xié)調(diào)水平，在某種程度上，這些影響因子不僅會(huì)影響縣域本身TES 系統(tǒng)耦合協(xié)調(diào)水平，而且也對(duì)相鄰城市TES 系統(tǒng)耦合協(xié)調(diào)發(fā)展水平產(chǎn)生影響。

通過對(duì)TES 系統(tǒng)耦合協(xié)調(diào)性時(shí)空格局與空間效應(yīng)研究的結(jié)論分析，提出以下建議。

（1）強(qiáng)力推進(jìn)河北省長(zhǎng)城沿線低耦合協(xié)調(diào)水平縣域旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益的協(xié)同發(fā)展，崇禮區(qū)、尚義縣、赤城縣、沽源縣和康?？h等縣（區(qū)）應(yīng)強(qiáng)化旅游產(chǎn)業(yè)對(duì)經(jīng)濟(jì)的帶動(dòng)作用，在不損害環(huán)境質(zhì)量的基礎(chǔ)上，實(shí)現(xiàn)社會(huì)效益的提高；雙灤區(qū)、安新縣、容城等縣域應(yīng)注重生態(tài)環(huán)境的保護(hù)，安新縣、容城縣、雄縣等縣域應(yīng)做好長(zhǎng)城本體及沿線的文化景觀風(fēng)貌的整體性保護(hù)和修復(fù)，對(duì)相關(guān)文化遺產(chǎn)進(jìn)行搶險(xiǎn)加固或修繕性保護(hù)，在保護(hù)基礎(chǔ)上合理開發(fā)利用長(zhǎng)城文化旅游資源，從而帶動(dòng)當(dāng)?shù)芈糜萎a(chǎn)業(yè)發(fā)展，促進(jìn)社會(huì)效益的提高，帶動(dòng)城鄉(xiāng)協(xié)同發(fā)展。

（2）提高 “協(xié)調(diào)縣域” 在旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益方面的輻射與帶動(dòng)的 “乘數(shù)效應(yīng)” 。遷安市、橋西區(qū)、遵化市等縣域旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益系統(tǒng)的協(xié)調(diào)發(fā)展水平較高，應(yīng)利用其豐富的旅游資源和較高的經(jīng)濟(jì)水平，通過 “擴(kuò)散效應(yīng)” 推動(dòng)資源和優(yōu)勢(shì)向其他縣域溢出和輻射，推動(dòng)形成縣域經(jīng)濟(jì)發(fā)展與長(zhǎng)城保護(hù)的聯(lián)動(dòng)機(jī)制，借鑒成熟的管理經(jīng)驗(yàn)，針對(duì)目前面臨的主要問題和較大的潛在威脅采取有針對(duì)性的舉措。

（3）多措并舉促進(jìn)旅游產(chǎn)業(yè)、生態(tài)環(huán)境和社會(huì)效益持續(xù)協(xié)調(diào)發(fā)展。不僅要提升長(zhǎng)城沿線縣域整體協(xié)調(diào)水平，也要促進(jìn)縣域協(xié)調(diào)發(fā)展的可持續(xù)性。因此，以山海關(guān)、金山嶺、大境門和崇禮為核心，借助長(zhǎng)城國(guó)家文化公園發(fā)展的契機(jī)，爭(zhēng)取資金等優(yōu)惠政策，吸引人、財(cái)、物和信息等資源，促進(jìn)長(zhǎng)城沿線公共服務(wù)廊道的建設(shè)，發(fā)揮極點(diǎn)帶動(dòng)作用，促進(jìn)長(zhǎng)城文化旅游帶的整體發(fā)展。