南水北調中線水源涵養(yǎng)區(qū)（河南段）生態(tài)功能與高質量發(fā)展耦合協(xié)調研究

魏杰 程翔 劉贊 邢昱 申進朝 朱紅雷

摘 要：為厘清近年來南水北調中線水源涵養(yǎng)區(qū)（河南段）生態(tài)功能與經濟社會高質量發(fā)展的關系，依據(jù)遙感資料與水資源、環(huán)境、社會經濟統(tǒng)計數(shù)據(jù)，參照《生態(tài)環(huán)境狀況評價技術規(guī)范》（ＨＪ １９２—２０１５），構建生態(tài)功能和高質量發(fā)展評價指標體系，采用耦合協(xié)調度模型分析２０１４—２０１８ 年水源涵養(yǎng)區(qū)（河南段）生態(tài)功能和高質量發(fā)展趨勢及耦合協(xié)調特征。結果表明：２０１４—２０１８ 年水源涵養(yǎng)區(qū)（河南段）６ 縣（市）原始生態(tài)功能指數(shù)均呈增大趨勢，增幅范圍為１．０４～８．０７，對于標準化后的生態(tài)功能指數(shù)，除淅川縣和鄧州市整體表現(xiàn)為增大外，其他區(qū)域均表現(xiàn)為略微減?。唬?縣（市）高質量發(fā)展極不均衡，鄧州市經濟基礎好、發(fā)展快，歷年高質量發(fā)展指數(shù)均高于０．６５０，盧氏縣發(fā)展平穩(wěn)但水平較低，２０１８ 年高質量發(fā)展指數(shù)最小，為０．０９４，與鄧州市差距明顯；６ 縣（市）生態(tài)功能與高質量發(fā)展耦合協(xié)調水平差別較大，西峽縣為優(yōu)質協(xié)調，欒川縣、內鄉(xiāng)縣和淅川縣為良好協(xié)調，鄧州市和盧氏縣為中度失調。

關鍵詞：水源涵養(yǎng)；生態(tài)功能；高質量發(fā)展；耦合協(xié)調度；南水北調；河南段

中圖分類號： Ｘ８２６ 文獻標志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２３．０７．０１９

引用格式：魏杰，程翔，劉贊，等．南水北調中線水源涵養(yǎng)區(qū)（ 河南段） 生態(tài)功能與高質量發(fā)展耦合協(xié)調研究［ Ｊ］．人民黃河，２０２３，４５（７）：１０３－１０８．

南水北調中線工程總干渠從河南省淅川縣陶岔渠首引水，跨越長江流域、淮河流域、黃河流域、海河流域，經過１２ 個地級市到達北京、天津，為北京、天津等地飲用水安全提供有力保障，對我國水資源優(yōu)化配置起到非常大的作用［１－２］ 。２０１５ 年河南省人民政府辦公廳出臺關于丹江口水庫（河南轄區(qū)）飲用水水源地保護區(qū)劃的規(guī)定，制定了保護措施［３－４］ 。２０２１ 年習近平總書記主持召開推進南水北調后續(xù)工程高質量發(fā)展座談會并發(fā)表重要講話，強調要加強南水北調工程沿線水資源保護，持續(xù)抓好輸水沿線區(qū)和受水區(qū)的污染防治和生態(tài)環(huán)境保護工作。這一重要講話為推進南水北調后續(xù)工程高質量發(fā)展指明了方向。經濟社會高質量發(fā)展與生態(tài)環(huán)境保護存在相互耦合的關系，經濟社會發(fā)展對生態(tài)環(huán)境保護具有脅迫作用，環(huán)境惡化對經濟社會發(fā)展有約束作用。目前，國內外學者從不同視角對經濟社會發(fā)展與生態(tài)環(huán)境保護之間的協(xié)調關系進行了研究，研究方法主要有系統(tǒng)動力學法［５－６］ 、環(huán)境庫茲涅茨曲線（ＥＫＣ）法［７－８］ 、ＰＶＡＲ 模型［９－１０］ 、多元線性回歸分析法［１１］ 、生態(tài)足跡法［１２－１３］ 和耦合協(xié)調度模型［１４－１５］ 等。其中耦合協(xié)調度模型可以定量表征２ 個及以上系統(tǒng)間配合得當、和諧一致的程度，既可反映各系統(tǒng)間相互作用的強弱程度，又可反映各系統(tǒng)是在高水平上相互促進還是低水平上相互制約［１６］ ，目前該模型在城市化與自然災害評價［１７］ 、人口空間分布與城市發(fā)展評價［１８］ 、生態(tài)系統(tǒng)服務與水資源評價等方面得到了廣泛應用。

