基于博弈论组合赋权-TOPSIS 模型的安康市水资源承载力综合评价

赵 浩，赵立磊，田 辉，2，丁廉超，王 凡，王郅睿，韩朝辉，龚文强，侯孝东

1.中国地质调查局 西安矿产资源调查中心，陕西 西安 710100；2.中国地质调查局 沈阳地质调查中心，辽宁 沈阳 110034

0 引言

水资源是人类赖以生存和发展的重要基础.随着经济社会的发展、人口规模的增大和城镇化进程的持续推进，人类对水资源的需求不断增加，水资源供需矛盾、水环境污染、水生态退化和地下水超采等突出问题已经成为制约水资源与经济社会可持续发展的重要因素［1-2］.水资源承载力是反映人水关系可持续发展的重要指标，科学合理地评价区域水资源承载力对于保障水资源可持续开发利用、促进经济社会发展和保护生态环境健康具有重要意义［3］.

水资源承载力的概念最早于20 世纪80 年代末由我国学者提出，此后大批学者对此展开了一系列的研究工作，从不同角度给出了水资源承载力的理解和定义，并提出了不同评价方法［1-8］.本文在施雅风、惠泱河、王建华、左其亭、Peng 等学者研究的基础上将水资源承载力定义归纳为：在某一具体的历史发展阶段，以一定的科技水平为依据，以经济社会可持续发展和水资源可持续利用为原则，以维护生态环境良性发展为条件，在水资源合理配置的前提下，一定区域上的水资源可支撑的社会、经济和生态环境协调发展的最大规模.关于水资源承载力的评价方法主要分为经验公式法、综合评价法和系统分析法3 类.经验公式法主要围绕水资源，而与经济社会和生态环境之间的系统联系不够，导致其应用较少.系统分析法虽能从水资源系统整体性、动态性和多目标性出发分析水资源-经济-社会-生态环境各个系统的内在联系，但由于计算方法复杂且不直观，难以推广应用.综合评价法由于数学理论应用比较深入，且充分考虑了水资源承载力影响因素以及水资源承载力主体与客体之间的耦合作用［1，6］，是当前应用比较多的评价方法.当前，综合评价法主要包括综合指标法［9］、主成分分析法［10］、模糊综合评判法［11］、TOPSIS 法［12-13］、VIKOR 法［14］、灰色关联法［15］和物源可拓模型法［16］等.这些方法在指标体系选择、评判标准等级构建以及评价计算上均有各自的优势和不足.其中，TOPSIS 法由于对样本数量无严格限制、信息失真小且使用便捷等特点被广泛应用于水资源承载力评价中.

基于此，本研究以水资源、社会、经济、生态环境4个子系统构建水资源承载力评价指标体系，运用博弈论组合赋权的TOPSIS 模型对安康市2011—2020 年水资源承载力进行评价，并结合障碍度模型分析影响水资源承载力的主要障碍因子，为全面掌握安康市水资源承载力状况，促进水资源可持续开发利用和管理，系统协调水资源、社会、经济和生态环境的发展关系提供科学依据.

1 研究区概况

安康市位于陕西省东南部，北依秦岭、南靠巴山，汉江横贯东西，地貌呈现南北高山夹峙、河谷盆地居中的特点.水资源以地表径流为主，区域含水岩组主要为变质岩类裂隙含水岩组和碳酸盐岩类裂隙岩溶含水岩组，同时分布少量的岩浆岩类裂隙含水岩组、碎屑岩类孔隙裂隙含水岩组和松散岩类孔隙含水岩组，各含水岩组富水性以中等富水-弱富水为主，仅在碳酸盐岩和第四系松散地层存在少量强富水地段（见图1）.全市总面积23 529 km2，年平均降水量1 050 mm，多年平均水资源总量和人均水资源量均位于陕西省前列.全市面积92.5%为山地，江河溪沟密布，地质结构复杂，气候环境多变，降水量时空分布不均，丰水期降水占全年总量70%以上，洪涝、泥石流等灾害频发，水资源开发利用率约为7.15%.作为南水北调中线工程核心水源区，安康市承担着“一泓清水永续北上”的使命和责任，由于在水资源整体性保护和整治力度上存在不足，水资源开发利用与经济社会发展格局不相匹配，在一定程度上制约了安康市的生态文明建设和可持续发展.

图1 安康市水文地质简图Fig.1 Hydrogeological sketch map of Ankang City

2 研究方法

2.1 评价指标体系及数据来源

以水资源-社会-经济-生态环境复合系统理论为基础，按照科学性、客观性和可操作性的原则系统构建水资源承载力评价指标体系.根据部分学者在陕西省水资源承载力评价［17］、汉江流域（陕西段）水生态承载力评估［18］以及汉江流域水资源承载力综合评价研究［19］的相关成果，结合安康市经济社会发展水平和水资源开发利用现状，全面分析影响水资源承载力的关键因素，选取了26 个指标构建安康市水资源承载力综合评价指标体系.针对不同指标在水资源承载力评价中的作用，将其分为正向指标和负向指标.正向指标即效益型指标，对承载力起促进作用；负向指标即成本型指标，对承载力起反作用.各指标内容、计算方法及属性类别见表1.

表1 水资源承载力综合评价指标体系Table 1 Comprehensive evaluation index system of water resources carrying capacity

研究数据来源于2011—2020 年《陕西省水资源公报》《安康市统计年鉴》《安康市国民经济和社会发展统计公报》等公开资料，各指标经整理计算获得.

2.2 组合权重计算

常用的权重计算方法分为主观赋权、客观赋权和组合赋权.主观赋权方法主要有层次分析法（AHP）［12］、古林法（KLEE）［20］和二项系数法［21］；客观赋权法主要有熵权法［11］、变异系数法［22］、CRITIC 法［23］和最大离差法［24］；组合赋权即对主、客观权重的两种或者三种进行组合计算获得综合权重，在一定程度上降低单独使用某种方法产生的误差，是目前比较常见的赋权方法，主要计算方法有乘积法［22，25］、遗传算法［23］、博弈论法［24］和代数和法［26］.本研究结合指标数据特性选择熵权法和CRITIC 法分别计算权重，运用博弈论获得各指标组合权重.熵权法能够体现指标数据之间的离散程度，而CRITIC 法能够反映指标自身的变化和指标间的相关性，两种赋权方法耦合能够更好地反映指标权重［27］.

2.2.1 熵权法

熵权法是根据各评价指标信息效用值的差异程度来确定权重系数的客观赋权法.指标信息熵越小，则变异程度越大，在评价中重要程度也就越高，因此权重也就越大.反之，信息熵越大，权重越小［25］.主要计算步骤如下：

1）构建原始评价指标矩阵X

设有m 个评价指标，n 个评价对象，xij为第i 个评价指标在第j 个评价对象上的对应数值，则原始评价矩阵为：

2）归一化、正向化处理

由于各指标类型不同、量纲不一，为了排除这些差异的影响，需要对各指标进行无量纲化处理.对于熵权法，宜采用极值法［28］进行处理，得到标准化矩阵Y.

式中：yij为第j 个评价对象在第i 个评价指标上的标准值.

3）计算第i 个评价指标的信息熵

4）计算权重w1i

2.2.2 改进的CRITIC 法

CRITIC 法（criteria importance through intercriteria correlation）是由Diakoulaki 提出的通过计算指标的变异性和冲突性来综合衡量指标权重的一种客观赋权法.变异性指同一指标针对不同对象取值的差异性，通常用标准差来度量.冲突性指各指标之间的相关性大小，通常用相关系数来表示.通用计算方法见文献［23］，该计算公式存在两个问题［29-30］：1）各指标量纲、数量级不同，用标准差反映差异性存在准确性低、误差大的缺陷；2）指标间的相关系数有可能存在负值，而绝对值相同的正相关与负相关其相关性相同.针对此问题，本研究采用文献［23］的方法进行改进：1）运用变异系数度量指标差异性；2）通过相关系数绝对值来计算指标之间的冲突性.设Ci代表第i 个评价指标的信息量，表达公式为：

式中：σi为第i 个评价指标的标准差，Xi为第i 个评价指标的均值，σi/Xi即为第i 个评价指标的变异系数，rij为第i 个评价指标与第j 个评价指标间的相关系数.则第i 个评价指标的权重w2i计算公式为：

2.2.3 博弈论组合赋权法

基于博弈论将熵权法和改进的CRITIC 法相结合求得组合权重w，计算方法如下［31］：

1）通过熵权法和CRITIC 法计算各指标权重，分别得到权重集w1和w2，设a1、a2为线性组合系数，将权重集w1和w2进行线性组合：

