王釔霏，許 朗

（南京農(nóng)業(yè)大學(xué) 經(jīng)濟(jì)管理學(xué)院，江蘇 南京 210095）

1 引言

高質(zhì)量發(fā)展是重要的時(shí)代課題，綠色發(fā)展是高質(zhì)量發(fā)展的重要標(biāo)志。 高質(zhì)量發(fā)展的驅(qū)動(dòng)因素具有復(fù)雜性特征［1］，經(jīng)濟(jì)建設(shè)、城鄉(xiāng)發(fā)展、生態(tài)保護(hù)、資源節(jié)約、產(chǎn)業(yè)布局等多方面因素都與其息息相關(guān)［2－4］。 在涉及多角度協(xié)同的黃河流域高質(zhì)量發(fā)展研究中，耦合協(xié)調(diào)度模型因能夠較好地反映多因素互動(dòng)關(guān)系和協(xié)同影響效應(yīng)，在相關(guān)研究中應(yīng)用廣泛［5－8］。 利用耦合協(xié)調(diào)度模型進(jìn)行測(cè)算前，需要對(duì)指標(biāo)賦權(quán)，常用的指標(biāo)賦權(quán)法有專家評(píng)估法、層次分析法和熵值法［9］，絕大部分文獻(xiàn)研究認(rèn)為熵值法可以利用所選指標(biāo)自身的統(tǒng)計(jì)性質(zhì)確定各自權(quán)重，指標(biāo)權(quán)重值的可信度較專家評(píng)估法和層次分析法更高，因此本文采用熵值法對(duì)指標(biāo)賦權(quán)［7，10－11］。

黃河流域?qū)儆谫Y源型缺水地區(qū)，水資源稀缺是制約黃河流域農(nóng)業(yè)綠色可持續(xù)發(fā)展和生態(tài)系統(tǒng)健康的關(guān)鍵因素之一［12－13］，同時(shí)水資源的過(guò)度開(kāi)發(fā)也對(duì)黃河流域生態(tài)環(huán)境產(chǎn)生了損害［14］，良好的農(nóng)業(yè)水資源管理可以同時(shí)提高經(jīng)濟(jì)和環(huán)境效益［15－16］。 已有研究表明，農(nóng)業(yè)節(jié)水是解決農(nóng)業(yè)生產(chǎn)用水和生態(tài)用水矛盾的關(guān)鍵手段［17］，發(fā)展節(jié)水農(nóng)業(yè)可對(duì)周邊生態(tài)環(huán)境的改善產(chǎn)生積極作用。 但也有研究發(fā)現(xiàn)，農(nóng)業(yè)節(jié)水會(huì)通過(guò)影響農(nóng)田水循環(huán)、綠洲穩(wěn)定性等路徑對(duì)農(nóng)區(qū)地表生態(tài)環(huán)境產(chǎn)生消極影響［18－19］。 此外，高效節(jié)水灌溉技術(shù)的應(yīng)用通過(guò)改變農(nóng)業(yè)生產(chǎn)方式對(duì)土壤理化性質(zhì)、地下水埋深和生物多樣性等生態(tài)環(huán)境要素產(chǎn)生影響［20－21］。 因此，進(jìn)一步探討農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)間的互動(dòng)演進(jìn)關(guān)系對(duì)推動(dòng)黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展具有重要意義。農(nóng)業(yè)用水效率是農(nóng)業(yè)節(jié)水的重要表現(xiàn)形式，本文將農(nóng)業(yè)灌溉用水技術(shù)效率作為指標(biāo)之一納入農(nóng)業(yè)節(jié)水系統(tǒng)，并采用數(shù)據(jù)包絡(luò)分析法（Data Envelopment Analysis， DEA）進(jìn)行測(cè)算。 通過(guò)梳理相關(guān)文獻(xiàn)，對(duì)于農(nóng)業(yè)水資源利用效率主要采用隨機(jī)前沿分析方法（Stochastic Frontier Analysis，SFA）和數(shù)據(jù)包絡(luò)分析法（DEA）進(jìn)行計(jì)算。 DEA 作為非參數(shù)方法， 利用線性規(guī)劃構(gòu)建邊界生產(chǎn)函數(shù)， 避免了SFA 方法對(duì)函數(shù)形式等的預(yù)先設(shè)定導(dǎo)致的參數(shù)估計(jì)有效性降低，因此本文選用DEA 方法對(duì)農(nóng)業(yè)灌溉用水技術(shù)效率進(jìn)行測(cè)算［22］。

