魏海亮，邵利民，王寧寧

(1．海軍大連艦艇學(xué)院研究生管理大隊，遼寧 大連116018;2．海軍大連艦艇學(xué)院軍事海洋系，遼寧 大連116018;3．91040部隊，山東 青島266321)

1 引言

海霧是指海洋與大氣相互作用特定條件下出現(xiàn)在海上或沿海地區(qū)上空低層大氣的一種凝結(jié)現(xiàn)象，是懸浮于大氣邊界層中大量水滴或冰晶使得水平能見度小于1 km的危險性天氣現(xiàn)象［1］。激光在海霧中傳輸時要受到折射、透射、散射和吸收等作用，其中散射和吸收是引起激光衰減的主要原因。由于激光的波長與霧滴的尺度相當(dāng)，因此海霧對激光的衰減值得關(guān)注。目前國內(nèi)外已經(jīng)建立了許多用于計算激光大氣透過率的大型數(shù)據(jù)庫和計算包，比較典型的如LOWTRAN7［2］、MODTRAN［3］等，但其分辨率較低，計算適用范圍有限，計算較為繁瑣;基于能見度的經(jīng)驗霧衰減計算模型在工程中應(yīng)用廣泛［4］，但模型僅考慮能見度和激光波長，沒有更多考慮霧的其他物理特征;基于Mie散射理論的衰減計算模型考慮了霧的滴譜分布等因素［5－6］，由于是以單次散射為基礎(chǔ)，當(dāng)霧能見度較小時該模型就不再適用。

海霧的電磁性質(zhì)是基于傳輸激光波長、環(huán)境溫度等參數(shù)的非線性函數(shù)，很難從理論上直接給出它們與介電常數(shù)的確切函數(shù)表達(dá)式。研究隨機(jī)介質(zhì)介電常數(shù)的一般做法是在特定溫度下采用測量方法得到經(jīng)驗公式，如已經(jīng)取得的雨滴、冰、雪和沙塵等介電常數(shù)的經(jīng)驗公式，這些方法反映了隨機(jī)介質(zhì)中粒子的介電特性，對研究激光的傳輸效應(yīng)有重要作用，但不能充分反映隨機(jī)介質(zhì)整體的介電特性，不能體現(xiàn)介質(zhì)物理特性對等效介電特性的影響。本文根據(jù)電介質(zhì)極化理論推導(dǎo)出了霧滴粒子的等效介電常數(shù)的一般計算模型，并進(jìn)一步得出傳輸衰減的計算公式。

2 改進(jìn)的等效介電常數(shù)模型

海霧發(fā)生時，實際是霧滴粒子按照一定的譜分布彌散在空氣中，形成一種宏觀近似均勻的隨機(jī)介質(zhì)。假設(shè)霧滴的等效介電常數(shù)為εeff，則霧區(qū)中的電位移矢量為［7］:

其中，εw代表霧滴粒子的介電常數(shù)，即水的介電常數(shù)。

利用上式計算霧滴粒子的等效介電常數(shù)最大的局限就是靜場假設(shè)，因此為了計算激光條件下的等效介電常數(shù)，可用下式粒子極化率代替［9］:

其中，f為激光頻率;fr為霧滴粒子的共振頻率。

以平流霧經(jīng)驗霧滴譜模型為例，仿真、計算了霧區(qū)的等效介電常數(shù)，水的介電常數(shù)用雙Debye公式計算，霧滴半徑取0．01～36μm，溫度為10℃。圖1為霧區(qū)等效介電常數(shù)仿真和計算結(jié)果。

圖1 霧區(qū)等效介電常數(shù)虛部

從圖1中可以看出海霧介質(zhì)的等效介電常數(shù)是頻率和能見度的非線性函數(shù);同時當(dāng)能見度變化時霧滴譜分布也隨之不同，因此海霧的霧滴尺寸及滴譜分布對霧介質(zhì)等效介電常數(shù)也有一定的影響;而溫度對等效介電常數(shù)的影響則體現(xiàn)在其對水的介電常數(shù)的作用上，不過溫度對霧滴粒子的電磁傳輸特性影響相對較?。?］。

利用等效介電常數(shù)計算激光霧中傳輸特征衰減(衰減率)的公式為［10］:

式中，c0為真空光速。

3 仿真與討論

3．1 海霧的霧滴譜特征

霧滴譜是霧中不同大小霧滴數(shù)密度的分布，是表征霧結(jié)構(gòu)的一個重要特征參量。被廣泛使用的霧滴譜模型是Khragian－Mazin分布模型，利用含水量與能見度的關(guān)系可以得到經(jīng)驗平流霧滴譜分布如下:

其中，V為能見度;r為霧滴半徑。從式(7)中可以看出，改變能見度的大小就可以獲得不同的霧滴譜分布。

由于不同氣象條件下海霧過程都有其各自特點，因此霧滴譜的研究也多是個例研究，下面是分別根據(jù)舟山地區(qū)和滬寧地區(qū)兩次海霧過程實測數(shù)據(jù)擬合得出的霧滴譜分布［11－12］:

圖2為海霧的滴譜分布，其中圖2(a)為不同能見度情況下平流霧經(jīng)驗霧滴譜分布，隨著海霧能見度的降低，霧滴譜寬度變寬，霧滴數(shù)濃度峰值也隨之向大霧滴方向變化;圖2(b)為當(dāng)V=0．5 km時，平流霧經(jīng)驗霧滴譜分布與兩個海霧個例霧滴譜分布的比較，由圖中的霧滴譜分布可見，經(jīng)驗霧滴譜和滬寧地區(qū)海霧個例的霧滴譜較窄，主要表現(xiàn)為r＜15μm的小霧滴;而舟山地區(qū)海霧個例的霧滴譜則相對較寬，大的霧滴明顯增多，霧滴數(shù)濃度峰值在r=20μm附近。

圖2 霧的滴譜分布

3．2 激光霧中傳輸衰減仿真

在工程應(yīng)用中，激光霧中傳輸衰減計算通常采用經(jīng)驗衰減模型如Kim模型、Naboulsi模型等［4］。圖3(a)、(b)分別以0．63μm和10．6μm波長為例，對等效介電常數(shù)衰減模型與經(jīng)驗衰減模型的仿真比較，其中霧滴譜采用平流霧經(jīng)驗霧滴譜模型。從圖中可以看出，各模型中激光霧區(qū)衰減率隨能見度變化趨勢相同。當(dāng)激光波長為0．63μm時，兩經(jīng)驗?zāi)Ｐ陀嬎憬Y(jié)果基本重合;在V＞2 km情況下激光波長不管是0．63μm還是10．6μm，等效介電常數(shù)衰減模型與經(jīng)驗?zāi)Ｐ皖A(yù)測結(jié)果都差別不大，隨著能見度的降低特別是能見度低于1 km時各模型計算的衰減率差別變大。這是因為經(jīng)驗?zāi)Ｐ驮谟嬎銜r僅考慮了海霧的能見度因素，而能見度是海霧的宏觀物理特征，當(dāng)能見度很低時霧滴粒子多散射效應(yīng)對衰減率的影響不容忽視。式(4)是根據(jù)電介質(zhì)極化理論導(dǎo)出了海霧發(fā)生時霧滴粒子的等效介電常數(shù)的一般計算模型，在推導(dǎo)過程中考慮了霧滴粒子的平均退極化場和單個粒子的退極化場分別對粒子外部場和粒子內(nèi)部場的影響，考慮了霧滴粒子物理特性對等效介電常數(shù)的影響。相比經(jīng)驗?zāi)Ｐ停锢頇C(jī)理和概念清晰，能夠反映海霧物理特征對等效介電特性的影響，能夠方便地計算具有不同物理特征霧滴粒子的等效介電特性。

圖3 激光霧區(qū)衰減率隨能見度的變化曲線

從圖4(a)、(b)顯示了不同霧滴譜分布對激光衰減率的影響，其中圖4(a)采用平流霧經(jīng)驗霧滴譜，不同能見度情況下霧滴譜分布變化很大，如圖2(a)所示;圖4(b)中平流霧經(jīng)驗霧滴譜中V=0．5 km，各霧滴譜分布的差異從圖2(b)中可以看出。從圖中可以看出霧滴譜分布對激光衰減率的大小影響非常明顯，各霧滴譜分布的衰減率差別很大，體現(xiàn)了激光衰減對海霧環(huán)境物理特征的依賴。

圖4 不同霧滴譜分布對激光衰減率的影響

4 結(jié)束語

