南水北调东线江苏段水质变化趋势分析.docxVIP

南水北调东线江苏段水质变化趋势分析 由于北京、天津、河北、山东和淮河流区的生态环境恶化，在这种背景下，实施了长江移民工程，缓解了北方水资源的短缺。南水北调东线江苏段南起长江, 北至微山湖, 输水主干线约617km, 设计总水头差约40.0 m, 沿途设13个梯级, 共22处枢纽、34座泵站, 逐级抽水入东平湖, 全程约300km。为确保北调水质达到Ⅲ类水标准, 2000-2012年期间, 东线江苏段治污工程项目共102项, 实际投资70.24亿元, 具体包括工业结构调整项目16项, 工业综合治理项目49项, 城镇污水处理及再生利用项目26项, 流域综合整治项目6项, 截污导流项目5项。本文通过应用季节性Kendall检验法, 对2001-2012年南水北调东线江苏段14个控制断面水质变化进行分析评价, 进一步了解南水北调东线江苏段通过治污工程水质的改善效果。 1 网络配置和分析项目 1.1 洪水断面位置。在转 南水北调东线江苏段输水通过三江营和高港2个引水口门从长江引水, 三江营为主要引水口, 高港在长江水位较低时补充引水。 南水北调东线江苏段输水干线上有6个水质监测断面, 由南往北依次为:江都西闸断面、戴楼断面、五叉河口断面、马陵翻水站断面、张楼断面和蔺家坝断面。除输水干线外, 南水北调东线支流水质控制断面有8个, 由南往北依次为:泰西断面、三垛西断面、老山乡、临淮乡、沙集西闸断面、单集桥断面、沙庄桥断面和李集桥断面。 1.2 水质指标分析 根据南水北调东线江苏段地表水污染特点, 确定溶解氧、高锰酸盐指数、五日生化需氧量、氨氮、石油类、挥发酚6项指标为分析项目, 采用江苏省环保局提供的2001-2012年水质监测资料对南水北调东线江苏段进行水质趋势分析。 2 双尾趋势检验t 季节性肯达尔检验是Mann-Kendall检验的一种推广。该检验的思路是用多年收集的数据, 分别计算各季节 (或月份) 的Mann-Kendall检验统计量S及其方差Var (S) , 再把各季节统计量相加, 计算总统计量。如果季节数和年数足够大, 可通过总统计量与标准正态表之间的比较来进行统计显著性趋势检验。 季节性肯达尔检验法原理是将每断面历年相同月或同一时期的水质资料进行比较, 如果后面的值 (时间上) 高于前面的值, 记为“+”号, 低于则记作“-”号。如果“+”的个数比“-”多, 则判为上升趋势;如果“-”个数比“+”多, 则划为下降趋势;如果“+”、“-”的个数分别为50%, 则判为无趋势。其适用于数据值有随季节性规律变化的多年趋势分析, 该方法分析计算公式如下:季节性肯达尔检验的零假设H0为随机变量, 与时间独立, 假定全年12个月的水质资料具有相同的概率分布。 设有n年p月的水质资料序列X为: 式中:X11, …, Xnp为月水质浓度监测值。 (1) 对于p月中第i月 (i≤p) 的情况 令第i月历年水质系列相比较 (后面的数与前面的数之差) 的正负号之和Si为 第i月内可以作比较的差值数据组个数mi为 式中:ni为第i月内水质系列中非漏测值个数。 在零假设下, 随机系列Si(i=1, 2, …, p) 近似地服从正态分布, 则Si的均值和方差如下: 均值:E (Si) =0 (2) 对于p月总体情况 在零假设下, p月的均值和方差为 式中:Si和Sk(i≠h) 都是独立随机变量的函数, 即Si=f (Xi) , Sh=f (Xh) , 其中Xi为i月历年的水质序列, Xh为h月历年的水质序列, 并且Xi∩Xh=ф。 因为Xi和Xh分别来自i月和h月的水质资料, 并且总体时间序列X的所有元素是独立的, 故协方差Cov (Si, Sh) =0。将其代入式 (4) , 得 当n年水质系列有y个数相同时, 同样有 当n≥10时, S服从正态分布, 且标准方差Z为 (3) 趋势检验 Kendall检验计量t定义为:, 由此在双尾趋势检验中, 如果|z|≤zα/2, 则接受零假设。这里FN(Zα/2=α/2) , FN为标准正态分布函数, 即: α为趋势检验的显著水平, 其值为: 取显著性水平α为0.1和0.01, 即当α≤0.01时, 说明检验具有高度显著性水平, 当0.01α≤0.1时, 说明检验是显著的。 在α计算结果满足上述两条件情况下, 当t为正时, 则说明具有显著 (或高度显著性) 上升趋势;当t为负时, 则说明具有显著 (或高度显著性) 下降趋势;当t为零时, 则无趋势。 3 水质变化 经过计算, 南水北调东线江苏段各控断面水质变化趋势分析结果见表1。 3.1 排水系统污染指标cod5 (1) 溶解氧 (DO) 。溶解氧在14个断面中, 12个断面呈上升趋势, 占85.7%;2个断面呈下降趋势, 占14.3%, 上升趋势断面多于下降