Typhoon wave field simulation method, device, equipment and medium

Abstract

The invention discloses a typhoon wave field simulation method, device and equipment and a medium, and relates to the technical field of typhoon wave simulation, and the method comprises the steps: employing a polar coordinate grid to track the typhoon center of a target typhoon, and the typhoon center is the polar coordinate origin of the polar coordinate grid; constructing a typhoon forced 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 forced field; constructing a typhoon wave field model by adopting a two-dimensional wave ray theory parameterized model; on the polar coordinate grid, solving the typhoon wave field model based on an initial field to determine the sea wave spectrum peak frequency, the sea wave main wave direction and the sea wave energy which change along with time on each wave ray; the initial field comprises the reference wind speed, the sea wave spectrum peak frequency, the sea wave main wave direction and the sea wave energy of each grid point at the initial moment, and the grid point is the central point of the corresponding grid. According to the invention, the forecasting speed of the typhoon wave field data is improved.

Description

Technical Field

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

[0002] A tropical cyclone is a low-pressure vortex that occurs between the tropical and subtropical ocean surface and the atmospheric system. It is characterized by strong winds, low central pressure, and high-intensity rainfall. It is one of the most dangerous and destructive natural disaster phenomena. In order to ensure navigation safety and prevent coastal disasters, it is necessary to conduct a rapid and effective assessment of the wave field under the control of the typhoon.

[0003] In the northern hemisphere, due to the typhoon's translation and rotation characteristics, 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 cause of 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 forecast of the wave field under the action of the typhoon is a challenging task. At present, the third-generation wave model can simulate the spatial distribution of the wave field under the typhoon more reliably. However, due to the low resolution of the wind field data commonly used in the driving model due to the limitation of computational cost, it is difficult to correctly distinguish the rapidly changing wind field and the large local wind speed gradient in the typhoon. Therefore, it is of great significance to establish a reliable typhoon wave field model and realize rapid forecasting. 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 forecasting speed of typhoon wave field data.

[0005] To achieve the above objectives, this application provides the following solutions: In a first aspect, the present application provides a typhoon wave field simulation method, comprising: 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.

[0006] In a second aspect, the present application provides a typhoon wave field simulation device, the typhoon wave field simulation device comprising: 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 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.

[0007] In a third aspect, the present application provides 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 steps of any one of the above-mentioned typhoon wave field simulation methods.

[0008] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any one of the above-mentioned typhoon wave field simulation methods.

[0009] According to the specific embodiments provided in this application, this application discloses the following technical effects: The present application provides a typhoon wave field simulation method, device, equipment and medium, which adopts a polar coordinate grid to track the typhoon center of a target typhoon, wherein the typhoon center is the polar coordinate origin of the polar coordinate grid, and the polar coordinate grid encrypts the high wind speed area, thereby reducing the amount of calculation while highlighting the key areas, and improving the calculation speed. A two-dimensional wave ray theory parameterized model is adopted to construct a typhoon wave field model, and the typhoon wave field model requires fewer parameters and less calculation, which further improves the calculation speed, thereby improving the forecasting speed of typhoon wave field data. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0011] Figure 1 A schematic diagram of a flow chart of a typhoon wave field simulation method provided in an embodiment of the present application; Figure 2 A schematic diagram of the principle of a typhoon wave field simulation method provided in one embodiment of the present application; Figure 3 A schematic diagram of a polar coordinate grid provided in an embodiment of the present application; Figure 4 A schematic diagram of a typhoon wave field simulation effect provided in an embodiment of the present application; Figure 5 A schematic diagram showing the effect of verifying a typhoon wave field model using actual wave data measured by an altimeter provided in an embodiment of the present application; Figure 6 A schematic diagram of the effect of using the Holland wind farm to drive the SWAN wave model to verify the typhoon wave field model provided in an embodiment of the present application; Figure 7 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0012] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0013] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0014] This application provides a typhoon wave field simulation method, such as Figure 1 and Figure 2 As shown, the typhoon wave field simulation method includes steps 101 to 104.

[0015] Step 101: Track the typhoon center of the target typhoon using a polar coordinate grid, where the typhoon center is the polar coordinate origin of the polar coordinate grid.

[0016] Step 102: 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.

[0017] Step 103: Use the two-dimensional wave ray theory parameterized model to construct a typhoon wave field model.

[0018] Step 104: 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.

