基于迭代算法的大气HONO 和NO2 开放光路宽带腔增强吸收光谱测量*

孟凡昊 秦敏 方武 段俊 唐科 张鹤露 邵豆 廖知堂 谢品华

1) (中国科学院合肥物质科学研究院,安徽光学精密机械研究所,环境光学与技术重点实验室,合肥 230031)

2) (中国科学技术大学,合肥 230026)

3) (皖西学院电气与光电工程学院,六安 237012)

气态亚硝酸(HONO)作为羟基(OH)自由基的重要前体物,在大气中浓度低、寿命短、易损耗且活性强,针对大气HONO 的高灵敏度测量具有一定的挑战.本文介绍了基于迭代算法的开放光路宽带腔增强吸收光谱(OP-BBCEAS)技术应用于大气HONO 和NO2 的测量.常规BBCEAS 技术通过将经滤膜过滤后的环境空气由泵压入/抽入光学腔内进行测量,尽管可以减小气溶胶消光对测量的影响,但针对一些活性组分的测量则需要考虑光学腔和采样造成的吸附损耗和二次生成等壁效应.本文采用OP-BBCEAS 技术,开放光路的测量模式避免了上述壁效应的影响,基于迭代反演算法通过多次迭代确定有效吸收光程,然后采用差分光学吸收光谱的光谱拟合方法对光谱中HONO 和NO2 的吸收进行定量,克服了气溶胶颗粒Mie 散射消光和光源波动的宽带变化影响.在轻度(PM2.5 <75 µg/m3)和中度(PM2.5 >75 µg/m3)不同气溶胶污染状况下测量了实际大气HONO 和NO2 浓度,并同时与常规封闭腔BBCEAS 系统开展了测量对比.不同PM2.5 污染程度下两台BBCEAS系统测量的HONO 和NO2 浓度均显著性相关(R2>0.99),HONO 和NO2 浓度的测量差异(HONO ≤ 4.0%,NO2 ≤6.5%)均小于系统测量误差(HONO:8.1%,NO2:7.5%),验证了迭代反演算法应用于OP-BBCEAS 系统实际大气测量的可行性.

1 引言

气态亚硝酸(HONO)和氮氧化物(NOx=NO +NO2)作为羟基(OH)自由基和臭氧(O3)的重要前体物,是大气活性氮(Nr)循环的重要组成部分,在对流层光化学过程中扮演着至关重要的角色[1,2].HONO 作为一种反应活性气体,在大气中的浓度较低且寿命较短,其具有易损耗、高溶解性的特点.因此,针对大气HONO 的快速、准确、定量测量具有一定的挑战.自1979 年Perner 和Platt[3]采用差分吸收光谱(DOAS)技术首次观测到HONO 以来,包括湿化学法和光谱法的多种测量技术被应用于大气HONO 的测量.湿化学法如长程吸光光度法(LOPAP)利用吸收液吸收气态HONO,通过测量吸收液中亚硝酸根离子或衍生物的含量,实现对于大气HONO 的定量测量,其具有较高的探测灵敏度(pptv 量级,1 pptv=10–12).然而,该方法易受到大气中NOy物种的化学干扰[4],且时间分辨率受限,需要经常标定和更换化学溶液.光谱法作为一种直接测量方法,通过测量HONO 在特定波段(紫外或红外)的特征吸收光谱定量HONO 浓度,不易受到化学干扰的影响,主要包括DOAS 技术[3]、可调谐激光吸收光谱(TDLAS)技术[5]、腔衰荡光谱(CRDS)技术[6]和宽带腔增强吸收光谱(BBCEAS)技术[7−10].其中BBCEAS 是基于CRDS 发展的高灵敏度探测技术,通过增加气体在光学腔内的有效吸收光程,实现对痕量气体的高灵敏度探测.该技术具有结构简单、精度高、适用于外场测量和多种气体同时探测等优点,已成功应用于多种痕量气体(NO2,N2O5,HONO,IO,CHOCHO,NO 等)的高灵敏度探测研究[11−16].

