《长江水质评价和预测的数学模型》由会员分享,可在线阅读,更多相关《长江水质评价和预测的数学模型(14页珍藏版)》请在金锄头文库上搜索。
1、- 1 -长江水质评价和预测的数学模型阮门富、沈晓燕、郑丽心摘 要本文以长江流域的水质污染状况为背景,首先通过对 个观测站点近两年多的水17质状况有关数据的统计,应用概率统计分析方法,对长江水质状况作出了定量综合评价,并分析各地区水质污染情况,发现南昌水质污染较为严重。其次,针对长江干流 7 个观测站近一年多的基本数据,引用了河流污染物中非守恒物质的净化过程满足的一级反应动力学规律,建立河流中污染物一维稳态衰减规律的微分方程模型,求解出了这一年多长江干流上 7 个观测站的排污量,分析得出高锰酸盐指数主要的污染源是在湖南岳阳、湖北宜昌、重庆朱沱、江苏南京等地,每天的高锰酸盐指数排放量分别为 39
2、74.6 吨、3047.4 吨、 2808.9 吨、2713.3 吨;氨氮排放的主要污染源为湖南岳阳、湖北宜昌、重庆朱沱、江西九江等地,每天的氨氮排放量分别为 384.6309 吨、275.7372 吨、243.8681 吨、221.1189 吨。再次,利用过去 10 年的主要统计数据,运用基于灰色系统的灰色预测方法,预测出未来 10 年长江的水质污染情况。并通过预测到的未来 10 年内有关长江干流水质情况,在要求 IV 类和 V 类水的比例控制在 以内,且没有劣 V 类水前提下,建立以未%20来 10 年的处理的污水总量最小为目标的规划模型,通过求解模型可得未来 10 年每年需要处理的污水量
3、分别为 43.4 吨、56.2 吨、60.28 吨、70.75 吨、。最后,对长江水质污染状况给出了解决长江水质污染问题的一些可行的建议和意见。关键词:水质类别;灰色预测;水质污染 - 2 -1 问题的提出水是人类赖以生存的资源,保护水资源就是保护我们自己,对于我国大江大河水资源的保护和治理是重中之重。 长江是我国第一、世界第三大的河流,长江的水质污染程度却日趋严重,已引起了相关政府部门、专家们的高度重视。为了揭示长江的污染情况,引起人们的注意,专家们对长江沿线 17 个观测站(地区)近两年多主要水质指标进行了考察并给出了相关的检测数据、19952004 年长江流域水质的主要统计数据,以及干流
4、上 7 个观测站近一年多的基本数据(站点距离、水流量和水流速) ,和国标(GB3838-2002) 给出的地表水环境质量标准中 4 个主要项目标准限值(其中、类为可饮用水) , 现 要求对长江近两年多的水质情况做出定量的综合评价,并分析各地区水质的污染状况;进一步研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的主要污染源地区;假如不采取更有效的治理措施,依照过去 10 年的主要统计数据,对长江未来水质污染的发展趋势做出预测分析,比如研究未来 10 年的情况;根据预测分析,如果未来 10 年内每年都要求长江干流的类和类水的比例控制在 20%以内,且没有劣类水,那么每年需要处理的污水量 ;
5、提出自己对解决长江水质污染问题的切实可行的建议和意见。 2 模型的假设2.1 长江干流的自然净化能力是近似均匀的。2.2 影响水质类别的 4 个主要项目标准限值中取其最差水质类别为该样本水质类别。2.3 主要污染物高锰酸盐指数和氨氮的降解系数取 0.2(单位:1/天) 。2.4 水质的类别主要以溶解氧、高锰酸盐指数和氨氮这三个项目为衡量标准。3 符号的约定:第 个观测站的污染物的浓度(单位: ) ;iCLmg/:主要污染物高锰酸盐指数和氨氮的降解系数(单位:1/天) ;K:第 个观测站的河流断面平均水流速(单位: ) ;iu s/:第 个观测站与第 个观测站之间的距离(单位: ) ;ix1i
6、k: 第 个观测站的水流量(单位: ) ;il sm/3:第 个观测站(地区)的污染物总量(单位: ) ;iMg:第 个观测站的污染物总量中来自第 个观测站(地区)的污染物的量(单位:i 1i) ;mg:第 个观测站(地区)的排污量(单位: ) ;i mg4 模型的建立4.1 问题一411 根据题意的分析可知,I、II、III 类为可饮用水,而其中 III 类水属于被轻- 3 -度污染的水质,IV、V、劣 V 类水都是受到较严重污染的水且污染程度逐类增加,按污染程度的不同把各类水分三大类对干流的 7 个观测站点和支流 10 个观测站点出现三大类水的情况进行统计,例如:2003 年 6 月出现
7、IV,V,劣 V 类的站点是 0 个,III 类水的站点为 1 个,而 I,II 类是 6 个站点。通过对给出的附件 3 进行统计得到下面的表1(长江 7 个干流点,分别依次为:四川攀枝花、重庆朱沱、湖北宜昌南津关、湖南岳阳城陵矶、江西九江河西水厂、安徽安庆皖河口、江苏南京林山,在 28 个月中每个月的水质类型分布情况)和表 2(长江 10 个支流点,依次分别为:四川乐山岷江大桥,四川宜宾凉姜沟,四川泸州沱江二桥,湖北丹江口胡家岭,湖南长沙新港,湖南岳阳岳阳楼,湖北武汉宗关,江西南昌滁槎,江西九石,江苏扬州三江营,在 28 个月中每个月的水质类型分布情况)2003,6 2003,7 2003,
8、8 2003,9 2003,10 2003,11 2003,12 2004,1 2004,2 2004,3IV,V,劣 V0 0 0 0 0 0 0 0 1 0III 1 2 0 1 1 0 1 1 2 0I,II 6 5 7 6 6 7 6 6 4 72004,4 2004,5 2004,6 2004,7 2004,8 2004,9 2004,10 2004,11 2004,12 2005,1IV,V,劣 V0 0 0 0 0 1 0 0 0 0III 0 1 0 2 1 1 0 0 1 2I,II 7 6 7 5 6 5 7 7 6 52005,2 2005,3 2005,4 2005,5
9、 2005,6 2005,7 2005,8 2005,9IV,V,劣 V0 0 0 0 0 0 0 0III 1 2 1 1 1 0 2 1I,II 6 5 6 6 6 7 5 6表 12003,6 2003,7 2003,8 2003,9 2003,10 2003,11 2003,12 2004,1 2004,2 2004,3IV,V,劣V2 5 3 2 2 3 5 5 5 5III 2 1 1 3 2 4 2 1 1 2I,II 6 4 6 5 6 3 3 4 4 32004,4 2004,5 2004,6 2004,7 2004,8 2004,9 2004,10 2004,11 2004
10、,12 2005,1IV,V 3 4 2 3 1 1 2 2 3 2月份分类月份分 类- 4 -,劣VIII 2 2 3 4 4 2 2 2 3 3I,II 5 4 5 3 5 7 6 6 4 52005,2 2005,3 2005,4 2005,5 2005,6 2005,7 2005,8 2005,9IV,V,劣V4 1 1 1 2 1 1 1III 0 3 3 4 1 3 3 2I,II 6 6 6 5 7 6 6 7表 2从表 1 和表 2 可以看出:长江干流的水污染不严重,只有两个月份存在 IV,V,劣 V 类污染水。III 类轻度污染的水所占比重也比较小。长江支流的水有一定程度的污
11、染,28 个月中,每个月至少有一个支流的观测站点存在 IV 类、 V 类和劣 V 类水,最多的时候有 5 个支流的观测站点存在 IV 类、V 类和劣 V 类水。且枯水期和平水期比丰水期存在 IV 类、V 类和劣 V 类水的观测站点多一些。下面对长江干流 7 个观测站的水质的 4 个项目指标进行分析。PH 值为无量纲量在此不讨论。对于高锰酸盐指数(CODMn)、氨氮(NH3-N) 、溶解氧(DO)中的每一项监测项目,把七个观测站点的每月每一项目的数据进行累加,然后求出每月每一项目指标的平均值,再求出对应的类别(I 类,II 类,III 类,IV 类,V 类,劣 V 类) 。每月的水质类别对应有
12、3 个值,根据附件 3 的水质分类标准,我们取各项目中水质类别最差的那一个作为长江水的水质类别。图象如图 1 所示: (图 1)从图 1 可以看出,长江干流大部分时间里的水质都为 II 类,有个别月份的水是 III类水。第 3 个月的水由第 2 个月的 II 类水变成 III 类水; 第 8 个月水也由第 7 个月的II 类水变成 III 类水,一直持续到第 9 个月;第 19 个月水则由第 18 个月的 II 类水变成III 类水,持续的时间有 6 个月之长,即一直持续到第 24 个月。由此可见,水质的净化时间变长了。综上可得,长江干流的水污染不严重(这主要是因为长江水量目前还比较大,掩盖了
13、长江水质污染问题的重要性) ,而长江支流的水则有一定程度的污染,长江干流的自然净化能力越来越差。- 5 -412 下面对各地区的水质的污染状况进行分析:图 2-4 为 17 个观测站的平均的水质类别和各站各自水质最差、最好的水质类别。用均值的方式求出各个点 3 个项目 28 个月的水质各自对应标准值的均值,求出对应的类别。我们取各项目中水质类别最差的那一个作为长江水的水质类别。跟上面类似,可以求出长江 17 个观测站的水质类别,如图 2-4 所示:图 2-4表 3(如下)为长江在近两年多共 28 个月,17 个观测站每个观测站的水质分布情况:1 2 3 4 5 6 7 8 9 10 11 12
14、 13 14 15 16 17IV,V,劣 V2 0 0 0 0 0 0 15 4 8 0 11 3 1 25 4 1III 3 8 3 10 0 1 1 9 4 4 0 14 10 9 3 5 7I,II 23 20 25 18 28 27 27 4 20 16 28 3 15 18 0 19 20表 3综合表 3 和图 2-4 可以得到各地区水质的污染状况:1没有污染(水质类别为 I,II 类)的观测站:江西九江河西水厂、湖北丹江口胡家岭轻度污染(水质类别为 III 类)的观测站:重庆朱沱、湖北宜昌南津关、湖南岳阳城陵矶、安徽安庆皖河口、江苏南京林山2污染较严重(水质类别为 IV,V,劣
15、V 类)的观测站:四川攀枝花、四川乐山岷江大桥、四川宜宾凉姜沟、四川泸州沱江二桥、湖南长沙新港、湖南岳阳岳阳楼、湖站点分类- 6 -北武汉宗关、江西南昌滁槎、江西九石、江苏扬州三江营(其中水质类别为劣V 类的观测站为:四川泸州沱江二桥、江西南昌滁槎)其中水质最差的观测站为江西南昌滁槎,因为江西南昌滁槎在 28 个月中水质类别最好的那一个月都有轻度污染。 4.2 问题二一个观测站(地区)的水质污染主要来自于本地区的排污和上游的污水。即:本地污染物总量本地区排污量来自上游的污染物于是得到:本地区排污量本地污染物总量-来自上游的污染物本地污染物总量我们已经通过附录 3 得到,对于来自上游的污染物,我
16、们有下面的关系:进入环境中的污染物可以分为两大类:守恒物质和非守恒物质。因为江河自身对污染物有一定的自然净化能力,即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低,所以河流污染物高锰酸盐指数和氨氮属于非守恒物质,从而河流污染物高锰酸盐指数和氨氮的净化过程满足一级反应动力学规律2,即 KCdt积分,得 te0其中, 为污染物浓度, 为上一观测站的污染物浓度, 为反应时间。C0Ct在流速为 ,距离为 的两个观测站之间,有 ,于是uxuxtKe0记长江干流上的七个观测站的污染物浓度的观测值为 ,( ,2,,7 分别代表iC1四川攀枝花、重庆朱沱、湖北宜昌南津关、湖南岳阳城陵矶、江西九江河西水厂、安徽安庆皖河口、江苏南京林山)。则这七个观测站中第 个观测站污染物的总量(溶质i的量)为 iM( ,2,,7)TlCii1i来自上一观测站(即第 个观测