1. Introduction

JPEE

Journal of Power and Energy Engineering

2327-588X

Scientific Research Publishing

10.4236/jpee.2022.108003

JPEE-119195

Articles

Engineering

Reliability Analysis of Wind Energy Generation System Using Stochastic Method

Godwin

Diamenu

¹^*Joseph

Cudjoe Attachie

¹Christian

Kwaku Amuzuvi

Renewable Energy Engineering Department, University of Mines and Technology (UMaT), Tarkwa, Ghana

Electrical and Electronic Engineering Department, University of Mines and Technology (UMaT), Tarkwa, Ghana

08082022

1008264423, May 202213, August 2022 16, August 2022

2014

This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

Operators of renewable energy systems (RESs) must always manage uncertainty to some extent to ensure the reliability and the security of the electric power supply source. The guiding principle in this regard is to ensure service reliability and quality by balancing load variations with the variable renewable energy (VRE) sources. If the power generated by these VRE sources is not properly managed in conjunction with the varying load, the power grid may fail to achieve the required balance. To ensure its reliable operation, reliability analysis is vital for wind energy generation system (WEGS). This paper evaluated and assessed the reliability of WEGS and a proposed varying load by first using a stochastic approach to model the WEGS and the proposed varying load after which power generation indices were used to evaluate and assess the performance of the model. The WEGS and the varying load were modelled separately after which the two were combined into one model. Full availability, partial availability, the expected energy not supplied (EENS) or loss of energy expectation (LOEE), the mean or average instantaneous electric power generation and mean instantaneous generation deficiency were the indices used for the evaluation of the WEGS. The results indicated that the electric power generation will meet the power demand during most of the transition states of the WEGS with the expectation that the variation in the load will not be at fast pace and in large quantum.

Renewable Energy Stochastic Reliability Availability Transition Rate Transition State

1. Introduction

Power generation from low carbon energy resources has significantly increased in recent times. Wind energy is an important component of future energy systems to meet growing energy demands. The demand for wind energy is growing rapidly all over the world. According to the Global Wind Energy Council (GWEC), by the year 2020, 350 GW of wind power capacity would have been installed. According to European Wind Energy Association (EWEA), an increment of about 320 GW is expected in European wind power capacity by the year 2030 [1]. With the push to adopt more renewable energy resources (RERs), the total installed capacity of wind power throughout the world has exceeded 742.689 GW to date [2].

Fossil fuels are non-renewable and their associated prices keep fluctuating on the World. The increasing environmental and climatic conditions associated with fossil fuel usage had moved the research focus from conventional electricity resources to RERs. RERs, such as wind, solar energy from the sun, tidal wave, biomass from plants, small hydro sources and geothermal power, are clean alternatives to fossil fuels. Among them, wind energy is one of the most promising and rapidly developing RERs in the world today. The main attractions of wind energy are its naturally free availability everywhere and its low environmental impact [3]. Wind velocity determines the output of wind turbine generators (WTGs). Wind farms’ production output is stochastic and completely different from that of conventional power generating units due to the fast changing and unpredictable characteristics of wind velocity and the random nature of WTG failures. Wind power generation’s large variations are mostly caused by real-time changes in meteorological conditions [4].

Large-scale deployment of variable renewable energy (VRE) is associated with reliability and availability concerns. Wind energy generation system (WEGS) reliability and availability are vital and necessary requirements for meeting electric power demand [5]. Reliability is the ability of a system to perform its required function without failure, for a given time interval under given conditions, whereas availability is the ability to be in a state to perform as required [6]. To ensure its reliable operation, reliability analysis is vital for WEGS. In the last two decades, reliability evaluation methods for electric power systems, including VRE sources, have been developed.

The Markov model gives a simple description of a component, which is normally handled with mathematical methods. In order to be able to utilise it, the components of the system must be able to be described as a Markov model. This means that the system should be represented as a system lacking memory of previous states with identifiable system states. The system’s reliability can then be evaluated using techniques based on the principles of the K state Markov process. For this purpose, the state probabilities can only be obtained by solving K differential equations either analytically or numerically [4].

The reliability analysis techniques based on the Markov chain process may be classified into two time domain modes: discrete and continuous, and reliability analysis in these two modes can be performed using the system state transition diagram. The Markov process had been utilised to derive the probability value of each working state, the failure state, and the time required for the system to reach a steady state in reliability assessment over the years [7] [8].

A lot of reliability analysis of wind turbine (WT) system has been performed using WT components failure rates and restoration rates data in the literature. However, in the literature, there is no large-scale model for the study and evaluation of WEGS and a varying load. In this study, by means of different states and transitions between them, this paper is an insight into the performance analysis of WEGS and its accompanying varying load presented through a Markov model. The WEGS and the varying load were modelled separately. The two models were then combined into one model which was used to evaluate the generation indices. WEGS availability, expected energy not supplied (EENS) among others was performed using WT transition states due to the varying wind velocity.

2. The Proposed Modelling Approach

The Markov process approach is widely used to model the states of continuous and discrete stochastic processes of complex systems based on which the system reliability is ascertained. Markov process is a stochastic process evolving under the Markov property. The future state of the Markov process depends on the current state, it is conditionally independent of the past states of the stochastic process. Given a state X_k, the state X_k_+i, where i > 0 are independent of state X_k_-j where j > 0 [4] [8]. The Markov model is a method based on states and the transitions between them. There are two main assumptions that are associated with Markov model:

• The system which is modelled by Markov process approach is memory-less. Thus, the future probability of each event is a function of the current state and it is independent from the former transitions to reach the new state.

• The probability of transition between states is assumed to be constant.

