长江水质的数学评价和预测.docVIP

长江水质的数学评价和预测 朱金 杨浩乐 摘要： 本文应用插值原理对长江水质现状进行综合评价, 分析了各地区水质的污染状况,并应用状态空间模型和改进的kalmam递归式对长江未来水质状况进行预测分析。得出：如果不及时加以治理，长江污染状况将会随着工业的发展和时间的推移迅速恶化；并根据计算结果提出解决长江水质的建议和意见。 关键词： 水质， 插值评判， 状态空间模型， Kalmam递归式 。 问题复述: 水是生命之源,是人类赖以生存的资源,保护水资源就是保护我们自己,是我们每个人应尽的责任.对于我国大江大河水资源的保护和治理应是重中之重。 长江是我国第一、世界第三大河流,是华夏5000年文明的发祥地。如今,长江水质污染程度日趋严重已经引起相关政府部门和专家们的高度重视。 根据国家环保局的政府网站〔/apps/〕2003年6月至2005年9月发布的长江沿线8个省、17个水文监测站的水质报告,长江水利委员会编辑出版的《长江年鉴》1995年至2004年长江流域水质报告表以及2004年4月至2005年4月发布的长江干流主要观测站点的基本数据表明,原本“一江清水向东流”的母亲河正在悄悄被污染,而且污染指数在急剧攀升。 2004年10月10日,全国政协与中国发展研究院联合掀起了“保护长江万里行”活动。11月13日,南方新闻中心〔〕刊登了《专家称若不及时拯救,长江生态10年濒临崩溃》的文章。12月14日,北京青年报发出了《拿什么拯救 癌变长江》呼声。 根据2002年6月国家重新修订颁布的《地表水环境质量标准》,长江水质中Ⅳ类,Ⅴ类,劣Ⅴ类的比例正在逐年增加。新标准采用单因子极值评价方式似乎过于严格。所谓“极值法”是指以其中最差的指标来评价水质等级,也就是说水质标准当中只要1项达不到Ⅴ类,那么整体水质就要被判为劣Ⅴ类。这种“一票否决”式的考评方式有一定的片面性,对水质的考评应该综合各项指标后得出结论,现在只考虑其中1项指标不是太科学。为此,应用插值原理,并采用国家水利部门提供的最新监测数据〔附件3〕,对长江水质进行多因子线性综合评价。 对于长江未来水质的预测,根据资料〔附件4〕中数据,通过分析实验,长江水质污染状态具有很强的随机性,因此这种状态在时间状态的变化下具有“无后效性”的特点。为此,本文采用了状态空间转移模型和改进了的Kalmam递归式进行了较为合理、准确的预测分析。 问题〔1〕:对长江近2年多的水质情况做出定量的综合评价,并分析各地区水质的污染状况。 1 插值评判原理 1.1定义:已知 n＋1,j=0,1…,n, 其中互不相同,不妨设,求任一插值点x*〔≠〕处的插值y*.可以看成是由某个函数y=g〔x〕产生的,g的解析表达式可能十分复杂,或不存在封闭形式,也可以未知.㈠ 1.2 求解的基本思路是：构造一个相对简单的函数y=f﹙x﹚,f通过全部节点,即f﹙﹚=﹙j=0,1,…,n﹚,f﹙x﹚y*=f﹙x*﹚.2 水质综合评价的SP插值模型 设根据水质评价标准产生的某次水样的标准水质等级及其水质指标分别为y﹙i﹚﹛x*﹙j﹚|j=1～m﹜,=1～nn、m分别为样本容量和水质指标数目。污染越严重,水质等级就越高,最低水质等级设为1、最高水质等级设为N。为消除各水质指标的量纲效应,使建模具有一般性,对水质指标进行标准化处理 = (1) 式中:Ex(j)、Sx(j)分别为原第j个水质标准{x*(i.j)|I=1～n}{x(n＋1,j)|j=1～m}n个样本内插研究水体的水质等级yc(n＋1), 达到最小。式(2)中: (2) 式中:di为第i个样本的水质指标与研究水体水质指标之间的距离:wi为权重,表示第i个样本对内插研究水体的水质等级yc(n＋1)b为待定参数,一般为大于1的常数,b取得越大,则在点{x(n＋1,j)|j=1～mx(n＋1,j)|j=1～m (3) 这就是所求的对应点{x(n＋1,j)|j=1～m(1)根据水质评价标准随机生成水质等级样本系列x(i.j)及y(i),i＝1～n,j＝1～mb进行优化估计。在样本系列中任取某样本I,由其它n-1个样本进行Shepard插值,得到相应于水质等级y(i)的插值记yc(i)。可通过求解如下优化问题来优化估计参数b (4) s.t 1≤b≤5 式(6)的区间是根据笔者的经验确定的。这是一个一维非线性优化问题,模拟生物进化过程中优胜劣汰规则与群体内部染色体信息交换机制的加速遗传算法(Accelerating Genetic Algorithm,简称AGA)[7]是一种通用的全局性优化方法,用它来求解该问题显得十分简便而有效。AGA的详细算法可参考文献[7]。(3)