郭宁

摘要 以山东省滨州市重点流域水质为例，采用RAGA-PPCE模型对水质进行综合评价。选取对水体影响较大的pH、化学需氧量（COD）、总磷（TP）、总氮（TN）、氨氮（NH3-N）、溶解氧（DO）6项水质评价因子，构造水质综合评价分级标准。结果表明：RAGA-PPCE 模型水质评价结果与模糊综合法评价结果基本相同；RAGA-PPCE模型评价结果客观、合理，能够有效应用于水质综合评价。

关键词 RAGA-PPCE模型；水质评价；投影寻踪

中图分类号 S181；X52 文献标识码 A 文章编号 0517-6611（2017）07-0053-03

Comprehensive Evaluation of Water Quality in Binzhou Basin Based on RAGA-PPCE Model

GUO Ning

（School of Civil Eengineering，Binzhou University，Binzhou，Shandong 256600）

Abstract Taking key river water quality in Binzhou City of Shandong Province as an example，comprehensive evaluation of water quality with RAGA-PPCE model was conducted.Selecting 6 water quality evaluation factors including pH，dissolved oxygen，total phosphorus，total nitrogen，ammonia nitrogen，dissolved oxygen，comprehensive evaluation of water quality classification standard was constructed.The results showed that RAGA-PPCE model water quality evaluation results and fuzzy comprehensive evaluation results are basically the same.The evaluation result of RAGA-PPCE model was objective and reasonable，and it can be applied to comprehensive evaluation of water quality.

Key words RAGA-PPCE model；Water quality evaluation；Projection pursuit

在城市生态水系建设中，实施严格的水资源管理及水体污染监测都离不開评价水质状况。现今的评价方法有神经网络方法[1]、模糊综合层次分析方法[2]、单因素分析法[3]、灰色关联分析法[4]、主成分分析法[5]、综合水质标识指数法[6]、综合污染指数模型[7]、Topsis算法[8]，这些算法都已在实际水质评价中运用并得到很好的效果。投影寻踪法（PP）[9]属于直接由样本数据驱动的探索性数据分析方法，它把高维数据通过某种组合投影到低维子空间上进行聚类，通过分析低维空间数据点的变化特征来研究高维数据结构特征。聚类就是按照一定的规则或者标准将对象聚为一个个集合，所得到的分布是聚类前未确知的。它以类内具有较大密集度（高内聚）而类间具有相对大的散开度（低耦合）为目标，寻找一维投影方向，并根据每个样本相应的综合投影值对样本进行综合评价分析。

基于实数编码的加速遗传算法求解模型（RAGA-PPCE）将标准遗传算法改为实数编码的形式，采用加速循环的方式使求解速度加快[10]。然而，目前采用该方法对整个流域水系水质的评价却鲜见报道。笔者以滨州市主要河流水质（杏花河、徒骇河、孝妇河、漳卫新河4个断面）为对象，采用RAGA-PPCE找出最佳投影方向，构造水质综合评价分级标准，以期为我国河流水质评价提供科学依据。

1 RAGA-PPCE模型的建立

1.1 RAGA-PPCE算法

首先构造指标函数[11-12]，即把N维数据变成a=（a1，a2，…，an），为投影方向的一维投影值zi：

zi=nj=1aj×xi，j

然后根据{zi/i=1，2，…，n} 的一维散布图进行分类，其中，a为单位长度向量。投影函数可构造为：

Q（a）=σz×ρz

式中，σz为投影值zi的标准差，ρz为投影值zi的密度[13-14]。

1.2 加速遗传算法求解投影寻踪模型的步骤

RAGA的求解步骤[15]：①采用二进制编码。②初始化群体。③适度评价，目标值是0的情况，定义适应度函数。④产生第1子代群体。⑤杂交产生第2代子群。⑥变异产生第3代子群。⑦继续迭代，重复进行演化。⑧加速循环。选出优秀个体对应变化区间重新迭代，如此循环。

2 应用实例

2.1 数据来源

滨州市主要水系有黄河水系、小清河水系和海河水系。水系内人口稠密，经济发达，流域内河流众多，水系发育，小清河水系有孝妇河、杏花河、小清河，海河水系有徒骇河、马颊河、漳卫星河等。滨州市多年平均水资源总量为30.67亿m3，过境河流水资源量11.95亿m3。流域面积5 788 km2。人均水资源占有量267 m3，属于人均水资源量小于500 m3的水危机地区。笔者选取小清河水系和海河水系为研究对象，选取影响较大的氨氮（NH3-N）、总磷（TP）、总氮（TN）、化學需氧量（COD）、溶解氧（DO）、pH作为水质综合评价因子。

2.2 样本生成及参数设置

2.2.1 样本生成。选取孝妇河断面袁家橋（长山）（1#）、马颊河断面李辛桥（大店）（2#）、漳卫新河断面小泊头桥（3#）、徒骇河断面申桥（4#）作为采样断面，每条河流断面取173组水质数据集，用于求取 RAGA-PPCE模型最佳投影方向 a。

2.2.2 参数设置。RAGA-PPCE算法参数设置为最大迭代次数T=200，种群规模N=50[16-17]。

