城市的扩张和工业的发展使得人类活动对城市河流的干扰和破坏不断加剧, 进而导致了水系衰退、水体黑臭等水生态环境问题^[1]。随着近些年来城市河流生态健康受到社会的广泛关注与重视, 一系列生态修复措施得以实施以期能够改善城市河流的健康状态^[2]。因此, 进行城市河流的健康评价, 一方面可以检验生态修复措施的效果, 同时也是进行城市河流管理与保护的重要前提, 对城市的可持续发展具有重要意义。

目前, 河流健康评价方法主要包括指示物种法和指标体系法。指示物种法采用鱼类、浮游动植物、底栖动物以及微生物等生物类群的物种完整性来监测河流生态系统健康, 适用于受人类活动影响较小的自然河流的健康评价^[3-5]。指标体系法则是根据河流系统特征构建评价指标体系, 采用数学模型确定河流健康等级, 在河流健康评价中得到了广泛的应用^[6-7]。城市河流作为人类活动和自然过程相互作用较为强烈的区域, 河流资源的过度开发利用以及污染物排放等压力使得城市河流健康状况不断恶化^[8], 而为了应对河流健康状况的恶化, 针对城市河流实施了大量的生态治理修复措施, 人类活动施加给河流的压力和采取的积极措施共同影响着城市河流的健康状态^[9]。然而, 目前针对城市河流构建的评价指标体系主要集中在水文、水质、生境等河流的自然属性以及人类对河流施加的压力方面^[10], 而忽视了人类采取的积极措施对河流生态健康的影响, 从而使得表征河流健康状态的信息考虑不够全面, 综合考虑人类正负面影响对河流健康状态的影响, 可提供更加全面准确的城市河流健康评价结果。

压力-状态-响应(PSR)模型于20世纪70年代由David J. Rapport和Tony Friend提出, 后由经济合作与发展组织(OECD)和联合国环境规划署(UNEP)推广发展^[11], 由于能够清晰的描绘外界压力、系统状态以及响应决策与行动三者之间的关系, 被广泛应用于生态安全、生态环境承载力等方面的评价^[12], 基于PSR模型筛选评价指标, 除了可以将不同地区城市河流健康评价体系的构建统一到同一框架下, 还能考虑人类积极措施对河流健康状态的影响, 有效提高城市河流健康综合评价准确性和合理性。彭涛等^[13]、Zhao等^[14]采用PSR模型构建城市河流健康评价指标体系, 利用集对分析方法进行了宜昌黄柏河和合肥十五里河的健康评价, 取得了较为全面准确的评价结果。

采用常规数学模型进行河流健康等级评定时, 其评价结果往往具有一定的模糊性^[15], 并且仅能针对当前指标监测数据进行现状评价, 难以根据评价结果判断河流健康状态的变化趋势。物元可拓模型最早由我国学者蔡文提出, 主要用于解决复杂的不相容问题^[16], 同时可以根据相邻等级综合关联度之间的距离判断河流健康水平的发展趋势, 将其与PSR模型相结合, 能够有效提升城市河流健康的评价准确性。城市河流往往流经城市的多个区域, 由于发展程度与功能定位的不同, 各区域内河段包含河流物理形态、土地利用类型以及河流水环境等要素在内的生态系统往往存在较大的差异^[17], 从而导致城市河流的健康状态呈现出明显的空间异质性, 基于河流沿程生态系统特征的差异, 进行河流分段健康评价, 可以有效识别各河段生态健康的关键制约因素。同时, 为避免河流尺度指标造成各河段之间差异的平均化, 影响评价结果, 选取的指标应集中在河段尺度和断面尺度。

本研究以城市河流鼓楼区西北护城河为对象, 根据河流形态特征及其所处区域功能定位差异划分河段, 基于PSR模型构建评价指标体系, 采用物元可拓模型进行河流健康评价并判断其变化趋势, 以期为城市河流的针对性生态治理与修复提供参考。

1 研究区域与数据 1.1 研究区域概况

