Interactive Simulation Of Power Systems: Etap Applications And Techniques

INTERACTIVE SIMULATION OF POWER SYSTEMS: ETAP APPLICATIONS AND TECHNIQUES Keith Brown Student Member, IEEE Ooeration Tichnolow. Inc. Iriine, California --

Herminio Abcede, Member, IEEE Flour Daniel Inc. Irvine, California

Farrokh Shokooh Senior Member, IEEE Operation Technoiogy, Inc. Irvine, California

Gary Donner Member IEEE. Texaco Refiding and harketmg Los Ange1es;California

Abstract

base including (optionally) load flow and/or short circuit results.

For many years electrical engineers have relied on the power (but not necessarily the convenience) of mainframe computers to analyze and design power systems. Recent advances in microcomputers have brought the power of mainframe computing to the desktop, paving the way for straight-forward, easy to use engineering applications which are designed especially for the personal computer.

ETAP, while capable of handling 1000 buses, contains a load schedule program which tracks up to 10,000,000 load items, and reports the voltage and short-circuit current at the terminals of each load item. This capability makes ETAP suitable for large industrial facilities, as well as utility systems. ETAP 5.5 includes the following programs :

System studies are an integral part of power system engineering and design. A structured computer program that uses technically correct models, employs a user-friendly interface, uses a common data base, and traps user errors is a powerful tool which greatly enhances the engineer's efficiency and productivity. ETAP is an engineering design and analysis program which satisfies these criteria. In addition, ETAP performs numerical calculations with tremendous speed, automatically applies industry accepted standards, and provides easy to follow output reports. This paper discusses solution techniques used by ETAP for various problems, as well as other features that make ETAP unique. Introduction The Electrical Transient Analyzer Program, commonly known as ETAP, began as a mainframe program and was rewritten for the PC as an interactive power system analysis and design tool. Significant effort has been made to insure that ETAP will provide the user with the highest possible technical accuracy, while maintaining a thoroughly user-friendly format. The common data base exemplifies the interactive nature of ETAP. Separate data editors for Bus, Branch, and Machine data allow the user to model the system in a common data base. User-edited libraries provide typical data which can be substituted into the data base upon request. When system studies are to be performed, ETAP automatically extracts the necessary parameters from the common data base. An interface with AutoCad is provided to automatically produce a detailed one-line diagram of the network modeled in the data

*

* *

* * * * * * * * *

One Line Diagram Load Flow Short Circuit Dynamic Stability Motor Acceleration Motor Starting Cable Derating Cable Pulling Ground Grid Design Induction Machine Parameter Estimation Induction Machine TorqueISlip Curve Load Schedule Common Data Base

A complete system analysis includes load flow, short circuit, and transient stability studies. Because of the unique data sharing capability of ETAP's programs, the data for an entire system needs to be entered only once with three user-friendly editors (Bus, Branch, and Machine). By not having to continually enter data, the user is free to concentrate his attention on analytical tasks, rather than bookkeeping. Errors in data input are thus reduced and more meaningful study results are obtained.

!

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

1

Often, when conducting system studies, the engineer does not have a complete set of input data, i.e. data is not available for certain machines, transformers, etc. The ETAP data editors are equipped with a large set of typical data for these situations. Based on information already entered (such as voltage rating), ETAP will automatically substitute typical data into the appropriate data fields upon request.

is an interface with AutoCad whereby the user can create a detailed one-line diagram without the assistance of a draftsman or CAD operator. The program automatically extracts data from the common data base and generates the one-line diagram, optionally indicating load flow and/or short circuit study results on the drawing.

As

indicated in Figure 1, ETAP data editors accept raw manufacturer data in standard NEMA form, as well as in per unit values, giving the user added flexibility. XFRW 13 10,000 YVA

New XFRW

-

7.500 YVA 13.80-2.40 k V

XFRY TRB 7.500 YVA 13.80-2.40 k V

12.0 L'R

14.0 X/R

Figure 2. using a new program often requires spending a considerable amount of time flipping through manuals to figure out what the computer is asking for in a certain entry field. On-line tutorials in the ETAP editors all but eliminate this problem. Each entry field has its own tutorial which is printed at the bottom of the screen. Another useful feature is the comprehensive error trapping that is continually occurring within the program while data is being entered. If the editor detects an unreasonable input, especially one which will result in meaningless study results, the user is immediately alerted, and given instructions on how to correct the mistake. There are actually four levels of error detection in ETAP. As described above, during data entry, ETAP checks each data input for validity. When exiting a data editor, ETAP again checks the data for errors. A third error check is made before any study is run, checking the validity of the data extracted from the common data base. Finally, run-time errors are trapped during program execution.

ETAP

contains an on-line power conversion calculator which can be invoked from the data editors by pressing a function key. Real and reactive power, voltage, current, and power factor can be calculated and substitute& directly into the editors. A

graphic representation of the network modeled in tha common data base can be generated One-Line Diagram. This program w i n g ET-'.

ETAP One-Line Diagram Load Flow

