非相干多普勒测风激光雷达鉴频算法

非相干多普勒测风激光雷达鉴频算法

刘延文;孙学金;张传亮;李绍辉

【摘要】风场对大气的运动状态和运动趋势具有很强的代表性,所以风场的测量精度对数值天气预报和气候研究都至关重要.激光雷达具有很高的时空分辨率,在近几十年发展迅速,在对地观测中的作用越来越大,尤其在大气风场的探测中得到了很大的应用.本文从探测原理、分类和技术等方面分别对多普勒测风激光雷达进行了介绍,并着重总结了非相干多普勒测风激光雷达的多普勒频移算法,可以为我国星载激光雷达的研究提供参考.

【期刊名称】《激光与红外》

【年(卷),期】2018(048)010

【总页数】10页(P1204-1213)

【关键词】风;激光雷达;多普勒频移;边缘技术;条纹成像技术

【作者】刘延文;孙学金;张传亮;李绍辉

【作者单位】国防科技大学气象海洋学院,江苏南京211101;国防科技大学气象海洋学院,江苏南京211101;国防科技大学气象海洋学院,江苏南京211101;国防科技大学气象海洋学院,江苏南京211101

【正文语种】中文

【中图分类】TN958.98;X831

1 引言

风对于大气能量循环、污染物扩散、水汽和气溶胶粒子的输送都有重要的影响,更是大气环流的根本动力,是数值天气预报和气候研究中最迫切需要的数据,并且局地风场在飞机起飞与着陆、火箭发射及军事等方面具具有重要的意义[1-4],但在当前全球观测系统中,对风的观测并不充分,和气压、温度、湿度等大气基本变量相比,风的观测较为薄弱[5]。当今世界气象组织主要通过无线电探空网、机载仪器和气象雷达来对风场进行探测,但无线电探空网分布不均匀,在南半球和洋面上观测数据有限[6],边界层以上的风速可以根据地转理论和气压测量值计算得到[8],但该方法只适用于中高纬地区,所以需要直接探测资料来对大气流动有更准确的认识。激光雷达作为近几十年来快速发展的探测仪器,被用于温度[8-9]、气溶胶浓度和光学厚度[10]的测量,对风场的测量也有时空分辨率和探测精度高、探测范围大、响应速度快的优点,星载激光雷达更是可以获得高精度的全球风场。多普勒测风激光雷达通过测量随风场运动的分子和气溶胶粒子对激光造成的多普勒频移来实现风场的测量,分子和气溶胶的后向散射光被激光雷达接收,通过干涉仪后光子被电荷耦合器(CCD)接收,利用CCD上光子的累积和分布计算得到多普勒频移,利用多普勒频移就可以得到风速,所以多普勒频移的计算精度决定了风速测量精度。对于边缘技术和条纹成像技术,都有多种多普勒频移的检测算法,尤其对于条纹成像技术而言,不同算法的计算精度、计算速度和计算稳定性都不同,所以需要根据不同的要求选择不同的算法计算多普勒频移。本文总结了非相干测风激光雷达的不同算法,介绍个各种算法的适用性,比较了不同算法的精度和误差。

2 多普勒测风激光雷达的探测原理

多普勒测风激光雷达发射激光脉冲并接收大气后向散射信号,大气分子和气溶胶粒子随风场运动会与激光产生多普勒效应,从而改变散射回波的频率,通过测量和计算得到多普勒频移,就可以计算出径向风速,其表达式为:

vd=2vLOS/λ

(1)

式中,vd为多普勒频移;vLOS为径向风速;λ为波长。从式(1)可以看出,多普勒频移

与目标物的运动速度及激光雷达的工作波长有关,当激光雷达的载荷确定后,多普勒

频移只与目标(分子或气溶胶)的运动速度有关。根据径向风速和激光的天顶角(地基)或天底角(空基和天基)便可计算出水平径向风速,其表达式为:

vHLOS=vLOS/sinθ

(2)

其中,vHLOS为水平径向速度;θ为天顶角或天底角。

3 激光雷达的鉴频方法

按照鉴频方式不同,测风激光雷达可以分为非相干激光雷达[11-13]和相干激光雷达[14-15]。相干技术首先利用光电二极管将光信号转化为电信号[16],然后出射信号与回波信号混频,拍频信号的频率就是回波信号的多普勒频移,主要用于大气中低层Mie散射回波信号的频移测量[17]；非相干技术利用光学鉴频器将频移转化为功率、强度或功率空间分布变化,通过测量功率、强度或功率空间分布的变化来反演风速[18],探测范围可到大气中高层。

相干雷达的探测灵敏度、探测精度及信噪比高[19],噪声功率小,对太阳背景光不敏

感[20],但对光学准直性要求高,接收视场角的失配会严重影响激光雷达的性能,对于短波长的多普勒激光雷达校准要求更为严格,且无法准确探测Rayleigh散射信号。非相干探测对光学系统和激光性能的要求相对较低,系统结构简单,技术成熟且容易

实现,采用多脉冲累积可以提高信噪比和探测精度,减小激光散斑的影响,对Rayleigh散射和Mie散射信号都可以进行探测,利用短波长的激光可以提高Rayleigh散射强度,在气溶胶浓度为零的情况下也可以工作[21]。

非相干探测包括边缘技术和条纹成像技术。边缘技术是利用滤波器将后向散射回波信号的频移转化为能量变化,通过测量能量变化计算多普勒频移；条纹技术则是利

用不同波长的光透过干涉仪后干涉条纹位置不同的原理来实现多普勒频移的测量,

可采用的干涉仪有F-P标准具和Fizeau干涉仪等。F-P标准具的干涉条纹呈环状,且条纹的粗细和间隔不一致,风速与条纹半径变化为非线性关系[22],需要阵列式探

测器动态追踪条纹变化,和线形像素探测器CCD不匹配,不利于测量[23-24],利用环状阴极倍增管可以解决此问题,但其量子效率较低[25],另一种解决方法是将环形条

纹转变为线条纹,再用CCD探测[26]。Fizeau干涉仪的干涉条纹呈线形,可以直接

通过CCD检测条纹重心的变化来实现多普勒频移的测量[23],这种鉴频系统光路简单,系统效率高,干涉仪口径要求低[27]。条纹技术一般用于测量Mie散射信号,而边缘技术对Rayleigh散射信号和Mie散射信号都可以测量,在敏感性上两种方法并

没有显著的不同[21],风速测量精度十分接近[28]。

3.1 边缘技术

边缘技术可利用F-P标准具、M-Z干涉仪[29-31]、Michelson干涉仪、光栅、各种原子和分子滤波器[32],如钠、钾、碘[33]、银蒸汽滤波器等,探测灵敏度依赖于

分子与气溶胶的后向散射率及风速大小,F-P标准具是非相干探测的主要鉴频器[34-37],其透过率曲线具有陡峭的边缘,入射光频率变化会使透射光的的强度发生明显变化。

3.1.1 单边缘技术

由于标准具透过率曲线半宽点对应频率处的斜率最大,所以出射激光的频率υE锁定在此处,微小的频移将使标准具的透过率发生明显变化。由于布朗运动,Rayleigh散射和Mie散射光谱的谱宽有很大差别,但都可用高斯函数来表示:

(3)

式中,σx为散射信号光谱的标准差,下标“x”表示不同散射类型,“R”表示Rayleigh散射,“M”表示Mie散射。分子散射光谱的标准差σR=(8kT/mλ2)1/2；