1. 引言
我国是一个水资源短缺的国家,水资源问题已成为制约社会经济发展的主要因素。就水污染现状而言,工业废物、农村面源污染和城市生活垃圾等成为水污染主要来源 [1] [2]。当未经处理或处理不当的工业废水和生活污水排入水中,其数量超过水体自净能力时,就会造成水体污染,对人体健康产生影响 [3] [4]。
近年来,我国的各大河流污染现象屡见不鲜。据报道,由于水土流失及矿产资源开发,鄱阳湖重金属潜在生态风险较高 [5];黄河兰州市区段河流每年都汇入数千万吨经处理或未经处理后的工业、生活废水 [6];淮河高锰酸盐指数浓度从1991年到1996年间呈明显上升趋势,1997年至2001年指标浓度呈下降趋势,2002年该指标浓度开始缓步上升 [7] 等。
河流在宿州市发展中扮演了重要角色,包括供水、城市景观等。因此吸引了很多学者的关注,一系列研究得以展开,如李奇等 [8] 利用内梅罗综合指数评价法对水质污染状况进行了评估,得出宿州各河流水体均受到不同程度污染,但较2004年相比,沱河芦岭桥断面有明显改善。黄淑玲等 [9] 通过对水文特征、污染物特征、水环境现状调查,同时结合水环境容量计算方法分析得出目前宿州市沱河水质良好,大部分河段存在剩余水环境容量。季定民等 [10] 通过采用单因子评价,基于直观图解、线性趋势分析和季节性Kendall检验法对水质污染物主要指标进行了分析,得出宿州市沱河段水质状况处于不断改善之中。
本文通过对沱河宿州段进行调查采样,分析河水中NH₃-N和TP的含量及其空间变化特征,并结合空间自相关分析和数理统计分析,以获取沱河宿州段的水质状况以及与城市影响有关的信息,从而为水资源保护提供参考。
2. 方法与材料
2.1. 研究区简介
宿州市位于安徽省北部,东与江苏省宿迁接壤,西与河南省商丘毗邻,南与安徽省蚌埠相连。宿州市辖区内有多条河流,本文研究对象为位于宿州市东部的沱河。沱河源自河南商丘附近,由永城县濉溪流入宿州市,西北至东南走向贯穿宿州市,全长275.13 km。宿州市沱河段是南沱河的上段,起始于沱河进水闸,终止于外环南路,全长12.7 km,水域面积44.7 km2,具有巨大的生态价值和社会价值。沱河流域历年来时有污染,政府斥巨资防治,已经取得了一定的成效。
2.2. 采样与测试
通过对沱河宿州段进行实地考察,确定其主要流经方向(西北流向东南)及输入口、输出口,并沿着沱河水流相反方向进行采样。样品采集在2018年4月完成,本着每隔50 m采集一个样品的原则,共采集了40个样品(具体分布见图1)。样品采集完成后,将其存放于干净、密封的容器中,并给它们标上序号,送进实验室进行处理。为了防止样品中存在杂质对仪器及测试产生影响,所有样品均在室温下进行抽滤(0.22 μm膜),除去杂物后放入仪器中测试。
2.3. 数据分析
本文着重研究了沱河水流中的NH3-N、TP含量。研究过程中对采集数据的进行统计分析,包括最小值、最大值、平均值、标准差、变异系数及正态分布检验p值。其中NH3-N、TP含量的空间分布利用Excel绘制,其空间分析用Geoda完成。
3. 结果与讨论
3.1. 含量特征
国标(GB 3838-2002)根据NH3-N、TP含量将地表水分为五类见表1:I类主要适用于水源地及国家自然保护区;II、III类主要适用于集中式生活饮用水;IV类主要适用于一般工业用水区及人体非直接接触的娱乐用水区;V类主要适用于农业用水区及一般景观要求用水。
本文研究样品的NH3-N、TP含量的基本统计结果见表2。从表2中可以看出,NH3-N含量介于0.057和0.661 mg/L之间,平均值为0.199 mg/L,而TP的含量介于0.013和0.064 mg/L之间,平均值为0.028 mg/L,与国家地表水环境质量标准相比(表1),从NH3-N含量来看,有39个样品具有二类及更好水质,仅有一个为三类水质,而从TP含量来看,40个样品都具有二类及更好水质。
变异系数(CV)是衡量资料中各观测值变异程度的一个统计量,计算方式为标准差与平均数的比值。在环境科学研究中,变异系数通常用于分析污染物含量是否受到明显的人为影响。研究表明,变异系数的大小代表着人为影响程度的强弱,当CV < 0.20时,表示受人为影响较低,代表低程度的人为贡献率;相反,当CV > 0.80时,表示人为影响较严重,代表着较高的人为贡献率 [11]。本文研究样品中NH3-N含量相对TP含量具有更高的CV值,分别为0.682和0.538 (表2),说明NH₃-N受人为影响较TP明显,这与上文根据水质分类推测的结果是相一致的。
正态分布检验p值可以用于分析数据是否符合正态分布。当p ≥ 0.05时,可以认为符合正态分布,间接表明河流污染受人为影响较低。反之,当p < 0.05时为非正态分布,河流污染受人为影响严重 [12]。宿州市沱河NH3-N和TP的p值均 < 0.01 (表2),说明可能受人为因素的影响,与上文推测的结果一致。
Table 1. Surface water environmental quality standards (GB 3838-2002) (mg/L)
表1. 地表水环境质量标准(GB 3838-2002) (mg/L)
Table 2. Statistical results of NH3-N and TP contents in Tuohe River (mg/L)
表2. 沱河水NH3-N、TP含量统计结果(mg/L)
3.2. 空间分布特征
污染物的空间分布可以为探讨污染物的来源以及人为影响提供直观信息 [13] [14]。因为本次研究采样的河流沿西北流向东南,采集路线与河流走向一致,因此本文采用了Excel软件绘制了沱河水中NH3-N和TP含量的线性分布图并示于图2,从中可以直观的观察出NH3-N和TP含量的空间分布特征。
从图2可以看出,沱河水中NH₃-N和TP含量的线状分布图具有相似的特征,其高含量区域主要集中于沱河下游。这种情况表明,宿州市沱河地表水中NH3-N和TP含量虽然污染程度上存在差异,但二者在来源上可能有相似的地方。这种推测与上文根据变异系数及正态检验p值得到的结果是相一致的。
比较沱河水NH3-N和TP含量的线性分布图能发现,NH3-N含量高的区域相对分散成两个区域,而TP含量高的区域则连成一片,这与区域的实际情况有关:对于沱河下游而言,其附近为田地及芦苇丛,可能有农作物肥料以及动植物腐烂的分解产物随水土的流失进入水体使NH3-N和TP含量增加 [15]。此外,对于沱河上游来说,附近出现排污口及大量的工业废屑,NH3-N含量增加可能与生活污水和某些含氮工业废水的排放有关 [16]。
Figure 2. Linear distribution of NH3-N (Upper figure) and TP (Figure below) contents (mg/L) in Tuohe River
图2. 沱河水中NH3-N (上图)和TP (下图)含量(mg/L)的线性分布
3.3. 空间自相关分析
空间自相关分析是对某一地理变量空间分布相邻位置的相关性进行检验的一种统计方法,它是通过检测一个位置上的变异是否依赖于邻近位置上的变异来判断该变异是否存在空间自相关性 [17] [18]。空间自相关在1948年由Moran等 [19] 提出。目前,已经被广泛应用于多个领域的研究中 [20] [21] [22] [23] [24]。空间自相关分析以相应统计量来表征地理变量的关联程度和集聚模式,统计量的正负表示空间自相关性的正负,统计量越大,表示空间自相关性越强 [25]。
按照Moran’I指数的分类,采集的样品可以划分为四个类别,分别为高–高、低–低、低–高和高–低,可以表示为某个样品与周边样品之间的含量上的相对关系 [25]。如高–高表示高含量的样品且其周边样品含量也很高,因此,这种样品点可以被称为热点;相反,低–低表示低含量的样品且其周边样品含量也很低,这种样品点被称为冰点;而低–高和高–低则表示异常点,可能与其他外在的因素影响有关。分析结果如图3所示。
Figure 3. The results of spatial autocorrelation analyses of NH3-N and TP contents (mg/L)
图3. NH3-N和TP含量空间自相关分析结果(mg/L)
1) 热点:从图3中可以看出,TP含量的高–高聚类样品主要出现沱河下游。TP含量的高–高聚类样品数为7,NH₃-N含量的高–高聚类样品数为0。从其聚集的情况来看,这种热点聚集的地方可以被称为热区,即沱河下游热区,可能与磷的面源污染有关。
2) 冰点:从图3中可以看出,NH3-N和TP含量的低–低聚类样品所出现的位置不完全相同,前者主要分布于图幅中间位置,而后者则主要分布于图幅中部及以上位置,对应的样品数均为3,且前者相对集中,后者相对分散。因为其含量偏低,所以低–低聚类的样品应该反映了自然环境的影响,而不是污染所致。
3) 异常聚类:这里包括低–高和高–低聚类的样品。从样品数及分布来看,NH3-N和TP的异常聚类样品数均为2个。分布比较分散,说明可能与点源污染有关。
3.4. 数理统计分析
相关研究表明,自然界中绝大多数物质的含量服从正态分布规律,如自然界中元素的背景值或者单一来源的污染物含量等。对于城市地表水而言,在没有受到外来因素的影响下,地表水中的元素含量也应该符合正态分布规律。换言之,如果地表水受到了一系列因素的扰动,就可能出现偏态的情况 [26]。
从NH3-N和TP含量的密度图来看(图4),均表现出双峰式的特征。因此,可以考虑NH3-N和TP在未受其他影响时是正态分布的(低含量部分),而异常部分(高含量)说明样品受到其他外来因素(如人为活动)的影响。
为了验证上述分析,我们还利用QQ图进行了分析(图5)。QQ图是一种由标准正态分布的分位数为横坐标,样本值为纵坐标的散点图,若图中的散点近似在一条直线附近,则可以看成符合正态分布,相反若散点波动范围较大,无法在一条直线附近则说明不符合正态分布 [27]。观察图5可以发现,NH3-N、TP的QQ图均分为两段,都不是近似在一条直线附近,因此NH3-N和TP的样品来源区域可能受到其他外在因素的影响。
Figure 4. Density plots of NH3-N and TP contents (mg/L)
图4. NH3-N和TP含量的密度图(mg/L)
Figure 5. QQ plots of NH3-N and TP contents (mg/L)
图5. NH3-N和TP含量的QQ图(mg/L)
4. 结论
通过对沱河宿州段水样品中NH3-N和TP含量的分析,获得了以下认识:
1) 与国家地表水环境质量标准相比,从NH3-N含量来看,有39个样品具有二类及更好水质,仅有一个为三类水质,从TP含量来看,40个样品都具有二类及更好水质。
2) 空间分布特征表明,具有高含量TP的样品主要分布在沱河下游,可能与面源污染有关,与利用空间自相关分析得出的认识相一致。
3) NH3-N和TP的密度图成双峰式,并且其QQ图分为两段,说明不符合正态分布,可能受到外在人为活动的影响。
基金项目
宿州学院地质资源与地质工程重点学科2018年度青年人才培养开放课题资助。