# Generalized Extreme Value Distribution Models for the Assessment of Seasonal Wind Energy Potential of Debuncha, Cameroon

Authors: Nkongho Ayuketang Arreyndip; Ebobenow Joseph
Journal: Journal of Renewable Energy
 (2016)

---

## Abstract

The method of generalized extreme value family of distributions (Weibull, Gumbel, and Frechet) is employed for the first time to assess the wind energy potential of Debuncha, South-West Cameroon, and to study the variation of energy over the seasons on this site. The 29-year (1983–2013) average daily wind speed data over Debuncha due to missing values in the years 1992 and 1994 is gotten from NASA satellite data through the RETScreen software tool provided by CANMET Canada. The data is partitioned into min-monthly, mean-monthly, and max-monthly data and fitted using maximum likelihood method to the two-parameter Weibull, Gumbel, and Frechet distributions for the purpose of determining the best fit to be used for assessing the wind energy potential on this site. The respective shape and scale parameters are estimated. By making use of theP values of the Kolmogorov-Smirnov statistic (K-S) and the standard error (s.e) analysis, the results show that the Frechet distribution best fits the min-monthly, mean-monthly, and max-monthly data compared to the Weibull and Gumbel distributions. Wind speed distributions and wind power densities of both the wet and dry seasons are compared. The results show that the wind power density of the wet season was higher than in the dry season. The wind speeds at this site seem quite low; maximum wind speeds are listed as between 3.1 and 4.2 m/s, which is below the cut-in wind speed of many modern turbines (6–10 m/s). However, we recommend the installation of low cut-in wind turbines like the Savonius or Aircon (10 KW) for stand-alone low energy need.


---

## Body

## 1. Introduction

Cameroon which is also in the list of fast growing economies in Africa has experienced a very fast growth rate in the past two decades with the Government investing heavily on industrialization. Statistics has shown that the number of heavy and light industries has doubled since the early 90s. The population statistics has significantly increased from 15 million inhabitants in the 80s to about 21 million inhabitants. These are the two main concepts that plunge an economy energy sector into serious challenges like energy shortages and distribution. The concept of energy shortage is a global issue as many world economies now have embarked on alternative renewable energy sources to meet their energy demand. Also faced with the crisis of global warming, the Paris summit with Cameroon represented by its head of state was to tackle this issue with most world economies accepting the 1.5°C rise in temperature and shifting their energy production from fossil fuels combustion to renewable energy enhancement. Cameroon is known for its enormous energy potentials. There is abundance of heavy water falls in the south to harness hydroelectric power, high solar intensities especially the northern regions for solar power enhancement, and wind power in its coastal cities for wind energy enhancement. These sectors are underdeveloped and the country solely relies on hydroelectric power to meet its energy demand with most of these rivers seasonal. Therefore it is worthwhile looking for other alternative sources of energy like wind and solar energy. To successfully harness wind power at a site, a wind energy assessment needs to be carried out which usually results in analysing the wind speed data at the particular site and also studying the topographic nature of the site.Many researchers have applied statistical techniques to assess the wind energy potential of a particular site or sites like the commonly used Weibull statistics and others like the lognormal and Rayleigh distributions. For example, Fagbenle et al. [1] applied the Weibull statistics and other statistical methods to assess the wind energy potentials of two sites in north-east Nigeria and found that Maiduguri was the better of the sites as compared to Potiskum in terms of monthly and seasonal variation of mean wind speed, but they both can be suitable for stand-alone and medium scale wind power generation. Oyedepo et al. [2] also applied the Weibull statistics to analyse wind speed data and wind energy potential in three selected locations (Enugu, Owerri, and Onitsha) in south-east Nigeria, where they found that Enugu was the best site with high wind energy potential. Nikolai Nawri et al. applied the downscaling simulations technique performed with the Weather Research and Forecasting (WRF) model to determine the large-scale wind energy potential of Iceland with local wind speed distributions represented by Weibull statistics. They found that, in addition to seasonal and spatial variability, differences in average wind speed and power density also exist for different wind directions. Kollu et al. [3] applied the three mixture probability density functions: Weibull-extreme value distribution (GEV), Weibull-lognormal, and GEV-lognormal which were not tried before to describe (model) wind speed characteristics. By making use of statistical parameters such as maximum error in the Kolmogorov-Smirnov test, root mean square error, Chi-square error, coefficient of determination, and power density error as judgement criteria to assess the fitness of the probability density functions, they found that the Weibull-GEV PDF is able to describe unimodal as well as bimodal wind distributions accurately whereas GEV-lognormal PDF is able to describe familiar bell-shaped unimodal distribution well. They also found that mixture probability functions are better alternatives to conventional Weibull, two-component mixture Weibull, gamma, and lognormal PDFs to describe wind speed characteristics.Most of the results in this area show that the Weibull distribution is the best distribution to mathematically model the wind energy potential of a particular site. So in this paper, our focus and analysis are based on the fact that, being a member of the generalized extreme value class of distributions also known as generalized extreme value type III distribution, it is worthwhile also to study the other class of distributions (Gumbel (type I), Frechet (type II)) in the same family to know if there is another distribution that can best be used to describe the wind energy potential of a site.


### 1.1. Geophysical Description of Debuncha

Debuncha, being known as the rainiest place in central and west Africa, is a small Cameroonian village on the western coast of the city of Limbe at the windward side of mount Cameroon. Some wind energy experts have assessed this site and it has been found suitable for the installation of an on-shore wind turbine because of its topographic structure and also because of its accessibility. So our main aim in this work is to select in the family of generalized extreme value class of distributions which distribution is best suited to use for wind power assessment and modelling for the feature installation of an on-shore wind turbine, study the conditions under which the two-parameter Weibull distribution which has so far been termed the best statistical distribution for wind power modelling and assessment can best be used, and finally to study the wind power variation over the seasons on this site.The paper is arranged in three sections; in Section2, we present the mathematical models of the underlying theory and statistically analyse the data set. Section 3 is dedicated to discussing the results of our findings and we end with a conclusion in Section 4.


## 1.1. Geophysical Description of Debuncha

Debuncha, being known as the rainiest place in central and west Africa, is a small Cameroonian village on the western coast of the city of Limbe at the windward side of mount Cameroon. Some wind energy experts have assessed this site and it has been found suitable for the installation of an on-shore wind turbine because of its topographic structure and also because of its accessibility. So our main aim in this work is to select in the family of generalized extreme value class of distributions which distribution is best suited to use for wind power assessment and modelling for the feature installation of an on-shore wind turbine, study the conditions under which the two-parameter Weibull distribution which has so far been termed the best statistical distribution for wind power modelling and assessment can best be used, and finally to study the wind power variation over the seasons on this site.The paper is arranged in three sections; in Section2, we present the mathematical models of the underlying theory and statistically analyse the data set. Section 3 is dedicated to discussing the results of our findings and we end with a conclusion in Section 4.


## 2. Data and Method

### 2.1. Raw Data Analysis

The data provided to us for analysis is NASA satellite daily average wind data over Debuncha on the western coast of Limbe, Cameroon, with longitude 9.2149, latitude 4.0242, and elevation of about 36 m above sea level. This data contains 29 (1983–2013) years of mean daily data with missing values in the years 1992 and 1994. Hence these years were considered as outliers. The data was sorted out into monthly means (mean-monthly), monthly minimums (min-monthly), and monthly maximums (max-monthly) for the 29 years.The second analysis deals with comparing the wind energy potentials of the two main seasons that exist in Cameroon which are the wet and dry seasons. The wet seasons usually run from May to September, while the dry season usually runs from November to March. In April and October there usually occur the short Springs and Autumn seasons, respectively. So they were not included in the analysis. The data is again partitioned into the five months’ dry season (November to March) and five months’ wet season (May to September) for analysis over the 29 years’ period. The yearly maximum wind speeds over this area for the 29 years’ period are shown in Figure1 with missing data in 1992 and 1994. The min-monthly, mean-monthly, and max-monthly wind distributions are also shown in Figures 2, 3, and 4.Figure 1
Annual maximum wind speed for Debuncha.Figure 2
Plot of monthly minimums wind speed distribution over Debuncha for 31 years’ period.Figure 3
Plot of monthly average wind speed distribution over Debuncha for 31 years’ period.Figure 4
Plot of monthly maximums wind speed distribution over Debuncha for 31 years’ period.


### 2.2. Mathematical Analysis

We will start this section by reviewing the cumulative distribution function of the three-parameter generalized extreme value (GEV) distribution given by(1)Fx=exp⁡-1+ξx-μσ-1ξ.Define, forx:1+ξx-μ/σ>0, -∞<μ<∞, σ>0, and -∞<ξ<∞, where μ is the location parameter, ξ is the shape parameter, and σ is the scale parameter [4, 5]. There are three classes of the generalized extreme value family of distribution. Their differences depend only on the value of the shape parameter ξ.We have the Gumbel distribution withξ=0 with probability density function, cumulative probability distribution, and quantiles given by(2)fx=1σexp⁡zexp⁡-exp⁡zFx=1-exp⁡-exp⁡zQx=μ-σlog⁡-log⁡x,where (3)z≡x-μσ.Forξ<0 we have the Weibull distribution with probability density function, cumulative probability distribution, and quantile distribution given by(4)fx=ξσxσξ-1exp⁡-xσξ.Fx=1-exp-xσξQx=σ-log⁡1-x1/ξForξ>0, we have the Freshet class of distribution with cumulative probability distribution given by(5)fx=ξσσxξ+1exp⁡σxξFx=exp⁡-σxξQx=σ-log⁡x1/ξ.Now using the Weibull distribution, the mean values of the wind speedxm and standard deviation ν defined in terms of the Weibull parameters ξ and σ are given as(6)xm=σΓ1+1ξ,where Γ is the gamma function, and(7)ν=σ2Γ1+2ξ-Γ1+1ξ2.


### 2.3. Mean Power and Energy Densities

When assessing the wind energy potential of a site, there are two wind speeds that are of interest and must therefore be taken seriously into consideration. These are the most probable wind speedvmp and the wind speed carrying maximum energy vEmax⁡. They are given by the expressions [6, 7](8)vmp=σξ-1ξ1/ξ,vEmax⁡=σξ+2ξ1/ξ.Theoretically, the mean power density is proportional to the cube of the mean velocity given by(9)PD=PvmA=12ρvm3.It can also be calculated from the Weibull probability density function given by(10)PD=PvA=12ρσ3Γ1+3ξand the energy density is given by(11)ED=12ρσ3Γ1+3ξT,where PD is the power density, ED the energy density, Pv wind power, A cross-sectional area of rotor, ρ the air density assumed here to take the value 1.225, c scale parameter, ξ shape parameter, T period, and Γ the gamma function [1, 6, 7].The data set over the operating conditions (min-monthly, mean-monthly, and max-monthly) is fitted into the GEV class of distribution in R using the R-package “fitdistrplus” using the maximum goodness of fit and the parameters are estimated using the method of maximum likelihood which is given by [8, 9](12)ξ=∑i=1nviξln⁡vi∑i=1nviξ-∑i=1nlnvin-1,where vi is the wind speed in time step i and n is the number of data points.The scale parameter is given by(13)σ=1n∑i=1nviξ1/ξ.


## 2.1. Raw Data Analysis

The data provided to us for analysis is NASA satellite daily average wind data over Debuncha on the western coast of Limbe, Cameroon, with longitude 9.2149, latitude 4.0242, and elevation of about 36 m above sea level. This data contains 29 (1983–2013) years of mean daily data with missing values in the years 1992 and 1994. Hence these years were considered as outliers. The data was sorted out into monthly means (mean-monthly), monthly minimums (min-monthly), and monthly maximums (max-monthly) for the 29 years.The second analysis deals with comparing the wind energy potentials of the two main seasons that exist in Cameroon which are the wet and dry seasons. The wet seasons usually run from May to September, while the dry season usually runs from November to March. In April and October there usually occur the short Springs and Autumn seasons, respectively. So they were not included in the analysis. The data is again partitioned into the five months’ dry season (November to March) and five months’ wet season (May to September) for analysis over the 29 years’ period. The yearly maximum wind speeds over this area for the 29 years’ period are shown in Figure1 with missing data in 1992 and 1994. The min-monthly, mean-monthly, and max-monthly wind distributions are also shown in Figures 2, 3, and 4.Figure 1
Annual maximum wind speed for Debuncha.Figure 2
Plot of monthly minimums wind speed distribution over Debuncha for 31 years’ period.Figure 3
Plot of monthly average wind speed distribution over Debuncha for 31 years’ period.Figure 4
Plot of monthly maximums wind speed distribution over Debuncha for 31 years’ period.


## 2.2. Mathematical Analysis

We will start this section by reviewing the cumulative distribution function of the three-parameter generalized extreme value (GEV) distribution given by(1)Fx=exp⁡-1+ξx-μσ-1ξ.Define, forx:1+ξx-μ/σ>0, -∞<μ<∞, σ>0, and -∞<ξ<∞, where μ is the location parameter, ξ is the shape parameter, and σ is the scale parameter [4, 5]. There are three classes of the generalized extreme value family of distribution. Their differences depend only on the value of the shape parameter ξ.We have the Gumbel distribution withξ=0 with probability density function, cumulative probability distribution, and quantiles given by(2)fx=1σexp⁡zexp⁡-exp⁡zFx=1-exp⁡-exp⁡zQx=μ-σlog⁡-log⁡x,where (3)z≡x-μσ.Forξ<0 we have the Weibull distribution with probability density function, cumulative probability distribution, and quantile distribution given by(4)fx=ξσxσξ-1exp⁡-xσξ.Fx=1-exp-xσξQx=σ-log⁡1-x1/ξForξ>0, we have the Freshet class of distribution with cumulative probability distribution given by(5)fx=ξσσxξ+1exp⁡σxξFx=exp⁡-σxξQx=σ-log⁡x1/ξ.Now using the Weibull distribution, the mean values of the wind speedxm and standard deviation ν defined in terms of the Weibull parameters ξ and σ are given as(6)xm=σΓ1+1ξ,where Γ is the gamma function, and(7)ν=σ2Γ1+2ξ-Γ1+1ξ2.


## 2.3. Mean Power and Energy Densities

When assessing the wind energy potential of a site, there are two wind speeds that are of interest and must therefore be taken seriously into consideration. These are the most probable wind speedvmp and the wind speed carrying maximum energy vEmax⁡. They are given by the expressions [6, 7](8)vmp=σξ-1ξ1/ξ,vEmax⁡=σξ+2ξ1/ξ.Theoretically, the mean power density is proportional to the cube of the mean velocity given by(9)PD=PvmA=12ρvm3.It can also be calculated from the Weibull probability density function given by(10)PD=PvA=12ρσ3Γ1+3ξand the energy density is given by(11)ED=12ρσ3Γ1+3ξT,where PD is the power density, ED the energy density, Pv wind power, A cross-sectional area of rotor, ρ the air density assumed here to take the value 1.225, c scale parameter, ξ shape parameter, T period, and Γ the gamma function [1, 6, 7].The data set over the operating conditions (min-monthly, mean-monthly, and max-monthly) is fitted into the GEV class of distribution in R using the R-package “fitdistrplus” using the maximum goodness of fit and the parameters are estimated using the method of maximum likelihood which is given by [8, 9](12)ξ=∑i=1nviξln⁡vi∑i=1nviξ-∑i=1nlnvin-1,where vi is the wind speed in time step i and n is the number of data points.The scale parameter is given by(13)σ=1n∑i=1nviξ1/ξ.


## 3. Results and Discussion

The results of the method applied here are presented in Figures2–9 and Table 1.Table 1
Data summary statistics and parameter estimates.
Distributions Min wind speed data Average wind speed data Max wind speed data Weibull k 2.76 3.13 3.29 s.e 0.343 0.12 0.125 c (m/s) 1.62 2.56 3.63 s.e 1.00 0.05 0.063 P values of K-S 0.175 0.182 0.172 Gumbel k 1.22 1.98 2.833 s.e 0.019 0.026 0.037 c (m/s) 0.35 0.47 0.665 s.e 0.016 0.021 0.03 P values of K-S 0.149 0.108 0.121 Frechet k 4.88 4.88 4.73 s.e 0.22 0.23 0.20 c (m/s) 1.93 1.93 2.76 s.e 0.022 0.022 0.033 P values of K-S 0.109 0.083 0.074 Descriptive statistics Mean (m/s) 1.44 2.29 3.26 s.d 0.54 0.74 1.0 Min (m/s) 0.7 1.4 1.9 Max (m/s) 3.5 4.62 6.4 Skew 1.39 1.36 1.26 Kurtosis 1.41 0.9 0.79 s.e 0.03 0.04 0.05Figure 5
Mean wind speed distributions over the dry and wet seasons.Figure 6
Mean wind power densities of both the dry and wet seasons.Figure 7
Monthly minimum wind speeds fitted to the GEV class of distributions.Figure 8
Monthly average wind speeds fitted to the GEV class of distributions.Figure 9
Monthly maximum wind speeds fitted to the GEV class of distributions.The parameters of the Weibull, Gumbel, and Freshet distributions are estimated using the method of maximum likelihood and are displayed in Table1.Figure1 shows the yearly maximum wind speed distributions for the 29 years’ period. We see that there is an increase in the wind speed from where it usually fluctuates from 3.1 m/s to 4.2 m/s from 1983 to 2006 to about 6.0 m/s starting from the year 2007 to 2013. This might be as a result of the changes in the pressure gradient and local weather conditions of the place.Figures2, 3, and 4 show that one of the most windy months is August which is usually at the heart of the wet season, and we also note from these figures that the month with the least wind speed is the month of December which is in the dry season and the highest year with a high average recorded wind speed is 2013.Figures5 and 6 display the mean wind speeds and wind power densities of the dry and wet seasons. Considering the mean data, the mean and maximum wind velocities in the dry season were calculated to be 2.61 m/s and 3.22 m/s, respectively, while those in the wet season are 3.36 m/s and 4.17 m/s, respectively. Hence we see that, during the wet season, there is a possibility of generating higher power from a wind turbine on this area than in the dry season.From Figures7, 8, and 9, we fit the monthly minimum, monthly average, and monthly maximums wind speed data into the generalized extreme value class of distributions. The parameter was estimated using the method of maximum likelihood and presented in Table 1. Taking into consideration the P values of the Kolmogorov-Smirnov test statistic (K-S) and the standard error (s.e) analysis, we see that the Frechet distribution best fits with the data in all operating conditions (min, mean, and max) followed by the Gumbel and finally the Weibull distribution. Here, we also see that the density,QQ, and cdf plots throughout the data ranges best fit with the Frechet distribution than the Weibull and Gumbel distributions. These figures also show that the monthly average wind speeds best fit the Weibull distribution as compared to the monthly minimums and the monthly maximums. This can also be confirmed by its lowest standard error and the P values.Table1 also presents an overview of the descriptive statistics of the data set used. Here, we see that average monthly wind speeds over this site range from 1.4 m/s to 4.62 m/s. Monthly average power density ranges from 5 W/m2 to 60 W/m2. These values are generally not encouraging for the installation of modern wind turbines.


## 4. Conclusion

To conclude, we have partitioned a 29 years’ NASA satellite daily average wind speed data gotten through the RETScreen software tool provided by CANMET Canada into monthly minimum, monthly average, and monthly maximum wind speeds and consequently fitted into the generalized extreme value class of distributions (Weibull, Frechet, and Gumbel) for the purpose of determining the best distribution for use to model the wind energy potential of Debuncha, Cameroon, and also to determine the condition under which the Weibull distribution which is so far the most used and recommended distribution for wind data analysis can be applied to the data. Using the Kolmogorov-Smirnov statistic and standard error analysis, the results here show that the Frechet distribution best fits the data in all operating conditions than the Gumbel and Weibull distributions. We have also seen that the monthly average wind speeds best fit the Weibull distribution as compared to the monthly minimums and monthly maximums. Hence, we conclude that when doing a Weibull analysis of the data on this site, the monthly average wind speeds should be used. The wind speed distributions of the dry and wet seasons together with the power densities over these seasons were analysed where we have seen that there is a greater wind speed and hence a higher power density in the wet season as compared to the dry season. The time series analysis of the yearly maximum wind speed data over the 29 years’ period as shown in Figure1 shows that there is an increasing trend in the wind speed for Debuncha from a maximum of 3.1 m/s to 4.2 m/s from 1983 to 2006 to above 6.0 m/s in 2013 which shows an increase in the wind power potential at this site. Hence from small to medium scale power generation, this may be suitable.Generally, the average wind speeds over this site are relatively low as compared to the cut-in wind speeds of most modern wind turbines which are usually between 6 and 10 m/s. Hence to generate wind energy on this site, a wind farm made of small wind turbines with low cut-in wind speed such as the AIRCON (HAWT 10 KW) or the Savonius wind turbine is recommended.Wind resource studies typically focus on high resolution data 10 min to 1-hour intervals. But this data and ground station data for the site are not available at the moment for analysis. Hence, in future, we will make sure we get in touch with high resolution data or ground station data from Cameroon meteorological department over this site. We will also be looking forward to designing more adaptive small wind turbines which can make use of the low wind speed over this site by concentrating on the aerodynamic structure, reduced load, and air-foil design with high lift and low Reynold number.



---

*Source: 9357812-2016-10-25.xml*