摘 要：水資源效率與韌性是可持續(xù)性研究的2個(gè)核心問題，分析中國水資源效率與韌性的協(xié)同發(fā)展關(guān)系對(duì)促進(jìn)中國水資源可持續(xù)利用具有重要意義。以2000—2022年全國30個(gè)省區(qū)為研究對(duì)象，運(yùn)用超效率SBM模型和三維矢量模型測(cè)度了中國水資源效率和韌性水平，并使用哈肯模型探究了兩者的協(xié)同演化特征。結(jié)果表明：研究期內(nèi)水資源效率呈波動(dòng)式上升—緩慢下降趨勢(shì)，2022年中國區(qū)域水資源效率均值為0. 47，整體上處于非有效區(qū)；研究期內(nèi)中國水資源韌性水平整體呈平穩(wěn)波動(dòng)—快速上升—緩慢上升趨勢(shì)，區(qū)域內(nèi)水資源韌性水平存在梯度差異，表現(xiàn)為東部地區(qū)>東北地區(qū)>中部地區(qū)>西部地區(qū)；研究期內(nèi)中國水資源韌性與效率協(xié)同水平整體上呈波動(dòng)式上升趨勢(shì)，水資源效率為協(xié)同演化的序參量，對(duì)整個(gè)系統(tǒng)的演化進(jìn)程起積極導(dǎo)向作用。

關(guān)鍵詞：水資源；韌性；效率；協(xié)同度

中圖分類號(hào)：TV213 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1001-9235（2024）09-0047-09

當(dāng)前，中國正面臨著水資源供需矛盾突出、水環(huán)境污染嚴(yán)重、水資源利用率低等問題，對(duì)中國水資源系統(tǒng)可持續(xù)發(fā)展帶來了挑戰(zhàn)［1］。為了應(yīng)對(duì)日益嚴(yán)峻的水資源挑戰(zhàn)，有必要深入研究水資源系統(tǒng)效率與韌性水平，以期提高其對(duì)外界環(huán)境變化的抵抗能力、適應(yīng)能力和恢復(fù)能力。以往研究多集中于水資源效率測(cè)度及其影響因素分析，但單純強(qiáng)調(diào)水資源效率可能會(huì)使管理者過度關(guān)注資源利用效率的提高，忽視系統(tǒng)內(nèi)部的復(fù)雜性和穩(wěn)定性。因此，為了更好地提高水資源系統(tǒng)穩(wěn)定性，增強(qiáng)其對(duì)外界環(huán)境的適應(yīng)能力，本文探究了中國區(qū)域水資源效率和韌性水平，并分析了兩者間的協(xié)同關(guān)系，對(duì)推動(dòng)中國水資源系統(tǒng)朝著可持續(xù)方向發(fā)展具有重要意義。

目前，國內(nèi)外對(duì)水資源利用效率的研究相對(duì)成熟，主要集中在城市用水效率［2-3］、工農(nóng)業(yè)用水效率［4-6］方面。相關(guān)學(xué)者主要采用數(shù)據(jù)包絡(luò)分析法（DEA）［6-7］、隨機(jī)前沿分析法（SFA）［8-9］、主成分分析法［10］等對(duì)水資源利用效率進(jìn)行分析探究。水資源效率指標(biāo)主要有單要素指標(biāo)［11-13］與全要素指標(biāo)［14-16］。單要素指標(biāo)體系常用于農(nóng)業(yè)領(lǐng)域，評(píng)價(jià)指標(biāo)主要有萬元DGP 用水量、萬元工業(yè)增加值用水量、畝均灌溉用水量等。全要素指標(biāo)體系主要從資源消耗、經(jīng)濟(jì)社會(huì)發(fā)展、生態(tài)環(huán)境等方面構(gòu)建。目前測(cè)算水資源效率使用最多的是DEA，這是一種計(jì)算多決策單元效率的非參數(shù)方法，傳統(tǒng)的DEA模型忽略了松弛變量對(duì)水資源效率的影響，容易高估DMU（決策單元）效率值，而超效率SBM可避免由于投入產(chǎn)出變量松弛性導(dǎo)致效率評(píng)價(jià)不準(zhǔn)的問題，且決策單元不受［0，1］的限制。本文以水資源、人力、資本為投入要素，經(jīng)濟(jì)發(fā)展、生態(tài)保護(hù)為產(chǎn)出要素，并選用超效率SBM模型對(duì)中國30個(gè)省區(qū)的水資源效率進(jìn)行測(cè)算。

