中国激光, 2019, 46 (7): 0711001, 网络出版: 2019-07-11   

燃烧流场波长调制光谱吸收模型的研究 下载: 989次

Absorption Model of Wavelength Modulation Spectroscopy in Combustion Flow Field
作者单位
1 中国科学院合肥物质科学研究院安徽光学精密机械研究所环境光学与技术重点实验室, 安徽 合肥 230031
2 中国科学技术大学, 安徽 合肥 230026
摘要
基于可调谐半导体吸收光谱的波长调制技术,建立了精确的吸收模型。通过两条已知吸收中心的吸收谱线,对标准具自由光谱范围进行标定,并利用更贴近激光器出光特性的描述模型,得到激光器频率-时间响应,结合实验室标定和HITEMP数据库的杂合吸收谱线参数,建立了可与实际吸收直接比较的精确模型,以诊断燃烧流场。本研究以H2O为目标分子,选取吸收中心为7185.60 cm -1和6807.83 cm -1两条吸收线,利用扣除背景的归一化二次谐波信号峰值反演流场温度,并在管式高温炉上进行实验验证,最高测量温度为1500 K,相对误差小于3.1%。吸收模型的准确性决定了所测流场参数的准确性,该模型可应用到更为复杂的燃烧流场环境,实现流场参数的精确测量。
Abstract
In this study, an accurate absorption model is established by using the wavelength modulation technology based on tunable diode laser absorption spectroscopy. The free spectral range of an etalon is calibrated by the two absorption lines whose transition line centers are known, and the frequency-time response of the laser is obtained using a description model close to the output characteristics of the laser. Along with the hybrid absorption line parameters calibrated in the laboratory and the HITEMP database, an accurate model that can be directly compared with the actual absorption is established for diagnosing the combustion flow field. Herein, H2O is the target molecule, and the transition line centers of the two selected lines are 7185.60 and 6807.83 cm -1. The flow field temperature can be obtained via the peak of the normalized second harmonic signal with background subtracted and verified in a high-temperature flow field generated by a tubular high-temperature furnace. The maximum temperature is 1500 K, and the relative measurement error is less than 3.1%. Further, the accuracy of the absorption model determines the accuracy of the measured flow field parameters. The method used to establish the model can be applied to the complex combustion flow fields for accurately measuring the flow field parameters.

1 引言

可调谐半导体吸收光谱(TDLAS)技术是一种非入侵式光学测量技术,主要分为直接吸收(DA)和波长调制(WMS)两种。该技术可以实现原位测量,具有高灵敏度、高分辨、高选择、信号高保真的特点[1-3],能够对温度、组分浓度、流速等多流场参数进行同时测量[4-5]。波长调制技术的本质是将吸收信息转移到高频部分,有效抑制低频噪声,不需要基线,并且对非吸收损耗不敏感,即解调得到的谐波信号对原始强度信号变弱不敏感。该技术已成为流场诊断的理想方法之一[6-7]

波长调制技术的研究以Wilson等[8]采用数值仿真方法获得谐波信号开始,Philippe等[9]在二极管激光器中注入正弦调制电流,发展了典型的波长调制技术,当时的测量需要在固定环境中事先标定。Sun[10]和Goldstein等[11]利用被一次谐波归一化的二次谐波信号来消除激光器本身的影响,实现了流场参数的测量。水汽是燃烧的主要产物之一,水分子在红外波段具有丰富的振-转光谱,从而成为燃烧流场诊断的理想目标分子之一。本文以H2O为目标分子,重点研究吸收模型的建立、频率的准确标定、线型函数的选取、谱线参数的选取等对所建吸收模型的影响,利用建立的吸收模型反演得到了管式高温炉温度,并从实验上验证了该模型的准确性。

2 基本原理

当一束频率为υ的激光通过待测气体时,其出射光强和入射光强满足朗伯-比尔定律, 即

T[υ(t)]=ItIoυ=exp(-αυ),(1)

式中:t为时间;T[υ(t)]为透过率;It为透射光强;Io为入射光强;αυ为吸光度。

吸收的本质是光子的能量与跃迁的能级差相等,气体分子的吸收并不只针对单一频率。展宽机制使光谱加宽,加宽分为非均匀加宽和均匀加宽,两种加宽机制的半峰全宽可表示为

ΔυD=7.1632×10-7υ0TM,(2)

ΔυC=2p·χQ2γP-Q,3

式中:υ0为跃迁中心频率,cm-1;T为温度,K;M为分子质量,g·mol-1;p为压强,atm,1 atm=101.325 kPa;χQ为组分Q的体积分数;γP-Q为组分PQ间的展宽系数。

实现波长调制技术常用的方法是在可调谐半导体激光器中注入低频锯齿扫描信号和高频正弦调制信号,实现对激光器出光频率和强度的调制。加入调制后的激光频率υ(t)和激光器的出光强度I0(t)可分别表示为

υ(t)=υ-(t)+acos(ωt),(4)

I0t=I-0t1+i0coscosωt+φ1+i2coscos2ωt+φ2,5

式中: υ-(t)为激光器的中心频率;ω为角频率(ω=f,f为频率);a为调制深度,cm-1; I-0(t)为激光器的中心光强;i0i2分别为归一化一次和二次调制幅度;φ1φ2分别为一次和二次调制的初相位;频率表达式中初相位取0。这些参数均依赖于激光器本身的调制特性。由于υ(t)为时间t的偶函数,因而透过率T[υ(t)]也为t的偶函数,透过率的傅里叶展开形式可以表示为

T[υ(t)]=k=0+Hkυ-(t),acos(kωt),(6)

式中:Hkk阶傅里叶展开系数。透射强度的表达式为

I(t)=I0(t)·T[υ(t)]=k=0+[Xkcos(kωt)+Yksin(kωt)],(7)

式中:XkYk分别为XY通道的k阶分量。

经过解调,得到扣除背景的1f归一化的2f信号,表达式为[10]

C2f/1f_bgsub=X2ftR1f-X2f0tR1f02+Y2ftR1f-Y2f0tR1f02,(8)

式中:X2fY2f为经锁相放大器解调得到的2fX分量和Y分量信号; X2f0Y2f0为2fX分量和Y分量的背景信号;R1fR1f0分别为双通道解调得到的1f信号和背景。

3 模型建立

建立吸收模型,即根据光谱参数及初始流场条件模拟得到谐波信号,并根据流场参数对谐波信号的依赖反演得到流场参数。模型建立流程如图1所示。波长调制技术反演流场参数主要有全线型拟合方法以及谐波信号峰值或峰值比方法。这两种方法均要求模拟的流场条件尽可能与实际条件对应,谐波信号与实测信号尽可能一致,故要求吸收模型所用谱线参数相对精确,对其频率的时间响应标定相对精确。

图 1. 模型建立流程图

Fig. 1. Flow chart of model building

下载图片 查看所有图片

首先测量激光器的调制参数,得到激光器的频率-时间响应,根据数据库并结合实验室标定结果得到杂合谱线参数,给定初始流场参数,得到谱线展宽、吸光度,进而得到透过率函数,根据朗伯-比尔定律,乘以测量的无吸收光强信号,得到吸收后的强度信号;乘以参考信号经滤波得到n阶谐波信号Cnf,进一步得到被1f归一化的2f信号C2f/1f。波长调制技术对流场参数反演的本质是:将吸收模型得到的谐波信号与吸收后的实测强度信号解调得到的谐波信号进行比较,即模型中流场初始参数随实际测量环境而变化,并保证模拟过程与实测过程尽可能一致。选择表1所示的谱线参数,得到如图2所示的吸光度和扣除背景的归一化谐波信号C2f/1f_bgsub(压强p=1 atm,水汽体积分数ΧH2O=1%,温度T=1500 K,光程L=10 cm,调制频率fMod=200 kHz)。谱线参数包括谱线的中心频率υ0、谱线强度S(T0)、空气展宽系数γair、空气展宽系数对温度的依赖系数nair、 自展宽系数γself、压力频移系数δ0

表 1. 所选谱线参数

Table 1. Selected spectral line parameters

υ0 /cm-1S(T0) /(cm-2·atm-1)γair /(cm-1·atm-1)γself /(cm-1·atm-1)nairδ0 /(cm-1·atm-1)
7185.61.9×10-20.0410.1980.65-0.013

查看所有表

图 2. 模拟结果。(a)吸光度;(b)谐波信号

Fig. 2. Simulated results. (a) Absorbance; (b) harmonic signal

下载图片 查看所有图片

3.1 激光器调制参数

激光器频率被调制的同时,其强度也被调制。开始研究调制光谱技术时,认为二者的相位差为常数,之后被证实其值会随调制频率、注入电流、工作温度发生变化。为了更好地描述频率-时间响应,将表达式写为

υ(t)=υ-(t)+k=12(ak+bkt+ckt2)cos(kωt+θk),(9)

式中:akbkck分别为k次信号成分的多项式系数;θk为初始相位。

利用标准具对频率进行标定时,自由光谱范围(FSR)针对不同工作波长的激光器进行标定,其方法是同时扫描所选择的两条已知中心波数的吸收线,得到适用于特定波长的自由光谱范围。图3是中心波长为1392 nm的激光器同时采集的吸收信号和标定信号。

经数学运算求得吸光度,本次扫描吸收中心分别为7185.3943 cm-1和7185.5970 cm-1的两条水汽吸收线。图4是吸光度曲线,根据吸收中心绝对波数差值,推导出所用标准具FSR为0.6916 cm-1

图 3. 原始信号。(a)吸收信号;(b)标定信号

Fig. 3. Raw signals. (a) Absorption signal; (b) calibration signal

下载图片 查看所有图片

图 4. 吸光度及吸收中心

Fig. 4. Absorbance and absorption center

下载图片 查看所有图片

根据以上方法得到精确的FSR,对激光器的频率-时间响应进行标定,标定的依据是激光器出光强度和频率同时被调制,频率变化落后于强度变化,强度变大时频率降低,并且在标定时,确保寻峰准确,不遗漏。利用(8)式得到精确的激光器频率-时间响应如图5所示。

图 5. 激光器的频率-时间响应

Fig. 5. Frequency-time response of laser

下载图片 查看所有图片

该激光器的频率-时间响应数学表达式为

υt=0.357-4384t+3.7105t2+-0.094-158.6t+1.7105t2×

cos(2π3.2105t-0.045)+(9.610-3-17.8t)×cos(4π3.2105t-1.456)(10)

3.2 谱线参数选取

谱线参数包括中心频率υ0T0(参考温度296 K)及该温度下的线强S(T0)、空气展宽系数γair0、自展宽系数γself0、压力频移系数δ0以及温度依赖系数n。某温度下的线强、展宽系数均影响该温度下的谐波信号,常用的分子吸收光谱有HITRAN和HITEMP,根据HITRAN2012和HITEMP两个数据库谱线参数模拟吸光度(压强p=0.2 atm,ΧH2O=100%,温度T=1500 K,光程L=31.8 cm),如图6所示。

图 6. 两个数据库模拟得到的吸光度

Fig. 6. Absorption simulated on two databases

下载图片 查看所有图片

图6中可以看出:HITRAN2012数据库中缺少了低态能级较高的吸收线,这些吸收在常温下很弱,对于燃烧流场,HITEMP数据库中的吸收更贴近实际吸收;HITEMP数据库中的谱线参数大多基于半经验模型和理论计算得到,所以模型中的谱线参数应尽可能选择实验室标定结果。吸收线型选择Viogt线型,该线型描述了多普勒展宽和压力展宽两种展宽机制的作用。利用表1中的谱线参数,分别模拟不同压强下考虑临近吸收和不考虑临近吸收的吸光度,结果如图7所示。

图 7. 不同压强下考虑临近吸收和不考虑临近吸收时模拟的吸光度

Fig. 7. Simulated absorbance at different pressures with and without considering near absorption

下载图片 查看所有图片

图7中给出了ΧH2O=5%、光程L=10 cm、温度T=1500 K条件下,压强分别为1 atm和10 atm下模拟得到的吸光度,可见,临近吸收对所选吸收线存在影响。同样条件下加上高频正弦调制信号(调制频率为200 kHz),对模拟的吸收信号通过数字解调方法得到各自扣除背景的2f/1f信号,如图8所示。

图 8. 不同压强下考虑临近吸收和不考虑临近吸收时模拟的谐波信号。(a) 1 atm;(b) 10 atm

Fig. 8. Harmonic signals simulatedat different pressures with and without considering near absorption. (a) 1 atm; (b) 10 atm

下载图片 查看所有图片

图8可知,无论是常压还是高压,临近吸收同样影响所选吸收线的谐波信号,并且压强越大,影响越大。波长调制技术即利用谐波信号对流场参数进行反演。实际测量得到的谐波信号必定受所选吸收线以外的吸收的影响,将吸收模型与测量信号得到的谐波信号进行直接比较时,要求模型中的谱线参数应尽可能完备。

4 实验

为验证建立的吸收模型的实用性,利用近红外水汽吸收,结合波长调制技术进行高温流场温度的反演实验。实验选取吸收中心为7185.60 cm-1和6807.83 cm-1的两条吸收线,两条吸收线的归一化二次谐波信号的峰值比随温度而单调变化[12],利用该关系反演流场温度。实验装置如图9所示。

当设定温度T=1000 K、压强为1.34 atm时,解调得到的7185.60 cm-1和6807.83 cm-1两条吸收线的谐波信号及峰值如图10所示。

两只激光器均为日本NEL公司生产的分布反馈式半导体激光器,中心波长分别为1392 nm和1469 nm,其锯齿扫描频率和正弦调制频率均为1.1 kHz和322 kHz。恒温池为科晶公司生产的型号为1800X的管式高温炉,炉膛中间有长度为31.8 cm的恒温流场,高温炉温度由比例积分微分(PID)温度模块控制,可根据需要设定温度。炉膛内的压强由德国浦发公司生产的型号为CCR-361的真空计测量。

图 9. 实验装置简图

Fig. 9. Diagram of experimental setup

下载图片 查看所有图片

图 10. 扣除背景的2f/1f信号。(a)吸收中心为7185.60 cm-1;(b)吸收中心为6807.83 cm-1

Fig. 10. Signals of 2f/1f with background subtracted. (a) Central frequency is 7185.60 cm-1; (b) central frequency is 6807.83 cm-1

下载图片 查看所有图片

实验过程设定7个温度段,如图11所知,每隔100 K设定一个温度阶梯,从900 K升到1500 K,每个温度阶梯连续采集100组数据,对其结果分段计算测量得到的平均温度、绝对误差、相对误差以及标准差,如表2所示。

测量温度与设定温度的最大绝对误差为27.7 K,最大相对误差为3.1%,最大标准差为12.7 K,吻合得较好,验证了模型的准确性。误差来源主要有4个方面:1)模型中光谱参数带来的误差;2)激光器调制参数测量带来的误差;3)压强、光程直接测量带来的误差;4)恒温区流场不稳定带来的误差。各温度阶梯测量结果的最大偏差出现在900 K,其原因主要是所选的中心波长为6807.83 cm-1的吸收线低态能级为3319.4 cm-1,相对较大,在900 K时线强较弱,对应的谐波信号相对较弱,导致误差较大。随着温度上升,误差降低。通过对高温炉温度的精确测量,验证了波长调制吸收光谱模型的准确性。

图 11. 设定温度和测量温度

Fig. 11. Set and measured temperatures

下载图片 查看所有图片

表 2. 温度分析

Table 2. Analysis of temperature

Set temperature /KAverage temperature /KAbsolute error /KRelative error /%Standard deviation /K
900927.727.23.112.7
10001002.02.00.212.0
11001112.612.61.112.7
12001222.522.51.95.7
13001323.323.31.84.6
14001379.1-20.9-1.54.6
15001478.2-21.8-1.59.8

查看所有表

5 结论

