南水北调中线水源涵养区（河南段）生态功能与高质量发展耦合协调研究

魏杰 程翔 刘赞 邢昱 申进朝 朱红雷

摘 要：为厘清近年来南水北调中线水源涵养区（河南段）生态功能与经济社会高质量发展的关系，依据遥感资料与水资源、环境、社会经济统计数据，参照《生态环境状况评价技术规范》（ＨＪ １９２—２０１５），构建生态功能和高质量发展评价指标体系，采用耦合协调度模型分析２０１４—２０１８ 年水源涵养区（河南段）生态功能和高质量发展趋势及耦合协调特征。结果表明：２０１４—２０１８ 年水源涵养区（河南段）６ 县（市）原始生态功能指数均呈增大趋势，增幅范围为１．０４～８．０７，对于标准化后的生态功能指数，除淅川县和邓州市整体表现为增大外，其他区域均表现为略微减小；６ 县（市）高质量发展极不均衡，邓州市经济基础好、发展快，历年高质量发展指数均高于０．６５０，卢氏县发展平稳但水平较低，２０１８ 年高质量发展指数最小，为０．０９４，与邓州市差距明显；６ 县（市）生态功能与高质量发展耦合协调水平差别较大，西峡县为优质协调，栾川县、内乡县和淅川县为良好协调，邓州市和卢氏县为中度失调。

关键词：水源涵养；生态功能；高质量发展；耦合协调度；南水北调；河南段

中图分类号： Ｘ８２６ 文献标志码：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２３．０７．０１９

引用格式：魏杰，程翔，刘赞，等．南水北调中线水源涵养区（ 河南段） 生态功能与高质量发展耦合协调研究［ Ｊ］．人民黄河，２０２３，４５（７）：１０３－１０８．

南水北调中线工程总干渠从河南省淅川县陶岔渠首引水，跨越长江流域、淮河流域、黄河流域、海河流域，经过１２ 个地级市到达北京、天津，为北京、天津等地饮用水安全提供有力保障，对我国水资源优化配置起到非常大的作用［１－２］ 。２０１５ 年河南省人民政府办公厅出台关于丹江口水库（河南辖区）饮用水水源地保护区划的规定，制定了保护措施［３－４］ 。２０２１ 年习近平总书记主持召开推进南水北调后续工程高质量发展座谈会并发表重要讲话，强调要加强南水北调工程沿线水资源保护，持续抓好输水沿线区和受水区的污染防治和生态环境保护工作。这一重要讲话为推进南水北调后续工程高质量发展指明了方向。经济社会高质量发展与生态环境保护存在相互耦合的关系，经济社会发展对生态环境保护具有胁迫作用，环境恶化对经济社会发展有约束作用。目前，国内外学者从不同视角对经济社会发展与生态环境保护之间的协调关系进行了研究，研究方法主要有系统动力学法［５－６］ 、环境库兹涅茨曲线（ＥＫＣ）法［７－８］ 、ＰＶＡＲ 模型［９－１０］ 、多元线性回归分析法［１１］ 、生态足迹法［１２－１３］ 和耦合协调度模型［１４－１５］ 等。其中耦合协调度模型可以定量表征２ 个及以上系统间配合得当、和谐一致的程度，既可反映各系统间相互作用的强弱程度，又可反映各系统是在高水平上相互促进还是低水平上相互制约［１６］ ，目前该模型在城市化与自然灾害评价［１７］ 、人口空间分布与城市发展评价［１８］ 、生态系统服务与水资源评价等方面得到了广泛应用。

本文依据遥感影像资料与水资源、环境、社会经济统计数据，参照《生态环境状况评价技术规范》（ＨＪ１９２—２０１５），构建生态功能和高质量发展评价指标体系，分析南水北调中线丹江口水源区（河南段）水源涵养生态功能和高质量发展趋势及特征，并对两者的耦合协调度进行研究，以期为丹江口水源区经济社会高质量发展和生态保护提供参考。

１ 研究区概况与研究方法

１．１ 研究区概况

南水北调中线水源涵养区（河南段）位于河南省西南部，西北高东南低，区域地貌类型复杂，以山地、丘陵为主。丹江口水库汇水区域（河南段）涉及南阳市淅川县、西峡县、内乡县和邓州市，三门峡市卢氏县，洛阳市栾川县（分布区域示意见图１），总面积７ ８１５ ｋｍ２，其中南阳市汇水区面积６ ３６１ ｋｍ２，占河南段汇水区面积的８１％。丹江口水库水源涵养区（河南段）气候属于亚热带季风气候，降水充沛，年平均气温１５．８ °Ｃ，多年平均降水量８０４．３ ｍｍ，汛期６—９ 月降水量占全年降水量的５８％。本文以邓州市、淅川县、西峡县、内乡县、栾川县和卢氏县及南水北调中线水源涵养区（河南段）全段为评价单元，构建评价指标体系，对６ 县（市）水源涵养生态功能与高质量发展耦合协调关系进行研究。

