一种近红外光谱测定污水中COD和BOD₅指标的方法

专利名称：一种近红外光谱测定污水中COD和BOD₅指标的方法
技术领域：
本发明属于污水处理领域，尤其是涉及一种近红外光谱测定污水中COD和BOD5指标的方法。
背景技术：
有机物污染是水污染的重要污染物之一。在水质监测中，常采用化学需氧量(COD)和生化需氧量(BOD5)作为重要的综合指标来评价水体中有机物的污染程度。目前国内外测定COD指标的方法主要有重铬酸钾法、库伦法和比色法；测定BOD5指标的方法主要有稀释与接种法、荧光法和生物传感器法。这些测定方法通常需要对样品进行预处理，操作复杂，有时需要用到有毒试剂，不仅危害身体健康，而且会造成二次污染。近年来，近红外光谱分析技术迅速发展。利用该技术分析样品具有快速高效、不破坏样品、不消耗化学试剂、无污染、便于在线分析等优点，因此在诸多领域得到了广泛应用。近红外光谱在水质监测方面的应用亦有报道:对污水样品的近红外光谱简单预处理后直接建立测定模型；采用优化算法对污水样品的近红外光谱进行优化，然后建立测定模型；通过测定光程的变化与组合，提高近红外光谱测定模型的灵敏度。但由于污水中的水分在近红外光谱区域有强吸收，会对污水中待测定指标的光谱特征造成强烈干扰，导致最终测定指标的误差较大。

发明内容
本发明的目的是针对上述现有技术中存在的不足，提供一种近红外光谱测定污水中COD和BOD5指标的方法。本发明为实现上述目的，所采用的技术方案为:一种近红外光谱测定污水中COD和BOD5指标的方法，包括下述步骤:
(1)取污水样品，先用常规方法测定污水样品的COD和BOD5的标准值，然后测定污水样品的近红外光谱，得到污水样品的原始近红外光谱数据；
(2)取纯水样品，在与污水样品相同的条件下测定纯水样品的近红外光谱，得到纯水样品的近红外光谱数据；
(3)对污水样品的原始近红外光谱数据进行预处理后，用凯纳德-斯通算法将污水样品划分为校正集和验证集，通过计算扣除污水样品中纯水的近红外光谱，得到污水样品的剩余近红外光谱数据；
(4)将校正集污水样品的剩余近红外光谱数据分别与用常规方法测定的COD和BOD5的标准值进行关联，用偏最小二乘法分别建立污水中COD和BOD5指标的测定模型；
(5)将验证集污水样品的剩余近红外光谱数据代入步骤(4)建立的测定模型得到COD和BOD5的测定值，并与用常规方法测定的COD和BOD5的标准值比较，对测定模型进行验证、评价；
(6)取未知污水样品，测定其近红外光谱，重复步骤(3)后，得到未知污水样品的剩余近红外光谱数据，代入测定模型，即可测得未知污水样品中COD和B0D5指标。本发明的一种近红外光谱测定污水中C0D和B0D5指标的方法，步骤(3)中所述的 扣除污水样品中纯水的近红外光谱的计算方法为
①构建与纯水样品的近红外光谱数据正交的转换矩阵Z，构建公式如下
z=i -wx (w’父胃厂丨父^’其中’^为纯水样品的近红外光谱矩阵，I为单位矩阵，r为 w的转置矩阵；
②根据转换矩阵z，计算污水样品的剩余近红外光谱数据，计算公式如下
R=ZXS，其中，R为污水样品的剩余近红外光谱矩阵，S为污水样品的原始近红外光谱 矩阵。本发明的一种近红外光谱测定污水中C0D和B0D5指标的方法，步骤(3)中所述的 预处理为采用平滑滤波器对污水样品的原始近红外光谱数据进行滤波。本发明的一种近红外光谱测定污水中C0D和B0D5指标的方法，所述的测定模型对 应的近红外光谱的光谱范围包括三个区域，分别为I区域ASOOcnT1 SSOOcnT1、II区域 7556cm 1 7162cm ^III区域6329cm 1 5470cm 丄。本发明的一种近红外光谱测定污水中C0D和B0D5指标的方法，采用正交信号校正 法，污水样品的近红外光谱特征信息S可表示为污水样品中有机污染物的光谱特征信息S w、水的光谱特征信息W、污水中高频噪音等杂质光谱特征信息三者的加和，即S=AW+S
通过对污水样品近红外光谱数据的滤波预处理可去除高频噪音等杂质(Ssft)对光 谱的干扰。用污水样品的近红外光谱数据与纯水的正交转换矩阵相乘后，可从污水的近红 外光谱特征信息中扣除纯水的近红外光谱特征信息，以此来消除水在近红外区域的强吸收 对C0D和B0D5指标测定造成的干扰，提高了最终建立的测定模型的预测能力和精确度。本发明的一种近红外光谱测定污水中C0D和B0D5指标的方法，对未知污水样品测 定其近红外光谱数据并进行数据处理后，通过建立的测定模型直接计算得到未知污水样品 中C0D和B0D5指标。测定时间快、无需任何化学试剂，具有简便、迅速、准确、易操作的优点。 用本发明方法监测水质，大大提高了水质预警的时效性，具有较强的实用性。


图1为污水样品的原始近红外光谱。图2为纯水样品的近红外光谱。图3为污水样品的剩余近红外光谱。图4为未知污水样品的C0D指标的测定值和标准值关系图。图5为未知污水样品的B0D5指标的测定值和标准值关系图。
具体实施例方式本发明的一种近红外光谱测定污水中C0D和B0D5指标的方法，包括下述步骤
(1)取污水样品，先用常规方法测定污水样品的C0D和B0D5的标准值，然后测定污水样 品的近红外光谱，得到污水样品的原始近红外光谱数据；
(2)取纯水样品，在与污水样品相同的条件下测定纯水样品的近红外光谱，得到纯水样 品的近红外光谱数据；(3)对污水样品的原始近红外光谱数据进行预处理后，用凯纳德-斯通算法将污水样品划分为校正集和验证集，通过计算扣除污水样品中纯水的近红外光谱，得到污水样品的剩余近红外光谱数据；
(4)将校正集污水样品的剩余近红外光谱数据分别与用常规方法测定的COD和BOD5的标准值进行关联，用偏最小二乘法分别建立污水中COD和BOD5指标的测定模型；
(5)将验证集污水样品的剩余近红外光谱数据代入步骤(4)建立的测定模型得到COD和BOD5的测定值，并与用常规方法测定的COD和BOD5的标准值比较，对测定模型进行验证、评价；
(6)取未知污水样品，测定其近红外光谱，重复步骤(3)后，得到未知污水样品的剩余近红外光谱数据，代入测定模型，即可测得未知污水样品中COD和BOD5指标。步骤(3)中所述的扣除污水样品中纯水的近红外光谱的计算方法为:
①构建与纯水样品的近红外光谱数据正交的转换矩阵Z，构建公式如下:
z=1-ffx (w’ Xwr1Xw'其中，w为纯水样品的近红外光谱矩阵，I为单位矩阵，r为w的转置矩阵；
②根据转换矩阵z，计算污水样品的剩余近红外光谱数据，计算公式如下:
R=ZXS,其中，R为污水样品的剩余近红外光谱矩阵，S为污水样品的原始近红外光谱矩阵。该计算方法能够有效消除水的光谱吸收在COD和BOD5指标测定过程中造成的干扰，提高测定模型的预测能力和精确度。步骤(3)中所述的预处理为采用平滑滤波器对污水样品的原始近红外光谱数据进行滤波。通过对光谱数据的滤波可消除高频噪音等杂质对污水样品近红外光谱数据的干扰，提高测定模型的精确度。所述的测定模型对应的近红外光谱的光谱范围包括三个区域，分别为:I区域:8800cm 1 8500cm \ II 区域:7556cm 1 7162cm 1、TTT区域:6329cm 1 5470cm、在所述的
1、I1、111区域内通过计算能够将水的近红外光谱造成的干扰降至最低，使最终测定模型的精确度达到最高。本发明所选污水样品数量越多，所建立的测定模型越准确、精度越高，但建立测定模型的工作量也越大，因此，实际操作中为减少工作量，一般选取尽可能涵盖所有可能预测值的样品数量，本发明优选的污水样品数量在60-120个。纯水样品不同样品之间无显著差异，优选数量为5-10个。本发明步骤(I)中所述的测定COD指标的常规方法参照GB11914-89 (水质-化学需氧量的测定-重铬酸盐法)进行测定；所述的测定BOD5指标的常规方法参照GB7488-87(稀释与接种法)进行测定。本发明步骤(3)中所述的用凯纳德-斯通算法将污水样品进行划分为校正集和验证集。凯纳德-斯通算法(kennard-stone算法的音译)的基本原理为:根据不同污水样品之间的最大欧氏距离选择最具有代表性分布的污水样品为校正集。根据凯纳德-斯通算法，第i个污水样品到第j个污水样品的欧氏距离的计算公式为:
权利要求
1.一种近红外光谱测定污水中COD和BOD5指标的方法，其特征在于:包括下述步骤: (1)取污水样品，先用常规方法测定污水样品的COD和BOD5的标准值，然后测定污水样品的近红外光谱，得到污水样品的原始近红外光谱数据； (2)取纯水样品，在与污水样品相同的条件下测定纯水样品的近红外光谱，得到纯水样品的近红外光谱数据； (3)对污水样品的原始近红外光谱数据进行预处理后，用凯纳德-斯通算法将污水样品划分为校正集和验证集，通过计算扣除污水样品中纯水的近红外光谱，得到污水样品的剩余近红外光谱数据； (4)将校正集污水样品的剩余近红外光谱数据分别与用常规方法测定的COD和BOD5的标准值进行关联，用偏最小二乘法分别建立污水中COD和BOD5指标的测定模型； (5)将验证集污水样品的剩余近红外光谱数据代入步骤(4)建立的测定模型得到COD和BOD5的测定值，并与用常规方法测定的COD和BOD5的标准值比较，对测定模型进行验证、评价； (6)取未知污水样品，测定其近红外光谱，重复步骤(3)后，得到未知污水样品的剩余近红外光谱数据，代入测定模型，即可测得未知污水样品中COD和BOD5指标。
2.根据权利要求1所述的一种近红外光谱测定污水中COD和BOD5指标的方法，其特征在于:步骤(3)中所述的扣除污水样品中纯水的近红外光谱的计算方法为: ①构建与纯水样品的近红外光谱数据正交的转换矩阵Z，构建公式如下: z=1-ffx (w’ Xwr1Xw'其中，w为纯水样品的近红外光谱矩阵，I为单位矩阵，r为W的转置矩阵； ②根据转换矩阵Z，计算污水样品的剩余近红外光谱数据，计算公式如下: R=ZXS,其中，R为污水样品的剩余近红外光谱矩阵，S为污水样品的原始近红外光谱矩阵。
3.根据权利要求1所述的一种近红外光谱测定污水中COD和BOD5指标的方法，其特征在于:步骤(3)中所述的预处理为采用平滑滤波器对污水样品的原始近红外光谱数据进行滤波。
4.根据权利要求1所述的一种近红外光谱测定污水中COD和BOD5指标的方法，其特征在于:所述的测定模型对应的近红外光谱的光谱范围包括三个区域，分别为:I区域:8800cm 1 8500cm \ II 区域:7556cm 1 7162cm 1JII区域:6329cm 1 5470cm 1O
全文摘要
本发明提供了一种近红外光谱测定污水中COD和BOD5指标的方法，具体步骤如下(1)取污水样品，常规方法测定COD和BOD5的标准值，然后测定污水样品的近红外光谱；(2)取纯水样品，测定近红外光谱；(3)对污水样品的光谱数据预处理后，将污水样品划分为校正集和验证集，通过计算得到污水样品的剩余近红外光谱数据；(4)将校正集污水样品的剩余近红外光谱数据与COD和BOD5的标准值进行关联，建立测定模型；(5)对测定模型进行验证；(6)取未知污水样品，测定近红外光谱，重复步骤(3)后，将剩余光谱数据代入测定模型，得到COD和BOD5指标。本发明的方法具有简便、迅速、准确、易操作的优点。
文档编号G01N21/35GK103175805SQ201310076050
公开日2013年6月26日 申请日期2013年3月11日 优先权日2013年3月11日
发明者楚广诣, 聂磊, 贾瑞宝, 孙韶华, 孔令华, 姜红, 乔常红, 王辉, 李涛, 胡力华 申请人:楚广诣