CN102565771B - 一种基于单站风廓线雷达的风切变识别和追踪方法 - Google Patents
一种基于单站风廓线雷达的风切变识别和追踪方法 Download PDFInfo
- Publication number
- CN102565771B CN102565771B CN 201010586775 CN201010586775A CN102565771B CN 102565771 B CN102565771 B CN 102565771B CN 201010586775 CN201010586775 CN 201010586775 CN 201010586775 A CN201010586775 A CN 201010586775A CN 102565771 B CN102565771 B CN 102565771B
- Authority
- CN
- China
- Prior art keywords
- wind
- shear
- data
- wind shear
- module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于单站风廓线雷达的风切变识别和追踪方法,其具体步骤为:第一步,搭建风廓线雷达风切变追踪系统,包括数据采集模块、数据预处理模块、风切变识别模块和风切变追踪模块;第二步,数据采集模块采集风廓线雷达探测数据;第三步,数据预处理模块对风廓线雷达探测数据进行质量控制;第四步,风切变识别模进行风切变的自动探测和识别;第五步,风切变追踪模块实现风切变的自动追踪和报警;至此,实现了基于单站风廓线雷达的风切变识别和追踪方法。本方法具有简单易行、识别准确率高、全自动运行等优点,可以获取水平风的垂直切变和垂直风的垂直切变信息,有效探测追踪站点上空风切变出现的起始时间、高度和强度。
Description
技术领域
本发明涉及一种风切变识别和追踪方法,特别是一种基于单站风廓线雷达的风切变识别和追踪方法。
背景技术
风廓线雷达是以大气湍流为探测目标、具有晴空探测能力的脉冲多普勒测速雷达,主要用于测量风随高度的变化。由于风廓线雷达采用遥感方式,可连续无人值守运行,且探测资料具有种类多、时空分辨率高、廓线形式等众多优点,在高空(非地面)风探测方面正发挥着越来越重要的作用。大量的观测资料表明,对飞行安全有影响的切变线在风廓线雷达连续的观测资料中有着清晰的表现。因此利用风廓线雷达实现对风切变的探测具有重要的意义。
目前常用的地基风切变探测方法一般是通过包括低空风切变预警系统(LLWAS)、多普勒天气雷达、多普勒激光雷达和集成低空风切变警报系统实现的,由于这些系统与风廓线雷达探测原理和方式上的差异,所用的风切变识别和追踪方法并不适用于风廓线雷达。
发明内容
本发明的目的在于提供一种基于单站风廓线雷达的风切变识别和追踪方法,解决过去的风切变识别和追踪方法不适用于风廓线雷达的问题。
一种基于单站风廓线雷达的风切变识别和追踪方法,具体步骤为:
第一步搭建风廓线雷达风切变识别和追踪系统
风廓线雷达风切变识别和追踪系统包括:数据采集模块、数据预处理模块、风切变识别模块和风切变追踪模块。数据采集模块用于采集风廓线雷达探测数据,数据预处理模块用于对风廓线雷达探测数据进行质量控制,风切变识别模块用于风切变的自动探测和识别,风切变追踪模块用于风切变的自动追踪和报警。
第二步数据采集模块采集风廓线雷达探测数据
数据采集模块采集风廓线雷达的探测数据作为风切变探测、识别、追踪的初始输入信息。由于风切变的生命周期短,数据采集模块输出的风廓线时间分辨率为10min。
第三步数据预处理模块对风廓线雷达探测数据进行质量控制
数据预处理模块识别并剔除探测数据中由于外界干扰等原因产生的虚假数据,避免在风切变计算时产生错误信息,提供虚假情报。
当风向小于0°或大于360°,或风速值为负值时,判断为无物理意义的数据;
当低空风速大于80m/s时,判断为不符合一般统计特征的数据;
当上下高度层及前一个时刻的同一高度层都为无效数据时,判断为孤立的数据;
对于不满足时空连续性的数据,其判断准则为:
首先对数据进行中值滤波。将水平风分解为u、v分量,用于检测的数据记为(ui、vi),与其空间和时间上紧邻的数据之中值(um、vm)进行比较,如果二者间的差超过指定门限(Tu、Tv),则将被检测的风数据(ui、vi)标识为无效。门限的确定公式为:
Tu=max(Tu1,T2) (1)
Tv=max(Tv1,T2) (2)
Tu1=0.2|um+ui| (3)
Tv1=0.2|vm+vi| (4)
T2=0.67(-6.127×10-8h2+0.0012h+7.383) (5)
式中,h为距地面高度。
然后对数据进行水平风的垂直切变检查。在边界层,除了极端天气状况,水平风的垂直切变小于0.2s-1,更强的风切变对应着更为极端的天气状况。在此,设定最大水平风的垂直切变门限为0.2s-1,超过门限的情况下,风廓线雷达探测所依赖的局地均匀各向同性假设得不到满足,则提供的数据不能反映真实大气情况。
第四步风切变识别模块进行风切变的自动探测和识别
风切变是大气中空间两点之间每单位距离内风速和风向的变化,为一个向量,单位为m/s/30m或者s-1,公式为:
Vshear=(V2-V1)/d21 (6)
其中,V2和V1为空间两点上的风矢量,d21为空间两点间的距离。
将风切变分解为水平风的水平切变、水平风的垂直切变、垂直风的水平切变和垂直风的垂直切变。单站风廓线雷达探测的是站点上空的垂直风廓线,用于探测水平风的垂直切变和垂直风的垂直切变。
将水平风分解为沿跑道的分量和垂直于跑道的分量,分别记为纵风和横风,分解公式为:
Vx=Vcos(D+180°-A) (7)
Vy=Vsin(D+180°-A) (8)
其中,Vx为纵风,代表沿跑道的风分量,Vy为横风,代表垂直于跑道的风分量,D为风向,A为航向。相应的,沿跑道的风切变和垂直于跑道的风切变的计算公式为:
Vshear_x=(V2x-V1x)/d21 (9)
Vshear_y=(V2y-V1y)/d21 (10)
垂直风的垂直切变计算公式为:
Vshear_z=(V2z-V1z)/d21 (11)
另外,对风切变无方向性要求时,水平风的垂直切变的计算公式为:
第五步风切变追踪模块实现风切变的自动追踪和报警
对于当前时刻风廓线,计算其每个高度层的风切变强度,当没有超过门限的风切变高度层时,判断其前一时刻是否有超过门限的风切变高度层,没有则无动作,有则解除上次追踪;当前时刻风廓线中有超过门限的风切变高度层时,判断其前一时刻是否有超过门限的风切变高度层,没有则开始一次新的追踪并记录风切变所发生的时间、高度层及其强度,有则根据风切变所发生的高度层判断是否为同一切变线,是则继续追踪,不是则开始新的追踪。在风切变追踪开始和解除时,自动给出报警信号。
至此,实现了基于单站风廓线雷达的风切变识别和追踪方法。
本方法简单易行、识别准确率高、全自动运行,可以获取水平风的垂直切变和垂直风的垂直切变信息,有效探测追踪站点上空风切变出现的起始时间、高度和强度。
具体实施方式
一种基于单站风廓线雷达的风切变识别和追踪方法,具体步骤为:
第一步搭建风廓线雷达风切变识别和追踪系统
风廓线雷达风切变识别和追踪系统包括:数据采集模块、数据预处理模块、风切变识别模块和风切变追踪模块。数据采集模块用于采集风廓线雷达探测数据,数据预处理模块用于对风廓线雷达探测数据进行质量控制,风切变识别模块用于风切变的自动探测和识别,风切变追踪模块用于风切变的自动追踪和报警。
第二步数据采集模块采集风廓线雷达探测数据
数据采集模块采集风廓线雷达的探测数据作为风切变探测、识别、追踪的初始输入信息。由于风切变的生命周期短,数据采集模块输出的风廓线时间分辨率为10min。
第三步数据预处理模块对风廓线雷达探测数据进行质量控制
数据预处理模块识别并剔除探测数据中由于外界干扰等原因产生的虚假数据,避免在风切变计算时产生错误信息,提供虚假情报。
当风向小于0°或大于360°,或风速值为负值时,判断为无物理意义的数据;
当低空风速大于80m/s时,判断为不符合一般统计特征的数据;
当上下高度层及前一个时刻的同一高度层都为无效数据时,判断为孤立的数据;
对于不满足时空连续性的数据,其判断准则为:
首先对数据进行中值滤波。将水平风分解为u、v分量,用于检测的数据记为(ui、vi),与其空间和时间上紧邻的数据之中值(um、vm)进行比较,如果二者间的差超过指定门限(Tu、Tv),则将被检测的风数据(ui、vi)标识为无效。门限的确定公式为:
Tu=max(Tu1,T2) (1)
Tv=max(Tv1,T2) (2)
Tu1=0.2|um+ui| (3)
Tv1=0.2|vm+vi| (4)
T2=0.67(-6.127×10-8h2+0.0012h+7.383) (5)
式中,h为距地面高度。
然后对数据进行水平风的垂直切变检查。在边界层,除了极端天气状况,水平风的垂直切变小于0.2s-1,更强的风切变对应着更为极端的天气状况。在此,设定最大水平风的垂直切变门限为0.2s-1,超过门限的情况下,风廓线雷达探测所依赖的局地均匀各向同性假设得不到满足,则提供的数据不能反映真实大气情况。
第四步风切变识别模块进行风切变的自动探测和识别
风切变是大气中空间两点之间每单位距离内风速和风向的变化,为一个向量,单位为m/s/30m或者s-1,公式为:
Vshear=(V2-V1)/d21 (6)
其中,V2和V1为空间两点上的风矢量,d21为空间两点间的距离。
将风切变分解为水平风的水平切变、水平风的垂直切变、垂直风的水平切变和垂直风的垂直切变。单站风廓线雷达探测的是站点上空的垂直风廓线,用于探测水平风的垂直切变和垂直风的垂直切变。
将水平风分解为沿跑道的分量和垂直于跑道的分量,分别记为纵风和横风,分解公式为:
Vx=Vcos(D+180°-A) (7)
Vy=Vsin(D+180°-A) (8)
其中,Vx为纵风,代表沿跑道的风分量,Vy为横风,代表垂直于跑道的风分量,D为风向,A为航向。相应的,沿跑道的风切变和垂直于跑道的风切变的计算公式为:
Vshear_x=(V2x-V1x)/d21 (9)
Vshear_y=(V2y-V1y)/d21 (10)
垂直风的垂直切变计算公式为:
Vshear_z=(V2z-V1z)/d21 (11)
另外,对风切变无方向性要求时,水平风的垂直切变的计算公式为:
第五步风切变追踪模块实现风切变的自动追踪和报警
对于当前时刻风廓线,计算其每个高度层的风切变强度,当没有超过门限的风切变高度层时,判断其前一时刻是否有超过门限的风切变高度层,没有则无动作,有则解除上次追踪;当前时刻风廓线中有超过门限的风切变高度层时,判断其前一时刻是否有超过门限的风切变高度层,没有则开始一次新的追踪并记录风切变所发生的时间、高度层及其强度,有则根据风切变所发生的高度层判断是否为同一切变线,是则继续追踪,不是则开始新的追踪。在风切变追踪开始和解除时,自动给出报警信号。
至此,实现了基于单站风廓线雷达的风切变识别和追踪方法。
Claims (1)
1.一种基于单站风廓线雷达的风切变识别和追踪方法,其特征在于具体步骤为:
第一步搭建风廓线雷达风切变识别和追踪系统
风廓线雷达风切变识别和追踪系统包括:数据采集模块、数据预处理模块、风切变识别模块和风切变追踪模块;数据采集模块用于采集风廓线雷达探测数据,数据预处理模块用于对风廓线雷达探测数据进行质量控制,风切变识别模块用于风切变的自动探测和识别,风切变追踪模块用于风切变的自动追踪和报警;
第二步数据采集模块采集风廓线雷达探测数据
数据采集模块采集风廓线雷达的探测数据作为风切变探测、识别、追踪的初始输入信息;由于风切变的生命周期短,数据采集模块输出的风廓线时间分辨率为10min;
第三步数据预处理模块对风廓线雷达探测数据进行质量控制
数据预处理模块识别并剔除探测数据中由于外界干扰原因产生的虚假数据,避免在风切变计算时产生错误信息,提供虚假情报;
当风向小于0°或大于360°,或风速值为负值时,判断为无物理意义的数据;
当低空风速大于80m/s时,判断为不符合一般统计特征的数据;
当上下高度层及前一个时刻的同一高度层都为无效数据时,判断为孤立的数据;
对于不满足时空连续性的数据,其判断准则为:
首先对数据进行中值滤波;将水平风分解为u、v分量,用于检测的数据记为(ui、vi),与其空间和时间上紧邻的数据之中值(um、vm)进行比较,如果二者间的差超过指定门限(Tu、Tv),则将被检测的风数据(ui、vi)标识为无效;门限的确定公式为:
Tu=max(Tu1,T2) (1)
Tv=max(Tv1,T2) (2)
Tu1=0.2|um+ui| (3)
Tv1=0.2|vm+vi| (4)
T2=0.67(-6.127×10-8h2+0.0012h+7.383) (5)
式中,h为距地面高度;
然后对数据进行水平风的垂直切变检查;在边界层,除了极端天气状况,水平风的垂直切变小于0.2s-1,更强的风切变对应着更为极端的天气状况;在此,设定最大水平风的垂直切变门限为0.2s-1,超过门限的情况下,风廓线雷达探测所依赖的局地均匀各向同性假设得不到满足,则提供的数据不能反映真实大气情况;
第四步风切变识别模块进行风切变的自动探测和识别
风切变是大气中空间两点之间每单位距离内风速和风向的变化,为一个向量,单位为m/s/30m或者s-1,公式为:
Vshear=(V2-V1)/d21 (6)
其中,V2和V1为空间两点上的风矢量,d21为空间两点间的距离;
将风切变分解为水平风的水平切变、水平风的垂直切变、垂直风的水平切变和垂直风的垂直切变;单站风廓线雷达探测的是站点上空的垂直风廓线,用于探测水平风的垂直切变和垂直风的垂直切变;
将水平风分解为沿跑道的分量和垂直于跑道的分量,分别记为纵风和横风,分解公式为:
Vx=Vcos(D+180°-A) (7)
Vy=Vsin(D+180°-A) (8)
其中,Vx为纵风,代表沿跑道的风分量,Vy为横风,代表垂直于跑道的风分量,D为风向,A为航向;相应的,沿跑道的风切变和垂直于跑道的风切变的计算公式为:
Vshear_x=(V2x-V1x)/d21 (9)
Vshear_y=(V2y-V1y)/d21 (10)
垂直风的垂直切变计算公式为:
Vshear_z=(V2z-V1z)/d21 (11)
另外,对风切变无方向性要求时,水平风的垂直切变的计算公式为:
第五步风切变追踪模块实现风切变的自动追踪和报警
对于当前时刻风廓线,计算其每个高度层的风切变强度,当没有超过门限的风切变高度层时,判断其前一时刻是否有超过门限的风切变高度层,没有则无动作,有则解除上次追踪;当前时刻风廓线中有超过门限的风切变高度层时,判断其前一时刻是否有超过门限的风切变高度层,没有则开始一次新的追踪并记录风切变所发生的时间、高度层及其强度,有则根据风切变所发生的高度层判断是否为同一切变线,是则继续追踪,不是则开始新的追踪;在风切变追踪开始和解除时,自动给出报警信号;
至此,实现了基于单站风廓线雷达的风切变识别和追踪方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010586775 CN102565771B (zh) | 2010-12-14 | 2010-12-14 | 一种基于单站风廓线雷达的风切变识别和追踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010586775 CN102565771B (zh) | 2010-12-14 | 2010-12-14 | 一种基于单站风廓线雷达的风切变识别和追踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102565771A CN102565771A (zh) | 2012-07-11 |
CN102565771B true CN102565771B (zh) | 2013-08-28 |
Family
ID=46411660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010586775 Active CN102565771B (zh) | 2010-12-14 | 2010-12-14 | 一种基于单站风廓线雷达的风切变识别和追踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102565771B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983975B (zh) * | 2014-06-05 | 2017-01-11 | 中国气象局气象探测中心 | 基于两种雷达的大气运动垂直速度检测方法及系统 |
CN105607063B (zh) * | 2016-01-05 | 2017-12-15 | 北京无线电测量研究所 | 一种机场低空风切变探测方法及系统 |
US10247170B2 (en) * | 2016-06-07 | 2019-04-02 | General Electric Company | System and method for controlling a dynamic system |
CN106772387B (zh) * | 2016-12-21 | 2019-06-14 | 中国航空工业集团公司雷华电子技术研究所 | 一种风切变识别方法 |
CN109324335B (zh) * | 2018-12-17 | 2023-10-31 | 北京无线电测量研究所 | 一种基于激光雷达识别风切变的方法和系统 |
CN111208534A (zh) * | 2020-01-20 | 2020-05-29 | 安徽四创电子股份有限公司 | 一种激光雷达与风廓线雷达联合探测识别风切变的方法 |
CN113486844B (zh) * | 2021-07-26 | 2022-04-15 | 中科三清科技有限公司 | 水平风切变位置判断方法、装置、电子设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251596A (zh) * | 2008-03-21 | 2008-08-27 | 北京理工大学 | 一种风廓线雷达非相参无线电干扰抑制方法 |
-
2010
- 2010-12-14 CN CN 201010586775 patent/CN102565771B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251596A (zh) * | 2008-03-21 | 2008-08-27 | 北京理工大学 | 一种风廓线雷达非相参无线电干扰抑制方法 |
Non-Patent Citations (2)
Title |
---|
胡明宝等.风廓线雷达中风切变分析方法的初步研究.《气象科学》.2010,第30卷(第4期),第510-514页. |
风廓线雷达中风切变分析方法的初步研究;胡明宝等;《气象科学》;20100831;第30卷(第4期);第510-514页 * |
Also Published As
Publication number | Publication date |
---|---|
CN102565771A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102565771B (zh) | 一种基于单站风廓线雷达的风切变识别和追踪方法 | |
CN102081801B (zh) | 多特征自适应融合船舶跟踪和航迹检测方法 | |
CN105467393B (zh) | 一种毫米波云雷达三维风场反演方法及系统 | |
CN103033808A (zh) | 基于雷达回波特征的铁路路障探测与告警方法及装置 | |
CN104345312B (zh) | 一种毫米波云雷达的数据融合方法及系统 | |
CN110806193A (zh) | 地铁隧道形变检测系统 | |
CN103177584B (zh) | 一种基于枚举试探的车速检测方法 | |
CN104407346A (zh) | 一种基于信息融合的移动式跑道异物监测方法 | |
CN216110866U (zh) | 一种基于rssi的隧道车辆事故检测及预警装置 | |
CN114419825A (zh) | 基于毫米波雷达和摄像头的高铁周界入侵监测装置及方法 | |
CN206312365U (zh) | 智能超高超宽检测仪 | |
CN207317740U (zh) | 一种车辆长度测量装置 | |
CN109615880A (zh) | 一种基于雷达图像处理的车流量测量方法 | |
CN101425128A (zh) | 基于激光传感器深度图像的人流检测方法 | |
CN113869196A (zh) | 一种基于激光点云数据多特征分析的车型分类方法及装置 | |
CN104199002B (zh) | 抑制风力发电旋叶虚警目标的雷达信息处理方法 | |
CN202995036U (zh) | 一种基于雷达回波特征的铁路路障探测与告警方法的系统 | |
CN115236673A (zh) | 一种大型车辆多雷达融合感知系统和方法 | |
CN206974384U (zh) | 一种车辆外廓的测量单元 | |
CN203501986U (zh) | 载重汽车宽高自动测量装置 | |
CN112147457B (zh) | 一种基于Hilbert-Huang变换的地下综合管廊电缆故障检测定位系统及方法 | |
CN102494633A (zh) | 用于防止船舶撞击桥梁的智能探测装置 | |
CN202432995U (zh) | 用于防止船舶撞击桥梁的智能探测装置 | |
Zhang et al. | Transmission line distance measurement with millimeter wave radar based on improved robust Kalman filter | |
CN114892556B (zh) | 一种新型智慧斑马线系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20191225 Address after: 214000 28 Weiming Road, Binhu District, Wuxi City, Jiangsu Province Patentee after: Aerospace new weather Technology Co., Ltd Address before: 100854 box 203, box 45, No. 142, Beijing Patentee before: No.23 Institute of No.20 Academy, CASIC |
|
TR01 | Transfer of patent right |