Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

reading / British practice / Vol D - 1990 (ocr) ELECTRICAL SYSTEM & EQUIPMENT

.pdf
Скачиваний:
9
Добавлен:
17.07.2023
Размер:
30.18 Mб
Скачать

Reliability evaluation of power systems

TABLE 2.2

Second order — maintenance outage overlapping forced outage (TYPE 1 and TYPE 2)

 

 

Event

Failure rate

Residence time

 

 

failure

X

 

 

 

sequence

r

--

 

 

 

TYPE I

 

 

 

 

 

 

(a)

I",

2

XP, 2r;

r 1 r2/(r; + r2)

 

(b)

2,

I

X5X 1 i.

riri/(ri +ri)

 

 

 

 

 

-

 

 

 

TYPE 2

 

 

 

 

 

 

(a)

1 0 ,

2

X;X2ri

r; t c /(r; + t c )

 

(b)

2,

I

XX1r5

r5, t c /(rir + tc)

 

 

 

 

 

 

 

Where: x n

= Failure rate of component 'n'

 

 

 

 

r n

Repair time of component 'n'

 

 

 

 

t c

= Switching time of alternative (NO) path

 

 

 

Components out on maintenance denoted thus (')

 

 

 

 

 

 

Indices for TYPE 1 and TYPE 2 events are:

U; =

 

a

a

=

= Ul/X5

a

a

 

 

For larger systems, due to the greatly increased number of components and branches, to do so manually would be very tedious and time consuming. However, part of the validation process for each version of the program, involves checking by inspection that the minimal path and cut set deduction is functioning correctly for at least one fairly large system.

Verification of the numerical algorithms is carried out during the program development by undertaking manual calculations, using the appropriate equations from Section 2.5.11 of this chapter, of busbar and system indices for a small system, and comparing the results with those obtained from using the computer program on the same system with the same component reliability data.

2.7 Typical applications

2.7.1 Example of the calculation and use of busbar indices

This section contains a specific example of the busbar indices evaluation techniques as they were used to answer a specific question relating the design of SESs for future nuclear power stations.

The higher voltage levels (down to 3.3 kV) of the SES shown in Fig 2.31 was the initial design proposal. It was put forward as part of an overall exercise to evolve a standard design of SES and grid connection interface arrangement, which would possibly be suitable for future nuclear power stations.

Depending on the exact siting of the power station, there can be a considerable variation in the reliability of the available grid connections.

The frequency of LOSP (loss of off-site power) contributes to the overall summated frequency of degraded core and hence forms part of the safety case.

The question raised was 'to what extent could the provision of additional on-site generation compensate for poor grid reliability?'

The procedure adopted was to assume a figure of 0.075 failures/year for the frequency of loss of all grid supplies for an average duration of 2 hours and to calculate the busbar indices for the II kV and 3.3 kV levels:

(a)For the basic system of Fig 2.31, which only has standby generators at the 3.3 kV level, designed to back up the SES.

(b)For the basic system, with additional standby generators on the 11 kV station boards.

113

Electrical system analysis

Chapter 2

 

 

TABLE 2.3

Third order overlapping forced outages (TYPE 1)

 

 

Event

Failure rate

Residence time

 

 

failure

 

 

 

 

sequence

X

r

 

 

 

 

 

 

 

 

 

 

 

TYPE 1

(a)

1,

2,

3

XIX2r1X3rir2/(ri + r2)

1ir2r3/(r0"2 + r2r3 + 1 3rI)

(b)

1,

3,

2

XiX3rIX2rir3/(ri + r3)

rir2r3/(rtr2 + r2r3 I- r3r1)

(c)

2,

3,

1

X2X3r2X1r2r3/(r2 + r3)

r 1 r2 r3 /(r 1 r2 -+ r2r3 + r3 r)

(d)

3,

2,

1

X3X2r3X1r3r2/(r3 + r2)

rir2r3/(rtr2 + r2r3 + r3r1)

(e)

2,

1,

3

X2X1r2X3r2ri/(r2 + IT)

rir2r3/(rir2 + r2r3 + r3r1)

(f)

3,

1,

2

X3k1r3X2r3r1/(r3 + ri)

rir2r3/(rt r2 + r2r3 + r3ri)

(g)

1 (2 & 3)

XIX23r1

r1r2r3/(r1r2 + r2r3 + r3r1)

(h)

(2 & 3), 1

X23X1r2r3/(r2 -I- r3)

rir2r3/(rir2 + r2r3 + r31- 1)

(i)

(I & 2 & 3)

X123

r tr2r3 /(r 1 r2 + r2r3 + r3r1)

 

 

 

 

 

 

Where: X n

X23

X123 rn

=Independent failure rate of component 'n' CMF rate for components 2 and 3

CMF rate for components 1, 2 and 3

=Repair time of component 'n'

Indices for TYPE 1 events are:

Xi = EX,

U1 = EXr, ri = U1/X1

a

a

 

 

(c)For the basic system, with additional standby generators on the 11 kV station and unit boards.

The following assumptions were made:

• The additional on-site generators provided at the 11 kV level would be standby generators and not run continuously..

As normal practice in respect of the safety case, no claim is made for on-site generation from the main units.

Description of computer produced system diagrams

In all diagrams, the following busbar identification applies:

S1

A dummy busbar to represent the

 

common failure rate/restoration

 

times for all the grid lines (irrespec-

 

tive number)

52, 53, S4

Individual sections of 400 kV busbar

and S5

 

S6 and S13 — Generator main connections

S7 and 514 — Unit 11 kV switchboards

S8 and S15 — Unit 3.3 kV switchboards

S10 and S12 — Station 3.3 kV switchboards.

Note that the actual method of producing the systems and data for each of the cases studied was to draw the complete system only for case (c) and complete the reliability data entry. Case studies (a) and (b) were then undertaken as sensitivity studies, using the 'outage' facility of the GRASP2 computer program. The plant outaged for each particular case study is shown by dotted lines on the appropriate results diagrams.

Results

Only two of the busbar indices are of interest to permit a comparison to be made between the system with no additional on-site generation (case (a)) and the other two systems (cases (b) and (c)):

The average failure rate (FR) is of interest, since every loss of supply to any of the four 11 kV busbars leads to a reactor trip.

114

 

 

 

 

 

 

 

 

Reliability evaluation of power systems

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TABLE 2.4

 

 

 

 

 

 

 

Third order overlapping forced outages (TYPE 2)

 

 

 

 

 

 

 

 

 

 

 

 

Event

 

Failure rate

 

Residence time

 

 

 

failure

 

 

 

 

 

 

 

sequence

 

X

 

r

 

-

 

 

 

 

 

 

TYPE 2

 

 

 

 

 

 

 

 

 

(a)

1,

2,

3

 

X iX2riX31- 11- 2 /(ri

+ 1- 2)

r1r2t c /(r 1 r2 I- r2tc + tcrl)

 

(b)

1,

3,

2

 

X 1 X3 1. 1 X2r1 r3/(r1

+ 1- 3)

ri te r3/(r 1 Lc + 10. 3 + r3ri )

 

(C)

2,

3,

I

 

X2X3r2X1r2r3/(r2 + r3)

t c r2r3/(t c r2 + r2r3 4- r3tc)

 

(d)

3,

2,

1

 

X3X2r3Xtr3r2/(r3 + r2)

ter2r3/(ter2 + r2r3 + 7'30

 

(e)

2,

1,

3

 

X2X1r2X3r2r1/(r2 + 1'1)

rir2t c /(rtr2 + r2te + teri)

 

(e)

3,

I,

2

 

X3X11- 3X2r3r1/(r3 + ri)

rit c r3/(r1 t c + to- 3 + r3ti)

 

(g)

1 (2 & 3)

 

X1X23r1

 

rticArt + lc)

 

(h)

(2 & 3), 1

 

X 23X i r2r3/(r2 + 1- 3)

t c r2r3/(t e r2 + r2r3 + r3t c ")

 

(i)

(1 & 2 & 3)

 

X123

 

t c

 

 

 

 

 

 

 

 

 

 

Where: X n

X23

X123

r n

tc

=Independent failure rate of component 'n'

=CMF rate for components 2 and 3

CMF rate for components I, 2 and 3 Repair time of component 'n'

= Switching time of alternative (NO) path

Indices for TYPE 2 events are:

X2 = E X, U2 =

Xr, r2 = U2/12

a

a

 

 

Average outage time (AOT) is of interest, since it gives an indication of the availability levels required of the SES to cater for the reactor post-trip (initiated by the electrical system) period.

The average failure rate for those events initiated by failures within the electrical system which cannot be recovered within 2 hours is also of interest. For the safety case, this indicates the improvement in the immediate reactor post-trip period that would be gained from providing the additional on-site generation.

The diagrammatic results for AOT and FR for case

(a) are shown in Figs 2.32 and 2.33 respectively. Those for cases (b) and (c) are shown in Figs 2.34 to 2.37.

Table 2.16 summarises both busbar indices and the (0-2 h) FR for each busbar for all three cases.

The (0-2 h) FR figures were obtained by manually sorting the list of failure events and combining only those events which had AOTs greater than 2 hours.

It can be seen from Table 2.16 that the effect of adding auxiliary standby generators at the 11 kV level would be to increase the average failure rates of the

11 kV and 3.3 kV boards, whilst only showing a marginal reduction in the corresponding average outage durations. A similar increase is shown in the (0-2 h) average failure rate.

Both these effects are predictable, since they are due mainly to the addition of the active failure rates of the 11 kV (NO) circuit-breaker associated with the additional standby generators, and these events are recoverable within the (0-2 h) period by switching/ isolation and the use of alternative sources of supply.

2.7.2 Example of the calculation and use of system indices

This section contains a specific example of the evaluation techniques for system indices as they were used to compare the performance of two alternative proposals for the design of a SES for the Coal-Fired Reference Design in terms of the availability of the main boiler/ turbine-generator unit.

It was also required to assess the cost of lost generation caused by failures within the SES over the

115

Electrical system analysis

Chapter 2

 

 

TABLE 2.5

Third order — maintenance outage overlapping forced outage (TYPE I)

 

 

Event

 

Failure rate

 

 

Residence time

 

 

failure

 

 

 

 

 

 

 

 

sequence

 

X

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TYPE I

 

(a)

1 ,

2,

3

 

XP, 2ri'X3rrr2/(q + r2)

 

gr2r3/(r; r2 + r2 r3 + r3rj')

 

(b)

1 0 ,

3,

2

 

XI'X3r1A2r 1'r3/(q + r3)

 

r 1 r2r3/(r1 r2 + r2r3 4- r3r;)

(0

2",

3,

1

 

X2 362,1*- 3/(r; + r3)

 

ririr3/(rir + rp. 3 + r3r/)

 

(d)

3",

2,

I

 

xp,21-ix it-P- 2/(ri ÷ r2)

 

r1r2r3 /(rir2 + r2q + 6'1'1)

(0

2",

1,

3

 

XpiriX3r?1/(ri + ri)

 

riri73/(rirl + r2 r3 + r3rj)

 

(f)

3",

1,

2

 

X3X1r;X2r3r1/(r3 + ri)

 

r1r2ri./(rir2 + r3r; + *I)

 

(g)

1" (2 & 3)

XcX23r1

 

 

rjr2r3/(rjr2 + r2r3 + r3rj)

 

 

 

 

 

 

 

 

 

Where:

X n = Independent failure rate of component 'n'

 

 

 

 

 

X23 = CMF rate for components 2 and 3

 

 

 

 

 

In

Repair time of component 'n'

 

 

 

 

 

 

 

Components out on maintenance denoted thus ("3

 

 

 

 

 

 

 

 

 

Indices for TYPE I events are:

 

 

 

 

 

 

X; = EX",

= EX'r",

= 14/Xr

 

 

 

 

 

 

a

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

designed life of the station for use in the cost benefit exercise to assess whether the capital cost of providing the interconnection could be justified.

The basic design of SES proposed, is shown in Fig 2.38 with two variants:

A system that is purely radial top-to-bottom, with

no switchboard interconnections at any voltage level.

A system as shown in Fig 2.38, but having normally open interconnectors between corresponding switchboards of each train at each level within the unit system and also at the 11 kV level between each unit and station switchboard.

System conditions

Since the exercise was of a comparative nature, all cables (except those used to interconnect boards in the one system) were excluded from the studies to simplify the analysis.

It was assumed for the purpose of the study that the station was in the CMR condition. Thus, for example, all coal mills and fans were modelled as being in operation, whereas only one sootblower and two oil pumps were included, the remaining sootblower and oil pump being on standby.

One complete unit system (Unit 1 Trains A and B) and half of the station system (Station Train 1) were represented in the computer model.

Figures 2.39 and 2.40 show the computer-produced diagrams for the radial system and Figs 2.41 to 2.43 show those for the interconnected system.

Table 2.17 shows the correlation between the com- puter-generated busbar numbering system and the individual system board names.

A grid failure rate of 0.039 failures/year was used for the 400 kV and 132 kV systems and failure rate data accumulated for generators of 500 MW and above were used for the 900 MW unit.

It was recognised that plant maintenance policy would play an important role in this assessment and a logical maintenance scheme was assumed, which included a certain amount of in-line maintenance.

To calculate the system indices, a derating factor was applied to each of the main unit boards (i.e., those at 3.3 kV and above) of each system. For example, it was assumed that loss of supply to 11 kV Unit Board 1A1 would result in a reduction in output power from the main generator of 50%. In addition, combinations of various unit boards were considered; for example, loss of both 3.3 kV mill boards was assumed to result in 100 0/o loss of unit output, whereas loss of 3.3 kV mill

116

 

 

 

 

 

 

 

 

 

Power system performance analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TABLE 2.6

 

 

 

 

 

Third order maintenance outage overlapping forced outage (TYPE 2)

 

 

 

 

 

 

 

 

 

 

 

 

Event

 

Failure rate

 

 

Residence time

 

 

failure

 

 

 

 

 

 

 

 

 

sequence

 

X

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TYPE 2

 

(a)

1",

2,

3

 

Xi'X2qX3ri'r2/(r1 + r2)

 

ilr2t c /(11'r2 + r21 c + tcr;')

 

(b)

I",

3,

2

 

>4X3r1X2ri'r3/(r1 + r3)

 

ri'r3t c /(11'r3 + qtc + r3lc)

 

(c)

2",

3,

I

 

X 2 X3riXtr5r3/(ri + r3)

 

r r3t c /(rD- 3 + *c + r3tc)

 

(d)

3",

2,

I

 

Xp, 2 rp+ 1 r3 r2 /(ri` + r2)

 

r3 r2t c /(r3 r2 + ri Lc + r2 lc)

 

(e)

2",

3,

1

 

X 2 X1r2 ?■ 3rIr1/(r

+ ri)

 

ri'rit c /(r5r1 + rp. c + rile)

(0

3",

1,

2

 

X3X1r3X2rPii(ri.

+r1)

 

ririt c /(riri + rit c + ritc)

 

(g)

1" (2 & 3)

 

Xj'X2311

 

 

rcte/(rj' + t c )

 

 

 

 

 

 

 

 

 

 

 

Where:

 

X n = Independent failure rate of component 'n'

 

 

 

 

 

X23 = CMF rate for components 2 and 3

 

 

 

 

 

 

 

 

r n = Repair time of component 'n'

 

 

 

 

 

 

 

 

t c = Switching time of alternative (NO) path

 

 

 

 

Components out on maintenance denoted thus (")

 

 

 

 

 

 

 

 

 

 

Indices for TYPE 2 events are:

 

 

 

 

 

 

 

 

= EX", U = EX"r", r = 1.1;Al

 

 

 

 

 

 

 

a

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

board IA and 11 kV fan board 1B would result in a loss of 50; output.

Loss of output due to maintenance of the main generator was not included in the system indices calculation, since this is usually a planned outage and is common to any auxiliary system. For the purpose of this analysis, loss of supply to unit boards at the 415 V level was not considered to cause loss (total or partial) of output, due to the difficulty of handling the large number of such boards and combinations of them with the GRASP computer program. However, since the analysis was comparative and the omission was made in both systems, the effects of the omission are negligible.

Results

The system indices calculated for the radial system ER and LOG are shown in Fig 2.44 and those for the interconnected system are shown in Fig 2.45.

The indices show that, although the frequency of encounter ER of the derated states is fractionally higher for the interconnected SES than for the radial system (0.28 encounters/year at 100% and 0.37 at 50% against 0.27 and 0.33, respectively), the loss of generation (LOG) is lower by approximately 25 GWh/year. With

replacement generation costs assessed at £25 per MWh, there would be a saving of £625 000 per unit per year or £37.5 million for a two-unit station over a 30 year life if the interconnected design of SES is used.

3 Power system performance analysis

3.1 Load flow analysis

3.1.1 Introduction

Load flow analysis is used to assist in planning, designing and operating power networks. The network or system is modelled in its steady state operating modes, then the voltage at each busbar or board and the power flow in each branch or circuit is calculated.

Particular attention is being paid here to the analysis of power station electrical systems. If a SES fails to provide sufficient power to motor drives or other equipment, then generating output may be lost completely or be significantly reduced. There are economic penalties associated with the loss of generation and there may be a hazard to plant safety. Thus, it is essential to have

117

Electrical system analysis

Chapter 2

 

 

TABLE 2.7

Active failure overlapping - with forced outages — second order (TYPE 3)

 

 

Event

 

Failure rate

 

Residence time

 

 

failure

 

 

 

 

 

 

 

Sequence

 

X

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TYPE 3

 

(a)

la, 2

 

)4X2s1

 

sir2/(s] + r2)

 

(b)

2, la

 

X 2X l r2

 

sI r2/(s1 + r2)

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

Where: X n = Independent failure rate of component

 

 

X na

= Active failure rate of component 'n'

 

 

rn

Repair time of component 'n'

 

 

 

 

sn = Switching/isolating time of component 'n'

 

 

 

 

 

 

 

Indices for TYPE 3 events are:

 

 

 

 

X3 = EX, U3 = EXr, r3 = U3/X3

 

 

 

 

a

a

 

 

 

 

 

 

 

 

 

 

TABLE 2.8

Maintenance outage overlapping active failure — second order (TYPE 3)

 

 

Event

 

 

Failure rate

 

Residence time

 

 

failure

 

 

 

 

 

 

 

 

sequence

 

 

X

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TYPE 3

 

(a)

2, 1 a

 

X 2" Xj r2"

 

 

sirl/(si + ri)

 

(b)

l", 2a

 

Ai'X2alii

 

 

r;s2/(r; + s2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Where: X n

Independent failure rate of component

 

 

X

 

failure rate of component 'n'

 

 

r n = Repair time of component 'n'

 

 

Sn = Switching/isolating time of component 'n'

 

 

Component on maintenance indicated thus (")

 

 

 

 

 

 

 

 

Indices for TYPE 3 events are:

 

 

 

 

 

= EX",

U

EX"r", r

= 115/Xi

 

 

a

 

a

 

 

 

 

 

 

 

 

 

 

 

 

118

 

 

 

 

 

 

 

 

Power system performance analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TABLE 2.9

 

 

 

 

 

 

Active failure overlapping forced outages third order (TYPE 3)

 

 

 

 

 

 

 

 

 

 

 

 

 

Event

 

Failure rate

 

Residence time

 

 

failure

 

 

 

 

 

 

 

 

sequence

X.

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TYPE 3

 

(a)

la,

2,

3

 

XjX2s1X3s1r2/(sl+r2)

 

sir2r3/(sir2 + r2r3 +r3st)

 

(b)

la,

3,2

 

XIX3s1X2r3/(si + 1- 3)

 

str2r3/(s1r2 +r2r3 +r3s1)

 

(c)

2,

3,

la

 

X2X3r2Xjr2r3/(r2 +r3)

 

sIr2r3/(str2 +r2r3 +r3s1)

 

(d)

3,

2,

la

 

X3X2r3Xjr3r2/(r3 +r2)

 

sir2r3/(str2 +r2r3 +r3s1)

 

(e)

2,

la, 3

 

X2Xjr2X3r2s1/(r2 +s1)

 

s1r2r3/(sir2 +r2r3+r3s1)

 

(f)

3,

la,

2

 

X3Xjr3X2r3s1/(r3+st)

 

51r2r3/(sir2+r2r3 +r3s1)

 

(g)

la (2 & 3)

XjX23s1

 

si r2r3/(si r2 + r2r3 + r3s1)

 

(h)

(2

& 3),

la

X23)4r2r3/(r2 +r3)

 

sir2r3/(sir2 +r2r3 +r3si)

 

 

 

 

 

 

 

 

 

 

Where:

X Tr = Independent failure rate of component 'n'

 

 

 

 

 

Xt = Active failure rate of component I

 

 

 

 

 

 

 

 

X23 = CMF rate for components 2 and 3

 

 

 

 

 

 

 

 

r n = Repair time of component 'n'

 

 

 

 

 

 

 

 

SI

= Switching/isolating time of component I

 

 

 

 

 

 

 

 

 

 

 

TABLE 2.10

Maintenance outage overlapping with active failure and forced outage — third order (TYPE 3)

 

Event

 

Failure rate

 

Residence time

 

failure

 

 

 

 

 

 

sequence

 

X

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TYPE 3

(a)

2", 3,

la

 

XiX3r2Xtrir3/(r1 +1- 3)

 

s1rr3/(siri +-sir3 +r?3)

(h)

2",

la,

3

 

XIXMX3r1s1/(r'i+st)

 

ft

 

 

sir 2r3/(sir -!- sir; +r?3)

(c)

3",

2,

I

 

X;X2eiXteir2/(ri + r2)

 

sIrS 1- 2/(sir"3 +slr2 +r"3r2)

(d)

3",

la, 2

 

X"3 Xf riX2rSsi/(1 +sI)

 

sir ir2/(siri +str2+r;r2)

 

 

 

 

 

 

 

 

Where:

X r, =

Independent failure rate of component 'n'

 

 

 

Xj

Active failure rate of component 1

 

 

 

 

 

 

= Repair time of component `ri'

 

 

 

 

 

si = Switching/isolating time of component I

 

 

 

Component on maintenance indicated thus (")

 

 

 

 

 

 

 

 

 

119

Electrical system analysis

 

 

 

 

 

 

 

 

 

Chapter 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

TABLE 2.11

 

 

 

 

 

 

 

Overlapping forced outages second order (TYPE .5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Residence time

(TYPE 5)

 

 

Event

 

Failure

 

 

 

 

 

 

 

 

 

 

 

Repair

 

Restoration by

 

Restoration

 

 

 

 

rate

 

rate

 

component

 

from LES

 

 

 

 

 

 

 

 

 

repair

 

 

 

 

 

 

 

 

A

 

 

ii

 

rA

 

rg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(a)

1, 2

 

XIX2r:

 

Al +n2

 

1 4t

 

rit c /(ri +tc)

 

(b)

2, 1

 

X2X1r2

 

1.4 I

+ n2

 

I/pi

 

r2tc/(r2-1-tc)

 

(c)

(I +2)

 

Al2

 

Al +n2

 

1 4t

 

t c

 

 

 

 

 

 

 

 

 

 

 

 

For each event: Average repair time r = rAP + rg (I —P)

 

 

 

 

 

Where:

P =

 

 

 

 

 

 

 

 

 

 

 

 

 

A n = Failure rate of component

 

 

 

 

 

 

 

r n

Repair time of component 'n'

 

 

 

 

 

 

 

t c = Switching and start-up time of LES

 

 

 

 

 

 

 

EL

Time limit of LES

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Indices for TYPE 5 events are:

 

 

 

 

 

 

 

 

 

X = Xa + Xb + A c

 

 

 

 

 

 

 

 

 

U

 

ra

Ab rb +

Ac rc

 

 

 

r = U/X

well designed power station electrical systems, so their operating voltage and frequency limits are defined; if the auxiliary plant is specified to be supplied within these defined voltage and frequency limits, then its performance will be acceptable.

3.1.2 Program construction [4,51

General remarks

The equations defining load flow in a power network can be stated and digital computation then used to solve them. A great deal of research has gone into finding the best methods of doing this. Special numerical techniques are used, the most popular being based on nodal methods of network analysis, using Newton-Raphson based algorithms to solve the network equations. Gauss-Seidel based algorithms and Z-matrix methods have also been developed. These methods are discussed in detail in [6].

A set of complex simultaneous equations is set up relating currents and voltages within the network. Although the basic equations are linear, non-linearities arise for various reasons. Hence the need for an iterative numerical solution.

There are many ways of describing an electrical network or power system mathematically: all must satisfy Kirchoff's voltage and current laws. Nodal analysis has emerged as the method best suited to digital computation.

Assuming a balanced three-phase network it is sufficient to use the positive phase sequence impedances (and hence admittances) of circuit components such as transmission lines and cables, transformers, series and shunt capacitors. These are usually represented in perunit (p.u.) notation as lumped linear complex impedances at rated system frequency. Transmission circuits may be represented by their equivalent r-networks.

In nodal analysis an electrical system comprising transformers, lines, cables, inductors, series capacitors, shunt capacitors and loads is translated into a set of 'branches' expressed as impedances or admittances. Admittance values are normally used, because there are considerable advantages during computation in using admittance instead of impedance values for network branches. The busbars and switchboards are separately translated into a set of 'nodes' which serve as link points for the branches. All voltage sources (generators) are replaced by current injection at the appropriate node.

120

V2 — V1 =

Power system performance analysis

 

 

 

 

 

 

 

TABLE

2.12

 

 

 

 

Maintenance outage overlapping forced outages — second order (TYPE j)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Residence time

(TYPE 5)

 

 

 

 

 

 

 

 

 

 

 

Event

 

Failure

Repair

 

Restoration by

 

Restoration

 

 

 

 

 

rate

rate

 

component

 

from LES

 

 

 

 

 

 

 

 

 

repair

 

 

 

 

 

 

 

 

 

 

 

 

 

r'E'3

 

 

 

 

 

 

 

 

 

 

 

 

 

(a)

I",

2

 

X "X2

"

A i' + n2

 

1//.4"

 

ri't c /(r; + t e )

 

 

 

 

 

l

rl

 

 

 

 

 

 

(b)

2",

1

 

X5X1r5

11 1 + n5

 

1/12"

 

1- t c /(r5 + tc)

 

 

 

 

 

 

 

 

 

For each event: average repair time r" = r A P" + r" (1 - P")

 

 

 

 

 

Where: P"

c _ ti.

 

 

 

 

 

 

 

 

X n

Failure rate of component 'n'

 

 

 

 

 

 

r n =

Repair time of component 'n'

 

 

 

 

 

 

tc

Switching and start-up time of LES

 

 

 

 

 

 

tL

Time limit of LES

 

 

 

 

 

 

Component out on maintenance indicated thus (")

 

 

 

 

 

 

 

 

 

 

 

 

 

Indices for TYPE 5 events are: X" = X; + X; U" = X r + X; r; r" =

These techniques are employed to reduce the amount of computation needed to obtain the solution to the load flow problem.

Nodal analysis

Figure 2.46 illustrates how a voltage source may be replaced by an equivalent current source.

In Fig 2.46 (a) we have voltage source E with admittance Y connected to node I having voltage to ground of V. It follows that V = E — I/Y and I =

as shown in Fig 2.46 (b). If we define II = YE and 12 = YV, this is shown in Fig 2.46 (c) and the voltage source has been replaced by a current source.

In Fig 2.46 (d), we have a voltage source E and an admittance Y between nodes 1 and 2, then

E i/Y. If now we define = —EY and 12 = + EY, this is shown in Fig 2.46 (e), and the voltage source has been replaced by two current sources.

To illustrate the formation of the basic equations used in nodal analysis, consider the simple network shown in Fig 2.47. Figure 2.47 (a) illustrates a simple power system which translates to the Fig 2.47 (b) network.

The network has a current injected into each node, each node is connected to at least one other node,

and at least one node is connected to ground (three nodes are connected to ground, in this example).

Using Kirchoff's first law at each node, we have the following equations:

Yl (v a — vb) +

Y4 v a = ia

Y2 (vb y e )

Yl(v a — vb) + Y5 vb = lb

Y3 (v c — vd) — Y2(vb — vc) = ic

Y6 Vd — Y3 (Vc

Vd) = id

expressing the admittance coefficients in matrix form, the voltages as a vector and the currents as a vector gives:

yi + Y4

—Y1

0

0

 

v a

 

la

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—Y1

Yl+Y2+Y5

—Y2

0

 

vb

 

ib

 

 

 

 

 

Vc

 

 

0

—Y2

Y2 +Y3

—Y3

 

 

 

 

 

 

 

 

 

 

 

0

0

—Y3

Y3+ Y6

 

vd

 

id

 

 

 

 

 

 

 

 

or, in short form

YV - 1

 

 

 

 

 

121

Electrical system analysis

Chapter 2

 

 

TABLE 2.13

Overlapping forced outages — third order (TYPE 5)

 

 

 

 

 

 

 

 

 

Residence time

(TYPE 5)

 

Ev en t

 

Failure rate

 

 

 

 

 

 

 

 

 

 

 

 

Restoration by

 

Restoration

 

 

 

 

 

 

 

 

 

component

 

from LES

 

 

 

 

 

 

 

 

 

repair

 

 

 

 

 

 

 

 

X

 

 

 

rA,

 

r13

 

 

 

 

 

 

 

 

 

 

 

 

(a)

1,

2,

3

 

XIX2riX3rir2/(ri

4-

r2)

 

Up

rir2t e /(rir2 + rite + r21 c )

(b)

1,

3,

2

 

Xik3r1X2rir3/(r1

+ r3)

 

1/A

rir3t e /(rir3 + rit e 4- r3( c )

(c)

2,

3,

1

 

X2X3r2Xtr2r3/(r2

-I- r3)

 

Iiii

r2r3tc/(r2r3 + r2tc + r3tc)

(d)

3,

2,

!

 

X3X2r3X1r2r3/(r2 + r3)

 

1/y

r2r3t c /(r2r3 + r2te + r3tc)

(e)

2,

1,3

 

X2X ir2X3r2r1/(r2

t

rt)

 

1/i,t

rir2t e /(rir2 + ritc + r2tc)

(f)

3,

1,

2

 

X3X I r3X2r3ri/(r3 + ri)

 

1/A

r3ritc/(r3ri + 1- 3te + rit c )

(g)

1,

2 + 3

 

X1X23r

 

 

 

I/A

rit c /(ri + t e )

(h)

2

+ 3, 1

 

X23X1r2r3/( - 2 + r3)

 

 

1/y

r2r3t e /(r2r3 + r2tc + r3tc)

 

 

 

 

 

 

 

For each event: average repair time r = rA P + rg (1—P)

 

 

 

 

 

 

 

 

Where:

=

 

+11 2 +143

 

 

 

 

 

 

 

 

P

=

e AtL

 

 

 

 

 

 

 

 

Xn

 

Failure rate of component

 

 

 

 

 

 

 

 

X23 = CMF rate for components 2 and 3

 

 

 

 

 

 

 

 

rn

 

Repair time of component 'n'

 

 

 

 

 

 

 

 

Le = Switching and startup time of LES

 

 

 

 

 

 

 

 

tL

 

Time limit of LES

 

 

 

 

 

 

 

 

 

 

 

Indices for TYPE 5 events are:

 

 

 

 

 

 

 

 

 

x5 = EX,

U5 = EXr, r = U5/X5

 

 

 

 

 

 

 

a

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

It will be seen that nodal admittance matrix Y can be written by inspection and that, for larger networks, the matrix will he highly sparse (i.e., contains a high proportion of zero elements). The properties of the matrix are:

It is a square of order (n x n), where n is the number of nodes in the network.

• It is symmetrical (Yik = Yki, general element •th row, k'h column).

It is complex, i.e., it may contain conductance and susceptance terms.

Each diagonal element Ykk is the sum of the admittances connected to node k, including the branches to ground.

Each off-diagonal element Y lk is the negative of the branch admittance between nodes i and k.

In large networks, many off-diagonal elements will be zero, i.e., the matrix is sparse.

The matrix equation YV = I defines the injected currents in terms of network admittances and voltages. Mathematically speaking, the set of equations YV = I may or may not have a solution, depending on the characteristics of matrix Y. No solution is possible if matrix Y is singular and thus has no inverse. This happens if the system has no ground connections. If matrix Y is not singular, the ease of solution depends on whether the matrix is well conditioned; if it is, computational rounding errors will not accumulate during the numerical process. Fortunately, the admit-

122