Disclosure of Invention
In view of this, the embodiments of the present disclosure provide a reactive wind shear false alarm suppression method that fuses weather radar information, so as to achieve the purpose of reducing the reactive wind shear false alarm rate.
The embodiment of the specification provides the following technical scheme that a reactive wind shear false alarm suppression method fusing weather radar information comprises the following steps:
step one, wind field inversion is carried out on wind field information;
Step two, introducing wind field information to judge turbulence, correcting integration time length and extracting a predicted wind shear F factor;
Step three, calculating a reactive wind shear F factor according to a reactive wind shear mode, and comparing the reactive wind shear F factor with an alarm threshold to obtain a wind shear alarm result;
Step four, calculating and comparing the reactive wind shear F factor with the F factor of the aircraft position point in wind field inversion, if the integral average difference value of the reactive wind shear F factor and the predictive wind shear F factor is larger than the tolerance, detecting the maneuver type of the reactive wind shear module and finishing the correction of the reactive wind shear F factor;
And fifthly, based on the time interval sensitivity characteristic, obtaining a corrected F factor by weighting and fusing the predicted wind shear F factor and the reactive wind shear F factor, and re-executing the third step and giving a wind shear alarm result.
Further, the first step includes:
Step 1.1, wind field information obtained from predictive wind shear is defined under a radar coordinate system taking an airplane position as an origin, wherein each node comprises horizontal wind speed, vertical wind speed, F factor and turbulence characteristic value information;
Step 1.2, determining a longitude and latitude range of predictive wind shear forward looking detection according to the position, heading information and detection distance of the aircraft, taking set three points in the predictive wind shear forward looking detection range, converting the set three points into longitude and latitude values by using an ink card support projection mode, and taking out the maximum and minimum longitude and latitude values as boundary values of the predictive wind shear forward looking detection range;
And 1.3, meshing the area in the spatial range of the predictive wind shear forward-looking detection, forming a longitude and latitude grid, and filling the predictive wind shear detection result into the longitude and latitude grid.
Further, the second step is to find corresponding position information in the longitude and latitude grid formed in the first step according to the position of the aircraft, check whether turbulence information exists at the position of the aircraft in the longitude and latitude grid, adjust the integral time of the reactive wind shear factor F according to the turbulence characteristic value, and extract the predictive wind shear factor F.
Further, the third step specifically includes:
step 3.1, the reactive wind shear F factor is according to the formula Performing a calculation in whichFor the horizontal component of wind along the flight path, W is the vertical component of wind, V a is airspeed, g is gravitational acceleration;
And 3.2, comparing the integral average value of the reactive wind shear F factor for a period of time with an alarm threshold and giving an alarm result.
Further, step four includes step 4.1, calculating a difference integral average of the reactive wind shear factor F and the predictive wind shear factor F, and comparing the difference integral average with a decision threshold to determine whether a maneuver effect exists.
Further, step four comprises step 4.2, if there is a maneuver effect, comprising the steps of:
Aiming at the change characteristics of flight parameters in a wind shear environment, qualitatively analyzing the change trend of the flight parameters under different maneuvering types;
according to the change trend of flight parameters under different maneuver types, a matrix A1= [ a 1 a2 a3 a4 a5 a6 ] of the occurrence probability of each maneuver is established, wherein a 1~a6 represents the occurrence probability of the maneuver in the current flight stage of the aircraft;
A matrix b1= [ B 1 b2 b3 b4b5 b6 ] is established based on the similarity of maneuver type to the flight parameters in the wind shear environment, where B 1~b6 represents the probability of occurrence of the maneuver identified based on the aircraft parameters.
Judging the occurrence factor lambda=f (A 1,B1) of each current maneuvering type of the airplane according to A1= [ a 1 a2 a3 a4 a5 a6 ] and B1= [ B 1 b2 b3 b4 b5 b6 ];
Obtaining an experience influence factor f of the current maneuver according to the search maneuver type library table;
obtaining a large maneuver influence factor according to the experience influence factor, the occurrence factor and the maneuver influence correction function Fun=fun (f, lambda);
the corrected reactive wind shear factor F 1 is obtained by F 1=Fun·F1.
Further, the fourth step further comprises the step of 4.3, if no maneuver effect exists, generating a wind shear alarm.
Further, the fifth step includes:
The method comprises the steps of carrying out weighted fusion on a predicted wind shear F factor and a reactive wind shear F factor through F mix(t)=kF1(t)+(1-k)F2 (t), wherein F 2 (t) is the predicted wind shear F factor, k= (t 1-t)/(t1-t0), t is the moment when the predicted wind shear obtains wind field information of the position, t 0 is the fusion calculation starting time, and t 1 is the starting time of a detection blind zone;
By passing through Performing integral calculation and correcting the F factor, wherein t 2 is the current moment;
And comparing and judging the corrected F factor with an alarm threshold and generating an alarm result.
Compared with the prior art, the wind field inversion method has the advantages that at least one technical scheme adopted by the embodiment of the specification at least comprises the steps of introducing a wind field detected by a weather radar into a near-ground warning system for wind field inversion, and arranging wind field data into a data form which is convenient to use for reactive wind shear in order to reduce turbulence in the reactive wind shear and wind shear false alarms caused by motor actions and improve the performance of the reactive wind shear alarms. And correcting the integration time of the current F factor calculation by using turbulence characteristics of the detected wind field information in the stage of reactive wind shear alarm processing, and reducing the F factor contribution quantity caused by disturbance as much as possible. In addition, the reactive wind shear alarm processing combines with the large maneuver influencing factors, and calculates different weight proportion prediction wind shear F factors and reactive wind shear F factors based on time interval sensitivity, so that the false alarm probability is effectively reduced.
Detailed Description
Embodiments of the present application will be described in detail below with reference to the accompanying drawings.
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
As shown in fig. 1 to 5, the embodiment of the invention provides a reactive wind shear false alarm suppression method for fusing weather radar information, which comprises the following steps:
step one, wind field inversion is carried out on wind field information;
Step two, introducing wind field information to judge turbulence, correcting integration time length and extracting a predicted wind shear F factor;
Step three, calculating a reactive wind shear F factor according to a reactive wind shear mode, and comparing the reactive wind shear F factor with an alarm threshold to obtain a wind shear alarm result;
Step four, calculating and comparing the reactive wind shear F factor with the F factor of the aircraft position point in wind field inversion, if the integral average difference value of the reactive wind shear F factor and the predictive wind shear F factor is larger than the tolerance, detecting the maneuver type of the reactive wind shear module and finishing the correction of the reactive wind shear F factor;
And fifthly, based on the time interval sensitivity characteristic, obtaining a corrected F factor by weighting and fusing the predicted wind shear F factor and the reactive wind shear F factor, and re-executing the third step and giving a wind shear alarm result.
In order to reduce turbulence in reactive wind shear and wind shear false alarms caused by maneuvering, and improve the performance of reactive wind shear alarms, a wind field detected by a weather radar is introduced into a near-ground alarm system, wind field inversion is performed, and wind field data are organized into a data form which is convenient to use for reactive wind shear. And correcting the integration time of the current F factor calculation by using turbulence characteristics of the detected wind field information in the stage of reactive wind shear alarm processing, and reducing the F factor contribution quantity caused by disturbance as much as possible. In addition, the reactive wind shear alarm processing combines with the large maneuver influencing factors, and calculates different weight proportion prediction wind shear F factors and reactive wind shear F factors based on time interval sensitivity, so that the false alarm probability is effectively reduced.
The first step comprises the following steps:
Step 1.1, wind field information obtained from predictive wind shear is defined under a radar coordinate system taking an airplane position as an origin, wherein each node comprises horizontal wind speed, vertical wind speed, F factor and turbulence characteristic value information;
Step 1.2, determining a longitude and latitude range of predictive wind shear forward looking detection according to the position, heading information and detection distance of the aircraft, taking set three points in the predictive wind shear forward looking detection range, converting the set three points into longitude and latitude values by using an ink card support projection mode, and taking out the maximum and minimum longitude and latitude values as boundary values of the predictive wind shear forward looking detection range;
And 1.3, meshing the area in the spatial range of the predictive wind shear forward-looking detection, forming a longitude and latitude grid, and filling the predictive wind shear detection result into the longitude and latitude grid.
And step two, searching corresponding position information in the longitude and latitude grid formed in the step one according to the position of the aircraft, checking whether turbulence information exists at the position of the aircraft in the longitude and latitude grid, adjusting the integration time of the reactive wind shear factor F according to the turbulence characteristic value, and extracting the predictive wind shear factor F.
The third step specifically comprises:
step 3.1, the reactive wind shear F factor is according to the formula Performing a calculation in whichFor the horizontal component of wind along the flight path, W is the vertical component of wind, V a is airspeed, g is gravitational acceleration;
And 3.2, comparing the integral average value of the reactive wind shear F factor for a period of time with an alarm threshold and giving an alarm result.
Step four, calculating a difference integral average of the reactive wind shear factor F and the predictive wind shear factor F, and comparing the difference integral average with a judgment threshold to determine whether the maneuvering effect exists.
Step 4.2, if there is a maneuver effect, comprising the steps of:
Aiming at the change characteristics of flight parameters in a wind shear environment, qualitatively analyzing the change trend of the flight parameters under different maneuvering types;
according to the change trend of flight parameters under different maneuver types, a matrix A1= [ a 1 a2 a3 a4 a5 a6 ] of the occurrence probability of each maneuver is established, wherein a 1~a6 represents the occurrence probability of the maneuver in the current flight stage of the aircraft;
A matrix b1= [ B 1 b2 b3 b4b5 b6 ] is established based on the similarity of maneuver type to the flight parameters in the wind shear environment, where B 1~b6 represents the probability of occurrence of the maneuver identified based on the aircraft parameters.
Judging the occurrence factor lambda=f (A 1,B1) of each current maneuvering type of the airplane according to A1= [ a 1 a2 a3 a4 a5 a6 ] and B1= [ B 1 b2 b3 b4 b5 b6 ];
Obtaining an experience influence factor f of the current maneuver according to the search maneuver type library table;
obtaining a large maneuver influence factor according to the experience influence factor, the occurrence factor and the maneuver influence correction function Fun=fun (f, lambda);
the corrected reactive wind shear factor F 1 is obtained by F 1=Fun·F1.
And 4.3, if the motor action influence does not exist, generating a wind shear alarm.
The fifth step comprises:
The method comprises the steps of carrying out weighted fusion on a predicted wind shear F factor and a reactive wind shear F factor through F mix(t)=kF1(t)+(1-k)F2 (t), wherein F 2 (t) is the predicted wind shear F factor, k= (t 1-t)/(t1-t0), t is the moment when the predicted wind shear obtains wind field information of the position, t 0 is the fusion calculation starting time, and t 1 is the starting time of a detection blind zone:
By passing through Performing integral calculation and correcting the F factor, wherein t 2 is the current moment;
And comparing and judging the corrected F factor with an alarm threshold and generating an alarm result.
Compared with the reactive wind shear warning method of the existing near-earth warning system, the wind shear warning method disclosed by the invention has the advantages that the wind field detection information of the predictive wind shear is fused, the measured turbulence characteristic information is utilized to adjust the judgment sensitivity, the F factor fluctuation component contributed by turbulence is reduced, and the false alarm is reduced.
According to the invention, a fusion mode of the predictive wind shear F factor and the reactive wind shear F factor is introduced, and when the difference value of the reactive wind shear F factor and the predictive wind shear F factor is larger than a threshold value, the reactive wind shear module carries out weighted fusion on the reactive F factor and the predictive F factor according to the maneuvering state of the aircraft, so that reactive wind shear false alarm caused by maneuvering is reduced.
Specific application examples of the embodiment of the invention are as follows:
step one, inversion of a predictive wind shear wind field.
As shown in fig. 5, wind field information obtained from predicted wind shear is defined under a radar coordinate system with an aircraft position as an origin, wherein each node contains horizontal wind speed, vertical wind speed, F-factor, and turbulence characteristic value information. And determining longitude and latitude ranges Lat_Min, lat_Max, lon_Min and Lon-Max of the predictive wind shear forward looking detection according to the aircraft position, the heading information and the detection distance. Taking three points A-for_Dis, 0, B0, for_Dis, C for_Dis,0, converting them into longitude and latitude values by using the ink card support projection mode, and taking out the maximum and minimum longitude and latitude values as Lat_Min, lat_Max, lon_Min and Lon-Max boundary values. And meshing the area in the space range of [ Lat_Min, lat_Max ], [ Lon_Min, lon_Max ], and filling the predicted wind shear detection result into the longitude and latitude grid.
And secondly, introducing wind field information to judge turbulence, correcting the integration time length and extracting a predicted wind shear F factor.
And (3) searching corresponding position information in the inversion wind field result matrix formed in the step (I) according to the position of the aircraft, checking whether turbulence information exists at the current position, and adjusting the integration time of the F factor according to the size of the turbulence characteristic value. The predicted wind shear F factor is extracted for subsequent calculation.
Step three, reactive wind shear standard F factor alarm judgment
The reactive wind shear F factor is calculated according to equation 1, whereFor the horizontal component of wind along the flight path, W is the vertical component of wind, V a is airspeed and g is gravitational acceleration.
The F-factor represents a wind field term of the aircraft's response to wind shear, defining a threshold for dangerous wind shear in the on-board wind shear warning device. And judging whether the alarm threshold value is exceeded or not by adopting the F factor integral average value for a period of time, and giving an alarm result.
And step four, judging whether the maneuvering influence exists by comparing the F factor of the reactive wind shear with the integral average of F factor differences of airplane position points in the inversion wind field, if so, correcting the maneuvering influence of the F factor of the reactive wind shear, otherwise, generating a reactive wind shear alarm.
Aiming at the change characteristics of flight parameters in wind shear environment, the common maneuver types are listed, the change trend of each main flight parameter and the flight characteristics of civil aircraft are qualitatively analyzed, and the comparison table of the change trend of the flight parameters under different maneuver types is listed in Table 1.
TABLE 1 flight parameter variation trend reference table under different maneuver types
The calculation flow of the large maneuver influence factor is shown in fig. 3, and the specific steps are as follows:
Firstly, determining the occurrence probability of each maneuver type under the current flight scene through an analytic hierarchy process (AHP method) according to the input flight parameters and maneuver type database, and establishing a matrix A1= [ a 1 a2 a3 a4 a5 a6 ] of the occurrence probability of each maneuver by combining the characteristics of a low-altitude wind shear wind field and scene recognition, and similarly, establishing a matrix B1= [ B 1 b2 b3 b4 b5 b6 ] according to the similarity of the maneuver type and the flight parameters under the wind shear environment. Next, according to probability and similarity matrix information of each maneuver, determining occurrence factor λ=f (a 1,B1) of each maneuver at present, and looking up maneuver influence factors through a maneuver type library to determine experience influence factor f. Then, according to the experience influence factor and the occurrence factor, a large maneuver influence factor is obtained by the maneuver influence correction function fun=fun (f, λ). Finally, the reactive wind shear factor F 1, i.e., F 1=Fun·F1, is corrected.
And fifthly, multi-source F factor fusion is carried out, and then calculation and alarm detection are carried out.
The calculation schematic diagram of the factor F is shown in fig. 4, the total fusion time is T, and T 2 is the current time. Due to the presence of radar detection dead zone, i.e. the time interval t 1,t2, predictive wind shear has no fusion detection data available, this interval f=f 1.
The multisource F factor fusion interval is [ t 0,t1 ]. Because wind field changes faster, the predictive wind shear F 2 factor has a time-dependent characteristic, i.e., the longer the time, the lower the confidence, so a time variable is introduced as a criterion for measuring the weight. Based on the time interval sensitivity characteristic, the multi-source F factors (reaction formula F 1 and predictive formula F 2) are weighted and fused according to formula (2), wherein k= (t 1-t)/(t1-t0), and t is the time when the predictive wind shear obtains the position wind field information.
Fmix(t)=kF1(t)+(1-k)F2(t) (2)
The delta is obtained by integrating and averaging the recalculated F factor, and the calculation formula is as follows:
And comparing the wind shear with the alarm threshold value, generating a reactive wind shear alarm when the calculation result is higher than the alarm threshold value, and otherwise, not alarming.
The foregoing description of the embodiments of the invention is not intended to limit the scope of the invention, so that the substitution of equivalent elements or equivalent variations and modifications within the scope of the invention shall fall within the scope of the patent. In addition, the technical characteristics and technical scheme, technical characteristics and technical scheme can be freely combined for use.