關(guān)鍵詞：生態(tài)韌性；經(jīng)濟(jì)高質(zhì)量發(fā)展；協(xié)同；哈肯模型；黃河流域

0引言

高質(zhì)量發(fā)展是全面建設(shè)社會(huì)主義現(xiàn)代化國(guó)家的首要任務(wù)，經(jīng)濟(jì)高質(zhì)量發(fā)展就是要實(shí)現(xiàn)質(zhì)的有效提升和量的合理增長(zhǎng)。黃河流域自然資源較為豐富，是我國(guó)重要的工業(yè)與糧食基地，經(jīng)濟(jì)發(fā)展主要依賴重工業(yè)與能源，長(zhǎng)期形成的經(jīng)濟(jì)發(fā)展伴隨“高能耗、高污染”狀況，使得原本脆弱、退化的生態(tài)系統(tǒng)承受極大壓力，生態(tài)系統(tǒng)的恢復(fù)、保護(hù)、治理難度極大，如何取得經(jīng)濟(jì)高質(zhì)量發(fā)展與生態(tài)韌性提升的雙贏是當(dāng)前亟須解決的重大問(wèn)題。

經(jīng)濟(jì)高質(zhì)量發(fā)展是目前學(xué)界研究的熱點(diǎn)之一，研究?jī)?nèi)容主要集中在三方面：一是經(jīng)濟(jì)高質(zhì)量發(fā)展的影響因素，如社會(huì)資本、人力資本、網(wǎng)絡(luò)基礎(chǔ)設(shè)施、金融科技、國(guó)家創(chuàng)新型城市建設(shè)等，采用的計(jì)量模型主要有多期雙重差分模型、動(dòng)態(tài)空間自回歸模型、門檻模型等；二是經(jīng)濟(jì)高質(zhì)量發(fā)展水平測(cè)度與評(píng)價(jià)，運(yùn)用的方法主要有主成分分析法、熵權(quán)法、綜合評(píng)價(jià)法等；三是經(jīng)濟(jì)高質(zhì)量發(fā)展的時(shí)代背景，如中國(guó)式現(xiàn)代化、“一帶一路”建設(shè)、數(shù)字經(jīng)濟(jì)、新征程、綠色發(fā)展等。生態(tài)韌性也是學(xué)界的研究熱點(diǎn)之一，有關(guān)學(xué)者主要聚焦于生態(tài)韌性概念及其評(píng)估開(kāi)展研究，多數(shù)學(xué)者把生態(tài)韌性定義為生態(tài)系統(tǒng)遭受外部干擾后恢復(fù)到原來(lái)狀態(tài)或形成新平衡狀態(tài)的能力，生態(tài)韌性評(píng)價(jià)方法與經(jīng)濟(jì)高質(zhì)量發(fā)展評(píng)價(jià)方法類似，評(píng)估框架包括PSR、抵抗力一適應(yīng)力一恢復(fù)力、自然一經(jīng)濟(jì)一社會(huì)等。有關(guān)學(xué)者開(kāi)展了經(jīng)濟(jì)發(fā)展與生態(tài)韌性關(guān)系的研究，如郭力等研究了數(shù)字經(jīng)濟(jì)發(fā)展對(duì)黃河流域生態(tài)韌性的影響、蔣文鑫等對(duì)江蘇省城市生態(tài)韌性與經(jīng)濟(jì)發(fā)展水平的耦合協(xié)調(diào)關(guān)系進(jìn)行了測(cè)算、趙領(lǐng)娣等測(cè)算了我國(guó)53個(gè)海岸帶城市經(jīng)濟(jì)發(fā)展與生態(tài)韌性的耦合協(xié)調(diào)度。

梳理相關(guān)文獻(xiàn)表明，學(xué)界對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展和生態(tài)韌性在概念界定、影響因素、測(cè)度與評(píng)價(jià)等方面的研究成果較多，對(duì)于二者關(guān)系的研究尚處于起步階段，鮮有對(duì)黃河流域生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同關(guān)系及其演化機(jī)制的定量研究。鑒于此，本文在生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同演化理論分析基礎(chǔ)上，分別構(gòu)建黃河流域生態(tài)韌性、經(jīng)濟(jì)高質(zhì)量發(fā)展能力評(píng)價(jià)指標(biāo)體系，采用熵權(quán)-TOPSIS法對(duì)黃河流域各省（區(qū)）生態(tài)韌性、經(jīng)濟(jì)高質(zhì)量發(fā)展能力進(jìn)行測(cè)度，運(yùn)用哈肯模型探究黃河流域生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同關(guān)系及其演化規(guī)律，以期為黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展的全面推進(jìn)及相關(guān)研究提供參考。