本文依據(jù)遙感影像資料與水資源、環(huán)境、社會經濟統(tǒng)計數(shù)據(jù)，參照《生態(tài)環(huán)境狀況評價技術規(guī)范》（ＨＪ１９２—２０１５），構建生態(tài)功能和高質量發(fā)展評價指標體系，分析南水北調中線丹江口水源區(qū)（河南段）水源涵養(yǎng)生態(tài)功能和高質量發(fā)展趨勢及特征，并對兩者的耦合協(xié)調度進行研究，以期為丹江口水源區(qū)經濟社會高質量發(fā)展和生態(tài)保護提供參考。

１ 研究區(qū)概況與研究方法

１．１ 研究區(qū)概況

南水北調中線水源涵養(yǎng)區(qū)（河南段）位于河南省西南部，西北高東南低，區(qū)域地貌類型復雜，以山地、丘陵為主。丹江口水庫匯水區(qū)域（河南段）涉及南陽市淅川縣、西峽縣、內鄉(xiāng)縣和鄧州市，三門峽市盧氏縣，洛陽市欒川縣（分布區(qū)域示意見圖１），總面積７ ８１５ ｋｍ２，其中南陽市匯水區(qū)面積６ ３６１ ｋｍ２，占河南段匯水區(qū)面積的８１％。丹江口水庫水源涵養(yǎng)區(qū)（河南段）氣候屬于亞熱帶季風氣候，降水充沛，年平均氣溫１５．８ °Ｃ，多年平均降水量８０４．３ ｍｍ，汛期６—９ 月降水量占全年降水量的５８％。本文以鄧州市、淅川縣、西峽縣、內鄉(xiāng)縣、欒川縣和盧氏縣及南水北調中線水源涵養(yǎng)區(qū)（河南段）全段為評價單元，構建評價指標體系，對６ 縣（市）水源涵養(yǎng)生態(tài)功能與高質量發(fā)展耦合協(xié)調關系進行研究。

１．２ 數(shù)據(jù)來源

遙感數(shù)據(jù)源自２０１４—２０１５ 年Ｌａｎｄｓａｔ８ ＯＬＩ 影像（空間分辨率為３０ ｍ）和２０１６—２０１８ 年高分一號、二號影像（空間分辨率統(tǒng)一重采樣至２ ｍ），通過人機交互解譯獲得土地利用類型數(shù)據(jù)；數(shù)字高程模型（ＤＥＭ）數(shù)據(jù)源自分辨率為３０ ｍ 的ＡＳＴＥＲ ＧＤＥＭ 模型；環(huán)境數(shù)據(jù)源自２０１５—２０１９ 年《河南省環(huán)境統(tǒng)計公報》；高質量發(fā)展數(shù)據(jù)源自２０１５—２０１９ 年《河南省統(tǒng)計年鑒》。

１．３ 生態(tài)功能與高質量發(fā)展耦合協(xié)調度模型

１．３．１ 評價指標體系構建

１）生態(tài)功能評價指標體系。本文參照《生態(tài)環(huán)境狀況評價技術規(guī)范》（ＨＪ １９２—２０１５），選擇評價指標并確定權重（見表１）。水源涵養(yǎng)區(qū)生態(tài)功能指數(shù)（ＦＥＩ）計算模型如下：ＦＥＩ ＝０．１５×水源涵養(yǎng)指數(shù)＋０．１２×受保護區(qū)域面積比×１００＋０．０９×林地覆蓋率＋０．０６×草地覆蓋率＋０．０９×水域濕地面積比＋０．０９×（１００－耕地與建設用地面積比） ＋０．１８×（１００－主要污染物排放強度）＋０．０４×污染源排放達標率×１００＋０．０４×城鎮(zhèn)污水集中處理率×１００＋０．０８×水質達標率×１００＋０．０４×空氣質量達標率×１００＋０．０２×集中式飲用水水源地水質達標率×１００。

