1. A typhoon vortex initial wind field optimization method is characterized by comprising the following steps:

respectively acquiring a streamline direction vector matrix sequence and a wind direction vector matrix sequence of a vortex region in a continuous period T of typhoon to be optimized;

and calculating a cross product value of the streamline direction vector matrix sequence and the wind direction vector matrix sequence, judging whether the cross product value is gradually reduced along with the time within the continuous period of time T, and if the cross product value is not gradually reduced along with the time within the continuous period of time T, adjusting the wind speed value within the continuous period of time T.

2. The method for optimizing the initial wind field of the typhoon vortex according to claim 1, wherein the method for obtaining the streamline direction vector matrix sequence of the vortex region in the typhoon to be optimized within a continuous period of time T comprises the following steps:

obtaining the vortex remote sensing images at each moment in the continuous period of time T, constructing a streamline direction vector matrix of a vortex region in the vortex remote sensing images at each moment, and constructing a streamline direction vector matrix sequence { D) in the continuous period of time T based on the streamline direction vector matrix of the vortex region in the vortex remote sensing images at each moment_iIn which D is_iIs any vortex remote sensing image p at any ith moment in the continuous period of time T_iThe streamline direction vector matrix of the inner vortex area, i ═ 1,2,3.

3. The method of optimizing an initial wind field of a typhoon vortex according to claim 2, wherein the streamline direction vector matrix D_iObtained by the following steps:

from the vortex remote sensing image p_iMiddle extracted vortex area m_iFor the vortex region m_iPerforming grid division, and respectively extracting the vortex areas m by taking the brightness temperature values as variables_iA flow field for each grid within;

judging whether each grid is on the corresponding streamline field, if so, assigning the streamline vector of the grid on the corresponding streamline field to be the direction of the corresponding streamline field, and if not, assigning the streamline vector of the grid not on the corresponding streamline field to be 0;

based on the vortex region m_iThe flow line field direction of each grid in the vortex region m is constructed_iStreamline direction vector matrix D of_i。

4. The method for optimizing the initial wind field of the typhoon vortex, according to claim 1, is characterized in that the method for obtaining the wind direction vector matrix sequence of the vortex area in the typhoon to be optimized within a continuous period of time T comprises the following steps:

acquiring vortex remote sensing images and corresponding initial wind field data at each moment in the continuous period of time T;

respectively constructing wind direction vector matrixes of vortex areas in the vortex remote sensing images at each moment based on initial wind field data corresponding to the vortex remote sensing images at each moment, and constructing a wind direction vector matrix sequence { W) within a continuous period of time T based on the wind direction vector matrixes of the vortex areas in the vortex remote sensing images at each moment_i}，W_iIs any vortex remote sensing image p at any ith moment in the continuous period of time T_iI is 1,2,3, L is the total number of time instants of a continuous period of time T.

5. The method for optimizing the initial wind field of the typhoon vortex according to claim 4, wherein any vortex remote sensing image p at any ith moment in the continuous period T_iWind direction vector matrix W of vortex region_iThe method comprises the following steps:

from the vortex remote sensing image p_iMiddle extracted vortex area m_iFor the vortex region m_iCarrying out grid division;

calculating a wind direction vector corresponding to each grid based on the initial wind field data corresponding to each grid;

based on the vortex region m_iWind direction vector of each grid in the vortex region m_iWind direction vector matrix W_i。

6. The method for optimizing the initial wind field of a typhoon vortex according to claim 1, wherein the wind speed value of the continuous period of time T is adjusted, and the method comprises the following steps:

and selecting all wind speed values corresponding to the moments when the slope of the cross product value is greater than 0 from the continuous period of time T as alternative adjustment wind speed values, sequentially selecting the alternative adjustment wind speed values from morning to evening according to the time sequence for one-by-one adjustment, updating the adjusted cross product value of the continuous period of time T after each adjustment, and stopping the adjustment of the wind speed value at the subsequent moment if the updated cross product value becomes smaller gradually along with the time in the continuous period of time T.

