ESTABILIDAD ANGULAR DE UN SISTEMA ELÉCTRICO DE POTENCIA

CONSIDERANDO LA INTEGRACIÓN DE REGULADORES DE VELOCIDAD

Y VOLTAJE
 UNIVERSIDAD POLITÉCNICA SALESIANA
SEDE QUITO

CARRERA:
INGENIERÍA ELÉCTRICA

Trabajo de titulación previo a la obtención del título de

INGENIERO ELÉCTRICO

TEMA:
ESTABILIDAD ANGULAR DE UN SISTEMA ELÉCTRICO DE POTENCIA
CONSIDERANDO LA INTEGRACIÓN DE REGULADORES DE VELOCIDAD
Y VOLTAJE

AUTOR:
DANIELA ANABEL ANTAMBA ÁVILA

DIRECTOR:
CARLOS ANDRÉS BARRERA SINGAÑA

Quito, Noviembre 2018

 . Datos de Catalogación Bibliográfica .
Daniela Anabel Antamba Ávila

ESTABILIDAD ANGULAR DE UN SISTEMA ELÉCTRICO DE POTENCIA

CONSIDERANDO LA INTEGRACIÓN DE REGULADORES DE VELOCIDAD
Y VOLTAJE

Universidad Politécnica Salesiana, Quito – Ecuador 2018

Ingeniería Eléctrica

Breve reseña histórica e información de contacto.

Daniela Anabel Antamba Ávila (Y’1995 – M’08). Realizó sus

estudios de nivel secundario en el Colegio “República del Ecuador” de
la ciudad de Otavalo. Egresado de Ingeniería Eléctrica de la
Universidad Politécnica Salesiana. Su trabajo se basa en el estudio de
estabilidad angular en Sistemas Eléctricos de Potencia.
dantamba@[Link]

Dirigido por:

Carlos Andrés Barrera Singaña (Y’1989). Realizó sus estudios

superiores en la Escuela Politécnica Nacional de Quito, donde se
graduó de Ingeniero Eléctrico en octubre de 2013. Además, cursó
estudios de posgrado en la Universidad de Zaragoza, España, donde
obtuvo el título en Máster en Energías Renovables y Eficiencia
Energética con Especialidad en Sistemas Eléctricos. Actualmente es
profesor ocasional a tiempo completo de la Universidad Politécnica
Salesiana. cbarrera@[Link]

Todos los derechos reservados:

Queda prohibida, salvo excepción prevista en la ley, cualquier forma de reproducción,

distribución, comunicación pública y transformación de esta obra para fines comerciales,
sin contar con la autorización de los titulares de propiedad intelectual. La infracción de
los derechos mencionados puede ser constitutiva de delito contra la propiedad
intelectual. Se permite la libre difusión de este texto con fines académicos o
investigativos por cualquier medio, con la debida notificación a los autores.

DERECHOS RESERVADOS
©2018 Universidad Politécnica Salesiana
QUITO - ECUADOR

I
 DECLARATORIA DE COAUTORÍA DEL DOCENTE TUTOR

Yo, Carlos Andrés Barrera Singaña declaro que bajo mi dirección y asesoría fue
desarrollado el trabajo de titulación “ESTABILIDAD ANGULAR DE UN SISTEMA
ELÉCTRICO DE POTENCIA CONSIDERANDO LA INTEGRACIÓN DE
REGULADORES DE VELOCIDAD Y VOLTAJE” realizado por Daniela Anabel
Antamba Ávila, obteniendo un producto que cumple con todos los requisitos estipulados
por la Universidad Politécnica Salesiana para ser considerados como trabajo final de
titulación.

Quito D.M., Noviembre de 2018

………………………………….
Carlos Andrés Barrera Singaña
C.C.: 0503503336

II
 CESIÓN DE DERECHOS DE AUTOR

Yo, Daniela Anabel Antamba Ávila, con documento de identificación N° 1003686357,

manifiesto mi voluntad y cedo a la Universidad Politécnica Salesiana la titularidad sobre
los derechos patrimoniales en virtud de que soy autora del trabajo de grado/titulación
intitulado: “ESTABILIDAD ANGULAR DE UN SISTEMA ELÉCTRICO DE POTENCIA
CONSIDERANDO LA INTEGRACIÓN DE REGULADORES DE VELOCIDAD Y
VOLTAJE”, mismo que ha sido desarrollado para optar por el título de: Ingeniero
Eléctrico, en la Universidad Politécnica Salesiana, quedando la Universidad facultada
para ejercer plenamente los derechos cedidos anteriormente. En aplicación a lo
determinado en la Ley de Propiedad Intelectual, en mi condición de autor me reservo los
derechos morales de la obra antes citada. En concordancia, suscribo este documento en
el momento que hago entrega del trabajo final en formato impreso y digital a la
Biblioteca de la Universidad Politécnica Salesiana.

Quito, D.M., Noviembre de 2018

………………………………….
Daniela Anabel Antamba Ávila
C.C.: 1003686357

III
1. ÍNDICE GENERAL

1 Introducción ............................................................................................................... 2
2 Marco teórico ............................................................................................................. 4
2.1 Estabilidad de Sistemas Eléctricos de Potencia ...................................................... 4
2.2 Estabilidad Angular ................................................................................................. 4
2.2.1 Relación entre estabilidad angular y estabilidad transitoria ......................... 4
2.3 Modelos del sistema eléctrico de potencia para estudios de estabilidad angular .... 5
2.3.1 Generadores sincrónicos .............................................................................. 5
2.3.2 Transformadores .......................................................................................... 6
2.3.3 Líneas de transmisión ................................................................................... 7
2.4 Reguladores de velocidad ....................................................................................... 7
2.5 Reguladores de voltaje ............................................................................................ 7
3 Métodos de solución de la ecuación de oscilación..................................................... 8
3.1 Método Runge-Kutta 4to orden .............................................................................. 8
3.2 Modelación matemática .......................................................................................... 9
3.3 Modelación de los sistemas reguladores ............................................................... 10
4 Implementación del modelo matemático ................................................................. 11
4.1 Sistema propuesto para estudio de estabilidad ...................................................... 12
4.2 Modelo para estudio de estabilidad ....................................................................... 12
5 Análisis de resultados ............................................................................................... 13
6 Conclusiones ............................................................................................................ 16
6.1 Trabajos futuros .................................................................................................... 17
7 Referencias ............................................................................................................... 17
7.1 Matriz de Estado del Arte ..................................................................................... 20
7.2 Resumen de Indicadores ....................................................................................... 22
8 Anexos ..................................................................................................................... 24

IV
2. ÍNDICE DE FIGURAS

Figura 1. Estudio de estabilidad angular en un SEP. ........................................................ 3