海霧是一種具有吸收或者散射特性的傳輸介質(zhì)，研究激光在海霧中傳輸?shù)乃p特性，對于了解海霧的光電遮蔽性能具有重要意義。根據(jù)電介質(zhì)極化理論推導(dǎo)出了激光霧區(qū)傳輸?shù)牡刃Ы殡姵?shù)衰減計算模型，在推導(dǎo)過程中考慮了霧滴粒子的平均退極化場和單個粒子的退極化場分別對粒子外部場和粒子內(nèi)部場的影響，并充分考慮了海霧環(huán)境本身的物理特征對激光衰減特性的影響。本文建立的模型物理機(jī)理和概念清晰，提供了一種計算激光霧中衰減的新方法。由于相關(guān)的海霧衰減數(shù)據(jù)稀缺，本文的研究工作局限于理論分析，今后可結(jié)合實測數(shù)據(jù)對模型進(jìn)一步驗證。

［1］ XU Xiaofeng，GU Jianfeng，LI Yongping．Marine meteorological disasters［M］．Beijing:Meteorological Press，2009:102－104．(in Chinese)許小峰，顧建峰，李永平．海洋氣象災(zāi)害［M］．北京:氣象出版社，2009:102－104．

［2］ F Filbet，G Russo．High order numerical methods for the space non－h(huán)omogeneous boltzmann equation［J］．Comput Physics，2003，186:457－480．

［3］ LIU Weicaho，QI Linlin，HE Hongrang，et al．Numerical study on atmospheric transmission of 1．06μm laser［J］．Laser＆Infrared，2011，41(5):520－524．(in Chinese)劉偉超，齊琳琳，何宏讓，等．1．06μm激光大氣透過率特性的數(shù)值計算研究［J］．激光與紅外，2011，41(5):520－524．

［4］ LI Di，CHEN Hui．Study on attenuation characteristic of laser propagation in rain and fog［J］．Electronic Design Engineering，2011，19(9):1－5．(in Chinese)李娣，陳輝．激光大氣傳輸?shù)挠觎F衰減特性研究［J］．電子設(shè)計工程，2011，19(9):1－5．

［5］ LI Xuebin，GONG Chunwen，LI Chao，et al．Size distribution and attenuation for infrared radiation of fog particles［J］．Laser＆Infrared，2009，39(7):743－745．(in Chinese)李學(xué)彬，宮純文，李超等．霧滴譜分布和霧對紅外的衰減［J］．激光與紅外，2009，39(7):743－745．

［6］ KE Xizheng，MA Dongdong，LIU Jiani．Study attenuation of laser transmission in fog［J］．The Journal of Light Scattering，2009，21(2):104－109．(in Chinese)柯熙政，馬冬冬，劉佳妮．激光在霧中傳輸?shù)乃p研究［J］．光散射學(xué)報，2009，21(2):104－109．

［7］ PENG Weicai，CHEN Kanghua．Calculation of effective electromagnetic parameters in random mixture media of magnetic iron fiber［J］．Rare Metal Materials and Engineering，2005，34(9):1047－1410．(in Chinese)彭偉才，陳康華．磁性纖維隨機(jī)混合媒質(zhì)等效電磁參數(shù)的計算［J］．稀有金屬材料與工程，2005，34(9):1047－1410．

［8］ P Lorrain，D Corson．Electromagnetic fields and waves［M］．Beijing:People's Education Press，1980:62－87．(in Chinese)P勞蘭，D R考森．電磁場與電磁波［M］．北京:人民教育出版社，1980:62－87．

［9］ ZHOU Xingli，XIE Yanan，YANG Zhengde．Effect study of clouds and fogdrops on microwave transmission characteristics［J］．Meteorological Science and Technology，2011，39(5):661－665．(in Chinese)周星里，謝亞楠，楊正得．云霧對電磁傳輸特性影響的研究［J］．氣象科技，2011，39(5):661－665．

［10］M Z Gao，T SYeo，P SKooi，et al．Rain attenuation calculation using dielectric mixture with deformed rain drops［C］．Geoscience and Remote Sensing Symposium，IGARSS'93．Better Understanding of Earth Environment．，International，1993:1250－1252．［11］YANG Zhongqiu，XU Shaozu，GENG Biao．Formation and microphysical structure features of spring sea fog in Zhoushan region［J］．Acta Oceanologica Sinica，1989，11(4):431－438．(in Chinese)楊中秋，許紹祖，耿驃．舟山地區(qū)春季海霧的形成和微物理結(jié)構(gòu)［J］．海洋學(xué)報，1989，11(4):431－438．

［12］HUANG Jianping，MEI Qingyin，JIN Yongcai．Microphysical structure features and evolution processes of radiation fog in Huning region［J］．Meteorological Monthly，1998，24(5):3－7．(in Chinese)黃建平，梅清銀，靳永才．滬寧地區(qū)輻射霧的微物理結(jié)構(gòu)及其演變［J］．氣象，1998，24(5):3－7．