长江水质的评价与预测.doc

1、13 组 聂本武（建模） 张丰宇（写作）. 长江水质的评价与预测摘要本文讨论如何设计对长江水质污染情况进行综合评价，对各个地区水质污染状况分析，并判断出污染物高锰酸盐和氨氮的主要污染源，以及对未来水质情况进行预测的模型，然后根据预测的情况对长江未来的水质情况采取切实可行的治理方案，并提出合理的建议与意见。根据题目附件中已有的数据和搜集的一些综合评价和预测模型，并根据实际情况作了适当的假设，对不同要求的题目建立了不同模型并进行了较为完整的求解。对于问题一：题目要求对长江水质污染情况做出定量的综合评价。根据题目要求建立了模糊综合评价模型（模型一）来评价长江水质。本文首先对附件 3 中 2003.6

2、2005.9 这两年多来 17 个观测站 28 个月的水质数据进行处理，分别求出各个观测站水质处于各类污染的隶属度，建立单因子模糊评价矩阵，结合评价指标的权系数向量，求出反映 17 个观测站水质状况的模糊综合评价矩阵，并进行归一化处理。评价结果为：长江全流域 I 类水质断面占 17.65%，II 类水断面 47.06%，III 类水断面 23.53%，IV 类水断面5.88%， V 类水断面 5.88%，并得到各地区的水质情况。对于问题二：题目要求判断出污染物高锰酸盐和氨氮的主要污染源。根据题目要求建立了稳态一维对流扩散水质模型（模型二） 。本文首先利用附件 3 中给出的相关数据，求出长江干流

3、 6 个江段高锰酸盐和氨氮的污染量，再结合支流的地理位置及支流观测站的污染浓度数据，分析相关图像。最后得出长江干流近一年多主要污染物高锰酸盐和氨氮的污染源均主要分布在：湖北宜昌至湖南岳阳江段、重庆朱沱至湖北宜昌江段以及四川乐山地区。对于问题三：题目要求预测未来 10 年的水质情况。根据题目要求建立了 GM(1,1)模型（模型三） 。本文首先利用灰色系统理论对长江未来水质污染的发展趋势做出预测，然后用 19962004 年的模拟值、残差对报告表进行检验。经检验可知预测结果合理。最后得出结论为：可饮用水逐年下降，十年后将低于 55.1%。对于问题四：题目要求根据预测结果并将 IV、V、VI 类水的

4、比例控制在一定的比例内，求出每年需要处理的污水量。根据题目要求利用简单的比例关系求出每年需要处理的污水量。最后得出结论为：每年需要处理的污水量逐年递增，最小污水处理量为24.0897 亿吨，最大污水处理量为 68.8667 亿吨。对于问题五：题目要求提出解决长江水质污染切实可行的方案。本文通过对模型计算结果的分析，本文从国家、政府、企业和个人等方面提出了若干行之有效的建议。关键词：水质的评价与预测 模糊综合评价模型 稳态一维对流扩散水质模 GM(1,1)模型 一、 问题重述水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。专家们呼吁：“以人为本

5、，建设文明和谐社会，改善人与自然的环境，减少污染。 ” 长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。2004 年 10 月，由全国政协与中国发展研究院联合组成“保护长江万里行 ”考察团，从长江上游宜宾到下游上海，对沿线 21 个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出“若不及时拯救，长江生态 10 年内将濒临崩溃” （附件） ，并发出了“拿什么拯救癌变长江”的呼唤（附件 2） 。附件 3 给出了长江沿线 17 个观测站（地区）近两年多主要水质指标的检测数据，以及干流上个观测站近一年多的基

6、本数据（站点距离、水流量和水流速） 。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净化能力的指标称为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于 0.10.5 之间，比如可以考虑取 0.2 (单位：1/ 天)。附件 4 是“19952004 年长江流域水质报告”给出的主要统计数据。下面的附表是国标(GB3838-2002) 给出的地表水环境质量标准中 4

7、 个主要项目标准限值，其中、类为可饮用水。请你们研究下列问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如不采取更有效的治理措施，依照过去 10 年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来 10 年的情况。（4）根据你的预测分析，如果未来 10 年内每年都要求长江干流的类和类水的比例控制在 20%以内，且没有劣类水,那么每年需要处理多少污水？ （5）你对解决长江水质污染问题有什么切实可行的建议和意见。附表: 地表水环境质量标准 （G B3

8、8382002） 中 4 个主要项目标准限值（单位：mg/L）序号分 类标准值项 目类 类 类 类 类 劣类1 溶解氧(DO) 7.5（或饱和率90%）6 5 3 2 02 高锰酸盐指数(CODMn) 2 4 6 10 15 3 氨氮（NH3-N） 0.15 0.5 1.0 1.5 2.0 4 PH 值（无量纲） 6-9二、 模型假设（1）假设长江干流的自然净化能力近似均匀，在任何情况下都保持恒定 0.2；（2）所有的数据都来自科学采集；（3）入河排污口水量与水质变化稳定；（4）长江各观测河段的水流水质状态稳定；（5）四川攀枝花上流的污染忽略。三、 符号说明符号 符号说明Q水流量（m 3/s）

9、C污染物浓度（mg/L ）W污染指标权重u水流速（m/s ）x河道长度（ km）t时间（ s）k降解系数 （1/ 天）0.2k四、 问题分析水质的评价与预测的难点在于某观测点的水质状况一直处于动态变化之中，且在一定程度上受到上游水质的影响。对于问题一：由于水体污染指标之间的相容性(不存在传递性、绝对化) , 且水质评价本身是一种实践性、时空性、技术性均很强的多属性、多指标决策。本文针对水质评价的特点、要求以及求解过程中的难点, 采用了直接确定多指标权重和有限样本隶属度的模糊分析法 1。建模时需要找出影响水质的各主要因素，确定评价因子集、评价集、隶属函数，然后通过计算各因素的权重和隶属度，得到综

10、合隶属度，确定水质级别。对于问题二：要研究分析长江干流近一年多高锰酸盐指数和氨氮的主要污染源地区，可考虑分两步的策略：第一步：确定七个观测站之间受污染最严重的江段。污染物在河流中符合一级反应动力学，可以建立稳定一维均匀河流水质模型，分别求出长江干流上六个江段高锰酸盐和氨氮的污染量（包含了该江段支流所汇入的污染、该江段沿岸产生的污染、上游污染经自然净化后剩余的污染） ；第二步：对于上一步中求得的污染严重的江段，进一步缩小范围确定主要污染源地区。再结合支流、湖泊的地理位置及支流观测站的污染浓度数据，作出相关图象，进行分析。对于问题三：根据过去 10 年的主要统计数据，要对未来 10 年的水质污染作

11、出科学的预测。本文首先对附件 4 的数据进行分析，然后建立 GM（1,1）模型来对未来十年长江水质污染的发展趋势做出预测，再用灰色系统理论应用软件对数据进行处理得到未来十年的长江水质报告表，最后用 1996 年到 2004 年的模拟值、残差对报告表进行检验，由此可知得到的水质报告表是比较合理的。对于问题四：题目要求根据预测结果并将 IV、V、VI 类水的比例控制在一定的比例内，求出每年需要处理的污水量。根据题目要求利用简单的比例关系求出每年需要处理的污水量。对于问题五：根据以上四问的结果及附录一和附录二的信息提出解决长江污染问题的合理性建议。主要从国家、政府、企业和个人等方面提出可行的建议。五

12、、 模型的建立与求解5.1 问题一的求解水质污染程度是一个模糊概念，在进行评价时很难给出确切的表达，因此模糊综合评价模型来处理问题一，评价结果比较合理、更加接近客观实际。（1）建立模糊评判矩阵记模糊评判矩阵为 ，其中 表示在第 i 个观测站测得得处于第 j 级()ijmnRr()ijijrax污染程度的隶属度，隶属度是通过对隶属函数的计算来确定的，隶属函数一般采用“降半梯形”的函数。由于劣类污染物指标较极端，我们根据标准将水质分成 5 级。以溶解氧（DO）为例，即溶解氧应有对应于 5 个级别的隶属函数。以 DO 的监测值为自变量 x,对第 j 级别的隶属度为（水质标准见附表） 10,6()()

