毕业论文设计--长江水质的评价和预测建模论文.doc

高教社杯全国大学生数学建模竞赛

承 诺 书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B
我们的参赛报名号为（如果赛区设置报名号的话）：
所属学校（请填写完整的全名）：
参赛队员 (打印并签名) ：1.
2.
3.
指导教师或指导教师组负责人 (打印并签名)：

日期： 年 月 日

赛区评阅编号（由赛区组委会评阅前进行编号）：
高教社杯全国大学生数学建模竞赛

编 号 专 用 页

赛区评阅编号（由赛区组委会评阅前进行编号）：

赛区评阅记录（可供赛区评阅时使用）：

评
阅
人

评
分

备
注

全国统一编号（由赛区组委会送交全国前编号）：

全国评阅编号（由全国组委会评阅前进行编号）：

长江水质的评价和预测

摘要
水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国江河水资源的保护和治理应是重中之重。本文主要研究了以下四个问题：长江水质的综合评价、主要污染源的确定、预测问题和污水处理问题。并以此对解决长江水质污染问题提出一些切实可行的建议和意见。
问题一，我们建立模型一：问题一，我们建立模型一：选择CODMn和NH3-N等三个指标，构造目标目标函数，应用数值方法结合蒙特卡洛法解决矩阵处理中穿插着非线性规划的问题，规划出各指标的权重系数，建立综合评估系统，得出2003年6月份至2005年9月份各地区水质的综合评估值。算得国标水质标准Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ类水的水质综合评估值分别为：5.0225 0.3120 0.2496 0.2079 0.1248 0.0832。详细结果见模型一求解部分，见表四（第8页）。
问题二，我们建立模型二：针对高锰酸盐和氨氮化合物，由于其在水中能自然降解，利用衰减理论计算各地区的排放量而得到长江干流近一年多高锰酸盐指数和氨氮的主要污染源。结论是：长江干流近一年多主要污染物高锰酸盐指数的污染源主要在湖南岳阳、湖北宜昌、江苏南京；主要污染物氨氮的污染源主要在湖北宜昌、湖南岳阳、重庆朱沱。
问题三，我们建立模型三：由附表所给的数据进行曲线似合，得到每类水所对应的函数关系，进而进行预测。结论是：2006年废水排放量为222亿吨，2010年的为429亿吨，2014年的为563亿吨，可见废水年排放量逐年上升，2006年至2015年的具体数据见表十（第17页）、表十一(第17页)和表十二（第18页)。
问题四，我们建立模型四：找出废水年排放量与各类别水河长、长江年总流量的一一对应关系，取污水密度1.7千克每立方米,在满足题目中所给的条件下，求出每年应处理的废水。2005年需处理的污水174.3亿吨，2009年的为266.7亿吨，2013年的为242.8亿吨，详细结果见表十七（第19页）
最后，我们结合以上四个问提出控制污染源、加强污水处理、整治污染河段和调整产业结构等四条切实可行的建议和意见。详见第20页。

关键字：综合评价；权重系数；蒙特卡洛法；预测；曲线拟合

1. 问题重述

水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。 长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。
现有长江沿线17个观测站（地区）近两年多主要水质指标的检测数据，以及干流上７个观测站近一年多的基本数据（站点距离、水流量和水流速）。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低，反映江河自然净化能力的指标称为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.1~0.5之间，比如可以考虑取0.2　(单位：1/天)。 “1995~2004年长江流域水质报告”给出了主要统计数据。以国标(GB3838-2002) 给出的《地表水环境质量标准》为标准。
现要研究以下五个问题：
（1）地表水中的4个主要项目标准限值见附表，依此标准对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。
（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?
（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。
（4）根据预测分析，如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,那么每年需要处理多少污水？
（5）对解决长江水质污染问题提出一些切实可行的建议和意见。

2.问题分析

2.1 问题1
要对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况，需要根据2003年6月份至2005年9月份长江流域主要城市水质检测报告，选择指标，建立合理的综合评价系统：按指标分别计算评估值，依据计算所得的评估值分析各地区水质的污染状况。主要问题为各指标权重的确定。

2.2 问题2
要研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要所在地，就要先确定长江干流从上游到下游各城市（向长江）的排污量。
根据附件3“长江干流主要观测站点的基本数据”中的部分数据，综合于表一：
表一：观测站点间的距离
观测站点
四川攀枝花
重庆朱沱
湖北宜昌
湖南岳阳
江西九江
安徽安庆
江苏南京
站点间距离
0
950
1728
2123
2623
2787
3251
注：（单位：km）。
可得地理分布。
主要问题为：如何依据2003年6月份至2005年9月份长江流域主要城市水质检测报告和长江干流主要观测站点的基本数据计算出近一年多每个地区主要污染物高锰酸盐指数和氨氮的排放量，干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要所在地即为排放量的前几名。水中氢离子浓度和溶氧量无需考虑。

2.3 问题3
假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。
要保证预测的准确性就需分时段对同一水类的河长进行预测，获得近两年的水文资料可进一步验证和修正模型。
需要解决的问题：预测模型。

2.4 问题4
标准即为未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水。根据问题三中的可靠预测，计算应处理的Ⅳ类、Ⅴ类水和劣Ⅴ类水的质量，这三个数据的和即为每年应处理的废水总量。
主要问题：每年长江总流量（单位：亿立方米）、废水排放总量 （单位：亿吨）与长江流域水质报告表有怎样的联系。

2.5 问题5
依据以上四个问题的分析结果，对解决长江水质污染问题提出切实可行的建议和意见。

3.模型假设与符号说明
3.1模型假设
（1）江水中的pH、DO、CODMn和NH3-N等四项指标相互独立。
（2）pH值对其水质的综合评估值无影响。
（3）对处理污水不影响长江的年总流量，长江的总河长。
（4）高锰酸盐指数和氨氮的降解系数取0.2　(单位：1/天)。
（5）不考虑特大自然水灾害，如98大洪水。
（6）对长江水质变化起决定性作用的是废水排放量，而与其它的的因素无关。
（7）取污水的密度为1.7千克每立方米。

3.2符号说明
y
表示第i个地方的综合评估值

表示第j种指标的的权重系数

表示第i个监测站水质的第j种指标值

表示第i个监测站水质的类别

表示所有研究对象第i种指标的最大值

表示个水质标准所给的第i类水的综合评估值
s1
表示所有相邻两类水的综合评估标准值的平方差的比重和

表示第i个监测站水质的类别值

表示干流观测站点i和i+1间的距离

表示干流观测站点i的水流速

表示干流观测站点i-1到i的水流时间

表示干流观测站点j中CODMn()和NH3-N()的浓度

表示干流观测站点j每秒排放CODMn()和NH3-N()的浓度

表示干流观测站点i的水流量

表示干流观测站点j每秒排放CODMn()和NH3-N()的质量

表示未来第j年需处理的ⅣⅤ类（i=1）和劣Ⅴ类(i=2)所占百分比

表示未来第j年需处理的ⅣⅤ类（i=1）和劣Ⅴ类(i=2)的河长

表示未来第j年水文年全流域的河长

表示未来第j年需处理的污水的河长

4.模型的建立与模型的求解

4.1 模型一：评价模型
4.1.1模型一的建立
为了解决问题一，我们选择DO值、CODMn值的倒数和NH3-N值的倒数等三项指标，建立综合评估值计算公式。将每个地区在每一月份的DO值、CODMn值的倒数和NH3-N值的倒数等三个指标值作为一个研究对象。按指标分别计算各个研究对象的综合评估值，即可对长江近两年多的水质情况做出定量的综合评价，并以此为依据分析各地区水质的污染状况。
江水中4个主要项目标准限值见表二：
表二：《地表水环境质量标准》中4个主要项目标准限值（单位：mg/L）

