18-Power F (Speed)
18-Power F (Speed)
sciences
Article
Variability of the Wind Turbine Power Curve
Mahesh M. Bandi 1,∗ and Jay Apt 2
1 Collective Interactions Unit, Okinawa Institute of Science and Technology, Onna, Okinawa 9040495, Japan
2 Department of Engineering & Public Policy and Tepper School of Business, Carnegie Mellon University,
5000 Forbes Avenue, Pittsburgh, PA 15213, USA; apt@cmu.edu
* Correspondence: bandi@oist.jp; Tel.: +81-98-966-1546
1. Introduction
The wind turbine power curve relates the speed of wind blowing past a turbine to the power
generated by the turbine. Wind plant operators forecast the power they expect to generate by feeding
wind speed forecasts from numerical weather models to these power curves [1,2]. The power curves
are supplied to operators by turbine manufacturers, who calibrate them under standards specified
by the International Electrotechnical Commission (IEC) [3]. The IEC standard considers the average
behaviour between the mean wind speed v (x being the time-average of a time-varying quantity x (t))
and mean power output P, and hence does not locally hold in time for instantaneous measurements.
Indeed, instantaneous values of wind speed v(t) and wind power P(t) exhibit significant scatter about
the mean profile (Figure 1). Several studies [4–8] have focused on the factors contributing to this
variability, including turbulent fluctuations, wind shear, directional shear, directional fluctuations, etc.,
with the aim of accurately modeling the “mean” profile of the wind turbine power curve [9–12].
In this article, we first establish that the wind turbine power curve has features that are sensitive
to local environmental factors. We then employ a combination of theory and empirical results to
self-consistently account for this dependence on local factors. Our proposed approach was especially
developed within the framework of IEC standard 614001-12-1 to ensure its easy adoption by wind
plant operators. In particular, we present a minimal parameter
q description of wind speed variability
q
σv = (v(t) − v) and wind power variability σP = ( P(t) − P)2 in terms of the mean wind speed
2
v. Our objective here is to remain faithful to the IEC standard [3] that prescribes the power curve with
mean wind speed v as its only parameter. We present empirical evidence that the standard deviation
in wind speed σv systematically varies with mean wind speed v. At least in some instances, this
monotonic variation follows algebraic scaling of the form σv = C × vα , where C is a constant and α is
a fractional power. This scaling form—which we attribute to residual signal correlations that remain
post-averaging—then affords a description of mean wind power P and its standard deviation σP in
terms of v alone.
Figure 1. Instantaneous power P(t) versus instantaneous wind speed v(t) (solid grey circles) and
time-averaged power P versus time-averaged wind speed v (solid red circles) for the Howard data
set (Table 1). Considerable scatter in P(t) versus v(t) occurs about the time-averaged power curve.
The scatter increases with mean speed v, as qualitatively shown with blue arrows at v = 5, 7, and 9 m/s.
Our analysis of wind data obtained from three different planetary locations (Table 1) reveals that
in instances when the above scaling form is satisfied, the power-law exponent α varies with geographic
location, and hence must reflect local environmental factors not captured by manufacturer-supplied
calibration power curves. Consequently, when these calibration power curves are applied to forecast
wind power, wind speed variability transforms into power forecast uncertainty. Since the variability
always increases in tandem with mean speed, the resulting forecast error is multiplicative, so it
substantially increases uncertainty in wind power forecasts. We conclude with a proposal that wind
plant operators should recalibrate turbine power curves at the plant location to properly account for
variability arising from local environmental factors. This will help to reduce the uncertainty of wind
power forecasts.
At first sight, the monotonic increase in wind speed variability σv with mean speed v seems
at odds with the general belief that the turbulence intensity I = σvv decreases with speed. In many
instances, one observes that fluctuations σv exhibit non-monotonic behaviour with steady increase
in mean speed v (and Reynolds number Re) [13]. In particular, frictional losses with the walls cause
loss of pressure head in water flow down a river or pipe. This in turn causes loss in mean speed
Appl. Sci. 2016, 6, 262 3 of 9
(v), but fluid shear with the wall can increase fluctuations and eventual transition to turbulence [14].
Atmospheric flows close to the Earth’s surface share close correspondence with pipe and river flows,
in that the Earth’s surface behaves like a rough wall, but the flow is unbounded from above, where
a second confining wall is absent. Fluctuations σv (and therefore turbulence intensity I) will vary
non-monotonically with mean speed v close to Earth’s surface. However, within the range of wind
speeds between the cut-in and cut-out speed of a turbine, this variation should exhibit monotonic
increase in σv with a concomitant rise in v.
Whereas the IEC standard considers only mean quantities, as we show below, both mean power
output (Equation (5)) and its standard deviation (Equation (7)) strongly depend upon wind speed
variability σv , in addition to mean wind speed v. Strong local environmental dependence of wind
speed fluctuations naturally affects both the mean profile of the power curve and its variability. When
not properly accounted for, this increases forecast uncertainty, which in turn adds costs to renewable
energy production [15,16]. Understanding the source of variability and utilising it appropriately
therefore brings tangible benefits to the global renewable energy community.
1
p(t) = ρAv(t)3 (1)
2
Here v(t) is the time-varying streamwise component of velocity blowing past the turbine. We
ignore directional fluctuations merely in the interest of keeping our arguments accessible (as also
recently explained by Hedevang [8]), but the subsequent analysis can be easily extended to include
directional fluctuations [4].
Performing Reynolds decomposition on the time-varying velocity v(t) = v + ṽ(t) into a
time-independent mean speed v and a time-varying fluctuation ṽ(t) whose long-time average ṽ(t) ≡ 0
by definition, and substituting it in Equation (1), we obtain:
1 1
p(t) = ρA[v + ṽ(t)]3 = ρA[v3 + 3v2 ṽ(t) + 3vṽ(t)2 + ṽ(t)3 ] (2)
2 2
Owing to the turbine experiencing a drag force, it does not extract all the energy available in the
wind. Indeed, the power P(t) generated by the turbine is given by:
Term I in RHS of Equation (4) represents the turbine’s mean power output under steady
(time-independent) wind speed conditions. Term II is zero, since ṽ(t) ≡ 0 by definition. ṽ(t)3
can take both positive and negative values, owing to its being an odd function of time t. Consequently,
Term IV on the RHS of Equation (4) takes a very small but non-zero value, and will be neglected
in the foregoing analysis. Term III captures wind speed variability, and ṽ(t)2 , being quadratic (an
even function in time t), is always positive. Indeed, ṽ(t)2 ≡ (v(t) − v)2 ≡ σ2v , the quantity of specific
interest to our present analysis. As a point of comparison, a 20 day long-time average for the Howard
dataset (Table 1) yielded the values v3 = 299.85 m3 /s3 for term I, 3vṽ(t)2 = 192.26 m3 /s3 for term
III, and v(˜t)3 = 23.06 m3 /s3 for term IV. We are therefore justified in keeping term III (64% of term I)
and dropping term IV (only 8% of term I). Dropping terms II and IV and re-arranging terms I and III,
Equation (4) can be re-expressed as:
h i h i
P = K v3 + 3vṽ(t)2 = Kv3 1 + 3I 2 (5)
where I = σvv is the turbulent intensity, or alternatively the coefficient of variation in wind speed.
The presence of a second fluctuation term in Equation (5) requires a comparison of its strength
relative to the mean term (v3 ). When the turbulent intensity is low (I 1), wind conditions
are close to steady, and the mean power follows the standard P ∼ v3 relation. If the fluctuation
magnitude is a significant fraction of mean speed (I . 1), then 3I 2 ∼ 1. Two scenarios must then be
considered: one where the standard deviation of wind speed σv remains constant or decreases with
increasing mean wind speed v, and the second where σv increases with v. A consideration of the two
scenarios becomes important, because a constant σv leads to additive variability, whereas σv increasing
with v leads to multiplicative variability, which in turn has implications for the error. Our analysis of
wind data presented below always shows an increase in σv with v, hence we discuss only the case of
multiplicative variability. q
The standard deviation of power σP ≡ ( P(t) − P)2 can also be computed as a function of
v and σv . Defining instantaneous power as P(t) = P + P̃(t), where P̃(t) is the fluctuation about
the time-independent mean P, whose long time average is zero by definition ( P̃(t) ≡ 0), and using
Equations (3) and (4):
An algebraic expansion of Equation (6) followed by a long-time average of individual terms yields:
Keeping only the leading order term and ignoring all higher order terms, we arrive at a first order
approximation for the standard deviation of wind power:
q q
σP = P̃(t)2 = 3Kv2 ṽ(t)2 = 3Kv2 σv (7)
From Equation (7), we see that σP scales linearly with σv , but is amplified by a factor of 3Kv2 .
More importantly, we see that the magnitude of wind power fluctuations scales quadratically with
mean wind speed (σP ∼ v2 ). We note that the functional form for σP in Equation (7) is expected to
apply between the cut-in and rated speeds.
In Figure 1, we plot the power curve (P vs. v, solid red circles) with the instantaneous power
against speed (P(t) vs. v(t), solid grey circles) overlaid on top of the power curve for a 2.05 MW
REPower MM92 turbine (cut-in speed: 3 m/s, rated speed: 11.2 m/s, cut-out speed: 24 m/s,
rated power: 2050 kW) located at a wind farm operated by EverPower Wind Holdings in Howard,
Appl. Sci. 2016, 6, 262 5 of 9
NY (Table 1). The solid blue arrows marking the scatter in instantaneous wind speed values at v = 5,
7, and 9 m/s qualitatively demonstrate the monotonic increase in wind speed variability with mean
speed; i.e., the case of multiplicative variability.
Before proceeding, we note a nuance concerning the time-average for P and v. IEC standard
61400-12-1 [3] requires that each value of v and P in the power curve be determined by averaging
a 10 min time record (sampled at 1 Hz or higher) of v(t) and P(t), respectively. This translates to
an average over 600 s × 1 Hz = 600 samples—at least—for each time-averaged value in the power
curve. Since the Howard data set (see Table 1) was sampled at 0.2 Hz, a pertinent question to ask is
whether 600 s × 0.2 Hz = 120 samples are sufficient to achieve statistical convergence of the mean.
We performed a bootstrap protocol [17] for the Howard time series to generate ten different time
series by randomly shuffling values of the originally measured time series for v(t) and p(t). We then
measured the mean for each bootstrap time series as a function of increasing sample size for wind
speed and wind power. Figure 2 shows how the mean wind speed (Figure 2a) and mean wind power
(Figure 2b) converge as a function of increasing sample size. Since each sample is measured 5 s
apart in our data set, 120 samples equate to a 10 min interval prescribed by the IEC standard. As is
evident from Figure 2, asymptotic statistical convergence of the mean values is achieved only at
about 400 samples. Following from this statistical convergence test, we constructed the power curve
in Figure 1 by averaging over 400 samples (i.e., a 33.3 min time record for each value of v and P).
We emphasise that the bootstrap protocol tests for statistical convergence in the absence of correlations.
The successive values in the measured time series are correlated; i.e., the nth value will lie within a
certain band relative to the (n − 1)th and (n + 1)th values, depending upon correlation strength and
correlation time. When this measured series is randomly shuffled to generate a synthetic series for the
bootstrap test, the correlations between successive values in the time series are lost. On the other hand,
when using the actual time series for v(t) and P(t) in constructing the power curve from averages,
such correlations persist and reveal themselves in the variability, as we discuss below.
a) b)
Figure 2. Statistical convergence of (a) mean wind speed; and (b) mean wind power versus number of
samples contributing to the average. Each trace represents a separate data set generated with a bootstrap
protocol by randomly shuffling values of the original Howard time series (Table 1). The averages
approach asymptotic convergence around 400 samples.
v and wind power P, respectively, for the Howard data set. The standard deviation of wind speed σv
(solid red circles in Figure 3a) scales algebraically relative to the mean wind speed v with a power-law
fit that follows the form σv = 0.25 × v0.73 (solid line in Figure 3a). The standard deviation of wind
power fluctuations σP (solid blue squares in Figure 3a) also exhibits algebraic scaling, albeit shallower
0.49
than wind speed, with a power-law fit of the form σP = 3.9 × P (dashed line in Figure 3a). Wind
speed standard deviation σv for the Big Bear Lake data set plotted in Figure 3b (red solid circles)
exhibits a similar algebraic scaling, but with a different power-law exponent (σv = 0.4 × v0.6 , solid
black line in Figure 3b) from the Howard data set. Finally, for the Atacama data set, although σv
exhibits a systematic increase with v (solid blue squares in Figure 3c), we find that the scatter is
too high to make a conclusive determination of power-law scaling. We do include a power-law fit
(dashed line in Figure 3c), but one should not derive any inferences from this fit. Be that as it
may, the fact that σv exhibits monotonic increase with v for all data sets establishes the scenario of
multiplicative variability.
Figure 3. (a) log–log scale: The standard deviation of wind speed fluctuations σv (red solid circles)
versus mean wind speed v and the standard deviation of wind power fluctuations σP (blue solid
squares) versus mean wind power P for the Howard data set. Power-law fits to the data provided
0.49
σv = 0.25 × v0.73 (solid line) and 3.9 × P (dashed line); (b) Whereas σv vs. v for Big Bear Lake
(red solid circles) exhibits power-law scaling with a fit value of 0.4 × v0.6 (solid line); and (c) the
Atacama data set (blue solid squares) reveals a monotonic increase in σv vs. v. Although we include a
power-law fit σv = 0.31 × v0.68 (dashed line) for illustrative purposes, the scatter in the data does not
permit one to place any confidence in the fit value of the exponent.
0.49
The observed scaling for power fluctuations (σP = 3.9 × P ) is significant given that
0.5
Calif et al. [18] have recently reported the same σP ∼ P scaling for a single turbine as well as
for wind farms at various planetary locations. They interpreted this scaling within the context of
the “Taylor Rule” [19] (not to be confused with Taylor’s Hypothesis in turbulence theory [20]), also
called “Fluctuation Scaling” [21] in physics, where the 1/2 scaling exponent forms one of two limiting
cases. However, based on additional empirical evidence for wind speed fluctuations and calculations
resulting in Equation (7), we proffer an alternative interpretation for this scaling.
Consider a time-varying signal x (t) with a signal correlation time τ0 . The mean x and standard
deviation σX become truly time-independent when the average is taken over several multiples of the
correlation time τ0 . Instead, if the averaging interval τ < τ0 , residual correlations persist in x τ and σxτ
(subscript τ now denotes the averaging interval) which vary with averaging time τ. Wind speed (and
therefore wind power) fluctuations which reflect atmospheric turbulence possess long time correlations
extending up to 24 h timescales [22–26], whereas the IEC standard [3] specifies a 10 min time average
for calculation of v and P. Even our bootstrap protocol specifically tests for statistical convergence in
the absence of temporal correlations; generating a new randomised time series from the original time
series destroys temporal correlations in the signal. Consequently, v and σv retain residual correlations,
Appl. Sci. 2016, 6, 262 7 of 9
which are revealed as systematic variation in σv relative to v (the same applies to σP and P). Such
relative scaling between moments of a distribution is well known via Extended Self Similarity (ESS)
scaling [27], and was fruitfully exploited to accurately estimate deviations from scalings predicted by
the Kolmogorov theory of turbulence [28]. Our observed scalings between the mean and standard
deviation of wind speed follow in the same spirit. One cannot escape the residual correlations—and
hence the power-law relationship between σv and v—unless the averaging time is increased from the
IEC-stipulated 10 min to the signal correlation time of order 24 h. A 24 h averaging time for each
point on the power curve calls for several months worth of data collection effort, and hence is clearly
impractical. A more tractable—and still quantitatively defensible—approach is to incorporate the
power-law scaling between σv and v to self-consistently account for residual correlations in the power
curve and its variability, as we show below in Equation (8).
We can now substitute the scaling form σv = C × vα in Equations (5) and (7) to obtain:
σ2
P = Kv 3
1 + 2v = Kv3 [1 + C2 v2(α−1) ] (8a)
v
σP = 3Kv2 σv = 3CKv(2+α) (8b)
planetary surface roughness. This hidden length scale—if identified through extensive observational
measurements—could potentially extend the phenomenology of flow within the shear boundary
layer beyond mean wind profiles to include fluctuations. Such a phenomenology would of course
mark an advancement in geophysical fluid dynamics, but would also benefit the wind engineering
community. High time-resolution wind data sets that extend over several days are few, and are prized
by researchers. The generation of such data sets in several locations would be a worthy goal for
national research agencies.
4. Summary
Before closing, we observe that collective interactions between turbines in a wind plant also
change the local flow field in the plant’s vicinity [30,31], in turn compromising plant efficiency [32,33].
Decomposing plant-induced variability from natural environmental factors is often complicated [30,34].
The methods we employ are agnostic to the source of variability, hence can account for both
natural (environmental) and artificial (plant-induced) sources of variability, while working within the
framework of IEC standard 61400-12-1 [3].
In summary, we have explained the functional dependence of P and σP in terms of v and σv .
From empirical evidence, we exploited systematic variation of σv with v to recast expressions for P and
σP solely in terms of v, thus keeping in accord with IEC standard 61400-12-1 [3]. Given the dependence
of the fluctuation scaling exponent α on local environmental factors, we propose that wind operators
should recalibrate their turbine power curves post-installation at the site to accurately account for the
site’s multiplicative variability. We believe that doing so, and including confidence intervals for wind
power variability about the mean power curve profiles, would reduce wind power forecast uncertainty.
Acknowledgments: Mahesh M. Bandi was supported by the Collective Interactions Unit at the Okinawa Institute
of Science and Technology Graduate University. Jay Apt received support from the Carnegie Mellon Electricity
Industry Center and from the Carnegie Mellon Climate and Energy Decision Making Center (CEDM) formed
through a cooperative agreement between the National Science Foundation and CMU (SES-0949710). We thank
Dhaval Balodia and EverPower Wind Holdings for access to the Howard wind and wind turbine power data.
Mahesh M. Bandi acknowledges helpful discussions with P. Chakraborty.
Author Contributions: M.M.B. and J.A. conceived and performed the research. M.M.B. wrote the manuscript
and J.A. edited it.
Conflicts of Interest: The authors declare no conflict of interest.
References
1. Landberg, L. Short-term prediction of the power production from wind farms. J. Wind Eng. Ind. Aerodyn.
1999, 80, 207–220.
2. Foley, A.M.; Leahy, P.G.; Marvuglia, A.; McKeogh, E.J. Current methods and advances in forecasting of wind
power generation. Renew. Energy 2012, 37, 1–8.
3. IEC. Wind Turbines–Part 12-1: Power Performance Measurements of Electricity Producing Wind Turbines;
IEC61400-12-1; International Electrotechnical Commission: Geneva, Switzerland, 2005.
4. Choukulkar, A.; Pichugina, Y.; Clack, C.T.M.; Calhoun, R.; Banta, R.; Brewer, A.; Hardesty, M. A new
formulation for rotor equivalent wind speed for wind resource assessment and wind power forecasting.
Wind Energy 2015, doi:10.1002/we.1929.
5. Wharton, S.; Lundquist, J.K. Atmospheric stability affects wind turbine power collection. Environ. Res. Lett.
2012, 7, doi:10.1088/1748-9326/7/1/014005.
6. Paiva, L.T.; Rodrigues, C.V.; Palma, J.M.L.M. Determining wind turbine power curves based on operating
conditions. Wind Energy 2014, 17, 1563, doi:10.1002/we.1651.
7. Sumner, J.; Masson, C. Influence of Atmospheric Stability on Wind Turbine Power Performance Curves.
J. Sol. Energy Eng. 2006, 128, 531–538.
8. Hedevang, E. Wind turbine power curves incorporating turbulence intensity. Wind Energy 2014, 17, 173,
doi:10.1002/we.1566.
Appl. Sci. 2016, 6, 262 9 of 9
9. Carrillo, C.; Montaño, A.F.O.; Cidràs, J.; Díaz-Dorado, E. Review of power curve modelling for wind turbines.
Renew. Sustain. Energy Rev. 2013, 21, 572–581.
10. Litong-Palima, M.; Bjerge, M.H.; Cutululis, N.A.; Hansen, L.H.; Sørensen, P. Modeling of the dynamics of
wind to power conversion including high wind speed behavior. Wind Energy 2015, doi:10.1002/we.1876.
11. Mücke, T.A.; Wächter, M.; Milan, P.; Peinke, J. Langevin power curve analysis for numerical wind
energy converter models with new insights on high frequency power performance. Wind Energy 2015, 18,
doi:10.1002/we.1799.
12. Gottschall, J.; Peinke, J. How to improve the estimation of power curves for wind turbines. Environ. Res. Lett.
2008, 3, doi:10.1088/1748-9326/3/1/015005.
13. Allen, J.J.; Shockling, M.; Kunkel, G.; Smits, A.J. Turbulent flow in smooth and rough pipes. Philos. Trans. R.
Soc. A 2007, 365, doi:10.1098/rsta.2006.1939.
14. Gioia, G.; Chakraborty, P. Turbulent friction in rough pipes and the energy spectrum of the phenomenological
theory. Phys. Rev. Lett. 2006, 96, doi:10.1103/PhysRevLett.96.044502.
15. Lueken, C.; Cohen, G.E.; Apt, J. Costs of Solar and Wind Power Variability for reducing CO2 Emissions.
Environ. Sci. Technol. 2012, 46, 9761–9767.
16. Katzenstein, W.; Apt, J. The cost of wind power variability. Energy Policy 2012, 51, 233–243.
17. Efron, B. Bootstrap methods: Another look at the jackknife. Ann. Stat. 1979, 7, 1–26.
18. Calif, R.; Schmitt, F.G. Taylor Law in Wind Energy Data. Resources 2015, 4, 787–795.
19. Taylor, L.R. Aggregation, Variance and the Mean. Nature 1961, 189, 732–735.
20. Taylor, G.I. The spectrum of turbulence. Proc. R. Soc. Lond. A 1938, 164, doi:10.1098/rspa.1938.0032.
21. Eisler, Z.; Bartos, I.; Kertész, J. Fluctuation scaling in complex systems: Taylor’s law and beyond. Adv. Phys.
2008, 57, 89–142.
22. Lovejoy, S.; Schertzer, D.; Stanway, J.D. Direct Evidence of Multifractal Atmospheric Cascades from Planetary
Scales down to 1 km. Phys. Rev. Lett. 2001, 86, 5200–5203.
23. Muzy, J.F.; Baïle, R.; Poggi, P. Intermittency of surface-layer wind velocity series in the mesoscale range.
Phys. Rev. E 2010, 81, doi:10.1103/PhysRevE.81.056308.
24. Apt, J. The spectrum of power from wind turbines. J. Power Sources 2007, 169, 369–374.
25. Calif, R.; Schmitt, R.; Huang, Y. Multifractal description of wind power fluctuations using arbitrary order
Hilbert spectral analysis. Physica A 2013, 392, 4106–4120.
26. Bel, G.; Connaughton, C.P.; Toots, M.; Bandi, M.M. Grid-scale Fluctuations Forecast Error in Wind Power.
New J. Phys. 2016, 18, doi:10.1088/1367-2630/18/2/023015.
27. Benzi, R.; Ciliberto, S.; Tripiccione, R.; Baudet, C.; Massaioli, F.; Succi, S. Extended self-similarity in turbulent
flows. Phys. Rev. E 1993, 48, R29, doi:10.1103/PhysRevE.48.R29.
28. Frisch, U. Turbulence: The legacy of A. N. Kolmogorov; Cambridge University Press: Cambridge, UK, 1995.
29. Garratt, J.R. The Atmospheric Boundary Layer; Cambridge University Press: Cambridge, UK, 1994.
30. Makridis, A.; Chick, J. Validation of a CFD model of wind turbine wakes with terrain effects. J. Wind Eng.
Ind. Aerodyn. 2013, 123, 12–29.
31. Castellani, F.; Astolfi, D.; Garinei, A.; Proietti, S.; Sdringola, P.; Terzi, L.; Desideri, U. How wind turbines
alignment to wind direction affects efficiency? A case study through SCADA data mining. Energy Procedia
2015, 75, 697–703.
32. González-Longatt, F.; Wall, P.; Terzija, V. Wake effect in wind farm performance: Steady-state and dynamic
behavior. Renew. Energy 2012, 39, 329–338.
33. Kim, S.H.; Shin, H.K.; Joo, Y.C.; Kim, K.H. A study of the wake effects on the wind characteristics and fatigue
loads for the turbines in a wind farm. Renew. Energy 2015, 74, 536–543.
34. Wu, Y.T.; Porté-Agel, F. Modeling turbine wakes and power losses within a wind farm using LES:
An application to the Horns Rev offshore wind farm. Renew. Energy 2015, 75, 945–955.
c 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC-BY) license (http://creativecommons.org/licenses/by/4.0/).