不同温度压力对浓度反演精度的定量分析*

王钰豪 刘建国 徐亮 刘文清宋庆利 金岭 徐寒杨

1) (中国科学技术大学环境科学与光电技术学院, 合肥 230026)

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

3) (黑龙江龙凤山大气成分本底国家野外科学观测研究站, 哈尔滨 150209)

自研傅里叶变换红外光谱仪在龙凤山大气本底站测量CO2, CH4 等温室气体.自研仪器的测量结果与符合世界气象组织标准的本底站仪器的测量结果进行对比, 结果表明: 自研仪器与本底站仪器测量的CO2 浓度值相关系数为0.9576, 均方根误差为18.6015.自研仪器使用标准温度、压力下的校准光谱反演浓度, 但测量气体的温度随着气温变化, 导致自研仪器反演浓度有误差.基于以上分析, 提取高分辨率透射分子吸收数据库参数计算吸收截面并结合仪器线形计算不同温度、压力下的校准光谱, 根据不同温度、压力下的校准光谱来校准反演浓度.校准后, 自研与本底站仪器测量的CO2 浓度值相关系数为0.9637, 均方根误差为6.7800.自研与本底站仪器测量的CO2 浓度值相关系数提高, 绝对误差减小, 说明校准算法提高了测量结果的精确度.

1 引 言

温室气体的观测数据是研究气候变化的重要数据基础, 观测技术的发展使得温室气体数据精度不断提高.1995 年, Keeling 等[1]使用非色散红外法(non-dispersive infrared, NDIR)对温室气体CO2进行长期观测.NDIR 和气相色谱法(gas chromatography, GC)具有精度高、稳定性好等优点, 成为测量温室气体的主要方法[2].随着技术的不断发展, 光谱学仪器自动化程度和测量精度不断提高.光腔衰荡光谱(cavity ring-down spectroscopy, CRDS)法、傅里叶变换红外光谱法(Fourier transform infrared spectroscopy, FTIR)仪器也逐渐应用于大气监测领域, 能够实现长期无人观测[3].CRDS 和FTIR 最初应用于抽取式仪器, 定点测量的CRDS 和FTIR 抽取式仪器能够提供高精度、多种类气体数据, 但定点测量方式下CRDS 和FTIR 仪器无法提供区域大气变化数据[4].CRDS和FTIR 仪器在机载、星载观测方式下能获取区域大气变化数据和进行高层大气研究[5,6].仪器在机载、星载观测过程需要校准, 通常校准的两种方式是标准气体校准和利用计算校准光谱校准[7].Griffith[8]在1996 年将校准光谱应用于FTIR, Esler等[9]2000 年研究了在1 cm–1FTIR 仪上对温室气体观测并利用校准光谱校准, 能达到世界气象组织(world meteorological organization, WMO)的测量标准.Hammer 等[10]和Griffith 等[11]对流速、温度、压力等环境因素进行控制, 评估了温度、压力、流量等测量因素的交叉敏感性对仪器测量精度的影响, 他们的研究有效地提高了实验室仪器反演的精度, 减小了使用校准光谱方法因环境变量带来的误差.但在实验中, 仪器无法在控制环境变量的状态下测量属于常态, 本文基于高分辨率透射分子吸收数据库( high-resolution transmission molecular absorption database, HITRAN)计算吸收截面并结合仪器线形产生不同温度、压力的校准光谱, 反演浓度过程使用不同温度、压力的校准光谱.

2 原理与实验

2.1 定量分析物理基础

如图1 所示, 红外光源发出的光经过拋物镜反射成平行光进入分束片(半透射半反射), 分别经过动镜、定镜, 再次汇聚在分束片上形成干涉光, 干涉光经过样品后到达探测器.光源发出的光(IO(v))与经过样品吸收后的光(I(v))进行对比, 由朗伯-比尔定律得透过率T(v)、吸光度A(v)[12]分别为:

图1 干涉仪Fig.1.Interferometer.

