城市生活污水处理过程三维荧光光谱在线监测分析方法

杨金强， 赵南京， 殷高方， 俞志敏， 甘婷婷，王 翔， 陈 敏， 冯 春

1. 中国科学院环境光学与技术重点实验室， 中国科学院安徽光学精密机械研究所， 安徽 合肥 230031 2. 合肥学院生物与环境工程系， 安徽 合肥 230601 3. 中国科学技术大学， 安徽 合肥 230026 4. 安徽省环境光学监测技术重点实验室， 安徽 合肥 230031

引 言

城市生活污水在组分、 含量、 结构的多样性和复杂性， 致使污水处理系统具有时变性、 时滞性、 扰动性以及非线性等特征[1]。 污水处理过程水质的传统检测方法是基于水体COD、 TN等指标分析， 得到的数据往往不能真实反应污水中各污染物成分与含量， 难以实现实时监测， 还会产生二次污染， 导致无法及时准确地反馈各处理流程物料平衡和控制效果[2]。 实时有效地监测生活污水水质状况， 是控制污水处理过程物料投送和处理效果评价的重要依据。

相关研究表明[3], 三维荧光光谱可以很好应用于水体溶解性有机物(DOM)的检测和识别， 寻峰法是实现三维荧光光谱分析常用方法[4]， 其基于荧光基团特定的激发/发射波长的荧光峰， 通过获取特定位置荧光信息对水体荧光基团进行表征。 对于组成较为复杂水体， 不同物质成分信息常常存在着较高的相关性， 使得光谱信息存在一定程度重叠， 导致不同荧光峰识别不够准确， 仅能选取部分特定荧光峰用于分析， 还会造成大量光谱信息冗余。 针对此问题[5]， 论文研究采用三维荧光光谱结合主成分分析法对城市生活污水光谱特征进行分析， 基于特定荧光基团位置的不同， 对荧光光谱进行分区， 求取各区域第一主成分区域值， 并利用其值表征污水中DOM的含量变化特征， 研究建立其与常规水质参数相关性， 实现城市生活污水处理过程水质的在线实时监测。

1 实验部分

1.1 样品

实验分析水样取自合肥某生活污水处理厂， 处理工艺为A2/O(Anaerobic-Anoxic-Oxi)。 采样点选取为： 进水管道口、 厌氧池出水口、 缺氧池出水口、 好氧池出水口、 厂区出水口。 水样采集时间： 2018年9月至2019年3月， 平均每周1次， 共采取进水管道口、 厌氧池出水口、 缺氧池出水口、 好氧池出水口、 厂区出水口各25个样本， 合计总样本数为125个。 现场使用聚氯乙烯瓶(PVC)收集原水样， 取回后在实验室沉淀1～2 h后采用0.45 μm醋酸纤维滤膜过滤， 去除难溶性杂质， 处理后的水样在5 h以内分别完成COD、 TN、 三维荧光光谱测量。

实验采用重铬酸钾指标法分析测定COD， 过硫酸钾氧化—紫外分光光度法测定TN， 三维荧光光谱采用荧光分光光度计(HITACHI F-7000)进行测量， 其参数设置为： 发射和激发波长扫描范围200～500 nm， 扫描狭缝宽度5 nm， 扫描速度12 000 nm·min-1， 扫描间隔10 nm， PMT电压为600 V， 空白样为超纯水。

1.2 方法

对三维荧光矩阵数据Eij进行主成分分析[6]，Eij为一个61维向量组， 假定二阶矩阵存在， 求取其协方差矩阵

cij=cov(XiYj)=E{[Xi-E(Xi)][Xj-E(Xj)]}

(1)

经线性变换得到全新变量Y1，Y2， …，Yn(i为激发波长选取点，j为发射波长上选取点)。

(2)

图1 污水处理各流程三维荧光光谱

满足以下两个条件

(3)

其中第i主成分Yi=tixi, 方差λi，i=1，2，3，…，n， 其中第一主成分为Y1， 其特征向量系数为A1j=(a11，a12，…，a1j)。 为验证方法的有效性引入主成分贡献率λ， 见式(4)

(4)

2 结果与讨论

2.1 光谱分区

污水处理各流程三维荧光光谱如图1(a)— 图1(e)所示。

从图中可以看出， 荧光光谱中主要包含T1， T2， FA和HA四个荧光峰。 结合文献可知[7-8]， 四个荧光峰依次为芳香性蛋白类物质T1、 微生物代谢产物T2、 腐殖酸类物质FA和富里酸类物质HA。 随处理流程各荧光峰强度不断变化， 其中荧光峰T1与T2在整个流程中降低明显， 好氧工段T1降幅最大， 无氧工段T2降幅最大； 荧光峰FA和HA在各工段均有降低， 但降低幅度较小。

根据各荧光峰位置， 对激发波长(EX)和发射波长(EM)所形成的荧光区域进行分区， 如图1中(f)所示。 将荧光光谱分成区域A(Ex： 250～300 nm, Em： 280～380 nm)、 区域B(Ex： 200～250 nm, Em： 280～380 nm)、 区域C(Ex： 250～500 nm, Em： 380～500 nm)和区域D(Ex： 200～250 nm, Em： 380～500 nm)四个区域， 各区域代表性物质依次为： 微生物代谢产物、 芳香性蛋白类物质、 类腐殖酸和富里酸类物质。

2.2 指示性特征分析

研究分析了25次采样100个样本成分贡献率λ， 对于样本数据Eij，λ均超过85%， 因此可以采用第一主成分Y1描述样本所有信息。 利用水体各荧光峰所在区域的不同， 控制选取区域范围中特征矩阵向量系数Aij， 分别求出区域A， B， C和D第一主成分荧光区域值如图2所示。

图2 污水处理流程中各物质主成分区域值

由图2可知， 随着污水处理流程各物质第一主成分区域值不断变化， 其中芳香性蛋白类物质(区域B)第一主成分区域值减少量达86%， 在好氧工段降幅最大； 微生物代谢产物(区域A)第一主成分区域值减少量到达81.3%， 在无氧工段降幅最大； 腐殖酸类 (区域C)和富里酸类物质(区域D)第一主成分区域值在各流程较为稳定。 结合相关研究[9-11]： 对于芳香性蛋白类物质， 在污水处理各流程物质含量都会降低， 尤其在好氧工段， 由于好氧类细菌大量吸收小分子蛋白释放无机物导致其降低明显； 微生物代谢类产物在无氧工段为厌氧类细菌提供物质和能量进行硝化反应， 被大量分解吸收导致其含量下降明显； 城市生活污水中的腐殖酸和富里酸类物质主要来源于人类代谢产物和土壤， 在生物处理阶段含量较为稳定。 可以看出， 各区域主要污染物含量变化与各区域第一主成分区域值变化规律相符。 通过余弦相似度方法对不同流程各区域主成分区域值与各区域主荧光峰峰值作相似性分析， 结果如表1所示。

表1 第一主成分区域值与各区域荧光峰相似度分析结果

Table 1 The first principal component region value and the fluorescence peak similarity analysis results of each region

相似度芳香性蛋白类物质区域值/%微生物代谢产物区域值/%富里酸类物质区域值/%腐殖酸类物质区域值/%进水口98.4399.7293.2193.17厌氧池98.7698.9591.4294.17缺氧池97.9699.6393.3492.97好氧池98.7799.1192.9191.57出水口98.2199.2291.8991.83

从表1可以看出， 四个物质荧光峰所在区域第一主成分区域值∑iY1与荧光峰峰值具有很高的相似度。 因此， 可以利用第一主成分区域值对污水处理各流程物质进行分析， 其不仅包含了区域峰值点信息， 还包含了不同激发发射波长下物质信息， 避免了有效光谱信息浪费。

2.3 处理过程指标相关性分析

对于城市生活污水， COD和TN是评价污水处理效果和有机物去除能力的重要指标[12-14]， 代表了城市生活污水水体有机物污染程度。 由相关研究可知， 类蛋白类物质含量与COD具有较好的相关性， 而TN主要来源于溶解性有机物， 水体溶解性有机物主要由芳香性蛋白类物质、 微生物代谢产物、 腐殖酸类物质和富里酸类物质等构成。 因此， 研究选取芳香性蛋白类物质主成分区域值(S1)以及各物质主成分区域值比Yp/Yf(芳香性类蛋白和微生物代谢产物第一主成分区域值之和与腐殖酸类和富里酸类物质主成分区域值之和的比值)与水体COD及TN作相关性分析， 并对比了各流程中芳香性蛋白类物质的主成分区域值及各物质主成分区域比值(Yp/Yf)与水体COD及TN在各流程变化趋势。 结果如图3和图4所示。

由图3和图4可知， 伴随着污水不同处理流程， 利用主成分分析得到的荧光参数变化趋势与常规水质指标变化趋势具有很高的一致性。 芳香性蛋白类物质第一主成分区域值(S1)与COD呈多项式相关， 相关系数可以达到97.63%， 微生物代谢产物主成分区域值之和以及各物质主成分区域比值(Yp/Yf)与水体TN呈线性相关， 其相关系数可以达到94.02%。

图3 水质参数变化趋势对比图

图4 S1, Yp/Yf与COD及TN相关性

3 结 论

采用三维荧光光谱结合主成分分析方法对城市生活污水处理厂各处理流程水样进行了研究， 将污水三维荧光光谱分为芳香性蛋白类物质、 微生物代谢产物、 腐殖酸类和富里酸类物质四个区域， 采用主成分分析法对四个光谱区域进行分析， 研究了各荧光光谱区主成分区域值与物质变化规律， 通过与水质常规指标相关分析， 表明水体芳香性蛋白类物质第一主成分区域值与COD具有良好相关性， 相关系数r为97.63%， 各物质主成分区域比值Yp/Yf与TN相关系数r达94.02%。 因此， 三维荧光光谱结合主成分分析可为城市生活污水处理各工艺流程水质监测、 工艺优化以及处理效果评估提供快速在线分析方法。