Load flow studies are essential in the planning and operation of electrical systems. The results obtained from a load flow study (in conjunction with short-circuit study results) are used to size capacitors, feeders, transformers, and current-limiting reactors. Whether designing a new system, or analyzing an existing one, factors such as voltage drop, load capacity, power factor constraints, steady-state stability limits, transformer tap settings, and generator excitation levels must be considered. The load flow study aids in the evaluation of these factors, and provides a convenient way of predicting the effects of system modification and expansion. The two most commonly used methods for calculating load flow are Gauss-Seidel and NewtonRaphson [1,2]. The Gauss-Seidel method starts with an initial guess of bus voltages. Using the scheduled power injected into the buses as well as the most immediate voltages of the other buses, the voltage at each bus is recalculated. These corrected voltages are then used to recalculate the voltages at all buses except the swing bus, continuing until all bus voltages have been evaluated, thus completing the first iteration. Iterations continue until all voltage corrections are less than a specified precision. A modification to this method is the Accelerated Gauss-Seidel method, which uses an acceleration factor tn a p e d up the convergence of the solution.

1931

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

Like Gauss-Seidel, the Newton-Raphson method is an iterative algorithm for solving a set of nonlinear simultaneous equations. At each iteration, the Jacobian matrix is used as a linear approximation for calculating a new set of bus voltages. This results in a quadratic convergence of the solution. However, this method is sensitive to the initial guess of bus voltages, as well as the level of system nonlinearity. To overcome the initial condition dependency, the Newton-Raphson method is preceded by a few iterations of Gauss-Seidel. Due to the nature of the Newton-Raphson method, extremely large humps in the admittance matrix could delay or even prevent solution convergence. Depending on the system configuration and impedances, the Newton-Raphson and Gauss-Seidel methods converge at different speeds. It is therefore advantageous to have both methods available. The load flow analysis package in ETAP provides both the NewtonRaphson and the Accelerated Gauss-Seidel methods, providing the flexibility needed to model loop and radial systems efficiently. Since some systems may be solved even faster by raising or lowering the acceleration factor, ETAP allows the user to adjust this factor. Other adjustable solution controls include the precision and the maximum number

of iterations. All system parameters are extracted from the ETAP common data base. The standard output report produced by ETAP lists all input parameters on separate pages: bus input parameters, linejcable parameters, and transformer and reactor parameters. The load flow results show bus voltages, loads, generation, and power flows, as well as transformer tap setting. Equipment exposed to undervoltage or overvoltage can be severely damaged and can prevent the efficient operation of machinery. Overloaded cables and transformers have considerably reduced life expectancies. When conducting a load flow study, these conditions are of paramount importance to the engineer. ETAP provides four optional summary reports to bring this information to the engineer's attention:

* * * *

Overvoltage Buses Undervoltage Buses Overloaded Branches Branch Losses and Voltage Drop

The loading in the shows a

user simply specifies the voltage or limits, and ETAP lists all violations respective summary reports. Figure 3 sample ETAP load flow output report.

Motor Load ---------__------_____MU

----_

Mvar

_____

* 1

Main A

Gen.

13.80

100.0

1.3

25.00

7.43

.oo

.oo

2

Rain B

Load 13.80

99.4

-.5

.OO

.00

.oo

.oo

*

.00

-1.97

13.80 100.0

1.6

25.00

8.27

.oo

.oo

.OO

.OO

4 Sub. 2-A Load '2.40 100.9

-.1

.00

.oo

2.45

1.52

.oo

.00

-1.6

.00

.00

3.34

2.07

.OO

.OO

.00

.00

3.82

.00

.00 .00

3

Rain C

Gen.

5 Sub. 2-B Load

2.40 100.2 2.40 100.7

-.5

Load

4.16

99.7

-.7

Load

4.16

99.4

-.5

.00 .00

9 Linc end Load 13.80

99.9

1.6

.00

6 Sub.6AlB Load 7 FCCU 8 FCCU

-A -B

Figure 3 .

2 13 1.70 2.15 2 100 -2.25 -1.62 2 100 -2.24 -1.62 2 5 2.61 1.99 2 8 .oo .oo 2 6 .19 1.08 3 14 3 9 3 100 3 100 3 6

3.76 1.95 7.80 7.84 3.65

2.29 1.29 1.61 1.62 1.46

-2.500

.o 17.1

-2.500

85.5 83.5 97.9 -.625 97.9 -.625 92.8 -2.500

1 -3.19 -1.71 .74 .19 5

862. 88.1 182. 96.9

5 5

4 -.74 -.17 2 -2.60 -1.90

182. 97.5 774. 80.7

2.500

2 -.19 -1.07 3 -3.63 -1.30

259. 17.2 922. 94.1

2.500 2.500

835. 100.0

2.37

.OO

6 6

6.00

.OO

.OO

.OO

7

1

-6.00

.OO

.oo .oo

.oo

.oo

8

2

.oo

.oo

.oo

.oo

.oo .oo

Load

184. 98. 333. 335. .%11

62.0 81.1 81.1 79.6

4 4

.OO

ETAP

115. 117. 116. 138. 0. 46.

9

v

3 -1.95 -1.28 10 1.95 1.28

0.

.o

98. 83.5 98.

Flow Output Report

1932

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

83.5

2.500