Figura 2. Sistema estable de un grupo de generadores [9]. ............................................... 4
Figura 3. Sistema inestable de un grupo de generadores [9].............................................. 4
Figura 4. Circuito equivalente de un generador sincrónico – Modelo clásico [6]. ............ 5
Figura 5. Transformador de 2 devanados con tap [20]. ..................................................... 6
Figura 6. Modelo π del transformador [20]. ....................................................................... 7
Figura 7. Modelo π de líneas de transmisión. .................................................................... 7
Figura 8. Reguladores de velocidad y voltaje según IEEE [24]. ....................................... 8
Figura 9. Método de Runge-Kutta 4to Orden. ................................................................... 9
Figura 10. Modelo clásico de generador sincrónico para estudios de estabilidad. .......... 10
Figura 11. Modelación para un regulador de velocidad básico [24]. ............................... 11
Figura 12. Modelación para un regulador de voltaje básico [24]. ................................... 11
Figura 13. Sistema IEEE de 9 barras [13]. ....................................................................... 12
Figura 14. Ángulos de los rotores de los generadores. Sin Reguladores. ........................ 13
Figura 15. Diferencia angular. Sin Reguladores. ............................................................. 13
Figura 16. Ángulos de los rotores de los generadores. Reguladores de Voltaje. ............. 14
Figura 17. Diferencia angular. Reguladores de Voltaje. .................................................. 14
Figura 18. Voltaje interno generado en los generadores con reguladores de voltaje. ...... 14
Figura 19. Ángulos de los rotores de los generadores. Reguladores de Velocidad. ........ 15
Figura 20. Diferencia angular. Reguladores de Velocidad. ............................................. 15
Figura 21. Potencias Mecánicas de los generadores con reguladores de Velocidad. ...... 15
Figura 22. Ángulos de los rotores de los generadores. Reguladores de Velocidad y
Voltaje. ............................................................................................................................. 15
Figura 23. Diferencia angular. Reguladores de Velocidad y Voltaje. ............................. 16
Figura 24. Voltaje interno generado en los generadores. Reguladores de Voltaje y
Velocidad. ........................................................................................................................ 16
Figura 25. Potencias Mecánicas de los generadores. Reguladores de Velocidad y
Voltaje. ............................................................................................................................. 16
Figura 26. Resumen e Indicador de la Temática - Estado del Arte.................................. 22
Figura 27. Indicador de la formulación del problema – Estado del Arte ......................... 22
Figura 28. Indicador de solución – Estado del Arte ......................................................... 23

V
3. ÍNDICE DE TABLAS

Tabla 1. Matriz de Estado del Arte .................................................................................. 20

Tabla 2. Parámetros de las Líneas y Cargas del Sistema IEEE de 9 Barras. ................... 24
Tabla 3. Parámetros de los Generadores del Sistema IEEE de 9 Barras.......................... 24

VI
 ESTABILIDAD ANGULAR DE UN SISTEMA
ELÉCTRICO DE POTENCIA CONSIDERANDO LA
INTEGRACIÓN DE REGULADORES DE
VELOCIDAD Y VOLTAJE

Resumen Abstract
La estabilidad angular es de gran Angle stability has a highly importance in
importancia en los Sistemas de Potencia Power Systems as their study and
ya que su estudio y entendimiento comprehension allows to maintain in
permiten mantener en condiciones normal conditions of working
normales de operación “sincronismo” a "synchronism" to the machines of the
las máquinas del sistema luego de verse system after being affected by a
afectado por una perturbación. En este disturbance. In this paper is presented a
trabajo se presenta un desarrollo el cual development that allows to observe the
permite observar el comportamiento de behavior of a power system undergone a
un sistema de potencia sometido a una disturbance in transient period of time,
perturbación en el período transitorio, including the generators, speed and voltage
incluyendo en los generadores regulators. The numerical method that
reguladores de velocidad y voltaje. El permit to solve the oscillation equations it
método numérico que permite resolver is from a multi-machine system named
el sistema de ecuaciones de oscilación fourth order Runge-Kutta, getting an
de un sistema multimáquina es el de appraisal of the angular behavior from the
Runge-Kutta de cuarto orden, synchronous machines with about to the
obteniendo una estimación del time.
comportamiento angular de las
maquinas sincrónicas con respecto al
tiempo. Keywords: Angular stability, Voltage
regulator, Speed regulator, Runge-Kutta.
Palabras Clave: Estabilidad Angular,
Regulación de Voltaje, Regulación de
Velocidad, Runge-Kutta.

1
1 Introducción con respecto al sistema, por
consecuencia los rotores de las máquinas
La estabilidad es uno de los casos más comienzan a oscilar produciendo
importantes de estudio cuando se refiere cambios en la aceleración de los
a los sistemas eléctricos de potencia generadores [9].
(SEP), ya que permite estudiar el Para el estudio de la estabilidad se
comportamiento del mismo ante una considera un generador sincrónico, el
perturbación física, garantizando su cual está ligado a conceptos mecánicos,
operación en condiciones normales, matemáticos y físicos. En donde, si se
confiabilidad, calidad y seguridad [1]. pretende garantizar dicha estabilidad se
Los sistemas eléctricos de potencia debe relacionar a los generadores con las
son vulnerables debido a la naturaleza ecuaciones de oscilación, que se
del evento de la falla, el lapso de tiempo producen por una perturbación transitoria
que dure la perturbación, provocando en el SEP.
inestabilidad en el sistema [2][3]. Existen varios métodos de integración
Al referirse a máquinas sincrónicas para resolver la ecuación de oscilación
interconectadas entre si se debe tener como se puede encontrar en [10]. En
presente que pueden perder su donde al ser métodos complejos y de alta
estabilidad debido al: tipo de falla, precisión es indispensable el uso de
localización y tiempo de despeje. Si la paquetes computacionales.
duración de la perturbación del sistema Los métodos de integración
es alta, mayor es la desviación de la determinan el comportamiento del
frecuencia y ángulo del rotor, implicando ángulo con respecto al tiempo,
inestabilidad [4][5]. obteniendo las curvas de oscilación de
La estabilidad angular está los generadores.
relacionada con el comportamiento de El método de integración escogido
los ángulos de los rotores ante una para el análisis de estabilidad es el
perturbación sea de gran o pequeña método de Runge-Kutta de 4to orden ya
magnitud [6]. Durante la perturbación, que presenta una mejor respuesta y
el ángulo de rotor del generador mayor precisión al momento de obtener
comienza a oscilar, si la falla no es tan las curvas de los ángulos de los rotores
severa el sistema llega a estabilizarse [11][12].
hasta que el ángulo del rotor regrese a su El presente documento propone un
estado original o un nuevo estado método para analizar la estabilidad
estable, mientras que, si la falla es angular, mediante la integración de
demasiado grande el ángulo del rotor reguladores de velocidad y voltaje, con el
oscila sin amortiguación, lo que conlleva detalle que se puede presentar el
a la inestabilidad del sistema [7]. desempeño de cada uno de los
Partiendo del estudio de estabilidad se generadores del sistema y sus
puede determinar algunas acciones características de mayor importancia
correctivas, entre las cuales se tiene: el para el estudio que son: ángulo del rotor,
tiempo critico de liberación de falla, el voltaje interno generado y las potencias
diseño de esquemas para pérdida de mecánica y eléctrica, lo cual permite
carga y el diseño de esquemas de disparo apreciar gráficamente si el generador
transferido [8]. llego a perder el sincronismo con los
Cuando un generador pierde el demás generadores del sistema.
sincronismo gira a una velocidad distinta

2
 Figura 1. Estudio de estabilidad angular en un SEP.

En la presente investigación se realiza Modelado matemático del método de

el análisis de estudio, tomando como Runge-Kutta y los reguladores de
referencia el sistema de 9 barras velocidad y voltaje. Sección IV:
propuesto por la IEEE [13]. En la Implementación del modelo matemático.
Figura 1 se indica de manera gráfica las Sección V: Se indica el análisis de
distintas temáticas asociadas a un estudio resultados, en base a la metodología
de estabilidad angular, indicando que las propuesta. Sección VI: Se realiza las
fallas se producen en el sistema de conclusiones y trabajos futuros que se
transmisión. pueden llevar a cabo. Y finalmente en la
La estructura del presente documento Sección VII: Se muestra las referencias y
se encuentra dividido de la siguiente anexos.
manera: Sección I: Se muestra la
Introducción y antecedentes relacionados
al estudio de estabilidad angular.
Sección II: Se presenta el marco teórico
relacionado al concepto de estabilidad en
el SEP y estabilidad angular. Sección III:

3
2 Marco teórico
2.1 Estabilidad de Sistemas
Eléctricos de Potencia
El Sistema de Potencia se encuentra en
cambios constantes, por ello es necesario
garantizar niveles óptimos de voltaje,
potencias activas, potencias reactivas,
frecuencia de operación y ángulos de
rotores de las máquinas. En donde a
pesar de sufrir una perturbación, esta no
debe implicar ningún cambio [10]. Figura 2. Sistema estable de un grupo de generadores
En los SEP se considera tres factores [9].

que afectan la estabilidad del sistema. La

estabilidad de voltaje que se refiere en Mientras tanto en la Figura 3 se puede
conservar niveles óptimos de voltaje observar el mismo sistema de
garantizando la operación del SEP, la generadores, pero los ángulos proyectan
estabilidad de frecuencia que trata de una separación considerable, indicando
mantener un balance entre generación- que el sistema se encuentra en
carga y la estabilidad angular tiene como inestabilidad, perdiendo sincronismo
finalidad mantener en sincronismo los [9][15].
generadores [3].

2.2 Estabilidad Angular

Es la capacidad que los generadores
sincrónicos tienen para retomar
condiciones normales de operación
después de ser sujetos a una perturbación
leve o severa, donde la máquina
sincrónica tiene la facilidad de mantener
o restaurar el equilibrio entre la potencia
eléctrica y la potencia mecánica. Si este
equilibrio no se lograse el sistema entra
en inestabilidad produciendo un aumento Figura 3. Sistema inestable de un grupo de generadores
de oscilaciones angulares llevando a la [9].
perdida de sincronismo [6][14].
En la Figura 2 se puede observar un 2.2.1 Relación entre estabilidad
sistema en estado estable luego de una angular y estabilidad transitoria
perturbación. El comportamiento de los Cuando existen grandes perturbaciones
generadores es oscilatorio en donde los en el SEP, se utiliza la estabilidad
ángulos no proyectan una separación transitoria. En donde, el tiempo para el
demasiado grande. estudio de estabilidad es de 3 a 5
segundos después de ocurrida la falla y
puede extenderse hasta los 10 segundos
según la severidad y la probabilidad de
ocurrencia de la perturbación, en donde

4
por lo general, la respuesta del sistema se Existen dos tipos de rotor, los cuales
encuentra principalmente en los ángulos dependen de la velocidad. Para una baja
de potencia de los generadores del velocidad se utiliza el rotor de polos
sistema [3][10][16]. salientes, ya que manejan un mayor
número de polos como son las turbinas
2.3 Modelos del sistema eléctrico hidráulicas y para las que operan a altas
de potencia para estudios de velocidades se utiliza el rotor cilíndrico,
estabilidad angular que manejan un menor número de polos
2.3.1 Generadores sincrónicos como son las turbinas a gas y a vapor.
Constituyen una de las fuentes Considerando un generador en estado
principales de energía eléctrica para los estable se puede representar como una
SEP. Donde se convierte la energía fuente de voltaje y una impedancia en
mecánica en energía eléctrica serie.
abasteciendo grandes cargas, siendo a su
vez un punto muy importante en la
estabilidad.
Las máquinas sincrónicas poseen dos
elementos esenciales. Uno de ellos es el
devanado de campo, que transporta la
corriente continua y produce un campo
magnético y el otro es el devanado de
Figura 4. Circuito equivalente de un generador
armadura, el cual toma estos campos sincrónico – Modelo clásico [6].
magnéticos induciendo voltajes alternos,
en los devanados correctamente
distribuidos a 120° de separación de Donde:
manera que, se tiene una rotación
Voltaje interno generado.
uniforme del campo magnético.
Para la producción de un torque Caída de voltaje en la
estable, los campos del rotor y el estator resistencia de armadura.
deben girar a la misma velocidad
“velocidad sincrónica”. Donde el número Caída de voltaje en la
de polos de campo, se encuentra reactancia propia de
expresado por medio de la velocidad armadura.
mecánica del rotor y la frecuencia
eléctrica del estator [10]. Caída de voltaje en la
reactancia mutua de
armadura.
(1)
Voltaje de fase del
Donde: generador.

Velocidad en [rev/min]
En la Figura 4 se indica el circuito
Frecuencia en [Hz] equivalente del generador para
cualquiera de las fases. La ecuación de
Número de polos.

5
manera fasorial se indica a continuación El momento de inercia también puede
[6]: ser llamado momento de inercia del
rotor, y se puede represar con la letra M,
(2) y este a su vez en términos de H
(constante de inercia).
Donde es la impedancia sincrónica
y se encuentra definida por:
(7)
(3)
Finalmente, como resultado de
reemplazar las ecuaciones (5), (6), (7) en
Mientras que, para el generador en
la ecuación (4) se obtiene la ecuación de
estado dinámico, se representa con el
oscilación [17][18].
torque de aceleración el cual se obtiene
mediante el producto del momento de
inercia de la máquina por la aceleración (8)
angular.
2.3.2 Transformadores
(4) El transformador tiene como función
principal elevar o disminuir el voltaje
Cuando un cuerpo puesto en con la finalidad de garantizar
movimiento trata de oponerse se lo condiciones de carga y mantener
define como momento de inercia. operativo el sistema. Considerando un
transformador ideal se tiene un modelo
de dos devanados que están enrollados
(5) alrededor de un núcleo magnético [19],
en dicho modelo se puede considerar
En donde m es la masa del cuerpo y R como despreciable la rama de
es el radio del cuerpo. Por otro lado, el magnetización siempre y cuando se
comportamiento del ángulo de potencia δ trabaje en estado estable, por la razón de
en el estudio de estabilidad angular que la corriente de excitación es
puede ser expresado con respecto a un insignificante en comparación a la
eje de referencia que rote a velocidad corriente de carga [20].
sincrónica.
(6)

