METHOD OF DETERMINING AVERAGE WIND SPEED BY MEANS OF A LiDAR SENSOR

Abstract

The present invention relates to a method of determining the average wind speed in a vertical plane by use of a LiDAR sensor (2), comprising performing measurements (MES), constructing a measurement model (MOD M) and a wind model (MOD V), then using an adaptive Kalman filter (KAL) to determine the wind speed (v), and determining the average wind speed in the vertical plane being considered (RAWS).

Claims

1-14. (canceled)

15. A method of determining the average wind speed in a vertical plane with a LiDAR sensor positioned on a wind turbine, comprising steps of: a) constructing a model of the LiDAR measurements; b) constructing a wind model accounting for spatial coherence and temporal coherence of wind speed; c) performing measurements with LiDAR sensor of wind amplitude and direction in at least one vertical wind measurement plane spaced a distance from the wind turbine; d) determining speed at predefined estimation points in space upstream from the wind turbine by use of an adaptive Kalman filter using the model of the LiDAR measurements, the wind model, and measurements; and e) determining the average wind speed in the vertical plane at the distance from the wind turbine by use of the wind speeds determined for the predefined estimation points within the at least one vertical plane considered.

16. A method as claimed in claim 15, wherein the model of the LiDAR measurements is written as: m.sub.j,x(k)=a.sub.jv.sub.j,x(k)+b.sub.jv.sub.j,y(k)+c.sub.jv.sub.j,x(k), with m being the measurement, x being the longitudinal direction, j being a measurement beam of said LiDAR sensor, m.sub.j,x being the measurement on measurement beam j at distance x, k being the discrete time, v being the wind speed, v.sub.j,x the longitudinal component of the wind speed for measurement beam j, v.sub.j,y being the transverse component of the wind speed for measurement beam j, v.sub.j,z being the vertical component of the wind speed for measurement beam j, and a.sub.j, b.sub.j, c.sub.j being constant measurement coefficients for measurement beam j.

17. A method as claimed in claim 15, wherein the spatial coherence of the wind model is a function of a transverse coherence, a vertical coherence and a longitudinal coherence.

18. A method as claimed in claim 16, wherein the spatial coherence of the wind model is a function of a transverse coherence, a vertical coherence and a longitudinal coherence.

19. A method as claimed in claim 17, wherein the transverse coherence is written as: v.sub.x,y.sub.1=f.sub.l(v.sub.x,y.sub.2, y.sub.1−y.sub.2), with x being the longitudinal component, y.sub.1 and y.sub.2 being two transverse positions having the same longitudinal and vertical values, v.sub.x,y1 being the longitudinal component of the wind speed at position y.sub.1, v.sub.x,y2 being the longitudinal component of the wind speed at position y.sub.2, and f.sub.t being a predefined function.

20. A method as claimed in claim 18, wherein the transverse coherence is written as: v.sub.x,y.sub.1=f.sub.t(v.sub.x,y.sub.2, y.sub.1−y.sub.2), with x being the longitudinal component, y.sub.1 and y.sub.2 being two transverse positions having the same longitudinal and vertical values, v.sub.x,y1 being the longitudinal component of the wind speed at position y.sub.1, v.sub.x,y2 being the longitudinal component of the wind speed at position y.sub.2, and f.sub.t being a predefined function.

21. A method as claimed in claim 17, wherein the vertical coherence is written as: $v_{x,z1}={v_{x,z2}\left(\frac{z_1}{z_2}\right)}^{\alpha },$ with x being a longitudinal component, z.sub.1 and z.sub.2 being vertical positions having the same longitudinal and transverse values, v.sub.x,z1 being the longitudinal component of the wind speed at position z.sub.1, v.sub.x,z2 being a longitudinal component of wind speed at position z.sub.2, and α being the coefficient of the power law.

22. A method as claimed in claim 17, wherein the longitudinal coherence is written as: v.sub.x,x.sub.1(k)=f.sub.l(v.sub.x,x.sub.2(k), x.sub.1−x.sub.2), with x being a longitudinal component, k being discrete time, x.sub.1 and x.sub.2 being longitudinal positions having identical transverse and vertical values, v.sub.x,x1 being a longitudinal component of the wind speed at position x.sub.1, v.sub.x,x2 the longitudinal component of the wind speed at position x.sub.2, and f.sub.l being a predefined function.

23. A method as claimed in claim 15, wherein the temporal coherence of wind model is written as: ω(k)=A.sub.sω(k−1), with k being discrete time, ω being a vector comprising first longitudinal components of wind speed at n predefined estimation points, then transverse components of the wind speed for the n predefined estimation points, A.sub.s being a constant matrix which is the autocorrelation function of the wind speed obtained by a Kaimal spectrum.

