METHOD FOR PROCESSING TELEMETRY DATA FOR ESTIMATING A WIND SPEED

Abstract

A method for processing telemetry data for estimating a wind speed. The method includes a hybridization by temporal combination, and/or by weighting, and/or by averaged projection.

Claims

1. A method for processing telemetry data for estimating a wind speed, said method comprising hybridization by temporal combination comprising: a step (A) of vector reconstruction of at least two components of an average wind speed vector over a time interval (Ω), called partition time interval, starting from successive projections, over time, of an instantaneous wind speed vector; a step (B) of scalar reconstruction of at least one average wind speed value (Vh.sub.ave) over a time interval (T), called reference time interval, starting from a number (T/Ω) of the at least two components of the average wind speed vector reconstructed in step A; 2Ω is less than or equal to T and T/Ω corresponds to the number of the at least two components of the average wind speed over the partition time interval Ω included in the reference time interval T.

2. The method according to claim 1 comprising: in step A, reconstruction, on the basis of equations (1) to (7), of the at least two components (U.sub.Ω, V.sub.Ω) or (V.sub.Ω, W.sub.Ω) or (U.sub.Ω, W.sub.Ω) among three components (U.sub.Ω, V.sub.Ω, W.sub.Ω) of the average wind speed vector over the partition time interval Ω; the component U.sub.Ω being the component of the average wind speed vector in a spatial direction (d1) extending in a spatial plane (p1) and the component V.sub.Ω being the component of the average wind speed vector in a spatial direction (d2) extending in the spatial plane p1 and the component W.sub.Ω being the component of the average wind speed vector in a spatial direction (d3) orthogonal to the plane p1: $\begin{array}{cc}\left[\mathrm{Math}1\right]& \\ U_{\Omega }=\frac{\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Ni}}\right)-\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Si}}\right)}{2.\sin \theta },\mathrm{or}& \mathrm{equation}1\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}2\right]& \\ U_{\Omega }=\frac{\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Ni}}\right)-\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Si}}\right)}{2.\sin \theta .\cos \left(\frac{\alpha }{2}\right)},\mathrm{and}/\mathrm{or}& \mathrm{equation}6\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}3\right]& \\ V_{\Omega }=\frac{\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Ei}}\right)-\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Wi}}\right)}{2.\sin \gamma },\mathrm{or}& \mathrm{equation}2\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}4\right]& \\ V_{\Omega }=\frac{\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Ei}}\right)-\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Wi}}\right)}{2.\sin \gamma .\cos \left(\frac{\alpha }{2}\right)},\mathrm{and}/\mathrm{or}& \mathrm{equation}7\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}5\right]& \\ W_{\Omega }=\frac{\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Ni}}\right)+\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Si}}\right)}{2.\cos \theta },\mathrm{or}& \mathrm{equation}3\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}6\right]& \\ W_{\Omega }=\frac{\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Ei}}\right)+\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Wi}}\right)}{2\cos \gamma },\mathrm{or}& \mathrm{equation}4\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}7\right]& \\ W_{\Omega }=\frac{1}{M}.{.Math.}_M\left(S_{\mathrm{Vi}}\right),& \mathrm{equation}5\end{array}$ in which i is an integer comprised between 1 and M corresponding to the successive projections S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi, and S.sub.vi , over time, of the instantaneous wind speed vector over the partition time interval Ω; S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi, and S.sub.vi are the respective projections of the instantaneous wind speed vector along, respectively, a first axis (a1), a second axis (a2), a third axis (a3), a fourth axis (a4) and a fifth axis (a5) merged with the direction d3, 0 is a non-zero angle formed between the axis a1 and a normal to the plane p1 and between the axis a2 and the normal to the plane p1 and γ is a non-zero angle formed between the axis a3 and the normal to the plane p1 and the axis a4 and the normal to the plane p1, the first and second axes a1 and a2 are included in a plane (p2), the third and fourth axes a3 and a4 are included in a plane (p3) and the planes p2 and p3 form a non-zero angle αbetween them; in step B, scalar reconstruction, on the basis of equations (8) to (10) and of the at least two components of the average wind speed vector reconstructed in step A, of at least one average wind speed value (Vh.sub.ave) over the reference time interval T in plane p1 or p2 or p3, respectively: $\begin{array}{cc}\left[\mathrm{Math}8\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.1}}=\frac{1}{Q}.{.Math.}_Q\left(\sqrt{{\left(U_{\Omega }\right)}^2+{\left(V_{\Omega }\right)}^2+2.U_{\Omega }.V_{\Omega }.\cos \alpha }\right),& \mathrm{equation}8\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}9\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.2}}=\frac{1}{Q}.{.Math.}_Q\left(\sqrt{{\left(U_{\Omega }\right)}^2+{\left(W_{\Omega }\right)}^2}\right),& \mathrm{equation}9\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}10\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.3}}=\frac{1}{Q}.{.Math.}_Q\left(\sqrt{{\left(V_{\Omega }\right)}^2+{\left(W_{\Omega }\right)}^2}\right),& \mathrm{equation}10\end{array}$ Q is an integer comprised between 1 and (T/Ω) corresponding to the number of the at least two components of the average wind speed over the partition time interval Ω included in the reference time interval T.

