CN111399084B

CN111399084B - High-altitude rapid flow extraction method based on three-dimensional wind field data

Info

Publication number: CN111399084B
Application number: CN202010136588.9A
Authority: CN
Inventors: 陈曦; 聂建强; 严勇杰; 田云钢; 肖雪飞; 刘岩; 莫海健
Original assignee: CETC 28 Research Institute
Current assignee: CETC 28 Research Institute
Priority date: 2020-03-02
Filing date: 2020-03-02
Publication date: 2021-12-24
Anticipated expiration: 2040-03-02
Also published as: CN111399084A

Abstract

The invention provides a high altitude torrent extraction method based on three-dimensional wind field data. According to the meteorological definition of the tropical east wind torrent, the information of the length, the width, the height, the torrent axis, the torrent core, the horizontal and vertical wind shear and the like of the east wind torrent are obtained. According to the meteorology definition of the western wind torrent, the information of the length, the width, the height, the size and the intensity position of the torrent shaft, the horizontal wind shear, the vertical wind shear and the like of each type of western wind torrent are obtained. The invention extracts the high-altitude torrent information which has important influence on high-altitude flight from the complex three-dimensional wind field by using an objective method, provides a uniform high-altitude flight meteorological situation for all parties, and can effectively ensure the safety and the economy of air transportation.

Description

High-altitude rapid flow extraction method based on three-dimensional wind field data

Technical Field

The invention belongs to the field of aeronautical weather, and particularly relates to a high-altitude torrent extraction method based on three-dimensional wind field data.

Background

Currently, high-altitude wind data for aeronautical meteorological service comes from a WAFS (world Area weather System) system provided by the international civil aviation organization and the world meteorological organization, which is three-dimensional wind field data consisting of transverse wind (V) and latitudinal wind (U) on dozens of different equal-pressure surfaces. However, this data is not convenient for being directly used by all relevant interests in the air transportation industry, and the position and attribute information of the high-altitude torrent which has important influence on the flight is difficult to judge from the three-dimensional wind field. The high altitude rush current is a long and narrow high speed airflow band with wind speed more than 30m/s in the height range of about 10km or so of the upper troposphere, and has strong horizontal and vertical wind shear. The length of the rush current can reach thousands of kilometers or more, the width can reach hundreds of kilometers to thousands of kilometers, and the thickness can reach several kilometers. The high altitude rush current has important influence on aviation flight in the aspects of economy and safety. When flying in a torrent, the ground speed of the airplane can be increased, the navigation time can be shortened, and the fuel can be saved. The ground speed of the airplane can be reduced by the counter-rush flight, the flight time is increased, and more fuel is consumed. In addition, because strong wind shear exists near the high-altitude torrent, a large drift current can be generated when the aircraft flies through the high-altitude torrent, the influence on navigation calculation and heading keeping is caused, and the jolt of the aircraft can be caused, so that the aviation safety accident is caused. Therefore, the intensity and the spatial distribution information of the high-altitude torrent are mastered, and the method has important significance for aviation flight.

According to related patents and documents at home and abroad, at present, the aviation meteorological service at home and abroad only provides three-dimensional wind field information, and an algorithm technology for extracting high altitude torrent intensity and position information from three-dimensional wind field data is not provided.

Disclosure of Invention

At present, an aeronautical meteorological service department provides high-altitude wind information service for aviation transportation by providing warp wind and weft wind data on different equal pressure surfaces.

The invention particularly provides a high-altitude torrent extraction method based on three-dimensional wind field data, which comprises the following steps of:

step 1, acquiring three-dimensional wind direction and wind speed information on a flight altitude layer, and then extracting a tropical eastern torrent range and a tropical western torrent range;

step 2: extracting properties of the Dongfeng torrent;

and step 3: and extracting the properties of the western wind torrent.

The step 1 comprises the following steps:

step 1-1: acquiring three-dimensional wind direction and wind speed information on a flight height layer;

step 1-2: extracting the tropical eastern wind torrent;

step 1-3: and extracting the western acute flow.

The step 1-1 comprises the following steps:

step 1-1-1, performing height layer conversion to convert the vertical coordinate system from an equal-pressure surface to an equal-height surface, wherein the conversion formula is a Laplace pressure formula

Wherein P is₁,z₁T is the pressure, height and temperature of the reference surface at zero degrees Celsius, P₂,z₂Respectively target pressure and height, P₁＝1013.24，z ₁0, t 0, and the above formula is simplified to Z₂＝18400×log(1013.24/P₂) Converting the constant pressure surface to the nearest flying height according to a simplified formula;

step 1-1-2, converting the warp wind and weft wind data into wind direction and wind speed data: calculating the wind direction and the wind speed according to the warp wind and weft wind data on each height layer, wherein the specific formula is as follows:

wherein W_sIs the wind speed, W_dAnd obtaining three-dimensional wind direction and speed information on the flight height layer through the calculation.

The step 1-2 comprises the following steps:

the method comprises the steps of reducing the extraction range of tropical east wind rapids to a height range of 9000-16000 m, locking a latitude range between 15 DEG S and 30 DEG N, locking a longitude range between 0 DEG and 140 DEG E, carrying out binary labeling processing on three-dimensional wind direction and wind speed information according to a wind speed threshold of 30m/S, marking grid points with wind speeds of more than or equal to 30m/S as 1, marking grid points with wind speeds of less than 30m/S as 0, carrying out binary labeling processing on the three-dimensional wind speed and wind speed information according to the wind directions, marking grid points with wind directions of more than 0 DEG and less than 180 DEG as 1, and marking grid points with wind directions of less than 0 as 0; and extracting all grid points with wind speed and wind direction labels equal to 1 according to the space range of the east wind torrent, namely the tropical east wind torrent.

The steps 1-3 comprise:

and carrying out binary labeling processing on the three-dimensional wind direction and wind speed information according to a wind speed threshold value of 30m/s, marking grid points with wind speeds of more than or equal to 30m/s as 1 and marking grid points with wind speeds of less than 30m/s as 0, carrying out binary labeling processing on the three-dimensional wind speed and wind speed information according to the wind direction, marking grid points with wind directions of more than 180 degrees and less than 360 degrees as 1, and marking grid points which are not in the range as 0. The altitude range is locked between 5 and 13.5km, and all grid points with wind speed and wind direction labels equal to 1 are extracted, i.e. all west wind rapids.

The step 2 comprises the following steps:

step 2-1: analyzing a longitudinal section; starting from 0-degree E to 140-degree E, traversing the meridional sections within the range of 15-degree S to 30-degree N one by taking the longitude resolution of three-dimensional wind direction and wind speed information as a step length, judging whether the east wind torrent exists on each section, if the latitude range of the east wind torrent of one section is less than 500km, rejecting the east wind torrent, and setting the east wind torrent attribute label of the section as a False value False; if the latitude range of the east wind torrent is more than or equal to 500km, setting the east wind torrent attribute label of the section to be a True value True; for each radial section with the east wind torrent attribute label of True, calculating the width and height of the radial section, the central position of the minimum circumscribed rectangle, the maximum wind speed value and the maximum wind speed value position information according to the minimum circumscribed rectangle of the east wind torrent, and calculating according to 4 lattice points around the lattice point of the maximum wind speed value to obtain the horizontal wind shear

And vertical wind shear

Wherein W_s(max)Is the maximum wind speed value, W_s(s)Is the wind speed value, W, of the south grid of the maximum wind speed grid point_s(n)Is the wind speed value, W, of the grid on the north side of the maximum wind speed grid point_s(u)Is the wind speed value, W, of the grid above the maximum wind speed grid point_s(d)Is the wind speed value of the grid below the maximum wind speed grid point;

step 2-2: acquiring integral attribute information of the torrent: the method comprises the steps that only one east wind torrent is provided, starting from 0 DEG E to 140 DEG E, recording the start longitude and the end longitude of a continuous uninterrupted longitude sequence with the east wind torrent attribute label of True value True, and calculating the length L of the east wind torrent to be 110km multiplied by delta theta according to the difference delta theta of the start longitude and the end longitude; for the width, height, central position, maximum wind speed value position and horizontal wind shear S of all longitudinal sections in the starting and stopping longitude range_hAnd vertical wind shear S_vAveraging, and calculating to obtain the relevant attribute values of the whole east wind torrent, including width, height, central point position, horizontal and vertical wind shear;

step 2-3: acquiring attribute information of the torrent axis: combining the maximum wind speed values on all the warp-wise sections within the start-stop longitude range of the east wind torrent obtained in the step 2-2 into an ordered array, and recording the maximum wind speed values and the position information of the maximum wind speed values in the ordered array, namely the torrent core information on the torrent shaft; the maximum point on the east wind torrent starting radial profile, the position point of the torrent core and the maximum point on the east wind torrent ending radial profile form a torrent axis of the tropical east wind;

step 2-4: acquiring torrent horizontal distribution information: according to the height of the east wind torrent and the position information of the central point obtained in the step 2-2, extracting the position profiles of different wind speed thresholds on each flight height layer influenced by the east wind torrent, and acquiring the horizontal space distribution information of the east wind torrent on different height layers, wherein the specific method comprises the following steps: the method comprises the steps of using a contour function in a matplotlib library based on python language to respectively calculate closed position contour sequences with thresholds of 30, 35, 40 and 45m/s and different wind speed grades, using real numbers in 0-Nx (Ny) as coordinates to carry out position recording, wherein Nx and Ny are grid total numbers in x and y directions respectively, and carrying out coordinate conversion according to start and stop longitudes on a plane grid to obtain a longitude and latitude coordinate sequence of a position contour, wherein the calculation formula is as follows:

wherein (X, Y) is a position point coordinate calculated by a contour function, (J, W) is a longitude and latitude coordinate of the position point after conversion, and J_sAnd J_eRespectively start and end longitude, W, of the plane grid_sAnd W_eRespectively the starting latitude and the stopping latitude of the plane grid.

In step 2-3, if the wind speed value of one point in the ordered array is larger than the wind speed values of two adjacent points, the wind speed value of the point is the maximum wind speed value.

The step 3 comprises the following steps:

step 3-1: classifying the western wind torrent; starting from 0 degrees, going around the earth from west to east, ending at 0 degrees, traversing warp sections in a range from 90 degrees S to 90 degrees N one by taking longitude resolution of three-dimensional wind direction and wind speed information as a step length, if the latitude range of one torrent on one warp section is less than 500km, eliminating the torrent on the warp section, then reserving the first 4 torrents with the largest width in the remaining torrents according to the width of each torrent, eliminating the remaining torrents, ensuring that the number of the torrents on each section is equal to or less than 4, and finally classifying each torrent on the warp section by using a west torrent classification matrix. The contents of the first row and the first column of a western wind torrent classification matrix of a certain meridional section are used as titles, the first row lists 4 types of western wind torrents, and the first column lists torrent sequences in the meridional section. The numerical value in the classification matrix table of the western wind torrent indicates the number of grid points in a certain torrent range on the meridional section, which fall in the latitude distribution range of the four types of torrents. Selecting the torrent with the maximum number of grid points in the range of the four torrents as the type of western torrent, and if the maximum number of grid points is multiplied by the spacing distance between grids to be less than 1000km, marking the attribute label of the type of western torrent of the section as a False value False; if the maximum number of grid points multiplied by the spacing distance between the grids is greater than or equal to 1000km, then such west wind torrent attribute label of the profile is marked as True value True. For each warp-wise section from 0 degrees to 0 degrees around the earth from west to east, performing torrent classification by using the west torrent classification matrix, and for the warp-wise sections with True-value attribute labels of each category of west torrent, calculating the width and height of the warp-wise section, the central position of the minimum circumscribed rectangle, the maximum wind speed value and the maximum wind speed value position, and horizontal wind shear and vertical wind shear information according to the minimum circumscribed rectangle of the torrent by using the method in the step 2-1;

step 3-2: acquiring overall attribute information of various western wind torrents:

aiming at each type of western wind torrent, selecting a section sequence of which the attribute label is continuously True value True, recording the start-stop longitude of each continuous section sequence, and calculating by using the method of the step 2-2 to obtain the relevant attribute value of each torrent of each type of western wind torrent, including length, width and height;

step 3-3: acquiring attribute information of various western wind torrent axes:

calculating to obtain the torrent axis attribute information of each type of western torrent by using the method in the step 2-3 aiming at each type of western torrent;

step 3-4: acquiring various horizontal distribution information of the western wind torrent:

and (3) obtaining all flight height layers within the vertical range of each type of western wind torrent according to the upper and lower height boundaries of each type of western wind torrent, and obtaining the horizontal distribution information of each type of western wind torrent on each flight height layer by using the method in the step 2-4.

