基于波数漂移修正算法的免标定固定点波长调制技术 下载: 979次
1 引言
波长调制光谱(WMS)技术[1-4]是一种有效的检测气体参数的方法,具有灵敏度高、低频噪声抑制能力强、无需确定基线等优点,适合弱吸收[5]、强噪声干扰(如燃烧)下气体参数的测量。波长调制光谱技术可通过测量不同浓度的标准气体来制定标定曲线,进而求得气体的浓度[6-7]。但是实际测量环境难以保证与标定环境一致,故导致测量结果不准确。为了克服这一困难,一些研究人员开发出了一系列免标定方法。Rieker等[8]从理论上推导了谐波表达关系式,测量了分布反馈式激光器(DFB)光强和波长的时间响应特性,实现了气体温度、浓度的长时间免标定测量。Sun等[9-10]提出了一种谐波拟合的简化分析模型,利用实测的零气光强结合数字锁相、低通滤波技术来仿真测量信号,实现了高温高压环境下H2O、CO和CO2等气体浓度的测量。
根据发生信号是否叠加扫描信号可以将波长调制光谱技术分为扫描波长调制技术和固定点波长调制技术。扫描波长调制技术极大地拓展了波长调制光谱技术的应用范围,被广泛地应用于气体参数的测量中[11-13]。然而扫描波长法相比较于固定点法,其二次谐波在频谱占据范围较宽,测量结果会受截止频率的影响。固定点法谐波频谱窄,更利于多谱线耦合测量。Faroop等[14]基于固定点波长调制技术开发了一套测温系统,结合光谱数据库仿真得到了两根谱线中心处谐波高度的比值与温度的关系,实现了激波管内温度的准确测量,但这种方法需要已知其他组分的精确的碰撞展宽系数等谱线参数,过分依赖于光谱数据库。Cai等[15]提出了一种准固定点波长调制方法,在调制信号基础上叠加幅值极小的扫描信号,采用峰值点法实现了高压下温度与CO2浓度的准确测量,但其本质属于扫描波长调制技术。在气体参数长时间的测量过程中,激光器在工业恶劣环境下长期运行会面临外界条件的影响,例如环境温度的周期性起伏、便携式设备移动过程中的振动,以及移动过程中带来的环境温度的突变、激光器控制电路的噪声等因素,都会造成激光器内部谐振器有效长度的改变,最终导致激光器出光漂移[16]。商用激光器系统在环境变化(20~60 ℃)时的波长漂移量可达78.1 pm,经过增加额外的温度补偿电路改进温控系统后波长漂移量可减小为4.1 pm[17]。然而,上述研究均未考虑激光器出光中心波数漂移对测量的影响。
目前,针对出光中心波数漂移对测量的影响进行波数漂移修正以实现免标定固定点波长调制技术准确测量的研究未见公开报道。本文重点研究了一种基于波数漂移修正算法的免标定固定点波长调制方法。与前人的研究相比,该方法利用测量得到的展宽参与气体浓度的计算,无须利用光谱数据库中过多的参数,减小了光谱参数不确定性对测量的影响。同时,该方法考虑了激光器出光中心漂移对测量的影响,可以实时计算出光中心,有效提高了气体浓度测量的准确性。选用中心波数为4958.97 cm-1的分布反馈式激光器搭建CO2测量系统,验证了算法的有效性。
2 测量原理
当激光器受到频率为
式中:
其中,
式中:
在波长调制技术中,为了保证较高的测量信噪比,一般采用二次谐波作为分析信号进行气体浓度的测量。在激光器出光的中心波数
其中,
式中:
待测气体的摩尔分数
式中:
3 波数漂移修正算法研究
当激光器出光中心波数发生漂移时,需要实时计算出光中心修正气体浓度的测量结果。在上述理论分析的基础上,提出了基于波数漂移修正算法的免标定固定点波长调制测量方法,该方法的计算流程如
1) 测量得到入射光强
2) 测量得到的标准具信号,提取出激光器时频响应离散点,通过余弦函数拟合得到调制深度
3) 在分析光谱吸收率频谱信息后,对
4) 结合半峰半宽
4 实验分析与讨论
4.1 谱线选择
搭建固定点波长调制CO2测量系统,采用该系统对修正算法进行实验验证。根据HITRAN2016[19]数据库可知CO2在波数4958.97 cm-1(波长为2016.55 nm)附近有两根非常近的吸收谱线(吸收谱线参数见
在常温(296 K)常压(101.325 kPa)条件下,对上述9种气体在所选谱线位置处的光谱吸收率进行仿真,结果如
此外,还需要考虑其他气体组分的干扰。以矿井风道测量环境为例,矿井内的主要气体有O2、CO2、CO、H2S、NO2、SO2、NH3、CH4以及H2O,气体浓度如
表 1. 所选CO2吸收谱线参数(HIRTAN2016)
Table 1. Parameters of selected CO2 spectroscopic absorption lines (HITRAN2016)
|
图 2. 谱线仿真结果。(a) CO2在4958.9674 cm-1附近的吸收谱线; (b) CO2吸收线型的HITRAN数据库仿真结果(仿真条件:x=5.02%, p=101.325 kPa, l=20 cm, T=296 K)及Lorenz线型函数拟合;(c)拟合残差
Fig. 2. Spectral simulation results. (a) CO2 spectroscopic absorption line near 4958.9674 cm-1; (b) CO2 absorption profile simulated using HITRAN database (conditions: x=5.02%, p=101.325 kPa, l=20 cm, T=296 K) with Lorenz-function fitting; (c) fitting residuals (difference between corresponding data points)
表 2. 矿井中各气体组分及其浓度
Table 2. Species and concentrations of gas in mine
|
图 3. 各组分气体在2016.55 nm波长附近的光谱吸收率仿真结果(l=20 cm)
Fig. 3. Simulation of spectral absorbance of gas species near wavelength of 2016.55 nm (l=20 cm)
4.2 实验装置搭建及结果分析
实验装置示意图如
图 4. 固定点波长调制CO2测量系统
Fig. 4. Experimental setup of fixed-point wavelength modulation spectroscopy system for measuring CO2 concentration
在激光器启动过程中,中心温度达到设定值具有一定的时延,出光中心波数随时间发生一定程度的变化。为方便进行研究,利用该过程模拟激光器长期运行时出光中心波数漂移的场景,进而验证修正算法的有效性。实验前用流量为1 L/min的高纯氮气(纯度为99.999%)对气体吸收池进行吹扫,之后切换气路通入流量为1 L/min的CO2标准气体,CO2的摩尔分数为5.02%,待气体充满吸收池后进行实验。记录了激光器从启动到稳定的全过程。由于激光器达到稳定状态的用时较短,取该过程中前15 s的数据进行分析。
假设在调制频率足够高时,调制深度
图 5. 标准具信号与拟合结果。(a)原始信号;(b)激光器频率响应离散点与拟合曲线;(c)拟合残差
Fig. 5. Etalon signal andfitting results. (a) Raw signal; (b) discrete points and fitting curve of laser frequency response; (c) fitting residuals
根据测量得到的透射光强与入射光强可以计算出光谱的吸收率曲线,原始光强如
图 6. 调制系数随时间的变化。(a)调制深度;(b)调制相位
Fig. 6. Variation of modulation coefficient with time. (a) Modulation depth; (b) modulation phase
图 7. 实验信号与光谱吸收率。(a)原始光强信号;(b)去噪前后的光谱吸收率
Fig. 7. Experimental signals and spectral absorbance. (a) Original laser intensity signals; (b) spectral absorbance before and after denoising
如
图 8. 测量结果随时间的变化。(a)出光中心波数;(b)二次谐波幅值
Fig. 8. Variation of measurement result with time. (a) Wavenumber of output laser center; (b) amplitude of second harmonic
CO2气体摩尔分数的测量结果如
图 9. 浓度的测量结果。(a)采用波数漂移修正算法前后计算得到的CO2气体的摩尔分数;(b)测量相对误差
Fig. 9. Concentration measurement results. (a) Calculation results of mole fraction of CO2 before and after wavelength drift-correction algorithm; (b) relative error of measurement
5 结论
针对激光器出光中心波数漂移导致测量不准确的问题,提出了波数漂移修正算法,实现了免标定固定点波长调制的准确测量。搭建实验装置对CO2气体浓度进行测量,测量结果证明了算法的可行性。与前人的研究相比,所提算法通过测量获取吸收谱线的展宽,在计算气体浓度过程中无需其他组分的碰撞展宽系数、温度依赖指数等参数,减小了对光谱数据库的依赖性,避免了光谱参数不确定性对测量的影响。所提算法通过实时计算激光器出光中心波数,克服了激光器出光中心波数漂移对测量的不良影响,具有更高的测量准确度。
[1] 陈家金, 王贵师, 刘锟, 等. 免标定波长调制吸收光谱技术用于乙炔探测的研究[J]. 光学学报, 2018, 38(9): 0930001.
[3] 熊涌泉, 周宾, 王一红, 等. 基于波长调制光谱技术的免标定单线测量法[J]. 光学学报, 2017, 37(9): 0930002.
[4] 汪步斌, 周宾, 王浩, 等. 基于三角波调制的免标定气体浓度测量方法研究[J]. 光学学报, 2017, 37(9): 0930003.
[5] 臧益鹏, 聂伟, 许振宇, 等. 基于可调谐二极管激光吸收光谱的痕量水汽测量[J]. 光学学报, 2018, 38(11): 1130004.
[6] Liang W K, Bi Y F, Zhou Q, et al. Developing CH4 detection limit at λ=1.654 μm by suppressing optical interference fringes in wavelength modulation spectroscopy[J]. Sensors and Actuators B: Chemical, 2018, 255: 2614-2620.
[7] Li J S, Reiffs A, Parchatka U, et al. In situ measurements of atmospheric CO and its correlation with NOx and O3 at a rural mountain site[J]. Metrology and Measurement Systems, 2015, 22(1): 25-38.
[16] 唐七星, 张玉钧, 陈东, 等. 调谐激光吸收光谱波长偏移修正算法研究[J]. 光谱学与光谱分析, 2018, 38(11): 3328-3333.
[17] 袁松, 阚瑞峰, 何亚柏, 等. 可调谐半导体激光吸收光谱中激光器温度补偿[J]. 中国激光, 2013, 40(5): 0515002.
[20] 张浩, 赵振洋, 周松. 15):[J]. . 论矿井空气成分及含量. 科技致富向导, 2012, 84: 172.
Zhang H, Zhao Z Y. 15):[J]. Zhou S. Discussion on mine air composition, content . Guide of Sci-tech Magazine, 2012, 84: 172.
Article Outline
祝仰坤, 周宾, 王一红, 连久翔, 王汗青, 陈慈伟. 基于波数漂移修正算法的免标定固定点波长调制技术[J]. 光学学报, 2019, 39(8): 0830001. Yangkun Zhu, Bin Zhou, Yihong Wang, Jiuxiang Lian, Hanqing Wang, Ciwei Chen. Research on Calibration-Free Fixed-Point Wavelength Modulation Spectroscopy Based on Wavenumber Drift-Correction Algorithm[J]. Acta Optica Sinica, 2019, 39(8): 0830001.