气体为均匀非色散吸光物质是(1c)式成立的前提.A(单位1)为气体定量分析光谱区域的吸光度,k(1/molecule × cm–2)为样品的吸收系数,c(mol/mol)为气体的浓度,l(m)为光程长度.利用校准光谱进行浓度反演的过程为: HITRAN 中光谱参数逐线积分计算得到标准温度、压力下的校准光谱, 实际测量的吸光度谱最小二乘拟合校准光谱得到气体浓度[13,14].

2.2 实验条件

自研FTIR 仪器在龙凤山大气本底站的观测现场, 如图2 所示.其中图2(a)为仪器主机在线测量状态, 图2(b)为仪器外界进气口.

外界空气从低层进气口经过前期气体处理装置(除湿、过滤)后, 以1.5 L/min 流速进入怀特池.自研FTIR 仪器的分辨率为1 cm–1[15].光谱采集波段为400—4000 cm–1, 64 次累加测量光谱平均得测量光谱[16].

2.3 实验算法

HITRAN 参数可以预测和模拟大气中光的传输和发射[17].假设能级i和能级j之间跃迁, 真空中光谱线跃迁的波数vij(cm–1)为

HITRAN 数据库提供标准温度下, 单位体积的单个分子的光谱线强度Sij(cm–1/(molecule × cm–2)).温度不为标准温度时, 光谱线强度需要修正, 温度T下Sij(T)为

图2 (a)仪器主机在线测量状态; (b)仪器外界进气口Fig.2.(a) Online measurement status of instrument host; (b) external air inlet of the instrument.

其中,Tref= 296 K,E′′为跃迁能级的低能态能量,配分函数来自于TIPS-2011 数据库,Q(T)为配分函数之和:

考虑分子的不规则运动和分子之间碰撞造成的影响, 分子的不规则运动造成的多普勒展宽的高斯线型为[18]

其中αD(T)高斯线型的半高宽,M为物质的量,NA为阿伏伽德罗常数,kb为玻尔兹曼常数.洛伦兹线型不仅与温度有关, 也与压力有关:

其中γ(p,T)为洛伦兹线型的半高宽.Tref,Pref为296 K 和1 atm.p和pself分别为总气体压力和单独气体的分压,nself为空气自展宽系数,δair为线宽偏移.Voigt 线型为高斯线型和洛伦兹线型的卷积:

光谱线强度由Voigt 线型展宽并在对应波数叠加得吸收系数:

无量纲的光学厚度为吸收系数乘以吸收分子的列数密度(浓度):

仪器线型(instrument line shape, ILS)矩型函数,傅里叶变换后为sinc 函数:

其中γ为仪器的分辨率.光学厚度卷积ILS 得到仿真吸光度谱为[19,20]

以上介绍了吸光度谱的算法过程, 以CO2在2290—2380 cm–1为例, 光谱线强度(红色)、吸收系数(蓝色)、吸光度(橘黄色)如图3 所示.

图3 (a) CO2 在2290—2380 cm–1 的光谱线强度Sij;(b) CO2 在2290—2380 cm–1 的吸收系数kij; (c) CO2 仿真吸光度谱AijFig.3.(a) the spectral line intensity Sij at 2290–2380 cm–1 of CO2; (b) absorption coefficient kij at 2290–2380 cm–1 of CO2; (c) simulated absorbance spectrum Aij of CO2.

2.4 温度、压力校正算法

计算CO2不同温度(297—306 K)下的吸收光谱与296 K 吸收光谱在拟合区域的残差.残差谱进行二次拟合, 进行温度、压力的校正.图4(a)为296 —306 K 下CO2在2290—2380 cm–1的吸收光谱, 297—305 K 吸收光谱与296 K 吸收光谱在2290—2380 cm–1的残差在图4(b)中, 残差与吸收光谱之间相差两个数量级, 有非常丰富的信息.

3 实验系统及测量结果

3.1 仪器状态验证

2019 年07 月23 日, 选取测量的光谱数据进行分析, 计算光谱一致性.