Has the advantages that: the high altitude rush current has important influence on the flight, the ground speed can be increased by the downwind flight, the fuel consumption is reduced, and the ground speed is reduced by the upwind flight, so that the fuel consumption is increased. In addition, the flying problem that the flying comfort and safety of the airplane are affected by airplane bumping and the like easily occurs when the airplane passes through the high altitude torrent axis to fly. The conventional aeronautical meteorological service product only contains three-dimensional warp wind and weft wind data and does not have high-altitude torrent information. The method calculates the high-altitude torrent information such as torrent intensity, torrent axis position, horizontal distribution of torrent layers and the like from the three-dimensional warp wind and weft wind data, provides accurate and visual high-altitude torrent information for each participant of the air transportation, provides a reference basis for planning an air route, and can effectively guarantee the safety and the economy of the air transportation.

Drawings

The foregoing and/or other advantages of the invention will become further apparent from the following detailed description of the invention when taken in conjunction with the accompanying drawings.

FIG. 1 is a flow chart of the method of the present invention.

Fig. 2 is a schematic diagram of extraction of east wind torrent meridional section information.

FIG. 3 is a schematic diagram of horizontal and vertical wind shear calculations.

Detailed Description

The invention mainly comprises three parts: torrent extraction, tropical east torrent attribute extraction and west torrent attribute extraction. Fig. 1 shows a flow chart of the algorithm of the present invention. Each part comprises the following specific steps:

1. torrent extraction

Step 1-1: obtaining three-dimensional wind direction and wind speed information on flight altitude layer

The extraction of high altitude torrent needs to convert the data of the warp wind V and the weft wind U on the isobaric surface into the data of the wind direction and the wind speed on the flight height layer. First, a height layer conversion is required to convert the vertical coordinate system from an iso-pressure surface to an iso-height surface. Converting the formula to a Laplace high formula

Wherein P is₁,z₁Is the pressure and height of the reference surface, and the reference surface for high-altitude flight is the standard sea level (the air pressure is 1013.25hPa) at the temperature of 0 ℃, so the formula can be simplified into Z₂＝18400×log(1013.24/P₂). The isobaric surface is converted to the flying height nearest to the isobaric surface according to a simplified formula, and the isobaric surface is shown in the table 1And mapping with flying height.

TABLE 1

Isobaric surface layer (Unit: hPa)	Flight altitude layer (Unit: hft)
		700	FL100
600	FL140
		500	FL180
450	FL210
		400	FL240
350	FL270
		300	FL300
275	FL320
		250	FL340
225	FL370
		200	FL390
175	FL420
		150	FL450
125	FL480
		100	FL520

Secondly, the warp and weft wind data need to be converted into wind direction and wind speed data. Calculating the wind direction and the wind speed according to the warp wind and weft wind data on each height layer, wherein the specific formula is as follows:

wherein W_sIs the wind speed, W_dIs the wind direction, U is the latitudinal wind, and V is the latitudinal wind. Through the calculation, three-dimensional wind direction and wind speed information on a flight height layer can be obtained.

Step 1-2: extracting tropical eastern wind torrent

The tropical east torrent is a strong east torrent zone on the high pressure south side of the upper subtropical zone of the troposphere, with a latitude in the range of approximately 5-20N, a longitude in the range of approximately 30-130E, and an altitude in the range of approximately 13-15 km. According to the meteorological definition, the invention reduces the extraction range of the tropical east wind rapids to the height range of 9000 to 16000m, locks the latitude range between 15 DEG S and 30 DEG N and locks the longitude range between 0 DEG E and 140 DEG E by considering the change of the spatial distribution. And carrying out binary labeling processing on the three-dimensional wind direction and wind speed information according to a wind speed threshold value of 30m/s, marking grid points with wind speeds of more than or equal to 30m/s as 1 and marking grid points with wind speeds of less than 30m/s as 0, carrying out binary labeling processing on the three-dimensional wind speed and wind speed information according to the wind direction, marking grid points with wind directions of more than 0 degrees and less than 180 degrees as 1, and marking grid points which are not in the range as 0. And extracting all grid points with wind speed and wind direction labels equal to 1 according to the space range of the east wind torrent, namely the tropical east wind torrent.

Step 1-3: extracting the acute stream of the west wind

