Method for determining the available power of a wind park, and corresponding wind park

11578696 · 2023-02-14

Assignee

Inventors

Cpc classification

International classification

Abstract

A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having rotor blades, the blade angle of which can be adjusted is provided. A wind farm which is set up to carry out the method for determining an available power is provided. The method comprises providing a shading matrix which determines at least one effective wind speed of each of the wind power installations in the wind farm as a function of at least one wind speed and wind direction and wind farm throttling using a park wake model. The method makes it possible to accurately determine an available power of a wind farm even when the wind farm is operated with throttled power.

Claims

1. A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having a plurality of rotor blades, wherein blade angles of the plurality of rotor blades are adjustable, the method comprising: providing or using a shading matrix to determine an effective wind speed of each of the plurality of wind power installations in the wind farm as a function of at least one: wind speed, wind direction, and wind farm throttling using a farm wake model; and determining the available power of the wind farm based on the effective wind speed of each of the plurality of wind power installations, wherein the method for determining the available power of the wind farm comprises: determining a first set of the plurality of wind power installations having freely inflowing wind for a particular wind direction; ascertaining an average wind speed and wind direction as an average value of the wind speeds and wind directions of the first set of the plurality of wind power installations with freely inflowing wind; ascertaining an effective wind speed of each wind power installation for each of the plurality of wind power installations without freely inflowing wind from the average wind speed and wind direction taking into account the shading matrix for an unthrottled wind farm; ascertaining an available power of each of the plurality of wind power installations based on their power curve and the ascertained effective wind speed; and ascertaining the available power of the wind farm by summing the available powers of each of the plurality of wind power installations.

2. The method as claimed in claim 1, further comprising determining the shading matrix from a maximum installation power.

3. The method as claimed in claim 1, wherein the farm wake model is adaptive and comprises a self-learning method for adapting to changes in environment and for correcting errors in the farm wake model.

4. The method as claimed in claim 3, wherein different shading matrices are provided based on different farm wake models, wherein one of the shading matrices is selected based on a parameter indicative of atmospheric stability or a parameter interpolated between different shading matrices.

5. The method as claimed in claim 1, further comprising low-pass filtering for correcting errors in the available power.

6. A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having a plurality of rotor blades, wherein blade angles of the plurality of rotor blades are adjustable, the method comprising: providing or using a shading matrix to determine an effective wind speed of each of the plurality of wind power installations in the wind farm as a function of at least one: wind speed, wind direction, and wind farm throttling using a farm wake model; and determining the available power of the wind farm based on the effective wind speed of each of the plurality of wind power installations, wherein the method for determining the available power of the wind farm comprises: determining a first set of the plurality of wind power installations having freely inflowing wind for a particular wind direction; determining a wind speed and wind direction for the first set of the plurality of wind power installations with freely inflowing wind; ascertaining effective wind speeds of a second set of the plurality of wind power installations without freely inflowing wind taking into account the shading matrix for an unthrottled wind farm for each of the determined wind speeds and wind directions of the first set of the plurality of wind power installations with freely inflowing wind; ascertaining an installation-specific average effective wind speed, for at least the second set of the plurality of wind power installations without freely inflowing wind, as an average value of the effective wind speeds ascertained for each wind power installation with freely inflowing wind; ascertaining an available power of each of the plurality of wind power installations based on their power curve and the average effective wind speed ascertained in an installation-specific manner; and ascertaining the available power of the wind farm by summing the available powers of each wind power installation.

7. The method as claimed in claim 6, further comprising providing the available power of the wind farm in real time.

8. The method as claimed in claim 6, wherein determining the available power is also based on losses of a network of the wind farm.

9. A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having a plurality of rotor blades, wherein blade angles of the plurality of rotor blades are adjustable, the method comprising: providing or using a shading matrix to determine an effective wind speed of each of the plurality of wind power installations in the wind farm as a function of at least one: wind speed, wind direction, and wind farm throttling using a farm wake model; and determining the available power of the wind farm based on the effective wind speed of each of the plurality of wind power installations, wherein the method for determining the available power of the wind farm comprises: ascertaining individual wind speeds and wind directions of the plurality of wind power installations; calculating modeled wind speeds and wind directions based on a set farm-effective wind speed and wind direction taking into account the shading matrix and the wind farm throttling; optimizing the set farm-effective wind speed and wind direction based on a comparison of individual wind speeds and wind directions and the modeled wind speeds and wind directions; ascertaining the effective wind speed of each wind power installation from the optimized set farm-effective wind speed and wind direction taking into account the shading matrix for an unthrottled wind farm; ascertaining an available power of each of the plurality of wind power installations based on their power curve and the ascertained effective wind speed; and ascertaining the available power of the wind farm by summing the available powers of each of the plurality of wind power installations.

10. The method as claimed in claim 9, further comprising storing the shading matrix for each wind power installation as a table based on the farm-effective wind speed and direction and a maximum farm power.

11. The method as claimed in claim 9, wherein determining the available power of the wind farm is also based on air density.

12. The method as claimed in claim 9, wherein an effective wind direction of each wind power installation is ascertained taking into account the shading matrix.

13. A wind farm comprising a plurality of wind power installations, wherein the wind farm includes a controller configured to carry out the method as claimed in claim 9.

14. A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having a plurality of rotor blades, wherein blade angles of the plurality of rotor blades are adjustable, the method comprising: providing or using a shading matrix to determine an effective wind speed of each of the plurality of wind power installations in the wind farm as a function of at least one: wind speed, wind direction, and wind farm throttling using a farm wake model; and determining the available power of the wind farm based on the effective wind speed of each of the plurality of wind power installations, wherein the method for determining the available power of the wind farm comprises: determining a wind speed and wind direction for each of the plurality of wind power installations; inverting the shading matrix for each of the plurality of wind power installations to estimate a steady farm-effective wind speed and wind direction based on the individual wind speed and wind direction of the respective wind power installation and the wind farm throttling; ascertaining the effective wind speed of each wind power installation from the respectively estimated steady farm-effective wind speed and wind direction taking into account the shading matrix for an unthrottled wind farm; ascertaining an available power of each of the plurality of wind power installations based on their power curve and the respective ascertained effective wind speed; and ascertaining the available power of the wind farm by summing the available powers of each of the plurality of wind power installations.

15. The method as claimed in claim 14, further comprising: averaging all estimated farm-effective wind speeds and wind directions; wherein ascertaining the effective wind speed of each wind power installation comprises ascertaining the effective wind speed of each wind power installation from the averaged farm-effective wind speed and wind direction taking into account the shading matrix for the unthrottled wind farm.

16. A wind power installation in a wind farm, wherein the wind power installation includes a controller configured to carry out the method as claimed in claim 14.

17. A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having a plurality of rotor blades, wherein blade angles of the plurality of rotor blades are adjustable, the method comprising: providing or using a shading matrix to determine an effective wind speed of each of the plurality of wind power installations in the wind farm as a function of at least one: wind speed, wind direction, and wind farm throttling using a farm wake model; determining the available power of the wind farm based on the effective wind speed of each of the plurality of wind power installations; and correcting the wind speed and wind direction based on spreading of wind field through the farm, wherein the correcting is carried out as a weighted sliding average value.

18. The method as claimed in claim 17, further comprising: determining an average wind speed and wind direction acting on a rotor of at least one wind power installation based on an operating point of the at least one wind power installation.

19. The method as claimed in claim 18, wherein determining the average wind speed comprises: interpolation from at least three-dimensional tables of the wind speed as a function of a rotational speed, an aerodynamic power, and a pitch angle; and iteration based on a family of characteristics as a function of a tip-speed ratio and a pitch angle.

20. The method as claimed in claim 17, wherein the sliding average value is adapted based on a flow time.

21. The method as claimed in claim 17, wherein determining the available power is also based on further limitations associated with one or more of technical, maintenance, reduced-sound operation and shadowing of the plurality of wind power installations in the wind farm.

22. A method for determining an available power of a wind farm, wherein the wind farm comprises a plurality of wind power installations with a rotor having a plurality of rotor blades, wherein blade angles of the plurality of rotor blades are adjustable, the method comprising: providing or using a shading matrix to determine an effective wind speed of each of the plurality of wind power installations in the wind farm as a function of at least one: wind speed, wind direction, and wind farm throttling using a farm wake model; and determining the available power of the wind farm based on the effective wind speed of each of the plurality of wind power installations, wherein the farm wake model is adaptive and comprises a self-learning method for adapting to changes in environment and for correcting errors in the farm wake model, wherein different shading matrices are provided based on different farm wake models, wherein one of the shading matrices is selected based on a parameter indicative of atmospheric stability or a parameter interpolated between different shading matrices, and wherein each of the different farm wake models are based on different wake decay constants.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

(1) Further configurations and advantages are described below with reference to the accompanying figures.

(2) FIG. 1 schematically shows a wind power installation,

(3) FIG. 2 schematically shows a wind farm,

(4) FIG. 3 schematically shows a first method for determining an available power of a wind farm,

(5) FIG. 4 schematically shows a second method for determining an available power of a wind farm,

(6) FIG. 5 schematically shows a third method for determining an available power of a wind farm,

(7) FIG. 6 schematically shows a fourth method for determining an available power of a wind farm,

(8) FIG. 7 schematically shows a fifth method for determining an available power of a wind farm,

(9) FIG. 8 schematically shows a sixth method for determining an available power of a wind farm, and

(10) FIG. 9 shows a block diagram for explaining the algorithm on which the determination of the equivalent wind speed is based.

