A method, device and electronic equipment for estimating power generation of a wave power plant

Abstract

The application provides a method, device and electronic equipment for estimating power generation of a wave power generation device, wherein the method comprises: calculating a wave height-period data distribution of a wave force time domain signal half wave of the wave power generation device; calculating a half wave energy corresponding to the wave force time domain signal half wave according to the wave height-period data distribution of the wave force time domain signal half wave; and calculating an average power of the wave force time domain signal according to the wave height-period data distribution of the wave force time domain signal half wave and the calculated half wave energy of the wave force time domain signal. The method, device and electronic equipment for estimating power generation of the wave power generation device provided by the application do not need to solve a large number of or large optimization problems of the wave force time domain signal, and do not need to calculate the whole process of system operation, thereby greatly reducing the calculation time.

Description

Technical Field

[0001] This application relates to the field of wave power generation technology, and more specifically, to a method, apparatus, and electronic equipment for estimating the power generation capacity of a wave power generation device. Background Technology

[0002] Currently, when designing wave power generation devices, it is necessary to determine the core parameters of the device (such as float shape, size, maximum float travel, motor and converter capacity, etc.) for the wave conditions of a specific sea area, so as to minimize the cost of the wave power generation device while maximizing its wave energy capture efficiency. The core issue here is how to estimate the wave energy capture efficiency of a wave power generation device with given parameters.

[0003] In related technologies, methods for estimating the average power of a wave power generation device with given parameters fall into two categories: one is time-domain simulation based on model predictive control, and the other is frequency-domain optimization based on pseudospectral methods. However, both methods require solving a large number of optimization problems, resulting in slow computation speeds and low processing efficiency for estimating the wave energy capture efficiency of a wave power generation device with given parameters. Summary of the Invention

[0004] To address the aforementioned problems, the purpose of this application is to provide a method, apparatus, and electronic device for estimating the power output of a wave power generation device.

[0005] In a first aspect, embodiments of this application provide a method for estimating the power generation capacity of a wave power generation device, including:

[0006] Calculate the wave height-period data distribution of the half-wave of the wave force time-domain signal of the wave power generator; wherein, the wave force time-domain signal is a sinusoidal time signal;

[0007] Based on the wave height-period data distribution of the wave force time domain signal half-wave, the half-wave energy corresponding to the wave force time domain signal half-wave is calculated.

[0008] The average power of the wave force time domain signal is calculated based on the wave height-period data distribution of the half-wave of the wave force time domain signal and the calculated half-wave energy of the wave force time domain signal.

[0009] Secondly, embodiments of this application also provide an apparatus for estimating the power generation capacity of a wave power generation device, comprising:

[0010] The first calculation module is used to calculate the wave height-period data distribution of the wave force time domain signal half-wave of the wave power generation device; wherein, the wave force time domain signal is a sinusoidal time signal;

[0011] The second calculation module is used to calculate the half-wave energy corresponding to the half-wave of the wave force time domain signal based on the wave height-period data distribution of the half-wave.

[0012] The third calculation module is used to calculate the average power of the wave force time domain signal based on the wave height-period data distribution of the half-wave of the wave force time domain signal and the calculated half-wave energy of the wave force time domain signal.

[0013] Thirdly, embodiments of this application also provide a computer-readable storage medium storing a computer program, which, when executed by a processor, performs the steps of the method described in the first aspect above.

[0014] Fourthly, embodiments of this application also provide an electronic device, the electronic device including a memory, a processor and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the processor using the steps of the method described in the first aspect above.

[0015] In the solutions provided by the first to fourth aspects of this application, the wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device is calculated; based on the wave height-period data distribution of the wave force time-domain signal half-wave, the half-wave energy corresponding to the half-wave of the wave force time-domain signal is calculated; based on the wave height-period data distribution of the wave force time-domain signal half-wave and the calculated half-wave energy of the wave force time-domain signal, the average power of the wave force time-domain signal is calculated. Compared with the related technologies that require solving a large number or large-scale optimization problems, this method divides the wave force time-domain signal of the wave power generation device into multiple half-waves and calculates the average power of the wave force time-domain signal by calculating the half-wave energy corresponding to each half-wave of the wave force time-domain signal separately. The entire calculation process does not require solving a large number or large-scale optimization problems of the wave force time-domain signal and does not require calculation of the entire system operation process, which greatly reduces the calculation time. At the same time, the proposed analytical calculation method is also faster than the numerical method, further accelerating the calculation efficiency of the average power.

[0016] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1 A flowchart is shown below illustrating a method for estimating the power output of a wave power generation device according to Embodiment 1 of this application;

[0019] Figure 2 This invention provides a schematic diagram of the structure of a device for estimating the power output of a wave power generation device according to Embodiment 2 of this application;

[0020] Figure 3 A schematic diagram of the structure of an electronic device provided in Embodiment 3 of this application is shown. Detailed Implementation

[0021] In the description of this application, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are used only for the convenience of describing this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this application.

[0022] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this application, "multiple" means two or more, unless otherwise explicitly specified.

[0023] In this application, unless otherwise expressly specified and limited, the terms "installation," "connection," "linking," and "fixing," etc., should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; and they can refer to the internal connection between two components. Those skilled in the art can understand the specific meaning of the above terms in this application according to the specific circumstances.

[0024] Currently, when designing wave power generation devices, it is necessary to determine the core parameters of the device (such as float shape, size, maximum float travel, motor and converter capacity, etc.) for the wave conditions of a specific sea area, so as to minimize the cost of the wave power generation device while maximizing its wave energy capture efficiency. The core issue here is how to estimate the wave energy capture efficiency of a wave power generation device with given parameters.

[0025] In related technologies, methods for estimating the average power of a wave power generation device with given parameters fall into two categories: one is time-domain simulation based on model predictive control, and the other is frequency-domain optimization based on pseudospectral methods. However, both methods require solving a large number of optimization problems, resulting in slow computation speeds and low processing efficiency for estimating the wave energy capture efficiency of a wave power generation device with given parameters.

[0026] Based on this, this embodiment proposes a method, apparatus, and electronic device for estimating the power generation of a wave power generation device. The method involves calculating the wave height-period data distribution of the half-wave of the wave force time-domain signal from the wave power generation device; calculating the half-wave energy corresponding to the half-wave of the wave force time-domain signal based on the wave height-period data distribution; and calculating the average power of the wave force time-domain signal based on the wave height-period data distribution of the half-wave of the wave force time-domain signal and the calculated half-wave energy. The wave force time-domain signal of the wave power generation device is divided into multiple half-waves, and the average power of the wave force time-domain signal is calculated by individually calculating the half-wave energy corresponding to each half-wave. The entire calculation process does not require solving a large number of optimization problems related to the wave force time-domain signal, nor does it require calculation of the entire system operation process, greatly reducing the calculation time. Simultaneously, the proposed analytical calculation method is faster than numerical methods, further accelerating the calculation efficiency of the average power.

[0027] In the following embodiments, the terms "half-wave energy corresponding to half-wave of the wave force time-domain signal" and "half-wave energy of the wave force time-domain signal" have the same meaning.

[0028] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and embodiments.

[0029] Example 1

[0030] This embodiment proposes a method for estimating the power output of a wave power generation device, with the server as the execution entity.

[0031] See Figure 1 The flowchart shown illustrates a method for estimating the power output of a wave power generation device. This embodiment proposes a method for estimating the power output of a wave power generation device, including the following specific steps:

[0032] Step 100: Calculate the wave height-period data distribution of the wave force time domain signal half-wave of the wave power generation device; wherein, the wave force time domain signal is a sinusoidal time signal.

[0033] In step 100 above, the wave height-period data distribution of the wave force time domain signal half-wave includes: the correspondence between the wave height and period of the wave force time domain signal half-wave, and the joint probability distribution of the wave height and period of the wave force time domain signal half-wave.

[0034] Given that the wave height-period data distribution of the wave force time-domain signal half-wave corresponds to the wave height-period of the wave force time-domain signal half-wave, in order to calculate the wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device, the following steps (1) to (7) can be performed:

[0035] (1) Obtain the wave excitation coefficient of the float, the energy spectrum of the wave, and the fundamental frequency of the wave as preset;

[0036] (2) Based on the fundamental frequency, obtain the frequency points of multiples of the fundamental frequency i, and use the frequency points of multiples of the fundamental frequency i to sample the energy spectrum to obtain the value of the energy spectrum at multiples of the fundamental frequency i; where i∈[0,1,2…n];

[0037] (3) The wave excitation coefficient of the float is sampled using the frequency point of the i-fold fundamental frequency to obtain the wave excitation coefficient of the float at the frequency point of the i-fold fundamental frequency; wherein, the wave excitation coefficient of the float at the frequency point of the i-fold fundamental frequency includes: the amplitude and the first phase at the frequency point of the i-fold fundamental frequency.

[0038] (4) Randomly assign the second phase to the frequency points of i multiples of the fundamental frequency between 0 and 2π;

[0039] (5) Calculate the wave force time-domain signal w(t) using the following formula:

[0040]

[0041] Where n represents the number of frequency points that are multiples of the fundamental frequency; |F e (jiω0)| represents the first phase; S(iω0) represents the value of the energy spectrum at a frequency point that is a multiple of the fundamental frequency; ω0 represents the fundamental frequency; φ i Indicates the second phase; t represents the time of the wave force time-domain signal; ∠F e (jiω0) represents the amplitude at a frequency point that is a multiple of the fundamental frequency;

[0042] (6) Detect the zero-crossing point of the wave force time domain signal to obtain the half-wave of each wave force time domain signal in the wave force time domain signal;