基于以上分析，本文借鑒物理學(xué)耦合模型，分析黃河流域9 個(gè)?。ㄗ灾螀^(qū)）農(nóng)業(yè)節(jié)水與生態(tài)環(huán)境耦合協(xié)調(diào)關(guān)系的時(shí)空演進(jìn)特征，通過(guò)對(duì)比兩系統(tǒng)各指標(biāo)得分差異，探究影響黃河流域農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境協(xié)調(diào)發(fā)展的因素，以期為推進(jìn)黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展戰(zhàn)略提供針對(duì)性政策建議。

2 研究數(shù)據(jù)與研究方法

2.1 研究區(qū)域概況

黃河發(fā)源于青藏高原，流經(jīng)青海、四川、甘肅、寧夏、內(nèi)蒙古、山西、陜西、河南、山東9 個(gè)?。ㄗ灾螀^(qū)），全長(zhǎng)5 464 km，是我國(guó)僅次于長(zhǎng)江的第二大河。 根據(jù)《2019 年水資源公報(bào)》統(tǒng)計(jì)數(shù)據(jù)，黃河流域農(nóng)業(yè)用水量占比達(dá)到66.8%。 黃河流域?qū)儆谫Y源性缺水地區(qū)，因此各地區(qū)大力實(shí)施農(nóng)業(yè)節(jié)水政策以緩解農(nóng)業(yè)用水窘境。 同時(shí)，黃河流域的生態(tài)環(huán)境狀況堪憂，存在濕地蓄水量減少、生態(tài)惡化、動(dòng)物棲息地喪失、魚鳥(niǎo)種群等生物多樣性減少等多種生態(tài)環(huán)境問(wèn)題。 水資源稀缺和生態(tài)環(huán)境脆弱阻礙了黃河流域高質(zhì)量發(fā)展，積極推動(dòng)農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境健康協(xié)同發(fā)展是促進(jìn)黃河流域農(nóng)業(yè)高質(zhì)量發(fā)展的關(guān)鍵舉措。

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

本文所用數(shù)據(jù)均為公開(kāi)統(tǒng)計(jì)數(shù)據(jù)，具體來(lái)源為《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》（2010—2019 年）、《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》（2010—2019 年）、《中國(guó)統(tǒng)計(jì)年鑒》（2010—2019 年），黃河流域各?。ㄗ灾螀^(qū)）《環(huán)境狀況公報(bào)》（2010—2019 年）、《統(tǒng)計(jì)年鑒》（2010—2019 年）、水資源公報(bào)（2010—2019 年），以及各?。ㄗ灾螀^(qū)）的環(huán)保、農(nóng)業(yè)、氣象等相關(guān)部門公布的統(tǒng)計(jì)數(shù)據(jù)，主要用于黃河流域各?。ㄗ灾螀^(qū)）農(nóng)業(yè)水資源利用效率以及農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合協(xié)調(diào)度計(jì)算。

2.3 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境綜合評(píng)價(jià)指標(biāo)體系構(gòu)建

綜合考慮評(píng)價(jià)指標(biāo)對(duì)系統(tǒng)的解釋程度，同時(shí)考慮數(shù)據(jù)來(lái)源的可靠性與可操作性，構(gòu)建黃河流域農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境協(xié)同發(fā)展綜合評(píng)價(jià)指標(biāo)體系。 參考賈鳳伶等［23］、楊春紅［24］的研究構(gòu)建農(nóng)業(yè)節(jié)水指標(biāo)體系，包括3 個(gè)一級(jí)指標(biāo)，10 個(gè)二級(jí)指標(biāo)；根據(jù)陳星霖［25］、何曉瑤［19］的研究構(gòu)建生態(tài)環(huán)境指標(biāo)體系，包括3 個(gè)一級(jí)指標(biāo)，12 個(gè)二級(jí)指標(biāo)，具體指標(biāo)設(shè)置見(jiàn)表1。

表1 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)綜合評(píng)價(jià)指標(biāo)體系

2.4 研究方法

2.4.1 DEA 模型

