Unl

UNIVERSIDAD AUTNOMA DE NUEVO LEN
FACULTAD DE INGENIERA MECNICA Y ELCTRICA

SUBDIRECCIN DE POSGRADO
ANLISIS DE MODOS DE OSCILACIN Y AMORTIGUAMIENTO

EN SISTEMAS ELCTRICOS DE POTENCIA
POR
JUAN CARLOS RAMREZ SAUCEDO
COMO REQUISITO PARA OBTENER EL GRADO DE MAESTRA
EN CIENCIAS DE LA INGENIERA ELCTRICA CON
ESPECIALIDAD EN POTENCIA
NOVIEMBRE 2014
UNIVERSIDAD AUTNOMA DE NUEVO LEN

FACULTAD DE INGENIERA MECNICA Y ELCTRICA
SUBDIRECCIN DE POSGRADO

POR
JUAN CARLOS RAMREZ SAUCEDO
COMO REQUISITO PARA OBTENER EL GRADO DE MAESTRA
EN CIENCIAS DE LA INGENIERA ELCTRICA CON
ESPECIALIDAD EN POTENCIA
NOVIEMBRE 2014
Dedicatoria
Dedico este trabajo a mis padres, Gregoria y Antonio por su sencillez y apoyo
incondicional de siempre.
A mi pequea Anna, fuente de alegra inagotable y a su madre Sonia, mi fiable
y creativa, compaera del camino.
A mis hermanas Juanita, Sol, Ofelia, Norma, Lul, Mayra y Ara y a mis
hermanos, Gil, Too e Isma.
iv
Agradecimientos
Mi agradecimiento al Dr. Florencio Aboytes Garca por su valioso apoyo y
asesora en el desarrollo de este trabajo y por compartir su conocimiento y
experiencia.
Al Dr. Jos Antonio de la O Serna por su asesora en temas especficos.
A los Doctores Manuel Antonio Andrade Soto, Ernesto Vzquez Martnez y al
M.C. Luis David Archiga Rodrguez, por sus valiosa asesora y comentarios
para enriquecer ste trabajo de tesis.
A quienes, dentro del programa doctoral, contribuyeron en mi formacin con
cordialidad y empeo, Dr. Salvador Acha Daza, Dr. Hctor J. Altuve Ferrer, Dr.
Joaqun Collado, Dr. Rogelio Ramrez Barradas, Dr. Fernando Snchez Tello,
Dr. Oscar L. Chacn Mondragn.
A mis compaeros del programa doctoral, compaerismo que se convirti en
amistad, Laura, Mara de los ngeles, Irma Rosario, David Alejandro, Julio
Alberto, Armando, Josu, Pedro Martn, Francisco Humberto, Gustavo y Mario.
Al CONACYT.
Al personal administrativo de la FIME, especialmente a la Sra. Lupita y al Ing.
Benito Garza.
A todas y todos, gracias.
Resumen
Publicacin No._______________
Juan Carlos Ramrez Saucedo
Universidad Autnoma de Nuevo Len
Facultad de Ingeniera Mecnica y Elctrica
Profesor Asesor: Dr. Manuel Antonio Andrade Soto
Noviembre de 2014
En este trabajo de tesis se analiza el origen y causas de las oscilaciones

electromecnicas que se presentan en Sistemas Elctricos de Potencia ante la
ocurrencia de pequeas perturbaciones. A partir de dicho anlisis, se obtiene
informacin sobre las frecuencias de oscilacin, que pueden presentarse en el
SEP, as como las magnitudes relativas y las formas caractersticas de
oscilacin de las variables del SEP, que participan en las oscilaciones.
Se analizan dos tcnicas, el anlisis modal y el anlisis de Fourier, mediante las
cuales se obtiene informacin valiosa que contribuye al conocimiento
anticipado, del comportamiento del SEP ante disturbios. Se describen los
fundamentos matemticos de la tcnica de valores caractersticos, mediante la
cual se obtienen los modos de oscilacin de diferentes SEPs de prueba con y
sin amortiguamiento. Se analiza e interpreta la relacin de los valores y
vi
vectores caractersticos y los factores de participacin, con las frecuencias y el

amortiguamiento de la oscilacin de las variables de estado del SEP tales como
la posicin y velocidad angular de rotores de mquinas y flujos de potencia
activa. Se comprueban estas relaciones mediante la excitacin de modos de
oscilacin o simulacin de fallas en programas de estabilidad transitoria.
Adems, se analiza el impacto de pequeas perturbaciones en otras variables
como el voltaje y la potencia reactiva encontrando diferencias significativas.
La tcnica del anlisis de Fourier, junto con el procesamiento de seales en
tiempo discreto, se emplean para obtener las componentes de frecuencia de
seales grabadas (respuesta no lineal del SEP), de disturbios ocurridos en
SEPs de prueba. Con estos resultados se obtienen tambin, los modos de
oscilacin naturales del sistema, aspectos de fase de la oscilacin de variables
y datos sobre, que modo de oscilacin tiene mayor influencia en el
comportamiento oscilatorio de alguna mquina en particular. Se presentan una
comparacin de los resultados obtenidos con ambas metodologas.
Finalmente, y como parte de la conclusin de este trabajo, se remarca la
utilidad de ambas tcnicas a ser usadas de forma complementaria, en el
anlisis de SEPs, contribuyendo a un mejor conocimiento del comportamiento
del SEP ante disturbios.
vii
ndice
Pgina
Dedicatoria. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iv
Agradecimientos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Resumen. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
vi
Captulo 1. Introduccin
1.1
Motivacin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.2
Revisin bibliogrfica
1.2
Estructura de la tesis. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.3
Aportaciones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Capitulo 2. Caracterizacin de modos de oscilacin

2.1
Introduccin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2
Sistemas lineales y ecuaciones de estado. . . . . . . . . . . . . . . . . . .
2.3
Anlisis modal . . . . . . . . . . . . . . . . . . . . . . . . . .
10
2.3.1
Modos de oscilacin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10
2.3.2
Eigenvalores o valores caractersticos. . . . . . . . . . . . . . . .
11
2.3.3
Eigenvectores o vectores caractersticos. . . . . . . . . . . . . .
12
2.4 La respuesta en el tiempo de un sistema lineal. . . . . . . . . . . . . .

2.4.1
14
Matriz de participacin. . . . . . . . . . . . . . . . . . . . . . . . . . . . .
16
2.5 Conclusiones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
16
Capitulo 3. Interpretacin de modos de oscilacin en sistemas

elctricos de potencia
3.1 Introduccin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
18
3.2 Frecuencias naturales no amortiguadas. . . . . . . . . . . . . . . . . . . . .
19
viii
ndice (continuacin)
3.2.1
Obtencin del modelo linealizado de las ecuaciones de

flujos de potencia. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2.2
19
El clculo de modos como un problema de valores

caractersticos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
23
3.3 Interpretacin de modos de oscilacin en sistemas elctricos de

potencia de prueba. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
25
3.4 Conclusiones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
33
Capitulo 4. Amortiguamiento de oscilaciones

4.1 Introduccin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
35
4.2 Coeficientes de sincronizacin y amortiguamiento en SEPs. . . .
36
4.3 Anlisis de coeficientes de sincronizacin y amortiguamiento en

un SEP bsico. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
37
4.4 Obtencin del modelo lineal de un SEP multimquina. . . . . . . . . .
42
4.5 Excitacin de modos de oscilacin. . . . . . . . . . . . . . . . . . . . . . . . .
45
4.6 Interpretacin de modos de oscilacin en SEPs de prueba y
46
anlisis de amortiguamiento. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.6.1
Interpretacin del amortiguamiento en los modos de

oscilacin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.7 Conclusiones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
46
53
Capitulo 5. Oscilaciones de voltaje, potencia activa y potencia

reactiva
5.1 Introduccin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
55
5.2 Ecuaciones de voltaje, potencia activa y potencia reactiva. . . . . .
56
5.2.1.
Ecuaciones de voltajes en terminales de generadores. . .
56
5.2.2.
Ecuaciones de potencia activa y potencia reactiva. . . . . .
58
5.3 Anlisis de oscilaciones de voltaje y potencia. . . . . . . . . . . . . . . .
60
ix
ndice (continuacin)
5.4 Ejemplos ilustrativos de oscilaciones de voltaje y potencia. . . . . .
62
5.5 Conclusiones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
72
Capitulo 6. Determinacin de modos de oscilacin mediante el

procesamiento de respuestas en el tiempo
6.1 Introduccin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6.2 Determinacin del espectro de frecuencias de seales de
respuesta en el tiempo en un sistema elctrico de potencia. . . . .
6.3 Aplicacin en sistemas elctricos de potencia de prueba. . . . . . .
6.3.1.
Determinacin de modos de oscilacin en un SEP de

prueba. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6.3.2.
74
75
76
76
Comparacin de resultados. . . . . . . . . . . . . . . . . . . . . . . .
81
6.4 Conclusiones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
85
Capitulo 7. Conclusiones y recomendaciones

7.1 Conclusiones finales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
87
7.2 Aportaciones. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
92
7.2 Recomendaciones para trabajos futuros. . . . . . . . . . . . . . . . . . . . .
93
Referencias. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
95
ndice de tablas. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
100
ndice de figuras. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
101
Apndice 1. Sistemas de prueba de casos de estudio. . . . . . . . . . . . . .
103
Apndice 2. Artculos presentados del trabajo de tesis. . . . . . . . . . . . .
108
CAPITULO 1
Introduccin
1.1. Motivacin
Los sistemas elctricos de potencia (SEPs) estn sujetos de manera
permanente a diversas perturbaciones debido principalmente a la variacin de
la carga, a operaciones normales, cambios en generacin, conexin y
desconexin de elementos, etc. En ocasiones estas perturbaciones involucran
variaciones de velocidad de rotores de generadores y estn asociadas con
oscilaciones, que generalmente decrecen con el tiempo, pero en algunas
circunstancias, pueden persistir despus de un disturbio incrementar de
magnitud, ocasionando el disparo de lneas que aslan parte del sistema con la
subsecuente prdida de carga.
En los ltimos aos se han detectado en muchos sistemas de potencia en el
mundo, la aparicin de oscilaciones espontneas de baja frecuencia, las cules
se mantienen durante periodos de tiempo relativamente grandes, debido a que
no existe amortiguamiento. En otros casos, ms graves an, las oscilaciones
han provocado la desconexin de lneas de interconexin por prdida de
sincronismo, debido al crecimiento de las oscilaciones en el tiempo.
Desde el punto de vista de anlisis es importante identificar el origen de las
oscilaciones, calcular las frecuencias naturales, entender la fsica del
amortiguamiento, evaluar el efecto de controles y proponer alternativas para el
control de oscilaciones. Desde el punto de vista prctico, en la operacin de
SEPs se requiere tener recomendaciones precisas para minimizar el problema
1
de oscilaciones y tener guas tcnicas para influir en el diseo de controles de

nuevas centrales generadoras. En este caso es importante el estudio y anlisis
del probable comportamiento del sistema, previo a la ocurrencia de las
oscilaciones ya que el tenerlas identificadas permitir una mejor comprensin
del problema as como una estrategia de prevencin al momento que se
presenten.
En esta tesis se describen y se aplican dos metodologas para el clculo de los
modos de oscilacin electromecnicos o frecuencias naturales del sistema
elctrico; la primera en base al conocimiento de los modelos de la red elctrica
con las caractersticas de sus elementos y la otra en base al anlisis de
registros de disturbios que se presentan en el SEP.
1.2. Revisin bibliogrfica

En diversos sistemas de potencia en el mundo se han presentado condiciones
de inestabilidad ante pequeas perturbaciones al presentarse oscilaciones
electromecnicas no amortiguadas [30,12]. Lo anterior ha motivado el desarrollo
de diferentes trabajos enfocados al anlisis y control de oscilaciones
electromecnicas. A este anlisis se le conoce tambin en la literatura como
anlisis de estabilidad de pequea seal. En [1], Byerly describe una base
terica para el anlisis de los modos normales y formas modales (anlisis
modal) de un sistema de potencia obteniendo datos para una mejor
comprensin del comportamiento dinmico. En [3,4,24] Prerz-Arriaga y
Verghese proponen un marco llamado Anlisis Modal Selectivo para analizar y
reducir los grandes modelos lineales invariantes en el tiempo que se originan en
el estudio de la Estabilidad Dinmica en sistemas de potencia. Mediante el, se
puede determinar de manera eficiente y con precisin, los valores propios y
vectores propios de los modos naturales de inters.
Aunque inicialmente se realizaron estudios de oscilaciones electromecnicas

mediante simulaciones dinmicas sobre el modelo no lineal del SEP, donde sus
elementos son modelados con la mayor cantidad de detalles posibles, se sabe
que el comportamiento de un SEP ante pequeas perturbaciones puede
compararse al comportamiento de sistemas lineales lo cual permite linealizar la
dinmica del SEP en un punto de operacin y la aplicacin de mtodos de
estudio de estabilidad de sistemas lineales. Los mtodos propuestos para al
anlisis de estabilidad de pequea seal son: a) anlisis modal, b) anlisis de
torque, c) mtodo de respuesta en frecuencia y d) mtodos de identificacin
siendo los mayor uso en el estudio de las oscilaciones, por su eficacia
comprobada en sistemas reales, el anlisis modal y los mtodos de
identificacin [30].
En [27] Hauer presenta resultados iniciales del anlisis Prony a seales de
respuesta de un sistema de potencia. El anlisis Prony realiza una extensin del
anlisis de Fourier, estimando directamente la frecuencia, el amortiguamiento,
la magnitud y la fase relativa de las componentes modales presentes en una
seal dada, obtenida generalmente como salida de un programa de estabilidad
transitoria. Este anlisis consiste en una estrategia matemtica que ajusta una
combinacin lineal de trminos exponenciales a la seal muestreada en el
tiempo.
Por sus caractersticas, el anlisis Prony, junto con los algoritmos SteiglitzMcBride y el Eigensystem Realization, son considerados mtodos de
identificacin. En [25], Sanchez-Gasca y Chow, realizan una comparacin
emprica de estos tres mtodos de identificacin utilizados en el anlisis de
oscilaciones electromecnicas.
La estabilidad de oscilaciones electromecnicas dbilmente amortiguadas entre
grupos de generadores es un requerimiento fundamental para la operacin
adecuada en un SEP. En [28] se utilizan el anlisis modal y un mtodo basado
3
en anlisis de residuos para la aplicacin de FACTS (Flexible Ac Transmision

Systems), especficamente se discute la aplicacin de Compensadores
Estticos de Vars (CEV) y Capacitores Serie Controlados por Tiristores
(CSCT),
para
amortiguar
oscilaciones
interrea
concluyendo
que
los
dispositivos FACTS mejoran el desempeo dinmico de un sistema de

potencia.
En publicaciones recientes se ha aplicado el anlisis modal para evaluar el
comportamiento dinmico del SEP ante pequeas perturbaciones permitiendo
la identificacin directa de las causas que contribuyen en forma negativa al
amortiguamiento asociado con los modos electromecnicos, as como la
ubicacin ms adecuada de controles que contribuyan a mejorar el
amortiguamiento de los mismos[29].
Un nuevo enfoque en la supervisin y control de oscilaciones electromecnicas
del sistema de potencia se plantea en [33], donde se ha identificado una
tendencia a emplear mediciones fasoriales sincronizadas, provistas por
unidades de medicin fasorial (PMU) en nuevos esquemas de supervisin y
control. Se propone una metodologa para la supervisin y control de
oscilaciones de baja frecuencia en tiempo real, basada en ndices de
controlabilidad y observabilidad. El mtodo para estimar la frecuencia y
amortiguamiento
de
las
oscilaciones
realiza
una
combinacin
de
la
transformada de wavelet discreta y el algoritmo Prony.
1.3. Estructura de la tesis

En el captulo 1 se presenta los motivos que dan origen a esta tesis, realiza una
revisin
bibliogrfica
con
los
antecedentes
relacionados
al
tema
de
investigacin, se presenta la estructura del contenido del trabajo y se incluyen

las aportaciones ms significativas. En el Captulo 2 se presenta la teora
fundamental para el anlisis de oscilaciones electromecnicas, se parte de la
4
tcnica de valores caractersticos y se interpreta el significado dichos valores y

sus vectores caractersticos asociados en relacin con las frecuencias de
oscilacin y formas caractersticas de las mismas.
En el capitulo 3 se aplica la teora del captulo 2 para el anlisis de oscilaciones
electromecnicas
no
amortiguadas
en
SEPs.
Se
presentan
algunas
consideraciones para obtener los modos de oscilacin naturales y formas

caractersticas de oscilacin del SEP mediante el clculo de valores y vectores
caractersticos as como las justificaciones para la linealizacin de las
ecuaciones de potencia en un punto de operacin. Se describe la interpretacin
de los resultados obtenidos mediante el clculo de los modos de oscilacin, en
variables como la posicin y velocidad angular de rotores y la potencia elctrica,
con algunos ejemplos ilustrativos en SEPs de prueba y se simulan algunos
modos con la ayuda de un programa de estabilidad transitoria.
El captulo 4 considera las ecuaciones generales que incluyen los pares de
sincronizacin y amortiguamiento en la dinmica de las maquinas generadoras.
Se analizan SEPs de prueba obteniendo los modos de oscilacin e
interpretando el efecto del amortiguamiento. Se realizan comparaciones
respecto al anlisis presentado en el captulo 3 para el caso no amortiguado.
As mismo se detalla el proceso matemtico de la excitacin de modos de
oscilacin en SEPs.
En el captulo 5 se analiza el impacto de pequeas perturbaciones en otras
variables que definen el comportamiento reactivo de los generadores y los
voltajes del sistema. Se desarrollan las ecuaciones para el voltaje en terminales
y la potencia reactiva generada, se obtienen expresiones en funcin de la
posicin angular de rotores para un sistema bsico, usando este desarrollo se
definen las regiones lineales y no-lineales del comportamiento de las unidades
en funcin del punto de operacin. Mediante el anlisis de casos ilustrativos se
muestra la deformacin de voltajes y potencias que resultan de la condicin
5
operativa y del movimiento de las posiciones angulares de rotores de unidades

generadoras. El anlisis se extiende a diversos modos de oscilacin y
confirmando el mismo efecto de no-linealidad.
El captulo 6 presenta una alternativa para el clculo de los modos naturales de

un sistema de potencia cuando no se conocen por completo los parmetros
fsicos de la red de manera que no pueden utilizarse las tcnicas del anlisis
modal. La alternativa propuesta requiere tener registros de datos del
comportamiento del sistema de potencia durante perturbaciones que se hayan
presentado en dicho sistema. Mediante el uso de tcnicas del anlisis de
seales y anlisis de Fourier se analiza la informacin contenida en los registros
para obtener informacin de los modos de oscilacin del sistema en cuestin.
Tambin en ste captulo se presentan ejemplos ilustrativos de la obtencin e
interpretacin de los modos naturales en SEPs de prueba. Se realizan
comparaciones con los resultados del anlisis modal presentado en captulos
anteriores.
Finalmente el captulo 7 presenta las conclusiones de la tesis y algunas
recomendaciones para trabajos futuros.
1.4. Aportaciones
Este trabajo presenta dos alternativas para el clculo de los modos de
oscilacin electromecnicos de un SEP. La primera de ellas consiste en una
aplicacin del anlisis de sistemas lineales al clculo de los modos de
oscilacin (anlisis modal) de un SEP demostrndose su validez para predecir
la respuesta natural ante pequeas perturbaciones. La segunda alternativa
presenta una metodologa para obtener los modos de oscilacin de un SEP a
partir de registros de disturbios.
Se describe en forma clara, la relacin entre los valores caractersticos,

eigenvalores, eigenvectores y factores de participacin con las frecuencias de
oscilacin electromecnicas del SEP, as como las magnitudes relativas y fases
de dichas oscilaciones.
Para oscilaciones amortiguadas se realiza la descripcin matemtica de la
respuesta del SEP en funcin de los eigenvalores y eigenvectores.
Se extiende el anlisis de pequeas perturbaciones al anlisis de otras
variables confirmando que los modos electromecnicos se reflejan en la
posicin y velocidad angular y generalmente tambin en la potencia activa y
flujos de potencia activa pero no necesariamente en el voltaje y potencia
reactiva.
Se obtienen los modos naturales de oscilacin en sistemas de potencia de
prueba a partir de registros de disturbios ocurridos en el sistema. En esta
metodologa es necesario considerar que una falla en particular que ocurra en
el SEP podra excitar solo algunos de los modos de oscilacin naturales y
necesariamente la respuesta no lineal estar compuesta por contribuciones de
dichos modos.
Ambas metodologas pueden combinarse y ser de gran utilidad al aplicarse en
el anlisis de sistemas de potencia, tanto para la identificacin de oscilaciones
que pudieran presentarse en el sistema como para influir en el diseo de
controles.
CAPITULO 2
Caracterizacin de Modos de
Oscilacin
2.1
Introduccin
Los sistemas elctricos de potencia (SEP) estn sujetos a diversas

perturbaciones, que pueden ser desde variaciones normales en la carga,
cambios en la potencia de generacin, conexin y desconexin de elementos,
hasta fallas en los sistemas de generacin y transmisin. Toda perturbacin
origina un desbalance entre los pares elctrico y mecnico de generadores y
turbinas respectivamente, que produce la oscilacin de los grupos turbinagenerador, que a su vez ocasiona oscilaciones en las diversas variables en el
sistema de potencia; tales como posiciones y velocidades angulares de rotores,
flujos de potencia, voltajes y corrientes, potencia generada entre otras.
En el SEP tambin existen diferentes fuentes de amortiguamiento que reducen
la amplitud y duracin de la oscilacin. Entre ellas los devanados
amortiguadores en generadores, el comportamiento de la potencia activa y
reactiva de las cargas con el voltaje, los controles en sistemas de excitacin y
gobernacin de velocidad. En el mejor de los casos, si la fuente de
amortiguamiento es efectiva las oscilaciones disminuyen su amplitud decrecen
con el tiempo. Sin embargo, en algunos casos, hay oscilaciones que persisten
en el tiempo e incluso oscilaciones que incrementan su magnitud, ocasionando
el disparo de lneas y la perdida de cargas en el sistema. Desde punto de vista
8
del anlisis dinmico de un sistema es de gran importancia identificar el origen

de oscilaciones para tratar de minimizar este problema.
En este captulo se presenta la teora fundamental para el anlisis de
oscilaciones electromecnicas, se parte de la tcnica de valores caractersticos
y se interpreta el significado dichos valores y sus vectores caractersticos
asociados en relacin con las frecuencias de oscilacin y formas caractersticas
de las mismas.
2.2
Sistemas lineales y ecuaciones de estado
Las ecuaciones dinmicas que rigen el comportamiento del sistema de potencia

son no-lineales. Sin embargo, para el anlisis de pequeas perturbaciones
pueden ser linealizadas alrededor de un punto de operacin. El modelo
linealizado puede ser expresado como un conjunto de ecuaciones diferenciales
de primer orden y a la representacin en variables de estado
x (t ) A(t ) x(t ) B(t ) (t )
y (t ) C (t ) x(t ) D(t ) (t )
(2-1)
donde x (t ) denota el vector de derivadas de las variables x(t) con respecto al

tiempo. El vector funcin del tiempo x(t) de orden (nx1) es el llamado vector de
estado y sus componentes x1(t), x2(t), ... xn(t), son las variables de estado del
sistema. La funcin u(t), la excitacin al sistema es un vector de orden (mx1),
y(t) es el vector de las seales de salida de orden (px1). A(t) (nxn) es la matriz
de estado, B(t) es la matriz (nxm) de coeficientes de las entradas, C(t) (pxn) es
la matriz de coeficientes de las salidas y D(t) (pxm) es la matriz de coeficientes
que relaciona directamente las variables de entrada y salida.
El sistema lineal ser invariante en el tiempo (LIT) si estas matrices de
coeficientes son constantes. Si una o ms de estas matrices varan con el
tiempo el sistema es variante en el tiempo (LVT).
9
En el anlisis de sistemas de potencia se considera un tiempo fijo t 0 inicial, y las

propiedades de la solucin x(t) de una ecuacin lineal de estado, para un
estado inicial dado x(t0)=x0 y una seal de entrada u(t) especificada, son de
inters para t t0.
Ecuaciones lineales de estado de la forma (2-1) pueden originarse de muchas
formas y en efecto, las ecuaciones diferenciales clsicas de n-simo orden
pueden llevarse a la forma de la ecuacin de estado.
2.3
Anlisis modal
2.3.1 Modos de oscilacin

