CN103048050A - 一种测量目标真实温度的辐射测温方法和仪器 - Google Patents

一种测量目标真实温度的辐射测温方法和仪器 Download PDF

Info

Publication number
CN103048050A
CN103048050A CN2012104177586A CN201210417758A CN103048050A CN 103048050 A CN103048050 A CN 103048050A CN 2012104177586 A CN2012104177586 A CN 2012104177586A CN 201210417758 A CN201210417758 A CN 201210417758A CN 103048050 A CN103048050 A CN 103048050A
Authority
CN
China
Prior art keywords
temperature
wavelength
target
warm
approximate
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.)
Granted
Application number
CN2012104177586A
Other languages
English (en)
Other versions
CN103048050B (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.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
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 Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201210417758.6A priority Critical patent/CN103048050B/zh
Publication of CN103048050A publication Critical patent/CN103048050A/zh
Application granted granted Critical
Publication of CN103048050B publication Critical patent/CN103048050B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radiation Pyrometers (AREA)

Abstract

本发明公开了一种测量目标真实温度的辐射测温方法,包括:通过三波长辐射测温仪获取三个波长对应的实测亮温Tsi、目标的近似真温T以及三个比色温度Tci,三个波长为等间隔;根据实测亮温Ts1、Ts2、Ts3以及近似真温T的数值,确定目标在三个波长下的发射率ελ1、ελ2、ελ3;根据计算出的目标在三个波长下的发射率值及近似真温T值,计算目标对应三个波长的比色温度Tci’;根据实测比色温度与计算比色温度,求得均方差之和∑ΔT2;利用改变(Δ=1)后的近似真温的值进一步计算发射率以及比色温度值,进一步求得实测比色温度与该计算比色温度的均方差,比较两个均方差的大小,直至均方差之和最小时,即找到实际真温值T真。利用本发明的方法,可测得目标的真实温度。

Description

一种测量目标真实温度的辐射测温方法和仪器
技术领域
本发明涉及利用辐射进行测温的技术领域,特别涉及一种在不接触目标的情况下测量目标真实温度的方法和仪器。
背景技术
热辐射温度计是通过测定目标辐射强度,来测定目标温度的温度计。它的理论基础是关于物体辐射强度和吸收本领关系的基尔霍夫(kirchhoff)定律(在局部热平衡时目标的吸收本领等于其发射率),以及关于黑体辐射的普朗克(Planck)定律。辐射测温有很多优点(例如,它没有测量上限限制,不扰动被测目标的温度场和热平衡;它测量精度高,动态响应好等等);已成为非接触测温的主要方法。目前世界各国已有很多型号的辐射温度计产品,并得到广发应用。各国高温计量基准和标准温度计都是辐射温度计(而且都是亮度温度计)。但是辐射温度计的最大缺点是,它测得的实际目标的温度不是目标的真实温度(简称真温),而是目标的表观温度。目前表观温度主要包括亮度温度、辐射温度和比色温度(或称颜色温度,简称色温)三种。对应这三种表观温度可以引伸出三种基本的辐射测温方法和仪表;即亮度法测温仪表,全辐射测温仪表(简称辐射温度计)和比色法测温仪表。其中亮度法测温仪表是辐射测温中最重要、最普及的仪表。它在引入有效波长概念后,在理论上十分严格,已成了各国高温计量基准和标准仪表,并在计量量值传递和工业应用方面长期起主导作用。
热辐射温度计在检定(校准)中,都是以黑体炉为依据,是对黑体进行分度处理的。但在测量实际目标(非黑体)温度时,得到的是目标的表观温度(亮度温度、辐射温度或比色温度)。要知道目标的真实温度,则要知道目标的发射率ε(λ.T)。目标的发射率不但是目标的温度和波长的函数,而且还与目标的材料性质、表面状态等因素有关。这样就对确定目标的真温带来较大困难和不确定性。近百年来,辐射测温的这一重大缺点一直困扰着测温界。20世纪80年代以后,多国科学家对多波长温度计进行了大量研究,企图通过应用多波长温度计等方法测得实际目标的真温。从三波长温度计开始进一步研制六波长温度计,甚至35个波长的温度计[1]等均研制出来。但到目前为止,这些多波长温度计都不能测定目标的真实温度。本专利提出一种特定三波长辐射温度计,建立求真温数学模型,并通过高速迭代运算方法,实时测定目标的真实温度的方法和仪器。
发明内容
本发明为了克服现有技术存在的上述缺陷,即针对不能测定目标真实温度的问题,提出一种特定三波长辐射测温仪器以及利用该仪器进行测温的方法,建立求真温数学模型,并通过高速迭代运算方法,实时测定目标的真实温度。
本发明提出的测量目标真实温度的辐射测温方法包括:步骤1,通过三波长辐射测温仪获取三个波长对应的实测亮温Tsi、目标的近似真温T以及三个比色温度Tci,其中i的值为1、2和3,三个波长为等间隔;步骤2,根据实测亮温Ts1、Ts2、Ts3以及近似真温T的数值,确定目标在三个波长下的发射率ελ1、ελ2、ελ3;步骤3,根据步骤2计算出的目标在三个波长下的发射率值及近似真温T值,计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’;步骤4,根据步骤1中的实测比色温度与步骤3中计算出的比色温度,求得均方差之和∑ΔT2如下:∑ΔT2=(Tc1-Tc1’)2+(Tc2-Tc2’)2+(Tc3-Tc3’)2;步骤5,令近似真温T’=T+1,并根据实测亮温Ts1、Ts2、Ts3按照步骤2的方式计算出目标的三个波长的发射率ελ1′、ελ2′、ελ3′;步骤6,按照步骤5计算得到的发射率ελ1′、ελ2′、ελ3′及T’=T+1的值,按照步骤3的方式计算出目标的三个比色温度值Tc1”、Tc2”、Tc3”;步骤7,根据计算比色温度Tc1”、Tc2”、Tc3”值与实测比色温度Tc1、Tc2、Tc3,求得均方差之和为:(∑ΔT2)’=(Tc1-Tc1”)2+(Tc2-Tc2”)2+(Tc3-Tc3”)2;步骤8,若(∑ΔT2)’<(∑ΔT2),则近似真温T’=T+1比T更接近实际真温T,则令近似真温T”=T’+1=T+2,继续重复步骤5-7,直至均方差之和最小时,即认为找出实际真温T值,若(∑ΔT2)’>(∑ΔT2),则近似真温T’=T+1已远离实际真温T,向反方向求实际真温值,即令T’=T-1,继续按步骤5-7算出∑ΔT2,直至均方差之和最小时,即找到实际真温值T
本发明还提供一种测量目标真实温度的辐射测温仪器,该测温仪器为三波长测温仪器,其包括:温度获取单元,用于获取三个波长对应的实测亮温Tsi、目标的近似真温T以及三个比色温度Tci,其中i的值为1、2和3,三个波长为等间隔;高速计算单元,其进一步包括:第一单元,用于根据获取的实测亮温Ts1、Ts2、Ts3以及近似真温T的数值,确定目标在三个波长下的发射率ελ1、ελ2、ελ3;第二单元,根据第一单元计算出的目标在三个波长下的发射率值及近似真温T值,计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’;第三单元,根据实测比色温度与第二单元求得的比色温度,求得均方差之和∑ΔT2如下:∑ΔT2=(Tc1-Tc1’)2+(Tc2-Tc2’)2+(Tc3-Tc3’)2;第四单元,其令近似真温T’=T+1,并根据实测亮温Ts1、Ts2、Ts3按照第一单元的方式计算出目标的三个波长的发射率ελ1′、ελ2′、ελ3′;第五单元,按照第四单元计算得到的发射率ελ1′、ελ2′、ελ3′及T’=T+1的值,按照第二单元的方式计算出目标的比色温度值Tc1”、Tc2”、Tc3”;第六单元,根据计算比色温度Tc1”、Tc2”、Tc3”值与实测比色温度Tc1、Tc2、Tc3,求得均方差之和为:(∑ΔT2)’=(Tc1-Tc1”)2+(Tc2-Tc2”)2+(Tc3-Tc3”)2;第七单元,若(∑ΔT2)’<(∑ΔT2),则近似真温T’=T+1比T更接近实际真温T,则令近似真温T”=T’+1=T+2,继续重复第四-六单元的处理过程,直至均方差之和最小时,即认为找出实际真温T值,若(∑ΔT2)’>(∑ΔT2),则近似真温T’=T+1已远离实际真温T,向反方向求实际真温值,即令T’=T-1,继续重复第四-六单元的处理过程算出∑ΔT2,直至均方差之和最小时,即找到实际真温值T
附图说明
图1为本发明测量目标真实温度的辐射测温仪器结构图;
图2为本发明测量目标真实温度的辐射测温方法的软件流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
图1为本发明测量目标真实温度的辐射测温仪器结构图。
本发明的测温仪器为特定波长的三波长(实为三个窄波段)温度计,三个特定波长按照以下原则进行选取:(1)三个波长的间隔为等距,且间隔不宜太大。间隔太大影响本专利唯一假定的正确性(见后)。如在近红外区,以0.1~0.3μm(微米)左右为宜。例如取λ1=0.8μm,λ2=0.9μm,λ3=1.0μm三个波长来测温。当然间隔也不宜太小,因间隔太小时波长间差异小,影响测量精度。(2)为了减少中间气体吸收对测量影响,上述三波长应避开空气中常见的水蒸气、CO2等气体分子吸收带。(3)为了提高测量精度,上述三波长应选择在常用光电元件光谱灵敏曲线中灵敏度较高的区域。例如上述选取0.8~1.0之间三个波长,即选取常用硅光电管灵敏度较高的区域。(4)为了数据处理方便,目标光辐射强度变化范围(对应为目标温度范围),应选择在光电元件线性区工作范围内。因此温度计的光电元件输出的光电流讯号比例于该波长时目标的辐射强度。(5)这台三波长温度计先用黑体炉分度;可分别得到三条亮度温度分度曲线(分别对应波长λ1、λ2、λ3三个波长),及三条比色温度分度曲线(分别对应波长λ1比λ2,波长λ1比λ3,λ2比λ3);以及近似真温分度曲线R=f(T)(详见后(5)及(6)式)。当用这台三波长温度计测量目标温度时,可同时得到目标的三个亮度温度,三个比色温度及一个近似真温T,共7个参数。然后通过本发明提供的数学模型,按图2所示软件运算流程,迅速计算出目标的真实温度T。其中原则(1)及(5)最为重要,它反应了本发明特点,其他三个原则与一般温度计相同。。
参照图1,本发明测量目标真实温度的辐射测温仪器(下面简称测温仪)包括物镜(可调焦)1、带中心小孔的反射镜2、瞄准反射镜3、目标成像小透镜4、由电机M带动的调制盘5、测量光电管6、发光二极管LED7、控制用同步光电管8、观察用目镜9。
测温仪的物镜1可调焦,将目标成像于带中心小孔的反射镜2上。由2反射的光经反射镜3再反射后,经目镜9成像供人眼观察瞄准用。反射镜2中心小孔的光透过反射镜2后,由目标成像小透镜4聚焦成像在测量光电管6上。因此中心小孔的直径决定测量点的直径;物镜1的焦距f除以中心小孔直径d,即f/d的数值为测温仪的重要指标,称为测温仪的距离系数。测温用的光线经过4后,需通过由调制电机M带动的调制盘。调制盘上装有三个窄带滤光片,分别让λ1、λ2、λ3三个波长的光分时通过调制盘,照射在测量光电管6上。调制盘上有一个小孔,当安装在调制盘左侧的发光二极管LED的光透过小孔,短时间照射在同步用光电管8上。在调制盘每转动一周时产生一个同步脉冲信号,同步信号用于信号处理时区分三个测量波长的信号。测量光电管6的信号经前置放大器及A/D转换后,与同步信号共同输入到存储器。
因存储器内事先(在用黑体炉分度测温仪时)已存储三个波长的三条亮度温度及三条比色温度的分度曲线及R=f(T)曲线(见后说明);故通过存储器可迅速查找出对应三个波长的亮温Ts1、Ts2、Ts3及比色温度Tc1、Tc2、Tc3及近似真温T的数值。将Ts1、Ts2、Ts3、Tc1、Tc2、Tc3及T的数值输入到高速计算单元,按图2所示迭代运算软件流程,可快速计算出目标的真实温度T的数值。采用现代快速计算技术,可以在一百毫秒内完成迭代运算,可以认为是实时测量真实温度。图1所示的测温仪的结构中,选通三个波长是通过调制盘5转动实现的。实际上选通三个波长还可以通过棱镜分光,以及析光镜等方法实现。实现选通三个波长的方法很多,它不影响本发明的实现与效果。
当用黑体炉标定测温仪时,对应某波长的光电元件输出的光电流值Iλ为:
I λ = k λ M ( λ . T ) = k λ C 1 λ - 5 [ e C 2 / λT - 1 ] - 1 - - - ( 1 )
式中:Iλ为对应波长λ测温仪的光电元件的光电流值;kλ为对应波长λ测温仪的光电元件光电转换系数(波长确定后为常数);M(λ.T)为普朗克公式,式中C1、C2分别为第一辐射常数和第二辐射常数;T为黑体炉温度;λ为波长。在大多数实用场合C2>>λT,用维恩公式可以足够精确地替代普朗克公式,则上述(1)式简化为:
I λ = k λ C 1 λ - 5 e - C 2 / λT - - - ( 1 ) ,
当用亮度温度计测量实际目标(非黑体)的温度时,测得的是目标的亮度温度Ts(简称亮温Ts)。而目标的真实温度(简称真温)T和它的亮温Ts有以下关系(由亮温定义得出):
M(λ.Ts)=ελM(λ.T)                  (2)
以普朗克公式代入上式,有(2)式的另一表达形式:
Ts = C 2 λ ln [ e c 2 / λ i T - 1 ϵ λ + 1 ] - - - ( 2 )
式中:T为目标的真实温度;Ts为目标的亮度温度;ελ为目标在波长λ时的发射率(在0~1之间)。由维恩公式代入(2)式,可得出到亮温和真温关系简化式为
Ts = C 2 T C 2 - λ T ln ϵλ - - - ( 2 ) ,
因为是三波长测温,因此(2)及(2)’式可以有三个波长λ1、λ2、λ3,对应可以用Ts1、Ts2、Ts3表示三个波长的亮度温度。
同理按比色温度定义,可分别测出三个比色温度,用维恩公式代入后可得到三个比色温度(详可参见李吉林、肖功弼。俞伦鹏编著,辐射测温和检定/校准技术,中国计量出版社,2009年6月,P85):
Figure BDA00002315200600081
上式中,TC1代表波长λ1对波长λ2的比色温度;TC2代表波长
λ1对波长λ3的比色温度;TC3代表波长λ2对波长λ3的比色温度。其他参数意义同(2)式。
对于绝大多数固体或液体(含金属或非金属以及金属氧化物等),它们的发射率随波长的变化都较平缓(可参见李吉林.肖功弼.俞伦鹏编著;《辐射测温和检定/校准技术》;中国计量出版社.2009年8月出版;P36-P37),也就是说绝大部分固体或液体的发射率,随波长的变化在小范围内是近似线性变化。辐射温度计大多用于测量固体或液体温度,很少用于测量气体的温度。
我们假定测量目标(非黑体)的发射率随波长变化,在λ1至λ3小范围内是线性变化。这是本发明推导的唯一假定,而这一假定对绝大多数固体和液体都是近似正确。根据这一近似正确的假定,可以得出当λ1、λ2、λ3的间隔为等距时有:
ελ1×ελ3≈ε2 λ2               (4)
式中:ελ1为波长λ1对应目标的发射率;ελ2为波长λ2对应目标的发射率;ελ3为波长λ3对应目标的发射率。
因为是三波长测温仪,它可以同时测出目标对应三个波长的光电流值,对应为
Figure BDA00002315200600082
Figure BDA00002315200600083
Figure BDA00002315200600084
它们可以有以下关系(参考(1)式及(2)式可得出:
I λ 1 × I λ 3 I 2 λ 2 = ϵ λ 1 · k λ 1 · M ( λ 1 . T ) × ϵ λ 3 · k 3 · M ( λ 3 . T ) [ ϵ λ 2 · k λ 2 · M ( λ 2 . T ) ] 2 - - - ( 5 )
式中:ελ1、ελ2、ελ3分别为对应目标对应三个波长的发射率;
Figure BDA00002315200600092
Figure BDA00002315200600093
分别为温度计的光电元件对应三个波长的光电转换系数(对应每个波长均为常数);M(λ1.T)、M(λ2.T)、M(λ3.T)分别为普朗克公式即温度为T时对应三个波长的辐射强度(T为目标的真温)。将近似式(4)代入(5)式,并令
Figure BDA00002315200600095
可得:
R = K λ 1 · M ( λ 1 . T ) × K λ 3 · M ( λ 3 . T ) K λ 2 2 [ M ( λ 2 . T ) ] 2 ≈ f ( T ) - - - ( 6 )
通过(6)式可以看出R值已与目标的光谱发射率ελ无关,它只是目标真实温度T的函数。而R值完全可以通过三波长温度计对黑体炉进行检定实际测出并存储在存储器中。而在用黑体炉检定时,(4)式及(6)式精确相等(因黑体的发射率为1)。在导出(6)式时,采用了本发明的唯一近似正确的假定,即目标的发射率在λ1至λ3小范围内是线性变化,并由此导出近似式(4)。实际上目标发射率在λ1至λ3小范围内不一定呈线性变化,而且(4)式是近似正确。这样按(6)测出的目标真温会有一定误差。因此我们称按(6)式测出的目标真温为目标的近似真温T。
令目标的实际真温用T表示。近似真温T和实际真温T会有一定误差。但近似真温T肯定在T的附近(即两者间差异不会太大,在常用的工业用温度范围内,一般在100度左右)。这样可按下述方法(数学模型),以近似真温T为起点,进行迭代运算寻找实际真温。
图2为本发明测量目标真实温度的辐射测温方法流程图,图中主要体现的是根据目标近似真温T值,寻找目标实际真温T的迭代运算过程,减退包括:
步骤1,根据实际测出的目标三个亮温Ts1、Ts2、Ts3数值以及按(6)式求出的目标近似真温T值,根据(2)式可算出目标在三个波长下的发射率如下:
ϵ λi = M ( λ i . T si ) M ( λ i . T ) = I λi I ( λ i . T ) - - - ( 7 )
式中:
Figure BDA00002315200600102
分别为波长λ1、λ2、λ3对应的目标发射率
Figure BDA00002315200600104
分别为波长λ1、λ2、λ3对应的目标(非黑体)光谱辐射强度所产生的光电流
Figure BDA00002315200600107
Figure BDA00002315200600109
;I(λi.T)分别为波长λ1、λ2、λ3三条标定曲线,近似真温T对应的光电流。(可以从标定曲线内存表中查出)。即温度为T时黑体辐射强度所产生的光电流I(λ1.T)、I(λ2.T)、I(λ3.T)。
这样已经初步确定了目标对应三个波长下的发射率数值(当然还有一定误差,因为近似真温有误差)。
步骤2,根据(7)式算出的目标在三个波长下的近似发射率值及近似真温T值,可按(3)式算出目标对应三个波长的比色温度值Tc1’、Tc2’、Tc3’如下(在Tc上方加’表示此比色温度是初次计算出的):
Figure BDA000023152006001010
式中:T为近似真温,已按(6)式实际测出
C2为普朗克第二辐射常数,C2=1.4388×104μm.K
λi为波长(共有λ1、λ2、λ3三个波长)
为波长λi时的目标发射率。已按(7)式初步算出近似值。
步骤3,根据(8)式算出的目标三个比色温度值和实际测出的目标的三个比色温度值的差为零,则说明本发明唯一假定及相关计算完全正确。即近似真温T等于实际真温T。实际上总有一些误差(含假定正确性误差,(4)式近似性,实际测量误差以及计算误差等)。可按下式求出三个波长比色温度均方差之和∑ΔT2如下:(即最小二乘法)
∑ΔT2=(Tc1-Tc1’)2+(Tc2-Tc2’)2+(Tc3-Tc3’)2          (9)
式中:Tc1、Tc2、Tc3分别为实测目标的三个比色温度值
Tc1’、Tc2’、Tc3’分别为按(8)式算出的目标的三个比色温度值显然∑ΔT2越小说明近似真温T和实际真温T之间误差越小,只要找到∑ΔT2最小值处对应的近似真温值T,即找到了实际真温值T
步骤4,在近似真温T值附近,寻找实际真温T过程如下:
步骤401,令近似真温T’=T+1,及实际测出的目标的三个亮温Ts1、Ts2、Ts3按(7)式算出目标三个波长的发射率如下,
ϵ ′ λi = M ( λ i . T si ) M ( λ ′ 1 . T ′ ) = I λi I ( λ i . T ′ ) - - - ( 7 ) ,
式中符号含义与(7)式相同,计算的发射率ε上加’表示是进一步近似推算出的发射率,以区别首次按(7)式算出的发射率。
步骤402,根据(7)’式算出的
Figure BDA00002315200600113
Figure BDA00002315200600114
Figure BDA00002315200600115
值及T’=T+1值,可按(3)算出目标的三个比色温度值如下:
Figure BDA00002315200600121
在Tci上加”表示第二次算出的色温,其他参数意义同(8)式。
步骤403,根据(8)’式算出的Tc1”、Tc2”、Tc3”值与实际测出的Tc1、Tc2、Tc3,按(9)式求出均方差之和为
(∑ΔT2)’=(Tc1-Tc1”)2+(Tc2-Tc2”)2+(Tc3-Tc3”)2       (9)’
若(∑ΔT2)’<(∑ΔT2),则说明近似真温T’=T+1比T更接近实际真温T。则完全可令近似真温T”=T’+1=T+2,继续按(7)’、(8)’、(9)’式算下去,直至均方差之和最小时,即认为以寻求出实际真温T值。
若(∑ΔT2)’>(∑ΔT2),则说明近似真温T’=T+1已远离实际真温T,应反方向寻求实际真温值。即令T’=T-1,继续按(7)’、(8)’、(9)’式算出∑ΔT2,直至均方差之和最小时,即找到实际真温值T。以上迭代运算软件流程方块图如图2所示。
另外,本发明提出的测量目标真实温度的辐射测温仪器用于执行上述测温方法,该测温仪器为三波长测温仪器,其包括:温度获取单元,用于获取三个波长对应的实测亮温Tsi、目标的近似真温T以及三个比色温度Tci,其中i的值为1、2和3,三个波长为等间隔;;高速计算单元,其进一步包括:第一单元,用于根据获取的实测亮温Ts1、Ts2、Ts3以及近似真温T的数值,确定目标在三个波长下的发射率ελ1、ελ2、ελ3;第二单元,根据第一单元计算出的目标在三个波长下的发射率值及近似真温T值,计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’;第三单元,根据实测比色温度与第二单元求得的计算比色温度,求得均方差之和∑ΔT2如下:∑ΔT2=(Tc1-Tc1’)2+(Tc2-Tc2’)2+(Tc3-Tc3’)2;第四单元,其令近似真温T’=T+1,并根据实测亮温Ts1、Ts2、Ts3按照第一单元的方式计算出目标的三个波长的发射率ελ1′、ελ2′、ελ3′;
第五单元,按照第四单元计算得到的发射率ελ1′、ελ2′、ελ3′及T’=T+1的值,按照第二单元的方式计算出目标的三个比色温度值Tc1”、Tc2”、Tc3”;第六单元,根据计算比色温度Tc1”、Tc2”、Tc3”值与实测比色温度Tc1、Tc2、Tc3,求得均方差之和为:(∑ΔT2)’=(Tc1-Tc1”)2+(Tc2-Tc2”)2+(Tc3-Tc3”)2;第七单元,若(∑ΔT2)’<(∑ΔT2),则近似真温T’=T+1比T更接近实际真温T,则令近似真温T”=T’+1=T+2,继续重复第四-六单元的处理过程,直至均方差之和最小时,即认为找出实际真温T值,若(∑ΔT2)’>(∑ΔT2),则近似真温T’=T+1已远离实际真温T,向反方向求实际真温值,即令T’=T-1,继续重复第四-六单元的处理过程算出∑ΔT2,直至均方差之和最小时,即找到实际真温值T
以上迭代运算采用现代快速计算技术,在一百毫秒以内可以完成,可以认为是实时测量真实温度。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种测量目标真实温度的辐射测温方法,该方法包括:
步骤1,通过三波长辐射测温仪获取三个波长对应的实测亮温Tsi、目标的近似真温T以及三个比色温度Tci,其中i的值为1、2和3,三个波长为等间隔;
步骤2,根据实测亮温Ts1、Ts2、Ts3以及近似真温T的数值,确定目标在三个波长下的发射率ελ1、ελ2、ελ3
步骤3,根据步骤2计算出的目标在三个波长下的发射率值及近似真温T值,计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’;
步骤4,根据步骤1中的实测比色温度与步骤3中计算出的比色温度,求得均方差之和∑ΔT2如下:∑ΔT2=(Tc1-Tc1’)2+(Tc2-Tc2’)2+(Tc3-Tc3’)2
步骤5,令近似真温T’=T+1,并根据实测亮温Ts1、Ts2、Ts3按照步骤2的方式计算出目标的三个波长的发射率ελ1′、ελ2′、ελ3′;
步骤6,按照步骤5计算得到的发射率ελ1′、ελ2′、ελ3′及T’=T+1的值,按照步骤3的方式计算出目标的三个比色温度值Tc1”、Tc2”、Tc3”;
步骤7,根据计算的比色温度Tc1”、Tc2”、Tc3”值与实测比色温度Tc1、Tc2、Tc3,求得均方差之和为:(∑ΔT2)’=(Tc1-Tc1”)2+(Tc2-Tc2”)2+(Tc3-Tc3”)2
步骤8,若(∑ΔT2)’<(∑ΔT2),则近似真温T’=T+1比T更接近实际真温T,则令近似真温T”=T’+1=T+2,继续重复步骤5-7,直至均方差之和最小时,即认为找出实际真温T值,
若(∑ΔT2)’>(∑ΔT2),则近似真温T’=T+1已远离实际真温T,向反方向求实际真温值,即令T’=T-1,继续按步骤5-7算出∑ΔT2,直至均方差之和最小时,即找到实际真温值T
2.根据权利要求1所述的方法,其特征在于,在近红外区,三个波长之间的间隔范围在0.1~0.3μm。
3.根据权利要求1所述的方法,其特征在于,在步骤2,发射率ελ1、ελ2、ελ3通过下式求得:
Figure FDA00002315200500021
其中,
Figure FDA00002315200500022
分别为波长λ1、λ2、λ3对应的目标发射率ελ1、ελ2、ελ3;Iλi分别为波长λ1、λ2、λ3对应的目标光谱辐射强度所产生的光电流
Figure FDA00002315200500023
Figure FDA00002315200500024
Figure FDA00002315200500025
;I(λi.T)分别为波长λ1、λ2、λ3三条标定曲线,近似真温T对应的光电流。
4.根据权利要求3所述的方法,其特征在于,在步骤3,通过下式来计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’:
Figure FDA00002315200500026
其中:T为近似真温,C2为普朗克第二辐射常数,λi为波长,
Figure FDA00002315200500027
为波长λi时的目标发射率。
5.一种测量目标真实温度的辐射测温仪器,该测温仪器为三波长测温仪器,其包括:
温度获取单元,用于获取三个波长对应的实测亮温Tsi、目标的近似真温T以及三个比色温度Tci,其中i的值为1、2和3,三个波长为等间隔;
高速计算单元,其进一步包括:
第一单元,用于根据获取的实测亮温Ts1、Ts2、Ts3以及近似真温T的数值,确定目标在三个波长下的发射率ελ1、ελ2、ελ3
第二单元,根据第一单元计算出的目标在三个波长下的发射率值及近似真温T值,计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’;
第三单元,根据实测比色温度与第二单元求得的比色温度,求得均方差之和∑ΔT2如下:
∑ΔT2=(Tc1-Tc1’)2+(Tc2-Tc2’)2+(Tc3-Tc3’)2
第四单元,其令近似真温T’=T+1,并根据实测亮温Ts1、Ts2、Ts3按照第一单元的方式计算出目标的三个波长的发射率ελ1′、ελ2′、ελ3′;
第五单元,按照第四单元计算得到的发射率ελ1′、ελ2′、ελ3′及T’=T+1的值,按照第二单元的方式计算出目标的比色温度值Tc1”、Tc2”、Tc3”;
第六单元,根据计算比色温度Tc1”、Tc2”、Tc3”值与实测比色温度Tc1、Tc2、Tc3,求得均方差之和为:(∑ΔT2)’=(Tc1-Tc1”)2+(Tc2-Tc2”)2+(Tc3-Tc3”)2
第七单元,若(∑ΔT2)’<(∑ΔT2),则近似真温T’=T+1比T更接近实际真温T,则令近似真温T”=T’+1=T+2,继续重复第四-六单元的处理过程,直至均方差之和最小时,即认为找出实际真温T值,若(∑ΔT2)’>(∑ΔT2),则近似真温T’=T+1已远离实际真温T,向反方向求实际真温值,即令T’=T-1,继续重复第四-六单元的处理过程算出∑ΔT2,直至均方差之和最小时,即找到实际真温值T。6、根据权利要求5所述的测温仪器,其特征在于,在近红外区,三个波长之间的间隔范围在0.1~0.3μm。
6.根据权利要求6所述的测温仪器,其特征在于,在第一单元中,发射率ελ1、ελ2、ελ3通过下式求得:
Figure FDA00002315200500031
其中,
Figure FDA00002315200500041
分别为波长λ1、λ2、λ3对应的目标发射率ελ1、ελ2、ελ3
Figure FDA00002315200500042
分别为波长λ1、λ2、λ3对应的目标光谱辐射强度所产生的光电流
Figure FDA00002315200500043
Figure FDA00002315200500044
;I(λi.T)分别为波长λ1、λ2、λ3三条标定曲线,近似真温T对应的光电流。
7.根据权利要求6所述的测温仪器,其特征在于,
第二计算单元通过下式来计算目标对应三个波长的比色温度Tc1’、Tc2’、Tc3’:
Figure FDA00002315200500046
其中:T为近似真温,C2为普朗克第二辐射常数,λi为波长,
Figure FDA00002315200500047
为波长λi时的目标发射率。
8.根据权利要求5-7任一项所述的测温仪器,其特征在于,该仪器进一步包括物镜、带中心小孔的反射镜、瞄准反射镜、目标成像小透镜、由电机M带动的调制盘、测量光电管、发光二极管LED、控制用同步光电管、观察用目镜。
9.根据权利要求8所述的测温仪器,其特征在于,其中测温仪的物镜可调焦,将目标成像于带中心小孔的反射镜上,由带中心小孔的反射镜反射的光经瞄准反射镜再反射后,经目镜成像供人眼观察瞄准用,带中心小孔的反射镜的中心小孔的光透过该反射镜后,经过目标成像小透镜后,通过由调制电机M带动的调制盘,调制盘上装有三个窄带滤光片,分别让λ1、λ2、λ3三个波长的光分时通过调制盘,照射在测量光电管上,调制盘上有一个小孔,当安装在调制盘左侧的发光二极管LED的光透过小孔时,短时间照射在同步用光电管上,在调制盘每转动一周时产生一个同步脉冲信号,同步脉冲信号用于信号处理时区分三个测量波长的信号,测量光电管的信号经前置放大器及A/D转换后,与同步脉冲信号共同输入到存储器。
CN201210417758.6A 2012-07-30 2012-10-26 一种测量目标真实温度的辐射测温方法和仪器 Active CN103048050B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210417758.6A CN103048050B (zh) 2012-07-30 2012-10-26 一种测量目标真实温度的辐射测温方法和仪器

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201210266196XA CN102749141A (zh) 2012-07-30 2012-07-30 一种测量目标真实温度的辐射测温方法和仪器
CN201210266196.X 2012-07-30
CN201210417758.6A CN103048050B (zh) 2012-07-30 2012-10-26 一种测量目标真实温度的辐射测温方法和仪器

Publications (2)

Publication Number Publication Date
CN103048050A true CN103048050A (zh) 2013-04-17
CN103048050B CN103048050B (zh) 2014-09-24

Family

ID=47029514

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201210266196XA Pending CN102749141A (zh) 2012-07-30 2012-07-30 一种测量目标真实温度的辐射测温方法和仪器
CN201210417758.6A Active CN103048050B (zh) 2012-07-30 2012-10-26 一种测量目标真实温度的辐射测温方法和仪器

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201210266196XA Pending CN102749141A (zh) 2012-07-30 2012-07-30 一种测量目标真实温度的辐射测温方法和仪器

Country Status (1)

Country Link
CN (2) CN102749141A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103411940A (zh) * 2013-08-08 2013-11-27 哈尔滨工业大学 基于发射光谱的防热材料催化特性检测方法与测试装置
CN105181169A (zh) * 2015-05-12 2015-12-23 中国科学院苏州生物医学工程技术研究所 温度测量方法、温度测量系统和温度获取装置
CN105509901A (zh) * 2015-12-29 2016-04-20 温州大学 发光体相对色温的高精度计量方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105333962B (zh) * 2014-06-18 2018-06-22 中微半导体设备(上海)有限公司 一种修正双波段测温误差的温度测量方法及系统
CN106066340B (zh) * 2016-05-24 2018-11-23 南京航空航天大学 一种在常温环境中测量固体目标红外特征的模拟实验方法
CN106768366A (zh) * 2017-01-12 2017-05-31 杭州电子科技大学 一种利用微型光谱仪测量高温的方法及装置
CN111044565A (zh) * 2019-12-30 2020-04-21 西安近代化学研究所 一种燃烧转爆轰过程真温多光谱测量方法
CN111562019B (zh) * 2020-04-08 2021-04-06 太原理工大学 一种多光谱辐射测温方法及系统
CN112220460A (zh) * 2020-09-25 2021-01-15 电子科技大学 一种基于微波辐射技术的体内温度测量贴片
CN113175998B (zh) * 2021-03-26 2022-09-09 合肥工业大学 一种基于比色测温的金属材料表面温度测量方法
CN114112061B (zh) * 2021-10-11 2024-02-02 黑龙江业博科技有限公司 一种晶体管式焊接电源焊头温度快速在线测量方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0294747A2 (en) * 1987-06-06 1988-12-14 Minolta Camera Kabushiki Kaisha Pyrometer
US6247918B1 (en) * 1998-12-16 2001-06-19 Forney Corporation Flame monitoring methods and apparatus
WO2004036149A2 (en) * 2002-09-19 2004-04-29 Bechtel Bwxt Idaho, Llc Optical steam quality measurement system and method
RU2003121076A (ru) * 2003-07-08 2005-01-10 Открытое акционерное общество "КАМАЗ" (RU) Способ измерения цветовой температуры металла
US20060076495A1 (en) * 2004-09-09 2006-04-13 Volkswagen Ag Device for determining the surface condition of a roadway
CN201867924U (zh) * 2010-11-15 2011-06-15 上海翼捷工业安防技术有限公司 三波长红外火焰探测器
JP2012013702A (ja) * 2010-06-30 2012-01-19 General Electric Co <Ge> 多次元温度データを生成するためのマルチスペクトルシステム及び方法
CN102426753A (zh) * 2011-11-21 2012-04-25 上海大学 红外多参数火气探测器

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0294747A2 (en) * 1987-06-06 1988-12-14 Minolta Camera Kabushiki Kaisha Pyrometer
US6247918B1 (en) * 1998-12-16 2001-06-19 Forney Corporation Flame monitoring methods and apparatus
WO2004036149A2 (en) * 2002-09-19 2004-04-29 Bechtel Bwxt Idaho, Llc Optical steam quality measurement system and method
RU2003121076A (ru) * 2003-07-08 2005-01-10 Открытое акционерное общество "КАМАЗ" (RU) Способ измерения цветовой температуры металла
US20060076495A1 (en) * 2004-09-09 2006-04-13 Volkswagen Ag Device for determining the surface condition of a roadway
JP2012013702A (ja) * 2010-06-30 2012-01-19 General Electric Co <Ge> 多次元温度データを生成するためのマルチスペクトルシステム及び方法
CN201867924U (zh) * 2010-11-15 2011-06-15 上海翼捷工业安防技术有限公司 三波长红外火焰探测器
CN102426753A (zh) * 2011-11-21 2012-04-25 上海大学 红外多参数火气探测器

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
D.R.SNELLING,ET AL: "Spectrally Resolved Measurement of Flame Radiation to Determine Soot Temperature and Concentration", 《AIAA JOURNAL》, vol. 40, no. 9, 30 September 2002 (2002-09-30), pages 1789 - 1795 *
ZHANG XIANGYU,ET AL: "An improved colorimetric method for visualization of 2-D, inhomogeneous temperature distribution in a gas fired industrial furnace by radiation image processing", 《PROCEEDINGS OF THE COMBUSTION INSTITUTE》, vol. 33, no. 2, 15 September 2010 (2010-09-15), pages 2755 - 2762, XP028164884, DOI: doi:10.1016/j.proci.2010.06.119 *
徐宝昌,等: "一种改进的比色测温方法研究", 《光电工程》, vol. 38, no. 4, 30 April 2011 (2011-04-30), pages 1 - 6 *
杨词银,等: "基于实时标校的目标红外辐射测量新方法", 《红外与毫米波学报》, vol. 30, no. 3, 30 June 2011 (2011-06-30), pages 284 - 288 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103411940A (zh) * 2013-08-08 2013-11-27 哈尔滨工业大学 基于发射光谱的防热材料催化特性检测方法与测试装置
CN103411940B (zh) * 2013-08-08 2014-06-18 哈尔滨工业大学 基于发射光谱的防热材料催化特性检测方法与测试装置
CN105181169A (zh) * 2015-05-12 2015-12-23 中国科学院苏州生物医学工程技术研究所 温度测量方法、温度测量系统和温度获取装置
CN105181169B (zh) * 2015-05-12 2019-01-15 中国科学院苏州生物医学工程技术研究所 温度测量方法、温度测量系统和温度获取装置
CN105509901A (zh) * 2015-12-29 2016-04-20 温州大学 发光体相对色温的高精度计量方法
CN105509901B (zh) * 2015-12-29 2018-11-02 温州大学 发光体相对色温的高精度计量方法

Also Published As

Publication number Publication date
CN102749141A (zh) 2012-10-24
CN103048050B (zh) 2014-09-24

Similar Documents

Publication Publication Date Title
CN103048050B (zh) 一种测量目标真实温度的辐射测温方法和仪器
CN103076107B (zh) 基于太赫兹脉冲测量的燃烧温度传感装置及方法
US7661876B2 (en) Infrared target temperature correction system and method
CN102564610B (zh) 一种基于能量法的半透明材料高温辐射率测量装置及扣除背景辐射的修正方法
CN101793563B (zh) 多波段红外辐射自动测量系统
CN103063312B (zh) 一种测量物体发射率的测量系统及方法
CN112378616B (zh) 基于波长调制吸收光谱的高速流场多参数测量系统及方法
CN103675019B (zh) 一种红外热像仪快速测量材料表面发射率的方法
CN102889934A (zh) 实时测量温度的方法
CN101907492A (zh) 一种基于物体发射光谱的温度测量方法
CN103743489A (zh) 基于标准面源黑体的红外辐射计标定方法
CN106979832A (zh) 一种光纤分光测温系统及其测温方法
CN201716111U (zh) 多波段红外辐射自动测量系统
CN101419095A (zh) 灰体辐射率的测定方法
CN104864977A (zh) 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法
CN108801496B (zh) 一种基于交叠吸收谱的路径温度直方图测量系统与方法
CN102080990B (zh) 一种四波段高温测量装置及方法
US20160349113A1 (en) Characterization of absolute spectral radiance of an unknown ir source
CN105806491A (zh) 一种三波长二维温度场测量装置及方法
CN201464053U (zh) 一种基于双ccd的四通道温度场测量装置
CN110160660A (zh) 基于光场相机的高温部件多光谱测温方法及系统
CN103344601B (zh) 基于傅立叶变换红外光谱分析仪的半透明材料吸收系数测量方法
CN106644099B (zh) 一种多光谱温度场测量装置
CN201892573U (zh) 一种近红外辐射温度计
CN202814557U (zh) 一种实时测量非透明物体温度和光谱发射率的装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant