针对大气环境质量评价中广泛存在的指标测量值精度不高的特点,给出基于改进的层次分析法确定指标权重的模糊综合评价方法,并对天津市环境空气质量进行综合评价。该方法借助单指标污染情况,参照环境质量等级基准值,建立三角模糊数形式的模糊判断矩阵。基于改进的模糊层次分析法确定指标权重的模糊综合评价方法避免了传统模糊层次分析对专家主观判断的依赖,且所得结果较为合理。 Aiming at the low precision of the index measurement value widely existing in the air environmental quality evaluation, a fuzzy comprehensive evaluation method based on the improved analytic hierarchy process to determine the index weight is given, and a comprehensive evaluation of the environmental air quality in Tianjin is carried out. This method builds a fuzzy judgment matrix in the form of triangular fuzzy numbers with the help of single-indicator pollution conditions and reference values of environmental quality grades. The fuzzy comprehensive evaluation method based on the improved fuzzy analytic hierarchy process to determine the index weight avoids the dependence of the traditional fuzzy analytic hierarchy process on the subjective judgment of experts, and the results obtained are more reasonable.
针对大气环境质量评价中广泛存在的指标测量值精度不高的特点,给出基于改进的层次分析法确定指标权重的模糊综合评价方法,并对天津市环境空气质量进行综合评价。该方法借助单指标污染情况,参照环境质量等级基准值,建立三角模糊数形式的模糊判断矩阵。基于改进的模糊层次分析法确定指标权重的模糊综合评价方法避免了传统模糊层次分析对专家主观判断的依赖,且所得结果较为合理。
环境质量评估,模糊层次分析法,三角模糊数,污染指数
Zhenhua Li, Chun Ding
School of Mathematics and Physics, North China Electric Power University, Beijing
Received: Nov. 17th, 2021; accepted: Dec. 21st, 2021; published: Dec. 28th, 2021
Aiming at the low precision of the index measurement value widely existing in the air environmental quality evaluation, a fuzzy comprehensive evaluation method based on the improved analytic hierarchy process to determine the index weight is given, and a comprehensive evaluation of the environmental air quality in Tianjin is carried out. This method builds a fuzzy judgment matrix in the form of triangular fuzzy numbers with the help of single-indicator pollution conditions and reference values of environmental quality grades. The fuzzy comprehensive evaluation method based on the improved fuzzy analytic hierarchy process to determine the index weight avoids the dependence of the traditional fuzzy analytic hierarchy process on the subjective judgment of experts, and the results obtained are more reasonable.
Keywords:Environmental Quality Assessment, Fuzzy Analytic Hierarchy Process, Triangular Fuzzy Number, Pollution Index
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
环境质量评价是识别环境质量问题并制定合理修复方案进行环境治理的基础。目前常见的环境质量评价方法有单因子评价法 [
设考虑n个评价指标 f 1 , f 2 , ⋯ , f n ,对m个评价对象 u 1 , u 2 , ⋯ , u n 的环境质量等级进行评价,设 c i j 是第i个评价对象在第j个指标上的测量值, i = 1 , 2 , ⋯ , m ; j = 1 , 2 , ⋯ , n 。根据《环境空气质量标准》(GB3095-2012)和《环境空气质量评价技术规范(试行)》(HJ663-2013)的规定,并考虑实际数据获取情况,本文选取大气环境质量评价指标集F = {PM2.5,PM10,SO2,NO2,CO,O3}。结合《环境空气质量标准》(GB3095-2012)和天津市各区的实际污染状况,本文将评价等级分为优(I)、良(II)、轻度污染(III)、中度污染(IV)和重度污染(V)等五个等级,并给出大气环境质量指数分级标准,见表1。
大气环境评价指标 | I | II | III | IV | V |
---|---|---|---|---|---|
PM2.5 | 15 | 35 | 75 | 120 | 165 |
PM10 | 40 | 70 | 140 | 210 | 320 |
SO2 | 20 | 60 | 100 | 140 | 180 |
NO2 | 20 | 40 | 60 | 80 | 100 |
CO | 2 | 4 | 6 | 8 | 10 |
O3 | 100 | 160 | 220 | 280 | 340 |
表1. 环境空气质量分级标准
注:除CO质量浓度为mg/m3,其余评价指标单位为斗ug/m3。
定义1 [
μ M = { x m − l − l m − l , l ≤ x ≤ m ; x m − u − u m − u , m < x ≤ u ; 0 , else . (1)
式中, l < m < u ,l和u分别是M的下界和上界,m为M的中值即最可能值。
设三角模糊数 M 1 = ( l 1 , m 1 , u 1 ) , M 2 = ( l 2 , m 2 , u 2 ) ,其运算如下:
1) M 1 + M 2 = ( l 1 + l 2 , m 1 + m 2 , u 1 + u 2 )
2) M 1 − M 2 = ( l 1 − l 2 , m 1 − m 2 , u 1 − u 2 )
3) M 1 M 2 = ( l 1 l 2 , m 1 m 2 , u 1 u 2 )
4) 1 M 1 = ( 1 u 1 , 1 m 1 , 1 l 1 )
针对大气环境质量评价问题,本文给出一种基于改进的模糊层次分析确定指标权重的模糊综合评价方法。具体实现过程可描述如下。
下面给出建立模糊判断矩阵的方法,其具体步骤包括:
Step 1单指标污染情况建模。为分析各评价指标相应污染物的污染情况,我们以各评价等级的基准平均值为参照,定义污染指数,即 d j = c j / s j , j = 1 , 2 , ⋯ , n 。其中 s j 为第j个指标的各级标准值的算数平均值, c j 是第j个指标的实测值。它反映评价对象在该指标下的环境质量水平。 d j 越大表明环境中该指标相应的主要污染物含量越高,其环境危害也就越严重。相应的,其在环境质量评价中的相对重要性也就应该越高。
Step 2确定判断矩阵中三角模糊数的峰值点。由单项指标的污染指数的比值确定三角模糊 a i j = ( a l i j , a m i j , a u i j ) 的中值 a m i j ,即
a m i j = d i / d j , i , j = 1 , 2 , ⋯ , n (2)
a m i j 越大表示指标 f i 相对于 f j 的重要程度越高。
Step 3确定三角模糊数的上下边界。设n个评价指标的平均污染水平为 d ¯ = 1 n ∑ j = 1 n d j 。基于该均值可
确定相对重要性的下界和上界
a l i j = max { 0 , d i / d j − d ¯ / d i } , a u i j = d i / d j + d ¯ / d i (3)
Step 4确定模糊判断矩阵。根据以上步骤, f i 与 f j 的相对重要性可建模为三角模糊数 a i j = ( a l i j , a m i j , a u i j ) , i , j = 1 , 2 , ⋯ , n 。对n个评价指标 f 1 , f 2 , ⋯ , f n 的相对重要性进行两两比较,得到三角模
糊数判断矩阵 A = ( a i j ) n × n 。
给出改进的模糊层析分析法确定指标权重的计算过程:
Step 1对模糊判断矩阵 A = ( a i j ) n × n 作行归一化处理,
a i ¯ = ∑ j = 1 n a i j / ∑ i = 1 n ∑ j = 1 n a i j , i = 1 , 2 , ⋯ , n (4)
Step 2计算三角模糊数排序的可能度矩阵。根据文献 [
度,建立可能度矩阵 ( p i j ) n × n 。
Step 3计算各指标的权重
W = ( w 1 , w 2 , ⋯ , w n ) = ( ∑ j = 1 n p 1 j / ∑ i = 1 n ∑ j = 1 n p i j , ∑ j = 1 n p 2 j / ∑ i = 1 n ∑ j = 1 n p i j , ⋯ , ∑ j = 1 n p n j / ∑ i = 1 n ∑ j = 1 n p i j ) (5)
综上,可给出基于改进模糊层次分析法确定指标权重的模糊综合评价方法具体实施步骤:
Step 1构建评价因素集。根据评价标准及实际情况,进行评价指标的选取,建立指标集F
F = { f 1 , f 2 , ⋯ , f n } = { PM 2 .5 , PM 10 , SO 2 , NO 2 , CO , O 3 } (6)
Step 2建立评价集。依据各评价指标的取值范围结合专家经验,建立评语集V
V = { v 1 , v 2 , v 3 , v 4 , v 5 } = { Ι , Ι Ι , Ι Ι Ι , Ι V , V } (7)
Step 3利用改进的模糊层次分析方法确定指标权重向量W。
Step 4确定各指标的等级隶属函数 r j k , j = 1 , 2 , ⋯ , n ; k = 1 , 2 , 3 , 4 , 5 。本文选用梯形分布隶属函数 [
r j 1 = { 1 , c j ≤ μ j 1 ; μ j 2 − c j μ j 2 − μ j 1 , μ j 1 < c j ≤ μ j 2 ; 0 , μ j 2 < c j . (8)
当质量级别 2 ≤ k ≤ 4 时
r j k = { c j − μ j ( k − 1 ) μ j k − μ j ( k − 1 ) , μ j ( k − 1 ) < c j ≤ μ j k ; μ j ( k + 1 ) − c j μ j ( k + 1 ) − μ j k , μ j k < c j ≤ μ j ( k + 1 ) ; 0 , c j ≤ μ j ( k − 1 ) or μ j ( k + 1 ) < c j . (9)
当质量级别 k = 5 时
r j 5 = { 0 , c j ≤ μ j 4 ; c j − μ j 4 μ j 5 − μ j 4 , μ j 4 < c j ≤ μ j 5 ; 1 , μ j 5 < c j . (10)
其中, c j 为第j个评价指标实际测量的值, μ j k 为第j个评价指标的5个评价等级的分界线(见表1)。确定隶属矩阵 R = ( r j k ) n × 5 。
Step 5模糊综合评价。常用的模糊合成运算的算子有 M ( ∧ , ∨ ) , M ( • , ⊕ ) ,其中 ∧ 为取小运算, ∨ 为取大运算, • 为乘法运算, ⊕ 为求和运算,本文采用 ( • , ⊕ ) 算子进行矩阵合成运算,并按最大隶属度原则确定评价等级。
B = W ∘ R = ( b 1 , b 2 , b 3 , b 4 , b 5 ) (11)
天津市主要大气污染物为NO2、PM10、SO、PM2.5、CO、O3,根据中国空气质量在线监测分析平台获取2020年天津月统计历史数据,如表2所示,对天津市2020年的空气质量进行评价,环境空气质量评价标准见表1。
月份 | PM2.5 | PM10 | SO2 | NO2 | CO | O3 |
---|---|---|---|---|---|---|
2020-1 | 101 | 105 | 12 | 60 | 1.448 | 50 |
2020-2 | 61 | 62 | 9 | 33 | 0.931 | 68 |
2020-3 | 43 | 67 | 8 | 37 | 0.7 | 89 |
2020-4 | 40 | 79 | 9 | 34 | 0.68 | 123 |
2020-5 | 35 | 63 | 7 | 32 | 0.732 | 129 |
---|---|---|---|---|---|---|
2020-6 | 38 | 76 | 8 | 29 | 0.83 | 171 |
2020-7 | 41 | 53 | 5 | 23 | 0.852 | 163 |
2020-8 | 34 | 45 | 5 | 26 | 0.942 | 136 |
2020-9 | 29 | 44 | 7 | 35 | 0.86 | 113 |
2020-10 | 54 | 80 | 10 | 52 | 0.852 | 78 |
2020-11 | 47 | 76 | 10 | 53 | 0.89 | 47 |
2020-12 | 53 | 78 | 12 | 56 | 0.984 | 38 |
表2. 天津市2020年1~12月统计历史数据
注:除CO质量浓度为mg/m3,其余评价指标单位为ug/m3。
根据评价指标间相对污染程度的比较,可得到三角模糊数判断矩阵,进而利用改进的模糊层次分析法确定各指标的权重。以2020年5月为例,其基于污染指数的相对重要性建模可得到的三角模糊数判断矩阵如下所示。
( ( 0.1635 , 1 , 1.8365 ) ( 0.2204 , 1.0569 , 1.8935 ) ( 5.2610 , 6.0976 , 6.9341 ) ( 0 , 0.8003 , 1.6369 ) ( 2.6621 , 3.4986 , 4.3351 ) ( 0 , 0.7279 , 1.5645 ) ( 0.0620 , 0.9462 , 1.8303 ) ( 0.1158 , 1 , 1.8842 ) ( 4.8851 , 5.7692 , 6.6534 ) ( 0 , 0.7572 , 1.6414 ) ( 2.4261 , 3.3102 , 4.1944 ) ( 0 , 0.6887 , 1.5729 ) ( 0 , 0.1640 , 5.2649 ) ( 0 , 0.1733 , 5.2742 ) ( 0 , 1 , 6.1009 ) ( 0 , 0.1313 , 5.2321 ) ( 0 , 0.5738 , 5.6747 ) ( 0 , 0.1194 , 5.2203 ) ( 0.5800 , 1.2495 , 1.9190 ) ( 0.6511 , 1.3206 , 1.9901 ) ( 6.9496 , 7.6190 , 8.2885 ) ( 0.3305 , 1 , 1.6695 ) ( 3.7021 , 4.3716 , 5.0411 ) ( 0.2401 , 0.9096 , 1.5791 ) ( 0 , 0.2858 , 3.2126 ) ( 0 , 0.3021 , 3.2288 ) ( 0 , 1.7429 , 4.6696 ) ( 0 , 0.2288 , 3.1555 ) ( 0 , 1 , 3.9267 ) ( 0 , 0.2081 , 3.1348 ) ( 0.7648 , 1.3738 , 1.9827 ) ( 0.8430 , 1.4519 , 2.0609 ) ( 7.7677 , 8.3766 , 8.9856 ) ( 0.4905 , 1.0994 , 1.7084 ) ( 4.1973 , 4.8063 , 5.4152 ) ( 0.3911 , 1 , 1.6089 ) )
借助改进的模糊层次分析方法,可得到天津市2020年5月份各指标的权重。
W = ( 0.1782 , 0.1203 , 0.1153 , 0.2416 , 0.0453 , 0.2994 )
重复上述过程,可获得的天津市2020年其他月份的指标权重。
根据4.3节给出的隶属函数确定模糊隶属矩阵。同样以2020年5月为例,获得的模糊评价矩阵如下
R = ( 0 0.2333 1 0.7667 0 0 0 0 0 0 1 0 0 0 0 0.4000 1 0.6000 0 0 0 0 0 0 0 0.5167 0.4833 0 0 0 )
同理可得天津市其他月份的模糊评价矩阵。
通过获得的指标权重和模糊评价矩阵,可计算隶属函数,进而可以得到天津市2020年空气质量评价等级,见表3。由表3可以看出改进模糊层次分析法与熵权法、中国空气质量在线监测分析平台公布的AQI指数所确定的评价等级基本一致。
月份 | 隶属度 | 评价结果 | ||||||
---|---|---|---|---|---|---|---|---|
优 | 良 | 轻度污染 | 中度污染 | 重度污染 | 改进模糊层次分析法 | 熵权法 | AQI | |
2020-1 | 0.2662 | 0.0910 | 0.4665 | 0.1762 | 0 | 轻度污染 | 优 | 轻度污染 |
2020-2 | 0.4082 | 0.3942 | 0.1977 | 0 | 0 | 优 | 优 | 良 |
2020-3 | 0.3413 | 0.6105 | 0.0484 | 0 | 0 | 良 | 优 | 良 |
2020-4 | 0.3917 | 0.5693 | 0.0391 | 0 | 0 | 良 | 良 | 良 |
2020-5 | 0.4400 | 0.5601 | 0 | 0 | 0 | 良 | 良 | 良 |
2020-6 | 0.2618 | 0.6537 | 0.0846 | 0 | 0 | 良 | 良 | 轻度污染 |
2020-7 | 0.3889 | 0.5608 | 0.0504 | 0 | 0 | 良 | 良 | 轻度污染 |
2020-8 | 0.563 | 0.4372 | 0 | 0 | 0 | 优 | 优 | 良 |
2020-9 | 0.5825 | 0.4175 | 0 | 0 | 0 | 优 | 优 | 良 |
2020-10 | 0.2728 | 0.4032 | 0.3241 | 0 | 0 | 良 | 优 | 良 |
2020-11 | 0.2656 | 0.4466 | 0.2878 | 0 | 0 | 良 | 优 | 良 |
2020-12 | 0.2606 | 0.3621 | 0.3774 | 0 | 0 | 轻度污染 | 优 | 良 |
表3. 2020年1~12月评价结果
熵权法是从整体数据出发所得的权重,不具有针对性。所得结果整体偏好。AQI所确定的环境空气质量等级比模糊层次分析法和熵权法所得的环境质量等级整体偏差,是因为AQI指数评价方法只关注每月首要污染物,而忽略整体污染物的作用。若某两个月的首要污染物相同且AQI指数相同,但其他指标差异较大,AQI确定的等级为同一等级,而本文方法所确定的等级是不同的。因此,改进模糊层次分析法的模糊综合评价法更能关注环境空气质量的整体情况。
1) 2020年天津市的大气环境质量整体处于较好的水平,2020年1月和12月的评价等级为“轻度污染”,其余月份均为“优”和“良”。
2) 基于污染指数改进的模糊层次分析方法不依赖专家经验,能够更合理地确定各评价指标在环境质量评价中的相对重要性。
3) 采用模糊综合评价不仅可以给出环境质量的等级评价,还同时给出更具体的等级隶属度信息,可为环境质量决策提供更多的数据支持。
李振华,丁 春. 基于改进模糊层次分析法的环境空气质量综合评价Comprehensive Evaluation of Ambient Air Quality Based on Enhanced Fuzzy Analytic Hierarchy Process[J]. 环境保护前沿, 2021, 11(06): 1238-1245. https://doi.org/10.12677/AEP.2021.116149