24. A method as claimed in claim 15, wherein the adaptive Kalman filter (KAL) is applied to equations as follows: v.sub.x(k)=A.sub.sv.sub.x(k−1)+η(k) and $\{\begin{array}{c}{v}_{x,y_1}\left(k\right)-f_t\left(v_{x,y_2}\left(k\right),y_1-y_2\right)-{ϵ}_t\left(x_1,x_2{y}_1{y}_2{z}_1,z_2,k\right)\\ v_{x,z_1}-{v_{x,z_2}\left(\frac{z_1}{z_2}\right)}^{\alpha }={ϵ}_v\left(x_1,x_2{y}_1{y}_2{z}_1,z_2,k\right)\\ v_{x,x_1}\left(k\right)-f_l\left(v_{x,x_2}\left(k\right),x_1-x_2\right)={ϵ}_l\left(x_1,x_2{y}_1{y}_2{z}_1,z_2,k\right)\\ m_{j,x}\left(k\right)=a_j{v}_{j,x}\left(k\right)+b_j{v}_{j,y}\left(k\right)+c_j{v}_{j,z}\left(k\right)+{ϵ}_m\left(k\right)\end{array},$ with k being the discrete time, v being wind speed, x being a longitudinal component, y.sub.1 and y.sub.2 being transverse positions having identical longitudinal and vertical values, x.sub.1 and x.sub.2 being longitudinal positions having identical transverse and vertical values, z.sub.1 and z.sub.2 being vertical positions having identical longitudinal and transverse values, v.sub.x,y1 being a longitudinal component of the wind speed at position y.sub.1, v.sub.x,y2 being a longitudinal component of the wind speed at position y.sub.2, f.sub.t being a predefined function, v.sub.x,x1 being a longitudinal component of the wind speed at position x.sub.1, v.sub.x,x2 being a longitudinal component of the wind speed at position x.sub.2, f.sub.l a predefined function, v.sub.x,z1 the longitudinal component of the wind speed at position z.sub.1, v.sub.x,z2 being a longitudinal component of the wind speed at position z.sub.2, α being a coefficient of a power law, j being a measurement beam of the LiDAR sensor, m.sub.j,x being a measurement on measurement beam j at distance x, v.sub.j,x being a longitudinal component of wind speed for measurement beam j, v.sub.j,y being a transverse component of wind speed for measurement beam j, v.sub.j,z being a vertical component of wind speed for measurement beam j, a.sub.j, b.sub.j, c.sub.j being constant measurement coefficients for measurement beam j, η being noise of the equation of state, ε.sub.t being transverse noise, ε.sub.v being vertical noise, ε.sub.l being longitudinal noise, ε.sub.m being measurement noise, A.sub.s being a constant matrix which is an autocorrelation function of wind speed obtained by a Kaimal spectrum.

25. A method as claimed in claim 15, wherein the wind speed is determined at different points using equations: $\{\begin{array}{c}\hat{w}\left(k|k-1\right)A_s\hat{w}\left(k-1|k-1\right)\\ P\left(k|k-1\right)=A_{s}P\left(k-1|k-1\right)A_s^T+Q\end{array}\mathrm{and}\{\begin{array}{c}K\left(k\right)=P\left(k|k-1\right){C_{\alpha }^T\left(C_{a}P\left(k|k-1\right)C_{\alpha }^T+R\right)}^{-1}\\ \hat{w}\left(k|k\right)=\hat{w}\left(k|k-1\right)+K\left(k\right)\left(y\left(k\right)-C_a\hat{w}\left(k|k-1\right)\right)\\ P\left(k|k\right)=\left(I-K\left(k\right)C_a\right)P\left(k|k-1\right)\end{array},$ with k being discrete time, ω being a vector comprising first longitudinal components of wind speed at n predefined estimation points, {circumflex over (ω)}(k|k−1) being an estimation of vector ω(k) given the measurements performed until time k−1, {circumflex over (ω)}(k|k) being an estimation of vector ω(k) given measurements performed until time k, P(k|k−1) a covariance matrix of vector ω(k) given the measurements performed until time k−1, P(k|k) being a covariance matrix of vector ω(k) given the measurements performed until time k, A.sub.s is a constant matrix which is the autocorrelation function of wind speed obtained by the Kaimal spectrum, Q and R being covariance matrices of noises ε(k) and η(k), and C.sub.a is obtained by linearizing output equations around {circumflex over (ω)}(k|k−1), y(k) being the measurements of the LiDAR sensor and I being an identity matrix.

26. A method as claimed in claim 15, wherein the wind speed is determined in the vertical plane at a distance from the wind turbine by use of the average of the longitudinal components of the wind speed at points belonging to the vertical plane, for wind speeds considered included in a projection of a surface swept by the rotor of the wind turbine in the vertical plane being considered.

27. A method of controlling a wind turbine, comprising: a) determining the average wind speed by use of the method as claimed in claim 15; and b) controlling the wind turbine according to the average wind speed.

28. A computer program product comprising code instructions which performs steps of the method as claimed in claim 15, when the program is executed on at least one of a control and diagnosis unit of the wind turbine.

29. A LiDAR sensor, comprising a processor which performs the method as claimed in claim 15.

30. A wind turbine comprising a LiDAR sensor as claimed in claim 29, with LiDAR sensor being positioned on a nacelle of the wind turbine or in a hub of the wind turbine.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Other features and advantages of the method according to the invention will be clear from reading the description hereafter of embodiments given by way of non-limitative example, with reference to the accompanying drawings wherein:

[0039] FIG. 1 illustrates the steps of the method of determining the average wind speed according to an embodiment of the invention:

[0040] FIG. 2 illustrates a wind turbine equipped with a LiDAR sensor according to an embodiment of the invention:

[0041] FIG. 3 illustrates, for a first example, the comparison between the average wind speed 100 m from the rotor of a wind turbine obtained with the method according to one embodiment of the invention and the reference average wind speed; and

[0042] FIG. 4 illustrates, for a second example, the comparison between the average wind speed 110 m from the rotor of a wind turbine obtained with the method according to one embodiment of the invention and the reference average wind speed.

DETAILED DESCRIPTION OF THE INVENTION

[0043] The present invention relates to a method of determining the average wind speed in a vertical plane, by a LiDAR sensor positioned on a wind turbine.

[0044] According to the invention, the LiDAR sensor measures the wind speed over at least one measurement plane upstream from the wind turbine. There are several types of LiDAR sensors, for example scanning LiDAR sensors, continuous wave or pulsed LiDAR sensors. Within the context of the invention, a pulsed LiDAR is preferably used. However, the other LiDAR technologies may also be used while remaining within the scope of the invention.

[0045] LiDAR sensors provide fast measurement. Therefore, using such a sensor enables fast continuous determination of the average wind speed. For example, the sampling rate of the LiDAR sensor can range between 1 and 5 Hz (or more in the future), and it can be 4 Hz. Furthermore, the LiDAR sensor allows obtaining information relative to the wind upstream from the turbine, which is related to the wind blowing towards the turbine. The LiDAR sensor can therefore be used for predicting the wind speed in the turbine rotor plane.

[0046] FIG. 2 schematically shows, by way of non-limitative example, a horizontal-axis wind turbine 1 equipped with a LiDAR sensor 2 for the method according to one embodiment of the invention. LiDAR sensor 2 is used to measure the wind speed at a given distance over measurement planes PM (only two measurement planes are shown). Knowing the wind measurement in advance a priori allows providing substantial information. This figure also shows axes x, y and z. The reference point of this coordinate system is the center of the rotor. Direction x is the longitudinal direction corresponding to the direction of the rotor axis, upstream from the wind turbine, which also corresponds to the measurement direction of LiDAR sensor 2. Direction y, which is perpendicular to direction x, is the lateral or transverse direction located in a horizontal plane (directions x, y form a horizontal plane). Direction z is the vertical direction (substantially corresponding to the direction of tower 4) pointing up, axis z is perpendicular to axes x and y. The rotor plane is indicated by the rectangle in dotted line PR which is defined by directions y, z for a zero value of x. Measurement planes PM are planes formed by directions y, z at a distance from rotor plane PR (for a non-zero value of x). Measurement planes PM are parallel to rotor plane PR.

[0047] Conventionally, a wind turbine 1 converts the kinetic energy of the wind into electrical or mechanical energy. For conversion of wind energy to electrical energy, the following elements are used:

[0048] a tower 4 allows a rotor (not shown) to be positioned at a sufficient height to enable motion thereof (necessary for horizontal-axis wind turbines) and allowing at least one of the rotor to be positioned at a height enabling it to be driven by stronger and more regular winds than at ground level 6. Tower 4 generally houses part of the electric and electrical components (modulator, control, multiplier, generator, etc.),

[0049] a nacelle 3 mounted at the top of tower 4, housing mechanical, pneumatic and some electrical and electronic components (not shown) necessary for operating the machine. Nacelle 3 can rotate to orient the machine in the right direction,

[0050] a rotor, fastened to the nacelle, comprising blades 7 (generally three) and the hub of the wind turbine. The rotor is driven by the wind energy and it is connected by a mechanical shaft, directly or indirectly (via a gearbox and mechanical shaft system), to an electrical machine (electric generator) (not shown) that converts the energy recovered to electrical energy. The rotor is potentially provided with control systems such as a variable-angle blades or aerodynamic brakes,

[0051] optionally a transmission made up of two shafts (mechanical shaft of the rotor and mechanical shaft of the electrical machine) connected by a transmission (gearbox) (not shown).

[0052] As is visible in FIG. 2, which is an example embodiment of a pulsed LiDAR sensor, the LiDAR sensor 2 comprises four measurement beams or axes (b1, b2, b3, b4). By way of non-limitative example, the method according to the invention also works with a LiDAR sensor comprising any number of beams. The LiDAR sensor performs a punctual measurement at each point of intersection of a measurement plane PM and a beam (b1, b2, b3, b4). These measurement points are represented by black circles in FIG. 2, for the first measurement plane PM with the measurement points being denoted by PT1, PT2, PT3 and PT4. Processing the measurements at these measurement points allows determination of the wind speed in measurement planes PM.

[0053] Preferably, LiDAR sensor 2 can be mounted on nacelle 3 of wind turbine 1 or in the hub of wind turbine 1 (that is at the front end of the nacelle in the wind direction).

[0054] According to the invention, the method of determining the average wind speed comprises the steps of:

[0055] 1) Construction of a LiDAR sensor measurement model

