A reliability model of wind farm located in mountainous land with complex terrain, which considers the cable and wind turbine (WT) failures, is proposed in this paper. Simple wake effect has been developed to be applied to the wind farm in mountainous land. The component failures in the wind farm like the cable and WT failures which contribute to the wind farm power output (WFPO) and reliability is investigated. Combing the wind speed distribution and the characteristic of wind turbine power output (WTPO), Monte Carlo simulation (MCS) is used to obtain the WFPO. Based on clustering algorithm the multistate model of a wind farm is proposed. The accuracy of the model is analyzed and then applied to IEEERTS 79 for adequacy assessment.
1. Introduction
Recent public concerns on environmental issues associated with conventional fossil resources has motivated an increased interest in the development and use of alternative sources. Wind energy has been widely recognized as an important alternative with the development in modern WT technology. Many countries, e.g. US, China, Germany and Denmark have integrated large scale wind energy into their power systems with incentive politics promoted which ensure the purchase of the wind energy produced in spite of its uncompetitive prices
[1

3]
.
However the stochastic WFPO, depending on the wind, which is intermittent, has become one major difficulty while planning or operation a power system with increasing wind energy penetration. The impact of wind farms (WFs) on composite power system reliability is a major concern, especially after a series of large scale accidents of wind farms disconnected with the power grid occurred in China since 2011
[4]
. The intermittent nature of wind speed, along with the probabilistic behavior of outage of WTs and other components, makes output power of wind farms completely stochastic and different from those of conventional units. Therefore, one of the complexities of integration of wind power in power systems can be seen in the reliability model of wind farms.
Traditionally, the analytic and MCS approaches are always used to model the wind farm
[5

11]
. Generally a wind farm can be simplified to an equivalent multistate unit
[7
,
10

11]
, which can be used for reliability evaluation of a power system with wind farms.
Analytic approach is always used when the wind farm is not large. Sayas combines the stochastic characteristics of wind speed with the operational information of the turbines, such as the failure and repair rates, representing the wind farm by a Markov process
[5]
. Leite utilizes the model for reliability studies, providing an annual estimate of energy production and calculate several performance indexes
[6]
. Then, based on Markov process Dobakhshari models a wind farm like a multistate conventional unit
[7]
.
Although the model is more credible and reliable, with the wind energy integrated into the power grid in large scale the problem becomes more complex and difficult, especially after incorporating wake effect and other relevant factors in a Markov process.
Roy Billinton proposed MCS to mimic the operation of a wind farm hourly, taking into account the fluctuating characteristics of wind speeds, the random failure of generating units and other recognized factors
[8]
. Based on the model, assessment of riskbased capacity is carried out for the power system including wind farms
[9]
. In 2006 Karki and Billinton, etc. presents a simplified multistate unit model for reliability evaluation of power systems
[10]
. And the method is further simplified by determining the reasonable multistate representation for a wind farm generation model in reliability evaluation.
Compared with the analytical one, MCS is more applicable to model the wind farm with numerous WTs and always used in an existing composite power grid. Thus, in this paper MCS is selected. The wind farm model is mainly determined by the random nature of the site resource, the characteristic of WTPO, and the unavailability of the WT expressed by the unit forced outage rate (FOR)
[8

11]
. However, two main aspects should be further investigated:

1. More factors should be taken into account due to the differences between wind farms.

2. How to achieve the tradeoff between accuracy and the number of states.
The existing models consider only the wake effect of a wind farm in a flat area. Based on the relatively simple Jensen model, a partial wake effect model in flat terrain between WTs is established
[12

14]
. When the wind passes over a mountainous land, complex terrain feature will change the wind speed. Wind speed over a mountain is much different from that on a flat area. Therefore wake effect model in flat terrain is not sufficient to accurately evaluate power outputs of WTs at different parts of a wind farm. The terrain influence of a mountainous land on wind speed has not been considered in an adequacy study of a wind farm. Besides, the internal connection line failures in the wind farm have been recognized as a contributing factor to wind power system reliability, after three accidents of a cable breakdown in wind farms and the impact of it on wind farm reliability is investigated. In this paper the number of states and the accuracy of the model is investigated, providing a reference to the tradeoff. The overall idea of the paper is illustrated in
Fig. 1
.
The overall idea of this paper
2. Wind Turbine Power Output Considering Multiple Wake Effect in Complex Terrain
 2.1 Wind speed model
The probability density function (PDF) of wind speed is important in numerous wind energy applications. A large number of studies have been published in scientific literature related to renewable energies that propose the use of a variety of PDFs to describe the wind speed frequency distributions, such as Weibull, Rayleigh,
r
^{2}
and so on
[15]
. Amongst the various distributions, the two parameter Weibull distribution has been the most widely used, accepted and recommended distribution on wind energy
[16]
. In this paper, two parameter Weibull distribution is used to simulate the wind speed probability distribution, its probability density function (PDF) is given by
[17]
:
where
c
is the scale factor and
k
is the shape factor. Seasonal patterns of wind speed which significantly affect the reliability indexes is taken into account
[6]
. The values of
c
and
k
can be estimated using maximum likelihood method based on measured wind speeds
[16]
,
[18]
.
Cumulative distribution function (CDF) of two parameter Weibull distribution is given by:
 2.2 Partial wake effect model in complex terrain
WTs generate electricity by tapping into the energy in the wind. Consequently, the wind leaving the turbine must have a lower energy content than that arriving in front of the turbine
[14]
. In this paper, a detail wake model, considering partial and multiple wakes inside a wind farm, has been developed.
The Jensen wake model is a single wake model, which demonstrates linear expansion of the wake diameter, ignoring the changes of the wind speed caused by the altitude difference between the WTs.
Fig. 2
shows an overview of the wake model.
The Jensen wake model
When completely located in the wake of the upstream WT, the wind speed at location
x
is expressed by
[14]
:
where
v
_{0}
is the original wind speed,
d_{F}
is the declining coefficient of wind speed in flat terrain,
C
_{T}
is the thrust coefficient,
r=r_{0}+k_{w}x, r
is the wake effect radius,
k_{w}
depends on the terrain and environment,
k_{w}
=0.04 for nature wind, otherwise
k_{w}
=0.08
[14]
.
However; the wake effect model in a more general terrain should be proposed for the wind farm located in mountainous land where the difference of altitude between WTs is not negligible. Besides the difference of altitude and the wind direction may probably lead to the downstream WT in partial wake. Combining the basic theories of Jensen model and Lissaman model
[19]
, a partial wake effect model of WTs at the different altitudes in complex terrain is developed in this paper.
In a wind farm configuration, the wind speed at WT
j
is affected not only by the upstream WT that is directly in front of it, WT
i
, but also by other upstream WTs
[14]
. At first, only the effect on WT
j
caused by WT
i
is studied in this section. In a wind farm, the relative position of WT
i
and WT
j
is assumed to be shown in
Fig. 3
. In the figure,
h_{i}
and
h_{j}
are the rotor altitudes of WT
i
and WT
j
,
r_{i}
and
r_{j}
are the rotor radiuses of WT
i
and WT
j
and
x_{ij}
is the distance between WT
i
and WT
j
. Then the wind speed of WT
j
affected by WT
i
can be calculated using the following expression:
Partial block wake effect model
where
v_{i}
and v
_{j(i)}
are the wind speeds at WT
i
and WT
j
;
α
is the coefficient of wind speed variation with the altitude, which is generally 1/7.
r_{i}
(
x_{ij}
) is the wake effect radius at WT
j
caused by WT
i
;
r_{i}
(
x_{ij}
)=
r_{i}+k_{w}x_{ij}
;
A_{rotor}
is the area covered by the rotor of downstream WT
i
,
is the area projected by downstream WT
j
in the wake region of upstream WT
i
, which can be calculated through the following expression:
where Δ
h
=
h_{j}−h_{i}
;
d_{ij}
is the horizontal spacing between WT
i
and WT
j
. When Δ
h
=
h_{i}−h_{j}
=0 and
A_{shad.ij}=A_{rotor}
, it turns to be the complete wake effect model in flat terrain expressed by Eq. (3).
 2.3 Multiple wake effect model in complex terrain
