第10卷第2期 2016年4月 祭琴祭祭 墨 口 供水技术 WATER IECHN0L0GY V01.10 No.2 Apr.2016 爨 研究论述 甏 丹江口水库水质的综合评价 谭浩强 , 何文杰 , 陆颖臣 , 黄廷林 , 韩宏大 , 张 南 (1.天津市自来水集团有限公司,天津300040;2.西安建筑科技大学环境与市政工程学院,西安710055) 摘要: 采用基于层次分析的模糊综合评价法对丹江口水库水质进行了等级评价,选取氨 氮、总氮、总磷、高锰酸盐指数、浊度、溶解性总固体、溶解氧、硬度、硫酸根、氯离子和pH等指标作 为评价因子,建立了丹江口水库水质综合评价指标体系。评价结果表明,在低浊期,丹江口水库水 质介于I类一Ⅱ类水质之间,在高浊期,丹江口水库水质介于Ⅱ类一Ⅲ类水质之间。 关键词: 水质评价;层次分析法;模糊综合评价法; 水源 中图分类号:TU991.2 文献标志码:A 文章编号:1673—9353(2016)o2—0001—07 doi:10.3969/j.issn.1673—9353.2016.02.001 Comprehensive evaluation of Danjiangkou Reservoir water quality Tan Haoqiang ,He Wenjie ,Lu Yingchen ,Huang Tinglin , Han Hongda 。 Zhang Nan (1.Tianjin Waterworks Group Co.,Ltd.,Tianfin 300040,China;2.School ofEnvironmental and Municipal Engineering,Xi'an University fArcohitecture and Technology,Xi’an 710055,China) Abstract:The evaluation of Danjiangkou Reservoir water quality was conducted by fuzzy comprehensive evaluation method based on analytic hierarchy process.Danjiangkou Reservoir water quality evaluation index system was established,including the principal parameters of water quality,such as ammonia nitrogen,total nitrogen,total phosphorus,permanganate index,turbidity,totl diassolved solids,dissolved oxygen,hardness,sulfate,chloride and pH.The evaluation results showed that Danjiangkou Reservoir water quality was between national standard case I and national standard case lI, in high turbidity period.it was between national standard case lI and national standard caseⅢ. Key words: water quality evaluation; analytic hierarchy process; fuzzy comprehensive evaluation method; water source 为了保障城市供水安全,在南水北调中线工程 通水之前,对丹江口水库水质进行了连续3年的跟 在此,笔者选用比较成熟的基于层次分析的模 糊综合评价法来确定丹江口水库水质等级 J。模 糊综合评价法能够客观反映水质级别的模糊性与连 踪监测分析,在此基础上需要对丹江口水库的水质 状况进行综合评价,对丹江口水库水质有直观和准 确的认识,指导南水北调中线受水区的净水厂水处 续性,将层次分析与模糊评判方法相结合,能够避免 水质评价的主观性影响,提高评价的客观性和准 确性 。 理工艺运行,便于供水企业决策分析。 基金项目:国家水体污染控制与治理科技重大专项(2012ZX07404-002) 第10卷第2期 供水技术 2016年4月 1基于层次分析的模糊综合评价模型 基于层次分析的模糊综合评价流程如下 。 级数。 V={V1, 2,V3, 4, 5} =①建立综合评价的指标体系,在此基础上,对 指标体系进行层次划分,细化各个层次,确保评价指 标划分的合理性。 {I类,Ⅱ类,Ⅲ类,Ⅳ类,V类} (1) 2.3评价指标权重的确定 (1)判断矩阵的构建 评价指标的梯阶层次结构建立之后,需要对同 一②设计评价集,对评价目标进行合理表达,对 评价等级进行科学划分。 层次的各评价因素进行两两重要性比较,标度是 ③采用层次分析法确定各评价指标权重,具 体包括判断矩阵的构建、各评价指标权重的确定和 判断矩阵的一致性检验。 ④确定各评价指标隶属度,形成隶属度矩阵。 用来衡量两个元素重要性的评价方法。在层次分析 法中采用T.L.Saaty的1—9标度区分方法,在标度 确定的基础上,构建判断矩阵(图1)。准则层的判 断矩阵A和指标层的判断矩阵A。、A 和A,,分别见 表2~表5。 all a12 … a1 ⑤按照模糊综合评价的计算步骤,采用模糊 合成算子的加权平均型,得出模糊综合评价集,采用 加权集对模糊综合评价集进行加权计算,得出最终 的总体评价结果。 A= a21 a22 … a2n : ● : ● 2水质综合评价的应用 2.1 评价指标体系的构建 an1 0 … 口M 图1判断矩阵 Fig.1 Judgement matrix 评价指标的科学性、合理性和代表性是构建综 合评价模型的第一步,必须全面、系统地反映出各个 影响因素,力求评价指标体系直观、准确、可靠和易 操作。经过对比、修改和筛选,设计出丹江口水库水 质综合评价指标体系。该评价指标体系为三层次递 阶结构,具体如表1所示。 表1 丹江口水库水质综合评价指标体系 Tab.1 Index system of water quality comprehensive A 1 表2判断矩阵A Tab.2 Judgement matrix A Ul aIi a12 aI 3 £,2 a21 a31 a22 a32 a23 a33 evaluation for Danjiangkou Reservoir 表3判断矩阵A, Tab.3 Judgement matrix A1 目标层 准则层 指标层 氨氮 指标类型 定量 A1 I1 11 1 Ul2 1 3 1 营养元素u。 总氮 l2 总磷 , 定量 定量 all 2 a12 2 a31 2 : aI1 3 a11 U21 3 a21 a31 3 631 高锰酸盐指数 水质评价 综合性 指标 浊度 : 溶解性总固体 , 溶解氧 总硬度 。 化学稳定 硫酸根 定量 定量 定量 定量 定量 定量 13 表4判断矩阵A Tab.4 Judgement matrix A2 A2 。 2 U24 £,2 3 n a21 n 2 a12 2 a22 2 a32 2 aj2 3 a23 2 a33 2 a12 4 a24 2 a34 2 性 氯离子 pH U34 定量 定量 U:4 a42 1 a42 2 a43 2 口44 2 2.2评价集的设计 根据丹江口水库水质综合评价指标体系中各评 表5判断矩阵A3 Tab.5 Judgement matrix A3 价指标的实际情况,将评价等级划分为5级(I类 水体、Ⅱ类水体、Ⅲ类水体、Ⅳ类水体、V类水体), 这也是多数模糊综合评价模型中最常采用的评价等 ·2· £ 2016年4月 谭浩强,等:丹江口水库水质的综合评价 = 续表5(Continue) 第l0卷第2期 当 >13时,CR<0.1,说明判断矩阵具有满意 的一致性,评价指标权重分配合理;否则,需要对判 A3 2 3 a21 3 a22 3 623 3 aN , o;1 3 a33 2 Ⅱ42 3 a33 3 643 3 04l 近 a33 4 044 3 断矩阵进行重新分配,直到CR<0.1,取得满意的一 致性。 (2)权重确定与判断矩阵的一致性检验 在建立判断矩阵之后,需要确定各评价指标权 重值,采用方根法进行权重向量w计算: 根据上述计算方法,最后得出权重向量。 准则层权重向量: W=[W1,W2,W3] 指标层权重向量: W1=[w w12,w13] (7) (8) (9) (10) (2) w2=[ 21,W22,W23,W24] w3=[W31,W32,W33,W34] (3) 2.4评价指标隶属度的确定 当判断矩阵取得满意的一致性时,其逻辑和权 (1)评价指标等级标准的划分 重分配表明是合理的。一致性检验的步骤如下。 通过前面的评价集可知,将评价结果进行等级 W=[W。,W2,…,W ] ①计算判断矩阵A的最大特征根A 。 = 毫 ㈩ 故采用模糊区间表示评价指标等级的分级标准见 式(11): 划分,分为5等级,评价指标等级标准也相应地分为 5等级。由于评价指标等级标准本身存在模糊性, 式中(AW) 表示向量AW的第i个元素。 ②计算一致性指标c,。 c,: (5) V={Vl,V2,V3,V4, 5}={(k0,k1),(kl,k2), (k ,k ),(k,,k ),(k ,k )} (11) 式中k1,k2,k3,k4分别表示 1与V2,V2与V3, 3 ③通过查表(表6),确定判断矩阵的平均随 表6平均随机一致性指标 Tab.6 Average random e.,onseng,u8 index 机一致性指标 。 与V ,V 与 等级之间的区间临界值;k。,k 分别表 示 等级区间的下限和 等级区间的上限。当评 价指标数值超出等级区间的上下限时,该评价指标 在距离最近等级区间的隶属度为1,其余等级区间 的隶属度均为0。 I R/l 0 J 0 I.0 58 0 90 1 J. I. 12 1 24 1 I. I. 32 1 41 I. I ④计算一致性比例C尺。 结合《地表水环境质量标准》(GB 3838— (6) CR=雨CI 2002)、《生活饮用水卫生标准》(GB 5749--2006)和 净水工艺特点要求,综上所述,各评价指标的标准等 级汇总详见表7。 当n<3时,判断矩阵永远具有完全一致性。 表7丹江口水库水质评价指标等级 指标 氨氮Ul1/(mg·L ) 类型 递增型 I类Vl 0~0.15 0 0.2 O~O.Ol O~2 0 3 Ⅱ类 2 0.15~O.5 0.2 O.5 2~4 3~lO 等级 Ⅲ类 0.5~1.O 0.5~1.0 4~6 10—15 Ⅳ类9 1.0~1.5 1.O~1.5 O.05~O.1 6~10 15~2O V类口5 >1.5 >1.5 >O.1 >10 >20 总氮Ul2/(mg·L ) 总磷U13/(mg·L ) 浊度Us2/NTU 递增型 递增型 递增型 0.O1~O.O25 0.025 O.05 高锰酸盐指数u2 /(mg·L ) 递增型 溶解性总固体U2,/(mg·L ) 递增型 溶解氧 /(mg·L ) 硬度U31/(mg·L ) 递减型 递增型 O一300 >7.5 0~l50 3oo~50o 6~7.5 150—300 50o~1 Ooo 5~6 300—450 1 Ooo~2 O0o 3~5 450~550 >2 Ooo <3 >550 ·3· 第10卷第2期 供水技术 2016年4月 续表7(Continue) 指标 类型 I类 1 硫酸根U32/(mg·L ) 氯离子U3 /(mg·L ) pH U34 Ⅱ类 2 50—150 50~150 8.O~8.2 等级 Ⅲ类 Ⅳ类 250—350 250~350 7.2~7.6 V类 5 >350 >350 <7.2 递增型 递增型 递减型 0~50 0~50 >8.2 150—250 150—250 7.6~8.0 U (2)隶属度的计算 评价指标的隶属度通过模糊分布法来计算确 Ui( 23 定。首先建立各评价指标对应的等级函数,然后根 据各评价指标的实际值对应的等级区间来推求计算 得出。由于评价指标氨氮、总氮、总磷、高锰酸盐指 数、浊度、溶解性总固体、硬度、硫酸根和氯离子均属 于递增型指标,溶解氧和pH均属于递减型指标,采 用改进的梯形隶属函数进行隶属度计算,可消除评 价等级临界值附近隶属度取值的“突变”现象。 对递增型指标,如氨氮、总氮、总磷、高锰酸盐指 0.5[1一 Ui-k3] ¨ 、: 0·5[1+ k3-//i] , < (15) 0.5[1+ Ui-k4 1J 0.5[1一 k34-Ui1-] 0 /Z <k34 < 数、浊度、溶解性总固体、硬度、硫酸根和氯离子,各 评价等级相对隶属函数计算如下。 o.s[·+ ] <kl/-L ̄s(“ )= 0.5[ 一 J ≤“ 一< l6) 。(“ )= 。.5[ 一 】 0 1≤ < 12(12) U i 10-51 1+jk}4 -式中k。,k:,k,和k 分别表示等级区间 与 :、 与 、 ,与 、V4与 的临界值; 。:, 和 分别表示 等级区间/32、 ,和v4的中间值。 对递减型指标,如溶解氧和pH,各评价等级相对 隶属函数计算如下: ≥ 1 o.5[ 一 】 ( )= < k1≤ <k12 0.5[ + ] 【-+ ] 一12≤M <Jj}2(13) o.51 1+ /x ( )= 0.5 1 k2-k2≤“z<k∞ ui 】 ≥ 23 o.5[ 一 】kl ̄>ui>k12(17) 0 ≤ l O 0 “ <k12 k】2≤M < 2 [ 一 】 o.s[ 一 卜 - 0·5[1+ 】 (14) /x .(u )= o.51 1.k2-…ui l k2≤“ <k23 : 8 , ( )= k23≤ <k3 o.5【l+ ui- k3 1J o.s[ + k3≤ <.]}34 0.5[ 一 k3-ui】 ≥ 34 叫1一 ] ·4· 2016年4月 谭浩强,等:丹江口水库水质的综合评价 第1O卷第2期 u > 12 k2一 。 k2-k12 一+ 一k2 + 一 k23一k2 k2≥ >k23 (/Ai)= k3一 k3一k23 一(19) k23≥ >k3 矗mk3≥ >k /.ti ̄/ ̄34 k3 k34一k3 Lu 0 >k23 0.5[1一 ks-/ ̄i ( )= ,> s R0·5【t+七Ui -一k3 】 ,≥ > (2o) 0.5[1+ k4-ui】 1一 】 O u >k >Ij} 旧 1鞋](28)  ̄.Lv5(U )= o.5【1一 0.5[L 1+uU i - k ̄ J> (21) 式中k。,k2,k,,k ,k毪,后∞和k的意义同前所述。 2.5模糊综合评价 在本文2.4部分基础上,计算得出营养元素 、综合性指标 和化学稳定性 的隶属度矩 阵,分别为 。、 :和 ,。 』:If_ R2= l水质等级I I 类lⅡ类IⅢ类lⅣ类l V 类l (23) 进行水质评价,具体水质情况见表9。 表9水质概况 Tab.9 Ovm'view ofwater quality 项目 R3= 低浊期 高浊期 (24) 氨氮/(mg·L ) 总氮/(mg·L ) 总磷/(mg·L ) 高锰酸盐指数/(mg·L ) 0.06 1.55 0.02 1.6O O.10 1.65 0.o4 2.2O ·5· 第10卷第2期 续表9(Continue) 供水技术 2016年4月 重分配合理。 高浊期 17 项目 浊度/NTU 低浊期 1.9 因此,准则层权重向量和指标层权重向量如下: W=[0.163 4,0.297 0,0.539 6] W1=[0.163 4,0.297 0,0.539 6] (32) (33) 溶解性总固体/(mg·L ) 溶解氧/(mg·L ) 总硬度/(mg·L ) 硫酸根/(mg·L ) 氯离子/(mg·L ) pH 160 9.5 140 33 5 8.10 135 6.5 120 25 3.5 7.95 v ̄2=[0.164 4,0.284 8,0.105 1,0.445 7] (34) w3=[0.285 7,0.142 9,0.285 7,0.285 7] (35) 表12指标层权重向量 计算 Tab.12 Calculation of index layer weight vector w2 A2 1 按照本文2.3部分的步骤,计算准则层和指标 层评价指标的权重,具体详见表l0~表13。 表10准则层权重向量w计算 Tab.10 Calculation of criteria layer weight vector W 1 2 1 2 1/2 1 。 U24 2 3 1 3 l/3 1/2 1/3 1 w2 0.164 4 0.284 8 0.105 1 O.445 7 A… A l l W A 丑x CR 2 1 2 3 1/2 l 2 1/3 1/2 1 0.163 4 0.297 0 0.539 6 3.009 2 0.007 9 4.070 9 0.026 2 3 l/2 1/3 U24 3 2 表13指标层权重向量w 计算 表11指标层权重向量 计算 Tab.13 Calculation of index layer weight,#eetor w3 Tab.11 Calculation of index layer weight vector wl Al 1l A3 £,3 。 2 2 1 2 2 3 l 4 1 w3 0.285 7 0.142 9 A… C尺 UI1 1 12 1/2 13 1/3 Wl 0.163 4 A Ⅱ 32 1 1/2 1 1 1/2 1/2 1 1 1 1 l2 2 3 1 2 1/2 1 0.297 0 0.539 6 3.009 2 0.007 9 4 ” 0 0.285 7 0.285 7 Ul3 4 权重向量W、W。、w2和 的判断矩阵CR均 按照本文2.4部分的步骤,计算各定量评价指 标的隶属度,结果具体详见表l4。 <0.1时,判断矩阵具有满意的一致性,评价指标权 表14隶属度计算w 计算 水质期 指标 氨氮mg/L 总氮/(mg·L ) 总磷/(mg·L。) 0.O6 1.55 0.02 I类 0.643 0 0 Ⅱ类 0.357 O 0.833 Ⅲ类 O 0 O.167 Ⅳ类 0 0.4l7 0 V类 0 0.583 0 高锰酸盐指数/(mg·L ) 浊度/NTU 1.60 1.9 0.643 0.620 0.357 0.380 0 0 0 0 0 0 低浊期 溶解性总固体/(mg·L ) 溶解氧/(mg·L ) 总硬度/(mg·L ) 硫酸根/(mg·L ) 氯离子/(mg·L ) pH l60 9.5 140 33 5 8.10 0.792 0.864 0.559 0.627 0.737 0 0.208 0.136 0.441 0.373 0.263 1.0oo 0 0 0 0 0 0 O 0 0 0 0 0 0 0 0 0 0 O 氨氮/(mg·L ) 总氮/(mg-L ) 0.10 1.65 0.591 0 0.409 0 0 0 0 0.313 0 0.687 高浊期 总磷/(mg·L ) 高锰酸盐指数/(mg·L ) 浊度/NTU 0.04 2.20 17 O 0.40o 0 0 0.600 0 0.900 0 0.10O 0.100 0 0.900 0 0 0 溶解性总固体/(mg·L ) 135 0.811 0.189 0 0 0 ·6· 2016年4月 谭浩强,等:丹江口水库水质的综合评价 续表14(Continue) 第lO卷第2期 水质期 指标 溶解氧/(mg·L ) 总硬度/(mg·L ) I类 6.5 120 Ⅱ类 0.833 0.357 Ⅲ类 0.167 O Ⅳ类 0 0 V类 0 0 0 0.643 高浊期 硫酸根/(mg·L ) 氯离子/(mg·L ) pH 25 3.5 7.95 0.667 0.741 O 0.333 0.259 0.375 0 0 0.625 0 0 0 0 0 0 根据本文2.5部分的计算步骤,得出指标层营 养元素 、综合性指标 、化学稳定性 和丹江 口水库水质评价的模糊综合评价集,具体详见表l5 和表16。 表15指标层的评价结果 Tab.15 Evaluation result ofindex layer 项目 低浊期 高浊期 营养元素 [0.105,0.508,0.090,0.124,0.173] [0.097,0.067,0.486,0.147,0.204] 综合性指标 [0.751,0.249,0,0,0] [0.151,0.490,0.103,0.256,0] 化学稳定性 [0.460,0.540,0,0,0] [0.491,0.331,0.178,0,0] 表16丹江口水库水质的评价结果 Tab.16 Evaluation result of water quality for Danjiangkou Reservoir 【 水质期 评价结果 水质期 评价结果 I 低浊期 [0.488,0.449,0.015,0.020,0.028] 高浊期 [0.326,0.335,0.206,0.100,0.033] 根据式(31)计算丹江口水库水质评价分类等 期,原水需要进入预沉池,对浊度进行去除,降低后 级,如表17所示。可以看出,低浊期的水体介于I 续净水工艺处理负荷,适当增加混凝剂投加量,调整 类一Ⅱ类水质之间,高浊期的水质介于Ⅱ类一Ⅲ类 滤池运行工况。 水质之间。 表17丹江口水库水质的评价分类等级 参考文献: Tab.17 Classification grade ofwater quality evaluation [1]谭浩强,何文杰,韩宏大,等.丹江口水库水质监测分 for Danjiangkou Reservoir 析[J].供水技术,2015,9(6):1—5. 水质期 等级 水质期 等级 [2] 魏炜,岳丹丹,潘俊.基于层次分析的模糊综合评价法 低浊期 1.65 高浊期 2.18 在水质评价中的应用[J].供水技术,2015,9(5):8一 3结论与建议 l2. 应用基于层析分析法的模糊综合评价模型对丹 [3] 黄源.T市供水管网浊度特性的研究及其风险评估的 江口水库水质进行了等级评价,低浊期的水体介于 应用[D].哈尔滨:哈尔滨工业大学,2013. I类~Ⅱ类水质之间,高浊期的水体介于Ⅱ类~Ⅲ 类水质之间。在进行多水源联合调度时,需要考虑 丹江口水库不同时期的水质情况,从而形成合理的 调水方案。在净水厂实际生产运行中,在水质低浊 期,可以采用超越管的方式,直接超越预沉池,进入 净水工艺前端,适当减少混凝剂投加量;在水质高浊 工 作E收-者m稿a简日il介:期tan:2谭男饮究h0ao1浩用工,q6i水作强a—博n0安g。( 2士1全—96,8@1保45 1一主障26要技.c)从术o,m 事研 ·7·