Wavefront water roll calculation method, device, equipment and storage medium

Abstract

The application provides a wave surface water roller calculation method, device, equipment and storage medium, including: acquiring topographic data, wave elements and flow field data of a nearshore area; calculating a water roller inclination angle according to the topographic data, wave elements and flow field data; calculating an energy transfer rate according to the topographic data, wave elements and flow field data; calculating water roller energy through a water roller evolution equation according to the wave elements, flow field data, water roller inclination angle and energy transfer rate; calculating water roller stress according to the water roller energy, and adding the water roller stress to a water dynamic model to simulate nearshore three-dimensional wave-induced current movement. The application does not need to artificially adjust or calibrate parameters in the water roller evolution equation, can calculate wave surface water roller energy and water roller stress according to the topographic data, wave elements and flow field data, and consider the influence of the wave surface water roller energy and water roller stress on nearshore water dynamics, so that the simulation accuracy of three-dimensional wave-induced current can be effectively improved.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of coastal dynamics, and particularly relates to a wave face water roll calculation method, device, equipment and storage medium. BACKGROUND

[0002] After wave breaking, a turbulent water body moving with the wave will be generated in the front end of the wave crest, which is called wave face water roll. Numerous studies have shown that water roll energy can have an important influence on nearshore wave-induced current movement. Accurate simulation and calculation of wave face water roll energy is a prerequisite and important basis for establishing a reasonable nearshore wave-induced current mathematical model.

[0003] Currently, there are two methods for calculating wave face water roll energy. One method is to calculate the water roll area by a simple empirical formula, and then to obtain the water roll energy. This method establishes a simple empirical relationship between the water roll area and the wave parameters, and it is difficult to reasonably reflect the generation, development and decline process of the water roll body after wave breaking. At the same time, the form and parameter values of the existing empirical formula are very different, and there is currently no widely applicable calculation formula. The other method is to establish a water roll evolution equation based on the conservation of water 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: the water roll inclination angle and the energy transfer rate. These two parameters have a great influence on the calculation results of water roll energy. Currently, most studies take these two parameters as model calibration parameters and take them as constants. On the one hand, this makes it difficult to meet the accuracy of the model at different locations, and on the other hand, when there is a lack of measured data, it is difficult to determine these two parameters, which restricts the universality of the water roll model.

[0004] In summary, there is currently a lack of a method and related device that is universally applicable and can accurately calculate wave face water roll. SUMMARY

[0005] Therefore, the present application aims to provide a wave face water roll calculation method, device, equipment and storage medium to solve the problems of low calculation accuracy and poor universality of the current water roll energy.

[0006] To achieve the above-mentioned purpose, the technical scheme of the present application is as follows:

[0007] In a first aspect, the present application provides a wave face water roll calculation method, comprising:

[0008] obtaining topographic data, wave elements and flow field data of a nearshore area;

[0009] calculating a water roll inclination angle based on the topographic data, wave elements and flow field data, wherein the water roll inclination angle is calculated based on the topographic water depth in the topographic data, the wave height and wave length in the wave elements, and the wave period average water level in the flow field data;

[0010] calculate the energy transfer rate according to the topographic data, the wave elements and the flow field data, wherein the energy transfer rate is calculated based on the topographic water depth and the topographic slope in the topographic data, the wave height and the wave length in the wave elements, and the wave period mean water level in the flow field data;

[0011] calculate the water roller energy according to the wave elements, the flow field data, the water roller inclination and the energy transfer rate, wherein the wave elements include the wave speed, the wave direction and the wave energy dissipation rate caused by wave breaking, and the flow field data include the wave period mean water level and the water depth mean flow speed;

[0012] add the water roller stress to the water dynamic model to simulate the nearshore three-dimensional wave-generated current motion according to the water roller energy.

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

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

[0015] a first calculation module configured to calculate a water roller inclination according to the topographic data, the wave elements and the flow field data, wherein the water roller inclination is calculated based on the topographic water depth in the topographic data, the wave height and the wave length in the wave elements, and the wave period mean water level in the flow field data;

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

[0017] a third calculation module configured to calculate a water roller energy according to the wave elements, the flow field data, the water roller inclination and the energy transfer rate, wherein the wave elements include the wave speed, the wave direction and the wave energy dissipation rate caused by wave breaking, and the flow field data include the wave period mean water level and the water depth mean flow speed;

[0018] a model output module configured to add the water roller stress to the water dynamic model to simulate the nearshore three-dimensional wave-generated current motion according to the water roller energy.

[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 capable of running on the processor, wherein the processor executes the program to implement the method according to the first aspect.

[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 for causing the computer to execute the method according to the first aspect.

[0021] Compared with the prior art, the wave face water roll calculation method, device, equipment and storage medium provided by the present application have the following beneficial effects:

[0022] The wave face water roll calculation method, device, equipment and storage medium provided by the present application do not need to artificially adjust or calibrate the parameters in the water roll evolution equation, can calculate the wave face water roll energy and water roll stress according to the terrain data, wave elements and flow field data, and consider the influence of the wave face water roll on the nearshore hydrodynamic force, thereby effectively improving the simulation accuracy of three-dimensional wave-generated flow. BRIEF DESCRIPTION OF DRAWINGS

[0023] The accompanying drawings, which form a part of this application, are included to provide a further understanding of the application and are incorporated in and constitute a part of this application. The embodiments of this application illustrated in the drawings, and their description, are presented to add generic scope to this application. In the drawings:

[0024] Figure 1 A wave face water roll calculation method flowchart according to an embodiment of the present application;

[0025] Figure 2 A comparison chart of the wave face water roll calculation method according to an embodiment of the present application and the wave-generated flow simulation results of the prior art based on the LIP11D experiment 1a group;

[0026] Figure 3 A comparison chart of the wave face water roll calculation method according to an embodiment of the present application and the wave-generated flow simulation results of the prior art based on the LIP11D experiment 1b group;

[0027] Figure 4 A comparison chart of the wave face water roll calculation method according to an embodiment of the present application and the wave-generated flow simulation results of the prior art based on the LIP11D experiment 1c group;

[0028] Figure 5 A wave face water roll calculation device structure schematic diagram according to an embodiment of the present application;

[0029] Figure 6 An electronic device hardware structure schematic diagram according to an embodiment of the present application. DETAILED DESCRIPTION

[0030] For the purposes of the present application, the technical solutions and advantages thereof are more clearly apparent, the following further describes the present application in detail with reference to specific examples and with reference to the accompanying drawings.

[0031] It should be noted that, unless otherwise defined, technical terms or scientific terms used in the embodiments of the present application should be understood as their common meanings to those having ordinary skills in the art to which the present application pertains. The terms "first", "second", and similar terms used in the embodiments of the present application do not denote any order, quantity, or importance, but are merely used to distinguish different components. The terms "include", "contain", and similar terms mean that the elements or objects before the terms encompass the elements or objects listed after the terms and their equivalents, and do not exclude other elements or objects. The terms "connect" or "connected" and similar terms do not mean physical or mechanical connections, but can include electrical connections, whether direct or indirect. The terms "upper", "lower", "left", "right", and the like merely denote relative positional relationships, which can change when the absolute positions of the described objects change.

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

[0033] Referring to FIG. 1, Figure 1 The embodiment provides a wave field calculation method, and specifically includes the following steps.

[0034] In step S101, topographic data, wave elements, and flow field data of a nearshore area are obtained.

[0035] Specifically, in the embodiment, the LIP11D flume experiment (the LIP11D flume experiment is a conventional experiment in the art, which will not be described further herein) is taken as an example to obtain topographic water depth and topographic slope data.

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

[0037] The general ocean model (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-generated flow experiment, so the FVCOM-SWAN coupled model is mainly used to simulate and calculate the wave-generated flow, and the wave period average water level and water depth average flow velocity of the nearshore area are obtained from the model results.

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

[0039] Specifically, in the present embodiment, the water roll inclination is calculated according to the terrain water depth, the wave height, the wave length and the wave period average water level obtained in step S101, and the calculation formula is as follows:

[0040] (1)

[0041] In the formula, the water roll inclination is represented by the root mean square wave height of irregular wave is represented by the wave length is represented by and the Ursell number is represented by

[0042] (2)

[0043] In the formula, the water depth containing water level change is represented by the calculation formula of which is as follows:

[0044] (3)

[0045] In the formula, the terrain water depth is represented by and the wave period average water level is represented by

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

[0047] Specifically, in the present embodiment, the energy transmission rate is calculated according to the terrain water depth, the terrain slope and the wave element obtained in step S101, and the calculation formula is as follows:

[0048] (4)

[0049] In the formula, the energy transmission rate is represented by the terrain slope is represented by and the wave number is represented by

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

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

[0052] (5)

[0053] In the formula, represents the water roller energy, represents the wave speed, represents the wave direction angle, and respectively represent , the water depth average flow velocity in the direction, represents time, represents the wave energy dissipation rate caused by wave breaking, represents the water roller energy dissipation rate, and the calculation formula is as follows:

[0054] (6)

[0055] In the formula, is the gravitational acceleration.

[0056] Step S105, the water roller stress is calculated according to the water roller energy, and the water roller stress is added to the water dynamic model to simulate the nearshore three-dimensional wave-induced current motion.

[0057] Specifically, in the embodiment, the water roller stress is further calculated according to the water roller energy calculated in step S104, and the calculation formula is as follows:

[0058] (7)

[0059] In the formula, , , and are the wave surface water roller stress terms; , are the components of the wave number in the , directions; is the vertical distribution function, and the calculation formula is as follows:

[0060] (8)

[0061] wherein, is the attenuation coefficient related to the vertical distribution, denotes the vertical coordinate.

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

[0063] (9)

[0064] (10)

[0065] wherein, , and are the velocity components in the , and directions, respectively; is the Coriolis force coefficient; and are the total density and the reference density, respectively; is the atmospheric pressure at the sea surface; is the vertical eddy viscosity coefficient; , are the horizontal momentum diffusion terms; , , , , and are the three-dimensional radiation stress terms; is the transformed vertical coordinate, and the formula for calculating the same is:

[0066] (11)

[0067] Specifically, it is added to the FVCOM model of the embodiment, so as to consider the effect of the wave face water roller on the water flow. In the embodiment, based on the FVCOM-SWAN coupling model, the wave face water roller calculation method proposed in the application is used to consider the influence of the wave face water roller on the nearshore hydrodynamic motion, and the wave-generated flow and its vertical structure are simulated and calculated.

[0068] Figures 2 to 4The wave-induced current calculated by the method of the embodiment is compared with the wave-induced current calculated by the constant water roller angle and energy transfer rate. As can be seen from the figure, the currently commonly used water roller angle and energy transfer rate value calculation method will produce obvious deviation in the calculation of nearshore wave-induced current at some positions, and the water roller calculation method proposed in the application effectively improves this, so that the model can more accurately simulate and calculate the nearshore three-dimensional wave-induced current structure.

[0069] The model accuracy is evaluated by using the statistical parameters of relative root mean square error (RRMSE) and Pearson correlation coefficient (r):

[0070] (12)

[0071] (13)

[0072] In the formula, and respectively represent the simulation and experimental measurement results of the model, and respectively represent the average values of the simulation and measurement results, and n is the sample number. The statistical parameter calculation results are shown in Table 1. It can be seen that the method described in the application provides more accurate simulation results for different experimental groups. Compared with the currently commonly used method, the RRMSE of the wave-induced current calculated by the application is obviously lower, and the correlation coefficient is closer to 1, and the application significantly improves the accuracy of the model in calculating the wave-induced current.

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

[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 roller calculation method proposed in the embodiment can reasonably describe the evolution of the nearshore wave surface water roller without adjusting parameters, thereby significantly improving the simulation accuracy of the nearshore wave-induced current mathematical model.

[0077] It should be noted that the above describes some embodiments of the application. Other embodiments are within the scope of the appended claims. In some cases, the acts or steps recited in the claims can be performed in an order other than that described above and still achieve desirable results. In addition, the processes depicted in the figures do not necessarily require the particular order shown or sequential order to achieve the desired results. In certain implementations, multitasking and parallel processing can be advantageous.

[0078] ​Based on the same inventive concept, embodiments of the present application also provide a wave-induced current calculation device corresponding to the method of any of the above embodiments.

[0079] As shown in Figure 5 the wave-induced current calculation device comprises:

[0080] The data acquisition module 11 is configured to acquire topographic data, wave element data and flow field data of the nearshore area.

[0081] The first calculation module 12 is configured to calculate a water roll inclination angle based on the topographic data, the wave element data and the flow field data, wherein the water roll inclination angle is calculated based on a topographic water depth in the topographic data, a wave height and a wave length in the wave element data, and a wave period average water level in the flow field data.

[0082] The second calculation module 13 is configured to calculate an energy transmission rate based on the topographic data, the wave element data and the flow field data, wherein the energy transmission rate is calculated based on a topographic water depth and a topographic slope in the topographic data, a wave height and a wave length in the wave element data, and a wave period average water level in the flow field data.

[0083] The third calculation module 14 is configured to calculate a water roll energy based on the wave element data, the flow field data, the water roll inclination angle and the energy transmission rate, and through a water roll evolution equation, wherein the wave element data comprises a wave speed, a wave direction and a wave energy dissipation rate caused by wave breaking, and the flow field data comprises a wave period average water level and a water depth average flow speed.

[0084] The model output module 15 is configured to calculate a water roll stress based on the water roll energy, and add the water roll stress to a hydrodynamic model to simulate a nearshore three-dimensional wave-induced current motion.

[0085] For the convenience of description, the above device is described in various modules based on functions. Of course, the functions of the modules can be implemented in one or more software and / or hardware when implementing the embodiments of the present application.

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

[0087] Based on the same inventive concept, embodiments of the present application also provide an electronic device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the method of any of the above embodiments when executing the program.

[0088] Figure 6A more specific electronic device hardware structure schematic diagram provided by the embodiment is shown, and the device can 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 through the bus 1050 for communication between each other inside 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, etc., for executing related programs to implement the technical solutions provided by the embodiments of the present specification.

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

[0091] The input / output interface 1030 is used to connect input / output modules to realize information input and output. The input / output modules can be configured as components in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. The input device can include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device can 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 in the figure) to realize the communication interaction between the device and other devices. The communication module can realize communication through a wired manner (such as USB, network cable, etc.) or through a wireless manner (such as mobile network, WIFI, Bluetooth, etc.).

[0093] The bus 1050 includes a channel to transmit information between various components (such as the processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040) of the device.

[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 can also include other components necessary for normal operation. In addition, those skilled in the art can understand that the above device can also only contain the components necessary to implement the embodiments of the present application, and does not necessarily contain 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 preceding embodiments, and has the beneficial effects of the corresponding method embodiments, which are not described here.

[0096] Based on the same inventive concept, corresponding to the method of any of the above embodiments, the present application also provides a non-transitory computer readable storage medium storing computer instructions for causing the computer to execute the method of any of the above embodiments.

[0097] The computer readable medium of the present embodiment includes permanent and non-permanent, removable and non-removable media, which can be realized by any method or technology to store information. The information can be computer readable instructions, data structures, program modules 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 technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic tape disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device.

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

[0099] Those skilled in the art should understand that the discussion of any of the above embodiments is only exemplary and is not intended to imply that the scope (including claims) of the present application is limited to these examples; under the idea of the present application, the above embodiments or technical features in different embodiments can also be combined, the steps can be implemented in any order, and there are many other changes of different aspects of the embodiments of the present application as described above. In order to be brief, they are not provided in detail.

[0100] Additionally, to simplify the description and discussion, and so as not to obscure the embodiments of the application being presented, the well-known functions or constructions of integrated circuit (IC) chips and other components can or can not be shown in the figures and will be omitted as not to unnecessarily obscure the embodiments of the application being presented. Moreover, the devices can be shown in block diagram form in order to avoid obscuring the embodiments of the application, and this also acknowledges the fact that the details in regard to the implementation of the block diagram devices are highly dependent on the platform within which the embodiments of the application are to be implemented (i.e., these details should be well within the purview of one of ordinary skill in the art). Where specific details are set forth in order to describe an illustrative embodiment of the application, it will be apparent to one of ordinary skill in the art that the embodiments of the application can be practiced without, or with variation of, these specific details. Thus, the description is to be considered as illustrative and not restrictive, and the scope of the application should be determined not with reference to the above description, but should be given to the appended claims.

[0101] While the application has been described in connection with specific embodiments thereof, many alternatives, modifications and variations will be apparent to those skilled in the art in light of the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) can use the embodiments discussed.

[0102] Embodiments of the application are intended to cover all such alternatives, modifications and variations as falling within the scope of the appended claims. Accordingly, any omission, modification, equivalent replacement, improvement, etc. made in the spirit and principle of the embodiments of the application should be included in the scope of protection of the application.

Claims

1. A wave-surface water rolling calculation method, characterized by, The method comprises: obtaining topographic data, wave elements and flow field data of a nearshore area; calculating a water roller inclination angle based on the topographic data, the wave elements and the flow field data, wherein the water roller inclination angle is calculated based on a topographic water depth in the topographic data, a wave height and a wave length in the wave elements, and a wave period average water level in the flow field data; the water roller inclination angle calculation formula is: ； wherein , ; wherein denotes the water roll angle, denotes the root mean square wave height of an irregular wave, denotes the wavelength, denotes the Ursell number, denotes the water depth including water level variation, denotes the topographic water depth, denotes the wave period mean water level; calculating an energy transfer rate based on the topographic data, the wave elements and the flow field data, wherein the energy transfer rate is calculated based on a topographic water depth and a topographic slope in the topographic data, a wave height and a wave length in the wave elements, and a wave period average water level in the flow field data; the energy transfer rate calculation formula is: ； wherein represents an energy transfer rate, represents a terrain slope, represents a wave number, represents a root mean square wave height of an irregular wave, represents a wavelength, represents a water depth including a water level change; calculating water roller energy based on the wave elements, the flow field data, the water roller inclination angle and the energy transfer rate through a water roller evolution equation, wherein the wave elements include a wave speed, a wave direction and a wave energy dissipation rate caused by wave breaking, and the flow field data includes a wave period average water level and a water depth average flow speed; calculating water roller stress based on the water roller energy, and adding the water roller stress to a water dynamic model to simulate nearshore three-dimensional wave-induced current movement.

2. The method of claim 1, wherein, the water roller energy calculation formula is: ； wherein, represents the water roller energy, represents the wave speed, represents the wave direction angle, and respectively represent , the average flow velocity in the direction of the water depth, represents time, represents the energy transfer rate, represents the wave energy dissipation rate caused by wave breaking, represents the water roller energy dissipation rate, and the calculation formula is: ； wherein is the gravitational acceleration, denotes the water roll angle.

3. The method of claim 1, wherein, the water roller stress calculation formula is: ； wherein represents the water roll energy, , , and represents the wave face water roll stress term, , respectively represent the components of the wave number in the , directions, represents the vertical distribution function, which is calculated as ； wherein denotes the attenuation coefficient related to the vertical distribution, denotes the topographic water depth, denotes the wave period mean water level, denotes the vertical coordinate.

4. The method of claim 3, wherein, adding a water roller stress term to a momentum equation of the water dynamic model to obtain a wave-induced current result, and the momentum equation formula is: ； ； wherein , ; where, , and represent velocity components in the directions of , and , represents the Coriolis force coefficient, and represent total density and reference density, respectively, represents atmospheric pressure at the sea surface, represents the vertical eddy viscosity coefficient, , represent horizontal momentum diffusion terms, , , , , and represent three-dimensional radiation stress terms, represents the transformed vertical coordinate, represents the water depth including the water level variation, represents the topographic water depth, represents the wave period-averaged water level.

5. The method of claim 1, wherein, The method further comprises: evaluating the water dynamic model accuracy through a relative root mean square error and a Pearson correlation coefficient, wherein the relative root mean square error formula is: ； ； wherein and denote the results of the model simulation and experimental measurement, respectively, and denote the average of the simulation and measurement results, respectively, and n denotes the number of samples.

6. A wave-surface water rolling calculation device characterized by comprising: The method comprises: a data acquisition module configured to obtain topographic data, wave elements and flow field data of a nearshore area; a first calculation module configured to calculate a water roller inclination angle based on the topographic data, the wave elements and the flow field data, wherein the water roller inclination angle is calculated based on a topographic water depth in the topographic data, a wave height and a wave length in the wave elements, and a wave period average water level in the flow field data; the water roller inclination angle calculation formula is: ； wherein , ; wherein denotes the water roll angle, denotes the root mean square wave height of irregular waves, denotes the wavelength, denotes the Ursell number, denotes the water depth including water level variations, denotes the topographic water depth, denotes the wave period mean water level; a second calculation module configured to calculate an energy transfer rate based on the topographic data, the wave elements and the flow field data, wherein the energy transfer rate is calculated based on a topographic water depth and a topographic slope in the topographic data, a wave height and a wave length in the wave elements, and a wave period average water level in the flow field data; the energy transfer rate calculation formula is: ； wherein represents the energy transfer rate, represents the terrain slope, represents the wave number, represents the root mean square wave height of irregular waves, represents the wavelength, represents the water depth including water level variation; a third calculation module configured to calculate water roller energy based on the wave elements, the flow field data, the water roller inclination angle and the energy transfer rate through a water roller evolution equation, wherein the wave elements include a wave speed, a wave direction and a wave energy dissipation rate caused by wave breaking, and the flow field data includes a wave period average water level and a water depth average flow speed; a model output module configured to calculate water roller stress based on the water roller energy, and add the water roller stress to a water dynamic model to simulate nearshore three-dimensional wave-induced current movement.

7. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the method of any one of claims 1-5 when executing the program.

8. A non-transitory computer-readable storage medium, comprising: wherein, The non-transitory computer readable storage medium stores computer instructions for causing a computer to perform the method of any one of claims 1-5.

Images