賴建梁 張加春 陳 晞 林 鍇

(1.福建省德化九仙山氣象站，福建 泉州 362500; 2.泉州市氣象局，福建 泉州 362000)

1 概述

寶光(民間稱為“佛光”，以下以佛光表述)是一種特殊的大氣物理奇觀，此現(xiàn)象系陽(yáng)光從人體背后投射到云霧層，光線經(jīng)過(guò)云滴粒子的衍射和反射作用而形成圓形彩圈，人影又同時(shí)投射到彩圈云霧層上，猶如“佛祖再現(xiàn)”，有幸置身佛光中的人們自然認(rèn)為己有佛緣，乃至對(duì)其趨之若鶩。由于佛光的出現(xiàn)要求陽(yáng)光、地形(人體位置)和云層(佛光位置)呈由上而下直線分布[1]，能滿足要求的云層大多出現(xiàn)在高山地帶，因此佛光大多出現(xiàn)在山岳，由此而成為山岳的一項(xiàng)重要而獨(dú)特的氣象旅游景觀而倍受旅游者追捧，如峨眉山、黃山等[1-4]。

關(guān)于佛光的成因，國(guó)內(nèi)外學(xué)者提出過(guò)多種學(xué)說(shuō)，有“復(fù)雜散射”學(xué)說(shuō)[5]、“先反射、后衍射”學(xué)說(shuō)[6]及“先衍射、后反射”學(xué)說(shuō)[4，7]，學(xué)者賴比星[8]等進(jìn)一步制作了“先衍射、后反射”的佛光成像原理圖，此得到了中國(guó)氣象科學(xué)研究院張紀(jì)淮研究員和中國(guó)科學(xué)院大氣物理所黃美元研究員等有關(guān)專(zhuān)家的高度評(píng)價(jià)，但該學(xué)者在文獻(xiàn)[9]中又提到佛光乃為太陽(yáng)光先受云滴粒子反射、再繞射而成?？梢?jiàn)各種寶光成因原理至今尚無(wú)公認(rèn)理論。但從中也不難看出，衍射是各種佛光形成原理中的一個(gè)共同的物理過(guò)程，因此相關(guān)的衍射理論為本研究基礎(chǔ)之一。

由于佛光多出現(xiàn)在山麓云層之中，此為測(cè)量佛光圈大小及其所在位置的云滴粒子大小等物理參數(shù)帶來(lái)困難，因此，相關(guān)的觀測(cè)數(shù)據(jù)也鮮見(jiàn)。相關(guān)文獻(xiàn)對(duì)于形成佛光的衍射云滴大小大多只是做了粒子小的定性描述，并未進(jìn)行量上的闡述，如文獻(xiàn)[10]指出了佛光系太陽(yáng)光穿過(guò)特別小的云霧滴時(shí)發(fā)生“衍射”而擴(kuò)散出多圈光環(huán)所致。根據(jù)九仙山氣象站觀測(cè)人員多年的觀測(cè)經(jīng)驗(yàn)，即使在滿足“陽(yáng)光、人體和云層”呈由上而下直線分布的條件下，并非所有的云層均可產(chǎn)生佛光現(xiàn)象，佛光現(xiàn)象常出現(xiàn)在流動(dòng)的云霧層或云層上端，該云層具有稀薄的特點(diǎn)，可見(jiàn)佛光現(xiàn)象與云層的云滴粒子大小、濃度等方面存在著一定的關(guān)系。本文根據(jù)2020年1月17日14時(shí)50分一次實(shí)測(cè)的佛光相關(guān)數(shù)據(jù)，利用流體動(dòng)力學(xué)和衍射原理，分別計(jì)算所形成佛光的云滴粒子大小并進(jìn)行比對(duì)分析，由此對(duì)形成佛光的云滴粒子大小進(jìn)行評(píng)估，為今后深入研究佛光形成的云物理過(guò)程及研判佛光出現(xiàn)的可能性提供分析依據(jù)，進(jìn)一步豐富利用天氣形勢(shì)定性預(yù)報(bào)佛光的技術(shù)方法[4]。

2 資料與研究方法

2.1 資料

德化縣九仙山氣象站位于閩中戴云山脈，海拔高度1653m，呈南北狹長(zhǎng)走向，上午、下午的太陽(yáng)運(yùn)行軌跡極易在狹長(zhǎng)山頂?shù)奈?、東向兩側(cè)山坡形成佛光現(xiàn)象。2020年1月17日下午，在測(cè)站下方出現(xiàn)稀薄流云，上方則為晴空狀態(tài)，14時(shí)50分，在測(cè)站東側(cè)下方盤(pán)山公路上出現(xiàn)了一次難得的近距離佛光現(xiàn)象(圖1)。經(jīng)現(xiàn)場(chǎng)測(cè)量，手機(jī)拍照點(diǎn)C與佛光中心點(diǎn)O的距離為46.4m，佛光半徑為3.2m，兩處的垂直高差CE為31.2m，15時(shí)的本站氣壓為83560Pa，氣溫8.5℃(281.5K)。

圖1 佛光實(shí)測(cè)數(shù)據(jù)示意圖

2.2 研究方法

佛光系由懸浮于大氣中的特定大小云滴粒子衍射太陽(yáng)光而成，由此即可通過(guò)流體動(dòng)力學(xué)原理求算大氣中的云滴粒子大小。此外，由于佛光是云滴粒子衍射太陽(yáng)光而出現(xiàn)的若干個(gè)圍繞一共同中心點(diǎn)的色圈[11]，色圈的大小與懸浮于大氣中的云滴粒子大小有關(guān)，每個(gè)色圈的半徑對(duì)應(yīng)相應(yīng)的衍射角，通過(guò)測(cè)量佛光圈徑大小等參數(shù)，也可得到相應(yīng)的衍射云滴粒子大小，二者可互為比對(duì)。

2.2.1 云滴粒子大小的動(dòng)力學(xué)計(jì)算

2.2.1.1 斯托克斯定律

在流體中運(yùn)動(dòng)的物體會(huì)受到流體的阻力作用，流體的粘滯阻力與速度的關(guān)系即為斯托克斯定律[12]，其表達(dá)式為：

f=-6πηνr

(1)

式中，f為在大氣中升降的云滴粒子所受到的空氣粘滯阻力，r是云滴粒子半徑，ν是云滴粒子相對(duì)于空氣的速度，η是空氣的粘滯系數(shù)。

2.2.1.2 云滴粒子的大小與下沉速度的關(guān)系

低層大氣中的水汽在不穩(wěn)定大氣的作用下上升而凝結(jié)為液態(tài)小水滴(即云滴粒子)，產(chǎn)生氣態(tài)到液態(tài)的質(zhì)變。在一定高度的云滴粒子，由于云滴粒子密度(即水密度)遠(yuǎn)大于空氣密度，使得小云滴粒子因所受的重力大于所受的浮力而下降，且速度增加，與此同時(shí)，空氣的粘滯阻力亦增加。隨著云滴粒子下降速度加大，向上的空氣粘滯阻力隨之增加，阻礙云滴粒子的下沉，當(dāng)重力、阻力和浮力達(dá)到平衡時(shí)(圖2)，則云滴粒子所受的外力合力為0而緩慢勻速升降，由于其速度很小，使得云滴粒子呈現(xiàn)懸浮狀態(tài)。

圖2 空氣中的云滴受力分析

r2=9ην/(2g(ρ-ρ0))

(2)

由于云滴粒子密度(103kg/m3)遠(yuǎn)大于空氣密度(1.29kg/m3)，則(2)式可簡(jiǎn)化為：

(3)

2.2.1.3 空氣粘滯系數(shù)η的計(jì)算

在溫度T<2000K時(shí)，空氣粘度η可用薩特蘭的冪次公式[13]進(jìn)行計(jì)算：

(4)

式中，η0為T(mén)0=15℃時(shí)的空氣粘滯系數(shù)，值為1.79×10-5Pa·s，對(duì)于空氣，指數(shù)n在90K

2.2.2 衍射原理的云粒子大小計(jì)算

2.2.2.1 衍射原理

陳世訓(xùn)等[11]指出，云滴粒子之間的間隙起著小圓孔的作用，光線通過(guò)云隙而向前發(fā)生衍射分光現(xiàn)象(圖略)，形成繞一共同中心點(diǎn)的多個(gè)衍射彩色圈，即佛光圈，每個(gè)衍射圈的色序由外向內(nèi)分別為紅、橙、黃、綠、青、藍(lán)、紫。設(shè)衍射角度為θ，則第k個(gè)暗紋條的位置可由下式確定：

sinθ=(k+0.22)λ/2r

(5)

式中，k=0，1，2，…。

2.2.2.2 光衍射的云滴粒子大小范圍

文獻(xiàn)[11]指出，產(chǎn)生衍射的云滴粒子尺度范圍可由如下光學(xué)衍射區(qū)的界限來(lái)確定：

3<2π(n-1)r/λ<20

(6)

其中r為云滴粒子半徑，水的折射率n=1.335。在該尺度范圍內(nèi)的云滴粒子可產(chǎn)生衍射現(xiàn)象，而尺度太大則發(fā)生折射，太小則發(fā)生散射。

3 觀測(cè)結(jié)果與討論

3.1 云滴粒子大小的計(jì)算結(jié)果

3.1.1 流體動(dòng)力學(xué)方法計(jì)算結(jié)果

由實(shí)測(cè)氣象要素?cái)?shù)據(jù)代入公式(4)，求得空氣的粘滯系數(shù)η=1.75×10-5Pa·s；再由公式(3)得到云滴粒子半徑r與終端速度ν之間的關(guān)系式為：

(7)

由上式得到云滴粒子大小r與終端速度ν的對(duì)應(yīng)關(guān)系結(jié)果，詳見(jiàn)表1。

表1 云滴粒子半徑與終端速度及紅光第一衍射角的對(duì)應(yīng)關(guān)系計(jì)算結(jié)果

3.1.2 衍射原理計(jì)算結(jié)果

設(shè)圖1中繞中心的紅光第一衍射圈的衍射角為θ，則：

sinθ≈tanθ=OA/AC=3.2/46.4=0.07

求得衍射角θ=arcsin0.07=4.0°。

由式(6)計(jì)算形成紅光第一衍射圈(k=1，λ=0.76μm)的云滴粒子大小：

sinθ=(1+0.22)×0.76/(2r)=0.07

由此得到云滴粒子半徑r=6.6μm。

3.1.3 太陽(yáng)光的衍射云滴粒子大小范圍

利用光學(xué)衍射區(qū)的界限公式(6)計(jì)算引起紅光(λ=0.76μm)發(fā)生衍射的云滴粒子尺度范圍為1.1～7.2μm，佛光中的其他顏色光色圈的衍射云滴粒子尺度范圍則小于紅光的衍射云滴粒子尺度范圍，如紫光(λ=0.38μm)的衍射云滴粒子尺度范圍為0.54～3.6μm。

3.2 討論

根據(jù)實(shí)測(cè)的佛光圈徑大小和與觀測(cè)者距離參數(shù)、氣象要素并分別利用衍射原理和流體動(dòng)力學(xué)原理計(jì)算形成佛光的衍射云滴粒子大小，討論如下：

①形成佛光的云滴粒子尺度范圍：大氣中的云滴粒子半徑通常在1～100μm之間[13]，由計(jì)算所得到的紅光衍射云滴粒子半徑尺度范圍(1.1～7.2μm)可知，并非所有的云系(云滴粒子)均可產(chǎn)生人眼可見(jiàn)的佛光現(xiàn)象，只有尺度較小的云滴粒子才有機(jī)會(huì)，因此佛光現(xiàn)象常出現(xiàn)在粒徑較小的稀薄云霧或云層上端。

②由實(shí)測(cè)佛光圈徑大小及位置數(shù)據(jù)并依衍射原理計(jì)算可得，本次形成佛光的衍射云滴粒子半徑為6.6μm，其滿足紅光的衍射云滴粒子半徑尺度范圍條件。

③經(jīng)對(duì)實(shí)測(cè)氣象要素進(jìn)行流體動(dòng)力學(xué)原理計(jì)算，該天氣條件下粒徑為6.6μm的云滴粒子的終端速度為0.54cm·s-1，符合懸浮狀態(tài)特征。

4 結(jié)論

根據(jù)本次近距離實(shí)測(cè)的佛光圈徑大小、與觀測(cè)者距離以及現(xiàn)場(chǎng)氣象要素，分別采用衍射原理和流體動(dòng)力學(xué)原理計(jì)算形成佛光的衍射云滴粒子半徑，初步得到以下結(jié)論：

①兩種方法計(jì)算形成佛光的云滴粒子大小的結(jié)果合乎實(shí)際情況，可相互比對(duì)驗(yàn)證；兩種估算方法簡(jiǎn)單合理，便于在實(shí)際中使用。

②通過(guò)實(shí)測(cè)的氣象要素和應(yīng)用流體動(dòng)力學(xué)原理并結(jié)合懸浮特征要求、光衍射的云滴粒子半徑尺度范圍條件綜合估算佛光衍射云滴粒子大小，此法有助于對(duì)距離較遠(yuǎn)等不便進(jìn)行佛光圈徑與距離實(shí)測(cè)的佛光衍射云滴粒子以及歷史佛光衍射云滴粒子的分析總結(jié)。

③實(shí)測(cè)氣象要素的流體動(dòng)力學(xué)估算方法，所得云滴粒子大小可供研判佛光出現(xiàn)的可能性，可為游客提供良好的觀賞服務(wù)，具有較好的氣象服務(wù)應(yīng)用價(jià)值。

④根據(jù)流體動(dòng)力學(xué)方法得到的佛光衍射云滴粒子大小，可進(jìn)一步依衍射原理推算衍射角，有助于分析形成佛光的云物理過(guò)程原理。

由于形成佛光的云物理過(guò)程機(jī)制尚無(wú)公認(rèn)的完整理論，通過(guò)借助佛光實(shí)測(cè)資料的分析，可幫助人們從量級(jí)上加深對(duì)形成佛光衍射云滴粒子的認(rèn)識(shí)，并對(duì)進(jìn)一步探索形成佛光的云物理過(guò)程原理和研判佛光出現(xiàn)可能性提供參考依據(jù)。目前一款探測(cè)范圍為0.1～1000μm的噴霧激光粒度儀正被用于九仙山氣象站，其今后所實(shí)測(cè)的佛光衍射云滴粒子大小可進(jìn)一步驗(yàn)證上述研究成果，并推動(dòng)佛光的深度研究。