西北护城河位于江苏省南京市鼓楼区中部, 全长约5.6 km, 汇水面积3.95 km²(图 1), 主要流经公园、景区以及居民住宅区等区域, 各河段所处环境存在显著的差异, 河流形态沿程呈现较大变化。为准确描述西北护城河各河段健康状况的差异, 考虑河流周边土地利用类型的区别以及河流形态的变化情况, 并结合城市地下管网走向, 将西北护城河划分为6个河段进行河流的分段健康综合评价。河段A起于小桃园泵站, 止于桃园桥, 全长0.55 km, 河宽10—15 m, 沿岸植被较多, 沿河无居民区；河段B起于桃园桥, 止于中山北路, 河段长1.3 km, 均宽100 m, 护岸为软质护坡, 功能定位为景观河道；河段C起于中山北路, 止于绣球公园闸, 为公园内湖泊, 水流流速缓慢, 人流量较大；河段D起于绣球公园闸, 止于卢龙桥, 河段长0.54 km, 均宽15 m, 河段蜿蜒曲折, 两岸主要为商业用地；河段E起于卢龙桥, 止于建宁西路桥, 河段长0.81 km, 均宽75 m, 右岸为景区, 左岸为道路；河段F起于建宁西路桥, 止于金川门泵站, 河段长1.24 km, 河段均宽35 m, 两岸均为居民住宅区, 采用硬质衬砌护岸, 河段内较多生态浮床及曝气设施。西北护城河的污染源主要为地表径流污染、生活污水以及本底污染, 水量补给多由闸坝控制, 其健康状态受人类活动影响大, 具有典型的城市河流特征^[18]。

1.2 数据资料

本研究于2019年对西北护城河6个河段布设12个采样点(图 1)进行每个季度水环境和水生态监测, 获取包括水质指标(COD_Mn、NH₃-N和TP等)、透明度、浮游植物生物量等水质水生态监测数据, 计算指标现状值时参考《地表水资源质量评价技术规程》^[19]中水质评价数据频次要求, 对四个季度的监测数据进行算术平均处理。城市河流污染负荷主要来自降雨携带、地表径流以及河道本底3个方面。降雨携带的年入河污染物负荷估算通过进行每个季度的降雨实验, 获得雨水中包括COD_Mn、NH₃-N以及TP在内的污染物浓度, 结合2019年水文年鉴中鼓楼区年降雨量数据, 估算污染物降雨入河量；地表径流污染物入河负荷估算通过各季度降雨试验测量雨水流经不同下垫面(绿地、硬质地面等)后靠近河道一侧的地表径流中各污染物浓度, 参考杨迪虎^[20]的研究确定不同下垫面面源污染物入河系数, 通过解译2019年9月高分1号遥感影像获得各河段汇水范围内不同类型下垫面面积, 估算污染物地表径流入河量；河道本底污染负荷根据水质监测数据及《地表水环境质量标准》^[21]进行估算；采取的生态修复措施及其工程量数据来自《西北护城河水质提升工程可行性研究报告》。

2 评价模型 2.1 评价指标体系

PSR模型包含人类活动给河流带来的压力、河流在人类压力作用下的状态以及人类为改善河流状态而采取的响应措施三个方面, 能够全面考虑人类对河流健康的正负面影响, 因此本研究从压力、状态、响应三个方面构建评价指标体系, 以进行全面准确的河流健康评价。

(1) 压力指标

随着社会对城市河流健康的关注与重视, 同时由于城市河流自净能力有限, 对于河流点源污染排放管控也越来越严格, 众多城市河流已进行了排污口的整治, 禁止污染物通过管道直接排入河流, 因此, 河流中影响水质水生态的污染物COD_Mn、NH₃-N以及TP主要来自面源冲刷、降雨携带以及河道本底累积的内源污染^[22], 考虑到各河段纳污量的差异, 采用各河段COD_Mn、NH₃-N以及TP三个污染物指标的年排放量与对应河段纳污量的比值表征污染负荷对河流水质水生态造成的压力。城市河流周边滨水建筑物的体量大小会直接影响到河岸带的宽度, 对河岸带的状态产生一定的胁迫, 而河岸带作为水生生态系统和陆地生态系统进行物质交换的关键区域, 能一定程度表征河流健康状况, 因此将滨水建筑物体量作为河岸带压力来源的指标。

(2) 状态指标

城市河流在各种人类活动压力的作用下, 其各方面的状态会受到一定的影响。其中污染物的排放量大小会直接影响到城市河流的水质健康状态, 使得城市河流的营养状态发生改变, 进而引起水生动植物的多样性的变化, 对河流的水质水生态产生影响, 水质状态选用透明度和水质等级两个指标表征, 水质等级的计算综合考虑COD_Mn、NH₃-N以及TP三个主要污染物参考《地表水环境质量标准》^[21]和《生态河湖状况评价规范》^[23]确定, 水生态状态采用浮游植物多样性指标表征；滨水建筑物体量会影响到河岸带的健康状态, 例如滨水建筑物体量较大的一般为商业住宅区, 河岸带宽度较窄, 而滨水建筑物体量较小的一般为公园绿地, 其河岸带宽度较宽, 因此采用在滨水建筑物作用下的河岸带宽度表征河岸带的状态。

(3) 响应指标

人类为了改善河流在压力作用下的状态而采取的响应措施对于河流的健康水平会产生一定的影响。城市河流相较于自然河流, 其河岸带的植被大都为人工种植和维护^[24], 同时部分河段修建了生态护坡以还原护岸的自然状态, 河岸带的植被种植情况及对应的岸坡自然状态共同影响着污染物的截留及水陆物质的交换, 因此将河岸带植被覆盖率和岸坡自然程度作为河岸带状态的响应指标；对于水质水生态状态而言, 往往会采取减少排污口和进行生态修复从而达到改善水质水生态的目的, 因此可以将排污口处理率及生态修复措施对污染物的削减情况作为水质水生态的响应指标。

根据上述分析选取指标, 基于PSR模型构建了包含压力、状态、响应三个方面共12指标的城市河流评价指标体系, 指标及其计算方式具体如表 1所示。

2.2 评价标准

目前对于城市河流的健康评价尚未制定统一的标准, 因此本研究结合国内外河流健康评价的研究成果以及城市河流的实际情况, 将目标层城市河流健康等级分为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ五个等级, 分别对应于健康、亚健康、中等、亚病态、病态五个健康水平^[26]。同时采取同样方式对准则层进行分级, 对于压力准则层, 分为低压力(Ⅰ)、较低压力(Ⅱ)、中等压力(Ⅲ)、较高压力(Ⅳ)以及高压力(Ⅴ)；对于状态准则层, 分为好状态(Ⅰ)、较好状态(Ⅱ)、中等状态(Ⅲ)、较差状态(Ⅳ)以及差状态(Ⅴ)；对于响应准则层, 分为高响应能力(Ⅰ)、较高响应能力(Ⅱ)、中等响应能力(Ⅲ)、较低响应能力(Ⅳ)以及低响应能力(Ⅴ)。各等级标准值参考相关文献、标准规范并结合西北护城河实际情况制定, 具体如表 2所示。

2.3 评价方法

当所给的条件不能达到要实现的目标时, 称为不相容问题^[30-31]。对于河流健康评价而言, 理论上指标选取的越多, 对于河流健康的描述就会更加全面, 但往往有些指标数据的获取会存在一定难度, 即当前所选的有限的指标并不能完全实现对河流健康状态的全面描述, 即河流健康评价是一个不相容问题。而物元可拓模型主要用于解决不相容的复杂问题^[15], 因此, 本文采用物元可拓模型进行河流健康等级评定, 并判断其演变趋势, 同时对河流准则层进行评价, 确定各准则层所处的等级。为减少权重判定的主观性影响, 采用熵权法进行权重计算。

(1) 河流健康物元矩阵、经典域及节域构建：物元是描述事物的名称、特征及量值3个基本元素的简称。将研究事物记为M, 本研究中指河流健康水平, M的特征记为C, 本研究中指各个评价指标, M关于C的量值记为V, 即为评价指标处于某一等级时对应的量值, 将有序三元组R=(M, C, V)称为物元矩阵。若研究事物M具有n个特征c₁, c₂, c₃, …, c_n 及其对应量值v₁, v₂, v₃, …, v_n, 相应物元矩阵R表示为：

(1)

经典域表示待评价对象的特征处于各个等级时对应的取值范围, 节域表示待评价事物对应特征的所有数据范围, 根据制定的评价标准即可得到相应经典域和节域矩阵^[32]。

(2) 关联函数计算：根据关联度计算获得指标所隶属的等级, K_j(v_i)表示第i个指标属于第j个等级的关联度, K_j(v_i) 越大, 则隶属于哪个等级, 计算公式如下：

(2)

其中,

(3)

式中, ρ(v_i, v_oji)和ρ(v_i, v_pi)分别为点v_i与经典域v_oji和节域v_pi的距离, v_i、v_oji、v_pi分别为待评价事物的量值、经典域的量值范围和节域的量值范围。

(3) 权重计算：熵权法根据指标自身所含信息量确定指标权重, 是相对客观的权重确定方法^[33], 采用熵权法计算得到第i个指标相对于目标层的权重w_i。

