基于近红外胶体量子点阵列和光谱重构算法检测多种化学物质的方法研究

王宿慧，张 旭，孙梽珅，杨 杰，郭腾霄，丁学全

国民核生化灾害防护国家重点实验室，北京 102205

引 言

危险化学品是生产建设中的常用物资材料，作为化工大国，我国源于危险化学品的事故时有发生，造成人民生命财产的重大损失。由于化学武器和化学战剂难以消除，且始终是人类共同的重大潜在威胁。危险化学品和化学战剂的检测是处理化学危害和实施防护的首要环节，具有重要和迫切的现实需求和实际意义。当前，在众多化学物质检测技术手段中，红外检测技术由于具有无需预处理、非破坏性、灵敏度高、检测速度快、准确性好等特点，广泛应用于化工、生物医学、航空航天、食品安全等领域。

量子点(quantum dot，QD)是一种半导体纳米材料。由于其具有量子限域效应，即量子点的发射峰会随着量子点尺寸的增大而红移。因此可以通过控制量子点的形状、结构和尺寸来改变量子点的发射和吸收特性。并通过制备分光器件，结合相应的光谱重构算法实现光谱信息获取。具有快速、准确、结构简单紧凑、光通量大、抗震性能好等优点[1-4]。

2015年Bao等[5]实现了微型胶体量子点光谱仪的制作，将胶体量子点材料与聚乙烯缩丁醛混合后作为滤波材料，并使用移液枪将混合物沉积在电荷耦合元件(charge-coupled device，CCD)阵列上，通过测量每个给定的胶体量子点滤波器的总透射强度，计算重构原始(入射)光谱。2019年Yang等[6]通过在探测器上直接生长金属纳米线并作为分光器件，由于光子能量吸收只发生在相应带隙的纳米线上，从而实现了不同位置纳米线上点的特异性吸收，研制出了微米级的微型光谱仪，通过电子探测光电流，对一系列点和线的响应函数进行预校准，就可以计算重建入射光信号。2020年Zhu等[7]通过在光电探测器上沉积钙钛矿量子点，并结合相应光谱重构算法实现了250～1 000 nm波段范围的检测。由于以上微型光谱仪检测波段范围有限，且均为量子点与探测器直接结合，属于被动探测，因此所要求的重构目标光谱必须为类似于太阳光、LED光、白炽灯光等输入光信号，在黑暗条件下化学物质的检测信号较弱。Yan等[8]将3种不同发射波长的PbSe量子点沉积在单个LED光源上，使用匹配气体吸收光谱的方法。虽然可以实现乙炔、甲烷及氮气在矿井中的检测，但是其适用环境较为单一，检测目标较为局限，普适性有待进一步提高。为了解决以上问题，研究中将发射峰不同的近红外胶体硫化铅量子点(PbS)、硫化铅/硫化镉(PbS/CdS)量子点材料与紫外固化胶混合后沉积在RGB点阵屏上，制备成近红外胶体量子点阵列，通过计算机编程实现对阵列的发光和显示控制，结合光谱重构算法实现了900～1 600 nm较宽波段范围内，气态乙醇(C2H6O)样品、液态化学战剂模拟剂甲基膦酸二甲酯(dimethyl methylphosphonate，DMMP)、二氯甲烷(CH2Cl2)样品的分析检测。结果表明，具有64条光谱通道的近红外量子点阵列结合光谱重构算法实现了危险化学品C2H6O气体样品的光谱重构，光谱分辨率为13 nm，重建光谱吸收特征峰与标准光谱信号的偏差为0.407%。具有144条光谱通道近红外量子点阵列结合光谱重构算法实现了DMMP、CH2Cl2液体样品的光谱重构，光谱分辨率为4.861 nm，重建光谱吸收特征峰与标准光谱信号的最小偏差仅为0.043%，完全可以实现多种目标物的检测识别。为复杂环境下危化品泄露事故的应急检测和化学战剂快速检测提供了有效的技术手段和理论支撑。

1 实验部分

1.1 仪器与试剂

主要仪器为近红外光谱仪(美国Ocean Optics公司)。PbS和PbS/CdS量子点材料的标称发射峰为950±50，1 050±50，1 150±50，1 172±20，1 250±20，1 250±50，1 300±20，1 350±50，1 400±20，1 450±20，1 450±50 nm等，均购自星紫新材料和星烁纳米；紫外固化胶为美国Norland公司的NOA61型；检测物包括C2H6O(AR)，DMMP(AR)，CH2Cl2(AR)均购自北京化工厂。

1.2 方法

1.2.1 近红外胶体量子点阵列的制备

将溶剂为甲苯的50 mg·mL-1的PbS和PbS/CdS胶体量子点溶液(具有不同标准发射峰)分别与紫外固化胶以1∶30，1∶50和1∶70的体积比混合，经搅拌和超声处理后制备成均匀的量子点墨水，使用移液枪将25 μL不同种类胶体量子点墨水依次沉积在RGB点阵模块上的各个灯点上，置于紫外灯下照射至完全固化，近红外胶体量子点阵列即制备完毕。

1.2.2 系统构建

构建了用于检测近红外胶体量子点阵列的光致发射光谱(photoluminescence spectroscopy，PL spectrum)测量系统[图1(a)]，该系统主要由电控位移台、支架、光纤、近红外光谱仪以及计算机构成。通过在发射光谱测量系统中加入气室(比色皿)、准直镜、聚焦镜构建了用于检测样品的目标物检测系统[图1(b)]。该系统的工作原理是，将近红外胶体量子点阵列固定在电控位移台上，将光纤固定于支架上，并对准所测近红外胶体量子点阵列的首个灯点，调节准直镜、聚焦镜的焦距确保系统光通量最大，通过对电控位移台编程，控制量子点阵列相对光纤做X和Y轴移动，以S形轨迹遍历所有近红外胶体量子点阵列中的灯点[图1(c)]，光纤另一端连接近红外光谱仪，每选定阵列中的一个灯点进行一次光谱采集，直到把阵列的全部发射光谱采集完毕。

图1 发射光谱测量系统、目标物检测系统及阵列扫描轨迹

2 结果与讨论

2.1 近红外量子点阵列发射光谱特性

当近红外胶体量子点阵列的激发光源和基底为RGB点阵显示模块时，由于该模块的颜色可通过编程控制，因此需要编写不同程序改变点阵屏显示模块的颜色，测量在不同激发光波长的条件下，近红外胶体量子点阵列的发射光谱，通过比较光谱发射强度选择最适合激发近红外胶体量子点阵列的激发光波长。RGB点阵模块用于描述颜色的模式为传统RGB色彩模式，即通过红(R)、绿(G)、蓝(B)三个颜色通道的变化及其相互之间的叠加得到不同颜色，本实验选取了7种不同叠加方案，当颜色数值为255时，发光强度值最大，选取颜色叠加方案如表1所示。

表1 颜色叠加方案

实验随机选择了8个点，每个点经7种不同激发光激发，共58条发射光谱[图2(a)]，不同光谱对比如图2(b)所示，图2中不同颜色的线分别对应相应颜色的激发光。

图2 不同颜色激发光下近红外胶体量子点阵列发射光谱

由图3可知，当RGB点阵显示模块颜色不同时，近红外胶体量子点阵列中相同灯点的发射光谱峰形一致，但强度不同，当RGB点阵显示模块为青色时，量子点经激发后的发光强度最大，白色次之，红色最弱。因此设置点阵模块RGB值为0 255 255，采集近红外胶体量子点阵列的发射光谱[图3(a)]，为了验证近红外胶体量子点阵列的稳定性，在RGB点阵显示模块点亮后的6h内测量了3次近红外胶体量子点阵列的发射光谱，柱状图中共有64个矩形，矩形的高度对应光谱通道3次测量的光强平均值，矩形的中心位置对应光谱通道发射光谱的特征峰，纵向误差棒表示每条光谱通道光强差异，横向误差棒表示每条光谱通道发射光谱特征峰位差异[图3(b)]。

图3 近红外胶体量子点阵列发射光谱及稳定性对比

由图3(b)中的误差棒可知，所测量的3次发射光谱峰位和强度基本吻合，光谱一致性较高，器件稳定性高。且光谱测量全程通过编程控制位移台移动，选择需要测量的量子点灯点，实时性较好。

2.2 目标物检测应用研究

2.2.1 近红外胶体量子点阵列光谱信号处理

为了提高标准光谱测量稳定性和一致性，采用经验模态分解方法(empirical mode decomposition，EMD)提取高频信号以减小随机噪声干扰。即将信号分解为一系列表征不同时间尺度的波动[本征模函数(intrinsic mode function，IMF)]和趋势项(残差)，选择带有特征信息的波动并经过进一步处理，再与趋势项进行重组[图4(a)]。所选原始光谱与EMD预处理后的近红外胶体量子点阵列发射光谱如图4(b, c)所示。

图4 EMD分解及近红外胶体量子点阵列发射光谱

由图4(b,c)可知，与原始光谱相比，经EMD预处理后的近红外胶体量子点阵列发射光谱的随机高频噪声剔除效果良好，标准光谱的随机性大大降低，一致性得到提高。

2.2.2 光谱重构算法研究

将制备完成的近红外胶体量子点阵列放置在扫描平台上，随着扫描系统的推进，检测器将采集近红外胶体量子点阵列每条光谱通道的发射光谱Pi(λ)(其中λ为波长，i=1, 2, …,nF，为光源编号)，当在含有样品室的目标物检测系统的光路中放入目标物后，再次开启扫描模式，采集被样品气体调制后的近红外胶体量子点阵列的发射光谱Ii(λ)，若将待测目标气体的光谱记为φ(λ)，则可以列出方程组，见式(1)

(1)

假设有nλ个离散变量λ，当nλ≤nF时方程组有解，当nλ=nF时，方程组有唯一解，当nλ>nF时方程组无解，但可以拟合出近似解，即可通过解方程得到φ(λ)，实现目标光谱重建。在重构算法流程设计的基础上，基于MATLAB编写了光谱重构程序。目标光谱重构的具体流程如图5所示。

图5 目标光谱重构流程图

2.2.3 目标物检测

光路调节好后，保持其他光学元件固定，将气室取出，将红绿蓝色彩模式(RGB)点阵模块点亮，采集近红外量子点阵列发射光谱。将气室固定，采集近红外量子点阵列透过空气室时的光谱。使用移液枪向气室中注射5微升97%的乙醇溶液，待其挥发后，采集近红外量子点阵列透过乙醇蒸汽样品后的光谱。经验模态分解法(EMD)分解后的发射光谱、透过气室光谱、透过样品光谱图分别如图6(a, b, c)所示。

由图6(c)可以看出来，经过C2H6O蒸汽的近红外胶体量子点阵列的部分发射光谱某一特定波长处出现了强度下降的现象，其为C2H6O的吸收现象。通过最小二乘法和光谱重构算法，拟合得到背景光谱和含样品光谱，再将背景光谱和含样品光谱带入吸收光谱的计算公式中得到C2H6O的吸收光谱。拟合得到的背景光谱和含样品光谱、C2H6O重建吸收光谱和使用光谱仪测得的C2H6O的标准反射光谱分别如图7(a, b, c)所示。

图6 近红外胶体量子点阵列发射光谱、经过气室光谱及经过C2H6O光谱

图7 重建光谱及标准光谱

由图7可知，64条光谱通道近红外量子点阵列结合光谱重构算法实现了微量乙醇气体样品的光谱重构，光谱分辨率为13 nm，重构的乙醇气体吸收光谱在1 390.652 nm处有特征吸收峰，使用近红外光谱仪测得的乙醇气体反射光谱的特征峰为1 396.34 nm。重构目标物的吸收光谱与标准近红外光谱仪测量的吸收光谱吸收特征峰偏差0.407%，基本可以实现目标物的识别。

为了进一步验证近红外胶体量子点阵列与光谱重构算法检测液体样品的可行性，制备了12×12面阵规模的具有144条光谱通道的近红外胶体量子点阵列，将气室更换为微量比色皿，采用同样的方法分别检测了战剂沙林、芥子气的模拟剂DMMP和CH2Cl2。发射光谱、背景光谱、样品经EMD分解后的光谱图分别如图8(a, b, c)(DMMP)，(d, e, f)(CH2Cl2)所示。重建光谱及标准光谱分别如图9(a, b, c)(DMMP)，(d, e, f)(CH2Cl2)所示。

图8 近红外胶体量子点阵列发射光谱、镜过气室光谱及经过样品光谱

由图9(a—f)分别可知，144条光谱通道近红外量子点阵列结合光谱重构算法实现了微量DMMP和CH2Cl2样品的光谱重构，光谱分辨率为4.861 nm，重构的DMMP吸收光谱在1 170和1 385.08 nm处有特征吸收峰，使用近红外光谱仪测得的DMMP吸收光谱的特征峰为1 171.826和1 384.487 nm，偏差分别为0.156%和0.043%。重构的CH2Cl2吸收光谱在1 153.02，1 368.1和1 419.04 nm处有特征吸收峰，使用近红外光谱仪测得的CH2Cl2吸收光谱的特征峰为1 154.234，1 367.153和1 417.513 nm，偏差分别为0.105%，0.069%和0.108%，完全可以实现目标物的识别，因此使用近红外胶体量子点阵列作为分光器件检测多种固态、液态化学物质是可行有效的。后续实验将通过增加有效光谱通道提高重构光谱分辨率，通过提高近红外胶体量子点阵列发光强度提高目标物检测信噪比。

图9 重建光谱及标准光谱

3 结 论

以危险化学品C2H6O、化学战剂沙林、芥子气的模拟剂DMMP、CH2Cl2为目标物，制备了近红外胶体量子点阵列，结合EMD分解减小光谱信号的高频随机噪声，基于最小二乘法建立相应光谱重构算法，实现了上述目标物的检测识别。结果表明，具有64条光谱通道的近红外胶体量子点阵列结合重构算法实现的重构光谱分辨率为13 nm，重构吸收光谱与标准吸收光谱的吸收特征峰偏差为0.407%。具有144条光谱通道的近红外胶体量子点阵列实现的重构光谱分辨率可以达到4.861 nm，与标准吸收光谱相比，其特征峰最小偏差仅为0.043%，完全可以实现900～1 600 nm波段范围内气体、液体目标物的检测识别。后续研究通过增加阵列数量可有效提升重构光谱的光谱分辨率；通过增加所选量子点材料，可以实现从紫外到红外波段范围内的光谱检测；通过优化检测光路和重构算法参数，提高目标物检测信噪比。本研究为复杂环境下危化品泄露事故的应急检测和化学战剂快速检测提供了有效的技术手段和理论支撑。