Open Access Library Journal
Vol.05 No.11(2018), Article ID:88922,14 pages

Modeling Wind Energy Using Copula

Zuhair Bahraoui1*, Fatima Bahraoui2, M. Amin Bahraoui3

1Department of Mathematics, ESTSB, University Chouaib Doukkali, El Jadida, Morocco

2Laboratory of Heat Transfer and Energetic, Faculty of Sciences and Techniques, Tangier, Morocco

3Department of Mathematics, Faculty of Science and Technology, Tangier, Morocco

Copyright © 2018 by authors and Open Access Library Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: October 13, 2018; Accepted: November 27, 2018; Published: November 30, 2018


In most studies related to wind energy, the quantity of the air density is considered constant, but actually, we know that it is variable and depending on others natural factors. We present a new procedure to estimate the wind density power energy by simulating the components of the air density. The procedure uses the copula theory and demonstrates that the estimated power energy is higher if the air density is not constant.

Subject Areas:

Applied Statistical Mathematics


Weibull, Nonparametric Estimation, Copula, Wind Speed

1. Introduction

Understanding the relationship between the components of wind power is fundamental to exploits wind energy, as well as the identification of suitable natural sites to perform the energy-efficient design and giving the economic rentability [1]. This involves that estimating wind energy in eolian installations can be considered as a challenge to predict the productive capacity energy. The most accurate formula to calculate the power in wind farms is to make use the wind density power energy

P = 1 2 ρ ( T , P ) S V 3 (1)

where ρ ( T , P ) is the air density. This term is dependent generally on the temperature T, and the pressure P, both of which vary with height, see for example [2]. S represent the useful surface of the wind sensor and v is the wind speed. In most studies, the air density is taken constant (for physical simplifications or for not having sufficient data), and it’s replacing with the standard air density, taking the average sea temperature 15˚C, and 1 atmospheric pressure that is 1.225 kg/m3. In this case, the relationship between the power and wind speed is reduced to a linear function.

In this paper, we use the copula theory to handle the structure dependence beyond the linear correlation between the naturals variables: The wind speed, temperature, and pressure. The air density is considered not constant and the power should be affected by the dependence which is the idea of this research.

The term of copula was introduced by [3] in the multivariate analysis concept, but it was not until 1986 when [4] and [5] focused more light on this type of function, giving an application in the finance field and exploiting the proprieties of the Archimedean copula. Applications have been extended in all fields [6] and of course for naturals sample type variable. One of the most salient properties of copulas theory is that it is possible to separate independently the function of the dependence from the marginals distribution function ones [7]. We use the same principle to treat the components of the air density, i.e. estimating an adequate copula for the wind speed and temperature and from another side an adequate copula for the wind speed and pressure. The marginals distribution function can be estimated independently.

[8] with a pioneering works find a suitable parametric copula for wind speed variables and its direction. Recently there is an important study relating to the dependence function in renewable wind energy. For example, [9] examines another type of copulas families to handle the dependency between wind speed and their directions [10] utilize a pair copula (conditional dependence) to analyze the correlation between the winds farms, while [11] use the extremes value theory and copulas function to determines the correlation between wind turbines that compose the wind farm. Others [12] examine the dependence between wind power production and electricity prices.

The remainder of this article is organized as: Section 2 describes the parametric and non-parametric distribution function to estimate the marginals probability distribution. The selected copulas used in the study are presented in Section 3. In Section 4 we explain a procedure to simulate wind speed density energy power. An application to real data is analyzed in the Section 5. Finally, we conclude.

2. Marginals under Study

Given X 1 , X 2 , , X n a n independent and identically distributed sample with a unknown distri-bution function (cdf) F and probability density function (pdf) f. A natural way to estimates the distribution function F, is to consider the empirical distribution function

F e m p = 1 n i = 1 n I ( X i x ) (2)

where I ( X i x ) is the indicator function of a set I ( X i x ) . The correction of this approximation

F ˜ ( x ) = n n + 1 F e m p ( x ) (3)

