Asymmetric parameterized wind field early warning method and system

Abstract

The present invention relates to the technical field of typhoon prevention, and in particular to an asymmetric parameterized wind field early warning method and system, comprising the following processes: reading an input file; initializing simulated typhoon wind field information and grid distribution information according to the input file information for use in subsequent cyclic simulation; performing cyclic calculations according to the obtained variable information; extracting and temporarily storing the calculation results; and drawing according to the calculation results. The method has the advantages of being able to specifically simulate the evolution process of a typhoon and the accompanying time series changes of air pressure, wind speed, wind direction, etc., and thus being able to be used to effectively assess the scope and impact of disasters such as storm surges and typhoon waves in coastal areas of my country, thereby rationally allocating disaster relief resources, optimizing the design standards of coastal defense projects, and accurately formulating disaster prevention and mitigation strategies, and providing guarantees for the economic development of coastal cities.

Description

Technical Field

[0001] The present invention relates to the technical field of typhoon prevention, and in particular to an asymmetric parameterized wind field early warning method and system. Background Art

[0002] The Northwest Pacific warm pool breeds nearly 30% of the world's tropical cyclones (TCs), making it the most active region for TC activity in the world, and the intensity of typhoons ranks first in the world. In recent years, typhoons have gradually increased in intensity, their frequency has continued to rise, and the destructive power of typhoon storm surge disasters has increased. Coastal areas are high-risk areas susceptible to typhoons, with relatively high levels of economic development and concentrated populations. Nearly 40% of the world's population is concentrated within 100 kilometers of the coast. Therefore, typhoon disasters often cause huge economic losses and casualties, and early warning and forecasting of typhoon events are particularly important.

[0003] The sea surface wind field during a typhoon is the most important factor in accurately simulating typhoon waves. Atmospheric dynamics numerical models, such as the mesoscale climate model WRF and the global climate model MPAS, are one of the ways to study typhoon wind fields and play an important role in atmospheric research and operational forecasts and hindcasts. This type of model is designed to solve atmospheric dynamics and thermodynamic equations based on specific boundary conditions and background wind fields, taking into account parameterized physical schemes such as radiation, surface, planetary boundary conditions, and cumulus, microphysics, turbulence, and diffusion, and output atmospheric variables such as wind, pressure, precipitation, and temperature. However, the dynamic numerical model has a huge amount of calculations, so the timeliness of the simulation is limited, which is not conducive to rapid forecasting and probabilistic research to simulate a large number of typhoon events. In addition, the simulation performance of this type of atmospheric model is easily affected by many factors such as physical schemes, initial background fields, and boundary conditions, which poses a great challenge to accurately analyzing the intensity and path of typhoons, especially typhoons with asymmetric structures. Therefore, it is necessary to propose a method that supports the setting of typhoon structure and is efficient and convenient in calculation to solve the above problems.

[0004] As a model that combines pure experience, pure theory, and semi-empirical and semi-theoretical methods, the parameterized model can take into account the physical parameter schemes during the typhoon period and set different typhoon center structures according to the characteristics of the typhoon to simulate the typhoon wind field. At the same time, it can also be combined with wind field data such as FNL (the final operational global analysis data set provided by the National Center for Environmental Prediction (NCEP) of the United States), ERA5 (the latest generation of reanalysis data set provided by the European Center for Medium-Range Weather Forecasts (ECMWF)), CCMP (a multi-platform cross-calibrated sea surface wind field data set developed by the U.S. Naval Research Laboratory), and JRA-55 (a 55-year reanalysis data set jointly produced by the Japan Meteorological Agency (JMA) and the Japan Agency for Marine-Earth Science and Technology (JAMSTEC)) to build a new synthetic typhoon wind field, further improving the accuracy of wind field simulation during the formation, development, and extinction of typhoons. In addition, since the parameterized model does not involve complex thermal and physical processes such as heat exchange, water vapor transport, and turbulence in the atmosphere during typhoon movement, the calculation is more convenient. It can support the simulation of probabilistic massive typhoon events while ensuring computational efficiency, and accurately characterize the main characteristics of the typhoon wind field, which is of great significance to engineering applications and disaster prevention and mitigation.

