灰色聚类决策在水环境评价中的应用

《灰色聚类决策在水环境评价中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、灰色聚类决策在水环境评价中的应用刘宜峰1　张治1　谢志毅2（1.宁波市鄞州区水利水电勘测设计院　浙江　宁波　315040　2.宁波市鄞州区大嵩水利管理站　浙江　宁波　315141）摘要：为了综合反映水环境实际质量，采用灰色聚类决策方法进行综合评价，并借助计算机解决繁琐的计算，使水质评价既简单又科学。关键词：水环境；灰色聚类决策；应用目前在水环境质量评价中应用最广的当属单因子极值法，即以我国现行的地面水环境质量标准为评价依据的方法。根据《地面水环境质量标准》要求，参加评价的指标值单项超标，即表明该项使用功能不能保证。该方法目标明确，评级直观，操作简单，宜于

2、接受，适用性很广，但其“一票否决”式的评价方式具有一定的片面性，不能综合反映水环境的实际情况。例如在宁波市鄞州区平原河网水质评价中，由于极个别指标超标，整体水环境被判为劣V类，致使人们不再放心使用身边的河网水，而一律改用自来水。这给严重缺水的宁波市带来更大的用水压力。灰色聚类决策问题是根据若干指标将一些研究对象进行分类的问题。水质评价是根据各种水质指标将实测水环境质量归为某一级别的问题，这一问题可以用灰色聚类决策方法很好地解决。这种方法以综合评价的方式充分考虑了所有指标的权重，即考虑各个指标对水环境的影响，能更切实际地反映问题。运用VB语言编程计算，可以

3、解决灰色聚类决策方法繁琐的计算过程，提高工作效率。本文着重介绍灰色聚类决策方法在水环境质量评价中应用的计算过程。1灰色聚类决策综合评价方法灰色聚类是以灰数的白化函数生成为基础，它将收集到的分散信息，按照灰类进行归纳，以判断聚类对象所属灰类。记i=1，2，…，n为聚类对象；j=1，2，…，m为聚类指标；k=1，2，…，K为聚类灰数，即灰类。dij为第i个聚类对象对于第j个聚类指标的样本值，D是以dij为元素的样本矩阵。fjk为第j个聚类指标的k灰类的白化权函数，fjk∈[0,1]，如图1-1所示。图1-1　白化权函数fjk图中λjk为fjk的阀值，有客观阀

4、值和相对阀值两种。作者简介：刘宜峰（1976-），男，助理工程师，大学本科，主要从事水利工程设计工作。张　治（1979-），男，助理工程师，大学本科，主要从事水利工程设计工作。谢志毅（1964-），男，助理工程师，大学专科，主要从事水利工程管理工作。7记ηjk为灰色聚类权，它表示第j种指标属于第k灰类的权重。关于ηjk的计算可按下述3种情况分别考虑：(1)在聚类指标的意义、量纲相同时(1－1)若聚类指标的意义、量纲不同，且不同指标的样本值在数量上相关很大时，不能按上式计算ηjk，否则将会引起评估偏差。(2)当聚类指标的单位不同及指标绝对值相差很大时，先进

5、行无量纲处理：(1－2)式中，Sjk为第j种指标的第k个灰类的灰数(标准值)，或取白化权函数图(图1-1)中的阀值λjk(即fjk(d)等于1的界值)。Sj为第j种指标的参照标准。灰色聚类权值ηjk为：(1－3)(3)聚类权ηjk是衡量多指标对同一灰类的权重。为了避免指标重要而指标值相对较小造成评估偏差的弊病，可采用层次分析法(AHP)、德尔菲法等，根据聚类指标的相对重要性而确定灰色聚类权值ηjk。设σik为灰色聚类系数，它反映第i个聚类对象隶属于第k灰类的程度，按下式计算：(1－4)式中fjk(djk)系由样本值djk查相应的白化权函数图或代入相应的白

6、化权函数计算式而求得的白化权函数值。灰色聚类决策矩阵σc:聚类行向量σi=(σi1，σi2，…，σik)若有σik＊满足(1－5)称聚类对象i属于灰类k*。也就是说，在聚类行向量σi中找出最大聚类系数，该最大聚类系数所属灰类即为聚类对象i所属灰类。2应用实例在水质监测值中，将采样点作为聚类对象，监测指标作为聚类指标，而水质标准值作为各灰类白化权函数值的阀值。如此确定各评价指标的各灰类的白化权函数，进行聚类决策计算，得出聚类对象所属灰类，即判断出各采样点的水质级别。7下面以宁波市鄞州区平原河网水环境质量评价为例，具体说明灰色聚类决策综合评价的方法与过程。在

7、七个采样点（辰交、石碶、望春、横溪水库出口、铜盆浦、潘火、五乡等七处分别编号为1～7号）对8项污染指标进行了监测，实测值列于表2－1。现将七个采样点作为聚类对象(i=1，2，…，7)，8种污染指标作为聚类指标(j=1，2，…，8)，dij为第i个采样点对于第j个污染指标的实测值。表2－1　　　　　　　鄞州区平原河网水质实测值（mg/L）测点i污染指标j1234567溶解氧4.846.324.259.767.098.226.48高锰酸盐指数4.315.317.222.383.297.632.11生化需氧量2.513.914.491.281.646.221.

8、84挥发酚0.0020.0020.0020.0020.0020.0020.002