序
号
  分 类
            标准值
  项 目
Ⅰ类
Ⅱ类
Ⅲ类
Ⅳ类
Ⅴ类
劣Ⅴ类
1
溶解氧(DO)    ≥
7.5
（或饱和率90%）
6
5
3
2
　 0
2
高锰酸盐指数(CODMn) ≤
2
4
6
10
15
∞
3
氨氮（NH3-N）  ≤
0.15
0.5
1.0
1.5
2.0
∞
4
PH值（无量纲）
6---9

分析该表可认为江水中高锰酸盐指数和氨氮越低越好，溶解氧（浓度）越高越好。
我们定义综合评估值：分别将各地区在不同时间段内的三项指标数据进行处理，然后乘上相应的权重系数作和得到综合评估值；综合评估值越高，水质越好。水的类别值越高，水质越差，综合评估值越低。即若，则。
设第i个研究对象的DO值、CODMn值的倒数和NH3-N值的倒数等三个指标值分别为、和，第i个研究对象关于这三个指标的线性综合评估值计算公式为：

权重系数的确定是首先要解决的问题。通过比较不同权重系数下的Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水中相邻两类水的综合评估标准值的平方差的比重和确定合理的一组权重系数：
步骤一：将原始数据处理成DO值、CODMn值的倒数和NH3-N值的倒数等三个指标值
以DO值、CODMn值的倒数和NH3-N值的倒数等三个指标值对表 给出的水质标准数据进行处理，结果如下：
7.5000 0.5000 6.6667
6.0000 0.2500 2.0000
5.0000 0.1667 1.0000
3.0000 0.1000 0.6667
2.0000 0.0667 0.5000
步骤二：定义Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ类水的水质综合评估标准值
我们定义Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ类水的水质综合评估标准值、、、、和的计算：

（、、为所有研究对象第i种指标的最大值，所所给研究对象的影响）

之后用取代得到（i=1，2……6）
（将其DO值、CODMn值的倒数和NH3-N值的倒数等三个指标值乘以对应的权重，并经一定处理）

步骤三：规划出最好的权重系数取值
由于Ⅰ、Ⅱ类水的比重占绝大部分，考虑Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水的比重对相邻两类水的综合评估标准值的平方差的和的影响。
对于任何一组权重系数，可得Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ水中相邻两类水的综合评估标准值的平方差的比重和，将该比重和作为目标函数s1，便将问题转化为一个非线性规划问题。
为能更好地解决问题，我们有必要对算法做一些分析和介绍：以上问题的解决过程中，大部分工作是矩阵处理，通过比较不同权重系数下的Ⅰ、Ⅱ、Ⅲ类水中相邻两类水的综合评估标准值的平方差的比重和确定最佳的一组权重系数的过程，则是非线性规划问题。从整体思路上讲，既有矩阵处理，又有非线性规划问题，若用LINGO求解非线性规划问题，过程将十分复杂；况且在本模型中这二者是无法分割开来求解的，于是我们采用数值方法利用MATLAB软件求解。
采用数值方法利用MATLAB软件求解时，又遇到另一个问题：如何取值才能使权重系数的比例满足任意性呢？能否让权重系数的比例既满足任意性，又是最简比呢？
用蒙特卡洛法可以很好地解决上述问题：对权重取任意的一组0-1间的随机数，求这三个数的和，用该和数分别去除那组0-1间的随机数。这样做是为了保证权值和为1，即化为最简形式。将这组值代入Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ类水的水质综合评估标准值、、、、和的计算公式，就能求得目标函数在既定权重系数下的取值。
步骤四：各研究对象的综合评估值的计算
取出某一研究对象，将其数据进行步骤一、步骤二的处理。若该研究对象是属于Ⅰ类的，则将上一步处理结果乘以；若是Ⅱ类的，将上一步处理结果乘以……
4.1.2 模型一的解答：
本模型的MATLAB程序见附录 程序一（第 页），由“f2”进行了十四次仿真，（必需经过多次仿真，以作比较，去最好的）仿真结果如下：
s1
第1次仿真结果： 40.3773 0.9912 0.0080 0.0008
第2次仿真结果： 40.3614 0.9905 0.0079 0.0016
第3次仿真结果： 40.3695 0.9651 0.0347 0.0002
第4次仿真结果： 40.3557 0.9544 0.0452 0.0005
第5次仿真结果： 40.3742 0.9665 0.0335 0.0000
第6次仿真结果： 40.3893 0.9912 0.0086 0.0002
第7次仿真结果： 40.3630 0.9587 0.0410 0.0003
第8次仿真结果： 40.3705 0.9746 0.0249 0.0005
第9次仿真结果：40.3779 0.9859 0.0135 0.0006
第10次仿真结果：40.3784 0.9870 0.0124 0.0006
第11次仿真结果：40.3729 0.9963 0.0025 0.0012
第12次仿真结果：40.3299 0.9731 0.0245 0.0024
第13次仿真结果：40.3808 0.9806 0.0191 0.0002
第14次仿真结果：40.3833 0.9823 0.0175 0.0002
由此，取仿真结果的第6组数据：
=0.9912 =0.0086 0.0002
将上述数据带入程序一“f3”（第 页）求得
国标水质标准Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、劣Ⅴ类水的水质综合评估值分别为：
5.0225 0.3120 0.2496 0.2079 0.1248 0.0832

