Next Article in Journal
An Optimal Reactive Power Control Strategy for a DFIG-Based Wind Farm to Damp the Sub-Synchronous Oscillation of a Power System
Next Article in Special Issue
Batch Growth of Chlorella Vulgaris CCALA 896 versus Semi-Continuous Regimen for Enhancing Oil-Rich Biomass Productivity
Previous Article in Journal
A Simulation Framework for Optimal Energy Storage Sizing
Previous Article in Special Issue
Biofuels from the Fresh Water Microalgae Chlorella vulgaris (FWM-CV) for Diesel Engines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Diagnosis of the Best Weibull Methods for Wind Power Assessment for Agricultural Applications

1
School of Engineering and Technology, Central Queensland University, Rockhampton, QLD 4702, Australia
2
National Centre for Engineering in Agriculture, Faculty of Engineering and Surveying, University of Southern Queensland, Toowoomba, QLD 4350, Australia
*
Author to whom correspondence should be addressed.
Energies 2014, 7(5), 3056-3085; https://doi.org/10.3390/en7053056
Submission received: 20 February 2014 / Revised: 22 April 2014 / Accepted: 23 April 2014 / Published: 2 May 2014
(This article belongs to the Special Issue Renewable Energy for Agriculture)

Abstract

: The best Weibull distribution methods for the assessment of wind energy potential at different altitudes in desired locations are statistically diagnosed in this study. Seven different methods, namely graphical method (GM), method of moments (MOM), standard deviation method (STDM), maximum likelihood method (MLM), power density method (PDM), modified maximum likelihood method (MMLM) and equivalent energy method (EEM) were used to estimate the Weibull parameters and six statistical tools, namely relative percentage of error, root mean square error (RMSE), mean percentage of error, mean absolute percentage of error, chi-square error and analysis of variance were used to precisely rank the methods. The statistical fittings of the measured and calculated wind speed data are assessed for justifying the performance of the methods. The capacity factor and total energy generated by a small model wind turbine is calculated by numerical integration using Trapezoidal sums and Simpson's rules. The results show that MOM and MLM are the most efficient methods for determining the value of k and c to fit Weibull distribution curves.

1. Introduction

Energy and environment are the twin major crises in the world [1]. Because of this, both developed and developing countries are becoming increasingly more interested in using pollution free, cost effective and renewable sources of energy [2]. As a renewable and alternative energy, wind is the most common and fastest-growing source of energy in the world [35]. The characteristics of wind energy are important in different aspects regarding wind energy exploitation [6,7]. Wind is highly variable, both in space and in time [8]. The importance of this variability becomes critical since it is amplified by the cubic relationship of the available power to the wind speed (p = 0.5ρAV3). The wind power production faces the fluctuation of the wind velocity [912]. Therefore, accurate knowledge about the wind characteristics is needed for planning, design and operation of wind turbines [1316]. For the proper assessment, the variability of the wind over time can be divided into three distinct time scales. Firstly, the large time scale variability describes the variations of the amount of wind from one year to another, or even over periods of decades or more. Secondly, the medium time scale covers periods up to a year. These seasonal variations of the wind are much more predictable. Finally, the short term time scale variability covers time scales of minutes to seconds, also well known by the term “turbulence” and which is of critical interest in the wind turbine design process [1723]. For more than half a century the Weibull distribution has attracted the attention of statisticians working on theory and methods as well as various fields of statistics [22,2427]. Hundreds of papers have been written on this distribution; however the research is still ongoing. Together with the normal, exponential distributions, the Weibull distribution is the most popular model in statistics [28,29]. It is of utmost interest to theory orientated statisticians because of its great number of special features, and to practitioners because of its ability to fit to data from various fields, ranging from life data to weather data or observations made in economics and business administration, in health, in physical and social science, in hydrology, in biology or in the engineering sciences [3036].

Research is ongoing worldwide on the Weibull distribution to find the most reliable methods for wind energy estimation. The main question is how precisely the values of the Weibull shape factor “k” and scale factor “c” can be determined [3739]. For this reason different scientists and engineers have developed different methods to find the Weibull parameters for wind energy assessment. That is why several methods are found in literature to estimate the Weibull factors [38,4043]. Recently, Mohammadi and Mostafaeipour [44] used two methods (STDM and PDM) for wind data assessment in Zarrineh, Iran. In 2012, Costa Rocha et al. [45] dealt with the analysis and comparison of seven numerical methods for the assessment of effectiveness in determining the parameters for the Weibull distribution, using wind data collected for Camocim and Paracuru cities in the northeast region of Brazil. Before that, Chang [41] made a statistical study to compare the performance of six numerical methods in estimating Weibull parameters for wind energy application. Seguro and Lambert [43] concluded that the maximum likelihood method (MLM) performs better than the popularly used graphical method (GM). Akdağ and Dinler [40], Azad and Saha [46], and Azad and Alam [47] reviewed three conventional methods and Akdag and Dinler [40] also proposed a new method, namely the energy pattern factor method, for estimating the Weibull parameters. Chu and Ke [48] and Bhattacharya and Bhattacharjee [49] examined the estimation comparison between the MLM and the least squares method [50]. They found that the least squares method significantly outperforms the MLM when sample size is the same. Jowder [51] used empirical and graphical methods to analyze the wind power density at 10, 30, and 60 m height in the Kingdom of Bahrain; two Weibull parameters were estimated and compared. From the analysis, it was found that the empirical methods provide more accurate prediction of average wind speed and power density than the GM. Odo et al. [52] used the Weibull distribution based model for prediction of wind energy potential in Enugu, Nigeria over a period of 13 years. Oyedepo et al. [53] also analyzed south-east Nigeria wind data that spans from 24 years to 37 years and was measured at 10 m height. Abbas et al. [54] statistically analyzed the wind speed data in Pakistan to determine the best fitting distribution of wind speed. For this purpose, two parameters Gamma, Weibull, Lognormal and Rayleigh distributions, and three parameters Burr and Frechet distributions were fitted to data and parameters for each distribution were estimated using the MLM [42,50,5561].

Seguro and Lambart [43] calculated the value of the Weibull parameters by three methods. They recommended that the MLM is useful for time series wind speed data and the modified maximum likelihood method (MMLM) is recommended for use with wind data in frequency distribution format. Philippopoulos and Deligiorgi [62] statistically simulated the wind speed data in Athens, Greece based on the Weibull and autoregressive-moving average (ARMA) method. They found that ARMA methods are superior in simulating the frequency distribution of wind speed. Karpa and Naess [63] also analyzed the extreme value statistics of wind speed by the average conditional exceedance rate (ACER) method. Morgan et al. [64] examined the probability distributions for offshore wind speed. They concluded that the two-parameter lognormal distribution performs best for estimating extreme wind speeds, but still gives estimates with significant error. Stathopoulos et al. [65] used both numerical and statistical models for wind power prediction. Zhou et al. [66] comprehensively evaluated wind speed distribution models for a case study of North Dakota country sites. Wind energy estimation and analysis of wind regions through the Weibull distribution methods is widely used nowadays [35,37,6780].

This study summarizes the results of the 10-min time series wind speed data measured at 20 m and 30 m height in three windy sites, namely Kuakata, Kutubdia and Sitakunda, located in Bangladesh. Table 1 shows the site information in which DL9210 anemometer used to measure wind speed at a 10-min interval. The anemometers have two sensors each of 32 KB memory and have one/two sockets. Statistical work was involved to find the best method of Weibull distribution with high efficiency. The assessment of wind energy potential has been done by using the best method of Weibull distribution. Section 2 offers a detailed outline of the methodology and associated theories for the statistical analysis. The results and discussions are presented in Section 3. Concluding remarks are presented in Section 4.

2. Outline of the Methodology and Associated Theories

To investigate the feasibility of the wind energy resource at any site, there are basically two ways at present to evaluate wind power. The first and the most accurate method to calculate wind power potential is based on measured values that are recorded at meteorological stations. The second method to assess wind power potential is by using probability distribution functions, namely the Rayleigh distribution, Chi-squared distribution, Normal distribution, Binomial distribution, Poisson distribution and Weibull distribution. In this study, the authors used only the Weibull distribution for wind power assessment as presented and discussed below.

2.1. Weibull Probability Density Function

The Weibull probability density function is a two-parameter function characterized by a dimensionless shape parameter (k) and scale parameter (c in m/s). These two parameters determine the wind speed for optimum performance of a wind conversion system as well as the speed range over which the device is likely to operate [52,8487] as given in Equation (1):

f ( v ) = d F ( v ) d v = k c ( v c ) k 1 × e ( v c ) k
where v, k and c are wind speed (m/s), shape factor (dimensionless) and scale factor (m/s), respectively.

2.2. Cumulative Distribution Function or the Weibull Function

Cumulative distribution function is the integration of the Weibull density function. It is the cumulative of relative frequency of each velocity interval [31]. The equation of the Weibull Function is given by:

F ( v ) = 0 v f ( v ) d v
o r F ( v ) = 1 e ( v c ) k

All these distributions are used to determine the probability of occurrence. The nature of the occurrence affects the shape of the probability curve, and in the case of the wind regime, the cumulative curve probability nature mostly fits to the Weibull Function. Several methods to estimate Weibull factors are found in the literature. Some of these methods are:

(1)

Graphical method (GM);

(2)

Method of moments (MOM);

(3)

Standard deviation method (STDM);

(4)

Maximum likelihood method (MLM);

(5)

Power density method (PDM);

(6)

Modified maximum likelihood method (MMLM);

(7)

Equivalent energy method (EEM).

2.2.1. GM

The graph is constructed in such a way that the cumulative Weibull distribution becomes a straight line, with the shape factor k as its slope. Taking the logarithm of both sides, the expression of Equation (3) can be rewritten as:

{ 1 F ( v ) } 1 = e ( v c ) k
or:
ln { 1 F ( v ) } = k ln v k ln c