２ 結果分析

２．１ 生態(tài)功能評價對生態(tài)功能指數(shù)１２ 個分指標進行標準化處理，加權計算得到２０１４—２０１８ 年６ 縣（市）以及南水北調中線水源涵養(yǎng)區(qū)（河南段）全段生態(tài)功能指數(shù)標準化值（見表４）。生態(tài)功能指數(shù)標準化值范圍為０．０８０ ～０．８４２，除淅川縣和鄧州市的生態(tài)功能指數(shù)標準化值整體表現(xiàn)為增大外，其他區(qū)域均表現(xiàn)為略微減小，鄧州市較其他５ 縣歷年指數(shù)均明顯偏低，欒川縣、盧氏縣、西峽縣和淅川縣指數(shù)較高，空間分異明顯。欒川縣、盧氏縣和西峽縣林草地覆蓋率均在８０．０％以上，淅川縣水域濕地面積比高達４４．０％，遠高于其他縣域５．０％的水平；內鄉(xiāng)縣處于中等水平，林草地覆蓋率達５７．０％；鄧州市為平原地區(qū)，林草地覆蓋率低于２．０％，耕地面積比高達７８．６％，建設用地面積比為１６．０％，遠高于整個研究區(qū)耕地面積比２７．０％、建設用地面積比５．８％的水平。

對未進行標準化處理的分指標值進行加權計算，得到原始生態(tài)功能指數(shù)（ＦＥＩ），見表５。按照ＦＥＩ 值和原始生態(tài)功能指數(shù)年際變化值（ΔＦＥＩ）劃分生態(tài)功能等級與生態(tài)功能變化等級：ＦＥＩ＜４０ 時生態(tài)功能等級為差，４０≤ＦＥＩ＜５０ 時生態(tài)功能等級為較差，５０≤ＦＥＩ＜６０時生態(tài)功能等級為一般，６０≤ＦＥＩ＜７０ 時生態(tài)功能等級為良，ＦＥＩ≥７０ 時生態(tài)功能等級為優(yōu)；ΔＦＥＩ≤－４ 時生態(tài)功能變化等級為顯著變差，－４＜ΔＦＥＩ≤－２ 時生態(tài)功能變化等級為明顯變差，－２＜ΔＦＥＩ≤－１ 時生態(tài)功能變化等級為略微變差，－１＜ΔＦＥＩ＜１ 時生態(tài)功能變化等級為無明顯變化，１≤ΔＦＥＩ＜２ 時生態(tài)功能變化等級為略微變好，２≤ΔＦＥＩ ＜４ 時生態(tài)功能變化等級為明顯變好，ΔＦＥＩ≥４ 時生態(tài)功能變化等級為顯著變好。得到生態(tài)功能變化分級空間分布圖，見圖２。由于２０１５—２０１６ 年南水北調中線水源涵養(yǎng)區(qū)６ 縣（市）生態(tài)功能變化等級均為無明顯變化，因此２０１５—２０１６ 年生態(tài)功能變化分級空間分布圖省略。２０１４ 年６ 縣（市）ＦＥＩ值范圍為３６．４９～７０．９３，２０１８ 年ＦＥＩ 值范圍為４０．９６～７１．９７，２０１４—２０１８ 年ＦＥＩ 的最小值和最大值均有所增大。２０１４ 年生態(tài)功能等級為優(yōu)、良、一般和差的縣（市）個數(shù)分別為１、３、１ 和１，這４ 種級別的縣（市）面積分別占南水北調中線水源涵養(yǎng)區(qū)（河南段）總面積的２０．２％、５２．５％、１３．５％和１３．８％。至２０１８ 年，優(yōu)、良和較差的縣（市）個數(shù)分別為２、３ 和１，差類型未出現(xiàn)，這３ 種級別的縣（市）面積分別占南水北調中線水源涵養(yǎng)區(qū)（河南段） 總面積的４１．７％、４４．５％和１３．８％。２０１８ 年等級為優(yōu)的縣（市）面積明顯增加，處于較差類型的縣（市）面積占比較小。

對原始生態(tài)功能指數(shù)２０１４—２０１８ 年變化情況進行分析，所有縣（市）的ΔＦＥＩ 值均為正值，即ＦＥＩ 值均呈增大趨勢，說明水源涵養(yǎng)功能明顯提高。尤其２０１５年河南省人民政府辦公廳出臺關于丹江口水庫（河南轄區(qū)）飲用水水源地保護區(qū)劃的規(guī)定，６ 縣（市）切實執(zhí)行該規(guī)定，加大污染減排和生態(tài)環(huán)境治理力度，城鎮(zhèn)污水集中處理率明顯提高，使得２０１６—２０１７ 年生態(tài)功能變化尤為明顯，盧氏縣、西峽縣、內鄉(xiāng)縣、鄧州市以及南水北調中線水源涵養(yǎng)區(qū)（河南段） 的ＦＥＩ 值升高?？傮w來看，２０１４—２０１８ 年６ 縣（市）原始生態(tài)功能指數(shù)均呈上升趨勢，盧氏縣、內鄉(xiāng)縣、淅川縣和鄧州市生態(tài)功能顯著變好，淅川縣原始生態(tài)功能指數(shù)增幅最大（為８．０７），欒川縣和西峽縣的生態(tài)功能略微變好，原始生態(tài)功能指數(shù)增幅較小，都為１．０４。

２．２ 高質量發(fā)展評價

對高質量發(fā)展評價指標體系的１２ 個分指標進行標準化處理，加權計算得到２０１４—２０１８ 年南水北調中線水源涵養(yǎng)區(qū)（河南段）６ 縣（市）高質量發(fā)展指數(shù)（見表６）?？梢钥闯?，鄧州市作為國務院確定的丹江口庫區(qū)區(qū)域中心城市、南陽市區(qū)域副中心城市，發(fā)展較快，生產總值高，歷年高質量發(fā)展指數(shù)均高于０．６５０；２０２０年之前盧氏縣屬于國家級扶貧開發(fā)重點縣，基本被崤山、熊耳山、伏牛山覆蓋，是河南省面積最大的縣，卻也是河南最貧窮的縣，經濟實力排名全省最后，發(fā)展平穩(wěn)但水平較低，２０１４—２０１８ 年高質量發(fā)展指數(shù)為０．０６６～０．０９８，與鄧州市的差距明顯；欒川縣和西峽縣發(fā)展水平較高，內鄉(xiāng)縣和淅川縣發(fā)展較不顯著。

２．３ 耦合協(xié)調綜合評價

通過式（３） ～ 式（５） 計算得到２０１４—２０１８ 年南水北調水源涵養(yǎng)區(qū)（河南段） 各縣（市） 生態(tài)功能與高質量發(fā)展耦合協(xié)調度指數(shù)，見表７。可以看出，２０１４—２０１８ 年南水北調中線水源涵養(yǎng)區(qū)（河南段）６縣（市）Ｄ 值為０．１９１ ～０．８５８，內鄉(xiāng)縣、盧氏縣和鄧州市Ｄ 值呈波動增長趨勢，其中鄧州市Ｄ 值增幅最大，為０．２３４。

各縣（市）生態(tài)功能與高質量發(fā)展耦合協(xié)調存在差異，西峽縣耦合協(xié)調類型為優(yōu)質協(xié)調，說明西峽縣生態(tài)功能與高質量發(fā)展相互作用處于良性階段，發(fā)展均衡。欒川縣、內鄉(xiāng)縣和淅川縣為良好協(xié)調，說明３ 縣生態(tài)功能與高質量發(fā)展耦合協(xié)調度較高，發(fā)展比較均衡，但生態(tài)功能與高質量發(fā)展均有提升空間。鄧州市整體為中度失調，從表５～表７ 可以看出，鄧州市經濟社會發(fā)展迅速，發(fā)展水平較其他５ 縣高，但生態(tài)功能相對較弱，兩者發(fā)展明顯失衡。盧氏縣整體為中度失調，生態(tài)功能作用明顯，但高質量發(fā)展指數(shù)較小，兩者耦合協(xié)調度低。

