长江水质的评价和预测 (全国一等奖)

长江水质的评价和预测

摘要

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

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

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

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

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

一、问题重述

自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. 预测不考虑突变因素，如洪水、干旱等。

三、符号约定

ij

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

ik X ： 第k 观测站第i 监测项目在28个月中的平均值 (173,2,1 =k ) ijk

Y ： 第k 观测站第i 监测项目对第j 类水的隶属度 ik

W ： 第k 观测站第i 个监测项目的权重。

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

ik

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

ik

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

k

： 长江干流的降解系数

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

k

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

n Y ： 未来第n 年需要处理的污水量(103,2,1 =n ) n p ：未来第n 年的废水排放总量

k t ：

分别表示未来10年每年第Ⅳ、Ⅴ和劣Ⅴ类水的百分比(321、、

=k )

四、模型建立与求解

(一)、问题(1)

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

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

2.模型建立

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

k k k R A B ⨯= (1-1)

3．模型求解

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

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

⎪⎪⎩

⎪⎪

⎨⎧≥<<--≤=2

21122

1

101i ik i ik i i i ik i i ik k

i S X S X S S S X S S X Y (1-2) 5,4,3,2=j 类水：

⎪⎪⎪⎪⎩⎪⎪⎪⎪

⎨

⎧<≤<<--<≤--=+---+++)

1()1()1(12)1()

1()1()1(,0

j i ik j i ik ij

ik j i i i j i ik j i ik ij ij j i ik j i ijk

S X S X S X S S S S X S X S S S X S Y (1-3)