3. The method according to claim 1, in which a value of the partition time interval Ω is constant or is modified during acquisition of the telemetry data, said value of the partition time interval Ω being a function of: the type of telemetry system from which the telemetry data are acquired, and/or the atmospheric conditions during acquisition of said telemetry data.

4. A method for processing telemetry data for estimating a wind speed, said method comprising hybridization by weighting comprising: a step (C) of vector reconstruction of at least two components of an instantaneous wind speed vector starting from projections of the instantaneous wind speed vector; a step (D) of vector reconstruction over a time interval (T), called reference time interval, of at least two components of an average wind speed vector starting from a number N of the at least two components, comprised over the reference time interval T, of the instantaneous wind speed vector reconstructed in step C; a step (E) of scalar reconstruction of at least one instantaneous wind speed value starting from the at least two components of the average wind speed vector reconstructed in step C; a step (F) of determining at least one average wind speed value starting from the at least one instantaneous wind speed value reconstructed in step E; a step (G) of determining at least one average wind speed value over the reference time interval T starting from the at least two components of the average wind speed vector reconstructed in step D; and a step (H) of determining at least one average wind speed value (Vh.sub.ave) over the time interval T by weighting of a sum of the at least one average wind speed value reconstructed in step F and of the at least one average wind speed value determined in step G.

5. The method according to claim 4 comprising: in step C, vector reconstruction, on the basis of the respective equations (11) to (17), of the at least two components (U.sub.i, V) or (V.sub.i, W.sub.i) or (U.sub.i, W.sub.i) among three components (U.sub.i, V.sub.i, W.sub.i) of the instantaneous wind speed vector; i is an integer comprised between 1 and N corresponding to the number of successive projections of the instantaneous wind speed vector over the reference time interval T, U, being the component of the instantaneous wind speed vector in a spatial direction (d1) extending in a spatial plane (p1) and the component V, being the component of the instantaneous wind speed vector in a spatial direction (d2) extending in the spatial plane p1 and the component W, being the component of the average wind speed vector in a spatial direction (d3) orthogonal to the plane p1: $\begin{array}{cc}\left[\mathrm{Math}11\right]& \\ U_i=\frac{S_{\mathrm{Ni}}-S_{\mathrm{Si}}}{2.\sin \theta },\mathrm{or}& \mathrm{equation}11\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}12\right]& \\ U_i=\frac{S_{\mathrm{Ni}}-S_{\mathrm{Si}}}{2.\sin \theta .\cos \left(\frac{\alpha }{2}\right)},\mathrm{and}/\mathrm{or}& \mathrm{equation}16\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}13\right]& \\ V_i=\frac{S_{\mathrm{Ei}}-S_{\mathrm{Wi}}}{2.\sin \gamma },\mathrm{or}& \mathrm{equation}12\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}14\right]& \\ V_i=\frac{S_{\mathrm{Ei}}-S_{\mathrm{Wi}}}{2.\sin \gamma .\cos \left(\frac{\alpha }{2}\right)},\mathrm{and}/\mathrm{or}& \mathrm{equation}17\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}15\right]& \\ W_i=\frac{S_{\mathrm{Ei}}+S_{\mathrm{Wi}}}{2.\cos \gamma },\mathrm{or}& \mathrm{equation}13\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}16\right]& \\ W_i=\frac{S_{\mathrm{Ni}}+S_{\mathrm{Si}}}{2.\cos \theta },\mathrm{or}& \mathrm{equation}14\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}17\right]& \\ W_i=S_{\mathrm{Vi}},& \mathrm{equation}15\end{array}$ in which S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi, and S.sub.vi are projections of the instantaneous wind speed vector along, respectively, a first axis (a1), a second axis (a2), a third axis (a3), a fourth axis (a4) and a fifth axis (a5) merged with the direction d3, 6 is a non-zero angle formed between the axis a1 and a normal to the plane p1 and between the axis a2 and the normal to the plane p1 and γ is a non-zero angle formed between the axis a3 and the normal to the plane p1 and the axis a4 and the normal to the plane p1, the first and second axes a1 and a2 are included in a plane (p2), the third and fourth axes a3 and a4 are included in a plane (p3) and the planes p2 and p3 form a non-zero angle α between them, in step D, vector reconstruction, on the basis of equations (18) to (20), of the at least two components (Uvect.sub.N, Vvect.sub.N) or (Vvect.sub.N, Wvect.sub.N) or (Uvect.sub.N, Wvect.sub.N) of the average wind speed vector over the reference time interval T; the component Uvect.sub.N being the component of the wind speed in the spatial direction d1, the component Vvect.sub.N being the component of the wind speed in the spatial direction d2 and the component Wvect.sub.N being the component of the wind speed in the spatial direction d3: $\begin{array}{cc}\left[\mathrm{Math}18\right]& \\ {\mathrm{Uvect}}_N=\frac{1}{N}.{.Math.}_N\left(U_i\right),\mathrm{and}/\mathrm{or}& \mathrm{equation}18\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}19\right]& \\ {\mathrm{Vvect}}_N=\frac{1}{N}.{.Math.}_N\left(V_i\right),\mathrm{and}/\mathrm{or}& \mathrm{equation}19\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}20\right]& \\ {\mathrm{Wvect}}_N=\frac{1}{N}.{.Math.}_N\left(W_i\right),& \mathrm{equation}20\end{array}$ in step E, scalar reconstruction, on the basis of equations (21) to (23), of the at least one value (Vscal.sub.i) of the instantaneous wind speed; Vscal.sub.i corresponding to a temporal series of the instantaneous wind speed value in plane p1 or p2 or p3, respectively: [Math 21]
Vscal.sub.i.1=√{square root over ((U.sub.i).sup.2+(V.sub.i).sup.2+2.U.sub.i.V.sub.i. cos α)},  equation 21 [Math 22]
Vscal.sub.i.2=√{square root over ((U.sub.i).sup.2+(W.sub.i).sup.2)},  equation 22 [Math 23]
Vscal.sub.i.3=√{square root over ((V.sub.i).sup.2+(W.sub.i).sup.2)},  equation 23, in step F, determination, on the basis of equations (24) to (26) and starting from the value Vscal.sub.i.1 or Vscal.sub.i.2 or Vscal.sub.i.3 of the instantaneous wind speed reconstructed in step E, of the at least one average wind speed value (Vhscal.sub.ave) in plane p1, p2 or p3 respectively over the reference time interval T: $\begin{array}{cc}\left[\mathrm{Math}24\right]& \\ {\mathrm{Vhscal}}_{\mathrm{ave}\text{.1}}=\frac{1}{N}.{.Math.}_N\left({\mathrm{Vscal}}_{i\text{.1}}\right),& \mathrm{equation}24\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}25\right]& \\ {\mathrm{Vhscal}}_{\mathrm{ave}\text{.2}}=\frac{1}{N}.{.Math.}_N\left({\mathrm{Vscal}}_{i\text{.2}}\right),& \mathrm{equation}25\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}26\right]& \\ {\mathrm{Vhscal}}_{\mathrm{ave}\text{.3}}=\frac{1}{N}.{.Math.}_N\left({\mathrm{Vscal}}_{i\text{.3}}\right),& \mathrm{equation}26\end{array}$ in step G, determination, on the basis of equations (27) to (29) and starting from the at least two reconstructed components of the average wind speed vector, of the at least one average wind speed value (Vvect.sub.ave) in plane p1, p2 or p3 respectively over the reference time interval T: $\begin{array}{cc}\left[\mathrm{Math}27\right]& \\ {\mathrm{Vhvect}}_{\mathrm{ave}\text{.1}}=\sqrt{{\left({\mathrm{Uvect}}_N\right)}^2+{\left({\mathrm{Vvect}}_N\right)}^2+2.{\mathrm{Uvect}}_N.{\mathrm{Vvect}}_N.\cos \alpha },& \mathrm{equation}27\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}28\right]& \\ {\mathrm{Vhvect}}_{\mathrm{ave}\text{.2}}=\sqrt{{\left({\mathrm{Uvect}}_N\right)}^2+{\left({\mathrm{Wvect}}_N\right)}^2},& \mathrm{equation}28\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}29\right]& \\ {\mathrm{Vhvect}}_{\mathrm{ave}\text{.3}}=\sqrt{{\left({\mathrm{Vvect}}_N\right)}^2+{\left({\mathrm{Wvect}}_N\right)}^2},& \mathrm{equation}29\end{array}$ in step H, calculation, on the basis of equations (30) to (32) and starting from the pairs of reconstructed wind speed values (Vhscal.sub.ave.1, Vhvect.sub.ave.1) or (Vhscal.sub.ave.2 and Vhvect.sub.ave.2) or (Vhscal.sub.ave.3, Vhvect.sub.ave.3), of at least one weighted average wind speed value (Vh.sub.ave) in plane p1, p2 or p3 respectively over the reference time interval T: [Math 30]
Vh.sub.ave.1=(1−P).Vhscal.sub.ave.1+P.Vhvect.sub.ave.1,  equation 30 [Math 31]
Vh.sub.ave.2=(1−P).Vhscal.sub.ave.2+P.Vhvect.sub.ave.2,  equation 31 [Math 32]
Vh.sub.ave.3=(1−P).Vhscal.sub.ave.3+P.Vhvect.sub.ave.3,  equation 32 in which P is a dimensionless weighting factor comprised between 0 and 1.

6. The method according to claim 4, in which the factor P is greater than 0.2 and/or less than 0.6.

7. The method according to claim 4, in which the value of the factor P is constant or is modified during acquisition of the telemetry data or when implementing the method, said value of the partition time interval Ω being a function of: the type of telemetry system from which the telemetry data are acquired, and/or the atmospheric conditions during acquisition of said telemetry data.

