Method for processing telemetry data to estimate wind speed

Abstract

The present invention relates to a method for processing telemetry data to estimate wind speed, the method comprising hybridization by temporal combination and / or hybridization by weighting and / or hybridization by averaged projection.

Description

Technical Field

[0001] The present invention relates to a method for estimating wind speed from telemetry data. The object of the present invention is, in particular, to determine the wind speed from telemetry data collected by a lidar, radar, or sodar system.

[0002] The present invention relates to a method based on the reconstruction of the mean wind speed based on data resulting from measurements made successively or continuously at different points in space by a lidar, radar, or sodar system.

Background Art

[0003] Methods for estimating wind speed by reconstructing data obtained from measurements made by a lidar, radar, or sodar system are known. In the state of the art, two alternative types of reconstruction have been found. The first method consists of considering the mean wind speed to be equal to the mean over a measurement time interval of one tenth of the norm of the instantaneous wind speed vector. This method is generally referred to as the "scalar method". The second method consists of considering the mean wind speed to be equal to the norm of the mean wind speed vector over a measurement interval of one tenth. This method is generally referred to as the "vector method".

[0004] When lidar measurements are made under conditions of strong atmospheric turbulence, i.e., conditions where there are significant variations in wind direction or wind speed over a ten-minute time interval, the estimates obtained by existing methods are significantly different from the actual wind speed. A cup anemometer is regarded as a standard tool for measuring wind speed.

[0005] In fact, under conditions of low atmospheric turbulence, i.e., conditions where the variations in wind direction or wind speed are small over a ten-minute time interval, the relative error between the wind speed estimated by the state-of-the-art method and the wind speed measured by a cup anemometer is less than 1%. However, under conditions of strong atmospheric turbulence, the relative error between the wind speed estimated by the state-of-the-art method and the wind speed measured by a cup anemometer can reach an absolute value of 4%.

[0006] It is known in the state of the art that a standard tool for determining wind speed is a cup anemometer. The determination of wind speed is closely related to the measuring device used and the transfer function used for the determination.

[0007] The present invention in particular proposes a method for estimating wind speed that enables overcoming the drawbacks of state-of-the-art methods, and / or proposes a method for estimating a wind speed such that the determined wind speed value is as close as possible to the wind speed value measured by a cup anemometer, and / or proposes a method for estimating a wind speed such that the error in the determined wind speed value is less than 1% under conditions of strong atmospheric turbulence and under conditions of low atmospheric turbulence, and / or proposes a method for estimating a wind speed that is not sensitive to the atmospheric conditions under which the measurement is made for the purpose. **Disclosure of the Invention**

[0008] For this purpose, a method for processing telemetry data to estimate wind speed has been proposed.

[0009] According to a first alternative form, the method starts from the continuous projection over time of the instantaneous wind speed vector and comprises a step (A) of vector reconstruction of at least two components of the mean wind speed vector over a time interval (Ω), called a split time interval, and starts from the number (T / Ω) of at least two components of the mean wind speed vector reconstructed in step A and comprises a step (B) of scalar reconstruction of at least one mean wind speed value (Vh ave ) over a time interval (T), called a reference time interval, including a hybridization by a temporal combination, where 2Ω is less than or equal to T and T / Ω corresponds to the number of at least two components of the mean wind speed over the split time interval Ω comprised in the reference time interval T.

[0010] The method of hybridization by temporal combination according to the first alternative form is In step A, based on formulas (1) to (7), at least two of the three components (U Ω , V Ω , W Ω ) of the mean wind speed vector over the divided time interval Ω (U Ω , V Ω ) or (V Ω , W Ω ) or (U Ω , W Ω ) are reconstructed, where the component U Ω is the component of the mean wind speed vector in the spatial direction (d1) extending in the spatial plane (p1), the component V Ω is the component of the mean wind speed vector in the spatial direction (d2) extending in the spatial plane p1, and the component W Ω is the component of the mean wind speed vector in the spatial direction (d3) orthogonal to the plane p1.

Number

Number

Number

Number

Number

Number

Number

[0011] Further, the hybridization method by temporal combination according to the first alternative form is In step B, based on equations (8) to (10) and at least two components of the mean wind speed vector reconfigured in step A, at least one mean wind speed value (Vh ave ) for a scalar reconstruction over a reference time interval T in each of the planes p1 or p2 or p3 [Number] [Number] [Number] Q is an integer included between 1 and (T / Ω) corresponding to the number of at least two components of the average wind speed over the divided time interval Ω included in the reference time interval T Step may be included.

