黃河流域生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)與文旅融合耦合協(xié)調(diào)發(fā)展研究

張嘉貝，王珂珂，施天樂(lè)，2，劉潤(rùn)歌

（1.鄭州大學(xué) 管理學(xué)院，河南 鄭州 450001；2.鄭州大學(xué) 黃河生態(tài)保護(hù)與區(qū)域協(xié)調(diào)發(fā)展研究院，河南 鄭州 450001；3.黃河科技學(xué)院 藝體學(xué)部，河南 鄭州 450006）

0 引言

黃河流域是中華文明的重要發(fā)祥地，孕育了河湟文化、關(guān)中文化、齊魯文化等瑰麗的黃河文化，分布著西安、洛陽(yáng)、開(kāi)封等多朝古都，為文旅融合發(fā)展提供了得天獨(dú)厚的現(xiàn)實(shí)條件。黃河流域也是中國(guó)重要的生態(tài)安全屏障，是中國(guó)人口活動(dòng)和經(jīng)濟(jì)發(fā)展的重要區(qū)域［1］。近年來(lái)，黃河流域生態(tài)環(huán)境脆弱、水土流失嚴(yán)重、水旱災(zāi)害頻發(fā)等生態(tài)問(wèn)題仍舊突出［2］，而作為文旅融合發(fā)展的自然資源基礎(chǔ)，生態(tài)保護(hù)一定是黃河流域高質(zhì)量發(fā)展的重要前提。如何推動(dòng)黃河流域生態(tài)及文旅協(xié)同發(fā)展現(xiàn)已成為社會(huì)各界關(guān)注的一項(xiàng)重大議題。依托寬帶中國(guó)，5G及工業(yè)互聯(lián)網(wǎng)的數(shù)字經(jīng)濟(jì)［3］已成為破解生態(tài)及文旅發(fā)展問(wèn)題的關(guān)鍵所在［4］。2019 年，黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展上升為重大國(guó)家戰(zhàn)略。此后，《黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展規(guī)劃綱要》對(duì)黃河流域生態(tài)優(yōu)先發(fā)展、數(shù)字經(jīng)濟(jì)發(fā)展、文旅融合發(fā)展等作出了戰(zhàn)略部署。因此深入探究黃河流域生態(tài)環(huán)境保護(hù)、數(shù)字經(jīng)濟(jì)及文旅融合之間的內(nèi)在關(guān)聯(lián)對(duì)黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展具有重要的現(xiàn)實(shí)意義。

學(xué)界對(duì)于生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)和文旅融合的研究成果較為豐富，但相關(guān)研究多聚焦于兩兩之間的關(guān)聯(lián)［5，6］。近年來(lái)，學(xué)界基于不同視角展開(kāi)了對(duì)數(shù)字經(jīng)濟(jì)和生態(tài)環(huán)境相關(guān)關(guān)系的研究。張彥博從生態(tài)規(guī)制的視角分析認(rèn)為數(shù)字經(jīng)濟(jì)的發(fā)展與互聯(lián)網(wǎng)的廣泛應(yīng)用為公眾提供了更多的參與機(jī)會(huì)，進(jìn)而極大改善了環(huán)境保護(hù)的現(xiàn)狀，并利于改變政府對(duì)企業(yè)污染排放的管控模式［7］；郭峰從污染減排的角度研究，表明互聯(lián)網(wǎng)的信息傳播能力和網(wǎng)絡(luò)經(jīng)濟(jì)的融合，為城市工業(yè)技術(shù)的發(fā)展帶來(lái)了巨大的潛力，提高了要素的產(chǎn)出效率，此外，利用數(shù)字技術(shù)，不僅可以有效降低城市地區(qū)的污染物排放量，而且可以突破產(chǎn)業(yè)數(shù)字融合的制度障礙，進(jìn)一步提高減排的效果［8］；彭代彥從全社會(huì)的角度分析，數(shù)字經(jīng)濟(jì)的發(fā)展不僅能改善公眾的生態(tài)意識(shí)，還可以深刻地影響其行為，進(jìn)而實(shí)現(xiàn)可持續(xù)發(fā)展目標(biāo)［9］。

伴隨著數(shù)字時(shí)代的到來(lái)，數(shù)字經(jīng)濟(jì)已成為文化旅游產(chǎn)業(yè)的一個(gè)重要組成部分，對(duì)文旅產(chǎn)業(yè)之間的融合深度［10］、融合路徑［11－13］、融合模式［14］等帶來(lái)重要變革。陳宇等以黑龍江文旅產(chǎn)業(yè)融合為背景，在數(shù)字經(jīng)濟(jì)的作用下提出平臺(tái)化推進(jìn)黑龍江文旅產(chǎn)業(yè)深度融合的建議［15］；陳斌等以四川南江縣文旅游資源普查與開(kāi)發(fā)為背景探索數(shù)字經(jīng)濟(jì)背景下文旅融合新路徑［16］；王雄青、張祝平、黃先開(kāi)等學(xué)者分別就紅色文化旅游、鄉(xiāng)村旅游、數(shù)字經(jīng)濟(jì)下的城市文化旅游以及新時(shí)代背景下的文旅融合高質(zhì)量發(fā)展提出了具體的策略和實(shí)施路徑，為推動(dòng)文旅產(chǎn)業(yè)發(fā)展提供了有益的建議［17－19］；李鳳亮、李麗云等學(xué)者深入探討了以低碳為基礎(chǔ)的文化旅游融合發(fā)展的可行性、可持續(xù)性以及可行的實(shí)施方案［20，21］。

文旅是促進(jìn)生態(tài)文明建設(shè)的重要力量。作為新興產(chǎn)業(yè)，文旅產(chǎn)業(yè)具有能耗低、污染少、排放少的特點(diǎn)，是綠色生態(tài)和環(huán)境友好型產(chǎn)業(yè)，這正契合了生態(tài)文明的綠色、環(huán)保、循環(huán)、永續(xù)發(fā)展理念［22］。多數(shù)學(xué)者對(duì)于兩者耦合路徑、協(xié)同發(fā)展以及體系構(gòu)建進(jìn)行了深入研究，其中蘇永波基于文旅與生態(tài)文明建設(shè)相互影響、相互交叉、相互滲透的天然的耦合性，從系統(tǒng)的角度，借助于主輔嵌入式耦合系統(tǒng)理論構(gòu)建耦合矩陣，分析耦合點(diǎn)探索耦合路徑［23］。楊曉東等構(gòu)建黃河流域文旅與生態(tài)保護(hù)同都市圈協(xié)同發(fā)展聯(lián)合體，進(jìn)而探索區(qū)域文旅合作發(fā)展新模式［24］。魏妮茜等以產(chǎn)業(yè)、產(chǎn)品、科技和服務(wù)為基礎(chǔ)，構(gòu)筑出一個(gè)具有高質(zhì)量發(fā)展的生態(tài)系統(tǒng)，深入研究了文化旅游與生態(tài)環(huán)境之間的關(guān)聯(lián)性及其實(shí)現(xiàn)途徑［25］。

當(dāng)前學(xué)界關(guān)于生態(tài)保護(hù)、文旅融合與經(jīng)濟(jì)增長(zhǎng)的研究多聚焦于單一要素或三者之間的兩種，如探究文旅融合發(fā)展的社會(huì)效益、經(jīng)濟(jì)可持續(xù)發(fā)展能力的測(cè)算、生態(tài)環(huán)境與旅游業(yè)發(fā)展的沖突，生態(tài)保護(hù)與經(jīng)濟(jì)發(fā)展的相互作用等，而將生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)與文旅融合三者納入統(tǒng)一研究框架，探究其耦合協(xié)調(diào)水平的實(shí)證研究鮮少。并且相關(guān)研究多從國(guó)家、省域及單一城市層面展開(kāi)，研究范圍有待進(jìn)一步細(xì)化，且少有基于重點(diǎn)流域進(jìn)行探索，忽略了流域地理單元生態(tài)保護(hù)－文旅融合－經(jīng)濟(jì)增長(zhǎng)的關(guān)聯(lián)性。因此，本文以黃河流域57 個(gè)地級(jí)市2012－2021 年面板數(shù)據(jù)為基礎(chǔ)，測(cè)算生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)、文旅融合（ecological protection，digital economy，and integration of culture and tourism，EDC）三系統(tǒng)的耦合協(xié)調(diào)度，探究三系統(tǒng)之間的空間相關(guān)性，并通過(guò)障礙因子模型進(jìn)一步探究，以期為黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展提供科學(xué)的理論參考。

1 研究區(qū)域、研究方法及數(shù)據(jù)來(lái)源

1.1 研究區(qū)域

黃河起源于青藏高原巴顏喀拉山脈北部，流經(jīng)9 個(gè)省份，最終匯入渤海，全長(zhǎng)5 464km，被譽(yù)為“母親河”，是中國(guó)重要的能源資源富集區(qū)、生態(tài)保護(hù)屏障和經(jīng)濟(jì)發(fā)展區(qū)，具有重要的經(jīng)濟(jì)和社會(huì)意義。為了確保研究范圍的完整性，本文以自然流域?yàn)榛A(chǔ)，結(jié)合地域研究其單元的完整性和黃河流域經(jīng)濟(jì)發(fā)展的直接關(guān)聯(lián)性，并對(duì)黃河流域進(jìn)行了詳細(xì)的劃分。本文以黃河流域的57 個(gè)地級(jí)市為研究對(duì)象，其中，上游地區(qū)包括包頭、銀川、中衛(wèi)等17 個(gè)地級(jí)市，中游地區(qū)包括臨汾、太原、銅川等25 個(gè)地級(jí)市，下游地區(qū)包括鄭州、濟(jì)南、德州等15 個(gè)地級(jí)市（圖2）。

圖1 研究區(qū)域Figure 1 Study area

圖2 2012—2021 年黃河流域生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)和文旅融合發(fā)展水平均值Figure 2 The integrated development mean value of ecological protection，digital economy，and integration of cultural tourism in the Yellow River Basin，2012－2021

1.2 研究方法

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

在耦合度模型的基礎(chǔ)上建立EDC 系統(tǒng)的耦合度模型，具體步驟及計(jì)算公式如下：第一步，利用EDC系統(tǒng)的耦合度模型計(jì)算耦合度Z：

式中：Z 為耦合度；E 為生態(tài)保護(hù)水平；D 為數(shù)字經(jīng)濟(jì)發(fā)展水平；C為文旅融合水平。

第二步，為更好體現(xiàn)生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)與文旅融合三系統(tǒng)間的協(xié)調(diào)發(fā)展程度，改進(jìn)耦合度模型建立EDC系統(tǒng)的耦合協(xié)調(diào)度模型，計(jì)算公式如下：

