CN103900723B - 利用辐射谱的全部信息测火焰实际温度的方法及测量系统 - Google Patents

利用辐射谱的全部信息测火焰实际温度的方法及测量系统 Download PDF

Info

Publication number
CN103900723B
CN103900723B CN201410169861.2A CN201410169861A CN103900723B CN 103900723 B CN103900723 B CN 103900723B CN 201410169861 A CN201410169861 A CN 201410169861A CN 103900723 B CN103900723 B CN 103900723B
Authority
CN
China
Prior art keywords
measured
temperature
digitlization
actual temperature
physical model
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
CN201410169861.2A
Other languages
English (en)
Other versions
CN103900723A (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.)
Tianjin Yitong Electrical Technology Development Group Ltd
Original Assignee
Tianjin Yitong Transmission & Distribution Electric Technology Co ltd
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 Tianjin Yitong Transmission & Distribution Electric Technology Co ltd filed Critical Tianjin Yitong Transmission & Distribution Electric Technology Co ltd
Priority to CN201410169861.2A priority Critical patent/CN103900723B/zh
Publication of CN103900723A publication Critical patent/CN103900723A/zh
Priority to US15/306,493 priority patent/US11359967B2/en
Priority to PCT/CN2015/074899 priority patent/WO2015161727A1/zh
Priority to EP15782663.7A priority patent/EP3136065B1/en
Application granted granted Critical
Publication of CN103900723B publication Critical patent/CN103900723B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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/0096Radiation pyrometry, e.g. infrared or optical thermometry for measuring wires, electrical contacts or electronic systems
    • 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/02Constructional details
    • G01J5/025Interfacing a pyrometer to an external device or network; User interface
    • 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/60Radiation pyrometry, e.g. infrared or optical thermometry using determination of colour temperature
    • G01J5/601Radiation pyrometry, e.g. infrared or optical thermometry using determination of colour temperature using spectral scanning
    • 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/80Calibration

Abstract

本发明公开了一种利用辐射谱的全部信息测火焰实际温度的方法及测量系统,本发明应用“能级结构修正”进行更加理论化的数据处理,可以利用辐射谱的全部信息。采用键盘的输入方式或数据传输的输入方式,获取到“能级结构修正参数”,并最终获取到被测火焰更加准确的实际温度值。本发明有效地克服了人们在采用“多光谱测温法”时,必须经过计算工作量大的数据处理进行“辐射率修正”,才能得到火焰实际温度的缺陷。目前已有的“多光谱测温法”只能利用辐射谱中的几个“单色辐射本领”信息;本发明能够利用辐射谱中的全部,一般有几千个“单色辐射本领”的信息。

Description

