Method, system, device and storage medium for generating streamline of wind field and visual wind shear product in airport area

Abstract

The present invention discloses a method, system, device and storage medium for generating airport area wind field streamline and visual wind shear products, belonging to the technical field of airport meteorological detection and analysis. The method includes: obtaining triangular partitions; obtaining the wind speed vector of any point within the triangular partitions; setting the starting point of the wind field streamline of the wind and obtaining the wind speed vector of the starting point; starting from each starting point of the triangular partitions, obtaining the atmospheric particle trajectories on the wind field streamline and drawing the wind field streamline; obtaining and marking the maximum wind shear direction and the maximum wind shear value; obtaining divergence and vorticity, calculating and drawing the normalized product of the displayed divergence and vorticity. The wind field streamline products constructed by the present invention, the calculated and drawn map of the maximum wind shear magnitude and direction in the partition, and the calculated and drawn map of the divergence vorticity product greatly improve the overall trend display degree, accuracy and visualization degree of the airport wind field information for aircraft takeoff and landing, and can provide better reference for the safety guarantee work of airport aircraft takeoff and landing.

Description

Technical Field

[0001] The present invention belongs to the technical field of airport meteorological detection and analysis, and in particular relates to a method, system, equipment and storage medium for generating airport regional wind field streamlines and visualized wind shear products. Background Art

[0002] Among the meteorological factors related to airport aviation safety, wind is most closely related to flight and is directly related to the safety of aircraft takeoff and landing. In particular, low-altitude wind shear and turbulence seriously endanger the safety of aviation activities, which can cause aircraft to land prematurely, run off the runway, and hit obstacles. Therefore, wind field information is the most concerned information in airport aviation meteorology. Wind field detection, wind field parameter solution, and wind shear identification have always been hot spots and difficulties that have been studied at home and abroad.

[0003] In the related technologies, the airport wind field information detection equipment mainly includes remote sensing equipment based on radar and field measurement equipment based on anemometers (wind poles). The advantage of remote sensing equipment is that it has a long detection distance, but it has the problem of limited resolution and being greatly affected by weather. The detection assumption of field measurement equipment such as wind poles is limited, but it can accurately measure the wind speed vector at the location, which is more suitable for the detection of low-altitude wind information in a small range such as airport runways. However, the layout of wind poles cannot hinder the takeoff and landing of aircraft. Therefore, it must be adapted to local conditions and laid out on both sides of the runway at a certain distance from the runway at a certain interval. Therefore, how to use the data detected by this relatively sparse wind speed sensor array to obtain wind field information over the entire runway becomes a key issue. Summary of the invention

[0004] The purpose of the present invention is to provide a method, system, device and storage medium for generating airport area wind field streamlines and visualized wind shear products, so as to solve the problem that the wind field information over the entire runway cannot be obtained.

[0005] The embodiment of the present application is implemented by providing a method for generating wind field streamlines and visualized wind shear results in an airport area, including: obtaining triangular partitions;

[0006] Assume a wind field model, and solve the parameters of the assumed wind field model to obtain the wind speed vector of any point within the triangular partition;

[0007] Set the starting point of the wind field streamline and obtain the wind speed vector at the starting point;

[0008] Starting from each starting point of the triangular partition, the atmospheric particle trajectories on the wind field streamlines are obtained and the wind field streamlines are drawn;

[0009] Obtain and mark the maximum wind shear direction and maximum wind shear value of the triangular partition;

[0010] Obtain the divergence and vorticity within the triangular partition, calculate and plot the normalized product of the divergence and vorticity.

[0011] In some embodiments, the method for obtaining the triangular partition includes: connecting the stations on both sides of the runway into triangular partitions that are sequentially connected, have no gaps, and do not overlap with each other according to the nearest principle.

[0012] In some embodiments, assume a wind field model, and the formula of the wind field model is as follows:

[0013] ;

[0014] In the formula, represents the coordinates of any point in the triangular partition, that is, ; represents the wind speed vector at any point , expressed as , and respectively represent the wind speed components in the direction; and respectively represent the unit vectors in the direction.

[0015] In some embodiments, in the formula of the wind field model, the wind speed components in the direction are respectively:

[0016] ;

[0017] ;

[0018] In the formula, represents the first-order derivative of the wind speed component in the direction with respect to the abscissa variable ; represents the first-order derivative of the wind speed component in the direction with respect to the ordinate variable ; represents the first-order derivative of the wind speed component in the direction with respect to the abscissa variable ; represents the first-order derivative of the wind speed component The wind speed component in the direction at the set coordinate origin initial value, denotes the wind speed component in the direction at the set coordinate origin initial value.

[0019] In some embodiments, , , , , , the solution methods for the parameters include:

[0020] Construct a matrix according to the 3 wind pole positions at the 3 vertices of each triangular partition , , :

[0021] ;

[0022] Then solve to obtain , , , , , parameters:

[0023] , ;

[0024] In the formula, , , are respectively the wind speed components in the direction measured by the wind poles at the 3 vertices, , , are respectively the wind speed components in the direction measured by the wind poles at the 3 vertices.

[0025] In some embodiments, starting from each starting point of each triangular partition, use an iterative algorithm to obtain the atmospheric particle trajectories on each wind field streamline and draw the wind field streamlines; and / or

[0026] Obtain and mark the maximum wind shear direction and maximum wind shear value of each triangular partition; and / or

[0027] Obtain the parameters related to wind shear in each triangular partition, which characterize the degree of wind field convergence / divergence and the degree of wind field vorticity: divergence and vorticity, calculate the normalized product of divergence and vorticity, and draw and display the normalized product of divergence and vorticity in each triangular partition with corresponding color scales.

[0028] In some embodiments, the iterative algorithm includes:

[0029] Starting from a point P on the edge of the triangular partition, i.e., let: , where is the starting point position coordinate for drawing the wind field streamline, is the position coordinate of point P;

[0030] Obtain the displacement of the atmospheric particle to get a new position: , , , where n represents the current point and the current moment, n + 1 represents the next point and the next moment, and n starts from 1, represents the position where the atmospheric particle is at the th moment, represents the direction wind speed component at the current position, represents the direction wind speed component at the current position, represents the time interval of the new position of the atmospheric particle, represents the direction displacement of the atmospheric particle from the th moment after passing through to reach the th moment, represents the direction displacement of the atmospheric particle from the th moment after passing through to reach the th moment, represents the position where the atmospheric particle reaches at the th moment;

[0031] Judge whether the trajectory of the atmospheric particle reaches the boundary of the triangular partition. If yes, stop the iteration. If no, return to the previous step and continue to loop; and / or

[0032] The definition of the wind shear magnitude is , and the expression of the wind shear magnitude is:

[0033] ;

[0034] Where represents the azimuth angle, ranges from , represents the direction wind speed component The first-order derivative with respect to the abscissa variable ​ Represents The wind speed component in the direction with respect to the ordinate variable The first derivative of, Represents The wind speed component in the direction with respect to the abscissa variable The first derivative of, Represents The wind speed component in the direction with respect to the ordinate variable The first derivative of;

[0035] Find the azimuth angle that makes maximum, which is the wind direction with the maximum wind shear. At this time, the wind shear value is the maximum wind shear value; and / or

[0036] The method for obtaining the normalized product of divergence and vorticity includes:

[0037] The divergence and vorticity of each triangular partition. For example, for the th triangular partition, the calculation formulas for divergence and vorticity are as follows:

[0038] Divergence ;

[0039] Vorticity ;

[0040] In the formula, Represents The wind speed component in the direction with respect to the abscissa variable The first derivative of, Represents The wind speed component in the direction with respect to the ordinate variable The first derivative of, Represents The wind speed component in the direction with respect to the abscissa variable The first derivative of, Represents The wind speed component in the direction with respect to the ordinate variable The first derivative of;

[0041] Calculate the product of the divergence and vorticity of each triangular partition For the th triangular partition:

[0042]

[0043] Normalize the product of the divergence and vorticity to obtain the normalized product of divergence and vorticity , for the th triangular partition:

[0044] ;

[0045] In the formula, represents the maximum value of the divergence vorticity product in all triangular partitions, is the total number of triangular partitions.

[0046] Correspondingly, the embodiment of the present application also provides an airport area wind field streamline and visualization wind shear result generation system, including:

[0047] A triangular partition module, used to obtain triangular partitions;

[0048] A wind speed vector module, used to assume a wind field model, solve the parameters of the assumed wind field model, and obtain the wind speed vector of any point within the triangular partition;

[0049] A starting point wind speed vector module, used to set the starting point of the wind field streamline of the wind and obtain the starting point wind speed vector;

[0050] A wind field streamline drawing module, used to start from each starting point of the triangular partition, obtain the atmospheric particle trajectory on the wind field streamline, and draw the wind field streamline;

[0051] A maximum wind shear direction and maximum wind shear value module, used to obtain and mark the maximum wind shear direction and maximum wind shear value of the triangular partition;

[0052] A module for calculating and drawing the normalized product of divergence and vorticity, used to obtain the divergence and vorticity within the triangular partition, and calculate and draw the normalized product of divergence and vorticity.

[0053] Correspondingly, the embodiment of the present application also provides a computer device, including a storage and a processor. When the computer program stored in the storage is executed by the processor, the processor is caused to execute the steps of the above method.

[0054] Correspondingly, the embodiment of the present application also provides a computer-readable storage medium, storing a computer program. When the computer program is executed by a processor, the processor is caused to execute the steps of the above method.

[0055] In summary, due to the adoption of the above technical solutions, the beneficial effects of the present invention are:

[0056] The present invention first obtains triangular partitions, then assumes a wind field model, and then solves for the wind field model parameters based on the data of wind poles (wind speed sensors). Based on the solved wind field model parameters, the wind field streamlines in the coverage area of the wind pole array are constructed; based on the solved wind field model parameters, the wind direction with the maximum wind shear and the maximum wind shear value in each different triangular partition are further obtained, and then the maximum wind shear directions in each area are depicted, and the maximum wind shear values in each triangular partition are marked, as the characterization of the absolute magnitude and direction of the maximum wind shear in each partition; based on the solved wind field model parameters, the divergence and vorticity of each non-overlapping partition, which are closely related to the wind shear, are further obtained, and the normalized divergence and vorticity product are obtained and marked and displayed in a graphical form, as the characterization of the relative magnitude of the average wind shear in each partition. The wind field streamlines constructed and the wind shear products solved and drawn can provide the overall wind field trend and the quantitative evaluation of wind shear required for the safety guarantee of aircraft takeoff and landing at the airport.