The above equation represents a relationship between ln(v) and −ln{1 − F(v)}. Therefore, the horizontal axis of this plot on the Weibull paper is v while the vertical axis is ln(1 − F(v))−1. The result is a straight line with slope k. For v = c, one finds F(v) = 1 − e−1 = 0.632 and t an estimation for the value of c, by drawing a horizontal line at F(v) = 0.632. The intersection point with the Weibull line gives the value of c.

2.2.2. MOM

The MOM is another technique commonly used in the field of parameter estimation. If the numbers v1, v2, …, vn represent a set of data, then an unbiased estimator for the n-th origin moment is given by:

m n ^ = 1 n i = 1 n v i n
where m n ^ stands for the estimation of mn. In a Weibull distribution, the nth moment readily follows from Equation (3). With an expression of the Gamma function Γ(x), the average wind speed can be expressed as a function of c and k. The integral found cannot be solved, however it can be reduced to a standard integral, the gamma function, as follows:
Γ ( x ) = 0 α y x 1 × e y d y
where y = ( v c ) k and v c = y x 1 ; x = 1 + 1 k and, after a few manipulations:
μ n = ( 1 c n ) n k Γ ( 1 + n k )

After a few manipulations:

μ n = C × Γ ( 1 + 1 k ) = 0.8525 + 0.0135 k + e [ 2 + 3 ( k 1 ) ]

This formula can easily be handled by pocket calculators in energy output calculations. The accuracy of the approximation is within 0.5% for 1.6 < k < 3.5.Then from Equation (8) we can find the first and the second moments as follows:

m 1 = μ n ^ = ( 1 c ) 1 k Γ ( 1 + 1 k )
and:
m 2 = μ n 2 ^ + σ n 2 ^ = ( 1 c ) 2 k [ Γ ( 1 + 2 k ) { Γ ( 1 + 1 k ) } 2 ]

When we divide m2 by the square of m1, we get an expression which is a function of the shape factor k only:

Γ ( 1 + 2 k ) Γ 2 ( 1 + 1 k ) Γ 2 ( 1 + 1 k )

On taking the square roots of the equation, we have the coefficient of variation (COV):

COV = Γ ( 1 + 2 k ) Γ 2 ( 1 + 1 k ) Γ ( 1 + 1 k )

In this case, this method can be used as an alternative to the MLM. The value of k and c can be easily determined by the following equations:

σ = c [ Γ ( 1 + 2 k ) Γ 2 ( 1 + 1 k ) ] 1 / 2
σ v ¯ = Γ ( 1 + 2 k ) [ Γ ( 1 + 1 k ) ] 2 1

After some calculation we can find:

k = ( 0.9874 σ v ¯ ) 1.0983

The Weibull scale factor can be calculated by:

v ¯ = c Γ ( 1 + 1 k )

2.2.3. STDM

In the STDM, the Weibull factors can be obtained as follows:

k = ( σ v ¯ ) 1.086
and:
c = v ¯ Γ ( 1 + 1 k )
where and σ are mean wind speed and standard deviation of wind speed for any specified periods of time respectively, and can be calculated [44,88] as follows:
v ¯ = 1 n i = 1 n v i

By determining the mean wind speed , the standard deviation σ of wind speed becomes:

σ = [ 0 α ( v v ¯ ) v d v ]
or:
σ = [ ( 1 n 1 i = 1 n ( v v ¯ ) 2 ) ] 1 / 2

One can find next an expression for σ in terms of k and c with = c × Γ(1 + 1/k) and also Γ(x) is the gamma function and is defined as:

Γ ( x ) = 0 α exp ( u ) u x 1 d x

2.2.4. MLM

Maximum likelihood estimation has been the most widely used method for estimating the parameters of the Weibull distribution. The commonly used procedure of MLM proposed by Cohe [50], Harter and Moore [55] and Gove [89] due to its very desirable properties. Let v1, v2, … vn be a random sample of size n drawn from a probability density function f(vi, θ), where θ is an unknown parameter. The likelihood function of this random sample is the joint density of the n random variables and is a function of the unknown parameter [89] given as:

L = i = 1 n f v i ( v i , θ )

Thus, Equation (24) is the likelihood function. The maximum likelihood estimator (MLE) of θ, say θ̂, is the value of θ that maximizes L or, equivalently, the logarithm of L. Often, but not always, the MLE of θ is a solution of (dlogL)/dθ = 0, where solutions that are not functions of the sample values v1, v2, …, vn are not admissible, nor are solutions which are not in the parameter space. Now, we are going to apply the MLE to estimate the Weibull parameters, namely the shape and the scale parameters [90]. Consider the Weibull probability density function given in Equation (3), and then the likelihood function will be represented as:

L ( v i , k , c ) = i = 1 n { ( k c ) ( v i c ) k 1 e x p [ ( v i c ) k ] }

On taking the logarithms of both sides of Equation (25), we obtain the estimating log-likelihood function:

ln ( L ) = n ln ( k ) n k ln ( c ) + ( k 1 ) i = 1 n ln ( v i ) i = 1 n ( v i c ) k

Differentiating Equation (26) with respect to k and c in turn and equating to zero, we have:

ln ( L ) c = n ( k c ) + ( k c ) i = 1 n ( v i c ) k = 0
ln ( L ) k = ( n k ) + i = 1 n ( v i c ) k i = 1 n ( v i c ) k l n ( v i c ) = 0

From Equation (27):

c = [ 1 n i = 1 n ( v i ) k ] 1 / k

When is obtained, then ĉ is can be determined. To solve by using the Newton-Raphson method as given below, let f(k) be the same as Equation (28) and taking the first differential of f(k), we have:

f ( k ) = ( n k 2 ) i = 1 n ( v i c ) k l n 2 ( v i c )

Substituting Equation (29) into Equation (28) gives:

f ( k ) = ( n k ) + i = 1 n [ ( v i ) { ( 1 n ) i = 1 n ( v i ) i k } ] i = 1 n [ ( v i ) k { ( 1 n ) i = 1 n ( v i ) i k } ] l n [ ( v i ) { ( 1 n ) i = 1 n ( v i ) i k } 1 / k ]

Substituting Equation (29) into Equation (30), we get:

f ( k ) = { ( n k ) + i = 1 n [ ( v i ) k { ( 1 n ) i = 1 n ( v i ) k } ] ln 2 [ ( v i ) { ( 1 n ) i = 1 n ( v i ) k } 1 / k ] }

Therefore, is obtained from the equation below by carefully choosing an initial value for ki and iterating the process until it converges:

k i + 1 = k i ( n k ) + i = 1 n [ v i { ( 1 n ) i = 1 n ( v i ) k } 1 / k ] { ( n k ) + i = 1 n [ ( v i ) k { ( 1 n ) i = 1 n ( v i ) k } ] ln 2 [ ( v i ) { ( 1 n ) i = 1 n ( v i ) k } 1 / k ] } i = 1 n [ ( v i ) k { ( 1 n ) i = 1 n ( v i ) k } ] ln [ ( v i ) { ( 1 n ) i = 1 n ( v i ) k } 1 / k ] { ( n k ) + i = 1 n [ ( v i ) k { ( 1 n ) i = 1 n ( v i ) k } ] ln 2 [ ( v i ) { ( 1 n ) i = 1 n ( v i ) k } 1 / k ] }

2.2.5. PDM

To obtain the shape factor and scale factor through this method, firstly the energy pattern factor is computed. The energy pattern factor usage is for turbine aerodynamic design. The energy pattern factor is related to the averaged data of wind speed and is defined as a ratio between mean of cubic wind speed to cube of mean wind speed. The energy pattern factor Epf is expressed as [40,45]:

E pf = Total amount of power available in the wind Power calculated by cubing the mean wind speed or E pf = 1 n i = 1 n v i 3 ( 1 n i = 1 n v i ) 3 = v 3 ¯ ( v ¯ ) 3 = Γ ( 1 + 3 k ) Γ 3 ( 1 + 1 k )

Once the energy pattern factor is calculated by using the above equation, the Weibull shape factor and scale factor can be estimated from the following formulas:

k = 1 + 3.69 E pf 2
c = v ¯ Γ ( 1 + 1 k )
where Epf is the energy pattern factor; and Γ is the gamma function.

2.2.6. MMLM

The MMLM can only be considered if the available data of wind speed are already in the shape of the Weibull distribution. The solution of the equations in the MLM requires some numerical iteration by the Newton-Raphson method [41]. The Weibull parameters are determined by the following equations:

k = [ i = 1 n v i k l n ( v i ) f ( v i ) i = 1 n v i k f ( v i ) i = 1 n l n ( v i ) f ( v i ) f ( v 0 ) ] 1
c = [ 1 f ( v 0 ) i = 1 n ( v i ) k f ( v i ) ] 1 / k
where vi is the wind speed central to bin i; n is the number of bins; f(vi) represents the Weibull frequency for the wind speed range within bin I; and f(v0) is the probability for wind speed to equal or exceed zero.

2.2.7. EEM

Consider a random sample of v1, v2, …, vn by relative frequency of occurrence in a given interval of wind speed. Then a random variable observation discrete Wv associated with wind speed can be obtained from:

W v = f v ( v )

This observation random is also related with the Weibull parameters k and c from the equation of the probability of occurrence W. W(v) is the observed frequency of the wind speed for the interval of (v − 1)V < v. The mathematical representation of W(v) is:

W ( v ) = Q ( v 1 ) Q ( v )
W ( v ) = { e ( v 1 c ) k e ( v c ) k }
where Q(v) and the probability of occurrence of wind speeds equal to or higher than v, given by Q(v) = 1 − F(v). Then the random observation Wv can be written using Equations (39) and (41) as:
W v = W ( v ) + ε v = { e ( v 1 c ) k e ( v c ) k } + ε v