韌性是指系統(tǒng)在受到外部干擾之后不斷調(diào)整自身，使其恢復(fù)到之前的狀態(tài)［17］。眾多學(xué)者將韌性理論應(yīng)用到經(jīng)濟(jì)學(xué)、城市規(guī)劃、地理學(xué)等學(xué)科領(lǐng)域。國外學(xué)者對(duì)水資源韌性的研究主要集中在韌性理論構(gòu)建［18］、策略研究及影響因素［19-20］等方面。Sandoval- Soils等［18］以基于可持續(xù)發(fā)展指數(shù)，從可靠性、恢復(fù)性和脆弱性角度構(gòu)建了水資源指標(biāo)體系，探究了美國和墨西哥的水資源可持續(xù)性狀態(tài)。Schlvter等［19］基于生態(tài)保護(hù)視角，將水環(huán)境流量納入徑流指標(biāo)體系中，評(píng)估了阿穆達(dá)里亞和三角洲濕地水資源韌性水平。國內(nèi)學(xué)者對(duì)水資源韌性研究較少，水資源系統(tǒng)韌性評(píng)價(jià)及機(jī)理分析［17，21-22］、水庫防洪韌性評(píng)價(jià)［23］是其主要研究方向。孫才志等［21］采用SBM-DEA模型及發(fā)展協(xié)調(diào)模型對(duì)全國31個(gè)省區(qū)的水資源效率和韌性進(jìn)行了測(cè)度，并探究了兩者的耦合協(xié)調(diào)關(guān)系。楊亞鋒等［17］從抵抗性、適應(yīng)性、恢復(fù)性角度構(gòu)建了水資源韌性指標(biāo)評(píng)價(jià)體系，使用可變勾股VIKOR決策模型對(duì)長(zhǎng)江經(jīng)濟(jì)帶11個(gè)省級(jí)行政區(qū)水資源韌性水平進(jìn)行了測(cè)度。在指標(biāo)體系的構(gòu)建上，相關(guān)學(xué)者主要基于經(jīng)濟(jì)、社會(huì)、資源、環(huán)境視角，從水資源抵抗性、恢復(fù)性、適應(yīng)性等方面去構(gòu)建指標(biāo)體系。水資源韌性測(cè)算方面通常采用熵權(quán)法、層次分析法、TOPSIS法、智能算法等從指標(biāo)體系整體角度進(jìn)行韌性評(píng)價(jià)。

當(dāng)前水資源韌性研究主要側(cè)重于評(píng)價(jià)分析水資源韌性本身，而對(duì)水資源韌性與效率協(xié)同關(guān)系的研究較少。鑒于此，本文首先運(yùn)用熵權(quán)-TOPSIS賦權(quán)法、三維矢量模型和超效率SBM模型測(cè)度中國30個(gè)省級(jí)行政單位（不包括西藏）的水資源效率及韌性發(fā)展水平，進(jìn)而引入哈肯模型探究?jī)烧叩膮f(xié)同關(guān)系，旨在為中國水資源利用和可持續(xù)發(fā)展政策提供依據(jù)和參考。

1 研究方法

1. 1 熵權(quán)法-TOPSIS賦權(quán)

熵權(quán)法作為一種客觀賦權(quán)法可以排除主觀因素，確保評(píng)價(jià)結(jié)果的客觀性。TOPSIS模型則是一種常用的排序方法，通過評(píng)價(jià)對(duì)象與最優(yōu)解之間的距離程度來確定各個(gè)評(píng)價(jià)對(duì)象的優(yōu)劣。將熵權(quán)法和TOPSIS模型結(jié)合使用，能夠使評(píng)價(jià)結(jié)果更加客觀準(zhǔn)確。因此，本文使用熵權(quán)-TOPSIS方法確定水資源韌性指標(biāo)的權(quán)重，詳細(xì)計(jì)算方法見參考文獻(xiàn)［24］。

1. 2 三維矢量模型

水資源韌性體系中抵抗性、恢復(fù)性和適應(yīng)性指標(biāo)涉及到水資源、水環(huán)境、社會(huì)經(jīng)濟(jì)等方面因素，三者表達(dá)功能和意義各不相同，且彼此有著復(fù)雜的非線性關(guān)系，以往通過子系統(tǒng)加權(quán)求得系統(tǒng)韌性的方法存在明顯的缺陷。參考趙自陽等［22］提出的三維空間矢量模型，本文構(gòu)造了水資源韌性抵抗性-恢復(fù)性-適應(yīng)性三維空間矢量模型（圖1）。

