TDLAS波长调制法中调制深度与高次谐波中心幅值关系的研究

引用本文

陈昊, 鞠昱, 韩立. TDLAS波长调制法中调制深度与高次谐波中心幅值关系的研究. 光谱学与光谱分析, 2021,41(12): 3676-3681
CHEN Hao, JU Yu, HAN Li. Research on the Relationship Between Modulation Depth and Center of High Order Harmonic in TDLAS Wavelength Modulation Method. Spectroscopy and Spectral Analysis, 2021,41(12): 3676-3681. 复制到剪切板

Doi: 10.3964/j.issn.1000-0593(2021)12-3676-06
Permissions

《光谱学与光谱分析》期刊社所有

TDLAS波长调制法中调制深度与高次谐波中心幅值关系的研究

陈昊^1,², 鞠昱³, 韩立¹

1.中国科学院电工研究所, 北京 100190

2.中国科学院大学, 北京 100049

3.北京航天易联科技发展有限公司, 北京 100176

作者简介: 陈昊, 1990年生, 中国科学院电工研究所博士研究生 e-mail: chenhao9@mail.iee.ac.cn

收稿日期: 2020-11-27 修订日期: 2021-02-19

基金项目: 国家自然科学基金项目(41705141)资助

摘要

可调谐半导体激光光谱技术(TDLAS)是近年来发展十分迅速的光谱检测技术, 相较于其他光谱检测技术, 它具有高灵敏度、高分辨率、实时监测、便携性好、小型化等优点, 在工业环保、医疗检测、气象监测等领域得到了广泛的应用。 TDLAS波长调制法中谐波信号易受气压影响, 经研究发现气压的影响是调制深度对谐波信号的影响, 基于TDLAS技术谐波法的原理, 研究了各次谐波与调制深度的关系, 通过计算四次谐波与二次谐波中心幅值比, 利用调制深度函数推算当前气压环境的调制深度, 调整调制频率幅度, 使得调制深度接近各次谐波最佳调制深度值, 使谐波信号信噪比最佳, 提高检测精度。实验通过国瑞智GRZ5031湿度发生器产生固定为1 000 ppm的水汽, 调节气阀控制密封箱内不同的气压环境, 采用TDLAS水汽检测系统获得了10.2~177.9 kPa气压条件下的二次谐波和四次谐波信号, 并进行了仿真与实验分析。仿真结果显示: 四次谐波与二次谐波中心幅值比的理论值和仿真值最大相对误差为-1.44%, 调制深度的理论值与仿真值最大相对误差为1.78%, 说明了仿真下的调制深度函数曲线与理论一致。实验结果显示: 根据调制深度函数推算调制深度值, 当 m=2.226 7时, 实测的二次谐波中心频率幅值达到最大值, 当 m=4.061 0时, 实测的四次谐波中心频率幅值达到最大值, 与理论结果一致; 在30.2 kPa< p<177.9 kPa时, 调制深度与气压乘积mp值相对误差较小, 最大相对误差不超过±3.2%, 说明了此气压条件下的 mp值波动不大, 通过调制深度函数推算的调制深度值与实际值近似, 验证了调制深度函数理论的准确性。

关键词: TDLAS; 气压; 调制深度; 高次谐波

中图分类号:O433.41 文献标识码:A

Research on the Relationship Between Modulation Depth and Center of High Order Harmonic in TDLAS Wavelength Modulation Method

CHEN Hao^1,², JU Yu³, HAN Li¹

1. Institute of Electrical Engineering, Chinese Academy of Sciences, Beijing 100190, China

2. University of Chinese Academy of Sciences, Beijing 100049, China

3. Beijing Aerospace Yilian Science and Technology Development Company, Beijing 100176, China

Abstract

