Method for processing telemetry data to estimate wind speed - Patent Application 20070122967

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. In particular, the object of the present invention is to determine wind speed from telemetry data collected by a lidar, radar, or SODA system.

[0002] The present invention relates to a method for reconstructing average wind speed based on data resulting from measurements taken sequentially or continuously at different points in space by a LiDAR, radar, or SODA system. [Background technology]

[0003] Methods for estimating wind speed are known by reconstructing data obtained from measurements performed by LiDAR, radar, or SODA systems. In the latest technology, two alternative types of reconstruction have been found. The first method consists of considering the average wind speed as equal to the average of the norm of the instantaneous wind speed vector over a 10-minute measurement time interval. This method is generally called the "scalar method". The second method consists of considering the average wind speed as equal to the norm of the average wind speed vector over a 10-minute measurement interval. This method is generally called the "vector method".

[0004] Under conditions of strong atmospheric turbulence, i.e., when LiDAR measurements are performed under conditions of significant fluctuation in wind direction or speed over 10-minute time intervals, estimates obtained by existing methods differ significantly from actual wind speeds. Cup anemometers are considered a standard tool for measuring wind speed.

[0005] In practice, under conditions of low atmospheric turbulence, i.e., when wind direction or speed fluctuates little over a 10-minute time interval, the relative error between wind speed estimated by state-of-the-art methods and wind speed measured by a cup anemometer is less than 1%. However, under conditions of strong atmospheric turbulence, the relative error between wind speed estimated by state-of-the-art methods and wind speed measured by a cup anemometer can reach an absolute value of 4%.

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

[0007] The invention, in particular, To propose a method for estimating wind speed that will enable overcoming the shortcomings of the latest technological methods, and / or To propose a method for estimating wind speeds that are as close as possible to the wind speeds measured by a cup anemometer, and / or To propose a method for estimating wind speed such that the error in the determined wind speed value is less than 1% under conditions of strong and low atmospheric turbulence, and / or To propose a method for estimating wind speed that is not sensitive to the atmospheric conditions under which measurements are taken. The purpose is to achieve this. [Disclosure of the Invention]

[0008] For this purpose, methods have been proposed for processing telemetry data to estimate wind speed.

[0009] According to the first alternative form, the method is: Step (A) involves vector reconstruction of at least two components of the average wind speed vector over time intervals (Ω) called divided time intervals, starting from a continuous projection of the instantaneous wind speed vector over time, Starting with the number of at least two components (T / Ω) of the average wind speed vector reconstructed in step A, at least one average wind speed value (Vh) over a time interval (T) called the reference time interval. ave Step (B) of the scalar reconstruction of ) and This includes hybridization by temporal combinations, 2Ω is less than or equal to T, and T / Ω corresponds to the number of at least two components of the average wind speed over the divided time interval Ω contained 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 components (U Ω , V Ω , W Ω ) of the three components of the average wind speed vector over the divided time interval Ω, namely (U Ω , V Ω ) or (V Ω , W Ω ) or (U Ω , W Ω ) are reconstructed. 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.

Number

Number

Number

Number

Number

Number

Number

[0011] Furthermore, the hybridization method using temporal combinations of the first alternative form is: In step B, based on equations (8) to (10) and at least two components of the average wind speed vector reconstructed in step A, at least one average wind speed value (Vh) over a reference time interval T in each of the planes p1, p2, or p3 is calculated. ave ) is a scalar reconstruction step,

number

number

number

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

[0013] According to the second alternative form, the method is: 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, Step (D) of vector reconstruction over a time interval (T) called the reference time interval, starting with a number N of at least two components of the instantaneous wind speed vector reconstructed in step C, which are included over the reference time interval T, and of at least two components of the mean wind speed vector, Step (E) is a 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, Step (F) determines at least one average wind speed value, starting from at least one instantaneous wind speed value reconstructed in step E. Step (G) determines at least one average wind speed value over a reference time interval T, starting from at least two components of the average wind speed vector reconstructed in step D. The weighting of 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 results in at least one average wind speed value (Vh) over time interval T. ave Step (H) to determine This includes weighted hybridization.

