Wave energy device array power estimation output method based on wave directional spectrum analysis

Abstract

The application discloses a wave energy device array power estimation output method based on wave directional spectrum analysis, first, the measured wave data of a target sea area is obtained, and a wave climate database containing effective wave height, energy period or peak period, frequency spectrum and direction distribution is established; then, a multidirectional wave directional spectrum is reconstructed according to a one-dimensional wave spectrum and a direction distribution function, and multiple combinations are constructed; further, a wave energy device array hydrodynamic state space model is established, the average capture power, power matrix and annual power generation of the array under different sea conditions are calculated in combination with PTO damping parameters and multidirectional irregular wave excitation; finally, taking the measured spectrum-measured direction distribution combination as a benchmark, the power estimation deviation under other combinations is calculated and corrected, and the array power estimation result facing the target sea area is output, so that the risk of overestimation or underestimation of power caused by unidirectional waves and idealized theoretical wave spectrum assumptions is reduced, and the power estimation reliability of the wave energy device array under real multidirectional sea conditions is improved.

Description

Technical Field

[0001] This invention relates to the field of marine renewable energy and wave energy power generation device array modeling and power estimation technology, specifically to a wave energy device array power estimation output method based on wave direction spectrum analysis. Background Technology

[0002] Wave energy, with its wide distribution, high energy density, and ability to be deployed in conjunction with offshore platforms or marine observation equipment, represents an important direction for the development of marine renewable energy. To enhance energy harvesting capacity per unit area of ​​sea, multiple wave energy generation devices are typically arranged in an array. While array arrangements allow for the sharing of mooring, power transmission, and maintenance resources, they also introduce radiation, diffraction, and shielding effects between devices, creating a strong coupling relationship between array power and wave incident direction, wave spectrum morphology, and array geometry.

[0003] In the early design and techno-economic evaluation of wave energy device arrays, theoretical wave spectrum and single incident direction are often used to quickly estimate the power of the device. This method is simple to implement, but it differs from real sea conditions. Real sea conditions are usually characterized by multi-directional irregular waves with both frequency and direction dispersion, and the direction distribution in different sea areas may exhibit asymmetric or non-Gaussian characteristics. If only the assumption of unidirectional waves is used, it is easy to amplify the favorable tectonic interference in the array, resulting in an overestimation of annual power generation and power benefits. If an over-diffuse theoretical direction distribution is used, the array power may also be underestimated.

[0004] For a single axisymmetric point-absorbing wave energy device, the directional sensitivity may not be obvious under ideal conditions; however, when multiple devices are arranged in linear, triangular, square, or other arrays, the overall geometry of the array exhibits significant differences in response to the incident direction. Therefore, it is difficult to reliably infer the annual power generation and the suitability of the array's layout under real sea conditions based solely on the performance of a single device or the array's performance under unidirectional wave conditions. Summary of the Invention

[0005] This invention provides a power estimation output method for wave energy device arrays based on wave direction spectrum analysis. The method is reasonably designed and can simultaneously consider wave power spectrum type, wave direction distribution, array hydrodynamic coupling effect, and measured sea state characteristics of the target sea area. It avoids power estimation deviations caused by using only theoretical wave spectrum and unidirectional wave assumptions, thus more accurately reflecting the power capture capability of wave energy device arrays under multidirectional irregular wave conditions in the target sea area. It provides a basis for wave energy device array engineering design, array layout, power prediction, and economic evaluation, thereby solving the problem that existing array power estimation methods, such as theoretical wave spectrum and unidirectional wave assumptions, are prone to estimation deviations and fail to reflect the real multidirectional wave environment.

[0006] The technical solution adopted by the present invention to solve the above-mentioned technical problems is as follows:

[0007] A method for estimating the output power of a wave energy device array based on wave direction spectrum analysis, the method comprising the following steps:

[0008] S1, determine the wave energy device array object to be estimated, determine the spatial arrangement and spacing parameters of each wave energy device, and obtain measured wave data of the target sea area; the wave energy device is a heave-type point absorption wave energy device, and the array is composed of multiple wave energy devices.

[0009] S2. Establish the hydrodynamic motion model and wave climate database corresponding to the wave energy device array. Incorporate the hydrostatic restoring force, radiation force, PTO force and wave excitation force on the wave energy device into the array dynamic equation. Referring to the hydrodynamic coupling between different wave energy devices in the array, perform state space approximation on the convolution integral term in the radiation force to construct a continuous state space model of the wave energy device array.

[0010] S3. Construct a multi-directional wave excitation model, combine the one-dimensional wave power spectrum with the wave direction distribution function to form a wave direction spectrum, and generate a multi-directional irregular wave field based on the wave direction spectrum.

[0011] S4. Obtain measured wave climate data for the target sea area. The measured wave climate data includes at least significant wave height, energy cycle or peak period, wave power spectrum and wave direction distribution information. Establish a sea state statistical matrix based on significant wave height and energy cycle.

[0012] S5. Determine the PTO damping parameters of the wave energy device according to the dominant frequency or peak period of each sea state, so that the wave energy device can obtain a higher average capture power under the corresponding sea state.

[0013] S6 compares the power estimation results under different combined operating conditions with the power estimation results under the benchmark operating condition, calculates the power gain or loss index, obtains the degree of influence of wave power spectrum type and direction distribution on the power estimation results of wave energy device array, and outputs the array power estimation results.

[0014] The equation of motion for any single floating body in the array can be expressed as:

[0015] ;

[0016] Among them, f *i b (t), f *i r (t), f *i pto (t), f *i e(t) represents the hydrostatic restoring force, radiation force, PTO force, and wave excitation force on the corresponding buoy, respectively;

[0017] The static restoring force can be expressed as:

[0018] ;

[0019] Where, k *i h x is the hydrostatic stiffness coefficient. *i (t) represents the undulating displacement of the pontoon;

[0020] The radiative force can be expressed as:

[0021] ;

[0022] Where, μ *i*j For the added mass at infinite frequency, This is the radiation impulse response function of the corresponding radiative force; The buoy's heave speed; It can be approximated in state-space form using MATLAB, thereby accelerating the computation;

[0023] The general state space can be obtained in the following form:

[0024] .

[0025] The wave energy device includes a cylindrical float and an underwater reaction column or reaction plate. The cylindrical float serves as a wave energy capture component, and the underwater reaction column or reaction plate serves as a stabilizing component, enabling the wave energy device to capture energy under heave freedom.

[0026] The spatial arrangement includes linear arrays, equilateral triangular arrays, and square arrays; the spacing between adjacent wave energy devices in the array is an integer multiple of the device's characteristic diameter;

[0027] The wave direction spectrum is represented as:

[0028] ;

[0029] Where S(ω) is the one-dimensional wave power spectrum, D(θ) is the wave direction distribution function, ω is the wave angular frequency, and θ is the wave propagation direction;

[0030] The wave direction distribution function satisfies the normalization condition, making the total energy of the wave direction spectrum consistent with the total energy of the one-dimensional wave power spectrum.

[0031] The combined operating conditions include theoretical wave power spectrum with unidirectional wave direction distribution, theoretical wave power spectrum with theoretical direction distribution, theoretical wave power spectrum with measured direction distribution, measured wave power spectrum with unidirectional wave direction distribution, measured wave power spectrum with theoretical direction distribution, and measured wave power spectrum with measured direction distribution; wherein, the measured wave power spectrum with measured direction distribution is used as the reference operating condition.

[0032] The power gain or loss metric is:

[0033] ;

[0034] Where G is the power gain or loss, P i To estimate the power under the combined operating conditions to be compared, P r This is the power estimate under the baseline operating condition;

[0035] The array power estimation results include the sea state average power matrix, annual power generation, and power estimation deviations under different wave power spectrum and direction distribution conditions.

