6、全国大学生数学建模竞赛A0807

A题之一（全国一等奖）

长江水质的评价和预测模型

参赛学校：广西师范大学 参赛学生：黄勇萍、覃荣存、陶胜达

指导教师：建模指导组

摘 要

本文研究长江水质的评价和预测，对长江近两年多的水质情况做出定量的综合评价，并分析各观测点的污染状况，分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源，依照过去10年的主要统计数据，对长江未来10年水质污染的发展趋势做出预测，未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水。我们得到以下结果：

问题1：我们利用文献[1]给出的算法算出干流中各断面的水质综合指数，然后比较各观测点在不同时间的水质综合指数，对长江近两年来的水质情况做出定量的综合评价，并分析各地区水质的污染状况。

问题2：我们利用文献[2]中零维水质模型的基本方程，

''，求解出第n段河道排污量QnCn?QnCn?PQn?1Cn?1dCn???kCdtVnVn'Pn?aV'n?QnC''n?(Qn?1?Qn)Cn?Qn?1Cn?1，从而判断出04年4月至05年4月污染源

主要位于湖北宜昌南津关至湖南岳阳城陵矶和四川攀枝花龙洞至重庆朱沱之间。

问题3：利用10年长江总流量、废水排放总量，及长江流域枯水期、丰水期、水文年，各类水所占比例，利用文[3]中给出的GM(1,1) 模型对长江未来水质污染的发展趋势做出预测分析。

问题4：枯水期第i类水占干流的比例为fij,将fij?fij(Aj,Wj)近似看成Aj,Wj,WjAj的

一次函数，利用最小二乘法确定常数项和系数。未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,计算出每年可排放污水最大值，从而算出未来10年内每年需处理的污水如下表 年 处理废水量 2005 102.38 2006 119.78 2007 128.79 2008 146.10 2009 161.38 2010 184.23 2011 184.22 2012 220.83 2013 253.46 2014 271.31 问题5：我们简单分析了长江水污染严重的原因，对改善长江流域水环境提出了一些建议。

关键词：水质 污染物 浓度 流量

21

一、 问题的重述

１． 问题的提出

水是人类赖以生存的资源，保护水资源就是保护人类自己。长江是我国第一、世界第三的淡水河，它与我国国民生活息息相关。随着中国工农业的发展，长江水系受到污染，探讨与研究长江水质的状况，评价和预测长江水质，对长江水质污染发展趋势的预测以及提出解决长江水质污染问题方案有着紧密的关系。

原题中附件３给出了长江沿线１７个观测站（地区）近两年多主要水质指标的检测数据，以及干流上７个观测站近一年多的基本数据（站点距离、水流量和水流速），附件4是“1995~2004年长江流域水质报告”给出的主要统计数据。下面的附表是国标(GB3838-2002) 给出的《地表水环境质量标准》中4个主要项目标准限值，其中Ⅰ、Ⅱ、Ⅲ类为可饮用水。

附表: 《地表水环境质量标准》（GB3838—2002）中4个主要项目标准限值 单位：mg/L 分 类 序 标准值 号 项 目 7.5 1 溶解氧(DO) ≥ （或饱和率90%） 2 3 4 高锰酸盐指数(CODMn) ≤ 氨氮（NH3-N） ≤ PH值（无量纲） 2 0.15 4 0.5 6 1.0 6---9 10 1.5 15 2.0 ∞ ∞ 6 5 3 2 0 Ⅰ类 Ⅱ类 Ⅲ类 Ⅳ类 Ⅴ类 劣Ⅴ类 通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。要求研究如下问题：

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

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

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

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

（5）提出对解决长江水质污染问题有什么切实可行的建议和意见。 下图为附件３中给出的１7个观测站的地理图形：

22

大渡河岷江汉江金沙江8沱江丹江口水库11176武汉516九江鄱阳湖南昌赣江15安庆南京7j9宜宾1102泸州重庆宜昌34 岳阳13洞庭湖长沙1214江都扬州长江口上海湘江图解：粗线表示干流，细线表示支流，数字1~17依次分别表示：16四川攀枝花安徽安庆皖河口2重庆朱沱3湖北宜昌南天关4 湖南岳阳城陵矶5江西九江河西水厂10四川泸州沱江二桥15江西南昌滁槎7江苏南京林山8四川乐山岷江大桥13湖南岳阳岳阳楼9四川宜宾凉姜沟14湖北武汉宗关11湖北丹江口胡家岭12湖南长沙新港16江西九江蛤蟆石17江苏扬州三江营