[0019] The present application adopts a polar coordinate grid to track the typhoon center of the target typhoon. The typhoon center is the polar coordinate origin of the polar coordinate grid. The polar coordinate grid will encrypt the high wind speed area, which not only highlights the key areas but also reduces the amount of calculation and improves the calculation speed. A two-dimensional wave ray theory parameterized model is adopted to construct a typhoon wave field model. The typhoon wave field model requires fewer parameters and less calculation, which further improves the calculation speed, thereby improving the forecast speed of typhoon wave field data.

[0020] During the movement of a typhoon, the range of typhoon waves is almost the same as the range of the typhoon wind field. 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 high wind speed areas, the use of polar coordinate grids will encrypt the area, highlighting key areas while reducing the amount of calculation.

[0021] In an exemplary embodiment, a typhoon wave field simulation method of the present application is described by taking Typhoon Lotus as an example. In the polar coordinate grid, the typhoon moving direction is the starting azimuth , 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, and the polar coordinate grid is as follows Figure 3 shown.

[0022] In an exemplary embodiment, although there are individual differences in the spatial distribution of the typhoon wind field, most of them can be approximately considered as translational vortices. The typhoon forcing field of the present application is the Holland typhoon wind field model (Holland wind field), which is a circularly symmetric wind field model based on the gradient wind equation. The expression of the Holland typhoon wind field model is: .

[0023] in, The radius from the typhoon center The gradient wind speed outside the atmospheric boundary layer, It is the ambient air pressure away from the typhoon, that is, the sea surface pressure 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: , dimensionless parameter defines the shape of the typhoon vortex, is the maximum wind speed radius of the target typhoon.

[0024] In an exemplary embodiment, to calculate the reference wind speed , gradient wind speed It is necessary to perform correction under the assumption of the roughness of the boundary layer profile surface, and determine the reference wind speed on the corresponding wave ray according to the typhoon forcing field, specifically including: using a correction formula to correct the gradient wind speed to obtain the reference wind speed, and the reference wind speed is the wind speed at a standard reference height of 10m.

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

[0026] The wind speed is spiraling toward the center of the typhoon, with an average inflow angle of ,in , , All are preset parameters.

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

[0028] The typhoon information provided by IBTrACS is once every 3 hours, so the Holland wind field that drives the model is also set to change once every 3 hours. During this period, the typhoon wind field only shifts in position, and the range and wind speed remain unchanged. The translation speed is calculated based on the typhoon position at two moments. The results of the Holland typhoon wind field model in the first 3 hours of operation are used as the initial field for the next 3 hours of operation.

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

[0030] The two-dimensional wave-ray theory parameterized model starts from the energy and momentum conservation equations and derives the self-similarity relationship between the wave spectral peak frequency, main wave direction and energy based on the wave growth relationship in a limited wind zone.

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

[0032] (1) , which represents the change of wave energy e.

[0033] (2) , which represents the change of the peak frequency of the spectrum. The nonlinear interaction of the waves reduces the peak frequency of the spectrum.

[0034] (3) , which reflects the wave direction changes.

[0035] (4) , which expresses the propagation of wave rays over time.

[0036] 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, indicating the change in spectrum peak frequency caused by wind energy input. represents the energy change caused by the vertical divergence of the group velocity, and are constants, is the peak wave number, It is the inverse of the wave age. and The step function involved is Indicates wind direction. Unit vector representing the direction of the group velocity.

[0037] Peak phase velocity =g / , the peak group velocity is 0.5 times the peak phase velocity, and g is the gravitational acceleration. Therefore, the two parameters of the peak group velocity and the peak frequency can be directly calculated.

[0038] The average group velocity and the peak group velocity can be obtained through the JONSWAP spectrum (Joint North Sea Wave Project Spectrum) The relationship between: , ; and represent the dimensionless forms of wind energy input and energy dissipation, respectively, , , dominates the evolution of wave energy, is the wind energy input in units, is the energy dissipated in units.

[0039] Since the wave spectrum generated by typhoon force conforms to the traditional limited wind zone JONSWAP spectrum, it has "spectral stability". and The coefficients can be derived using the JONSWAP spectrum.

[0040] In an exemplary embodiment, the initial field is In the initial field, the main wave direction of the ocean wave is the direction corresponding to the reference wind speed, and the ocean wave spectrum peak frequency and the ocean wave energy are obtained according to the wind-wave growth relationship in finite wind time.

[0041] 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: , .

[0042] 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. , , and All of them can be obtained from the wind-wave growth relationship in limited wind areas.

[0043] In one exemplary embodiment, during the differential equation solution, the shorter wind time , get the initial time The peak frequency and energy of the spectrum, the main wave direction of the sea wave at the initial moment is the same as the wind direction at that location. Since the typhoon wave field model is constructed using a two-dimensional wave ray theory parameterized model on a polar coordinate grid, the field result 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 a grid. The wave parameter of the wave ray with the largest wavelength is selected as the value at the grid point. The typhoon wave field model uses the fourth-order Runge-Kutta method to solve the differential equations on a polar coordinate grid, and the model time step is set to 300s.