8. The method according to claim 4, comprising estimation of the fluctuations a of the wind speed over the reference time interval T according to equation (33): $\begin{array}{cc}\left[\mathrm{Math}33\right]& \\ \sigma =c.\sqrt{\frac{.Math.{\mathrm{Vhscal}}_{\mathrm{ave}}-{\mathrm{Vhvect}}_{\mathrm{ave}}.Math.}{{\mathrm{Vh}}_{\mathrm{ave}}}}& \mathrm{equation}33\end{array}$ in which c is a positive number and a is a zero or positive dimensionless number.

9. A method for processing telemetry data for estimating a wind speed, said method comprising an averaged projection comprising: a step (I) of vector reconstruction of at least two components of an instantaneous wind speed vector starting from projections of the instantaneous wind speed vector; and a step (J) of determining at least one average wind speed value over the time interval T by projection, over the time interval T, of the at least two components of the instantaneous wind speed vector reconstructed in step I.

10. The method according to claim 9 comprising: in step I, vector reconstruction, on the basis of equations (34) to (40), of the at least two components (U.sub.i, V.sub.i) or (V.sub.i, W.sub.i) or (U.sub.i, W.sub.i) among three components (U.sub.i, V.sub.i, W.sub.i) of the instantaneous wind speed vector; i is an integer comprised between 1 and N corresponding to the number of successive projections of the instantaneous wind speed vector over a time interval (T) called reference time interval, U, being the component of the instantaneous wind speed vector in a spatial direction (d1) extending in a spatial plane (p1), the component V, being the component of the instantaneous wind speed vector in a spatial direction (d2) extending in the spatial plane p1 and the component W, being the component of the average wind speed vector in a spatial direction (d3) orthogonal to the plane p1: $\begin{array}{cc}\left[\mathrm{Math}34\right]& \\ U_i=\frac{S_{\mathrm{Ni}}-S_{\mathrm{Si}}}{2.\sin \theta },\mathrm{or}& \mathrm{equation}34\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}35\right]& \\ U_i=\frac{S_{\mathrm{Ni}}-S_{\mathrm{Si}}}{2.\sin \theta .\cos \left(\frac{\alpha }{2}\right)},\mathrm{and}/\mathrm{or}& \mathrm{equation}39\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}36\right]& \\ V_i=\frac{S_{\mathrm{Ei}}-S_{\mathrm{Wi}}}{2.\sin \gamma },\mathrm{or}& \mathrm{equation}35\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}37\right]& \\ V_i=\frac{S_{\mathrm{Ei}}-S_{\mathrm{Wi}}}{2.\sin \gamma .\cos \left(\frac{\alpha }{2}\right)},\mathrm{and}/\mathrm{or}& \mathrm{equation}40\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}38\right]& \\ W_i=\frac{S_{\mathrm{Ei}}+S_{\mathrm{Wi}}}{2.\cos \gamma },\mathrm{or}& \mathrm{equation}36\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}39\right]& \\ W_i=\frac{S_{\mathrm{Ni}}+S_{\mathrm{Si}}}{2.\cos \theta },\mathrm{or}& \mathrm{equation}37\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}40\right]& \\ W_i=S_{\mathrm{Vi}},& \mathrm{equation}38\end{array}$ in which S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi, and Svi are projections of the instantaneous wind speed vector along, respectively, a first axis (a1), a second axis (a2), a third axis (a3), a fourth axis (a4) and a fifth axis (a5) merged with the direction d3, 6 is a non-zero angle formed between the axis a1 and a normal to the plane p1 and between the axis a2 and the normal to the plane p1 and γ is a non-zero angle formed between the axis a3 and the normal to the plane p1 and the axis a4 and the normal to the plane p1, the first and second axes a1 and a2 are included in a plane (p2), the third and fourth axes a3 and a4 are included in a plane (p3) and the planes p2 and p3 form a non-zero angle α between them, in step J, determination, on the basis of equations (41) to (42) and starting from the at least two reconstructed components of the instantaneous wind speed vector, of the at least one average wind speed value (Vh.sub.ave) in plane p1, p2 or p3 respectively over the reference time interval T: $\begin{array}{cc}\left[\mathrm{Math}41\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.1}}=\frac{1}{N-1}.Math.{.Math.}_{N-1}\left(\frac{\begin{array}{c}{U}_{i+1}.U_i+V_{i+1}.V_i+U_{i+1}.V_i.\cos \alpha +\\ V_{i+1}.U_i.\cos \alpha \end{array}}{\sqrt{{\left(U_i\right)}^2+{\left(V_i\right)}^2+2.U_i.V_i.\cos \alpha }}\right),& \mathrm{equation}41\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}42\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.2}}=\frac{1}{N-1}.Math.{.Math.}_{N-1}\left(\frac{U_{i+1}.U_i+W_{i+1}.W_i}{\sqrt{{\left(U_i\right)}^2+{\left(W_i\right)}^2}}\right),& \mathrm{equation}42\end{array}$ $\begin{array}{cc}\left[\mathrm{Math}43\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.3}}=\frac{1}{N-1}.Math.{.Math.}_{N-1}\left(\frac{V_{i+1}.V_i+W_{i+1}.W_i}{\sqrt{{\left(V_i\right)}^2+{\left(W_i\right)}^2}}\right),& \mathrm{equation}43\end{array}$