1理論分析與研究方法

1.1生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同演化理論分析

由德國(guó)物理學(xué)家哈肯創(chuàng)立的協(xié)同學(xué)，是研究復(fù)雜系統(tǒng)自組織演化的學(xué)科，其認(rèn)為相互協(xié)同的兩個(gè)或多個(gè)子系統(tǒng)構(gòu)成的復(fù)合系統(tǒng)，具有開(kāi)放性、自適應(yīng)性、非平衡性及非線性等特征。黃河流域生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展的協(xié)同，就構(gòu)成了這樣的復(fù)合系統(tǒng)（稱為生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)）。其中：系統(tǒng)開(kāi)放性是指各子系統(tǒng)間、系統(tǒng)與外界可以進(jìn)行物質(zhì)、能量與信息交換；系統(tǒng)自適應(yīng)性是指在外界發(fā)生變化使信息、物質(zhì)和能量大量涌人系統(tǒng)時(shí)，各子系統(tǒng)及其組成部分能夠整合、利用資源并推動(dòng)系統(tǒng)發(fā)展，主動(dòng)適應(yīng)外界的變化：系統(tǒng)非平衡性是指各子系統(tǒng)及其組成要素等存在時(shí)空差異，如黃河流域各?。▍^(qū)）之間存在差異、不同時(shí)期也存在差異；系統(tǒng)非線性是指各子系統(tǒng)相互作用及其效應(yīng)是復(fù)雜多變的。

在生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)中，生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展能力為狀態(tài)參量，其中：生態(tài)韌性是保障經(jīng)濟(jì)高質(zhì)量發(fā)展的物質(zhì)基礎(chǔ)，當(dāng)生態(tài)韌性增強(qiáng)時(shí)，可通過(guò)抵抗力、適應(yīng)力和恢復(fù)力推動(dòng)經(jīng)濟(jì)持續(xù)發(fā)展（擴(kuò)大經(jīng)濟(jì)規(guī)模、提高經(jīng)濟(jì)效率、完善經(jīng)濟(jì)結(jié)構(gòu)等），從而提高黃河流域經(jīng)濟(jì)高質(zhì)量發(fā)展能力：經(jīng)濟(jì)高質(zhì)量發(fā)展通過(guò)正反饋機(jī)制（經(jīng)濟(jì)支撐、環(huán)境規(guī)制等）為生態(tài)韌性增強(qiáng)（完善生態(tài)本底、提高治理能力、推動(dòng)生態(tài)修復(fù)等）創(chuàng)造條件，當(dāng)經(jīng)濟(jì)發(fā)展到一定規(guī)模時(shí)對(duì)生態(tài)韌性提出更高的要求。綜上，生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)中的兩個(gè)狀態(tài)參量受正負(fù)反饋機(jī)制的作用，呈現(xiàn)出不同速度、不同幅度的隨機(jī)波動(dòng)，其中之一通過(guò)自組織行為形成控制復(fù)合系統(tǒng)結(jié)構(gòu)和動(dòng)態(tài)演化方向的序參量；復(fù)合系統(tǒng)在序參量驅(qū)動(dòng)下，不斷跳躍結(jié)構(gòu)突變的臨界點(diǎn)，形成新的有序結(jié)構(gòu)，其有序程度越高協(xié)同程度就越高，如圖1所示。

基于黃河流域生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同演化的自組織特征，本文利用哈肯模型對(duì)兩個(gè)狀態(tài)參量的協(xié)同演化規(guī)律進(jìn)行分析，識(shí)別生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)的序參量，揭示協(xié)同演化機(jī)制，確定復(fù)合系統(tǒng)有序發(fā)展的平衡點(diǎn)。

1.2評(píng)價(jià)指標(biāo)體系構(gòu)建

1）生態(tài)韌性評(píng)價(jià)指標(biāo)體系。依據(jù)《城市生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)技術(shù)指南》（GB/T 43236-2023）、借鑒有關(guān)學(xué)者的研究，從系統(tǒng)韌性演化（在受到外界擾動(dòng)前進(jìn)行抵抗、受到擾動(dòng)時(shí)進(jìn)行適應(yīng)、擾動(dòng)結(jié)束后進(jìn)行恢復(fù)）視角，遵循合理性、全面性和科學(xué)性原則，按照抵抗力、適應(yīng)力、恢復(fù)力3個(gè)維度構(gòu)建生態(tài)韌性評(píng)價(jià)指標(biāo)體系，見(jiàn)表1。其中：抵抗力是指生態(tài)系統(tǒng)抵抗外界干擾的能力，與生態(tài)系統(tǒng)自身屬性和受到的擾動(dòng)情況有關(guān)，用自身屬性、水污染程度、大氣污染程度、固體污染程度、垃圾污染程度5個(gè)指標(biāo)來(lái)表征：適應(yīng)力表現(xiàn)為生態(tài)系統(tǒng)受到干擾時(shí)社會(huì)主體采取多種措施使生態(tài)系統(tǒng)實(shí)現(xiàn)對(duì)干擾的適應(yīng)，側(cè)重于生態(tài)系統(tǒng)中的社會(huì)屬性且與抵抗力中的污染指標(biāo)對(duì)應(yīng)，用污水治理能力、固體治理能力、廢氣治理能力、垃圾治理能力、生態(tài)建設(shè)投入5個(gè)指標(biāo)來(lái)表征：恢復(fù)力為擾動(dòng)結(jié)束后生態(tài)系統(tǒng)恢復(fù)到原有狀態(tài)或形成新的穩(wěn)定狀態(tài)的能力，用綠化修復(fù)程度、公園修復(fù)程度、道路修復(fù)程度、排水管道修復(fù)程度、供水修復(fù)程度5個(gè)指標(biāo)來(lái)表征。

2）經(jīng)濟(jì)高質(zhì)量發(fā)展能力評(píng)價(jià)指標(biāo)體系。根據(jù)黨的十九大報(bào)告對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展的要求、借鑒有關(guān)學(xué)者的研究，按照經(jīng)濟(jì)規(guī)模、效率、結(jié)構(gòu)、持續(xù)性4個(gè)維度構(gòu)建經(jīng)濟(jì)高質(zhì)量發(fā)展能力評(píng)價(jià)指標(biāo)體系，見(jiàn)表2。其中：規(guī)模指經(jīng)濟(jì)體量的大小，包括經(jīng)濟(jì)量、消費(fèi)量、財(cái)政收入3個(gè)指標(biāo)，為使各?。▍^(qū)）指標(biāo)值具有可比性，用相對(duì)（人均）指標(biāo)來(lái)表征；效率指經(jīng)濟(jì)轉(zhuǎn)化效率，包括金融效率、勞動(dòng)效率、固定資產(chǎn)投資效率3個(gè)指標(biāo)；結(jié)構(gòu)包含產(chǎn)業(yè)結(jié)構(gòu)、城鄉(xiāng)結(jié)構(gòu)、經(jīng)濟(jì)開(kāi)放結(jié)構(gòu)、消費(fèi)結(jié)構(gòu)4個(gè)指標(biāo)：持續(xù)性指經(jīng)濟(jì)發(fā)展能力，包括信息化程度、創(chuàng)新能力、整體活力3個(gè)指標(biāo)。

3）數(shù)據(jù)來(lái)源。本文以黃河流域九?。▍^(qū)）為研究單元、以2012-2022年為研究期，各指標(biāo)數(shù)據(jù)來(lái)源于歷年《中國(guó)統(tǒng)計(jì)年鑒》《中國(guó)城市建設(shè)統(tǒng)計(jì)年鑒》及各?。▍^(qū)）《統(tǒng)計(jì)年鑒》、統(tǒng)計(jì)局網(wǎng)站等，缺失數(shù)據(jù)采用線性插值法補(bǔ)全。

1.3研究方法

1.3.1生態(tài)韌性和經(jīng)濟(jì)高質(zhì)量發(fā)展能力測(cè)度

