基于波长调制-直接吸收光谱方法的CO分子1567 nm处谱线参数高精度标定*

王振 杜艳君 丁艳军 彭志敏

(清华大学能源与动力工程系, 电力系统与发电设备控制与仿真国家重点实验室, 北京100084)

(2019 年 12 月 9日收到; 2019 年 12 月 19 日收到修改稿)

直接吸收光谱(DAS)可直接测量分子吸收率函数, 并通过拟合吸收率函数确定待测气体参数. 波长调制-直接吸收光谱(WM-DAS)在DAS基础上, 结合了波长调制光谱(WMS)中谐波分析思想, 利用傅里叶变换复现吸收率函数, 可有效提高吸收率函数的测量精度. 本文利用WM-DAS方法结合长光程气体吸收池, 在室温低压条件下, 对CO分子1567 nm处R5—R11近红外弱吸收谱线吸收率函数进行了精确复现, 其拟合残差标准差低至3 × 10–5, 随后根据测得的吸收率函数对谱线的碰撞展宽、Dicke收敛以及速度依赖的碰撞展宽系数等光谱参数进行了高精度标定, 并将其与高灵敏度的连续波腔衰荡光谱(CW-CRDS)测量结果进行了比较, 实验结果表明该方法与CW-CRDS测量结果具有高度一致性, 更具有系统简单、测量速度快、对环境要求低等优点.

1 引言

气体分子吸收谱线的物理常数(如碰撞展宽系数[1], 线强度[2]等)是测量气体浓度、温度等不可或缺的参数, 它可通过拟合测得的分子吸收率函数得到, 因此, 吸收率函数的测量精度决定了谱线物理常数的标定精度. 目前, 常用的气体分子吸收率函数测量方法主要有傅里叶变换红外吸收光谱[3,4]和可调谐激光二极管吸收光谱(TDLAS)[5−7], 其中TDLAS采用窄带激光扫描分子吸收谱线, 具有波长选择性强、测量速度快、灵敏度高等优点.TDLAS中常用的直接吸收光谱(DAS)[8−11]可直接测量吸收率函数, 通过拟合吸收率函数确定待测气体参数或谱线光谱常数, 其物理概念清晰、操作简单, 在强吸收谱线光谱常数标定和气体参数测量中得到广泛应用.

尽管DAS在气体参数测量和谱线参数标定等方面具有操作简单、物理概念清晰等优点, 但诸如“振动噪声”、“暗噪声”、“光噪声”和“比例噪声”等限制了DAS测量精度的进一步提高[12,13], 难以用于弱吸收条件(如气体浓度低、谱线强度弱等)气体参数或谱线光谱常数的高精度测量. 为此, 科研工作者常采用延长有效吸收光程的方法以提高气体分子对激光的吸收作用[14−19], 进而提高吸收率函数的测量信噪比. 基于几何光学的Herriott[14−17]和White[18,19]等多次反射池通过增加激光反射次数延长吸收光程, 如反射100次即可将10 cm长吸收池的有效吸收光程提高到10 m左右, 在痕量气体监测、弱吸收谱线参数标定等领域得到广泛应用. 与此同时, 为进一步提高吸收率函数的测量精度, 科研工作者进行了诸多研究, 如采用稳定光强的方法以减小光强波动对吸收率函数测量的影响[20],或采用基于物理光学的腔增强吸收光谱[21−23]或连续波腔衰荡光谱(CW-CRDS)[24−26]将有效吸收光程提高到数km以上, 但该类方法对环境要求高,系统操作复杂.

与上述延长吸收光程、稳定光强等研究方法不同, Du等[27]将波长调制光谱(WMS)谐波分析思想引入到DAS中, 通过傅里叶变换中特征频谱复现分子吸收率函数, 即波长调制-直接吸收光谱法(WM-DAS). 该方法融合了 DAS (免标定、可测量吸收率函数)和WMS (高信噪比、抗干扰能力强)的优点, 将吸收率函数测量信噪比提高了约1个量级, 其拟合残差标准差低至约 1 × 10–5量级[28,29].该方法可高精度标定分子谱线的光谱常数, 如碰撞展宽系数、温度指数、线强度等, 即使在弱吸收条件下仍能测得吸收率函数的精细结构, 可实现高精度的Rautian[30]和Galatry[31]线型函数中的Dicke收敛系数、Speed-dependent Voigt[32]线型函数中的速度依赖的碰撞展宽系数等参数的精确测量.

