CN114877951B - 一种燃烧温度场及水蒸气含量同步测量装置及方法 - Google Patents

一种燃烧温度场及水蒸气含量同步测量装置及方法 Download PDF

Info

Publication number
CN114877951B
CN114877951B CN202210710542.2A CN202210710542A CN114877951B CN 114877951 B CN114877951 B CN 114877951B CN 202210710542 A CN202210710542 A CN 202210710542A CN 114877951 B CN114877951 B CN 114877951B
Authority
CN
China
Prior art keywords
terahertz wave
light intensity
dimensional light
intensity information
temperature field
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
CN202210710542.2A
Other languages
English (en)
Other versions
CN114877951A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN202210710542.2A priority Critical patent/CN114877951B/zh
Publication of CN114877951A publication Critical patent/CN114877951A/zh
Application granted granted Critical
Publication of CN114877951B publication Critical patent/CN114877951B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00Measuring or testing not otherwise provided for
    • G01D21/02Measuring two or more variables by means not covered by a single other subclass
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/0014Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation from gases, flames
    • G01J5/0018Flames, plasma or welding
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/10Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors
    • G01J5/34Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors using capacitors, e.g. pyroelectric capacitors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3504Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3581Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3504Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
    • G01N2021/354Hygrometry of gases
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Power Engineering (AREA)
  • Plasma & Fusion (AREA)
  • Toxicology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本申请提供了一种燃烧温度场及水蒸气含量同步测量装置及方法,包括:第一测量模块、第二测量模块、第三测量模块和数据处理模块;第一测量模块用于产生未穿过火焰的参考太赫兹波以及穿过火焰的两束测量太赫兹波,将参考太赫兹波和一束测量太赫兹波分别投射到第二测量模块和第三测量模块,采集另一束的二维光强信息;第二测量模块用于将参考太赫兹波分为两束,将一束投射到第三测量模块,采集另一束的二维光强信息;第三测量模块用于生成太赫兹波和参考太赫兹波的干涉条纹图像;数据处理模块用于对两个二维光强信息和干涉条纹图像进行处理,得到火焰的温度场分布和水蒸气浓度分布。本申请实现了火焰温度和水蒸气含量的同步测量。

Description

一种燃烧温度场及水蒸气含量同步测量装置及方法
技术领域
本申请涉及光学燃烧测量技术领域,尤其是涉及一种燃烧温度场及水蒸气含量同步测量装置及方法。
背景技术
温度和组分测量是燃烧诊断的重要监测参数,其中,温度是表征热释放最直接的物理量,而水蒸气作为燃烧中最主要的产物,其浓度是表征燃烧效率的重要参量,因此燃烧温度和水蒸气浓度测量对于燃烧状态的判断、预测和诊断具有非常重要的意义。
传统的热电偶或烟气分析仪等接触式测量手段虽然简单、直观,但存在感应元件易于损坏、被测场易受干扰等难题。随着激光技术的发展,灵敏度很高、测量响应速度快的激光光谱技术成为燃烧诊断的重要手段之一,其中激光干涉测温是基于光波的相位移特性,通过燃烧场的干涉图来反映温度场的空间分布;可调谐二极管激光吸收光谱(TDLAS)技术基于燃烧粒子的吸收光谱获取具燃烧温度或组分浓度信息。但现有的激光干涉或TDLAS技术还不能满足对火焰温度和水蒸气组分浓度的同时在线测量的需求。
发明内容
有鉴于此,本申请提供了一种燃烧温度场及水蒸气含量同步测量装置及方法,以解决上述技术问题。
第一方面,本申请实施例提供了一种燃烧温度场及水蒸气含量同步测量装置,包括:第一测量模块、第二测量模块、第三测量模块和数据处理模块;
所述第一测量模块,用于产生未穿过火焰样品的参考太赫兹波以及穿过火焰样品的两束测量太赫兹波,将参考太赫兹波投射到第二测量模块,将一束测量太赫兹波投射到第三测量模块,采集另一束测量太赫兹波的二维光强信息,发送至数据处理模块;
所述第二测量模块,用于将参考太赫兹波分为两束,将一束参考太赫兹波投射到所述第三测量模块,采集另一束参考太赫兹波的二维光强信息,发送至数据处理模块;
所述第三测量模块,用于实现测量太赫兹波和参考太赫兹波的干涉,生成干涉条纹图像,发送至数据处理模块;
所述数据处理模块,用于对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布。
进一步,所述第一测量模块包括:自由电子激光器,光隔离器,第一离轴抛物面镜,第二离轴抛物面镜,第一分束镜,第二分束镜和第一热释电探测器;所述火焰样品位于第一分束镜和第二分束镜之间;
所述自由电子激光器,用于发射连续太赫兹波;
所述光隔离器的中心轴与所述自由电子激光器发射的太赫兹波平行,用于隔离反射的太赫兹回波,使太赫兹波单方向传播;
所述第一离轴抛物面镜位于太赫兹波的光路上,所述第二离轴抛物面镜与所述第一离轴抛物面镜相对放置;所述第一离轴抛物面镜和第二离轴抛物面镜对太赫兹波进行反射准直,实现对太赫兹波的扩束;
所述第一分束镜位于第二离轴抛物面镜反射的光路上,用于将扩束后的太赫兹波分为两束:一束透射的测量太赫兹波穿过火焰样品到达第二分束镜;另一束反射的参考太赫兹波进入所述第二测量模块;
所述第二分束镜用于将测量太赫兹波分为两束,一束透射的测量太赫兹波进入第一热释电探测器,另一束反射的测量太赫兹波进入第三测量模块;
所述第一热释电探测器为阵列式探测器,用于采集透射的测量太赫兹波的二维光强信息,发送至数据处理模块。
进一步,所述第一测量模块还包括:速度调制器,用于调制火焰样品的速度。
进一步,所述第二测量模块包括:第三分束镜和第二热释电探测器;
所述第三分束镜位于参考太赫兹波的光路上,用于将参考太赫兹波分为两束,一束透射的参考太赫兹波进入第二热释电探测器,另一束反射的参考太赫兹波进入第三测量模块;
所述第二热释电探测器为阵列式探测器,用于测量透射的参考太赫兹波的二维光强信息,发送至数据处理模块。
进一步,所述第三测量模块包括:第四分束镜、第三离轴抛物面镜、第四离轴抛物面镜和CCD相机;
所述第四分束镜位于第二分束镜的反射的测量太赫兹波光路,和第三分束镜的反射的参考太赫兹波的光路上,用于实现测量太赫兹波和参考太赫兹波的干涉;
所述第四离轴抛物面镜与第三离轴抛物面镜相对放置,所述第三离轴抛物面镜和第四离轴抛物面镜对干涉后的太赫兹波进行反射准直,实现干涉后的太赫兹波的缩束,将缩束的干涉后的太赫兹波投射至所述CCD相机;
所述CCD相机,用于对缩束的干涉后的太赫兹波进行成像,将干涉条纹图像发送至数据处理模块。
进一步,所述数据处理模块具体用于:
获取测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像;
对干涉条纹图像进行处理,得到火焰样品的温度场分布;
对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布。
进一步,对干涉条纹图像进行处理,得到火焰样品的温度场分布,得到火焰样品的温度场分布;包括:
对干涉条纹图像进行滤波处理,去除燃烧过程中杂散光对干涉条纹图像的噪声干扰,得到去噪后的干涉条纹图像;
对去噪后的干涉条纹图像进行图像二值化处理,得到二值化图像;
采用数学形态学对二值化图像进行细化,提取出干涉条纹中心线,标记干涉条纹级数N;
基于干涉条纹级数,计算火焰样品的折射率分布n(x,y):
Figure BDA0003706714260000041
其中,nref为空气折射率,λ为太赫兹波波长,L为太赫兹波通过火焰区域宽度;
根据火焰样品的折射率分布n(x,y)、理想气体定律及Gladstone-Dale关系,计算火焰样品区域的温度场分布T(x,y):
Figure BDA0003706714260000042
其中,Tref是参考温度。
进一步,对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布;包括:
构建多项式形式的配分函数Q(T),多项式系数由HITRAN数据库获得,T为温度;
温度T下水蒸气吸收谱线的谱线强度S(T)为:
Figure BDA0003706714260000051
其中,S(T0)为预设温度T0的谱线强度,k为玻尔兹曼常数,h为普朗克常数,c为光速,E″为分子低能态的能量,υ0为吸收谱线中心频率。
基于测量太赫兹波的二维光强I(x,y)、参考太赫兹波的二维光强I0(x,y)、火焰样品折射率分布n(x,y)和温度场分布T(x,y),利用理想气体定律及Gladstone-Dale关系计算火焰中水蒸气的浓度分布X(x,y):
Figure BDA0003706714260000052
其中,G为混合气体的Gladstone-Dale系数,R为理想气体常数,S(T(x,y))是温度T(x,y)的水蒸气吸收谱线的谱线强度,
Figure BDA0003706714260000053
是洛伦兹型线型函数:
Figure BDA0003706714260000054
其中,Δν为碰撞线宽;υ为频率。
第二方面,本申请实施例提供了一种燃烧温度场及水蒸气含量同步测量方法,应用于本申请实施例的装置,包括:
打开所述自由电子激光器,所述自由电子激光器发出连续太赫兹波;
所述第一热释电探测器采集测量太赫兹波的二维光强信息,发送至数据处理模块;所述第二热释电探测器采集参考太赫兹波的二维光强信息,发送至数据处理模块;所述CCD相机生成干涉条纹图像,发送至数据处理模块;
所述数据处理模块对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布。
进一步,所述数据处理模块对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布;包括:
获取测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像;
对干涉条纹图像进行处理,得到火焰样品的温度场分布;
对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布。
本申请实现了火焰温度和水蒸气含量的同步测量。
附图说明
为了更清楚地说明本申请具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本申请的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的燃烧温度场及水蒸气含量同步测量装置的结构示意图;
图2为本申请实施例提供的燃烧温度及水蒸气含量同步测量方法流程图。
图标:
101-自由电子激光器,102-光隔离器,
103-第一离轴抛物面镜,104-第二离轴抛物面镜,
105-第三离轴抛物面镜,106-第四离轴抛物面镜,
107-第一分束镜,108-第二分束镜,109-第三分束镜,
110-第四分束镜,111-火焰样品,112-第一热释电探测器,
113-第二热释电探测器,114-CCD相机,115-电子设备。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。
因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
首先对本申请实施例的设计思想进行简单介绍。
太赫兹波是一种处于红外和微波间的过渡波段,对于很多介电材料和非极性的液体有良好的穿透性,在安全监测、无损检测、药品检测等领域有广泛应用。作为目前电磁空间唯一尚未完全开发利用的频谱资源,太赫兹频段涵盖了大多数燃烧气体分子的纯转动跃迁频率,具有对气体分子温度及含量变化敏感的特点;太赫兹波波长较长,在非均匀物质中有较少的散射和偏折;同时太赫兹辐射对红外辐射噪声不敏感;这些独特的优势使得太赫兹波在燃烧诊断领域展现了一定的应用前景。
针对现有的激光干涉或TDLAS技术还不能满足对火焰温度和水蒸气组分浓度的同时在线测量的需求,基于太赫兹技术,本申请提供了一种燃烧温度场及水蒸气含量同步测量装置,利用分束镜替换现有马赫-曾德干涉光路中的反射镜,以同时获取燃烧场的干涉图像及透射光谱;通过干涉图像确定相位信息,进而重构火焰的温度场,基于燃烧的气体吸收谱线获取燃烧产物中水蒸气浓度信息,以此解决现有技术中存在的对于火焰温度和水蒸气组分浓度的同步在线测量困难的问题。
在介绍了本申请实施例的应用场景和设计思想之后,下面对本申请实施例提供的技术方案进行说明。
如图1所示,本申请实施例提供了一种燃烧温度场及水蒸气含量同步测量装置,包括:第一测量模块、第二测量模块、第三测量模块和数据处理模块;所述数据处理模块设置在电子设备115上;
所述第一测量模块,用于产生未穿过火焰样品的参考太赫兹波以及穿过火焰样品的两束测量太赫兹波,将参考太赫兹波投射到第二测量模块,将一束测量太赫兹波投射到第三测量模块,采集另一束测量太赫兹波的二维光强信息,发送至数据处理模块;
在本实施例中,所述第一测量模块包括:自由电子激光器101,光隔离器102,第一离轴抛物面镜103,第二离轴抛物面镜104,第一分束镜107,第二分束镜108和第一热释电探测器112;所述火焰样品111位于第一分束镜107和第二分束镜108之间;所述第一测量模块还包括:速度调制器,用于调制火焰样品111的速度。
所述自由电子激光器101,用于发射连续太赫兹波;
所述光隔离器102的中心轴与自由电子激光器101发射的太赫兹波平行,用于隔离反射的太赫兹回波,使太赫兹波单方向传播;
所述第一离轴抛物面镜103位于太赫兹波的光路上,所述第二离轴抛物面镜104与第一离轴抛物面镜103相对放置;所述第一离轴抛物面镜103和第二离轴抛物面镜104对太赫兹波进行反射准直,实现对太赫兹波的扩束;
所述第一分束镜107位于第二离轴抛物面镜104反射的光路上,用于将扩束后的太赫兹波分为两束:一束透射的测量太赫兹波穿过火焰样品111到达第二分束镜108;另一束反射的参考太赫兹波进入所述第二测量模块;
所述第二分束镜108用于将测量太赫兹波分为两束,一束透射的测量太赫兹波进入第一热释电探测器112,另一束反射的测量太赫兹波进入第三测量模块;
所述第一热释电探测器112为阵列式探测器,用于采集透射的测量太赫兹波的二维光强信息,发送至数据处理模块。
所述第二测量模块,用于将参考太赫兹波分为两束,将一束参考太赫兹波投射到所述第三测量模块,采集另一束参考太赫兹波的二维光强信息,发送至数据处理模块;
在本实施例中,所述第二测量模块包括:第三分束镜109和第二热释电探测器113;
所述第三分束镜109位于参考太赫兹波的光路上,用于将参考太赫兹波分为两束,一束透射的参考太赫兹波进入第二热释电探测器113,另一束反射的参考太赫兹波进入第三测量模块;
所述第二热释电探测器113为阵列式探测器,用于测量透射的参考太赫兹波的二维光强信息,发送至数据处理模块。
所述第三测量模块,用于实现测量太赫兹波和参考太赫兹波的干涉,生成干涉条纹图像,发送至数据处理模块;
所述第三测量模块包括:第四分束镜110、第三离轴抛物面镜105、第四离轴抛物面镜106和CCD相机114;
在本实施例中,所述第四分束镜110位于第二分束镜108的反射的测量太赫兹波光路,和第三分束镜109反射的参考太赫兹波的光路上,用于实现测量太赫兹波和参考太赫兹波的干涉;
所述第四离轴抛物面镜106与第三离轴抛物面镜105相对放置,共同对干涉后的太赫兹波进行反射准直,实现干涉后的太赫兹波的缩束,将缩束的干涉后的太赫兹波投射至所述CCD相机114;
所述CCD相机114,用于对缩束的干涉后的太赫兹波进行成像,将干涉条纹图像发送至数据处理模块。
所述数据处理模块,用于获取测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像;对干涉条纹图像进行处理,得到火焰样品的温度场分布;对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布。
其中,对干涉条纹图像进行处理,得到火焰样品的温度场分布,得到火焰样品的温度场分布;包括:
基于偏微分方程的方法对干涉条纹图像进行滤波处理,以滤掉燃烧过程中杂散光对干涉条纹图像的噪声干扰,得到去噪后的干涉条纹图像;
对去噪后的干涉条纹图像进行图像二值化处理,得到二值化图像;
在本实施例中,采用具有精度较高和自动化程度高的导数符号二值图法进行图像二值化。
采用数学形态学对二值化图像进行细化,提取出干涉条纹中心线,标记干涉条纹级数N;
基于干涉条纹级数,得到测量太赫兹波与参考太赫兹波间相位差
Figure BDA0003706714260000101
计算火焰样品的折射率分布n(x,y):
Figure BDA0003706714260000111
其中,nref为空气折射率,λ为太赫兹波波长,L为太赫兹波通过火焰区域宽度;
根据火焰样品的折射率分布n(x,y)、理想气体定律及Gladstone-Dale关系,计算火焰样品区域的温度场分布T(x,y):
Figure BDA0003706714260000112
其中,Tref是参考温度。
其中,对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布;包括:
构建多项式形式的配分函数Q(T),多项式系数由HITRAN数据库获得,T为温度;
温度T下水蒸气吸收谱线的谱线强度S(T)为:
Figure BDA0003706714260000113
其中,S(T0)为预设温度T0的谱线强度,k为玻尔兹曼常数,h为普朗克常数,c为光速,E″为分子低能态的能量,υ0为吸收谱线中心频率。
基于测量太赫兹波的二维光强I(x,y)、参考太赫兹波的二维光强I0(x,y)、火焰样品折射率分布n(x,y)和温度场分布T(x,y),利用理想气体定律及Gladstone-Dale关系计算火焰中水蒸气的浓度分布X(x,y):
Figure BDA0003706714260000114
其中,G为混合气体的Gladstone-Dale系数,R为理想气体常数,S(T(x,y))是温度T(x,y)的水蒸气吸收谱线的谱线强度,
Figure BDA0003706714260000121
是洛伦兹型线型函数:
Figure BDA0003706714260000122
其中,Δν为碰撞线宽;υ为频率。
基于上述实施例,本申请实施例提供了一种燃烧温度场及水蒸气含量同步测量方法,包括如下步骤:
步骤201:打开所述自由电子激光器,所述自由电子激光器发出连续太赫兹波;
步骤202:所述第一热释电探测器采集测量太赫兹波的二维光强信息,发送至数据处理模块;所述第二热释电探测器采集另一束参考太赫兹波的二维光强信息,发送至数据处理模块;所述CCD相机生成干涉条纹图像,发送至数据处理模块;
步骤203:所述数据处理模块对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布。
在本实施例中,该步骤具体包括:
获取测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像;
对干涉条纹图像进行处理,得到火焰样品的温度场分布;
对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布。
此外,尽管在附图中以特定顺序描述了本申请方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
尽管已描述了本申请的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
最后应说明的是:以上各实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述各实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。