從耦合協(xié)調發(fā)展類型來看，西峽縣生態(tài)功能與高質量發(fā)展耦合協(xié)調發(fā)展類型為同步型，生態(tài)功能與高質量發(fā)展相互促進作用明顯。其他縣（市）均為不均衡發(fā)展，內鄉(xiāng)縣、欒川縣和淅川縣為生態(tài)功能主導型，盧氏縣為高質量發(fā)展滯后型，鄧州市２０１４—２０１７ 年為生態(tài)功能損益型、２０１８ 年為生態(tài)功能滯后型。這充分證明近年來國家生態(tài)補償相關政策的實施，極大促進了水源涵養(yǎng)區(qū)部分縣（市）高質量發(fā)展，提升了生態(tài)環(huán)境質量。但生態(tài)功能與高質量發(fā)展耦合協(xié)調發(fā)展類型分級明顯，針對所處階段，建議欒川縣、淅川縣、盧氏縣在提升水源涵養(yǎng)功能的同時，嚴格產業(yè)準入機制，積極調整產業(yè)結構，從傳統(tǒng)依靠資源的制作業(yè)向創(chuàng)新驅動的新產業(yè)轉型，發(fā)展綠色經濟。鄧州市應創(chuàng)新環(huán)境保護和治理模式，發(fā)展經濟的同時不破壞生態(tài)的穩(wěn)定性和完整性，控制工業(yè)和建設用地規(guī)模，大力實施生態(tài)工程，以生態(tài)功能高水平促進經濟社會高質量發(fā)展。

３ 結論與建議

３．１ 結論

本文以南水北調中線水源涵養(yǎng)區(qū)（河南段）為研究區(qū)，構建生態(tài)功能和高質量發(fā)展評價指標體系，建立耦合協(xié)調度模型，分析生態(tài)功能與高質量發(fā)展之間的耦合協(xié)調情況，得出結論如下。

１）通過分析２０１４—２０１８ 年南水北調中線水源涵養(yǎng)區(qū)（河南段）６ 縣（市）原始生態(tài)功能指數(shù)變化情況，發(fā)現(xiàn)６ 縣（市）指數(shù)均呈上升趨勢，增幅范圍為１．０４～８．０７。其中：盧氏縣、內鄉(xiāng)縣、淅川縣和鄧州市生態(tài)功能顯著變好，淅川縣原始生態(tài)功能指數(shù)增幅最大，為８．０７，欒川縣和西峽縣生態(tài)功能略微變好。然而，對于標準化后的生態(tài)功能指數(shù)，２０１４—２０１８ 年除淅川縣和鄧州市表現(xiàn)為增大外，其他縣域均表現(xiàn)為略微減小。

２）分析２０１４—２０１８ 年南水北調中線水源涵養(yǎng)區(qū)（河南段）６ 縣（市）高質量發(fā)展指數(shù)變化情況，發(fā)現(xiàn)各縣（市）發(fā)展極不均衡，鄧州市經濟基礎好、發(fā)展快，歷年高質量發(fā)展指數(shù)均高于０．６５０；而盧氏縣高質量發(fā)展指數(shù)最小，２０１４—２０１８ 年高質量發(fā)展指數(shù)為０．０６６ ～０．０９８。

３）２０１４—２０１８ 年南水北調中線水源涵養(yǎng)區(qū)（河南段）各縣（市）生態(tài)功能與高質量發(fā)展耦合協(xié)調類型存在差異，西峽縣耦合協(xié)調類型為優(yōu)質協(xié)調，欒川縣、內鄉(xiāng)縣和淅川縣為良好協(xié)調，鄧州市和盧氏縣整體為中度失調。南水北調中線水源涵養(yǎng)區(qū)（河南段）整體耦合協(xié)調發(fā)展水平較高，但不同縣（市）耦合協(xié)調發(fā)展不均衡，高質量發(fā)展有較大提升空間，鄧州市水源涵養(yǎng)生態(tài)功能質量需穩(wěn)步提高。

