基于模糊综合评价法的大伙房水库上游水质评价及预测

来源：https://www.shuizhifenxi.com/ 作者：余氯检测仪 时间：2019-08-01

　　摘要：运用模糊综合评价法及层次分析法对辽宁省浑河大伙房水库上游北杂木、古楼及台沟断面丰、枯水期的水质进行评价，同时将水质评价结果(水质模糊综合评价指数)与对应上游来水量进行相关性分析，建立来水量与水质模糊综合评价指数的非线性回归模型，并利用该模型进行水质预测。结果表明：非线性回归模型具有较高的拟合效果，应用该模型可以对上游各来水断面的水质类别进行预测。

　　关键词：水质评价;水质预测;模糊综合评价法;层次分析法;非线性回归模型;大伙房水库上游;水质模糊综合评价指数

　　大伙房水库输水工程是解决辽宁省中部沈阳、鞍山、抚顺、营口、辽阳、盘锦及南部大连七城市水资源短缺问题的主要供水水源地。大伙房水库水质级别直接影响辽宁各大城市居民饮用水安全，所以对近年来大伙房水库入库水质进行评价及预测水质显得十分必要。

　　水体本身是一个多元的复杂体系，影响水质评价的物理、化学及生物因素具有不确定性和模糊性，常用的评价方法可概括为单因子指数评价法和综合评价法两大类，其中综合评价法又包括模糊综合评价法、模糊模式识别法、灰色系统综合评价法、灰关联分析法、人工神经网络法、主成分分析法、综合水质标识指数法及标准类别指数评价法[1-9]。水质预测模型大致可概括为数理统计预测方法、灰色系统预测法、神经网络模型预测法、水质模拟模型预测法及混沌理论预测法等[10-14]。

　　本文运用模糊综合评价法及层次分析法对大伙房水库上游水质进行综合评价，建立来水量与水质模糊综合评价指数的非线性回归模型，可用于突发洪水影响下的水质预测，进行水质污染预警。

　　1模糊综合评价法

　　模糊综合评价法的基本思路是：首先运用隶属函数对各单项污染因子分别进行评价，其集合构成一个模糊矩阵R;再考虑各单项污染因子对水质影响的程度大小给予不同的权重W总，进而构成一个权重集;最后通过复合运算，得出水质模糊综合评价指数FCL，藉此判别水质级别。

　　1.1单因子隶属度函数

　　表征水体质量的各项因子的监测值可构成一个因子集{xi}，而根据不同用途和特征划分的水质级别标准又构成一个标准集{Sij}。

　　对于第1级水质，隶属度函数为

　　yi1=1xi≤Si1

　　xi-Si2Si1-Si2Si1

　　0xi≥Si2(1)

　　式中：i为第i种因子，(i=1，2，…，n);xi为第i种因子的实测浓度;Si1为第1级水质标准中规定的第i种因子标准浓度值;Si2为第2级水质标准中规定的第i种因子的标准浓度值。

　　对于第2级至(m-1)级水质，其隶属度函数为

　　yij=1xi≤Sij 　　xi-Si，(j-1)Sij-Si，(j-1)Si，(j-1)

　　xi-Si，(j+1)Sij-Si，(j+1)Sij

　　0xi≥Si，(j+1)或xi≤Si，(j-1)

　　j=2，3，…，(m-1)(2)

　　对于第末级m水质，隶属度函数为

　　yim=1x≥Sim

　　xi-Si(m-1)Sim-Si(m-1)Si(m-1)

　　0xi≤Si(m-1)(3)

　　式中：xi为第i种因子实测值;Sim为第i种因子的第m级评价标准。

　　1.2模糊矩阵的建立

　　利用以上隶属函数，按照水质分级标准，可以计算出各项因子分别对各级水质指标的隶属度，从而构成模糊矩阵R如下：

　　R=(yij)n×m=y11…y1m

　　yn1…ynm(4)

　　1.3确定权重

　　权重是基于污染物对水质影响的大小来确定的，准确地计算各污染指标权重是评价水质的重点。运用层次分析法确定权重可有效避免主观臆断。

　　(1)建立层次结构模型。将模型构成划分为目标层A、准则层B和方案层C三个层次。

　　(2)构造判断矩阵。通过对准则层相互比较，确定各准则对于目标的标度值，即构造目标层判断矩阵A-B;对方案层相互比较，确定各方案对于准则层的标度值，即构造准则层判断矩阵B-C。在层次分析法中，引进矩阵判断标度(1至9标度法)。

　　(3)层次单排序。根据矩阵理论计算得出同一层次各因素对于上一层次对应因素的相对重要性权重值W，它是本层次中所有元素对上一层次而言进行重要性排序的基础。层次单排序可归结为求判断矩阵的特征值和特征向量问题。即对于判断矩阵A，计算满足：AW=λmaxW的特征根与特征向量。