[0057] The near-surface wind field trend and wind shear information over the airport are crucial for the safety of aircraft takeoff and landing. The wind field information obtained by the existing technologies is insufficient in terms of the overall trend display, visualization degree, and fineness, etc., and thus cannot provide sufficiently accurate and intuitive wind field information for pilots and air traffic control staff (such as tower controllers). The wind field streamline product constructed by the present invention can provide the overall wind field trend near the surface over the airport. The calculated and drawn map of the maximum wind shear magnitude and direction in each partition can provide the accurate value and direction information of the maximum wind shear in each partition. The calculated and drawn map of the divergence vorticity product can provide the relative magnitude comparison of the average wind shear degree in each partition. These products greatly improve the overall trend display, accuracy, and visualization degree of the airport wind field information for aircraft takeoff and landing, and can provide better reference for the safety guarantee work of aircraft takeoff and landing at the airport. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 It is a flowchart of the method for generating the wind field streamline and the visualized wind shear product in the airport area provided by the embodiment of the present invention;

[0059] Figure 2 It is a wind pole array and runway diagram provided by the embodiment of the present invention;

[0060] Figure 3 It is a diagram of the triangular partition division in the detection area provided by the embodiment of the present invention;

[0061] Figure 4 It is a single triangular partition diagram provided by the embodiment of the present invention;

[0062] Figure 5 It is a streamline starting point diagram of a single triangular partition provided by the embodiment of the present invention;

[0063] Figure 6It is the layout diagram of the wind pole array at a certain airport provided by the embodiment of the present invention;

[0064] Figure 7 It is the triangular partition diagram of the wind pole detection area at a certain airport provided by the embodiment of the present invention;

[0065] Figure 8 It is the wind speed vector diagram at the positions of each wind pole measured by the wind pole array at a certain airport provided by the embodiment of the present invention;

[0066] Figure 9 It is the wind field streamline diagram drawn for a certain airport provided by the embodiment of the present invention;

[0067] Figure 10 It is the diagram of the maximum wind shear direction and wind shear value marked for a certain airport provided by the embodiment of the present invention;

[0068] Figure 11 It is the normalized divergence vorticity product diagram drawn for a certain airport provided by the embodiment of the present invention. Specific embodiments

[0069] In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0070] The technical solution of this application is as follows:

[0071] As Figure 1 shown, in the first aspect, the embodiment of this application provides a method for generating wind field streamlines and visual wind shear products in the airport area, including:

[0072] S01. Obtain a triangular partition;

[0073] S02. Assume a wind field model, solve the parameters of the assumed wind field model, and obtain the wind speed vector at any point within the triangular partition;

[0074] S03. Set the starting point of the wind field streamline of the wind and obtain the wind speed vector at the starting point;

[0075] S04. Starting from each starting point of the triangular partition, obtain the atmospheric particle trajectory on the wind field streamline and draw the wind field streamline;

[0076] S05. Obtain and mark the maximum wind shear direction and the maximum wind shear value of the triangular partition;

[0077] S06. Obtain the divergence and vorticity within the triangular partition, calculate and draw the normalized product of the displayed divergence and vorticity.

[0078] The present invention first obtains triangular partitions, then assumes a wind field model, and then solves for the wind field model parameters based on the data of wind poles (wind speed sensors). Based on the solved wind field model parameters, the wind field streamline of the area covered by the wind pole array is constructed; based on the solved wind field model parameters, the wind direction with the largest wind shear and the maximum wind shear value in each different triangular partition are further obtained, and then the direction of the maximum wind shear in each area is depicted, and the maximum wind shear value of each triangular partition is marked, as a representation of the absolute magnitude and direction of the maximum wind shear in each partition; based on the solved wind field model parameters, the divergence and vorticity of each non-overlapping partition, which are closely related to the wind shear, are further obtained, and the normalized divergence and vorticity product are obtained and marked and displayed in a graphical form, as a representation of the relative magnitude of the average wind shear in each partition. The wind field streamline constructed and the wind shear products solved and drawn can provide the overall wind field trend and quantitative assessment of wind shear required for the safety guarantee of aircraft takeoff and landing at the airport.

[0079] The near-surface wind field trend and wind shear information over the airport are crucial for the safety of aircraft takeoff and landing. The wind field information obtained by the prior art is insufficient in terms of the overall trend display, visualization degree, and fineness, etc., so that it cannot provide sufficiently accurate and intuitive wind field information for pilots and air traffic control staff (such as tower controllers). The wind field streamline product constructed by the present invention can provide the overall wind field trend near the surface over the airport. The calculated and drawn map of the maximum wind shear magnitude and direction in each partition can provide the accurate value and direction information of the maximum wind shear in each partition. The calculated and drawn map of the divergence-vorticity product can provide the relative magnitude comparison of the average wind shear degree in each partition. These products greatly improve the overall trend display, accuracy, and visualization degree of the airport wind field information for aircraft takeoff and landing, and can provide better reference for the safety guarantee work of aircraft takeoff and landing at the airport.

[0080] In S01:

[0081] In some embodiments, the method for obtaining triangular partitions includes: connecting the stations on both sides of the runway in accordance with the nearest principle to form triangular partitions that are sequentially connected, have no gaps, and do not overlap with each other.

[0082] Exemplarily, as Figure 2 shown, 6 wind pole arrays with wind speed sensors are arranged around the aircraft runway.

[0083] Exemplarily, as Figure 3 shown, connect the stations on both sides of the runway in Figure 2 to form triangular partitions that are sequentially connected, have no gaps, and do not overlap with each other, as Figure 3 shown.

[0084] In S02:

[0085] It can be understood that the wind speed vector at any point within the triangular partition can be obtained based on the solved wind field model parameters.

[0086] Furthermore, assume a wind field model, and the formula of the wind field model is as follows:

[0087] ;

[0088] In the formula, represents the coordinates of any point in the triangular partition, that is, ; represents the wind speed vector at any point and is expressed as , and respectively represent the wind speed components in the direction and the direction; and respectively represent the unit vectors in the direction and the direction.

[0089] It can be understood that the wind at each point in the triangular partition is a vector.

[0090] As shown in Figure 4 , furthermore, each triangular partition has a total of three vertices, and the number of wind poles at each vertex is 1. The positions of the three wind poles are , , . The wind speed vectors measured by the wind poles at the three vertex positions are respectively denoted as , , , , , are respectively the wind speed components in the direction measured by the wind poles at the three vertices, , , are respectively the wind speed components in the direction measured by the wind poles at the three vertices.

[0091] It can be understood that the wind speed vector at the position can be expressed as or .

[0092] It can also be understood that , , contain 6 known horizontal and vertical wind speed components. Therefore, the wind field is first assumed to be a model with 6 unknown parameters.

[0093] Furthermore, in the formula of the wind field model, the wind speed components in the direction and the

[0094] direction are respectively:

[0095] ;

[0096] In the formula, represents the first derivative of the wind speed component in the direction with respect to the abscissa variable ; represents the first derivative of the wind speed component in the direction with respect to the ordinate variable ; represents the first derivative of the wind speed component in the direction with respect to the abscissa variable ; represents the first derivative of the wind speed component in the direction with respect to the ordinate variable ; represents the initial value of the wind speed component in the direction at the set coordinate origin ; represents the initial value of the wind speed component in the direction at the set coordinate origin ;

[0097] It can be understood that , , , , , are six parameters to be solved.

[0098] It can be understood that the wind speed component in the direction is the horizontal wind speed;

[0099] Furthermore, , , , , , The solution methods for the parameters include:

[0100] ​​​​​​According to the positions of the three wind vanes at the three vertices of each triangular partition , , construct a matrix:

[0101] ;

[0102] Then solve to obtain , , , , , parameters:

[0103] , ;

[0104] In the formula, , , are the wind speed components in the direction measured by the wind vanes at the three vertices respectively, , , are the wind speed components in the direction measured by the wind vanes at the three vertices respectively.

[0105] It can be understood that so far, the parameters of the wind field model described by the formulas and are solved. Substitute the formulas and back into , then the wind speed vector at any point within the triangular partition can be obtained from this wind field model.

[0106] In the above S03:

[0107] Exemplarily, as shown in Figure 5 , equally spaced points are selected on each side of the triangular partition as the starting points of the wind field streamlines, and then the wind speed vectors of each starting point are obtained. Taking point P as an example, assuming its coordinates are , substitute this coordinate into to obtain the wind speed vector of point P.

[0108] In the above S04:

[0109] Furthermore, starting from each starting point of each triangular partition, the atmospheric particle trajectories on each wind field streamline are obtained using an iterative algorithm, and the wind field streamlines are drawn.

[0110] Even further, the iterative algorithm includes:

[0111] S041. Starting from a point P on the edge of the triangular partition, i.e., let: , where is the starting point position coordinate for drawing the wind field streamline, is the position coordinate of point P;

[0112] S042. Obtain the displacement of the atmospheric particle to get a new position: , , , where n represents the current point and the current moment, n + 1 represents the next point and the next moment, n starts from 1, represents the position where the atmospheric particle is at the th moment, represents the -direction wind speed component at the current position, represents the -direction wind speed component at the current position, represents the time interval of the new position of the atmospheric particle, represents the -direction displacement of the atmospheric particle from the th moment to the th moment during this period, represents the -direction displacement of the atmospheric particle from the th moment to the th moment to the th moment during this period, represents the -direction displacement of the atmospheric particle from the th moment to the

[0113] S043. Determine whether the trajectory of the atmospheric particle reaches the boundary of the triangular partition. If "yes", the iteration stops; if "no", return to the previous step and continue the loop.

[0114] In the above S05:

[0115] Furthermore, obtain and mark the maximum wind shear direction and the maximum wind shear value of each triangular partition.

[0116] Even further, the definition of the wind shear magnitude is , and the expression for the wind shear magnitude is:

[0117] ;

[0118] where represents the azimuth angle, has a value range of , represents Wind speed component in the first derivative with respect to the abscissa variable , represents wind speed component in the first derivative with respect to the ordinate variable , represents wind speed component in the first derivative with respect to the abscissa variable , represents wind speed component in the first derivative with respect to the ordinate variable ;

[0119] Find the azimuth angle that makes the largest, which is the wind direction with the maximum wind shear. At this time, the wind shear value is the maximum wind shear value.

[0120] Exemplarily, draw arrowed lines at the centers of each triangular partition. The arrow direction represents the maximum wind shear direction of the triangular partition, and the line length is proportional to the magnitude of the maximum wind shear. At the same time, mark the maximum wind shear value of the triangular partition next to the arrow line, with the unit of 1 / second.

[0121] In the said S06:

[0122] Furthermore, obtain the parameters related to wind shear in each triangular partition, which characterize the degree of wind field convergence / divergence and the degree of wind field vorticity: divergence and vorticity, and calculate the normalized product of divergence and vorticity. Plot and display the normalized product of divergence and vorticity in each triangular partition with corresponding color scales.

[0123] It can be understood that the normalized product value of divergence and vorticity corresponds to different colors.

[0124] Even further, the method for obtaining the normalized product of divergence and vorticity includes:

[0125] The divergence and vorticity of each triangular partition. For example, for the th triangular partition, the calculation formulas for divergence and vorticity are as follows:

[0126] Divergence ;

[0127] Vorticity ;

[0128] In the formula, represents wind speed component in the first derivative with respect to the abscissa variable , represents Wind speed component in the first derivative of the ordinate variable, represents wind speed component in the first derivative of the abscissa variable, represents wind speed component in the first derivative of the ordinate variable;

[0129] Calculate the product of the divergence and vorticity of each triangular partition , for the th triangular partition:

[0130] ;

[0131] Normalize the product of the divergence and vorticity to obtain the normalized divergence and vorticity product , for the th triangular partition:

[0132] ;

[0133] In the formula, represents the maximum value of the divergence vorticity product in all triangular partitions, is the total number of triangular partitions.

[0134] It can be understood that in each triangular partition, the normalized product of the divergence and vorticity calculated above is plotted and displayed with corresponding color scales. ranges from 0 to 1, characterizing the average wind shear degree of each partition. Different values correspond to different colors.

[0135] In a second aspect, an embodiment of the present application provides a system for generating airport area wind field streamlines and visual wind shear products, including:

[0136] A triangular partition module for obtaining triangular partitions;

[0137] A wind speed vector module for assuming a wind field model, solving the parameters of the assumed wind field model, and obtaining the wind speed vector of any point within the triangular partition;

[0138] A starting point wind speed vector module for setting the starting point of the wind field streamline of the wind and obtaining the starting point wind speed vector;

[0139] A module for plotting wind field streamlines for obtaining the atmospheric particle trajectories on the wind field streamline starting from each starting point of the triangular partition and plotting the wind field streamlines;

[0140] The maximum wind shear direction and maximum wind shear value module is used to obtain and mark the maximum wind shear direction and maximum wind shear value of the triangular partition;

[0141] The module for calculating, plotting, and displaying the normalized product of divergence and vorticity is used to obtain the divergence and vorticity within the triangular partition, calculate, plot, and display the normalized product of divergence and vorticity.

[0142] In the triangular partition module:

[0143] In some embodiments, the method for obtaining the triangular partition includes: connecting the stations on both sides of the runway into triangular partitions that are sequentially connected, have no gaps, and do not overlap according to the nearest principle.

[0144] Exemplarily, as Figure 2 shown, 6 wind pole arrays with wind speed sensors are arranged around the aircraft runway.

[0145] Exemplarily, as Figure 3 shown, connect the stations on both sides of the runway in Figure 2 into triangular partitions that are sequentially connected, have no gaps, and do not overlap, as Figure 3 shown.

[0146] In the wind speed vector module:

[0147] It can be understood that the wind speed vector at any point within the triangular partition can be obtained based on the solved wind field model parameters.

[0148] Furthermore, assume a wind field model, and the formula of the wind field model is as follows:

[0149] ;

[0150] In the formula, represents the coordinates of any point in the triangular partition, that is, ; represents the wind speed vector at any point , expressed as , and respectively represent direction and direction wind speed components; and respectively represent direction and direction unit vectors.

[0151] It can be understood that the wind at each point in the triangular partition is a vector.

[0152] As Figure 4As shown, further, each triangular partition has a total of three vertices, and the number of wind vanes at each vertex is 1. The positions of the three wind vanes are , , . The wind speed vectors measured by the wind vanes at the three vertex positions are respectively denoted as , , . , , are respectively the wind speed components in the direction measured by the wind vanes at the three vertices. , , are respectively the wind speed components in the direction measured by the wind vanes at the three vertices.

[0153] It can be understood that the wind speed vector at the position can be expressed as or .

[0154] Furthermore, it can be understood that , , contain six known horizontal and vertical wind speed components. Therefore, the wind field is first assumed to be a model with six unknown parameters.

[0155] Further, in the formula of the wind field model, the wind speed components in the direction and

[0156] ;

[0157] ;

[0158] In the formula, represents the first derivative of the wind speed component in the direction with respect to the abscissa variable . represents the first derivative of the wind speed component in the direction with respect to the ordinate variable . represents the first derivative of the wind speed component in the direction with respect to the abscissa variable . represents the first derivative of the wind speed component in the direction with respect to the ordinate variable . represents The wind speed component in the direction at the set coordinate origin initial value at represents the wind speed component in the direction at the set coordinate origin initial value at.

