2005高教社杯全国大学生数学建模竞赛题目A 长江水质的评价和预测

长江水质的评价和预测

摘要

文章在已有数据的基础上，建立了水质依靠流量、流速和降解系数的数学模型，找出了污染源的所在地。建立一元线性回归模型，对后十年污水治理做出了预测。利用Matlab，C语言程序进行求解。得出了有关结论。

针对问题一，根据03、04年长江流域水质报告表，对长江近两年多的水质情况做出定量的综合评价。对每一个地区在近两年的28个月中的水质情况进行统计，找出该地区污染的种类及该种类污染出现的频率。以此为依据分析各地区水质的污染状况。

针对问题二，根据主要污染物在各个观测点的观测数据，建立了水质依靠流量、流速和降解系数的数学模型。对长江干流沿岸各个地段的排污量进行统计，找出了主要污染源所在地区：长江中游湖北宜昌至湖南岳阳段。

针对问题三，根据各个年份废水排放量总量，采用一元线性回归模型找出废水排放量总量与年份之间的关系。根据水文年支流和干流的相关数据和各年长江总流量和废水排放量，得出长江总流量中废水排放量的比例，利用matlab对1995-2004年长江总流量中废水排放量的比例拟合（不考虑1998年特大洪水），对未来十年废水排放量占长江总流量比例进行预测。从预测结果中，发现污水百分比呈逐年上升的趋势（从2005年的3.27%到2014年的6.24%），由此说明长江污水的处理迫在眉睫。

针对第四问，依照过去10年的Ⅳ类、Ⅴ类水和Ⅵ水的统计数据，通过数据拟合构建了一元线性回归模型、预测的未来十年Ⅳ类、Ⅴ类水和Ⅵ水占长江总水量的百分比。引入流量的概念，得到长江的总水量HS。由治理污水的标准建立分段函数。从而求出未来十年每年需要处理的污水量。

针对第五问，提出了解决长江水质污染问题的从四方面着手的方案：沿江工厂的整治，民众意识的唤醒，上游植被的保护，以及法律的硬性要求。

一问题重述

2004年10月“保护长江万里行”考察团，对长江沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面。为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”，并发出了“拿什么拯救癌变长江”的呼唤。

通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。江河自身对污染物都有一定的自然净化能力，反映江河自然净化能力的指标称为降解系数。长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.1~0.5之间，比如可以考虑取0.2(单位：1/天)。

请根据，长江沿线17个观测站近两年多主要水质指标的检测数据，干流上７个观测站近一年多的基本数据（站点距离、水流量和水流速），“1995~2004年长江流域水质报告”给出的主要统计数据，以及《地表水环境质量标准》。

研究以下问题：

（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。

（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?

（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。

（4）根据你的预测分析，如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,那么每年需要处理多少污水？

（5）你对解决长江水质污染问题有什么切实可行的建议和意见。

二问题分析

针对问题一，要分两步来解决。第一步，我们根据03、04年长江流域水质报告表，对长江近两年多的水质情况做出定量的综合评价。第二步，分析近两年多长江各地区水质的污染状况。为此我们先对每一个地区在近两年的28个月中的水质情况进行统计，再找出该地区污染的种类及该种类污染出现的频率。

第二个问题分为两个部分：一是研究、分析长江干流近一年多主要污染物高锰酸盐污染源主要在哪些地区。二是研究、分析长江干流近一年多主要污染物氨氮的污染源主要在哪些地区。解决这两部分问题所用的方法是一样的。要想找出污染源，就要对长江干流沿岸各个地段的排污量进行统计。排污量高的地段就是主要污染源所在地区。

对于问题三：首先筛选取出水文年支流和干流的相关数据和各年长江总流量和废水排放量，对1995-2004年长江总流量中废水排放量的比例拟合，对每年废水排放量占长江总流量比例进行预测。通过预测出的数据分析得出废水排放量占长江总流量比例变化情况，由此对未来十年长江干流和支流污染量的污染量的变化趋势作出总体的分析说明。

问题四要求我们预测分析如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,每年需要处理多少污水。流量除以流速可以得到过水断面的面积。干流处有七个检测点，求出每个检测点的过水断面，求平均值。此值即为长江的平均过水断面。用过水断面的平均值乘以平均河长可以得到长江水的总体积。为了求解需要处理的污水量，我们还要求出需要处理的污水占总水量的百分比。此时问题可分两类进行讨论。第一类，Ⅳ加Ⅴ水所占长江水的百分比未超过20%。此时为了使长江水质达到题目中的要求，只需对劣Ⅵ水进行处理。所需处理的污水量即为长江中劣Ⅵ水的总量，用劣Ⅵ水占长江总水量的百分比乘以长江的总水量即可得到。第二类，Ⅳ、Ⅴ水所占长江水的百分比超过20%。此时为了使长江水质达到题目中的要求需要做两方面的处理。一是将Ⅳ、Ⅴ水进行处理，使其在长江水中所占百分比降至20%。二是对劣Ⅵ水进行处理。所需处理的污水总量即为所需处理的Ⅳ、Ⅴ水的量与长江中劣Ⅵ水的总量之和。用Ⅳ、Ⅴ水占长江总水量的百分比之和，减去20%再乘以长江的总水量，即可得到所需处理的Ⅳ、Ⅴ水的量。用劣Ⅵ水占长江总水量的百分比乘以长江的总水量即可得到长江中劣Ⅵ水的总量。用一元线性回归模型对Ⅰ、Ⅱ、Ⅲ类水所占百分比与年污水排放量，Ⅰ、Ⅱ、Ⅲ类水所占百分比与万倍年污水浓度，Ⅰ、Ⅱ、Ⅲ类水所占百分比与年污水浓度的倒数，分别进行拟合。选取拟合度最高的模型来预测未来十年长江Ⅰ、Ⅱ、Ⅲ类水所占百分比。用一元线性回归模型对Ⅳ、Ⅴ类水所占百分比与年污水排放量，Ⅳ、Ⅴ类水所占百分比与年污水占长江总水量百分比，Ⅳ、Ⅴ类水所占百分比与年污水占长江总水量百分比的倒数，分别进行拟合。选取拟合度最高的模型来预测未来十年长江Ⅳ、Ⅴ水所占百分比。

针对第五问，提出了解决长江水质污染问题的从四方面着手的方案：沿江工厂的整治，民众意识的唤醒，上游植被的保护，以及法律的硬性要求。