DEA 模型是研究投入產(chǎn)出效率的重要分析工具。DEA 方法具體可以分為投入主導(dǎo)或產(chǎn)出主導(dǎo)兩種類型，投入主導(dǎo)的DEA 方法是為了測(cè)算生產(chǎn)單元相對(duì)于給定產(chǎn)出水平下最小可能投入的效率；產(chǎn)出主導(dǎo)的DEA 方法是為了度量實(shí)際產(chǎn)出與給定投入水平的最大可能產(chǎn)出差距。 本文主要考察水資源作為農(nóng)業(yè)投入要素的利用效率，因此采用投入主導(dǎo)的DEA 方法。 基于可操性、可獲得性原則，結(jié)合我國(guó)農(nóng)業(yè)水資源利用現(xiàn)狀，選取農(nóng)作物播種總面積、農(nóng)機(jī)總動(dòng)力、化肥施用量、農(nóng)業(yè)從業(yè)人數(shù)和農(nóng)業(yè)灌溉用水量作為投入指標(biāo)，選取農(nóng)業(yè)產(chǎn)值作為產(chǎn)出指標(biāo)。

此外，DEA 模型通常分為兩種模型：一種為規(guī)模報(bào)酬不變（CCR）模型，一種為規(guī)模報(bào)酬可變（BCC）模型。 因各省（自治區(qū)）在農(nóng)業(yè)投入產(chǎn)出方面并不恒定，為了更加貼近實(shí)際情況，本文采用基于規(guī)模報(bào)酬可變模型對(duì)黃河流域各省（自治區(qū)）農(nóng)業(yè)水資源利用效率進(jìn)行測(cè)算。 假設(shè)有n個(gè)決策單元，一個(gè)?。ㄗ灾螀^(qū)）為一個(gè)決策單元，每個(gè)決策單元DMUj具有m種投入（記作X1j，X2j，…，Xmj），s種產(chǎn)出（記作Y1j，Y2j，…，Ysj），得到DEA 模型表達(dá)式為

式中：S＋為產(chǎn)出松弛變量；S－為投入松弛變量；ε為非阿基米德無(wú)窮小量；E為m × j階和s × j階單位矩陣；λj為各?。ㄗ灾螀^(qū)）投入和產(chǎn)出的權(quán)向量；X為m ×j階投入矩陣；Y為s ×j階產(chǎn)出矩陣；Xj、Yj分別為第j個(gè)決策單元的投入、產(chǎn)出量；θ為決策單元的有效值。 根據(jù)DEA 法對(duì)黃河流域9 個(gè)?。ㄗ灾螀^(qū)）農(nóng)業(yè)灌溉用水技術(shù)效率的計(jì)算結(jié)果見(jiàn)表2。

表2 2010—2019 年黃河流域各?。ㄗ灾螀^(qū)）農(nóng)業(yè)灌溉用水技術(shù)效率

2.4.2 數(shù)據(jù)標(biāo)準(zhǔn)化和熵值法權(quán)重計(jì)算

統(tǒng)計(jì)口徑、指標(biāo)單位等差異使得農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境指標(biāo)存在量綱差異，同時(shí)不同指標(biāo)對(duì)系統(tǒng)存在正向和逆向兩種作用，為便于對(duì)比分析計(jì)算結(jié)果，采用極值法對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，正向指標(biāo)標(biāo)準(zhǔn)化公式為

負(fù)向指標(biāo)標(biāo)準(zhǔn)化公式為

式中：i、j分別為第i行、第j列；x′ij為標(biāo)準(zhǔn)化數(shù)據(jù)；xmax和xmin分別為某個(gè)指標(biāo)的最大值和最小值。 由于后期需要進(jìn)行熵值法測(cè)算，因此要求所有指標(biāo)值在［0，1］范圍內(nèi)。 參考相關(guān)文獻(xiàn)做法，本文將標(biāo)準(zhǔn)化處理后的數(shù)據(jù)均向右平移0. 001 個(gè)單位。

采用主觀賦值法確定指標(biāo)權(quán)重時(shí)容易產(chǎn)生較大偏差，為降低主觀因素干擾，本文采用較為常用的熵值法計(jì)算各指標(biāo)權(quán)重，具體步驟如下。

（1）數(shù)據(jù)均值化處理。 其表達(dá)式為