It is quite possible that one WT is affected by a multiple wake effect in a wind farm containing a large number of WTs. Therefore the multiple wake effect model in complex terrain should be taken into account. This model assumes that the kinetic energy deficit of interacting wakes is equal to the sum of the energy deficits of the individual wakes
[20]
. Thus, the velocity at the intersection of several wakes can be calculated using the following expression:
where
v_{j0}
is the wind speed at WT
j
without any wake effect;
N
is the total number of upwind influencing turbines;
v_{eq.j}
is the wind speed at WT
j
such that all the wakes are taken into account;
v_{j}
is the wind speed at WT
j
affected by the individual WT
i
.
 2.4 Wind turbine power output model
The WTPO is a function of the wind speed. The randomness of wind speed determines the WTPO. The relationship between the WTPO and the wind speed is nonlinear due to the combined effects of aeroturbine and generator characteristics
[21]
. However, the assumption of the linear characteristic of the WTPO with the wind speed simplifies the calculations without roughly errors. Ignoring the differences between WT types, the commonly used linear relationship between the WTPO and the wind speed can be calculated using the following expression
[22]
:
where
P_{r}
,
v_{ci}
,
v_{r}
and
v_{co}
are the rated power output, the cutin wind speed, the rated wind speed, and the cutout wind speed of the WT, respectively;
k
_{1}
=
P_{r}
/(
v_{r}−v_{ci}
);
k
_{2}
=−
k
_{1}
v_{ci}
.
3. Wind Farm Reliability Model
 3.1 Wind farm element reliability model: Connection lines and wind turbines
Generally, four main parts: WTs, package transformers, collection lines, substations constitute the wind farm. In present techniques, the impact of connection line reliability to wind farms has not been considered in wind farms reliability modeling. The lines, which are long with high failure rates and are used for collecting and delivering energy, have a great impact on the wind farm. A series of large scale accidents of wind farms disconnected with the power grid occurred in China since 2011 just due to the breakdown of cable
[4]
.
 3.2 The reliability model of the connection line
The model of cable is investigated in this section as an example, but the model is also applicable to common connection lines like overhead line. Electric power network topology in the wind farm comprises 3 general forms: chain, ring and star
[23]
. Twostate model is applied to each cable, representing the operation (up) and outage (down) state. The Markov chain of cable is illustrated in
Fig. 4
, where
λ_{cable}
,
μ_{cable}
are failure rate and repair rate.
Markov chain model of the cable
The probability calculated from the Markov chain above is as follows:
Only longterm state probabilities are of interest, the normal and failure state probabilities of each component are expressed as follows:
 3.3 The Reliability of the Wind Turbine
Operation of the WT is affected by the connecting cable. Some outage of the WTs is due to the cable failure and the equivalent outage of WTs increase. The model of the WTs is illustrated below.
Markov chain model of the WT
The common cause outage indicates the WT outage due to the cable failures. The WT outage caused by the failure itself is alone to be separated from the common cause outage. The probability of WT being available
can be calculated below:
Where
is the equivalent availability of the
i
th WT considering the common cause outage,
P_{WTi}
is the availability of the
i
th WT,
F_{i}
is a failure set which makes the
i
th WT unavailable due to the common cause of cable,
f
represents a failure that belongs to
F_{i}
.
Fig. 6
below shows 10 WTs connected to the bus in chain form.
Chain form of the WTs
The equivalent availability of each turbine is demonstrated below in
Fig. 7
.
Equivalent availability of the WTs
As it is seen in the figure above, the equivalent availability decreases with the number of cables used to connect into the grid increases.
4. Wind Farm MultiState Probabilistic Model Based On KMeans Clustering
 4.1 Wind farm power output
The MCS is proposed to mimic the WFPO. The step of the calculation is as follows:

1) Input the location coordinates of each WT.