(4) 综合关联度计算：根据最大关联度原则, 确定研究事物评价等级, 待评价事物M_x(x=1, 2, 3, …, n)关于等级j的关联度K_j(M_x)为：

(4)

式中, m表示指标数量。若K_jx=max K_j(M_x), 则待评价事物M_x属于等级j。同时根据K_jx与K_(j-1)x及K_(j+1)x的距离大小, 判断当前位于等级j的待评价事物M_x的等级变化趋势。

3 结果分析 3.1 河段指标层现状值及对应等级

利用2019年数据资料, 根据表 1中各指标计算方法得到指标现状值, 结合表 2中的评价标准, 将各指标现状值代入式(2), 获得指标关于各个等级的关联度, 从而确定各指标所处的等级, 指标现状值及其所处等级如表 3所示。从表 3中可以看出, 尽管选取的城市河流的尺度较小, 但同一指标的等级在不同河段仍然存在明显的差异。各河段处于Ⅲ级以上(Ⅰ、Ⅱ)的指标数量占比分别为58.3%、66.7%、50%、25%、33%以及16.7%, 总体而言, 河段A, B、C各指标状态相对较好, 而河段D、F则相对较差。

3.2 指标权重分配

根据各河段指标现状值, 采用熵权法计算获得城市河流健康评价体系指标权重, 结果如表 4所示。从表 4中可以看出, 准则层权重表现为压力指标>状态指标>响应指标, 各指标相对于目标层的权重介于0.049—0.140, 权重较大的指标主要为COD_Mn、NH₃-N、TP排放量以及河流植被覆盖率, 权重均达到了0.12以上, 而其他指标的权重则相对较小, 根据熵权法的内涵, 表明COD_Mn、NH₃-N、TP排放量、河流植被覆盖率四个指标在不同河段现状值的离散程度大于其他指标, 即上述四个指标在不同河段的差异性相对更大, 对综合评价结果的影响程度也更大。

3.3 河段准则层及综合等级

将指标相对于准则层和目标层的权重以及关于各等级的关联度代入式(4), 计算得到各河段准则层以及河流健康水平关于各等级的关联度, 关联度结果如表 5所示, 各河段准则层等级及河流综合等级如图 2所示。从表 5和图 2中可以看出, 在河流健康水平综合等级方面, 河段A、B处于亚健康等级, 河段C、E处于中等等级, 而河段D、F则处于亚病态等级, 同时, 根据河流健康水平所处等级的关联度与其相邻水平关联度之间的距离, 可以看出河段C、D、E的健康水平均与更优等级的关联度距离相近, 表明该河段健康水平有进一步提升的趋势, 而河段A、E、F的健康水平则存在恶化的趋势, 同时, 各河段健康水平的变化趋势与该河段准则层的变化趋势也呈现出了较高的一致性。

4 讨论 4.1 河流健康水平的空间异质性

城市河流相较于自然河流, 除了需要保有部分自然生态功能外, 还要发挥一定的社会服务功能^[33-35], 其通常会流经城市的不同功能区域, 如公园绿地、居民区、商业区等, 这使得城市河流沿程周边的土地利用类型呈现出显著的差异性, 影响各区域河段的污染负荷输入量及滨水景观格局。同时, 为满足不同区域内河流相应服务功能的发挥, 人类往往会根据区域功能定位对河流的物理形态进行一定程度的人为干预, 如截弯取直、河道拓宽、岸坡衬砌^[29]等, 使得城市河流自身的形态结构沿程发生较大的变化, 这两方面原因使得城市河流在不同河段的生态系统特征呈现一定的差异, 进而导致河流健康水平的空间异质性。本研究中, 河段A、B作为城市景观河段, 主要处在公园绿地内, 河岸带自然程度较好, 能够截留一部分降雨径流冲刷过程中的污染负荷, 降低入河污染负荷量, 其健康水平相对较高; 而采用硬质护岸的河段(如河段D、F), 往往处于人类活动较为频繁的商业区或居民区, 存在较大的面源污染, 在降雨期间污染物易直接冲刷入河, 造成较大的外源污染, 使得河段压力水平较高, 难以维持较好的状态, 整体健康水平较低; 即使针对这些河段采取了一定的生态措施(曝气复氧、生态浮床等), 但不注意控制外部点面源污染的输入, 使得污染负荷入河量远高于生态修复措施对河段污染物的削减量(如河段F), 其健康水平仍然难以得到进一步改善, 导致显著的河段健康水平差异。

