一种重金属浓酸液体中柠檬铬黄浓度的检测方法
【专利摘要】本发明公开了一种重金属浓酸液体中柠檬铬黄浓度的检测方法，该方法以不同柠檬铬黄浓度的重金属浓酸液体作为测试样本，获取各个测试样本在设定波数范围内以硅片为衬底时的拉曼光谱，并根据所有测试样本的拉曼光谱，分别采用连续投影算法和多元线性回归分析法确定定标波数，根据各个测试样本的柠檬铬黄浓度，以及各个定标波数处的峰强与520cm-1处的峰强的比值构建浓度-强度比的第一线性回归模型作为定标模型，利用定标模型测量得到待测样本中柠檬铬黄的浓度。利用Raman测试进行分析，操作简单，不需要进行繁琐、耗时的样品制备过程，同时避免了其他来源的干扰进而保证了测试的柠檬铬黄的准确性，且以硅作为对比大大提高了测试的准确性。
【专利说明】一种重金属浓酸液体中柠檬铬黄浓度的检测方法
【技术领域】
[0001]本发明涉及柠檬铬黄浓度检测【技术领域】，具体涉及一种重金属浓酸液体中柠檬铬黄浓度的检测方法。
【背景技术】
[0002]柠檬铬黄是铅铬黄的其中一种，是一种用作油性合成树脂涂料、印刷油墨、水彩和油彩的颜料，色纸、橡胶和塑料制品的无机着色颜料，由于其具有完善的颜料应用性能、相对低廉的价格和完整的色泽范围，因此得到了广泛的应用。其主要化学成分为铬酸铅，铬酸铅对人体的危害是巨大的，可引起贫血、肾损害、铅中毒、皮炎、湿疹、铬鼻病和皮肤溃疡等，国际癌症研究中心(IARC)已将“铬和某些铬化合物”列入对人类致癌的化学物质。而每生产I吨铅铬黄颜料大约排放出120-150吨废水，废水中一般含有超过国家排放标准5-10倍以上的铅、铬离子及其化合物的悬浮物。废水的治理主要是通过调节液体的PH值，使铅、铬离子反应生成沉淀，以达到去除的效果。
[0003]目前对液体中柠檬铬黄的检测主要是通过测定液体中铅、铬等重金属来评定的，主要的检测方法主要有:原子吸收光谱法、电感耦合等离子体法、原子荧光光谱法和溶出伏安法等。
[0004]原子吸收光谱法是基于蒸气相中被测元素的基态原子对其原子共振辐射的吸收强度来测定试样中被测元素含量的一种方法。该法的优点是选择性强、灵敏度高、分析范围广，但是在多元素检测时不能同时分析，难熔元素的检测灵敏度差，对于基体复杂的样品分析，尚存某些干扰问题需要解决。
[0005]电感耦合等离子体法主要包括电感耦合等离子体原子发射光谱(ICP-AES)法和电感耦合等离子体质谱(ICP-MS)法。ICP-AES是高频感应电流产生的高温将反应气加热、电离，利用元素发出的特征谱线进行测定，它的灵敏度高，干扰小，线性宽，可同时或顺序测定多种金属元素；电感藕合等离子体质谱(ICP-MS)分析技术是将电感藕合等离子体与质谱联用，利用电感藕合等离子体使样品汽化，将待测金属分离出来，从而进人质谱进行测定，通过离子荷质比进行无机元素的定性分析、半定量分析、定量分析，同时进行多种元素及同位素的测定，具有比原子吸收法更低的检测限，是痕量元素分析领域中最先进的方法，但价格昂贵，易受污染。
[0006]原子荧光光谱法(AFS)的原理是原子蒸气吸收一定波长的光辐射而被激发，受激原子随后通过激发过程发射出一定波长的光辐射，在一定的实验条件下，其辐射强度与原子浓度成正比。原子荧光光谱法具有灵敏度高，选择性强，试样量少和方法简单等特点；但其不足之处是应用范围还不够广泛。
[0007]溶出伏安法又称反向溶出极谱法，这种方法是使被测的物质，在待测离子极谱分析产生极限电流的电位下电解一定的时间，然后改变电极的电位，使富集在该电极上的物质重新溶出，根据溶出过程中所得到的伏安曲线来进行定量分析。该方法的灵敏度很高，故在超纯物质分析中具有实用价值，但是影响溶出电流的因素有很多，如富集时间、搅拌速度和电位扫描速率等。
[0008]以上方法都是通过鉴定溶液中重金属铅和铬的存在，进而推断出液体中残留的柠檬铬黄含量，但是在液体处理过程中，无法排除铅、铬的其他来源。所以，单靠重金属铅和铬的检测无法确定液体中铅、铬一定来源于柠檬铬黄。并且用以上方法检测时需要用到大量的试剂进行前处理，过程繁琐，无法做到快速检测。此外，目前对柠檬铬黄废水处理过程中，每一流程里残余柠檬铬黄的监测还鲜见报道。
【发明内容】
[0009]针对现有技术的不足，本发明提供了一种重金属浓酸液体中柠檬铬黄浓度的检测方法。
[0010]一种重金属浓酸液体中柠檬铬黄浓度的检测方法，包括:
[0011](I)以不同柠檬铬黄浓度的重金属浓酸液体作为测试样本，获取各个测试样本在设定波数范围内以硅片为衬底时的拉曼光谱；
[0012](2)根据所有测试样本的拉曼光谱，分别采用连续投影算法提取若干组特征峰，每组特征峰的数量不同，采用连续投影算法时以520CHT1处的列向量作为初始投影向量；
[0013](3)利用多元线性回归分析法确定各组特征峰的验证均方根误差，选择验证均方根误差最小的一组特征峰作为特征指纹峰，并以特征指纹峰作为定标波数，根据各个测试样本的柠檬铬黄浓度，以及相应拉曼光谱中各个定标波数处的峰强与520cm—1处的峰强的比值构建浓度-强度比的第一线性回归模型作为定标模型；
[0014]所述的线性回归模型为:
[0015]Y= 148.288-0.304 λ ^0.222 λ 2-0.08422 λ 3+0.353 λ 4-0.339 λ 5+0.323 λ 6-0.444λ7+0.161λ8-0.03072 λ 9，
[0016]其中，λ!、λ 2、λ 3、λ 4、λ 5、λ 6、λ 7、λ 8 和 λ 9 分别为 1988cm \ 1564cm \979cm \913cm \732cm \678cm \570cm \395cm 1 和 118cm 1 处的峰强与 520cm 1 处的峰强的比值；
[0017](4)获取待测样本在设定波数范围内以硅片为衬底时的拉曼光谱，计算该拉曼光谱中各个定标波数处的峰强和520cm—1处的峰强的比值，并代入定标模型计算得到待测样本中柠檬铬黄的浓度。
[0018]拉曼光谱是基于拉曼效应建立起来的分子结构表征技术，起源于晶体或分子振动(和点阵振动)与转动，拉曼谱线的位置、强度和线宽可提供分子振动、转动方面的信息，可以据此实现分子中某些化学键和官能团的“指纹鉴别”。拉曼光谱作为分子水平的测试手段，易于实现复杂混合物体系的成分鉴定分析。与其他依靠检测属铅和铬元素的方法比较，利用拉曼检测能够保证测试的铅和铬元素来源于柠檬铬黄，进而保证了测试的柠檬铬黄的准确性，避免了其他来源的干扰。
[0019]本发明的重金属浓酸液体多指的是重金属的H2SO4溶液。
[0020]本发明中的硅片多采用单晶硅片，且与测试样本接触面为抛光面，有利于增强520cm_1处的拉曼振动。
[0021]本发明中在连续投影算法时，将硅的特征峰(520CHT1处的峰)的列向量作为初始投影向量，确保了在处理大数据样本时，结果的唯一性，同时也大大加快了数据的处理速度。另一方面，在建模过程中，挑选了不影响溶液结构性质的衬底硅片，以其特征峰做为参照，将测试样本的各个特征指纹峰的强度与520cm—1处的峰强的比值构建定标模型，可以实现拉曼光谱的半定量检测，大大提高了测试的准确性；
[0022]多元线性回归分析是用来研究一个因变量与一组自变量之间的依存关系的，步骤
(3)中根据线性回归分析的结果，选择均方根误差最小的一组特征峰作为特征指纹峰计算待测样本中柠檬铬黄的浓度，能够提高测量结果准确度。
[0023]所述步骤(I)包括如下步骤:
[0024](1-1)将硅片置入容器底部后，向容器中注入测试样本；
[0025](1-2)将注有测试样本的容器放置在显微拉曼光谱仪的载物台上测试该测试样本的拉曼光谱。
[0026]获取以硅片为衬底的Raman光谱(拉曼光谱)时，可以直接将样本均匀涂抹在硅片上，然后将均匀涂抹的硅片放置在显微拉曼光谱仪的载物台上测试样本的拉曼光谱。但是由于液体具有流动性，且所需要的测试样本量是微量的，液体表面存在着张力，直接涂抹无法确保样品表面的平整，易对实验造成影响。其次，采用涂抹时，很难保每次涂抹的测试样本的量刚好相等，从而存在测试误差。本发明中利用容器盛放测试样本，便于对测试样本进行定量，也能够使表面平整，有利于减小因测试条件引起的测试误差。
[0027]本发明中为保证量相同，每次均将容器注满，然后利用刮板沿容器顶面去除多余的液体。
[0028]通常采用圆柱形容器，相应的，所述硅片为圆形，且硅片直径比容器的内径小I?2mm ο
[0029]进行Raman测试时，为保证能够釆集到娃衬底的Raman振动，尽量使娃片能够覆盖整个容器底部，且在测试时尽量不要扫描靠近容器边缘的点。若技术条件允许时，可直接将Si片焊接于容器内的底部，或采用硅材料的容器。
[0030]本发明中Raman测试的测试条件如下:测试激光波长为532nm，测试激光功率为25mv，曝光时间为ls，曝光次数为I次，采集孔径为20 μ m，物镜为20倍，扫描点数为30。
[0031]作为优选,测试样本的数量为50?150。
[0032]单独通过某一个测试样本的Raman光谱很难准确的确定柠檬铬黄的特征指纹峰，本发明中通过对大样本进行统计学分析，能够准确的找出柠檬铬黄振动相关的特征指纹峰。通常样本数越多，特征指纹峰判定越准确，但是这样会导致计算量大，效率低。因此测试样本的数量需要根据实际情况考虑，不可太高，也不能太低。另外，为便于实现所有测试样本可分为若干组(通常为5?7组)，每一组的朽1檬铬黄浓度相同,不同组之间的朽1檬铬黄浓度进行梯度设置。
[0033]作为优选，设定波数范围为112.3?2717.2cm^0
[0034]重金属浓酸液体中与柠檬铬黄相关的振动峰分布在该波数范围内，因此通过设定该112.3?2717.2CHT1的波数范围扫描。
[0035]所述步骤(2)各组特征峰中特征峰的个数为5?15。
[0036]提取得到的特征峰的组数根据实际情况设定，不同组中特征峰的个数互不相同，兼顾到模型的建模效率以及模型的准确性，将特征峰的个数设定为5?15，每组包含的特征峰个数均不同。
[0037]所述步骤(3)中构建定标模型后还包括更新定标波数和定标模型，更新后的定标模型为:
[0038]Y = 142.115-0.211 λ「0.08067 λ 3+0.343 λ 4_0.317 λ 5+0.334 λ 6_0.332 λ 7。
[0039]更新方法如下:
[0040]对所述的第一线性回归模型进行方差分析筛选定标波数，并将筛选后的定标波数作为最终的定标波数，根据所有测试样本的柠檬铬黄浓度，以及筛选后的各个定标波数处的峰强和520CHT1处的峰强的比值构建浓度-强度比的第二线性回归模型作为定标模型。
[0041]通过对第一线性回归模型进行方差分析对特征指纹峰进行筛选，确定得到最终的定标波数，进一步提高定标模型的准确性，进而减小最终测试的浓度值与实际浓度值之间的偏差。
[0042]本发明中未作特殊说明，Y均表示测试样本中柠檬铬黄的含量，其单位为mg/ml。
[0043]与现有技术相比，本发明具有如下优点:
[0044](I)利用Raman测试进行分析，利用拉曼检测能够保证测试的铅和铬元素来源于柠檬铬黄，进而保证了测试的柠檬铬黄的准确性，避免了其他来源的干扰；
[0045](2)利用硅片作为衬底，一方面，由于本实验的检测对象是液体，用硅片作为衬底，可以方便实验中的聚焦；另一方面，娃片的信号峰单一，主要在52001^1处,对测试样本的信号干扰少；同时，以测试样本的各个特征指纹峰的强度与520cm—1处的峰强的比值构建定标模型，可以实现拉曼光谱的半定量检测，大大提高了测试的准确性；
[0046](3)操作简单，避免了传统柠檬铬黄含量测量的提取、消解等繁琐、耗时的样品制备过程，为快速有效地实时监测柠檬铬黄废水处理过程中残余柠檬铬黄的浓度提供有效手段，具有良好的应用前景。
【具体实施方式】
[0047]下面将结合具体实施例和对比例对本发明进行详细描述。
[0048]分别称取lg、0.8g、0.6g、0.4g、0.2g朽1檬铬黄于50ml烧杯中，用移液管移入1mlH2SO4浓溶液，用玻璃棒充分搅拌混匀，静置过夜后，配制成5个浓度梯度的液体样本。每个浓度水平取15个样本，5个浓度梯度的合计得到75个测试样本。
[0049]实施例1
[0050]一种重金属浓酸液体中柠檬铬黄浓度的检测方法，包括:
[0051](I)以不同柠檬铬黄浓度的重金属浓酸液体作为测试样本，获取各个测试样本在2717.2CHT1?112.3cm^波数范围内以硅片为衬底时的拉曼光谱。
[0052]本实施例中采用96孔平底培养皿(孔径:6.4mm，底面积:0.32cm2，容积为0.36ml)作为容器，一个孔即作为一个容器。将每孔的底部首先放置一块直径为5_，厚度为0.5mm的圆形硅片，用移液管移取0.4ml液体于培养皿孔中，用刮板沿着培养皿边缘将多余的液体去除，然后将培养皿放置在玻璃片上，待用。
[0053]将盛有测试样本的容器放置在显微拉曼光谱仪的载物台上测试各个测试样本的拉曼光谱。每个样本的测试条件相同，均如下:
[0054]测试激光波长为532nm,测试激光功率为25mv,曝光时间为Is,曝光次数为I次，米集孔径为20 μ m,物镜为20倍，扫描点数为30。
[0055](2)根据所有测试样本的拉曼光谱，分别采用连续投影算法提取若干组特征峰，每组特征峰的数量不同，采用连续投影算法时以520CHT1处的列向量作为初始投影向量；
[0056]根据测试样本的个数M(本实施例中M = 75)和波数的个数K组成大小为MXK的光谱矩阵X，光谱矩阵X中的元素Xu为第i个测试样本在第j个波数处的拉曼峰的强度值，j = I个波数对应的波数最大，依次向后减小。
[0057]连续投影算法(SPA算法)是一种前向循环选择方法，它从光谱矩阵X中任意一个波长对应的一列组作为投影向量开始，每次循环，计算该投影向量在未选入的波长对应的向量上的投影，再将投影向量的模最大的波长引入到波长组合，直到循环N次。每一次新选入的波长，都与前一个线性关系最小。
[0058]本实施例中一共需要提取11组特征峰，各组特征峰包含的特征峰的个数分别为
5、6、……15。对于每一组特征峰均采用SPA方法进行提取，记每组特征峰的个数为N，则SPA方法包括如下步骤:
[0059](2-1)初始化:n = I (第一次迭代)，在光谱矩阵X中520CHT1对应的元素组成一个列向量作为投影向量，即初始投影向量，记为Xkfe)(即j = k(o)，且第k(0)波数为520CHT1)；
[0060](2-2)集合 S 定义为:S = {j, I ^ j ^ K, j g{k(0),...,k(n-l)H，即还没有被选择进波长链的列向量，其中，k(n-l)表示第η次迭代选出来的最大投影所在的列向量，根据公式:
[0061 ] Pxj = Xj- (x/xkfo-D) Xkto-D (X1 k(n-1) Xk(n-1))，
[0062]分别计算Xj在S集合中各个波数对应的列向量上的投影，并根据公式:
[0063]k (n) = arg (max | | Pxj | , j e S)
[0064]确定使投影最大的j的取值，并记为k(n)，其中I |PXj| I为投影向量在Xj上的投影的模；
[0065](2-3)如果n〈N，则令n = n+1，并返回步骤(2-1)，并以Xk(n)作为初始投影向量，否则停止，并以每次最大投影对应的波数作为特征峰所在的位置，进而得到除520cm—1外，另外包含N个特征峰的一组特征峰。
[0066](3)分别对每一组N个特征峰建立线性回归模型，使用多元线性回归分析法来判断所建模型的优劣，选出最小的RMSEP的一组特征峰作为特征指纹峰。以选择的特征指纹峰作为定标波数，根据各个测试样本的柠檬铬黄浓度，以及相应拉曼光谱中各个定标波数处的峰强与520CHT1处的峰强的比值构建浓度-强度比的第一线性回归模型作为定标模型；
[0067]本实施例中得到的第一线性回归模型为:
[0068]Y = 148.288-0.304 λ ^0.222 λ 2~0.08422 λ 3+0.353 λ 4-0.339 λ 5+0.323 λ 6-0.444λ7+0.161λ8-0.03072 λ 9
[0069]其中，λ!、λ 2、入 3、入 4、入 5、入 6、入 7、入 s 和入 9 分别为 1988cm \ 1564cm \979cm \913cm \732cm \678cm \570cm \395cm 1 和 118cm 1 处的峰强与 520cm 1 处的峰强的比值；
[0070](4)获取待测样本在设定波数范围内以硅片为衬底时的拉曼光谱，计算该拉曼光谱中各个定标波数处的峰强和520cm—1处的峰强的比值，并代入定标模型计算得到待测样本中柠檬铬黄的浓度。
[0071]利用定标模型对这25个待测样本的柠檬铬黄浓度的预测结果如表I所示，模型的相关系数为0.943405，均方根误差为12.11951。说明该模型能够实现浓H2SO4溶液中柠檬铬黄浓度的有效检测。
[0072]表I
[0073]
真实值；1页泖1值真实值I预泖H直真实值预泖H直
(mg/ml)(mg/ml)(mg/ml) (mg/ml) (mg/ml) (mg/ml)
100109.1896061.070 2019.183
10097.7456059.809 20 20.109
100108.1756065.734 20 23.869
100109.6716059.993 20 23.794
100103.3316070.284 2017.941
8082.5394042.992
8077.5474046.000
8088.7714037.571
8074.7044045.059
8072.4824045.771
[0074]实施例2
[0075]与实施例1相同，所不同的是步骤(3)中构建定标模型后还包括更新定标波数和定标模型，更新方法如下:
[0076]对第一线性回归模型进行方差分析筛选定标波数，并将筛选后的定标波数作为最终的定标波数，并根据所有测试样本的柠檬铬黄浓度，以及拉曼光谱中筛选后的各个定标波数处的峰强和520CHT1处的峰强的比值构建浓度-强度比的第二线性回归模型作为更新后的定标模型。
[0077]对定标模型进行方差分析，结果如表2所示
[0078]表2
[0079]
【权利要求】
1.一种重金属浓酸液体中朽1檬铬黄浓度的检测方法，其特征在于，包括: (1)以不同柠檬铬黄浓度的重金属浓酸液体作为测试样本，获取各个测试样本在设定波数范围内以硅片为衬底时的拉曼光谱； (2)根据所有测试样本的拉曼光谱，分别采用连续投影算法提取若干组特征峰，每组特征峰的数量不同，采用连续投影算法时以520cm—1处的列向量作为初始投影向量； (3)利用多元线性回归分析法确定各组特征峰的验证均方根误差，选择验证均方根误差最小的一组特征峰作为特征指纹峰，并以特征指纹峰作为定标波数，根据各个测试样本的柠檬铬黄浓度，以及相应拉曼光谱中各个定标波数处的峰强与520cm—1处的峰强的比值构建浓度-强度比的第一线性回归模型作为定标模型； 所述的线性回归模型为:
Y= 148.288-0.304 λ ^0.222 λ 2-0.08422 λ 3+0.353 λ 4-0.339 λ 5+0.323 λ 6-0.444 λ 7+0.161 λ8-0.03072 λ 9，
其中，λ ρ λ2、λ 3、λ 4、λ 5、λ 6、λ 7、λ8 和 λ9 分别为 1988cm_1> 1564cm_1>979cm_1>913cm \732cm \678cm \570cm \395cm 1 和 118cm 1 处的峰强与 520cm 1 处的峰强的比值； (4)获取待测样本在设定波数范围内以硅片为衬底时的拉曼光谱，计算该拉曼光谱中各个定标波数处的峰强和520cm—1处的峰强的比值，并代入定标模型计算得到待测样本中柠檬铬黄的浓度。
2.如权利要求1所述的重金属浓酸液体中柠檬铬黄浓度的检测方法，其特征在于，所述步骤(I)包括如下步骤: (1-1)将硅片置入容器底部后，向容器中注入测试样本； (1-2)将注有测试样本的容器放置在显微拉曼光谱仪的载物台上测试该测试样本的拉曼光谱。
3.如权利要求2所述的重金属浓酸液体中柠檬铬黄浓度的检测方法，其特征在于，所述硅片为圆形，且硅片直径比容器的内径小I?2mm。
4.如权利要求1所述的重金属浓酸液体中柠檬铬黄浓度的检测方法，其特征在于，测试样本的数量为50?150。
5.如权利要求1所述的重金属浓酸液体中柠檬铬黄浓度的检测方法，其特征在于，设定波数范围为112.3?2717.2CHT1。
6.如权利要求1所述的重金属浓酸液体中柠檬铬黄浓度的检测方法，其特征在于，所述步骤(2)各组特征峰中特征峰的个数为5?15。
7.如权利要求1?6中任意一项权利要求所述的重金属浓酸液体中朽1檬铬黄浓度的检测方法，其特征在于，所述步骤(3)中构建定标模型后还包括更新定标波数和定标模型，更新后的定标模型为:
Y= 142.115-0.211 λ「0.08067 λ 3+0.343 λ 4_0.317 λ 5+0.334 λ 6_0.332 λ 7。
8.如权利要求7所述的重金属浓酸液体中柠檬铬黄浓度的检测方法，其特征在于，更新方法如下: 对所述的第一线性回归模型进行方差分析筛选定标波数，并将筛选后的定标波数作为最终的定标波数，根据所有测试样本的柠檬铬黄浓度，以及筛选后的各个定标波数处的峰强和520CHT1处的峰强的比值构建浓度-强度比的第二线性回归模型作为定标模型。
【文档编号】G01N21/65GK104132931SQ201410364517
【公开日】2014年11月5日 申请日期:2014年7月28日 优先权日:2014年7月28日
【发明者】李晓丽, 孙婵骏, 何勇 申请人:浙江大学