Short-circuit Short-circuit studies are performed for an electric power system in order for the engito determine an adequate protection neer scheme for the system and to coordinate the operation of protective devices during fault conditions. Almost instantaneously, shortcircuit currents produce powerful magnetic forces and intense heat in the power system, which can result in considerable damage to affected equipment if not promptly interrupted. These values of short-circuit currents must be determined by the engineer to ensure that the short-circuit ratings of all equipment are adequate to handle the currents available at their locations. In addition, the engineer can use data about these currents in conjunction with timecurrent characteristics of protective devices to determine the required relay settings and fuse ratings to achieve coordinated operation during faults. A properly coordinated protection system will generally result in fast and selective isolation of faulted equipment with minimum equipment outage. Calculation of short-circuit currents for industrial power systems tend to be more complex because of the mixture of sources contributing currents to the fault. In a typical modern industrial system the basic sources of fault currents are the utility, the in-plant generation, and synchronous and induction motors. These sources contribute additional exponentially decaying currents which make fault current magnitudes at Various locations time dependent. ANSI, IEEE, and IEC provide recommendations [3,4,5] as to how the decay of the ac and dc components of fault currents from various sources should be considered in fault calculations. The ETAP Short-circuit program is a userfriendly, fully interactive program for calculating rms values of fault currents corresponding to different time periods after fault occurrence. Specifically, it calculates the 112 cycle, 1.5 to 4 cycle and 30 cycle m s symmetrical fault currents for threephase, line-to-ground, line-to-line and line-to-line-to-ground faults.

For calculating device duties,

provides two methods for determining X/R ratios. The first method, which is in full compliance with ANSI and IEEE Standards, finds the equivalent resistance and reactance of the entire system, resulting in a single valued X/R ratio for a given fault location. The second method differs from the first in that it combines the individual branch current contributions (each with a separate X/R ratio) into a total value. This method, assumes the individual branch unique to ET-, contributions to be exponentially decaying currents, with time constants corresponding X/R ratios of the to the equivalent individual branches. For each time period, these exponential currents are evaluated and summed to get the total fault duty at a bus. The IEEE method was designed for simple manual calculations. As a result, it yields conservative results to insure safety. With the power of today's computers, ETAP is able to track individual fault contributions and their rates of decay. ETAp

By doing so, a more accurate representation of the power system is achieved. The user has the option of choosing either method. ETAP

provides two methods for selecting the ratio of individual machines. The first method, Variable X/R Ratio, calculates the X/R ratio using the subtransient reactances for 112 cycle currents, and transient reactances for 1.5 to 4 cycle currents. The second method, Fixed X/R Ratio, calculates X/R ratio using the subtransient the reactance for both 112 cycle and 1.5 to 4 cycle fault currents. This method is consistent with the IEEE adjustment tables which consider the X/R ratio a measure of the time constant of the exponentially decaying dc component during a fault at the machine terminals. Note that the fixed X/R ratio does not imply a varying armature resistance, but is included merely for consistency with the way the multiplying factor tables have been formed. X/R

For branches, the impedances of transformers and current-limiting reactors are adjusted if tolerances are specified. ETAP automatically uses the lower impedance values for conservative results. The reactance representation that ETAP uses for the different sources in calculating the 112, 1.5 to 4 and 30 cycle fault currents are summarized in Table 1, where XI1 represents utility short-circuit reactance and machine subtransient reactance, and represents machine transient reactance. This reactance selection is automatically performed by ETAP.

Contributions From

2 Cycle

Cycle

Cycle

Synchronous Generators

X '1

X"

X'

Synchronous Motors

X"

X'

Induction Motors

X"

X'

utilities

Table 1. Recommended reactance representation for utilities and machines in short-circuit calculations.

Table 2 presents typical reactances that the program will automatically use for induction motors of varying sizes if the user decides to use this feature of ETAP [ 6 ] . As in the other programs, the user can access typical data for use by the short-circuit program by pressing a function key during data entry. Table 2 presents a locked-rotor current (LRC) of 650% for medium voltage and 600% for low voltage motors. Using these LRC's, the corresponding subtransient and transient reactances for the given HP size and speed ratings are reported based on ANSI and/or IEEE recommendations.

1933

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

The short-circuit results are printed in a format that can easily be understood by even the novice user. Selection of appropriate breaker duties and bus bracing for each location in the power system is directly obtained from the printout. ETAP applies all appropriate ANSI, IEEE, IEC, and UL standards and procedures for device duty selection. For relay setting and coordination, fault contributions from all branches in the network to a given fault location ?ay be calculated and printed. This calculation is facilitated by an efficient, fast, and reliable algorithm using the bus impedance matrix with nodal injection method. Table 3 summarizes the application of the ETAP Short-circuit program output data.

@ 41800 rpm

X'=3.0

X"

Table 2. Typical parameters recommended for induction motors for use in short-circuit calculations (if actual data is unavailable). Except for a LRC of 650 % for MV motors, all parameters are based on ANSIIIEEE standards.

112 Cycle Currents

Mv