The first hypothesis says that “The energy density is a parameter that helps in the determination of parameters of the Weibull distribution for applications in wind energy”. The related factor part deterministic must meet the following conditions: (a) be variable with random value expected value equal to 0: E(ε) = 0; (b) be variable random variance with constant: v(ε) = σ2; and (c) the occurrences of ε are non-correlated: COVi, εj). To ensure the condition of equivalence initially proposed in the hypothesis, it is the equality between Equation (4) and v m 3 = 1 n i = 1 n v i 3. The equation resulting from this transaction expresses the parameter c as a function of speed cubed average of observations and the parameter k as:

c = ( v m 3 Γ ( 1 + 3 k ) ) 1 / 3

Substituting Equation (43) into Equation (42) gives:

W v = { e [ ( v 1 ) { Γ ( 1 + 3 k ) } 1 / 3 ( v m 3 ) 1 / 3 ] k e [ ( v ) { Γ ( 1 + 3 k ) } 1 / 3 ( v m 3 ) 1 / 3 ] k } + ε v

The estimate of the parameter k may be obtained from an estimator of least squares given by the expression [45]:

i = 1 n [ W v i e [ ( v i 1 ) { Γ ( 1 + 3 k ) } 1 / 3 ( v m 3 ) 1 / 3 ] k + e [ ( v i ) { Γ ( 1 + 3 k ) } 1 / 3 ( v m 3 ) 1 / 3 ] k ] 2 = i = 1 n ( ε v i ) 2
where Wvi is the observed frequency of the wind speed; n is the number of intervals of the histogram of speed; vi is the value of the upper limit of the i-th speed interval; v m 3 is the mean of the cubic wind speed; and εvi is the error of the approximation. Once the value of the parameter k is calculated, the value of c can be obtained directly from Equation (43).

2.3. COV

The COV is defined as the ratio between mean standard deviation to mean wind speed expressed as a percentage. It demonstrates the mutability of wind speed and can be expressed as [29]:

COV ( % ) = σ v ¯ × 100

2.4. Wind Speed Varies with Altitude

Wind energy is indirect solar energy because it is generated by the temperature difference between the equator and the poles which drives the thermal system by solar radiation. It is known that wind speed varies with altitude, however, wind blows relatively slowly at low altitude and wind speed then increases with altitude. Different relationships are found in the literature to calculate wind speed at any height [9193]. The Weibull factors used for these calculations must be obtained from the best possible method of Weibull distribution. For this purpose, the Weibull factors are initially calculated at desired heights, then wind speed and wind power are obtained [94,95]. The calculation procedures use the following relationships:

η = 0.37 0.0881 ln c 10
k h = k 10 [ 1 0.0881 ln ( z h z ref ) ]
c h = c 10 ( z h z ref ) n
v h = c h Γ ( 1 + 1 k h )
P ¯ h = 1 2 ρ v ¯ h 3 = 1 2 ρ v ¯ h 3 Γ ( 1 + 3 k h ) Γ 3 ( 1 + 3 k h ) = 1 2 ρ c h 3 Γ ( 1 + 3 k h )
E D = 1 2 ρ c h 3 Γ ( 1 + 3 k h ) × T

The mean energy density over a period of time, T, is the product of mean power density and the time period. k10 and c10 are the shape factor and scale factor at a height of 10 m and η is the power law coefficient. kh, ch, vh, zh and Ph are shape factor, scale factor, wind speed and wind power at the desired height, respectively. The zref is the reference height. ρ is the air density and for standard conditions (i.e., at sea level with temperature of 15 °C and pressure of 1 atmosphere) and is equal to 1.225 kg/m3 and v is the wind speed (m/s).In this paper the authors used the Weibull factors to find wind speed and wind power at a height of 50 m.

2.5. Statistical Error Analysis/Goodness of Fit

To find the best method for the analysis, some statistical parameters were used to analyze the efficiency of the above mentioned methods. The following tests were used to achieve this goal:

(a)

zRelative percentage of error (RPE)

RPE = ( x i , w y i , m y i , m ) × 100 %

(b)

Root mean square error (RMSE)

RMSE = [ 1 N i = 1 N ( y i , m x i , w ) 2 ] 1 / 2

(c)

Mean percentage error (MPE)

MPE = 1 N i = 1 N ( x i , w y i , m y i , m ) × 100 %

(d)

MAPE = 1 N i = 1 N | x i , w y i , m y i , m | × 100 %

(e)

Chi-square error

χ 2 = i = 1 N ( y i , m x i , w ) 2 x i , w

(f)

Analysis of variance or efficiency of the method

R 2 = i = 1 N ( y i , m z i , v ¯ ) 2 i = 1 N ( y i , m x i , w ) 2 i = 1 N ( y i , m z i , v ¯ ) 2
where N is the number of observations; yi,m is the frequency of observation or i-th calculated value from measured data; xi,w is the frequency of Weibull or i-th calculated value from the Weibull distribution; zi,v̄ is the mean of i-th calculated value from measured data. RPE shows the percentage deviation between the calculated values from the Weibull distribution and the calculated values from measured data. MPE shows the average of percentage deviation between the calculated values from the Weibull distribution and the calculated values from measured data, and MAPE shows the absolute average of percentage deviation between the calculated values from the Weibull distribution and the calculated values from measured data. Best results are obtained when these values are close to zero. R2 determines the linear relationship between the calculated values from the Weibull distribution and the calculated values from measured data. The ideal value of R2 is equal to 1 [30,44,45,96,97].

3. Results and Discussion

In this statistical analysis, data from three wind monitoring stations were used to diagnose the best method of the Weibull distribution. The most important results of this analysis, based on hourly, monthly, seasonal and annual figures, are presented. Table 2 shows a lot of information about the sites based on monthly mean as well as annual mean wind speed, number of observations in hours, maximum and minimum wind speed, standard deviation, and COV/turbulence intensity at 20 m and 30 m height.

A sample data frequency distribution and cumulative frequency distribution has been presented in Table 3.

The hourly variation of wind speed at the selected altitude has been presented in Figure 1. Wind speed is higher during the daytime than the nighttime at every site. The COV/turbulence with the time series is at its minimum during the higher windy periods.

The hourly mean COV is presented in Figure 2. An important issue which has been clearly shown in Figure 2 is that turbulence is comparatively low at 30 m height than at 20 m height. Therefore, at the higher level the wind velocity stream is more uniform, i.e., the COV is lower.

From the sample frequency distribution (Table 3), it can be clearly seen that more than 50% of the frequency is between 5 m/s and 13 m/s of wind speed at every station. Table 3 can be used to clearly identify the total number of hours at certain wind speed available in a month. Anyone can find monthly wind speed scenario at a glance from the Table. The results are similar for other months at each station. Figures 1 and 2, shows the hourly mean variations of wind speed and coefficient of variation for the period of 24 h as a mean of a year. It has been shown that the wind speed is low until 8:00 AM to 10:00 AM, and after that it increases during the day until the maximum value is gained from 2:00 PM to 4:00 PM, and afterwards it decrease again until the end of the day. The same wind characteristics have been found at 20 m and 30 m height at every station. The maximum and minimum values of hourly mean wind speeds at 30 m height are 4.95 m/s at Station-III and 3.10 m/s at Station-II. The variation of wind speed at Station-I is lower than the others. The COV is lower when the wind speeds become a maximum. Its values are very high at nighttime. From Table 3 and Figures 1 and 2, it is clear that 30 m height wind speed data at Station-I shows better results than the other sites. To find the best Weibull method, the statistical analysis of the Weibull distribution at 20 m and 30 m height is presented in Tables 4, 5, 6, 7, 8 and 9.

In the statistical analysis, seven methods were used to determine the shape parameter k and scale parameter c of the Weibull distribution. For comparison of these seven methods to each other and to find out the efficiency of the methods, six statistical tools were used, i.e., relative percentage error (RPE), RMSE, MPE, MAPE, chi-square error (χ2), and analysis of variance or efficiency of the method (R2). It is important to note that only one column is required to rank the methods, since the above six criteria all gave the same relative results. For a more precisely diagnosis, the authors used these six statistical tools to rank the methods. The statistical analysis results for the seven numerical methods and the six statistical test results are shown in Tables 4, 5, 6, 7, 8 and 9 for the three stations at 20 m and 30 m height, respectively. The results show that the method of moments (MOM) and MLM give better results than other methods, where the value of k and c becomes almost the same by these two methods. For MLM, RPE and RMSE becomes zero and the value of MPE becomes negative, because, by this method, the calculated value becomes greater than the measured one. But the most important statistical test gives a chi-square error of χ2 = 0.0010 and the efficiency of the method is R2 = 0.9997, where the best results are obtained when these values are close to zero and unity, respectively. For Station-II, in Table 7, MLM and PDM methods show better performance than others. The common factor is that the MLM gives the better performance at every site. Therefore, it can be said that the MLM is the fisrt and MOM is the second most efficient method for wind data assessment at 30 m height.

From Table 4, it can be clearly seen that MOM (χ2 = 0.0001, R2 = 0.9999), MLM (χ2 = 0.0005, R2 = 0.9998) and MMLM (χ2 = 0.0002, R2 = 0.9999) give very close results and show better performance than other methods. In Table 6, MOM and PDM methods show better results than others. In the case of Station-III, MLM and PDM methods show better results than others in Table 8. Summarizing these results, it can be easily said that MOM and MLM methods were applicable for every height at any location wind data assessment. The PDM method is better for wind data assessment at the lower height, but not perfect for the higher height.

Table 10 summaries the test results of the seven methods and ranking of the methods according to their performance and efficiency in wind data assessment. The rankings were done by considering minimum error and maximum efficiency according to first to seventh positions respectively. Regarding this test, six statistical tools have been used and considered four decimal places of each value by numerical iteration methods. One test is enough to rank the methods, but for more precise analysis, the authors used more tools which helped to verify the discussion about the best method. In this statistical research work, it has been found that the MOM achieved the first position and the MLM took the second position in rank. Although the PDM got the third position, this method has better performance for low height wind data assessment. At increased height, PDM was a less efficient method than others. But MOM and MLM methods were applicable at any altitude with minimum error and maximum efficiency. Our first goal has been satisfied by the above statistical analysis where we identified the best methods to determine the Weibull distribution. Another goal is to select the best wind site by using these best methods, which has been analyzed below.

The procedure mentioned in this study is not only applicable in case study sits, it can be applied in any climatic conditions at any site in any countries in the world. For example, Mohammadi and Mostafaeipour [44] used STDM and PDM for wind turbine utilization in Zarriuneh at Kurdistan mountainous province in Iran. Chang [41] used moment method, empirical method, GM, MLM, MMLM and energy pattern factor method at three wind farms (Dayuan, Hengchun and Penghu) experiencing different weather conditions in Taiwan. Costa Rocha et al. [45] used EEM, moment method, MLM, etc. for wind energy generation in coastal area of the State of Ceara, located in the northeast region of Brazil. Dorvlo [98] used moment method, Chi-square method, empirical method and regression method for estimation of wind energy in four weather stations (Marmul, Masirah, Sur and Thumrait) in Oman. Many authors have already used the methods at different geographical locations for wind energy estimation [30,43,99]. The paper presented more generalized form of the Weibull distribution methods which is validated by widely acceptable error correction methods for practical application. This procedure is applicable at any geographical location; at any weather condition (i.e., summer, winter, spring etc.) at any altitude in any country in the world. The main outcome of the study is the development of the processes for identifying best methods for wind power generation which is applicable to any location.

Table 11 presents the variation of wind speed during different months of the year and the annual mean values. From the average values, the variation of the maximum and minimum wind speed occurred in April to September and October to March respectively at every station. At Station-I, the maximum and minimum wind speed occurred in June and October with values of 5.07 m/s and 2.57 m/s, respectively. The COV and standard deviation at this site varies between 21.01%–35.68% and 0.54–2.13, respectively. On the other hand, Station-III has wind speeds similar to Station-I, but the COV and standard deviation varies between 15.93%–129.53% and 0.71–13.07, respectively. In August and September, the COV is 50.1% and 129.53%, which means that these months have some irregular wind behavior, i.e., gusty wind exists during these months. Station-II has wind behavior (wind speed) performance lower than other stations. Therefore, Station-I has shown better wind characteristics than the other Stations, hence it is the selected wind site. The mean available power and energy in the wind at 30 m height at Station-I is analyzed below.

The monthly mean available power is analyzed using both MOM and MLM methods in Table 12. From this table, it can be seen that the available power from April to September was above 100 W/m2 and the maximum and minimum power occurred in August, 179.36 W/m2, and October, 11.8 W/m2, respectively. Only two statistical tools (RMSE and Chi-square error) have been used for comparing the MOM and MLM methods. Here, the MLM method gives better performance with minimum error. Using the MLM method, the Weibull shape factor and scale factor vary between 2.52–5.65 and 2.79–6.69 m/s. respectively.

Therefore, the six months from April to September show more potential wind generated power than other months in the year. In this research work, the Weibull parameters at heights of 30 m and 20 m were determined.

Table 13 presents the seasonal analysis of wind speed, Weibull parameters and available power at 20 m and 30 m height. It also presents the extrapolated values at 50 m. The months in each season can be classified as winter (November, December and January), spring (February, March and April), summer (May, June and July), and autumn (August, September and October).

The mean velocity, shape factor, scale factor, etc. at 50 m were extrapolated using Equations (44)(48) and these extrapolated results are presented in Table 14.

Table 15 also includes yearly or annual mean values. According to this table, the maximum and minimum values of wind speed and wind power have been observed in summer and spring with the values of 5.26 m/s and 4.03 m/s and the power is 108.24 W/m2 and 48.16 W/m2 respectively at 30 m height. In the literatures, Elliot et al. [100], Yu et al. [101], Ilinca et al. [102] and Zhou et al. [103] classified the wind power in seven categories which are shown in Table 14. Considering these wind power classes, it has been found that Station-I has poor wind power at every height (i.e., 20 m, 30 m and 50 m). For a more realistic analysis, the available energy in the wind of this site was determined by using numerical iteration methods and is presented in Table 15.

The energy estimation of wind regimes by the Weibull based approach has been presented in Table 15. It also shows the most frequent wind speed, velocity contributing the maximum energy, energy density and total energy intensity. From the above table, it can be clearly seen that the energy density and total energy intensity are higher in summer than during autumn seasons and also are a minimum in spring. Therefore, summers have more wind potential than any other seasons. To find out the energy generated by a wind turbine, a small model wind turbine (NACA 4418 turbine power curve) with rated power of 20 kW or 0.02 MW, rated speed of 8 m/s and hub height of 30 m has been used for this calculation. Values of capacity factor and total energy output were obtained using numerical integration methods of Trapezoidal Sums and Simpson's 1/3 Rule and the values were verified to four decimal places. From Table 16, it can be seen that, in summer, 13.2938 MW·h of energy can be extracted by the suggested model wind turbine. This means that only 147.67 kW/day of power can be generated by the wind turbine in summer and the value becomes 50% less in the rainy season. Therefore, it can be finally concluded that the analytical results show that Station-I does not have the necessary available wind power potential for large turbines but has sufficient wind power for small wind turbine both for electricity generation and pumping water for irrigation.

4. Conclusions

