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 . 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
where 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 . 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  in the multivariate analysis concept, but it was not until 1986 when  and  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  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 . 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.
 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,  examines another type of copulas families to handle the dependency between wind speed and their directions  utilize a pair copula (conditional dependence) to analyze the correlation between the winds farms, while  use the extremes value theory and copulas function to determines the correlation between wind turbines that compose the wind farm. Others  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 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
where is the indicator function of a set . The correction of this approximation
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 . This transformation is closed to the uniform distribution and its similar to the so-called the pseudo-observation introduced by . The second nonparametric approximation is using the classical kernel estimator (CKE) 
where is some known kernel distribution function  , like Epanechnikov or the Gaussian distribution function, the most used function used for this approximation, see for an efficiency comparison . 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, , where and sd is the standard deviation of the sample and R is the sample interquartile range .
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,  .
The CDF of Weibull distribution is: .
The pdf density function
The mean and variance are, respectively:
is the Euler’s Gamma function. To estimate the parameters
, exists various methods like the mean squared method or the moment method .
In this paper we use the maximum likelihood method to estimate the parameters . This method is easy simple tools and consists to maximize the product density:
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 with uniform marginal distribution functions . If FX, and FY are marginals distribution of variables bivariate , then from the Sklar theorem  for each bivariate distribution function H there exist a hidden cdf function called copula C such that . For more proprieties on copula one can consult  and .
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:
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 be a bivariate random vector with marginal probability distribution functions (pdfs) and . Also, let and two bounded non-constant function such that:
The joint bivariate pdf introduced by  is defined as:
and the associated copula distribution function is:
The density is:
where and 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. 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),  as:
where and and and
To give a range of the parameter η, we use the result giving by  when the support of and is not only belong in i.e. if the support of is contained in and that of is contained in where a, b, c and d are finite real numbers, then:
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 . 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, , and then one can handle the negative dependence.
The Frank copula is defined as:
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 . From (1) and dividing by the surface area we have:
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  a
The wind power density (11) can be calculated in two way. The first way is considering the air density constant. In this case the mean power produced until an observation z is:
where f is the (pdf) of the wind speed. When the value we have the average of the power energy. The second case if the air density , is not constant. For n registration of the data , the mean wind power energy can be calculated until the observation zk.
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  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 :
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
3) Set , where denotes aquasi-inverse of . The quasi-inverse is:
4) The desired first pair is 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 , where u3 is a uniform variable related to pressure variable.
6) Taking the inverse , of the marginal distribution function used in (3), (for the inverse method one can consult  ), 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  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.
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 the nonparametric estimation of copula density can be expressed:
where k is the kernel pdf function and is one of the smoothed parameters for copula estimation .
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  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 . This method consists firstly to estimate the distribution function of the marginal (CDF), and calculating the pair using (5). The second stage is to minimize the log-likelihood function for a copula, that is .
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, .
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  , the CIC tools was introduced by  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.
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.
 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.
 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.
 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.
 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.
 Bolancé, C., Bahraoui, Z. and Artis, M. (2014) Quantifying the Risk Using Copulate with Nonparametric Marginals. Insurance: Mathematics and Economics, 58, 46-56. https://doi.org/10.1016/j.insmatheco.2014.06.008
 Pobockova, I. and Sedliackova, Z. (2014) Comparison of Four Methods for Estimating the Weibull Distribution Parameters. Applied Mathematical Sciences, 8, 4137-4149.
 Lee, M.L. (1996) Properties and Applications of the Sarmanov Family of Bivariate Distributions. Communications in Statistics-Theory and Methods, 25, 1207-1222.