懸空波導(dǎo)的激光雷達(dá)探測能力個例研究

馬征征 王曉賓

（中國電波傳播研究所 電波環(huán)境特性及模化技術(shù)重點實驗室，山東 青島266107）

引 言

激光雷達(dá)作為一種遙感探測設(shè)備，能夠始長期通過測量對流層大氣溫度和濕度來間接獲取大氣折射率.早在20世紀(jì)70年代，美國賓州州立大學(xué)已開致力于研究可用于測量電波環(huán)境的拉曼激光雷達(dá)［1］，其第三代和第五代系統(tǒng)于1990至2002年間參與了十幾項行動計劃.同時科研人員還將其第五代激光雷達(dá)系統(tǒng)和無線電探空儀的探測結(jié)果進(jìn)行了比對［2］.比對結(jié)果表明，激光雷達(dá)系統(tǒng)具有精確、連續(xù)、即時、實時、無需額外支持（如氣瓶、氣球）、不怕風(fēng)以及無準(zhǔn)備時間等諸多優(yōu)勢，但其問題是垂直分辨率較低且白天噪聲大.另外，英國奇爾波頓天文臺于2002年研制了一套用于探測大氣波導(dǎo)的紫外拉曼激光雷達(dá)系統(tǒng)［3］.21世紀(jì)初，瑞士聯(lián)邦理工學(xué)院也研制了探測大氣溫度和濕度的移動型拉曼激光雷達(dá)系統(tǒng)［4］.

相比傳統(tǒng)手段，激光雷達(dá)雖然具有連續(xù)、即時和遙感等諸多優(yōu)點，但其目前在某些方面仍有待加強(qiáng)和改進(jìn).例如，受光路結(jié)構(gòu)限制，激光雷達(dá)通常具有一定的探測起始距離，這一定程度上限制了貼地波導(dǎo)的探測.同時，激光雷達(dá)的探測原理表明，其接收信號強(qiáng)度大致與探測距離的平方成反比關(guān)系.因此，相比貼地波導(dǎo)，處于更高高度的懸空波導(dǎo)對探測系統(tǒng)性能提出了較高要求.另外，鑒于大氣波導(dǎo)較窄的層高（約十幾米至幾十米），對激光雷達(dá)的空間分辨率也提出了一定要求.

本文將針對一例通過無線電探空儀實測的懸空波導(dǎo)，來仿真考察激光雷達(dá)接收信號泊松起伏所引入誤差對其探測能力的影響.對于其它能夠影響探測精度的因素，如雷達(dá)定標(biāo)、擬合式和經(jīng)驗公式引入、彈性信號和背景噪聲混入等，這里不作研究［5］.

1 探測原理

采用使用最廣泛的拉曼激光雷達(dá)［6］.其通過振-轉(zhuǎn)拉曼散射信號測量大氣濕度，通過純轉(zhuǎn)動拉曼散射信號測量大氣溫度.拉曼激光雷達(dá)方程為

式中：系數(shù)K＝AN0Δrη，是與距離無關(guān)的量；體后向散射系數(shù)β（λ0，λX，r）＝n（r）σ（λ0，λX）；消光系數(shù)α＝αa＋αm；X表示不同的物質(zhì)種類；其他各量的含義見表1.

表1 公式（1）所使用的物理量

其中，噪聲NB可以通過擬合消除.

對于濕度探測，通常采用對H2O和N2分子兩個通道測量結(jié)果的比值來確定，有

在完成各個模塊的FPGA硬件實現(xiàn)后,先要通過Modelsim進(jìn)行功能仿真,然后進(jìn)行板級測試,驗證模塊功能是否正確。整套系統(tǒng)軟硬件測試平臺如圖10所示。

式中：e為水汽分壓；PN2為氮氣分壓，均隨高度分布.

對于溫度探測，通常采用兩個不同波長（或波段內(nèi)）純轉(zhuǎn)動拉曼散射信號的測量結(jié)果來反演.由于兩個通道波長相對接近，消光系數(shù)近似相同.兩通道比值為