（2）求取指標(biāo)熵值信息。 其表達(dá)式為

（3）求取指標(biāo)權(quán)重與綜合指數(shù)。 公式為

式中：Xij為數(shù)據(jù)標(biāo)準(zhǔn)化后比值；ej為熵值；k為對(duì)象數(shù)量；zj為權(quán)重；Iij為變量得分。

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

耦合度表示系統(tǒng)相互聯(lián)系程度，協(xié)調(diào)度表示系統(tǒng)協(xié)調(diào)發(fā)展程度［26－28］。 本研究構(gòu)建的農(nóng)業(yè)節(jié)水與生態(tài)環(huán)境兩系統(tǒng)耦合協(xié)調(diào)度表達(dá)式為

式中：C為農(nóng)業(yè)節(jié)水系統(tǒng)與生態(tài)環(huán)境系統(tǒng)的耦合度；Iaws為農(nóng)業(yè)節(jié)水綜合指數(shù)；Ieco為經(jīng)濟(jì)發(fā)展綜合指數(shù)；T為農(nóng)業(yè)節(jié)水與生態(tài)環(huán)境綜合協(xié)調(diào)指數(shù)；α、β為待定系數(shù)，鑒于兩個(gè)系統(tǒng)發(fā)展同等重要，故權(quán)重相等，即α和β取值均為0.5；i和s分別為農(nóng)業(yè)節(jié)水要素和生態(tài)環(huán)境要素的編號(hào)，r、g分別為對(duì)應(yīng)要素的個(gè)數(shù)；wi和es分別為農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)各要素的標(biāo)準(zhǔn)化指標(biāo)值；ai和bs分別為農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)各要素的權(quán)重，采用熵值法計(jì)算；D為生態(tài)環(huán)境與經(jīng)濟(jì)發(fā)展的耦合協(xié)調(diào)度。

參考王少劍等［29］、田深圳等［30］、李玉紅［31］的研究成果，結(jié)合研究區(qū)實(shí)際情況，將耦合與協(xié)調(diào)程度劃分為不同階段，見(jiàn)表3。

表3 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境耦合協(xié)調(diào)度階段劃分

3 結(jié)果分析

3.1 黃河流域農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)時(shí)空演變分析

圖1 中各指標(biāo)變化趨勢(shì)表明，2010—2019 年農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)呈現(xiàn)波動(dòng)上升的趨勢(shì)，從0.293 上升到0.376。 從農(nóng)業(yè)節(jié)水系統(tǒng)一級(jí)指標(biāo)時(shí)序變化特征看，2010—2019 年節(jié)水灌溉指數(shù)和農(nóng)業(yè)用水效益指數(shù)呈現(xiàn)增長(zhǎng)趨勢(shì)，農(nóng)業(yè)用水指數(shù)呈現(xiàn)小幅度下降趨勢(shì)。其中：農(nóng)業(yè)用水指數(shù)雖有所下降，但下降幅度較小，從2010 年的0.097 降至2019 年的0.092；節(jié)水灌溉指數(shù)增長(zhǎng)較為明顯，從2010 年的0.088 上漲至2019 年的0.130，在2018 年達(dá)到研究期內(nèi)最高峰，為0.155；農(nóng)業(yè)用水效益指數(shù)增長(zhǎng)速度最快，2010—2019 年從0.088上升至0.154，且在研究期內(nèi)始終保持上升趨勢(shì)。

圖1 黃河流域農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)時(shí)序演變情況

由圖2 可知，農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)呈現(xiàn)“東高西低”省際空間格局，山東農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)最高（為0.511），最低的是青海（僅為0.158），自然資源稟賦和經(jīng)濟(jì)社會(huì)發(fā)展水平可能是形成農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)“東高西低”省際空間格局的主要驅(qū)動(dòng)因素。從研究期內(nèi)各省（自治區(qū)）農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)時(shí)序變化情況分析，山東和寧夏農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)增長(zhǎng)最快，2010—2019 年兩地農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)分別上漲0.155 和0.141。 通過(guò)對(duì)農(nóng)業(yè)節(jié)水系統(tǒng)一級(jí)指標(biāo)的省際差異分析，農(nóng)業(yè)用水效益指數(shù)也呈現(xiàn)“東高西低”的省際空間格局，東部地區(qū)經(jīng)濟(jì)發(fā)展水平較高，在農(nóng)業(yè)節(jié)水設(shè)施建設(shè)和農(nóng)業(yè)節(jié)水技術(shù)推廣上投入力度較大，單位產(chǎn)值或產(chǎn)量的需水量與西部地區(qū)相比較少［32］；節(jié)水灌溉指數(shù)在內(nèi)蒙古、山西、甘肅、寧夏等地呈現(xiàn)較高水平，這與當(dāng)?shù)厮Y源稀缺程度較高、農(nóng)民水資源稀缺意識(shí)較強(qiáng)有關(guān)［33］。