7. The method for optimizing the initial wind field of the typhoon vortex according to the claim 6, wherein the adjustment of the wind speed value at any j time comprises the following steps:

streamline direction vector matrix D for j time_jSeparately calculating the streamline direction vector matrix D_jThe direction difference between the streamline direction vector and the wind direction vector of the grid with each streamline direction vector not being 0;

and sorting the direction difference of the grids with the streamline direction vectors not being 0 from large to small, and selecting grid wind speed values with the preset proportion and the front sorting to adjust.

8. The method for optimizing an initial wind field of a typhoon vortex according to claim 7, wherein the following steps are performed for any grid of grid wind speed values to be adjusted:

respectively calculating a wind speed value u (k, f) in the north-south direction of the grid and a wind speed value v (k, f) in the east-west direction of the grid, and adjusting the wind speed value u (k, f) in the north-south direction of the grid and the wind speed value v (k, f) in the east-west direction of the grid according to the wind direction of the grid until the sum of the wind speed adjustment variances between the wind speed value u (k, f) in the north-south direction and the wind speed value v (k, f) in the east-west direction of the grid and the surrounding grid is minimum;

wherein, the Luv calculation formula is as follows:

luv is a sum of wind speed adjustment variances, α is a first coordinate difference between the grid and a surrounding grid, β is a second coordinate difference between the grid and the surrounding grid, u (k, f) is a wind speed value in the north-south direction of the grid, u (k + α, f + β) is a wind speed value in the north-south direction of the grid surrounding the grid, v (k, f) is a wind speed value in the east-west direction of the grid, and v (k + α, f + β) is a wind speed value in the east-west direction of the grid surrounding the grid.

9. A computer system comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any of the preceding claims 1 to 8 are carried out by the processor when the computer program is executed by the processor.

10. A computer storage medium having a computer program stored thereon, wherein the program is adapted to perform the steps of the method of any one of claims 1 to 8 when executed by a processor.

Background

Landing typhoons are destructive to power grid striking, so that hardware fittings are damaged, towers are broken, large-area power failure is caused, and great harm is caused. The typhoon position and the typhoon strength can be accurately predicted, and the method has important value for preventing typhoon of the power grid.

Typhoon prediction is carried out, and firstly, the initial position and the strength of the typhoon are determined. Since typhoons tend to be generated on the ocean, observation data is scarce. At present, the positioning calculation of typhoon mainly depends on satellite monitoring data, but the satellite monitoring usually has better monitoring capability on vortex with obvious typhoon characteristics, but the positioning calculation deviation of the vortex with weak typhoon characteristics such as unobvious typhoon eyes is large, and the typhoon is directly predicted by directly adopting a numerical simulation method, and the later-period prediction effect is poor due to lack of observation data and large calculation deviation.

Therefore, the low accuracy of the existing numerical simulation typhoon prediction method has become a technical problem to be solved urgently by those skilled in the art.

Disclosure of Invention

The invention provides a typhoon vortex initial wind field optimization method, a typhoon vortex initial wind field optimization system and a storage medium, which are used for solving the technical problem that the accuracy of the typhoon prediction method of numerical simulation is low.

In order to solve the technical problems, the technical scheme provided by the invention is as follows:

a typhoon vortex initial wind field optimization method comprises the following steps:

respectively acquiring a streamline direction vector matrix sequence and a wind direction vector matrix sequence of a vortex region in a continuous period T of typhoon to be optimized;

and calculating the cross product value of the streamline direction vector matrix sequence and the wind direction vector matrix sequence, judging whether the cross product value is gradually reduced along with the time within a continuous period of time T, and if the cross product value is not gradually reduced along with the time within the continuous period of time T, adjusting the wind speed value within the continuous period of time T.

Preferably, the method for acquiring the streamline direction vector matrix sequence of the vortex region in the typhoon to be optimized within a continuous period of time T comprises the following steps:

obtaining vortex remote sensing images at each moment in a continuous period of time T, constructing a streamline direction vector matrix of a vortex region in the vortex remote sensing images at each moment, and constructing a streamline direction vector matrix sequence { D) in the continuous period of time T based on the streamline direction vector matrix of the vortex region in the vortex remote sensing images at each moment_iIn which D is_iIs any vortex remote sensing image p at any ith moment in a continuous period of time T_iThe streamline direction vector matrix of the inner vortex area, i ═ 1,2,3.

Preferably, the streamline direction vector matrix D_iObtained by the following steps:

remote sensing of images p from vortices_iMiddle extracted vortex area m_iFor the vortex region m_iPerforming grid division, and respectively extracting vortex areas m by taking the brightness temperature value as a variable_iA flow field for each grid within;

judging whether each grid is on the corresponding streamline field, if so, assigning the streamline vector of the grid on the corresponding streamline field to be the direction of the corresponding streamline field, and if not, assigning the streamline vector of the grid not on the corresponding streamline field to be 0;

based on the vortex region m_iFlow field direction of each grid in the vortex region m_iStreamline direction vector matrix D of_i。

Preferably, the method for acquiring the wind direction vector matrix sequence of the vortex region in the typhoon to be optimized within a continuous period of time T comprises the following steps:

acquiring vortex remote sensing images and corresponding initial wind field data at each moment in a continuous period of time T;

respectively constructing a wind direction vector matrix of an inner vortex region of the vortex remote sensing image at each moment based on initial wind field data corresponding to the vortex remote sensing image at each moment, and constructing a wind direction vector matrix sequence { W) within a continuous period of time T based on the wind direction vector matrix of the inner vortex region of the vortex remote sensing image at each moment_i}，W_iIs any vortex remote sensing image p at any ith moment in a continuous period of time T_iI is 1,2,3, L is the total number of time instants of a continuous period of time T.

Preferably, any vortex remote sensing image p at any ith moment in continuous period T_iWind direction vector matrix W of vortex region_iThe method comprises the following steps:

remote sensing of images p from vortices_iMiddle extracted vortex area m_iFor the vortex region m_iCarrying out grid division;

calculating a wind direction vector corresponding to each grid based on the initial wind field data corresponding to each grid;

based on the vortex region m_iWind direction vector of each grid in the vortex region m_iWind direction vector matrix W_i。

Preferably, the wind speed value for a continuous period of time T is adjusted, including the following steps:

selecting all wind speed values corresponding to the moments when the slope of the cross product value is greater than 0 from a continuous period of time T as alternative adjustment wind speed values, sequentially selecting the alternative adjustment wind speed values from morning to evening according to the time sequence for one-by-one adjustment, updating the cross product value of the adjusted continuous period of time T after each adjustment, and stopping the adjustment of the wind speed value at the subsequent moment if the updated cross product value is gradually reduced along with the time in the continuous period of time T.

Preferably, the adjustment of the wind speed value at any j time comprises the following steps:

streamline direction vector matrix D for j time_jSeparately calculating a streamline direction vector matrix D_jThe direction difference between the streamline direction vector and the wind direction vector of the grid with each streamline direction vector not being 0;

and sorting the direction difference of the grids with the streamline direction vectors not being 0 from large to small, and selecting grid wind speed values with the preset proportion and the front sorting to adjust.

Preferably, for any grid of grid wind speed values to be adjusted, the following steps are performed:

respectively calculating a wind speed value u (k, f) in the north-south direction of the grid and a wind speed value v (k, f) in the east-west direction of the grid, and adjusting the wind speed value u (k, f) in the north-south direction of the grid and the wind speed value v (k, f) in the east-west direction of the grid according to the wind direction of the grid until the sum of the wind speed adjustment variances between the wind speed value u (k, f) in the north-south direction and the wind speed value v (k, f) in the east-west direction of the grid and the surrounding grid is minimum;

wherein, the Luv calculation formula is as follows:

luv is a sum of wind speed adjustment variances, α is a first coordinate difference between the grid and a surrounding grid, β is a second coordinate difference between the grid and the surrounding grid, u (k, f) is a wind speed value in the north-south direction of the grid, u (k + α, f + β) is a wind speed value in the north-south direction of the grid surrounding the grid, v (k, f) is a wind speed value in the east-west direction of the grid, and v (k + α, f + β) is a wind speed value in the east-west direction of the grid surrounding the grid.

A computer system comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the steps of the method being performed when the computer program is executed by the processor.

A computer storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the above method.

The invention has the following beneficial effects:

1. according to the typhoon vortex initial wind field optimization method, system and storage medium, a streamline direction vector matrix sequence and a wind direction vector matrix sequence of a vortex region in a continuous period T of typhoon to be optimized are obtained respectively; calculating a cross product value of the streamline direction vector matrix sequence and the wind direction vector matrix sequence, judging whether the cross product value is gradually reduced along with the time within a continuous period of time T, and if the cross product value is not gradually reduced along with the time within the continuous period of time T, adjusting the wind speed value within the continuous period of time T; compared with the prior art, the typhoon prediction method based on numerical simulation is constrained by extracting typhoon spiral characteristics (namely the streamline direction vector matrix sequence of the vortex region) of the typhoon, so that the optimized calculation of the typhoon initial field is realized, the uncertainty of typhoon prediction is reduced, and the typhoon prediction precision is improved.

2. In the optimal selection method, the brightness temperature data is skillfully adopted by the streamline calculation method, the spiral characteristic of the typhoon initial vortex is extracted, the spiral characteristic is taken as the constraint, the initial wind field of the typhoon numerical mode is corrected, and the method is convenient and easy to operate.

3. The method has good universality, and can be used for different regions of the world where typhoons are easy to occur to initialize and calculate the typhoons.

4. By adopting the technology of the invention, the automatic adjustment calculation of the initial vortex of the typhoon wind field can be realized, the calculation speed is high, the correction can be carried out for many times, and the typhoon prediction accuracy is further improved.

In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The present invention will be described in further detail below with reference to the accompanying drawings.

Drawings

The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:

FIG. 1 is a flow chart of a method for optimizing an initial wind field of a typhoon vortex in a preferred embodiment of the invention;

FIG. 2 is a schematic view of the center of the typhoon and the helix vector (W) in the preferred embodiment of the present invention (the length of the arrow is meaningless, and the direction of the arrow represents the helix direction);

fig. 3 is a schematic wind direction diagram of the typhoon center and helix versus position in the numerical pattern in the preferred embodiment of the invention (the length of the arrow indicates the magnitude of the full wind speed, and the direction of the arrow indicates the wind direction of the numerical pattern versus position of the helix).

Detailed Description

The embodiments of the invention will be described in detail below with reference to the drawings, but the invention can be implemented in many different ways as defined and covered by the claims.

The first embodiment is as follows:

as shown in fig. 1, this implementation discloses a method for optimizing an initial wind field of a typhoon vortex, which includes the following steps:

respectively acquiring a streamline direction vector matrix sequence and a wind direction vector matrix sequence of a vortex region in a continuous period T of typhoon to be optimized;

and calculating the cross product value of the streamline direction vector matrix sequence and the wind direction vector matrix sequence, judging whether the cross product value is gradually reduced along with the time within a continuous period of time T, and if the cross product value is not gradually reduced along with the time within the continuous period of time T, adjusting the wind speed value within the continuous period of time T.

In addition, in the embodiment, a computer system is also disclosed, which includes a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the computer program, the steps of the method are implemented.

Furthermore, in the present embodiment, a computer storage medium is also disclosed, on which a computer program is stored, which when executed by a processor implements the steps in the above method.

