陜穎穎，唐健雄，劉雨靖

(湖南師范大學旅游學院，中國 長沙 410081)

城鄉(xiāng)發(fā)展不平衡與鄉(xiāng)村發(fā)展不充分一直是我國面臨的重大問題,為了解決這一問題，國家相繼提出了新型城鎮(zhèn)化和鄉(xiāng)村振興兩大戰(zhàn)略。在城鄉(xiāng)發(fā)展過程中，城鎮(zhèn)化與鄉(xiāng)村振興產(chǎn)生了共生效應，涉及自然、政治、社會、經(jīng)濟、文化等方面[1]?！?019年新型城鎮(zhèn)化建設重點任務》明確提出要促進“人口就地就近城鎮(zhèn)化”。就地城鎮(zhèn)作為新型城鎮(zhèn)化的一種重要模式也為鄉(xiāng)村振興帶來了全新的解決思路，是促進鄉(xiāng)村振興和新型城鎮(zhèn)化發(fā)展的重要手段。旅游引導的城鎮(zhèn)化(簡稱“旅游城鎮(zhèn)化”，下同)改變了地區(qū)原有經(jīng)濟模式、地域景觀、生活方式和社會結構等，推動區(qū)域非農(nóng)化發(fā)展向小城鎮(zhèn)(市)演化，具有就地城鎮(zhèn)化特點[2]，推動了張家界、鳳凰等地區(qū)脫貧，實現(xiàn)了當?shù)氐泥l(xiāng)村振興[3]。2018年習近平總書記提出：鄉(xiāng)村振興是否成功，關鍵在于農(nóng)村居民是否獲得真正的幸福感[4]。旅游城鎮(zhèn)化作為一種以旅游業(yè)帶動城鎮(zhèn)化發(fā)展的特殊城市化模型，能夠借用旅游產(chǎn)業(yè)的關聯(lián)作用，有效提升農(nóng)村居民的收入水平、受教育程度、居住條件和消費層次等，縮短城鄉(xiāng)差距，提高居民的幸福感[5]。但是目前有不少旅游小鎮(zhèn)陷入發(fā)展困境，面臨產(chǎn)業(yè)規(guī)模小、經(jīng)營低效、品質(zhì)不高等問題，沒有給城鄉(xiāng)統(tǒng)籌帶來積極效果，反而對鄉(xiāng)村振興產(chǎn)生了負面影響。這表明，旅游城鎮(zhèn)化發(fā)展與鄉(xiāng)村振興之間存在復雜的非線性關系，我們需要重新審視旅游城鎮(zhèn)化發(fā)展過程，剖析它與鄉(xiāng)村振興的關系。因此，從旅游城鎮(zhèn)化的視角來分析其與鄉(xiāng)村振興協(xié)調(diào)發(fā)展的模式，理應成為學界共同關注的話題。

旅游城鎮(zhèn)化自2000年提出以來逐漸受到政府和學界的廣泛關注，并且從理論和實證層面證實了旅游是驅動城鎮(zhèn)化的重要力量之一[6]。Mullins 最早提出有關旅游城鎮(zhèn)化概念，認為旅游城鎮(zhèn)化是一種建立在享樂的銷售和消費基礎上的城市化模式[7]。國內(nèi)學者認為旅游城鎮(zhèn)化不僅是中國城鎮(zhèn)化多元道路的一種重要模式，而且旅游城鎮(zhèn)化的發(fā)展能夠對當?shù)亟?jīng)濟、環(huán)境、社會、文化產(chǎn)生一定的影響[8，9]。傳統(tǒng)城鎮(zhèn)化是資本為了經(jīng)濟增長而盲目擴張的城鎮(zhèn)化，是一種缺乏可持續(xù)性的城鎮(zhèn)化，給居住地帶來了環(huán)境污染、耕地浪費以及“農(nóng)村空心化”等一系列問題。區(qū)別于傳統(tǒng)城鎮(zhèn)化，旅游城鎮(zhèn)化表現(xiàn)出旅游產(chǎn)業(yè)聚集、城市經(jīng)濟結構轉變、人居環(huán)境的改善、基礎設施完善等特征，契合中國目前的“以人為核心”的新型城鎮(zhèn)化發(fā)展方向，是新型城鎮(zhèn)化的重要抓手[10]。鄉(xiāng)村振興是當今中國鄉(xiāng)村發(fā)展一個重大戰(zhàn)略指標，如何推動鄉(xiāng)村振興合理、快速、有效地發(fā)展，始終是學術界關注的熱點。自2017年以來，鄉(xiāng)村振興戰(zhàn)略著眼于從根本解決農(nóng)民、農(nóng)村、農(nóng)業(yè)“三農(nóng)”問題，補齊鄉(xiāng)村發(fā)展短板，促進城鄉(xiāng)經(jīng)濟社會均衡發(fā)展和鄉(xiāng)村充分發(fā)展。

那么旅游城鎮(zhèn)化與鄉(xiāng)村振興是如何協(xié)調(diào)發(fā)展呢？國內(nèi)學者大都集中在對新型城鎮(zhèn)化與鄉(xiāng)村振興的發(fā)展的研究上。龐新軍[11]等人認為，城鎮(zhèn)化的發(fā)展過程必然促進鄉(xiāng)村的發(fā)展。新型城鎮(zhèn)化戰(zhàn)略能夠為鄉(xiāng)村振興制造條件，提供人力、物力、財力方面的流通，從而實現(xiàn)城鄉(xiāng)可持續(xù)和健康發(fā)展。徐維祥[12]等人從鄉(xiāng)村振興與新型城鎮(zhèn)化的耦合協(xié)調(diào)入手，認為湖南省鄉(xiāng)村振興與新型城鎮(zhèn)化耦合協(xié)調(diào)呈現(xiàn)向上轉移的態(tài)勢。鄉(xiāng)村振興的目標是使居民獲得幸福感、滿足感，縮小城鄉(xiāng)差距，以實現(xiàn)共同富裕。婁伶俐[13]認為幸福感是對個人或家庭的總體綜合認知的感覺體現(xiàn)，是衡量人們生活質(zhì)量和社會福利的重要指標之一。麻學鋒[14]等通過研究張家界市旅游城鎮(zhèn)化與城市居民幸福的機制，發(fā)現(xiàn)城鎮(zhèn)化集聚水平的提高對城鎮(zhèn)居民幸福水平的提升具有正向促進作用。由此可見，旅游城鎮(zhèn)化作為一種特殊城鎮(zhèn)化模式，對鄉(xiāng)村經(jīng)濟、文化、社會、生態(tài)有很強的關聯(lián)帶動作用，是推動鄉(xiāng)村振興和城鄉(xiāng)融合發(fā)展的重要動力[15]。

綜上所述，旅游城鎮(zhèn)化與鄉(xiāng)村振興之間并不是一種簡單的線性演化關系。在研究方法上，傳統(tǒng)的耦合協(xié)調(diào)模型只能片面研究兩者之間的關系，缺乏對旅游城鎮(zhèn)化與鄉(xiāng)村振興發(fā)展趨勢的預測，并且不能揭示其非線性機制的演化趨勢。在研究內(nèi)容方面，旅游城鎮(zhèn)化和鄉(xiāng)村振興兩個獨立體系的研究成果相對豐富，但對旅游城鎮(zhèn)化和鄉(xiāng)村振興兩個體系之間的協(xié)調(diào)發(fā)展關系、時空演變規(guī)律以及驅動機制研究還有待探索。基于此，本文選取張家界市為研究區(qū)域，采取SD模型仿真的方法，建立旅游城鎮(zhèn)化與鄉(xiāng)村振興的存量流量圖，依據(jù)張家界“十四五”規(guī)劃進行參數(shù)調(diào)節(jié)，深入探討旅游城鎮(zhèn)化與鄉(xiāng)村振興的作用機制，對不同情景下旅游城鎮(zhèn)化與鄉(xiāng)村振興的關系機制和演化趨勢進行仿真模擬，旨在揭示旅游城鎮(zhèn)化發(fā)展背景下鄉(xiāng)村振興的影響機理，并為張家界市提出未來鄉(xiāng)村振興發(fā)展的政策意見。

1 研究區(qū)域、系統(tǒng)結構與數(shù)據(jù)來源

1.1 研究區(qū)域概況

武陵山片區(qū)作為典型的全國集中連片特困區(qū)域，是我國扶貧的重點區(qū)域之一，片區(qū)內(nèi)擁有豐富的自然資源，因此旅游產(chǎn)業(yè)成為該區(qū)域經(jīng)濟發(fā)展的重要驅動力。張家界市位于湖南省西北部，處于武陵山脈腹地，現(xiàn)轄兩區(qū)兩縣(武陵源區(qū)、永定區(qū)、桑植縣以及慈利縣)，是我國內(nèi)陸一個偏遠狹小的城市。1982年，張家界建立第一個國家森林公園，以旅游立市，不斷發(fā)展成為國家重點旅游城市之一，在2020中國旅游城市排行榜排名第46名。近年來，隨著旅游產(chǎn)業(yè)的不斷發(fā)展，帶動張家界市經(jīng)濟水平不斷升高。截止2020年底，張家界市旅游收入占全市GDP約69.57%，城鎮(zhèn)居民人均可支配收入達到27 884元，農(nóng)村人均可支配收入11 538元，貧困發(fā)生率逐漸降低，旅游業(yè)已成為推動當?shù)孛撠氈赂?，實現(xiàn)鄉(xiāng)村振興的重要途徑。在旅游產(chǎn)業(yè)的驅動下，城鎮(zhèn)化水平達到51.65%，成為旅游產(chǎn)業(yè)推動城鎮(zhèn)化水平發(fā)展的典型城市?；诖?，本文以張家界市為研究案例地，具有一定的典型性和代表性，以期能夠為張家界以及片區(qū)內(nèi)同類型旅游地區(qū)鄉(xiāng)村振興的發(fā)展提供參考。

1.2 旅游城鎮(zhèn)化與鄉(xiāng)村振興協(xié)調(diào)發(fā)展的作用機理