由此得到如下结论：
①国标水质标准及其综合评估值，见下表三：
表三 ：国标水质标准及其综合评估值

序
号
分 类
标准值
项 目
Ⅰ类
Ⅱ类
Ⅲ类
Ⅳ类
Ⅴ类
劣Ⅴ类
0
综合评估值 ≥
5.0225
0.3120
0.2496
0.2079
0.1248
0.0832
1
溶解氧(DO)  ≥
7.5
（或饱和率90%）
6
5
3
2
0
2
高锰酸盐指数(CODMn) ≤
2
4
6
10
15
∞
3
氨氮（NH3-N） ≤
0.15
0.5
1.0
1.5
2.0
∞
4
PH值（无量纲）
6---9

②各地区在2003年7月份至2005年9月份间每个月份的综合评估值详细结果如下表四：
表四：各地区各月综合评估值
点位名称
200306
200307
200308
200309
200310
200311
200312
四川攀枝花
0.2233
0.2114
4.2636
4.6265
0.2807
4.7838
0.3066
重庆朱沱
0.2745
0.2788
0.2792
0.2802
0.2945
0.3112
0.3396
湖北宜昌南津关
0.2038
0.2059
0.2171
0.3459
0.3884
2.9366
0.3144
湖南岳阳城陵矶
0.2112
0.2471
0.2585
0.2729
0.209
0.2804
0.2176
江西九江河西水厂
0.2022
0.2041
0.2246
0.2426
0.2501
0.282
0.2807
安徽安庆皖河口
0.2135
0.2196
0.2307
0.1801
0.2497
0.2563
0.283
江苏南京林山
0.2253
0.208
0.2112
0.216
0.225
0.2745
0.2921
四川乐山岷江大桥
0.0914
0.1007
0.0716
0.1375
0.1519
0.1053
0.1394
四川宜宾凉姜沟
0.2491
0.1757
0.2829
0.2977
0.3297
0.2688
0.1383
四川泸州沱江二桥
0.0875
0.1022
0.2337
0.2018
0.1845
0.1648
0.0827
湖北丹江口胡家岭
5.3602
5.1554
0.2948
0.2322
0.2749
0.2854
0.333
湖南长沙新港
0.1684
0.2262
0.0944
0.1749
0.1184
0.1603
0.1556
湖南岳阳岳阳楼
0.2045
0.2119
0.2716
0.2872
0.3254
0.2654
0.2179
湖北武汉宗关
0.21
0.1171
0.1533
0.2063
0.2422
0.2247
0.3221
江西南昌滁槎
0.1354
0.1055
0.0466
0.0352
0.054
0.0803
0.0676
江西九江蛤蟆石
0.2243
0.1985
0.2177
0.2331
0.2859
0.2328
0.2636
江苏扬州三江营
0.2253
0.0978
0.1965
0.1505
2.6103
0.3189
0.3655
点位名称
200401
200402
200403
200404
200405
200406
200407
四川攀枝花龙洞
0.309
0.2068
0.2928
0.2951
0.2373
0.2871
0.2763
重庆朱沱

0.333
0.2371
0.2768
0.2383
0.2537
0.2565
0.252
湖北宜昌南津关
0.2667
0.2735
0.3786
0.2481
0.3185
0.2598
0.2553
湖南岳阳城陵矶
0.2151
0.2101
0.2787
0.2774
0.3753
0.2641
0.2208
江西九江河西水厂
0.2736
0.2308
0.2814
0.2298
0.2301
0.237
0.2246
安徽安庆皖河口
0.2813
0.2853
0.2347
0.2468
0.2128
0.2084
0.1973
江苏南京林山
0.2874
3.2629
0.2792
0.2393
0.2113
0.2021
0.1537
四川乐山岷江大桥
0.0862
0.0634
0.0437
0.0822
0.0972
0.1519
0.1159
四川宜宾凉姜沟
0.105
0.368
0.1879
0.1517
0.2406
0.2866
0.2078
四川泸州沱江二桥
0.0364
0.0559
0.038
0.034
0.1077
0.2086
0.195
湖北丹江口胡家岭
0.2909
4.2573
4.8189
4.5177
4.8517
5.308
4.8931
湖南长沙新港
0.1463
0.1195
0.1123
0.1318
0.1501
0.1832
0.1988
湖南岳阳岳阳楼
0.2328
0.2496
0.141
0.2732
0.2872
0.177
0.1709
湖北武汉宗关
0.3231
0.2725
0.2804
0.2491
0.2615
0.1501
0.153
江西南昌滁槎
0.0624
0.0077
0.0508
0.135
0.121
0.0441
0.1209
江西九江蛤蟆石
0.2456
0.1935
0.2294
0.2203
0.1522
0.2237
0.2237
江苏扬州三江营
0.3722
0.3186
0.2981
0.2435
0.1736
0.2158
0.2435
点位名称
200408
200409
200410
200411
200412
200501
200502
四川攀枝花龙洞
0.3627
0.2049
5.1116
0.305
5.4133
5.3619
4.9318
重庆朱沱
0.2529
0.2198
0.2994
0.3227
0.2716
0.2898
0.2793
湖北宜昌南津关
0.2266
0.2468
0.2742
2.9366
0.2961
0.3004
0.3394
湖南岳阳城陵矶
0.2695
0.2794
2.6103
0.2478
0.2765
0.2179
0.315
江西九江河西水厂
0.2051
0.209
0.2462
0.2778
0.2886
0.3688
0.3656
安徽安庆皖河口
1.9578
0.2211
0.1979
0.2204
4.5461
0.2954
0.3085
江苏南京林山
0.2075
0.21
0.2006
0.2062
5.0743
0.3623
0.3427
四川乐山岷江大桥
0.1509
0.1504
0.2272
0.2354
0.1663
0.0836
0.0733
四川宜宾凉姜沟
0.271
0.2775
4.7074
0.3006
0.3043
0.369
0.333
四川泸州沱江二桥
4.264
4.2037
4.9668
0.2912
0.2533
0.2873
0.2713
湖北丹江口胡家岭
3.9402
0.2411
4.77
6.1453
0.4146
6.0432
63.0211
湖南长沙新港
0.1868
0.2364
0.2146
0.1804
0.1941
0.2663
2.3917
湖南岳阳岳阳楼
0.2187
0.2621
0.2318
0.2195
0.3071
0.2923
0.3196
湖北武汉宗关
0.2236
0.2243
0.2491
0.2863
0.2321
0.2404
0.3068
江西南昌滁槎
0.0637
0.1081
0.0798
0.064
0.0809
0.0626
0.0579
江西九江蛤蟆石
0.2227
0.226
0.1729
0.1588
0.2009
0.3655
0.3689
江苏扬州三江营
0.1744
0.1767
0.2654
0.2967
0.2819
0.3656
0.3753
点位名称
200503
200504
200505
200506
200507
200508
200509
四川攀枝花龙洞
0.3135
4.7474
0.2736
0.2873
0.2813
0.2915
0.2774
重庆朱沱
0.2715
0.2242
0.1943
0.3297
0.2664
0.2706
0.2677
湖北宜昌南津关
0.3134
0.2684
0.2573
0.2318
2.2841
0.1698
0.2125
湖南岳阳城陵矶
0.2793
0.2442
3.589
0.2793
0.284
0.2044
0.2923
江西九江河西水厂
5.2078
4.05
0.2422
0.2194
0.2344
0.2374
0.2244
安徽安庆皖河口
0.2892
0.2824
0.2171
0.2605
0.2324
0.2344
0.2325
江苏南京林山
5.522
0.225
0.1984
0.2221
0.2269
0.2012
0.1985
四川乐山岷江大桥
0.1373
0.0785
0.1621
0.206
0.1365
0.1665
0.2155
四川宜宾凉姜沟
0.29
0.2522
0.2328
0.2569
0.2332
0.2491
0.2674
四川泸州沱江二桥
0.2507
0.3156
0.2709
0.2579
0.2191
0.2188
0.2885
湖北丹江口胡家岭
6.4633
4.8435
4.1369
4.2448
4.0087
4.1921
4.3501
湖南长沙新港
0.308
0.1945
0.147
0.227
0.1684
0.1497
0.1585
湖南岳阳岳阳楼
0.2031
0.2192
0.1605
0.2487
0.2592
0.2526
0.2221
湖北武汉宗关
0.3036
0.271
0.1543
0.0953
0.141
0.2031
0.2181
江西南昌滁槎
0.1346
0.1617
0.1043
0.0958
0.0736
0.0842
0.0739
江西九江蛤蟆石
0.346
0.2729
0.2153
0.2064
0.2083
0.2351
0.2413
江苏扬州三江营
0.372
0.3623
0.1985
0.143
0.2041
0.1442
0.2366