目前,BBCEAS 技术中的光学腔主要采用密封腔体的结构设计,环境大气通过泵采样至光学腔内进行测量,对于大气自由基(如NO3)或易损耗(如HONO)活性气体的测量,需要标定采样损耗和腔壁吸附的影响[11,15].区别于前者,开放光路宽带腔增强吸收光谱(OP-BBCEAS)技术可以避免采样损耗和腔壁吸附的影响,其已经成功应用于大气模拟舱痕量气体如NO3和气溶胶的研究[17−20].此外,Wu 等[21]报道了基于常规BBCEAS 光谱反演算法的OP-BBCEAS 系统应用于实验室清洁空气HONO 和NO2的测量研究,系统积分时间90 s下对HONO 和NO2探测灵敏度(1σ)分别为0.43×10–9和1×10–9.由于常规BBCEAS 反演算法依赖于光强的绝对稳定,观测中发现高气溶胶浓度会影响光谱拟合结果,相较于无气溶胶影响时的光谱拟合,其光谱拟合残差较大[22].Horbanski 等[23]首次提出了基于迭代算法的BBCEAS 技术,利用DOAS反演消除了宽带变化的影响.模拟显示在有效吸收光程减少约80%的情况下,该算法仍然能够准确反演被测痕量气体的浓度.Tang 等[24]验证了该算法对光强波动的不敏感性,并与常规BBCEAS 反演结果进行了对比,不同算法的反演结果显著性相关(HONO:R2=0.94.NO2:R2=0.99).然而,上述算法仅在采用封闭腔结构的BBCEAS 系统上得到了验证,并未开展过实际大气的应用研究.

本文介绍了基于迭代算法的大气HONO 和NO2开放光路宽带腔增强吸收光谱定量测量技术,开放光路的测量模式避免了光学腔和采样管表面损耗和二次生成的影响.基于迭代反演算法,通过多次迭代确定有效吸收光程,然后利用DOAS 反演HONO 和NO2浓度,消除了气溶胶颗粒Mie 散射和光源波动的宽带变化影响.通过对比不同算法的反演结果并开展基于迭代算法的OP-BBCEAS 系统与常规封闭腔BBCEAS 系统的HONO 和NO2测量对比实验,验证了迭代算法应用于OP-BBC EAS 系统测量的可行性.

2 BBCEAS 系统和原理

2.1 BBCEAS 测量原理

BBCEAS 技术是基于痕量气体对光辐射的特征吸收实现对痕量气体的定性和定量测量,其实质是通过光在有限长光学腔内的多次反射增加有效吸收光程,提高对痕量气体的探测灵敏度.由于腔增强吸收光谱技术的有效吸收光程很长,常规封闭腔BBCEAS 技术通常会采用聚四氟乙烯(PTFE)过滤膜过滤气溶胶等颗粒物,消除气溶胶颗粒Mie散射消光的影响.然而,对于OP-BBCEAS 技术而言,不仅需要考虑气体分子吸收和瑞利(Rayleigh)散射的影响,还需要考虑气溶胶颗粒Mie散射消光的影响.常规BBCEAS 反演算法依赖于光源的绝对稳定,通过测量绝对光强的变化确定吸收系数α(λ),利用最小二乘算法拟合被测气体的吸收截面和测量的吸收系数α(λ) 反演被测气体的浓度:

式 中cHONO,cNO2和cO2-O2分别为HONO,NO2和O2-O2的浓度;σHONO(λ) ,σNO2(λ)和σO2-O2(λ) 分别为HONO[25],NO2[26]和O2-O2[27]高分辨率截面与仪器函数卷积后获得的吸收截面;多项式项a1,a2和a3为背景基线中的宽带变化(气溶胶颗粒Mie散射等).由于该算法利用绝对光强的变化反演被测气体的浓度,对环境变化十分敏感,环境空气中气溶胶消光随时间的变化和光源波动均会影响被测气体的浓度反演.