研究了基于波长调制技术的吸收模型精确建立方法,吸收模型建立的准确性直接影响流场参数反演的准确性。以H2O为目标分子,重点讨论激光器频率-时间响应精确标定,提出了一种简单实用的FSR标定方法。另外,在燃烧流场应用领域,谱线参数相对复杂,谱线参数参考具有更多“热线”吸收的HITEMP数据库并对其标定。利用建立的吸收模型测量了稳定高温流场温度,测量的最大相对误差为3.1%,验证了模型的准确性,有利于波长调制技术在燃烧流场诊断方面的应用。

参考文献

[1] Sun K, Chao X, Sur R, et al. Analysis of calibration-free wavelength-scanned wavelength modulation spectroscopy for practical gas sensing using tunable diode lasers[J]. Measurement Science and Technology, 2013, 24(12): 125203.

[2] 阚瑞峰, 夏晖晖, 许振宇, 等. 激光吸收光谱流场诊断技术应用研究与进展[J]. 中国激光, 2018, 45(9): 0911005.

    Kan R F, Xia H H, Xu Z Y, et al. Research and progress of flow field diagnosis based on laser absorption spectroscopy[J]. Chinese Journal of Lasers, 2018, 45(9): 0911005.

[3] Goldenstein CS. Wavelength-modulation spectroscopy for determination of gas properties in hostile environments[D]. Stanford: Stanford University, 2014.

[4] Sur R, Sun K, Jeffries J B, et al. Scanned-wavelength-modulation-spectroscopy sensor for CO, CO2, CH4 and H2O in a high-pressure engineering-scale transport-reactor coal gasifier[J]. Fuel, 2015, 150: 102-111.

[5] Goldenstein C S, Spearrin R M, Jeffries J B, et al. Infrared laser-absorption sensing for combustion gases[J]. Progress in Energy and Combustion Science, 2017, 60: 132-176.

[6] Goldenstein C S, Schultz I A, Spearrin R M, et al. Scanned-wavelength-modulation spectroscopy near 2.5 μm for H2O and temperature in a hydrocarbon-fueled scramjet combustor[J]. Applied Physics B, 2014, 116(3): 717-727.

[7] Sun K, Sur R, Chao X, et al. TDL absorption sensors for gas temperature and concentrations in a high-pressure entrained-flow coal gasifier[J]. Proceedings of the Combustion Institute, 2013, 34(2): 3593-3601.

[8] Wilson G V H. Signal amplitudes of modulation-broadened Lorentz and Gaussian magnetic resonances[J]. Journal of Applied Physics, 1965, 36(11): 3505-3506.

[9] Philippe L C, Hanson R K. Laser diode wavelength-modulation spectroscopy for simultaneous measurement of temperature, pressure, and velocity in shock-heated oxygen flows[J]. Applied Optics, 1993, 32(30): 6090-6103.

[10] SunK. Utilization of multiple harmonics of wavelength modulation absorption spectroscopy for practical gas sensing[D]. Stanford: Stanford University, 2013.

[11] Goldenstein C S, Schultz I A, Jeffries J B, et al. Two-color absorption spectroscopy strategy for measuring the column density and path average temperature of the absorbing species in nonuniform gases[J]. Applied Optics, 2013, 52(33): 7950-7962.

[12] 孙鹏帅, 张志荣, 夏滑, 等. 基于波长调制技术的温度实时测量方法研究[J]. 光学学报, 2015, 35(2): 0230001.

    Sun P S, Zhang Z R, Xia H, et al. Study on real-time temperature measurement based on wavelength modulation technology[J]. Acta Optica Sinica, 2015, 35(2): 0230001.

张步强, 许振宇, 刘建国, 夏晖晖, 聂伟, 袁峰, 阚瑞峰. 燃烧流场波长调制光谱吸收模型的研究[J]. 中国激光, 2019, 46(7): 0711001. Buqiang Zhang, Zhenyu Xu, Jianguo Liu, Huihui Xia, Wei Nie, Feng Yuan, Ruifeng Kan. Absorption Model of Wavelength Modulation Spectroscopy in Combustion Flow Field[J]. Chinese Journal of Lasers, 2019, 46(7): 0711001.

本文已被 5 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!