13、/1.,7.,75DOxx2()(.)/,6.50,.DOxx或3(3)/2,(),6DOx或4()(5)/2,30,DOxx或51()(),3DOx将各监测断面的监测数据代入前面确定的隶属函数中，就可以计算其隶属度，进而建立每个断面的单因子模糊评价矩阵。（2）确定评价指标的权系数向量当对多个（m 个）目标进行综合模糊评价时，还要对各个目标分别加权。权重是衡量因子集中某一因子对水质污染程度影响相对大小的量,权重系数越大,则该因子对水质的影响程度越大。设第 i 个目标权系数为 Wi，则可得权系数向量 A=(W1,W2,Wm)，满足。这里仅需要考虑 DO、CODMn，NH 3-N 这三种污染物对水

14、质的影响，1,0iiiW对它们赋予不同的权重 Wi。我们以污染物的超标情况决定权重，各因素的监测值相对于水质标准的超标越大，对污染的贡献越大，从而权重越大，可用下面的公式求权重系数：对于 CODMn 和 NH3-N 指标： ! iiCIS对于 DO 指标，因为 DO 与其它因素性质相反，实测 DO 浓度大，说明水质污染不严 2重，水质好。所以 DO 的权重赋值取 Ci/Si 的倒数，即：iiI其中，I i（i=1,2,3）表示第 i 种污染物的权重；C i 表示第 i 种污染物的浓度实测均值；Si 为第 i 种污染物的六个级别浓度标准限值的均值。为了进行模糊复合运算，还必须对各因子权重进行归一

15、化处理，即： 31iiiIW根据各观测站中三种污染因素的实测浓度 Ci 和平均浓度 Si 得到了因素权重分配集X 的归一化处理结果（见表 1） 。表 1 三种因素权重的归一化处理结果观测站 溶解氧 DO 高锰酸盐指数 CODMn 氨氮 NH3-N1 0.5034 0.3247 0.17192 0.4658 0.2490 0.28523 0.4623 0.3230 0.21474 0.3929 0.3760 0.23115 0.5563 0.2982 0.14556 0.5254 0.2891 0.18557 0.6072 0.2754 0.11748 0.3452 0.2919 0.36289

16、 0.3994 0.2818 0.318810 0.3555 0.2354 0.409111 0.5893 0.3058 0.104812 0.3535 0.1789 0.467613 0.3807 0.3723 0.247014 0.4948 0.3541 0.151015 0.1429 0.0569 0.800216 0.4288 0.3694 0.201817 0.4566 0.3216 0.2218（3）求模糊评价矩阵 B利用矩阵的模糊乘法得到模糊综合评判矩阵 B：121212 212(,)(,)nm nmnrrARWbrr 归一化处理后得： 121(,)nmiiibbB同理可以得出剩

17、余 16 个观测站水质数据经归一化处理后所得的模糊综合评判矩阵B2B17（见表 2） 。表 2 17 个观测站水质数据归一化处理结果分类观测站点I 类 II 类 III 类 IV 类 V 类 所属类别四川攀枝花（干流） 0.3432 0.3182 0.2214 0.1172 0 1重庆朱沱 （干流） 0.4495 0.2753 0.2753 0 0 1湖北宜昌（干流） 0.3706 0.3706 0.2589 0 0 2湖南岳阳（干流） 0.3382 0.3382 0.3236 0 0 2江西九江（干流） 0.5000 0.5000 0 0 0 2安徽安庆（干流） 0.3921 0.3921

18、0.2157 0 0 2江苏南京（干流） 0.4550 0.4550 0.0899 0 0 2四川乐山（支流） 0.2000 0.2000 0.2000 0.2000 0.2000 4四川宜宾（支流） 0.2276 0.2276 0.1816 0.1816 0.1816 2四川泸州（支流） 0.2000 0.2000 0.2000 0.2000 0.2000 2湖北丹江口（支流） 0.6341 0.3659 0 0 0 1湖南长沙（支流） 0.2500 0.2500 0.2500 0.2500 0 3湖南岳阳（支流） 0.2514 0.2514 0.2514 0.2459 0 3湖北武汉（支流

19、） 0.2757 0.2757 0.2757 0.1728 0 3江西南昌（支流） 0.0427 0.2393 0.2393 0.2393 0.2393 5江西九江（支流） 0.2686 0.2686 0.2314 0.2314 0 2江苏扬州（支流） 0.2810 0.2810 0.2810 0.1569 0 3（4）评价结果通过对以上数据处理，最终得出长江全流域水质（见表 3） 。表 3 长江全流域水质分类水质情况 1 类 2 类 3 类 4 类 5 类水质比例 17.65% 47.06% 23.53% 5.88% 5.88%将表 3 数据绘制为饼状图（见图 1）图 1 长江全流域水质5.

