Issue |
Renew. Energy Environ. Sustain.
Volume 1, 2016
|
|
---|---|---|
Article Number | 22 | |
Number of page(s) | 4 | |
DOI | https://doi.org/10.1051/rees/2016025 | |
Published online | 13 June 2016 |
Research Article
Comparison of different methods for evaluation of wind turbine power production based on wind measurements
1 Institute of Physical Energetics, 21 Aizkraukles Str., Riga LV-1006, Latvia
2 Ventspils University College, 101 Inzenieru Str., Ventspils LV-3601, Latvia
⁎ e-mail: elmag@inbox.lv
Investigations of the wind shear up to the height of 200 (m) on the Latvian coast of the Baltic Sea have been carried out using a Pentalum SpiDAR laser measuring complex. Based on wind speeds measurements for three levels – 30, 40 and 50 (m), assessment of the operational efficiency of the wind turbines for heights 100, 140 and 180 (m) have been performed. For comparison, this analysis involves five different approaches: the Rayleigh frequency distribution, three different Weibull frequency distributions and method based on approximation of the cubic wind speed. Results are compared with measurements on the corresponding heights.
© V. Bezrukovs et al., published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The methods of analysis are based on extrapolation of the results obtained by measurements of the wind speed and standard deviation (SD), and at modelling the probability density curve of wind speed function F(V) using Rayleigh's models and Weibull's parameters k (shape factor) and c (scale factor) [1].
As an example, it is taken the operating characteristics of a 2300 (kW) Enercon E-82 wind turbines (WTs) [2] and have shown to what extent the precision of the calculation results obtained for the WT operational efficiency at different heights depends on the choice of the model. The proper choice of the model would enable calculation of the annual electricity production (AEP) by generators, making it possible to estimate the efficiency of WT offered for the use of a particular WT type.
2 The results of wind shear investigation
Figure 1 shows the map of disposition of a laser measuring complex Pentalum SpiDAR in the territory of Latvia on the Baltic Sea coast [3]. The measuring complex is installed in the territory of the Ventspils International Radio-Astronomy Centre at the height of 13.2 (m) above ground level and at a distance of 2.5 (km) from the seashore on a platform with coordinates N 57°33′12.00″ and E 21°51′16.00″.
The Pentalum SpiDAR performs 5 (s) measurements of wind speed and direction for ten height levels: 30, 40, 50, 80, 100, 120, 140, 160, 180 and 200 (m). Apart from that, the records are made for pressure, humidity, and temperature of the air. Web system is developed for retrieving data from SpiDAR ftp server.
Figure 2 displays the curves of average monthly wind speed values, Vavg, which are calculated by averaging the bins, i.e. the measured wind speed values V accumulated within 10 min in the time period T from 01.03.2014 to 01.03.2015 using the laser measuring complex Pentalum SpiDAR at the heights 30, 50, 100, 140 and 180 (m) above ground.
Figure 3 demonstrates the distribution curves of the average wind speed Vavg and the average cubic wind speed Vavg cub (m/s) values along with the SD value for the wind shear up to the height of 200 (m). These values have been calculated based on analysis of the data accumulated during 12 months. The distribution of the obtained wind shear values is well approximated by the expression:(1)where
-
hr is the height of wind speed measurement (m),
-
h is the height of WT location (m),
-
α is the approximation coefficient.
Using the force function for approximation of the measurement results, expression (1) will assume the form [1]:(2)
The distribution curve of SD values is approximated using the exponential function which for the measurement results obtained will be:(3)
To approximate the values of cubic wind speed the force function described by the equation [2] can be used:(4)
For adequate comparison of the modelling methods and their further analysis, at determination of expressions for approximating functions (2), (3) and (4), the values corresponding to the three heights (30, 40 and 50 (m)) are used. In these cases, the regularities describing the distributions of wind speed Vavg and SD values give a good coincidence with the measurement results up to the height of 200 (m). However, for the cubic wind speed Vavg cub the divergence of the approximating curve calculated using expression (4) from the real values increases with height.
![]() |
Fig. 1 Laser measuring complex Pentalum SpiDAR installed on the Baltic Sea coast at a height of 13.2 (m) above ground level and at a distance of 2.5 (km) from the sea coast. |
![]() |
Fig. 2 Distribution curves of the average wind speed values Vavg (m/s) obtained using the laser measuring complex Pentalum SpiDAR for the heights 30, 50, 100, 140 and 180 (m) above ground in the measurement time T. |
![]() |
Fig. 3 Approximation curves for the measured values of wind speed Vavg, Vavg cub (m/s), and for the standard deviation, SD, depending on height h by three points for heights of 30, 40 and 50 (m) above ground. |
3 Comparison of methods for calculation AEP of WT
The measuring complex Pentalum SpiDAR makes it possible to obtain a realistic picture for wind shear up to 200 (m); therefore, estimation of the model-related discrepancies in the calculated values of AEP by WTs located at different heights can be performed.
The algorithms of calculating the forecasted AEP for the heights of 80, 100, 140, 180 and 200 (m) could systemized as follows:
- –
I (Weibull) – AEP determination by modelling the Weibull frequency distribution curves based on average wind speeds Vavg and SD for corresponding heights.
- –
II (Vavg cub) – AEP determination by the value of average cubic wind speed Vavg cub calculated using approximating function (4).
- –
III (Rayleigh) – AEP determination using approximating function (2) and calculating the Rayleigh wind speed frequency distribution.
- –
IV (Weibull) – AEP determination using approximating functions (2) and (3) in order to find k- and c-factor values.
- –
V (Weibull) – AEP determination by approximating k- and c-factor values of Weibull frequency distribution function based on calculated k and c values at heights 30, 40 and 50 (m).
The modelling of annual energy production using models I ÷ V is performed taking the power curve P(V) (kW) of a 2300 kW wind turbine of the E82 type. Figure 4 shows the WT power curve, Weibull and Rayleigh curves of wind speed frequency distributions for 100 (m) height. The value of average cubic wind speed for this height Vavg cub = 7.57 (m/s) and for this cubic wind speed corresponding power of generator is 694 (kW).
Analysis of curves in Figure 4 shows that the methods for AEP calculation considered above are not mutually connected by measurement parameters. This makes it possible to compare the modelling results with each other as well as with the calculation results obtained from measurements.
The total electricity production W (GWh) to which the wind stream energy can be converted in a definite period of time is determined by the equality:(5)where the values of function P(Vi) are corresponding to the power curve for the generator, while those of function F(Vi) – to the Weibull's wind speed frequency distribution curve at the height of rotor hub and respective wind speeds Vi (m/s).
Table 1 shows the results for forecast AEP (GWh) calculated using the data obtained during wind speed measurements as well as the results for forecast AEP according to variants I ÷ V at heights of 80, 100, 140, 180 and 200 (m) above ground for the WT type E82. To determine the relative uncertainty of the results obtained at AEP calculation for different models the below shown expression is used:(6)where Wmeasur is the forecasted AEP value (GW) calculated from the wind speed measurements, Wmod is the forecasted AEP value (GW) calculated from the results of modelling (models I ÷ V).
The results given in Table 1 are presented graphically in Figure 5, where it is seen that the relative deviation of AEP calculation WΔ is decreasing with the height h increasing in the Rayleigh model while in the Weibull model this value depends but slightly on the height changes. As seen the results of WΔ calculation for first and second models vary with height from +0.4 to −17.6 (%).
At the same time, the relative error of the AEP calculation based on the average cubic wind speed is varying in the limits from −3.1 to +6.2 (%). This peculiarity of the models’ behaviour can be employed for increasing the precision of AEP calculation due to summing up and averaging predicted values.
![]() |
Fig. 4 Curves of wind speed frequency distributions F(V) (%), by Weibull and Rayleigh, the power curve P(V) (kW) of WT type E82 for the height of 100 (m) and the corresponding value of average cubic wind speed Vavg cub. |
Results of AEP calculation for a WT type E82 obtained using models: I (Weibull), II (Vavg cub), III (Rayleigh), IV (Weibull), V (Weibull) and based on the wind speed measurements at heights 80, 100, 140, 180 and 200 (m).
![]() |
Fig. 5 Relative deflection WΔ (%) of AEP calculation based on models I ÷ V for the heights 80, 100, 140, 180 and 200 (m). |
4 Conclusions
Comparison of AEP calculation methods based on wind speed measurements on 30, 40 and 60 (m) height and wind shear extrapolation up to 200 (m) height shows that application of Rayleigh and Weibull models gives a result with relative deflection up to −17 (%). The method of application of the value of average cubic wind speed gives relative deflection up to 6.2 (%).
Acknowledgments
This paper was supported by European Social Fund project “Assessment of wind energy potential in Latvia and environmental impact from wind energy installations”, No. 2014/0010/1DP/1.1.1.2.0/13/APIA/VIAA/033.
References
- J. Manwell, J. McGowan, A. Rogers, Wind energy explained: theory, design and application ( John Wiley & Sons Ltd., 2009), 689 [Google Scholar]
- Enercon wind energy converters: product overview 2010, Technical information. Version: July 2010: www.enercon.de/p/downloads/EN_Productoverview_0710.pdf [Google Scholar]
- V. Bezrukovs, V. Bezrukovs, A. Zacepins, Comparative efficiency of wind turbines with different heights of rotor hubs: performance evaluation for Latvia, J. Phys.: Conf. Ser. 524012113 (2014), 9 p. DOI: 10.1088/1742-6596/524/1/012113 [Google Scholar]
Cite this article as: Valerijs Bezrukovs, Aleksejs Zacepins, Vladislavs Bezrukovs, Vitalijs Komasilovs, Comparison of different methods for evaluation of wind turbine power production based on wind measurements, Renew. Energy Environ. Sustain. 1, 22 (2016)
All Tables
Results of AEP calculation for a WT type E82 obtained using models: I (Weibull), II (Vavg cub), III (Rayleigh), IV (Weibull), V (Weibull) and based on the wind speed measurements at heights 80, 100, 140, 180 and 200 (m).
All Figures
![]() |
Fig. 1 Laser measuring complex Pentalum SpiDAR installed on the Baltic Sea coast at a height of 13.2 (m) above ground level and at a distance of 2.5 (km) from the sea coast. |
In the text |
![]() |
Fig. 2 Distribution curves of the average wind speed values Vavg (m/s) obtained using the laser measuring complex Pentalum SpiDAR for the heights 30, 50, 100, 140 and 180 (m) above ground in the measurement time T. |
In the text |
![]() |
Fig. 3 Approximation curves for the measured values of wind speed Vavg, Vavg cub (m/s), and for the standard deviation, SD, depending on height h by three points for heights of 30, 40 and 50 (m) above ground. |
In the text |
![]() |
Fig. 4 Curves of wind speed frequency distributions F(V) (%), by Weibull and Rayleigh, the power curve P(V) (kW) of WT type E82 for the height of 100 (m) and the corresponding value of average cubic wind speed Vavg cub. |
In the text |
![]() |
Fig. 5 Relative deflection WΔ (%) of AEP calculation based on models I ÷ V for the heights 80, 100, 140, 180 and 200 (m). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.