考虑到WM-DAS的优点, 本文利用高信噪比的WM-DAS方法结合长光程Herriott池, 在室温低压条件下, 对CO分子1567 nm处R5—R11近红外弱吸收谱线 (约 2 × 10–23cm–1/(mol·cm–2))的吸收率函数进行了高精度复现, 同时根据测得的吸收率函数标定了谱线碰撞展宽、Dicke收敛和速度依赖的碰撞展宽等系数, 并将测量的吸收率函数、标定的光谱参数及其不确定度与高灵敏度的CW-CRDS测量结果进行了比较.

2 实验系统

实验系统如图1所示, 虚线方框内为WMDAS[27−29]系统, 其他部分与 CW-CRDS[24−26]系统共用. 两系统的气体池通过聚四氟管线联通, 待测气体经过干燥和过滤后进入气体池中, 两气体池内的气体压力、温度、浓度均相同. 两系统共用一个中心波长1567 nm的分布反馈半导体激光器(电流、温度调谐范围分别为 18—100 mA, 10—40 ℃,功率约7 mW), 激光束通过光隔离器以减少对激光器的光反馈. 输出激光分为两束, 一束通过声光调制器进入衰荡腔, 另一束进入Herriott池[14−17],并通过程序控制激光控制器的外调制开关, 以实现WM-DAS和CW-CRDS分时测量. 激光绝对波长和相对波长分别采用波长计(Bristol 671A)和干涉仪(Thorlabs, 自由光谱区30 GHz)进行标定. Herriott池由一对间隔约 1 m、曲率半径 2 m的镀银反射镜(反射率约98%)组成, 总光程约为120 m. 激光光束在Herriott池内经过多次反射后由入射孔出射, 经探测器接收并通过高速采集卡采集.

图1 WM-DAS与 CW-CRDS的系统原理图 (LC, 激光电流 和 温 度 控 制 器 ; FI, 光 纤 隔 离 器 ; AOM, 声 光 调 制 器 ;APD, 雪崩光电二极管; PD, 光电二极管; DDG, 数字延迟发生器; PZT, 压电换能器; WM, 波长计)Fig. 1. System schematic diagram of WM-DAS and CWCRDS. LC, laser current and temperature controller; FI,fiber isolator; AOM, acousto-optic modulator; APD, avalanche photodiode; PD, photodiode; DDG, digital delay generator; PZT, piezoelectric transducer; WM, wavelength meter.

衰荡腔(腔长约0.5 m)由一对高反射率(反射率高于99.99%)镜片组成. 通过PZT扫描腔长,使腔长扫描范围大于一个自由光谱区, 以确保任意波长的激光均可耦合进腔内. 腔另一端的出射光由光电探测器接收, 当探测器达到预设触发电平时,由数字延迟发生器发送脉冲信号, 使声光调制器失去能量, 从而关闭进入腔内的激光以形成单指数衰减信号. 采用高速数据采集卡同时采集脉冲信号和单指数衰减信号, 并利用计算机对实验数据实时处理, 快速拟合[24]得到衰荡时间. 采用步进式扫描激光电流来改变激光波长, 从而得到随波长改变的衰荡时间t(v), 采集多个扫描周期并平均以提高信噪比. 随后, 通过公式将衰荡时间t(v)转换为吸收系数k(v), 其中c为光速,k0表示镜面反射率、散射、吸收等导致的损耗[24−26],再对1/(ct(v))拟合即可得到谱线参数.

3 波长调制-直接吸收光谱

WM-DAS考虑了激光输出光强和波长的非线性, 建立了激光波长与光强之间的关系[27]. 在WM-DAS 方法中, 定义x= cos(wt+h), 其中t为扫描时间,w为正弦调制信号的角频率,h为基倍频初始相位角, 激光光强I和激光波长v可表示为

式中,v0为激光中心波长,a1和a2为调制深度,Ak和Bk是k次特征频率的傅里叶系数的实部和虚部. 根据(1)和(2)式可重构出光强I与波长v的关系, 其中‘－’和‘+’分别代表光强上升沿和下降沿. 根据Beer-Lambert定律可得到吸收率函数如下[27−29]:

式中,It(v)和I0(v)分别为重构的透射光强和入射光强.P(atm)为气体压力 (1 atm = 1.01325 ×105Pa),S(cm–2·atm–1)为谱线强度,T(K)为气体温度,X为气体摩尔分数,L(cm)为光程,j(v)(cm)为线型函数, 可以用 Voigt[27−30], Rautian[31],Galatry[32], Speed-dependent Voigt[33,34]等线型来描述.

其中Voigt线型(VP)可以表示为[27−30]

式 中W(x,y)为 误 差 函 数[33],这 里 ,, 其中gs是待测气体分子的自身碰撞展宽系数,g0是背景气体分子对待测气体分子的碰撞展宽系数,p1是待测气体分压,p2是背景气分压;gD是多普勒半宽,其中M是相对分子质量,T为气体温度. 当背景气体的分子质量大于待测气体分子质量时,j(v)可采用Rautian线型(RP)描述[31]:

与 GP 和 RP 相比, Speed-dependent Voigt线型(SDVP)[33,34]的精度更高, 其表达式为

式中w(z)和z±表达式如下:

4 实验结果与分析

实验中采用 WM-DAS对 CO分子 1567 nm附近的R5—R11谱线进行了测量. Herriott池有效光程120 m, 气体温度、压力和CO浓度分别为288 K, 18 kPa 和 0.1% (背景气 N2). 激光扫描频率、扫描范围分别为 1 kHz, 0.4 cm–1, 单次实验共采集100个正弦波周期(用时0.1 s), 同时采集相应的干涉仪信号(•)进行激光波长标定, 如图2所示, 其中蓝色曲线为透射光强信号, 红色曲线为波长标定结果. 蕴含气体吸收率函数信息的透射光强信号傅里叶系数为Ak和Bk, 将Ak和Bk及通过干涉仪标定的激光波长系数a1,a2和h等参数代入(1)—(3)式中即可复现吸收率函数.

图2 测量的100个正弦波周期的激光光强及激光相对波长标定结果(FSR为自由光谱范围), 以及蕴含气体吸收率函数信息的透射光强傅里叶变换(FFT)系数Fig. 2. Measured transmitted intensities of 100 periods of sinusoidal waves and fitted frequency (FSR, free spectral range), and fast Fourier transform (FFT) coefficients Ak and Bk of transmitted light intensity.

图3 (a)展示了WM-DAS方法对R10谱线的测量结果, 为便于与 CW-CRDS比较, 将 WMDAS所测吸收率(a)转换为吸收系数(k). CWCRDS方法采用步进式扫描激光波长, 间隔约0.002 cm–1, 扫描范围约 0.4 cm–1, 共扫描 100次,用时约20 min. 图3(b)展示了CW-CRDS测量结果, 由于测量的吸收系数仅与衰荡时间相关而与光强无关[24−26], 残差波动小, 测量数据更加平滑, 可以清晰看出VP拟合残差中“w”形的精细结构, 其原因在于VP未考虑Dicke收敛效应[31,32], 观察到该精细结构也说明了CW-CRDS具有很高的灵敏度. 与此相比, WM-DAS通过提取周期性正弦信号的整数倍特征频谱来复现气体吸收光谱, 可有效减小或消除振动、电磁等多种噪声干扰[12,13],VP拟合残差中也可清晰地观察到“w”形的精细结构, 这验证了本文 WM-DAS测量结果与CWCRDS相一致. 与 VP 不同, RP 考虑了 Dicke收敛效应, 拟合时可以消除“w”形残差. 从 RP拟合结果可知, WM-DAS和CW-CRDS两种方法测得的吸收率函数拟合残差标准差相差约2.5倍, 相应的SNR(1s)相差约2.5倍, 但WM-DAS的测量时间 (0.1 s)远小于 CW-CRDS (20 min), 在大量的谱线参数标定以及工业现场的气体快速监测中更有优势.

图3 测量的 CO 光谱及其最佳 Voigt和 Rautian 线型拟合结果(XCO, c, t和sSD分别为CO浓度、光速、衰荡时间和残差的标准差) (a) WM-DAS; (b) CW-CRDSFig. 3. Measured absorption function of CO and the best fits of Voigt and Rautian profile (XCO, c, t, and sSD represent the CO concentration, light velocity, ring down time and standard deviation of the residual, respectively): (a) WMDAS; (b) CW-CRDS.