[0056] 2) Construction of a wind model

[0057] 3) Wind measurement

[0058] 4) Determination of the wind speed

[0059] 5) Determination of the average wind speed.

[0060] Steps 3), 4) and 5) are carried out in real time. Steps 1) and 2) can be carried out offline and prior to the real-time steps, and they can be performed in this order, in the reverse order or simultaneously. All the steps are described in detail in the rest of the description.

[0061] FIG. 1 schematically illustrates, by way of non-limitative example, the steps of the method according to an embodiment of the invention. The method determines the average wind speed in a vertical plane by use of a LiDAR sensor arranged on a wind turbine. A first step constructs offline a wind model MOD V and a measurement model MOD M. Then, the amplitude and the direction of the wind MES are measured in real time by means of the LiDAR sensor. The wind speed v is then determined in real time at various points by use of an adaptive Kalman filter KAL, which uses wind model MOD V, measurement model MOD M and measurements MES. Finally, the average wind speed RAWS is determined from wind speed v at different points.

[0062] 1) Construction of a LiDAR Sensor Measurement Model

[0063] This step constructs a model of the LiDAR sensor measurements. It is a model relating the components of the wind speed to the measurement signal from the LiDAR sensor.

[0064] According to one embodiment of the invention, the LiDAR sensor measurement model can be written as follows: m.sub.j,x(k)=a.sub.jv.sub.j,x(k)+b.sub.jv.sub.j,y(k)+c.sub.jv.sub.j,z(k), with m being the measurement, x being the longitudinal direction, j being a measurement beam of the LiDAR sensor, m.sub.j,x being the measurement on measurement beam j at distance x, k being the discrete time, v the wind speed, v.sub.j,x the longitudinal component of the wind speed for measurement beam j, v.sub.j,y being the transverse component of the wind speed for measurement beam j, v.sub.j,z being the vertical component of the wind speed for measurement beam j, a.sub.j, b.sub.j, c.sub.j being constant measurement coefficients for measurement beam j. Measurements coefficients a.sub.j, b.sub.j, c.sub.j depend only on the beam angles of the LiDAR sensor and are not dependent on the measurement distances. The measurement coefficients a.sub.j, b.sub.j, c.sub.j can be data provided by the LiDAR sensor manufacturer.

[0065] 2) Construction of a Wind Model

[0066] This step constructs a wind model, which accounts for the spatial coherence and the temporal coherence which define the wind speed and its components at any point in space according to various parameters, notably time and the position in space (therefore according to the coordinates of the point considered in the (x, y, z) system). In other words, a wind model meeting the spatial coherence constraints and the temporal coherence constraints is constructed. These spatial and temporal coherences allow the wind model to be representative of the wind and to provide precise determination of the wind speed.

[0067] According to one implementation of the invention, the wind model can determine the longitudinal and transverse components of the wind speed. Alternatively, the wind model can determine the three components of the wind speed.

[0068] According to one embodiment of the invention, the spatial coherence used in the wind model can depend on a transverse coherence, a longitudinal coherence and a vertical coherence. The representativity of the wind model is thus improved.

[0069] For this embodiment, the transverse coherence can be written by the equation: v.sub.x,y.sub.1=f.sub.t(v.sub.x,y.sub.2, y.sub.1−y.sub.2), with x being the longitudinal component, y.sub.1 being and y.sub.2 being two transverse positions having the same longitudinal (x.sub.1=x.sub.2=x) and vertical (z.sub.1=z.sub.2=z) values, v.sub.x,y1 being the longitudinal component of the wind speed at position y.sub.1, v.sub.x,y2 being the longitudinal component of the wind speed at position y.sub.2, f.sub.t being a known predefined function. Thus, the longitudinal component of the wind speed at point y.sub.1 depends on the longitudinal component of the wind speed at point y.sub.2 and on the distance between points y.sub.1 and y.sub.2. According to an example embodiment, predefined function f.sub.t can be an exponential function.

