黑体红外波段辐射亮度响应的通用公式
崔双龙, 孙博君, 孙晓刚*
哈尔滨工业大学仪器科学与工程学院, 黑龙江 哈尔滨 150001
*通讯联系人 e-mail: sxg@hit.edu.cn

作者简介: 崔双龙, 1988年生, 哈尔滨工业大学仪器科学与工程学院博士研究生 e-mail: cui.shuanglong@qq.com

摘要

以红外测温技术为背景, 着重研究了f(T)=λ1λ2RλL(T)dλCTn模型, 即某段波长范围内, 黑体辐射在探测器上引起的响应, 称其为黑体波段辐射亮度响应。 对于不同的探测器, 不同的波长区间, 不同的温度范围, 有不同的 C n n值难以准确获取, 多数研究者试验时使用Inagaki及Okamoto在1996年提出的三个固定波段模型, 不能很好的扩展到任意波段的探测器。 通过使用“维恩近似公式”代替普朗克公式, 从理论上推导出 f( T)的解析式, 得到了黑体波段辐射亮度响应的通用公式, 从而能够通过理论计算的方式, 求取任意波段内黑体波段辐射亮度响应 f( T)。 使用黑体波段辐射亮度响应的通用公式进行了两项仿真工作。 一是将通用公式在全波段内进行积分, 得到解析式 Mbb=5.238 5×10-8 T4, 并与斯蒂芬玻尔兹曼定律对比。 通用公式求得的系数 σ'=5.238 5×10-8与斯蒂芬玻尔兹曼常数 σ=5.667 9×10-8差值为0.429 4×10-8。 二是使用通用公式计算出8~13 μm波段内黑体波段辐射亮度响应 f( T), 并绘图与Inagaki及Okamoto文章中的拟合结果 f( T)≈0.136 σT4 .09进行对比, 结果曲线基本一致。 两项仿真说明了通用公式的正确性, 在此基础上, 进一步进行实验验证。 以实验室内面源辐射体为目标, 根据所提出的通用公式, 计算被测目标的发射率 ε, 并将之与目标发射率参考值 ε0对比。 面源辐射体实验结果:ε̅0=0.92为参考值,ε̅=0.93为测量值, 发射率误差为0.01。 实验误差较小, 说明所提出的通用公式可用于红外测温的工程实践中。 通用公式与原模型 f( T)≈ CTn相比, 最大的优势在于可以在任意波段内, 不需考虑温度分区, 通过理论计算的方式, 求取黑体波段辐射亮度响应 f( T), 具有通用性。 黑体波段辐射亮度响应的通用公式进一步完善了红外测温技术的基础理论。

关键词: 红外测温; 黑体波段辐射亮度响应; f( T)≈ CTn;; n
中图分类号:TN219 文献标志码:A
Universal Formula of Blackbody Waveband Radiation Brightness Response in the Infrared Temperature Measurement Technology
CUI Shuang-long, SUN Bo-jun, SUN Xiao-gang*
School of Instrumental Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
*Corresponding author
Abstract

Based on the infrared temperature measurement technology, this paper focuses on the modelf(T)=λ1λ2RλL(T)dλCTn. This model represents the response of blackbody radiation brightness on the detector in a certain wavelength range, which is called the blackbody waveband radiation brightness response in this paper. For different detectors, different wavelength ranges, different temperature ranges, there are different C and n. The value of n is difficult to obtain accurately. Most researchers used the three fitting results of C and n proposed by Inagaki T and Okamoto Y in 1996, and the results could not be well extended to any waveband detectors. In this paper, by using the Wien Approximation Formula to replace Planck’s formula, the analytical formula of f( T) is theoretically derived which is a universal formula of f( T), so that we can obtain the blackbody waveband radiation brightness response f( T) in any waveband through theoretical calculations. The universal formula is applied in simulations. Simulation 1: Integrate the universal formula in the whole waveband to obtain an analytical formula Mbb=5.238 5×10-8 T4, and compare it with Stephen Boltzmann’s law. The difference between the coefficient σ'=5.238 5×10-8 obtained by the universal formula and the Stefan Boltzmann constant σ=5.667 9×10-8 is 0.429 4×10-8. Simulation 2: Acquire the spectral responsivity of the HgCdTe detector with an effective wavelength of 8~13 μm from reference [2], calculate the blackbody waveband radiation brightness response f( T) of the detector in the 8~13 μm band, plot the result to compared with f( T)≈0.136× σT4. 09 shown in reference [2], and the results are basically the same. Two simulations illustrate the correctness of the analytical formula. On this basis, further experimental verification is performed. A surface source radiator in the laboratory was applied as the target. Through the universal formula, the emissivity of target ε can be calculated and compared with the reference value of true emissivity of target ε0. Experimental results of the surface source radiator:ε̅0=0.92 is the reference value of true target emissivity,ε̅0=0.93 is measured value, the error of emissivity is 0.01. The small error indicates that the model proposed in this paper can be used in the engineering practice of infrared temperature measurement technology. In this paper, a universal formula that can replace the model f( T)≈ CTn is proposed. Compared with the original model, the greatest advantage of the universal formula is that it can be applied in any waveband without considering the temperature partition and it is a universality to calculate the blackbody waveband radiation brightness response by theoretical calculation. The universal formula further perfects the basic theory of infrared temperature measurement technology.