３．２ 建議

對于南水北調中線水源涵養(yǎng)區(qū)，國家有明確的產業(yè)發(fā)展要求以及嚴格的產業(yè)準入機制。為彌補國家政策帶來的經濟損失，提出建議：逐年增加橫向生態(tài)補償資金，大力支持發(fā)展第三產業(yè)；經濟社會高質量發(fā)展水平較高的受水區(qū)積極補償水源地在生態(tài)環(huán)境保護中支出的資金成本，加大援助力度；積極開展生態(tài)功能提升、生態(tài)保護和補償資金使用的科學評估，為水源涵養(yǎng)區(qū)的生態(tài)補償提供科學依據(jù)，切實促進水源地生態(tài)環(huán)境質量提高、鄉(xiāng)村振興和經濟社會高質量發(fā)展。

參考文獻：

［１］ 劉占坤，路祥翼，陳軍，等．南水北調中線工程沿線地區(qū)用水空間競合關系研究［Ｊ］．中國農村水利水電，２０２１（６）：１９１－１９６．

［２］ 仲志余，劉國強，吳澤宇．南水北調中線工程水量調度實踐及分析［Ｊ］．南水北調與水利科技，２０１８，１６（１）：９５－９９，１４３．

［３］ 邢領航，范治暉，陳進，等．南水北調中線水源地飲用水源保護區(qū)劃分研究［Ｊ］．人民長江，２０１５，４６（１９）：８３－８９．

［４］ 黃文清，谷樹忠．南水北調中線水源地保護法制保障能力分析［Ｊ］．人民長江，２０１７，４８（２３）：１８－２１，３９．

［５］ 陳華君，褚鈺，付景保．南水北調中線水源區(qū)生態(tài)產業(yè)與環(huán)境耦合發(fā)展情景分析［Ｊ］．水資源保護，２０２２，３８（６）：１９４－２０１．

［６］ ＷＵ Ｄ Ｓ，ＮＩＮＧ Ｓ．Ｄｙｎａｍｉｃ Ａｓｓｅｓｓｍｅｎｔ ｏｆ Ｕｒｂａｎ Ｅｃｏｎｏｍｙ?Ｅｎｖｉｒｏｎｍｅｎｔ?Ｅｎｅｒｇｙ Ｓｙｓｔｅｍ Ｕｓｉｎｇ Ｓｙｓｔｅｍ Ｄｙｎａｍｉｃｓ Ｍｏｄｅｌ：Ａ Ｃａｓｅ Ｓｔｕｄｙ ｉｎ Ｂｅｉｊｉｎｇ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｒｅｓｅａｒｃｈ，２０１８，１６４：７０－８４．

［７］ 張陳俊，章恒全．新環(huán)境庫茲涅次曲線：工業(yè)用水與經濟增長的關系［Ｊ］．中國人口·資源與環(huán)境，２０１４，２４（５）：１１６－１２３．

［８］ ＫＯＯＮＤＨＡＲ Ｍ Ａ，ＳＨＡＨＢＡＺ Ｍ，ＭＥＭＯＮ Ｋ Ａ，ｅｔ ａｌ．Ａ Ｖｉ?ｓｕａｌｉｚａｔｉｏｎ Ｒｅｖｉｅｗ Ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ Ｌａｓｔ Ｔｗｏ Ｄｅｃａｄｅｓ ｆｏｒ Ｅｎ?ｖｉｒｏｎｍｅｎｔａｌ Ｋｕｚｎｅｔｓ Ｃｕｒｖｅ “ ＥＫＣ” Ｂａｓｅｄ ｏｎ Ｃｏ?ＣｉｔａｔｉｏｎＡｎａｌｙｓｉｓ Ｔｈｅｏｒｙ ａｎｄ Ｐａｔｈｆｉｎｄｅｒ Ｎｅｔｗｏｒｋ Ｓｃａｌｉｎｇ Ａｌｇｏｒｉｔｈｍｓ［Ｊ］．Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｃｉｅｎｃｅ ａｎｄ Ｐｏｌｌｕｔｉｏｎ Ｒｅｓｅａｒｃｈ，２０２１，２８（１３）：１６６９０－１６７０６．

