Open access peer-reviewed chapter

# Power System Reliability: Mathematical Models and Applications

By Rabah Medjoudj, Hassiba Bediaf and Djamil Aissani

Submitted: December 13th 2016Reviewed: October 25th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71926

## 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
• 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 Sand Das the supply and the load demand respectively, they compute the reliability of a multistate system (MSS) as:

R=PrSDOrR=1LOLPE1

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=1j=1MTjj=1MPrSDjTjE2

where the operation period T is divided into M intervals, and each interval has duration Tjand 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 iis 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 tjduring which the load demand exceeds Cjcan 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=CjPLCj=jpjtj100E3

where pjis 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, runits have failed from a total of nidentical and independent installed units in the generating system, the probability pjis given by:

pj=nrA¯r1A¯nrE4

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 Atequals reliability Rt. In repairable systems,Atwill be equal to or greater thanRt. 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 Dand may be defined as [8]:

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

LOLP=1EAOrEA=1LOLPE6

#### 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.

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λkE7

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

EDLC=k=1nλkTkE8

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

EENS=L.EDLCE9

whereλk,Tkare failure rate and failure duration of an item kand Lis 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 1aas a minor maintenance, is characterized by an improvement factor m1. A second action is considered and denoted 2band 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 takingm1and m2equal 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=tmtb0tmhtdtE10

and

MDT=tatb0tmhtdt,E11

where,tm, taand tbare respectively, the preventive maintenancePMinterval, the PMand corrective maintenance CMtimes on replacement; the operational availability is defined as:

A=tmtb0tmhtdttm+taE12

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

tm+tahtm0tmhtdt=tatbE13

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 SubsystemnDistributionParametersdksdn0.05Decision
MT/LV Transformer17Weibullβ=2.459579;
η=5.909754×106
0.22640.308Not rejected
(Tr)Exponentialλ=0.0077854670.3136Rejected

### Table 1.

Distribution functions parameters estimation.

The obtained results show that, based on the Kolmogorov Smirnov KStest [15], dksis lower thandn0.05, the Weibull distribution is not rejected; however, with the exponential law, dksis greater thandn0.05, the hypothesis is not accepted. In Table 2 are gathered reliability indices, where both the feeder failure frequency Fiand the transformer failure rate hare 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 SubsystemMUT(hours)MDT(hours)MTBF(hours)AFi(1/year)h(1/year)
MT/LV Transformer87358.339.9787368.320.99980.100240.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.

Let X12, X13, X1, X2, X3be the random variables representing the duration of the operation until failure, the duration of the operation until maintenance, the duration of the operation (stateS1), the duration of the interruption (stateS2) and the duration of the maintenance (stateS3), respectively. The estimated parameters of the random variables following the Weibull distributions are listed in Table 3.

VariableNDistributionParametersdksdn0.05Decision
X1227Weibullβ=1.0644;
η=3.2827×106
0.13560.25438Not rejected
X1327Weibullβ=1.6689;
η=7.2589×105
0.21730.25438Not rejected
X227Weibullβ=0.6764;
η=196.00159
0.12190.25438Not rejected
X327Weibullβ=1.03894;
η=123.6005
0.12190.25438Not 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 PMstage 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,jtE14

where R0,jis the initial reliability of the jthstage and RV,jtis the reliability degradation of surviving parts on this stage. Considering periodical PMwhich interval is tm, the reliability of surviving parts is defined as:

RV,jt=R1m1(tj1tmE15

With: j1tmtjtmand m1,0<m11is the improvement factor of action 1a.

To model the reliability of systems followingPM, the effects of various actions on R0,jand RV,jmust be evaluated.

R0,j=Rf,j1=R0,j1.RtmE16

where R0,j1, Rf,j1indicate the initial and final reliability values of the system on the j1thstage.

Action 1bcan 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 1bcan be expressed as:

R0,j=Rf,j1+m2R0Rf,j1E17

where R0denotes the initial reliability of the new system.

The system reliability is expressed as:

Rjt=R0,jexptj1t_m/m_1ηβE18

where R0,jis the initial reliability of the jthstage and m1is the improvement factor of action 1a. The benefit of component maintenance on the jthstage is defined as [14]:

Bi,k=tjRi,j+1tdttjRi,jtdtCi,kE19

where i,kdenote, respectively, the ithsubsystem 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=Ttb,mi=1ntj1tjhi,jtdtT+inti,k,aE20

where nis the number of components or subsystems and ti,k,ais the time of the PMactions 1a, 1band 2pand T, the cycle time. In the following are described the different types of PMactions 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.Creplacement,

• Action 2p: transformer replacement.

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

βηtadaystbdaysm1m2C1a$C1b$C2p$2.454103.993.5280.800.9060015008600 ### Table 4. The useful parameters for benefit evaluation. The obtained results are: (tmI.C.C= 715 days, tmTr= 2415 days); however, the maintenance interval for the system is Tm=min7152415=715days. 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: 1. *: no maintenance is needed, 2. 0: nothing to do (after an inspection), 3. 1: action 1ais carried out, 4. 2: action 1bis carried out, 5. 3: action 2pis carried out, 6. Rj+1Tm: The instantaneous reliability at j+1×Tm. StageAction proposedRj+1TmBenefit$Action Retained
Rth=0.8Rth=0.9Rth=0.95Rth=0.8Rth=0.8Rth=0.95Rth=0.8Rth=0.8Rth=0.95
11a0.92720.92720.92720.4378002
1b0.8285
2p0.1601
21a0.81550.81550.92342.09440.6371012
1b1.76180.9393
2p0.32540.1827
31a0.66180.82860.92303.02841.25260.6601122
1b2.58591.80100.9394
2p0.47100.33260.1844
41a0.72870.92200.92221.75240.6830202
1b2.14570.9393
2p0.39450.1861
51a0.92000.80560.92221.53050.7058012
1b1.08510.9393
2p0.34290.1877
61a0.80380.80320.92193.14441.55990.7514112
1b2.67120.34300.9392
2p0.48810.31880.1910
71a0.64550.87210.92151.72711.33590.7741212
1b2.17350.60470.9392
2p0.40160.36870.1927
81a0.71810.83760.92111.72711.33590.7741212
1b2.17350.60470.9392
2p0.40160.36870.1927
91a0.91900.79160.92071.33030.7967012
1b0.92710.9392
2p0.43000.1944
101a0.80240.73850.92031.30840.8193012
1b1.21770.9391
2p0.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.95and actions (1a, 1b)benefits are 0.43780.8285,respectively. The action 1bis 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.

### 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Ω3corresponding 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,andY3tcombined with the shock process denoted Dtas given in Figure 3. The sets of states are represented by Ω1=M1,,11,01, which corresponds to the degradation 1 withM1+1states, Ω2=M2,,12,02, corresponds to the degradation 2 with M2+1, and Ω3=M3,,13,03corresponds to the degradation 3 with M3+1states. 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,02and Ω3=M3,,13,03,and their corresponding intervals are given as follows:

0<Y1tWMState M1.

WM<Y1tWM1State M11.

.

.

.

W2<Y1tW1State 11.

G1=W1<Y1tState 01.

0<Y2tAMState M2.

AM<Y2tAM1State M12.

.

.

.

A2<Y2tA1State 12.

G2=A1<Y2tState 02.

0<Y3tZMState M3.

ZM<Y3tZM1State M13.

.

.

.

Z2<Y3tZ1State 13.

G3=Z1<Y3tState 03.

The states’ space of the system is defined by Ωμ=M10Fwith M+2states. 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Ω1and the degradation process 2 is at state j2Ω2and the process of degradation 3 is in the statek3Ω3. It is assumed that at the present time, the system is not at fault condition (catastrophic state F). Thus, state Fcan 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.10E21

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

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

The matrix Hrepresented in Figure 4 gives information about the resulting states space and component of M+1elements corresponding to each space leaving by the functionf. The line at the top of the matrix Hrepresents 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 Hrepresents the degradation process 3. The elements of the matrix Hrepresent the states f(i1,j2,k3)=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>G3where 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 M11M12M13or 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 stateF. The same process is repeated at each degradation stage with the exception of the states 010203

Assumptions:

1. The system occupies M+2states, where 0and Fare the states of failure, state iis 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 Dtis 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>Gifor 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 Giand Sare critical level of degradation and shocks, respectively.

The reliability function is defined as follows:

RMt=Pstate1=i=1MPit=PY1G1Y2G2Y3G3DtSE22

The system will fail if any of the degradation rates exceeds the critical level Gi; i=1,2,3or 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 Yitexceeds a level Gi; i=1,2,3or the process of damages cumulates shocks, and D (t) exceeds the level S. It is assumed that the state spaces associated toY1,Y2and Y3are Ω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.

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.

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.

## How to cite and reference

### Cite this chapter Copy to clipboard

Rabah Medjoudj, Hassiba Bediaf and Djamil Aissani (December 20th 2017). Power System Reliability: Mathematical Models and Applications, System Reliability, Constantin Volosencu, IntechOpen, DOI: 10.5772/intechopen.71926. Available from:

### Related Content

Next chapter

#### Techno-Economic Feasibility Study of Autonomous Hybrid AC/DC Microgrid System

By Atanda K. Raji

First chapter

#### Microassembly Using Water Drop

By Taksehi Mizuno

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books