崔雙龍，孫博君，孫曉剛

哈爾濱工業(yè)大學(xué)儀器科學(xué)與工程學(xué)院，黑龍江 哈爾濱 150001

引 言

溫度是一個(gè)重要的物理參數(shù)。溫度測(cè)量一直是人們研究的熱點(diǎn)。紅外測(cè)溫技術(shù)是一種常用的光學(xué)測(cè)溫手段，具有響應(yīng)速度快、不破壞被測(cè)目標(biāo)溫度場(chǎng)等優(yōu)點(diǎn)。如何將探測(cè)器測(cè)得的能量更準(zhǔn)確的轉(zhuǎn)化為被測(cè)目標(biāo)的真實(shí)溫度一直是一個(gè)重要的研究方向，存在著許多未解的問(wèn)題，黑體輻射在某波段探測(cè)器上引起響應(yīng)的理論計(jì)算也是其中之一。

紅外測(cè)溫中，一直使用f(T)≈CTn模型來(lái)描述黑體波段輻射亮度響應(yīng)，該模型中n值的準(zhǔn)確性直接影響測(cè)溫結(jié)果的準(zhǔn)確性。對(duì)于不同的探測(cè)器，不同的波長(zhǎng)區(qū)間，不同的溫度范圍內(nèi)，有不同的C及n[1]。n值難以準(zhǔn)確獲取[2]，多數(shù)研究者都在使用Inagaki及Okamoto在1996年提出的三個(gè)擬合結(jié)果[3]，還不能很好的擴(kuò)展到任意波段的探測(cè)器。本文通過(guò)理論推導(dǎo)，提出了一種可以替代f(T)≈CTn模型的通用公式，通用公式與原模型f(T)≈CTn相比，最大的優(yōu)勢(shì)在于可以在任意波段內(nèi)，不需考慮溫度分區(qū)，直接計(jì)算求取黑體波段輻射亮度響應(yīng)f(T)，具有通用性。黑體波段輻射亮度響應(yīng)的通用公式進(jìn)一步完善了紅外測(cè)溫技術(shù)的基礎(chǔ)理論。

1 紅外測(cè)溫原理

紅外測(cè)溫中，探測(cè)器接收到的輻射能量包括三部分： 目標(biāo)自身發(fā)出的輻射，目標(biāo)反射環(huán)境的輻射和大氣輻射。測(cè)得的物體表面光譜輻射亮度為[4]

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

(1)

式(1)中，Lλ表示光譜輻射亮度，τaλ為光譜大氣透過(guò)率，ελ為被測(cè)目標(biāo)的光譜發(fā)射率，Ts為被測(cè)目標(biāo)表面溫度，ρλ為被測(cè)目標(biāo)的反射率，Tu為被測(cè)目標(biāo)周?chē)h(huán)境溫度，εaλ為大氣的光譜發(fā)射率，Ta為大氣溫度。在紅外測(cè)溫技術(shù)使用的窄波段范圍內(nèi)，ελ，τaλ，ρλ，εaλ可近似認(rèn)為與λ無(wú)關(guān)。當(dāng)被測(cè)目標(biāo)為朗伯體時(shí)，公式變化為[4]

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

(2)

式(2)中，Tr為探測(cè)器測(cè)得物體的輻射溫度。f(T)為某段波長(zhǎng)范圍內(nèi)，黑體輻射亮度在探測(cè)器上引起的響應(yīng)，本文稱(chēng)其為黑體波段輻射亮度響應(yīng)。f(T)可近似為f(T)≈CTn。

(3)

式(3)中，λ2及λ1為探測(cè)器響應(yīng)波段的上下限，Rλ為探測(cè)器的光譜響應(yīng)度，Mbλ(T)為黑體在波長(zhǎng)λ處溫度T時(shí)的輻射出射度，由普朗克公式計(jì)算。若被測(cè)目標(biāo)為灰體，即ε=α=1-ρ，對(duì)大氣可認(rèn)為εa=αa=1-τa，通過(guò)將f(T)≈CTn代入可得[4]

(4)