[0070] For this embodiment, the vertical coherence can be written by the equation as:

[00004]$v_{x,z1}={v_{x,z2}\left(\frac{z_1}{z_2}\right)}^{\alpha }$

with x being the longitudinal component, z.sub.1 and z.sub.2 being vertical positions having the same longitudinal (x.sub.1=x.sub.2=x) and transverse (y.sub.1=y.sub.2=y) values, v.sub.x,z1 being the longitudinal component of the wind speed at position z.sub.1, v.sub.x,z2 the longitudinal component of the wind speed at position z.sub.2, α the coefficient of the power law. For this equation, the reference framework of the height z is defined with respect to the base of the wind turbine tower (and not at the LiDAR sensor). Thus, the longitudinal component of the wind speed at point z.sub.1 depends on the longitudinal component of the wind speed at point z.sub.2 and on the ratio between the heights of points z.sub.1 and z.sub.2. Coefficient α of the power law can be chosen to be constant, or it can be estimated using LiDAR sensor measurements, for example according to the method described in the patent application FR-19/06,569.

[0071] For this embodiment, the longitudinal coherence can be written by the equation: v.sub.x,x.sub.1(k)=f.sub.l(v.sub.x,x.sub.2(k), x.sub.1−x.sub.2), with x being the longitudinal component, x.sub.1 and x.sub.2 being longitudinal positions having the same transverse (y.sub.1=y.sub.2=y) and vertical (z.sub.1=z.sub.2=z) values, v.sub.x,x1 being the longitudinal component of the wind speed at position x.sub.1, v.sub.x,x2 being the longitudinal component of the wind speed at position x.sub.2, f.sub.l being a known predefined function. Thus, the longitudinal component of the wind speed at point x.sub.1 depends on the longitudinal component of the wind speed at point x.sub.2 and on the distance between points x.sub.1 and x.sub.2. According to an example embodiment, predefined function f.sub.l can be an exponential function.

[0072] The temporal coherence is understood to be the variation with time of the wind speed components in a single position, that is for the same values x, y and z. In other words, the temporal coherence can be formulated as a relation between the wind speed components between two consecutive discrete time intervals, denoted by k and k−1.

[0073] According to an implementation of the invention, one known temporal coherence is obtained using the Kaimal spectrum that can be defined by:

[00005]$S_t\left(f\right)={\sigma }_t^2\frac{4\frac{L_t}{U}}{{\left(1+6f\frac{L_t}{U}\right)}^{\frac{5}{3}}},$

with f being the frequency in Hertz, t being the component of the wind speed (t can therefore correspond to x, y or z), S.sub.t being the Kaimal spectrum of component t of the wind speed, U being the average wind speed at the height of the wind turbine rotor, L.sub.t being the integral scale parameter of component t of the wind speed and σ.sub.t being the variance determined by the wind turbulence intensity. Indeed, the Kaimal spectrum allows determination of a discrete transfer function that can relate a wind value at time k to a wind value at time k−1.

[0074] For the embodiment where only the longitudinal and transverse components of the wind speed are determined, ω can be a vector of dimensions 2n, which first comprises the longitudinal components of the wind speed for the n points being considered, then the transverse components of the wind speed for the n points are considered. To illustrate this vector ω in a simple case, when a first point has longitudinal and transverse wind speed components v.sub.x1, v.sub.y1, and a second point has longitudinal and transverse wind speed components v.sub.x2, v.sub.y2, vector ω is written as:

ω=(v.sub.x1v.sub.x2v.sub.y1v.sub.y2).sup.T.

[0075] Using this notation and noting that the Kaimal spectrum is the Fourier transform of the autocorrelation function of the wind speed, the following equation can be written for the temporal coherence: ω(k)=A.sub.gφ(k−1), with A.sub.s being a constant matrix which is the autocorrelation function of the wind speed obtained by a Kaimal spectrum. Matrix A.sub.s can be obtained from the Kaimal spectrum formula as defined above. Thus, this equation gives the connection between wind speed w at time k and wind speed ω at time k−1.

[0076] Alternatively, for the temporal coherence, the von Karman spectrum or any similar representation can be used. 3) Wind measurement

[0077] In this step, the wind amplitude and direction are continuously measured in at least one measurement plane distant from the wind turbine, by the LiDAR sensor. This measurement corresponds to the signal received by the LiDAR sensor in response to the signal emitted by the LiDAR sensor. Indeed, by interferometry and Doppler effect, part of the laser signal emitted by the LiDAR sensor is reflected by the air molecules at the measurement point and also by the aerosols (suspended dust and microparticles).