When calculating 112 cycle and 1.5 to 4 cycle fault currents for ac high-voltage circuit breakers, synchronous generators are modeled with the same reactance (subtransient). Note that this results in the same generator rms symmetrical contribution for both 112 cycle and 1.5 to 4 cycle fault currents, i.e., no ac decay. However, the symmetrical interrupting fault current for each bus is calculated as a combination of the utility, generator and motor contributions. For this reason, the program determines from the network configuration if a generator contribution to a given fault location is local or remote based on the amount of impedance separating the generator from the fault location. Therefore, for a given fault, the percentage of local and remote contributions is determined to calculate the NACD (No AC Decay) ratio. For the remote component, a multiplying factor with no ac decrement is applied while for the local component a factor with both ac and dc decrement is used.

Three-phase f a u l t a t bus nuker:

1

Main A

*.***********.****

,

13.80

CI Closing & Latching Capability

1.5-4 Cycle Currents zapability

LV Cf Interrupting Capability

Fuse

Interrupting Capability

SWGR/ Bus Bracing ICC Relab Instantaneous Settings

---

-

Figure 4 contains a sample ETAP short-circuit output report indicating fault contributions one bus level away from the faulted bus.

kV

bus

(

Prefault Voltage =

100.00 X )

------ ----No. Name No. --__- _--___-_- _ _ _ _ s=ii==

12 100 100 4 7 101

5A Cable Utility Utility Sub. 2-A FCCU - A Gen. #1

Manentary Duty:

1 1 1 1 1

' 1

.124 .238 .241 .071 .o60 .129

.94 82.14 82.14 23.30 20.3s 100.00

Spmetrical = 34.848 Aspmetrical = 58.077 Asyrmetrical = 55.756 CB Capability =

Interrupting Duty:

77.000

-1.890 -9.310 -9.310 -1.372 -1.418 -11.536

15.2 39.1 38.6 19.4 23.8

1.894 9.313 9.313 1.374 1.419

89.5

11.537

M

rms W rms M rms

M

(

MF = 1.667, X/R = = 1.60 a t X/R =

( MF

53.3 frm separate R B X networks ) 25.0 )

rms

syrmetrical = 33.936 Aspmetrical = 46.238

M rms M rms

( (

NAU) r a t i o = .66 ) MF = 1.362, X I R = 54.7

= 40.200

W rm6

(

5 cycle CB,

CB Capability

Figure 4.

Overcurrent Settings

Table 3 . Device duty and usage of fault current contributions.

To BUS

Contritution F r m Bus

30 Cycle Currents

total

frm separate R B X networks ) current r a t i n g bnsis )

ETAP sample short-circuit report. 1934

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

mansient (Dynamic) Stability Over the years, the size and complexity of utility power systems have increased tremendously with the additions of large power plants and long distance high-voltage ties between utility power pools. Similarly, industrial power systems, particularly oil refineries, have grown in size and complexity with the popularity of in-plant cogeneration systems and large synchronous and induction motor drives. As these power systems grow bigger and become more complex, the need to predict system behavior more accurately, both dynamic and steady-state, has also increased. In this paper! the words transient and dynamic stabilities are used interchangeably since both are commonly used to indicate studies which include the machine swing equations.

For a typical industrial power system consisting of induction and synchronous motors, in-plant synchronous generators, and utility interties, a power system engineer using ETAP investigates two types of system instabilities. The first type of transient stability study investigates the system behavior during the first cycle of the machine power angle swing. This study is commonly known as the first swing analysis. First swing analyses can cover disturbances such as faults, generator rejection, abrupt dropping of a large load, and loss of utility ties. A simulation time of 1 second is typically adequate for this type of study. The second type of stability analysis is commonly referred to as multi-swing analysis and requires more sophisticated simulation models and numerical solution techniques to accurately predict system performance over a relatively longer simulation time. A

good stability analysis program employs technically accurate simulation models for various power system components, namely, generators, generator excitation systems, prime mover speed-governing systems, synchronous motors, induction motors, static loads and system equivalents. The ETAP stability program includes all of the above simulation models. Moreover, it uses efficient numerical techniques to solve the differential equations of the simulation models and the algebraic equations for the network. All relevant system parameters are extracted from the common data base, and the load flow program is automatically run prior to the stability solution. The ETAP stability program 1s fully menu-driven. An input menu screen with an on-line tutorial is provided for each power system component to be included in the study. Running the program and accessing the output are as simple as pressing a series of function keys. Although the ETAP stability program is designed for industrial power systems, it could also be used For most utility power systems, except for those with HVDC interties. The simulation models and the required parameters of the ETAP stability program are deecribed in the following paragraphs.

For

synchronoum machines, implemented. The simpler

two modelm model does

are not

include the effects of damper windings, while the more complex model does [ 7 ] . Machine parameters including subtransient, transient and synchronous reactances and time constants are used in the more sophisticated model. As such, this model is suitable for studies that concentrate on damping of machine oscillations. The effect of magnetic saturation is considered for both models. Total inertia constant for the generator and prime mover is also a required input. For excitation systems, different models are provided from which a user can select the type that best matches the actual exciter to be used [ 8 ] . These include: IEEE Type 1 Continuously Acting Regulator With Rotating Rotating Exciter System; IEEE Type 2 Rectifier Exciter With Static Regulator Static System With System; and IEEE Type 3 Terminal Potential.

-

-

-