［９］ 李茜，胡昊，羅海江，等．我國經濟增長與環(huán)境污染雙向作用關系研究：基于ＰＶＡＲ 模型的區(qū)域差異分析［Ｊ］．環(huán)境科學學報，２０１５，３５（６）：１８７５－１８８６．

［１０］ ＣＨＡＲＦＥＤＤＩＮＥ Ｌ，ＫＡＨＩＡ Ｍ．Ｉｍｐａｃｔ ｏｆ Ｒｅｎｅｗａｂｌｅ ＥｎｅｒｇｙＣｏｎｓｕｍｐｔｉｏｎ ａｎｄ Ｆｉｎａｎｃｉａｌ Ｄｅｖｅｌｏｐｍｅｎｔ ｏｎ ＣＯ２ Ｅｍｉｓｓｉｏｎｓａｎｄ Ｅｃｏｎｏｍｉｃ Ｇｒｏｗｔｈ ｉｎ ｔｈｅ ＭＥＮＡ Ｒｅｇｉｏｎ：Ａ Ｐａｎｅｌ ＶｅｃｔｏｒＡｕｔｏｒｅｇｒｅｓｓｉｖｅ （ＰＶＡＲ） Ａｎａｌｙｓｉｓ［Ｊ］．Ｒｅｎｅｗａｂｌｅ Ｅｎｅｒｇｙ，２０１９，１３９：１９８－２１３．

［１１］ 杜鑫，許東，付曉，等．遼河流域遼寧段水環(huán)境演變與流域經濟發(fā)展的關系［Ｊ］．生態(tài)學報，２０１５，３５（６）：１９５５－１９６０．

［１２］ ＤＥＳＴＥＫ Ｍ Ａ，ＳＩＮＨＡ Ａ．Ｒｅｎｅｗａｂｌｅ，Ｎｏｎ?Ｒｅｎｅｗａｂｌｅ ＥｎｅｒｇｙＣｏｎｓｕｍｐｔｉｏｎ，Ｅｃｏｎｏｍｉｃ Ｇｒｏｗｔｈ，Ｔｒａｄｅ Ｏｐｅｎｎｅｓｓ ａｎｄ Ｅｃｏｌｏｇ?ｉｃａｌ Ｆｏｏｔｐｒｉｎｔ：Ｅｖｉｄｅｎｃｅ ｆｒｏｍ Ｏｒｇａｎｉｓａｔｉｏｎ ｆｏｒ Ｅｃｏｎｏｍｉｃ Ｃｏ?Ｏｐｅｒａｔｉｏｎ ａｎｄ Ｄｅｖｅｌｏｐｍｅｎｔ Ｃｏｕｎｔｒｉｅｓ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ ＣｌｅａｎｅｒＰｒｏｄｕｃｔｉｏｎ，２０２０，２４２：１１８５３７．

［１３］ ＤＯＧＡＮ Ｅ，ＵＬＵＣＡＫ Ｒ，ＫＯＣＡＫ Ｅ，ｅｔ ａｌ．Ｔｈｅ Ｕｓｅ ｏｆ Ｅｃｏ?ｌｏｇｉｃａｌ Ｆｏｏｔｐｒｉｎｔ ｉｎ Ｅｓｔｉｍａｔｉｎｇ ｔｈｅ Ｅｎｖｉｒｏｎｍｅｎｔａｌ ＫｕｚｎｅｔｓＣｕｒｖｅ Ｈｙｐｏｔｈｅｓｉｓ ｆｏｒ ＢＲＩＣＳＴ ｂｙ Ｃｏｎｓｉｄｅｒｉｎｇ Ｃｒｏｓｓ?ＳｅｃｔｉｏｎＤｅｐｅｎｄｅｎｃｅ ａｎｄ Ｈｅｔｅｒｏｇｅｎｅｉｔｙ ［Ｊ］． Ｓｃｉｅｎｃｅ ｏｆ Ｔｈｅ ＴｏｔａｌＥｎｖｉｒｏｎｍｅｎｔ，２０２０，７２３：１３８０６３．