[0078] According to an implementation of the invention, the measurement planes can be at a longitudinal distance (along axis x of FIG. 2) from the rotor plane preferably ranging between 50 and 400 m. It is thus possible to determine the evolution of the wind speed over a long distance upstream from the wind turbine, which also allows the accuracy of the average wind speed determination to be improved.

[0079] According to an embodiment of the invention, the wind speed measurement can be performed in several measurement planes (whose measurement distances are not imposed by the method according to the invention) to facilitate wind speed determination, which allows the user of the LiDAR sensor to freely parametrize the LiDAR sensor.

[0080] For the embodiment using a pulsed LiDAR, the measurements are obtained successively at the measurement points illustrated in FIG. 2, starting with beam b1, then beam b2, . . . and finally beam b4. An interesting characteristic of this system is that it allows measurement of the projection of the wind speed at several distances, simultaneously, for a given beam. It is thus possible to obtain, for example, 10 successive distances between 50 m and 400 m, at a sampling rate of the LiDAR sensor. At each sampling time, only the measurements of the selected current beam are refreshed.

[0081] 4) Determination of the Wind Speed

[0082] This step determines the wind speed at various points of the space upstream from the wind turbine, by use of an adaptive Kalman filter using the wind model constructed in step 2, the LiDAR sensor measurement model constructed in step 1 and the measurements performed in step 3. The various wind speed determination points are predefined estimation points. Application of the Kalman filter allows obtaining a state observer. The adaptive Kalman filter enables adaptation of the noise covariance matrix according to the wind speed. Thus, the filter is efficient over a wide wind speed range. Moreover, the adaptive Kalman filter is robust wind speed variations.

[0083] It should be observed that a state observer or a state estimator is, in automation and systems theory, an extension of a model represented as a state representation. When the state of the system is not measurable, an observer allowing the state to be reconstructed from a model is constructed.

[0084] For an embodiment using the equations illustrated in step 2, the following state model can be written, with the equation of state:

v.sub.x(k)=A.sub.sv.sub.x(k−1)+η(k)

and the output equations:

[00006]$\{\begin{array}{c}{v}_{x,y_1}\left(k\right)-f_t\left(v_{x,y_2}\left(k\right),y_1-y_2\right)-{ϵ}_t\left(x_1,x_2{y}_1{y}_2{z}_1,z_2,k\right)\\ v_{x,z_1}-{v_{x,z_2}\left(\frac{z_1}{z_2}\right)}^{\alpha }={ϵ}_v\left(x_1,x_2{y}_1{y}_2{z}_1,z_2,k\right)\\ v_{x,x_1}\left(k\right)-f_l\left(v_{x,x_2}\left(k\right),x_1-x_2\right)={ϵ}_l\left(x_1,x_2{y}_1{y}_2{z}_1,z_2,k\right)\\ m_{j,x}\left(k\right)=a_j{v}_{j,x}\left(k\right)+b_j{v}_{j,y}\left(k\right)+c_j{v}_{j,z}\left(k\right)+{ϵ}_m\left(k\right)\end{array},$

with η being the noise of the equation of state, ε.sub.t being the transverse noise, ε.sub.v being the vertical noise, ε.sub.l being the longitudinal noise and ε.sub.m being the measurement noise.

[0085] Thus, the problem of estimating vector ω(k) becomes a state estimation problem, which does not require imposing the position of the measurement planes of the LiDAR sensor. One way of estimating the unknown state vector ω(k), which can take into account the information on noises η(k) and ε(k), applies the algorithm of the adaptive Kalman filter, with the following notation:

[00007]$\mathrm{.Math.}=\left(\begin{array}{c}{\mathrm{.Math.}}_t\\ {\mathrm{.Math.}}_v\\ {\mathrm{.Math.}}_l\\ {\mathrm{.Math.}}_m\end{array}\right).$

Indeed, the adaptive Kalman filter provides the solution to the optimization problem:

[00008]$\underset{w\left(k\right)}{\min }{J}_r\left(k\right)$

with

[00009]$J_r\left(k\right)=\left\{{\left(w\left(0\right)-\stackrel{\_}{w}\left(0\right)\right)}^T{P}_0^{-1}\left(w\left(0\right)-\stackrel{\_}{w}\left(0\right)\right)+\underset{j=1}{\overset{k}{.Math.}}\left({\eta \left(j-1\right)}^T{Q}^{-1}\eta \left(j-1\right)+{ϵ\left(j\right)}^T{R}^{-1}ϵ\left(j\right)\right)\right\}$

where P.sub.0, Q and R are adjustment matrices of suitable dimensions, and ω(0) is the average value of the initial state ω(0).

[0086] In order to solve this optimization problem by use of the adaptive Kalman filter, the following hypotheses can be made, notably for a mathematical interpretation of P.sub.0, Q and R: [0087] ω(0) is a random vector uncorrelated with noises η(k) and ε(k) [0088] ω(0) has a known average ω(0) with P.sub.0 as the covariance matrix, that is:

P.sub.0=E[ω(0)−ω(0))(ω(0)−ω(0)).sup.T] [0089] η(k) and ε(k) are zero-mean uncorrelated white noise processes with covariance matrices Q and R respectively, i.e.:

[00010]$E\left[\eta \left(k\right){\eta \left(j\right)}^T\right]=\{\begin{array}{c}Q\mathrm{if}k=j\\ 0\mathrm{if}k\ne j\end{array}E\left[ϵ\left(k\right){ϵ\left(j\right)}^T\right]=\{\begin{array}{c}R\mathrm{if}k=j\\ 0\mathrm{if}k\ne j\end{array}E\left[ϵ\left(k\right){\eta \left(j\right)}^T\right]=0\mathrm{for}\mathrm{all}k,j$

[0090] This last hypothesis implies that Q and R are symmetric positive semidefinite matrices.

[0091] Furthermore, given that, in the state model, noises ε.sub.l, ε.sub.v, and ε.sub.t depend on measurement distances x.sub.1, x.sub.2, y.sub.1, y.sub.2, z.sub.1, z.sub.2, covariance matrix R is adapted according to the measurement distances. According to one embodiment, R can be a polynomial function of the measurement distances. Alternatively, R can be obtained from a map, a neural network, etc.

[0092] The following notations can be adopted: [0093] {circumflex over (ω)}(k|k−1) is the estimation of vector ω(k) given the measurements performed until time k−1; [0094] {circumflex over (ω)}(k|k) is the estimation of vector ω(k) given the measurements performed until time k; [0095] P(k|k−1) is the covariance matrix of vector ω(k) given the measurements performed until time k−1; and [0096] P(k|k) is the covariance matrix of vector ω(k) given the measurements performed until time k.

[0097] Then, the algorithm of the adaptive Kalman filter is used to determine the wind speed at various points, using the following equations:

[0098] On the one hand, a temporal update:

[00011]$\{\begin{array}{c}\hat{w}\left(k|k-1\right)A_s\hat{w}\left(k-1|k-1\right)\\ P\left(k|k-1\right)=A_{s}P\left(k-1|k-1\right)A_s^T+Q\end{array}$

[0099] On the other hand, a measurement update:

[00012]$\{\begin{array}{c}K\left(k\right)=P\left(k|k-1\right){C_{\alpha }^T\left(C_{a}P\left(k|k-1\right)C_{\alpha }^T+R\right)}^{-1}\\ \hat{w}\left(k|k\right)=\hat{w}\left(k|k-1\right)+K\left(k\right)\left(y\left(k\right)-C_a\hat{w}\left(k|k-1\right)\right)\\ P\left(k|k\right)=\left(I-K\left(k\right)C_a\right)P\left(k|k-1\right)\end{array}$

with C.sub.a obtained by linearizing the output equations of the state model around {circumflex over (ω)}(k|k−1), y(k) the measurements of the LiDAR sensor and I the identity matrix.

[0100] Thus, these steps allow determination of the vector ω, which comprises the components of the wind speed at several different points. In other words, these steps allow determination of the components of the wind speed at several different points.

[0101] 5) Determination of the Average Wind Speed

[0102] This step determines the average wind speed in a vertical plane at a distance upstream from the wind turbine (the distance is defined by means of the longitudinal direction) by use of the wind speeds determined in step 4, in particular the wind speeds in the vertical plane are considered.

[0103] According to one embodiment, the average wind speed can be the average of the longitudinal components of the wind speed in the plane being considered.

[0104] According to a preferred embodiment of the invention, the average wind speed can be the average of the longitudinal components of the wind speed in the plane considered, which considers only the values of the wind speed in a surface corresponding to the surface swept by the rotor of the wind turbine. In other words, the surface swept by the rotor of the wind turbine (a circle of radius corresponding to the length of the wind turbine blades at nacelle height) is projected onto the vertical plane being considered, and the wind speeds are averaged for the points of the vertical plane belonging to this projection. This average speed is generally referred to as RAWS (Rotor Average Wind Speed) and commonly used for at least one of control, diagnosis, and monitoring of a wind turbine.

[0105] The present invention also relates to a method of controlling a wind turbine equipped with a LiDAR sensor. The following steps are carried out for this method:

[0106] determining the average wind speed by use of the method of determining the average wind speed according to any one of the variants described above; and

[0107] controlling the wind turbine according to the average wind speed that is determined.

[0108] Precise real-time determination of the average wind speed allows suitable wind turbine control in terms of minimization of the effects on the turbine structure and maximization of the recovered power. Indeed, through this control, the LiDAR sensor allows reduction of the loads on the structure, the blades and the tower representing 54% of the cost. Using a LiDAR sensor therefore allows optimizing the wind turbine structure and thus reducing the costs and maintenance.

[0109] The method can further comprise an intermediate step that determines the average wind speed in the rotor plane of the wind turbine from the average wind speed determined by the method. The wind movement time between the vertical plane and the rotor plane can therefore be taken into account (it can be calculated notably by considering Taylor's frozen turbulence hypothesis), it is also possible to account for the induction phenomenon between the vertical plane and the rotor plane (by use of an induction factor for example). The induction factor reflects the wind deceleration upstream from the wind turbine related to the presence of the wind turbine blades. The wind turbine is then controlled according to the average wind speed in the rotor plane.

[0110] According to an implementation of the invention, at least one of the inclination angle of the blades and the electrical recovery torque of the wind turbine generator can be controlled according to the wind speed. Other types of regulation devices can be used.

[0111] According to an embodiment of the invention, at least one of the inclination angle of the blades and electrical recovery torque can be determined by use of wind turbine maps according to the wind speed at the rotor. For example, the control method described in patent application FR-2,976,630 A1 which corresponds to US patent application 2012-0,321,463 can be applied.

[0112] Furthermore, the invention relates to a computer program product comprising code instructions designed to carry out the steps of one of the methods described above (method of determining the wind speed in the rotor plane, control method). The program can be executed on a processor of the LiDAR sensor or any similar processor linked to the LiDAR sensor or to the wind turbine.

[0113] According to an aspect, the present invention also relates to a LiDAR sensor for a wind turbine, comprising a processor configured to implement one of the methods described above (method of determining the average wind speed, control method).

[0114] According to an implementation of the invention, the LiDAR sensor can be a scanning LiDAR sensor, a continuous wave LiDAR sensor or a pulsed LiDAR sensor. The LiDAR sensor is preferably a pulsed LiDAR sensor.

[0115] The invention also relates to a wind turbine, notably an offshore (at sea) or an onshore (on land) wind turbine equipped with a LiDAR sensor as described above. According to an embodiment of the invention, the LiDAR sensor can be arranged on the nacelle of the wind turbine or in the hub of the turbine (at the end of the nacelle of the wind turbine). The LiDAR sensor is so oriented to measure the wind upstream from the turbine (i.e. before the wind turbine and along the longitudinal axis thereof, designated by axis x in FIG. 2). According to an embodiment, the wind turbine can be identical to the wind turbine illustrated in FIG. 2.

[0116] For the embodiment of the control method, the wind turbine can comprise a control, for example for control of the pitch angle of at least one blade of the wind turbine or of the electrical torque, for implementing the control method according to the invention.

[0117] It is clear that the invention is not limited to the embodiments of the methods described above by way of example and that it encompasses any variant embodiment.

Example

[0118] The features and advantages of the method according to the invention will be clear from reading the example hereafter.

[0119] For this example, the wind is simulated by a simulator, as well as the LiDAR sensor measurements, and the average wind speed is determined by the method according to an embodiment of the invention. This embodiment of the invention uses the spatial and temporal coherence equations described, and it determines the average longitudinal component of the wind speed in a vertical plane.

[0120] According to a first configuration, the measurement plane distances are: [50, 70, 90, 100, 120, 140, 160, 180, 190, 200] meters.

[0121] FIG. 3 illustrates a comparison of the average wind speed RAWS in m/s as a function of time T for a 100 m distance between the rotor plane and the vertical plane. In this figure, the curve in dotted line represents reference curve REF and the curve in solid line represents the average wind speed curve EST obtained by the method according to the invention. It is noted that the two curves are nearly superposed. Therefore, the method according to the invention enables precise determination of the average wind speed.

[0122] According to a second configuration, the measurement plane distances are: [50, 80, 90, 110, 130, 150, 170, 180, 190, 200] meters.

[0123] FIG. 4 illustrates a comparison of the average wind speed RAWS in m/s as a function of time T for a 110 m distance between the rotor plane and the vertical plane. In this figure, the curve in dotted line represents reference curve REF and the curve in solid line represents the average wind speed curve EST obtained by the method according to the invention. It is noted that the two curves are nearly superposed. Therefore, the method according to the invention enables precise determination of the average wind speed.

[0124] These two curves also allow showing that the method is accurate whatever the distance considered, without imposing any measurement plane distance.