A three state Markov chain system transition diagram is given in Figure 1. State ‘0’ represents the normal working state, and states “1” and “2” respectively represents partial and abnormal working states under different subsystem failure conditions.

For the three-state repairable system in Figure 1, the system state function is given as follows:

W ( t ) = { 0 , normal working state 1 , partial working state 2 , faulty working state

where, λ is failure rate of the system, μ is the repair rate of the system, Δt is the represents a very short time interval.

Based on the state transition diagram given in Figure 1, the state transition matrix equation in the time interval [t, t + Δt] for the system for each state represented by an r-by-r transition matrix based on the transitions between the states is derived as follows:

[ P 0 ( t + Δ t ) P 1 ( t + Δ t ) P 2 ( t + Δ t ) ] = [ 1 − ( λ 0 , 1 + λ 0 , 2 ) μ 1 , 0 μ 2 , 0 λ 0 , 1 1 − ( μ 1 , 0 + λ 2 , 1 ) μ 2 , 1 λ 0 , 2 λ 2 , 1 1 − ( μ 2 , 0 + μ 2 , 1 ) ] [ P 0 ( t ) P 1 ( t ) P 2 ( t ) ] (1)

By taking the limit as ∆t → 0 in Equation (1), the differential equations of the system state probabilities are obtained as follows:

{ d p 0 d t = limit Δ → 0 p 0 ( t + Δ t ) Δ t = 1 − ( λ 0 , 1 + λ 0 , 2 ) p 0 ( t ) + μ 1 , 0 p 1 ( t ) + μ 2 , 0 p 2 ( t ) d p 1 d t = limit Δ → 0 p 1 ( t + Δ t ) Δ t = λ 0 , 1 p 0 ( t ) + 1 − ( μ 1 , 0 + λ 2 , 1 ) p 1 ( t ) + μ 2 , 1 p 2 ( t ) d p 2 d t = limit Δ → 0 p 2 ( t + Δ t ) Δ t = λ 0 , 2 p 0 ( t ) + λ 2 , 1 p 1 ( t ) + 1 − ( μ 2 , 0 + μ 2 , 1 ) p 2 ( t ) (2)

where, the subscripts indicate the state the transition emanates from, and the next state it is going to terminate. Where, P₀ is the probability of the system being in normal working state, P₁ and P₂ is the system being in partial and abnormal working states respectively.

For the four-state repairable system in Figure 2, its state function is given as follows.

W ( t ) = { 0 , normal working state ( fullyavailable ) 1 , partial working state ( partially available ) 2 , degraded state ( somehow availabletoanextent ) 3 , faulty state ( not available )

Figure 2 is a four state transition diagram for a Markov chain process. The state transition matrix equation in the time interval [t, t + Δt] for the system for each state represented by an r-by-r transition matrix based on the transitions between the states is derived as follows:

[ P 0 ( t + Δ t ) P 1 ( t + Δ t ) P 2 ( t + Δ t ) P 3 ( t + Δ t ) ] = [ 1 − ( λ 0 , 1 + λ 0 , 2 ) μ 1 , 0 μ 2 , 0 0 λ 0 , 1 1 − ( μ 1 , 0 + λ 1 , 3 ) 0 μ 3 , 1 λ 0 , 2 0 1 − μ 2 , 0 0 0 λ 1 , 3 0 1 − μ 3 , 1 ] [ P 0 ( t ) P 1 ( t ) P 2 ( t ) P 3 ( t ) ] (3)

where, the subscripts indicate the state the transition emanates from, and the next state it is going to terminate, λ is the failure rate of the system, μ is the repair rate of the system, Δt is very short time interval, P₀ is the probability of the system being in normal working state, P₁ is the probability of the system being in partial working state, P₂ is the probability of the system being in degraded state, P₃ is the probability of the system being in abnormal working state.

By taking the limit as ∆t → 0 in Equation (3), the differential equations of the system state probabilities is obtained as follows:

{ d p 0 d t = limit Δ t → 0 p 0 ( t + Δ t ) Δ t = 1 − ( λ 0 , 1 + λ 0 , 2 ) p 0 ( t ) + μ 1 , 0 p 1 ( t ) + μ 2 , 0 p 2 ( t ) + 0 p 3 ( t ) d p 1 d t = limit Δ t → 0 p 1 ( t + Δ t ) Δ t = λ 0 , 1 p 0 ( t ) + 1 − ( μ 1 , 0 + λ 1 , 3 ) p 1 ( t ) + 0 p 2 ( t ) + μ 3 , 1 p 3 ( t ) d p 2 d t = limit Δ t → 0 p 2 ( t + Δ t ) Δ t = λ 0 , 2 p 0 ( t ) + 0 p 1 ( t ) + 1 − μ 2 , 0 p 2 ( t ) + 0 p 3 ( t ) d p 3 d t = limit Δ t → 0 p 3 ( t + Δ t ) Δ t = 0 p 0 ( t ) + λ 1 , 3 p 1 ( t ) + 0 p 2 ( t ) + 1 − μ 3 , 1 p 3 ( t ) (4)

3. Case Study

To validate the proposed model, a four-state WEGS model with a four-state load model is proposed. The proposed modelling framework and the reliability analysis procedure is given by Figure 3.

3.1. WEGS Model

The model for the WEGS is made up of wind farms comprising a number of WTs. It is a four-state and eight-transition states model where a certain amount

of electric power is generated at each state to feed the electric load. The states are based on the WT power output curve given in Figure 4, which establishes a relation between the wind speed at its hub height and the power that can be extracted from the wind and injected into the power grid. It is very useful for evaluating the WT performance [9]. Even though it has cut-in, rated and cut-out speeds, another speed range is introduced between the cut-in speed and the rated speed which is known as the mid-range speed resulting in the four-state transition states instead of the original three-state transition states. The transition state diagram based on the WT output power curve is given in Figure 5. Table 1 gives the wind speed and the assumed electric power generated during each state. The analysis assumed that no electric power is produced when the wind speed is below the cut-in speed.

The assumed transition rates among the various states of the WEGS is given by the transition rate matrix, M defined as:

M = ( m i j ) = ( 0.247 0.284 0.302 0.167 0.276 0.249 0.235 0.240 0.348 0.324 0.233 0.095 0.175 0.121 0.302 0.402 ) (5)

The assumed transition rate matrix is chosen in such a way that every row sum up to one which is in line with how the transition of the electric power demand occurs in the electric power and electric power demand combined model given by Figure 7, where an arrow indicates the direction of the electric power generation transition states (row wise) and electric power demand transition states (column wise) respectively.

3.2. Model for the Electric Power Demand

In reliability analysis, the model for the proposed load is very necessary since it indicates the electric power demand that must be met by the generation source.

Table 1 Electric power produced by the WEGS during each state

S/No.	Wind Speed	Transition States	Electric Power Generated (MW)
1	3.50 - 6.50 m/s	State 0	10
2	6.50 - 8.50 m/s	State 1	120
3	8.5 - 14.5 m/s	State 2	200
4	14.5 - 18.5 m/s	State 3	40

The transition can occur between two states at any particular time. An electric power system must match generation and load in real time. The electric power demand that must be met by the generation source fluctuates within certain time interval. Hence, in modelling the electric power demand or the load, uncertainty or randomness must be taken into consideration. The state transition diagram for the varying electric load is given in Figure 6 indicating four different states comprising different load capacities as given in Table 2. The electric power load is assumed to be varying between the stated minimum and maximum values.

Table 2 Varying power demand states for reliability analysis

S/No.	Transition States	Electric Power Demand (MW)
1	State 0	50
2	State 1	100
3	State 2	140
4	State 3	160

The assumed transition rates among the various states of the electric power demand is given by the transition rate matrix, N defined as:

N = ( n i j ) = ( 0.18 0.19 0.32 0.23 0.24 0.31 0.29 0.30 0.25 0.28 0.21 0.27 0.33 0.22 0.18 0.20 ) (6)

The assumed transition rate matrix is chosen in such a way that every column sum up to one which is in line with how the transition of the electric power demand occurs in the combined electric power generation and electric power demand model given by Figure 7.

3.3. Combined WEGS and Electric Power Demand Model

The objective of WEGS is to deliver the generation capacity to the load, it is possible to combine the equivalent transition diagram of the generation subsystem with the equivalent transition diagram of the electric power load subsystem. The states transition diagram for the WEGS and an electric power demand combined is given in Figure 7.

For a combined WEGS and the electric power demand, the output of the WEGS is represented by a number of an electric power generation states and corresponding transition rate matrix is given by M = [m_ij]. That of the electric power demand is represented by b number of transition states with the corresponding transition rate matrix given by N = [n_ij]. With the number of transition states for the WEGS being a, and that of the electric power demand being b, the total number of transition states for the WEGS and that of the electric power load combined is a x b transition states [10] [11]. The transition states for the

WEGS is given by g ∈ ( g 0 , g 1 , g 2 , g 3 ) and that of the electric power demand is also given by d ∈ ( d 0 , d 1 , d 2 , d 3 ) . For this case study, the number of transition states of the WEGS and that of the electric power demand together gives sixteen (16) transition states. The state transition diagram for the combined WEGS and the electric power demand model is given in Figure 7. The horizontal transition states represent the power generation states of the wind farms in increasing order from left to right columns and the vertical transition states represents the demand states of the load. It is a state-transition-rate diagram as opposed to a state-transition diagram because it shows the rate at which the process moves from state to state and not the probability of moving from one state to another [12]. More than one transition cannot occur at a time. The combined electric power generation and electric power demand model has 16 transition states, which are categorized into two different classes, states (1, 2, 5, 6, 9, 10, 13, 14) that fully meets the electric power demand irrespective of the state in which it is, and those states (0, 3, 4, 7, 8, 11, 12. 15) that partially meets the electric power demand. The power generation and power demand at the various transition states is given in Table 3.

A graph of electric power generation and power demand at the various transition states against the various probabilities at those transition states is given by Figure 8.

The probabilities, p_i(t), 0 ≤ i ≤ 15 of the combined WEGS and electric power demand of the WEGS at the various states at the time, t + ∆ are:

Table 3 Power generation and power demand at the various transition states

States	Probabilities, P(t)	Electric Power Generation (MW)	Electric Power Demand (MW)
0	p 0 ( t )	10	50
1	p 1 ( t )	10	100
2	p 2 ( t )	10	140
3	p 3 ( t )	10	160
4	p 4 ( t )	120	50
5	p 5 ( t )	120	100
6	p 6 ( t )	120	140
7	p 7 ( t )	120	160
8	p 8 ( t )	200	50
9	p 9 ( t )	200	100
10	p 10 ( t )	200	140
11	p 11 ( t )	200	160
12	p 12 ( t )	40	50
13	p 13 ( t )	40	100
14	p 14 ( t )	40	140
15	p 15 ( t )	40	160

[ P 0 ( t + Δ t ) ] = [ ( 1 − g 01 + g 02 + g 03 + d 01 + d 02 + d 03 ) Δ t P 0 ( t ) + g 10 Δ t P 1 ( t ) + g 20 Δ t P 2 ( t ) + g 30 Δ t P 3 ( t ) + d 10 Δ t P 4 ( t ) + d 20 Δ t P 8 ( t ) + d 30 Δ t P 12 ( t ) ] [ P 0 ( t ) ] (7)

[ P 1 ( t + Δ t ) ] = [ g 01 Δ t P 0 ( t ) + ( 1 − g 10 + g 12 + g 13 + d 01 + d 02 + d 03 ) Δ t P 1 ( t ) + g 21 Δ P 2 ( t ) + g 31 Δ t P 3 ( t ) + d 10 Δ t P 5 ( t ) + d 20 Δ t P 9 ( t ) + d 30 Δ t P 13 ( t ) ] [ P 1 ( t ) ] (8)

[ P 2 ( t + Δ t ) ] = [ g 12 Δ t P 0 ( t ) + g 02 P 1 ( t ) Δ t + ( 1 − g 20 + g 21 + g 23 + d 01 + d 02 + d 03 ) Δ t P 2 ( t ) + g 32 Δ t P 3 ( t ) + d 10 Δ t P 6 ( t ) + d 20 Δ t P 10 ( t ) + d 30 Δ t P 14 ( t ) ] [ P 2 ( t ) ] (9)

[ P 3 ( t + Δ t ) ] = [ g 03 Δ t P 0 ( t ) + g 13 Δ t P 1 ( t ) + g 13 Δ t P 2 ( t ) ( 1 − g 30 + g 31 + g 32 + d 01 + d 02 + d 03 ) Δ t P 3 ( t ) + P 7 ( t ) d 10 Δ t + P 11 ( t ) d 02 Δ t + P 15 ( t ) d 30 Δ t ] [ P 3 ( t ) ] (10)

[ P 4 ( t + Δ t ) ] = [ ( 1 − g 01 + g 02 + g 03 + d 10 + d 12 + d 13 ) Δ t P 4 ( t ) + g 10 Δ t P 5 ( t ) + g 20 Δ t P 6 ( t ) + g 30 Δ t P 7 ( t ) + d 01 Δ t P 0 ( t ) + d 21 Δ t P 8 ( t ) + d 31 Δ t P 12 ( t ) ] [ P 4 ( t ) ] (11)

[ P 5 ( t + Δ t ) ] = [ ( 1 − g 10 + g 12 + g 13 + d 10 + d 12 + d 13 ) Δ t P 5 ( t ) + g 01 Δ t P 4 ( t ) + g 21 Δ t P 6 ( t ) + g 31 Δ t P 7 ( t ) + d 01 Δ t P 1 ( t ) + d 21 Δ t P 9 ( t ) + d 31 Δ t P 13 ( t ) ] [ P 5 ( t ) ] (12)

[ P 6 ( t + Δ t ) ] = [ ( 1 − g 20 + g 21 + g 23 + d 10 + d 12 + d 13 ) Δ t P 6 ( t ) + g 02 Δ t P 4 ( t ) + g 12 Δ t P 5 ( t ) + g 32 Δ t P 7 ( t ) + d 01 Δ t P 2 ( t ) + d 21 Δ t P 10 ( t ) + d 31 Δ t P 14 ( t ) ] [ P 6 ( t ) ] (13)

[ P 7 ( t + Δ t ) ] = [ ( 1 − g 30 + g 31 + g 32 + d 10 + d 12 + d 13 ) Δ t P 7 ( t ) + g 03 Δ t P 4 ( t ) + g 13 Δ t P 5 ( t ) + g 23 Δ t P 6 ( t ) + d 01 Δ t P 3 ( t ) + d 21 Δ t P 11 ( t ) + d 31 Δ t P 13 ( t ) ] [ P 7 ( t ) ] (14)

[ P 8 ( t + Δ t ) ] = [ ( 1 − g 01 + g 02 + g 03 + d 20 + d 21 + d 23 ) Δ t P 8 ( t ) + g 10 Δ t P 9 ( t ) + g 20 Δ t P 10 ( t ) + g 30 Δ t P 11 ( t ) + d 02 Δ t P 0 ( t ) + d 12 Δ t P 4 ( t ) + d 32 Δ t P 12 ( t ) ] [ P 8 ( t ) ] (15)

[ P 9 ( t + Δ t ) ] = [ g 01 Δ t P 8 ( t ) + ( 1 − g 10 + g 12 + g 13 + d 20 + d 23 + d 21 ) Δ t P 9 ( t ) + g 21 Δ t P 10 ( t ) + g 31 Δ t P 11 ( t ) + d 02 Δ t P 1 ( t ) + d 12 Δ t P 5 ( t ) + d 32 Δ t P 13 ( t ) ] [ P 9 ( t ) ] (16)

[ P 10 ( t + Δ t ) ] = [ g 02 Δ t P 8 ( t ) + g 12 Δ t P 9 ( t ) + ( 1 − g 20 + g 21 + g 23 + d 20 + d 21 + d 23 ) Δ t P 10 ( t ) + g 32 Δ t P 11 ( t ) + d 02 Δ t P 2 ( t ) + d 12 Δ t P 6 ( t ) + d 32 Δ t P 14 ( t ) ] [ P 10 ( t ) ] (17)

[ P 11 ( t + Δ t ) ] = [ g 03 Δ t P 8 ( t ) + g 13 Δ t P 9 ( t ) + g 23 Δ t P 10 ( t ) + ( 1 − g 30 + g 31 + g 32 + d 10 + d 21 + d 23 ) Δ t P 11 ( t ) + d 01 Δ t P 3 ( t ) + d 12 Δ t P 7 ( t ) + d 32 Δ t P 15 ( t ) ] [ P 11 ( t ) ] (18)

[ P 12 ( t + Δ t ) ] = [ ( 1 − g 01 + g 02 + g 03 + d 30 + d 31 + d 32 ) Δ t P 12 ( t ) + P 13 ( t ) g 10 Δ t + P 14 ( t ) g 20 Δ t + g 30 Δ t P 15 ( t ) + d 03 Δ t P 0 ( t ) + d 13 Δ t P 4 ( t ) + d 23 Δ t P 8 ( t ) ] [ P 12 ( t ) ] (19)

[ P 13 ( t + Δ t ) ] = [ g 01 Δ t P 12 ( t ) + ( 1 − g 10 + g 12 + g 13 + d 30 + d 31 + d 32 ) Δ t P 13 ( t ) + g 21 Δ t P 14 ( t ) + g 31 Δ t P 15 ( t ) + d 03 Δ t P 1 ( t ) + d 13 Δ t P 5 ( t ) + d 23 Δ t P 9 ( t ) ] [ P 13 ( t ) ] (20)

[ P 14 ( t + Δ t ) ] = [ g 02 Δ t P 12 ( t ) + g 12 Δ t P 13 ( t ) + ( 1 − g 20 + g 21 + g 23 + d 30 + d 31 + d 32 ) Δ t P 14 ( t ) + P 15 ( t ) g 32 Δ t + d 03 Δ t P 2 ( t ) + d 13 Δ t P 6 ( t ) + d 23 Δ t P 10 ( t ) ] [ P 14 ( t ) ] (21)

[ P 15 ( t + Δ t ) ] = [ g 03 Δ t P 12 ( t ) + g 13 Δ t P 13 ( t ) + g 23 Δ t P 14 ( t ) + ( 1 − g 30 + g 31 + g 32 + d 30 + d 31 + d 32 ) Δ t P 15 ( t ) + d 03 Δ t P 3 ( t ) + d 13 Δ t P 7 ( t ) + d 23 Δ t P 11 ( t ) ] [ P 15 ( t ) ] (22)

Taking the limit as ∆t → 0 in Equation (7) to Equation (22), and by using the Laplace transform as well as the stated assumptions for each, Equation (7) to Equation (22) for the state probabilities, p_i(t), 0 ≤ i ≤ 15 are obtained as follows:

For Equation (7), P₍₀₎ = 1, P₍₁₎ = P₍₂₎ = P₍₃₎ = P₍₄₎ = P₍₈₎ = P₍₁₂₎ = 0;

P 0 ( t ) = e − ( g 01 + g 02 + g 03 + d 01 + d 02 + d 03 ) t (23)

For Equation (8), P₍₁₎ = 1, P₍₀₎ = P₍₂₎ = P₍₃₎ = P₍₅₎ = P₍₉₎ = P₍₁₃₎ = 0;

P 1 ( t ) = e − ( g 10 + g 12 + g 13 + d 01 + d 02 + d 03 ) t (24)

For Equation (9), P₍₂₎ = 1, P₍₀₎ = P₍₁₎ = P₍₃₎ = P₍₆₎ = P₍₁₀₎ = P₍₁₄₎ = 0;

P 2 ( t ) = e − ( g 20 + g 21 + g 23 + d 01 + d 02 + d 03 ) t (25)

For Equation (10), P₍₃₎ = 1, P₍₀₎ = P₍₁₎ = P₍₂₎ = P₍₇₎ = P₍₁₁₎ = P₍₁₅₎ = 0;

P 3 ( t ) = e − ( g 30 + g 31 + g 32 + d 01 + d 02 + d 03 ) t (26)

For Equation (11), P₍₄₎ = 1, P₍₀₎ = P₍₅₎ = P₍₆₎ = P₍₇₎ = P₍₈₎ = P₍₁₂₎ = 0;

P 4 ( t ) = e − ( g 01 + g 02 + g 03 + d 10 + d 12 + d 13 ) t (27)

For Equation (12), P₍₅₎ = 1, P₍₁₎ = P₍₄₎ = P₍₆₎ = P₍₇₎ = P₍₉₎ = P₍₁₃₎ = 0;

P 5 ( t ) = e − ( g 10 + g 12 + g 13 + d 10 + d 12 + d 13 ) t (28)

For Equation (13), P₍₆₎ = 1, P₍₂₎ = P₍₄₎ = P₍₅₎ = P₍₇₎ = P₍₁₀₎ = P₍₁₄₎ = 0;

P 6 ( t ) = e − ( g 20 + g 21 + g 23 + d 10 + d 12 + d 13 ) t (29)

For Equation (14), P₍₇₎ = 1, P₍₃₎ = P₍₄₎ = P₍₅₎ = P₍₆₎ = P₍₁₁₎ = P₍₁₃₎ = 0;

P 7 ( t ) = e − ( g 30 + g 31 + g 32 + d 10 + d 12 + d 13 ) t (30)

For Equation (15), P₍₈₎ = 1, P₍₀₎ = P₍₄₎ = P₍₉₎ = P₍₁₀₎ = P₍₁₁₎ = P₍₁₂₎ = 0;

P 8 ( t ) = e − ( g 01 + g 02 + g 03 + d 20 + d 21 + d 23 ) t (31)

For Equation (16), P₍₉₎ = 1, P₍₁₎ = P₍₅₎ = P₍₈₎ = P₍₁₀₎ = P₍₁₁₎ = P₍₁₃₎ = 0;

P 9 ( t ) = e − ( g 10 + g 12 + g 13 + d 20 + d 23 + d 21 ) t (32)

For Equation (17), P₍₁₀₎ = 1, P₍₂₎ = P₍₆₎ = P₍₈₎ = P₍₉₎ = P₍₁₁₎ = P₍₁₄₎ = 0;

P 10 ( t ) = e − ( g 20 + g 21 + g 23 + d 20 + d 21 + d 23 ) t (33)

For Equation (18), P₍₁₁₎ = 1, P₍₃₎ = P₍₇₎ = P₍₈₎ = P₍₉₎ = P₍₁₀₎ = P₍₁₅₎ = 0;

P 11 ( t ) = e − ( g 30 + g 31 + g 32 + d 10 + d 21 + d 23 ) t (34)

For Equation (19), P₍₁₂₎ = 1, P₍₀₎ = P₍₄₎ = P₍₈₎ = P₍₁₃₎ = P₍₁₄₎ = P₍₁₅₎ = 0;

P 12 ( t ) = e − ( g 01 + g 02 + g 03 + d 30 + d 31 + d 32 ) t (35)

For Equation (20), P₍₁₃₎ = 1, P₍₁₎ = P₍₅₎ = P₍₉₎ = P₍₁₂₎ = P₍₁₄₎ = P₍₁₅₎ = 0;

P 13 ( t ) = e − ( g 10 + g 12 + g 13 + d 30 + d 31 + d 32 ) t (36)

For Equation (21), P₍₁₄₎ = 1, P₍₂₎ = P₍₆₎ = P₍₁₀₎ = P₍₁₂₎ = P₍₁₃₎ = P₍₁₅₎ = 0;

P 14 ( t ) = e − ( g 20 + g 21 + g 23 + d 30 + d 31 + d 32 ) t (37)

For Equation (22), P₍₁₅₎ = 1, P₍₃₎ = P₍₇₎ = P₍₁₁₎ = P₍₁₂₎ = P₍₁₃₎ = P₍₁₄₎ = 0;

P 15 ( t ) = e − ( g 30 + g 31 + g 32 + d 30 + d 31 + d 32 ) t (38)

The availability, A(t) of the WEGS at any time, t is the sum of all the probabilities, thus, Equation (23) to Equation (38) which is given by Equation (39) [8] [13]:

A ( t ) = ∑ i = 1 n P i ( t ) , n = 16 (39)

Substituting the probabilities into Equation (39), results in Equation (40) as follows:

A ( t ) = e − ( g 01 + g 02 + g 03 + d 01 + d 02 + d 03 ) t + e − ( g 10 + g 12 + g 13 + d 01 + d 02 + d 03 ) t + e − ( g 30 + g 31 + g 32 + d 01 + d 02 + d 03 ) t + e − ( g 01 + g 02 + g 03 + d 10 + d 12 + d 13 ) t + e − ( g 10 + g 12 + g 13 + d 10 + d 12 + d 13 ) t + e − ( g 20 + g 21 + g 23 + d 10 + d 12 + d 13 ) t + e − ( g 30 + g 31 + g 32 + d 10 + d 12 + d 13 ) t + e − ( g 01 + g 02 + g 03 + d 20 + d 21 + d 23 ) t + e − ( g 10 + g 12 + g 13 + d 20 + d 23 + d 21 ) t + e − ( g 20 + g 21 + g 23 + d 20 + d 21 + d 23 ) t + e − ( g 30 + g 31 + g 32 + d 10 + d 21 + d 23 ) t + e − ( g 01 + g 02 + g 03 + d 30 + d 31 + d 32 ) t + e − ( g 10 + g 12 + g 13 + d 30 + d 31 + d 32 ) t + e − ( g 20 + g 21 + g 23 + d 30 + d 31 + d 32 ) t + e − ( g 30 + g 31 + g 32 + d 30 + d 31 + d 32 ) t (40)

The full availability, A f ( t ) of the WEGS during which the electric power demand is fully met occurs at states 1, 2, 5, 6, 9, 10, 13, and 14. Therefore, the full availability, A f ( t ) is given by Equation (41):

A f ( t ) = P 1 + P 2 + P 5 + P 6 + P 9 + P 10 + P 13 + P 14 = e − ( g 01 + g 02 + g 03 + d 01 + d 02 + d 03 ) t + e − ( g 10 + g 12 + g 13 + d 01 + d 02 + d 03 ) t + e − ( g 10 + g 12 + g 13 + d 10 + d 12 + d 13 ) t + e − ( g 20 + g 21 + g 23 + d 10 + d 12 + d 13 ) t + e − ( g 10 + g 12 + g 13 + d 20 + d 23 + d 21 ) t + e − ( g 20 + g 21 + g 23 + d 20 + d 21 + d 23 ) t + e − ( g 10 + g 12 + g 13 + d 30 + d 31 + d 32 ) t + e − ( g 20 + g 21 + g 23 + d 30 + d 31 + d 32 ) t (41)

The partial availability, A p ( t ) of the WEGS during which the electric power demand is not fully met occurs at states 0, 3, 4, 7, 8, 11, 12, and 15. Therefore, the partial availability, A p ( t ) is given by Equation (42):

A p ( t ) = P 0 + P 3 + P 4 + P 7 + P 8 + P 11 + P 12 + P 15 = e − ( g 01 + g 02 + g 03 + d 01 + d 02 + d 03 ) t + e − ( g 30 + g 31 + g 32 + d 01 + d 02 + d 03 ) t + e − ( g 01 + g 02 + g 03 + d 10 + d 12 + d 13 ) t + e − ( g 30 + g 31 + g 32 + d 10 + d 12 + d 13 ) t + e − ( g 01 + g 02 + g 03 + d 20 + d 21 + d 23 ) t + e − ( g 30 + g 31 + g 32 + d 10 + d 21 + d 23 ) t + e − ( g 01 + g 02 + g 03 + d 30 + d 31 + d 32 ) t + e − ( g 30 + g 31 + g 32 + d 30 + d 31 + d 32 ) t (42)

The assumed transition rates were used to compute the possible probabilities at the various transition states during an interval of one hour up to the fourth hour which is given by Table 4 using Equation (23) to Equation (39). A transition state probability graph based on Table 4 for the WEGS is given by Figure 9.

Table 4 Probabiliites of the combined WEGS and the power demands at the various states and at different times

States	Probabilities, P(t)	Time, T₁ = 1	Time, T₂ = 2	Time, T₃ = 3	Time, T₄ = 4
0	p 0 ( t )	0.224	0.449	0.673	0.898
1	p 1 ( t )	0.225	0.450	0.675	0.900
2	p 2 ( t )	0.221	0.443	0.664	0.885
3	p 3 ( t )	0.262	0.524	0.786	1.049
4	p 4 ( t )	0.205	0.410	0.615	0.821
5	p 5 ( t )	0.206	0.411	0.617	0.822
6	p 6 ( t )	0.202	0.405	0.607	0.809
7	p 7 ( t )	0.240	0.479	0.719	0.958
8	p 8 ( t )	0.211	0.423	0.634	0.846
9	p 9 ( t )	0.212	0.424	0.635	0.847
10	p 10 ( t )	0.208	0.417	0.625	0.834
11	p 11 ( t )	0.249	0.499	0.748	0.997
12	p 12 ( t )	0.227	0.453	0.680	0.907
13	p 13 ( t )	0.227	0.454	0.682	0.909
14	p 14 ( t )	0.224	0.447	0.671	0.894
15	p 15 ( t )	0.265	0.530	0.794	1.059

The mean or average instantaneous electric power generation, G_t by the WEGS at time, t based on Table 3 is given as follows [14] [15]:

G t = ∑ i = 0 15 g i p i ( t ) = 10 p 0 ( t ) + 10 p 1 ( t ) + 10 p 2 ( t ) + 10 p 3 ( t ) + 120 p 4 ( t ) + 120 p 5 ( t ) + 120 p 6 ( t ) + 120 p 7 ( t ) + 200 p 8 ( t ) + 200 p 9 ( t ) + 200 p 10 ( t ) + 200 p 11 ( t ) + 40 p 12 ( t ) + 40 p 13 ( t ) + 40 p 14 ( t ) + 70 p 15 ( t ) (43)

Assuming the average electric power demand is 90 MW (which is the average of the four electric power demands at the various transition states), the WEGS mean instantaneous generation deficiency, G_def is obtained as follows [14] [15]:

G d e f = ∑ i = 0 3 p i ( t ) max ( d a v − g i , 0 ) = 80 p 0 ( t ) + 30 p 1 ( t ) + 50 p 3 ( t ) (44)

where, G_def is the generation deficiency, d_av is average electric power demand, p_i(t) is the probability at a particular state, g_i(t) is the electric power generation at a particular state. A graph of mean instantaneous generation G_t and mean instantaneous generation deficiency, G_def is given by Figure 10. The expected energy not supplied (EENS) or loss of energy expectation (LOEE) also known expected unserved energy (EUE) is a measure of the amount of electricity demand (in MWh in a given year) expected to be lost when demand exceeds the available generation. Therefore, the expected energy not supplied (EENS) or expected unserved energy (EUE) to the consumers during the 24-hour operation of the WEGS is as follows [14] [15]:

EENS = ∫ 0 24 D d e f d t = ∫ 0 24 80 d t + 30 d t + 50 d t = 3840000 kWh (45)

The annual expected energy not supplied (EENS) is computed as follows [14] [15]:

EENS = ∫ 0 8760 D d e f d t = ∫ 0 8760 80 d t + 30 d t + 50 d t = 1401600 kWh (46)

The surplus electric power generation, G_surplus by the WEGS at time, t is given by [15]:

G s u r p l u s = ∑ i = 0 3 p i ( t ) max ( g i − d a v , 0 ) = 110 p 2 ( t ) (47)

G s u r p l u s = ∫ 0 24 110 d t = 2640000 kWh (48)

Based on Table 1 and Table 2 for the assumed power generation and load demand at the various transition states, the expected demand not supplied (EDNS) which is the difference between the expected power to be generated and the expected power demand by the consumers, is obtained at each transition state as given in Table 5 [16]. Graphs based on the various transition states are depicted by Figures 11-14 respectively.

Table 5 Power generation, power demand and EDNS

States	Generation	Demand	EDNS
0	10	50	40
	120	100	x
	200	140	x
	40	160	120
1	10	100	90
	120	50	x
	200	160	x
	40	140	100
2	10	140	130
	120	160	40
	200	50	x
	40	100	60
3	10	160	150
	120	140	20
	200	100	x
	40	50	10

4. Conclusion

Reliability analysis of WEGS and a proposed varying load was carried out using various electric power generation related indices to ascertain how a large-scale wind farm comprising WTs will behave under varying wind speeds and varying load conditions. It was obvious from the analysis that during certain transition states, the electric power demand will be more than that of the electric power generation. And during certain transition states, the electric power generation will exceed that of the power demand. It was also clear from the analysis that the probability that the WEGS system will be in a particular state is fairly high indicating that the electric power demand during that particular state is likely to be met. On the other hand, some transition states probabilities which is an indication that the electric power load during those transition states might not be met. The mean instantaneous generation graph and mean generation deficiency graph have shown that during the early hours of WEGS operation, there will be some generation deficiencies. However, as time progresses, the mean instantaneous power generation will exceed the mean electric power demand. The uniqueness of the study is that it brought to bear the behavior of a large-scale WEGS under different wind speeds, and a varying load conditions it is supplying. Due to the volatility of the price of fossil fuel and instability around the globe, this study is intended to examine how a future WEGS and other renewable energy resources for power generation can be operated on a large scale in the same manner like that of the conventional power generation system without excessive instability situations. A further study is proposed to ascertain how the model will behave when the model assumptions are varied and the transition rates as well.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Diamenu, G., Attachie, J.C. and Amuzuvi, C.K. (2022) Reliability Analysis of Wind Energy Generation System Using Stochastic Method. Journal of Power and Energy Engineering, 10, 26-44. https://doi.org/10.4236/jpee.2022.108003

References1

Kim, D. and Hur, J. (2017) Stochastic Prediction of Wind Generating Resources Using the Enhanced Ensemble Model for Jeju Island’s Wind Farms in South Korea. Sustainability, 9, Article 817. https://doi.org/10.3390/su9050817

Su, H. and Chen, L. (2021) Stochastic Model Analysis of Preventive Maintenance based on Wind Turbine Gearbox Space-Time Evolution. Journal of Applied Science and Engineering, 255, 721-732.

Zhu, C. and Li, Y. (2018) Stability Control and Reliable Performance of Wind Turbines. IntecOpen, London, 169-175.

Ding, Y., Singh, C., Goel, L., Ostergaard, J. and Wang, P. (2014) Short-Term and Medium-Term Reliability Evaluation for Power Systems with High Penetration of Wind Power. IEEE Transactions on Sustainable Energy, 5, 896-906. https://doi.org/10.1109/TSTE.2014.2313017

Yusuf, I., Ismail, A L. and Ali, U.A. (2020) Reliability Analysis of Communication Network System with Redundant Relay Station under Partial and Complete Failure. Journal of Mathematics and Computer Science, 10, 863-880.

Ozturk, S., Fthenakis, V.M. and Faulstich, S. (2018) Failure Modes, Effects and Criticality Analysis for Wind Turbines Considering Climatic Regions and Comparing Geared and Direct Drive Wind Turbines. Energies, 11, Article 2317. https://doi.org/10.3390/en11092317

Hou, K., Jia, H., Xu, X., Liu, Z. and Jiang, Y. (2016) A Continuous Time Markov Chain Based Sequential Analytical Approach for Composite Power System Reliability Assessment. IEEE Transactions on Power Systems, 31, 738-748, https://doi.org/10.1109/TPWRS.2015.2392103

Zuo, W. and Li, K. (2021) Three-state Markov Chain Based Reliability Analysis of Complex Traction Power Supply Systems. 5th IEEE International Conference on System Reliability and Safety, Palermo, 24-26 November 2021, 74-79. https://doi.org/10.1109/ICSRS53853.2021.9660623

Xiao, Z., Zhao, Q., Yang, X. and Zhu, A. (2020) A Power Performance Online Assessment Method of a Wind Turbine Based on the Probabilistic Area Metric. Applied Sciences, 10, Article 3268. https://doi.org/10.3390/app10093268

Roy, A. and Chatterjee, K. (2018) Reliability Analysis of a Multi-State Wind Farm Using Markov Process. Safety and Reliability, 37, 3-24. https://doi.org/10.1080/09617353.2017.1418811

Singh, C., Jirutitijaroen, P. and Mitra, J. (2019) Electric Power Grid Reliability Evaluation—Models and Methods. John Wiley and Sons, Inc., Hoboken, New Jersey, USA, 117-160. https://doi.org/10.1002/9781119536772

Ibe, O. (2013) Markov Processes for Stochastic Modelling. 2nd Edition,&#32 Elsevier, London. https://doi.org/10.1016/B978-0-12-407795-9.00015-3

Kalaiarasi, S., Anita, A.M. and Geethanjalii, R. (2017) Analysis of System Reliability Using Markov Technique. Global Journal of Pure and Applied Mathematics, 13, 5265-5273.

Levitina, G., Xing, L. and Huang, H.Z. (2019) Dynamic Availability and Performance Deficiency of Common Bus Systems with Imperfectly Repairable Components. Reliability Engineering and System Safety, 189, 58-66. https://doi.org/10.1016/j.ress.2019.04.007

Lisnianskia, A., Frenkel, I. and Khvatskin, L. (2017) On Sensitivity Analysis of Aging Multi-State System by Using LZ-Transform. Reliability Engineering and System Safety, 166, 99-108. https://doi.org/10.1016/j.ress.2016.12.001

Al-Shaalan, A.M. (2019) Reliability Evaluation of Power Systems. IntecOpen, London, United Kingdom. https://doi.org/10.5772/intechopen.85571

States	Generation	Demand	EDNS
0	10	50	40
	120	100	x
	200	140	x
	40	160	120
1	10	100	90
	120	50	x
	200	160	x
	40	140	100
2	10	140	130
	120	160	40
	200	50	x
	40	100	60
3	10	160	150
	120	140	20
	200	100	x
	40	50	10

States	Generation	Demand	EDNS
0	10	50	40
	120	100	x
	200	140	x
	40	160	120
1	10	100	90
	120	50	x
	200	160	x
	40	140	100
2	10	140	130
	120	160	40
	200	50	x
	40	100	60
3	10	160	150
	120	140	20
	200	100	x
	40	50	10

States	Generation	Demand	EDNS
0	10	50	40
	120	100	x
	200	140	x
	40	160	120
1	10	100	90
	120	50	x
	200	160	x
	40	140	100
2	10	140	130
	120	160	40
	200	50	x
	40	100	60
3	10	160	150
	120	140	20
	200	100	x
	40	50	10