[发明专利]一种清洁河流的断面综合评价方法

权利要求书

1.一种清洁河流的断面综合评价方法，包括以下步骤：

(1)构建清洁河流断面综合评价要素及指标：

清洁河流断面评价要素包括水文水质特征Y₁、生物特征Y₂、功能特征Y₃与地貌特征Y₄；

水文水质特征Y₁包括流速X₁、单位m/s；

河流生态需水保证率X₂，记为W_E，单位％,计算公式为：

式(1)中W_E为生态需水保证率，M_E为流量大于或等于生态环境需水流量的天数，M_T为总天数；

浊度X₃、单位NTU；

溶解氧X₄、单位mg/L；

化学需氧量X₅、单位mg/L；

氨氮X₆、单位mg/L；

总磷X₇、单位mg/L；

铅X₈、单位mg/L；

六价铬X₉、单位mg/L；

有机质X₁₀、单位mg/kg；

底泥潜在生态危害指数X₁₁，记为RI，计算公式为：

RI＝∑E_rⁱ＝∑T_rⁱ×C_rⁱ 式(2)；

式(2)中：RI表示底泥潜在生态危害指数，E_rⁱ表示某单个重金属i的潜在生态危害指数；T_rⁱ为第i种重金属的毒性系数，反映重金属的毒性水平及水体对重金属污染的敏感程度；C_rⁱ(C_rⁱ＝C_sⁱ/C_nⁱ)表示某单个重金属i的污染指数，C_sⁱ为底泥重金属i的浓度实测值，C_nⁱ为底泥重金属i的背景参比值；

X₁～X₁₁构成水文水质特征评价指标；

生物特征Y₂包括浮游动物香农多样性指数X₁₂、浮游植物香农多样性指数X₁₃、固着藻类香农多样性指数X₁₄，底栖动物种数X₁₅，X₁₂、X₁₃、X₁₄的香农多样性指数记为H,计算公式为：

H＝-∑p_ilnp_i 式(3)；

式(3)中：p_i表示第i种物种个数占总物种个数的比例；

X₁₂～X₁₅构成生物特征评价指标；

功能特征Y₃包括景观效应X₁₆，面源污染强度X₁₇、单位kg/km²，污水处理率X₁₈、单位％，污水处理率X₁₈计算公式为：

X₁₆～X₁₈构成功能特征评价指标；

地貌特征Y₄包括纵向连通性指标X₁₉和河岸带稳定性指标X₂₀，纵向连通性指标X₁₉的计算公式为：

X₁₉～X₂₀构成地貌特征评价指标；

(2)对河流断面的现场采样和监测：

利用分析仪器进行现场监测、对河流断面进行物理化学指标测定和生物鉴定，物理指标包括流速、河宽、水深，化学指标包括水体污染物指标和底泥污染物指标，水体污染物指标为浊度、溶解氧、化学需氧量、氨氮、总磷、铅、六价铬；底泥污染物指标包括Cu、Pb、Zn、Cr、Cd、Hg、As和有机质，表征地貌数据的纵向连通性采用ArcGIS软件进行数据提取计算获得，用水样采集器、底泥采样器、激光测距仪、卷尺、便携式流速仪、便携式溶解氧测定仪、多参数水质分析仪、便携式冰箱、电子称对河流现场采样，对河流现场采集的样品用原子吸收分光光度计、原子荧光光度计、电感耦合等离子体发射光谱仪进行水质分析，对浮游动物、浮游植物、固着藻类采用生物显微镜鉴定，根据现场测定结果，计算得到清洁河流断面评价指标的实测值；

(3)数据的标准化处理：

由于不同评价指标间的量纲不同，为了消除量纲，提高不同指标数据的可比性，需要对n个评价断面，20个评价指标的实测值r_ij进行标准化处理；根据各评价指标的内涵，评价指标分为“效益型”和“成本型”两类指标，对于“效益型”指标属性值越大越好；而“成本型”指标为属性值越小越好，为了对评价指标建立惩罚机制，采用极差标准化法对各指标数据进行标准化处理，对于效益型和成本型两类指标实测值的标准化结果均用x_ij表示，具体计算公式如下：

对于效益型评价指标，实测值标准化计算公式为：

对于成本型评价指标，实测值标准化计算公式为：

式(6)、式(7)中max{r_ij}、min{r_ij}分别为指标实测值的最大值和最小值；根据以上的标准化结果可知，对于效益型指标，r_ij越大，水质越好，x_ij越小，其指标情况越好；对于“成本型”评价指标，r_ij越小水质越好，x_ij越小，其指标情况越好；

(4)评价指标数据转换矩阵的建立：

河流有n个评价断面O₁,O₂,…,O_n，20个评价指标X₁,X₂,…,X₂₀构成评价指标系统，x_ij(i＝1,2,…,n；j＝1,2,…,20)为被评价断面O_i关于指标X_j实测值的规范化值，记为N＝{1,2,…,n}(3≤n≤100),M＝{1,2,…,20}，评价断面O_i的指标值集合为Z_i＝{x_i1,x_i2,…,x_i20}，将清洁河流断面各评价指标分级限值作为标准断面处理，建立标准评价断面数据矩阵，对清洁河流断面分为不清洁、较不清洁、较清洁、清洁、很清洁5级；

