Tomo I
Tomo I
Modelamiento y análisis
A mis hijos Alfonso y Miguel, que son mi alegría, e inspiración, a mi esposa Eliana,
mi complemento, a la memoria de mi padre Alfonso, quien fuera mi guía y amigo,
a la memoria de mi madre Ana Emilia y a mis hermanos.
Ramón.
A mi esposa Suany, a mis hijos Laura, David y Juan Pablo que son mi fuente de
inspiración.
Antonio.
Mauricio.
Contenido
i
1.14 Ecuaciones de potencia inyectada . . . . . . . . . . . . . . . . . . . . . . . . 28
1.15 Flujo de potencia por una línea . . . . . . . . . . . . . . . . . . . . . . . . . . 29
1.15.1 FP en una línea modelada como Zserie . . . . . . . . . . . . . . . . . . 30
1.16 Pérdidas en SEP multinodales . . . . . . . . . . . . . . . . . . . . . . . . . . 31
ii
4 Operación de los sistemas de potencia 111
4.1 Estado de la operación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
4.2 Caracteristicas de la carga . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
4.2.1 Variación de la carga con la tensión y la frecuencia . . . . . . . . . . . 114
4.3 Balance de Potencia Vs frecuencia . . . . . . . . . . . . . . . . . . . . . . . . 115
4.4 Balance de potencia Vs tensión . . . . . . . . . . . . . . . . . . . . . . . . . . 117
iii
iv
Índice de figuras
v
2.19 Curva típica de demanda en Colombia . . . . . . . . . . . . . . . . . . . . . . 56
2.20 Modelo de la línea de transmisión . . . . . . . . . . . . . . . . . . . . . . . . 59
2.21 Red eléctrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
2.22 Representación del transformador con impedancia dividida . . . . . . . . . . . 68
2.23 Impedancia del transformador referida al secundario . . . . . . . . . . . . . . 68
2.24 Impedancia del transformador referida al primario . . . . . . . . . . . . . . . . 69
2.25 Modelo del transformador en p.u. visto en el secundario . . . . . . . . . . . . . 69
2.26 Modelo del transformador en p.u. visto en el primario . . . . . . . . . . . . . . 70
2.27 Modelo del transformador en p.u. con el tap en la posición nominal . . . . . . 70
2.28 Sistema eléctrico básico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
2.29 Modelo en impedancias p.u. del generador-transfomador . . . . . . . . . . . . 73
2.30 Ejercicio de valores p.u. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
2.31 Sistema eléctrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
2.32 Z p.u. en el transformador 3φ . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
2.33 Modelo equivalente del transformador Y − Y . . . . . . . . . . . . . . . . . . 76
2.34 Modelo equivalente del transformador Y − ∆ . . . . . . . . . . . . . . . . . . 76
2.35 Red radial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
2.36 Sistema en malla abierta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
2.37 Sistema en malla cerrada con transformadores en la malla . . . . . . . . . . . . 80
2.38 Sistema en malla cerrada sin transformadores en la malla . . . . . . . . . . . . 81
2.39 Circuito de malla cerrada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
2.40 Circuito de malla cerrada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
2.41 Transformador monofásico con tres arrollamientos . . . . . . . . . . . . . . . 84
2.42 Modelo del transformador de tres arrollamientos . . . . . . . . . . . . . . . . . 84
2.43 Diagrama en p.u. del transformador tridevanado . . . . . . . . . . . . . . . . . 86
2.44 Sistema con transformador tridevanado . . . . . . . . . . . . . . . . . . . . . . 86
2.45 Representación del sistema que contiene un transformador tridevanado . . . . . 87
2.46 Diagrama de una red eléctrica de 7 nodos . . . . . . . . . . . . . . . . . . . . 87
2.47 Admitancia y taps en el secundario . . . . . . . . . . . . . . . . . . . . . . . . 89
2.48 Modelo del transformador con taps . . . . . . . . . . . . . . . . . . . . . . . . 90
2.49 Admitancia al lado opuesto del cambiador de taps . . . . . . . . . . . . . . . . 90
2.50 Modelo del transformador con taps . . . . . . . . . . . . . . . . . . . . . . . . 91
2.51 Modelo del transformador a) . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
2.52 Modelo del transformador b) . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
2.53 Modelo de transformadores en paralelo . . . . . . . . . . . . . . . . . . . . . 93
2.54 Transformador regulante para el control de la magnitud de tensión . . . . . . . 94
2.55 Modelo del transformador desfasador . . . . . . . . . . . . . . . . . . . . . . 95
2.56 Autotransformador elevador . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
2.57 Representación simbólica del autotransformador . . . . . . . . . . . . . . . . . 96
2.58 Ensayo de cortocircuito en el autotransformador . . . . . . . . . . . . . . . . . 97
2.59 Circuito equivalente del autotransformador . . . . . . . . . . . . . . . . . . . . 97
2.60 Representación de la operación del autotransformador . . . . . . . . . . . . . . 98
2.61 Ejemplo propuesto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
2.62 Circuito equivalente del regulador de tensión . . . . . . . . . . . . . . . . . . 101
2.63 Circuito equivalente del regulador de tensión . . . . . . . . . . . . . . . . . . 101
vi
3.1 Red eléctrica de 4 nodos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
3.2 Red eléctrica de 5 nodos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
3.3 Red eléctrica sin ordenamiento nodal . . . . . . . . . . . . . . . . . . . . . . . 108
3.4 Red eléctrica con ordenamiento nodal . . . . . . . . . . . . . . . . . . . . . . 109
vii
viii
Prefacio
La red de transmisión de energía eléctrica es considerada como una de las infraestructuras más
importantes de un país. Su función es esencial para transportar la energía eléctrica desde los
grandes centros de producción a los grandes centros de consumo. La red eléctrica cumple
funciones técnicas, sociales y económicas.
Uno de los problemas más importantes y más complejos en el ámbito de la Ingeniería
Eléctrica es el inherente a la formalización y formulación matemática de los conceptos
asociados a la determinación de las tensiones y los flujos de potencia resultantes en la red
eléctrica de los sistemas eléctricos de potencia. Este problema se conoce en Ingeniería como
el problema de flujo de carga. Su importancia se refleja en el continuo interés por parte de los
investigadores de presentar nuevas formulaciones matemáticas así como de técnicas de solución
para este problema.
El problema de flujo de carga es el aspecto central de los análisis de sistemas eléctricos de
potencia y su formulación y solución es una tarea obligatoria en los estudios de planeamiento
a largo plazo de sistemas eléctricos, en los estudios de seguridad (confiabilidad, contingencias,
cortocircuito, estabilidad) y en los estudios que se realizan de manera permanente para la
programación de la operación diaria de un sistema eléctrico y para su control.
A medida que los sistemas eléctricos evolucionan por causa de las nuevas tecnologías,
este problema debe ser revisado y adecuado a las nuevas circunstancias. Sin embargo, la base
teórica fundamental sigue residiendo en los conceptos y formulaciones del problema tradicional
de flujo de carga.
Esta obra está centrada principalmente en el análisis de los sistemas eléctricos de potencia
(SEP) a partir de la solución del problema de Flujo de Carga (FC) considerando una amplia
variedad de modelos y representaciones de los diferentes componentes eléctricos de la red. La
representación matricial de la red es un aspecto esencial en este texto, de forma que facilite la
implementación computacional de métodos de solución de sistemas de gran tamaño. También
este aspecto es fundamental para el análisis de sistemas con estructuras matemáticas de gran
complejidad.
Generalmente, una obra con el alcance que aquí se expone es presentada en único
tomo. Sin embargo, por restricciones editoriales referentes al número máximo de páginas por
publicación, se ha optado por presentar el contenido completo de la obra en dos tomos. El
primero aborda aspectos sobre el modelamiento del probema de FC tradicional. El segundo
tomo aborda el problema de FC óptimo, el cual proporciona al estudiante algunas herramientas
básicas de optimización matemática aplicables al planeamiento y operación del SEP.
El propósito fundamental de este libro es presentar el modelamiento del problema de FC
y de sus principales fenómenos asociados. También se presentan las bases para desarrollar
análisis de sistemas eléctricos de potencia, así como cubrir las principales técnicas matemáticas
ix
involucradas en la solución de los mismos. El texto está dirigido a estudiantes de pregrado y de
posgrado de Ingeniería Eléctrica, y ciencias afines y ha sido estructurado a partir de las notas
de clase del curso “análisis de sistemas de potencia” ofrecido en el programa de Ingeniería
Eléctrica de la Universidad Tecnológica de Pereira (UTP).
Este libro está organizado en seis capítulos, que comprenden un abordaje del análisis de
sistemas de potencia desde los aspectos básicos e introductorios, pasando por una formulación
tradicional del problema de flujo de carga, hasta una formulación que involucra el flujo de carga
linealizado. El primer capítulo presenta una descripción de los elementos y conceptos básicos
de los Sistemas Eléctricos de Potencia (SEP) que permiten al lector construir una taxonomía
de la red y sus componentes para que, posteriormente en el capítulos dos, se pueda realizar la
representación y el modelamiento matemático de los sistemas a través de sistemas de ecuaciones
que lo representen. El capítulo tres complementa la construcción del modelo matemático a
través de la representación matricial de la red.
Los capítulos 4, 5 y 6 cubren los principales aspectos operativos de los SEP, diferentes
técnicas de solución para el problema de FC no lineal y su versión linealizada. Cada capítulo
incluye ejemplos y ejercicios que ayudan al lector a interiorizar los conceptos que hemos
considerado más relevantes.
Para comprender el contenido de este texto se requiere que el lector conozca los conceptos
fundamentales de la teoría de circuitos, las bases teóricas del algebra lineal y conceptos
asociados a la máquina síncrona operando como generador.
Queremos agradecer a los colegas, estudiantes y profesores, quienes a lo largo del tiempo
han ayudado a enriquecer el contenido de este documento. La principal motivación para
sacar adelante este proyecto es que sirva como punto de referencia y consulta para las nuevas
generaciones de estudiantes y profesionales, quienes a su vez, con toda seguridad, nos ayudarán
y motivarán a incorporar nuevos y mejores contenidos a este trabajo.
x
xi
ELEMENTOS BÁSICOS PARA EL
1
ANÁLISIS DE LOS SISTEMAS DE
POTENCIA
1
Capítulo 1. Elementos básicos de los SEP
• Controlar, supervisar y gobernar en forma particular cada uno de los elementos del
sistema de generación, transmisión, subtransmisión y distribución, y controlar la
operación global del sistema para obtener en cada momento las mejores condiciones,
tanto en estado estable como durante periodos transitorios.
2
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
3
Capítulo 1. Elementos básicos de los SEP
• Un problema local se puede transformar en un problema que afecta la red como un todo,
por ejemplo, problemas de estabilidad y colapsos (blackout).
Efecto En esta franja de tiempo, la energía eléctrica suministrada a la carga adicional es entregada
transitorio por el circuito más próximo a la carga (circuito propio). Este efecto puede ser perceptible
local o no, dependiendo del tamaño de la carga: una bombilla de baja potencia no muestra
un efecto perceptible, un motor de gran potencia puede producir una caída de tensión
observable en los demás dispositivos eléctricos que se encuentran en operación en su
proximidad. En todos los casos aparecerá un efecto transitorio local, así sea pequeño. En
consecuencia, en esta franja de tiempo se puede decir que la nueva carga toma una parte
de la energía almacenada (por ejemplo, en forma magnética) del circuito adyacente.
• Transitorio electromecánico (10−1 s): Después del impacto inicial, existe una respuesta
mecánica del sistema. La energía consumida por la nueva carga es solicitada a los rotores
Pérdida de
energía de las unidades generadoras (y turbinas) que se caracterizan porque poseen elementos
cinética almacenadores de energía (masas propias y volantes de inercia). La energía almacenada
4
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
5
Capítulo 1. Elementos básicos de los SEP
tanto, puede requerirse una última etapa que puede ocurrir con un tiempo de respuesta
considerablemente mayor al de las etapas anteriores. Se trata del redespacho económico
de generación considerando restricciones de seguridad operativa (minimización de
riesgos de colapso). En este proceso de redespacho se consideran los costos de
combustibles de plantas térmicas, las eficiencias de las unidades generadoras y la
disponibilidad de agua de los embalses. El nuevo programa de generación será el que
finalmente definirá el o las unidades generadoras del sistema que atenderán la demanda
de la nueva carga conectada.
Como se puede observar, la energía que requiere la nueva carga es tomada de diferentes
lugares, conforme transcurre el tiempo, y finalmente será un generador o un grupo de
generadores los que asuman el valor de su demanda.
1.4.1 Turbina
Es un dispositivo mecánico en el cual, la energía cinética proveniente de fuentes hidráulicas,
térmicas o eólicas, es transformada en energía mecánica rotacional en su eje.
Gran parte de la energía cinética transportada por un fluido (aire, agua, aceite, etc) es
Energía transferida al eje de una turbina cuando el fluido golpea sus álabes, provocando su rotación. No
cinética
de un
toda la energía cinética del fluido es transformada en energía cinética en el eje de la turbina,
fluido una parte de ella se convierte en calor, debido a la fricción existente entre las partes móviles,
otra parte continua asociada al fluido que sale de la turbina y otra parte se convierte en energía
cinética tranferida a masas de inercia y al aire circundante. Una vez que la turbina ha tomado
una parte importante de la energía del fluido y la ha disponibilizado en su eje, otro eje de acople
unido al eje de la turbina en un extremo y al eje del generador en el otro extremo, transporta
esta energía hasta el eje del generador eléctrico. Las turbinas tienen formas y tamaños variados
dependiendo del fluido que usan, de la cantidad de energía que deben captar del energético
primario, y de la velocidad y la presión del fluido.
1.4.2 Generador
Es una máquina eléctrica constituida por una parte estática y otra rotativa, acopladas entre sí
tanto eléctrica como magnéticamente. El generador transforma la energía mecánica rotacional
entregada por una turbina hidráulica o de vapor, a través de un eje común, en energía eléctrica
disponible en un arreglo de conductores eléctricos, generalmente en configuración trifásica, la
cual es inyectada a una red de transmisión.
Los rotores de los generadores se construyen de polos salientes o polos lisos (cilíndricos),
lo cual implica, para máquinas de igual capacidad, diferentes valores de inercias.
6
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
El generador entrega en sus bornes tanto potencia eléctrica activa como potencia eléctrica
reactiva. La magnitud de ambas potencias puede modificarse entre unos límites mínimos y
Reguladores
máximos. Fundamentalmente el generador entrega una cantidad preestablecida de energía de tensión
activa al sistema y ajusta su potencia reactiva de tal forma que su tensión terminal se encuentre y velocidad
en un cierto valor requerido por el sistema. La estabilización del valor de la tensión terminal se
logra con la acción del denominado regulador de tensión. La estabilización de la frecuencia de
la tensión en 50 Hz o en 60 Hz (según el sistema) se logra a través de la acción del denominado
regulador de velocidad.
Por razones económicas (minimización de pérdidas) la transmisión hacia sitios de
consumo que se encuentran a gran distancia se realiza generalmente a tensiones elevadas
Transformador
(230kV, 345kV, 500kV o 750kV), sin embargo, los generadores no pueden operar con esos elevador
niveles de tensión debido a limitaciones físicas y de nivel de aislamiento eléctrico. Con las
tecnologías convencionales, los generadores operan típicamente en el rango de 10kV a 30kV.
En consecuencia, cuando los generadores se encuentran muy retirados de los centros de carga
resulta necesario utilizar transformadores elevadores, que tienen como función transformar los
niveles de tensión de generación en niveles de tensión de transmisión, con la consecuente
reducción de los niveles de corriente, y por lo tanto, de los niveles de pérdidas óhmicas. Una
innovación tecnológica en desarrollo es el denominado generador-transformador (powerformer,
ABB) que puede producir tensiones a los niveles de transmisión, eliminando así la necesidad
de los transformadores elevadores. Estas máquinas utilizan ranuras profundas en los estatores
y en ellas se insertan conductores convencionales de alta tensión, que permiten la operación sin
causar problemas de aislamiento.
1.4.3 Cargas
Entre todos los componentes del sistema, posiblemente los que ofrecen mayores dificultades
para su modelamiento son las cargas.
Contrario a lo que ocurre con, por ejemplo, generadores y líneas de transmisión, que son
proyectadas y tienen un comportamiento altamente previsible, las cargas representan agregados
Modelos
de consumidores con características muy diversas vistos desde las subestaciones de distribución. de carga.
Además de la diversidad de los elementos que las componen, su variabilidad en el tiempo es un Ver sección
factor adicional que dificulta su modelamiento, aunque a este nivel de tensión esta variabilidad 2.5, pág. 54
resulte relativamente previsible.
La manera más usual de modelar las cargas consiste en representarlas a través de valores
constantes de potencias activas y reactivas (modelo de potencia constante). Otros modelos
posibles son el de corriente constante y el de impedancia constante. Combinaciones ponderadas
de los tres modelos también pueden ser utilizadas.
En la mayoría de los casos prácticos, las cargas son de tipo inductivo, debido a los efectos
de los motores de inducción y de los reactores utilizados en iluminación. Por lo tanto, en
general, las cargas consumen potencia activa y solicitan al sistema potencia reactiva. Una regla
usada en algunas ocasiones, a niveles de transmisión y subtransmisión, sugiere que la potencia
reactiva solicitada por la carga (que a nivel del sistema de potencia es el consumo visto desde
las subestaciones de distribución) no debe superar la mitad del valor de la potencia activa.
7
Capítulo 1. Elementos básicos de los SEP
1.4.4 Subestación
Es un centro de interconexión que permite la conformación de una red eléctrica radial,
enmallada al mismo nivel de tensión o enmallada a diferente nivel de tensión. Las subestaciones
favorecen la distribución de la energía y la modificación de los niveles de tensión transmitidos.
Una subestación puede ser del tipo elevador o reductor.
En las subestaciones elevadoras, la potencia entregada por los generadores a tensiones
medias entre (10 - 30 kV) son transformadas a niveles de tensión muy superiores para permitir su
transmisión económica a las áreas de consumo. En las subestaciones reductoras (115 - 500 kV)
se reciben las potencias provenientes de uno o varios centros de generación, para que una vez
transformados a nivel de tensión apropiado puedan distribuirse a los usuarios o sub-transmitirse
a otras regiones.
Componentes de las subestaciones.
• Transformador de potencia.
Permite convertir niveles de tensión altos en bajos y viceversa. Además y según el
diseño cumplen funciones de regulación bien sea de la magnitud de la tensión primaria o
secundaria, afectando el flujo de potencia reactiva o bien, regulando los flujos de potencia
activa y reactiva. La facilidad de alterar los niveles de tensión a través de transformadores
es uno de los mayores atractivos del uso de la corriente alterna y esto justifica su amplia
utilización.
Los transformadores de potencia se caracterizan por presentar pérdidas de potencia
activa y reactiva, y por poseer un valor de reactancia inductiva generalmente superior
al valor de su resistencia (típicamente, el valor de la reactancia es 20 a 50 veces el
valor de su resistencia). Las excepciones, en general, se presentan en los modelos de
transformadores tridevanados, en los cuales pueden aparecer valores o relaciones poco
usuales. (pueden aparecer, por ejemplo, reactancias negativas en estos modelos). De otro
lado, los transformadores reguladores permiten la modificación de su tap, lo que afecta
la relación de transformación y su impedancia. En los transformadores más comunes, los
denominados transformadores en fase, el tap es una cantidad real (escalar). El efecto del
tap de un transformador depende del tipo de elementos que se encuentran conectados a
lado y lado de este elemento. Si el transformador está conectado entre dos nodos de una
red radial y en uno de ellos existe un generador con regulador de tensión y en el otro una
carga, al alterar la posición del tap, la magnitud de la tensión en el lado del generador
permanecerá prácticamente inalterada, mientras que la magnitud de la tensión en el lado
de la carga va a presentar un cambio. Cuando el transformador en fase está conectado a
una red de tal forma que hace parte de una malla o trayectoria cerrada, la modificación
del tap afectará generalmente la distribución del flujo de potencia reactiva en la malla.
• Interruptores de potencia.
En condiciones normales de trabajo, estos elementos permiten abrir o cerrar un circuito
eléctrico para interrumpir cuando sea conveniente o necesario el flujo de potencia
en determinadas partes del sistema. En condiciones anormales de trabajo abren
automáticamente para valores determinados de sobrecarga o cortocircuito aislando el
elemento de funcionamiento anormal.
8
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Los seccionadores.
Aislan partes vivas del sistema eléctrico de partes que temporalmente deben estar fuera
de servicio. Los seccionadores no se diseñan para interrumpir corrientes de falla; ellos se
pueden maniobrar sólo cuando los circuitos donde están actuando no presentan paso de
corriente.
Operación
• Pararrayos. en vacío
Controlan sobretensiones de carácter atmosférico o producidas por maniobras (apertura
o cierre de circuitos energizados), evitando la propagación de oscilaciones que ponen en
peligro la vida de los equipos.
Descargador de
• Transformadores de potencial y de corriente.
sobretensión
Permiten realizar mediciones de los valores eficaces o magnitudes fasoriales de los DSP
parámetros correspondientes a la tensión y a la corriente y sirven además para emprender
tareas de protección al alimentar los relés.
TP’s y
• Barrajes.
TC’s
El barraje es la parte física del sistema que se asocia a los nodos de la red. Las barras
o barrajes, en realidad, son conductores con resistencia eléctrica despreciable, por lo
menos cuando se les compara con las impedancias de las líneas y transformadores, lo
que justifica su representación como nodos eléctricos. Las barras están localizadas, en
general, en las denominadas subestaciones del sistema de potencia, y están constituidas
por interconexiones de de barras conductoras conectadas entre sí a través de interruptores
o seccionadores. En cada subestación pueden existir diferentes niveles de tensión
(generalmente 2 o 3) según la salida de los transformadores. A cada nivel de tensión se
asocia un barraje, al cual confluyen (o de donde parten) diferentes alimentadores trayendo
(llevando) potencia de (hacia) otros centros de generación o consumo. Para garantizar la
continuidad de la operación de la subestación, el barraje crespondiente a cada nivel de
tensión puede duplicarse o si se requiere puede seccionarse en tramos.
9
Capítulo 1. Elementos básicos de los SEP
10
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
electricidad, los despachos de generación diaria se realizan basados en ofertas económicas que
realizan varios agentes pertenecientes al sistema, con capacidad de generación, interesados en
salir programados en dichos despachos. Normalmente las ofertas son ordenadas de menor a
mayor precio y son programadas en orden creciente de precio hasta atender la demanda del
sistema. Lógicamente la capacidad de generación efectiva del sistema siempre es mayor que la
demanda.
Los sistemas eléctricos han evolucionado desde los esquemas verticalmente integrados
hacia los esquemas con mercado de electricidad, ya que los esquemas sin mercado normalmente
involucran inversiones innecesarias y riesgos fiscales elevados para los aparatos económicos de
los países. Los esquemas con mercado permiten involucrar aspectos de competitividad en el
sector eléctrico y su objetivo es mitigar el riesgo proyectado a los consumidores finales tanto
desde el punto de vista de la continuidad del servicio como de los precios que se debe pagar por
la energía eléctrica.
A diferencia de los esquemas verticales, los esquemas con mercado de electricidad tienen
agentes propietarios del sistema de generación, transmisión y distribución, los cuales realizan
inversiones que deben ser recuperadas en el largo plazo a través de la operación diaria de
sus elementos. Otro aspecto importante es que los agentes compiten entre sí a través de sus
ofertas y el precio de la energía eléctrica resulta de la interacción oferta-demanda, emulando
un mercado tradicional. Desafortunadamente las redes eléctricas existentes no generan siempre
condiciones de acceso abierto para todos los competidores, es decir, un agente con una buena
oferta económica puede no salir despachado por causa de las restricciones del sistema de
transmisión.
En general, los esquemas con mercado pueden no resultar muy eficientes si no se realiza
una adecuación de la red de transmisión y de la regulación existente, con reglas claras y estables
en el tiempo. De otro lado, un esquema de mercado exitoso en un país puede resultar ineficiente
en otro país, aunque tengan características similares. Para que un esquema desregulado sea
eficiente debe ser diseñado para las características propias de dicho sistema. Este aspecto
da origen a una línea muy interesante de investigación denominada arquitectura de mercados
eléctricos.
El tema de mercados eléctricos involucra aspectos técnicos muy importantes que pasan
obligatoriamente por los análisis de sistemas de potencia y la solución de problemas de flujo de
carga, aspecto que será tratado con profundidad en este texto.
11
Capítulo 1. Elementos básicos de los SEP
energía primaría son alteradas y convertidas en energía eléctrica que puede ser transportada y
suministrada a quienes demandan este tipo de energía.
En los sistemas eléctricos de potencia tradicionalmente se ha hecho uso de las energías
primarias provenientes de los recursos fósiles, también denominadas energías no renovables, y
se ha utilizado la energía cinética presente en el agua que fluye de manera natural por los ríos o
las corrientes marinas, o que fluye de manera controlada por canales o ductos construidos por
el hombre. La energía primaria proveniente del agua fluyente hace parte de las denominadas
energías renovables.
En un sistema eléctrico de potencia, las denominadas Plantas de generación térmicas
e hidráulicas son las encargadas de transformar grandes cantidades de energía primaria
proveniente del gas, del carbón o del agua fluyente en energía eléctrica, suficiente para abastecer
grandes conglomerados de consumidores agrupados en poblaciones o para atender la demanda
de grandes industrias.
12
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Gas natural
• Carbón
• Biomasa
• Basuras o desperdicios
• Alcoholes combustibles
• Sol
• Energía Solar
• Energía Geotérmica
• Energía de Biomasa
• Energía Eólica
13
Capítulo 1. Elementos básicos de los SEP
• Energía Maremotriz
• Energía de aceites vegetales combustibles
• Energía de alcoholes combustibles
G C
Generador Transformador
Líneas de
transmisión
Nodo de
generación Nodo de
carga
14
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
SG G
C
S C-K
S DC
Figura 1.4. Representación de un nodo de carga
15
Capítulo 1. Elementos básicos de los SEP
1.7.3 Topología
Se define así al estado de conexión momentánea de los elementos físicos del sistema que
integran una “configuración” o estructura física.
n
X
• Balance nodal: SG i − SD i = Sij ; i = 1, 2, . . . k
j=1
m
X R
X
• Balance global: SG i − SD j − SL = 0
i=1 j=1
G1 T1 LT1 T2
LT2 D2
16
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
G1 T1 T2 G2
SD1 SD2
Doble
circuito
SD
17
Capítulo 1. Elementos básicos de los SEP
SG1 S D1
SG 2 SD2
Red Interconectada
M M
S Gm S Dn
• Uso racional a lo largo del tiempo de toda la energía primaria existente en el sistema
integrado.
18
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
G1 T1 T2 G2
Línea de transmisión
Área 1 Área 2
diferentes y mejores estrategias de coordinación entre los ORs para lograr una operación segura
y económica del sistema global. Algunos factores básicos pueden dificultar estas labores
de coordinación. Por ejemplo, para obtener el punto de operación de un área particular,
generalmente, el OR debe conocer todos los datos de red del sistema global. Actualmente,
con la existencia de mercados competitivos, no siempre es posible disponer de esta información
debido a políticas de confidencialidad propias de cada OR. Por otro lado, tener acceso a esta
información implica manejar grandes volúmenes de información y consecuentemente lidiar con
un mayor número de problemas técnicos de comunicación y con tiempos computacionales que
podrían llegar a ser prohibitivos.
El concepto de descomposición en regiones o esquema descentralizado surge como
una solución que proporciona flexibilidad, velocidad de cómputo y robustez matemática a los
Esquema
procesos de coordinación en un Sistema Eléctrico de Potencia (SEP) multi-áreas. Este concepto descentra-
permite una operación independiente de cada OR sujeta a un proceso de coordinación global lizado
que busca una operación óptima integrada. Los SEP multi-área representan adecuadamente
lo que sucede en los mercados eléctricos desregulados, en donde la operación del sistema
es fuertemente afectada por la competencia de los participantes del mercado y donde la
información de la red es clasificada. La idea básica consiste en dividir el problema global en
sub-problemas regionales asociados a cada área en donde cada OR resuelve su propio problema
de flujo de potencia óptimo (FPO). Posteriormente, un proceso iterativo de coordinación global
intercambia alguna información estratégica entre las áreas a fin de alcanzar la misma solución
de un FPO tradicional centralizado.
Las principales razones para adoptar el esquema descentralizado son:
19
Capítulo 1. Elementos básicos de los SEP
• Es posible considerar situaciones reales en que las diferentes áreas de un SEP manejan
diferentes tipos de restricciones operativas. En este esquema cada OR resuelve su propio
problema de flujo de potencia.
S Gi
S i1
S Di Si 2
M
S ik
Figura 1.12. Representación de la potencia neta inyectada en un nodo
20
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Para el nodo i en cuestión se define la potencia neta inyectada del sistema externo a la
diferencia entre la potencia generada en ese nodo y la potencia allí demandada.
En la formulación j y k sólo incluyen los nodos que tienen unión física con el nodo i.
Así, en el nodo i debe existir un balance de potencia del sistema interno y externo, de la
siguiente forma:
Balance
nodal
l
X t
X
SGi − SDi = Sij − Ski (1.3)
j=1 k=1
l
X t
X
SGi − SDi − Sij + Ski = 0 (1.4)
j=1 k=1
SGi − SDi = SN i
⇒ (PGi − PDi ) + j (QGi − QDi ) = PN i + jQN i (1.5)
l
X t
X
PGi − PDi − Pij − Pki = 0 (1.6)
j=1 k=1
l
X t
X
QGi − QDi − Qij − Qki = 0 (1.7)
j=1 k=1
En el nodo i puede darse el caso particular de no tener una de las componentes básicas:
generación o carga, o ninguna de éstas.
21
Capítulo 1. Elementos básicos de los SEP
• Variables no controladas o de perturbación PDi , QDi son variables con las cuales no es
posible establecer control y presentan un comportamiento impredecible a nivel de usuario
final, sin embargo a nivel de transmisión su comportamiento sí es predecible. Para su
estimación se utilizan diferentes técnicas de proyección que dependen del horizonte de
tiempo (corto, mediano o largo plazo).
• Variables de control PGi , QGi son variables que permiten operar y controlar el sistema
eléctrico. Su rango de operación está fijado entre unos límites mínimos y máximos.
Los valores de operación son determinados para garantizar la calidad, seguridad y
confiabilidad del sistema eléctrico.
• Variables de estado | Vi |, θi
Estas últimas son aquellas variables dependientes (afectadas por las variables de control
y de perturbación) con cuyo conocimiento se puede estimar el estado del sistema, es decir las
condiciones de operación en un determinado momento: flujos de potencia, corrientes a través
de las líneas, pérdidas, etc.
22
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
En consecuencia, por cada nodo sólo podrán establecerse dos ecuaciones, una para la
potencia activa y otra para potencia reactiva:
l
X t
X
PN i = PGi − PDi = fi (Vi , Y ) = Pij − Pki (1.8)
j=1 k=1
l
X t
X
QN i = QGi − QDi = fi (Vi , Y ) = Qij − Qki (1.9)
j=1 k=1
donde: l = Nodos a los cuales llega el flujo. t = Nodos de los cuales sale el flujo.
En las ecuaciones anteriores se observa que son formulaciones escritas en función de las
variables nodales y parámetros de red: PGi , QGi , PDi , QDi , | Vi |, θi .
Por tanto, por cada nodo se establecen:
• Nodo donde se especifica PGi y QGi , denominado nodo (P, Q) o de potencia neta
inyectada.
23
Capítulo 1. Elementos básicos de los SEP
• Límites de tensión:
Vmin ≤ V ≤ Vmax
24
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Para realizar analisis del sistema eléctrico se han planteado modelos matemáticos en
función de variables y parámetros nodales. En un principio se propone un modelo lineal, cuyas
variables son la corriente nodal inyectada y la tensión nodal.
I i1
i 2
Ii2
M
I ii I in
n
Combinando la ecuación 1.10 primera ley de Kirchhoff y la ecuación 1.11 segunda ley
de Kirchhoff y donde yik corresponde a la admitancia serie de la línea entre los nodos ik y yio
como la sumatoria de todas las admitancias que conectan al nodo i con la referencia, la cual
incluye efectos shunt de las líneas, condensadores y reactores, se obtiene:
25
Capítulo 1. Elementos básicos de los SEP
yio Vo + yi1 V1 ... + yin Vn − (yio + yi1 ... + yin )Vi + Iii = 0
−yio Vo − yi1 V1 ... + (yio + yi1 .. + yin )Vi ... − yin Vn = Iii
n
X
Siendo el valor de los términos: Yii = yij y Yij = −yij .
j=0
n
X
La ecuación anterior representada en forma algebraica es: yik Vk = Iik
k=0
26
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
n
X
– Elementos de la diagonal: Ykk = yik
i=0
– Elementos fuera de la diagonal: Yki = −yki i: Nodo con el cual el nodo k tiene
conexión física.
I N = YN · V N
YN−1 IN = YN−1 · YN · VN
VN = YN−1 · IN
ZN = YN−1
27
Capítulo 1. Elementos básicos de los SEP
V1 Z11 · · · · · · Z1n I11
.. ··· ··· .
..
.
.. = ··· ··· .
. ..
··· ···
Vn Zn1 · · · · · · Znn Inn
La notación matricial es la siguiente: Vbus = Zbus · Ibus
1
Yik = −yik =
zik
Esto es verdad para la Ybus siempre y cuando no existan acoples mutuos entre las líneas.
En la matriz admitancia nodal, el elemento (i − k) corresponde a la admitancia física entre
dichos nodos.
1
Zik = zik = falso para (Zbus )
yik
En la matriz impedancia nodal, el elemento (i-k) corresponde a la impedancia eléctrica
(de transferencia) entre dichos nodos.
PGi + jQGi
1
PDi + jQDi 2
M
n
Figura 1.14. Representación del balance de potencia nodal
28
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
SN i = SGi − SDi
Potencia neta inyectada en el sistema interno.
SN i = Vi Iii∗
La corriente inyectada en el nodo i es definida como:
n
X
Iii = Yik Vk
K=0
n
!
X
QN i = QGi − QDi = Imag Vi Yik∗ Vk∗
K=0
Nota: Estas ecuaciones son expresadas en componentes polar o rectangular y son del tipo
algebraicas no lineales. Para su solución se requieren métodos numéricos especiales.
∗
Sik = Pik + JQik = Vi Iik
Cuando el flujo de potencia se calcula en el sentido de nodo k al nodo i se tiene:
∗
Ski = Pki + JQki = Vk Iki
Las pérdidas de potencia en la línea se calculan como la suma de las potencias vistas en
los nodos i y k.
29
Capítulo 1. Elementos básicos de los SEP
k l
Pkl Plk
Qkl Qlk
Figura 1.15. Flujo de potencia que fluye por la sección serie en el modelo de la línea
1.15.1 Cálculo del flujo de potencia por una línea representada por su
impedancia serie
El modelo matemático de la línea de transmisión utilizado para el estudio de flujo de potencia
es el equivalente pi (Π) de parámetros concentrados, constituido por una impedancia serie y
dos elementos paralelos conectados en los extremos de la misma. El primer modelo estudiado
considera la línea corta, en este se desprecian los elementos paralelos y la línea se representa a
través de los parámetros serie.
Representación de los parámetros serie como impedancia: zkl = rkl + jxkl (Resistencia y
reactancia).
Representación de los parámetros serie como admitancia: ykl = gkl + jbkl (Conductancia
y susceptancia).
Para efectos del modelo matemático la línea es representada a través de la conductancia y
susceptancia serie, su cálculo en función de los parámetros impedancia es el siguiente:
rkl
gkl = 2
rkl + x2kl
xkl
bkl = − 2
rkl + x2kl
k l
Rkl + jX kl
Vk Vl
I kl
Las tensiones nodales en los nodos que conectan la línea, medidos respecto al sistema de
referencia se escriben como:
30
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
serie
IkL = IkL = ykL (V̄k − V̄L )
Efectuando reemplazos en la admitancia y tensiones nodales se obtiene:
∗
SkL = V̄k∗ IkL = PkL − JQkL
PkL − JQkL = Vk e−Jθk (gkL + JbkL )(Vk eJθk − VL eJθL )
El flujo de potencia entre k y l expresado en su componente real (P) e imaginaria (Q) es
el siguiente:
31
Capítulo 1. Elementos básicos de los SEP
La pérdida técnica total es igual a la suma de las pérdidas individuales que se producen
en los elementos del sistema. Este valor también puede ser calculado realizando el balance de
potencia. La potencia de pérdidas resulta entonces igual a la potencia total generada menos la
potencia total demandada:
m
X R
X
SL = SGi − SDj
i=1 j=1
n
X n
X
SL = SN i = Vi Iii∗
i=1 i=1
SL = VNT IN∗
Estas pérdidas pueden ser expresadas en función de la matriz impedancia de barra y de
las corrientes nodales o en función de la matriz admitancia de barra y de las tensiones nodales.
SL = INT ZN IN∗
Ejercicio
Dado el siguiente sistema de cinco nodos y seis líneas, y los valores de sus tensiones
nodales, se pide determinar:
• La potencia nodal.
32
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
2 5
1 4
Con el fin de verificar los valores obtenidos, se presentan las potencias inyectadas en los
nodos:
SG1 = 303 MW ; 159 MVAr
SD2 = 200 MW ; 100 MVAr
SD3 = 150 MW ; 70 MVAr
SG4 = 160 MW ; 153 MVAr
SD5 = 90 MW ; 60 MVAr
33
Capítulo 1. Elementos básicos de los SEP
34
REPRESENTACIÓN DE LOS
2
SISTEMAS ELÉCTRICOS
35
Capítulo 2. Representación de los SEP
Υ
Υ ∆ Υ Υ ∆
Figura 2.1. Representación de un sistema usando diagrama unifilar
36
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
respecto a unos valores base. En este último caso se dice que el sistema se encuentra en por
unidad.
Interruptor de potencia
Fusible
Barra
Transformador dos
devanados
Transformador tres
devanados
Seccionador
Conexión de un transformador de
tres devanados ∆ Υ Υ
Línea de transmisión
Carga estática
Capacitor
Reactor inductivo
37
Capítulo 2. Representación de los SEP
comportamiento del sistema completo. El modelaje es una de las áreas del análisis de sistemas
eléctricos y se define como la representación matemática de los elementos de la red. La
complejidad o simplicidad del modelo depende del grado de exactitud requerido, así como
el tipo de estudio a realizar. En este capítulo se presentan los modelos que se utilizan para
representar los elementos básicos del sistema de potencia en estudios de estado estacionario.
Los compensadores capacitivos e inductivos se representan por su valor de reactancia o por el
valor de la potencia reactiva que suministran o consumen a la tensión de alimentación. Para
los diferentes elementos del sistema de potencia existen modelos más detallados utilizados en
otra clase de estudios, como por ejemplo los modelos de los elementos usados en estudios
de distorsión armónica, los cuales pueden ser consultados en los textos especializados de las
diferentes temáticas.
Xd"
PG Q G VG
38
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
puede entregar potencia reactiva (condensador síncrono) o absorber potencia reactiva (reactor
síncrono).
En cualquier caso, dado que la máquina es un elemento físico real del sistema de potencia,
se debe tener en cuenta que su capacidad es límitada y que sus potencias activa Pk y reactiva
Qk deben satisfacer las siguientes restricciones:
Si los límites Pkmin , Pkmax , Qkmin y Qkmax de las restricciones fueran fijos e
independientes entre sí, estos valores permitirían establecer una región de operación factible
Límites
para el generador que puede representarse geométricamente como una región rectangular que ideales de
encierra todas las combinaciones de valores de potencia activa y de potencia reactiva que pueden generación
ser programadas durante la operación y que no violan los límites superiores ni inferiores de
potencia activa y reactiva. Los puntos que pertenecen al borde del rectángulo también son
factibles porque las restricciones incluyen la igualdad. La figura 2.5 muestra geométricamente
esta región factible. Se puede observar que la potencia activa que el generador puede entregar
se encuentra en la parte positiva, es decir, no se considera que un generador absorba potencia
activa del sistema en condiciones normales de operación. También se observa que la potencia
reactiva límite que puede suministrar es mayor que la reactiva límite que puede absorber. De
manera aproximada, Qkmin es 0,4 veces Qkmax y de signo opuesto.
Q (MVAR)
P P
P(MW )
39
Capítulo 2. Representación de los SEP
Q (MVAR)
Región factible
P(MW )
En los problemas de flujo de carga normalmente se especifican las tensiones que los
generadores deben suministrar en sus bornes. Para cada valor de tensión, que es posible
establecer en terminales de un generador síncrono, existe un valor de potencia reactiva asociada
al generador. En consecuencia, una vez establecido el valor deseado de tensión en bornes
del generador, es necesario calcular el valor asociado de potencia reactiva y verificar si el
punto (Pk ,Qk ) pertenece a la región factible establecida a través de la curva de capacidad
40
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
del generador. Esto quiere decir que se deben cumplir los límites mínimos y máximos de
generación de potencia reactiva asociados al nivel de generación actual de potencia activa.
Puede verificarse en la figura 2.6 que para un cierto valor de potencia activa programada, varían
los límites mínimos y máximos de generación de potencia reactiva.
A continuación se presenta la forma de establecer los límites de la curva de capacidad de
un generador síncrono y su significado, para condiciones de operación en régimen de estado
estacionario.
Estator Estator
A A
C C
B B
Rotor Rotor
B C B C
A A
41
Capítulo 2. Representación de los SEP
Snom
Iamax = √ (2.4)
3 × Vnom
También la potencia aparente S puede expresarse en función de la potencia activa P y
reactiva Q del generador a través de la expresión:
p
S= P 2 + Q2 (2.5)
Al hacer una representación geométrica de los puntos de operación (P ,Q) que no exceden
la corriente de armadura límite de un generador particular se obtienen los puntos encerrados
en una circunferencia de radio Snom centrada en el origen de los ejes P y Q, figura 2.8, y que
satisfacen la desigualdad:
p
P 2 + Q2 ≤ Snom (2.6)
La potencia activa que la máquina le entrega al sistema de potencia está dada por:
√
P = 3Vt × Ia cosθ
42
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Q (MVAR)
Qk
S
Pk
P(MW )
Snom (MVA)
Aunque los generadores síncronos pueden operar temporalmente como motores, por
ejemplo durante contingencias, lo normal es que exclusivamente operen como generadores.
A aquellas máquinas que no se les permite la operación como motor en ningún instante del
tiempo se les agrega la restricción: P ≥ 0. La operación exclusiva como generador (potencia
activa entregada) modifica la región factible de la curva de capacidad, como se muestra en la
figura 2.9.
43
Capítulo 2. Representación de los SEP
Q (MVAR)
S (MVA)
P(MW )
Región factible
como un generador que suministra potencia reactiva al sistema, el producto Eg .cosδ debe ser
mayor a la magnitud de la tensión en terminales Vt y el ángulo de la tensión inducida (δ),
debe ser también mayor que el ángulo de la tensión en terminales. La primera condición:
|Eg .cosδ| > |Vt | garantiza que el flujo de potencia reactiva va desde el generador hacia el
sistema. La segunda condición: δ > 0 garantiza que la potencia activa fluye desde el generador
hacia el sistema. Para que la operación de la máquina sea estable, además se debe garantizar que
el ángulo δ sea menor de 90o eléctricos. En consecuencia, el diagrama fasorial para la operación
como generador sobrexcitado sigue la forma presentada en la figura 2.11.
P
Ia Q
j Xs
If
Eg Vt
ω
Figura 2.10. Circuito equivalente del generador síncrono
44
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Eg∠δ
Eg∠δ
jXs Ia
δ
Vt ∠ 0 o
δ
Vt ∠ 0 o
a) b) Ia
∗
o (Eg ∠δ) − Vt
S = (Vt ∠0 = Vt )Ia∗ (2.7)
jXs
2
j
S = Vt · (Eg ∠− δ) − Vt · (2.8)
Xs
Dado que la potencia activa suministrada por el generador corresponde a la parte real de
la expresión anterior y la potencia reactiva corresponde a la parte imaginaria, estas potencias
pueden calcularse a través de:
Vt · E g Vt · E g
P =− · sen(−δ) = · senδ (2.9)
Xs Xs
2
Vt · E g Vt
Q= · cosδ − (2.10)
Xs Xs
Reescribiendo estas expresiones se tiene:
)
V ·E
P = tXs g · Sen(δ) 2 2
Vt2
2
Vt Vt ·Eg
⇒ VtX·Es g = P 2 + Q + Xs
(2.11)
Q + Xs = Xs · Cos(δ)
2 2
V2 Vt · E g
P + Q+ t
2
= (2.12)
Xs Xs
Se puede observar que los valores de P y Q que permiten la operación del generador
sin producir sobrecalentamiento del devanado de excitación corresponden a todos los puntos
V2
pertenecientes al interior de una circunferencia de radio VtX·Es g con centro en el punto (0, − Xts ).
La figura 2.12 muestra esta región.
Al superponer la región factible establecida a partir de la corriente de armadura máxima,
figura 2.9, y la región factible obtenida con la condición de corriente de excitación máxima,
figura 2.12, se obtiene la región resultante mostrada en la figura 2.13.
45
Capítulo 2. Representación de los SEP
Vt E g
Xs
δ P
2
Vt
-X
s
En la figura 2.13 los puntos dentro de la circunferencia de radio Snom centrada en el origen
(0, 0) representan la región de valores de P y Q que pueden programarse en el generador y que
no originan corrientes de armadura superiores al valor máximo permitido. Los puntos dentro de
V2
la circunferencia de radio mayor VtX·Es g y con centro en (0, − Xts ) representan los valores de P y Q
que al programarse no llevan a la corriente de excitación mas allá de su límite permitido. Puede
observarse que existen valores de potencia activa y reactiva que no producen sobrecalentamiento
del devanado de armadura pero que requerirían de corrientes de excitación superiores a su valor
máximo. Estos son los puntos del círculo de radio Snom que no pertenecen al círculo de radio
Vt ·Eg
Xs
. La región sombreada es la intersección de los círculos y representa al conjunto de valores
de P y de Q que simuntáneamente cumplen con el valor límite de la corriente de armadura y de
la corriente de campo.
46
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Q (MVAR)
Límite de corriente
de armadura
Límite de corriente
Snom de excitación
P(MW )
Región factible
2
Vt
-
Xs
Vt E g
Xs
47
Capítulo 2. Representación de los SEP
Eg∠δ Vt ∠ 0 o
Vt ∠ 0 o
δ
δ
δ
Vt ∠ 0 o Eg∠δ
a) b) Eg∠δ c)
48
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
de agua disponible, una reducción del caudal produce una reducción de la energía mecánica y
por lo tanto de la energía eléctrica que podrá ser obtenida. Este límite aparece en la curva de
capacidad del generador como una línea recta perpendicular al eje de potencia activa y que pasa
por el valor de potencia máxima. Ver figura 2.15.
49
Capítulo 2. Representación de los SEP
Q (MVAR)
P(MW )
Un generador sincrónico debe satisfacer sus especificaciones de placa entre el 95% y el 105% de
su tensión nominal (Norma ANSI C50,12-1982). La figura 2.16 muestra la porción sobrexcitada
del diagrama de cargabilidad de un generador síncrono de 65,78 MVA, 13,8 kV y factor de
potencia 0,95, cuando se varía la tensión entre 12,42 kV (90% del valor nominal) y 14,49 kV
(105% del nominal). Puede observarse que la potencia reactiva básicamente queda limitada por
el calentamiento del campo: a mayor tensión terminal el límite de potencia reactiva disminuye.
Cuantitativamente se observa que el incremento del 5% en la tensión terminal reduce en
cerca del 14% el límite de entrega de potencia reactiva. La reducción del 10% en la tensión
incrementa en más del 20% la capacidad de reactiva. En la zona restringida por la corriente de
armadura se aprecia que la potencia disminuye con la tensión, sin embargo, la limitación se ve
50
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
MVAR
MW
Figura 2.16. Efecto de la tensión en la curva de capacidad de la máquina síncrona
La presión del gas utilizado como refrigerante (generalmente hidrógeno) tiene marcada
influencia sobre los límites de generación de potencia reactiva. La figura 2.17 muestra las
curvas de capacidad de un turboalternador con especificaciones de placa: 83,4 MVA, factor de
potencia 0,85 y una presión de 1,05 bares (la curva más interna).
En esta figura se puede observar que al aumentar en dos atmósferas la presión del gas
refrigerante se puede aumentar en un 17% la capacidad de potencia reactiva.
La temperatura del agua o del aire de refrigeración en los conductos de entrada también es
utilizada por los fabricantes como variable para establecer los márgenes de generación de
potencia reactiva. Las curvas son similares a las presentadas en la figura 2.17 con la diferencia
que en este caso dependen de la temperatura en grados centígrados. Al aumentar la temperatura
se reduce la región factible de la misma manera que ocurre cuando la presión del gas refrigerante
disminuye.
51
Capítulo 2. Representación de los SEP
40
20
0 20 40 60 80
MW
Figura 2.17. Efecto de la presión del gas refrigerante en la capacidad de la máquina síncrona
(RP + RS ) J ( X P + XS )
52
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Las líneas medias adoptan generalmente el modelo Π o T. Las longitudes de una línea
de transmisión media en relación al nivel de tensión se presentan en la tabla 2.2.
Modelo de
línea media
Tensión de la Línea (kV) Longitud Máxima (km)
V < 115 80 ≤ L ≤ 200
115 ≤ V < 500 40 ≤ L ≤ 200
V ≥ 500 20 ≤ L ≤ 100
En las líneas medias el modelo circuital considera el efecto capacitivo. Para el modelo Π
la resistencia total se conecta en serie con la reactancia total y a lado y lado de esta impedancia
serie se conecta una admitancia paralelo a tierra que modela el efecto capacitivo de la línea.
53
Capítulo 2. Representación de los SEP
Las líneas largas pueden adoptar un modelo Π idéntico al de las líneas medias,
corrigiendo los valores de zserie y yshunt con las siguientes expresiones:
Modelo de
línea larga
senh (γ · l)
z =z· (2.14)
γ·l
tangh γ·l
2
y =y· γ·l
(2.15)
2
donde:
l: longitud de la línea.
√
γ: constante de propagación γ = y·z .
y: admitancia shunt por unidad de longitud.
z: impedancia serie por unidad de longitud.
Desde el punto de vista de los estudios de flujo de carga, el modelo adoptado es el de línea
corta y línea media sin efectuar correcciones a la impedancia serie y admitancia paralelo. Desde
el punto de vista de cortocircuito, se efectúan algunas simplificaciones, así: se desprecian las
admitancias paralelo y se desprecia la resistencia serie de las líneas de alta tensión. En sistemas
de distribución el valor de resistencia es significativo frente a la reactancia y por lo tanto debe
ser considerado.
Las líneas de transmisión conducen la energía generada en las centrales eléctricas hasta
los centros de consumo, usando niveles de alta tensión del orden de los cientos de Kilovoltios.
En general, la transmisión de energía se realiza utilizando sistemas de corriente alterna
trifásica, que se representan usando los modelos presentados anteriormente, pero también puede
realizarse usando sistemas de corriente continua. Los sistemas de transmisión de alta tensión
en corriente continua: HVDC, desempeñan un papel importante cuando se utilizan de manera
complementaria con los sistemas de corriente alterna. Para largas distancias, la transmisión
HVDC se convierte en una alternativa interesante. Los circuitos de corriente continua ofrecen
mejores posibilidades para controlar los flujos de potencia, lo que no ocurre con los sistemas
de corriente alterna. Dado que la generación y la distribución se realiza en corriente alterna, los
sistemas de corriente continua requieren de rectificadores e inversores para su interconexión con
el sistema de corriente alterna. Históricamente, los costos asociados a estos equipos han sido la
principal limitante de la utilización de los sistemas de HVDC. En este texto sólo se considera la
transmisión usando corriente alterna trifásica.
54
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Domiciliarias o residenciales.
• Comerciales.
• Industriales.
• Son altamente previsibles en magnitud, para cada franja horaria (para niveles de
transmisión y subtransmisión).
• La variación de la carga es generalmente muy lenta si se les compara con las constantes
de tiempo del sistema de energía eléctrica. En cada momento se asume operación en
régimen permanente.
• Una componente del consumo depende del nivel de tensión existente en el nodo de
carga (impedancia constante) y otra componente del consumo es fundamentalmente
independiente del nivel de tensión existente (potencia constante). La parte del consumo
que depende de la tensión es la carga del tipo calefacción y la que no depende de la tensión
es la carga motorizada.
55
Capítulo 2. Representación de los SEP
• La carga típica es simétrica, es decir, se asume que todos los equipos son trifásicos y que
el consumo en las fases es equilibrado (balanceado). Lo anterior implica que las cargas
monofásicas o que se encuentran conectadas a solo dos fases o a dos fases y neutro,
se encuentran distribuidas de tal forma que el sistema resultante es balanceado. Debe
tenerse en cuenta que, desde el punto de vista estadístico, un número elevado de cargas
asimétricas conectadas arbitrariamente al sistema tienden a producir una disribución
equilibrada de corrientes en las tres fases.
• Presenta una curva de demanda a nivel local y global con una clara diferenciación entre
el consumo mínimo, medio y máximo. En Colombia esta curva de demanda presenta dos
picos como se muestra en la figura 2.19.
4 8 12 16 20 24
T
En estudios de flujo de carga para sistemas de potencia es común que las cargas sean
representadas como potencia constante y en el estudio de redes de distribución es común que
Diferentes
tipos de las cargas sean representadas como potencia, corriente y/o como impedancia constante, o la
estudios combinación de estas.
En estudios de cortocircuito normalmente las cargas no son consideradas. En sistemas
con altos niveles de corto circuito, la corriente que fluye por las líneas debido al corto es muy
superior a la corriente producida por la carga normal. En sistemas con bajo nivel de corto
circuito, la corriente que fluye por las líneas debido a la carga resulta ser significativa respecto
a la de cortocircuito, y por lo tanto debe ser considerada.
En otra clase de estudios, como en análisis armónico, el modelamiento de las cargas es
mucho más complejo y resulta definitivo. En general, según el tipo de análisis que se realiza se
adoptan modelos más simples o más detallados para describir el comportamiento de las cargas.
Los modelos de carga comunmente utilizados para análisis de sistemas eléctricos en
estado estacionario son:
Modelo de inyección de potencia constante:
Este modelo representa generalmente grandes consumos vistos desde las subestaciones.
Potencia
constante Los valores de potencia activa y reactiva que representan la carga son el resultado de la suma de
los consumos de usuarios residenciales, comerciales e industriales que comparten su conexión
con una barra de la subestación. Los valores de P y Q se obtienen a partir de mediciones en
la subestación y datos estadísticos y se representan por curvas de demanda. En este modelo, P
56
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
P = Kp · V pv ; Q = Kq · V qv (2.21)
Si además se supone que las cargas no dependen de la tensión, pv = qv = 0, se obtiene el
modelo de inyección de potencia constante.
P = Kp ; Q = Kq (2.22)
57
Capítulo 2. Representación de los SEP
Si se hace pv = qv = 1 y pf = qf = 0 se obtiene:
P = Kp · V Q = Kq · V (2.23)
Que es equivalente a:
P Q
= Kp = I; = Kq = I (2.24)
V V
Que representa el modelo de corriente constante.
Finalmente si pv = qv = 2 y pf = qf = 0, se obtiene el modelo de admitancia:
P = Kp · V 2 ; Q = Kq · V 2 (2.25)
De aquí se deduce que:
P Q
= Kp = Y ; = Kq = Y (2.26)
V2 V2
A continuación se presentan algunos valores típicos de parámetros de carga.
Carga pv qv pf qf
Motor de inducción a media carga 0,2 1,6 1,5 -0,3
Motor de inducción a plena carga 0,1 0,6 2,8 1,8
Lámpara de filamento 1,7 0 0 0
Lámpara Fluorescente 1,2 3,0 -1,0 2,8
Horno de reducción 1,9 2,1 -0,5 0
Planta de aluminio 1,8 2,2 -0,3 0,6
Calefactor 2,0 0 0 0
Estas características se pueden combinar para representar una carga acumulada conectada
en un nodo particular.
Modelos de elementos en paralelo
En estudios de flujo de carga, los efectos capacitivos de las líneas de transmisión se
Modelo para representan como elementos paralelo. También los modelos de transformadores con taps, los
Elemento en
Paralelo compensadores capacitivos y los compensadores inductivos se modelan usando un equivalente
con elementos paralelo. En estos casos el componente paralelo se modela como inyección
de potencia o por su valor de admitancia. En estudios de corto circuito, generalmente son
despreciados o en algunos casos representados por su valor de admitancia.
58
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
k m
ikm Zkm = Rkm + J Xkm imk
Vk ∠θ k Vm ∠θ m
sh sh
Jbkm Jbkm
59
Capítulo 2. Representación de los SEP
∗
Skm = Pkm − jQkm = Ek∗ ikm (2.29)
∗ −jθk
Skm = Vk e ykm (Vk ejθk − Vm ejθm ) + jbsh
km Vk e
jθk
(2.30)
∗
2
Skm = jbsh 2
km Vk + ykm Vk − Vk Vm e
−j(θk −θm )
(2.31)
Se asume que: θkm = θk - θm
∗
Skm = jbsh 2 2
km Vk + (gkm + jbkm )(Vk − Vk Vm Cosθkm + jVk Vm Senθkm ) (2.32)
∗
Skm = Vk2 gkm − Vk Vm gkm Cosθkm − Vk Vm bkm Senθkm
− j −(bkm + bsh 2
km )Vk + Vk Vm bkm Cosθkm − Vk Vm gkm Senθkm (2.33)
El cálculo del flujo de potencia activa (Pmk ) y reactiva (Qmk ) de la barra m hacia la barra
k es el siguiente:
∗ ∗
Smk = Pmk − jQmk = Em imk (2.36)
∗
Smk = Vm e−jθm ykm (Vm ejθm − Vk ejθk ) + jbsh
km V m e jθm
(2.37)
∗
2
Smk = jbsh 2
km Vm + ykm Vm − Vk Vm e
j(θk −θm )
(2.38)
∗
Smk = jbsh 2 2
km Vm + (gkm + jbkm )(Vm − Vk Vm Cosθkm − jVk Vm Senθkm ) (2.39)
∗
Smk = Vm2 gkm − Vk Vm gkm Cosθkm + Vk Vm bkm Senθkm
− j −(bkm + bsh 2
km )Vm + Vk Vm bkm Cosθkm + Vk Vm gkm Senθkm (2.40)
60
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Qe = −bsh 2 2 2 2
km (Vk + Vm ) − bkm (Vk + Vm − 2Vk Vm Cosθkm ) (2.44)
Qe = −bsh
km (Vk
2
+ Vm2 ) − bkm | Ek − Em | 2
(2.45)
-bsh 2 2
km (Vk + Vm ): Generación de potencia reactiva en el elemento shunt. (Pérdidas negativas,
ya que bsh
km mayor que 0).
61
Capítulo 2. Representación de los SEP
• Cuando la tensión de un nodo o una potencia que fluye por una línea de transmisión
se expresa en por unidad, determinar si este valor se encuentra dentro de los límites
permisibles es más directo que si se utilizan los valores reales. Los valores en p.u. tienen
una interpretación similar a la de cantidades en porcentaje. Por ejemplo, si la tensión de
un nodo se escala respecto a la tensión nominal de dicho nodo, cuando la tensión tiene un
valor igual a la tensión nominal, se dice que la tensión está al 100% o equivalentemente
en 1 p.u. Dado que los niveles máximos y mínimos permitidos están asociados a valores
porcentuales, por ejemplo, si la tensión se considera normal cuando se encuentra entre el
95% y el 105% del valor nominal, esto equivale a decir que debe encontrarse entre 0,95
y 1,05 p.u. En consecuencia, si la tensión calculada en un nodo resulta ser de 0,97 p.u.
se dice que se encuentra dentro del rango permitido y para esto no se requiere conocer su
valor real en kV.
• Los valores de impedancias de equipos del mismo tipo (por ejemplo máquinas síncronas
de polos salientes) con valores nominales muy diferentes tienen valores de impedancia
en ohmios también muy diferentes, sin embargo, si sus impedancias se especifican en
p.u. estos valores resultan muy similares. Esta característica es útil para determinar si los
valores que se están utilizando en el análisis son apropiados y también sirven para asignar
valores típicos de impedancia a un equipo del cual no se tiene información. Diferentes
características de muchos aparatos eléctricos pueden ser estimados por comparación de
otros expresados en p.u.
62
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Z cc ZL Z cc
63
Capítulo 2. Representación de los SEP
Cantidad Unidad en VA y Ω en VA y V
(V A)1/2 (V A)
I A Ω1/2 V
V V Ω .(V A)1/2
1/2 V
S VA (VA) (VA)
V2
Z Ω Ω (V A)
S = Potencia base
V = Tensión base
Z = Impedancia base
I = Corriente base
Si se especifican V y S , entonces la corriente base y la impedancia base, en sistemas
monofásicos, se calculan usando las siguientes expresiones:
V =Z I
S =V I
V2
Z = Impedancia base en Ω = S
LN
1φ
S 1φ
I = Corriente base en A = V LN
La tensión base V siempre representa una magnitud sin ángulo en kV, sin embargo, la
potencia base puede ser una potencia aparente, por ejemplo, la potencia en MVA del equipo
64
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
de mayor capacidad del sistema, o puede ser una potencia activa en MW, por ejemplo, la
potencia activa límite asociada a la línea de transmisión de mayor capacidad, o puede ser una
potencia reactiva, por ejemplo, la potencia reactiva nominal de un compensador reactivo. En
consecuencia, la potencia base puede tener como unidades MVA, MW ó MVAR.
Potencia base ⇒ S =V ×I
Potencia base ⇒ P
Potencia base ⇒ Q
Observaciones sobre la impedancia base:
a. Los valores base corresponden a cantidades sin ángulo, los valores físicos son valores
complejos: magnitud y ángulo o compuestos por una parte real y una imaginaria. En
consecuencia la impedancia en p.u tiene parte real y parte imaginaria. Por esta razón, al
dividir una cantidad compleja por su valor base, se altera su magnitud pero no su ángulo.
Por ejemplo, si se define una impedancia Z como:
Z = R + jX(Ω)
Por lo tanto,
R(p.u.) = R(Ω)/Z (Ω)
X(p.u.) = X(Ω)/Z (Ω)
donde
Z∠θ
Zp.u. ∠θ =
p.u.
Z
La impedancia en p.u. tiene el mismo valor de ángulo que la impedancia en Ω. En
concordancia con lo anterior, se puede afirmar que:
Z =X =R
c. Los valores en p.u. se obtienen dividiendo el valor físico por el valor base. Cálculo de
la impedancia p.u.:
Z S
Zp.u. = p.u. = Z p.u.
Z (V )2
65
Capítulo 2. Representación de los SEP
MVA
Zp.u. = Z p.u.
(kV )2
S
Sp.u. = p.u.
S
I
Ip.u. = p.u.
I
V
Vp.u. = p.u.
V
Similarmente la potencia es representada en sus componentes activa y reactiva y expresada
en p.u. como sigue:
P Q
Pp.u. = y Qp.u. =
S S
La potencia es calculada como el producto de la tensión por la corriente y expresada de la
siguiente manera:
I¯ = I¯p.u. ∗ I (A)
V̄ = V̄p.u. ∗ V (V )
Z̄ = Z̄p.u. ∗ Z (Ω)
R = Rp.u. ∗ Z (Ω)
S = Sp.u. ∗ S (VA)
P = Pp.u. ∗ S (W )
66
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Z MVA base1
Zp.u. base1 = =Z 2
(2.46)
Zbase1 (kV base1 )
2
(kV base1 )
Z = Zp.u. base1 (2.47)
(MVA base1 )
Nota: La palabra base1 se refiere a la base del elemento que será incorporado a la base de
datos de la red, base2 se refiere a la base de la red y con la cual tendrá que quedar el elemento
que será incorporado al sistema.
Para referirlo a la nueva base del sistema se tiene:
Z MVA base2
Zp.u. base2 = =Z 2
(2.48)
Zbase2 (kV base2 )
Comparando los valores base de tensión y potencia del sistema con los valores base
utilizados por el fabricante del equipo se pueden presentar varias situaciones:
67
Capítulo 2. Representación de los SEP
• Formas de referir la impedancia del modelo del transformador mostrado en la figura 2.22
del primario al secundario (figura 2.23) o del secundario al primario (figura 2.24).
IP ZP ZS IS
nP / nS
+ + - +
VP VS
- -
Figura 2.22. Representación del transformador con impedancia dividida
IP ZS IS
+ + -
+ +
VP
VP VS
a
- - -
a :1
Figura 2.23. Impedancia del transformador referida al secundario
68
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
IP ZP = ZS × a2 IS
+ + - +
+
aVS VS
VP
-
- 1: a
-
Figura 2.24. Impedancia del transformador referida al primario
La tensión base del secundario del transformador será entonces una variable dependiente,
determinada a partir de la base del primario y de la relación de transformación nominal (ver
figura 2.25).
+ -
+ +
VP VS
VPbase VSbase
- -
1 :1
Figura 2.25. Modelo del transformador en p.u. visto en el secundario
69
Capítulo 2. Representación de los SEP
VPbase S VSbase
IP × Z P × base 2 Is ×
Sbase VPbase Sbase
+ +
VP VS
VPbase VSbase
- -
1:1
Figura 2.26. Modelo del transformador en p.u. visto en el primario
Nótese que se está considerando el caso particular en que el tap del transformador está en
la posición nominal o sea a = a
Se observa que las impedancias referidas al primario o al secundario del transformador
son iguales, cuando son medidas en p.u.
ZP × P ZS × P
ZP p.u. = ZSp.u. ⇒ 2
=
(VP ) (VS )2
ZP ZS
⇒ =
(VP )2 (VS )2
Lo anterior puede ser verificado con la siguiente relación:
2
2 VP
Zp = a ZS = ZS
VS
+ +
VP VS
VPbase VSbase
- -
Figura 2.27. Modelo del transformador en p.u. con el tap en la posición nominal
Ejemplo 2.1.
70
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
X P = X 1 + a2 X 2
X 1 + a2 X 2 X 1 + a2 X 2
XP p.u. = = MVA
ZP (VP )2
X1 a2 X 2 + X 1
XS = X2 + =
a2 a2
a2 X 2 + X 1 a2 X 2 + X 1
XS p.u. = = MVA
a 2 ZS (VP /VS )2 (VS )2
X 1 + a2 X 2
XS p.u. = MVA
(VP )2
Ejemplo 2.2.
Generador Transformador
Capacidad (MVA) 75 40
Voltaje (kV) 20 20/34,5
Reactancia (p.u.) 0,09 0,20
71
Capítulo 2. Representación de los SEP
(20)2
Generador:Z (Ω) = 0, 09 = 0,48 Ω
75
(20)2
Transformador Z (Ω) = 0,20 = 2,0 Ω (En el primario del
40
transformador)
(34,5)2
Z(Ω) = 0, 2 = 5.95 Ω (En el secundario del
40
tranformador)
Nota: Como se observa, las impedancias en p.u. vistas desde el primario y secundario del
transformador son iguales.A pesar de esto, es importante recordar que los valores base con
los cuales se calcularon dichas impedancias son diferentes, ya que se cuenta con valores
base en el primario y secundario por tal razon los valores reales de las impedancias en Ω
vistas desde cada lado son diferentes.
Ejercicio 2.1. Convertir a p.u. la información del sistema de energía eléctrica presentado en la
figura 2.30.
Datos del sistema:
T1 : Rel. de transf. = 13,8/115 kV, Capacidad = 10 MVA Reactancia X1 = 0,1 p.u.
T2 : Rel. de transf. = 115/34,5 kV, Capacidad = 10 MVA Reactancia X2 = 0,08 p.u.
Zc = 300 Ω
Se selecciona V = 13,8 kV en la región A y una potencia base S = 10 MVA.
Ejercicio 2.2. Dado el sistema de energía eléctrica mostrado en la figura 2.31 determinar:
72
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
T1 T2
A B C ZC
• Valores en p.u. de los parámetros del sistema con base en las nuevas base.
73
Capítulo 2. Representación de los SEP
1 T1 2 3 4
T2
L1
G1 C1
V (L−L)
V (L−n) = √
3
S (3ϕ)
S (1ϕ) =
3
Sin embargo en un sistema trifásico se emplean como valores base la tensión línea - línea
y como potencia base la potencia trifásica, tal como se describe a continuación:
S (3ϕ)
I =√
3V (L−L)
MVA (3φ)
I =√ (kA)
3 kV (L−L)
2
[kV (L−L) ]
Z = Ω =R =X
MVA (3ϕ)
74
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
1 MVA (3ϕ)
Y = = 2
℧=G =B
Z [kV (L−L) ]
Conexión ∆ − Y
x2
x1 x1 x1
3
x2 x2
x1
75
Capítulo 2. Representación de los SEP
+ + +
VP
3 ZC ZC
VP - VC
ZC
- -
b b
√
V p V p / 3
ab = b = √
Vs Vsb / 3
Por lo tanto, es indiferente emplear tensiones de línea o tensiones de fase.
La impedancia de la carga reflejada al primario es dada por: ZcP = (ab )2 Zc
+ + +
VP
3 ZC ZC
VC
VP - VC
3
ZC
- -
Si se supone una relación de espiras tal que las tensiones de línea en el primario y en el
secundario son iguales, este transformador es equivalente al de conexión Y − Y .
Sin embargo, para este caso, la relación del número de espiras entre primario y secundario
es:
√
VbP / 3
a=
VbS
76
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
VbP
aeq =
VbS
Lo anterior corresponde a la relación entre las tensiones de línea (fase-fase) nominales del
primario y el secundario.
De esta forma es posible tratar los transformadores Y − Y e Y − ∆ utilizando el mismo
modelo, es decir, para el caso de un transformador con conexión Y − ∆ en lugar de pensar en
la relación de espiras que realmente existe, se trabaja con la relación equivalente conforme fué
descrito antes.
Nota: Al trabajar en valores por unidad, las cargas trifásicas se representan en Y
equivalente.
• V = VL−N y S =S 3φ /3
• V = VL−L y S =S 3φ
Esta última es la opción más usual. Lo anterior basado en que un sistema 3φ equilibrado
se puede resolver como uno monofásico, puesto que no circula corriente por el neutro.
1 4
G1
M1
2 3
T1 T2
L1
G2
C1
77
Capítulo 2. Representación de los SEP
Carga C1 Resistencia: RC = 2, 0 Ω
b
Línea de Transmisión L1 : Potencia nominal: SL1 = 100 MVA
b
Tensión nominal: VL1 = 115 kV
Reactancia: XL1 = 10% = 0, 1 p.u.
78
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
G1 T1 T2
L1
∆ Υ Υ ∆
G2 T3 T4
L 2 L3 C1
∆ Υ Υ ∆
Υ ∆ Υ ∆
T5 T6
C2
Figura 2.36. Sistema en malla abierta
En el caso de la figura 2.38, a pesar de ser enmallada también puede ser tratada de la
misma forma que el ejercicio anterior, ya que no presenta transformadores en la malla.
La configuración que se observa en la figura 2.37 presenta una dificultad adicional, aqui la
malla es cerrada a través de transformadores. Estos elementos dependiendo de la conexión y/o
localización de los taps pueden introducir un desfasamiento angular y/o elevación (disminución)
del nivel de tensión nodal. Por lo anterior, el cerramiento de mallas a través de transformadores
podrá introducir complejidad adicional en el cálculo de los valores p.u.
Es importante aclarar que las situaciones que involucran los transformadores ∆ − Y
presentados en este caso son de alguna forma artificiales y son presentados con el fin de ilustrar
las dificultades para definir el sistema p.u. en sistemas enmallados.
También se aclara que esta situación es poco probable en la práctica.
ii) Conexión Y − ∆:
79
Capítulo 2. Representación de los SEP
G1 T1 T2
L1
∆ Υ Υ ∆
G2 T3 T4
L 2 L3 C1
∆ Υ Υ ∆
Υ ∆ Υ ∆
T5 T6
C2
Figura 2.37. Sistema en malla cerrada con transformadores en la malla
80
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
G1 T1 T2
L1
∆ Υ Υ ∆ C1a
G2 T3 T4
L 2 L3
C1b
∆ Υ Υ ∆
Υ ∆ Υ ∆
T5 T6
C2a C2b
Figura 2.38. Sistema en malla cerrada sin transformadores en la malla
81
Capítulo 2. Representación de los SEP
T1 L1 T2
T3 T4 C1a C1b
G1 G2 L2 L3
T5 T6
C2a C2b
Figura 2.39. Circuito de malla cerrada sin correción de magnitud de voltaje y ángulo a través
de un transformador desfasador
82
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
T1 L1 T2
T3
A.T .
T4 C1a C1b
G1 G2 L2 L3
T5 T6
C2a C2b
Figura 2.40. Circuito de malla cerrada con correción de magnitud de voltaje y ángulo a través
de un transformador desfasador
Los cinco parámetros que definen el modelo del transformador de tres arrollamientos son:
ZP , ZS , Zt , a s y a t
Las pruebas que son necesarias para la obtención del modelo son:
Zp + Zs = Zps
Zp + Zt = Zpt
(abs )2 Zs + (abs )2 Zt = Zst
83
Capítulo 2. Representación de los SEP
Ip Is
+ +
Vp Vs
- -
It
- Vt +
1: a s +
Zs
Vs
-
+ Zp
Vp 1: a t +
- Zt
Vt
1
Zp = (Zps + Zpt − Zst /(abs )2 )
2
84
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
1
Zs = (Zps − Zpt + Zst /(abs )2 )
2
1
Zt = (−Zps + Zpt + Zst /(abs )2 )
2
En las expresiones anteriores los valores de impedancia de Zps y Zpt están referidos a Ω
del primario y Zst a Ω del secundario.
Cuando Zps , Zpt y Zst son dadas en p.u., las expresiones son dadas como:
1
Zp = (Zps + Zpt − Zst )p.u.
2
1
Zs = (Zps − Zpt + Zst )p.u.
2
1
Zt = (−Zps + Zpt + Zst )p.u.
2
Ejemplo 2.3.
10
Zst = 6% = 8%
7.5
En por unidad, respecto a la base del sistema se tiene:
85
Capítulo 2. Representación de los SEP
1
Zp = J (0,07 + 0,09 - 0,08) = J0,04 p.u.
2
1
Zs = J (0,07 - 0,09 + 0,08) = J0,03 p.u.
2
1
Zt = J (-0,07 + 0,09 + 0,08) = J0,05 p.u.
2
El diagrama en p.u. del transformador tridevanado es presentado en la figura 2.43:
Ficticio
S
1 2 3 4
Carga
5
Banco de
condensadores
86
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
F
4
1 2 3
Ejercicio 2.3. Determinar el diagrama de impedancias en por unidad del sistema de energía
eléctrica mostrado en la figura 2.46.
L1
1 2
T2 T3
L2
3 4
7
T1 T4
Cond.
5 6
C1
G1
Figura 2.46. Diagrama de una red eléctrica de 7 nodos
87
Capítulo 2. Representación de los SEP
a. Unidades monofásicas:
b. Transformador trifásico:
T4 (3ϕ) Transformador de tres devanados.
Relación de transformación: 220/110/44 kV.
Conexión: Y - Y - ∆.
Impedancias en por unidad:
Zps = 0, 18 p.u. (220 kV, 180 MVA)
Zpt = 0, 28 p.u. (220 kV, 180 MVA)
Zst = 0, 16 p.u. (110 kV, 100 MVA)
Generador:
′′
Tensión Factor de Capacidad React.X (%)
(kV) pot. (MW) Subtrans.
G1 13, 8 0, 85 125 25
Carga:
88
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
I1 It I2
y
+ +
+
V1 Vt V2
− −
89
Capítulo 2. Representación de los SEP
k m
ikm ty imk
Vk ∠θ k Vm ∠θ m
(t 2
−t )y (1 − t ) y
+
+ +
V1 Vt V2
- -
-
90
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
i2
a
= −(V1 − aV2 )Y
i2 = −a(V1 − aV2 )Y
k m
ikm y/ t imk
Vk ∠θ k Vm ∠θ m
(1−1/ t) y (1/ t 2
− 1/ t )y
91
Capítulo 2. Representación de los SEP
100
Zt1 p.u. = 0, 118 = J 0, 472 p.u.
25
• Transformador 2 cambio de base:
Se supone que el voltaje base con el cual se obtuvieron los valores en p.u. corresponden
al primario. Así solo se requiere cambio de base de potencia.
100
Zt2 p.u. = 0, 064 = J 0, 32 p.u.
20
k m
ikm y p.u.(1,09) imk
Vk ∠θ k Vm ∠θ m
Figura 2.51. Modelo del transformador 2 que simula la relación de transformación 115/37, 6
kV
92
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
k m
ikm - J 3,406 imk
Vk ∠θ k Vm ∠θ m
- J 0,30656 J 0,28125
T1 J 0 , 472
I1 I2
1 1,09 T 2 J 0,320
+ +
Vp Vs
-
-
93
Capítulo 2. Representación de los SEP
Transformador serie
A
Pkl = PF l y Qkl = QF l
94
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
k l
f
1 : e jθ kl
Pkl Pfl
Qkl Transformador ideal Qfl
VF = Vk a ; PF l = PKl
θF = θk + φkl ; QF l = Qkl
Sustituyendo estos valores en las expresiones PF l y QF l , se obtiene:
Pkl = (Vk a)2 gkl − (Vk a)Vl gkl Cos(Θkl + φkl ) − (Vk a)Vl bkl Sen(Θkl + φkl )
Qkl = −(Vk a)2 bkl − (Vk a)Vl gkl Sen(Θkl + φkl ) + (Vk a)Vl bke Cos(Θkl + φkl )
95
Capítulo 2. Representación de los SEP
+
+ + +
VPauto VP NP NS VS
− − − VSauto
−
Figura 2.56. Autotransformador elevador
+ +
VS
−
VSauto
+
+
VPauto
VP
−
− −
NP
La relación de transformación para el transformador común es: r1φ = NS
.
La relación del autotranformador es: rauto = NPN+N
P
S
VP NP 1
rauto = = = NS
=
VP + VS NP + NS 1+ N P
1 1 r1φ
1 = r1φ +1 =
1 + r1φ r1φ + 1
r1φ
96
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
is
Vs
Como se observa, los puntos a y b se encuentran al mismo potencial, por lo tanto la fuente
de tensión VS está conectada directamente al devanado secundario. Así, ZSPauto = ViSS (igual
que en el caso del transformador común).
En la figura 2.59 se presenta el circuito equivalente del autotransformador.
Z SP
+ +
VPauto NP NP + NS VSauto
− −
97
Capítulo 2. Representación de los SEP
iSauto = iS
+ +
VS
NS
iPauto −
VSauto = VP + VS
+ iP
+
VPauto = VP
VP
NP
−
− −
La expresión (2.49), dada en función de las variables del transformador común (VPauto =
VP , VSauto = VP + VS , iPauto = iP + iS , iSauto = iS ), puede ser presentada de la siguiente forma:
98
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
ZSPauto ZSP
Zautop.u. = 2 = 2
(VPbase +V )
NP
S V
NS Sbase
+VSbase
base
NP
NP
NS
+1 S1φbase +1 S1φbase
NS
2 !2
NS
NP +NS
ZP S 1 ZP S
Zautop.u. = VP2
= NP VP2
base
NS
+1 base
NP NP
NS
+1 S1φbase NS
+1 S1φbase
ZP S ZSP
Zautop.u. =
NP
2 =
NP
2
NS
+1 VP2 NS
+1 VS2
base base
NP NP
NS
+1 S1φ base NS
+1 S1φ base
ZP S ZSP
Zautop.u. = VP2
= VS2
NP NP
NS
+1 base
S1φbase NS
+1 base
S1φbase
ZP Sp.u. ZSPp.u.
Zautop.u. = = (2.53)
NP NP
NS
+1 NS
+1
Así, en el transformador monofásico se tiene que:
Ejemplo 2.4.
99
Capítulo 2. Representación de los SEP
iSauto = iS
+ +
NS 33 kV
iPauto −
VSauto = VP + VS
+ iP
+
VPauto = VP
NP 67 kV
−
− −
100
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
NS
i′p
VS
+ iP VS′
VP
NP
− − −
P S
tY
(1 − t ) Y
(t 2
− t )Y
A partir del análisis de la figura 2.62, al aplicar la primera y segunda Ley de Kirchhoff y
la Ley de Ohm, se tiene lo siguiente:
it = (Vs′ − Vs )Y (2.54)
is = −(Vs′ − Vs )Y (2.55)
N p ip = N s it (2.56)
i′p = ip + it ⇒ ip = i′p − it (2.57)
101
Capítulo 2. Representación de los SEP
Np (i′p − it ) = Ns it
Np i′p − Np it = Ns it
Np i′p = (Np + Ns )it
i′p Np + Ns
=
it Np
Respecto a la relación de voltajes:
Vp Np
=
Vs′ Np + Ns
Np +Ns
Por definición t = Np
, por lo tanto:
i′p
it = (2.58)
t
′
Vs = t Vp (2.59)
i′p
it = (Vs′− Vs )Y ⇒ = (t Vp − Vs )Y ⇒ i′p = (t2 Vp − tVs )Y
t
′
is = −(Vs − Vs )Y ⇒ is = −(t Vp − Vs )Y ⇒ is = (−tVp + Vs )Y
Estas expresiones pueden escribirse matricialmente, obteniéndose una representación
matemática, la cual obedece a una relación de variables nodales así: ibus = [Ybus ] Vbus .
′ 2
ip t Y −tY Vp
= (2.60)
is −tY Y Vs
Finalmente, el modelo del regulador en representación primitiva es mostrado en la figura
2.63
102
REPRESENTACIÓN MATRICIAL DE
3
LAS REDES ELÉCTRICAS
3.1 Introducción
Las redes de energía eléctrica son definidas por medio de modelos matemáticos, que pueden ser
de diferente tipo, dependiendo de la aplicación y exactitud requerida en los resultados. Desde el
punto de vista de la formulación matemática los modelos pueden ser algebraicos, lineales o no
lineales, o integrodiferenciables lineales o no lineales. Los modelos algebraicos lineales pueden
ser representados por medio de arreglos matriciales, expresados en función de variables nodales
como la tensión y las corrientes inyectadas. Dependiendo del tipo de estudio que se realiza, el
modelo puede ser más o menos detallado, es decir, puede contener expresiones más complejas o
más simples. También es posible considerar modelos reducidos con menos variables o modelos
completos que consideren todas las variables del problema.
Una forma de representar matemáticamente las redes eléctricas usando modelos lineales
es aquella que relaciona las variables nodales de tensión y de corriente en las cuales la corriente
se expresa en función de la tensión o al contrario. En el primer caso, ibus = Ybus × Vbus
expresa la corriente nodal o inyectada en función de las tensiones nodales, en el segundo caso,
Vbus = Zbus ×ibus expresa las tensiones nodales en función de las corrientes nodales o corrientes
inyectadas.
La representación de la red se lleva a cabo a través de las matrices admitancia e
impedancia nodal.
La matriz [Ybus ] representa la conectividad física entre los nodos del sistema, en cambio
la matriz [Zbus ] representa la conectividad eléctrica entre los nodos del sistema. Para que el
elemento Yi−j de la matriz [Ybus ] sea diferente de cero, debe existir conexión física directa entre
los nodos i − j a través de un elemento de circuito. Si no existe conexión física directa entre
estos nodos el elemento Yi−j es igual a cero, sin embargo, el elemento Zi−j de la matriz [Zbus ]
103
Capítulo 3. Representación matricial
puede ser diferente de cero si existe una ruta, así no sea directa, que interconecte estos dos
nodos.
Algoritmos usados en la formación de las matrices YBus y ZBus :
– Zbus = [Ybus ]−1 (inversa que usa técnicas de matrices dispersas y factorización
triangular)
– Algoritmo basado en inyecciones de corriente, (ver sección 3.10).
104
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
IA = YAA VA + YAB VB (3.2)
0 = YBA VA + YBB VB (3.3)
Ejemplo 3.1.
1 2
3 4
105
Capítulo 3. Representación matricial
equivalente −1
YAA = YAA − YAB YBB YBA
Y14 [ Y144 ] [ Y41 Y42 Y43 ]
−1
YAB YBB YBA = Y24
Y34
Y14 Y41 Y14 Y42 Y14 Y43
Y44 Y44 Y44
−1 Y24 Y41 Y24 Y42 Y24 Y43
YAB YBB YBA = Y44 Y44 Y44
Y34 Y41 Y34 Y42 Y34 Y43
Y44 Y44 Y44
Y14 Y41
Y11 − Y44
Y12 − Y14Y44Y42 Y13 − Y14 Y43
Y44
equivalente Y24 Y41
YAA = Y21 − Y44
Y22 − Y24Y44Y42 Y23 − Y24 Y43
Y44
Y34 Y41
Y31 − Y44
Y32 − Y34Y44Y42 Y33 − Y34 Y43
Y44
Ejercicio 3.1. Dado el sistema mostrado en la figura 3.2,se pide construir la matriz Ybus y
eliminar de esta el nodo 5. Usar para esto el método de reducción de Kron.
Línea R X Y /2
1−2 0, 0316 0, 2114 0, 105
1−3 0, 0152 0, 0944 0, 072
2−3 0, 0121 0, 0700 0, 082
2−4 0, 0140 0, 1250 0, 096
2−5 0, 0230 0, 055 0, 043
3−4 0, 0246 0, 150 0, 175
4−5 0, 0102 0, 170 0, 202
106
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
1 2 5
3 4
El ordenamiento de las ecuaciones nodales tienen por objetivo presentar el orden en que las
ecuaciones deberán ser solucionadas a fin de minimizar el número de nuevos elementos que
puedan resultar al momento de efectuar los procedimientos de solución de dichas ecuaciones.
Las ecuaciones mencionadas anteriormente son del tipo AX = b, en estas la variable X
corresponde a la incógnita y resuelta a través de remplazos sucesivos de las variables de una
ecuación en las otras y que es descrito como sigue: X = A−1 b.
En el proceso de inversión de la matriz A algunas posiciones con valor cero pasarán
a tomar un valor diferente de cero y son denominados elementos de relleno. Con el fin de
minimizar el número de estos elementos, las filas y columnas de la matriz A son ordenadas
usando algún algoritmo de ordenamiento nodal.
Los primeros esquemas de ordenamiento para aplicaciones en ingeniería eléctrica fueron
propuestos por W. Tinney y denominados por él mismo como Tinney I, II y III respectivamente.
En el primero, el orden es asignado según el número de conexiones fisicas al nodo y su orden
es de menor a mayor número. En el segundo esquema de ordenamiento, es asignado a la lista
de ordenamiento un solo nodo a la vez, y cada vez que esto sucede es simulada la inversa de la
matriz y la misma se va reduciendo en tamaño de uno. En el tercer esquema de ordenamiento se
efectua un estudio exhaustivo simulando en cada caso la primera posición para todos los nodos.
Despues de seleccionado un nodo, se repite el proceso con los nodos no ordenados, de esta
forma el tamaño de las matrices se reduce en uno.
Ejemplo 3.2.
107
Capítulo 3. Representación matricial
2 3
4 1 5
1 x x x x x
2
x x 0 x 0
3
x 0 x 0 0
4 x x 0 x 0
5 x 0 0 0 x
Al simular la inversa de esta matriz se observa que en algunas posiciones cambia de valor
cero a uno diferente.
En la siguiente matriz se muestran las posiciones marcadas con (*) las cuales durante la
simulación de la inversa pasan de cero a un valor diferente de este.
1 x x x x x
2 x x ∗ x ∗
3 x ∗ x ∗ ∗
4 x x ∗ x ∗
5 x ∗ ∗ ∗ x
Caso 2: En la figura 3.4 se muestra el sistema después de ser ordenado con esquema
denominado tinney I.
Este sistema ha sido ordenado usando un esquema de ordenamiento básico. Como se
observa se trata de un sistema simple, levemente enmallado, por lo que cualquier tipo de
ordenamiento nodal seguramente encuentra un orden en el cual no se presenten elementos de
relleno.
La estructura física representada en una matriz es como sigue:
1 x x
2 x x
3 x x x
4 x x x
5 x x x x x
Invertida la matriz anterior no se observan cambios en los elementos iguales a cero.
108
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
3 2
4 5 1
1 x x
2
x x
3
x x x
4 x x x
5 x x x x x
Como conclusión se tiene que, después de ordenados los nodos y simulada la inversa de
dicha matriz, no se observa alteración en los elementos que son iguales a cero. Sin embargo,
cabe destacar, que este no es el caso general de las redes eléctricas en los cuales a pesar de que
se efectúe un ordenamiento óptimo aparecen algunos elementos de relleno. Lo que se trata en
este caso es entonces de minimizar el número de éstos elementos.
Un tipo de estructura especial son las redes radiales, en las cuales despues de
efectuado el ordenamiento no se presentan cambios en los elementos iguales a cero.
La conclusión anterior es aplicable a cualquier tipo de red cuya estructura sea radial.
Otros temas que son importantes dentro de las técnicas de representación matricial
utilizadas en el análisis de sistemas eléctricos de potencia comprenden: factorización
triangular, la descomposición de Cholesky, el algoritmo de construcción de la matriz Ybus ,
cálculo de la matriz Zbus mediante el proceso de factorización LDU, correccióon de la Ybus
por efectos mutuos, entre otros. En Anexo 1 del tomo II de esta obra se analizan estos y otros
aspectos, de forma que el lector interesado pueda profundizar en lo relacionado con técnicas de
representación matricial.
109
Capítulo 3. Representación matricial
110
OPERACIÓN DE LOS SISTEMAS DE
4
POTENCIA
Normal
Control preventivo
Restaurar cargas
E, I E, I
Restauración Alerta
Control de
Resincronización emergencia
E, I E, I
Colapso Emergencia
En más del 99% del tiempo el sistema opera en el estado normal. En este estado tanto la
frecuencia como la tensión presentan los valores preestablecidos para funcionamiento normal.
Estos valores de frecuencia y tensión resultan de mantener un balance entre la potencia activa y
111
Capítulo 4. Operación de los sistemas de potencia
reactiva demandadas por las cargas y la potencia activa y reactiva abastecida por las fuentes de
generación.
La igualdad entre generación y demanda es un prerrequisito fundamental para que el
sistema opere normalmente, indicado por el símbolo (E). El segundo símbolo (I), indica que
ciertas inigualdades pueden presentarse en el estado normal.
El estado normal también puede ser caracterizado por un cierto nivel de seguridad que
requiere un cierto margen de generación en la forma de reserva rodante.
Si el control preventivo falla o si ocurre un fallo severo, el sistema en ocasiones podrá
entrar en un estado de emergencia. Este estado de emergencia podrá ser alcanzado directamente
Estado
de
del estado normal o a través del estado de alerta.
emergencia Si las acciones de control de emergencia fallan, el sistema entra en colapso,
desconectando partes del sistema y conformando sistemas aislados. Algunos de estos sistemas
aislados pueden contener suficiente generación para abastecer las cargas (Monticelli and García,
Colapso
1999).
Si muchos de los generadores salen de operación, los que queden operando son
sobrecargados, entonces el sistema podrá entrar en un colapso total o blackout.
Blackout
La cadena de eventos que lleva el sistema de un estado normal al colapso puede ocurrir
en unos pocos segundos o en varios minutos. El proceso de restauración es mucho más lento,
que puede tomar algunas horas y en ocasiones hasta dias.
Algunos de los problemas a estudiar en la operación de los sistemas de potencia en estado
normal son los siguientes:
• Relación MW - Frecuencia
• Despacho económico
• Motores en general
• Equipos de calentamiento
• Equipos electrónicos
• Equipo de iluminación
112
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Desde el punto de vista del sistema de energía eléctrica, las cargas pueden ser separadas
en tres (3) grupos:
• Cargas domiciliarias
• Cargas industriales
4 8 12 16 20 24
T
Figura 4.2. Curva de carga típica
Estas cargas presentan características muy diferentes con relación al tamaño, simetría (1φ
o 3φ), constancia de la carga y el período de funcionamiento.
En un sistema eléctrico bien proyectado, las cargas presentan las siguientes características:
• Las variaciones de carga generalmente son muy lentas cuando son comparadas con las
constantes de tiempo del sistema de energía eléctrica; así este es considerado como
operando en régimen permanente.
• La carga típica consume potencia reactiva; esto es, la carga típica es inductiva por la
presencia de los motores eléctricos.
113
Capítulo 4. Operación de los sistemas de potencia
Ek∗ Ek2 Vk2
S = Ek = = (4.2)
Z∗ Z∗ (R − jωL)
Vk2 (R + jωL) RVk2 ωLVk2
S= . = 2 + j (4.3)
(R − jωL) (R + jωL) R + (ωL)2 R2 + (ωL)2
dado que S = P + jQ, entonces:
RVk2 ωLVk2
Pk = ; Qk = (4.4)
R2 + (ω.L)2 R2 + (ω.L)2
Como ω = 2π.f, entonces:
RVk2 ωLVK2
Pk = ; Qk = (4.5)
R2 + (2π.f.L)2 R2 + (2π.f.L)2
Observaciones:
• Pk y Qk crecen o decrecen con el cuadrado de la tensión aplicada Vk .
• Cuando la frecuencia f aumenta, Pk disminuye, y Qk lo hace en menor proporción.
En general para todas las cargas existen relaciones del tipo:
114
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Observaciones:
Para el caso de un conjunto de cargas, simultáneamente no es posible obtener relaciones
analíticas. En este caso son obtenidas empíricamente valores típicos para las derivadas
parciales, de acuerdo con los tipos de carga (cargas compuestas).
f = 60+0, 02Hz.
Si el equilibrio generación total - demanda total se pierde por causa de una reducción de
la demanda del sistema, aparecen los siguientes efectos instantáneos:
115
Capítulo 4. Operación de los sistemas de potencia
Turbina
ω Generador Carga
• Los generadores siguen suministrando la energía que entregaban antes del cambio en
la carga pero la carga consume menos energía. La energía generada que ya no se
está consumiendo es almacenada en forma mecánica en las volantes de inercia de los
generadores y asume la forma de energía cinética.
• Los generadores siguen suministrando la energía que entregaban antes del cambio en
la carga pero la carga consume más energía. La energía requerida por la carga que no
116
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Observación importante:
Dentro de los aspectos relacionados con la calidad en el suministro de la energía eléctrica
se encuentra el control de la frecuencia en su valor nominal. En consecuencia, un aumento o una
disminución de la frecuencia más allá del umbral de variación permitido por la regulación debe
producir una respuesta de los actuadores de sistema con el propósito de reducir o aumentar la
energía primaria que se entrega a los generadores y así ajustar la energía generada para que sea
igual a la demandada por el sistema. En este proceso, el aspecto fundamental está asociado con
el control de las válvulas que permiten más o menos entrada de agua a las turbinas hidráulicas y
con el control de las válvulas que permiten mayor o menor paso de combustible en las centrales
térmicas. Todo esto es adecuadamente coordinado por un sistema de control.
117
Capítulo 4. Operación de los sistemas de potencia
k m
Ek = Vk ∠0 Em = Vm∠δ
P + jQ
Figura 4.5. Sistema simplificado de dos barras
Simplificaciones:
• Vk es mantenido constante
Ek i∗km = P + jQ (4.11)
De la ecuación anterior se obtiene:
P + jQ P + jQ P − jQ
i∗km = = o
⇒ ikm = (4.12)
Ek Vk ∠0 Vk
Reemplazando las ecuaciones anteriores:
o P − jQ
Em = Vk ∠0 − (jXkm ) (4.13)
Vk
Xkm Xkm
E m = Vk − Q −j P (4.14)
Vk Vk
De la gráfica se observa que:
118
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
X km
Q
Em Vk
E m′ X km
P
Vk
• Con carga pico (máxima), aumentan las necesidades de potencia activa y reactiva;
entonces, la tendencia de Vm es a disminuir y por lo tanto, se debe colocar capacitores
sincronos en m para generar la potencia reactiva necesaria.
• Con carga baja (mínima), el efecto capacitivo de la línea de transmisión genera potencia
reactiva mayor que la necesaria en la carga, y en este caso el flujo de potencia reactiva
puede ser invertido, lo que lleva a un aumento de tensión; entonces puede ser necesario
instalar consumidores de potencia reactiva en la barra m, esto es, se deben colocar
reactores en paralelo.
Dado:
E2 = V2 ∠0o = 1∠0o y E1 = V1 ∠θ1 − θ2 = 1∠δ
119
Capítulo 4. Operación de los sistemas de potencia
G1 G2
1 2
P12 = 10 p.u. y 1
b12 = − X112 = − 0.03 = − 100
3
120
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Ejercicio:
Ejercicio 4.1. Otra forma de resolver el ejemplo anterior (4.1) es utilizando las ecuaciones de
flujo de potencia activa y reactiva a través de una línea, descritas en la sección 1.15 (pág. 29).
Resolver el ejemplo utilizando estas ecuaciones.
121
Capítulo 4. Operación de los sistemas de potencia
122
SOLUCIÓN AL PROBLEMA DE FLUJO
5
DE POTENCIA
123
Capítulo 5. Solución al problema de FP
Los elementos de la matriz [Ybus ] pueden ser descritos en forma polar o rectangular, así:
El modelo matemático del sistema con el cual se resuelve el problema de flujo de carga se
expresa en función de las tensiones nodales, las cuales pueden ser representadas en componentes
polares o rectangulares, siendo la más usada en los escritos técnicos la representación polar.
n
X
I~i = Y~i1 V~1 + Y~i2 V~2 + · · · + Y~in V~n = Y~in V~n (5.4)
n=1
124
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
i i
PGi Pi1 1 1
QGi Qi1
Pi 2 2 Qi 2 2
Pij M j
M
Q ij j
M M
PDi QDi
1 1 1 1
Pi Programada Pi Calculada Qi Programada Qi Calculada
125
Capítulo 5. Solución al problema de FP
n
X
Pk + JQk = V~k Y~km
∗ ~∗
Vm (5.14)
m=1
n
X
Pk + JQk = Vk eJθk Ykm e−Jϕkm Vm e−Jθm (5.15)
m=1
Xn
Pk + JQk = |Vk ||Ykm ||Vm |∠θk − θm − ϕkm (5.16)
m=1
n
X
Pk = |Vk ||Ykm ||Vm |Cos(θk − θm − ϕkm ) (5.17)
m=1
Xn
Qk = |Vk ||Ykm ||Vm |Sen(θk − θm − ϕkm ) (5.18)
m=1
n
X
Pk + JQk = V~k ∗ ~∗
Y~km Vm (5.19)
m=1
n
X
Pk + JQk = Vk eJθk Ykm e−Jϕkm Vm e−Jθm (5.20)
m=1
Xn
Pk + JQk = Vk Vm (Cosθkm + JSenθkm )(Gkm − JBkm ) (5.21)
m=1
n
X
Pk = Vk Vm (Gkm Cosθkm + Bkm Senθkm ) (5.22)
m=1
Xn
Qk = V k Vm (Gkm Senθkm − Bkm Cosθkm ) (5.23)
m=1
126
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Por cada nodo son expresadas dos ecuaciones como se presenta a seguir:
n
X
Pk = Vk Vm (Gkm Cosθkm + Bkm Senθkm ) (5.25)
m=1
Xn
Qk = V k Vm (Gkm Senθkm − Bkm Cosθkm ) (5.26)
m=1
127
Capítulo 5. Solución al problema de FP
Fortalezas: de fácil implementación, requiere poca memoria, las respuestas obtenidas son
proximas a las reales, por tratarce de un modelo no lineal. Debilidades: requiere de gran
tiempo de cómputo y diverge fácilmente por malos condicionamientos.
Entre los grandes aportantes al flujo de carga resuelto usando el método de Newton
Rahson están Tinney y Walker quienes presentaron la solución a dicho problema, con
la aplicación de las matrices dispersas, esquemas de ordenamiento nodal y algoritmos
eficientes de programación.
- También existen otros métodos alternativos que han ayudado en la discusión del tema,
como es el caso del flujo de carga fasorial Iwamoto and Tamura (1981).
128
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
129
Capítulo 5. Solución al problema de FP
F1 (X1 X2 · · · , Xn ) = Y 1
.. .. (5.31)
. .
Fn (X1 , X2 · · · , Xn ) = Yn
130
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
F1 (X1 , · · · , Xn ) = Y1
.. .. .. (5.32)
. . .
Fn (X1 , · · · , Xn ) = Yn
F̄ (X̄) = Ȳ (5.33)
• Soluciones múltiples
• Para su solución pueden ser usados, métodos exactos o métodos iterativos. Para este caso
se utilizan los segundos. (En estos, la solución final depende del estado inicial, así el
sistema generalmente converge a una solución cercana al estado inicial. Esta condición
es favorablemente aprovechada en la solución de problemas en redes eléctricas, en las
cuales la magnitud de la tensión está próxima a 1 p.u. y el ángulo en radianes cercano a
cero).
• Newton Raphson
• No lineales
1. k = 0
131
Capítulo 5. Solución al problema de FP
F( X )
El proceso diverge
Y=X
El proceso converge
X4 X3 X2 X1 X0
X
Figura 5.2. Proceso de convergencia del método de Gauss
132
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Con gran poder de convergencia (el número de iteraciones es independiente del número
de ecuaciones a resolver).
Sistema unidimensional
Sea F(x) = 0, una Función no lineal (una ecuación con una incógnita). La solución de este
problema consiste en determinar el valor de x para el cual la función cumple la condición de
igualdad, esto es F(x) = 0
El problema es resuelto linealizando la función alrededor de un punto e iterando al rededor
del mismo. La linealización de la función anterior se explica a través de la expansión en serie
de Taylor así:
2
k k k ′ k ∆X k k ′′ k
g(X + ∆X ) = g(X ) + g (X )(∆X ) + g (X ) + ...... (5.36)
2!
Se desprecian los términos con potencias superiores a 2, quedando entonces la siguiente
expresión:
′
g(X k + ∆X k ) = g(X k ) + g (X k )(∆X k ) = 0 (5.37)
1. v = 0
133
Capítulo 5. Solución al problema de FP
X V +1 = X V + ∆X V
′
Siendo la corrección de la variable de estado: ∆X V = −g(X V )/g (X V )
6. v = v + 1 , regresar a 2.
Una variante al método es obtenida considerando constante la derivada durante todo el
proceso.
F( X )
Pendiente
Error
X0 X
∆X 4 ∆X 3
∆X 2 ∆X 1
Sistema n-dimensional
El modelo matemático es expresado de la siguiente forma:
ḡ(X̄) = 0 (5.39)
134
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Ejemplo 5.1.
∂g2 ∂g2
g2 (X1 , X2 ) ≈ g2 (X1V , X2V ) + ∆X1V + ∆X2V
∂X1 ∂X2
Algoritmo para solucionar el sistema de ecuaciones g (X) = 0 por Newton.
135
Capítulo 5. Solución al problema de FP
Ejercicio:
Ejercicio 5.1. Resolver el siguiente sistema de ecuaciones no Lineales de segundo orden
empleando los métodos de Newton, Gauss y Gauss-Seidel y analizar los resultados.
x21 − x22 + x1 + 1 = 0
2x1 x2 + x2 = 0
PG , QG
PD , QD V ariables nodales
|V |, θ
De las seis variables nodales mencionadas anteriormente, dos son conocidas, estas son la
demanda de potencia activa y reactiva: PD y QD
Así las ecuaciones nodales se plantean en función de las siguientes cuatro variables
nodales
−Potencia inyectada (SG − SD )
· activa
PG , QG · reactiva
|V |, θ
−Tension
· magnitud
· ángulo
- Variables de control: PG y QG
- Variables de estado: |V | y θ
136
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Generalmente como nodo slack es seleccionado el nodo de mayor generación. Este nodo
generador deberá tener la suficiente cantidad de reactivos para el control de la magnitud
de la tensión.
Se supone conocida la potencia activa generada en todos los nodos generadores, excepto
en el nodo denominado slack, que corresponde al nodo generador de mayor capacidad;
además las demandas en los nodos son conocidas. La potencia generada en el sistema
excepto en el nodo slack es:
n
X
Potencia generada conocida PkG = PG (5.46)
i=1
i 6= slack
n
X
Potencia demandada PkD = PD (5.47)
i=1
Al reescribir la ecuación 5.45, el balance global del sistema cumple con la siguiente
ecuación:
PG + PkG − PkD = P pérdidas (5.48)
slack
De la ecuación anterior, el nodo slack debe generar una potencia para garantizar el balance
global del sistema:
PG = −PkG + PkD + P pédidas
slack
137
Capítulo 5. Solución al problema de FP
Por lo tanto la potencia generada en el nodo slack deberá ser la suficiente para cubrir
la diferencia entre (−PkG + PkD ) más el valor de las pérdidas. Las pérdidas son una
variable no conocida, por lo tanto, la variable PG tendrá que ser libre para que se
slack
pueda cumplir la ecuación anterior.
Al ser la magnitud de la tensión una variable especificada, el QG sera una variable
slack
libre (no especificada). Los nodos generadores incluyendo el nodo slack se encargan
de establecer el balance global del sistema, tanto en potencia activa como en potencia
reactiva, cubriendo los faltantes (sobrantes) que presente el sistema.
De esta manera, en el nodo slack las variable quedan clasificadas así:
|V |, θ ⇒ Variables especificadas (el ángulo θ se asume a 0o y será referencia angular para
los demás nodos del sistema).
P, Q ⇒ Variables libres.
A este tipo de barra pertenecen aquellos nodos que tengan elementos con capacidad para
controlar la magnitud de la tensión.
Barras candidatas a ser denominadas de tensión o voltaje controlado:
– Generadores síncronos
– Motores síncronos
– Compensadores sincrónicos
– Compensadores estáticos activos (controlados por tiristores)
138
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Barras de carga (P-Q) donde la magnitud de la tensión se requiere que este dentro de unos
límites prefijados: Vmin ≤ V ≤ Vmax . Este tipo de barra permite determinar la cantidad
de potencia reactiva requerida en la red a fin de mantener la tensión dentro del límite
establecido.
En el problema de flujo de potencia este tipo de nodo es operado de la siguiente manera:
Cuando uno de los límites de tensión de barra es sobrepasado, el nodo cambia de tipo
(P-Q) a (P - V), Q queda entonces libre y se fija V. Resuelto el flujo se calcula la diferencia
(Qespecif icado − Qcalculado ), que determina la cantidad de reactivos (L-C) que deberán ser
inyectados.
n
X
Pk − JQk = V~k∗ Y~km V~m (5.50)
m=1
n
X
Pk − JQk = V~k∗ Y~kk V~k + V~k∗ Y~km V~m (5.51)
m=1
m 6= k
n
Pk − JQk X
= Y~kk V~k + Y~km V~m (5.52)
V~ ∗
k
m=1
m 6= k
De la ecuación anterior, se despeja la tensión del nodo k:
n
1 P − JQ X
k k
V~k = ~ ~
− Y V
km m (5.53)
Y~kk ~
Vk ∗
m=1
m 6= k
En los nodos tipo tensión controlada se calcula la potencia reactiva generada:
n
P
Qk = V k Vm (Gkm Senθkm − Bkm Cosθkm )
m=1 (5.54)
Siendo (Qk = QGk − QDk )
En estos nodos se deberá satisfacer que:
QGk min ≤ QGk ≤ QGk max (5.55)
139
Capítulo 5. Solución al problema de FP
Ejercicio:
Ejercicio 5.2. Dada la red de energía eléctrica de 3 nodos mostrada en la figura 5.4, determinar
los valores de tensiones nodales, flujos de potencia por las líneas y pérdidas. Emplear para su
solución el método de Gauss Seidel.
Línea Rp.u. Xp.u. Y /2p.u.
1 − 2 0, 01 0, 06 0, 005
Pbase = 100 MVA
1 − 3 0, 03 0, 20 0, 01
2 − 3 0, 04 0, 25 0, 015
Datos nodales
N odo PG QG PD QD Tensión
(M W ) (MVAR) (M W ) (MVAR) (p.u.)
1 150 80 60 30 1, 02
2 30 15 20 10 1, 01
3 − − 45 15 −
140
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Información nodal
V=V+1
K=1 (Contador de nodos)
Si
¿Es K el nodo
K=K+1
Slack?
No
Si
¿Es K el nodo
Calcular
Generador?
n
No Qk = Vk ∑Vm (Gkm Senθ km − BkmCosθ km )
m =1
Calcular
K=K+1 Qk min ≤ Qk ≤ Qkmáx
1 Pk − jQ k n Límite violado
Vk u = − ∑ Ykm V m
Ykk V ∗
k m =1, ( m ≠ k )
Qk = Qkmáx ó Qk = Qk min
¿Se calcularon
tensiones en todos los
nodos?
l =0
l = l +1
Vlu
Vlu = Vl0 Especificado
Vlu
l ≥ nodos de
generación
141
Capítulo 5. Solución al problema de FP
i
SGi S i1 1
Si 2 2
M
S ij j
M
SDi
1 1
Sistema externo Sistema interno
En estado estacionario la potencia inyectada por el sistema externo deberá ser igual a la
inyectada por el sistema interno, se debe cumplir entonces que:
n
P
Sni = Vi eJθi Yim e−Jϕim Vm e−Jθm
m=1 (5.58)
(Fórmula que establece el balance nodal)
El modelo matemático está constituido por 2n ecuaciones nodales, estas son clasificadas
según el tipo de nodo como ecuaciones del subsistema 1 y subsistema 2. Las ecuaciones del
subsistema 1 establecen el modelo matemático del problema de flujo de carga y son escritas
como sigue:
142
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Las ecuaciones del subsistema 2 son el complemento de las 2n ecuaciones nodales y son
usadas en el cálculo de los resultados complementarios.
X ¯k+1 = X¯k + ∆X
¯k
• Las potencias netas inyectadas dependen de su propia tensión, así como tambien de las
tensiones de los otros nodos con los cuales dicho nodo está interconectado y parámetros
de la red.
• Si los nodos i y k no están interconectados no aparece término Yik en la matriz Ybus por
lo tanto el término (i − k) no incide en la sumatoria de la potencia nodal.
143
Capítulo 5. Solución al problema de FP
Ejemplo 5.2.
1 2 3
∆P1 δ P1 δ P1 δ P1 δ P1
∆θ1
δθ1 δθ 2 δ V1 δ V2
∆P2 δ P2 δ P2 δ P2 δ P2 δ P2 δ P2
∆θ 2
δθ1 δθ 2 δθ3 δ V1 δ V2 δ V3
∆P3 δ P3 δ P3 δ P3 δ P3 ∆θ3
δθ 2 δθ3 δ V2 δ V3
=
∆Q δ Q1 δ Q1 δ Q1 δ Q1
∆V
1 δθ δθ 2 δ V1 δ V2 1
1
∆Q2 δ Q2 δ Q2 δ Q2 δ Q2 δ Q2 δ Q2
∆V2
δθ δθ 2 δθ3 δ V1 δ V2 δ V3
1
∆Q3 δ Q3 δ Q3 δ Q3 δ Q3 ∆V3
δθ 2 δθ3 δ V2 δ V3
Como se observa de la figura anterior, las ecuaciones nodales están descritas en función de
(P, Q, V y θ). El modelo matemático del problema de flujo de carga toma la siguiente estructura
matricial.
∆P H N ∆θ NPQ + NPV
∆ Q = J L ∆V
NPQ
Los errores de potencia (∆P y ∆Q) son calculados como la diferencia entre la potencia
inyectada por el sistema externo (especificado) e interno (calculado) de la siguiente manera:
∆P = P especificado − P calculado
144
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
∆Q = Qespecificado − Qcalculado
El modelo de n ecuaciones linealizadas para el problema de flujo de carga es escrito de la
siguiente manera:
∆P δP δP ∆θ
δθ δV
δQ δQ
∆Q ∆V
δθ δV
Las ecuaciones podrán ser planteadas en forma polar ó rectangular, dependiendo de como
se describa la variable de estado (tensión).
145
Capítulo 5. Solución al problema de FP
flujo de carga desacoplado, flujo de carga óptimo entre otros y en la mayoría de estos se parte
del hecho de que el modelo se plantea en componentes polares.
n
X
Pk + JQk = Vk eJθk Ykm e−Jϕkm Vm e−Jθm (5.64)
m=1
δQk
Deducción de los términos de la diagonal: ( δP
δθk
k
y δθk
)
n
X
Pk + JQk = Vk2 Ykk e−Jϕkk + Vk eJθk Ykm e−Jϕkm Vm−Jθm e−Jθm
m=1
m 6= k
n
δPk δQk X
+J = JVk eJθk Ykm e−Jϕkm Vm e−Jθm
δθk δθk
m=1
m 6= k
n
δPk δQk X
+J = −JVk2 Ykk e−Jϕkk + J(Vk eJθk Ykm e−Jϕkm Vm e−Jθm )
δθk δθk
m=1
| {z }
Pk +JQk
δPk δQk
+J = −JVk2 Ykk e−Jϕkk + J(Pk + JQk )
δθk δθk
δPk δQk
+J = −JVk2 (Gkk − JBkk ) + J(Pk + JQk )
δθk δθk
δPk
= −Vk2 Bkk − Qk (5.65)
δθk
δQk
= −Vk2 Gkk + Pk (5.66)
δθk
δQk
Deducción de los términos de la diagonal: ( δP k
δVk k
V y V )
δVk k
146
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
n
δPk δQk X
+J = 2Vk Ykk eJϕkk + eJθk Ykm e−Jϕkm Vm e−Jθm
δVk δVk
m=1
m 6= k
n
δPk δQk X
Vk + J Vk = 2Vk2 Ykk e−Jϕkk + Vk eJθk Ykm e−Jϕkm Vm e−Jθm
δVk δVk
m=1
m 6= k
n
δPk δQk 2 −Jϕkk Jθk
X
Vk + J Vk = Vk Ykk e + Vk e Ykm e−Jϕkm Vm e−Jθm
δVk δVk
m=1
δPk δQk
Vk + J Vk = Vk2 Ykk e−Jϕkk + Pk + JQk
δVk δVk
δPk δQk
Vk + J Vk = Vk2 (Gkk − JBkk ) + Pk + JQk
δVk δVk
Partiendo en parte real e imaginaria se tiene:
δPk
Vk = Vk2 Gkk + Pk (5.67)
δVk
δQk
Vk = −Vk2 Bkk + Qk (5.68)
δVk
δPk δQk
Deducción de términos fuera de la diagonal: ( δθ m
y δθm
)
n
X
Jθk
Pk + JQk = Vk e Ykm e−Jϕkm Vm e−Jθm
m=1
δPk + JδQk
= −JVk eJθk Ykm e−Jϕkm Vm∗ e−Jθm
δθm
δPk + JδQk
= −J(ek + Jfk )(Gkm − JBkm )(em − Jfm )
δθm
Partiendo en parte real e imaginaria se obtiene:
δPk
= −ek em Bkm − ek fm Gkm + fk em Gkm − fk fm Bkm = Hkm (5.69)
δθm
δQk
= −ek em Gkm + ek fm Bkm − fk em Bkm − fk fm Gkm = Jkm (5.70)
δθm
δPk δQk
Deducción de los términos fuera de la diagonal: ( δV m
Vm y V )
δVm m
δPk δQk
+J = Vk eJθk Ykm e−Jϕkm e−Jθm
δVm δVm
δPk δQk
Vm + JVm = Vk eJθk Ykm e−Jϕkm Vm e−Jθm
δVm δVm
147
Capítulo 5. Solución al problema de FP
El balance de potencia nodal establece que QGi − QDi = Qi . Entonces QGi = Qi + QDi .
Con este valor de potencia reactiva generada, se verifica el cumplimiento de los valores
límites, presentando algunas de las siguintes situaciones:
b) Caso contrario, el nodo cambia a tipo PQ. Identificar el límite violado y especificar en
este nodo los valores de PG y QG .
En la siguiente iteración, son verificados de nuevo los nodos generadores, así como
también, los nodos que al inicio del proceso eran PV y que despues de la verificación fueron
cambiados a nodos PQ.
Para la verificación de estos últimos, se usa el concepto V ≈ KQ (tensión es proporcional
a los reactivos inyectados).
La verificación de estos nodos se inicia con la comparación del Vicalculado y Vicalculado ,
donde i es el nodo a verificar. Como el objetivo es conocer si el nodo que ahora es tipo PQ,
puede regresar a nodo tipo PV y la tensión de este nodo es V esp. , se pueden presentar algunas
de las siguientes situaciones:
esp.
a) Vi = Vicalc. − ∆Vi
148
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
esp.
b) Vi = Vicalc. + ∆Vi
Donde ±∆Vi es la magnitud de la tensión requerida para lograr el V esp. .
La decisión de regresar de nuevo a nodo PV se toma con base en las dos situaciones de
tensión calculadas anteriormente y el límite de generación reactiva fijado, Qmin max
Gi o QGi , así:
Dada la red mostrada en la figura 5.8 de 3 nodos y 3 líneas, una demanda de 720 MW y
410 MVAR y una generación de 800 MW, solucionarlo por el método de Newton Raphson .
1 2
G1
L1
D1 D2
L2 L3
G3 D3
149
Capítulo 5. Solución al problema de FP
Para la solución del problema de flujo de carga usando el método de Newton Raphson se
organiza y cálcula la siguiente información:
Y11 Y12 Y13
Ybus = Y21 Y22 Y23
Y31 Y32 Y33
Condiciones iniciales: V̄10 = 1, 02∠0, V20 = 1, 0∠0, V30 = 1, 02∠0
Errores de potencia activa y reactiva
[∆P |∆Q] = ∆P21 ∆P31 ∆Q12
Después de calcular los errores de balance en el nodo, de potencia activa y reactiva, se
efectuá chequeo por convergencia
Matriz jacobiana
H22 H23 N22
J = H32 H33
N22
J22 J23 L22
Corrección de la variable de estado
Corrección de las variables de estado = ∆θ21 ∆θ31 ∆V21
Se actualizan las variables de estado en la primera iteración:
150
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
NEWTON
MÉTODO GAUSS NEWTON NEWTON DESAC. WARD Y LINEAL
SEIDEL ACOPL. DESAC. RÁPIDO HALE
TENSIONES NODALES
1 1,02 1,02 1,02 1,02 1,02 1,02
2 0,82530 0,82395 0,82376 0,82425 0,82373 1
3 1,02 1,02 1,02 1,02 1,02 1,02
ÁNGULOS [Grados]
1 0 0 0 0 0
2 -14,08927 -14,18044 14,18525 -14,17450 -14,17018 -12,47039
3 -0,45146 -0,48744 -0,48822 -0,48807 -0,48363 -0,07168
GENERACIÓN
POT, ACTIVA [MW]
1 335,81031 337,16553 337,2086 337,14800 336,99700 320
2 0,00 0,00 0,00 0,00 0,00 0,00
3 400,17283 400,02119 400,03910 399,954 400 401,39002
POT, REACT,[Mvar]
1 182,6472 183,78625 183,4653 183,192 183,4840
2 0,00 0,00 0,00 0,00 0,00
3 292,22033 293,53217 293,69550 293,274 293697
DEMANDA
POT. ACTIVA [MW]
1 200 200 200 200 200 200
2 300 300 300 300 300 300
3 220 220 220 220 220 220
FLUJOS DE POT.
POT. ACTIVA [MW]
1-2 126,88280 127,52616 127,5262 127,49602 127,43322 118,60998
1-3 8,92751 9,63937 9,63937 9,65194 9,56409 1,39002
2-3 -179,09535 -179,55640 -179,55640 -179,51834 -179,45680
POT. REACT.[Mvar]
1-2 92,88280 93,24392 93,2439 93,07546 93,35540
1-3 -9,78492 -9,88180 -9,88180 9,88351 9,87158
2-3 -95,814 -96,311 -96,311 -96,17627 -96,461
FLUJOS DE POT.
POT. ACTIVA [MW]
2-1 -119,93251 -120,48967 -120,4897 -120,47064 -120,39685
3-1 -8,91736 -9,62753 -9,62753 -9,64007 -9,55244
3-2 -189,09019 189,64873 189,64873 189,59362 189,55244 181,39002
POT. REACT[Mvar]
2-1 -63,271 -63,474 93,4745 -63,38658 -63,58296
3-1 -7,02733 -7,098880 -7,09880 -7,09688 -7,11030
3-2 149,4276 150,6309 150,6309 150,3711 150,8076
PÉRD. EN LÍNEAS
POT. ACTIVA [MW]
1-2 6,95029 7,03648 7,03648 7,02539 7,03637
1-3 0,01016 0,01184 0,01184 0,01187 0,01165
2-3 9,99484 10,09233 10,09233 10,07528 10,09564
POT. REACT.[Mvar]
1-2 29,16028 29,76946 29,76946 29,68888 29,77244
1-3 -16,99226 16,98060 16,98060 16,98039 16,98188
2-3 53,61364 54,31919 54,31919 54,19485 54,34588
151
Capítulo 5. Solución al problema de FP
Parámetros de la red
(líneas y
transformadores)
Información
nodal
Calcular
∆Pi (i = 1, nd ) ; ∆Qi (i = 1, nd )
Si
Comprobar
, ∆Q ≤ ξ
Max ∆Pque i i
Solución Pare
No
∆P H N ∆θ
∆Q = J L ∆ V / V
−1
∆θ H N ∆P
∆ V / V = J L ∆ Q
∆V k
V k +1 = V k 1 + k
V
θ k +1 = θ k + ∆θ k
k = k +1
Figura 5.9. Algoritmo para la solución del problema de flujo de potencia por Newton Raphson
152
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Ejemplo 5.4.
1 2
3 4
8, 9852 − 44, 8360i −3, 8156 + 19, 0781i −5, 1696 + 25, 8478i 0
−3, 8156 + 19, 0781i 8, 9852 − 44, 8360i 0 −5, 1696 + 25, 8478i
Ybus =
−5, 1696 + 25, 8478i
0 8, 1933 − 40, 8638i −3, 0237 + 15, 1185i
0 −5, 1696 + 25, 8478i −3, 0237 + 15, 1185i 8, 1933 − 40, 8638i
Condiciones iniciales:
Primera iteración
153
Capítulo 5. Solución al problema de FP
Errores de potencia
[∆P |∆Q] = −1, 5966 −1, 9395 2, 2129 −0, 4465 −0, 8345
45, 4429 0 −26, 3648 8, 8818 0
0 41, 2687 −15, 4209 0 8, 1328
Jacobiano =
−26, 3648 −15, 4209 41, 7857 −5, 2730 −3, 0842
−9, 0886 0 5, 2730 44, 2290 0
0 −8, 2537 3, 0842 0 40, 4590
V1 V2 V3 V4
V = 1, 0000 0, 9834 0, 9710 1, 0200
Anguloso = 0 −0, 9309 −1, 7879 1, 5438
Segunda iteración
Errores de potencia
[∆P |∆Q] = −0, 0323 −0, 0645 0, 0359 −0, 0342 −0, 0620
44, 3749 0 −25, 6778 7, 1397 0
0 39, 7018 −14, 7736 0 5, 9619
Jacobiano =
−26, 1256 −15, 1217 41, 2473 −4, 1296 −2, 1828
−10, 3562 0 6, 2998 43, 0530 0
0 −9, 6597 3, 8597 0 38, 4642
44, 3270 0 −25, 6508 7, 0969 0
0 39, 6096 −14, 7398 0 5, 8755
Jacobiano =
−26, 1026 −15, 0938 41, 1963 −4, 1183 −2, 1655
−10, 33721 0 6, 3048 42, 9755 0
0 −9, 6932 3, 8683 0 38, 3186
154
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
V1 V2 V3 V4
|V | = 1, 0000 0, 9824 0, 9690 1, 0200
Anguloso = 0 −0, 9760 −1, 8720 1, 5231
Resultados obtenidos con ecuaciones del subproblema 2
Para el cálculo se usa una potencia base 100 MVA
Potencias activas inyectadas:
P1 = 136, 7948 Mw
P2 = −169, 9966 Mw
P3 = −199, 9854 Mw
P4 = 237, 9958 Mw
Potencias Reactivas inyectadas:
Q1 = 83, 4977 MVAR
Q2 = −105, 3459 MVAR
Q3 = −123, 9235 MVAR
Q4 = 131, 8393 MVAR
Flujos de Potencia:
P12 = 38,6883 Mw Q12 = 22,2969 Mvar
P13 = 98,1065 Mw Q13 = 61,2008 Mvar
P24 = -131,5350 Mw Q24 = -74,1109 Mvar
P34 = -102,9101 Mw Q34 = -60,3650 Mvar
P21 = -38,4616 Mw Q21 = -31,2349 Mvar
P31 = -97,0753 Mw Q31 = -63,5585 Mvar
P42 = 133,2504 Mw Q42 = 74,9166 Mvar
P43 = 104,7453 Mw Q43 = 56,9227 Mvar
155
Capítulo 5. Solución al problema de FP
i j
1:t
+ +
+
Vi t Vi Vj
- - -
i j
ty
Vi t (t − 1) y (1 − t ) y Vj
156
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
tnuevo
i = tviejo
i + ci (Vireg − Vi ).(t valor del tap) (5.79)
Vs
VS Vs ×VpBase
La relación de transformación es: tip.u. = Base
Vp = Vp ×VsBase
Vp
Base
Vireg (Tensión que se desea mantener), Vi (tensión calculada)
ci = [−0.5 ⇔ −1, 0] rango de valores típicos empleados.
En cada iteración se calcula un nuevo valor de t el cual se reemplaza en el modelo Π del
transformador para posteriormente incluirlo en la matriz Ybus , la que posteriormente se emplea
en el cálculo de los errores de las potencias inyectadas y de la matriz jacobiana. El flujo de
carga converge cuando se cumple el criterio de convergencia.
El tap podrá ser controlado automáticamente en el flujo de carga y para lograrlo se requiere
incluir la variable tap en el modelo del problema de flujo de carga. Lo anterior se logra
efectuando en el nodo en el que se desea controlar la tensión un cambio de variable de estado,
así, la corrección de la tensión asociada al nodo donde se desea controlar, se sustituye por la
variable corrección de tap. Además de lo anterior, en la matriz jacobiana la columna asociada
a la tensión donde esta será controlada, se sustituye por la derivada de las potencias inyectadas
activa y reactiva, respecto al tap. En la figura 5.13 se presenta un sistema de cinco nodos que
contiene 2 transformadores reguladores.
1 2 3
6 4 5
El modelo matemático del flujo de carga incluidos los dos transformadores de la figura
5.13 se presenta a seguir. En este el nodo 1 es tipo P V y el nodo 6 es tipo V θ
157
Capítulo 5. Solución al problema de FP
θ1 θ2 θ3 θ4 θ5 t2 t3 V4 V5
δ∆P1 δ∆P1 δ∆P1
∆P1 0 0 0 0 0 0
δθ1 δθ 4 δ V4
δ∆P2 δ∆P2 δ∆P2 δ∆P2
∆P 0 0 0 0 0
2 δθ 2 δθ 4 δ t2 δ V4
δ∆P3 δ∆P3 δ∆P3 δ∆P3
∆P 0 0 0 0 0
3 δθ3 δθ5 δ t3 δ V5
δ∆P4
∆P δ∆P4 δ∆P4
0
δ∆P4 δ∆P4 δ∆P4
0
δ∆P4
4 δθ1 δθ 2 δθ 4 δθ5 δ t2 δ V4 δ V5
∆P = 0 δ∆P5 δ∆P5 δ∆P5 δ∆P5 δ∆P5 δ∆P5
0 0
5 δθ3 δθ 4 δθ5 δ t3 δ V4 δ V5
∆Q2 0 δ∆Q2 δ∆Q2 δ∆Q2 δ∆Q2
0 0 0 0
δθ 2 δθ 4 δ t2 δ V4
∆Q3 0 δ∆Q3 δ∆Q3 δ∆Q3 δ∆Q3
0 0 0 0
δθ3 δθ5 δ t3 δ V5
∆Q4 δ∆Q4 δ∆Q4 δ∆Q4 δ∆Q4 δ∆Q4 δ∆Q4 δ∆Q4
0 0
δθ δθ 2 δθ 4 δθ5 δ t2 δ V4 δ V5
1
∆Q5 δ∆Q5 δ∆Q5 δ∆Q5 δ∆Q5 δ∆Q5 δ∆Q5
0 0
δθ3 δθ 4 δθ5
0
δ t3 δ V4 δ V5
Variables: θ1 θ2 θ3 θ4 θ5 t2 t3 V4 V5
Los términos de la matriz Jacobiana que tienen relación con los nodos en los cuales se
controla la tensión usando los taps de los transformadores se calculan así:
Los términos de la matriz jacobiana que relacionan la potencia activa y reactiva con la
tensión en magnitud y ángulo son:
158
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• La diferencia angular en radianes θik entre nodos adyacentes es pequenã por lo tanto:
Senθik ≃ θik ; Cosθik ≃ 1, 0.
• Submatriz J: las conductancias Gik son multiplicadas por valores de coseno y las
susceptancias Bik por valores seno. Por lo tanto, se obtienen valores relativamente
pequeños.
Con este método se obtienen los mismos resultados respecto al acoplado. A pesar de que el
número de iteraciones es mayor, estas ocurren en un menor tiempo de computo. Las matrices
H y L son aproximadamente la cuarta parte del tamaño de la matriz jacobiana del método
acoplado. Por lo tanto la inversa de estas matrices requiere de menos tiempo de cálculo respecto
de la matriz jacobiana del método acoplado.
Características de las matrices H y L:
• Están en función de las tensiones nodales, por lo cual cambian en cada iteración
159
Capítulo 5. Solución al problema de FP
δ∆Q(θk , V k )
∆Q(θk , ∆θk , V k ) = ∆Q(θk , V k ) + .∆θk
δθ
k k
Se observa que aparece el término: δ∆Q(θ
δθ k
,V )
, que a pesar de no ser la matriz J ( δQ
δθ
), tiene
parte de la información contenida en esta. Esto quiere decir que el problema preserva en parte
el acoplamiento del subproblema Q − θ, contenido en la matriz J.
Ejemplo 5.5.
Resolver el sistema de dos barras mostrado en la figura 5.15 por el método flujo de carga
desacoplado de Newton.
Datos del sistema
Barra 1: (V-θ)
Entonces: V1 = 1, 0 p.u θ1 = 0, 0 : incógnitas : P1 , Q1
Barra 2: (P-Q)
Entonces: P2 = −0, 30 p.u. Q2 = 0, 07 p.u. incógnitas : V2 , θ2
Tolerancia de convergencia: ξ = 0, 002 p.u.
Parámetros del sistema (línea de transmisión):
R1−2 = 0,2 p.u. , X1−2 = 1,0 p.u. , y/2 = 0,02 p.u.
(Sbase = 100 MVA)
Expresiones requeridas en el cálculo
P2 = 0, 1923V22 − 0, 1923V2 Cosθ2 + 0, 9615V2 Senθ2
Q2 = 0, 9415V22 − 0, 1923V2 Senθ2 − 0, 9615V2 Cosθ2
∆P2 = −0, 30 − P2
∆Q2 = 0, 07 − Q2
H22 = −Q2 + 0, 9415V22
L22 = V12 [Q2 + 0, 9415V22 ]
Proceso iterativo de solución
• Condiciones iniciales
160
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Información nodal
Calcular ∆ PN = PG − PD − PN (V , Y )
Comprobar si Comprobar si
∆PN < ξ ∆QN < ξ
no no
∆PN = [ H ]∆θ
k k
θ k +1 = θ k + ∆θ k
Solución
k1 = k1 + 1
Calcular ∆QN = QG − QD − QN (V , Y )
Comprobar Comprobar
∆QN < ξ si ∆PN < ξ si
no
k
kno∆ V
∆QN = [ L ] no
Vk
V k +1 = V k + ∆V k
k2 = k2 +1
Figura 5.14. Algoritmo para la solución del problema de flujo de potencia por Newton
desacoplado
161
Capítulo 5. Solución al problema de FP
k m
ikm Zkm = Rkm + J Xkm imk
Vk ∠θ k Vm ∠θ m
sh sh
y/ 2= Jbkm y/ 2= Jbkm
162
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Tercera iteración (P − θ)
ii) P2 (V22 , θ22 ) = 0, 1923(V22 )2 − 0, 1923V22 Cosθ22 + 0, 9615V22 Senθ22
P2 (V22 , θ22 ) = −0, 29887
∆P2 (V22 , θ22 ) = −0, 30 − (−0, 29887) = −0, 00113
163
Capítulo 5. Solución al problema de FP
Observación:
Si fuera efectuada una iteración adicional, los valores obtenidos de θ2 y V2 son los
siguientes:
θ2 = −0, 33186 (rad.)
V2 = 0, 97522 (p.u.)
• Pequeña diferencia angular expresada en radianes entre nodos adyacentes, por lo tanto:
Sen(θk − θj ) ≈ 0
Cos(θk − θj ) ≈ 1, 0
164
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Valores de las tensiones en los nodos del sistema son próximas a 1 p.u.
′
Hkm = Vm (Gkm Senθkm − Bkm Cosθkm ) (5.86)
′ Qk
Hkk = − − Vk Bkk (5.87)
Vk
′
Lkm = Gkm Senθkm − Bkm Cosθkm (5.88)
′ Qk
Lkk = 2 − Bkk (5.89)
Vk
De esta forma el modelo del problema de flujo de carga desacoplado rápido asume la
siguiente forma:
∆P ′
= [H ]∆θ (5.90)
V
∆Q ′
= [L ]∆V (5.91)
V
Con base en las características propias del sistema de transmisión se asumen las siguientes
aproximaciones:
• Cosθkm ≈ 1; Senθkm ≈ 0
165
Capítulo 5. Solución al problema de FP
Llegando así a matrices constantes que dependen solamente de los parámetros del sistema.
′ ′ ′ ′′
La estructura de las matrices [H ] = [B ] y [L ] = [B ] es semejante a la de la matriz Ybus , con
las siguientes diferencias:
′ ′
• B , es la parte imaginaria de la matriz Ybus , multiplicada por un signo menos. En B no
aparece la línea y columna de la barra V θ
′′ ′
• B , al igual que B , es la parte imaginaria de la matriz Ybus , multiplicada por un signo
′′
menos. En B ademas de no aparecer las líneas y columnas de las barras V θ, tampoco
aparecen las de las barras PV.
Finalmente, las ecuaciones del método desacoplado rápido son escritas de la siguiente forma:
∆P ′
= [B ]∆θ (5.98)
V
∆Q ′′
= [B ]∆V (5.99)
V
′ ′′
[B ] y [B ] son matrices reales, dispersas y tiene la misma estructura de [H] y [L]
respectivamente. Por ser matrices simétricas, solamente la triangular superior es almacenada.
′ ′′
Algunas consideraciones de las matrices B y B son:
′
• en [B ] se omiten efectos de flujos MVAR, por ejemplo reactancia shunt y el modelo del
transformador en fase para control de tensión.
′′
• de [B ] se omiten efectos del trsnformador desfasador
El método desacoplado rápido presenta un mejor desempeño cuando son despreciadas las
resistencias serie en la formación de B ó de B , entonces bkm = X1km , obteniendo las versiones:
′ ′′
([B],[B]), ([X],[B]) y ([B],[X]). En la primera se lleva en cuenta la resistencia serie de las líneas
′ ′′
de transmisión, tanto en la fomación de [B ] como en [B ], en la segunda la resistencia se
166
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
′ ′′
desprecia en la fomación de [B ] y en la tercera se desprecia en la formación de [B ]. De esta
forma las matrices se conforman de la siguinte manera:
Si se desprecia la resistencia serie de las líneas de transmisión
1
Bkm = − Xkm
X 1
B⇒
Bkk = Xkm
mΩk
167
Capítulo 5. Solución al problema de FP
Información nodal
Calcular ∆ PN = PG − PD − PN (V , Y )
Comprobar si Comprobar si
∆PN < ξ ∆QN < ξ
no
no
∆ PN
= [ B ′ ]∆ θ k
Vk
θ k +1 = θ k + ∆θ k
Solución
k1 = k1 + 1
Calcular ∆QN = QG − QD − QN (V , Y )
Comprobar si Comprobar si
∆QN < ξ ∆PN < ξ
no
no
∆QN
= [ B′′]∆V k
Vk
V k +1 = V k + ∆V k
k2 = k2 +1
Figura 5.16. Algoritmo para la solución del problema de flujo de potencia por Newton
desacoplado rápido
168
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
′′
- Dimensión de B : NPQ
• Condiciones iniciales
• Primera iteración (P − θ)
ii) P2 (V20 , θ20 ) = 0, 1923(V20 )2 − 0, 1923V20 Cosθ20 + 0, 9615V20 Senθ20 P2 (V20 , θ20 ) =
0, 0 ⇒ ∆P2 (V20 , θ20 ) = −0, 30 − 0 = −0, 30
iii) |∆P2 (V20 , θ20 )| < ξ ? [NO]
∆P2 (V20 ,θ20 ) ′ −0,30
iv) V20
= B ∆θ20 ⇒ ∆θ20 = 1∗1
= −0, 30
v) θ21 = θ20 + ∆θ20 = 0 − 0, 30 = −0, 30
θ21 = −0, 30
vi) p = 0 + 1 = 1
vii) Kq = 1 ⇒ ir al paso (x)
• Segunda iteración (P − θ)
169
Capítulo 5. Solución al problema de FP
vii) kq = 1 ⇒ ira(x)
• Tercera iteración (P − θ)
170
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Observación:
El error es reducido después de una iteración adicional se obtiene:
θ2 = −0, 33193
V2 = 0, 97518
Entonces:
" #
′′ ′
X X
Bkk = Lkk = − bkm − bsh
k + bsh
km (5.100)
m∈Ωk m∈Ωk
171
Capítulo 5. Solución al problema de FP
′′ Qk
• Análisis Modificado llevando en cuenta que el término de la diagonal Bkk = −Bkk + Vk2
X
Qk = V k Vm (Gkm Senθkm − Bkm Cosθkm )
m∈k
X
Qk = −Bkk Vk2 + Vk Vm (Gkm Senθkm − Bkm Cosθkm )
m∈Ωk
′ Qk
Lkk = −Bkk +
Vk2
" !#
P P
L′ kk = −2 bkm + bsh k + bsh
km +
P m∈Ωk m∈Ωk (5.101)
1
Vk
Vm (−gkm Senθkm + bkm Cosθkm )
m∈Ωk
" #
P
L′ kk = −2 bsh k + bsh
km −
m∈Ωk (5.102)
P n o
bkm + bkm 1 − VVmk Cosθkm + VVmk gkm Senθkm
m∈Ωk
• Vk ≈ Vm ≈ 1
′
Entonces los términos de Lkk se reducen a:
172
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
" #
′
X X
Lkk = −2 bsh
k + bsh
km − bkm (5.103)
m∈Ωk m∈Ωk
" #
′′ ′
X X
Bkk = Lkk = −2 bsh
k + bsh
km − bkm (5.104)
m∈Ωk m∈Ωk
Polares:
k k k
∆PN i δPN i /δθi δPN i /δVi ∆θi
= (5.106)
∆QN i δQN i /δθi δθN i /δVi ∆Vi
173
Capítulo 5. Solución al problema de FP
Información nodal
k =0 (Contador de iteraciones)
k = k +1
i = i +1
δPt δPt
δθ δVt
t
δQt δQt
δθ t δVt
θ t k +1 = θ t k + ∆θ t k , Vt = Vt + ∆Vt
k +1 k k
Tensión actualizada:
Si
i < nodos
No
ep p ξ y eq p ξ Resultados complementarios
No
Figura 5.17. Algoritmo para la solución del problema de flujo de potencia Ward y Hale
174
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
slack no se plantean ecuaciones, para los nodos P V , solo se plantea la ecuación (P − θ), con
la que se resuelte el θ en dicho nodo y para los nodos de carga se plantean las 2 ecuaciones que
resuelven la variable de estado (V, θ) en dicho nodo.
Lkk = Q
′
k
Vk2
− Bkk
y los fuera de la diagonal
′
Nkm = Cosθkm Gkm + Senθkm Bkm
′
Lkm = −Cosθkm Bkm + Senθkm Gkm
′
Hkm = −Vm Cosθkm Bkm + Vm Senθkm Gkm
′
Jkm = −Vm Cosθkm Gkm − Vm Senθkm Bkm
k
′
Para el caso de Lkk , subproblema Q − V , no aplica este tipo de simplificación y por lo
tanto es determinado el término Lkk = Q
′
k
V2
− Bkk
k
175
Capítulo 5. Solución al problema de FP
Qk 1 X
= −B kk + Vm (−gkm Senθkm + bkm Cosθkm )
Vk2 Vk m∈Ω
k
Qk X X
=(− b km − bsh sh
km − bk )+
Vk2 m∈Ωk m∈Ωk
1 X
Vm (−gkm Senθkm + bkm Cosθkm )
Vk m∈Ω
k
Qk X Vm
2
= − bkm (1 − bkm Cosθkm )+
Vk m∈Ωk
Vk
1 X X
Vm (−gkm Senθkm ) − bsh sh
km − bk
Vk m∈Ω m∈Ω
k k
Qk X Vm X
= − b km (1 − b km ) − bsh sh
km − bk
Vk2 m∈Ω
V k
m∈Ω
k k
Finalmente se obtiene:
Qk X
= − bsh sh
km − bk
Vk2 m∈Ω k
Por lo tanto:
′
X
Lkk = Qk /Vk2 − Bkk = −( bsh sh
km + bk ) − Bkk
m∈ Ωk
176
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
xtemp. = [A−1 ]a
• Paso 2: determinar y
−1
y = [Dequiv. ]bequiv.
• Paso 3: determinar x
177
Capítulo 5. Solución al problema de FP
∆P2 −b g ∆θ2
=
∆Q2 −g −b ∆V2
[Link] Version BX
Para la solución se aplica la eliminación gaussiana en bloque, en la cual el primer bloque de
ecuaciones es reemplazado en el segundo, de la siguiente manera:
∆P2 −b g ∆θ2
g = g2
∆Q2 - b ∆P2 0 −b − b
∆V2
Donde:
g δ∆Q2
∆Q2 − ∆P2 = ∆Q2 (V2 , θ2 + ∆θ2 ) = ∆Q2 (V2 , θ2 ) + × ∆θ2
b δθ2
g g
∆Q2 − ∆P2 = ∆Q2 (V2 , θ2 ) − g × ∆θ2 = ∆Q2 (V2 , θ2 ) − × ∆P2
b b
El término del jacobiano:
g2 −x R 2 + x2 R2
−b − =− 2 − ×
b R + x2 −x (R2 + x2 )2
g2 x R2 (x2 + R2 ) 1
−b − = 2 + = =
b R + x2 (R2 + x2 )x (R2 + x2 )x x
El modelo anterior es entonces reescrito de la siguiente manera:
∆P2 −b g ∆θ2
= 1
∆Q2equiv. 0 x
∆V2
Este modelo que es equivalente al anterior es resuelto usando un algoritmo de 3 pasos así:
• Paso 1:
∆P2 = −b∆θ2 + g∆V2
∆θ2 = −b−1 ∆P2 + b−1 g∆V2
Temporalmente se desprecia el segundo término y se calcula: ∆θ2temp. = −b−1 ∆P2
178
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
• Paso2:
Calcular: ∆V2 = x∆Q2equiv. = x(∆Q2 − gb−1 ∆P2 ) = x∆Q2 (V2 , θ2 + ∆θ2temp. )
• Paso 3:
Calcular:
∆θ2 = −b−1 ∆P2 + b−1 g∆V2
∆θ2 = ∆θ2temp. + b−1 g∆V2
Esta versión que es denominada BX, puede ser resuelta en 2 o 3 pasos, dependiendo de
las caracteristicas de la red eléctrica estudiada.
En redes con R << x, es decir cuando los valores de la conductancia (g) son pequeñas en
relación a la susceptancia (b), caracterísitico de las redes de transmisión de alta tensión, puede
ser resuelto usando un algoritmo que usa los dos primeros pasos y la corrección obtenida con el
tercer paso sería entonces poco significativa.
En redes de distribución, donde el valor de la resistencia es comparable o mayor que el de
la reactancia, quiere decir, el valor de la conductancia es mayor o igual que la susceptancia, se
recomienda usar un algoritmo que incluya los 3 pasos descritos anteriormente.
[Link] Version XB
Como en la versión BX, se aplica un procedimiento de eliminación gaussiana en bloque. En este
procedimiento de eliminación, el segundo bloque de ecuaciones es remplazado en el primero,
resultando el siguiente modelo equivalente.
2
∆P2 + gb ∆Q2 −b − gb 0 ∆θ2
=
∆Q2 −g −b ∆V2
1
∆P2equiv. 0 ∆θ2
= x
∆Q2 −g −b ∆V2
Esta versión es resuelta usando el siguiente algoritmo de 3 pasos:
• Paso 1:
∆Q2 = −g∆θ2 - b∆V2
∆V2 = − gb ∆θ2 - b−1 ∆Q2
Temporalmente se desprecia el primer término y se calcula: ∆V2temp. = −b−1 ∆Q2
• Paso2:
Calcular: ∆θ2 = x∆P2equiv. = x(∆P2 + gb−1 ∆Q2 ) = x∆P2 (V2 + ∆V2temp. , θ2 )
• Paso 3:
Calcular:
∆V2 = − gb ∆θ2 - b−1 ∆Q2
∆V2 = − gb ∆θ2 + ∆V2temp.
Como en el caso anterior, esta versión puede ser resuelta en 2 o 3 pasos, dependiendo de
las características de la red.
179
Capítulo 5. Solución al problema de FP
180
FLUJO DE CARGA LINEALIZADO
6
6.1 Introducción
El flujo de carga linealizado, denominado también flujo de carga DC o simplemente FCDC, es
una forma aproximada de resolver ecuaciones de flujo de carga no lineal. En la deducción del
flujo de carga DC se considera la estrecha relación existente entre el flujo de potencia activa
en las líneas de transmisión y transformadores, y la abertura angular entre los extremos de las
líneas o los transformadores.
Pkm ⇒ δ = θk − θm
• La aproximación es válida solamente para sistemas de transmisión; esto es, para niveles
elevados de tensión. La aproximación es mejor cuanto mayor fuese el nivel de tensión.
181
Capítulo 6. Flujo de carga linealizado
• El modelo de flujo de carga C.C. no tiene en cuenta las magnitudes de las tensiones
nodales y el flujo de potencia reactiva. Desde este punto de vista asume que el problema
reactivo ya está resuelto.
Vk ≈ Vm ≈ 1p.u. (6.5)
Senθkm ≈ θkm (en radianes) (6.6)
−xkm 1
bkm = 2 2
≈− (6.7)
rkm + xkm xkm
θkm θk − θm
Pkm = = (6.8)
xkm xkm
La ecuación 6.8 tiene la forma de la Ley de Ohm aplicada a una resistencia, tal como se
presenta en la figura 6.1
El flujo Pkm de un Transformador en fase es:
182
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
ikm Pkm
Vk θk
Vk − Vm θk − θm
ikm =
rkm
rkm ↔ Pkm =
X km
X km
Vm θm
Figura 6.1. Comparación de modelos de corriente continua y D.C.
183
Capítulo 6. Flujo de carga linealizado
Pkm tiene entonces dos componentes, una variable que depende del estado del sistema
θkm /xkm y otra fija e igual a ϕkm /xkm
θk θm
k θp =θk +ψkm m
Pkm =
(θkm+ψkm)
Xkm
Esquema equivalente
θk θm
k m
Xkm
ψ km ψ km
−
X km X km
Figura 6.2. Modelo del Transformador Desfasador Puro
Donde la parte del flujo invariante ϕkm /xkm aparece como una carga adicional en la barra
k, y una generación adicional en la barra m.
184
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
k
PGk Pk1 1
Pk 2 2
Pkn ⋮ n
⋮
PDk
Figura 6.3. Ecuación de balance de potencia nodal
θk θk θk θk
Pk = − xθk1
1
− θ2
xk2
− θ3
xk3
− ··· − θn
xkn
+ xk1
+ xk2
+ xk3
+ ··· + xkn
Pk = − xθk1
1
− θ2
xk2
··· + 1
xk1
+ 1
xk2
+ ··· + xkn
1
θk · · · − θn
xkn
θk
Pk = − xθk1
1
− θ2
xk2
··· + xkk
··· − θn
xkn
185
Capítulo 6. Flujo de carga linealizado
Siendo que:
1
′
Bkm = −x−1
km = − (6.21)
xkm
′
X X 1
Bkk = x−1
km = (6.22)
m∈Ωk m∈Ωk
xkm
′
La matriz B en 6.20 es singular y se debe eliminar una de las ecuaciones y adoptar para
la barra correspondiente como referencia angular θk = 0
La matriz es singular, ya que las pérdidas son despreciadas y la suma de las componentes
de P es nula; esto es, la inyección de potencia en una barra cualquiera puede ser obtenida a
partir de la suma algebraica de las demás. Así, ahora se transforma en un sistema no singular
de dimensión nb − 1. La solución de este sistema algebraico entrega los ángulos de las nb − 1
barras restantes.
Si el sistema eléctrico tiene transformadores en fase o desfasadores, entonces el sistema
de ecuaciones 6.20 continúa válido con las siguientes observaciones:
′
• La formación de la matriz B cuando existen transformadores en fase o desfasadores, es
realizada exactamente igual que en el caso de las líneas de transmisión.
186
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Ejemplo 6.1.
Dado el siguiente sistema de Tres (3) Barras, calcular los flujos de potencia activa por las
líneas.
1 2
P13 P23
P3 = -1.0
θ13 0−(−3/8) 3
P13 = x13
= 1/2
= 4
⇒ P13 = 0, 75
θ23 −1/4−(−3/8) 1
P23 = x23
= 1/2
= 4
⇒ P23 = 0, 25
187
Capítulo 6. Flujo de carga linealizado
X X
Pk = (1 − Cosθkm )gkm + x−1
km Senθkm (6.28)
m∈Ωk m∈Ωk
188
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Comparando las ecuaciones 6.29 y 6.18, se observa que la diferencia está en el término:
1X 2
gkm θkm (6.30)
2 m∈Ω
k
Que representa las pérdidas aproximadas de las líneas de transmisión que salen de la
barra k. Este hecho puede ser verificado analizando las pérdidas de transmisión en la línea k-m.
Ejemplo 6.2.
• Resolver el sistema 6.33 con los valores de pérdida obtenidas en el paso anterior. Con
los θ conocidos, se calculan los flujos en las líneas.
189
Capítulo 6. Flujo de carga linealizado
1
15
g12 = 2
15
2 = 26
( 151 ) +( 13 )
1
5
g13 = 1 2
10
2 = 13
( ) ( 12 )
10
+
1
5
g23 = 1 2
10
1 2
= 13
( ) ( )
10
+ 2
2
15 1
2
Pe12 = g12 θ12 = 26
0− − = 26 15
× 26
= 0, 036
4
2
5 3
2
Pe13 = g13 θ13 = 13
0− − = 135 × 9
× 64
= 0, 054
8
2
5 1 3
2
Pe23 = g23 θ23 = 13
−4 − − = 13 ×5 64 = 0, 006
8
1 2
P1 = 1.5 P2 = -0.5
-0.045 -0.021
P3 = -1.0 -0.030
190
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
Ejemplo 6.3.
En el sistema de 3 barras anterior, sustituir la línea 1-2 por un desfasador con las siguientes
características:
ϕ12 = 5◦ , x12 = 1/3
Recalcular los flujos en las líneas: ϕ12 = 5◦ = 0, 0873rad
ϕ12 0, 0873
P12 = = 1 = 0, 2618
x12 3
1 2
-0.2618 0.2618
P3 = -1.0
′ θ12 0,1846
P12 = X12
= 1 = 0, 5538 p.u.
3
′ ϕ12
P12 = P12 + X12
= 0, 5538 + 0, 2618 = 0, 8156 p.u.
θ13 0,3423
P13 = X13
= 1 = 0, 6846 p.u.
2
191
Capítulo 6. Flujo de carga linealizado
θ23 −0,1846−(−0,3423)
P23 = X23
= 1 = 0, 3154 p.u.
2
1 2
0.8156
P1 = 1.5 P2 = -0.5
0.6844 0.3156
P3 = -1.0
Figura 6.7. Sistema de 3 nodos mostrando los flujos (incluido el transformador desfasador)
192
Bibliografía
Abadie, J. and Carpentier, J.: (1969). Generalization of the wolf reduced gradient method
to the case of nonlinear constraints., Recent Advances in Mathematical Programming. in
Optimization, R. Fletcher (ed), New York: Academic Press. pp. 37–47.
Bazaraa, M., Sherali, H. and Shetty, C.: 1993, Nonlinear Programming: Theory and
Algorithms, John Wiley & Sons.
Calderón, L. F. and Barón, A.: 1988, Introducción al análisis de sistemas de potencia en estado
estacionario, Publicaciones facultad de Ingeniería Universidad Nacional Bogotá.
Castillo, E., Conejo, A. J., Pedregal, P., García, R. and Alguacil, N.: 2002, Building and Solving
Mathematical Programming Models in Engineering and Science, John Wiley & Sons.
Dommel, H. and Tinney, W.: (1968). Optimal power flow solutions, IEEE Trans. Power
Apparat. Syst. 87, 1866–1876.
Granada, M., Rider, M., Mantovani, J. and Shahidehpour, M.: (2008). Multi-areas optimal
reactive power flow.
Iwamoto, S. and Tamura, Y.: (1981). A load flow calculation method for ill-conditioned power
systems, IEEE Trans. Power App. Syst. 100(4), 1736–1743.
Kindermann, G.: 2003, Curto-Circuito, Florianópolis: Edição do Autor.
UFSCEEL-LABPLAN.
Lebow, W. M., Mehra, R. K., Nadira, R., Rouhani, R. and Usoro, P.: (1984). Optimization of
reactive volt-amperes (var) sources in system planning, EPRI Report El-3729, Project 2109-1
1.
Monticelli, A.: 1999, State Estimation in Electric Power Systems a Generalized Approach,
Kluwer Academic Publishers. Massachusetts, USA.
Monticelli, A. and García, A.: 1999, Introducao a sistemas de energía eléctrica, editora da
Unicamp, Brasil.
Monticelli, A. J.: 1983, Fluxo de Carga em Redes de Energia Elétrica, Editora Edgard Blücher
LTDA.
Ness, J. E. V.: (1959). Iteration methods for digital load flow studies, IEEE Transactions on
Power Apparatus and Systems 78, 583–588.
Renegar, J.: (1988). A polynomial-time algorithm, based on newton method, for linear
programming, Math. Program 40, 59–93.
Rosenbrock, H.: (1960). Automatic method for finding the greatest or least value of a function,
Computer J pp. 175–184.
Sasson, A., Trevino, C. and Aboytes, F.: (1971). Improved newton’s load flow through a
minimization technique, IEEE Trans. Power App. Syst. 90, 1974–1981.
Sun, D., Ashley, B., Beuler, B., Hughes, A. and Tinney, W.: (1984). Optimal power flow by
newton approach, IEEE Transactions on Power Apparatus and Systems 103(10), 2864–2880.
193
Bibliografía
Sun, D. I., Demaree, K. D. and Brewer, B.: (1990). Application and adaptation of newton
for optimal power flow., In Application of Optimization Methods for Economy/Security
Functions in Power System Operations. An IEEE Tutorial pp. 14–90.
Tinney, W. F. and Walker, J. W.: (1967). Direct solutions of sparse network equations by
optimally ordered triangular factorizations., Proceedings of the IEEE 55, 1801–1809.
Tinney, W. and Hart, C.: (1967). Power flow solution by newton’s method, IEEE Transactions
on Power Apparatus and Systems 86, 1449–1456.
Ward, J. and Hale, W.: (1956). Digital computer solution of power flow problems, AIEE
Transactions Power App. Syst. 75, 398–404.
Wolfe, P.: (1963). Methods of nonlinear programming, Recent Advances in Mathematical
Programming. in R.L. Graves and P. Wolfe (eds), McGraw-Hill .
Zollenkopf, K.: (1971). Bi-factorization-basic computation algorithm and programming
techniques, Large sparse sets of linear equations, edited by Reid, J.K., N. York, Academic
Press pp. 75–97.
194
Índice general
195
Índice general
196
Ramón A. Gallego R., Antonio H. Escobar Z., Mauricio Granada E.
197
Índice general
de Control, 18
de Estado, 18
dependientes, 18
No controladas, 17
Nodales, 17
198