In this work, statistical diagnosis of the best Weibull distribution methods for wind data analysis is presented. By using the available wind data, the values of shape factor k and scale factor c were determined using seven methods and were then investigated as to how efficiently the methods can estimate the Weibull factors with minimum error. To satisfy the main objectives of this work, six statistical tools were used to find the best method of Weibull distribution. It is important to note that any one of these statistical tools, namely RPE, RMSE, MPE, MAPE, chi-square error (χ2), and analysis of variance or efficiency of the method (R2) is good enough to rank the methods, however, analysis using all of them was done to rank more precisely. The results show that the MOM and MLM are the most efficient methods for determining the value of k and c to fit the Weibull distribution curves. The PDM is more efficient for low altitude wind data but is not efficient for higher altitude wind data. MOM and MLM methods are more efficient with less error and are applicable for any altitude. Other methods such as MMLM, EEM, GM and STDM are the least efficient methods to fit the Weibull distribution curves for the assessment of wind speed data. Another objective of this work was to find the best wind site using the best Weibull distribution methods and calculate available wind power. Monthly mean wind speed was found to be relatively higher in Station-I than that of the other sites. As a result this is our selected site. The MLM method has shown better results than MOM in calculating monthly mean wind power at Station-I. The value of shape factor, scale factor, wind speed and power was determined at 50 m height using extrapolation of numerical equations to satisfy the wind power classes as discussed in this paper. The poor class wind power has been found in this site in each altitude. Furthermore, energy density and total energy intensity per unit area has been analyzed by numerical iteration methods. Finally, the energy extracted by a small model wind turbine has been analyzed by using numerical integration methods of Trapezoidal Sums and Simpson's 1/3 Rule. This study offers a new pathway on how to evaluate feasible locations for wind energy assessment which is applicable at any windy sites in any country in the world.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baños, R.; Manzano-Agugliaro, F.; Montoya, F.G.; Gil, C.; Alcayde, A.; Gómez, J. Optimization methods applied to renewable and sustainable energy: A review. Renew. Sustain. Energy Rev. 2011, 15, 1753–1766. [Google Scholar]
  2. Cruz-Peragon, F.; Palomar, J.M.; Casanova, P.J.; Dorado, M.P.; Manzano-Agugliaro, F. Characterization of solar flat plate collectors. Renew. Sustain. Energy Rev. 2012, 16, 1709–1720. [Google Scholar]
  3. Acker, T.L.; Williams, S.K.; Duque, E.P.; Brummels, G.; Buechler, J. Wind resource assessment in the state of Arizona: Inventory, capacity factor, and cost. Renew. Energy 2007, 32, 1453–1466. [Google Scholar]
  4. Munteanu, I.; Bratcu, A.I.; Cutululis, N.A.; Ceanga, E. Optimal Control of Wind Energy Systems: Towards A Global Approach; Springer: London, UK, 2008. [Google Scholar]
  5. Montoya, F.G.; Aguilera, M.J.; Manzano-Agugliaro, F. Renewable energy production in Spain: A review. Renew. Sustain. Energy Rev. 2014, 33, 509–531. [Google Scholar]
  6. Hernández-Escobedo, Q.; Manzano-Agugliaro, F.; Zapata-Sierra, A. The wind power of Mexico. Renew. Sustain. Energy Rev. 2010, 14, 2830–2840. [Google Scholar]
  7. Hernández-Escobedo, Q.; Saldaña-Flores, R.; Rodríguez-García, E.R.; Manzano-Agugliaro, F. Wind energy resource in Northern Mexico. Renew. Sustain. Energy Rev. 2014, 32, 890–914. [Google Scholar]
  8. Hernandez-Escobedo, Q.; Manzano-Agugliaro, F.; Gazquez-Parra, J.A.; Zapata-Sierra, A. Is the wind a periodical phenomenon? The case of Mexico. Renew. Sustain. Energy Rev. 2011, 15, 721–728. [Google Scholar]
  9. Vincent, C.L.; Giebel, G.; Pinson, P.; Madsen, H. Resolving nonstationary spectral information in wind speed time series using the Hilbert-Huang transform. J. Appl. Meteorol. Climatol. 2010, 49, 253–267. [Google Scholar]
  10. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar]
  11. Chen, Z.; Gomez, S.A.; Mccormick, M. A Fuzzy Logic Controlled Power Electronic System for Variable Speed Wind Energy Conversion Systems. Proceedings of the 8th International Conference on Power Electronics and Variable Speed Drives, London, UK, 18–19 September 2000; pp. 114–119.
  12. Ahmed Shata, A.; Hanitsch, R. Electricity generation and wind potential assessment at Hurghada, Egypt. Renew. Energy 2008, 33, 141–148. [Google Scholar]
  13. Akhmatov, V. Influence of wind direction on intense power fluctuations in large offshore windfarms in the North Sea. Wind Eng. 2007, 31, 59–64. [Google Scholar]
  14. Sorensen, P.; Cutululis, N.A.; Vigueras-Rodríguez, A.; Jensen, L.E.; Hjerrild, J.; Donovan, M.H.; Madsen, H. Power fluctuations from large wind farms. IEEE Trans. Power Syst. 2007, 22, 958–965. [Google Scholar]
  15. Chen, Z.; Blaabjerg, F. Wind farm—A power source in future power systems. Renew. Sustain. Energy Rev. 2009, 13, 1288–1300. [Google Scholar]
  16. Hsieh, C.-H.; Dai, C.-F. The analysis of offshore islands wind characteristics in Taiwan by Hilbert–Huang transform. J. Wind Eng. Ind. Aerodyn. 2012, 107, 160–168. [Google Scholar]
  17. Belu, R.; Koracin, D. Wind characteristics and wind energy potential in western Nevada. Renew. Energy 2009, 34, 2246–2251. [Google Scholar]
  18. Chang, T.P. Wind energy assessment incorporating particle swarm optimization method. Energy. Convers. Manag. 2011, 52, 1630–1637. [Google Scholar]
  19. Adaramola, M.; Paul, S.; Oyedepo, S. Assessment of electricity generation and energy cost of wind energy conversion systems in north-central Nigeria. Energy. Convers. Manag. 2011, 52, 3363–3368. [Google Scholar]
  20. Ohunakin, O.S.; Akinnawonu, O.O. Assessment of wind energy potential and the economics of wind power generation in Jos, Plateau State, Nigeria. Energy. Sustain. Dev. 2012, 16, 78–83. [Google Scholar]
  21. Ohunakin, O.S. Assessment of wind energy resources for electricity generation using WECS in North-Central region, Nigeria. Renew. Sustain. Energy Rev. 2011, 15, 1968–1976. [Google Scholar]
  22. Akpinar, E.; Akpinar, S. Statistical analysis of wind energy potential on the basis of the Weibull and Rayleigh distributions for Agin-Elazig, Turkey. Proc. Inst. Mech. Eng. Part A J. Power Energy 2004, 218, 557–565. [Google Scholar]
  23. Ahmed, A.S. Wind energy as a potential generation source at Ras Benas, Egypt. Renew. Sustain. Energy Rev. 2010, 14, 2167–2173. [Google Scholar]
  24. Fadare, D. A statistical analysis of wind energy potential in Ibadan, Nigeria, based on Weibull distribution function. Pac. J. Sci. Technol. 2008, 9, 110–119. [Google Scholar]
  25. Akpinar, E.K. A statistical investigation of wind energy potential. Energy Sources Part A 2006, 28, 807–820. [Google Scholar]
  26. Ahmed Shata, A.; Hanitsch, R. Evaluation of wind energy potential and electricity generation on the coast of Mediterranean Sea in Egypt. Renew. Energy 2006, 31, 1183–1202. [Google Scholar]
  27. Akdağ, S.A.; Güler, Ö. Evaluation of wind energy investment interest and electricity generation cost analysis for Turkey. Appl. Energy 2010, 87, 2574–2580. [Google Scholar]
  28. Azad, A. Statistical Weibull distribution analysis for wind power of the two dimensional ridge areas. Int. J. Adv. Renew. Energy Res. 2013, 1, 8–14. [Google Scholar]
  29. Ahmed, S.; Mahammed, H. A statistical analysis of wind power density based on the Weibull and Ralyeigh models of “Penjwen Region” Sulaimani/Iraq. Jordan J. Mech. Ind. Eng. 2012, 6, 135–140. [Google Scholar]
  30. Celik, A.N. A statistical analysis of wind power density based on the Weibull and Rayleigh models at the southern region of Turkey. Renew. Energy 2004, 29, 593–604. [Google Scholar]
  31. Kavak Akpinar, E.; Akpinar, S. A statistical analysis of wind speed data used in installation of wind energy conversion systems. Energy Convers. Manag. 2005, 46, 515–532. [Google Scholar]
  32. Safari, B.; Gasore, J. A statistical investigation of wind characteristics and wind energy potential based on the Weibull and Rayleigh models in Rwanda. Renew. Energy 2010, 35, 2874–2880. [Google Scholar]
  33. Kose, R.; Ozgur, M.A.; Erbas, O.; Tugcu, A. The analysis of wind data and wind energy potential in Kutahya, Turkey. Renew. Sustain. Energy Rev. 2004, 8, 277–288. [Google Scholar]
  34. Zaharim, A.; Najid, S.K.; Razali, A.M.; Sopian, K. Analyzing Malaysian Wind Speed Data Using Statistical Distribution. Proceedings of the 4th IASME/WSEAS International Conference on Energy and Environment, Cambridge, UK, 24–26 February 2009; pp. 360–370.
  35. Islam, M.; Saidur, R.; Rahim, N. Assessment of wind energy potentiality at Kudat and Labuan, Malaysia using Weibull distribution function. Energy 2011, 36, 985–992. [Google Scholar]
  36. Daoo, V.; Panchal, N.; Sunny, F.; Sitaraman, V.; Krishnamoorthy, T. Assessment of wind energy potential of Trombay, Mumbai (19.1 N; 72.8 E), India. Energy Convers. Manag. 1998, 39, 1351–1356. [Google Scholar]
  37. Celik, A. Weibull representative compressed wind speed data for energy and performance calculations of wind energy systems. Energy Convers. Manag. 2003, 44, 3057–3072. [Google Scholar]
  38. Ulgen, K.; Hepbasli, A. Determination of Weibull parameters for wind energy analysis of Izmir, Turkey. Int. J. Energy Res. 2002, 26, 495–506. [Google Scholar]
  39. Langlois, R. Estimation of Weibull parameters. J. Mater. Sci. Lett. 1991, 10, 1049–1051. [Google Scholar]
  40. Akdağ, S.A.; Dinler, A. A new method to estimate Weibull parameters for wind energy applications. Energy Convers. Manag. 2009, 50, 1761–1766. [Google Scholar]
  41. Chang, T.P. Performance comparison of six numerical methods in estimating Weibull parameters for wind energy application. Appl. Energy 2011, 88, 272–282. [Google Scholar]
  42. Minami, F.; Brückner-Foit, A.; Munz, D.; Trolldenier, B. Estimation procedure for the Weibull parameters used in the local approach. Int. J. Fract. 1992, 54, 197–210. [Google Scholar]
  43. Seguro, J.V.; Lambert, T.W. Modern estimation of the parameters of the Weibull wind speed distribution for wind energy analysis. J. Wind Eng. Ind. Aerodyn. 2000, 85, 75–84. [Google Scholar]
  44. Mohammadi, K.; Mostafaeipour, A. Using different methods for comprehensive study of wind turbine utilization in Zarrineh, Iran. Energy Convers. Manag. 2013, 65, 463–470. [Google Scholar]
  45. Costa Rocha, P.A.; de Sousa, R.C.; de Andrade, C.F.; da Silva, M.E.V. Comparison of seven numerical methods for determining Weibull parameters for wind energy generation in the northeast region of Brazil. Appl. Energy 2012, 89, 395–400. [Google Scholar]
  46. Azad, A.; Saha, M. Weibull's analysis of wind power potential at coastal sites in Kuakata, Bangladesh. Int. J. Energy Mach. 2011, 4, 36–45. [Google Scholar]
  47. Azad, A.; Alam, M. A Statistical tools for clear energy: Weibull's distribution for potentiality analysis of wind energy. Int. J. Adv. Renew. Energy Res. 2012, 1, 240–247. [Google Scholar]
  48. Chu, Y.-K.; Ke, J.-C. Computation approaches for parameter estimation of Weibull distribution. Math. Comput. Appl. 2012, 17, 39–47. [Google Scholar]
  49. Bhattacharya, P.; Bhattacharjee, R. A study on Weibull distribution for estimating the parameters. Wind Eng. 2009, 33, 469–476. [Google Scholar]
  50. Cohen, A.C. Maximum likelihood estimation in the Weibull distribution based on complete and on censored samples. Technometrics 1965, 7, 579–588. [Google Scholar]
  51. Jowder, F.A. Wind power analysis and site matching of wind turbine generators in Kingdom of Bahrain. Appl. Energy 2009, 86, 538–545. [Google Scholar]
  52. Odo, F.; Offiah, S.; Ugwuoke, P. Weibull distribution-based model for prediction of wind potential in Enugu, Nigeria. Pelagia Res. Libr. Adv. Appl. Sci. Res. 2012, 3, 1202–1208. [Google Scholar]
  53. Oyedepo, S.O.; Adaramola, M.S.; Paul, S.S. Analysis of wind speed data and wind energy potential in three selected locations in south-east Nigeria. Int. J Energy Environ. Eng. 2012, 3. [Google Scholar] [CrossRef]
  54. Abbas, K.; Alamgir.; Khan, S.A.; Ali, A.; Khan, D.M.; Khalil, U. Statistical analysis of wind speed data in Pakistan. World Appl. Sci. J. 2012, 18, 1533–1539. [Google Scholar]
  55. Harter, H.L.; Moore, A.H. Maximum-likelihood estimation of the parameters of gamma and Weibull populations from complete and from censored samples. Technometrics 1965, 7, 639–643. [Google Scholar]
  56. Odell, P.M.; Anderson, K.M.; D'agostino, R.B. Maximum likelihood estimation for interval-censored data using a Weibull-based accelerated failure time model. Biometrics 1992, 48, 951–959. [Google Scholar]
  57. Choi, S.; Wette, R. Maximum likelihood estimation of the parameters of the gamma distribution and their bias. Technometrics 1969, 11, 683–690. [Google Scholar]
  58. Cacciari, M.; Mazzanti, G.; Montanari, G. Comparison of maximum likelihood unbiasing methods for the estimation of the Weibull parameters. IEEE Trans. Dielectr. Electr. Insul. 1996, 3, 18–27. [Google Scholar]
  59. Zanakis, S.H.; Kyparisis, J. A review of maximum likelihood estimation methods for the three-parameter Weibull distribution. J. Stat. Comput. Simul. 1986, 25, 53–73. [Google Scholar]
  60. Lemon, G.H. Maximum likelihood estimation for the three parameter Weibull distribution based on censored samples. Technometrics 1975, 17, 247–254. [Google Scholar]
  61. Jiang, S.; Kececioglu, D. Maximum likelihood estimates, from censored data, for mixed-Weibull distributions. IEEE Trans. Reliab. 1992, 41, 248–255. [Google Scholar]
  62. Philippopoulos, K.; Deligiorgi, D. Statistical simulation of wind speed in Athens, Greece based on Weibull and ARMA models. Int. J. Energy Environ. 2009, 3, 151–158. [Google Scholar]
  63. Karpa, O.; Naess, A. Extreme value statistics of wind speed data by the ACER method. J. Wind Eng. Ind. Aerodyn. 2013, 112, 1–10. [Google Scholar]
  64. Morgan, E.C.; Lackner, M.; Vogel, R.M.; Baise, L.G. Probability distributions for offshore wind speeds. Energy Convers. Manag. 2011, 52, 15–26. [Google Scholar]
  65. Stathopoulos, C.; Kaperoni, A.; Galanis, G.; Kallos, G. Wind power prediction based on numerical and statistical models. J. Wind Eng. Ind. Aerodyn. 2013, 112, 25–38. [Google Scholar]
  66. Zhou, J.; Erdem, E.; Li, G.; Shi, J. Comprehensive evaluation of wind speed distribution models: A case study for North Dakota sites. Energy Convers. Manag. 2010, 51, 1449–1458. [Google Scholar]
  67. Sardar Maran, P.; Ponnusamy, R. Wind power density estimation using meteorological tower data. Int. J. Renew. Sustain Energy. 2013, 2, 110–114. [Google Scholar]
  68. Bivona, S.; Bonanno, G.; Burlon, R.; Gurrera, D.; Leone, C. Stochastic models for wind speed forecasting. Energy Convers. Manag. 2011, 52, 1157–1165. [Google Scholar]
  69. Daniel, A.; Chen, A. Stochastic simulation and forecasting of hourly average wind speed sequences in Jamaica. Sol. Energy. 1991, 46, 1–11. [Google Scholar]
  70. Kamal, L.; Jafri, Y.Z. Time series models to simulate and forecast hourly averaged wind speed in Quetta, Pakistan. Sol. Energy 1997, 61, 23–32. [Google Scholar]
  71. Cadenas, E.; Rivera, W. Wind speed forecasting in the south coast of Oaxaca, Mexico. Renew. Energy 2007, 32, 2116–2128. [Google Scholar]
  72. Alexiadis, M.; Dokopoulos, P.; Sahsamanoglou, H. Wind speed and power forecasting based on spatial correlation models. IEEE Trans. Energy Convers. 1999, 14, 836–842. [Google Scholar]
  73. Keyhani, A.; Ghasemi-Varnamkhasti, M.; Khanali, M.; Abbaszadeh, R. An assessment of wind energy potential as a power generation source in the capital of Iran, Tehran. Energy 2010, 35, 188–201. [Google Scholar]
  74. Mirhosseini, M.; Sharifi, F.; Sedaghat, A. Assessing the wind energy potential locations in province of Semnan in Iran. Renew. Sustain. Energy Rev. 2011, 15, 449–459. [Google Scholar]
  75. Alamdari, P.; Nematollahi, O.; Mirhosseini, M. Assessment of wind energy in Iran: A review. Renew. Sustain. Energy Rev. 2012, 16, 836–860. [Google Scholar]
  76. Jaramillo, O.; Saldaña, R.; Miranda, U. Wind power potential of Baja California Sur, México. Renew. Energy 2004, 29, 2087–2100. [Google Scholar]
  77. Mathew, S.; Pandey, K.; Kumar, V.A. Analysis of wind regimes for energy estimation. Renew. Energy 2002, 25, 381–399. [Google Scholar]
  78. Weisser, D. A wind energy analysis of Grenada: An estimation using the “Weibull” density function. Renew. Energy 2003, 28, 1803–1812. [Google Scholar]
  79. Chang, T.-J.; Wu, Y.-T.; Hsu, H.-Y.; Chu, C.-R.; Liao, C.vM. Assessment of wind characteristics and wind turbine characteristics in Taiwan. Renew. Energy 2003, 28, 851–871. [Google Scholar]
  80. Akpinar, S.; Akpinar, E.K. Estimation of wind energy potential using finite mixture distribution models. Energy Convers. Manag. 2009, 50, 877–884. [Google Scholar]
  81. Exif. Kutubdia in Bangladesh, 2008. Available online: https://www.flickr.com/photos/rahmanmm/2370901630/meta/ (accessed on 15 April 2014).
  82. Khadem, S.K.; Hussain, M. A pre-feasibility study of wind resources in Kutubdia Island, Bangladesh. Renew. Energy 2006, 31, 2329–2341. [Google Scholar]
  83. Rahaman, S. A Field Report Ongeoenvironmental and Biological Aspects of Sitakunda; Jahangirnagar University: Dhaka, Bangladesh, 2010. [Google Scholar]
  84. Ramírez, P.; Carta, J.A. The use of wind probability distributions derived from the maximum entropy principle in the analysis of wind energy. A case study. Energy Convers. Manag. 2006, 47, 2564–2577. [Google Scholar]
  85. Kantar, Y.M.; Usta, I. Analysis of wind speed distributions: Wind distribution function derived from minimum cross entropy principles as better alternative to Weibull function. Energy Convers. Manag. 2008, 49, 962–973. [Google Scholar]
  86. Weibull, W. A statistical distribution function of wide applicability. J. Appl. Mech. 1951, 18, 293–297. [Google Scholar]
  87. Xiao, Y.; Li, Q.; Li, Z.; Chow, Y.; Li, G. Probability distributions of extreme wind speed and its occurrence interval. Eng. Struct. 2006, 28, 1173–1181. [Google Scholar]
  88. Manwell, J.F.; Mcgowan, J.G.; Rogers, A.L. Wind Energy Explained: Theory, Design and Application; John Wiley & Sons Ltd.: London, UK, 2002; p. 577. [Google Scholar]
  89. Gove, J.H. Moment and maximum likelihood estimators for Weibull distributions under length-and area-biased sampling. Environ. Ecol. Stat. 2003, 10, 455–467. [Google Scholar]
  90. Guure, C.B.; Ibrahim, N.A. Bayesian analysis of the survival function and failure rate of Weibull distribution with censored data. Math. Probl. Eng. 2012, 2012. [Google Scholar] [CrossRef]
  91. Sedefian, L. On the vertical extrapolation of mean wind power density. J. Appl. Meteorol. 1980, 19, 488–493. [Google Scholar]
  92. Justel, A.; Peña, D.; Zamar, R. A multivariate Kolmogorov-Smirnov test of goodness of fit. Stat. Probab. Lett. 1997, 35, 251–259. [Google Scholar]
  93. Watson, G.S. Goodness-of-fit tests on a circle. Biometrika 1961, 48, 109–114. [Google Scholar]
  94. Genest, C.; Rémillard, B.; Beaudoin, D. Goodness-of-fit tests for copulas: A review and a power study. Insur. Math. Econ. 2009, 44, 199–213. [Google Scholar]
  95. Balderrama, J.A.; Masters, F.J.; Gurley, K.R. Peak factor estimation in hurricane surface winds. J. Wind Eng. Ind. Aerodyn. 2012, 102, 1–13. [Google Scholar]
  96. Justus, C.; Mikhail, A. Height variation of wind speed and wind distributions statistics. Geophys. Res. Lett. 1976, 3, 261–264. [Google Scholar]
  97. Justus, C.; Hargraves, W.; Mikhail, A.; Graber, D. Methods for estimating wind speed frequency distributions. J. Appl. Meteorol. 1978, 17, 350–353. [Google Scholar]
  98. Dorvlo, A.S.S. Estimating wind speed distribution. Energy Convers. Manag. 2002, 43, 2311–2318. [Google Scholar]
  99. Lun, I.Y.F.; Lam, J.C. A study of Weibull parameters using long-term wind observations. Renew. Energy 2000, 20, 145–153. [Google Scholar]
  100. Elliott, D.L.; Schwartz, M.N. Wind Energy Potential in the United States. Proceedings of the Project Energy '93: Real Energy Technologies—Environmentally Responsible—Ready for Today, Independence, MO, USA, 21–23 June 1993.
  101. Yu, X.; Qu, H. Wind power in China—Opportunity goes with challenge. Renew. Sustain. Energy Rev. 2010, 14, 2232–2237. [Google Scholar]
  102. Ilinca, A.; Mccarthy, E.; Chaumel, J.-L.; Rétiveau, J.-L. Wind potential assessment of Quebec Province. Renew. Energy 2003, 28, 1881–1897. [Google Scholar]
  103. Zhou, W.; Yang, H.; Fang, Z. Wind power potential and characteristic analysis of the Pearl River Delta region, China. Renew. Energy 2006, 31, 739–753. [Google Scholar]

Nomenclature

P

Total power, W/m2

A

Area, m2

Pe

Practically extractable power, W/m2

k

Dimensionless shape parameter

c

Scale parameter (m/s)

f(v)

Weibull probability density function

F(v)

Cumulative distribution function

C

Constant

L

Likelihood function

COV

Coefficient of variation

vi

Random sample of wind speed central to bin i

n

Number of sample or bin

f(vi)

Weibull frequency for wind speed ranging within bin i

f(v≥0)

Probability for wind speed ≥ 0

Wvi

Observed frequency of the wind speed

v m 3

Mean of the cubic wind speed, m/s

εvi

Error of the approximation

zref

Reference height, m

kh

Weibull shape factor at desired height

ch

Weibull scale factor at desired height, m/s

vh

Wind speed at desired height, m/s

Ph

Power at desired height, W/m2

N

Total number of observations

yi,m

i-th calculated value from measured data

xi,w

i-th calculated value from the Weibull distribution

ziv̄

Mean of i-th calculated value from measured data

R2

Analysis of variance

Γ

Gamma function

σ

Standard deviation of wind speed, m/s

ρ

Air density, kg/m3

θ

Unknown parameter for maximum likelihood function

η

Power law coefficient

χ2

Chi-square error

Figure 1. Hourly mean wind speed at 20 m and 30 m height for the selected sites.
Figure 1. Hourly mean wind speed at 20 m and 30 m height for the selected sites.
Energies 07 03056f1 1024
Figure 2. Hourly mean COV of turbulence at 20 m and 30 m height.
Figure 2. Hourly mean COV of turbulence at 20 m and 30 m height.
Energies 07 03056f2 1024
Table 1. Site name, location for wind speed measuring at 20 m and 30 m height above surface.
Table 1. Site name, location for wind speed measuring at 20 m and 30 m height above surface.
Wind station nameSite nameLatitude (°) NLongitude (°) EHeight above the sea levelReference
Station-IKuakata21°54.76′90°08.24′3 m[81]
Station-IIKutubdia21°54.71′91°52.43′0–4 m[82]
Station-IIISitakunda22°35.68′91°42.52′9 m[83]
Table 2. Monthly mean wind speed (m/s) data for selected wind stations at 20 m and 30 m height. COV: coefficient of variation.
Table 2. Monthly mean wind speed (m/s) data for selected wind stations at 20 m and 30 m height. COV: coefficient of variation.
WindMonth of the yearYearly mean

JanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
Station-IH = 30 mN (h)7446727447207447207447447207447207448760
vmean (m/s)3.113.573.544.984.975.075.745.974.952.572.983.224.22
vmax (m/s)5.085.795.528.268.287.559.469.299.523.55.777.457.12
vmin (m/s)1.662.252.442.522.322.483.232.341.971.51.311.262.11
σ (m/s)0.780.860.821.541.681.561.612.132.140.540.761.391.32
COV (%)25.0824.0923.1630.9233.830.7728.0535.6843.2321.0125.543.1730.37

H = 20 mN (h)7446727447207447207447447207447207448760
vmean (m/s)2.152.622.653.763.813.824.444.753.811.762.122.553.19
vmax (m/s)3.654.144.296.036.255.588.158.217.332.774.97.25.71
vmin (m/s)1.221.681.781.881.691.772.711.721.530.990.90.971.57
σ (m/s)0.520.690.661.11.241.171.371.831.730.450.651.41.07
COV (%)24.1926.3424.9129.2632.5530.6330.8638.5345.4125.5730.6654.932.82

Station-IIH = 30 mN (h)7446727447207447207447447207447207448760
vmean (m/s)3.112.83.153.734.244.535.424.743.392.252.962.653.58
vmax (m/s)5.265.694.816.366.957.357.587.587.483.314.154.375.91
vmin (m/s)1.641.291.672.222.352.062.821.991.291.251.71.731.83
σ (m/s)0.951.330.691.161.391.521.171.561.740.590.50.641.1
COV (%)30.5547.521.931.132.7833.5521.5932.9151.3326.2216.8924.1530.87

H = 20 mN (h)7446727447207447207447447207447207448760
vmean (m/s)1.691.791.92.733.33.814.633.952.581.261.451.112.52
vmax (m/s)3.514.473.75.16.316.317.266.797.132.412.682.424.84
vmin (m/s)0.70.640.761.491.391.382.241.610.870.560.590.51.06
σ (m/s)0.711.250.611.081.421.411.181.451.840.460.430.441.02
COV (%)42.0169.8332.1139.5643.0337.0125.4936.7171.3236.5129.6639.6441.91

Station-IIIH = 30 mN (h)7446727447207447207447447207447207448760
vmean (m/s)2.723.173.364.354.464.815.365.1510.092.252.72.324.23
vmax (m/s)4.166.215.147.367.558.167.5913.7639.23.483.693.679.16
vmin (m/s)1.161.241.962.092.632.142.912.311.491.231.861.111.84
σ (m/s)0.841.430.821.421.41.791.142.5813.070.480.430.712.18
COV (%)30.8845.1124.432.6431.3937.2121.2750.1129.5321.3315.9330.639.2

H = 20 mN (h)7446727447207447207447447207447207448760
vmean (m/s)1.962.472.643.723.844.264.794.875.971.431.821.663.29
vmax (m/s)3.365.614.356.626.867.346.8914.3518.582.62.672.76.83
vmin (m/s)0.660.791.331.652.011.562.311.731.280.760.930.861.32
σ (m/s)0.731.450.821.421.421.791.092.885.330.420.360.521.52
COV (%)37.2458.731.0638.1736.9842.0222.7659.1489.2829.3719.7831.3341.32
Table 3. Frequency distribution and cumulative frequency distribution for the selected sites in June.
Table 3. Frequency distribution and cumulative frequency distribution for the selected sites in June.
Wind speed (m/s)Station-IStation-IIStation-III



Frequency (%)Cumulative frequency (%)Frequency (%)Cumulative frequency (%)Frequency (%)Cumulative frequency (%)
0–10.4030.4030.2690.2690.1340.134
1–22.1512.5541.8822.1510.8060.94
2–35.2427.7963.6295.785.6456.585
3–49.00516.80114.24720.02717.33923.924
4–517.74234.54319.08639.11324.05947.983
5–620.4354.97324.32863.44118.41466.397
6–718.01172.98418.2881.72117.74284.139
7–812.63485.61811.55993.288.19992.338
8–98.60294.225.24298.5223.89896.236
9–104.03298.2521.34499.8662.01698.252
10–111.61399.865099.8661.34499.596
11–120.1341000.1341000.26999.865
12–13----0.134100
Table 4. The Weibull distribution analysis for Station-I at 20 m height. GM: graphical method; MOM: method of moments; STDM: standard deviation method; MLM: maximum likelihood method; PDM: power density method; MMLM: Modified maximum likelihood method; and EEM: equivalent energy method.
Table 4. The Weibull distribution analysis for Station-I at 20 m height. GM: graphical method; MOM: method of moments; STDM: standard deviation method; MLM: maximum likelihood method; PDM: power density method; MMLM: Modified maximum likelihood method; and EEM: equivalent energy method.
Statistical methodsThe Weibull parametersStatistical test efficiency


k (-)c (m/s)RPE (%)RMSEMPE (%)MAPE (%)χ2R2
GM3.633.550.07850.05000.82834.25710.14320.9577
MOM3.583.54−0.05230.0408−0.03940.03940.00010.9999
STDM3.994.4626.67360.921927.982627.9822.25930.1542
MLM3.53.550.00000.0000−0.04380.37000.00050.9998
PDM3.113.53−0.78450.1581−0.66490.66480.02850.9917
MMLM3.543.550.07840.05000.07230.18850.00020.9999
EEM3.373.54−0.28760.09570.15072.46870.04510.9866

Note: Dimensionless shape factor (k) and scale factor (c) in m/s, relative percentage error (RPE), root mean square error (RMSE), mean percentage error (MPE), mean absolute percentage error (MAPE), chi-square error (χ2), analysis of variance or efficiency of the method (R2).

Table 5. The Weibull distribution analysis for Station-I at 30 m height.
Table 5. The Weibull distribution analysis for Station-I at 30 m height.
Statistical methodsThe Weibull parametersStatistical test efficiency


k (-)c (m/s)RPE (%)RMSEMPE (%)MAPE (%)χ2R2
GM2.714.790.88800.19300.62501.70400.03500.9900
MOM3.874.68−0.55300.1530−0.54200.58900.00230.9990
STDM2.644.781.73700.27102.33103.71800.23800.9320
MLM3.784.680.00000.0000−0.05600.34400.00100.9997
PDM2.584.66−1.20400.2250−1.18701.23900.01500.9960
MMLM3.184.67−0.80900.1850−0.92100.95400.00640.9980
EEM3.264.730.49300.14400.62901.62200.02900.9920
Table 6. The Weibull distribution analysis for Station-II at 20 m height.
Table 6. The Weibull distribution analysis for Station-II at 20 m height.
Statistical methodsThe Weibull parametersStatistical test efficiency


k (-)c (m/s)RPE (%)RMSEMPE (%)MAPE (%)χ2R2
GM3.112.70−3.87420.3120−4.15464.56440.13360.9787
MOM2.822.810.06620.04080.04220.04220.00020.9999
STDM3.43.7434.20530.927840.003540.0032.84530.3654
MLM2.842.93.04630.27686.78237.08910.30630.9475
PDM2.652.820.13240.05770.11700.26720.00030.9999
MMLM2.842.851.58940.20003.48233.66710.07950.9865
EEM2.882.76−1.82120.2140−1.93772.19660.03370.9945
Table 7. The Weibull distribution analysis for Station-II at 30 m height.
Table 7. The Weibull distribution analysis for Station-II at 30 m height.
Statistical methodsThe Weibull parametersStatistical test efficiency


