陳 立，于 偉

（曲阜師范大學(xué)地理與旅游學(xué)院，山東 日照 276800）

中國(guó)區(qū)域經(jīng)濟(jì)的長(zhǎng)期快速增長(zhǎng)帶來(lái)了水資源供給不足與利用粗放并存、水環(huán)境污染、水生態(tài)失衡等一系列問(wèn)題[1]。在加快推進(jìn)經(jīng)濟(jì)高質(zhì)量發(fā)展時(shí)代背景下，立足人多水少的基本國(guó)情，協(xié)調(diào)好水資源利用效率和經(jīng)濟(jì)發(fā)展水平之間的關(guān)系成為亟待解決的重要現(xiàn)問(wèn)題。山東省作為沿海經(jīng)濟(jì)大省，地處黃河流域下游的東部季風(fēng)區(qū)，降水年際變化大，是全國(guó)旱災(zāi)頻發(fā)地之一，全省人均水資源量?jī)H為全國(guó)平均水平的1/7，2/3 的城市屬于缺水型城市。水資源短缺嚴(yán)重制約了山東省經(jīng)濟(jì)的綠色轉(zhuǎn)型，作為黃河流域的龍頭省份，如何協(xié)調(diào)水資源硬約束與區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展的關(guān)系是亟待研究的科學(xué)問(wèn)題[2]。

伴隨發(fā)展中國(guó)家城鎮(zhèn)化的快速推進(jìn)，水資源與經(jīng)濟(jì)協(xié)調(diào)發(fā)展成為國(guó)際社會(huì)和國(guó)內(nèi)外學(xué)者高度關(guān)注的問(wèn)題[3]。1990 年，聯(lián)合國(guó)發(fā)布的《水與可持續(xù)發(fā)展準(zhǔn)則：原理與政策方案》論述了水資源開(kāi)發(fā)與經(jīng)濟(jì)發(fā)展之間的緊密關(guān)系。1992 年，《里約環(huán)境與發(fā)展宣言》首次從可持續(xù)視角探討了水資源與經(jīng)濟(jì)協(xié)調(diào)發(fā)展的問(wèn)題。國(guó)內(nèi)外學(xué)者對(duì)水資源和經(jīng)濟(jì)系統(tǒng)之間關(guān)系的研究議題主要集中在水資源與經(jīng)濟(jì)社會(huì)發(fā)展要素的時(shí)空匹配[4-6]、水資源與經(jīng)濟(jì)發(fā)展的耦合協(xié)調(diào)與仿真模擬[7-15]、水資源利用與經(jīng)濟(jì)發(fā)展的脫鉤分析[16-19]、水資源和經(jīng)濟(jì)協(xié)同發(fā)展的路徑[20-21]等方面。近年來(lái)，隨著黃河流域生態(tài)保護(hù)與高質(zhì)量發(fā)展上升為國(guó)家戰(zhàn)略，水資源和區(qū)域經(jīng)濟(jì)的耦合協(xié)調(diào)度評(píng)價(jià)成為研究熱點(diǎn)[2，21-27]。國(guó)內(nèi)外學(xué)者針對(duì)水資源和區(qū)域經(jīng)濟(jì)關(guān)系開(kāi)展了豐富的案例研究，取得了大量有益的成果，但仍存在以下不足：第一，已有案例研究多采用單一的耦合協(xié)調(diào)度模型衡量二者的關(guān)系，難以反映水資源與區(qū)域經(jīng)濟(jì)系統(tǒng)耦合協(xié)調(diào)的動(dòng)態(tài)變化。第二，學(xué)者大多關(guān)注省域、流域等大尺度地域單元，對(duì)縣域等中小尺度區(qū)域關(guān)注不足，縣域作為我國(guó)經(jīng)濟(jì)發(fā)展的重要載體，水資源利用效率不高問(wèn)題往往比較典型，值得深入研究。第三，以往研究多從區(qū)域整體層面分析動(dòng)態(tài)耦合協(xié)調(diào)情況，忽略了區(qū)域經(jīng)濟(jì)發(fā)展和水資源稟賦的差異，難以刻畫(huà)區(qū)域內(nèi)部的動(dòng)態(tài)耦合情況及其演化歷程。

