劉連偉,楊淼淼,樊宏杰,許振領(lǐng),姚 梅
(中國人民解放軍63892部隊(duì),河南洛陽471003)
·紅外技術(shù)及應(yīng)用·
一種表面發(fā)射率的測量方法研究
劉連偉,楊淼淼,樊宏杰,許振領(lǐng),姚 梅
(中國人民解放軍63892部隊(duì),河南洛陽471003)
大型實(shí)驗(yàn)設(shè)備由于表面溫度難以控制而使其表面發(fā)射率精確測量較為困難。在分析表面發(fā)射率測量誤差主要來源的基礎(chǔ)上,提出一種表面發(fā)射率工程測量方法。分析并推導(dǎo)了到達(dá)探測器上的輻射量,利用兩發(fā)射率已知的設(shè)備,實(shí)現(xiàn)了環(huán)境輻射量的測量,從而消除了環(huán)境輻射對(duì)表面發(fā)射率測量的影響,有效提高了測量精度。設(shè)計(jì)了實(shí)驗(yàn)方法并進(jìn)行了發(fā)射率測量實(shí)驗(yàn),實(shí)驗(yàn)結(jié)果表明未考慮環(huán)境輻射的測量方法與該方法的計(jì)算結(jié)果相對(duì)誤差可達(dá)20.37%,驗(yàn)證了該測量方法的可行性和準(zhǔn)確性。
發(fā)射率;紅外熱像儀;環(huán)境輻射
表面發(fā)射率是目標(biāo)紅外輻射特性建模與仿真的重要參數(shù),發(fā)射率測量精度直接影響目標(biāo)紅外輻射特性計(jì)算的精度。目前發(fā)射率測量實(shí)驗(yàn)方法包括量熱法、反射率法、能量法和多波長法等[1-2],這些方法大多需要對(duì)被測量表面進(jìn)行精確控溫或?qū)y量環(huán)境要求較高,然而許多設(shè)備如車輛等體積龐大、結(jié)構(gòu)復(fù)雜,難以精確控制其表面溫度,一般的利用紅外熱像儀或紅外測溫儀測量此類設(shè)備表面發(fā)射率的工程方法多是根據(jù)熱像儀定標(biāo)系數(shù)計(jì)算物體表面輻射亮度,再結(jié)合物體表面溫度反推其發(fā)射率[3-6],但此種方法沒有考慮物體反射環(huán)境輻射對(duì)測量結(jié)果的影響,文獻(xiàn)[7]中對(duì)影響發(fā)射率測量精度的因素進(jìn)行了簡要分析,但并未給出適用于紅外熱像儀的發(fā)射率測量方法。
本文提出一種基于紅外熱像儀的表面波段發(fā)射率工程測量方法,理論分析并推導(dǎo)了考慮環(huán)境輻射的發(fā)射率計(jì)算公式,提高了發(fā)射率測量精度。
發(fā)射率測量的工程方法一般將被測表面近似為朗伯體,本文提出的方法遵循這一近似;文中所述發(fā)射率測量實(shí)驗(yàn)熱像儀與被測物體距離小于10 m,故忽略大氣衰減對(duì)測量結(jié)果的影響[7]。
2.1 傳統(tǒng)測量方法
測量表面發(fā)射率的傳統(tǒng)工程方法為利用紅外熱像儀測量并求出被測設(shè)備紅外輻射亮度,然后根據(jù)被測設(shè)備溫度求出表面發(fā)射率,具體計(jì)算步驟為:
(1)對(duì)紅外熱像儀進(jìn)行定標(biāo),求出熱像儀響應(yīng)系數(shù)A、B;
(2)利用熱像儀測量被測物體輻射灰度Hs,根據(jù)熱像儀響應(yīng)公式計(jì)算目標(biāo)輻射亮度Ls:
其中,λ1,λ2為紅外熱像儀響應(yīng)起始波長和終止波長;c1=3.74×104W·μm4·cm-2;c2=1.44×104μm·K。
2.2 誤差來源分析
傳統(tǒng)工程方法測量物體表面發(fā)射率主要誤差在于沒有考慮環(huán)境輻射對(duì)測量各因素的影響,具體可分為以下四個(gè)方面。
2.2.1 對(duì)響應(yīng)系數(shù)A、B的影響
測量試驗(yàn)中,一般把黑體與被測設(shè)備并排擺放以保證熱像儀對(duì)黑體和被測設(shè)備同時(shí)清晰成像,此時(shí)熱像儀接收的黑體輻射包括黑體自身輻射和黑體反射環(huán)境輻射兩部分,則熱像儀的響應(yīng)率A和背景電平B的計(jì)算公式變?yōu)椋ㄒ詢牲c(diǎn)定標(biāo)為例):
(3)測量被測物體表面溫度Ts,利用普朗克公式計(jì)算其表面發(fā)射率εs:
式中,下標(biāo)m、n,代表定標(biāo)黑體設(shè)定不同的溫度對(duì)應(yīng)的熱像儀測量灰度值H和黑體輻射亮度值L;εb為黑體發(fā)射率。定義Lrm=Erm/π為環(huán)境對(duì)黑體輻射照度的等效輻射亮度,Erm為環(huán)境輻射對(duì)黑體的輻射照度。
環(huán)境輻射對(duì)響應(yīng)系數(shù)A、B的影響表現(xiàn)為直接計(jì)算出的背景電平B′值中包括黑體反射環(huán)境輻射。當(dāng)環(huán)境中無強(qiáng)輻射源時(shí),由于黑體發(fā)射率很高,反射環(huán)境輻射對(duì)B值的影響較小,但數(shù)據(jù)處理時(shí)應(yīng)避免B值的使用。
2.2.2 對(duì)輻射亮度測量結(jié)果的影響
利用熱像儀測量被測物體表面輻射時(shí),熱像儀接收的能量同樣包括物體自身輻射和物體反射環(huán)境輻射兩部分,則公式(1)重新表示為:
其中,Ls為物體自身的紅外輻射亮度;Lrs為環(huán)境對(duì)物體的等效輻射亮度。
環(huán)境輻射對(duì)物體表面輻射亮度測量的影響表現(xiàn)為熱像儀接收輻射中始終同時(shí)包括物體自身輻射和物體反射環(huán)境輻射兩部分。由于物體表面發(fā)射率較低,以傳統(tǒng)工程方法測得的表面發(fā)射率與實(shí)際表面發(fā)射率誤差較大,圖1和圖2為長波紅外波段表面發(fā)射率測量誤差與實(shí)際表面發(fā)射率和環(huán)境等效輻射亮度之間的關(guān)系,圖中縱軸為傳統(tǒng)測量方法計(jì)算的發(fā)射率與實(shí)際發(fā)射率之比。
圖1 發(fā)射率測量誤差隨表面發(fā)射率變化關(guān)系Fig.1 Measurement error versus surface emissivity
圖2 發(fā)射率測量誤差隨環(huán)境等效輻射亮度變化關(guān)系Fig.2 Measurement error versus environmental equivalent Radiance
圖1 為當(dāng)物體表面溫度為293 K、環(huán)境等效輻射亮度為25 W/(m2·sr)時(shí),傳統(tǒng)工程測量方法計(jì)算的發(fā)射率誤差與表面實(shí)際發(fā)射率的關(guān)系,從圖中可以看出當(dāng)表面發(fā)射率從0.3變化至0.9時(shí),傳統(tǒng)工程測量方法計(jì)算的發(fā)射率誤差從約170%變化至約10%,表面發(fā)射率越低,傳統(tǒng)工程測量方法的誤差越大。圖2為當(dāng)物體表面溫度為293 K、表面發(fā)射率為0.8時(shí),傳統(tǒng)工程測量方法計(jì)算的發(fā)射率誤差與環(huán)境等效輻射亮度的關(guān)系,從圖中可以看出當(dāng)環(huán)境等效輻射亮度從10W/(m2·sr)變化至40 W/(m2·sr)時(shí),傳統(tǒng)工程測量方法計(jì)算的發(fā)射率誤差從約8%變化至約29%,環(huán)境等效輻射亮度越小,傳統(tǒng)工程測量方法的誤差越小。
2.2.3 漫反射與鏡面反射的影響
工程上用熱像儀測量物體表面發(fā)射率時(shí)一般假設(shè)表面為朗伯體,但實(shí)際表面的反射同時(shí)包括漫反射和鏡面反射兩種,對(duì)不透明表面:
式中,ρ為表面反射率;ε為表面發(fā)射率;ρ1為漫反射率;ρ2為鏡面反射率。
漫反射與鏡面反射的差異可用圖3表示,環(huán)境整體對(duì)表面的等效輻射亮度為L0,某一方向環(huán)境輻射亮度為L0d,則對(duì)應(yīng)方向熱像儀接收的反射環(huán)境輻射為:
圖3 漫反射與鏡面反射示意圖Fig.3 Schematic diagram of diffuse and specular reflection
由于大型設(shè)備一般處于露天等較為復(fù)雜且不可控的紅外環(huán)境中,環(huán)境輻射的空間差異性會(huì)通過表面漫反射和鏡面反射的區(qū)別對(duì)發(fā)射率測量產(chǎn)生影響,實(shí)驗(yàn)布局的不同會(huì)得到不同的測量結(jié)果。
2.2.4 熱像儀響應(yīng)漂移的影響
實(shí)驗(yàn)過程中,由于熱像儀自身的特性其響應(yīng)特性會(huì)隨時(shí)間而漂移,某型氧化釩非制冷型熱像儀開機(jī)35min后對(duì)同一溫度的黑體測得灰度變化如圖4所示。
圖4 熱像儀響應(yīng)漂移測試數(shù)據(jù)Fig.4 Response drift of thermal infrared imager test data
從圖中可以看出,4 min內(nèi)熱像儀對(duì)同一溫度黑體響應(yīng)灰度從5320變化至5275,灰度變化達(dá)45,熱像儀響應(yīng)率A和背景電平B也發(fā)生改變。熱像儀響應(yīng)特性的變化會(huì)在測量數(shù)據(jù)處理時(shí)帶來誤差。
傳統(tǒng)的利用紅外熱像儀測量表面發(fā)射率時(shí),由于熱像儀接收輻射中始終同時(shí)包括物體自身輻射和反射環(huán)境輻射且不能區(qū)分測量兩者的大小,導(dǎo)致測量結(jié)果存在誤差。若利用兩發(fā)射率已知的設(shè)備,可測量環(huán)境對(duì)被測物體的等效輻射亮度,并通過改進(jìn)實(shí)驗(yàn)方法,消除2.2.3和2.2.4節(jié)中誤差來源,提高測量精度。
3.1 鏡面反射近似
當(dāng)環(huán)境輻射的空間分布較為均勻時(shí),可近似認(rèn)為環(huán)境整體等效輻射亮度與任一方向的環(huán)境輻亮度相等,公式(6)可表示為:
一般的,光學(xué)暗室內(nèi)因墻壁地面均涂有消光漆其發(fā)射率相同,若條件限制也可選擇無窗的庫房,當(dāng)室內(nèi)溫度穩(wěn)定時(shí)可認(rèn)為各方向輻射亮度相同。
3.2 環(huán)境等效輻射亮度測量
利用發(fā)射率已知的面源黑體(發(fā)射率為εb)和標(biāo)準(zhǔn)漫反射板(發(fā)射率為εd),將兩者緊密擺放并使其表面法線方向一致,則可認(rèn)為環(huán)境對(duì)兩者的等效輻射亮度相同,為Le。對(duì)面源黑體,熱像儀響應(yīng)灰度為:
對(duì)標(biāo)準(zhǔn)漫反射板,熱像儀響應(yīng)灰度為:
由式(8)和式(9),則環(huán)境等效輻射亮度計(jì)算表達(dá)式為:
其中,Lb和Ld為面源黑體和標(biāo)準(zhǔn)漫反射板輻射亮度,根據(jù)各自表面溫度及發(fā)射率計(jì)算。
3.3 設(shè)備表面發(fā)射率測量方法一
若將設(shè)備與面源黑體、標(biāo)準(zhǔn)漫反射板緊密擺放且使其表面法線方向一致,則認(rèn)為環(huán)境對(duì)三者等效輻射亮度相同。熱像儀對(duì)被測設(shè)備響應(yīng)灰度為:
式中,εs為設(shè)備發(fā)射率,為待求參數(shù)。由于設(shè)備輻射亮度Ls中同樣含有待求參數(shù)εs,因此以設(shè)備輻射出射度表示:
式中,Es為設(shè)備表面輻射出射度,可根據(jù)普朗克公式計(jì)算。
由式(8)-式(11),并帶入式(12),可得設(shè)備表面發(fā)射率計(jì)算公式:
式中,環(huán)境等效輻射亮度Le根據(jù)公式(10)計(jì)算。
3.4 設(shè)備表面發(fā)射率測量方法二
由于熱像儀響應(yīng)率A和背景電平B隨時(shí)間變化引起測量誤差,結(jié)合公式(10)和公式(13)表達(dá)形式,當(dāng)調(diào)整面源黑體溫度使其灰度分別與標(biāo)準(zhǔn)漫反射板、設(shè)備灰度相同時(shí),可在表達(dá)式中約去響應(yīng)率A。類似公式(12)形式,以輻射出射度代替輻射亮度,式(10)、(13)可表示為:
式中,Ed為標(biāo)準(zhǔn)漫反射板輻射出射度;Es為設(shè)備輻射出射度;Eb1為當(dāng)面源黑體與標(biāo)準(zhǔn)漫反射板灰度相同時(shí)面源黑體的輻射出射度;Eb2為當(dāng)面源黑體與設(shè)備灰度相同是黑體輻射出射度。
由公式(15)中可知該方法不需要進(jìn)行熱像儀響應(yīng)標(biāo)定,可避免熱像儀響應(yīng)線性近似帶來的誤差[8]。
4.1 實(shí)驗(yàn)設(shè)備與布局
實(shí)驗(yàn)用熱像儀為氧化釩非制冷型熱像儀,像素?cái)?shù)為320×240,響應(yīng)波段為7.5~13.5μm,最小焦距為5m;面源黑體為SR-800型,控溫精度優(yōu)于0.008 K,溫度均勻性優(yōu)于0.015 K,表面發(fā)射率為0.97;標(biāo)準(zhǔn)漫反射板在7.5~13.5μm波段的平均發(fā)射率為0.68。
實(shí)驗(yàn)布局如圖5所示。
圖5 測量實(shí)驗(yàn)布局圖Fig.5 Sketch map of experiment layout
測量試驗(yàn)選擇在紅外環(huán)境簡單的庫房內(nèi)進(jìn)行,窗戶做遮光處理,室內(nèi)各處溫度與發(fā)射率基本一致,室內(nèi)各方向輻射亮度近似相等;面源黑體、標(biāo)準(zhǔn)漫反射板與設(shè)備緊密擺放,并使其表面法線方向相同;熱像儀距設(shè)備9 m,調(diào)整其焦距使其對(duì)設(shè)備清晰成像;實(shí)驗(yàn)中使用測溫儀器測量環(huán)境溫度、標(biāo)準(zhǔn)漫反射溫度和設(shè)備表面溫度。
4.2 實(shí)驗(yàn)方法
實(shí)驗(yàn)過程中發(fā)現(xiàn)熱像儀響應(yīng)漂移明顯,因此采用3.4節(jié)的方法進(jìn)行測量實(shí)驗(yàn),實(shí)驗(yàn)步驟如下:
(1)熱像儀開機(jī)并穩(wěn)定約10 min,調(diào)整熱像儀焦距使其對(duì)設(shè)備、面源黑體和標(biāo)準(zhǔn)漫反射板清晰成像,調(diào)整完畢后不再調(diào)整熱像儀焦距;
(2)利用熱像儀的中心像素點(diǎn)P測量標(biāo)準(zhǔn)漫反射板灰度H1;
(3)調(diào)整熱像儀角度使像素點(diǎn)P對(duì)準(zhǔn)面源黑體,調(diào)整面源黑體溫度使P點(diǎn)灰度為H1;
(4)重復(fù)步驟(2)、(3)使二者灰度完全相同,記錄此時(shí)面源黑體溫度與標(biāo)準(zhǔn)漫反射板溫度;
(5)利用熱像儀像素點(diǎn)P測量設(shè)備灰度H2;
(6)調(diào)整熱像儀角度使像素點(diǎn)P對(duì)準(zhǔn)面源黑體,調(diào)整面源黑體溫度使P點(diǎn)灰度為H2;
(7)重復(fù)步驟(5)、(6)使二者灰度完全相同,記錄此時(shí)面源黑體溫度與設(shè)備溫度;
(8)重復(fù)(5)~(7)測量設(shè)備其他部位。
(9)實(shí)驗(yàn)中實(shí)時(shí)監(jiān)測并記錄環(huán)境溫度。
進(jìn)行了兩組實(shí)驗(yàn),測量了設(shè)備表面4個(gè)不同位置的發(fā)射率,測量圖像如圖6所示,測量數(shù)據(jù)及處理結(jié)果如表1所示,傳統(tǒng)發(fā)射率測量值指以傳統(tǒng)測量方法處理實(shí)驗(yàn)數(shù)據(jù)(公式(2))得到的發(fā)射率,改進(jìn)發(fā)射率測量值指以改進(jìn)測量方法處理實(shí)驗(yàn)數(shù)據(jù)(公式(15))得到的發(fā)射率。實(shí)驗(yàn)中熱像儀距設(shè)備9m,大氣透過率近似為1,大氣路徑輻射亮度近似為0。
圖6 實(shí)驗(yàn)測量圖像Fig.6 Measuring Images
表面位置1與表面位置2涂有相同的漆層,因此其發(fā)射率測量結(jié)果相同;表面位置3由于漆層被破壞,發(fā)射率要高于位置1、2;表面位置4由于其結(jié)構(gòu)類似于腔式黑體,因此其發(fā)射率較高,與面源黑體的發(fā)射率相同。
由表1可知,由于環(huán)境輻射影響,以傳統(tǒng)測量方法處理實(shí)驗(yàn)數(shù)據(jù),得到的發(fā)射率要普遍高于改進(jìn)測量方法得到的發(fā)射率,發(fā)射率越低,相對(duì)誤差越大,變化規(guī)律與2.2.2節(jié)中理論分析結(jié)果相一致,驗(yàn)證了測量結(jié)果的準(zhǔn)確性。
傳統(tǒng)測量方法和改進(jìn)測量方法處理得到的表面位置1和表面位置2的發(fā)射率大小關(guān)系恰好相反,其原因在于實(shí)驗(yàn)環(huán)境等效輻射亮度達(dá)43.07 W/(m2·sr),高于設(shè)備的輻射亮度,因此發(fā)射率越低,設(shè)備反射環(huán)境輻射越大,熱像儀接收輻射量越高,以傳統(tǒng)測量方法測得的發(fā)射率也越大,因此其相對(duì)誤差也越大。
表1 實(shí)驗(yàn)測量數(shù)據(jù)及處理結(jié)果Tab.1 Testing data and processing results
本文主要考慮環(huán)境輻射對(duì)發(fā)射率傳統(tǒng)工程測量方法的影響,提出環(huán)境等效輻射亮度概念,理論分析了傳統(tǒng)工程測量方法誤差與發(fā)射率和環(huán)境等效輻射照度的大小和變化關(guān)系。在此基礎(chǔ)上,提出了改進(jìn)的發(fā)射率工程測量方法,并理論推導(dǎo)了兩種測量方法的發(fā)射率測量公式。設(shè)計(jì)并組織了發(fā)射率測量實(shí)驗(yàn),實(shí)驗(yàn)結(jié)果表明,傳統(tǒng)測量方法與改進(jìn)測量方法得到的發(fā)射率相對(duì)誤差可達(dá)20.37%,相對(duì)誤差變化規(guī)律與理論分析結(jié)果一致,驗(yàn)證了該方法的準(zhǔn)確性。
該方法相對(duì)于傳統(tǒng)工程測量方法而言提高了發(fā)射率測量精度,尤其對(duì)發(fā)射率偏低的設(shè)備,測量精度提升較為明顯;實(shí)驗(yàn)所需配套設(shè)備易于獲取,實(shí)驗(yàn)步驟簡單,具有很強(qiáng)的可操作性,因此具有良好的實(shí)際應(yīng)用價(jià)值。
該方法為基于熱像儀的表面波段發(fā)射率測量方法,在測量波段內(nèi)對(duì)光譜值進(jìn)行了平均化處理,包括熱像儀的光譜響應(yīng)特性、標(biāo)準(zhǔn)漫反射板的光譜發(fā)射率等,因此測量結(jié)果會(huì)存在誤差,此外,實(shí)驗(yàn)中近似環(huán)境各空間方向輻射亮度相同,也會(huì)對(duì)測量結(jié)果造成一定誤差。
[1] Zhang Jianqi,F(xiàn)ang Xiaoping.Infrared physics[M].Xi’an:Xidian University Press,2004.(in Chinese)張建奇,方小平.紅外物理[M].西安:西安電子科技大學(xué)出版社,2004.
[2] Dai Jingmin,Song Yang,Wang Zongwei.Review of spectral emissivitymeasurement[J].Infrared and Laser Engineering,2009,38(4):710-715.(in Chinese)戴景民,宋揚(yáng),王宗偉.光譜發(fā)射率測量技術(shù)[J].紅外與激光工程,2009,38(4):710-715.
[3] Luo Xiaochun,Sun Jiyin.Method for measuring obiect surface emissivitv[J].Automatic Measurement andd Control,2007,26(8):72-74.(in Chinese)羅曉春,孫繼銀.物體表面發(fā)射率的測量方法[J].自動(dòng)測量與控制,2007,26(8):72-74.
[4] Wang Xianbing,Yang Shizhi,Qiao Yanli,et al.Field measurementmethods for thermal infrared spectral emissivity of terrestrial surfacemeterial[J].Journal of Atomspheric and Environmental Optics,2006,1(2):105-111.(in Chinese)王先兵,楊世植,喬延利,等.地面目標(biāo)熱紅外光譜發(fā)射率的野外測量方法研究[J].大氣與環(huán)境光學(xué)學(xué)報(bào),2006,1(2):105-111.
[5] Xie Pinhua,Wei Qinnong,Liu Jianguo,et al.A method of measuring infrared emissivity of paint coatings by bidirectional reflectance[J].Laser&Infrared,1999,29(1):42-45.(in Chinese)謝品華,魏慶農(nóng),劉建國,等.用雙向反射法測量漆層的紅外發(fā)射率[J].激光與紅外,1999,29(1):42-45.
[6] Huang Longxiang,Shen Xiangheng,Song Jiangtao.Measure target wideband emissivity with thermal imager[J]. Laser&Infrared,2009,39(2):159-161.(in Chinese)黃龍祥,沈湘衡,宋江濤.基于熱像儀的物體波段發(fā)射率的測量[J].激光與紅外,2009,39(2):159-161.
[7] Yang Li,Kou Wei,Liu Huikai,et al.Surface emissivity measurement and error analysis using infrared thermography[J].Larse&Infrared,2002,32(1):43-45.(in Chinese)楊立,寇蔚,劉慧開,等.熱像儀測量物體表面輻射率及誤差分析[J].激光與紅外,2002,32(1):43-45.
[8] Wu Junhui,Zhu Zhen,Li Dongqiang,et al.Error calculation of IR radiance measurement with thermal cameras caused by data processing based on waveband[J].Infrared Technology,2005,27(1):23-28.(in Chinese)吳軍輝,朱震,李東強(qiáng),等.基于波段處理而引起的熱像儀輻射測量誤差計(jì)算[J].紅外技術(shù),2005,27(1):23-28.
M ethod for surface em issivity measurement
LIU Lian-wei,YANG Miao-miao,F(xiàn)AN Hong-jie,XU Zhen-ling,YAO Mei
(Unit63892 of PLA,Luoyang 471003,China)
It is difficult to accurately measure the large-scale equipment surface emissivity,because it is difficult to control the temperature of the equipment.A method of emissivity measurement is proposed based on analyzing the main source of emissivitymeasurementerrors.The radiation that reaches the detector is analyzed and derived,environmental IR radiationmeasurements are achieved by using two apparatuswith known emissivity,thereby the environmental radiation on the emissivity measurements effects is eliminated.Experimental method is designed and emissivity measurements experiment ismade,analysis results show that,compared with themethod ofno considering the environmental radiation,there error is reduced by 20.37%.Themethod has feasibility and accuracy for engineeringmeasurement of surface emissivity.
emissivity;thermal infrared imager;environmental IR radiation
TN21
A
10.3969/j.issn.1001-5078.2014.02.0
1001-5078(2014)02-0152-06
總裝預(yù)研基金(No.51303040203)資助。
劉連偉(1985-),男,助理工程師,本科,主要從事紅外對(duì)抗仿真實(shí)驗(yàn)技術(shù)方面的研究。E-mail:yyoohh@126.com
2013-06-28