[0044] Taking Typhoon Lotus (2009) as an example, the simulation results are as follows: Figure 4 As shown, Figure 4 Parts (a), (b), (c), (d), (e) and (f) are the typhoon wave field simulation effects at 0h, 6h, 12h, 18h, 24h and 30h respectively. Figure 4 middle is the effective wave height of typhoon waves. The typhoon wave field model fully demonstrates the spatial structure of the typhoon wave field and the development process of the waves. The wave field of the typhoon wave field model has the following characteristics: (1) It is structurally asymmetric. Due to the rotation and translation characteristics of the typhoon, the group velocity resonance phenomenon caused makes the typhoon wave field have obvious asymmetry. (2) The 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 waves, the wind and waves in this area have the characteristics of high energy and high group velocity, and continue to diverge outward to the low wind speed area to form surge waves; in the low wind speed area of ​​the typhoon, the surge waves propagating from the high wind speed area are mixed with the local wind waves. The surge waves, as the dominant waves, are the main components of the mixed wave field in this area. In addition, the typhoon wave field model was verified using the actual wave data measured by the altimeter, such as Figure 5 As shown, Figure 5 Part (a) shows the wind speed comparison effect. Figure 5 The middle part (b) shows the contrast effect of effective wave height. Figure 5 Part (c) is a plot of the correlation coefficient R between the satellite significant wave height and the typhoon wave field model. Overall, it is consistent with the altimeter data. The Holland wind field, which is the same as the typhoon wave field model, is used to drive the SWAN (Simulating WAves Nearshore) model. The results are as follows: Figure 6As shown, Figure 6 Part (a) shows the comparison between satellite data, two-dimensional wave rays and SWAN model, and part (b) shows the correlation coefficient between the effective wave height of satellite data and the effective wave height of SWAN model. The results of typhoon wave field model are similar to SWAN, which proves that the method has high credibility. At the same time, the amount of calculation is much lower than that of SWAN model, which shows that the typhoon wave field model can simulate the typhoon wave field quickly and effectively. The typhoon wave field model has a short running time and can be run on a personal computer. It can be used as a substitute for the third-generation wave spectrum model when computing resources are insufficient.

[0045] This application adopts the above-mentioned typhoon wave field simulation method based on two-dimensional wave ray theory. The parameters required for calculation are only the central air pressure, maximum wind speed radius and position information of the typhoon at each moment, so as to predict and simulate the time and space distribution of typhoon wave energy, main wave direction and spectral peak frequency.

[0046] In an exemplary embodiment, 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, specifically including: using the fourth-order Runge-Kutta method to solve the typhoon wave field model to determine the wave spectrum peak frequency, wave main wave direction and wave energy that change with time on each wave ray.

[0047] Based on the same inventive concept, the embodiment of the present application also provides a typhoon wave field simulation device for implementing the typhoon wave field simulation method involved above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above method, so the specific limitations in one or more typhoon wave field simulation device embodiments provided below can refer to the limitations of the typhoon wave field simulation method above, and will not be repeated here.

[0048] In an exemplary embodiment, the present application provides a typhoon wave field simulation device comprising: The module for determining the relationship between the typhoon center and the polar coordinate grid is used to track the typhoon center of the target typhoon using the polar coordinate grid, wherein the typhoon center is the polar coordinate origin of the polar coordinate grid.

[0049] The typhoon forcing field construction module is used 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.

[0050] The typhoon wave field model construction module is used to construct a typhoon wave field model using a two-dimensional wave ray theory parameterized model.

[0051] 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 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.

[0052] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 7 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. 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. 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 an external device. 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, a typhoon wave field simulation method is implemented.

[0053] Those skilled in the art will understand that Figure 7 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components. In an exemplary embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

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

[0055] 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 used 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 must comply with relevant regulations.

[0056] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and 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-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0057] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a data processing logic of a programmable logic device, etc., but is not limited thereto.

[0058] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0059] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this 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.

2. The typhoon wave field simulation method according to claim 1, characterized in that: 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.

3. 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.

4. The typhoon wave field simulation method according to claim 3, 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.

5. The typhoon wave field simulation method according to claim 1, characterized in that: 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. Unit vector representing the direction of the group velocity.

6. 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.

7. 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.

8. 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.

9. 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 7.

10. 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 7 is implemented.

Images