Tunable Semiconductor Laser Spectroscopy (TDLAS) has been a rapidly developing spectral detection technology in recent years. Compared with other spectral detection technologies, TDLAS has the advantages of high sensitivity, high resolution, real-time monitoring, good portability and miniaturization, and has been widely used in the fields of industrial environmental protection, medical detection, meteorological monitoring and so on. The harmonic signal in the TDLAS wavelength modulation method is susceptible to air pressure. It is found that the influence of air pressure is the influence of modulation depth on the harmonic signal. Based on the principle of the harmonic method in TDLAS technology, the relationship between each harmonic and modulation depth is studied. The modulation depth of the current air pressure environment is calculated by calculating the center amplitude ratio of the fourth harmonic to the second harmonic. To adjust the amplitude of the modulation frequency so that the modulation depth is close to the optimal modulation depth value of each subharmonic, and the signal-to-noise ratio of the harmonic signal is optimal to improve the detection accuracy. The second and fourth harmonic signals under 10.2~177.9 kPa pressure were extracted by the TDLAS water vapor detection system. The simulation and experimental analysis were carried out. The simulation results show that the maximum relative error of the center amplitude ratio of the fourth harmonic to the second harmonic is -1.44%, and the maximum relative error of the modulation depth simulation to the theoretical value is 1.78%. The experimental results show that the modulation depth value is calculated based on the modulation depth function. When m=2.226 7, the measured central frequency amplitude of the second harmonic reaches the maximum value. When m=4.061 0, the measured central frequency amplitude of the fourth harmonic reaches the maximum value, which is consistent with the theoretical results. When 30.2 kPa< p<177.9 kPa, the relative error between the modulation depth and the pressure product MP value is small. The maximum relative error is not more than ±3.2%. It shows that the MP value under this pressure does not fluctuate much. The modulation depth value calculated by the modulation depth function approximates the actual value, which verifies the accuracy of the modulation depth function theory.

Key words: TDLAS; Pressure; Modulation depth; High harmonic

文章图片

引言

可调谐半导体激光光谱技术(tunable diode laser absorption spectroscopy, TDLAS)是利用气体分子对激光选频吸收, 计算入射光与出射光的光功率变化, 实现对待测气体浓度的定量检测。近年来大量学者对TDLAS技术进行了研究, 相较于其他光谱检测技术, 它具有高灵敏度、高分辨率、实时监测、便携性好、小型化等优点, 在工业环保、医疗检测、气象监测等领域得到了广泛的应用^{[1, 2, 3]}。

当气体浓度较低或气体吸收峰较弱的情况下, 吸收产生的光强变化微弱, 直接法(direct absorption spectroscopy, DAS)测量无法获得理想的信噪比, 需要采用谐波法(wavelength modulation spectroscopy, WMS)进行检测, 它是利用锁相放大器解调特定频率的谐波进行检测, 计算各次谐波峰值或者峰值比表征气体浓度, 具有更低的检测下限和灵敏度^{[4, 5]}。各次谐波信号与调制深度m相关, 当m≈ 2.2时, 二次谐波的峰值达到最大; 当m≈ 4.1时, 四次谐波的峰值达到最大^[6]。当气压发生剧烈变化时, 气体半高宽随气压变化^[7], 而调制深度为调制频率幅度与气体半高宽的比值, 因此调制深度也受气压影响, 导致各次谐波的幅值和信噪比都受到影响, 测量结果容易产生误差。

如何确定当前气压环境下的调制深度值就显得至关重要, 基于此本文研究了TDLAS技术中谐波法的原理, 推导了在中心频率位置下的一次谐波、二次谐波、三次谐波以及四次谐波幅值强度表达式, 四次谐波与二次谐波中心幅值强度比只与调制深度值相关, 建立了调制深度函数, 通过调制深度函数来推算当前环境下的调制深度值, 再调整调制频率幅度, 使调制深度m尽可能接近2.2, 来保证二次谐波信号的幅值保持在最大值, 提高信噪比和测量精度。

1 TDLAS调制深度函数

TDLAS技术理论基础是Beer-Lambert定律, 它描述了一束特定频率的激光进入气体样品前后的光强变化, 如式(1)所示。

$I_{t} = I_{0} \exp [- α (ν) cL]$ (1)

式(1)中, I_t为穿过待测气体后的透射光光强, I₀为入射光强, α (ν )为吸收系数, c为待测气体的浓度, L为光吸收路径长度。

对于谐波法, 光源的驱动电流叠加了一个低频的小信号, 激光器的出光频率和光强都会产生调制。

$ν = \overset{̅}{ν} + a \sin ωt$ (2)

$I_{0} = {\overset{̅}{I}}_{0} + Δ i \sin ωt$ (3)