圖2 黃河流域農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)空間分布情況

3.2 黃河流域生態(tài)環(huán)境綜合發(fā)展指數(shù)時(shí)空測(cè)度分析

如圖3 所示，2010—2019 年黃河流域生態(tài)環(huán)境綜合發(fā)展指數(shù)呈現(xiàn)上升趨勢(shì)，從2010 年的0.226 上升至2019 年的0.256，其中2014 年達(dá)到研究期間峰值，為0.276。 2013 年習(xí)近平總書記發(fā)表“綠水青山就是金山銀山”的重要講話，為黃河流域生態(tài)文明建設(shè)指明了方向，推動(dòng)了黃河流域生態(tài)環(huán)境發(fā)展。 從一級(jí)指標(biāo)變化趨勢(shì)看，2010—2019 年生態(tài)彈性力、環(huán)境承載力和生態(tài)發(fā)展力均呈現(xiàn)上升趨勢(shì)，但初始水平較低，增長(zhǎng)速度較慢，導(dǎo)致生態(tài)環(huán)境綜合發(fā)展指數(shù)并未明顯增長(zhǎng)。

圖3 生態(tài)環(huán)境綜合發(fā)展指數(shù)時(shí)序演變情況

根據(jù)圖4 可知，黃河流域生態(tài)環(huán)境綜合發(fā)展指數(shù)呈現(xiàn)“東低西高”的省際空間格局和以內(nèi)蒙古、青海為核心的“雙核領(lǐng)先”的流域空間格局，生態(tài)環(huán)境綜合發(fā)展指數(shù)極差為0.207，東西部地區(qū)生態(tài)環(huán)境差異較大。黃河流域東部地區(qū)生態(tài)環(huán)境綜合發(fā)展指數(shù)普遍低于西部地區(qū)，其中：山東生態(tài)環(huán)境綜合發(fā)展指數(shù)最低，為0.147；內(nèi)蒙古生態(tài)環(huán)境綜合發(fā)展指數(shù)最高，為0.354，其次是青海，為0.344。 從生態(tài)環(huán)境系統(tǒng)一級(jí)指標(biāo)空間差異上看，內(nèi)蒙古和青海在環(huán)境承載力和生態(tài)發(fā)展力上領(lǐng)先于流域內(nèi)其他?。ㄗ灾螀^(qū)），在二者的雙重作用下內(nèi)蒙古和青海的生態(tài)環(huán)境綜合發(fā)展指數(shù)與其他?。ㄗ灾螀^(qū)）相比也具有較明顯的優(yōu)勢(shì)，因此黃河流域生態(tài)環(huán)境綜合發(fā)展指數(shù)呈現(xiàn)出內(nèi)蒙古、青海“雙核領(lǐng)先”的空間分布特征。 其主要原因是東部地區(qū)城市人口活動(dòng)和經(jīng)濟(jì)發(fā)展較為活躍，資源過(guò)度開(kāi)發(fā)和人類活動(dòng)制約了生態(tài)環(huán)境建設(shè)，生態(tài)環(huán)境綜合發(fā)展水平相對(duì)較低。

圖4 生態(tài)環(huán)境綜合發(fā)展指數(shù)空間分布情況

3.3 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合度時(shí)空演變分析

圖5 表明， 2010—2019 年黃河流域農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)長(zhǎng)期處于高水平耦合階段，耦合度基本維持在0.95 左右，說(shuō)明黃河流域的農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境兩系統(tǒng)關(guān)聯(lián)較為密切。

圖5 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合度時(shí)序演變情況

