郑海明，朱小朋，冯帅帅，贾桂红

华北电力大学机械工程系，河北 保定 071003

引 言

苯-甲苯-二甲苯(Benzene,Toluene and Xylene，BTX)是大气挥发性有机物(volatile organic compounds，VOCs)的主要组成成分。在公路两侧空气的VOCs中，BTX含量达到46%。此外，工业溶剂的使用、汽油的挥发、家庭室内装修、化工厂泄漏等也是BTX的主要来源，工人长期暴露在苯系物的环境中致癌风险将极大提高[1-3]。如何准确有效地监测BTX浓度很有必要。

目前，针对BTX的监测方法主要有：GC-MS、离子化检测器、FTIR、激光诱导荧光、NDIR和DOAS技术等[4-6]。GC-MS、离子化检测器、FTIR等价格较贵，需对样品进行预处理，适用于实验室标准分析，而工业现场的高污染以及复杂的环境条件对这几种方法的测量精度影响严重。差分吸收光谱(differential optical absorption spectroscopy,DOAS)可避免这种污染的影响[7-8]，能够实现在线定性定量分析。

本实验研究利用紫外光源对以上苯系物质连续激发，并结合差分吸收光谱技术实现对单组分苯系物以及BTX混合气体的同时在线监测。

1 差分吸收光谱技术监测原理

差分吸收光谱是一种利用气体物质对吸收谱线的差分吸收特征来反演浓度的方法，目前已经发展成为一种监测大气环境污染物的重要手段[9]。DOAS的核心基础是修正后的Lambert-Beer定律，如式(1)所示

(1)

式(1)中：I0(λ)为光源发出的原始光谱强度；I(λ)为光谱仪接收到的测量光谱强度；σi(λ)为气体的标准吸收截面；L为气体吸收光程；εMie(λ)和εRay(λ)分别为米散射系数和瑞利散射系数。

(2)

(3)

式(3)中OD′表示差分吸光度。

2 实验部分

设计了一套由标准BTX液体制备标准气态BTX的装置，如图1所示。实验系统采用7ISU75系列三光栅单色仪，分辨率0.01 nm，光源采用7ILX150UV紫外增强型氙灯光源作为发射端，其连续光谱覆盖整个紫外波段，在250～280 nm范围光源具有良好的稳定性。CS200气体质量流量计其精确度达到0.35%。混合气室为316L不锈钢制成，光程长度为400 mm。标准液体为北京计量院生产的单组分苯、甲苯、二甲苯、甲醇(色谱纯)，其纯度均达到99.99%以上。

图1 实验系统1:氮气源；2:PTFE管；3:CS200流量控制阀；4:BTX标准液；5:硅胶软管；6:BT100-2J精密蠕动泵；7:雾化器； 8:混合气室；9:加热装置；10:过滤器；11:紫外增强型氙灯光源；12:吸收池；13:废气流出端；14:数据采集系统； 15:PC数据处理系统；16:废液回收罐；17:壳体Fig.1 Experimental system1:Nitrogen source;2:PTFE tube;3:CS200 flow control valve;4:BTX standard fluid;5:Silicone hose;6:BT100-2J precision peristaltic pump;7:Atomizer;8:Mixing gas chamber;9:Heating device;10:Filter;11:Ultraviolet enhanced xenon light source;12:Absorption cell;13:Exhaust gas outflow end;14:Data acquisition system;15:PC data processing system;16:Waste liquid recovery tank;17:Shell

3 结果与讨论

为了更好地研究BTX混合气体的浓度监测，首先以单组分苯反演在单一物质条件下浓度测量情况，然后利用上述实验装置对BTX混合气体组分进行实验与分析。

3.1 单组分苯实验结果

针对单组分气态苯特征吸收，选取252.5～253.8 nm作为研究波段，当充入单组分苯标气时，基于差分吸收光谱法的基本原理，其反演浓度的关键在于获取差分吸收截面及差分吸光度。在实际的测量过程中，系统及环境噪声会叠加在吸收光谱上，影响测量精度。传统处理吸收光谱噪声的方法常采用多项式平滑滤波，会使吸收截面上的细节信息或高频分量丢失。小波变换具有良好的时频局域化特性，能通过伸缩和平移对信号进行多分辨率分析并可聚焦到信号的任意细节。

小波去噪的基本原理[10]：给定一个信号，首先将信号展开成小波的平移和伸缩之和，然后把欲舍弃项的系数去掉或加以适当的修改，由修改后的小波系数做逆向小波变换得到去噪后的信号。去噪效果采用信噪比公式评价，如式(4)

(4)

图2 平滑去噪图Fig.2 Schematic of smooth denoising

图3 小波变换去噪图Fig.3 Schematic of Wavelet transform denoising

吸收截面是反映气体对特征光谱吸收能力强弱的常量，它与物质的种类、温度和压力有关，与浓度无关。研究中通过实验获取吸收截面，计算的基础为Lambert-Beer定律，在一定压力和温度下获取入射光谱和衰减光谱[11-13]，由式(5)计算得到

(5)

由于标准吸收截面在实验室中测定，其分辨率较高，权威HITRAN数据库中标准吸收截面给出的只是BTX在少数几个温度和压强下的吸收截面，对于实际条件下的温度和压力不能够完全匹配。为了能够更好地得到反演气体的浓度，提高测量精度，通过测定在实验室环境下的温度、压力来获取吸收截面。如图4为本研究获取的单组分苯吸收截面与HITRAN数据库当中的苯标准吸收截面对比。

由图4可知：根据式(5)获取的吸收截面与理论标准吸收截面图相比较，整体变化趋势一致，最高点处与标准吸收截面有差异。分析认为是由于在实际的测量当中，其温度、压强并不是标准状态且外界环境噪声等存在干扰，这种差异从另一方面也反映出采用实验获取吸收截面的灵活性。若直接采用HITRAN数据库中吸收截面值将造成由于温度与压强变化导致的浓度反演误差。为了能够与实际监测环境相符合，采用实验室获取的吸收截面作为标准吸收截面。

图4 苯吸收截面对比图Fig.4 Comparison of benzene absorption cross-section

研究中分别采用二项式系数滤波、多项式拟合和Savitzky-Golay滤波法去除慢变部分得到差分吸收部分，如图5所示为三种滤波方法的对比图。

图5 滤波效果对比图Fig.5 Comparison of filtering effect

由图5中可以看出多项式拟合与Savitzky-Golay滤波拟合相比二项式拟合更能提取快变特征，通过对甲苯浓度值的反演得到Savitzky-Golay滤波法多项式拟合误差更低，去除慢变部分效果更好，故选取了Savitzky-Golay滤波方法去除慢变部分。将差分吸光度与差分吸收截面最小二乘法拟合得到反演浓度。

通过多组实验发现不同浓度下的苯与横坐标包围面积与浓度之间有一定的比例关系，如图6所示。为探究这种比例关系，以单组分苯为例探究该比例关系。实验中配置了7.314，9.136，13.494，15.884，19.884，23.606和30.764 mg·m-3等7组不同浓度值的标准气态苯，为能够更好地表述积分面积与浓度之间的关系，采用多项式拟合该数据点如图7所示，最终建立关系如式(6)

图6 气态苯吸光度Fig.6 Absorbance of gaseous benzene

图7 苯拟合曲线Fig.7 Benzene fitting curve

y=2×10-3x3-0.0306x2+2.058x-32.82

(6)

式(6)中，x为苯积分面积，y为苯的浓度值。

为验证式(6)的合理性，研究中再次配置了3种不同浓度的气态苯，测量两种物质的吸光度积分面积并分别代入式(6)，另一方面与采用传统的最小二乘法反演浓度作比较。如图8所示为气态苯的差分吸光度，实验结果如表1所示。结果显示：利用积分面积法测量浓度误差均在10%以内，最小二乘法测量出的浓度值其误差均在5%以下，通过对比两种算法反演误差可得知最小二乘法更加稳定，精度更高。分析原因：(1) 单组分气态苯的吸光度出现单一峰值的现象，利用多项式拟合去除慢变化过程当中，拟合效果不理想。(2) 积分面积法能避免因光谱仪器设备老化出现的光谱偏移造成的计算误差。根据HJ/T76—2007标准中准确度的规定可知积分面积法的计算精度在规定范围内，能够满足固定污染源等工业领域对苯系物的实时在线监测要求。

图8 苯差分吸光度曲线Fig.8 Differential absorbance curve of benzene

表1 气态苯积分面积法和最小二乘法反演浓度结果Table 1 Inversion results of gaseous benzene based on integrated area method and least square method

3.2 BTX混合组分实验结果与分析

在吸收池内充入不同浓度的BTX混合气体进行实验和分析，苯在253和259 nm处有特征吸收峰，甲苯在267 nm处有较强吸收峰，二甲苯在250～275 nm处均有吸收峰，故选取了255～275 nm波段作为研究波段，该波段可将BTX的特征吸收均包含其内，如图9(a)；由于苯、甲苯在270～275 nm之间没有明显特征吸收，故先反演出二甲苯的浓度，利用二甲苯浓度值反演出在255～275 nm波段的二甲苯的差分吸光度如图9(b)，在该波段内利用总差分吸光度去除通过反演得到的二甲苯的差分吸光度便可得到在此波段内苯、甲苯混合差分吸光度，在波段265～270 nm波段甲苯有特征吸收而苯没有特征吸收，故可认为此波段的差分吸光度仅为甲苯的差分吸光度如图9(c)，通过最小二乘法反演得到甲苯浓度值进而反演得到255～275 nm波段的甲苯差分吸光度，通过剩余差分吸光度去除获得苯的差分吸光度如图9(d)，进而可得苯的浓度值。

图9 BTX混合浓度反演过程(a)：BTX差分吸光度；(b)：BTX差分吸光度与二甲苯差分吸光度；(c)：苯、甲苯差分吸光度与甲苯差分吸光度；(d)：苯差分吸光度Fig.9 BTX mixed concentration inversion process(a)：BTX differential absorbance；(b)：BTX differential absorbance and xylene differential absorbance； (c)：Differential absorbance of benzene and toluene；(d) Differential absorbance of benzene

利用上述装置在同等条件下配置如表2的BTX混合气体。由于最小二乘法反演浓度更为稳定，误差更小，故混合气体浓度反演采用最小二乘法反演其浓度值，反演得到如表3所示结果。从表3中可以看出，第二组数据的各项误差均处于极小值，BTX浓度的反演结果比其他组好。

表2 BTX混合气体配比(mg·m-3)Table 2 BTX mixed gas ratio (mg·m-3)

表3 不同浓度的混合气体反演结果(mg·m-3)Table 3 Inversion results of different concentrations mixed gases (mg·m-3)

4 结 论

对单组分气态苯进行了浓度反演实验，并在此基础上进一步对BTX混合气体做了研究，最终实验结果表明：积分面积法能够应用于特征吸收为单峰值波形的物质，但最小二乘法相比积分面积法测量误差更小，稳定性更好；利用反演差分吸光度的方法对混合BTX混合物的浓度反演二甲苯测量误差均在2%以下，但对甲苯和苯的测量误差逐渐增大，苯的反演最大误差达到了9.07%，苯的测量精度受到了二甲苯、甲苯测量精度的影响以及苯特征吸收波段特点的影响。针对BTX混合气体反演精度的提高需要探究更优良的算法或更加精密的光谱测量设备。