20、2 问题二的求解问题二要研究分析长江干流近一年多高锰酸盐指数和氨氮的主要污染源地区。由问题出发，污染源地区的确定分两步走：（1）确定七个观测站之间受污染最严重的江段I 模型二的建立建立稳态一维河流水质模型 1由于江河自身对污染物都有一定的自然净化能力，所以污染物在江河中迁移时，随着迁移距离的增大，浓度会不断下降。其次水流速度对污染物的浓度也有明显的影响。例如：水流速度很大时，许多污染物到下游时还不及降解，浓度变化缓慢。基于上述分析，可利用污染物质在河流中运动变化的基本模型“稳态一维对流扩散模型” 2：（1）212cctxxuEs其中：c 为整个断面的平均浓度；u 为流速;E 为离散系数;x 为

21、河道长度 ;t 为时间;s1、 s2 分别反映内部降解变化和外部影响。由于长江水流状态一般为稳定、均匀流动状态，与江长相比，水面宽度很小，可认为污染物浓度在横断面上分布比较均匀，即对同一地点进行水质采样在无外界环境的影响下，可认为污染物浓度不随时间变化。其次，因为长江水流状态稳定、均匀，可认为污染物浓度随迁移距离线性变化，且不受外部影响。s 1 与江水的自然降解和污染物浓度有关。所以一维对流扩散的水质模型可简化得到这里的稳定一维河流水质模型：（2） cxduk其中：x 为河段长度； u 为河段平均流速； c 为污染物浓度； k 为降解系数。对（2）式分离变量积分得：（3）0xkuce其中：c

22、0 为初始位置的污染物浓度。利用（3）式，建立反映六江段污染物总量的模型 2考查干流上相邻两点间的污染情况，即在考虑自然降解的情况下这段流域上支流和干流的污染状况。干流上各江段污染物总量可以表达为：（4）(1)(,(1)()iixkuiiiiLQCe其中： 表示干流上第 i+1 个观察站（地区）的污染量；Q 表示水流量；C(1)i表示污染物浓度； 表示干流第 i 个观察站（地区）的污染经过自然降解(1)(,()iixkuie到下一个观测站（地区）时的剩余污染量。II 模型二的求解由于所给数据为 13 个月的各观测点（干流上）的离散数据，考虑用 13 个月的污染物总合来表征各江段的污染状况。利用

23、（4）式分别对高锰酸盐指数和氨氮进行计算，计算结果见表 4。表 4 长江六个江段的污染总量（单位：万吨）江段 1-2 2-3 3-4 4-5 5-6 6-7CODMn 109.5481 118.8481 155.0080 100.6016 56.5220 105.8197江段 1-2 2-3 3-4 4-5 5-6 6-7NH3-N 9.5109 10.7537 15.0006 8.6236 4.0140 0.5338由上表的数据可知：对于高锰酸盐（CODMn）和氨氮（NH 3-N）两种污染，湖北宜昌南津关至湖南岳阳城陵矶江段的污染最严重，其次为重庆朱沱至湖北宜昌南津关，和四川攀枝花重庆朱沱。

24、（2）缩小范围确定污染源上文已经确定了长江上高锰酸盐（CODMn）和氨氮（ NH3-N）两种污染最严重的江段。为了确定具体的污染源地区，我们采用作图法：两湖区进出口处污染比较 1CODMn 浓度比较（见图 2）：图 2 两湖进出口 CODMn 污染浓度比较结论：a 洞庭湖出口处 CODMn 浓度增大，从类水变成了类水。但比照地表水环境质量标准 （GB3838-2002）可知，其仍属于可饮用水范围。可见洞庭湖湖畔工厂对洞庭湖水质中的 CODMn 指标有清度污染。由上文结果可分析得知 污染最严重的湖北宜昌南津关至湖南岳阳城陵矶江段的 CODMn 主要污染源为该江段沿岸地区，而不是洞庭湖。b 鄱阳湖

25、出入口的水质按照国家质量标准，均属于类，为可饮用水。即鄱阳湖湖畔基本无 CODMn 污染。NH3-N 浓度比较（见图 3）：图 3 两湖进出口 NH3-N 污染浓度比较结论：a 洞庭湖出口处 NH3-N 浓度减小，从类水变成了类水。但比照地表水环境质量标准 （GB3838-2002）可知，其仍属于可饮用水范围。可见洞庭湖湖畔工厂对洞庭湖水质基本没有 NH3-N 污染，湖水湖畔入口流域的 NH3-N 污染起到了净化作用。由上文结果可分析得知污染最严重的湖北宜昌南津关至湖南岳阳城陵矶江段的 NH3-N 主要污染源仍为该江段沿岸地区，而不是洞庭湖。b 鄱阳湖对上游 NH3-N 污染的净化作用更加明显

26、。分析可知鄱阳湖入口江西南昌滁槎河段为 NH3-N 的主要污染源。岷江流域两观测点污染比较 2CODMn 浓度比较(见图 4)：图 4 四川乐山与四川宜宾 CODMn 污染浓度比较结论：近一年中，四川乐山的 CODMn 污染浓度普遍高于四川宜宾，其中第 1、11 月属于污染状态。NH3-N 浓度比较：注：“交汇前”指四川乐山观测站， “交汇后”指四川宜宾观测站图 5 四川乐山与四川宜宾 NH3-N 污染浓度比较结论：近一年中，四川乐山的 NH3-N 污染浓度普遍高于四川宜宾，且第 1、4、10、11 月这四个月份均属于污染状态。综上所述：长江干流近一年多高锰酸盐（CODMn）和氨氮（ NH3-

27、N）的主要污染源是：湖北宜昌至湖南岳阳江段沿岸、岷江流域的四川乐山地区、重庆朱沱至湖北宜昌江段沿岸。5.3 问题三的求解题目要求对未来长江的水质情况进行预测。由于只有过去十年的统计数据，信息不完整，所以用对信息质量要求不高的灰色系统分析法 3来作预测，建立 GM（1,1）模型。由于 1998 年长江流域遭遇罕见的洪水灾害，数据变化很大，因此我们先将长江 1998 年的数据作为灾变数据剔除。记 ，其中 表示第 i 年长江流域水文(1),2,()xxn ()xi年可饮用水（、类）的河长比例。（1）模型的建立令 为 GM（1,1）建模序列，表示灰色导数 1 (0)x(0)()(0)(0)1,2,xx

28、n其中： (0),3xkk令 为 的 AGO 序列 2 (1)()(1)()(1)(1)0(1)()1,2,kmxxn令 为 的均值（MEAN）序列，表示白化背景值 3 (1)z()x(1)(1)(1)()0.5.2,3,zkxkzn则得到 GM（1,1）的灰微分方程模型为 (0)(1)xkab其中： (1)(0) (1)(0)222(1)()(1)()2(1)(1)(0)2 2()()2nnnkkknnnnkkkkzzxkazzzxkbz经变换得到目标函数： (0)(1)xkbazk约束条件为： (1)(0) (1)(0)222(1)()(1)()2(1)(1)(0)2 2()()2.nnn

29、kkknnnnkkkkzxzxkazstzzxkbz（2）模型的求解利用 MATLAB 编程求出当 k=11,1220 时的 1 (0)xk因为 ，变换后得到 2 (0)()xk(0)2,1,3)xkk最后预测出长江流域未来十年水文年可饮用水河长比例（见图 6） 。图 6 未来十年可饮用水河长比例预测值进行残差检验，结果见表 5：表 5 长江全流域可饮用水的河长比例实测值与预测值比较（单位：%）年份 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004预测 93.1 83.31 81.42 79.57 77.77 76.00 74.27 72.59

30、 70.94 69.33实测 93.1 85.3 80.7 88.4 80.2 74 73.7 76.7 77 68残差 0 -1.97 0.72 -8.82 -2.43 2.00 0.57 -4.11 -6.56 1.33相对误差 0 -2.33 0.89 9.99 3.03 2.70 0.78 5.36 -8.46 1.95从上表可以看出，除 1998 年因为长江发生罕见自然灾害，预测值误差较大外。其余相对误差都小于 9%，可以认为使用 GM（1,1）模型预测效果是令人满意的。将长江未来 10 年水质的预测结果用表 6 表示：表 6 对长江流域未来十年可饮用水河长比例的预测（单位：%）年份

31、 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014比例 67.55 66.22 64.71 63.24 61.81 60.40 59.03 57.69 56.38 55.10由上表分析得出：未来十年内，长江流域可饮用水的河长比例逐年下降。到 2014 年，全流域可饮用水的河长比例下降至 55.10%，即近一半为非饮用水河段。污染如此严重，长江生态濒临崩溃。5.4 问题四的求解根据问题三的预测，长江水体的污染将日益严重，急待解决的是如何使未来 10 年内每年的干流、类水的比例控制在 20%以内，且没有劣类水。利用解决问题三的方法，基于过去 10

32、年的统计数据，预测未来十年（2005 i2014）内的排污量。令每年污水年均排放量为 Pi； 1令每年全流域非饮用水类的河长比例为 qi； 2令每年干流、类水和劣类水的河长比例分别为 ri、t i； 3由于污水是非饮用水类，即每年的污水主要分布在全流域、类和劣类水的河段上。在假定干流和支流的污染河段上，单位长度河段的污水含量大致相同的前提下，利用简单的比例关系可粗略地预测未来 10 年内每年的污水处理量。基于过去 10 年统计数据，分别以干流河长和全流域河长的均值作为未来 10 年的评价河长（L 干 ：干流评价河长；L 全 ：全流域评价河长） ，其数学计算式可描述如下：(1)(2)()()(0

33、5214)i iii iiiiiPLqLqi全 全 全其中：Li 表示第年的污水处理量；L(1)i 表示第年需要治理的、类河长，L (1)i=L 干 *(ri-20%)，r i 表示第 i 年、类河长的预测比例；L(2)i 表示第年需要治理的类河长，L (2)i=L 干 *ti，t i 表示第 i 年类河长的预测比例；Pi 表示第 i 年污水年均排放量的预测值；qi 表示第 i 年全流域非饮用水类的预测比例。由此解得每年需要处理的污水量见表 7：表 7 每年需要处理的污水量（单位：亿吨）年份 2005 2006 2007 2008 2009污水年处理量 24.0897 28.6064 33.0

34、651 37.6167 42.3310年份 2010 2011 2012 2013 2014污水年处理量 47.2051 52.2562 57.5597 63.0846 68.8667从上表可以看出 2005 年到 2014 年需要处理的污水量呈现逐年上升趋势，与污水排放量呈现逐年上升趋势一致，故认为，求得结果比较合理。5.5 问题五的求解从问题一到问题四的求解过程可知 1995 年到 2004 年废水排放量逐年增加，并且经过预测可知 2005 年到 2014 年废水排放量也在逐年增加。由数据分析可知，年废水排放量增加速度有增无减，按此趋势发展下去，将会导致长江无饮用水可用。综上分析可知：近年

35、来长江水质情况已经严重恶化，对人类的生活产生了严重威胁，由此可知解决长江水污染问题迫在眉睫。有鉴于此，本文就国家、政府、企业和个人等方面提出若干可行的建议和对策：（1）加快国家保护长江的立法进程。一方面制定促进循环经济发展的政策和法律法规；另一方面，还需要将建议加大对违法排污行为的处罚力度，要真正的做到“谁污染，谁治理” ，把对长江的生态环境的合理开发纳入沿江城市政府官员的政绩考核体系，对严重破坏生态环境并造成生态恶果的地区，应执行官员任用的一票否决制。（2）强化政府污水治理和监督管理的职能。政府要对每个污染源必须进行有效的治理,各排污单位要根据其废水的特点,选择相应的污水处理工艺和设备对污水

36、进行治理,达到排放标准后才能排放。关闭或责令停业整改污染排放已经超标的造纸、皮革、钢铁、化工等重污染企业。对未治理或治理后未达标的工业废水一律不准向湖泊水域直接排放。（3）加大企业处理工业废水的资金投入。随着中国社会的快速发展和城市化进程的加快，大型企业的污水排放量大幅增加，而企业的污水处理能力明显不足。除了污水处理能力有限外，由于资金投入不足，导致企业现有的污水处理工艺达不到国家新的标准要求。为解决财政能力有限的难题，取得社会效益和经济效益 “双赢”的效果，企业要加大污水处理项目建设资金投入（4）提高全民的环保意识。努力唤醒并提高全民的环保意识，大力普及环保教育，加强环保宣传。使“环境保护，

37、人人有责”的观念深入人心。六、 模型的评价与优化6.1 模型的优点（1）对于问题一建立的模糊评价模型能够合理的反映出各地区的水质状况，同时利用综合指标定量的反映出地区水质的优劣状况和长江干流的总体水质状况。（2）对于问题二建立的模型找到的污染源区域与事实的情况比较符合，说明我们建立的模型较合理。（3）问题三建立的灰色预测模型利用有限的数据预测出未来年份的长江各污染物的发展趋势，具有利用少量数据预测有限的信息的优点。（4）我们所建立的模型对题目中所给的数据进行了合理的处理，充分地利用了数据源提供的信息，基本上解决了“长江水质的评价和预测”问题。（5）本文在研究各种污染产生的影响时，将数据、表格和

38、图示相结合，使结果一目了然。6.2 模型的缺点（1）模型二中忽略了各污染物内部反应与相互作用项，如生物化学中的生长与降解变化。（2）由于长江河道的复杂性和周围环境的影响使得预测存在着一定的误差。6.3 模型的优化可以针对实际情况对模型做进一步讨论，首先对江河自净系数，在问题(2)中所建立的模型中，我们认为长江干流的自净系数的均匀且相等的，实际情况下自净系数是变化的而且可通过经典的斯特里特菲尔普斯(Streeter-Phelps)公式对不同江段的自净系数求解，这样可以在对长江水质测量分析等应用中提高测量精度；另外还可以增加对长江流域支流的分析考虑，一般长支流的流量对长江干流影响不大，可实际上长江

39、有些支流却是比较重要的污染源；最后可以对长江水质预测的进一步考虑，在处理问题(3)的过程中，我们曾采用灰色预测的方法对过去 10 年的数据进行处理时发现结果的杂乱性。实际应用中，我们可以通过积累更多，更有效的历史数据或者研究更科学，更有效的预测算法保证预测的精度。七、 模型的推广通过对题目的解读不难发现这是一类综合评价和对未来情况预测的问题。本文建立的是一个对长江水质的评价和预测的模型，仔细分析建立的模型，不难发现：这个模型不仅适用于各类水的水质评价与预测问题，它对于与之相关的综合评价和预测的问题都可以起到指导作用。本题的求解是一个典型的综合评价与预测的问题，模型的使用范围非常广泛，在各种水的

40、水质评价和预测，各种商品未来销售情况的预测，学生综合成绩的评价等领域都有着积极的指导意义。总之，凡是涉及到对事物的综合评价和预测的问题，都可以用本模型进行解决。参考文献1谢季坚,刘承平 .模糊数学方法及其应用.武汉:华中科技大学出版社,2000.2徐祖信等. 基于数学模型的苏州河上游和支流水质对干流的影响分析.水动力学研究与发展,A 辑第 19 卷第 6 期,733_743,2004.3邓聚龙.灰预测与灰决策.武汉:华中科技大学出版社,2000.4郑阿奇.MATLAB 实用教程 .北京:电子工业出版社,2004.5宋来忠,王志明 .数学建模与实验.北京:科学出版社,2005.附录附录一：（见题

41、目）附录二：（见题目）附录三：（见题目）附录四：（见题目）附录五：将题目所给附件 2 中长江流域主要城市水质检测报告中的数据以矩阵形式保存到MATLAB 中，数据名为 d,文件名为 data.其中矩阵行数对应题目附件中数据的行数，1-6列分别代表 PH、DO、CODMn、NH3-N 浓度和本月、上月水质等级。附录六：建立单因子模糊评判矩阵load datam,n=size(d);do=zeros(17,5);mn=zeros(17,5);nh3=zeros(17,5);for j=1:17,for i=1:28,x1=d(i-1)*17+j,2);x2=d(i-1)*17+j,3);x3=d(

42、i-1)*17+j,4); if x17.5,do(j,1)=do(j,1)+1;elseif x1=6endif x30.15endendend附录七：确定评价指标的权系数向量d=xlsread(d1.xls,sheet1,A1: G476)z=zeros(17,3);for i=1:17,for k=1:28,z(i,:)=z(i,:)+d(k-1)*17+i,2:4);endendx=z./28;s(1)=mean(7.5 6 5 3 2);s(2)=mean(2 4 6 10 15);s(3)=mean(0.15 0.5 1.0 1.5 2);a=zeros(17,3)a(:,1)=s

43、(1)./x(:,1);a(:,2)=x(:,2)./s(2);a(:,3)=x(:,3)./s(3);for j=1:3for i=1:17w(i,j)=a(i,j)/sum(a(i,:)end end w(i,j)附录八：求综合模糊评判矩阵load datafor i=1:17,R=do(i,:);mn(i,:);nh3(i,:);A=w(i,:);for j=1:5,for k=1:3,tj(k)=min(A(k),R(k,j);endB(i,j)=max(tj);endendfor i=1:17,for j=1:5,C(i,j)=B(i,j)/sum(B(i,:);endenddisp

44、(综合模糊评判矩阵:)C附录九：分别求出长江干流各江段污染物总量load datav=3.7 2.1 0.9 0.9 1.0 1.1 1.23.7 1.9 0.8 0.9 1.1 1.1 1.23.9 2.1 1.2 1.3 1.5 1.5 1.64.1 2.3 1.4 1.5 1.5 1.6 1.73.8 2.1 1.4 1.4 1.5 1.7 1.75.1 4.8 1.7 1.9 2.1 3.4 3.43.1 2.3 1.5 1.6 1.6 1.7 1.92.7 1.9 0.7 0.8 0.9 0.9 1.03.1 1.5 0.7 0.8 0.8 0.8 0.92.1 1.5 0.5 0.

45、6 0.7 0.7 0.82.0 1.0 0.4 0.6 0.7 0.7 0.81.9 0.9 0.4 0.6 0.8 0.8 0.92.1 1.2 0.4 0.5 0.7 0.8 0.8;for i=1:6,for j=1:13,jun_v(i,j)=(v(j,i+1)+v(j,i)/2;endendjun_v=jun_v;s=950 (1728-950) (2123-1728) (2623-2123) (2787-2623) (3251-2787);s=s*1000for i=1:13,t(i,:)=s./(jun_v(i,:);endt=t./3600./24;k=3600*24*30;

46、L_codmn=zeros(1,7);L_nh3=zeros(1,7);a=0.2Q=3690 13800 21000 25600 28100 29500 298003720 13100 19800 20500 29800 34000 345004010 14200 20300 22600 29500 32100 331004660 16400 22700 24100 27000 31900 321003740 10600 24000 25900 32100 33400 351006280 47600 53500 53800 72800 74200 810003260 16200 19100

47、22300 24800 31000 384001500 8170 10600 12000 14600 17000 19600951 6550 7400 10700 13200 14100 14900712 4020 4570 8190 10900 12300 14400612 3603 4510 7980 10300 13700 15100623 4740 5180 7040 14300 21400 21500642 3650 5400 7240 15100 20200 22100;L_codmn(1)=sum(codmn(:,1).*Q(:,1)*k;L_nh3(1)=sum(nh3(:,1

48、).*Q(:,1)*k;for i=2:7,L_codmn(i)=k*sum(codmn(:,i).*Q(:,i)-codmn(:,i-1).*Q(:,i-1).*exp(-a.*t(:,i-1);L_nh3(i)=k*sum(nh3(:,i).*Q(:,i)-nh3(:,i-1).*Q(:,i-1).*exp(-a.*t(:,i-1);endL_codmn=L_codmn/1000/1000/10000L_nh3=L_nh3/1000/1000/10000disp(各江段分别排放CODMn总量（单位：万吨） )L_codmndisp(各江段分别排放NH3-N总量（单位：万吨） )L_nh3附录十：利用灰度预测对长江全流域未来10年可饮用水河长比例预测n=9;g0=25.8+42.6+24.7,15.3+20.2+49.8,12.2+24.9+43.6,5.2+39.8+35.2,5.6+32.8+35.6,5.9+33.1+34.7,4.4+44+28.3,4.7+41.5+31.3,1.2+26.9+39.9;x0=sqrt(g0);x1=zeros(1,n);x1(1)=x0(1);for i=2:n,x1(i)=x1(i-1)+x0(i);endfor i=1:n-1,