CN109583593A - A low-altitude wind shear identification method based on automatic weather station - Google Patents

A low-altitude wind shear identification method based on automatic weather station Download PDF

Info

Publication number
CN109583593A
CN109583593A CN201811285720.1A CN201811285720A CN109583593A CN 109583593 A CN109583593 A CN 109583593A CN 201811285720 A CN201811285720 A CN 201811285720A CN 109583593 A CN109583593 A CN 109583593A
Authority
CN
China
Prior art keywords
wind
data
shear
wind shear
automatic weather
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.)
Granted
Application number
CN201811285720.1A
Other languages
Chinese (zh)
Other versions
CN109583593B (en
Inventor
张靖
陈少应
刘淑昕
雷雅慧
余小强
刘志鹏
戚颖
仇逸菲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Create Electronics Co ltd
Original Assignee
Sun Create Electronics Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sun Create Electronics Co ltd filed Critical Sun Create Electronics Co ltd
Priority to CN201811285720.1A priority Critical patent/CN109583593B/en
Publication of CN109583593A publication Critical patent/CN109583593A/en
Application granted granted Critical
Publication of CN109583593B publication Critical patent/CN109583593B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00Indicating or recording presence, absence, or direction, of movement
    • G01P13/02Indicating direction only, e.g. by weather vane
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft

Landscapes

  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

The invention relates to a low-altitude wind shear identification method based on an automatic weather station, which comprises the following steps: acquiring historical data of each wind element in a monitoring area; performing quality control and time synchronization on historical data to obtain a data list; selecting effective data of actually generated wind shear, and calculating the shear value of each wind element; scaling a standard value according to the distance between every two stations of the automatic meteorological station, and comparing the standard value with a wind shear judgment threshold value to obtain a label for judging whether wind shear occurs; using each wind element as input, using a tag of whether wind shear occurs as a result to form a data set, and using a SMOTE algorithm to balance the number of the two tag types; training by using a machine learning method, and outputting a recognition model; and reading the latest data of each automatic weather station, and inputting the latest data into the identification model to obtain an identification result. The method has the advantages of high identification precision and good local adaptability, can effectively identify the horizontal wind shear near the ground of the runway, and improves the guarantee efficiency of meteorological service.

Description

A kind of low-level wind shear recognition methods based on automatic weather station
Technical field
The present invention relates to Data of Automatic Weather analysis fields, in particular relate to a kind of low latitude based on automatic weather station Wind shear recognition methods.
Background technique
Wind shear is a kind of common meteor, refer in sustained height or different height short distance wind direction (or) wind speed Quick variation.According to the difference of wind vector, wind shear can be divided into the horizontal shear of horizontal wind, horizontal vertical wind shear With vertical wind shear.It is referred to as low wherein at 600 meters of flying height or less, mainly in the wind shear of takeoff and landing stage generation Empty wind shear.It is counted according to World Meteorological Organization and International Civil Aviation Organization, low-level wind shear is to takeoff and landing stage flight safety Maximum weather phenomenon is threatened, many great air crashes are caused by low-level wind shear.And in aircraft landing, for resolution Wind field situation highly below is the most important thing of observation.
Currently, relying primarily on space wind field in relevant equipment search coverage for the detection identification of wind shear both at home and abroad Distribution and situation of change.Automatic weather station is generally deployed in runway two sides, close automatically horizontal nearby for intuitively obtaining runway Wind field, therefore in traditional wind shear recognition methods, it can use Data of Automatic Weather and calculate the horizontal shear of surface wind.Its In, more mature efficient system is the LLWAS system in the U.S., but its to ground survey wind devices deployment quantity and position have compared with High requirement, and it is insufficient to the specific aim of local landform, climatic characteristic, device characteristics.Current many medium and small airports are in construction Wind devices are surveyed on the ground for being not equipped with the quantity as required by LLWAS system and installation site, and LLWAS system wind shear algorithm is in local Adaptability it is also invalidated reliable, the judgement of forecaster is relied primarily in practical business, accuracy and coverage still have promotion Space.
Summary of the invention
According to problems of the prior art, the present invention provides a kind of, and the low-level wind shear based on automatic weather station is known Other method has the advantages that accuracy of identification is high, local adaptability is good, can carry out to runway horizontal wind shear near the ground effective Identification, improve the guarantee efficiency of Meteorological Services.
The invention adopts the following technical scheme:
A kind of low-level wind shear recognition methods based on automatic weather station, which comprises the steps of:
S1 obtains the historical data of all collected each wind elements of automatic weather station in monitoring region;
S2 carries out quality control to historical data with rejecting abnormalities data, obtains the valid data of each wind element, and according to The data moment of setting carries out time synchronization to valid data, obtains the corresponding data list of valid data;
S3 chooses the valid data in the time range that front and back is set at the time of actually occurring wind shear in data list, The shear at all automatic weather stations being calculated in the monitoring region corresponding each wind element each data moment between standing two-by-two Value;
S4, by it is described monitoring region in automatic weather station stand two-by-two between apart from bi-directional scaling standard value, will scale Standard value afterwards obtains whether each data moment occurs wind shear compared with wind shear discrimination threshold, and provides each Whether the data moment occurs the label of wind shear;
S5 is constituted as a result using each wind element of each time data as input, the label for whether occurring wind shear Data set, and SMOTE algorithm EDS extended data set is used to make two kinds of tag types quantity balance that wind shear whether occurs, it obtains Data set after balance;
S6 is trained the data set after balance using machine learning method, chooses the best training result of recall ratio, And corresponding model is exported, using this model as identification model;
S7 reads the latest data of each automatic weather station in monitoring region, using the method for step S2 to newest number in real time According to being pre-processed, and will treated that data input the identification model, and then obtain recognition result.
Preferably, in step S1, the wind element includes instantaneous wind speed, instantaneous wind direction, gustiness, fitful wind wind direction, two Minute mean wind speed, two minutes mean wind directions, two minutes maximum wind velocities, two minutes maximum wind directions, ten minutes mean wind speeds, very Clock mean wind direction, ten minutes maximum wind velocities, ten minutes maximum wind directions.
It is further preferred that carrying out quality control to historical data includes data normalization inspection processing, gas in step S2 Wait educational circles's limit value inspection processing, the processing of extreme value range check, internal consistency inspection processing, time consistency inspection processing;Institute It states data normalization inspection processing to refer to and delete the repeated data occurred in historical data, and corresponding to historical data Filename, the date, the time, automatic weather station site information carry out improve and uniform format;The climatology boundary value inspection Processing refers to the local historical climate feature according to locating for Meteorological Field relevant criterion and monitoring region, goes through to what can not be occurred History data are rejected;The extreme value range check processing refers to the history monthly average wind speed for choosing each automatic weather station acquisition For basic unit, calculate the standard deviation of the annual monthly average wind speed, using this month history maximum wind velocity plus twice of standard deviation as Maximum cuts twice of standard deviation as minimum using this month history minimum windspeed, and rejecting is more than that maximum value is less than minimum Historical data;The internal consistency inspection processing specifically refers to carry out inspection processing to the physical link each wind element, i.e., Then wind speed when according to ten minutes maximum wind velocity >=two minute maximum wind velocities, two minutes maximum wind velocity >=instantaneous wind speeds, wind direction being 0 ° These judgment criterias of≤0.2m/s, rejecting abnormalities data;Time consistency inspection processing specifically refer to it is a certain current when It is checked with the change rate of five minutes corresponding historical datas before this moment at quarter, upper limit threshold is greater than to change rate Corresponding historical data is rejected, while to a certain current time and one hour corresponding historical data before this moment Change rate checked, to change rate be less than lower threshold corresponding historical data reject, the upper limit threshold and Lower threshold local historical climate feature according to locating for monitoring region determines.
Still more preferably, in step S2, it is specific that time synchronization is carried out to valid data according to the data moment of setting Refer to the data list that valid data are established to valid data using per half a minute as the data moment, i.e., by each of valid data Entry collects in that the hithermost moment half a minute corresponding data list of its original moment.
Still more preferably, in step S3, front and back each 12 at the time of actually occurring wind shear in data list is chosen Hour in valid data, according in valid data instantaneous wind speed and instantaneous wind direction obtain instantaneous vector wind, according to fitful wind wind Speed and fitful wind wind direction obtain vector fitful wind, obtain two minutes mean vectors according to two minutes mean wind speeds and two minutes mean wind directions Wind obtains ten minutes average vector winds according to ten minutes mean wind speeds and ten minutes mean wind directions, calculates in the monitoring region All automatic weather stations stand two-by-two between the phasor difference of instantaneous vector wind, the phasor difference of vector fitful wind, two minutes average vector winds Phasor difference, ten minutes average vector winds phasor difference, and then obtain all automatic weather stations stand two-by-two between each data moment Shear value.
Still more preferably, in step S4, according to the international judgment criteria of wind shear, i.e. between two automatic weather stations Spacing when being 4km, when the wind shear value between two automatic weather stations is greater than 15 section, then be judged to that wind shear has occurred, International judgment criteria is floated downward 10% to 20%, in this, as wind shear discrimination threshold.
Still more preferably, in step S4, the distance between the automatic weather station in the monitoring region is stood two-by-two is pressed Scaling standard value, by the standard value after scaling compared with wind shear discrimination threshold, if the instantaneous wind at a certain data moment Shear value is greater than wind shear discrimination threshold or battle array wind shear value is greater than wind shear discrimination threshold, average wind shear in two minutes Value is greater than wind shear discrimination threshold and instantaneous wind shear value and wind shear discrimination threshold difference are within 20%, then when this data Have occurred wind shear quarter, it is on the contrary then wind shear does not occur;And whether wind shear is occurred to each data moment and provides corresponding mark Label, that is, have occurred the label that wind shear provides "Yes", the label that wind shear provides "No" do not occur.
Still more preferably, in step S5, two kinds of tag types quantity balance that wind shear whether occurs is specific Refer to that data that label is "Yes" and label are the data bulk difference of "No" less than 5%.
Still more preferably, step S6 specifically refers to that support vector machines, logic is respectively adopted using machine learning method It returns and three kinds of models of decision tree, the data set after input balances is trained, and fold into capable intersection using 10 in training process Verifying chooses the best training result of recall ratio, and exports corresponding model, using this model as identification model.
The advantages and beneficial effects of the present invention are:
1) present invention first acquisition monitoring region in automatic weather station each wind element historical data, to historical data into The control of row quality simultaneously carries out time synchronization to valid data according to the data moment of setting, obtains the corresponding data column of valid data Table;The valid data in the time range that front and back is set at the time of actually occurring wind shear in data list are chosen, are calculated The shear value at each wind element each data moment;Again by it is described monitoring region in automatic weather station stand two-by-two between distance by than Example scaling standard value, by the standard value after scaling compared with wind shear discrimination threshold, and then obtains whether each data moment is sent out Raw wind shear, and provide the label whether each data moment occur wind shear;Using each wind element of each time data as Input, the label for whether occurring wind shear constitute data set as a result, and use SMOTE algorithm EDS extended data set make whether Two kinds of tag types quantity balance of wind shear, the data set after being balanced occurs;Then, using machine learning method to flat Data set after weighing apparatus is trained, and chooses the best training result of recall ratio, and export corresponding model, using this model as knowledge Other model;Finally, the latest data of each automatic weather station in monitoring region, the method pair controlled using the quality are read in real time Latest data is pre-processed, and by treated, data input the identification model, and then obtain recognition result.The present invention adopts It is analyzed and is trained with based on Data of Automatic Weather, effective identification can be made to low latitude horizontal wind shear, compare traditional base In the wind shear identification technology of Data of Automatic Weather, the history number of each wind element is analyzed in innovation by the way of machine learning According to containing special to the monitoring local landform in region, climate characteristic, each site location feature and each measuring device data Point etc. be relatively difficult to directly give amendment or it is existing study considering for not yet specific influence factor, so that larger improve Low level wind The accuracy of shear identification, provides strong support for Meteorological safeguard and service.
Detailed description of the invention
Fig. 1 is the flow chart of method of the invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other Embodiment shall fall within the protection scope of the present invention.
As shown in Figure 1, a kind of low-level wind shear recognition methods based on automatic weather station, includes the following steps:
S1 obtains the historical data of all collected each wind elements of automatic weather station in monitoring region;
Specifically, the wind element includes instantaneous wind speed, instantaneous wind direction, gustiness, fitful wind wind direction, two minutes average winds Speed, two minutes mean wind directions, two minutes maximum wind velocities, two minutes maximum wind directions, ten minutes mean wind speeds, ten minutes average winds To, ten minutes maximum wind velocities, ten minutes maximum wind directions.
S2 carries out quality control to historical data with rejecting abnormalities data, obtains the valid data of each wind element, and according to The data moment of setting carries out time synchronization to valid data, obtains the corresponding data list of valid data;
Specifically, carrying out quality control to historical data includes data normalization inspection processing, climatology boundary value inspection Processing, the processing of extreme value range check, internal consistency inspection processing, time consistency inspection processing;The data normalization inspection Reason is investigated and prosecuted to refer to and deletes the repeated data occurred in historical data, and to the corresponding filename of historical data, the date, when Between, the site information of automatic weather station improve and uniform format;The climatology boundary value inspection processing refers to according to gas The local historical climate feature as locating for industry relevant criterion and monitoring region, picks the historical data that can not occur It removes;The extreme value range check processing refers to that choosing the history monthly average wind speed of each automatic weather station acquisition is basic unit, The standard deviation for calculating the annual monthly average wind speed, using this month history maximum wind velocity plus twice of standard deviation as maximum, with this The moon, history minimum windspeed cut twice of standard deviation as minimum, and rejecting is more than historical data of the maximum value less than minimum;Institute It states internal consistency inspection processing to specifically refer to carry out inspection processing to the physical link each wind element, i.e., most according to ten minutes When big wind speed >=two minute maximum wind velocity, two minutes maximum wind velocity >=instantaneous wind speeds, wind direction are 0 ° then wind speed≤0.2m/s these Judgment criteria, rejecting abnormalities data;Time consistency inspection processing specifically refer to a certain current time and this moment it The change rate of five minutes preceding corresponding historical datas is checked, the corresponding history number of upper limit threshold is greater than to change rate It is carried out according to being rejected, while to the change rate at a certain current time and one hour corresponding historical data before this moment Check, the corresponding historical data for being less than lower threshold to change rate is rejected, the upper limit threshold and lower threshold according to Local historical climate feature locating for region is monitored to determine.
Time synchronization is carried out to valid data according to the data moment of setting to specifically refer to valid data with per half a minute The data list of valid data is established for the data moment, i.e., by each entry in valid data collect its original moment near In that close moment half a minute corresponding data list.
S3 chooses the valid data in the time range that front and back is set at the time of actually occurring wind shear in data list, The shear at all automatic weather stations being calculated in the monitoring region corresponding each wind element each data moment between standing two-by-two Value;
Specifically, the valid data before and after choosing at the time of actually occurring wind shear in data list in each 12 hours, According in valid data instantaneous wind speed and instantaneous wind direction obtain instantaneous vector wind, sweared according to gustiness and fitful wind wind direction Fitful wind is measured, obtains two minutes average vector winds according to two minutes mean wind speeds and two minutes mean wind directions, it is average according to ten minutes Wind speed and ten minutes mean wind directions obtain ten minutes average vector winds, calculate all automatic weather stations in the monitoring region two-by-two The phasor difference of instantaneous vector wind between standing, the phasor difference of vector fitful wind, the phasor difference of two minutes average vector winds, ten minutes it is average The phasor difference of vector wind, so obtain all automatic weather stations stand two-by-two between each data moment shear value.
S4, by it is described monitoring region in automatic weather station stand two-by-two between apart from bi-directional scaling standard value, will scale Standard value afterwards obtains whether each data moment occurs wind shear compared with wind shear discrimination threshold, and provides each Whether the data moment occurs the label of wind shear;
Specifically, according to the international judgment criteria of wind shear, i.e. when spacing between two automatic weather stations is 4km, two When wind shear value between a automatic weather station is greater than 15 section, then it is judged to that wind shear has occurred, it will be under international judgment criteria Floating 10% to 20%, in this, as wind shear discrimination threshold.
By it is described monitoring region in automatic weather station stand two-by-two between apart from bi-directional scaling standard value, after scaling Standard value is compared with wind shear discrimination threshold, if the instantaneous wind shear value at a certain data moment is greater than wind shear discrimination threshold, Or battle array wind shear value is greater than wind shear discrimination threshold, two minutes average wind shear values are greater than wind shear discrimination threshold and instantaneous Wind shear value and wind shear discrimination threshold difference are within 20%, then this data moment has occurred wind shear, on the contrary then do not send out Raw wind shear;And whether wind shear is occurred to each data moment and provides corresponding label, that is, wind shear has occurred and provides "Yes" Label, the label that wind shear provides "No" does not occur.
S5 is constituted as a result using each wind element of each time data as input, the label for whether occurring wind shear Data set, and SMOTE algorithm EDS extended data set is used to make two kinds of tag types quantity balance that wind shear whether occurs, it obtains Data set after balance;
Specifically, two kinds of tag types quantity balance that wind shear whether occurs specifically refers to the number that label is "Yes" According to label be "No" data bulk difference less than 5%.
S6 is trained the data set after balance using machine learning method, chooses the best training result of recall ratio, And corresponding model is exported, using this model as identification model;
Specifically, three kinds of support vector machines, logistic regression and decision tree models is respectively adopted using machine learning method, Data set after input balance is trained, and folds into row cross validation using 10 in training process, and it is best to choose recall ratio Training result, and corresponding model is exported, using this model as identification model.
S7 reads the latest data of each automatic weather station in monitoring region, using the method for step S2 to newest number in real time According to being pre-processed, and will treated that data input the identification model, and then obtain recognition result.
Method of the invention is illustrated below with reference to embodiment.
8 automatic meteorological station equipments are now deployed with around Hefei City new bridge International airport runway, first to be located at runway north and south Illustrate the specific embodiment of the method for the present invention for the automatic weather station A and automatic weather station B at both ends.
It is about 2.7km at a distance between two station automatic weather station A and automatic weather station B, acquires automatic weather station A and automatic gas As the historical data that station two station B is built a station so far certainly, the method according to the invention first pre-processes historical data, rejects not Trust data, and historical data being synchronized in time, by data available naturalization to synchronization to utilize.
Concern report actually occurs 12 hours before and after wind shear moment periods, extracts the wind element letter wherein paid close attention to Breath, including instantaneous wind speed, wind direction, gustiness, wind direction, two minutes mean wind speeds, wind directions, ten minutes wind speed, wind directions.With the world The universal standard (4km range inscribe becomes larger in 15 sections (about 7.7m/s)) is zoomed in and out by distance, and floats downward 10%, i.e. threshold value W= 7.7*(2.7/4)*0.9≈4.6m/s.The shear value of each wind element at two stations is calculated, and using above-mentioned threshold value W as automatic meteorological Stand A and two station automatic weather station B each wind element shear discrimination standard.
During differentiation, in conjunction with forecaster's experience, instantaneous wind shear value, fitful wind are successively considered from big to small by weight Shear value, two minutes average wind shear values, the shear value of ten minutes average winds, instantaneous wind shear value sentences greater than wind shear Other threshold value or battle array wind shear value are greater than wind shear discrimination threshold, two minutes average wind shear values are greater than wind shear and differentiate threshold Value and instantaneously wind shear value and wind shear discrimination threshold difference, then labeled as wind shear occurs, otherwise mark within 20% For wind shear does not occur.Front and back 12 hours automatic weather station A and automatic gas at the time of actually occurring wind shear to report As the wind shear label of each wind element in station two station B and label constitutes data set.
Using SMOTE algorithm EDS extended data set, it will expand to and be labeled as labeled as the data entry that wind shear does not occur The data entry of wind shear occurs unanimously or relatively (the two quantity difference controls within 5%).
The data set expanded is inputted into support vector machines, three kinds of models of logistic regression and decision tree are trained, use 10 folding cross validations take recall ratio as preferential, output optimal result model, using this model as identification model.
The real time data for acquiring automatic weather station A and two station automatic weather station B, carries out the pre- of method same as described above to it After processing, trained identification model is inputted, and then obtain recognition result.
In conclusion the present invention provides a kind of low-level wind shear recognition methods based on automatic weather station, has and knows The advantage that other precision is high, local adaptability is good, can effectively identify runway horizontal wind shear near the ground, improve gas As the guarantee efficiency of service.

Claims (9)

1.一种基于自动气象站的低空风切变识别方法,其特征在于,包括如下步骤:1. a low-altitude wind shear identification method based on automatic weather station, is characterized in that, comprises the steps: S1,获取监测区域内所有自动气象站采集到的各风要素的历史数据;S1, obtain the historical data of each wind element collected by all automatic weather stations in the monitoring area; S2,对历史数据进行质量控制以剔除异常数据,得到各风要素的有效数据,并根据设定的数据时刻对有效数据进行时间同步,得到有效数据对应的数据列表;S2, perform quality control on the historical data to eliminate abnormal data, obtain valid data of each wind element, and perform time synchronization on the valid data according to the set data time to obtain a data list corresponding to the valid data; S3,选取数据列表中实际发生风切变的时刻前后设定的时间范围内的有效数据,计算得到该监测区域内的所有自动气象站两两站间对应的各风要素每个数据时刻的切变值;S3, select the valid data in the set time range before and after the time when the wind shear actually occurs in the data list, and calculate the shear value of each wind element corresponding to each data moment between all the automatic weather stations in the monitoring area. variable value; S4,将所述监测区域内的自动气象站两两站间的距离按比例缩放标准值,将缩放后的标准值与风切变判别阈值比较,进而得到每个数据时刻是否发生风切变,并给出每个数据时刻是否发生风切变的标签;S4, scaling the standard value of the distance between two automatic weather stations in the monitoring area in proportion, and comparing the scaled standard value with the wind shear judgment threshold, and then obtaining whether wind shear occurs at each data moment, And give the label of whether wind shear occurs at each data moment; S5,将每个时刻数据的各风要素作为输入、是否发生风切变的标签作为结果构成数据集,并使用SMOTE算法扩充数据集使得是否发生风切变的两种标签类型数量平衡,得到平衡后的数据集;S5, take each wind element of the data at each moment as the input and the label of whether wind shear occurs as the result to form a data set, and use the SMOTE algorithm to expand the data set to balance the number of the two types of labels whether the wind shear occurs, and the balance is obtained. the latter dataset; S6,利用机器学习方法对平衡后的数据集进行训练,选取查全率最好的训练结果,并输出相应的模型,以此模型作为识别模型;S6, use the machine learning method to train the balanced data set, select the training result with the best recall rate, and output the corresponding model, and use this model as the recognition model; S7,实时读取监测区域的各自动气象站的最新数据,采用步骤S2的方法对最新数据进行预处理,并将处理后的数据输入所述识别模型,进而得到识别结果。S7 , read the latest data of each automatic weather station in the monitoring area in real time, use the method of step S2 to preprocess the latest data, and input the processed data into the recognition model to obtain the recognition result. 2.根据权利要求1所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S1中,所述风要素包括瞬时风速、瞬时风向、阵风风速、阵风风向、两分钟平均风速、两分钟平均风向、两分钟最大风速、两分钟最大风向、十分钟平均风速、十分钟平均风向、十分钟最大风速、十分钟最大风向。2. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 1, is characterized in that: in step S1, described wind element comprises instantaneous wind speed, instantaneous wind direction, gust wind speed, gust wind direction, two minutes Average wind speed, 2-minute average wind direction, 2-minute maximum wind speed, 2-minute maximum wind direction, 10-minute average wind speed, 10-minute average wind direction, 10-minute maximum wind speed, and 10-minute maximum wind direction. 3.根据权利要求2所述的一种基于自动气象站的低空风切变识别方法,其特征在于;步骤S2中,对历史数据进行质量控制包括数据标准化检查处理、气候学界限值检查处理、极值范围检查处理、内部一致性检查处理、时间一致性检查处理;所述数据标准化检查处理是指对历史数据中出现的重复数据进行删除,并对历史数据相应的文件名、日期、时间、自动气象站的站点信息进行完善和格式统一;所述气候学界限值检查处理是指根据气象行业相关标准和监测区域所处的当地历史气候特征,对不可能出现的历史数据进行剔除;所述极值范围检查处理是指选取每个自动气象站采集的历史月平均风速为基本单位,计算每年该月平均风速的标准差,以该月历史最大风速加上两倍标准差作为极大值,以该月历史最小风速减掉两倍标准差作为极小值,剔除超过极大值、小于极小值的历史数据;所述内部一致性检查处理具体是指对各风要素间的物理联系进行检查处理,即根据十分钟最大风速≥两分钟最大风速、两分钟最大风速≥瞬时风速、风向为0°时则风速≤0.2m/s这些判断标准,剔除异常数据;所述时间一致性检查处理具体是指对某一当前时刻与此时刻之前的五分钟的相应的历史数据的变化率进行检查,对变化率大于上限阈值的相应的历史数据进行剔除,同时对某一当前时刻与此时刻之前的一小时的相应的历史数据的变化率进行检查,对变化率小于下限阈值的相应的历史数据进行剔除,所述上限阈值和下限阈值根据监测区域所处的当地历史气候特征确定。3. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 2, is characterized in that; in step S2, carry out quality control to historical data including data standardization check processing, climatology limit value check processing, Extremum range check processing, internal consistency check processing, and time consistency check processing; the data standardization check processing refers to deleting the duplicate data appearing in the historical data, and checking the corresponding file name, date, time, The site information of the automatic weather station is improved and the format is unified; the climatological limit value check processing refers to the elimination of impossible historical data according to the relevant standards of the meteorological industry and the local historical climate characteristics of the monitoring area; the The extreme value range check processing refers to selecting the historical monthly average wind speed collected by each automatic weather station as the basic unit, calculating the standard deviation of the monthly average wind speed each year, and taking the historical maximum wind speed of the month plus twice the standard deviation as the maximum value, Taking the historical minimum wind speed of the month minus twice the standard deviation as the minimum value, and excluding the historical data that exceeds the maximum value and is less than the minimum value; the internal consistency check process specifically refers to the physical connection between the wind elements. Inspection and processing, that is, according to the judgment criteria of the maximum wind speed in ten minutes ≥ the maximum wind speed in two minutes, the maximum wind speed in two minutes ≥ instantaneous wind speed, and the wind speed when the wind direction is 0°, the abnormal data is excluded. Specifically, it refers to checking the rate of change of the corresponding historical data at a certain current moment and five minutes before this moment, eliminating the corresponding historical data whose rate of change is greater than the upper limit threshold, and at the same time checking the rate of change between a certain current moment and five minutes before this moment. The rate of change of the corresponding historical data for one hour is checked, and the corresponding historical data whose rate of change is less than the lower threshold is eliminated. The upper threshold and the lower threshold are determined according to the local historical climate characteristics where the monitoring area is located. 4.根据权利要求3所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S2中,根据设定的数据时刻对有效数据进行时间同步具体是指将有效数据以每半分钟为数据时刻建立有效数据的数据列表,即将有效数据中的每个条目归集到其原始时刻最靠近的那个半分钟时刻对应的数据列表内。4. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 3, it is characterized in that: in step S2, according to the data moment of setting time synchronization to valid data specifically refers to valid data to A data list of valid data is established for the data moment every half-minute, that is, each entry in the valid data is collected into the data list corresponding to the half-minute moment closest to its original moment. 5.根据权利要求4所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S3中,选取数据列表中实际发生风切变的时刻前后各十二小时内的有效数据,根据有效数据中的瞬时风速和瞬时风向得到瞬时矢量风,根据阵风风速和阵风风向得到矢量阵风,根据两分钟平均风速和两分钟平均风向得到两分钟平均矢量风,根据十分钟平均风速和十分钟平均风向得到十分钟平均矢量风,计算该监测区域内的所有自动气象站两两站间的瞬时矢量风的矢量差、矢量阵风的矢量差、两分钟平均矢量风的矢量差、十分钟平均矢量风的矢量差,进而得到所有自动气象站两两站间每个数据时刻的切变值。5. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 4, is characterized in that: in step S3, selects the effective time within each twelve hours before and after the moment when wind shear actually occurs in the data list. data, obtain the instantaneous vector wind according to the instantaneous wind speed and instantaneous wind direction in the valid data, obtain the vector gust according to the gust wind speed and gust wind direction, obtain the two-minute average vector wind according to the two-minute average wind speed and two-minute average wind direction, and obtain the two-minute average vector wind according to the ten-minute average wind speed and The ten-minute average wind direction is obtained from the ten-minute average vector wind, and the vector difference of the instantaneous vector wind, the vector difference of the vector gust, the vector difference of the two-minute average vector wind, and the ten-minute average vector wind between all automatic weather stations in the monitoring area are calculated. The vector difference of the average vector wind, and then the shear value of each data moment between the two stations of all automatic weather stations is obtained. 6.根据权利要求5所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S4中,根据风切变的国际判断标准,即两个自动气象站之间的间距为4km时,两个自动气象站之间的风的切变值大于15节时,则判定为发生了风切变,将国际判断标准下浮10%至20%,以此作为风切变判别阈值。6. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 5 is characterized in that: in step S4, according to the international judgment standard of wind shear, namely the distance between two automatic weather stations When it is 4km and the wind shear value between the two automatic weather stations is greater than 15 knots, it is judged that wind shear has occurred, and the international judgment standard is lowered by 10% to 20% as the wind shear judgment threshold . 7.根据权利要求6所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S4中,将所述监测区域内的自动气象站两两站间的距离按比例缩放标准值,将缩放后的标准值与风切变判别阈值比较,若某一数据时刻的瞬时风的切变值大于风切变判别阈值,或阵风的切变值大于风切变判别阈值、两分钟平均风的切变值大于风切变判别阈值且瞬时风的切变值与风切变判别阈值差值在20%以内,则此数据时刻发生了风切变,反之则未发生风切变;并对每个数据时刻是否发生风切变给出相应的标签,即发生了风切变给出“是”的标签,未发生风切变给出“否”的标签。7. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 6, is characterized in that: in step S4, the distance between two automatic weather stations in described monitoring area is scaled proportionally Standard value, compare the scaled standard value with the wind shear discrimination threshold, if the instantaneous wind shear value at a certain data moment is greater than the wind shear discrimination threshold, or the gust shear value is greater than the wind shear discrimination threshold, two If the shear value of the minute average wind is greater than the wind shear judgment threshold, and the difference between the instantaneous wind shear value and the wind shear judgment threshold is within 20%, then the wind shear occurs at the moment of this data, otherwise there is no wind shear. ; and give a corresponding label to whether wind shear occurs at each data moment, that is, if wind shear occurs, a "Yes" label is given, and no wind shear occurs, a "No" label is given. 8.根据权利要求7所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S5中,所述是否发生风切变的两种标签类型数量平衡具体是指标签为“是”的数据与标签为“否”的数据数量差小于5%。8. A low-altitude wind shear identification method based on an automatic weather station according to claim 7, characterized in that: in step S5, the quantity balance of the two types of labels whether the wind shear occurs or not specifically refers to the label as The number of "yes" data differs from those labeled "no" by less than 5%. 9.根据权利要求8所述的一种基于自动气象站的低空风切变识别方法,其特征在于:步骤S6具体是指利用机器学习方法,分别采用支持向量机、逻辑回归和决策树三种模型,输入平衡后的数据集进行训练,且训练过程中采用10折进行交叉验证,选取查全率最好的训练结果,并输出相应的模型,以此模型作为识别模型。9. a kind of low-altitude wind shear identification method based on automatic weather station according to claim 8, is characterized in that: step S6 specifically refers to utilizing machine learning method, adopts three kinds of support vector machine, logistic regression and decision tree respectively Model, input the balanced data set for training, and use 10-fold cross-validation in the training process, select the training result with the best recall rate, and output the corresponding model, which is used as the recognition model.
CN201811285720.1A 2018-10-31 2018-10-31 Low-altitude wind shear identification method based on automatic meteorological station Active CN109583593B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811285720.1A CN109583593B (en) 2018-10-31 2018-10-31 Low-altitude wind shear identification method based on automatic meteorological station

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811285720.1A CN109583593B (en) 2018-10-31 2018-10-31 Low-altitude wind shear identification method based on automatic meteorological station

Publications (2)

Publication Number Publication Date
CN109583593A true CN109583593A (en) 2019-04-05
CN109583593B CN109583593B (en) 2020-11-13

Family

ID=65920911

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811285720.1A Active CN109583593B (en) 2018-10-31 2018-10-31 Low-altitude wind shear identification method based on automatic meteorological station

Country Status (1)

Country Link
CN (1) CN109583593B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113962462A (en) * 2021-10-25 2022-01-21 中国科学院空天信息创新研究院 Wind field stability prediction method and system based on convolutional neural network
CN115169236A (en) * 2022-07-22 2022-10-11 黄河水利委员会水文局 Automatic identification method for cut line based on combination of data mining and weather analysis
CN115169133A (en) * 2022-07-19 2022-10-11 中国消防救援学院 A kind of wind farm monitoring method and system
CN115830380A (en) * 2022-12-06 2023-03-21 中科三清科技有限公司 Wind shear line identification method and device based on artificial intelligence, storage medium and terminal
CN115861811A (en) * 2022-12-06 2023-03-28 中科三清科技有限公司 Wind shear region identification method and device, storage medium and terminal

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105607063A (en) * 2016-01-05 2016-05-25 北京无线电测量研究所 Detection method and system for low-altitude wind shear at airport
CN106772387A (en) * 2016-12-21 2017-05-31 中国航空工业集团公司雷华电子技术研究所 A kind of wind shear recognition methods

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105607063A (en) * 2016-01-05 2016-05-25 北京无线电测量研究所 Detection method and system for low-altitude wind shear at airport
CN106772387A (en) * 2016-12-21 2017-05-31 中国航空工业集团公司雷华电子技术研究所 A kind of wind shear recognition methods

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JAMES N. K. LIU ET AL.: "Chaotic Oscillatory-Based Neural Network for Wind Shear and Turbulence Forecast With LiDAR Data", 《IEEE》 *
蒋立辉: "基于仿真雷达图像的低空风切变类型识别研究", 《激光与红外》 *
迟继峰 等: "复杂地形多测风塔综合地貌及风切变拟合修正的风资源评估方法研究", 《华电技术》 *
郭忠立: "双流机场风切变预警系统方案设计及LLWAS算法研究", 《气象水文海洋仪器》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113962462A (en) * 2021-10-25 2022-01-21 中国科学院空天信息创新研究院 Wind field stability prediction method and system based on convolutional neural network
CN115169133A (en) * 2022-07-19 2022-10-11 中国消防救援学院 A kind of wind farm monitoring method and system
CN115169133B (en) * 2022-07-19 2025-05-30 中国消防救援学院 A wind farm monitoring method and system
CN115169236A (en) * 2022-07-22 2022-10-11 黄河水利委员会水文局 Automatic identification method for cut line based on combination of data mining and weather analysis
CN115830380A (en) * 2022-12-06 2023-03-21 中科三清科技有限公司 Wind shear line identification method and device based on artificial intelligence, storage medium and terminal
CN115861811A (en) * 2022-12-06 2023-03-28 中科三清科技有限公司 Wind shear region identification method and device, storage medium and terminal

Also Published As

Publication number Publication date
CN109583593B (en) 2020-11-13

Similar Documents

Publication Publication Date Title
CN109583593A (en) A low-altitude wind shear identification method based on automatic weather station
Rasmussen et al. Verification of the origins of rotation in tornadoes experiment: VORTEX
Hon Predicting low-level wind shear using 200-m-resolution NWP at the Hong Kong International Airport
CN115755220B (en) Airport gust forecasting and correcting method based on combination of numerical simulation and deep learning
Bramberger et al. Vertically propagating mountain waves—A hazard for high-flying aircraft?
Tafferner et al. ADWICE: Advanced diagnosis and warning system for aircraft icing environments
CN109407177A (en) Dense fog identifying system and methods for using them based on machine learning and conventional meteorological observation
Hagelin et al. Nowcasting with the AROME model: First results from the high-resolution AROME airport
Pauley et al. Assimilation of in-situ observations
Hu et al. Objective verification of Clear-Air Turbulence (CAT) diagnostic performance in China using in situ aircraft observation
Eilts et al. Oklahoma downbursts and their asymmetry
Díaz-Fernández et al. Mountain Waves Analysis in the Vicinity of the Madrid‐Barajas Airport Using the WRF Model
Gurke et al. The development of the wake vortices warning system for Frankfurt airport: Theory and implementation
Rubnich et al. An algorithm to identify robust convective weather avoidance polygons in en route airspace
Bartok et al. Data mining for fog prediction and low clouds detection
CN111967653B (en) Method for constructing airport runway wind forecast model and forecast method and system
Fonseca et al. Wind forecasts for rocket and balloon launches at the Esrange space center using the WRF model
Chan et al. LIDAR ground-based velocity track display analyses and surface observations of a vortex shedding event observed at the Hong Kong International Airport on April 11, 2011
Dean et al. Volcanic ash transport and dispersion models
CN203164444U (en) A device forecasting heavy rainfall by using the selection of boundary points and the determination of an aggregation area
Mandel An early look at the development of an unmanned automated surface aviation weather observation system
KR102746587B1 (en) Method of producing meteorological data for digital twin and system for producing meteorological data for digital twin
CN118981616A (en) Aviation flight weather analysis method and system
Xiao et al. Threshold analysis method for aircraft avoiding convective weather
CN207488522U (en) A kind of measure thunder and lightning group and the associated system of meteorologic factor

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant