Wave surface water rolling calculation method, device and equipment and storage medium

Abstract

The invention provides a wave surface water rolling calculation method and device, equipment and a storage medium. The wave surface water rolling calculation method comprises the steps that terrain data, wave elements and flow field data of a near-shore area are acquired; calculating a water roll inclination angle according to the topographic data, the wave elements and the flow field data; according to the topographic data, the wave elements and the flow field data, calculating to obtain an energy transmissibility; according to the wave elements, the flow field data, the water roll inclination angle and the energy transmission rate, water roll energy is calculated through a water roll evolution equation; and calculating water rolling stress according to the water rolling energy, and adding the water rolling stress into the hydrodynamic model to simulate the near-shore three-dimensional wave-generated flow motion. According to the method, parameters in the water rolling evolution equation do not need to be manually adjusted or calibrated, the wave surface water rolling energy and the water rolling stress can be calculated according to the topographic data, the wave elements and the flow field data, the influence of the wave surface water rolling energy and the water rolling stress on coastal water power is considered, and the simulation precision of the three-dimensional wave-generated flow can be effectively improved.

Description

Technical Field

[0001] The present application belongs to the field of coastal dynamics technology, and in particular, relates to a wave surface water roll calculation method, device, equipment and storage medium. Background Art

[0002] When waves break, turbulent water will be generated at the front of the wave crest, which moves with the waves. This is called wave roll. Many studies have shown that the energy of water roll has an important impact on the movement of nearshore wave-induced currents. Accurate simulation and calculation of the energy of wave roll is a prerequisite and important basis for establishing a reasonable mathematical model of nearshore wave-induced currents.

[0003] At present, there are two main methods for calculating the water roll energy of wave surface. One method is to calculate the water roll area through a simple empirical formula, and then obtain the water roll energy. This method establishes a simple empirical relationship between the water roll area and the wave parameters, which is difficult to reasonably reflect the generation, development and decay process of the water roll body after the wave breaks; at the same time, the form and parameter values ​​of the existing empirical formulas are also very different, and there is currently no widely applicable calculation formula. Another method is to establish a water roll evolution equation based on the conservation of water body energy. This method can more reasonably reflect the change process of water roll energy after wave breaking. However, there are two parameters in the water roll evolution equation that need to be determined before calculation: water roll inclination angle and energy transfer rate. These two parameters have a great influence on the calculation results of water roll energy. Most of the current studies use these two parameters as model calibration parameters and take them as constants. On the one hand, it is difficult to meet the accuracy of the model at different locations. At the same time, when there is a lack of measured data, these two parameters are difficult to determine, which restricts the versatility of the water roll model.

[0004] In summary, there is currently a lack of universal methods and related devices that can accurately calculate wave surface water roll. Summary of the invention

[0005] In view of this, the present application aims to propose a wave surface water roll calculation method, device, equipment and storage medium to solve the current problems of low accuracy and poor versatility in water roll energy calculation.

[0006] To achieve the above purpose, the technical solution of this application is implemented as follows:

[0007] In a first aspect, the present application provides a method for calculating wave surface water roll, comprising:

[0008] Obtain terrain data, wave elements and flow field data in nearshore areas;

[0009] The water roll angle is calculated based on the terrain data, wave elements and flow field data, wherein the water roll angle is calculated based on the terrain water depth in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data;

[0010] The energy transfer rate is calculated based on the terrain data, wave elements and flow field data, wherein the energy transfer rate is calculated based on the terrain water depth and terrain slope in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data;

[0011] The water roll energy is calculated according to the wave elements, flow field data, water roll inclination angle and energy transfer rate through the water roll evolution equation, wherein the wave elements include wave speed, wave direction and wave energy dissipation rate caused by wave breaking, and the flow field data includes the average water level of the wave period and the average flow velocity of the water depth;

[0012] The water roll stress is calculated based on the water roll energy, and the water roll stress is added to the hydrodynamic model to simulate the three-dimensional wave-induced current movement near the shore.

