长江水质的评价和预测（全国一等奖） - 图文

长江水质的评价和预测

摘要

本文用模糊数学的方法，通过计算各评价因子的隶属度和权重，得到了长江近两年多的水质情况的综合评价结论：Ⅰ类水比例为25%、Ⅱ类水比例为23%、Ⅲ类水比例为20%、Ⅳ类水比例小于1%、Ⅴ类水比例为30%、劣Ⅴ类水比例小于2%，如下面饼图，其中可饮用水比例为68%，不可饮用水比例为32%。结果显示不可饮用水的比例很大，可以说明长江污染情况已经相当严重。

对于问题(2)，我们通过建立反映长江水质的一维稳态微分方程模型，并求解得到各观测站浓度的计算公式，用Matlab编程计算，计算结果显示,高锰酸盐污染源主要在：湖北宜昌南津关和湖南岳阳城陵矶。氨氮污染源主要在：重庆朱沱和湖南岳阳城陵矶。

对于问题(3)，根据近10年的水文年数据建立灰色系统预测模型,得到了未来10年长江全流域、干流、支流河长百分比的值，据此画出相应的走势图，由此确定水质污染的发展趋势，我们的结论是：长江未来10年的污染会越来越严重。

对于问题(4)，我们首先建立排污量的灰色系统预测模型，得出未来10年的排污总量，根据长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水的要求，建立了每年需要处理污水量的计算公式，得到了未来10年每年需要处理的污水量，见下表(单位：亿吨)：

年份 需处理的污水量 年份 需处理的污水量 2005 19.5390 2006 132.6777 2007 143.4669 2008 193.9729 2009 225.4419 2010 276.0293 2011 316.5134 2012 367.8958 2013 413.7736 2014 463.8682

关键词：模糊数学 隶属度 权重 微分方程 灰色系统

1

一、问题重述

自2004年10月“保护长江万里行”行动发起后，考察团对沿线21个重点城市做了实地考察，认识到了母亲河长江受到了严重的污染，为此，专家提出了拯救长江的呼唤，给出了下面这些有待解决的问题。

(1)对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。

(2)研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?

(3)假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。

(4)根据你的预测分析，如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,那么每年需要处理多少污水？

(5)你对解决长江水质污染问题有什么切实可行的建议和意见。

题目附件中给出了解决上述问题的各类数据。附件3给出了长江沿线17个观测站(地区)近两年多主要水质指标的检测数据，以及干流上7个观测站近一年多的基本数据(站点距离、水流量和水流速);附件4是“1995~2004年长江流域水质报告”给出的主要统计数据;附表是《地表水环境质量标准》中4个主要项目标准限值，其中Ⅰ、Ⅱ、Ⅲ类为可饮用水。污染物高锰酸盐指数和氨氮的降解系数，考虑取0.2 (单位：1/天)。

已知条件：通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水；污染物都有一定的自然净化能力(指标称为降解系数)；自然净化能力可以认为是近似均匀的。

二、模型假设

1．污染物排放入长江后迅速混合在水中。

2. 把长江认为是一维的，不考虑河宽，水深，横断面。

3. 主要污染物高锰酸盐指数和氨氮的降解系数取为常数0.2。 4. 一个地区的污染只来自于上游的污水和本地区的排污。 5. 预测不考虑突变因素，如洪水、干旱等。

Sij

三、符号约定

： 第i监测项目第j类水的标准限值 (i?1,2,3,4,j?1,2,3,4,5,6)

： 第k观测站第i监测项目在28个月中的平均值 (k?1,2,3?17) ： 第k观测站第i监测项目对第j类水的隶属度 ： 第k观测站第i个监测项目的权重。

XikYijkWikBkRkcik： 第k观测站模糊数学方法综合评价的结果 Ak： 第k观测站各评价因子的权重向量 ： 第k观测站隶属度的模糊关系矩阵

： 第k观测站第i个监测项目在近18个月中的平均浓度 (mg/L) c0 ： 各污染物的初始浓度

Cik： 第k观测站第i个监测项目浓度的计算值

k ： 长江干流的降解系数

x ： 长江干流相邻观测站间的距离(m)

2

uk： 第k观测站在近12个月中水流平均流速(m/s)

Yn： tk：

未来第n年需要处理的污水量(n?1,2,3?10) pn：未来第n年的废水排放总量

2、3) 分别表示未来10年每年第Ⅳ、Ⅴ和劣Ⅴ类水的百分比(k?1、

四、模型建立与求解

(一)、问题(1)

此问用模糊数学法对长江水质情况作定量的综合评价。 1.模糊数学法基本原理：

把普通集合理论中的非“0”则“1”的绝对隶属函数用[0，1]来刻画。 本文选用监测项目集合U?{ PH值（无量纲）、溶解氧(DO) 、高锰酸盐指数(CODMn)、 氨氮（NH3-N）}4项作为评价因子。由于PH值标准限值是一个区间范围，无法用模糊数学法对其做定量分析，而且附件3的数据显示，各观测点的PH值都位于区间[6，9]内，所以此小题只对后3项评价因子做综合评价。

2.模型建立

模糊数学的综合评价是通过模糊关系矩阵Rk和评价因子的权重矩阵Ak的复合运算进行评价，实际上是对各项评价因子进行加权合成，所以可以得到模型：

Bk?Ak?Rk (1-1)

3．模型求解

3.1 计算各评价因子对各类水的隶属度

用线形隶属函数确定各评价因子对各类水的隶属度的计算公式如下： j?1类水：

?1??S?Xik??i2?Si2?Si1?0??Si(j?1)?Xik??Si(j?1)?Sij?X?S?iki(j?1)???Si2?Si1????0Xik?Si1Si1?Xik?Si2 (1-2) Xik?Si2 Yi1kj?2,3,4,5类水：

Sij?Xik?Si(j?1) YijkSi(j?1)?Xik?Sij (1-3)

Xik?Si(j?1),Xik?Si(j?1)

3

j?6类水：

Yi6k?1??X?Si5??ik?Si6?Si5?0?Xik?Si6Si5?Xik?Si6 (1-4) Xik?Si5

由上述公式(1-2)，(1-3)，(1-4)可计算求出各观测站中各监测项目对于各类水的隶属度，它们都是1行6列的向量,将各观测站的向量按行放到一个矩阵中，得到各观测站的模糊关系矩阵Rk，以第1个观测站为例：

0 0 0 0 0 0????R1?0 0 0 0 0 0 ???? 0.9061 0.0939 0 0 0 0 ??用同样的方法我们可以求得另外16个观测站的模糊关系矩阵Rk。

3.2 各评价因子权重的计算

本文给出的4项评价因子中，溶解氧(DO)在水中含量越高表明水质越好，而高锰酸盐指数(CODMn)和氨氮（NH3-N）在水中含量越低表明水质越好，所以我们取以下的权重公式：

溶解氧(DO)权重公式为：

6?SW2k?j?12jmX (1-5)

2k高锰酸盐指数(CODMn)、 氨氮（NH3-N）权重公式为：

Wik?mX6ik (1-6)

?Sj?1ij其中m为水质的类数，即m?6，且i?3,4。

用上面公式(1-5), (1-6)计算，得到各观测站各评价因子的权重和可能出现大于1的情况，但模糊数学运算只允许在[0，1]区 …… 此处隐藏：2959字，全部文档内容请下载后查看。喜欢就下载吧 ……