式(2)和式(3)中, ν 为扫描频率, $\overset{̅}{ν}$ 为中心扫描频率, a为调制频率幅度, ${\overset{̅}{I}}_{0}$ 为中心扫描频率下的光强, Δ i为光强调制放大系数, ω 为调制频率。采用谐波法对弱信号进行检测, 则α (ν )cL≪1, 因此式(1)可以近似如式(4)

$I_{t} = ({\overset{̅}{I}}_{0} + Δ i \sin ω t) (1 - α (ν) cL)$ (4)

当采用洛伦兹线型表示吸收系数α (ν )时

$\begin{array}{l} α (ν) = SNf (ν) = SN \frac{γ}{π ((\overset{̅}{ν} + a \sin ωt - ν_{0})^{2} + γ^{2})} \\ = \frac{SN}{π γ} \frac{1}{({(\frac{\overset{̅}{ν} + α \sin ω t - ν_{0}}{γ})}^{2} + 1)} \end{array}$ (5)

式(5)中, ν ₀为气体吸收峰中心频率, S为气体吸收线强, N为标准状况下(STP)的气体摩尔分子数, γ 为待测气体吸收谱半高宽, f(ν )为洛伦兹线型函数, 调制深度m=a/γ , 令a₀=SN/π γ , 即气体吸收峰中心位置(ν =ν ₀)的吸收系数, 则式(5)简化为

$α (ν) = a_{0} \frac{1}{({(\frac{\overset{̅}{ν} - ν_{0}}{γ} + m \sin ω t)}^{2} + 1)}$ (6)

对于气体吸收峰中心位置, 即ν =ν ₀, 式(6)代入式(4)得

$\begin{array}{l} I_{t} = ({\overset{̅}{I}}_{0} + Δ i \sin ω t) (1 - a_{0} \frac{1}{({(m \sin ω t)}^{2} + 1)} cL) \\ = {\overset{̅}{I}}_{0} + Δ i \sin ω t - \frac{{\overset{̅}{I}}_{0}}{({(m \sin ω t)}^{2} + 1)} a_{0} cL - \\ \frac{Δ i \sin ω t}{({(m \sin ωt)}^{2} + 1)} a_{0} cL \end{array}$ (7)

对式(7)按照ω 进行傅里叶级数展开, 得到气体吸收峰中心位置的各次谐波强度, 其中一次谐波到四次谐波幅值强度表达式如式(8)— 式(11)

$I_{1 f} = Δ i - \frac{2 \sqrt[]{1 + m^{2}} - 2}{m^{2} \sqrt[]{1 + m^{2}}} Δ i a_{0} cL$ (8)

