高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

机动型激光测风雷达倾斜风场修正算法研究

李策 刘俊伟 赵培娥 周杰 谢日华 罗雄 周鼎富

引用本文:
Citation:

机动型激光测风雷达倾斜风场修正算法研究

    作者简介: 李策(1990-), 男, 硕士研究生, 现主要从事激光雷达技术的研究.
    通讯作者: 周鼎富, df_xh_zhou@sina.com
  • 中图分类号: TN958.98

Correction method of tilt wind field of mobile wind lidar

    Corresponding author: ZHOU Dingfu, df_xh_zhou@sina.com
  • CLC number: TN958.98

  • 摘要: 为了解决激光测风雷达处于倾斜放置的情况下测出的风场信息与实际风场之间存在的误差,分析了激光雷达的测量原理、放置情况以及存在的问题,采用空间坐标变换的方法对雷达倾斜状态进行理论计算并给出了修正计算公式。对比了水平放置和横向倾斜7.2°放置的激光雷达的实测数据,发现修正前由于倾斜放置导致雷达所测垂直气流数据间产生的误差较大,加入修正算法进行仿真后,数据曲线吻合良好,修正效果明显。结果表明,空间坐标变换的修正算法可使激光雷达在处于非特定位置时仍能正确地反映实际风场信息,从而验证了该修正算法的正确性和可行性。
  • Figure 1.  Theory of VAD scanning

    Figure 2.  Theory of DBS scanning

    Figure 3.  Contrast before and after lidar tilt

    Figure 4.  Schematic diagram of lidar tilt

    Figure 5.  Transformation of wind field between the horizontal plane and the inclined plane

    Figure 6.  Wind speed decomposition in the inclined plane

    Figure 7.  Wind speed decomposition in the horizontal plane

    Figure 8.  Coordinate system after origin coincidence

    Figure 9.  Comparison between wind speeds measured with No.1 and No.2 lidars on 100m level

    Figure 10.  Comparison between wind directions measured with No.1 and No.2 lidars on 100m level

    Figure 11.  Comparison between vertical air flow measured with No.1 and No.2 lidars on 50m level

    Figure 12.  Comparison between vertical air flow measured with No.1 and No.2 lidars on 100m level

    Figure 13.  Comparison between vertical air flow measured with No.1 and No.2 lidars on 150m level

    Figure 14.  Compare between vertical air-flow measured with No.1 and No.2 lidars on 200m level

  • [1]

    LAI D, CHEN Y, ZHOU D F, et al. Beam scanning of lidar and the simulation of the improved VAD inversion methods[J]. Laser Technology, 2008, 32(6): 584-586(in Chinese). 
    [2]

    LHERMITTE R M, ATLAS D. Precipitation motion by pulse Doppler[C]//American Meteorological Society. 9th Weather Radar Confe-rence.New York, USA: American Meteor Society, 1961: 218-223.
    [3]

    BROWNING K A, WEXLER R. The determination of kinematic properties of a wind field using Doppler radar[ J]. Journal of Applied Meteorology, 1968, 7(1):105-113. doi: 10.1175/1520-0450(1968)007<0105:TDOKPO>2.0.CO;2
    [4]

    LI L, XIE Y F, DONG G Y, et al. Wind field inversion technique for scanning wind lidar[J]. Chinese Optics, 2013, 6(2):251-258(in Chinese). 
    [5]

    NEWSOM R K, LIGON D, CALHOON R, et al. Retrieval of micro scale wind and temperature fields from single-and dual-Doppler lidar data[J]. Journal of Applied Meteorology and Climatology, 2005, 44(9): 1324-1344. doi: 10.1175/JAM2280.1
    [6]

    BANAKH V A, WERNER C. Retrieval of the turbulence parameters from spatial statistics of Doppler estimates[J]. Proceedings of the SPIE, 2004, 5575: 55-66. doi: 10.1117/12.566458
    [7]

    JIANG H Y, GE R Sh. A new retrieval technique for single Doppler radar[J]. Quarterly Journal of Applied Meteorology, 1997, 8(2):219-223. 
    [8]

    QIU C J, XU Q. A simple adjoint method of wind analysis for single-Doppler data[J]. Journal of Atmospheric and Oceanic Technology, 1992, 9(5): 588-598. doi: 10.1175/1520-0426(1992)009<0588:ASAMOW>2.0.CO;2
    [9]

    WAN R, TANG D Z, ZHANG P, et al. VAD eiementary analysis of nonlinear wind field[J]. Scientia Meteorologica Sinica, 2003, 23(3):314-324(in Chinese).
    [10]

    BAI J, ZHOU X B, WANG L K, et al. Error analysis of doppler radar expanded VAD technique[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2001, 37(1):48-54(in Chinese). 
    [11]

    REN P, WU Y L, ZHOU D F, et al. Numerical optimization in VAD invertion technique for wind lidar[J]. Laser Technology, 2009, 33(6):664-669(in Chinese). 
    [12]

    XUE S Q, YANG H P, GAO Y C. Comparison between VAD retrieved wind and radiosonde data[J].Scientia Meteorologica Sinica, 2010, 30(6):856-861(in Chinese). 
    [13]

    XU M Q, WANG B, XU H Z. 3D-radar altimetry definition analysis and error correction[J].Modern Defence Technology, 2013, 41(2):123-127(in Chinese). 
  • [1] 薛灏张全李鹏 . 半导体激光器远场光强特性分析. 激光技术, 2004, 28(6): 661-663.
    [2] 董吉辉胡企铨 . 空间碎片的激光雷达探测. 激光技术, 2007, 31(2): 185-187,191.
    [3] 亢甲杰张福民曲兴华 . 激光雷达坐标测量系统的测角误差分析. 激光技术, 2016, 40(6): 834-839. doi: 10.7510/jgjs.issn.1001-3806.2016.06.013
    [4] 王玉诏 . 一种用于激光雷达信号的自适应分段平滑算法. 激光技术, 2014, 38(5): 599-602. doi: 10.7510/jgjs.issn.1001-3806.2014.05.005
    [5] 陈文革黄铁侠卢益民 . 机载海洋激光雷达发展综述. 激光技术, 1998, 22(3): 147-152.
    [6] 刘长久杨华军赖燔 . 空间交会激光雷达信息测量技术. 激光技术, 2006, 30(6): 608-610,613.
    [7] 马泳冀航梁琨林宏 . 载波调制激光雷达技术在海洋探测中的应用. 激光技术, 2008, 32(4): 346-349.
    [8] 吕潮峰杜正春 . 基于激光雷达测量技术的集装箱卡车定位系统. 激光技术, 2007, 31(6): 596-599.
    [9] 汪少林曹开法胡顺星魏合理 . 对激光雷达几何因子的分析与测量. 激光技术, 2008, 32(2): 147-150.
    [10] 陶小红胡以华雷武虎蔡晓春 . 经验模式分解用于激光雷达大气回波处理. 激光技术, 2008, 32(6): 590-592,595.
    [11] 张小红杨敏郭亮 . 合成孔径激光雷达研究进展. 激光技术, 2011, 35(2): 255-259. doi: 10.3969/j.issn.1001-3806.2011.02.031
    [12] 唐磊吴海滨孙东松舒志峰 . 瑞利多普勒激光雷达F-P标准具的设计与校准分析. 激光技术, 2014, 38(2): 283-288. doi: 10.7510/jgjs.issn.1001-3806.2014.02.029
    [13] 王杰陈亚峰刘秋武杨杰黄见胡顺星 . 测污激光雷达探测气溶胶消光特性和Ångström指数. 激光技术, 2018, 42(6): 727-732. doi: 10.7510/jgjs.issn.1001-3806.2018.06.001
    [14] 吴俊杰王耀辉徐足音任佳莉张博义 . 基于多普勒激光雷达的机场边界层高度研究. 激光技术, 2023, 47(6): 778-785. doi: 10.7510/jgjs.issn.1001-3806.2023.06.008
    [15] 梁晓峰杨泽后王顺艳陈涌陈春利李晓锋李晶周鼎富 . 基于差分吸收激光雷达有毒有害气体遥测进展. 激光技术, 2021, 45(1): 53-60. doi: 10.7510/jgjs.issn.1001-3806.2021.01.010
    [16] 阳琴陈孝林曾诚徐诗月杨俊敏张志坚杨峰高剑波 . 基于DPMZM的微波光子倍频激光雷达仿真分析. 激光技术, 2023, 47(6): 729-735. doi: 10.7510/jgjs.issn.1001-3806.2023.06.001
    [17] 王春勇谢俊卞保民李振华 . 提高激光雷达测距精度的最佳信号门限比. 激光技术, 2007, 31(4): 408-411.
    [18] 李兵杨赟秀李潇呙长冬寇先果孔繁林袁鎏郑博仁 . APD阵列及其成像激光雷达系统的研究进展. 激光技术, 2023, 47(3): 310-316. doi: 10.7510/jgjs.issn.1001-3806.2023.03.004
    [19] 张春光张玉钧韩道文刘文清陈臻懿 . 机动车颗粒物的激光雷达监测. 激光技术, 2009, 33(2): 130-133.
    [20] 张飞飞夏海云孙东松 . 1.55μm相干激光雷达系统的硬目标速度探测. 激光技术, 2012, 36(5): 602-606. doi: 10.3969/j.issn.1001-3806.2012.05.007
  • 加载中
图(14)
计量
  • 文章访问数:  5401
  • HTML全文浏览量:  3851
  • PDF下载量:  17
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-07-28
  • 录用日期:  2016-09-30
  • 刊出日期:  2017-05-25

机动型激光测风雷达倾斜风场修正算法研究

    通讯作者: 周鼎富, df_xh_zhou@sina.com
    作者简介: 李策(1990-), 男, 硕士研究生, 现主要从事激光雷达技术的研究
  • 1. 西南技术物理研究所, 成都 610041
  • 2. 中国人民解放军 93169部队, 四平 136000

摘要: 为了解决激光测风雷达处于倾斜放置的情况下测出的风场信息与实际风场之间存在的误差,分析了激光雷达的测量原理、放置情况以及存在的问题,采用空间坐标变换的方法对雷达倾斜状态进行理论计算并给出了修正计算公式。对比了水平放置和横向倾斜7.2°放置的激光雷达的实测数据,发现修正前由于倾斜放置导致雷达所测垂直气流数据间产生的误差较大,加入修正算法进行仿真后,数据曲线吻合良好,修正效果明显。结果表明,空间坐标变换的修正算法可使激光雷达在处于非特定位置时仍能正确地反映实际风场信息,从而验证了该修正算法的正确性和可行性。

English Abstract

    • 风是大气气象重要要素之一,对人类的日常工作及工农业生产、军事活动等都有重要的影响,对风场的准确测量,是气象观测的一项重要工作。传统风场测量仪器有微波气象雷达、声雷达、探空仪等设备,而激光相干测风雷达与传统测风仪器相比,具有高精度、高时空分辨率等特点,被广泛地应用于低空大气风场的测量。

      激光测风雷达的工作原理是将激光发射到目标空域,沿光束传播途径上运动的气溶胶粒子会对激光脉冲产生后向散射效应和多普勒频移。测风雷达通过接收带有多普勒频移的散射光回波,对其进行光电转换和信号处理,并从中检测出多普勒频移,最终推算出激光视线方向上的风矢量信息。而要得到所要求的3维风场信息,则需再对风场进行反演运算。因此,径向风速的反演算法尤为重要。在风场反演方法中,速度方位显示法[1](velocity azimuth dispaly,VAD)和多普勒光束摆动(Dopplor beam swinging,DBS)四波束合成解算法因为运算量小且计算简单等特点被广泛应用。

      激光多普勒测风雷达探测到的风矢量为径向速度分量,当仰角固定时,径向速度在以雷达为中心的等距离圈上随方位大致呈简谐曲线分布,简谐曲线的振幅和位相反映了雷达扫描范围内的风速和风向,因此,要求激光测风雷达工作时须严格水平安装。另一方面,激光测风雷达因其体积小、机动性强等特点,在车载或其它机动性设备中被广泛应用,因此其机动性受地形因素影响很大,导致雷达在使用过程中难以保持水平,这样测出的风场数据则与真实值之间存在误差。本文中通过引入空间几何学的方法,利用空间坐标变换,对雷达倾斜后的风场误差进行修正。

    • 由于激光测风雷达直接探测得到的径向风矢量为沿扫描光束方向的风场分量。因此,需要选取适当的方法对径向风速进行反演。风场反演算法在国内外被普遍应用,在20世纪60年代速度方位显示方法由LHERMITTE提出,之后由BROWNING等人将其完善[3]。到20世纪70年代,在局部地区均匀风的条件下有人提出了速度体积处理技术,但此方法需要进行体扫,运算量大,分辨率也只有10m。到21世纪,由NEWSOM和BANAKH等人提出并改进了DBS四波束合成计算的方法[4-6],此方法至今仍被广泛运用。而在国内风场反演方法发展也十分迅速,JIANG[7]和QIU[8]等人分别提出了涡度-散度方法和简单伴随函数方法,对风场反演技术的发展提供了支持。

      现阶段风场反演算法以速度方位显示和DBS四波束合成计算两种算法居多,两种算法均基于假定风场稳定的情况下对风场进行多波束扫描,再通过反演计算得到高精度的风速风向信息,并计算出垂直气流信息。

    • VAD风场反演算法的基本原理是假定在均匀风和线性风的条件下,利用雷达在某一固定仰角上作360°的全方位圆周扫描,获取不同方位上的径向风速值后, 把一系列的径向风速拟合为方位角的余弦函数,最后求得扫描区域内的平均风场。所拟合的3个参量代表 3维风场信息[9-12]。激光光束的扫描方式如图 1所示。

      Figure 1.  Theory of VAD scanning

      设风矢量为(u, v, w),在雷达直角坐标系中:u沿着x轴,v沿着y轴,w沿着z轴,圆锥扫描半角γ是激光束指向与z正轴的夹角,扫描方位角θx坐标正轴为0°起点。

      选取正东方向为x轴正向,正北方向为y轴正向,测得的径向风速为vr,水平风速为vh和垂直气流为w,则有:

      $ {v_{\rm{r}}} = {v_{\rm{h}}}\sin \gamma + w\cos \gamma $

      (1)

      式中,γ为圆锥扫描半角。

      根据扫描角θ,水平风速vh可由水平分量uw合成,则得到公式:

      $ {v_{\rm{h}}} = u\cos \theta + v\sin \theta $

      (2)

      将(2)式带入(1)式,可得:

      $ {v_{\rm{r}}} = \left( {u\cos \theta + v\sin \theta } \right)\sin \gamma + w\cos \gamma $

      (3)

      令$a = \sqrt {{u^2} + {v^2}} $,由(3)式得:

      $ {v_{\rm{r}}} = b\cos \left( {a - \theta } \right) + c $

      (4)

      式中,$b = \sqrt {{u^2} + {v^2}} \sin \gamma , c = w\cos \gamma $。

      综上推理,VAD扫描方式中,每一个高度层都对应一个余弦曲线,余弦曲线的振幅、相位、偏距等参量均包含了水平风速、风向以及垂直气流等信息。

    • DBS原理即用四波束合成计算方法进行风场反演。DBS四波束扫描过程中,假设风场矢量保持不变,如图 2所示。设风矢量为(uvw),在雷达直角坐标系中:u沿着x轴,v沿着y轴,w沿着z轴,圆锥扫描半角γ是激光束指向与z正轴的夹角,扫描方位角θx坐标正轴为0°起点。

      Figure 2.  Theory of DBS scanning

      径向风速与风向矢量的关系可表示为:

      $ {v_{\rm{r}}} = u\cos \theta \sin \gamma + v\sin \theta \sin \gamma + w\cos \gamma $

      (5)

      设4个方向测得的径向风速分别为vr, 1vr, 2vr, 3vr, 4,对应的方位角分别为θ1θ2θ3θ4,则有:

      $ \left\{ \begin{array}{l} {v_{{\rm{r,1}}}} = u\cos {\theta _1}\sin \gamma + v\sin {\theta _1}\sin \gamma + w\cos \gamma \\ {v_{{\rm{r,2}}}} = u\cos {\theta _2}\sin \gamma + v\sin {\theta _2}\sin \gamma + w\cos \gamma \\ {v_{{\rm{r,3}}}} = u\cos {\theta _3}\sin \gamma + v\sin {\theta _3}\sin \gamma + w\cos \gamma \\ {v_{{\rm{r,4}}}} = u\cos {\theta _4}\sin \gamma + v\sin {\theta _4}\sin \gamma + w\cos \gamma \end{array} \right. $

      (6)

      对于激光雷达,4个扫描方位角与N,E,S,W这4个方向对准,即θ1θ2θ3θ4,分别为0°,90°,180°,270°。

      则(6)式可简化为:

      $ \left\{ \begin{array}{l} {v_{{\rm{r,E}}}} = {v_{{\rm{r,1}}}} = u\sin \gamma + w\cos \gamma \\ {v_{{\rm{r,N}}}} = {v_{{\rm{r,2}}}} = v\sin \gamma + w\cos \gamma \\ {v_{{\rm{r,W}}}} = {v_{{\rm{r,3}}}} = - u\sin \gamma + w\cos \gamma \\ {v_{{\rm{r,S}}}} = {v_{{\rm{r,4}}}} = - v\sin \gamma + w\cos \gamma \end{array} \right. $

      (7)

      进而推导出:

      $ \left\{ \begin{array}{l} u = \frac{{{v_{{\rm{r,E}}}} - {v_{{\rm{r,W}}}}}}{{2\sin \gamma }}\\ v = \frac{{{v_{{\rm{r,N}}}} - {v_{{\rm{r,S}}}}}}{{2\sin \gamma }}\\ w = \frac{{{v_{{\rm{r,E}}}} + {v_{{\rm{r,N}}}} + {v_{{\rm{r,W}}}} + {v_{{\rm{r,S}}}}}}{{4\cos \gamma }} \end{array} \right. $

      (8)
    • 可以看出,VAD和DBS原理均在雷达水平放置的情况下来进行风速测量和风场反演。而当雷达随其放置载体(例如车载雷达)经过各种地形时,会遇到上坡、下坡、转弯等各类状况,使得测风雷达不能保证时刻处于水平状态进行测量,导致风场反演结果与真实情况不符。如图 3所示,其中箭头为假定水平风和垂直气流。

      Figure 3.  Contrast before and after lidar tilt

    • 假定雷达在车体上放置时与车体底盘平行,则由车体倾斜而引起的雷达相对于水平地面的倾斜方向与车体一致。雷达倾斜示意图如图 4所示。

      Figure 4.  Schematic diagram of lidar tilt

      当车体沿轴向方向倾斜,用β表示,β>0°表示朝前倾斜,β < 0°表示朝后倾斜;车体垂直于轴向方向倾斜(横向倾斜),用α表示,α>0°表示朝右边倾斜,α < 0°表示朝左边倾斜。

      由于雷达倾斜的角度与车体倾斜角度相一致,因此雷达倾斜情况同样为:横向倾斜α,轴向倾斜为β,由此相关计算如下。

      图 4所示,设立方体ABCD-A1B1C1D为雷达,AD方向为雷达轴向,ABCD所处平面为水平地面,雷达D处为车头一侧,车体经横向倾斜α,轴向倾斜β以后,使雷达到达GHFC位置,于是:∠ABE=∠DCF=α,∠EFH=∠BCG=β

      设立方体边长为1,由于BRFC为平行四边形,所以EFAA1,又由于平面BRFC绕着轴CF旋转β角到达位置GHFC,易知GHFC为平行四边形。

      延长HFAD的延长线交于点P,连接CP,过点HHKCP,连接AK,则∠AKH为平面GHFC与平面ABCD的二面角。

      因为正方体ABCD-A1B1C1D棱长为1,所以有AE=ABtanα=tanα,经计算,可求出:

      $ \left\{ \begin{array}{l} AH = \tan \alpha + \tan \beta \\ AK = \frac{{\tan \alpha + \tan \beta }}{{\sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } }} \end{array} \right. $

      (9)

      因此,车体倾斜的角度为:

      $ \angle AKH = \arctan \sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } $

      (10)

      同理可得,倾斜方向与车体轴向的夹角为:

      $ \angle PAK = \arccos \frac{{\tan \beta }}{{\sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } }} $

      (11)

      由于雷达倾斜具有一定的方向性,通过分析可知:

      (1) 当α≥0°, β>0°时:

      $ \left\{ \begin{array}{l} \angle AKH = \arctan \sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } \\ \angle PAK = \arccos \frac{{\tan \beta }}{{\sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } }} \end{array} \right. $

      (12)

      (2) 当α≤0°, β>0°时:

      $ \left\{ \begin{array}{l} \angle AKH = \arctan \sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } \\ \angle PAK = 2{\rm{ \mathsf{ π} }} - \arccos \frac{{\tan \beta }}{{\sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } }} \end{array} \right. $

      (13)

      (3) 当α≤0°, β<0°时:

      $ \left\{ \begin{array}{l} \angle AKH = \arctan \sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } \\ \angle PAK = 2{\rm{ \mathsf{ π} }} - \arccos \frac{{\tan \beta }}{{\sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } }} \end{array} \right. $

      (14)

      (4) 当α≥0°, β<0°时:

      $ \left\{ \begin{array}{l} \angle AKH = \arctan \sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } \\ \angle PAK = \arccos \frac{{\tan \beta }}{{\sqrt {{{\tan }^2}\alpha + {{\tan }^2}\beta } }} \end{array} \right. $

      (15)

      (5) 当α=0°, β=0°时:

      $ \left\{ \begin{array}{l} \angle AKH = 0\\ \angle PAK = 0 \end{array} \right. $

      (16)
    • 当雷达与正北方向成γ角,同时自身倾斜角度为θ时,设原坐标系中的风矢量为(u, v, w),倾斜后新坐标系中的风矢量为(vx, vy, vz),${v_{\rm{h}}} = \sqrt {v_x^2 + v_y^2} $为倾斜以后测得的水平风速值的大小,ξ为倾斜后测得的风向的值,vz为倾斜后测得的垂直气流,如图 5所示。A为倾斜前的水平面,B为倾斜以后的面,OP为倾斜方向,OPOPA平面内的投影。

      Figure 5.  Transformation of wind field between the horizontal plane and the inclined plane

      vx-O-vyu-O′-v都为地理坐标系(原点与雷达重合,vy轴和v轴指向正北方向,vx轴和u轴指向正东方向,vz轴和w轴指向垂直向下方向)。

      则在vx-O-vy平面内,以OPx轴,OQy轴建立新的直角坐标系,将风矢量为[vx, vy, vz]重新分解,如图 6所示。

      Figure 6.  Wind speed decomposition in the inclined plane

      (1) 当ξγ时,风矢量与倾斜方向的夹角为π-(ξ-γ),水平风速在OP轴和OQ轴上的分量分别为:

      $ \left\{ \begin{array}{l} OP = {v_{\rm{h}}}\cos \left( {{\rm{ \mathsf{ π} }} + \gamma - \xi } \right) = - {v_{\rm{h}}}\cos \left( {\gamma - \xi } \right)\\ OQ = {v_{\rm{h}}}\sin \left( {{\rm{ \mathsf{ π} }} + \gamma - \xi } \right) = {v_{\rm{h}}}\sin \left( {\xi - \gamma } \right) \end{array} \right. $

      (17)

      (2) 同理,当ξ < γ时,水平风速在OP轴和OQ轴上的分量分别为:

      $ \left\{ \begin{array}{l} OP = {v_{\rm{h}}}\cos \left( {{\rm{ \mathsf{ π} }} + \gamma - \xi } \right) = - {v_{\rm{h}}}\cos \left( {\gamma - \xi } \right)\\ OQ = {v_{\rm{h}}}\sin \left( {{\rm{ \mathsf{ π} }} + \gamma - \xi } \right) = {v_{\rm{h}}}\sin \left( {\xi - \gamma } \right) \end{array} \right. $

      (18)

      可见,无论ξγ的关系如何,水平风速在OP轴和OQ轴上的分量计算结果相同,即:

      $ \left\{ \begin{array}{l} OP = - {v_{\rm{h}}}\cos \left( {\xi - \gamma } \right)\\ OQ = {v_{\rm{h}}}\sin \left( {\gamma - \xi } \right) \end{array} \right. $

      (19)

      同理,在u-O′-v平面内,如图 7所示。以OPx轴、OQ′为y轴建立新的直角坐标系,则由图 4可知,∠OPO′=θOPv轴正向的夹角为γ

      Figure 7.  Wind speed decomposition in the horizontal plane

      经计算可得:

      $ \left\{ \begin{array}{l} u = O'P\sin \gamma - O'Q'\cos \gamma \\ v = O'P\cos \gamma + O'Q'\sin \gamma \end{array} \right. $

      (20)

      图 7中坐标系POQ-vz与坐标系POQ′-w相互叠加,使原点重合,得两坐标系的关系如图 8所示。其中,M代表的是原来的OO′,MS代表OPMK代表OPMD代表OQ′和OQ,可得到:

      $ \left\{ \begin{array}{l} MS = MK\cos \theta - {v_z}\sin \theta \\ MD = MD\\ w = MK\sin \theta + {v_z}\cos \theta \end{array} \right. $

      (21)

      Figure 8.  Coordinate system after origin coincidence

      进而由(20)式和(21)式可推导出:

      $ \left\{ \begin{array}{l} u = - {v_{\rm{h}}}\cos \left( {\xi - \gamma } \right)\cos \theta \sin \gamma - {v_z}\sin \theta \cdot \\ \;\;\;\;\sin \gamma - {v_{\rm{h}}}\sin \left( {\xi - \gamma } \right)\cos \gamma \\ v = - {v_{\rm{h}}}\cos \left( {\xi - \gamma } \right)\cos \theta \cos \gamma - {v_z}\sin \theta \cdot \\ \;\;\;\;\cos \gamma + {v_{\rm{h}}}\sin \left( {\xi - \gamma } \right)\sin \gamma \\ w = - {v_{\rm{h}}}\cos \left( {\xi - \gamma } \right)\sin \theta + {v_z}\cos \theta \end{array} \right. $

      (22)

      对于风向计算,设倾斜前风向为ξxk为中间变量,则相关计算如下:

      $ k = \left\{ \begin{array}{l} \alpha \tan \left( {v{u^{ - 1}}} \right),\left( {u \ge 0,v > 0} \right)\\ \alpha \tan \left( {v{u^{ - 1}}} \right) + {\rm{ \mathsf{ π} }},\left( {u \le 0,v > 0} \right)\\ \alpha \tan \left( {v{u^{ - 1}}} \right) + {\rm{ \mathsf{ π} }},\left( {u \le 0,v < 0} \right)\\ \alpha \tan \left( {v{u^{ - 1}}} \right) + 2{\rm{ \mathsf{ π} }},\left( {u \ge 0,v < 0} \right)\\ 0,\left( {u = 0,v = 0} \right) \end{array} \right. $

      (23)

      从而推导出:

      $ {\xi _x} = \left\{ \begin{array}{l} 3{\rm{ \mathsf{ π} }}/2 - k,\left( {0 \le k \le {\rm{ \mathsf{ π} }}/2} \right)\\ 3{\rm{ \mathsf{ π} }}/2 - k,\left( {{\rm{ \mathsf{ π} }}/2 \le k \le {\rm{ \mathsf{ π} }}} \right)\\ 3{\rm{ \mathsf{ π} }}/2 - k,\left( {{\rm{ \mathsf{ π} }} \le k \le 3{\rm{ \mathsf{ π} /2}}} \right)\\ 7{\rm{ \mathsf{ π} }}/2 - k,\left( {3{\rm{ \mathsf{ π} }}/2 \le k \le 2{\rm{ \mathsf{ π} }}} \right) \end{array} \right. $

      (24)

      最后通过ξx×180/π将风向由弧度换算为角度便得出风向信息。

      在雷达随车体横向倾斜α、而轴向倾斜β后,雷达沿着与车头方向夹角为∠PAK的方向倾斜∠AKH。如果已知雷达所在车体轴向与正北方向的夹角为δ的话,则倾斜方向与正北方向的夹角为γ=δ+∠PAK,(vh, ξ, vz)分别为倾斜后测得的水平风速、风向、垂直气流,则可推出真实的水平风速为:

      $ {V_x} = \sqrt {{u^2} + {v^2}} $

      (25)

      修正后垂直气流(即倾斜以前的垂直气流)为:

      $ \begin{array}{*{20}{c}} {w = {v_z}\cos \angle AKH - {v_{\rm{h}}}\cos \left( {\xi - \gamma } \right) \cdot }\\ {\sin \angle AKH = {v_z}\cos \angle AKH - }\\ {{v_{\rm{h}}}\cos \left( {\xi - \delta - \angle PAK} \right)\sin \angle AKH} \end{array} $

      (26)

      在雷达倾斜的情况下(倾斜角度在10°以内),对于垂直气流测量的影响较大,而对水平风速风向的影响并不大[13](将在下一节的实验例证中运用数据对比曲线来证明),面对此类问题,作者提出运用空间坐标变换原理,将雷达的倾斜状况进行整体修正,计算量小且修正效果明显。

    • 为了验证该方法的正确性和有效性,本文中通过对实际风场的测风试验对该方法进行验证。在试验中,选用两台规格及参量相同的激光测风雷达在几个典型高度层上进行对比试验,试验中大气环境条件较好,能见度大于20km,气温0℃~10℃,湿度较低,试验主要内容如下:(1)1号雷达水平放置,2号雷达倾斜放置(横向倾斜7.2°),两台雷达均不添加修正算法,同时工作,输出数据并绘制曲线,证明倾斜后的垂直气流相差很大,而水平风速、风向相差不大; (2)为了达到实时对比,将前面得到的2号雷达的数据提取出来,添加修正算法,将修正后数据与1号雷达数据进行对比,证明修正算法起到了很好的修正作用,其中,2号雷达的放置情况为横向倾斜-7.2°。

      首先,提取100m高度层的风速、风向数据,并绘制曲线进行对比说明,如图 9图 10所示。

      Figure 9.  Comparison between wind speeds measured with No.1 and No.2 lidars on 100m level

      Figure 10.  Comparison between wind directions measured with No.1 and No.2 lidars on 100m level

      图 9图 10中可以看出,在修正前,由于倾斜放置而导致的风速、风向测量误差并不大。而由修正前后的曲线对比可以看出,修正后的风速、风向与修正前非常接近,相当于从理论上证明了雷达倾斜放置情况(10°以内)对于风速、风向的测量结果影响较小。

      其次,提取50m,100m,150m及200m共4个典型高度层的垂直气流数据进行修正前与修正后的实验对比,并根据数据绘制出曲线,如图 11~图 14所示。

      Figure 11.  Comparison between vertical air flow measured with No.1 and No.2 lidars on 50m level

      Figure 12.  Comparison between vertical air flow measured with No.1 and No.2 lidars on 100m level

      Figure 13.  Comparison between vertical air flow measured with No.1 and No.2 lidars on 150m level

      Figure 14.  Compare between vertical air-flow measured with No.1 and No.2 lidars on 200m level

      从图中可以看出,修正前由于倾斜放置导致雷达所测垂直气流数据间产生的误差较大,加入修正算法进行仿真后,数据曲线吻合良好,修正效果明显。

    • 经过实验对比,首先得出在雷达处于倾斜放置状态下进行测量时,对于风速及风向的影响程度较小,而对垂直气流的影响程度很大。其次,在进行垂直气流的误差修正过程中可以看到,当雷达未加修正算法时,1号与2号雷达垂直气流曲线差异较大,而在加入修正算法后,两台雷达所显示的垂直气流曲线相当吻合,修正效果十分明显。从而证明了该算法对于垂直气流的修正补偿是有效可行的。值得注意的是,本算法的验证条件是基于倾斜角度在10°以内、风场相对稳定(均匀风和线性风场)的条件下进行修正的,而当风场波动较大或倾斜角度过大时,该算法的有效性还有待继续考证。

参考文献 (13)

目录

    /

    返回文章
    返回