純轉(zhuǎn)動拉曼散射截面與溫度有關(guān)，比值中包含了溫度信息，即Q＝Q（T）.再通過與理論計算的或標(biāo)定的Q-T關(guān)系比對，即可得到溫度.更具體的反演過程以及散射截面和消光系數(shù)的計算可參考早期的工作［5］.作為簡單情況，設(shè)定溫度探測中兩個通道均探測純轉(zhuǎn)動拉曼散射的單支譜線.具體設(shè)置為氮氣的第4和第14支反-斯托克斯轉(zhuǎn)動能級譜，這兩根譜線的散射截面與溫度成相反變化關(guān)系.

相比濕度和溫度，氣壓是相對穩(wěn)定的量，這里直接使用無線電探空儀的測量結(jié)果.為了將地球表面處理成平面，通常使用經(jīng)地球曲率修正的大氣修正折射指數(shù)M（單位M-unit）.對于30～100GHz范圍內(nèi)的電磁波可計算得到大氣修正折射指數(shù)為

式中：P、T、e分別是大氣壓、溫度和水汽分壓，單位依次是hPa、K和hPa；h是高度；Re是地球平均半徑.

對于實際情況，無論是濕度測量還是溫度測量，各通道獲得的信號（光子數(shù)）中包含了起伏噪聲.根據(jù)激光雷達(dá)的探測原理，該噪聲服從泊松分布，有

其將會在反演結(jié)果中引入誤差.式中，N＊（r）是無信號起伏時的接收光子數(shù).

2 仿真結(jié)果

仿真使用2009年6月17日約9點在廣東省汕頭市南澳島通過高空氣球探測得到的一例懸空波導(dǎo)，圖1給出了該懸空波導(dǎo)的大氣修正折射指數(shù)高度剖面.

圖1 2009年6月17日廣東省汕頭市南澳島高空氣球探測的大氣修正折射指數(shù)高度剖面

探空儀時間分辨率為2s，近似勻速上升.2.5 km高度以下對應(yīng)的高度分辨率為10.4m.其中1 km高度以下的平均高度分辨率為10.7m；而1km高度以上的平均高度分辨率為10.3m.

該次實驗中測量到了兩個大氣波導(dǎo)結(jié)構(gòu)，均屬于較弱的例子.其參數(shù)如表2所示.這里以較高高度的大氣波導(dǎo)為例.

表2 無線電探空儀測量的大氣波導(dǎo)參數(shù)

對于激光雷達(dá)，考慮一套中型系統(tǒng).采用美國Continuum公司的Surelite系列激光器.具體參數(shù)為：發(fā)射波長為355nm，重復(fù)頻率為20Hz，脈沖能量為80mJ；發(fā)射累積時間為10min.根據(jù)公式

計算得到發(fā)射總光子數(shù)（～1021）.對于中型系統(tǒng)，接收望遠(yuǎn)鏡直徑設(shè)為40cm.對于大氣波導(dǎo)幾十米的層高，高度門為200ns即15m.雷達(dá)效率簡單設(shè)定為0.01［7］.

按照公式（5）在激光雷達(dá)測量濕度和溫度的四個通道接收信號中參入泊松分布的信號起伏，并對反演結(jié)果使用了寬度為30m的平滑.圖2給出了一次加入信號起伏后的探測結(jié)果，可以看到激光雷達(dá)成功地捕獲了該懸空波導(dǎo).圖3給出了另一次加入信號起伏后的探測結(jié)果.很明顯，由于信號起伏的影響，本次探測未能辨識出該懸空波導(dǎo).

圖2 一次激光雷達(dá)仿真探測結(jié)果（成功捕獲到了該懸空波導(dǎo)）

基于蒙特卡洛（Monte-Carlo，MC）方法引入泊松分布的信號起伏，考慮10 000次重復(fù)試驗，得到激光雷達(dá)對該懸空波導(dǎo)的辨識概率，以及波導(dǎo)厚度、波導(dǎo)強(qiáng)度測量結(jié)果，如表3所示.

圖3 另一次激光雷達(dá)仿真探測結(jié)果（未能辨識出該懸空波導(dǎo)）

表3 激光雷達(dá)仿真測量的大氣波導(dǎo)參數(shù)

表3結(jié)果表明，在時空分辨率為10min和15 m時，激光雷達(dá)對該懸空波導(dǎo)的辨識概率較高，達(dá)到了92%.層高的測量值也與真實值十分接近.同時，表3還進(jìn)一步列出了其它時空分辨率下的仿真結(jié)果.意料之中的，當(dāng)探測累積時間減少一半時，激光雷達(dá)的辨識概率降低至88%.很明顯這是因為隨著累積時間的減少，接收信號的光子數(shù)同時減少，從而造成信號起伏相對增強(qiáng).另外，空間分辨能力提升1倍（至7.5m）對辨識概率基本沒有影響.這是因為空間分辨能力提升的同時單一高度門內(nèi)的信號光子數(shù)也會減少，從而抵消了空間分辨率能力提升帶來的益處.此外，表3結(jié)果表明，各時空分辨率情況下，波導(dǎo)強(qiáng)度的測量結(jié)果與真實值均有一定程度的偏差.波導(dǎo)強(qiáng)度測量精度應(yīng)受限于激光雷達(dá)系統(tǒng)的整體性能.當(dāng)然，無論是激光雷達(dá)功率孔徑積的增大，還是系統(tǒng)效率的提升，均能提高激光雷達(dá)系統(tǒng)對大氣波導(dǎo)的探測能力.

3 結(jié) 論

我們仿真考察了一臺中型拉曼激光雷達(dá)系統(tǒng)在三種時空分辨率下對一例實測較弱懸空波導(dǎo)的探測能力，包括探測的辨識概率、波導(dǎo)厚度和強(qiáng)度.

仿真結(jié)果表明，對于該例較弱的懸空波導(dǎo)，中型激光雷達(dá)系統(tǒng)對其具有較高的捕獲能力，辨識概率能達(dá)到90%左右，且辨識概率隨著探測累積時間的增加而提升.同時在三種不同時空分辨率下，雷達(dá)系統(tǒng)對波導(dǎo)層厚度的測量結(jié)果均十分接近真實值.相比而言，雷達(dá)對波導(dǎo)強(qiáng)度的測量值誤差要稍大一些，但也是可以接受的.

總的來說，利用激光雷達(dá)探測懸空波導(dǎo)能夠具有可觀的預(yù)期.進(jìn)一步地，我們期望以此工作的個例研究為契機(jī)，開展激光雷達(dá)對大氣波導(dǎo)探測能力的系統(tǒng)研究.當(dāng)然，這還需要進(jìn)一步的理論和實驗工作.

［1］PHILBRICK C R.EM ＆EO properties of the lower atmosphere［C］／／BACIMO 2003Conference.Monterey，2003.

［2］HARRISOND E.Comparison of Lidar and Mini-rawin Sonde Profiles［D］.Monterey：Naval Postgraduate School，1998.

［3］AGNEW J L.Lidar Measurement of Tropospheric Radio Refractivity［R］.Radiocommunications Agency Project，2002.

［4］FROIDEVAUX M，SERIKOV I，BURGOS S，et al.A new lidar for water vapor and temperature measurements in the atmospheric boundary layer［J］.Asia Flux Newsletter Issue，2009，28：13-17.

［5］馬征征，王曉賓，康士峰，等.純轉(zhuǎn)動拉曼測溫激光雷達(dá)濾光器參數(shù)設(shè)計［J］.電波科學(xué)學(xué)報，2012，27（5）：937-942.MA Zhengzheng，WANG Xiaobin，KANG Shifeng，et al.Designing filter parameters of pure rotational raman lidar for detecting duct［J］.Chinese Journal of Radio Science 2012，27（5）：937-942.（in Chinese）

［6］WEITKAMP C.Lidar：Range-resolved Optical Remote Sensing of the Atmosphere［M］.Springer，2005：102.

［7］POVEY A C.Development of an optimal estimation retrieval scheme for a Raman lidar system［R／OL］.［2013-09-12］（2011-12-08）http：／／www.atm.ox.ac.uk／group／eodg／1st-year-reports／adam-povey.pdf.