Induction Motor Load Dynamics: Impact on Voltage Recovery ...
Induction Motor Load Dynamics: Impact on Voltage Recovery ... Induction Motor Load Dynamics: Impact on Voltage Recovery ...
2 representation. Section V presents some preliminary results with an example test system that comprises induction motor loads. Finally, section VI concludes the paper. II. PROBLEM STATEMENT The problem of transient voltage sags during .disturbances and recovery after the disturbance has been removed is quite well known. The importance of the problem has been well identified; its significance is increasing especially in modern restructured power systems that may frequently operate close to their limits under heavy loading conditions. Furthermore, the increased number of voltage-sensitive loads and the requirements for improved power system reliability and power quality are imposing more strict criteria for the voltage recovery after severe disturbances. It is well known that slow voltage recovery phenomena have secondary effects such as operation of protective relays, electric load disruption, motor stalling, etc. Many sensitive loads may have stricter settings of protective equipment and therefore will trip faster in the presence of slow voltage recovery resulting in loss of load with severe economic consequences. A typical situation of voltage recovery following a disturbance is illustrated in Fig. 1. Note there is a fault during which the voltage collapses to a certain value. When the fault clears, the voltage recovers quickly to another level and then slowly will build up to the normal voltage. The last period of slow recovery is mostly affected by the load dynamics and especially induction motor behavior. The objective of the paper is to present a method that can be used to study voltage recovery events after a disturbance. More specifically the problem is stated as follows: Assume a power system with dynamic loads, like, for example, induction motors. A fault occurs at some place in the system and it is cleared by the protection devices after some period of time. The objective it to study the voltage recovery after the disturbance has been cleared at the buses where dynamic or other sensitive loads are connected and also determine how these loads affect the recovery process. This paper proposes a hybrid approach to the study of voltage recovery that is based on static load flow techniques taking also into account the essential dynamic features of the load. This approach provides a more realistic tool compared to traditional load flow, avoiding however the full scale transient simulation which requires detailed system and load dynamic models. III. SINGLE-PHASE QUADRATIZED POWER FLOW (SPQPF) WITH INDUCTION MOTOR REPRESENTATION A. Overview of Single Phase Quadratized Power Flow The proposed system modeling is based on the single phase quadratized power flow. The idea of this power flow model is to have a set of power flow equations of degree no greater than two, i.e. have a set of linear or quadratic equations. This can be achieved without making any sort of approximations, so the power system model is an exact model. Since most of the equations turn out to be linear and the degree of nonlinearity for the nonlinear equations is restricted to at most two, this results in improved convergence characteristics in the iterative solution and therefore improved execution speed at no expense in the accuracy of the solution. Voltage (pu) 1.00 0.95 0.90 0.85 0.80 0.75 0.70 0.65 0.60 Fault Fault Cleared
3 Application of the connectivity constraints (Kirchoff’s current law) at each bus yields the quadratized power flow equations for the whole system: X 0 Y X X 0 T T F1 X F2 X b G( X ) , (2) where X : system state vector, Y : linear term coefficient matrix (admittance matrix), F : quadratic term coefficient matrix, i b : driving vector. The solution to the quadratic equations is obtained using the Newton-Raphson iterative method: 1 X X J X ) G( X ) (3) n n1 ( n1 n1 where n : iteration step, J ( ) 1 : Jacobian matrix at iteration n 1. X n The Iterative procedure terminates when the norm of the QPF equations is less than a defined tolerance. Therefore, the SPQPF equations G ( X ) 0 comprise a different mathematical system of nonlinear algebraic equations compared to the traditional power flow equations. The state vector consists of the real and imaginary part of the voltage at each bus and of additional internal state variables for each device. Some of these internal states are the additional variables introduced for the quadratization of the equations. The system G ( X ) 0 consists of the current balance equations at each bus, plus additional internal equations for each one of the nonlinear devices that exist in the system. Most of the equations are linear equations. All the nonlinear equations are of order at most quadratic. B.
- Page 1: 1 Induction <stron
- Page 5 and 6: 5 quadratic in the torque equilibri
- Page 7 and 8: 7 Motor active pow
2<br />
representati<strong>on</strong>. Secti<strong>on</strong> V presents some preliminary results<br />
with an example test system that comprises inducti<strong>on</strong> motor<br />
loads. Finally, secti<strong>on</strong> VI c<strong>on</strong>cludes the paper.<br />
II. PROBLEM STATEMENT<br />
The problem of transient voltage sags during .disturbances<br />
and recovery after the disturbance has been removed is quite<br />
well known. The importance of the problem has been well<br />
identified; its significance is increasing especially in modern<br />
restructured power systems that may frequently operate close<br />
to their limits under heavy loading c<strong>on</strong>diti<strong>on</strong>s. Furthermore,<br />
the increased number of voltage-sensitive loads and the<br />
requirements for improved power system reliability and power<br />
quality are imposing more strict criteria for the voltage<br />
recovery after severe disturbances. It is well known that slow<br />
voltage recovery phenomena have sec<strong>on</strong>dary effects such as<br />
operati<strong>on</strong> of protective relays, electric load disrupti<strong>on</strong>, motor<br />
stalling, etc. Many sensitive loads may have stricter settings of<br />
protective equipment and therefore will trip faster in the<br />
presence of slow voltage recovery resulting in loss of load<br />
with severe ec<strong>on</strong>omic c<strong>on</strong>sequences. A typical situati<strong>on</strong> of<br />
voltage recovery following a disturbance is illustrated in Fig.<br />
1. Note there is a fault during which the voltage collapses to a<br />
certain value. When the fault clears, the voltage recovers<br />
quickly to another level and then slowly will build up to the<br />
normal voltage. The last period of slow recovery is mostly<br />
affected by the load dynamics and especially inducti<strong>on</strong> motor<br />
behavior.<br />
The objective of the paper is to present a method that can<br />
be used to study voltage recovery events after a disturbance.<br />
More specifically the problem is stated as follows: Assume a<br />
power system with dynamic loads, like, for example,<br />
inducti<strong>on</strong> motors. A fault occurs at some place in the system<br />
and it is cleared by the protecti<strong>on</strong> devices after some period of<br />
time. The objective it to study the voltage recovery after the<br />
disturbance has been cleared at the buses where dynamic or<br />
other sensitive loads are c<strong>on</strong>nected and also determine how<br />
these loads affect the recovery process.<br />
This paper proposes a hybrid approach to the study of<br />
voltage recovery that is based <strong>on</strong> static load flow techniques<br />
taking also into account the essential dynamic features of the<br />
load. This approach provides a more realistic tool compared to<br />
traditi<strong>on</strong>al load flow, avoiding however the full scale transient<br />
simulati<strong>on</strong> which requires detailed system and load dynamic<br />
models.<br />
III. SINGLE-PHASE QUADRATIZED POWER FLOW (SPQPF)<br />
WITH INDUCTION MOTOR REPRESENTATION<br />
A. Overview of Single Phase Quadratized Power Flow<br />
The proposed system modeling is based <strong>on</strong> the single phase<br />
quadratized power flow. The idea of this power flow model is<br />
to have a set of power flow equati<strong>on</strong>s of degree no greater<br />
than two, i.e. have a set of linear or quadratic equati<strong>on</strong>s. This<br />
can be achieved without making any sort of approximati<strong>on</strong>s,<br />
so the power system model is an exact model. Since most of<br />
the equati<strong>on</strong>s turn out to be linear and the degree of<br />
n<strong>on</strong>linearity for the n<strong>on</strong>linear equati<strong>on</strong>s is restricted to at most<br />
two, this results in improved c<strong>on</strong>vergence characteristics in<br />
the iterative soluti<strong>on</strong> and therefore improved executi<strong>on</strong> speed<br />
at no expense in the accuracy of the soluti<strong>on</strong>.<br />
<strong>Voltage</strong> (pu)<br />
1.00<br />
0.95<br />
0.90<br />
0.85<br />
0.80<br />
0.75<br />
0.70<br />
0.65<br />
0.60<br />
Fault<br />
Fault Cleared<br />
<str<strong>on</strong>g>Motor</str<strong>on</strong>g>s will trip<br />
if voltage sags<br />
for too l<strong>on</strong>g<br />
-1.00 -0.50 0.00 0.50 1.00 1.50 2.00<br />
Sec<strong>on</strong>ds<br />
Fig. 1. Possible behavior of voltage recovery during and after a disturbance.<br />
The first step in expressing the power flow equati<strong>on</strong>s in<br />
quadratic form is to avoid the trig<strong>on</strong>ometric n<strong>on</strong>linearities.<br />
This can be achieved by utilizing rectangular coordinates<br />
instead of the traditi<strong>on</strong>ally used polar coordinates for<br />
expressing the voltage phasors. Therefore, the system states<br />
are not the voltage magnitudes and angles, but instead the real<br />
and imaginary parts of the voltage phasors. This results in a<br />
set of polynomial equati<strong>on</strong>s. If the degree of n<strong>on</strong>linearity of<br />
these equati<strong>on</strong>s is more than two, then quadratizati<strong>on</strong> of the<br />
equati<strong>on</strong>s can be achieved by introducing additi<strong>on</strong>al state<br />
variables. Note that the quadratizati<strong>on</strong> is performed without<br />
any approximati<strong>on</strong>s, and the resulting quadratic model is an<br />
exact model.<br />
The system modeling is performed <strong>on</strong> the device level, i.e.<br />
a set of quadratic equati<strong>on</strong>s is used to represent the model of<br />
each device. A generalized comp<strong>on</strong>ent model is used,<br />
representing every device, which c<strong>on</strong>sists of the current<br />
equati<strong>on</strong>s of each device, which relate the current through the<br />
device to the states of the device, al<strong>on</strong>g with additi<strong>on</strong>al<br />
internal equati<strong>on</strong>s that model the operati<strong>on</strong> of the device. The<br />
general form of the model, for any comp<strong>on</strong>ent k , is as in (1)<br />
k<br />
i<br />
<br />
Y<br />
0<br />
<br />
k<br />
x<br />
k<br />
x<br />
<br />
x<br />
<br />
<br />
<br />
k T<br />
k T<br />
k k<br />
F x <br />
1<br />
<br />
k k<br />
F x b , (1)<br />
2<br />
<br />
<br />
<br />
k<br />
i is the current through the comp<strong>on</strong>ent,<br />
x k is the<br />
where<br />
vector of the comp<strong>on</strong>ent states and b k the driving vector for<br />
each comp<strong>on</strong>ent. Matrix Y k models the linear part of the<br />
comp<strong>on</strong>ent and matrices F the n<strong>on</strong>linear (quadratic) part.<br />
k<br />
i