[0159] It can be understood that , , , , , are 6 parameters to be solved.

[0160] It can be understood that the wind speed component in the direction is the horizontal wind speed;

[0161] Furthermore, , , , , , The solution methods for the parameters include:

[0162] According to the 3 wind vane positions at the 3 vertices of each triangular partition , , construct a matrix:

[0163] ;

[0164] Then the solution obtains , , , , , parameters:

[0165] , ;

[0166] In the formula, , , are the wind speed components in the direction measured by the wind vanes at the 3 vertices respectively, , , are the wind speed components in the direction measured by the wind vanes at the 3 vertices respectively.

[0167] It can be understood that so far, the formula and the wind field model parameters described are solved, and substituting the formulas and back into , the wind speed vector at any point within the triangular partition can be obtained from this wind field model.

[0168] In the starting point wind speed vector module described:

[0169] Exemplarily, as shown in Figure 5 , points are selected at equal intervals on each side of the triangular partition as the starting points of the wind field streamlines, and then the wind speed vectors of each starting point are obtained. Taking point P as an example, assuming its coordinates are , substituting these coordinates into to obtain the wind speed vector of point P.

[0170] In the wind field streamline drawing module described:

[0171] Furthermore, starting from each starting point in each triangular partition, the atmospheric particle trajectories on each wind field streamline are obtained using an iterative algorithm, and the wind field streamlines are drawn.

[0172] Even further, the iterative algorithm includes:

[0173] Starting from a point P on the edge of the triangular partition, i.e., let: , where is the position coordinate of the starting point for drawing the wind field streamline, is the position coordinate of point P;

[0174] Calculate the displacement of the atmospheric particle to obtain a new position: , , , where n represents the current point and the current time, n + 1 represents the next point and the next time, n starts from 1, represents the position where the atmospheric particle is located at the th time, represents the wind speed component in the direction at the current position, represents the wind speed component in the direction at the current position, represents the time interval of the new position of the atmospheric particle, represents the displacement in the direction made by the atmospheric particle from the th time to the th time during this period, direction, represents the displacement in the direction made by the atmospheric particle from the The displacement of atmospheric particles during the period up to the th moment in the direction, represents the position where the atmospheric particles reach at the th moment;

[0175] Judge whether the trajectory of the atmospheric particles reaches the boundary of the triangular partition. If "yes", the iteration stops; if "no", return to the previous step and continue the loop.

[0176] In the said maximum wind shear direction and maximum wind shear value module:

[0177] Furthermore, calculate and mark the maximum wind shear direction and maximum wind shear value of each triangular partition.

[0178] Even further, the definition of the wind shear magnitude is , and the expression of the wind shear magnitude is:

[0179] ;

[0180] In the formula, represents the azimuth angle, has a value range of , represents the wind speed component in the first-order derivative with respect to the abscissa variable , represents the wind speed component in the first-order derivative with respect to the ordinate variable , represents the wind speed component in the first-order derivative with respect to the abscissa variable , represents the wind speed component in the first-order derivative with respect to the ordinate variable ;

[0181] Find the azimuth angle that makes the largest, which is the wind direction with the maximum wind shear. At this time, the wind shear value is the maximum wind shear value.

[0182] Exemplarily, draw arrowed lines at the centers of each triangular partition. The arrow direction represents the maximum wind shear direction of the triangular partition, and the line length is proportional to the magnitude of the maximum wind shear. At the same time, mark the maximum wind shear value of the triangular partition next to the arrow line, with the unit of 1 / second.

[0183] In the said calculating and drawing the display of the normalized product of divergence and vorticity module:

[0184] Furthermore, the parameters closely related to wind shear in each triangular partition, which characterize the degree of wind field convergence and divergence and the degree of wind field vortex, are obtained: divergence and vorticity, and the normalized product of divergence and vorticity is calculated. The normalized product of divergence and vorticity is plotted and displayed in each triangular partition using the corresponding color scale.

[0185] It can be understood that the normalized product values ​​of divergence and vorticity correspond to different colors.

[0186] Furthermore, the method for obtaining the normalized product of divergence and vorticity includes:

[0187] The divergence and vorticity of each triangular partition, such as The calculation formulas of divergence and vorticity are as follows:

[0188] Divergence ;

[0189] Vorticity ;

[0190] In the formula, express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of ;

[0191] Find the product of the divergence and vorticity of each triangular partition , for Triangular partitions:

[0192] ;

[0193] The product of divergence and vorticity is normalized to obtain the normalized divergence and vorticity product , for Triangular partitions:

[0194] ;

[0195] In the formula, represents the maximum value of the divergence vorticity product among all triangular partitions, is the total number of triangular partitions.

[0196] It can be understood that in each triangular partition, the normalized product of the divergence and vorticity calculated previously is plotted and displayed using the corresponding color scale. ranges from 0 to 1, representing the average wind shear degree of each partition. Different values correspond to different colors.

