Sea area wave surface height prediction method after intervention of wave energy device and related device

Abstract

The invention discloses a sea area wave surface height prediction method after intervention of a wave energy device and a related device, and relates to the field of ocean energy utilization, and the method comprises the steps: building a numerical pool and setting parameters such as boundary conditions based on the main body size of the wave energy device, and carrying out the comparative analysis of an ideal wave numerical pool and a theoretical wave sine curve; importing the modeled wave energy device, setting a dynamic domain, a static domain and a boundary layer grid, calculating the conversion efficiency of the wave energy device, comparing the conversion efficiency with the conversion efficiency of the device obtained by the test until a relative error range is met, and calculating the wave energy device based on the final verified wave pool numerical model containing the device. Wave surface average change data in a period is obtained, a wave surface height distribution function variable factor is determined through a fitting method, wave surface change functions before and after the device are obtained, and then the wave surface heights before and after the device can be rapidly obtained. According to the method, the wave height change of the surrounding sea area after the wave energy device is arranged can be quickly calculated.

Description

TECHNICAL FIELD

[0001] The application relates to the technical field of ocean energy utilization, in particular to a method for predicting wave height in a sea area after intervention of a wave energy device and a related device. BACKGROUND

[0002] After the wave energy device is intervened in the sea area, the energy in the sea area is extracted, and the original wave power propagation process is changed. The wave height is lifted in front of the device, and attenuation occurs behind the device, that is, the "tail wave" phenomenon occurs, and the input energy of the rear device changes. With large-scale deployment of wave energy devices in the sea area, the overall wave field of the sea area changes compared to the original state, which will sharply affect the overall wave energy resource distribution. However, the existing wave energy resource evaluation methods do not consider the changes in the wave field caused by the large-scale deployment of wave energy devices, and the changes in resources caused by the changes in the wave surface are not considered. At the same time, since the wave surface propagation is always changing with the period, how to establish a model to convert the instantaneous changes at each time into a unified average change is also a challenge. If the traditional numerical simulation method is used, complete modeling of all devices is required, which not only consumes a large amount of computing resources, but also takes too long to calculate, making it difficult to obtain the wave field height changes in time, and unable to quickly obtain the changes in the overall wave energy resources, and also limiting the rapid optimization ability of the device array arrangement.

[0003] In summary, there is an urgent need to develop a method for predicting the wave surface and the tail wave surface after the wave energy device is intervened in the sea area. SUMMARY

[0004] The purpose of the present application is to provide a method for predicting the wave height in the sea area after intervention of a wave energy device, which can quickly calculate the wave height changes around the wave energy device.

[0005] To achieve the above-mentioned purpose, the present application provides the following solutions.

[0006] In a first aspect, the present application provides a method for predicting the wave height in the sea area after intervention of a wave energy device, which comprises the following steps.

[0007] Based on the main size of the wave energy device, a numerical tank is established and boundary conditions and numerical tank parameters are set to obtain an ideal wave numerical tank.

[0008] The ideal wave numerical tank is compared and analyzed with a theoretical wave sine curve to obtain an ideal wave tank numerical model after verification.

[0009] After modeling the wave energy device, the ideal wave pool numerical model is imported, the dynamic domain, static domain and boundary layer grid are set, the conversion efficiency of the wave energy device is calculated, and the conversion efficiency of the device obtained by the test is compared until the relative error range is met, and the final verified wave pool numerical model containing the device is obtained.

[0010] Based on the final verified wave pool numerical model containing the device, a period is divided into several equidistant time points, the wave surface change position at each time point is counted and averaged, and the average change data of the wave surface in front of the device and the wave surface behind the device in a period is obtained.

[0011] Based on the average change data of the wave surface in front of the device and the wave surface behind the device in a period, the least square fitting method is used to determine the wave surface height distribution function variable factor, that is, the wave height distribution function in front of the device based on the natural exponential type and the tail wave height distribution function behind the device based on the oscillation attenuation are obtained.

[0012] Based on the wave height distribution function in front of the device based on the natural exponential type and the tail wave height distribution function behind the device based on the oscillation attenuation, the wave surface height in front of and behind the wave energy device is calculated.

[0013] In a second aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored on the memory and capable of running on the processor, and the processor executes the computer program to realize the wave energy device intervention sea area wave surface height prediction method of the first aspect.

[0014] In a third aspect, the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the wave energy device intervention sea area wave surface height prediction method of the first aspect.

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

[0016] The present application provides a wave energy device intervention sea area wave surface height prediction method and related device, which comprises: based on the main size of the wave energy device, a numerical wave pool is established, boundary conditions and numerical wave pool parameters are set, an ideal wave numerical wave pool is obtained, and the initial conditions are matched with the water flow in and out, boundary constraints and other conditions in the actual wave environment.

[0017] The ideal wave numerical tank is compared with the theoretical wave sine curve, and the verified ideal wave tank numerical model is obtained. Through this comparison, it can be effectively tested whether the previously set numerical tank parameters and boundary condition settings are reasonable. If the two are highly consistent, it means that the ideal wave numerical tank is reliable, and then the verified ideal wave tank numerical model is obtained, ensuring the accuracy of the subsequent simulation of wave energy device conversion efficiency and wave surface changes.

[0018] After modeling the wave energy device and importing it into the ideal wave tank numerical model, the dynamic domain, static domain, and boundary layer grid are set, and the wave energy device conversion efficiency is calculated. The device conversion efficiency obtained by the test is compared until the relative error range is met, and the final verified wave tank numerical model containing the device is obtained. This can distinguish different properties in the simulation, accurately simulate the motion state of the wave energy device under the action of waves in the dynamic domain, and simulate the normal wave propagation without interference from the device motion in the static domain, providing a reasonable and effective model basis for further analysis of the interaction between the wave energy device and the wave. By continuously optimizing the grid settings and boundary conditions of these key areas until the preset conversion efficiency or the relative error range is met, the model can more accurately simulate the interaction between the wave energy device and the wave, providing a reliable model basis for subsequent prediction of wave height.

[0019] Based on the final verified wave tank numerical model containing the device, a period is divided into several equidistant time points, the wave surface change position at each time point is counted and averaged, and the average change data of the wave surface in front of and behind the device in a period is obtained. This approach can effectively eliminate the randomness and contingency of wave surface changes at a single time point, and obtain a more representative result that reflects the overall trend of wave surface changes, providing a stable and reliable data basis for subsequent fitting of wave surface change functions.

[0020] Based on the average change data of the wave surface in front of and behind the device in a period, the least square fitting method is used to determine the wave height distribution function variable factor, i.e., the wave height distribution function in front of the device based on the natural exponential type and the tail wave height distribution function behind the device based on the oscillation decay are obtained. The function can succinctly and accurately describe the change law of the wave surface with time and space in the form of mathematical expression, and convert the complex process of wave surface changes into an analyzable and calculable mathematical model, providing a mathematical tool for calculating the wave height in front of and behind the wave energy device, and greatly improving the accuracy and efficiency of predicting the wave height.

[0021] Based on the natural exponential-based wave height distribution function before the device and the oscillation attenuation-based tail wave height distribution function after the device, the wave surface height before and after the wave energy device is calculated. These specific wave surface height data intuitively and accurately reflect the actual state of the wave surface after the wave energy device is intervened, and the precise prediction of the wave surface height after the wave energy device is intervened is realized, thereby providing strong data support for accurate assessment of wave energy resources and optimization of wave energy device array.

[0022] In summary, the method provides important support and guarantee for wave energy resource assessment and array optimization research. By accurately predicting the wave surface height after the wave energy device is intervened, the prediction accuracy of wave energy resources and power generation can be improved, providing strong support for wave energy device array optimization and development and utilization, and helping wave energy to become a more reliable clean energy. At the same time, it can help designers fully consider the influence of waves when designing wave energy devices, optimize structural design, and improve the reliability and safety of wave energy devices. BRIEF DESCRIPTION OF DRAWINGS

[0023] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0024] Figure 1 The flowchart of a wave surface height prediction method after a wave energy device is intervened according to an embodiment of the present application is shown.

[0025] Figure 2 The overall framework diagram of a wave surface height prediction method after a wave energy device is intervened according to an embodiment of the present application is shown.

[0026] Figure 3 The flowchart of a least squares fitting method according to an embodiment of the present application is shown.

[0027] Figure 4 The structure diagram of a wave energy device according to an embodiment of the present application is shown.

[0028] Figure 5 The ideal wave numerical tank diagram according to an embodiment of the present application is shown.

[0029] Figure 6 The diagram for comparison between the extracted wave surface position and the theoretical wave sine curve according to an embodiment of the present application is shown.

[0030] Figure 7 The overall diagram of numerical tank grid division according to an embodiment of the present application is shown.

[0031] Figure 8 A numerical pool initial position diagram provided by an embodiment of the present application.

[0032] Figure 9 A boundary layer setting diagram provided by an embodiment of the present application.

[0033] Figure 10 A power conversion efficiency verification result diagram provided by an embodiment of the present application.

[0034] Figure 11 A wave surface diagram of an initial state provided by an embodiment of the present application.

[0035] Figure 12 A T / 2 moment simulation curve and a theoretical wave height comparison diagram provided by an embodiment of the present application.

[0036] Figure 13 A T / 2 moment wave surface relative change diagram in front of the device provided by an embodiment of the present application.

[0037] Figure 14 A T / 2 moment wave surface relative change diagram behind the device provided by an embodiment of the present application.

[0038] Figure 15 A wave surface change diagram in front of the device after period averaging provided by an embodiment of the present application.

[0039] Figure 16 A wave surface change diagram behind the device after period averaging provided by an embodiment of the present application.

[0040] Figure 17 A wave surface distribution function diagram in front of the device provided by an embodiment of the present application.

[0041] Figure 18 A wave surface distribution function diagram behind the device provided by an embodiment of the present application.

[0042] Figure 19 A structure diagram of a computer device provided by an embodiment of the present application. DETAILED DESCRIPTION

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