[0014] Preferably, the hybridization method by weighting using a second alternative form is: In step C, based on equations (11) to (17), the three components of the instantaneous wind speed vector (U i , V i , W i ) at least two of the components (U i , V i ) or (V i , W i ) or (U i , W i A step of vector reconstruction of ), where i is an integer between 1 and N corresponding to the number of continuous projections of instantaneous wind speed vectors over a reference time interval T, and U i This is a component of the instantaneous wind velocity vector in the spatial direction (d1) extending within the spatial plane (p1), and component V i This is the component of the instantaneous wind velocity vector in the spatial direction (d2) extending within the spatial plane p1, and component W i This is the component of the average wind speed vector in the spatial direction (d3) perpendicular to the plane p1,

number

number

number

number

number

number

number

[0015] Furthermore, the second alternative method of hybridization by temporal combination is: In step D, based on equations (18) to (20), at least two components of the average wind speed vector over the reference time interval T (Uvect N Vvect N ) or (Vvect N Wvect N ) or (Uvect N Wvect N A step in vector reconstruction of the component Uvect N This is the component of the wind velocity in the spatial direction d1, and the component Vvect N This is the component of wind velocity in the spatial direction d2, and the component Wvect N This is the component of wind velocity in the spatial direction d3,

number

number

number

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, and more preferably equal to 0.33.

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

[0018] The second alternative hybridization method, weighted by a different form, may include estimating the wind speed variation σ over a reference time interval T using equation (33).

number

[0019] According to the third alternative form, the method is: 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, Step (J) determines at least one mean wind speed value over time interval T by projecting at least two components of the instantaneous wind speed vector reconstructed in step I over time interval T. Includes averaged projections.

[0020] Preferably, the hybridization method by averaged projection using a third alternative form is: In step I, based on equations (34) to (40), the three components of the instantaneous wind speed vector (U i , V i , W i ) at least two of the components (U i , V i ) or (V i , W i ) or (U i , W i A step in the vector reconstruction of ), where i is an integer between 1 and N that corresponds to the number of continuous projections of instantaneous wind speed vectors over a time interval (T) called the reference time interval, U i This is a component of the instantaneous wind velocity vector in the spatial direction (d1) extending within the spatial plane (p1), and component Vi This is the component of the instantaneous wind velocity vector in the spatial direction (d2) extending within the spatial plane p1, and component W i This is the component of the average wind speed vector in the spatial direction (d3) perpendicular to the plane p1,

number

number

number

number

number

number

number

[0021] Furthermore, a third alternative form of hybridization using temporal combinations is: In step J, based on equations (41) to (42), at least one average wind speed value (Vh) in planes p1, p2, or p3 over a reference time interval T, starting from at least two reconstructed components of the instantaneous wind speed vector. ave ) is a decision-making step,

number

number

number

[0022] The hybridization method by averaged projection using a third alternative form may include the estimation of wind direction (dir) in the plane p1 according to equation (44).

number

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

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

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

[0026] According to the present invention, a computer program is also proposed that, when the program is executed by a computer, includes instructions that cause the computer to perform a method in any one of the first, second, and / or third alternative forms.

[0027] According to the present invention, When executed by a computer, it includes instructions that cause the computer to perform one of the first, second, and / or third alternative forms of the method, and / or The computer program according to the present invention is recorded. Recording media are also proposed. [Brief explanation of the drawing]

[0028] Other advantages and features of the present invention will become apparent from the detailed description of embodiments and embodiments, which are not limiting, as well as from the accompanying drawings below. [Figure 1] This is a schematic diagram of a perspective view of an optical system for acquiring telemetry data used to carry out the method according to the present invention. [Figure 2]This is a functional diagram of a first alternative embodiment of the method according to the present invention. [Figure 3] This is a functional diagram of a second alternative embodiment of the method according to the present invention. [Figure 4] This is a functional diagram of a third alternative embodiment of the method according to the present invention. [Modes for carrying out the invention]

