《毕业设计论文-数学建模论文长江水质的评价与预测.doc》由会员分享，可在线阅读，更多相关《毕业设计论文-数学建模论文长江水质的评价与预测.doc（29页珍藏版）》请在金锄头文库上搜索。
1、高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的
2、话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：长江水质评价预测模型摘要本文通过对各监测点水质指标数据的分析，建立了水质综合评价和预测模型，得到了水质现状以及各类水质和废水排放量的发展趋势，并以此为基础预测出未来十年需处理的
3、废水量。对水质进行综合评价时，以监测平均值作为水质总体情况的评价指标。将所有监测值按照枯水期，丰水期，平水期分为三类，求出每个时期各指标的平均值，再确定出各个时期的权重，最终得到监测平均值。结果显示，50%以上的水处于类，存在6%左右的劣类水。就各干流和支流的水质情况而言，支流的污染情况更为严重。为寻找长江干流氨氮和高锰酸盐污染源，引入差分方程，描述上游流入当地的污染量和当地排放的污染量之和等于当地污染物的监测值这一关系。再考虑降解作用，求解出各个监测点的当地排污量，得到氨氮和高锰酸盐污染物的主要污染源为重庆朱沱至湖南岳阳城陵矶这一段水域的结论。为预测出未来十年水质的发展趋势，本文通过建立灰色
4、预测模型，预测出可饮用水，+类水，劣类水这三项的百分比以及每年的废水排放量。用后残差检验预测精度，发现精度偏低。于是采用吐故纳新的思想处理首尾数据，以修正模型，得到较为满意的水质预测值。预测显示，十年之后长江流域的可饮用水比例将不足50%，而全流域的劣类水将上升到50.3%，废水排放量呈逐年递增趋势。预测要使水质达标未来十年每年所需处理的废水量，先根据灰色预测模型所得的三类水质百分比和废水排放量，然后利用SPSS拟合出废水量关于上述三类水质百分比的多元线性函数。接着按照所规定的水质比例控制标准，得到新的三类水质百分比，代入拟合函数求得废水排放量，进而求得所需处理的废水量。结果表明，需处理的废水
5、量以较大幅度逐年递增，到2014年将达到65.3737亿吨。最后，对所建模型进行评价和改进，并且就水质评价和预测的结果给出了具体的建议。关键词：长江水质 综合评价 差分方程 灰色预测 废水处理 一、问题重述长江水质的污染程度日趋严重，2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。通常认为一个观测站的水质污染主要来自当地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，长江干流的自然净化能力可以认为是近似均匀的。现有长江沿线17个
6、观测站（地区）近两年多主要水质指标的检测数据，干流上个观测站近一年多的基本数据（站点距离、水流量和水流速）以及“19952004年长江流域水质报告”给出的主要统计数据。我们需要研究如下问题：(1)对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。(2)研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?(3)假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。(4)根据以上的预测分析，如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,那么每年
7、需要处理多少污水？ (5)对解决长江水质污染问题 提出实际可行的建议和意见。二、问题分析对长江水质的评价和预测，可通过分析各污染物的当前浓度和在一定时期内的变化规律，给出各监测站点水质所处于的平均等级，并找出主要污染源，达到评价的目的；再根据过去十年的数据，建立预测模型，用适当的函数对各类水质的比例变化以及废水排放量进行预测，在此基础上，求解出控制类和类水的比例在20%以内且不含劣类水需处理的废水量。在进行水质的综合评价时，考虑到枯水期，丰水期和平水期的水流量有较大的差异，不能将所有的月份各污染物浓度相加平均。于是可先将所有统计数据按上述三个时期进行分类，求出每个时期各污染物的平均浓度指标，再
8、设置不同的权重将三个时期的浓度综合，最后得到各个监测站点的水质等级。在就高猛酸盐和氨氮的浓度寻找长江干流的主要污染源时，关键是确定污染源的认定指标。监测到的污染物浓度值是上游和当地排污的综合结果，不足以用来衡量一个地区是否为污染源，因此需要找出当地的实际排污量，用此作为污染源的认定指标。鉴于观测点之间是逐级影响的关系，可建立差分方程模型对问题进行求解。在对未来十年的水质情况进行预测时，发现各年数据之间有一定的变化趋势，但并不明显，若直接进行函数的拟合会使结果显得粗糙，于是可考虑用灰色预测的原理建立模型，最后得到未来十年的预测数据。根据预测求解未来十年每年需处理的废水量时，首先得明确废水的概念，
9、考虑到、 类是可饮用水，那么废水至少为类水，又因为江水对废水中的污染物浓度起到了一个稀释的作用，所以废水本身的污染物浓度一定大于与江水混合后的浓度。不过对劣类水的各项指标的浓度难以确定，所以考虑换一种角度解决问题。可以通过找出每年的废水排放量与各类水质百分比之间的关系，用函数表示出来，然后再通过各类水质需达到的百分比要求反过来计算废水量，两次废水量之间的差值便是需要处理的废水量。三、基本假设1.假设处于同一水文期的各月水流量相等。2.假设两个站点之间的各支流流量相同。3.假设两干流监测点之间的支流流量为两干流站点的流量之差取平均。4.假设攀枝花为所研究水域的源头，污染物浓度仅由当地排污决定。5
10、.假设污染物不沉淀，且与江水均匀混合6.假设干流和支流的降解系数取为0.27.假设处理后的废水变为可饮用水。四、变量说明:第个监测站点第项水指标的平均浓度(mg/l)：第个监测站点第项水指标在第类水文期的平均浓度(mg/l) ：第类水文期水指标浓度所占的权重：干流监测站点的污染物监测值：支流监测站点的污染物监测值：自上游流入监测点的污染物量：相邻两干流监测站点和之间河段的平均流速（km/天）：监测站点和之间的距离（km）：污染物从监测站点流入的时间：降解系数（g/天）水文期：指枯水期，丰水期，平水期水质指标：附表中地表水环境质量标准列出的四个监测项：溶解氧、高锰酸盐指数、氨氮、PH值五、模型建
11、立和求解5.1水质综合评价模型对题目中所给的17个监测站点逐一考虑。对某一站点来说，共有28次监测值，按照枯水期，丰水期和平水期的月份划归，将28组数据分为三类。在假设处于同一水文期的各月水流量相等的情况下，可将同一水文期的四项主要污染物浓度相加求出平均值。最后再确定出三个时期的权重，加权平均便可求得这一站点在两年里的各项平均水质指标。某站点的平均浓度的表达式如下：其中，表示第个监测站点第个指标的平均浓度，表示第个监测站点第个指标在第类水文期的平均浓度，表示第类水文期指标浓度所占的权重。i取1,2,3分别对应枯水期、丰水期、平水期，m取117代表17个监测点,n取14代表4项水质指标。 分类结
12、果如下：枯水期：2004年14月，2005年14月；(共8个月)丰水期：2003年610月，2004年510月，2006年59月；(共16个月)平水期：2003年1112月，2004年1112月；(共4个月)接着确定各水文期所占的权重，求出各个水文期在一年内的平均水流量，将水流量归一化便可得到各个水期的权重。在此我们假设处于同一水文期的各月水流量相等。得出的三个水文期的权重分别为：枯水期：0.2322；丰水期：0.5544；平水期：0.2134。求解结果如下：表1各监测站点各指标的平均值 （单位：mg/L） 全年平均值地点PHDOCODMnNH3-N水质类别四川攀枝花8.24629.17162
13、.45320.16842重庆朱沱7.91898.98972.05260.33072湖北宜昌南津关7.76318.52662.8530.26052湖南岳阳城陵矶7.83148.65043.76580.32982江西九江河西水厂7.43067.76872.4270.16062安徽安庆皖河口7.45257.44692.54710.22712江苏南京林山7.66317.47682.11930.12792四川乐山岷江大桥7.52155.63515.20540.88343四川宜宾凉姜沟8.06598.99332.76370.4342四川泸州沱江二桥7.68136.86923.27640.74733湖北丹江口胡家岭7.86229.37581.96770.09371湖南长沙新港7.09917.10632.42120.93583湖南岳阳岳阳楼7.72578.31134.17510.38033湖北武汉宗关7.9447.47853.38170.19632江西南昌滁槎7.1765.90332.44325.01176江西九江蛤蟆石7.6393