旅游城鎮(zhèn)化是一種以旅游業(yè)帶動城鎮(zhèn)化發(fā)展的新型城鎮(zhèn)化模式，由旅游業(yè)、城鎮(zhèn)化與生態(tài)環(huán)境三個系統(tǒng)構成[16]。旅游產(chǎn)業(yè)是鄉(xiāng)村地區(qū)實現(xiàn)產(chǎn)業(yè)興旺的關鍵驅動力，具有強大的產(chǎn)業(yè)聯(lián)動作用，能夠產(chǎn)生消費力和投資動力，促進區(qū)域經(jīng)濟的發(fā)展，形成以旅游產(chǎn)業(yè)為主導的泛旅游產(chǎn)業(yè)集群，從而帶動人力、物力、財力、資本集聚，促進城鎮(zhèn)化水平的提高[17]。十九大報告將鄉(xiāng)村振興分為產(chǎn)業(yè)興旺、生態(tài)宜居、鄉(xiāng)風文明、治理有效和生活富裕5個層面，旅游城鎮(zhèn)化的發(fā)展使貧困地區(qū)經(jīng)濟、生態(tài)、社會和文化等方面得到改善，能夠推動當?shù)禺a(chǎn)業(yè)結構發(fā)生改變、人居生活環(huán)境得到改善，農(nóng)村基礎設施完善，居民消費水平提高，鄉(xiāng)村氛圍逐漸融洽，居民幸福感得到提升，從而實現(xiàn)鄉(xiāng)村振興的發(fā)展(圖1)。

圖1 旅游城鎮(zhèn)化驅動鄉(xiāng)村振興發(fā)展的作用機制Fig. 1 Mechanism of tourism urbanization driving rural revitalization and development

1.3 旅游城鎮(zhèn)化與鄉(xiāng)村振興協(xié)調(diào)發(fā)展系統(tǒng)因果反饋

系統(tǒng)動力學(System Dynamics，SD)是Forrester在1956年創(chuàng)立的，是用系統(tǒng)、整體的視角，結合定性與定量分析的手段，借助計算機模擬、分析、解決復雜多系統(tǒng)問題的方法。系統(tǒng)動力學能夠根據(jù)系統(tǒng)的結構與功能，探討系統(tǒng)之間的因果反饋關系，尋找該系統(tǒng)問題產(chǎn)生的根源[18]。在旅游城鎮(zhèn)化與鄉(xiāng)村振興的研究領域中，國內(nèi)學者主要將SD仿真應用于旅游城鎮(zhèn)化響應機制分析[16]、旅游城鎮(zhèn)化的研究[19]、工業(yè)與旅游業(yè)推動鄉(xiāng)村振興發(fā)展的機制[20]以及鄉(xiāng)村旅游[21]等方面。

本文主要涉及旅游城鎮(zhèn)化與鄉(xiāng)村振興兩個子系統(tǒng)，張家界作為“旅游強市”，是以旅游產(chǎn)業(yè)帶動城鎮(zhèn)化水平發(fā)展，故旅游城鎮(zhèn)化從人口、經(jīng)濟、旅游3個子系統(tǒng)選取解釋變量，選取旅游總收入、旅游總收入占GDP的比重、旅游人數(shù)、旅游業(yè)投資、城鎮(zhèn)化率等相關變量來構建旅游城鎮(zhèn)化的因果循環(huán)圖[22](圖2)。同時依據(jù)鄉(xiāng)村振興的五個維度以及張家界市鄉(xiāng)村的發(fā)展現(xiàn)狀，從產(chǎn)業(yè)興旺中選取農(nóng)業(yè)生產(chǎn)總值、第一產(chǎn)業(yè)總值、農(nóng)業(yè)機械化水平等指標來表征產(chǎn)業(yè)發(fā)展的水平，從鄉(xiāng)風文明中選取居民受教育水平、農(nóng)村普通在校學生數(shù)等指標來表征農(nóng)村居民的文化素養(yǎng)，從生活富裕中選取農(nóng)村居民人均可支配收入、居民工資性收入來表征居民的收入水平，從生態(tài)宜居中選用化肥使用量、農(nóng)藥使用量以及森林覆蓋率等指標來表征農(nóng)村居民自然環(huán)境宜居狀況，從治理有效中選取村居民委員會的個數(shù)、通自來水村的個數(shù)等指標來表征農(nóng)村治理水平的效果，從而構建鄉(xiāng)村振興的因果循環(huán)圖[23](圖3)。

圖2 旅游城鎮(zhèn)化子系統(tǒng)Fig. 2 Tourism urbanization subsystem

圖3 鄉(xiāng)村振興子系統(tǒng)Fig. 3 Rural revitalization subsystem

1.4 旅游城鎮(zhèn)化與鄉(xiāng)村振興系統(tǒng)結構分析

本文構造的旅游城鎮(zhèn)化與鄉(xiāng)村振興的系統(tǒng)包含產(chǎn)業(yè)子系統(tǒng)、社會文化子系統(tǒng)、生態(tài)環(huán)境子系統(tǒng)和居民幸福子系統(tǒng)4個子系統(tǒng)模塊(圖4)。

圖4 旅游城鎮(zhèn)化驅動鄉(xiāng)村振興發(fā)展的存量流量圖Fig. 4 Stock flow diagram of rural revitalization and development driven by tourism urbanization