11. The method according to claim 2 comprising estimation of a wind direction (dir) in plane p1 according to equation (44): $\begin{array}{cc}\left[\mathrm{Math}44\right]& \\ \mathrm{Dir}={\tan }^{-1}\left(\frac{\mathrm{Vrec}}{\mathrm{Urec}}\right),& \mathrm{equation}44\end{array}$ in which tan.sup.−1 is the arc tangent function, the estimated wind direction is an angular value between the wind direction and the direction d1 and in which Vrec and Urec are each: a scalar value of a component of the wind speed in plane p1 over the reference time interval T, or an average vector speed of a component of the wind speed in plane p1 over the reference time interval T.

12. The method according to claim 1, comprising a step of measurement of the projections S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi and S.sub.Vi of the instantaneous wind speed vector by means of at least one measuring laser beam extending, respectively, along a first axis a1, a second axis a2, a third axis a3, a fourth axis a4 and a fifth axis a5.

13. The method according to claim 1 implemented by computer.

14. A data processing device comprising means arranged and/or programmed and/or configured for implementing the method according to claim 1.

15. A computer program comprising instructions which, when the program is executed by a computer, lead the latter to implement the method according to claim 1.

16. A recording medium: comprising instructions which, when they are executed by a computer, lead to implementation of the method according to claim 1.

17. A recording medium: comprising instructions which, when they are executed by a computer, on which the computer according to claim 15 is recorded.

Description

DESCRIPTION OF THE FIGURES

[0063] Other advantages and features of the invention will become apparent on reading the detailed description of implementations and embodiments, which are in no way limitative, and the following attached drawings:

[0064] FIG. 1 shows a diagrammatic representation of a slantwise view of an optical system for the acquisition of telemetry data used for implementing the method according to the invention,

[0065] FIG. 2 is a functional diagram of a first alternative of the method according to the invention,

[0066] FIG. 3 is a functional diagram of a second alternative of the method according to the invention,

[0067] FIG. 4 is a functional diagram of a third alternative of the method according to the invention.

DESCRIPTION OF THE EMBODIMENTS

[0068] As the embodiments described hereinafter are in no way limitative, it is possible, in particular, to consider variants of the invention comprising only a selection of the characteristics described, in isolation from the other characteristics described (even if this selection is isolated within a phrase comprising these other characteristics), if this selection of characteristics is sufficient to confer a technical advantage or to differentiate the invention with respect to the state of the prior art. This selection comprises at least one, preferably functional, characteristic without structural details, or with only a part of the structural details if this part alone is sufficient to confer a technical advantage or to differentiate the invention with respect to the state of the prior art.

[0069] FIG. 1 illustrates an example of an optical system 1 for acquiring telemetry data. According to the example in FIG. 1, the optical system 1 emits five optical measurement beams, each extending along a different axis a1, a2, a3, a4 and a5. As a non-limitative example, this system may be a LIDAR of the Lidar type with continuous measurement technology or Lidar with pulsed measurement technology. In this example, it is envisaged to estimate the average wind speed in plane p1. By extension of the method to an additional spatial dimension, a person skilled in the art will also be able to estimate the wind speed in space. According to the example in FIG. 1, axis a5 is vertical, axis a1 is inclined by an angle θ, which here is equal to 28°, relative to axis a5 towards the magnetic north, axis a2 is inclined by an angle θ relative to axis a5 towards the south, axis a3 is inclined by an angle γ, which here is equal to 28°, relative to axis a5 towards the east and axis a4 is inclined by an angle γ relative to axis a5 towards the west. The planes p2 and p3 form an angle α, which here is equal to 90°. In this example, the angle θ formed by axis a1 and axis a2 relative to axis a5 and the angle γ formed by axis a3 and axis a4 relative to axis a5 are identical. A person skilled in the art will also be able to adapt the method according to the invention to the case when these angles are different.

[0070] In metrology, wind may be characterized by its direction and its force or magnitude. In practice, wind is defined by a wind vector comprising three components (U, V, W), generally U represents the component of the wind vector along the axis from north to south, V represents the component of the wind vector along the axis from east to west and W represents the component of the wind vector on the axis normal to the surface of the earth at the point of measurement. This wind vector is measured by measuring the velocity of displacement of particles along each of the beams. The instantaneous values measured along each of the beams are projected components S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi and S.sub.vi of the wind vector. For this example, the system delivers 5 measurements S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi and S.sub.vi every 4 seconds. Thus, a measurement is available about every 0.8 seconds. Over a time interval Ω, there are therefore M, which is equal to division of Q (in seconds) by 4, sets of projected components S.sub.Ni S.sub.Si, S.sub.Ei, S.sub.Wi and S.sub.vi. It is then necessary to reconstruct the components (U, V, W) of the wind vector from the measured instantaneous projections S.sub.Ni, S.sub.Si, S.sub.Ei, S.sub.Wi and S.sub.Vi. In practice, the method for processing telemetry data for estimating a wind speed according to the invention may be implemented on the data measured in real time or on stored data such as measured stored data, statistical data or data that have not been measured (for example data from simulation).