[0005] Based on this, this case is brought forward. Summary of the invention

[0006] One of the purposes of the present invention is to provide an asymmetric parameterized wind field early warning method, which can specifically simulate the evolution process of a typhoon and the accompanying time series changes in air pressure, wind speed, wind direction, etc., and can then be used to effectively assess the scope and impact of disasters such as storm surges and typhoon waves in my country's coastal areas.

[0007] To achieve the above object, the technical solution of the present invention is as follows:

[0008] An asymmetric parameterized wind field early warning method, comprising:

[0009] S1. Read the input file, including the following steps:

[0010] S11. Read the start time, end time, simulation time step, simulation year, read the number of slice data, and initialize the slice data time, typhoon center longitude and latitude, typhoon range determination parameters, maximum sustained wind speed, typhoon center pressure, maximum wind speed radius, and shape parameters according to the number of slice data;

[0011] S12. Read the grid longitude and latitude interval and grid longitude and latitude range;

[0012] S13. Reading angle and distance conversion method parameters, background pressure, atmospheric density, and wind speed adjustment parameters;

[0013] S14. Read output point information;

[0014] S15. Read wind speed correction parameters and output file name;

[0015] S2. Initialize the simulated typhoon wind field information and grid distribution information according to the input file information obtained in step S1 for use in subsequent cyclic simulations, including the following steps:

[0016] S21. Initialize typhoon time difference sequence information and calculate the distance and direction of adjacent typhoon centers;

[0017] S22. Calculate the moving speed of the typhoon center;

[0018] S23. Initialize the grid distribution and variable information, and calculate the number of cycles;

[0019] S3. Perform loop calculation based on the variable information obtained in steps S1 and S2;

[0020] S4. extract and temporarily store the calculation results;

[0021] S5. Draw according to the calculation results.

[0022] Furthermore, the calculation method of step S21 is as follows:

[0023] First, according to all the typhoon slice data obtained in step S1, the time difference between each slice data time and the start time is calculated in seconds, thereby obtaining the time difference sequence information of the typhoon slice;

[0024] Calculate the distance and direction of the typhoon center of adjacent slice data. The distance calculation formula is as follows:

[0025]

[0026] in Indicates the center of the typhoon To the center of the typhoon , φ1 and φ2 represent the longitudes of the two points. represents latitude, R represents the radius of the earth in nmi, arctan(r,z) represents the angle of the right triangle with adjacent side length z and opposite side length r, where

[0027]

[0028]

[0029] The direction calculation formula is as follows:

[0030]

[0031] in Indicates the center of the typhoon To the center of the typhoon The azimuth of the direction is calculated in clockwise direction with due north as 0°. arctan(y,x) represents the angle of the right triangle with adjacent side length x and opposite side length y. rad2dg=0.017453293 represents the radian corresponding to each angle, which is used to convert angles to radians.

[0032] Furthermore, according to the distances and time differences between adjacent typhoon centers calculated in step S21, the moving speed of the typhoon center can be obtained, and the formula is as follows:

[0033]

[0034] in Indicates the center of the typhoon To the center of the typhoon The moving speed, Indicates that the typhoon center is at point The time difference between the time point and the initial time point, Indicates that the typhoon center is at point The time difference between the time point and the initial time point, so Indicates the typhoon center from point Move to time.

[0035] Further, step S23 initializes grid information and related variables according to the simulation information read in step S1, and calculates the number of cycles, wherein the grid information includes the number of grid longitude and latitude grid points and the longitude and latitude division vectors, and the related variables include east-west wind speed, north-south wind speed, air pressure, latitude vector, and longitude vector;

[0036] The number of cycles is calculated as follows:

[0037]

[0038] Where num represents the number of loops. Indicates the last typhoon center point The time difference between the time point and the initial time point, dt represents the time interval of calculation.