Claims (8)

1.一种燃烧温度场及水蒸气含量同步测量装置,其特征在于,包括:第一测量模块、第二测量模块、第三测量模块和数据处理模块;
所述第一测量模块,用于产生未穿过火焰样品的参考太赫兹波以及穿过火焰样品的两束测量太赫兹波,将参考太赫兹波投射到第二测量模块,将一束测量太赫兹波投射到第三测量模块,采集另一束测量太赫兹波的二维光强信息,发送至数据处理模块;
所述第二测量模块,用于将参考太赫兹波分为两束,将一束参考太赫兹波投射到所述第三测量模块,采集另一束参考太赫兹波的二维光强信息,发送至数据处理模块;
所述第三测量模块,用于实现测量太赫兹波和参考太赫兹波的干涉,生成干涉条纹图像,发送至数据处理模块;
所述数据处理模块,用于对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布;
所述数据处理模块具体用于:
获取测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像;
对干涉条纹图像进行处理,得到火焰样品的温度场分布;
对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布;
其中,对干涉条纹图像进行处理,得到火焰样品的温度场分布,包括:对干涉条纹图像进行滤波处理,去除燃烧过程中杂散光对干涉条纹图像的噪声干扰,得到去噪后的干涉条纹图像;
对去噪后的干涉条纹图像进行图像二值化处理,得到二值化图像;
采用数学形态学对二值化图像进行细化,提取出干涉条纹中心线,标记干涉条纹级数N;
基于干涉条纹级数,计算火焰样品的折射率分布n(x,y):
Figure FDA0004181446260000021
其中,nref为空气折射率,λ为太赫兹波波长,L为太赫兹波通过火焰区域宽度;
根据火焰样品的折射率分布n(x,y)、理想气体定律及Gladstone-Dale关系,计算火焰样品区域的温度场分布T(x,y):
Figure FDA0004181446260000022
其中,Tref是参考温度。
2.根据权利要求1所述的燃烧温度场及水蒸气含量同步测量装置,其特征在于,所述第一测量模块包括:自由电子激光器,光隔离器,第一离轴抛物面镜,第二离轴抛物面镜,第一分束镜,第二分束镜和第一热释电探测器;所述火焰样品位于第一分束镜和第二分束镜之间;
所述自由电子激光器,用于发射连续太赫兹波;
所述光隔离器的中心轴与所述自由电子激光器发射的太赫兹波平行,用于隔离反射的太赫兹回波,使太赫兹波单方向传播;
所述第一离轴抛物面镜位于太赫兹波的光路上,所述第二离轴抛物面镜与所述第一离轴抛物面镜相对放置;所述第一离轴抛物面镜和第二离轴抛物面镜对太赫兹波进行反射准直,实现对太赫兹波的扩束;
所述第一分束镜位于第二离轴抛物面镜反射的光路上,用于将扩束后的太赫兹波分为两束:一束透射的测量太赫兹波穿过火焰样品到达第二分束镜;另一束反射的参考太赫兹波进入所述第二测量模块;
所述第二分束镜用于将测量太赫兹波分为两束,一束透射的测量太赫兹波进入第一热释电探测器,另一束反射的测量太赫兹波进入第三测量模块;
所述第一热释电探测器为阵列式探测器,用于采集透射的测量太赫兹波的二维光强信息,发送至数据处理模块。
3.根据权利要求2所述的燃烧温度场及水蒸气含量同步测量装置,其特征在于,所述第一测量模块还包括:速度调制器,用于调制火焰样品的速度。
4.根据权利要求2所述的燃烧温度场及水蒸气含量同步测量装置,其特征在于,所述第二测量模块包括:第三分束镜和第二热释电探测器;
所述第三分束镜位于参考太赫兹波的光路上,用于将参考太赫兹波分为两束,一束透射的参考太赫兹波进入第二热释电探测器,另一束反射的参考太赫兹波进入第三测量模块;
所述第二热释电探测器为阵列式探测器,用于测量透射的参考太赫兹波的二维光强信息,发送至数据处理模块。
5.根据权利要求4所述的燃烧温度场及水蒸气含量同步测量装置,其特征在于,所述第三测量模块包括:第四分束镜、第三离轴抛物面镜、第四离轴抛物面镜和CCD相机;
所述第四分束镜位于第二分束镜的反射的测量太赫兹波光路,和第三分束镜的反射的参考太赫兹波的光路上,用于实现测量太赫兹波和参考太赫兹波的干涉;
所述第四离轴抛物面镜与第三离轴抛物面镜相对放置,所述第三离轴抛物面镜和第四离轴抛物面镜对干涉后的太赫兹波进行反射准直,实现干涉后的太赫兹波的缩束,将缩束的干涉后的太赫兹波投射至所述CCD相机;
所述CCD相机,用于对缩束的干涉后的太赫兹波进行成像,将干涉条纹图像发送至数据处理模块。
6.根据权利要求5所述的燃烧温度场及水蒸气含量同步测量装置,其特征在于,对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布;包括:
构建多项式形式的配分函数Q(T),多项式系数由HITRAN数据库获得,T为温度;
温度T下水蒸气吸收谱线的谱线强度S(T)为:
Figure FDA0004181446260000041
其中,S(T0)为预设温度T0的谱线强度,k为玻尔兹曼常数,h为普朗克常数,c为光速,E″为分子低能态的能量,υ0为吸收谱线中心频率;
基于测量太赫兹波的二维光强I(x,y)、参考太赫兹波的二维光强Io(x,y)、火焰样品折射率分布n(x,y)和温度场分布T(x,y),利用理想气体定律及Gladstone-Dale关系计算火焰中水蒸气的浓度分布X(x,y):
Figure FDA0004181446260000042
其中,G为混合气体的Gladstone-Dale系数,R为理想气体常数,S(T(x,y))是温度T(x,y)的水蒸气吸收谱线的谱线强度,
Figure FDA0004181446260000043
是洛伦兹型线型函数:
Figure FDA0004181446260000044
其中,Δv为碰撞线宽;υ为频率。
7.一种燃烧温度场及水蒸气含量同步测量方法,应用于权利要求5-6任一项所述的燃烧温度场及水蒸气含量同步测量装置,其特征在于,包括:
打开所述自由电子激光器,所述自由电子激光器发出连续太赫兹波;
所述第一热释电探测器采集测量太赫兹波的二维光强信息,发送至数据处理模块;所述第二热释电探测器采集参考太赫兹波的二维光强信息,发送至数据处理模块;所述CCD相机生成干涉条纹图像,发送至数据处理模块;
所述数据处理模块对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布。
8.根据权利要求7所述的燃烧温度场及水蒸气含量同步测量方法,其特征在于,所述数据处理模块对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像进行处理,得到火焰样品的温度场分布和水蒸气浓度分布;包括:
获取测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和干涉条纹图像;
对干涉条纹图像进行处理,得到火焰样品的温度场分布;
对测量太赫兹波的二维光强信息、参考太赫兹波的二维光强信息和火焰样品的温度场分布进行处理,得到火焰样品的水蒸气浓度分布。
CN202210710542.2A 2022-06-22 2022-06-22 一种燃烧温度场及水蒸气含量同步测量装置及方法 Active CN114877951B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210710542.2A CN114877951B (zh) 2022-06-22 2022-06-22 一种燃烧温度场及水蒸气含量同步测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210710542.2A CN114877951B (zh) 2022-06-22 2022-06-22 一种燃烧温度场及水蒸气含量同步测量装置及方法