熵權(quán)法能夠充分利用各指標(biāo)原始數(shù)據(jù)包含的信息、減少主觀賦權(quán)法造成的偏差；TOPSIS法主要應(yīng)用于多指標(biāo)綜合評(píng)估，用各樣本值與最優(yōu)、最劣理想解的差距計(jì)算貼進(jìn)度（即與最優(yōu)理想解的貼進(jìn)度），具體方法參見(jiàn)文獻(xiàn)。本研究采用熵權(quán)法計(jì)算生態(tài)韌性評(píng)價(jià)指標(biāo)和經(jīng)濟(jì)高質(zhì)量發(fā)展能力評(píng)價(jià)指標(biāo)的權(quán)重，采用TOPSIS法計(jì)算黃河流域各?。▍^(qū)）歷年生態(tài)韌性貼進(jìn)度、經(jīng)濟(jì)高質(zhì)量發(fā)展能力貼進(jìn)度，并將其分別作為生態(tài)韌性、經(jīng)濟(jì)高質(zhì)量發(fā)展能力（J）。

1.3.2生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同水平測(cè)度

哈肯模型將系統(tǒng)狀態(tài)參量分為快變量和慢變量?jī)深?，通過(guò)構(gòu)建方程計(jì)算快變量的系統(tǒng)不平衡點(diǎn)，并運(yùn)用熱學(xué)中的絕熱近似原理消去快變量，得到慢變量方程和演化方程組，進(jìn)而識(shí)別系統(tǒng)發(fā)展的序參量，并研判系統(tǒng)自組織過(guò)程中的演化機(jī)制與演化趨勢(shì)。把哈肯模型用于研究黃河流域生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同演化機(jī)制，具有科學(xué)性和可行性。

物理方程是針對(duì)連續(xù)型隨機(jī)變量設(shè)定的，為了將其用于生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)研究，對(duì)式（1）、式（2）分別按式（7）、式（8）進(jìn)行離散化處理后，進(jìn)行序參量、勢(shì)函數(shù)、協(xié)同穩(wěn)定點(diǎn)、樣本點(diǎn)與協(xié)同穩(wěn)定點(diǎn)距離等的分析計(jì)算。

對(duì)樣本點(diǎn)與協(xié)同穩(wěn)定點(diǎn)的距離按下式進(jìn)行正向歸一化處理，作為衡量生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)協(xié)同水平的指標(biāo)：

3實(shí)證研究

3.1生態(tài)韌性、經(jīng)濟(jì)高質(zhì)量發(fā)展能力測(cè)算結(jié)果

2012-2022年黃河流域各省（區(qū)）生態(tài)韌性、經(jīng)濟(jì)高質(zhì)量發(fā)展能力測(cè)算結(jié)果分別見(jiàn)表3、表4。

由表3、表4可知：2012-2022年，黃河流域生態(tài)韌性和經(jīng)濟(jì)高質(zhì)量發(fā)展能力在整體上保持穩(wěn)步上升的協(xié)同發(fā)展趨勢(shì)，生態(tài)韌性與新經(jīng)濟(jì)高質(zhì)量發(fā)展能力相比上升速度更快、上升幅度更大，經(jīng)濟(jì)高質(zhì)量發(fā)展能力在2020-2022年發(fā)生小幅波動(dòng)；黃河流域九?。▍^(qū)）生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展能力均表現(xiàn)出上升趨勢(shì)，生態(tài)韌性上升幅度均比經(jīng)濟(jì)高質(zhì)量發(fā)展能力的上升幅度大。

3.2生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同演化實(shí)證分析

1）序參量識(shí)別。基于2012-2022年黃河流域九?。▍^(qū)）生態(tài)韌性和經(jīng)濟(jì)高質(zhì)量發(fā)展能力測(cè)算結(jié)果，分別假設(shè)生態(tài)韌性和經(jīng)濟(jì)高質(zhì)量發(fā)展能力為序參量，運(yùn)用經(jīng)濟(jì)計(jì)量分析軟件EViews按式（7）和式（8）進(jìn)行回歸，結(jié)果見(jiàn)表5。假設(shè)生態(tài)韌性為序參量時(shí)，控制參量y1=-0. 029 7，y2=-0.0227，a=0.123 0，6=0.0054，絕熱近似假設(shè)（原假設(shè)）不成立；假設(shè)經(jīng)濟(jì)高質(zhì)量發(fā)展能力為序參量時(shí)，控制參量，且Y2gt;0，絕熱近似假設(shè)（原假設(shè)）成立，即黃河流域經(jīng)濟(jì)高質(zhì)量發(fā)展能力為序參量。