鑒于此，本文基于生態(tài)脆弱區(qū)黃河三角洲的分縣數(shù)據(jù)，綜合利用耦合協(xié)調(diào)度模型和剪刀差方法探討水資源和區(qū)域經(jīng)濟(jì)系統(tǒng)之間的耦合協(xié)調(diào)情況及其動(dòng)態(tài)變化，為黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展國(guó)家戰(zhàn)略的順利施提供科學(xué)依據(jù)。

1 研究區(qū)域概況

黃河三角洲是黃河攜帶大量泥沙在渤海凹陷處沉積形成的沖積平原。2009 年，國(guó)務(wù)院通過(guò)《黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃》，標(biāo)志著該區(qū)域的保護(hù)開(kāi)發(fā)上升為國(guó)家戰(zhàn)略。依據(jù)該規(guī)劃確定黃河三角洲的研究范圍，包括東營(yíng)、濱州兩個(gè)地級(jí)市所轄全部縣域單元，德州慶云縣和樂(lè)陵市，濰坊壽光市、昌邑市和寒亭區(qū)，煙臺(tái)萊州市和淄博高青縣，共計(jì)19 個(gè)縣域單元。

2021 年，黃河三角洲GDP 合計(jì)7752.47 億元，占山東省GDP 的9.32%；面積2.65×104km2，占山東省陸域面積的16.78%。人均用水量約為395 m3，遠(yuǎn)高于山東省206 m3的平均水平?？梢?jiàn)，黃河三角洲的經(jīng)濟(jì)發(fā)展水平不高，但水資源消耗較高。黃河三角洲生態(tài)環(huán)境較為脆弱，降水不均、旱災(zāi)頻發(fā)、用水短缺現(xiàn)象比較嚴(yán)重，水資源已成為該區(qū)域經(jīng)濟(jì)發(fā)展的剛性約束。以生態(tài)脆弱區(qū)黃河三角洲為研究對(duì)象，探究該區(qū)域水資源和區(qū)域經(jīng)濟(jì)的協(xié)調(diào)發(fā)展動(dòng)態(tài)演化及其影響因素，對(duì)于指引黃河流域的生態(tài)保護(hù)和高質(zhì)量發(fā)展具有參考價(jià)值。

2 指標(biāo)體系與研究方法

2.1 指標(biāo)體系

結(jié)合黃河三角洲的特點(diǎn)，初步選取水資源利用水平與區(qū)域經(jīng)濟(jì)發(fā)展水平的評(píng)價(jià)指標(biāo)體系，在征詢專家意見(jiàn)的基礎(chǔ)上進(jìn)行修改，形成本文的評(píng)價(jià)指標(biāo)，采用熵值法確定指標(biāo)權(quán)重（表1）。水資源系統(tǒng)包括用水負(fù)荷、利用效率、用水結(jié)構(gòu)3 個(gè)層面的7 個(gè)指標(biāo)；區(qū)域經(jīng)濟(jì)系統(tǒng)包括發(fā)展規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)、居民生活和生態(tài)效應(yīng)4 個(gè)層面的10 個(gè)指標(biāo)。其中，水結(jié)構(gòu)均衡度參考易晶晶等[28]設(shè)計(jì)的方法進(jìn)行計(jì)算，反映用水結(jié)構(gòu)的均衡程度，值越大表明單一用水類別優(yōu)勢(shì)度越小，用水系統(tǒng)越趨于良性。水資源供給模數(shù)指標(biāo)是總供水量與土地面積的比值。產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化指數(shù)參考付凌暉[29]的計(jì)算方式，按照三次產(chǎn)業(yè)的比重進(jìn)行計(jì)算。

表1 水資源利用與經(jīng)濟(jì)發(fā)展評(píng)價(jià)指標(biāo)體系

2.2 研究方法

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

耦合協(xié)調(diào)模型指兩個(gè)或者多個(gè)系統(tǒng)之間的相互影響和協(xié)調(diào)發(fā)展的程度，反映水資源利用和經(jīng)濟(jì)系統(tǒng)之間的協(xié)調(diào)狀態(tài)。計(jì)算公式如下：

式中，D為耦合協(xié)調(diào)度，C為耦合關(guān)聯(lián)度，T為綜合評(píng)價(jià)指數(shù)，f（x）為區(qū)域經(jīng)濟(jì)發(fā)展水平評(píng)價(jià)值，f（y）為水資源系統(tǒng)評(píng)價(jià)值，α、β分別為經(jīng)濟(jì)系統(tǒng)和水資源利用的貢獻(xiàn)度，本文認(rèn)為二者同等重要，即α=β=0.5。參照已有研究成果[24]，黃河三角洲水資源利用和經(jīng)濟(jì)系統(tǒng)之間耦合協(xié)調(diào)發(fā)展階段及類型判定見(jiàn)表2。

表2 耦合協(xié)調(diào)度與耦合等級(jí)

2.2.2 相對(duì)發(fā)展度模型

耦合協(xié)調(diào)模型反映區(qū)域經(jīng)濟(jì)和水資源之間的協(xié)調(diào)發(fā)展程度，但無(wú)法度量二者之間的發(fā)展差距。因此，引入相對(duì)發(fā)展度模型體現(xiàn)二者之間的超前或滯后發(fā)展程度。公式如下：

式中，P為相對(duì)發(fā)展度。根據(jù)文獻(xiàn)[22]，當(dāng)P＞1.2 時(shí)，為區(qū)域經(jīng)濟(jì)發(fā)展超前型；當(dāng)0.8＜P≤1.2，為區(qū)域經(jīng)濟(jì)發(fā)展和水資源同步發(fā)展型；當(dāng)P≤0.8 時(shí)，為區(qū)域經(jīng)濟(jì)發(fā)展滯后型。

2.2.3 剪刀差方法

考慮水資源系統(tǒng)和區(qū)域經(jīng)濟(jì)系統(tǒng)之間的關(guān)系處于動(dòng)態(tài)變化中，利用剪刀差方法[23]分析黃河三角洲兩個(gè)系統(tǒng)耦合關(guān)系的動(dòng)態(tài)演化。計(jì)算公式如下：

式中，F(xiàn)′（x）代表區(qū)域經(jīng)濟(jì)系統(tǒng)的變化趨勢(shì)，F(xiàn)′（y）代表水資源系統(tǒng)的變化趨勢(shì)，α是給定時(shí)刻兩系統(tǒng)變化曲線的夾角值，α越大說(shuō)明兩系統(tǒng)變化趨勢(shì)的差異越大。

2.3 數(shù)據(jù)來(lái)源

水資源數(shù)據(jù)主要來(lái)源于2011—2020 年各地市的水資源公報(bào)，部分公報(bào)中沒(méi)有水資源數(shù)據(jù)的地區(qū)，在人民政府官方網(wǎng)站“政民互動(dòng)”板塊發(fā)起數(shù)據(jù)收集請(qǐng)求，水利局通過(guò)電子郵件回復(fù)并發(fā)送了相關(guān)數(shù)據(jù)。經(jīng)濟(jì)數(shù)據(jù)來(lái)源于2012—2021 年?yáng)|營(yíng)、濱州等相關(guān)地市的統(tǒng)計(jì)年鑒，以及各縣域單元相應(yīng)的國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。PM2.5利用美國(guó)哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)與應(yīng)用中心提供的全球PM2.5的年均濃度數(shù)據(jù)整理計(jì)算得到（數(shù)據(jù)獲取網(wǎng)址：

https://sedac.ciesin.columbia.edu/data/set/sdei-globalannual-gwr-pm2-5-modis-misr-seawifs-aod）。各縣域單元2011—2020 年的碳排放數(shù)據(jù)在Chen et al.[30]基礎(chǔ)上通過(guò)遙感影像數(shù)據(jù)反演得出。

3 結(jié)果分析

3.1 區(qū)域整體的耦合協(xié)調(diào)度及其動(dòng)態(tài)變化