Keyword: Infrared temperature measurement; Blackbody waveband radiation brightness response; f( T)≈ CTn; Value of n
引 言

温度是一个重要的物理参数。 温度测量一直是人们研究的热点。 红外测温技术是一种常用的光学测温手段, 具有响应速度快、 不破坏被测目标温度场等优点。 如何将探测器测得的能量更准确的转化为被测目标的真实温度一直是一个重要的研究方向, 存在着许多未解的问题, 黑体辐射在某波段探测器上引起响应的理论计算也是其中之一。

红外测温中, 一直使用f(T)≈ CTn模型来描述黑体波段辐射亮度响应, 该模型中n值的准确性直接影响测温结果的准确性。 对于不同的探测器, 不同的波长区间, 不同的温度范围内, 有不同的Cn[1]n值难以准确获取[2], 多数研究者都在使用Inagaki及Okamoto在1996年提出的三个拟合结果[3], 还不能很好的扩展到任意波段的探测器。 本文通过理论推导, 提出了一种可以替代f(T)≈ CTn模型的通用公式, 通用公式与原模型f(T)≈ CTn相比, 最大的优势在于可以在任意波段内, 不需考虑温度分区, 直接计算求取黑体波段辐射亮度响应f(T), 具有通用性。 黑体波段辐射亮度响应的通用公式进一步完善了红外测温技术的基础理论。

1 红外测温原理

红外测温中, 探测器接收到的辐射能量包括三部分: 目标自身发出的辐射, 目标反射环境的辐射和大气辐射。 测得的物体表面光谱辐射亮度为[4]

Lλ=τaλελLλ(Ts)+τaλρλLλ(Tu)+εaλLλ(Ta)(1)

式(1)中, Lλ 表示光谱辐射亮度, τ aλ 为光谱大气透过率, ε λ 为被测目标的光谱发射率, Ts为被测目标表面温度, ρ λ 为被测目标的反射率, Tu为被测目标周围环境温度, ε aλ 为大气的光谱发射率, Ta为大气温度。 在红外测温技术使用的窄波段范围内, ε λ , τ aλ , ρ λ , ε aλ 可近似认为与λ 无关。 当被测目标为朗伯体时, 公式变化为[4]

f(Tr)=τaεf(Ts)+τaρf(Tu)+εaf(Ta)(2)

式(2)中, Tr为探测器测得物体的辐射温度。 f(T)为某段波长范围内, 黑体辐射亮度在探测器上引起的响应, 本文称其为黑体波段辐射亮度响应。 f(T)可近似为f(T)≈ CTn

f(T)=λ1λ2RλLbλ(T)dλ=λ1λ2RλM(T)πdλCTn(3)

式(3)中, λ 2λ 1为探测器响应波段的上下限, Rλ 为探测器的光谱响应度, Mbλ (T)为黑体在波长λ 处温度T时的辐射出射度, 由普朗克公式计算。 若被测目标为灰体, 即ε =α =1-ρ , 对大气可认为ε a=α a=1-τ a, 通过将f(T)≈ CTn代入可得[4]

Ts=1ε1τaTrn-(1-ε)Tun-εaτaTan1n(4)