[0197] In a third aspect, the present application provides a computer device, including a storage and a processor. When the computer program stored in the storage is executed by the processor, the processor is caused to execute the steps of the method for generating the streamline of the wind field in the airport area and the visualized wind shear product as described above.

[0198] Among them, the computer device can be a desktop computer, a notebook, a palm computer, a cloud server, and other computing devices. The computer device can perform human-computer interaction with the user through a keyboard, a mouse, a remote control, a touchpad, a voice control device, and other means.

[0199] The memory includes at least one type of readable storage medium, and the readable storage medium includes flash memory, a hard disk, a multimedia card, a card-type memory (such as an SD or D interface display memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a magnetic disk, an optical disk, etc. In some embodiments, the memory can be an internal storage unit of the computer device, such as the hard disk or memory of the computer device. In other embodiments, the memory can also be an external storage device of the computer device, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) card, a Flash Card, etc. equipped on the computer device. Of course, the memory can also include both the internal storage unit and the external storage device of the computer device. In this embodiment, the memory is commonly used to store the operating system installed on the computer device and various application software, such as the program code of the method for generating the streamline of the wind field in the airport area and the visualized wind shear product. In addition, the memory can also be used to temporarily store various data that have been output or will be output.

[0200] In some embodiments, the processor may be a central processing unit (CPU), a controller, a microcontroller, a microprocessor, or other data processing chips. The processor is generally used to control the overall operation of the computer device. In this embodiment, the processor is used to run the program code stored in the memory or process data, such as running the program code of the method for generating the streamline of the wind field in the airport area and the visualized wind shear product.

[0201] Fourthly, the present application provides a computer-readable storage medium storing a computer program, which, when executed by a processor, causes the processor to execute the steps of the method for generating the streamline of the wind field in the airport area and the visualized wind shear product as described above.

[0202] Wherein, the computer-readable storage medium stores an interface display program, which can be executed by at least one processor to cause the at least one processor to execute the steps of the method for generating the streamline of the wind field in the airport area and the visualized wind shear product as described above.

[0203] Application Example

[0204] This application example adopts the layout of the wind pole array with wind speed sensors at a certain domestic airport. As Figure 6 , taking the position of wind pole 1 as the coordinate origin, the coordinates (in meters) of the 9 wind poles are respectively: R1 = [0, 0], R2 = [0, 1721.8], R3 = [0, 3370.3], R4 = [2188.9, -2294.7], R5 = [-2590.6, -2294.7], R6 = [-2590.6, 565.6], R7 = [-2590.6, 6315.0], R8 = [0, 6315.0], R9 = [2188.9, 6315.0].

[0205] Triangulate the entire area covered by the wind poles into triangular partitions, as Figure 7 shown.

[0206] The wind speed vectors measured on wind poles 1 - 9 are as shown by the red arrows in Figure 8 . The specific wind speed (in meters per second) vectors measured by each wind pole at a certain moment are:

[0207] = [-9.8524, 2.6399], = [-5.3623, -4.4995],

[0208] = [-3.6773, -7.2171], =[-4.9790, 4.4831],

[0209] =[-2.8531, -0.9270], =[-2.0138, -2.4868],

[0210] =[1.6771, -1.2638], =[-12.4195, -6.8842],

[0211] =[-1.5045, -1.9965];

[0212] Obtain the wind speed vector of any point within the triangular partition.

[0213] Obtain the wind speed vectors of the starting points of each triangular partition.

[0214] Based on the iterative algorithm of the present invention, obtain the atmospheric particle trajectories, draw streamlines, and reflect the overall trend of the wind field in the airport area. The results are as Figure 9 shown.

[0215] Obtain the maximum wind shear direction and maximum wind shear value of each triangular partition, draw and mark them, as Figure 10 shown.

[0216] Take the normalized divergence vorticity product of each triangular partition, mark the values in each triangular partition, and draw them with corresponding color scales to characterize the average wind shear degree of each partition, as Figure 11 shown; Figure 11 Among them, the color scale ranges from dark blue to yellow, and different colors represent different normalized divergence vorticity product values.

[0217] The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent replacements, and improvements made within the spirit and principles of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for generating wind field streamlines and visualized wind shear products in an airport area, characterized in that: include: Get triangle partitions; Assume a wind field model, solve the assumed wind field model parameters, and obtain the wind speed vector of any point within the triangular partition; Set the starting point of the wind field streamline and obtain the wind speed vector at the starting point; Starting from each starting point of each triangular partition, an iterative algorithm is used to obtain the atmospheric particle trajectory on each wind field streamline and draw the wind field streamline. The iterative algorithm includes: Starting from a point P on the edge of the triangular partition, let: , where The coordinates of the starting point for drawing the wind field streamlines. is the position coordinate of point P; Calculate the displacement of atmospheric particles and obtain the new position: , , , where n represents the current point and the current time, n+1 represents the next point and the next time, and n starts from 1. Indicates The position of the atmospheric particles at a given moment, Indicates the current location Directional wind speed component, Indicates the current location Directional wind speed component, The time interval representing the new position of the atmospheric particle, Indicates that from The moment begins to pass Arrive at The atmospheric particles make Direction displacement, Indicates that from The moment begins to pass Arrive at The atmospheric particles make Direction displacement, Indicates The position where the atmospheric particles reach at a certain moment; Determine whether the atmospheric particle trajectory reaches the boundary of the triangular partition. If yes, the iteration stops; if not, return to the previous step to continue the cycle; Obtain and mark the maximum wind shear direction and maximum wind shear value of the triangular partition; Find the divergence and vorticity within a triangular partition, and calculate and plot the normalized product of divergence and vorticity.

