A typhoon wave field simulation method, device, equipment and medium

Abstract

The present application discloses a typhoon wave field simulation method, device, equipment and medium, relating to the technical field of typhoon wave simulation. The method includes: tracking the typhoon center of a target typhoon by using a polar coordinate grid, where the typhoon center is the pole origin of the polar coordinate grid; constructing a typhoon forcing field for each wave ray in the polar coordinate grid, and determining the reference wind speed on the corresponding wave ray according to the typhoon forcing field; constructing a typhoon wave field model by using a two-dimensional wave ray theory parameterization model; on the polar coordinate grid, solving the typhoon wave field model based on the initial field to determine the wave spectrum peak frequency, the main wave direction of the sea wave and the sea wave energy changing with time on each wave ray; the initial field includes the reference wind speed, the wave spectrum peak frequency, the main wave direction of the sea wave and the sea wave energy at each grid point at the initial moment, and the grid point is the center point of the corresponding grid. The present application improves the prediction speed of typhoon wave field data.

Description

Technical Field

[0001] The present application relates to the technical field of typhoon wave simulation, and particularly to a typhoon wave field simulation method, device, equipment and medium. Background Art

[0002] A tropical cyclone is a low-pressure vortex occurring between the ocean surface and the atmospheric system in the tropics and subtropics. It is characterized by strong wind speed, low central pressure, and is accompanied by high-intensity rainfall. It is one of the most dangerous and destructive natural disaster phenomena. To ensure navigation safety and prevent coastal disasters, it is necessary to quickly and effectively evaluate the wave field under typhoon control.

[0003] In the Northern Hemisphere, due to the translation and rotation characteristics of typhoons, the waves in the right sector of the typhoon will move with the typhoon under the action of high wind speed and last for a period of time (typhoons rotate counterclockwise in the Northern Hemisphere and clockwise in the Southern Hemisphere). The resulting typhoon waves have extremely strong energy, and the reason for this phenomenon is called group velocity resonance. Due to the rapid changes in the typhoon wind field in space and time and the influence of the group velocity resonance effect, the simulation and prediction of the wave field under typhoon action is a challenging task. Currently, the third-generation wave model can relatively reliably simulate the spatial distribution of the wave field under typhoon. However, due to the wind field data commonly used in the driving model, which has a low resolution due to computational cost limitations, it is difficult to correctly distinguish the rapidly changing wind field and the large local wind speed gradient within the typhoon. Therefore, it is of great significance to establish a reliable typhoon wave field model and achieve rapid prediction. Summary of the Invention

[0004] The purpose of the present application is to provide a typhoon wave field simulation method, device, equipment and medium, which improves the prediction speed of typhoon wave field data.

[0005] To achieve the above purpose, the present application provides the following solutions:

[0006] In the first aspect, the present application provides a typhoon wave field simulation method, including:

[0007] Tracking the typhoon center of the target typhoon using a polar coordinate grid, where the typhoon center is the pole origin of the polar coordinate grid;

[0008] Constructing a typhoon forcing field for each wave ray in the polar coordinate grid, and determining the reference wind speed on the corresponding wave ray according to the typhoon forcing field;

[0009] Constructing a typhoon wave field model using a two-dimensional wave ray theory parameterization model;

[0010] On the polar coordinate grid, the typhoon wave field model is solved based on the initial field to determine the time-varying wave spectral peak frequency, the main wave direction of the ocean wave, and the ocean wave energy on each wave ray; the initial field includes the reference wind speed, the wave spectral peak frequency, the main wave direction of the ocean wave, and the ocean wave energy at each grid point at the initial moment, and the grid point is the center point of the corresponding grid.

[0011] In a second aspect, the present application provides a typhoon wave field simulation device, which includes:

[0012] A relationship determination module between the typhoon center and the polar coordinate grid, configured to track the typhoon center of the target typhoon using the polar coordinate grid, and the typhoon center is the polar coordinate origin of the polar coordinate grid;

[0013] A typhoon forcing field construction module, configured to construct a typhoon forcing field for each wave ray in the polar coordinate grid, and determine the reference wind speed on the corresponding wave ray according to the typhoon forcing field;

[0014] A typhoon wave field model construction module, configured to construct a typhoon wave field model using a two-dimensional wave ray theory parameterization model;

[0015] A solution module, configured to solve the typhoon wave field model on the polar coordinate grid based on the initial field to determine the time-varying wave spectral peak frequency, the main wave direction of the ocean wave, and the ocean wave energy on each wave ray; the initial field includes the reference wind speed, the wave spectral peak frequency, the main wave direction of the ocean wave, and the ocean wave energy at each grid point at the initial moment, and the grid point is the center point of the corresponding grid.

[0016] In a third aspect, the present application provides a computer device, including: a memory, a processor, and a computer program stored on the memory and executable on the processor, and the processor executes the computer program to implement the steps of the typhoon wave field simulation method described in any one of the above.

[0017] In a fourth aspect, the present application provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the typhoon wave field simulation method described in any one of the above are implemented.

[0018] According to the specific embodiments provided by the present application, the following technical effects are disclosed by the present application:

[0019] The present application provides a typhoon wave field simulation method, device, equipment and medium. The typhoon center of the target typhoon is tracked by using a polar coordinate grid, and the typhoon center is the pole origin of the polar coordinate grid. The polar coordinate grid encrypts the high wind speed area, highlighting the key area while reducing the calculation amount and improving the calculation speed. The typhoon wave field model is constructed by using a two-dimensional wave ray theory parameterization model. The typhoon wave field model requires fewer parameters and less calculation amount, further improving the calculation speed, and thus improving the forecasting speed of typhoon wave field data. Brief Description of the Drawings

[0020] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings required in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0021] Figure 1 It is a schematic flow chart of a typhoon wave field simulation method provided by an embodiment of the present application;

[0022] Figure 2 It is a schematic principle diagram of a typhoon wave field simulation method provided by an embodiment of the present application;

[0023] Figure 3 It is a schematic diagram of a polar coordinate grid provided by an embodiment of the present application;

[0024] Figure 4 It is a schematic diagram of the simulation effect of a typhoon wave field provided by an embodiment of the present application;

[0025] Figure 5 It is a schematic diagram of the effect of verifying the typhoon wave field model by using the measured sea wave data of an altimeter provided by an embodiment of the present application;

[0026] Figure 6 It is a schematic diagram of the effect of verifying the typhoon wave field model by using the Holland wind field to drive the SWAN sea wave model provided by an embodiment of the present application;

[0027] Figure 7 It is a schematic diagram of the structure of a computer device provided by an embodiment of the present application. Detailed Description of the Embodiments

[0028] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present application without creative efforts shall fall within the protection scope of the present application.

[0029] To make the above objects, features, and advantages of the present application more obvious and understandable, the present application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0030] The present application provides a typhoon wave field simulation method, as Figure 1 and Figure 2 shown, the typhoon wave field simulation method includes steps 101 to 104.

[0031] Step 101: Use a polar coordinate grid to track the typhoon center of the target typhoon, and the typhoon center is the pole origin of the polar coordinate grid.

[0032] Step 102: Construct a typhoon forcing field for each wave ray in the polar coordinate grid, and determine the reference wind speed on the corresponding wave ray according to the typhoon forcing field.

[0033] Step 103: Use a two-dimensional wave ray theory parameterization model to construct a typhoon wave field model.

[0034] Step 104: On the polar coordinate grid, solve the typhoon wave field model based on the initial field to determine the time-varying wave spectrum peak frequency, main wave direction, and wave energy on each wave ray; the initial field includes the reference wind speed, wave spectrum peak frequency, main wave direction, and wave energy at each grid point at the initial moment, and the grid point is the center point of the corresponding grid.

[0035] The present application uses a polar coordinate grid to track the typhoon center of the target typhoon. The typhoon center is the pole origin of the polar coordinate grid. The polar coordinate grid will encrypt the high wind speed area, highlighting the key area while reducing the calculation amount and improving the calculation speed. A two-dimensional wave ray theory parameterization model is used to construct a typhoon wave field model. The typhoon wave field model requires fewer parameters and less calculation amount, further improving the calculation speed, thereby improving the prediction speed of typhoon wave field data.

