小议长江水质的评价和预测

摘要：本文旨在研究长江过去1O年的水质变化情况，以此来预测未来1O年长江的水质情况，并对长江的治理提供一些具有可行性的建议和意见。在第一问中，我们从模糊数学的角度出发提出综合污染指数的概念，并运用该指数分析各地区的水质污染情况，得出结论，04年较O3年恶化，O5年有所好转；在第二问中．我们将7个长江干流观测点的水质报告表和基本数据表相结合，用每秒流过观测点的水中所含污染物的量减去上一个观测点的水中所古污染物经过白净后残余的量，即为两个观测点之间污染物增加的量，得出：高锰酸盐和氨氮的虽主要污染源在湖南岳阳城及其上游地区；在第三问中，我们建立了针对各类水所占评价河长百分比的一次累加拟合模型和时间序列模型，得出结论：1O年之后，I、Il类水都已不存在，III类水仅存O.62％；在第四问中，我们采用多元线性回归，得出长江流量，废水排放总量，Ⅳ类，V类，Ⅵ类水占水文年全流域长度半分比之间的线性关系，以第三问中预测数据为基础，以未来l0年总共要处理的污水为目标，建立线性规划并求解。

关键词:归一化法；自净系数；一次累加拟合模型；时间序列法；多元线性回归模型

1基本假设

1）假设主要污染物高锰酸盐指数和氨氮的自身降解系数为0.2（单位：1/天）；2）假设观测点j和观测点j＋1之间江水的流速是所测两点j和j＋1的速度的平均值；3）假设忽略观测点j和观测点j＋1之间增加的污染物的自净；4）假设2005～2014年的评价河长与2004年的评价河长相等；5）假设在2005～2014年间不会发生大旱大涝等自然灾害性天气；6）假设长江水的总流量，废水排放总量,IV类、V类、劣V类三种类别的水在水文年全流域的百分比之间成线性关系。

2问题分析

问题一中数据较多，我们可以用统计的方法，求出每年各类水所占的一个百分比，以此来说明长江水质的一个整体变化趋势。针对各地区而言，我们可以把影响水质的四个主要因素进行加权，从而求出一个综合污染指数，以此为标准来判断各地区水质的污染状况。在问题二中，由于江水在流动过程中对污染物可以进行一定的自身降解，所以水质最差的地区不一定就是污染源，如果上游污染严重，那么经过积累作用，到达下游的时候水质也会变得很差。仔细分析附表四，发现每年评价的河长差别很大，不利于我们建模求解。我们认为不妨从各类水所占百分比入手，预测出未来10年各类水所占的百分比即可。问题四的求解需要建立在问题三的基础上，由于题目要求没有劣V类水的存在，因此很容易得出劣V类水要全部进行处理。

3模型的建立与求解

3.1问题一

首先，我们从整体着手，对进两年多来长江水质进行整体的定量分析，统计每年长江水质按照水质标准分类而成的各种水类其绝对数值和相对百分比，得出长江水质04年较03年均比较差，但05年水质恶化现象得到了一定的控制，其水质明显好于04年。两年来，长江水质呈现出先恶化再好转的趋势。然后，对于各地区水质的污染状况，我们综合考虑对评判水质的四个因子：溶解氧、高锰酸盐指数、氨氮和PH值，采用模糊数学的方法对四者进行加权，从而得到一个综合污染指数W以这个指数来定量说明各地区水质的污染情况。我们以I类水：溶氧量7.5mg/L，高锰酸盐2mg/L，氨氮0.15mg/L，ph值取6～9的中值7.5为标准提出了一个综合污染指数的概念

很明显，按照综合污染指数函数的定义式，综合污染指数越高，水质越差。由以上表格（表格略）分析得出各地区水质污染状况，总体来讲，17个观测点04年水质比03年水质明显恶化，05年水质有较大程度的好转。

3.2问题二

在问题二中，我们把两观测点间汇入支流的主要污染物高锰酸盐指数和氨氮的自身降解忽略不计，并把水质主要污染物高锰酸盐指数和氨氮的自身降解系数（以下简称自净系数）假设为0.2（单位1/天）。我们根据附件三中2004.4-2005.4的长江主干水质监测报告测算出2004.4-2005.4间长江水流过两相邻观测点所需时间：第j点到j＋1点距离为……。以两个观测站点的速度平均值……。

作为两站点间水流速度，时间即为由每个点的水流量我们可以计算出在各观测点每秒流过的水中所含有的CODMn的量为：CODMn的浓度流量，即水流经过自净到达下一观测站点时，每秒流过的水中所含有的CODMn的量为：CODMn的浓度流量自净系数的天数次幂，即两站点间增加的CODMn的量

我们在用进行数值计算得过程中，增加量有部分为负值，而实际当中负值是不可能存在的，因此将负值全部修正为0。2004年4月～2005年4月13个月CODMn的平均增加值

同理，根据上面的推导过程，我们也可以得到NH3-N的平均增加值。