[0012] The value of the divided time interval Ω may be constant or modified during the acquisition of the telemetry data, and said value of the divided time interval Ω is a function of the type of telemetry system from which the telemetry data was acquired, and / or the atmospheric conditions during the acquisition of said telemetry data .

[0013] According to a second alternative form, the method is a step (C) of vector reconstruction of at least two components of the instantaneous wind speed vector starting from the projection of the instantaneous wind speed vector, and a step (D) of vector reconstruction over a time interval (T) called the reference time interval of at least two components of the average wind speed vector starting from the number N of at least two components of the instantaneous wind speed vector reconstructed in step C included over the reference time interval T, and a step (E) of scalar reconstruction of at least one instantaneous wind speed value starting from at least two components of the average wind speed vector reconstructed in step C, and a step (F) of determining at least one average wind speed value starting from at least one instantaneous wind speed value reconstructed in step E, and a step (G) of determining at least one average wind speed value over the reference time interval T starting from 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 ave ) over the time interval T by weighting the sum of at least one average wind speed value reconstructed in step F and at least one average wind speed value determined in step G and including hybridization by weighting.

[0014] Preferably, the method of hybridization by weighting according to the second alternative form is In step C, based on each of the formulas (11) to (17), at least two of the three components (U i , V i , W i ) of the instantaneous wind speed vector, namely, a vector reconstruction step of (U i , V i ) or (V i , W i ) or (U i , W i ), where i is an integer included between 1 and N corresponding to the number of consecutive projections of the instantaneous wind speed vector over the reference time interval T, U i is the component of the instantaneous wind speed vector in the spatial direction (d1) extending in the spatial plane (p1), the component V i is the component of the instantaneous wind speed vector in the spatial direction (d2) extending in the spatial plane p1, and the component W i is the component of the mean wind speed vector in the spatial direction (d3) orthogonal to the plane p1.

Number

Number

Number

Number

Number

Number

Number

[0015] Further, the method of hybridization by temporal combination according to the second alternative form is In step D, based on equations (18) - (20), at least two components (Uvect N , Vvect N ) or (Vvect N , Wvect N ) or (Uvect N , Wvect N ) of the average wind speed vector over the reference time interval T are re - constructed in a vector, where the component Uvect N is the component of the wind speed in the spatial direction d1, the component Vvect N is the component of the wind speed in the spatial direction d2, and the component Wvect N is the component of the wind speed in the spatial direction d3. [Number] [Number] [Number] a step that is; In step E, based on equations (21) to (23), a step of scalar reconstruction of at least one value (Vscal i ) of the instantaneous wind speed, where Vscal i corresponds to the time series of the instantaneous wind speed values in each of the planes p1, p2, or p3,

Number

Number

Number

Number

Number

Number

Number

Number

Number

Number

Number

Number

[0016] The coefficient P may be 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.

[0017] The value of the coefficient P may be constant or modified during the acquisition of the telemetry data or when implementing the method, and the value of the division time interval Ω is The type of the telemetry system from which the telemetry data was obtained, and / or the atmospheric conditions during the acquisition of the telemetry data is a function of.

[0018] The method of hybridization by weighting according to a second alternative form may include estimating the variation σ of the wind speed over the reference time interval T according to Equation (33), [Number] where c is a positive number and σ is a zero or positive dimensionless number.

[0019] According to a third alternative form, the method step (I) of vector reconstruction of at least two components of the instantaneous wind speed vector starting from the projection of the instantaneous wind speed vector, and step (J) of determining at least one average wind speed value over the time interval T by the projection of at least two components of the instantaneous wind speed vector reconstructed in step (I) over the time interval T including an averaging projection.

[0020] Preferably, the method of hybridization by averaging projection according to the third alternative form in step (I), based on Equations (34) to (40), at least two of the three components (U i , V i , W i ) of the instantaneous wind speed vector (U i , V i ) or (V i , W i ) or (U i , W i ) vector reconstruction step, where i is an integer included between 1 and N corresponding to the number of consecutive projections of the instantaneous wind speed vector over a time interval (T) called the reference time interval, U i is the component of the instantaneous wind speed vector in the spatial direction (d1) extending in the spatial plane (p1), and the component Vi is the component of the instantaneous wind speed vector in the spatial direction (d2) extending within the spatial plane p1, and the component W i is the component of the average wind speed vector in the spatial direction (d3) orthogonal to the plane p1,