According to the typhoon vortex initial wind field optimization method, system and storage medium, a streamline direction vector matrix sequence and a wind direction vector matrix sequence of a vortex region in a continuous period T of typhoon to be optimized are obtained respectively; calculating a cross product value of the streamline direction vector matrix sequence and the wind direction vector matrix sequence, judging whether the cross product value is gradually reduced along with the time within a continuous period of time T, and if the cross product value is not gradually reduced along with the time within the continuous period of time T, adjusting the wind speed value within the continuous period of time T; compared with the prior art, the typhoon prediction method based on numerical simulation is constrained by extracting typhoon spiral characteristics (namely the streamline direction vector matrix sequence of the vortex region) of the typhoon, so that the optimized calculation of the typhoon initial field is realized, the uncertainty of typhoon prediction is reduced, and the typhoon prediction precision is improved.

Example two:

the second embodiment is the preferred embodiment of the first embodiment, and the difference between the first embodiment and the second embodiment is that the specific steps of the typhoon vortex initial wind field optimization method are refined:

in the embodiment, a typhoon vortex initial wind field optimization method based on vortex reconstruction is disclosed, the method extracts typhoon spiral features according to monitoring data of synchronous satellites at adjacent moments, uses the spiral features as constraints, optimizes the typhoon initial vortex wind field, puts the optimized and corrected initial field into a numerical mode again to start predictive computation, and comprises the following specific steps:

(1) synchronous satellite typhoon monitoring data screening

Setting a monitoring range of initial typhoon according to the monitoring range of the synchronous satellite, judging whether vortex motion exists in the monitoring range of the typhoon, and continuously collecting vortex satellite remote sensing images with the time length of T when vortex motion is found, wherein T is more than or equal to 3 hours.

(2) Calculation of typhoon spiral characteristics of geostationary satellite

Selecting a vortex area M according to the synchronous satellite vortex remote sensing image collected in the step (1), and dividing the vortex area M into M multiplied by n grids, wherein the grid division is consistent with the adopted typhoon numerical mode initial field. And each grid takes the brightness temperature value as a variable and adopts ArcGis software to extract the streamline of the area M. And (5) repeating the step (2) to obtain the flow line field H of each satellite remote sensing image.

(3) Streamline direction vector matrix calculation

And (3) judging whether each grid is on the streamline in the area M according to the streamline field obtained in the step (2). When the grid is on the flow line, the flow line direction is taken as the value of the point, and when the grid is not on the flow line, the value of the point is 0. Traversing all grid points in the region M to obtain a streamline direction vector matrix D_iObtaining the streamline direction vector matrix at all the moments to obtain a vector matrix sequence (D)_i}。

Wherein, the concrete steps of the steps (2) and (3) are as follows:

(2.1) constructing a streamline direction vector matrix of a vortex area in the vortex remote sensing image at each moment;

wherein, any vortex remote sensing image p at any ith moment in continuous period T_iStreamline direction vector matrix D of inner vortex area_iObtained by the following steps:

remote sensing of images p from vortices_iMiddle extracted vortex area m_iFor the vortex region m_iPerforming grid division, and respectively extracting vortex areas m by taking the brightness temperature value as a variable_iA flow line field of each grid, wherein i is 1,2,3, and L is a total number of time instants of a continuous period of time T;

judging whether each grid is on the corresponding streamline field, if so, assigning the streamline vector of the grid on the corresponding streamline field to be the direction of the corresponding streamline field, and if not, assigning the streamline vector of the grid not on the corresponding streamline field to be 0;

based on the vortex region m_iFlow field direction of each grid in the vortex region m_iStreamline direction vector matrix D of_i。

(2.2) constructing the streamline direction in a continuous period of time T based on the streamline direction vector matrix of the vortex region in the vortex remote sensing image at each momentVector matrix sequence { D_iIn which D is_iIs any vortex remote sensing image p at any ith moment in a continuous period of time T_iThe streamline direction vector matrix of the inner vortex area, i ═ 1,2,3.

(4) Typhoon numerical mode wind field wind direction calculation

Collecting typhoon numerical value mode initial wind field data at corresponding time according to the satellite observation time determined in the step (1), calculating a wind direction vector in each grid, and obtaining an initial field wind direction vector matrix W_iObtaining the initial field wind direction vector matrix W at all the time_iObtaining a field wind direction vector matrix sequence { W_i}。

Wherein, the step (4) comprises the following steps:

(4.1) acquiring initial wind field data corresponding to the vortex remote sensing images at each moment in a continuous period of time T;

(4.2) respectively constructing a wind direction vector matrix of the vortex area in the vortex remote sensing image at each moment based on the initial wind field data corresponding to the vortex remote sensing image at each moment:

calculating a wind direction vector corresponding to each grid based on the initial wind field data corresponding to each grid;

based on the vortex region m_iWind direction vector of each grid in the vortex region m_iWind direction vector matrix W_i。

(4.3) constructing a wind direction vector matrix sequence { W) in a continuous period of time T based on the wind direction vector matrix of the vortex region in the vortex remote sensing image at each moment_i}，W_iIs any vortex remote sensing image p at any ith moment in a continuous period of time T_iI is 1,2,3, L is the total number of time instants of a continuous period of time T.

(5) Vortex initial field wind field constraint judgment

Obtaining a streamline vector matrix D at all moments according to the step (3)_iAnd the numerical mode field vector matrix W obtained in the step (4)_iCalculating the cross product of the two values, and recording as S_iAnd for all S_iAnd sorting is carried out.

S_i＝D_i×W_i

When S is_iThe numerical mode is gradually close to the real observation along with the gradual reduction of the time, the adjustment is not needed, and the typhoon numerical mode calculation can be directly carried out.

When S is_iGradually adjusting each moment one by one along with the gradual non-reduction of the time, wherein the step specifically comprises the following steps:

selecting all wind speed values corresponding to the moments with the cross product value slopes larger than 0 (namely the cross product value at the moment is larger than the cross product value at the previous moment) from a continuous period of time T as alternative adjustment wind speed values, sequentially selecting the alternative adjustment wind speed values from the morning to the evening according to the time sequence for carrying out adjustment one by one, updating the adjusted cross product value of the continuous period of time T after each adjustment, and stopping the adjustment of the wind speed values at the subsequent moments if the updated cross product value is gradually reduced along with the time in the continuous period of time T.

(6) Initial wind field calculation correction

When the step (5) judges that the adjustment is needed, executing the step, wherein the specific adjustment process is as follows:

comparing streamline vector matrix D of time j to be adjusted_jAnd a numerical mode field vector matrix W_jTo D, pair_jThe grid point (existence direction vector) d (k, f) of which the middle is not 0 is calculated_jThe direction difference Δ p of the direction vector w (k, f) of the middle corresponding position;

sorting the delta p, and selecting the first 30% with larger value for adjustment.

Selecting a numerical value to be adjusted to calculate the wind speed u (k, f) in the north-south direction and the wind speed v (k, f) in the east-west direction in the grid, adjusting two wind speed values according to the wind direction of the corresponding streamline grid, wherein the adjustment process meets the minimum of variance and Luv; wherein, the Luv calculation formula is as follows:

and (5) repeating the step (6) to complete the adjustment of all grids at the current moment j, and repeating the step (5) to judge whether the adjustment is needed at the next moment.

(7) Typhoon numerical prediction calculation

And (5) after the adjustment of the step (5) and the step (6) is completed, carrying out typhoon numerical calculation to obtain the typhoon forecast at the future moment.

Example three:

the third embodiment is a preferred embodiment of the second embodiment, and is different from the second embodiment in that the third embodiment is combined with a specific typhoon scene, and comprises the following steps:

(1) synchronous satellite typhoon monitoring data screening

Setting the monitoring range of the initial typhoon to be 0-10 degrees N and 120-130 degrees E according to the monitoring range of the synchronous satellite, judging that vortex motion exists in the monitoring range of the typhoon, and continuously collecting remote sensing images of the sunflower 8 satellite for 5 hours.

(2) Calculation of typhoon spiral characteristics of geostationary satellite

Selecting a vortex area M according to the synchronous satellite vortex remote sensing image collected in the step (1), and dividing the vortex area M into 10 multiplied by 10 grids, wherein the grid division is consistent with the adopted typhoon numerical mode initial field. And each grid takes the brightness temperature value as a variable and adopts ArcGis software to extract the streamline of the area M. And (5) repeating the step (2) to obtain the flow line field H of each satellite remote sensing image.

(3) Streamline direction vector matrix calculation

And (3) judging whether each grid is on the streamline in the area M according to the streamline field obtained in the step (2). When the grid is on the flow line, the flow line direction is taken as the value of the point, and when the grid is not on the flow line, the value of the point is 0. Traversing all grid points in the region M to obtain a streamline direction vector matrix D shown in figure 2, and obtaining a vector matrix sequence { D ] at all moments_i}。

(4) Typhoon numerical mode wind field wind direction calculation

Collecting typhoon numerical value mode initial time of corresponding time according to satellite observation time determined in step (1)Starting wind field data, calculating wind direction vectors in each grid to obtain an initial field wind direction vector matrix W shown in figure 3, and obtaining a W vector matrix sequence { W_i}。

(5) Vortex initial field wind field constraint calculation

Obtaining a streamline vector matrix D at all moments according to the step (3)_iAnd the numerical mode field vector matrix W obtained in the step (4)_iCalculating the cross product of the two values, and recording as S_iAnd for all S_iAnd sorting is carried out.

S_i＝D_i×W_i

When S is_iThe numerical mode is gradually close to the real observation along with the gradual reduction of the time, the adjustment is not needed, and the typhoon numerical mode calculation can be directly carried out.

When S is_iAnd gradually adjusting each moment one by one without reducing the time.

(6) Initial field assimilation correction

When the step (5) judges that the adjustment is needed, executing the step, wherein the specific adjustment process is as follows:

comparing streamline vector matrix D of time j to be adjusted_jAnd a numerical mode field vector matrix W_jTo D, pair_jThe grid point (existence direction vector) d (k, f) of which the middle is not 0 is calculated_jDirection difference Δ p of direction vector w (k, f) of middle corresponding position

Sorting the deltap, and selecting the first 30% with larger value for adjustment.

Selecting a numerical value to be adjusted to calculate the wind speed u (k, f) in the north-south direction and the wind speed v (k, f) in the east-west direction in the grid, adjusting two wind speed values according to the wind direction of the corresponding streamline grid, wherein the adjustment process meets the minimum of variance and Luv

Luv is the sum of wind speed adjustment variances, alpha is a first coordinate difference value between the grid and a surrounding grid, beta is a second coordinate difference value between the grid and the surrounding grid, u (k, f) is a wind speed value in the north-south direction of the grid, u (k + alpha, f + beta) is a wind speed value in the north-south direction of the grid around the grid, v (k, f) is a wind speed value in the east-west direction of the grid, and v (k + alpha, f + beta) is a wind speed value in the east-west direction of the grid around the grid

And (5) repeating the step (6) to complete the adjustment of all grids at the current moment j, and repeating the step (5) to judge whether the adjustment is needed at the next moment.

(7) Typhoon number calculation

And (5) after the adjustment of the step (5) and the step (6) is completed, carrying out typhoon numerical calculation to obtain the typhoon forecast at the future moment.

In summary, according to the typhoon vortex initial wind field optimization method, the typhoon vortex initial wind field optimization system and the storage medium, a streamline direction vector matrix sequence and a wind direction vector matrix sequence of a vortex region in a continuous period of time T of typhoon to be optimized are respectively obtained; calculating a cross product value of the streamline direction vector matrix sequence and the wind direction vector matrix sequence, judging whether the cross product value is gradually reduced along with the time within a continuous period of time T, and if the cross product value is not gradually reduced along with the time within the continuous period of time T, adjusting the wind speed value within the continuous period of time T; compared with the prior art, the typhoon prediction method based on numerical simulation is constrained by extracting typhoon spiral characteristics (namely the streamline direction vector matrix sequence of the vortex region) of the typhoon, so that the optimized calculation of the typhoon initial field is realized, the uncertainty of typhoon prediction is reduced, and the typhoon prediction precision is improved.

The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.