[0036] This invention employs the aforementioned structure, establishing a hydrodynamic motion model corresponding to the wave energy device array. The static restoring force, radiation force, PTO force, and wave excitation force experienced by the wave energy devices are incorporated into the array's dynamic equations. Referring to the hydrodynamic coupling between different wave energy devices in the array, a state-space approximation is performed on the convolution integral term in the radiation force, constructing a continuous state-space model of the wave energy device array. By constructing a multi-directional wave excitation model, a one-dimensional wave power spectrum is combined with the wave direction distribution function to form a wave direction spectrum. By comparing the power estimation results under different combined operating conditions with the power estimation results under the baseline operating condition, power gain or loss indices are calculated to obtain the degree of influence of the wave power spectrum type and direction distribution on the power estimation results of the wave energy device array, and the array power estimation results are output. This method has the advantages of being accurate, reliable, simple, and practical. Attached Figure Description

[0037] Figure 1 This is a schematic diagram of the technical route of the present invention.

[0038] Figure 2 This is a reconstruction diagram of the wave direction spectrum of the present invention.

[0039] Figure 3 This is a diagram showing the estimated output effect of the wave energy device array of the present invention. Detailed Implementation

[0040] To clearly illustrate the technical features of this solution, the invention will be described in detail below through specific implementation methods and in conjunction with the accompanying drawings.

[0041] like Figures 1-3As shown, the wave energy device array power estimation output method based on wave direction spectrum analysis includes the following steps:

[0042] S1, determine the wave energy device array object to be estimated, determine the spatial arrangement and spacing parameters of each wave energy device, and obtain measured wave data of the target sea area; the wave energy device is a heave-type point absorption wave energy device, and the array is composed of multiple wave energy devices.

[0043] S2. Establish the hydrodynamic motion model and wave climate database corresponding to the wave energy device array. Incorporate the hydrostatic restoring force, radiation force, PTO force and wave excitation force on the wave energy device into the array dynamic equation. Referring to the hydrodynamic coupling between different wave energy devices in the array, perform state space approximation on the convolution integral term in the radiation force to construct a continuous state space model of the wave energy device array.

[0044] S3. Construct a multi-directional wave excitation model, combine the one-dimensional wave power spectrum with the wave direction distribution function to form a wave direction spectrum, and generate a multi-directional irregular wave field based on the wave direction spectrum.

[0045] S4. Obtain measured wave climate data for the target sea area. The measured wave climate data includes at least significant wave height, energy cycle or peak period, wave power spectrum and wave direction distribution information. Establish a sea state statistical matrix based on significant wave height and energy cycle.

[0046] S5. Determine the PTO damping parameters of the wave energy device according to the dominant frequency or peak period of each sea state, so that the wave energy device can obtain a higher average capture power under the corresponding sea state.

[0047] S6 compares the power estimation results under different combined operating conditions with the power estimation results under the benchmark operating condition, calculates the power gain or loss index, obtains the degree of influence of wave power spectrum type and direction distribution on the power estimation results of wave energy device array, and outputs the array power estimation results.

[0048] The equation of motion for any single floating body in the array can be expressed as:

[0049] ;

[0050] Among them, f *i b (t), f *i r (t), f *i pto (t), f *i e (t) represents the hydrostatic restoring force, radiation force, PTO force, and wave excitation force on the corresponding buoy, respectively;

[0051] The static restoring force can be expressed as:

[0052] ;

[0053] Where, k *i h x is the hydrostatic stiffness coefficient. *i (t) represents the undulating displacement of the pontoon;

[0054] The radiative force can be expressed as:

[0055] ;

[0056] Where, μ *i*j For the added mass at infinite frequency, This is the radiation impulse response function of the corresponding radiative force; The buoy's heave speed; It can be approximated in state-space form using MATLAB, thereby accelerating the computation;

[0057] The general state space can be obtained in the following form:

[0058] .

[0059] The wave energy device includes a cylindrical float and an underwater reaction column or reaction plate. The cylindrical float serves as a wave energy capture component, and the underwater reaction column or reaction plate serves as a stabilizing component, enabling the wave energy device to capture energy under heave freedom.

[0060] The spatial arrangement includes linear arrays, equilateral triangular arrays, and square arrays; the spacing between adjacent wave energy devices in the array is an integer multiple of the device's characteristic diameter;

[0061] The wave direction spectrum is represented as:

[0062] ;

[0063] Where S(ω) is the one-dimensional wave power spectrum, D(θ) is the wave direction distribution function, ω is the wave angular frequency, and θ is the wave propagation direction;

[0064] The wave direction distribution function satisfies the normalization condition, making the total energy of the wave direction spectrum consistent with the total energy of the one-dimensional wave power spectrum.

[0065] The combined operating conditions include theoretical wave power spectrum with unidirectional wave direction distribution, theoretical wave power spectrum with theoretical direction distribution, theoretical wave power spectrum with measured direction distribution, measured wave power spectrum with unidirectional wave direction distribution, measured wave power spectrum with theoretical direction distribution, and measured wave power spectrum with measured direction distribution; wherein, the measured wave power spectrum with measured direction distribution is used as the reference operating condition.

[0066] The power gain or loss metric is:

[0067] ;

[0068] Where G is the power gain or loss, P i To estimate the power under the combined operating conditions to be compared, P r This is the power estimate under the baseline operating condition;

[0069] The array power estimation results include the sea state average power matrix, annual power generation, and power estimation deviations under different wave power spectrum and direction distribution conditions.

[0070] The working principle of the wave energy device array power estimation output method based on wave direction spectrum analysis in this invention embodiment is as follows: it can simultaneously consider wave power spectrum type, wave direction distribution, array hydrodynamic coupling effect, and measured sea state characteristics of the target sea area, avoiding power estimation deviation caused by only using theoretical wave spectrum and unidirectional wave assumption, thereby more accurately reflecting the power capture capability of wave energy device array under multidirectional irregular wave conditions in the target sea area. It provides a basis for wave energy device array engineering design, array layout, power prediction, and economic evaluation, so as to solve the problem that theoretical wave spectrum and unidirectional wave assumption in existing array power estimation are prone to estimation deviation and difficult to reflect the real multidirectional wave environment, and improve the reliability of wave energy device array power estimation under real multidirectional sea conditions.

[0071] In the overall scheme, the estimation output method includes the following steps:

[0072] The wave energy device array to be estimated is determined, the spatial arrangement and spacing parameters of each wave energy device are determined, and measured wave data of the target sea area are obtained; the wave energy device is a heave-type point absorption wave energy device, and the array is composed of multiple wave energy devices.

[0073] Establish a hydrodynamic motion model and wave climate database corresponding to the wave energy device array. Incorporate the hydrostatic restoring force, radiation force, PTO force and wave excitation force on the wave energy device into the array dynamic equation. Referring to the hydrodynamic coupling between different wave energy devices in the array, perform state space approximation on the convolution integral term in the radiation force to construct a continuous state space model of the wave energy device array.

[0074] A multi-directional wave excitation model is constructed, which combines the one-dimensional wave power spectrum with the wave direction distribution function to form a wave direction spectrum, and generates a multi-directional irregular wave field based on the wave direction spectrum.

[0075] Acquire measured wave climate data for the target sea area. The measured wave climate data includes at least significant wave height, energy period or peak period, wave power spectrum and wave direction distribution information. Establish a sea state statistical matrix based on significant wave height and energy period.

[0076] The PTO damping parameters of the wave energy device are determined according to the dominant frequency or peak period of each sea state, so that the wave energy device can obtain a higher average capture power under the corresponding sea state.

[0077] The power estimation results under different combined operating conditions are compared with the power estimation results under the benchmark operating condition. The power gain or loss index is calculated to obtain the degree of influence of wave power spectrum type and direction distribution on the power estimation results of wave energy device array, and the array power estimation results are output.

[0078] Specifically, in this embodiment of the invention, linear arrays, equilateral triangular arrays, and square arrays are used as candidate array arrangements. The capture power of different array arrangements under the same sea conditions is estimated, and a hydrodynamic motion model of the wave energy device array is established. The heave motion equation of any single floating body in the array can be expressed as:

[0079] ;

[0080] Among them, f *i b (t), f *i r (t), f *i pto (t), f *i e (t) represents the hydrostatic restoring force, radiation force, PTO force, and wave excitation force on the corresponding buoy, respectively;

[0081] The static restoring force can be expressed as:

[0082] ;

[0083] Where, k *i h x is the hydrostatic stiffness coefficient. *i (t) represents the undulating displacement of the pontoon;

[0084] The radiative force can be expressed as:

[0085] ;

[0086] Where, μ *i*j For the added mass at infinite frequency, This is the radiation impulse response function of the corresponding radiative force; The buoy's heave speed; It can be approximated in state-space form using MATLAB, thereby accelerating the computation;

[0087] The general state space can be obtained in the following form:

[0088] .

[0089] Where q is the radiative force state variable, and A r B r C r D r It is the state-space approximation matrix of the radiation force.

[0090] In this invention, PTO power can be simplified into a speed proportional damping system controlled by the following factors:

[0091] ;

[0092] Among them, f i and S i For the float and reaction plate of the wave energy device, C pto is the damping coefficient.

[0093] Furthermore, by integrating hydrodynamics and PTO forces into a system using a multibody hydrodynamic coupling framework, a unified state-space representation of the two-body point absorber array is derived, with the specific equations as follows:

[0094] ;

[0095] in:

[0096] ;

[0097] ;

[0098] M is the combined inertia matrix, K h is the hydrostatic regenerative coefficient.

[0099] In a continuous-time array, the overall dynamics can be discretized under the assumption of piecewise constant excitation, resulting in a discrete state space, expressed as:

[0100] ;

[0101] in:

[0102] ;

[0103] The wave direction spectrum is represented as:

[0104] ;

[0105] Wherein, S(ω) is the one-dimensional wave power spectrum, D(θ) is the wave direction distribution function, ω is the wave angular frequency, and θ is the wave propagation direction; the wave direction spectrum can provide a complete statistical characterization of the sea state and quantify the frequency-related energy distribution and angular propagation characteristics.

[0106] Furthermore, D(θ) can be expressed as:

[0107] ;

[0108] The wave energy device includes a cylindrical float and an underwater reaction column or reaction plate. The cylindrical float serves as a wave energy capture component, and the underwater reaction column or reaction plate serves as a stabilizing component, enabling the wave energy device to capture energy under heave freedom.

[0109] The frequency domain dynamics of a cylindrical float undergoing heave motion can be expressed as:

[0110] ;

[0111] After merging and simplifying, we get:

[0112] ;

[0113] in:

[0114] ;

[0115] The magnitude of the complex heave velocity is obtained by performing determinant analysis on the linear system in the formula:

[0116] ;

[0117] The complex amplitude of the relative velocity between the float and the floating body can be expressed as:

[0118] ;

[0119] After simplification, we get:

[0120] ;

[0121] in:

[0122] ;

[0123] Therefore, the instantaneous power of the wave energy device can be calculated as follows:

[0124] ;

[0125] In the formula, † is a complex conjugate operator, and through time integration, these conjugate terms satisfy:

[0126] ;

[0127] Therefore, the time-averaged power of a wave energy device can be expressed as:

[0128] ;

[0129] Therefore, we get:

[0130] ;

[0131] To obtain optimal PTO damping, the average power must satisfy:

[0132] ;

[0133] The optimal damping can be derived as follows:

[0134] ;

[0135] Substituting this into the average power formula, we can obtain the maximum capture power:

[0136] ;

[0137] A wave energy device array comprises multiple wave energy device units. The spatial arrangement and spacing parameters of each wave energy device unit are determined as needed. The spatial arrangement can be a linear array, an equilateral triangular array, a square array, or other regular or irregular array forms. The spacing between adjacent wave energy device units can be an integer multiple of the device's characteristic diameter, preferably three times the device's characteristic diameter, to balance sea area utilization and array hydrodynamic interaction.

[0138] Furthermore, the displacement, velocity, and radiation force state variables of each wave energy device unit in the array are combined into an array state vector to establish a continuous state-space model of the wave energy device array. This model can simultaneously reflect the hydrodynamic coupling between the float and the underwater reactor inside a single device, as well as the radiation, diffraction, and shielding effects between different wave energy device units in the array.

[0139] Subsequently, a multi-directional wave excitation model was constructed. Irregular waves in actual sea areas do not propagate in a single direction, but rather have distribution characteristics in both frequency and direction.

[0140] Subsequently, measured wave climate data for the target sea area were acquired. This measured wave climate data includes at least the significant wave height H. s Energy cycle T e Or peak period Tp The system collects wave power spectrum and wave direction distribution information; it also performs zoned statistical analysis on the measured wave climate data based on significant wave height and energy cycle, establishing a sea state statistical matrix. Specifically, the significant wave height is divided into several wave height intervals, and the energy cycle is divided into several cycle intervals, and statistical analysis is performed on each wave power spectrum and wave direction distribution. s and T e The frequency or probability of occurrence of sea state intervals is used to obtain the sea state distribution characteristics of the target sea area.

[0141] For each sea state interval, the corresponding measured wave power spectrum and measured wave direction distribution are extracted as the basic input for subsequent array power estimation. The sea state statistical matrix reflects the long-term wave resource distribution in the target sea area and provides probability weights for estimating the array's annual power generation.

[0142] Finally, we constructed combined operating conditions with different wave power spectra and different directional distribution functions, and compared the power estimation results under different combined operating conditions with the power estimation results under the benchmark operating condition.

[0143] For the combined operating conditions, there are theoretical wave power spectrum-unidirectional wave direction distribution, theoretical wave power spectrum-theoretical direction distribution, theoretical wave power spectrum-measured direction distribution, measured wave power spectrum-unidirectional wave direction distribution, measured wave power spectrum-theoretical direction distribution, and measured wave power spectrum-measured direction distribution; wherein, the measured wave power spectrum-measured direction distribution is used as the reference operating condition.

[0144] Calculate the array power estimate for each operating combination. Compare the power estimate under the compared operating combination with the power estimate under the baseline operating combination, and calculate the power gain or loss index.

[0145] ;

[0146] Where G is the power gain or loss, P i To estimate the power under the combined operating conditions to be compared, P r The power estimate is the value under the baseline operating condition; by calculating the power gain or loss index, the influence of different wave power spectrum types and different directional distribution forms on the power estimation results of the wave energy device array is obtained.

[0147] Furthermore, the output results include the array sea state average power matrix, annual power generation, power estimation deviation under different combinations of operating conditions, the influence of wave power spectrum type on power estimation, the influence of directional distribution on power estimation, and the final array power estimation results.

[0148] It is worth noting that the output estimation method in this invention can avoid the power estimation bias caused by using only theoretical spectrum and unidirectional wave assumptions. It can more accurately reflect the power capture capability of wave energy device arrays under multidirectional irregular wave conditions in the target sea area, and provide a basis for wave energy device array engineering design, array layout, power prediction and economic evaluation.

[0149] In summary, the wave energy device array power estimation output method based on wave direction spectrum analysis in this embodiment of the invention can simultaneously consider wave power spectrum type, wave direction distribution, array hydrodynamic coupling effect, and measured sea state characteristics of the target sea area. It avoids power estimation deviation caused by using only theoretical spectrum and unidirectional wave assumption, thus more accurately reflecting the power capture capability of wave energy device array under multidirectional irregular wave conditions in the target sea area. This provides a basis for wave energy device array engineering design, array layout, power prediction, and economic evaluation, solving the problem that theoretical spectrum and unidirectional wave assumption in existing array power estimation are prone to estimation deviation and are difficult to reflect the real multidirectional wave environment, thereby improving the reliability of wave energy device array power estimation under real multidirectional sea conditions.

[0150] The above specific embodiments should not be construed as limiting the scope of protection of the present invention. For those skilled in the art, any alternative improvements or modifications made to the embodiments of the present invention shall fall within the scope of protection of the present invention.

[0151] Any aspects of this invention not described in detail are well-known to those skilled in the art.

Claims

1. A method for estimating the output power of a wave energy device array based on wave direction spectrum analysis, characterized in that, The output estimation method includes the following steps: S1, determine the wave energy device array object to be estimated, determine the spatial arrangement and spacing parameters of each wave energy device, and obtain measured wave data of the target sea area; the wave energy device is a heave-type point absorption wave energy device, and the array is composed of multiple wave energy devices. S2. Establish the hydrodynamic motion model and wave climate database corresponding to the wave energy device array. Incorporate the hydrostatic restoring force, radiation force, PTO force and wave excitation force on the wave energy device into the array dynamic equation. Referring to the hydrodynamic coupling between different wave energy devices in the array, perform state space approximation on the convolution integral term in the radiation force to construct a continuous state space model of the wave energy device array. S3. Construct a multi-directional wave excitation model, combine the one-dimensional wave power spectrum with the wave direction distribution function to form a wave direction spectrum, and generate a multi-directional irregular wave field based on the wave direction spectrum. S4. Obtain measured wave climate data for the target sea area. The measured wave climate data includes at least significant wave height, energy cycle or peak period, wave power spectrum and wave direction distribution information. Establish a sea state statistical matrix based on significant wave height and energy cycle. S5. Determine the PTO damping parameters of the wave energy device according to the dominant frequency or peak period of each sea state, so that the wave energy device can obtain a higher average capture power under the corresponding sea state. S6 compares the power estimation results under different combined operating conditions with the power estimation results under the benchmark operating condition, calculates the power gain or loss index, obtains the degree of influence of wave power spectrum type and direction distribution on the power estimation results of wave energy device array, and outputs the array power estimation results.

2. The wave energy device array power estimation output method based on wave direction spectrum analysis according to claim 1, characterized in that, The equation of motion for any single floating body in the array can be expressed as: ； Among them, f *i b (t), f *i r (t), f *i pto (t), f *i e (t) represents the hydrostatic restoring force, radiation force, PTO force, and wave excitation force on the corresponding buoy, respectively; The static restoring force can be expressed as: ； Where, k *i h x is the hydrostatic stiffness coefficient. *i (t) represents the undulating displacement of the pontoon; The radiative force can be expressed as: ； Where, μ *i*j For the added mass at infinite frequency, This is the radiation impulse response function of the corresponding radiative force; The buoy's heave speed; It can be approximated in state-space form using MATLAB, thereby accelerating the computation; The general state space can be obtained in the following form: 。 3. The wave energy device array power estimation output method based on wave direction spectrum analysis according to claim 1, characterized in that: The wave energy device includes a cylindrical float and an underwater reaction column or reaction plate. The cylindrical float serves as a wave energy capture component, and the underwater reaction column or reaction plate serves as a stabilizing component, enabling the wave energy device to capture energy under heave freedom.

4. The wave energy device array power estimation output method based on wave direction spectrum analysis according to claim 1, characterized in that: The spatial arrangement includes linear arrays, equilateral triangular arrays, and square arrays; the spacing between adjacent wave energy devices in the array is an integer multiple of the device's characteristic diameter; The wave direction spectrum is represented as: ； Where S(ω) is the one-dimensional wave power spectrum, D(θ) is the wave direction distribution function, ω is the wave angular frequency, and θ is the wave propagation direction; The wave direction distribution function satisfies the normalization condition, making the total energy of the wave direction spectrum consistent with the total energy of the one-dimensional wave power spectrum.

5. The wave energy device array power estimation output method based on wave direction spectrum analysis according to claim 1, characterized in that: The combined operating conditions include theoretical wave power spectrum with unidirectional wave direction distribution, theoretical wave power spectrum with theoretical direction distribution, theoretical wave power spectrum with measured direction distribution, measured wave power spectrum with unidirectional wave direction distribution, measured wave power spectrum with theoretical direction distribution, and measured wave power spectrum with measured direction distribution; wherein, the measured wave power spectrum with measured direction distribution is used as the reference operating condition.

6. The wave energy device array power estimation output method based on wave direction spectrum analysis according to claim 1, characterized in that: The power gain or loss metric is: ； Where G is the power gain or loss, P i To estimate the power under the combined operating conditions to be compared, P r This is the power estimate under the baseline operating condition; The array power estimation results include the sea state average power matrix, annual power generation, and power estimation deviations under different wave power spectrum and direction distribution conditions.

Images