原點(diǎn)O 代表系統(tǒng)韌性最差/最不理想的狀態(tài)。Xt、Yt、Zt分別為恢復(fù)性、適應(yīng)性、抵抗性在t 時(shí)刻的狀態(tài)，P（Xt，Yt，Zt ）是水資源系統(tǒng)在t 時(shí)刻的韌性狀態(tài)，PP'⊥OS。其中，Xt為t 年每個(gè)指標(biāo)權(quán)重與指標(biāo)數(shù)據(jù)歸一化值的乘積之和，Yt、Zt計(jì)算方法與Xt相同。OP的長(zhǎng)度（LOP）是水資源韌性的發(fā)展度，見式（1）；θ 為水資源韌性的協(xié)調(diào)度，是指水資源實(shí)際發(fā)展軌跡與理想軌跡之間的擬合度量，見式（2）；L*OP 是最終水資源韌性指數(shù)，本文定義為韌性發(fā)展度在理想軌跡路徑OS 上的投影，見式（3）。

1. 3 超效率SBM模型

超效率SBM模型將投入產(chǎn)出的松弛變量直接納入目標(biāo)函數(shù)，有效避免了松弛變量對(duì)測(cè)度值的影響，解決了決策單元最大閾值為1的問題，有助于進(jìn)一步比較決策單元，計(jì)算公式參閱文獻(xiàn)［25］。

1. 4 哈肯模型

協(xié)同學(xué)理論的哈肯模型用于描述具有自組織特征的系統(tǒng)演化趨勢(shì)，其通過“絕熱近似”假設(shè)識(shí)別系統(tǒng)內(nèi)部序參量，呈現(xiàn)內(nèi)部要素從無序到有序的演化過程［26］。在水資源效率與韌性系統(tǒng)中，系統(tǒng)內(nèi)部存在自組織性。水資源韌性是效率發(fā)展的基礎(chǔ)條件，較高的水資源韌性通過優(yōu)良的資源稟賦條件、完善的基礎(chǔ)設(shè)施、科學(xué)的節(jié)水政策等提高了水資源利用效率。水資源利用效率的提高可以降低用水成本、減輕對(duì)生態(tài)系統(tǒng)的壓力、減少水資源沖突，從而間接提升水資源韌性水平。基于水資源效率與韌性的自組織特征，本文利用哈肯模型對(duì)系統(tǒng)內(nèi)部演化過程進(jìn)行了分析，具體的模型構(gòu)建如下。

a）提出假設(shè)，構(gòu)建演化方程。本文分別用q1和q2表示水資源韌性和效率的狀態(tài)變量，假設(shè)q1是序參量，2 個(gè)變量的系統(tǒng)演化所產(chǎn)生的影響表示為式（4）、（5）：

q1 = -λ1q1 - aq1q2 （4）

q2 = -γ2q2 + bq21（5）

式中：q1、q2為狀態(tài)變量；γ1、γ2分別為2個(gè)子系統(tǒng)的阻尼系數(shù)；a、b 為q1和q2的相互作用程度。

b）演化方程求解。判斷方程是否滿足“絕熱近似條件”：|γ2|≥|γ2|，且γ2>0。若滿足，則q1為系統(tǒng)發(fā)展的序參量，否則轉(zhuǎn)至①，從而求得序參量。

c）若滿足以上絕熱近似條件，可令q2=0，得到序參量演化方程，見式（6）。對(duì)q1的相反積分求得系統(tǒng)勢(shì)函數(shù)，見式（7）。

d）協(xié)同得分值求解。中國水資源效率與韌性的穩(wěn)定點(diǎn)由式（7）等于0求得［q*，v（q）］，在系統(tǒng)中任意q 點(diǎn)與其到均衡點(diǎn)的距離用來評(píng)價(jià)水資VaSjcANiMsz6j9q2Ci19qQ==源韌性與效率的狀態(tài)，即協(xié)同值d。

由于哈肯模型最初主要應(yīng)用于物理學(xué)領(lǐng)域，物理方程中多為連續(xù)隨機(jī)變量，所以在社會(huì)經(jīng)濟(jì)數(shù)據(jù)中要進(jìn)行離散化處理，得到水資源韌性（WES）與效率（EF）的演化方程，其中t 為時(shí)間變量：

WES（t） = （1 - γ1 ）WES（t - 1） - aWES（t - 1） EF（t - 1） （9）

EF = （1 - γ2 ）EF（t - 1） + bWES2（t - 1） （10）

2 指標(biāo)體系構(gòu)建與數(shù)據(jù)來源

2. 1 水資源效率和韌性指標(biāo)的建立