Donde:
Posición angular con
respecto al eje de referencia.

Velocidad sincrónica del Figura 5. Transformador de 2 devanados con tap [20].

generador. Donde:
Valor en por unidad del
Desplazamiento angular
cambiador de tap.
desde el eje de referencia.

6
 Resistencia y reactancia de
primario y secundarios (10)
referidos al lado de bajo
voltaje.
(11)
Al modelar el transformador mediante
el modelo π con cambiador de taps, se
tiene el siguiente circuito [21]:

Figura 7. Modelo π de líneas de transmisión.

2.4 Reguladores de velocidad

Para contralar la velocidad de los
generadores se utiliza un mecanismo el
cual es responsable de controlar las
Figura 6. Modelo π del transformador [20].
válvulas de las turbinas de vapor ó la
posición de las compuertas en las
Donde y es la admitancia del modelo turbinas hidráulicas, ajustándose el
y se expresa de la siguiente manera: torque mecánico. Esto ocurre en
condiciones normales de operación y
(9) durante perturbaciones [10][23].

2.5 Reguladores de voltaje

2.3.3 Líneas de transmisión
En el sistema de excitación del
Son elementos que transportan la energía
generador, el regulador de voltaje es la
eléctrica desde la generación a la carga
parte principal, ya que es capaz de
siendo un factor importante para la
detectar cambios en el voltaje de salida,
estabilidad del sistema y mantener
tomándose las medidas correspondientes
niveles óptimos de voltaje. Una línea de
para su restauración en caso de existir
transmisión se caracteriza por tres
una variación no deseada en el SEP.
parámetros:
Independientemente de la velocidad de
excitación de la respuesta, este no se verá
Resistencia en [Ω].
afectado hasta que el regulador de voltaje
se lo indique. En el caso que se tenga un
Reactancia serie en [Ω].
regulador lento, este influenciara a la
estabilidad del sistema, por ello es
Suceptancia paralelo en [S].
importante indicar los parámetros de
control proporcional para que el
Estos parámetros en función de la
regulador aumente los límites de
longitud de la línea de transmisión dan
estabilidad del generador y por
lugar a los distintos modelos de una línea
consecuencia las del sistema.
de transmisión, donde están presentes la
En la Figura 8, se presenta el
impedancia serie y la admitancia paralelo
diagrama de bloques propuesto por la
(Figura 7) [22].
IEEE, el cual contiene los dos
reguladores antes mencionados y como

7
estos se integran a un generador
sincrónico [24].
(12)
Figura 8. Reguladores de velocidad y voltaje según IEEE [24].

3 Métodos de solución de la ( ) (13)

ecuación de oscilación
3.1 Método Runge-Kutta 4to ( ) (14)
orden
El método Runge-Kutta, es un método
explicito, en donde el valor de la variable (15)
dependiente x en cualquier valor de t se
calcula a partir del conocimiento de los [ ] (16)
valores de x de los pasos de tiempo
anteriores. En donde x, es el ángulo de
los generadores y t es el tiempo de Donde:
duración de la perturbación. Este método Pendiente inicio del intervalo
es adecuado implementarlo en el caso de
estudio por la razón de que tiene una alta Primera pendiente en la mitad
precisión [10][11]. del intervalo usando k1.
Runge-Kutta se construye a partir del
método de Taylor, en donde Runge- Segunda pendiente en la
Kutta tiene la ventaja que elimina la mitad del intervalo usando k2.
necesidad de calcular y evaluar derivadas
de , haciendo uso de proyecciones Pendiente al final del
de punto medio [25][26]. La solución intervalo usando k3.
propuesta es equivalente a la aplicación
de las fórmulas de Taylor de cuarto La interpretación grafica de las
orden y se puede representar con base en pendientes estimadas para modelo de
la regla de 1/3 de Simpson de la Runge-Kutta, se lo puede apreciar en la
siguiente manera [12]: Figura 9.

8
 tiene los elementos fuera de la diagonal
[9][27].
̅ (23)

̅ (24)

[ ] (25)

Una vez obtenida la matriz de

admitancia se procede a calcular la
misma en condiciones previas a la falla,
Figura 9. Método de Runge-Kutta 4to Orden. durante la falla y luego que la falla haya
ocurrido y .
En cuanto a este documento, las
ecuaciones de cuarto orden de Con la obtención de la matriz de
Runge-Kutta son: admitancias, se puede determinar las
matrices reducidas (n x n) empleando el
(17) método de Kron (26) [28][29]:
 Se proporcionan barras
adicionales para los voltajes
( ) (18) internos del generador, agregando
la simulación de una impedancia
de falla de ser necesario.
( ) (19)
 Todos los elementos de
impedancias se deben modificar
(20) en admitancias.

[ ] (21) [ ] (26)

3.2 Modelación matemática [ ] (27)

Para un estudio de estabilidad angular es
necesario conocer como parámetros En donde n se usa para mencionar las
iniciales. los elementos que conforman el barras de los generadores mientras que r
sistema, obteniendo la matriz de se usa para los nodos restantes. Para
admitancias conocida como “Ybarra”. obtener el voltaje interno de los
generadores E, se debe partir de los datos
̅ ̅ ̅ (22) de flujo de potencia.

Donde ̅ contiene los elementos de la

matriz de admitancias, ̅ tiene los
elementos diagonales de la matriz y ̅

9
 Al igualar los componentes reales e
imaginarios, se puede obtener la
siguiente ecuación, dada por:

∑
(32)
Figura 10. Modelo clásico de generador sincrónico
para estudios de estabilidad. Donde:
Magnitud del voltaje interno
En términos de la potencia, se derivan
del generador i.
las ecuaciones (28) y (29):
̅ ̅ Magnitud del voltaje interno
(28)
del generador j.

(29) Módulo de la matriz reducida

ij.
Determinándose finalmente el voltaje
interno del generador [9]: Ángulo de la matriz reducida
ij.

( ) ( ) (30) Ángulo del voltaje interno

del generador i.

Donde: Ángulo del voltaje interno

Magnitud del voltaje interno del generador j.
del generador.
En la modelación es necesario incluir
Ángulo del voltaje interno del la ecuación de oscilación ya que es la
generador. responsable del movimiento del rotor y
relaciona el torque de inercia con los
Voltaje de la barra de torque mecánicos y eléctricos [28] [29].
generación. En donde, y son condiciones
iniciales del ángulo del voltaje interno
Reactancia transitoria del generado y la velocidad respectivamente.
generador.
3.3 Modelación de los sistemas
Potencia activa de generación.
reguladores
Potencia reactiva de Para una máquina sincrónica se
generación. considera dos casos de regulación: un
regulador de voltaje y un regulador de
La potencia en la red en una barra i, se velocidad. El regulador de velocidad al
establece como la salida de potencia sufrir cambios en produce cambios en
eléctrica de la máquina i, definido como: el torque mecánico, el cual depende de la
caída de velocidad además de la función
̅ ̅ (31) de transferencia del regulador de
velocidad y la fuente de energía primaria

10
(turbina), como se muestra en la Figura
11. [ ] (34)

Donde:

Voltaje de referencia.

Voltaje inicial.
Figura 11. Modelación para un regulador de velocidad
básico [24]. Variación de voltaje.

[ ] (33) Voltaje interno generado.

Ganancia.
Donde:
Constante de tiempo.
Velocidad angular inicial.
4 Implementación del modelo
Velocidad angular final.
matemático
Desplazamiento angular. Los sistemas eléctricos de potencia están
sometidos comúnmente a perturbaciones,
Potencia mecánica de pequeña o gran magnitud que puede
causar la inestabilidad angular. El no
Ganancia. poder reintegrar el sistema al estado
inicial de operación conlleva a
Constante de tiempo consecuencias tales como: el aumento de
la oscilación angular de varios
Los reguladores de voltaje tienen la generadores perdiendo el sincronismo
facilidad de detectar cambios en el con respecto a los otros generadores del
voltaje de salida. El diagrama de bloques sistema, disminución en la transferencia
presenta como modelo una relación en el de potencia e incluso un apagón total,
dominio , entre el cambio del voltaje de conocido como Blackout.
excitación y voltaje entre los terminales Se considera un estudio de estabilidad
del generador, como se muestra en la angular mediante la ecuación de
Figura 12. oscilación y el método de Runge-Kutta,
el cual permite conocer el
comportamiento de los ángulos y las
potencias, cuando se agrega reguladores
de velocidad y voltaje a las maquinas
sincrónicas.
Para completar el estudio de
estabilidad angular son integradas las
gráficas del comportamiento de las
Figura 12. Modelación para un regulador de voltaje oscilaciones de ángulo de las maquinas
básico [24].
con respecto a valores iniciales,

11
comparando así los resultados obtenidos. Algoritmo Estabilidad Angular
Se considera el software MATLAB para Paso 1:
Ingreso de Datos del sistema.
la realización del modelo de estabilidad Paso 2:
angular. Construcción de la Yb antes, durante y después de la
perturbación.
4.1 Sistema propuesto para Paso 3:
Matrices Ynn, Ynr, Yrn, Yrr a partir del Paso 2.
estudio de estabilidad Paso 4:
Para realizar el análisis referente a la Reducción de la Yb obtenidas en el Paso 3.
Y Reducida = Ynn – (Ynr * Yrr -1 * Yrn)
estabilidad angular se toma como Y Reducida: Matriz de admitancias reducidas.
referencia el sistema de 9 barras Paso 5:
propuesto por la IEEE indicado en la E
Figura 13. E: Voltaje interno generado en cada generador
Paso 6:
Condiciones iniciales para:
Wsinc: Generadores
V, θ: Regulación de voltaje
ΔPm, ΔE: Regulación de velocidad.
Paso 7: Runge-Kutta
Constantes Km. Ke, Tm y Te
Tiempo de análisis = Time/dt
t < tprefalla, Yb Reducida Prefalla
tpre < t < tdespeje, Yb Reducida Falla
t > tdespeje, Yb Reducida Posfalla
Paso 8: Potencia Eléctrica
Pe(i): Pot. Eléctrica para n generadores
Paso 9: Proyección de Punto Medio Para t = 0
Inicio del intervalo.
w1(i), d1(i), ΔPm1(i), ΔPe1(i), ΔE1(i)
Proyección de Punto Medio Para t = dt/2
Segunda Proyección
w2(i), d2(i), ΔPm2(i), ΔPe2(i), ΔE2(i)
Proyección de Punto Medio Para t = dt
Tercera Proyección.
w3(i), d3(i), ΔPm3(i), ΔPe3(i), ΔE3(i)
Proyección de Punto Medio Para t = 0.5
Final del intervalo.
w4(i), d4(i), ΔPm4(i), ΔPe4(i), ΔE4(i)
Paso 10: Formulación de las Ec’s Diferenciales
wr(i) → Vel. Angular, dr(i) → Angulo del Gen,
ΔPmr(i) → Pot. Mecánica, ΔPer(i) → Pot. Eléctrica,
Er(i) → Voltaje interno generado.
Paso 11:
Solución de las Ec’s Diferenciales obtenidas en el Paso
10.
Paso 12: Regulador de Velocidad
Si Reg-Vel = True {Pm(i) = Pm(i) + Pmr(i)}
Figura 13. Sistema IEEE de 9 barras [13]. Paso 13: Regulador de Voltaje
Si Reg-Vol = True {E(i) = E(i) + Er(i)}
4.2 Modelo para estudio de Paso 14: Valores Finales Aplicando Reguladores
di_gen: Angulo de los generadores
estabilidad wi_gen: Velocidad angular
A continuación, se muestra el algoritmo Ei_gen: Voltaje interno generado
Pm: Potencia Mecánica
implementado para resolver la ecuación Pe: Potencia Eléctrica
de oscilación mediante el modelo de Paso 15: Graficas
Runge-Kutta en un SEP multimáquina. Figura 1: [Tiempo, Angulo Generador]
Figura 2: [Tiempo, Angulo Final – Angulo Inicial]

12
5 Análisis de resultados
En el análisis de resultados se propone
determinar a partir de los ángulos del
rotor de los generadores, la estabilidad
transitoria con cuatro casos propuestos:
sin utilizar reguladores, utilizando
reguladores de voltaje, utilizando
reguladores de velocidad y utilizando
ambos reguladores.
En el planteamiento del problema se
tiene una falla trifásica en la línea 5-7
cercana a la barra 7. Dicha falla tiene una Figura 14. Ángulos de los rotores de los generadores.
duración de 5 ciclos (0.083 segundos). Sin Reguladores.
En donde los generadores se encuentran En donde se puede apreciar que todos
representados por el modelo clásico y las los generadores comienzan a oscilar y el
cargas por impedancias constantes. incremento angular inicia con los
Como resultado de cada caso se generadores 2 y 3 produciendo una
obtiene una gráfica de los ángulos de desaceleración. También se puede mirar
cada generador y otra gráfica con la que el grupo de los generadores alcanza
diferencia de los ángulos, donde se toma un ángulo máximo de 279.2° para un
como referencia al generador 1. tiempo de 2 segundos.
La diferencia de los ángulos permite En la Figura 15., se presenta la
determinar si el sistema es estable o diferencia de ángulos y se puede
inestable. Si la diferencia angular alcanza observar un ángulo máximo alcanzado de
valores máximos y luego decrecen el 80.68° en un tiempo de 0.43 segundos, lo
sistema es estable, mientras que si la que indica que el generador 2 alcanzo
diferencia angular se incrementa mayor oscilación.
indefinidamente el sistema es inestable.
Además de las gráficas mencionadas
anteriormente, se tiene las gráficas del
comportamiento para el voltaje interno
generado y las potencias mecánicas de
cada generador, los cuales determinan la
correcta implementación de los
reguladores con sus respectivas
ganancias y constantes de tiempo.
Los resultados obtenidos del estudio
de estabilidad angular sin aplicar
reguladores se indican en la Figura 14.

Figura 15. Diferencia angular. Sin Reguladores.

Para los resultados obtenidos del

estudio de estabilidad angular empleando
solo reguladores de voltaje, se incorpora

13
el modelo descrito en la ecuación (34),
con los siguientes valores:

Ganancia
Constante de tiempo
La Figura 16 se aprecia el
comportamiento de los ángulos de los
rotores. En donde se puede ver que todos
los generadores comienzan a oscilar y el
incremento angular inicia con los
generadores 2 y 3 produciendo una Figura 17. Diferencia angular. Reguladores de Voltaje.
desaceleración. También se puede ver
que el grupo de los generadores alcanzo En la Figura 18 se puede observar
un ángulo máximo de 277.6° para un que, al utilizar los reguladores de voltaje,
tiempo de 2 segundos. el voltaje interno generado de cada
generador sufre una variación como es
de esperarse luego de la falla. En donde,
estos valores son muy pequeños que no
afectan considerablemente al sistema.

Figura 16. Ángulos de los rotores de los generadores.

Reguladores de Voltaje.

En la Figura 17 de la diferencia de
ángulos, se puede observar un ángulo
Figura 18. Voltaje interno generado en los generadores
máximo alcanzado de 80.61° en un con reguladores de voltaje.
tiempo de 0.43 segundos, lo que indica
que el generador 2 alcanzo mayor Para los resultados obtenidos del
oscilación, siendo de un valor menor que estudio de estabilidad angular utilizando
el logrado sin reguladores. solo reguladores de velocidad, se
incorpora el modelo descrito en la
ecuación (33) con los siguientes valores:

Ganancia
Constante de tiempo
La Figura 19 muestra la actuación de
los ángulos de los rotores. En donde se
puede visualizar que todos los

14
generadores comienzan una oscilación,
pero con una escala muy pequeña en
comparación a la de los casos anteriores.

Figura 21. Potencias Mecánicas de los generadores con

reguladores de Velocidad.

En el último caso se utilizó las

mismas ganancias y constantes de tiempo
Figura 19. Ángulos de los rotores de los generadores.
indicados en los casos de reguladores de
Reguladores de Velocidad. velocidad y voltaje, en donde se obtiene
el siguiente comportamiento de los
En la Figura 20 existe la diferencia de ángulos de los rotores ilustrados en la
ángulos y se puede observar un ángulo Figura 22.
máximo alcanzado de 19.53° en un
tiempo de 0.07 segundos, lo que indica
que el generador 2 alcanzo mayor
oscilación.

Figura 22. Ángulos de los rotores de los generadores.

Reguladores de Velocidad y Voltaje.

En donde se puede ver que todos los

Figura 20. Diferencia angular. Reguladores de generadores comienzan oscilaciones
Velocidad. angulares con los mismos valores del
Para las potencias mecánicas de los caso de estudio de solo reguladores de
generadores, se observa en la Figura 21 velocidad. Sin embargo, los valores de
que los generadores 2 y 3 comienzan a los voltajes internos generadores son
oscilar. Sin embargo, estos valores se diferentes ya que también se incluyen los
amortiguan en un tiempo de 1.2 reguladores de voltaje.
segundos gracias a los reguladores de En la Figura 23, la diferencia de
velocidad implementados. ángulos se puede observar un ángulo
máximo alcanzado de 19.53° en un
tiempo de 0.07 segundos, lo que indica

15
que el generador 2 alcanzo mayor
oscilación.

Figura 25. Potencias Mecánicas de los generadores.

Reguladores de Velocidad y Voltaje.

Al tener todos los casos de estudio, se

Figura 23. Diferencia angular. Reguladores de puede observar que los generadores que
Velocidad y Voltaje.
sufrieron mayor oscilación fueron los
Como es de esperarse en la Figura 24 que se encontraron más cercanos a la
los voltajes internos generados no falla y no perdieron sincronismo ya que
presentan mayor variación en ninguna diferencia de ángulos excedió el
comparación al caso de estudio sin valor de ángulo máximo de 90°.
utilizar reguladores ya que para este caso Por otro lado, al aplicar los
se incluyen los reguladores de voltaje reguladores se puede apreciar una
que llegan a mantener en niveles reducción considerable en las
adecuados los valores de E. oscilaciones de los ángulos de los
rotores, mejorando la estabilidad del
sistema eléctrico de potencia.

6 Conclusiones
Para la solución de la ecuación de
oscilación se puede utilizar diversos
métodos numéricos. Sin embargo, el
método de Runge-Kutta permite obtener
un mejor análisis de las curvas de
oscilación que se plantean en el
problema de estabilidad con una mejor
precisión.
Figura 24. Voltaje interno generado en los
generadores. Reguladores de Voltaje y Velocidad. Los reguladores tanto de velocidad
como de voltaje, al tener correctamente
Para las potencias mecánicas de los parametrizados los valores de ganancias
generadores, se observa en la Figura 25 y constantes de tiempo, permiten
que los generadores 2 y 3 comienzan a garantizar la estabilidad del sistema ante
oscilar. Donde los valores se amortiguan una perturbación.
en un tiempo de 1.2 segundos, siendo los El modelo implementado en el
mismos que el caso de estudio de solo presente documento para estabilidad
reguladores de velocidad. angular puede ser realizado para
cualquier SEP, donde los generadores

16
deben ser representados por el modelo Structural Vulnerability
clásico y las cargas por impedancias Assessment Based on an
constantes. Improved Maximum Flow
La estabilidad y sus fenómenos Approach,” IEEE Trans. Smart
asociados cuando son estudiados son de Grid, vol. 9, no. 2, pp. 777–785,
gran importancia ya que permiten tener Mar. 2018.
un mayor grado de confiabilidad en la [3] P. Kundur ; J. Paserba ; V.
operación del SEP. Ajjarapu ; G. Andersson ; A.
Los generadores sincrónicos al estar Bose ; C. Cañizares ; N.
sometidos a constantes cambios físicos Hatziargyriou ; D. Hill ; A.
como mecánicos pueden perder su Stankovic ; C. Taylor ; T. Van
estabilidad por ello es necesario Cutsem ; V. Vittal, “Definition
incorporar los reguladores de voltaje y and Classification of Power
velocidad ya que con ellos los valores de System Stability IEEE/CIGRE
voltaje y potencia no se ven ampliamente Joint Task Force on Stability
afectados. Terms and Definitions,” IEEE
Trans. Power Syst., vol. 19, no. 3,
6.1 Trabajos futuros pp. 1387–1401, Aug. 2004.
Los resultados obtenidos a partir de este [4] M.-E. Hamedani-golshan, A.
documento pueden ser tomados para Davoudi, and A. Bidram,
desarrollar futuras investigaciones como: “Loading constraints for first
 Incorporar modelos de swing stability margin
reguladores en los cuales sean enhancement of distributed
incluidos turbinas y excitatrices generation,” IET Gener. Transm.
en los reguladores de velocidad y Distrib., vol. 6, no. 12, pp. 1292–
voltaje respectivamente. 1300, Dec. 2012.
 Incrementar el PSS (Power [5] Xiaoyu Wang and W. Freitas,
System Stabilizer) en el regulador “Influence of Voltage Positive
de voltaje. Feedback Anti-Islanding Scheme
 Optimización de los parámetros on Inverter-Based Distributed
de los reguladores ganancias y Generator Stability,” IEEE Trans.
constantes de tiempo, definiendo Power Deliv., vol. 24, no. 2, pp.
una función objetivo adecuada. 972–973, Apr. 2009.
[6] C. Barrera and J. Játiva, Análisis
de Estabilidad Angular en un
7 Referencias Sistema de Laboratorio, vol. 33,
[1] N. Zhou, D. Trudnowski, J. W. no. 1. [Verlag nicht ermittelbar],
Pierre, S. Sarawgi, and N. Bhatt, 2014.
“An algorithm for removing [7] M. F. Trebilcock, F. Santamaría,
trends from power-system and J. A. Alarcón, “ANÁLISIS
oscillation data,” in 2008 IEEE DE ESTABILIDAD
Power and Energy Society TRANSITORIA EN UN
General Meeting - Conversion SISTEMA INDUSTRIAL CON
and Delivery of Electrical Energy GENERACIÓN PROPIA
in the 21st Century, 2008, pp. 1–7. INTERCONECTADO CON EL
[2] J. Fang, C. Su, Z. Chen, H. Sun, SISTEMA DE POTENCIA,” Inf.
and P. Lund, “Power System tecnológica, vol. 25, no. 1, pp. 77–

17
 84, 2014. Dekker, 2002.
[8] W. S. Zimmermann, S. Hopp, M. [17] U. G. Vaidya, R. N. Banavar, and
Bondeur, and Der-Nien Chen, N. M. Singh, “A chaotic
“Transient stability study of the phenomenon in the damped power
Hsin Yu Co-Generation plant in swing equation,” in Proceedings
Hsin-Chu science based industrial of the 38th IEEE Conference on
park in Taiwan,” in 2000 IEEE Decision and Control (Cat.
Power Engineering Society Winter No.99CH36304), vol. 5, pp. 4650–
Meeting. Conference Proceedings 4655.
(Cat. No.00CH37077), vol. 1, pp. [18] Jun Zhou and Y. Ohsawa,
452–457. “Improved Swing Equation and Its
[9] P. Anderson and A. Fouad, Power Properties in Synchronous
System Control And Stability. . Generators,” IEEE Trans. Circuits
[10] P. Kundur, N. Balu, and M. Syst. I Regul. Pap., vol. 56, no. 1,
Lauby, Power System Stability pp. 200–209, Jan. 2009.
and Control. 1994. [19] D. J. Glover and M. . Sarma,
[11] E. Hernández, “Metodología Para Power System Analysis & Design,
El Análisis En Estado SI. 2012.
Estacionario De Sistemas [20] C. Javier Jiménez Salazar,
Eléctricos No Lineales,” IPN, “Incorporación de Modelos de
Mexico, 2004. Carga en Estado Estable Para el
[12] E. Cruz and R. Lizama, “Software Cálculo de Flujos de Potencia,”
Para Integración de Funciones y Escuela Politécnica Nacional,
Solución de Ecuaciones 2018.
Diferenciales Por Métodos [21] O. Manoochehri, A. Asoodeh, and
Numéricos,” PONTIFICIA K. Forooraghi, “π Model Dual-
UNIVERSIDAD CATÓLICA Band Impedance Transformer for
DEL PERÚ ESCUELA DE Unequal Complex Impedance
GRADUADOS, 2006. Loads,” IEEE Microw. Wirel.
[13] S. A. Shaikh, K. Kumar, A. R. Components Lett., vol. 25, no. 4,
Solangi, S. Kumar, and A. A. pp. 238–240, Apr. 2015.
Soomro, “Short Circuit Analysis [22] L. H-Restrepo Gladys Caicedo
&amp; Over current Relaying Delgado Ferley Castro-Aranda,
Coordination of IEEE 9-Bus “Transmission line models for
System,” in 2018 5th International electromagnetic transients in
Multi-Topic ICT Conference power systems,” 2008.
(IMTIC), 2018, pp. 1–6. [23] IEEE Power Engineering Society.,
[14] C. Gallardo, D. Andagoya, G. Power Development and Power
Salazar, and H. Arcos, Análisis de Generation Committee., Institute
Estabilidad Angular del Sistema of Electrical and Electronics
Eléctrico Ecuatoriano con Alta Engineers., and IEEE-SA
Penetración Eólica, vol. 33, no. 1. Standards Board., IEEE
[Verlag nicht ermittelbar], 2014. recommended practice for
[15] E. W. Kimbark, Power system preparation of equipment
stability. IEEE Press, 1995. specifications for speed-governing
[16] R. Natarajan, Computer - Aided of hydraulic turbines intended to
Power System Analysis. Marcel drive electric generators. Institute

18
 of Electrical and Electronics
Engineers, 2007.
[24] IEEE, 1207-2011 IEEE Guide for
the Application of Turbine
Governing Systems for
Hydroelectric Generating Units.
2011.
[25] R. L. Burden, J. D. Faires, and A.
M. Burden, Numerical analysis. .
[26] M. Syai’in and Kuo Lung Lian,
“Microgrid power flow using
Homotopic and Runge-Kutta
Method,” in 2015 IEEE 2nd
International Future Energy
Electronics Conference (IFEEC),
2015, pp. 1–8.
[27] C. A. Vásquez Barba, “Óptima
ubicación de fuentes de
generación distribuida basado en
el perfil de tensión de sistemas,”
2016.
[28] A. M. Kettner and M. Paolone,
“On the Properties of the Power
Systems Nodal Admittance
Matrix,” IEEE Trans. Power Syst.,
vol. 33, no. 1, pp. 1130–1131, Jan.
2018.
[29] F. Dorfler and F. Bullo, “Kron
Reduction of Graphs With
Applications to Electrical
Networks,” IEEE Trans. Circuits
Syst. I Regul. Pap., vol. 60, no. 1,
pp. 150–163, Jan. 2013.