[0039] Furthermore, the step S3 includes step S31: performing linear interpolation on some variables according to the current cycle time point, and the some variables include: the longitude and latitude of the typhoon center at the current time, the maximum sustained wind speed, the central air pressure, the central air pressure change rate, the moving speed of the typhoon center at the current moment, the moving direction, the maximum wind speed radius, the shape parameters and the background air pressure.

[0040] Further, the step S3 includes step S32: calculating the distances from all grid points to the typhoon center point;

[0041] When the angle and distance conversion method parameter value is 0, the calculation method is:

[0042]

[0043] in Represents the i-th point To the center of the typhoon distance;

[0044] When the angle and distance conversion method parameter value is 1, the distances from all grid points to the typhoon center are calculated according to the calculation method in step S21.

[0045] Further, the step S3 includes step S33: calculating the relevant parameters, including: index parameters x and b, sea surface temperature T s , surface pressure P at the maximum wind speed radius s , vapor pressure q at 90% relative humidity s , virtual sea surface temperature T vs , air density ρ air ;

[0046] The calculation formulas for the relevant parameters are as follows:

[0047]

[0048] T vs =(T s +273.15)(1+0.61q s );

[0049]

[0050] where p c Indicates the central air pressure, p atm represents the background pressure, Indicates the rate of change of central air pressure, V o represents the moving speed of the typhoon center at the current moment. These parameters are calculated in step S31. Represents the latitude of the i-th point.

[0051] Further, the step S3 includes step S34: calculating the transition wind field and the annular wind field;

[0052] The calculation method of the transition wind field is as follows:

[0053]

[0054] in Represents the i-th point The transition wind field value, R mrepresents the maximum wind speed radius, which is calculated in step S31;

[0055] The annular wind field calculation method is as follows:

[0056]

[0057] in Represents the i-th point The annular wind field value, V R Indicates the maximum sustained wind speed at the current moment, which is calculated in step S31.

[0058] Further, the step S3 includes step S35: calculating grid point results, where the grid point results to be calculated include east-west wind speed, north-south wind speed, and air pressure;

[0059] Assume that the longitude and latitude coordinates of the grid point in the i-th row and j-th column are The distance to the center of the typhoon is d i,j , tr represents the typhoon range determination parameter, which can be obtained in step S1 and is calculated as follows:

[0060] When i,j / R m When ≤tr:

[0061]

[0062] When i,j / R m >tr:

[0063]

[0064] in Respectively represent the grid points in the i-th row and j-th column The east-west wind speed, north-south wind speed and air pressure; α is the wind speed adjustment parameter, which is obtained in step S1; u c Indicates the moving speed of the typhoon center at the current moment, θ c represents the moving direction of the typhoon center at the current moment, which can be calculated by step S31; φ is a known coefficient, and the calculation method is as follows:

[0065]

[0066] Thus, the east-west wind speed, north-south wind speed and air pressure results of all grid points are obtained.

[0067] A second object of the present invention is to provide a system based on the above-mentioned asymmetric parameterized wind field early warning method, comprising:

[0068] The database module is responsible for adding, deleting, modifying and querying all results of wind farm events;

[0069] The wind farm event module is responsible for constructing a wind farm event process and controlling the operation logic of the entire system;

[0070] Input processing module, responsible for reading input files and inputting wind farm event information into the system;

[0071] The calculation module is responsible for starting the asymmetric parameterized wind field warning system for calculation;

[0072] The output processing module is responsible for extracting and temporarily storing the calculation results after the model calculation is completed;

[0073] The analysis module is responsible for drawing according to the calculation results of the output processing module.

[0074] The advantages of the present invention are that it can simulate the evolution process of a typhoon and the accompanying time series changes in air pressure, wind speed, wind direction, etc., and can be used to effectively evaluate the scope and impact of disasters such as storm surges and typhoon waves in coastal areas of my country, so as to rationally allocate disaster relief resources, optimize the design standards of coastal defense projects, and accurately formulate disaster prevention and mitigation strategies to provide guarantees for the economic development of coastal cities. In the field of scientific research, the typhoon landing process is a hot topic in current typhoon research because it involves multi-phase coupled dynamic processes of sea, air and land and multi-scale interaction mechanisms. In this process, the accurate location of the landing site, the attenuation of typhoon intensity after approaching the coast and landing, the disappearance of the typhoon core on land or its return to the ocean to continue to operate after landing, and other scientific issues all rely on the support of a large amount of observational data and reanalysis data. BRIEF DESCRIPTION OF THE DRAWINGS

[0075] Figure 1 Schematic diagram of the relationship between modules of the asymmetric parameterized wind farm warning system in the embodiment;

[0076] Figure 2 It is a schematic diagram of the operation logic of the asymmetric parameterized wind farm early warning system in the embodiment;

[0077] Figure 3 Schematic diagram of the execution process of the calculation module in the embodiment;

[0078] Figure 4 Schematic diagram of the plane wind field distribution at 6 moments in the wind field simulation of Typhoon Ampil in the embodiment. DETAILED DESCRIPTION

[0079] The present invention is further described in detail below in conjunction with the embodiments. It should be understood that the directions or positional relationships indicated by the terms "upper", "lower", "front", "back", "left", "right", "top", "bottom", "inside" and "outside" herein are based on the directions or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific direction, be constructed and operated in a specific direction, and therefore cannot be understood as limiting the present invention.

[0080] like Figure 1 As shown, the asymmetric parameterized wind farm early warning system constructed in this embodiment mainly includes six modules, namely: database module, wind farm event module, input processing module, calculation module, output processing module and analysis module. Figure 2 The operating logic between the modules in the system is demonstrated. The following is a detailed introduction to the six modules.

[0081] The database module is responsible for adding, deleting, modifying and querying all results of wind farm events. Its main functions are as follows:

[0082] 1) Compare and query wind farm events in the database

[0083] After inputting wind farm change information into the system, the system will first pass the information to the database module, which will compare the information with historical wind farm events, find similar items, and output the stored data and results; if there is no similar item, it will output 0.

[0084] 2) Add new wind field event data and results

[0085] If there is no similar item in the database for the input wind field, the result will be obtained through system simulation calculation; the result will be added to the database module as new event information.

[0086] 3) Modify the specified wind field event data and results

[0087] For some historical simulation event data and results, there may be areas that need to be updated and adjusted, which can be achieved through this function.

[0088] 4) Delete the specified wind farm event data and results

[0089] Some historical simulation event data and results may need to be deleted, which can be achieved through this function.

[0090] The wind farm event module is responsible for constructing a wind farm event process and controlling the operation logic of the entire system. Its main functions are as follows:

[0091] 1) Accept input information and query it in the database module

[0092] The project module first accepts wind farm input information and system operation related information, and compares it with the wind farm event database based on the wind farm information. If there are similar projects, the project module directly extracts the data and results of similar projects in the database.

[0093] 2) Construct input processing module, execution module, output processing module and result analysis module

[0094] If there is no similar project in the database, it is necessary to calculate the wind field event process generated by the input wind field information; at this time, the project module will start to build the input processing module, execution module, output processing module and analysis module for subsequent calculations.

[0095] 3) Start running wind field event calculation

[0096] If there is no similar project in the database, after the project module is built, different modules will be started in sequence for calculation.

[0097] 4) Store the calculated project data and results in the wind farm event database

[0098] When the wind farm event process calculation is completed, the project module is responsible for storing the calculation results in the output processing module, as well as the pictures and animations in the analysis module into the wind farm event database.

[0099] The input processing module is responsible for reading the input file and inputting wind farm event information into the system. Its main functions are as follows:

[0100] 1) Read and organize basic characteristic information of wind field events

[0101] The basic characteristic information of the wind field includes: the number of slice data along the typhoon path, background air pressure, atmospheric density, wind speed adjustment parameters, and slice data; the slice data includes: slice data time, longitude and latitude of the typhoon center, maximum sustained wind speed, typhoon center air pressure, maximum wind speed radius, shape parameters, and input site information;

[0102] 2) Read and organize wind field event simulation information

[0103] The simulation information includes: simulation start time, simulation end time, time step (in seconds), simulation year, latitude and longitude resolution, latitude and longitude range, angle and distance conversion method parameters, typhoon range determination parameters, wind speed correction parameters, output point information, and output file name; the output point information includes: the number of output points and the latitude and longitude of the output points.

[0104] The execution process of the input processing module includes the following steps:

[0105] S1: Reading the input file includes the following process:

[0106] S11: Read the start time, end time, simulation time step (in seconds), simulation year, and the number of slice data; initialize the slice data time, typhoon center longitude and latitude, typhoon range determination parameters, maximum sustained wind speed, typhoon center pressure, maximum wind speed radius, and shape parameters according to the number of slice data;

[0107] S12: Read the grid longitude and latitude interval and the grid longitude and latitude range;

[0108] S13: Reading angle and distance conversion method parameters, background air pressure, atmospheric density, and wind speed adjustment parameters;

[0109] S14: Read output point information;

[0110] S15: Read wind speed correction parameters and output file name.

[0111] It should be noted that steps S11 to S15 must be performed in sequence. Once the sequence is incorrect or an error occurs in the intermediate reading, the system will stop directly.

[0112] The calculation module is responsible for starting the asymmetric parameterized wind field warning system for calculation. Its main functions are as follows:

[0113] 1) Set the number of computing nodes and cores

[0114] The module automatically identifies idle nodes and available cores of a computer or cluster and uses all available cores. In addition, users can also specify the nodes and cores to be used.

[0115] 2) Write and execute the compilation script

[0116] For each model run, the execution module will delete the previous compilation results and recompile to ensure that it is not affected by historical results;

[0117] 3) Write and execute the task startup script

[0118] Model calculation relies on the PBS job management system. The execution module is responsible for setting up and writing relevant scripts and starting the task for calculation.

[0119] Figure 3 The calculation module execution process is shown, which includes the following steps:

[0120] S2: Initialize wind field information and grid information;

[0121] S3: Circular calculation results.

[0122] In step S2, according to the input file information obtained in step S1, the simulated typhoon wind field information and grid distribution information are initialized for use in subsequent cyclic simulations. Step S2 includes the following steps:

[0123] S21: Initialize typhoon time difference sequence information and calculate the distance and direction of adjacent typhoon centers;

[0124] S22: Calculate the moving speed of the typhoon center;

[0125] S23: Initialize the grid distribution and related variables, and calculate the number of cycles.

[0126] The calculation method of step S21 is as follows:

[0127] First, based on all the typhoon slice data obtained in step S1, calculate the time difference between each slice data time and the start time in seconds, thereby obtaining the time difference sequence information of the typhoon slice. Calculate the distance and direction of the typhoon center of adjacent slice data, and the distance calculation formula is as follows:

[0128]

[0129] in Indicates the center of the typhoon To the center of the typhoon The distance between the two points, φ1 and φ2 represent the longitudes of the two points. Indicates latitude. R indicates the radius of the earth, in nmi. Arctan(r,z) indicates the angle of the right triangle with adjacent side length z and opposite side length r.

[0130]

[0131] The direction calculation formula is as follows:

[0132]

[0133] in Indicates the center of the typhoon To the center of the typhoon The azimuth of the direction is calculated in clockwise direction with due north as 0°. arctan(y,x) represents the angle of the right triangle with adjacent side length x and opposite side length y. rad2dg=0.017453293 represents the radian corresponding to each angle, which is used to convert angles to radians.

[0134] Step S22 obtains the moving speed of the typhoon center based on the distance and time difference between adjacent typhoon centers calculated in step S21, and the formula is as follows:

[0135]

[0136] in Indicates the center of the typhoon To the center of the typhoon The moving speed, Indicates that the typhoon center is at point The time difference between the time point and the initial time point, Indicates that the typhoon center is at point The time difference between the time point and the initial time point, so Indicates the typhoon center from point Move to time.

[0137] Step S23 initializes the grid information and related variables according to the simulation information read in step S1, and calculates the number of cycles. The grid information includes: the number of grid longitude and latitude grid points and the longitude and latitude division vectors; the related variables include: east-west wind speed, north-south wind speed, air pressure, latitude vector and longitude vector.

[0138] The number of cycles is calculated as follows:

[0139]

[0140] Where num represents the number of loops. Indicates the last typhoon center point The time difference between the time point and the initial time point, dt represents the time interval of calculation.

[0141] Step S3 performs a loop calculation based on the variable information obtained in steps S1 and S2, and performs the following steps in each loop:

[0142] S31: perform linear interpolation on some variables according to the current cycle time point;

[0143] S32: Calculate the distances from all grid points to the typhoon center;

[0144] S33: Calculate relevant parameters;

[0145] S34: Calculate the transition wind field and the annular wind field;

[0146] S35: Calculate grid point results.

[0147] Some of the variables in step S31 include: the longitude and latitude of the typhoon center at the current time, the maximum sustained wind speed, the central air pressure, the central air pressure change rate, the moving speed of the typhoon center at the current moment, the moving direction, the maximum wind speed radius, the shape parameters and the background air pressure.

[0148] In step S32, there are two algorithms for calculating the distance from all grid points to the typhoon center. The calculation method to be selected is determined according to the input parameter angle and the distance conversion method parameter. If the parameter value is 0, the calculation method is:

[0149]

[0150] Represents the i-th point To the center of the typhoon distance; if the parameter value is 1, the distance from all grid points to the typhoon center is calculated according to the calculation method in step S21.

[0151] In step S33, the relevant parameters calculated include: index parameters x and b, sea surface temperature T s , surface pressure P at the maximum wind speed radius s , vapor pressure q at 90% relative humidity s , virtual sea surface temperature T vs , air density ρ air The calculation formulas for the relevant parameters are as follows:

[0152]

[0153] where p c Indicates the central air pressure, p atm represents the background pressure, Indicates the rate of change of central air pressure, V o represents the moving speed of the typhoon center at the current moment. These parameters are calculated in step S31. Represents the latitude of the i-th point.

[0154] In step S34, the method for calculating the transition wind field is as follows:

[0155]

[0156] in Represents the i-th point The transition wind field value, R m represents the maximum wind speed radius, which is calculated in step S31. The calculation method of the annular wind field is as follows:

[0157]

[0158] in Represents the i-th point The annular wind field value, V R Indicates the maximum sustained wind speed at the current moment, which is calculated in step S31.

[0159] In step S35, the grid point results to be calculated include: east-west wind speed, north-south wind speed and air pressure. For the convenience of display, assume that the longitude and latitude coordinates of the grid point in the i-th row and j-th column of the grid are The distance to the center of the typhoon is d i,j , the calculation method is as follows:

[0160] tr represents the typhoon range determination parameter, which can be obtained in step S1. i,j / R m When ≤tr:

[0161]

[0162] When i,j / R m >tr:

[0163]

[0164] in Respectively represent the grid points in the i-th row and j-th column The east-west wind speed, north-south wind speed and air pressure are obtained, and α is the wind speed adjustment parameter, which is obtained in step S1. c Indicates the moving speed of the typhoon center at the current moment, θ c Indicates the moving direction of the typhoon center at the current moment, which can be calculated in step S31. γ is a known coefficient, and the calculation method is as follows:

[0165]

[0166] Thus, the east-west wind speed, north-south wind speed and air pressure results of all grid points are obtained.

[0167] The output processing module is responsible for extracting and temporarily storing the calculation results after the model calculation is completed.