2. The method for generating airport area wind field streamlines and visualized wind shear products according to claim 1, characterized in that: The method for obtaining the triangular partition includes: according to the nearest principle, connecting the stations on both sides of the runway into triangular partitions that are connected in sequence, have no gaps and do not overlap each other.

3. The method for generating airport area wind field streamlines and visualized wind shear products according to claim 2, characterized in that: Assuming the wind field model, the formula of the wind field model is as follows: ； In the formula, Represents any point in a triangle partition The coordinates of ; Indicates any point The wind speed vector at is expressed as , and Respectively Direction and Directional wind speed component; and Respectively Direction and The unit vector of the direction.

4. The method for generating airport area wind field streamlines and visualized wind shear products according to claim 3, characterized in that: In the formula of the wind field model, Direction and The wind speed components in the directions are: ； ； In the formula, express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of express Directional wind speed component At the set coordinate origin The initial value of express Directional wind speed component At the set coordinate origin The initial value of .

5. The method for generating airport area wind field streamlines and visualized wind shear products according to claim 4, characterized in that: , , , , , The methods for solving the parameters include: According to the 3 wind pole positions at the 3 vertices of each triangle partition , , Construct the matrix: ； Then the solution is , , , , , parameter: , ； In the formula, , , They are measured by the wind poles at the three vertices. Directional wind speed component, , , They are measured by the wind poles at the three vertices. Directional component of wind speed.

6. The method for generating airport area wind field streamlines and visualized wind shear products according to claim 2, characterized in that: Obtain and annotate the maximum wind shear direction and maximum wind shear value for each triangular subarea; and / or Obtain the parameters closely related to wind shear in each triangular partition, which characterize the degree of wind field convergence and divergence and the degree of wind field vortex: divergence and vorticity, and calculate the normalized product of divergence and vorticity. Draw and display the normalized product of divergence and vorticity in each triangular partition using the corresponding color scale.

7. The method for generating airport area wind field streamlines and visualized wind shear products according to claim 6, characterized in that: The wind shear magnitude is defined as , the wind shear magnitude is expressed as: ； In the formula, represents the azimuth, The value range is , express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of ; Find out The maximum azimuth is the wind direction with the maximum wind shear, and the wind shear value at this time is the maximum wind shear value; and / or The methods for obtaining the normalized product of divergence and vorticity include: The divergence and vorticity of each triangular partition, such as The calculation formulas of divergence and vorticity are as follows: Divergence ; Vorticity ; In the formula, express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of express Directional wind speed component For the horizontal axis variable The first derivative of express Directional wind speed component For the vertical axis variable The first derivative of ; Find the product of the divergence and vorticity of each triangular partition , for Triangular partitions: ； The product of divergence and vorticity is normalized to obtain the normalized divergence and vorticity product , for Triangular partitions: ； In the formula, represents the maximum value of the divergence vorticity product among all triangular partitions, is the total number of triangle partitions.

8. Airport area wind field streamline and visual wind shear product generation system, characterized by: include: A triangle partitioning module, used to obtain triangle partitioning; The wind speed vector module is used to assume a wind field model and solve the assumed wind field model parameters to obtain the wind speed vector of any point within the triangular partition; The starting point wind speed vector module is used to set the starting point of the wind field streamline and obtain the starting point wind speed vector; The module for drawing wind field streamlines starts from each starting point of each triangular partition, uses an iterative algorithm to obtain the atmospheric particle trajectory on each wind field streamline, and draws the wind field streamlines. The iterative algorithm includes: Starting from a point P on the edge of the triangular partition, let: , where The coordinates of the starting point for drawing the wind field streamlines. is the position coordinate of point P; Calculate the displacement of atmospheric particles and obtain the new position: , , , where n represents the current point and the current time, n+1 represents the next point and the next time, and n starts from 1. Indicates The position of the atmospheric particles at a given moment, Indicates the current location Directional wind speed component, Indicates the current location Directional wind speed component, The time interval representing the new position of the atmospheric particle, Indicates that from The moment begins to pass Arrive at The atmospheric particles make Direction displacement, Indicates that from The moment begins to pass Arrive at The atmospheric particles make Direction displacement, Indicates The position where the atmospheric particles reach at a certain moment; Determine whether the atmospheric particle trajectory reaches the boundary of the triangular partition. If yes, the iteration stops. If not, return to the previous step to continue the cycle. The maximum wind shear direction and maximum wind shear value module is used to obtain and mark the maximum wind shear direction and maximum wind shear value of the triangular partition; The module for calculating and drawing the normalized product of divergence and vorticity is used to obtain the divergence and vorticity within the triangular partition, and calculate and draw the normalized product of divergence and vorticity.

9. Computer device, characterized in that The invention comprises a storage and a processor, wherein the storage stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that A computer program is stored, and when the computer program is executed by a processor, the processor is caused to perform the steps of the method according to any one of claims 1 to 7.

Images