二、符合说明：

为干流上第 j个观测点断面的综合水质指数； 为干流上第 j个观测点断面综合水质类别

pij为j端面i项污染指标的分指数

为i项污染指标的权重

记第 n段河段为从干流上第n个观测站到第n＋1个观测站之间的河段

Sn为第 n段河段长度

tn表示水流经过第 n段河段所花的时间 vn表示水流经过干流上第 n个观测点的速度 an为第n段中河水加速度 k为污染物的降解系数

Pn第 n段河段污染物每天流入量（g/天） Qn为干流第n个观测站的流量

Cn为某种污染物（高锰酸盐指数或氨氮）的浓度

' Qn第 n段河段支流的流量

'第 n段河段支流该种污染物的浓度 CnAj为长江第j年总流量 Wj为第j年废水排放总量

fij枯水期第i类水（Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ，劣Ⅴ类分别记为1，2，3，4，

23

5，6），占干流的比例

ej每年需要处理的污水量

本文对问题(1),(2),(3)，（4）分别建立模型，来解答这四个问题。

问题1、

1、问题的分析

根据附件3的定性分析，我们初步知道各个观测点每月的水质类别，以及我们要监测的各主要监测项目的单位含量，依据这些数据，我们对长江水质近两年多的水质情况做出定量分析。文献[1]给出了水质综合指数的一种算法，我们利用[1]给出的算法算出干流中各断面的水质综合指数，然后比较各观测点在不同时间的水质综合指数，以及观测点与观测点之间在相同时间内的水质综合指数，再对长江近两年来的水质情况做出定量的综合评价，并分析各地区水质的污染状况。

2、模型的建立及求解

在对长江近两年多的水质情况做出定量的综合评价时，我们首先需要知道这两年多时间里各断面在各个时间段的综合水质指数。为此，我们采用文献[1]给出的模型：

式中，为j断面的综合水质指数；为断面综合水质类别（Ⅰ，Ⅱ，Ⅲ，Ⅳ，Ⅴ及劣Ⅴ类分别对应1，2，3，4，5，6）；Q为系数，它的作用是满足：

Q??i?1nwipij?w?1

i pij为j端面i项污染指标的分指数（为浓度值与标准值的比值，这里取Ⅰ类标准值）；为i项污染指标的权重。本小问中取溶解氧(DO)指数、高锰酸盐指数(CODMn)、氨氮（NH3-N）、PH值（无量纲）对应的i分别为1，2，3，4。由于溶解氧(DO)指数越大越好，实际计算时取p1j＝10－对应观测点的溶解氧(DO)指数。

这样综合水质指数就由证书和小数两部分构成：整数部分代表端面水质功能类别，使综合水质指数与水质类别统一起来；小数部分为各污染指数的加权平均，以区别同类水质的优劣，提高了分辨率。

的算法如下：wi?si1 si5这种定义权重的方法是以污染物超标倍数对河流水质的贡献率大小为依据的。例如，我们取附表中的高锰酸盐和氨氮的标准值来说明该算法： 高锰酸盐的Ⅰ类标准值为2，Ⅴ类标准值为15，就得到wi高?si1=2/15 si5氨氮的Ⅰ类标准值为0.15，Ⅴ类标准值为2.0，就得到wi氨?

24

qsi1=0.15/2.0 si5公式中Q是经验系数，要考虑到小数既不能过小，又要具有较高的分辨率。Q值选取如下：

当水质类别为Ⅰ类时，Q=1； 当水质类别为Ⅱ类时，Q=0.47； 当水质类别为Ⅲ类时，Q=0.45； 当水质类别为Ⅳ类时，Q=0.41； 当水质类别为Ⅴ类时，Q=0.18； 当水质类别为劣Ⅴ类时，Q=0.17。

当水质类别为Ⅴ类时，可能出现指数过大，Q=0.17仍不能使

Q?wipij?wi?1，则规定

Q?wipij?w?0.999