Publications (2)

Publication Number Publication Date
CN114877951A CN114877951A (zh) 2022-08-09
CN114877951B true CN114877951B (zh) 2023-07-04

Family

ID=82681769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210710542.2A Active CN114877951B (zh) 2022-06-22 2022-06-22 一种燃烧温度场及水蒸气含量同步测量装置及方法

Country Status (1)

Country Link
CN (1) CN114877951B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013024842A (ja) * 2011-07-26 2013-02-04 Hitachi High-Tech Control Systems Corp テラヘルツ波イメージング装置
CN103919530A (zh) * 2014-04-21 2014-07-16 首都师范大学 一种增强生物组织太赫兹波成像信号强度的***和方法
JP6266719B1 (ja) * 2016-08-31 2018-01-24 フェムトディプロイメンツ株式会社 テラヘルツ時間領域分光装置
CN109781656B (zh) * 2018-12-27 2020-12-04 深圳市华讯方舟太赫兹科技有限公司 基于太赫兹的水蒸气检测***及检测方法
CN211235536U (zh) * 2019-12-31 2020-08-11 鞍山亚赛电磁设备有限公司 一种燃烧场温度与气体组分浓度测试仪
CN111077109B (zh) * 2020-01-16 2021-06-11 北京航空航天大学 一种基于双光梳光谱技术的温度和浓度测量***和方法
CN114061783B (zh) * 2022-01-18 2022-04-26 北京航空航天大学 基于马赫曾德干涉的高频热释放率脉动场测量装置及方法