2.3 构造水质评价分级标准及结果

根据地表水环境质量国家标准（GB 3838—2002）及地表水水域环境功能和保护目标，按功能高低依次划分为Ⅰ～Ⅴ类水，所选指标对应的各类水环境质量标准见表1。

由表2可知，孝妇河袁家桥（长山）断面 Ⅴ 类水质以上占60.08%，马颊河李辛桥（大店）断面 Ⅴ 类水质以上占68.74%，漳卫新河小泊头桥断面 Ⅴ 类水质以上占44.49%，徒骇河申桥断面 Ⅴ 类水质占95.08%，前3个断面的水质都没有 Ⅰ 和 Ⅱ 类。徒骇河申桥断面水质最好，漳卫新河小泊头桥断面水质最差。a=（0.045 0，0.215 4，0.299 6，0.884 0，0.124 1，0.255 0）为最佳投影方向。

图1～4是河流水质评价指标分布情况，图中5条直线是5类划分标准（从上而下分别是 Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ），可以看到水质的分类并非固定不变，而是呈动态变化的，数据跨过夏、冬两季，随着冬季到来河水流量减小，水质变差，从图中可以明显看出，随时间的推移水质越来越差。4个断面的水质都是冬季最差。

2.4 评价结果

用MATLAB软件计算出投影值，见表3（由于篇幅关系每个断面从173个中任选出10个指标）。由表1可知，RAGA-PPCE模型水质评价结果与模糊综合法结果相同。从最佳投影方向来看，各评价指标均为正值，且相差不大，说明各指标投影方向一致，对于水质评价等级的影响相同。通过定量计算，得到以下评价：4条河在 173个数据水质均有波动，但1#、2#、3#均未达到过 Ⅰ、Ⅱ 类水标准，最高等级是 Ⅲ 类水，且各水质之间差异较为明显。4#水质在173个数据内的两极分化比较明显，水质波动范围较大，总体水质最好；漳卫新河小泊头桥断面水质有54%处于 Ⅴ 类水以下，3#达到劣 Ⅴ 类水质次数最多，相比于其他河水其水质条件最差。

3 小结

该结果表明，RAGA-PPCE模型评价结果客观、合理，能够有效应用于水质综合评价。利用 RAGA-PPCE模型计算最佳投影方向a，不但提高了该模型的评价精度，而且为模型最佳投影方向 a 的选取提供了一种新的途径和方法。

参考文献

[1] 赵小勇，付强，邢贞相，等.投影寻踪模型的改进及其在生态农业建设综合评价中的应用[J].农业工程学报，2006，22（5）：222-225.

[2] 谭永明，孙秀玲.基于加速遗传算法与投影寻踪的水质评价模型[J].水电能源科学，2008，26（6）：42-44.

[3] 王洪梅，卢文喜，辛光，等.灰色聚类法在地表水水质评价中的应用[J].节水灌溉，2007（5）：20-22.

[4] 王艳琼，白秀琴.基于BP神经网络模型的水质评价及预测[J].武汉工业学院学报，2007，26（1）：64-67.

[5] 陈兴，程吉林，刘芳.BP神经网络用于水质评价的参数确定[J].水利与建筑工程学报，2007，5（1）：12-15.

[6] 曾永，樊引琴，王丽伟，等.水质模糊综合评价法与单因子指数评价法比较[J].人民黄河，2007，29（2）：45.

[7] 沃飞，陈效民，吴华山，等.灰色聚类法对太湖地区农村地下水水质的评价[J].安全与环境学报，2006，6（4）：38-41.

[8] 徐健，吴玮，黄天寅，等.改进的模糊综合评价法在同里古镇水质评价中的应用[J].河海大学学报（自然科学版），2014，42（2）：143-149.

[9] 徐兵兵，张妙仙，王肖肖.改进的模糊层次分析法在南苕溪临安段水质评价中的应用[J].环境科学学报，2011，31（9）：2066-2072.

[10] 申剑，史淑娟，周扬，等.基于改进灰色关联分析法的丹江口流域地表水环境质量评价[J].中国环境监测，2014，30（5）：41-46.

[11] 盧文喜，祝廷成.应用人工神经网络评价长春南湖水的营养状态[J].地理科学，1999，19（5）：462-465.

[12] 张欣莉，丁晶，李祚泳，等.投影寻踪新算法在水质评价模型中的应用[J].中国环境科学，2000，20（2）：187-189.

[13] 付新峰，谷晓伟，刘晓岩，等.乌梁素海生态功能定位初步分析[J].人民黄河，2008，30（10）：61-62.

[14] 管延海，李强，柴成繁.模糊数学方法在天津市地下水水质评价中的应用[J].地下水，2008，30（2）：27-28.

[15] 周惠成，董四辉.基于投影寻踪的水质评价模型[J].水文，2005，25（4）：14-17.

[16] 汪明武，金菊良，李丽.基于实码加速遗传算法的投影寻踪方法在砂土液化势评价中的应用[J].岩石力学与工程学报，2004，23（4）：643.

[17] 杨晓华，杨志峰，郦建强.水质综合评价的遗传投影寻踪插值模型[J].环境工程，2004，22（3）：69-71.