Method for joint detection and identification of wind shear by using laser radar and wind profile radar
Technical Field
The invention relates to the field of meteorological detection and laser radar application, in particular to a method for identifying wind shear by combined detection of a laser radar and a wind profile radar.
Background
Wind shear is an atmospheric phenomenon, the wind direction, the wind speed vary in the horizontal and/or vertical distance in the air. Wind shear can be divided into horizontal shear of horizontal wind, vertical shear of horizontal wind, and vertical shear of vertical wind, the horizontal shear of horizontal wind is the change of wind direction and/or wind speed in horizontal distance, the vertical shear of horizontal wind is the change of wind direction and/or wind speed in vertical distance, the shear of vertical wind is the change of vertical wind (i.e. lift airflow) in horizontal or track direction, and the down-wash airflow is a form of shear of vertical wind, and is a strong down-wash airflow, and the down-wash airflow with small range and strong intensity is called micro-down-wash airflow. According to the statistics of the world meteorological organization and the international civil aviation organization, wind shear is the weather phenomenon which has the greatest threat to the flight safety of the airplane in the take-off and landing stage, and a plurality of major air accident accidents are caused by wind shear.
At present, the existing detection wind shear in China is single-station detection through a single radar, the whole airport runway is difficult to cover, the obtained data is limited, and the detection accuracy is low.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method for jointly detecting and identifying wind shear by using a laser radar and a wind profile radar, which can enlarge the detection coverage and improve the detection accuracy.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for detecting and identifying wind shear by combining a laser radar and a wind profile radar comprises the following steps:
s1, acquiring wind field data above the wind profile radar through the wind profile radar, and acquiring wind field data above the laser radar through the laser radar, wherein the wind field data comprise wind speed, height position of the wind speed and time of the wind speed;
s2, acquiring wind shear data over an airport runway through wind field data over a laser radar and wind field data over a wind profile radar, wherein the wind shear data comprise a wind shear value, wind shear generation time and a wind shear height position; the wind shear value is the absolute value of the difference between the wind speed acquired by the wind profile radar and the wind speed acquired by the laser radar at the same height;
and S3, if the wind shear value is larger than the standard value, the alarm system of the airport prompts alarm information, and the alarm information comprises the height position and the time when the wind shear value is larger than the standard value.
Further, the value of wind shear is Fi:
Fi=|Vi*cos(Di+θ)-vi*cos(di+θ)|
Wherein, FiIs the value of wind shear at the height i position, ViThe wind speed at the position of the wind profile radar altitude i is DiIs the wind direction v at the position i of the altitude above the wind profile radariIs the wind speed at the laser radar altitude i position, diThe wind direction is the wind direction at the position of the laser radar altitude i, namely the included angle between the wind direction and the horizontal plane, and theta is a constant.
Further, in the present invention,
wherein v is1i、v2i、v3i、v4iRespectively a first beam and a second beam of a lidarPartial velocity of wind, v, obtained at the height i position for each beam, third beam, fourth beam1i、v2i、v3i、v4iThe included angle with the horizontal direction is α, and the offset is the set calibration angle.
Further preferably, the first beam and the third beam are symmetrically arranged with respect to the vertical direction, and the second beam and the fourth beam are symmetrically arranged with respect to the vertical direction.
Further preferably, the first beam, the second beam, the third beam and the fourth beam are uniformly distributed on the circumference formed by the laser radar scanning.
Further preferably, at the same height i, the lidar acquires the wind speed v at the height i above the lidariTime of tLaserThe wind profile radar acquires the wind speed V at the altitude iiTime of tWind profileIf t isLaserAnd tWind profileIf the absolute value of the difference is less than the set threshold value, the wind speed v is savediOtherwise, the wind speed v is deletediAnd filling wind speed at adjacent time and/or adjacent position, wherein the wind speed after filling is used as the wind speed v 'at the laser radar overhead height i'i。
It is further preferred that the first and second liquid crystal compositions,
wherein, hi-1Is the height at the i-1 position, hiHeight at position i, hi+1Is the height at the i +1 position, hi+1Greater than hi,hiGreater than hi-1,xi+1Is hi+1Altitude wind velocity vi+1Abscissa, y, in a rectangular coordinate systemi+1Is hi+1Altitude wind velocity vi+1Ordinate, x, in a rectangular coordinate systemi-1Is hi-1Altitude wind velocity vi-1Abscissa, y, in a rectangular coordinate systemi-1Is hi-1Altitude wind velocity vi-1A vertical coordinate in a rectangular coordinate system; wherein, at a height hi+1At the wind of the first beam acquisitionHas a partial velocity v1(i+1)The partial velocity v of the wind picked up by the second beam2(i+1)The partial velocity v of the wind picked up by the third beam3(i+1)The partial velocity v of the wind picked up by the fourth beam4(i+1)Velocity v of windi+1Is v is1(i+1)、v2(i+1)、v3(i+1)And v4(i+1)The resultant speed of (c); at a height hi-1At the first beam, the partial velocity of the wind acquired by the first beam is v1(i-1)The partial velocity v of the wind picked up by the second beam2(i-1)The partial velocity v of the wind picked up by the third beam3(i-1)The partial velocity v of the wind picked up by the fourth beam4(i-1)Velocity v of windi-1Is v is1(i-1)、v2(i-1)、v3(i-1)And v4(i-1)The resultant velocity of (c).
Furthermore, the laser radar and the profile radar are respectively positioned at two ends of the length direction of the runway.
Further preferably, the adjacent time is 5s
The invention has the following beneficial effects:
(1) the method utilizes data provided by laser radar vertical profile mode detection and original wind profile radar data to perform data analysis and calculation, judge whether wind shear occurs and perform wind field data and alarm interface display.
The laser radar and the wind profile radar are simultaneously used for detecting, on one hand, the detection range can be enlarged, the air above the runway is covered, and more detection data are obtained, so that the accuracy of the obtained wind shear is improved. On the other hand, the laser radar and the profile radar are mutually used as references, whether the acquired data are accurate or not can be judged, the detection accuracy is further improved, and the occurrence of airplane accidents is reduced.
(2) The laser radar data and the wind profile radar data acquired in real time are identified in a data conversion calculation mode in a three-dimensional coordinate system, and the accuracy is high. The wind shear recognition of different height positions above the airport runway is carried out by using different detection radars, so that the limitation problem of a single type of radar is solved, and the problem that the traditional low-altitude wind shear early warning system based on the ground station cannot detect a high-altitude wind field is also solved.
(3) The detection method can acquire the wind shear of different heights above the runway so as to know the height of the fan to be safe.
(4) And when the time difference value of the wind speeds acquired by the two radars is greater than the threshold value, eliminating the wind speed, and filling the wind profile radar data and the laser radars at close time or adjacent positions, so that the accuracy of subsequent wind shear calculation is improved, invalid data are eliminated, the data quality is improved, and the accuracy of wind shear identification and the integrity of alarm information are improved.
(5) The laser radar has high resolution, improves the detection accuracy of horizontal wind vertical shear through data layered calculation, focuses on the overhead area of the landing point of the airplane, and ensures the flight safety of the airplane during landing.
Drawings
FIG. 1 is a schematic view of a radar beam scan of the present invention;
FIG. 2 is a flow chart of the present invention;
FIG. 3 is a schematic illustration of the distribution of the lidar and the wind profile radar of the present invention on an airstrip;
FIG. 4 is a schematic diagram of the detection result of the present invention on an alarm system.
Detailed Description
The technical scheme of the invention is clearly and completely described below by combining the embodiment and the attached drawings of the specification. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Examples
A method for detecting and identifying wind shear by combining a laser radar and a wind profile radar is disclosed, as shown in FIG. 2, and comprises the following steps:
s1, as shown in FIG. 3, the wind profile radar and the laser radar are distributed at two ends of the runway, the runway is arranged in the north and south, and the airplane flies to the laser radar from the wind profile radar and also can fly to the wind profile radar from the laser radar. Two identical radars, such as two lidar, may be employed. In this embodiment, two different types of radars, namely, a wind profile radar and a laser radar, are used, which solves the problems of insufficient detection of other devices and simultaneous detection of horizontal shear and horizontal wind vertical shear. And the detection of the gravity point region in accordance with the ICAO international standard. Wind field data above the wind profile radar is obtained through the wind profile radar, wind field data above the laser radar is obtained through the laser radar, the wind field data comprise wind speed, the height position where the wind speed is located and the time when the wind speed occurs, and the obtained wind speed is shown in figure 4. The method comprises the following specific steps:
the lidar is deployed at the north end of the runway, and periodically performs scanning of a five-beam system by adopting a vertical profile scanning mode, as shown in fig. 1, wherein four beams form an included angle with the horizontal plane. The wind profile radar is deployed at the south end of the runway, and periodically scans to acquire detection data in the vertical direction. For example: acquiring a group of data of a laser radar vertical profile mode every 85 seconds, wherein the group of data comprises a plurality of radial data in the vertical direction, and one radial data comprises information such as wind direction and wind speed on layers with different heights; and acquiring a group of wind profile radar profile data every 150 seconds or so, wherein the data comprises information such as wind direction and wind speed on layers with different heights. Radial, i.e. along the direction of the beam.
Wherein v is1i、v2i、v3i、v4iThe method comprises the steps of obtaining the wind partial speed of a first wave beam, a second wave beam, a third wave beam and a fourth wave beam of the laser radar at a height i position respectively, establishing the wind partial speed on a polar coordinate, wherein the wind speed is a vector, two wind speeds are subtracted by 2, the absolute value of the sum of the two wind speeds is obtained, the first wave beam and the third wave beam are symmetrically arranged in the vertical direction, and the second wave beam and the fourth wave beam are symmetrically arranged in the vertical directionThe first beam, the second beam, the third beam and the fourth beam are symmetrically arranged in the vertical direction and are uniformly distributed on the circumference formed by scanning the laser radar, and v is1i、v2i、v3i、v4iThe angle to the horizontal is α, the offset is the set calibration angle, and the offset is 30 °.
At the position of the same height i, the laser radar acquires the wind speed v at the position of the height i above the laser radariTime of tLaserThe wind profile radar acquires the wind speed V at the altitude iiTime of tWind profileIf t isLaserAnd tWind profileIf the absolute value of the difference is less than the set threshold value, the wind speed v is savediOtherwise, the wind speed v is deletediAnd filling wind speed at adjacent time and/or adjacent position, wherein the wind speed after filling is used as the wind speed v 'at the laser radar overhead height i'i. In this embodiment, the threshold is set to 60s, and the wind speed V is obtained at the height i at the adjacent timeiThe front 5s and the back 5 s.
Wherein, hi-1Is the height at the i-1 position, hiHeight at position i, hi+1Is the height at the i +1 position, hi+1Greater than hi,hiGreater than hi-1,xi+1Is hi+1Altitude wind velocity vi+1Abscissa, y, in a rectangular coordinate systemi+1Is hi+1Altitude wind velocity vi+1Ordinate, x, in a rectangular coordinate systemi-1Is hi-1Altitude wind velocity vi-1Abscissa, y, in a rectangular coordinate systemi-1Is hi-1Altitude wind velocity vi-1A vertical coordinate in a rectangular coordinate system; wherein, at a height hi+1At the first beam, the partial velocity of the wind acquired by the first beam is v1(i+1)The partial velocity v of the wind picked up by the second beam2(i+1)The partial velocity v of the wind picked up by the third beam3(i+1)The partial velocity v of the wind picked up by the fourth beam4(i+1)Velocity v of windi+1Is v is1(i+1)、v2(i+1)、v3(i+1)And v4(i+1)The resultant speed of (c); at a height hi-1At the first beam, the partial velocity of the wind acquired by the first beam is v1(i-1)The partial velocity v of the wind picked up by the second beam2(i-1)The partial velocity v of the wind picked up by the third beam3(i-1)The partial velocity v of the wind picked up by the fourth beam4(i-1)Velocity v of windi-1Is v is1(i-1)、v2(i-1)、v3(i-1)And v4(i-1)Resultant velocity, wind speed vi+1Partial velocity and wind velocity vi-1Are different, so that the wind speed vi+1And wind speed vi-1And they are not identical in their coordinates on the coordinate system.
S2, acquiring wind shear data over an airport runway through wind field data over a laser radar and wind field data over a wind profile radar, wherein the wind shear data comprise a wind shear value, wind shear generation time and a wind shear height position; the value of the wind shear is the absolute value of the difference between the wind speed acquired by the wind profile radar and the wind speed acquired by the laser radar on the same height.
Value of wind shear Fi:
Fi=|Vi*cos(Di+θ)-vi*cos(di+θ)|
Wherein, FiIs the value of wind shear at the height i position, ViThe wind speed at the position of the wind profile radar altitude i is DiIs the wind direction v at the position i of the altitude above the wind profile radariIs the wind speed at the laser radar altitude i position, diThe wind direction is the wind direction at the position of the laser radar altitude i, namely the included angle between the wind direction and the horizontal plane, theta is a constant, namely the north-offset angle of the north side of the runway, and in the embodiment, theta is 30 degrees.
If in step S1, if tLaserAnd tWind profileThe absolute value of the difference is greater than a set threshold value, the value of wind shear FiV in the calculation formulaiIs replaced by v'i。
S3, if the value of wind shear is larger than the standard value, the alarm system of the airport prompts the alarm information, otherwise, the alarm information is not prompted, as shown in fig. 4, the DBS mode alarm displays the alarm information, and the alarm information includes the height position and time at which the value of wind shear is larger than the standard value. In this example, the standard value was 7.5 m/s.