紫外光谱法检测COD时波长和低温影响的探究
李鑫1, 苏成志1,2,*, 于丹阳1, 盛宇博1, 常川1, 石磊1, 姜吉光1
1. 长春理工大学机电工程学院, 吉林 长春 130022
2. 长春理工大学人工智能研究所, 吉林 长春 130022
*通讯联系人 e-mail: Chengzhi_su@126.com

作者简介: 李 鑫, 1994年生, 长春理工大学机电工程学院硕士研究生 e-mail: cust_xin_li@126.com

摘要

COD代表了水体受还原性物质污染的程度。 相对于采用传统方法检测COD, 存在检测时间长且操作复杂等缺点, 紫外光谱法以其检测速度快, 无需使用化学试剂等特点成为了主流的检测方法。 基于朗伯-比尔定律, 以邻苯二甲酸氢钾粉末配置的COD标准溶液为研究对象, 针对低温环境下利用紫外光谱法检测COD精度的问题, 分别对COD的最佳检测波长和温度对COD检测值的影响进行研究。 同时选择长春市某地区地表水为研究对象, 验证COD最佳检测波长在实际水样中的适用性及温度补偿模型的准确度。 在研究检测波长对COD检测值的影响时, 选用256, 266, 276, 286和296 nm共5个波长对样本进行回归分析, 它们的吸光度分别为 A256, A266, A276, A286, A296, 将吸光度 A与COD标准溶液浓度值进行线性回归, 通过拟合结果得出: 276, 286和296 nm处模型具有代表性, 且在286 nm处拟合效果最好, 296 nm次之, 最后为276 nm, 其中286 nm处相关系数 r为0.994 6, 决定系数 R2为0.989 4, 波长为296 nm处和方差SSE=0.011 4, 预测均方根误差RMSE=0.037 7, 但其决定系数 R2较低, 可见在286 nm处COD检测值与吸光度具有最高的相关性, 又探究了标准温度(20 ℃)下8 mg·L-1的COD实际水样和标准水样的光谱吸收情况, 得出286 nm同样适用于实际水样的检测, 可见286 nm处为最佳检测波长。 在研究温度对COD检测值的影响时, 采集不同温度下COD实际水样与标准水样的紫外吸收光谱, 经过分析得出: COD实际水样中紫外光谱吸收度会随温度升高而增大。 为了减弱在COD测量中温度的影响, 根据最小二乘原则, 建立温度补偿模型。 利用实际水样验证温度补偿模型的准确度, 同时进行误差分析, 分析结果表明: COD的实际值与补偿后值的最大相对误差为6.38%, 最小相对误差为0.63%, 且多数相对误差集中在4%, 由此可见, COD温度补偿模型补偿精度高, 效果良好。 结果表明: COD检测选取的最佳波长及温度补偿模型可有效的提高COD低温检测精度。

关键词: 化学需氧量; 低温补偿模型; 紫外光谱法; 吸光度
中图分类号:O657.3 文献标志码:A
Study on the Influence of Wavelength and Low Temperature on COD Detection by Ultraviolet Spectroscopy
LI Xin1, SU Cheng-zhi1,2,*, YU Dan-yang1, SHENG Yu-bo1, CHANG Chuan1, SHI Lei1, JIANG Ji-guang1
1. College of Mechanical and Electric Engineering, Changchun University of Science and Technology, Changchun 130022, China
2. Institute of Artificial Intelligence, Changchun University of Science and Technology, Changchun 130022, China
*Corresponding author
Abstract

COD represents the degree of water pollution by reducing substances, compared with the traditional method to detect COD, the detection time is long and the operation is complicated. Ultraviolet spectroscopy has become a mainstream detection method due to its fast detection speed and no need for chemical reagents. Based on the Lambert Beer law, using potassium hydrogen phthalate powder to prepare standard solution as an experimental object, aiming at the detection accuracy of COD ultraviolet spectrum at low temperature, the optimal detection wavelength of COD and the influence of temperature on the detection value of COD were studied respectively. At the same time, the surface water in a certain area of Changchun City is selected as the research object to verify the applicability of the best COD detection wavelength in the actual water sample and the accuracy of the temperature compensation model. When studying the influence of detection wavelength on COD detection value, choose A256, A266, A276, A286, A296 five wavelengths to regression analysis samples, including A256, A266, A276, A286, A296 is A wavelength of 256, 266, 276, 286 and 296 nm absorbance, the absorbance of A linear regression with the COD standard solution, can be seen from the fitting data 276, 286 and 296 nm model representative, 286 nm in fitting out the best effect, 296 nm, Finally, it is 276 nm, of which the correlation coefficient r of 286 nm is 0.994 6 and the determination coefficient R2 is 0.989 4, SSE=0.011 4 and RMSE=0.037 7 at 296 nm, but the determinant coefficient R2 is low. It can be seen that 286 nm is the highest correlation and the smallest error. The results show that 286 nm is also suitable for the detection of actual water samples, and 286 nm is the best detection wavelength. When studying the influence of temperature on COD detection value, UV absorption spectra of COD water samples and standard water samples were collected at different temperatures. The results show that the UV absorption of COD solution increases with the increase of temperature. After thoroughly studying the spectral absorption of the actual and standard water samples with the same concentration of COD at the standard temperature (20 ℃), in order to eliminate the influence of temperature on COD measurement, a temperature compensation model was established by the least square method. The accuracy of the temperature compensation model is verified by the actual water sample, and the error analysis is carried out at the same time. The results show that the maximum relative error between the actual value of COD and the compensated value is 6.38%, the minimum relative error is 0.63%, and most of the relative errors are concentrated in 4%, which shows that the fitting effect of the model is good. Thus, COD temperature compensation model has high compensation accuracy and good effect. Finally, the conclusion is drawn that the best wavelength and temperature compensation model selected for COD detection can effectively improve the accuracy of COD low temperature detection.

Keyword: Chemical oxygen demand; Low temperature compensation model; Ultraviolet spectrometry; Absorbance
引 言

化学需氧量(chemical oxygen demand, COD)[1, 2]是指在一定条件下氧化水样中还原性有机物质所消耗的氧化剂的量。 在水质检测中, 还原性物质对水体污染的程度是重要检测指标之一[3]。 目前国内外检测COD的常用方法有重铬酸钾氧化法、 高锰酸盐指数法、 超声消解法、 微波消解法、 化学发光法和电化学方法等[4, 5], 但这些方法都存在需要取样到实验室进行测量, 操作要求水平高, 成本高且对环境污染严重, 测量范围窄等问题[6]。 相对于以上检测方法存在的缺点, 紫外光谱法以其操作快捷简便, 对环境无污染等优点实现了对水质原位的实时测量, 也因此成为了光谱法中最具代表性的方法[7, 8]。 在利用紫外光谱法检测COD时, 大部分学者都是仅对COD标准水样进行探究, 例如杨虹等用邻苯二甲酸氢钾配置的标准水样中的COD进行探究[9], 只针对COD标准水样进行探究往往不具有普遍应用价值。 因此如何利用紫外光谱法更加精确的在标准水样的基础上探究实际水样中的COD值成为目前各国学者探索问题之一。

波长的选择是影响紫外光谱法检测COD精度的一个重要因素, 当前很多学者都是对某一波段进行研究, 徐熠刚[10]等利用254 nm波长的单波长紫外光对水中COD检测进行探究, 杨云开等[11]为了消除水体浑浊度对紫外光吸收造成的影响在波长546 nm处检测COD。 但在实际水样中物质成分复杂, 部分无机物、 金属离子在254 nm处同样有吸光特性, 使测量结果不准确, 因此本文对不同波段的紫外光进行研究并选出最佳波长。

温度也是影响COD检测的重要因素之一[12], 大部分学者只是得到温度对COD检测存在影响的结论, 如周昆鹏等[13]探究了温度对COD检测的影响, 但并未建立温度补偿模型, 很难应用到实际检测中。 因此建立一种能将低温环境中COD检测值实时补偿到标准温度的温度补偿模型是很有必要的。

采用标准水样进行探究, 对不同浓度标准水样的吸光度进行回归分析, 通过拟合结果确定COD最佳检测波长, 并根据最小二乘原则, 在最佳检测波长下建立温度补偿模型, 并以实际水样为研究对象, 验证了最佳波长的适用性及温度补偿模型的准确性。

1 实验部分
1.1 仪器与样品

实验仪器是由厦门市奥谱天成科技有限公司生产的闪烁氙灯、 可调光程式光纤传感探头和全光谱光纤光谱仪、 中控器组成, 通过用光纤传感系统来测得溶液的吸收光谱, 该系统包括光源、 可调光程式光纤传感探头、 光谱仪、 光纤、 中控器、 显示器组成, 测试总波长范围为200~1 000 nm, 分辨率2.3 nm, 信噪比> 500:1, 系统如图1所示。

图1 实验仪器实物图Fig.1 Physical diagran of experimental divice