基于迭代算法的OP-BBCEAS 技术,是BBC EAS 技术和DOAS 反演算法的结合[23,24].通过BBC EAS 技术增加有效吸收光程提高对痕量气体的探测灵敏度,利用DOAS 反演算法不受宽带变化影响的特点,使用窄带差分结构反演被测气体浓度,消除了气溶胶颗粒Mie 散射消光和光源波动宽带变化的影响,其光学厚度(DCE(λ))的定义如下:

式中Itot和Itot0分别为有气体吸收和无气体吸收时的透射光强,为有效吸收光程,ci和σi(λ)分别为第i种气体的浓度和吸收截面,εR(λ) 为Rayleigh 散射和Mie 散射的宽带消光.当已知有效吸收光程时,就可以利用DOAS 反演算法反演被测气体的浓度.然而,由于BBCEAS 技术的有效吸收光程强烈依赖于波长,其有效吸收光程不等于“空腔”时的光程

因此,BBCEAS 技术利用DOAS 反演算法反演被测气体浓度时需要修正有效吸收光程.采用Horbanski 等[23]提出的迭代算法修正有效吸收光程,通过多次迭代计算BBCEAS 的有效吸收光程,然后利用DOAS 进行光谱拟合获取HONO 和NO2的浓度信息.该算法具有DOAS 反演算法的优点,可以通过高通滤波的方式消除气溶胶颗粒Mie 散射消光和光源波动宽带变化的影响,而不会影响窄带吸收结构的拟合,实现了开放光路测量模式下HONO 和NO2浓度的准确测量.Horbanski 等[23]和Tang 等[24]已经详细介绍了迭代算法计算有效吸收光程的过程,有效吸收光程可表示为

由(4)式和(5)式可知,若已知光学厚度,则可以修正有效吸收光程.基于迭代算法假定HONO 和NO2的浓度已知,根据(2)式可计算光学厚度[23,24].通过多次迭代确定真实的有效吸收光程,从而获得最终的HONO 和NO2反演浓度,其迭代的停止条件为两次迭代反演的浓度差异小于光谱拟合的误差.

2.2 实验系统

不同于常规封闭腔BBCEAS 系统,OP-BBC EAS 系统实际测量时无光学腔,仅在仪器标定时将光学腔放置在两片高反镜之间,实际测量时将其移除.因此,OP-BBCEAS 系统避免了因光学腔和采样造成的吸附损耗和二次生成的影响.本文采用高反镜间距为33 cm 的OP-BBCEAS 系统来验证基于迭代算法的OP-BBCEAS 技术应用于大气HONO 和NO2测量的可行性.

OP-BBCEAS 系统示意图如图1 所示,系统主要由LED 光源、透镜、高反镜、离轴抛物面镜和光谱仪等组成.LED(LZ1-00 UV00)光源中心波长为365 nm,半峰全宽13 nm,光功率可达到1680 mW.通过PID 算法控制与LED 芯片相连的半导体制冷片,利用热敏电阻温度探头进行实时温度反馈,实现LED 光源恒温(20 ± 0.1) ℃控制.光源发出的光通过消色差透镜(f=50 mm)准直后,耦合到由两片高反镜(d=25 mm)组成的间距为33 cm的开放腔中.开放腔测量状态下,环境空气经PTFE过滤膜(0.2 µm)后由隔膜泵(KNF)抽取至两端高反镜前维持高反镜前端正压,以降低环境空气中气溶胶颗粒物对高反镜造成的污染.透射光由90°离轴抛物面镜(f=25.4 mm)聚焦耦合进入光纤,经光纤传输至光谱仪(QE65000,Ocean Optics)采集光谱信号,通过电脑分析采集的光谱信号获得被测痕量气体的浓度信息.

图1 OP-BBCEAS 系统示意图Fig.1.Schematic diagram of the OP-BBCEAS instrument setup.

3 实验结果与分析

3.1 有效吸收光程