Also Published As

Publication number Publication date
CN114877951A (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
US9182282B2 (en) Multi-analyte optical computing system
US8213006B2 (en) Multi-analyte optical computing system
US7990538B2 (en) Signal processing for optical computing system
US8213012B2 (en) Stability for optical computing system
RU2400715C2 (ru) Способ градуировки спектрометра
US7911605B2 (en) Multivariate optical elements for optical analysis system
US7834999B2 (en) Optical analysis system and optical train
JP6053523B2 (ja) ファブリーペローフーリエ変換分光計
US8862445B2 (en) Selecting spectral elements and components for optical analysis systems
US8039803B2 (en) Phase contrast imaging method and apparatus
US20070165210A1 (en) High-density channels detecting device
CN110927121A (zh) 一种基于白光干涉光谱的相位型spr检测装置及方法
Luo et al. 4D surface shape measurement system with high spectral resolution and great depth accuracy
US5406377A (en) Spectroscopic imaging system using a pulsed electromagnetic radiation source and an interferometer
Rauscher et al. Dual-channel mid-infrared sensor based on tunable Fabry-Pérot filters for fluid monitoring applications
CN114877951B (zh) 一种燃烧温度场及水蒸气含量同步测量装置及方法
CN110658144B (zh) 一种气体成像装置
JPH0414298B2 (zh)
US8223339B1 (en) Scanning spectroscopy modulation for frequency domain spectral analysis
Tasumi Introduction to infrared spectroscopy
Delwiche Basics of spectroscopic analysis
JP2009019893A (ja) センシング方法及びセンシング装置
US20240142372A1 (en) Adulteration and authenticity analysis method of organic substances and materials by terahertz spectroscopy
WO1996000887A1 (en) An improved optical sensor and method
Vickers et al. On-line monitoring by Raman spectroscopy: instrument control and calibration

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
GR01 Patent grant
GR01 Patent grant