改进的灰色聚类法在地下水质量评价中的应用

2012-11-27李亚松张兆吉费宇红孟素花郭春艳

水资源保护 2012年5期

李亚松，张兆吉，费宇红，钱永，孟素花，郭春艳

(中国地质科学院水文地质环境地质研究所，河北石家庄 050061)

地下水不仅是供给人们日常用水的主要来源，同时也是复杂的生态环境系统中敏感的组成因子之一，地下水的变化往往会影响生态环境系统的天然平衡状态[1]。鉴于地下水在生物圈和地质环境中的特殊地位，对地下水质量进行全面合理的评价显得尤为重要。

国内学者前期进行的区域地下水质量评价，其评价指标体系多为常规指标，如溶解性总固体、总硬度、耗氧量和三氮等[2-5]，评价结果只代表地下水的一般质量。随着地下水有机污染工作的开展，越来越多的有机指标可以被检测出来，目前，全世界已在水中测定出2221种有机化学污染物，其中765种存在于自来水中，20种已被确认为致癌物，23种为可疑致癌物，18种为促致癌物，56种为诱变物质[6]。经调查发现，不同的有机化合物还对人体各器官及系统有不同的危害，如:三氯乙烯刺激皮肤，对肺、肝、肾等内脏器官有损害;苯刺激皮肤、眼睛，对神经系统、免疫功能、胃肠有损害[7]。很多有机污染物在地下水环境中很难通过自然降解过程去除，很可能会长期存在并发生累积，所以将有机毒理指标引入评价指标体系就显得尤为重要。已经有学者就内梅罗指数法和模糊数学法等在综合指标评价中的问题进行了研究[8-10]，但是对于同样应用较为广泛的灰色聚类法，则鲜有报道。

灰色聚类法是目前地下水质量评价工作中较为常用的方法，但当评价指标体系中痕量指标较多时，其评价结果容易失真。笔者通过修正其聚类权的计算方法将原有方法加以改进，得到改进后的灰色聚类评价法评价结果同模糊数学综合评价法相乘相加算子的评价结果有较好地一致性，可以较好地反映区域地下水质量状况。

1 研究区概况

研究区位于海河流域子牙河水系滹沱河冲洪积扇，主要包括石家庄市以及下属6区17个市县。研究区面积约6000 km2，地势平坦，地面标高由山前的100 m降至扇缘的50 m左右，地面坡降0.25% ～0.1%。研究区属暖温带半湿润半干旱大陆型季风气候区，1956—2002年多年平均气温13.3℃，平均降水量为531.4 mm。降水量具有年内和年际变化不均匀的特点，6—8月份降水量占全年降水量的70%～80%。包气带多为层状非均质黄色、黄褐色亚砂土、粉细砂、中粗砂及亚黏土组成，其下伏含水层为上更新统－全新统砂卵石层[11]。研究区天然条件下地下水主要接受大气降水入渗补给，主要排泄方式为人工开采。

研究区是全国著名的工业中心及产粮基地。自20世纪80年代以来，工矿企业不断增多，并且其大都建立在透水性好、水量充足的河流两岸，河床多为砂性土，大量工业及生活废水通过无防渗的沟渠排入河流，使得污染物连续渗漏，以直接或间接方式进入地下含水层，造成地下水，特别是浅层地下水的污染。此外，随着农业生产的发展、人口的增长和人民生活水平的提高，对水资源的需求量越来越大，地下水的持续过量开采已经在局部地区形成了水位降落漏斗，地下水位的不断下降引发了一系列的环境地质问题，其中代表性的是水质变差、地面沉降等等[12]。

2 案例研究

2.1 评价数据来源

研究以“滹沱河冲洪积扇地下水污染调查评价”项目为依托，按照《地下水污染调查评价规范》中相关要求[13]，在太行山前滹沱河冲洪积扇平原区共采集地下水样品508个，样品测试单位为国土资源部水文地质环境地质专业测试中心，平行样品测试单位为上海澳实分析检测公司，所有样品的采集和测试均进行了严格的质量控制。

2.2 参评指标体系及评价标准

综合对比实验室测试项目及《地下水污染调查评价规范》，确定参评指标体系为50项(表1)，评价标准为1993年制订、2007年修订的GB/T14848－2007《地下水质量标准》[14]。

表1 评价指标体系

2.3 评价方法原理

灰色聚类评价法的基本思路就是以灰数的白化函数生成为基础，将收集的聚类对象观测值的分散信息，按照灰类进行归纳，计算水体水质中各因子的实测浓度与各级水质标准的灰度(关联度)，然后根据灰度大小确定水体水质的级别。对处于同类水质的不同水体可通过其与该类标准水体的灰度大小进行优劣比较。

2.3.1 评价流程

a.记i=1，2，…，n为聚类样本，即各个测点或断面;记j=1，2，…，m为聚类指标，即各评价指标;记k=1，2，…，p为灰类别，即按某一标准的分级。

将滹沱河冲洪积扇508个采样点作为聚类样本(即 i=1，2，…，508)，把50项评价指标作为聚类指标(即 j=1，2，…，50)，GB/T14848—2007 按《地下水质量标准》可以分为5个灰类(即k=1，2，…，5)。

b.给出聚类白化数Xij。聚类白化数Xij为各样品的测试结果，本文将各个采样点上各项检测指标的实测值作聚类白化数Xij。

c.根据选取的分级标准，构造出白化函数。

对Ⅰ级水的白化函数，即k=1:

对Ⅱ级至Ⅳ级水的隶属函数，即k=2，3，4:

对Ⅴ级水的隶属函数，即k=5:

式中:Yijk为i水样中第j个污染因子对k级标准的白化函数;Cjk为j污染因子k级标准限值。

当Xij给定，可以用以上白化函数求出i水样中第j个污染因子分别对各级水的灰类，然后按聚类指标所属类别，确定不同的白化函数。

d.聚类权的确定。考虑量纲的影响，应对原始数据进行无量纲化处理，可以使用平均标准方法或背景值方法，即各因子的污染分级标准值分别除以其对应的平均标准值或Ⅱ类水标准值。以前一种处理方法为例，设处理后的第j个污染因子对应k级标准的处理结果为Rjk，则聚类权的计算公式为

式中:Njk为第j个污染因子第k个灰类的权重。

e.聚类系数的计算。聚类系数是通过灰数白化函数的生成而得到的，它反映了聚类样本对灰类的亲疏程度，其计算公式为

式中:Fjk为第i个样本关于第k个灰类的聚类系数。

将所得出的聚类系数进行归一化处理，按最大隶属原则，最大值所对应的灰类即为此样本所属类别[15]。

2.3.2 评价结果分析

评价后发现，评价结果全部为Ⅰ类水，显然不符合客观实际。主要是由于评价指标中大部分为痕量指标，用传统方法计算聚类权值，造成了聚类系数在Ⅰ类水处聚集，导致评价结果失真。

2.3.3 改进聚类权

鉴于上述情况，需使得聚类权在相同指标不同灰类及相同灰类不同指标间都可以有梯度关系，故对聚类权的计算方法修正如下。

a.计算相同级别对于不同指标的权重。首先将各指标的单位统一化，假定指标i的j级标准值为Kij，相同级别不同指标的权重为Mij:

b.计算相同指标不同级别的权重。由于有机指标在Ⅱ类和Ⅲ类标准值间跨度很大，故将参评指标分为有机指标和无机指标两类，分别按下式进行运算得出各指标不同级别的均权值Nij。

式中:n为无机指标(有机指标)个数。

c.将两权重复合运算，归一化后得到聚类权值Pij:

2.3.4 修正结果分析

利用修正后的聚类权计算聚类系数后得出，评价结果失真的现象得以消除，计算各水质级别样品数占样品总数的比例可知，地下水质量主要为Ⅱ类水和Ⅲ类水，超Ⅲ类水占样品总数的29.6%(表2)。由评价结果分布图(图1)可以看出，超Ⅲ类水点主要分布在城区、城镇周围区、排污河道两侧及引污水灌溉农田地区，符合野外污染源荷载调查对该区地下水质量状况的预判，因此利用改进的灰色聚类法可以对区域地下水质量做出较为客观的评价。

对主要超标因子进行分析可知，影响区域地下水质量的指标主要为溶解性总固体、总硬度、硝酸盐氮等无机指标，可以得出，工农业发展及人类活动加剧是造成本区水质变化的主要原因，可以得出，工农业发展及人类活动加剧是造成本区水质变化的主要原因，为保证本区地下水资源的可持续利用，应在一定程度内限制工业用水量、规范工业排污量、提高农业工作中农药、化肥的使用效率等。

表2 评价结果对比 %

3 结论

由评价结果分布图(图1)可以得出，区内绝大多数地区水质状况良好，影响区域地下水质量的指标主要为溶解性总固体、总硬度、硝酸盐氮等无机指标，水质较差区主要分布在城区、城镇周围区、排污河道两侧、地表污染水体分布区及引污水灌溉农田地区。可以得出，工农业发展及人类活动加剧是造成本区水质变化的主要原因，为保证本区地下水资源的可持续利用，应在一定程度内限制工业用水量、规范工业排污量、提高农业工作中农药、化肥的使用效率等。