[0044] The application can quickly and accurately predict the wave height change and influence range of the wave field after the single wave device is intervened, thereby achieving the rapid prediction of the wave field change after the wave energy device is arranged in a large scale, improving the accuracy of wave energy resource assessment and total power generation calculation of the wave energy array field, supporting the optimization of wave energy device arrangement array, and reasonably arranging to improve the total power generation.

[0045] In order to make the above-mentioned purposes, features and advantages of the application more obvious and easy to understand, the application will be further described in detail below in combination with the drawings and specific embodiments.

[0046] In one exemplary embodiment, as shown in Figure 1 and Figure 2 , a wave energy device intervention sea area wave height prediction method is provided, which comprises the following steps.

[0047] S1: Based on the main size of the wave energy device, a numerical tank is established and boundary conditions and numerical tank parameters are set to obtain an ideal wave numerical tank; the main size includes height, length and width; the boundary conditions are that the left side of the tank is defined as an inlet Inlet, which is set as a velocity inlet (Velocity Inlet), the upper and lower sides are defined as Up and Down respectively, which are set as velocity inlets (Velocity Inlet) and no-slip walls (No slip Wall) respectively, and the right side outlet Outlet is defined as a pressure outlet (Pressure Outlet); the numerical tank parameters include turbulence model, multiphase flow, wave, gravity effect and time solving method; the ideal wave numerical tank is a tank established based on a Cartesian coordinate system.

[0048] S2: Comparing the ideal wave numerical tank with the theoretical wave sine curve, an ideal wave tank numerical model after verification is obtained.

[0049] S3: After modeling the wave energy device, it is introduced into the ideal wave tank numerical model, the dynamic domain, static domain and boundary layer grid are set, the wave energy device conversion efficiency is calculated, compared with the device conversion efficiency obtained by experiment, until the relative error range is met, and the final verified wave tank numerical model containing the device is obtained.

[0050] S4: Based on the final verified wave tank numerical model containing the device, a period is divided into several equidistant time points, the wave surface change position at each time point is counted and averaged to obtain the average change data of the wave surface before and after the device in a period.

[0051] S5: Based on the average change data of the wave surface in front of the device and the wave surface behind the device in one period, the least square fitting method is used to determine the variable factor of the wave height distribution function, that is, the natural exponential type wave height distribution function in front of the device and the oscillation attenuation type wave height distribution function behind the device are obtained.

[0052] S6: Based on the natural exponential type wave height distribution function in front of the device and the oscillation attenuation type wave height distribution function behind the device, the wave height in front of and behind the wave energy device is calculated.

[0053] By implementing the above steps S1 to S6, based on the numerical simulation method, the ideal wave numerical pool is constructed, the boundary conditions and numerical pool parameters are reasonably set, the wave energy device is modeled in the ideal wave numerical pool, the wave energy device is verified by the regular theoretical wave surface, and the wave height change around the wave energy device is calculated by comparing the change of the regular theoretical wave surface.

[0054] The overall scheme of the present application is shown as follows.

[0055] I. First, an ideal wave pool is established.

[0056] First, the size of the wave energy device to be modeled is obtained, mainly the height, length and width. The ideal wave numerical pool matching the device is established in the numerical model, the pool is established based on the Cartesian coordinate system, the positive direction of the vertical direction Z points downward, X the positive direction points to the right, Y the direction should follow the right-hand rule.

[0057] Set reasonable boundary conditions, the left side of the pool is defined as Inlet, set as velocity inlet, the upper and lower sides are defined as Up and Down, respectively, set as velocity inlet and no-slip wall, the right side outlet Outlet is defined as pressure outlet.

[0058] Hexahedral structured grid is adopted to establish the ideal wave numerical pool without device, and the grid is encrypted at the water-air interface. The turbulent flow model of the wave numerical pool adopts Reynolds time-averaged turbulent flow model, and the VOF multiphase flow is set. Other specific settings are shown in Table 1.

[0059] Table 1 Numerical pool parameters and setting types

[0060] The Reynolds time-averaged turbulent flow model can select k-e or SST k-w model, the wave adopts Fifth Order VOF wave, and the simulated wave height and wave period are set.

[0061] II. Verify the numerical wave tank simulation effect.

[0062] The simulated wave tank is compared with the theoretical wave sine curve to verify the model, and then the ideal wave tank numerical model after verification is obtained.

[0063] III. Model the wave energy device.

[0064] The wave energy device is modeled in three-dimensional modeling software, and the model is appropriately simplified to retain the main part of the energy harvesting structure. The moment of inertia, mass, center of gravity and other parameters are obtained. The device is imported into the numerical model, and the water-air interface is encrypted. At the same time, the boundary layer encryption method is used to encrypt the edge of the device.

[0065] IV. Verify the numerical model.

[0066] The wave energy device is set as a 6-DOF (6-DOF) motion model, and the moment of inertia, center of gravity and other parameters of the wave energy device are set. The calculation domain is divided into static and dynamic domains. The numerical model is verified by calculating the conversion efficiency of the wave energy device.

