InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » "System Reliability", book edited by Constantin Volosencu, ISBN 978-953-51-3706-1, Print ISBN 978-953-51-3705-4, Published: December 20, 2017 under CC BY 3.0 license. © The Author(s).

# Power System Reliability: Mathematical Models and Applications

By Rabah Medjoudj, Hassiba Bediaf and Djamil Aissani
DOI: 10.5772/intechopen.71926

Article top

## Overview

Figure 1. Superimposition of the system available capacity and the load model.

Figure 2. Three states diagram.

Figure 3. Diagram of transition states of a system subject to four failure processes.

Figure 4. Three-dimensional matrix of system’ states, where I 1 ∈ Ω 1 , J 2 ∈ Ω 2 , K 3 ∈ Ω 3 , and L ∈ Ω .

Figure 5. The three-dimensional matrix of the studied system.

Figure 6. Probabilities of the states and the system reliability changing.

# Power System Reliability: Mathematical Models and Applications

Rabah Medjoudj1, Hassiba Bediaf2 and Djamil Aissani3
Show details

## Abstract

This chapter deals with power systems reliability including technical, economical, and decisional aspects. Knowing that almost 90% of failures occur in the distribution systems, great interest was dedicated to this part of the system, and the first work was oriented to reliability indices defined as objectives to attempt and as performance measures in the electricity market. Some works deal with the managers’ behavior, and the customers reactions are modeled using economic criteria in uncertain future and inspired from game theory. When studying components, degradation models were introduced and combined with the effects of socks to study the reliability changing during system operation. In some works, the correlation between maintenance policies and reliability aspects was highlighted. In a recent work, considering the importance of new technologies integration and renewable energy insertion to power systems, it was revealed that reliability aspects and energy sustainability are two fundamental issues of progress in a given society.

Keywords: power systems reliability, distribution functions, degradation modeling, maintenance, decision-making

## 1. Introduction

In general way, power system reliability addresses the issues of service interruption and power supply loss. In several cases, it is defined as an objective to attempt in terms of indices directly related to the customer. Typical reliability index values for US utilities are SAIFI, SAIDI, and CAIDI. Over time, they become standard values for evaluating the reliability of electrical systems and used in several publications. Medjoudj et al. [1], in their recent publication, defined other indices as reliability subcriteria in their decision-making attributes, giving the best model of a smart energy grid. These indices are discussed in Section 2 with an application to a real case study. After this classical definition, some works have integrated data analysis and processing taking into account the calculation of distribution parameters such as those used for Weibull and Weibull-Markov processes. These works were vulgarized with applications in power system reliability by Van Casteren et al. [2] and Medjoudj et al. [3]. The most techniques used in power system reliability optimization and processing is the reliability centered preventive maintenance (RCM). Several publications have highlighted that in most cases of multicomponent systems, the maintenance actions arrive very early without any effects on the system or very late, that is, the need of curative maintenance with its negative consequences. This issue is also treated in the case of a differentiated service of reliability in the case of customers with different requirements of the reliability level. These two concepts are discussed in Section 3. Recent publications have highlighted the interest of combining reliability attributes for maintenance actions in the case of degrading systems and components [4]; however, a novel work developed recently at the LaMOS unit research dealing with multiple degradation processes is applied to power switchgear and is discussed in Section 4. The discussions and the conclusion highlighting the place of reliability in a power energy smart grid are given in Section 5.

## 2. Power system reliability indices

In the beginning, the methods used were classical to evaluate reliability indices of distribution systems such as failure frequency, mean failure times, mean time between failure and energy not supplied. These indices help decision makers to define technical and management measures to perform systems. After that was introduced the notion of loss of load probability (LOLP), which has a lot of applications in load modeling and electrical parameters dimensioning. It is significant for any power enterprise to analyze customer satisfaction. A variety of indices have been developed to measure reliability and its cost in power systems area such as loss of load probability (LOLP), loss of load expectation (LOLE), expected frequency of load curtailment (EFLC), expected duration of load curtailment (EDLC), expected duration of a curtailment (EDC), and expected energy not supplied (EENS) [1].

### 2.1. Loss of load probability

LOLP is an expected value, sometimes calculated on the basis of the peak hourly load of each day and sometimes on each hour’s load (24 in a day). Moreover, in the beginning, LOLP is used to characterize the adequacy of generation to serve the load on the bulk power system; it does not directly model the reliability of the transmission and distribution system where the majority of outages actually occur [5]. Nourelfeth and Ait Kadi [6] have recently made that the LOLP is usually used to estimate the reliability index. Considering S and D as the supply and the load demand respectively, they compute the reliability of a multistate system (MSS) as:

 R=PrS≥DOrR=1−LOLP (1)

Using the well-known formulation of the LOLP given in several publications and discussed in the following section, they generalized the MSS reliability index R

 R=1∑j=1MTj∑j=1MPrS≥DjTj (2)

where the operation period T is divided into M intervals, and each interval has duration Tj and a required demand level Dj . In the same context, as advancement in reliability applications, Taboada et al. [7] have generalized the use of LOLP as a reliability index. Using a series-parallel system, they calculate the availability of each part of the power including transmission and distribution system using LOLP model.

#### 2.1.1. Loss of load probability formulation

The generation system reliability calculations are based mainly on two analytical methods, which differ by the load model used. The first method is the calculation of the LOLP; there the load is given by the load duration curve. The second method is the frequency and duration approach, by which, besides the probability, the frequency and duration of load levels higher than the generation capacity can be determined. The loss of load probability method associates each value of MW outage with a certain cumulative probability, thus producing a capacity outage table. The expected load loss of the system is obtained from the capacity outage table, and a daily load peak variation curve is derived from the daily load curves. The LOLP, number of days on which capacity is insufficient, is obtained by adding the probability that the amount of capacity on forced outage, on day i is greater than or equal to the reserve on day i , for all days of the period being studied.

For a system state where the remaining generating capacity is Cj , the percentage of time tj during which the load demand exceeds Cj can be determined from the load curve L . The overall probability that the load demand will not be met is called the loss of load probability and is given by the Eq. [ 4 ]:

 LOLP=∑jPC=CjPL≻Cj=∑jpjtj100 (3)

where pj is the probability associated to the number of the failed generating units at time tj , and it is formulated through the following development.

In addition to complete failures, generating units may experience partial failure, when they continue to operate but at reduced capacity levels. They are also taken out of service from time to time for preventive maintenance. Using a simple two-state model for the operation of a unit, its failure probability is given depending on its failure outage rate (FOR), which can be assumed as the unit steady-state unavailability denoted A¯ . If at time tj , r units have failed from a total of n identical and independent installed units in the generating system, the probability pj is given by:

 pj=nrA¯r1−A¯n−r (4)

A case of unequal size of the units can appear.

### 2.1.2.1. First purpose

It is well known that availability is a measure of success used primarily for repairable systems. For nonrepairable systems, availability At equals reliability Rt . In repairable systems, At will be equal to or greater than Rt . In the optimistic case, the availability is greater than the reliability. Following the Levitin and Lisnianski development for a multistate generating system (MSGS), the availability expectation is the function of demand D and may be defined as [8]:

The index 1EA is often used and treated as loss of load probability and can be written as:

 LOLP=1−EAOrEA=1−LOLP (6)

### 2.1.2.2. Second purpose

This purpose highlights the correlation between the system reliability, the energy availability, and the loss of load probability. To understand this correlation, we consider a multistate generating repairable system MSGS connected to a load L, and on a given period of time, we draw two curves representing the evolutions of the system available capacity (SAC) and the hourly system load (HSL), respectively, as shown in Figure 1. Depending on the states of generating units (up or down) that involve partial or total failure of a simple unit or of several units, the appearance of dips in the same curve reflects units’ breakdowns, and the resumption to the initial level of capacity indicates that repairs were made. One of the most reliability indices that concerns more the utility than the customer is the energy not supplied (ENS) given by the dashed lines under the curve. Their corresponding time intervals denote durations, where the consumption exceeded the production, and therefore, we have loss of load. The decreasing level of system reliability is highlighted by degrading state, corresponding to each decreasing in the SAC curve behavior.

### Figure 1.

Superimposition of the system available capacity and the load model.

The generating system failures can occur in two ways: either through unit failures or through load increases. There is a loss of load when the demand is greater than the supply. However, there is a loss of supply when a failure occurs in the upstream of the load point. It is important to retain this difference. In this context, two relevant questions are to be asked and the answers were given in [9, 10] with a case study application, such as: What happens when load increases? What is the consequence of generating system failure?

### 2.2. Frequency and duration indices

Almost every electricity utility computes reliability indices on an annual basis. The most important reliability indices involving decision-making criteria are given as follows [1]:

The Expected Frequency of Load Curtailment in (fault/yr):

 EFLC=∑k=1nλk (7)

The Expected Duration of Load Curtailment in (hrs/yr):

 EDLC=∑k=1nλkTk (8)

The Expected Energy Not Supplied in (kWh/yr):

 EENS=L.EDLC (9)

where λk,Tk are failure rate and failure duration of an item k and L is the load curtailed at a considered load point, respectively. Application is done for a part of the distribution system of Algiers city (Algeria). Considering the electrical characteristics (network topology, section length, power value at load points and the fault search method) and reliability parameters mentioned earlier, the overall system reliability indices are computed.

### 2.3. Reliability indices improvement

To improve the reliability level, technical and organizational measures are considered during system planning and operation. The actions currently carried out are as follows: intensifying the operations of maintenance; networks reorganization, looping and meshing systems, and automation of networks. In [11], some options are added such as load transfer between feeders, undergrounding circuits, and replacement of aging equipment. From a practical standpoint, this application allows to highlight the goodness of each measure to the system performances by a simple comparison of reliability indices. The results of reliability indices improvement are published in [1].

## 3. Interruptions modeling and reliability service differentiation

In several studies dealing with electrical distribution system reliability, the objective often sought by the energy distributor is the balance between the required reliability level and its cost. In the following, we develop two important points of view in the reliability of the electrical systems: that relating to modeling and that of the differentiation of electricity prices according to the level of reliability required with a minimal guaranteed reliability level for customers without any prior requirement [12].

### 3.1. Interruption modeling using the Weibull-Markov process

In the last decade, a novel vision of interruption modeling in power systems was developed and consists of the Weibull-Markov process. The purpose is to model the failure and operating data according to Weibull distribution proprieties, while retaining those assigned to the Markov model where the system occupies discrete states. This process was initially developed by Van Castaren [2] and was applied successfully by Pivatolo [13] and Medjoudj et al. [3]. Applications were made to highlight maintenance policies gathered on three types of actions: namely nondestructive action which does not improve reliability level but slows the system degradation. This action, denoted 1a as a minor maintenance, is characterized by an improvement factor m1 . A second action is considered and denoted 2b and can touch some of the components of a system up to their replacement. To this action is associated an improvement factor m2 , and the maintenance is a major one. The third and final proposed maintenance action is on the renewal of equipment, and it is assumed to be perfect, and after its implementation, the system is assumed as good as new, and it is denoted 2p . From a practical standpoint, this action is highlighted by taking m1 and m2 equal to the unity. This notion is introduced by Tsai et al. [14] for a mechatronic system and applied for power systems by Medjoudj et al. [3]. In this part of the section, we introduce the concept of the differentiated reliability with an application to the case of an electrical MV/LV substation. Starting from the expression of reliability function expressed by a desired threshold, the need of performing preventive maintenance action at time is decided regarding the behavior of this function at the coming stage of maintenance. Then, the choice of the type of action to perform is dictated by the value of the maximum benefit brought by this action. Threshold reliability is allocated to the opposite risk of system failure occurrence.

#### 3.1.1. Reliability data analysis

Considering the formulations of mean up time and mean down time of an item is given respectively by:

 MUT=tm−tb∫0tmhtdt (10)

and

 MDT=ta−tb∫0tmhtdt, (11)

where, tm , ta and tb are respectively, the preventive maintenance PM interval, the PM and corrective maintenance CM times on replacement; the operational availability is defined as:

 A=tm−tb∫0tmhtdttm+ta (12)

Subsequently, the PM interval for maximizing the availability can be derived by differentiating Eq. 9 to time tm , such as dAdt=0 , and the differential result is:

 tm+tahtm−∫0tmhtdt=tatb (13)

For data treatment and statistical processing, forced and planned outages are collected over 17 years of system operation continuously at the national company of electricity and gas center (SONELGAZ) of Bejaia city, Algeria. For an MV/LV transformer which is a critical item of an electrical substation, the estimated parameters and the adequate probability distribution functions are listed in Table 1.

Component or Subsystem n DistributionParameters dks dn0.05 Decision
MT/LV Transformer17Weibull β=2.459579;
η=5.909754×106
0.2264 0.308 Not rejected
(Tr)Exponential λ=0.007785467 0.3136 Rejected

### Table 1.

Distribution functions parameters estimation.

The obtained results show that, based on the Kolmogorov Smirnov KS test [15], dks is lower than dn0.05 , the Weibull distribution is not rejected; however, with the exponential law, dks is greater than dn0.05 , the hypothesis is not accepted. In Table 2 are gathered reliability indices, where both the feeder failure frequency Fi and the transformer failure rate h are added. The nonacceptation of the exponential distribution is comforted by the review results of reference [16], where the authors state that the exponential law, usually used to describe failures, is not always 100% suitable for electricity distribution systems.

Component or Subsystem MUT (hours) MDT (hours) MTBF (hours) A Fi (1/year) h (1/year)
MT/LV Transformer 87358.33 9.97 87368.32 0.9998 0.10024 0.00778

### Table 2.

Reliability indices of the power transformer.

In this study, it is assumed that the substation failures are due either to the transformer or to the internal cable connector failures. A three states diagram (working, failure and maintenance) is dressed for the life cycle modeling of the substation as shown in Figure 2.

### Figure 2.

Three states diagram.

Let X12 , X13 , X1 , X2 , X3 be the random variables representing the duration of the operation until failure, the duration of the operation until maintenance, the duration of the operation (state S1 ), the duration of the interruption (state S2 ) and the duration of the maintenance (state S3 ), respectively. The estimated parameters of the random variables following the Weibull distributions are listed in Table 3.

Variable N DistributionParameters dks dn0.05 Decision
X12 27Weibull β=1.0644;
η=3.2827×106
0.1356 0.25438 Not rejected
X13 27Weibull β=1.6689;
η=7.2589×105
0.2173 0.25438 Not rejected
X2 27Weibull β=0.6764;
η=196.00159
0.1219 0.25438 Not rejected
X3 27Weibull β=1.03894;
η=123.6005
0.1219 0.25438 Not rejected

### Table 3.

Parameters estimation of the random variables following Weibull distributions.

#### 3.1.2. Reliability under preventive maintenance

The improvement of maintenance to reliability is developed using two factors, and the selection of the action to do for the components on every PM stage is decided by maximizing system benefit in maintenance. Depending on the percent of the survival parts of system when it is maintained, the reliability function is:

 Rjt=R0,j.RV,jt (14)