i利用该模型算出四川攀枝花、重庆朱沱、湖北宜昌南、湖南岳阳城、江西九江河、安徽安庆皖、江苏南京林7个观测点的水质综合指数表，并根据相应的数据画出对应的图形如下：

图一四川攀枝花综合水质

从图一中可看到该观测点水质污染状况：Ⅰ类水占28.57%，Ⅱ类水占53.57%， Ⅲ 类水占10.71%，Ⅳ类水占7.14%。由表可看出同类水中水质仍有差别,例如：03年6月和04年4月同是Ⅱ类水，但是前者综合指数为2.0489，后者为2.2237，相比可知前者水质更好。

25

图二重庆朱沱综合水质

从图中可看到该观测点水质污染状况：Ⅰ类水占0%，Ⅱ类水占71.43%, Ⅲ 类水占28.57%，Ⅳ类水占0%.由表可看出同类水中水质仍有差别,例如：03年6月和7月同是Ⅱ类水，但是前者综合指数为2.3256，后者为2.2049，相比可知后者水质更好。

图三湖北宜昌南津关综合水质

从图中可看到该观测点水质污染状况：Ⅰ类水占0%,Ⅱ类水占89.29%, Ⅲ 类水占

26

10.71%,Ⅳ类水占0%.由表可看出同类水中水质仍有差别,例如：04年4月和5月同是Ⅱ类水，但是前者综合指数为2.4162，后者为2.2568，相比可知后者水质更好。

图四

从图中可看到该观测点水质污染状况：Ⅰ类水占0%,Ⅱ类水占64.29%, Ⅲ 类水占35.71%,Ⅳ类水占0%.由表可看出同类水中水质仍有差别,例如：05年4月和5月同是Ⅱ类水，但是前者综合指数为2.3602，后者为2.2899，相比可知后者水质更好。

图五江西九江河西水厂综合水质

27

从图中可看到该观测点水质污染状况：Ⅰ类水占7.14%,Ⅱ类水占92.86%, Ⅲ 类水占0%,Ⅳ类水占0%.由表可看出同类水中水质仍有差别,例如：03年6月和04年4月同是Ⅱ类水，但是前者综合指数为2.1718，后者为2.3094，相比可知前者水质更好。

图六安徽安庆皖河口综合水质

从图中可看到该观测点水质污染状况：Ⅰ类水占3.57%,Ⅱ类水占92.86%, Ⅲ 类水占3.57%,Ⅳ类水占0%.由表可看出同类水中水质仍有差别,例如：03年6月和04年4月同是Ⅱ类水，但是前者综合指数为2.3151，后者为2.285，相比可知后者水质更好。

图七江苏南京林山综合水质

28

从图中可看到该观测点水质污染状况：Ⅰ类水占7.14%,Ⅱ类水占89.29%, Ⅲ 类水占3.57%,Ⅳ类水占0%.由表可看出同类水中水质仍有差别,例如：03年7月和04年5月同是Ⅱ类水，但是前者综合指数为2.2218，后者为2.1658，相比可知后者水质更好。 4、模型结果分析：我们采用文献[1]给出的水质综合指数的算法，该算法能比较准确的反映水质的真实状况，比直接确定水的类别能较为准确的反映水质的真实状况。 5、模型改进：Q值的选取、wi的确定有待改进。 问题二：

符号说明：

记第 n段河段为从第n个观测站到第n＋1个观测站之间的河段

Sn为第 n段河段长度

tn表示水流经过第 n段河段所花的时间 vn表示水流经过干流上第 n个观测点的速度 an为第n段中河水加速度 k为污染物的降解系数

Pn第 n段河段污染物每天流入量（g/天） Qn为干流第n个观测站的流量

Cn为某种污染物（高锰酸盐指数或氨氮）的浓度

' Qn第 n段河段支流的流量

'第 n段河段支流该种污染物的浓度 Cn 问题的分析：已知干流上每个观测点断面水流量、流速、污染物浓度，要判断是否干流上两个观测点之间是否有污染物排放。 模型假设：

1、污染物的降解系数k为常数；

2、第 n段河段中江水流速均匀改变，即an为常数； 3、污染物的排放是均匀排放。 模型建立：

记第 n段河段为从干流上第n个观测站到第n＋1个观测站之间的河段，Sn为第 n段河段长度，设tn表示水流经过第 n段河段所花的时间，在第n段中河水加速度为an，