The western wind torrent refers to a strong western wind central area of a prevailing western wind zone near the top of a convection layer, the horizontal width of the strong western wind central area is approximately 800-1000km, the vertical thickness of the strong western wind central area is approximately 6-10km, the length of the strong western wind central area is ten thousand kilometers, and even the strong western wind central area surrounds the earth for a circle. According to the difference of geographical position, the west wind torrent can be divided into subtropical west wind torrent and polar front torrent (temperate west wind torrent), and the south and north hemispheres are all existed. The approximate latitudes of the subtropical western currents range from 15 ° n(s) to 40 ° n(s), and the latitudes of the extreme currents range from 40 ° n(s) to 70 ° n(s). And carrying out binary labeling processing on the three-dimensional wind direction and wind speed information according to a wind speed threshold value of 30m/s, marking grid points with wind speeds of more than or equal to 30m/s as 1 and marking grid points with wind speeds of less than 30m/s as 0, carrying out binary labeling processing on the three-dimensional wind direction and wind speed information according to the wind direction, marking grid points with wind directions of more than 180 degrees and less than 360 degrees as 1, and marking grid points which are not in the range as 0. The altitude range is locked between 5 and 13.5km, and all grid points with wind speed and wind direction labels equal to 1 are extracted, i.e. all west wind rapids.

2. Dongfeng torrent attribute extraction

Step 2-1: meridional section analysis

From 0 ° E to 140 ° E, the longitudinal sections (as shown in fig. 2) in the range of 15 ° S to 30 ° N are traversed one by one with the longitude resolution of the three-dimensional wind direction and wind speed information as a step, and the presence or absence of an east wind torrent is determined for each section. If the latitude range (width) of the east wind torrent is less than 500km, the torrent is rejected, and the east wind torrent attribute label of the section is set to False; if the range of the altitude of the torrent is more than or equal to 500km,the east torrent attribute tag of this profile is set to True. For each meridional section with the east wind torrent attribute label of True, calculating the width and height of the meridional section, the central position of the minimum circumscribed rectangle, the maximum wind speed value and the maximum wind speed value position information according to the minimum circumscribed rectangle (shown in figure 2) of the east wind torrent, and calculating to obtain the horizontal wind shear according to 4 grid points (shown in figure 3) around the maximum wind speed value grid point

And vertical wind shear

Wherein W_s(max)Is the maximum wind speed value, W_s(s)Is the wind speed value, W, of the south grid of the maximum wind speed grid point_s(n)Is the wind speed value, W, of the grid on the north side of the maximum wind speed grid point_s(u)Is the wind speed value, W, of the grid above the maximum wind speed grid point_s(d)Is the wind speed value of the grid below the maximum wind speed grid point.

Step 2-2: obtaining the integral attribute information of the torrent

And recording the start longitude and the end longitude of a continuous uninterrupted longitude sequence with the east wind torrent attribute label of True from 0 DEG E to 140 DEG E, and calculating the length L of the east wind torrent to be 110km multiplied by delta theta according to the difference delta theta of the start longitude and the end longitude. For the width, height, central position, maximum wind speed value position and horizontal wind shear S of all longitudinal sections in the starting and stopping longitude range_hAnd vertical wind shear S_vAnd (4) averaging, and calculating to obtain the width, height, central point position, horizontal and vertical wind shear and other related attribute values of the whole east wind torrent.

Step 2-3: obtaining attribute information of torrent axis

The torrent axis is a long axis of the torrent center, is a connecting line of maximum wind speed points on each longitudinal vertical section in the torrent area, and is provided with one or more maximum value centers of wind speed, namely torrent cores, on the axis of the long and narrow torrent zone. Since wind speed and wind shear near the fast axis are large and have important influence on flight safety, the intensity and position information of the fast axis needs to be extracted. And (3) combining the maximum wind speed values on all the warp-wise sections in the east wind torrent start-stop longitude range obtained in the step (2-2) into an ordered array, and recording the maximum wind speed value (the wind speed value is greater than the wind speed of two adjacent points in the ordered array) and the position information of the maximum wind speed value in the ordered array, namely obtaining the torrent core strength and the position information. The maximum point on the east wind torrent starting radial profile, the position point of the torrent core and the maximum point on the east wind torrent ending radial profile form a torrent axis of the tropical east wind.

Step 2-4: obtaining torrent horizontal distribution information

And (3) extracting the position profiles of different wind speed thresholds on each flight height layer influenced by the east wind torrent according to the height and the central position information of the east wind torrent obtained in the step (2-2), and acquiring the horizontal spatial distribution information of the east wind torrent on different height layers. Specifically, closed position contour sequences with different wind speed levels, such as the threshold values of 30, 35, 40, 45m/s and the like, are calculated by using a contour function in a matplotlib library based on python language. Because the contour function uses real numbers in 0-Nx (Ny) as coordinates to record the position, where Nx and Ny are total number of grids in x and y directions, coordinate conversion is required to be performed according to start and stop longitudes and latitudes on a planar grid, and a longitude and latitude coordinate sequence of the position contour is obtained. The calculation formula is as follows:

