CN105547469B - 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 - Google Patents
基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 Download PDFInfo
- Publication number
- CN105547469B CN105547469B CN201510889073.5A CN201510889073A CN105547469B CN 105547469 B CN105547469 B CN 105547469B CN 201510889073 A CN201510889073 A CN 201510889073A CN 105547469 B CN105547469 B CN 105547469B
- Authority
- CN
- China
- Prior art keywords
- flame
- mrow
- radiation
- msub
- objective function
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 40
- 230000003287 optical effect Effects 0.000 title abstract 4
- 230000005855 radiation Effects 0.000 claims abstract description 132
- 230000003595 spectral effect Effects 0.000 claims abstract description 68
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 230000000704 physical effect Effects 0.000 claims description 32
- 230000005540 biological transmission Effects 0.000 claims description 18
- 230000001052 transient effect Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 5
- 238000000295 emission spectrum Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000004141 dimensional analysis Methods 0.000 abstract 1
- 238000009529 body temperature measurement Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000002485 combustion reaction Methods 0.000 description 3
- 230000005457 Black-body radiation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003889 chemical engineering Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000009851 ferrous metallurgy Methods 0.000 description 1
- 238000002309 gasification Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 229910002059 quaternary alloy Inorganic materials 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J1/00—Photometry, e.g. photographic exposure meter
- G01J1/42—Photometry, e.g. photographic exposure meter using electric radiation detectors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Radiation Pyrometers (AREA)
Abstract
基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,涉及高温热辐射测量技术领域。本发明是为了适应火焰光场探测中对于泛尺度分析的需求。本发明所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,利用具有微透镜阵列的光场相机获取高温火焰不同方向上的辐射强度信息,通过将脉冲激光照射到火焰上,由探测器接收到的介质边界上的出射光谱辐射强度通过逆问题求解得到介质的光谱辐射特性参数,结合火焰自身的出射辐射强度重建出高温火焰的三维温度场。通过本发明的仿真计算,可实现对高温火焰的辐射特性参数和温度场的重建,并为光场相机的标定、测量等工作提供理论基础。
Description
技术领域
本发明属于高温热辐射测量技术领域。
背景技术
高温燃烧现象广泛存在于航空航天、能源动力、钢铁冶金、化工等领域,如火箭发动机、燃气轮机、内燃机、电站锅炉、煤气化反应器等高温设备中。火焰温度场的测量对于燃烧诊断具有十分重要的意义。常用的温度场测量方法主要包括接触式测温和非接触式测温两类,其中非接触式测温方法因为其不干扰流场、测量精度高、实时连续测量等优点而得到了越来越广泛的应用。光场成像技术是一种典型的非接触式测温技术,可通过一次成像获取三维全场辐射数字化信息。
对于非接触式测温技术,其关键问题在于建立高效、准确、稳定的介质温度场重建算法。参与性介质温度场反演属于典型的不适定问题,尤其是当介质的辐射特性参数未知时,需要同时反演介质的温度场和辐射特性参数,这给问题求解带来了很大困难。
利用反演的方法测量半透明介质的辐射物性,通常需要利用半透明材料表面的扩散光信号,目前对扩散光信号的测量主要有三种方法:时域测量、频域测量和稳态测量。时域测量由于能够提供较多的介质信息,具有很高的信噪比和灵敏度,在半透明介质物性测量中具有广阔的应用前景。
火焰光场探测中的热辐射传输过程,会涉及到空间、时间、角度、频率四个变量的多个数量级变化,如:对于采用微透镜阵列技术的光场相机而言,由于被探测火焰与相机镜头之间巨大的几何尺度差异,通过微透镜进行探测角度离散后的探测像元所代表的探测角度范围与相机自身的探测角度范围又存在着两个数量级以上的差距,此时的泛尺度分析就十分必要。
发明内容
本发明是为了适应火焰光场探测中对于泛尺度分析的需求,现提供基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法。
基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,该方法是基于以下装置实现的,所属装置包括:激光控制器1、激光头2、光电探测器4、数据采集处理系统5和微透镜阵列光场相机6;
激光控制器1的激光控制信号输出端连接激光头2的激光控制信号输入端,光电探测器4均匀分布在火焰3周围,光电探测器4的光电探测信号输出端连接数据采集处理系统5的光电探测信号输入端,微透镜阵列光场相机6的断层三维辐射场信号输出端连接数据采集处理系统5的断层三维辐射场信号输入端;
所述方法包括以下步骤:
步骤一:开启激光控制器1,使激光头2发射出的脉冲激光入射到火焰3上,然后使激光头2顺时针旋转,且每次旋转的角度均为30°,直至激光头2旋转回初始位置,激光头2每旋转一次则利用光电探测器4采集一次火焰3透射或反射出的激光,然后将光电探测器4获得的所有光电探测信号发送至数据采集处理系统5中,数据采集处理系统5分别对其获得的光电探测信号进行处理,获得火焰3边界的出射光谱辐射强度值Imi,i为激光头2旋转次数,i=1,2,...,12;
步骤二:关闭激光控制器1,保持火焰3稳定,利用微透镜阵列光场相机6从单角度成像获取火焰3的断层三维辐射场信号,然后将该断层三维辐射场信号发送至数据采集处理系统5,数据采集处理系统5对断层三维辐射场信号进行处理,获得火焰3自身的出射光谱辐射强度测量值In;
步骤三:将步骤一获得的火焰3边界的出射光谱辐射强度值Imi分别与步骤二获得的火焰3自身的出射光谱辐射强度测量值In作差,获得消除了自辐射的光谱辐射强度测量值Mi;
步骤四:根据逆问题算法假设火焰光谱辐射物性参数为k,将该火焰光谱辐射物性参数k代入瞬态辐射传输方程,并对该瞬态辐射传输方程求解获得消除了自辐射的光谱辐射强度估计值Ei;
步骤五:以步骤三获得的光谱辐射强度测量值Mi与步骤四获得的光谱辐射强度估计值Ei的最小二乘作为火焰光谱辐射物性参数目标函数方程,并求解该火焰光谱辐射物性参数目标函数方程的目标函数值,然后将该目标函数值与火焰光谱辐射物性参数的给定目标函数阈值进行比较,判断目标函数值是否小于给定目标函数阈值,是则步骤四中假设的火焰光谱辐射物性参数k为真实的火焰光谱辐射物性参数,然后执行步骤六,否则返回步骤四;
步骤六:利用真实的火焰光谱辐射物性参数和逆问题算法假设火焰温度场为T,将假设火焰温度场T代入火焰3的温度场辐射传输方程,并对火焰3的温度场辐射传输方程求解获得火焰3自身的出射辐射强度估计值I'n;
步骤七:利用步骤二获得的出射光谱辐射强度测量值In和步骤六获得的出射辐射强度估计值I'n构造温度场目标函数方程,并求解该温度场目标函数方程的温度场目标函数值,然后将该温度场目标函数值与温度场给定目标函数阈值进行比较,判断该温度场目标函数值是否小于温度场给定目标函数阈值,是则步骤六中假设的火焰温度场T为火焰3的真实温度场,完成火焰温度泛尺度光场探测,否则返回步骤六。
本发明涉及一种四元体系的泛尺度分析方法,具体地说是一种基于微透镜阵列与脉冲激光的空间——时间——角度——频率四种变量有多个数量级变化时的高温火焰光谱辐射特性参数及温度场的重建方法。本发明所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,主要包括:高温火焰辐射测量、高温火焰辐射特性参数计算、温度场重建等环节。利用具有微透镜阵列的光场相机获取高温火焰不同方向上的辐射强度信息,通过将脉冲激光照射到火焰上,由探测器接收到的介质边界上的出射光谱辐射强度通过逆问题求解得到介质的光谱辐射特性参数,结合火焰自身的出射辐射强度重建出高温火焰的三维温度场。通过本发明的仿真计算,可实现对高温火焰的辐射特性参数和温度场的重建,并为光场相机的标定、测量等工作提供理论基础。
附图说明
图1为实现步骤一时所采用的硬件的结构示意图;
图2为实现步骤二时所采用的硬件的结构示意图;
图3为基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法的流程图。
具体实施方式
具体实施方式一:参照图1、图2和图3具体说明本实施方式,本实施方式所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,该方法是基于以下装置实现的,所属装置包括:激光控制器1、激光头2、光电探测器4、数据采集处理系统5和微透镜阵列光场相机6;
激光控制器1的激光控制信号输出端连接激光头2的激光控制信号输入端,光电探测器4均匀分布在火焰3周围,光电探测器4的光电探测信号输出端连接数据采集处理系统5的光电探测信号输入端,微透镜阵列光场相机6的断层三维辐射场信号输出端连接数据采集处理系统5的断层三维辐射场信号输入端;
所述方法包括以下步骤:
步骤一:开启激光控制器1,使激光头2发射出的脉冲激光入射到火焰3上,然后使激光头2顺时针旋转,且每次旋转的角度均为30°,直至激光头2旋转回初始位置,激光头2每旋转一次则利用光电探测器4采集一次火焰3透射或反射出的激光,然后将光电探测器4获得的所有光电探测信号发送至数据采集处理系统5中,数据采集处理系统5分别对其获得的光电探测信号进行处理,获得火焰3边界的出射光谱辐射强度值Imi,i为激光头2旋转次数,i=1,2,...,12;
步骤二:关闭激光控制器1,保持火焰3稳定,利用微透镜阵列光场相机6从单角度成像获取火焰3的断层三维辐射场信号,然后将该断层三维辐射场信号发送至数据采集处理系统5,数据采集处理系统5对断层三维辐射场信号进行处理,获得火焰3自身的出射光谱辐射强度测量值In;
步骤三:将步骤一获得的火焰3边界的出射光谱辐射强度值Imi分别与步骤二获得的火焰3自身的出射光谱辐射强度测量值In作差,获得消除了自辐射的光谱辐射强度测量值Mi;
步骤四:根据逆问题算法假设火焰光谱辐射物性参数为k,将该火焰光谱辐射物性参数k代入瞬态辐射传输方程,并对该瞬态辐射传输方程求解获得消除了自辐射的光谱辐射强度估计值Ei;
步骤五:以步骤三获得的光谱辐射强度测量值Mi与步骤四获得的光谱辐射强度估计值Ei的最小二乘作为火焰光谱辐射物性参数目标函数方程,并求解该火焰光谱辐射物性参数目标函数方程的目标函数值,然后将该目标函数值与火焰光谱辐射物性参数的给定目标函数阈值进行比较,判断目标函数值是否小于给定目标函数阈值,是则步骤四中假设的火焰光谱辐射物性参数k为真实的火焰光谱辐射物性参数,然后执行步骤六,否则返回步骤四;
步骤六:利用真实的火焰光谱辐射物性参数和逆问题算法假设火焰温度场为T,将假设火焰温度场T代入火焰3的温度场辐射传输方程,并对火焰3的温度场辐射传输方程求解获得火焰3自身的出射辐射强度估计值I'n;
步骤七:利用步骤二获得的出射光谱辐射强度测量值In和步骤六获得的出射辐射强度估计值I'n构造温度场目标函数方程,并求解该温度场目标函数方程的温度场目标函数值,然后将该温度场目标函数值与温度场给定目标函数阈值进行比较,判断该温度场目标函数值是否小于温度场给定目标函数阈值,是则步骤六中假设的火焰温度场T为火焰3的真实温度场,完成火焰温度泛尺度光场探测,否则返回步骤六。
本实施方式中,火焰3为高温火焰,其温度范围在1000~3000℃之间。
通过激光控制器1触发激光脉冲,照射到火焰3表面,激光在火焰3中被吸收和散射后会从火焰边界透射或反射出来,利用布置在火焰3周围的光电探测器4接收出射信号,信号经电缆进入数据采集处理系统5,最终获得火焰3边界不同方向的出射光谱辐射强度值Imi。
步骤三中,消除了自辐射的光谱辐射强度测量值Mi为:激光穿透火焰3后消除了火焰3自身辐射的光谱辐射强度值。
步骤四中,逆问题算法采用共轭梯度法。消除了自辐射的光谱辐射强度估计值Ei为:在激光照射的情况下火焰边界的出射光谱辐射强度估计值与无激光照射时火焰本身的出射光谱辐射强度估计值的差值,由于估计火焰边界的出射光谱辐射强度时也需要配合估计出激光头2旋转过程的所有值,因此本步骤中所估计的差值也为i个。
具体实施方式二:本实施方式是对具体实施方式一所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法作进一步说明,本实施方式中,步骤四中,假设的火焰光谱辐射物性参数k包括:火焰3的吸收系数κa和火焰3的散射系数κs,火焰3的衰减系数κe为火焰3的吸收系数κa和火焰3的散射系数κs之和。
具体实施方式三:本实施方式是对具体实施方式二所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法作进一步说明,本实施方式中,步骤四中利用离散坐标法对瞬态辐射传输方程求解,瞬态辐射传输方程的表达式如下:
其中,t为辐射时刻,r为辐射位置,Ω为辐射入射方向,Ω′为辐射散射方向,c为火焰中的光速,Φ(Ω′,Ω)为火焰3的散射相函数。
本实施方式中,根据上述参数的定义可知Ei(r,Ω,t)为在t时刻、辐射位置为r、辐射入射方向为Ω处的辐射强度。
具体实施方式四:本实施方式是对具体实施方式一所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法作进一步说明,本实施方式中,步骤五中火焰光谱辐射物性参数目标函数方程的表达式如下:
其中,Fobj表示火焰光谱辐射物性参数目标函数值,Nm表示估计的辐射强度值的个数。
具体实施方式五:本实施方式是对具体实施方式一所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法作进一步说明,本实施方式中,步骤六中火焰3的辐射传输方程的表达式如下:
其中,I(s,s)表示位置为s、方向为s处的辐射强度,Ib(s)表示黑体辐射强度,可表示为其中σ为黑体辐射常数,T为火焰温度场,κa表示火焰3的吸收系数,κs表示火焰3的散射系数,κe表示火焰3的衰减系数,Φ(si,s)表示散射相函数,S(s,s)表示广义源项。
本实施方式中,广义源项S(s,s)由有限体积法求得。
若要求解沿BA方向的辐射强度,则可以将辐射传输方程改写为沿BA方向的形式,如下所示:
则sBA表示辐射方向BA方向,表示BA方向的广义源项。
火焰3自身的出射辐射强度估计值I'n能够通过对辐射方程进行分段积分来求得,该出射辐射强度估计值I'n为任意探测方向上从边界到探测器的出射辐射强度,如下式所示:
其中,表示第m个网格的广义源项,△sm表示第m个网格的微元长度,m表示网格总数,l表示网格序号,△sl表示第l个网格的微元长度,表示第i个网格的广义源项。
具体实施方式六:本实施方式是对具体实施方式一所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法作进一步说明,本实施方式中,步骤七中温度场目标函数的表达式如下:
其中,Fobj(T)表示温度场目标函数值。
Claims (5)
1.基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,该方法是基于以下装置实现的,所属装置包括:激光控制器(1)、激光头(2)、光电探测器(4)、数据采集处理系统(5)和微透镜阵列光场相机(6);
激光控制器(1)的激光控制信号输出端连接激光头(2)的激光控制信号输入端,光电探测器(4)均匀分布在火焰(3)周围,光电探测器(4)的光电探测信号输出端连接数据采集处理系统(5)的光电探测信号输入端,微透镜阵列光场相机(6)的断层三维辐射场信号输出端连接数据采集处理系统(5)的断层三维辐射场信号输入端;
其特征在于,所述方法包括以下步骤:
步骤一:开启激光控制器(1),使激光头(2)发射出的脉冲激光入射到火焰(3)上,然后使激光头(2)顺时针旋转,且每次旋转的角度均为30°,直至激光头(2)旋转回初始位置,激光头(2)每旋转一次则利用光电探测器(4)采集一次火焰(3)透射或反射出的激光,然后将光电探测器(4)获得的所有光电探测信号发送至数据采集处理系统(5)中,数据采集处理系统(5)分别对其获得的光电探测信号进行处理,获得火焰(3)边界的出射光谱辐射强度值Imi,i为激光头(2)旋转次数,i=1,2,...,12;
步骤二:关闭激光控制器(1),保持火焰(3)稳定,利用微透镜阵列光场相机(6)从单角度成像获取火焰(3)的断层三维辐射场信号,然后将该断层三维辐射场信号发送至数据采集处理系统(5),数据采集处理系统(5)对断层三维辐射场信号进行处理,获得火焰(3)自身的出射光谱辐射强度测量值In;
步骤三:将步骤一获得的火焰(3)边界的出射光谱辐射强度值Imi分别与步骤二获得的火焰(3)自身的出射光谱辐射强度测量值In作差,获得消除了自辐射的光谱辐射强度测量值Mi;
步骤四:根据逆问题算法假设火焰光谱辐射物性参数为k,将该火焰光谱辐射物性参数k代入瞬态辐射传输方程,并对该瞬态辐射传输方程求解获得消除了自辐射的光谱辐射强度估计值Ei;
步骤五:以步骤三获得的光谱辐射强度测量值Mi与步骤四获得的光谱辐射强度估计值Ei的最小二乘作为火焰光谱辐射物性参数目标函数方程,并求解该火焰光谱辐射物性参数目标函数方程的目标函数值,然后将该目标函数值与火焰光谱辐射物性参数的给定目标函数阈值进行比较,判断目标函数值是否小于给定目标函数阈值,是则步骤四中假设的火焰光谱辐射物性参数k为真实的火焰光谱辐射物性参数,然后执行步骤六,否则返回步骤四;
步骤六:利用真实的火焰光谱辐射物性参数和逆问题算法假设火焰温度场为T,将假设火焰温度场T代入火焰(3)的温度场辐射传输方程,并对火焰(3)的温度场辐射传输方程求解获得火焰(3)自身的出射辐射强度估计值I'n;
步骤七:利用步骤二获得的出射光谱辐射强度测量值In和步骤六获得的出射辐射强度估计值I'n构造温度场目标函数方程,并求解该温度场目标函数方程的温度场目标函数值,然后将该温度场目标函数值与温度场给定目标函数阈值进行比较,判断该温度场目标函数值是否小于温度场给定目标函数阈值,是则步骤六中假设的火焰温度场T为火焰(3)的真实温度场,完成火焰温度泛尺度光场探测,否则返回步骤六。
2.根据权利要求1所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,其特征在于,步骤四中,假设的火焰光谱辐射物性参数k包括:火焰(3)的吸收系数κa和火焰(3)的散射系数κs,火焰(3)的衰减系数κe为火焰(3)的吸收系数κa和火焰(3)的散射系数κs之和。
3.根据权利要求2所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,其特征在于,步骤四中利用离散坐标法对瞬态辐射传输方程求解,瞬态辐射传输方程的表达式如下:
<mrow>
<mfrac>
<mn>1</mn>
<mi>c</mi>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>&Omega;</mi>
<mo>&CenterDot;</mo>
<mo>&dtri;</mo>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&kappa;</mi>
<mi>e</mi>
</msub>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<msub>
<mi>&kappa;</mi>
<mi>s</mi>
</msub>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<msub>
<mo>&Integral;</mo>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
</msub>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msup>
<mi>&Omega;</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>&Omega;</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>d&Omega;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
其中,t为辐射时刻,r为辐射位置,Ω为辐射入射方向,Ω′为辐射散射方向,c为火焰中的光速,Φ(Ω′,Ω)为火焰(3)的散射相函数。
4.根据权利要求1所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,其特征在于,步骤五中火焰光谱辐射物性参数目标函数方程的表达式如下:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>m</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>m</mi>
</msub>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
其中,Fobj表示火焰光谱辐射物性参数目标函数值,Nm表示估计的辐射强度值的个数。
5.根据权利要求1所述的基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法,其特征在于,步骤七中温度场目标函数的表达式如下:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>I</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
其中,Fobj(T)表示温度场目标函数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510889073.5A CN105547469B (zh) | 2015-12-04 | 2015-12-04 | 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510889073.5A CN105547469B (zh) | 2015-12-04 | 2015-12-04 | 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105547469A CN105547469A (zh) | 2016-05-04 |
CN105547469B true CN105547469B (zh) | 2018-01-30 |
Family
ID=55826820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510889073.5A Active CN105547469B (zh) | 2015-12-04 | 2015-12-04 | 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105547469B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106018286B (zh) * | 2016-05-18 | 2018-12-11 | 哈尔滨工业大学 | 基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法 |
CN106023082B (zh) * | 2016-05-24 | 2019-01-08 | 哈尔滨工业大学 | 基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置及其重建方法 |
CN108983273B (zh) * | 2018-08-22 | 2021-09-24 | 中国辐射防护研究院 | 一种人体内部三维剂量分布的实时测量系统及方法 |
CN109632718B (zh) * | 2019-01-25 | 2021-05-04 | 哈尔滨工业大学 | 基于时频光信息融合的弥散介质光学参数场测量装置及方法 |
CN109974893B (zh) * | 2019-03-27 | 2020-11-03 | 东南大学 | 一种梯度折射率火焰三维温度场测量方法 |
CN110445563B (zh) * | 2019-09-09 | 2021-12-24 | 上海无线电设备研究所 | 一种微波电磁信号的火焰衰减测试方法 |
CN111795746B (zh) * | 2020-06-09 | 2021-09-10 | 哈尔滨工业大学 | 基于主被动光学层析融合探测的火焰多参数场协同测量方法 |
CN114388073B (zh) * | 2022-01-14 | 2024-08-30 | 哈尔滨工业大学 | 基于参数水平集与调频激光的火焰多参数瞬态测量装置及方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8908087B2 (en) * | 2010-11-16 | 2014-12-09 | Nikon Corporation | Multiband camera, and multiband image capturing method |
CN102519391B (zh) * | 2011-12-28 | 2013-11-20 | 东南大学 | 基于弱饱和二维图像的物体表面三维立体图像重建方法 |
JP2014095688A (ja) * | 2012-10-09 | 2014-05-22 | Ricoh Co Ltd | 撮像装置及び撮像システム |
CN104819774A (zh) * | 2015-04-22 | 2015-08-05 | 哈尔滨工业大学 | 一种基于微透镜阵列的火焰光场探测泛尺度分析方法 |
CN104897632A (zh) * | 2015-06-01 | 2015-09-09 | 哈尔滨工业大学 | 一种基于扫描式平面激光诱导荧光成像系统测量瞬态燃烧场oh基浓度三维空间分布的方法 |
-
2015
- 2015-12-04 CN CN201510889073.5A patent/CN105547469B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105547469A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105547469B (zh) | 基于微透镜阵列与脉冲激光的火焰温度泛尺度光场探测方法 | |
Huang et al. | Simultaneous reconstruction of 3D temperature distribution and radiative properties of participating media based on the multi-spectral light-field imaging technique | |
CN105547485B (zh) | 基于微透镜阵列与调制激光的火焰温度泛尺度光场探测方法 | |
CN105571741B (zh) | 基于微透镜阵列与连续激光的火焰温度泛尺度光场探测方法 | |
CN106018286B (zh) | 基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法 | |
Aldave et al. | Review of thermal imaging systems in composite defect detection | |
CN103557886B (zh) | 一种基于超光谱的高温气体二维瞬态温度场、浓度场分布同时测量的方法 | |
Murray-Bruce et al. | Estimating localized sources of diffusion fields using spatiotemporal sensor measurements | |
Sun et al. | A novel calibration method of focused light field camera for 3-D reconstruction of flame temperature | |
CN111795746B (zh) | 基于主被动光学层析融合探测的火焰多参数场协同测量方法 | |
Niu et al. | Temperature field reconstruction of 3D luminous flames based on light field tomography theory | |
CN104729430B (zh) | 一种塔式太阳能热发电用定日镜面形检测方法 | |
Yu et al. | Rapid tomographic reconstruction based on machine learning for time-resolved combustion diagnostics | |
CN108362733B (zh) | 基于锁相热波与光学层析相结合的半透明材料光热特性分布测量方法 | |
Caliot et al. | Validation of a Monte Carlo integral formulation applied to solar facility simulations and use of sensitivities | |
Huang et al. | Application of Landweber method for three-dimensional temperature field reconstruction based on the light-field imaging technique | |
Liu et al. | Volumetric imaging of flame refractive index, density, and temperature using background-oriented Schlieren tomography | |
Shi et al. | Three-dimensional temperature reconstruction of diffusion flame from the light-field convolution imaging by the focused plenoptic camera | |
Liu et al. | Investigation and optimization of sampling characteristics of light field camera for flame temperature measurement | |
Niu et al. | Three-dimensional inhomogeneous temperature tomography of confined-space flame coupled with wall radiation effect by instantaneous light field | |
Mingyuan et al. | An efficient regulation approach for tomographic reconstruction in combustion diagnostics based on TDLAS method | |
Li et al. | Effect of nonuniform radiation properties on flame temperature reconstruction based on light field imaging | |
Dai et al. | 3-D soot temperature and volume fraction reconstruction of afterburner flame via deep learning algorithms | |
Li et al. | Joint method for reconstructing three-dimensional temperature of flame using Lucy-Richardson and nearest neighbor filtering using light-field imaging | |
Wei et al. | 3-D reconstruction algorithm of flame based on inversion calculation of thermal radiation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190617 Address after: 150000 Heilongjiang Harbin Dalian economic and Trade Zone, the North Road and Xingkai Road intersection Patentee after: HIT ROBOT GROUP Co.,Ltd. Address before: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin Patentee before: Harbin Institute of Technology |
|
TR01 | Transfer of patent right | ||
PP01 | Preservation of patent right |
Effective date of registration: 20240626 Granted publication date: 20180130 |
|
PP01 | Preservation of patent right |