经分析，显然可得出以下结论：1)、高锰酸盐得主要来源为：湖南岳阳城陵矶及其上游地区，湖北宜昌南津关及其上游地区，重庆朱沱及其上游地区。2)、氨氮得主要来源为：湖南岳阳城陵矶及其上游，湖北宜昌南津关及其上游，重庆朱沱及其上游地区。即高锰酸盐和氨氮均主要来自于：湖南岳阳城陵矶及其上游地区。

3.3问题三

模型一一次累加拟合模型

观察已知数据可知各数据具有很大的随机性和波动性，所以我们采用了百分比的一次累加序列来进行多项式的拟合，并以拟合出来的多项式为依据，计算未来10年的各类水占评价河长的百分比：我们以2004年的评价河长为准，为固定值：全流域评价河长：39412km，干流评价河长：6341km，支流评价河长为：33071km。

首先，我们对数据序列按以下规则进行一次累加，得到一次累加数据序列

累加后的数据序列成为了一个单调递增序列，克服了原始数据序列的随机性和波动性。然后我们对一次累加数据序列进行多项式拟合，得出的表达式

然后将待预测的值带入，即得到在时刻的取值，根据以下式子即可得出待预测量的值：

由于从I到劣V类的水是单独预测的，故不能保证其总和为100％，因此我们应求得其第y年的和，用每类水所占评价河长百分比数据在第y年的和中所占百分比来修正：

再根据2004年的评价河长，即可得到各类水所占的河长。

模型二：基于时序分析法的水质总体趋势预测

将1995～2004年水文年全流域水质分类比例列为矩阵X，令：

得到新的矩阵如下:

预测方程式如下：，对数变换得。

令：，，则，，。

将上表有关数据带入计算得R＝1.9337，T＝1.8635，，。

预测2006年长江水域可饮用水得比例为

预测误差指标选择均方根误差RMSE；令预测误差：，则

预测值区间

式中：N－1为自由度，本次预测自由度为10－1＝9。为概率度，为避免与时间t混淆，因此将概率度由t换为。本次预测样本数未超过20项，为小样本，用t分布概率度，查t分布表，按80％信度计，自由度为9，t值为1.38，则：

基于时间序列法研究得出的结论可知，长江水质虽然有时某一年的情况较前一年可能会有所好转，但是它的总体趋势是在逐年恶化。其中好转的情况体现在随机误差较大这一特征之上。这与前面所提到的一次累加序列模型相统一的，同时也是对一次累加序列模型不能体现水质情况随机性的一个有力补充。

3.4问题四

本问要求我们计算如果未来10年内每年长江干流的IV类和V类水的比例控制在20%以内，且没有劣V类水,那么每年需要处理的污水量。即以控制污水总量的方式来控制IV类、V类、劣V类水的比例，因此我们需要知道污水总量与IV类、V类、劣V类水类水所占比例之间的关系。在此，我们建立了多元线性回归模型来求解。

从整体考虑，在过去的10年中，长江水的总流量，废水排放总量,IV类、V类、劣V类三种类别的水在水文年全流域的百分比之间的关系可以用以下式子来表示：

其中：

代入以上数据，关系式变为：

在2005～2014年间，我们每年都要处理一定量得污水，上式变为

变换得：

因为要处理得污水应尽可能得少，故建立以下目标函数

把以上数据代入规划模型，将模型化简后得：

目标函数：

目标函数及约束化简为：

线性规划求解得：1390.621亿吨。

将z输出整理即得出表中的数据，然后根据下式即可求出每年应处理得污水量。

3.5问题五

通过对附件一、二的认真研读，以及查阅相关的资料我们从政策层面方面提出以下两点建议：1、加强宣传力度唤醒民众正确认识长江的现状。2、立即制定相应的法律法规。加大对违法排污行为的处罚力度。在具体治理长江方面的几点意见：1、加强对中小支流的治理。2、加强主要污染城市的污水处理力度，污水达不到标准，则由政府出面买单。3、对污染严重的企业实行关停并转，集中建设实力强、能治理好污染的大企业。4、大力发展低污染农业，控制农药化肥的使用量。5、提高警惕。坚决防止江水治了又污。

4模型的评价

在本文中，我们没有使用单一的模型对问题进行分析求解，而是综合运用了多种数学模型，并且增加了一定的检验环节。这样就增加了我们所求得数据的合理性，同时也增加了论文的实际参考价值。在运用归一化法综合评价长江水质情况时，我们很巧妙的将多个指标转化为单一指标，这样就很好的解决了多指标难以比较的难题。在解决问题三时，通过一次累加序列拟合法和时间序列法的有机结合，解决了长江水质情况随机性和波动性较大从而难以准确预测的难题。

在改进方面我们认为：在问题一中模糊数学确定的污染指数可以改进为用模糊神经网络进行判断的指数。可以更好的反映污染程度。以及不同点的关系；问题二中，可以改进为灰色评价，以更好的反映污染物的来源；问题三中，我们用以预测的一次累加拟合模型虽然在数据点以内误差较小，但在数据点以外，进行预测时，误差将是不可知的，可以改进为使用时问序列法、灰色系统理论、BP神经网络进行预测，或者将其三者有机的结合在一起，通过优化组合形成最优化权值组合法，这样既可以达到较高的精度要求又可以避免应用大量的参数，达到最小误差预测的目的。