2) Generate the wind speed and wind direction.

3) Obtain the state of WTs and cables according to MCS.

4) Calculate the equivalent wind speed at each location of WT, especially the WT in the state of outage does not impact the WTs downstream.

5) Calculate all the WTPO according to the equivalent wind speed.

6) Calculate the WFPO.

7) Determine whether it is the end of the simulation time. If it is, end the simulation and generate the output. If not, return to step 2).
 4.2 MultiState probabilistic model based on clustering
In order to apply time series WFPO simulated hourly to power system reliability evaluation, the multistate probabilistic model is obtained through clustering the WFPO data.
Clustering is a process of partitioning or grouping a set of data objects into a number of clusters such that similar patterns are assigned to one cluster
[24]
.
K
means clustering is one of the simplest unsupervised learning algorithms that solve the wellknown clustering problem. In this paper,
K
means clustering will be used to partition a time series of WFPO into
K
different clusters and the wind farm will be equivalent to a
K
state conventional unit.
The
K
means clustering attempts to find the cluster centers (
c
_{1}
,
c
_{2}
, … ,
c
_{K}
) such that the sum of the squared distances(called Distortion) of each data point (
x_{i}
) to its nearest cluster center (
c
_{k}
) is minimized, as shown by following expression:
where
d
(
x_{i}
,
c_{k}
) is the Minkowski distance function and its expression is given by:
where
p
is a positive integer;
m
is the dimension of vector
c_{k}
.
Typically
d
(
x_{i}
,
c_{k}
) is chosen as the Euclidean distance, which is a special case (
p
=2) of the Minkowski distance.
The
K
means clustering includes the following steps
[24]
:

1) Fix the number of clustersKand initialize the center locations of each cluster, (c1(0),c2(0),
,cK(0)). Each cluster center is amdimensional vector i.e.,ci(0)={ci1(0),ci2(0), …,cim(0)}.

2) Start the iterative procedure. Set the iteration count t =1.

3) Calculate the distance measure,dki(t−1), betweenkth cluster center andith data set (data point inmspace). The distance measure used is Euclidean distance as given by Eq. (15).

4) Assign each data objectxito its nearest cluster centerck.

5) Update each cluster centerck(t)as the mean of allxithat have been assigned as closest to it as given by Eq. (16).
where
n_{k}
is the number of data items belonging to
k
th cluster.

6) Calculate the DistortionDisas given by Eq. (13). This function depicts the sum of all intra cluster distances, which is lower the better.

7) If the value ofDishas converged, return the final cluster centers (c1(t),c2(t), …,cK(t)). Else sett=t+1 and go to (3).
5. Example Analysis
 5.1 Relevant data of wind farm
In this paper, a wind farm containing 2MW WTs with a total of 160MW installed (the main technical parameters of the WT are given in
Table 1
[25]
.) is selected as an example to do the relevant research and its specific layout is shown in
Fig. 8
. The topography of the wind farm is complex and the rotor altitude of each WT varies. The wind speed and wind direction of the wind farm are measured by the anemometer tower, which is at the altitude of 70
m
. Besides, the data next to the WT in
Fig. 8
is the rotor altitude of the WT.
The main technical parameters of asynchronous wind turbine
The main technical parameters of asynchronous wind turbine
The layout of wind turbines in wind farm
Due to the randomness of wind direction, wind rose is used to describe the wind direction. And the wind rose of the selected wind farm is shown in
Fig. 9
. In this paper, twoparameter Weibull distribution is used to simulate the wind speed. Thus, assuming the values of
c
and
k
change monthly, the values of
c
and
k
from Capo Vado wind farm in
[26]
are used for further analysis in this paper and the values of
c
and
k
are presented in
Table 2
.
Wind rose of the wind farm
Monthly Weibull parameters for Capo Vado site
Monthly Weibull parameters for Capo Vado site
 5.2 Modeling of wind farm power output probability distribution
Using a 50year time series of wind speed and wind direction simulated hourly, a time series of power output of the wind farm (seen in
Fig. 8
) is calculated with a multiple wake effect model in complex terrain taken into consideration. With
K
means clustering used to divide the time series of power output of the wind farm into
K
(
K
= 6,7,8,9,10) different clusters, multistate WFPO probability distribution models is shown in
Table 3
.
Multistate wind farm power output probability distribution models
Multistate wind farm power output probability distribution models
Table 3
shows that the wind farm mentioned in Section 5.1 is equivalent to a (6~10)state conventional unit. And the maximum of WFPO increases with the increase of the state number; however, its corresponding probability decreases.
In order to study the accuracies of using (6~10)state WFPO probability distribution models to simulate the WFPO, the error curves of comparing the simulated power output with the actual power output of the wind farm are plotted in
Fig. 10
.
Error curves of multistate wind farm power output probability distribution models
From
Fig. 10
, it can be noted that the errors between the simulated power output and the actual power output of the wind farm are volatile, but the overall errors are all less than 1% when the simulation times is about 50% of the total simulation times (438,000 times, 50 years, simulating hourly). Compared with the mean and standard deviation of the time series of the actual WFPO calculated in
Fig. 10
when the simulation times is 438,000, the errors of the means and standard deviations of the simulated time series of the WFPO are shown in
Table 4
.
Errors of the mean and standard deviation of every multistate model
Errors of the mean and standard deviation of every multistate model
Combining with
Fig. 10
(when the simulation times is 438,000, the errors of (6~10)state models: 7>8>6>10>9) and
Table 4
, the means reflect the accuracy rate of the multistate models and the standard deviations are related to the volatility of the multistate models to simulate WFPO. Because of the probability in the process of simulation, the accuracy rates of the multistate models do not increase with the state number, but do the volatility of the multistate models.
Due to the uncertainties of the simulation accuracies of the (6~10)state models, boxplot is used to describe the volatility of the error of every multistate model in simulating WFPO (seen in
Fig. 11
). The boxplot is obtained through using every multistate model to simulate WFPO for 100 simulations (the time of every simulation is 50 years).
Boxplot of every multistate model
As it can be seen from
Fig. 11
, the average error of every multistate model in simulating WFPO varies from 0.1% to 0.2% and the maximum error is not more than 0.5%. Therefore, every multistate model can ensure its reliability and accuracy when simulating WFPO. In other words, the multistate models can effectively simulate WFPO.
 5.3 Wind farm adequacy assessment