式中：C為各市間三者的協(xié)調(diào)指數(shù)；T 為耦合協(xié)調(diào)度；α、β、ε分別為生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)、文旅融合的系統(tǒng)權(quán)重。參考已有的成果［29］及數(shù)據(jù)分布情況，本文認(rèn)為生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)、文旅融合同等重要，故令α ＝β ＝ε ＝1／3，并將耦合度及耦合協(xié)調(diào)度均分為10 個(gè)階段（表2）。

1.2.2 空間相關(guān)性分析

莫蘭指數(shù)（Moran′s I）是度量相關(guān)變量空間相關(guān)性的一個(gè)重要指標(biāo)［30，31］，本文通過(guò)測(cè)算全局莫蘭指數(shù)（Global Moran′s I）與局部莫蘭指數(shù)（Local Moran′s I）來(lái)探究黃河流域EDC 的空間相關(guān)性及空間集聚特征。全局莫蘭指數(shù)的計(jì)算公式如下：

式中：n為城市數(shù)；S2為樣本方差；Wij為空間權(quán)重；xi和xj為i 城市和j 城市的耦合協(xié)調(diào)度為均值。I∈［－1，1］，I ＞0 表示正向自相關(guān)，即表現(xiàn)為高值與高值相鄰、低值與低值相鄰；I ＜0 表示負(fù)向自相關(guān)，即表現(xiàn)為高值與低值相鄰；I ＝0 則表示耦合協(xié)調(diào)度在空間上隨機(jī)分布。

雖然全局莫蘭指數(shù)能夠揭示黃河流域EDC 系統(tǒng)的耦合協(xié)調(diào)度是否存在空間自相關(guān)性，但是并不能揭示各地市耦合協(xié)調(diào)度的空間集聚方式。因此，利用局部莫蘭指數(shù)進(jìn)行進(jìn)一步分析，計(jì)算公式如下：

式中：S2、Wij和的含義同式（4）。若Ii＞0，表示局部區(qū)域3 個(gè)系統(tǒng)耦合協(xié)調(diào)度相似，集聚模式表現(xiàn)為“高高”集聚或“低低”集聚；若Ii＜0，表示局部區(qū)域3 個(gè)系統(tǒng)耦合協(xié)調(diào)度不相似，集聚模式表現(xiàn)為“低高”集聚或“高低”集聚；若Ii＝0，表示局部區(qū)域耦合協(xié)調(diào)度隨機(jī)分布。

1.2.3 趨勢(shì)分析

利用趨勢(shì)分析技術(shù)，結(jié)合優(yōu)秀的多項(xiàng)式擬合，將區(qū)域散點(diǎn)內(nèi)插，形成一個(gè)完整的趨勢(shì)面，這樣就能夠有效地半定量地研靠大規(guī)模的地理數(shù)據(jù)，并且能夠更加精確地捕捉地理要素的空間分布特征和變化趨勢(shì)［32］。通過(guò)ArcGIS軟件中的趨勢(shì)分析探究黃河流域EDC系統(tǒng)的空間分異特征，計(jì)算公式如下：

1.2.4 Kriging插值法

基于變異函數(shù)理論和結(jié)構(gòu)分析的Kriging 插值法是一種有效的無(wú)偏最優(yōu)估計(jì)方法，可以有效地估計(jì)目標(biāo)區(qū)域內(nèi)的變量，從而達(dá)到最優(yōu)的結(jié)果［33－35］。對(duì)黃河流域地級(jí)市耦合協(xié)調(diào)度進(jìn)行Kriging 插值分析，可以直觀反映各市三要素耦合協(xié)調(diào)程度及變化趨勢(shì)，計(jì)算公式如下：

1.2.5 障礙度模型

協(xié)調(diào)發(fā)展是各子系統(tǒng)及相關(guān)要素相互適應(yīng)、相互依存的發(fā)展過(guò)程，這個(gè)過(guò)程表現(xiàn)為從低級(jí)到高級(jí)、從簡(jiǎn)單到復(fù)雜、從混亂到有序的動(dòng)態(tài)演變過(guò)程。因此，耦合協(xié)調(diào)的水平與各子系統(tǒng)及其相關(guān)要素發(fā)展水平密切相關(guān)。本文采用障礙度模型［36］探究影響黃河流域EDC系統(tǒng)耦合協(xié)調(diào)的各系統(tǒng)的主要障礙因素，計(jì)算公式如下：

式中：Dij為指標(biāo)偏離度，即單項(xiàng)指標(biāo)與耦合協(xié)調(diào)目標(biāo)之間的差距；oij和Oij分別為指標(biāo)層和要素層的障礙度，它們的值越大，說(shuō)明它們對(duì)黃河流域EDC系統(tǒng)耦合協(xié)調(diào)的阻礙作用越大。

1.3 指標(biāo)體系構(gòu)建及數(shù)據(jù)來(lái)源

