长江水质的评价和预测大学生数学建模竞赛.doc

长江水质的评价和预测 摘要 本文通过对长江近几年的相关数据的分析，建立了长江目前水质情况的评价模型、长江流域主要污染物的污染源模型，预测长江未来10年的水质变化模型，并给出了治理长江污染的建议。 针对问题一：首先：利用MATLAB建立长江流域17个城市水质主要监测项目（PH、DO、CODMN、NH3-N）图表模型；然后：分析长江流域17个城市主要监测项目，对各地区水质情况分别做出分析；最后：对长江流域水质污染情况给出评价。 针对问题二：首先：由“相邻观测站之间的水流时间=相邻观测站距离÷水流速度”得出水由上一观测站到下一观测站的水流时间；然后：由“高锰酸盐量（或氨氮的量）=相邻观测站之间水流量×高锰酸盐量（或氨氮的量）指数”分别得出单位时间内高锰酸盐量和氨氮的量；再由“下游剩余高锰酸盐量（或氨氮的量）=高锰酸盐量（或氨氮的量）×（1-高锰酸盐量（或氨氮的量）×相邻观测站之间水流时间段内的高锰酸盐量（或氨氮的量）的降解系数）”得出上游高锰酸盐和氨氮经过一定时间的降解到下游时剩余高锰酸盐量和氨氮的量；最后：求得长江流域各观测站的月平均排放高锰酸盐的量和氨氮的量，利用MATLAB绘制高锰酸盐的量和氨氮的量各站点的月排放平均值的柱状图，通过比较得出锰酸盐和氨氮的主要污染源。 针对问题三：首先：由“废水排放浓度=废水排放总量÷长江总流量”求得长江前10年各年废水排放浓度；然后：利用MATLAB拟制长江废水排放浓度函数；最后：对长江未来10年各年废水排放浓度做出预测，绘制前10年的废水排放浓度及未来10年的拟合废水排放浓度曲线图。 针对问题四：首先：由“需处理废水量=长江总流量×（IV类水百分比+V类水百分比）×80%+长江总流量×劣V类水百分比”求得长江前10年各年需处理的污水量；然后：利用MATLAB拟制需处理污水量函数；最后：对长江未来10年各年需处理污水量做出预测，绘制前10年的需处理污水量及前后20年的拟合需处理废水量曲线图。 针对问题五：在上述4问的基础上，对长江水污染的情况做出可行性的建议。 关键词：图表模型 拟制函数 拟合曲线 MATLAB 柱状图 一、问题重述 水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。长江是我国第一、世界第三大河流，近几年来，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。 在长江沿线设有观测站以监测长江的水质，现有长江沿线17个观测站（地区）近28个月4个主要水质指标的检测数据，以及干流上7个观测站近一年多的基本数据（站点距离、水流量和水流速）（见附件3）。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。一般说来，污染物在水环境中可通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。该降解系数为江河自然净化能力指标。事实上，长江干流的自然净化能力可以认为是近似均匀的。给出“1995~2004年长江流域水质报告”的主要统计数据（见附件4）及国家规定4个主要水质指标的标准限值 三、模型假设 假设降解作用在整个流域中是平稳进行的，取整个流域的降解系数恒为0.2（单位：1/天）。 假设一个观测站（地区）的水质污染主要来自于本地区自身的排污和上游污水。 考虑污染源就在观测站点当地，即有多少个观测站点就有一个包含相应个数的污染源候选地的候选集。 忽略各污染物内部反应与相互作用项，如生物化学中的生长与降解变化。 未来十年的污染速率是一样的。 四、符号说明 x 年份（将年份看作数字，如“1995-2004”为“1-20”） y 废水排放浓度 水在第个时间段从第个观测站到个观测站点的时间（） 第个观测站在第个时间段内高锰酸盐的量 第个观测站在第个时间段内氨氮的量 第个观测站在第个时间段内水的流量 第个观测站在第个时间段内水的流速 第个观测站在第个时间段内高锰酸盐的指数 第个观测站在第个时间段内氨氮的指数 第个观测点距初始站点的距离 第个观测站在时间段内上游高锰酸盐经过降解到下游时剩余的量 第个观测站在时间段内上游氨氮经过降解到下游时剩余的量 第个观测站在时间段自身排放高锰酸盐的量 第个观测站在时间段自身排放氨氮的量 降解系数 长江总流量 需处理污水量 Ⅳ类水百分比 Ⅴ类水百分比 劣Ⅴ类水百分比 废水排放总量 五、模型的建立与求解 5.1、问题一 由附件3，分别建立长江流域17个城市在2003.6-2005.9（28个月）的主要监测项目图表模型。 图表说明： 横轴： 5.1.1四川攀枝花2003.6--2005.9长江水质监测图（附录1.Fun1） 可