2011—2020 年，黃河三角洲水資源綜合利用水平和經(jīng)濟(jì)發(fā)展水平均呈現(xiàn)動(dòng)態(tài)提升的態(tài)勢(shì)。水資源系統(tǒng)中，用水負(fù)荷所占權(quán)重最大，說(shuō)明水資源利用水平主要受制于用水負(fù)荷，尤其是農(nóng)業(yè)用水量、工業(yè)用水量的增加會(huì)加劇水資源系統(tǒng)的壓力，導(dǎo)致水資源系統(tǒng)綜合指數(shù)的下降。區(qū)域經(jīng)濟(jì)系統(tǒng)中，居民生活和經(jīng)濟(jì)發(fā)展規(guī)模的權(quán)重較大，在一定程度上說(shuō)明區(qū)域經(jīng)濟(jì)發(fā)展水平主要受居民生活質(zhì)量和經(jīng)濟(jì)發(fā)展規(guī)模的影響。

耦合協(xié)調(diào)度D計(jì)算結(jié)果見(jiàn)表3。2011—2020 年，黃河三角洲水資源-經(jīng)濟(jì)系統(tǒng)的協(xié)調(diào)狀態(tài)經(jīng)歷了勉強(qiáng)協(xié)調(diào)和初級(jí)協(xié)調(diào)兩個(gè)階段，系統(tǒng)的耦合類型由低水平磨合向高水平協(xié)調(diào)狀態(tài)過(guò)渡。具體可分為兩個(gè)階段，2011—2017 年，黃河三角洲經(jīng)濟(jì)發(fā)展水平明顯提高，與水資源利用水平評(píng)估值的差距逐步縮小，推動(dòng)了二者協(xié)調(diào)水平的提升，但總體的協(xié)調(diào)狀態(tài)仍為勉強(qiáng)協(xié)調(diào)。2018 年之后，由于經(jīng)濟(jì)進(jìn)一步優(yōu)化發(fā)展提升，與水資源利用系統(tǒng)的耦合協(xié)調(diào)度穩(wěn)定在0.6 以上，協(xié)調(diào)狀態(tài)演化至初級(jí)協(xié)調(diào)階段。2017 年之前，經(jīng)濟(jì)發(fā)展滯后制約了黃河三角洲水資源-經(jīng)濟(jì)系統(tǒng)耦合協(xié)調(diào)狀態(tài)的提升，之后，經(jīng)濟(jì)發(fā)展水平與水資源利用水平差距縮小，系統(tǒng)向初級(jí)協(xié)調(diào)狀態(tài)演化。

表3 黃河三角洲水資源-經(jīng)濟(jì)系統(tǒng)耦合協(xié)調(diào)演化趨勢(shì)

為分析黃河三角洲水資源-經(jīng)濟(jì)系統(tǒng)的動(dòng)態(tài)演化，首先對(duì)系統(tǒng)的綜合評(píng)價(jià)值進(jìn)行曲線擬合，利用擬合方程進(jìn)行求導(dǎo)，得出系統(tǒng)基于年份的變化速率，據(jù)此計(jì)算演化速率的剪刀差值（α），揭示兩個(gè)系統(tǒng)的動(dòng)態(tài)耦合特征（圖1）。2011—2020 年，黃河三角洲水資源-經(jīng)濟(jì)系統(tǒng)演化速率的剪刀差呈現(xiàn)先上升后下降再上升的態(tài)勢(shì)。2011—2014 年，剪刀差由0.046°升至0.728°，該時(shí)段水資源發(fā)展綜合指數(shù)和經(jīng)濟(jì)發(fā)展指數(shù)均呈穩(wěn)步上升態(tài)勢(shì)，但經(jīng)濟(jì)系統(tǒng)上升速率相對(duì)較快，反映了黃河三角洲水資源和經(jīng)濟(jì)系統(tǒng)之間發(fā)展水平差異較大，互相之間無(wú)明顯的限制。2015—2017 年，剪刀差值逐步降至0.046°，表明水資源-經(jīng)濟(jì)系統(tǒng)之間的差距不斷縮小。該時(shí)期水資源系統(tǒng)發(fā)展水平相對(duì)較高，但增長(zhǎng)速率明顯下降；區(qū)域經(jīng)濟(jì)發(fā)展質(zhì)量相對(duì)不高，但變化速率呈現(xiàn)增長(zhǎng)趨勢(shì)，說(shuō)明該時(shí)期區(qū)域經(jīng)濟(jì)對(duì)水資源的壓力逐步增大，正在以損害水資源系統(tǒng)為代價(jià)推動(dòng)經(jīng)濟(jì)增長(zhǎng)。2017 年以來(lái)，剪刀差值迅速增大至2.366°，表明系統(tǒng)的演化速率差異明顯。該時(shí)段水資源利用和經(jīng)濟(jì)發(fā)展均經(jīng)歷了回調(diào)后的快速提升，水資源回升幅度稍大，說(shuō)明水資源支撐經(jīng)濟(jì)發(fā)展的能力有所增強(qiáng)，水資源和區(qū)域經(jīng)濟(jì)開(kāi)始走向協(xié)調(diào)發(fā)展。