光谱扫描软件为系统自带软件。 系统将闪烁氙灯的光源通过光纤发射到可调光程式光纤传感探头的反射镜上, 全光谱光谱仪将反射镜反射回来的光进行数据采集, 得到COD吸光度数据传输到显示平台中, 从而获得全波段的溶液光谱图。

实验样品为依据国标规定的重铬酸钾法所配比的COD标准溶液与采集长春市某地区地表水的实际水样, 该实际水样浓度经山东青岛某环境研究院采用重铬酸钾国标法检测, 检测结果显示: COD的值为8 mg· L-1, 该方法检测上限700 mg· L-1, 检测下限4 mg· L-1, 相对标准偏差为1.3%~11%。

1.2 光谱的采集

为确定在紫外光谱法下检测COD的最佳波长, 配制浓度为20, 50, 125, 160, 250和350 mg· L-1的COD标准溶液, 并采集其在全波长下的吸收光谱图, 见图2。

图2 不同浓度COD标准水样光谱图Fig.2 COD spectrum of standard water sample with different concentration

为验证最佳波长在实际水样中的适用性, 分别对8 mg· L-1 COD实际水样和8 mg· L-1 COD标准水样进行光谱采集, 如图3所示。

图3 COD实际水样与标准水样(20 ℃)光谱图Fig.3 COD spectrum of actual water sample and standard water sample at 20 ℃

为了研究温度对使用紫外光谱法测量水中COD时测量值是否受到影响, 分别选取浓度为10, 20, 30, 40, 50, 60, 70, 80, 90和100 mg· L-1的COD标准水样, 在温度为0, 5, 10, 15, 20, 25和30 ℃时的COD检测值, 见表2, 以COD浓度为80 mg· L-1的标准水样随温度变化的光谱图为例, 见图4。

图4 80 mg· L-1 COD标准水样随温度变化的光谱图Fig.4 80 mg· L-1 COD spectrogram of standard water sample with different temperature

研究实际水样中温度对COD测量的影响, 采集了COD实际水样在温度为0~30 ℃不同温度下的紫外吸收光谱, 图5为实际水样COD随温度变化的光谱图。

图5 COD实际水样随温度变化的光谱图Fig.5 COD spectrogram of actual water sample with different temperature

2 结果与讨论
2.1 光谱最佳波长的选取

根据朗伯-比尔定律可知, 使用全光谱光谱仪可以测出若干组已知浓度的溶液样品的吸光度。 以此来计算得出COD溶液浓度与吸光度的匹配曲线图, 进而计算得出未知COD溶液的浓度。

采集不同浓度的COD标准水样的紫外吸收光谱, 分别对波长为256, 266, 276, 286和296 nm处的吸光度值进行采集, 依据所得数据建立一元线性回归模型, 并进行线性拟合, 建模结果见表1

表1 COD标准水样在不同波长下统计结果 Table 1 Modeling results of COD standard water sample at different wavelengths
表2 不同浓度COD标准水样随温度变化的吸光度 Table 2 Absorbance of different concentration COD standard water sample with different temperature

一般R2> 0.9时, 模型为优, 0.8< R2≤ 0.9为良, 0.6< R2≤ 0.8为一般, 其余为差, 由表1可知, 286 nm波长处拟合曲线的相关系数r为0.994 6, 决定系数R2为0.989 4, 拟合效果最好, 通过对比和方差(SSE)与均方根误差(RMSE)的情况来作为判断回归方程拟合优劣的另一依据, 其中和方差(SSE)和均方根误差(RMSE)越接近于0, 说明模型选择和拟合越好, 数据预测也越成功, 和方差(SSE)与均方根误差(RMSE)公式见式(1)和式(2), 通过表1可以得出在波长为286和296 nm处的和方差与均方根误差最小, 拟合效果好。 可见COD浓度在286 nm处吸收的最好, 有很强的参考价值, 为最佳检测波长。

SSE=i=1n(yi-y˙i)2(1)

RMSE=1ni=1n(yi-y˙i)2(2)

式(1)和式(2)中yi为模型的预测值, y˙i为标准值, n为采集样本数。

采集浓度为8 mg· L-1标准温度下实际水样与标准水样的紫外吸收光谱, 如图3所示, COD实际水样与标准水样光谱图整体趋势基本一致, 在大部分波段实际水样吸光度大于标准水样, 这是因为实际水样中存在的粘土、 微生物和Cl-, N O2-等离子在紫外区产生吸光度, 实际水样在240 nm处出现最大吸收峰, 次峰出现在245 nm处, 波谷出现在253 nm处, 在275 nm附近再次出现波峰, 接近于标准水样波峰的280 nm, 可见, 在286 nm附近, 在实际水样与标准水样中均出现了吸光度的峰值, 综合上述分析, 故应选用为286 nm作为最佳波长进行实际光谱的探究。

2.2 温度对COD检测精度的影响

由2.1节可得COD最佳检测波长为λ =286 nm, 故对温度研究的实验应在λ =286 nm处进行, 考虑该方法适用于东北地区野外环境, 因此本文重点研究在低温(0~20 ℃)情况下, 温度对测量COD值的影响。 实验以80 mg· L-1COD标准水样为研究对象, 根据图4显示, 在COD浓度一定时, 温度不同, COD标准液的吸光度存在明显差异, 这种现象的原因是温度升高导致COD标准液中水分子动能增大, 原子外层的电子能量随之发生变化, 从而影响COD标准液对紫外光的吸收。

验证光谱法下温度对实际水样COD的影响, 采集COD实际水样在温度为0~30 ℃不同温度下的紫外吸收光谱。 见图5, 从图5可以看出, 温度对其检测值的影响比较大, 随着温度的增加, COD吸光度曲线整体上移, 与图4光谱图上升趋势一致, 故应建立COD温度补偿的模型, 校正温度对紫外光谱法检测COD的影响。 采用实际水样对温度补偿模型进行误差分析。

由图4和图5可知, 在同一浓度中, 温度升高, 吸光度随之增长, 上升趋势一致。 为此采用二元三次多项式表示吸光度值与温度值x、 COD浓度值y的函数关系, 故设

z(x, y)=a0+a1x+a2y+a3x3+a4xy+a5y2+a6x3+a7x2y+a8xy2+a9y3(3)

S( a0* , a1* , …, a9* )是关于未知数 a0* , a1* , …, a9* 的目标函数。 具体见式(4)

S(a0* , a1* , , a9* )=i=1n[(a0* +a1* xi+a2* yi+a3* xi2+a4* xiyi+a5* yi2+a6* xi3+a7* xi2yi+a8* xiyi2+a9* yi3)-zi](4)

其中 a0* , a1* , …, a9* 是式(3)中a0, a1, …, a9的最优解。

分别对 a0* , a1* , …, a9* 导数

Sak* =0 (k=0, 1, , 9)(5)

解得式(5)为

(a0* , a1* , a2* , a3* , a4* , a5* , a6* , a7* , a8* , a9* )T=A-1b(6)

式(6)中

A=i=1n1i=1nyi3i=1nyi3i=1nyi610×10(7)

b=i=1n[zixiziyizixi2zixiyiziyi2zixi3zixi2yizixiyi2ziyi3zi]T(8)

解得式(6)为

a0* =-0.65 64, a1* =0.002 254, a2* =0.012 17, a3* =-6.63× 10-5, a4* =-4.243× 10-5, a5* =1.284× 10-4, a6* =1.158× 10-6, a7* =1.091× 10-6, a8* =4.321× 10-7, a9* =8.566× 10-7

a0* , a1* , …, a9* 代入式(3), 得到所求拟合曲线方程为

$z(x,y)=-0.65 64+0.02 254x+0.012 17y-6.63\times 10^{-5}x^{3}-4.243\times 10^{-5}xy+1.284\times 10^{-4}y^{2}+\\1.158\times 10^{-6}x^{3}+1.091\times 10^{-6}x^{2}y+4.321\times 10^{-7}xy^{2}+8.566\times 10^{-7}y^{3}$(9)

经计算得方程(9)相关系数r为0.995 6, 决定系数R2为0.994 9, 决定系数大于0.9, 由此可见该模型拟合效果良好, 从整体上较好的反映了离散数据的变化规律。 方程(9)经软件仿真后建立的COD温度补偿模型见图6。

图6 COD温度补偿模型图Fig.6 COD temperature compensation model

为了验证COD补偿模型的准确度, 以COD实际水样为研究对象, 依次取水样在0~30 ℃范围内的不同个温度值, 采集该实际水样的COD补偿前值, 将补偿前值经模型补偿后得到表3, 将表3中的所得到的数据进行取点, 具体见图7, 由图7可以看出实际值与补偿后COD值基本吻合, 拟合模型效果较好。

表3 8 mg· L-1实际水样吸光度随温度变化数据表 Table 3 Absorbance of 8 mg· L-1 COD actual water sample with different temperature