[0071] Preferably, the instantaneous wind speed vector, on the basis of which the method for processing telemetry data for estimating a wind speed according to the invention is implemented, is measured by telemetry, for example by a LIDAR. Preferably, the method for processing telemetry data for estimating a wind speed according to the invention is implemented on the data relating to the instantaneous wind speed vector that is measured by telemetry, for example by a LIDAR.

[0072] A typical case of use of wind measurement is measurement of the power available for producing wind turbine energy, in this case the measurement interval T is typically a time interval of 10 min, which makes it possible to isolate the energy produced by the wind turbine. This interval is called the reference time interval T. The method for processing telemetry data according to the invention makes it possible to estimate the average wind speed over this reference interval.

[0073] According to a particular embodiment of the first alternative of the method according to the invention, the method of hybridization by temporal combination comprises: [0074] in step A, reconstruction, on the basis of equations 1 and 2 and the two components U.sub.Ω, V.sub.Ω of the wind vector along the north/south and east/west axis respectively, of the average wind speed vector over the partition time interval Ω; the component U.sub.Ω being the component of the average wind speed vector in a spatial direction d1 corresponding to the north/south axis extending in a spatial plane p1 corresponding to the plane tangential to the surface of the earth at the level of the measurement point and the component V.sub.Ω being the component of the average wind speed vector in a spatial direction d2 corresponding to the east/west axis extending in the spatial plane p1:

[00011]$\begin{array}{cc}\left[\mathrm{Math}1\right]& \\ U_{\Omega }=\frac{\frac{1}{M}.Math.{.Math.}_M\left(S_{\mathrm{Ni}}\right)-\frac{1}{M}.Math.{.Math.}_M\left(S_{\mathrm{Si}}\right)}{2.Math.\sin \theta },\mathrm{and}/\mathrm{or}& \mathrm{equation}1\end{array}$$\begin{array}{cc}\left[\mathrm{Math}3\right]& \\ V_{\Omega }=\frac{\frac{1}{M}.Math.{.Math.}_M\left(S_{\mathrm{Ei}}\right)-\frac{1}{M}.Math.{.Math.}_M\left(S_{\mathrm{Wi}}\right)}{2.Math.\sin \gamma },\mathrm{and}/\mathrm{or}& \mathrm{equation}2\end{array}$

in which i is an integer comprised between 1 and M, which here is equal to 25, since the partition time interval is 60 seconds and there is a measurement of the 5 projections of the instantaneous wind speed vector every 4 seconds, [0075] in step B, scalar reconstruction, on the basis of equation 6 and the two components of the average wind speed vector reconstructed in step A, U.sub.Ω and V.sub.Ω, of the average value V.sub.have of the horizontal wind speed over the reference time interval T in plane p1:

[00012]$\begin{array}{cc}\left[\mathrm{Math}8\right]{\mathrm{Vh}}_{\mathrm{ave}\text{.1}}=\frac{1}{Q}.Math.{.Math.}_Q(\sqrt{{\left(U_{\Omega }\right)}^2+{\left(V_{\Omega }\right)}^2+2.U_{\Omega }.V_{\Omega }.\cos \left(\alpha =90\right))}=\sqrt{{\left(U_{\Omega }\right)}^2+{\left(V_{\Omega }\right)}^2},& \mathrm{equation}8\end{array}$

Q is equal to T/Ω and corresponds to the number of the components U.sub.Ω and V.sub.Ω of the average wind speed over the partition time interval Ω included in the reference time interval T.

[0076] In practice, the value of the partition time interval Ω is constant or is modified during acquisition of the telemetry data, said value of the partition time interval Ω being a function of: [0077] the type of telemetry system from which the telemetry data are acquired, and/or [0078] the atmospheric conditions during acquisition of said telemetry data.

[0079] The value of Q may be adapted to the amplitude of variation of the direction and of the speed of the horizontal wind, indicated for example by calculating the standard deviation of the direction and the horizontal wind speed, or to the value of the estimated average wind speed.

[0080] FIG. 3 shows the functional schematic diagram of the method for processing telemetry data for estimating a wind speed according to a second alternative according to the invention. According to the second alternative, the method comprises hybridization by weighting comprising: [0081] a step C of vector reconstruction of at least two components of an instantaneous wind speed vector from projections of the instantaneous wind speed vector, [0082] a step D of vector reconstruction over a time interval T, called reference time interval, of at least two components of an average wind speed vector starting from N, which for a time interval of 10 minutes is 600 seconds and for acquisition of a set of projection of the instantaneous speed vector every 4 seconds is equal to 150, of the two reconstructed components, included over the reference time interval T, of the instantaneous wind speed vector, [0083] a step E of scalar reconstruction of at least one average wind speed value over the time interval T starting from the at least two reconstructed components of the average wind speed vector, [0084] a step F of determination of at least one value of the norm of the instantaneous wind speed starting from projections of the instantaneous wind speed vector, [0085] a step G of determination of at least one average value of the norm of the wind speed over the reference time interval T starting from at least one value of the norm of the reconstructed instantaneous wind speed, [0086] a step (H) of determination of at least one average value Vh.sub.ave of the wind speed over the time interval T by weighting of a sum of the at least one average wind speed value reconstructed in step E and of the at least one average wind speed value determined in step G.

[0087] According to a particular embodiment of the second alternative of the method according to the invention, the method of hybridization by weighting comprises: [0088] in step C, reconstruction, on the basis of the respective equations 11 and 12, of at least two components (U.sub.i, V) of an instantaneous wind speed vector along the north/south and east/west axis, respectively; i is an integer comprised between 1 and N corresponding to the number of successive projections of the instantaneous wind speed vector over a time interval T called reference time interval, U, being the component of the instantaneous wind speed vector in a spatial direction d1 corresponding to the north/south axis extending in a spatial plane p1 corresponding to the plane tangential to the surface of the earth at the level of the measurement point and the component V, being the component of the instantaneous wind speed vector in a spatial direction d2 corresponding to the east/west axis extending in the spatial plane p1:

[00013]$\begin{array}{cc}\left[\mathrm{Math}11\right]& \\ U_i=\frac{S_{\mathrm{Ni}}-S_{\mathrm{Si}}}{2.Math.\sin \theta },\mathrm{and}/\mathrm{or}& \mathrm{equation}11\end{array}$$\begin{array}{cc}\left[\mathrm{Math}13\right]& \\ V_i=\frac{S_{\mathrm{Ei}}-S_{\mathrm{Wi}}}{2.Math.\sin \gamma },\mathrm{and}/\mathrm{or}& \mathrm{equation}12\end{array}$ [0089] in step D, reconstruction, on the basis of equations 18 and 19, of the two components (Uvect.sub.N, Vvect.sub.N) of an average wind speed vector over the reference time interval T along the north/south and east/west axis, respectively; the component Uvect.sub.N being the component of the wind speed in the spatial direction d1 and the component Vvect.sub.N being the component of the wind speed in the spatial direction d2:

[00014]$\begin{array}{cc}\left[\mathrm{Math}18\right]& \\ {\mathrm{Uvect}}_N=\frac{1}{N}.Math.{.Math.}_N\left(U_i\right),\mathrm{and}& \mathrm{equation}18\end{array}$$\begin{array}{cc}[\mathrm{Math}19& \\ {\mathrm{Vvect}}_N=\frac{1}{N}.Math.{.Math.}_N\left(\mathrm{Vi}\right),\mathrm{and}& \mathrm{equation}19\end{array}$ [0090] in step E, reconstruction, on the basis of equation 21, of at least one scalar value (Vscal.sub.i) of the instantaneous wind speed; Vscal.sub.i corresponding to a temporal series of the scalar value of the instantaneous wind speed respectively in plane p1:

[00015]$\begin{array}{cc}\left[\mathrm{Math}21\right]{\mathrm{Vscal}}_{i\text{.1}}=\sqrt{{\left(U_i\right)}^2+{\left(V_i\right)}^2+2.U_i.V_i.\cos \left(\alpha =90\right)}=\sqrt{{\left(U_i\right)}^2+{\left(V_i\right)}^2},& \mathrm{equation}21\end{array}$ [0091] in step F, calculation, on the basis of equation 24 and starting from the scalar value Vscal.sub.i..sub.1 of the reconstructed instantaneous wind speed, of the average value Vhscal.sub.ave of the norm of the wind speed respectively in plane p1 over the reference time interval T

[00016]$\begin{array}{cc}\left[\mathrm{Math}24\right]& \\ {\mathrm{Vhscal}}_{\mathrm{ave}\text{.1}}=\frac{1}{N}.{.Math.}_N\left({\mathrm{Vscal}}_{i\text{.1}}\right),& \mathrm{equation}24\end{array}$ [0092] in step G, calculation, on the basis of equation 27 and starting from the two reconstructed components of the average wind speed vector, of the average wind speed value Vvect.sub.ave in plane p1 over the reference time interval T:

[00017]$\begin{array}{cc}\left[\mathrm{Math}27\right]& \\ {\mathrm{Vhvect}}_{\mathrm{ave}\text{.1}}=\sqrt{{\left({\mathrm{Uvect}}_N\right)}^2+{\left({\mathrm{Vvect}}_N\right)}^2+2.{\mathrm{Uvect}}_N.{\mathrm{Vvect}}_N.\cos \left(\alpha =90\right)}=\sqrt{{\left({\mathrm{Uvect}}_N\right)}^2+{\left({\mathrm{Vvect}}_N\right)}^2},& \mathrm{equation}27\end{array}$ [0093] in step H, calculation, on the basis of equation 30 and starting from the pairs of reconstructed wind speed values (Vhscal.sub.ave.1, Vhvect.sub.ave.1), of the weighted average wind speed value V.sub.have in plane p1 over the reference time interval T:

[Math 30]

[0094]
Vh.sub.ave1=(1−P).Vhscal.sub.ave1+P.Vhvect.sub.ave1,  equation 30

in which P is a dimensionless weighting factor comprised between 0 and 1.

[0095] The optimum value of P depends on: [0096] the type of telemetry system from which the telemetry data are acquired, and/or [0097] the atmospheric conditions during acquisition of said telemetry data.
The factor P is greater than 0.2 and/or less than 0.6, preferably greater than 0.3 and/or less than 0.5, more preferably equal to 0.33. Under the standard atmospheric conditions and for the telemetry system with the configuration presented in FIG. 1, the number that makes it possible to obtain the best estimates is approximately 0.33.

[0098] The method comprises estimation of the fluctuations a of the wind speed over the reference time interval T according to equation 33:

[00018]$\begin{array}{cc}\left[\mathrm{Math}33\right]& \\ \sigma =c.\sqrt{\frac{.Math.{\mathrm{Vhscal}}_{\mathrm{ave}}-{\mathrm{Vhvect}}_{\mathrm{ave}}.Math.}{{\mathrm{Vh}}_{\mathrm{ave}}}},& \mathrm{equation}33\end{array}$

in which it is a positive number and a is a zero or positive dimensionless number. This estimation is an approximation of the value of the standard deviation of the horizontal speed and of the direction that makes it possible to classify the measured wind flow in categories to be defined as high turbulence or low turbulence.

[0099] FIG. 4 shows the functional schematic diagram of the method for processing telemetry data for estimating a wind speed according to a third alternative according to the invention. According to the third alternative, the method comprises hybridization by averaged projection comprising: [0100] a step I of vector reconstruction of at least two components of an instantaneous wind speed vector starting from projections of the instantaneous wind speed vector, [0101] a step J of determination of at least one average value V.sub.have of the wind speed over the time interval T by projection, over the time interval T, of the at least two components of the instantaneous wind speed vector reconstructed in step I.

[0102] According to a particular embodiment of the third alternative of the method according to the invention, the method of hybridization by averaged projection comprises: [0103] in step I, vector reconstruction, on the basis of equations 34 and 35, of the two components (U.sub.i, V) of the instantaneous wind speed vector along the north/south and east/west axis, respectively; i is an integer comprised between 1 and N corresponding to the number of successive projections of the instantaneous wind speed vector over a time interval T called reference time interval, U, being the component of the instantaneous wind speed vector in a spatial direction d1 corresponding to the north/south axis extending in a spatial plane p1 corresponding to the plane at the surface of the earth at the level of the measurement point and the component V, being the component of the instantaneous wind speed vector in a spatial direction d2 corresponding to the east/west axis extending in the spatial plane p1:

[00019]$\begin{array}{cc}\left[\mathrm{Math}34\right]& \\ U_i=\frac{S_{\mathrm{Ni}}-S_{\mathrm{Si}}}{2.\sin \theta },\mathrm{and}/\mathrm{or}& \mathrm{equation}34\end{array}$$\begin{array}{cc}\left[\mathrm{Math}37\right]& \\ V_i=\frac{S_{\mathrm{Ei}}-S_{\mathrm{Wi}}}{2.\sin \gamma },\mathrm{and}/\mathrm{or}& \mathrm{equation}35\end{array}$ [0104] in step J, determination, on the basis of equation 41 and starting from the two reconstructed components of the instantaneous wind speed vector, of the average value V.sub.have of the wind speed respectively in plane p1 over the reference time interval T:

[00020]$\begin{array}{cc}\left[\mathrm{Math}41\right]& \\ {\mathrm{Vh}}_{\mathrm{ave}\text{.1}}=\frac{1}{N-1}.Math.{.Math.}_{N-1}\left(\frac{\begin{array}{c}{U}_{i+1}.U_i+V_{i+1}.V_i+U_{i+1}.V_i.\cos \left(\alpha =90\right)+\\ V_{i+1}.U_i.\cos \left(\alpha =90\right)\end{array}}{\sqrt{{\left(U_i\right)}^2+{\left(V_i\right)}^2+2.U_i.V_i.\cos \left(\alpha =90\right)}}\right)=\frac{1}{N-1}.Math.{.Math.}_{N-1}\left(\frac{U_{i+1}.U_i+V_{i+1}.V_i}{\sqrt{{\left(U_i\right)}^2+{\left(V_i\right)}^2}}\right),& \mathrm{equation}41\end{array}$

[0105] Of course, the invention is not limited to the examples that have just been described, and numerous adjustments may be made to these examples without departing from the scope of the invention. Thus, it is conceivable to combine variants or steps of the embodiments described above.

[0106] Moreover, the different characteristics, forms, variants and embodiments of the invention can be combined with one another in various combinations, provided that they are not incompatible or mutually exclusive.