水資源效率投入指標(biāo)選用用水總量、資本存量和全行業(yè)從業(yè)人員，本文采用“永續(xù)盤存法”估算資本存量，以張軍等［27］計(jì)算的2000年的資本存量為起點(diǎn)。期望產(chǎn)出指標(biāo)選用地區(qū)生產(chǎn)總值GDP，以2000年為基期通過各省區(qū)的GDP指數(shù)進(jìn)行處理，以消除價(jià)格影響。非期望產(chǎn)出指標(biāo)選用灰水足跡和化學(xué)需氧量（COD）排放量，灰水足跡計(jì)算方式參考文獻(xiàn)［28］。投入和產(chǎn)出指標(biāo)的選取均參考相關(guān)文獻(xiàn)［14，29-30］。水資源韌性指標(biāo)分別從抵抗性、恢復(fù)性和適應(yīng)性3個(gè)方面選?。?7，21-22，31］，具體指標(biāo)與等級(jí)劃分見表1。

2. 2 數(shù)據(jù)來源

本文數(shù)據(jù)主要來源于2000—2022年《中國統(tǒng)計(jì)年鑒》《中國環(huán)境統(tǒng)計(jì)年鑒》《中國城鄉(xiāng)建設(shè)統(tǒng)計(jì)年鑒》和各省《水資源公報(bào)》，部分缺失數(shù)據(jù)通過線性插值法補(bǔ)充。

3 結(jié)果分析

3. 1 水資源效率和韌性測(cè)度結(jié)果分析

本文采用DEArun軟件中的超效率SBM模型對(duì)2000—2022年中國30個(gè)省區(qū)（市）的水資源效率進(jìn)行測(cè)算，選用熵權(quán)-TOPSIS法和三維矢量模型測(cè)算了研究區(qū)內(nèi)水資源韌性水平，具體計(jì)算結(jié)果見表2。表2中東部、中部、西部及東北地區(qū)劃分是依據(jù)國家統(tǒng)計(jì)局2021年發(fā)布的標(biāo)準(zhǔn)。根據(jù)計(jì)算結(jié)果繪制了分區(qū)域效率和韌性均值變化趨勢(shì)（圖2）。

從表2可知，中國水資源效率整體上呈緩慢上升趨勢(shì)，均值由0. 26（2000 年）增長(zhǎng)至0. 49（2022年）。2011—2019年效率水平呈明顯上升趨勢(shì)，說明隨著“十二五”和“十三五”規(guī)劃所倡導(dǎo)的綠色、低碳經(jīng)濟(jì)理念，經(jīng)濟(jì)向著高質(zhì)量發(fā)展方向轉(zhuǎn)變，水資源利用由粗放型向集約型轉(zhuǎn)變?nèi)〉昧溯^好效果。2020—2022年整體呈下降趨勢(shì)，原因可能是2020年新冠肺炎疫情的暴發(fā)，經(jīng)濟(jì)生產(chǎn)活動(dòng)的減少導(dǎo)致各地生產(chǎn)總值放緩，從而影響了水資源利用效率。

由圖2a可知，水資源效率水平整體上表現(xiàn)為東部地區(qū)>東北地區(qū)>中部地區(qū)>西部地區(qū)。東部地區(qū)產(chǎn)業(yè)結(jié)構(gòu)多樣化且擁有更多的資金和技術(shù)支持，因此水資源利用效率最高。2011—2022年東北和中部地區(qū)的水資源效率均值與中西部地區(qū)相當(dāng)，造成這一差異的可能原因在于傳統(tǒng)產(chǎn)業(yè)的式微、人口外流以及經(jīng)濟(jì)增長(zhǎng)乏力等因素導(dǎo)致了東北和中部地區(qū)水資源利用效率增長(zhǎng)速度受限，這也顯示了優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)必要性和緊迫性。從區(qū)域內(nèi)部來看，高于全國效率均值水平的有13個(gè)省份，效率均值排前五名的是上海、北京、廣東、浙江、山東，主要分布在東部地區(qū)且均屬于水資源高效區(qū)。效率均值處于0. 3～0. 4的有福建、湖北、湖南、江西、四川、遼寧等省份，主要分布在中部地區(qū)。這些地區(qū)由于受東部地區(qū)經(jīng)濟(jì)、技術(shù)輻射影響，水資源利用效率相對(duì)較高。山西、內(nèi)蒙、河北等地效率水平低于東部和中部地區(qū)，這些省多以高耗水的重工業(yè)如煤炭、鋼鐵為支柱產(chǎn)業(yè)，通常對(duì)水資源的需求量較大，且在生產(chǎn)過程中對(duì)水的回收和再利用率較低，水資源效率不高。水資源效率均值處在末尾的省份有甘肅、青海、寧夏和新疆。這些省份位于西部地區(qū)，區(qū)域內(nèi)因缺水嚴(yán)重、農(nóng)業(yè)用水占比較大但農(nóng)業(yè)所產(chǎn)生的經(jīng)濟(jì)效益不高以及節(jié)水灌溉技術(shù)應(yīng)用不足等原因?qū)е滤Y源效率整體偏低。

表2 表明中國水資源韌性均值由0. 35（2000年）增長(zhǎng)至0. 58（2022年），平均值為0. 45，說明中國水資源韌性水平處于不斷提升的狀態(tài)，但整體處于發(fā)展階段，還有較大的上升空間。從圖2b 可以看出，2010年水資源韌性整體呈下降趨勢(shì)，主要原因是該年份多數(shù)省份發(fā)生旱災(zāi)，降水稀少，導(dǎo)致地缺墑、人缺水。2011—2016年的水資源韌性值由0. 41上升至0. 53，上升幅度相對(duì)較大，這主要得益于“十二五”規(guī)劃的實(shí)施，政府在生態(tài)環(huán)境保護(hù)、水資源管理政策、水利基礎(chǔ)設(shè)施等方面推行了相關(guān)的利好措施，提高了水資源的抵抗性和適應(yīng)性。2017—2022年水資源韌性發(fā)展相對(duì)平穩(wěn)，這可能是該時(shí)期水利基礎(chǔ)設(shè)施趨于完善、節(jié)水措施有效推行、污水處理率顯著提升等原因，從而導(dǎo)致水資源韌性值的年變化率減小。

由圖2b可知，2000—2022年水資源韌性均值由高到低依次為東部地區(qū)、東北地區(qū)、中部地區(qū)和西部地區(qū)。中部地區(qū)和東北地區(qū)的水資源韌性在2016—2018年呈現(xiàn)輕微下降趨勢(shì)，這可能是由于區(qū)域內(nèi)產(chǎn)業(yè)結(jié)構(gòu)多次調(diào)整升級(jí)，水資源的抵抗能力和恢復(fù)能力不穩(wěn)定，導(dǎo)致該時(shí)期水資源韌性水平略有下降。區(qū)域內(nèi)部韌性水平存在梯度差異。水資源韌性均值在0. 55以上的?。ㄖ陛犑校┦钦憬?、上海、廣東、江西和四川，屬于第一梯隊(duì)。均值在0. 45以上且排名前六的省份為江蘇、福建、山東、湖南、湖北和吉林，屬于第二梯隊(duì)。這些省份主要集中在東部和中部地區(qū)，區(qū)域內(nèi)由于降水和地表水資源豐富、水利基礎(chǔ)設(shè)施相對(duì)完善、環(huán)保資金充足等原因，使得水資源恢復(fù)性和適應(yīng)性水平相對(duì)較高。水資源韌性均值在0. 4以下的省（自治區(qū)）有甘肅、青海、寧夏、貴州、新疆，屬于第三梯隊(duì)，這些省份主要集中在西部地區(qū)，由于地表水資源較為匱乏、水資源利用率不高、生態(tài)環(huán)境脆弱等原因?qū)е滤Y源恢復(fù)能力和適應(yīng)能力相對(duì)較差。

3. 2 水資源效率與韌性協(xié)同性分析

基于哈肯模型基本原理，以水資源效率（δW）和水資源韌性（RW）為變量，通過Eviews軟件得到模型方程，結(jié)果見表3。

演化方程為q1=-0. 18q1+0. 671q13，勢(shì)函數(shù)為v=0. 009q12-0. 168q14。令q=0，得出演化方程3個(gè)解q?=0；q?=1. 232；q=-1. 232，由于水資源效率與韌性均為正值，故只考慮q>0的部分，將解為正的值代入勢(shì)函數(shù)得出v=-0. 373，均衡點(diǎn)為（1. 232，-0. 373），在系統(tǒng)中任意q 點(diǎn)與其到均衡點(diǎn)的距離用來評(píng)價(jià)水資源韌性與效率的狀態(tài)。距離公式表示為：d =根號(hào)下（q - 1. 232）2 + （v + 0. 373）2。d 越大代表區(qū)域水資源韌性與效率協(xié)調(diào)程度越低，需要對(duì)其進(jìn)行逆向化處理，最終得到兩者的協(xié)同得分（表4）。

從表4可以得出，中國水資源效率與韌性的協(xié)同度由0. 24（2000年）上升至2022 年的0. 48（2022年），整體協(xié)同水平不高。2000—2010年水資源效率與韌性協(xié)同水平處于平穩(wěn)發(fā)展?fàn)顟B(tài)；2011—2019年兩者協(xié)同水平上升幅度增大。主要由于西部大開發(fā)戰(zhàn)略、東北振興戰(zhàn)略的實(shí)施推動(dòng)了西部和東北地區(qū)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整和升級(jí)，提高了水資源利用效率，加強(qiáng)了水資源的韌性。同時(shí)，南水北調(diào)工程中線和東線工程的竣工，改善了水資源的供需格局，很大程度上解決了中部和東部地區(qū)缺水的問題，促進(jìn)了水資源的均衡配置和高效利用。協(xié)同度在2000—2021年呈下降趨勢(shì)，原因可能是2020年之后經(jīng)濟(jì)發(fā)展處于停滯狀態(tài)，地區(qū)生產(chǎn)總值增長(zhǎng)放緩，政府對(duì)水資源管理和利用的滯后優(yōu)化，從而影響了兩者的協(xié)同水平。

2000—2022 年水資源效率與韌性協(xié)同均值由高到低依次為東部地區(qū)、東北地區(qū)、中部地區(qū)、西部地區(qū)。東部地區(qū)擁有豐富的水資源儲(chǔ)備，水利基礎(chǔ)設(shè)施較為完善，經(jīng)濟(jì)發(fā)展迅速，因此在水資源韌性與效率協(xié)作水平方面最高。西部地區(qū)多以農(nóng)業(yè)和能源產(chǎn)業(yè)為主，水資源需求量大且多數(shù)省份缺水嚴(yán)重，導(dǎo)致水資源供需矛盾突出，水資源效率與韌性的協(xié)同水平最低。區(qū)域內(nèi)部協(xié)同水平同樣存在梯度差異。上海（0. 64）、福建（0. 59）、浙江（0. 57）、廣東（0. 59）、山東（0. 58）等省份協(xié)同水平屬于第一梯隊(duì)，主要分布在東部地區(qū)。這些地區(qū)經(jīng)濟(jì)發(fā)展水平及水利基礎(chǔ)設(shè)施建設(shè)能夠支持高密度人口的用水需求，水資源效率和韌性發(fā)展水平均高。四川（0. 55）、遼寧（0. 53）、湖北（0. 51）、江西（0. 49）、云南（0. 48）等省份屬于第二梯隊(duì)，主要分布在中西部偏南地區(qū)，這些地區(qū)受東部地區(qū)經(jīng)濟(jì)、技術(shù)輻射作用明顯，產(chǎn)業(yè)集聚性更強(qiáng)，對(duì)資源的利用率更高，因此協(xié)同水平僅次于東部地區(qū)。河南（0. 45）、河北（0. 43）、內(nèi)蒙（0. 42）、山西（0. 41）、黑龍江（0. 41）處于中間水平，這些省份多以資源型城市為主，工業(yè)耗水量大、水資源回收利用率低、水環(huán)境污染嚴(yán)重等問題突出，致使水資源效率與韌性的協(xié)同關(guān)系受損。陜西（0. 4）、甘肅（0. 39）、青海（0. 38）、寧夏（0. 38）等省份屬于第四梯隊(duì)，多分布于西北地區(qū)。這些省份農(nóng)業(yè)用水量大但所產(chǎn)生的經(jīng)濟(jì)效益不高，水資源效率水平較低，水資源恢復(fù)能力和適應(yīng)能力較差，因此整體效率與韌性協(xié)同度不高。

4 結(jié)論與建議

4. 1 結(jié)論

a）研究期間中國水資源效率整體上呈波動(dòng)式上升（2000—2019年）—緩慢下降（2020—2022年）趨勢(shì)，2009—2019年上升速度較快。全國30個(gè)省份的水資源均值均小于1，整體上處于非有效區(qū)。

b）研究期間中國水資源韌性水平呈平穩(wěn)波動(dòng)（2000—2010年）—快速上升（2011—2016年）—緩慢上升（2017—2022年）趨勢(shì)，韌性均值為0. 45，整體上處于較低水平，還有較大上升空間。

c）研究期間中國水資源效率與韌性協(xié)同度呈緩慢上升趨勢(shì)，由0. 24（2000年）上升至0. 48（2022年），整體上處于中度協(xié)同水平。水資源效率作為協(xié)同演化的序參量，對(duì)整個(gè)演化過程起積極導(dǎo)向作用。