tnvn?1?vn则有an?,设在第n个观测站时t?0,v(t)?vn?ant，因此有?(vn?ant)dt?Sn,

0tn从而有(vnt?tn?v?v12tnant)|0?Sn,即 n?1ntn?Sn,求出水流经过第 n段河段所花的时间 222Sn。

vn?1?vn 29

长江干流第n个观测站的流量为Qn，流速为vn，某种污染物（高锰酸盐指数或氨

'氮）的浓度为Cn，该种污染物的降解系数为k,若在n个观测站处支流的流量为Qn，该'种污染物的浓度为Cn，污染源的该种污染物每天流入量为Pn，在第n段流量满足方程

'Qn?1?Qn?Qn （1） ''QCQnCn?QnCn?Pn该种污染物浓度增加量为，减少量为n?1n?1?kC，浓度变

VnVn化满足方程：

'' （2） QnCn?QnCn?PQCdCn??n?1n?1?kCdtVnVn其中Vn为第n个观测站到第n+1个观测站的水体体积。（2）也可参见文献[2]第三章第二节零维水质模型的基本方程。

Q?Qn?1tn表示，由（1）我们近似的用Vn?n，（2）得到 2

'QnCn?(Qn?1?Qn)Cn?PQCdCn ??n?1n?1?kC （3）dtVnVn不妨设支流和污染物从第 n段河段中间流入，可以想象为中间设了一个观测点，将第 n段河段分为两段，记中间虚拟的观测点为n'，设河水流经第一段、第二段的时间分别为tn,1,tn,2，则tn,1?tn,2段满足

dC??kC （4） dt?tn，第一段由于没有支流和污染物注入，河水在第一

到第n'个观测点后，该种污染物的浓度为

C?Cne''n?ktn,1，

'QnC''n?(Qn?1?Qn)Cn?PQCdCn??n?1n?1?kC （5） dtV'nV'n其中V'n为第二段水体的体积，V'n?Qn?1tn,2，即有

'QnC''n?(Qn?1?Qn)Cn?PdCn?Qn?1Cn?1??kC （6） dtV'n'QnC''n?(Qn?1?Qn)Cn?Pn?Qn?1Cn?1为简便起见记?a，

V'n 30

dC?a?kC （7） dt解得C(t)?a?me?kt，其中m为常数。 ka?kt''t?0时，C(0)=?m?Cn?Cnen,1 （8）

kt?tn,2时，C(tn,2)=由（8）得

a?kt?men,2?Cn?1 （9） km?Cne代入（9）得，

?ktn,1a－k （10）

aa?ktn,1?kt?（Cne－）en,2?Cn?1 （11） kka?ktn,2?ktn（1－e)?Ce?Cn?1 （12） 化简得 nk求得

k(Cn?1?Cne?ktn)a=?kt1－en,2，

t其中n,2由

?tn,20(vn?1?vn?1?vn1t)dt?Sn确定，即

tn2vn?1tn,2?1vn?1?vn21tn,2?Sn。

2tn2'由a的表达式得QnC''n?(Qn?1?Qn)Cn?Pn?Qn?1Cn?1?aV'n，从而可求出

'Pn?aV'n?QnC''n?(Qn?1?Qn)Cn?Qn?1Cn?1。 （13）

若Pn?0，则有污染源，否则就没有。

模型求解：实际计算时取k?0.2, 计算得所有的Pn值，见附件2。判断出04年4月至05年4月污染源主要位于湖北宜昌南津关至湖南岳阳城陵矶和四川攀枝花龙洞至重庆朱沱之间。

模型检验与结果分析：我们将求得的Pn值代入（3）式，求解方程（3），由（10）求出

m，由（9）确定C(tn,2)，比较C31

n?1与C(tn,2)，发现误差非常小。

模型评价与改进：由于没有提供各支流流入长江入口的流量、流速及污染物浓度和流入点距离观测站的距离，本模型假设支流与污染物均从每个河段的中间注入，如能提供上述数据，则模型更完善。 问题三、