本文緊扣黃河流域生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)與文旅融合協(xié)調(diào)發(fā)展的作用機(jī)制，構(gòu)建黃河流域生態(tài)—數(shù)字經(jīng)濟(jì)—文旅系統(tǒng)耦合協(xié)調(diào)評(píng)估指標(biāo)體系（表1）。生態(tài)保護(hù)評(píng)估指標(biāo)體系的構(gòu)建遵循科學(xué)性、代表性和可操作性原則，具體篩選依據(jù)主要是生物多樣性、生態(tài)系統(tǒng)組成結(jié)構(gòu)、生態(tài)系統(tǒng)服務(wù)功能以及生態(tài)系統(tǒng)的主要威脅。數(shù)字經(jīng)濟(jì)評(píng)估指標(biāo)體系的構(gòu)建遵循體系科學(xué)、指標(biāo)適用、設(shè)計(jì)合理原則，具體篩選依據(jù)主要是高質(zhì)量發(fā)展的要求、數(shù)字經(jīng)濟(jì)發(fā)展規(guī)律、地區(qū)發(fā)展的實(shí)際以及數(shù)據(jù)來(lái)源的可行性。通過(guò)建立文旅評(píng)估指標(biāo)體系，可以準(zhǔn)確地衡量黃河流域文旅融合發(fā)展水平，更好地了解兩大產(chǎn)業(yè)之間的層次結(jié)構(gòu)關(guān)系，該指標(biāo)體系綜合考慮文化產(chǎn)業(yè)和旅游產(chǎn)業(yè)的發(fā)展情況，并根據(jù)實(shí)際情況篩選出最佳的評(píng)價(jià)指標(biāo)，從而更加準(zhǔn)確地反映兩大產(chǎn)業(yè)的融合發(fā)展水平。

表1 黃河流域EDC系統(tǒng)耦合協(xié)調(diào)評(píng)估指標(biāo)體系Table 1 Evaluation index system for coupling coordination of EDC system in the Yellow River Basin

表2 耦合協(xié)調(diào)度等級(jí)劃分標(biāo)準(zhǔn)__Table 2 Classification criteria for coupling coordination degree levels

遵循以上指標(biāo)選取依據(jù)，生態(tài)保護(hù)指標(biāo)選取主要借鑒朱婧文的研究成果［26］，從生態(tài)基礎(chǔ)、生態(tài)壓力與環(huán)境污染3 個(gè)維度選取指標(biāo)反映黃河流域生態(tài)保護(hù)發(fā)展情況。數(shù)字經(jīng)濟(jì)增長(zhǎng)指標(biāo)選取主要借鑒李蕾的成果［27］，從數(shù)字化基礎(chǔ)、數(shù)字化產(chǎn)業(yè)及數(shù)字普惠金融3 個(gè)方面選取指標(biāo)衡量黃河流域經(jīng)濟(jì)增長(zhǎng)情況。文旅融合指標(biāo)選取主要參考唐承財(cái)?shù)难芯砍晒?8］，從文旅資源、文旅效益、文旅支撐3 個(gè)方面選取指標(biāo)反映黃河流域文旅融合發(fā)展水平。最終構(gòu)建包括生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)、文旅融合3 個(gè)子系統(tǒng)層9個(gè)要素層24 個(gè)指標(biāo)層的黃河流域EDC系統(tǒng)耦合協(xié)調(diào)評(píng)估指標(biāo)體系，指標(biāo)層權(quán)重由熵權(quán)法計(jì)算所得，要素層權(quán)重由對(duì)應(yīng)指標(biāo)層權(quán)重加總所得。

研究樣本中各指標(biāo)數(shù)據(jù)主要來(lái)源于2012—2021 年各市（州）統(tǒng)計(jì)年鑒、各市（州）國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)、各市（州）生態(tài)環(huán)境局發(fā)布的《環(huán)境狀況公報(bào)》、《中國(guó)能源統(tǒng)計(jì)年鑒》、《中國(guó)城市統(tǒng)計(jì)年鑒》以及各省份統(tǒng)計(jì)年鑒等相關(guān)文件。部分?jǐn)?shù)據(jù)缺失采用均值替換法和插值法補(bǔ)齊。

3 結(jié)果及分析

3.1 均值分析

2012—2021 年間，黃河流域生態(tài)保護(hù)、數(shù)字經(jīng)濟(jì)、文旅融合3 個(gè)系統(tǒng)的發(fā)展水平及綜合發(fā)展水平均呈上升趨勢(shì)（圖3）?？傮w來(lái)看，3 個(gè)系統(tǒng)的綜合發(fā)展指數(shù)由2012 年的0.28 增加至2021 年的0.35，總體發(fā)展較為平穩(wěn)。其中，2019—2021 年綜合發(fā)展指數(shù)因文旅融合降低而略有降低，2021 年綜合發(fā)展指數(shù)有較大回升。從各個(gè)系統(tǒng)來(lái)看，生態(tài)保護(hù)指數(shù)呈逐年上升趨勢(shì)，生態(tài)保護(hù)指數(shù)2019—2021 年年均增幅大于2012—2019 年，這是因?yàn)?019 年黃河流域的生態(tài)保護(hù)與高質(zhì)量發(fā)展被列入了國(guó)家的重要戰(zhàn)略，黃河流域各地市對(duì)生態(tài)保護(hù)投資力度加大，減污降碳力度增強(qiáng)，從而推動(dòng)了經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展；數(shù)字經(jīng)濟(jì)指數(shù)由2012 年的0.11 增加至2021 年的0.26，整體呈現(xiàn)上升趨勢(shì)，2012—2013 年增幅較小，其數(shù)值始終低于0.13，2014 年后增幅較為明顯，與我國(guó)實(shí)施“寬帶中國(guó)”政策密不可分，2017 年后，流域數(shù)字經(jīng)濟(jì)呈現(xiàn)較為穩(wěn)定的增長(zhǎng)態(tài)勢(shì)；文旅融合指數(shù)由2012 年的0.12 增加至2021 年的0.14，與另外兩個(gè)系統(tǒng)相比增幅較小但年度間波動(dòng)較大，受疫情沖擊及國(guó)際環(huán)境變化影響，2019—2021 年黃河流域各地市國(guó)內(nèi)外游客接待數(shù)量及旅游總收入明顯減少，2021 年文旅融合指數(shù)達(dá)到近6 年最低水平。

