system: [SYSTEM]

user: 2. 研究方法 2.1. 评价体系 2.1.1. 评价体系指标 建立河流健康评价指标体系的基础是科学地通过选择适当指标，反映和衡量河流所处的状态，从而为寻求河流的可持续发展提供方向。考虑到评价指标的跨学科跨领域和种类、项目繁多的特点 [ 13 ]，本文筛选了4类代表性的指标，分别反映河流的水文水资源、水质、水生态和社会服务方面，旨在尽可能涵盖河流完整的功能。 1) 水文水资源 为了维持河流生态系统的生态环境功能，河流必须要能够提供足够的水量。本文选取水文学中生态流量的指标，计算其满足程度来表征河流径流量满足其生态系统需求的程度。生态流量满足度是河流流量能够满足生态流量的天数与总天数的比值，计算公式如下： α i j = D i j D = ∑ sgn ( Q i j k − Q i j ) D (1) sgn ( Q i j k − Q i j ) = { 1       Q i j k > Q i j 0       Q i j k < Q i j (2) 式中，α ij 表示第i年第j月的生态流量满足度；D ij 表示第i年第j月生态流量的满足天数(d)；D表示第i年第j月的总天数(d)；Q ijk 表示第i年第j月第k日的河道流量(m³/s)；Q ij 表示第i年第j月的生态流量(m³/s)。 2) 水质 保障河流健康不仅需要提供足够的水量，也需要保证合格的水质。本文根据《地表水环境质量标准》(GB 3838-2002)，选取了以下四个评价项目：DO水质状况、有机污染水质状况、无机污染水质状况和重金属污染状况。其中有机污染水质状况包含高锰酸盐指数、化学需氧量、五日生化需氧量、氨氮4个水质指标，无机污染水质状况主要考虑总磷指标，重金属污染状况包含铜、锌、硒、砷、汞、镉、六价铬、铅8个水质指标。 对于只包含一个指标的DO水质状况和无机污染水质状况，其项目值采用相应水质指标评估年的平均浓度。对于包含多个指标的有机污染水质状况和重金属污染状况，将它们各自包含的水质指标分别进行赋分，见表1和表2。根据监测数据，每个水质指标将其全年12个月的月均浓度按照汛期和非汛期分别进行平均，再分别确定汛期与非汛期的赋分，取其最低赋分作为此水质指标的赋分，最后取所包含的各个指标的赋分平均值作为该项目的赋分 [ 13 ]。 Table 1 The score distribution of organic pollutio 水质指标 指标浓度(mg/L) 高锰酸盐指数 ≤2 ≤4 ≤6 ≤10 ≤15 化学需氧量 ≤10 ≤15 ≤20 ≤30 ≤40 五日生化需氧量 ≤2 ≤3 ≤4 ≤6 ≤10 氨氮 ≤0.15 ≤0.5 ≤1 ≤1.5 ≤2 赋分 100 75 50 25 0 表1. 有机污染水质状况水质指标赋分 Table 2 The score distribution of heavy metal pollutio 水质指标 指标浓度(mg/L) 铜 ≤0.01 ≤0.2575 ≤0.505 ≤0.7525 ≤1 锌 ≤0.05 ≤0.5375 ≤1.025 ≤1.5125 ≤2 硒 ≤0.01 ≤0.0125 ≤0.015 ≤0.0175 ≤0.02 砷 ≤0.05 ≤0.0625 ≤0.075 ≤0.0875 ≤0.1 汞 ≤0.0005 ≤0.000625 ≤0.00075 ≤0.000875 ≤0.001 镉 ≤0.001 ≤0.00325 ≤0.0055 ≤0.00775 ≤0.01 铬(六价) ≤0.01 ≤0.0325 ≤0.055 ≤0.0775 ≤0.1 铅 ≤0.01 ≤0.0325 ≤0.055 ≤0.0775 ≤0.1 赋分 100 75 50 25 0 表2. 重金属污染状况水质指标赋分 3) 水生态 水生态健康调查评价中，大多数指标体系只考虑了水生动物，而评价一条河流时，这样的考虑过于笼统和片面。因此本文将水生动物分为底栖动物和浮游动物，另外将浮游植物也考虑其中，采用Shannon Wiener指数来评价它们的多样性 [ 14 ]，计算公式如下： H = − ∑ i = 1 s ( p i ) ( ln p i ) (3) 式中，S表示总种数；p i 表示第i种类的个体数占总个体数的比例。 4) 社会服务 河流健康是自然价值与人类社会服务价值的统一体，健康不仅意味着自然功能的完整性，还要确保其社会功能的正常发挥 [ 12 ]。以往的研究往往只从水文、水质和水生态的角度去分析河流健康，忽略了水的“社会–自然”二元特性。人类活动对河流健康的影响越来越大，经济需水和生态需水的矛盾越来越突出 [ 15 ]。本文充分考虑河流完整的功能，在评价指标中加入体现社会服务功能的通航水位保证率和水资源开发利用率两个指标。 河道航运功能的保证能促进货物的顺利流通，有利于社会经济的发展。河道航运功能的保证采用通航水位保证率来衡量，计算公式如下： β i j = D i j D = ∑ sgn ( Z i j k − Z i j ) D (4) sgn ( Z i j k − Z i j ) = { 1       Z i j k > Z i j 0       Z i j k < Z i j (5) 式中，β ij 表示第i年第j月的通航水位保证率；D ij 表示第i年第j月通航水位得到满足的天数(d)；D表示第i年第j月的总天数(d)；Z ijk 表示第i年第j月第k日的河道水位(m)；Z ij 表示第i年第j月的通航水位(m)。 水资源的开发利用是改造和利用自然的一个方面，其目的是发展社会经济，产生社会效益。水资源开发利用的内容很广，包括农业灌溉、工业用水和生活用水等等。本文用水资源开发利用率来衡量水资源开发利用程度，计算公式如下： R u = W u W t (6) 式中，R u 为水资源开发利用率；W u 为已开发利用的水资源总量；W t 为流域的水资源总量。 2.1.2. 评价体系建立 结合上述河流四个方面的特性，选取的河流健康综合评价体系见表3。 Table 3 River health comprehensive evaluation syste 准则层 指标层 计算说明 水文水资源 生态流量满足度 计算评估年的日径流能够满足选取断面生态流量的程度。 水质 DO水质状况 根据评估年的监测数据，对12个月的溶解氧浓度取平均。 有机污染水质状况 计算方法详见2.1.1。 无机污染水质状况 根据评估年的监测数据，对12个月的总磷浓度取平均。 重金属污染状况 计算方法详见2.1.1。 水生态 底栖动物生物多样性 采用底栖动物Shannon Wiener指数。 浮游动物生物多样性 采用浮游动物Shannon Wiener指数。 浮游植物生物多样性 采用浮游植物Shannon Wiener指数。 社会服务 通航水位保证率 统计评估年日水位大于等于通航水位的天数百分比。 水资源开发利用率 评估年水资源开发利用量／多年平均水资源量×100%。 表3. 河流健康评价指标体系 2.1.3. 评价标准 评价指标的标准参照国际及国家已经正式发布的标准和通过实际调查与理论分析综合确定的指标评价标准值 [ 13 ]。结合各指标特点，对每项指标划定5级，分别为I (理想状况)、II (健康)、III (亚健康)、IV (不健康)、V (病态)。指标层各项指标评价标准见表4。 Table 4 The criterion of each selected indicato 准则层 指标层 I II III IV V 水文水资源 生态流量满足度 80 60 40 20 0 水质 DO水质状况 7.5 6 5 3 2 有机污染水质状况 100 75 50 25 0 无机污染水质状况 0.02 0.1 0.2 0.3 0.4 重金属污染状况 100 75 50 25 0 水生态 底栖动物生物多样性 4 3 2 1 0 浮游动物生物多样性 4 3 2 1 0 浮游植物生物多样性 4 3 2 1 0 社会服务 通航水位保证率 95 85 70 60 0 水资源开发利用率 10 30 40 50 60 表4. 指标层各项指标评价标准 2.2. 模糊物元模型 2.2.1. 模糊物元 给定事物的名称N，它关于特征C有量值为V，以有序三元R(N，C，V)组作为描述事物的基本元，简称物元。如果物元模型中量值V具有模糊性，便称它为模糊物元。如果m个事物的n维物元组合在一起，使构成m个事物n维复合物元R nm 。若将R nm 的量值改写为模糊物元量值，称为m个事物n维复合模糊物元，记作 R n m = [ M 1 ⋯ M m C 1 v 11 ⋯ v 1 m ⋮ ⋮ ⋮ C n v n 1 ⋯ v n m ] (7) 式中，M j 为第j个事物( j = 1 ， 2 ， ⋯ ， m )；C i 为第i个特征( i = 1 ， 2 ， ⋯ ， n )；v ij 为第j个事物第i个特征对应的模糊量值。 2.2.2. 从优隶属度模糊物元 各单项评价指标相应的模糊量值从属于标准方案中各对应评价指标相应的模糊量值的隶属程度，称为从优隶属度。根据从优隶属度计算公式对所有指标的原始值标准化，以消除量纲的影响。 由此可得到从优隶属度模糊物元 R ˜ n m ： R ˜ n m = [ M 1 ⋯ M m C 1 x 11 ⋯ x 1 m ⋮ ⋮ ⋮ C n x n 1 ⋯ x n m ] (8) 2.3. 熵权法确定指标权重 本文采用熵权法来确定河流健康评价中各个评价指标的权重。熵权法源于信息科学中的信息熵，其值用于反映信息的无序化程度，进而可以度量信息量的大小。若某项指标所含有的信息量越大，则表示该指标在评价体系中的作用也就越大。此方法可以避免权重赋值过程中主观性的影响。