CN107036983B - 基于调制光谱实现非均匀燃烧场气体参数定量测量的方法 - Google Patents

基于调制光谱实现非均匀燃烧场气体参数定量测量的方法 Download PDF

Info

Publication number
CN107036983B
CN107036983B CN201710142698.4A CN201710142698A CN107036983B CN 107036983 B CN107036983 B CN 107036983B CN 201710142698 A CN201710142698 A CN 201710142698A CN 107036983 B CN107036983 B CN 107036983B
Authority
CN
China
Prior art keywords
line
gas parameter
line width
temperature
wavelength modulation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710142698.4A
Other languages
English (en)
Other versions
CN107036983A (zh
Inventor
洪延姬
王广宇
金星
屈东胜
宋俊玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Original Assignee
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Peoples Liberation Army Strategic Support Force Aerospace Engineering University filed Critical Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority to CN201710142698.4A priority Critical patent/CN107036983B/zh
Publication of CN107036983A publication Critical patent/CN107036983A/zh
Application granted granted Critical
Publication of CN107036983B publication Critical patent/CN107036983B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/171—Systems in which incident light is modified in accordance with the properties of the material investigated with calorimetric detection, e.g. with thermal lens detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/1717—Systems in which incident light is modified in accordance with the properties of the material investigated with a modulation of one or more physical properties of the sample during the optical investigation, e.g. electro-reflectance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/1717—Systems in which incident light is modified in accordance with the properties of the material investigated with a modulation of one or more physical properties of the sample during the optical investigation, e.g. electro-reflectance
    • G01N2021/1725—Modulation of properties by light, e.g. photoreflectance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N2021/3125—Measuring the absorption by excited molecules
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/06—Illumination; Optics
    • G01N2201/061—Sources
    • G01N2201/06113—Coherent sources; lasers

Abstract

本发明基于调制光谱技术，提供了一种实现非均匀燃烧场气体参数定量测量的方法，该方法包括波长调制光谱测量步骤、波长调制光谱拟合步骤、非均匀流场中气体参数计算步骤。其具体内容是选择在被测温度范围内线强度随温度线性变化的吸收谱线；通过实验测得透射光强信号，获得实验测得的谐波信号；以积分吸光度、高斯线宽和洛伦兹线宽为自由变量，基于最小二乘方法拟合谱线的谐波信号；计算非均匀燃烧场中的气体参数。

Description

基于调制光谱实现非均匀燃烧场气体参数定量测量的方法
技术领域
本发明属于光谱测量技术领域，是针对非均匀燃烧场气体参数定量测量开发的一种方法。
背景技术
针对燃烧场开始研究时，大部分的燃烧场都是非均匀的，非均匀流场中气体参数的定量测量对于分析燃烧效率以及研究燃烧过程具有十分重要的意义。热电偶、压力传感器等常用的测量手段具有明显的局限性，如侵入式测量干扰流场、容易产生激波、测量精度有限等。本发明涉及的波长调制光谱(wavelength modulation spectroscopy,WMS)技术属于光谱测量技术，能够在不干扰流场的情况下实现流场内气体参数的快速测量，在燃烧场诊断等领域具有广阔的应用前景。
WMS技术是利用吸收光谱信号实现气体参数的测量，在瞬态、强震动等恶劣环境中优势十分明显。利用WMS实现均匀燃烧场内气体参数的定量测量问题已经得到解决，但由于吸收光谱的光线积分效应，如何实现非均匀燃烧场内气体参数的定量测量成为了一个难题。
下述的文献也涉及到了燃烧场气体参数测量实现方法的相关内容。
1、美国斯坦福大学L.S.Chang,C.L.Strand等人(“Supersonic Mass-FluxMeasurements via Tunable Diode Laser Absorption and Nonuniform FlowModeling”，AIAA Journal,Vol.49,No.12,2011,pp.2783-2791.),发现了采用波长调制光谱测量时，边界层对测量结果影响很大，选择对边界层不敏感的谱线对，能够有效减少边界层的影响。
2、美国斯坦福大学Liu X，J.B.Jeffries等人(“Measurement of Non-uniformTemperature Distributions Using Line-of-Sight Absorption Spectroscopy”，AIAAJournal,Vol.45,No.2,2007,pp.411-419.),提出了采用多波长方法实现非均匀燃烧场气体参数的测量，研究了吸收谱线数目和温度区间长度对测量结果的影响。然而多波长方法的缺点十分明显，多波长方法主要有形状拟合法和温度分区法，形状拟合法需要根据经验预先给出流场内温度和组分浓度的分布形状，在实际中的应用十分有限而温度分区法只能获得温度和组分浓度沿光路的概率密度函数，无法获得在光路上的具体分布。
3、英国剑桥大学的Cai W和C.F.Kaminski等人(“Multiplexed absorptiontomography with calibration-free wavelength modulation spectroscopy”，AppliedPhysics Letters,No.104,2014,pp.154106.),提出了基于波长调制光谱方法实现非均匀流场区域的二维重建研究。该方法将吸收光谱技术与图像重建技术相结合，能够实现气体的温度和组分浓度的二维分布测量。这种方法需要大量的探测和接收装置，才能保证重建结果质量，对测量系统的硬件设备要求较高。另外该方法以基于搜索的重建算法为核心，通过大量的数值仿真过程，才能获得较好的重建结果，具有空间分辨率低、计算周期长等显著缺点。
4、清华大学车璐等人(“Calibration-Free Wavelength Modulation For GasSensing In Tunable Diode Laser Absorption Spectroscopy”，Applied Physical B,No.117,2014,pp.1211–1219.)提出了基于免标定波长调制光谱方法的气体参数测量。分析了流场中气体压强和吸收分子组分浓度的定量测量过程。然而该方法只适用于均匀流场中，无法定量测量非均匀流场中气体参数。
本发明实现了非均匀燃烧场中气体参数的定量测量，突破了理论和技术上的难题，形成了一套切实可行的方法。
发明内容
本发明的目的是解决非均匀燃烧场中气体参数的定量测量问题，提出了一种基于调制光谱实现非均匀燃烧场气体参数定量测量的方法。从理论上分析了实现非均匀流场气体参数定量测量存在的问题及解决的方法并通过数值仿真验证了该方法的可能性。基于该方法能够实现非均匀燃烧场中气体参数的非侵入式定量测量，对于研究燃烧过程、提高燃烧效率以及减少污染物排放等意义重大。尤其在发动机流场诊断领域，利用该技术能够实现发动机内部流场的温度等参数的定量测量，为发动机的故障诊断提供可靠的依据。
本发明提供了一种基于调制光谱实现非均匀燃烧场气体参数定量测量的方法，包括：
步骤一、波长调制光谱测量步骤；
步骤二、波长调制光谱拟合步骤；
步骤三、非均匀流场中气体参数计算步骤；
其特征在于：
步骤一的波长调制光谱测量步骤包括：(1)预估非均匀流场中气体温度的变化范围，选择具有合适低跃迁态能级的两条谱线，使得所选谱线的线强度在该温度范围内具有随温度近似线性变化的特性；(2)通过实验测量所选谱线穿过非均匀流场后的透射光强信号随时间的变化；(3)创建数值锁相放大器，解调实验测得的透射光强信号，得到两条谱线的谐波信号；
步骤二的波长调制光谱拟合步骤包括：(1)计算非均匀流场中吸收谱线的积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC的变化范围；(2)以积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC为自变量，仿真所选谱线的透射光强信号随时间的变化；(3)创建数值锁相放大器，解调仿真的透射光强信号，得到仿真的谐波信号；(4)基于最小二乘方法拟合实验测得的谐波信号，获得两条谱线穿过非均匀流场后的积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC；
步骤三的非均匀流场中气体参数计算步骤包括：(1)利用线性函数拟合气体温度变化范围内所选谱线的线强度随温度的变化；(2)利用两条谱线的积分吸光度A计算非均匀流场内的气体温度和吸收分子的分压。
本发明基于调制光谱方法，在简化数据处理的同时实现了非均匀燃烧场气体参数的定量测量，优点如下：
(1)解决了非均匀燃烧场气体参数定量测量问题；选择在温度变化范围内具有近似线性变化特性的线强度的吸收谱线开展测量，简化了谱线的线强度随温度的变化趋势，结合实验和仿真的谐波信号，通过光谱拟合法实现了气体参数的定量测量。
(2)开发了基于最小二乘方法实现谐波信号拟合的方法；以积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC为自变量，实现了基于最小二乘方法实现谐波信号拟合的方法。
附图说明
图1为本发明的基于调制光谱实现非均匀燃烧场气体参数定量测量的流程图；
图2为本发明的实验测得的透射光强信号随时间的变化；
图3为本发明的基于最小二乘方法的拟合谐波信号流程图；
图4为本发明的仿真和实验测得的谐波信号对比图；
图5为本发明的所选谱线的线强度随温度线性变化的示意图；
具体实施方案
结合附图和实施案例对本发明的基于调制光谱实现非均匀燃烧场气体参数定量测量的实现方法作进一步详细描述。图1给出了本发明的基于调制光谱实现非均匀燃烧场气体参数定量测量的流程图。
本发明的实现方案如下：本发明方法分三个步骤，即波长调制光谱测量步骤、波长调制光谱拟合步骤、非均匀流场中气体参数计算步骤；
本发明基于调制光谱实现非均匀燃烧场气体参数定量测量的实现方法，实施步骤如下：
步骤一、波长调制光谱测量步骤
基于调制光谱实现非均匀燃烧场气体参数定量测量的实施方法第一步，是波长调制光谱测量步骤，共分为3个任务，实施方案如下：
①、选择合适的吸收谱线
预估被测流场的气体温度变化范围，根据谱线的低跃迁态能级选择两条在该温度变化范围内具有线性变化特性的谱线。
②、采集波长调制光谱测量信号
调节激光器的电流和温度，使得激光波长在所选两条谱线附近。激光器输出的光经合束器后变成一束光。激光穿过非均匀流场后由探测器接收，而后读取探测到的光电信号。该信号是经过非均匀燃烧场的透射光强信号，在吸收谱线中心频率附近产生了吸收，可作为了实验测得的调制光谱信号。实验测得的透射光强信号随时间的变化如图2所示。
③、解调调制光谱信号
利用锁相放大技术解调实验测得的透射光强信号，获得所选谱线谐波信号。两条谱线的调制频率不同，锁相放大时通过调制频率区别两条谱线的谐波信号。
步骤二、波长调制光谱拟合步骤
①、计算谱线积分吸光度的变化范围
根据预估的燃烧场气体温度、压强和吸收分子组分浓度的变化范围，计算所选谱线的积分吸光度A的变化范围，计算公式为
A＝PXS(T)L (1)
式中，P为气体压强，X为吸收分子的组分浓度，S(T)为吸收谱线的线强度，L为吸收光程。
②、计算谱线高斯线宽和洛伦兹线宽的变化范围
根据流场内气体参数(温度、气体压力和吸收分子组分浓度)变化范围，计算所选谱线高斯线宽ΔvD和洛伦兹线宽ΔvC的变化范围，计算公式为
式中，v0是谱线的中心频率，M是吸收分子的摩尔质量，γi(T0)是谱线的展宽系数，T0是参考温度，ni是相应的温度指数。
③、最小二乘方法拟合谐波信号
谱线的谐波信号是关于谱线积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC的函数，以实验测得的谱线的谐波信号为基准，采用最小二乘方法拟合实验测得的谐波信号，获得所选谱线的积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC。
最小二乘拟合步骤为：将积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC的变化范围离散成若干个点，以这三个变量变化范围内的离散点为自由变量，仿真基于不同离散点的所选谱线的所有谐波信号，计算所有谐波信号与实验测得的谐波信号所有点的平方差之和(SSE)，判断是否满足条件(SSE≤ε，ε为很小的数，可设为0.01)，满足，拟合结束，获得A，ΔvD和ΔvC，不满足，更新A，ΔvD和ΔvC，重新拟合直至满足收敛条件。基于最小二乘方法的拟合谐波信号流程见图3。仿真和实验测得的谐波信号对比图见图4。
步骤三、非均匀流场中气体参数计算步骤；
①拟合温度变化范围内所选谱线的线强度随温度的变化
所选谱线的线强度在温度变化范围内具有随温度线性变化的特性，所选谱线的线强度随温度线性变化的示意图见图5。利用线性函数拟合谱线的线强度在该温度范围内的变化特性，拟合公式为
S1(T)＝m1T+b1，S2(T)＝m2T+b2 (3)
式中，S1(T)和S2(T)分别为拟合后两条谱线的线强度，m1、b1、m2和b2分别为常数，可通过线性拟合获得。
②、计算计算非均匀流场内的气体温度和吸收分子的分压。
利用最小二乘方法获得的两条谱线穿过非均匀流场的积分吸光度(A1和A2)计算非均匀流场内的气体温度，计算公式为
式中，温度T的物理含义是吸收分子分压权重的温度积分平均值，其表达式为
利用计算获得的气体温度T，可计算吸收分子的分压，计算公式为
式中，分压PL的物理含义是光路上吸收分子分压的平均值。

Claims (1)

1.基于调制光谱实现非均匀燃烧场气体参数定量测量的方法，其实现步骤为：
步骤一、波长调制光谱测量步骤；
步骤二、波长调制光谱拟合步骤；
步骤三、非均匀流场中气体参数计算步骤；
其特征在于：
步骤一的波长调制光谱测量步骤包括：
(1)预估非均匀流场中气体温度的变化范围，选择具有合适低跃迁态能级的两条谱线，使得所选谱线的线强度在该温度范围内具有随温度近似线性变化的特性；
(2)通过实验测量所选谱线穿过非均匀流场后的透射光强信号随时间的变化；
(3)创建数值锁相放大器，解调实验测得的透射光强信号，得到两条谱线的谐波信号；
步骤二的波长调制光谱拟合步骤包括：
(1)计算非均匀流场中吸收谱线的积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC的变化范围；
(2)以积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC为自变量，仿真所选谱线的透射光强信号随时间的变化；
(3)创建数值锁相放大器，解调仿真的透射光强信号，得到仿真的谐波信号；
(4)基于最小二乘方法拟合实验测得的谐波信号，获得两条谱线穿过非均匀流场后的积分吸光度A，高斯线宽ΔvD和洛伦兹线宽ΔvC；
步骤三的非均匀流场中气体参数计算步骤包括：
(1)利用线性函数拟合气体温度变化范围内所选谱线的线强度随温度的变化；
(2)利用两条谱线的积分吸光度A计算非均匀流场内的气体温度和吸收分子的分压。
CN201710142698.4A 2017-03-10 2017-03-10 基于调制光谱实现非均匀燃烧场气体参数定量测量的方法 Active CN107036983B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710142698.4A CN107036983B (zh) 2017-03-10 2017-03-10 基于调制光谱实现非均匀燃烧场气体参数定量测量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710142698.4A CN107036983B (zh) 2017-03-10 2017-03-10 基于调制光谱实现非均匀燃烧场气体参数定量测量的方法

Publications (2)

Publication Number Publication Date
CN107036983A CN107036983A (zh) 2017-08-11
CN107036983B true CN107036983B (zh) 2018-09-25

Family

ID=59534330

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710142698.4A Active CN107036983B (zh) 2017-03-10 2017-03-10 基于调制光谱实现非均匀燃烧场气体参数定量测量的方法

Country Status (1)

Country Link
CN (1) CN107036983B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112304897B (zh) * 2020-09-14 2021-09-14 中国人民解放军战略支援部队航天工程大学 一种用于燃烧场二维重建的光谱选择方法和系统
CN112255186B (zh) * 2020-10-30 2021-09-28 中国人民解放军战略支援部队航天工程大学 一种非均匀流场积分吸光度的计算方法和系统
CN112798550B (zh) * 2021-04-14 2021-07-13 四川大学 一种宽测量范围的激光吸收光谱燃烧诊断方法
CN113435104B (zh) * 2021-05-21 2024-02-23 中国人民解放军战略支援部队航天工程大学 一种非均匀流场积分吸光度的计算方法和系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216426A (zh) * 2008-01-16 2008-07-09 浙江大学 基于扩展卡尔曼滤波理论的气体状态定量分析仪
CN103323115B (zh) * 2013-06-24 2015-02-11 清华大学 基于波长调制的气体吸收谱线线宽和线型系数的测量方法
CN103868885A (zh) * 2014-03-27 2014-06-18 清华大学 基于复合多次谐波的气体浓度在线测量方法

Also Published As

Publication number Publication date
CN107036983A (zh) 2017-08-11

Similar Documents

Publication Publication Date Title
CN107036983B (zh) 基于调制光谱实现非均匀燃烧场气体参数定量测量的方法
Bao et al. Relative entropy regularized TDLAS tomography for robust temperature imaging
Goldenstein et al. Scanned-wavelength-modulation spectroscopy near 2.5 μm for H 2 O and temperature in a hydrocarbon-fueled scramjet combustor
CN104903703B (zh) 气体吸收分光装置以及气体吸收分光方法
Wang et al. Simultaneous measurement of 2-dimensional H2O concentration and temperature distribution in premixed methane/air flame using TDLAS-based tomography technology
Tuzson et al. Field intercomparison of two optical analyzers for CH 4 eddy covariance flux measurements
CN106949985B (zh) 一种基于混合光抽运的碱金属气室内部温度的精密测量方法
CN107063553B (zh) 一种利用波长调制光谱测量气体压强和组分浓度的装置和方法
CN109100325A (zh) 一种基于光谱吸收率二次谐波特征提取的气体浓度测量方法
CN108760681A (zh) 一种基于波形分解的路径平均温度测量系统与方法
CN105388120B (zh) 基于wmrf模型的免定标波长调制光谱气体检测方法
CN108709871B (zh) 一种基于tdlas的气体浓度检测方法
CN108981953A (zh) 一种基于干涉调制原理的激光吸收光谱测温方法和系统
CN113447458B (zh) 一种基于激光吸收阻抗谱的气体温度与浓度参数测量方法
CN107389606A (zh) 一种基于可调谐半导体激光吸收谱的光程长度分析方法
CN107219190A (zh) 一种基于三角波调制的免标定气体参数测量方法
CN106289712B (zh) 一种内流道阻力测量方法
Xia et al. Probing greenhouse gases in turbulent atmosphere by long-range open-path wavelength modulation spectroscopy
Xu et al. Noise immune TDLAS temperature measurement through spectrum shifting by using a Mach–Zehnder interferometer
Lan et al. Multi-harmonic measurements of line shape under low absorption conditions
Li et al. Robust cepstral analysis at variable wavelength scan depth for narrowband tunable laser absorption spectroscopy
CN107941467A (zh) 直接获取分布反馈半导体激光源电流调制波长响应的方法
Jiang et al. TDLAS-WMS second harmonic detection based on spectral analysis
Zhang et al. Simultaneous detection of multiple gas concentrations with multi-frequency wavelength modulation spectroscopy
Du et al. Measurement of the absolute absorbance based on wavelength modulation spectroscopy

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 101416 3380 mailbox No. 86, equipment college, Bayi Road, Huairou District, Beijing, No. 86

Applicant after: China People's Liberation Army Strategic Support Unit Space Engineering University

Address before: 101416 3380 mailbox No. 86, equipment college, Bayi Road, Huairou District, Beijing, No. 86

Applicant before: PLA 'S Equipment College

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant