Nedic Tech Rep 2002
Nedic Tech Rep 2002
Nedic Tech Rep 2002
Dusko Nedic
Transformers
Transformers are used to transfer power between different voltage levels or to regulate real or reactive flow through a particular transmission corridor. Most transformers come equipped with taps on the windings to adjust either the voltage transformation or the reactive flow through the transformer. Such transformers are called either load-tap-changing (LTC) transformers or on load tap-changing (OLTC) transformers. Another type of transformer is known as a phase-shifting transformer (or phase shifter). Phase-shifting transformers, which are less common than LTC transformers, vary the angle of the phase shift across the transformer in order to control the MW power flow through the transformer. This type of control can be very useful in controlling the flow of real power through a transmission system. The emphasis of this document is the modelling of different types of transformers in AC load flow. The first section is a short introduction in different types of transformer considered in such AC load flow, afterwards a comprehensive branch model and the expressions for line flows trough transformer are given. The summary of the previous work based on several references was given in the section Review of the Previous work. The succeeding sections contain the basic ideas, formulations, adjustment interactions and algorithm for tap adjustments based on AC sensitivity analysis. Small examples are given in the last two sections for the sake of practicality.
Transformers with fixed taps operate at the given off-nominal turns ratio and phase shift, and will remain fixed at those values during the entire power flow solution process unless manually changed by the user.
When on automatic voltage control, the transformer taps automatically change to keep the voltage at the regulated bus (usually one of the terminal buses of the transformer) within a voltage range between the minimum voltage and maximum voltage values. When on automatic reactive power control, the transformer taps automatically change to keep the reactive power flow through the transformer (measured at the from bus) within a user-specified range. When a transformer is on phase shift control, the transformer phase shift angle automatically changes to keep the MW flow through the transformer (measured at the regulated bus end) between the minimum and maximum flow values (with flow into the transformer assumed positive).
Ii 1:t y V'i
Ij
Vi
Vj
Figure 1. Transformer equivalent circuit The voltage and current ratio can be then defined as follows:
Vi : Vi' = 1 : t and I i : I i' = t * : 1 due to Vi* I i = Vi'* I i'
(T1)
Vi -is the complex voltage at the i end of the line i-j, Vi' -is the complex voltage behind the ideal transformer,
V j -is the complex voltage at the j end of the line i-j,
I i -is the complex current at the i end of the line i-j, I i' -is the complex current behind the ideal transformer,
I j -is the complex current at the j end of the line i-j,
The transformer equivalent circuit shown in Fig. 1. can be transformed to an equivalent circuit using the following equations:
I i = t * I i' = t * (Vi' V j ) y = t * (tVi V j ) y = t 2Vi y t *V j y ,
or in a matrix form:
I i t 2 y t * y Vi = y V j I j ty
(T2)
Ii
yt * /(yt)
Ij
yt *(t-1) Vi
y (1-t*) Vj
Figure 2. Comprehensive branch model for8 transformers Based on equation (T2) a comprehensive branch model is shown in Fig. 2. It should be noted that only phase shifter transformer has y ij y ji , while for all others types of transformer t * = t and consequently y ij = y ji . Besides, this branch model assumes that the transformer admittance is behind the off nominal side of transformer. Some other branch models are given in [1].
(T3)
Vi = Vi e ji , Vi = Vi , i = Vi ,
V j = V je
j j
(T4) (T5)
5
, V j = V j , j = V j ,
t = te j , t = t , = t ,
y = ye j = g + jb, = arctan b g.
(T6) (T7)
Substituting the complex variables with the polar coordinates given in equations (T4-T7), equation (T3) can be rewritten as:
j (i j +)
),
(T8)
{ } {
(T9)
{ } {
(T10)
S ji = V j2 y V j*Vi ty
or, in terms of real and reactive power as:
(T11)
{ } {
(T12)
{ } {
(T13)
(T14)
where x is the vector of independent variables (for example active and reactive power at a PQ bus) and y is the vector of dependent variables (for example voltage and angle for a PQ bus). This system of non-linear
equations can be linearised at a particular point x 0 by expanding it into a Taylors series and retaining only the first- order terms. Therefore,
g ( x, u, p) x
x0
x =
g ( x , u, p) u
x0
u ,
(T15)
where the vector of independent variables y is split into the vector of controllable variables u (for example tap changer values) and the vector of fixed parameters p . The changes in any desired function f ( x , y ) with respect to the single parameter change u i can be found from the total differential as follows [6]:
f = f ( x, y ) x
x0
x +
f ( x , y ) u
x0
g ( x , u, p) x
x0
x0
g ( x, u, p) ui
x0
x0
ui +
f ( x, u, p) ui
x0
ui
f f = ui x
J 1
x0
g ui
f ui
x0
(T16)
Using equation (T16) the change of control variable u i can be adjusted to achievestrike the desired change in function f if the sensitivity value
f is known. This sensitivity analysis with respect to a single u i
parameter change can be implemented in tap adjustments, selecting appropriate desired functions. Thus, in the case of the voltage-controlled transformers the desired function is the voltage at the controlled bus. For reactive power control transformers the desired function is the reactive power flow at the from bus. In both cases, the controllable variable u i is the tap changer value. The calculation of sensitivity shown in the previous equation requires the following calculation steps: 1. calculation of
f , x
f f and depends on the desired function choice as it will be shown separately for u i x
The calculation of
and j is used to control the voltage at bus j within a specified range. Therefore, the desired function is f = V j and the control variable is u i = t ij . Then, the calculation required by the steps 1,3, and 4 (see the previous section) is:
f = [0 0 L 1 0] = e Tj , where the non-zero entry refers to the position of the voltage V j x
in the state vector (vector of independent variables) of the Newton-Raphson iteration procedure.
0 M 0 PB i M t 2 PB 2tVi g + ViV j y cos (i j j ViV j y cos( j i ) t 0 g 0 = = M ui M QBi 2tV 2b + V V y sin ( ) i i j i j t ViV j y sin ( j i ) QB j M t 0 M 0 ) ,
(T17)
where PBi , PB j , QBi , QB j are the real and reactive balance equation at buses i and j, respectively. These balance equations give the total injection of real and reactive power into a bus, summing generation, load and line flows to/from the bus ( PBi / j = PiG j Pi L j Piinjj / / / and
f =0 u i
other hand, the calculation of the inverse Jacobian is time demanding, especially for large systems. This calculation can be avoided using the illustration of the
f x
x0
J 1
x0
g ui
x0
shown in Fig. 3. It can be seen that such structure has only one non-zero value in the vector
beneficial in the sense that the calculation of the inverse Jacobian matrix can be avoided. In essence, in Fig.3
we can observe that only one row of the inverse Jacobian matrix is required. If the dimension of the Jacobian matrix is N, than having found the lower (L ) and the upper (U ) triangular sub-matrices, an inverse matrix calculation would require N forward-backward substitutions ( J * J 1 = I LU * [ j1
j2 L jN ] = [e1 e2 L e N ] ). In each of these substitutions ( LU * ji = ei , i = 1, N ) one column ( ji , ) of the inverse
Jacobian matrix is determined. However, this time consuming calculation of each column to determine only one row of the inverse Jacobian can be avoided using the equation ( J 1 ) T = ( J T ) 1 . Therefore, instead of N forward-backward substitutions only one substitution ( U T LT * j row = e j ) will be required to determine the row j
j row that corresponds to the position of the independent variable V j (see Figure 3). The next simplification in j
the sensitivity calculation is related to the LU decomposition. At each iteration step in a Newton-Raphson procedure, LU decomposition is needed. Using the assumption that the change of the Jacobian matrix is not so dramatic between two subsequent steps in a Newton-Rapshon iteration procedure, already obtained L and U sub-matrices can be used for the calculation of the corresponding inverse Jacobian row. This assumption is even more sensible if one takes into account that the tap adjustments will take place only if the maximal mismatch in the iteration procedure is relatively small. A complete algorithm of the suggested tap adjustments is shown in Fig. 4. Special attention is paid to acceleration factors in order to prevent excessive tap movements or to avoid oscillations in the iteration procedure. These oscillations might appear if a tap change of a transformer cause unwanted change of another voltage controlled by another transformer. These adjustment interactions are summarised in [7], which is in the authors opinion is the best paper ever written on this topic. Therefore, the effect of a tap change can be sometime mitigated/aggravated by another tap change. The acceleration factors are equal to 1.0 if the normal tap change is required, and less than 1.0 for excessive tap changes. It should be noted that the impact of simultaneous tap adjustments on divergence has not been fully investigated. However, the testing has shown that a correct choice of the acceleration factors makes this algorithm very robust in terms of convergence. The tap adjustments based on this approach have been successfully tested on IEEE 14 bus system and the NGC power system (1100 buses, 1700 branches). An illustrative example is given at the end of this document, emphasising some crucial points of the suggested AC sensitivity approach.
Vj position in the state vector
og ou
i position in the state vector j position in the state vector
ej
-1
jj row
10
no
yes
Calculate the sensitivity for each transformer if its regulated voltage is not within the limits. If all the regulated voltages are within the limits, return to Newton-Raphson iteration procedure Using these sensitivites determine the new tap positions and update the admittances. An accelerating factor should be determined for each tap change in order to prevent some large tap movements. Return to Newton-Raphson iteration procedure.
11
(T18)
(T19)
There are a few changes with respect to the voltage control transformer sensitivity calculation. The vector
f has now four non-zero entries, and consequently the calculation of four rows of the inverse Jacobian is x
x0
J 1
x0
g ui
x0
og ou
i position in the state vector
J-1
Vi position in the state vector Vj position in the state vector
P load(MW) 0.0 95 55
Q load(MW) 0.0 43 45
Figure 5 A small power system It can be seen in Table 2 that the following scenarios were analysed: 1. LTC transformer between buses 3 and 2 controls the voltage at the bus number 3 within the interval 0.96 0.98. 2. LTC transformer between buses 3 and 2 controls the voltage at the bus number 2 within the interval 0.98 0.985.
V3 f = = 0.32023 is calculated in the second iteration of the N-R t 32 t 32
iteration procedure, because the controlled voltage V3 = 0.950782 was out of the specified range, and correction of
V3wanted = 0.00921767 was
14
t 32 = 0.028744 . After the tap position had been changed, a new iteration of the N-R procedure was carried
change is slightly larger than the wanted one, which is the consequence of the implemented linearisation. In the next iteration, the voltage at bus number 3 was V3 = 0.9599439 , which is still out of the specified range, but within an acceptable tolerance (0.002).
V2 f = = 0.318302 was calculated in the second iteration of the Nt 32 t 32
R iteration procedure, because the controlled voltage V 2 = 0.950174 was out of the specified range, and correction of V 2wanted = 0.0298259 was wanted. However, the sensitivity is now positive and consequently a positive tap change of t 32 = 0.0937 is required. After the tap position had been changed, a new iteration of the N-R procedure was carried out and the voltage at bus number 2 was changed to
V 2 = 0.98180 . In the
next iteration, the voltage at bus number 2 was V 2 = 0.9785541 , which is still out of the specified range, but within an acceptable tolerance (0.002).
It can be seen in the modified Table 2 that the following scenario was analysed: Reactive power control transformer between buses 3 and 2 controls the reactive power flow at bus 3 within the range 0.2 0.3.
f Q32 = = 2.580067 was calculated in the second iteration of the N-R t32 t32
iteration procedure, because the controlled reactive power flow Q32 = 0.00637 was out of the specified
wanted = 0.20637 was wanted. The required tap movement was t32 = 0.07998 and range, and correction of Q32
2 0.026756 0.0021399 0.021884 x 3 x wanted = 0.00175 = = 0.026923 t32 V2 0.3366308 0.026911 V3 0.336474
The change of independent variables calculated in the nest iteration of the N-R iteration procedure, after the tap change t32 took place was:
0.00229 0.00144 x = 0.02821 0.02544
and the reactive power flow at the from bus was improved to Q32 = 0.225 . No other change in the tap position was required during the iteration procedure.
16
References:
1. L.V. Barboza, H. H. Zurn, R. Salgado: Load Tap Change Transformers A Modelling Reminder, IEEE Power Engineering Review, February 2001. 2. R.N. Allan, C. Arruda: LTC Transformer and MVAR violations in the Fast Decoupled Load Flow, IEEE transaction, PAS-101,1982, pp. 3328-3332. 3. S. K. Chang, V. Brandwajn:Adjusted Solutions in Fast Decoupled Load Flow, IEEE Transaction on Power Systems, Vol. 3, No. 2, May 1998. 4. N. M. Peterson, W. S. Meyer:Automatic Adjustment of Transformer and phase Shifter Taps in the Newton Power flow, IEEE Transaction, PAS, Vol PAS-90, No. 1, January/February 1971. 5. B.Stott:Review of Load-Flow Calculations Methods, IEEE Proceedings, Vol. 62, pp 916-929, july 1974. 6. W.F. Tinney, H.W. Dommel: Steady State sensitivity analysis, 4th PSCC Grenoble, 1972. 7. S. K. Chang, V. Brandwajn: Solving the Adjustment Interactions in Fast Decoupled Load Flow,IEEE Transaction on Power Systems, Vol. 6, No. 2, May, 1991.
17