Los modos de oscilacin naturales de un sistema elctrico describen la
respuesta libre, natural o no forzada del mismo. El objetivo es predecir el
funcionamiento futuro del sistema en el dominio del tiempo mediante el anlisis
en el dominio de la frecuencia compleja.
Para realizar el anlisis modal, se parte de la ecuacin de estado (2-1),
considerando que el sistema es LIT y se aplica la transformada de Laplace de
donde se obtiene una representacin alternativa para la ecuacin de estado en
el dominio de la frecuencia
X ( s) SI A x(0) SI A bU ( s)
1
donde X ( s) x(t )e st
Y ( s) c X ( s)
(2-2)
En general, la respuesta natural de un sistema LIT con un estado inicial x(t 0)=x0
en el tiempo t0=0 y en ausencia de nuevas entradas obedecer
X (s) SI A x0
1
(2-3)
10
La matriz [SI-A] es algunas veces llamada la matriz caracterstica de A, su

determinante
a(s) det SI A s - 1 s - 2 ...s - n
(2-4)
es conocido como el polinomio caracterstico de A, y sus races {i} son

llamados los valores caractersticos de A [13]. Los eigenvalores pueden ser
reales o complejos. Si A es real los eigenvalores necesariamente ocurren en
pares complejos conjugados. Para el sistema de potencia la matriz de estado A
tiene coeficientes reales.
Es claro, de la ecuacin (2-1), que estos eigenvalores describen la respuesta
natural (u=0) de la realizacin {A,b,c} (la descripcin en variables de estado), y
por lo tanto ellos son con frecuencia llamados los modos o frecuencias
naturales de oscilacin.
2.3.2 Eigenvalores o valores caractersticos
Los eigenvalores de A estn dados por los valores escalares de para los
cuales existe una solucin con pi diferente de cero
Ap p
donde:
(2-5)
A es una matriz nxn con coeficientes reales
p es un vector nx1
Para encontrar los eigenvalores la ecuacin (2-5) puede ser escrita como
A I p 0
(2-6)
para una solucin no trivial
11
det A I 0
(2-7)
Al encontrar las n soluciones de , obtenemos los eigenvalores de A.

2.3.3 Eigenvectores o vectores caractersticos
Para algn eigenvalor i, el vector columna pi el cual satisface la ecuacin (2-5)
es conocido como el eigenvector derecho de A asociado con el eigenvalor i
Api i pi
(2-8)
i 1,2,...n
el eigenvector pi tiene la forma

p1i
p
pi 2i

pni
(2-9)
En forma similar, el vector rengln de la forma

qiT qi1 qi 2 qin
(2-10)
que satisface la ecuacin

qiT A i qiT
(2-11)
i 1,2,...n
es llamado el eigenvector izquierdo de A asociado con el eigenvector i.

Para eigenvalores diferentes, los vectores caractersticos derecho e izquierdo
correspondientes, forman un conjunto ortogonal, es decir
12
qTj pi ki
i j
i j
(2-12)
En vista de que los eigenvectores estn determinados con algn multiplicador

escalar, es una prctica comn normalizar estos vectores de tal forma que para
i = j tenemos
(2-13)
ki qiT pi 1
Siguiendo con el anlisis de las propiedades de la matriz A, se definen las

siguientes matrices:
sea P pi [ p1 p2 ... pn ]
(2-14)
una matriz nxn con los vectores caractersticos derechos de A y sea
Q qiT
q1T

qnT

(2-15)
una matriz nxn formada por los vectores caractersticos izquierdos de A.

En trminos de estas matrices y considerando las ecuaciones (2-8) y (2-13), se
puede escribir
AP P
(2-16)
QP I , Q P1
(2-17)
y de (2-16) y (2-17)
13
QAP
(2-18)
donde la matriz diag1, 2 ,n
(2-19)
es la matriz de valores caractersticos de A.
2.4
La respuesta en el tiempo de un sistema lineal
La solucin en el dominio del tiempo del sistema puede ser analizada por medio
de la transformada inversa de Laplace. Arreglando la ecuacin (2-3) tenemos
X (s) SI A x0 PP1SI PP1 APP1

1
x0
(2-20)
que tambin puede escribirse como
X (s) P(SI P 1 AP) P1
(2-21)
x0
y de (2-17) y (2-18)
X (s) P(SI )1 P1x0
(2-22)
como P y P-1 son las matrices con los vectores caractersticos derechos e
izquierdos, y (SI-) es una matriz diagonal, se tiene
X ( s)
i 1
pi (qiT x0 )
s i
(2-23)
donde al aplicar la transformada inversa de Laplace se obtiene la respuesta en

el tiempo del sistema en funcin de eigenvalores y eigenvectores.
14
x(t )
q x e p
n
i 1
T
i
(2-24)
Ntese que qTi x0 es una cantidad escalar.
Un modo de oscilacin en particular est excitado si la respuesta solo incluye

un valor caracterstico i.
x(t ) ei t c
(2-25)
donde c es un vector constante proporcional al eigenvector pi.

En la ecuacin (2-24) se observa que el eigenvector izquierdo qTi determina la
influencia de las condiciones iniciales x0 sobre el modo i, mientras que el
eigenvector derecho pi determina la actividad de cada variable de estado en, el
i-esimo modo, es decir la forma del modo. Adems se observa que la respuesta
en el tiempo del sistema en general est formada por una combinacin de todos
los modos de oscilacin del sistema.
Si bien los eigenvectores dan informacin importante de la magnitud y la forma
del modo, la contribucin de cada uno es dependiente de la escala, es decir de
las unidades en que se miden las variables de estado, por lo que la informacin
es til cuando las variables de estado se expresan en las mismas unidades. Los
eigenvectores son particularmente valiosos para determinar los aspectos de
fase y la forma del modo.
En general la respuesta en el tiempo del sistema es de carcter no lineal, est
dada por la sumatoria de varias funciones senoidales de diferente fase y
frecuencia, correspondientes a los modos de oscilacin del sistema. Sin
embargo, en algunos casos un modo de oscilacin puede dominar el
15
comportamiento oscilatorio del sistema debido a su mayor participacin en ese

modo.
2.4.1 Matriz de participacin
Un problema al usar los eigenvectores derechos e izquierdos para identificar las
relaciones entre los estados y los modos es el escalamiento de los elementos
de los eigenvectores. La matriz de participacin es comnmente utilizada para
evitar este problema de escalamiento. Considerando que el elemento pki mide la
actividad de la variable de estado xk en el i-esimo modo y qki da un factor de
peso de esta actividad en el modo, entonces es posible definir
Pki qki pik
(2-26)
como la matriz de participacin y los factores de participacin Pki son nmeros

adimensionales y miden la participacin neta del k-esimo estado en el i-esimo
modo y viceversa. Los factores de participacin, en vista de la normalizacin de
los eigenvectores, tienen la particularidad de que, ya sea asociados con algn
modo o con alguna variable de estado, son igual a la unidad, es decir
n
P
i 1 ki
2.5
1 y
P
k 1 ki
(2-27)
Conclusiones
La respuesta en el tiempo de un sistema dinmico puede describirse en funcin

de los modos de oscilacin del sistema. El anlisis modal, el cual se realiza a
partir de definir el estado del sistema en una condicin inicial, comprende el
clculo de valores y vectores caractersticos y a partir de estos ltimos se
calculan los factores de participacin.
Los valores caractersticos contienen informacin referente a las frecuencias
naturales de oscilacin, los vectores caractersticos contienen informacin sobre
16
la fase y participacin relativa de las variables de estado en cada modo y los

factores de participacin miden la participacin neta de cada variable de estado
en cada modo.
En un sistema con n modos de oscilacin y ante alguna perturbacin arbitraria,
uno o ms de ellos son excitados. Si se desea estimular solo un modo es
suficiente ajustar el vector de condiciones iniciales x0, con un vector
proporcional al eigenvector correspondiente al modo de inters.
17
CAPITULO 3
Interpretacin de los modos de
oscilacin en sistemas elctricos de
potencia
3.1
Introduccin
Los modos de oscilacin en SEPs definen la forma natural en que los rotores
de los generadores oscilan al perturbarse el equilibrio del sistema. Ellos estn
asociados a las inercias de los generadores (ecuacin de oscilacin) y a los
pares mecnicos y elctricos de turbinas y generadores. Si bien, las ecuaciones
que rigen el comportamiento dinmico del SEP son en general no-lineales, la
teora presentada en el captulo 2 para el anlisis de Sistemas Lineales, se
aplica en este captulo para el anlisis de oscilaciones electromecnicas en
SEPs.
En este anlisis se deja fuera, por lo pronto, el efecto del amortiguamiento en
las oscilaciones ya que se analiza en el captulo 4. Los generadores son
representados por el modelo clsico con encadenamientos de flujo constante y
se realizan algunas consideraciones adicionales para obtener los modos de
oscilacin naturales y formas caractersticas de oscilacin del SEP mediante el
clculo de valores y vectores caractersticos. Se incluyen algunos ejemplos
ilustrativos donde se interpreta los resultados del anlisis modal y con la ayuda
18
de un programa de simulacin se observan los modos de oscilacin en SEPs

de prueba.
3.2
Frecuencias naturales no amortiguadas
Las frecuencias naturales no amortiguadas que existen en un sistema son

descritas por un sistema de ecuaciones diferenciales de la forma
(3-1)
x Ax 0
para llegar a esta forma las ecuaciones de aceleracin que describen la

dinmica de las mquinas, se escriben sin el trmino del amortiguamiento [7].
2 H i
i PTi PEi
w0
(3-2)
Se considera que en el momento de la perturbacin, la potencia de la turbina

permanece constante, el modelo clsico (un voltaje constante detrs de una
reactancia transitoria) se emplea para representar a los generadores y el
sistema es linealizado en un punto de operacin.
3.2.1 Obtencin del modelo linealizado de las ecuaciones de flujos de
potencia
La linealizacin del sistema se logra considerando que la posicin angular de
rotores est dada por un valor inicial o de estado estable ms una desviacin
incremental
i i 0 i
(3-3)
al obtener la segunda derivada considerando que i 0 es cero, i i ,

entonces la ecuacin de oscilacin puede ser escrita como:
19
2 H i
i PTi PEi
w0
(3-4)
siendo H la constante de inercia del rotor del grupo turbina generador.

Si las cargas se modelan como una impedancia constante y la red se reduce
hasta los nodos internos de los generadores, la generacin elctrica para cada
mquina puede ahora ser escrita como:
PEi Ei2Gii
E E B sen
Ei E j Gij cos i j
j 1, j i
ij
(3-5)
j 1, j i
la ecuacin (3-5) est en funcin de las magnitudes de voltajes internos de

generadores, ngulos de rotores y admitancias de transferencia; Gij y Bij son las
componentes real e imaginaria de la matriz de admitancias reducida y G ii es la
componente real del trmino diagonal.
Esta ecuacin no-lineal para la potencia elctrica generada se linealiza
utilizando el argumento en (3-3) de tal modo que
i j i 0 i j 0 j
(3-6)
y adems con las aproximaciones siguientes
) sen
cos(
cos( i j ) cos i 0 j 0 i j sen( i 0 j 0 )

sen( i j
i0
j0
i0
j0 )
(3-7)
se obtiene la ecuacin linealizada para la aceleracin del rotor, la cual puede

escribirse como:
20
2 H i
i PTi Ei2Gii
w0
E E B sen(
n
ij
i0
E E G cos(
n
ij
i0
j 0 ) i j sen( i 0 j 0 )
j 1 j i
j 0 ) i j cos( i 0 j 0 )
(3-8)
j 1 j i
como en estado estable no debe existir aceleracin, se cumple que P Ti=PEi0.

por tanto (3-8) se simplifica a
n
2 H i
Ei E j Gij i j sen( i 0 j 0 )
i
w0
j 1 j i
E E B
i
ij
j cos( i 0 j 0 )
j 1 j i
(3-9)
Si se separan los trminos i, la ecuacin (3-9) se expresa como:
n
2 H i
Ei E j Gij sen( i 0 j 0 ) i
Ei E j Gij sen( i 0 j 0 ) j
i
w0
j 1 j i
j 1 j i
(3-10)
n
Ei E j Bij cos( i 0 j 0 ) j
Ei E j Bij cos( i 0 j 0 ) i
j 1 j i
j 1 j i
Usando (3-10) para cada generador, las ecuaciones diferenciales para el

sistema pueden ser escritas en forma matricial
w0
g b i 2H i i
(3-11)
donde
g ij Ei E j Gij sen( i 0 j 0 )
g ii
n
j 1, j i
(3-12)
E i E j Gij sen( i 0 j 0 )
n
j 1, j i
21
g ij
(3-13)
bij Ei E j Bij cos( i 0 j 0 )
bii
n
j 1, j i
(3-14)
E i E j Bij cos( i 0 j 0 )
n
j 1, j i
bij
(3-15)
Definiendo la constante Mi=2Hi/w0 y dividiendo cada ecuacin por su

correspondiente Mi, se obtiene
A
(3-16)
donde los elementos de la matriz A estn dados por
A aij g ij bij / M i
(3-17)
La matriz de estado A contiene informacin condensada de datos del sistema

como las reactancias de la red, el punto de operacin, inercias de generadores,
modelo de cargas, etc.
Partiendo de la ecuacin (3-16) y considerando la relacin entre la posicin
angular, la velocidad angular y la aceleracin angular el modelo del sistema
puede ser representado por un conjunto de ecuaciones diferenciales de la
forma
i 0
I
i
A i
0 i
( 3-18)
Este sistema tiene la forma general x Ax por lo que algunas caractersticas de

la estabilidad del sistema pueden ser analizadas a partir de los eigenvectores
de la matriz de estado A. Las tcnicas para el anlisis de sistemas lineales son
ahora vlidas para un sistema de potencia con el modelo descrito en (3-18).
22
En el anlisis realizado se ha simplificado el modelo y no se ha incluido el

amortiguamiento de las oscilaciones. En el caso general, el problema es mas
complejo cuando se incluye en el modelo los elementos que contribuyen al
amortiguamiento a las oscilaciones. En tal caso se debe incluir el detalle
necesario para lograr el control de las oscilaciones. Sin embargo, generalmente
se considera, y con justificacin, que la frecuencia de una oscilacin es
primeramente determinada por las inercias de generadores, por el nivel de
generacin, por las reactancias de la red y un promedio de las diferencias
angulares entre mquinas y a cuyo anlisis se ha dedicado gran parte de este
capitulo. En el capitulo 5 se analiza el efecto del amortiguamiento en el clculo
de las frecuencias naturales de un SEP.
3.2.2 El clculo de modos de oscilacin como un problema de valores
caractersticos
Se presenta una forma alterna [1] para el clculo de las frecuencias naturales
del sistema, donde las oscilaciones incrementales del rotor son escritas de tal
forma que:
sen( t )
(3-19)
donde y son vectores de dimensin n, el calculo de modos de oscilacin y

sus formas caractersticas puede ser expresado como un problema de
eigenvalores.
De (3-19) tenemos que
cos( t )
(3-20)
2 sen( t )
(3-21)
substituyendo (3-19) y (3-21) en (3-16)

23
Asen( t ) 2 sen( t )
(3-22)
Si (3-22) es vlida para todo tiempo,
A I 0
(3-23)
debe satisfacerse. Si los valores caractersticos i y los vectores caractersticos

i cumplen lo siguiente,
A i I i 0
(3-24)
las frecuencias naturales estn dadas por

i i
(3-25)
Los vectores caractersticos correspondientes i indican las amplitudes relativas

de las oscilaciones del rotor de los generadores para una frecuencia dada
(forma caracterstica de oscilacin), como puede verse en (3-19) y (3-20) para
la posicin y velocidad angular.
3.3
Interpretacin de modos de oscilacin en sistemas elctricos de

potencia de prueba
Ejemplo 3.1. Sistema Mquina-Bus Infinito.

Se tiene un generador de 300 MW de capacidad nominal, se analiza los modos
del sistema para las condiciones de generacin mnima, media y mxima.
A) Generacin mnima.
La figura 3.1 muestra el sistema de prueba a analizar. El generador G1
entrega 100MW al sistema.
24
Zeq
G1
Bus
Infinito
100Mw
Figura 3.1. Sistema de prueba de Mquina-Bus infinito.
En la tabla 3-1 se presentan resultados del anlisis modal del sistema, se tiene
un modo de oscilacin de frecuencia 1.3569Hz (modo 2). Se tendr una
oscilacin sin amortiguamiento. Los factores de participacin nos indican que la
mayor amplitud en la oscilacin es del generador G1 mientras que el bus infinito
se mantendr sin oscilar. Este factor es de 0.0005 para la velocidad y posicin
angular en el bus infinito y en G1 es de 0.4995.
RESULTADOS DEL ANLISIS MODAL

Modos
Mag
Fase
0.0000i
8.5254i
8.5254
90
1.3569
DATOS DE LA OSCILACIN
D = Cte. de Amortiguamiento
W = Frecuencia angular (rad/s)
f = Frecuencia (Hz)
FORMA CARACTERISTICA DE OSCILACIN

(MAGNITUD Y FASE)
modo 1
FACTORES DE PARTICIPACIN
modo 1
modo 2
w BI
180
0.001
0.4995
0.0005
w1
180
180
0.0005
0.4995
d BI
180
0.001
-90
0.4995
0.0005
d1
180
90
0.0005
0.4995
VAR_GEN
modo 2
DESCRIPCIN DEL CASO: Sistema Mquina_Bus Infinito

Generacin Mnima. PG1=100MW
Tabla 3-1 Anlisis modal del sistema mquina-bus infinito en condiciones de

generacin mnima.
25
Adems la fase del modo nos indica que existe un defasamiento de 90 grados
entre la velocidad y posicin angular para cada mquina.
Velocidad angular y Posicin Angular
378
BI
G1
w (r/s)
377.5
377
376.5
376
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
d(grad)
20
15
10
5
0
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 3.2. A) Oscilacin de la frecuencia. B) Oscilacin de la posicin angular.
En la figura 3-2, se muestran los grficos obtenidos al estimular el modo de

oscilacin del sistema en un programa de estabilidad transitoria. La figura 3.2A,
muestra la oscilacin de la velocidad angular en G1, mientas que en el Bus
infinito no hay cambio, esto se observa tambin para la posicin angular en la
figura 3.2B.
B) Generacin Media.
Modos
Mag
Fase
0
0
0.0000i
8.2776i
0
8.2776
0
90
D = Cte. de Amortiguamiento
0
1.3174 W = Frecuencia angular (rad/s)
f = Frecuencia (Hz)
FORMA CARACTERISTICA DE OSCILACIN

(MAGNITUD Y FASE)
modo 1
VAR_GEN
w BI
w2
d BI
d2
1
1
1
1
180
180
180
180
modo 2
0.001
1
0.001
1
0
180
-90
90
modo 1
modo 2
0.4995
0.0005
0.4995
0.0005
0.0005
0.4995
0.0005
0.4995
DESCRIPCIN DEL CASO: Sistema Mquina_Bus Infinito

Generacin media. PG1=200Mw
Tabla 3-2 Anlisis modal del sistema mquina-bus infinito en condiciones de

generacin media.
26
Para las condiciones de generacin media G1 est entregando al sistema de

transmisin 200 MW. La tabla 3-2 muestra la informacin modal. Se observa
que ocurri una pequea disminucin en la frecuencia de oscilacin, ahora
f=1.3174Hz.
Los dems datos son iguales a la condicin anterior. Esta nueva condicin de

378
w (r/s)
377.5
377
376.5
376
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
40
d(grad)
30
20
BI
G1
10
0
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 3.3. A) Oscilacin de la velocidad angular. B) Oscilacin de la

posicin angular.
operacin del sistema se observa en la figura 3.3B en donde la posicin angular

de G1 oscila entre los 30 grados.
C) Generacin Mxima.
Para este caso la frecuencia de oscilacin disminuye aun ms a un valor de
1.2440Hz, los aspectos de fase y participacin de mquinas permanecen sin
cambio.
27
Ejemplo 3.2. Sistema de dos mquinas finitas

La figura 3.4 muestra un sistema de potencia de dos mquinas unidas por un
enlace de 230 KV. Se realiza un anlisis modal para la condicin de operacin
mostrada.
4
G1
G2
50 Mw
100 Mw
50 Mw
Figura 3.4. Sistema de prueba de dos mquinas.
En la tabla 3-3 vemos que ante alguna perturbacin en el sistema, las variables
elctricas tendrn una frecuencia de oscilacin de 1.5838 Hz y la participacin
de ambos generadores en la oscilacin ser prcticamente igual. Las mquinas
oscilan una en contra de la otra, esto lo apreciamos en la fase del modo 1 y
grficamente en la figura 3.5A.
Modos
1
2
0
0

Magnitud
Fase ()
f
9.9515i
9.9515
90
1.5838 D = Cte. de Amortiguamiento
W = Frecuencia angular (rad/s)
0.0000i
0
0
0
W
f = Frecuencia (Hz)
FORMA CARACTERISTICA DE OSCILACIN (MAGNITUD Y FASE)
VAR_GEN
w1
w2
d1
d2
modo 1
1
0
0.9781
180
1
-90
0.9781
90
DESCRIPCIN DEL CASO
Sistema de dos Mquinas Finitas
modo 2
1
1
1
1
0
0
0
0
modo 1
modo 2
0.2528
0.2472
0.2528
0.2472
0.2472
0.2528
0.2472
0.2528
Tabla 3-3 Anlisis modal de un sistema de prueba de dos mquinas.
28
En la figura 3.5 se presentan los resultados grficos obtenidos en el programa

de estabilidad transitoria al estimular el modo de oscilacin.
Se observa tambin que la potencia elctrica entregada por los generadores al
sistema oscila tambin a la frecuencia del modo.
Posicin angular de rotores
20
15
d(grad)
10
5
0
-5
-10
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
Potencia Elctrica Generada

150
Pe (Mw)
100
G1
G2
50
0
-50
-100
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 3.5. Excitacin del modo de oscilacin en el sistema de prueba de dos

mquinas A) Posicin angular. B) Potencia Elctrica
De acuerdo al anlisis de los modos de oscilacin electromecnicos del SEP,

para el caso de modelar el sistema de prueba con dos generadores, se tiene un
solo modo de oscilacin en vista de que se ha tomado uno de los generadores
como referencia y ante cualquier tipo de perturbacin (1) se tendr una oscilacin
de una misma frecuencia. A manera de comprobacin de lo anterior, se realiza
un ejercicio de observar la respuesta del sistema ante la ocurrencia de una falla.
Se simula una falla trifsica en uno de los nodos del sistema con el fin de
observar y confirmar la frecuencia de oscilacin de las variables. La condicin
de operacin es la mostrada en el sistema de prueba de la figura 3.4.
(1) Valido estrictamente para perturbaciones pequeas
29
En la figura 3.6 se tienen los resultados grficos de la respuesta del sistema

ante la falla. El periodo de la oscilacin (1) es de 0.635s dando una frecuencia de
oscilacin de 1.575 Hz que es aproximadamente igual a la calculada en el
anlisis modal de 1.58Hz.
20
15
G1
G2
10
d(grad)
5
0
-5
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.4
1.6
1.8

150
100
Pe (Mw)
50
T=0.635 seg
0
-50
-100
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
Figura 3.6. Oscilacin de la A) posicin angular y B) potencia elctrica al ocurrir una falla.
Ejemplo 3.3. Sistema multimquina.
G1
G4
G5
G2
G3
Figura 3.7. Sistema de prueba para el caso

multimquina.
(1) El tiempo en segundos de un ciclo de la oscilacin.
30
La figura 3.7 muestra un sistema de potencia multimquina. En el se tienen

cinco generadores enlazados mediante lneas de transmisin de caractersticas
diferentes, en cuyo modelo se han incluido los efectos resistivos y reactivos.
En la tabla 3-4 se muestran datos para los modos de oscilacin del sistema y la
forma caracterstica de oscilacin de la posicin angular. Se observa que G1
oscila en contra de las dems mquinas y es esta mquina la que mayor
participacin tiene en la oscilacin, seguida por las mquinas G5, G2, G4 y G3
respectivamente.

FORMA CARACTERISTICA DE OSCILACIN EN MAGNITUD Y FASE
G1
G2
G3
G4
G5
Modo 1
f1=2.279 Hz
1
90
0.1124
-90
0.136
-90
0.1086
-90
0.6523
-90
Modo 2
f2 = 0 Hz
1
-90
1
-90
1
-90
1
-90
1
-90
Modo 3
f3 = 1.4772 Hz
0.4229
180
0.3366
0
0.4667
0
0.3184
0
1
180
Modo 4
f4=1.1179 Hz
0.0425
90
0.3192
-90
1
90
0.2496
-90
0.0598
90
Modo 5
f 5 = 1.0504 Hz
0.0129
-90
0.6974
-90
0.0772
-90
1
90
0.0173
-90
0.0009
0.0848
0.3738
0.0389
0.0016
0
0.1964
0.0011
0.3024
0.0001
0.3508
G1
0.0074
G2
G3
0.0048
G4
0.0052
0.1319
G5
DESCRIPCIN DEL CASO
Sistema multimquina.
0.0902
0.1502
0.0674
0.1125
0.0797
0.058
0.0612
0.0529
0.041
0.2868
Nota: Valores correspondientes a la oscilacin de la posicin angular
Tabla 3-4 Anlisis modal de un sistema de prueba multimquina.
En la figura 3.8A y 3.8B se muestran la oscilacin de la posicin angular y la

potencia elctrica al excitar el modo 1, confirmando los resultados mostrados en
la tabla 3-4.
31

4
d(grad)
-2
-4
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
Pe (Mw)
50
G1
G2
G3
G4
G5
-50
-100
0
0.1
0.2
0.3
0.4
0.5
t(seg)
0.6
0.7
0.8
0.9
Figura 3.8. Oscilacin de la A) posicin angular y B) potencia elctrica para el modo1 (f=2.28Hz).
Para los modos 3, 4 y 5 se grafica la posicin angular para cada generador en

las figuras 3.9A, 3.9B y 3.9C, respectivamente.
Posicin Angular de rotores
G1
G2
G3
G4
G5
d (grad)
-5
0
0.5
1.5
0.5
1.5
0.5
1.5
d (grad)
5
0
-5
d(grad)
5
0
-5
t (seg)
Figuras 3.9 A, 3.9B y 3.9C. Modos de oscilacin en un sistema multimquina

A) Modo 3, f=1.47 Hz B)Modo 4, f=1.11Hz C)Modo 5, f=1.05Hz
32
3.4 Conclusiones
El anlisis modal es de gran utilidad en el anlisis del comportamiento del
sistema de potencia ya que es posible conocer las frecuencias naturales de
oscilacin que pueden presentarse al ocurrir una perturbacin. Si bien al ocurrir
una falla en el sistema, la respuesta oscilatoria no es una frecuencia pura sino
que se trata de una combinacin de las diferentes frecuencias naturales, en
ocasiones un modo puede dominar el comportamiento del sistema de potencia
o parte de l.
En forma individual se obtiene informacin sobre la fase y participacin relativa
de cada variable de estado en los diferentes modos del sistema, esto es, se
puede predecir la forma caracterstica de oscilacin as como aspectos de
amplitudes de oscilacin de las variables.
Por otra parte, utilizando la informacin obtenida a partir de este anlisis, por
ejemplo la informacin de los vectores caractersticos, como informacin de
entrada en un programa de estabilidad transitoria, es posible observar el
comportamiento oscilatorio de las unidades generadoras, la forma caracterstica
de la oscilacin as como el nivel de participacin de cada unidad generadora,
haciendo de esta manera ms eficiente un estudio de estabilidad.
Adems, usando el anlisis modal, se puede determinar el efecto de la
condicin inicial del sistema en la frecuencia de oscilacin al ocurrir una
perturbacin. En los ejemplos mostrados se observa que la frecuencia de
oscilacin es afectada por la condicin inicial del sistema, es decir el punto de
operacin. En los ejemplos mostrados se observa que a mayor cargabilidad del
sistema se tiene un punto de operacin ms crtico y una frecuencia natural de
oscilacin menor. En la prctica el hecho de tener una frecuencia de oscilacin
menor hace que esta sea ms difcil de detectar o bien que se requiera de un
mayor tiempo de anlisis y capacidad de memoria del equipo de simulacin y/o
registro.
33
El estado del sistema de potencia se define en base a las posiciones relativas

de rotores y a la seleccin de una referencia para las desviaciones angulares.
Por esta razn para un SEP de n generadores es comn hablar de n-1 modos
de oscilacin electromecnicos, ya que en el clculo de los n valores
caractersticos, uno de ellos (f=0 Hz) corresponde al modo no oscilatorio o de
estado estable.
Se confirma que un modo de oscilacin en particular puede ser excitado y en el
SEP esto se logra modificando las condiciones iniciales de mquinas con un
vector proporcional al vector caracterstico derecho correspondiente al modo de
oscilacin que se desea observar.
Si bien, ya hemos mencionado que las ecuaciones que rigen el comportamiento
dinmico del SEP son no-lineales, ellas se linealizan en un punto de operacin
con lo que esta formulacin es til para anlisis de pequeas perturbaciones
donde las variables tienen movimientos reducidos.
34
CAPITULO 4
Amortiguamiento de oscilaciones
4.1
Introduccin
En el anlisis realizado en el captulo 3 se han considerado las frecuencias

naturales de oscilacin no amortiguadas lo cual nos ha permitido simplificar de
alguna manera, el anlisis de oscilaciones electromecnicas, sin embargo,
como ya se ha mencionado, en el sistema elctrico de potencia, existen
diferentes factores que tratan de reducir la amplitud y duracin de la oscilacin.
El efecto de estos factores, entre ellos, los devanados amortiguadores de
mquinas, el devanado de campo, el regulador automtico de voltaje, el
gobernador de velocidad y el comportamiento de cargas, en el amortiguamiento
de oscilaciones es mejor comprendido mediante los conceptos de pares de
sincronizacin y pares de amortiguamiento.
En este captulo se consideran las ecuaciones generales, que incluyen los
pares de sincronizacin y amortiguamiento, en la dinmica de las maquinas
generadoras. Se analizan SEPs de prueba obteniendo los modos de oscilacin
naturales e interpretando los efectos del amortiguamiento. Se realizan
comparaciones respecto del anlisis presentado en el captulo 3 para el caso no
amortiguado. As mismo se detalla el proceso matemtico de la excitacin de
modos de oscilacin en SEPs, se presentan ejemplos ilustrativos de ello y se
presentan conclusiones al respecto.
35
4.2
Coeficientes de sincronizacin y amortiguamiento en SEPs
A una frecuencia dada la mquina desarrolla ambos: par de sincronizacin y par

de amortiguamiento y ellos definen el comportamiento dinmico de la mquina
ante disturbios; el par de sincronizacin se encuentra en fase con el ngulo del
rotor de la mquina, depende del punto de operacin del sistema y est
relacionado con la fuerza restauradora para mantener el equilibrio mientras que
el par de amortiguamiento se encuentra en fase con la velocidad angular de la
mquina y tiene que ver con el crecimiento o reduccin de la amplitud y
duracin de la oscilacin.
La condicin para que el sistema sea estable es que ambos pares sean
positivos, de lo contrario el sistema ser inestable. En orden de considerar los
efectos de los pares de sincronizacin y pares de amortiguamiento, la ecuacin
de oscilacin que describe la dinmica de las mquinas se escribe en su forma
ms general
2 H D

Pm Pe
w0
w0
donde
(4-1)
d
w0
dt
(4-2)
H es la constante de inercia
D = constante de amortiguamiento
En la ecuacin (4-1), la constante de amortiguamiento D incluye en forma
aproximada el amortiguamiento natural del sistema, es decir la friccin, el efecto
de devanados amortiguadores, y el comportamiento de cargas. Por ejemplo un
coeficiente
de
1-3
pu
es
comnmente
usado
para
representar
el
amortiguamiento debido a la friccin del aire de la turbina y el efecto de cargas.

Un valor mucho ms alto, arriba de 25 pu, es usado algunas veces como una
36
representacin
del
amortiguamiento,
cuando
fuentes
importantes
de
amortiguamiento han sido omitidas del modelo de la mquina [8].
4.3
Anlisis de los coeficientes de sincronizacin y amortiguamiento en

un SEP bsico
Para el sistema de la figura 4-1 el cambio en la potencia acelerante dado por el

trmino de la derecha en la ecuacin (4-1), considerando que la potencia
mecnica no cambia al ocurrir
un disturbio, puede escribirse en funcin de las componentes del par
X3
X1
E1
X2
E2
Figura 4-1 Sistema de prueba de dos mquinas.
desarrollado por la mquina:
Pm Pe K1 K 2
donde K1
(4-3)
E1E2
cos 0
X eq
(4-4)
substituyendo (4-3) en (4-1) y arreglando la ecuacin
2 H D

K 2 K1 0
w0
w0
(4-5)
K1 y K2 son los coeficientes de sincronizacin y amortiguamiento de la mquina

y sus controles respectivamente. En el modelo empleado en nuestro anlisis no
37
se incluye a los controles de la mquina, consideramos el coeficiente K 2=0 y

usaremos
la
constante
de
amortiguamiento
D,
para
representar
el
amortiguamiento existente en el sistema. La ecuacin (4-5) puede escribirse

como
2 H D

K1 0
w0
w0
(4-6)
La ecuacin linealizada de aceleracin ecuacin (4-6) es vlida ante pequeas

perturbaciones, sin embargo, nos permite ilustrar la dinmica del sistema
alrededor de un punto de operacin.
Se encuentra que la ecuacin caracterstica de (4-6) es
D
Kw
1 0 0
2H
2H
(4-7)
Las races de esta ecuacin son de la forma
1, 2 jwd
K1w0
D2
D
j
4H
2 H 16 H 2
(4-8)
de lo cual se desprende que
D
4H
wd
(4-9)
K1w0
D2
2H
16H 2
(4-10)
La parte real de la ecuacin (4-8) define el amortiguamiento y la parte

imaginaria define la frecuencia de oscilacin amortiguada, ecuaciones (4-9) y
(4-10) respectivamente.
38
La ecuacin (4-6) es de la forma general

2 2 wn wn2 0
(4-11)
en donde la frecuencia natural no amortiguada es
wn
K 1 w0
2H
(4-12)
y la razn de amortiguamiento queda definida como
1 D
2 2 Hwn
(4-13)
De la ecuacin (4-8) se observa que al introducir el amortiguamiento (D), este

no afecta en gran medida la frecuencia de oscilacin, y cuando no existe
amortiguamiento (D=0), se tiene que wd = wn. Las expresiones para el
amortiguamiento y la frecuencia de oscilacin amortiguada wd pueden
simplificarse a:
wn
(4-14)
wd wn 1 2
(4-15)
En la figura 4-2 se muestra la ubicacin de las races en el plano complejo. El

ngulo de amortiguamiento es
1 2
cos1 tan1
(4-16)
39
el cual queda definido por los valores de amortiguamiento y frecuencia de

oscilacin.
Para que el sistema sea estable las races deben localizarse en el semiplano
izquierdo del plano complejo.
jwd
Figura 4-2 . Ubicacin de las races en el plano complejo.
De este modo para tener amortiguamiento positivo en las oscilaciones la

constante de amortiguamiento D debe ser mayor que cero; si D es negativo se
tiene un crecimiento en las oscilaciones. Para D=0 se tienen oscilaciones que
se mantienen en el tiempo, en estos ltimos dos casos se tiene inestabilidad
producida por los pares de amortiguamiento. Otro tipo de inestabilidad es
producida por los pares de sincronizacin. Esta se tiene cuando el coeficiente
de sincronizacin K1 es negativo, en cuyo caso no se tiene oscilacin sino un
crecimiento sostenido en la diferencia angular.
En la figura 4-3 se ilustra grficamente lo anterior. Se grafica, para diferentes
condiciones de operacin, la posicin angular de rotores en un sistema de
potencia bsico.
40
En la figura 4-3A se observa claramente que la condicin necesaria para que el

sistema sea estable es la existencia de pares de sincronizacin y
amortiguamiento positivos.
A falta de cualquiera de ellos el sistema es
inestable, figuras 4-3B a 4-3F.
1000
500
d (g ra d )
d (g ra d )
20
10
A ) D> 0, K 1> 0
B ) D> 0, K 1< 0
-5 0 0
-1 0 0 0
-1 0
0
100
0.5
1.5
1000
C) D< 0, K 1> 0
500
d (g ra d )
d (g ra d )
50
0
GE N1
0
-5 0
GE N2
-5 0 0
-1 0 0 0
-1 0 0
0
0.5
1.5
1.5
1000
d (g ra d )
20
d (g ra d )
D) D< 0, K 1< 0
10
E ) D= 0, K 1> 0
0
F ) D= 0, K 1< 0
500
0
-5 0 0
-1 0
-1 0 0 0
0
t (s e g )
0.5
1
t (s e g )
Figura 4-3 Efecto de los coeficientes de sincronizacin K1 y

amortiguamiento D en la estabilidad.
En este apartado hemos demostrado, para un sistema bsico, que cuando

existe amortiguamiento en el sistema, se tienen valores caractersticos los
cules dan informacin tanto de la frecuencia natural de oscilacin como de la
presencia de amortiguamiento en el sistema, ecuacin 4-8, donde la parte real
define el amortiguamiento del modo de oscilacin y la parte imaginara define la
frecuencia de oscilacin. Se observa adems que la frecuencia natural de
41
oscilacin no se ve afectada en gran medida por la presencia de

amortiguamiento.
4.4
Obtencin del modelo lineal de un SEP multimquina
Considerando que la derivada del cambio en la posicin angular es el cambio

en la velocidad angular ecuacin (4-2), la ecuacin (4-1) puede escribirse para
cada mquina como:
D
2H i
i
i Pmi Pei
w0
w0
i i
donde
i = 1,2,...n,
(4-17)
(4-18)
n = numero de unidades generadoras
En este caso de tener presente el trmino de amortiguamiento, las mismas

consideraciones
hechas en el apartado 3.2 son consideradas; en estado
estable Pmi = Pei, se linealiza el flujo de potencia en un punto de operacin y la

potencia mecnica no cambia al ocurrir una perturbacin. La ecuacin 3-10 se
transforma a:
2H i
Di
i
w0
w0
Ei E j Bij i j cos( i 0 j 0 )
j 1 j i
E E G
i
ij
j sen( i 0 j 0 )
j 1 j i
(4-19)
en la ecuacin anterior solo se ha agregado el termino que contiene el
amortiguamiento D a la ecuacin 3-10. Factorizando la desviacin incremental
de ngulos obtenemos
42
2H i
Di
i

w0
w0
E E B
n
ij
cos( i 0 j 0 ) Ei E j Gij sen( i 0 j 0 ) i j
j 1 j i
(4-20)
Para una condicin inicial de operacin del sistema de potencia, i0 son
conocidos por lo que el termino entre parntesis cuadrados en ecuacin (4-20),
es una constante. Esta ecuacin puede simplificarse a:
2Hi
Di
i

w0
w0
K
1ij
(4-21)
j 1 j i
Donde
K1ij Ei E j Bij cos( i 0 j 0 ) Ei E j Gij sen( i 0 j 0 )
(4-22)
K1ij es el coeficiente de sincronizacin entre las mquinas i y j y representa el

cambio en la potencia elctrica de la mquina i debido al cambio en la posicin
angular entre las mquinas i y j con los dems ngulos constantes.
Escribiendo la ecuacin (4-21) para cada generador, el sistema de ecuaciones
diferenciales de primer orden puede ser escrito en forma matricial
i K S i K D i
(4-23)
que tambin puede ser escrita como
i K D
I
i
KS
O
i

i
(4-24)
43
KS es una matriz que contiene informacin de los coeficientes de sincronizacin

e inercias de mquinas y KD es una matriz que contiene informacin del
amortiguamiento en el sistema. As mismo es el vector de velocidades
angulares y es el vector de posiciones angulares de rotores. Los pares de
sincronizacin estn en fase con la desviacin de la posicin angular y los
pares de amortiguamiento estn en fase con la desviacin de velocidad angular
ecuacin (4-23). I es una matriz identidad y O es una matriz de ceros.
Hasta ahora hemos llegado al modelo linealizado para representar la dinmica
del SEP ecuacin (4-24), el cual tiene forma general de un conjunto de
ecuaciones diferenciales de primer orden de la forma
x Ax
(4-25)
donde x representa el vector de primeras derivadas de las variables de estado

para la velocidad y posicin angular, respectivamente, contenidas en el vector
x. La matriz
K
A D
I
KS
O
(4-26)
contiene la informacin del estado del sistema y est formada por submatrices
como se ve en (4-23) y (4-26).
Ahora
algunas caractersticas de la estabilidad del sistema puede ser
analizadas a partir de los eigenvectores de la matriz de estado A y todos los

argumentos para el anlisis de sistemas lineales, dados en el captulo dos, son
ahora vlidos para el sistema de potencia cuyo modelo ha sido representado en
(4-24) y en el cual se han incluido ambos coeficientes, de sincronizacin y
amortiguamiento.
44
4.5
Excitacin de modos de oscilacin
Segn lo establecimos en el captulo 2, la respuesta en el tiempo de un sistema

como el de la ecuacin (4-25), puede ser expresada en funcin de los
elementos de los vectores caractersticos y los eigenvalores de la matriz de
estado A; rescribiendo la ecuacin 2-24
x(t )
q x e p
n
i 1
T
i
(4-27)
se observa adems que la respuesta dinmica del sistema, ante alguna

perturbacin, es una combinacin lineal de los modos naturales del sistema.
Para un sistema con n modos de oscilacin, la respuesta en el tiempo puede
ser encontrada a partir de (4-27) como:
x(t ) p1 q1T x0 e 1 t p2 q2T x0 e 2 t p3 q3T x0 e 3 t ... pn qnT x0 e n t
(4-28)
y de las propiedades de los vectores caractersticos (definidas en la seccin

2.3.3) sabemos que
qTj pi 1
0
i j
(4-29)
i j
Usando este argumento y observando la ecuacin (4-28), es fcil identificar que

un modo de oscilacin en particular, digamos el i-simo, puede ser excitado si
el vector de condiciones iniciales Xo es modificado con un vector proporcional al
vector caracterstico derecho Pi correspondiente al modo de inters, esto es,
haciendo Xo=P1 para excitar el modo 1 en (4-28), tendremos una seal de
respuesta de la forma
45
x(t ) p1 q1T p1 e 1 t p2 q2T p1 e 2 t p3 q3T p1 e 3 t ... pn qnT p1 e n t
(4-30)
de lo cual se desprende que la respuesta en el tiempo para las variables del

sistema ser de la forma
x1 (t ) p11 e 1 t
x (t )
p 21 e 1 t
2
x(t )
t
x n (t ) p n1 e 1
(4-31)
es decir, todas las variables de estado del sistema tienen una frecuencia de
oscilacin nica, la correspondiente a uno de los modos de oscilacin, el modo
uno en este caso particular.
4.6
Interpretacin de modos de oscilacin en SEPs de Prueba y anlisis

de amortiguamiento
4.6.1 Interpretacin del amortiguamiento en los modos de oscilacin

En el caso de tener presencia de amortiguamiento al calcular los valores
caractersticos del sistema representado en la matriz de estado A, ecuacin
4-26 son de la forma de la ecuacin 4-8, es decir
i i jwd i
( 4-32)
los cules, como era de esperarse, tienen una parte real y una parte imaginara;
la parte real define el amortiguamiento del modo el cual se ve reflejado en todas
las variables de estado y la parte imaginaria define la frecuencia de oscilacin
amortiguada. Ya que la respuesta en el tiempo de un sistema al estimular uno
de los modos de oscilacin digamos el modo 1 en particular, est descrita por la
46
ecuacin 4-31 y considerando una de las variables de estado se tiene que para
esa variable en particular la respuesta en el tiempo es de la forma:
xi (t ) pi1 e1t pi1 e( 1 j1 ) t pi1 e 1 t e j1 t

donde:
( 4-33)
pi1 contiene informacin de la magnitud y fase para la
respuesta de xi(t) al excitar el modo 1

Esta importante expresin, define una respuesta oscilatoria de frecuencia
definida por la parte imaginaria del valor caracterstico, de amplitud y fase
definidas por la magnitud y fase de pi1 y una envolvente de amortiguamiento
definida por la parte real del valor caracterstico. Cada elemento del vector
caracterstico puede expresarse de la forma
~
pki Aki ki A
ki
( 4-34)
~ tal como lo hacemos en circuitos elctricos

en donde hemos definido el fasor A
como la parte que contiene informacin suficiente sobre la magnitud y fase de

una respuesta oscilatoria[10] en este caso se trata de una respuesta oscilatoria
amortiguada.
La ecuacin 4-33 puede ser escrita ahora como
~ e1t A e j e 1 t e j1 t A e 1 t e j (1t )
xi (t ) A
( 4-35)
Esta es una expresin que define la respuesta en el tiempo de una de las

variables de estado del sistema, la variable i en este caso, al excitar uno de los
modos de oscilacin, el modo 1 en este caso particular; la Amplitud es A, la fase
es , la frecuencia es y el coeficiente de amortiguamiento es .
47
A continuacin se presentan algunos ejemplos ilustrativos del amortiguamiento

en SEPs de prueba.
Ejemplo 4.1. Sistema de dos mquinas finitas (vase tambin el ejemplo
3.2).
Descripcin de la condicin Inicial. La mquina G2 entrega 50MW sin perdidas
a una carga y 50MW a la lnea de transmisin. Ambas mquinas cuentan con
amortiguamiento.
G1
G2
50 Mw
50 Mw
100 Mw
Figura 4-4 Sistema de prueba de dos mquinas con amortiguamiento en

mquinas
La tabla 4-1 muestra los resultados obtenidos del anlisis modal. Se observa
que la frecuencia del modo es de 1.5833Hz, (compare con el caso sin
amortiguamiento presentado en el ejemplo 3.2 donde la frecuencia es de
1.5838Hz).
Las mquinas oscilan una en contra de la otra y la mquina G1 tiene una
participacin ligeramente mayor. Se tiene un coeficiente de amortiguamiento D
= -0.5 por lo que se tendr una oscilacin amortiguada en el caso de que el
sistema sea perturbado.
48
Modos
1
2
-0.25
-0.5

Magnitud
Fase ()
f
D = Coef. de Amortiguamiento
9.9484i
9.9515 91.4395
0
0.0000i
0.5
180
1.5833 W = Frecuencia angular (rad/s)
W
f = Frecuencia (Hz)
FORMA CARACTERISTICA DE OSCILACIN (MAGNITUD Y FASE)
modo 1
VAR_GEN
modo 2
modo 1
modo 2
w1
w2
dI
d2
1
1.4395
1
180
0.2529
0.4945
0.2473
0.5055
0.9781 -178.561
1
180
1
-90
1
0
0.2529
0
0.9781
90
1
0
0.2473
0
DESCRIPCIN DEL CASO
Sistema de Dos Mquinas.
Generacin PG1=-50MW, PG2=100MW y carga de 50MW en G2
Tabla 4-1 Anlisis modal de un sistema de Dos Mquinas con Amortiguamiento
En la figura 4-5 se muestra grficamente el resultado de estimular el modo de

oscilacin en un programa de estabilidad transitoria, en donde se observa una
oscilacin amortiguada.
378
w (r/s)
377.5
377
376.5
376
-0.5
0.5
1.5
2.5
1.5
2.5
t(seg)
d(grad)
20
10
0
-10
-0.5
0.5
1
t(seg)
Figura 4-5 A) Oscilacin de la velocidad y B) Oscilacin de la posicin angular.
49
Ejemplo 4.2. Sistema multimquina del ejemplo 3.3.

El sistema de la figura 4-6 se encuentra operando en estado estable. Se han
incluido efectos de amortiguamiento en los modelos de las mquinas G1 y G5.
Se realiza un anlisis de los modos de oscilacin electromecnicos obteniendo
la informacin que se muestra en la tabla 4-2.
G1
G4
G5
G2
G3
Figura 4-6 Sistema de Prueba Multimquina.
En la Tabla 4-2, se observa que los valores caractersticos son nmeros

complejos con informacin del amortiguamiento y frecuencias naturales de
oscilacin. El amortiguamiento presente en el sistema, se refleja en los modos
1, y 3 pero para los modos 4 y 5 no se tiene amortiguamiento o este es muy
pequeo por lo que en el caso de presentarse una perturbacin que excite
alguno de estos modos oscilatorios, se tendr un problema de inestabilidad.
Para el primer modo de oscilacin (modo1), el anlisis modal realizado nos
indica que la mquina G1 oscilar en contra de las dems unidades
generadoras. Se observa que la mayor participacin en la oscilacin la tiene G1
seguido de G5, G2, G4 y G3 respectivamente. En la tabla 4-2, se observan
tambin los aspectos de fase de la oscilacin. A partir de la tabla 4-2, una
50
interpretacin similar a la descrita para el modo 1, puede tambin hacerse para

cada uno de los modos del sistema.
FORMA CARACTERISTICA DE OSCILACIN EN MAGNITUD Y FASE DE LA POSICIN ANGULAR
Modo 1
Modo 2
Modo 3
Modo 4
Modo 5
-0.4167 +14.3120i
-0.3013 + 0.0000i
-0.3187 + 9.2673i
-0.0022 + 7.0240i
-0.0001 + 6.5998i
f 1 = 2.2778 Hz
f 2 = 0.0000 Hz
f 3 = 1.4749 Hz
f 4 = 1.1179 Hz
f 5 = 1.0504 Hz
G1
38.88
0.9989
0.00
0.4245
-55.68
0.0425
-103.37
0.0129
G2
0.1124
-145.38
0.9964
0.00
0.3376
116.62
0.3187
77.40
0.6975
165.04
166.17
G3
0.1359
-145.57
0.9967
0.00
0.4668
115.02
-101.29
0.0772
166.08
G4
0.1086
-145.34
0.9964
0.00
0.3194
116.84
0.2492
77.46
-14.07
G5
0.6523
-139.98
0.00
-59.32
0.0596
-107.14
0.0173
161.78
G1
0.3511
0.1814
0.0586
0.0009
G2
0.0074
0.3003
0.0616
0.0847
0.1964
G3
0.0048
0.1349
0.0529
0.3742
0.0011
G4
0.0052
0.2249
0.0413
0.0388
0.3024
G5
0.132
0.1605
0.2871
0.0016
0.0001
DESCRIPCIN DEL CASO
Anlisis de Amortiguamiento en un Sistema multimquina.
Tabla 4-2 Informacin obtenida en el Anlisis Modal para un sistema Multimquina.

Se observan valores caractersticos con parte real e imaginaria indicacin de la
presencia de amortiguamiento en el sistema
En el ejemplo 3.3 se obtuvieron las frecuencias naturales siguientes, de este

sistema de prueba, pero sin incluir los efectos del amortiguamiento:
f 1 = 2.2790 Hz
f 2 = 0.0000 Hz
f 3 = 1.4772 Hz
f 4 = 1.1179 Hz
f 5 = 1.0504 Hz
Para el caso de incluir los efectos de amortiguamiento en algunas mquinas, las

frecuencias naturales son:
f 1 = 2.2778 Hz
f 2 = 0.0000 Hz
f 3 = 1.4749 Hz
f 4 = 1.1179 Hz
f 5 = 1.0504 Hz
Comparando ambos resultados anteriores, se confirma que la presencia de

amortiguamiento en el sistema no afecta en gran medida el clculo de las
frecuencias naturales de oscilacin.
51
d(grad)
G3
1
0
G4
-1
G1
G2
G3
G4
G5
G5
G1
G2
-2
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0.9
Pe (Mw)
20
0
G3
G1
G2
G3
G4
G5
G4
-20
G2
-40
-60
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
t(seg)
Figura 4-7 Excitacin del modo 1 de 2.2778Hz. En la figura se muestra la
oscilacin de A) la posicin angular y B) la potencia elctrica.
Posteriormente se aplica una perturbacin en el sistema de prueba en un

programa de simulacin de estabilidad transitoria para excitar un modo de
oscilacin. El resultado grfico se muestra en la figura 4-7 para el modo 1. Al
observar este resultado grfico se confirma lo que habamos obtenido a partir
del anlisis modal en cuanto a la forma caracterstica de oscilacin, niveles de
participacin de mquinas y amortiguamiento del modo. Cabe mencionar que
los aspectos de fase no se observan en la respuesta del sistema ante la
perturbacin aplicada en el Programa de Simulacin de Estabilidad Transitoria.
En la figura 4-8 se muestra la excitacin del modo 5, uno de los modos sin
amortiguamiento, con la finalidad de observar la inestabilidad. Se observa
tambin que debido a que se trata de un modo de ms baja frecuencia, se
requiere de un tiempo de simulacin mayor para observar el proceso oscilatorio.
52
Asimismo a la hora de detectar estas bajas frecuencias se requiere de mayor

nmero de registros, en el caso de monitorear el comportamiento de un
sistema.
d(grad)
G2
G4
G1
G2
G3
G4
G5
G3
0
G1
-5
G5
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.4
1.6
1.8
Pe (Mw)
40
G2
G4
G5
20
0
G3
-20
G1
-40
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
Figura 4-8 Excitacin de modo 5, f=1.05Hz. Se observa que en este modo, el

amortiguamiento presente en el sistema no tiene efecto alguno. Se tienen oscilaciones en
forma permanente.
4.7
Conclusiones
Se modela el sistema de potencia incluyendo los efectos de amortiguamiento

para contemplar el efecto de diversos factores que contribuyen a reducir la
amplitud y tiempo de una oscilacin. Ya que se trata de un modelo linealizado,
esta formulacin es til para anlisis de perturbaciones pequeas y una vez
definido el estado del sistema en funcin de coeficientes de sincronizacin y
amortiguamiento se realiza el anlisis modal.
53
A partir del anlisis modal, el clculo de los valores y vectores caractersticos,

se obtiene informacin valiosa la cual es til para conocer las frecuencias
naturales que conforman el sistema de manera de tener una idea clara de la
respuesta del sistema ante alguna perturbacin previo a la ocurrencia de la
misma. Esta informacin es de gran utilidad antes de iniciar un estudio de
estabilidad ya que es posible usar esta informacin para conocer el
comportamiento oscilatorio de los modos del sistema, los aspectos de amplitud
y fase de la oscilacin de variables, si se cuenta con amortiguamiento para las
diferentes frecuencias naturales, etc. Estos resultados ayudan a tener una
mejor comprensin de los anlisis que se realizan en el SEP.
En este captulo hemos demostrado que cuando existe amortiguamiento en el
sistema, se tienen valores caractersticos los cules dan informacin tanto de la
frecuencia natural de oscilacin como de la presencia de amortiguamiento en el
sistema, donde la parte real define el amortiguamiento del modo de oscilacin y
la parte imaginara define la frecuencia de oscilacin. Tambin hemos mostrado
que para un modo de oscilacin, el amortiguamiento se ve reflejado en todas
las variables de estado, e incluso en otras variables como la potencia elctrica.
Se ha confirmado que la frecuencia natural de oscilacin no se ve afectada en
gran
medida
por la
presencia
de
amortiguamiento
y que
el tener
amortiguamiento positivo necesariamente llevar al sistema a una nueva

condicin de estado estable.
54
CAPITULO 5
Oscilaciones de voltaje, potencia
activa y potencia reactiva
5.1
Introduccin
El anlisis de la estabilidad en un sistema elctrico es una tarea importante en

la planificacin y la operacin del sistema. Los SEPs estn sujetos de manera
permanente a diversas perturbaciones debidos principalmente a la variacin de
la carga, a operaciones normales, cambios en generacin, conexin y
desconexin de elementos y tambin debido a fallas en generadores y lneas de
transmisin. El sistema para que sea estable debe producir simultneamente
fuerzas de sincronizacin que permitan mantener las mquinas unidas
operando en paralelo y fuerzas de amortiguamiento que lleven a una operacin
en estado estable despus de cualquier perturbacin.
Generalmente cuando se perturba un sistema de potencia se produce la
excitacin de algunos modos de oscilacin electromecnicos que se observan
en la dinmica de la posicin angular y en la potencia activa de las mquinas
que participan en esos modos, tambin se detectan en la dinmica de algunos
flujos de potencia activa en el sistema de transmisin.
Las tcnicas de anlisis lineal se han utilizado tradicionalmente para
caracterizar las oscilaciones electromecnicas en los sistemas de potencia,
identificando los generadores que participan en cada modo de oscilacin y
55
tambin su nivel de participacin, siendo este anlisis vlido estrictamente ante

pequeas perturbaciones donde las variables de inters tienen movimientos
reducidos.
El objetivo de este captulo es analizar el impacto de estas pequeas
perturbaciones en otras variables que definen el comportamiento reactivo de los
generadores y los voltajes del sistema. En la primera parte del captulo se
desarrollan las ecuaciones para el voltaje en terminales y la potencia reactiva
generada, se obtienen expresiones en funcin de la posicin angular de rotores
para un sistema bsico, usando este desarrollo se definen las regiones lineales
y no-lineales del comportamiento de las unidades en funcin del punto de
operacin.
Mediante el anlisis de casos ilustrativos se muestra la deformacin de voltajes
y potencias que resultan de la condicin operativa y del movimiento de las
posiciones angulares de rotores de unidades generadoras. El anlisis se
extiende a diversos modos de oscilacin y se confirma el mismo efecto de nolinealidad para un sistema con dos modos de oscilacin. Se incluyen
conclusiones referentes al anlisis y los desarrollos presentados en el trabajo.
5.2
5.2.1
Ecuaciones de voltaje, potencia activa y potencia reactiva

Ecuaciones de voltajes en terminales de generadores
Para desarrollar las ecuaciones de los voltajes terminales se utiliza el sistema

V1
V2
X3
X1
E11
X2
Figura 5.1 Sistema de prueba de

E22
56
dos mquinas.
de prueba mostrado en la figura 1.

La corriente en el sistema est dada por:
1
E1 E 2
j( X 1 X 2 X 3 )
( 5-1)
aplicando la ley de Kirchhoff, se obtienen los voltajes en terminales de las

mquinas, V1 y V2:
V1 E1
( 5-2)
jX 1 I
V2 E1 j X1 X 3 I
1
( 5-3)
substituyendo la ecuacin (5-1), y arreglando las ecuaciones resulta:

1
V1
E1 ( X 2 X 3 ) E2 ( X 1 )
X1 X 2 X 3
V2
E1 ( X 2 ) E2 ( X 1 X 3 )
X1 X 2 X 3
( 5-4)
( 5-5)
y si en (5-4) y (5-5) se aplica la identidad de Euler, los voltajes pueden

expresarse como:
X2 X3
X1
V1
E2 cos( 2 )
E1 cos(1 )
X1 X 2 X 3
X1 X 2 X 3
X 2 X 3
X1
j
E2 sen( 2 )
E1sen(1 )
X1 X 2 X 3
X 1 X 2 X 3
57
( 5-6)
X1 X 3
X2
V2
E1 cos( 1 )
E2 cos( 2 )
X1 X 2 X 3
X1 X 2 X 3
X1 X 3
X2
j
E2 sen( 2 )
E1sen( 1 )
X1 X 2 X 3
X 1 X 2 X 3
( 5-7)
En forma alterna, partiendo de las ecuaciones (5-4) y (5-5), la magnitud de los

voltajes tambin puede expresarse de la manera siguiente:
2
X2 X3 2
X1
2
V1
E2
E1
X1 X 2 X 3
X1 X 2 X 3
2 X1X 2 X 3
E E cos( 1 2 )
X 1 X 2 X 3 2 1 2
2
( 5-8)
2 X1 X 3
X2
2
V2
E1
E2
X1 X 2 X 3
X1 X 2 X 3
2 X 2 X1 X 3
E E cos( 1 2 )
X 1 X 2 X 3 2 1 2
2
( 5-9)
En las ecuaciones (5-6) y (5-7) el voltaje se expresa en su parte real e

imaginaria, con funciones coseno y seno de los ngulos de las mquinas,
respectivamente; en cambio en (5-8) y (5-9) el voltaje est dado por un trmino
constante ms otro que depende de una funcin coseno de la diferencia de
ngulos.
5.2.2
Ecuaciones de potencia activa y potencia reactiva
En el sistema de la figura 5.1, la potencia de generacin para la mquina 1 est

dada por:
( 5-10)
S1 V1 I * P1 jQ1
o bien
( 5-11)
S1 V1 I P1 jQ1
58
Substituyendo en (5-11) la ecuacin (5-1) para la corriente y el conjugado de la

ecuacin (5-4) para el voltaje en terminales de la mquina 1, se obtiene:
E 1 ( X 2 X 3 ) E 2 2 ( X 1 ) E1 1 E 2 2
P1 jQ1 1
(X1 X 2 X 3 )
j ( X 1 X 2 X 3 )
( 5-12)
donde al realizar las operaciones algebraicas pertinentes se obtienen las

ecuaciones para la potencia de generacin, activa y reactiva de la mquina 1:
P1
E1 E 2
sen( 12 )
X1 X 2 X 3
( 5-13)
2
X 2 X 3 2
X1
E
Q1
E
1
2
2 2
X 1 X 2 X 3
X 1 X 2 X 3
X X 2 X 3
E E cos( 12 )
1
2 1 2
X 1 X 2 X 3
( 5-14)
donde 12 1 2
( 5-15)
de manera similar para la mquina 2, las ecuaciones quedan:
P2
E1 E 2
sen( 12 )
X1 X 2 X 3
( 5-16)
2
X1 X 3 2
X2
Q2
E
E
2
2
2 1
X 1 X 2 X 3
X 1 X 2 X 3
X X1 X 3
2
E E cos( 12 )
2 1 2
X 1 X 2 X 3
( 5-17)
Las potencias activas, ecuaciones (5-13) y (5-16), son una funcin senoidal de
la diferencia de los ngulos de mquinas y tienen una amplitud determinada por
59
la magnitud de los voltajes internos de mquinas y reactancias del sistema. Las

potencias reactivas estn dadas por un trmino constante sumado de otro
trmino que depende de una funcin coseno de la diferencia de los ngulos.
5.3
Anlisis de oscilaciones de voltaje y potencia
En el captulo 3 demostramos, con algunas consideraciones, que en un sistema

de
potencia
de
generadores
existen
(n-1)
modos
de
oscilacin
electromecnicos, los cuales pueden ser calculados usando una representacin

clsica de mquinas. Los modos de oscilacin estn directamente relacionados
con las oscilaciones de variables como velocidades angulares (w), ngulos de
rotores () y potencias activas, sin embargo estas oscilaciones producen
oscilaciones en todas las variables del sistema de potencia, entre ellas el voltaje
y la potencia reactiva. El objetivo de este captulo es el anlisis de este
comportamiento.
En el anlisis lineal de modos de oscilacin se considera una variacin senoidal
de la posicin angular de las mquinas a una velocidad angular de oscilacin w,
es decir:
Asen(wt )
( 5-18)
de esta forma si cada mquina que participa en la oscilacin tiene una variacin
angular como en la ecuacin (5-18), la potencia activa ser una funcin no
lineal ya que en las ecuaciones (5-13) y (5-16) se observa que depende del
seno de la diferencia de dos senoidales, sin embargo si el movimiento de los
ngulos se realiza en la regin ms lineal de la funcin seno (figura 5.2), es
decir alrededor de cero grados, la oscilacin de la potencia activa sigue un
comportamiento similar a la posicin angular. De esta forma el modo de
oscilacin es observable tambin en la potencia activa generada.
60
En el caso del voltaje en terminales de generadores, estos se pueden expresar

como la suma de un trmino constante que depende de las magnitudes de los
voltajes internos de mquinas y las reactancias de la red y un trmino que
depende del coseno de la diferencia de los ngulos, los voltajes tendrn en
general un comportamiento no-lineal ante la variacin de los ngulos de rotores.
f(t)=cos(t)
f(t)=sen(t)
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
-0.2
-0.2
REGION
-0.4
LINEAL
-0.4
REGION
-0.6
-0.6
LINEAL
-0.8
-0.8
-1
-200
-1
-100
100
200
-200
-100
0
t(grados)
100
200
t(grados)
Figura 5.2 Regin lineal de las funciones coseno y seno.
Para explicar este comportamiento, en la figura 5.2 se marcan las regiones

lineales de las funciones coseno y seno, la regin lineal en la funcin coseno se
tiene para valores entre 40 y 140 grados, es decir para diferencias angulares
grandes. En cambio para valores cercanos a cero grados el comportamiento es
altamente no-lineal. Como el comportamiento del voltaje y la potencia reactiva
dependen del coseno de una diferencia angular, la respuesta en el tiempo
seguir la oscilacin de ngulos solo si se opera en la regin lineal, en el caso
contrario se
producir una distorsin causada por la no linealidad de la
caracterstica.
Otro aspecto relevante en el anlisis es la forma de la oscilacin de los voltajes
en terminales para ambas mquinas, las ecuaciones (5-8) y (5-9) muestran que
dependen del coseno de la misma diferencia angular, lo cual indica que ambos
61
voltajes en terminales tendrn la misma forma de oscilacin, en cambio la

amplitud puede variar ya que depende de la combinacin, en forma diferente
para cada mquina, de los voltajes internos y reactancias de la red.
A su vez la potencia reactiva de generacin puede ser positiva o negativa en la
condicin inicial ya que el trmino constante en las ecuaciones (5-14) y (4-17)
se compone de una suma algebraica de otros dos miembros, uno positivo y uno
negativo. Adems, depende de la misma diferencia angular que el voltaje en
terminales, solo que el trmino que da la amplitud de oscilacin generalmente
es negativo, lo que har que la oscilacin de la potencia reactiva sea opuesta a
la del voltaje. La conclusin del anlisis realizado es que los voltajes y las
potencias reactivas generadas no oscilan de la misma forma que los ngulos de
las mquinas. A continuacin se presentan algunos casos de prueba para
ilustrar el comportamiento de las variables del sistema ante una perturbacin.
5.4
Ejemplos ilustrativos de oscilaciones de voltaje y potencia
Ejemplo 5.1. La figura 5.3 muestra el sistema en su condicin de operacin. El

flujo de potencia por la lnea de transmisin se encuentra en su nivel mximo
con lo que la diferencia angular es grande. Se pretende observar la linealidad
de voltaje en terminales y potencia reactiva al mismo tiempo la no linealidad en
la potencia activa.
G1
2
G2
300 Mw
300 Mw
Figura 5.3 Sistema de prueba de dos mquinas en condiciones de flujo mximo de potencia.
En la figura 5.4A y 5.4B se grafican respectivamente la posicin y la diferencia

angular entre mquinas, esta ltima oscilando entre 55 y 70.
62
En la figura 5.4C se observa una oscilacin lineal de una frecuencia igual a la

del modo de oscilacin de la posicin angular para el voltaje en terminales. Lo
mismo se observa en la figura 5.4D para la potencia reactiva.
Diferencia angular de rotores
80
40
60
20
d(grad)
d(grad)

60
40
20
-20
0
-40
0
0.5
1.5
0.5
1.5
G1
G2
Voltaje en terminales de Mquinas
Potencia reactiva generada
1.06
120
1.04
Q(Mvar)
Vt(pu)
1.02
1
100
80
0.98
60
0.96
0
0.5
1
t(seg)
1.5
0.5
1
t(seg)
1.5
Potencia electrica generada

400
300
Pe (Mw)
200
100
E
0
-100
-200
-300
-400
0.5
1
t(seg)
1.5
Figura 5.4 Comportamientos de variables del sistema de dos mquinas ante una perturbacin.
En la figura 5.4E se observa la no linealidad en la oscilacin de la potencia

activa generada donde la curva est un poco deformada. Esto se debe a que la
regin lineal de funcin seno se tiene entre los 50 y 50 grados.
Ejemplo 5.2. En este caso el sistema se encuentra descargado al momento de

la perturbacin, vea la figura 5.5. El objetivo es operar el sistema en la regin
63
ms no-lineal para observar el comportamiento de voltajes y potencias

reactivas.
2
G2
G1
0
Figura 5.5 Sistema de prueba operando en la regin no-lineal.
En la figuras 5.6A y 5.6B se observa el comportamiento de los ngulos de

rotores y la diferencia angular en el momento en que se perturba el sistema, las
6
4
2
d(grad)
d(grad)
0
-2
-2
-4
-6
-4
0
0.5
1.5
0.5
1.5
G1
G2
1.02
10
5
Q(Mvar)
Vt(pu)
1.015
1.01
1.005
D
0
-5
-10
0
0.5
1
t(seg)
1.5
0.5
1
t(seg)
1.5
Figura 5.6. A) Posicin angular B) Diferencia angular C) y D) Comportamiento no lineal del

Voltaje y la Potencia Reactiva respectivamente.
64
mquinas oscilan una contra la otra, la diferencia angular se mueve alrededor

de 6 grados. Como resultado de esta no linealidad, en la figura 5.6C se
observa una doble frecuencia en la oscilacin del voltaje de las mquinas. En la
figura 5.6D se muestra el comportamiento de la potencia reactiva que tiene una
forma de oscilacin similar pero opuesta al voltaje.
Ejemplo 5.3. En este caso los generadores alimentan a una carga, el objetivo
es analizar el cambio en la condicin de operacin y observar el
comportamiento resultante.
2
G2
G1
100
200
100
Figura 5.7 Efecto en el cambio en la condicin de operacin en el

sistema de potencia de prueba.
En este caso el sistema est en un punto de operacin donde la diferencia

angular se mueve dentro de la regin no-lineal y est variando entre -14 y 3
grados, figura 5.8B.
Se observa una mayor distorsin para el voltaje en terminales de la mquina 2
(figura 5.8C), ya que la posicin angular de esta mquina est oscilando en la
regin ms no-lineal.
En la figura 5.8D, se observa que en la oscilacin de la potencia reactiva
tambin existe distorsin, aunque, el efecto de la no linealidad es menor que en
el voltaje como se aprecia tambin en las ecuaciones (5-14) y (5-17).
65
15
10
d(grad)
d(grad)
B
5
-5
-10
0.5
1.5
-15
0.5
G1
G2
1.5

30
1.02
1.015
20
Q(Mvar)
1.01
Vt(pu)
-5
1.005
D
10
1
0
0.995
0.99
0
0.5
1
t(seg)
1.5
-10
0.5
1
t(seg)
1.5
Figura 5.8. Efecto de la condicin de operacin en las variables del sistema.
Ejemplo 5.4. En el sistema de tres mquinas de la figura 5.9 se estimulan los

modos de oscilacin por separado para observar el comportamiento de voltajes
y potencias reactivas en cada modo.
1
2
G2
G1
3
G3
Figura 5.9 Sistema de potencia de prueba de tres generadores.
En este caso el sistema est descargado, el objetivo es confirmar el

comportamiento de voltajes y potencias reactivas para esta regin no-lineal.
66
En la figura 5.10A y 5.10B se observa la oscilacin de los ngulos para el

primer modo de oscilacin (f=2.52 Hz). El efecto no-lineal se ve reflejado en la
figura 5.11A y 5.11B, donde se observa una doble frecuencia en la oscilacin de
voltaje y potencia reactiva.
Para el segundo modo de oscilacin (1.63 Hz) se confirma tambin que la
frecuencia de oscilacin del voltaje y la potencia reactiva tambin es el doble de
la oscilacin de los ngulos, figuras 5.12 y 5.13.
Modo 1, f = 2.5214 Hz.

Diferencia angular
4
Difd (grad)
2
0
-2
-4
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.8
G1
G2
G3
d(grad)
-2
-4
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
Figura 5.10 (A y B). Diferencia y posicin angular de mquinas.
67
Modo 1, f = 2.5214 Hz.

1.02
1.0199
Vt(pu)
1.0198
1.0197
1.0196
1.0195
1.0194
0
0.2
0.4
0.6
0.8
1.2
1.4
1.8 G1 2
G2
G3
1.6

0.4
0.3
Q(Mvar)
0.2
0.1
0
-0.1
-0.2
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 5.11 A) Voltaje y B) Potencia reactiva de generadores
Modo 2, f = 1.6314Hz
Diferencia angular
10
Difd (grad)
-5
-10
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.8

G1
G2
G3
d(grad)
-5
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
Figura 5.12 (A y B). Diferencia y posicin angular de mquinas.
68
Modo 2, f = 1.6314Hz
1.02
Vt(pu)
1.0195
1.019
1.0185
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.8
G1
G2
G3

1
Q(Mvar)
0.8
0.6
0.4
0.2
0
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
Figura 5.13 A) Voltaje y B) Potencia reactiva de generadores
Ejemplo 5.5. En este caso el sistema de tres mquinas suministra la potencia

demandada en dos nodos de carga. El objetivo es operar el sistema en una
regin lineal.
Modo 1 f= 2.5240 Hz
Diferencia angular
0
Difd (grad)
-5
-10
-15
-20
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
d(grad)
G1
G2
G3
-5
-10
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 5.14 A) Diferencia y B) Posicin angular de mquinas.
69
Modo 1 f= 2.5240 Hz

1.022
Vt(pu)
1.021
1.02
1.019
1.018
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.4
1.6
1.8

8
7
G1
G2
G3
Q(Mvar)
6
5
4
3
2
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
Figura 5.15 A) Voltaje y B) Potencia reactiva de generadores.
En las figuras 5.14 y 5.16 se observa la respuesta en el tiempo para los ngulos
de rotores ante el disturbio; los voltajes y potencias reactivas oscilan a la
frecuencia de la posicin angular (figuras 5.15 y 5.17), sin embargo hay un
grado de distorsin dada por la no linealidad del coseno.
Para ambos modos de oscilacin se cumplen estas observaciones, adems que
la oscilacin de cada potencia reactiva generada es opuesta al voltaje en
terminales de esa mquina.
70
Modo 2 f=1.6146 Hz
Diferencia angular
0
Difd (grad)
-5
-10
-15
-20
-25
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
G1
G2
G3

10
d(grad)
5
0
-5
-10
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 5.16(A y B). Diferencia y posicin angular de mquinas.

1.03
1.02
G1
G2
G3
1.015
1.01
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.4
1.6
1.8

10
5
Q(Mvar)
Vt(pu)
1.025
-5
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.8
Figura 5.17.A) Voltaje y B) Potencia reactiva de generadores
71
5.5
Conclusiones
Comportamiento de voltajes en terminales

La condicin inicial para los ngulos de los voltajes internos de mquinas
depende del punto de operacin del sistema, al ocurrir una perturbacin se
tendrn oscilaciones en estas variables del sistema. Si la perturbacin es
pequea los ngulos de los voltajes internos de mquinas siguen una variacin
senoidal alrededor de su condicin inicial y oscilan en forma opuesta para el
caso de dos mquinas. Los voltajes en terminales dependen directamente de la
variacin de estos ngulos y por lo tanto la forma de oscilacin es tambin
funcin del punto de operacin.
En general el comportamiento dinmico de los voltajes es de naturaleza nolineal debido a que la oscilacin es descrita por una suma de senos y cosenos
de funciones senoidales. Adems, los voltajes de los generadores oscilan en la
misma direccin ya que dependen principalmente de funciones coseno, en las
cuales, sin importar si crece o decrece la diferencia angular, el resultado
siempre ser positivo.
La distorsin de los voltajes en terminales de mquinas depender del ngulo
del voltaje interno de las mquinas, si se opera en la regin no-lineal la
distorsin produce otras frecuencias de oscilacin en los voltajes, si se opera en
la regin lineal la distorsin ser menor. La conclusin es que los voltajes no
siguen el comportamiento del modo de oscilacin de las posiciones angulares.
Comportamiento de la potencia reactiva
La oscilacin de la potencia reactiva generada es en general no-lineal ya que
est dada por un valor constante sumado o restado de una funcin coseno de la
diferencia de dos funciones senoidales. Sin embargo, esta oscilacin puede
tener un comportamiento bastante lineal, lo que depende directamente del
comportamiento de la diferencia angular.
72
El tener la posibilidad de que el trmino constante pueda ser positivo o negativo

indica que la mquina en la condicin inicial puede estar entregando o
recibiendo potencia reactiva, el trmino que produce la oscilacin puede a su
vez ser positivo o negativo, lo cual describe la forma de oscilacin.
Se debe resaltar que no es importante si la diferencia angular es positiva o
negativa ya que el coseno es una funcin par, lo que s es relevante en la
oscilacin es que si existe cambio de signo en la diferencia angular, que se
tiene cuando el sistema opera en una regin no-lineal, se tiene distorsin en el
comportamiento de la potencia reactiva. Si bien, tanto la potencia reactiva como
el voltaje en terminales
dependen del coseno de la diferencia angular, la
oscilacin de la potencia reactiva es menos sensible debido a la amplitud de

ste trmino (ecuaciones (5-8), (5-9), (5-14) y (5-17)).
La potencia activa es en general una funcin no-lineal cuya oscilacin es
descrita por una funcin seno de la diferencia de los ngulos de mquinas, sin
embargo, se ha demostrado e ilustrado que la oscilacin de la potencia activa
sigue el modo de oscilacin de la posicin angular, esto, en el caso tener
diferencias angulares pequeas. La explicacin es que se opera en la regin
lineal de la funcin seno.
Cuando la diferencia angular aumenta conforme aumenta la condicin de carga
del sistema la oscilacin de la potencia activa es no-lineal tenindose
frecuencias adicionales a la del modo de oscilacin de la posicin angular. La
explicacin es que se opera en la regin no-lineal de la funcin seno.
73
CAPITULO 6
Determinacin de modos de
oscilacin mediante el
procesamiento de respuestas en el
tiempo
6.1
Introduccin
En captulos anteriores hemos comentado que la respuesta en el tiempo de un

sistema lineal ante una perturbacin es una combinacin lineal de todas las
frecuencias naturales de dicho sistema. Por medio del anlisis modal hemos
calculado estas frecuencias naturales (modos de oscilacin) as como los
vectores caractersticos a partir de haber definido el estado del sistema en una
condicin inicial en funcin de parmetros de generadores, caractersticas de la
red elctrica, modelos de cargas, etc.
Adems ha sido posible, para un sistema elctrico de potencia de prueba, lograr
la excitacin de un solo modo de oscilacin y hemos podido observar por
separado, las oscilaciones de las variables de estado del sistema a una sola
frecuencia natural de oscilacin utilizando un programa de simulacin.
En este captulo se presenta una alternativa para el clculo de los modos
naturales de un sistema elctrico de potencia cuando no se conocen por
74
completo los parmetros fsicos de la red de manera que no es sencillo

establecer, con esta base, el estado del sistema en alguna condicin inicial de
operacin y por tanto no pueden utilizarse las tcnicas del anlisis modal.
La alternativa propuesta consiste en la obtencin de los modos naturales de
oscilacin de sistemas elctricos a partir de registros de datos. Esta alternativa
requiere tener registros de datos del comportamiento del sistema de potencia,
para un intervalo de tiempo determinado, obtenidos durante perturbaciones que
se hayan presentado en dicho sistema. Mediante el uso de tcnicas apropiadas,
se procesa la informacin contenida en los registros, para obtener los modos de
oscilacin del sistema en cuestin.
Adems de lo antes descrito, se presentan tambin en ste captulo, ejemplos
ilustrativos de la obtencin e interpretacin de los modos naturales en SEPs de
prueba. Los registros del comportamiento del SEP de prueba, se obtienen con
base en la simulacin de fallas en dicho sistema, con la ayuda de un programa
de estabilidad transitoria. Posteriormente, mediante el adecuado procesamiento
de la informacin contenida en tales registros, de la respuesta en el tiempo del
sistema ante la falla, se calculan los modos naturales del SEP de prueba,
aspectos de fase y factores de participacin y finalmente los resultados
obtenidos se comparan contra los resultados obtenidos en captulos anteriores
a partir de las tcnicas del anlisis modal obteniendo conclusiones relevantes.
6.2
Determinacin del espectro de frecuencias de seales de respuesta

en el tiempo en un sistema elctrico de potencia
La respuesta en el tiempo de un sistema ante una perturbacin es la

superposicin de componentes de diferentes frecuencias. Si se tiene la
respuesta en el tiempo de un sistema, a partir de dicha informacin es posible
determinar el espectro de frecuencias de alguna seal as como los aspectos de
fase. Para ello pueden utilizase tcnicas apropiadas de muestreo y anlisis de
75
seales as como la transformada de Fourier en tiempo discreto aplicada a la

seal de respuesta [21, 16]. Para aplicar esta teora a un SEP debe contarse
con los registros necesarios, procesar adecuadamente las seales, realizar los
clculos necesarios y realizar la interpretacin de los resultados lo cual forma
parte del objetivo de este captulo.
6.3
Aplicacin en sistemas elctricos de potencia de prueba
En el anlisis de SEPs con frecuencia se utilizan esquemas de monitoreo y

medicin de parmetros elctricos (registradores de eventos que incluyen
Unidades de Medicin Fasorial, PMUs por su nombre en ingles) que faciliten el
estudio del comportamiento del sistema ante disturbios. A partir de registros
apropiados y considerando los datos registrados como muestras de una seal
discreta, es posible analizar la seal de respuesta en el tiempo, de modo de
obtener informacin de la magnitud y fase de las frecuencias naturales que
componen a dicha seal.
En el ejemplo siguiente, mediante la aplicacin de tcnicas de procesamiento
digital de seales, tcnicas de muestreo as como la transformada de Fourier en
tiempo discreto se obtienen los modos naturales de oscilacin de un SEP de
prueba.
6.3.1 Determinacin de modos de oscilacin en un SEP de Prueba

Ejemplo 6.1. En este ejemplo, se utiliza el mismo SEP de prueba utilizado en el
ejemplo 3.3 que corresponde a un sistema multimquina de cinco generadores
el cual se muestra en la figura 6-1. Con la ayuda de un programa de estabilidad
transitoria, se modela el sistema de prueba en el cual se registra el
comportamiento de la velocidad angular de rotores de maquinas al ocurrir una
perturbacin. A partir de esa informacin se realiza el anlisis de las seales
76
muestreadas para obtener las frecuencias que componen la respuesta de cada

generador.
El sistema de la figura 6.1 se encuentra operando en estado estable cuando se
presenta una falla trifsica, en uno de los nodos del sistema, que se libera en un
ciclo. Se registra el comportamiento de variables en diversos puntos del sistema
y en la figura 6-2 se muestra grficamente la respuesta oscilatoria de las
mquinas la cual, como es de esperarse, es una respuesta no lineal.
G1
G4
G5
G2
G3
Figura 6-1 Sistema de Prueba Multimquina.
Las seales presentan una componente de directa que es necesario eliminar

previo al procesado de las muestras. Posteriormente se elige una ventana de
muestreo que permita observar mejor las componentes frecuenciales tratando
de evitar que puedan confundirse con los lbulos laterales que se obtienen
durante el procesamiento de la seal.
Con estas y otras aplicaciones del anlisis de seales se obtienen los
resultados que se presentan enseguida, utilizando las tcnicas del anlisis de
Fourier para el clculo del espectro de frecuencias, de las seales de respuesta
del sistema en cuestin.
77
Velocidad angular
376.99
G1
G2
G3
G4
G5
376.98
w (rad/s)
376.97
376.96
376.95
376.94
376.93
376.92
376.91
0
0.5
1.5
t(seg)
2.5
Figura 6-2 Oscilacin de velocidades de mquinas al ocurrir un disturbio.
En la figura 6.3 se muestra, para cada mquina, el espectro de frecuencias

obtenido a partir de la respuesta no lineal de la velocidad de rotores. Se puede
observar las componentes principales de la respuesta de las mquinas.
Tambin se observa que hay frecuencias de oscilacin que se presentan en
varias maquinas por ejemplo, aunque en diferentes proporciones, la frecuencia
de 2.27 Hz est presente en todos los generadores.
Algunas frecuencias tienen una mayor magnitud lo cual indica mayor influencia
de tal frecuencia de oscilacin en una mquina en particular por ejemplo en la
respuesta de G1, la frecuencia de 2.27 Hz tiene la mayor magnitud y domina el
comportamiento oscilatorio de la maquina 1 lo cual tambin se observa en la
figura 6.2 donde dicha respuesta tiene un periodo de 0.45seg. Se observa
tambin como las frecuencias naturales de oscilacin que componen a este
sistema de prueba participan en forma diferente en cada generador.
78
Figura 6-3 Espectro de Frecuencias de las Oscilaciones de mquinas.
La figura 6.4 muestra el espectro obtenido para la mquina 1 donde se observa

tanto la magnitud como la fase de las frecuencias que componen la respuesta
de la mquina ante la falla.
79
Figura 6-4 Amplitud y

Fase de componentes de
la
mquina 1. Para
f=2.27Hz la fase es de 173.1, para f=1.47Hz la
fase es 8 y para
f=1.11Hz la fase es de 169.5.
En la tabla 6-1 se muestran las frecuencias, en magnitud y fase, que componen

la respuesta en el tiempo de cada generador. Se observa cual de dichas
frecuencias tiene mayor influencia en el comportamiento de una mquina en
particular, esto es, la componente de mayor magnitud, siendo 2.27 Hz en G1,
1.04 Hz en G2, 1.11 Hz en G3, 1.05 Hz en G4 y 2.27 Hz en G5.
G1
G2
G3
Frec
Mag
Fase
Frec
Mag
Fase
Frec
Mag
Fase
2.268
1.473
1.110
4.975
0.335
0.228
-173.1
8.0
-169.5
1.040
1.120
2.268
1.475
1.791
1.576
0.553
0.266
15.1
-11.4
6.9
178.7
1.113
2.268
1.473
5.383
0.665
0.361
-176.2
6.9
-171.7
G4
G5
Frec
Mag
Fase
Frec
Mag
Fase
1.053
2.268
1.473
3.055
0.530
0.249
166.9
6.9
-171.7
2.268
1.473
1.110
3.281
0.778
0.318
6.9
8.0
-169.3
Tabla 6-1 Espectro de frecuencias de la respuesta oscilatoria de los generadores G1,

G2, G3, G4 y G5 ante la ocurrencia de un disturbio.
80
Es importante el hecho de que los resultados de la tabla 6-1, pueden

compararse entre todas las mquinas, por ejemplo, la frecuencia de 1.11 Hz es
la que presenta la mayor amplitud de oscilacin ante el disturbio con una
magnitud relativa de 5.83 en el generador G3.Estos resultados son
particularmente tiles para efectos de comparacin de magnitudes de
oscilacin. Se observa tambin que son 4 las frecuencias que, en mayor o
menor medida, contribuyen a la respuesta no lineal de todas las maquinas del
sistema y ellas son las frecuencias naturales del sistema, tal como se aprecia
en la tabla 6-3.
6.3.2 Comparacin de resultados

El SEP de prueba utilizado en el ejemplo 6.1 es el mismo que se utiliz en el
ejemplo 3.3. De este sistema se conocen sus parmetros de red y en el captulo
3 se calcularon, a partir del anlisis modal, sus frecuencias naturales. A
continuacin se realiza una comparacin de los resultados obtenidos al analizar
el comportamiento del sistema ante el disturbio y su relacin con las frecuencias
naturales de oscilacin obtenidas a partir del anlisis modal.
En la tabla 6-2 la informacin de la tabla 6-1 se ha reacomodado y normalizado
para fines de comparacin con los resultados obtenidos en el clculo de modos
del sistema en el ejemplo 3.3 los cuales se muestran en la tabla 6-3.
Los resultados de las tablas 6-2 y 6-3, muestran la magnitud y fase de la
participacin relativa de cada mquina en cada modo de oscilacin. De este
modo se puede observar cual generador tiene mayor participacin en una
frecuencia dada, por ejemplo G1 en 2.27 Hz, G5 en 1.47 Hz, G3 en 1.11 Hz y
G4 en 1.05 Hz.
81
Tambin se puede observar en la tabla 6-2, la forma caracterstica de oscilacin

para algn modo en particular, por ejemplo en el modo de 2.27 Hz, G1 oscila en
contra de los dems generadores, tal como se aprecia tambin en la tabla 6-3.
RESULTADOS DEL ANLISIS ESPECTRAL (FORMA CARACTERISTICA DE OSCILACIN EN

MAGNITUD Y FASE)
f = 2.27 Hz
f = 1.47 Hz
f =1.11Hz
G1
1.000
-173.13
0.430
8.01
0.042
-169.50
G2
0.111
6.87
0.342
178.67
0.293
-11.40
G3
0.134
6.88
0.464
-171.7
1.000
-176.23
G4
0.107
6.86
0.319
-171.75
G5
0.659
6.87
1.000
8.00
0.059
f = 1.05 Hz
0.586
15.1
1.000
166.86
-169.3
DESCRIPCIN DEL CASO: Sistema multimquina.

Valores correspondientes a la oscilacin de la velocidad angular
Tabla 6-2 Espectro de frecuencias de la respuesta oscilatoria de los generadores G1, G2, G3,
G4 y G5 ante la ocurrencia de un disturbio. Las frecuencias se han ordenado para indicar en
cuales generadores est presente cada frecuencia. Por ejemplo la frecuencia de 1.05 Hz solo
est presente en G2 y G4.

FORMA CARACTERISTICA DE OSCILACIN EN MAGNITUD Y FASE
Modo 1
Modo 2
f1=2.279 Hz
f2 = 0 Hz
G1
1
90
1
-90
G2 0.112
-90
1
-90
G3 0.136
-90
1
-90
G4 0.109
-90
1
-90
G5 0.652
-90
1
-90
DESCRIPCIN DEL CASO
Modo 3
f3 = 1.4772 Hz
0.4229
180
0.3366
0
0.4667
0
0.3184
0
1
180
Modo 4
f4=1.1179 Hz
0.0425 90
0.3192 -90
1
90
0.2496 -90
0.0598 90
Modo 5
f 5 = 1.0504 Hz
0.0129
-90
0.6974
-90
0.0772
-90
1
90
0.0173
-90
Sistema multimquina.
Tabla 6-3 Clculo de modos de oscilacin del ejemplo 3.3
82
Al comparar ambos resultados mostrados en las tablas 6-2 y 6-3 observamos

que, el disturbio ocurrido en el sistema provoc una respuesta oscilatoria en las
mquinas, compuesta por las frecuencias naturales del sistema en cuestin. En
la tabla 6-2 se observa que, con ligeras variaciones, la participacin relativa de
las maquinas ocurri de acuerdo a lo que pudiera haberse previsto de la tabla
6-3.
A partir de esta comparacin tambin es importante mencionar que la matriz de
participacin utilizada en el anlisis lineal para mostrar resultados de
participacin neta de cada variable de estado en cada modo de oscilacin
puede compararse con el espectro de frecuencias mostrado en la tabla 6-1. En
la tabla 6-4 se muestran los factores de participacin obtenidos en el anlisis
lineal.
Factores de Participacin. Sistema multimquina

f1=2.27 f3=1.47 f4=1.11 f5=1.05
G1
0.349
0.060
0.001
0.000
G2
0.007
0.063
0.088
0.192
G3
0.005
0.052
0.374
0.002
G4
0.005
0.041
0.036
0.306
G5
0.134
0.285
0.002
0.000
Tabla 6-4 Factores de participacin para la oscilacin de la velocidad
angular de mquinas. Resultados del anlisis lineal.
Leyendo en forma horizontal, observamos en orden de magnitud que 2.27 Hz,

1.47 Hz, 1.11 Hz son las frecuencias que participaran en la oscilacin de G1.
De forma similar, en la tabla 6-1 los resultados muestran que la respuesta de
G1 ante el disturbio se compone de las frecuencias 2.27 Hz, 1.47 Hz y 1.11 Hz
en ese orden de magnitud.
83
Si continuamos con la comparacin pero utilizando la informacin de los

factores de participacin, ahora presentada en la tabla 6-5, donde solamente se
cambia el arreglo de la tabla 6-4, podemos comparar directamente con los
resultados del anlisis espectral (tabla 6-1). En ambas tablas 6-1 y 6-5, se
muestran las frecuencias que participan en la respuesta de cada generador.
Factores de participacin. Sistema multimquina.
G1
Frec
2.27
1.47
1.11
1.05
Mag
0.349
0.060
0.001
0.000
G2
G3
G4
Frec
Mag
Frec
Mag
Frec
Mag
1.04
0.192
1.11
0.374
1.05
0.306
1.12
0.088
1.47
0.052
1.47
0.041
1.48
0.063
2.27
0.005
1.11
0.036
2.27
0.007
1.05
0.002
2.27
0.005
G5
Frec
1.47
2.27
1.11
1.05
Mag
0.285
0.134
0.002
0.000
Tabla 6-5 Factores de participacin para la oscilacin de la velocidad

angular de mquinas. Resultados del anlisis lineal. Cambia el arreglo de
la informacin de la Tabla 6-4.
Adems con este arreglo es sencillo observar que, junto con los aspectos de
fase ya presentados (tabla 6-3) es posible predecir el comportamiento no lineal
de la respuesta oscilatoria de las maquinas.
Para el resto de los generadores no se cumple del todo este orden sin embargo
es conveniente mencionar que no se espera que los resultados obtenidos sean
exactamente iguales, pero si muy aproximados, ya que en el caso presentado
en el ejemplo 6.1, los registros analizados corresponden a la respuesta no lineal
del sistema ante una falla trifsica liberada en un ciclo. Es conveniente
mencionar que el sistema de potencia est sujeto a diferentes disturbios
pudiendo cada uno de ellos excitar en forma diferente los modos oscilatorios.
Sin embargo se ha demostrado que es posible, con est tcnica, obtener las
frecuencias naturales de oscilacin de un sistema, as como, con una muy
buena aproximacin, los aspectos de magnitud relativa y fase de la
participacin de las unidades generadoras en la oscilacin.
84
6.4
Conclusiones
En este captulo se han calculado las frecuencias que componen la respuesta

no lineal de un SEP de prueba a partir de registros del comportamiento de
variables de unidades generadoras ante un disturbio ocurrido en el SEP.
Se han comparado los resultados aqu obtenidos con los resultados obtenidos
en el captulo 3 mediante las tcnicas del anlisis modal, comprobando que la
respuesta no lineal de las unidades generadoras del SEP de prueba, est
compuesta por las frecuencias naturales del sistema.
Con el procesamiento de las seales de respuesta de las mquinas, ha sido
posible obtener con muy buena aproximacin, los aspectos de magnitud y fase
de las frecuencias naturales del sistema. Tambin ha se ha podido observar,
cual modo de oscilacin tiene mayor influencia en el comportamiento de una
mquina en particular, informacin que puede ser de utilidad en la identificacin
y solucin de problemas de oscilaciones pobremente amortiguadas. Tambin se
pudo observar cual generador tiene mayor participacin en una frecuencia
dada.
Con los resultados obtenidos se confirma que, con el conocimiento de los
modos naturales del sistema, es posible predecir con magnitudes relativas y
con buen acercamiento el comportamiento del sistema ante disturbios.
Se demuestra que, con la adecuada seleccin de registros, es posible obtener
las frecuencias naturales de oscilacin de un sistema, mediante el adecuado
procesamiento de las seales de la respuesta en el tiempo obtenidas a partir de
disturbios ocurridos en dicho sistema.
Si en un sistema de potencia se piensa que hay problemas de oscilaciones de
baja frecuencia se puede implementar el monitoreo del mismo para estudiar su
comportamiento ante disturbios; a partir de dicho comportamiento se pueden
85
identificar las frecuencias causa del problema para posteriormente proporcionar

en el sistema el amortiguamiento de las mismas.
86
CAPITULO 7
Conclusiones y recomendaciones
7.1
Conclusiones finales
El anlisis de los modos de oscilacin de Sistemas Elctricos de Potencia

(SEP) (anlisis modal), se realiza a partir de definir el estado del sistema en una
condicin inicial de operacin y comprende el clculo de valores y vectores
caractersticos y a partir de estos ltimos, se calculan los factores de
participacin.
En los resultados obtenidos en SEPs de prueba se comprueba que los valores
caractersticos contienen informacin referente a las frecuencias naturales de
oscilacin y del amortiguamiento del sistema, que los vectores caractersticos
contienen informacin tanto de la influencia de la condicin inicial en los modos
como de la fase y participacin relativa de las variables de estado en cada
modo y que los factores de participacin miden la participacin neta de cada
variable de estado en cada modo.
Mediante el anlisis modal aplicado a un SEP es posible predecir las
frecuencias naturales de oscilacin que pueden presentarse al ocurrir una
perturbacin. Los resultados muestran que ante una perturbacin, la respuesta
oscilatoria de las variables del sistema ser una combinacin de las diferentes
frecuencias naturales de dicho sistema.
87
Con la ayuda de un programa de estabilidad transitoria se realiza la

estimulacin de modos de oscilacin en particular en el SEP. Se observa el
comportamiento oscilatorio de unidades generadoras, forma caracterstica de
oscilacin y nivel de participacin de cada unidad, mediante el uso de la
informacin del anlisis modal como informacin de entrada al programa de
simulacin;
esta
informacin
permite
predecir
en
cierta
medida
el
comportamiento oscilatorio del sistema ante disturbios.

Mediante el anlisis modal, se determina que la condicin inicial del sistema
tiene un efecto en la frecuencia de oscilacin al ocurrir una perturbacin. En los
resultados obtenidos se observa que a mayor cargabilidad del sistema se tiene
un punto de operacin ms crtico y una frecuencia natural de oscilacin menor.
En la prctica el hecho de tener una frecuencia de oscilacin menor hace que
esta sea ms difcil de detectar y que se requiera de un mayor tiempo de
anlisis y capacidad de memoria del equipo de simulacin y/o registro.
El estado del SEP se define en base a las posiciones relativas de rotores y a la
seleccin de una referencia para las desviaciones angulares. Por esta razn
para un SEP de n generadores es comn hablar de n-1 modos de oscilacin
electromecnicos, ya que en el clculo de los n valores caractersticos, uno de
ellos (f=0 Hz) corresponde al modo no oscilatorio o de estado estable.
Se confirma que un modo de oscilacin en particular puede ser excitado y en
las simulaciones de estabilidad transitoria del SEP esto se logra modificando las
condiciones iniciales de mquinas con un vector proporcional al vector
caracterstico derecho correspondiente al modo de oscilacin que se desea
observar. En SEPs reales los modos oscilatorios son excitados por las
perturbaciones que se presentan.
Las ecuaciones que rigen el comportamiento dinmico del SEP son no-lineales,
por lo que en la tcnica del anlisis modal, ellas se linealizan en un punto de
88
operacin por lo que esta formulacin es til para anlisis de pequeas

perturbaciones donde las variables tienen movimientos reducidos.
Se modela el sistema de potencia incluyendo los efectos de amortiguamiento
para contemplar el efecto de diversos factores que contribuyen a reducir la
amplitud y tiempo de una oscilacin. Por tratarse de un modelo linealizado, esta
formulacin es til para anlisis de perturbaciones pequeas y una vez definido
el estado del sistema en funcin de coeficientes de sincronizacin y
amortiguamiento se realiza el anlisis modal.
Se demuestra que cuando existe amortiguamiento en el sistema, se tienen
valores caractersticos los cules dan informacin tanto de la frecuencia natural
de oscilacin como de la presencia de amortiguamiento en el sistema, donde la
parte real define el amortiguamiento del modo de oscilacin y la parte
imaginara define la frecuencia de oscilacin. Tambin se muestra que para un
modo de oscilacin amortiguado, el amortiguamiento se ve reflejado en todas
las variables de estado, e incluso en otras variables como la potencia elctrica.
Se confirma que la frecuencia natural de oscilacin no se ve afectada en gran
medida por la presencia de amortiguamiento y que el tener amortiguamiento
positivo necesariamente llevar al sistema a una nueva condicin de estado
estable.
La condicin inicial para los ngulos de los voltajes internos de mquinas
depende del punto de operacin del sistema y al ocurrir una perturbacin se
tienen oscilaciones en estos ngulos. Si la perturbacin es pequea los ngulos
de los voltajes internos de mquinas siguen una variacin senoidal alrededor de
su condicin inicial. Los voltajes en terminales dependen directamente de la
variacin de estos ngulos y por lo tanto la forma de oscilacin es tambin
funcin del punto de operacin.
89
En general el comportamiento dinmico de los voltajes es de naturaleza nolineal debido a que la oscilacin es descrita por una suma de senos y cosenos
de funciones senoidales. Adems, los voltajes de los generadores oscilan en la
misma direccin ya que dependen principalmente de funciones coseno, en las
cuales, sin importar si crece o decrece la diferencia angular, el resultado
siempre ser positivo.
La distorsin de los voltajes en terminales de mquinas depender del ngulo
del voltaje interno de las mquinas, si se opera en la regin no-lineal la
distorsin produce otras frecuencias de oscilacin en los voltajes, si se opera en
la regin lineal la distorsin ser menor. La conclusin es que los voltajes no
siguen el comportamiento del modo de oscilacin de las posiciones angulares.
La oscilacin de la potencia reactiva generada es en general no-lineal ya que
est dada por un valor constante sumado o restado de una funcin coseno de la
diferencia de dos funciones senoidales. Sin embargo, esta oscilacin puede
tener un comportamiento bastante lineal, lo que depende directamente del
comportamiento de la diferencia angular.
El tener la posibilidad de que el trmino constante pueda ser positivo o
negativo indica que la mquina en la condicin inicial puede estar entregando o
recibiendo potencia reactiva, el trmino que produce la oscilacin puede a su
vez ser positivo o negativo, lo cual describe la forma de oscilacin. Se debe
resaltar que no es importante si la diferencia angular es positiva o negativa ya
que el coseno es una funcin par, lo que si es relevante en la oscilacin es que
si existe cambio de signo en la diferencia angular, que se tiene cuando el
sistema opera en una regin no-lineal, se tiene distorsin en el comportamiento
de la potencia reactiva. Si bien, tanto la potencia reactiva como el voltaje en
terminales dependen del coseno de la diferencia angular, la oscilacin de la
potencia reactiva es menos sensible debido a la amplitud de ste trmino
(ecuaciones (5-8), (5-9), (5-14) y (5-17)).
90
La potencia activa es en general una funcin no-lineal cuya oscilacin es

descrita por una funcin seno de la diferencia de los ngulos de mquinas, sin
embargo, se ha demostrado e ilustrado que la oscilacin de la potencia activa
sigue el modo de oscilacin de la posicin angular, esto, en el caso tener
diferencias angulares pequeas. La explicacin es que se opera en la regin
lineal de la funcin seno. Cuando la diferencia angular aumenta conforme
aumenta la condicin de carga del sistema la oscilacin de la potencia activa es
no-lineal tenindose frecuencias adicionales a la del modo de oscilacin de la
posicin angular. La explicacin es que se opera en la regin no-lineal de la
funcin seno.
Se calculan las frecuencias que componen la respuesta no lineal de un SEP de
prueba a partir de registros del comportamiento de variables de unidades
generadoras ante un disturbio ocurrido en el SEP. Se comparan los resultados
aqu obtenidos, con los resultados obtenidos mediante las tcnicas del anlisis
modal, comprobando que la respuesta no lineal de las unidades generadoras
del SEP de prueba, est compuesta por las frecuencias naturales del sistema.
Con el procesamiento de las seales de respuesta de las mquinas ante una
falla, ha sido posible obtener con muy buena aproximacin, los aspectos de
magnitud y fase de las frecuencias naturales del sistema. Tambin ha se ha
podido observar, cual modo de oscilacin tiene mayor influencia en el
comportamiento de una mquina en particular, informacin que puede ser de
utilidad en la identificacin y solucin de problemas de oscilaciones pobremente
amortiguadas.
Con los resultados obtenidos a partir de registros de disturbios, se confirma
que, con el conocimiento de los modos naturales del sistema, es posible
predecir con magnitudes relativas y con buen acercamiento el comportamiento
del sistema ante disturbios.
91
Se demuestra que, con la adecuada seleccin de registros, es posible obtener

las frecuencias naturales de oscilacin de un sistema, mediante el adecuado
procesamiento de las seales de la respuesta en el tiempo obtenidas a partir de
disturbios ocurridos en dicho sistema.
Si en un sistema de potencia se piensa que hay problemas de oscilaciones de
baja frecuencia se puede implementar el monitoreo del mismo para estudiar su
comportamiento ante disturbios; a partir de dicho comportamiento se pueden
identificar las frecuencias causa del problema para posteriormente proporcionar
en el sistema el amortiguamiento de las mismas.
7.2
Aportaciones
Este trabajo presenta una aplicacin del anlisis de sistemas lineales al clculo
de los modos de oscilacin de un Sistema Elctrico de Potencia. Se obtiene un
modelo linealizado para el SEP y se demuestra su validez para predecir la
respuesta natural del SEP ante pequeas perturbaciones mediante el anlisis
de los modos de oscilacin del sistema.
Se presenta, de forma clara, la relacin entre los valores caractersticos,

eigenvalores, eigenvectores y factores de participacin con las frecuencias
naturales de oscilaciones electromecnicas del SEP, las magnitudes relativas y
fase de dichas oscilaciones.
Se extiende el anlisis de pequeas perturbaciones al anlisis de los modos de

oscilacin electromecnicos a otras variables como el voltaje y potencia reactiva
confirmando que dichos modos se reflejan en la posicin y velocidad angular y
generalmente tambin en la potencia activa y flujos de potencia activa pero no
necesariamente en el voltaje y potencia reactiva.
92
Para el caso de oscilaciones amortiguadas se logra la descripcin matemtica

de la respuesta del SEP en funcin de los eigenvalores y eigenvectores en
donde la parte real e imaginaria de los eigenvalores describen el
amortiguamiento y la frecuencia respectivamente. Se demuestra y se ilustra que
el amortiguamiento afecta en menor medida la frecuencia de oscilacin y
tambin que para que el sistema sea estable se requiere la presencia de
amortiguamiento positivo.
En forma alterna, en este trabajo de tesis se presenta una metodologa para la

obtencin de los modos naturales de oscilacin para un sistema de potencia a
partir de registros de disturbios ocurridos en el sistema con la observacin de
que una falla en particular podr excitar algunos de ellos. Esta metodologa es
particularmente til para el caso donde no se conoce por completo los
parmetros de la red elctrica de forma que no puede emplearse la tcnica del
anlisis modal para el clculo de los modos. La combinacin de ambas
metodologas puede ser de gran utilidad al ser aplicadas en el anlisis de
sistemas de potencia para la identificacin de oscilaciones que pudieran
presentarse en el sistema ocasionando la prdida de la estabilidad del sistema
o bien el disparo de elementos importantes del sistema con la consecuente
prdida de carga y/o generacin.
7.3
Recomendaciones para trabajos futuros
Una vez calculados los modos de oscilacin de un sistema, la informacin

obtenida puede ser de gran utilidad usndola de la manera adecuada.
En un sistema del cual se conocen sus parmetros y se puede modelar, es
posible, mediante tcnicas del anlisis de seales, determinar cuales modos de
oscilacin pueden ser excitados ante diferentes disturbios que pudieran
presentarse en el sistema. Primeramente en un programa de estabilidad
93
transitoria se podra simular la respuesta del sistema ante diferentes disturbios,

a partir de esa respuesta obtener los modos que pueden ser excitados.
Ya que un sistema elctrico es dinmico y continuamente variable en cuanto a
sus condiciones de cargabilidad esto podra significar una gran carga de
trabajo, sin embargo el analista del sistema de potencia puede seleccionar
algunas condiciones de operacin que sean las ms significativas del
comportamiento del sistema o en algunas condiciones especiales por ejemplo la
hora de mxima demanda, el mes de mxima demanda, etc. y de esta forma
obtener ayuda valiosa en el mejoramiento de la seguridad operativa del
sistema, as como en el mejor conocimiento anticipado del comportamiento del
sistema ante algunos disturbios.
Una vez conocidos los modos que pueden ser excitados, se puede identificar
los de mayor participacin y simular el efecto que tendra en el sistema el lograr
amortiguar dichos modos. Posteriormente, y una vez confirmado mediante
simulaciones de estabilidad transitoria, que es til el amortiguar positivamente
dichos modos de oscilacin, se tendra que analizar la forma y lugar
convenientes en donde se podra proveer al sistema, de las condiciones de
amortiguamiento para los modos de oscilacin de mayor inters. Lograr esto
seguramente sera un gran avance en la prevencin de comportamientos no
deseados que pueden ocurrir en el SEP tales como las oscilaciones por baja
frecuencia que pueden llegar a provocar la desconexin de elementos como
lneas de transmisin, generadores o cargas del SEP.
94
Referencias
[1] R. T. Byerly, D. E. Sherman, D. K. McLain, Normal modes and mode shapes
applied to dynamic stability analysis, IEEE Transactions on PAS vol. 94,
March/April 1975.
[2] F. Aboytes, F. Snchez, Oscilaciones electromecnicas en sistemas
elctricos de potencia, IEEE, RVP-88, Acapulco, Guerrero, Agosto 1988.
[3] I. J. Perez-Arriaga, G. C. Verghese, F. C. Schweppe, Selective modal
analysis with applications to electric power systems, Part I: Heuristic
introduction, IEEE Transactions on PAS, vol. 101, Septiembre/1982.
[4] G. C. Verghese, I. J. Perez-Arriaga, F. C. Schweppe, Selective modal
analysis with applications to electric power systems, Part II: The dynamic
stability problem, IEEE Transactions on PAS, vol. 101, Septiembre/1982.
[5] P. Kundur, G. J. Rogers, D. Y. Wong, L. G. Lauby, A comprehensive
computer program package for small signal stability analysis of power systems,
IEEE Transactions on Power Systems, November/1990.
[6] Florencio Aboytes, Estabilidad Transitoria en Sistemas Elctricos de
Potencia, Comisin Federal de Electricidad, Mxico, D.F., Julio 1990.
[7] Florencio Aboytes, Estabilidad Dinmica en Sistemas Elctricos de
Potencia, Comisin Federal de Electricidad, Monterrey, Nuevo Len, Octubre
1991.
95
[8] P. M. Anderson, A. A. Fouad, Power System Control and Stability, Vol. I,

Iowa, U.S.A. 1977.
[9] E. W. Kimbark, Power System Stability Vol. III, Wiley, New York, 1964.
[10] Prabha Kundur, Power System Stability and Control, McGraw-Hill,
California, 1994.
[11] Florencio Aboytes, Anlisis de Sistemas de Potencia Monterrey N. L.,
1990.
[12]
Fernando
Snchez
Tello,
Anlisis
Control
de
Oscilaciones
Electromecnicas en Sistemas Elctricos de Potencia, Tesis, UANL, Monterrey

N.L., 1997.
[13] Thomas Kailath, Linear Systems, Prentice-Hall, Inc., Englewoods Clifs,
N.J., 1980.
[14] Florencio Aboytes, Control de Voltaje en Sistemas Elctricos de Potencia,
Comisin Federal de Electricidad, UANL, Monterrey, N.L. Mxico, Mayo, 1991.
[15] J. F. Hauer, Introduction to the 1989 IEEE/PES Symposium on
eigenanalysis
and
frequency
domain
methods
for
system
dynamic
performance, Portland, Oregon.

[16] Alan V. Oppenheim, Alan S. Willsky Seales y Sistemas, 2 Edicin,
Prentice Hall Hispanoamericana, S.A., Mxico, 1998.
[17] David J. DeFatta, Joseph G. Lucas, William S. Hodgkiss Digital Signal
Processing, Chong Moh Offset Printing Pte. Ltd., Singapore, 1988.
96
[18] David J. DeFatta, Joseph G. Lucas, William S. Hodgkiss Matlab for

windows, Chong Moh Offset Printing Pte. Ltd., Singapore, 1988.
[19] J. F. Hauer, The use of prony analysis to determine modal content and
equivalent models for measured power system response, Portland, Oregon.
[20] Juan C. Ramrez, Florencio Aboytes G. Anlisis de oscilaciones de
Voltajes y Potencias Reactivas en Sistemas Elctricos de Potencia, IEEE,
RVP-2000, Gro. Mxico.
[21] Alan V. Openheim, Ronald W. Schafer, Jhon R. Back, Tratamiento de
Seales en tiempo Discreto, 2 Edicin, Prentice Hall, Madrid, 2000.
[22] F. Aboytes, F. Snchez, R. Ramrez, Simulador Interactivo de Sistemas de
Potencia, SISP-PC Manual Del Usuario, CFE, CENACE, Mayo 1993.
[23] Juan C. Ramrez, Alfonso Lpez, Luis D. Archiga, Determinacin de
Modos de Oscilacin a partir de registros de disturbios en Sistemas Elctricos
de Potencia, CFE, CENACE, IEEE, RVP-2009, Gro. Mxico.
[24] I. J. Prez - Arriaga, G. C. Verghese, F. Luis Pagola, F. C. Schweppe,
Selective modal analysis in power systems, American Control Conference,
San Francisco, EE.UU. June/1983.
[25] J. J. Sanchez-Gasca, J. H. Chow, Performance comparison of three
identification methods for the analysis of electromechanical oscillations, IEEE
Transactions on power systems, August/1999.
[26] D. Y. Wong, G. J. Rogers, B. Porreta, P. Kundur, Eygenvalue analysis of
very large power systems, IEEE Transactions on power systems, May/1988.
97
[27] J. F. Hauer, C. J. Demeure, L. L. Sharf, Initial results in prony analysis of

power systems response signals IEEE Transactions on Power Systems, vol. 5,
February/1990.
[28] Jess H. Hernndez L, A. R. Messina, Utilizacin de dispositivos de control
en la red de transmisin para amortiguar oscilaciones de potencia, Comit
Nacional de CIGRE, Junio/2001.
[29] J. G Caldern-Guizar, Anlisis modal de oscilaciones electromecnicas en
sistemas elctricos de potencia, IIE, Mxico, Agosto/2007.
[30] J. L Rueda, D. G. Colom, Mtodos para el anlisis de estabilidad de
pequea seal en sistemas elctricos de potencia, XII ERIAC del CIGR,
Brasil, Mayo/2007.
[31] C. Sharma, P. Shing. The effect of load types on power systems
oscilations, IEEE/PES transmission and distribution conference and exposition,
Chicago, Il. 2008.
[32] R. Castellanos B., A. R. Messina, H. Sarmiento U., G. Pampin.
Assessment of remedial control schemes for damping transient oscilations in
the Mexican system, Conference on Power Systems, Winnipeg, Canada,
CIGR, October/2008.
[33] C. A. Jurez, D. G. Colom, Un nuevo enfoque en la Supervisin y Control
de Oscilaciones Electromecnicas del Sistema de Potencia, IEEE, El Salvador,
Noviembre/2011.
[34] J. Turunen, J. Thambirajah, M. Larsson, B. C. Pal, N. F. Thornhill, L. C.
Haarla, W. W. Hung, A. M. Carter, T. Rauhala. Comparison of three
98
electromechanical oscillation damping estimation methods, IEEE Transactions

on power systems, November/2011.
[35] M. A. Andrade, A. R. Messina, C. A. Rivera, D. Olguin. Identification of
instantaneous attributes of torsional shaft signals using Hilbert transform, IEEE
Transactions on power systems, August/2014.
[36] A. R. Messina, M. A. Andrade, J. H. Hernndez, R. Betancourt. Analysis
and characterizations of power system nonlinear oscillations using Hilbert
spectral analysis, The Open Electrical and Electronic Engineering Journal,
June/2007.
99
ndice de tablas
Tabla
3-1
Pgina
Anlisis modal del sistema mquina-bus infinito en condiciones
de generacin mnima. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3-2
25
Anlisis modal del sistema mquina-bus infinito en condiciones

de generacin media. . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . .
26
3-3
Anlisis modal de un sistema de prueba de dos mquinas. . . .
28
3-4
Anlisis modal de un sistema de prueba multimquina . . . . . . .
31
4-1
Anlisis modal de un sistema de dos mquinas con

amortiguamiento. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4-2
Informacin obtenida en el anlisis modal para un sistema

multimquina. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6-1
49
51
Espectro de frecuencias de la respuesta oscilatoria de los

generadores G1, G2, G3, G4 y G5 ante la ocurrencia de un
disturbio. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6-2
80
Espectro de frecuencias de la respuesta oscilatoria de los

generadores G1, G2, G3, G4 y G5 ante la ocurrencia de un
disturbio. Las frecuencias se han ordenado para indicar en
cuales generadores est presente cada frecuencia . . . . . . . . . .
82
6-3
Clculo de modos de oscilacin del ejemplo 3.3 . . . . . . . . . . . .
82
6-4
Factores de participacin para la oscilacin de la velocidad

angular de mquinas. Resultados del anlisis lineal. . . . . . . . . .
6-5
83
Factores de participacin para la oscilacin de la velocidad

angular de mquinas. Resultados del anlisis lineal. Cambia el
arreglo de la informacin de la tabla 6-4. . . . . . . . . . . . . . . . . . .
100
84
ndice de figuras
Figura
Pgina
3.1
Sistema de prueba de Mquina-Bus infinito. . . . . . . . . . . . . . .
3.2
A) Oscilacin de la frecuencia. B) Oscilacin de la posicin

angular. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.3
A) Oscilacin de la velocidad angular. B) Oscilacin de la

posicin angular. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.4
Sistema de prueba de dos mquinas. . . . . . . . . . . . . . . . . . . . .
3.5
Excitacin del modo de oscilacin en el sistema de prueba de

dos mquinas A) Posicin angular. B) Potencia Elctrica . . . . .
3.6
25
26
27
28
29
Oscilacin de la A) posicin angular y B) potencia elctrica al

ocurrir una falla. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
30
3.7
Sistema de prueba para el caso multimquina. . . . . . . . . . . . . .
30
3.8
Oscilacin de la A) posicin angular y B) potencia elctrica

para el modo1 (f=2.28Hz). . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
32
3.9
A, B y C. Modos de oscilacin en un sistema multimquina . . .
32
4.1
37
4.2
Regin lineal de las funciones coseno y seno . . . . . . . . . . . . . .
40
4.3
Efecto de los coeficientes de sincronizacin K1 y

amortiguamiento D en la estabilidad. . . . . . . . . . . . . . . . . . . . .
41
4.4
Sistema de prueba de dos mquinas con amortiguamiento. . . .
48
4.5
A) Oscilacin de la velocidad y B) Oscilacin de la posicin

angular. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
49
4.6
Sistema de Prueba Multimquina. . . . . . . . . . . . . . . . . . . . . . .
50
4.7
Excitacin del modo 1 de 2.2778Hz. . . . . . . . . . . . . . . . . . . . . .
52
4.8
Excitacin de modo 5, f=1.05Hz. . . . . . . . . . . . . . . . . . . . . . . . .
53
101
ndice de figuras (continuacin)

5.1
56
5.2
Regin lineal de las funciones coseno y seno . . . . . . . . . . . . . .
61
5.3
Sistema de prueba de dos mquinas en condiciones de flujo

mximo de potencia. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.4
62
Comportamientos de variables del sistema de dos mquinas

ante una perturbacin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
63
5.5
Sistema de prueba operando en la regin no-lineal. . . . . . . . . .
64
5.6
A) Posicin angular, B) Diferencia angular, C) y D)

Comportamiento no lineal del Voltaje y la Potencia Reactiva
respectivamente. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.7
Efecto en el cambio en la condicin de operacin en el

Sistema de Potencia de prueba . . . . . . . . . . . . . . . . . . . . . . . .
5.8
64
65
Efecto de la condicin de operacin en las variables del

sistema. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
66
5.9
Sistema de Potencia de prueba de tres generadores . . . . . . .
66
5.10
(A y B). Diferencia y posicin angular de mquinas. . . . . . . . .
67
5.11
A) Voltaje y B) Potencia reactiva de generadores. . . . . . . . . . .
68
5.12
68
5.13
69
5.14
A) Diferencia y B) Posicin angular de mquinas. . . . . . . . . . .
69
5.15
70
5.16
71
5.17
71
6.1
Sistema de Prueba Multimquina. . . . . . . . . . . . . . . . . . . . . . .
77
6.2
Oscilacin de velocidades de mquinas al ocurrir un disturbio.
78
6.3
Espectro de Frecuencias de las Oscilaciones de mquinas . . .
79
6.4
Amplitud y Fase de componentes de la mquina 1. . . . . . . . . .
80
102
APNDICE 1
Sistemas de prueba utilizados para los casos de
estudio
103
SISTEMA MAQUINA-BUS INFINITO:
Datos de Nodos
1Nodo
2_Pl
3_Ql
4_Tc
5_Tn
6SBase
7KVline
Comentario:
0.0
0.0
100
20
Nodo de BT G1
0.0
0.0
100
20
Nodo de BT G2
0.0
0.0
100
230
Nodo de carga
0.0
0.0
100
230
Nodo de carga
0.0
0.0
100
230
Nodo de carga
Datos de Generadores:
1Nodo
Comentario:
2PE_Gen 3Vt_Esp
5.0
Bus Infinito Gen_1
1.0
1.02
Generador 2
Datos de Lneas y Transformadores:

1_Ns
2_Nll
3_r
4_x
5_y/2
6_R.T.
7_ID
8SBase
9Vline
Comentario:
0.04
100
230
Trasf. 1
0.04
100
230
Trasf. 2
0.01206 0.08242 0.08
100
230
L.T.
0.02412 0.16484 0.164
100
230
L.T.
0.01206 0.08242 0.08
100
230
L.T.
104
SISTEMA DE DOS MAQUINAS FINITAS CON PC=50MW:

Datos de Nodos
1Nodo
2_Pl
3_Ql
4_Tc
5_Tn
6SBase
7KVline
Comentario:
0.0
0.0
100
20
Nodo BT G1
0.0
0.0
100
20
Nodo BT G2
0.0
0.0
100
230
Nodo de carga
0.5
0.0
-1
100
230
Nodo de carga
1Nodo
Comentario:
2PE_Gen 3Vt_Esp
5.0
Gen_1
1.0
1.02
Gen_2

1_Ns
2_Nll
3_r
4_x
5_y/2
6_R.T.
7_ID
8SBase
9Vline
Comentario:
0.040
100
230
Trasf. 1
0.040
100
230
Trasf. 2
0.100
100
230
L.T.
105
SISTEMA MULTIMAQUINA DE 5 GENERADORES.
Datos de Nodos
1Nodo
2_Pl
3_Ql
4_Tc
5_Tn
6SBase
0.0
0.0
100
20
Nodo BT G1
0.0
0.0
100
20
Nodo BT G2
0.0
0.0
100
20
Nodo BT G3
0.0
0.0
100
20
Nodo BT G4
0.0
0.0
100
20
Nodo BT G5
0.0
0.0
100
230
Nodo de carga
0.0
0.0
100
230
Nodo de carga
0.0
0.0
100
230
Nodo de carga
0.0
0.0
100
230
Nodo de carga
10
0.0
0.0
100
230
Nodo de carga
7KVline Comentario:
Comentario:
1Nodo 2PE_Gen 3Vt_Esp

1
Compensador.
0.05
1.02
Gen_2
0.05
1.02
Gen_3
0.05
1.02
Gen_4
0.05
1.02
Gen_5

1_Ns
2_Nll
3_r
4_x
5_yn2
6_R.T.
7_ID
8SBase
9Vline
Comentario:
0.033
100
230
Trasf. 1
0.033
100
230
Trasf. 2
0.033
100
230
Trasf. 3
0.033
100
230
Trasf. 4
10
0.033
100
230
Trasf. 5
0.0146
0.0973 0.092
100
230
L.T.
0.0434
0.2891 0.275
100
230
L.T.
0.0293
0.1987 0.191
100
230
L.T.
0.0035
0.0233 0.025
100
230
L.T.
106
SISTEMA DE 3 GENERADORES.
Datos de Nodos
1Nodo
2_Pl
3_Ql
4_Tc
5_Tn
6SBase
7KVline
Comentario:
0.0
0.0
100
20
Nodo BT G1
0.0
0.0
100
20
Nodo BT G2
0.0
0.0
100
20
Nodo BT G3
0.0
0.0
100
230
Nodo de carga
0.0
0.0
-1
100
230
Nodo de carga
0.0
0.0
-1
100
230
Nodo de carga
Comentario:
1Nodo 2PE_Gen 3Vt_Esp

1
3.0
1.02
Compensador.
0.0
1.02
Gen_2
0.0
1.02
Gen_3
1_Ns
2_Nll
3_r
4_x
5_yn2
6_R.T.
7_ID
8SBase
9Vline
Comentario:
0.010
100
230
Trasf. 1
0.010
100
230
Trasf. 2
0.010
100
230
Trasf. 3
0.260
100
230
L.T.
0.060
100
230
L.T.
0.260
100
230
L.T.
107
APNDICE 2
Artculos presentados referentes al trabajo de
tesis
108
ANLISIS DE OSCILACIONES DE VOLTAJES Y POTENCIAS REACTIVAS EN

SISTEMAS ELCTRICOS DE POTENCIA
F. Aboytes *
J.C. Ramrez*
Unidad de Ingeniera Avanzada

Centro Nacional de Control de Energa
Comisin Federal de Electricidad
faboytes@cfe.gob.mx
*Programa Doctoral en Ingeniera Elctrica

Facultad de Ingeniera Mecnica y Elctrica
Universidad Autnoma De Nuevo Len
jramirez@die.fime.uanl.mx
RESUMEN
Se presenta el anlisis del comportamiento de voltajes
terminales y potencias reactivas de generadores al estimular
una oscilacin electromecnica en un sistema de potencia
bsico. Se detalla el comportamiento no-lineal de estas
variables y se ilustran casos crticos en sistemas de prueba.
En el trabajo se presentan conclusiones referentes al efecto
de la posicin del rotor del generador en la dinmica de los
voltajes en terminales y la potencia reactiva de generadores
para los modos de oscilacin lineales.
INTRODUCCIN
El anlisis de la estabilidad en un sistema elctrico es una
tarea importante en la planificacin y la operacin de un
sistema. Los sistemas de potencia estn sujetos de manera
permanente
a
diversas
perturbaciones
debido
principalmente a la variacin de la carga, a operaciones
normales, cambios en generacin, conexin y desconexin
de elementos y tambin debido a fallas en generadores y
lneas de transmisin. El sistema para que sea estable debe
producir simultneamente fuerzas de sincronizacin que
permitan mantener las mquinas unidas operando en
paralelo y fuerzas de amortiguamiento que lleven a una
operacin en estado estable despus de cualquier
perturbacin.
Generalmente cuando se perturba un sistema de potencia se
produce la excitacin de algunos modos de oscilacin
electromecnicos que se observan en la dinmica de la
posicin angular y en la potencia activa de las mquinas
que participan en esos modos, tambin se detectan en la
dinmica de algunos flujos de potencia activa en el sistema
de transmisin.
Las tcnicas de anlisis lineal se han utilizado

tradicionalmente para caracterizar las oscilaciones
electromecnicas en los sistemas de potencia, identificando
los generadores que participan en cada modo de oscilacin y
tambin su nivel de participacin, este anlisis es vlido
estrictamente ante pequeas perturbaciones donde las
variables de inters tienen movimientos reducidos.
El objetivo de este trabajo es analizar el impacto de estas
pequeas perturbaciones en otras variables que definen el
comportamiento reactivo de los generadores y los voltajes
del sistema. En la primera parte del artculo se desarrollan
las ecuaciones para el voltaje en terminales y la potencia
reactiva generada, se obtienen expresiones en funcin de la
posicin angular de rotores para un sistema bsico, usando
este desarrollo se definen las regiones lineales y no-lineales
del comportamiento de las unidades en funcin del punto de
operacin.
Mediante el anlisis de casos ilustrativos se muestra la
deformacin de voltajes y potencias reactivas resultado de la
condicin operativa y del movimiento de las posiciones
angulares de rotores de unidades generadoras. El anlisis se
extiende a diversos modos de oscilacin y se confirma el
mismo efecto de no-linealidad para un sistema con dos
modos de oscilacin. Se incluyen conclusiones referentes al
anlisis y los desarrollos presentados en el trabajo.
ECUACIONES DE VOLTAJE
Para desarrollar las ecuaciones de los voltajes terminales se
utiliza el sistema de prueba mostrado en la figura 1.
X3
V1
X1
X2
E11
V2
E22
Figura 1. Sistema de prueba de dos mquinas

La corriente en el sistema est dada por:
1
E1 E 2
I
j( X 1 X 2 X 3 )
109
(1)
que depende de una funcin coseno de la diferencia de

ngulos.
aplicando la ley de Kirchhoff, se obtienen los voltajes en

terminales de las mquinas, V1 y V2:
1
ECUACIONES DE POTENCIA REACTIVA
V1 E1 jX 1 I
(2A, 2B)
En el sistema de la figura 1, la potencia de generacin para

la mquina 1 est dada por:
V2 E1 j X 1 X 3 I
1
substituyendo la ecuacin 1, y arreglando las ecuaciones

resulta:
1
E ( X 2 X 3 ) E2 ( X 1 )
V1 1
X1 X 2 X 3
1
V2
(6)
o bien
*
*
S1 V1 I P1 jQ1
(7)
Substituyendo en (7) la ecuacin (1) para la corriente y el

conjugado de la ecuacin (3) para el voltaje en terminales
de la mquina 1, se obtiene:
(3A, 3B)
E1 ( X 2 ) E 2 ( X 1 X 3 )
X1 X 2 X 3
y si en (3A) y (3B) se aplica la identidad de Euler, los

voltajes pueden expresarse como:
X2 X3
X1
V1
E2 cos( 2 )
E1 cos(1 )
X
X
X
X
X
X
2
3
2
3
1
1
X 2 X 3
X1
j
E2 sen( 2 )
E1sen(1 )
X
X
X
X
X
X
2
3
2
3
1
1
(4A)
X1 X 3
X2
V2
E1 cos( 1 )
E2 cos( 2 )
X
X
X
2
3
1
X1 X 2 X 3
X1 X 3
X2
j
E2 sen( 2 )
E1sen( 1 )
X
X
X
X
X
X
2
3
2
3
1
(4B)
E 1 ( X 2 X 3 ) E 2 2 ( X 1 ) E1 1 E 2 2
P1 jQ1 1
(X1 X 2 X 3 )
j ( X 1 X 2 X 3 )
(8)
donde al realizar las operaciones algebraicas pertinentes se
obtienen las ecuaciones para la potencia de generacin,
activa y reactiva de la mquina 1:
P1
E1 E 2
sen( 12 )
X1 X 2 X 3
X 2 X 3 2
2
X1
E
E
Q1
2 1
2 2
X 1 X 2 X 3
X 1 X 2 X 3
X X 2 X 3
E E cos( 12 )
1
2 1 2
X 1 X 2 X 3
donde 12 1 2
En forma alterna, partiendo de (3), la magnitud de los

voltajes tambin puede expresarse de la manera siguiente:
2
S1 V1 I * P1 jQ1
de manera similar para

quedan:
(10)
(11)
la mquina 2, las ecuaciones
X2 X3 2
X1
2
V1
E1
E2
X1 X 2 X 3
X1 X 2 X 3
2 X1X 2 X 3
E E cos( 1 2 )
X 1 X 2 X 3 2 1 2
2
(9)
P2
(5A)
X1 X 3 2
2
X2
Q2
E
E
2 1
2 2
X 1 X 2 X 3
X 1 X 2 X 3
X X1 X 3
2
E E cos( 12 )
2 1 2
X 1 X 2 X 3
2 X1 X 3
X2
2
2
V2
E2
E1
X
X
X
X
X
X
2
3
2
3
1
1
2 X 2 X1 X 3
E E cos( 1 2 )
X 1 X 2 X 3 2 1 2
E1 E 2
sen( 12 )
X1 X 2 X 3
(5B)
(12)
(13)
Las potencias activas, ecuaciones (9) y (12), son una

funcin senoidal de la diferencia de los ngulos de mquinas
y tienen una amplitud determinada. Las potencias reactivas
estn dadas por un trmino constante sumado de otro
trmino que depende de una funcin coseno de la diferencia
de los ngulos.
En la ecuacin (4) el voltaje se expresa en su parte real e

imaginaria, con funciones coseno y seno de los ngulos de
las mquinas, respectivamente; en cambio en la ecuacin
(5) el voltaje est dado por un trmino constante, ms otro
110
ANLISIS DE OSCILACIONES
ELECTROMECNICAS
f(t)=cos(t)
Se sabe que en un sistema de potencia de n generadores

existen (n-1) modos de oscilacin electromecnicos, los
cuales pueden ser calculados usando una representacin
clsica de mquinas y algunas consideraciones [1,2]. Los
modos de oscilacin estn directamente relacionados con
las oscilaciones de variables como velocidades angulares
(), ngulos de rotores () y potencias activas, sin embargo
estas oscilaciones producen oscilaciones en todas las
variables del sistema de potencia, entre ellas el voltaje y la
potencia reactiva. El objetivo de esta seccin es el analizar
este comportamiento.
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
-0.2
-0.2
-0.4
-0.6
-0.8
-0.8
-1
-200
REGION
LINEAL
-0.4
REGION
LINEAL
-0.6
-1
-100
100
200
-200
-100
0
t(grados)
100
200
t(grados)
Figura 2. Funciones coseno y seno
misma forma de oscilacin, en cambio la amplitud puede

variar ya que depende de la combinacin, en forma diferente
para cada mquina, de los voltajes internos y reactancias de
la red, (ecuacin (5)).
En el anlisis lineal de modos de oscilacin se considera

una variacin senoidal de la posicin angular de las
mquinas a una velocidad angular de oscilacin w, es decir:
A sen(wt )
f(t)=sen(t)
A su vez la potencia reactiva de generacin puede ser

positiva o negativa en la condicin inicial ya que el trmino
constante en las ecuaciones (10) y (13) se compone de una
suma algebraica de otros dos miembros, uno positivo y uno
negativo. Adems, depende de la misma diferencia angular
que el voltaje en terminales, solo que el trmino que da la
amplitud de oscilacin generalmente es negativo, lo que
har que la oscilacin de la potencia reactiva sea opuesta a
la del voltaje. La conclusin del anlisis realizado es que los
voltajes y las potencias reactivas generadas no oscilan de la
misma forma que los ngulos de las mquinas.
(14)
de esta forma si cada mquina que participa en la oscilacin

tiene una variacin angular como en (14), la potencia activa
ser una funcin no lineal ya que en las ecuaciones (9) y
(12) se observa que depende del seno de la diferencia de
dos senoidales, sin embargo si el movimiento de los
ngulos se realiza en la regin ms lineal de la funcin seno
(figura 2), es decir alrededor de cero grados, la oscilacin
de la potencia activa sigue un comportamiento similar a la
posicin angular. De esta forma el modo de oscilacin es
observable tambin en la potencia activa generada.
A continuacin se presentan algunos casos de prueba para

ilustrar el comportamiento de las variables del sistema ante
la aplicacin de una perturbacin.
En el caso del voltaje en terminales de generadores, estos se

pueden expresar como la suma de un trmino constante que
depende de las magnitudes de los voltajes internos de
mquinas y las reactancias de la red y un trmino que
depende del coseno de la diferencia de los ngulos, los
voltajes tendrn en general un comportamiento no-lineal
ante la variacin de los ngulos de rotores.
Caso 1. En este caso el sistema se encuentra descargado al

momento de la perturbacin. El objetivo es operar el sistema
en la regin ms no-lineal para observar el comportamiento
de voltajes y potencias reactivas.
1
Para explicar este comportamiento, en la figura 2 se marcan

las regiones lineales de las funciones coseno y seno, la
regin lineal en la funcin coseno se tiene para valores
entre 30 y 150 grados, es decir para diferencias angulares
grandes. En cambio para valores cercanos a cero grados el
comportamiento es altamente no-lineal. Como el
comportamiento del voltaje y la potencia reactiva dependen
del coseno de una diferencia angular, la respuesta en el
tiempo seguir la oscilacin de ngulos solo si se opera en
la regin lineal, en el caso contrario se producir una
distorsin causada por la no linealidad de la caracterstica.
2
G2
G1
0
Figura 3. Sistema de prueba para el caso 1
En la figuras 4A y 4B se observa el comportamiento de la

diferencia angular y los ngulos de rotores en el momento
en que se perturba el sistema, las mquinas oscilan una
contra la otra, la diferencia angular est alrededor de 3
grados.
Otro aspecto relevante en el anlisis es la forma de la

oscilacin de los voltajes en terminales para ambas
mquinas, las ecuaciones (5A) y (5B) muestran que
dependen del coseno de la misma diferencia angular, lo
cual indica que ambos voltajes en terminales tendrn la
Como resultado de esta no linealidad, en la figura 5A se

observa una doble frecuencia en la oscilacin del voltaje de
la mquina 1. La mquina 2 tiene un comportamiento
similar. En la figura 5B se muestra el comportamiento de la
111
potencia reactiva que tiene un forma de oscilacin similar

pero opuesta al voltaje.
Se observa una mayor distorsin para el voltaje en

terminales de la mquina 2 (figura 8A), ya que la posicin
angular de esta mquina est oscilando en la regin nolineal.
Diferencia angular
4
Difd (grad)
2
0
En la figura 8B, se observa que en la oscilacin de la

potencia reactiva tambin existe distorsin, aunque, de
acuerdo a las ecuaciones (10) y (13), el efecto de la no
linealidad es menor que en el voltaje.
-2
-4
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8

2
d(grad)
1
G1
G2
-1
1.02
1.015
-2
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
1.01
Vt(pu)
Figura 4(A y B). Diferencia y posicin angular de mquinas
1.005
1
0.995
0.99
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8

1.0005
G1
G2
30
Q(Mvar)
Vt(pu)
20
10
-10
0.9995
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 8. A)Voltaje y B)Potencia reactiva de generadores
-10.5
Q(Mvar)
0.2
G1
G2
-11
Caso 3. En este caso se pretende operar al sistema en la

regin lineal para determinar el impacto en la oscilacin de
voltajes y potencia reactiva.
-11.5
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 5. A) Voltaje y B)Potencia reactiva de generadores
Caso 2. En este caso los generadores alimentan a una carga,

el objetivo es analizar el cambio en la condicin de
operacin y observar el comportamiento resultante.
1
100
50
G2
En la figura 10A se observa la diferencia angular oscilando

en valores cercanos a los 15 grados, en (11A) y (11B) se
muestra la oscilacin de voltajes y potencias reactivas, se
observa un poco de distorsin, pero oscilan a la misma
frecuencia del modo de oscilacin del ngulo. La
explicacin es que se opera cerca de la regin lineal.
100
200
2
G2
50
G1
100
G1
En este caso el sistema est en un punto de operacin donde

la diferencia angular se mueve dentro de la regin no-lineal
y est variando entre -14 y 3 grados, figura 7.
Diferencia angular
Diferencia angular
5
20
15
Difd (grad)
Difd (grad)
-5
10
5
-10
-15
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.4
1.6
1.8

20
15
15
10
d(grad)
d(grad)
10
5
0
-5
G1
G2
5
0
G1
G2
-5
-10
-10
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.8
Figura 10 (A y B). Diferencia y posicin angular de

mquinas
Figura 7 (A y B). Diferencia y posicin angular de mquinas
112
1.03
1.02
1.0199
1.0198
1.01
Vt(pu)
Vt(pu)
1.02
1.0197
1.0196
1.0195
0.99
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.0194
0
G1
G2
0.2
0.4
0.6
0.8
30
1.2
1.4
1.6
G1
G2
G3
1.8

0.4
20
0.2
Q(Mvar)
Q(Mvar)
0.3
10
0
0.1
0
-10
-0.1
-20
0
0.5
1.5
t(seg)
2.5
-0.2
Figura 11. A)Voltaje y B)Potencia reactiva de

generadores
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Para el segundo modo de oscilacin (1.63 hz) se confirma

tambin que la frecuencia de oscilacin del voltaje y la
potencia reactiva tambin es el doble de la oscilacin de los
ngulos, figuras 15 y 16.
En el sistema de tres mquinas de la figura 12 se

estimularn los modos de oscilacin por separado para
observar el comportamiento de voltajes y potencias
reactivas en cada modo.
Modo 2 f= 1.6314Hz
Diferencia angular
10
Difd (grad)
2
G2
G1
-5
-10
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8

G1
G2
G3
d(grad)
3
G3
Figura 12. Sistema de prueba de tres mquinas
-5
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Caso 4. En este caso el sistema est descargado, el objetivo

es confirmar el comportamiento de voltajes y potencias
reactivas para esta regin no-lineal.

1.02
Vt(pu)
1.0195
1.019
1.0185
0
0.2
0.4
0.6
Modo 1 f= 2.5214 Hz
0.8
1.2
1.4
1.6
1.8
G1
G2
G3

1
Diferencia angular
0.8
Q(Mvar)
Difd (grad)
0.6
0.4
0.2
-2
0
0
-4
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
G1
G2
G3
d(grad)
Caso 5. En este caso el sistema de tres mquinas suministra

la potencia demandada en dos nodos de carga. El objetivo es
operar el sistema en una regin lineal.
-2
-4
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 13(A y B). Diferencia y posicin angular de mquinas
En las figuras 17 y 19 se observa la respuesta en el tiempo

para los ngulos de rotores ante el disturbio; los voltajes y
potencias reactivas oscilan a la frecuencia de la posicin
angular (figuras 18 y 20), sin embargo hay un grado de
distorsin dada por la no linealidad del coseno.
En la figura 13A y 13B se observa la oscilacin de los

ngulos para el primer modo de oscilacin (2.52 hz), el
efecto no-lineal se ve reflejado en la figura 14A y 14B, se
observa una doble frecuencia en la oscilacin de voltajes y
potencias reactivas.
113
CONCLUSIONES
Modo 1 f= 2.5240 Hz
Diferencia angular
0
Comportamiento de voltajes en terminales
Difd (grad)
-5
-10
La condicin inicial para los ngulos de los voltajes internos

de mquinas depende del punto de operacin del sistema, al
ocurrir una perturbacin se tendrn oscilaciones en estas
variables del sistema. Si la perturbacin es pequea los
ngulos de los voltajes internos de mquinas siguen una
variacin senoidal alrededor de su condicin inicial y
oscilan en forma opuesta para el caso de dos mquinas. Los
voltajes en terminales dependen directamente de la
variacin de estos ngulos y por lo tanto la forma de
oscilacin es tambin funcin del punto de operacin.
-15
-20
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
d(grad)
G1
G2
G3
-5
-10
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Para ambos modos de oscilacin se cumplen estas

observaciones, adems que la oscilacin de cada potencia
reactiva generada es opuesta al voltaje en terminales de esa
mquina.
En general el comportamiento dinmico de los voltajes es de

naturaleza no-lineal debido a que la oscilacin es descrita
por una suma de senos y cosenos de funciones senoidales.
Adems, los voltajes de los generadores oscilan en la misma
direccin ya dependen principalmente de funciones coseno,
en las cuales, sin importar si crece o decrece la diferencia
angular, el resultado siempre ser positivo.

1.022
Vt(pu)
1.021
1.02
1.019
La distorsin de los voltajes en terminales de mquinas

depender del ngulo del voltaje interno de las mquinas, si
se opera en la regin no-lineal la distorsin produce otras
frecuencias de oscilacin en los voltajes, si se opera en la
regin lineal la distorsin ser menor. La conclusin es que
los voltajes no siguen el comportamiento del modo de
oscilacin de las posiciones angulares.
1.018
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
1.4
1.6
1.8

8
7
G1
G2
G3
Q(Mvar)
6
5
4
3
2
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
Comportamiento de la potencia reactiva
Modo 2 f= 1.6146 Hz
Diferencia angular
La oscilacin de la potencia reactiva generada es en general

no-lineal ya que est dada por un valor constante sumado o
restado de una funcin coseno de la diferencia de dos
funciones senoidales. Sin embargo, esta oscilacin puede
tener un comportamiento bastante lineal, lo que depende
directamente del comportamiento de la diferencia angular.
Difd (grad)
-5
-10
-15
-20
-25
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
G1
G2
G3

10
d(grad)
El tener la posibilidad de que el trmino constante pueda

ser positivo o negativo indica que la mquina en la
condicin inicial puede estar entregando o recibiendo
potencia reactiva, el trmino que produce la oscilacin
puede a su vez ser positivo o negativo, lo cual describe la
forma de oscilacin. Se debe resaltar que no es importante si
la diferencia angular es positiva o negativa ya que el coseno
es una funcin par, lo que si es relevante en la oscilacin es
que si existe cambio de signo en la diferencia angular, que
se tiene cuando el sistema opera en una regin no-lineal, se
tiene distorsin en el comportamiento de la potencia
reactiva. Si bien, tanto la potencia reactiva como el voltaje
en terminales dependen del coseno de la diferencia angular,
la oscilacin de la potencia reactiva es menos sensible
debido a la amplitud de ste trmino (ecuaciones (5), (10) y
(13)).
0
-5
-10
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8

1.03
Vt(pu)
1.025
1.02
G1
G2
G3
1.015
1.01
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8

10
Q(Mvar)
-5
0
0.2
0.4
0.6
0.8
1
t(seg)
1.2
1.4
1.6
1.8
Figura 20.A)Voltaje y B)Potencia reactiva de generadores
114
Science and Technology de Londres, Inglaterra en 1974.

Desde 1975 trabaja en Comisin Federal de Electricidad
donde ha ocupado diversos puestos tcnicos en la Gerencia de
Estudios y la Gerencia del Centro Nacional de Control de
Energa (CENACE). Actualmente es Jefe de la Unidad de
Ingeniera Avanzada del CENACE.
BIBLIOGRAFIA
[1] R. T. Byerly, D. E. Sherman, D. K. McLain, Normal
modes and mode shapes applied to dynamic stability
analysis, IEEE Transactions on PAS vol. 94, March/April
1975.
Adicionalmente, desde 1975 ha impartido cursos de posgrado

y dirigido trabajos de investigacin en la Seccin de
Graduados del ESIME-IPN, la Divisin de Estudios de
Posgrado de la UNAM, y en el Programa de Doctorado de la
Facultad de Ingeniera Mecnica y Elctrica de la
Universidad Autnoma de Nuevo Len (FIME-UANL).
Actualmente es Presidente del Comit Doctoral de FIMEUANL.
[2] F. Aboytes, F. Sanchez, Oscilaciones electromecnicas

en sistemas elctricos de potencia, RVP-88, Acapulco,
Guerrero, Agosto 1988.
[3] I. J. Perez-Arriaga, G. C. Verghese, F. C. Schweppe,
Selective modal analysis with applications to electric
power systems, Part I: Heuristic introduction, IEEE
Transactions on PAS, vol. 101, Septiembre/1982.
Es Fellow Member del IEEE, Acadmico de Nmero de la

Academia Mexicana de Ingeniera, miembro de CIGRE,
Presidente del Consejo Mxico del IEEE, y es Investigador
Nacional del Sistema Nacional de Investigadores.
[4] G. C. Verghese, I. J. Perez-Arriaga, F. C. Schweppe,

Selective modal analysis with applications to electric
power systems, Part II: The dynamic stability problem,
IEEE Transactions on PAS, vol. 101, Septiembre/1982.
[5] G. J. Roger, P. Kundur, Small signal stability of power
systems, Ontario, Canada.
Juan C. Ramrez, naci en Zacatecas en 1974, donde

realizo sus estudios primarios. Desde 1987 vive en
Monterrey y en 1997 se gradu de Ingeniero Mecnico
Electricista en la Facultad de Ingeniera Mecnica y
Elctrica de la Universidad Autnoma de Nuevo Len
donde actualmente realiza estudios de maestra en el
Programa Doctoral en Ingeniera Elctrica. Ha trabajado
para la iniciativa privada en el rea de proyectos y asesora
en la rama electromecnica y a la fecha es catedrtico de
tiempo parcial en la Universidad del Norte en la divisin de
Ingeniera Elctrica. Su rea de inters se centra en el
anlisis de oscilaciones electromecnicas en Sistemas
Elctricos de Potencia.
[6] Florencio Aboytes, Estabilidad Transitoria en Sistemas

Elctricos de Potencia, Comisin Federal de Electricidad,
Mxico, D.F., Julio 1990.
[7] Florencio Aboytes, Estabilidad Dinmica en Sistemas
Elctricos de Potencia, Comisin Federal de Electricidad,
Monterrey, Nuevo Len, Octubre 1991.
Dr. Florencio Aboytes Garca. Ingeniero Mecnico
Electricista y Maestro en Ciencias del Instituto Tecnolgico y
de Estudios Superiores de Monterrey en 1968 y 1969.
Doctorado en Ingeniera Elctrica en el Imperial College of
115
DETERMINACIN DE MODOS DE OSCILACIN A PARTIR DE REGISTROS DE

DISTURBIOS EN SISTEMAS ELCTRICOS DE POTENCIA
Juan Carlos Ramrez Saucedo
CFE - rea de Control Noreste
juan.ramirez01@cfe.gob.mx
Alfonso Lpez Delgado

alfonso.lopez@cfe.gob.mx
Luis David Archiga Rodrguez

luis.arechiga@cfe.gob.mx
En el anlisis de SEPs es comn el uso de esquemas de

monitoreo y medicin de parmetros elctricos
(registradores de eventos que incluyen Unidades de
Medicin Fasorial, PMUs por su nombre en ingles) que
faciliten el estudio del comportamiento del sistema ante
disturbios. A partir de registros apropiados y
considerando los datos registrados como muestras de una
seal discreta, es posible analizar la seal de respuesta en
el tiempo de modo de obtener informacin de la
magnitud y fase de las frecuencias que componen a dicha
seal.
Resumen
Se presenta una alternativa para el clculo de los modos
de oscilacin naturales de un sistema elctrico de
potencia (SEP) para cuando no se conocen por completo
los parmetros fsicos de la red, de manera que no es
sencillo establecer, con esta base, el estado del sistema en
alguna condicin inicial de operacin y por tanto no
pueden utilizarse las tcnicas convencionales del anlisis
modal.
La alternativa propuesta requiere tener registros de datos
del comportamiento del SEP, para un intervalo de tiempo
determinado, durante perturbaciones que se hayan
presentado en dicho sistema. Mediante el uso de tcnicas
apropiadas se analiza la informacin contenida en los
registros para obtener informacin de los modos de
oscilacin del sistema en cuestin.
Desde el punto de vista de la estabilidad del Sistema de

Potencia, es conveniente identificar las oscilaciones de
baja frecuencia que se pueden presentar ya que podran
causar desde una pequea variacin en los parmetros
elctricos hasta prdida de algunos elementos del SEP.
Este trabajo presenta ejemplos ilustrativos de la
obtencin e interpretacin de los modos naturales en
sistemas elctricos de potencia de prueba. Los registros
del comportamiento del sistema se obtienen mediante la
simulacin de fallas en dicho sistema. Posteriormente,
con la informacin contenida en tales registros, de la
respuesta en el tiempo del sistema ante la falla, se
calculan los modos naturales del sistema de prueba y se
comparan los resultados con los obtenidos mediante
tcnicas del anlisis modal.
Palabras clave: modos de oscilacin, anlisis lineal,

espectro de frecuencias, estabilidad, oscilaciones
electromecnicas, oscilaciones de baja frecuencia.
I. Introduccin
La respuesta en el tiempo de un sistema ante una
perturbacin es la superposicin de componentes de
diferentes frecuencias. Si se tiene la respuesta en el
tiempo de un sistema, a partir de dicha informacin es
posible determinar el espectro de frecuencias de alguna
seal as como los aspectos de fase.
II. Desarrollo
En el caso que se presenta se modela un SEP de prueba y
con la ayuda de un programa de Estabilidad Transitoria
se registra el comportamiento de la velocidad angular de
rotores de maquinas al ocurrir una perturbacin en uno
de los nodos del sistema. A partir de esa informacin se
realiza el anlisis de las seales muestreadas para obtener
las frecuencias que componen la respuesta de cada
generador mediante el clculo del espectro de frecuencias
de dichas seales.
Para ello pueden utilizase tcnicas apropiadas de

muestreo y anlisis de seales, la Transformada de
Fourier en Tiempo Discreto, Anlisis Prony, aplicadas a
la seal de respuesta [1, 2, 5]. Para aplicar esta teora a
un SEP debe contarse con los registros necesarios,
procesar adecuadamente las seales y realizar la
interpretacin de los resultados lo cual forma parte del
objetivo de este artculo.
En forma alterna, a partir de los parmetros de la red

elctrica, se presentan los resultados obtenidos a partir de
las tcnicas del anlisis lineal con las cuales se obtienen
los modos de oscilacin electromecnicos del sistema en
cuestin.
116
Con ambos resultados, los obtenidos a partir de los datos

del disturbio y los obtenidos a partir del anlisis lineal se
realizan comparaciones para observar si es posible
obtener las frecuencias naturales de oscilacin para un
sistema del cual se desconoce el modelo de la red y se
describen las observaciones del comportamiento del
sistema ante el disturbio y su relacin con las frecuencias
naturales de oscilacin.
tratando de evitar que puedan confundirse con los lbulos

laterales que se obtienen durante el procesamiento de la
seal. Con estas consideraciones, se obtiene los
resultados que siguen utilizando las tcnicas del anlisis
de Fourier.
En la figura 3 se muestra el espectro de frecuencias
obtenido a partir de los datos registrados de la velocidad
de cada maquina. En ella se pueden observar las
componentes principales de la respuesta de cada mquina.
Algunas frecuencias de oscilacin se presentan en varias
maquinas por ejemplo la frecuencia de 2.27hz est
presente en todos los generadores en diferentes
proporciones.
El sistema de la figura 1 se encuentra operando en estado

estable cuando se presenta una falla trifsica que se libera
en 1 ciclo.
G4
G1
G2
G5
G3
Figura 1 Sistema de Prueba Multimquina.

Se registra el comportamiento en diversos puntos del
sistema y en la figura 2 se muestra grficamente la
respuesta oscilatoria de la velocidad de las mquinas que
se obtiene una vez que la falla es liberada.
Velocidad angular
376.99
G1
G2
G3
G4
G5
376.98
w (rad/s)
376.97
376.96
Figura 3 Espectro de Frecuencias de las

Oscilaciones de mquinas.
376.95
376.94
Algunas frecuencias tienen una mayor magnitud lo cual

indica mayor influencia en una mquina en particular por
ejemplo en la respuesta de G1, la frecuencia de 2.27hz
tiene la mayor magnitud y domina el comportamiento
oscilatorio de la maquina G1, lo cual tambin se observa
en la figura 2, donde dicha respuesta tiene un periodo de
0.45seg. Se observa tambin como las frecuencias
naturales de oscilacin que componen a este sistema de
prueba participan en forma diferente en cada generador.
376.93
376.92
376.91
0
0.5
1.5
t(seg)
2.5
Figura 2 Oscilacin de velocidades de mquinas ante el

disturbio.
Ntese la no linealidad de la respuesta oscilatoria, G1

presenta una oscilacin casi lineal pero con variacin en
amplitud. G2 oscila de forma no lineal y con amplitud
decreciente. G3 oscila permanentemente con forma no
lineal. G4 oscila de forma altamente no lineal y G5 con
menor amplitud tiene una respuesta tambin no lineal.
La figura 4 muestra el espectro obtenido para la mquina

G1 donde se observa tanto la magnitud como la fase de
las frecuencias que componen la respuesta de la mquina
ante la falla. En la tabla 1 se muestran las frecuencias en
magnitud y fase que componen la respuesta de cada
generador donde se observa cual de dichas frecuencias
tiene mayor influencia en el comportamiento de una
mquina en particular, esto es, la componente de mayor
Las seales presentan una componente de directa que es

necesario eliminar previo al procesado de las muestras.
Posteriormente se elige una ventana de muestreo que
permita observar mejor las componentes frecuenciales
117
muestran la magnitud y fase de la participacin relativa

de cada mquina en cada modo de oscilacin.
magnitud, 2.27hz en G1, 1.04hz en G2, 1.11hz en G3,

1.05hz en G4 y 2.27hz en G5.
RESULTADOS DEL ANLISIS ESPECTRAL (FORMA CARACTERISTICA DE OSCILACIN EN

MAGNITUD Y FASE)
G1
G2
G3
G4
G5
f = 2.27 Hz
1.000 -173.13
0.111
6.87
0.134
6.88
0.107
6.86
0.659
6.87
f = 1.47 Hz
0.430
8.01
0.342
178.67
0.464
-171.7
0.319
-171.75
1.000
8.00
f =1.11Hz
0.042
-169.50
0.293
-11.40
1.000
-176.23
0.059
f = 1.05 Hz
0.586
15.1
1.000
166.86
-169.3
DESCRIPCIN DEL CASO: Sistema multimquina. Valores correspondientes a la oscilacin

de la velocidad angular
Tabla 2 Espectro de frecuencias de la respuesta oscilatoria de

los generadores G1, G2, G3, G4 y G5 ante la ocurrencia de un
disturbio. Las frecuencias se han ordenado para indicar en
cuales generadores est presente cada frecuencia. Por ejemplo
la frecuencia de 1.05hz solo est presente en G2 y G4.
De este modo se puede observar cual generador tiene

mayor participacin en una frecuencia dada, por ejemplo
G1 en 2.27hz, G5 en 1.47hz, G3 en 1.11hz y G4 en
1.05hz as como la participacin relativa de cada
mquina, sin embargo la mayor utilidad de esta
distribucin de los datos es la observacin de la forma
caracterstica de oscilacin de cada modo en funcin de
la fase; en 2.27hz la maquina G1 oscila en contra del
resto de las maquinas, en 1.47hz G1 y G5 oscilan en
contra de G2, G3 y G4, etc. Esto puede observarse tanto
en la tabla 2 como en la tabla 3.
Figura 4 Amplitud y Fase de componentes de la

mquina 1. Para f=2.27hz la fase es de -173.1,
para f=1.47hz la fase es 8 y para f=1.11hz la fase
es de -169.5.
Es importante el hecho de que estos resultados pueden

compararse entre todas las maquinas, por ejemplo, la
mayor amplitud de oscilacin ante el disturbio se
presenta en G3 con una magnitud relativa de 5.38 para la
frecuencia de 1.11hz. La tabla 1 es sobre todo til para
efectos de comparacin de magnitudes de oscilacin.
FREC
G1
MAG
FASE
FREC
G2
MAG
FASE
FREC
G3
MAG
FASE
2.268
4.975
-173.1
1.040
1.791
15.1
1.113
5.383
-176.2
1.473
0.335
8.0
1.120
1.576
-11.4
2.268
0.665
6.9
1.110
0.228
-169.5
2.268
0.553
6.9
1.473
0.361
-171.7
1.475
0.266
178.7

F OR M A C AR AC TER IS TIC A DE OS C ILAC IN EN M AGNITUD Y F AS E
Modo 1
f1=2.279 Hz
G1
1
90
G2 0.1124
-90
G3
0.136
-90
G4 0.1086
-90
G5 0.6523
-90
DESCRIPCIN DEL CASO
Sistem a m ultim quina.
FREC
G4
MAG
FREC
G5
MAG
FASE
1.053
3.055
166.9
FASE
2.268
3.281
6.9
2.268
1.473
0.530
0.249
6.9
1.473
0.778
8.0
-171.7
1.110
0.318
-169.3
Modo 2
f2 = 0 Hz
1
-90
1
-90
1
-90
1
-90
1
-90
Modo 3
f3 = 1.4772 Hz
0.4229
180
0.3366
0
0.4667
0
0.3184
0
1
180
Modo 4
f4=1.1179 Hz
0.0425
90
0.3192
-90
1
90
0.2496
-90
0.0598
90
Modo 5
f 5 = 1.0504 Hz
0.0129
-90
0.6974
-90
0.0772
-90
1
90
0.0173
-90
Tabla 3 Clculo de modos de oscilacin

del sistema. Anlisis lineal.
A partir de la tabla 2 se observa que, con ligeras

variaciones, la participacin relativa de las maquinas
ocurri de acuerdo a lo que pudiera haberse previsto del
anlisis modal presentado en la tabla 3. Al comparar
ambos resultados del anlisis de la respuesta ante el
disturbio mostrados en la tabla 2 y los mostrados en la
tabla 3 observamos que, el disturbio ocurrido en el
sistema provoc una respuesta oscilatoria en las
mquinas, compuesta por las frecuencias naturales del
sistema en cuestin. Conclusiones similares se obtienen
tambin acerca de la forma caracterstica de oscilacin.
Tabla 1 Espectro de frecuencias de la respuesta

oscilatoria de los generadores G1, G2, G3, G4 y
G5 ante la ocurrencia de un disturbio.
Se observa tambin que son 4 las frecuencias que, en

mayor o menor medida, contribuyen a la respuesta de
todas las maquinas del sistema; 2.27hz, 1.47hz, 1.11hz y
1.05hz y ellas son las frecuencias naturales del sistema,
como se aprecia en la tabla 3.
En la tabla 2 la informacin de la tabla 1 se ha
reacomodado y normalizado solo para fines de
comparacin con los resultados obtenidos en el clculo
de modos de oscilacin electromecnicos del sistema los
cuales se muestran en la Tabla 3. Los resultados
A partir de esta comparacin tambin es importante

mencionar que la matriz de participacin utilizada en el
anlisis lineal para mostrar resultados de participacin
neta de cada variable de estado en cada modo de
oscilacin puede compararse con el espectro de
118
frecuencias mostrado en la tabla 1. En la Tabla 4 se

muestran los factores de participacin obtenidos en el
anlisis lineal.
G1
G2
G3
G4
G5
III. CONCLUSIONES
f1=2.27 f3=1.47 f4=1.11 f5=1.05
0.349
0.060
0.001
0.000
0.007
0.063
0.088
0.192
0.005
0.052
0.374
0.002
0.005
0.041
0.036
0.306
0.134
0.285
0.002
0.000
1.
Con la metodologa presentada se han calculado

los modos de oscilacin de un sistema de
potencia de prueba a partir de conocer la
respuesta de los maquinas ante un disturbio
ocurrido en el sistema.
2.
Tambin con el anlisis de las seales de

respuesta en el tiempo, ha sido posible obtener,
con una muy buena aproximacin, los aspectos
de magnitud y fase de las frecuencias naturales
o modos de oscilacin electromecnicos del
sistema.
3.
Ha sido posible observar cual modo de

oscilacin tiene mayor influencia en el
comportamiento de un generador en particular y
tambin se pudo observar cual generador tiene
mayor participacin en una frecuencia dada.
4.
Con los resultados aqu obtenidos se confirma

que, con el conocimiento de los modos naturales
del sistema, es posible predecir con magnitudes
relativas y con buen acercamiento el
comportamiento del sistema ante disturbios.
5.
Se demuestra que, con la adecuada seleccin de

registros, es posible obtener las frecuencias
naturales de oscilacin de un sistema mediante
el anlisis de seales de la respuesta en el
tiempo de disturbios ocurridos en el sistema.
6.
Si en un sistema de potencia se piensa que hay

problemas de oscilaciones de baja frecuencia, se
puede implementar el monitoreo del mismo para
estudiar su comportamiento ante disturbios; a
partir de dicho comportamiento se pueden
identificar las frecuencias causa del problema
para posteriormente proporcionar en el sistema
el amortiguamiento adecuado.
7.
Esta herramienta puede ser de gran utilidad en la

supervisin del SEP.
Descripcin del caso: Sistema multimquina
Valores correspondientes a la oscilacin

de la velocidad angular
Tabla 4 Factores de participacin para la oscilacin de la
velocidad angular de mquinas. Resultados del anlisis
lineal.
Leyendo en forma horizontal, observamos en orden de

magnitud que 2.27hz, 1.47hz, 1.11hz son las frecuencias
que participaran en la oscilacin de G1. En la Tabla 1
este orden es 2.27hz, 1.47hz y 1.11hz.
Si continuamos con la comparacin pero utilizando la
informacin ahora presentada en la tabla 5 donde
solamente se cambia el arreglo de la tabla 4, podemos
comparar directamente con los resultados del anlisis
espectral (tabla 1). Adems con este arreglo es sencillo
observar que junto con los aspectos de fase ya
presentados (tabla 3) es posible predecir el
comportamiento no lineal de la respuesta oscilatoria de
las maquinas.
Factores de participacin. Sistema multimquina.
G1
Frec
2.27
1.47
1.11
1.05
Mag
0.349
0.060
0.001
0.000
G2
G3
G4
Frec
Mag
Frec
Mag
Frec
Mag
1.04
0.192
1.11
0.374
1.05
0.306
1.12
0.088
1.47
0.052
1.47
0.041
1.48
0.063
2.27
0.005
1.11
0.036
2.27
0.007
1.05
0.002
2.27
0.005
G5
Frec
1.47
2.27
1.11
1.05
Mag
0.285
0.134
0.002
0.000
Tabla 5 Factores de participacin para la oscilacin de la

velocidad angular de mquinas. Resultados del anlisis lineal.
Para el resto de los generadores no se cumple del todo

este orden sin embargo es conveniente mencionar que no
se espera que los resultados obtenidos sean exactamente
iguales, pero si muy aproximados, ya que en el caso
presentado, los registros analizados corresponden a una
falla trifsica liberada en un ciclo pero en el sistema
existen diferentes disturbios que pueden ocurrir pudiendo
excitar en forma diferente los modos oscilatorios. Sin
embargo se ha mostrado que es posible, con est tcnica,
obtener las frecuencias naturales de oscilacin de un
sistema, as como, con una muy buena aproximacin, los
aspectos de magnitud relativa y fase de la participacin
de las unidades generadoras en la oscilacin.
IV. BIBLIOGRAFA
119
1.
Alan V. Openheim, Ronald W. Schafer, Jhon R.

Back, Tratamiento de Seales en tiempo
Discreto, 2 Edicin, Prentice Hall, Madrid,
2000.
2.
Alan V. Oppenheim, Alan S. Willsky Seales y

Sistemas,
2
Edicin,
Prentice
Hall
Hispanoamericana, S.A., Mxico, 1998.
3.
F. Aboytes, Juan C. Ramrez, Anlisis de

oscilaciones de voltajes y potencias reactivas en
sistemas elctricos de potencia, RVP-2000,
Acapulco, Guerrero.
4.
J.F. Hauer, C.J. Demeure, L.L. Scharf, Initial

Results in Prony analysis of power system
response signals, IEEE 1990.
5.
F. Aboytes, F. Sanchez, Oscilaciones

electromecnicas en sistemas elctricos de
potencia, RVP-88, Acapulco, Guerrero, Agosto
1988.
6.
R. T. Byerly, D. E. Sherman, D. K. McLain,

Normal modes and mode shapes applied to
dynamic stability analysis, IEEE Transactions on
PAS vol. 94, March/April 1975.
7.
Florencio Aboytes, Estabilidad Transitoria en

Sistemas Elctricos de Potencia, Comisin
Federal de Electricidad, Mxico, D.F., Julio 1990.
8.
Juan C. Ramrez, Tesis de Maestra. Anlisis de

Modos de Oscilacin y Amortiguamiento en
Sistemas Elctricos de Potencia. DIE-FIMEUANL.
M.C. Luis David Archiga Rodrguez,

naci en Monterrey NL. Se gradu de
Ingeniero Mecnico Electricista en la
Facultad de Ingeniera Mecnica y
Elctrica (FIME) de la Universidad
Autnoma de Nuevo Len (UANL) en
1985, realiz sus estudios de maestra en el programa
Doctoral de Ingeniera Elctrica de la FIME de la UANL
de 1986 a 1991, en la actualidad trabaja en el ACNE de
la CFE, desempeando el puesto de Subgerente de
Transacciones Comerciales.
V. CURRICULUM VITAE
Ing. Juan Carlos Ramrez Saucedo.
Naci en Zacatecas en 1974, y desde 1986
radica en Nuevo Len. Actualmente trabaja
para la CFE en el rea de Control Noreste
y desde 1999 al 2008 fue catedrtico en la
Universidad del Norte en Monterrey, N. L.
Realiz sus estudios de Maestra en el Programa Doctoral
de la Facultad de Ingeniera Mecnica y Elctrica en la
Universidad Autnoma de Nuevo Len. Tambin ha
trabajado durante unos 6 aos para la industria privada,
entre otros aspectos, en el anlisis de calidad de energa y
en el diseo de Transformadores Industriales.
M.C. Alfonso Lpez Delgado: Naci en
Tres Valles Ver. en 1966. Se gradu como
Ingeniero Electromecnico en 1990 en el
Instituto Tecnolgico de Tuxtepec Oax.
Obtuvo el grado de Maestra en Ingeniera
Elctrica en 1999 en el Programa Doctoral
de la Facultad de Ingeniera Mecnica y Elctrica en la
Universidad Autnoma de Nuevo Len. Actualmente
labora en el rea de Control Noreste en Monterrey N. L.
120

Unl

Cargado por

Copyright:

Formatos disponibles

Unl

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Unl

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD AUTNOMA DE NUEVO LEN

FACULTAD DE INGENIERA MECNICA Y ELCTRICA

ANLISIS DE MODOS DE OSCILACIN Y AMORTIGUAMIENTO

UNIVERSIDAD AUTNOMA DE NUEVO LEN

ANLISIS DE MODOS DE OSCILACIN Y AMORTIGUAMIENTO

En este trabajo de tesis se analiza el origen y causas de las oscilaciones

vectores caractersticos y los factores de participacin, con las frecuencias y el

Capitulo 2. Caracterizacin de modos de oscilacin

Sistemas lineales y ecuaciones de estado. . . . . . . . . . . . . . . . . . .

Eigenvalores o valores caractersticos. . . . . . . . . . . . . . . .

Eigenvectores o vectores caractersticos. . . . . . . . . . . . . .

2.4 La respuesta en el tiempo de un sistema lineal. . . . . . . . . . . . . .

Capitulo 3. Interpretacin de modos de oscilacin en sistemas

3.2 Frecuencias naturales no amortiguadas. . . . . . . . . . . . . . . . . . . . .

Obtencin del modelo linealizado de las ecuaciones de

El clculo de modos como un problema de valores

3.3 Interpretacin de modos de oscilacin en sistemas elctricos de

Capitulo 4. Amortiguamiento de oscilaciones

4.2 Coeficientes de sincronizacin y amortiguamiento en SEPs. . . .

4.3 Anlisis de coeficientes de sincronizacin y amortiguamiento en

4.4 Obtencin del modelo lineal de un SEP multimquina. . . . . . . . . .

4.5 Excitacin de modos de oscilacin. . . . . . . . . . . . . . . . . . . . . . . . .

4.6 Interpretacin de modos de oscilacin en SEPs de prueba y

Interpretacin del amortiguamiento en los modos de

Capitulo 5. Oscilaciones de voltaje, potencia activa y potencia

5.2 Ecuaciones de voltaje, potencia activa y potencia reactiva. . . . . .

Ecuaciones de voltajes en terminales de generadores. . .

Ecuaciones de potencia activa y potencia reactiva. . . . . .

5.3 Anlisis de oscilaciones de voltaje y potencia. . . . . . . . . . . . . . . .

Capitulo 6. Determinacin de modos de oscilacin mediante el

Determinacin de modos de oscilacin en un SEP de

Capitulo 7. Conclusiones y recomendaciones

7.2 Recomendaciones para trabajos futuros. . . . . . . . . . . . . . . . . . . . .

Apndice 1. Sistemas de prueba de casos de estudio. . . . . . . . . . . . . .

Apndice 2. Artculos presentados del trabajo de tesis. . . . . . . . . . . . .

de oscilaciones y tener guas tcnicas para influir en el diseo de controles de

1.2. Revisin bibliogrfica

Aunque inicialmente se realizaron estudios de oscilaciones electromecnicas

en anlisis de residuos para la aplicacin de FACTS (Flexible Ac Transmision

dispositivos FACTS mejoran el desempeo dinmico de un sistema de

transformada de wavelet discreta y el algoritmo Prony.

1.3. Estructura de la tesis

investigacin, se presenta la estructura del contenido del trabajo y se incluyen

tcnica de valores caractersticos y se interpreta el significado dichos valores y

consideraciones para obtener los modos de oscilacin naturales y formas

operativa y del movimiento de las posiciones angulares de rotores de unidades

El captulo 6 presenta una alternativa para el clculo de los modos naturales de

Se describe en forma clara, la relacin entre los valores caractersticos,

Los sistemas elctricos de potencia (SEP) estn sujetos a diversas

del anlisis dinmico de un sistema es de gran importancia identificar el origen

Sistemas lineales y ecuaciones de estado

Las ecuaciones dinmicas que rigen el comportamiento del sistema de potencia

donde x (t ) denota el vector de derivadas de las variables x(t) con respecto al

En el anlisis de sistemas de potencia se considera un tiempo fijo t 0 inicial, y las

2.3.1 Modos de oscilacin

La matriz [SI-A] es algunas veces llamada la matriz caracterstica de A, su

es conocido como el polinomio caracterstico de A, y sus races {i} son

para una solucin no trivial