　　求出矩阵的最大特征值后，还必须进行矩阵一致性检验。评价判断矩阵的一致性检验指标为CI=λmax-nn-1(n为判断矩阵的阶数)，将CI与平均一致性指标RI进行比较，即CR=CIRI，当CR<0.1时，认为判断矩阵满足一致性，否则就要对判断矩阵进行调整。1至10阶矩阵的RI值见表1。

　　(4)层次总排序。层次总排序是计算同一层次所有因素对于最高层相对重要性数值，这一过程是从最高层向最低层逐层进行的。采用层次分析法确定权重W总(表2)。

　　1.4模糊综合评价

　　将权重向量W总和模糊矩阵R合成，得到模糊综合评价结果矩阵：

　　B=W总・R(5)

　　构造水质标准类别矩阵S=[1，2，3，4，5]，则计算模糊综合评价指数为

　　FCL=B・S(6)

　　2基于模糊综合评价法的浑河大伙房水库上

　　游段水质评价

　　针对浑河上游北杂木、苏子河古楼和社河台沟3个断面，选取7种实测污染因子作为水质评价集，即评价集U={溶解氧，氨氮，高锰酸盐指数，BOD5，总氮，总磷，粪大肠菌群}。以《地表水环境质量标准》(GB 3838-2002)评价标准为依据建立评价标准集V={Ⅰ，Ⅱ，Ⅲ，Ⅳ，Ⅴ}，采用基于层次分析法的模糊综合评价法进行水质评价。

　　2.1建立模糊矩阵

　　应用式(1)-式(4)计算得到各断面的模糊矩阵。以2011年为例，北杂木、古楼和台沟断面模糊矩阵R分别为

　　2.2评价指标的权重

　　根据层次分析法将单排序分为A、B、C三个层次(表3)。

　　采用层次分析法确定污染因子权重，计算结果见表4。

　　总排序一致性检验：CI=∑3i=1bi(CI)i=0.0004，CR=0.000 3<0.1，即总排序满足一致性。

　　因素集U={溶解氧，氨氮，高锰酸盐指数，BOD5，总氮，总磷，粪大肠菌群}的权重集W总={0.035，0.217，0.088，0.111，0.011，0.230，0.308}。

　　2.3模糊综合评价

　　采用2004年-2012年水质监测数据，运用模糊综合评价法计算浑河大伙房水库上游北杂木、古楼及台沟断面全年及丰、枯水期的水质模糊综合评价指数FCL及水质级别，结果见表5。

　　由表5可见，北杂木断面全年水质类别均为Ⅱ类，丰水期水质类别除2011年为Ⅰ类以外其他年份均为Ⅱ类，枯水期水质类别除2012年为Ⅲ类以外其他年份均为Ⅱ类，进一步分析发现该断面当年枯水期总磷污染平均浓度为0.28 mg/L，超过国家地表水Ⅲ类水质评价标准80%。古楼断面全年水质类别均在Ⅱ类水质以上，自2008年起水质类别上升为Ⅰ类(2010年除外)，水质发展趋势良好，进一步分析发现该断面2010年丰水期时氨氮污染浓度为0.71 mg/L，超过国家地表水Ⅱ类水质标准42%;高锰酸盐指数污染浓度为5.65 mg/L，超过国家地表水Ⅱ类水质标准41.25%;总氮污染浓度为3.03 mg/L，超过国家地表水Ⅴ类水质标准51.5%，且这三种污染均大于古楼断面历年监测浓度。台沟断面全年水质类别为Ⅱ类，2008年转为Ⅰ类，水质总体趋势良好，丰水期在2009年由Ⅱ类转为Ⅰ类，枯水期水质类别均为Ⅱ类。

　　3水质模糊评价指数预测

　　本文通过建立来水量与水质模糊综合评价指数的非线性回归模型，对水质模糊综合评价指数进行预测。

　　3.1浑河大伙房水库上游段丰、枯水期划分

　　根据浑河大伙房水库上游段2004年-2012年实测来水量，计算丰、枯水期多年平均月入库水量(图1)。由图1可见，大伙房水库上游丰水期为4月-9月，期间入库水量大于50 m3/s;枯水期为12月-2月，期间入库水量小于15 m3/s;平水期为3月、10月及11月，期间入库水量大于20 m3/s且小于30 m3/s[15]。 　　3.2非线性预测模型的建立

　　利用丰水期、枯水期日平均来水量与水质模糊综合评价指数计算值，建立二者的非线性回归模型(表6)。

　　3.3非线性预测模型回归分析