［１４］ 孫斌，徐渭，薛建春，等．黃河流域城市群城鎮(zhèn)化與生態(tài)保護耦合協(xié)調研究［Ｊ］．人民黃河，２０２２，４４（６）：１６－２１．

［１５］ ＳＨＩ Ｔ，ＹＡＮＧ Ｓ Ｙ，ＺＨＡＮＧ Ｗ，ｅｔ ａｌ．Ｃｏｕｐｌｉｎｇ ＣｏｏｒｄｉｎａｔｉｏｎＤｅｇｒｅｅ Ｍｅａｓｕｒｅｍｅｎｔ ａｎｄ Ｓｐａｔｉｏｔｅｍｐｏｒａｌ Ｈｅｔｅｒｏｇｅｎｅｉｔｙ Ｂｅ?ｔｗｅｅｎ Ｅｃｏｎｏｍｉｃ Ｄｅｖｅｌｏｐｍｅｎｔ ａｎｄ Ｅｃｏｌｏｇｉｃａｌ Ｅｎｖｉｒｏｎｍｅｎｔ：Ｅｍｐｉｒｉｃａｌ Ｅｖｉｄｅｎｃｅ ｆｒｏｍ Ｔｒｏｐｉｃａｌ ａｎｄ Ｓｕｂｔｒｏｐｉｃａｌ Ｒｅｇｉｏｎｓｏｆ Ｃｈｉｎａ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｃｌｅａｎｅｒ Ｐｒｏｄｕｃｔｉｏｎ，２０２０，２４４：１１８７３９．

［１６］ 王成，唐寧．重慶市鄉(xiāng)村三生空間功能耦合協(xié)調的時空特征與格局演化［Ｊ］．地理研究，２０１８，３７（６）：１１００－１１１４．

［１７］ ＺＨＡＮＧ Ｚ Ｘ，ＬＩ Ｙ．Ｃｏｕｐｌｉｎｇ Ｃｏｏｒｄｉｎａｔｉｏｎ ａｎｄ ＳｐａｔｉｏｔｅｍｐｏｒａｌＤｙｎａｍｉｃ Ｅｖｏｌｕｔｉｏｎ Ｂｅｔｗｅｅｎ Ｕｒｂａｎｉｚａｔｉｏｎ ａｎｄ ＧｅｏｌｏｇｉｃａｌＨａｚａｒｄｓ：Ａ Ｃａｓｅ Ｓｔｕｄｙ ｆｒｏｍ Ｃｈｉｎａ［Ｊ］．Ｓｃｉｅｎｃｅ ｏｆ Ｔｈｅ ＴｏｔａｌＥｎｖｉｒｏｎｍｅｎｔ，２０２０，７２８：１３８８２５．

［１８］ ＨＡＮ Ｈ Ｌ，ＬＩ Ｈ．Ｃｏｕｐｌｉｎｇ Ｃｏｏｒｄｉｎａｔｉｏｎ Ｅｖａｌｕａｔｉｏｎ ＢｅｔｗｅｅｎＰｏｐｕｌａｔｉｏｎ ａｎｄ Ｌａｎｄ Ｕｒｂａｎｉｚａｔｉｏｎ ｉｎ Ｈａ?Ｃｈａｎｇ Ｕｒｂａｎ Ａｇ?ｇｌｏｍｅｒａｔｉｏｎ［Ｊ］．Ｓｕｓｔａｉｎａｂｉｌｉｔｙ，２０２０，１２（１）：３５７．

［１９］ 唐曉靈，杜莉．基于引力模型的區(qū)域經濟發(fā)展與生態(tài)環(huán)境耦合協(xié)調研究：以陜西省為例［Ｊ］．生態(tài)經濟，２０２０，３６（７）：１６４－１６９．

［２０］ 逯進，常虹，汪運波．中國區(qū)域能源、經濟與環(huán)境耦合的動態(tài)演化［Ｊ］．中國人口·資源與環(huán)境，２０１７，２７（２）：６０－６８．

【責任編輯 栗 銘】