k (-)c (m/s)RPE (%)RMSEMPE (%)MAPE (%)χ2R2
GM4.173.88−1.81500.2550−1.51982.74440.06100.9790
MOM3.993.94−0.58200.1440−0.59300.64900.00740.9974
STDM4.474.8322.64400.900423.863023.8631.90200.1780
MLM3.923.96−0.09300.0580−0.09200.52200.00200.9990
PDM3.253.99−0.25600.0960−0.25600.25600.00040.9998
MMLM3.593.98−6.51600.4830−8.31608.52702.34200.2280
EEM4.234.3911.19400.633011.771011.7710.51800.7970
Table 8. The Weibull distribution analysis for Station-III at 20 m height.
Table 8. The Weibull distribution analysis for Station-III at 20 m height.
Statistical methodsThe Weibull parametersStatistical test efficiency


k (-)c (m/s)RPE (%)RMSEMPE (%)MAPE (%)χ2R2
GM3.043.660.02540.0289−0.88254.49010.15680.979
MOM3.143.62−0.25360.0913−0.17110.17110.00310.9996
STDM3.434.9535.12551.074334.535334.5355.70010.0331
MLM3.133.650.17750.07640.08260.44530.00110.9998
PDM2.773.65−0.05070.0408−0.03110.03110.00010.9999
MMLM3.143.64−1.16660.1958−0.82051.40340.09240.9877
EEM2.913.660.00000.0000−0.4492.27020.0410.9945
Table 9. The Weibull distribution analysis for Station-III at 30m height.
Table 9. The Weibull distribution analysis for Station-III at 30m height.
Statistical methodsThe Weibull parametersStatistical test efficiency


k (-)c (m/s)RPE (%)RMSEMPE (%)MAPE (%)χ2R2
GM2.14.71−0.59100.1580−0.52502.86300.06400.9950
MOM3.764.450.07900.05800.13200.16500.00100.9999
STDM2.33.73−21.5400.9540−11.41013.92020.6800.3470
MLM3.664.5−0.47300.1414−0.30950.53400.01400.9990
PDM2.324.741.49800.2520−0.07303.46300.30800.9740
MMLM2.994.620.15800.0820−0.45101.97600.05800.9950
EEM3.034.09−11.8450.7080−6.05107.43105.96900.5632
Table 10. Ranking of the methods by statistical test results.
Table 10. Ranking of the methods by statistical test results.
Statistical methodsStation-IStation-IIStation-IIIDiscussion



20 m30 m20 m30 m20 m30 m
GMSixthSixthFifthFourthSixthFourth-
MOMFirstSecondFirstThirdThirdFirstThe First choice
STDMSeventhSeventhSeventhSeventhSeventhSeventh-
MLMThirdFirstSixthSecondSecondSecondThe second choice
PDMFourthFourthSecondFirstFirstFifthThe third choice
MMLMSecondThirdFourthSixthFifthThird-
EEMFifthFifthThirdFifthFourthSixth-
Table 11. Monthly mean wind speed, power law coefficient (η), the Weibull shape factor (k) and scale factor (c) at 20 m, 30 m and 50 m height.
Table 11. Monthly mean wind speed, power law coefficient (η), the Weibull shape factor (k) and scale factor (c) at 20 m, 30 m and 50 m height.
MonthPower law coefficient ηMeasured valueExtrapolated value


At 20 m heightAt 30 m heightAt 50 m height

v20k20c20v30k30c30v50k50c50
January0.2622.154.332.353.114.423.403.554.633.89
February0.252.624.182.893.574.413.914.054.624.44
March0.2512.654.122.913.544.413.874.024.624.4
April0.2193.763.864.164.983.585.535.573.756.18
May0.2193.813.454.254.973.305.555.583.466.21
June0.2183.823.954.235.073.845.635.694.026.29
July0.2074.443.414.935.743.826.346.394.007.05
August0.2034.752.915.345.973.176.696.663.327.42
September0.2183.812.44.314.952.525.605.562.646.26
October0.281.764.361.942.575.652.792.985.923.22
November0.2662.123.062.342.983.703.263.373.873.73
December0.2562.551.992.893.222.493.633.682.614.14
Table 12. Monthly mean power based on measured data and calculated data by statistical methods.
Table 12. Monthly mean power based on measured data and calculated data by statistical methods.
MonthPower (W/m2)MOMMLM


k (-)c (m/s)P (W/m2)RMSEχ2k (-)c (m/s)P (W/m2)RMSEχ2
January18.754.503.4121.841.75780.99044.423.421.671.70880.9724
February27.874.713.9032.632.18174.413.9133.242.3173
March26.714.923.8631.402.16564.413.8731.962.2913
April73.113.585.5397.574.94573.585.5397.574.9457
May71.473.255.55101.465.47633.35.55101.045.4378
June75.873.65.63102.855.19423.845.63101.215.0339
July110.113.986.34143.455.77413.826.34144.75.8813
August123.883.066.68181.117.56513.176.69179.367.4485
September70.612.485.58118.126.89272.525.6117.256.8293
October10.055.472.7811.701.28455.652.7911.81.3229
November16.214.423.2719.381.78043.73.2619.821.9000
December20.812.483.6332.463.41322.493.6332.373.4000
Table 13. Seasonal mean wind speed in (m/s), the Weibull shape factor (-), scale factor (m/s) and power density (W/m2) at 30 m and 50 m height.
Table 13. Seasonal mean wind speed in (m/s), the Weibull shape factor (-), scale factor (m/s) and power density (W/m2) at 30 m and 50 m height.
SeasonsηMeasured dataExtrapolated data


At 20 m heightAt 30 m heightAt 50 m height

v20k20c20P20v30k30c30P30v50k50c50P50
Winter0.243.193.53.5526.214.223.784.759.14.763.965.2783.6
Spring0.243.014.053.3220.24.034.134.448.24.554.335.0168.7
Summer0.214.023.64.4748.645.263.655.8108.245.893.836.52157.2
Autumn0.233.443.223.8634.14.53.785.172.35.073.965.63100.67
Table 14. Wind power classification.
Table 14. Wind power classification.
Power classPotentialPower density and wind speed at 10 m (33 ft)Power density and wind speed at 30 m (98 ft)Power density and wind speed at 50 m (164 ft)



Power (W/m2)Speed (m/s)Power (W/m2)Speed (m/s)Power (W/m2)Speed (m/s)
Class-IPoorP10 ≤ 100≤4.4P30 ≤ 160≤5.1P50 ≤ 200≤5.6
Class-IIMarginalP10 ≤ 150≤5.1P30 ≤ 240≤6.0P50 ≤ 300≤6.0
Class-IIIModerateP10 ≤ 200≤5.6P30 ≤ 320≤6.5P50 ≤ 400≤7.0
Class-IVGoodP10 ≤ 250≤6.0P30 ≤ 400≤7.0P50 ≤ 500≤7.5
Class-VVery goodP10 ≤ 300≤6.4P30 ≤ 480≤7.5P50 ≤ 600≤8.0
Class-VIExcellentP10 ≤ 400≤7.0P30 ≤ 640≤8.2P50 ≤ 800≤8.8
Class-VIIExcellentP10 ≤ 1000≤9.4P30 ≤ 1600≤11.0P50 ≤ 2000≤11.9
Table 15. The most frequent wind velocity (VFmax), velocity contributing the maximum energy (VEmax), energy density (ED) and total energy intensity (EDT).
Table 15. The most frequent wind velocity (VFmax), velocity contributing the maximum energy (VEmax), energy density (ED) and total energy intensity (EDT).
SeasonsVFmax (m/s)VEmax (m/s)Energy density, ED (kW/m2)Total energy intensity, EDT (kW/m2)
Winter4.315.240.0586129.408
Spring4.154.880.0492105.082
Summer5.356.580.1148253.531
Autumn4.645.630.0728160.667

Yearly mean4.315.240.0586513.412
Table 16. Energy generated by the wind turbine and the capacity factor by numerical integration. The turbine rated power is 20 kW or 0.02 MW, rated speed is 8 m/s and hub height at 30 m.
Table 16. Energy generated by the wind turbine and the capacity factor by numerical integration. The turbine rated power is 20 kW or 0.02 MW, rated speed is 8 m/s and hub height at 30 m.
SeasonsNumerical integration

Trapezoidal sumsSimpson's rule

Capacity factor, CFTotal energy output, ET (MW·h)Capacity factor, CFTotal energy output, ET, (MW·h)
Winter0.12385.46890.12385.4689
Spring0.09203.92830.09203.9283
Summer0.301013.29370.301013.2938
Autumn0.16987.50050.16987.5005

Yearly mean0.123821.69730.123821.6974

Share and Cite

MDPI and ACS Style

Azad, A.K.; Rasul, M.G.; Yusaf, T. Statistical Diagnosis of the Best Weibull Methods for Wind Power Assessment for Agricultural Applications. Energies 2014, 7, 3056-3085. https://doi.org/10.3390/en7053056

AMA Style

Azad AK, Rasul MG, Yusaf T. Statistical Diagnosis of the Best Weibull Methods for Wind Power Assessment for Agricultural Applications. Energies. 2014; 7(5):3056-3085. https://doi.org/10.3390/en7053056

Chicago/Turabian Style

Azad, Abul Kalam, Mohammad Golam Rasul, and Talal Yusaf. 2014. "Statistical Diagnosis of the Best Weibull Methods for Wind Power Assessment for Agricultural Applications" Energies 7, no. 5: 3056-3085. https://doi.org/10.3390/en7053056

APA Style

Azad, A. K., Rasul, M. G., & Yusaf, T. (2014). Statistical Diagnosis of the Best Weibull Methods for Wind Power Assessment for Agricultural Applications. Energies, 7(5), 3056-3085. https://doi.org/10.3390/en7053056

Article Metrics

Back to TopTop