问题的分析：已知过去10年长江总流量、废水排放总量，及长江流域枯水期、丰水期、水文年，各类水所占比例，要求对长江未来水质污染的发展趋势做出预测分析，研究未来10年的10年长江总流量、废水排放总量，及长江流域枯水期、丰水期、水文年，各类水所占比例。

符号说明：

记G?0??t?为一组随时间变化的河流水质实测数据序列

模型假设：假设文中X?1?具有近似指数变化规律。

文[3]指出灰色系统中应用最为广泛的GM (1, 1) 模型自20 世纪80 年代邓聚龙教授提出以来, 由于其所需样本数据少, 计算简便等优点, 已广泛应用于社会、经济、生态、农业等各个领域, 尤其在小样本、贫信息不确定系统和缺乏数据的情况下, 也得到了成功的应用, 这就决定了灰色系统中的GM (1, 1) 模型在预测、决策等领域占有重要的地位。

我们采用GM (1, 1) 模型来预测未来10年长江总流量，废水排放总量，每年枯水期、丰水期、水文年各种水质的百分比。下面GM (1, 1) 模型见文献[3]。

设G?0??t?为一组随时间变化的河流水质实测数据序列，采用自然对数变换对G?0??t??0???0??0?G?G1,G???2?,?0??进行光滑处理，得到x，即：

,G?0??n???x?0??t??lnG?t??0?

x?0??x?0??1?,x?0??2?,?,x?0??n?

?由附件４给出的过去１０年对长江水质的调查监测的主要统计数据，研究发现长江水质受季节因素影响，而在同一季节这种影响波动又具有相似的波动趋势。我们针对这些特征建立GM(1,1)模型数学模型对长江水质进行分析与预测。

模型的建立：

000由模型假设中可知原始数据序列为x???x???1?,x???2?,?,x?0??n?,对x?0?进行一次

?累加生成，得到生成序列为：

x?x1?1???1??1?,x?2?,?1?,x?1??n??, 其中

x?1??t???x?0??i?,???????i?1,2,i?1t,n

假设X??具有近似指数变化规律，则灰色微分方程的白化方程初值问题为：

??1?dx1??ax???u,??dt??1? （14） 0????x?1??x?1?,

32

（1）式中的a,u为待定参数

将上式离散化，微分方程变为差分方程，得到均值生成序列为

Z?1??t??Z?1??2?,Z?1??3?,?,Z?1??n?

?从而得到GM（1，1）灰差分方程如下：

X?0??t??aZ?1??t??u. （15）

按照最小二乘准则采用最小二乘法可求得参数

a??a,u???BTB?BTYN ， （16）