如圖6 所示，從空間維度上看各?。ㄗ灾螀^(qū)）農(nóng)業(yè)節(jié)水與生態(tài)環(huán)境耦合度無(wú)較大差異，均處于“高水平耦合”階段。 地域分布大致呈現(xiàn)“中間高，兩邊低”的特征，內(nèi)蒙古、四川和甘肅耦合度均保持在0.990 以上，與其他省份相比處于相對(duì)領(lǐng)先的水平，耦合度較低的?。ㄗ灾螀^(qū)）集中在黃河下游地區(qū)。 黃河下游泥沙堆積嚴(yán)重，導(dǎo)致灌溉用水困難，阻礙了農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境保護(hù)的相互作用。 2010—2019 年絕大部分省（自治區(qū)）兩系統(tǒng)耦合度在較高水平上呈現(xiàn)小幅度波動(dòng)上升趨勢(shì)。 山東省耦合度水平有明顯的下降趨勢(shì)，兩系統(tǒng)耦合度從2010 年的0.892 下降至2014 年的0.776，從“高水平耦合”下滑至“磨合”，2014—2019 年又逐步回升至“高水平耦合”，但2019 年山東耦合度水平較2010 年仍有所下降。

圖6 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合度時(shí)空演變情況

3.4 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境耦合協(xié)調(diào)度時(shí)空演變分析

由圖7 可知，2010—2019 年農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合協(xié)調(diào)度總體呈現(xiàn)由“基本失調(diào)”階段向“基本協(xié)調(diào)”階段過(guò)渡的增長(zhǎng)趨勢(shì)，但整體水平較低，高水平耦合協(xié)調(diào)發(fā)展趨勢(shì)不明顯；2010—2019 年各省（自治區(qū)）耦合協(xié)調(diào)度極差與方差均有所下降，說(shuō)明各?。ㄗ灾螀^(qū)）耦合協(xié)調(diào)度差距減小，整體耦合協(xié)調(diào)發(fā)展趨于穩(wěn)定；兩系統(tǒng)耦合協(xié)調(diào)度時(shí)序演進(jìn)特征以2014 年為分水嶺，2010—2013 年系統(tǒng)耦合協(xié)調(diào)度平均值在“基本失調(diào)”階段保持逐年上升趨勢(shì)，2014 年系統(tǒng)耦合協(xié)調(diào)度進(jìn)入“基本協(xié)調(diào)”階段，但2014—2019 年系統(tǒng)耦合協(xié)調(diào)度始終處于生態(tài)環(huán)境滯后的“基本協(xié)調(diào)”階段，說(shuō)明該階段農(nóng)業(yè)節(jié)水水平提升仍是促進(jìn)兩系統(tǒng)協(xié)調(diào)發(fā)展的主要驅(qū)動(dòng)力，進(jìn)一步改善黃河流域生態(tài)環(huán)境能夠有效提升兩系統(tǒng)協(xié)調(diào)發(fā)展水平。

圖7 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合協(xié)調(diào)度時(shí)序演變情況

黃河流域9 個(gè)?。ㄗ灾螀^(qū)）系統(tǒng)耦合協(xié)調(diào)度平均值介于0.436 ～0.623 之間，處于“基本失調(diào)”和“基本協(xié)調(diào)”兩個(gè)階段，整體呈現(xiàn)“東高西低”的空間分布特征（見(jiàn)圖8）。 其中：甘肅和寧夏均處于“基本失調(diào)”階段，青海處于“基本失調(diào)－農(nóng)業(yè)節(jié)水受阻”階段，農(nóng)業(yè)節(jié)水是制約青海高質(zhì)量發(fā)展的主要因素之一；其他各?。ㄗ灾螀^(qū)）均處于“基本協(xié)調(diào)”階段，其中山西、山東、河南和陜西處于生態(tài)環(huán)境滯后的“基本協(xié)調(diào)”階段，東部地區(qū)生態(tài)環(huán)境問(wèn)題是黃河流域高質(zhì)量發(fā)展進(jìn)程中需要解決的重點(diǎn)問(wèn)題。 絕大部分?。ㄗ灾螀^(qū)）耦合協(xié)調(diào)度呈現(xiàn)上升趨勢(shì)，其中：甘肅和寧夏耦合協(xié)調(diào)度增長(zhǎng)趨勢(shì)較為明顯，2010—2019 年兩地耦合協(xié)調(diào)度均呈現(xiàn)從“基本失調(diào)”階段 向“基本協(xié)調(diào)”階段邁進(jìn)的向好趨勢(shì)；寧夏耦合協(xié)調(diào)度增長(zhǎng)最快，從2010 年的0.381 增長(zhǎng)至2019 年的0.488，2018 年達(dá)到研究期內(nèi)峰值，為0.533；青海省2010—2019 年始終處于“基本失調(diào)”階段，但耦合協(xié)調(diào)度有所增長(zhǎng)；其他各?。ㄗ灾螀^(qū)）耦合協(xié)調(diào)度均在“基本協(xié)調(diào)”階段保持小幅度增長(zhǎng)趨勢(shì)。