[0013] In the second aspect, based on the same inventive concept, the present application also provides a wave surface water roll calculation device, comprising:

[0014] A data acquisition module is configured to acquire terrain data, wave elements and flow field data of a nearshore area;

[0015] A first calculation module is configured to calculate a water roll angle according to the terrain data, wave elements and flow field data, wherein the water roll angle is calculated based on the terrain water depth in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data;

[0016] a second calculation module, configured to calculate an energy transfer rate according to the terrain data, the wave elements and the flow field data, wherein the energy transfer rate is calculated based on the terrain water depth and terrain slope in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data;

[0017] A third calculation module is configured to calculate water roll energy through a water roll evolution equation according to the wave elements, flow field data, water roll inclination angle and energy transfer rate, wherein the wave elements include wave speed, wave direction and wave energy dissipation rate caused by wave breaking, and the flow field data includes the average water level of the wave period and the average flow velocity of the water depth;

[0018] The model output module is configured to calculate the water roll stress according to the water roll energy and add the water roll stress to the hydrodynamic model to simulate the three-dimensional wave-induced current movement near the shore.

[0019] In a third aspect, based on the same inventive concept, the present application also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method described in the first aspect is implemented.

[0020] In a fourth aspect, based on the same inventive concept, the present application also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method described in the first aspect.

[0021] Compared with the prior art, the wave surface water roll calculation method, device, equipment and storage medium described in this application have the following beneficial effects:

[0022] The wave surface water roll calculation method, device, equipment and storage medium described in the present application do not require manual adjustment or calibration of parameters in the water roll evolution equation. The wave surface water roll energy and water roll stress can be calculated based on terrain data, wave elements and flow field data, and their influence on nearshore hydrodynamics can be taken into account, which can effectively improve the simulation accuracy of three-dimensional wave-induced flow. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] The drawings constituting a part of the present application are used to provide a further understanding of the present application. The illustrative embodiments and descriptions of the present application are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0024] Figure 1 This is a flow chart of a method for calculating wave surface water roll described in an embodiment of the present application;

[0025] Figure 2 It is a comparison diagram of the wave surface water roll calculation method described in the embodiment of the present application and the wave-induced flow simulation results of the prior art based on the LIP11D experiment group 1a;

[0026] Figure 3 It is a comparison diagram of the wave surface water roll calculation method described in the embodiment of the present application and the wave-induced flow simulation results of the prior art based on the LIP11D experiment group 1b;

[0027] Figure 4 It is a comparison diagram of the wave surface water roll calculation method described in the embodiment of the present application and the wave-induced flow simulation results of the prior art based on the LIP11D experiment group 1c;

[0028] Figure 5 This is a schematic diagram of the structure of a wave surface water roll calculation device described in an embodiment of the present application;

[0029] Figure 6 This is a schematic diagram of the hardware structure of the electronic device described in the embodiment of the present application. DETAILED DESCRIPTION

[0030] In order to make the objectives, technical solutions and advantages of the present application more clearly understood, the present application is further described in detail below in combination with specific embodiments and with reference to the accompanying drawings.

[0031] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present application should be the usual meanings understood by people with ordinary skills in the field to which the present application belongs. The "first", "second" and similar words used in the embodiments of the present application do not represent any order, quantity or importance, but are only used to distinguish different components. "Including" or "comprising" and similar words mean that the elements or objects appearing in front of the word cover the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0032] The embodiments of the present application are described in detail below with reference to the accompanying drawings.

[0033] See also Figure 1 As shown, this embodiment provides a method for calculating wave surface water roll, which specifically includes the following steps:

[0034] Step S101, obtaining terrain data, wave elements and flow field data of the nearshore area.

[0035] Specifically, in this embodiment, the LIP11D water tank experiment (the LIP11D water tank experiment is a conventional experiment in the art and will not be further described here) is used as an example to obtain terrain water depth and terrain slope data.