4. 2 建議

a）西部地區(qū)水資源相對(duì)短缺，因此要積極優(yōu)化用水方式，加強(qiáng)對(duì)節(jié)水灌溉技術(shù)的研發(fā)和應(yīng)用。中部和東北地區(qū)要積極發(fā)展低耗水、節(jié)水產(chǎn)業(yè)，加大技術(shù)研發(fā)投入，提升區(qū)域自主創(chuàng)新能力，促進(jìn)新一輪產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級(jí)。

b）各地要嚴(yán)格按照“雙控”行動(dòng)的要求，合理制定可監(jiān)測(cè)、可量化的節(jié)水指標(biāo)，從源頭上減少水資源使用量。政府可以建立產(chǎn)業(yè)節(jié)水型企業(yè)標(biāo)桿，使其在整個(gè)產(chǎn)業(yè)鏈上形成集聚輻射效應(yīng)，促進(jìn)上下游產(chǎn)業(yè)完善用水規(guī)制，整體提升區(qū)域水資源用水效率水平。

c）各地要加快推動(dòng)跨區(qū)域合作，促進(jìn)不同地區(qū)間水資源信息交流、技術(shù)合作和經(jīng)驗(yàn)共享，優(yōu)化資源共享模式和方案，提升水資源分配方面的銜接度和協(xié)調(diào)度，實(shí)現(xiàn)水資源效率與韌性的共同發(fā)展。

參考文獻(xiàn)：

［1］ 鄭德鳳，王明，安志英. 水-能源-糧食紐帶系統(tǒng)研究脈絡(luò)梳理與分析［J］. 水利經(jīng)濟(jì)，2022，40（6）：8-15，102.

［2］ 劉浩然，黨素珍，連炎清. 黃河流域國家中心城市水資源利用效率分析［J］. 人民黃河，2023，45（11）：86-90，94.

［3］ 張明斗，翁愛華. 長(zhǎng)江經(jīng)濟(jì)帶城市水資源利用效率的空間關(guān)聯(lián)網(wǎng)絡(luò)及形成機(jī)制［J］. 地理學(xué)報(bào)，2022，77（9）： 2353-2573.

［4］ 江激宇，劉嘉銘，張士云. 淮河流域谷物生產(chǎn)水足跡及用水效率研究［J］. 農(nóng)業(yè)資源與環(huán)境學(xué)報(bào)，2024，41（2）：371-382.

［5］ 黨麗娟. 黃河流域工業(yè)用水效率提升策略研究［J］. 人民黃河，2023，45（1）：76-81.

［6］ 陳杰，許朗. 基于面板三階段DEA-Malmquist模型的中國農(nóng)業(yè)綠色水資源利用效率研究［J］. 地理科學(xué)，2023，43（4）：709-718.

［7］ 張義敏，蔡旭東，陳雅芬，等. 基于DEA的水資源利用效率評(píng)價(jià)及節(jié)水潛力分析［J］. 人民珠江，2015，36（6）： 121-124.

［8］ 劉維哲，常明，王西琴. 基于隨機(jī)前沿的灌溉用水效率及影響因素研究：以陜西關(guān)中地區(qū)小麥為例［J］. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2018，26（9）：1407-1414.

［9］ 任玉芬，方文穎，歐陽志云，等. 基于面板數(shù)據(jù)的我國城市水資源水環(huán)境隨機(jī)前沿面分析［J］. 環(huán)境科學(xué)學(xué)報(bào)，2020，40（7）：2638-2643.

［10］ 高媛媛，許新宜，王紅瑞，等. 中國水資源利用效率評(píng)估模型構(gòu)建及應(yīng)用［J］. 系統(tǒng)工程理論與實(shí)踐，2013，33（3）： 776-784.

［11］ 劉家宏，邵薇薇，王浩，等. 水資源利用效率頻譜分析方法及應(yīng)用［J］. 水利水運(yùn)工程學(xué)報(bào)，2019（6）：132-138.

［12］ 宋國君，何偉. 中國城市水資源利用效率標(biāo)桿研究［J］. 資源科學(xué)，2014，36（12）：2569-2577.

［13］ 唐明，周涵杰，許文濤，等. 區(qū)域用水效率綜合評(píng)價(jià)：新方法研究及其應(yīng)用［J］. 節(jié)水灌溉，2022（5）：89-96.

［14］ 趙良仕，孫才志，鄭德鳳. 中國省際水資源利用效率與空間溢出效應(yīng)測(cè)度［J］. 地理學(xué)報(bào)，2014，69（1）：121-133.

