馬舒慶徐振飛毛節(jié)泰劉達(dá)新張春波楊 玲甄小瓊

1）（中國氣象局氣象探測中心，北京100081）2）（成都信息工程學(xué)院電子工程學(xué)院，成都610225）

3）（北京大學(xué)物理學(xué)院，北京100080）4）（北京市觀象臺(tái)，北京100089）

以黑體為目標(biāo)的能見度參考標(biāo)準(zhǔn)試驗(yàn)研究

馬舒慶1）*徐振飛2）毛節(jié)泰3）劉達(dá)新1）張春波4）楊 玲2）甄小瓊2）

1）（中國氣象局氣象探測中心，北京100081）2）（成都信息工程學(xué)院電子工程學(xué)院，成都610225）

3）（北京大學(xué)物理學(xué)院，北京100080）4）（北京市觀象臺(tái)，北京100089）

黑體為目標(biāo)的能見度測量是從科西米德定律出發(fā)，通過工業(yè)相機(jī)對黑體目標(biāo)物進(jìn)行拍攝，建立黑體和背景天空的數(shù)學(xué)模型，求得能見度。分析表明：當(dāng)能見度小于30 km，黑體黑度和CCD（change-coupled device）工業(yè)相機(jī)的面非均勻性對該方法測得的能見度造成的誤差約為3.7%。該方法與前向散射能見度儀測量進(jìn)行對比，能見度變化趨勢上具有一致性；但在高能見度時(shí)，黑體能見度測量值大于前向散射能見度測量值，低能見度時(shí)，黑體測量值小于前向散射能見度測量值。

黑體；黑度；能見度；面均勻性

引 言

氣象能見度［1-2］是與日常生活和生產(chǎn)密切相關(guān)的氣象要素［3-4］。如何實(shí)現(xiàn)氣象能見度自動(dòng)化觀測，國內(nèi)外做了大量工作［5-6］，開發(fā)了多種氣象能見度觀測儀［7-8］，能見度觀測儀器的研發(fā)和試驗(yàn)考核都需要能見度參考標(biāo)準(zhǔn)來驗(yàn)證。1988年10月—1989年5月WMO（世界氣象組織）在英國舉行了第1次國際能見度比對試驗(yàn)［9］，來自18個(gè)廠家的25個(gè)能見度儀器參加了此次試驗(yàn)，其中包括1臺(tái)后向散射能見度儀、多臺(tái)前向散射能見度儀和多種基線的透射式能見度儀。從中選擇一部分儀器的測量結(jié)果的中間值作為比對的參考值。選擇用于計(jì)算參考值的能見度儀時(shí)，需要考慮其量程范圍以及其觀測性能。與當(dāng)?shù)氐娜斯び^測進(jìn)行比對后給出每臺(tái)能見度儀器的不確定度。英國BIRAL公司采用FAA（美國航空協(xié)會(huì)）認(rèn)可的標(biāo)準(zhǔn)透射式能見度儀作為參考［10］，對VPF-710，VPF-730能見度儀器進(jìn)行為期兩年的大氣消光系數(shù)對比試驗(yàn)，其中VPF-710采用前向散射式技術(shù)，而VPF-730采用后向散射式技術(shù)。美國Environment Sensors公司的Crosby［11］指出，AES（加拿大大氣環(huán)境局）用芬蘭Vaisala的大氣透射表來對比英國BIRAL公司的能見度傳感器。德國GMS（German Meteorological Service）的Waas［12］用大氣透射表與前向散射型能見度儀進(jìn)行比對試驗(yàn)。毛節(jié)泰等［13］利用望遠(yuǎn)光度計(jì)，根據(jù)能見度定義測量大氣能見度。謝興生等［14］提出數(shù)字?jǐn)z像法測量氣象能見度原理和方法。呂偉濤等［15］對于長方體結(jié)構(gòu)的黑體進(jìn)行了理論分析和計(jì)算，給出了黑體近軸向的半球方向反射率與其特性參數(shù)之間的關(guān)系。王緬等［16］提出了利用透射能見度儀標(biāo)定和漫反射板校準(zhǔn)前向散射式能見度儀的方法。呂偉濤等［17］分析了雙亮度差方法測量白天氣象能見度時(shí)的誤差。周建榮等［18］也將大氣透射表觀測結(jié)果與人工觀測結(jié)果進(jìn)行了比對分析試驗(yàn)。程紹榮等［19］提出采用漸變率光衰減鏡頭的方法對大氣能見度儀進(jìn)行標(biāo)定。

從國內(nèi)外開展的能見度對比試驗(yàn)可以看出，能見度對比試驗(yàn)的參考標(biāo)準(zhǔn)主要是透射式能見度儀、多種儀器的分析計(jì)算值以及人工觀測。透射式能見度儀量程較小，多種儀器分析計(jì)算值和人工觀測值的不確定性均為建立能見度參考標(biāo)準(zhǔn)的難點(diǎn)所在。本文提出了從能見度定義出發(fā)的黑體能見度參考標(biāo)準(zhǔn)，該標(biāo)準(zhǔn)具有客觀和準(zhǔn)確的優(yōu)點(diǎn)。

1 以黑體為目標(biāo)的能見度參考標(biāo)準(zhǔn)構(gòu)成

圖1給出了黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)。采用工業(yè)相機(jī)獲取黑體和天空的圖像。工業(yè)相機(jī)鏡頭朝南，在相機(jī)南面680 m處選擇窗口朝北的房屋，將其窗口內(nèi)部做消光處理，使其作為黑體。相機(jī)每6 s拍攝1次黑體的圖片。取圖像中天空灰度值和圖像中黑體灰度值，換算成亮度值D∞和D′L。

圖1 以黑體為目標(biāo)的能見度參考標(biāo)準(zhǔn)示意圖Fig.1 Schematic diagram of visibility reference standard blackbody targets

當(dāng)氣柱光學(xué)特性均勻，那么

其中，D∞是天空視亮度，DL是黑體視亮度，D是氣柱亮度，DL0是黑體亮度，k是消光系數(shù)，L是黑體目標(biāo)距離，s是黑體黑度。按氣象能見度的概念，文獻(xiàn)［2］中給出水平方向黑體及其周邊水平天空構(gòu)成的視亮度對比小于閾值0.05，其相應(yīng)的距離為氣象能見距離，即氣象能見度

假定工業(yè)相機(jī)拍攝圖片中天空處為基準(zhǔn)，那么DL=D′L×A，其中D′L為直接用工業(yè)相機(jī)灰度值轉(zhuǎn)換得到的黑體視亮度，A=1-u是考慮了工業(yè)相機(jī)非均勻性的訂正參數(shù)，u是非均勻性參數(shù)。將DL=D′L×A帶入式（6），得

對R進(jìn)行微分得到

其中，

黑體黑度帶來的能見度誤差和CCD面非均勻性帶來的能見度誤差可能會(huì)符號(hào)相反，互相抵消一部分，但最大誤差是兩者疊加，因此，在能見度不超過30 km時(shí)最大相對誤差約為3.7%。

2 黑 體

2.1 黑體結(jié)構(gòu)

由于在觀測場外架設(shè)一個(gè)視角（相對在觀測場中的相機(jī)）足夠大的物體非常困難，因此租用在相機(jī)視場內(nèi)的房屋。房屋只留一個(gè)400 mm×400 mm窗口，內(nèi)部空間為2500 mm×2500 mm×2600 mm，采用能吸收光線的黑布覆蓋室內(nèi)6面墻及天花板和地面。從400 mm×400 mm窗口進(jìn)入的光線被黑布吸收和散射，經(jīng)過多次散射和吸收，使從窗口射出的光線極其微弱。

2.2 白卡測量

使用FieldSpec Pro型光譜儀測得白卡反射率為0.7332，測量方法如下：①將標(biāo)準(zhǔn)白卡和待測白卡如圖2放置，使之正面向陽光；②標(biāo)準(zhǔn)卡測量（探頭垂直向下，離標(biāo)準(zhǔn)卡約25 cm）2次，如此測量待測白卡20次，再測量標(biāo)準(zhǔn)卡2次，由于邊上有污點(diǎn)，所以探測面積始終保持在中心1/2面積內(nèi)；③對試驗(yàn)數(shù)據(jù)進(jìn)行計(jì)算分析，按式（11）得到待測卡的反射率。

式（11）中，λ是波長，Rλ是目標(biāo)反射率，n是目標(biāo)測量次數(shù)，F(xiàn)iλ是來自目標(biāo)的反射輻射，F(xiàn)0λ是來自參考板的平均反射輻射，ρλ是測量時(shí)刻參考板的反射比因子，它由實(shí)驗(yàn)室測量的參考板BRF雙向反射比因子插值獲得。

圖2 白卡反射率測量Fig.2 White card reflectivity measurement

2.3 黑度和黑度測量

理論上可見光的黑體是一個(gè)亮度為零的物體，但實(shí)際加工的黑體不可能亮度為零，因此，定義加工的黑體的黑度為出射光與入射光之比。假定黑體的入射光光量為L0，出射光光量為L1，黑體的黑度為

對房屋黑體的黑度測量方法是在窗口下沿放一塊標(biāo)準(zhǔn)白卡，相機(jī)鏡頭距離窗口600 mm，相機(jī)可同時(shí)拍到白卡和窗口（圖3），從相機(jī)拍攝到的圖片讀取白卡和窗口的亮度數(shù)據(jù)，將窗口的亮度作為出射光亮度，白卡的亮度除以0.7332作為入射光亮度，兩者的比值為黑度，實(shí)測房屋的黑度數(shù)據(jù)為0.0018。

標(biāo)定方法使用近焦廣角段拍攝，作為實(shí)際情況中長焦小角度拍攝的衡量標(biāo)準(zhǔn)。

圖3 黑體的黑度測量示意圖Fig.3 Schematic diagram of blackbody target measurement

3 相機(jī)標(biāo)定

在中國科學(xué)院安徽光學(xué)精密機(jī)械研究所采用積分球進(jìn)行工業(yè)相機(jī)非均勻性和線性標(biāo)定。積分球主要性能指標(biāo)：①出光口輸出面非均勻性為0.19%，②角度非均勻性為水平±50°最大偏差0.0154，垂直最大偏差0.0139，③積分球非穩(wěn)定性為0.80%（9 h）。

3.1 工業(yè)圖像非均勻性標(biāo)定

將相機(jī)固定在積分球前的支架上，其光軸垂直于出光口面，相機(jī)距積分球出光口30 cm（圖4）。

圖4 CCD拍攝平臺(tái)Fig.4 The CCD camera platform

采樣50次取平均值。非均勻性為

式（13）中，u為非均勻性，ˉV為各有效點(diǎn)測量的平均值，Si為各有效點(diǎn)測量值的標(biāo)準(zhǔn)偏差。計(jì)算標(biāo)準(zhǔn)偏差

式（14）中，x為樣本值，ˉx為樣本平均值，n為樣本數(shù)量。

圖5a為相機(jī)帶長焦鏡頭拍攝照片，圖5b為分析圖，可以看出相機(jī)拍攝照片中部值較大，向四周遞減。相機(jī)中部1/4面積，面非均勻性為0.24%。

3.2 工業(yè)相機(jī)線性標(biāo)定

將積分球的光源在24盞、16盞、8盞、4盞的情況下，光圈固定，快門也固定為250μs，得到灰度與亮度的對應(yīng)關(guān)系數(shù)據(jù)（表1）和數(shù)據(jù)曲線（圖6）。

圖5 帶長焦鏡頭CCD拍攝圖（a）及三維示意圖（b）Fig.5 CCD with telephoto lens shooting figure（a）and 3D sketch（b）

表1 積分球不同燈盞亮度灰度數(shù)據(jù)Table 1 Integral ball lamp brightness of different gray scale data

圖6 亮度-灰度轉(zhuǎn)換曲線Fig.6 Brightness gray transformation curve

4 試驗(yàn)結(jié)果

由如圖1所示的黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)與前向散射能見度儀進(jìn)行同步對比，前向散射能見度儀是凱邁公司CJY-2C前向散射式能見度儀，當(dāng)能見度在10 km以下，CJY-2C相對誤差小于10%，CJY-2C與黑體能見度測量系統(tǒng)的工業(yè)相機(jī)相距15 m。

圖7給出了2011年11月6日兩種能見度測量設(shè)備測量數(shù)據(jù)的時(shí)間曲線（時(shí)間為北京時(shí)，下同）。黑體測量值小于前向散射能見度測量值，平均相對偏差為-0.26，標(biāo)準(zhǔn)偏差為0.04。

圖8給出了2011年11月8日兩種能見度測量設(shè)備測量數(shù)據(jù)的時(shí)間曲線。15：02以前黑體測量值大于前散能見度測量值，平均相對偏差為0.07，標(biāo)準(zhǔn)偏差為0.09。15：02以后黑體測量值小于前散能見度測量值，平均相對偏差為-0.16，標(biāo)準(zhǔn)偏差為0.04。