2）根据博弈论思想，对（10）式中的线性组合系数a1、a2进行优化，以组合权重w 和w1、w2离差最小化为目标确定目标函数：

3）根据矩阵微分性质，确定（11）式最优化一阶导数条件的线性微分方程组为：

4）由（12）式计算得到优化组合系数a1、a2，归一化处理得到a1*=a1/（a1＋a2），a2*=a2/（a1＋a2），则基于博弈论组合赋权的综合权重w 为：

2.3 TOPSIS 模型

TOPSIS 法（technique for order preference by similarity to an ideal solution）即“逼近理想解排序法”，是由C.L.Hwang 和K.Yoon 于1981 年首次提出的一种多目标决策分析方法.该方法能够充分利用原始数据信息，精确地反映各评价方案之间的差距，对数据分布及样本含量没有严格限制，数据计算简单易行，是一种常见的组内综合评价方法.通过构造评价对象与正、负理想解之间的距离来判断评价对象的优劣程度.具体计算步骤如下［27］：

1）原始数据同趋势化、归一化.将负向指标进行倒数处理，然后采用公式（14）计算归一化值，即由原始矩阵X 得到无量纲决策矩阵Z＝（zij）m×n

2）构造加权决策矩阵V

3）计算正负理想解Ai＋、Ai-

4）计算各评价对象与正负理想解的距离Adj＋、Adj-

5）计算贴近度Gj

式中：Gj越接近于1，证明该对象与正理想解越靠近，评价结果越优，反之越差.

2.4 障碍度模型

通过障碍度模型进行障碍因子诊断，是在水资源承载力评价的基础上进一步识别影响水资源承载力的主要障碍因素，有助于认识区域水资源开发利用现状，提高水资源承载力水平.某项指标的障碍度越大，说明该指标对水资源承载力的阻碍作用越强，反之越弱.计算时引入因子贡献度（Fi）、指标偏离度（Iij）和障碍度（Pij），公式如下［32］：

式中：wi为各指标所占权重；Ti为各子系统所占的权重；zij为标准化处理后的指标值.

3 计算结果

3.1 权重计算结果

根据2.2 节熵权法和改进的CRITIC 法分别计算权重集w1和w2，通过博弈论法计算出优化组合系数a1*=0.3688，a2*=0.6312，进而计算各指标的组合权重w.分配结果列于表2.

表2 博弈论组合赋权的指标权重分配结果Table 2 Index weight distribution results of game theory combination weighting

3.2 贴近度Gj 计算结果

根据各指标组合权重和TOPSIS 模型评价方法及步骤，对2011—2020 年安康市水资源承载力及各子系统进行综合评价，并计算水资源承载力复合系统与各子系统评价值的相关性，结果见表3.

表3 安康市2011—2020 年水资源承载力评价结果Table 3 Evaluation results of water resources carrying capacity in Ankang City during 2011-2020

4 评价结果分析

4.1 复合系统水资源承载力分析

由图2 可知，安康市水资源承载力变化趋势呈“W”型，贴近度2019 年最高，2013 年最低.总体上可分为3 个阶段：第一阶段为2011—2013 年，贴近度值从0.4793 降低至0.3292，下降31.32%，下降幅度较大，主要表现为水资源、社会和生态环境3 个子系统承载力评价值均呈下降趋势，尤其是水资源子系统，降幅高达96.55%；第二阶段为2013—2017 年，贴近度值从0.3292 上升到0.6133，增长86.30%，增长幅度较大，主要表现为4 个子系统评价值均呈上升趋势，水资源、社会、经济和生态环境子系统分别增长2460.06%、6.67%、77.93%、16.05%；第三阶段为2017—2020 年，贴近度呈波动态势，虽然2019 年贴近度值最高，但仍处于下降趋势，主要表现为水资源、经济和生态环境3个子系统评价值均呈波动下降趋势，尤其是水资源和生态环境子系统，分别下降46.90%和27.08%.从表3中的均值和相关系数分析可知，经济子系统对水资源承载力的贡献最大，其次为水资源子系统、社会子系统和生态环境子系统.

图2 安康市2011—2020 年水资源承载力评价值Fig.2 Evaluation values of water resources carrying capacity in Ankang City during 2011-2020

综合来看，在2011—2020 年整个评价期内，安康市水资源承载力总体变化不大，呈波动状态，承载力水平处于中等，表明水资源仍有较大开发利用空间.

4.2 各子系统水资源承载力分析

根据表3 和图2 评价结果，对各子系统承载力进行分析.

4.2.1 水资源子系统

水资源子系统承载力波动幅度很大，总体呈下降趋势，与复合系统水资源承载力变化趋势基本一致.2011—2013 年贴近度由0.9934 骤降至0.0343，降幅高达96.55%，2013—2020 年贴近度呈波动上升趋势，但在2016、2018、2020 年均较前一年有明显的下降.从统计数据分析可知，水资源子系统贴近度值与产水模数、人均水资源量变化趋势保持一致.2011—2013 年间，产水模数从66.17×104m3/km2减少到28.64×104m3/km2，人均水资源量从5 907.63 m3/人下降到2 576.15 m3/人，分别下降56.72%和56.39%.该时段内人口密度从112.01 人/ km2降低到111.15 人/ km2，说明在常住总人口变化不大的情况下，水资源总量的减少，使得水资源供需矛盾加大，从而导致水资源承载力的骤降.

4.2.2 社会子系统

社会子系统承载力呈波动上升趋势，贴近度从2011 年的0.3118 增长到2020 年的0.6952，增幅达122.96%.综合来看存在2 个上升阶段和3 个下降阶段：上升阶段为2013—2015 年、2016—2019 年，前者贴近度值从0.2818 增长至0.5525，后者从0.3005 增长至0.7139，分别增长96.06%和137.57%.下降阶段为2011—2013 年、2015—2016 年、2019—2020 年，贴近度值分别从0.3118 下降至0.2818，0.5525 下降至0.3006，0.7139 下降至0.6952.从原始数据来看，上升段主要是由于人口密度、人口自然增长率和人均日生活用水量的降低以及城镇化率的升高综合作用的结果.下降段2011—2013 年人口自然增长率从2.57‰增长到2.85‰，相比人口密度、城镇化率和城镇人均日生活用水量变化幅度大，导致承载力水平略微降低.而2015—2016 年和2019—2020 年两个阶段主要是由于人口自然增长和城镇化进程的加快，人均日生活用水量快速升高，尤其是2020 年从115.91 L 增长到125.56 L，增幅达8.33%，用水量的增加导致承载力水平下降.表明在城镇化进程加快的过程中大量人口进入城市，随着用水量的增大，政府部门对水资源的节约管控和群众对水资源的节约意识需要加强和提升.

4.2.3 经济子系统

经济子系统承载力呈上升趋势，贴近度值从2011年的0.2635 增长到2020 年的0.6051，增长129.64%.从原始数据分析可知，主要原因是2011—2020 年人均国内生产总值（GDP）和第三产业比重的增长以及万元GDP 用水量、万元工业增加值用水量、万元农业增加值用水量大幅降低综合作用的结果.其中人均GDP 从15 107 元增长到43 666 元，第三产业比重从40.00%增至45.28%，万元GDP 用水量从184.4 m3/万元下降到71.4 m3/万元，万元工业增加值用水量从69.76 m3/万元下降到24.55 m3/万元，万元农业增加值用水量从935.66 m3/万元下降到405.35 m3/万元，变化幅度分别为＋189.04%、＋5.28%、-61.28%、-64.81%、-56.68%.

此外，2014、2020 年贴近度值较上一年有所降低，2014 年降低的主要原因为农田灌溉亩均用水量较上一年增长达59.34%，对承载力水平有较大影响，而2020 年主要是因为人均GDP 降低以及万元GDP 用水量和万元工业增加值用水量增大所致.

综上可知，随着经济规模的快速增长和产业结构的调整，以及水资源利用率的快速提升，安康市经济子系统承载力能力已经达到一定水平，但仍需进一步优化产业结构，降低水资源量的消耗，使经济发展保持良好的增长态势.

4.2.4 生态环境子系统

生态环境子系统承载力变化趋势总体呈“V”型，可分为2011—2015 年的缓慢下降阶段和2015—2020年的波动上升阶段.2011—2015 年贴近度值从0.4956下降至0.3322，降幅32.97%.从原始数据来看，生态环境子系统承载力下降的主要原因是生态环境用水率、污水排放量、单位耕地施肥量的增加和水土流失治理面积的减少.计算可知，2015 年生态环境用水率、污水排放量和单位耕地施肥量较2011 年分别增长130.77%、35.93%、16.75%，可见生态环境用水量的增加加大了对水资源的负载，虽然此阶段城市污水日处理能力和城市污水处理率有了较大提升，但仍不足以抵消污水排放量对水资源承载力的影响程度.2015—2020 年贴近度值从0.3322 上升至0.3760，属于波动上升阶段.分析数据可知，森林覆盖率和城市绿化覆盖率的增加、城市污水日处理能力和城市污水处理率的提升以及化学需氧量、氨氮排放量的大幅降低是承载力上升的主要原因.其中，森林覆盖率和城市绿化覆盖率分别达到67.95%和40.15%，污水日处理能力达到21.47×104t/d，城市污水处理率为96.69%.而化学需氧量排放量和氨氮排放量分别较2015 年下降16.38%和46.32%.

总体上，从生态环境用水率、污水日处理能力和城市污水处理率的快速增长以及化学需氧量和氨氮排放量的大幅度下降，可见安康市在经济发展过程中对环境保护的重视程度越来越高，尤其是“十三五”期间对水污染方面的整治力度有了显著增强.然而，从图2可以看出，2020 年贴近度值较2019 年有一定的下降.从原始数据分析可知，2020 年人均GDP 较2019 年降低7.00%，而化学需氧量排放量和氨氮排放量分别上升55.35%和13.94%.虽然经济处于负增长状态，但部分污染物排放量反而有所增长，说明在污染防治方面有所松懈，应该引起重视.

4.3 水资源承载力障碍度因素诊断

4.3.1 各子系统主要障碍因素分析

利用公式（19）—（21）计算2011—2020 年各指标障碍度，可得各子系统承载力主要障碍因子及变化趋势，如图3 所示.

图3 2011—2020 年安康市水资源承载力各指标障碍度Fig.3 Obstacle degrees of indexes for water resources carrying capacity in Ankang City during 2011-2020

从图3a 可知，水资源子系统主要障碍因子为产水模数、人均水资源量和水资源开发利用率，且三者障碍度总体呈上升趋势.从原始数据分析可知，2011年以来，产水模数从66.17×104m3/km2降低至46.47×104m3/km2，人均水资源量从5 907.63 m3/人降低至4 359.51 m3/人，二者呈波动下降趋势，而水资源开发利用率从4.72%上升至7.15%.水资源子系统中承载力受气候变化和人类活动双重影响，气候变化作为不可抗因素时，人类活动对水资源的需求在逐步加大，因此，应结合水资源量进行合理配置和开发利用.

从图3b 可知，社会子系统各指标障碍度总体小于1.0%，普遍较小，承载力主要障碍因子为人口自然增长率和城镇人均日生活用水量，而人口密度和城镇化率影响相对较小.

从图3c 可知，经济子系统各指标障碍度总体呈下降趋势，主要障碍因子为万元工业增加值用水量、人均GDP、万元GDP 用水量和农田灌溉亩均用水量.从原始数据分析可知，万元工业增加值用水量从69.76 m3/万元下降至24.55 m3/万元，人均GDP 从15 107 元增长至43 666 元，万元GDP 用水量从184.4 m3/万元下降至71.4 m3/万元，农田灌溉亩均用水量从449.0 m3/亩增长至701.1 m3/亩.可以看出，随着经济快速发展和人民生活水平的提高，安康市逐步将高能耗的工矿企业转向绿色低碳的农牧渔产业，单位GDP 用水量大幅下降，而农田灌溉亩均用水量却在大幅上涨.因此在农业灌溉上要采用节约集约用水技术，提升节水意识，降低农用灌溉用水量.

从图3d 可知，生态环境子系统各指标障碍度普遍较高，其主要障碍因子为氨氮排放量、化学需氧量排放量、城市污水日处理能力、单位耕地施肥量、生态环境用水率和污水排放量.除生态环境用水率和污水排放量的障碍度逐步增高外，其余指标障碍度均呈下降趋势.从原始数据分析可知，“十二五”到“十三五”期间，随着生态环境保护意识的增强和环保力度的加大，化学需氧量排放量和氨氮排放量分别降低22.66%和50.47%，而生态环境用水率却从0.39%上升至1.69%，说明人们对生态环境的重视越来越强，生态环境对水资源的需求在逐步加大.虽然污水排放总量在快速增长，但污水处理率已经达到96.69%，排入环境中的实际排污量在不断降低.整体来看，生态环境子系统各指标障碍度对水资源承载力作用较强，应继续加大环境保护力度，不断降低污染物排放量，继续提升污水日处理能力和处理率，减少生态环境污染，提升水资源承载能力.

4.3.2 复合系统主要障碍因素分析

根据各指标障碍度计算可得各子系统相对复合系统障碍度（图4），并选取障碍度前10 位的指标作为主要障碍因子列于表4.

表4 2011—2020 年安康市水资源承载力主要障碍因子及障碍度Table 4 Main obstacle factors and obstacle degrees of water resources carrying capacity in Ankang City during 2011-2020

图4 2011—2020 年安康市水资源承载力各子系统障碍度Fig.4 Obstacle degrees of subsystems for water resources carrying capacity in Ankang City during 2011-2020

由图4 可知，2011—2020 年评价期内，生态环境子系统障碍度最高，其次为经济子系统、水资源子系统和社会子系统.其中，生态环境和水资源子系统障碍度分别升高3.28%和20.16%，经济和社会子系统障碍度分别下降18.27%和4.71%.

从表4 分析可知，安康市水资源承载力复合系统主要障碍因子为氨氮排放量、化学需氧量排放量、生态环境用水率、城市污水日处理能力、污水排放量和万元工业增加值用水量，其次为单位耕地施肥量、产水模数和人均水资源量，其他指标障碍度相对较小.

总体变化可分为两个阶段：第一个阶段是2011—2015 年，障碍度以氨氮排放量、化学需氧量排放量、城市污水日处理能力和万元工业增加值用水量为主，其中氨氮排放量和化学需氧量排放量障碍度基本不变，万元工业增加值用水量和城市污水日处理能力障碍度在逐渐降低，而生态环境用水率和污水排放量的障碍度在快速升高.从原始数据可知，该时段万元工业增加值用水量从69.76 m3/万元大幅降低至20.48 m3/万元，城市污水日处理能力从6×104t/d 提升至16×104t/d，污水排放量从1 078.90×104t/d 增长到1 466.56×104t/d，生态环境用水率也从0.39%提升至0.90%.第二阶段是2016—2020 年，生态环境用水率的障碍度提升至首位，其次为氨氮排放量、化学需氧量排放量、污水排放量和城市污水日处理能力.从数据可知，到2020 年，环境用水率已增长至1.69%，而氨氮排放量已经从2015 年的4 538 t 降低至2 436 t，化学需氧量排放量也从34 529 t 降低至28 872 t，但污水排放量已经增长至3292×104t.

以上数据表明，“十二五”期间，随着经济的快速发展，安康市开始重视工业用水的控制，在逐步优化工业产业结构，向节水型社会发展，同时也加大了生态环境的管理与建设，但在氨氮、化学需氧量排放量上还未进行有效的控制.而在“十三五”期间，随着生态文明建设的加快推进，建设资源节约型和环境友好型社会以及山水林田湖生命共同体等理念的提升，安康市作为南水北调中线工程水源地核心区对生态环境保护越来越重视.从各子系统障碍度和主要指标障碍度综合来看，应进一步统筹协调水资源、经济、社会和生态环境复合系统，加大环境保护力度，不断优化产业结构，通过技术革新进一步降低水资源消耗量和污染物排放量，使水资源承载力水平在经济社会发展的过程中逐步提升.

5 结论

本研究以水资源-社会-经济-生态环境复合系统理论为基础构建安康市水资源承载力评价指标体系，采用博弈论组合赋权法，运用TOPSIS 模型和障碍度模型对安康市2011—2020 年水资源承载力进行综合评价，主要结论如下：

1）从承载力综合水平来看，2011—2020 年安康市水资源承载力变化趋势呈“W”型，综合评价值从0.4793 下降至0.4740，总体变化不大，水资源仍有较大开发利用空间.其中，经济子系统对水资源承载力的贡献最大，其次为水资源子系统、社会子系统和生态环境子系统.

2）从子系统承载力来看，2011—2020 年水资源和生态环境子系统承载力总体呈下降趋势，社会和经济子系统承载力呈上升趋势.在2020 年，各子系统承载力均有所下降.

3）障碍度分析表明，生态环境子系统障碍度最高，其次为经济子系统和水资源子系统，社会子系统障碍度最低.其中，氨氮排放量、化学需氧量排放量、生态环境用水率、城市污水日处理能力和万元工业增加值用水量是主要的障碍因子.