标准化处理的评价矩阵为x_Ⅰj、x_Ⅱj、x_Ⅲj、x_Ⅳj 4个清洁河流标准断面O_Ⅰ、O_Ⅱ、O_Ⅲ、O_Ⅳ的Ⅰ、Ⅱ、Ⅲ、Ⅳ类标准上限值标准化处理后的数据，则标准断面O_Ⅰ的指标值集合为Z_I＝{x_I1,x_I2,...,x_I20}，记为N^*＝{Ⅰ、Ⅱ、Ⅲ、Ⅳ}，因此，N扩展为(N∪N^*)；

评价指标数据转换矩阵A为：

其中，x_ij(i＝1,2,…,n,Ⅰ，…，Ⅳ；j＝1,2,…,20)；

(5)自主优势量矩阵的建立：

根据式(8)：

式(9)：

λ_ij＝μα_ij+υβ_ij i∈(N∪N^*),j∈M 式(9)

其中：μ为竞争性目标偏爱系数，υ为发展性目标惩罚系数，计算出自主优势量λ_ij(i∈(N∪N^*),j∈M)，得到自主优势量矩阵B＝[λ_ij]_(n+4)×20；

λ_ij为评价断面O_i(i∈(N∪N^*))关于指标X_j(j∈M)的自主优势量，其值越大说明评价断面自主优势量越大，μ,υ∈[0,1],μ+υ＝1；k表示第k个列优势评价断面，p表示第k个列优势评价断面下第p个指标；α_ij,β_ij为评价断面O_i(i∈(N∪N^*))在指标X_j(j∈M)上的列优势量及行优势量，列优势量α_ij(i∈(N∪N^*),j∈M)反映了评价断面O_i的第j项指标与其他(n+4)-1个评价断面整体之间改善的差异，行优势量β_ij(i∈(N∪N^*),j∈M)反映了评价断面O_i的第j项指标与其他m-1项指标整体的优势状况差异，评价指标m，1≤m≤20；

自主优势量矩阵B＝[λ_ij]_(n+4)×20，评价断面O_i的自主优势量向量为λ_i＝(λ_i1,λ_i2,…,λ_i20)(i∈(N∪N^*))，B是根据式(8)与式(9)对评价指标数据转换矩阵A上的信息转换而成，其λ_ij(i∈(N∪N^*),j∈M)表征出评价断面O_i在指标X_j上的优势程度，因而，对应的指标分别是评价断面O_i的最优势及最劣势指标；

(6)计算得到指标优势序组：

对自主优势量向量λ_i＝(λ_i1,λ_i2,…,λ_i20)(i∈(N∪N^*))中各分量由大到小降序排列，得到一个有序组，将{1,2,…,n,Ⅰ,…，Ⅳ}上的元素指定给有序(n+4)元组上的对应元素，记λ_ij对应的序值为d_ij，称d_ij为评价断面O_i在指标X_j上的优势数；

对评价断面O_i的优势数集合{d_ij|j∈M}中元素按由小到大升序排序，得到评价断面O_i的指标优势序组，记为D_i,i∈(N∪N^*)；

(7)计算评价指标的最终优势权向量：

最终优势权向量计算公式为：

其中，u为连续递增函数；1≤m≤20；

b(b≥0)为可变参数；

当b＝0时，ω_ij＝1/m，此时，ω_i完全均衡各指标的作用，指标自身优势削弱程度最弱；当b→+∞时，ω_i1→1；此时，ω_i2,ω_i3,…,ω_im→0，ω_i1最大化削弱了评价断面O_i的自身指标优势，优势弱化能力达到最大；

为关于指标X_j(j∈M)的优势权向量，其中i∈(N∪N^*)；将中分量按照被评价断面O_i优势数d_ij升序排列的结果，即为指标优势序组的优势权向量；

清洁河流断面评价指标分为多个评价要素体系，需要计算重新排列后的评价要素优势权向量与相应的评价要素下的指标优势权向量并得到评价指标的最终优势权向量

(8)计算评价断面O_i综合评价值；

清洁河流综合评价值为：

其中，λ_ij为评价断面O_i关于评价指标X_j的自主优势量，为评价断面O_i关于评价指标X_j的优势权向量，

(9)清洁河流评价断面的优劣排序：

通过比较综合评价值y_i(i∈(N∪N^*))得到各评价断面O_i的相对优劣；根据n个评价断面的y_i值与标准断面y_Ⅰ、y_Ⅱ、y_Ⅲ、y_Ⅳ值的相对情况，对n个评价断面O₁,O₂,…,O_n进行评价等级划分，得到清洁河流的评价断面O_i的优劣排序，综合评价值y_i(i∈(N∪N^*))越小，其综合评价结果越好，从而实现对清洁河流的断面综合评价。