UNIVERSIDAD AUTÓNOMA DE NUEVO LEÓN

FACULTAD DE INGENIERÍA MECÁNICA Y ELÉCTRICA

DIVISIÓN DE ESTUDIOS DE POSGRADO

iLEKfUttOlM
VERTUTE

A N A L I S I S DE RESONANCIA SUBSINCRONA Y
CONTRAM! E DIDÀS

EN OPCION AL GRADO DE
MAESTRO EN CIENCIAS
DE LA INGENIERÍA ELÉCTRICA C
ESPECIALIDAD EN POTENCIA

PRESENTA:
ALEJANDRO MENDOZA VIVEROS

SN NICOLAS DE LOS GARZA, N, L. MAYO DEL 2004

T M
2 5 8 5 3
. M2
F IME
2 0 0 4

. M 4 7
1 0 2 0 1 5 0 0 4 8
UNIVERSIDAD AUTONOMA DE NUEVO LEON

FACULTAD DE INGENIERÍA MECÁNICA Y ELÉCTRICA

DIVISIÓN DE ESTUDIOS DE POSGRADO

ANALISIS DE RESONANCIA SUBSINCRONA Y

CONTRAM EDIDAS

TESIS

EN OPCION AL GRADO DE
MAESTRO EN CIENCIAS
DE LA INGENIERÍA ELÉCTRICA CON
ESPECIALIDAD EN POTENCIA

PRESENTA:
ALEJANDRO MENDOZA VIVEROS

SN NICOLÁS'DE LOS GARZA, N. L. MAYO DEL 20Û4.

 Universidad A u t ó n o m a de N u e v o L e ó n

F a c u l t a d de Ingeniería M e c á n i c a y E l é c t r i c a

Subdirección de Estudios de Postgrado

Los miembros de comité de tesis recomendamos que la tesis A N A L I S I S DE

RESONANCIA SUBSÍNCRONA Y CONTRAMEDIDAS, realizada por el alumno

Alejandro Mendoza Viveros, matrícula 1 1 1 5 0 4 6 sea aceptada para su defensa como opción al

grado de Maestro en Ciencias de la Ingeniería Eléctrica con especialidad en Potencia.

El Comité de Tesis

Asesor
Dr. Fernando Sánchez Tello
 AGRADECIMIENTOS

Primeramente a Dios por permitirme culminar esta etapa tan importante en

mi vida profesional y personal. A mis padres Dionisio Mendoza Vázquez y
Yolanda Viveros Vivaldo por haberme formado con tanto empeño y cariño. A
ti Leonardo por ser un buen amigo y un excelente hermano. Gracias por las
palabras de apoyo en los momentos difíciles.

A ti Ana parte fundamental de mi Vida, por todo el cariño, el apoyo,

comprensión y paciencia que me haz brindado siempre, gracias y mil
disculpas por todos esos fines de semana que te hice pasar en el postgrado

A los Ingenieros Lino Coria Cisneros y José Luis Lemus Duarte

A todos mis maestros y compañeros del postgrado por haberme compartido

sus conocimientos y amistad

Al Dr. Fernando Sánchez Tello por su apoyo y amistad a lo largo de mis

estudios de postrado, así como sus valiosos consejos que fueron sin duda
claves para la culminación del trabajo. Gracias Dr. por su paciencia.

A los Doctores Ernesto Vázquez Martínez y Rogelio Ramírez Barradas por

sus valiosos comentarios y observaciones

A los Ingenieros Emilio Luna Quiroz y Francisco Pérez Ponce de León por
haberme facilitado los medios para la impresión de este trabajo

Al CONACYT por el apoyo brindado en mis estudios de postgrado

Dedico esta tesis a todos aquellos que a lo largo de mi vida me han formado
y que tan generosamente me han transmitido sus conocimientos. Gracias por
confiar en mí

/
 RESUMEN

Publicación No

Alejandro Mendoza Viveros, M.C. en Ingeniería Eléctrica

Universidad Autónoma de Nuevo León, 2003

Profesor Asesor: Dr. Fernando Sánchez Tello

El trabajo presenta un análisis conceptual del problema de resonancia subsíncrona y

estudia una interacción de diversos elementos de control con el sistema mecánico del

generador. Se ilustra el efecto relevante de algunos parámetros del sistema en el

amortiguamiento de oscilaciones torsionales.

El problema de resonancia subsíncrona se aborda desde el estudio de sus causas

primarias hasta el diseño de una solución general problema. El material presentado debería

ayudar a aquellos involucrados en la planificación y diseño de sistemas eléctricos a entender

la naturaleza y alcance del problema de resonancia subsíncrona para el mejor

aprovechamiento de la red de trasmisión y sus controles.

Se ilustra el detalle del modelado utilizado para el estudio de resonancia subsíncrona

y se muestra su utilización con un sistema de prueba típico del IEEE.

Se estudian los efectos de la compensación serie paralelo de la red de transmisión

con los modos d oscilación mecánicos del generador. Con el mismo fin se evalúa el impacto

del sistema de excitación del generador , el compensador estático de VAR's y el estabilizador

de potencia y el dispositivo de compensación serie controlada por tiristores (TCSC por sus

siglas en ingles ).

Para el TCSC, se propone una técnica de control modal, la cual permite estabilizar y

mejorar el amortiguamiento de cualquier frecuencia torsional en el sistema, para un rango

amplio de compensación serie en el sistema.

En el análisis de todos los casos presentados se emplean las técnicas de análisis

modal y de eigenvalores.
 INDICE
AGRADECIMIENTOS V

RESUMEN vi

ÍNDICE vii

CAPÍTULO I INTRODUCCIÓN 1

1.1 MOTIVACIÓN 1

1 . 2 REVISIÓN HISTÓRICA 2

1 . 3 CONTENIDO DEL TRABAJO 3

1 . 4 OBJETIVOS 4

CAPÍTULO II CONCEPTOS BÁSICOS Y M O D E L A D O DE ELEMENTOS.5

2 . 1 INTRODUCCIÓN 5

2 . 2 RESONANCIA 6

2 . 3 COMPENSACIÓN SERIE EN LÍNEAS DE TRANSMISIÓN 6

2 . 4 MODULACIÓN DE FRECUENCIA 8

2.4.1 efecto de generador de inducción 12

2.4.2 interacciones torsionaies 13
2.4.3 pares transitorios 13
2.5 MÉTODOS DE ANÁLISIS 14

2.5.1 análisis de eigenvalores 14

2.5.2 barrido de frecuencia 15
2.5.3 análisis de respuesta en la frecuencia 15
2.5.4 análisis aproximado en el dominio de la frecuencia 15

2.5.5 análisis en el dominio del tiempo 15

2 . 6 MODELADO DEL SISTEMA ELÉCTRICO 16

2.6.1 el modelo de la máquina 16

2.6.2 el modelo de la red 23

2.6.3 el sistema de excitación 24

2.6.4 el estabilizador 24
[Link] el filtro torsional 26

2.7 MODELADO PARA EL SISTEMA MECÁNICO 27

CAPÍTULO III EVALUACIÓN DEL EFECTO DE GENERADOR DE

INDUCCIÓN Y LA RESONANCIA S U B S Í N C R O N A 30

3 . 1 INTRODUCCIÓN 30

3.2FORMACIÓN DE LA MATRIZ DE ESTADOS DEL SISTEMA 31

3.2.1 el modelo incremental de la red 31

3.2.2 efecto de! rotor de la máquina 32
3.2.3 el acoplamiento entre el sistema eléctrico y el sistema mecánico 32
3 . 3 ANÁLISIS DE OSCILACIONES AUTOEXCITADAS 34

3.3.1 evaluación de! efecto de generador de inducción 34

3.3.2 análisis de interacciones torsionales 41
[Link] forma incremental de las ecuaciones de! sistema mecánico 41

3 . 4 EL P S S Y SU INTERACCIÓN CON EL SISTEMA MECÁNICO 44

CAPÍTULO IV EL COMPENSADOR ESTÁTICO DE VARS Y SU

INFLUENCIA SOBRE LOS MODOS TORSIONALES 47

4 . 1 INTRODUCCIÓN 47

4 . 2 REGULACIÓN DE VOLTAJE EN EL PUNTO MEDIO DE LA LÍNEA 48

4 . 3 INFLUENCIA DE LA COMPENSACIÓN PARALELO SOBRE LA CARACTERÍSTICA DE TRANSMISIÓN DE

POTENCIA 49

4 . 4 ANÁLISIS DE ESTABILIDAD PARA LOS MODOS TORSIONALES DEL GRUPO TURBINA-GENERADOR .. 5 0

4.4.1 el modelo de la red sin compensación 50

4.4.2 modelado de la línea con un CEV conectado en el punto medio 51
4.4.3 obtención del modelo en espacio de estado del CEV 55
4.4.4 análisis de valores característicos 56

CAPÍTULO V EL TCSC Y SU INTERACCIÓN CON EL SISTEMA

MECÁNICO 58

5 . 1 INTRODUCCIÓN 58

5 . 2 MODELO DEL T C S C EN ESTADO ESTACIONARIO 59

5 . 3 IMPEDANCIA FUNDAMENTAL DEL T C S C 69

5 . 4 PRINCIPIOS BÁSICOS DE OPERACIÓN 72

5.4.1 operación en corto circuito 73

5.4.2 modo de operación en bloqueo 73

5.4.3 modo de operación vernier 74
5.4.4 problemas de diseño 74

vi i 1
 5 . 5 VENTAJAS Y APLICACIONES 74

5.5.1 control de flujo de potencia activa y reactiva 75

5.5.2 mejoramiento de la estabilidad transitoria 75
5.5.3 control de flujo no deseado ( loop flow ) 75

5.5.4 amortiguamiento de oscilaciones subsíncronas 75

5 . 6 ESTRATEGIA DE CONTROL CONVENCIONAL 76

5 . 7 ESTRATEGIA DE CONTROL MODAL 76

5.7.1 cálculo de componentes modales 77

5 . 8 DESCOMPOSICIÓN DE LA SEÑAL UTILIZANDO FILTROS ACTIVOS 78

5 . 9 MODELO MATEMÁTICO PARA EL T C S C PARA ESTUDIOS DE RESONANCIA SUBSÍNCRONA 79

5.9.1 aproximación de la característica de reactancia de! TCSC 79

5.9.2 Obtención de las ecuaciones ilnealizadas para el TCSC 81

5 . 1 0 ANÁLISIS DE EIGENVALORES 81

CAPÍTULO VI OTROS ESQUEMAS DE A M O R T I G U A M I E N T O PARA

OSCILACIONES SUBSÍNCRONAS ....85

6 . 1 INTRODUCCIÓN 85

6 . 2 ESQUEMA DE AMORTIGUAMIENTO N G H - S S R 86

6 . 3 E L T C B R ( THYRISTOR CONTROLLED BRAKING RESISTOR ) 89

6.3.1 amortiguamiento de oscilaciones subsíncronas 90

6 . 4 EL ESTABILIZADOR DINÁMICO 92

CAPITULO VII CONCLUSIONES Y RECOMENDACIONES 95

7 . 1 INTRODUCCIÓN 95

7 . 2 CONCLUSIONES Y RECOMENDACIONES 96

7 . 3 APORTACIONES 97

7 . 4 TRABAJOS FUTUROS 97

APENDICE A 99

APENDICE B 101

APENDICE C 102

REFERENCIAS 104

INDICE DE FIGURAS 106

INDICE DE TABLAS 108

 CAPÍTULO I

INTRODUCCION

1.1 MOTIVACIÓN