[0168] The analysis module is responsible for drawing according to the calculation results of the output processing module. Its main functions are as follows:

[0169] 1) Draw the wind field information in the determined area, including the occurrence and movement of typhoons;

[0170] 2) Draw the changes in air pressure at any location in the area;

[0171] 3) Draw the wind speed changes at any location in the area.

[0172] For ease of understanding, this embodiment shows the simulation results of Typhoon Ampil in July 2018. Figure 4It shows that Typhoon Ampil experienced a process of wind circle range from large to small in the entire time series of its formation, development and extinction, and the wind speed reached the maximum before and after landing. Typhoon Ampil showed a strong asymmetric relationship at the beginning. When it approached the land on July 21, the shape of the wind circle was approximately circular, but the center was still off the center of the circle. Then, the wind speed contour lines at 02:00 and 13:00 on July 22 showed that due to the weakening effect of the terrain on the wind field, the wind speed after landing was significantly weakened, and the high wind speed area was reduced, showing a strong asymmetry, which is consistent with the change trend of the maximum wind speed radius.

[0173] The above embodiments are only used to explain the concept of the present invention, rather than to limit the protection of the rights of the present invention. Any non-substantial changes to the present invention using this concept shall fall within the protection scope of the present invention.

Claims

1. An asymmetric parameterized wind field early warning method, characterized in that: include: S1. Read the input file, including the following steps: S11. Read the start time, end time, simulation time step, simulation year, read the number of slice data, and initialize the slice data time, typhoon center longitude and latitude, typhoon range determination parameters, maximum sustained wind speed, typhoon center pressure, maximum wind speed radius, and shape parameters according to the number of slice data; S12. Read the grid longitude and latitude interval and grid longitude and latitude range; S13. Reading angle and distance conversion method parameters, background pressure, atmospheric density, and wind speed adjustment parameters; S14. Read output point information; S15. Read wind speed correction parameters and output file name; S2. Initialize the simulated typhoon wind field information and grid distribution information according to the input file information obtained in step S1 for use in subsequent cyclic simulations, including the following steps: S21. Initialize typhoon time difference sequence information and calculate the distance and direction of adjacent typhoon centers; S22. Calculate the moving speed of the typhoon center; S23. Initialize the grid distribution and variable information, and calculate the number of cycles; S3. Perform loop calculation based on the variable information obtained in steps S1 and S2; S4. extract and temporarily store the calculation results; S5. Draw according to the calculation results.

2. The asymmetric parameterized wind field early warning method according to claim 1, characterized in that: The calculation method of step S21 is as follows: First, according to all the typhoon slice data obtained in step S1, the time difference between each slice data time and the start time is calculated in seconds, thereby obtaining the time difference sequence information of the typhoon slice; Calculate the distance and direction of the typhoon center of adjacent slice data. The distance calculation formula is as follows: in Indicates the center of the typhoon To the center of the typhoon , φ1 and φ2 represent the longitudes of the two points. represents latitude, R represents the radius of the earth in nmi, arctan(r,z) represents the angle of the right triangle with adjacent side length z and opposite side length r, where The direction calculation formula is as follows: in Indicates the center of the typhoon To the center of the typhoon The azimuth in radians is calculated clockwise with due north as 0°. arctan(y,x) represents the angle of the right triangle with adjacent side length x and opposite side length y. rad2dg=0.017453293 represents the radian corresponding to each angle, which is used to convert angles to radians.

3. The asymmetric parameterized wind field early warning method according to claim 2, characterized in that: According to the distance and time difference between adjacent typhoon centers calculated in step S21, the moving speed of the typhoon center can be obtained, and the formula is as follows: in Indicates the center of the typhoon To the center of the typhoon The moving speed, Indicates that the typhoon center is at point The time difference between the time point and the initial time point, Indicates that the typhoon center is at point The time difference between the time point and the initial time point, so Indicates the typhoon center from point Move to time.