为了进一步验证WM-DAS的测量精度, 利用该方法对不同压力下CO分子R5—R11谱线的光谱参数进行了测量, 并与CW-CRDS测量结果进行对比, 其中, 碰撞展宽(gc)及速度依赖的碰撞展宽 (g2)选用 SDVP[33,34]的拟合结果, Dicke收敛(b)选用 RP[31,32]的拟合结果. 气体温度为 288 K,温度波动小于 0.3 K, CO浓度为 0.5% (背景气为N2), 压力范围 1—25 kPa. 图4(a)和图4(b)给出了不同压力下, 两种方法对CO分子R10谱线的gc及g2和b测量结果, 测量数据的拟合线性度分别达到 0.9998 (gc)和 0.98 (g2和b), 这说明两方法测量数据有较好的一致性. 在上述参数中,b和g2数值远小于碰撞展宽gc, 对光谱信噪比有较高的要求, 因而其测量误差略大于gc. 事实上, 虽然WM-DAS的 SNR低于CW-CRDS的 SNR, 但是WM-DAS测量速度更快测量点更多, 每个测量点仅需 1 s (信号采集和数据处理), 测量过程受环境温度、气压波动影响更小, 因而b和g2的线性拟合结果与CW-CRDS相一致.

图4 不同压力下测得的光谱参数 (WM-DAS为红色, CW-CRDS为黑色) (a) gc (圆); (b) b (正方形), g2 (三角)Fig. 4. Measured spectral parameters for various pressures (WM-DAS (red), CW-CRDS (black)): (a) gc (dot); (b) b (square), g2(triangle).

表1 WM-DAS和 CW-CRDS 测量的光谱参数及其不确定度Table 1. Measured spectroscopic parameters and uncertainties.

表1列出了WM-DAS和CW-CRDS两种方法采用四种线型 (VP, GP, RP, SDVP)拟合得到的碰撞展宽系数g0(T0), Dicke收敛系数b0(T0)和速度依赖的碰撞展宽系数g2(T0). 测量不确定度主要来源于气体温度 (< 0.2 K)、压力 (< 30 Pa)、浓度 (< 50 ppm)以及吸收率函数拟合误差, 经计算得到g0(T0)总的不确定度范围约为0—1% (CWCRDS)和 0—1.5% (WM-DAS). 其中, 两种方法采用VP拟合的g0(T0)的相对误差均不超过1%,采用精度更高的RP/GP/SDVP拟合的g0(T0)的相对误差仅略大于1%, 主要是由Dicke收敛系数及速度依赖的碰撞展宽系数的拟合不确定引起的.由于b0(T0)和g2(T0)值比g0(T0)小, 对光谱信噪比要求更高, 因此CW-CRDS测得的b0(T0)和g2(T0)的不确定度范围为5%—15%, 与文献[35]的测量结果0.7%—14%相接近. 相比之下, WM-DAS通过提取特征频率复原吸收率函数, 并通过快速测量大量的数据点来降低不确定度,b0(T0)和g2(T0)不确定度范围达到10%—20%.

为了进一步评价WM-DAS和CW-CRDS两种方法的测量下限, 在相同条件下, 分别对WMDAS和CW-CRDS进行了Allan方差[37]分析, 实验结果如图5所示. 其中CW-CRDS采用固定中心波长的测量方式, 积分时间约25 s时(固定波长下的衰荡时间采集速度约为0.01 s, 25 s相当于平均2500次), 基于CW-CRDS的CO检测限可到30 ppb这与文献 [38]检测限相一致, 这验证了本文CW-CRDS测量结果的可靠性. 与此相比, 积分时间约 100 s时, 基于 WM-DAS的 CO 检测限可到 80 ppb, 略高于 CW-CRDS, 对应的吸收系数约 2 × 10–10cm–1.

图5 两种方法测量的 Allan 方差Fig. 5. Allan variance measured by the two methods.

5 结论

本文利用免标定、高信噪比的WM-DAS方法结合长光程Herriott池, 在低压常温条件下, 对CO分子1567 nm附近7条近红外弱吸收谱线(R5—R11)的吸收率函数进行了测量, 光谱拟合残差标准差低至 3 × 10–5. 随后, 采用不同的线型函数对测得的吸收率函数进行拟合得到了谱线光谱参数及其测量不确定度, 并将其与高灵敏度的CW-CRDS测量结果进行比较. 实验结果表明:CW-CRDS方法测得的吸收率函数信噪比约为长光程WM-DAS方法的2.5倍, 两种方法测得的谱线物理参数具有高度一致性, 其中VP线型拟合的碰撞展宽系数的相对误差小于1 %. 与此同时, 基于WM-DAS方法的CO的检测限达到80 ppb, 对应的吸收系数低至 2 × 10–10cm–1, 略高于 CWCRDS方法, 但其测量速度更快, 并且具有系统简单、成本低、可长期稳定运行的优点, 预期可为弱吸收谱线的测量提供新的测量方法.