EXAMEN DE MESA
CARRERA: INGENIERIA ELECTRICA
ASIGNATURA:SISTEMAS DE POTENCIA II
DOCENTE: RENE F. SOLIZ MARTINEZ
PERIODO: II-2010
Fecha: 9 - Sep - 2010
1. Calcular lo parametros en por unidad Base 100 MVA del siguiente sistema
25 Puntos
U1 = 115 kV
U2 = 10 kV Un = 115 kV Un = 115 kV
Sng=70 MVA Sn = 75 MVA L = 82 km Sn3f = 12 MVAR
R1g = 0.007 pu Uk% = 8 % R1 = 0.13 ohm/km
X1g = 0.2 pu Pcc = 120 kW X1 = 0.41 ohm/km
tap = 0.8 pu B1 = 2.7 micro S/km
Para el transformador, tomar en cuenta:
2. Describa los pasos del algoritmo de Newton - Raphson para el calculo de flujos de carga
20 Puntos
3. Realizar los siguientes calculos para la determinación de los flujos de carga del siguiente sistema 55 Puntos
LINEA R X B/2 Y linea
pu pu pu pu
LINEA 1-2 0.01450 j 0.09078 j 0.10529 1.715721 +j* -10.741596
LINEA 2-3 0.00893 j 0.05587 j 0.06479 2.789577 +j* -17.45282
LINEA 3-1 0.01302 j 0.07419 j 0.08240 2.294808 +j* -13.076177
BARRA Carga Vector solucion
pu pu < radianes
Barra 1 0.6+j*0.32 1.05<0
Barra 2 1.15+j*0.67 0.959539<-0.105915
Barra 3 1.8+j*1.23 0.948095<-0.112381
a. Hallar el Ybus 15 puntos
b. Hallar las pérdidas de potencia activa en la línea 1-3 10 puntos
c. Hallar las pérdidas de potencia reactiva en la línea 1-3 10 puntos
d. Hallar la potencia activa entregada por el generador Slack 10 puntos
e. Hallar la potencia reactiva entregada por el generador Slack 10 puntos
100 Puntos
Total
Generador
Zg 0.007+j*0.2 en Base 70 MVA
0.01 +j* 0.28571 en base 100 MVA 4
Transformador
R% 0.16
Z% 8
X% 7.99839984
Zt % 0.16+j*7.9984 en Base 75 MVA
Zt % 0.21333 +j* 10.66453
en base 100 MVA
Zt pu 0.00213 +j* 0.10665en base 100 MVA 4
para el equivalente en susceptancias
Yik 0.18719 +j* -9.37273en base 100 MVA
Kti 0.3125
Ktk -0.25
Yp ik 0.23399 +j* -11.71591 4
Yshu i 0.0585 +j* -2.92898
Yshu k -0.0468 +j* 2.34318
Linea
Zesp 0.13+j*0.41
Ztotal 10.79 +j* 34.03
Zbase 132.25
Zt pu 0.081588 +j* 0.257316
en base 100 MVA 4
Bt 221.4
Bt pu 0.02928
Bt/2 pu 0.01464 en base 100 MVA 4
Capacitor
Sn3f 0+j*-12 MVA
Vlinea 115 KV
Sn3f conj 0+j*12
Zcap f 0 +j* -1102.08333 ohm
Zbase 132.25
Zcap pu 0 +j* -8.333333 en base 100 MVA
Ycap pu 0 +j* 0.12 en base 100 MVA 5
ALGORITMO DE NEWTON RAPHSON
1. Lectura de Datos
2. Formar el YBUS
3. definir el tipo de Barras SLACK, PV, PQ
4. Calcular potencias Especificadas
Pesp = P generacio − Pc arg a
Q esp = Q generacio − Q c arg a
5. Identificar los datos y las incógnitas (V1< fi1, etc..)
6. Adoptar un valor para el vector de estado (flat start)
7. Calcular las potencias de inyección
N
Pkcalc = ∑Vk *Vm[Gkm * Cos(θk − θm) + Bkm * Sen(θk − θm)]
m =1
N
Qkcalc = ∑Vk *Vm[Gkm * Sen(θk − θm) − Bkm * Cos(θk − θm)]
m=1
8. Con 4 y 7 calcular el error en la inyección
∆ P = P esp − Pk calc
∆Q = Q esp − Qk calc
9. Si el error es menor a Epsilon (error máximo aceptable) ir a 15, caso contrario continuar
10. Formar matriz Jacobiano Reducido
11. Invertir el Jacobiano Reducido
12. Hallar la corrección para el vector de estado (delta VE), resolviendo la ecuación matricial
−1
∆θ ∂∂Piθj ∂Pi
∂Vj Vj ∆P
∆V = ∂Qi ∂Qi *
Vj ∆Q
∂θj ∂Vj
13. Corregir el vector de estado
θ k
= θ k −1
+ ∆θ
k −1
V k
= V + ∆V
14. Ir a 7
15. Calcular Inyecciones finales de Generación (barras Slack y PV)
16. Calcular Flujos de Carga en líneas y transformadores
17. Calculo de pérdidas
LINEA Z B/2 Y
pu pu 1 2 3
LINEA 1-2 0.0145+j*0.09078 0+j*0.10529 1.715721 +j* -10.741596 1 4.0105+j*-23.6301 -1.7157 +j* 10.7416 -2.2948 +j* 13.0762
LINEA 2-3 0.00893+j*0.05587 0+j*0.06479 2.789577 +j* -17.45282 2 4.5053+j*-28.0243 -2.7896 +j* 17.4528
LINEA 3-1 0.01302+j*0.07419 0+j*0.0824 2.294808 +j* -13.076177 3 5.0844+j*-30.3818
BARRA Carga Vector solucion
pu pu < radianes
Barra 1 0.6+j*0.32 1.05<0
Barra 2 1.15+j*0.67 0.959539<-0.105915
Barra 3 1.8+j*1.23 0.948095<-0.112381
V fi rad fi deg
1 1.05 0.00000 0
2 0.95954 -0.10592 -6.068768966
3 0.9481 -0.11238 -6.438899702
fi k -fi m V1 mag fi 1 mag V i mag fi i mag
0 1.05 0.00000 1.05 0.00000 G11 4.0105 B11 -23.6301
0.10592 1.05 0.00000 0.95954 -0.10592 G12 -1.7157 B12 10.7416
0.11238 1.05 0.00000 0.9481 -0.11238 G13 -2.2948 B13 13.0762
4.4215763 V1 V1 ( G11 Cos (f1 - f1) + B11 seno( f1-f1))
-0.574749 V1 V2 ( G12 Cos (f1 - f2) + B12 seno( f1-f2))
-0.810254 V1 V3 ( G13 Cos (f1 - f3) + B13 seno( f1-f3)
Piny 1 3.0365737 Piny 1 = Pesp 1 = P gen - P carga
Pcarga 1 0.6
Pgen 1 3.6365737
Qiny 1 1.9163316
Qcarga 1 0.32
Qgen 1 2.2363316
[ ]
Pij = −Gij * Vi 2 − Vi * Vj * Cos(δi − δj ) + Vi * Vj * Bij * Sen(δi − δj ) + Vi 2 *
Gij '
2
fi k -fi m V1 mag fi 1 mag V i mag fi i mag
0.11238 1.05 0.00000 0.9481 -0.11238 G13 -2.2948 B13 13.0762
-0.11238 0.9481 -0.11238 1.05 0.00000 G31 -2.2948 B31 13.0762
Pij de 1 a 3 1.719763293 0.259942632 1.459820661
de 3 a 1 -1.667113973 -0.207293312 -1.459820661
perdidas 0.05264932
Qij = −Vi 2 *
Bij ' sh _ linea
2
[ ]
+ Bij * Vi 2 − Vi * Vj * Cos(δi − δj ) + Vi * Vj * Gij * Sen(δi − δj )
Qij de 1 a 3 1.134235852
de 3 a 1 -0.999128712
perdidas 0.135107139
ALGORITMO DE NEWTON RAPHSON
1. Lectura de Datos
2. Formar el YBUS
3. definir el tipo de Barras SLACK, PV, PQ
4. Calcular potencias Especificadas
a. P esp = P generacio − P c arg a
b. Qesp = Q generacio − Qc arg a
5. Identificar los datos y las incógnitas (V1< fi1, etc..)
6. Adoptar un valor para el vector de estado (flat start)
7. Calcular las potencias de inyección
N
a. Pk calc = ∑Vk *Vm[Gkm * Cos(θk − θm) + Bkm * Sen(θk − θm)]
m =1
N
b. Qk calc = ∑Vk *Vm[Gkm * Sen(θk − θm) − Bkm * Cos(θk − θm)]
m =1
8. Con 4 y 7 calcular el error en la inyección
∆P = Pesp − Pk calc
a.
∆Q = Qesp − Qk calc
9. Si el error es menor a Epsilon (error máximo aceptable) ir a 15, caso contrario
continuar
10. Formar matriz Jacobiano Reducido
11. Invertir el Jacobiano Reducido
12. Hallar la corrección para el vector de estado (delta VE), resolviendo la ecuación
matricial
−1
∆θ ∂∂θj ∂Pi
Vj ∆P
Pi
∂Vj
a. = ∂Qi ∂Qi *
∆V ∂θj ∂Vj Vj ∆Q
13. Corregir el vector de estado
a. θ k = θ k −1 + ∆θ
b. V k = V k −1 + ∆V
14. Ir a 7
15. Calcular Inyecciones finales de Generación (barras Slack y PV)
16. Calcular Flujos de Carga en líneas y transformadores
17. Calculo de pérdidas
[ ]
Pij = −Gij * Vi 2 − Vi * Vj * Cos(δi − δj ) + Vi * Vj * Bij * Sen (δi − δj ) + Vi 2 *
Gij ' sh _ linea
2
Qij = −Vi 2 *
Bij ' sh _ linea
2
[ ]
+ Bij * Vi 2 − Vi * Vj * Cos(δi − δj ) + Vi * Vj * Gij * Sen (δi − δj )