$I_{2 f} = - \frac{2 (2 + m^{2} + 2 \sqrt[]{1 + m^{2}}}{m^{2} \sqrt[]{1 + m^{2}}} {\overset{̅}{I}}_{0} a_{0} cL$ (9)

$I_{3 f} = \frac{2 (4 + m^{2} - \frac{1}{\sqrt[]{1 + m^{2}}} - 3 \sqrt[]{1 + m^{2}})}{m^{4}} Δ i a_{0} cL$ (10)

$I_{4 f} = - \frac{2 (8 + m^{4} - 8 \sqrt[]{1 + m^{2}} - 4 m^{2} (- 2 + \sqrt[]{1 + m^{2}}))}{m^{4} \sqrt[]{1 + m^{2}}} \cdot {\overset{̅}{I}}_{0} a_{0} cL$ (11)

利用四次谐波式(11)与二次谐波式(9)幅值强度进行比值计算, 得到一个关于调制深度m的表达式。

$\frac{I_{4 f}}{I_{2 f}} = f (m) = \frac{8 + m^{4} - 8 \sqrt[]{1 + m^{2}} - 4 m^{2} (- 2 + \sqrt[]{1 + m^{2}})}{m^{2} (2 + m^{2} - 2 \sqrt[]{1 + m^{2}})}$ (12)

中心频率位置下的四次谐波与二次谐波比值仅与调制深度m相关, 与浓度、温度、入射光强都不相关, 根据式(12)可以得到当前气压环境下的调制深度值m。待测气体浓度的表征方式是以二次谐波峰值(即吸收峰中心位置, ν =ν ₀)与标准浓度值进行拟合表示, 根据式(9)可知, 当m≈ 2.2时, 二次谐波幅值强度达到最大, 信噪比最好, 通过计算四次谐波与二次谐波的比值得到当前气压环境下的调制深度m, 调整调制频率幅度a, 让调制深度m尽可能接近2.2, 以此获得幅值最强、信噪比最好的二次谐波信号, 提高测量准确度。

2 仿真分析

2.1 调制深度函数曲线

根据式(9)、式(11)以及式(12)可以得到在中心频率处的, 二次谐波中心幅值I₂_f(m)、四次谐波中心幅值I₄_f(m)以及两者幅值强度比值I₄_f/I₂_f(调制深度函数f(m))关于调制深度m的函数曲线, 如图1所示。

	Figure Option View Download New Window
	图1 二次谐波中心幅值、四次谐波中心幅值以及调制深度函数曲线Fig.1 Second harmonic center amplitude, fourth harmonic center amplitude and modulation depth function curve

图1(a)是二次谐波幅值I₂_f随调制深度m变化的曲线, I₂_f(m)不是单调的, 当m≈ 2.2时, 二次谐波幅值达到最大值; 图1(b)是四次谐波幅值I₄_f随调制深度m变化的曲线, I₄_f(m)也不是单调的, 当m≈ 4.1时, 四次谐波幅值达到最大值; 图1(c)是调制深度函数曲线f(m), 它是单调递增的函数。实测中为了获得更好的二次谐波信噪比, 应该尽可能将调制深度控制在2.2附近, 而I₂_f(m)与I₄_f(m)都是非单调的, 存在一个最大值, 通过I₂_f(m)与I₄_f(m)的曲线变化来判断当前调制深度就显得尤为困难, 但是调制深度函数f(m)是单调递增, 并且只与调制深度m相关, 因此根据调制深度函数f(m), 计算四次谐波与二次谐波幅值强度比值来反推当前环境下的调制深度是可行的, 并具有唯一性。

2.2 仿真验证

为了验证调制深度函数f(m)反推调制深度的可行性, 选择水汽中心频率为7 306.752 1 cm^-1作为吸收峰, 吸收强度为1.8× 10^-20 cm^-1· (molec· cm^-2)^-1, 气体浓度设定为1 000 ppm, 吸收光程为30 cm, 模拟了气压从0.3~1.5 atm时二次谐波以及四次谐波信号, 并对其进行了分析。

图2(a)是气压0.3~1.5 atm时的二次谐波模拟信号, 二次谐波中心幅值在p=0.9 atm时达到最大, 此时调制深度理论值为2.378 8, 二次谐波中心幅值最大值位置与图1中I₂_f(m)曲线结果相吻合; 图2(b)是气压0.3~1.5 atm时的四次谐波模拟信号, 四次谐波中心幅值在p=0.5 atm时达到最大, 此时调制深度理论值为4.281 9, 四次谐波中心幅值最大值位置与图1中I₄_f(m)曲线结果相吻合。表1计算了各个气压下的调制深度m和四次谐波与二次谐波中心幅值比I₄_f/I₂_f的理论值和仿真值, 随着气压的升高, 调制深度理论值逐渐减小, 四次谐波与二次谐波中心幅值比的理论值和仿真值最大相对误差为-1.44%, 反推的调制深度的理论值与仿真值最大相对误差为1.78%, 说明了仿真下的实际调制深度函数f(m)与理论函数式(12)接近, 符合理论推导。

	Figure Option View Download New Window
	图2 不同气压下二次谐波与四次谐波的仿真信号Fig.2 Simulation signals of the second and fourth harmonics under different air pressures

表1 不同气压下调制深度和四次谐波二次谐波中心幅值比的理论值和实际值 Table 1 Theoretical and practical values of modulation depth m and harmonic center amplitude ratio under different air pressure

3 实验部分

仿真结果验证了调制深度函数f(m)的仿真结果与理论值一致, 下面进一步用实验进行论证。实验装置示意图如图3所示, TDLAS水汽测试系统选用中心频率为7 306.79 cm^-1的DFB半导体激光器作为光源, 气室探头的吸收光程为30 cm, 采用锁相放大的方式对二次谐波和四次谐波信号进行提取。国瑞智GRZ5013湿度发生器产生低浓度水汽气体, 浓度设定为1 000 ppm, 自制密封箱内放入气室探头并接入压力传感器, 做好密封处理, 实验前将密封箱内气体抽空, 随后将1 000 ppm的水汽气体通过气体阀门缓慢通入密封箱内, 通过控制阀门来调节密封箱内的气压, 记录不同气压条件下的二次谐波与四次谐波信号。

	Figure Option View Download New Window
	图3 实验装置示意图Fig.3 Schematic diagram of experimental device

4 结果与讨论

根据所搭建的实验装置, 通过国瑞智GRZ5031湿度发生器浓度产生1 000 ppm的水汽, 调节减压阀门控制密封箱的气压范围在10.4~177.9 kPa, 相敏检波器输入二倍频和四倍频的参考信号, 调整滤波器和放大器参数, 得到了不同气压下的二次谐波和四次谐波信号, 如图4所示。

	Figure Option View Download New Window
	图4 不同气压下的二次谐波和四次谐波实测信号Fig.4 Measured second and fourth harmonic signals at different air pressures

图4(a)是不同气压下的二次谐波信号, 图4(b)是不同气压下的四次谐波信号, 计算了中心频率下的二次谐波和四次谐波信号的幅值以及幅值比I₄_f/I₂_f, 利用式(12)调制深度函数f(m)推算了各个气压下的调制深度值, 如表2所示。

表2 不同气压下的二次谐波与四次谐波的中心幅值、幅值比以及调制深度 Table 2 Central amplitude, amplitude ratio and modulation depth m of the second and fourth harmonics under different air pressures

根据表2中的中心频率下的幅值比以及对应的调制深度m值, 可以得到二次谐波I₂_f(m)和四次谐波I₄_f(m)曲线, 如图5所示。当m=2.226 7时, 实测的二次谐波中心频率幅值达到最大值, 而m=4.061 0时, 实测的四次谐波中心频率幅值达到最大值, 这与图1理论结果的m≈ 2.2时二次谐波幅值达到最大值以及m≈ 4.1时四次谐波幅值达到最大值有很好的吻合。图5中实测的二次谐波I₂_f(m)和四次谐波I₄_f(m)的函数曲线与图1的理论曲线变化趋势相同, 说明了通过调制深度函数f(m)推算的调制深度值m接近实际环境下的调制深度值, 并未存在较大误差, 论证了理论的准确性。

	Figure Option View Download New Window
	图5 实测的二次谐波和四次谐波中心幅值曲线Fig.5 Curves of measured second harmonic and fourth harmonic center amplitude

为了进一步验证理论的准确性, 计算调制深度m与气压p的乘积值可以来判断结果的准确性, 气体半高宽γ 表达式如式(13)所示

$γ (p, T) = \frac{p}{p_{0}} {(\frac{T_{0}}{T})}^{n} [(1 - c) γ_{air} (p_{0}, T_{0}) + c γ_{self} (p_{0}, T_{0})]$ (13)

式(13)中, n为温度系数, p₀和T₀分别为标准气压和标准温度, γ _air是空气吸收谱线半高宽, γ _self是待测气体吸收谱线半高宽。根据式(13)可知气体半高宽γ 在温度和气体浓度不变的情况下与气压p是正相关的, 即γ =ε p, ε 为正相关系数, 而调制深度m=a/γ , 调制频率幅度a=mγ =ε mp。实验中湿度发生器产生固定浓度为1 000 ppm的水汽并且保持密封箱内温度恒定, 因此实测的调制深度m与气压p的乘积应该为定值, 即mp=a/ε , 因此计算不同气压下mp值可以验证推算的调制深度是否准确。图6为不同气压条件下关于mp平均值的相对误差曲线。

	Figure Option View Download New Window
	图6 不同气压条件下关于mp平均值的相对误差曲线Fig.6 Relative error curve of mean value of mp under different pressure conditions

根据图6所示, mp的平均值为109.46, 各个气压条件下的mp值的相对误差显示, 在10.4 kPa< p< 30.2 kPa时, 相对误差较大, 最大相对误差为-18%, 低压下相对误差较大的原因可能是真空抽气时并不能完全将密封箱内的气体抽完, 开始通气时通入的气体有一部分是残留在气管内的气体, 造成了一定的干扰; 在30.2 kPa< p< 177.9 kPa时, 相对误差较小, 最大相对误差不超过± 3.2%, 说明了各个气压条件下的mp值变化不大, 通过调制深度函数f(m)推算的调制深度值m与实际环境下的值接近, 验证了深度函数f(m)理论的准确性。

5 结论

TDLAS波长调制法中谐波信号受气压影响较大, 经研究发现气压的影响是调制深度对谐波信号的影响, 本文研究了各次谐波与调制深度的关系, 通过计算四次谐波与二次谐波中心幅值比, 利用调制深度函数推算当前气压环境的调制深度, 调整调制频率幅度, 使得调制深度接近各次谐波最佳调制深度值, 谐波信号信噪比最佳, 为验证结论进行了仿真与实验。仿真结果显示: 四次谐波与二次谐波中心幅值比的理论值和仿真值的最大相对误差为-1.44%, 调制深度理论值与仿真值的最大相对误差为1.78%, 说明了仿真下的调制深度函数曲线与理论函数式(12)接近, 符合理论推导。实验结果显示: 根据调制深度函数推算的调制深度值, 当m=2.226 7时, 实测的二次谐波中心频率幅值达到最大值, 而m=4.061 0时, 实测的四次谐波中心频率幅值达到最大值, 与仿真结果一致; 在30.2 kPa< p< 177.9 kPa时, 调制深度与气压乘积mp值相对误差较小, 最大相对误差不超过± 3.2%, 说明了此气压条件下的mp值波动不大, 通过调制深度函数f(m)推算的调制深度值m与实际值接近, 验证了调制深度函数f(m)理论的准确性。

参考文献

文献选项

[1]	CHEN Hao, JU Yu, HAN Li, et al(陈昊, 鞠昱, 韩立, 等). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2020, 40(10): 3015. [本文引用:1]
[2]	LI Ming-xing, CHEN Bing, RUAN Jun, et al(李明星, 陈兵, 阮俊, 等). Optics and Precision Engineering(光学精密工程), 2020, 28(7): 1424. [本文引用:1]
[3]	Lee Jungwun, Bong Cheolwoo, Yoo Jihyung, et al. Optic Express, 2020, 28(14): 21121. [本文引用:1]
[4]	YAO Lu, LIU Wen-qing, LIU Jian-guo, et al(姚路, 刘文清, 刘建国, 等). Chinese Journal of Lasers(中国激光), 2015, 42(2): 0215003. [本文引用:1]
[5]	LI Min, BO Feng-ming(李珉, 柏逢明). Laser Journal(激光杂志), 2015, 36(10): 75. [本文引用:1]
[6]	QU Dong-sheng, HONG Yan-ji, WANG Guang-yu, et al(屈东胜, 洪延姬, 王广宇, 等). Journal of Infrared and Millimeter Waves(红外与毫米波学报), 2016, 35(4): 470. [本文引用:1]
[7]	ZHANG Ke-ke, LIU Shi-xuan, CHEN Shi-zhe, et al(张可可, 刘世萱, 陈世哲, 等). Instrument Technique and Sensor(仪表技术与传感器), 2016, (1): 53. [本文引用:1]

2020

0.0

... 近年来大量学者对TDLAS技术进行了研究, 相较于其他光谱检测技术, 它具有高灵敏度、高分辨率、实时监测、便携性好、小型化等优点, 在工业环保、医疗检测、气象监测等领域得到了广泛的应用^[1,2,3] ...

2020

0.0

2020

0.0

2015

0.0

... 当气体浓度较低或气体吸收峰较弱的情况下, 吸收产生的光强变化微弱, 直接法(direct absorption spectroscopy, DAS)测量无法获得理想的信噪比, 需要采用谐波法(wavelength modulation spectroscopy, WMS)进行检测, 它是利用锁相放大器解调特定频率的谐波进行检测, 计算各次谐波峰值或者峰值比表征气体浓度, 具有更低的检测下限和灵敏度^[4,5] ...

2015

0.0

2016

0.0

... 1时, 四次谐波的峰值达到最大^[6] ...

2016

0.0

... 当气压发生剧烈变化时, 气体半高宽随气压变化^[7], 而调制深度为调制频率幅度与气体半高宽的比值, 因此调制深度也受气压影响, 导致各次谐波的幅值和信噪比都受到影响, 测量结果容易产生误差 ...