由于光学腔内气体的吸收光程和高反镜镜片反射率相关,需要标定镜片反射率R(λ) 才能准确反演被测气体的浓度.镜片反射率的标定可以通过在高反镜之间增加光学腔后,测量已知浓度的气体吸收或散射引起的光强变化确定,图2 示意了采用高纯氮气(N2,99.999%)和高纯氦气(He,99.999%)的Rayleigh 散射差异性标定高反镜镜片反射率的标定结果.图中显示了高反镜镜片反射率随波长的变化曲线,黑线和红线分别为腔内充满高纯氮气和高纯氦气时的光谱信号,计算的HONO 吸收峰368 nm 处高反镜镜片反射率为0.99983,空腔时的有效吸收光程为1.94 km.

图2 氮气(N2)谱(黑色)、氦气(He)谱(红色)和镜面反射率曲线(蓝色)Fig.2.Nitrogen (N2) spectrum (black line),helium (He)spectrum (red line) and the derived curve of mirror reflectivity (blue line).

不同于常规封闭腔BBCEAS 系统,开放光路的测量模式会受到环境空气中气溶胶颗粒Mie 散射消光的影响,如图3(a)所示.图中显示了环境空气中气溶胶颗粒Mie 散射消光对透射光强的影响,在有和无气溶胶过滤装置时,透射光谱强度变化明显,其峰值光强下降了约21%,估算368 nm 处的气溶胶颗粒Mie 散射消光系数约为 3×10–6cm–1,与污染城市环境中的气溶胶消光水平相当[28].考虑到气溶胶颗粒Mie 散射消光会减少测量时的有效吸收光程,根据有效吸收光程(Leff(λ))表达式[29]:

式中αRay(λ) 和αMie(λ)分别为Rayleigh 散射和Mie散射消光系数;d为腔长;为被测气体总的吸收(10–8cm–1量级),其与高反镜镜片反射率引起的消光(10–6cm–1量级)相比可忽略不计.计算了空腔、腔内充满N2和实际大气测量时的有效吸收光程,如图3(b)所示.相较于空腔时的有效吸收光程,气溶胶颗粒Mie 散射消光显著减少了测量时的有效吸收光程,368 nm 处的最大有效吸收光程仅为1.33 km.

图3 环境空气中气溶胶颗粒Mie 散射对透射光谱强度和有效吸收光程的影响 (a)气溶胶颗粒Mie 散射对透射光谱强度的影响;(b)气溶胶颗粒Mie 散射对有效吸收光程的影响Fig.3.Influence of Mie scattering of aerosol particles in ambient air on transmission spectral intensity and effective absorption optical path:(a) Influence of Mie scattering of aerosol particles on transmission spectral intensity;(b) influence of Mie scattering of aerosol particles on effective absorption optical path.

3.2 HONO 和NO2 的浓度反演

由于BBCEAS 技术采用宽带光源(LED)覆盖了HONO,NO2和O2-O2的特征吸收波段,综合考虑LED 光谱范围、高反镜的高反区域、HONO和NO2吸收峰等因素,选择拟合波段为362.4—389 nm.该拟合波段除HONO,NO2和O2-O2吸收外,其他大气分子(CH2O,O3,BrO,IO 和CHOC HO 等)也有结构性吸收.然而,考虑到其吸收较弱,实际大气浓度远低于系统的探测限,在光谱拟合中可以忽略.大气HONO 和NO2的迭代反演实例如图4(a)所示,拟合得到的HONO 和NO2浓度分别为(1.81±0.11)×10–9和(9.77±0.13)×10–9.光谱拟合残差的标准偏差为1.75×10–4,拟合残差无明显结构,光谱拟合效果较好.

由于常规BBCEAS 反演算法是利用绝对光强的变化反演HONO 和NO2的浓度,气溶胶消光会影响光谱拟合结果,导致光谱拟合结果较差且拟合残差较大,如图4(b)所示.相较于无气溶胶消光影响时的光谱拟合(拟合残差的标准偏差为6.30×10–10cm–1),存在气溶胶消光影响时光谱拟合残差的标准偏差为7.10×10–9cm–1,增加了约1 个数量级.相反,基于迭代反演的OP-BBCEAS 算法,其利用了DOAS 反演算法的优势消除了气溶胶颗粒Mie 散射消光和光源波动宽带变化的影响,在有和无气溶胶消光影响时光谱拟合无明显差异,积分时间为60 s 时HONO 和NO2的探测灵敏度(1σ)分别为95×10–12和270×10–12.

图4 基于迭代反演的OP-BBCEAS 算法和常规BBCEAS 反演算法存在气溶胶消光影响时的大气HONO 和NO2 反演实例 (a)基于迭代反演的OP-BBCEAS 算法存在气溶胶消光影响时的光谱拟合结果,拟合残差的标准偏差为1.75×10–4;(b)常规BBCEAS 反演算法存在气溶胶消光影响时的光谱拟合结果,拟合残差的标准偏差为7.10×10–9 cm–1Fig.4.Examples of HONO and NO2 retrieval of iterative retrieval algorithm and conventional retrieval algorithm with the influence of aerosol extinction:(a) Spectral fitting results of iterative retrieval algorithm with the influence aerosol extinction,the standard deviation of fit residual is 1.75×10–4;(b) spectral fitting results of conventional retrieval algorithm with the influence aerosol extinction,the standard deviation of fit residual is 7.10×10–9 cm–1.

由高斯误差传播可以确定系统的测量误差.OP-BBCEAS 系统的测量误差主要是由高反镜镜片反射率的标定误差、标准吸收截面的测量误差和光谱拟合反演误差组成.其中高反镜镜片反射率的标定误差为5%,文献报道的HONO[25]和NO2[26]标准吸收截面的测量误差分别为5%和4%,光谱拟合反演的误差为4%.根据误差传递函数,HONO和NO2的总测量误差分别为8.1%和7.5%.

3.3 基于迭代算法的OP-BBCEAS 系统与常规封闭腔BBCEAS 系统测量结果对比

为验证基于迭代算法的OP-BBCEAS 系统测量准确性,开展了OP-BBCEAS 系统和常规封闭腔BBCEAS 系统HONO 和NO2测量对比实验,常规封闭腔BBCEAS 系统的详细描述见文献[11].两台系统均放置于合肥市科学岛安徽光学精密机械研究所综合实验楼6 楼(31.89°N,117.17°E),距地面约20 m.常规封闭腔BBCEAS 系统采样管长度约为2 m,采样口放置于OP-BBCEAS 系统附近.考虑到实际测量时OP-BBCEAS 系统采用过滤空气作为吹扫保护气流,为验证过滤空气作为吹扫保护气流是否能够保证高反镜镜片反射率的稳定,标定了实际大气测量前后的高反镜镜片反射率曲线,如图5 所示.环境空气测量前后的高反镜镜片反射率变化差异<0.1%,表明采用过滤空气作为吹扫保护气流能够保证实际测量时高反镜镜片反射率的稳定.

图5 原始镜面反射率曲线(黑色)和实际大气测量后镜面反射率曲线(红色)Fig.5.Initial mirror reflectivity curve (black line) and mirror reflectivity curve after the atmospheric measurements(red line).

基于迭代算法的OP-BBCEAS 系统和常规封闭腔BBCEAS 系统测量不同PM2.5浓度下HONO和NO2浓度的时间序列,如图6 所示.测量期间,PM2.5浓度变化范围为2—149 µg/m3,HONO 浓度变化范围为0.77×10–9—2.80×10–9,平均浓度为(1.75±0.49)×10–9;NO2浓度变化范围为5.20×10–9—25.24×10–9,平均浓度为(13.76 ± 5.12)×10–9.图7(a)和图7(b)分别为轻度(PM2.5<75 µg/m3)和中度(PM2.5>75 µg/m3)不同气溶胶污染状况下两台BBCEAS 系统测量HONO 和NO2浓度的相关性,不同PM2.5污染程度下两台BBCEAS 系统测量的HONO 和NO2浓度均显著性相关.PM2.5浓度小于75 µg/m3时两台BBCEAS 系统测量的HONO 和NO2浓度的相关性系数R2分别为0.998和0.999,HONO 和NO2浓度的测量差异分别为2.4%和6.3%.PM2.5浓度大于75 µg/m3时两台BB CEAS 系统测量的HONO 和NO2浓度的R2分别为0.999 和0.999,HONO 和NO2浓度的测量差异均为4.0%.两台BBCEAS 系统测量的HONO 和NO2浓度差异均在系统的测量误差(HONO:8.1%.NO2:7.5%)范围内,其可能是由于光学腔/采样损耗或测量空气团的差异造成的.

图7 基于迭代算法的OP-BBCEAS 系统和常规封闭腔BBCEAS 系统在不同PM2.5 浓度下测量HONO 和NO2 浓度的相关性(a)轻度(PM2.5<75 µg/m3)和中度(PM2.5>75 µg/m3)气溶胶污染状况下两台BBCEAS 系统测量HONO 浓度的相关性;(b)轻度(PM2.5<75 µg/m3)和中度(PM2.5>75 µg/m3)气溶胶污染状况下两台BBCEAS 系统测量NO2 浓度的相关性Fig.7.Correlation of HONO and NO2 concentrations measured by OP-BBCEAS system based on iterative algorithm and conventional close-path BBCEAS system at different PM2.5 concentrations:(a) The correlation between HONO concentration measured by two BBCEAS instruments in light (PM2.5<75 µg/m3) and moderate (PM2.5>75 µg/m3) aerosol loading;(b) the correlation between NO2 measured by two BBCEAS instruments in light (PM2.5<75 µg/m3) and moderate (PM2.5>75 µg/m3) aerosol loading.

4 结 论

本文介绍了基于迭代算法的开放光路宽带腔增强吸收光谱技术应用于大气HONO 和NO2的定量测量.开放光路的测量模式避免了因光学腔和采样造成的吸附损耗和二次生成的影响,采用过滤空气作为吹扫保护气流,能够保证开放腔测量状态下高反镜镜片反射率的稳定.开放光路测量时,气溶胶颗粒Mie 散射消光显著减少了测量时的有效光吸收光程,相较于空腔时368 nm 处的最大有效吸收光程1.94 km,最大有效吸收光程减少为1.33 km.常规BBCEAS 反演算法受气溶胶消光影响,存在气溶胶消光时的光谱拟合较差,拟合残差的标准偏差较无气溶胶消光影响时增大约1 个数量级.采用迭代反演算法反演HONO 和NO2浓度,通过多次迭代计算有效吸收光程,利用DOAS反演HONO 和NO2浓度,消除了气溶胶颗粒Mie散射消光和光源波动宽带变化的影响,在有和无气溶胶消光影响时其光谱拟合无明显差异,在积分时间为60 s 时HONO 和NO2的探测灵敏度(1σ)分别为95×10–12和270×10–12.开展轻度(PM2.5<75 µg/m3)和中度(PM2.5>75 µg/m3)不同气溶胶污染状况下基于迭代算法的OP-BBCEAS 系统大气HONO 和NO2的测量,并且与常规封闭腔BB CEAS 系统测量结果进行了对比.对比结果显示,不同PM2.5污染程度下两台BBCEAS 系统测量的HONO 和NO2浓度均显著性相关(R2>0.99),HONO 和NO2浓度的测量差异(HONO ≤ 4.0%,NO2≤ 6.5%)均小于系统测量误差,验证了迭代算法应用于OP-BBCEAS 系统实际大气HONO和NO2测量的可行性.未来将利用基于迭代算法的OP-BBCEAS 系统开展大气自由基及其活性前体物的测量,探究自由基及其活性前体物对大气氧化性的影响.