is considered in the empirical estimation because using directly (2) to estimate the copula function can cause the boundary problems. That is the copula distribution function implanted to identify the dependence structure may not integrate one, for a discussion one can consult [13]. This transformation is closed to the uniform distribution and its similar to the so-called the pseudo-observation introduced by [4]. The second nonparametric approximation is using the classical kernel estimator (CKE) [14]

F ^ ( x ) = 1 n i = 1 n K ( x X i h ) (4)

where K ( x ) = x k ( t ) d t is some known kernel distribution function [14] , like Epanechnikov or the Gaussian distribution function, the most used function used for this approximation, see for an efficiency comparison [15]. The parameter h, that satisfied generally the condition h → 0 when n → ∞ , controlled the smoothness of the estimated function.

In this paper, we use the specific value, h = 3.572 σ n 1 / 3 , where σ = min ( s d , R / 1.34 ) and sd is the standard deviation of the sample and R is the sample interquartile range [14].

Finally, to fit wind speed distribution function we add a Weibull distribution function. Used generally in reliability field and failure times in physical systems, this distribution function is recommended by the International Standard IEC 61400-12 and is widely used in the context of eolian analysis of a wind speed, [1] [2].

The CDF of Weibull distribution is: F k , λ = 1 e ( x λ ) k .

The pdf density function

f k , λ = { ( k λ ) ( x λ ) k 1 e ( k λ ) k s i x 0 0 s i x < 0

The mean and variance are, respectively:

E ( X ) = λ Γ ( 1 + 1 k ) and V a r ( X ) = λ 2 [ Γ ( 1 + 2 k ) + Γ 2 ( 1 + 1 k ) ]

such that, Γ ( x ) = 0 + t x 1 e t d t is the Euler’s Gamma function. To estimate the parameters
k , λ , exists various methods like the mean squared method or the moment method [16].

In this paper we use the maximum likelihood method to estimate the parameters [17]. This method is easy simple tools and consists to maximize the product density:

L V ( θ ) = i = 1 n f ( X i , θ ) (5)

where consist the set parameters to be estimated

3. Dependence e Model

To reveals the structure dependence between the pairs variable, (Wind speed, Temperature) and (Wind speed, Pressure), three parametric copula have been considered. A copula is a distribution function defined in the cubic interval [ 0 , 1 ] 2 with uniform marginal distribution functions U [ 0 , 1 ] . If FX, and FY are marginals distribution of variables bivariate ( X , Y ) , then from the Sklar theorem [3] for each bivariate distribution function H there exist a hidden cdf function called copula C such that H ( x , y ) = C ( F X ( x ) , F Y ( y ) ) . For more proprieties on copula one can consult [18] and [19].

The first copula considered is the Gaussian copula which pertains to the implicit copulas. Copula associated with elliptical distribution and represent a symmetrical dependency. In addition, they become important whether we are analyzing the right or left tail of the distribution function.

3.1. Gaussian Copula

If we denote by ρ be the linear correlation coefficient between two random variables X and Y, the Gaussian copula with parameter ρ is expressed:

C ρ ( u , v ) = Φ ρ ( Φ 1 ( u ) , Φ 1 ( v ) ) = Φ 1 ( u ) Φ 1 ( v ) 1 2 π 1 ρ 2 exp ( 2 ρ s t s 2 t 2 2 ( 1 ρ 2 ) ) d s d t

where Φ ρ is the two-dimensional standard Normal distribution function with correlation coefficient equal to ρ , and Φ is the standard Normal one-di- mensional distribution.

3.2. Sarmanov Copula

The second copula family considered is the Sarmanov copula. The range of the subfamilies is infinite due to its way of constructing a copula. One can find a special sub-copula associated to each marginal distribution function. Let ( X , Y ) be a bivariate random vector with marginal probability distribution functions (pdfs) f X and f Y . Also, let ψ 1 and ψ 2 two bounded non-constant function such that:

+ f X ( t ) ψ 1 ( t ) d t = 0 , + f Y ( t ) ψ 2 ( t ) d t = 0

The joint bivariate pdf introduced by [20] is defined as:

h ( x , y ) = f X ( x ) f Y ( y ) ( 1 + η ψ 1 ( x ) ψ 2 (y))

and the associated copula distribution function is:

C ( u , v ) = u v + η 0 u 0 v ψ 1 ( F X 1 ( t ) ) ψ 2 ( F Y 1 ( s ) ) d t d s (6)

The density is:

c ( u , v ) = 1 + η ψ 1 ( F X 1 ( u ) ) ψ 2 ( F Y 1 ( v ) ) (7)

where F X and F Y are the cumulative distribution functions (cdf’s) of X and Y, respectively. Parameter η is a real number that satisfies the condition for all x and y. 1 + η ψ 1 ( x ) ψ 2 ( y ) 0 for all x and y.

Note that when = 0, X and Y are independent. This parameter is related to the correlation between X and Y (if it exists), [21] as:

C o r r ( X , Y ) = η v 1 v 2 σ 1 σ 2 (8)

where ψ 1 ( x ) = x μ X and ψ 2 ( x ) = x μ X and μ X = E ( X ) and

μ Y = E ( Y ) (9)

To give a range of the parameter η, we use the result giving by [13] when the support of f X and f Y is not only belong in [ 0 , 1 ] i.e. if the support of f X is contained in [ a , b ] and that of f Y is contained in [ c , d ] where a, b, c and d are finite real numbers, then:

max ( ( b a ) ( d c ) ( μ X a ) ( μ Y b ) , ( b a ) ( d c ) ( b μ X ) ( d μ Y ) ) ( b a ) ( d c ) η min ( ( b a ) ( d c ) ( μ X a ) ( d μ Y ) , ( b a ) ( d c ) ( b μ X ) ( μ Y c ) ) (10)

3.3. Frank Copula

The ultimate copula is the Frank copula. This copula belongs to the so-called Archimedean copula, a family of dependence function used for here nice analytical proprieties [5]. The main characteristic of Frank copula is that it does not present dependence in the extremes, but in the center. The parameter of dependence can take a large range, θ ( , 0 ] [ 0 , + ) , and then one can handle the negative dependence.

The Frank copula is defined as:

C θ ( u , v ) = 1 θ ln ( 1 ( 1 e θ u ) ( 1 e θ v ) 1 e θ ) .

4. Simulate the Wind Density Power Energy from a Copula

In this section, we describe a procedure to simulate the statistical behavior of the wind power density using a Monte Carlo method [22]. From (1) and dividing by the surface area we have:

P w = 1 / 2 ρ ( T , P ) V 3 (11)

This term can represent the kinetic energy per unit area related to the wind. Now employing the ideal gas law, one expressed the air density [1] a

ρ ( T , p ) = 1.225 [ 288.15 T ] [ P 1013.3 ] (12)

The wind power density (11) can be calculated in two way. The first way is considering the air density ρ ( T , p ) constant. In this case the mean power produced until an observation z is:

A P w ( z ) = ρ ( T , p ) 2 0 z f ( v ) v 3 d v (13)

where f is the (pdf) of the wind speed. When the value z we have the average of the power energy. The second case if the air density ρ ( T , P ) , is not constant. For n registration of the data ρ ( T , p ) = ( ρ ( T 1 , P 1 ) , ρ ( T 2 , P 2 ) , , ρ ( T n , P n ) ) , the mean wind power energy can be calculated until the observation zk.

M P w ( z k ) = 1 2 i = 1 n ρ ( T k , P k ) z k 3 (14)

Now to simulate the wind power energy density we start by simulating the wind speed variable coupled with the temperature and the pressure.

The same procedure describes in [19] to generate a pair of copula is used. To generate a two-dimensional random variable we do serve the conditional distribution of the random vector ( U , V ) :

P ( V U | U = u ) = C u (v)

where C u ( v ) = lim Δ u 0 + C ( u + Δ u , v ) C ( u , v ) Δ u = C ( u , v ) u .

The following algorithm simulate the wind density power energy:

1) Start by fixing a copula C, of wind speed and temperature.