圖3 2012、2021 年黃河流域EDC趨勢(shì)面Figure 3 EDC trend surface of the Yellow River Basin in 2012 and 2021

3.2 耦合協(xié)調(diào)度趨勢(shì)及時(shí)空分析

運(yùn)用耦合協(xié)調(diào)度模型測(cè)算黃河流域EDC 系統(tǒng)的耦合協(xié)調(diào)度，借助ArcGIS 軟件分別繪制2012、2021 年黃河流域EDC 系統(tǒng)的趨勢(shì)面和時(shí)空演變圖（圖3）。其中，趨勢(shì)面的X軸表示東西方向，Y軸表示南北方向，Z軸表示耦合協(xié)調(diào)度。

從整體時(shí)空演化特征來(lái)看（圖4），總體上黃河流域EDC呈現(xiàn)出“東高西低、南高北低”的空間分布特征。從東西方向來(lái)看，2012 年?yáng)|西方向的趨勢(shì)線具有明顯的“東高西低”的特征，說(shuō)明黃河流域東部地區(qū)EDC系統(tǒng)耦合協(xié)調(diào)發(fā)展?fàn)顟B(tài)較佳；2021 年，東西趨勢(shì)線表現(xiàn)出“東西高中部低”的“U”型特征，說(shuō)明黃河流域中西部地區(qū)耦合協(xié)調(diào)度為黃河流域最低，同時(shí)西側(cè)趨勢(shì)線的坡度小于東側(cè)，說(shuō)明西部地區(qū)的耦合協(xié)調(diào)度差異小于東側(cè)。從南北方向來(lái)看，2012—2021 年南北趨勢(shì)線經(jīng)歷了從“凹”到“凸”的演變，整體表現(xiàn)為“南高北低”，說(shuō)明黃河流域南部EDC系統(tǒng)耦合協(xié)調(diào)水平高于北部地區(qū)，并且隨著時(shí)間推移南北地區(qū)耦合協(xié)調(diào)度差距增大。原因在于：首先，黃河中西部地區(qū)雖然擁有多樣的自然景觀與優(yōu)美的生態(tài)環(huán)境，但部分地市仍存在較為粗放的數(shù)字經(jīng)濟(jì)模式，數(shù)字經(jīng)濟(jì)水平相對(duì)較弱，省會(huì)及中心城市三系統(tǒng)發(fā)展水平與耦合協(xié)調(diào)水平不及東部，對(duì)周圍地市的溢出效應(yīng)及輻射帶動(dòng)作用不夠明顯；其次，黃河下游，即黃河的東南部地區(qū)是中國(guó)古文明的搖籃，孕育了仰韶文化、龍山文化等瑰麗的黃河文化，擁有眾多歷史悠久的古都城市，具有豐厚的文化底蘊(yùn)和文旅資源，文旅融合發(fā)展水平較高；最后，黃河?xùn)|部地區(qū)人口基數(shù)大，高校及科研院所眾多，較中西部地區(qū)具有更加充足的人才儲(chǔ)備及豐碩的科研成果，經(jīng)濟(jì)高質(zhì)量發(fā)展基礎(chǔ)更為堅(jiān)實(shí)。

圖4 2012、2021 年黃河流域EDC耦合系統(tǒng)的Kriging插值圖Figure 4 Kriging interpolation diagram map of EDC coupling system in the Yellow River Basin in 2012 and 2021