T?1其中

)?x?0??2????Z(1?2?1????(1)??0??x?3????Z?3?1??????? ， YN?? B???. ?????)??0????Z(1?n?1???x?n??Z?1??t???x?1??t?1???1???x?1??t?,????t?2,3,,n,?????0,1?, （17）

取?=0.5，求出参数a、u，解（1）式得初值问题的解

u??at?1u?01x???t???x???1???e??? （18）

a?a?对x?1??t?作一次累减生成，还原到原始数据，得

x?0??0??t????x?1???u??a?t?1?a,? t=2~n. (19) ??1?e?ea?模型求解：实际算出未来10年长江总流量，废水排放总量，每年枯水期、丰水期、水文年各类水质的百分比，见附件4。

模型的检验：我们还原到原始数据，发现预测的1995－2004年各项数据与实际已知数据差距较小，结果令人满意。 模型评价与改进：正如文[3]所指出的GM (1, 1) 模型具有所需样本数据少, 计算简便等优点。

问题4、

问题分析：要求未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水。

符号说明：

Aj为长江第j年总流量 Wj为第j年废水排放总量

，fij枯水期第i类水（Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ，劣Ⅴ类分别记为1，2，3，4，5，6）占干流的比例

33

ej每年需要处理的污水量

W模型假设：枯水期第i类水占干流的比例为fij近似为Aj，Wj，j的一次函数。

Aj经观察可看出枯水期长江干流的Ⅳ类和Ⅴ类水及劣Ⅴ类水占干流的比例一般高于丰水期占干流的比例，要求未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,为简便起见，我们要求枯水期长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水。

模型建立及求解：长江第j年总流量为Aj，废水排放总量为Wj，枯水期第i类水（Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ，劣Ⅴ类分别记为1，2，3，4，5，6），占干流的比例为

fij,则fij为Aj，Wj的函数，即fij?fij(Aj,Wj).

fij?fij(Aj,Wj)?bij?cijAj?dijWj?eij利用最小二乘法，确定bij,cij,dij,eij见附件5。

WjAj， （20）

要求未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,即要求

f4j?f5j?0.2，f6j?0，设每年实际排放的污水为W'j，则

W'jAj?b5j?c5jAj?d5jW'j?e5jW'jAj?0.2b4j?c4jAj?d4jW'j?e4j，b6j?c6jAj?d6jW'j?e6jW'jAj?0。我们可确定第j年可排放的污水最大值

W'j。

设每年需要处理的污水为ej，则实际排放的污水为W'j?Wj?ej，即

2013 253.46 2014 271.31 ej?Wj?W'j，计算得未来10年内需处理的污水如下表：

年 处理废水量 2005 102.38 2006 119.78 2007 128.79 2008 146.10 2009 161.38 2010 184.23 2011 184.22 2012 220.83 模型的检验：我们将算出的每年排放污水W'j用W'j代替（20）中Wj，算出

?Wj?ej，

f4j?f5j?0.2，f6j?0。

34

模型评价与改进：本模型只要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，

Wj且没有劣Ⅴ类水。并假设枯水期第i类水占干流的比例为fij近似为Aj，Wj，的

Aj一次函数。

问题（5）

根据题目提供的附件1和附件2的叙述信息以及对问题（1）、（2）、（3）、（4）的探讨与回答，可以清楚的看到，长江水污染已经达到了相当严重的程度，长江沿岸废水的排放如不加以控制或处理，5年后长江流域水质Ⅳ类和Ⅴ类水的比例将达到30%，劣Ⅴ类水将达到10％以上。这将极大影响到长江流域人们的生存。导致长江水质严重变坏的，有以下几个最根本的主要的缘故：

１． 民众、领导环境保护意识相当差，忧患意识低，且大多只顾追求眼前利益而

置长远利益于不顾。 ２． 工业的发展，尤其是重污染企业的建立与发展，使得工业三废与生活污物直

接且过多的排放到长江。

针对以上几点，我们对处理长江水质污染问题提出以下的建议与意见：

１．首先，必须唤起民众、尤其是沿江城市的民众的关注，加强他们的环保意识，

使他们清楚的知道，环境保护不只是国家的事情，也是每个人应该做且必须做的事情。

２．长江水质严重变坏最直接的原因就是长江沿岸有许多操作不规范的重污染企业，没有对污物进行处理就直接排放入河中，许多生活垃圾也直接倒进河内，针对此现象，应制定由目前条块分割的管理方式逐步过渡到集开发、利用和保护管理于一体的企业化管理体制。根据水体功能，制定合理的水质目标和相应的地方水环境质量标准、污染物排放标准；推行排放标准的年限制；在不同经济水平地区实行有区别的水环境管理制度。逐步推行总量控制和排污许可证制度，有偿使用环境容量。应用市场机制，制订实行水的有偿使用，制定合理的水资源价格政策、排污交易政策、配套法规和标准。工业水污染防治要优先发展全生产过程控制技术，城市和乡镇水污染防治要规划先行，坚持区域综合治理和流域综合防治、污染集中控制的方向，建设示范工程。

参考文献：

[1] 梁德华，蒋火华. 河流水质综合评价方法的统一和改进. 中国环境监测, 18（2）:63~66, 2002. [2]金腊华，徐峰俊，水环境数值模拟与可视化技术，化学工业出版社，2004年1月

[3] 何文章,宋国乡,吴爱弟. 估计GM(1,1)模型中参数的一族算法. 系统工程理论与实践, 1:69~75，

2005.

[3] 李如忠,汪家权,钱家忠. 基于灰色动态模型群法的客流水质预测研究. 水土保护通报,

22(4):10~12 ,2002.

[4] 李如忠,王超. 灰色动态模型群法在河流水质预测中的应用初探. 中国农村水利水

电.1:76~78,2003.

评委简评

1、文章各项目齐全，顺序安排合理，格式规范。文字表达较清晰，写作有条理，对每一问题的处理基本做到有始有终。

35

本文来源：https://www.bwwdw.com/article/ikrd.html