[0036] The simulating waves nearshore (SWAN) model is used to simulate the wave field of the LIP11D tank experiment. Wave elements such as nearshore wave height, wavelength, wave direction, wave speed, and wave energy dissipation rate caused by wave breaking are obtained from the model results.

[0037] The finite-volume community ocean model (FVCOM) based on the finite volume method is used to simulate the flow field. The LIP11D flume experiment is a wave-induced flow experiment, so this embodiment mainly uses the FVCOM-SWAN coupling model to simulate and calculate the wave-induced flow, and obtains the nearshore wave period average water level and depth average flow velocity from the model results.

[0038] Step S102, calculating the water roll angle according to the terrain data, wave elements and flow field data, wherein the water roll angle is calculated based on the terrain water depth in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data.

[0039] Specifically, in this embodiment, the water roll angle is calculated according to the terrain water depth, wave height, wavelength and wave period average water level obtained in step S101, and the calculation formula is:

[0040] (1)

[0041] In the formula, represents the water roll angle, represents the RMS wave height of the irregular wave, represents the wavelength, represents the Ursell number, which is calculated as:

[0042] (2)

[0043] In the formula, It represents the water depth including water level changes, and its calculation formula is:

[0044] (3)

[0045] In the formula, Indicates topographic water depth; Indicates the average water level over a wave period.

[0046] Step S103, calculating the energy transfer rate according to the terrain data, wave elements and flow field data, wherein the energy transfer rate is calculated based on the terrain water depth and terrain slope in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data.

[0047] Specifically, in this embodiment, the energy transfer rate is calculated according to the terrain water depth, terrain slope and wave elements obtained in step S101, and the calculation formula is:

[0048] (4)

[0049] In the formula, represents the energy transfer rate, represents the terrain slope, Indicates the wave number.

[0050] Step S104, according to the wave elements, flow field data, water roll inclination angle and energy transfer rate, the water roll energy is calculated by the water roll evolution equation, wherein the wave elements include wave speed, wave direction and wave energy dissipation rate caused by wave breaking, and the flow field data includes the average water level of the wave period and the average flow velocity of the water depth.

[0051] Specifically, in this embodiment, according to the wave speed, wave direction and wave energy dissipation rate caused by wave breaking obtained in step S101, the average water level and average flow velocity of the water depth during the wave period, the water roll inclination angle calculated in step S102, and the energy transfer rate calculated in step S103, the water roll energy is calculated by the water roll evolution equation, and the calculation formula is:

[0052] (5)

[0053] In the formula, Indicates the boiling energy of water. represents the wave speed, represents the wave direction angle, and Respectively , The average velocity of water depth in the direction, Indicates time, represents the wave energy dissipation rate caused by wave breaking, It represents the water rolling energy dissipation rate, and its calculation formula is:

[0054] (6)

[0055] In the formula, is the acceleration due to gravity.

[0056] Step S105: Calculate the water rolling stress based on the water rolling energy, and add the water rolling stress to the hydrodynamic model to simulate the three-dimensional wave-induced current movement near the shore.

[0057] Specifically, in this embodiment, the water rolling stress is further calculated based on the water rolling energy calculated in step S104, and the calculation formula is:

[0058] (7)

[0059] In the formula, , , and is the wave surface water rolling stress term; , The wave number is , Directional weight; is the vertical distribution function, and its calculation formula is:

[0060] (8)

[0061] In the formula, is the attenuation coefficient related to the vertical distribution, Represents the vertical coordinate.

[0062] The water rolling stress term calculated by formula (7) is , , and It can be added to the momentum equation of the hydrodynamic model, which is:

[0063] (9)

[0064] (10)

[0065] In the formula, , and They are , and Directional velocity component; is the Coriolis force coefficient; and are the total density and reference density, respectively; is the atmospheric pressure at the sea surface; is the vertical eddy viscosity coefficient; , is the horizontal momentum diffusion term; , , , , and is the three-dimensional radiation stress term; is the transformed vertical coordinate, and its calculation formula is:

[0066] (11)

[0067] Specifically, it is added to the FVCOM model of this embodiment to consider the effect of wave surface water roll on water flow. In this embodiment, based on the FVCOM-SWAN coupling model, the wave surface water roll calculation method proposed in this application is adopted to consider the effect of wave surface water roll on nearshore hydrodynamic motion, and simulate and calculate wave-induced flow and its vertical structure.

[0068] Figures 2 to 4The results of wave-induced flow calculated by the method described in this embodiment are compared with those calculated by constant water roll angle and energy transfer rate. As can be seen from the figure, the commonly used water roll angle and energy transfer rate value determination method will produce obvious deviations in some positions when calculating nearshore wave-induced flow. The water roll calculation method proposed in this application makes effective improvements to this, so that the model can more accurately simulate and calculate the nearshore three-dimensional wave-induced flow structure.

[0069] The statistical parameters relative root mean square error (RRMSE) and Pearson correlation coefficient (r) were used to evaluate the model accuracy:

[0070] (12)

[0071] (13)

[0072] In the formula, and represent the results of model simulation and experimental measurement, respectively. and represent the average values ​​of simulation and measurement results respectively, and n is the number of samples.

[0073] The calculation results of statistical parameters are shown in Table 1. It can be seen that for different experimental groups, the method of the present invention provides more accurate simulation results. Compared with the currently commonly used methods, the RRMSE of the wave-induced flow results calculated by the present invention is significantly lower, and the correlation coefficient is closer to 1. The present invention significantly improves the accuracy of the model in calculating wave-induced flow.

[0074] Table 1 Relative root mean square error (RRMSE) and Pearson correlation coefficient (r) of different methods

[0075]

[0076] The above results show that the wave surface water roll calculation method proposed in this embodiment can reasonably describe the evolution of nearshore wave surface water roll without considering adjustment parameters, thereby significantly improving the simulation accuracy of the nearshore wave-induced flow mathematical model.

[0077] It should be noted that the above describes some embodiments of the present application. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the above embodiments and still achieve the desired results. In addition, the processes depicted in the accompanying drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0078] Based on the same inventive concept, corresponding to any of the above-mentioned embodiments and methods, an embodiment of the present application further provides a wave surface water roll calculation device.

[0079] like Figure 5 As shown, the wave surface water roll calculation device comprises:

[0080] The data acquisition module 11 is configured to acquire terrain data, wave elements and flow field data of the nearshore area;

[0081] A first calculation module 12 is configured to calculate a water roll angle according to terrain data, wave elements and flow field data, wherein the water roll angle is calculated based on the terrain water depth in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data;

[0082] The second calculation module 13 is configured to calculate the energy transfer rate according to the terrain data, the wave elements and the flow field data, wherein the energy transfer rate is calculated based on the terrain water depth and terrain slope in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data;

[0083] The third calculation module 14 is configured to calculate the water roll energy according to the wave elements, flow field data, water roll inclination angle and energy transfer rate through the water roll evolution equation, wherein the wave elements include wave speed, wave direction and wave energy dissipation rate caused by wave breaking, and the flow field data includes the average water level of the wave period and the average flow velocity of the water depth;

[0084] The model output module 15 is configured to calculate the water roll stress according to the water roll energy, and add the water roll stress to the hydrodynamic model to simulate the three-dimensional wave-induced current movement near the shore.

[0085] For the convenience of description, the above devices are described in terms of functions and are divided into various modules. Of course, when implementing the embodiments of the present application, the functions of each module can be implemented in the same or multiple software and / or hardware.

[0086] The device of the above embodiment is used to implement the corresponding method in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be described in detail here.

[0087] Based on the same inventive concept, corresponding to any of the above-mentioned embodiments and methods, an embodiment of the present application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method described in any of the above embodiments is implemented.

[0088] Figure 6A more specific schematic diagram of the hardware structure of an electronic device provided in this embodiment is shown, and the device may include: a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. The processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040 are connected to each other through the bus 1050 in the device.

[0089] The processor 1010 can be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0090] The memory 1020 may be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 1020 may store an operating system and other application programs. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program codes are stored in the memory 1020 and are called and executed by the processor 1010.

[0091] The input / output interface 1030 is used to connect the input / output module to realize information input and output. The input / output module can be configured in the device as a component (not shown in the figure), or it can be externally connected to the device to provide corresponding functions. The input device may include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device may include a display, a speaker, a vibrator, an indicator light, etc.

[0092] The communication interface 1040 is used to connect a communication module (not shown) to realize communication interaction between the device and other devices. The communication module can realize communication through wired means (such as USB, network cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0093] The bus 1050 includes a path that transmits information between the various components of the device (eg, the processor 1010 , the memory 1020 , the input / output interface 1030 , and the communication interface 1040 ).

[0094] It should be noted that, although the above device only shows the processor 1010, the memory 1020, the input / output interface 1030, the communication interface 1040 and the bus 1050, in the specific implementation process, the device may also include other components necessary for normal operation. In addition, it can be understood by those skilled in the art that the above device may also only include the components necessary for implementing the embodiments of the present specification, and does not necessarily include all the components shown in the figure.

[0095] The electronic device of the above embodiment is used to implement the corresponding method in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be described in detail here.

[0096] Based on the same inventive concept, corresponding to any of the above-mentioned embodiments, the present application also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method described in any of the above embodiments.

[0097] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, read-only compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, tape disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device.

[0098] The computer instructions stored in the storage medium of the above embodiments are used to enable the computer to execute the method described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0099] A person skilled in the art should understand that the discussion of any of the above embodiments is merely illustrative and is not intended to imply that the scope of the present application (including the claims) is limited to these examples. In line with the concept of the present application, the technical features in the above embodiments or different embodiments may be combined, the steps may be implemented in any order, and there are many other variations of the different aspects of the embodiments of the present application as described above, which are not provided in detail for the sake of simplicity.

[0100] In addition, to simplify the description and discussion, and in order not to make the embodiments of the present application difficult to understand, the known power / ground connections to the integrated circuit (IC) chip and other components may or may not be shown in the provided drawings. In addition, the device may be shown in the form of a block diagram to avoid making the embodiments of the present application difficult to understand, and this also takes into account the fact that the details of the implementation of these block diagram devices are highly dependent on the platform on which the embodiments of the present application are to be implemented (that is, these details should be fully within the scope of understanding of those skilled in the art). Where specific details (e.g., circuits) are set forth to describe exemplary embodiments of the present application, it is obvious to those skilled in the art that the embodiments of the present application can be implemented without these specific details or with changes in these specific details. Therefore, these descriptions should be considered illustrative rather than restrictive.

[0101] Although the present application has been described in conjunction with specific embodiments of the present application, many alternatives, modifications and variations of these embodiments will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may use the discussed embodiments.

[0102] The embodiments of the present application are intended to cover all such substitutions, modifications and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of the present application should be included in the scope of protection of the present application.

Claims

1. A method for calculating wave surface water roll, characterized in that: include: Obtain terrain data, wave elements and flow field data in nearshore areas; The water roll angle is calculated based on the terrain data, wave elements and flow field data, wherein the water roll angle is calculated based on the terrain water depth in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data; The energy transfer rate is calculated based on the terrain data, wave elements and flow field data, wherein the energy transfer rate is calculated based on the terrain water depth and terrain slope in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data; The water roll energy is calculated according to the wave elements, flow field data, water roll inclination angle and energy transfer rate through the water roll evolution equation, wherein the wave elements include wave speed, wave direction and wave energy dissipation rate caused by wave breaking, and the flow field data includes the average water level of the wave period and the average flow velocity of the water depth; The water roll stress is calculated based on the water roll energy, and the water roll stress is added to the hydrodynamic model to simulate the three-dimensional wave-induced current movement near the shore.

2. The method according to claim 1, characterized in that: The water roll angle calculation formula is: ； in, , ; In the formula, represents the water roll angle, represents the RMS wave height of the irregular wave, represents the wavelength, Represents Ursell number, Indicates the water depth including water level changes, Indicates the terrain water depth, Represents the average water level during the wave period.

3. The method according to claim 1, characterized in that The energy transfer rate calculation formula is: ； In the formula, represents the energy transfer rate, represents the terrain slope, represents the wave number, represents the RMS wave height of the irregular wave, represents the wavelength, Indicates the water depth including water level changes.

4. The method according to claim 1, characterized in that: The water boiling energy calculation formula is: ； In the formula, Indicates the water boiling energy, represents the wave speed, represents the wave direction angle, and Respectively , The average velocity of water depth in the direction, Indicates time, represents the energy transfer rate, represents the wave energy dissipation rate caused by wave breaking, It represents the water rolling energy dissipation rate, and its calculation formula is: ； In the formula, is the acceleration due to gravity, Indicates the water roll angle.

5. The method according to claim 1, characterized in that The water rolling stress calculation formula is: ； In the formula, Indicates the water boiling energy, , , and represents the wave surface water rolling stress term, , Respectively represent the wave number in , The weight in direction, represents the vertical distribution function, and its calculation formula is: ； In the formula, represents the attenuation coefficient related to the vertical distribution, Indicates the terrain water depth, represents the average water level during the wave period, Represents the vertical coordinate.

6. The method according to claim 5, characterized in that The water rolling stress term is added to the momentum equation of the hydrodynamic model to obtain the wave-induced flow result. The momentum equation formula is: ； ； in, , ; In the formula, , and Respectively , and The velocity component in the direction, represents the Coriolis force coefficient, and denote the total density and reference density respectively, is the atmospheric pressure at the sea surface, represents the vertical eddy viscosity coefficient, , represents the horizontal momentum diffusion term, , , , , and represents the three-dimensional radiation stress term, represents the transformed vertical coordinate, Indicates the water depth including water level changes, Indicates the terrain water depth, Represents the average water level during the wave period.

7. The method according to claim 1, characterized in that Also includes: The accuracy of the hydrodynamic model is evaluated by the relative root mean square error and the Pearson correlation coefficient, wherein the relative root mean square error formula is: ； ； In the formula, and represent the results of model simulation and experimental measurement, respectively. and They represent the average values ​​of simulation and measurement results respectively, and n represents the number of samples.

8. A wave surface water roll calculation device, characterized in that: include: A data acquisition module is configured to acquire terrain data, wave elements and flow field data of a nearshore area; A first calculation module is configured to calculate a water roll angle according to the terrain data, wave elements and flow field data, wherein the water roll angle is calculated based on the terrain water depth in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data; a second calculation module, configured to calculate an energy transfer rate according to the terrain data, the wave elements and the flow field data, wherein the energy transfer rate is calculated based on the terrain water depth and terrain slope in the terrain data, the wave height and wavelength in the wave elements, and the wave period average water level in the flow field data; A third calculation module is configured to calculate water roll energy through a water roll evolution equation according to the wave elements, flow field data, water roll inclination angle and energy transfer rate, wherein the wave elements include wave speed, wave direction and wave energy dissipation rate caused by wave breaking, and the flow field data includes the average water level of the wave period and the average flow velocity of the water depth; The model output module is configured to calculate the water roll stress according to the water roll energy and add the water roll stress to the hydrodynamic model to simulate the three-dimensional wave-induced current movement near the shore.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the method according to any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium, characterized in that: in, The non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to cause a computer to execute the method according to any one of claims 1 to 7.

Images