长江水质的评的价和预测24628.docVIP

长江水质的评价和预测 一．摘要：本文在参考一些数据,文献的基础上对长江水质以及变化趋势综合分析并建模，对母亲河的水质做出了一个客观的评价并对水质的变化趋势做出了预测，针对问题一，水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。专家们呼吁：“以人为本，建设文明和谐社会，改善人与自然的环境，减少污染。”长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。 根据给定的统计数据预测长江未来（10年）水质的污染变化趋势。并在保证比例控制的条件下，做出每年合理的污水处理计划，建立相应的模型，以达到使第IV类与V类水保持在一定比例之内并满足没有劣V类水的目的。解决长江水质污染的问题已迫在眉睫，最后根据长江水质污染的现实问题提出解决问题的可行性建议和意见。 附表: 《地表水环境质量标准GB3838—2002）项目标准限值 单位：mg/L 序 号 ? 分 类标准值 Ⅰ类 Ⅱ类 Ⅲ类 Ⅳ类 Ⅴ类 劣Ⅴ类 1 溶解氧(DO)??????????≥ 7.5 （或饱和率90%） 6 5 3 2 2 高锰酸盐指数(CODMn) ≤ 2 4 6 10 15 ∞ 3 氨氮（NH3-N）????? ≤ 0.15 0.5 1.0 1.5 2.0 4 PH值（无量纲） 6---9 问题一 四．建模过程 1）符号设置： 溶解氧的浓度的倒数值（1/DO） 高锰酸盐的指数浓度（CODMn） 氨氮浓度（NH3－N） 2）模型假设： 不考虑江水的蒸发量、渗透量以及降雨量对河水总量的影响 2、各种污染物流入长江干流后都能与江水充分混合。 3、入河排污口水量与水质变化稳定 4、四川攀枝花上流的污染忽略。 5、检测数据真实可靠，具有代表性。 6、长江干流及支流的自然净化能力、同一条支流上观测站的水流量是均匀 的 ，每一个观测站的检测范围相同。 7、高锰酸盐指数和氨氮的降解系数为0.2（单位：1/天）。 8、污染物排到水流中后，其量的变化只与水的自然净化有关。 9、一个观测站的水质污染主要来自于本地区的排污和上游的污水。 10、废水是均匀排放的，即每一个时刻排放的废水量相同。 11、所有废水的排放量都未经是处理，全为直排。 3）模型的建立与求解： （一）问题一：基于主成分分析的水质综合评价模型 1、主成分分析原理及计算过程 （1）将样本数据标准化，得到标准化后的数据矩阵为 X=,其中：=，， ，为原始样本数据。 （2）计算相关系数矩阵 R= (3)计算R的特征方程的m个非负的特征值. (4)计算对应特征值的相应的特征向量，m个特征向量成 的特征向量矩阵 (5)由特征向量C组成个新因子，有. (6)选择p(pm)个主分量,这p个主分量的累积方差贡献率大于等于85%以上,即 . (7)根据计算p个主成分的值。 （8）计算样本综合评价值，式中。 2、长江水质的主成分分析评价 我们选用3个指标对长江水质进行评价，这三个指标分别是：溶解氧的浓度的倒数值（），高锰酸盐的指数浓度（）， 氨氮浓度（）。样本数据如下表一所示： （三）问题三：基于灰度预测与时间序列分析法的未来水质污染变化趋势模型 根据问题分析，我们将干流不同时期可饮用水河长所占比例作为评价长江污染情况。 1 GM（1，1）模型的建立与计算过程 时间序列有X有10 个观察值，X={ X(1), X(2) ... X(10)，通过累加生成新序列X={ X(1), X(2) ... X(10)， 则GM（1，1）模型相应的微分方程为： +x= 其中：α 称为发展灰数；μ 称为内生控制灰数。 2、设为待估参数向量，，可利用最小二乘法求解。解得： 求解微分方程，即可得预测模型： ， 利用dsp进行求解，得到灰度预测下的1995 年到2014 年后的结果为 Y={ Y(1), Y(2) ... Y(20)}={93.1， 79.78， 78.246， 76.741， 75.264， 73.817 ， 72.397， 71.004， 69.638， 68.298， 66.985， 65.696， 64.432， 63.1， 61.977 ， 60.785，59.616， 58.469， 57.344，56.241} 得到数据如下图： 图 全流域枯水期实测值