1 Matriz admitancia Y

BUS
Los componentes de un sistema de potencia se encuentran interconectados el ectricamente lo que
hace que un cambio de un par ametro de un elemento o de la potencia generada o demandada
en un nodo, produzca una variaci on de las corrientes, potencias o tensiones existentes en otras
partes del sistema. Un efecto que inicialmente es localizado puede generar efectos globales, y su
intensidad depende, entre otras cosas, de la distancia el ectrica existente entre los elementos. La
relaci on de interdependencia entre los diferentes elementos del sistema de potencia puede ser
adecuadamente caracterizada por la matriz admitancia: Y
BUS
o la matriz impedancia Z
BUS
. A
continuaci on se presenta la forma de obtener la matriz Y
BUS
a partir de las relaciones circuitales
que se pueden plantear en el sistema de potencia.
1.1. Determinaci on de las relaciones entre corrientes y ten-
siones nodales en un sistema interconectado
Todas las tensiones nodales del sistema pueden ser agrupadas en un vector siguiendo la
enumeraci on arbitraria que el analista haya hecho de los nodos. Las corrientes se dividen en dos
grandes grupos. El primer grupo lo conforman las corrientes asociadas a los generadores y a las
cargas del sistema. Estas corrientes se denominan corrientes inyectadas y, por convenci on,
siempre se toman entrando al nodo donde se encuentra conectado el generador o la carga. Estas
corrientes presentan una particularidad: provienen de la parte externa del sistema que representa
la parte no est atica o variable. El segundo grupo de corrientes lo conforman las corrientes que
uyen por el interior del sistema de potencia, es decir, a trav es de las lneas, de los transfor-
madores, y de los elementos pasivos. Estas corrientes uyen por la denominada parte interna
del sistema o simplemente por el sistema interno. La gura 1 muestra un sistema de potencia
con 4 nodos, 3 cargas y dos generadores, en donde el sistema se ha separado en la parte interna
y la parte externa.
En resumen, el sistema de potencia se divide en dos sistemas: el sistema externo y el sistema
interno. Al primer grupo pertencen los generadores y las cargas del sistema. El sistema interno lo
conforman las lneas, transformadores y dem as elementos pasivos que permanecen inalterados
durante la operaci on del sistema. Los nodos permiten la interconexi on entre el sistema interno
y el sistema externo.
sistema externo
sistema interno
Figura 1: Representaci on de un sistema de potencia a trav es del sistema interno y externo.
Cuando se formula matem aticamente el sistema solo se escriben las ecuaciones que relacionan
las corrientes inyectadas del sistema, las tensiones nodales y los par ametros de los elementos
del sistema interno escritos en forma de admitancias. Las corrientes que uyen por el interior
del sistema no se plantean explcitamente. Las relaciones existentes entre las tensiones nodales,
las corrientes netas inyectadas en cada nodo y las admitancias de los elementos del sistema, se
determinan aplicando a la red la primera y la segunda ley de Kirchhoff.
La gura 2 muestra el sistema de potencia de 4 nodos de la gura 1 modicado. En la nueva
representaci on, los generadores y las cargas, que conforman el sistema externo, se han reem-
plazado por las corrientes inyectadas I
1
, I
2
, I
3
e I
4
. Como se puede observar, la corriente I
1
representa el efecto neto del generador conectado al nodo 1 y la carga conectada al mismo no-
do. En el nodo 2, I
2
es la corriente inyectada por el generador conectado al nodo 2. En los nodos
3 y 4, las corrientes inyectadas I
3
e I
4
representan corrientes de las cargas conectadas a dichos
nodos.
La gura 2 tambi en presenta la estructura que asumir a el sistema de potencia cuando se repre-
sente matem aticamente. Las corrientes inyectadas representan el sistema externo, la matriz
Y
BUS
representa al sistema interno y las tensiones nodales V
1
, V
2
, V
3
y V
4
, asociadas a los
nodos, permiten unir estos dos sistemas.
sistema interno:
[Y]BUS
I
4
I
I
I
1
2
3
1
2
3
V
4
V
V
V
Figura 2: Representaci on de un sistema de potencia por corrientes inyectadas, tensiones
nodales y [Y ]
bus
.
Para llevar el sistema de potencia a una representaci on matem atica, se plantea inicialmente la
ley de Kirchhoff de corrientes a todos los nodos del sistema. Para esto es necesario escribir los
par ametros del sistema interno en forma de admitancias. Para un elemento con impedancia z
ij
conectado entre los nodos i y j, la admitancia se calcula como y
ij
= 1/z
ij
. La gura 3 muestra
las lneas de transmisi on del sistema de la gura 1 representadas a trav es de sus admitancias. A
partir de esta representaci on se pueden plantear las ecuaciones de corriente en los nodos.
Suma de corrientes que salen del nodo 1 = suma de corrientes que entran al nodo 1:
y
12
(V
1
V
2
) + y
13
(V
1
V
3
) = I
1
(1)
Suma de corrientes que salen del nodo 2 = suma de corrientes que entran al nodo 2:
y
12
(V
2
V
1
) + y
24
(V
2
V
4
) = I
2
(2)
Suma de corrientes que salen del nodo 3 = suma de corrientes que entran al nodo 3:
y
13
(V
3
V
1
) + y
34
(V
3
V
4
) = I
3
(3)
Suma de corrientes que salen del nodo 4 = suma de corrientes que entran al nodo 4:
y
34
(V
4
V
3
) + y
24
(V
4
V
2
) = I
4
(4)
I
4
I
I
I
1
2
3
1
2
3
V
4
V
V
V
34
y
y
y
y
12
13
24
Figura 3: Representaci on de par ametros de las lneas usando admitancias.
En las expresiones anteriores, las corrientes que circulan por el sistema interno han sido expre-
sadas en funci on de las tensiones nodales y las admitancias de los elementos usando la relaci on:
Corriente en la lnea ij: I
ij
= y
ij
(V
i
V
j
)
Esto quiere decir que, para el ejemplo anterior, las corrientes del sistema interno no aparecer an
explicitamente escritas en las ecuaciones y en su lugar aparecen las tensiones nodales y los
par ametros de las lneas.
Las expresiones (1), (2), (3) y (4) pueden reescribirse de la siguiente manera:
Suma de corrientes que salen del nodo 1 = suma de corrientes que entran al nodo 1:
(y
12
+ y
13
)V
1
y
12
V
2
y
13
V
3
= I
1
(5)
Suma de corrientes que salen del nodo 2 = suma de corrientes que entran al nodo 2:
y
12
V
1
+ (y
12
+ y
24
)V
2
y
24
V
4
= I
2
(6)
Suma de corrientes que salen del nodo 3 = suma de corrientes que entran al nodo 3:
y
13
V
1
+ (y
13
+ y
34
)V
3
y
34
V
4
= I
3
(7)
Suma de corrientes que salen del nodo 4 = suma de corrientes que entran al nodo 4:
y
24
V
2
y
34
V
3
+ (y
24
+ y
34
)V
4
= I
4
(8)
Reescribiendo las expresiones (5), (6), (7) y (8) en forma de un producto de dos vectores se
tiene:
I
1
= [ (y
12
+ y
13
) y
12
y
13
0 ]

V
1
V
2
V
3
V
4

I
2
= [ y
12
(y
12
+ y
24
) 0 y
24
]

V
1
V
2
V
3
V
4

I
3
= [ y
13
0 (y
13
+ y
34
) y
34
]

V
1
V
2
V
3
V
4

I
4
= [ 0 y
24
y
34
(y
24
+ y
34
) ]

V
1
V
2
V
3
V
4

Estas ecuaciones pueden reescribirse de manera compacta en forma matricial, de la siguiente

forma:

I
1
I
2
I
3
I
4

(y
12
+ y
13
) y
12
y
13
0
y
12
(y
12
+ y
24
) 0 y
24
y
13
0 (y
13
+ y
34
) y
34
0 y
24
y
34
(y
24
+ y
34
)

V
1
V
2
V
3
V
4

(9)
En esta representaci on matricial, el vector de corrientes se denomina vector de corrientes in-
yectadas, el vector de tensiones se denomina vector de tensiones nodales o variables de estado,
y la matriz conformada por las admitancias de las lneas (representadas con letras min usculas)
se denomina matriz admitancia o Y
BUS
. En forma general, los elementos de la matriz Y
BUS
se
denominan con letras may usculas y cada elemento toma los ndices de la la y la columna a
la que pertenecen. En consecuencia, el elemento Y
ij
representa el elemento de la matriz Y
BUS
situado en la la i y en la columna j. Para el ejemplo anterior se tiene:

I
1
I
2
I
3
I
4

Y
11
Y
12
Y
13
Y
14
Y
21
Y
22
Y
23
Y
24
Y
31
Y
32
Y
33
Y
34
Y
41
Y
42
Y
43
Y
44

V
1
V
2
V
3
V
4

(10)
Al comparar las expresiones (9) y (10) se concluye que:

Y
11
Y
12
Y
13
Y
14
Y
21
Y
22
Y
23
Y
24
Y
31
Y
32
Y
33
Y
34
Y
41
Y
42
Y
43
Y
44

(y
12
+ y
13
) y
12
y
13
0
y
12
(y
12
+ y
24
) 0 y
24
y
13
0 (y
13
+ y
34
) y
34
0 y
24
y
34
(y
24
+ y
34
)

Se puede concluir que: 1) Los elementos de la diagonal de la matriz Y

BUS
son positivos y se
obtienen sumando las admitancias de los elementos que llegan a cada nodo; 2) Los elementos
por fuera de la diagonal son negativos y corresponden a las admitancias de los elementos que
unen los nodos. Si no existe un elemento uniendo dos nodos entonces all hay una impedancia
innita (circuito abierto) o una admitancia cero. Escrito en forma matem atica se tiene que:
Y
ii
=
n

j=1
y
ij
Y
ij
= y
ij
En general, para un sistema de n nodos se tiene:

I
1
.
.
.
I
i
.
.
.
I
n

Y
11
... Y
1i
... Y
1n
.
.
.
.
.
.
.
.
.
.
.
.
Y
i1
... Y
ii
... Y
in
.
.
.
.
.
.
.
.
.
.
.
.
Y
n1
... Y
ni
... Y
nn

V
1
.
.
.
V
i
.
.
.
V
n

Cuya notaci on es la siguiente:

[I] = [Y
BUS
][V ]
En el ejemplo anterior vamos a suponer que los valores de las impedancias de las lneas en p.u.
corresponden a los siguientes valores:
z
12
= (0,12 + j0,41)p.u.; z
13
= (0,15 + j0,54)p.u.; z
24
= (0,13 + j0,46)p.u. y z
34
= (0,11 +
j0,39)p.u. El sistema de potencia puede representarse de la forma mostrada en la gura 4.
I
4
I
I
I
1
2
3
1
2
3
V
4
V
V
V
z 12 = (0.12 + j 0.41) p.u.
z

2
4
=

(
0
.
1
3

+

j

0
.
4
6
)

p
.
u
.

Figura 4: Sistema de potencia en funci on de impedancias.
Al calcular la inversa de las impedancias se obtienen las siguientes admitancias:
y
12
=
1
z
12
=
1
(0,12+j0,41)p.u.
= (0,6575 j2,2466)p.u.
y
13
=
1
z
13
=
1
(0,15+j0,54)p.u.
= (0,4776 j1,7192)p.u.
y
24
=
1
z
24
=
1
(0,13+j0,46)p.u.
= (0,5689 j2,0131)p.u.
y
34
=
1
z
34
=
1
(0,11+j0,39)p.u.
= (0,6699 j2,3752)p.u.
El sistema en funci on de admitancias corresponde al de la gura 5.
En consecuencia, los elementos de la matriz Y
BUS
son los siguientes:
Y
11
es la suma de las admitancias que llegan al nodo 1, por lo tanto:
Y
11
= y
12
+ y
13
= (0,6575 j2,2466)p.u. + (0,4776 j1,7192)p.u.
Y
11
= (1,1351 j3,9658)p.u.
I
4
I
I
I
1
2
3

y 12 = (0.6575 - j 2.2466)p.u.
y
2
4
=

(
0
.
5
6
8
9

-
j

2
.
0
1
3
1
)
p
.
u
.

Figura 5: Sistema de potencia en funci on de admitancias.
Y
12
es la admitancia del elemento que une al nodo 1 y el nodo 2, con signo negativo.
Y
12
= y
12
= (0,6575 j2,2466)p.u. = (0,6575 + j2,2466)p.u.
Y
13
es la admitancia del elemento que une al nodo 1 y el nodo 3, con signo negativo.
Y
13
= y
13
= (0,4776 j1,7192)p.u. = (0,4776 + j1,7192)p.u.
Y
14
es la admitancia del elemento que une al nodo 1 y el nodo 4, con signo negativo.
Y
14
= y
14
= 0 p.u. no existe una lnea que una el nodo 1 y el nodo 4.
Y
21
es la admitancia del elemento que une al nodo 2 y el nodo 1, con signo negativo.
Y
21
= y
21
= (0,6575 + j2,2466)p.u. y
21
es la misma admitancia y
12
Y
22
es la suma de las admitancias que llegan al nodo 2.
Y
22
= y
12
+ y
24
= (0,6575 j2,2466)p.u. + (0,5689 j2,0131)p.u.
Y
22
= (1,2264 4,2597)p.u.
Y
23
es la admitancia del elemento que une al nodo 2 y el nodo 3, con signo negativo.
Y
23
= 0 p.u. no existe una lnea que una el nodo 2 y el nodo 3.
Y
24
es la admitancia del elemento que une al nodo 2 y el nodo 4, con signo negativo.
Y
24
= y
24
= (0,5689 + j2,0131)p.u.
Y
31
es la admitancia del elemento que une al nodo 3 y el nodo 1, con signo negativo.
Y
31
= y
13
= (0,4776 + j1,7192)p.u. y
31
es la misma admitancia y
13
Y
32
es la admitancia del elemento que une al nodo 3 y el nodo 2, con signo negativo.
Y
32
= y
32
= 0 p.u. no existe una lnea que una el nodo 3 y el nodo 2.
Y
33
es la suma de las admitancias que llegan al nodo 3.
Y
33
= y
13
+ y
34
= (0,4776 j1,7192)p.u. + (0,6699 j2,3752)p.u.
Y
33
= (1,1475 j4,0944)p.u.
Y
34
es la admitancia del elemento que une al nodo 3 y el nodo 4, con signo negativo.
Y
34
= y
34
= (0,6699 + j2,3752)p.u.
Y
41
es la admitancia del elemento que une al nodo 4 y el nodo 1, con signo negativo.
Y
41
= y
41
= 0 p.u. no existe una lnea que una el nodo 4 y el nodo 1.
Y
42
es la admitancia del elemento que une al nodo 4 y el nodo 2, con signo negativo.
Y
42
= y
42
= (0,5689 + j2,0131)p.u.
Y
43
es la admitancia del elemento que une al nodo 4 y el nodo 3, con signo negativo.
Y
43
= y
43
= (0,6699 + j2,3752)p.u.
Y
44
es la suma de las admitancias que llegan al nodo 4.
Y
44
= y
24
+ y
34
= (0,5689 j2,0131)p.u. + (0,6699 j2,3752)p.u.
Y
44
= (1,2388 j4,3883)p.u.
La matriz Y
BUS
asume entonces la siguiente forma:

(1,1351 j3,9658) (0,6575 + j2,2466) (0,4776 + j1,7192) 0

(0,6575 + j2,2466) (1,2264 4,2597) 0 (0,5689 + j2,0131)
(0,4776 + j1,7192) 0 (1,1475 j4,0944) (0,6699 + j2,3752)
0 (0,5689 + j2,0131) (0,6699 + j2,3752) (1,2388 j4,3883)

La matriz Y
BUS
en formato de magnitud y angulo asume la siguiente forma (todos los valores
se encuentran en p.u.):

4,125074,0277 2,3408106,3139 1,7843105,5241 00

2,3408106,3139 4,4328 73,9377 00 2,0920105,7808
1,7843105,5241 00 4,252174,3441 2,4678105,7512
00 2,0920105,7808 2,4678105,7512 4,559874,2353

1.2. Caractersticas especiales de la matriz Y

bus
Es una matriz cuadrada de tama no (n n) siendo n el n umero de nodos.
Es una matriz sim etrica. esto quiere decir que el elemento Y
ij
es igual al elemento Y
ji
.
Los elementos de la diagonal son positivos y se pueden calcular sumando las admitancias
de los elementos conectados a los nodos. Por ejemplo, el elemento Y
ii
es la suma de las
admitancias de los elementos conectados al nodo i. Por lo tanto:
Y
ii
=
n

j=1
y
ij
Los elementos por fuera de la diagonal Y
ij
de la matriz Y
BUS
son el negativo de la admi-
tancia del elemento y
ij
. Por lo tanto:
Y
ij
= y
ij
La matriz es diagonalmente dominante, es decir, la magnitud de los elementos de la dia-
gonal es num ericamente mayor que la magnitud de los elementos fuera de la diagonal.
Esto puede observarse en el ejemplo resuelto.
En sistemas de potencia donde la reactancia supera en tres o m as veces el valor de la
resistencia, el angulo de los elementos de la diagonal resultan negativos y de magnitud
menor a 90 grados. Los angulos de los elementos fuera de la diagonal resultan positivos
y superiores pero cercanos a 90 grados.
En sistemas reales de gran tama no la matriz Y
Bus
es dispersa, es decir, contiene muchos
ceros. Esta condici on ocurre porque normalmente s olo existe conexi on directa entre nodos
geogr acamente vecinos y no existe conexi on directa entre nodos geogr acamente dis-
tantes. En un sistema el ectrico de la vida real del orden de 100 nodos, un nodo promedio
puede estar interconectado con 3 o 4 nodos, y unos pocos nodos pueden tener conexi on
con 8 o 9 nodos. La gura 6 muestra el disgrama unilar del sistema el ectrico colombiano
2012 reducido a 93 nodos y 155 corredores y que es utilizado como sistema de prueba en
investigaciones de planeamiento de sistemas de transmisi on de energa el ectrica. Puede
observarse que el nodo 27 tiene conexi on con 8 nodos, mientras los nodos 51 y 52 solo
presentan conexi on con 2 nodos. Las lneas punteadas son conexiones futuras mientras las
lneas llenas representan corredores existentes. Si se observa el esquema se puede notar
que en general los nodos se conectan unicamente con algunos (pocos) nodos vecinos.
50
81
45
54
48 63
47
53
46
56
51 52
43
42
41
40
38
57
84
49
39
88
68 86
37
61
12
66
22
76
82
60
62
16
20
23
21
15
55
75
34 70
69
13
14
74
64
30
73
72
27
89
32
65
36
35
19
2
9
4
5
77 79
93
67
87
78
6
10
1
92
11
85
83
25
28
29
80
26
44
91
90
3
7
71
59
33
31
24
8
18
17
58
Figura 6: Sistema Colombiano de 230KV y 500KV, de 93 nodos y 155 corredores.
Se deja al lector la tarea de determinar el n umero de conexiones que tiene cada uno de los 93
nodos de este sistema y sacar una conclusi on para el sistema el ectrico colombiano 2012.
1.3. Forma alternativa de construir la matriz Y
BUS
.
La matriz Y
BUS
tambi en puede construirse a trav es de la multiplicaci on de tres matrices, que en
principio son m as sencillas de construir. Seg un esta forma alternativa, la matriz Y
BUS
se puede
obtener a trav es del siguiente producto:
Y
BUS
= A Y
primitiva
A
T
Donde:
A es la matriz incidencia nodo-rama o tambien denominada matriz incidencia nodo-
elemento. Es una matriz cuyo n umero de las corresponde al n umero de nodos del sistema
y el n umero de columnas al n umero de elementos que contiene el sistema interno.
A
T
es la transpuesta de la matriz A, es decir, una matriz donde las las de A se convierten
en las columnas de A
T
.
Y
primitiva
es una matriz que solo contiene elementos en la diagonal, si no existen acoples
mutuos entre elementos. La cantidad de t erminos en la diagonal de esta matriz es igual
a n umero de elementos que contenga el sistema interno y cada posici on en la diagonal
asume el valor de la admitancia de cada elemento.
Para el sistema de potencia de la gura 5, vamos a asumir que la lnea de transmisi on que une
los nodos 3 y 4 sale de servicio por razones de mantenimiento. El sistema asume entonces la
forma mostrada en la gura 7. Se desea encontrar la matriz Y
BUS
del sistema de potencia para
la nueva conguraci on.
Para el sistema de la gura 7 se dene las siguiente matriz primitiva en p.u.:
1-2 1-3 2-4
1 2
1 3
2 4

0,6575 j2,2466 0 0
0 0,4776 j1,7192 0
0 0 0,5689 j2,0131

Como se puede observar, los elementos del sistema, que en este caso son lneas de transmisi on
son nombrados de acuerdo a los nodos entre los que se encuentran conectados. Por lo tanto,
el elemento 1-2 representa la lnea que une los nodos 1 y 2. El elemento 1-3 es la lnea que
une los nodos 1 y 3; y el elemento 2-4 es la lnea que une los nodos 2 y 4. Observese tambi en
Figura 7: Sistema Colombiano de 230KV y 500KV, de 93 nodos y 155 corredores.
que el elemento de la diagonal donde se intersecta la la 1-3 con la columna 1-3, contiene la
admitancia de la lnea de transmisi on que une los nodos 1 y 3, y los dem as elementos de la la
1-3 y de la columna 1-3 son ceros.
A continuaci on construimos la matriz incidencia nodo-rama que denominamos A. Para esto
construimos una matriz con tantas las como nodos tenga el sistema, nombrados en orden
ascendente, y tantas columnas como elementos tenga el sistema. Dado que el sistema de la
gura 7 tiene 4 nodos (1, 2, 3 y 4) y tres lneas de transmisi on (1-2, 1-3 y 2-4), esta matriz resulta
de cuatro las y tres columnas (4x3). A continuaci on se muestra la matriz A correspondiente a
este sistema de potencia. Debe tenerse cuidado de colocar los elementos en las columnas de A
en el mismo orden en que fueron colocadas en la matriz Y
primitiva
.
1-2 1-3 2-4
A =
1
2
3
4

1 1 0
1 0 1
0 1 0
0 0 1

Cada columna de la matriz incidencia nodo-rama o nodo-elemento s olo debe contener una posici on con
valor 1 y una posici on con valor -1. Las dem as posiciones deben ser ceros. Cada columna est a asociada
a un elemento del sistema as: la columna 1 est a asociada al elemento 1-2, la columna 2 al elemento 1-3
y la columna 3 al elemento 2-4. Cada columna debe tener un 1 en la la correspondiente al nodo inicial
del elemento y un -1 en la la correspondiente al nodo nal del elemento. Por lo tanto, en la columna 1,
que corresponde al elemento 1-2, debe haber un 1 en la la 1 y un -1 en la la 2. En la columna 2, que
corresponde al elemento 1-3, debe haber un 1 en la la 1 y un -1 en la la 3. Finalmente, en la columna
3, que corresponde al elemento 2-4, debe haber un 1 en la la 2 y un -1 en la la 4.
Ahora realizamos la operaci on: Y
BUS
= A Y
primitiva
A
T
.

1 1 0
1 0 1
0 1 0
0 0 1

0,6575 j2,2466 0 0
0 0,4776 j1,7192 0
0 0 0,5689 j2,0131

1 1 0 0
1 0 1 0
0 1 0 1

Que resulta igual a:

1 1 0
1 0 1
0 1 0
0 0 1

0,657 j2,246 0,657 + j2,246 0 0

0,477 j1,719 0 0,477 + j1,719 0
0 0,568 j2,013 0 0,568 + j2,013

Finalmente se llega a:

1,1351 j3,9658 0,6575 + j2,2466 0,4776 + j1,7192 0

0,6575 + j2,2466 1,2264 4,2597 0 0,5689 + j2,0131
0,4776 + j1,7192 0 0,4776 j1,7192 0
0 0,5689 + j2,0131 0 0,568 j2,013

Que corresponde a la matriz Y

BUS
del sistema de potencia analizado.