Induction Motor Load Dynamics: Impact on Voltage Recovery ...
Induction Motor Load Dynamics: Impact on Voltage Recovery ...
Induction Motor Load Dynamics: Impact on Voltage Recovery ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
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