CN108593115B - 粉尘环境中高精度红外测温方法、装置、设备及存储介质 - Google Patents
粉尘环境中高精度红外测温方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN108593115B CN108593115B CN201810234561.6A CN201810234561A CN108593115B CN 108593115 B CN108593115 B CN 108593115B CN 201810234561 A CN201810234561 A CN 201810234561A CN 108593115 B CN108593115 B CN 108593115B
- Authority
- CN
- China
- Prior art keywords
- dust
- temperature
- testee
- atmosphere
- infrared
- 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
- 239000000428 dust Substances 0.000 title claims abstract description 302
- 238000005259 measurement Methods 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000003860 storage Methods 0.000 title claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 7
- 230000001052 transient effect Effects 0.000 claims abstract description 5
- 230000002452 interceptive effect Effects 0.000 claims abstract description 4
- 230000005855 radiation Effects 0.000 claims description 114
- 230000003595 spectral effect Effects 0.000 claims description 18
- 230000003287 optical effect Effects 0.000 claims description 15
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 12
- 238000010521 absorption reaction Methods 0.000 claims description 12
- 239000002245 particle Substances 0.000 claims description 10
- 238000002310 reflectometry Methods 0.000 claims description 9
- 238000009529 body temperature measurement Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 7
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 6
- 239000001569 carbon dioxide Substances 0.000 claims description 6
- 239000000843 powder Substances 0.000 claims description 6
- 230000008033 biological extinction Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000004891 communication Methods 0.000 claims description 2
- 239000006185 dispersion Substances 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- 229910000831 Steel Inorganic materials 0.000 description 6
- 239000010959 steel Substances 0.000 description 6
- 238000001931 thermography Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000003723 Smelting Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000003595 mist Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 241000209035 Ilex Species 0.000 description 1
- 229910000661 Mercury cadmium telluride Inorganic materials 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001447 compensatory effect Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- WPYVAWXEWQSOGY-UHFFFAOYSA-N indium antimonide Chemical compound [Sb]#[In] WPYVAWXEWQSOGY-UHFFFAOYSA-N 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000004571 lime Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation 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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- 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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
-
- 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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明实施例提供了一种粉尘环境中高精度红外测温方法,包括:根据粉尘透射率模型、粉尘作用下的大气透射率模型及被测物体表面的真实温度模型,获取粉尘环境中红外测温补偿模型;采用所述粉尘环境中红外测温补偿模型,测量并得到粉尘环境中的高精度温度数据。本发明实施例还提供了一种主动交互设备、装置及非暂态可读存储介质,用来实现所述方法。本发明可以获取高精度的红外测温的温度数据。
Description
技术领域
本发明实施例涉及温度检测技术领域,尤其涉及一种粉尘环境中高精度红外测温方法、装置、设备及存储介质。
背景技术
温度是钢铁流程中至关重要的参数之一,在钢铁制造过程中检测并调控温度可以避免产品缺陷从而保证产品质量。在温度检测设备中,非制冷焦平面红外热像仪是一种常用的测温设备,具有非接触、非侵入和实时连续测温等优点。然而,红外热像仪的测温结果容易受到各种各样的环境因素的影响,尤其是在工业过程中,现场环境复杂且多变,存在高温、高粉尘、多振动等扰动,对红外热像仪的测温结果造成较大的误差,限制了红外热像仪的广泛应用。
针对环境因素对于红外测温精确度的影响,可以通过温度补偿的方式来提高红外测温的精确度,补偿方式可以分为两大类:硬件补偿和软件补偿。硬件补偿的方式包括设计补偿电路或者自动调节红外镜头,这些补偿手段主要是为了克服红外设备工作温度对测温结果的影响。软件补偿的方式通常利用数值拟合、神经网络等算法,定量预测需要补偿的温度值,主要解决测温距离环境温度对红外测温结果的影响。上述补偿方法只是针对一些影响因素进行补偿,难以解决其他环境因素的干扰。在钢铁冶炼过程中,红外热像仪通常安装在某一固定位置且工作温度变化不大,影响红外测温准确度的主要因素是生产现场中存在的粉尘,因此,如何克服粉尘对测温准确度的影响,在粉尘环境下仍然能够测量并得到高精度的红外测温数据,就成为业界亟待解决的问题。
发明内容
针对现有技术存在的上述问题,本发明实施例提供了一种粉尘环境中高精度红外测温方法、装置、设备及存储介质。
一方面,本发明实施例提供了一种粉尘环境中高精度红外测温方法,包括:根据粉尘透射率模型、粉尘作用下的大气透射率模型及被测物体表面的真实温度模型,建立粉尘环境中红外测温补偿模型;采用所述粉尘环境中红外测温补偿模型,测量并得到粉尘环境中的高精度温度数据。
另一方面,本发明实施例提供了一种主动交互设备及一种非暂态可读存储介质。所述一种主动交互设备包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中:所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行所述一种粉尘环境中高精度红外测温方法。所述一种非暂态可读存储介质存储程序指令,用于执行所述一种粉尘环境中高精度红外测温方法。
最后,本发明实施例提供了一种粉尘环境中高精度红外测温装置,其特征在于,包括:粉尘辐射能获取模块、红外辐射能计算模块、被测物体真实温度计算模块及红外测温补偿模型建立模块;
所述粉尘辐射能获取模块,用于获取红外热像仪接收到的粉尘的辐射能并发送至所述红外辐射能计算模块;
所述红外辐射能计算模块,用于根据所述粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能并发送至所述被测物体真实温度计算模块;
所述被测物体真实温度计算模块,用于根据所述红外辐射能,计算粉尘环境下红外热像仪测得被测物体的真实温度并发送至所述红外测温补偿模型建立模块;
所述红外测温补偿模型建立模块,用于根据所述粉尘环境下红外热像仪测得被测物体的真实温度,并建立粉尘环境下的红外测温补偿模型。
本发明实施例提供了一种粉尘环境中高精度红外测温方法、装置、设备及存储介质,通过在粉尘环境下建立红外测温的补偿模型,可以获取高精度的红外测温的温度数据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例中粉尘环境中高精度红外测温方法的整体流程图;
图2是本发明实施例中建立粉尘环境中红外测温补偿模型示意图;
图3是现有技术中粉尘环境下红外热像仪测温结果示意图;
图4是本发明实施例中粉尘环境下补偿后的红外测温结果示意图;
图5是本发明实施例中对测温结果补偿后的红外热像仪误差示意图;
图6是本发明实施例的硬件设备工作示意图;
图7是本发明实施例中粉尘环境中高精度红外测温装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供了一种粉尘环境中高精度红外测温方法、装置、设备及存储介质。参见图1,图1是本发明实施例中粉尘环境中高精度红外测温方法的整体流程图,包括:
S101:根据粉尘透射率模型、粉尘作用下的大气透射率模型及被测物体表面的真实温度模型,建立粉尘环境中红外测温补偿模型。
其中,所述粉尘透射率模型、粉尘作用下的大气透射率模型及被测物体表面的真实温度模型的实现方案如下:
(1)建立粉尘透射率模型并计算粉尘作用下的大气透射率
红外辐射在大气中的传输能力可以用透射率来表征。被测目标的红外辐射在光路上传播时会衰减,主要是大气中存在的水蒸气、二氧化碳等的吸收以及云雾等微粒的散射。当被测目标和红外镜头的光路中存在粉尘时,会严重影响到红外测温的精度。光谱透射率τ(λ)可以基于布盖尔-朗伯定律描述:
τ(λ)=I(λ)/I0(λ)=exp[-k(v,n)·c·l] (1)
式中,τ(λ)为粉尘透过率,I(λ)为透射辐射强度,I0(λ)为入射辐射强度,k(v,n)为粉尘的质量消光系数,c为粉尘的质量浓度,l粉尘的厚度,也就是红外辐射通道长度。粉尘对不同波长的红外光的消减也符合Lambert-Beer定律。v=2πrf/c,f为入射辐射频率,r为粉尘微粒的半径,c为入射辐射在粉尘中的传输速率。
由式(1)可知,当质量浓度、质量消光系数、辐射通道长度确定后,就可以计算出粉尘的透过率。由于粉尘的反射系数可以忽略不计,因此,质量消光系数主要由吸收系数和散射系数组成,即:
k(v,n)=α(v,n)+δ(v,n) (2)
式中α(v,n)为粉尘的吸收系数,δ(v,n)为粉尘的散射系数。n=nr-ini为粉尘的复数折射率,nr表征粉尘散射部分,ni表征粉尘吸收部分。
粉尘微粒的粒径总有一个分布,典型的有对数分布、正态分布、威布尔分布和洛森-莱姆莱尔分布等。本文实验中的高炉粉尘微粒符合对数分布:
式中,r为粉尘微粒半径;rg为粉尘微粒的集合平均半径;σ为粉尘半径的标准差;β为拟合参量。
对于吸收系数α(v,n)和散射系数δ(v,n),可按照下面公式计算:
式中,Qa(v,n)为粉尘的吸收因子;Qs(v,n)为粉尘的散射因子;Qa(v,n)、Qs(v,n)与粉尘的衰减因子Qe(v,n)之间存在如下关系:
Qe(v,n)=Qa(v,n)+Qs(v,n) (5)
Qa(v,n)和Qe(v,n)可按常用的Vander Hulst近似计算:
式中,k=2v(nr-1),b=arctan[ni/(ni-1)]。
联立式(2)~(6)可得粉尘的质量消光系数k(v,n),再将k(v,n)带入式(1)中,即可求得粉尘的透射率。
在钢铁冶炼工厂中,不存在云雾等气象环境,因此,在计算钢铁冶炼工业现场中的大气透射率时,忽略云雾等微粒的影响,只考虑大气中影响红外辐射的三个最重要的因素:水蒸气、二氧化碳和粉尘。当被测目标与红外镜头之间的光路上不存在粉尘时,影响大气透射率的因素主要是水蒸气和二氧化碳,大气透射率的计算公式如式(7)所示。
式中和分别为水蒸气和二氧化碳单独作用下的透射率。当被测目标与红外镜头之间的光路上存在粉尘时,影响大气透射率的因素主要是水蒸气、二氧化碳和粉尘。红外辐射的大气透射率是这三个因素共同作用的结果,因此,大气透射率应由式(8)计算。
式中τdust表示粉尘透射率,可以由所建粉尘透射率模型求得。
红外热像仪一般工作在长红外波段,将其工作波段记为[λ1,λ2]。为便于后续计算,用表示该工作波段内的平均大气透射率。将该波段平均分为n个长度为dλ的小波段,n为正整数,即
利用等面积原则可计算平均大气透射率,具体见公式(9)
在后续的红外测温机理模型推导过程中,使用的是平均大气透射率。
(2)建立红外测温模型
假设光谱吸收率αλ、光谱反射率ρλ和光谱透射率τλ。一般来说,在一定温度下,入射到一个物体表面的红外辐射,一部分能量被吸收,一部分能量被反射,还有一部分能量经过物体透射过去,即会发生吸收、反射、透射三种物理现象。假设投射到物体表面上的辐射能量为1,按能量守恒定理,有光谱吸收率αλ、光谱反射率ρλ和光谱透射率τλ之和为1,即
αλ+ρλ+τλ=1 (10)
默认被测物体是不透明物体,故其透过率τλ=0,根据式(10)可知
ρλ=1-αλ (11)
根据基尔霍夫定律,在同样的温度下,各种不同物体对相同波长的单色辐射出射度与单色吸收比之比值都相等,并等于该温度下黑体对同一波长的单色辐射出射度。即
ελ=αλ (12)
所以光谱反射率为
ρλ=1-ελ (13)
对于大气,认为其反射率ρa=0,则大气的吸收率可以表示为
αa=1-τa (14)
同理,根据基尔霍夫定律可得,εa=αa,
εa=1-τa (15)
设T0被测物体表面真实温度,Tr红外热像仪的辐射温度,Ta大气温度,Te环境温度,ε被测物体表面发射率,εa为大气的发射率,τa大气透射率,Wλ表示红外热像仪接收到的辐射能,W0被测物体自身辐射能,We被测物体反射周围环境的辐射能,Wa大气辐射能。
图2描述了光路上不存在粉尘时红外热像仪接收到的红外辐射能量。用Wλ表示红外热像仪接收到的红外辐射能量,它由三个部分构成,来自被测物体自身辐射能、被测物体反射的周围环境的辐射能和大气的辐射能,见公式(16)。
Wλ=ετaW0+ρλτaWe+εaWa (16)
设A0为红外热像仪最小空间张角所对应的目标的可视面积,d为被测目标到红外热像仪之间的距离,S为红外热像仪透镜的面积,R为红外探测器的光谱响应度,表示红外探测器把红外辐射能转变为电信号的能力。
则与辐射能量相对应的红外热像仪信号电压可以表示为
VS=SA0d-2{τa[ε∫△λW0(λ,T)Rdλ+ρλ∫△λWe(λ,T)Rdλ]+εa∫△λWa(λ,T)Rdλ} (17)
通常情况下A0d-2为常数。
令k=SA0d-2,f(t)=∫△λW(λ,T)Rdλ,则
VS=k{τa[εf(T0)+ρλf(Te)]+εaf(Ta)} (18)
又因为ρλ=1-ελ,εa=1-τa,所以式(18)可以改写为式(19)。
VS=k{τa[εf(T0)+(1-ελ)f(Te)]+(1-τa)f(Ta)} (19)
令VS/k=f(Tr),则式(19)可变换为式(20)
f(Tr)=τa[εf(T0)+(1-ελ)f(Te)]+(1-τa)f(Ta) (20)
当被测物体为黑体,即ελ=1,且大气透射率τa=1时,则
f(Tr)=f(T0) (21)
即红外热像仪测得的辐射温度等于被测物体的真实温度。
当被测物体不是黑体,即ελ<1,且大气透射率τa<1时,则
f(Tr)≠f(T0) (22)
红外热像仪测得的辐射温度不等于被测物体的真实温度。
根据普朗克定律可知,
当探测器的光谱响应度R不随波长λ变化时,对式(23)积分可得到
f(T)=cTn (24)
式中c为常值,n取决于不同的探测元材料,对于HgCdTe(8~13μm)探测器,n=4.09,对于InSb(2~5μm)探测器,n=8.68。
将式(24)代入式(20),可得辐射温度的计算公式,见式(25)。
Tr n=τa[εT0 n+(1-ελ)Te n]+(1-τa)Ta n (25)
对式(25)进行变换,可得被测物体表面的真实温度计算公式(26)。
S102:采用所述粉尘环境中红外测温补偿模型,测量并得到粉尘环境中的高精度温度数据。
在钢铁冶炼过程应用红外热像仪检测温度时,当红外热像仪的安装参数固定后,影响红外测温精度的主要是粉尘,而粉尘的存在导致大气透射率τa变小,因此,需要分析大气透射率对测温精确度的影响。图3描述了光路上存在粉尘时红外热像仪接收到的红外辐射能量。
先对式(26)进行对数变换,求其关于大气透射率τa的偏微分,可得当大气透射率变化时,造成的红外测温误差。
设τd是光路上存在粉尘时的大气透射率,Trd是粉尘影响下的红外热像仪的辐射温度,T0d是被测物体的真实温度。
其中,所述粉尘环境中红外测温补偿模型根据粉尘环境下红外热像仪接收到的红外辐射能、粉尘环境下红外热像仪测得的物体表面温度及粉尘透射率,对被测物体的表面温度进行修正补偿,得到高精度红外测温数据
参见图2,图2是本发明实施例中建立粉尘环境中红外测温补偿模型示意图,包括:
S201:获取红外热像仪接收到的粉尘的辐射能,包括:
获取粉尘透射率、粉尘影响下的大气透射率及粉尘辐射能;
根据所述粉尘透射率、粉尘影响下的大气透射率及粉尘辐射能,获取红外热像仪接收到的粉尘的辐射能。
当红外热像仪与被测对象之间的光路存在粉尘时,粉尘除了影响大气透射率外,粉尘自身也有一定温度,因此,红外热像仪也会接收到来自粉尘的辐射能。设粉尘的温度为Td,粉尘的发射率为εdust,粉尘的影响下的大气透射率为τd,粉尘的辐射能为Wd,则红外热像仪接收到的粉尘的辐射能可以表示为
Wλd=εdustτdWd (28)
S202:根据所述粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能,包括:获取被测物体自身辐射能、被测物体反射的周围环境的辐射能及大气的辐射能;
根据所述被测物体自身辐射能、被测物体反射的周围环境的辐射能、大气的辐射能及粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能。
设Wλ表示红外热像仪接收到的红外辐射能量,τd为粉尘影响下的大气透射率。当红外热像仪与被测对象的光路间存在粉尘时,它由四个部分构成,来自被测物体自身辐射能W0、被测物体反射的周围环境的辐射能We、大气的辐射能Wa以及粉尘的辐射能Wd,原始的计算红外辐射能的式(16)应修正为式(29)。
Wλ=ετdW0+εdustτdWd+ρτdWe+εaWa (29)
式中ε为被测物体表面发射率,εa为大气的发射率,ρ为被测对象的光谱反射率,τd粉尘影响下的大气透射率。
S203:根据所述红外辐射能,计算粉尘环境下的被测物体的真实温度,包括:
采用电压信号描述所述红外辐射能,计算粉尘环境下的被测物体的真实温度。所述红外辐射能包括:被测对象的表面温度、粉尘自身温度、大气温度及周围环境温度。
设A0为红外热像仪最小空间张角所对应的目标的可视面积,d为被测目标到红外热像仪之间的距离,S为红外热像仪透镜的面积,R为红外探测器的光谱响应度,表示红外探测器把红外辐射能转变为电信号的能力。根据红外辐射能的计算公式(29),与红外热像仪接收到的辐射能量相对应的红外热像仪信号电压VS可以表示为
式中通常情况下A0d-2为常数。令k=SA0d-2,f(t)=∫△λW(λ,T)Rdλ,则
VS=k{τd[εf(T0)+εdustf(Td)+ρf(Te)]+εaf(Ta)} (31)
又因为ρ=1-ε,εa=1-τa,所以式(31)可以改写为式(32)。
VS=k{τd[εf(T0)+εdustf(Td)+(1-ε)f(Te)]+(1-τa)f(Ta)} (32)
令VS/k=f(Trd),则式(32)可变换为式(33)
f(Trd)=τd[εf(T0)+εdustf(Td)+(1-ελ)f(Te)]+(1-τa)f(Ta) (33)
式(33)描述了存在粉尘时红外辐射温度与被测对象的表面温度、粉尘自身温度、大气温度、周围环境温度之间的关系。
设αdust为粉尘的光谱吸收率,ρdust为粉尘的光谱反射率,τdust为粉尘的光谱透射率。由公式(10)可知,αdust+ρdust+τdust=1,同时可以忽略粉尘的光谱反射率ρdust,可以算得粉尘的光谱吸收率为
αdust=1-τdust
又因为εdust=αdust,所以εdust=1-τdust,再根据普朗克定律式(23)和式(24),同理可推导出式(34)。
Trd n=τd[εT0d n+(1-τdust)Td n+(1-ε)Te n]+(1-τa)Ta n (34)
求解式(34)可知,被测物体的真实温度T0d应该表示为:
式中Trd为红外辐射温度,T0d为粉尘影响下的被测对象的表面真实温度,Td为粉尘的温度,Te为周围环境温度,Ta为大气温度,ε为被测对象的表面发射率,τd为粉尘影响下的大气透射率,τdust为粉尘透射率,n取决于不同的探测元材料。
S204:根据所述真实温度,建立粉尘环境下的红外测温补偿模型,包括:
获取粉尘影响下的红外热像仪的辐射温度;
根据所述被测物体的真实温度及所述粉尘影响下的红外热像仪的辐射温度,建立粉尘环境下的红外测温补偿模型。
红外热像仪在测温应用中,当被测物体与红外镜头之间存在粉尘时,依旧采用原始的大气透射率,按式(26)计算被测物体的温度,即
式中T′0表示当光路上存在粉尘时红外热像仪测得的被测物体的真实温度。
由式(35)和(36)可知,当利用红外热像仪检测被测物体的温度时,粉尘导致的大气透射率的改变会造成较大的测温误差,因此,需要对原始的红外测温结果进行补偿。
红外热像仪的辐射温度反映了红外镜头接收到红外辐射能量,而粉尘的存在导致大气透射率减小,使红外镜头接收到红外辐射能量减小。根据粉尘影响下红外热像仪测得的温度T′0,由式(36)可推导出粉尘影响下的红外热像仪的辐射温度。
Trd n=τa[ε0T0'n+(1-ε0)Te n]+(1-τa)Ta n (37)
在研究粉尘对红外测温结果的影响时,认为环境温度Te、大气温度Ta和被测物体的发射率ε是已知的,并且当红外热像仪选定后,n就是常数。当测量光路上存在粉尘时,应按照式(35)计算被测物体的真实温度。因此,将式(37)代入式(35)中,可得
由式(7)和(8)可知,
故式(39)可化简为
式中T′0表示当光路上存在粉尘时红外热像仪测得的被测物体的表面真实温度,T0d为粉尘影响下的被测对象的表面真实温度
式(40)便是红外热像仪测温的补偿公式。当被测物体与红外镜头之间存在粉尘时,根据红外热像仪的测温结果,采用式(40)对粉尘影响下的红外测温结果进行补偿。
当被测对象与红外测温镜头的光路间存在粉尘时,若已知被测物体的发射率ε、环境温度Te、大气温度Ta和红外热像仪的检测温度T′0,且τdust可根据粉尘透射率模型求得,依据式(40)可以实时地对红外热像仪的测温结果进行补偿,从机理上克服粉尘造成的误差,提高红外热像仪的测温精确度。
参见图3,图3是现有技术中粉尘环境下红外热像仪测温结果示意图,包括:
粉尘环境下温度301、无粉尘环境下温度302及黑体温度303。由图中可见,未采用本发明技术方案的前提下,测量到的温度误差较大。
参见图4,图4是本发明实施例中粉尘环境下补偿后的红外测温结果示意图,包括:
无粉尘温度401、粉尘环境下补偿后温度402及黑体温度403。由图中可见,通过采用本发明的技术方案,可以显著缩小测量到的温度误差。
参见图5,图5是本发明实施例中对测温结果补偿后的红外热像仪误差示意图,包括:
有粉尘误差501、补偿后误差502及无粉尘误差503。由图中可见,通过采用本发明的技术方案,补偿后误差502与无粉尘误差503十分接近,在初始阶段甚至出现补偿后误差502比无粉尘误差503更低的情况。
参见图6,图6是本发明实施例的硬件设备工作示意图,所述硬件设备包括:一种粉尘环境中高精度红外测温设备601、处理器602及存储介质603。
粉尘环境中高精度红外测温设备601:所述一种粉尘环境中高精度红外测温设备601实现所述一种粉尘环境中高精度红外测温方法。
处理器602:所述处理器602加载并执行所述存储介质603中的指令及数据用于实现所述的一种粉尘环境中高精度红外测温方法。
存储介质603:所述存储介质603存储指令及数据;所述存储介质603用于实现所述的一种粉尘环境中高精度红外测温方法。
参见图7,图7是本发明实施例中粉尘环境中高精度红外测温装置示意图,包括:
粉尘辐射能获取模块701、红外辐射能计算模块702、被测物体真实温度计算模块703及红外测温补偿模型建立模块704;
所述粉尘辐射能获取模块701,用于获取红外热像仪接收到的粉尘的辐射能并发送至所述红外辐射能计算模块702;
所述红外辐射能计算模块702,用于根据所述粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能并发送至所述被测物体真实温度计算模块703;
所述被测物体真实温度计算模块703,用于根据所述红外辐射能,计算粉尘环境下红外热像仪测得被测物体的真实温度并发送至所述红外测温补偿模型建立模块704;
所述红外测温补偿模型建立模块704,用于根据所述粉尘环境下红外热像仪测得被测物体的真实温度,并建立粉尘环境下的红外测温补偿模型。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (6)
1.一种粉尘环境中高精度红外测温方法,其特征在于,包括:
根据粉尘透射率模型、粉尘作用下的大气透射率模型及被测物体表面的真实温度模型,获取粉尘环境中红外测温补偿模型;
采用所述粉尘环境中红外测温补偿模型,测量并得到粉尘环境中的高精度温度数据;
其中,所述粉尘透射率模型具体为:
τ(λ)=I(λ)/I0(λ)=exp[-k(v,n)·c·l];
式中,τ(λ)为粉尘透过率,I(λ)为透射辐射强度,I0(λ)为入射辐射强度,k(v,n)为粉尘的质量消光系数,c为粉尘的质量浓度,l粉尘的厚度;k(v,n)=α(v,n)+δ(v,n),α(v,n)为粉尘的吸收系数,δ(v,n)为粉尘的散射系数,n=nr-ini为粉尘的复数折射率,nr表征粉尘散射部分,ni表征粉尘吸收部分;v=2πrf/c,f为入射辐射频率,r为粉尘微粒的半径,c为入射辐射在粉尘中的传输速率;
式中,Qa(v,n)为粉尘的吸收因子;Qs(v,n)为粉尘的散射因子;
式中,r为粉尘微粒半径;rg为粉尘微粒的集合平均半径;σ为粉尘半径的标准差;β为拟合参量;
式中,k=2v(nr-1),b=arctan[ni/(ni-1)],Qe(v,n)=Qa(v,n)+Qs(v,n);
所述大气透射率模型具体为:
式中和分别为水蒸气和二氧化碳单独作用下的透射率,τdust表示粉尘透射率;
所述被测物体表面的真实温度模型具体为:
式中,T0被测物体表面真实温度,Tr红外热像仪的辐射温度,Ta大气温度,Te环境温度,ε被测物体表面发射率,τa大气透射率,n取决于不同的探测元材料;
所述获取粉尘环境中红外测温补偿模型,具体包括:求解所述被测物体表面的真实温度模型关于大气透射率τa的偏微分,确定当大气透射率变化时,造成的红外测温误差;
所述采用所述粉尘环境中红外测温补偿模型,测量并得到粉尘环境中的高精度温度数据,具体包括:
所述粉尘环境中红外测温补偿模型根据粉尘环境下红外热像仪接收到的红外辐射能、粉尘环境下红外热像仪测得的物体表面温度及粉尘透射率,对被测物体的表面温度进行修正补偿,得到高精度红外测温数据。
2.根据权利要求1所述的方法,其特征在于,所述根据粉尘透射率模型、粉尘作用下的大气透射率模型及被测物体表面的真实温度模型,获取粉尘环境中红外测温补偿模型,包括:
获取红外热像仪接收到的粉尘的辐射能;
根据所述粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能;
根据所述红外辐射能,计算粉尘环境下红外热像仪测得被测物体的真实温度;
根据所述真实温度,建立粉尘环境下的红外测温补偿模型;
其中,所述获取红外热像仪接收到的粉尘的辐射能,包括:
获取粉尘透射率、粉尘影响下的大气透射率及粉尘辐射能;
根据所述粉尘透射率、粉尘影响下的大气透射率及粉尘辐射能,获取红外热像仪接收到的粉尘的辐射能,具体通过如下公式实现:
Wλd=εdustτdWd
式中,εdust为粉尘的发射率,τd为粉尘的影响下的大气透射率,Wd为粉尘的辐射能,Wλd为红外热像仪接收到的粉尘的辐射能;
所述根据所述粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能,包括:
获取被测物体自身辐射能、被测物体反射的周围环境的辐射能及大气的辐射能;
根据所述被测物体自身辐射能、被测物体反射的周围环境的辐射能、大气的辐射能及粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能,具体通过如下公式实现:
Wλ=ετdW0+εdustτdWd+ρτdWe+εaWa
式中,ε为被测物体表面发射率,εa为大气的发射率,ρ为被测对象的光谱反射率,τd为粉尘影响下的大气透射率,W0为被测物体自身辐射能、We为被测物体反射的周围环境的辐射能、Wa为大气的辐射能,Wd为粉尘的辐射能;
所述根据所述红外辐射能,计算粉尘环境下的被测物体的真实温度,包括:
采用电压信号描述所述红外辐射能,计算粉尘环境下的被测物体的真实温度,具体方法为:
VS=SA0d-2{τd[ε∫△λW0(λ,T)Rdλ+εdust∫△λWd(λ,T)Rdλ+ρ∫△λWe(λ,T)Rdλ]+εa∫△λWa(λ,T)Rdλ}
式中,A0为红外热像仪最小空间张角所对应的目标的可视面积,d为被测目标到红外热像仪之间的距离,S为红外热像仪透镜的面积,R为红外探测器的光谱响应度,W0(λ,T)为波长λ温度T时被测物体自身辐射能,We(λ,T)为波长λ温度T时被测物体反射的周围环境的辐射能,Wd(λ,T)为波长λ温度T时粉尘的辐射能,Wa(λ,T)为波长λ温度T时大气的辐射能;
所述根据所述真实温度,建立粉尘环境下的红外测温补偿模型,具体包括:
当被测物体与红外镜头之间存在粉尘时,根据红外热像仪的测温结果,采用如下公式对粉尘影响下的红外测温结果进行补偿:
式中,T0'表示当光路上存在粉尘时红外热像仪测得的被测物体的表面真实温度,T0d为粉尘影响下的被测对象的表面真实温度,Td为粉尘的温度。
3.根据权利要求1或2所述的方法,其特征在于,所述红外辐射能包括:被测对象的表面红外辐射能、粉尘的红外辐射能、大气温度及周围环境的辐射能。
4.一种主动交互设备,其特征在于,包括:
至少一个处理器;以及
与所述处理器通信连接的至少一个存储器,其中:
所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至3任一所述方法。
5.一种非暂态可读存储介质,其特征在于,所述非暂态可读存储介质存储程序指令,所述程序指令用于执行如权利要求1至3任一所述方法。
6.一种粉尘环境中高精度红外测温装置,其特征在于,包括:
粉尘辐射能获取模块,用于获取红外热像仪接收到的粉尘的辐射能并发送至红外辐射能计算模块;
红外辐射能计算模块,用于根据所述粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能并发送至被测物体真实温度计算模块;
被测物体真实温度计算模块,用于根据所述红外辐射能,计算粉尘环境下红外热像仪测得被测物体的真实温度并发送至红外测温补偿模型建立模块;
红外测温补偿模型建立模块,用于根据所述粉尘环境下红外热像仪测得被测物体的真实温度,并建立粉尘环境下的红外测温补偿模型;
其中,所述粉尘辐射能获取模块具体用于:
获取粉尘透射率、粉尘影响下的大气透射率及粉尘辐射能;
根据所述粉尘透射率、粉尘影响下的大气透射率及粉尘辐射能,获取红外热像仪接收到的粉尘的辐射能,具体通过如下公式实现:
Wλd=εdustτdWd
式中,εdust为粉尘的发射率,τd为粉尘的影响下的大气透射率,Wd为粉尘的辐射能,Wλd为红外热像仪接收到的粉尘的辐射能;
所述红外辐射能计算模块具体用于:
获取被测物体自身辐射能、被测物体反射的周围环境的辐射能及大气的辐射能;
根据所述被测物体自身辐射能、被测物体反射的周围环境的辐射能、大气的辐射能及粉尘的辐射能,计算粉尘环境下红外热像仪接收到的红外辐射能,具体通过如下公式实现:
Wλ=ετdW0+εdustτdWd+ρτdWe+εaWa
式中,ε为被测物体表面发射率,εa为大气的发射率,ρ为被测对象的光谱反射率,τd为粉尘影响下的大气透射率,W0为被测物体自身辐射能、We为被测物体反射的周围环境的辐射能、Wa为大气的辐射能,Wd为粉尘的辐射能;
所述被测物体真实温度计算模块具体用于:
采用电压信号描述所述红外辐射能,计算粉尘环境下的被测物体的真实温度,具体方法为:
VS=SA0d-2{τd[ε∫△λW0(λ,T)Rdλ+εdust∫△λWd(λ,T)Rdλ+ρ∫△λWe(λ,T)Rdλ]+εa∫△λWa(λ,T)Rdλ}
式中,A0为红外热像仪最小空间张角所对应的目标的可视面积,d为被测目标到红外热像仪之间的距离,S为红外热像仪透镜的面积,R为红外探测器的光谱响应度,W0(λ,T)为波长λ温度T时被测物体自身辐射能,We(λ,T)为波长λ温度T时被测物体反射的周围环境的辐射能,Wd(λ,T)为波长λ温度T时粉尘的辐射能,Wa(λ,T)为波长λ温度T时大气的辐射能;
所述红外测温补偿模型建立模块具体用于:
当被测物体与红外镜头之间存在粉尘时,根据红外热像仪的测温结果,采用如下公式对粉尘影响下的红外测温结果进行补偿:
式中,T0'表示当光路上存在粉尘时红外热像仪测得的被测物体的表面真实温度,T0d为粉尘影响下的被测对象的表面真实温度,τdust为粉尘透射率,n=nr-ini为粉尘的复数折射率,nr表征粉尘散射部分,ni表征粉尘吸收部分;Te环境温度;Td为粉尘的温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810234561.6A CN108593115B (zh) | 2018-03-21 | 2018-03-21 | 粉尘环境中高精度红外测温方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810234561.6A CN108593115B (zh) | 2018-03-21 | 2018-03-21 | 粉尘环境中高精度红外测温方法、装置、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108593115A CN108593115A (zh) | 2018-09-28 |
CN108593115B true CN108593115B (zh) | 2019-10-11 |
Family
ID=63626948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810234561.6A Active CN108593115B (zh) | 2018-03-21 | 2018-03-21 | 粉尘环境中高精度红外测温方法、装置、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108593115B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109238473B (zh) * | 2018-11-02 | 2020-08-11 | 汕头大学 | 一种红外图像阵列传感器灰尘干扰进行校准的装置与方法 |
CN111047211B (zh) * | 2019-12-23 | 2023-12-19 | 河北山宝环境工程有限责任公司 | 粉尘数据综合分析系统及分析方法 |
CN111272296B (zh) * | 2020-02-19 | 2020-12-01 | 中南大学 | 一种降低光路中粉尘对红外测温影响的校正方法及系统 |
CN111238650B (zh) * | 2020-03-11 | 2021-04-13 | 中南大学 | 一种基于特殊红外光谱的高炉铁水测温方法及系统 |
CN111521270A (zh) * | 2020-04-23 | 2020-08-11 | 烟台艾睿光电科技有限公司 | 一种体温筛查报警系统及其工作方法 |
CN111751003B (zh) * | 2020-06-10 | 2022-12-13 | 四川省东宇信息技术有限责任公司 | 一种热像仪温度修正系统、方法及热像仪 |
CN111811662B (zh) * | 2020-07-21 | 2023-08-08 | 绍兴晨璞网络科技有限公司 | 一种高温烤箱用非接触红外热成像装置及其热力补偿方法 |
CN113721529A (zh) * | 2021-09-03 | 2021-11-30 | 成都理工大学 | 一种基于fpga的数据采集系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879121B (zh) * | 2012-09-13 | 2014-10-22 | 北京百纳威尔科技有限公司 | 测量温度的方法及通讯设备 |
CN102853916B (zh) * | 2012-09-26 | 2014-08-27 | 中国神华能源股份有限公司 | 一种用于对煤垛表面进行远距离红外测温的方法及系统 |
CN103604504B (zh) * | 2013-10-15 | 2016-06-29 | 中国人民解放军海军工程大学 | 一种红外辐射精确测温方法 |
CN103971337A (zh) * | 2014-04-29 | 2014-08-06 | 杭州电子科技大学 | 一种利用大气传输特征的红外图像去雾霾方法 |
CN106441586B (zh) * | 2016-08-23 | 2019-07-05 | 上海华乘智能设备有限公司 | 基于温湿度自动补偿的红外带电检测装置及其补偿方法 |
-
2018
- 2018-03-21 CN CN201810234561.6A patent/CN108593115B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108593115A (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108593115B (zh) | 粉尘环境中高精度红外测温方法、装置、设备及存储介质 | |
Wang et al. | Compensation method for infrared temperature measurement of explosive fireball | |
CN209247174U (zh) | 自定标热像检测仪 | |
Huang et al. | A new wavelength selection criterion for two-color pyrometer interfered with participating media | |
Li et al. | A new CMTF evaluation model for dynamic target in photoelectric imaging system | |
Meca et al. | Calculation and optimisation of the maximum uncertainty in infrared temperature measurements taken in conditions of high uncertainty in the emissivity and environment radiation values | |
JPH02134522A (ja) | H↓2o吸収補正型放射温度計 | |
Shanks et al. | Stokes resolved differential temperature: an important metric of polarimetric precision in the long-wave infrared | |
US20190137386A1 (en) | Gas observation method | |
Kruczek | Conditions for use of long-wave infrared camera to measure the temperature of the sky | |
Guoqing et al. | Nonlinear atmospheric correction based on neural network for infrared target radiometry | |
Semenov et al. | Optico-electronic dust meter | |
Li et al. | Gas imaging detectivity model combining leakage spot size and range | |
Liu et al. | Research on evaluation index of stray light suppression ability of star sensor based on signal-to-noise ratio | |
Hao et al. | A method for testing the infrared spectral radiation of aircraft engine exhaust plume | |
Ko et al. | Atmospheric characterization on the Kennedy Space Center Shuttle Landing Facility | |
Rogers et al. | Reporting NETD: why measurement techniques matter | |
Wang et al. | Research of suppressing stray light of cold shield with different black coating | |
Zhang et al. | High temperature imaging dynamic measurement techniques under harsh experimental environment | |
Sentenac et al. | Evaluation of a charge-coupled-device-based video sensor for aircraft cargo surveillance | |
Huang et al. | Research on accurate non-contact temperature measurement method for telescope mirror | |
Liu et al. | Applicability analysis of operating range model of airborne infrared detection system | |
Li | Detection distance calculation model of flying target and atmospheric influence analysis | |
JP6969542B2 (ja) | 温度計測システム及び温度計測方法 | |
Chen et al. | Measurement of atmospheric particle number density fluctuations based on moiré deflectometry |
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 |