Number

Number

Number

Number

Number

Number

Number

[0021] Furthermore, the hybridization method by temporal combination according to the third alternative form is In step J, starting from at least two reconstructed components of the instantaneous wind speed vector, based on formulas (41) to (42), at least one average wind speed value (Vh ave ) in each of the planes p1, p2, or p3 over the reference time interval T, which is a step of determining

Number

Number

Number

[0022] The hybridization method by averaging projection according to the third alternative form may include the estimation of the wind direction (dir) in the plane p1 by formula (44),

Number

[0023] The method according to any one of the first, second, and / or third alternative forms may include the step of measuring the projections S Ni , S Si , S Ei , S Wi , and S Vi of the instantaneous wind speed vector by at least one measurement laser beam extending along each of the respective axes a1, a2, a3, a4, and a5.

[0024] The method according to any one of the first, second, and / or third alternative forms may be implemented by a computer.

[0025] According to the present invention, there is also proposed a data processing device comprising means arranged and / or programmed and / or configured to implement the method according to any one of the first, second, and / or third alternative forms.

[0026] According to the present invention, there is also proposed a computer program comprising instructions which, when executed by a computer, cause the computer to implement the method according to any one of the first, second, and / or third alternative forms.

[0027] According to the present invention, when executed by a computer, comprising instructions which cause the computer to implement the method according to any one of the first, second, and / or third alternative forms, and / or a recording medium on which is recorded a computer program according to the present invention, is also proposed. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] Other advantages and features of the present invention will become apparent from the following detailed description of non-limiting embodiments and embodiments, as well as from the following attached drawings.

Figure 1

Figure 2

Figure 3

Figure 4

Embodiments for Carrying Out the Invention

[0029] The embodiments described below are in no way limiting. In particular, when this selection of characteristics is sufficient to provide a technical advantage or to distinguish the present invention with respect to the state of the prior art, (even if this selection is separated within a phrase containing these other characteristics), variants of the present invention may be considered to include only the selection of the characteristics described, separately from the other characteristics described. This selection may include at least one, preferably functional, characteristic without structural details or with only some of the structural details, if this selection alone is sufficient to provide a technical advantage or to distinguish the present invention with respect to the state of the prior art.

[0030] FIG. 1 shows an example of an optical system 1 for acquiring telemetry data. According to the example of FIG. 1, the optical system 1 emits five optical measurement beams each extending along different axes a1, a2, a3, a4, and a5. As a non-limiting example, this system may be a Lidar type LIDAR with continuous measurement technology or a Lidar with pulse measurement technology. In this example, it is assumed that the average wind speed in the plane p1 is estimated. By extending this method to additional spatial dimensions, one skilled in the art can also estimate the wind speed in space. According to the example of FIG. 1, the axis a5 is vertical, the axis a1 is inclined by an angle θ (here equal to 28°) with respect to the axis a5 towards magnetic north, the axis a2 is inclined by an angle θ with respect to the axis a5 towards the south, the axis a3 is inclined by an angle γ (here equal to 28°) with respect to the axis a5 towards the east, and the axis a4 is inclined by an angle γ with respect to the axis a5 towards the west. The planes p2 and p3 form an angle α here equal to 90°. In this example, the angle θ formed by the axes a1 and a2 with respect to the axis a5, and the angle γ formed by the axes a3 and a4 with respect to the axis a5 are the same. One skilled in the art can also adapt the method according to the present invention when these angles are different.

[0031] In measurement, wind can be characterized by its direction and its force or magnitude. In practice, wind is defined by a wind vector containing 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 perpendicular to the ground surface at the measurement point. This wind vector is measured by measuring the displacement velocity of particles along each of the beams. The instantaneous values measured along each of the beams are the projection components S Ni 、S Si 、S Ei 、S Wi 、and S Vi . In this example, the system measures five measurement values S Ni 、S Si 、S Ei 、S Wi, and S Vi are delivered. Thus, measurements are available approximately every 0.8 seconds. Thus, over the time interval Ω, M projection components S equal to the value obtained by dividing Ω (in seconds) by 4 Ni , S Si , S Ei , S Wi , and S Vi exist. Then, the components (U, V, W) of the wind vector need to be reconstructed from the measured instantaneous projections S Ni , S Si , S Ei , S Wi , and S Vi . In practice, the method for processing telemetry data according to the present invention to estimate wind speed can be implemented for real-time measured data, or stored data such as measured stored data, statistical data, or unmeasured data (e.g., data from simulations).

