全文HTML
--> --> -->实验中所采用的IBBCEAS系统如图1所示. 峰值波长约为448 nm的高功率蓝光LED (LZ1-04B2P5, LedEngin)发出的光通过透镜耦合进入光学腔. LED的温度由温度传感器(PT1000)测量, 并通过热电制冷器(TEC)将温度控制在25 ℃ ± 0.1 ℃, 以减少温度波动对光强的影响. 光学腔两侧由两个直径为1英寸的高反镜(0.99982@436.1 nm)组成. 光学腔一端的透射光通过带通滤光片消除杂散光, 经过90°离轴抛物面反射镜后聚焦至光纤(光纤直径Φ为600 μm, 数值孔径NA为0.22). 将光纤的另一端连接到光谱仪(Ocean Optics, QE65000). QE65000光谱仪中的电荷耦合器件(CCD)制冷至–10.0 ℃用于减小暗电流的影响. 光信号通过光谱仪的光栅色散后, 被探测器接收转化为电信号, 经处理后传输到计算机, 进行下一步处理. 本系统主要通过3个质量流量计分别对一路采样气路和两路吹扫气路进行控制. 采样气路和吹扫气路流速分别维持在1和0.2 L/min. 采样管、各种接头和光学腔均采用PFA材质, 采样口前端采用孔径为2 μm的聚四氟乙烯膜过滤器, 用于去除气溶胶以确保整个实验过程中高反镜不被污染.
图 1 基于蓝光LED的非相干宽带腔增强吸收光谱系统组成示意图
Figure1. Schematic of the incoherent broadband cavity-enhanced absorption spectrometer based on blue LED.
3.1.IO拟合波段选择及系统镜片反射率标定
LED发出的光经过增强腔后, 其被光谱仪采集接收后的光谱半高宽主要分布在430—455 nm, IO对应该波段的主要吸收峰分布在436.1, 444.8及454.8 nm处, 如图2所示. 在IO[19]主要的吸收波段同时还存在NO2[20]、水汽[21]和O4等成分的吸收, 采用IBBCEAS系统对IO定量时, 需要考虑高反镜的高反区域与IO的主要吸收峰波段能够较好匹配, 后续实验采用435—465 nm波段作为IO的光谱反演拟合波段.图 2 (a)镜片反射率标定曲线, 黑线为镜片反射率曲线, 红线为氦气谱, 蓝线为氮气谱; (b) 435—465 nm波段主要吸收成分IO, NO2, H2O, O4的吸收截面
Figure2. (a) Reflectivity calibration of the mirror reflectivity. The blue and red curves represent intensity spectrum when the cavity is filled with N2 and He, respectively. The black line represents derived mirror reflectivity curve. (b) Cross sections of NO2 (sky blue line), H2O (grey line), IO (red line) and O4 (black line) in the 435–465 nm band.
由(1)式可知, 镜片反射率是吸收系数公式参数之一. 吸收系数的计算和镜片反射率相关, 因此需要标定镜片反射率来准确获得待测气体的浓度. 利用氮气和氦气之间的瑞利散射间的差异性标定高反镜的镜片反射率[22], 表达式为
2
3.2.有效腔长标定
实际测量时, 采用高纯氮吹扫保护高反镜以避免后者与空气接触时被污染, 氮气吹扫导致系统有效腔长小于高反镜间的距离并影响测量结果, 因此需要对有效腔长进行标定. 我们采用Liang等[23]描述的方法对有效腔长进行标定, 通过在光学腔中通入高纯氧气, 测量打开吹扫和关闭吹扫时的光谱, 将O4吸收截面[24]与光谱仪的仪器函数卷积后, 利用约477 nm处的O4吸收进行光谱反演获得O4浓度. 通入高纯氧流速为1 L/min, 氮气吹扫流速为0.2 L/min, 后续实验过程均保持此吹扫流速. 由于通入的气体是高纯氧气, 因此在开启及关闭氮气吹扫时测量的O4浓度应是唯一的, O4浓度差异是吹扫导致的. 根据(3)式标定系统的有效腔长:图 3 依次关闭和开启N2吹扫时测量的O2浓度时间序列图
Figure3. Temporal variation of O2 concentration with and without N2 purge gas.
2
3.3.系统稳定性及探测灵敏度分析
系统的稳定性会影响其探测灵敏度, 理想的稳定系统可以通过延长积分时间提高探测灵敏度, 然而实际测量过程会受到系统漂移和其他噪声的影响, 导致系统仅在一定时间范围内保持稳定.通过Allan方差评估系统的稳定性并预测最佳采集时间, 利用标准偏差确定系统的探测灵敏度[25]. Allan方差可以通过下式计算:
图 4 IBBCEAS系统的评估 (a)和(b)分别是测量N2背景下光谱反演的NO2和IO的浓度时间序列; (c)和(d)分别是NO2和IO的艾伦方差及标准偏差随平均时间的变化曲线
Figure4. Evaluation of the performance of IBBCEAS instrument. Panels (a) and (b) are the time series of NO2 and IO concentrations with 3 s acquisition time when the cavity is filled with N2. Panels (c) and (d) are the variation curves of Allan and standard deviation plots for NO2 and IO with mean time, respectively
系统 | 时间分 辨率 | 检测限(2σ) | 参考文献 |
LP-DOAS | 60 s | 1.25 pptv | Commane等[2] (2011) |
MAX-DOAS | 60 s | 1.3 × 1013 molecule·cm–2 | Coburn等[6] (2011) |
LIF | 300 s | 0.6 pptv | Gravestocket等[8] (2010) |
CRDS | 30 s | 10 pptv | Wada等[9] (2007) |
ML-CEAS | 300 s | 20 ppqv | Grilli等[10] (2012) |
IBBCEAS | 60 s | 30 pptv | Vaughan等[14] (2008) |
IBBCEAS | 60 s | 4.4 pptv | Ashu-Ayem等[12] (2012) |
IBBCEAS | 22 min | 0.6 pptv | Thalman等[13] (2010) |
IBBCEAS | 60 s | 1.9 pptv | This work |
表1相关IO测量仪器检测限和时间分辨率对比
Table1.Comparison of detection limit and time resolution of correlated IO measuring instruments.
2
3.4.NO2验证实验
IO主要分布在海洋边界层, 而城市边界层检测到IO的浓度较低, 低于系统检测限; 而NO2和IO在435—465 nm波段都存在吸收峰, 因此采用NO2作为验证气体, 通过测量实际大气中的NO2来验证系统的性能. 利用本系统与中心波长为368 nm的IBBCEAS系统[26]及LP-DOAS系统开展NO2的对比实验, 三套系统放置于合肥市科学岛综合实验楼6层, LP-DOAS系统采用开放光路的测量方式, 其角反射镜位于350 m外行政楼楼顶6层, 两套IBBCEAS系统采样管设置在同一位置, LP-DOAS系统放置在IBBCEAS系统的隔壁实验室. 三套系统持续测量3天的NO2浓度时间序列图见图5. 图6(a)为两套IBBCEAS系统的测量数据相关性分析图, 相关系数R2为0.993, 斜率为0.98, 截距为0.02. 图6(b)为中心波长436 nm的IBBCEAS系统(本文系统)与LP-DOAS系统的测量数据相关性分析图, 相关系数R2为0.959, 斜率为1.01, 截距为–0.27. 两套IBBCEAS系统测量结果离散程度优于中心波长为436 nm的IBBCEAS系统(本文系统)与LP-DOAS系统的测量结果离散程度, 其原因可能是IBBCEAS系统是闭路点式采样测量大气中的NO2浓度, 而LP-DOAS系统是开放式光路测量大气NO2的平均浓度, 以及两种方法测量的空气团可能不一致.图 5 三套系统测量的NO2浓度时间序列图, 黑线为中心波长为436 nm的IBBCEAS系统的NO2浓度时间序列, 红线为中心波长为368 nm的IBBCEAS系统的NO2浓度时间序列, 蓝线为LP-DOAS系统的NO2浓度时间序列
Figure5. Comparison of NO2 concentration measured by three different instruments. The black, red and blue dotted lines denote the NO2 concentrations measured by IBBCEAS (center wavelength: 368 nm and 436 nm) and LP-DOAS, respectively.
图 6 (a) IBBCEAS (中心波长368 nm)和IBBCEAS (中心波长436 nm)测量的NO2浓度的相关性; (b) IBBCEAS (中心波长436 nm)和LP-DOAS测量的NO2浓度的相关性
Figure6. (a) Correlation analysis of NO2 concentrations measured by two IBBCEAS instruments (center wavelength: 368 nm and 436 nm); (b) correlation of NO2 concentrations measured by IBBCEAS instrument (center wavelength: 436 nm) and LP-DOAS.
2
3.5.IO损耗标定及线性实验
通过鼓泡法将溶于碘化钾(KI)溶液的碘带出并将其光解后与臭氧反应[16]产生IO样气. 分别将0.031 g的I2和0.839 g的KI溶解于去离子水中, 采用200 mL容量瓶定容配制碘与碘化钾混合溶液(I2/KI/H2O). 将溶于KI溶液的碘光解, 利用氮气作为载气将光解后的碘带出, 将氧气通入臭氧发生装置制备臭氧, 通过臭氧与光解后的碘反应制备IO. 图7为实验室IO的测量示意图.图 7 实验室IO的测量示意图
Figure7. Schematic diagram of laboratory IO measurement system.
由于IO的反应活性较高, 为减小采样损耗, 实验中通过增加采样流速以减少待测气体在空腔中的停留时间. 采样管、各种接头和光学腔均采用PFA材质, PFA的化学惰性可以减少壁效应对测量结果的影响. 在交替使用3 m和10 m采样管下, 以1.2 L/min的流速将实验室制备的IO通入IBBCEAS系统进行检测, 以研究采样管对IO损耗的影响. 获得浓度时间序列如图8所示, 持续测量近2 h的实验结果表明, 3 m采样管内的IO平均浓度约为604 pptv, 其浓度变化约为1.67%, 浓度基本保持稳定, 通过采用3 m采样管时测量的IO浓度为基准评估10 m采样管内的IO损耗, 约为1.98%. 通过改变载气N2的流速, 进而改变IO浓度进行线性实验. 设置臭氧流速恒定为1.2 L/min, 鼓泡氮气流速分别为0.02, 0.04, 0.08, 0.16, 0.32 L/min, 通过测量到的第一组IO浓度为基准, 改变氮气流速进而改变IO配比浓度. 图9(a)为IO线性实验的浓度时间序列图. 整个线性实验的测定时间持续近75 min. 配比的IO浓度范围(39—530 pptv)基本覆盖文献中报道的海洋边界层IO浓度范围. 图9(b)为测量IO的浓度与配比的浓度混合比的相关性, 相关系数R2为0.99.
图 8 IO采样损耗标定, 黑点表示采用10 m PFA采样管测得的IO, 红点表示采用3 m PFA采样管测得的IO
Figure8. Measurements of IO loss in the sampling tube, black dots correspond to the IO measured with the 10 m PFA sample tube, red dots correspond to the IO measured with the 3 m PFA sample tube.
图 9 (a)线性实验下IO的浓度时间序列图; (b)测量IO的浓度与配比浓度的相关性
Figure9. (a) Different concentrations of IO measured by IBBCEAS; (b) the correlation analysis between the average of these concentration gradients and the normalized mixing ratio.
2
3.6.海带排放实验
外场观测已证实海带会排放碘[27], 实验中利用本系统测量海带排放的碘与臭氧反应生成IO. 通过高速重复脉冲(FRR)光诱导单周转快相荧光动力学测量技术测量海带的活性, FRR是基于叶绿素荧光的植物光合速率的快速检测方法, 其原理是将光合作用能流过程分段, 通过分析不同诱导模式下叶绿素荧光动力学曲线, 分段获得主导光合作用能流效率的光合参数; 在此基础上, 根据生物膜能流过程, 建立基于叶绿素荧光的浮游植物光合速率的定量分析方法, 实现植物光合作用状态和生长潜能的实时快速检测[28]. 根据快相荧光模型, 当激发时间趋于零时的荧光产率即为最小荧光值F0[28]. 通过对诱导曲线的前m个荧光曲线数据和激发启动时间t秒采用线性最小二乘算法获得F0[28]. 最大可变荧光Fv即为Fm与F0的差值[28]. 综上即可计算出来光化学量子效率参数Fv/Fm[28]. 通过测量Fv/Fm表征海带的活性, Fv/Fm值越高代表海带活性越强. 海带与海水置于三口烧瓶中, 利用氙灯使用非汇聚的方式对准三口烧瓶, 臭氧通过鼓泡法通入烧瓶, 鼓出的气体通入到IBBCEAS系统进行测量. 不同浓度下的IO反演实例以及拟合残差见图10. 图11(a)和图11(b)分别为IO和NO2的浓度时间序列图, 阴影部分面积表示测量浓度的2σ拟合误差. 分别在实验前后测量海带的Fv/Fm, 实验前均值为0.59, 实验后均值为0.002. 实验前30 min未开启臭氧发生器及氙灯, 用氧气作载气, 结果显示没有IO生成. 在30 min时开启氙灯, IO浓度依然在基线附近. 在60 min时打开臭氧发生器, 发现有IO生成. 测量到的IO最高值为50 pptv, 达到最高值后IO呈下降趋势. 整个海带排放实验的测定时间持续约3 h, 其中发现在1—3 h中IO浓度呈指数型下降, 其原因可能是随实验进行海带的活性不断降低, 碘排放速率下降.图 10 随时间变化的不同浓度的IO反演实例, 最底层是反应时间为137.5 min的光谱反演的拟合残差
Figure10. Example of retrieved absorption spectra of different IO concentrations varying with time, the fitting residuals of spectral retrieved with reaction time of 137.5 min is showed in the bottom of the figure.
图 11 (a) IO的浓度时间序列图, 红色三角形代表海带的Fv/Fm值; (b) NO2的浓度时间序列图, 阴影部分面积表示测量浓度的2σ拟合误差
Figure11. (a) Time series of IO concentrations measured by the IBBCEAS. The red triangle denote the Fv/Fm of luminaria; (b) time series of NO2 concentrations, and the shadow area represent the 2σ deviation error.