(1)產(chǎn)業(yè)發(fā)展子系統(tǒng) 在鄉(xiāng)村振興的發(fā)展下，旅游產(chǎn)業(yè)可以推動旅游資源由城市流向農(nóng)村或者是貧困地區(qū)，是促進農(nóng)村地區(qū)減貧的重要手段，實現(xiàn)人人共享成果、滿足人民美好生活的需要[24]。快速發(fā)展旅游產(chǎn)業(yè)不僅能夠促進區(qū)域經(jīng)濟水平的發(fā)展，加快城鎮(zhèn)化的進程，還能夠提供大量就業(yè)崗位，提高農(nóng)村就業(yè)人口，促進居民人均可支配收入增加，農(nóng)民的生活水平也會隨之提升。因此本文主要選取農(nóng)業(yè)生產(chǎn)總值、旅游固定資產(chǎn)投資額、GDP、游客數(shù)量、旅游總收入、第三產(chǎn)業(yè)從業(yè)人員、A級以上景區(qū)數(shù)量等變量來反映旅游子系統(tǒng)。其中，農(nóng)業(yè)生產(chǎn)總值、旅游人數(shù)、GDP為水平變量，以農(nóng)業(yè)生產(chǎn)總值的增加量、旅游人數(shù)增加量、GDP的增加量為速率變量，旅游產(chǎn)業(yè)投資額、旅游就業(yè)人數(shù)、農(nóng)業(yè)投資、第一產(chǎn)業(yè)產(chǎn)值、旅游總收入、旅游總收入增長指數(shù)等為輔助變量，并與其他子系統(tǒng)連接，來構建旅游產(chǎn)業(yè)子系統(tǒng)。

(2)社會文化子系統(tǒng) 城鎮(zhèn)化是人口非農(nóng)化、經(jīng)濟非農(nóng)化和空間利用非農(nóng)化的區(qū)域變化過程[25]。曹翔認為農(nóng)村人口城鎮(zhèn)化對提高居民消費具有重要推動作用[26]，因此促進農(nóng)村城鎮(zhèn)化的發(fā)展對于提高鄉(xiāng)村振興水平起到一定的作用。一方面鄉(xiāng)村地區(qū)通過發(fā)展旅游產(chǎn)業(yè)，帶動第三產(chǎn)業(yè)GDP的增長，居民人均可支配收入增加，居民消費結構發(fā)生變化，從而使居民受教育水平和醫(yī)療水平逐步提高。另一方面，隨著居民生活水平提高，居民外出旅游意識加強，從而又帶動旅游業(yè)經(jīng)濟的增長。同時隨著城鎮(zhèn)化進程加快，居民生產(chǎn)、生活受到經(jīng)濟快速發(fā)展以及旅游資源、生態(tài)環(huán)境等多方面的影響，鄉(xiāng)村人口外遷數(shù)量增加，人才流失嚴重。因此選取城鎮(zhèn)化率和農(nóng)村居民人均可支配收入作為水平變量，以城鎮(zhèn)化增長率和可支配收入增加量作為速率變量，以農(nóng)村普通高校在讀學生、鄉(xiāng)村人口的數(shù)量、財政支出等作為輔助變量，來構建社會文化子系統(tǒng)。

(3)生態(tài)環(huán)境子系統(tǒng) 生態(tài)環(huán)境對城鎮(zhèn)化具有支撐作用[27]，一方面生態(tài)環(huán)境的承載力決定了旅游目的地的可持續(xù)發(fā)展，另一方面生態(tài)環(huán)境需要消耗旅游城鎮(zhèn)化過程中所產(chǎn)生的各種廢棄物。因此旅游城鎮(zhèn)化的發(fā)展要以生態(tài)環(huán)境的承載力為限，超出限度就會導致居民的生活水平下降，居民幸福感降低[28]，人均耕地面積減少，鄉(xiāng)村振興水平下降，居民人口流失嚴重。因此選取耕地面積、A級景區(qū)為水平變量，以耕地面積增加量、耕地面積減少量、A級景區(qū)增加量為速率變量，以環(huán)保投資支出、化肥使用強度、農(nóng)藥使用強度為輔助變量，并引入環(huán)境污染變化指數(shù)變量綜合反映環(huán)境子系統(tǒng)。

(4)居民幸福子系統(tǒng) 居民幸福指數(shù)是一種衡量居民生活幸福的指數(shù)、它涵蓋了經(jīng)濟發(fā)展指標、收入消費指標、生存環(huán)境指標以及精神狀況指標，能夠很好地反映鄉(xiāng)村振興的發(fā)展狀況。目前關于居民幸福感測算已經(jīng)有了豐富的成果。程國棟、徐中民等人首先提出建立國民幸福生活核算體系的構想[29]。郭玲玲等人在原有公式的基礎上引入了環(huán)境污染指數(shù)，對公式進行了修正[30]。本文結合前人的研究，基于旅游城鎮(zhèn)化的發(fā)展特點，在公式中增加城鎮(zhèn)化的發(fā)展水平指數(shù)，運用客觀指標構建農(nóng)村居民幸福感指數(shù)，客觀反映農(nóng)村居民幸福感受到外界因素的影響。

(1)

式中，H表示居民幸福指數(shù)，W表示旅游收入增長指數(shù)，U為城鎮(zhèn)化水平，G為基尼系數(shù)，J為失業(yè)率，E為環(huán)境污染指數(shù)，S為旅游風險指數(shù)，I為居民消費價格指數(shù)。

1.5 主要變量說明與數(shù)據(jù)收集

系統(tǒng)模擬的空間邊界為張家界市，仿真時間為2005—2035年，仿真步長為1年，以2005年指標數(shù)據(jù)為初始值，以方案1數(shù)據(jù)作為參照基數(shù)，同時以2005—2020年的歷史數(shù)據(jù)為模型參數(shù)的確定依據(jù)?；跀?shù)據(jù)的可獲得性、可衡量性和科學性原則，模型中原始數(shù)據(jù)參數(shù)全部來源于《湖南省統(tǒng)計年鑒》、《張家界市統(tǒng)計年鑒》以及張家界統(tǒng)計公報。根據(jù)系統(tǒng)各變量的特性以及相互之間的關系,主要采取算術平均法、回歸分析、表函數(shù)、趨勢法等進行模型賦值(表1)。

表1 部分變量關系式Table 1 Partial variable relation

2 仿真結果與討論

2.1 模型的歷史性檢驗

在對張家界市旅游城鎮(zhèn)化驅動鄉(xiāng)村系統(tǒng)的機制進行建模的過程中，結合張家界市旅游業(yè)本身的特征，并參考旅游發(fā)展及其他產(chǎn)業(yè)發(fā)展的系統(tǒng)動力學模型，擬合了實際系統(tǒng)的主要特點，所選數(shù)據(jù)真實且具有現(xiàn)實意義，系統(tǒng)模型與實際結構相吻合。為確保模型的科學性、有效性和精確性，從系統(tǒng)中選取農(nóng)業(yè)生產(chǎn)總值、GDP、農(nóng)村居民人均可支配收入、旅游總收入4個具有代表性的變量進行歷史性檢驗，時間為2005—2020年，由于數(shù)據(jù)量較大，從中截取2015—2020年的數(shù)據(jù)進行檢驗。相關研究認為檢驗誤差在-10%～15%變動均可以接受，檢驗結果表明：誤差介于-0.000 89～0.004 3，模型的精確度較高，可用于預測發(fā)展和政策模擬(表2)。

表2 系統(tǒng)動力模型的檢驗結果Table 2 Test results of system dynamic model

2.2 參數(shù)設置

為了探究不同政策下旅游城鎮(zhèn)化與鄉(xiāng)村振興的協(xié)調(diào)發(fā)展變化，選取旅游業(yè)固定資產(chǎn)投資系數(shù)、農(nóng)業(yè)固定資產(chǎn)投資系數(shù)、第三產(chǎn)業(yè)產(chǎn)值比重、環(huán)保投資占比、城鎮(zhèn)化水平增長率等指標作為調(diào)控參數(shù)，參數(shù)的調(diào)整根據(jù)張家界市2005—2020年的發(fā)展趨勢和張家界市“十四五”規(guī)劃，對未來10年的調(diào)控策略進行多情景仿真研究。需要說明的是，旅游業(yè)、農(nóng)業(yè)投資系數(shù)越大，證明政府對旅游業(yè)以及農(nóng)業(yè)的支持程度越高；第三產(chǎn)業(yè)產(chǎn)值比重反映了旅游產(chǎn)業(yè)產(chǎn)值占國內(nèi)生產(chǎn)總值的比重，側面反映了旅游業(yè)帶來的經(jīng)濟增量；環(huán)境保護力度越大，越有利于改善鄉(xiāng)村的生態(tài)環(huán)境，從而促進鄉(xiāng)村振興的發(fā)展；城鎮(zhèn)化水平增長率表明了旅游業(yè)促進城鎮(zhèn)化水平的發(fā)展程度。通過分析比較各模擬方案的現(xiàn)狀與趨勢，得出張家界市各方案的優(yōu)勢和劣勢，通過考察旅游總收入、農(nóng)村人均可支配收入、耕地面積、環(huán)境污染指數(shù)、居民幸福指數(shù)等表量的動態(tài)變化，對旅游城鎮(zhèn)化與鄉(xiāng)村振興進行仿真模擬研究(表3)。

表3 系統(tǒng)動力模型的情景設置及調(diào)控參數(shù)變化Table 3 Scenario setting and regulation parameter change of system dynamic model

2.3 仿真結果

2.3.1 模擬方案比較

方案一：自然增長模式

保持系統(tǒng)參數(shù)不變，按照現(xiàn)有參數(shù)進行模擬(圖5)：(1)經(jīng)濟增長明顯，農(nóng)村居民人均可支配收入呈現(xiàn)明顯上升趨勢，由2020年的11 538元/人增長到2035年的31 453元/人，旅游總收入在2020—2030年呈現(xiàn)迅速增長的趨勢，而后在2030—2035年增長緩慢，旅游總收入到2035年達到1 857億元。由于產(chǎn)業(yè)結構未發(fā)生改變，農(nóng)業(yè)生產(chǎn)總值呈現(xiàn)緩慢上升，到2035年增長至90.2億元。(2)環(huán)境污染指數(shù)在2020—2035年呈現(xiàn)上升趨勢，隨著疫情防控的有效開展，人們出現(xiàn)報復性的旅游消費，導致游客數(shù)量增多，環(huán)境污染加重。而且按照現(xiàn)有的耕地保護措施，人均耕地面積在2020—2035年將呈現(xiàn)持續(xù)性下降趨勢，到2035年人均耕地面積與2020年相比減少了0.532 hm2。(3)居民幸福指數(shù)在2020—2035年緩慢增長，由此可以看出居民幸福指數(shù)是由旅游收入增長指數(shù)、環(huán)境污染、居民消費水平、就業(yè)人數(shù)等多種因素共同作用的結果。

綜上所述，在自然增長模式下，旅游業(yè)的發(fā)展雖然能夠帶來一定的經(jīng)濟增長，但是隨著城鎮(zhèn)化進程的增加會對環(huán)境污染造成明顯的影響。

方案二：高速增長模式

與自然增長模式相比，高速發(fā)展模式下，城鎮(zhèn)化發(fā)展速度將更快，從而使鄉(xiāng)村經(jīng)濟發(fā)展力度更大(圖5)。(1)經(jīng)濟迅速增長，農(nóng)村居民可支配收入、旅游總收入在4個方案中發(fā)展最快、數(shù)值也最高。農(nóng)村居民可支配收入和旅游總收入在2035年分別達到33 234元/人和4 123億元。農(nóng)業(yè)生產(chǎn)總值開始緩慢上升，到2035年增加至100億元。(2)環(huán)境污染指數(shù)從2020—2035年迅速上升，是4個方案中上升最快的。由此可以看出伴隨旅游城鎮(zhèn)化迅速發(fā)展，農(nóng)村人才流失，科技水平下降，農(nóng)藥和化肥使用量增多，并且超出了環(huán)境治理能力。人均耕地面積在4個方案中下降最快，到2035年人均耕地面積僅為0.23畝/人。(3)居民幸福指數(shù)在2020—2035年之間開始急劇下降，說明片面發(fā)展旅游城鎮(zhèn)化，不利于居民幸福感的提升，同時也會使居民參與旅游的意識下降。

綜上所述，高速增長模式只是單純追求經(jīng)濟的迅速增長，雖然能短暫地提升居民的收入水平，但與此而來的是環(huán)境的嚴重污染，居民的外遷率提高，鄉(xiāng)村人口減少，不利于鄉(xiāng)村振興的發(fā)展。

方案三：低速增長模式

此方案是一種相對保守的發(fā)展模式，重在對農(nóng)村環(huán)境進行合理的保護，使農(nóng)村居民擁有一個宜居的生活環(huán)境(圖5)。(1)經(jīng)濟增長速度緩慢，旅游總收入和農(nóng)業(yè)生產(chǎn)總值在四個方案中是最低的。受旅游固定資產(chǎn)投資比重和農(nóng)業(yè)投資率的減少，經(jīng)濟增長被抑制。到2035年，旅游總收入和農(nóng)村居民人均可支配收入311億元和9 604元/人，農(nóng)業(yè)生產(chǎn)總值為77億元，位列四種發(fā)展模式之尾。(2)環(huán)境污染指數(shù)，隨著環(huán)境保護投資系數(shù)的增加，環(huán)境污染指數(shù)在2020—2035年呈現(xiàn)下降趨勢，人均耕地面積開始呈現(xiàn)上升趨勢，與高速發(fā)展和自然增長模式相比，耕地利用呈現(xiàn)出良性發(fā)展的趨勢。(3)居民幸福指數(shù)隨著環(huán)境污染指數(shù)的下降呈現(xiàn)出先上升(2020—2030年)后下降(2030—2035年)的趨勢，超越其他3種模式，證明居民的幸福指數(shù)與生活環(huán)境息息相關。

綜上所述，低速增長模式雖然使生態(tài)環(huán)境得到保護，但是經(jīng)濟增長指數(shù)卻逐漸下降，居民的收入水平降低，在一定程度上抑制了當?shù)芈糜纬擎?zhèn)化的發(fā)展。

方案四：協(xié)調(diào)發(fā)展模式

針對上述3種發(fā)展模式存在的弊端，提出了基于產(chǎn)業(yè)發(fā)展-生態(tài)環(huán)境-社會文化-居民幸福協(xié)調(diào)發(fā)展的模式，并進行參數(shù)調(diào)整(圖5)。(1)經(jīng)濟增長較快，并且呈現(xiàn)上升趨勢。到2035年，農(nóng)村居民人均可支配收入和旅游總收入32 312元/人和3 924億元，僅次于高速發(fā)展模式。農(nóng)業(yè)生產(chǎn)總值是4個方案里最高的，到2035年達到153億元。(2)與其他3種模式不同的是，協(xié)調(diào)發(fā)展模式在發(fā)展經(jīng)濟的同時扭轉了環(huán)境污染指數(shù)，環(huán)境污染指數(shù)的下降僅次于低速發(fā)展模式。在此模式下，耕地面積呈現(xiàn)緩慢減少的趨勢，人均耕地面積開始迅速上升。(3)居民幸福指數(shù)在2020年開始明顯上升，僅次于低速發(fā)展模式。

圖5 不同情景下旅游城鎮(zhèn)化驅動鄉(xiāng)村振興的仿真結果分析Fig. 5 Simulation results of rural revitalization driven by tourism urbanization under different scenarios

綜上所述，協(xié)調(diào)發(fā)展模式既能使經(jīng)濟得到很好的發(fā)展，同時也能避免環(huán)境污染，提高居民幸福指數(shù)，從而促進鄉(xiāng)村振興的可持續(xù)發(fā)展。

2.3.2 不同方案下的鄉(xiāng)村振興發(fā)展水平 由于鄉(xiāng)村振興是一個龐大且復雜的系統(tǒng)，涉及產(chǎn)業(yè)、經(jīng)濟、人口、社會、環(huán)境等領域，影響因素頗多[31]。因此，本文通過調(diào)整旅游城鎮(zhèn)化發(fā)展水平的參數(shù)來模擬不同情景下鄉(xiāng)村振興發(fā)展水平。由圖6可知，不同方案下的鄉(xiāng)村振興發(fā)展水平存在明顯的差異。自然增長模式、低速增長模式和協(xié)調(diào)增長模式下，鄉(xiāng)村振興發(fā)展指數(shù)呈現(xiàn)上升趨勢，其中協(xié)調(diào)發(fā)展模式下，鄉(xiāng)村振興發(fā)展指數(shù)上升最快，表明鄉(xiāng)村振興的發(fā)展要綜合考慮產(chǎn)業(yè)興旺、生態(tài)宜居、生活富裕、鄉(xiāng)風文明和治理有效這5個層面，缺一不可。在高速增長模式下，由于片面追求經(jīng)濟的增長，鄉(xiāng)村振興指數(shù)呈現(xiàn)先上升(2020—2025年)再下降(2025—2035年)的趨勢。自2017年以來我國開始大力發(fā)展鄉(xiāng)村振興，使鄉(xiāng)村經(jīng)濟開始得到發(fā)展。而后隨著疫情的爆發(fā)，旅游業(yè)發(fā)展停滯不前。在后疫情時代(2020—2030年)，旅游業(yè)開始緩慢恢復，促使張家界的鄉(xiāng)村振興指數(shù)開始緩慢上升，但隨之而來帶來的是生態(tài)環(huán)境遭到破壞，農(nóng)村居民外遷、農(nóng)村“空心化”等問題，故2030年之后鄉(xiāng)村振興發(fā)展指數(shù)開始迅速下跌。在協(xié)調(diào)發(fā)展模式下，由于綜合考慮產(chǎn)業(yè)經(jīng)濟、生態(tài)環(huán)境、居民幸福等因素，鄉(xiāng)村振興指數(shù)在2035年達到0.768，上升至良好的水平狀態(tài)，說明協(xié)調(diào)發(fā)展模式下可以更好提升張家界鄉(xiāng)村振興。

圖6 不同情景下鄉(xiāng)村振興發(fā)展指數(shù)Fig. 6 Rural revitalization and development index under different scenarios

3 研究結論與政策建議

3.1 研究結論

本文運用系統(tǒng)動力學的方法，對張家界市旅游城鎮(zhèn)化與鄉(xiāng)村振興協(xié)調(diào)發(fā)展進行研究。首先構建旅游城鎮(zhèn)化與鄉(xiāng)村振興的SD模型，包括產(chǎn)業(yè)經(jīng)濟子系統(tǒng)、社會文化子系統(tǒng)、生態(tài)環(huán)境子系統(tǒng)以及居民幸福4個子系統(tǒng)。其次，設計了自然發(fā)展、高速發(fā)展、低速發(fā)展以及協(xié)調(diào)發(fā)展的4種不同方案，最后進行仿真模擬并且對4種方案進行了比較分析。結果表明：

(1)本文構建的旅游城鎮(zhèn)化與鄉(xiāng)村振興的系統(tǒng)動力學模型具有很好的模擬性和穩(wěn)定性，說明該模型能夠反映真實的情況。而且從系統(tǒng)的角度剖析旅游城鎮(zhèn)化與鄉(xiāng)村振興的演化規(guī)律，分析其發(fā)展特點和趨勢，可以在一定程度上避免做出片面的決策。通過仿真模擬發(fā)現(xiàn)，張家界地區(qū)各個子系統(tǒng)在仿真模擬期間總體保持了增長態(tài)勢，但增長速度有所差異，在不同子系統(tǒng)中不同因素的發(fā)展趨勢也存在差異。對比4種模擬方案發(fā)現(xiàn)，方案四的仿真效果最優(yōu)，此方案兼顧產(chǎn)業(yè)發(fā)展-社會文化-生態(tài)環(huán)境-居民幸福4個子系統(tǒng)之間的協(xié)同發(fā)展，為張家界市鄉(xiāng)村振興的發(fā)展提供了最優(yōu)的發(fā)展模式。

(2)通過對張家界市鄉(xiāng)村振興與旅游城鎮(zhèn)化的協(xié)調(diào)發(fā)展仿真研究，認為張家界市旅游城鎮(zhèn)化與鄉(xiāng)村振興協(xié)調(diào)發(fā)展應綜合考慮產(chǎn)業(yè)發(fā)展、生態(tài)環(huán)境、社會文化以及居民幸福等要素，缺一不可。由于旅游城鎮(zhèn)化的發(fā)展對鄉(xiāng)村振興的影響表現(xiàn)為旅游城鎮(zhèn)化發(fā)展所帶來的多種外部因素，既有居民收入水平、教育水平、醫(yī)療水平的提高，同時也會帶來環(huán)境污染、貧富差距增大、公共安全風險等問題。而且在系統(tǒng)結構上，旅游城鎮(zhèn)化與鄉(xiāng)村振興的協(xié)調(diào)發(fā)展由產(chǎn)業(yè)經(jīng)濟、社會文化、生態(tài)環(huán)境、以及居民幸福等各個系統(tǒng)通過一個或多個變量相互作用和關聯(lián)，形成具有特定功能的有機整體，因此制定相關政策時應綜合考慮旅游城鎮(zhèn)化與鄉(xiāng)村振興內(nèi)部和外部各種因素，才能實現(xiàn)鄉(xiāng)村的可持續(xù)發(fā)展。

3.2 政策建議

根據(jù)研究結果，本文基于協(xié)調(diào)發(fā)展模式，優(yōu)化旅游城鎮(zhèn)化的發(fā)展結構，提升農(nóng)村居民幸福感，從而促進鄉(xiāng)村振興的可持續(xù)發(fā)展。

(1)政府應加大環(huán)境保護治理支出占比，提高環(huán)境治理能力，正如習近平總書記所說：“地球是全人類賴以生存的唯一家園”。人與自然是命運共同體，我們要在發(fā)展旅游產(chǎn)業(yè)的過程中保護生態(tài)環(huán)境，在生態(tài)環(huán)境保護中發(fā)展旅游產(chǎn)業(yè)，才能使居民真正獲得幸福感。尤其隨著城鎮(zhèn)化水平的提高，我們更要在注重經(jīng)濟效應的同時關注鄉(xiāng)村環(huán)境與生態(tài)的發(fā)展情況，建設生態(tài)宜居、生活富裕、鄉(xiāng)風文明的美麗鄉(xiāng)村。

(2)由于旅游行業(yè)本身具有脆弱性和敏感性，張家界市過度依靠旅游業(yè)的發(fā)展，抗風險能力較弱，2020年1月因疫情防控張家界暫時關閉了所有的景區(qū)，國內(nèi)游客接待數(shù)量由2019年的7 912.3萬人次下跌到2020年的4 949.2萬人次，旅游總收入由905.6億元下跌至568.95億元，給張家界經(jīng)濟發(fā)展帶來了巨大的損失。因此政府要促進旅游產(chǎn)業(yè)轉型發(fā)展，注重產(chǎn)業(yè)結構的轉型升級，使農(nóng)業(yè)與一二三產(chǎn)業(yè)融合發(fā)展，推動“工業(yè)+旅游”深度融合，堅持城鄉(xiāng)融合發(fā)展，發(fā)揮自身優(yōu)勢，將紅色、民族、民俗的特色旅游文化與當?shù)氐牡匦蔚孛步Y合起來，創(chuàng)造文旅融合的特色旅游產(chǎn)品，同時完善旅游公共衛(wèi)生安全機制，提高旅游衛(wèi)生安全應急能力，促進旅游產(chǎn)業(yè)的可持續(xù)發(fā)展和高質(zhì)量發(fā)展。

此外，本研究雖具有理論意義和實踐價值，但也存在一定的局限性。一方面鄉(xiāng)村振興和旅游城鎮(zhèn)化是兩個錯綜復雜的系統(tǒng)，指標體系繁雜，指標選取具有不確定性，比如鄉(xiāng)村振興的五個維度的所有指標并不能完全包含在內(nèi)，只是選取了部分有代表性的指標，而且模型主要依靠產(chǎn)業(yè)發(fā)展、生態(tài)環(huán)境、社會文化、居民幸福這4個系統(tǒng)來考慮鄉(xiāng)村振興與旅游城鎮(zhèn)化的協(xié)調(diào)發(fā)展，實際上并不只限于這幾種因素。另一方面，旅游城鎮(zhèn)化驅動鄉(xiāng)村振興的發(fā)展是一個時變系統(tǒng)，同一地區(qū)處于不同發(fā)展階段的驅動機制系統(tǒng)的演變可能會發(fā)生變化，會受到各種不同政策和突發(fā)情況的影響，本文在研究中未考慮突發(fā)性疫情對旅游業(yè)的影響，因此不能檢驗該系統(tǒng)動力模型的普適性和科學性，需要持續(xù)關注旅游城鎮(zhèn)化與鄉(xiāng)村振興協(xié)調(diào)發(fā)展的動態(tài)趨勢，這也是未來研究中要考慮的因素。