長江水质的评价和预测模型-wzh.docVIP

长江水质的评价和预测模型 摘要 本文研究了近年来长江流域水质监测指标和主要的统计数据，建立了长江水质定量的综合评价模型、干流各地区污染物排放量的计算模型和未来十年水质的预测模型等三个数学模型。 根据各种污染物对水质等级评定的贡献不同，利用加权平均的方法，使水质评定的等级由整数和小数两部分组成；建立了长江水质定量的综合评价模型，并对整个长江流域两年多来水质的变化情况以及长江沿线各地区水质的污染状况作了分析、比较，得出的主要结论是：长江流域平均水质在冬季较差，在夏季较好；干流各地区水质较为平均，且比支流的水质要好。 根据污染物在水流中满足自然降解的微分方程和一个地区的水质污染主要来自于本地区的排污和上游的污水的关系，建立了长江干流沿岸各地区主要污染物排放量的计算模型，并找到了长江干流高锰酸盐、氨氮的主要污染源是湖南岳阳城陵矶湖北宜昌南津关i,j): 第j个断面（观测站）的第I项监测数据（i=1,2,3,4时分别对应 PH值， 溶解氧浓度，高锰酸盐指数浓度，氨氮浓度） x（i）：第i年废水排放总量 y（i）：第i年长江总流量 Z（i，j）：第i年干流中j类水所占的百分比 3. 模型的建立与求解 3.1 模型一：长江水质定量的综合评价模型 断面综合水质指数公式 （1） （2）公式1中r为系数，它的作用是使 表1 r 的经验值 水质类别 Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ 劣Ⅴ r的经验值 1.00 0.47 0.45 0.41 0.18 0.17 若上述r仍不能使 ，则的取值为0.999 表2 《地表水环境质量标准GB3838—2002）项目标准限值单位：mg/L 序 号 ? 分 类标准值 Ⅰ类 Ⅱ类 Ⅲ类 Ⅳ类 Ⅴ类 劣Ⅴ类 1 PH值（无量纲） 6---9 ??????????≥ 7.5 （或饱和率90%） 6 5 3 2 ≤ 2 4 6 10 15 ∞ 4 氨氮（NH3-N）????? ≤ 0.15 0.5 1.0 1.5 2.0 （3）公式1中为第种污染指数的权重，它的计算方法是： 设Ⅰ类水标准限值为，Ⅴ类水标准限值为。，的数值从表 2查出 令，取， 特别的，（因为）， 取 （4）公式1中： （5）各种污染物的权重计算结果如下： 表3 各种污染物的权重 污染物 PH值 溶解氧 高锰酸盐指数 氨氮 权重 0.0497 0.1450 0.2899 0.5154 （6）根据公式1，用Matlab计算17个主要观测站（地区）水质定量的综合评价。（计算程序见附录二的程序1） 3.2 模型二： 长江干流主要污染物排放量计算模型 （1）干流污染物自然降解方程 设污染物浓度，为降解系数， 得微分方程： ，初始条件： 当取时， 解得 （2） （2）补充假设： a） 设2004年4月至2005年5月分别编号为第1—13月 b） 为第个月水流从地区到地区所用时间 c） 为第个月，第个地区的排污量 d） 为第个月，第个地区的水流量 （3）根据2.1基本假设（2），相邻两个观测站的水流速取这两个站点的平均值，利用Matlab计算时间，计算结果见附录一的表12，计算程序见附录二的程序2 （4）排污量计算公式： ， （3） 特别的，我们认为位于最上游的四川攀枝花的排污量。 （5）根据公式3，用Matlab计算高锰酸盐和氨氮排放量，计算程序见附录二的程序3，计算结果见附录一的表10，表11 3.3 模型三 ： 未来十年长江水质发展趋势预测模型[3] （1）废水排放总量预测 图1 废水排放总量x（i）与年份i的关系 根据图1，建立如下的二次回归模型： （4） 用Matlab计算得2005-2014年废水排放量：（计算程序见附录二的程序4） 表 4 2005-2014年废水排放量（取整数） 年份 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 污水排放量（亿吨） 308 331 356 382 410 440 471 504 538 574 （2）长江总流量预测 图2 长江总流量y 对年份 i 的散点图 分析图2，我们发现总流量y有明显的周期性，即每四年出现一次流量极大值。（比如1998，2002都是流量的极大值年份。）按此规律，并使用1995-2004年的数据，我们