2) Generate two independent random variables u1 and z from a Uniform distribution U ( 0 , 1 )

3) Set u 2 = C u 1 [ 1 ] ( z ) , where C u 1 [ 1 ] denotes aquasi-inverse of C u 1 . The quasi-inverse is:

C u 1 [ 1 ] ( z ) = { inf { x | C u 1 z } if z = 0 C u 1 1 if z ( 0 , 1 ) inf { x | C u 1 t } if z = 1

4) The desired first pair is ( u 1 , u 2 ) where u2 is a uniform variable related to the temperature variable.

5) Fixing the random variable u1, and considering now the copula of wind speed and pressure, we repeated the procedure to give a pair ( u 1 , u 3 ) , where u3 is a uniform variable related to pressure variable.

6) Taking the inverse F 1 ( u i ) , i = 1 , 2 , 3 , of the marginal distribution function used in (3), (for the inverse method one can consult [22] ), the pairs variable (Wind speed, temperature, pressure) coupled and conserving the same structure of the dependence in sample data.

7) Replacing this terms in (11).

5. Data and Result

The data represent the registration of wind speed, temperature, and pressure collected in the region Hrarza, situated in the north Morocco kingdom. Near on the straits of Gibraltar and surrounded by two seas, the Mediterranean and the Atlantic, this region suffers a gusty wind. The annual mean wind speed exceeds 6 m/s. For comparison with another region of the kingdom, we can consult [23] in the north of Morrocan. The Hrarza station covers the registration at the beginning of 2000. We use only the last years of the registration covering 365-day maximal wind speed and their correspondent temperature and pressure.

The first lecture of the graphical behavior of this variable Figure 1 is that the temperature variables are more predictable than the wind speed and the pressure. Probably due to the seasonal comportment. In Table 1 we give the descriptive statistics of the data. A simple reading for wind speed variable indicates that the distribution can be right skewness. The Jarque-Bera test of normality indicates that either variable is distributed Normally.

As we noted in Section 2 the marginal adjustment is resumed in Figure 2 and Figure 3. For the temperature and the pressure variables we have used the nonparametric estimation and for the wind speed, we have added the Weibull adjustment. We note that the nonparametric estimation of the distribution gives a well fit for the three variables. To visualize the relationship between the pair (Wind speed, Temperature), (Wind speed, Pressure) we make use the regression

Figure 1. Daily wind speed, temperature and pressure.

Figure 2. Distribution and density of wind speed. (a) Density of the wind speed; (b) CDF of the wind speed.

Figure 3. Density function of temperature and pressure. (a) Density of temperature; (b) Density of pressure.

Table 1. Descriptive statistics.

and a scatter plot, Figure 4 and Figure 5 which indicate a negative correlation if one considers a linear regression. On the other hand, the Kplot, [24] a useful

Figure 4. Regression of wind speed and temperature. (a) Wind vs Temperature; (b) Regression line.

Figure 5. Regression of wind speed and pressure. (a) Wind vs pressure; (b) Regression line.

tool to identify the dependency structure between two variables confirm this relation, i.e. the point plotted is under the diagonal line wish indicate a negative dependency.

To visualize the intensity of the dependence in the tail and in the center the nonparametric estimation of copula density, Figure 6 & Figure 7 is also considered. For an independent and identical variables ( U i , V i ) , i = 1 , 2 , , n the nonparametric estimation of copula density can be expressed:

c ^ ( u , v ) = 1 n i = 1 n k ( U i u h * ) k ( V i v h * ) , (15)

where k is the kernel pdf function and h * is one of the smoothed parameters for copula estimation [25].

Figure 6. Nonparametric density copula. (a) Wind and temperature; (b) Wind and pressure.

Figure 7. K-Plot for dependent data. (a) K-plot wind and temperature; (b) K-plot wind and pressure.