[0067] V. Wave surface change analysis.

[0068] After the numerical model is verified, the liquid surface change caused by the wave energy device is analyzed. The theoretical height generated by the theoretical wave surface sine function is compared with the wave surface after the device is involved, and the period is divided into M equidistant time points. The change of the wave surface at each time point is calculated, and the average change is calculated. According to this, the period-averaged wave surface height change can be obtained.

[0069] VI. Fit the wave surface height change function.

[0070] The wave surface change is divided into two parts before and after the device, and the least square method and other fitting methods are used for fitting to obtain the wave height distribution function before and after the device. According to this, the wave surface change before and after the device can be quickly calculated.

[0071] As an optional implementation, in step S1, the ideal wave numerical tank size is determined. First, the height and length of the wave energy device to be modeled are obtained. The ideal wave numerical tank is established, and the overall height of the tank is at least 5 times the height of the device, and the width is at least 5 times the width of the device.

[0072] .

[0073] .

[0074] wherein, Height of the ideal wave numerical tank; Width of the ideal wave numerical tank; Height of the wave energy device; Width of the wave energy device.

[0075] In terms of length, the overall length is at least 31 times the length of the device: at least 10 times the length in front of the device and at least 20 times the length behind the device.

[0076] .

[0077] .

[0078] .

[0079] wherein, Length of the ideal wave numerical tank; Length in front of the wave energy device; Length behind the wave energy device; Length of the wave energy device.

[0080] According to the above dimensions, an ideal wave numerical tank is established. First, the coordinate system is based on the Cartesian coordinate system, and the vertical direction Z is in the positive direction, pointing downward, X the positive direction points to the right, Y the direction should follow the right-hand rule. X The point where =0 is set as the origin of the wave energy device to be placed.

[0081] Boundary conditions are set, the left side of the tank is defined as Inlet, set as a velocity inlet, the top and bottom are defined as Up and Down, respectively, set as a velocity inlet and a no-slip wall, respectively, and the right side outlet Outlet is defined as a pressure outlet, as shown in Table 2.

[0082] Table 2 Boundary and setting type

[0083] The numerical model (ideal wave numerical tank) is meshed using hexahedral structured mesh, and the overall mesh size should be selected as 0.1 times the smaller value of the device height or length. First, determine which one is smaller, the height or the length, and the overall mesh size of the tank should be based on: 0.1 times the value.

[0084] .

[0085] .

[0086] wherein, This is the smaller of the height and length of the wave energy device; The overall grid size of the ideal wave numerical water tank.

[0087] Furthermore, the water-air interface is densified, and the minimum grid size is shown below.

[0088] .

[0089] in, The minimum grid size for an ideal wave numerical water tank.

[0090] The turbulence model, multiphase flow, wave and time solution methods of the numerical model are set as shown in Table 1. The Reynolds time-averaged turbulence model is selected for the turbulence model, the VOF multiphase flow model is used for the multiphase flow, and the fifth-order wave model is used for the wave. Considering the effect of gravity, the implicit unsteady calculation format is selected for the time calculation.

[0091] The wave model is set as Fifth Order VOF Wave. For the specific theory of this wave model, please refer to Fenton JD's paper "A Fifth-Order Stokes Theory For Steady Waves", Journal of Waterway Port Coastal and Ocean Engineering, 1985, 111(2): 216-234; this application will not elaborate further.

[0092] The Reynolds-averaged turbulence model can be selected from the Ke model, the Kw model, or the SST Kw model. Regarding the waterline, a minimum distance of 3.5 should be set from the bottom of the pool. The height is the reference position of the water-air interface (waterline).

[0093] Set the wave height of the regular wave to H=0.2. The wave period is T, and the minimum computation time step is set to T / 1000. After calculating the model for 10T, the model is considered stable. After extracting 10T... x The wavefront position at =0 was continuously sampled for 5T time intervals, and the wavefront undulations were dimensionless. Let the waterline position be... , amplitude A=0.5 H; The position of the liquid surface at a certain moment is The relative change in liquid level can be calculated by subtracting the baseline position from the water level rise and then dividing by the amplitude; this gives the result. The range of variation is [-1, 1].

[0094] As an optional implementation, in step S2, the ideal wave numerical tank is compared with the theoretical wave sine curve, and the ideal wave tank numerical model after verification is obtained, specifically including the following steps.

[0095] S21: Set the wave height, wave period and minimum calculation time step of the regular wave.

[0096] S22: Extract the wave surface position at =0 after the preset wave period x , continuously extract several wave periods, and carry out dimensionless treatment on the wave surface fluctuation.

[0097] S23: Compare the extracted wave surface position with the theoretical wave sine curve to obtain the difference between the theoretical value and the actual value.

[0098] S24: Determine whether the difference exceeds the preset threshold.

[0099] S25: If yes, the ideal wave numerical tank fails the verification, and the numerical tank parameters are reset.

[0100] S26: If no, the ideal wave numerical tank passes the verification, and the ideal wave tank numerical model after verification is obtained.

[0101] In step S2, the extracted wave surface position is compared with the theoretical sine value to obtain the difference between the theoretical value and the actual value; the theoretical sine value can be generated according to the following formula.

[0102] .

[0103] Wherein, A is the amplitude, m; is the angular frequency, rad / s; ; is the initial phase offset, rad; t is the time, s; T is the wave period of the current working condition.

[0104] According to the comparison result, the model boundary condition and the like can be adjusted to obtain the numerical tank (ideal wave tank numerical model after verification) best matching the theoretical value.

[0105] As an optional implementation, in step S3, wave energy device three-dimensional modeling simplification is performed, Solidworks or other three-dimensional modeling software is used to appropriately simplify the wave energy device, the main structure is retained, and the moment of inertia, mass, center of gravity and other parameters of the device are obtained.

[0106] The simplified model is imported into the ideal wave pool numerical model, the surrounding area of the wave energy device is set as a dynamic domain, the remaining area is set as a static domain, and the wave energy device motion model is set as a 6 degree of freedom (6DOF) motion mode. The moment of inertia, mass, center of gravity and other parameters of the wave energy device are input into the model. The remaining waterline position, boundary conditions and foundation grid division are consistent with the verified ideal wave pool.

[0107] As an optional implementation, in step S3, the wave energy device is modeled and imported into the ideal wave pool numerical model, the dynamic domain, static domain and boundary layer grid are set, and the wave energy device conversion efficiency is calculated. The calculated conversion efficiency is compared with the conversion efficiency obtained by the test, and the relative error range is satisfied to obtain the final verified wave pool numerical model containing the device, which includes the following steps.

[0108] S31: The wave energy device is modeled and imported into the ideal wave pool numerical model, wherein the surrounding area of the wave energy device is set as a dynamic domain, and the other areas are set as a static domain.

[0109] S32: The dynamic domain is encrypted, and the boundary layer grid is set around the wave energy device to calculate the first layer boundary layer grid height.

[0110] S33: Based on the first layer boundary layer grid height, the boundary layer thickness and growth rate are set to obtain the final grid division numerical pool.

[0111] S34: Based on the final grid division numerical pool, the wave energy device conversion efficiency is calculated.

[0112] S35: When the calculated wave energy device conversion efficiency and the conversion efficiency obtained by the test cannot satisfy the error requirement, the boundary layer grid or boundary condition parameter of the dynamic domain and the surrounding area of the device is adjusted until the relative error range is satisfied, and the final verified wave pool numerical model containing the device is obtained. Specifically, the surrounding water area of the wave energy device is encrypted, prism grid is used, and boundary layer grid is set around the device to capture the complex wave changes around the device. The first layer boundary layer grid height calculation method is as follows.

[0113] 1) First, calculate the Reynolds number.

[0114] .

[0115] wherein, Reynolds number; reference velocity; reference length; hydrodynamic viscosity.

[0116] The hydrodynamic viscosity is: For wave energy device, in the general sea environment, the reference velocity The value can be taken as follows.

[0117] .

[0118] Where, A is the amplitude, 0.5 H , is the frequency, ; the reference length can be taken as the device length value.

[0119] 2) Then calculate the friction coefficient: the friction coefficient can be estimated by an empirical formula.

[0120] .

[0121] Where, is the friction coefficient.

[0122] 3) Then calculate the friction velocity: the friction velocity estimation method is shown as follows.

[0123] .

[0124] Where, is the friction velocity.

[0125] 4) Finally, solve the first layer boundary layer grid height.

[0126] .

[0127] Where, is the first layer boundary layer grid height; is the wall dimensionless distance. According to the calculation requirement, the wall dimensionless distance value can generally be taken as 1~20; finally, the first layer grid height of the boundary layer .

[0128] As an optional implementation, the boundary layer thickness and growth rate are set, which can generally be set to 5~20 layers, and the growth rate can be set to 1.1~1.3.

[0129] As an optional implementation, through the model output result, first calculate the device time-averaged output power, and finally verify the model according to the wave energy device conversion efficiency calculation result.

[0130] The device time-averaged output power calculation formula is as follows.

[0131] .

[0132] Where, is the time-averaged output power of the device, kW; is the area enclosed by the device output power duration curve and the time axis, kW·s; t is the total calculation time.

[0133] The conversion efficiency of the wave energy device is the ratio of the power obtained by the device to the wave input power. The wave input power can be calculated by the regular wave incident power. The regular wave incident power can be obtained by multiplying the regular wave power density and the device upwind width.

[0134] The calculation expression of the regular wave power density is as follows.

[0135] .

[0136] wherein, J is the regular wave power density, kW / m; is the seawater density, ; g is the gravitational acceleration, ; H is the regular wave height, m, T is the wave period, s.

[0137] The calculation formula of the incident power of the regular wave is as follows.

[0138] .

[0139] wherein, is the incident power of the regular wave, kW; B is the upwind width, m.

[0140] Finally, the calculation formula of the conversion efficiency of the wave energy device is as follows.

[0141] .

[0142] wherein, is the conversion efficiency of the wave energy device.

[0143] As an optional implementation, if the conversion efficiency verification effect is not good (for example, the threshold value can be set to within 15% of the relative error), the device surrounding grid division form needs to be adjusted until the verification is good.

[0144] As an optional implementation, in step S4, after the accurate numerical model is verified, the wave surface change around the wave energy device can be analyzed. In order to better analyze the periodic change of the wave surface, one period is divided into M time points (M≥6), and the wave surface position at each time point is recorded as , x is the X-direction distance ，t respectively, M time points in one period.

[0145] The wave surface change is divided into before and after the device, and the wave surface change is converted into the change of the relative theoretical value. The specific method is as follows.

[0146] First, the wave surface of the theoretical value at each time is obtained, and the theoretical liquid surface at each position at each time is , and the wave surface relative position is calculated according to the following formula.

[0147] .

[0148] When is negative, it means that the wave surface height after the device intervention is lower than the theoretical value, and when is positive, it means that the wave surface height after the device intervention is higher than the theoretical value.

[0149] The water surface change is converted into the relative change, and the wave surface change at each time is extracted, and finally the changes of M times (M≥6) are averaged, and then the average change of the wave surface in a period can be obtained.

[0150] According to the regular wave power density formula, the wave energy density is proportional to the square of the wave height, that is, the larger the wave height, the more energy is contained and transmitted in the water body at this time.

[0151] Therefore, the wave surface change needs to be weighted and averaged. When the wave height is high, the wave energy contained in the water body is larger, and the interaction between the device and the wave is most intense, and the wave surface change when the energy is large needs to be considered.

[0152] The following weighted average method based on spatial change is used to calculate the weight of each position x at each time. The theoretical wave height at this position is as follows.

[0153] .

[0154] The above formula shows that is X = x the theoretical wave height at position at time , and the square is used to reflect the proportional relationship between energy and wave height. The absolute value is used to make the weight always positive, and the wave crest and the wave trough have equal importance.

[0155] This method can reflect the different contributions of different spatial positions to the average wave surface change due to different wave phases during wave propagation, and especially emphasizes the disturbance effect of high-energy areas (wave crest / wave trough).

[0156] Finally, the weight corresponding to M times is equal to 1, that is, .

[0157] Then the weighted average wave face variation is calculated: the wave face variation after period average, for each position x , using its exclusive weight to calculate the weighted average.

[0158] .

[0159] Wherein, is the average variation data of the wave face in front of the device and the wave face behind the device within a period; is the wave face relative position at the position at the moment x ; is the weight at the position at the moment x ; M is the total moment.

[0160] Optionally, the wave face variation is averaged by using an average function, i.e. a general average method is selected, and a weighted average method is not used.

[0161] .

[0162] According to the above formula, the wave face relative variation in front of the device is obtained, ; and the wave face relative variation behind the device is obtained, .

[0163] As an optional implementation, in step S5, in order to further obtain the wave face variation function in front of the device and behind the device, a least square fitting method (threshold value is that the determination coefficient is at least 75%) is used to determine the wave face variation function variable factor in front of the device and behind the device (the fitting process is shown in Figure 3 ). When the wave propagates to the device, there is reflection and energy accumulation, and both effects cause the wave face to have a lifting phenomenon. And this phenomenon is enhanced as it approaches the device, showing an exponential growth law.

[0164] Therefore, the natural exponential type wave height distribution function in front of the device can be obtained.

[0165] .

[0166] Wherein, is the natural exponential type wave height distribution function in front of the device; is an amplitude factor; is an accumulation factor; C is a wave-approaching offset factor; x is X the direction distance. The function can reflect the wave height lifting in front of the device, and the wave face relative variation tends to zero in the far distance away from the device, i.e. close to the theoretical wave face.

[0167] The least square fitting method is used to determine 、 and C , and the incident wave height distribution function can be obtained.

[0168] The wave energy device extracts energy from the sea, and the total energy decays. There is a phenomenon of wave height reduction on the wave surface behind the device. According to the principle of energy conservation, and because the wave energy density is proportional to the square of the wave height ( E ∝ ), the reduction of energy will inevitably lead to the reduction of the transmitted wave height. The change of the wave surface behind the device can be expressed as a function of oscillation decay with the propagation distance.

[0169] Therefore, the tail wave height distribution function behind the device based on oscillation decay can be obtained.

[0170] .

[0171] wherein, is the tail wave height distribution function behind the device based on oscillation decay; is the tail wave amplitude factor; is the frequency, ; is the phase shift factor, is the tail wave shift factor; is the tail wave decay factor. This function can indicate that the wave surface decays behind the device, and slowly approaches the theoretical wave surface away from the device.

[0172] The least square fitting method is used to determine 、 、 and , and the tail wave height distribution function can be obtained.

[0173] .

[0174] wherein, is the relative change function of the wave surface before and after the device.

[0175] Finally, in step S6, the prediction expression of the wave surface height before and after the device can be uniformly expressed as follows.

[0176] .

[0177] wherein, is the wave surface height before and after the wave energy device; is the theoretical wave surface height.

[0178] Through the method, the wave surface height distribution at any position around the wave energy device can be quickly calculated without re-performing numerical simulation, thereby providing quantitative basis for array layout optimization and device interference analysis.

[0179] The method realizes the rapid prediction of the wave surface height before and after the wave energy device by converting the complex instantaneous wave surface change into periodic average change and further fitting into a continuous function. Compared with the traditional numerical simulation method, the method has the advantages of high calculation efficiency, strong applicability, and easy engineering application, and is especially suitable for wave field interference analysis, resource assessment and array optimization after the array layout of the wave energy device.

[0180] The application will be described below with specific embodiments.

[0181] The wave energy device to be studied is shown in Figure 4 , with a height of 5 m, a width of 7 m, and a length of 1 m; an ideal numerical tank is established according to the above dimensions (as shown in Figure 5 ): the boundary conditions and other parameters are set according to the description above, the regular wave height is set to 1 m, and the wave period is set to 4 s. The calculation step is set to 0.004 s.

[0182] The water surface change of the position X =0 for 5 consecutive periods is extracted and compared with the ideal wave sine curve: the ideal sine curve is , in order to facilitate comparison, the period is dimensionless; the comparison result is shown in Figure 6 .

[0183] The correlation coefficient reaches 98.5% for 5 consecutive periods, proving the feasibility of the numerical model. After verifying the ideal wave tank, the wave energy device is introduced into the numerical model, and the overall grid division diagram of the numerical tank is obtained, as shown in Figure 7 ; the initial position of the numerical tank is shown in Figure 8 ; the boundary layer setting is shown in Figure 9 ; and the conversion efficiency verification result is shown in Figure 10 .

[0184] The relative error of the highest point conversion efficiency is 5.71%, and the overall trend is close, so it is considered that the numerical model is verified, and the simulation result is reliable.

[0185] Now start the wave surface analysis, which is divided into 6 time points. Figure 11 is the initial state of the wave surface: (asterisk represents the position of the wave energy device); Figure 12 is the comparison between the simulation curve and the theoretical wave height at T / 2; Figure 13 is the relative change of the wave surface in front of the device at T / 2; Figure 14 is the relative change of the wave surface after the device at T / 2.

[0186] The wave front changes at 6 time points were averaged, and the wave front changes before the device after the period average were obtained as shown in Figure 15 ; the wave front changes after the device after the period average were as shown in Figure 16 .

[0187] The wave front height distribution function variable factors before and after the device were determined by using the least square method, and the wave front height distribution functions before and after the device were obtained as shown below.

[0188] .

[0189] .

[0190] The wave front height distribution function before the device is shown in Figure 17 , and the tail wave height distribution function after the device is shown in Figure 18 . The results show that the wave function is very consistent with the wave value of the numerical simulation result, the root mean square error RMSE is 0.0063; the determination coefficient is 84.4%, the root mean square error RMSE is 0.00283; the determination coefficient is 76.6%.

[0191] The relative change functions of the wave front before and after the device can be expressed as follows.

[0192] .

[0193] Finally, the prediction expression of the wave front height before and after the device can be expressed as follows.

[0194] .

[0195] In an exemplary embodiment, a computer device is provided, which can be a server or a terminal, and its internal structure diagram can be as shown in Figure 19As shown in the figure. The computer device includes a processor, a memory, an Input / Output (I / O) interface and a communication interface. Among them, the processor, the memory and the input / output interface are connected through the system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capability. 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 the main size of the wave energy device. The input / output interface of the computer device is used to exchange information between the processor and the external device. The communication interface of the computer device is used to communicate with the terminal outside through the network connection. The computer program is executed by the processor to realize a wave energy device intervention post-sea wave height prediction method.

[0196] Those skilled in the art can understand that, Figure 19 The structure shown in the figure is only a block diagram of part of the structure related to the scheme of the present application, and does not constitute a limitation on the computer device to which the scheme of the present application is applied. The specific computer device can include more or fewer components than those shown in the figure, or combine certain components, or have a different component arrangement.

[0197] In one exemplary embodiment, a computer device is also provided, including a memory and a processor, the memory storing a computer program, and the processor executing the computer program to realize the above-mentioned method embodiments.

[0198] In one exemplary embodiment, a computer readable storage medium is provided, storing a computer program, which is executed by a processor to realize the above-mentioned method embodiments.

[0199] The technical features of the above embodiments can be combined in any way. In order to make the description concise, not all possible combinations of the technical features in the above embodiments are described, but as long as the combination of the technical features does not exist contradictory, it should be considered as the scope of the present application.

[0200] The principles and implementation modes of the present application are described by applying specific examples herein. The above embodiment is only used to help understand the method and its core idea of the present application; at the same time, for those skilled in the art, according to the idea of the present application, the specific implementation mode and application range will be changed. In view of the above, the content of the specification should not be understood as a limitation of the present application.

Claims

1. A method for predicting wave height in a post-intervention sea area of a wave energy device, characterized in that, The method comprises the following steps: Based on the size of the wave energy device, a numerical wave tank is established, and boundary conditions and numerical wave tank parameters are set to obtain an ideal wave numerical wave tank; The ideal wave numerical wave tank is compared with the theoretical wave sine curve to obtain the ideal wave numerical wave tank after verification; The wave energy device is modeled and imported into the ideal wave numerical wave tank, the dynamic domain, the static domain and the boundary layer grid are set, and the conversion efficiency of the wave energy device is calculated, which is compared with the conversion efficiency obtained by the test, until the relative error range is met, and the wave numerical wave tank containing the device after final verification is obtained; Based on the wave numerical wave tank containing the device after final verification, a period is divided into several equidistant time points, the wave surface change position at each time point is counted and averaged to obtain the average change data of the wave surface in front of and behind the device in a period; Based on the average change data of the wave surface in front of and behind the device in a period, the least square fitting method is used to determine the wave surface height distribution function variable factor, that is, the wave height distribution function in front of the device based on the natural exponential type and the tail wave height distribution function behind the device based on the oscillation attenuation are obtained; Based on the wave height distribution function in front of and behind the device, the wave height in front of and behind the wave energy device is calculated.

2. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The ideal wave numerical wave tank is compared with the theoretical wave sine curve to obtain the ideal wave numerical wave tank after verification, which specifically comprises: The wave height, wave period and minimum calculation time step of the regular wave are set; extracting a preset wave period x the wave surface position at =0, continuously extracting several wave periods, and normalizing the wave surface fluctuation; The extracted wave surface position is compared with the theoretical wave sine curve to obtain the difference between the theoretical value and the actual value; It is judged whether the difference exceeds the preset threshold; If yes, the ideal wave numerical wave tank does not pass the verification, and the numerical wave tank parameters are reset; If no, the ideal wave numerical wave tank passes the verification, and the ideal wave numerical wave tank after verification is obtained.

3. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The wave energy device is modeled and imported into the ideal wave numerical wave tank, the dynamic domain, the static domain and the boundary layer grid are set, and the conversion efficiency of the wave energy device is calculated, which is compared with the conversion efficiency obtained by the test, until the relative error range is met, and the wave numerical wave tank containing the device after final verification is obtained, which specifically comprises: The wave energy device is modeled and imported into the ideal wave numerical wave tank, wherein the periphery of the wave energy device is set as the dynamic domain, and other areas are set as the static domain; The dynamic domain is encrypted, and the boundary layer grid is set around the wave energy device to calculate the first layer boundary layer grid height; Based on the first layer boundary layer grid height, the boundary layer thickness and the growth rate are set to obtain the final grid division numerical wave tank; Based on the final grid division numerical wave tank, the conversion efficiency of the wave energy device is calculated; When the calculated conversion efficiency of the wave energy device and the conversion efficiency obtained by the test cannot meet the error requirement, the boundary layer grid or the boundary condition parameter of the dynamic domain and the area around the device is adjusted until the relative error range is met, and the wave numerical wave tank containing the device after final verification is obtained.

4. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The calculation formula of the wave energy device conversion efficiency is: ； ； ； ； wherein, is the wave energy device conversion efficiency; is the time-averaged output power of the device; is the incident power of regular waves; B is the oncoming wave width; J is the regular wave power density; is the seawater density; g is the gravitational acceleration; H is the regular wave height, T is the wave period; is the area enclosed by the device output power duration curve and the time axis; t is the total calculation time.

5. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The formula for calculating the average change of the wave surface in front of and behind the device in one period is: or ; wherein, is the average change of the wavefront before and after the device in a cycle; is the time instant; x is the relative position of the wavefront at the position at the time instant; is the weight at the position at the time instant; x M is the total time instant.​​​ 6. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The expression of the natural exponential type wave height distribution function in front of the device is: ； wherein, is a natural exponential based wavefront height distribution function prior to the device; is an amplitude factor; is a packing factor; C is a wavefront offset factor; x is X directional distance.

7. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The expression of the oscillation decay type wave height distribution function behind the device is: ； wherein, is an apparatus-based coda height distribution function based on oscillatory decay; is a coda amplitude factor; is a frequency; is a phase shift factor, is a coda shift factor; is a coda decay factor; x is X direction distance.

8. The wave energy device post-intervention sea area wave face height prediction method of claim 1, wherein, The formula for calculating the wave surface height in front of and behind the wave energy device is: ； ； wherein, is the wave height in front of the wave energy device; is the theoretical wave height; is the relative change function including the wave height in front of the device and the wave height behind the device; is the wave height distribution function in front of the device based on a natural exponential type; is the wave height distribution function behind the device based on an oscillation decay.

9. A computer device comprising: The memory, the processor and the computer program stored in the memory and capable of running on the processor, characterized in that the processor executes the computer program to realize the wave surface height prediction method after the wave energy device is intervened in the sea area according to any one of claims 1-8.

10. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program is executed by the processor to realize the wave surface height prediction method after the wave energy device is intervened in the sea area according to any one of claims 1-8.

Images