Download - Brazil WindPower

Accurate Long Term Wind Resource Assessment
through Multivariate Analysis
Dario Patanè, Mario Benso, Carlos Hernández,
Fernando de la Blanca, Cristóbal López
Ereda
Paseo del Marqués de Monistról 11, Madrid, Spain
[email protected], [email protected], [email protected]
[email protected], [email protected]
ABSTRACT
Estimating the long term variability of wind data is a key factor in the precision for wind
resource assessment. The standard approach, i.e. the Measure-Correlate-Predict method (MCP),
consists in extrapolating from a historical wind speed time series, a long term series for the target
mast located at the site under evaluation.
Data from masts with long measurement periods and numerical weather models data provide
a large amount of potential inputs for MCP. In order to process all available data, a novel
Multivariate MCP (MMCP) methodology was developed. This technique takes efficiently into
account all inputs and extracts the maximum amount of information relevant for the local wind
climate of the target mast.
The multivariate MCP was tested against the standard MCP using a database of 13 long term
masts located in complex terrains in different regions of Spain and a substantial improvement in the
long term wind resource assessment was found.
Keywords: wind resource assessment, long term analysis
1. Introduction
Fluctuations of wind resource are one of the major sources of uncertainty on the energy
production estimate of a potential wind park, thus strongly influencing its financial risk [1]. To deal
with this issue, it is customary to resort to MCP analysis [2, 3]. It consists first in identifying among
the available long term reference wind data the series with the highest correlation with the target
one in the concurrent period. Second, the relationship between the target and reference data is
modeled by means of a single variable regression and eventually a long term series for the target
site is extrapolated. Usually a linear regression is employed to fit the data binned into wind
direction sectors. A variety of MCP algorithms have been proposed in the literature [4-11] to
improve the linear MCP using nonlinear techniques or exploring the joint probability distribution of
data.
Usually long term data from masts located close to the target site or weather numerical
models data, such as NCAR/NCEP reanalysis data are considered [12]. In the latter case several
wind and other variables time series relative to a grid of geographical coordinates and height levels
are available. Every variable, even those not highly correlated to the target mast time series, might
provide valuable pieces of information about the target local wind climate that are usually discarded
in the standard MCP analysis. In order to exploit all available data, a multivariate generalization of
the MCP approach was developed. In the following section we describe the novel MMCP
methodology. In order to compare the latter with the standard technique, long term data of a set of
meteorological masts located in Spain were studied (Section 3). The goodness of the MMCP against
the standard MCP is eventually tested using a set of statistical parameters in Section 4.
2. Method
The MMCP methodology consists first in identifying the available long term masts and the
numerical weather model nodes closest to the target mast. The selected variables may be highly
correlated among themselves and some of them may provide no information about the target local
wind climate. In order to take into account these issues, the MMCP is built as a two steps process
(see Fig. 1). First the input variables are analyzed by means of a feature selection algorithm. In this
step variables bringing little or redundant information about the target local wind climate are
discarded. Second, a multivariate regression similar to principal component analysis is performed,
based on the cross-correlation matrix of the selected variables and the target wind time series [13].
As a result a long term wind speed time series at the target site is extrapolated.
Fig.1 Schematic structure of the MMCP method.
3. Validation
In order to compare MCP and MMCP, a set of meteorological masts located in different
regions of Spain with long measurement periods were studied [14]. A higher number of masts were
initially taken into account, but after being scrutinized to reject anomalous data such as gross trends
or step changes (see Fig. 2) a final number of 13 was selected. In particular, masts close to obstacles
(such as buildings or trees) were discarded. The final set of masts considered is shown in Fig. 3 and
their characteristics are shown in Table 1.
Fig. 2 Anomalous trend of wind velocity V and Weibull shape parameter K for a rejected mast. The
running averages shown in the plot are rescaled with respect to their mean values.
Fig. 3 Locations of the selected masts
Table 1. Details of sites used in the analysis.
Id
Mast
1 Gorramendi
2
Aralar
3
Trinidad
4
Carrascal
5
Tafalla
6 Carcastillo
7
Bardenas
8 LomaNegra
9
Corrubedo
10 CIS Ferrol
11
Muralla
#
UTMAvailability
Available
UTM-y Region
x
%
Years
626330 4785520 Nav
85
16
584586 4756273 Nav
76
16
583286 4740665 Nav
86
16
609766 4726608 Nav
93
16
608897 4708597 Nav
87
16
626578 4692352 Nav
92
16
617067 4673673 Nav
93
16
634369 4658997 Nav
91
16
497771 471638
Gal
97
9.7
560578 4815885 Gal
99
9.7
518427 4732806 Gal
93
8.5
Mean
Speed m/s
7.97
7.32
7.22
6.4
3.26
2.93
5.45
7.27
4.18
3.1
7.05
12
13
Ancares 669916 4743224
Marroxo 623277 4703623
Gal
Gal
94
99
8.8
9.7
4.92
2.34
For each valid mast, NCAR/NCEP reanalysis data relative to the closest geographical nodes
were considered as historical data set. For the MCP technique the reanalysis wind speed time series
with the highest correlation with the mast was considered as input. On the other hand, for the
MMCP analysis, all time series relative to different height levels and geographical coordinates were
used as inputs. Mast data relative to a short time window (a ‘reference period’) were used to
construct a relationship between the mast and the NCAR data. This relationship was then applied to
the long term NCAR data to extrapolate a mast wind speed time series for the entire measurement
period. The estimated series was then compared with the real one by evaluating the parameters
described below. The process was repeated sliding the time window along the measurement period
and averaging the results.
4. Results
In order to compare the goodness of the MMCP regression with respect to the MCP one, the
following parameters were calculated: the Root Mean Squared Error (RMSE) of the daily mean
wind speed and the coefficient of determination R2 between the measured and extrapolated time
series. Moreover, the wind speed distribution obtained from the extrapolated time series was
analyzed by means of the RMSE of the shape factor K of the Weibull distribution. Finally, the
RMSE of the mean wind speed (<V>) with respect to the entire measurement period was addressed,
too.
The results obtained for a meteorological mast by varying the length of the reference period
are shown in Fig 4. Results relative to the entire set of analyzed masts for a reference period of one
year are summarized in Fig. 5. A substantial improvement of the MMCP results with respect to the
standard MCP was found for all parameters considered, especially for the shape factor estimation.
Fig. 4 Results obtained for mast #1. RMSEs of daily wind speed and Weibull K parameter (upper
panel) and long term wind speed (lower panel) as a function of the reference period length used for
the regression. RMSEs are relative to the corresponding long term mean values. Also the correlation
between extrapolated and measured long term time series is shown.
Fig. 5 Results obtained averaging all one year reference periods for the masts considered.
5. Conclusions
Summarizing, a novel multivariate approach to long term resource assessment was proposed.
The MMCP regression (as indicated by RMSE_V and R2) was found on average to be about 16%
more accurate with respect to the standard MCP; this led to an increase of the quality of the
estimation of the long term wind resource about 19% better than the MCP one (RMSE_<V>).
Moreover, a remarkable improvement of the estimation of wind speed distribution (about 51%
RMSE_K) was found.
MMCP can be exploited also within measurement campaigns in complex terrains. In this
case, in fact, the measurement period of an additional mast can be extended extrapolating data from
close masts. Moreover the multivariate approach proposed could be a valuable tool to perform
statistical downscaling of weather numerical model wind forecasts to wind park sites.
MMCP
reveals to be an effective approach in the improvement of precision of long term estimation of wind
resource and hence a tool to reduce the wind resource uncertainty and reduce the risk of financing
wind farms.
REFERENCES
[1] Raftery P, Tindal AJ, Wallenstein M, Johns J, Warren B and Dias Vaz F. Understanding the
risks of financing wind farms. EWEC, 1999.
[2] Mortensen NG. Wind Measurements for Wind Energy Applications: A Review. BWEA, 1994.
[3] Baily B. Wind Resource Assessment Handbook. AWS Scientific, Inc. Report, 1997.
[4] Riedel V, Strack M. Robust approximation of functional relationships between meteorological
data: Alternative measure-correlate-predict algorithms. Proc. EWEA, 2001.
[5] Landberg L, Mortenson NG. A comparison of physical and statistical methods for estimating
the wind resource at a site. BWEA 1993.
[6] Clive P. Non-linearity in MCP with Weibull Distributed Wind Speeds. Wind Eng, 2008; 32, 3.
[7] Mortimer AA. A new correlation/prediction method for potential wind farm sites. Proc. BWEA,
1994.
[8] King C, Hurley B. The SpeedSort, DynaSort and Scatter wind correlation methods, Wind Eng
2005; 29, 3, 217-241.
[9] Garcia-Rojo R. Algorithm for the estimation of the long-term wind climate at a meteorological
mast using a joint probabilistic approach. Wind Eng, 2004; 28, 2, 213-224.
[10]
King C, Hurley B. The Moulded Site Data (MSD) wind correlation method; description and
assessment. Wind Eng 2004; 28, 6, 649-666.
[11]
Rogers AL, Rogers JW, Manwell JF, Comparison of the Performance of Four MeasureCorrelate-Predict Algorithms, JWEIA 2005; 93, 3, 243-264.
[12]
The NCEP/NCAR Reanalysis Project, NOAA/ESRL Physical Sciences Division,
http://www.cdc.noaa.gov/cdc/reanalysis/reanalysis.shtml
[13]
Izenman A J, Modern Multivariate Statistical Techniques: Regression, Classification, and
Manifold Learning, Springer Texts in Statistics, 2008.
[14]
Daily wind speed and direction time series of the analyzed masts were provided by the
Consellería de Medio Ambiente, Xunta de Galicia, and the Departamento de Desarrollo Rural,
Gobierno de Navarra.