[0029] The embodiments described below are not limiting in any way, and in particular, if the selection of this characteristic is sufficient to provide a technical advantage or to distinguish the invention with respect to the state of the prior art, it is possible to consider that modifications of the invention include only the selection of the characteristic described, apart from the other characteristics described (even if this selection is separated within a phrase containing these other characteristics). If this selection is sufficient to provide a technical advantage or to distinguish the invention with respect to the state of the prior art, then the selection includes at least one, preferably functional, characteristic having no structural detail or only a portion of a structural detail.

[0030] Figure 1 shows an example of an optical system 1 for acquiring telemetry data. According to the example of Figure 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 Figure 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 fact, wind is defined by a wind vector that includes 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 takes five measurement values S Ni 、S Si 、S Ei 、S Wi, and S Vi It delivers the data. Therefore, measurements are available approximately every 0.8 seconds. Thus, over a time interval Ω, there are M projected components S equal to the value obtained by dividing Ω (in seconds) by 4. Ni S Si S Ei S Wi , and S Vi A set exists. Next, the measured instantaneous projection S Ni S Si S Ei S Wi , and S Vi It is necessary to reconstruct the components of the wind vector (U, V, W). In practice, the method for processing telemetry data to estimate wind speed according to the present invention can be applied to data measured in real time, or to 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 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 to estimate wind speed is performed on data relating to instantaneous wind speed vectors measured by telemetry, for example by LIDAR.

[0033] A typical use of wind measurements is to measure the power available to generate energy for a wind turbine, in which case the measurement interval T is typically a 10-minute time interval, allowing for the isolation of 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 allows for the estimation of the average wind speed over this reference interval. According to a particular 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, equations 1 and 2, as well as the two components U of the wind vector along the north / south axis and east / west axis, respectively Ω , V Ω Based on this, a step of reconstructing the average wind speed vector over a divided time interval Ω, wherein the component U Ω The component of the average wind speed vector is the spatial direction d1 corresponding to the north / south axis extending within the spatial plane p1 corresponding to the plane tangential to the ground surface at the level of the measurement point, and component V Ω This is the component of the average wind speed vector in the spatial direction d2 corresponding to the east / west axis extending within the spatial plane p1,

number

number

number

[0034] In practice, the value of the division time interval Ω is constant or modified during the acquisition of telemetry data, and the said value of the division time interval Ω is The type of telemetry system from which the telemetry data was acquired, and / or Atmospheric conditions during the acquisition of the aforementioned telemetry data It is a function of .

[0035] The value of Ω can be adapted to the amplitude of the variation in the direction and speed of the horizontal wind, which can be shown, for example, by calculating the standard deviation of the direction and horizontal wind speed, or to the value of the estimated mean wind speed. Figure 3 shows a schematic diagram of the function of a method for processing telemetry data and estimating wind speed according to a second alternative embodiment of the present invention. According to the second alternative embodiment, the method is: 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, Step D of vector reconstruction over a time interval called the reference time interval T, starting with N of two reconstructed components of the instantaneous wind speed vector included over a reference time interval T (a 10-minute time interval is 600 seconds, and obtaining a set of projections of the instantaneous velocity vector every 4 seconds is equal to 150), Step E of scalar reconstruction of at least one average wind speed value over a time interval T, starting from at least two reconstructed components of the average wind speed vector, Step F of determining at least one value of the norm of instantaneous wind speed, starting from the projection of the instantaneous wind speed vector, Step G of determining at least one mean value of the wind speed norm over a reference time interval T, starting from at least one value of the reconstructed instantaneous wind speed norm, The weighting of the sum of at least one mean wind speed value reconstructed in step E and at least one mean wind speed value determined in step G gives at least one mean wind speed Vh over time interval T. ave The decision steps (H) and This includes weighted hybridization.

[0036] According to a specific embodiment of a second alternative form of the method according to the present invention, the weighted hybridization method is: In step C, based on equations 11 and 12 respectively, at least two components (U) of the instantaneous wind speed vector along the north / south axis and the east / west axis are obtained. i , V i ) is a reconstruction step where i is an integer between 1 and N that corresponds to the number of continuous projections of instantaneous wind speed vectors over a time interval T called the reference time interval, U i This is the 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 plane tangential to the ground surface at the level of the measurement point, and component V i This is the component of the instantaneous wind velocity vector in the spatial direction d2, which corresponds to the east / west axis extending within the spatial plane p1.

number

number

number

number

number

number

number

number

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

[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, and more preferably equal to 0.33. Under standard atmospheric conditions, for the telemetry system configuration shown in Figure 1, the number that allows for obtaining the best estimate is approximately 0.33.

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

number

[0040] Figure 4 shows a schematic diagram of the function of a method for processing telemetry data and estimating wind speed according to a third alternative embodiment of the present invention. According to the third alternative embodiment, the method is: 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, By projecting at least two components of the instantaneous wind speed vector reconstructed in step I over a time interval T, at least one mean value Vh of the wind speed over the time interval T is obtained. ave Step J of the decision This includes hybridization by averaging projection. According to a specific embodiment of a third alternative form of the method according to the present invention, the hybridization method by averaging projection is: In step I, based on equations 34 and 35, the two components of the instantaneous wind speed vector along the north / south axis and the east / west axis are calculated (U i , V i A step in the vector reconstruction of ), where i is an integer between 1 and N, corresponding to the number of continuous projections of instantaneous wind speed vectors over a time interval T called the reference time interval, U iThis is the 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 plane of the ground surface at the level of the measurement point, and component V i This is the component of the instantaneous wind velocity vector in the spatial direction d2, which corresponds to the east / west axis extending within the spatial plane p1.

number

number

number

[0041] Of course, the present invention is not limited to the examples described herein, and many modifications can be made to these examples without departing from the scope of the invention. Therefore, modifications or combinations of the embodiments described above are conceivable.

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

Claims

1. 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 the mean wind speed vector over a time interval (Ω) called the split time interval, starting from successive projections over time of the instantaneous wind speed vector; Starting from the number (T / Ω) of the at least two components of the mean wind speed vector reconstructed in step A, at least one mean wind speed value (Vh ａｖｅ ) and a step (B) of scalar reconstruction of Hybridization by time combination, comprising: 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 divided time interval Ω included in the reference time interval T. method.

2. In step A, the three components (U Ω , V Ω , W Ω The at least two components (U Ω , V Ω ) or (V Ω , W Ω ) or (U Ω , W Ω ) a step of reconstructing said component U Ω is the component of the mean wind speed vector in a spatial direction (d1) extending in a spatial plane (p1), and the component V Ω is the component of the average wind speed vector in a spatial direction (d2) extending in the spatial plane p1, and the component W Ω is the component of the average wind speed vector in a spatial direction (d3) perpendicular to the plane p1, [0010] [0025] [0030] [0045] [0050] [006] [0070] where i is the successive projection S of the instantaneous wind speed vector over the divided time interval Ω Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S ｖｉ is an integer between 1 and M corresponding to S Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S ｖｉ are the respective 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 merged with the direction d3, θ is a non-zero angle formed between the axis a1 and the normal to the plane p1 and between the axis a2 and the normal to the plane p1, γ is a non-zero angle formed between the axis a3 and the normal to the plane p1 and between the axis a4 and the normal to the plane p1, the first and second axes a1 and a2 are contained in a plane (p2), the third and fourth axes a3 and a4 are contained in a plane (p3), and the planes p2 and p3 form a non-zero angle α therebetween. Steps and In step B, based on equations (8) to (10) and the at least two components of the average wind speed vector reconstructed in step A, at least one average wind speed value (Vh ａｖｅ ) a step of scalar reconstruction of [0080] [0090] [0010] Q is an integer between 1 and (T / Ω) corresponding to the number of the at least two components of the average wind speed over the divided time interval Ω included in the reference time interval T. Steps and The method of claim 1 , comprising:

3. The value of the division time interval Ω is constant or modified during acquisition of the telemetry data, and the value of the division time interval Ω is the type of telemetry system from which the telemetry data was obtained; and / or Atmospheric conditions during acquisition of said telemetry data The method of claim 1 or 2, wherein the function is

4. 1. 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 the instantaneous wind speed vector, starting from a projection of said instantaneous wind speed vector; A step (D) of vector reconstruction of at least two components of the mean wind speed vector over a time interval (T) called the reference time interval, starting from the number N of the at least two components of the instantaneous wind speed vector reconstructed in step C included over the 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 mean 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; (G) 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 reconstructed in step F and the at least one average wind speed value determined in step G, at least one average wind speed value (Vh ａｖｅ (H) determining A method comprising weighted hybridization comprising:

5. In step C, the three components of the instantaneous wind speed vector (U ｉ , V ｉ , W ｉ At least two components (U ｉ , V ｉ ) or (V ｉ , W ｉ ) or (U ｉ , W ｉ ) a vector reconstruction step, where 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, and U ｉ is the component of the instantaneous wind speed vector in a spatial direction (d1) extending in a spatial plane (p1), and the component V ｉ is the component of the instantaneous wind speed vector in a spatial direction (d2) extending in the spatial plane p1, and the component W ｉ is the component of the average wind speed vector in a spatial direction (d3) perpendicular to the plane p1, ##EQU00011## ##EQU00012## ##EQU00013## ##EQU14## ##EQU00015## ##EQU00016## ##EQU00017## In the formula, S Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S Ｖｉ are the 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 merged with the direction d3, θ is a non-zero angle formed between the axis a1 and the normal to the plane p1 and between the axis a2 and the normal to the plane p1, γ is a non-zero angle formed between the axis a3 and the normal to the plane p1 and between the axis a4 and the normal to the plane p1, the first and second axes a1 and a2 are contained in a plane (p2), the third and fourth axes a3 and a4 are contained in a plane (p3), and the planes p2 and p3 form a non-zero angle α between them. Steps and In step D, at least two components (U) of the average wind speed vector over the reference time interval T are calculated based on equations (18) to (20). Ｎ , Vvect Ｎ ) or (Vvect Ｎ , Wvect Ｎ ) or (Uvect Ｎ , Wvect Ｎ ) vector reconstruction step, Ｎ is the component of the wind speed in the spatial direction d1, and the component V Ｎ is the component of the wind speed in the spatial direction d2, and the component W Ｎ is the component of the wind speed in the spatial direction d3, [0018] [0019] [0020] and In step E, at least one value of the instantaneous wind speed (Vscal ｉ ) a scalar reconstruction step of Vscal ｉ corresponds to the time series of instantaneous wind speed values ​​in plane p1, p2 or p3, respectively; ##EQU00021## [0022] [0023] and In step F, the instantaneous wind speed value Vscal reconstructed in step E is calculated based on equations (24) to (26). ｉ．１ Or Vscal ｉ．２ Or Vscal ｉ．３ At least one average wind speed value (Vhscal) in the planes p1, p2, or p3, respectively, over said reference time interval T, starting from ａｖｅ ) determining the ##EQU00024## [0025] [0026] and In step G, starting from the at least two reconstructed components of the mean wind speed vector, at least one mean wind speed value (Vvect ａｖｅ ) determining the [0027] [0028] [0029] and In step H, the reconstructed wind speed value (Vhscal ａｖｅ．１ , Vhvect ａｖｅ．１ ) or (Vhscal ａｖｅ．２ and Vhvect ａｖｅ．２ ) or (Vhscal ａｖｅ．３ , Vhvect ａｖｅ．３ ) in the planes p1, p2, or p3, respectively, over said reference time interval T, ａｖｅ ) [0030] [0031] [0032] where P is a dimensionless weighting factor between 0 and 1 Steps and The method of claim 4 , comprising:

6. The method of 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 modified during acquisition of the telemetry data or when performing the method, and the value of the division time interval Ω is the type of telemetry system from which the telemetry data was obtained; and / or Atmospheric conditions during acquisition of said telemetry data The method according to any one of claims 5 to 6, wherein the value is a function of

8. Estimation of the variation σ of the wind speed over the reference time interval T according to equation (33), [Equation 33] where c is a positive number and σ is 0 or a positive dimensionless number. The method according to any one of claims 4 to 7.

9. 1. A method for processing telemetry data to estimate wind speed, the method comprising: The average projection includes: The three components of the instantaneous wind speed vector (U ｉ , V ｉ , W ｉ At least two components (U ｉ , V ｉ ) or (V ｉ , W ｉ ) or (U ｉ , W ｉ ) vector reconstruction step I, where i is an integer between 1 and N corresponding to the number of successive projections of the instantaneous wind speed vector over a time interval (T) called the reference time interval, and U ｉ is the component of the instantaneous wind speed vector in a spatial direction (d1) extending in a spatial plane (p1), and the component V ｉ is the component of the instantaneous wind speed vector in a spatial direction (d2) extending in the spatial plane p1, and the component W ｉ is the component of the average wind speed vector in a spatial direction (d3) perpendicular to the plane p1, [0034] [Equation 35] [0036] [Equation 37] [Equation 38] [0039] [0040] In the formula, S Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S Ｖｉ are the 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 merged with the direction d3, θ is a non-zero angle formed between the axis a1 and the normal to the plane p1 and between the axis a2 and the normal to the plane p1, γ is a non-zero angle formed between the axis a3 and the normal to the plane p1 and between the axis a4 and the normal to the plane p1, the first and second axes a1 and a2 are contained in a plane (p2), the third and fourth axes a3 and a4 are contained in a plane (p3), and the planes p2 and p3 form a non-zero angle α between them. Step I; At least one average wind speed value (Vh) in the plane p1, p2, or p3, respectively, over the reference time interval T, is calculated by projecting the at least two components of the instantaneous wind speed vector reconstructed in step I over the time interval T based on equations (41) to (43). ａｖｅ A step J of determining [0041] [0042] [0043] Step J, A method comprising:

10. including the estimation of the wind direction (dir) at plane p1 according to equation (44), [0044] In the formula, tan －１ is an arctangent function, the wind direction Dir estimated from Equation 44 is an angle value between the wind direction and the direction d1, and Vrec and Urec are respectively a scalar value of the component of the wind speed in plane p1 over the reference time interval T; or Average vector velocity of the components of the wind speed in plane p1 over the reference time interval T That is, 10. The method of any one of claims 2, 5 and 9.

11. The projection S of the instantaneous wind speed vector is measured by at least one measurement laser beam extending along a first axis a1, a second axis a2, a third axis a3, a fourth axis a4, and a fifth axis a5, respectively. Ｎｉ , S Ｓｉ , S Ｅｉ , S Ｗｉ , and S Ｖｉ The method of any one of claims 2, 5 and 9, comprising measuring:

12. A method according to any one of claims 1 to 11 which is computer-implemented.

13. A data processing device comprising means arranged and / or programmed and / or configured to carry out 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 carry out a method according to any one of claims 1 to 12.

15. and / or comprising instructions which, when executed by a computer, cause the method according to any one of claims 1 to 12 to be carried out; A computer program according to claim 14, Recording medium.

Images