圖7 2011年11月6日黑體能見度數(shù)據(jù)Fig.7 Blackbody visibility data on 6 November 2011

圖8 2011年11月8日黑體能見度數(shù)據(jù)Fig.8 Blackbody visibility data on 8 November 2011

由圖7和圖8可以看出，黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)和前向散射能見度儀測量能見度，在能見度變化的趨勢上具有一致性，但在高能見度（大于8 km）時(shí)，黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)測量值大于前向散射能見度測量值，低能見度（小于3 km）時(shí)黑體測量值小于前向散射能見度測量值。這種特征可能是因?yàn)榍跋蛏⑸淠芤姸葍x原理上不能測量大氣中氣溶膠的吸收作用，當(dāng)能見度高時(shí)，氣溶膠較少，對光的吸收作用可以忽略，前向散射能見度儀測量比較準(zhǔn)確。當(dāng)能見度低時(shí)，氣溶膠較多，對光的吸收作用較大，而前向散射能見度儀只測量散射衰減，不能測量吸收衰減，測量能見度值可能偏高，導(dǎo)致比黑體能見度測量值大。這只是一個(gè)初步分析，證實(shí)該結(jié)論還需同步觀測氣溶膠溶度和光學(xué)特性進(jìn)行深入研究。

5 小 結(jié)

1）以黑體為目標(biāo)黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)由圖像傳感器和黑體目標(biāo)組成。結(jié)構(gòu)簡單，自動(dòng)化運(yùn)行。

2）黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)的測量誤差來自黑體的黑度和CCD工業(yè)相機(jī)的非均勻性。分析可知，能見度在30 km以下，黑體的黑度和CCD工業(yè)相機(jī)的非均勻性帶來的誤差約為3.7%。

3）開展黑體能見度測量系統(tǒng)與前向散射能見度儀對比觀測試驗(yàn)，統(tǒng)計(jì)了兩者變化趨勢，相關(guān)系數(shù)為0.3～0.6。黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)與前向散射能見度儀測量在能見度變化趨勢上具有一致性，但高能見度時(shí)（大于8 km），黑體能見度參考標(biāo)準(zhǔn)系統(tǒng)測量值大于前向散射能見度儀測值，低能見度時(shí)（小于3 km），黑體能見度測量值小于前向散射能見度儀測值。

［1］ 中國氣象局.地面氣象觀測規(guī)范.北京：氣象出版社，2003：21-27.

［2］ WMO.Guide to Meteorological Instrumentsand Methods of Observation（No.8）.2010：230-233.

［3］ 黃健，吳兌，黃敏輝，等.1954—2004年珠江三角洲大氣能見度變化趨勢.應(yīng)用氣象學(xué)報(bào)，2008，19（1）：61-70.

［4］ 袁成松，卞光輝，馮民學(xué).高速公路上低能見度的監(jiān)測與預(yù)報(bào).氣象，2003，29（11）：36-40.

［5］ 林云，孫向明，張小麗，等.深圳市大氣能見度與細(xì)粒子濃度統(tǒng)計(jì)模型.應(yīng)用氣象學(xué)報(bào)，2009，20（2）：252-256.

［6］ 劉西川，高太長，劉磊，等.降水現(xiàn)象對大氣消光系數(shù)和能見度的影響.應(yīng)用氣象學(xué)報(bào)，2010，21（4）：433-441.

［7］ 曾書兒，王改利.能見度的觀測及其儀器.應(yīng)用氣象學(xué)報(bào)，1999，10（2）：207-212.

［8］ 馬舒慶，吳可軍，陳冬冬.天氣現(xiàn)象自動(dòng)化觀測系統(tǒng)設(shè)計(jì).氣象，2011，37（9）：1166-1172.

［9］ WMO.CIMO International Visibility Intercomparison.1989.

［10］ 馬遠(yuǎn).英國BIRAL HSS高精度能見度/天氣現(xiàn)象儀.氣象水文海洋儀器，2005，2（7）：47-50.

［11］ Crosby D.Visibility Sensor Accuracy：What’s Realistic？The 12th Symposium on Meteorological Observations and Instrumentation 2003 American Meteorological Society Annual Meeting，2003.

［12］ Waas S.Field Test of Forword Scatter Visibility Sensors at German Airports.TECH-2008，2008.