4.2 模型二：流量模型
4.2.1模型二的建立
因长江源头的水可以考虑是由冰雪融化而得,故对长江干流起始点四川攀枝花按Ⅰ类水来确定是否是污染源，且取该Ⅰ类水主要污染物高锰酸盐指数和氨氮的浓度分别为2mg/L和0.15mg/L。考虑自然降解对河道污染的缓解作用，根据污染物浓度（与污染标准比较）来判断该地区是否为某种主污染物的主要污染源，且取降解系数取0.2　(单位：1/天)。将长江干流近一年多（04.04—05.04）主要污染物高锰酸盐指数和氨氮以及水流量和水流速综合于“附录表十八”。
各观测站点间水流速度取相邻两观测站点的平均速度，即为：
故各观测站点间水流时间为：
() 其中
由降解系数（取为0.2），运用衰减理论，可得干流各观测站点每秒排放CODMn()和NH3-N()的浓度为：

其中，取
若 则有
否则
运用“附录表十八”中的数据，用Matlab程序（见附录程序二）求解可得各观测站点的每秒排污污染物CODMn和NH3-N的质量，见表五：

表五：观测站点一年多每秒排污污染物的质量
观测站点
四川攀枝花
重庆朱沱
湖北宜昌
湖南岳阳
江西九江
安徽安庆
江苏南京
2004.04
CODMn
0.1107
3.4680
2.6946
6.6351
5.7673
2.5448
1.5268
NH3-N
0
0.2010
0.4515
0.4431
0.5711
0.1169
0
2004.05
CODMn
0.8556
0.5368
4.9567
5.5280
6.2888
3.7084
2.3159
NH3-N
0
0.3155
0.4733
0.4708
0.3935
0.3411
0.0063
2004.06
CODMn
0.2005
2.6931
5.9120
4.1129
3.5651
0
4.1756
NH3-N
0
0.2305
0.5821
0.4513
0.0320
0.1827
0
2004.07
CODMn
0.1864
3.5836
4.7356
6.3059
2.2256
0
3.1384
NH3-N
0
0.2647
0.2253
0.6768
0.0211
0.3134
0
2004.08
CODMn
1.4212
0
5.9174
6.4857
2.1132
0.4507
2.6715
NH3-N
0.3179
0
0.3266
0.6172
0.3356
0
0
2004.09
CODMn
2.575
3.256
5.505
10.604
0
0
9.520
NH3-N
0
0.501
1.090
0.729
0.675
0
0
2004.10
CODMn
0
2.0697
6.0056
3.5326
0.8394
1.9535
3.1520
NH3-N
0
0.3204
0.2867
0.5372
0.0729
0.0714
0
2004.11
CODMn
0.1200
0.8468
1.5621
2.5349
2.3811
1.4039
2.3066
NH3-N
0
0.2037
0.1934
0.3581
0.2062
0
0.0314
2004.12
CODMn
0
0.4909
1.6212
3.0860
2.1164
0.4871
2.3615
NH3-N
0
0.3160
0.0577
0.2878
0.0769
0.0978
0.0981
2005.01
CODMn
0
0.3992
0.7886
3.9334
1.8536
1.5861
2.8424
NH3-N
0
0.2099
0.0675
0.3158
0.1661
0.1050
0.2585
2005.02
CODMn
0
0.5854
0.8560
2.3461
2.7410
1.3800
2.6481
NH3-N
0
0.1876
0.0671
0.3245
0.0732
0.1494
0.4131
2005.03
CODMn
0
0.8102
1.0431
2.6942
1.6473
3.5094
2.4525
NH3-N
0.0069
0.2346
0.0857
0.2783
0.0624
0.3536
0
2005.04
CODMn
0
0.6392
1.2084
1.9195
2.5111
2.3492
2.1785
NH3-N
0
0.2217
0.0426
0.2943
0.1377
0.1422
0.3679
注：（单位：10kg/s）

4.22 模型二解答：
用图表方法将上表五的数据整理为图一和图二

图一：观测站点排放物CODMn随时间的变化

图二：观测站点排放物随时间的变化
干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要所在地即为其排放量的前几名。由于我们只研究干流七个地区的情况，污染主要所在地不妨取三个。由图一和图二可知：
长江干流近一年多主要污染物高锰酸盐指数的污染源主要在：湖南岳阳、湖北宜昌、江苏南京。
长江干流近一年多主要污染物氨氮的污染源主要在：湖北宜昌、湖南岳阳、重庆朱沱。
4.3 模型三：预测模型
4.3.1模型三的建立
由1995年~2004年长江流域水质报告，应分时分类地对河长进行预测，即如将1995年枯水期全流域评价河长、1996年枯水期全流域评价河长……2004年枯水期全流域评价河长等十个数据作处理，来预测2005年枯水期全流域评价河长、2006年枯水期全流域评价河长……2014年枯水期全流域评价河长等十个数据。
下面我们来分析预测未来长江水质的变化，通过对长江十年来的流量统计可知，除98 年外，长江每年的流量是一个比较稳定的值，大约在9 千亿立方米到