长江水质评价及预测模型的建立与分析.pdf

《长江水质评价及预测模型的建立与分析.pdf》由会员分享，可在线阅读，更多相关《长江水质评价及预测模型的建立与分析.pdf（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 1 长江水质评价及预测模型的建立与分析 摘要：水资源污染日趋严重，本文仅就长江流域若干观测点的最近 10 年的观测数据进行建模分析，探讨水资源污染的预测与控制的问题。对于问题一、我们首先采用模糊综合评判法进行评价，我们以溶解氧、高锰酸盐指数、氨氮组成的评价参数集为 U=u 1,u 2,u 3 ,水质等级集为 v 1,v 2,v 3,v 4,v 5,v 6 ，用隶属度来刻画长江沿岸 1 7个测点水质在三个参数方面属于不同水质级别的情况，由观测值和隶属度确定模糊评判矩阵，由最大隶属原则，求出 2

2、8 个时点平均值得到 1 7 个测点的 2 年多的平均水质情况，其中江西南昌滁槎污染最为严重（测点 1 5），其次是四川乐山岷江大桥（测点 8）。我们还采用了另外一种一般的综合指数评价方法对水质情况作综合评价。选取等水的 4 项指标值 P H0=7，D O0=5，C O D0=6，N H0=1 为评价标准。先计算分指数 K，然后再把分指数加权（）求和得到综合指标 W=Kj。因为 1 7 个测点分别有 2 8 个时点的监测数据，我们把这 2 8 个月按年分为枯水期、丰水期和平水期。分类标准参照题中附件 4的说明（3）。然后用平均的方法求出按年按期分类的综合指标矩阵。用 E x c e l 分别画

3、出2 0 0 3、2 0 0 4、2 0 0 5 年各测点的水质综合评价柱状图（见附件中图 1、图 2、图 3）。其中等水的量化指标为 0.1，值越低水质越好。最后得出：江西南昌滁槎污染最为严重（测点 1 5），其次是四川乐山岷江大桥（测点 8）。可以看出两种方法评价的结果是一致的。对于问题二、考虑到江水的流动性和自降解性，测点浓度不能与污染源等同起来，所以我们把每个监测值看作是由上游对它的影响和本地污染两部分组成的。所以我们考虑水流速度、测点间距计算出降解次数，进而计算出一个测点对下一个测点的影响度，通过减去来自上游的影响进而求出每个测点的本地污染度。这样就可以通过比较本地污染度找出主要污染

4、源了。具体计算出主要污染源有：江西南昌滁槎、四川泸州沱江二桥、四川乐山岷江大桥、湖南长沙新港、湖南岳阳城陵矶、重庆朱沱、湖南岳阳楼等地区。对于问题三、问题要求预测未来十年的水质污染趋势，我们用预测未来 10 年中六 种水的比例来评价水质污染趋势。具体用灰色系统方法进行预测。又由于每年的总评价河长不一样，所以不能直接预测，但是比例是不会变的，我们把 10 年的六种水的比例用同一个总评价河长（我们假定为 10000）转化。再进行预测，得出未来 10 年的六种水的河长，然后再转化为比例值。可以看出未来 10 年长江劣类水越来多，污染越来越严重。对于问题四、我们把 6 种水重新分为 3 类：清洁水（I

5、、和）、污染水（和）、劣质水（）。用灰色系统方法预测未来 1 0 年原 6 种水的比例，折算出新 3 类水的比例，建立废水率与 3类水比例的逐步回归模型，确定未来 1 0年每一年废水率与各种水比例之间的回归方程。然后根据回归方程得到在控制状态下（和类不超过 2 0%、劣类没有）的废水率，然后与不加控制的预测值进行比较，得到需要处理的废水量（亿吨）。2005 年 2006 年 2007 年 2008 年 2009 年 2010 年 2011 年 2012 年 2013 年 2014 年 79.4 97.5 118.9 141.9 166.1 192.6 231.3 261.3 294.8 330

6、.1 对于问题五、根据我们所做模型的分析和预测结果，以及目前我国治理长江水质污 染的现状给出：合理制定工业发展规划、污水处理设施的利用、控制排污总量、加大执法力度等四点建议和意见。关键字：水质评价 模糊综合评价 综合指数评价 灰色系统 逐步回归 2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 2 正文：一、问题的提出“黄河脏了，长江可以救黄河，长江完了，拿什么救长江？”。我们的生命线长江正在倍受煎熬，保护长江、保护水资源就是保护我们自己。而如今长江污染日趋严重。那么长江现在到底现况如何？污染源在哪里？如果再不采取有效措施，长江的命运将会如何

7、？最主要的污染物污水的排放对长江产生什么样的影响？怎么样规划、采取怎么样的措施才能使长江在保持“生命力”的前提下达到环境与经济和谐发展？等等这些都是我们亟待解决的问题。二、问题的分析 2.1 对水质情况的定量综合评价 水质的评价是对水中各种物质的浓度进行测定，分析是否超过某个标准的过程。题目中只给出了 4 种物质作为评价指标。通过查阅资料我们知道一般分析水质情况的方法有指数评价法、模糊评价法、灰色评价法、物元分析法、人工神经网络评价法等诸多方法1。我们在这里先用模糊评价法作出评价，然后用比较简单常用的指数评价法进行验证。2.2 寻找主要污染源 污染源是指污染的源头，即污染物的排放点。但是题目只

8、给出了各个测点的污染物的浓度值，我们很容易理解一个测点的监测值高不一定说明该点就是主要污染源。因为江水是流动的，一个测点的数据应该包括两个部分，即来自上游的污染物的影响和本地的污染物的浓度。各测点本地污染物的浓度才是我们考察的目标。那么上游污染物对下测点的影响度该怎么计算呢？这就要考虑两测点间距离、水流速度及污染物降解与上游污染物影响度之间的关系。2.3 水质污染趋势的预测 预测需要有历史数据支持，根据题目所给数据的特点：过去 10 年每一年中 6 类水的比例，以及水流量总量及排污量等不确定因素。我们可以灰色系统方法进行预测。且用一年中 6 类水的比例也可充分表明水质状况。2.4 水质污染的控

9、制 在 2.3 的预测中，我们知道 6 类水的比例反映水质污染的程度，且题目也是要求我们通过对污水的处理达到对 6 类水的比例进行控制的目的，即 6 类水的比例是我们的控制目标，污水量是我们的控制点。所以我们要先找到污水量与 6 类水比例之间的关系，我们采用多元回归对污水排放比例和 6 类水的比例进行逐步回归分析，这样我们可以得到一个回归方程，它就是我们要找排污比例与 6 类水比例之间的关系式，通过它我们就可以在给定控制标准下的污水排放上限，然后跟我们预测的不加控制的排污量比较就得到了需要处理的污水量。三、模型的假设（1）假设溶解氧(D O)浓度越高水质越好，不考虑过含氧情况。（2）假设各监测

10、指标之间无相互作用。（3）假设我们研究的长江是一条平直的河流。（4）假设所给数据真实可靠。（5）假设水质状况只与题目给我们的 4 个项目有关，不考虑其他项目 2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 3 四、号的定义与说明 4.1 模糊评判模型符号定义 符号 符号说明 Ai j 第 i 个参数在第 j 级别上的标准值 Xi,k 第 i 个参数在某一级别上的监测值的第 k 个 li,j i 参数监测值介于 Ai,j-1到 Ai,j之间的个数 Li 第 i 个参数污染物监测值的个数 Zi,j i参数污染物监测值介于 Ai,j-1到 Ai,

11、j之间的 li,j个监测值的平均值 ri,j,j zi,j对第 j 级水质的隶属度 wi 计算得到的 i 参数的权重 pi,j i 参数而言，介于 Ai,j-1到 Ai,j之间的监测值发生在j 水质下的概率 ai 归一化处理后 i 参数的权重 qi,j i 参数发在 j 水质下的模糊概率 qj 水域水体出现 j 级水质的模糊综合概率 4.2 综合指数评价模型符号定义 符号 符号说明 W 水质综合指标 K 分指数（监测数据与评价标准之比）评价指标的权重 Cj 全流域 17 个测点 COD 的监测值 Nj 全流域 17 个测点 NH3 的监测值 Pj 全流域 17 个测点 PH 的监测值 Dj 全

12、流域 17 个测点 DO 的监测值 Czj 支流 10 个测点 COD 的监测值 Nzj 支流 10 个测点 NH3 的监测值 Cgj 干流 7 个测点 COD 的监测值 Ngj 干流 7 个测点 NH3 的监测值 五、模型的建立与求解 5.1 长江水质的综合评价 5.1.1 模糊综合评判模型 2 根据水域情况的质量标准我们把水污染监测浓度看成是一个离散的随机变量，用概率统计方法进行统计可以得到水域属于某个标准的概率，因为可以拟定不同的水域标准，评价参数集为 U=u1,u2,u3,水质分级集为v1,v2,v3,v4,v5,v6,其中 u1,u2,u3 分别表示为溶解氧，高锰酸盐指数，氨氮（N

13、H 3-N），因为 P H值对水域影响不大，所以对其不2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 4 予考虑，v1-v6 分别表示为类到劣类，设 i 参数污染物监测值共有 Li个，其中介于 Ai,j-1到 Ai,j之间的监测值有 li,j个，高锰酸盐指数，氨氮（N H 3-N），的监测值为(i=2,3)li,1 监测值 Xi,k Ai,1 li,j=li,j Ai,j-1Xi,kAi,j(j=2,3,4,5,6)li,n+1 Xi,k Ai,n(k=1,2,Li)溶解氧的监测值为(i=1)：li,1 监测值 Xi,kAi,1 li,j=

14、li,j-1 Ai,j-1Xi,kAi,j(j=2,n)li,n Xi,kAi,n 那么，对于 i 参数而言，介于 Ai,j-1到 Ai,j之间的监测值 Xi,k发生在 Ai,j-1下或Ai,j下 的概率为 pi,j=li,j/Li 我们用隶属度来刻画水质分级界限：已知水质等级标准为Ai,j,i参数污染物监测值Xi,k介于Ai,j-1到Ai,j之间的li,j个监测值的平均值为 zi,j=(Xi,k)/li,j(k=1,2,li,j),Xi,k为第 i个参数污染物位于 Ai,j-1到 Ai,j之间的监测值的第 k 个。则对 i 参数而言，zi,j对第 j 级水质和j-1 级水质的隶属度 ri,j

15、,j和 ri,j,j-1,可由下列各式求得(其中 i=2,3)，即：ri,1,1=1 zi,1Ai,j ri,j,j-1=Ai,j-zi,jAi,j-Ai,j-1 Ai,j-1zi,jAi,j ri,n+1,n=0 ri,j,j=zi,j一 Ai,j/Ai,j一 Ai,j-1 Ai,j-1zi,jAi,j (zi,n+1=0)ri,j,j-1=0 zi,j=0 ri,n+1,n=1 ri,1,1=0 zi,1=0 ri,j,j=0 (zi,n+1 Ai,n 而 i=1 时对于溶解氧的隶属度的求法与上面方法相反 对于评价参数的权重的确定：对于溶解氧权重按如下确定 w1=(x0-x1)/(x0-s

16、1)，而高锰酸盐指数，氨氮的权重分别为 wi=xi/si，其中 xi-第 i种污染物的实测浓度算术平均值，x 0-溶解氧在某条件下的饱和浓度（标准浓度），si-第 i种污染物各级标准的算术平均值。对其进行归一化处理得到 ai=wi/wi（i=1,2,3）3 个参数构成权重矩阵即为 A=（a1,a2,a3）三个指标的权重（见表 1）：2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 5 表 1 指标 溶解氧 高锰酸盐指数氨氮 a 0.2523 0.6854 0.053 i 参数发生各级水质的模糊概率为qi,j=ri,j,j*pi,j+ri,j+

17、1,j*pi,j+1(i=1,2,3,j=1,2,3,4,5,6),则模糊概率关系矩阵为q36,于是每个单指标模糊概率评价矩阵为 q16=A*q36（见附表1）,其中qj=ai*qi,j，即为水域水体出现 j 级水质的模糊综合概率，并按最大隶属原则求得qj判断属于第几个级别。在我们已知条件中我们对 28 个月中 17 个长江流域主要城市的水质监测项目求出平均值按以上说明做出模糊综合评判，由 matlab 程序求出 17 个城市 28 个月的平均水质情况，得到长江沿岸17 个观测点近两年多的水质平均情况（见表 2）：表 2 测点 1 2 3 4 5 6 7 8 9 1 0 1 1 1 2 1 3

18、 1 4 1 5 1 6 1 7 水质等级 IV 劣 5.1.2 综合指数模型3 为了验证我们评价的客观性，我们现在用综合指数评价模型对模糊综合评价模型进行验证。首先求分指数 K 4 个指标的 K 分别为：Pj/P H0 ；Dj/D O0 ；Cj/C O D0；Nj/N H0 其中选取第三类等级水的 4项标准 P H0=7，D O0=5，C O D0=6，N H0=1 为评价标准。由于 D O 是逆指标，P H 是非单调指标（不是值越大表示水质越好或越坏）。所以综合指标计算公式如下：W=10PHPj+jDDO0+0CODCj+0NHNj （不加权的）求权数 权重的确定是参照附件 3 中的数据，

19、通过统计 2 8 组数据中“主要污染指标中”各指标出现次数，依据次数越多，权重越大的思想计算各权重（见表 3）表 3 P H D O0 C O D0 N H0 1.1%2 5.9%1 6.5%5 6.5%W=110PHPj+2jDDO0+30CODCj+40NHNj （加权）因为有 28 个月，这样我们就得到了一个 1728 的综合指标矩阵，数据过于庞大，2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 6 不便直观分析，所以我们按照原题附件 4 中说明（3）的划分标准把 28 个月划分为 2003 年丰水期、2003 年平水期；2004 年

20、枯水期、2004 年丰水期、2004 年平水期；2005年枯水期、2005 年丰水期、2005 年平水期。每个时期按所包含月份求平均，这样就大大缩小了数据量。用直观的柱形图表示如下图（按年评价图见附件图1图3）：水质综合指数评价直观图00.0 50.10.1 50.20.2 50.3四川攀枝花龙洞重庆朱沱?湖北宜昌南津关?湖南岳阳城陵矶?江西九江河西水厂?安徽安庆皖河口?江苏南京林山?四川乐山岷江大桥?四川宜宾凉姜沟?四川泸州沱江二桥?湖北丹江口胡家岭?湖南长沙新港?湖南岳阳岳阳楼?湖北武汉宗关?江西南昌滁槎?江西九江蛤蟆石?枯水期平均丰水期平均平水期平均 图 1 把题目中分类标准中的前 5

21、类标准作为 5 个样本，也按上述综合指标的计算方法计算的出 6 个等级水质的量化标准（见表 4）表 4 劣 0.0 3 0.0 6 0.1 0.1 6 0.2 2 0.2 2 从图 1 中不难看出各个测点在 3 个时期的水质等级，比如：江西南昌滁槎的水质最差，在枯水期和平水期已经达到了劣类标准。在丰水期也已达到类标准。然后我们再把 3 个时期的平均指标再综合，得出长江最近一两年的水质综合评价指标（见图 2）3 个水期综合评价00.0 30.0 60.0 90.1 20.1 50.1 80.2 10.2 40.2 70.31234567891 0 1 1 1 2 1 3 1 4 1 5 1 6

22、1 71 7 个观测点 图 2 根据上表中的量化标准可以得出水质等级评价（见表 5）2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 7 表 5 测点 1 2 3 4 5 6 7 8 9 1 0 1 1 1 2 1 3 1 4 1 5 1 6 1 7 水质等级 劣 比较综合指数模型和模糊综合评判模型得出的结论，发现两者基本一致，所以可以判定我们模型是客观有效的 5.2 主要污染物 C O D 和 N H 3-N 的污染源分析 5.2.1 因为支流互不干扰且不受干流影响，所以支流可以先分离出来单独考虑。我们先找出长江支流10个测点从2004年

23、4月到2005年 4月13个月中COD 和NH两项指标的监测数据，然后求平均得到支流每个测点 COD 和 NH 的监测值 Czj和 Nzj（见附件图4），又由于C O D和 N H在 4项指标中的权重为 1 6.5%和 5 6.5%，所以它们两者之间的权重可以重新划分为 0.2 3 和 0.7 7。我们在把 Czj和 Nzj分别进行归一，然后再加权平均得到综合指标 Wj。Wj=0.23Czj/Czj+0.77 Nzj/Nzj 用柱形图表示（见下图 3）支流测点本地污染综合评价W j00.0 50.10.1 50.20.2 50.30.3 5测点1 234567891 0、图 3 我们取 Wj=

24、0.05 为比较标准。可以看出污染较严重的测点依次是：测点 8（江西南昌滁槎）、测点 3（四川泸州沱江二桥）、测点 1（四川乐山岷江大桥）、测点 5（湖南长沙新港）、测点 6（湖南岳阳楼）。因为各个支流之间是互不干扰的，所以我们可以初步判断这5 个地区为主要污染源。5.2.2 然后我们再来考虑干流上的污染源。我们假设同一个测点上的监测值在几天之内是不变的，也就是说，可以用同一时刻上上游的监测值经过几天降解后的剩余两表示其现时刻对下游的影响。对干流各个测点间降解天数（次数）m 的计算：各测点水流速度 Vj（1 3 个月平均）。（j=1 7）（单位由秒换算成天）测点之间水流平均速度：211,+=j

25、jjjVVV （j=1 6）各测点间距离：Sj，j+1；2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 8 所以 mj=Sj+1/1,+jjV 根据表中数据算出干流各测点 C O D 和 N H 1 3 个月的平均监测值 Cgj 和 Ngj 因为没有数据，我们在这里不考虑第一个测点上游对它的影响，认为它的监测值就等于本地污染。从第 2 点开始计算本地污染（降解）：Cg*j=Cgj-Cgj-1jm)2.01(Ngj*=Ngj-Ngj-1jm)2.01(所以本地污染 C O D 和 N H 的综合评价指标 Wj为（先把 Cg*j 和 Ngj*进

26、行归一，然后加权平均）：Wj=0.23Cg*j/Cg*j+0.77 Ngj*/Ngj*用柱形图表示（见图 4）干流测点本地污染综合评价W j00.0 50.10.1 50.20.2 51234567测点 图 4 因为干流水质整体情况优于支流，所以在确定干流污染源时不能用支流的比较标准。我们取 Wj=0.15。得出干流污染源为：测点 4（湖南岳阳城陵矶）、测点 2（重庆朱沱）。5.3 对长江未来水质污染的发展趋势的预测分析 由于河流水质状况存在许多的影响因素，而题中所给数据又十分有限，系统存在很强的不明确的因素，于是我们考虑用灰色系统法对其进行分析和预测。，灰色系统 G M（1，1）模型是依据系

27、统中已知的多种因素的综合资料，将此资料的时间序列按微分方程拟合去逼近上述时间序列所描述的动态过程，进而外推，达到预测的目的。这种拟合得到的模型是时间序列的一阶微分方程，因此，简记为 G M(1，1)模型。建立 G M(1，1)模型 给定原始时间 1 9 9 5-2 0 0 4 年资料列：x(0)=（x(0)(1)，x(0)(2)，x(0)(3)，x(0)(4)，x(0)(5)，x(0)(6)，x(0)(7)，x(0)(8)，x(0)(9)，x(0)(1 0)，）=（2 5 6 4,1 5 3 3,1 2 2 5,1 2 5 6,1 5 7 1,1 7 0 5,2 0 1 3,1 5 7 2,1

28、 8 2 4,4 7 3），对 x(0)坐 A G O 生成，有 x(1)=A G O x(0)，x(1)(k)=k1m)0(),m(x则 x(1)=（x(1)(1)，x(1)(2)，x(1)(3)，x(1)(4)，x(1)(5)，x(1)(6)，x(1)(7)，x(1)(8)，x(1)(9)，x(1)(1 0)，x(1)(1 1)，x(1)(1 2)）2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 9=（2 5 6 4,4 0 9 7,5 3 2 2,6 5 7 8,8 1 4 9,9 8 5 4,1 1 8 6 7,1 3 4 3 9,

29、1 5 2 6 3,1 5 7 2 6），对上述 x(0)的 G M（1，1）参数 a,b，按下述算式辨识：=(BTB)-1BTyN 基于 x(0)0(x与 x(1)，有=+=124434815.20484815.16836915.13594315.1078918673015.6858115.489151293081)2()1(5.01)2()1(5.01)2()1(5.01)2()1(5.01)2()1(5.01)2()1(5.01)2()1(5.01)2()1(5.01)2()1(5.01)10(1)9(1)8(1)7(1)6(1)5(1)4(1)3(1)2()1()1()1()1()1(

30、)1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1()1(xxxxxxxxxxxxxxxxxxzzzzzzzzzB yN=x(0)(2)，x(0)(3)，x(0)(4)，x(0)(5)，x(0)(6)，x(0)(7)，x(0)(8)，x(0)(9)，x(0)(1 0)T=（2 5 6 4,1 5 3 3,1 2 2 5,1 2 5 6,1 5 7 1,1 7 0 5,2 0 1 3,1 5 7 2,1 8 2 4,4 7 3）T。将 B,yN代入辨识算式，有.=baa=(BTB)-1BTyN=4115.640.14699

31、9-，a=0.1 4 6 9 9 9 b=4 1 1 5.6 4 得 G M(1，1)模型为 灰微分方程 x(0)(k)0.1 4 6 9 9 9 z(1)(k)=4 1 1 5.6 4；白化方程 dtdx)1(0.1 4 6 9 9 9 z(1)(k)=4 1 1 5.6 4；白化方程的时间响应式 =2 5 4 3 3.5 7te1062105.02 7 9 9 7.5 7，)1()1(),()1()1()1()0()1()1()0(xxtxtxtx=+=+=2 5 6 4 得还原方程：)1()0(+tx=1 5 1 5te-0.146999 预测结果（见表 6）abeabxtxat+=+)

32、1()1()0()1(2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 10 表 6 注：预测值=)k(x)k(x)0()0(，?=%)k(x)k(x)k(x)0()0()0(因 1998 年发洪水，数据做了相应调整 为相对误差（单位%）把预测值转换成比例（见表 7）：表 7 年份 I 类 类 类 类 类 劣类 2005 0.0 2 4 6 4 6 0.3 4 5 1 5 8 0.3 1 1 2 8 7 0.1 1 0 5 4 9 0.0 5 8 3 9 1 0.1 4 9 9 6 8 2006 0.0 2 0 6 7 4 0.3 4 0 1

33、 0 2 0.3 0 0 4 5 5 0.1 0 8 5 0 4 0.0 6 0 1 8 3 0.1 7 0 0 8 2 2007 0.0 1 7 2 8 4 0.3 3 3 8 1 8 0.2 8 8 8 8 7 0.1 0 6 0 8 3 0.0 6 1 8 0.1 9 2 1 2 8 2008 0.0 1 4 3 9 9 0.3 2 6 3 0 9 0.2 7 6 6 2 3 0.1 0 3 2 8 3 0.0 6 3 2 1 9 0.2 1 6 1 6 6 2009 0.0 1 1 9 3 5 0.3 1 7 6 2 9 0.2 6 3 7 6 7 0.1 0 0 1 3 6 0.0

34、6 4 3 6 0.2 4 2 1 7 3 2010 0.0 0 9 8 5 3 0.3 0 7 7 8 6 0.2 5 0 3 7 2 0.0 9 6 6 4 8 0.0 6 5 2 4 9 0.2 7 0 0 9 3 2011 0.0 0 8 0 9 5 0.2 9 6 8 5 2 0.2 3 6 5 4 6 0.0 9 2 8 5 1 0.0 6 5 8 3 5 0.2 9 9 8 2 2 2012 0.0 0 6 6 1 7 0.2 8 4 9 1 6 0.2 2 2 3 9 5 0.0 8 8 7 6 4 0.0 6 6 1 0 3 0.3 3 1 2 0 6 2013 0.0 0

35、5 3 8 0.2 7 2 0 7 9 0.2 0 8 0 4 1 0.0 8 4 4 2 9 0.0 6 6 0 4 0.3 6 4 0 3 1 2014 0.0 0 4 3 6 2 0.2 5 8 4 7 7 0.1 9 3 6 0 3 0.0 7 9 8 9 0.0 6 5 6 3 5 0.3 9 8 0 3 3 未来 10 年长江水质污染情况趋势用折线图表示（见图 5）年分 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 实际值 2564 1533 1225 1256 1571 1705 2013 1572 1824 473 预测值

36、2564 1544 1523 1240 1483 1463 1843 1424 1705 387 I类 0 0.007 0.243-0.012-0.056-0.141-0.084-0.094-0.065-0.181 实际值 4225 2030 2582 2645 12113 9944 11289 15574 15964 10602 预测值 4225 2303 2911 2847 7931 9187 10641 12326 14278 12539 类 0 0.134 0.127 0.076-0.345-0.076-0.057-0.208-0.105 0.182 实际值 2452 4998 435

37、7 5786 10738 10792 11855 10022 12048 15725 预测值 2452 5541 5287 6134 9095 9186 10423 11827 13420 15227 类 0 0.108 0.213 0.060-0.153-0.148-0.120 0.180 0.113-0.031 实际值 387 974 1329 905 2888 5037 4784 3544 2458 5833 预测值 387 869 1554 855 2601 4701 3463 3400 2610 5318 类 0-0.107 0.169-0.055-0.099-0.066-0.276

38、-0.040 0.061-0.088 实际值 297 191 259.8 186 1897 1320 1876 1136 2246 2325 预测值 297 175 300 205 1523 1240 1503 1821 2207 2675 类 0-0.083 0.154 0.102-0.197-0.060-0.198 0.602-0.017 0.150 实际值 0 311 339.7 180 1239 1608 2329 3538 3973 4454 预测值 0 326 334 112 1482 1976 2635 3351 4683 5244 劣类 0 0.0 4 8-0.0 1-0.3

39、7 0.1 9 6 0.2 2 8 0.1 3 1-0.0 5 0.1 7 8 0.1 7 7 2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 11 未来1 0 年长江水质（六类水比例）趋势00.10.20.30.40.52 0 0 5 年2 0 0 6 年2 0 0 7 年2 0 0 8 年2 0 0 9 年2 0 1 0 年2 0 1 1 年2 0 1 2 年2 0 1 3 年2 0 1 4 年I 类类类类类劣类 图 5 由图可以轻易看出，劣类水在未来 1 0 年中上升趋势很快，类水上升趋势较慢，其它类水成下降趋势。也就是说如果不加控制

40、，长江水质污染将会越来越严重。5.4 对长江水资源污染的控制分析 在问题三中我们是对长江全流域的水质情况做的预测。现在我们用同样的方法用过去 1 0 年长江干流的数据预测未来 1 0 的 6 类水的比例，以及总流量和废水排放量。然后我们重新把水分为 3 大类：清洁水：I 类、类和类 污染水：类和类 劣质水：劣类 各类比等于包含原类别的加和，所以 3 类水的比例，以及总流量和废水排放量（见表 8）表 8 年份 清洁水 x 1 污染水 x2 劣质水 x3 总流量 T 废水量 F 废水比例 f 2005 0.5 9 2 0 2 2 0.2 2 4 2 0 6 0.1 8 3 7 7 2 9611 3

41、03 0.0 3 1 5 2 6 2006 0.5 3 0 9 7 4 0.2 3 2 8 7 0.2 3 6 1 5 7 9541 322 0.0 3 3 7 4 9 2007 0.4 6 6 0 2 0.2 3 7 3 2 5 0.2 9 6 6 5 4 9471 343 0.0 3 6 2 1 6 2008 0.3 9 9 4 3 0.2 3 6 9 2 6 0.3 6 3 6 4 3 9402 365 0.0 3 8 8 2 2 2009 0.3 3 3 9 4 2 0.2 3 1 4 8 5 0.4 3 4 5 7 3 9334 388 0.0 4 1 5 6 8 2010 0.2

42、7 2 3 0 3 0.2 2 1 3 6 3 0.5 0 6 3 3 4 9266 413 0.0 4 4 5 7 2 2011 0.2 1 6 7 7 9 0.2 0 7 4 0 8 0.5 7 5 8 1 3 9198 440 0.0 4 7 8 3 6 2012 0.1 6 8 8 1 0.1 9 0 7 7 9 0.6 4 0 4 1 2 9131 468 0.0 5 1 2 5 4 2013 0.1 2 8 9 1 0 8 0.1 7 2 7 0 5 0.6 9 8 3 8 4 9065 499 0.0 5 5 0 4 7 2014 0.0 9 6 8 0 9 8 0.1 5 4

43、2 9 1 0.7 4 8 9 0 8999 531 0.0 5 9 0 0 7 然后根据废水比例 f 对三类水比例 x 1、x 2、x 3 进行逐步回归分析，得出逐步回归方程（见表 9 第 2 列）。按照题意，如果要控制污染水的比例在 0.2 以下，劣质水为 0，所以将 x 1=0.8、x 2=0.2、x 3=0 代入回归方程可计算出废水允许最大排放量 F*（见表 9 第 4 列）。那么需要处理的废水量F=F-F*，（见表 9 第 5 列）。2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 12 表 9 年份 逐步回归方程 F*（亿吨）F（

44、亿吨）2 0 0 5 年 f=-0.0 2 7 3 x 1+0.0 4 5 1 8 p0.0 2 0 3 2 2 3.6 2 8 7 7 9.3 7 1 2 5 2 0 0 6 年 f=-0.0 3 0 4 x 1+0.0 4 7 8 5 p0.0 0 3 5 2 2 4.4 9 9 7 9 7.5 0 0 2 7 2 0 0 7 年 f=-0.0 3 2 4 5 x 1+0.0 4 9 6 2 p0.0 0 0 5 2 2 4.0 8 3 9 1 1 8.9 1 6 1 2 0 0 8 年 f=-0.0 3 3 9 1 x 1+0.0 5 0 8 6 p0.0 0 0 1 2 2 3.1 2

45、 8 3 1 4 1.8 7 1 7 2 0 0 9 年 f=-0.0 3 4 9 8 x 1+0.0 5 1 7 6 p0.0 0 0 1 2 2 1.9 2 5 2 1 6 6.0 7 4 8 2 0 1 0 年 f=-0.0 3 5 8 9 x 1+0.0 5 2 5 p0.0 0 0 1 2 2 0.4 1 9 6 1 9 2.5 8 0 4 2 0 1 1 年 f=-0.0 2 2 1 4 x 1+0.0 2 0 3 2 x 3+0.0 4 0 4 0 p0.0 0 0 4 2 0 8.6 8 4 2 2 3 1.3 1 5 8 2 0 1 2 年 f=-0.0 2 1 7 4 x

46、1+0.0 2 1 1 5 x 3+0.0 4 0 0 3 p0.0 0 0 2 2 0 6.7 0 7 6 2 6 1.2 9 2 4 2 0 1 3 年 f=-0.0 2 0 7 3 x 1+0.0 2 3 0 4 x 3+0.0 3 9 1 1 p0.0 0 0 1 2 0 4.1 9 8 2 2 9 4.8 0 1 8 2 0 1 4 年 f=-0.0 1 8 8 4 x 2+0.0 4 5 1 1 x 3+0.0 1 8 5 6 p0.0 0 0 1 2 0 0.9 2 9 7 3 3 0.0 7 0 3 5.5 一些建议和意见 1、合理制定工业发展规划 人口在增长，社会在发展。绝对

47、阻止工业的发展是不现实的，我们只能在控制工业发展规模一定的情况下，根据长江目前污染状况，合理规划工业的地理位置，使长江的生态生命力（降解能力）达到最大。比如：在江西南昌滁槎地区，污染已经很严重，工业发展规模就应该相对减少。2、污水处理设施的利用 中国目前有相当一部分的污水处理场处于停止状态，没有发挥作用，主要原因是处理费用高，企业在追求经济利益的目的下以牺牲环境为代价。这就要求环保部门加大管理力度，不定期进行检查。加大处罚力度。每个工业企业建立自己的处理厂，严格执行许可证排污制度。3、控制排污总量 一个地区的污水排放不仅要控制污染物浓度，还要控制总量，每个地区的总量上限根据水质评测给出，总的排

48、放量不能超过上限，否则就要关闭一些企业。4、加大执法力度 环境保护并非新鲜事物，法规已很健全，关键在于执法不严。如果能够提高执法人员的素质，加强执法力度，（对违法乱纪现象予以严办！）还我山河指日可待。六、模型的评价 模糊综合评判法所做的对长江沿岸17个城市28个月的平均值所做的综合评价由于没有考虑季节对它们的影响，时间方面具有一定的局限性，同时由于题目给的评判水等级的项目只有 4 个，且 PH 值对其影响不大，所做出的评判难免有一定的误差，但通过用总体说来该模型还是比较成功的。通过用原理简单、常用的综合指数评价模型进行验证。发现结果吻合的很好，说明我们的模糊综合评价法是客观有效的。逐步回归模型

49、通过连续使用预测数据不断回归得到 15 个不同的回归方程，使所得到的结果更有说服力，避免了单一回归方程计算未来 10 年的废水率而带来的误差。2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学)蔡虎，卢亚军。许超 指导教师 汪晓银 13 参考文献 1 史晓新，夏军，水环境质量评价灰色模式识别模型及应用，http:/ 2 马建华，季凡，水质评价的模糊概率综合评价法，水文，第 3 期，第 21-24 页(1994)3 陆渝蓉，地球水环境学，南京，南京大学出版社 1993 年 附件：图1：00.10.20.30.40.5四川攀枝花龙洞重庆朱沱?湖北宜昌南津关?湖南岳阳城陵矶?江西九江河西

50、水厂?安徽安庆皖河口?江苏南京林山?四川乐山岷江大桥?四川宜宾凉姜沟?四川泸州沱江二桥?湖北丹江口胡家岭?湖南长沙新港?湖南岳阳岳阳楼?湖北武汉宗关?江西南昌滁槎?江西九江蛤蟆石?江苏扬州三江营?2 0 0 3 F2 0 0 3 P 图2：00.10.20.30.40.5四川攀枝花龙洞重庆朱沱?湖北宜昌南津关?湖南岳阳城陵矶?江西九江河西水厂?安徽安庆皖河口?江苏南京林山?四川乐山岷江大桥?四川宜宾凉姜沟?四川泸州沱江二桥?湖北丹江口胡家岭?湖南长沙新港?湖南岳阳岳阳楼?湖北武汉宗关?江西南昌滁槎?江西九江蛤蟆石?江苏扬州三江营?2 0 0 4 K2 0 0 4 F2 0 0 4 P 图3：0