１．２ 数据来源

遥感数据源自２０１４—２０１５ 年Ｌａｎｄｓａｔ８ ＯＬＩ 影像（空间分辨率为３０ ｍ）和２０１６—２０１８ 年高分一号、二号影像（空间分辨率统一重采样至２ ｍ），通过人机交互解译获得土地利用类型数据；数字高程模型（ＤＥＭ）数据源自分辨率为３０ ｍ 的ＡＳＴＥＲ ＧＤＥＭ 模型；环境数据源自２０１５—２０１９ 年《河南省环境统计公报》；高质量发展数据源自２０１５—２０１９ 年《河南省统计年鉴》。

１．３ 生态功能与高质量发展耦合协调度模型

１．３．１ 评价指标体系构建

１）生态功能评价指标体系。本文参照《生态环境状况评价技术规范》（ＨＪ １９２—２０１５），选择评价指标并确定权重（见表１）。水源涵养区生态功能指数（ＦＥＩ）计算模型如下：ＦＥＩ ＝０．１５×水源涵养指数＋０．１２×受保护区域面积比×１００＋０．０９×林地覆盖率＋０．０６×草地覆盖率＋０．０９×水域湿地面积比＋０．０９×（１００－耕地与建设用地面积比） ＋０．１８×（１００－主要污染物排放强度）＋０．０４×污染源排放达标率×１００＋０．０４×城镇污水集中处理率×１００＋０．０８×水質达标率×１００＋０．０４×空气质量达标率×１００＋０．０２×集中式饮用水水源地水质达标率×１００。

２ 结果分析

２．１ 生态功能评价对生态功能指数１２ 个分指标进行标准化处理，加权计算得到２０１４—２０１８ 年６ 县（市）以及南水北调中线水源涵养区（河南段）全段生态功能指数标准化值（见表４）。生态功能指数标准化值范围为０．０８０ ～０．８４２，除淅川县和邓州市的生态功能指数标准化值整体表现为增大外，其他区域均表现为略微减小，邓州市较其他５ 县历年指数均明显偏低，栾川县、卢氏县、西峡县和淅川县指数较高，空间分异明显。栾川县、卢氏县和西峡县林草地覆盖率均在８０．０％以上，淅川县水域湿地面积比高达４４．０％，远高于其他县域５．０％的水平；内乡县处于中等水平，林草地覆盖率达５７．０％；邓州市为平原地区，林草地覆盖率低于２．０％，耕地面积比高达７８．６％，建设用地面积比为１６．０％，远高于整个研究区耕地面积比２７．０％、建设用地面积比５．８％的水平。

对未进行标准化处理的分指标值进行加权计算，得到原始生态功能指数（ＦＥＩ），见表５。按照ＦＥＩ 值和原始生态功能指数年际变化值（ΔＦＥＩ）划分生态功能等级与生态功能变化等级：ＦＥＩ＜４０ 时生态功能等级为差，４０≤ＦＥＩ＜５０ 时生态功能等级为较差，５０≤ＦＥＩ＜６０时生态功能等级为一般，６０≤ＦＥＩ＜７０ 时生态功能等级为良，ＦＥＩ≥７０ 时生态功能等级为优；ΔＦＥＩ≤－４ 时生态功能变化等级为显著变差，－４＜ΔＦＥＩ≤－２ 时生态功能变化等级为明显变差，－２＜ΔＦＥＩ≤－１ 时生态功能变化等级为略微变差，－１＜ΔＦＥＩ＜１ 时生态功能变化等级为无明显变化，１≤ΔＦＥＩ＜２ 时生态功能变化等级为略微变好，２≤ΔＦＥＩ ＜４ 时生态功能变化等级为明显变好，ΔＦＥＩ≥４ 时生态功能变化等级为显著变好。得到生态功能变化分级空间分布图，见图２。由于２０１５—２０１６ 年南水北调中线水源涵养区６ 县（市）生态功能变化等级均为无明显变化，因此２０１５—２０１６ 年生态功能变化分级空间分布图省略。２０１４ 年６ 县（市）ＦＥＩ值范围为３６．４９～７０．９３，２０１８ 年ＦＥＩ 值范围为４０．９６～７１．９７，２０１４—２０１８ 年ＦＥＩ 的最小值和最大值均有所增大。２０１４ 年生态功能等级为优、良、一般和差的县（市）个数分别为１、３、１ 和１，这４ 种级别的县（市）面积分别占南水北调中线水源涵养区（河南段）总面积的２０．２％、５２．５％、１３．５％和１３．８％。至２０１８ 年，优、良和较差的县（市）个数分别为２、３ 和１，差类型未出现，这３ 种级别的县（市）面积分别占南水北调中线水源涵养区（河南段） 总面积的４１．７％、４４．５％和１３．８％。２０１８ 年等级为优的县（市）面积明显增加，处于较差类型的县（市）面积占比较小。

