測風激光雷達信號載噪比和有效性算法研究

羅 甜，胡芳仁，朱海龍

(1.南京郵電大學 光電工程學院，江蘇 南京 210046;2.南京先進激光技術研究院先進全固態激光技術研發中心，江蘇 南京 210038)

?

測風激光雷達信號載噪比和有效性算法研究

羅 甜1，胡芳仁1，朱海龍2

(1.南京郵電大學 光電工程學院，江蘇 南京 210046;2.南京先進激光技術研究院先進全固態激光技術研發中心，江蘇 南京 210038)

信號的載噪比和有效性是衡量測風激光雷達系統性能的重要指標。在不同的去噪方法基礎上提出了信號載噪比的新算法，并通過與國外測風激光雷達進行外場對比實驗驗證了該載噪比算法的正確性以及該算法對不同去噪方法的適用性。同時,提出了判別信號有效性的新算法，并通過實驗數據對比分析得出，載噪比和有效性新算法在判別信號有效性上具有很高的一致性，從而證明了該有效性算法的正確性。此外，與載噪比相比，有效性算法能進一步縮小信號的有效范圍，提高有效精度，從而能大大提高反演出的風場信息的準確性。

測風激光雷達；載噪比；有效性；信號分析

0 引言

多普勒測風激光雷達具有測量精度高、時空分辨率高及響應速度快等優點[1],廣泛應用于大氣湍流探測、風切變測量[2]、大氣風廓線測量[3]和機場測風等，其在航空航天和軍事應用等領域都有著十分重要的作用。發展相干多普勒測風激光雷達進行風場探測已成為國內外激光雷達的一個重要發展方向。近年來，隨著測風激光雷達技術的不斷發展和完善[4-5]，國內外出現了多種類型的測風激光雷達，這些新型的測風雷達系統都是針對硬件進行優化或對探測技術進行改進[6-8]等，少有針對雷達系統回波信號的載噪比算法進行簡化和信號有效性進行精確判別的研究。縱觀國內外，空難事故頻繁發生，其中很大一部分原因是測風激光雷達回波信號反演出的風場信息可靠性差。如何利用新的科學手段及時判別風場信息的可靠性，將是我們能及時做出相應應急措施，從而有效地降低空難事故發生概率的關鍵。

本文針對如何判別測風激光雷達回波信號的可靠性,提出了判別測風激光雷達回波信號有效性的新算法。由于回波信號載噪比也能在一定程度上判別回波信號的有效性，同時它也是衡量測風激光雷達系統性能的重要指標，因此本文還提出了回波信號載噪比的新算法。

文中所用數據是由南京先進激光技術研究院所研制的相干多普勒測風激光雷達于2015年5月13日在青島即墨所測得的。該測風激光雷達發射信號脈寬為180 ns,重頻為10 kHz，中頻為120 MHz。

1 回波信號去噪分析

由于測風激光雷達回波信號包含大量的噪聲，因此在進行載噪比計算和有效性判別之前需要對測風激光雷達回波信號進行去噪。目前常用的去噪方法有：減最后距離門回波信號去噪法[9]和減當前距離門擬合基底噪聲去噪法[10]。減最后距離門回波信號去噪法認為相干多普勒測風激光雷達最后一個距離門即激光雷達最高探測距離內的回波信號中沒有氣象信號，全部是噪聲，所以可以將最后一個距離門的回波信號作為基底噪聲進行去噪。晴空時，認為測風激光雷達最后距離門內的回波信號全是噪聲信號是基本合理的，但在其他天氣(如降水天氣)時，最后距離門的回波信號包含有不可忽略的氣象信息。這時若還以最后距離門回波信號作為噪聲功率標準，就會給譜參數估計帶來較大的誤差。減當前距離門擬合基底噪聲去噪法是通過多項式最小二乘法對每個距離門內的功率譜基線進行擬合，最后將擬合出來的曲線作為當前距離門的噪聲進行去噪。該算法的優勢是其噪聲功率譜是針對每個距離門單獨計算得到的，可以在一定程度上減少由于基底噪聲自身的變化抖動帶來的不利影響。

圖1為640 m和1 540 m高空處，沒有經過任何去噪處理的回波信號與分別經過上述兩種不同的去噪方法后所得的回波信號的功率譜對比圖。

圖1 640 m和1 540 m處不去噪和兩種不同方法去噪后的的回波信號功率譜對比圖

圖中實線為未經去噪處理的回波信號的功率譜線，虛線為減去最后距離門去噪的回波信號的功率譜線，點線為減基底擬合噪聲后的功率譜線。同時,線上的圓點分別表示在指定的頻率區間內功率譜峰值點。由圖1(a)和圖1(b)可以明顯看出，未經去噪處理的回波信號功率譜在160 MHz之后開始出現大幅上升趨勢。這種上升趨勢是由測風激光雷達探測器自身噪聲所引起的。經去噪后的頻譜線與未去噪處理的頻譜線相比更為平滑，且在160 MHz之后沒有出現幅值大幅上升的現象，這就說明上述兩種去噪方法都能在一定程度上抑制激光雷達系統的探測器所帶來的不利影響。同時，從圖1(b)可以看出，高空探測時，去噪和不去噪的回波信號功率譜峰值并不在同一頻率點。經去噪后的信號功率峰值頻率更接近中頻120 MHz,說明去噪后所得到的峰值頻率更接近真實的峰值頻率。所以在反演風場信息之前應先對回波信號進行去噪處理。

圖2為上述兩種不同去噪方法對同一回波信號處理后的功率譜對比圖。

圖2 640 m和1 540 m處減最后距離門去噪和減當前距離門擬合噪聲去噪的功率譜對比圖

圖2中實線為減最后距離門去噪的功率譜線，虛線為減當前距離門擬合基底噪聲去噪的功率譜線。從圖2(a)可以看出，在低空探測時，這兩種去噪方法具有很好的一致性，即功率譜峰值點相同。但圖2(b)顯示出高空探測時，不同的去噪處理所得到的功率譜的峰值點不在同一點，即得到的多普勒頻移量不同，最終會導致反演出的風場信息不同。

圖3是最后距離門噪聲和640 m當前距離門擬合基底噪聲的噪聲對比圖。

圖3 最后距離門噪聲和640 m當前距離門擬合噪聲對比圖

從圖3中可以看出最后距離門噪聲和640 m當前距離門擬合基底噪聲并不完全相同，640 m當前距離門的擬合基底噪聲與最后距離門噪聲相比更為平滑，這有利于將氣象信號的功率譜峰值從噪聲中區分出來。

綜上所述，兩種不同的去噪方法在低空探測時具有相同的去噪效果，但在高空探測時存在一定的差異，即會得到不同的風場信息。而通過上述理論分析可知減當前距離門擬合基底去噪更為合理和精確，但該去噪方法計算復雜，需要對每一個距離門進行基底擬合，運算時間長，而減最后距離門去噪不需要這個過程，這將大大提高運算速度，節省時間。因此可以根據測風激光雷達系統的不同探測距離和運算時間來合理選擇不同的去噪方法。

2 載噪比的算法設計及其驗證

激光雷達回波信號和噪聲之間的比值與激光雷達的探測范圍和精度密切相關，是評估激光雷達系統性能的一個重要指標。載噪比[8,10-11]的理論公式為：

(1)

式中，η為系統效率，Ex為激光器發射能量，β為氣溶膠后向散射系數，λ為波長，T為傳輸R處的大氣透射率，D為望遠鏡孔徑，h為普朗克常量，B為探測電路帶寬，R為探測距離。可以看出上述載噪比理論公式包含很多的參量。

本文提出的載噪比算法不需要上述載噪比理論公式中的眾多參量，其基本思路是在經去噪后的回波信號功率譜的基礎上，根據雷達發射信號的脈寬T，取以峰值頻率為中心的一段連續的2/T區間，認為該2/T區間內功率譜曲線的面積為氣象信號的面積S1，剩下的面積為噪聲的面積S2，由此得出初步去噪后的信號載噪比CNR=S1/S2。

如圖4所示，經過去噪后的回波信號功率譜圖中，兩條虛線是根據峰值頻率(圓點所在頻率)左右各頻移1/T=5.5 MHz對應的頻率。依據上述載噪比的計算思路可認為兩條虛線之間的面積是氣象信號的面積，而剩余部分為噪聲的面積，由此可計算出載噪比。

圖4 氣象信號面積示意圖

為了驗證該算法的正確性，通過與國外測風激光雷達做對比實驗的方法來進行驗證分析。將兩臺測風激光雷達放置于同一場地，統一指向，使其外界條件達到一致，同時測量大氣風場信息，最后將本文提出的新算法算得的載噪比與國外測風激光雷達自身處理得到的載噪比作對比分析來驗證該算法是否正確。

圖5為兩臺測風激光雷達于2015年5月13日在即墨所測的在40 m處連續17個小時(1 020分鐘)的數據對比圖。雷達1為南京先進激光技術研究院所研制的測風激光雷達，雷達2為國外測風激光雷達。由于兩個激光雷達系統的載噪比算法不同，因此在分析載噪比對比圖時，只需看兩者的變化趨勢是否一致，即對比是否是在同一點上升、下降即可。因為在相同的時刻出現相同的變化趨勢就說明這兩臺測風激光雷達對變化的風場探測到了相同的結果。

圖5 載噪比對比圖

如圖5所示，兩條曲線的變化趨勢是一致的，在相同的時刻出現拐點，即說明兩臺測風激光雷達對變化的風場探測到的結果相同，這也就驗證了本文提出的載噪比算法的正確性。

由于該載噪比算法是在回波信號經過一定的去噪處理基礎上算得的，因此，在驗證了本文提出的載噪比算法的正確性的基礎上，還需驗證該載噪比算法對上述兩種不同的去噪方法是否均適用。

圖6是在上述兩種不同的去噪方法基礎上所得到的載噪比對比圖，如圖所示，測風激光雷達系統設有200個距離門，每個距離門為30 m。圖中兩條曲線的變化趨勢基本一致，即說明該載噪比的算法對上述兩種去噪方法均適用。同時，從圖中可以看出，從第14距離門到第50距離門左右，其載噪比呈明顯的下降趨勢，這符合式(1)中載噪比與探測距離成反比的趨勢。而第50距離門之后的載噪比在總體趨勢上是相對平穩的，這與式(1)是矛盾的，即說明該回波信號是不可靠的，超出了測風激光雷達系統的探測范圍。因此，可利用載噪比的這種變化趨勢來初步判別信號的有效性。

圖6 載噪比隨探測距離變化的對比圖

3 信號的有效性判別

本文除了用上述載噪比的變化趨勢來判別信號的有效性外，還提出了另一種方法來判別信號的有效性，即：通過最大概似法對其去噪后的功率譜幅值進行均值估計，然后將該功率譜的峰值m與最大概似法估得的均值n作比值c=m/n，根據比值c的變化趨勢來判斷信號的有效性。

圖7為減最后距離門去噪時c的變化趨勢,從圖中可以看出c在第10到40距離門范圍內基本穩定，當探測距離超過第40距離門之后,c值有較大的波動。分析其原因在于第10距離門之前和第40距離門之后，通過最大概似法估算出的均值n為負值，而若此時的回波信號是有效信號，則它的功率譜幅值必定都大于零，因為只有當噪聲將氣象信號完全淹沒時才會出現均值n小于零的情況，而此時的信號是不可靠的，即無效的。由此可根據比值c的這種變化趨勢來判斷信號的有效性。

圖7 減最后距離門去噪時,比值c隨距離門的變化趨勢圖

與上文中直接用載噪比來判別信號的有效性相比，用峰值與均值的比值來判別信號的有效性能縮小有效探測距離的范圍，這有利于進一步提高最后反演得到的風速信息的可靠性。

圖8為減當前距離門擬合基底噪聲去噪時, 其功率譜最大值m與最大概似法估得的均值n的比值c的變化趨勢圖。

圖8 減擬合基底噪聲去噪時,比值c隨距離門的變化趨勢圖

由圖8可以看出，從第11距離門到第40距離門，其比值c一直是緩慢增加的，40距離門之后基本趨于平穩，若將c緩慢增加的區間看作是信號有效區間，則其有效區間與圖8中的有效區間是完全一致的，由此可見減擬合基底噪聲去噪時也可根據c的變化趨勢來判別信號的有效性。

綜上所述，減最后距離門去噪時，比值c為正且趨于平穩的連續區間為信號的有效區間；減擬合基底去噪時，比值c呈上升的連續區間為信號的有效區間。

4 結論

本文對于目前存在的兩種去噪方法進行了對比分析，結果顯示在低空探測時，這兩種去噪方法在多普勒頻移量上具有很好的一致性，但在高空探測時存在一定的差異。經理論分析得出，減當前距離門擬合基底去噪更為精確，但該方法計算復雜，耗時長。因此應針對測風激光雷達系統所需的探測范圍、精度和數據處理運行時間來合理選擇去噪方法。文中通過與國外測風激光雷達進行外場對比實驗，驗證了文中提出的載噪比算法的正確性，同時還對比分析驗證了該算法對兩種不同的去噪方法均適用。

在用載噪比初步判別信號有效區間的基礎上，本文通過實驗數據對比分析驗證了文中提出的信號有效性算法的正確性，同時該算法能進一步縮小信號的有效區間，從而提高了信號的有效精度，進而能大大提高其反演出的風場信息的可靠性。

[1] 周小林,孫東松,鐘志慶,等.多普勒測風激光雷達研究進展[J].大氣與環境光學學報,2007，2(3):161-168.

[2] TARG R,STEAKLEY B C,HAWLEY J G,et al.Coherent lidar airborne wind sensor II: flight-test results at 2 and 10 μm[J].Applied Optics,1996,35(36):7117-7127.

[3] Diao Weifeng,Zhang Xin,Liu Jiqiao,et al.All fiber pulsed coherent lidar development for wind profiles measurements in boundary layers[J].China Optical Letter,2014,12(7):072801.

[4] 王旭，何佩琨，毛二可.一種脈沖多普勒雷達數字信號處理機的設計[J].電子技術應用，2008，34(2):39-41.

[5] 陳勇，姚新宇，潘玉林,等.雷達實時仿真中的脈沖壓縮技術研究[J].電子技術應用,2012,38(1):118-121.

[6] Shangguan Mingjia,Xia Haiyun,Wang Chong,et al.All-fiber upconversion high spectral resolution wind lidar using a Fabry-Perot interferometer[J].Optics Express,2016,24(17):19322-19336.

[7] Yao Yuan, Li Feng, Chen Lian,et al.Note: a sub-sampling technique for frequency locking in Doppler wind lidar[J].Review of Scientific Instruments,2016,87(5):056107.

[8] 濮莉莉,周煜,孫劍鋒,等.激光雷達部分相干外差探測接收特性研究[J].光學學報,2011,31(12):260-268.

[9] 何平,李柏,吳蕾,等.確定風廓線雷達功率譜噪聲功率方法[J].應用氣象學報,2013,24(3):297-303.

[10] 步志超.相干測風激光雷達系統設計及數據處理算法研究[D].北京:北京理工大學,2014.

[11] 劉立人.合成孔徑激光成像雷達(VI)：時空散斑效應和外差探測信噪比[J].光學學報,2009,29(8):2326-2332.

Research on the algorithms of carrier to noise ratio and effectiveness of wind lidar signal

Luo Tian1, Hu Fangren1, Zhu Hailong2

(1. School of Optoelectronic Engineering, Nangjing University of Posts and Telecommunications, Nanjing 210046, China;2. Advanced Solid-State Laser Technology Research and Development Center, Nanjing Institute of Advanced Laser Technology, Nanjing 210038, China)

The carrier to noise ratio and the effectiveness of signal are important index for the performance of wind lidar. In our paper, we propose two new algorithms to analyze the carrier to noise ratio and the effectiveness of the signal, respectively. The out-field comparative experiments results indicate that the proposed carrier to noise ratio algorithm is correct and it is feasible to different denoising methods. Also, the experimental results show that there are good consistencies between the carrier to noise ratio and the effectiveness algorithms in judging the effectiveness of the signal. Meanwhile, compared with the carrier to noise ratio algorithm, the new effectiveness algorithm can reduce the effective range of the lidar signal, thus improving the accuracy of the effectiveness of the signal and wind information.

wind lidar; carrier to noise ratio; effectiveness; signal analysis

TN956

A

10.19358/j.issn.1674- 7720.2017.12.003

羅甜，胡芳仁，朱海龍.測風激光雷達信號載噪比和有效性算法研究[J].微型機與應用，2017,36(12)：8-12.

2016-12-26)

羅甜(1992-)，通信作者，女，碩士研究生，主要研究方向：測風激光雷達信號處理。E-mail:18795958495@163.com。

胡芳仁(1973-)，男，博士，教授，主要研究方向：微納米尺度的新型光電集成器件。

朱海龍(1988-)，男，碩士研究生，高級工程師，主要研究方向：測風激光雷達系統。