The IEEE Type 1 excitation system is representative of the majority of modern systems now in service and presently being supplied. As variations of the Type 1 excitation system, Type 2 and Type 3 have similarities to Type 1. The Type 2 system applies to units with the main input to the damping loop provided from the regulator output rather than the exciter output. Although the regulator transfer function of the Type 3 is similar to Type 1, the regulator output is combined with a signal which represents the self-excitation from the generator terminals. This feature makes Type 3 different from Type 1. Three models are provided €or governorturbine systems to represent the most commonly used turbines and governing systems [ 9 ] . Type ST is used for steam turbine, Type GT for gas turbine, while Type GP models a general turbine prime mover application. The mathematical derivations of these models are consistent with the IEEE Committee Report on governors and turbines. Induction machines are represented dynamically by their equivalent circuit as obtained either from locked rotor test data or from manufacturer's performance curves. ETAP provides two induction machine models. The sinpler model (Figure 5) uses a single-rotor equivalent circuit with an internal voltage source behind the machine locked-rotor reactance and stator resistance. The slip of the machine is a variable and is calculated from the acceleration equation. Changes in rotor flux linkages are considered in the model. The second model includes the effects of rotor deep bars and double cage machines. Electrical dynamics in the rotor circuit are, therefore, represented more accurately with this model. Other parameters that are required for both models include the combined inertia constant of motor and load and the mechanical load torque-speed characteristic. In cases where the circuit parameters are not known, the parameters can be estimated from manufacturer's performance curves, namely, torque, current, and power factor versus speed curves, using the ETAP induction motor parameter estimation program.

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

out#

4 8.98B N

out mtor load:

P= 9-

ETAP provides a total of six time events for specifying the occurrence times of the network changes. The first time event is defaulted at time T=O.O while all the other time events are user-specified. The seventh time event corresponds to the total simulation time. Up to a maximum of 4 2 network changes can be specified per time event, which is generally adequate for simulation of most commonly occurring system disturbances.

czm P= 48)

2,458 N

1.serullR

out uoltage: Unag. =1W.93x uang. = -8.86 aeg.

9.W PaRr Factor llachlne I l24 llachlne I125

The available models for synchronous motors are identical to those for synchronous generators, however, the excitation system and speed-governing system are not represented. Field excitation is considered constant during the simulation. The required subtransient and transient parameters for the more sophisticated model allow the damping of machine oscillations to be investigated more closely as the system frequency changes during a transient. Other loads are modeled either as constant kVA or constant impedance loads. In the constant kVA model, the load volt amperes is held constant only within a specified range. When the voltage goes out of this range, the model automatically reverts to constant impedance. With constant impedance representation, the load real and reactive power vary as a function of voltage squared. In the ETAP stability program, motor representation can be switched from dynamic models to constant kVA and vice versa for a particular simulation by pressing a single function key. If the detailed response of a given motor is not important for the case being studied, then the constant kVA representation is used for that particular run. In general this results in less computing time. In the ETAP stability program, the swing bus is automatically represented as a constant voltage of constant frequency behind a fixed impedance corresponding to the fault availability at the bus. This is done by the program without user intervention. With this representation, the swing bus acts as a sink or source of power as the system oscillates after a disturbance. Generally, the presence of a swing bus in a given simulation results in a better transient stability margin for the system. Simulation of system disturbances in the ETAP stability program is carried out by inputting a series of switching events in the dynamic simulation menu. This menu is used to specify the network changes such as short circuit, load/generator rejection, circuit switching, and/or impact loading situations. The menus provide for entry of the type of switching along with their time of occurrences. The specific network changes that are allowed in ETAP are branch addition and deletion, and bus revision (swing, generator, load, delete, fault, and accelerate)

A print/plot menu is provided in the ETAP stability program for use in selecting the variables to be included in the program output report. The calculated values of these variables are printed out at an interval corresponding to the specified time step or multiples thereof. A plot of these values as a time function may also be included in the output if the user so specifies. The variables that the user can specify for inclusion in the program output report are:

* * * * * *

Synchronous machine power angles Generator exciter voltages Bus voltages and frequencies Induction motor slips Branch power flows Apparent impedances

The various study categories for which the ETAP stability program may be used are: Generator Parameter Specifications which may involve determining generator reactances, generator-turbine combined inertia, and unit transformer reactance to achieve a certain stability margin for the generator as part of a system under specified disturbances; Generator Control which may cover specific areas, namely, excitation response time, excitation voltage ceiling, and generator dynamic breaking for investigation relative to system transient performance; System Control which may deal wfth controlled separation of various islands of a utility system and also load shedding to prevent system collapse during an overload condition of an ascended subsystem; Plant Equipment which may cover various conditions relating to the starting of large motors and their stability during fault disturbances in the system; Transmission System Design which could address items such as line design, voltage level, requirements for intermediate switching stations, voltage support, series compensation and shunt compensation relative to transient stability performance of the system; Transmission Line Protection which may study breaker and relaying speeds, breaker schemes, automatic reclosing and transfer tripping relative to critical fault clearing time.

1936

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

user to accelerate motors under a variety of circumstances.

Motor Acceleration In large industrial power systems and power plants, acceleration or re-acceleration of large motor drives is always necessary. It is therefore crucial f o r the design engineer to predict accurately whether or not a large motor would be able to successfully start and to determine how long it would take to reach the steady state condition. If analysis indicates that operating problems would arise, the engineer can then identify necessary changes or adjustments to the system during its design stage. I f problems are not accurately predicted, then practical solutions at the later stage of plant testing may be limited and uneconomical. A computer program that is based on technically accurate simulation models for motor acceleration and re-acceleration analysis is therefore very useful to an engineer to help investigate all the W h a t ift1 scenarios early in the design of the system.

The required data €or this study are described previously under the Transient (Dynamic) Stability Program. The motor to be started is modeled dynamically, preferably using the induction motor double cage model. This is generally an accurate representation for any ac starting motor since the rotor reactance and resistance generally change with the machine speed. The torque-speed characteristics of the driven load are specified by indicating its shape and the machine loading at the steady state condition. ETAP uses a curve fitting program to form a 5th order polynomial for the load torque. For no-load starting, the final value of the load torque can be either zero or a small value to reflect motor friction and windage losses. ETAP can provide a compplete load flow report at each time event, i.e., at the first time event (T=O.O before starting), second time event, and so on including the end of simulation time. As described earlier in this paper, plots and print outs of the different motor variables as functions of time can be requested from the program printlplot menu to obtain a complete picture of the motor acceleration. A sample motor acceleration curve generated by ETAP is shown in Figure 6.

The ETAP stability program has the capability to simulate motor acceleration or reacceleration dynamically. Motors to be started are specified in the dynamic simulation menu as machines to be accelerated in a specific time event. In addition, other network changes can be specified at the instant of motor starting. This feature allows the

~ n d .motor Slip (X) 104 (3500 HP )

-

124 (1250 HP )

-

vs

125 (700

X

Time (Secwds) HP 1

-

+

33.3

.oooo .0500

.3500

.4000 .4500 .5500 .6000 .6500 .7000 .?so0

+

.woo

+ .a500 +

.woo .I

*

*

+

.E500 +

-X

1.2500 1.3000 1.3500

*

1.6500 +x

r.moo

t*

1.8ooo

it

1.7500 +x 1.8500 +x 1.9000 +x 1.9500

*I

Figure 6.

Sample Motor Acceleration Output Plot

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

Cable m a t i n g Electrical power systems are only as reliable as the cables that network them together. A complete analysis of a power system should therefore include an analysis of the cable ampacities. This analysis is complicated since the ampacity of a conductor varies with different installation configurations. Ampacity is defined as "the current in amperes a conductor can carry continuously under the conditions of use (conditions of the surrounding medium in which the cables are installed) without exceeding its temperature rating." Therefore, a cable ampacity study is the calculation of the temperature rise of the conductors in a cable system under steady state conditions. The National Electrical Code (NEC) [lo] now accepts the Neher-McGrath method [11] of determining the ampacities of conductors. NEC gives ampacity tables for a few specific configurations, but since most practical applications have yore complex configurations than those given in the tables, the engineer must turn to a computer program which calculates the temperatures of each conductor in the installation, and verifies that they do not exceed the rated maximum temperature. ETAP Cable Derating Program (ETAP-CD) is based on the Neher-McGrath method, and not onl-y calculates temperature rise in cables, but also contains options for sizing cables, and optimizing ampacities.

ETAP-CD contains many extra features tailored especially for the design engineer. User-friendly data entry menus as well as a cable library database make the cable derating process straight-forward and hassle free. Since the design engineer starts from scratch, simply finding the temperature rise of cables in a given installation is not sufficient. This is why ETAP-CD includes special options for sizing cables and optimizing ampacities. Finding ampacities which optimize the amount of current that can flow through a given underground installation would usually be done by trial and error, but with the Find Ampacities option, the engineer simply enters the maximum rated temperature for the cables in the duct bank, and ETAP-CD finds ampacities which will bring each conductor to the specified temperature, thus maximizing the current through the entire installation. Alternatively, it may be desirable to run each unique cable size at the same ampacity. ETAP-CD computes these ampacities just as easily. If the loading requirements are known for each cable, and the cable sizes have not been determined, selecting the Size Cables option tells ETAP-CD to choose adequate cable sizes from the library. User-friendliness is an important aspect that is often overlooked in engineering applications. Special care has been taken to make data entry as simple as possible in ETAP-CD. A separate menu was designed to aid in the entry of repetitive data which also calculates the spacing of the cables for direct substitution into the appropriate entry field. popup menus allow the user to select a run option (see Figure 7).

When an electrical current flows through a cable, it generates heat. The type of cable and its location in the installation are two of the factors which determine how many compon nts of heat generation are present, such as I R losses, sheath losses, dielectric losses, and so on. The heat flows from these sources through a series of thermal resistances to the surrounding environment [12]. The cable operating temperature is directly related to the amount of heat generated and the effective thermal resistance through which it flows. The Neher-McGrath method involves the application of a series of thermal equivalents of Ohm's and Kirchoff's laws to a thermal circuit. This circuit includes several parallel paths with heat entering at several points. The steady state Conductor temperature can be determined by the temperature differential created across the thermal resistances as the heat flows to the ambient temperature.

5

Important characteristics of the installation that determine thermal resistance between the cable and the ultimate surrounding environment include the type of conduit (steel, PVC, etc.), the depth of the underground installation, the thermal resistivity of the earth surrounding the installation, and the conductor insulation thermal resistivity. Factors that determine the external heat sources include the proximity of nearby cables and the temperature of the surrounding environment. The heat generated in each cable depends on the cable attributes, such as the type of conductor and the thermal resistivity of the insulation, as well as the dielectric losses.

-

ill

I

1 1

I

2 2 2 3 3

3

4 4 4

-

Clrarlt uller)*ulber or Iu=der Mber

lun

Figure 7.

ETAP-CD Cable Location Menu

The

output report can be printed by pressing single function key. All input parameters (including any data extracted from the library) are printed, followed by a structured output which graphically shows the location of the cables in the installation, making the report easy to use. Figure 8 shows a sample duct bank system, with the output report shown in Figure 9.

a

1938

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

Induction Machine Parmeter - t h a t i o n

Available circuit parameters of induction machines are usually sufficient for shortcircuit analysis only. Many steady state and transient studies require more sophisticated models which include the resistances and reactances representing the stator, rotor, and magnetizing branches. Additional parameters used in some models include the effects of saturation, hysteresis, eddy currents, deep rotor bars, and double rotor cages. Figure 8. Sample 3x5 duct bank installation COLUINS:

IOU 1

RW 2

2

--__---

3

4

5

350 138.5

350 146.9 75.0

350 176.3 75.0

350 110.3 75.0

149.3

350 124.8 75.0

350 157.6 75.0

_ _ _ _ _ _ _ _----------.-

Cable: Mp.: Tanp.:

500 209.8 75.0

500 175.9 75.0

Cable:

500 173.8

99.0

75.0

500 129.4 75.0

Mp.:

500 185.3

500 148.4

350 115.1

hap.:

75.0

75.0

75.0

Ap.: 7anp.: Cable:

ROU 3

1

_----__

75.0 350

75.0

An equivalent circuit representing the steady state behavior of a polyphase induction motor is shown in Figure 10 where the deep rotor bars and double cage effects are included in the rotor reactance and resistance as a function of slip and full-load values.

350 75.0

Figure 9. Output report from ETAP-CD using the Find Ampacities (Uniform Temperatures) option.

Figure 10. Equivalent circuit representing polyphase induction motors.

Cable pulling

In the design of underground distribution 8y8temSI the limiting factors imposed by cable pulling design criteria must be considered. For most underground installations, the number and location of splices, electrical manholes, pull boxes and other pull points in the system are known to be dependent upon the maximum allowable pulling lengths of cables. Therefore, in order to achieve an optimum design that will be economical and will not cause damage or distortion to the cables being pulled into the conduits, the practical application of cable pulling theory and design criteria is necessary. The ETAP Cable Pulling program, ETAP-CP, was developed using widely accepted cable pulling equations and criteria. It is user-friendly, fully interactive, and menu driven, and is certainly a handy tool for a cable run designer. The user can model three dimensional configurations by specifying lengths and slopes of pull segments, angles and radii of bends, coefficients of friction, etc. Cable libraries that can be modified and expanded by the user are provided. These include physical and electrical data for various cables of different voltage ratings. The ETAp Cable Pulling program is used to determine tensions and sidewall pressures that a cable is subjected to when pulled into conduits. A point-by-point calculation method ie performed at every conduit bend and pull point. The program is applicable to simple as well as complex conduit configurations where both the forward and reverse pulling tensions as well as sidewall pressures are calculated to determine the best direction of pull. Additionally, it checks for percent fill, jam ratio and other NEC and cable requirements. The program automatically flags any installation problems.

The task of obtaining machine design data can be difficult, if not impossible for some machines. Manufacturers do not always supply complete data, while owners may misplace data that has been supplied. For most small machines, the only available data appears on their nameplates. Additional data available for some larger motors are the steady-state performance curves, which include the torque, current, and power factor versus slip. The Parameter Estimation Program (ETAP-PE), assumes that the nameplate data along with two points from the motor performance curves are available. These points, which include the locked rotor torquelcurrent and breakdown torque, may be obtained from actual machine tests or from curves supplied by the manufacturer. Indeed, these points are often the only reliable information available from the performance curves since the rest of the curve might have been provided merely as an indication of the machine's generic behavior. The algorithm for estimating the equivalent circuit parameters [13] proceeds by first determining the parameters of a simplified equivalent circuit. Next the method utilizes exact machine equations to obtain a more accurate estimation. A sensitivity analysis of the performance characteristics with respect to the circuit parameters is then performed to obtain a higher degree of accuracy. illustrate the application of ETAP-PE, a three-phase, 1000 HP, 4.16kV, 60Hz induction motor was selected. This particular example 20 iterations of sensitivity required adjustments until the performance characteristics calculated from the estimated Circuit parameters were within 1 percent of the given values, as indicated in the computer output shown in Figure 11. TO

1939

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

Inprt Data ( 60 Hz ):

ETAP-GRD

provides

options

which

allow

the

engineer to optimize the design and to deterMax. Locked Rotor _ * _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ Torque _ _ _ _ __ _ _ _ . _ _ _ _mine ____ ---- - - economically sound implementathe most X T X PF X Eff. X X I X Slip X PF tion of a given configuration. For a fixed Rated ( F u l l Load)