where R0,j is the initial reliability of the jth stage and RV,jt is the reliability degradation of surviving parts on this stage. Considering periodical PM which interval is tm , the reliability of surviving parts is defined as:

 RV,jt=R1m1(t−j−1tm (15)

With: j1tmtjtm and m1,0<m11 is the improvement factor of action 1a .

To model the reliability of systems following PM , the effects of various actions on R0,j and RV,j must be evaluated.

 R0,j=Rf,j−1=R0,j−1.Rtm (16)

where R0,j1 , Rf,j1 indicate the initial and final reliability values of the system on the j1th stage.

Action 1b can improve the surviving parts of the system and also recover the failed parts. Generally, the impact of this action on the failed parts can be measured by an improvement factor m2 , which is also set between 0 and 1 representing the restored level except the surviving parts. According to the definition, the initial reliability on the action 1b can be expressed as:

 R0,j=Rf,j−1+m2R0−Rf,j−1 (17)

where R0 denotes the initial reliability of the new system.

The system reliability is expressed as:

 Rjt=R0,jexpt−j−1t_m/m_1ηβ (18)

where R0,j is the initial reliability of the jth stage and m1 is the improvement factor of action 1a . The benefit of component maintenance on the jth stage is defined as [14]:

 Bi,k=∫tj∞Ri,j+1tdt−∫tj∞Ri,jtdtCi,k (19)

where i,k denote, respectively, the ith subsystem or component and the maintenance action considered and Ci,k , the action cost. The advantageous one will correspond to the maximum of the benefit, that is, Bi=MaxBi,k . Once the action of maintenance is defined and retained, the availability of the system at any stage is processed as:

 As,j=T−tb,m∑i=1n∫tj−1tjhi,jtdtT+∑inti,k,a (20)

where n is the number of components or subsystems and ti,k,a is the time of the PM actions 1a , 1b and 2p and T , the cycle time. In the following are described the different types of PM actions in the case of the power transformer.

• Action 1a : cleaning, lubricating, tightening and oil-level verification,

• Action 1b : oil and internal cable connectors’ I.C.C replacement,

• Action 2p : transformer replacement.

Parameters are needed to compute the benefits such as the distribution function parameters βη PM and CM times ( ta and tb ), maintenance actions costs ( C1a , C1b , and C2P ) listed in Table 4, and the threshold value of reliability ( Rth = 0.8).

β η tadays tbdays m1 m2 C1a $C1b$ C2p $2.45 4103.99 3.5 28 0.80 0.90 600 1500 8600 ### Table 4. The useful parameters for benefit evaluation. The obtained results are: ( tm I.C.C = 715 days, tm Tr = 2415 days); however, the maintenance interval for the system is Tm=min7152415=715 days. The maintenance action to retain is based on the maximum benefit value, and the results at different maintenance stages are listed in Table 5 using the following notations: • *: no maintenance is needed, • 0: nothing to do (after an inspection), • 1: action 1a is carried out, • 2: action 1b is carried out, • 3: action 2p is carried out, • Rj+1Tm : The instantaneous reliability at j+1×Tm . StageAction proposed Rj+1Tm Benefit$Action Retained
Rth=0.8 Rth=0.9 Rth=0.95 Rth=0.8 Rth=0.8 Rth=0.95 Rth=0.8 Rth=0.8 Rth=0.95
1 1a 0.92720.92720.9272 0.4378002
1b 0.8285
2p 0.1601
2 1a 0.81550.81550.9234 2.09440.6371012
1b 1.76180.9393
2p 0.32540.1827
3 1a 0.66180.82860.92303.02841.25260.6601122
1b 2.58591.80100.9394
2p 0.47100.33260.1844
4 1a 0.72870.92200.92221.7524 0.6830202
1b 2.1457 0.9393
2p 0.3945 0.1861
5 1a 0.92000.80560.9222 1.53050.7058012
1b 1.08510.9393
2p 0.34290.1877
6 1a 0.80380.80320.92193.14441.55990.7514112
1b 2.67120.34300.9392
2p 0.48810.31880.1910
7 1a 0.64550.87210.92151.72711.33590.7741212
1b 2.17350.60470.9392
2p 0.40160.36870.1927
8 1a 0.71810.83760.92111.72711.33590.7741212
1b 2.17350.60470.9392
2p 0.40160.36870.1927
9 1a 0.91900.79160.9207 1.33030.7967012
1b 0.92710.9392
2p 0.43000.1944
10 1a 0.80240.73850.9203 1.30840.8193012
1b 1.21770.9391
2p 0.48670.1960

### Table 5.

Maintenance plan depending on reliability thresholds.

The results listed in Table 5 can be interpreted as follows: at every maintenance stage, verify for each component if its reliability for the coming stage is greater or equal to the Rth .

• If the condition is realized, the decision is doing nothing. For the example of the transformer, at j=1 , the reliability is Rj+1×Tm=R2×Tm=0.9272>Rth=0.80 ; however, no maintenance is needed for the stage j=1 .

• If no, compute the benefit for each action proposed and choose the maximum value. For example, for the threshold value Rth=0.95 , at the first stage, j=1 , the reliability at the coming stage of maintenance Rj+1×Tm=R2×Tm=0.9272<Rcrit=0.95 and actions (1a , 1b) benefits are 0.43780.8285, respectively. The action 1b is retained looking at the maximum value of the benefit.

The risk management is highlighted by thresholds of reliability. Depending on the reliability level reached, or fixed a priori, maintenance operations can be decided. The objectives are the determination of maintenance frequencies on an item and consequently their costs. It will be remarked that a high level of reliability is required (i.e., the risk of failure is minimized), the maintenance frequency increases, and subsequently, the cost increases.

## 4. Competing failure processes of oil circuit breaker

The components constituting a high voltage oil circuit breaker (HVOCB) are subject to various degradations, namely the aging of the insulating oil in the arc’s extinguishing chamber, the contacts wear out and the sharp breakdown of bus bars supports. In this section, we have modeled the behavior of this item as subject to a competing three degradation processes by using the Markov state diagram as given in Figure 3. The states were defined using thresholds of degradation parameters. To the degradation processes was associated a shock process highlighting the effects of short circuit arrivals on the HVOCB when defaults occur at the downstream feeder. The novelty in this work is outlined by the use of three-dimensional matrix to show the possible states, where the HVOCB can sojourn.

#### Figure 3.

Diagram of transition states of a system subject to four failure processes.

### 4.1. Case of three degradation processes modeling

We consider that the processes of degradation are modeled using continuous probability functions, and the operating condition of the system is characterized by a number of states which space is noted by Ωμ .

Following the Li and Pham theory [17], we consider the three state spaces Ω1 , Ω2 ,and Ω3 corresponding to the degradation processes Y1t,Y2t,andY3t, respectively. After obtaining the state spaces Ω1 , Ω2 , and Ω3 , we develop a methodology to establish a relationship between the states of the system Ωμ , the set of degradation states and catastrophic state due to shocks arrivals Ω1Ω2Ω3F.

The study deals with the three processes of degradation Y1t,Y2t,andY3t combined with the shock process denoted Dt as given in Figure 3. The sets of states are represented by Ω1=M1,,11,01 , which corresponds to the degradation 1 with M1+1 states, Ω2=M2,,12,02 , corresponds to the degradation 2 with M2+1 , and Ω3=M3,,13,03 corresponds to the degradation 3 with M3+1 states. Figure 3 shows the transition between states of a system submitted to four failure processes.

The equivalence relations between degradation states Ω1=M1,,11,01 ; Ω2=M2,,12,02 and Ω3=M3,,13,03, and their corresponding intervals are given as follows:

• Degradation Process 1:

0<Y1tWM State M1.

WM<Y1tWM1 State M11.

.

.

.

W2<Y1tW1 State 11.

G1=W1<Y1t State 01.

• Degradation Process 2:

0<Y2tAM State M2.

AM<Y2tAM1 State M12.

.

.

.

A2<Y2tA1 State 12.

G2=A1<Y2t State 02.

• Degradation Process 3:

0<Y3tZM State M3.

ZM<Y3tZM1 State M13.

.

.

.

Z2<Y3tZ1 State 13.

G3=Z1<Y3t State 03.

The states’ space of the system is defined by Ωμ=M10F with M+2 states. In this part, we develop a function that generates the relation between the states’ space of the system Ωμ and degradation states’ space Ω1Ω2Ω3F . For example, at a given time t, it is assumed that the degradation process 1 is at the state i1Ω1 and the degradation process 2 is at state j2Ω2 and the process of degradation 3 is in the state k3Ω3 . It is assumed that at the present time, the system is not at fault condition (catastrophic state F ). Thus, state F can be ignored for the moment, we must therefore seek a function relationship between Ω and Ω1Ω2Ω3 . Instead of Ωμ and Ω1Ω2Ω3F . The operation can be described by the mathematical function formulated as follows:

 f:R=Ω1×Ω2×Ω3→Ω=M….10 (21)

where R=Ω1×Ω2×Ω3=i1j2k3i1Ω1j2Ω2k3Ω3

The function f is defined by: Ω1×Ω2×Ω3fH.

The matrix H represented in Figure 4 gives information about the resulting states space and component of M+1 elements corresponding to each space leaving by the function f . The line at the top of the matrix H represents the states of the degradation process 1, the right column of the matrix represents the states of the degradation process 2, and the top page of the matrix H represents the degradation process 3. The elements of the matrix H represent the states f(i1,j2, k3)=L.

### Figure 4.

Three-dimensional matrix of system’ states, where I1Ω1 , J2Ω2 , K3Ω3 , and LΩ .

We note that in the matrix H, some elements are zeros and it can be assumed that when degradation 1 is in a certain advanced state I101<I1<M1 , the degradation 2 is also in a certain weak state I202<I2<M2 , and the degradation 3 is also in a certain weak state I303<I3<M3 , it is considered as a failure condition. We also notice that fM1M2M3=M , and initially, the system is in a perfect state.

We define time until failure by: T=inf{t:Y1t>G1,Y2t>G2,Y3t>G3 where D>S} .

It is important to know that the life of the system depends on a single process among the three degradations and of that of the shock. However, the system failure is caused by the process that occurs first, exceeding its critical value corresponding to the level which can bring the system back to failure.

### 4.2. System case study modeling and application

Initially, the system is considered in good states of operation (M1,M2 , and M3) . It can pass first, to the degradation states M11M12M13 or to the state of catastrophic failure (state F ), due to random shock. When the system reaches the first state of degradation, it can either remain in this state or go to the second degradation state M21M22M23 , or it passes to state F . The same process is repeated at each degradation stage with the exception of the states 010203

Assumptions:

1. The system occupies M+2 states, where 0 and F are the states of failure, state i is a degradation state, 1<i<M ;

2. No repair or maintenance is carried out on the system;

3. Yit; i=1,2,3 ; is a not decreasing and not negative function. With respect to time t , it corresponds to an irreversible accumulation of damage;

4. Yit; i=1,2,3; and Dt is statistically independent implying that the state of one process will have no effect on the other state;

5. At time t=0 , the system is at state M ;

6. The system may fail due to;

• Degradation process if: Yit>Gi for i=1,2,3 ;

• Random shock process (the system passes to the condition of the catastrophic failure state F ), if: Dt=t=1NtXt>S ; with Gi and S are critical level of degradation and shocks, respectively.

The reliability function is defined as follows:

 RMt=Pstate≥1=∑i=1MPit=PY1≤G1Y2≤G2Y3≤G3Dt≤S (22)

The system will fail if any of the degradation rates exceeds the critical level Gi ; i=1,2,3 or the process of shock also exceeds the critical level S .

The system subject to three displacement processes is defined by:

1. The process of increasing degradation representing the wear of the contacts of the circuit breaker is denoted by Y1t ;

2. The process of increasing degradation representing the aging of oil insulating circuit of the circuit breaker is denoted by Y2t ;

3. The degradation process of bus bars supports is denoted by Y3t ;

4. A random process of cumulative shock damages is given by Dt=t=1NtXt .

We obtain a system with four competing degradation processes. For the wear of the contacts: Y1t , M1=2 ; for the aging of the oils Y2t : M2=2 ; and for the degradation of the supports: Y3t , M3=3 ; the system fails if the process of degradation Yit exceeds a level Gi ; i=1,2,3 or the process of damages cumulates shocks, and D (t) exceeds the level S. It is assumed that the state spaces associated to Y1 , Y2 and Y3 are Ω1=211101 , Ω2=221202 , Ω3=33231303, respectively. Consequently, the space of the system is defined as follows: Ω=3210F . Thus, the function f is defined as being the Cartesian product of three sets following f:R=Ω1×Ω2×Ω3Ω=321, and the result is illustrated by the three-dimensional matrix given in Figure 5.

### Figure 5.

The three-dimensional matrix of the studied system.

The implementation of the abovementioned models under the MATLAB software which has given the results of the probabilities of the sojourns in different states and the system reliability changing is shown in Figure 6. The system (oil circuit breaker, bus bars) is in good condition for 3 years (1100 days) with a probability greater than 0.9. After this period, the latter decreases exponentially to reach the zero value after 4 years (1490 days) without any maintenance actions.

### Figure 6.

Probabilities of the states and the system reliability changing.

The evolution of the probability of the system in degradation state 2 is complementary to that of the probability of the degradation state 3. Indeed, during 3 years of operation, the probability of being in state 2 is zero. Then, it increases exponentially to reach its maximum value of 0.41 up to 4 years (1280 days). As this state is transient, its probability function decreases to 0 after 5 years. For the system reliability, we note that during the first 8 years, it is expected to decrease by 20% due to the random shock process that governs the system during this period.

## 5. Conclusion and discussions

Recently researchers in electrical systems have proposed differentiated electricity service based on reliability and have shown some inconveniences to apply it into a real case. In the same location are connected both consumers with high reliability requirements, with an agreement to pay more and others who are not concerned. Because the technical measure proposed is to add a reliable feeder, the differentiation is not quite possible. Our proposition consists of organizational measure and is oriented to maintenance actions on MV/LV substations, which are directly connected to the end users of the network. The differentiated service is directly related to the reliability of the substation where the improvement is a function of maintenance actions and the frequency of interventions. For statistical considerations and for interruption (forced or scheduled outages) modeling, we have applied the Weibull-Markov approach rather than the Markov method, which is usually used for the case of electrical systems. It has been proven that it is possible to maintain in another way than the classical one based on systematic preventive maintenance. In this chapter, it is shown that the maintenance is decided on the reliability level and benefit bases. Another critical component of electrical substation is studied using competing failure processes and consists a circuit breaker. The reliability aspects are formulated in the bases of oil aging, contacts wear and bus bars support degradation. Investigations conducted by Pham in a theoretical framework have been applied successfully to complement system such as electrical system. The models applied on simple numerical examples have been validated by application to a real case engineering. During system operation, the results analysis of the network current state allows to the decision maker to reach better information and target the equipment that reduces the performances of the system and practicing suitable maintenance actions. Recent studies in energy sustainability and smart energy grid have revealed that reliability is the main criterion taken into account by decision makers in electricity market behavior and a performance index for electric utilities classifications.

## Acknowledgements and tributes

In 1988, a group of professors in mathematics and engineering (A. Aissani, D. Aissani, K.D. Haim, A. Boubakeur, and A. Ouabdeslam) had organized the national conference named MFSI at the University of Bejaia, in Algeria, where during 2 days, the notion of reliability was vulgarized in the field of engineering. From that manifestation was born the group of work in power system reliability FSE2. Over than 1500 various works were conducted in engineer, master and doctorate theses dealing with all the aspects of power systems reliability and with a large cooperation with other universities and various manufactures and services. A great number of applications were done around the power systems including production, transportation and distribution parts. In recent years, a lot of novelties were developed compared to what is done over the word, such as the Weibull-Markov modeling in data analysis, nonparametric distributions in switching components behavior, Box and Jenkins models in blackouts forecasting and reliability aspects in smart grids development and multicriteria optimization. The results were valorized in a great number of international conference proceedings and in valuable international journals. This chapter dealing with power system reliability constitutes an interesting opportunity to express our acknowledgments and tributes to these pioneers of reliability in Algeria for what they have given for research.

## References

1 - Medjoudj R, Aissani D, Haim KD. Power customer satisfaction and profitability analysis using multi-criteria decision making methods. International Journal of Electrical Power and Energy Systems (Elsevier). 2013;45(1):331-339. DOI: 10.1016/j.ijepes.2012.08.062
2 - Van Castaren J. Reliability assessment in electrical power systems: The Weibull-Markov stochastic model. IEEE Transactions on Industry Applications. 2000;36(6)
3 - Medjoudj R, Aissani D, Boubakeur A, Haim KD. Interruption modeling in electrical power distribution systems using Weibull-Markov model. Proceedings of the Institution of Mechanical Engineers (IMEchE), Part O: Journal of Risk and Reliability. 01-06-2009;223:145-157. DOI: 10.1243/1748006XJRR215
4 - Iberraken F, Medjoudj R, Medjoudj R, Aissani D. Combining reliability attributes to maintenance policies to improve high-voltage oil circuit breaker performances in the case of competing risks. Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability 1748006X15578572. First published on 02-04-2015. DOI: 10.1177/1748006X15578572
5 - Kuech JD, Kirby BJ, Overholt PN and Markel LC. Measurement practices for reliability and power quality. U.S Department of energy. ORNL/TM-2004/91
6 - Nourelfeth M, Ait-Kadi D. Optimization of series-parallel multi-state systems under maintenance policies. Reliability Engineering and System Safety. 2007;92:1620-1626
7 - Taboada HA, Espiritu JF, Coit DW. Design allocation of multistate series-parallel systems for power systems planning: A multiple objective evolutionary approach. Proceedings IMEchE, Part O: Journal of Risk and Reliability. 2008:381-391
8 - Levitin G, Lisnianski A. Multi-state system reliability analysis and optimization (UGF and GAA). In: Pham H, editor. Chapter 4, Handbook of Reliability Engineering. Springer; 2003
9 - Medjoudj R. Reliability Aspects of Electrical Systems: Semi-Markovian Modeling and Maintenance Optimization. Doctorate thesis: University of Bejaia, Algeria; 2009
10 - Iberraken F, Medjoudj R, Aissani D. Mathematical models to support the issue of electrical blackouts in the context of smart grids. Lecture Notes on Information Theory, Engineering and Technology Publishing. December 2014;2(4). DOI: 10.12720/lnit.2.4.295-301 (LNIT,ISSN: 2301-3788, www.lnit.org), 2014
11 - Endrenyi J. Reliability Modeling in Electric Power Systems. Toranto: Wiley & Interscience; 1978
12 - Brown ER, Marshall MW. The Cost of Reliability. ABB consulting report. 2001
13 - Pievatolo A. The downtime distribution after a failure of a system with multistate independent components. IMATI-CNR Technical report, Presented at the Mathematical Methods on Reliability conference, MMR, Glasgow, 2007
14 - Tsai YT, Wang KS, Tsai LC. A study of availability centred preventive maintenance for multi-component systems. Reliability Engineering and System Safety. 2004;84:261-269
15 - Dhilon BS. Design Reliability, Fundamentals and Applications. USA: CRC press LLC; 1999
16 - Billinton R, Goel R. An analytical approach to evaluate probability distributions associated with the reliability indices of electric distribution systems. IEEE Transactions on Power Delivery. 1986;1(3):245-251
17 - Li W, Pham H. Reliability modeling of multi-state degraded systems with multi-competing failures and random shocks. IEEE Transactions on Reliability. 2005;54(2):297-230