基于吸收光谱和激光干涉技术的气体压力测量方法研究

张博涵, 杨 军, 黄乾坤, 谢兴娟

航空工业北京长城计量测试技术研究所, 北京 100095

引 言

在气体压力检测方面，目前主要还是利用压力传感器进行定点测量，然而随着激光技术的发展和传感器自身限制日渐明显，越来越多的研究机构将研究重点投入到光学技术在压力非接触测量领域的应用当中，其中包括激光干涉技术和TDLAS技术。

TDLAS技术是利用波长可调谐的激光器扫描出气体特征吸收谱线，再进行相关处理，来实现对气压的测量。Buchholz等根据水汽在7 299.4 cm-1附近吸收峰进行了15～80 kPa的环境气压测量；中国科学院安徽光机所陈祥等利用TDLAS技术进行了30～100 kPa水汽气压测量。激光干涉技术是通过测量得到的干涉条纹信号得到光程差的信号，进而根据压力和光程间的关系得到压力值。美国国家标准与技术研究院(NIST)进行了通过激光干涉法来实现气体压力和温度的测量[1-2]；中国科学技术大学合肥微尺度物质科学国家研究中心的许玉蓉[3]等设计了一种通过法布里-珀罗腔实现氩气折射率精密的测量方法，并进行了气体压力的反演。但无论哪种光学技术，在气体压力测量过程中都要面临温度耦合问题，因此如何解决气体压力温度解耦问题一直是研究工作中的重点。

在以往的工作中对TDLAS技术和激光干涉技术测压进行了相关的研究，为了实现气体压力温度的解耦，此工作内容是以纯CO2为研究对象，进行气体压力组合技术的测量，建立有关积分吸光度和折射率融合的气体压力温度解算模型，这种研究思路可为今后的气压光学测量温度补偿方法提供参考价值。

1 原理分析

TDLAS技术直接吸收法是依据朗伯-比尔(Beer-Lambert)定律[4]，通过积分吸光度实现气体压力的测量，其公式如式(1)

(1)

式(1)中，P是目标气体的压强；I0是入射光的强度；It是吸收后出射光的强度；X是目标气体的浓度；S(T)是线强度函数；L是物理光程；A是吸收谱线的积分吸光度。

其中谱线强度S(T)是表达气体吸收谱线强弱的函数[5]

(2)

式(2)中，T0是参考温度，一般为296 K；T是实验中的温度；Q(T)是气体吸收的配分函数；E″是低能级的能量；ν0是吸收谱线的中心频率。

配分函数Q(T)通常采用的是近似的三次多项式来表达

Q(T)=a+bT+cT2+dT3

(3)

当实验温度低于500 K时，参数a，b，c，d分别取值为-1.361 7，0.948 99，-6.925 9×10-4，2.597 2×10-6。

基于折射率的气体压力激光干涉测量技术依据洛伦兹-洛伦茨(Lorentz-Lorenz)方程

(4)

式(4)中，n为目标气体的折射率；ε0为真空介电常数；N是单位体积内的气体分子数，α为极化率。

又单位体积内气体分子数N可表达为

(5)

式(5)中，ρ为目标气体的密度，Mmol为目标气体的摩尔质量；NA为阿伏伽德罗常数。

假设实验环境满足理想气体状态条件，则有

(6)

式(6)中，ρ为气体密度；R为理想气体常数。

由此，根据式(4)，式(5)，式(6)可得折射率-压力关系模型

(7)

根据以上分析，可以得知基于积分吸光度的压力测量模型和基于折射率的压力测量模型都与气体压力和温度有关。而对于光谱法测压模型中温度参数复杂性，无法简单地进行方程联立求解，因此选择对含有温度变量的谱线强度函数进行三次多项式拟合，考虑到配分函数参数取值条件，选择对CO2气体在273～500 K温度条件下，在4 988.655 cm-1处的吸收谱线强度函数进行拟合。得到的拟合结果如图1所示。

图1 线强度函数拟合结果

拟合得到的方程为

S(T)=-4.058×10-9T3+5.568×10-6T2-

2.6×10-3T+0.4197

(8)

拟合结果R2值为0.9995，证明与原曲线吻合度极高。

因此基于吸收光谱和折射率融合的压力温度解算方程组如式(9)

(9)

式(9)中，a1，b1，c1，d1分别为拟合系数-4.058×10-9，5.568×10-6，-2.6×10-3，0.419 7。

2 实验部分

实验系统如图2，主要分为压力控制单元、光谱测量单元和激光干涉测量单元。

图2 气体压力测量系统示意图

其中压力控制单元采用的是纯度为99.99%的高纯二氧化碳，在室温条件下通过高精度压力控制器(瑞力普 PC7510A)控制实验气室内的压力分别稳定在101.33 kPa(1个大气压)、151.99 kPa(1.5个大气压)、202.65 kPa(2个大气压)、303.98 kPa(3个大气压)四个压力点后进行实验测量。温度检测单元选择的是热电偶式传感器(Autonics TC4SP)进行温度检测。

光谱测量单元采用的是中心波长为2 004 nm的(LD-PDINC，PL-DFB-2004-A-A81-PA)DFB激光器，控制发射光覆盖CO2在4 988.655和4 987.308 cm-1处的吸收峰，出射光经过光电探测器(Thorlabs，PDA10D2)接收，在经过装有NI采集卡(PXI-6115)的PXI设备采集数据。

激光干涉单元采用的是德国SIOS激光干涉仪(SP 2000 DI)，波长为632.8 nm，出射光经过气室，通过反射镜原路返回，搭配相应的上位机软件，采集光程差数据。

3 结果与讨论

3.1 光谱测量结果

利用PXI设备的数据采集卡，采集到的不同压力下CO2吸收谱线原始信号如图3。

图3 不同压力下的吸收信号

从图3中可以看出，在激光扫描周期内有两个吸收波段，通过查询HITRAN数据库和波长计检测并且利用线性插值的方式确定频率横坐标，从而得到具有两个吸收峰的吸光度-波数曲线图。

如图4所示，CO2气体在4 987.308 cm-1处的吸收峰曲线并不完整，为减小因拟合造成的误差，选择气体在4 988.655 cm-1处吸收峰进行洛伦兹函数拟合和积分吸光度值计算，结果如图5。

图4 不同压力下的吸光度曲线

图5 不同压力下的积分吸光度

3.2 光程差测量结果

通过激光干涉仪采集到的不同压力下的光程变化数据如图6(a)所示，然后根据测量得到的光程差数据，并以真空条件下折射率为1，得到不同压力下折射率变化情况

(10)

计算结果如图6(b)所示。从图6中可以看出，随着压力的增大，气体折射率和光程呈正比例关系在增加。

图6 不同压力下的测量结果

3.3 气体压力温度解算结果

根据3.1节和3.2节得到的积分吸光度和折射率值，结合式(9)，得到气体测量过程中的压力值和温度值。其结果如图7和图8所示。

图7 压力测量结果

图8 温度测量结果

从图中可以看出，压力和温度解算最大误差在303.98 kPa测量点，压力测量相对误差为3.61%，温度绝对误差为7.66 K；压力解算最小相对误差在101.33 kPa测量点，为0.50%；温度解算最小绝对误差在202.65 kPa测量点，为0.78 K；该结果验证了理论模型在接近常温环境中的可行性。

4 结 论

通过分析以朗伯-比尔定律为基础的TDLAS直接吸收法气体压力测量原理，和基于折射率的激光干涉测压原理, 提出一种组合技术测量方法，根据光谱技术测量得到的吸光度曲线和激光干涉技术得到的光程差，建立气体吸收积分吸光度和折射率融合后的压力温度测量模型。实验结果表明：在室温下，压力测量平均相对误差为1.99%，温度测量平均绝对误差为3.29 K，测量结果与参考结果具有良好的一致性。进一步探索高温情形下的测量实验将是下一步的主要工作。