wherein (X, Y) is the coordinate of a certain position point calculated by a contour function, (J, W) is the longitude and latitude coordinate of the position point after conversion, and J_sAnd J_eStart-to-end longitude, W, of a planar grid_sAnd W_eThe starting and stopping latitude of the plane grid.

3. Western wind torrent attribute extraction

Step 3-1: classification of western wind torrent

Traversing warp sections in a range from 90 degrees S to 90 degrees N one by taking longitude resolution of three-dimensional wind direction and wind speed information as a step length from 0 degrees to 0 degrees from one week around the earth from west to east, if the latitude range (width) of a certain torrent on a certain warp section is less than 500km, rejecting the torrent on the warp section, then reserving the first 4 torrents with the largest width in the remaining torrents according to the width of each torrent, rejecting the remaining torrents, ensuring that the number of the torrents on each section is equal to or less than 4, and finally classifying each torrent on the warp section by using a west torrent classification matrix. The western torrent classification matrix is shown in table 2, and the numerical values in the table represent the number of grid points in a torrent range on a certain longitudinal section in the latitude distribution range of four types of torrents. The four rushes select the rushes with the largest number of grid points in the range as the western-style rushes, if the largest number of grid points multiplied by the spacing distance between grids (1 latitude equal to about 110km) is less than 1000km, as shown by the southern hemisphere extreme rushes in table 2, the attribute label of such western-style rushes of the profile is denoted as False; if the maximum number of grid points multiplied by the separation distance between the grids is greater than or equal to 1000km, such west wind torrent attribute label of the section is True, for example, in table 2 the north hemisphere polar torrent attribute label on the section is True, torrent 1 is north hemisphere polar torrent, torrent 2 is north hemisphere subtropic west wind torrent on the section, and torrent 3 is south hemisphere subtropic west wind torrent on the section. And (3) carrying out torrent classification on each meridional section from 0 degrees to 0 degrees around the earth from west to east by using the west torrent classification matrix, and calculating the width, the height, the central position of the minimum circumscribed rectangle, the maximum wind speed value and the maximum wind speed value position of the meridional section, and horizontal wind shear and vertical wind shear information according to the minimum circumscribed rectangle of the torrent by using the method of the step 2-1 for the meridional section with the west torrent attribute label of True.

TABLE 2

Step 3-2: obtaining overall attribute information of various western wind torrents

And (2) selecting a profile sequence with continuous True attribute labels for each type of western wind torrent (for each type of western wind torrent, the number of the torrent may not be unique, so that more than one continuous profile sequence is possible), recording start and stop longitudes of each continuous profile sequence, and calculating by using the method in the step 2-2 to obtain relevant attribute values such as length, width, height and the like of each torrent of each type of western wind torrent.

Step 3-3: obtaining attribute information of various western wind torrent shafts

And (4) calculating to obtain the torrent axis attribute information of each type of the western torrent by using the method in the step (2-3) aiming at each type of the western torrent.

Step 3-4: obtaining various western wind torrent horizontal distribution information

At present, aeronautical weather service products at home and abroad only contain three-dimensional warp wind and weft wind data and do not have high-altitude torrent information. The method calculates the high-altitude torrent information such as torrent intensity, torrent axis position, horizontal distribution of torrent layers and the like from the three-dimensional warp wind and weft wind data, provides accurate and visual high-altitude torrent information for each participant of the air transportation, provides a reference basis for planning an air route, and can effectively guarantee the safety and the economy of the air transportation.

The invention provides a high-altitude torrent extraction method based on three-dimensional wind field data, and a plurality of methods and ways for specifically implementing the technical scheme, the above description is only a preferred embodiment of the invention, and it should be noted that, for those skilled in the art, a plurality of improvements and decorations can be made without departing from the principle of the invention, and these improvements and decorations should also be regarded as the protection scope of the invention. All the components not specified in the present embodiment can be realized by the prior art.

Claims

1. A high-altitude torrent extraction method based on three-dimensional wind field data is characterized by comprising the following steps:

step 2: extracting properties of the Dongfeng torrent;

and step 3: extracting the properties of the western wind torrent;

the step 1 comprises the following steps:

step 1-2: extracting the tropical eastern wind torrent;

step 1-3: extracting the western wind torrent;

the step 1-1 comprises the following steps:

Wherein P is₁,z₁T is the pressure, height and temperature of the reference surface at zero degrees Celsius, P₂、z₂Target pressure and height, respectively, the above formula is simplified to Z₂＝18400×log(1013.24/P₂) Converting the constant pressure surface to the nearest flying height according to a simplified formula;

wherein W_sIs the wind speed, W_dThe three-dimensional wind direction and speed information on the flight height layer is obtained through the calculation;

the step 1-2 comprises the following steps:

the method comprises the steps of narrowing an extraction range of tropical east wind torrent to a height range of 9000-16000 m, locking a latitude range between 15 DEG S and 30 DEG N, locking a longitude range between 0 DEG and 140 DEG E, carrying out binary labeling processing on three-dimensional wind direction and wind speed information according to a wind speed threshold of 30m/S, marking grid points with wind speeds of more than or equal to 30m/S as 1, marking grid points with wind speeds of less than 30m/S as 0, carrying out binary labeling processing on the three-dimensional wind direction and wind speed information according to the wind directions, marking grid points with wind directions of more than 0 DEG and less than 180 DEG as 1, and marking grid points with wind directions of not in the range as 0; extracting all grid points with wind speed and wind direction labels equal to 1 according to the space range of the east wind torrent, namely the tropical east wind torrent;

the steps 1-3 comprise:

according to a wind speed threshold value of 30m/s, carrying out binary labeling processing on the three-dimensional wind direction and wind speed information, marking grid points with wind speeds of more than or equal to 30m/s as 1 and grid points with wind speeds of less than 30m/s as 0, then carrying out binary labeling processing on the three-dimensional wind direction and wind speed information according to the wind directions, marking grid points with wind directions of more than 180 degrees and less than 360 degrees as 1, and marking grid points which are not in the range as 0; locking the height range between 5 and 13.5km, and extracting all grid points with wind speed and wind direction labels equal to 1, namely all western wind torrent;

the step 2 comprises the following steps:

And vertical wind shear

wherein (X, Y) is a position point coordinate calculated by a contour function, (J, W) is a longitude and latitude coordinate of the position point after conversion, and J_sAnd J_eRespectively start and end longitude, W, of the plane grid_sAnd W_eRespectively the starting latitude and the stopping latitude of the plane grid;

in the step 2-3, if the wind speed value of one point in the ordered array is larger than the wind speed values of two adjacent points, the wind speed value of the point is the maximum wind speed value;

the step 3 comprises the following steps:

step 3-1: classifying the western wind torrent; starting from 0 degrees, going around the earth from west to east, going to 0 degrees, traversing warp sections in a range from 90 degrees S to 90 degrees N one by taking longitude resolution of three-dimensional wind direction and wind speed information as a step length, if the latitude range of one torrent on one warp section is less than 500km, eliminating the torrent on the warp section, then reserving the first 4 torrents with the largest width in the remaining torrents according to the width of each torrent, eliminating the remaining torrents, ensuring that the number of the torrents on each section is equal to or less than 4, and finally classifying each torrent on the warp section by using a west torrent classification matrix:

the contents of a first row and a first column of a western wind torrent classification matrix of a meridional section are used as titles, the first row lists 4 types of western wind torrents, the first column lists torrent sequences in the meridional section, numerical values in a western wind torrent classification matrix table represent the number of grid points in a torrent range on the meridional section, which fall in a latitude distribution range of four types of torrents, the four types of torrents select the torrent with the largest number of grid points in the range as the type of western wind torrent, and if the number of the largest grid points is multiplied by the spacing distance between grids and is less than 1000km, the attribute label of the type of western wind torrent of the section is marked as a False value False; if the maximum grid point number multiplied by the spacing distance between the grids is greater than or equal to 1000km, the attribute label of the western wind torrent of the section is marked as True value True; for each warp-wise section from 0 degrees to 0 degrees around the earth from west to east, performing torrent classification by using the west torrent classification matrix, and for the warp-wise sections with True-value attribute labels of each category of west torrent, calculating the width and height of the warp-wise section, the central position of the minimum circumscribed rectangle, the maximum wind speed value and the maximum wind speed value position, and horizontal wind shear and vertical wind shear information according to the minimum circumscribed rectangle of the torrent by using the method in the step 2-1;

step 3-3: acquiring attribute information of various western wind torrent axes: