Modelo de Flujo de Carga Desbalanceado para El Análisis de Redes de Distribución Con Generación Distribuida PDF
Modelo de Flujo de Carga Desbalanceado para El Análisis de Redes de Distribución Con Generación Distribuida PDF
Modelo de Flujo de Carga Desbalanceado para El Análisis de Redes de Distribución Con Generación Distribuida PDF
Facultad de Ingeniería
Escuela de Ingeniería Mecánica Eléctrica
FACULTAD DE INGENIERÍA
TRABAJO DE GRADUACIÓN
AL CONFERÍRSELE EL TÍTULO DE
I
3. MODELO DE CARGAS .......................................................................... 33
3.1. Cargas conectadas en Y .......................................................... 34
3.1.1. Cargas de potencia constante ................................. 35
3.1.2. Cargas de impedancia constante ............................ 36
3.1.3. Cargas de corriente constante ................................ 37
3.2. Cargas conectadas en Delta .................................................... 38
3.2.1. Cargas de potencia constante ................................. 39
3.2.2. Cargas de voltaje constante .................................... 39
3.2.3. Cargas de corriente constante ................................ 40
3.3. Bancos de capacitores ............................................................. 41
3.4. Cargas distribuidas................................................................... 43
II
4.2.2.2. Reguladores de voltaje conectados
en estrella ............................................ 97
4.2.2.3. Reguladores de voltaje conectados
en delta .............................................. 100
4.2.2.4. Reguladores de voltaje conectados
en delta abierta .................................. 102
III
7.7. Validación del modelo de generadores .................................. 154
IV
ÍNDICE DE ILUSTRACIONES
FIGURAS
V
16. Diagrama fasorial de la conexión Delta-Y aterrizada Step Up
Connection ............................................................................................ 57
17. Conexión Delta-Y aterrizada Step Down............................................... 62
18. Diagrama fasorial de la conexión Delta-Y aterrizada Step Down
Connection ............................................................................................ 63
19. Conexión Y- Delta Step Up ................................................................... 65
20. Diagrama fasorial de la conexión Y-Delta Step Up ............................... 66
21. Conexión Y-D Step Down Connection .................................................. 71
22. Diagrama fasorial de voltajes en la conexión Y-D Step Down .............. 72
23. Conexión (izquierda) y diagrama fasorial (derecha) de la conexión
de transformadores Y-Y ........................................................................ 74
24. Conexión (izquierda) y diagrama fasorial (derecha) de la conexión
de transformadores Delta-Delta ............................................................ 75
25. Conexión (izquierda) y diagrama fasorial (derecha) de la conexión
Y abierta- Delta abierta Step Up Connection usando las fases A y B
en el lado del primario ........................................................................... 78
26. Conexión (izquierda) y diagrama fasorial (derecha) de la conexión
Y abierta- Delta abierta Step Down Connection usando las fases A
y B en el lado del primario..................................................................... 82
27. Esquema del circuito para regulador por etapas ................................... 87
28. Circuito equivalente (izquierda) y circuito equivalente abreviado
(derecha) de un regulador de voltaje en su conexión tipo A en
posición aumentadora ........................................................................... 88
29. Circuito equivalente (izquierda) y circuito equivalente abreviado
(derecha) de un regulador de voltaje en su conexión tipo B en
posición aumentadora ........................................................................... 89
30. Esquema básico de un compensador de caída de línea ...................... 91
31. Reguladores de voltaje conectados en estrella..................................... 97
32. Reguladores de voltaje conectados en delta ...................................... 100
VI
33. Reguladores de voltaje conectados en delta abierta .......................... 103
34. Uso del multiport based method para crear la matriz sensitiva .......... 113
35. Algoritmo para el cálculo iterativo de inyección de potencia reactiva
de los generadores tipo PV ................................................................ 118
36. Reordenamiento de los nodos según el número de capas ................ 123
37. Representación de las matrices generalizadas para cualquier
elemento entre dos nodos consecuentes ........................................... 127
38. Algoritmo de solución de flujo de carga en redes de distribución ....... 130
39. Comportamiento de la convergencia en valores por unidad del
voltaje en el nodo 730 ........................................................................ 156
40. Comportamiento de la convergencia en valores por unidad del
voltaje en el nodo 730 ........................................................................ 157
41. Comportamiento de la convergencia en valores por unidad del
voltaje en el nodo 734 ........................................................................ 157
42. Comportamiento de la convergencia en valores por unidad del
voltaje en el nodo 720 ........................................................................ 158
43. Comportamiento de la convergencia de voltaje en los nodos 730,
734 y 720 del caso IEEE 37 Node test feeder .................................... 159
44. Comparación de los perfiles de voltaje de la fase A del caso
IEEE37 node test feeder con tres escenarios diferentes .................... 162
45. Comparación de los perfiles de voltaje de la fase B del caso
IEEE37 node test feeder con tres escenarios direferentes ................ 163
46. Comparación de los perfiles de voltaje de la fase C del caso
IEEE37 node test feeder con tres escenarios direferentes ................ 164
VII
TABLAS
VIII
XVI. Comparación entre los voltajes y las corrientes un regulador de
voltaje de tipo B conectado en estrella con una regulación
independiente en las fases ................................................................. 145
XVII. Comparación entre los voltajes y las corrientes un regulador de
voltaje de tipo B conectado en delta abierta entre las fases AB y CB
con una regulación independiente en las fases.................................. 146
XVIII. Comparación de los perfiles de voltaje del caso IEEE37 node test
feeder obtenidos con el programa y los presentados en el caso
de prueba ........................................................................................... 149
XIX. Comparación entre las corrientes obtenidas con el programa y las
presentadas en el caso de prueba IEEE37 node test feeder ............. 150
XX. Flujo de potencia activa y reactiva caso IEEE37 node test feeder ..... 151
XXI. Consumo de potencia aparente de las cargas caso IEEE37 node
test feeder .......................................................................................... 152
XXII. Comparación entre las pérdidas halladas con el programa y las
presentadas en el caso de prueba IEEE37 node test feeder ............. 153
XXIII. Comparación de los perfiles de voltaje caso IEEE37 node test
feeder con tres escenarios diferentes ................................................ 161
IX
X
LISTA DE SÍMBOLOS
Símbolo Significado
Y Admitancia en derivación.
Ancho del blindaje en un conductor subterráneo de
tipo tape shield.
Ángulo entre un par de líneas desde el conductor i
hacia su imagen hasta la imagen del conductor j.
Angulo de factor de potencia.
Ángulo de voltaje.
A Área transversal de un conductor.
Campo eléctrico.
Capacitancia.
Carga eléctrica.
Coeficiente de potencial.
Coeficiente de resistividad con la temperatura.
Conductividad.
Corriente eléctrica.
Corriente eficaz oRMS (Raíz cuadrada media).
Corriente entre las fases x y y de una conexión en
delta, donde x y y pueden ser cualquier
combinación de fases.
Corriente de línea.
d Diámetro.
Distancia.
XI
Distancia entre los conductores i y j en pies.
Distancia media geométrica mutua.
Distancia media geométrica mutua entre los
conductores i y j.
Distancia media geométrica propia.
Enlaces de flujo.
Flujo magnético.
Frecuencia angular del sistema en radianes por
segundo.
Frecuencia en Hertz.
Fuerza contra electromotriz.
Z Impedancia.
Impedancia de capacitores, es una reactancia
capacitiva.
Impedancia propia del conductor i en /milla.
Impedancia mutua entre los conductores i y j en
/milla.
Inductancia.
Inyección de corriente reactiva por capacitores.
Inyección de corriente de un generador.
Longitud.
Matriz de admitancias en derivación.
Matriz de capacitancias de fase.
Matriz de coeficientes de potencia del fase.
Matriz de coeficientes de potencial primitivos.
Matriz de coeficientes de potencial primitivos entre
conductores.
XII
Matriz de coeficientes de potencial primitivos entre
conductores de fase y neutrales.
Matriz de coeficientes de potencial primitivos entre
neutrales.
Matriz de coeficientes de potencial primitivos entre
neutrales y conductores de fase.
Matriz de impedancias de fase.
Matriz de impedancias primitivas.
Matriz de impedancias primitivas entre conductores.
Matriz de impedancias primitivas entre conductores
de fase y neutrales.
Matriz de impedancias primitivas entre neutrales.
XIII
Matriz generalizada que relaciona el vector de
voltajes de un nodo m con el vector de corrientes de
el nodo siguiente n.
Matriz generalizada que relaciona el vector de
voltajes de un nodo m con el vector de voltajes de el
nodo siguiente n.
Matriz generalizada que relaciona el vector de
voltajes de un nodo n con el vector de corrientes de
su nodo anterior m.
Matriz generalizada que relaciona el vector de
voltajes de un nodo n con el vector de voltajes de su
nodo anterior m.
Matriz sensitiva.
Número de espiras del devanado x de un
transformador, x puede ser uno o dos.
Número de hilos neutrales en un conductor
subterráneo de tipo cocentric neutral.
Número de iteración.
Número de vueltas de una bobina.
Operador .
Operador .
Operador /milla.
Operador .
µ Permeabilidad.
Permitividad.
Potencia aparente.
Profundidad de penetración o piel.
XIV
Radio.
Radio del conductor i en pies.
Radio medio geométrico.
Radio medio geométrico en pies.
X Reactancia.
Relación de transformación.
Relación de transformación para reguladores de
voltaje por etapas.
Resistencia efectiva en AC del conductor i en
/milla.
R Resistencia efectiva en corriente alterna.
Resistencia en corriente directa.
Resistividad.
Resistividad de la tierra en -metro.
T Temperatura.
Vector de corrientes de línea para el nodo m.
Vector de voltajes de línea a neutro para el nodo m.
Voltaje o potencial eléctrico.
XV
XVI
GLOSARIO
XVII
Arranque plano Se asume que todos los nodos de un sistema
poseen su voltaje nominal en el inicio de la solución
de flujo de carga.
XVIII
Cargas de potencia Cargas que mantienen su consumo de potencia
constante constante sin importar los cambios en el voltaje del
sistema.
XIX
Convergencia Criterio usado para considerar aceptable la solución
de un método numérico, comúnmente se usa como
criterio de convergencia la diferencia entre los
resultados de dos iteraciones consecuentes.
XX
Distancia media Raiz -esima de los productos de las distancias de
geométrica propia cada hilo del conductor hacia sí mismo y a cada uno
(DMG propia) de los otros hilos, también llamada radio medio
geométrico (RMG).
XXI
Equivalente de Circuito equivalente de una red que permite
Thevenin relacionar el voltaje y la corriente presentes en un
punto determinado.
XXII
Generación Generación de escala media o pequeña conectada
distribuida directamente al sistema de distribución.
XXIII
Histéresis Fenómeno creado en el núcleo ferromagnético del
transformador por el cambio continuo de los
dominios magnéticos del mismo y un retraso en su
reorientación.
XXIV
Ley de Lenz Oposición de voltajes inducidos a los cambios en la
corriente que los induce.
XXV
Matriz de admitancias Matriz que contiene el las admitancias en derivación
en derivación equivalentes para un segmento de línea.
XXVI
Método de Método desarrollado por Fortescue el cual establece
componentes que los voltajes y corrientes de un sistema trifásico
simétricas desbalanceado pueden ser representados por un
grupo de tres sistemas trifásicos balanceados
llamados componentes simétricos.
XXVII
Permeabilidad Constante física que relaciona el campo magnético
en un objeto con el material.
XXVIII
Red radial Red de distribución con alimentación únicamente
por uno de sus extremos y la transmisión de energía
a la carga de forma radial desde el alimentador.
XXIX
Resistencia Oposición de un material al flujo de corriente
eléctrica.
XXX
Tolerancia Valor máximo aceptado para la magnitud del error
usado como criterio de convergencia en la solución
de métodos numéricos.
XXXI
XXXII
RESUMEN
XXXIII
Para la construcción del modelo se clasifican los componentes según su
conexión con los nodos; en serie o en derivación (también llamados en
paralelo). Los elementos en serie requieren para su análisis el uso de matrices
generalizadas para relacionar los voltajes y corrientes de dos nodos
consecuentes y así poder considerar cada elemento de la red como un bloque,
al que entran y salen corrientes.
XXXIV
OBJETIVOS
General
Específicos
XXXV
XXXVI
INTRODUCCIÓN
XXXVII
generadores conectados a estas redes, por lo que se desconocen el
comportamiento y la respuesta del sistema al agregar una cantidad
relativamente grande de generación desde la red de distribución.
1
KERSTING, William. “Distribution System Modeling and Analysis”.
XXXVIII
1. SISTEMAS DE DISTRIBUCIÓN
1.1. Generalidades
1
Los tipos más comunes de redes eléctricas son:
2
solar y eólica, que pueden conectarse directamente en los puntos de consumo,
lo que cambia radicalmente la estructura de los sistemas de distribución,
transformándolos de un sistema clásicamente radial a uno en anillo cerrado con
más de una fuente de alimentación.
3
Seccionadores: se encargan de dividir una red eléctrica en partes, para
coordinar la operación de alimentación de energía eléctrica. Se usan
también para aislar áreas en estado de falla.
4
para los equipos que evitan que la descarga llegue a las líneas y equipos
de la subestación y pararrayos a los equipos que drenan a tierra la
corriente de una descarga electroatmosférica.
5
1.2. Matrices generalizadas
6
del otro lado salen corrientes a otro nivel de voltaje (nodo N). La matriz [a]
relaciona el vector de voltajes del nodo m con el vector de voltajes del nodo n,
la matriz [b] relaciona el vector de voltajes del nodo m con el vector de
corrientes que salen del nodo n, la matriz [c] relaciona el vector de corrientes
que entran al nodo m con el vector de voltajes del nodo n y la matriz [d]
relacionan el vector de corrientes que entran al nodo m con el vector de
corrientes que salen del nodo n (1-1).
(1-1)
También es posible relacionar las variables del nodo n (aguas abajo) con
las variables del nodo m:
(1-2)
Donde [A] y [B] son matrices que relacionan el vector de voltajes del nodo
n con los vectores de voltajes y corrientes del nodo m, [C] y [D] son matrices
que relacionan el vector de corrientes del nodo n con las variables del nodo m.
7
actualizar los valores del voltaje considerando las corrientes calculadas en la
iteración anterior.
(1-3)
Donde [A] y [B] son matrices que relacionan el vector de voltajes del nodo
n con el vector de voltajes del nodo m y el vector de corrientes que entran al
nodo n. En este caso no se incluyen las relaciones para hallar las corrientes
entrantes al nodo n, ya que se usan las corrientes calculadas en la iteración
anterior, en el capítulo 6 se explica el algoritmo de solución de flujo de carga
que aclara este hecho.
8
2. MODELO DE LÍNEA DE DISTRIBUCIÓN
9
Las líneas de distribución no siempre son líneas transpuestas, lo que
acentúa el desbalance entre las fases, mientras que las líneas de
transmisión son casi siempre transpuestas.
Estas diferencias hacen evidente que se debe usar un modelo más exacto
para líneas de distribución, considerando las propiedades de la línea.
10
2.1.1. Tipos de conductores
11
2.1.2. Resistencia
(2-1)
(2-2)
(2-3)
(2-4)
12
La distribución uniforme de la corriente en la sección transversal de un
conductor, solamente se presenta en corriente directa, en cambio para corriente
alterna conforme aumenta la frecuencia se acentúa más la no uniformidad de la
corriente en la sección, a este fenómeno se llama efecto piel y en conductores
circulares la densidad de corriente se incrementa desde el interior hacia la
superficie. Los flujos alternos de corriente inducen mayores voltajes en el
interior del conductor que en la superficie, y por la Ley de Lenz los voltajes
inducidos se oponen a los cambios de corriente que los producen, creando
voltajes inducidos más elevados en el interior del conductor causando una
mayor densidad de corriente en la superficie.
(2-5)
2
El desarrollo matemático de la profundidad del efecto piel se encuentra en HYAT, William, Teoría
electromagnética. p. 369-376.
13
de sus conductores. Para fines prácticos se puede usar un valor de resistencia
efectiva que toma en cuenta todos los factores explicados anteriormente:
(2-6)
2.1.3. Inductancia
(2-7)
(2-8)
14
Todos los hilos de un conductor tienen el mismo radio y una fracción de la
corriente total I/n.
(2-9)
(2-10)
(2-11)
15
Figura 3. Configuración de conductores representativa para las
ecuaciones (2-10) y (2-11)
16
John Carson4 desarrolló a partir de la ecuación (2-9) y el método de las
imágenes ecuaciones para calcular la impedancia propia (2-12) y mutua (2-13)
de las líneas para cualquier configuración de fases, incluyendo el retorno de la
corriente por tierra, considerando a esta como un conductor sólido, uniforme,
infinito, con una superficie plana y una resistividad constante.
(2-12)
(2-13)
Donde:
4
CARSON, John, Wave propagation in overhead wires with ground return, Bell System Technical Journal.
Vol. 5. p. 65-80.
17
(2-14)
(2-15)
(2-16)
(2-17)
18
Para el modelo de la línea se usan las ecuaciones modificadas de
Carson5, en las que se usa únicamente el primer término de P y Q y se
sustituyen los valores de -metro, y Hertz:
(2-18)
(2-19)
5
Una explicación amplia del desarrollo de las ecuaciones modificadas de Carson puede ser consultada en:
CARSON, John, Wave propagation in overhead wires with ground return, Bell System Technical Journal.
Vol. 5. p. 65-80.
.
19
Esta formulación considera que el retorno de la corriente y la
concatenación de flujos magnéticos entre fases y entre fases y neutro. Es
posible utilizarlas para crear un modelo de línea que considera los desbalances
de corrientes y su efecto en la línea, ya que al cambiar la corriente en
cualquiera de las fases cambia su campo magnético asociado y esto afecta a
las otras fases del sistema. Este efecto se manifiesta como una impedancia
entre fases, tal como se muestra en la figura 4.
20
Con las impedancias de un segmento de línea es posible construir una
matriz de impedancias primitivas de la siguiente forma:
(2-20)
(2-21)
21
(2-22)
(2-23)
El radio medio geométrico del neutral, para los conductores con neutrales
concéntricos es:
(2-24)
(2-25)
22
Donde es la resistencia de un hilo neutral.
(2-26)
(2-27)
23
El radio medio geométrico es el radio del círculo que pasa en medio del
blindaje (figura. 6):
(2-28)
2.1.6.1. Conductancia
El efecto corona ocurre cuando la intensidad del campo eléctrico hace que
el aire se ionice eléctricamente y se vuelva conductor. Las pérdidas debidas a la
conductancia en derivación comúnmente se desprecian, debido a que su efecto
es mucho menor al generado por la resistencia en serie, pero principalmente
debido a la dificultad para su medición, ya que es muy variable al depender de
las condiciones meteorológicas entre otras.
24
2.1.6.2. Capacitancia
(2-29)
(2-30)
(2-31)
25
(2-32)
(2-33)
(2-34)
(2-35)
26
Nuevamente se usa el método de las imágenes, para considerar el efecto
que crea el suelo sobre el campo eléctrico de los conductores; considerando
que el suelo es un conductor uniforme, perfecto y con una longitud infinita, en
este caso se considera que los conductores tienen una carga definida para un
instante dado y que las imágenes contienen cargas iguales a las de los
conductores pero con signos opuestos. Con este método se puede formar
primero una matriz de coeficientes de potencial primitivos (2-36) y (2-37), para
después aplicar la reducción de Kron (2-38) y tener una matriz de coeficientes
de potencial útil para el modelo (2-39), de forma análoga a la formación de la
matriz de impedancias de fase.
(2-36)
O de forma compacta
(2-37)
(2-38)
(2-39)
27
(2-40)
(2-41)
(2-42)
Para el caso del conductor subterráneo con blindaje tipo cinta se puede
considerar como un conductor con una cantidad infinita de hilos neutrales, por
lo que al aplicar (2-42) el segundo término del denominador se aproxima a cero:
(2-43)
28
En el caso de conductores subterráneos se debe tener en cuenta que la
permitividad del medio no es aproximadamente igual a 1, sino que es la
permitividad del material aislante entre el conductor de fase y el neutral. Para el
caso de las capacitancias mutuas entre fases es evidente que se vuelven nulas,
ya que el campo eléctrico se queda confinado al interior de los neutrales, por lo
que los términos fuera de la diagonal de la matriz de capacitancias son cero.
(2-44)
29
Figura 8. Modelo de línea de distribución
(2-45)
(2-46)
(2-47)
6
El desarrollo de estas ecuaciones está ampliamente explicado en: KERSTING, William. Distribution
System Modeling and Analysis. p. 125-129. Y para líneas de transmisión balanceadas en: GLOVER,
Duncan, SARMA, Mulukutla y OVERBYE Thomas. Power System Analisis and Design. p. 145-192.
.
30
(2-48)
(2-49)
(2-50)
Donde:
(2-51)
31
32
3. MODELO DE CARGAS
33
Potencia activa y reactiva constante
Corriente constante
Impedancia constante
Cualquier combinación de las anteriores
34
Figura 9. Cargas conectadas en Y
35
(3-1)
(3-2)
36
Para obtener el consumo de corriente en este caso, es necesario usar la
ley de Ohm, que es aplicable únicamente para el caso de cargas con
impedancia constante:
(3-3)
(3-4)
37
3.2. Cargas conectadas en Delta
(3-5)
38
3.2.1. Cargas de potencia constante
(3-6)
(3-7)
39
Para encontrar el consumo de corriente en un momento dado hay que
valerse nuevamente de la ley de Ohm, y usar la impedancia entre fases
encontrada con las ecuaciones (3-7) y el voltaje aplicado entre fases:
(3-8)
(3-9)
40
3.3. Bancos de capacitores
41
(3-10)
(3-11)
(3-12)
42
Y la corriente de los capacitores conectados en delta es:
(3-13)
43
Figura 11. Modelo equivalente de cargas distribuidas
44
Tabla I. Tipos de cargas
45
46
4. MODELO DE TRANSFORMADORES Y REGULADORES DE
VOLTAJE
4.1. Transformadores
47
Los factores que deben tomarse en cuenta para el modelado de
transformadores son:
48
Figura 12. Modelo equivalente exacto del transformador monofásico
(4-1)
(4-2)
49
(4-3)
(4-4)
De forma matricial:
(4-5)
(4-6)
(4-7)
(4-8)
(4-9)
50
(4-10)
(4-11)
(4-12)
(4-13)
7
KERSTING, William. Distribution System Modeling and Analysis. p. 147-150.
51
Figura 13. Modelo aproximado del transformador monofásico
(4-14)
(4-15)
(4-16)
(4-17)
(4-18)
(4-19)
Al comparar las ecuaciones (4-6)-(4-9), (4-12) y (4-13) con las (4-14)- (4-
19) se observa fácilmente que las últimas son mucho más sencillas, pero
52
conllevan un error de aproximadamente 0,1 por ciento, lo que puede ser
importante en sistemas grandes o situaciones que requieran mucha exactitud
debido a que se trabaja con magnitudes de miles de voltios aún hablando de
sistemas de distribución, que poseen voltajes relativamente bajos.
Delta-Y aterrizada
Y-Delta
Y aterrizada-Y aterrizada
Delta-Delta
Y abierta-Delta abierta
53
corriente generada por el desbalance se queda circulando dentro de la delta. Se
debe tener especial cuidado con el desfasamiento generado por los diferentes
tipos de conexión de los transformadores, ya que cada conexión introduce un
desfase diferente en la onda de voltaje y corriente del secundario respecto al
primario en un múltiplo de 30 grados.
54
Figura 14. Bloque representativo de una conexión cualquiera de
transformadores trifásicos
55
estas presentan. Como se mencionó en la sección anterior es posible conectar
los transformadores de manera que el voltaje y corriente del primario se
adelanten o atrasen respecto de las conexiones del secundario.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders.http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html. Consulta. 15 de septiembre de
2011.
56
Figura 16. Diagrama fasorial de la conexión Delta-Y aterrizada Step Up
Connection
(4-20)
(4-21)
57
Ahora, haciendo uso del método de las componentes simétricas
desarrollado por Fortescue se puede relacionar el voltaje de línea a neutro del
lado primario con el voltaje de línea a neutro del lado del secundario de la
siguiente forma:
(4-22)
(4-23)
Donde:
(4-24)
58
(4-25)
(4-26)
(4-27)
(4-28)
(4-29)
59
De la ecuación (4-29) se hace evidente la obtención de las dos primeras
matrices generalizadas:
(4-30)
(4-31)
Ahora, para obtener los voltajes del lado secundario en función de las
magnitudes del lado primario se debe relacionar los voltajes de línea-línea con
los de línea a neutro del lado primario:
(4-32)
(4-33)
(4-34)
60
De esta forma se encuentran las matrices generalizadas [A] y [B]:
(4-35)
(4-36)
(4-37)
(4-38)
(4-39)
61
(4-40)
(4-41)
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders.http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html. Consulta. 15 de septiembre de
2011.
62
Figura 18. Diagrama fasorial de la conexión Delta-Y aterrizada Step
Down Connection
(4-42)
8
El proceso para la obtención de estas matrices en la conexión Step down está ampliamente explicado en:
KERSTING, William. Distribution System Modeling and Analysis. p. 201-208.
63
(4-43)
(4-44)
(4-45)
(4-46)
(4-47)
64
Figura 19. Conexión Y- Delta Step Up
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders.http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html. Consulta. 15 de septiembre de
2011.
65
Figura 20. Diagrama fasorial de la conexión Y-Delta Step Up
(4-48)
66
(4-49)
(4-50)
Reordenando la ecuación 4-48 y despejando para los voltajes del lado del
secundario:
(4-51)
Relacionando las corrientes dentro de la delta del secundario con las que
salen del transformador se puede formar la siguiente relación:
(4-52)
La matriz que 3x3 en la ecuación (4-52) es una matriz singular, por lo que
no tiene una matriz inversa, debido a esto es necesario modificar esta ecuación.
(4-53)
67
De esta manera se puede incluir otra ecuación independiente que
relaciona las corrientes dentro y fuera de la delta:
(4-54)
(4-55)
Se pueden relacionar las corrientes que entran a los devanados con las
que salen de los mismos:
(4-56)
(4-57)
(4-58)
68
(4-59)
(4-60)
(4-61)
(4-62)
(4-63)
(4-64)
69
(4-65)
(4-66)
(4-67)
Sustituyendo (4-67) en (4-66) y despejando para los voltajes del lado del
secundario:
(4-68)
(4-69)
70
(4-70)
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders.http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html. Consulta. 15 de septiembre de
2011.
9
Las matrices generalizadas para la conexión Y-Delta Step down connection están ampliamente
desarrolladas en: KERSTING, William. Distribution System Modeling and Analysis. p. 212-218.
71
Figura 22. Diagrama fasorial de voltajes en la conexión Y-D Step Down
(4-71)
(4-72)
(4-73)
72
(4-74)
(4-75)
(4-76)
73
Figura 23. Conexión (izquierda) y diagrama fasorial (derecha) de la
conexión de transformadores Y-Y
(4-77)
(4-78)
(4-79)
(4-80)
10
Las matrices generalizadas para la conexión Y-Y están ampliamente desarrolladas en: KERSTING,
William. Distribution System Modeling and Analysis. p. 222-224.
74
(4-81)
(4-82)
75
Las matrices generalizadas para la conexión Delta-Delta11 son:
(4-83)
(4-84)
(4-85)
(4-86)
(4-87)
(4-88)
11
Las matrices generalizadas para la conexión Delta-Delta están ampliamente desarrolladas en:
KERSTING, William. Distribution System Modeling and Analysis. p. 224-231.
76
Las ecuaciones (4-83) y (4-88) se expresan como un producto indicado
porque no se pueden reducir a una forma más simple.
Esta conexión puede ser usada para alimentar una combinación de cargas
monofásicas y trifásicas, pero es más común en la alimentación de una carga
monofásica pesada y una carga trifásica ligera. Esta conexión muchas veces se
usa en casos de emergencia cuando la conexión Y-Delta pierde un
transformador y se requiere seguir alimentando la carga a una capacidad de
carga reducida, en caso de perder una fase la potencia total del banco se
reduce a 57,7 por ciento de la capacidad nominal y el factor de utilización se
reduce a 86,6 por ciento para los dos transformadores restantes.
77
Únicamente se desarrollan las matrices generalizadas para la conexión de
transformadores Y abierta- Delta abierta conectada en las fases A y B, las
demás conexiones se obtienen de forma equivalente.
(4-89)
78
(4-90)
(4-91)
(4-92)
(4-93)
(4-94)
(4-95)
(4-96)
(4-97)
79
No se considera la corriente en la fase a del lado del secundario por que
los voltajes y del lado del primario se pueden expresar en función
únicamente de las corrientes e . Los voltajes de línea a línea del secundario
pueden ser sustituidos por los voltajes de línea a neutro:
(4-98)
(4-99)
(4-100)
La relación entre las corrientes del lado del primario y el secundario es:
(4-101)
(4-102)
(4-103)
80
Usando ahora (4-92), sustituyendo el voltaje Vca que debe ser igual a la
suma del negativo de los dos voltajes restantes y considerando una
combinación de las ecuaciones (4-93), (4-94) y (4-96), para sustituir las
corrientes de línea a línea con las de línea a neutro:
(4-104)
Ahora usando la ecuación (4-27) se pueden encontrar los voltajes del lado
del secundario en función de las magnitudes del secundario:
(4-105)
(4-106)
(4-107)
81
matrices para la conexión Y abierta-Delta abierta Step Down12 (4-108) a (4-113)
mostrada en la figura 26.
(4-108)
(4-109)
(4-110)
12
El proceso para la obtención de estas matrices en la conexión Step down está ampliamente explicado
en: KERSTING, William. Distribution System Modeling and Analysis. p. 236-240.
82
(4-111)
(4-112)
(4-113)
Donde:
(4-114)
(4-115)
(4-116)
83
(4-117)
84
4.2. Reguladores de voltaje
85
Los autotransformadores están formados por dos devanados, un
devanado en serie que es equivalente al devanado secundario y el devanado
común que es equivalente al devanado primario.
86
Figura 27. Esquema del circuito para regulador por etapas
13
De acuerdo al Standard ANSI/IEEE C57.15-1986.
87
al voltaje del devanado común, ya sea para aumentar o disminuir el voltaje a la
salida regulador.
88
Figura 29. Circuito equivalente (izquierda) y circuito equivalente
abreviado (derecha) de un regulador de voltaje en su
conexión tipo B en posición aumentadora
14
KERSTING, William. Distribution System Modeling and Analysis. p. 168.
89
Tabla III. Relación entre voltajes y corrientes de entrada con los de
salida para las diferentes conexiones de reguladores de
voltaje
Tipo B
(4-118)
90
El número de taps lo determina el compensador de caída de línea (figura
30), que detecta cambios en la corriente representativos de las caídas de
voltaje en la línea por medio de un CT (Transformador de corriente), un PT
(transformador de potencial) y una impedancia de compensación que es
equivalente a la impedancia total que debe existir entre el regulador y el centro
de carga. El compensador detecta si el voltaje en el centro de carga se
encuentra dentro de valores aceptables, si no es así activa un motor encargado
de cambiar la posición de taps en el transformador para ajustar el voltaje
siempre con el propósito de mantener el voltaje en el centro de carga lo más
cercano al nominal.
91
El compensador de caída de línea detecta la diferencia entre el voltaje de
la línea (referido al voltaje usado en la carga) y la caída de voltaje equivalente
representado por la caída de voltaje dada en R’+j X’. Si el voltaje en la bobina
del relé es diferente al que debe existir (que comúnmente es equivalente al de
la carga) el relé cambia los taps del transformador y los ajusta a una posición en
la cual el voltaje en la carga sea aproximadamente igual al requerido.
(4-119)
(4-120)
(4-121)
92
Ahora se transforma la corriente de la línea a una corriente equivalente
referida al lado del compensador:
(4-122)
(4-123)
(4-124)
Esta caída de voltaje representada como R’+jX’ (ver figura 30), es usada
por el circuito compensador para cuantificar las caídas de voltaje hasta el centro
de carga, por lo que dependiendo de la magnitud de R’+jX’, así será la cantidad
de taps que el circuito cambiará para ajustar el voltaje.
Conexión monofásica
Dos reguladores conectados en Y abierta
Tres reguladores conectados en estrella
93
Dos reguladores conectados en delta abierta
Tres reguladores conectados en delta cerrada
(4-125)
(4-126)
(4-127)
(4-128)
94
(4-129)
(4-130)
(4-131)
] (4-132)
95
(4-133)
(4-134)
(4-135)
(4-136)
96
4.2.2.2. Reguladores de voltaje conectados en
estrella
97
En la conexión mostrada (figura 31) cada regulador está conectado en su
posición aumentadora y tiene su propio circuito compensador de caída de línea.
(4-137)
(4-138)
(4-139)
(4-140)
(4-141)
(4-142)
98
Las matrices generalizadas para una conexión de reguladores trifásicos en
Y formada por reguladores con una conexión tipo B son:
(4-143)
(4-144)
(4-145)
(4-146)
(4-147)
(4-148)
99
4.2.2.3. Reguladores de voltaje conectados en delta
100
En la conexión mostrada (figura 32) cada regulador está conectado en su
posición aumentadora y tiene su propio circuito compensador de caída de línea.
(4-149)
(4-150)
(4-151)
(4-152)
(4-153)
(4-154)
(4-155)
101
(4-156)
(4-157)
(4-158)
(4-159)
(4-160)
102
Figura 33. Reguladores de voltaje conectados en delta abierta
(4-161)
(4-162)
(4-163)
103
(4-164)
(4-165)
(4-166)
(4-167)
(4-168)
(4-169)
(4-170)
(4-171)
104
5. MODELO DE GENERACIÓN DISTRIBUIDA
Se llaman tecnologías con recursos renovables a las que son usadas para
la generación de energía eléctrica usando fuentes como:
105
Eólica: energía producida por el caudal de viento.
Geotérmica: energía producida por medio del calor natural de la tierra, que
puede extraerse del vapor, agua, gases, excluidos los hidrocarburos, o a
través de fluidos inyectados para este fin.
Confiabilidad de la red.
106
Seguridad del personal operativo.
Por lo que se requiere que los proyectos sean analizados desde el punto
de vista de la red, considerando el desbalance inherente de las líneas de
distribución y las características distintivas de la red de distribución, para evitar
que se genere cualquiera de estos problemas.
Para conectar los GD’s a la red de distribución según norma técnica para
la conexión, operación, control y comercialización de la Comisión Nacional de
Energía Eléctrica, se requiere de protecciones y líneas de transmisión entre el
generador y la red capaces de soportar la magnitud de los nuevos flujos de
potencia y dependiendo del tipo de generación transformadores. Los relés de
107
protección deberán ser aplicados tanto al generador como a la interconexión
eléctrica y deberán ser coordinados con la empresa eléctrica distribuidora.
108
baterías. Esto permite una inyección de potencia desbalanceada, lo que hasta
cierto punto puede beneficiar al sistema de distribución al minimizar el
desbalance entre las líneas, pero puede perjudicar la calidad de la energía en
ese punto al introducir armónicos en el sistema.
109
mantienen relación entre sí, de esta forma si se conoce la magnitud y ángulo de
una de ellas es posible ubicar las demás respecto a esta.
Los generadores del tipo PQ son generadores que aportan una cantidad
específica de potencia activa y reactiva sin importar el nivel de voltaje, debido a
esto se pueden modelar como cargas negativas inyectando corriente a la red,
ya que no poseen ninguna restricción adicional.
(5-1)
110
Donde δ es el ángulo de voltajes de línea a neutro y θ es el ángulo de
factor de potencia. El signo menos indica que se está inyectando corriente a la
red, mientras que para el caso de las cargas se consideraba el consumo con
signo positivo.
(5-2)
(5-3)
111
5.2.3. Modelo de generadores tipo PV
112
Figura 34. Uso del multiport based method para crear la matriz
sensitiva
(5-4)
113
siendo los elementos dentro de la diagonal la suma de matrices de impedancias
de los elementos entre el generador i y el alimentador y los elementos fuera de
la diagonal son la suma de las matrices de impedancias de los elementos en
común entre los generadores i y j. El vector de corrientes de tamaño (nx1) que
multiplica a la matriz sensitiva contiene n subvectores de (3x1), que
corresponden a los incrementos de la corriente en las tres fases de un nodo. El
vector de tamaño (nx1) de diferencias de voltaje contiene n subvectores de
(3x1) que alojan los incrementos del voltaje en las tres fases de cada
generador.
(5-5)
(5-6)
114
Para la incorporación de los generadores al modelo propuesto se debe
agregar un proceso iterativo independiente del flujo de carga:
(5-7)
(5-8)
115
(5-9)
116
(5-10)
(5-11)
(5-12)
117
que circula en las líneas de los nodos anteriores inyectando corriente desde el
punto donde está conectado. El algoritmo para el cálculo de la inyección de
potencia reactiva para generadores de tipo PV se presenta a continuación:
118
6. MODELO PARA EL ANÁLISIS DE FLUJO DE CARGA
15
Estos métodos tradicionales de flujo de carga están ampliamente desarrollados en: GRAINGER, John y
STEVENSON, William. Análisis de sistemas de potencia.
119
integran una red de distribución para una representación más exacta del
modelo.
120
Tabla V. Clasificación de los elementos de un sistema de distribución
Switch (combinación de
Carga puntual
interruptor y seccionador)
Cargas distribuidas
121
La formulación de este método surge, básicamente de la aplicaciones de
las leyes de Kirchoff en dos barridos consecutivos a la red, su origen es el
método de red en escalera para el análisis de circuitos eléctricos lineales (por
eso se le llama también ladder iterative technique), aunque debido a la
presencia de cargas de potencia constante y generadores no es posible aplicar
directamente las leyes desarrolladas para sistemas lineales, lo que lo vuelve un
método iterativo para el análisis de redes no lineales, de forma que cada
corriente y voltaje calculado con estas ecuaciones es una aproximación lineal a
los valores reales.
16
SHIRMOHAMMADI, HONG, SEMLYEN y LUO, A compensation based power flow method for weakly
meshed distribution and transmission networks, IEEE Transactions on Power Systems. p. 753-759
122
manera secuencial el análisis de la red, ya que posibilita trabajar por capas,
independientemente de la configuración del sistema.
123
Teniendo el voltaje del alimentador se asume un arranque plano para
todos los nodos del sistema y el algoritmo de solución sigue los pasos
mencionados anteriormente:
(6-1)
(6-2)
124
La matriz de admitancias es una matriz equivalente de las admitancias de
los capacitores y las admitancias de la carga de impedancia constante
según el modelo del capítulo 3.
(6-3)
(6-4)
125
corriente total que entra a cada nodo, siendo esta la suma de las
corrientes de todas las líneas que salen del nodo más el consumo de
corriente es ese nodo (ecuación (6-5)).
(6-5)
Después de calcular las corrientes para todos los nodos del último nivel se
avanza al siguiente nivel (capa 7 en la figura 36), esto se hace con las
matrices generalizadas desarrolladas para todos los elementos en serie
ecuación (6-6), entiéndase elemento en serie transformadores, segmentos
de línea, reguladores de voltaje, switchs y cargas distribuidas. De esta
forma se pueden calcular los voltajes de los nodos siguientes (capa 7 de la
figura 36) y la corriente que sale de cada uno de los nodos de la capa 7
hacia los nodos de la capa 8 (véase figura 36).
Para este caso se usan las matrices generalizadas [a], [b], [c] y [d] que
dependen del elemento en serie que se esté modelando; capítulo 2 para
los segmentos de líneas y capítulo 4 para transformadores y reguladores
126
de voltaje, estas relacionan los voltajes y corrientes del nodo M con los
voltajes y corrientes del nodo N (figura 37) de la siguiente forma:
(6-6)
127
Se calculan nuevamente las corrientes que salen de todos los nodos del
nivel actual (nivel o capa 7 según la figura 36), considerando todos los
componentes en paralelo, y las corrientes que salen de los nodos por
medio de los componentes en serie hacia otros nodos, estos componentes
pueden estar de manera individual o como una combinación de
componentes.
(6-7)
128
Forward sweep: se actualizan los voltajes en un barrido desde el
alimentador hasta los nodos de la capa final. El voltaje se actualiza usando
el calculado para el nodo en la capa superior.
El forward sweep inicia con el voltaje nominal del alimentador y usa las
corrientes calculadas en el backward sweep para actualizar los voltajes de
la red. Esto se hace con las matrices generalizadas [A] y [B] desarrolladas
en los capítulos 2 y 4:
(6-8)
129
Figura 38. Algoritmo de solución de flujo de carga en redes de
distribución
130
Las principales ventajas de este método ante los métodos tradicionales
para el análisis de flujo de carga son:
131
132
7. VALIDACIÓN DE LOS MODELOS
17
El algoritmo para la solución de redes ligeramente malladas se presenta en: LUO Y SEMLYEN, Efficient
load flow for large weakly meshed networks. p. 1309-1316. Y en: SHIRMOHAMMADI, HONG, SEMLYEN y
LUO, A compensation based power flow method for weakly meshed distribution and transmission
networks. p. 753-759.
133
el caso de redes ligeramente malladas que se puede incorporar de una manera
sencilla debido a su similitud con el modelo de generadores.
18
El modelo para transformadores con admitancias en derivación se encuentra en: WANG, Zhuding,
CHEN, Fen y LI, Jingui, Implementing Transformer Nodal Admittance Matrices Into Backward/Forward
Sweep-Based Power Flow Analysis for Unbalanced Radial Distribution Systems. p. 1831-1836.
134
El modelo de generadores permite la incorporación de cualquier cantidad
de generadores, pero se debe tener en cuenta que la eficiencia numérica del
programa se ve afectada directamente por la cantidad de generadores que sean
agregados.
135
7.2.1. Validación del modelo de línea
136
Tabla VI. Comparación de las matrices de impedancias y admitancias
para el caso IEEE34 node test feeder
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
137
Tabla VII. Comparación de las matrices de impedancias y admitancias
para el caso IEEE37 node test feeder
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
138
7.3. Validación del modelo de cargas
Los resultados presentados por el IEEE4 node test feeder muestran dos
tipos de desfasamiento Step Up y Step Down, para las conexiones Y aterrizada-
139
Y aterrizada, y Delta-Delta, pero como se explica en el capítulo 4 estas
conexiones únicamente puede presentar dos desfases: de 120 y -120 grados, lo
cual no cumple con el American Standard Thirty-Degree, en el cual se
especifica que solo las conexiones Y-delta y Delta-Y pueden tener desfases de
múltiplos de 30 grados (llamados comúnmente Step Down y Step Up),
dependiendo de si el voltaje del primario se adelanta al del secundario o
viceversa, pero este desfasamiento no es posible con las conexiones Y
aterrizada-Y aterrizada y Delta- Delta.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
140
Tabla IX. Comparación entre los voltajes y las corrientes para la
conexión Delta-Y aterrizada Step Down Connection
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
141
Tabla XI. Comparación entre los voltajes y las corrientes para la
conexión Y-Delta Step Down Connection
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
142
Tabla XIII. Comparación entre los voltajes y las corrientes para la
conexión Delta-Y Step Up Connection
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
143
Tabla XV. Comparación entre los voltajes y las corrientes para la
conexión Y abierta - Delta abierta Step Up Connection
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
144
Tabla XVI. Comparación entre los voltajes y las corrientes un regulador
de voltaje de tipo B conectado en estrella con una regulación
independiente en las fases
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
En la tabla XVI los datos más importantes son los del número de taps, ya
que son los que determinan cuál será el nivel de la regulación del voltaje de
acuerdo a la ecuación (4-117)
145
fases AB y CB la regulación del voltaje CA se realiza como una combinación de
los otros dos voltajes.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
146
Seguidamente: se calcula la caída de voltaje equivalente entre el
regulador y el centro de carga con los valores de R’ y X’, la corriente de
salida del regulador y la corriente nominal del primario del CT, la corriente
del secundario no es importante para los cálculos, ya que se elimina en el
proceso:
147
valores menores a 0,5 a su valor inmediato inferior. De esta forma el valor
resultante de taps es:
Tabla XVIII. Comparación de los perfiles de voltaje del caso IEEE37 node
test feeder obtenidos con el programa y los presentados en el
caso de prueba
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
149
qué tipo de elemento es, los números enteros son segmentos de línea, donde el
número corresponde a la configuración de línea usada, los números 0,5, 1,5 y
2,5 identifican si el segmento de red tiene un switch, un transformador o un
regulador de voltaje respectivamente.
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
150
Tabla XX. Flujo de potencia activa y reactiva caso IEEE37 node test
feeder
151
Tabla XXI. Consumo de potencia aparente de las cargas caso IEEE37
node test feeder
152
Tabla XXII. Comparación entre las pérdidas halladas con el programa y
las presentadas en el caso de prueba IEEE37 node test feeder
Fuente: IEEE Power Engineering Society Distribution, system Analysis Subcommite, Radial Test
Feeders. Consulta. 15 de septiembre de 2011.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders.html.
153
En las tablas XX y XXI se muestra la comparación entre los flujos de
potencia calculados con el programa y los presentados en el caso de prueba
IEEE 37 node test feeder. En la tabla XXII se muestra la comparación entre las
pérdidas en el sistema halladas con el programa y las presentadas en el IEEE
37 node test feeder.
154
se realizaron dos pruebas al caso IEEE 37 node test feeder, una con un
generador conectado en el nodo 730, y otra con tres generadores conectados
en los nodos 720, 727 y 730.
Los generadores usados son de tipo PV por que el modelo del tipo PQ es
igual al modelo de las cargas de potencia constante y se usó un límite de
potencia reactiva alto en los generadores para demostrar la capacidad
reguladora de estos. La conexión de los generadores es en estrella para regular
los voltajes de línea a neutro. Los generadores son probados con el modelo de
regulación independiente en las fases.
155
Figura 39. Comportamiento de la convergencia en valores por unidad
del voltaje en el nodo 730
0.008
0.006
0.004 Fase A
0.002 Fase B
0 Fase C
-0.002 0 2 4 6 8 10 12 14
-0.004
-0.006
Número de iteraciones
156
Figura 40. Comportamiento de la convergencia en valores por unidad
del voltaje en el nodo 730
0
0 5 10 15 20
-0.01 Fase A
-0.02 Fase B
Fase C
-0.03
-0.04
Número de Iteraciones
0
0 5 10 15 20
-0.01
Fase A
-0.02
Fase B
-0.03 Fase C
-0.04
-0.05
Número de iteraciones
157
Figura 42. Comportamiento de la convergencia en valores por unidad
del voltaje en el nodo 720
0
Diferencia de voltajes Vpu
0 5 10 15 20
-0.005
-0.01 Fase A
-0.015 Fase B
Fase C
-0.02
-0.025
-0.03
Número de iteraciones
158
Figura 43. Comportamiento de la convergencia de voltaje en los nodos
730, 734 y 720 del caso IEEE 37 Node test feeder
159
En la tabla XXIII se muestran los perfiles de voltaje del caso IEEE 37 node
test feeder, para cuatro situaciones diferentes:
160
Tabla XXIII. Comparación de los perfiles de voltaje caso IEEE37 node test
feeder con tres escenarios diferentes
161
En las figuras 44, 45 y 46 se muestran la comparación de los perfiles de
voltaje de cada fase para las situaciones detalladas anteriormente.
1.000
GEN 730
0.990
0.980
0.970 Gen
730,734,7
0.960 20
0.950
0.940
0 10 20 30 40
NODO
162
Figura 45. Comparación de los perfiles de voltaje de la fase B del caso
IEEE37 node test feeder con tres escenarios direferentes
1.030
PERFILES DE VOLTAJE FASE B
Original
1.020
1.010
RV
VOLTAJE Vpu
1.000
Gen 730
0.990
0.980
Gen 720,
730, 734
0.970
0.960
0 10 20 30 40
NODO
163
Figura 46. Comparación de los perfiles de voltaje de la fase C del caso
IEEE37 node test feeder con tres escenarios direferentes
1.020
PERFILES DE VOLTAJE FASE C
1.010
Original
1.000
RV
VOLTAJE Vpu
0.990
Gen 730
0.980
0.950
0 10 20 30 40
NODO
164
CONCLUSIONES
165
6. El método iterativo Forward and Backward sweeps posee muy buenas
características de convergencia, aunque estas se reducen
considerablemente, conforme aumenta el número de generadores en el
sistema.
166
RECOMENDACIONES
167
6. A las personas interesadas en el desarrollo de software para el análisis
de flujo de carga se les propone usar de manera equilibrada la relación
entre exactitud y velocidad de cómputo, considerando la aplicación que
se le piense dar al software.
168
BIBLIOGRAFÍA
169
7. GLOVER, Duncan; SARMA, Mulukutla; OVERBYE, Thomas. Power
system analisis and design. 4a ed. USA: Thomson Learning, 2008.
752 p. ISBN: 0-534-54884-9.
170
14. KHUSHALANI, Sarika, Development of power flow with distributed
generators and reconfiguration for restoration of unbalanced
distribution systems. Director: Noel Schulz. Department of Electrical
and Computer Engineering, Mississippi State University, 2006.
107 p.
16. LUO, G. X; SEMLYEN, A. Efficient load flow for large weakly meshed
networks. Toronto: IEEE Transactions On Power Systems, 1990. p.
1309-1316.
171
20. TENG. Modeling distributed generators in three phase distribution load
flow. Institution of Engineering and Technology Generation,
Transmission and Distribution. 2008, vol. 2, num. 3. p.330-340.
22. WANG, Zhuding; CHEN, Fen; LI, Jingui, Implementing Transformer Nodal
Admittance Matrices Into Backward/Forward Sweep-Based Power
Flow Analysis for Unbalanced Radial Distribution Systems. USA:
IEEE Transactions On Power Systems, 2004. p. 1831-1836.
172