圖1 黃河三角洲水資源-經(jīng)濟(jì)系統(tǒng)的剪刀差

3.2 黃河三角洲水資源和區(qū)域經(jīng)濟(jì)分區(qū)縣動(dòng)態(tài)耦合分析

3.2.1 黃河三角洲各縣域耦合協(xié)調(diào)度的時(shí)序變化

研究期內(nèi)各縣域耦合協(xié)調(diào)度均有不同程度的提升，但區(qū)域差異明顯（圖2），東部沿萊州灣的地區(qū)耦合協(xié)調(diào)等級(jí)明顯高于西部地區(qū)，如東營(yíng)區(qū)、壽光市、濱城區(qū)等地。2011 年，黃河三角洲大部分縣域處于瀕臨失調(diào)和勉強(qiáng)協(xié)調(diào)等級(jí)，陽(yáng)信縣、惠民縣、高青縣處于極度失調(diào)階段，主要原因是以上地區(qū)的經(jīng)濟(jì)發(fā)展水平偏低，影響了系統(tǒng)的協(xié)調(diào)發(fā)展。相比之下，東營(yíng)區(qū)、壽光市、萊州市處于初級(jí)協(xié)調(diào)階段，這些地區(qū)經(jīng)濟(jì)發(fā)展質(zhì)量較好，水資源利用效率也較高，因此具有較好的耦合協(xié)調(diào)狀況。至2014年，除鄒平市、高青縣外的縣域耦合協(xié)調(diào)度均有所上升，原因是2012 年山東省政府貫徹落《行最嚴(yán)格水資源管理制度的意見(jiàn)》，推行《山東省用水總量控制管理辦法》，引導(dǎo)了各縣域用水方式和經(jīng)濟(jì)發(fā)展方式的轉(zhuǎn)變，促進(jìn)了系統(tǒng)協(xié)調(diào)發(fā)展，這種上升態(tài)勢(shì)延續(xù)至2017 年。之后，各縣域發(fā)展分異顯著，東營(yíng)區(qū)、濱城區(qū)、壽光市、陽(yáng)信縣等地的耦合協(xié)調(diào)度繼續(xù)大幅提升，墾利區(qū)、利津縣、廣饒縣等地出現(xiàn)了小幅度的下降。以墾利區(qū)為例，當(dāng)?shù)伛詈蠀f(xié)調(diào)度等級(jí)下降的直接原因是人均生活用水量下降明顯，導(dǎo)致水資源利用效率回調(diào)，與經(jīng)濟(jì)發(fā)展水平的差距拉大。

圖2 黃河三角洲各縣域耦合協(xié)調(diào)度的變化

3.2.2 縣域單元水資源與區(qū)域經(jīng)濟(jì)的動(dòng)態(tài)耦合

借助MATLAB（2021）計(jì)算研究區(qū)水資源系統(tǒng)和區(qū)域經(jīng)濟(jì)系統(tǒng)的剪刀差值（圖3）。根據(jù)剪刀差值的變化特征，可將各縣域動(dòng)態(tài)耦合情況分為以下4類。

圖3 黃河三角洲各縣域水資源-經(jīng)濟(jì)系統(tǒng)動(dòng)態(tài)耦合