图4 (a) CO2 在2290—2380 cm–1 的仿真吸光度谱Aij; (b) 297—306 K 下仿真吸光度谱相比于296 K 仿真吸收光谱在拟合波段的残差Fig.4.(a) The simulated absorbance spectrum Aij under 286–306 K of CO2; (b) the residual of the simulated absorbance spectrum at 297–306 K compared to the simulated absorption spectrum at 296 K in the fitted band.

图5 (a)光谱数据(1000—3500 cm–1); (b)光谱数据的一致性Fig.5.(a) Spectral data; (b) consistency of spectral data.

STD 为标准偏差(standard deviation),Sn–1,Sn为相邻时间测量得到的两条光谱在2100—2200 cm–1的光谱数据, 计算出的值越接近于1, 说明仪器状态越稳定.2019 年07 月23 日光谱数据如图5(a)所示, 对该光谱数据的一致性进行分析, 结果如图5(b)所示.一致性结果表明仪器状态正常, 测量光谱有效.

3.2 观测结果的分析

用均方根 (root mean square, RMS)误差来衡量CO2绝对浓度的差异.xi表示FTIR 测量的CO2浓度,yi表示相同测量时间CRDS 仪器测量的CO2浓度.公式如下:

均方根误差越小说明两者测量时间内浓度越接近.分析2019 年07 月23 日到2019 年08 月01 日测量的CO2浓度数据, 在图6(a)中使用同一坐标轴显示了CO2浓度.计算了FTIR 与CRDS 仪器测量的CO2浓度的相关系数, 结果如图6(b)所示.FTIR 与CRDS 仪器测量的CO2浓度值的相关系数为0.9576, RMS 误差为18.6015.经过温度、压力校准算法后, FTIR 与CRDS 仪器测量的CO2浓度值相关系数为0.9637, 均方根误差为6.7800.校正后CO2浓度对比结果如图6(c)所示, 线性拟合FTIR 与CRDS 仪器测量CO2浓度值的结果如图6(d)所示.

4 讨 论

图6 (a) CO2 浓度趋势对比; (b) CO2 浓度值线性拟合; (c)校准后CO2 浓度趋势对比; (d)校准后CO2 浓度值线性拟合Fig.6.(a) CO2 concentration trend comparison; (b) CO2 concentration value linear fitting; (c) CO2 concentration trend comparison after calibration; (d) linear fitting of CO2 concentration value after calibration.

观测时, 气体温度、压力并不能保持标准状态.基于标准温度、压力的校准光谱拟合吸光度光谱反演浓度会导致反演浓度误差.图4(b)中给出了不同温度的CO2吸收光谱相对于296 K 下CO2吸收光谱的误差.实际吸收光谱比标准吸收光谱高,则反演得出的CO2浓度会低于实际的CO2浓度.CO2297 K 的吸收光谱与296 K 的吸收光谱形成的残差谱相比于296 K 的吸收光谱的光谱误差总和为1.8%, 反演出的浓度误差为0.41%, 所以计算出不同温度、压力吸收光谱用来校准反演浓度.

5 结 论

自研FTIR 仪器与CRDS 仪器测量CO2的浓度值的相关系数能达到0.9637, RMS 误差为6.7800.以上结果说明了自研FTIR 仪器的可靠性, 本次实验为自研FTIR 仪器的研究工作提供了重要的数据基础.

校准前, FTIR 仪器与CRDS 仪器测量的CO2浓度值的相关系数为0.9576, RMS 误差为18.6015.温度、压力校准后, FTIR 仪器与CRDS 仪器测量的CO2浓度值的相关系数为0.9861, RMS 误差为6.7800.以上结果说明温度、压力校正算法一定程度上可以提高使用校准光谱的仪器测量气体反演的精确度.

计算不同温度、压力下的校准光谱得出的校正系数只能一定程度上校正反演浓度.反演浓度的波段不同会造成校准光谱不同, 那么反演效果可能并不相同.反演波段的识别和选择、反演算法的选择、硬件控制温度、压力, 在使用温度, 压力校准,可以使仪器测量精度进一步提高.