式(4)即为常用的红外测温公式。 根据文献[3], 273~473 K范围内, 有效波长2~5 μ m的InSb探测器f(T)≈ 1.12× 10-13× σ T8.68, 有效波长6~9 μ m的HgCdTe探测器f(T)≈ 5.16× 10-5× σ T5.53, 有效波长8~13 μ m的HgCdTe探测器f(T)≈ 0.136× σ T4.09, 其中σ =5.6679× 10-8[W· (m2· K4)-1]为斯蒂芬玻尔兹曼常数。

文献[1, 2]指出, f(T)≈ CTn对于不同的探测器, 不同的波长区间, 不同的温度范围内, 有不同的Cn。 且难以确定。 针对此现象, 对Planck公式的近似公式“ 维恩近似公式” 进行积分, 通过理论推导给出了一种可以替代f(T)≈ CTn模型的通用公式, 该通用公式可通过探测器的光谱响应范围及温度, 计算黑体波段辐射亮度响应。

2 f(T)推导

根据式(3)对Planck公式进行积分, 未能得到有效的解析解, 因而退一步对“ 维恩近似公式” 进行积分。 对于普朗克公式当 ec2/λT≫1时, 可近似为式(5)[5]

MbλT=c1λ51ec2λT-1c1λ5e-c2λT(5)

其中c1c2为普朗克辐射常数, c1=2π hc3=3.741 5× 108 (W· μ m4· m-2), c2=hc/KB=1.438 79× 104 (μ m· K)。 式(5)即为维恩近似公式。 实际应用中, 在误差可接受的情况下, 维恩近似公式可以替代普朗克公式使用。 Rλ 一般可拟合为λ 的二次多项式, 即Rλ =2++c, 则黑体波段辐射亮度响应为

f(T)=λ1λ21πc1λ5e-c2λT(aλ2++c)dλ(6)

x=c2/λ T

f(T)=c1c24πT4-c(x3+3x2+6x+6)ex-ac22(1+x)+bc2T(x2+2x+2)T2exx1x2(7)

特别的, 当Rλ =c为常数, 即式(7)中a=b=0时

f(T)=c1c24πT4Rλ{[-e-x(x3+3x2+6x+6)]|x1x2}(8)

式(7)即为本工作提出的计算黑体波段辐射亮度响应的通用公式, 可计算任意波段内的f(T), 将f(T)带入式(2)后可通过二分法或迭代法计算被测目标温度Ts

3 通用公式仿真
3.1 Rλ 为常数时全波段仿真

斯蒂芬玻尔兹曼定律是将普朗克公式在全波段积分所得, 表达式为Mbb=σ T4, 其中σ 为斯蒂芬波尔兹曼常数, σ =5.667 9× 10-8[W· (m2· K4)-1]。 黑体波段辐射亮度响应公式按全波段计算, 即x1=0, x2=∞ 得到的结果理论上是与斯蒂芬玻尔兹曼定律近似的, 可用于仿真验证。 全波段积分得

Mbb=5.2385×10-8T4(9)

σ '=5.238 5× 10-8[W· (m2· K4)-1]与σ =5.667 9× 10-8[W· (m2· K4)-1]接近, 差值为0.429 4× 10-8[W· (m2· K4)-1], 其偏差来自于维恩公式对普朗克公式的近似。 将斯蒂芬玻尔兹曼定律与式(9)绘图对比如图1。

图1 通用公式与斯蒂芬玻尔兹曼公式对比Fig.1 Comparison between the universal formula and the Stephen Boltzmann formula

3.2 Rλ 为二次多项式时波段辐射仿真

文献[3]给出了三种探测器在273~473 K内f(T)≈ CTn的拟合结果, 取其中有效波长8~13 μ m的HgCdTe探测器进行仿真对比。 文献中探测器响应曲线及拟合结果如图2。

图2 文献[3]中探测器响应曲线及拟合结果Fig.2 Detector response curves and fitting results in [3]

根据式(7)进行求解并将之与文中给出的f(T)≈ 0.136× σ T4.09进行对比, 如图3。 从图中可以看出通用公式计算结果与参考文献拟合结果基本一致。

图3 通用公式与f(T)=CTn对比Fig.3 Comparison between the universal formula and f(T)=CTn

3.3 仿真小结