（1）剪刀差呈現(xiàn)先減后增的“V 型”特征，表明兩系統(tǒng)的演化速率差異先逐漸縮小，后逐步增大，東營(yíng)區(qū)、墾利區(qū)、濱城區(qū)、惠民縣、寒亭區(qū)、壽光市、昌邑市、樂(lè)陵市、慶云縣屬于該類。以壽光市為例，2011—2016 年，壽光市經(jīng)濟(jì)快速增長(zhǎng)，帶來(lái)了水資源消耗、水環(huán)境污染等問(wèn)題，經(jīng)濟(jì)增長(zhǎng)對(duì)水資源利用造成了一定的壓力，限制了水資源-區(qū)域經(jīng)濟(jì)系統(tǒng)的協(xié)調(diào)發(fā)展。2016 年年底，南水北調(diào)壽光段陸續(xù)建設(shè)完成，有效緩解了水資源供給不足的問(wèn)題。2017 年，壽光市開(kāi)展節(jié)水農(nóng)業(yè)技術(shù)改造，提高了水資源利用效率。與此同時(shí)，壽光市全面落2016 年12 月中共中央辦公廳、國(guó)務(wù)院辦公廳印發(fā)《關(guān)于全面落河長(zhǎng)制的意見(jiàn)》，水污染現(xiàn)象得到治理。因此，2017 年以來(lái)，水資源指數(shù)由負(fù)轉(zhuǎn)正，且增速不斷增大，而經(jīng)濟(jì)系統(tǒng)逐步轉(zhuǎn)入高質(zhì)量發(fā)展階段，增長(zhǎng)速率逐漸下降，引起2017—2020 年兩系統(tǒng)剪刀差值迅速上升，演化速率差異增大。

（2）剪刀差呈現(xiàn)“雙減”“雙增”的“W 型”特征，表明兩系統(tǒng)之間的演化速率差異先逐步縮小，后逐步增大，而后重復(fù)該過(guò)程，河口區(qū)、廣饒縣、鄒平市、無(wú)棣縣、陽(yáng)信縣、博興縣、萊州市、高青縣屬于該類。以萊州市為例，2011—2012 年，萊州市經(jīng)濟(jì)發(fā)展水平出現(xiàn)波動(dòng)，水資源指數(shù)增長(zhǎng)速率減緩，經(jīng)濟(jì)系統(tǒng)開(kāi)始對(duì)水資源系統(tǒng)產(chǎn)生負(fù)面影響，兩系統(tǒng)演化速率差異減小。2013—2015 年，經(jīng)濟(jì)指數(shù)持續(xù)上升，水資源指數(shù)持續(xù)下降，以水環(huán)境質(zhì)量降低為代價(jià)換取經(jīng)濟(jì)發(fā)展，兩系統(tǒng)演變速率差異擴(kuò)大。2016—2018 年，萊州市政府完善了農(nóng)田水利基礎(chǔ)設(shè)施、建設(shè)了水系連通補(bǔ)水工程，水資源指數(shù)下降速率減緩，兩系統(tǒng)演化速率差異又逐步減小。2019—2020 年，隨著萊州市政府繼續(xù)加強(qiáng)水資源保護(hù)，推進(jìn)重點(diǎn)水利工程建設(shè)，水資源指數(shù)由下降轉(zhuǎn)化為上升。此外，由于萊州市政府去產(chǎn)能、加快新舊動(dòng)能轉(zhuǎn)換，經(jīng)濟(jì)系統(tǒng)出現(xiàn)一定程度的下降，兩系統(tǒng)演化速率差異增大。

（3）剪刀差呈現(xiàn)“N 型”的波動(dòng)上升特征，表明兩系統(tǒng)之間的演化速率差異先變大后變小，而后演化速率差異又逐步變大，沾化區(qū)屬于此類。2011—2016 年，沾化區(qū)水資源指數(shù)呈現(xiàn)先降后再降的周期性波動(dòng)，區(qū)域經(jīng)濟(jì)系統(tǒng)先增后減。說(shuō)明該時(shí)期水資源和區(qū)域經(jīng)濟(jì)相互影響較小，互相之間限制較小。2017 年，沾化區(qū)推進(jìn)企業(yè)轉(zhuǎn)型升級(jí)，經(jīng)濟(jì)指數(shù)不斷提升但增速下降，同時(shí)該區(qū)加快推進(jìn)重點(diǎn)水利工程建設(shè)，水資源系統(tǒng)持續(xù)優(yōu)化且水資源指數(shù)增速提升，因此兩系統(tǒng)演化速率差異增大。