The package kdecopula [26] is used for R available language and environment for statistical computing and graphics. We observe that the associated copula for both, (Wind speed, Temperature) and (Wind speed, Pressure) doesn’t have extremes dependance in the tail, but a large dependence in the center and tend to have symmetry form. The elliptical or an Archimedean copula can be a good fit for this type of dependence. In the Table 2 we summarize the estimate parameters for a proposed copula accompanied with the CIC criteria value to select the adequate statistics model. The method used in estimating the parameters of the copula is the inference from marginal, (IFM) introduced by [18]. This method consists firstly to estimate the distribution function of the marginal (CDF), and calculating the pair ( F ^ X ( X i ) , F ^ X ( Y i ) ) = ( u ^ i , v ^ i ) , i = 1 , 2 , , n using (5). The second stage is to minimize the log-likelihood function for a copula, that is L V C ( θ ) = i = 1 n C θ ( u ^ i , v ^ i ) .

The (IFM) method has an advantage that they avoid the excess time optimization. Using a global estimation of the log-likelihood function considering both the marginals and a copula doesn’t guarantee the existence of the minimum, [27].

Table 2. Results of the fit of the copulas with nonparametric and semiparametric marginals.

The CIC criteria give the best selection of copula among another copula. Like the Akaike information criterion [28] , the CIC tools was introduced by [29] and it is designed when the dependence relation is considered. Analyzing the result of the Table 2 we can say that the Sarmanov copula gives the best fit for the dependency between wind speed and temperature and more specify using a Weibull distribution for wind speed and the empirical (CDF) for a temperature as marginals distribution function. In another hand, Frank copula give the best fit for the wind speed and Pressure and one can only use the empirical estimation for the marginals to give the best model.

A simulation of daily energy is illustrated in Figure 8(a) using the procedure describes above in Section 4. We can generate then the best model that reflect the dependence structure of the variables coupling with their estimated marginal probability function. With the red lines, we have the wind density power (11) simulated when the air density is constant and with black lines if the air density is no constant. We note that there exists a notable difference between the two cases, and generally, the power energy is higher if the air density not constant. Finally, we plot, Figure 8(b), the mean of the wind density power energy simulated both in the case (13) and (14). The variation of mean in the case of the density is not constant is hight, stabilized at the total mean 252.181. Using the Weibull pdf function the mean of the density power energy is equal 214.8 less.

6. Conclusion

In this work, we have analyzed the statistical behavior of the wind energy in

Figure 8. Power energy. (a) Power from sarmanov (Weibull, Emp) and frank (Weibull, Emp) copula; (b) Simulate aggregate power.

north Morocco, from the point the view of the copula. Our objective is to give a most accurate simulation method for a density power energy in the eolian park. Fitting the marginal probability component of the wind density variable, we have noted that the nonparametric approach can give the best fit. To capture the negative dependence between the variables we have used two type of copula. The Archimedean and the elliptical copula. We have incorporated Sarmanov copula wish has never been used in this type of data and giving a suitable model. The procedure introduced to simulate the wind density power energy can generalize to the other field related to renewable energies since the density of the air is always present.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper and there has been no significant financial support for this work that could have influenced its outcome.

Cite this paper

Bahraoui, Z., Bahraoui, F. and Bahraoui, M.A. (2018) Modeling Wind Energy Using Copula. Open Access Library Journal, 5: e4984.


  1. 1. Burton, T., Jenkins, N., Sharpe, D. and Bossanyi, E. (2001) Wind Energy Handbook. John Wiley and Sons, Inc., Chichester.

  2. 2. Manwell, J., McGowan, J. and Rogers, A.L. (2010) Wind Energy Explained: Theory, Design. John Wiley and Sons.

  3. 3. Sklar, A. (1959) Fonctions de répartition a′n dimensions et leurs marges. Publications de Institut de Statistique de l’université de Paris, 8, 229-231.

  4. 4. Genest, C. and MacKay, J. (1986) The Joy of Copulas: Bivariate Distributions with Uniform Marginals. The American Statistician, 40, 280-283.

  5. 5. Genest, C. (1987) Frank’s Family of Bivariate Distributions. Biometrika, 74, 549-555.

  6. 6. McNeil, A.J., Frey, R. and Embrechts, P. (2005) Quantitative Risk Management: Concepts, Techniques, and Tools. Princeton University Press, Princeton.

  7. 7. Genest, C. and Favre, A.C. (2007) Everything You Always Wanted to Know about Copula Modeling but Were Afraid to Ask. Journal of Hydrologic Engineering, 12, 347-368.

  8. 8. Coles, S. and Walshaw, D. (1994) Directional Modeling of Extreme Wind Speeds. Applied Statistics, 43, 139-157.

  9. 9. Soukissiana, T.H. and Karathanasia, F.E. (2017) On the Selection of Bivariate Parametric Models for Wind Data. Applied Energy, 188, 229-236.

  10. 10. Cao, J. and Yan, Z. (2017) Probabilistic Optimal Power Flow Considering Dependences of, Wind Speed among Wind Farms by the Pair-Copula Method. International of Electrical Power and Energy Systems, 84, 296-307.

  11. 11. D’Amicoa, G., Petronib, F. and Pratticoc, F. (2015) Wind Speed Prediction for Wind Farm Applications by Extreme Value Theory and Copulas. Journal of Wind Engineering and Industrial Aerodynamics, 145, 229-236.

  12. 12. Pircalabu, A., Hvolbya, T., Jung, J. and Hog, H. (2017) Joint Price and Volumetric Risk in Wind Power Trading: A Copula Approach. Energy Economics, 62, 139-154.

  13. 13. Bolancé, C., Bahraoui, Z. and Artis, M. (2014) Quantifying the Risk Using Copulate with Nonparametric Marginals. Insurance: Mathematics and Economics, 58, 46-56.

  14. 14. Silverman, B.W. (1986) Density Estimation for Statistics and Data Analysis. Chapman and Hall, CRC Finance Series, London.

  15. 15. Wand, M.P. and Jones, M.C. (1995) Kernel Smoothing. Chapman and Hall, London.

  16. 16. Pobockova, I. and Sedliackova, Z. (2014) Comparison of Four Methods for Estimating the Weibull Distribution Parameters. Applied Mathematical Sciences, 8, 4137-4149.

  17. 17. Cohen, A.C. (1965) Maximum Likelihood Estimation in the Weibull Distribution Based on Complete and on Censored Samples. Technometrics, 7, 579-588.

  18. 18. Joe, H. (1997) Multivariate Models and Dependence Concept. Chapman and Hall, London.

  19. 19. Nelsen, R.B. (2006) An Introduction to Copulas. 2nd Edition, Springer, Portland.

  20. 20. Sarmanov, O.V. (1966) Generalized Normal Correlation and Two-Dimensional Fréchet. Soviet Mathematics. Doklady, 25, 1207-1222.

  21. 21. Lee, M.L. (1996) Properties and Applications of the Sarmanov Family of Bivariate Distributions. Communications in Statistics-Theory and Methods, 25, 1207-1222.

  22. 22. Devroye, L. (1986) Non-Uniform Random Variate Generation. Springer, New York.

  23. 23. Nfaoui, H., Buret, J.J. and Sayigh, A.A.M. (1998) Wind Characteristics and Wind Energy Potential in Morocco. Solar Energy, 63, 51-60.

  24. 24. Genest, C. and Boies, J.C. (2003) Detecting Dependence with Kendall Plots. Journal of the American Statistical Association, 57, 275-284.

  25. 25. Nagler, T. (2014) Kernel Methods for Vine Copula Estimation. Universi at Munchen, München.

  26. 26. Nagler, T. (2017) Kdecopula: An R Package for the Kernel Estimation of Bivariate Copula Densities.

  27. 27. Cherubini, U., Luciano, E. and Vecchiato, W. (2004) Copula Methods in Finance. Wiley, Chichester.

  28. 28. Akaike, H. (1973) Information Theory and an Extension of the Maximum Likelihood Principle. 2nd International Symposium on Information Theory, Tsahkadsor, 2-8 September 1971, 267-281.

  29. 29. Gr?nnrberg, S. and Hjort, N.L. (2014) The Copula Information Criteria. Scandinavian Journal of Statistics, 41, 436-459.