Los sistemas eléctricos de potencia (SEP's) experimentan una gran variedad de

oscilaciones las cuales son excitadas durante perturbaciones y maniobras operativas de
rutina. La mayoría de estas oscilaciones desaparecen debido al amortiguamiento positivo
neto del sistema. Sin embargo, bajo algunas circunstancias, una oscilación específica puede
alcanzar una magnitud considerable, que puede crecer o mantenerse por un período
prolongado de tiempo, provocando un daño al aislamiento de las líneas de transmisión, una
disminución en el tiempo de vida útil de las partes mecánicas de los equipos, un daño
permanente a los generadores ó la pérdida de estabilidad del sistema.

Dados los parámetros de las líneas de transmisión y el equipo, las oscilaciones del
sistema que involucran solo componentes eléctricos pasivos, presentan frecuencias
resonantes substancialmente mayores a la frecuencia nominal del sistema. Sin embargo,
cuando las oscilaciones involucran el equipo mecánico giratorio y el eléctrico, acoplados
estos a través del flujo magnético, aparecen frecuencias menores a la frecuencia nominal.
Las oscilaciones que involucran la masa e inercia del grupo turbina-generador presentan
frecuencias de oscilación electromecánicas en el rango de 0.1 hz hasta varios hertz.

Las oscilaciones subsíncronas con frecuencias menores a la frecuencia síncrona, se

encuentran en el rango de 10-50Hz y resultan de oscilaciones mecánicas entre la masa de
turbinas individuales y el generador acoplado en la misma flecha, generalmente de gran
longitud. Estas oscilaciones mecánicas, están eléctricamente acopladas con el sistema a
través del rotor de la máquina. El problema de resonancia subsíncrona comenzó a estudiarse
a partir del daño físico que se presentó en la flecha del grupo turbina-generador de la planta
Mojave de la Southern California Edison [1].
 Por otra parte, se ha reconocido que todos los controles activos de alta velocidad en
un SEP tales como líneas en corriente directa de alto voltaje (HVDC), los controladores
FACTS y el control de excitación de los generadores tienen la capacidad para mitigar el
problema [10]. Sin embargo, una mala sintonización de estos dispositivos puede causar una
disminución en la vida útil de los generadores o en el caso más drástico, la pérdida del
mismo.

1.2 REVISIÓN HISTORICA

Analíticamente el problema de resonancia subsíncrona es complejo, existen muchos

trabajos publicados orientados hacia el tratamiento del problema de resonancia subsíncrona,
sin embargo es difícil encontrar una fuente que describa con detalle los aspectos básicos del
problema.

La primera experiencia de resonancia subsíncrona que se tiene registrada es de la

década de los 70's en una planta al oeste de los Estados Unidos, en donde se tuvo pérdida
total de la máquina. A partir de ese momento, se empezó a trabajar en la identificación del
problema. Este trabajo fue publicado en 1976 por Hall M.C. y Dodges, D.A. [1],
Posteriormente se publicó el modelo de prueba propuesto por el IEEE para estudios de
resonancia subsíncrona en 1977 [2]. Para 1980 se publicó un artículo en el que se trataron
por primera vez las interacciones torsionales con los enlaces de corriente directa (HVDC)
conectados a turbo generadores [4], En 1982 I.M . Cannay [5] propone una nueva
aproximación para el amortiguamiento y la interacción torsional, lo que aportó nuevas
herramientas para el análisis de amortiguamiento en el sistema mecánico y eléctrico de la
red, así como las posibles frecuencias naturales de la flecha con posibilidad de inestabilidad.
A lo largo de este periodo y posterior a él, se han publicado infinidad de artículos sobre el
tema.

En años recientes destacan dos publicaciones relacionadas al tema de resonancia

subsíncrona, en 1990 [7] P.M. Anderson, [Link] y [Link] Ness publican un libro que
trata a detalle el modelado de elementos orientado al estudio de resonancia subsíncrona en
el dominio de la frecuencia. En 1999 el Dr. K.R. Padlyar hace nuevas aportaciones a la
bibliografía incluyendo en su libro modelos matemáticos para el análisis de dispositivos
FACTS, y enlaces de corriente directa, mediante una técnica de eigenvalores [10].
1.3 C O N T E N I D O DEL TRABAJO

El trabajo desarrolla un tratamiento conceptual del tema de resonancia subsíncrona y

está dividido en 7 capítulos. A lo largo de cada capítulo se cubren diversos aspectos, tales
como modelado, análisis, evaluación y propuestas de solución al problema. Inicialmente se
¡lustra el proceso de modelado realizando un análisis matemático detallado de cada
elemento involucrado. La evaluación del problema de Resonancia Subsíncrona se realiza
utilizando técnicas de eigenvalores. El análisis se respalda por un buen número de ejemplos
ilustrativos.

La distribución de los temas a tratar está hecha de la siguiente manera:

En el capítulo 2 se revisan conceptos básicos relacionados con el problema de

resonancia subsíncrona. Se obtienen los modelos en espacio de estado para el generador, la
flecha del grupo turbina-generador, el excitador y el estabilizador.

El capítulo 3 presenta la linealización de los modelos matemáticos desarrollados en el

capítulo 2, se plantea la matriz de estados del sistema y se estudia el fenómeno de
generador de inducción analíticamente.

El capítulo 4 presenta el comportamiento del compensador estático de VAR's (CEV) y

su influencia sobre los modos torsionales de la flecha. Se presentan los modelos
matemáticos utilizados para la red de transmisión sin compensación y con un CEV
conectado en el punto medio de la línea de transmisión.

El capítulo 5 presenta el detalle de modelado de la compensación serie controlada por

tiristores (TCSC, por sus siglas en ingles) y un análisis de su comportamiento y aplicación
para mejorar al estabilidad del sistema y amortiguamiento de oscilaciones torsionales. Se
muestra la aplicación de la técnica de control modal en la solución del problema de
resonancia subsíncrona.

El capítulo 6 informa de diversas técnicas actualmente empleadas en el

amortiguamiento a las oscilaciones torsionales.
 El capítulo 7 presenta las conclusiones del trabajo, así como las contribuciones y
posibles trabajos que pueden realizarse tomando como base los conceptos y modelos aquí
presentados.

Finalmente Los parámetros para el sistema utilizado se presentan en el apéndice A, el

apéndice B muestra el diagrama fasorial correspondiente y las relaciones utilizadas para el
cálculo de condiciones iniciales utilizadas en este trabajo. El apéndice C muestra una
metodología para la sintonización del control para el TCSC

1.4 OBJETIVOS

• Presentar un análisis conceptual sobre el problema de resonancia subsíncrona, de

una manera clara mediante ejemplos numéricos y el análisis de eigenvalores.

• Ilustrar en detalle el modelado del sistema mecánico, de la red así como de cada uno
de los diferentes componentes requeridos para el estudio de resonancia subsíncrona.

• Presentar una metodología, que ayude a la solución del problema de auto excitación

asociado con la red, utilizando el sistema de excitación de la máquina.

• Presentar de una manera clara y detallada un modelo lineal para el TCSC utilizado

para el cálculo de eigenvalores en estudios de resonancia subsíncrona.

• Realizar un análisis conceptual detallado del TCSC en un sistema máquina bus

infinito.

• Presentar una técnica de control modal para la desintonización del fenómeno de

resonancia subsíncrona logrando una solución global del problema de resonancia
subsíncrona, independientemente del nivel de compensación serie del sistema.

• Aplicación de técnicas de análisis modal y de eigenvalores orientadas al estudio de

resonancia subsíncrona.

/
 CAPÍTULO II

CONCEPTOS BÁSICOS Y MODELADO DE ELEMENTOS

2.1 INTRODUCCIÓN

Debido a la creciente demanda de energía eléctrica y la dificultad que se tiene

actualmente para la construcción de nuevas líneas de transmisión, la compensación serie en
los sistemas de potencia modernos ha surgido como una solución económica al problema de
capacidad de transmisión de potencia. Sin embargo, existen problemas derivados
caracterizados por oscilaciones de baja frecuencia que tienen como consecuencia la
interacción entre el sistema eléctrico y el sistema mecánico del grupo turbina generador.
Este efecto, que es conocido desde la década de los 70's y que se ha referido en la literatura
como resonancia subsíncrona [1], puede causar daño permanente a la flecha dependiendo
de la magnitud de las oscilaciones.

Existen dos problemas potenciales relacionados con el problema de oscilaciones de

baja frecuencia o resonancia subsíncrona. El primero de ellos se debe a la autoexcitación
provocada por las corrientes circulantes en la red de frecuencia subsíncrona y el segundo
debido a pares transitorios provocados por eventos en la red. El primero es un problema
escencialmente de estado estable, y el segundo un problema dinámico.

En este capítulo se revisan algunos conceptos básicos relacionados con el problema

de resonancia subsíncrona, se obtienen los modelos matemáticos para la máquina, el
excitador, la red, el excitador y el estabilizador del sistema de potencia así como su
representación en espacio de estado, Los modelos para el Compensador Estático de VARs
( CEV ) y la compensación serie controlada por tiristores ( TCSC por sus siglas en ingles) se
presentan en los capítulos 4 y 5 respectivamente.
2.2 RESONANCIA

La resonancia es una característica de los sistemas físicos, la cual puede o no ser

deseable, dependiendo del propósito del sistema.

En sistemas eléctricos, la resonancia se presenta en circuitos que tienen tanto

inductores como capacitores y aún cuando el teorema de máxima transferencia de potencia
dice que esta se obtiene cuando se llega a la condición resonante, esta condición de
operación debe evitarse ya que puede desencadenar daños irreversibles en equipos
conectados al sistema.

La real academia de la lengua da la siguientes definiciones para el fenómeno de resonancia:

Resonancia .- Aumento elevado de la amplitud de una oscilación bajo la influencia de

una acción periódica de una frecuencia próxima.

Resonancia Eléctrica.- fenómeno producido en un circuito oscilante cuando es

alimentado por una tensión alterna de frecuencia próxima a su frecuencia de oscilación.

2.3 LA COMPENSACIÓN SERIE DE LÍNEAS DE T R A N S M I S I Ó N

En un sistema de transmisión no compensado, como resultado de las fallas y otros

disturbios, aparecen componentes de cd que desplazan las corrientes de los devanados de
estator del generador. El resultado de este desplazamiento es la aparición de componentes
de par de entrehierro a una frecuencia de deslizamiento de 60 Hz. En el caso de fallas
desbalanceadas, las componentes de secuencia negativa, resultan en pares de entrehierro
de 120 Hz, por lo que es necesario evitar frecuencias que se encuentren cercanas a los 60 y
120 Hz.

En el caso de los sistemas compensados en serie, la situación es diferente, ya que a

diferencia de la componente transitoria de cd en la corriente de falla en sistemas no

compensados, puede presentarse una componente transitoria de corriente alterna de

frecuencia igual a la frecuencia natural (con) derivada de la interacción entre los elementos
/

capacitivos e inductivos de la red ( LC ), cuya expresión esta dada por :

Si se multiplica y se divide la expresión anterior por co0,

a
con = . ° = = con I^ rad/seg

t = = t Hz (2.2)

Para el circuito de la Figura 2,1,

f0 = es la frecuencia síncrona en Hz, Xc = Reactancia de Dispersión del

Oí. = 2x f rad/seg Transformador

X, = Reactancia Inductiva de la Línea

Xe=X +xt+xL
Xc = Reactancia Capacitiva de
X" = Reactancia Subtransitoria
Compensación
de la Máquina

En sistemas compensados en serie pueden presentarse oscilaciones subsíncronas

amortiguadas negativamente. Estas oscilaciones pueden presentarse como una condición de
autoexcitación o bien como una interacción entre el sistema mecánico y el sistema eléctrico,

1 X,L

/ "V

Figura 2.1.- Sistema de Transmisión Compensado en Serie

2.4 MODULACIÓN DE FRECUENCIA

La respuesta en el tiempo para un circuito sobreamortiguado en el que se tienen

inductores, resistencias y capacitores, está definida en la ecuación (2.3):

/(f) = /<[/lsin(cy1f + + Se"5"*' S\n(o)2t+ y/2)] (2.3)

En la ecuación (2.3) se tienen presentes dos frecuencias diferentes. La primera de

ellas es una componente de frecuencia igual a la nominal del sistema, mientras que la
segunda componente, corresponde a una señal senoidal amortiguada, cuya frecuencia
depende exclusivamente de los elementos del sistema.

Lo anterior se puede comprobar, mediante el análisis de una red R-L-C simple

conectada en serie. Si se supone que el voltaje de alimentación para el circuito descrito
anteriormente tiene la forma:

v(t) = V21/sen (a^ + e) (2.4)

Este voltaje puede transformarse al marco de referencia de la frecuencia de la

siguiente manera:

_ s [sen (0)1 + co. feos (0)1

V(s) = y/2V-± ^ y= (2.5)
S + ÍUj

Por otra parte, la impedancia del circuito en el marco de referencia de la frecuencia

esta dada por:

Z(s) = R + sL + ~ (2.6)
sC

de (2.5) y (2.6) puede obtener la expresión para la corriente:

 SV(S)
V(S) __ —
(2.7)
[
' Z (si „ + sL +
R . —
1 _s2 .+R.— S +1
2
SC L LC

La ecuación (2.7) puede representarse de la siguiente manera:

sV(s) sl/(s)

CÙ-,

En donde se puede definir:

03n = na
yTt' frecuencia tural no amortiguada

R IC re ac
C ' '°n am
ortiguamiento

a = -fój n El grado de amortiguamiento

á>2 = ú)nyJl~Cz Frecuencia amortiguada

Sustituyendo (2.5) en (2.8) y se resuelve la ecuación para la corriente se tiene:

i(t) = /c[/4sin(íü1t + y/x) + Be'<úi2t sin(<a2t +- (2.9)

La expresión anterior es idéntica que la ecuación (2.3). Para la ecuación (2.9) se tiene
lo siguiente:

y¡2co2CVcos(0)
K~
\[CÚ\ - o>l) + 4f(ú\cü'2
A = co2yjl + a2 col

1 - 2a£co2 + a¿<x>l
B = co.

r \
-2dCo)¡
- tan
\
2£o>¡ + cú\ - co

En las expresiones anteriores puede apreciarse las dos frecuencias relacionadas con la

corriente.

o)^ Frecuencia del sistema

<x>2 Frecuencia relacionada con los elementos de la red

Se puede definir a <x>2 de la siguiente manera:

1 (AL- R2C
\

i ¡= —
co2 = dú)
n
2L , C

Si se supone que el resultado anterior es la solución de la fase a perteneciente a una

red trifásica balanceada, entonces, la solución para las fases b y c será la misma, pero con

coeficientes diferentes en la componente de respuesta transitoria. Estas componentes se

verán afectadas al ser multiplicadas por las funciones sen{9) y cos(0), con sus respectivos

desplazamientos de ±120 grados para el caso de las fases b y c al aplicarles la

transformación de Park. En este caso, la componente transitoria de la corriente debida a los

elementos almacenadores de energía se transforma y aparecen términos de la forma:

cos(é i )sen((y 2 t)

Donde 0 es redefinido de la manera siguiente;

e = coxt + C +
La definición anterior permite mostrar que:

eos(0)sen(co 2 t) = — \sen{0 + co2t) - sen{0 - <v2t)]

2
(2.10)
= -{cos[(® 1 +ffl 2 )t + fl-[(®1 - f f l 2 ) t + f]}

Finalmente se puede concluir que para un circuito como el de la Figura 2.1, las

componentes transitorias de las corrientes circulantes en la red, se ven reflejadas en el rotor

de la máquina (esto como resultado de la transformación de Park) como corrientes

circulantes cuya frecuencia es la suma y la diferencia de las frecuencias a>xy co2, lo cual

provoca la aparición de pares subsíncronos y supersíncronos, de los cuales, los primeros son

de particular interés en estudios de resonancia subsíncrona, debido a que las frecuencia

naturales de oscilación de la flecha de la máquina se encuentran en este rango de

frecuencias. La consecuencia de estas corrientes (de frecuencia subsíncrona), es la inyección

de energía a la masa de la flecha. Esta energía, eventualmente puede presentar oscilaciones

muy cercanas a las frecuencias naturales resonantes de la red, desencadenando un

movimiento oscilatorio que resulta en grandes esfuerzos para la flecha del generador que

pueden ser poco amortiguados. Si estas corrientes se mantienen por un periodo de tiempo

prolongado, pueden causar fatiga al material de la flecha y aumentan las posibilidades de

falla mecánica en el generador.

El IEEE define al fenómeno de resonancia subsíncrona de la siguiente manera:

" Resonancia subsíncrona es la condición del sistema eléctrico de potencia en la que

existe un intercambio de energía entre la red eléctrica y la turbina de un generador en una o

más de las frecuencias naturales del sistema, por debajo de la frecuencia síncrona del

sistema".

Como ya se ha dicho, el fenómeno de resonancia subsíncrona puede estudiarse como

oscilaciones auto-excitadas o pares transitorios. En el estudio de oscilaciones autoexcitadas

se consideran dos casos. El efecto de generador de inducción, en donde la auto excitación

ocurre en el modo de la red y el estudio de interacciones torsionales, caracterizado por la

interacción entre el sistema eléctrico y el sistema mecánico del generador.

2.4.1 El efecto de generador de inducción

El efecto de generador de inducción es provocado por la auto excitación del sistema

eléctrico como resultado de una f.m.m. producida por las corrientes de frecuencia

subsíncrona.

El deslizamiento para la máquina síncrona operando como generador de inducción

esta dado por la expresión:

fn

En donde:

fn = Velocidad del rotor

f0 = Velocidad Síncrona

Debido a que fn < f0 el deslizamiento es negativo, la resistencia del rotor a frecuencia

subsíncrona es vista como una resistencia negativa desde las terminales de la máquina. Si la

sumatoria entre esta resistencia "negativa" y la resistencia externa del sistema ( compuesta

por la resistencia de armadura, la resistencia del transformador y la resistencia de la línea

de transmisión) es negativa, entonces las corrientes de frecuencia subsíncrona que circulan

por la red dan paso a una condición de operación inestable. Este fenómeno es conocido en la

literatura como el efecto de generador de inducción, debido a que la máquina se comporta

como un generador de inducción.

^ Xs rr Xr

Xm Í (1 -S)
^ < se desprecia ^ '-rr
A

r =r +
" > s

Figura 2.2.- Circuito Equivalente Simplificado para una Máquina

Síncrona como Función del Deslizamiento S
 La Figura 2.2 muestra un circuito simplificado de una máquina síncrona similar al

utilizado para motores de inducción. El efecto de saliencia se ha despreciado, los efectos del

circuito de rotor se han concentrado y se representan mediante una resistencia r r y una

reactancia de dispersión X r . Únicamente se representan los efectos del deslizamiento S y de

la frecuencia subsíncrona fn.

2.4.2 Interacciones torsionales

Las interacciones torsionales en un sistema eléctrico de potencia aparecen cuando las

oscilaciones del rotor del generador de frecuencia fn inducen componentes de voltaje a una

frecuencia fern dada por:

Si la componente subsíncrona fem se encuentra cerca de la frecuencia eléctrica de la

red fer, aparecen pares subsíncronos producidos por las componentes de voltaje

subsíncrono. Estos pares pueden ser sostenidos o crecientes tal que superan el

amortiguamiento mecánico del sistema giratorio provocando una condición de auto

excitación.

A la condición de operación en donde interactúan el sistema mecánico del generador

y el sistema eléctrico se le conoce como interacciones torsionales,

2.4.3 Pares transitorios

Las maniobras en la red como la liberación de alguna falla en el sistema pueden

provocar pares eléctricos transitorios.

En un par eléctrico transitorio existen muchas componentes con comportamientos

exponenciales de decaimiento así como pares oscilatorios de múltiples frecuencias que van

desde las subsíncronas hasta la frecuencia nominal del sistema. En sistemas de transmisión

sin compensación serie, estos transitorios resultan en respuestas de tipo exponencial, cuya

constante de decaimiento depende de la relación inductancia resistencia del sistema. Sin

embargo, en redes de transmisión compensadas en serie las corrientes transitorias tienen

una forma similar a la que se muestra en la ec (2.9) y el número de frecuencias oscilatorias

depende de la inductancia, la resistencia y la capacitancia de la red de transmisión.

Una red de transmisión con muchos capacitores serie instalados tendrá muchas

frecuencias subsíncronas y si alguna de estas frecuencias de la red se encuentra lo

suficientemente cerca o llega a coincidir con alguna de las frecuencias naturales del grupo

turbina generador se pueden presentar pares muy grandes debido a que son directamente

proporcionales a la magnitud de la corriente de oscilación.

La presencia de estos pares disminuye la vida útil de la flecha como resultado de la

fatiga que experimenta el material. La fatiga es un proceso de cambio estructural

permanente que se presenta en materiales sometidos a condiciones de esfuerzo y tensión

en alguno o algunos puntos. Estos esfuerzos pueden desencadenar en crujido o fractura

completa del material después de un número suficiente de eventos en que los pares tienen

magnitudes elevadas.

La fatiga acumulada puede desencadenar en un umbral donde la probabilidad de la

iniciación de crujido en áreas de esfuerzo elevado sea significante. Una vez iniciado el

crujido, este puede ser tal que resulte en daño irreversible y puede presentarse hasta la

ruptura de la flecha.

2.5 METODOS DE ANÁLISIS

Los métodos de análisis en el dominio de la frecuencia pueden ser usados para

facilitar el análisis de sistemas muy complejos. Estos sistemas conservan las características

resonantes de las redes que representan. También existen métodos en el dominio del

tiempo que pueden ser utilizados para el estudio de resonancia subsíncrona. Esta sección

describe brevemente algunos de estos métodos.

2.5.1 Análisis de eigenvalores [19]

Este método da información relacionada las diferentes frecuencias involucradas en el

modelo, así como el amortiguamiento de cada frecuencia del sistema. Para este tipo de

análisis se debe modelar la red y los generadores mediante un sistema de ecuaciones

diferenciales linealizadas puestas en la forma:

 X = Ax + Bu

Los eigenvalores se obtienen a partir de la matriz de estados (A) del sistema. La parte

imaginaria de los eigenvalores representa la frecuencia de oscilación del modo, mientras que

la parte real hace referencia al amortiguamiento de esta oscilación. Si la parte real de todos

los eigenvalores es positiva, el sistema es inherentemente inestable.

El análisis de eigenvalores es el método adoptado en este trabajo.

2.5.2 Barrido de frecuencia [19]

Esta técnica puede ser particularmente efectiva en el estudio del efecto de generador

de inducción y consiste en calcular la resistencia y la reactancia equivalente de la red vista

por el devanado de estator de un generador para diferentes valores de frecuencia. El método

da información acerca de las frecuencias naturales del sistema y la tendencia que existe

hacia la auto excitación y los problemas de resonancia subsíncrona (SSR por sus siglas en

ingles). Esta técnica es de mucha ayuda en análisis preliminares de SSR.

2.5.3 Análisis de respuesta en la frecuencia [19]

La estabilidad de oscilaciones subsíncronas puede ser estudiada mediante el uso de

diagramas de Nyquist multidimensionales. Este tipo de análisis puede manejar modelos

detallados y sistemas de gran tamaño.

2.5.4 Análisis aproximado en el dominio de la frecuencia [19]

Analiza la estabilidad individual de cada modo y está limitado a una detección

aproximada de SSR.

2.5.5 Análisis en el dominio del tiempo

Para este tipo de análisis se utilizan programas en el dominio del tiempo para

calcular la respuesta transitoria, como el programa para transitorios electromagnéticos

EMTP. En este caso se tiene una representación trifásica del sistema eléctrico, lo que
permite un modelo detallado del equipo, incluyendo los efectos no lineales. Esto facilita el

análisis de pares transitorios debidos a SSR, sin embargo el uso de estas técnicas de análisis

están restringidas a sistemas pequeños.

2.6 M O D E L A D O DEL SISTEMA ELÉCTRICO

2.6.1 El modelo de la máquina

En esta sección se desarrollan las ecuaciones para una máquina de rotor sólido. Este

modelo es utilizado para el estudio del fenómeno de auto inducción e interacciones

torsionales.

La Figura 2.3 muestra un modelo para una máquina síncrona con tres devanados en

el estator y tres devanados de rotor. En donde fü es el devanado de campo del rotor, Kd el

devanado amortiguador en el eje de directo y KqeI devanado en el eje de cuadratura. Las

ecuaciones de voltaje que representan a la máquina síncrona como un generador están

dadas por las ecuaciones (2.11) y (2.12):

>

Figura 2.3.-Máquina Síncrona

^abe ^s^abe +
^ ^abe (2.11) V ^ V + ^ v , (2.12)

Los encadenamientos de flujo representados por AgOC y Aqdr en las ecuaciones (2.11)
y (2.12) están dados por la expresión:

¿ahe L, l
abc
(2.13)
L, _ V

En donde:

Ls = Matriz de inductancias propias del estator.

Lsr = Matriz de inductancias mutuas entre el rotor y el estator.
Lr = Matriz de inductancias propias de rotor.

Las matrices Ls, Lsr y Lr están definidas de la siguiente manera:

-¿•a c o s 2 ^ r
-I]
- L g eos 2er -Ls c o s 2 |
" i i
/ \

eos 2 ,cos^r -
k =
" 3, -Le eos 2 (er + n) (2.14)
-Lb

> eos 2
V R
3,
eos 2 [er + L
a -Lb c o s 2 |

L
¡kq
c
°s< Lsfdsener L5kasener
2,t 2n
L = - — | Lsfüsen | 0, - — | L^sen \dr (2.15)
3
2n 2n
Lskqcos\dr+~ | L¡fdsen\er+ — j Lskdsen\0r+ —

+
L/kq t-mkq
L = 0 l-lfd +
Lmfd fdkd (2.16)
o -fOkó L/kú
+

En donde :

N r
N. x2
LA - nnjrlat U =
XMor,a 2
v 2 y
Una descripción más detallada de las inductancias anteriores puede encontrarse en [16]. Si
se definen las inductancias magnetizantes como:

U =3 (l*-Lb) Lma=^(LA+LB)

Estas pueden utilizarse para referir los parámetros del rotor al estator:

[Link] 2
LSkQ
V
\

N,
,
s y
3 \t~mq (2.17) ~mkq
>v

N 2
,3)
~mq (2.20)

¡sfd = (2.18) fd] ( 1) (2.21)

N
s) l 3J
M* \2
t-skd ~ - I L 'md (2.19) {Nkû '2
(2.22)
KNs s .3
J

De la misma manera:

ÜL (2.23) (2.25)
h
=

¡ Y j

V
J = v
j (2.24)
KNJ

En donde el subíndice j corresponde a los subíndices q y d ( definidos para las

cantidades de rotor ) según sea el caso.

La ecuación para los encadenamientos de flujos con las ecuaciones de rotor referidas

a estator se pueden expresar como:

''abes -i sbes
(2.26)
1
qdr !<a L 'qdr

en donde:
 eos dr LmdsenOr Lmasendr
27C_
U - L„qcos[er-f Laseri Lmasen (2.27)
3 / 3 ,
2n 2n'
^ c o s f ^ l L^sen | 6r + — | i^se/7 I + —

L = 0 md (2.28)
0 ~md Likd +

Las ecuaciones de voltaje referidas a los devanados de estator se pueden expresar de

la manera siguiente:

rs + pLs pLsr
- /
3 bes
(2.29)
Kr r'r+pír qdr

En donde p denota una derivada.

Aplicando la transformación de Park a las ecuaciones de voltaje y corriente, se puede

escribir:

K*o = TVztc íqdo ~ 7~idtjc

En donde T y 7 " e s t á n definidas como:

f
2^ 2n
eos i? eos eos 9 +
3 ; V 3 COS0 seno 1
2n 2 k\ f 2n
7- = ! sene sen^ú sen e+ r- 1 - cos^fl - — sen i9 1
3 - T [ ~3~ 3 ) { 3
1 1 i „ 2?r
— sen f i9 + — 1
2 2 2 3 J l 3
(2.30) (2.31)

Realizando las operaciones indicadas se llega a los siguientes resultados:

Kdo = ~rivdo + P\do + (2.32)

Las ecuaciones para el rotor permanecen sin cambios y se puede escribir:

 gdr- r'qdr
ri. PK (2.33)

Para los encadenamientos de flujo tenemos:

_ TLJ-1 TLS
''•qdo -i qdo
(2.34)
1
Qdr _ | (Or L [Link]

Si se toma 6 en la transformación de Park como el ángulo de referencia del rotor

(0 = 0r), se obtiene una simplificación en las matrices para las inductancias mutuas y

propias tal que los valores para dichas expresiones dentro de las ecuaciones serán
constantes, por lo que se puede definir:

+ L
<r,q 0 0" L/nq 0 0 "
1
TL.T' = 0 L
ls + 0 (2.35) Kr = 0 i-md (2.36)
0 0 0 0 0

0 0

1(4 )(ry 0 0 (2.37)

0 0

En forma expandida:

(2.38)
K = -Li*qikq + L mq (-i q + ¡ kq ) (2.41)

= -l-Jd + Lma(-ia +/w) (2.39)

•4 = - 4 4 + L m Â-' ! a + 4 + L ) (2-42)
¿o = " V . (2.40)
Kd = - C J L + L ma (-¡ d + L + L ) (2.43)

Considerando que y = co0X las ecuaciones (2.38) -(2.43) se pueden expresar en

' 1x
función de las reactancias y si se iguala X - V en las ecuaciones (2.32) y (2.33) se
coo y

tiene:

V
ç = ~rsíq Wd + — PVq (2-44) Vkç = Wkq +
jrP^Kq (2.47)

. a>r 1
= -rjds'd - — rq t rra (2.45) v
m = rfdhd +
Py/fd (2.48)
Û)„ CÜ„ CÜ„

^o = 'Ve + — PV0 (2.46) V

Kd = rM>kd + P^Kd (2.49)
0i„ m
¥q = -XJq+Xmq(-iq+ikq) (2.50) ¥kq = - X ^ + Xma(-iq + ikq) (2.53)

¥d = ~X,sia + Xma(-ia + +4 ) (2.51) v4 = - X ' J f d + X m q (-/ d + 4 + 4 ) (2.54)

^ =-V 0 (2.52) = -X,JM + + 4 + '«) (2-55)

De la misma manera se puede definir el conjunto de reactancias siguiente:

Xq - Xls + Xmq (2.56) Xfd - Xm + Xmd (2.59)

xü=xls + xma (2.57)

xka - Xlkd + Xmd (2.60)

Xfq - Xlkq + Xmq (2.58)

En donde:

xmd(xd-xls xmq [ xq - x„
Xifd - (2.61) Xíkg - (2.63)
xmd-(xd-xl5 xmq-(xq-x,

xd -X,£)XmaXlfd
Xfra - r (2.62)
XifdXm(J - ( Xü - X (s W Xmd + xlfd

La resistencia para los diferentes devanados, se puede calcular a partir de las

constantes de tiempo de la máquina de la siguiente manera;

Xikq + Xmq (2.64) X md X iftf

T
(2.66)
<°o qo ÚJ„T.
í
o do Xm(¡ + Xlfd j
X,fd + Xmd (2.65)
®oTdo

Sustituyendo las expresiones para los encadenamientos de flujo dadas por las ec's

(2.50)-(2.55) en las ees (2.44) - (2.49) se llega a las expresiones:

v.q -
<'a
s'q + ™ i-Xdid +x j f d + XmJkd ) — p (~Xqiq + Xmqi'kcl (2.67)
0K ÓJ„

'rJd - (~Xai0 + xmJk0 ) + — p (~Xüid + x j f d +x j k d ) (2.68)

co. íü„

v. = -rJc P(-X¿0) (2.69)

^ = rkqikq +—p{XJka
' kg kq- Xmamq
i, q v
*¿ = rJxd +
—P{Xküi\a ~ Xmdid + Xmdifd
OK

(2.70) (2.72)

Vf* = + — P (*«'« " + -^Â*

(2.71)
En forma compacta se puede escribir:

v = Lpi + Go)r¡ + Ri (2.73)

En donde:

-X* 0 Xmd Xmü 0 0 Xq 0 0 - X mq

0 0 0 Xmq -X* 0
i Xmd Xmd 0
L = — ~Xmd 0 X* Xmd 0 (2.74) G = — 0 0 0 0 0 (2.75)
co„
~Xmd 0 Xma Xkd 0 0 0 0 0 0
0 ~Xmq 0 0 * * 0 0 0 0 0

~rs 0 0 0 0 v vq 0 o
0 -r s 0 0 0
R= 0 0 r
fd 0 0 (2.76)
0 0 0 0 L /fd 'ka •kq
0 0 0 0

En las expresiones anteriores se han omitido los superindices en las variables para las
cantidades de rotor y se da por entendido que todas las cantidades están referidas al
estator. En lo sucesivo se omitirán todos los superindices cuando se haga referencia a las
cantidades de rotor.
2.6.2 El modelo de la red

/ v-Y Y \

Figura 2 . 4 . - D i a g r a m a Simplificado d e una Línea con C o m p e n s a c i ó n S e r i e

El circuito mostrado en la Figura 2.4 se puede representar por el conjunto de

ecuaciones siguientes:

1/(0-1/ {t) = Rl{t) + Lpl(t) (2.77) p(|/e(t) -É(0)= (2.78)

Asumiendo que no existen corrientes de secuencia cero , la transformación de Park

puede ser reducida a:

l{t) = Iocos{9) + Iasen{9) V(t) = Vgcos(d) + Vdsen(0)

Sustituyendo estas expresiones en las ecuación (2.77) se obtienen las expresiones de

caída de voltaje en la reactancia inductiva y la resistencia en sus componentes qd.

K - - (2.79) (vd - vcd) = ríü +—XLp¡d (2.80)

¿y„ ÉÜ„ co,.

Similarmente para el elemento capacitivo se tiene:

±p{vcd - Ed) = XJa +^(vcg-Eq) 4

(2.81) ±p{vcq -Eq) =X i L
- Ed
co, CO ' ' ®„ 0)„
(2.82)

Las ecuaciones (2.79) - (2.82) se pueden expresar de forma matricial.

- ^CD rL 0 1 XL 0 P'a 0 -XL

+ — - £
(2.83)
0 r, Cù~ 0 XL Piq co„ X, 0

~p(vcd-Ed) X o" V (2.84)

ry.,
2.6.3 Ei sistema de excitación

+
x
1 + sTc

V,

SKF
1 + ST-

Figura 2.5.-Diagrama de Bloques para el Sistema de Excitación

Con base al diagrama de bloques mostrado en la Figura 2.5, se plantean las

ecuaciones para el sistema de excitación:

SEfd=—[KE{Vs-VF~Vc)~Efa] (2.85) M i Eh-^tV,

1
E
TETF 1
B1 F r.F
(2.86)

En forma matricial:

KE ' ' Ke '

Te
"

-fa TE rE Te
\ ~
Efü
Kt + l/. (2.87)
Kc ±+k£Kf kekf kekf
TJr TJf y T
[ eTF . [ T£Tf J
2.6.4 El estabilizador

Las oscilaciones de baja frecuencia del rotor se encuentran en el rango de 0.2 a 2.0

Hz y se presentan debido a la oscilación de los rotores como resultado de algún disturbio. En

las líneas de transmisión pueden apreciarse en forma de oscilaciones de potencia que

pueden limitar la transferencia de potencia y amenazar la seguridad del sistema eléctrico.

Una solución rentable de proporcionar amortiguamiento a estas oscilaciones es el

Estabilizador del Sistema de Potencia ( PSS por sus siglas en ingles ). El PSS tiene una señal

de entrada (U) que se deriva de la velocidad del rotor, frecuencia del nodo, potencia

eléctrica o la combinación de estas variables. Su salida ( O P ) s e suma a la entrada del AVR

con el propósito de proporcionar amortiguamiento adicional sin que se afecte el par de

sincronización de las frecuencias de oscilación críticas.

 sTa 1 + S7; 1 + sT3 OP
1 + sT^ 1 + ST2 TT7T7
(a)

l
i+sT tó - v V > -í:
TTÍT: '

KS TI T3

(b)
Figura 2.6.- (a) Diagrama de Bloques para el Estabilizador; (b) Diagrama Equivalente

La Figura 2.6 (a) muestra un diagrama de bloques para el estabilizador. Un análisis

matemático del diagrama equivalente presentado en la Figura 2.6 (b) conduce a las
siguientes ecuaciones.

sX2 = i - l xi + 1-ZL KU - X,
(2.88) 7, vn y 7
1 7. 7.
(2.89)

TJ3
T2T4
rt
T2J
xl + i-L. X, ZL-ZL KrU - X, (2.90)
v VT2 TAJ

En forma matricial:

K,
0 0
7, 7
1 1
0 X2 1/ (2.91)
vT-2 7 T2 7,2
X,
HZL 7 r i 1 Ks73
M
v W l 4T
T4 74 y
[Link] El filtro torsional

Se ha observado que las vibraciones torsionales pueden presentarse en unidades

térmicas de gran tamaño que utilizan señales de estabilización basadas en la medición de la

velocidad de la flecha, por lo que se ha propuesto el uso de filtros torsionales diseñados

apropiadamente para este propósito. Al mismo tiempo el diseño del filtro debe asegurar un

margen de seguridad suficiente para interacciones torsionales, y un cambio de fase mínimo

a bajas frecuencias para un comportamiento óptimo del estabilizador.

El modelo en espacio de estado para el filtro torsional puede obtenerse a partir del diagrama

de bloques presentado en la Figura 2.7 y está dado por las ecuaciones:

s X = (OP - 2frnX - conV5) (2.92) - (2.93)

En forma matricial:

X <¡>n °>n OP
0 (2.94)
y. co_ o

OP
s + 2Cans

(a)

OP a>n A>N
>r
V Y
A. s S

2C <

(b)

F i g u r a 2 . 7 . - ( a ) D i a g r a m a d e B l o q u e s p a r a el Filtro T o r s i o n a l ; ( b ) D i a g r a m a Equivalente
2.7 MODELADO PARA EL SISTEMA MECÁNICO

Los rotores de las unidades de generación térmica son unidades muy complejas que
contienen secciones de flecha maquinadas de varios tamaños, así como diversos
acoplamientos. Las diferentes secciones de turbina contienen discos, alabes y otros
componentes adicionales. El generador incluye ranuras de bobina y anillos deslizantes, por
lo que este tipo de rotores tienen un número de vibraciones torsionales muy elevado. Estas
vibraciones pueden estar por arriba o debajo de la frecuencia del sistema, y para algunos
estudios se puede tener la necesidad de modelar el sistema de la flecha con bastante
precisión. Sin embargo, en el estudio de interacciones torsionales se puede utilizar un
modelo concentrado como el que se muestra en la Figura 2.8.

AP "" PM ' BPA < BPB GEN eXC

— . x y y y .y

Figura 2 . 8 . - M o d e l o de Flecha para Estudios de Resonancia Subsíncrona

La Figura 2.9 muestra el modelo masa-resorte correspondiente a la representación de

seis segmentos de flecha mostrada en la Figura 2.8. El modelo presentado considera cuatro
segmentos para representación de la turbina, dos de baja presión, uno de presión media y
uno mas que representa a la sección de alta presión, el segmento que representa al rotor
del generador y el excitador.

'l2 — D23 D34 - D45

-y- * M ¿, -VVY- M,3 -YVV M.4 "Y YY' Me

:> v v

"12
K 23 K 34 ___ K 45 K 56
D
i - D3 D4 " D5 - D6

F i g u r a 2 . 9 . - M o d e l o M a s a - R e s o r t e p a r a la F l e c h a d e l G r u p o Turbina-Generador

El modelo masa resorte para la flecha del generador mostrado en la figura 2.9 supone
que la flecha está formada por varios segmentos conectados entre sí por flechas de rigidez
finita representadas por la constante de resorte K . Se consideran el amortiguamiento D
entre secciones adyacentes de flecha y cada una de las secciones a la referencia. La
ecuación de movimiento para las n secciones de flecha representadas en la Figura 2.9 es:
 H9 + DO + K$ = T (2.95)
(o,
o

En donde 0 representa el vector angular de desplazamiento o torcedura (twists) de

cada una de las secciones de flecha con respecto a un punto de referencia común en grados

eléctricos. Por lo que 0 es positivo en la dirección de rotación de la flecha. 6 se mide con

respecto a la primera sección de flecha por lo tanto = 0 en la ecuación (2.95). La

frecuencia angular base en rad/seg es co0 .T es el vector de par externo aplicado a cada una

de las diferentes secciones de flecha. H es una matriz diagonal, es decir

H = diag(H1, H2, H3,...Hn), cuyos elementos son las constantes de inercia de cada una de las

secciones de flecha. Las constantes de inercia de las diferentes masas están dadas con

respecto a la misma base de volt-amp del generador,

Utilizando el ángulo de desplazamiento de la primera sección de flecha como

referencia para el desplazamiento angular de las demás secciones de flecha, la matriz de
rigidez K, para el conjunto de masas, es una matriz tridiagonal cuya forma estandarizada se
muestra en (2.96), en donde k0es la constante de resorte de la porción de flecha conectada

entre las secciones i y j en P.U. de par/rad.

K = (2.96)

D (2.97)
 La matriz D está compuesta de los coeficientes de amortiguamiento, los cuales
representan el amortiguamiento torsional entre masas. Similar a la matriz K, la matriz D
también resulta en una matriz tridiagonal como la que se muestra en:

En donde d10es el amortiguamiento entre la i-ésima sección de flecha y la referencia

rotatoria síncrona, y d¡j es el amortiguamiento entre al i-ésima y la j-ésima sección de

flecha. En la práctica los coeficientes de amortiguamiento son difíciles de determinar con

precisión y sus valores se determinan normalmente en el dominio modal ( mediante pruebas
de frecuencia del generador ) en donde las componentes de amortiguamiento modal son
más fáciles de determinar.

La forma expandida para las ecuaciones, en el caso de seis secciones de flecha, se

muestran a continuación:

M I P K ) - -(D12 + D U H + Dl2a2 - K12S, + K12S2 (2.98)

— pS1 =Qlí- 1 (2.99)

M2p((o2) = -(D]1+D23+D22)ú)2+D]2ü>t+Di3a>3+Kl2ói-(Kl2 + K2i)S2 + K11S¡ (2.100)

— pS2=o,-l (2.101)

M^p(^) = ^{D 23 +D M +D ü )m i +D a ío 1 +D^a) A +K 23 S 2 -(K 23 +K ii )S 3 +K 3i ó A (2.102)

— pS3=co3-l (2.103)

K P M - Í O ^ + D , , + +¿V>3-(K¡A+KI5)Ó4 +K„SI+K4,S5 (2.104)

-pó4=co4-1 (2.105)

+ + + (2.106)

— PS5=6)5-1 (2.107)

Mf¡P{ms) = -{Dw¡+D66)a>, + Dxa>s + -Kib5b (2.108)

— pó6 =<»6 - 1 (2.109)

 CAPÍTULO III

EVALUACIÓN DEL EFECTO DE GENERADOR DE

INDUCCIÓN Y LA RESONANCIA SUBSÍNCRONA

3.1 INTRODUCCIÓN

En el capítulo 2 se dieron algunas definiciones importantes relacionadas con la

resonancia subsíncrona, se derivaron las ecuaciones, así como la representación en espacio
de estado para los diferentes elementos involucrados en el estudio. En este capítulo se
linealizan los modelos matemáticos descritos anteriormente, se presenta la matriz de
estados utilizada para el estudio del efecto de generador de inducción, el cual es estudiado
mediante un ejemplo numérico. Los valores característicos calculados de la matriz de
estados del sistema son presentados para la comprobación de los resultados obtenidos
analíticamente. Se muestra el límite de estabilidad para este fenómeno y se analiza el
comportamiento del amortiguamiento del modo subsíncrono de la red ante la variación de
las resistencias de rotor de la máquina.

El análisis de interacciones torsionales se realiza observando el comportamiento de

los diferentes modos naturales de oscilación de la flecha cuando la frecuencia del modo
subsíncrono de la red se mueve a lo largo del rango en donde se encuentran ubicadas las
frecuencias de estos modos. Se presenta el vector de variables de estado utilizado para este
estudio, la región de resonancia completa para el modo torsional de la flecha con frecuencia
de 160 rad/seg así como las diferentes regiones de inestabilidad de cada uno de los modos
para el sistema de la Figura 2.4 cuyos datos se presentan en el apendice A.

Al final del capítulo se analiza la influencia del estabilizador de la máquina sobre los
modos torsionales. Se demuestra que este problema puede solucionarse mediante la
selección adecuada de un filtro pasabajas.
3.2 FORMACIÓN DE LA MATRIZ DE ESTADOS DEL SISTEMA

Los modelos en espacio de estado para los diferentes elementos del sistema eléctrico

y mecánico considerados en el estudio de resonancia subsíncrona pueden ser puestos en

forma incremental suponiendo que:

X = Xo+AX (3.1)

U = U0+MJ (3.2)

La suposición anterior permite la linealizacion del modelo de tal forma que para el

análisis de pequeñas señales se tiene:

Ax = AAx + BAu (3.3)

3.2.1 El modelo incremental de la red

De acuerdo con la suposición (3.1) y (3.2), las expresiones (2.67) y (2.68) se

pueden expresar en forma incremental de la siguiente manera:

Al/, = rLAÍq + — XLSAiq + XLAid + XLIÜOA0) + AV — SAVCQ = XcAiq - AVcdAv + EQOAS - AVCÜ

(3.4)
(3.6)

Al/fl = rLA¡a + ^XLSAia + XLAÍq + XJ^a + A^ J ^ ^ = ^ + ^ ^ + ^ + ^

Ú)n
(3.5)
(3.7)

AVq = -rsAiq - XaAid - XaIü0AÚJ + XmdAifd + XmaIfaoAco + XmdAikd - — XqSA¡q + -~XmqStikq (3.8)

Al/, = -r s A/ d + XqAiq + XqIqoAv-XmqAikq ~—XdsAia + — XmdsAifd + -LXmdSAikd (3.9)

®0

Igualando (3.4) con (3.8) y (3.5) con (3.9) se obtiene el modelo linealizado para la

red de transmisión incluyendo los efectos del estator de la máquina, (ec's (3.10) y (3.11)).

— (X, + X t ) s A i a + — XmdsAifd + — XmdsAikd =(rL +rs)Aid -ÍXL +Xq)Ai + XmqAika

v
o>~ * (3.10)
-(XL+XqIqo)Aa> + AVa
 .±(Xq+XL)s&iq +^-XmdsAikq =(rL +rs)A/, +(XL + Xd)Aid -XmdAifd -Xm0Aikd
^ (3.11)
+ [(XL+Xd)ldo-XmdIFd0]Aa> + AVi

Las ecuaciones (3.8) y (3.9) representan los efectos de los devanados de estator de

la máquina.

3.2.2 Efectos del rotor de la máquina

Los efectos para el rotor de la máquina pueden incluirse al considerar las expresiones

de voltaje para cada uno de estos devanados. Estas ecuaciones se obtienen de (2.67) a

(2,69). Las expresiones linealizadas para estos voltajes están dadas en las ecuaciones

siguientes:

• — XmdsAid+—XfdsA!fd + — XmásMkü = -rfdAifd (3.12)

a
CO o <¿o o

— XmdsAid +—XmdsAifd +-^XkdsAJkd ~ -rkdAikd (3.13)

^o °>0

~Xms&iq + — XkqsMkq = -rkqAikq (3.14)

(,J OJ
0 0

3.2.3 El acoplamiento entre el sistema eléctrico y el sistema mecánico

La ecuación de oscilación proporciona el acoplamiento entre el sistema eléctrico y el

sistema mecánico. Esta ecuación puede representarse en forma incremental de la siguiente

manera:

MsA0 = ATm-ATe (3,15)

— SAS = A® (3.16)

El par electromagnético en el marco de referencia qdo puede expresarse como;

T
e = iQ ~ Vq 'a
Sustituyendo y á y dados en las expresiones (2.43) y (2.44) se obtiene:

ATe = (Xq - Xd)lqoAld+[(Xq - Xd)ld0 + XmdIfüo]Mq + XmJqoAIfa + XmdIqoAlkd - XmqIdoAlkq (3.17)

Finalmente si se sustituye (3,17) en (3.15) se obtiene la expresión incremental para

la velocidad.

Ms&o) = Tm- {(Xd - Xq)lqoMd + [(Xd - Xq)ldo - Xmdlfd] Aig - XmdIqoAifd - XmdIqoAikd + XmgId0Aikq)

(3.18)

— sAS = Acó (3.19)

Finalmente, el conjunto de ecuaciones puede ser puesto en la forma:

Ex = Fx + Hu (3.20)

En donde E, F y H son matrices. Si se premultiplica ambos lados de la ecuación

(3.20) por E 1, se obtiene la representación en espacio de estado. Para el estudio del efecto
de generador de inducción, el modelo de la flecha del grupo turbina-generador puede ser
excluido, y las matrices mencionadas anteriormente tienen la forma mostrada en (3.22) a
(3.26). Colocando (3.20) en la misma forma que (3.3).

X = E^Fx + E'lHU = Ax + Bu (3.21)

En donde:

A = E 'F

B = E-'H

0 Xma Xmd 0 0 0 0 0
0 0 0 Xmq 0 0 0 0
0 Xmd 0 0 0 0 0
0 Xma 0 0 0 0 0
0 ~Xmq 0 0 0 0 0 0 (3.22)
Ú3„ * *
0 0 0 0 0 co0M 0 0 0
0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0 1
 (0+1) 0 0 Xmq 0 1 0
{XL + Xq)lao
(r,+rL) ~X<T>d ~Xma 0 0 0 1
~Xma Ifdo
0 0 ~rfd 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
[xü-xq)id0
~XmdIqo ~XmdIqo X mqfdo 0 0 0 0
y ¡vdT fdo
0 0 0 0 0 1 0 0 0
0 0 0 0 0 1
Edo
0 0 0 0 -1 0

0 0 ' fd
0 0 0 0 0 0
Hf = X m<J (3.24)
0 0 0 0 0 1 0 0 0

xc =[Md A iq A iFd A ikd A ikq A cog ASg AVcd AL/,,] (3.25)

u = ^fc (3.26)
A R„

3.3 A N A L I S I S DE O S C I L A C I O N E S A U T O E X C I T A D A S

3.3.1 Evaluación del efecto de generador de inducción

Como ya se ha mencionado, la compensación serie de líneas de transmisión provoca

la aparición de corrientes circulantes en la red de frecuencia natural (f n ). Esta frecuencia

está dada por la expresión (2,1) y reescrita en (3,27).

=ü)
oJy- (3.27)

En donde:
Xc = reactancia capacitiva

Xe = X + x7 + XL

y-

El circuito equivalente para cantidades subsícronas de la máquina en los ejes dq se

muestra en la Figura 3.1 [3].

f >« f í -y f
" £
A
x
>"i £ x» L x f ?«
mtf f Y0 y,
v..
[*2_
S

(a) (b)

F i g u r a 3 . 1 . - ( a ) C i r c u i t o E q u i v a l e n t e P a r a el E j e d ; ( b ) C i r c u i t o E q u i v a l e n t e P a r a el E j e q

En la ecuación (3,27) Xe >Xe, por lo que con < co0. Esto provoca que la máquina

síncrona, tenga un comportamiento similar al de un generador de inducción inyectando

energía de secuencia positiva a una frecuencia menor que la del sistema. La cantidad de

energía disipada en el estator dependerá del valor de la resistencia del circuito del rotor y

del deslizamiento. En otras palabras, la cantidad de energía absorbida dependerá de la

resistencia del estator y del circuito externo, ya que dependiendo de la resistencia total del

circuito, se puede presentar o no una situación de auto excitación en el sistema eléctrico.

La modulación en la frecuencia del sistema se hace más grande a medida que el

grado de compensación aumenta. El deslizamiento se hace más grande y la energía
generada inductivamente aumenta. El resultado es una disminución en el amortiguamiento
del modo de oscilación subsíncrono de la red, el cual puede llegar a ser negativo originando
así una condición inestable en el sistema.
Ejemplo

Un sistema máquina-bus infinito se encuentra compensado en serie al 50% de la

línea. Si la impedancia de la línea de transmisión es Z = 0.004 + j0.56 p . u c a l c u l a r la

frecuencia subsíncrona de la red como consecuencia de la compensación serie y los circuitos

equivalentes para la máquina en los ejes q y d. Determine si el sistema es estable para esta

condición de operación. Utilice los parámetros mostrados en la tabla 3.1,

Tabla 3.1.- Parámetros de la Máquina Tabla 3.2.- Valores Caráeteristieos para el

Análisis de Generador de
X te =0.130 X'q = 0.228 V - 0.850 Inducción

X d = 1.790 Xd =0.135 ^ - 0.032

X, =1.710 X¡ =0.200 <„=0.050 MODQ 1

: .1.
X" = 0.169 -1.47 ± 589.791 -1.49 ± 610.10Í MR +
4 - 4.300
-0.12 ± 164.12Í 0.15 ± 143.77i MR"
-1.26 ± 9.85Í -1.38 ± 10.30i ME
-51.03 -50.86

-1.40 -1.39

-1.00 -1.00

-5.53 + 19.79Í -5.75 ± 19.43Í MEXC

SOLUCION:

Primeramente se encuentran las reactancias y resistencias correspondientes a los ejes

qd. De las ecuaciones (2.56) a (2.66) se tiene:

mq" X-X.
q s ¡s = 1 . 7 1 0 - 0 . 1 3 0 = 1.580~p.u

xmd = Xd-Xls = 1 . 7 9 0 - 0 . 1 3 0 = 1.660 p.u

1
MR+ Modo de la red supersíncrono, MR- Modo de la red subsíncrono, ME Modo
electromecánico, MEXC Modo del excitador
 _ ~ x is) _ 1.660(0,169-0.130) _ Q4
= = = P
Xfflrf 1.660-(0.169-0.130) ' "

~ X„ XmaXIFd

x t f d x m ú ~{x a - x , A ( x m a + x l f d )

(0.135-0.130)(1.660) (0.04)
= 0.0057 p.u
(0.04)(1.660) - (0,135 - 0.130)(1.660 + 0.04)

X mq (X q - 1.580(0.228-0.130)
Xlkc, - = 0.104 p.u
Xmq-(x'q-xls 1.580-(0.228-0.130)

r>i = ( X „ + X M ) __ 0 - 0 4 , 1 . 6 6 0 = a 0 0 1
CO~T
o do 1
377(4.300)

XmgX,fd 1.660(0.04)
rM = X ikd 0.0057 + = 0.0037 p.u
XMR¡ + XLFÜ J 377(0.032) 1.660 + 0.04

Q
^ = ZjSLlXm. = : 1 0 4 + 1 - 5 8 0 = 0.0053 p.u
'*<7
377(0,850)

La reactancia subtransitoria se puede encontrar de la siguiente manera:

x - 0.135 + 0.200 = 0 -u

La frecuencia natural:

0.28
= 214 rad /seg (3.28)
'xü +x c
+XL ,169 + 0.14 + 0.56

La componente subsíncrona de la frecuencia:

CÚ. = d0 - 03n = 3 7 7 - 2 1 4 = 163 rad ¡ seg (3.29)

0 bien,

fer=fo-fn = 6 0 - 3 4 = 26 Hz (3.30)

El deslizamiento de la máquina está dado por:

¿ - ^ 3 4 - 6 0
(3.31)
t 34

La Figura 3.2 muestra los circuitos equivalentes para la máquina en los ejes dq
respectivamente. La resistencia y reactancia equivalentes del rotor, se pueden calcular
como,

r d +rq
X = (3.32) r = (3.33)

la resistencia total del sistema

r d + rq 0.0036 + 0.0061
1to«, = r , + r L + = 0.0015 + 0 . 0 0 4 - = 0.000652 p.u

>0.0738 >0.0738

>0.0227 J >0-0032 >0.0590 «

j0.9423 >0.8969 ¿ "
-0.0013 > -0.0049 -0.0070 >

V
• '"YX
>0.0772 j0.1292

-0.0036 -0.0061

(a) (b)

Figura 3.2.- Circuitos equivalentes para el estudio del efecto de Generador de Inducción;
(a) Circuito Equivalente eje d; (b) Circuito Equivalente eje q
 Debido a que la resistencia total es positiva, el sistema resulta ser estable para esta

condición de compensación serie. Por otra parte si se incrementa en 10% la compensación

serie de la línea, la frecuencia subsíncrona de la red se encuentra en 143.77 rad/seg y la

resistencia total es de -0,000544 p.u, por lo que el sistema es inestable. La Tabla 3.2

muestra los valores característicos para cada una de las condiciones de compensación

descritas anteriormente.

Si se considera el estabilizador del sistema de potencia, la condición de inestabilidad

que se presenta al compensar el sistema al 60%, puede ser superada, tomando como señal
de retroalimentación la potencia de la máquina. La Tabla 3.3, muestra los eigenvalores para
el sistema compensado al 60%, considerando el estabilizador de la máquina y señal de
potencia activa para la estabilización del sistema.

Tabla 3.3.-Eigenvalores para el Sistema Compensado al 60%

Considerando el Estabilizador de la Máquina.

Como puede apreciarse en la Tabla 3.3, el efecto de generador de inducción puede

solucionarse, utilizando el estabilizador de la máquina, mediante la retroalimentación de una
señal de control adecuada. Los coeficientes utilizados para los bloques de atraso y adelanto
en el caso estudiado son de: T l = 0.05, T2= 0.01, K=0.12 y Tw= 10.

El límite de estabilidad del sistema por efecto de generador de inducción se muestra

en la Figura 3.3. Se puede observar como al Incrementar la compensación serie de la línea
 de transmisión, la resistencia externa debe ser más grande para poder compensar el
efecto de la resistencia negativa del rotor de tal forma que el sistema no sea autoexcitado.

0 45

0.05-
0 0.001 0 002 0.003 0 004 0.005 0 006 0.007 0.008
R e s i s t e n c i a e x t e r n a en p.u

Figura 3.3.- Límite de Estabilidad por Efecto de Generador de Inducción

El efecto de la variación de las resistencias del circuito de rotor sobre el

amortiguamiento del modo subsíncrono se muestra en la Figura 3.4. Para este caso se tomo
un valor para resistencia externa de 0.01 p.u con el sistema compensado al 80% de la
línea. Puede observarse que el efecto de ir incrementando gradualmente el valor en
cualquiera de las resistencias del rotor se ve reflejado directamente sobre el
amortiguamiento del modo subsíncrono de la red. Esto se debe al efecto negativo que tienen
las resistencias del rotor en el fenómeno de generador de inducción para el caso mostrado
en la Figura 3.4, la variación de rkges el caso más drástico.

0.02 0.01

0.018 0.009

0.016
0.Û0Ô
0.007
0.014

0.006
0.012
0.005
0.01
0.004
0.00«
0.003
0.006
0.002

0.094
0.001
0.002
-1 -0.8 -0.6 -0.4 -0.2 0 0,2 0.4
Parte rea! de! m o d o subsíncrono

Figura 3.4,- Efecto de la Variación de las Resistencias de Rotor sobre el

Amortiguamiento del Modo Subsíncrono
3.3.2 Análisis de interacciones torsionales.

Otro problema relacionado con la compensación serie de líneas de transmisión es la

interacción entre el sistema mecánico y el sistema eléctrico. Este fenómeno es conocido en
la literatura como interacciones torsionales y su estudio requiere el modelado del sistema
mecánico, y deben linealizarse las ecuaciones que representan a la flecha del grupo turbina-
generador.

[Link] Forma incremental para las ecuaciones del sistema mecánico

Aplicando (3.1) y (3.2) a las expresiones (2.98) a (2.109) se obtiene el conjunto de

ecuaciones linealizadas para la representación del sistema mecánico en el modelo del

sistema. Estas ecuaciones están representadas por las expresiones siguientes:

M^i&co^ = "(D12 + Dn)A<ox + D12AOJ2 - K12AS1 + K12AS2 (3.34)

— SA^ = A ^ (3.35)

M2S{ACO2) = -(C 12 + D23 + D22)ACO2 + D^Aco^ + D23Aco3 + Kí2ASx - (Kn + K23)A52 + K23AS3 (3.36)

~sAÓ2=Aa2 (3.37)

M3S(A&3) = ~(D23 + D3A + D32)Aa>3 + D23Aoj2 + D34AÚ)4 + K23AS2 - (K23 + K3A)A83 + Ki4A¿> A (3.38)

~SA$3=ACÚ3 (3.39)

MaS(ACÚ4) - (D34 + D4S + D44)Aa>4 + D45A<v5 + D3aA(O2 - (K34 + K 4 5 )A¿ 4 + K3AA53 + KaiA8s (3.40)

(3.41)

MsS {ACO5) = -(D 45 + D55 + D 5 5 )A« 5 + D45A<Y4 + D56ACÜ6 + K45AS4 - (K45 + K56)ASs + K56AS6 (3.42)

— SA£5 = A cos (3.43)

M6S(AG>6) = -(D 5 6 + D66)ACO6 + D56AÚ)5 + K56AS5 - K56AS6 (3.44)

— SA¿>, - Acó, (3.45)

 Las ecuaciones (3.34) a (3.45) pueden ser puestas en la forma (3.20) y

posteriormente convertirse a la forma (3.21) de tal manera que puedan ser incluidas en el

modelo. El vector de estados considerando las ecuaciones para la representación de la flecha

esta dado por:

AXT = [AÍ d Ai q Ai fd Mkd Mkq AcOs ASs AV^ AVcq Aú)í A3Í A ^ AS2 A®, A<?3
(3.46)
AÚ)4 A A Ú > 6 A¿ 6 ]

Tabla 3.4.- Valores Carácteristicos para el Análisis de Resonancia Subsíncrona

0.G5 p.u p.u

r
2
í ' MODO Xt Xc~2G% Xc=30% , ' ' Xc F 4 0 %

MR + -4.77 + 466.81i -4.82 ± 504.091 -4.85 + 532.69i -4.87 ± 556.80i |

MR - -4.37 + 287.08i -4.23 ± 249.7 Oi -4.15 ± 220.811 -4.50 + 197.63i

M5 -0.01 + 298.18i -0.01 ± 298.18i -0.01 ± 298.18i -0.01 + 298.18i

M4 -0.13 ± 203.04i -0.12 ± 203.08i -0.05 + 203.28i 0.44 + 202.06i

M3 -0.03 ± 160.63i -0.03 + 160.65i -0.02 + 160.68i -0.01 ± 160.741

M2 -0.04 ± 127.03i -0.04 ± 127.03i -0.04 + 127.04i -0.04 ± 127.05i

MI -0.04 + 99.16Í -0.04 ± 99.21Í -0.03 ± 99.27Í -0.03 ± 99.34Í

ME -0.33 ± 8.72Í -0.54 ± 9.06Í -0.59 ± 9.42i -0.65 ± 9.81Í

ÍO% m ^ }
U. w'v
-4.89 ± 578.04¡ -4.91 ± 597.23i -4.93 ± 614.89i -4.94 + 631.32i -4.95 + 646.76i

-3.86 ± 175.58i -4.37 ± 157.12i -3.44 ± 138.64i -3.41 + 121.921 -5.52 + 104.58Í

-0.01 ± 298.181 -0.01 ± 298.18i -0.01 + 298.18i -0.01 + 298.18i -0.01 ± 298.18i

-0.11 ± 202.72i -0.13 ± 202.82i -0.13 ± 202.86i -0.13 ± 202.88i -0.14 + 202.89i

0.06 ± 160.95Í 0.78 ± 159.93Í 0.01 ± 160.34i -0.02 ± 160.45i -0.03 + 160.49i

-0.04 ± 127.06i -0.03 ± 127.091 0.01 ± 127.19i 0.16 + 126.711 -0.03 ± 126.92Í

-0.03 ± 99.44Í -0.02 ± 99.59Í 0.01 ± 99.84Í 0.12 ± 100.43i 2.74 ± [Link]

-0.72 ± 10.23Í -0.79 ± 10.69Í -0.88 ± 11.20¡ -0.98 ± 11.781 -1.11 ± 12.45Í

2
MR+ Modo de la red supersíncrono, MR- Modo de la red subsíncrono, M I Modo torsional 1,
M2 Modo torsional 2, M3 Modo torsional 3, M4 Modo torsional 4, M5 Modo torsional 5, ME
Modo electromecánico
 La Tabla 3.4 muestra los valores característicos para diferentes grados de
compensación serie en el sistema de prueba de la Figura 2.4. El análisis detallado de estos
resultados muestra que a medida que la frecuencia del modo de la red disminuye como
resultado de la compensación serie de la línea, los modos de frecuencias cercanas son
desestabilizados. El amortiguamiento negativo del modo se hace más crítico cuanto más
próximas se encuentren estas frecuencias, mientras que el modo de la red incrementa su
amortiguamiento.

Se puede observar como el modo torsional de frecuencia 298 rad/seg no cambia su

amortiguamiento y se mantiene estable aun cuando la frecuencia del modo de la red se

encuentre muy próxima a él. Esto se debe a que su inercia modal es muy elevada.

La Figura 3.5, presenta una región de resonancia subsíncrona correspondiente al

modo torsional de frecuencia 160 rad/seg para el sistema de prueba utilizado, el cual

corresponde al primer modelo de prueba propuesto por el IEEE para estudios de resonancia

subsíncrona [14]. Puede apreciarse como para valores elevados en la resistencia del sistema

la región llega a cerrarse. En la Tabla 3.4, se puede observar que este modo es el más

afectado por la compensación de la línea.

En la misma Figura 3.5 se presenta el detalle para una resistencia externa máxima de
0.02 p.u, debido a que la resistencia para sistemas de potencia reales no es muy elevada.

0.«

0.4

c!. 0.35
c
«
03
I Región de resonare» subsbicrona0.4
iu a25
es
8 0.2 0.3
a

O.
° 0.15 0.2

0.1
0.1

0.05
0.01 0.02
O 0,1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0,9
Resistencia externa en p.u

Figura 3.5.- Región de Inestabilidad por Resonancia Subsíncrona

 Regiones similares a la mostrada en la Figura 3.5 pueden encontrarse para todos los
modos. La Figura 3.6 muestra las diferentes zonas de inestabilidad correspondientes a cada
uno de los modos torsionales para el primer modelo de prueba del IEEE.

0.45

g Región de inestabilidad modo 99 rad/seg

g Región de inestabilidad modo 127 rad/seg

g Región de inestabilidad modo 160 rad/seg

^ R e g i ó n de inestabilidad modo 203 rad/seg

E Estable

I Inestable

0.05
0 0.0020.0040.0060.008 0.01 0 . 0 1 2 0 . 0 1 4 0 . 0 1 6 0 . 0 1 8 0.02

Resistencia externa en p.u

Figura 3 . 6 . - Límites de Estabilidad por Efecto de G e n e r a d o r de

Inducción y Resonancia Subsíncrona.

3.4 El PSS Y SU INTERACCIÓN CON EL SISTEMA MECÁNICO

El PSS proporciona adelanto de fase a frecuencias que se encuentran por debajo del
modo torsional de más baja frecuencia. Esto puede resultar en posible desestabilización de
uno o más modos torsionales debido a una reducción en el margen de ganancia provocado
por el PSS. Este problema puede solucionarse con la instalación de un filtro típicamente
pasabaja de segundo orden, La función de transferencia para el filtro se presenta en la
sección 2,6.4.1 y se presenta nuevamente en (3.47)

FILTRO = (3.47)
ss + 2£cons + nCÜ\
2

£ y <¡)n son parámetros que deben ser elegidos por el diseñador. La selección de

requiere el conocimiento de la frecuencia del modo torsional de más baja frecuencia.

Para el caso estudiado se ha seleccionado £ = 0.5 y = 20 ; el diagrama de Bode

para la función de transferencia dada en (3.47) se muestra en la Figura 3.7, en donde puede

apreciarse como la frecuencia de corte se encuentra alrededor de los 10 rad/seg, la primera

 Tabla 3 . 5 . - V a l o r e s Carácteristicos para el Análisis del PSS

C o n PSS con señal Con . PSS . con

Con e x c i t a t í ó r Cori PSS con
de velocidad sin Seftaí d e potencia
MODO-. Sin P S S - / • Filtro torsionaí
Filtro torsionaí sin FHtro torsionaí

MR + -4.77 £ 466.811 -4.77 ± 466.81Í -4.77 ± 466.81i -4.77 ± 466.811

MR" -4.37 ± 287.07i -4.37 ± 287.07Í -4.37 ± 287.07Í -4.37 ± 287.07i

MI -0.01 ± 298.18i -0.01 ± 298.18Í -0.01 ± 298.181 -0.01 + 298.18¡

M2 -0.13 ± 203.041 -0.12 ± 203.03i -0.13 ± 203.04i -0.13 + 203.04i

M3 -0.03 ± 160.64Í -0.02 + 160.64i -0.03 ± 160.64i -0.03 + 160.64i

M4 -0.04 ± 127.03i -0.04 ± 127.03i -0.04 ± 127.03i -0.04 + 127.03i

M5 -0.04 ± 99.21Í 0.01 + 99.19i -0.04 ± 99.21Í -0.04 ± 9 9 . 2 Ü

ME 0.53 ± 9.70i -1.01 ± 9.70i -0.06 + 7.19Î -1.01 + 9.70i

EXCITADOR -7.55 ± 16.39Í -4.90 ± 16.48Í -9.23 ± 15.65Í -4.85 ± 15.961

-47.96 -45.67 -47.91 -48.35

-2.82 -2.81 -2.91 -2.87

-1.00 -1.00 -1.00 -1.00

-28.53 -23.31 -0.07

-16.29 -17.40 -7.97 ± 4.51i

-0.100 -0.10

FILTRO -2.43 ± 12.95i

 CAPÍTULO IV

EL COMPENSADOR ESTÁTICO DE VARS Y SU

INFLUENCIA SOBRE LOS MODOS TORSIONALES

4.1 INTRODUCCIÓN

La compensación paralelo en los sistemas eléctricos de potencia es utilizada

para e! control de voltaje en la red de transmisión. Estos equipos tienen un impacto

2 2
significativo sobre las pérdidas por I r e I X, así como en la característica de

transmisión de potencia, Lo anterior hace posible incrementar la capacidad del sistema

de transmisión y al mismo tiempo tener márgenes de estabilidad transitoria más

grandes. Estos efectos contribuyen al mejoramiento de las condiciones operativas del

sistema logrando así un sistema más rentable y de mayor eficiencia.

El compensador estático de VARs ( CEV ) es un elemento de sistema que

presenta características dinámicas muy interesantes y por lo tanto una alternativa para

la solución a los problemas de soporte de voltaje y estabilidad. Desde su aparición a

mediados de la década de los setenta, ha venido revolucionando la operación de los

sistemas eléctricos de todo el mundo y hoy en día es visto como la primera

generación de controladores FACTS.

En este capítulo se analiza el comportamiento del compensador estático de

VARs y la influencia que tiene sobre los modos de oscilación de la flecha del grupo

turbina-generador. Se desarrollan los modelos matemáticos utilizados para la red de

transmisión sin compensación, y con un CEV en el punto medio de la línea. El modelo

en espacio de estado para el CEV también es presentado. Al final del capítulo se

analiza el comportamiento del sistema para diferentes condiciones de operación y se

hace una comparación del comportamiento de los diferentes modos.

4.2 REGULACIÓN DE VOLTAJE EN EL P U N T O MEDIO D E LA

LÍNEA

La Figura 4.1 (a) muestra un sistema simple de dos máquinas en donde se

conecta un compensador estático de VARs en el punto medio de la línea de

transmisión. El compensador se representa con una fuente de voltaje d e la misma

amplitud que los voltajes de envío y recepción conectada en el punto m e d i o marcado

como Vm. El efecto de conectar un CEV en este punto es el de dividir la l í n e a en dos

X
secciones independientes de reactancia — cada uno.

La primera línea conecta los puntos de envío y medio, mientras q u e la segunda

enlaza el punto medio con el extremo receptor. Las relaciones entre VE, V m / VR y las

corrientes en cada sección de línea Iem e Imr se muestran en el diagrama f a s o r i a l de la

Figura 4.1 (b):

1 V £
2 2
(a) (b)

Figura 4.1.- (a) Sistema Simple de dos Máquinas Compensado en P a r a l e l o

(b) Diagrama Fasorial

Las expresiones para voltajes, corrientes y potencia del circuito de l a figura

4.1(a) se pueden derivar del diagrama fasorial de la Figura 4.1 (b),

Vem = Vfflr = l / c o s | (4.1) I em = I mr = I = (4.2)

En donde:
La potencia transmitida esta dada por,

s_
P
> = 'V
ArY>Iem
em i Jmr = Jm eos | = l/7 eos | (4.3)
4

(4.4)

Similarmente,

(4.5)

4.3 LA INFLUENCIA DE LA COMPENSACIÓN PARALELO SOBRE

LA CARACTERÍSTICA DE T R A N S M I S I Ó N DE POTENCIA

En la sección anterior se vio como la compensación paralelo puede incrementar

significativamente la potencia transferible máxima. Este aumento en la característica
de transferencia de potencia, permite al sistema de potencia incrementar su capacidad
de transmisión, al mismo tiempo que se consigue un margen de estabilidad transitoria
más amplio, por lo que el sistema puede mantener el sincronismo ante posibles
oscilaciones del sistema más severas, La Figura 4.2 muestra la curva de transferencia
de potencia P - S para el sistema de la Figura 4.1 (a).

P *

P„n
V2
PT

s
K
Figura 4.2.- Curva P-5
 En la Figura 4.2, el sistema sin compensación está representado por la curva A,

Mientras que el sistema compensado está caracterizado por la curva B. Es evidente

que la compensación paralelo en el punto medio de la línea incrementa la característica

de transmisión de potencia, así como el margen de estabilidad transitoria, que para el

caso del sistema no compensado está representada por el área A, mientras que para el

caso compensado está representada por A+B.

4.4 ANÁLISIS DE ESTABILIDAD PARA LOS MODOS

T O R S I O N A L E S DEL GRUPO T U R B I N A - G E N E R A D O R

4.4.1 El modelo de la red sin compensación

Antes de derivar las ecuaciones para el sistema de transmisión con

compensación paralelo, se derivarán las ecuaciones para la red de transmisión no

compensada, con la finalidad de poder hacer una comparación de los modos

torsionales entre ambos casos.

La Figura 4.3 muestra un circuito RL en serie, el cual representa la línea de

transmisión para el sistema máquina bus infinito.

V rL LL E

;" A \ o r ^ o o

Figura 4.3.-Línea de T r a n s m i s i ó n no C o m p e n s a d a

El circuito mostrado en la Figura 4.3, puede ser representado por la siguiente

ecuación:

V(t)-E{t) = rí{t) + Lpi(t) (4.6)

La ecuación anterior puede ser puesta en función de la reactancia y

representarse en el marco de referencia qdo de la siguiente manera:

Va-Eg= rLiq + —XLpiq + ~XJd (4.7) Vd-Ed = rjd + — XLpid - — XLiq (4.8)
^ CÜN C0-. CON

Las expresiones (4.7) y (4.8) pueden transformarse ai dominio de Laplace y ser

puestas en forma incremental mediante la sustitución de (3.1) y (3.2):

AVq = rLAi + — XLsAÍ + XLAíd + XLId0AOJ - Ed AS (4.9)

¿y„

AVd = rLAÍü + XLA¡q - XLIQ0AC0 + EqoAS (4.10)

Finalmente si se igualan (3.8) con (4.9) y (3.9) con (4.10) se obtiene el modelo

de la red con los efectos de estator incluidos.

— {Xd+X1+X2) sAi + — XmqsAikq = (rs +r1+r2) Aiq - (Xq + X1+X2) Aid - XmdAifd
a
°>o o
- + [(X, + + X 2 ) Ido - XmdIfü0 ] Aa> - Ed0AS

~—-(Xd + X1+X2)sAÍd + — XmdsAifd + -i-X dsAikd = (r + r, + r2) Aid

<O0 ' CÚ0 o)0 (4.12)
- (X, 4- X , + X 2 ) Aiq + XwqAÍkq - (Xq + X , + X 2 ) IqQAco + 5 q o A i

El vector de estados para el sistema sin compensación es:

A iq A id A ifd A ikd A ikq Acog A Sg] (4.13)

4.4.2 Modelado de la Línea con un CEV conectado en el punto medio

Anteriormente se mostró la influencia de la compensación paralelo sobre la

característica de transferencia de potencia. En esta sección se desarrollan las

ecuaciones utilizadas para la evaluación del efecto que tiene el compensador estático

de VARs con los modos torsionales del grupo turbina-generador. Para este propósito se

hará referencia a la figura siguiente:

•150048
 V r. Vm r, k
'Y'Y^y'A

Figura 4 . 4 . - S i s t e m a de T r a n s m i s i ó n con C o m p e n s a c i ó n Shunt

El análisis de! circuito mostrado en la Figura 4.4, permite plantear las

ecuaciones siguientes;

(4.14)

Vcev-0 = -B-% (4.16)

Ke, - E = r2i2 + L2pi2 (4.15)

Transformando las ecuaciones anteriores al marco de referencia qdo se tiene:

Vq - Kevq = rJq + kP', +

<i a (4.17) V^ - Erf = r2i2ú + L2pi2a - coL2i2q (4.20)

1
^ - Vce* = + kPi* - aL-Jo (4.18)
x
' Vcevg = -6"
^ /'lq (4.21)

(4.19) (4.22)

Aplicando LCK en el nodo marcado como 1/ :

>2 = (' + O

Sustituyendo la expresión anterior en (4.19) y (4.20)

K««, = r2 (>'„ - ilq) + L2p(iq + / l a ) + COL2 (id + / l t f ) + Eq (4.23)

 = r2 {'* + iu¡) + LzP(Í* + - WL, (iq + iia) + (4.24)

De (4.21) y (4.22)

'i, = C4-25) /tó - - B V ^ (4.26)

Sustituyendo (4.25) y (4.26) en (4,23) y (4.24) respectivamente se obtiene:

^ = r2 (iq - BV^) + L2p (iq - BV^q) - coL, (/, - BVcevd) + Eq (4.27)

^ = " BVcevd) + L2p{id _ BVceva)- vL, (iq - BVcevq) + Ed (4.28)

Si se sustituyen (4.27) y (4.28) en (4.17) y (4.18)

V* = + r2)iq+(L1+L7)piq+«>(L1+L^) - r2BKe.q ~ LzPBV^ + [Link]^ + Eq (4.29)

+ +
= (1 + O ' , (¿i " «(¿i + M ' , - ^BV^ - ¡-2pBVcevd + a ^ B l ^ + (4.30)

Del diagrama fasorial de la Figura 4.1 (b):

Eq - Eeos s

E„ = EsenS

Sustituyendo las expresiones anteriores para Eqy Ed en (4.29) y (4,30) se

obtienen las ecuaciones de la red con un compensador Estático de VARs en el punto

medio de la línea.

Vq-{ri+r2)iq+{L1+L2)piq+o,{L1+L2)id-r2BVcevg-L2pBVcevq+ü)L2BVcevd+EcosS (4.31)

K, = (r! + r2)ia +{L1+L¿)p¡d-(o{Lí+ L2)iq-r2BVce„d-L2pBVcevd+a>L2BVcevq + EsenS (4.32)

Las ecuaciones de la red con los efectos del estator de la máquina incluidos, en
función de las reactancias las ecuaciones (4.31) y (4.32) pueden expresarse en la
forma:
K, = (ri + r
2)'q + — + X2)piq + — {X1+X2)id-r2BVcevq -±x2pBVcevQ
(4.33)
+ ~X2BVcevü+EcosS

= ( r i + r2)id + — ( X 1 + X2)pid - ^-{X, + X2)iq - r2BVcevd -±X2pBVcevd

(4.34)
+ ^-X2BVcevq+EsenS
co„

Transformando las expresiones (4.33) y (4.34) al dominio de la frecuencia y

aplicando (3.1) se obtienen las expresiones linealizadas:

A ^ = (rt + r2)AiQ + [X, - X2)Aid + [{X, - X2)Ido + X2VcevüoB0]Acó - r2VcevqAB - EdA8

1 1 (4.35)
+ ±-{Xl+X2)sMq-fX2VcevqosAB

AV.^ir^r^Ai.-iX^X^A^+l-iX^X^I^+X.V^B^Aa^-r.V^AB^E^S
1 1 (4.36)
+ —(X1+X2)sAid-—X2VcevdosAB

Igualando (3.8) y (3.9) con (4.35) y (4.36) respectivamente se obtienen las

ecuaciones de la red con efectos de estator incluidos.

- — (xd + x, + X2) sAid + — XmdsAifd +1- X^sAi^ +-A- X2V^üosAB =

(r; + r 2 )A/ t í - ( X ^ +X2)A/q (4.37)

+ (X2Vcevq0 ~ r2VcñVd0)AB + EqoAó

-j^(Xq +X1+ X2)sAiq -i--^-XmqsAikq +^X2VcevqosAB =

(rs + rx + r2)AÍQ ~{Xa + X, - X2)Aid + XmdAifd ~ XmdAikd (4.38)

+ [ - ( X d + X , - X2)lqo + X2Vcevd0B0] Acó + [X2Vcev(!o ~ r2Vcevü0)AB - EdoAS

4.4.3 Obtención del modelo en espacio de estado del CEV

El modelo utilizado para la represtación del compensador estático de VARs se muestra

en la Figura 4.5

V «(i + sr,)
(1 + sT 2 ) (1 + sT)

(a)

v+ i T
>
Tz 1 + ST
lt ST}

(b)

Figura 4.5.- (a) Diagrama de Bloques para el CEV; (b) Diagrama de

Bloques Equivalente

La Figura 4.5 (a) muestra un diagrama de bloques para el compensador estático

de VARs, En la Figura 4.5 (b) se muestra el diagrama de bloques equivalente, de

donde se derivan las siguientes ecuaciones:

sX = V X
T2~T22J r T2T T2T T

(4.39) (4.40)

Utilizando (3.1) y (3.2) se pueden expresar las ecuaciones (4,39) y (4.40) de la

siguiente manera:

1 fC KT KT 1
SAX = Mref + i - i Al/ — — AX SAB = ^ A X + —¿-Al/, REF -~±AV-±AB
\TZ ' T
2 J T } ' T 2 J r, T TT T2T - T

(4.41) (4.42)

En forma matricial:

1 "1 TX- "Ti 1"

0
AX R2 "AX" T2 TI
+ T}'T2 A\/ (4.43)
AS K 1 _A B KT,
KTI
. T T. L Ti T T2T _
El vector de estados con el compensador estático de VARs incluido,

[M d tiq A/w Mkd Mkq Acag ASg AX AS] (4.44)

4.4.4 Análisis de valores característicos

La Tabla 4.1 muestra los eigenvalores para diferentes condiciones de operación del

sistema de prueba de la Figura 2.4.

Tabla 4.1.-Valores Característicos para el Análisis del CEV en las Interacciones Torsionales

MR -9.23 ± 376.96Í
ill§i
-9.23 ± 376.961 -9.11 ± 377.04i

M5 -0.01 ± 298.18Í -0.01 ± 298.18Í -0.01 ± 298.18i

M4 -0.13 ± 202.99Í -0.10 ± 202.99Í -0.13 ± 202.97Í

M3 -0.03 ± 160.62i -0.00 ± 160.62i -0.01 ± 160.61Í

M2 -0.04 ± 127.02i -0.03 ± 127.03Í -0.03 ± 127.02i

MI -0.04 ± 99.14Í 0.15 ± 99.20i 0.14 ± 99.20i

ME 0.73 ± 9.20i -1.01 ± 9.20i -1.01 ± 8.98i

EXCITADOR -7.68 ± 16. l O i -5.73 ± 16.85Í -5.85 ± 16.49Í

-47.72 -43.85 -43.88

-3.01 -3.00 -3.25

-1.00 -1.00 -1.00

-0.10 -0.10

-80.05 ± l l . S l i -87.15 ± [Link]

-20.00

I -80.00

1
MR Modo de la red, M i Modo torsional 1 M2 Modo torsional 2 M3 Modo torsional 4 M4
Modo torsional 4 M5 Modo torional 5 ME Modo electromecánico
La Tabla 4.1 muestra los valores característicos con diferentes elementos de control en
el sistema de prueba de la Figura 2.4. En general puede decirse que el compensador
estático de VARs tiene una acción amortiguadora para los modos torsionales. En el
caso particular del modo de frecuencia 298 rad/seg, el amortiguamiento se mantiene.
Esto se debe a que la sección de flecha correspondiente a este modo tiene una inercia
muy grande, por lo que no se ve afectado. El amortiguamiento negativo del modo de
frecuencia igual a 99 rad/seg, se presenta debido a que a la ausencia del filtro
torsional en el PSS
 CAPÍTULO V

EL TCSC Y SU INTERACCION CON EL SISTEMA

MECÁNICO

5.1 INTRODUCCIÓN

A lo largo de capítulos anteriores se ha evaluado la Interacción que existe entre el

sistema mecánico y el sistema eléctrico debido a la compensación serie, así como a la
presencia de algunos elementos c o m o es el caso del estabilizador del sistema de potencia.
Se ha visto como a medida que la compensación serie se incrementa, la componente
subsíncrona de la frecuencia del sistema eléctrico se va a acercando a cada uno de los
diferentes modos naturales de oscilación de la flecha, provocando una situación resonante
en el sistema.

La presencia de dispositivos FACTS en los sistemas eléctricos de potencia modernos

ha contribuido enormemente a la mejora en la operación de estos sistemas. El capacitor
serie controlado por tiristores (CSCT ó mas comúnmente TCSC por sus siglas en ingles),
tiene características de diseño, que permiten tener un control sobre el efecto de resonancia
subsíncrona en el sistema eléctrico, permitiendo con esto mejorar la transferencia de
potencia de una manera segura para el sistema.

Este capítulo presenta los conceptos básicos así como el desarrollo de ecuaciones
para la consideración del TCSC en estudios de resonancia subsíncrona. El TCSC es
considerado un elemento de ajuste rápido y eficiente a la reactancia de la red, el cual está
compuesto de un reactor controlado por tiristores "TCR" conectado en paralelo con un
capacitor fijo.
5.2 MODELO DEL TCSC EN ESTADO ESTACIONARIO

La obtención de las ecuaciones de voltaje y corriente en el TCSC se logra mediante la

solución de las ecuaciones diferenciales en el dominio de Laplace para los diferentes
periodos de conducción y no conducción de los tiristores. La Figura 5.1 muestra un circuito
LC paralelo, en donde los tiristores están representados mediante un interruptor, el cual
permanece abierto o cerrado según sea su condición de operación.

Si se supone que la corriente de línea es senoidal y se toma como referencia el cruce

por cero de la corriente para el disparo de los tiristores se tiene:

/(f) = Isen(cot) (5.1)

Asumiendo que el primer tiristor ( simbolizado por el Interruptor en la figura 5.1 ) se

dispara en el tiempo t = tx, se tiene

' ' C W c C M í O (5.2)

En donde:

ic (t) = corriente en el capacitor

iL (t) = Corriente en el inductor

Durante el primer periodo de conducción del tiristor, el voltaje en el inductor y el

capacitor son iguales, es decir:
 t ti

L
<=K
J

Figura 5.1.-Circuito Eléctrico para un TCSC

Aplicando LVK al circuito de la figura 5.1 se obtiene:

L
'dJL {t) =
~h l { t ) d t +
^ " ^ (5,3)

En donde ^(tjes el voltaje del capacitor al momento del encendido del primer

tiristor y U(t - t¡_) el operador de retardo.

Expresando las ecuaciones (5.2) y (5.3) en el dominio de Laplace:

í ( s ) = / c (s) + í t ( s ) (5.4)

sut(s) = ^ + ie-^(tl) (5.5)

Debido a que el periodo de conducción inicia a partir de t^se debe obtener J(s)a

partir de ^, por lo que 7(s) está dada por:

/(s) = i { / ( t ) U ( t - t O }

I (s) = Ie'st' [ s e n ( ^ í ) s + o) e o s (aiÍ! ) ] (5.6)

Despejando I c (s) de (5.5)

J c (s) = s a LCT t (s)-Ce-*.v;(t 1 ) (5.7)

Sustituyendo (5.7) y (5.6) en (5.4) y resolviendo para I L (s) tenemos:

IL (s) = a>20Isen « ) e"st> — ^ — + eos « ) e"st'

(5.8)

(S +6>t)

en donde:

1
Cú„ --
sÍLC

La expresión en el tiempo para la corriente en el inductor IL(t) se obtiene aplicando

la transformada inversa de Laplace a la ecuación (5.8); para cada término se tiene:

i 1
r1 -1 2 —Y e"s£l = — senfat -
s + coi co„

2^\cos{cot-a>tí)-co${co0t-co0tí)\
(s2 + « 2 )(s 2 +ü>l) — OJ

1 -2 (eot - mtJ - cosen (co0t - coDf,)]

(s2+«,2)(s2+íy2) j co0co(co2~co2)

Por lo tanto

IL (t) = A [sen (ffltj eos (co ( t - 0 ) + eos ( o t j sen (©(t-t,))] - ¿ s e n « ) eos (a>0 ( t - t j )
- B e o s « ) sen (a>0 (t -t^ + DV. (t^sen (a>0 (t - t,))

(5.9)
En donde:

CO. - (ù

CÚ0 -co2

D = co0C

Agrupando:

IL (t) - Asen (mt) - Asen (m^ ) eos (o>0 (t - t j ) - B eos (ait^ ) sen (co0 (t - ) )
(5.10)
+ D\/c(t1)5en(íyo(t-f1))

Agrupando (5.10) se obtiene:

i(t) = Asen(ú)t) + /q cos(w 0 t) + k2sen(co0t) (5.11)

En donde:

kx = B eos (a) sen (ka) - Asen (a) eos (ka)-DVc(tl)sen(ka)

K2 = DVC (tx) eos (ka) - Asen (a) sen (ka) - B eos (a) e o s (Ara)

(0

a = cotx - ángulo de disparo del tiristor

La expresión (5.11) es la ecuación general para la corriente que circula a través de

los tiristores. El voltaje en el capacitor está dado e n la ecuación (5.12)

Vc(t) = L±i(t) (5.12)

Sustituyendo la expresión para i(t) dada en la ecuación (5.11) se obtiene:

Vc (t) = AXL eos (<yt) + AkXLsen (a) sen (co0t - ka)- BkXL eos (a) eos (o0t - ka)
<5-13)
+ DkXLVc (tjcos^f - ka)

O bien en forma reducida,

Ve (t) = AXL eos(eat) - kX^senfat) + kXLK2 cos(co0t) (5.14)

En donde:

X L = col—

Las ecuaciones (5.11) y (5.14) son las ecuaciones generales para la corriente en los
tiristores y el voltaje en el capacitor en el primer periodo de conducción, esto es a partir del
cierre del interruptor mostrado en la figura en t = ^ 1.

Para el caso en el que los tiristores no están en conducción, se tiene que la corriente
del inductor es cero y la corriente del capacitor es igual a la corriente de línea, por lo tanto:

ic (t) = i(t) = Isen(a>t) (5.15)

Y el voltaje en el capacitor está dado por la expresión:

Vc (f) - ^ \lsen(a)t)dt + Vc (0) (5.16)

Co

En donde Vc (0) es la condición inicial para el voltaje en el capacitor. Si se supone un

tiempo de inicio en t = 0 y se resuelve (5.16) se obtiene:

1
Las ecuaciones ( 5 . 1 1 ) y ( 5 . 1 4 ) son válidas en el intervalo t e [ t j , t¿ ff ] . Para el s e g u n d o i n t e r v a l o de conducción

de los tiristores, las ecuaciones son las mismas, la diferencia está en el t i e m p o d e cierre d e l i n t e r r u p t o r , el cual para

el segundo intervalo de conducción será: t = t2 = tx + , por lo que el s e g u n d o intervalo se define para

t e [t2, J , en donde t¿ ff y son el t i e m p o de a p a g a d o del primer y s e g u n d o t i r i s t o r r e s p e c t i v a m e n t e .

 (5.17)
= -/Xe[cos(fflt)-l] + l/ c (0)

La expresión (5.17) se cumple para el período ¿yte[0,a]. Para el período cot e [cyt^,^ + a ]

se tiene:

CrJt

(5.18)

l/ c (tj f f )se puede obtener de la ecuación (5.14), sustituyendo t = t10ff. De manera

similar se obtiene la expresión para el intervalo de cot e [¿at^, 2;r] .

(5.19)

En donde Vc (t* ff )se calcula a partir de (5.14) sustituyendo t = t\ñy cambiando por t2.

Las ecuaciones anteriores determinan el comportamiento transitorio del TCSC durante

un ciclo. El dispositivo alcanza su estado estable cuando el pico máximo de la corriente de

línea coincide con el pico mínimo de la corriente que circula por los tiristores. Esta también

es simétrica con respecto al cruce por cero del voltaje en el capacitor. El voltaje en el

capacitor en estado estable puede calcularse al momento del disparo. Sustituyendo cot -

en la ecuación (5.14) se tiene:

(5.20)
y
 Despejando ^(í^ se obtiene el voltaje en el capacitor en estado estable al

momento del disparo:

l/ c (t 1 )=-|cos(a)tan (5.21)

Si se sustituye (5.21) en (5.11) se obtiene la corriente de los tiristores en estado estable.

/ x sen (a) í n
i (t) = Asen (eot) - A f ' ^ eos\co D t-^k (5.22)
eos

La ecuación (5.22) se cumple en el intervalo t +tc~], en donde ta es el tiempo

de conducción del primer tiristor en estado estable. El segundo pulso de disparo se da 180°

después del primer pulso de disparo, produciendo un flujo de corriente en el segundo tiristor

con la polaridad opuesta a la de (5.22). La corriente del segundo tiristor esta dada por:

sen(a)
i(t) = Asen(cot) + A ^cosl co0t ~~-k (5.23)
eos fclf-a

Sustituyendo las ecuaciones (5.22) y (5.23) en (5.12) se obtienen las expresiones

para el voltaje en el capacitor.

VC (t) - AXL eos (cot) + AkXL ^ sen | WJ- _ * k (5.24)

eos k\l-a
2

sen(a) i 3/r .
Ve(t) = AXL eos (cot)-AkXL seni co0t~~k (5.25)
K
cosi kI — - a

Despejando i c (t) de la ecuación (5.2) se tiene:

Vc ( t j s e encuentra en las expresiones para kx y k2 que aparecen en la expresión (5.20)

 (5.26)

Estas expresiones describen el comportamiento dei TCSC en estado estacionario en el

intervalo en que no existe conducción por parte de ninguno de los tiristores. Por otra parte,

el voltaje en estado estable cuando los tiristores se encuentran abiertos está dado por:

=¿ í ¡sen^dt + V^) (5.27)

2t-c¡

Donde l/ c (t o f f )es el voltaje del capacitor en el momento de la conmutación de los

tiristores. De (5.25) se tiene:

2 TI - A
Ve (toff) = Vc — I = AXl eos (a) - AkXLsen(a) tan (5.28)
a> v U „

Sustituyendo (5.28) en (5.27) y resolviendo para Vc(t)se obtiene:

K (f) = - ¿ C O S
(*0|L- + v
c for) - -IX c [cos(©0 - cos(a)] + Ve [toff) (5.29)

La expresión anterior se cumple en el intervalo cot <

= [0,«] y cot\2x - a, 2/r]. Para

cot e [re - a,7t + a] se tiene:

f
1 ,
=- \ lsen(vt)dt + Vc(toff) (5.30)

En donde l/c (f o i r )se calcula de (5.24) y está dado por:

!t - a f /
vc M = ^c [ — | = -AXL eos (a) + AkXLsen(a)tan (5.31)
v U
Y de (5.30) se tiene que:

= + ^ ( t o f f ) = - / X c [ c o s ( ® t ) + cos(a)] + l/(;(t0(f) (5.32)

La expresión anterior se cumple en el intervalo a>t <k\_n - a,n + a~\. Definiendo

V0 = Vc (toff) y comparando (5.28) con (5.31) puede observarse que:

K M = -ve(toff) = vt¡

Finalmente las expresiones para el voltaje en el capacitor para los periodos de

conducción en estado estable están dadas por:

Vc (T) - -IXc [COS(ÚÍÍ) - eos (A)] + V0 (5.33)

Vc (t) - -IXe [eos (a*) + eos (a)] -VQ (5.34)

La ecuación (5.33) es válida para los periodos en donde cot e [0,ar] y cot e [2n - a, 2tÍ\,

mientras que la ecuación (5.34) es válida en el intervalo cot € - a, n + a~\.

Tabla 5.1.- Corriente en el Inductor en Estado Estable Durante un Ciclo

it(t) =0
cot e [0, a)

, , i . sen (a) / \
,t (t) = Asen (cot) - y - eos co0t — k
^ 2 J cat e [a, n - o)
eos k I — - a 1
l U J,
it (t) = 0 cot - a, X + a)

=
't ( 0 Asen (tafr) - A cot e [/r + a, 2K - a)

it(t) =0 cot e \2n ~ a, 2-n)

 Las ecuaciones desarrolladas anteriormente describen el comportamiento del T C S C en

estado estable durante un ciclo de operación, y se muestran de manera resumida en las

Tablas 5.1 y 5.2.

Tabla 5.2.- Voltaje en el Capacitor en Estado Estable Durante un Ciclo

Ve )- [cos M - eos («)] + l/0 cot € [0, a)

sen(a)
Vc (t) = AXL eos (cot) + AkXL -sen[á>0t-jkj
cot e [«, n - a)
cos

Vc (t) = -ImXc [cos(fflt) + cos (a)] _ l/0 cot e \n - a, n + a)

sen (a)
Vc ( t ) = AXl cos (cot) AkXL sen^ai0t-~kj
cot e [tt + a, 2n - a)
cos
H f - I

Vc (t) = -ImXc [cos( f f l t) - cos (a)] + Va cot € [2n - a, 1n)

4.0
— a = 30°

X
3.0

i»
a)
o
tf
2.0

1.0
/ -a-15°

\f \
< a = 0*

c
<0
v -1.0

o
U
-2.0
X
-3.0

-4.0
0 0.005 0.01 0.01S 0.005 0.01 0.015
Tiempo Tiempo

Figura 5.2.- Voltaje en el Capacitor y Corriente en el Inductor

en la Región Inductiva

Las Figuras 5.2 y 5.3 muestran las formas de onda típicas para el voltaje del capacitor

y corrientes en el inductor en la región inductiva y capacitiva. Puede observarse que el grado

de distorsión de la onda de voltaje es mayor cuando el T C S C es operado en la región

inductiva. Otra observación importante es el incremento en la amplitud del voltaje del

capacitor conforme el ángulo de disparo se acerca al punto de resonancia (ver punto de

resonancia en la Figura 5.4), Para el modo de operación de bloqueo (a = 0 o )el voltaje es

senoidal, así como en el modo de operación continuo (a = 90°). La polaridad del voltaje en

el capacitor depende de la zona de operación.

En el comportamiento de la corriente puede apreciarse algo similar a lo observado en

el voltaje del capacitor. El valor del ángulo de disparo determina la dirección de la corriente
a través del reactor controlado y del capacitor, las cuales son opuestas y la magnitud de las
corrientes se incrementa conforme la operación del dispositivo se acerca al punto de
resonancia.

Tiempo Tiempo

Figura 5.3.- Voltaje en el Capacitor y Corriente en el Inductor

en la Región Capacitiva

5.3 IMPEDANCIA F U N D A M E N T A L DEL T C S C

Debido a que la corriente que circula por los tiristores contiene distorsión armónica,

se debe realizar el análisis de Fourier para la obtención de una expresión a frecuencia

fundamental. Así,

i°L{t) = I?sen(at) (5.35)

En donde:

(t)es la corriente de los tiristores a frecuencia fundamental.

I°es el coeficiente de Fourier a frecuencia fundamental dado por la expresión:

1° = 7 í 'l (t)sen(cot)dá>t (5.36)

En donde:

Úi

En las figuras 5.2 y 5.3 puede apreciarse que la onda de corriente de los tiristores

iL (t) cuenta con simetría de media onda impar, por lo que la componente fundamental se

obtiene evaluando (5.36) a partir de (5.22)

sen (<*)
Asen (cot) - A \ eos| ú > 0 t - ^ k sen (cot) dcot
eos

(5.37)

4 k tan sen2 (a) + (k2 + l) sen (2a)

n-2a~\ A
= A
K k2 -1 7T

La impedancia fundamental del TCSC a frecuencia fundamental está definida mediante:

o
VTCSC
70
TCSC
_ no
TCSC J
yyo
TCSC
(5.38)

En donde

es e
-ese ' voltaje fasorial del TCSC a frecuencia fundamental

I es la corriente fasorial de la línea

 Debido a que el voltaje a través del TCSC es Igual al voltaje del capacitor, la ecuación

(5.38) puede escribirse de la siguiente manera:

-j*c % (5.39)

I - íí
-TCSC

Z°rcsc = -J*c + JXC

I

z°csc =-jXc + jXc-f (5.40)

1
En donde X . =
coC

I°c es la corriente fasorial del capacitor a frecuencia fundamental

1°L es la corriente fasorial en el reactor controlado por tirostores ( TCR ) a frecuencia

fundamental.

Sustituyendo (5.37) en (5.40) se obtiene la impedancia para el TCSC a frecuencia

fundamental en estado estable, la cual está compuesta únicamente de la parte reactiva que

depende del ángulo de disparo a , y de las reactancias capacitiva Xc e inductiva XLdel

dispositivo:

4k tanl k\ sen2 (a) + (k2 +l) sen (2a)

A (K -2a
¿rcsc - -JXC + JXC n
J Im{k2-1) 7T

Ak tan k\a-£ sen2 (a) + (k2 +1) sen (2a)

=
-Jxc + J{xc + xiC ) | — j + j
71

(5.41)

En donde:

Xc — = X c + XLC
 2
y
A
*
X IC
_ ^

Xc
Uk2-i)= XL

LC
xc-x L

La Figura 5.4 muestra la curva característica para la impedancia del TCSC a

frecuencia fundamental como función del ángulo de disparo. Los parámetros para el TCSC
utilizados en la obtención de la curva son los siguientes:

- 10 mH, C = 212 pF

150
—
« ico Punto de resonancia
Región Inductiva j
50
£
(J
0
1
i -50
^ Región Capacitiva
-100

-150
0 10 20 30 40 50 60 70 80 90
Anguto de deparo ( grados)

Figura 5.4.-Característica de Impedancia del TCSC a Frecuencia Fundamental

5.4 PRINCIPIOS BÁSICOS DE OPERACIÓN

La Figura 5.5 muestra un esquema básico para un TCSC compuesto por un solo
módulo. En la práctica pueden ser dos o mas módulos, dependiendo de las necesidades que
se tengan. Para reducir el costo, el TCSC puede ser utilizado junto con un capacitor fijo.
Cada módulo puede ser operado en corto circuito,en bloqueo y vernier.

Figura 5.5.-Esquema Básico para el TCSC

5.4.1 Operación en corto circuito

Para lograr este modo de operación, las válvulas del tiristor son disparadas cada 180°

(en ambas direcciones), con esto se logra que la corriente en el reactor sea continua y

senoidal. De esta manera se obtiene un flujo de corriente más grande en el reactor que en el

capacitor, por lo que este modo de operación es utilizado principalmente para proteger al

capacitor de posibles sobrevoltajes.

A este modo de operación también se le conoce como modo "TSR" (Thyristor

Switched Reactor). El diagrama esquemático para este modo de operación se presenta en la

Figura 5.6a

5.4.2 Modo de operación en bloqueo

Para lograr esta condición de operación, se deben mantener bloqueadas las válvulas

de los tiristores, lo cual provoca que la reactancia equivalente del TCSC sea igual a la del

capacitor fijo. Este modo de operación es equivalente a tener un capacitor fijo. El diagrama

esquemático para este modo de operación se presenta en la Figura 5.6b.

.'Ü-YY' A
K
> II
(a) En cortocircuito
— ' m ' <
N
(c) Vernier

•
(b) En Bloqueo

Figura 5.6.-Modos de Operación del T C S C

5.4.3 Modo de operación Vernier

Para lograr este modo de operación, los tiristores deben ser operados en la región

(cmm < a < 180°), de tal forma que la conducción se presenta en una parte del ciclo. El valor

efectivo de la reactancia del TCSC en la región capacitiva aumenta conforme a se va

reduciendo y alcanza un valor máximo para a = am]n. El valor de la reactancia del T C S C

puede llegar a ser hasta tres veces Xc. a,™se encuentra alrededor del valor de a

correspondiente a la resonancia paralelo entre el "TCR" y el capacitor fijo (a frecuencia

fundamental) 2 . El diagrama esquemático para este modo de operación se presenta en la

Figura 5.6c

5.4.4 Problemas de diseño

La selección del tamaño del reactor controlado y del capacitor fijo tiene efectos

importantes en el comportamiento del TCSC. Por razones económicas, la reactancia del

capacitor puede ser lo suficientemente pequeña, tal que la caída de voltaje a través de este

elemento en condiciones de operación normales se encuentre por debajo de un valor

determinado, el cual está definido por los requerimientos de diseño del sistema. Una vez que

la reactancia del capacitor fijo se ha seleccionado, la frecuencia natural del circuito LC

(co0)es directamente dependiente de la inductancia del reactor controlado, por lo que

valores mas bajos de L resultan en valores de frecuencia natural y picos de corriente más

elevados. Mientras que valores más grandes para el reactor provocan una característica de

reactancia más plana.

5.5 VENTAJAS Y APLICACIONES

La capacidad de cambiar el nivel de compensación serie rápidamente mediante la

variación del ángulo de disparo del tiristor, hacen al TCSC muy atractivo en muchas

aplicaciones. Algunas de estas aplicaciones se presentan a continuación.

2
Generalmente el modo de operación Vernier es utilizado únicamente en la región capacitiva
5.5.1 Control de flujo de potencia activa y reactiva.
j

En una línea de transmisión compensada en serie, el nivel de compensación puede

ser utilizado para el control del flujo de potencia activa y reactiva a través de la línea. Esto

se debe a que el TCSC es un compensador serie controlable, por lo,que puede utilizarse

para este propósito.

5.5.2 Mejoramiento de la estabilidad transitoria

El criterio de áreas iguales puede utilizarse para mostrar la capacidad de este

elemento para el mejoramiento de la estabilidad transitoria. Un T C S C puede ajustarse

dinámicamente al nivel de compensación serie total de la línea, incrementándose la energía

de aceleración o desaceleración de una máquina de una manera controlada obteniéndose así

un mejoramiento en la estabilidad transitoria del sistema.

5.5.3 Control de flujo no deseado ( Loop Flow )

Algunas veces en los sistemas de potencia se desea incrementar la cantidad de

potencia a través de rutas definidas. Esto puede resultar en un incremento de potencia no

programado a través de otras rutas paralelas. Este flujo de potencia no deseado es conocido

como "Loop Flow". El control de reactancia continuo y rápido, permite controlar la cantidad

de potencia activa y reactiva que fluye a través de una línea de transmisión por lo tanto el

flujo de potencia no planeado (loop flow) a través de rutas imprevistas puede ser

minimizado localizando adecuadamente el TCSC.

5.5.4 Amortiguamiento de oscilaciones subsíncronas

La frecuencia de oscilaciones torsionales están en el rango de las frecuencias

subsíncronas ( 5-55 Hz ) . Un TCSC puede variar la impedancia de la línea en el rango de

frecuencias subsíncronas. Por lo tanto, basado en una estrategia de control predeterminada,

es posible modular el flujo de potencia activa y reactiva de la línea para mitigar la oscilación

torsional de la flecha.
5.6 ESTRATEGIA DE CONTROL CONVENCIONAL

La estrategia de control convencional para la mitigación de resonancia subsíncrona es

utilizar una señal de error que contenga componentes para todos los modos torsionales
inestables y aplicar esta señal a través de un lazo de control. La Figura 5.7 muestra un
esquema del controlador utilizado para este propósito. La señal utilizada para este propósito
puede ser la velocidad del rotor del generador, la frecuencia o la potencia de salida. En los
resultados presentados en la sección 5.9 se utilizó como señal de retroalimentación la señal
de velocidad para el rotor de la máquina.

El objetivo del esquema es dar amortiguamiento a todos los modos de oscilación de la

flecha, con un margen de fase adecuado para todas las [Link] de operación del
sistema. Esto dificulta la utilización de un solo bloque de control para este propósito, ya que
no es fácil lograr un amortiguamiento para todos los modos de la flecha en todos los niveles
de compensación de línea deseados, por lo que esta técnica no siempre puede ser utilizada
para dar amortiguamiento a todos los modos torsionales del sistema. Esto se debe a que la
acción de control para dar amortiguamiento a algún modo en particular, pudiera excitar
negativamente otro modo de oscilación, lo que hace difícil la sintonización.

retardo de t i e m p o
e n el p r o c e s o d e d i s p a r o
señal de _
A a
referencia ~ -N 1
- «w ^ ->

Medición

Figura 5.7.-Diagrama de Bloques para el Esquema de Control Convencional

5.7 ESTRATEGIA DE CONTROL MODAL

La idea principal en esta estrategia de control es la descomposición de la señal

medida en sus componentes modales y utilizar cada una de estas señales como señal de
retroalimentación en un esquema de control similar al visto en la sección anterior pero de
manera independiente. El proceso de descomposición de esta señal, puede realizarse
mediante manipulación matemática o mediante la utilización de filtros activos.
5.7.1 Cálculo de componentes modales.

Las desviaciones de velocidad de la flecha del grupo turbina-generador Acó pueden

ser utilizadas para obtener las velocidades modales correspondientes. Estas componentes

pueden encontrarse a partir de la ecuación (5.42)

[Affl] = [ 0 ] [ A f f l J (5.42)

En donde Q esta formada por los eigenvectores de [tt Invirtiendo la ecuación

(5,42) se tiene:

[Afií„] = [e] _ 1 [Affl] (5.43)

En donde:

[AO>] = [ACÚhp ACÚ ¡p AALPA ACOlpb ACOG ACO^CF

[A^M] = [A£»ML A¿YM2 ACONÚ AÜ)M4 AO>MS &Ú)M6F

[Q] es u n a matriz de 6 x 6

De la ecuación (5.43) se obtienen las velocidades modales Aco m . El i esimoelemento

del vector [A¿y m ], está relacionado a la velocidad de la masa mediante la ecuación (5.44)

A Ú)¡ = an AÚ)hp + a¡2 Acú1P + anAú)LPA + a¡4AcoLPB+ai5Acoa+a¡6Aú)EX (5.44)

En donde los elementos a^son los elementos de la matriz \Q\ Las desviaciones de

la velocidad de la flecha pueden ser medidas instantáneamente, y como consecuencia, cada

una de las velocidades modales Acom puede ser usada como señal de retroalimentación para

mitigar las oscilaciones de la flecha. En la Figura 5.8 se muestra el diagrama de bloques

para el esquema de control modal.

3
Las matrices H y K se definieron en la sección 2.7
 H, (s)

A¿>„
retardo de tiempo
-> h2(S)
\ en el proceso de disparo
+
\ Aa
1 + sT

¿a>m>
> H,(S)

A cu.
H„ (s)

F i g u r a 5 . 8 . - D i a g r a m a d e B l o q u e s p a r a el E s q u e m a d e C o n t r o l Modal

5.8 DESCOMPOSICIÓN DE LA SEÑAL DE CONTROL UTILIZANDO

FILTROS ACTIVOS

Para este propósito se utiliza un filtro que tiene como salidas las velocidades modales

obtenidas a partir de la señal medida Acó. Un filtro pasa altas es utilizado para eliminar

señales de baja frecuencia como el modo ds oscilación entre áreas. La señal de salida del

filtro pasa altas, se introduce a un grupo de filtros notch para separar la señal de velocidad

modal requerida. Esta aproximación es válida debido a que las frecuencias modales son

conocidas de manera muy precisa, ya que estas no se ven afectadas por la variación en los

parámetros del sistema y las condiciones de operación. Por lo tanto se pueden diseñar filtros

notch para eliminar señales desconocidas.

5.9 MODELO MATEMÁTICO DEL TCSC PARA ESTUDIOS DE

RESONANCIA SUBSÍNCRONA

5.9.1 Aproximación de la característica de reactancia del TCSC

Para la aproximación de la curva característica del T C S C en estudios de resonancia

subsíncrona, se supondrá un generador entregando el 90% de su capacidad en MVA a un

factor de potencia de 0.99 en atraso y un voltaje en terminales de 1.0 p.u.

Se desea compensar el sistema al 68%, esto deberá repartirse entre un banco de

capacitores fijo y el T C S C de la siguiente manera:

a) El 76.5% de la compensación serie será mediante un banco de capacitores fijo. Es

decir X c = 0.29 p.u

b) El 23.5% restante estará dado por el TCSC, lo cual representa una reactancia

equivalente del elemento igual a Xrcsc = 0.0892 p.u. Esto puede lograrse con un

ángulo de di sparo de los tiristores de a 0 = 51.4°

El sistema utilizado para el análisis de valores característicos y su efecto en el fenómeno de

resonancia subsíncrona se muestra en la Figura 5.9

x. R XT
i BUS
v V -Ú INFINITO

Figura 5.9.- Modelo de prueba

La aproximación de la curva característica en la región capacitiva mostrada en la

Figura 5.4, se puede lograr utilizando mínimos cuadrados de segundo orden:

X ( a ) = 2.31+4.55 a -2.31 a 2 p.u. (5.45)

 La curva alrededor del punto de estado estable correspondiente a la ecuación (5.46)
se muestra en la Figura 5.10. Se puede observar que para variaciones no mayores de ±2°,
la aproximación se asemeja bastante a la curva real.

-0.075 0.0891
0.08
-0.0691
-0.0891
-0.085
-0.0891
-0.09
-0.0891
-0.09S -0.0892
-0.1 -0.0892
curva real
-0.105 0.0892 ajuste

-0.11
-0.0892
-0.0892
-0.115
SO 52 SI.38 Sl.< 51.42 51.44
Angulo de disparo Angulo de lloaro

Figura 5.10.- Aproximación de Segundo Orden

Ya que el objetivo es amortiguar las oscilaciones debidas a pequeñas perturbaciones,

la curva puede representarse mediante una aproximación de primer orden alrededor del

punto de estado estable. Está aproximación esta dada por la ecuación:

x(a) = x(a0) + x' (a0)Aa

X(a) = -0.0892 + 0.4080 A a (5.47)

-0.07

-0.07S

-0.08

£V -0.08 S
V
ta
>
9 0.09
%
10 0.095
V
c
8 -0.1
ta
<2
-0.105

~0.11

-0.115
48 49 50 51 52 S3 54
Ángulo de disparo

Figura 5.11.- Aproximación de Primer Orden

 En la Figura 5.11, puede observarse como j la aproximación de primer orden, resulta

ser bastante precisa para valores no muy alejados del punto de estado estable

correspondiente a X r c s c = 0.0892 p.u y un ángulo de disparo de a 0 = 51.4°.

5.9.2 Obtención de las ecuaciones linealizadas para el TCSC

En el sistema de la Figura 5.9, la caída de voltaje a través del T C S C está dada por la

expresión;

= —^/(t) (5.48)

Transformando la ecuación (5.48) al dominio de Laplace y poniendo las ecuaciones en

el marco de referencia qdo se tiene:

— «OW - Ed) = XC («)/„ + «>(VTCSCQ - EQ) (5.49)

a
o

~5(vrcscq -Eq) = Xc [a)iq - ® (Vrcscd -Ed) (5.50)

Finalmente se linealizan las expresiones (5.49) y (5.50) para obtener:

~SAVTCSCÜ = Xc {<x0)tid + Xc {a0)Id0^a + ® 0 Al/ rcsa; + to0EtíoAS + VTCSCgoAa (5.51)

SAV
TCSCQ = XC (« e ) AIQ + x'c (A C ) IqoAa - A>0AVTCSCÜ + (o0EqoAS - VTCSCd0A(Ú (5.52)

5.10 ANALISIS DE EIGENVALORES

La Tabla 5.3 muestra los resultados obtenidos con el control convencional. Puede

observarse como al ir incrementando la compensación en el sistema, se requieren nuevos

ajustes que permitan estabilizar el sistema. Es importante observar que no es posible

estabilizar el sistema en un rango de compensación muy grande con el mismo control,

debido a que se presentan condiciones de inestabilidad en otros modos, por lo que la

sintonización del control es complicada.

 Utilizando una estrategia de control modal,; se obtuvieron los resultados mostrados en
la Tabla 5.4, Cada uno de los compensadores que forman el lazo de control modal tiene una
función de transferencia de primer orden dada por la expresión:

/ / ( 5 ) = * ! ± 5 i ' ( 5 . 5 3 )
V
' 1 + T2S

Los parámetros seleccionados para cada uno de los bloques de adelanto-atrasó se

muestran en la Tabla 5.5. Estos fueron seleccionados mediante la solución del sistema de

ecuaciones ( para cada uno de los bloques ) presentado en el apéndice C

Tabla 5.5.- Parámetros para el Control Modal

FRECUENCIA
MODO * . K '• * • Ti(seg) ' . - T?(seg)
RAD/SEG

M2 . 203 20 0.0042 0.0064

M3 160 -20 0.0048 0,0066

M4 127 "• 20 0.006 0.01

M6 , 99 -30 0.008 0.0172

En el caso de los resultados presentados en la Tabla 5.4, puede observarse como al ir

incrementando la compensación en el sistema, los parámetros para cada uno de los bloques
de adelanto-atraso que componen el lazo de control permanecen constantes. Es importante
recalcar como a diferencia de la estrategia de control convencional, cuando se utiliza la
estrategia de control modal, se puede eliminar efectivamente el fenómeno de resonancia
subsíncrona en todo el rango de compensación.
 Tabla 5.3.- Valores Carácteristicos del Sistema Utilizando la Estrategia de Control Convencional
)

Efecto del Control Convencional en el Amortiguamiento de Oscilaciones

Torsíonales
20% 30% 40% 50%
Porcentaje
Xc= 20.35% Xc= . 46.9% Xc= 60.17% Xc= 68.14%
1
compensación
XTCSC« 79,64% XTCSC- 53.1% XTCSC= 39.82% , XTCSC=31.85% .
Tl=0.08 T1=Q.G8
Sin .
T2=0,0034 T2=0.0034 T2=0.9Q21
Modo Contra*
K=0.52 K=0.52 K=G.55

MR+ -4.82 ± 511.51i -4.86 ± 541.74i -4.89 ± 567.26i -4.91 ± 589.72i

HR- -4.19 ± 242.27i -3.98 ± 212.26i -3.82 ± 186.72i -3.62 ± 164.51Í

M1 -0.23 ± 298.18i -0.23 ± 298.18i -0.23 ± 298.18i -0.23 ± 298-lSi

M2 -0.03 ± 203.08i -0.04 ± 202.95i -0.03 ± 202.91i -0.06 ± [Link]

M3 -0.05 ± 160.64i -0.04 ± 160.61Í -0.03 + 160.64i -0.03 ± 160.39i

M4 -0.04 + 127.03i -0.04 ± 127.02i -0.04 + 127.03i -0.03 ± 127.02¡

M5 -0.07 ± 99.18i -0.03 ± 9 9 . 1 Ü -0.02 ± 99.161 0.05 ± 99.09Í

ME -1.00 ± 8.721 -0.88 ± 9.02i -0.95 ± 9.37Í -0.87 ± 9.74Í

MR 0.00 ± [Link] 0.00 ± [Link] 0.00 ± [Link] 0.00 ± [Link]

60% 70% 80%: 90%

73.45% Xc- 7?.24% ; Xc= 80.08% Xc= . 82.30%
XTC5C=26.54% XTCSC= 22.75% XTCSC^ 13-91% XTCsc= 17.69%
T I =#.08 Tl-0,08
T2=0.0021 T2=0.002l :
INESTABLE INESTABLE
K=0.55 K=0.55

-4.93 ± 610.05i -4.94 ± 628.73i -4.95 ± 646.15i -4.96 ± 662.48Í

-3.38 ± 143.95i -3.23 ± 125.27i -4.54 ± 105.36i -4.79 ± 94.05Í

-0.23 ± 298.18Í -0.23 ± [Link] -0.23 ± 298.18i -0.23 ± 298.18i

-0.05 ± 202.98i -0.05 ± 202.97Í -0.04 ± 202.90Í -0.04 ± 202.91Í

-0.09 ± 160.60i -0.06 ± 160.59Í -0.05 + 1 6 0 . 5 l i -0,05 ± 160.53i

-0.02 ± 127.02i -0.01 ± 126.90Í -0.03 ± 126.93Í -0.04 ± 126.97Í

0.09 + 99.161 0.19 + 99.30i 1.73 ± 101.74Í 2.43 ± 96.53Í

-0.95 ± 10.18Í -1.05 ± 10.69Í -1.70 ± 11.35Í -1.92 ± 12.08Í

0.00 ± [Link] 0.00 ± [Link] -0.00 ± [Link] 0.00 ± [Link]

 Tabla 5.4.- Valores Característicos del Sistema Utilizando la Estrategia de Control Modal

Amortiguamiento de Oscilaciones Torsionales Mediante el Control Modal 1

20% 30% 40% 50% ¡I
Porcentaje
Xc= 20,35% Xc= 46.9% Xc—. 60.17% Xc= 68.14% I
compensación
X TC sc= 79.64% Xrcsc- 53.1% XTCSC- 39.82% XTCSC-31,85% 1

MR+ -4.82 ± 511.51Í -4.86 ± 541.68Í -4.89 ± 567.21Í -4.91 ± 589.701

MR- -4.19 ± 242.27i -3.86 ± 213.27i -3.69 ± 187.32i -3.01 ± 166.33Í

M1 -0.23 ± 298.18i -0.23 ± 298.18i -0.23 ± 298.18i -0.23 ± 298.18i

M2 -0.03 ± [Link] -0.16 ± 202.51Í -0.12 ± 203.18Í -0.06 ± 203.031

M3 -0.05 + 160.64i -0.01 ± 160.30i -0.02 ± 1 6 0 . l l i -0.53 ± 158.891

M4 -0.04 ± 127.03i -0.02 ± 126.79i -0.02 ± 126.72Í -0.02 ± 126.59i

M5 -0.07 ± 99.18¡ -0.03 ± 98.79Í -0.02 ± 98.79Í -0.02 ± 98.771

ME -1.00 ± 8.72i -1.09 ± 9.03i -1.18 ± 9.38i -1.28 + 9.78i

MR 0.00 ± [Link] 0.00 + [Link] 0.00 ± [Link] 0.00 ± [Link]

60% 70% 80% 90%.

Xc~ 73.45% Xc= 77.24% Xc= 80.08% Xc= 82.30%

X T C S C=26.54% X T tsc= 22.75% XTCSC" 1 9 . 9 1 % XJCSC" 17.69%

-4.92 ± 610.031 -4.94 ± 628.72i -4.95 ± 6 4 6 . l l i -4.96 ± 662.45i

-2.96 ± 144.17Í -1.62 ± 122.56i -2.27 ± 107.52i -1.76 ± 90.051

-0.23 ± 298.18i -0.23 ± 298.181 -0.23 ± 298.18i -0.23 ± 298.18i

-0.05 ± 202.99i -0.04 ± 202.97Í -0.04 ± 202.96i -0.04 ± 202.96Í

-0.26 ± 161.15» -0.11 ± 160.84¡ -0.08 ± 160.74i -0.06 ± 160.69i

-0.08 ± 126.20i -1.25 ± 129.48Í -0.16 ± 127.59i -0.07 ± 127.3H

-0.01 ± 98.73Í -0.01 ± 98.64Í -0.09 ± 98.21Í -0.22 ± 99.57Í

-1.40 ± 10.22Í -1.55 ± 10.74Í -1.73 ± 11.34Í -1.96 ± 12.07Í

-0.00 ± [Link] 0.00 ± [Link] 0.00 ± [Link] -0.00 ± [Link]

• -

NOTA: Los modos marcados en color azul indican el modo sobre el que está actuando el

control
 CAPITULO VI

OTROS ESQUEMAS DE AMORTIGUAMIENTO PARA

OSCILACIONES SUBSÍNCRONAS

6.1 INTRODUCCIÓN

A lo largo del trabajo se ha visto como el problema de resonancia subsíncrona está

relacionado con la compensación serie, sin embargo existen algunos puntos importantes que
deben tenerse en consideración.

Primeramente se debe decir que los puntos resonantes son repentinos y los modos
están significantemente separados. Además, las líneas de transmisión compensadas en serie
tienen frecuencias resonantes específicas, las cuales pueden ser localizadas entre los modos.

Otro aspecto importante es que cuando se tienen muchas plantas conectadas a un

sistema, estas tienen frecuencias resonantes diferentes, por lo que tienden a amortiguarse
unas con otras. El caso extremo esta representado por una planta en un sistema radial
compensado en serie.

Existen diferentes maneras de proveer amortiguamiento para reducir pares

transitorios. Muchos de los dispositivos FACTS que existen en el mercado, tienen funciones
adicionales para dar amortiguamiento a estos pares subsíncronos, por lo que este capítulo
presenta algunos tipos de controladores utilizados para el amortiguamiento de oscilaciones
subsíncronas.
6.2 ESQUEMA DE AMORTIGUAMIENTO NGH-SSR

El nombre del esquema NGH, se debe a las iniciales de su inventor ( Narain G.

Hingorani [9]. Los objetivos de este esquema de amortiguamiento son los siguientes:

• Minimizar el par eléctrico subsíncrono, los pares mecánicos y por consiguiente la

torcedura de la flecha.

• Evitar la aparición de oscilaciones, debidas a resonancia subsíncrona de estado

estable.

• Suprimir el de Offset de! capacitor serie, resultado de una falla, liberación de falla,

recierre y otros disturbios, ya que estos alimentan los pares eléctricos subsíncronos.

• Proteger el capacitor serie contra sobrevoltajes

• Reducir los esfuerzos del capacitor, incluyendo sobrevoltajes y el nivel de corriente de

descarga. Eliminar las oscilaciones asociadas con la descarga del capacitor durante el
bypass, por lo tanto reducir el costo del capacitor.

El principio básico de funcionamiento está basado en las siguientes observaciones:

mayor que menor que mayor que

8.33 ms 8.33 ms 8.33 ms
t» T
8.33 ms l t^ 8 . 3 3 ms f

P u l s o s de d i s p a r o

(a)
 (b)

Figura 6.1.- (a) Señal senoidal combinada con una señal de cd (b) Señal
senoidal combinada con una señal de frecuencia subsíncrona

En la Figura 6.1 (a), el voltaje de 60 Hz se combina con un voltaje de CD. Por otra
parte, en la Figura 6.1 (b) la misma señal de voltaje se vuelve a combinar, pero en esta
ocasión con un voltaje de frecuencia subsíncrona. En ambos casos, puede observarse como
algunos semiciclos son más grandes que el semiciclo nominal correspondiente a la
frecuencia de 60 Hz ( 8.33 ms ). Cualquier combinación, con voltajes de CD, voltaje
subsíncrono o voltaje asociado con cualquier oscilación relacionada con estabilidad de baja
frecuencia, tendrá como resultado algunos semiciclos más grandes que los semiciclos de
frecuencia nominal.

El esquema de amortiguamiento para oscilaciones subsíncronas (NGH) consiste

básicamente de un conjunto de tiristores en anti-paralelo conectados en serie con una
pequeña impedancia, que puede ser una resistencia, una inductancia o una combinación de
los dos. Este arreglo se conecta en paralelo con un banco de capacitores serie. En la Figura
6.2 se muestra un diagrama unifilar simplificado para el esquema NGH.
 C a p a c i t o r serie

-
Línea ! • 'yyr^ r ^ | L'nea

K Pulsos
-<— Control ^ Ajuste
de disparo I de tiempo
de disparo

Voltaje
del c a p a c i t o r

Figura 6.2.- Diagrama Básico para el esquema de amortiguamiento NGH

La hipótesis sustentada por el esquema es la siguiente:

El desbalance de carga existente en el capacitor, se transfiere a la inductancia del

sistema, lo que produce oscilaciones. Si la carga, producto del desbalance se elimina del
capacitor serie, el sistema será desintonizado efectivamente para cualquier otra frecuencia
diferente a la nominal, por lo que el control del disparo de los tiristores es diseñado tal que
cuando se detecta el cruce por cero del voltaje del capacitor se mide la duración del
semiciclo siguiente. Tan rápido como el semiciclo excede el tiempo de ajuste ( 8.33 ms para
sistemas de 60 Hz ), se dispara el tiristor correspondiente, para descargar el capacitor y
llevar su corriente a cero. En cada cruce por cero del voltaje del capacitor, se inicia una
nueva medición del tiempo de duración del semiciclo, con la finalidad de descargar el
capacitor en el momento en el que el semiciclo exceda el tiempo de ajuste.

Puede observarse que la operación del dispositivo NGH es independiente para cada
fase, además no requiere de ninguna señal proveniente del generador o detección especifica
de ninguna frecuencia subsíncrona en particular, así como tampoco la presencia de voltaje
de polarización inversa de CD en el capacitor. El control puede localizarse en la plataforma
del capacitor, o a nivel de tierra.

Información más detallada sobre este esquema de amortiguamiento puede

encontrarse en la referencia [18].

6.3 EL T C B R ( THYRISTOR-CONTROLED BRAKING RESISTOR )

El TCBR es una resistencia conectada en paralelo a través de un par de tiristores

conectados en antiparalelo. Cada pierna de un TCBR trifásico, es operado medio ciclo
encendido y medio ciclo apagado. Se emplea para mejorar la estabilidad transitoria del
sistema de potencia y las oscilaciones subsíncronas, mediante la reducción de la energía
disponible total para la aceleración, y por tanto, las desviaciones de velocidad de una
máquina durante una falla.

El TCBR es utilizado en:

• La prevención de inestabilidad transitoria, durante el primer ciclo de la oscilación del

sistema de potencia, mediante la toma de energía que de otra manera sería utilizada
para la aceleración o des-aceleración del generador.

• Mejorar el amortiguamiento para prevenir inestabilidad dinámica debida a

oscilaciones de baja frecuencia entre sistemas interconectados.

• Amortiguamiento de resonancia subsíncrona ( SSR ), debida a la compensación serie

de líneas de transmisión.

• Reducir y amortiguar rápidamente los pares de flecha subsíncronos, por lo que las
posibilidades de un recierre de alta velocidad exitoso aumentan en líneas cercanas a
la planta de generación.

• Facilitar la sincronización de un generador. Un generador fuera de sincronización

puede generar pares de flecha más severos que los producidos por una falla trifásica
en el generador.
 ¡^ Sistema
Generador Línea
f | de potencia
-t^y
v y

Filtro —>- A m >- ^

* TCBR

Aû)r.

Figura 6.3.- Diagrama básico para el esquema TCBR

Las funciones listadas anteriormente maximizan el valor de los recursos de

transmisión y generación estableciendo niveles apropiados y control inteligente. La mejor

localización para un TCBR es cercana al generador que tenga la necesidad de apertura

durante condiciones de inestabilidad transitoria. El esquema para el TCBR se muestra en la

Figura 6.3.

6.3.1 Amortiguamiento de Oscilaciones Subsíncronas

La Figura 6.4 muestra las formas de onda para el ángulo de potencia y desviación de

velocidad correspondientes a una oscilación de potencia de baja frecuencia, incluyendo la

potencia para el TCBR. La desviación de velocidad es el parámetro de control básico, ya que

el T C B R es insertado cuando la desviación de velocidad filtrada es positiva, y es sacado de

operación cuando esta desviación es negativa.

Durante la perturbación, la potencia del TCBR puede ser pequeña o cero,

dependiendo del voltaje de terminales durante el disturbio. Una vez que la perturbación

desaparece, el voltaje se incrementa y consecuentemente la potencia en el TCBR.

Ya que el TCBR es básicamente una resistencia y únicamente consume potencia, esto

hace que el TCBR funcione como un freno cuando la velocidad de la máquina se incrementa.

Por otra parte, el TCBR no puede contribuir en el incremento de la velocidad cuando esta cae

por debajo de su valor nominal. Esto es un rasgo característico del TCBR. Así mismo es

conveniente colocarlo en terminales de la máquina, con la finalidad de obtener una

respuesta rápida durante aceleración transitoria. En la Figura 6.4a se muestra una señal

típica de desviación de velocidad no amortiguada. La señal mostrada, tiene una componente

significante de aproximadamente 2Hz correspondiente a la frecuencia de oscilación del

sistema de potencia, así como frecuencias mecánicas torsionales cuya componente

dominante es de 20Hz y en menor proporción de 33Hz y 43 Hz.

Si el objetivo es controlar la oscilación del sistema de potencia, la medición de la

desviación de velocidad del generador de la planta debe pasarse a través de un filtro pasa
bajas para separar la componente de baja frecuencia y de esta manera controlar las
oscilaciones de baja frecuencia.

1 1 \
1 4-4-i. • -
•1¡
—

1 i ; I j
iiii
^——
i—
i l •S
—i-,
í 1 1 -H- L
A. r n

( a ) desviación de velocidad

Pr *

cot
( b ) Pr, Acó B a j a f r e c u e n c i a Aco

P. A

Cút
( C ) Pr Aa> C o m b i n a d a

Pr "

( d ) Pr Ao> S u b s í n c r o n a

Figura 6.4.- Pulsos de potencia en el TCBR utilizando señales de entrada (a) desviación
de velocidad A® ; (b) Pr con desviación de velocidad A© de baja frecuencia;
(c) Pr Con desviación de velocidad combinada;(d)P r Con desviación de velo-
cidad subsíncrona

Si por el contrario, lo que se desea es controlar los pares subsíncronos o utilizar

recierres de alta velocidad sin que se presenten pares torsionales elevados en la flecha del
generador, la medición de la desviación de velocidad se debe pasar a través de un filtro pasa
altas diseñado para obtener la señal para la desviación subsíncrona de la velocidad
completa. La secuencia de entrada y salida para el TCBR en el caso de desviaciones de
velocidad subsíncronas, se muestra en la Figura 6.4d. Esta estrategia resulta en instantes de
inserción del TCBR más cortos correspondientes a las oscilaciones subsíncronas.

La Figura 6.4c muestra los instantes de inserción del TCBR para la desviación de
velocidad cuando esta no se pasa por el filtro. Es evidente que esta señal tiene como
objetivo amortiguar las oscilaciones de baja frecuencia hasta que estas son lo
suficientemente pequeñas como para cruzar frecuentemente el eje de referencia.

La desviación de velocidad es esencial para la correcta explotación del TCBR. Esta

puede obtenerse a partir de técnicas de medición disponibles como la medición directa de la
desviación de velocidad del generador, o bien derivarse a partir del voltaje de salida del
generador o el voltaje en el punto de conexión del TCBR.

6.4 EL ESTABILIZADOR DINÁMICO

El diagrama esquemático para este dispositivo se muestra en la Figura 6.5 y está

compuesto por una reactancia modulada conectada en paralelo a través de un
transformador al nodo de un generador. La modulación de la reactancia se hace mediante el
disparo de tiristores, que se disparan en respuesta a una señal basada en la desviación de
velocidad del generador.

La suceptancia en paralelo del TCR está aada por la relación:

6 ( C T ) = ™ H ( 6 1 )

Donde o-es el ángulo de conducción y XL es la reactancia del reactor en paralelo. El

ángulo de conducción a está relacionado al ángulo de disparo a por la relación:

a = 2(?r - <z)

El ángulo de disparo a es modulado alrededor de un punto de trabajo en respuesta a

la oscilación del rotor del generador. Cuando no se presentan oscilaciones del rotor, el

compensador dinámico se comporta como una carga reactiva continua.

 El diseño del estabilizador dinámico es tal que genera la cantidad de corriente
necesaria y de fase apropiada para compensar las corrientes de frecuencia subsíncrona en la
armadura del generador debidas a la resonancia de la red.

G enerador
X\ Transformador
V .

- y. A.A Transformador
Estabilizador

>o

Control del
ángulo de disparo
» A

Figura 6.5.- Diagrama básico para el estabilizador dinámico

Algo interesante que se puede observar acerca de la operación del estabilizador

dinámico es que, aún bajo las peores condiciones resonantes, resulta una pequeña
modulación en la corriente del reactor para corrientes de oscilación grandes en la armadura
del generador debida a la amplificación por la resonancia. La amplificación puede ser
máxima si el estabilizador está localizado en un punto en el sistema donde toda la reactancia
inductiva está en un solo lado del estabilizador, mientras que toda la reactancia capacitiva
está en otro lado.

Las ventajas del estabilizador dinámico son las siguientes:

• Los requerimientos de mantenimiento se esperan que sean mínimos

• La operación no es sensible a variaciones en la frecuencia del sistema o temperatura

ambiente.

Las desventajas son:

• No provee protección contra generador de inducción o problemas de pares
transitorios. El amortiguamiento únicamente es proporcionado cuando la oscilación
del rotor está presente.

• El estabilizador introduce corrientes armónicas que pueden ser minimizadas

mediante:

a) Conexión de los TCR en delta

b) Poner dos devanados secundarios, uno conectado en delta y otro conectado

en estrella para eliminar la 5 a y 7 a armónica.

 CAPÍTULO VII

CONCLUSIONES Y RECOMENDACIONES

7.1 INTRODUCCIÓN

El trabajo presenta conceptualmente el problema de resonancia subsíncrona, e ilustra

las consecuencias catastróficas que este puede desencadenar en un sistema de potencia

debido a corrientes de frecuencia subsíncrona circulantes en la en la red de transmisión,

Se ha podido observar mediante el análisis de eigenvalores como la compensación

serie de líneas de transmisión induce modos de frecuencia subsíncrona, de valores
peligrosamente cerca a los modos de oscilación naturales de la flecha. Adlcionalmente se ha
demostrado como eventualmente el modo de la red (subsíncrono) puede inestabilizarse, y
dar paso con esto al fenómeno conocido como "el efecto de generador de inducción".

Por otra parte, se ha estudiado la interacción que existe entre el sistema mecánico y

el sistema eléctrico, cuando se encuentran conectados al sistema diversos elementos de

control, como es el caso del sistema de estabilización o el de los compensadores serie o

paralelo.

Para el caso particular del TCSC, el trabajo presenta una solución completa al

problema de resonancia subsíncrona mediante el control modal.

En este capítulo, se presentan de manera resumida algunos de los puntos más

importantes derivados de este trabajo.

7.2 CONCLUSIONES Y R E C O M E N D A C I O N E S

• La compensación serie de líneas de transmisión, puede provocar el fenómeno de

resonancia subsíncrona, dependiendo del grado de compensación serie en el sistema.

• Para una condición específica de resistencia en el sistema, se puede presentar el

fenómeno conocido como "efecto de generador de inducción". Sin embargo, el límite
de estabilidad por este fenómeno se encuentra muy elevado.

• La presencia del estabilizador de potencia puede excitar algunos modos naturales de

oscilación de la flecha, y dar paso a una condición de resonancia subsíncrona.

• La compensación paralelo en el punto medio de la línea de transmisión, puede

utilizarse como una medida alternativa para incrementar la capacidad de
transferencia de potencia del sistema, sin el riesgo de un problema de resonancia.

• El efecto de generador de inducción, provocado por la autoexcitación del modo de

oscilación de la red, puede ser resuelto a través del sistema de excitación de la
máquina.

• La aplicación de técnicas de control convencional para la solución del problema de

resonancia subsíncrona, es poco efectiva debido a la complejidad en el ajuste de los
parámetros del control.

• La aplicación de técnicas de control modal en la solución global del problema de

resonancia subsíncrona es efectivo y simple.

• Es importante y recomendable realizar estudios de interacciones torsionales en

sistemas en los que se pretenda incrementar la capacidad de transferencia aplicando
técnicas de compensación serie en líneas de transmisión.

• En estudios de resonancia subsíncrona la consideración del amortiguamiento

mecánico es de mucha importancia para la obtención de resultados confiables.
7.3 APORTACIONES

• El trabajo presenta un análisis conceptual sobre el problema de resonancia

subsíncrona, el cual se presenta de manera clara mediante ejemplos numéricos y el

análisis de eigenvalores.

• Se ilustra en detalle el modelado del sistema mecánico, de la red así como de cada
uno de los diferentes componentes requeridos para el estudio de resonancia
subsíncrona.

• Se presenta una solución al problema de auto excitación asociado con la red,

utilizando el sistema de excitación de la máquina.

• Se presenta de una manera clara y detallada un modelo lineal para el TCSC utilizado
para el cálculo de eigenvalores en estudios de resonancia. La relevancia, de este
modelado radica en su simplicidad para el entendimiento conceptual del problema de
resonancia subsíncrona.

• Se realiza un análisis conceptual detallado del TCSC en un sistema máquina bus

infinito.

• Se presenta una técnica de controi modal para la desintonización del fenómeno de

resonancia subsíncrona logrando una solución global del problema de resonancia

subsíncrona, independientemente del nivel de compensación serie del sistema.

7.4 TRABAJOS FUTUROS

• Es importante la realización de estudios de resonancia subsíncrona aplicando técnicas

de simulación en el tiempo, con la finalidad de poder observar formas de onda y otros

detalles que pueden perderse en el análisis en el dominio de la frecuencia.

• Se debe mejorar, y optimizar el algoritmo desarrollado en este trabajo para el

amortiguamiento de oscilaciones torsionales. La mejora debe orientarse a su

codificación en algún lenguaje de programación, así como la validación del mismo

mediante pruebas de campo y/o resultados históricos registrados.

• Se debe implementar un algoritmo multimáquina que permita estudiar el fenómeno a
nivel de sistema.

• Debe analizarse la interacción del sistema mecánico de cada generador con otros
dispositivos FACTS al ser colocados en la red de transmisión.
 APENDICE A
PARÁMETROS PARA EL SISTEMA UTILIZADO

Al MODELO DE PRUEBA

El primer modelo de prueba para estudios de resonancia subsíncrona del IEEE [2], se

ha venido ocupando a lo largo del trabajo. La Figura A . l muestra un diagrama de este

modelo.

% BUS
^ l >-
vV R > INFINITO
F
A

Figura A l . - Primer m o d e l o de p r u e b a para estudios de resonancia subsíncrona

A2 CONSTANTES DE TIEMPO

Las constantes de corto circuito y circuito abierto se especifican en la Tabla A . l

TABLA A l . - Constantes de la máquina

V V . U : ;;;: C O N S T A N T E S D Ê T I E
CIRCUITO A^IERTO^ .. V A L O R : I- -¿I C M R T ^ C M I M R P ' -

4.300 0.4000

0.032 0.0259

0.850 0.1073

,, - > R9, 0.050 0.0463

 A3 IMPEDANCIAS

Las impedancias normalizadas para la red x y para la máquina, se presentan en las

Tablas A2 y A3 respectivamente:

T A B L A A 2 . - I m p e d a n c i a s d e la r e d TABLA A3.- Impedancias de la máquina

-pA^AmujE^: ;:SEP¡jEÍGfAr
iifrfnjii
i i l l i l
: - * 0.02 0.50

, ûî- 0.14 0.14

0.50 1.56
- «•.r:::.-
0.06 0.06

0.35 0.35
....i,,;. » .-.-r

A4 INERCIAS Y CONSTANTES DE RESORTE

La Tabla A4 muestra los valores correspondientes a cada sección de flecha del modelo

utilizado.

TABLA A4.- Parámetros para la flecha del generador

t - .^sgGgo^^lf-Éqlfií,

I F í
131; \ *.?.} n : .í- í r V : ;1 í-£¿ iP lililí ifíiíi31 " ^^ "-•••..•ií'--':-^.
ï v'.'u.— '"'::!:^
V •• • t
« ••i'i.i^..
HP J . . . . 0.092897 19.303
» , ..1 .. • • .> : "V?
0.155589 34.929

, T -I IPA 0.858670 52.038

Î ... '• : i
LP& Î 0.884215 ^fl^Pg-iGE^? i : 70.858

Generador r 0.868495 2.82

-• • * • • Í. r; -
-, » -. %
Excitador » 0.0342165

1 Las impedancias se encuentran normalizadas a la base del generador 892.4 M V A

 APENDICE B

CONDICIONES INICIALES

I x.

Eje d
K /
/ tr,
/

nx.-xc)

I r.

' I

Eje q

F i g u r a B l . - D i a g r a m a f a s o r i a i p a r a e! c á l c u l o d e c o n d i c i o n e s iniciales

B1 Cálculo de condiciones iniciales

Ed = E sen ((5)

Eq = Ecos(S)

V^ = Ed + Xciq

Ka = Eq - Xcid

Vd =V sen (5-0)

VQ =Vcos(S-<¿>)

Id = i sen (S -9)

Iq = ¡cos(5 -e)

Ifa=(Vq + XdId + raIq)l(X3d)

 APENDICE C
ECUACIONES PARA EL AJUSTE DEL TCSC

La función de transferencia para el sistema de estudio esta dada por:

lap/ace
X = Ax + Bu ->> sx = Ax(s) + Bu(s)

(si - A)x(s) = Bu(s)

x(s) = (si - yA)_1Bu(s)

Laplace
Y - Cx Y(s) = Cx(s)

Y(s) = C(sI -A)1Bu(s)

= C(sl - A) 1B - H(s)

u{s)

si se evalúa un punto cualquiera s= a + jco entonces se tiene que:

H(s) = H(CT + ja) = a + jb

U 1 + s7
"i _ _ K y
s
1 + ST2

Figura Cl.* Lazo de control para el esquema de control modal

Ahora bien si:

~ =K y se evalúa H(a + j&) en el diagrama de bloques de la Figura C1 se tiene qL

'2
 1 + (cr + jco) T2K
Ks = a + jb
1 + (cr + jai) T2

Al desarrollar algebraicamente la expresión anterior se tiene que:

a + r2o-a - coT2b - KS - CTT2KKS = 0

coT2a + b + T2crb - COT2KKS = 0

El cual es un sistema de ecuaciones de primer orden con dos incógnitas ( K s y 7 2 ), por lo que

la solución al sistema de ecuaciones dependerá de los valores asiganados a o-,¿yy la

relación entre y J2.

 REFERENCIAS

[1] Hall, M.C. and Hodges, D.A. "Experience with 500 KV sub synchronous resonance
and resulting turbine generator shaft damage. Al Mohave generating station", IEEE PES
Winter meting , Publication CH1066-0-PWR, (76): pp. 22-30, 1976.

[2] " First Benchmark Model for Computer Simulation of Subsynchronous Resonance"
IEEE tracsactions on Power Apparatus and Systems, Vol. PAS-96,5 Sept/Oct 1977

[3] Hisham Mohammad AM Hamdan "Multivariable Frequency Response Approach

Applied To Power System Dynamics and Control" Thesis 1977

[4] Bahrman, M., Larsen, E.V., Piwko, R.J., Patel H.S:, "Experience with HVDC-Turbine

generator torsional interaction at square Buttle", IEEE Transactions, PAS-96(6):1777-

1786. 1980.

[5] L.M. Cannay "A Novel Approach to the Torsional Interaction and Electrical
Damping of the Synchronous Machine" Part I: Theory, IEEE Trans, on PAS, Oct. 1982,

pp. 3630-3638.

[6] I.M. Cannay "A Novel Approach to the torsional Interaction And Electrical Damping

Of The Synchronous Machine" Part II: Application to an arbitrary Network Theory,

IEEE Trans, on PAS, Oct. 1982, pp.3639-3647.

[7] Anderson - Agrawal - Van Ness "Subsynchronous Resonance in Power

Systems'!EEE PRESS 1990.

[8] Lennart Ángquist, Gunnar Ingestrom, Hans-Áke Jónsson " Dinamical Performance

of TCSC Schemes" ABB Power systems CIGRÉ 1996.

[9] Chee Mun Ong " Dynamic Simulation of Electric Machinery" Prentice Hall 1997

[10] K.R. Padiyar, "Análisis of Subsynchronous Resonance in Power Systems "

Academic Publishers, 1999
[11] Antonio Ricardo de Mattos Tenorio, Carlos A. Gama " Resolvendo Problemas de

Ressonáncia Subsíncrona através do TCSC e do Dispositivo NGH-Análise de estratégias

de controle" XV Seminario Nacional de Producao e Transmissao de Energía Elétrica

octubro 1999 Iguacu - Paraná - Brasil.

[12] Ricardo Jorge Davalos Marín "Modelado y Análisis Dinámico del TCSC y su

Aplicación para Mejorar el Comportamiento Electromecánico en Redes Eléctricas". Tesis

de maestría CINVESTAV- Unidad Guadalajara 2000

n
[13] Carlos Alberto Rivera Salamanca Análisis de Interacciones torsionales
Subsíncronas en sistemas de potencia con compensación serie y paralelo". IPN tesis
Doctorado 2000.

[14] E. Fitzgerald, Charles Kingsley, Stephen D. Umans "Máquinas Eléctricas" 5a Ed.

McGraw-Hill

[15] Stephen J. Chapman " Máquinas Eléctricas" McGraw-Hill.

[16] Paul C. Krause - Oleg Wasynczuk - Scott D. Sudhoff "Analysis of Electric

Machinery " IEEE PRESS

[17] Prabha K u n d u r " Power System Stability and Control * McGraw-Hill, Inc.

[18] Narain G. Hingorani - Laszlo GyuGyi " Understanding Facts"

[19] "Damping Subsynchronous Oscillations Using Thyrlstor-Controlled Series

CapacitorJ' thesis University of Toronto.

[20] Alejandro Mendoza, Fernando Sánchez "Análisis y modelado de elementos para

estudios de Resonancia Subsíncrona" UANL RVP 2002

 INDICE DE FIGURAS

FIGURA 2 . 1 SISTEMA DE TRANSMISIÓN COMPENSADO EN SERIE 7

FIGURA 2 . 2 CIRCUITO EQUIVALENTE PARA UNA MÁQUINA SÍNCRONA COMO FUNCIÓN DEL

DESLIZAMIENTO S 12

FIGURA 2 . 3 MÁQUINA SÍNCRONA 16

FIGURA 2 . 4 DIAGRAMA SIMPLIFICADO DE UNA LÍNEA COMPENSADA EN SERIE 23

FIGURA 2 . 5 DIAGRAMA DE BLOQUES PARA EL SISTEMA DE EXCITACIÓN 24

FIGURA 2 . 6 a DIAGRAMA DE BLOQUES PARA EL ESTABILIZADOR 25

FIGURA 2 . 6 b DIAGRAMA EQUIVALENTE 25

FIGURA 2 . 7 a DIAGRAMA DE BLOQUES PARA EL FILTRO TORSIONAL 26

FIGURA 2 . 7 b D I A G G R A M A EQUIVALENTE 26

FIGURA 2 . 8 MODELO DE FLECHA PARA ESTUDIOS DE RESONANCOA SUBSÍNCRONA 27

FIGURA 2 . 9 MODELO MASA-RESORTE PARA LA FLECHA DEL GRUPO TURBINA-GENERADOR 27

FIGURA 3 . 1 a C I R Q U I T O EQUIVALENTE PARA EL EJE d 35

FIGURA 3 . l b C I R C U I T O EQUIVALENTE PARA EL EJE q 35

FIGURA 3 . 2 CIRCUITOS EQUIVALENTES PARA EL ESTUDIO DE GENERADOR DE I N D U C C I Ó N 38

FIGURA 3 . 3 LÍMITE DE ESTABILIDAD POR EFECTO DE GENERADOR DE INDUCCIÓN 40

FIGURA 3 . 4 EFECTO DE LA VARIACIÓN DE LAS RESISTENCIAS DE ROTOR SOBRE EL

AMORTIGUAMIENTO DEL MODO SUBSÍNCRONO 40

FIGURA 3 . 5 REGIÓN DE INESTABILIDAD POR RESONANCIA SUBSÍNCRONA 43

FIGURA 3 . 6 LÍMIRTES DE ESTABILIDAD POR EFECTO DE GENERADOR DE I N D U C C I Ó N Y RESONANCIA

SUBSÍNCRONA 44

FIGURA 3 . 7 DIAGRAMA DE BODE PARA LA FUNCIÓN DE TRANSFERENCIA DEL FILTRO 45

FIGURA 4 . l a SISTEMA SIMPLE DE DOS MÁQUINAS COMPENSADO EN PARALELO 48

FIGURA 4 . l b D I A G R A M A FASORIAL 48

FIGURA 4 . 2 CURVA P- 8 49

FIGURA 4 . 3 LÍNEA DE TRANSMISIÓN NO COMPENSADA 50

FIGURA 4 . 4 SISTEMA DE TRANSMISIÓN CON COMPENSACIÓN SHUNT 52

FIGURA 4 . 5 a DIAGRAMA DE BLOQUES PARA EL CEV 55

FIGURA 4 . 5 b D I A G R A M A DE BLOQUES EQUIVALENTE 55

FIGURA 5 . 1 CIRCUITO ELÉCTRICO PARA UN TCSC 60

FIGURA 5 . 2 VOLTAJE EN EL CAPACITOR Y CORRIENTE EN EL INDUCTOR EN LA REGIÓN I N D U C T I V A . . . 68

FIGURA 5 . 3 VOLTAJE EN EL CAPACITOR Y CORRIENTE EN EL INDUCTOR EN LA REGIÓN CAPACITIVA .. 6 9

FIGURA 5 . 4 CARACTERISTICA DE IMPEDANCIA DEL TCSC A FRECUENCIA FUNDAMENTAL 72

FIGURA 5 . 5 ESQUEMA BÁSICO PARA EL TCSC 72

FIGURA 5 . 6 MODOS DE OPERACIÓN DEL TCSC 73

FIGURA 5 . 7 DIAGRAMA DE BLOQUES PARA EL ESQUEMA DE CONTROL CONVENCIONAL 76

FIGURA 5 . 8 DIAGRAMA DE BLOQUES PARA EL ESQUEMA DE CONTROL MODAL 78

FIGURA 5 . 9 MODELO DE PRUEBA 79

FIGURA 5.10APROXIMACIÓN DE SEGUNDO ORDEN 80

FIGURA 5 . 1 1 APROXIMACIÓN DE PRIMER ORDEN 80

FIGURA 6. l a SEÑAL SENOIDAL COMBINADA CON UNA SEÑAL DE CD 86

FIGURA 6 . 1 bSEÑAL SENOIDAL COMBINADA CON UNA SEÑAL DE FRECUENCIA SUBSÍNCRONA 87

FIGURA 6 . 2 DIAGRAMA BÁSICO PARA EL ESQUEMA DE AMORTIGUAMIENTO NGH 88

FIGURA 6 . 3 DIAGRAMA BÁSICO PARA EL ESQUEMA DEL TCBR 90

FIGURA 6 . 4 PULSOS DE POTENCIA EN EL TCBR UTILIZANDO SEÑALES DE ENTRADA 91

FIGURA 6 . 5 DIAGRAMA BÁSICO PARA EL ESTABILIZADOR DINÁMICO 93

FIGURA A l PRIMER MODELO DE PRUEBA DE IEEE PARA ESTUDIOS DE RESONANCIA SUBSÍNCRONA.. 9 9

FIGURA B l DIAGRAMA FAS O RIA L PARA EL CÁLCULO DE CONDICIONES INICIALES 101

FIGURA C L LAZO DE CONTROL PARA EL ESQUEMA DE CONTROL MODAL 102

 INDICE DE TABLAS

TABLA 3 . 1 PARÁMETROS DE LA MÁQUINA 36

TABLA 3 . 2 VALORES CARÁCTERISTIC05 PARA EL ANÁLISIS DE GENERADOR DE INDUCCIÓN 36

TABLA 3 . 3 EIGENVALORES PARA EL SISTEMA COMPENSADO AL 6 0 % CONSIDERANDO EL

ESTABILIZADOR DE LA MÁQUINA 39

TABLA 3 . 4 VALORES CARACTERISTICOS PARA EL ANÁLISIS DE RESONANCIA SUBSÍNCRONA 42

TABLA 3 . 5 VALORES CARACTERISTICOS PARA EL ANÁLISIS DEL PSS 46

TABLA 4 . 1 VALORES CARACTERISTICOS PARA EL ANÁLISIS DEL CEV EN LAS INTERACCIONES

TORSIONALES 56

TABLA 5 . 1 CORRENTE EN EL INDUCTOR EN ESTADO ESTABLE DURANTE UN CICLO 67

TABLA 5 . 2 VOLTAJE EN EL CAPACITOR DURANTE UN CICLO EN ESTADO ESTABLE 68

TABLA 5 . 3 VALORES CARÁCTERISTICOS PARA EL SISTEMA UTILIZANDO LA ESTRATEGIA DE

CONTROL CONVENCIONAL 83

TABLA 5 . 4 VALORES CARACTERISTICOS PARA EL SISTEMA UTILIZANDO LA E5TRATEGIA DE

CONTROL MODAL 84

TABLA 5 . 5 PARÁMETROS PARA EL CONTROL MODAL 82

TABLA A l CONSTANTES DE MÁQUINA 99

TABLA A 2 IMPEDANC1AS DE LA RED 100

TABLA A3 IMPEDANCIAS DE LA MÁQUINA 100

TABLA A 4 PARÁMETROS PARA LA FLECHA DEL GENERADOR 100