［13］ 毛節(jié)泰，李建國.氣象能見度與望遠(yuǎn)光度計(jì).大氣科學(xué)，1984，2（6）：170-177.

［14］ 謝興生，陶善昌，周秀驥.數(shù)字?jǐn)z像法測量氣象能見度.科學(xué)通報(bào)，1999，44（1）：97-100.

［15］ 呂偉濤，陶善昌.數(shù)字?jǐn)z像能見度觀測系統(tǒng)中實(shí)用黑體技術(shù)的應(yīng)用.應(yīng)用氣象學(xué)報(bào)，2003，12（6）：692-699.

［16］ 王緬，劉文清，陸亦懷，等.氣溶膠前向散射大氣能見度測量系統(tǒng)傳遞系數(shù)的標(biāo)定及校準(zhǔn)方法.光學(xué)技術(shù)，2008，34（3）：334-337.

［17］ 呂偉濤，陶善昌，譚涌波.雙亮度差方法測量白天氣象能見度時(shí)的誤差分析.應(yīng)用氣象學(xué)報(bào)，2005，16（5）：619-628.

［18］ 周建榮，楊銳.能見度自動(dòng)觀測系統(tǒng)與人工觀測的對比分析.中國科技論文在線，2008，1（13）：1436-1439.

［19］ 程紹榮，呂軍，楊寧，等.大氣能見度儀的定標(biāo)系統(tǒng)及其定標(biāo)方法：中國，CN101661155A.2010-03-03.

Experimental Research on Visibility Reference Standard for Blackbody Targets

Ma Shuqing1）Xu Zhenfei2）Mao Jietai3）Liu Daxin1）Zhang Chunbo4）Yang Ling2）Zhen Xiaoqiong2）

1）（Meteorological Observation Center of CMA，Beijing 100081）
2）（Department of Electronic Engineering，Chengdu University of Information Technology，Chengdu 610225）
3）（School of Physics，Peking University，Bejing 100871）4）（Beijing Weather Observatory，Beijing 100089）

The reference standard system of blackbody visibility is composed of industrial camera，blackbody and industrial personal computer.With the industrial camera lens facing south，shooting blackbody and sky background，the chosen house is 680 meters away from the camera and the window of the house faces north.A series of extinction processing are finished to ensure the house as a blackbody.The camera takes pictures of the blackbody every 6 seconds.The mathematical model of the reference standard system of blackbody is established according to the definition of meteorological visibility and the Cosimi De’s law，which shows the relationship between the meteorological visibility and sky/blackbody luminosity，and the relationship between the meteorological visibility and the blackbody blackness or the industrial camera image uniformity at the same time.Using the mathematical model，the error caused by the blackbody blackness and image uniformity of the camera is analyzed.When visibility is less than 30 km，the error caused by blackbody blackness and CCD industrial camera’s uniformity of the measured system is about 3.7%.Blackbody blackness is measured by photographic method，putting a standard white card under the blackbody window，and the camera lens is 600 mm away from the window，shooting the white card and window at the same time.The brightness of the window is taken as the light-emission luminance，and the brightness white card divided by 0.7332 is taken as the incident light intensity.The ratio between them is the brightness，and the blackness of measurements for the house is 0.0018.Integrating sphere is used for uniformity and linear calibration of the camera.Compared with the forward scattering visibility instrument measurement and the reference standard system of blackbody visibility，the value of the visibility is consistent，but is greater under high visibility condition，and the value of the blackbody visibility measurement is lower under low visibility condition.Perhaps the cause is that the principle of the forward scattering visibility instrument cannot measure atmospheric aerosol absorption.Under low visibility condition，the aerosol optical absorption effect is larger，and the forward scattering visibility instrument cannot measure the absorption attenuation and measurements may be higher.To confirm this conclusion，further solubility and synchronous observation of aerosol optical properties are needed.

blackbody；black degree；visibility；surface uniformity

馬舒慶，徐振飛，毛節(jié)泰，等.以黑體為目標(biāo)的能見度參考標(biāo)準(zhǔn)試驗(yàn)研究.應(yīng)用氣象學(xué)報(bào)，2014，25（2）：129-134.

2013-05-23收到，2013-10-28收到再改稿。

公益性行業(yè)（氣象）科研專項(xiàng)（GYHY201106047）

*email：msqaoc＠cma.gov.cn