[0043] (7) Extract the wave height and period of each wave force time domain signal half wave respectively, and obtain the correspondence between the wave height and period of each wave force time domain signal half wave.

[0044] In step (1) above, the wave excitation coefficient of the float, the energy spectrum of the wave, and the pre-set fundamental frequency of the wave are pre-cached in the server.

[0045] In step (3) above, the wave excitation coefficient of the float at a frequency point that is a multiple of the fundamental frequency is a complex number. Therefore, the wave excitation coefficient of the float at a frequency point that is a multiple of the fundamental frequency includes the real part and the imaginary part of the complex number, respectively, for the amplitude and the first phase at the frequency point that is a multiple of the fundamental frequency.

[0046] In step (5) above, see Figure 2 The diagram shows a wave force time-domain signal w(t), where w(t) is a sinusoidal time signal.

[0047] In step (6) above, the zero-crossing point of the wave force time domain signal is detected, and the half-wave signal between adjacent zero-crossing points is determined as a half-wave of the wave force time domain signal in a wave force time domain signal.

[0048] In step (7) above, see Figure 3 The diagram illustrates the correspondence between wave height and period of N sets of wave force time-domain signal half-waves. The wave height-period data distribution of the wave force time-domain signal half-waves is defined as follows: N sets of wave height-period correspondences for wave force time-domain signal half-waves. The wave height-period correspondence of each set of wave force time-domain signal half-waves in the N sets of wave force time-domain signal half-wave correspondences is expressed as (W... r D r ), r∈1,…,N.

[0049] Of course, in addition to the process of calculating the correspondence between wave height and period of half wave of wave force time domain signal in steps (1) to (7) above, this embodiment can also use any existing scheme to calculate the correspondence between wave height and period of half wave of wave force time domain signal, which will not be elaborated here.

[0050] When the wave height-period data distribution of the wave force time-domain signal half-wave is the joint probability distribution of the wave force time-domain signal half-wave, in order to calculate the wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device, the following steps (11) to (12) can be performed:

[0051] (11) Obtain the energy spectrum of the wave and the wave excitation coefficient of the float;

[0052] (12) The energy spectrum of the waves and the wave excitation coefficient of the float were processed using the ocean wave data analysis software (WAFO) to obtain the joint probability distribution of the wave height-period of the half-wave of the wave force time domain signal.

[0053] Alternatively, when the wave height-period data distribution of the wave force time-domain signal half-wave is a joint probability distribution of the wave force time-domain signal half-wave, the calculation of the wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device includes the following steps (21) to (23):

[0054] (21) Obtain the energy spectrum of the wave and the wave excitation coefficient of the float;

[0055] (22) Calculate the energy spectrum characteristic parameters based on the energy spectrum of the waves and the wave excitation coefficient of the float;

[0056] (23) Based on the energy spectrum characteristic parameters, the LH probability distribution model is obtained, and the joint probability distribution of wave height-period of the wave force time domain signal half wave is calculated according to the LH probability distribution model.

[0057] In step (12) above, the specific process of using WAFO software to process the energy spectrum of the wave and the wave excitation coefficient of the float to obtain the joint probability distribution of the wave height-period of the half-wave of the wave force time domain signal is existing technology and will not be described in detail here.

[0058] In steps (22) to (23) above, the energy spectrum characteristic parameters m0, m1, and m2 are calculated as follows:

[0059]

[0060] Then calculate the parameter v and the coefficient c. LH as follows:

[0061]

[0062]

[0063] Therefore, the LH probability distribution model is shown in the following formula:

[0064]

[0065] Where ω represents the angular frequency; S(ω) represents the energy spectrum; F e(jω) represents the buoyancy excitation force coefficient; p(W,D) represents the joint probability distribution of wave height-period of the half-wave of the wave force time domain signal; W represents wave height; D represents period.

[0066] The three cases described above for calculating the wave height-period data distribution of the half-wave of the wave force time-domain signal from the wave power generation device all assume that an accurate wave energy spectrum can be obtained. However, only the effective wave height H can be obtained. s and peak period T p In such cases, a wave spectrum model is needed to obtain the energy spectrum of the wave.

[0067] Specifically, in order to obtain only the effective wave height H of the wave s and peak period T p In this case, to obtain the energy spectrum of the wave, the following steps (1) to (2) can be performed:

[0068] (1) Obtain the effective wave height and peak period of the wave;

[0069] (2) Input the effective wave height and peak period of the wave into the spectrum model for modeling to obtain the energy spectrum of the wave.

[0070] In step (1) above, the effective wave height and peak period of the wave are pre-stored in the server.

[0071] In step (2) above, the spectrum model includes, but is not limited to, the JONSWAP model and the Bretschneider model.

[0072] For example, by inputting the effective wave height and peak period of the wave into the Bretschneider model for modeling, the energy spectrum of the wave is obtained as follows:

[0073]

[0074] Where, ω p ω is the peak frequency. p =2π / T p .

[0075] The specific process of inputting the effective wave height and peak period of the wave into the spectrum model to obtain the energy spectrum of the wave is existing technology and will not be described in detail here.

[0076] Using the information from steps (1) to (2) above, only the effective wave height H can be obtained. s and peak period T pIn this case, after obtaining the energy spectrum of the wave, we can return to the three cases mentioned above for calculating the wave height-period data distribution of the half-wave of the wave force time domain signal of the wave power generator, and calculate the wave height-period data distribution of the half-wave of the wave force time domain signal of the wave power generator.

[0077] After calculating the wave height-period data distribution of the wave force time domain signal half-wave through the above step 100, the half-wave energy corresponding to the wave force time domain signal half-wave is calculated through the following step 102.

[0078] Step 102: Calculate the half-wave energy corresponding to the half-wave of the wave force time domain signal based on the wave height-period data distribution of the half-wave.

[0079] In step 102 above, the half-wave energy corresponding to the half-wave of the wave force time-domain signal can be calculated in the following two ways.

[0080] In the first case, in order to calculate the half-wave energy corresponding to the half-wave of the wave force time-domain signal, the following steps (1) to (3) can be performed:

[0081] (1) Obtain the maximum displacement of the wave power generation device;

[0082] (2) For the wave height-period data distribution of any wave force time-domain signal half-wave, the constraint coefficient α is calculated using the following formula:

[0083]

[0084] Among them, R eq Z represents the damping coefficient that matches the period recorded in the wave height-period data distribution of the wave force time-domain signal half-wave for calculating the current half-wave energy; m Indicates the maximum displacement of the wave generator; W represents the wave height; D represents the period.

[0085] (3) The half-wave energy of the wave force time-domain signal is calculated using the following formula:

[0086]

[0087] Where E(W, D) represents the half-wave energy of the wave force time-domain signal.

[0088] In step (1) above, the maximum displacement of the wave generator is pre-cached in the server.

[0089] In step (2) above, the server has a pre-stored correspondence between period and damping coefficient; then, in order to obtain the damping coefficient, the server queries the correspondence between period and damping coefficient to find the damping coefficient that matches the period recorded in the wave height-period data distribution of the wave force time domain signal half-wave for calculating the current half-wave energy.

[0090] Here, when the wave height-period data distribution of the wave force time-domain signal half-wave is N sets of wave height-period correspondences for the wave force time-domain signal half-wave, the wave height-period correspondence (W) for each set of wave force time-domain signal half-wave is... r D r The half-wave energy of the wave force time-domain signal calculated is expressed as E(W). r D r ), r∈1,...,N.

[0091] Given the wave height-period data distribution of a half-wave in the wave force time-domain signal as p(W, D), p(W, D) is essentially a binary function, represented as a two-dimensional array (table) in the server. The meaning of p(W, D) is the probability (density) that for any given half-wave, its wave height is W and its period is D. Therefore, given p(W, D), the server can substitute W and D within the selected probability range from p(W, D) into the formula for calculating the constraint coefficient α to calculate the constraint coefficient α.

[0092] In the second case, the calculation of the half-wave energy of the wave force time-domain signal by utilizing the wave height-period data distribution of the half-wave of the wave force time-domain signal further includes the following steps (1) to (2):

[0093] (1) Based on the wave height and period recorded in the wave height-period data distribution of the wave force time domain signal half wave, a regular wave is obtained whose wave height is a preset multiple of the wave height of the wave force time domain signal half wave and whose period is a preset multiple of the period of the wave force time domain signal half wave.

[0094] (2) Input the obtained regular wave into WecOptTool to calculate the half-wave energy of the wave force time domain signal.

[0095] In step (1) above, the preset quantity multiple is 2.

[0096] The specific process of obtaining a regular wave whose wave height is a preset multiple of the wave height of the wave force time domain signal half-wave and whose period is a preset multiple of the period of the wave force time domain signal half-wave, based on the wave height and period recorded in the wave height-period data distribution of the wave force time domain signal half-wave, is existing technology and will not be elaborated here.

[0097] Given p(W,D), a regular wave can be obtained by selecting W and D within the probability range of p(W,D), where the wave height is a preset multiple of the wave height of the half-wave of the wave force time domain signal and the period is a preset multiple of the period of the half-wave of the wave force time domain signal.

[0098] In step (2) above, the obtained regular wave is input into WecOptTool to calculate the half-wave energy of the wave force time domain signal.

[0099] Step 104: Calculate the average power of the wave force time domain signal based on the wave height-period data distribution of the half-wave of the wave force time domain signal and the calculated half-wave energy of the wave force time domain signal.

[0100] Specifically, in order to calculate the average power of the wave force time-domain signal, step 104 above can be performed by executing the following steps (1) to (2):

[0101] (1) When the wave height-period data distribution of the half-wave of the wave force time-domain signal adopts the correspondence between the wave height and period of the half-wave of the wave force time-domain signal, the average power P of the wave force time-domain signal is calculated by the following formula:

[0102]

[0103] Where N represents the number of wave height-period correspondences for the half-wave of the wave force time-domain signal; E(W r D r D represents the half-wave energy of the wave force time-domain signal calculated from the correspondence between wave height and period of the r-th group of wave force time-domain signals; r The period represents the correspondence between the wave height and period of the r-th group of wave force time-domain signals.

[0104] (2) When the wave height-period data distribution of the half-wave of the wave force time-domain signal adopts the joint probability distribution of the wave height-period of the half-wave of the wave force time-domain signal, the average power P of the wave force time-domain signal is calculated by the following formula:

[0105]

[0106] Where p(W,D) represents the joint probability distribution of wave height-period of half-wave of wave force time-domain signal; E(W,D) represents the half-wave energy of wave force time-domain signal calculated from the joint probability distribution of wave height-period of half-wave of wave force time-domain signal.

[0107] In step (1) above, This represents the total energy of N half-waves in the wave force time-domain signal. The total time of N half-waves in the wave force time domain signal is represented by the sum of the two values, which gives the average power P.

[0108] In step (2) above, the result of the integral of ∫∫p(W,D)E(W,D)dWdD represents the total energy of the N half waves in the wave force time domain signal, and the result of the integral of ∫∫p(W,D)DdWdD represents the total time of the N half waves in the wave force time domain signal. The difference between the two is the average power P.

[0109] In summary, this embodiment proposes a method for estimating the power output of a wave power generation device. It calculates the wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device. Based on this data distribution, the half-wave energy corresponding to each half-wave is calculated. Finally, based on the wave height-period data distribution and the calculated half-wave energy, the average power of the wave force time-domain signal is calculated. Compared to related technologies that require solving numerous or large-scale optimization problems, this method divides the wave force time-domain signal of the wave power generation device into multiple half-waves and calculates the average power of the wave force time-domain signal by individually calculating the half-wave energy corresponding to each half-wave. The entire calculation process does not require solving numerous or large-scale optimization problems related to the wave force time-domain signal and does not require calculations throughout the entire system operation process, significantly reducing computation time. Furthermore, the proposed analytical calculation method is faster than numerical methods, further accelerating the calculation efficiency of the average power.

[0110] Example 2

[0111] This embodiment proposes an apparatus for estimating the power output of a wave power generation device, used to execute the method for estimating the power output of a wave power generation device proposed in Embodiment 1 above.

[0112] See Figure 2 The diagram shows a structural schematic of a device for estimating the power output of a wave power generation device. This embodiment proposes a device for estimating the power output of a wave power generation device, comprising:

[0113] The first calculation module 200 is used to calculate the wave height-period data distribution of the wave force time domain signal half-wave of the wave power generation device; wherein, the wave force time domain signal is a sinusoidal time signal;

[0114] The second calculation module 202 is used to calculate the half-wave energy corresponding to the half-wave of the wave force time domain signal based on the wave height-period data distribution of the half-wave.

[0115] The third calculation module 204 is used to calculate the average power of the wave force time domain signal based on the wave height-period data distribution of the half-wave of the wave force time domain signal and the calculated half-wave energy of the wave force time domain signal.

[0116] In one embodiment, the wave height-period data distribution of the wave force time-domain signal half-wave includes: the correspondence between the wave height and period of the wave force time-domain signal half-wave.

[0117] Specifically, the first calculation module is used for:

[0118] Obtain the wave excitation coefficient of the float, the energy spectrum of the wave, and the pre-set fundamental frequency of the wave;

[0119] Based on the fundamental frequency, frequency points of multiples of the fundamental frequency are obtained, and the energy spectrum is sampled using the frequency points of multiples of the fundamental frequency to obtain the value of the energy spectrum at multiples of the fundamental frequency; where i∈[0,1,2…n].

[0120] The wave excitation coefficient of the float is sampled using the frequency point of the i-th multiple of the fundamental frequency to obtain the wave excitation coefficient of the float at the frequency point of the i-th multiple of the fundamental frequency; wherein, the wave excitation coefficient of the float at the frequency point of the i-th multiple of the fundamental frequency includes: the amplitude and the first phase at the frequency point of the i-th multiple of the fundamental frequency.

[0121] Between 0 and 2π, a second phase is randomly assigned to frequency points that are multiples of the fundamental frequency (i).

[0122] The wave force time-domain signal w(t) is calculated using the following formula:

[0123]

[0124] Where n represents the number of frequency points that are multiples of the fundamental frequency; |F e (jiω0)| represents the first phase; S(iω0) represents the value of the energy spectrum at a frequency point that is a multiple of the fundamental frequency; ω0 represents the fundamental frequency; φ i Indicates the second phase; t represents the time of the wave force time-domain signal; ∠F e (jiω0) represents the amplitude at a frequency point that is a multiple of the fundamental frequency;

[0125] The zero-crossing point of the wave force time-domain signal is detected to obtain the half-wave of each wave force time-domain signal in the wave force time-domain signal;

[0126] The wave height and period of each wave force time domain signal half wave are extracted to obtain the correspondence between wave height and period of each wave force time domain signal half wave.

[0127] In summary, this embodiment proposes a device for estimating the power output of a wave power generation device. It calculates the wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device. Based on this data distribution, the half-wave energy corresponding to each half-wave is calculated. Finally, based on the wave height-period data distribution and the calculated half-wave energy, the average power of the wave force time-domain signal is calculated. Compared to related technologies that require solving numerous or large-scale optimization problems, this method divides the wave force time-domain signal of the wave power generation device into multiple half-waves and calculates the average power of the wave force time-domain signal by individually calculating the half-wave energy corresponding to each half-wave. The entire calculation process does not require solving numerous or large-scale optimization problems related to the wave force time-domain signal and does not require calculations throughout the entire system operation process, significantly reducing computation time. Furthermore, the proposed analytical calculation method is faster than numerical methods, further accelerating the calculation efficiency of the average power.

[0128] Example 3

[0129] This embodiment proposes a computer-readable storage medium storing a computer program. When the computer program is run by a processor, it executes the steps of the method for estimating the power generation of a wave power generation device as described in Embodiment 1 above. For a detailed implementation, please refer to Method Embodiment 1, which will not be repeated here.

[0130] In addition, see Figure 3 The diagram shows the structure of an electronic device. This embodiment also proposes an electronic device, which includes a bus 51, a processor 52, a transceiver 53, a bus interface 54, a memory 55, and a user interface 56. The electronic device includes a memory 55.

[0131] In this embodiment, the electronic device further includes: one or more programs stored in the memory 55 and executable on the processor 52, configured to be executed by the processor to perform the one or more programs for the following steps (1) to (3):

[0132] (1) Calculate the wave height-period data distribution of the wave force time domain signal half-wave of the wave power generation device; wherein, the wave force time domain signal is a sinusoidal time signal;

[0133] (2) Based on the wave height-period data distribution of the wave force time domain signal half-wave, the half-wave energy corresponding to the wave force time domain signal half-wave is calculated.

[0134] (3) The wave height-period data distribution of the half-wave of the wave force time domain signal and the half-wave energy of the wave force time domain signal are calculated, and the average power of the wave force time domain signal is calculated.

[0135] Transceiver 53 is used to receive and send data under the control of processor 52.

[0136] The bus architecture (represented by bus 51) can include any number of interconnected buses and bridges, linking various circuits including one or more processors represented by processor 52 and memory represented by memory 55. Bus 51 can also link various other circuits such as peripheral devices, voltage regulators, and power management circuits, which are well known in the art and therefore will not be further described in this embodiment. Bus interface 54 provides an interface between bus 51 and transceiver 53. Transceiver 53 can be a single element or multiple elements, such as multiple receivers and transmitters, providing a unit for communicating with various other devices over a transmission medium. For example, transceiver 53 receives external data from other devices. Transceiver 53 is used to transmit data processed by processor 52 to other devices. Depending on the nature of the computing system, a user interface 56 may also be provided, such as a keypad, display, speaker, microphone, or joystick.

[0137] Processor 52 is responsible for managing bus 51 and general processing, such as running general-purpose operating system 551 as described above. Memory 55 can be used to store data used by processor 52 during operation.

[0138] Optionally, the processor 52 may be, but is not limited to, a central processing unit, a microcontroller, a microprocessor, or a programmable logic device.

[0139] It is understood that the memory 55 in the embodiments of this application can be volatile memory or non-volatile memory, or may include both volatile and non-volatile memory. The non-volatile memory can be read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), or flash memory. The volatile memory can be random access memory (RAM), which is used as an external cache. By way of example, but not limitation, many forms of RAM are available, such as Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced Synchronous DRAM (ESDRAM), Synchlink DRAM (SLDRAM), and Direct Rambus RAM (DRRAM). The memory 55 of the systems and methods described in this embodiment is intended to include, but is not limited to, these and any other suitable types of memory.

[0140] In some implementations, memory 55 stores elements such as executable modules or data structures, or subsets thereof, or extended sets thereof: operating system 551 and application programs 552.

[0141] The operating system 551 includes various system programs, such as the framework layer, core library layer, and driver layer, used to implement various basic business functions and handle hardware-based tasks. The application program 552 includes various applications, such as a media player and a browser, used to implement various application functions. Programs implementing the methods of the embodiments of this application can be included in the application program 552.

[0142] In summary, this embodiment proposes a computer-readable storage medium and electronic device that calculates the wave height-period data distribution of a half-wave of the wave force time-domain signal from a wave power generation device. Based on this data distribution, the half-wave energy corresponding to the half-wave of the wave force time-domain signal is calculated. Furthermore, based on the wave height-period data distribution and the calculated half-wave energy, the average power of the wave force time-domain signal is calculated. Compared to related technologies that require solving numerous or large-scale optimization problems, this method divides the wave force time-domain signal from the wave power generation device into multiple half-waves and calculates the average power of the wave force time-domain signal by individually calculating the half-wave energy corresponding to each half-wave. The entire calculation process does not require solving numerous or large-scale optimization problems related to the wave force time-domain signal and does not require calculation of the entire system operation process, significantly reducing computation time. Simultaneously, the proposed analytical calculation method is faster than numerical methods, further accelerating the calculation efficiency of the average power.

[0143] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for estimating the power output of a wave power generation device, characterized in that, include: Calculate the wave height-period data distribution of the half-wave of the wave force time-domain signal of the wave power generator; wherein, the wave force time-domain signal is a sinusoidal time signal; Based on the wave height-period data distribution of the wave force time domain signal half-wave, the half-wave energy corresponding to the wave force time domain signal half-wave is calculated. Based on the wave height-period data distribution of the half-wave of the wave force time-domain signal and the calculated half-wave energy of the wave force time-domain signal, the average power of the wave force time-domain signal is calculated. The step of calculating the half-wave energy corresponding to the half-wave of the wave force time-domain signal based on the wave height-period data distribution of the half-wave includes: Obtain the maximum displacement of the wave power generation device; For the wave height-period data distribution of any half-wave of a wave force time-domain signal, the constraint coefficients are determined using the following formula. Perform the calculation: ； in, The damping coefficient represents the period matching recorded in the wave height-period data distribution of the wave force time-domain signal half-wave for which the current half-wave energy is calculated; This indicates the maximum displacement of the wave generator; Represents wave height; Indicates period; The half-wave energy of the wave force time-domain signal is calculated using the following formula: ； in, This represents the half-wave energy of the wave force time-domain signal.

2. The method according to claim 1, characterized in that, The wave height-period data distribution of the wave force time domain signal half-wave includes: the correspondence between the wave height and period of the wave force time domain signal half-wave. The wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device includes: Obtain the wave excitation coefficient of the float, the energy spectrum of the wave, and the pre-set fundamental frequency of the wave; Based on the fundamental frequency, frequency points of multiples of the fundamental frequency are obtained, and the energy spectrum is sampled using the frequency points of multiples of the fundamental frequency to obtain the value of the energy spectrum at multiples of the fundamental frequency; where i∈[0,1,2…n]. The wave excitation coefficient of the float is sampled using the frequency point of the i-th multiple of the fundamental frequency to obtain the wave excitation coefficient of the float at the frequency point of the i-th multiple of the fundamental frequency; wherein, the wave excitation coefficient of the float at the frequency point of the i-th multiple of the fundamental frequency includes: the amplitude and the first phase at the frequency point of the i-th multiple of the fundamental frequency. Between 0 and 2π, a second phase is randomly assigned to frequency points that are multiples of the fundamental frequency (i). The wave force time-domain signal is calculated using the following formula. : ； in, This represents the number of frequency points that are multiples of the fundamental frequency; Indicates the first phase; This represents the numerical value of the energy spectrum at frequency points that are multiples of the fundamental frequency (i). Indicates the fundamental frequency; Indicates the second phase; This represents the time domain signal of wave force; This represents the amplitude at a frequency point that is a multiple of the fundamental frequency (i). The zero-crossing point of the wave force time-domain signal is detected to obtain the half-wave of each wave force time-domain signal in the wave force time-domain signal; The wave height and period of each wave force time domain signal half wave are extracted to obtain the correspondence between wave height and period of each wave force time domain signal half wave.

3. The method according to claim 2, characterized in that, The wave height-period data distribution of the wave force time domain signal half-wave includes: the joint probability distribution of the wave height-period of the wave force time domain signal half-wave; The wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device includes: Obtain the energy spectrum of the waves and the wave excitation coefficient of the buoy; The energy spectrum of the wave and the wave excitation coefficient of the float were processed using WAFO software to obtain the joint probability distribution of the wave height-period of the half-wave of the wave force time domain signal. or, The wave height-period data distribution of the wave force time-domain signal half-wave of the wave power generation device includes: Obtain the energy spectrum of the waves and the wave excitation coefficient of the buoy; Calculate the characteristic parameters of the energy spectrum based on the energy spectrum of the wave and the wave excitation coefficient of the buoy; Based on the energy spectrum characteristic parameters, the LH probability distribution model is obtained, and the joint probability distribution of wave height-period of the half-wave of the wave force time domain signal is calculated according to the LH probability distribution model.

4. The method according to claim 3, characterized in that, The wave height-period data distribution of the half-wave of the wave force time-domain signal is defined as follows: The correspondence between wave height and period of half-wave in the time-domain signal of wave force. The correspondence between wave height and period for each half-wave of a group of wave force time-domain signals is expressed as follows: , , The half-wave energy of the wave force time-domain signal, calculated from the correspondence between wave height and period of each group of wave force time-domain signals, is expressed as: , ; The step of calculating the average power of the wave force time-domain signal based on the wave height-period data distribution of the half-wave of the wave force time-domain signal and the calculated half-wave energy of the wave force time-domain signal includes: When the wave height-period data distribution of the half-wave of the wave force time-domain signal adopts the correspondence between wave height and period of the half-wave of the wave force time-domain signal, the average power of the wave force time-domain signal is calculated by the following formula. : ； in, The number of correspondences between wave height and period for half-waves of wave force time-domain signals; Indicates the first The half-wave energy of the wave force time-domain signal is calculated from the correspondence between wave height and period of the half-wave of the wave force time-domain signal. Indicates the first The period in the correspondence between wave height and period of a half-wave of a group of wave force time-domain signals; When the wave height-period data distribution of the half-wave of the wave force time-domain signal is adopted using the joint probability distribution of wave height-period of the half-wave of the wave force time-domain signal, the average power of the wave force time-domain signal is calculated by the following formula. : ； in, This represents the joint probability distribution of wave height and period of a half-wave of a wave force time-domain signal. The half-wave energy of the wave force time-domain signal is calculated when the joint probability distribution of wave height and period of the half-wave of the wave force time-domain signal is expressed.

5. A device for estimating the power output of a wave power generation device, characterized in that, include: The first calculation module is used to calculate the wave height-period data distribution of the half-wave of the wave force time-domain signal of the wave power generation device; wherein, the wave force time-domain signal is a sinusoidal time signal; The second calculation module is used to calculate the half-wave energy corresponding to the half-wave of the wave force time domain signal based on the wave height-period data distribution of the half-wave. The third calculation module is used to calculate the average power of the wave force time domain signal based on the wave height-period data distribution of the half-wave of the wave force time domain signal and the calculated half-wave energy of the wave force time domain signal. The step of calculating the half-wave energy corresponding to the half-wave of the wave force time-domain signal based on the wave height-period data distribution of the half-wave includes: Obtain the maximum displacement of the wave power generation device; For the wave height-period data distribution of any half-wave of a wave force time-domain signal, the constraint coefficients are determined using the following formula. Perform the calculation: ； in, The damping coefficient represents the period matching recorded in the wave height-period data distribution of the wave force time-domain signal half-wave for which the current half-wave energy is calculated; This indicates the maximum displacement of the wave generator; Represents wave height; Indicates period; The half-wave energy of the wave force time-domain signal is calculated using the following formula: ； in, This represents the half-wave energy of the wave force time-domain signal.

6. The apparatus according to claim 5, characterized in that, The wave height-period data distribution of the wave force time domain signal half-wave includes: the correspondence between the wave height and period of the wave force time domain signal half-wave. The first calculation module is specifically used for: Obtain the wave excitation coefficient of the float, the energy spectrum of the wave, and the pre-set fundamental frequency of the wave; Based on the fundamental frequency, frequency points of multiples of the fundamental frequency are obtained, and the energy spectrum is sampled using the frequency points of multiples of the fundamental frequency to obtain the value of the energy spectrum at multiples of the fundamental frequency; where i∈[0,1,2…n]. The wave excitation coefficient of the float is sampled using the frequency point of the i-th multiple of the fundamental frequency to obtain the wave excitation coefficient of the float at the frequency point of the i-th multiple of the fundamental frequency; wherein, the wave excitation coefficient of the float at the frequency point of the i-th multiple of the fundamental frequency includes: the amplitude and the first phase at the frequency point of the i-th multiple of the fundamental frequency. Between 0 and 2π, a second phase is randomly assigned to frequency points that are multiples of the fundamental frequency (i). The wave force time-domain signal is calculated using the following formula. : ； in, This represents the number of frequency points that are multiples of the fundamental frequency; Indicates the first phase; This represents the numerical value of the energy spectrum at frequency points that are multiples of the fundamental frequency (i). Indicates the fundamental frequency; Indicates the second phase; This represents the time domain signal of wave force; This represents the amplitude at a frequency point that is a multiple of the fundamental frequency (i). The zero-crossing point of the wave force time-domain signal is detected to obtain the half-wave of each wave force time-domain signal in the wave force time-domain signal; The wave height and period of each wave force time domain signal half wave are extracted to obtain the correspondence between wave height and period of each wave force time domain signal half wave.

7. A computer-readable storage medium storing a computer program thereon, characterized in that, When the computer program is run by the processor, it performs the steps of the method described in any one of claims 1-4.

8. An electronic device, characterized in that, The electronic device includes a memory, a processor, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the processor of the steps of the method according to any one of claims 1-4.

Images