[0036] During the movement of the typhoon, the range of typhoon waves is almost the same as the range of the typhoon wind field. And since the wave ray propagation trajectory is a ray, and the dominant waves in a large part of the area under the action of the typhoon are generated by the high wind speed area, using a polar coordinate grid will encrypt this area, highlighting the key area while reducing the calculation amount.

[0037] In an exemplary embodiment, taking Typhoon Linfa as an example, a typhoon wave field simulation method of the present application is described. In the polar coordinate grid, the typhoon moving direction is the starting azimuth angle , the grid azimuth angle resolution is , the radial spacing of the grid is , the grid radial distribution starts from and ends at . is the radial spacing of the grid, is the radius difference between the inner circle and the outer circle where the grid is located, is the radius of the outer circle where the grid is located, is the radius of the smallest circle in the polar coordinate grid, is the radius of the largest circle in the polar coordinate grid, , , is the maximum wind speed radius of the target typhoon. The polar coordinate grid is as shown in Figure 3 .

[0038] In an exemplary embodiment, although there are individual differences in the spatial distribution of the typhoon wind field, it can generally be approximated as a translational vortex. The typhoon forcing field of the present application is the Holland typhoon wind field model (Holland wind field). This Holland typhoon wind field model is a circular symmetric wind field model based on the gradient wind equation. The expression of the Holland typhoon wind field model is: .

[0039] Among them, is the gradient wind speed outside the atmospheric boundary layer at a radius from the typhoon center, is the environmental air pressure far from the typhoon, that is, the sea surface air pressure without typhoon, is the air pressure at the typhoon center, is the air density, is the Coriolis parameter, and are both dimensionless parameters, and have the relationship with the maximum wind speed radius of the target typhoon as . The dimensionless parameter defines the shape of the typhoon vortex, is the maximum wind speed radius of the target typhoon.

[0040] In an exemplary embodiment, in order to calculate the reference wind speed , the gradient wind speed Correction needs to be carried out under the assumption of the roughness of the boundary layer profile surface, and the reference wind speed on the corresponding wave ray is determined according to the typhoon forcing field, specifically including: correcting the gradient wind speed by using a correction formula to obtain the reference wind speed, and the reference wind speed is the wind speed at the standard reference height of 10m.

[0041] The correction formula is expressed as: ; where is the reference wind speed, is a constant.

[0042] The wind speed is in a spiral shape pointing to the typhoon center, and the average inflow angle where , , are all preset parameters.

[0043] The IBTrACS (International Best Track Archive for Climate Stewardship) typhoon best track dataset is used to generate the typhoon forcing field. IBTrACS provides data from twelve global institutions or historical databases and is currently the most complete global historical typhoon dataset. This dataset includes typhoon information such as the longitude and latitude of the typhoon center every three hours, the minimum central pressure, the maximum wind speed near the center, and the radius of the maximum wind speed. In this embodiment, the typhoon information provided by its latest version v04 from the National Oceanic and Atmospheric Administration (NOAA) of the United States is used to generate the wind field driving the model.

[0044] The typhoon information provided by IBTrACS is once every 3h, so the Holland wind field driving the model is also set to change once every 3h. During this period, the typhoon wind field only has a translation in position, and the range and wind speed remain unchanged. The translation speed is calculated based on the typhoon positions at two moments. The result of the first 3h operation of the Holland typhoon wind field model is used as the initial field for the next 3h operation to continue running.

[0045] According to the frequency of the provided typhoon information, the Holland wind field driving the model also makes corresponding changes. During the same wind field driving period, the typhoon wind field only has a translation in position, and the range and wind speed remain unchanged. The translation speed is calculated based on the typhoon positions at two moments.

[0046] The two-dimensional wave ray theory parameterization model starts from the energy and momentum conservation equations and derives the self-similar relationships of the spectral peak frequency, main wave direction, and energy of the ocean waves based on the ocean wave growth relationship under a limited wind area.

[0047] In an exemplary embodiment, the typhoon wave field model is expressed as the following four differential equations.

[0048] (1) , this equation represents the change in the wave energy e.

[0049] (2) , this equation represents the change in the spectral peak frequency. The non - linear interaction of ocean waves causes the spectral peak frequency to decrease.

[0050] (3) , this equation reflects the change in the wave direction .

[0051] (4) , this equation represents the propagation of the wave ray over time.

[0052] Among them, t represents time, is the mean group velocity, e represents the wave energy, is the spectral peak frequency of the ocean wave, is the wind energy input, and the wind energy input is determined according to the reference wind speed, is the energy dissipation, is the spectral peak group velocity, is a constant, is the acceleration due to gravity, is the main wave direction of the ocean wave, represents the grid - point position, is a function of the wave age, representing the change in the spectral peak frequency caused by the wind energy input, represents the energy change caused by the vertical divergence of the group velocity, and are both constants, is the spectral peak wave number, represents a step function related to the reciprocal of the wave age and , represents the wind direction, represents the unit vector in the direction of the group velocity.

[0053] The spectral peak phase velocity =g / . The spectral peak group velocity is 0.5 times the spectral peak phase velocity, and g is the acceleration due to gravity. Therefore, the two parameters of the spectral peak group velocity and the spectral peak frequency can be directly calculated.

[0054] The relationship between the mean group velocity and the spectral peak group velocity can be obtained through the JONSWAP spectrum (Joint North Sea Wave Project Spectrum): ; , ; and Represent the dimensionless forms of wind energy input and energy dissipation respectively, , , which dominate the evolution of wave energy, is the dimensional wind energy input, is the dimensional energy dissipation.

[0055] Since the wave spectrum generated by typhoon forcing conforms to the traditional finite fetch JONSWAP spectrum type and has "spectral type stability", and and other coefficients can all be derived using the JONSWAP spectrum.

[0056] In an exemplary embodiment, the initial field is . In the initial field, the main wave direction of the waves is the direction corresponding to the reference wind speed, and the wave spectral peak frequency and the wave energy are obtained according to the wind wave growth relationship in finite wind duration.

[0057] The wave spectral peak frequency and the wave energy are expressed according to the wind wave growth relationship in finite wind duration as: , .

[0058] Wherein, is the dimensionless wave spectral peak frequency, is the dimensionless wave energy, is the dimensionless wind duration, t represents time, u represents wind speed, g represents gravitational acceleration, , , and are all constants. , , and can all be obtained from the wind wave growth relationship in finite fetch.

[0059] In an exemplary embodiment, during the solution process of the differential equation, with a relatively short wind duration , the spectral peak frequency and energy at the initial moment are obtained. The main wave direction of the waves at the initial moment is the same as the wind direction at that place. Since a two-dimensional wave ray theory parameterization model is adopted on the polar coordinate grid to construct the typhoon wave field model, the result of the field is composed of a large number of wave rays. During the operation of the typhoon wave field model, the wave rays will cross and overlap, and there may be multiple values in one grid. The wave parameters of the wave ray with the largest wavelength are selected as the values at that grid point. The typhoon wave field model uses the fourth-order Runge-Kutta method to solve the differential equation under the polar coordinate grid, and the model time step is set to 300 s.

[0060] Taking the typhoon "Lekima" (2009) as an example, the method of this application is used for simulation, and the results are asFigure 4 As shown Figure 4 In parts (a), (b), (c), (d), (e) and (f), the simulated effects of typhoon wave fields at times 0h, 6h, 12h, 18h, 24h and 30h are shown in sequence. Figure 4 In is the significant wave height of the typhoon wave. The typhoon wave field model fully demonstrates the spatial structure of the typhoon wave field and the development process of the sea waves. The typhoon wave field model has the following characteristics: (1) It has asymmetry in structure. Due to the rotation and translation characteristics of the typhoon, the group velocity resonance phenomenon causes the typhoon wave field to have obvious asymmetry. (2) The sea waves generated in the high wind speed area of the typhoon dominate the entire wave field. In the high wind speed area of the typhoon, due to the high wind energy input obtained by the sea waves, the wind waves in this area have the characteristics of high energy and high group velocity, and continuously diverge outward to the low wind speed area to form swell; in the low wind speed area of the typhoon, the swell propagated from the high wind speed area is mixed with the local wind waves, and the swell, as the dominant wave, is the main component of the mixed wave field in this area. In addition, the typhoon wave field model was verified using the measured sea wave data of the altimeter. As Figure 5 shown Figure 5 In part (a), the wind speed comparison effect is shown. Figure 5 In part (b), the significant wave height comparison effect is shown. Figure 5 In part (c), the effect diagram of the correlation coefficient R between the satellite significant wave height and the typhoon wave field model is shown. Generally, it is in good agreement with the altimeter data. Using the same Holland wind field as the typhoon wave field model to drive the SWAN (Simulating WAves Nearshore) model, the results are as Figure 6 shown Figure 6 In part (a), the comparison effect of satellite data, two-dimensional wave rays and the SWAN model is shown. In part (b), the correlation coefficient effect between the satellite data significant wave height and the SWAN model significant wave height is shown. The results of the typhoon wave field model are similar to those of the SWAN, proving that this method has high credibility. At the same time, the computational amount is much lower than that of the SWAN model, indicating that this typhoon wave field model can quickly and effectively simulate the typhoon wave field. The typhoon wave field model has a short running time and can be run on a personal computer. It can be used as an alternative to the third-generation sea wave spectrum model when the computing resources are insufficient.

[0061] This application adopts the above-mentioned typhoon wave field simulation method based on the two-dimensional wave ray theory. The parameters required for calculation are only the central pressure, maximum wind speed radius and position information of the typhoon at each moment, and the spatio-temporal distributions of the sea wave energy, main wave direction and spectral peak frequency of the typhoon wave can be predicted and simulated.

[0062] In an exemplary embodiment, on the polar coordinate grid, the typhoon wave field model is solved based on the initial field to determine the time-varying wave spectral peak frequency, the main wave direction of the ocean waves, and the ocean wave energy on each wave ray. Specifically, it includes: using the fourth-order Runge-Kutta method to solve the typhoon wave field model to determine the time-varying wave spectral peak frequency, the main wave direction of the ocean waves, and the ocean wave energy on each wave ray.

[0063] Based on the same inventive concept, an embodiment of the present application further provides a typhoon wave field simulation device for implementing the above-mentioned typhoon wave field simulation method. The solution provided by this device to solve the problem is similar to the solution described in the above method. Therefore, the specific limitations in one or more embodiments of the typhoon wave field simulation device provided below can refer to the limitations on the typhoon wave field simulation method in the above text, and will not be repeated here.

[0064] In an exemplary embodiment, the present application provides a typhoon wave field simulation device, including:

[0065] A typhoon center and polar coordinate grid relationship determination module, configured to use the polar coordinate grid to track the typhoon center of the target typhoon, where the typhoon center is the polar coordinate origin of the polar coordinate grid.

[0066] A typhoon forcing field construction module, configured to construct a typhoon forcing field for each wave ray in the polar coordinate grid, and determine the reference wind speed on the corresponding wave ray according to the typhoon forcing field.

[0067] A typhoon wave field model construction module, configured to construct a typhoon wave field model using a two-dimensional wave ray theory parameterization model.

[0068] A solution module, configured to solve the typhoon wave field model based on the initial field on the polar coordinate grid to determine the time-varying wave spectral peak frequency, the main wave direction of the ocean waves, and the ocean wave energy on each wave ray; the initial field includes the reference wind speed, wave spectral peak frequency, main wave direction of the ocean waves, and ocean wave energy at each grid point at the initial moment, and the grid point is the center point of the corresponding grid.