4.2 河流健康评价方法适用性

自然河流中, 不同生境的生物群落组成往往具有明显差异性^[36], 但城市河流在人类活动的干扰下, 其水文特征、渠化、底泥等环境因子的改变使得河流生境的复杂性和异质性降低, 不同生境中生物群落组成的相似性明显增加^[37], 虽然可以基于指示生物物种确定河流健康水平, 但难以辨别河流健康水平的空间差异性。因此, 本研究根据西北护城河周边的土地利用类型(例如公园绿地各居民小区的分界处)、河流物理形态(例如河宽即河流形状突变处)以及城市管网分布进行河段的划分；基于PSR框架筛选河段尺度指标构建了城市河流健康综合评价体系, 考虑了人类活动对河流健康状态的正负面影响, 指标分类明确, 能够较为全面的描述城市河流的综合健康状态, 适用于人类活动和自然过程相互作用强烈的城市河流。利用常规的模糊综合评价法^[38]、集对分析法^[39]以及灰色关联度分析法^[40]等方法确定河流健康等级时, 仅能获得单次评价结果, 无法进行河流未来健康状态的预测, 而物元可拓模型根据河流健康水平隶属于某一等级的综合关联度与其相邻等级综合关联度的距离, 除了可以确定当前河流所处的健康状态, 还可以预测下一阶段河流健康水平以及各项指标的发展趋势。为管理者下一阶段城市河流的针对性管理和保护提供科学的指导。

4.3 河段健康水平的变化趋势及改善措施

根据河段准则层及健康水平关联度与相邻水平关联度的距离, 可以看出河段B、C、D的健康水平有向更优等级发展的趋势, 而河段A、E、F的健康水平有恶化的趋势。尽管河段D有向更优等级发展的趋势, 但其现状健康水平接近病态, 宜针对河段D纳污量低、污染物易沉积的情况, 进行定期的河底清淤, 减少内源污染。河段A整体上有恶化的趋势, 但当前健康水平较高, 为亚健康水平, 保持或进一步降低当前压力输入水平及提高响应级别即可使河段A维持在理想的水平。合理的生态治理措施可以缓解河流的压力, 改善河流的状态, 最终达到提高河流健康水平的目的。从表 5中可以看出, 河段C、F当前压力均处于较高及以上水平, 但河段C由于采取了较高的响应措施, 河段状态水平维持在中等, 最终健康水平也为中等。因此, 针对河段F应重点排查沿河排污口设置情况, 坚决取缔非法排污口, 减小点源污染, 同时建立完善的雨污分流系统, 避免降雨径流过程将大量污染物冲刷入河, 增加河段压力。河段E的健康水平呈现出由当前中等水平向亚病态过渡的趋势, 虽然响应措施向更优发展, 但河段压力水平却在恶化, 结合准则层权重来看, 压力准则层权重最大, 对河流健康水平的影响也最大, 应加强生态治理措施, 设立合理的河道控制红线, 在控制红线内适当拓宽河岸带, 多设置人工生态护岸, 增加水-空气-植被-土壤之间的交互作用, 为动植物提供栖息生境, 增加生物多样性, 以更高的响应能力抵消压力增大对河流健康的影响, 改善河段F健康水平。同时, 严抓对应河流各级河长的考核工作, 保证河湖长制的长效性, 也是促进城市河流生态健康改善与恢复的重要措施。

5 结论

(1) 根据PSR模型构建的西北护城河健康评价体系包含压力指标、状态指标和响应指标3个方面, 涵盖COD_Mn排放量、岸坡自然程度、水质等级等12个指标, 3类指标的权重排序为压力指标>状态指标>响应指标。

(2) 西北护城河6个河段的健康水平存在明显的差异, 河段A、B为亚健康水平、河段C、E为中等水平, 河段D、F为亚病态水平；各河段生态系统组成要素的区别与其健康水平的差异有较大的关联性, 流经公园绿地的河段健康水平较高, 流经居民区的河段健康水平相对较低, 从趋势变化上来看, 河段B、C、D的健康水平有向更优等级发展的趋势, 而河段A、E、F的健康水平有恶化的趋势。

(3) 本研究基于PSR模型构建评价指标体系, 结合熵权法和物元可拓模型, 能够对不同河段的外界压力程度、系统状态以及响应决策能力进行分级评估, 考虑了人类修复措施对河流健康水平的影响, 能较为准确全面的评价各河段当前健康水平并判断其变化趋势。