［15］ 趙良仕，冷明祥，孫才志. 中國省際全要素水資源綠色效率及節(jié)水減排潛力分析［J］. 地域研究與開發(fā)，2021，40（4）：123-128.

［16］ 郭炳南，馮雨，張浩. 長(zhǎng)江經(jīng)濟(jì)帶工業(yè)綠色全要素水資源效率的區(qū)域差異與隨機(jī)收斂［J］. 生態(tài)經(jīng)濟(jì)，2023，39（2）：165-172.

［17］ 楊亞鋒，王紅瑞，鞏書鑫，等. 可變勾股模糊VIKOR水資源系統(tǒng)韌性評(píng)價(jià)調(diào)控模型及應(yīng)用［J］. 水利學(xué)報(bào)，2021，52（6）：633-646.

［18］ SANDOVAL-SOLIS S， MCKINNEY D C， LOUCKS D P.Sustainability index for water resources planning and management［J］. Journal of Water Resources Planning and Management，2011，137（5）：381-390.

［19］ SCHLVTER M， KHASANKHANOVA G， TALSKIKH V， et al.Enhancing resilience to water flow uncertainty by integrating environmental flows into water management in the Amudarya River， Central Asia［J］. Global and Planetary Change，2013，110：114-129.

［20］ SWEYA L N， WILKINSON S J E I. A tool for measuring environmental resilience to floods in Tanzania water supply systems ［J］. Ecological Indicators，2020，112. DOI：10. 1016/j. ecolind.2020. 106165.

［21］ 孫才志，孟程程. 中國區(qū)域水資源系統(tǒng)韌性與效率的發(fā)展協(xié)調(diào)關(guān)系評(píng)價(jià)［J］. 地理科學(xué)，2020，40（12）：2094-2104.

［22］ 趙自陽，王紅瑞，張力，等. 長(zhǎng)江經(jīng)濟(jì)帶水資源-水環(huán)境-社會(huì)經(jīng)濟(jì)復(fù)雜系統(tǒng)韌性調(diào)控模型及應(yīng)用［J］. 水科學(xué)進(jìn)展，2022，33（5）：705-717.

［23］ 張曉琦，劉攀，陳進(jìn)，等. 基于條件風(fēng)險(xiǎn)價(jià)值理論的水庫群防洪庫容協(xié)同作用［J］. 水科學(xué)進(jìn)展，2022，33（2）： 298-305.

［24］ 王松茂，牛金蘭. 黃河流域城市生態(tài)韌性時(shí)空演變及其影響因素［J］. 生態(tài)學(xué)報(bào)，2023，43（20）：8309-8320.

［25］ 張倩歌，楊茂. 基于超效率SBM模型的黃河水資源利用效率評(píng)價(jià)［J］. 人民黃河，2022，44（9）：111-115.

［26］ 王松茂，牛金蘭. 山東省旅游經(jīng)濟(jì)與城市生態(tài)韌性協(xié)同演化研究［J］. 地理學(xué)報(bào)，2023，78（10）：2591-2608.

［27］ 張軍，吳桂英，張吉鵬. 中國省際物質(zhì)資本存量估算：1952—2000［J］. 經(jīng)濟(jì)研究，2004（10）：35-44.

［28］ 韓琴，孫才志，鄒瑋. 1998-2012年中國省際灰水足跡效率測(cè)度與驅(qū)動(dòng)模式分析［J］. 資源科學(xué)，2016，38（6）：1179-1191.

［29］ 田澤，方琪，陳一飛. 我國三大流域水資源利用效率時(shí)空演變特征與影響因素分析［J］. 水利經(jīng)濟(jì)，2023，41（6）：4-11，110.

［30］ 左其亭，張志卓，馬軍霞. 黃河流域水資源利用水平與經(jīng)濟(jì)社會(huì)發(fā)展的關(guān)系［J］. 中國人口·資源與環(huán)境，2021，31（10）：29-38.

［31］ 孫宇，劉維忠，盛洋. 基于PSR模型的新疆水資源經(jīng)濟(jì)生態(tài)韌性時(shí)空差異及影響因素分析［J］. 干旱區(qū)地理，2023，46（12）：2017-2028.

（責(zé)任編輯：程茜）

基金項(xiàng)目：遼寧省經(jīng)濟(jì)社會(huì)發(fā)展研究課題（2023lslybkt-080）；教育部人文社會(huì)科學(xué)重點(diǎn)研究基地“海洋可持續(xù)發(fā)展”學(xué)科交叉專向研究課題（XKJC2023004）