利用辐射谱的全部信息测火焰实际温度的方法及测量系统
技术领域
本发明涉及仪器仪表领域的辐射谱测温方法及设备领域,尤其涉及一种利用辐射谱的全部信息测火焰实际温度的方法及测量系统。
背景技术
火焰温度测量技术目前是国内外研究的热门课题。它在节能减排、航空航天、战略和战术武器的开发等领域都具有重要的意义和应用价值。火焰温度的测量方法很多,主要有接触法、成像法、激光光谱法、声波法和辐射法等。从国际实用温标的角度来看,无论是用于温度标准的复现还是实际的测量,在银凝固点(1234.0开)以上的温区,辐射测温方法都是首选方法。辐射法分为“发射-吸收光谱法”和“高温计法”。“发射-吸收光谱法”的缺点是:结构庞大、造价昂贵、技术复杂、容易破坏被测火焰的组分和热平衡状态等。“高温计法”的缺点是:需要进行“辐射率修正”后才能知道火焰的真实温度。目前获得物体真温和“辐射率”数据最有效的方法是“多光谱测温法”,它是利用多个光谱下的火焰辐射亮度信息,经过数据处理进行“辐射率修正”得到火焰的真实温度。
根据黑体辐射理论,人们定义了三种光学温度。即1、根据黑体的“辐射总能量”来测量的“辐射温度”;2、根据黑体“光谱中能量的分布”来测量的“色温度”;3、根据黑体在“某一波段内的辐射能量”来测量的“亮温度”。对于理想黑体,上述光学温度与实际温度是相同的;但对于实际物体,上述光学温度与实际温度并不相同。
目前公知的,获得火焰真温和“辐射率”数据最有效的方法是“多光谱测温法”,它是利用多个(实际上只有几个)光谱下的火焰辐射亮度信息,经过数据处理进行“辐射率修正”得到火焰的真实温度。因为“辐射率修正”需要的计算工作量大,只能利用有限几个光谱信息,所以这种方法的测温准确度难以进一步提高。“多光谱测温法”属于根据黑体的“光谱中能量的分布”来测量的“色温度”法。
发明内容
本发明提供了一种利用辐射谱的全部信息测火焰实际温度的方法及测量系统,本发明提高了测温的准确度,降低了计算的复杂度,详见下文描述:
一种利用辐射谱的全部信息测火焰实际温度的方法,所述方法包括以下步骤:
(1)通过处于校准状态的实际温度测量系统中的数字化光谱仪,测出待测对象在任意一已知温度下的发射本领随波长λ变化的第一数字化分布谱,求出峰值波长λm
(2)将第一数字化分布谱输入到安装有第一物理模型的单片机或PC机上进行数据处理,用最小二乘法对第一数字化分布谱作曲线拟合,获取反映能级结构的参数X(1)、X(2);
(3)获取待测对象的维恩位移定律修正公式,并作为第二物理模型输入到实际温度测量系统中的单片机或PC机上;
(4)将验证后的反映能级结构参数的X(2))输入到所述第二物理模型中,实现对实际温度测量系统的校准;
(5)通过处于测温状态的实际温度测量系统中的光谱仪,测出待测对象在待测温度下的发射本领随波长λ变化的第二数字化分布谱,获得峰值波长λm的具体数值;
(6)将步骤(5)获得的峰值波长λm,输入处于测温状态的实际温度测量系统内的PC机或单片机,根据所述第二物理模型进行运算处理,获取待测对象的温度值;
(7)将待测对象的温度值通过显示器予以显示。
所述第一物理模型具体为:
在适用于理想黑体的普朗克公式中,添加X(1)、X(2)两个能级结构修正参数后形成的公式:
E ( λ , T ) = X ( 1 ) C 1 λ - 5 ( e X ( 2 ) C 2 λ T - 1 ) - 1
其中,E(λ,T)为非黑体发射的光谱辐射通量密度,单位为Wcm-2·μm-1;C1=3.741844×10-12W·cm-2,称为第一辐射常数;C2=1.438832544cm·K,称为第二辐射常数;λ为辐射波长,单位为μm;T为待测对象的任一已知固定温度,单位为K;X(1)、X(2)为能级结构参数,是在普朗克公式中添加的两个待定参数。
所述第二物理模型具体为:
λmT=X(2)×2.898×10-3(米*开)。
所述实际温度测量系统具有校准和测温两种状态,校准是为了获取待测对象反映能级结构的参数X(1)、X(2)的具体数值;测温是为了获取待测对象的温度值。
一种利用辐射谱的全部信息测火焰实际温度的测量系统,所述测量系统包括:
数字化光谱仪,具有校准和测温两种功能状态,当处于校准状态时,用于测出待测对象在任意一已知温度下的发射本领随波长λ变化的第一数字化分布谱,求出峰值波长λm;当处于测温状态时,用于测出待测对象在待测温度T下的发射本领随波长λ变化的第二数字化分布谱;
第一输入装置,用于将第一数字化分布谱输入到安装有第一物理模型的第一数据处理装置;
第一数据处理装置,采用最小二乘法对第一数字化分布谱作曲线拟合,获取反映能级结构的参数X(1)、X(2);获取待测对象的维恩位移定律修正公式,维恩位移定律修正公式作为第二物理模型;
第二输入装置,用于将第二物理模型输入到第二数据处理装置;
校准装置,用于将验证后的反映能级结构参数的X(2))输入到第二物理模型中,实现对数字化光谱仪的校准;
第三输入装置,用于将第二数字化分布谱输入到第二数据处理装置;
第二数据处理装置,用于由第二数字化分布谱获得峰值波长λm的具体数值;根据峰值波长λm、第二物理模型进行运算处理,获取待测对象的温度值;
显示装置,用于显示待测对象的温度值。
所述数字化光谱仪具有校准和测温两种状态,校准是为了获取待测对象反映能级结构的参数X(1)、X(2)的具体数值;测温是为了获取待测对象的温度值。
本发明提供的技术方案的有益效果是:目前已有的“比色法”,因为“辐射率修正”计算工作量大,一般只能选择有限几个单波长的“单色辐射本领”进行数据处理,有的研究还用到人工神经网络技术。本发明将只能适用于理想黑体的维恩位移定律,通过“能级结构修正”,改造成为能够适用于特定非黑体待测对象的“广义维恩位移定律”。理论化程度高、计算工作量小(测温时只需做简单的除法),可以将光谱仪接收到的多达几千个单波长“单色辐射本领”,全部用于数据处理,精度大为提高。此外,本发明除了具有一般“色温度法”抗环境干扰能力强的优点外,还有理论清晰、计算简单、测温准确度高、校准程序只需要在一个温度点进行等优势。上述各项优势都是利用自然规律本身实现的。因此,本发明提供的技术方案及仪器,可以推广到整个辐射测温领域。
附图说明
图1为一种利用辐射谱的全部信息测火焰实际温度的方法的流程图;
图2为一种利用辐射谱的全部信息测火焰实际温度的测量系统的结构示意图;
图3为理想二极管灯丝温度为1800K时,发射本领E(λ,T)按波长λ分布曲线,If(A)=0.55 X(1)=688.151507034063的示意图;
图4为理想二极管灯丝温度为1880K时,发射本领E(λ,T)按波长λ分布曲线,If(A)=0.60 X(1)=565.451265451300的示意图;
图5为理想二极管灯丝温度为1960K时,发射本领E(λ,T)按波长λ分布曲线,If(A)=0.65 X(1)=319.697031070087的示意图;
图6为理想二极管灯丝温度为2040K时,发射本领E(λ,T)按波长λ分布曲线,If(A)=0.70 X(1)=405.677545908357的示意图;
图7为理想二极管灯丝温度为2120K时,发射本领E(λ,T)按波长λ分布曲线,If(A)=0.75 X(1)=322.483340509147的示意图;
图8为理想二极管灯丝温度为2200K时,发射本领E(λ,T)按波长λ分布曲线,If(A)=0.80 X(1)=256.636032423250的示意图;
图9为理想二极管灯丝发射本领E(λ,T)按波长λ分布曲线的峰值波长λm,随温度T增高向短波方向移动。两者间的关系由“广义维恩位移定律”λmT=X(2)×2.898×10-3确定。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。
本发明将只能适用于理想黑体的维恩位移定律,通过“能级结构修正”,改造成为能够适用于特定非黑体待测对象的“广义维恩位移定律”。对于该非黑体待测对象,其“色温度”与实际温度相同,具体说明如下:
维恩位移定律的内容是:在理想黑体的,发射本领E(λ,T)按波长λ分布曲线上,E(λ,T)有一最大值称为峰值,即最大的单色发射本领。相应于这最大值的波长,用λm表示,叫做峰值波长。对于每一温度,有一条发射本领E(λ,T)按波长λ分布曲线。随着温度T的增高,λm向短波方向移动,两者间的关系由以下公式确定:
λmT=2.898×10-3(米*开) (1)
称为维恩位移定律。
说明:引用“物理学基本常数”计算更准确的常数值,公式(1)应改写为:λmT=2.8979507428634305162188796020618×10-3(米*开)。
由光谱仪测得峰值波长λm后,即可由公式(1)计算出温度T来。但这样算出的温度,只是待测对象的“色温度”(即假设待测对象是黑体时的温度),不是待测对象的实际温度。
下面具体阐明本发明怎样将只能适用于理想黑体的维恩位移定律,通过“能级结构修正”改造成为能够适用于特定非黑体待测对象的“广义维恩位移定律”。达到对于该非黑体待测对象其“色温度”与实际温度相同的目的,详见下文描述:
101:通过处于校准状态的实际温度测量系统中的数字化光谱仪测出待测对象在任意一已知温度下的发射本领随波长λ变化的第一数字化分布谱,并求出其中的峰值波长λm
实际应用时,数字化光谱仪的可测波长范围应将峰值波长λm包含在内;利用公式(1)可以估算是否满足上述要求。其中,任意一已知温度可由“标准仪器”测得,例如:标准光电高温计、标准光学高温计、标准铂铑30-铂铑6。
102:将第一数字化分布谱输入到安装有第一物理模型的单片机或PC机上进行数据处理,用最小二乘法对第一数字化分布谱作曲线拟合,获取反映能级结构的参数X(1)、X(2);
第一物理模型是指:在只能适用于理想黑体的普朗克公式中,添加X(1)、X(2)两个“能级结构修正参数”后形成的公式:
E ( λ , T ) = X ( 1 ) C 1 λ - 5 ( e X ( 2 ) C 2 λ T - 1 ) - 1 - - - ( 2 )
其中,E(λ,T)为非黑体发射的光谱辐射通量密度,单位为Wcm-2·μm-1;C1=3.741844×10-12W·cm-2,称为第一辐射常数;C2=1.438832544cm·K,称为第二辐射常数;λ为辐射波长,单位为μm;T为待测对象的任一已知固定温度,单位为K。X(1)、X(2)为“能级结构修正参数”,是在普朗克公式中添加的两个待定参数。
说明:C1、C2更准确的常数值,引用“物理学基本常数”计算应为:
C1=3.7418439644240131737525825×10-12W·cm-2
C2=1.4388325438316932513026737224237cm·K
为了提高测量的准确度,以上步骤101-102也可以在多个已知温度点进行,然后将得到的多个X(2)值平均。
103:获取待测非黑体的维恩位移定律修正公式,并作为第二物理模型输入到实际温度测量系统中的单片机或PC机上;
将只能适用于黑体的维恩位移定律(公式3)改写为能适用于待测非黑体的维恩位移定律修正公式,即公式(4):
λmT=2.898×10-3(米*开) (3)
λmT=X(2)×2.898×10-3(米*开) (4)
其中,维恩位移定律修正公式(4)的推证过程如下:
在公式(2)中引入
d X = - ( X ( 2 ) C 2 λ 2 T ) d λ = - ( TX 2 X ( 2 ) C 2 ) d λ
公式(2)变为
E ( λ , T ) = C 1 T 5 X 5 C 2 5 ( e x - 1 )
为了求出上式中极大值的位置,需要对它求微分并令其等于零。
d E ( λ , T ) d x = C 1 T 5 { ( e x - 1 ) 5 X 4 - X 5 e x } C 2 5 ( e x - 1 ) 2 = 0
由此得出满足上式的条件为:
5eX-XeX-5=解0此方程得Xm=4.965即或λmT=X(2)×1.44cm·K/4.965=X(2)×2.898×10-3(米*开),此即公式(4)。
104:将验证后的反映能级结构参数的X(2))输入到第二物理模型中,实现对实际温度测量系统的校准;
其中,对反映能级结构参数的X(2)的验证具体为:将步骤101中获得的λm和已知温度,以及步骤102获得的X(2)代入公式(4),验证正确性。当误差较小时,将验证后的反映能级结构参数的X(2))输入到第二物理模型中,实现对实际温度测量系统的校准;否则,即误差较大时,回到步骤101,重新进行校正过程。
另外,如果被测对象是透明火焰(辐射线状或带状光谱),在采集光谱样本时需要在火焰内人工加入烟粒,使其变为发光火焰(辐射连续光谱)。
105:通过处于测温状态的实际温度测量系统中的光谱仪,测出待测对象在待测温度T下的单色发射本领E(λ,T)随波长λ变化的第二数字化分布谱,从而获得峰值波长λm的具体数值;
其中,峰值波长λm应包含在光谱仪的可测波长范围之内,以便能直接测出它的数值;利用公式(1)可以估算是否满足上述要求。
106:将步骤105获得的峰值波长λm,输入处于测温状态的实际温度测量系统内的PC机或单片机,根据第二物理模型进行运算处理,获取被测对象的温度值;
107:将被测对象的温度值通过显示器予以显示。
其中,实际温度测量系统有确定被测对象反映能级结构参数、以及确定被测对象温度的两种功能。
一种利用辐射谱的全部信息测火焰实际温度的测量系统,参见图2,该系统包括:
数字化光谱仪,具有校准和测温两种功能状态,当处于校准状态时,用于测出待测对象在任意一已知温度下的发射本领随波长λ变化的第一数字化分布谱,求出峰值波长λm;当处于测温状态时,用于测出待测对象在待测温度T下的发射本领随波长λ变化的第二数字化分布谱;
第一输入装置,用于将第一数字化分布谱输入到安装有第一物理模型的第一数据处理装置;
第一数据处理装置,采用最小二乘法对第一数字化分布谱作曲线拟合,获取反映能级结构的参数X(1)、X(2);获取待测对象的维恩位移定律修正公式,维恩位移定律修正公式作为第二物理模型;
第二输入装置,用于将第二物理模型输入到第二数据处理装置;
校准装置,用于将验证后的反映能级结构参数的X(2))输入到第二物理模型中,实现对数字化光谱仪的校准;
第三输入装置,用于将第二数字化分布谱输入到第二数据处理装置;
第二数据处理装置,用于由第二数字化分布谱获得峰值波长λm的具体数值;根据峰值波长λm、第二物理模型进行运算处理,获取待测对象的温度值;
显示装置,用于显示待测对象的温度值。
本发明实施例对上述各部件不做限制,只要能完成上述功能的器件均可,例如:第一数据处理装置和第二数据处理装置可以为单片机或PC机;第一输入装置、第二输入装置和第三输入装置可以为USB接口或WIFI等无线传输方式;显示装置可以为显示器等。
由优选实例可以看出,X(2)和λmT基本上不随温度T变化,其微小变化为实验误差所致。这一结果直接证明了我们的理论预期:由黑体模型导出的普朗克公式,其函数形式可以适用于任何热辐射物体,各种热辐射体的区别可以由该函数形式中“反映能级结构的参数”不同来描述。因此,黑体辐射理论中,包括维恩位移定律在内的所有结论都可以推广到其它热辐射体。这是本发明的理论和实验基础。
本实验充分地揭示了自然规律的本来面目,暴露了“辐射率修正”方法的缺陷。现代量子论认为,热(光)辐射的本质是微观粒子量子跃迁。历史的经验表明,这种涉及物质内部微观过程的问题,经典理论(“辐射率修正”方法的理论基础)已经显得无能为力,只有用现代量子理论(本发明的理论基础)才能解决。
本优选实例采用的非黑体热辐射源为“理想二极管”灯丝。它是0.075毫米的细长钨丝,可以近似地视为“无限长直导线”。钨丝的“辐射率”为0.44,显然属于非黑体。理想二极管的灯丝温度可以由电流控制。厂家已事先设计并测出了理想二极管灯丝电流与温度之间的关系,如表1所示。表1中所列温度,是灯丝的实际温度。
表1 理想二极管灯丝电流与温度关系
If(A) 0.50 0.55 0.60 0.65 0.70 0.75 0.80
T(×103k) 1.72 1.80 1.88 1.96 2.04 2.12 2.20
表2 测温范围为1800开到2000开的实验数据
表3 测温范围为2000开到2200开的实验数据
由于条件限制,目前提供的实验数据中,λm值是由实测数据外推得到的,误差较大。但是,如果根据以上数据分两档来测温度,其测温准确度亦可在目前规定的“工作计量器具”允许误差范围之内。可比较仪器是“比色测温仪”,其允许的误差范围是0.01到0.025。
提供的实验数据计算处理如下:
第一档测温范围为1800开到2000开。
λmT的平均值为2427.191*10-6(米*开),因而采用的公式为:λmT=2427.191*10-6(米*开);由此计算得:
T(0.55A)=2427.191*10-6/1.32755*10-6=1828(开);绝对误差为:1828-1800=28(开)。相对误差为:28/1800=0.016。
T(0.60A)=2427.191*10-6/1.2903*10-6=1881(开);绝对误差为:1881-1880=1(开)。相对误差为:1/1880=5.3*10-4=0.00053。
T(0.65A)=2427.191*10-6/1.2583*10-6=1929(开);绝对误差为:1929-1960=-31(开)。相对误差为:-31/1960=-0.0158。
最大绝对误差为31开;相对误差为0.016,已在目前规定的“工作计量器具”允许误差范围之内。
第二档测温范围为2000开到2200开。
λmT的平均值为2514.938*10-6(米*开),因而采用的公式为:λmT=2514.938*10-6(米*开);由此计算得:
T(0.70A)=2514.938*10-6/1.2348*10-6=2037(开);绝对误差为:2037-2040=-3(开)。相对误差为:3/1960=0.0015。
T(0.75A)=2514.938*10-6/1.18485*10-6=2123(开);绝对误差为:2123-2120=3(开)。相对误差为:3/2120=0.0014。
T(0.80A)=2514.938*10-6/1.1427*10-6=2201(开);绝对误差为:2201-2200=1(开)。相对误差为:1/2200=0.0005。
最大绝对误差为3开;相对误差为0.0015,已经在目前规定的“计量标准器具”允许误差(2—15开)范围之内。
上述实验结果中,第一档误差较大;第二档误差较小。原因是:本实验所用光谱仪只能测1微米以下波长的光。第一档温度较低,峰值波长较长,为1.2583-1.32755微米。因此,实验数据中实测数较少、外推数较多,结果误差较大。第二档温度较高,峰值波长较短,为1.1427-1.2348微米。因此,实验数据中实测数较多、外推数较少,结果误差较小。
可以预测,如果规定“对于数字化光谱仪的要求是:其可测波长范围应将峰值波长λm包含在内”,那么“利用辐射谱的全部信息测火焰实际温度的测量系统”的测温准确度更会有所提高,完全可以在目前规定的辐射测温仪“计量标准器具”允许的误差(2—15开)范围之内。
说明:如果不分档,采用的公式应为:
λmT=X(2)*b=0.8550*2.897*10-3=2.455*10-3(米*开)
图3-图9为实验数据拟合曲线:粗线为测得的谱线,细线为拟合得到的曲线。横轴的单位是m。通过上述图形可以看出本发明取得了较好的实验效果,验证了本发明的可行性。
本发明实施例对各器件的型号除做特殊说明的以外,其他器件的型号不做限制,只要能完成上述功能的器件均可。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种利用辐射谱的全部信息测火焰实际温度的方法,其特征在于,所述方法包括以下步骤:
(1)通过处于校准状态的实际温度测量系统中的数字化光谱仪,测出待测对象在任意一已知温度下的发射本领随波长λ变化的第一数字化分布谱,求出峰值波长λm
(2)将第一数字化分布谱输入到安装有第一物理模型的单片机或PC机上进行数据处理,用最小二乘法对第一数字化分布谱作曲线拟合,获取反映能级结构的参数X(1)、X(2);
(3)获取待测对象的维恩位移定律修正公式,并作为第二物理模型输入到实际温度测量系统中的单片机或PC机上;
(4)将验证后的反映能级结构参数的X(2)输入到所述第二物理模型中,实现对实际温度测量系统的校准;
(5)通过处于测温状态的实际温度测量系统中的数字化光谱仪,测出待测对象在待测温度下的发射本领随波长λ变化的第二数字化分布谱,获得峰值波长λm的具体数值;
(6)将步骤(5)获得的峰值波长λm,输入处于测温状态的实际温度测量系统内的PC机或单片机,根据所述第二物理模型进行运算处理,获取待测对象的温度值;
(7)将待测对象的温度值通过显示器予以显示;
其中,所述第一物理模型具体为:
在适用于理想黑体的普朗克公式中,添加X(1)、X(2)两个能级结构修正参数后形成的公式:
E ( λ , T ) = X ( 1 ) C 1 λ - 5 ( e X ( 2 ) C 2 λ T - 1 ) - 1
其中,E(λ,T)为非黑体发射的光谱辐射通量密度,单位为Wcm-2·μm-1;C1=3.741844×10-12W·cm-2,称为第一辐射常数;C2=1.438832544cm·K,称为第二辐射常数;λ为辐射波长,单位为μm;T为待测对象的任一已知固定温度,单位为K;X(1)、X(2)为能级结构参数,是在普朗克公式中添加的两个待定参数。
2.根据权利要求1所述的一种利用辐射谱的全部信息测火焰实际温度的方法,其特征在于,所述第二物理模型具体为:
λmT=X(2)×2.898×10-3(米*开)。
3.根据权利要求1所述的一种利用辐射谱的全部信息测火焰实际温度的方法,其特征在于,所述实际温度测量系统具有校准和测温两种状态,校准是为了获取待测对象反映能级结构的参数X(1)、X(2)的具体数值;测温是为了获取待测对象的温度值。
4.一种利用辐射谱的全部信息测火焰实际温度的测量系统,其特征在于,所述测量系统包括:
数字化光谱仪,具有校准和测温两种功能状态,当处于校准状态时,用于测出待测对象在任意一已知温度下的发射本领随波长λ变化的第一数字化分布谱,求出峰值波长λm;当处于测温状态时,用于测出待测对象在待测温度T下的发射本领随波长λ变化的第二数字化分布谱;
第一输入装置,用于将第一数字化分布谱输入到安装有第一物理模型的第一数据处理装置;
第一数据处理装置,采用最小二乘法对第一数字化分布谱作曲线拟合,获取反映能级结构的参数X(1)、X(2);获取待测对象的维恩位移定律修正公式,维恩位移定律修正公式作为第二物理模型;
第二输入装置,用于将第二物理模型输入到第二数据处理装置;
校准装置,用于将验证后的反映能级结构参数的X(2)输入到第二物理模型中,实现对数字化光谱仪的校准;
第三输入装置,用于将第二数字化分布谱输入到第二数据处理装置;
第二数据处理装置,用于由第二数字化分布谱获得峰值波长λm的具体数值;根据峰值波长λm、第二物理模型进行运算处理,获取待测对象的温度值;
显示装置,用于显示待测对象的温度值;
其中,所述第一物理模型具体为:
在适用于理想黑体的普朗克公式中,添加X(1)、X(2)两个能级结构修正参数后形成的公式:
E ( λ , T ) = X ( 1 ) C 1 λ - 5 ( e X ( 2 ) C 2 λ T - 1 ) - 1
其中,E(λ,T)为非黑体发射的光谱辐射通量密度,单位为Wcm-2·μm-1;C1=3.741844×10-12W·cm-2,称为第一辐射常数;C2=1.438832544cm·K,称为第二辐射常数;λ为辐射波长,单位为μm;T为待测对象的任一已知固定温度,单位为K;X(1)、X(2)为能级结构参数,是在普朗克公式中添加的两个待定参数。
5.根据权利要求4所述的一种利用辐射谱的全部信息测火焰实际温度的测量系统,其特征在于,所述数字化光谱仪具有校准和测温两种状态,
校准是为了获取待测对象反映能级结构的参数X(1)、X(2)的具体数值;
测温是为了获取待测对象的温度值。
CN201410169861.2A 2014-04-25 2014-04-25 利用辐射谱的全部信息测火焰实际温度的方法及测量系统 Active CN103900723B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201410169861.2A CN103900723B (zh) 2014-04-25 2014-04-25 利用辐射谱的全部信息测火焰实际温度的方法及测量系统
US15/306,493 US11359967B2 (en) 2014-04-25 2015-03-23 Method for measuring actual temperature of flame by using all information of radiation spectrum and measurement system thereof
PCT/CN2015/074899 WO2015161727A1 (zh) 2014-04-25 2015-03-23 利用辐射谱的全部信息测火焰实际温度的方法及测量系统
EP15782663.7A EP3136065B1 (en) 2014-04-25 2015-03-23 Method for measuring actual temperature of flame using all information about radiation spectrum and measurement system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410169861.2A CN103900723B (zh) 2014-04-25 2014-04-25 利用辐射谱的全部信息测火焰实际温度的方法及测量系统

Publications (2)

Publication Number Publication Date
CN103900723A CN103900723A (zh) 2014-07-02
CN103900723B true CN103900723B (zh) 2016-08-24

Family

ID=50992189

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410169861.2A Active CN103900723B (zh) 2014-04-25 2014-04-25 利用辐射谱的全部信息测火焰实际温度的方法及测量系统

Country Status (4)

Country Link
US (1) US11359967B2 (zh)
EP (1) EP3136065B1 (zh)
CN (1) CN103900723B (zh)
WO (1) WO2015161727A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103900723B (zh) * 2014-04-25 2016-08-24 天津送变电易通电力科技有限公司 利用辐射谱的全部信息测火焰实际温度的方法及测量系统
CN104864977B (zh) * 2014-12-17 2018-02-06 西北工业大学 一种考虑多波长光谱辐射的火箭发动机燃气温度测试方法
CN108398191A (zh) * 2018-01-29 2018-08-14 南京理工大学 基于光纤光谱仪的高精度色温测量方法及系统
CN108470418A (zh) * 2018-04-02 2018-08-31 深圳汇创联合自动化控制有限公司 一种预警准确的建筑物火灾预警系统
CN111044565A (zh) * 2019-12-30 2020-04-21 西安近代化学研究所 一种燃烧转爆轰过程真温多光谱测量方法
CN112629694B (zh) * 2020-11-26 2023-01-13 中国大唐集团科学技术研究院有限公司火力发电技术研究院 一种燃煤电站炉膛的温度检测方法
CN112834051B (zh) * 2020-12-31 2023-01-24 华北电力大学 一种基于多光谱相机的辐射测温法
CN113401360B (zh) * 2021-06-16 2023-03-10 电子科技大学 一种基于多波段光学辐射测温的航空发动机涡轮盘温度测量装置
CN114046883B (zh) * 2021-10-27 2023-12-22 东南大学 火焰温度及碳烟浓度分布同时重建方法、装置及存储介质
CN117405262B (zh) * 2023-12-15 2024-02-23 常州泰斯科电子有限公司 一种温度测试仪的多点温度采集方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0335224A2 (en) * 1988-03-29 1989-10-04 Nippon Steel Corporation Radiation thermometry
CN101000264A (zh) * 2007-01-15 2007-07-18 哈尔滨工业大学 光谱极值测温方法
US20100246631A1 (en) * 2009-03-24 2010-09-30 K-Space Associates, Inc. Blackbody fitting for temperature determination
CN102374902A (zh) * 2010-08-11 2012-03-14 曹柏林 一种提高辐射温度计测温准确度的量子论修正方法
CN103123281A (zh) * 2011-11-17 2013-05-29 财团法人工业技术研究院 非接触式温度测量方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5132922A (en) 1989-01-12 1992-07-21 Massachusetts Institute Of Technology Emissivity independent multiwavelength pyrometer
US5265036A (en) * 1991-02-25 1993-11-23 United Technologies Corporation Turbine pyrometer system for correction of combustor fireball interference
EP1006345A1 (en) * 1998-12-04 2000-06-07 Jurij K Lingart Method of determining the actual temperature of real bodies
CN1162648C (zh) * 2001-09-29 2004-08-18 上海理工大学 一种火焰监测诊断测量方法及装置
CN101907492A (zh) * 2010-07-16 2010-12-08 浙江工业大学 一种基于物体发射光谱的温度测量方法
CN103644972B (zh) 2013-11-29 2017-04-12 西北核技术研究所 一种多光谱辐射测温装置及方法
CN103900723B (zh) * 2014-04-25 2016-08-24 天津送变电易通电力科技有限公司 利用辐射谱的全部信息测火焰实际温度的方法及测量系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0335224A2 (en) * 1988-03-29 1989-10-04 Nippon Steel Corporation Radiation thermometry
CN101000264A (zh) * 2007-01-15 2007-07-18 哈尔滨工业大学 光谱极值测温方法
US20100246631A1 (en) * 2009-03-24 2010-09-30 K-Space Associates, Inc. Blackbody fitting for temperature determination
CN102374902A (zh) * 2010-08-11 2012-03-14 曹柏林 一种提高辐射温度计测温准确度的量子论修正方法
CN103123281A (zh) * 2011-11-17 2013-05-29 财团法人工业技术研究院 非接触式温度测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Determination of the true temperature of emitted radiation bodies from generalized Wien"s displacement law;A I Fisenko等;《J.Phys.D:Appl.Phys.》;19991130;第32卷(第22期);第2882~2885页 *
固体火箭发动机羽焰测温方法研究;范传新等;《宇航计测技术》;20060831;第26卷(第4期);第18~22页 *

Also Published As

Publication number Publication date
CN103900723A (zh) 2014-07-02
EP3136065B1 (en) 2020-03-04
US11359967B2 (en) 2022-06-14
WO2015161727A1 (zh) 2015-10-29
US20170045398A1 (en) 2017-02-16
EP3136065A1 (en) 2017-03-01
EP3136065A4 (en) 2018-10-03

Similar Documents

Publication Publication Date Title
CN103900723B (zh) 利用辐射谱的全部信息测火焰实际温度的方法及测量系统
Christidis et al. Dramatically increasing chance of extremely hot summers since the 2003 European heatwave
CN107687899B (zh) 一种红外测温方法及系统
CN104571217B (zh) 航天器热真空试验温度控制方法及系统
CN108801474A (zh) 一种四光谱涡轮叶片辐射测温方法
Lukovic et al. LED-based Vis-NIR spectrally tunable light source-the optimization algorithm
Chehbouni et al. Directional effect on radiative surface temperature measurements over a semiarid grassland site
CN102374902A (zh) 一种提高辐射温度计测温准确度的量子论修正方法
CN107271047A (zh) 一种不均匀温度场的红外辐射能测试平台及测试方法
CN103336995B (zh) 一种百万千瓦级光伏发电基地实时测光网络的构建方法
CN101907492A (zh) 一种基于物体发射光谱的温度测量方法
CN103646175A (zh) 一种目标光谱辐射亮度的计算方法
CN105043555B (zh) 一种计算光谱发射率和真实温度的方法
Garrido et al. Characterization of the KTH high-flux solar simulator combining three measurement methods
Jin et al. Understanding the heterogeneity of soil moisture and evapotranspiration using multiscale observations from satellites, airborne sensors, and a ground-based observation matrix
CN104280119A (zh) 一种双列对消红外光谱仪的定标系统
Yu et al. Development of an absorption-corrected method for 3D computed tomography of chemiluminescence
CN106896079B (zh) 耐高温镍基合金材料的光谱发射率建模方法和测量系统
CN103542934B (zh) 一种光谱响应度校准方法及其装置
Ni et al. Estimation of solar-induced fluorescence using the canopy reflectance index
CN203551111U (zh) 一种光谱响应度校准装置
CN109598059A (zh) 一种基于代理模型的热防护系统优化设计方法及设计系统
CN106768382A (zh) 一种基于辅助光源的辐射温度测量装置及方法
CN103674247A (zh) 一种丰度值可控的草地混合光谱获取实验方法
CN100561150C (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
ASS Succession or assignment of patent right

Owner name: TIANJIN ELECTRIC POWER TRANSMISSION AND TRANSFORMA

Free format text: FORMER OWNER: LIU YOUXIANG

Effective date: 20141223

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 300102 NANKAI, TIANJIN TO: 300456 TANGGU, TIANJIN

TA01 Transfer of patent application right

Effective date of registration: 20141223

Address after: 300456 2-CO22, room 1, outer ring north road, Tianjin Airport Economic Zone

Applicant after: Tianjin Electric Power Transmission and Transformation Yitong Electric Power Technology Co., Ltd.

Address before: 300102, No. 58, Changjiang Road, Nankai District, Tianjin, No. 6

Applicant before: Liu Youxiang

C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20171211

Address after: No. 6, No. 58, Changjiang Road, Nankai District, Tianjin

Patentee after: Tianjin Yitong Electrical Technology Development Group Limited

Address before: 300456 2-CO22, room 1, outer ring north road, Tianjin Airport Economic Zone

Patentee before: Tianjin Electric Power Transmission and Transformation Yitong Electric Power Technology Co., Ltd.