DETAILED DESCRIPTION

(11) FIG. 1 shows a schematic illustration of a wind power installation 100. The wind power installation 100 has a tower 102 and a nacelle 104 on the tower 102. An aerodynamic rotor 106 having three rotor blades 108 and a spinner 110 is provided on the nacelle 104. The aerodynamic rotor 106 is caused to rotate by the wind during operation of the wind power installation and therefore also rotates an electrodynamic rotor of a generator which is directly coupled to the aerodynamic rotor 106. The electrical generator is arranged in the nacelle 104 and generates electrical energy.

(12) FIG. 2 shows a wind farm 112 having, by way of example, three wind power installations 100 which may be identical or different. The three wind power installations 100 are therefore representative of fundamentally any desired number of wind power installations in a wind farm 112. The wind power installations 100 provide their power, namely the generated current in particular, via an electrical farm network 114. In this case, the currents or powers respectively generated by the individual wind power installations 100 are added, and a transformer 116 is usually provided and steps up the voltage in the farm in order to then feed it into the supply network 120 at the feed-in point 118 which is also generally referred to as the PCC. FIG. 2 is only a simplified illustration of a wind farm 112 which does not show a controller, for example, even though a controller is naturally present. The farm network 114 may also be configured differently, for example, by virtue of a transformer also being present at the output of each wind power installation 100, for example, to name just one other exemplary embodiment.

(13) The wind farm 112 can be operated with reduced power, for example in the event of a network-related power reduction. In this case, one or more or all of the wind power installations 100 operate(s) with a power which is below the power possible for the wind situation. In order to provide a P.sub.available signal, wind power installations must also be able to ascertain their available power during such reduced-power operation. As mentioned, wind power installations 100 in the wake see a higher wind speed during the power reduction than in the non-reduced case, which results in an overestimation of the available farm power. The aim of the method for determining an available power of a wind farm 112 is therefore to obtain the available power of a wind farm, as the sum of the powers of the wind power installations, taking into account the farm effect, on the basis of, for example, a shearing force, a terrain model and a farm plan for the given throttled farm power, wind speed and wind direction.

(14) This consideration is effected, as described with reference to the subsequent figures, using a shading matrix 200 which determines at least one effective wind speed of each of the wind power installations 100 in the wind farm 112 as a function of at least one wind speed and wind direction and wind farm throttling using a farm wake model. Various methods for implementing the shading matrix 200 are described in detail below with reference to FIGS. 3 to 8.

(15) FIG. 3 schematically shows a first method a) for determining the available power of the wind farm 112 using the shading matrix 200. The borders with rounded corners in the figures schematically combine a hierarchical or structural responsibility of the wind farm, individual installation controller and farm controller or farm computer relating to individual method steps or calculations.

(16) The wind farm 112 shown in FIG. 3 and also in the subsequent FIGS. 4 to 8 comprises nine wind power installations 100 which are arranged in a strict matrix arrangement. It goes without saying that this is only a very simple example of an actual wind farm 112 and the teaching can likewise be applied to any desired wind farms. The wind farm 112 symbolically has inflowing wind 130 from the left, with the result that three wind power installations 132 have freely inflowing wind. The further wind power installations 134 are in the wake of the wind power installations 132 with freely inflowing wind and experience a reduced wind speed.

(17) In a first step 310, the installations 132 with freely inflowing wind are selected or determined from the wind power installations 100. The determination in step 310 is carried out on the basis of the wind direction.

(18) In step 320, an undisrupted wind speed and wind direction of the wind power installations 132 with freely inflowing wind are then determined. The undisrupted wind speed and wind direction correspond to the wind speed which is not influenced by the rotor. It is preferably ascertained using the wind estimator which is described in detail with reference to FIG. 9. Alternatively, the undisrupted wind speed can also be determined by means of anemometers and the like fitted to the wind power installation with suitable corrections.

(19) The ascertained undisrupted wind speeds and wind directions custom character are averaged in a step 330 in order to obtain an average wind speed and wind direction custom character In FIG. 3, the wind speed and wind direction are represented as a wind vector, for example. The average value is represented as the sum of the wind vectors divided by the number n.sub.frei of wind power installations 132 with freely inflowing wind. It should be noted that this formation of an average value is only one example and other functions can also be used in a similar manner to form the average value. Together with the average wind speed and wind direction custom character a power during the non-limited operation of the wind farm P.sub.nom, is used in a step 340 to obtain effective wind speeds v.sub.1 . . . n of each wind power installation taking into account the shading matrix 200. Non-limited operation is present if the wind farm is operated with the permitted maximum power and not with a lower, throttled power P.sub.gedr, for example on account of network requirements.

(20) As an alternative to all wind power installations, the wind power installations 132 with freely inflowing wind can be completely or partially excluded from step 340 and the available power of the wind power installations 132 with freely inflowing wind can be ascertained directly on the basis of the power curve for the wind speed determined in step 320, that is to say the wind speed determined in step 320 corresponds to the effective wind speed for the wind power installations with freely inflowing wind in this embodiment.

(21) In a step 350, the effective wind speed v.sub.1 . . . n of each wind power installation is corrected on account of spreading of the wind field through the farm. A flow time through the farm is determined, for example, as the quotient of the extent of the wind farm 112 in the inflow direction to the wind speed, for example the estimated wind speed. Other possibilities for correction in step 350 are also possible. The correction in step 350 corresponds to a farm-layout-dependent flow delay through the wind farm.

(22) The corrected effective wind speeds are then transformed in step 360 into installation-dependent power values on the basis of the power curves P.sub.LK, are summed and are transferred to a farm-available power P.sub.verfg,Park.

(23) FIG. 4 schematically shows a second method b) for determining an available power of a wind farm 112. Like in the first method a), the wind power installations 132 with freely inflowing wind are first of all determined according to the second method b) in step 310 and their wind speed and wind direction are determined in step 320.

(24) An effective wind speed v.sub.1 . . . n,1 . . . n.sub.frei is determined for each of the determined wind speeds and wind directions of the wind power installations with freely inflowing wind in step 430 for each wind power installation 100 using a shading matrix 200 for an unthrottled wind farm. Therefore, for each of the n wind power installations 100 in the wind farm 112 of the number of wind power installations with freely inflowing wind n.sub.frei, effective wind speeds v.sub.i,1 . . . n.sub.frei are accordingly determined in step 430.

(25) For each of the wind power installations, an average value of the effective wind speeds for this wind power installation is then determined in step 440 and the available power of each of the wind power installations is ascertained on the basis of their power curve for the average effective wind speed ascertained in an installation-specific manner.

(26) In an alternative embodiment, the wind power installations 132 with freely inflowing wind can be completely or partially excluded from steps 430 and 440 and the available power of the wind power installations 132 with freely inflowing wind can be directly determined on the basis of the power curve for the wind speed determined in step 320.

(27) The ascertained available powers P.sub.1 . . . n are corrected in a step 450, which corresponds to a farm-layout-dependent flow delay through the wind farm, as in step 350. The difference is only that, in the second method b), the available powers P.sub.1 . . . n are corrected instead of the ascertained effective wind speeds v.sub.1 . . . n.

(28) Finally, the available powers are summed in step 460 in order to obtain the available power of the wind farm 112.

(29) FIG. 5 schematically shows a third method c) for determining an available power of the wind farm 112. In contrast to the first method a) and the second method b), all wind power installations 100 in the wind farm 112 are selected in the third method c) in step 510 and their wind speed and wind direction are determined in step 520. The method for ascertaining the undisrupted wind speeds and wind directions v.sub.Schätz,1 . . . n and step 520 can correspond to the method in step 320, that is to say can be carried out, for example, with the aid of the wind estimator which is described with reference to FIG. 9.

(30) The third method c) is an iteration method which is aimed at optimally correlating the wind speed and direction for all wind power installations. For this purpose, an effective farm wind speed custom character is set and model speeds of each wind power installation custom character are determined in step 530 for a throttled power P.sub.gedr on the basis of the shading matrix 200. The modeled wind speeds and wind directions custom character resulting from the farm speed custom character which has been set are correlated with the undisrupted individual wind speeds and wind directions custom character and a measure of the correlation of both variables is determined in step 540. In the example in FIG. 5, a quadratic distance is determined from both, which distance is intended to be as short as possible in order to achieve an optimum correlation. For this purpose, the farm wind speed custom character which has been set is adapted until the correlation achieves a particular quality. All iterative optimization algorithms can be used as an alternative to the least squares method shown.

(31) The optimized farm-effective wind speed custom character is then determined in step 550 taking into account the shading matrix 200 and, for an unthrottled wind farm with power P.sub.nom, the effective wind speed v.sub.i of each wind power installation is determined.

(32) The available power of the farm P.sub.verfg,Park can then be determined and summed in step 560 using the power curves and the ascertained effective wind speeds of each of the wind power installations.

(33) FIG. 6 schematically shows a fourth method d) for ascertaining the available powers P.sub.verfg,Park of the wind farm 112. Steps 510 and 520 correspond to the steps shown in FIG. 5.

(34) The individual wind speeds and wind directions custom character are supplied, in step 630, to an inverted shading matrix 200′ in order to estimate a steady farm-effective wind speed and wind direction v.sub.Park,1 . . . n together with the throttled power P.sub.gedr. The shading matrix 200 is inverted in an installation-specific manner in order to arrive at the shading matrix 200′.

(35) The steady farm-effective wind speed custom character estimated for each wind power installation is averaged in step 640 and the averaged farm-effective wind speed custom character is used in step 650 to determine effective wind speeds v.sub.1 . . . n of each wind power installation using the shading matrix 200 for an unthrottled wind farm with a maximum power P.sub.nom. In step 660, the effective wind speeds v.sub.1 . . . n are then used to determine the available power of the wind farm P.sub.verfg,Park using the respective power curves and subsequent summation.

(36) FIG. 7 schematically shows a fifth method e) for ascertaining the available power P.sub.verfg,Park of the wind farm 112. In contrast to the first to fourth methods, a large part of the method e) can be implemented on a controller of the wind power installation 100 since no operating information relating to the other wind power installations is required. More specifically, all steps combined under 710 can be implemented on the individual wind power installation. Each wind power installation ascertains an available power P.sub.verfg,i on the basis of ascertained individual wind speeds and wind directions v.sub.Schätz,i taking into account the shading acting on the individual installation, which is summarized as the aerodynamic farm efficiency in bracket 750. In step 720, the installation-specific estimation includes an estimation using the inverted shading matrix 200′ with throttled power P.sub.gedr in order to estimate a farm-effective wind speed custom character This is then used in step 730 to infer the installation-effective wind speed v.sub.i using the shading matrix 200 and the technically possible power or the nominal power P.sub.nom. The installation-effective wind speed v.sub.i is then used in step 740 to estimate the available power P.sub.verfg,i using the power curve. The farm controller 760 must then include only the signals of the available power P.sub.verfg,i which are transmitted by the individual wind power installations 100, in step 770. The shading matrix 200 and the inverted shading matrix 200′ can be provided in an installation-specific manner in this method e) and can be implemented on controllers of the respective wind power installations 100.

(37) FIG. 8 shows a further method f) which is a variation of the method e) and, instead of the installation-specific speeds v.sub.Schätz,i, takes into account the shading on the basis of the power values. In the method f) too, the wind farm controller 760 must only sum the individual signals of available power P.sub.verfg,i provided by the installations in step 770. The remaining steps can be implemented on the individual wind power installation in a similar manner to the method e). In step 810, an available power is calculated on the basis of the working point and optionally on the basis of sensor signals from the wind power installation. The available power P.sub.verfg which is individually calculated for all installations is used, together with the determined wind direction and the pitch angle, in step 820 to take into account the shading effects and to provide an available power P.sub.verfg,i corrected by shading effects.

(38) The method f) combines individual steps of the method e) in one step 820; in particular, the inverse shading matrix 200′, which is interpolated on the basis of the throttled power, and the shading matrix 200 of the unthrottled maximum technically possible supply are combined in such a manner that the intermediate step using the wind speed and direction valid throughout the farm in method e) can be dispensed with on the installation or the farm computer and is required only for the (one-off) calculation of the matrices 200, 200′.

(39) For the shading matrix, a combination to form an aerodynamic efficiency, which only implicitly depends on the wind speed, is therefore carried out in step 820. In addition to the pitch angle α, the minimum pitch angle α.sub.min, for which the wind power installation would be operated in the non-power-reduced case, that is to say with the maximum technically possible supply, is provided.

(40) As already stated at the outset, a wind speed determined by a nacelle anemometer is regularly disrupted by rotor effects. Furthermore, the measured values from a nacelle anemometer have a good correlation with the actual power values from a wind power installation only in the case of averaging over long periods of 10 minutes and more and are therefore not particularly well suited to determining the currently available power, as can also be gathered from reference [1]. In addition, the anemometer measurement is greatly dependent on a power reduction and further influences which cannot be used to correct signals or can be used only with difficulty.

(41) A second possibility, as likewise stated at the outset, is therefore to calculate the available power or the individual undisrupted wind speed from the current wind power installation power (for example the power of a wind power installation with full-scale converter which is fed into the DC intermediate circuit) and the blade angle of the wind power installation, which is described below with reference to FIG. 9.

(42) If the power of a wind power installation is limited below the available power, the supply of mechanical power must be reduced. This can be carried out by increasing the blade angle. In this case, the reduction in the aerodynamic power coefficient c.sub.P(λ,α) with an increase in the blade angle α generally depends on the tip-speed ratio λ.

(43) If a wind power installation is now continuously operated at the same tip-speed ratio, in particular at the optimum tip-speed ratio λ.sub.opt, even in the case of varying wind speeds as a result of known power or rotational speed control, a fixed power reduction results from the increase in the blade angle with respect to a particular optimum blade angle α.sub.opt. If the power-reducing effect is determined for each blade angle, the available power can be calculated from the current wind power installation power P.sub.ist and the power reduction coefficient k(α)=1−c.sub.p(λ.sub.opt,α)/c.sub.p(λ.sub.opt,α.sub.opt) according to the following formula:

(44) P verf u .Math. gbar = P ist 1 - k ( α )

(45) However, a constant tip-speed ratio cannot be maintained under all operating conditions. Deviations from the optimum may occur, in particular in the case of severe curtailment. On account of the inversely proportional relationship between the power reduction coefficient and the estimated available power, small deviations in the assumed power reduction coefficient may already result here in errors when calculating the available power, in particular in the case of severe curtailment.

(46) In order to solve this, an algorithm is now proposed which can estimate the available power on the basis of the wind speed and therefore at all operating points of a wind power installation as reliably as possible.

(47) FIG. 9 therefore shows a block diagram 1 of a proposed algorithm and therefore of a proposed method for determining an equivalent wind speed v.sub.Wind which can be used as an undisrupted wind speed v.sub.Schätz according to the methods described above.

(48) Two observers which observe or estimate other values which cannot be measured on the basis of particular measurement variables are used in said algorithm.

(49) The power observer 2 illustrated in FIG. 9 determines the available internal power, as the power P.sub.ava,DC available at the DC voltage intermediate circuit, from the power P.sub.DC fed into the DC voltage intermediate circuit and measured and the rotor rotational speed ω.

(50) The variables of rotational speed (ω) and power in the DC intermediate circuit (P.sub.DC) are therefore captured by means of measurement. The rotational speed (ω) of the generator corresponds to the rotational speed of the rotor if there is no transmission. The available internal power P.sub.ava,DC is determined or calculated therefrom with the aid of the state observer 2, which available internal power can also be referred to as the available aerodynamic power based on the intermediate circuit ({circumflex over (P)}.sub.ava,DC), that is to say it corresponds to the aerodynamic power minus the generator losses and possibly minus further losses such as friction or current conversion losses.

(51) The power observer may be based on the following simplified system description which can also be referred to as an acceleration equation of a rotating single-mass system:

(52) J ω ˙ = T m e c h - T e l = 1 ω ( P m e c h - P e l )

(53) In this formula, J denotes the moment of inertia of the rotating masses of the wind power installation, ω denotes the rotational speed of the wind power installation, T.sub.mech and P.sub.mech denote the torque and power, respectively, arising at the shaft due to the wind and T.sub.el/P.sub.el denote the electromechanical torque and the power, respectively, of the generator. Here a dot once again denotes derivatives of variables with respect to time. If a generator power that exceeds the mechanical power is taken out over a period of time, the wind power installation consequently decelerates.

(54) For the power observer 2 and proceeding from the acceleration equation of a rotating single-mass system, above, it is possible to set up the following state space model of the wind power installation, having the state variables of rotational speed (ω) and mechanical torque (T.sub.mech):

(55) [ ω ˙ T ˙ m e c h ] = [ 0 1 J 0 0 ] [ ω T m e c h ] + [ - 1 J 0 ] T el

(56) Of the variables used here, it is not possible to measure the mechanical torque T.sub.mech and therefore the latter should be calculated from the measurement data by way of a state observer. Since, moreover, the rotational speed signal often is only measured with a low resolution and with a low sampling rate, a state observation is proposed for this value, too. A suitable observer structure can be formulated as follows, with reference being made to reference [4] for further explanation:

(57) [ ω ˆ . T ˆ . m e c h ] = [ - k ω 1 J - k T 0 ] [ ω ˆ T ˆ m e c h ] + [ k ω - 1 J k T 0 ] [ ω T e l ]

(58) Here, in contrast to the measured input variables of rotational speed w and electrical torque

(59) T e l = P D C ω ,
observed variables are denoted by {circumflex over ( )}. P.sub.DC is the electrical power fed into the DC voltage intermediate circuit. The two parameters k.sub.ω, and k.sub.T influence the dynamic behavior and, in the case of a time-discrete implementation, the stability of the state observer, too, and have to be chosen taking these aspects into account.

(60) The aerodynamic power {circumflex over (P)}.sub.ava,DC is calculated by way of the product of rotational speed and torque. Then, this corresponds to the determined available internal power P.sub.ava,DC.

(61) Input variables for determining the available internal power P.sub.ava,DC are the measured electrical power P.sub.DC that is fed into the DC intermediate circuit and the measured wind power installation rotational speed ω. Then, the available electrical power, in relation to the DC intermediate circuit, P.sub.ava,DC is the output value. This observer is important since the wind power installation, under certain circumstances, is not at a steady operating point at which the power P.sub.DC corresponds to the currently available wind power, said wind power installation instead almost permanently accelerating or decelerating during normal operation and, for example, possibly deviating from a normal rotational speed-power characteristic for network support, in which energy is taken from kinetic energy of the rotation.

(62) The generator losses P.sub.vGen are added to the available internal power P.sub.ava,DC that is related to the DC intermediate circuit. These losses can either be determined by measurements during operation or can be gathered from a characteristic for the respective generator working point (optionally taking account of temperature measured values). Additionally, if relevant, friction and current conversion losses can also be added thereto. The result of this addition is the available generator input power of the wind power installation, referred to in a simplifying manner as P.sub.mech. For this calculation, the generator losses expected for the maximum available input power should be used instead of the current generator losses.

(63) This power and the measured wind power installation rotational speed Co and the current blade angle α are used as input values for the second observer, the wind observer 4. The latter observes the wind speed v.sub.Wind that leads to the available power at the current blade angle α. In principle, this relates to a solution v to the equation

(64) P m e c h = 1 2 ρ A v 3 C p ( λ , α ) ,
where P.sub.mech is assumed to be known from the power observer 2, ρ denotes the air density, A denotes the rotor area and C.sub.p (λ,α) denotes the power coefficient depending on the tip-speed ratio λ and the blade angle α. This is a nonlinear equation, in which the wind speed v is included directly and indirectly via the tip-speed ratio

(65) λ = ω R v
(ω in turn denotes the wind power installation rotational speed in rad/s and R denotes the radius of the aerodynamic rotor). Different solutions thereto are known from reference [2].

(66) On account of the nonlinear properties of the equation, a nonlinear approach to determining the wind speed as the basis for the wind observer was found to be advantageous and reference is made in this respect to references [2, 3], listed under the phrase immersion and invariance.

(67) However, in contrast to references [2] and [3], it is not the measured electrical wind power installation power that is used as an input variable for the wind observer 4 in the solution proposed here, but instead the available internal power that is determined or calculated by way of the power observer 2. This is important for the correct, at least desired functionality of the wind observer 4 in the transient operation of the wind power installation (i.e., in the case of strong acceleration or deceleration, for example for network support, in which energy is taken from kinetic energy of the rotation). Preferably, the wind observer 4 also additionally outputs the available power Pava (α.sub.opt,λ.sub.opt) of the wind power installation. This is the power that could be produced by the wind power installation at the equivalent wind speed v.sub.Wind, determined by the wind observer 4, if the optimum blade angle α.sub.opt and the optimum tip-speed ratio λ.sub.opt were set. By way of example, this can be calculated on the basis of the equivalent wind speed v.sub.Wind and taking account of known relationships between wind and produced power in the specified conditions (optimum blade angle α.sub.opt and optimum tip-speed ratio λ.sub.opt). By way of example, an appropriate characteristic can be stored to this end. Additionally, it is proposed to output the current tip-speed ratio λ which can be calculated from the wind speed, rotational speed and rotor diameter. The stated available power P.sub.verf. (α.sub.opt,λ.sub.opt) would be the available power not corrected by farm effects. Whereas the methods a) to e), for example, are based directly on the determined equivalent wind speed v.sub.Wind as the input variable for determining the farm-available power, the method f) can be based on the available power P.sub.verf(α.sub.opt,λ.sub.opt).

(68) The proposed solution can be used for applications of positive and negative balancing energy with wind energy, both primary balancing energy and secondary balancing energy. Here, the improved accuracy, in particular, can be important. The solution can be used as a replacement for a second nacelle anemometer.

(69) By correctly determining the undisrupted wind speed, at least determining the latter with high accuracy, it is possible to calculate the maximum available power of the wind power installation in many operating states. This allows the correct, at least relatively accurate, determination of, e.g., negative balancing power provided, where a wind power installation or a wind farm has to reduce its power by a certain value below the available power. Furthermore, the requirement for a second wind meter, for example for satisfying requirements directed at safety, may be solved by means of software instead of hardware.

(70) The proposed solution is at least an improvement over previously known solutions. The core problem of determining the wind speed lies in solving the aerodynamic power equation:

(71) P m e c h = 1 2 ρ A v 3 C p ( λ , α ) , ( 1 )
where ρ denotes the air density, A denotes the rotor area, which can synonymously also be referred to as the rotor blade area, v denotes the wind speed and C.sub.p(λ, α) denotes the power coefficient as a function of the tip-speed ratio λ and the blade angle α. Even though, strictly speaking, this equation denotes the mechanical power of the rotor at the generator rotor, it can be determined approximately from the measured electrical power and the estimated or calculated generator losses and can therefore be assumed to be known, albeit as a variable afflicted by noise. The difficulty in solving this equation for the unknown v lies in the twofold input of the variables into the equation, both directly (v.sup.3) and via the tip-speed ratio

(72) 0 λ = ω m R v
(ω.sub.m once again denotes the wind power installation rotational speed in rad/s).

(73) A numerically efficient and, at the same time, accurate solution according to reference [2] can be obtained by way of an observer according to the “immersion and invariance” principle, as shown in reference [3]. To this end, the following system of state equations is implemented:

(74) v ˆ . w I = γ [ 1 J P m e c h ω m - Φ ( ω m , v ˆ w I + γ ω m ) ] , ( 2 ) v ˆ w = v ˆ w I + γ ω m , ( 3 )
where

(75) Φ ( ω m , v w ) = 1 2 ρ A J v w 3 ω m C p ( r ω m v w , α ) ,
J denotes the moment of inertia of the wind power installation rotor, P.sub.mech denotes the available DC power of the generator, estimated by the power observer, plus the generator losses, γ>0 is an adjustable adaptation parameter of the filter, ρ denotes the air pressure, A denotes the rotor area and r denotes the rotor radius of the wind power installation. The derivative of the variable custom character.sub.w.sup.I with respect to time is denoted by a dot.

(76) The output values of the wind observer are calculated by numerically solving this system of state equations. Here, the variable custom character.sub.w calculated in equation (3) forms the estimated wind speed v.sub.Wind. The estimated tip-speed ratio λ then emerges from the measured rotational speed and the estimated wind speed according to the aforementioned equation. The available power emerges from the solution of equation (1) while taking account of losses arising in the wind power installation at the estimated optimum working point and of possible technical or operational restrictions, which may result in a limitation of this available power.

(77) A decisive modification of the wind observer described here in relation to the solution described in reference [3] lies in the use of an available power, estimated by a separate observer, for the calculation of custom character.sub.w.sup.I (equation (2)). By contrast, the measured electrical power or the electrical torque is used directly in reference [3], under the assumption that the wind power installation in the partial load range is always operated near its optimum operating point. However, in practice, there is a significant deviation between the current electrical power and the power achievable at the optimum operating point, or the power available in the wind, on account of acceleration and deceleration processes in the case of a variable wind speed and, in part, also as desired on account of certain requirements such as network frequency support, in which energy is taken from kinetic energy of the rotation, and so a separate determination of the available power is proposed prior to solving the aerodynamic power equation by way of the wind observer for the purposes of accurately determining the wind speed at all operating points.

(78) A further deviation in relation to the solution presented in reference [3] consists of taking account of the current blade angle α when determining the power coefficient C.sub.P(λ,α). This is of decisive importance for using the wind observer to determine the maximum possible power since the blade angle, as a rule, is greater than the optimum blade angle in the curtailed operation of a wind power installation. Consequently, for determining the maximum available power, the wind speed estimated by the wind observer and the calculation of the power coefficient for this wind speed and the optimum blade angle α.sub.opt and an optimum tip-speed ratio λ.sub.opt, i.e., C.sub.p(λ.sub.opt,α.sub.opt) and the solution of the aerodynamic power equation (1) with the estimated wind speed and the optimum C.sub.p(λ.sub.opt,α.sub.opt) are advantageous.

(79) An alternative but complicated model is known from reference [5], for example.

(80) Consequently, an improvement in the accuracy of determining the available power in the case of strong curtailment is created.

(81) Particular advantages of the solution shown with reference to FIG. 9 may consist in determining the available power at, where possible, all working points at which the wind power installation is in operation, taking account of operating states in which the wind power installation deviates from its steady working point, inter alia when providing network services, determining the equivalent wind speed in the rotor plane without additional sensor systems and, moreover determining the rotational speed and acceleration of the wind power installation in a highly accurate manner in the case of an appropriate implementation.

(82) In joint consideration with the methods described in FIG. 3 to FIG. 8, particularly good, simple and accurate determination of the available farm power can be achieved in combination.

REFERENCES

(83) [1] K. E. Johnson, L. Y. Pao, M. J. Balas, L. J. Fingersh: Control of Variable-Speed Wind Turbines—Standard and Adaptive Techniques for Maximizing Energy Capture. IEEE Control Systems Magazine, June 2006, pages 70-81. [2] M. N. Soltani, T. Knudsen, M. Svenstrup, R. Wisniewski, P. Brath, R. Ortega, K. Johnson: Estimation of Rotor Effective Wind Speed: A Comparison. IEEE Transactions on Control Systems Technology, 21(4), July 2013, pages 1155-1167. [3] R. Ortega, F. Macilla-David, F. Jaramillo: A Globally Convergent Wind Speed Estimator for Windmill Systems. In Proc. 2011 50.sup.th IEEE Conference on Decision and Control, Orlando, Fla., USA, December 2011, pages 6079-6084. [4] C. M. Verrelli, A. Savoia, M. Mengoni, R. Marino, P. Tomei, L. Zarri: On-line Identification of Winding Resistances and Load Torque in Induction Machines. IEEE Transactions on Control Systems Technology, vol. 22(4), July 2014. [5] T. Göçmen, G. Giebel, P.-E. Réthoré, J. P. Murcia Leon, J. R. Kristoffersen: Uncertainty Quantification of the Real-Time Reserves for Offshore Wind Power Plants, WindEurope Summit, September 2016