______ _ _ _ _ _

_ _ H_ P_ _ _ 1000.0

kV

Poles

4.00

2

______ .___ _ _ __ _ _ __ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 89.0

1.50

Calculated Full Load Current =

126.0

F u l l Load Speed

= 3546.0

F u l l Load Torque

=

M a x i m Torque

96.0

220.0

450.0

35.0

15.0

A~ps

RPM

1481.3

ft-Lb

= 3258.8

ft-lb

number of rods, the program finds the optimum number of conductors. Based on the user-specified costs of conductors and rods (including installation costs), ETAP-GRD finds the optimum number of conductors and rods for the configuration. Both of these options find solutions which limit the step and touch potentials to acceptable values. The program also checks the ground conductors and rods for their capability to carry the ground fault current.

Estimated Parameters ( X ) :

Conclusion Engineers depend on continual development of software to keep in step with fast-paced technological advancement. ETAP's designers are professional electrical engineers who are dedicated to researching and implementing new methods as well as developing new applications which will enhance the productivity and capability of ETAP users.

X"

= =

XOC

= 357.73

X Open-circuit reactance

X/R

= =

15.22 13.81

X/R r a t i o (X"/Rs) a t f u l l - l o a d a t locked-rotor

Tdo' =

0.634 0.601

Rotor open-circuit time constant (sec.) a t f u l l - l o a d a t locked-rotor

0.053

Rotor resistance cage factor [ Rr = ( 1 + S*Kr)'Rr,s=o

22.85 21.80

=

Kr

=

-0.093

WX

X subtransient reactance a t f u l l - l o a d a t locked-rotor ( = Xlr )

Rotor reactance cage factor

I

t Xr = (1 + S*Kx)*Xr,s=o I

(a)

"1"1 1 -

Today's electrical engineer is limited only by the software available to him. ETAP is a technically accurate, user-friendly software package designed with the working engineer in mind. From the common data base input editors, to the clearly structured output reports, ETAP has pushed forward the frontier of desktop power system analysis software. References 1. Glenn W. Stagg and Ahmed H. El-Abiad, Computer Methods in Power System Analysis, McGraw-Hill, 1968.

x

1

C -

0

U

R P

R R

U

E-

E -

N

2. J. Arrillaga, C.P. Arnold, and B.J. Harker, Computer Modelling of Electrical Power Systems, John Wiley & Sons, 1983. 3. ANSIfIEEE C37.010-1979, IEEE Application Guide for AC High-Voltage Circuit Breakers Rated on a Symmetrical Current Basis.

1 -

4. ANSIfIEEE C37.50-1979, Guide for Calculation of Fault Currents for Application of AC High-Voltage Circuit Breakers Rated on a Total Current Basis.

Figure 11.

ETAP-PE example output report. Ground Grid Design

When a short-circuit occurs, large currents must be passed to the earth without endangering personnel in the vicinity. Based on ANSIfIEEE Std 80-1986 [14], the Ground Grid program, ETAP-GRD, calculates the required number of parallel ground conductors and ground rods for safety during a ground fault in the system. For a given grid configuration with a specified number of conductors and rods, the program calculates and reports the step and touch potentials and flags unsafe conditions.

5. ANSIfIEEE C37.13-1981, IEEE Standard for Low-Voltage AC Power circuit Breakers Used in Enclosures.

6. Walter C. Huening, Jr. , "Calculating short-circuit currents with contributions from Induction Motorsa1,IEEE Trans. Ind. App. Vol. IA-18, No. 2, MarfApr 1982. 7. Edward Wilson Kimbark, Power Stability: Synchronous Machines, Publications, 1956.

system Dover

on computer 8. IEEE committee report modelling of excitation systems, "Excitation system models for power system stability analysisIng IEEE Trans. Power App. Sys., Vol. PAS-100, No. 2, Feb. 1981.

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.

9. IEEE committee report for governors and turbines, "Dynamic models for steam and hydro turbines in power system studies," IEEE Trans. Power App. Sys., Vol. PAS-92, pp. 1904-15, Jul/Dec 1973. 10. National 70-1990.

Electrical

Code

(NEC), NFPA

11. J.H. Neher and M.H. McGrath, "The calculation of the temperature rise and load capability of cable systems, AIEE Trans. Power App. Sys. (pt. 111), Vol. 76, pp. 752-72, Oct. 1957. 12. Farrokh Shokooh and Harold M. Knutson, "Ampacity derating for underground cables,'* IEEE I K P S Conference in Baltimore, Maryland, May 1988. 13. Somchai Ansuj, Farrokh Shokooh, and Roland Shinzinger, "Parameter estimation for, induction machines based on sensitivity analysis," IEEE Trans. Ind. Appl., vol. IAS-25, pp. 1035-40, Nov/Dec 1989. 14. ANSI/IEEE Std 80-1986, IEEE Safety in AC Substation Grounding.

Guide for

1941

Authorized licensed use limited to: UNIVERSITY OF WESTERN ONTARIO. Downloaded on March 9, 2009 at 20:20 from IEEE Xplore. Restrictions apply.