[0032] Preferably, the instantaneous wind speed vector based on the implementation of the method for processing telemetry data according to the present invention to estimate wind speed is measured by telemetry, for example, by LIDAR. Preferably, the method for processing telemetry data according to the present invention is implemented for data regarding the instantaneous wind speed vector measured by telemetry, for example, by LIDAR.

[0033] A typical case of using wind measurements is the measurement of the power available to generate wind turbine energy, in which case the measurement interval T is typically a 10-minute time interval, which enables the separation of the energy generated by the wind turbine. This interval is called the reference time interval T. The method for processing telemetry data according to the present invention enables the estimation of the average wind speed over this reference interval. According to a specific embodiment of the first alternative form of the method according to the present invention, the method of hybridization by temporal combination is In step A, two components U Ω and V Ω of the wind vector along the north / south axis and the east / west axis respectively, based on which, a step of reconstructing the average wind speed vector over the divided time interval Ω, where the component U Ω is the component of the average wind speed vector in the spatial direction d1 corresponding to the north / south axis extending in the spatial plane p1 corresponding to the plane tangent to the ground at the level of the measurement point, and the component V Ω is the component of the average wind speed vector in the spatial direction d2 corresponding to the east / west axis extending in the spatial plane p1,

Number

Number

Number

[0034] Actually, the value of the divided time interval Ω is constant or corrected during the acquisition of the telemetry data, and the said value of the divided time interval Ω is The type of the telemetry system from which the telemetry data was acquired, and / or the atmospheric conditions during the acquisition of the telemetry data is a function of.

[0035] The value of Ω can be adapted to, for example, the amplitude of the variations in the direction and speed of the horizontal wind, as indicated by calculating the standard deviation of the direction and horizontal wind speed, or the value of the estimated mean wind speed. FIG. 3 shows a functional schematic diagram of a method for processing telemetry data to estimate wind speed according to a second alternative form of the present invention. According to the second alternative form, the method comprises step C of vector reconstruction of at least two components of the instantaneous wind speed vector from the projection of the instantaneous wind speed vector, and step D of vector reconstruction over a time interval T called the reference time interval of at least two components of the mean wind speed vector, starting from N of the two reconstructed components of the instantaneous wind speed vector included over the reference time interval T (the time interval of one tenth is 600 seconds, and the acquisition of a set of projections of the instantaneous speed vector every 4 seconds is equal to 150), and step E of scalar reconstruction of at least one mean wind speed value over a time interval T, starting from at least two reconstructed components of the mean wind speed vector, and step F of determination of at least one value of the norm of the instantaneous wind speed, starting from the projection of the instantaneous wind speed vector, and step G of determination of at least one average value of the norm of the wind speed over a reference time interval T, starting from at least one value of the norm of the reconstructed instantaneous wind speed, and step (H) of determination of at least one average value Vh of the wind speed over a time interval T by weighting the sum of at least one average wind speed value reconstructed in step E and at least one average wind speed value determined in step G ave and includes weighted hybridization.

[0036] According to a particular embodiment of the second alternative form of the method according to the present invention, the method of weighted hybridization is In step C, based on each of Equation 11 and Equation 12, at least two components (U i , V i ) of the instantaneous wind speed vector along the north / south axis and the east / west axis respectively, where i is an integer between 1 and N corresponding to the number of consecutive projections of the instantaneous wind speed vector over a time interval T called the reference time interval, U i is the component of the instantaneous wind speed vector in the spatial direction d1 corresponding to the north / south axis extending in the spatial plane p1 corresponding to the plane tangent to the ground at the level of the measurement point, and the component V i is the component of the instantaneous wind speed vector in the spatial direction d2 corresponding to the east / west axis extending in the spatial plane p1,

Number

Number

Number

Number

[0037] The optimal value of P depends on the type of the telemetry system from which the telemetry data was acquired, and / or the atmospheric conditions during the acquisition of the telemetry data and.

[0038] The coefficient 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. In the case of the telemetry system of the configuration presented in FIG. 1 under standard atmospheric conditions, the number that enables obtaining the best estimate is approximately 0.33.

[0039] The method includes an estimation of the variation σ of the wind speed over a reference time interval T according to Equation 33,

Number

[0040] FIG. 4 shows a functional schematic diagram of a method for processing telemetry data to estimate the wind speed according to a third alternative form of the present invention. According to the third alternative form, the method starts with the projection of the instantaneous wind speed vector, a step I of vector reconstruction of at least two components of the instantaneous wind speed vector, and at least one average value Vh of the wind speed over a time interval T by the projection of at least two components of the instantaneous wind speed vector reconstructed in step I over a time interval T ave a step J of determination of and includes hybridization by averaging projection. According to a specific embodiment of the third alternative form of the method according to the present invention, the method of hybridization by averaging projection in step I, based on Equation 34 and Equation 35, two components (U i , V i ) of the instantaneous wind speed vector along the north / south axis and the east / west axis, respectively, a step of vector reconstruction, where i is an integer included between 1 and N corresponding to the number of consecutive projections of the instantaneous wind speed vector over a time interval T called the reference time interval, and U iis a component of the instantaneous wind speed vector in the spatial direction d1 corresponding to the north / south axis extending within the spatial plane p1 corresponding to the ground plane at the level of the measurement point, and the component V i is a component of the instantaneous wind speed vector in the spatial direction d2 corresponding to the east / west axis extending within the spatial plane p1, [Number] [Number] is a step, and In step J, based on Equation 41, starting from the two reconstructed components of the instantaneous wind speed vector, the average value Vh of the respective wind speeds in the plane p1 over the reference time interval T ave is a step of determination, [Number] is a step and includes.

[0041] Of course, the present invention is not limited to the examples described here, and many adjustments can be made to these examples without departing from the scope of the present invention. Therefore, it is conceivable to combine the modifications or steps of the above-described embodiments.

[0042] Furthermore, different characteristics, forms, modifications, and embodiments of the present invention can be combined with each other in various combinations as long as they are not incompatible or mutually exclusive.

Claims

1. A method for processing telemetry data to estimate wind speed, the method comprising: a step (A) of vector reconstruction of at least two components of an average wind speed vector over a time interval (Ω), called a split time interval, starting from a continuous projection of an instantaneous wind speed vector over time; Starting from the number (T / Ω) of said at least two components of said mean wind speed vector reconfigured in step A, at least one mean wind speed value (Vh ａｖｅ ) of scalar reconstruction over a time interval (T) called the reference time interval (B) and including hybridization by temporal combination, where 2Ω is less than or equal to T, and T / Ω corresponds to the number of said at least two components of said average wind speed over said split time interval Ω included in said reference time interval T; a method.

2. In step A, based on formulas (1) to (7), at least two of the three components (U Ω , V Ω , W Ω ) of the average wind speed vector over the division time interval Ω are reconstructed, where the component U Ω is the component of the average wind speed vector in the spatial direction (d1) extending in the spatial plane (p1), the component V Ω is the component of the average wind speed vector in the spatial direction (d2) extending in the spatial plane p1, and the component W Ω is the component of the average wind speed vector in the spatial direction (d3) orthogonal to the plane p1. The reconstruction step is one of (U Ω , V Ω ) or (V Ω , W Ω ) or (U Ω , W Ω ). 【Number 1】 【Number 2】 【Number 3】 【Number 4】 [Number 5] 【Number 6】 【Number 7】 wherein i is an integer included between 1 and M corresponding to S, S, S, S, S, and S over the divided time interval Ω of the instantaneous wind speed vector over time Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S ｖｉ ; S Ｎｉ 、 S Ｓｉ 、 S Ｅｉ 、 S Ｗｉ 、 and S ｖｉ are the respective projections of the instantaneous wind speed vector along the first axis (a1), second axis (a2), third axis (a3), fourth axis (a4), and fifth axis (a5) each fused with the direction d3, θ is a non-zero angle formed between the axis a1 and the normal of the plane p1 and between the axis a2 and the normal of the plane p1, γ is a non-zero angle formed between the axis a3 and the normal of the plane p1 and between the axis a4 and the normal of 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 α therebetween a step; In step B, based on the at least two components of the mean wind speed vector reconfigured in equations (8) to (10) and step A, at least one mean wind speed value (Vh ａｖｅ ) over the reference time interval T in each of planes p1, p2, or p3, which is a step of scalar reconstruction 【Number 8】 【Number 9】 【Number 10】 where Q is an integer included between 1 and (T / Ω) corresponding to the number of said at least two components of said average wind speed over said split time interval Ω included in said reference time interval T a step and The method according to claim 1, comprising.

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

4. A method for processing telemetry data to estimate wind speed, the method comprising: a step (C) of vector reconstruction of at least two components of said instantaneous wind speed vector, starting from a projection of the instantaneous wind speed vector; a step (D) of vector reconstruction of at least two components of an average wind speed vector over a time interval (T), called a reference time interval, starting from the number N of said at least two components of said instantaneous wind speed vector reconstructed in step C and included over said reference time interval T; a step (E) of scalar reconstruction of at least one instantaneous wind speed value, starting from said at least two components of said average wind speed vector reconstructed in step C; a step (F) of determining at least one average wind speed value, starting from said at least one instantaneous wind speed value reconstructed in step E; a step (G) of determining at least one average wind speed value over said reference time interval T, starting from said at least two components of said average wind speed vector reconstructed in step D; By weighting the sum of the at least one average wind speed value reconfigured in step F and the at least one average wind speed value determined in step G, at least one average wind speed value (Vh ａｖｅ ) over the time interval T is determined (step H) and including hybridization by weighting.

5. In step C, based on each of the formulas (11) to (17), three components (U ｉ , V ｉ , W ｉ ) of the instantaneous wind speed vector, at least two components (U ｉ , V ｉ ) or (V ｉ , W ｉ ) or (U ｉ , W ｉ ) are subjected to vector reconstruction. Here, i is an integer included between 1 and N corresponding to the number of consecutive projections of the instantaneous wind speed vector over the reference time interval T. U ｉ is the component of the instantaneous wind speed vector in the spatial direction (d1) extending in the spatial plane (p1). The component V ｉ is the component of the instantaneous wind speed vector in the spatial direction (d2) extending in the spatial plane p1. The component W ｉ is the component of the mean wind speed vector in the spatial direction (d3) orthogonal to the plane p1. 【Number 11】 【Number 12】 【Number 13】 【Number 14】 【Number 15】 【Number 16】 【Number 17】 wherein S Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S Ｖｉ are projections of the instantaneous wind speed vector along a first axis (a1), a second axis (a2), a third axis (a3), a fourth axis (a4), and a fifth axis (a5) respectively, each of which is fused with the direction d3, θ is a non-zero angle formed between the axis a1 and the normal of the plane p1 and between the axis a2 and the normal of the plane p1, γ is a non-zero angle formed between the axis a3 and the normal of the plane p1 and between the axis a4 and the normal of 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 α therebetween a step; In step D, based on formulas (18) to (20), at least two components (Uvect Ｎ , Vvect Ｎ ) or (Vvect Ｎ , Wvect Ｎ ) or (Uvect Ｎ , Wvect Ｎ ) of the average wind speed vector over the reference time interval T are subjected to vector reconstruction, where the component Uvect Ｎ is the component of the wind speed in the spatial direction d1, the component Vvect Ｎ is the component of the wind speed in the spatial direction d2, and the component Wvect Ｎ is the component of the wind speed in the spatial direction d3, 【Number 18】 【Number 19】 【Number 20】 a step that is; In step E, based on equations (21) to (23), a step of scalar reconstruction of at least one value (Vscal ｉ ) of the instantaneous wind speed, where Vscal ｉ corresponds to the time series of the instantaneous wind speed values in planes p1, p2, or p3, respectively, 【Number 21】 【Number 22】 【Number 23】 a step that is; In step F, based on formulas (24) to (26), the value Vscal of the instantaneous wind speed reconfigured in step E ｉ．１ or Vscal ｉ．２ or Vscal ｉ．３ starting from, the step of determining at least one average wind speed value (Vhscal ａｖｅ ) in each of the planes p1, p2, or p3 over the reference time interval T, 【24 Points】 【Number 25】 [Number 26] a step that is; In step G, starting from at least two reconfigured components of the mean wind speed vector, based on equations (27) to (29), at least one mean wind speed value (Vvect ａｖｅ ) in each of the planes p1, p2, or p3 over the reference time interval T is determined, 【Number 27】 【Number 28】 【Number 29】 a step that is; In step H, based on equations (30) to (32), starting from a pair of the reconstructed wind speed values (Vhscal ａｖｅ．１ , Vhvect ａｖｅ．１ ), or (Vhscal ａｖｅ．２ and Vhvect ａｖｅ．２ ), or (Vhscal ａｖｅ．３ , Vhvect ａｖｅ．３ ), the step of calculating at least one weighted average wind speed value (Vh ａｖｅ ) in each of the planes p1, p2, or p3 over the reference time interval T, 【30 numbers】 【Number 31】 【Number 32】 where P is a dimensionless weighting coefficient included between 0 and 1 step and The method according to claim 4, comprising

6. The method according to claim 5, wherein the coefficient P is greater than 0.2 and / or less than 0.

6.

7. The value of the coefficient P is constant or corrected during the acquisition of the telemetry data or when the method is implemented, and the value of the division time interval Ω is the type of the telemetry system from which the telemetry data is acquired, and / or the atmospheric conditions during the acquisition of the telemetry data The method according to any one of claims 5 to 6, which is a function of

8. including the estimation of the variation σ of the wind speed over the reference time interval T according to equation (33), 【Number 33】 where c is a positive number and σ is a zero or positive dimensionless number, The method according to any one of claims 4 to 7.

9. A method for processing telemetry data to estimate wind speed, the method comprising including an averaging projection, the averaging projection comprising Derived from the projection of the instantaneous wind speed vector and from equations (34) to (40), three components (U ｉ , V ｉ , W ｉ ) of the instantaneous wind speed vector, and at least two of these components (U ｉ , V ｉ ) or (V ｉ , W ｉ ) or (U ｉ , W ｉ ) in step I of vector reconstruction, where i is an integer included between 1 and N corresponding to the number of consecutive projections of the instantaneous wind speed vector over a time interval (T) called the reference time interval, U ｉ is said component of the instantaneous wind speed vector in a spatial direction (d1) extending in a spatial plane (p1), said component V ｉ is said component of the instantaneous wind speed vector in a spatial direction (d2) extending in said spatial plane p1, and said component W ｉ is said component of the mean wind speed vector in a spatial direction (d3) orthogonal to said plane p1, 【Number 34】 【Number 35】 【Number 36】 【Number 37】 【Number 38】 【Number 39】 【Number 40】 wherein S Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S Ｖｉ are projections of the instantaneous wind speed vector along the first axis (a1), the second axis (a2), the third axis (a3), the fourth axis (a4), and the fifth axis (a5) respectively, which are fused with the direction d3, θ is a non-zero angle formed between the axis a1 and the normal of the plane p1 and between the axis a2 and the normal of the plane p1, γ is a non-zero angle formed between the axis a3 and the normal of the plane p1 and between the axis a4 and the normal of 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 α therebetween step I, and Step J of determining at least one average wind speed value (Vh ａｖｅ ) in each of the planes p1, p2, or p3 over the reference time interval T by projection based on equations (41) to (43) over the time interval T of the at least two components of the instantaneous wind speed vector reconfigured in step I 【Number 41】 【Number 42】 【Number 43】 step J which is The method comprising

10. including the estimation of the wind direction (dir) in the plane p1 according to equation (44), 【Number 44】 where tan －１ is the arctangent function, the wind direction Dir estimated from Equation 44 is the angular value between the wind direction and the direction d1, and Vrec and Urec are each the scalar value of the component of the wind speed in the plane p1 over the reference time interval T, or the average vector velocity of the component of the wind speed in the plane p1 over the reference time interval T being The method according to any one of claims 2, 5, and 9.

11. The projection S of the instantaneous wind speed vector by at least one measurement 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 Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S Ｖｉ The method according to any one of claims 2, 5, and 9, comprising the step of measuring

12. The method according to any one of claims 1 to 11, implemented by a computer.

13. A data processing device comprising means arranged and / or programmed and / or configured to implement the method according to any one of claims 1 to 12.

14. A computer program comprising instructions which, when executed by a computer, cause the computer to implement the method according to any one of claims 1 to 12.

15. including instructions which, when executed by a computer, cause the computer to implement the method according to any one of claims 1 to 12, and / or on which the computer program according to claim 14 is recorded, A recording medium.

Images