High-precision wind shear recognition algorithm and device based on wind speed filtering and storage medium

Abstract

The invention discloses a high-precision wind shear recognition algorithm based on wind speed filtering. The algorithm comprises the following steps: S1, collecting wind field data by a wind measurement radar; the wind field data comprises radial wind speeds and radar pitch angles # imgabs0 # S2, and filtering the radial wind speeds of different distance doors; s3, performing inversion based on the filtered radial wind speed and the radar pitch angle to obtain a horizontal wind speed and a horizontal wind direction; s4, calculating wind shear factors of different distance doors according to the filtered horizontal wind speed and horizontal wind direction; and S5, combining the wind shear factor with a wind shear standard, and carrying out wind shear identification judgment, compared with the prior art, the radial wind speed of each distance door is filtered, and the calculation precision of the wind shear factor is greatly improved.

Description

Technical Field

[0001] The present invention relates to the field of wind shear warning, and in particular to a high-precision wind shear identification algorithm, device and storage medium based on wind speed filtering. Background Art

[0002] Wind shear is an atmospheric phenomenon, which is a sudden change in wind speed in the horizontal and vertical directions. Wind shear is an important factor leading to flight accidents, especially low-altitude wind shear.

[0003] Wind shear can be divided into horizontal wind shear, vertical wind shear and vertical wind shear according to wind direction. Some vertical shear is mainly judged based on wind speed and direction. In the process of measuring vertical shear with three-dimensional laser wind radar, a large amount of radial wind speed needs to be collected. These radial wind speeds will contain a large amount of invalid data or low-precision data. The existing data filtering method is relatively simple, and the low-precision data filtering is incomplete, which greatly affects the calculation of wind shear factor, resulting in inaccurate calculation results and prone to false alarms and delayed reports. Summary of the invention

[0004] The purpose of the present invention is to address the problem that the wind speed collected in the prior art is incompletely filtered, resulting in low calculation accuracy of the wind shear factor. To address this problem, a high-precision wind shear identification algorithm, device and storage medium based on wind speed filtering are proposed.

[0005] In order to achieve the above object, the present invention adopts the following technical solution:

[0006] A high-precision wind shear identification algorithm based on wind speed filtering includes the following steps:

[0007] S1. Basic wind radar collects wind field data; wind field data includes radial wind speed and radar pitch angle

[0008] S2, filtering radial wind speeds at different distances from doors;

[0009] S3, obtaining the horizontal wind speed and horizontal wind direction based on the filtered radial wind speed and radar pitch angle inversion;

[0010] S4. Calculate the wind shear factors of different distance gates according to the filtered horizontal wind speed and horizontal wind direction;

[0011] S5. The wind shear factor is combined with the wind shear standard to identify and judge wind shear.

[0012] As a further preferred embodiment of the present invention, the wind field data includes radial wind speeds of five beams in one period, which are UrV, UrN, UrW, UrE, and UrS respectively.

[0013] As a further preferred embodiment of the present invention, filtering radial wind speeds at doors at different distances includes:

[0014] S21. Assuming the radial wind speed of UrV is Vfifth, UrS and UrN are decomposed into:

[0015]

[0016] Among them, V x1 is the horizontal component of UrS; V x2 is the horizontal component of UrN;

[0017] S22. Calculate the north-south component error △V x ;

[0018]

[0019] S23, UrE and UrW are decomposed into:

[0020]

[0021] Among them, V y1 is the horizontal component of UrW; V y2 is the horizontal component of UrE;

[0022] S24. Calculate the east-west component error △V y ;

[0023] S25, based on the requirements of wind speed accuracy less than 0.1m / s and wind direction accuracy less than 1°, we get

[0024]

[0025] Where threshold = △V x =△V y ;

[0026] S26, calculate the east-west component error △V y and the north-south component error △V x The error range of

[0027]

[0028] S27. Filter the east-west component and the north-south component according to the error range.

[0029] As a further preferred embodiment of the present invention, the inversion based on the radial wind speed and the radar pitch angle to obtain the horizontal wind speed and the horizontal wind direction includes:

[0030] S31, based on UrS and UrN decomposition to obtain the north-south component V x ;

[0031] S32, based on UrE and UrW decomposition to obtain the east-west component V y ;

[0032] S33. Get the zenith velocity V based on UrS, UrN, UrE and UrW z ;

[0033] S34. According to the east-west component, the north-south component and the zenith component, the horizontal wind speed and the horizontal wind direction are inverted.

[0034] A high-precision wind shear identification device based on wind speed filtering, the device comprising one or more processors and a memory, the memory being used to store one or more programs; when the one or more programs are executed by the processor, the processor is enabled to implement a wind shear identification method based on laser radar.

[0035] A computer storage medium stores at least one program instruction, and the at least one program instruction is used to be loaded and executed by a processor to implement a wind shear identification method based on laser radar.

[0036] The high-precision wind shear identification algorithm, device and storage medium based on wind speed filtering proposed in the present invention have the following beneficial effects compared with the prior art:

[0037] 1. The present invention greatly improves the calculation accuracy of the wind shear factor by filtering the radial wind speed of each range gate;

[0038] 2. The radial wind speed of one cycle and five beams obtained by radar is inverted, and large-scale wind speed inversion at different range gates can be performed;

[0039] 3. The wind shear factor can be used to obtain the vertical wind shear conditions at different distance gates, providing accurate alarms for aircraft takeoff and landing. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 This is a schematic diagram of the radar of the present invention obtaining radial wind speed with one cycle and five beams.

[0041] Figure 2 This is the front view of the horizontal wind speed inversion by the DBS mode of the present invention;

[0042] Figure 3 is the horizontal wind speed V of the synthetic true vector wind h , schematic diagram of horizontal wind direction β0;

[0043] Figure 4 This is the wind shear inversion diagram obtained by the present invention through DBS. DETAILED DESCRIPTION

[0044] The present invention is described in detail below with reference to the accompanying drawings and specific embodiments.

[0045] Combination Figure 1-4 , a high-precision wind shear identification algorithm based on wind speed filtering, comprising the following steps:

[0046] S1. The basic wind radar collects wind field data; the wind field data includes radial wind speed with 5 beams in one cycle, namely UrV, UrN, UrW, UrE, and UrS; the wind field data also includes the radar pitch angle like Figure 1 As shown,.

[0047] S2. Filter the radial wind speed at different distances from the door.

[0048] The specific steps include:

[0049] S21. Assuming the radial wind speed of UrV is Vfifth, UrS and UrN are decomposed into:

[0050]

[0051] Among them, V x1 is the horizontal component of UrS; V x2 is the horizontal component of UrN;

[0052] S22. Calculate the north-south component error △V x ;

[0053]

[0054] S23, UrE and UrW are decomposed into:

[0055]

[0056] Among them, V y1 is the horizontal component of UrW; V y2 is the horizontal component of UrE.

[0057] S24. Calculate the east-west component error △V y .

[0058] S25, based on the requirements of wind speed accuracy less than 0.1m / s and wind direction accuracy less than 1°, we get

[0059]

[0060] Where threshold = △V x =△V y .

[0061] S26, calculate the east-west component error △V y and the north-south component error △V x The error range of

[0062]

[0063] S27, according to the error range of the east-west component △V y and the north-south component △V x to filter.

[0064] S3. Horizontal wind speed and horizontal wind direction are obtained based on the filtered radial wind speed and radar elevation angle.

[0065] The following steps are involved:

[0066] The radial wind speeds UrN and UrS collected in the north-south direction are back-projected to the north-south direction (x direction) to obtain the north-south component V of the true wind vector x , the radial wind speeds UrN and UrS collected in the north-south direction are back-projected to the zenith direction (z direction) to obtain the zenith component Vz of the true wind vector.

[0067] S31, based on UrS and UrN decomposition to obtain the north-south component V x ;

[0068]

[0069] get

[0070] S32, based on UrE and UrW decomposition to obtain the east-west component V y ;

[0071]

[0072] get

[0073] S33. Get the zenith velocity V based on UrS, UrN, UrE and UrW z ;

[0074] S34, according to the east-west component V x , north-south component V y and the zenithal component V z , the horizontal wind speed V is obtained by inversion h and horizontal wind direction β0,

[0075] get

[0076]

[0077] S4. Calculate the wind shear factor P of different distance gates according to the filtered horizontal wind speed and horizontal wind direction.

[0078]

[0079] Where i is the range gate number; V h is the horizontal wind speed of the current range gate, β0 is the horizontal wind direction of the current range gate; when i=1, V h1 is the horizontal wind speed of the range gate one layer above the current range gate, β1 is the horizontal wind direction of the range gate one layer above the current range gate; when i = 2, V h2 is the horizontal wind speed of the next range gate, β2 is the horizontal wind direction of the next range gate; the wind shear factor P is calculated. i , P1 is the wind shear factor formed between the current range gate and the previous range gate, and P2 is the wind shear factor formed between the current range gate and the next range gate.

[0080] V h1 、V h2 , β1 and β2 are calculated in the same way as V h The calculation method is the same as β.

[0081] Wind shear factor P = max(P1, P2);

[0082] Compare the wind shear factor P1 formed between the current range gate and the previous range gate with the wind shear factor P2 formed between the current range gate and the next range gate, and select the largest one as the wind shear factor P of the current range gate.

[0083] S5. The wind shear factor is combined with the wind shear standard to identify and judge wind shear.

[0084] Based on the ICAO wind shear standard, the current wind shear warning range is determined according to the wind shear factor P of the current range gate.

[0085] Table 1 shows the conversion of ICAO wind shear standards

[0086]

[0087] Issue a wind shear warning based on the current wind shear warning range.

[0088] When the wind shear factor P is within the range of 0-2.56, it means that the current vertical wind shear is slight and no warning is issued;

[0089] When the wind shear factor P is within the range of 2.57-4.11, it indicates that the current vertical wind shear is moderate, and a moderate wind shear warning is issued;

[0090] When the wind shear factor P is within 4.12-6.17, it means that there is strong vertical wind shear and a strong wind shear warning is issued;

[0091] When the wind shear factor P is above 6.18, it means that there is severe vertical wind shear and a strong wind shear warning is issued.

[0092] A high-precision wind shear identification device based on wind speed filtering, the device comprising one or more processors and a memory, the memory being used to store one or more programs; when the one or more programs are executed by the processor, the processor is enabled to implement a wind shear identification method based on laser radar.

[0093] A computer storage medium stores at least one program instruction, and the at least one program instruction is used to be loaded and executed by a processor to implement a wind shear identification method based on laser radar.

[0094] The above shows and describes the basic principles, main features and advantages of the present invention. Those skilled in the art should understand that the above embodiments do not limit the present invention in any form, and any technical solution obtained by equivalent replacement or equivalent transformation falls within the protection scope of the present invention.

Claims

1. A high-precision wind shear identification algorithm based on wind speed filtering, characterized in that: The following steps are involved: S1. Basic wind radar collects wind field data; wind field data includes radial wind speed and radar pitch angle S2, filtering radial wind speeds at different distances from doors; S3, obtaining the horizontal wind speed and horizontal wind direction based on the filtered radial wind speed and radar pitch angle inversion; S4. Calculate the wind shear factors of different distance gates according to the filtered horizontal wind speed and horizontal wind direction; S5. The wind shear factor is combined with the wind shear standard to identify and judge wind shear.

2. The high-precision wind shear identification algorithm based on wind speed filtering according to claim 1 is characterized in that: The wind field data includes radial wind speeds of five beams in one period, namely UrV, UrN, UrW, UrE, and UrS.

3. The high-precision wind shear identification algorithm based on wind speed filtering according to claim 1 is characterized in that: The filtering of radial wind speeds at doors at different distances includes: S21. Assuming the radial wind speed of UrV is Vfifth, UrS and UrN are decomposed into: Among them, V x1 is the horizontal component of UrS; V x2 is the horizontal component of UrN; S22. Calculate the north-south component error △V x ; S23, UrE and UrW are decomposed into: Among them, V y1 is the horizontal component of UrW; V y2 is the horizontal component of UrE; S24. Calculate the east-west component error △V y ; S25, based on the requirements of wind speed accuracy less than 0.1m / s and wind direction accuracy less than 1°, we get Where threshold = △V x =△V y ; S26, calculate the east-west component error △V y and the north-south component error △V x The error range of S27. Filter the east-west component and the north-south component according to the error range.

4. The high-precision wind shear identification algorithm based on wind speed filtering according to claim 3 is characterized in that: The inversion based on radial wind speed and radar pitch angle to obtain horizontal wind speed and horizontal wind direction includes: S31, based on UrS and UrN decomposition to obtain the north-south component V x ; S32, based on UrE and UrW decomposition to obtain the east-west component V y ; S33. Based on UrS, UrN, UrE and UrW, the zenith velocity V is obtained z ; S34. According to the east-west component, the north-south component and the zenith component, the horizontal wind speed and the horizontal wind direction are inverted.

5. A high-precision wind shear identification device based on wind speed filtering, characterized in that: The device includes one or more processors and a memory, wherein the memory is used to store one or more programs; when the one or more programs are executed by the processor, the processor is enabled to implement a wind shear identification method based on laser radar.

6. A computer storage medium, characterized in that: The computer storage medium stores at least one program instruction, and the at least one program instruction is used to be loaded and executed by a processor to implement a wind shear identification method based on laser radar.

Images