通过将通用公式与斯蒂芬玻尔兹曼定律及文献[3]中拟合结果的对比, 说明了通用公式的正确性。 通用公式源于维恩近似公式, 适用于 ec2/λT≫1的情况。 在误差允许的情况下, 通过式(7)可计算任意波段内黑体波段辐射亮度响应, 简单易行更具通用性。

4 通用公式实验验证
4.1 实验方法

为进一步验证通用公式的有效性, 设计实验进行验证。 实验方法为: 选取测量目标, 调节其表面温度至某一定值且稳定, 通过铂电阻测定其表面温度Ts及环境温度Tu, 分别记录。 使用热像仪对其进行测量, 将热像仪发射率设置为1, 测量目标温度, 记录示值为Tr。 根据Ts, TuTr可计算目标发射率ε (近距离测量, 大气透过率τ 近似为1; 被测目标近似为灰体, 即ρ =1-ε )[1]

ε=f(Tr)-f(Tu)f(Ts)-f(Tu)(10)

调节热像仪的发射率, 使其示值温度T=Ts, 记此时发射率ε 0为目标发射率参考值[6, 7]。 若ε ε 0, 则通用公式可以用于红外测温技术的工程实践之中。

4.2 面源辐射体实验数据

以实验室内面源辐射体为目标进行实验, 已知该辐射体发射率高于0.9, 按4.1所述实验方法进行实验, 连续测量22组数据, 列于表1

表1 面源辐射体实验数据 Table 1 Surface source radiator experimental data

根据参考值ε 0的测量数据, 有效数字取小数点后2位。 通过拉伊达准则进行粗大误差判别, 最后将有效数据求取均值, ε̅0=0.92作为参考值, ε̅=0.93作为测量值, 发射率误差为0.01。 面源辐射体实验装置如图4。

图4 实验装置Fig.4 Experimental device

通过实验, 验证了通用公式的正确性。 实验结果 ε̅ε̅0, 说明黑体波段辐射亮度响应的通用公式可以用于工程实践。

5 结 论

红外测温技术中, 多数研究者通过模型f(T)≈ CTn来描述黑体波段辐射亮度响应, 式中n值准确性直接影响测温结果的准确性, 且难以准确获取。 提出了一种黑体波段辐射亮度响应的通用公式。 并通过仿真说明了通用公式的正确性。 进一步通过实验证明了, 通用公式可以替代f(T)≈ CTn模型, 应用于工程实践。 通用公式与原模型相比, 最大的优势是可以在任意波段内, 不需考虑温度分区, 通过理论计算的方式求取f(T), 具有通用性。 黑体波段辐射亮度响应的通用公式进一步完善了红外测温技术的基础理论。

参考文献
[1] LI Yan-feng, ZHANG Zhi-jie, ZHAO Chen-yang, et al(李岩峰, 张志杰, 赵晨阳, ). Chinese Journal of Sensors and Actuators(传感技术学报), 2017, 30(9): 1348. [本文引用:3]
[2] BAI Jing-chen, YU Qing-bo, HU Xian-zhong, et al(白敬晨, 于庆波, 胡贤忠, ). Journal of Northeastern University·Natural Science(东北大学学报·自然科学版), 2013, 34(12): 1747. [本文引用:2]
[3] Inagaki T, Okamoto Y. NDT & E International, 1996, 29(6): 363. [本文引用:4]
[4] LI Huan, LI Yun-hong, HUI Shi-ke, et al(李欢, 李云红, 惠史科, ). Computer and Digital Engineering(计算机与数字工程), 2017, (11): 223. [本文引用:3]
[5] SHI Xiao-guang, HUAN Ke-wei, GAO Lan-lan(石晓光, 宦克为, 高兰兰). Infrared Physics(红外物理). Hangzhou: Zhejiang University Press(杭州: 浙江大学出版社), 2013. [本文引用:1]
[6] LI Wen-jun, XU Yong-da, ZHENG Yong-jun(李文军, 徐永达, 郑永军). China Measurement & Test(中国测试), 2017, 43(6): 12. [本文引用:1]
[7] WANG Zhong-ren, LI Dong-yang, RUAN Ban-chao, et al(王中任, 黎冬阳, 阮班超, ). Machinery Design & Manufacture(机械设计与制造), 2016, 6: 226. [本文引用:1]