為打破行政區(qū)域界限，本文借助ArcGIS10.7 軟件繪制Kriging插值圖以便更直觀地對(duì)局部時(shí)空演化特征進(jìn)行分析（圖4）。2012 年，黃河流域EDC存在明顯的地區(qū)差異，各協(xié)調(diào)水平的城市數(shù)量分布近似于正態(tài)分布。西安、鄭州兩個(gè)省會(huì)城市憑借其文旅資源豐富、交通發(fā)達(dá)便利、基礎(chǔ)設(shè)施較完善等優(yōu)勢(shì)，其耦合協(xié)調(diào)度高于其他地市，分別達(dá)到0.65和0.64，率先進(jìn)入初級(jí)協(xié)調(diào)水平；淄博、濟(jì)南、太原、東營(yíng)、洛陽(yáng)、呼和浩特和濱州7 個(gè)地市緊隨其后，為勉強(qiáng)失調(diào)水平；超過(guò)50%的城市處于瀕臨失調(diào)水平，主要分布在內(nèi)蒙古和黃河流域中下游地區(qū)；位于黃河流域中上游的16 個(gè)城市EDC 較低，處于輕度失調(diào)水平。2021 年，黃河流域EDC顯著提高，有34 個(gè)地市的EDC實(shí)現(xiàn)了跨協(xié)調(diào)等級(jí)的提升，所有地市均達(dá)到瀕臨失調(diào)及以上水平。西安以0.70 的耦合協(xié)調(diào)度成為57 個(gè)地市中唯一1 個(gè)達(dá)到中級(jí)協(xié)調(diào)的地市，這與其豐富的文旅資源、發(fā)達(dá)的交通網(wǎng)絡(luò)、持續(xù)改善的生態(tài)環(huán)境與平穩(wěn)堅(jiān)實(shí)的經(jīng)濟(jì)基礎(chǔ)密不可分；其次，黃河下游的東營(yíng)、濟(jì)南、淄博、鄭州、洛陽(yáng)、濱州躋身初級(jí)協(xié)調(diào)水平；勉強(qiáng)失調(diào)水平的城市增加至15個(gè)，呈鑲嵌式分布特征，除太原、呼和浩特、蘭州、銀川這4 個(gè)省會(huì)城市外，其余多分布在鄭州與濟(jì)南周圍，由此可見(jiàn)，黃河流域下游EDC普遍較高；瀕臨失調(diào)的地市占比為61.4%，連片分布于黃河流域中上游，其中，甘肅與寧夏的地市耦合協(xié)調(diào)水平常年位于黃河流域57個(gè)地市排名的尾部，說(shuō)明這兩個(gè)省區(qū)未來(lái)應(yīng)緊抓生態(tài)文明建設(shè)，大力發(fā)展數(shù)字經(jīng)濟(jì)和文旅市場(chǎng)，盡快縮小與流域其他地市的耦合協(xié)調(diào)水平差距。

3.3 空間相關(guān)性分析

首先，運(yùn)用Global Moran′s I分析黃河流域EDC系統(tǒng)耦合協(xié)調(diào)水平的空間效應(yīng)，通過(guò)GeoDa 軟件對(duì)面板數(shù)據(jù)進(jìn)行計(jì)算得到黃河流域EDC 協(xié)調(diào)度Global Moran′s I，并采用MCMC 方法驗(yàn)證其顯著性狀況（表3）。結(jié)果顯示，黃河流域EDC的全局莫蘭指數(shù)隨時(shí)間推移而動(dòng)態(tài)變化且均在5%的水平上顯著為正，表明黃河流域EDC系統(tǒng)的耦合協(xié)調(diào)水平與鄰近城市的發(fā)展存在著正向關(guān)聯(lián)，Global Moran′s I 指數(shù)的變化趨勢(shì)表明，隨著距離城市的增加，EDC 系統(tǒng)的耦合協(xié)調(diào)水平也在不斷提高，這反映了流域內(nèi)地區(qū)之間的聯(lián)系日益緊密，協(xié)調(diào)發(fā)展的空間聚集特征也越來(lái)越明顯。

表3 黃河流域EDC系統(tǒng)耦合協(xié)調(diào)度的全局莫蘭指數(shù)Table 3 Global Moran′s I index of coupling coordination degree of EDC system in the Yellow River Basin

為進(jìn)一步探究黃河流域每個(gè)地市EDC 系統(tǒng)耦合協(xié)調(diào)發(fā)展與周邊城市的關(guān)系，根據(jù)Local Moran′s I計(jì)算結(jié)果繪制2012、2021 年的黃河流域EDC 的LISA聚類圖（圖6）。

由圖5 可知，黃河流域57 個(gè)地市間EDC 的集聚程度主要以“高高”集聚與“低低”集聚兩種類型為主導(dǎo)，具有較高的空間穩(wěn)定性，說(shuō)明黃河流域EDC系統(tǒng)耦合協(xié)調(diào)發(fā)展水平存在一定的兩極分化?！案吒摺奔蹍^(qū)域主要分布在黃河下游地區(qū)，2012—2021年“高高”集聚區(qū)的地市變動(dòng)較小，除開(kāi)封市由“高高”集聚變?yōu)椤安伙@著”外，山東的濟(jì)南、淄博、東營(yíng)、泰安和濱州以及河南鄭州6 個(gè)地市一直處于“高高”集聚區(qū)域。“低低”集聚區(qū)集中于黃河的中上游，2012 年以甘肅和寧夏兩省區(qū)的地市為主，2021 年蘭州EDC 與周圍地市相比較高，成為“高低”集聚區(qū)，而“低低”集聚區(qū)有向黃河中部遷移的趨勢(shì)，數(shù)量上由2012 年的9 個(gè)地市減少至2021 年的5 個(gè)地市，其中慶陽(yáng)、白銀、定西和吳忠脫離“低低”集聚區(qū)?！暗透摺奔叟c“高低”集聚區(qū)具有一定的空間流動(dòng)性且較為分散。2012 年的聊城形成了“低高”集聚區(qū)，說(shuō)明其EDC 與周邊地市相比仍有一定差距，而2021 年并未出現(xiàn)“低高”集聚區(qū)；中游大多數(shù)地市未呈現(xiàn)出明顯的集聚特征，2012 年的蘭州和銀川與2021 年的太原、呼和浩特和西安形成了“高低”集聚區(qū)，說(shuō)明這5 個(gè)地市的EDC 雖然較高，但對(duì)周圍城市輻射帶動(dòng)能力較弱。

圖5 2012、2021 年黃河流域EDC的LISA聚類圖Figure 5 LISA Cluster Map of EDC in the Yellow River Basin in 2021

3.4 障礙因子分析

為更細(xì)致地探究影響EDC 系統(tǒng)耦合協(xié)調(diào)水平的因素，借助障礙度模型探究影響2012—2021 年黃河流域EDC系統(tǒng)耦合協(xié)調(diào)發(fā)展差異的指標(biāo)層障礙因子。由于涉及的具體指標(biāo)與年份較多，因此研究采用2012—2021 年各指標(biāo)的障礙度均值代表各指標(biāo)最終的障礙度并進(jìn)行排序，本文將各地市各年份障礙度大于5 且排名前3 的障礙因子稱為其主要障礙因子（表4）。

表4 黃河流域EDC系統(tǒng)指標(biāo)層主要障礙因子Table 4 Main obstacle factors of the EDC system indicator layer in the Yellow River Basin

從生態(tài)保護(hù)系統(tǒng)來(lái)看，黃河流域所轄地市主要障礙因子均包括人均公園綠地面積（X12）和空氣質(zhì)量?jī)?yōu)良率（X13），說(shuō)明以人均公園綠地面積與空氣質(zhì)量?jī)?yōu)良率為代表的城市生態(tài)環(huán)境建設(shè)情況對(duì)黃河流域生態(tài)保護(hù)的阻礙程度較高，其可能的原因是黃河流域工業(yè)化快速發(fā)展及建設(shè)用地不斷擴(kuò)張?jiān)谝欢ǔ潭壬舷拗屏顺鞘猩鷳B(tài)環(huán)境建設(shè)；同時(shí)隨著城鎮(zhèn)化的不斷發(fā)展以及陸續(xù)開(kāi)放“二孩政策”與“三孩政策”使得城市人口密度增加，城市生態(tài)環(huán)境承載力降低，從而對(duì)生態(tài)保護(hù)發(fā)展產(chǎn)生制約作用。除此之外，有38%的地市第三障礙因子為化肥施用量（X22），這些地市均位于黃淮海平原和汾渭平原等重要的農(nóng)業(yè)生產(chǎn)區(qū)，這提醒黃河流域各地市應(yīng)注重提升農(nóng)田質(zhì)量，推動(dòng)綠色農(nóng)田建設(shè)及農(nóng)業(yè)可持續(xù)發(fā)展；而內(nèi)蒙古和山西的大部分地市第三障礙因子均為萬(wàn)元GDP 能耗（X21），這表明黃河流域生態(tài)保護(hù)系統(tǒng)的障礙因子存在區(qū)域異質(zhì)性。

從數(shù)字經(jīng)濟(jì)系統(tǒng)來(lái)看，有64%的地市主要障礙因子及其排序?yàn)槊壳嘶ヂ?lián)網(wǎng)寬帶接入用戶數(shù)（Y12）＞電信業(yè)務(wù)總量占GDP 比例（Y23）＞互聯(lián)網(wǎng)普及率（Y13），每千人互聯(lián)網(wǎng)寬帶接入用戶數(shù)與互聯(lián)網(wǎng)普及率反映了數(shù)字化程度，電信業(yè)務(wù)總量占GDP比例反映了黃河流域經(jīng)濟(jì)發(fā)展趨勢(shì)，由此可知，黃河流域各地市面臨著數(shù)字化基礎(chǔ)薄弱以及社會(huì)經(jīng)濟(jì)發(fā)展不平衡不充分的尷尬窘境，嚴(yán)重制約了數(shù)字經(jīng)濟(jì)的發(fā)展。

從文旅融合系統(tǒng)來(lái)看，有53 個(gè)地市的第一障礙因子為公路客運(yùn)量（Z33），障礙度均值為21.06%，表明黃河流域的公路交通建設(shè)限制了黃河流域文旅產(chǎn)業(yè)的發(fā)展；90%的地市主要障礙因子還包括A 級(jí)以上景區(qū)密度（Z12）和國(guó)內(nèi)外游客接待數(shù)（Z22），這兩個(gè)指標(biāo)一定程度上反映了當(dāng)?shù)匚穆觅Y源豐富程度和城市吸引力，說(shuō)明雖然黃河流域是中華文明的發(fā)源地，自然景觀多樣，文旅資源豐富，但是當(dāng)前黃河流域各地市仍然缺乏對(duì)當(dāng)?shù)匚幕纳钊肜斫庖约皩?duì)內(nèi)涵的挖掘，文旅宣傳力度亟需加強(qiáng)。

4 結(jié)論及政策建議

4.1 結(jié)論

本文測(cè)算了黃河流域EDC耦合協(xié)調(diào)水平，并對(duì)其時(shí)空演化特征、空間集聚特征及障礙因子開(kāi)展了研究。主要結(jié)論如下：①時(shí)間維度上，2012—2021年黃河流域EDC各子系統(tǒng)均保持上升趨勢(shì)，三者的耦合協(xié)調(diào)水平有顯著提高，實(shí)現(xiàn)了由輕度失調(diào)向?yàn)l臨失調(diào)的轉(zhuǎn)變，各地市呈現(xiàn)出向更高協(xié)調(diào)水平集聚的發(fā)展趨勢(shì)。②空間維度上，黃河流域EDC協(xié)調(diào)度呈現(xiàn)出“東高西低、南高北低”的空間分布特征，東西向表現(xiàn)為東部地區(qū)地市間耦合協(xié)調(diào)差異大于西部地區(qū)，南北向表現(xiàn)為10 年間南北地區(qū)耦合協(xié)調(diào)度差異有所增大。③從空間相關(guān)性的分析來(lái)看，黃河流域EDC耦合協(xié)調(diào)度具有正向關(guān)聯(lián)性和空間集聚性，集聚程度主要以“高高”集聚與“低低”集聚兩種類型為主導(dǎo)，具有較高的空間穩(wěn)定性。④從障礙因子的分析來(lái)看，影響黃河流域生態(tài)保護(hù)的主要障礙因子是人均公園綠地面積、空氣質(zhì)量?jī)?yōu)良率和化肥施用量，生態(tài)保護(hù)障礙因子具有明顯的地區(qū)差異；影響黃河流域數(shù)字經(jīng)濟(jì)的主要障礙因子是每千人互聯(lián)網(wǎng)寬帶接入用戶數(shù)、電信業(yè)務(wù)總量占GDP比例以及互聯(lián)網(wǎng)普及率；影響黃河流域文旅融合的主要障礙因子是公路客運(yùn)量、A 級(jí)以上景區(qū)密度和國(guó)內(nèi)外游客接待數(shù)。

4.2 政策建議

黃河流域各市應(yīng)在繼續(xù)以生態(tài)環(huán)境為著力點(diǎn)，依托數(shù)字經(jīng)濟(jì)，在實(shí)現(xiàn)文旅融合發(fā)展的基礎(chǔ)上，加大對(duì)數(shù)字化經(jīng)濟(jì)方面的支持，實(shí)施生態(tài)環(huán)境、數(shù)字經(jīng)濟(jì)、文旅的有效融合，深入推進(jìn)數(shù)字化生態(tài)、文旅創(chuàng)新體系的建設(shè)。同時(shí)政府要根據(jù)我國(guó)的基本國(guó)情以及協(xié)調(diào)發(fā)展?fàn)顟B(tài)，完善出臺(tái)一系列實(shí)際可行、高效的協(xié)調(diào)發(fā)展規(guī)劃，并構(gòu)建完整的生態(tài)、數(shù)字經(jīng)濟(jì)、文旅耦合協(xié)調(diào)發(fā)展的統(tǒng)籌機(jī)制，進(jìn)一步推動(dòng)生態(tài)環(huán)境、數(shù)字經(jīng)濟(jì)、文旅的耦合協(xié)調(diào)。堅(jiān)持“提質(zhì)增效、重點(diǎn)發(fā)展”的原則繼續(xù)助力經(jīng)濟(jì)欠發(fā)達(dá)地區(qū)與較發(fā)達(dá)地區(qū)的協(xié)調(diào)發(fā)展。

采取差異化對(duì)策實(shí)現(xiàn)空間協(xié)同發(fā)展，打造區(qū)域協(xié)調(diào)發(fā)展體系?；诘赜騼?yōu)勢(shì)，東部地區(qū)應(yīng)積極推動(dòng)發(fā)展，不斷鞏固成果，并充分利用資源、技術(shù)和經(jīng)驗(yàn)；西部地區(qū)的協(xié)調(diào)度較低，要補(bǔ)齊短板，加速推進(jìn)傳統(tǒng)產(chǎn)業(yè)的轉(zhuǎn)型升級(jí)，促進(jìn)新興產(chǎn)業(yè)的發(fā)展，同時(shí)結(jié)合地域特色，推進(jìn)“生態(tài)＋”的發(fā)展模式，依托數(shù)字經(jīng)濟(jì)打造具有影響力的黃河文化生態(tài)走廊和黃河文化旅游目的地品牌，從而帶動(dòng)黃河流域生態(tài)環(huán)境、數(shù)字經(jīng)濟(jì)和文旅的整體協(xié)調(diào)發(fā)展。耦合協(xié)調(diào)度高的地區(qū)應(yīng)該充分利用其溢出效應(yīng)，并將其輻射到周邊地區(qū)。

疏通耦合協(xié)調(diào)發(fā)展堵點(diǎn)，提高耦合協(xié)調(diào)發(fā)展水平?；趯?duì)黃河流域生態(tài)—數(shù)字經(jīng)濟(jì)—文旅系統(tǒng)耦合協(xié)調(diào)發(fā)展障礙因子的分析，為推動(dòng)黃河流域高質(zhì)量發(fā)展，首先要提升沿黃城市生態(tài)建設(shè)水平，以綠色發(fā)展為引領(lǐng)，加快產(chǎn)業(yè)布局優(yōu)化與結(jié)構(gòu)調(diào)整，推動(dòng)產(chǎn)業(yè)全面綠色低碳轉(zhuǎn)型。其次要推進(jìn)多領(lǐng)域協(xié)同創(chuàng)新，進(jìn)一步突出黃河流域科技創(chuàng)新聯(lián)盟作用，建設(shè)豫魯沿黃科創(chuàng)大走廊，爭(zhēng)創(chuàng)綜合性國(guó)家科學(xué)中心，營(yíng)造優(yōu)良創(chuàng)新生態(tài)；最后要提升對(duì)外開(kāi)放水平，依托自貿(mào)試驗(yàn)區(qū)、跨境電商綜試區(qū)、市場(chǎng)采購(gòu)貿(mào)易方式試點(diǎn)等，先行先試、聯(lián)動(dòng)創(chuàng)新，持續(xù)優(yōu)化營(yíng)商環(huán)境。