（4）剪刀差呈現(xiàn)先增后減的“倒U 型”特征，表明兩系統(tǒng)之間的演變速率差異先逐步增大至峰值，而后又逐步減小，利津縣屬于此類。2011—2017 年，當(dāng)?shù)亟?jīng)濟(jì)增長(zhǎng)速率持續(xù)下降但仍保持正增長(zhǎng)，水資源指數(shù)持續(xù)負(fù)增長(zhǎng)，說(shuō)明區(qū)域經(jīng)濟(jì)增長(zhǎng)壓力持續(xù)增加，正以破壞水資源支撐能力為代價(jià)來(lái)發(fā)展經(jīng)濟(jì)，兩系統(tǒng)演化速率差異持續(xù)擴(kuò)大。2018—2020 年，經(jīng)濟(jì)系統(tǒng)綜合評(píng)價(jià)值轉(zhuǎn)變?yōu)樨?fù)增長(zhǎng)，水資源指數(shù)下降且降幅增大，兩指數(shù)均為負(fù)增長(zhǎng)，因此兩系統(tǒng)演化速率差異減小。

4 結(jié)束語(yǔ)

在建立評(píng)價(jià)指標(biāo)體系的基礎(chǔ)上，利用耦合協(xié)調(diào)度模型和剪刀差方法分析黃河三角洲19 個(gè)縣域單元的水資源-經(jīng)濟(jì)系統(tǒng)的動(dòng)態(tài)耦合協(xié)調(diào)關(guān)系。研究發(fā)現(xiàn)，黃河三角洲水資源綜合利用水平和經(jīng)濟(jì)發(fā)展水平均呈現(xiàn)動(dòng)態(tài)提升的態(tài)勢(shì)，二者的耦合協(xié)調(diào)狀態(tài)經(jīng)歷了勉強(qiáng)協(xié)調(diào)和初級(jí)協(xié)調(diào)兩個(gè)階段，耦合類型由低水平磨合向高水平協(xié)調(diào)狀態(tài)過(guò)渡。水資源與區(qū)域經(jīng)濟(jì)系統(tǒng)剪刀差值呈現(xiàn)先上升后下降再上升的態(tài)勢(shì)。尤其是2017 年后，水資源利用和經(jīng)濟(jì)發(fā)展均經(jīng)歷了回調(diào)后的快速提升，水資源回升幅度稍大，說(shuō)明水資源支撐經(jīng)濟(jì)發(fā)展的能力有所增強(qiáng)，水資源和區(qū)域經(jīng)濟(jì)初步走向協(xié)調(diào)發(fā)展。

黃河三角洲水資源、區(qū)域經(jīng)濟(jì)系統(tǒng)綜合評(píng)價(jià)值存在一定差距，耦合協(xié)調(diào)類型差異較大，耦合協(xié)調(diào)水平空間分布不均衡。19 個(gè)縣域單元剪刀差可以分為“V 型”“W 型”“N 型”“倒U 型”4 類，以“V 型”“W 型”為主。表明各縣域水資源與區(qū)域經(jīng)濟(jì)耦合協(xié)調(diào)動(dòng)態(tài)演化的差異明顯，但大部分水資源對(duì)區(qū)域經(jīng)濟(jì)的約束減小，二者協(xié)同演化態(tài)勢(shì)增強(qiáng)。

本文在縣域尺度上分析了水資源和區(qū)域經(jīng)濟(jì)耦合協(xié)調(diào)動(dòng)態(tài)變化，未來(lái)可進(jìn)一步細(xì)化至鄉(xiāng)鎮(zhèn)街道尺度刻畫(huà)二者耦合協(xié)調(diào)更為細(xì)致的變化過(guò)程。另外，在研究變化規(guī)律的基礎(chǔ)上，需要加強(qiáng)對(duì)黃河三角洲水資源—經(jīng)濟(jì)系統(tǒng)耦合協(xié)調(diào)的動(dòng)態(tài)模擬預(yù)測(cè)，據(jù)此制定調(diào)控策略，探索協(xié)調(diào)發(fā)展路徑，為黃河三角洲的生態(tài)保護(hù)和高質(zhì)量發(fā)展提供科學(xué)依據(jù)。