In this paper, the IEEERTS 79
[27]
Test System with the wind farm (in Section 5.1) integrated into the 16th node is applied to carry out the adequacy assessment. And the multistate models (seen in
Table 3
) are used to simulate the WFPO to investigate the impact of multistate models on the IEEEBTS 79 generating system adequacy. The reliability indexes LOLP (Lost of Load Probability) and EDNS (Expected Demand Not Supplied) are calculated for the corresponding analysis.
Fig. 12
shows the results of LOLP and EDNS calculated using (6~10)state models to simulate WFPO. In
Fig. 12
, the LOLP and EDNS both fluctuate slightly with the state number, basically in accord with the reliability impact of the actual WFPO on the original test system.
LOLP and EDNS of every multistate model
Fig. 13
shows the values of EDNS change with the various annual peak loads ranging from 0.9p.u. to 1.1p.u. (the peak load of the test system is 2850 MW) with (6~10)state models integrated into the test system. It can be seen from
Fig. 13
that the values of EDNS are relatively close, using a model with six or more states to represent the wind farm. The differences between the EDNS using the actual WFPO and the EDNSs using (6~10)state models to represent the wind farm are relatively small but increase with the various annual peak loads ranging from 0.9p.u. to 1.1p.u..
Curves of EDNS with the various annual peak loads
From above analysis, the multistate models can be used to represent wind farm in reliability evaluation of a power system with wind farms.
In order to analyze the accuracies of using multistate models to represent wind farm in reliability evaluation of a power system with wind farms, the boxplots of LOLP error and EDNS error are obtained through 50 simulations (the time of every simulation is 50 years) (seen in
Fig. 14
). The errors of LOLP or EDNS are calculated between the one calculated using the actual WFPO and the ones calculated using (6~10)state models to represent the wind farm in the reliability evaluation of a power system with wind farms.
Boxplots of LOLP error and EDNS error
Fig. 14
shows that the average error of LOLP is volatile (varies between 1% and 2%) and the maximum error is not more than 5%; the average error of EDNS decreases with the increasing state number (varies between 2% and 3%) and the maximum error is not more than 7%. Therefore, the multistate models can effectively represent wind farm in reliability evaluation of a power system with wind farms and ensure the accuracy.
6. Conclusion
This paper has introduced a modeling approach for WFPO probability distribution based on clustering algorithm. Firstly, combining the basic theories of Jensen model and Lissaman model, a partial wake effect model of WTs at the different altitudes in complex terrain is established. A wind farm reliability model, especially considering the outage of wind farm internal connection lines is set up in this paper. Finally, after an introduction of
K
means clustering in detail, a wind farm can be equivalent to a multistate conventional unit with
K
means clustering as a method.
Taking an installed capacity of 160 MW wind farm as an example, it turns out that the multistate models can effectively simulate the WFPO. The studies on the IEEERTS 79 Test System, with reliability indexes LOLP and EDNS calculated, indicate that the multistate models can also effectively represent wind farm in reliability evaluation of a power system with wind farms and ensure the accuracy.
Acknowledgements
This work was supported by a grant from the National High Technology Research and Development Program of China (863 Program) (No.2013BAA02B00).
BIO
Wenxia Liu She received B.S. degree in Radio technology from Nanjing University of Science and Technology in 1990, and obtained M.S. degree and PhD degree in Electric power system and its automation from Northeast of China Electric Power University and North China Electric Power University in 1995 and 2009 respectively. She is currently associate professor in the school of Electrical and Electronic Engineering, North China Electric Power University, Beijing, China. Her research interests are risk assessment in power system, power system dispatch and management automation
Qi Chen He received B.S degree in electrical engineering from North China Electric Power University in 2013. He is currently a M.S. degree candidate in school of electrical engineering in North China Electric Power University. His major interest lies in the reliability evaluation in power system.
Yuying Zhang She received the B.S. degree in electrical engineering from the Southwest Jiaotong University, Chengdu, China, in 2012 and currently working towards the M.S. degree in electrical engineering at North China Electric Power University. Her research interests include active distribution system planning, investment evaluation of smart grid, and reliability theories.
Guobing Qiu He received B.S degree and M.S. degree in electrical engineering from North China Electric Power University in 2011 and 2014. He is currently working in Jinhua electric power corporation. His major interest lies in the reliability modeling in power system.
Chenghui Lin He was born in 1983. He received his M.S. degree in Sichuan University. He is working in Guizhou Electric Power Research Institute. His major interest lies in the electric power system and its automation.
Sperling K.
,
Hvellund F.
,
Mathiesen BV.
2010
“Evaluation of wind power planning in Denmarktowards an integrated perspective”
Energy
35
(12)
5443 
5454
DOI : 10.1016/j.energy.2010.06.039
Zhang L.
,
Qiu W.
2012
“Consideration of Largescale Wind Farm Disconnection from Grid.”
Electric Power Construction
33
(3)
11 
14
Castro Sayas F.
,
Allan R. N.
1996
“Generation availability assessment of wind farms”
IET
Generation, Transmission and Distribution, IEE Proceedings
143
(5)
Dobakhshari A.S.
,
FotuhiFiruzabad M.
2009
“A reliability model of large wind farms for power system adequacy studies”
Energy Conversion, IEEE Transactions on
24
(3)
792 
801
DOI : 10.1109/TEC.2009.2025332
Andréa P. Leite
,
Borges C. L.
,
Falcao D. M.
2006
“Probabilistic wind farms generation model for reliability studies applied to Brazilian sites”
IEEE Trans. Power Systems
21
(4)
1493 
1501
Billinton R.
,
Chen H.
,
Ghajar R.
1996
“A sequential simulation technique for adequacy evaluation of generating systems including wind energy”
Energy Conversion, IEEE Transactions on
11
(4)
728 
734
DOI : 10.1109/60.556371
Billinton R.
,
Chen H.
1998
“Assessment of riskbased capacity benefit factors associated with wind energy conversion systems”
Power Systems, IEEE Transactions on
13
(3)
1191 
1196
Karki R.
,
Hu P.
,
Billinton R.
2006
“A simplified wind power generation model for reliability evaluation”
Energy conversion, IEEE Transactions on
21
(2)
533 
540
DOI : 10.1109/TEC.2006.874233
Billinton R.
,
Gao Y.
2008
“Multistate wind energy conversion system models for adequacy assessment of generating systems incorporating wind energy”
Energy Conversion, IEEE Transactions on
23
(1)
163 
170
DOI : 10.1109/TEC.2006.882415
WANG H.
,
BAI X.
2012
“Assessment on generation adequacy and capacity credit for largescale wind farm”
Power System Technology
6
036 
Jensen N. O.
1983
A note on wind generator interaction
GonzálezLongatt F.
,
Wall P.
,
Terzija V.
2012
“Wake effect in wind farm performance: Steadystate and dynamic behavior”
Renewable Energy
39
(1)
329 
338
DOI : 10.1016/j.renene.2011.08.053
Carta J.A.
,
Ramírez P.
,
Velázquez S.
2009
“A review of wind speed probability distributions used in wind energy analysis: Case studies in the Canary Islands”
Renewable and Sustainable Energy Reviews
13
(5)
933 
955
DOI : 10.1016/j.rser.2008.05.005
Akdağ S. A.
,
Dinler A.
2009
“A new method to estimate Weibull parameters for wind energy applications”
Energy conversion and management
50
(7)
1761 
1766
DOI : 10.1016/j.enconman.2009.03.020
Hetzer John
,
Yu D. C.
,
Bhattarai K.
2008
“An economic dispatch model incorporating wind power”
Energy Conversion, IEEE Transactions on
23
(2)
603 
611
DOI : 10.1109/TEC.2007.914171
ZHENG R.
2008
“Modeling of LargeScale Wind Farms in the Probabilistic Power Flow Analysis Considering Wake Effects [J]”
Journal of Xi’an Jiaotong University
12
019 
de and Mikel P. G.
2012
“Power generation efficiency analysis of offshore wind farms connected to a SLPC (single large power converter) operated with variable frequencies considering wake effects”
Energy
37
(1)
455 
468
DOI : 10.1016/j.energy.2011.11.010
Wen J.
,
Zheng Y.
,
Donghan F.
2009
“A review on reliability assessment for wind power”
Renewable and Sustainable Energy Reviews
13
(9)
2485 
2494
DOI : 10.1016/j.rser.2009.06.006
Wu L. W.
,
Zhang J. H.
,
Liu R.X
2012
“A wind farm reliability model considering wind turbine faults and its application”
Automation of Electric Power Systems
16
(36)
31 
35
Jin J.
,
Ai Q.
,
Xi L.
2007
“Internal electrical wiring system of offshore wind farms”
East China Electric Power
35
(10)
20 
23
Kalyani S.
,
Swarup K. S.
2011
“Particle swarm optimization based Kmeans clustering approach for security assessment in power systems”
Expert Systems with Applications
38
(9)
10839 
10846
DOI : 10.1016/j.eswa.2011.02.086
Bie Z.
2009
“Studies on models and algorithms of the power system probabilistic production simulation integrated with wind farm”
IEEE
Power & Energy Society General Meeting, 2009. PES’09. IEEE
1 
7
2010
“Monthly and seasonal assessment of wind energy characteristics at four monitored locations in Liguria region (Italy)”
Renewable and Sustainable Energy Reviews
14
(7)
1959 
1968
DOI : 10.1016/j.rser.2010.04.015
Subcommittee P. M.
1979
“IEEE reliability test system”
Power Apparatus and Systems, IEEE Transactions on
6
2047 
2054