4. The asymmetric parameterized wind field early warning method according to claim 3, characterized in that: Step S23 initializes grid information and related variables according to the simulation information read in step S1, and calculates the number of cycles, wherein the grid information includes the number of grid longitude and latitude grid points and the longitude and latitude division vectors, and the related variables include east-west wind speed, north-south wind speed, air pressure, latitude vector, and longitude vector; The number of cycles is calculated as follows: Where num represents the number of loops. Indicates the last typhoon center point The time difference between the time point and the initial time point, dt represents the time interval of calculation.

5. The asymmetric parameterized wind field early warning method according to claim 1, characterized in that: The step S3 includes step S31: performing linear interpolation on some variables according to the current cycle time point, and some variables include: the longitude and latitude of the typhoon center at the current time, the maximum sustained wind speed, the central air pressure, the central air pressure change rate, the moving speed of the typhoon center at the current moment, the moving direction, the maximum wind speed radius, the shape parameters and the background air pressure.

6. The asymmetric parameterized wind field early warning method according to claim 5, characterized in that: The step S3 includes step S32: calculating the distances from all grid points to the typhoon center point; When the angle and distance conversion method parameter value is 0, the calculation method is: in Represents the i-th point To the center of the typhoon distance; When the angle and distance conversion method parameter value is 1, the distances from all grid points to the typhoon center are calculated according to the calculation method in step S21.

7. The asymmetric parameterized wind field early warning method according to claim 6, characterized in that: The step S3 includes step S33: calculating the relevant parameters, including: index parameters x and b, sea surface temperature T s , surface pressure P at the maximum wind speed radius s , vapor pressure q at 90% relative humidity s , virtual sea surface temperature T vs , air density ρ air ; The calculation formulas for the relevant parameters are as follows: T vs =(T s +273.15)(1+0.61q s ); where p c Indicates the central air pressure, p atm represents the background pressure, Indicates the rate of change of central air pressure, V o represents the moving speed of the typhoon center at the current moment. These parameters are calculated in step S31. Represents the latitude of the i-th point.

8. The asymmetric parameterized wind field early warning method according to claim 7, characterized in that: The step S3 includes step S34: calculating the transition wind field and the annular wind field; The calculation method of the transition wind field is as follows: in Represents the i-th point The transition wind field value, R m represents the maximum wind speed radius, which is calculated in step S31; The annular wind field calculation method is as follows: in Represents the i-th point The annular wind field value, V R Indicates the maximum sustained wind speed at the current moment, which is calculated in step S31.

9. The asymmetric parameterized wind field early warning method according to claim 8, characterized in that: The step S3 includes step S35: calculating grid point results, where the grid point results to be calculated include east-west wind speed, north-south wind speed, and air pressure; Assume that the longitude and latitude coordinates of the grid point in the i-th row and j-th column are The distance to the center of the typhoon is d i,j , tr represents the typhoon range determination parameter, which can be obtained in step S1 and is calculated as follows: When i,j / R m When ≤tr: When i,j / R m >tr: in Respectively represent the grid points in the i-th row and j-th column The east-west wind speed, north-south wind speed and air pressure; α is the wind speed adjustment parameter, which is obtained in step S1; u c Indicates the moving speed of the typhoon center at the current moment, θ c represents the moving direction of the typhoon center at the current moment, which can be calculated by step S31; γ is a known coefficient, and the calculation method is as follows: Thus, the east-west wind speed, north-south wind speed and air pressure results of all grid points are obtained.

10. A system based on the above asymmetric parameterized wind field early warning method, characterized in that: include: The database module is responsible for adding, deleting, modifying and querying all results of wind farm events; The wind farm event module is responsible for constructing a wind farm event process and controlling the operation logic of the entire system; Input processing module, responsible for reading input files and inputting wind farm event information into the system; The calculation module is responsible for starting the asymmetric parameterized wind field warning system for calculation; The output processing module is responsible for extracting and temporarily storing the calculation results after the model calculation is completed; The analysis module is responsible for drawing according to the calculation results of the output processing module.

Images