　　应用SPSS软件对表6各非线性预测模型进行回归分析及非线性相关分析，结果见表6、表7。可见，除丰水期台沟断面回归模型的R2为0.0817外，其他各时期各断面R2均呈现较显著回归水平，说明除丰水期台沟断面外各时期各断面非线性回归方程与评价离散值拟合度较高。从表7可见，来水量与古楼断面丰、枯水期水质模糊综合评价指数均呈显著正相关(分别为p=0.807，p=0.762)，与北杂木丰、枯水期及台沟断面枯水期水质模糊综合评价指数均呈正相关(分别为p=0.253，p=0.577，p=0.62)，与台沟断面丰水期水质模糊综合评价指数呈负相关(p=-0.197)，负相关的结果是由于评价指数值与拟合曲线差值较大造成的，并不影响对台沟断面丰水期水质评价指数进行预测。

　　3.4水质预测

　　应用非线性预测模型对各断面丰水期、枯水期水质模糊综合评价指数进行预测，从结果(表8、表9)可见，预测值与计算值差值百分比均在20%以内[16]，两者比较吻合。北杂木断面、古楼断面、台沟断面丰水期预测指数分别为1.76、1.67、1.66，水质评价等级均为Ⅱ类;枯水期预测指数分别为3.32、1.59、2.5，水质评价等级分别为Ⅲ、Ⅱ、Ⅲ类。

　　4结论

　　(1)北杂木、古楼、台沟断面全年水质类别均在Ⅱ类以上，丰水期各断面水质类别均在Ⅱ类以上，枯水期除2012年北杂木断面水质类别为Ⅲ类以外，其他各断面各年均在Ⅱ类以上。

　　(2)通过建立来水量与水质模糊综合评价指数的非线性回归方程，发现除台沟断面丰水期水质模糊综合评价指数与来水量相关性较差以外，其他各断面各时期均呈正相关。

　　(3)北杂木、古楼、台沟断面丰水期水质预测评价等级均为Ⅱ类;枯水期水质预测评价等级分别为Ⅲ、Ⅱ、Ⅲ类。

　　需要指出的是，应用本文的模糊综合评价法进行水质预测时，数据系列越长回归系数R2越接近于1，预测结果越精确。

　　参考文献(References)：

　　[1]甘霖，张强.指数评价法在川北某村地下水水质评价中的应用[J].职业与健康，2009，25(24)：2670-2672.(GAN Lin，ZHANG Qiang.Application of index evaluation method in evaluation of ground water quality in a certain village of Northern Sichuan[J].Career and health，2009，25(24)：2670-2672.(in Chinese))

　　[2]曾永，樊引琴，王丽伟.水质模糊综合评价法与单因子指数评价法比较[J].人民黄河，2007，29(2)：64-65.(ZENG Yong，FAN Yin-qin，WANG Li-wei.Water quality fuzzy comprehensive evaluation method compared with single factor index evaluation method[J].The people of the Yellow River，2007，29(2)：64-65.(in Chinese))

　　[3]徐力刚，叶昌，张奇.基于模糊模式识别的地下水水质综合评价研究[J].水文地质工程地质，2011，38(5)：7-12.(XU Li-gang，YE Chang，ZHANG Qi.Application of fuzzy pattern recognition for the comprehensive assessment of groundwater quality[J].Hydrogeology Engineering geology，2011，38(5)：7-12.(in Chinese))

　　[4]李晴新，朱琳，陈中智.灰色系统法评价近海海洋生态系统健康[J].南开大学学报，2010，42(1)：39-43.(LI Qing-xin，ZHU Lin，CHEN Zhong- zhi.Ecosystem health assessment of marine sublittoral ecosystem by grey system method[J].Journal of Nankai University，2010，42(1)：39-43.(in Chinese))

　　[5]赵天燕，贺方升，侯棋棕.一种改进的灰关联分析法在水质评价中的应用[J].安全与环境学报，2005，5(1)：26-29.(ZHAO Tian-yan，HE Fang-sheng，HOU Qi-zong.Application of improved gray relation analysis method for water quality assessment in aquatic system[J].Journal of safety and the environment，2005，5(1)：26-29.(in Chinese))

　　[6]向娜.基于神经网络和人工蜂群算法的水质评价和预测研究[D].广州：华南理工大学，2012.(XIANG Na.Exploration for water quality assessment and prediction based on neural networks and artificial bee colony algorithm[D].Guangzhou：South China university，2012.(in Chinese)) 　　[7]姚焕玫，黄仁涛，刘洋，等.主成分分析法在太湖水质富营养化评价中的应用[J].桂林工学院学报，2005，25(2)：248-251.(YAO Huan-mei，HUANG Ren-tao，LIU Yang，et al.Principal component analysis of water high eutrophication evaluation in Taihu Lake[J].Journal of Guilin institute of technology，2005，25(2)：248-251.(in Chinese))

　　[8]胡成，苏丹.综合水质标识指数法在浑河水质评价中的应用[J].生态环境学报，2011，20(1)：186-192.(HU Chen，SU Dan.Application of comprehensive water quality identification index in water quality assessment of Hun River[J].Journal of ecological environment，2011，20(1)：186-192.(in Chinese))

　　[9]潘术香，李莲芳，张宝莉.利用标准类别指数评价法评价北京主要水系河流水质[J].农业环境科学学报，2004，23(3)：560-564.(PAN Shu- xiang，LI Lian-fang，ZHANG Bao-li.Utillizing standard category index method to assessthe water quality of main rivers of Beijing[J].Agricultural Journal of Environmental Science，2004，23(3)：560-564.(in Chinese))

　　[10]韦璐，江敏，余根鼎，等.凡纳滨对虾养殖塘叶绿素a与水质因子主成分多元线性回归分析[J].中国水产科学，2012，19(4)：620-625.(WEI Lu，JIANG Min，YU Gen-ding.Exploring the influence of water quality parameters on chlorophylla in Litopenaeus vannamei culture ponds using the method of multiplelinear regression of principal components[J].Agricultural journal of environmental science，2012，19(4)：620-625.(in Chinese))

　　[11]魏智宽，蒋世云，李少旦.灰色理论在龙江突发镉污染水质预测中的应用[J].水资源与水工程学报，2013，24(3)：135-141.(WEI Zhi -kuan，JIANG Shi-yun，LI Shao-dan.Application of grey theory in water quality prediction of sudden cadmium pollution in Longjiang River[J].Journal of water resources and water engineering，2013，24(3)：135-141.(in Chinese))

　　[12]梁楠.基于人工神经网络的水质预测及MATLAB实现[D].西安：长安大学，2007.(LIANG Nan.Research on artificial neural networks for water quality forecasting and realization on the MATLAB[D].Xi′an：Chang′an University，2007.(in Chinese))

　　[8]陈坚，娄红岩，李娟.大口径自由侧翻式拍门结构和撞击力计算[J].灌溉排水学报，2008，27(4)：44-47.(CHEN Jian，LOU Hong-yan，LI Juan，et al.Structure and impact force calculating of the Side-turn-over flap value of large diameter[J].Journal of Irrigation and Drainage，2008，27(4)：44-47.(in Chinese))

　　[9]王志祥，李端明，李娜，等.泵站不同类型拍门的开启角和撞击力分析[J].中国农村水利水电， 2012(7)：87-90.(WANG Zhi-xiang，LI Duan-ming，LI Na，et al.Opening angle and impact force analysis of different types flap valve in pumping station[J].China Rural Water and Hydropower，2012(7)：87-90.(in Chinese))

　　[10]Anthony Spiteri Staines.Modeling and analysis of a cruise control system[J].International Journal of Electrical and Computer Engineering，2008，3(10)：652-656.

　　[11]Demmin DM，Hallman Eric.Power take-off(PTO) safety[R].Department of Agricultural and Biological Engineering，Cornall University，New York，2009.

　　[12]杨帆，周济人，刘超.泵装置拍门水力损失数值模拟与试验[J].农业机械学报，2011，42(9)：108-112.(YANG Fan，ZHOU Ji-ren，LIU Chao.Numerical simulation and experiment of hydraulic loss of pump unit flap valve[J].Transactions of the Chinese Society of Agricultural Machinery，2011，42(9)：108-112.(in Chinese))

　　[13]成立，刘超.大型泵站水力稳定性探讨[J].南水北调与水利科技，2009，7(2)：75-77.(CHENG Li，LIU Chao.A discussion on the hydraulic stability of large pumping station[J].South-to-North Water Transfers and Water Science & Technology，2009，7(2)：75-77.(in Chinese))

　　[14]何志坚，叶锋.侧向式拍门流态分析及能效估算[J].陕西水利，2013：119-120.(HE Zhi-jian，YE Feng.Flow analysis and efficiency estimation of lateral flap valve[J].Shanxi Water Resources，2013：119-120.(in Chinese))

　　[15]朱红耕，冯汉民.拍门水头损失系数试验研究[J].排灌机械，1994(3)：49-51.(ZHU Hong-gen，FENG Han-min.Experimental study on hydraulic loss coefficient of flap valve[J].Drainage and Irrigation Machinery，1994(3)：49-51.(in Chinese))

——本文由丰临科技整理发布，内容供参考，如有侵权，请联系删除，谢谢！上海丰临科技有限公司为你提供浊度仪(浊度计)、在线浊度仪、余氯仪、余氯分析仪、工业在线pH计、cod测定仪、pH计等多种水质检测仪，水质分析仪，欢迎您前来选购，丰临科技竭诚为您服务！