[0069] In an exemplary embodiment, a computer device is provided. The computer device can be a server or a terminal, and its internal structure diagram can be as Figure 7As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O), and a communication interface. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store typhoon wave field simulation data. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, it implements a typhoon wave field simulation method.

[0070] Those skilled in the art can understand that Figure 7 the structure shown in is only a block diagram of some structures related to the solution of this application, and does not constitute a limitation on the computer device to which the solution of this application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have different component arrangements. In an exemplary embodiment, a computer device is provided, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

[0071] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by the processor, the steps in the above method embodiments are implemented.

[0072] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use, and processing of relevant data need to comply with relevant regulations.

[0073] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, database, or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memories. Non-volatile memories can include read-only memory (ROM), magnetic tapes, floppy disks, flash memories, optical memories, high-density embedded non-volatile memories, resistive random access memories (ReRAM), magnetoresistive random access memories (MRAM), ferroelectric random access memories (FRAM), phase change memories (PCM), graphene memories, etc. Volatile memories can include random access memory (RAM) or external cache memories, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0074] The databases involved in the embodiments provided in the present application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., without limitation. The processors involved in the embodiments provided in the present application can be general-purpose processors, central processors, graphics processors, digital signal processors, data processing logics of programmable logics, etc., without limitation.

[0075] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, not all possible combinations of the technical features in the above embodiments are described. However, as long as these technical feature combinations do not conflict, they should all be considered to be within the scope described in this specification.

[0076] Specific examples are used in this article to elaborate on the principles and implementation manners of the present application. The descriptions of the above embodiments are only used to help understand the method and its core idea of the present application; at the same time, for those of ordinary skill in the art, according to the idea of the present application, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present application.

Claims

1. A typhoon wave field simulation method, characterized in that: The typhoon wave field simulation method comprises: Tracking the typhoon center of the target typhoon using a polar coordinate grid, wherein the typhoon center is the polar coordinate origin of the polar coordinate grid; constructing a typhoon forcing field for each wave ray in the polar coordinate grid, and determining a reference wind speed on the corresponding wave ray according to the typhoon forcing field; A typhoon wave field model is constructed using a two-dimensional wave ray theory parameterized model; On the polar coordinate grid, the typhoon wave field model is solved based on the initial field to determine the wave spectrum peak frequency, wave main wave direction and wave energy that change with time on each wave ray; the initial field includes the reference wind speed, wave spectrum peak frequency, wave main wave direction and wave energy of each grid point at the initial moment, and the grid point is the center point of the corresponding grid; during the operation of the typhoon wave field model, the wave rays will cross and overlap, and when there are multiple wave parameters in a grid, the wave parameter of the wave ray with the largest wavelength is selected as the wave parameter at the grid point; the wave parameters include wave spectrum peak frequency, wave main wave direction and wave energy; The typhoon wave field model is expressed as: ； ； ； ； Where t represents time, is the average group velocity, e represents the wave energy, is the peak frequency of the wave spectrum, is wind energy input, the wind energy input is determined according to a reference wind speed, For energy dissipation, is the peak group velocity, is a constant, is the acceleration due to gravity, is the main wave direction of the ocean wave, represents the grid point location, is a function of wave age, represents the energy change caused by the vertical divergence of the group velocity, is a constant, is the peak wave number, It is the inverse of the wave age. and The step function involved is Indicates wind direction. The unit vector representing the direction of the group velocity; In the polar coordinate grid, the grid azimuth resolution is , the radial spacing of the grid is , the radial distribution of the grid is from Start to Finish, is the radial spacing of the grid, is the radius difference between the inner circle and the outer circle where the grid is located, is the radius of the outer circle where the grid is located, is the radius of the smallest circle in the polar grid, is the radius of the largest circle in the polar grid, , , is the maximum wind speed radius of the target typhoon.

2. The typhoon wave field simulation method according to claim 1, characterized in that: The typhoon forcing field is the Holland typhoon wind field model, and the expression of the Holland typhoon wind field model is: ； in, The radius from the typhoon center The gradient wind speed outside the atmospheric boundary layer, is the ambient air pressure on the sea surface when there is no typhoon. is the air pressure at the center of the typhoon, is the air density, is the Coriolis parameter, and are dimensionless parameters, and The relationship with the maximum wind speed radius of the target typhoon is: , is the maximum wind speed radius of the target typhoon.

3. The typhoon wave field simulation method according to claim 2, characterized in that: Determining the reference wind speed on the corresponding wave ray according to the typhoon forcing field specifically includes: The gradient wind speed is corrected by using a correction formula to obtain the reference wind speed; the reference wind speed is the wind speed at a reference height of 10 m; The correction formula is expressed as: ;in, is the reference wind speed, is a constant.

4. The typhoon wave field simulation method according to claim 1, characterized in that: In the initial field, the main wave direction of the sea wave is the direction corresponding to the reference wind speed, and the sea wave spectrum peak frequency and the sea wave energy are obtained according to the wind-wave growth relationship in finite wind time; The wave spectrum peak frequency and the wave energy are expressed as follows according to the relationship between wind and wave growth in finite wind conditions: , ; in, is the dimensionless peak frequency of the wave spectrum, is the dimensionless wave energy, For dimensionless wind, t represents time, u represents wind speed, and g represents gravitational acceleration. , , and are all constants.

5. The typhoon wave field simulation method according to claim 1, characterized in that: On the polar coordinate grid, solving the typhoon wave field model based on the initial field to determine the wave spectrum peak frequency, wave main wave direction and wave energy that vary with time on each wave ray specifically includes: The fourth-order Runge-Kutta method is used to solve the typhoon wave field model to determine the wave spectrum peak frequency, wave main wave direction and wave energy that vary with time on each wave ray.

6. A typhoon wave field simulation device, characterized in that: The typhoon wave field simulation device comprises: A module for determining the relationship between a typhoon center and a polar coordinate grid, for tracking the typhoon center of a target typhoon using a polar coordinate grid, wherein the typhoon center is the polar coordinate origin of the polar coordinate grid; A typhoon forcing field construction module, used to construct a typhoon forcing field for each wave ray in the polar coordinate grid, and determine a reference wind speed on the corresponding wave ray according to the typhoon forcing field; Typhoon wave field model construction module, used to construct a typhoon wave field model using a two-dimensional wave ray theory parameterized model; A solution module is used to solve the typhoon wave field model on the polar coordinate grid based on the initial field to determine the wave spectrum peak frequency, wave main wave direction and wave energy that change with time on each wave ray; the initial field includes the wind speed, wave spectrum peak frequency, wave main wave direction and wave energy of each grid point at the initial moment, and the grid point is the center point of the corresponding grid; during the operation of the typhoon wave field model, the wave rays will cross and overlap, and when there are multiple wave parameters in a grid, the wave parameter of the wave ray with the largest wavelength is selected as the wave parameter at the grid point; the wave parameters include wave spectrum peak frequency, wave main wave direction and wave energy; The typhoon wave field model is expressed as: ； ； ； ； Where t represents time, is the average group velocity, e represents the wave energy, is the peak frequency of the wave spectrum, is wind energy input, the wind energy input is determined according to a reference wind speed, For energy dissipation, is the peak group velocity, is a constant, is the acceleration due to gravity, is the main wave direction of the ocean wave, represents the grid point location, is a function of wave age, represents the energy change caused by the vertical divergence of the group velocity, is a constant, is the peak wave number, It is the inverse of the wave age. and The step function involved is Indicates wind direction. The unit vector representing the direction of the group velocity; In the polar coordinate grid, the grid azimuth resolution is , the radial spacing of the grid is , the radial distribution of the grid is from Start to Finish, is the radial spacing of the grid, is the radius difference between the inner circle and the outer circle where the grid is located, is the radius of the outer circle where the grid is located, is the radius of the smallest circle in the polar grid, is the radius of the largest circle in the polar grid, , , is the maximum wind speed radius of the target typhoon.

7. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the typhoon wave field simulation method according to any one of claims 1 to 5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the typhoon wave field simulation method according to any one of claims 1 to 5 is implemented.

Images