19
7.1 Matriz de Estado del Arte
Tabla 1. Matriz de Estado del Arte

ESTABILIDAD ANGULAR DE UN SISTEMA ELÉCTRICO DE POTENCIA CONSIDERANDO LA INTEGRACIÓN DE REGULADORES DE

VELOCIDAD Y VOLTAJE

FORMULACIÓN DEL
PROPUESTAS PARA
PROBLEMA RESTRICCIONES DEL SOLUCIÓN
DATOS TEMÁTICA RESOLVER EL
FUNCIONES PROBLEMA PROPUESTA
PROBLEMA
OBJETIVO

Modelos de reguladores de velocidad y

Inestabilidad de sistemas de potencia

Constantes de tiempo de reguladores

Modelación de sistemas de potencia

Estabilidad sistemas multimáquina

Sistemas de control de velocidad y

Regulación de velocidad y voltaje

Integración de reguladores de
Modelación de reguladores de

Ganancias de reguladores
TÍTULO DEL ARTÍCULO

Regulador de Velocidad

Estabilidad Transitoria
Regulador de Voltaje
Estabilidad Angular

Estabilidad Angular

Métodos numéricos

velocidad y voltaje

velocidad y voltaje
Runge-Kutta

Runge-Kutta
voltaje

voltaje
CITAS
ITEM

AÑO

An algorithm for removing trends from power-system

1 2008 35                
oscillation data
 
2 2016
Power System Structural Vulnerability Assessment Based on
11               
an Improved Maximum Flow Approach
  
Loading constraints for first swing stability margin
3 2012 8               
enhancement of distributed generation   
Influence of Voltage Positive Feedback Anti-Islanding
4 2009 20                 
Scheme on Inverter-Based Distributed Generator Stability

Análisis de Estabilidad Transitoria en un Sistema Industrial
5 2014 con Generación Propia Interconectado con el Sistema de 7                  
Potencia

6 2000
Transient stability study of the Hsin Yu Co-Generation plant
7                 
in Hsin-Chu science based industrial park in Taiwan

7 2008 Power System Control and Stability 5202                 

8 1994 Power System Stability and Control. 21327                 

9 2010 Numerical analysis 292                 

10 2015
Microgrid power flow using Homotopic and Runge-Kutta
Method
14                  
Short Circuit Analysis &amp; Over current Relaying
11 2018
Coordination of IEEE 9-Bus System
0                  
12 2014
Análisis de Estabilidad Angular del Sistema Eléctrico
Ecuatoriano con Alta Penetración Eólica
0                  
13 2014 Análisis de Estabilidad Angular en un Sistema de Laboratorio 1                
 

20
14 1995 Power system stability 758                

15 2002 Computer - Aided Power System Analysis. 74   
              
16 2004 Definition and Classification of Power System Stability 2417   
              
17 2012 Power System Analysis & Design 1388   
             

18 2018
Incorporación de Modelos de Carga en Estado Estable Para el
Cálculo de Flujos de Potencia
0                  
π Model Dual-Band Impedance Transformer for Unequal
19 2015
Complex Impedance Loads
23                  
Transmission line models for electromagnetic transients in
20 2008
power systems
12                  
recommended practice for preparation of equipment
21 2007 specifications for speed-governing of hydraulic turbines 2                 
intended to drive electric generators 
Óptima ubicación de fuentes de generación distribuida basado
22 2016 0              
en el perfil de tensión de sistemas
   
23 2004
Metodología Para El Análisis En Estado Estacionario De
0              
Sistemas Eléctricos No Lineales
   
24 2006
Software Para Integración de Funciones y Solución de
0                 
Ecuaciones Diferenciales Por Métodos Numéricos

On the Properties of the Power Systems Nodal Admittance
25 2017 5                 
Matrix

26 2013
Kron Reduction of Graphs with Applications to Electrical
Networks
336                  
27 2009
Improved Swing Equation and Its Properties in Synchronous
Generators
23                  
28 1999 A chaotic phenomenon in the damped power swing equation 7                 

CANTIDAD: 14 2 1 5 9 7 6 1 10 12 4 4 3 4 1 4 2 1

21
7.2 Resumen de Indicadores

TEMÁTICA
Estabilidad Angular
15

10

Estabilidad Transitoria 5 Regulador de Voltaje

Runge-Kutta Regulador de Velocidad

Figura 26. Resumen e Indicador de la Temática - Estado del Arte

FORMULACIÓN DEL PROBLEMA

12

10

0
Modelación de sistemas Estabilidad Angular Sistemas de control de Estabilidad sistemas
de potencia velocidad y voltaje multimaquina

Figura 27. Indicador de la formulación del problema – Estado del Arte

22
Figura 28. Indicador de solución – Estado del Arte

23
8 Anexos
Tabla 2. Parámetros de las Líneas y Cargas del Sistema IEEE de 9 Barras.

Parámetros de la Líneas de Transformación de Carga a admitancia

Transmisión constante
DESDE HASTA R [p.u.] X [p.u] Y/2 [p.u.] G [p.u.] B [p.u.]
4 5 0.0100 0.0850 0.0880 0.000000000 0.000000000
4 6 0.0170 0.0920 0.0790 0.000000000 0.000000000
5 7 0.0320 0.1610 0.1530 1.260976262 -0.504390505
6 9 0.0390 0.1700 0.1790 0.877611554 -0.292537185
7 8 0.0085 0.0720 0.0745 0.000000000 0.000000000
8 9 0.0119 0.1008 0.1045 0.968959833 -0.339060000
1 4 0.0000 0.1184 0.0000 0.000000000 0.000000000
2 7 0.0000 0.1823 0.0000 0.000000000 0.000000000
3 9 0.0000 0.2399 0.0000 0.000000000 0.000000000

Tabla 3. Parámetros de los Generadores del Sistema IEEE de 9 Barras.

Datos de los Generadores

H [MW*s] Pm [p.u] Xd' [p.u.] θ [grad] V [p.u.] P [p.u.] Q [p.u.]
23.64 0.7163577 0.0608 0 1.040 0.7163577 0.27024337
6.4 1.629955043 0.1198 9.280018 1.025 1.62995504 0.06647461
3.01 0.849961462 0.1813 4.664742 1.025 0.84996146 -0.10872716

24