对原始生态功能指数２０１４—２０１８ 年变化情况进行分析，所有县（市）的ΔＦＥＩ 值均为正值，即ＦＥＩ 值均呈增大趋势，说明水源涵养功能明显提高。尤其２０１５年河南省人民政府办公厅出台关于丹江口水库（河南辖区）饮用水水源地保护区划的规定，６ 县（市）切实执行该规定，加大污染减排和生态环境治理力度，城镇污水集中处理率明显提高，使得２０１６—２０１７ 年生态功能变化尤为明显，卢氏县、西峡县、内乡县、邓州市以及南水北调中线水源涵养区（河南段） 的ＦＥＩ 值升高。总体来看，２０１４—２０１８ 年６ 县（市）原始生态功能指数均呈上升趋势，卢氏县、内乡县、淅川县和邓州市生态功能显著变好，淅川县原始生态功能指数增幅最大（为８．０７），栾川县和西峡县的生态功能略微变好，原始生态功能指数增幅较小，都为１．０４。

２．２ 高质量发展评价

对高质量发展评价指标体系的１２ 个分指标进行标准化处理，加权计算得到２０１４—２０１８ 年南水北调中线水源涵养区（河南段）６ 县（市）高质量发展指数（见表６）。可以看出，邓州市作为国务院确定的丹江口库区区域中心城市、南阳市區域副中心城市，发展较快，生产总值高，历年高质量发展指数均高于０．６５０；２０２０年之前卢氏县属于国家级扶贫开发重点县，基本被崤山、熊耳山、伏牛山覆盖，是河南省面积最大的县，却也是河南最贫穷的县，经济实力排名全省最后，发展平稳但水平较低，２０１４—２０１８ 年高质量发展指数为０．０６６～０．０９８，与邓州市的差距明显；栾川县和西峡县发展水平较高，内乡县和淅川县发展较不显著。

２．３ 耦合协调综合评价

通过式（３） ～ 式（５） 计算得到２０１４—２０１８ 年南水北调水源涵养区（河南段） 各县（市） 生态功能与高质量发展耦合协调度指数，见表７。可以看出，２０１４—２０１８ 年南水北调中线水源涵养区（河南段）６县（市）Ｄ 值为０．１９１ ～０．８５８，内乡县、卢氏县和邓州市Ｄ 值呈波动增长趋势，其中邓州市Ｄ 值增幅最大，为０．２３４。

各县（市）生态功能与高质量发展耦合协调存在差异，西峡县耦合协调类型为优质协调，说明西峡县生态功能与高质量发展相互作用处于良性阶段，发展均衡。栾川县、内乡县和淅川县为良好协调，说明３ 县生态功能与高质量发展耦合协调度较高，发展比较均衡，但生态功能与高质量发展均有提升空间。邓州市整体为中度失调，从表５～表７ 可以看出，邓州市经济社会发展迅速，发展水平较其他５ 县高，但生态功能相对较弱，两者发展明显失衡。卢氏县整体为中度失调，生态功能作用明显，但高质量发展指数较小，两者耦合协调度低。

从耦合协调发展类型来看，西峡县生态功能与高质量发展耦合协调发展类型为同步型，生态功能与高质量发展相互促进作用明显。其他县（市）均为不均衡发展，内乡县、栾川县和淅川县为生态功能主导型，卢氏县为高质量发展滞后型，邓州市２０１４—２０１７ 年为生态功能损益型、２０１８ 年为生态功能滞后型。这充分证明近年来国家生态补偿相关政策的实施，极大促进了水源涵养区部分县（市）高质量发展，提升了生态环境质量。但生态功能与高质量发展耦合协调发展类型分级明显，针对所处阶段，建议栾川县、淅川县、卢氏县在提升水源涵养功能的同时，严格产业准入机制，积极调整产业结构，从传统依靠资源的制作业向创新驱动的新产业转型，发展绿色经济。邓州市应创新环境保护和治理模式，发展经济的同时不破坏生态的稳定性和完整性，控制工业和建设用地规模，大力实施生态工程，以生态功能高水平促进经济社会高质量发展。

３ 结论与建议

３．１ 结论

本文以南水北调中线水源涵养区（河南段）为研究区，构建生态功能和高质量发展评价指标体系，建立耦合协调度模型，分析生态功能与高质量发展之间的耦合协调情况，得出结论如下。

１）通过分析２０１４—２０１８ 年南水北调中线水源涵养区（河南段）６ 县（市）原始生态功能指数变化情况，发现６ 县（市）指数均呈上升趋势，增幅范围为１．０４～８．０７。其中：卢氏县、内乡县、淅川县和邓州市生态功能显著变好，淅川县原始生态功能指数增幅最大，为８．０７，栾川县和西峡县生态功能略微变好。然而，对于标准化后的生态功能指数，２０１４—２０１８ 年除淅川县和邓州市表现为增大外，其他县域均表现为略微减小。

２）分析２０１４—２０１８ 年南水北调中线水源涵养区（河南段）６ 县（市）高质量发展指数变化情况，发现各县（市）发展极不均衡，邓州市经济基础好、发展快，历年高质量发展指数均高于０．６５０；而卢氏县高质量发展指数最小，２０１４—２０１８ 年高质量发展指数为０．０６６ ～０．０９８。

３）２０１４—２０１８ 年南水北调中线水源涵养区（河南段）各县（市）生态功能与高质量发展耦合协调类型存在差异，西峡县耦合协调类型为优质协调，栾川县、内乡县和淅川县为良好协调，邓州市和卢氏县整体为中度失调。南水北调中线水源涵养区（河南段）整体耦合协调发展水平较高，但不同县（市）耦合协调发展不均衡，高质量发展有较大提升空间，邓州市水源涵养生态功能质量需稳步提高。

３．２ 建议

对于南水北调中线水源涵养区，国家有明确的产业发展要求以及严格的产业准入机制。为弥补国家政策带来的经济损失，提出建议：逐年增加横向生态补偿资金，大力支持发展第三产业；经济社会高质量发展水平较高的受水区积极补偿水源地在生态环境保护中支出的资金成本，加大援助力度；积极开展生态功能提升、生态保护和补偿资金使用的科学评估，为水源涵养区的生态补偿提供科学依据，切实促进水源地生态环境质量提高、乡村振兴和经济社会高质量发展。

参考文献：

［１］ 刘占坤，路祥翼，陈军，等．南水北调中线工程沿线地区用水空间竞合关系研究［Ｊ］．中国农村水利水电，２０２１（６）：１９１－１９６．

［２］ 仲志余，刘国强，吴泽宇．南水北调中线工程水量调度实践及分析［Ｊ］．南水北调与水利科技，２０１８，１６（１）：９５－９９，１４３．

［３］ 邢领航，范治晖，陈进，等．南水北调中线水源地饮用水源保护区划分研究［Ｊ］．人民长江，２０１５，４６（１９）：８３－８９．

［４］ 黄文清，谷树忠．南水北调中线水源地保护法制保障能力分析［Ｊ］．人民长江，２０１７，４８（２３）：１８－２１，３９．

［５］ 陈华君，褚钰，付景保．南水北调中线水源区生态产业与环境耦合发展情景分析［Ｊ］．水资源保护，２０２２，３８（６）：１９４－２０１．

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

［７］ 张陈俊，章恒全．新环境库兹涅次曲线：工业用水与经济增长的关系［Ｊ］．中国人口·资源与环境，２０１４，２４（５）：１１６－１２３．

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

［９］ 李茜，胡昊，罗海江，等．我国经济增长与环境污染双向作用关系研究：基于ＰＶＡＲ 模型的区域差异分析［Ｊ］．环境科学学报，２０１５，３５（６）：１８７５－１８８６．

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

［１１］ 杜鑫，許东，付晓，等．辽河流域辽宁段水环境演变与流域经济发展的关系［Ｊ］．生态学报，２０１５，３５（６）：１９５５－１９６０．

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

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

［１４］ 孙斌，徐渭，薛建春，等．黄河流域城市群城镇化与生态保护耦合协调研究［Ｊ］．人民黄河，２０２２，４４（６）：１６－２１．

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

［１６］ 王成，唐宁．重庆市乡村三生空间功能耦合协调的时空特征与格局演化［Ｊ］．地理研究，２０１８，３７（６）：１１００－１１１４．

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

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

［１９］ 唐晓灵，杜莉．基于引力模型的区域经济发展与生态环境耦合协调研究：以陕西省为例［Ｊ］．生态经济，２０２０，３６（７）：１６４－１６９．

［２０］ 逯进，常虹，汪运波．中国区域能源、经济与环境耦合的动态演化［Ｊ］．中国人口·资源与环境，２０１７，２７（２）：６０－６８．

【责任编辑 栗 铭】