图7 COD补偿前、 补偿后与实际值比较图Fig.7 Comparison of COD before and after compensation and actual COD

由于相对误差能反映补偿模型的可靠性, 因此对实际值与补偿后的COD值进行相对误差分析, 相对误差值越小, 模型的补偿效果越好, 表达式如式(10)

Re=y˙i-yiyi×100%(10)

式(10)中, y˙i为补偿后的COD值, yi为实际值, n为采集样本数。 表3列出了每组数据的相对误差, 从表中可以得出, 最大相对误差为6.38%, 最小相对误差为0.63%, 大部分均集中在4%以下, 说明COD温度补偿模型可以比较准确地拟合样本, 个别误差较大的原因可能是光谱测量时水样的微动以及COD测量时烧杯表面的清洁不均所致。

该模型的检测精度为± 4%, 检测结果在误差允许范围之内, 这表示本文所建立的温度补偿模型能够在使用紫外光谱法检测水中COD时, 减弱温度对测量值的影响。

3 结 论

针对利用紫外光谱法低温环境下检测COD精度的问题, 以邻苯二甲酸氢钾粉末配置的COD标准溶液和采集长春市某地区地表水的实际水样为研究对象, 对COD的最佳检测波长和温度对COD检测值的影响进行研究。

实验选取286 nm作为COD检测的最佳波长, 验证了最佳波长对检测野外实际水样中COD的值同样适用。 根据最小二乘原则, 建立吸光度、 温度、 COD浓度值三者之间的COD温度补偿模型。 采用实际水样为研究对象, 同时进行温度补偿模型的误差分析, 结果表明: COD的实际值与补偿后值基本吻合, 最大相对误差为6.38%, 最小相对误差为0.63%, 且多数相对误差集中在4%。 由此可见, COD温度补偿模型补偿值在误差范围内, 补偿精度高, 效果良好。

本研究可对在野外低温环境下使用紫外光谱法检测COD提供有效依据。

参考文献
[1] WANG Li-qun, ZHANG Ming, LIN Feng-mei, et al(王立群, 张明, 林丰妹, ). Chemical Analysis and Meterage(化学分析计量), 2018, 27(3): 113. [本文引用:1]
[2] SONG Jian-jun, ZHAO Ling(宋建军, 赵凌). Transducer and Microsystem Technologies(传感器与微系统), 2018, 37(5): 30. [本文引用:1]
[3] Li J, Tong Y, Guan L, et al. Optik, 2018, 174: 591. [本文引用:1]
[4] LIU Fei, DONG Da-ming, ZHAO Xian-de, et al(刘飞, 董大明, 赵贤德, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2017, 37(9): 2724. [本文引用:1]
[5] YANG Xiao-min, SU Wen-hong, HUANG Xu-bing, et al(杨晓敏, 苏文洪, 黄旭冰, ). City and Town Water Supply(城镇供水), 2019, (1): 34. [本文引用:1]
[6] Hadiyanto, Silviana, PAdetya N, et al. Journal of Physics: Conference Series, 2019, 1217: 012052. [本文引用:1]
[7] YANG Xiao-rong, JIANG Tao, CHENG Ying, et al(杨孝容, 江滔, 成英, ). Journal of Environment and Health(环境与健康杂志), 2018, 35(6): 551. [本文引用:1]
[8] LI Meng-fei, YAO Meng(李萌飞, 姚梦). Journal of Yangtz University·Natural Science Edition(长江大学学报·自科版), 2018, 15(13). [本文引用:1]
[9] YANG Hong, ZHAO Li(杨虹, 赵莉). Tianjin Science & Technology(天津科技), 2015, 42(6): 19. [本文引用:1]
[10] XU Yi-gang, LI Qing, WU Yi, et al(徐熠刚, 李青, 吴轶, ). Computer Measurement & Control(计算机测量与控制), 2017, 25(11): 311 [本文引用:1]
[11] YANG Yun-kai, HE Sheng-hui, YUAN De-fang, et al(杨云开, 何胜辉, 元德仿, ). Intelligent City(智能城市), 2019, 5(16): 140. [本文引用:1]
[12] LI Wen, JIN Xu, ZHANG Zhi-yong, et al(李文, 金旭, 张志永, ). Laser & Optoelectronics Progress(激光与光电子学进展), 2019, 56(13): 131201. [本文引用:1]
[13] Zhou Kunpeng, Bi Weihong, Zheng Qihang, et al. Optoelectronics Letters, 2016, 12(6): 461. [本文引用:1]