非分光多波长调制吸收光谱测量方法 下载: 656次
1 引言
脉冲爆轰发动机(PDE)是一种新型发动机,对其工作过程中产生的燃气流场组分温度浓度等指标进行分析,能够更好地研究其工作机理,为评估工作性能提供数据支持[1-5]。在PDE高温高压的工作环境下,常规的接触式测量技术难以适用,可调谐半导体激光吸收光谱(TDLAS)技术作为非接触气体测试技术,具有灵敏度高、分辨率高、响应速度快和环境适应性强等优点,适用于爆轰流场的在线检测[6-8]。
目前已经有许多学者将TDLAS技术应用于不同环境下的气体测量,主要分为直接吸收法与波长调制法两种。阚瑞峰等[9-10]选取了1391.7 nm中心处单根谱线采取直接吸收法,搭建双光路测量系统实现了风洞中H2O速度测量,并且还应用直接吸收法对平焰炉中的H2O温度进行测量;张雅琪等[11]应用直接吸收法,选取7153.7208 cm-1/7153.7488 cm-1和7154.353 cm-1/7154.3538 cm-1两组谱线对火焰中H2O浓度进行测量;吕晓静等[12-14]应用直接吸收法,通过测量CO、H2O等气体组分的浓度、温度以及速度参数,对气液两相脉冲爆轰发动机的工作性能进行研究;Rieker等[15-16]在高温条件下应用波长调制技术选用4条谱线对H2O进行温度测量;Li等[17]通过一次谐波标定二次谐波方法,实现了冲击管中H2O的浓度测量。
波分复用技术可在同一根光纤中同时传输多波长光信号,相比时分复用技术能保证更高的测量效率[18-20]。但波分复用技术后续需要通过光栅等手段将不同波长的信号分离,且波长相近时可能会导致相互之间产生干扰[21-23],测量系统相对复杂,对测量环境要求较高。
本文提出一种基于未分光多谱线解调的波长调制吸收光谱爆轰燃气测量方法,通过解调合成透射信号,无需通过分光装置来实现吸收谱线的分离,简化了测试系统,适用于复杂现场测量,具有重要研究意义和应用前景。
2 波长调制分离解调方法
激光吸收光谱技术以分子吸收光谱理论为基础,一束频率连续变化的平行光穿透气体介质后,气体会对特定频率的光产生不同程度的吸收[24-26],激光强度会相应减弱,该过程遵循比尔朗伯定律[15]:
式中:P(l)、X(l)和T(l)分别为总静压、摩尔分数和温度;L为光程长度;
式中:参考温度T0为296 K,该温度下谱线强度可通过HITRAN光谱数据库得到;Q(T)为分子内部分割函数,可通过多项式拟合求近似值;h是普朗克常数;k为玻尔兹曼常数;c为光的传播速度;
分布式反馈(DFB)激光器在高频电流驱动下激光强度可以描述为[6]
式中,a2、a1和a0分别是激光强度扫描波长非线性系数、线性系数和常数。波长调制技术是通过在锯齿扫描信号上叠加高频正弦信号驱动激光器实现的。激光器出射频率为
式中:v为中心频率;a为调制深度;f为调制频率。出射激光强度表示为
定义透射光强与入射光强的比值
傅里叶级数展开为
式中:
根据比尔朗伯定律计算得到透射光强乘以参考信号获得二次谐波的X分量和Y分量为
式中,
式中:G为探测器增益系数;H1、H3为奇函数,在中心处为0。H2值为
二次谐波峰值与气体浓度成正比关系。定义调制系数m,其大小为两倍的调制幅度除以线宽
式中,M=1-x2+m2,在中心处x=0,此时有二次谐波峰值P2f:
通过查询HITRAN数据库中两条谱线的加宽系数可知,在一定浓度范围内,线宽近似相等,m也近似相等,可用二次谐波峰值比代替谱线强度比:
HITRAN数据库中谱线强度随温度变化曲线如
图 1. 谱线强度/谱线强度比与温度的对应关系(7185.60 cm-1/7444.35 cm-1)。(a)谱线强度;(b)谱线强度比
Fig. 1. Correspondence between line strength/line strength ratio and temperature (7185.60 cm-1/7444.35 cm-1). (a) Line strength; (b) line strength ratio
因单一激光器波长扫描范围较小,且气体吸收谱线并不接近,所以需要通过波分复用技术将不同波段的激光耦合进一根光纤内,经过气体吸收区后再将不同波段分开,目前主要通过光栅等手段对合波后的激光进行分光处理。本文不再通过借助后续光路设计来实现分光,对经过吸收的透射信号解调得到二次谐波信号,调节吸收谱线中心位波长位置,尽可能减少不同谱线之间叠加带来的影响,实现二次谐波信号的分离。
3 仿真分析
3.1 非分光多谱线波长调制吸收光谱特征
首先对非分光多谱线波长调制吸收光谱特征进行分析。选取7185.60 cm-1/7444.35 cm-1两条H2O特征吸收谱线进行扫描,H2O温度设定为1273 K,压力为2 atm,摩尔分数为0.1,测量路径长度为10 cm。分别对两个波段激光器加载锯齿波信号,对两条H2O特征吸收谱线进行扫描,波长扫描范围分别为7443.5~7447.5 cm-1与7182.5~7186.5 cm-1,谱线间距离为2.25 cm-1,锯齿波扫描频率为1 Hz。分别加载正弦调制信号,正弦波调制频率为1 kHz,调制系数为2。通过仿真获取的一个完整扫描周期内的调制吸收光谱信号如
对采集得到的调制激光吸收光谱信号进行解调制,解调制结果如
图 3. 7185.60 cm-1/7444.35 cm-1合波透射信号解调后二次谐波信号值
Fig. 3. Second harmonic signal after demodulation of multiplexed transmission signals of 7185.60 cm-1/7444.35 cm-1
3.2 调制系数与吸收谱线相对位置对于测量的影响
为了研究调制系数m对于温度测量的影响。在仿真中保持其他测量工况参数(T=1273 K,P=2 atm,X=0.1,L=10 cm)不变,通过改变调制幅度来改变调制系数,进而考察调制系统对于温度测量精度的影响。仿真结果如
图 4. 不同调制系数m对计算结果的影响。(a)二次谐波信号;(b)温度计算相对误差
Fig. 4. Influence of modulation coefficient m on simulation results. (a) Second harmonic signal; (b) temperature relative error
为了探究两条谱线相对位置对测量信号的影响,不改变环境工况条件,只改变激光器初始波长位置,从而调整两条谱线相对位置,分析谱线叠加的影响。仿真结果如
图 5. 谱线相对位置对仿真结果的影响
Fig. 5. Influence of relative position of spectral line on simulation result
3.3 不同调制频率对于非分光吸收光谱测量的影响
采用非分光波长调制吸收光谱技术进行测量时,两条谱线之间的相对位置对于测量结果具有重要影响。当测量环境压力较高或谱线对位置较为接近时,谱线间会产生明显干扰。为了消除这种干扰,可以对两条谱线加载不同的调制频率,从而使得测量在不同的调制频段内同步进行。
对7185.60 cm-1/7444.35 cm-1两条H2O特征吸收谱线在相对位置距离为0.2 cm-1时分别采用不同调制频率(f1为800 Hz,f2为1 kHz)进行扫描,测量得到非分光时域信号及其频谱如
图 6. 不同调制频率信号叠加结果。(a)时域图;(b)频域图
Fig. 6. Different modulation frequency signal superposition results. (a) Time domain diagram; (b) frequency domain diagram
对非分光吸收光谱时域信号的频谱进行细致分析。首先,在调制信号基频及其倍频位置处出现明显的吸收信号,这与调制激光吸收光谱理论相吻合。调制吸收光谱的解调制过程就是将特定倍频信号中的吸收光谱信息进行提取;其次,在调制信号基频及其倍频位置处信号脉宽较窄(小于100 Hz),属于窄脉冲信号。因此,当两条谱线调制频率设定为不同值时,其调制信号基频及其倍频均不重叠,可以将其中两条谱线不同倍频处的谐波信号准确还原。
针对
图 7. 两条谱线相同调制频率与不同调制频率的解调结果。(a)调制频率均为1 kHz;(b)7185.60 cm-1谱线调制频率为800 Hz;(c)7444.35 cm-1谱线调制频率为1 kHz
Fig. 7. Demodulation results of two spectral lines with same modulation frequency and different modulation frequencies. (a) Modulation frequencies both 1 kHz; (b) line 7185.60 cm-1 with modulation frequency of 800 Hz; (c) line 7444.35 cm-1 with modulation frequency of 1 kHz
值得注意的是,选取的两条吸收谱线调制频率在频谱图中基频与倍频信号间相互重叠时,会导致解调时产生更高阶谐波信号,信号间出现干扰;选取的两条谱线调制频率过近时,会导致解调制时噪声较大,信噪比较低,谐波信号失真,温度测量误差较大。因此两者调制频率不宜选取相近及倍数关系的组合。
4 针对气液两相爆轰流场测量的仿真分析与结果讨论
脉冲爆轰发动机工作时爆轰燃气温度与压力变化范围大,可以用于验证本文波长调制气体测量方法。本文以气液两相脉冲爆轰发动机管口处爆轰燃气为研究对象,采用仿真方法分析不同环境工况条件下的测量效果。
仿真中脉冲爆轰发动机长1.5 m,内径80 mm,分别以汽油和氧气为燃料和氧化剂,工作时初始温度设定为296 K,初始压力设定为1 atm,燃料填充系数设定为1。建立考虑粘性影响下的气液两相爆轰流场轴对称控制方程,并采用时空守恒元/求解元(CE/SE)方法进行求解。计算中设定燃料氧化剂化学反应释放能量被气体吸收,忽略发动机壁面与外界的能量交换。通过仿真计算得到脉冲爆轰发动机爆轰过程管内流场变化全过程,典型时刻的爆轰流场速度场及压力场如
图 8. CE/SE方法仿真得到的管口速度和压力图。(a)速度图;(b)压力图
Fig. 8. Speed and pressure graphs obtained by CE/SE method simulation. (a) Speed graph; (b) pressure graph
针对脉冲爆轰发动机管口处爆轰燃气发展变化过程进行仿真测量。测试系统方案如
表 1. 不同时刻通过仿真得到的温度和计算得到的温度比较
Table 1. Comparison of temperatures obtained through simulation and calculated temperatures at different time
|
爆轰波到达测量位置后,管内压力由爆轰波阵面压力值迅速下降,2.2 ms时压力达到最大为32 atm;在3.5 ms后接近常压。与之相比,爆轰燃气温度变化较为缓慢,经过2.8 ms后温度达到最大值2199 K;在爆轰波过后5 ms后降低至1100 K以下,随后温度下降趋势减缓;在20 ms后到达600 K左右并趋于稳定。研究中选取爆轰过程的几个典型时刻,对爆轰燃气温度进行了模拟测量。
图 10. 典型时刻二次谐波解调结果。(a)调制频率均为500 kHz解调结果;(b)2.2 ms时刻7444.35 cm-1谱线调制频率500 kHz解调结果;(c)2.2 ms时刻7185.60 cm-1谱线调制频率400 kHz解调结果
Fig. 10. Second harmonic demodulation results at typical moments. (a) Modulation frequencies both 500 kHz demodulation result; (b) line 7444.35 cm-1 with modulation frequency 500 kHz demodulation result at 2.2 ms; (c) line 7185.60 cm-1 with modulation frequency 400 kHz demodulation result at 2.2 ms
根据本文方法计算得到温度值与实时温度对比如
5 结论
主要研究了可适用于复杂环境测量的非分光多谱线波长调制气体测量方法。在传统波长调制测量方法基础上,通过控制激光器温度和输出电流,使得叠加后透射激光信号中吸收光谱谱线间存在合适的相对距离,进而无需分光即可得到多条调制谱线的二次谐波信号。
利用7185.60 cm-1/7444.35 cm-1谱线对,分别对不同调制系数和谱线间相对距离等工况下的二次谐波信号进行了分析。结果表明,调制系数m在2.2附近时最为合适,调制系数过大或者过小都将引起较大的测量误差。同时,两条谱线相对位置大于两倍线宽时,非分光解调信号对于测量无影响。为了进一步解决在特殊工况条件下谱线间相互叠加干扰对于测量的影响,分别对两条测量谱线加载不同调制频率。理论分析了不同调制频率下调制吸收光谱信号频谱,研究发现调制信号基频及其倍频位置处出现明显的吸收信号,信号脉宽小于100 Hz;当两条谱线调制频率不同时,两条谱线调制信号的基频及其倍频均不重叠,因此可以将两条谱线二次谐波信号完整提取。
为验证本文波长调制气体测量方法,以气液两相爆轰发动机爆轰过程为研究对象,通过采用CE/SE方法获得爆轰过程中燃气参数变化过程,并以此建立模型进行求解。结果表明,非分光多谱线波长调制吸收光谱测量方法可以准确测量爆轰过程燃气温度变化特征,爆轰过程燃气温度测量误差小于9.63%。该方法适用于基于波长调制的多谱线在线测量领域,无需搭建分光装置即可获取测量结果,优化了测量系统,更加适用于在现场进行实时测量,有广阔的工程应用前景。
[1] 周帏, 翁春生. PDE出口爆轰波与射流诱发水中流场变化规律的数值仿真[J]. 水下无人系统学报, 2017, 25(3): 167-173.
Zhou W, Weng C S. Numerical simulation of the underwater flow field induced by detonation wave and jet from pulse detonation engine outlet[J]. Journal of Unmanned Undersea Systems, 2017, 25(3): 167-173.
[2] 黄孝龙, 李宁, 翁春生. 多管PDE爆轰噪声传播过程物理特性实验研究[J]. 声学学报, 2021, 46(3): 415-422.
Huang X L, Li N, Weng C S. Experimental study on the physical characteristics of noise propagation process of multi-tube pulse detonation engine[J]. Acta Acustica, 2021, 46(3): 415-422.
[3] 蒋弢, 翁春生. 进气温度和当量比对脉冲爆轰发动机工作过程影响[J]. 航空动力学报, 2015, 30(1): 29-37.
Jiang T, Weng C S. Effect of inlet temperature and equivalence ratio on working process of pulse detonation engine[J]. Journal of Aerospace Power, 2015, 30(1): 29-37.
[4] 翁春生, 王杰, 白桥栋, 等. 脉冲爆轰发动机进气压力对爆轰影响的实验研究[J]. 弹道学报, 2008, 20(3): 1-4.
Weng C S, Wang J, Bai Q D, et al. Experimental study on effect of inlet pressure on detonation in pulse detonation engine[J]. Journal of Ballistics, 2008, 20(3): 1-4.
[5] 许桂阳, 翁春生, 李宁. 脉冲爆轰发动机近场爆轰噪声特性研究[J]. 南京理工大学学报, 2014, 38(6): 781-786.
Xu G Y, Weng C S, Li N. Detonation noise characteristics of pulse detonation engine in near-field[J]. Journal of Nanjing University of Science and Technology, 2014, 38(6): 781-786.
[6] 吕晓静, 李宁, 徐恩华, 等. 基于TDLAS技术的脉冲爆轰发动机管内燃气CO的测试研究[J]. 海峡科技与产业, 2018(6): 77-79.
Lü X J, Li N, Xu E H, et al. Research on carbon monoxide test of gas in pulse detonation engine tube based on TDLAS technology[J]. Technology and Industry Across the Straits, 2018(6): 77-79.
[7] 吕晓静, 李宁, 徐恩华, 等. 基于TDLAS技术结合多普勒测速的脉冲爆轰发动机冲量测试[J]. 科技风, 2018(32): 154-155.
Lü X J, Li N, Xu E H, et al. Impulse test of pulse detonation engine based on TDLAS technology and Doppler velocity measurement[J]. Technological Wind, 2018(32): 154-155.
[8] 吕晓静, 李宁, 翁春生. 基于激光吸收光谱技术的PDE管内燃气测试[J]. 激光杂志, 2018, 39(10): 137-142.
Lü X J, Li N, Weng C S. Exhaust diagnostics in PDE tube based on absorption spectrum technology[J]. Laser Journal, 2018, 39(10): 137-142.
[9] 阚瑞峰, 夏晖晖, 许振宇, 等. 激光吸收光谱流场诊断技术应用研究与进展[J]. 中国激光, 2018, 45(9): 0911005.
[10] 彭于权, 阚瑞峰, 许振宇, 等. 基于可调谐半导体激光吸收光谱技术的甲烷/空气预混平焰炉温度测量[J]. 大气与环境光学学报, 2019, 14(3): 228-234.
[11] 张雅琪. 可调谐半导体激光吸收光谱技术用于火焰中CO和H2O的浓度测量研究[D]. 杭州: 浙江大学, 2019.
ZhangY Q. Concentration measurement of CO and H2O in flame using tunable diode laser absorption spectroscopy[D]. Hangzhou: Zhejiang University, 2019.
[13] 吕晓静, 李宁, 翁春生. 基于多谱线吸收光谱技术的脉冲爆轰发动机管外流场测试研究[J]. 光谱学与光谱分析, 2016, 36(3): 624-630.
[14] 吕晓静, 李宁, 翁春生. 基于双光路吸收光谱技术的气液两相爆轰燃气诊断技术研究[J]. 光谱学与光谱分析, 2014, 34(3): 582-586.
[15] Rieker G B, Jeffries J B, Hanson R K. Calibration-free wavelength-modulation spectroscopy for measurements of gas temperature and concentration in harsh environments[J]. Applied Optics, 2009, 48(29): 5546-5560.
[16] Rieker G B, Liu X, Li H, et al. Measurements of near-IR water vapor absorption at high pressure and temperature[J]. Applied Physics B, 2007, 87(1): 169-178.
[17] Li H J, Rieker G B, Liu X, et al. Extension of wavelength-modulation spectroscopy to large modulation depth for diode laser absorption measurements in high-pressure gases[J]. Applied Optics, 2006, 45(5): 1052-1061.
[18] Liu X, Jeffries J B, Hanson R K, et al. Development of a tunable diode laser sensor for measurements of gas turbine exhaust temperature[J]. Applied Physics B, 2006, 82(3): 469-478.
[19] Liu J T C, Jeffries J B, Hanson R K. Large-modulation-depth 2f spectroscopy with diode lasers for rapid temperature and species measurements in gases with blended and broadened spectra[J]. Applied Optics, 2004, 43(35): 6500-6509.
[20] Liu J T C, Jeffries J B, Hanson R K. Wavelength modulation absorption spectroscopy with 2f detection using multiplexed diode lasers for rapid temperature measurements in gaseous flows[J]. Applied Physics B, 2004, 78(3/4): 503-511.
[21] Mattison D W, Oehlschlaeger M A, Morris C I, et al. Evaluation of pulse detonation engine modeling using laser-based temperature and OH concentration measurements[J]. Proceedings of the Combustion Institute, 2005, 30(2): 2799-2807.
[22] Mattison D W, Jeffries J B, Hanson R K, et al. In-cylinder gas temperature and water concentration measurements in HCCI engines using a multiplexed-wavelength diode-laser system: sensor development and initial demonstration[J]. Proceedings of the Combustion Institute, 2007, 31(1): 791-798.
[23] Sanders S T, Mattison D W, Jeffries J B, et al. Rapid temperature tuning of a 1.4-mum diode laser with application to high-pressure H2O absorption spectroscopy[J]. Optics Letters, 2001, 26(20): 1568-1570.
[24] 侯月, 黄克谨, 于冠一, 等. 基于红外TDLAS技术的高精度CO2同位素检测系统的研制[J]. 红外与激光工程, 2021(4): 20200083.
Hou Y, Huang K J, Yu G Y, et al. Development on high precision CO2 isotope measurement system based on infrared TDLAS technology[J]. Infrared and Laser Engineering, 2021(4): 20200083.
[25] 王国水, 郭奥, 刘晓楠, 等. TDLAS气体检测系统仿真与影响因素分析[J]. 光谱学与光谱分析, 2021, 41(10): 3262-3268.
[26] 吕文静, 李红莲, 李文铎, 等. TDLAS技术调制参量的优化及实验研究[J]. 激光技术, 2021(3): 336-343.
Lü W J, Li H L, Li W D, et al. Optimization and experimental research on modulation parameters of TDLAS technology[J]. Laser Technology, 2021(3): 336-343.
[27] Reid J, Labrie D. Second-harmonic detection with tunable diode lasers: comparison of experiment and theory[J]. Applied Physics B, 1981, 26(3): 203-210.
[28] Kang Y, Li N, Weng C S, et al. Acoustic characteristics of pulse detonation engine with ellipsoidal reflector[J]. Chinese Physics B, 2018, 27(10): 104703.
Article Outline
石逸凡, 李宁, 黄孝龙, 康杨, 李灿, 翁春生. 非分光多波长调制吸收光谱测量方法[J]. 激光与光电子学进展, 2023, 60(1): 0107001. Yifan Shi, Ning Li, Xiaolong Huang, Yang Kang, Can Li, Chunsheng Weng. Non-Dispersive Multi-Wavelength Modulation Absorption Spectrum Measurement Method[J]. Laser & Optoelectronics Progress, 2023, 60(1): 0107001.