圖8 農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合協(xié)調(diào)度時(shí)空演變情況

4 結(jié)論

本文以黃河流域9 個(gè)?。ㄗ灾螀^(qū)）的數(shù)據(jù)樣本為例，利用耦合協(xié)調(diào)度模型對(duì)農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境兩系統(tǒng)協(xié)調(diào)發(fā)展的時(shí)空演進(jìn)情況進(jìn)行測(cè)算與分析，得出以下結(jié)論。

（1）2010—2019 年，黃河流域農(nóng)業(yè)用水綜合發(fā)展指數(shù)保持逐年上升趨勢(shì)，農(nóng)業(yè)節(jié)水成效較弱是制約農(nóng)業(yè)用水綜合發(fā)展指數(shù)增長(zhǎng)的關(guān)鍵因素；在空間維度上，農(nóng)業(yè)節(jié)水綜合發(fā)展指數(shù)呈現(xiàn)“東高西低”省際空間格局，加大對(duì)西部地區(qū)水利設(shè)施建設(shè)投入、制定合理有效的農(nóng)業(yè)節(jié)水激勵(lì)機(jī)制是促進(jìn)東西部地區(qū)農(nóng)業(yè)節(jié)水平衡發(fā)展的重要措施。

（2）2010—2019 年，黃河流域生態(tài)環(huán)境綜合發(fā)展指數(shù)呈現(xiàn)上升趨勢(shì)，但其初始水平較低且增長(zhǎng)緩慢，生態(tài)環(huán)境治理成效不顯著，應(yīng)通過(guò)增加草地面積、濕地面積和提升水土流失治理率，著力提升生態(tài)發(fā)展力，促進(jìn)生態(tài)系統(tǒng)健康發(fā)展。 在空間維度上，黃河流域生態(tài)環(huán)境綜合發(fā)展指數(shù)呈現(xiàn)“東低西高”的省際空間格局，其中內(nèi)蒙古和青海生態(tài)環(huán)境綜合發(fā)展指數(shù)較其他省優(yōu)勢(shì)顯著，形成了以內(nèi)蒙古、青海為核心的“雙核領(lǐng)先”的流域空間格局。

（3）2010—2019 年，黃河流域農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)耦合度長(zhǎng)期處于高水平耦合階段，兩系統(tǒng)間相互作用的關(guān)聯(lián)度較強(qiáng)。 在空間維度上，耦合度水平呈現(xiàn)“中間高，兩邊低”的地域分布特征，黃河下游各省農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境關(guān)聯(lián)度較弱。

（4）2010—2019 年，農(nóng)業(yè)用水和生態(tài)環(huán)境系統(tǒng)耦合協(xié)調(diào)度總體呈現(xiàn)由“基本失調(diào)”階段向“基本協(xié)調(diào)”階段過(guò)渡的增長(zhǎng)趨勢(shì)，但高水平耦合協(xié)調(diào)發(fā)展趨勢(shì)不明顯，黃河流域農(nóng)業(yè)節(jié)水和生態(tài)環(huán)境系統(tǒng)整體呈現(xiàn)“高耦合－低協(xié)調(diào)”的分離特征；東部地區(qū)兩系統(tǒng)耦合協(xié)調(diào)度與西部地區(qū)相比較高，但生態(tài)環(huán)境發(fā)展滯后狀況較西部地區(qū)明顯，西部地區(qū)農(nóng)業(yè)節(jié)水對(duì)系統(tǒng)耦合協(xié)調(diào)度發(fā)展的制約作用較為顯著。