Sistema 3 Barras Marti

Sea el Presente Sistema

FiguraN1: Diagrama unifilar del Sistema Jose Marti para el estudio de Flujo de Potencia

Donde los Datos del Circuitos se presentan a continuacin:

A continuacin se llevaran a cabo la serie de tareas que corresponden al estudio de

Flujo de Potencia del Sistema compuesto por tres barras del Jos Marti

1.- Construccin de manera manual de la matriz admitancia Ybus del Sistema de

potencia, expresada en el Sistema por Unidad .

Primer Paso: Determinar las Bases del Sistemas.

Vbase=230KV

Sbase=100MVA

Zbase=(230kV)2/(100MVA)=529

Segundo paso: Determinar las Impedancias del sistema.

A) Lneas de Transmisin.

Lnea de transmisin de conexin Barra 1-2:

Datos conocidos de la lnea de transmisin
Longitud.=120Km; R(/Km)=0.08 ;X( /Km)=0.477 ;b(nF/Km)=9.053(nF/Km)
Entonces se procede a determinar la impedancia de la lnea:

Zlinea = 120 Km * 0.08

+ 0.477
= 9.6 + 57.24 j = (58.039480.4793)
Km
Km

Se sabe que:
b(nF/Km)=9.053(nF/Km) , por lo que

b(nF)=9.053(nF/Km)*120Km=1.0863F

La reactancia Capacitiva de la lnea ser:

XcLT =

1
= 2.44172 K
2 * 60 Kz *1.08636 F

Linea de transmisin de conexin Barra 2-3:

Datos conocidos de la linea de transmisin
long.=180Km

R(/Km)=0.08 ;

X( /Km)=0.477 ; b(nF/Km)=9.053(nF/Km)

Entonces se procede a determinar la impedancia de la linea:

+ 0.477
Zlinea = 180 Km * 0.08
= 14.4 + 85.86 j = (87.059280.4793)
Km
Km

Se sabe que la suceptancia capacitva de la Lnea es:

b(nF/Km)=9.053(nF/Km) , entonces por la longitud de la linea es:

b(nF)=9.053(nF/Km)*180Km=1.62954F
Por lo que la reactancia capacitiva de la lnea de transmisin es:
XcLT =

1
= 1.62781K
2 * 60 Kz *1.62954 F

Lnea de transmisin de conexin Barra 1-3:

Datos conocidos de la lnea de transmisin
longitud=80Km

R(/Km)=0.05 X( /Km)=0.477 b(nF/Km)=9.053(nF/Km)

Entonces se procede a determinar la impedancia de la lnea:

Zlinea = 80 Km * 0.05
+ 0.477
= 6.4 + 38.16 j = (38.69380.4793)
Km
Km

Se sabe que la suceptancia capacitva de la Lnea es:

b(nF/Km)=9.053(nF/Km), entonces por la longitud de la lnea es:
b(nF)=9.053(nF/Km)*80Km=724.24nF
Por lo que la reactancia capacitiva de la lnea de transmisin es:
XcLT =

1
= 3.66257 K
2 * 60 Kz * 724.24nF

Tercer paso: Determinacin de las impedancias del sistema en por unidad.

Lnea de transmisin de conexin Barra 1-2:

9.6 + 57.24 j
3
3
Z LT 1_ 2 =
= 18.1474 x10 + 108.204 x10 j
529

2.44172 K
XjLT 1_ 2 =
= 4.61572 P.U .
529

1
XjLT 1_ 2

Linea de transmisin de conexin Barra 2-3:

1
= 216.651x10 3
4.61572 P.U .

 14.4 + 85.86 j
3
3
Z LT 2 _ 3 =
= 27.2212 x10 + 162.306 x10 j
529

1.62781K
Xj2 _ 3 =
= 3.07715 P.U .
529

1
XjLT 2 _ 3

1
= 324.976 x10 3
3.07715 P.U .

Lnea de transmisin de conexin Barra 1-3:

6.4 + 38.16 j
3
3
Z LT 1_ 3 =
= 12.0983 x10 + 72.1361x10 j
529
3.66257 K
Xj1_ 3 =
= 6.92358 P.U .
529

1
XjLT 1_ 3

1
= 144.434 x10 3 .
6.92358 P.U .

Cuarto Paso: Diagrama de reactancias del sistema de transmisin.

A continuacin se onstuye el Diagrama contemplando las reactancias del sistema de

transmisin:
1/2(Xc12)
1/2(Xc13)

Barra1

1/2(Xc13)

Z13

1/2(Xc12)

Z12

Barra3

1/2(Xc23)

Barra2

Z23

1/2(Xc23)

FiguraN2: Diagrama de reactancias del sistema de transmisin

Quinto Paso: Calculo de los elementos de la matriz admitancia de Barra.

Por ser un sistema que contempla tres barras la Matriz admitancia ser de la forma:

Y 11 Y 12 Y 13
[Y ] = Y 21 Y 22 Y 23
Y 31 Y 32 Y 33
Donde cada elemento contempla el inverso de la impedancia de los tramos de las lneas
de transmisin que se muestran en la Figura N2 .Se procede entonces a determinar los
elementos que componen la matriz:

Elementos de la diagonal principal:

YC 12 YC 13
+
=
2
2
(1.50758 8.98895 j ) + (2.26137 13.4834 j ) + (108.325 x10 3 j ) + (72.217 x10 3 j ) = (3.76894983 22.291820973 j )

Y 11 = Y 12 + Y 13 +

YC 12 YC 23
+
=
2
2
(1.50758 8.98895 j ) + (1.00505 5.99263 j ) + (108.325 x10 3 ) + (162.488 x10 3 j ) =
= 2.49503 14,9815 j

Y 22 = Y 12 + Y 23 +

YC 13 YC 23
=
+
2
2
(2.26137 13.4834 j ) + (1.00505 5.99263 j ) + (72.217 x10 3 ) + (162.488 x10 3 j ) =
= 3.2664 19.24134 j

Y 33 = Y 13 + Y 23 +

Elementos Fuera de la Diagonal Principal.

Y 12 = Y 12 = (1.507579935222 8.9889453637642 j )
Y 13 = (2.2613699028335 13.483418045647 j )
Y 23 = (1.0050532901483 5.9926302425096 j )

Entonces la Matriz admittancia de barra sera:

(3.76894983 22.291820973 j )
(1.507579935222 8.9889453637642 j ) (2.2613699028335 13.483418045647 j )
(1.507579935222 8.9889453637642 j )
(2.495016 14.98157 j )
(1.0050532901483 5.9926302425096 j ) = Ybus
(2.2613699028335 13.483418045647 j ) (1.0050532901483 5.9926302425096 j )
(3.2662231 19.24134 j )

2.- Implementacin del Programa PowerWorld y contraste del resultado del

matriz de admitancia Ybus obtenido mediante el software de programacin

FiguraN3: Interfaz del software de simulacin PowerWorld mostrando la resolucin del sistema de
potencia en estudio mostrado en la Figura N1

FiguraN4: Interfaz del software de simulacin PowerWorld mostrando el resultado de la matriz

admitancia del sistema que se muestra en la figura N3

En la Figura N4 se muestra el resultado de la matriz admitancia de barra del

sistema de potencia en estudio mediante la simulacin en el software PowerWorld.
Cada elemento que compone la matriz admitancia en esta solucin son bastante

aproximado al que es presentado mediante la aplicacin de los mtodos de algebra

compleja y manejo de favores estudiados durante el curso de Sistemas de Potencia
II. Cabe destacar, y como se muestra en la figura N4, la precisin en cada resultado
de cada elemento que compone la matriz admitancia e de solo Dos decimales, por lo
que una comparacin absoluta no pudiera realizarse.

3.- Planteamiento de las ecuaciones de balance de Potencia en cada barra del

Sistema de Potencia del Jos Marti en estudio.

Primer Paso: Se establecen los tipos de barras existentes en el sistema:

Tabla N1: Establecimiento de los tipo de barra existentes en el sistema en el estudio

Bus
1
2
3

X
X

P
X

Q
X
X

Type
Slack
PV
Load

Segundo paso: se establecen los valores en por unidad de las variables

establecidas en cada barra del sistema

Elementos en p.u de la Barra 2

.
- P2 =

(Pg 2 Pcarga2) = (200MW 450MW ) = 2.5 pu

- Q2 =

(Qg 2 Ql 2) = (Qg calc 200MW )

Sbase

Sbase

100 MVA

100 MVA

en donde se tiene que :

Qg calc = IImV 2 * [Link] = IIm V 2 * [Y 21xV 1 + Y 22 xV 2 + Y 23 xV 3]

i =1

-V2 =

234.6 Kv
= 1.02
[Link]

Elementos en p.u de la Barra 3

- P3 =

(Pg 3 Pload ) = (250MW ) = 2.5 pu

- Q3 =

(Qg 3 Qload 3) = ( 100MVA) = 1 pu

Sbase

Sbase

100 MVA

100 MVA

Tercer paso: Se plantea las ecuaciones de balance de potencia las cuales son de
la forma:
n

Pk jQk = V * [Link]
i =1

Ecuacin de Balance de potencia para la Barra 2

n

P 2 jQ 2 = V 2* Y [Link] , entonces sustituyendo los

i =1

elementos conocidos se tiene que:

Esta era de la forma

( 2.5) Qg calc 200MW = conj (V 2)[(1,.5075 + 8.98892 j )V 1 + (2.49503 14.98157 j )V 2 + (1.00505 + 5.9926 j )V 3]
100 MVA

En donde se sabe que:

Qg calc = Im V 2 * [(1.50758 8.98895 j ) xV 1 + (2.49503 14.98157 j ) xV 2 + (1.00505 5.99263 j ) xV 3]

Ecuacin de Balance de Potencia para la Barra 3

n

P3 jQ3 = V 3* Y [Link] , entonces sustituyendo los

i =1

elementos conocidos se tiene que:

Esta sera de la forma

( 2.5 + 1.0) = conj (V 3)[(2.26137 + 13.4834 j ) xV 1 + (1.00505 + 5.99263 j ) xV 2 + (3.2264 19.24134 j ) xV 3]

4.- Resolucin de manera Manual de las 3 primeras iteraciones del flujo de

potencia del sistema, Empleando el mtodo de Gauss Seidel y comparacin de los
resultados con los obtenidos en el software PowerWorld.
Se procedera entonces a la determinacin del flujo de potencia del sistema
estudiado:

Primer Paso: Determinacin de la condiciones iniciales en el proceso de

iteracion.

Una vez clasificadas las barras del sistema de potencia, mostradas en la Tabla N1,
entonces se establecen sus respectivas condiciones iniciales:

Barra PV 20 = 0
Barra PQ V 3 = 1, 30 = 0

Segundo Paso: Planteamiento de las Ecuaciones iterativas

Ecuacin iterativa de la Barra 2

V2( K +1) =

V2( K +1) =

1 P2esp jQ2K +1 21
( K +1)

2
.
m
Vm
Y [Link] ( K )
Y

*k
V2
Y 22
m =1
m = 2 +1

2.5 jQ 2
1
(1.50757 + 8.9889 j )V1( K +1) (1.00505 + 5.9926 j )V3K

(* K )
2.49503 14.9815 j V2

Donde:

)]

( K +1)
Q 2 = Qcalc
2 pu = Im V2(*K ) Y 21 * Vm ( K +1) + Y 22 * V2( K ) + Y 23 * V3( K ) 2 pu

)]

Q 2 = Im V2(*K ) (1.50757 + 8.9889 j ) * Vm( K +1) + (2.49503 14.9815 j ) * V2( K ) + (1.00505 + 5.9926 j ) * V3( K ) 2 pu

Ecuacin iterativa de la Barra 3

( K +1)
3

n
1 P3 jQ3 31
( K +1)
(K )
=

Y
3
m
.
Vm
Y
3
m
.
Vm

Y 33 V3*k
m =1
m =3+1

Tercer Paso: proceso numrico iterativo para la determinacin de las

condiciones del sistema.

Debido a que el generador conectado a la barra 2 del sistema estudiado, posee limites de
generacin de potencia reactiva los cuales son Qmin=-100 MVAr y Qmax=230MVAr
se verificara que la potencia reactiva en cada iteracin se mantenga entre estos limites
de potencia ya si asegurar que el voltaje en la barra se mantenga igual para cada
interaccin conservando la caracterstica de barra PV.
Los Limites de potencia reactiva en por unidad seran:
Q min p.u =

100 MVAr
= 1 p.u
100 MVAr

Q max p.u =

250100 MVAr
= 2.3 p.u
100 MVAr

Iteracin 1 para K=0

Q 21 = IIM [(1.020)((1.50757 + 8.9889 j ) * (1.020) + (2.49503 14.9815 j ) * (1.020) + (1.00505 + 5.9926 j ) * (10) )] 2 pu
Q 21 = 2.06112 pu

Verificacin de Generacion de Potencia Reactiva

La potencia reactiva generada es: Qgen = Qcalc = 61.1248 x10 3 jp.u
Analizando el intervalo se tiene que Q min < Qgen < Q max , por lo que la barra
mantiene la clasificacin de ser una barra PV.
Se sigue con el proceso iterativo:
V2(1) =

2.5 1 j

1
(1.50757 + 8.9889 j )(1.040) (1.00505 + 5.9926 j )(1.020) = 1.02 8.9307
2.49503 14.9815 j (1.00)

V3(1) =

2.5 1 j

1
(1.50757 + 8.988 j )(1.040) + (1.00505 + 5.9926 J )(1.02 8.9307)
3.264 19.24134 j (1.00)

V3(1) = (0.984825 9.88616)

Iteracin 2 para K=1

(1.50757 + 8.9889 j ) * (1.040) +

2 pu
Q 2 2 = IIM (1.02 8.9307)*
(
2
.
49503
14
.
9815
j
)
*
(
1
.
02
8
.
9307
)
(
1
.
00505
5
.
9926
j
)
*
(
0
.
984825
9
.
88615
)

Q 2 2 = 1.6205 pu
Verificacin de Generacion de Potencia Reactiva
La potencia reactiva generada es:

Qgen = 0.379501 p.u

Analizando el intervalo se tiene que Q min < Qgen < Q max , por lo que la barra
mantiene la clasificacin de ser una barra PV.
Se sigue con el proceso iterativo:

( 2)
2

2.5 1 j
(1.50757 + 8.9889 j )(1.040) +
1

(1.08.9307)
= 1.02 12.214
=

2.49503 14.9815 j
+ (1.00505 + 5.9926 j )(0.984825 9.88616)

( 2)
3

2.5 1 j

(1.50757 + 8.988 j )(1.040) +

1

(1.00)
=

3.264 19.24134 j
(1.00505 + 5.9926 J )(1.02 12.2413)

V3(1) = 0.962229 10.4461

Iteracin 3 para K=2

(1.50757 + 8.9889 j ) * (1.040) +

2 p
Q 23 = IIM (1.02 12.2143)*
(
2
.
49503
14
.
9815
)
*
(
1
.
02
12
.
2143
)
(
1
.
00505
5
.
9926
)
*
(
0
.
962229
10
.
4461
j
j
)

Q 23 = 1.2428 pu
Verificacin de Generacin de Potencia Reactiva
La potencia reactiva generada es: Qgen = 0.75719 p.u
Analizando el intervalo se tiene que Q min < Qgen < Q max , por lo que la barra
mantiene la clasificacin de ser una barra PV.
Se sigue con el proceso iterativo:

( 3)
2

2.5 1.24281 j
(1.50757 + 8.9889 j )(1.040) +
1

(1.02 12.2143)
=

2.49503 14.9815 j
+ (1.00505 + 5.9926 j )(0.962229 10.4461)

V2(3) = 1.02 12.2753

( 3)
3

2.5 1 j

(1.50757 + 8.988 j )(1.040) +

1
*

(0.962229 10.4461)
=

3.264 19.24134 j
(1.00505 + 5.9926 J )(1.02 12.2753)

V3(3) = 0.959334 10.5916

Comparacin luego de tres iteraciones mediante el mtodo numrico Gauss-Seidel

con los resultados obtenidos mediante el Software de simulacin PowerWorld
A continuacin se muestra la interfaz del simulador PowerWorld mostrando el resultado
del flujo de potencia del sistema en estudio:

FiguraN5: Interfaz del software de simulacin PowerWorld mostrando el resultado de la simulacin del
flujo de potencia del sistema mostrado e la figura N3.

Los Resultados mostrados por el programa difieren de los obtenidos despus de

tres iteraciones mediante el mtodo de Gauss Seidel. Esta diferencia es atribuida al
siguiente hecho: Existe la posibilidad que para un mayor numero de iteracin de
iteraciones, la potencia reactiva de generacin en la barra 2 exceda en alguna dimensin
a los limites establecidos de potencia reactiva, por lo que la barra perdera la condicin
de barra PV y el voltaje de esta barra variara.

5.- Procedimiento a resolver la primera iteracin del flujo de potencia del sistema
empleando el mtodo de Newton-Rapshon. Para la primera iteracin, comparar el
Jacobiano numrico obtenido con el que se logra con el programa Powerworld.

Las ecuaciones de potencia activa y reactiva , por barra ,en el sistema de potencia
estudiado tendrn la forma:

Pi =| Vi |2 .Yii cos ii + | [Link] | . cos(in + n i )

n =i
n i

Qi = | Vi |2 . | Yii | .sen(ii ) + | [Link] | .sen(in + n i )

n =i
n i

Derivando cada ecuacin de potencia activa y reactiva respecto a cada elemento que la
compone, se llega a la expresin de la matriz jacobiana del sistema. La exprsion
iterativa tendr la forma:

P 2
2
P 3
2
Q 2
2
Q 3
2

P 2
3
P 3
3
Q 2
3
Q 3
3

P 2
|V 2 |
|V 2 |
P 3
|V 2 |
|V 2 |
Q 2
|V 2 |
|V 2 |
Q 3
|V 2 |
|V 2 |

P 2
|V 3 |
|V 3 |
P 3
|V 3 |
|V 3 |
*
Q 2
|V 3 |
|V 3 |
Q 3
|V 3 |
|V 3 |

P 2
3
| V 2 | P 3
=
|V 2 |
Q 2
| V 3 | Q 3
|V 3 |

En donde la expresin de la matriz Jacobina es:

J11

J21

P 2
2
P 3
2
Q 2
2
Q 3
2

P 2
3
P 3
3
Q 2
3
Q 3
3

P 2
|V 2 |
|V 2 |
P 3
|V 2 |
|V 2 |
Q 2
|V 2 |
|V 2 |
Q 3
|V 2 |
|V 2 |

P 2
|V 3 |
|V 3 |
P 3
|V 3 |
|V 3 |
Q 2
|V 3 |
|V 3 |
Q 3
|V 3 |
|V 3 |

J12

J22

Se mantienen entonces la condiciones iniciales de las variables del sistema segn la

clasificacin de barras hechas en la TablaN1

Barra PV 20 = 0
Barra PQ V 3 = 1, 30 = 0

Primer Paso: Determinacin de los elementos de la Matriz Jacobiana para

condiion de itracion K=0.

Segn la subdivisin realizada a la matriz cada elemento es:

Elementos de la Jacobiana J11

Fuera de la Diagonal

P 2
= | V 2.V 3.Y 23 | .sen( 23 + 3 2) =
3

P 2
= (1.02)(1)(6.07633) sen(1.73696rad + 0rad 0rad ) = 6.11248
3

Elementos de la Diagonal

P 2
= | V 2.V 1.Y 21 | .sen( 21 + 1 2) =
2
P 2
= (1.02)(1.04)(9.11449) sen(1.73696rad + 0rad 0rad ) + (1.02)(1)(6.07633) sen(1.73696 + 0 + 0) = 15.648
2

P 3
= | V 3.V 1.Y 31 | .sen( 31 + 1 3)+ | V 3.V 2.Y 32 | .sen( 32 + 2 3) =
3
P 3
= ((1)(1.04)(13.4986)) sen(1.6828 + 0 0) + ((1)(1.04)(6.07633)) sen(1.73696 + 0 0) =
3
P 3
= 20.1352
3

Elementos del Jacobiano J21

Fuera de la Diagonal

Qi
= | [Link] | . cos(ij + j i )
j

Q 2
= | V 2.V 3.Y 23 | . cos( 23 + 3 2) = (1.02)(1)(6.07633) cos(1.73696rad + 0rad 0rad ) =
3

Q 2
= 1.02513
3

Q 3
= | V 3.V 2.Y 32 | . cos( 32 + 2 3) = (1.02)(1)(6.07633) cos(1.73696rad ) = 1.02513
2

En la Diagonal

Qi N
= | [Link] | . cos(in + n i )
j n =i
n i

Q 2 N
= | V [Link].Y 2n | . cos( 2n + n 2) =
2 n =i
n i

| V 2.V 1.Y 21 | . cos( 21 + 1 2)+ | V 2.V 3.Y 23 | . cos( 23 + 3 2) =

(1.02)(1.04)(9.1144) cos(1.73696rad ) + (1.02)(1)(6.07633) cos(1.73696rad ) = 2.62437

Q 3
= (V 3.V 1.Y 31). cos( 31 + 1 3) + (V 3.V 2.Y 32) cos( 23 + 2 3)
3
Q 3
= (1.00)(1.04)(13.6717) cos(1.73696rad ) + (1)(1.02)(6.07633) cos(1.73696rad ) = 3.37695
3

Resolucin de la Jacobiana J12

Fuera de la Diagonal Vj

V3

P 2
Q 2
=
= 1.02513
V 3
3

V2

P 2
Q 3
=
= 1.02513
V 2
2

Vj

Qi
P1
=
2 j
Vj

En la Diagonal Principal:

Pi Q 3
2
=
+ (2) Vi Yii cos(ii ) =
Vj 3

V3

P3 Q 3
2
80.36
=
+ (2) V 3 Y 33 cos( 33) = 3.37695 + (2)(1.0) 2 (19.5166) cos
= 9.098
180
V 3 3

V2

P 2 Q 2
2
+ (2) V 2 Y 22 cos( 22) = 2.62437 + 2 | 1.02 | 2 . | 15.1879 | cos(1.40577 rad )
=
V 2 2

V2

P 2 Q 2
2
=
+ (2) V 2 Y 22 cos( 22) = 2.56731
V 2 2

Resolucin de la Jacobiana J22

Fuera de la diagonal

V3

Q 2 P 2
= 6.11248
=
V 3 3

V2

Q 3 P 3
= 6.11248
=
V 2 2

En la Diagonal

Vj

Q1 Pi
=
Vj j

Vj

Q1 Pi
2 Vi Yiisen(ii )
=
Vj j

V3

Q 3 P 3
2
80.36
=
2 V 3 Y 33sen( 33) = 20.1352 2(10) 2 *19.5166 sen
= 18.3475
V 3 3
180

V2

Q 2 P 2
2
2 V 2 Y 22sen( 22) = 15.648 2(1.02) 2 *[Link]( 1.40577 ) = 15.5256
=
V 2 2

Entonces construyendo los elementos de la matriz, el Jacobiano numrico para la

primera iteracin ser:

6.11248 2.56731 1.02513

15.648
6.11248 20.1352 1.02513
9.098
= Jacob
2.62437 1.02513
15.5256 6.11248
1.02513 3.37695 6.11248 18.3475

A continuacin se presenta el jacobiano numrico obtenido mediante el software

de simulacin PowerWorld:

FiguraN6: Interfaz del software de simulacin PowerWorld mostrando el resultado de la simulacin del
Jacobiano numrico del sistema de la figura N3.

Entonces comparando el resultado de la simulacin respecto a la obtenida

mediante la aplicacin a mano del mtodo de Newton-Rapshon, se observa una
aproximacin en los valores. Las diferencias radican nuevamente en que existe la
posibilidad que el proceso iterativo la potencia reactiva de generacin en la barra 2
supere los limites que este tiene establecido.

6.- Implementacin del software de programacin WinFlu para la comparacin

con los resultados del programa PowerWorld.

FiguraN6: Interfaz del software de simulacin WinFlu mostrando la resolucin del sistema de potencia
en estudio mostrado en la Figura N1

FiguraN7: Reporte de simulacin presentando los resultados del flujo de potencia del sistema mostrado
en Figura N1

Las diferencias que se presentan en los resultados de la simulacin por el

Software PowerWorld contemplaran diferentes tipos de configuraciones del programa,
as como tambin, la cantidad de datos nominales que requiere cada programa para
correr las simulaciones respectivas.