對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展能力演化方程進(jìn)行變換，得到如下勢(shì)函數(shù)：

3）生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同演化狀況。上述黃河流域生態(tài)韌性和經(jīng)濟(jì)高質(zhì)量發(fā)展能力測(cè)算結(jié)果及哈肯模型回歸結(jié)果表明，黃河流域生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展具有協(xié)同演化特征，在演化過(guò)程中經(jīng)濟(jì)高質(zhì)量發(fā)展為序參量，即經(jīng)濟(jì)高質(zhì)量發(fā)展發(fā)揮主導(dǎo)作用。從哈肯模型控制參量看，a= 0.031 2gt;0，說(shuō)明生態(tài)韌性對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展具有抑制作用，大力提高生態(tài)韌性不利于經(jīng)濟(jì)高質(zhì)量發(fā)展，但a數(shù)值較小，說(shuō)明其抑制作用較?。篵=0.4696gt;0，表明經(jīng)濟(jì)高質(zhì)量發(fā)展對(duì)生態(tài)韌性具有協(xié)同提升效應(yīng)，且6數(shù)值較大，表明經(jīng)濟(jì)高質(zhì)量發(fā)展對(duì)生態(tài)韌性提升效果明顯；Y1=-0.0360lt;0，說(shuō)明黃河流域生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)有序演化過(guò)程中經(jīng)濟(jì)高質(zhì)量發(fā)展已經(jīng)對(duì)生態(tài)韌性形成正反饋機(jī)制，但Y1絕對(duì)值較小，生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)尚處于低級(jí)有序階段；y2=0.1604gt;0，說(shuō)明黃河流域生態(tài)韌性對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展具有負(fù)反饋機(jī)制。

4）生態(tài)韌性與經(jīng)濟(jì)高質(zhì)量發(fā)展協(xié)同水平。把上述生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)協(xié)同穩(wěn)定點(diǎn)（0.6278，-0.0229）代入式（9），測(cè)算2012-2022年黃河流域九?。▍^(qū)）生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)協(xié)同水平（S值），結(jié)果見(jiàn)表6。依據(jù)S值，把生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)協(xié)同水平分為低（S=[0，0.333]）、中（S=（0.333， 0.667]）、高（S=（0.667，1]）3級(jí)，2012-2022年各級(jí)協(xié)同水平的省份數(shù)量統(tǒng)計(jì)見(jiàn)表7。

由表6、表7可知，黃河流域生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)整體協(xié)同水平（流域均值）在2012-2019年呈穩(wěn)步上升趨勢(shì)、在2020-2022年出現(xiàn)小幅波動(dòng)：研究期低協(xié)同水平的省份數(shù)量由8個(gè)減少為1個(gè)，中、高協(xié)同水平的省份數(shù)量分別由1、0個(gè)增加為5、3個(gè)，各?。▍^(qū)）協(xié)同水平的差異逐步縮小，其中山東一直維持中、高水平，而青海、甘肅、寧夏、山西一直在中、低水平徘徊。

4結(jié)論

1）黃河流域整體生態(tài)韌性在研究期高于經(jīng)濟(jì)高質(zhì)量發(fā)展能力且呈現(xiàn)穩(wěn)步上升趨勢(shì)，經(jīng)濟(jì)高質(zhì)量發(fā)展能力在2012-2019年穩(wěn)步上升、在2020-2022年出現(xiàn)小幅波動(dòng)。

2）黃河流域生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)尚處于低級(jí)有序階段，在協(xié)同演化過(guò)程中，經(jīng)濟(jì)高質(zhì)量發(fā)展能力為序參量、主導(dǎo)二者協(xié)同發(fā)展演化的路徑與方向并對(duì)生態(tài)韌性具有協(xié)同提升效應(yīng)，而生態(tài)韌性對(duì)經(jīng)濟(jì)高質(zhì)量發(fā)展有一定抑制作用但其擬制作用較小。

3）黃河流域生態(tài)韌性一經(jīng)濟(jì)高質(zhì)量發(fā)展復(fù)合系統(tǒng)協(xié)同水平在研究期呈現(xiàn)穩(wěn)步上升趨勢(shì)，在空間上呈現(xiàn)東高西低的格局，流域九?。▍^(qū)）協(xié)同水平在2018年以前差異較大、2018年以后差異縮小。