式(4)即為常用的紅外測(cè)溫公式。根據(jù)文獻(xiàn)[3]，273～473 K范圍內(nèi)，有效波長(zhǎng)2～5 μm的InSb探測(cè)器f(T)≈1.12×10-13×σT8.68，有效波長(zhǎng)6～9 μm的HgCdTe探測(cè)器f(T)≈5.16×10-5×σT5.53，有效波長(zhǎng)8～13 μm的HgCdTe探測(cè)器f(T)≈0.136×σT4.09，其中σ=5.6679×10-8[W·(m2·K4)-1]為斯蒂芬玻爾茲曼常數(shù)。

文獻(xiàn)[1-2]指出，f(T)≈CTn對(duì)于不同的探測(cè)器，不同的波長(zhǎng)區(qū)間，不同的溫度范圍內(nèi)，有不同的C及n。且難以確定。針對(duì)此現(xiàn)象，對(duì)Planck公式的近似公式“維恩近似公式”進(jìn)行積分，通過(guò)理論推導(dǎo)給出了一種可以替代f(T)≈CTn模型的通用公式，該通用公式可通過(guò)探測(cè)器的光譜響應(yīng)范圍及溫度，計(jì)算黑體波段輻射亮度響應(yīng)。

2 f(T)推導(dǎo)

根據(jù)式(3)對(duì)Planck公式進(jìn)行積分，未能得到有效的解析解，因而退一步對(duì)“維恩近似公式”進(jìn)行積分。對(duì)于普朗克公式當(dāng)ec2/λT?1時(shí)，可近似為式(5)[5]

(5)

其中c1和c2為普朗克輻射常數(shù)，c1=2πhc3=3.741 5×108(W·μm4·m-2)，c2=hc/KB=1.438 79×104(μm·K)。式(5)即為維恩近似公式。實(shí)際應(yīng)用中，在誤差可接受的情況下，維恩近似公式可以替代普朗克公式使用。Rλ一般可擬合為λ的二次多項(xiàng)式，即Rλ=aλ2+bλ+c，則黑體波段輻射亮度響應(yīng)為

(6)

設(shè)x=c2/λT則

(7)

特別的，當(dāng)Rλ=c為常數(shù)，即式(7)中a=b=0時(shí)

(8)

式(7)即為本工作提出的計(jì)算黑體波段輻射亮度響應(yīng)的通用公式，可計(jì)算任意波段內(nèi)的f(T)，將f(T)帶入式(2)后可通過(guò)二分法或迭代法計(jì)算被測(cè)目標(biāo)溫度Ts。

3 通用公式仿真

3.1 Rλ為常數(shù)時(shí)全波段仿真

斯蒂芬玻爾茲曼定律是將普朗克公式在全波段積分所得，表達(dá)式為Mbb=σT4，其中σ為斯蒂芬波爾茲曼常數(shù)，σ=5.667 9×10-8[W·(m2·K4)-1]。黑體波段輻射亮度響應(yīng)公式按全波段計(jì)算，即x1=0，x2=∞得到的結(jié)果理論上是與斯蒂芬玻爾茲曼定律近似的，可用于仿真驗(yàn)證。全波段積分得

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]，其偏差來(lái)自于維恩公式對(duì)普朗克公式的近似。將斯蒂芬玻爾茲曼定律與式(9)繪圖對(duì)比如圖1。

圖1 通用公式與斯蒂芬玻爾茲曼公式對(duì)比Fig.1 Comparison between the universal formula and the Stephen Boltzmann formula

3.2 Rλ為二次多項(xiàng)式時(shí)波段輻射仿真

文獻(xiàn)[3]給出了三種探測(cè)器在273～473 K內(nèi)f(T)≈CTn的擬合結(jié)果，取其中有效波長(zhǎng)8～13 μm的HgCdTe探測(cè)器進(jìn)行仿真對(duì)比。文獻(xiàn)中探測(cè)器響應(yīng)曲線(xiàn)及擬合結(jié)果如圖2。

根據(jù)式(7)進(jìn)行求解并將之與文中給出的f(T)≈0.136×σT4.09進(jìn)行對(duì)比，如圖3。從圖中可以看出通用公式計(jì)算結(jié)果與參考文獻(xiàn)擬合結(jié)果基本一致。

圖2 文獻(xiàn)[3]中探測(cè)器響應(yīng)曲線(xiàn)及擬合結(jié)果Fig.2 Detector response curves and fitting results in [3]

圖3 通用公式與f(T)=CTn對(duì)比Fig.3 Comparison between the universal formula and f(T)=CTn

3.3 仿真小結(jié)

通過(guò)將通用公式與斯蒂芬玻爾茲曼定律及文獻(xiàn)[3]中擬合結(jié)果的對(duì)比，說(shuō)明了通用公式的正確性。通用公式源于維恩近似公式，適用于ec2/λT?1的情況。在誤差允許的情況下，通過(guò)式(7)可計(jì)算任意波段內(nèi)黑體波段輻射亮度響應(yīng)，簡(jiǎn)單易行更具通用性。

4 通用公式實(shí)驗(yàn)驗(yàn)證

4.1 實(shí)驗(yàn)方法

為進(jìn)一步驗(yàn)證通用公式的有效性，設(shè)計(jì)實(shí)驗(yàn)進(jìn)行驗(yàn)證。實(shí)驗(yàn)方法為： 選取測(cè)量目標(biāo)，調(diào)節(jié)其表面溫度至某一定值且穩(wěn)定，通過(guò)鉑電阻測(cè)定其表面溫度Ts及環(huán)境溫度Tu，分別記錄。使用熱像儀對(duì)其進(jìn)行測(cè)量，將熱像儀發(fā)射率設(shè)置為1，測(cè)量目標(biāo)溫度，記錄示值為T(mén)r。根據(jù)Ts，Tu及Tr可計(jì)算目標(biāo)發(fā)射率ε(近距離測(cè)量，大氣透過(guò)率τ近似為1； 被測(cè)目標(biāo)近似為灰體，即ρ=1-ε)[1]

(10)

調(diào)節(jié)熱像儀的發(fā)射率，使其示值溫度T=Ts，記此時(shí)發(fā)射率ε0為目標(biāo)發(fā)射率參考值[6-7]。若ε≈ε0，則通用公式可以用于紅外測(cè)溫技術(shù)的工程實(shí)踐之中。

4.2 面源輻射體實(shí)驗(yàn)數(shù)據(jù)

以實(shí)驗(yàn)室內(nèi)面源輻射體為目標(biāo)進(jìn)行實(shí)驗(yàn)，已知該輻射體發(fā)射率高于0.9，按4.1所述實(shí)驗(yàn)方法進(jìn)行實(shí)驗(yàn)，連續(xù)測(cè)量22組數(shù)據(jù)，列于表1。

表1 面源輻射體實(shí)驗(yàn)數(shù)據(jù)Table 1 Surface source radiator experimental data

圖4 實(shí)驗(yàn)裝置Fig.4 Experimental device

5 結(jié) 論

紅外測(cè)溫技術(shù)中，多數(shù)研究者通過(guò)模型f(T)≈CTn來(lái)描述黑體波段輻射亮度響應(yīng)，式中n值準(zhǔn)確性直接影響測(cè)溫結(jié)果的準(zhǔn)確性，且難以準(zhǔn)確獲取。提出了一種黑體波段輻射亮度響應(yīng)的通用公式。并通過(guò)仿真說(shuō)明了通用公式的正確性。進(jìn)一步通過(guò)實(shí)驗(yàn)證明了，通用公式可以替代f(T)≈CTn模型，應(yīng)用于工程實(shí)踐。通用公式與原模型相比，最大的優(yōu)勢(shì)是可以在任意波段內(nèi)，不需考慮溫度分區(qū)，通過(guò)理論計(jì)算的方式求取f(T)，具有通用性。黑體波段輻射亮度響應(yīng)的通用公式進(jìn)一步完善了紅外測(cè)溫技術(shù)的基礎(chǔ)理論。