PDF

COMPENSADOR ESTTICO SNCRONO ( STATCOM ):
MODELADO Y SIMULACIN
EDGAR ENRIQUE SIERRA RODRGUEZ
UNIVERSIDAD INDUSTRIAL DE SANTANDER
FACULTAD DE INGENIERAS FSICO - MECNICAS
ESCUELA DE INGENIERAS ELCTRICA, ELECTRNICA Y
TELECOMUNICACIONES
BUCARAMANGA
2006
COMPENSADOR ESTTICO SNCRONO ( STATCOM ):
MODELADO Y SIMULACIN
EDGAR ENRIQUE SIERRA RODRGUEZ
Trabajo de grado para optar al ttulo de Ingeniero Electricista
Director
HERMANN RAL VARGAS TORRES
Doctor Ingeniero Electricista
UNIVERSIDAD INDUSTRIAL DE SANTANDER
FACULTAD DE INGENIERAS FSICO - MECNICAS
ESCUELA DE INGENIERAS ELCTRICA, ELECTRNICA Y
TELECOMUNICACIONES
BUCARAMANGA
2006
A Dios por darme la fortaleza
A la memoria de Denis Rodrguez
A mi pap
A mi hija Denis
A Alba Yaneth con su apoyo incondicional
A mis hermanos Claudia y Juan Carlos
AGRADECIMIENTOS
Quiero extender el mas grande agradecimiento a las personas que siempre estuvieron
presentes cuando mas las necesitaba , a mi hermana Claudia Sierra y su esposo Marcelo
Vergara.
A mis Tias Denis, Euderis y Aleida, siempre atentas.
A mis primos Pedro, Alma, Yelitza, y a todos los que en este proceso aportaron con sus
palabras de nimo y hermandad.
A mis amigos Guido Gonzales, Guido Polo, Dayana, Rafael Peate.
Al Dr. Hermann Ral Vargas, director del proyecto, quien me propuso la idea de vin-
cularme y entrar al mundo del ATP/ATPDraw y por darme la orientacin necesaria.
A los Ingenieros Carlos Rodrguez y Lizeth Meja, cuyas orientacines fuern bastante
precisas, importantes y oportunas.
RESUMEN
TTULO: COMPESADORESTTICOSNCRONO: MODELADOYSIMULACIN
*
AUTOR: EDGAR ENRIQUE SIERRA RODRGUEZ
**
PALABRAS CLAVES: Sistemas de Potencia, FACTS, STATCOM, ATP/EMTP.
DESCRIPCIN:
El rpido desarrollo de la tecnologa de la electrnica de potencia proporciona oportuni-
dades de desarrollar nuevos equipamientos para mejorar el desempeo de los sistemas
de potencia actuales. Durante las ltimas dcadas han sido propuestos e implementados
varios dispositivos de control llamados Flexible AC Transmission Systems (FACTS).
Los dispositivos FACTS se pueden emplear para el control del ujo de potencia, con-
trol de lazos de ujo, regulacin de tensin, mejoramiento de la estabilidad transitoria
y amortiguamiento de las oscilaciones de potencia.
En este trabajo se describe la teora y la tcnica de modelado de un dispositivo para
los Sistemas de Transmisin Flexible en AC (FACTS), llamado, Compensador Esttico
Sncrono (STATCOM), utilizando un programa de simulacin para transitorios electro-
magnticos (ATP/EMTP). El STATCOM es una fuente inversora de tensin de estado
slido acoplado con un transformador, conectado en paralelo con una lnea de trans-
misin. El STATCOM inyecta una corriente casi senoidal, de magnitud variable, en el
punto de conexin. Esta corriente inyectada esta casi en cuadratura con la tensin de la
lnea, mientras emula una reactancia inductiva o capacitiva en el punto de conexon con
la lnea de transmisin. La funcionalidad del STATCOM se verica por la regulacin
del ujo de corriente reactiva a travs de l. Esto es til para regular la tensin de la
lnea.
Tambin se presenta un modelo reducido del STATCOM, implementado bajo una he-
rramienta de simulacin de sistemas de potencia ante grandes perturbaciones (HAGP).
*
Proyecto de Grado
**
Facultad de Ingenieras Fsico-Mecnicas. Escuela de Ingenieras Elctrica, Elctronica y de Tele-
comunicaciones. Hermann Ral Vargas Torres.
ABSTRACT
TITLE: STATIC SYNCHRONOUS COMPENSATOR: MODELLING AND SIMU-
LATION
*
AUTHOR: EDGAR ENRIQUE SIERRA RODRGUEZ
**
KEY WORDS: Power Systems, FACTS, STATCOM, ATP/EMTP.
DESCRIPTION:
The rapid development of power electronics technology provides opportunities to de-
velop new power equipment to improve the performance of the actual power systems.
During the last decade, a number of control devices called Flexible AC Transmission
Systems (FACTS) technology have been proposed and implemented. FACTS devices
can be used for power ow control, loop-ow control, voltage regulation, enhancement
of transient stability and damping of power oscillations.
This work describes the theory and the modelling technique of a Flexible Alternating
Currents Systems (FACTS) device, namely, STATic synchronous COMpensator (STAT-
COM) using un electromagnetic transients program (ATP/EMTP) simulation package.
The STATCOM, a solid-state voltage source inverter coupled with a transformer, is tied
to a transmission line. A STATCOM injects an almost sinusoidal current, of variable
magnitude, at the point of connection. This injected current is almost in quadrature with
the line voltage, thereby emulating an inductive or a capacitive reactance at the point
of connection with the transmission line. The functionality of the STATCOM model is
veried by regulating the reactive current ow through it. This is useful for regulating
the line voltage.
A reduced model of the STATCOM is also presented, implemented under a tool of
simulation of systems of power before big interferences (HAGP).
*
Degree Project
**
School of Electrical Engineering. Hermann Ral Vargas Torres.
ndice general
1. Introduccin 1
1.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2. Historia de los dispositivos FACTS . . . . . . . . . . . . . . . . . . . . 2
1.3. Tipos de dispositivos FACTS . . . . . . . . . . . . . . . . . . . . . . . 3
1.4. Historia del STATCOM . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.5. Motivacin del trabajo . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.6. Objetivos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2. El STATCOM 12
2.1. Funcionamiento del STATCOM . . . . . . . . . . . . . . . . . . . . . 13
2.2. Componentes del STATCOM . . . . . . . . . . . . . . . . . . . . . . . 14
2.2.1. Inversores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2.1.1. Funcin . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2.1.2. Vlvulas utilizadas . . . . . . . . . . . . . . . . . . . 15
2.2.1.3. Topologas . . . . . . . . . . . . . . . . . . . . . . . 16
2.2.1.4. Frecuencia de conmutacin . . . . . . . . . . . . . . 18
2.2.2. Capacitor del lado de corriente continua . . . . . . . . . . . . . 18
2.2.2.1. Funcin . . . . . . . . . . . . . . . . . . . . . . . . 18
2.2.2.2. Topologas . . . . . . . . . . . . . . . . . . . . . . . 18
2.2.3. Transformadores . . . . . . . . . . . . . . . . . . . . . . . . . 19
I
NDICE GENERAL II
2.2.3.1. Transformadores en zig-zag . . . . . . . . . . . . . . 20
2.2.3.2. Transformadores en estrella - delta . . . . . . . . . . 24
3. Anlisis detallado del STATCOM de 6 y 12 pulsos 26
3.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2. Fuente inversora de Tensin (FIT) . . . . . . . . . . . . . . . . . . . . 27
3.3. Anlisis armnico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
3.4. Compensador esttico sncrono (STATCOM) basado en un inversor de
seis pulsos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.4.1. Intercambio de potencia reactiva . . . . . . . . . . . . . . . . . 33
3.4.1.1. Anlisis de las seales de corrientes AC . . . . . . . 34
3.4.1.2. Corriente en el capacitor . . . . . . . . . . . . . . . 38
3.4.1.3. Tensin DC en el capacitor . . . . . . . . . . . . . . 39
3.4.2. Intercambio de potencia activa y reactiva . . . . . . . . . . . . 41
3.5. Conversor de 12 pulsos . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.5.1. Seales de corriente AC . . . . . . . . . . . . . . . . . . . . . 47
4. Control vectorial de inversores trifsicos 49
4.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
4.2. Representacin de seales elctricas en el espacio vectorial . . . . . . . 51
4.2.1. Sistemas trifsicos a cuatro hilos . . . . . . . . . . . . . . . . . 51
4.2.2. Sistema trifsico de tres hilos . . . . . . . . . . . . . . . . . . 54
4.3. Estados de conmutacin en inversores VSI . . . . . . . . . . . . . . . . 56
4.4. Modulacin de vectores espaciales . . . . . . . . . . . . . . . . . . . . 61
5. Modelado del STATCOM 63
5.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.2. Representacin vectorial de cantidades trifsicas instantneas . . . . . . 64
NDICE GENERAL III
5.3. Ecuaciones del circuito equivalente . . . . . . . . . . . . . . . . . . . 67
5.4. Sistemas de control del STATCOM . . . . . . . . . . . . . . . . . . . . 68
5.4.1. Tcnica de control por intermedio de la modulacin SVPWM . 69
5.4.1.1. Implementacin digital del modelo del STATCOM . . 70
5.4.2. Tcnica de control de fase . . . . . . . . . . . . . . . . . . . . 73
5.4.2.1. Componentes bsicos del sistema de control . . . . . 74
5.4.2.2. PLL, phase locked loop . . . . . . . . . . . . . . . . 74
5.4.2.3. Control de reactivos . . . . . . . . . . . . . . . . . . 75
5.4.2.4. Generacin de la seal de desfasaje . . . . . . . . . . 77
6. Modelo reducido del STATCOM 78
6.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
6.2. Modelos del STATCOM . . . . . . . . . . . . . . . . . . . . . . . . . 79
6.2.1. Modelo de estabilidad transitoria . . . . . . . . . . . . . . . . . 79
6.2.2. Modelo en estado estable . . . . . . . . . . . . . . . . . . . . . 81
6.2.2.1. Sistema de 3 mquinas y 9 nodos . . . . . . . . . . . 82
6.3. Anlisis dinmico de estabilidad . . . . . . . . . . . . . . . . . . . . . 85
7. Resultados de simulaciones 90
7.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.2. Ejemplos implementados en ATP/EMTP y HAGP . . . . . . . . . . . . 91
7.2.1. Caso1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.2.2. Caso 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
7.2.2.1. Modelo del STATCOM en el software HAGP . . . . 100
7.2.2.2. Tcnica de control de fase . . . . . . . . . . . . . . . 101
7.2.2.3. Esfuerzos en los elementos de conmutacin . . . . . 104
7.2.3. Caso 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
7.2.3.1. Inversor . . . . . . . . . . . . . . . . . . . . . . . . 105
NDICE GENERAL IV
7.2.3.2. Tranformadores . . . . . . . . . . . . . . . . . . . . 105
7.2.3.3. Resultados de la respuesta dinmica del STATCOM . 109
7.2.3.4. Esfuerzos en los elementos de conmutacin . . . . . 111
8. Conclusiones y sugerencias 115
8.1. Aportes originales del trabajo . . . . . . . . . . . . . . . . . . . . . . . 116
8.2. Sugerencias para posteriores desarrollos . . . . . . . . . . . . . . . . . 116
Bibliografa 118
A. Ecuaciones diferenciales para incluir en el HAGP 122
A.1. Solucin en estado estable . . . . . . . . . . . . . . . . . . . . . . . . 122
A.2. Simulacin dinmica . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
A.3. Ecuaciones dinmicas del modelo del STATCOM . . . . . . . . . . . . 127
B. Gua de usuario del HAGP 129
B.1. Introduccin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
B.2. Instalacin y Ejecucon . . . . . . . . . . . . . . . . . . . . . . . . . . 129
B.3. Datos de entrada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
C. Modelos y libreras en ATP 132
C.1. Caso de simulacin 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . 132
ndice de guras
1.1. Forma bsica del Statcom . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2. FACTS- primera generacin -paralelo . . . . . . . . . . . . . . . . . . 4
1.3. FACTS- primera generacin -serie . . . . . . . . . . . . . . . . . . . . 5
1.4. FACTS- segunda generacin -paralelo-STATCOM . . . . . . . . . . . . 5
1.5. FACTS- Segunda generacin- serie-SSSC . . . . . . . . . . . . . . . . 6
1.6. FACTS-Tercera generacin-UPFC . . . . . . . . . . . . . . . . . . . . 6
1.7. FACTS- Cuarta generacin-IPFC . . . . . . . . . . . . . . . . . . . . . 7
1.8. Caracteristicas V vs I del SVC y el STATCOM . . . . . . . . . . . . . 7
2.1. Diagrama simplicado de dos fuentes de tensin conectados por una
impedancia X
L
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2. Diagramas fasoriales del sistema elctrico y el STATCOM . . . . . . . 14
2.3. Inversor de seis pulsos . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.4. Inversor trifsico compuesto por 3 inversores puente H monofsicos . . 16
2.5. Inversor de tres niveles NPC . . . . . . . . . . . . . . . . . . . . . . . 17
2.6. Topologas de capacitores del lado de corriente continua . . . . . . . . 19
2.7. Topologas con capacitores individuales por cada inversor . . . . . . . . 20
2.8. STATCOM de 12 pulsos utilizando transformadores zig-zag . . . . . . 20
2.9. Esquemas de conexin de los transformadores zig-zag . . . . . . . . . 21
V
NDICE DE FIGURAS VI
2.12. Diagrama de un STATCOM quasi-24 pulsos (esquema de conexiones
estrella - delta) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.1. Fuente inversora de tensin . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2. Seales de disparo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.3. Formas de onda de las tensiones lnea a lnea . . . . . . . . . . . . . . 28
3.4. Formas de onda de las tensiones lnea a neutro . . . . . . . . . . . . . . 29
3.5. Circuito equivalente para la secuencia 1-5-6 . . . . . . . . . . . . . . . 30
3.8. STATCOM con inversor de seis pulsos . . . . . . . . . . . . . . . . . . 34
3.9. Tensiones de fase del compensador y del sistema AC . . . . . . . . . . 35
3.10. Corrientes entre el sistema AC y el compensador . . . . . . . . . . . . 38
3.11. Tensin del capacitor a) generando potencia reactiva; b) absorbiendo
potencia reactiva . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.12. Corrientes del capacitor a) generando potencia reactiva; b) absorbiendo
potencia reactiva . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.13. Tensiones de fase del compensador y del sistema AC . . . . . . . . . . 42
3.14. a) v
ab
(t) y v
abY
(t)
2
; b) Tensin del conversor de 12 pulsos . . . . . . . 46
3.15. STATCOM de 12 pulsos . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.1. Tensiones de fase de un sistema trifsico . . . . . . . . . . . . . . . . . 52
4.2. Evolucin temporal de las tensiones de un sistema trifsico . . . . . . . 53
4.3. Evolucin temporal del vector suma . . . . . . . . . . . . . . . . . . . 53
4.4. Relacin entre las tensiones de fase y tensiones de lnea . . . . . . . . . 55
4.5. Inversor trifsico formado por tres semipuentes . . . . . . . . . . . . . 57
4.6. Estados de los interruptores . . . . . . . . . . . . . . . . . . . . . . . . 58
NDICE DE FIGURAS VII
4.7. Posicin de los diferentes vectores . . . . . . . . . . . . . . . . . . . . 60
4.8. Descomposicin del vector intermedio
V
re f
situado en el sector 1 . . . . 61
5.1. Marco de referencia . . . . . . . . . . . . . . . . . . . . . . . . 65
5.2. Marco de referencia d q . . . . . . . . . . . . . . . . . . . . . . . . . 66
5.3. Esquema general del STATCOM . . . . . . . . . . . . . . . . . . . . . 67
5.4. Esquema general del sistema de control del STATCOM . . . . . . . . . 70
5.5. Controladores adicionales para las seales de referencia del sistema de
control . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
5.6. Modelos del controlador desacoplados . . . . . . . . . . . . . . . . . . 72
5.7. Controlador de corrientes . . . . . . . . . . . . . . . . . . . . . . . . . 72
5.8. Esquema de control del STATCOM para el control de fase . . . . . . . 73
5.9. Esquema del PLL, phase locked loop. . . . . . . . . . . . . . . . . . . 75
5.10. Esquema del PLL+ seal de desfasaje. . . . . . . . . . . . . . . . . . . 76
5.11. Controlador P-I de potencia reactiva . . . . . . . . . . . . . . . . . . . 77
6.1. Modelo reducido del STATCOM con control de tensin PWM . . . . . 80
6.2. Tcnica de control PWM para control de tensin . . . . . . . . . . . . 81
6.3. Sistema 9 barras y 3 mquinas de la WSSC . . . . . . . . . . . . . . . 83
6.4. Diagrama de ujo simplicado del HAGP . . . . . . . . . . . . . . . . 89
7.1. Sistema de potencia para el caso de simulacin 1 . . . . . . . . . . . . 92
7.2. Sistema de control para el caso de simulacin 1 . . . . . . . . . . . . . 93
7.3. Seales de referencia del sistema de control para el caso de simulacin 1 93
7.4. Tensines del sistema AC y a la salida del conversor . . . . . . . . . . . 94
7.5. Armnicos en la tensin de salida del compensador . . . . . . . . . . . 95
7.6. Componentes en cuadratura de las corrientes de la carga y del compen-
sador . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
NDICE DE FIGURAS VIII
7.7. Componente en cuadratura de la corriente suministrada por la fuente de
tensin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
7.8. Compensacin dinmica del factor de potencia de la carga . . . . . . . 96
7.9. Dinmica de la tensin del capacitor en el lado DC del compensador . . 96
7.10. Conguracin de un STATCOM de 12 pulsos . . . . . . . . . . . . . . 97
7.11. Sistema de potencia para el caso de simulacin 2 . . . . . . . . . . . . 99
7.12. Regulador Automtico de Tensin, AVR . . . . . . . . . . . . . . . . . 99
7.13. Sistema de control de fase . . . . . . . . . . . . . . . . . . . . . . . . 99
7.14. Magnitud de tensin en la barra 8 . . . . . . . . . . . . . . . . . . . . . 100
7.15. ngulo del rotor del generador . . . . . . . . . . . . . . . . . . . . . . 101
7.16. Inclusin del STATCOM en barra 8 a los 0,5 seg . . . . . . . . . . . . . 102
7.17. Forma de onda de la tensin a la salida del inversor . . . . . . . . . . . 102
7.18. Espectro armnico de tensin en la salida del inversor . . . . . . . . . 102
7.19. ngulo del Generador . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
7.20. ngulo de desfase del control . . . . . . . . . . . . . . . . . . . . . . 103
7.21. Corriente en un dispositivo de conmutacin . . . . . . . . . . . . . . . 104
7.22. Frecuencia en radianes por segundo . . . . . . . . . . . . . . . . . . . 107
7.23. Seal de sincronismo del PLL . . . . . . . . . . . . . . . . . . . . . . 107
7.24. Potencia imaginaria instantnea inicial . . . . . . . . . . . . . . . . . . 107
7.25. Potencia activa instantnea inicial . . . . . . . . . . . . . . . . . . . . 108
7.26. Tensin del capacitor . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
7.27. Tensin del sistema V
s
, tensin 48 pulsos y corriente de compensacin Ic109
7.28. Potencia imaginaria instantnea . . . . . . . . . . . . . . . . . . . . . 110
7.29. Potencia activa instantnea . . . . . . . . . . . . . . . . . . . . . . . . 110
7.30. Tensin del capacitor . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
7.31. Tensin del sistema V
s
, tensin 48 pulsos y corriente de compensacin Ic111
NDICE DE FIGURAS IX
7.32. Corriente que circula por uno de los elementos de conmutacin, en
modo inductivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
7.33. Corriente que circula por uno de los elementos de conmutacin, en
modo capacitivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
7.34. Corriente que circula por uno de los elementos de conmutacin, con
referencia de potencia reactiva cero . . . . . . . . . . . . . . . . . . . . 112
7.35. Tensiones de los elementos de conmutacin, en modo inductivo . . . . 113
7.36. Tensiones de los elementos de conmutacin, en modo capacitivo . . . . 113
7.37. Tensiones de los elementos de conmutacin, con referencia de potencia
reactiva cero . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
A.1. Circuito equivalente del STATCOM . . . . . . . . . . . . . . . . . . . 123
A.2. Tcnica de control PWM para control de tensin . . . . . . . . . . . . 123
C.3. Caso de simulacin 3 - parte 1 . . . . . . . . . . . . . . . . . . . . . . 155
Lista de tablas
2.1. Valores de x y de y en funcin de . . . . . . . . . . . . . . . . . . . . 22
4.1. Estados de conmutacin . . . . . . . . . . . . . . . . . . . . . . . . . 59
6.1. Datos iniciales de los nodos . . . . . . . . . . . . . . . . . . . . . . . 83
6.3. Datos de las lneas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
6.5. Tensiones nodales sin la inclusin del STATCOM . . . . . . . . . . . . 84
6.7. Flujos de potencia por las lneas sin la inclusin del STATCOM . . . . 84
6.9. Valores nodales con la inclusin del STATCOM en el nodo 5 . . . . . . 86
6.11. Flujos de potencia por las lneas con la inclusin del STATCOM . . . . 86
7.1. Parmetros del sistema de referencia . . . . . . . . . . . . . . . . . . . 92
7.2. Parmetros de la maquina sncrona, con bases de 200MVA y 13.8kV . . 98
7.3. Parmetros del sistema de potencia en p.u., con potencia base de 240
MVA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
7.4. Parmetros del STATCOM en p.u., con potencia base de 150 MVA y
12 kV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
7.5. Secuencia de eventos de la simulacin caso 2 . . . . . . . . . . . . . . 100
7.6. Secuencia de eventos de la simulacin . . . . . . . . . . . . . . . . . . 106
X
Lista de Smbolos
A
D
Particin del jacobiano correspondiente a derivadas
de ecuaciones dinmicas con respecto variables de es-
tado
B
D
de ecuaciones dinmicas con respecto a tensiones
C
D
de ecuaciones de red con respecto a variables de esta-
do
F Vector de ecuaciones dinmicas
G Vector de ecuaciones de red
J Jacobiano
Y Matriz de admitancias de la red del sistema. Particin
del jacobiano correspondiente a derivadas de ecuacio-
nes de red con respecto a tensiones
i Vector de corrientes del sistema
v Vector de tensiones del sistema
x Vector de estado del sistema
V Fasor de tensin de barra
I Fasor de corriente de barra

H Constante de inercia
I Magnitud de fasor de corriente de barra
I
I
Componente imaginaria de fasor de corriente de barra
Y Admitancia
I
R
Componente real de fasor de corriente de barra
K
v
Ganancia del controlador de ndice de modulacin
XI
Lista de Smbolos XII
L Inductancia
P Potencia activa
Q Potencia reactiva
R Resistencia
V Magnitud de fasor de tensin
V
m
Valor pico de onda de tensin
V
I
Componente imaginaria de fasor de tensin de barra
V
DC
Tensin en el lado DC del STATCOM
V
R
Componente real de fasor de tensin de barra
X Reactancia
i Corriente instantnea
x Variable cualquiera
r, s, t (subndice) Denota cantidades de las fases de tensiones
d (subndice) Denota cantidades de eje directo
q (subndice) Denota cantidades del eje en cuadratura
t (subndice) Denota cantidades en terminales del generador
max(subndice) Denota cantidades de valor mximo
mn(subndice) Denota cantidades de valor mnimo
re f (subndice) Denota cantidades de referencia
R(subndice) Denota el valor de la componente real
I(subndice) Denota el valor de la componente imaginaria
k (superndice) Denota ksima iteracin del Mtodo de Newton
T Matriz de transformacin de marco de referencia
Incremento
Sumatorio
pi
ngulo de desfasaje de onda de tensin
G
c
Conductancia representativa de las
prdidas del STATCOM
K Constante
v
m1
, v
m2
v
d1
, v
d2
Variables de estado del STATCOM
B Suceptancia
Velocidad o frecuencia angular
Captulo 1
Introduccin
En ste captulo se presenta una breve resea histrica sobre los dispositivos FACTS.
Tambin, se citan ejemplos de equipamientos FACTS. Adems de la evolucin histrica
del STATCOM, se presenta la motivacin y los objetivos de sta tesis.
1
1.1 Introduccin 2
1.1. Introduccin
La industria de los semiconductores de potencia ha producido dispositivos de conmu-
tacin que pueden operar con lmites de tensin y corriente del orden de los kV y kA
respectivamente. Esto permite la aplicacin de la electrnica de potencia en equipos
de industrias, y en los sistemas de distribucin y transmisin de energa. Adems tales
equipos deben soportar fallas de decenas a centenas de MVAs.
La electrnica de potencia puede auxiliar el funcionamiento de los dispositivos utiliza-
dos en la industria, sustituyendo los metodos convencionales, proporcionando formas
mas rpidas de control, mas ecientes y econmicas de utilizacin de la energa.
En los sistemas de potencia, la utilizacin de la electrnica de potencia tiene sus ini-
cios en los proyectos de transmisin de corriente continua a alta tensin, mas conoci-
dos como HVDC, High Voltage Direct Current. stos sistemas utilizan conversores de
corriente alterna a corriente directa y vice-versa controlados por tiristores.
El control de los elementos semiconductores permite la construccin de equipamien-
tos capaces de realizar alteraciones rpidas de las caractersticas del sistema elctrico.
stos dispositivos operan dentro del concepto conocido como sistemas de tranmisin
exibles en corriente alterna FACTS(Flexible Alternating Current Transmission Sys-
tem), creado por Hingorani [N.G.Hingorani, 1988; N.G.Hingorani, 1990].
El STATCOM, sigla de Static Synchronous Compensator, es un dispositivo FACTS,
compuesto por inversores, transformadores, capacitores y controles, como se muestra
en la gura 1.1, el cual se implementa para controlar potencia reactiva, tanto capacitiva
como inductiva. ste dispositivo puede ser utilizado para hacer el control del factor de
potencia, regulacin de tensin en el punto de conexin y para mejorar la estabilidad
dinmica del sistema elctrico [Gyugyi, 1994].
El STATCOM a pesar de ser un dispositivo relativamente nuevo, ya posee una exten-
sa literatura. No obstante es importante saber y constatar con modelos adecuados los
alcances que pueda tener el STATCOM al incluirlo en un sistema de potencia.
1.2. Historia de los dispositivos FACTS
Los equipamientos bajo el concepto de FACTS comenzaron a ser implementados a
nales de la decada de los 60 por diversos grupos, y bajo el nombre de FACTS, Flexible
1.3 Tipos de dispositivos FACTS 3
Control del
Statcom
Inversor
Mediciones de corrientes
Mediciones
de tension
Capacitor
Transformador de
Acople
STATCOM
Figura 1.1: Forma bsica del Statcom
Alternating Current Transmission System, empez a existir a partir de 1988 cuando
Hingorani public sus artculos[N.G.Hingorani, 1988; N.G.Hingorani, 1990]. En esta
poca Hingorani vice-presidente del EPRI(Energy Power Research Institute), diriga
un grupo de trabajo que investigaba sobre el empleo de la electrnica de potencia en
los sistemas de potencia.
El concepto de FACTS agrupa a un conjunto de nuevos equipos de electrnica de po-
tencia que permiten mayor exibilidad al control de los sistemas elctricos. En este
caso entiendase exibilidad, como la capacidad de alteracin rpida y continua de los
parmetros que controlan la dinmica de funcionamiento de los sistemas elctricos.
1.3. Tipos de dispositivos FACTS
Los dispositivos FACTS se dividen en dos grandes grupos: los compensadores conec-
tados en paralelo y los compensadores conectados en serie.
A travs de su historia los dispositivos FACTS se han clasicado por generaciones. La
primera generacin son dispositivos FACTS que emplean tiristores y elementos semi-
conductores descritos en [N. Mohan, 1995].
Ejemplos de dispositivos conectados en paralelo con la red del sistema de potencia
pertenecientes a la primera generacin son: el TCR(Thyristor Controlled Reactor), el
TSC(Thyristor Switched Capacitor) y el SVC(Static Var Compensator). Diagramas de
estos dispositivos, se muestran en la gura 1.2.
Ejemplos de dispositivos conectados en serie con la red y pertenecientes a la primera
generacin son: el TSSC(Thyristor Switched Series Capacitor) y el TCSC (Thyristor
Controlled Series Capacitor) cuyos esquemas se muestran en la gura1.3.
La segunda generacin est compuesta por dispositivos que utilizan elementos semi-
conductores tales como los IGBTs (Insulated Gate Bipolar Transistor) y GTOs (Gate
Turn Off Thyristor), los cuales se describen en [N. Mohan, 1995].
Dispositivos FACTS que utilizan estos elementos de conmutacin como ejemplos de
compensadores paralelos son el STATCOM (Static Syncronnous Shunt Compensator),
cuyo esquema est mostrado en la gura1.4, y compensadores conectados en serie
como el SSSC(Static Syncronnous Series Compensator) cuyo esquema est mostrado
en la gura1.5.
La tercera generacin es compuesta por la integracin de dos dispositivos serie y parale-
lo en una misma lnea de transmisin. Un resultado de esto es el UPFC (Unied Power
Flow Controller), el cual, es un dispositivo combinado de un SSSC y un STATCOM,
conforme es mostrado en la gura 1.6.
La cuarta generacin es conformada por la integracin de dispositivos serie y paralelo
pero en diferentes lineas de transmisin que da lugar al nombre IPFC (Interline Power
Flow Controller) mostrado en la gura 1.7
TCR TSC
SVC
TCR TSC
Figura 1.2: FACTS- primera generacin -paralelo
TSSC
TCSC
C
L/2 L/2
Figura 1.3: FACTS- primera generacin -serie
Conversor
Capacitor
Control del
Statcom
PT
CT
Transformador
de Acople
Figura 1.4: FACTS- segunda generacin -paralelo-STATCOM
Conversor
Capacitor
Control del
SSSC
PT
CT
Transformador
de Acople
Figura 1.5: FACTS- Segunda generacin- serie-SSSC
Capacitor
de enlace
Control del
Statcom
Control del
SSSC
Figura 1.6: FACTS-Tercera generacin-UPFC
1.4 Historia del STATCOM 7
Conversor
serie
Conversor
Paralelo
Capacitor
Figura 1.7: FACTS- Cuarta generacin-IPFC
1.4. Historia del STATCOM
El desarrollo del STATCOM comenz en 1976, resultado de un trabajo de investigacion
japones ligado a las compaias Kansay Electric Co. Ltda y Mitsubishi Electric Corpo-
ration [Y. Sumi, 1981] . Este equipamiento utilizaba bloques de inversores compuestos
por tiristores y circuitos auxiliares complejos para la conmutacin forzada. Tena como
objetivo aumentar la capacidad de corriente de compensacin, en contraste con las li-
mitantes del SVC como muestra la gura 1.8 y explicada en las referencias [Narain
G. Hingorani, 2000].
Caracteristica V vs I
SVC
Caracteristica V vs I
STATCOM
Operacion
Transitoria
Operacion
Transitoria
I Ind
I cap I Ind I cap
Operacion
Continua
Operacion
Continua
Vbus
I I
Figura 1.8: Caracteristicas V vs I del SVC y el STATCOM
Tambin el surgimiento del STATCOM prometa para el futuro grandes reducciones
en cuanto al tamao del equipamiento total del dispositivo, eliminando los grandes
reactores y capacitores que poseen los SVCs.
En 1980 el primer prototipo del STATCOM (en esta poca este nombre no era utili-
zado) fue contrudo y probado [Y. Sumi, 1981]. Este fu dimensionado para una po-
tencia de 20MVA a una tensin de 77kV y utilizaba tiristores con circuitos auxiliares
de conmutacin forzada. El objetivo de ste dispositivo era mejorar la dinmica del
sistema elctrico mejorando su estabilidad, a travs de la compensacin de reactivos.
Para esta poca todava las potencias que se manejaban eran relativamente bajas, inclu-
yendo los problemas de complejidad y de conabilidad de los circuitos auxiliares de
conmutacin forzada y de los bajos lmites de corrientes y tensin de los dispositivos
auto-conmutados tipo GTO.
Los lmites existentes de dispositivos de conmutacion en los aos 70s atrasaron en una
dcada el desarrollo de todos los compensadores (SVC, STATCOM, etc.) de estado
slido.
Solamente el rpido avance de los dispositivos auto-conmutados, con mayores capa-
cidades de tensin y corrientes, ocurrido en la dcada de los 80s, fu posible el re-
surgimiento del inters por la construccin del STATCOM. En 1990 haba dispositivos
GTOs de 4,5kV para la tensin de bloqueo y 2,5kA para la corriente de interrupcin,
en el 2000 eran citados nuevos lmites de tensin de bloqueo de 6,5kV y corrientes de
interrupcin de 6,6kA.
Ejemplos de estos nuevos dispositivos de conmutacin son los IEGTs (Injection En-
hancement IGBT,(2,5-6,5kV/2kA), HiGT (High-Conductivity IGBT,(3,3kV/5kA)), y pa-
quetes compactos de IGBTs (4,5kV/2kA), y los GCTs (Gate Commutated Thyristor )
o los IGCTs (Integrated Gate Commutated Thyristor).
Las pruebas del STATCOM de 20 MVA [Y. Sumi, 1981], fueron buenas y el paso a
seguir era la construccin de un dispositivo de mayor potencia.
En 1990, un grupo de investigadores japoneses anunciaba estar bastante avanzado en
el desarrollo de un STATCOM de 80MVA a una tensin de 154kV [Edwards, 1988a].
A travs de pruebas con un modelo reducido de 2 kVA, 200V, 48 pulsos, fueron veri-
cados los sistemas de control y su principio de funcionamiento. Con un modelo parcial
de 12 pulsos fueron vericados y probados los circuitos snubber, como tambin los
dispositivos de conmutacin que fueron empleados.
Al mismo tiempo en Estados Unidos un grupo de investigadores de EPRI y ESERCO
(Empire State Energy Research Corporation), presentaba sus resultados obtenidos de
un modelo de 1 MVA probado en campo [Edwards, 1988b].
En 1992 el primer STATCOM entr en operacin en la subestacin de Inuyama, per-
teneciente a Kansay Electric Corporation Co., Ltd. [S.Mori, 1993]. ste STATCOM
de 48 pulsos posea una potencia de 80 MVA a una tensin de 154 kV y utilizaba
transformadores especiales con conguraciones en zig-zag. Fueron utilizados GTOs
de 4,5kV/3kA en los inversores.
En 1995, un STATCOM americano de 100 MVA relatado por [Narain G. Hingorani,
2000], fu construido en la subestacin de Sullivan, perteneciente a TVA, Tenneesse
Valley Authority. ste STATCOM fu implementado con una estructura de 48 pulsos
para neutralizar armnicos. Esta estructura utilizaba GTOs de 4,4kV/4kA, cinco uni-
dades de GTOs fueron conectados en serie asociados con diodos antiparalelos. En el
lado de corriente continua se coloc un capacitor de 65 kJ y 2948 F, con una tension
nominal de 6,6 kV, este STATCOM estaba conectado a una lnea de 161 kV.
Otro proyecto de construccin fu un UPFC de 320MVA por AEP American Electric
Power, instalado en la subestacin de Inez, Kentucky, relatado en 1998 por [C.Schauder,
1998]. Este proyecto en la fase inicial era un STATCOM de 160MVA conectado a
una lnea de 138 kV. En la segunda etapa otro STATCOM fue insertado de la misma
capacidad que el anterior, permitiendo tener un control de 320MVA. Pero despus,
ste ultimo STATCOM fu modicado por un SSSC, y as obtener la composicin de
un UPFC.
En las referencias de IEEE (Institute of Electrical and Electronics Engineer), el nombre
STATCOM (Static Sycronnous Compensator) ya est establecido, pero los nombres
como SVG(Static Var Generator), ASVC (Advanced Static Var Compensator) y ASVG
(Advanced Static Var Generator), todavia permanecen en Japn y Europa.
Algunos de estos nombres existen como variaciones del nombre STATCOM: el STAT-
CON, derivado de static condenser, D-STATCOM, Distribution-STATCOM fabricado
para sistemas de distribucin y un PWM STATCOM Pulse Width Modulation STAT-
COM, con inversores accionados a altas frecuencias de conmutacin y comportamiento
similar al de un ltro activo.
1.5 Motivacin del trabajo 10
1.5. Motivacin del trabajo
En la actualidad como un nmero de investigaciones del dispositivo STATCOM est
aumentando, es muy importante saber que tanta informacin es valedera para nuestros
futuros estudios en cuanto a la apropiacin de esta tecnologa (FACTS).
El requerimiento de este trabajo por intermedio de la Universidad Industrial de San-
tander de encontrarse a la vanguardia de la investigacin, ha dado origen al presente
proyecto, el cual busca modelar y simular el STATCOM. La investigacin del STAT-
COM permitir la formacin de profesionales con conocimientos en los dispositivos
FACTS y adems de sus mejoras en la insercin en sistemas de energa elctrica.
En el pasado en cuanto al diseo de redes, stas se diseaban de una manera sobre-
dimensionada, para satisfacer las demandas locales. Ahora, con al aumento de la de-
manda de energa y de las interconexiones nacionales e internacionales, las redes de
transmisin se han quedado rezagadas en su capacidad de transporte. La construccin
de nuevos es cada vez mas difcil, debido a sus costos, problemas medioambientales
y de servidumbres. Por tanto, se hace necesario tener unos criterios de eciencia para
un mejor aprovechamiento de las redes, y que mejor tener dentro de esos criterios de
eciencia el uso de tecnologas avanzadas como lo son los dispositivos FACTS.
Con la realizacin de este proyecto por medio de la herramienta de simulacin ATP y
su interfaz grca ATPDraw, se pretende tener un buen recurso informtico que permita
tener una respuesta aceptable de un sistema de potencia junto con el STATCOM.
Finalmente, y no menos importante es el deseo del investigador y del director, pro-
fundizar sus conocimientos en sta rea de investigacin, lo cual, permitir realizar
posteriores investigaciones en ste campo.
1.6. Objetivos
El objetivo fundamental de ste trabajo es la apropiacin de la tecnologa sobre
los sistemas de transmisin exible en corriente alterna FACTS, por intermedio
del modelado y simulacin de uno de sus ms representativos dispositivos como
lo es el STATCOM.
1.6 Objetivos 11
Realizar un modelo detallado del STATCOM para su posterior inclusin en sis-
temas de potencia balanceados, a travs de la implementacin con modelos digi-
tales creados con la herramienta de simulacin ATP-ATPDraw .
Aplicar el modelo digital implementado en ATP, a sistemas electricos de diver-
so tamao, para vericar mejoras en la respuesta dinmica con la inclusin del
STATCOM.
Realizar un modelo reducido del STATCOM para su posterior inclusin en la
herramienta de simulacin de sistemas de potencia ante grandes perturbaciones:
HAGP.
Captulo 2
El STATCOM
En este captulo se describir el principio de funcionamiento del STATCOM como
compensador de potencia reactiva inductiva y capacitiva. Tambin se describen sus
componentes mas importantes y sus topologas ms comunes para su implementacin.
Al nal se hace una breve comparacin entre algunas de dichas topologas.
12
2.1 Funcionamiento del STATCOM 13
2.1. Funcionamiento del STATCOM
El principio de funcionamiento del STATCOM puede ser descrito a travs del diagrama
mostrado en la gura 2.1, que es igual a una representacin del STATCOM en conjunto
con un sistema elctrico de potencia. Dicho sistema, se considera por su equivalente
Thevenin, al igual que el STATCOM, sin la representacin de los armnicos.
V
S V
I
X
L
i
L
V
L
=V
S
V
I
Figura 2.1: Diagrama simplicado de dos fuentes de tensin conectados por una impe-
dancia X
L
En este esquema (g 2.1), el sistema elctrico y el STATCOM, se representan respecti-
vamente por las tensiones Vs y Vi, adems de una reactancia equivalente representando
el transformador de acople.
Las potencia activas y reactivas entre las dos fuentes de tensin en rgimen permanente,
se describen a travs de:
P =
V
s
V
i
X
L
sin (2.1)
Q =
V
2
s
X
L
V
s
V
i
X
L
cos (2.2)
Considerando las ecuaciones (2.1) y (2.2) y los diagramas fasoriales de la gura 2.2, se
destacan las siguientes situaciones:
Cuando la tensin V
i
est adelantada de la tensin V
s
, existe una potencia activa
saliendo del STATCOM hacia la red. P<0, gura 2.2 a).
2.2 Componentes del STATCOM 14
(a) (b) (c)
(d) (e)
V
I V
I
V
I
V
I
V
I
V
S
V
S
V
S
V
S
V
S
V
L
V
L
V
L
V
L
i i
i
i
Figura 2.2: Diagramas fasoriales del sistema elctrico y el STATCOM
Cuando la tensin V
i
est atrasada de la tensin V
s
, exite una potencia activa
saliendo de la red hacia el STATCOM. P>0, gura 2.2 b).
Cuando la tensin V
i
est en fase con la tensin V
s
, y adems |V
i
|=|V
s
|, no existe
ujo de potencia activa ni reactiva, gura 2.2 c).
Cuando la tensin V
i
s
, y adems |V
i
| <|V
s
|, no
existe ujo de potencia activa, pero existe potencia reactiva inductiva, Q>0, el
sistema v al STATCOM como si fuera un inductor, ver gura 2.2 d).
Cuando la tensin V
i
s
, y adems |V
i
| >|V
s
|, no
existe ujo de potencia activa, pero existe potencia reactiva capacitiva, Q<0, el
sistema v al STATCOM como si fuera un capacitor, ver gura 2.2 e).
Vericndose as que una fuente de tensin con capacidad de control de fase puede
direccionar el ujo de potencia activa, y con control de amplitud, puede controlar la
potencia reactiva en aquel punto de conexin.
2.2. Componentes del STATCOM
Un STATCOM tiene su funcionamiento basado en las situaciones expuestas anterior-
mente.
En este caso el STATCOM es una fuente de tensin controlada en amplitud y fase, e
implementada a travs de inversores. Estos inversores poseen en su lado de corriente
continua capacitores para simular una fuente de tensin DC. El capacitor en el lado de
corriente continua tambien sirve como sistema almacenador de energa. Los transfor-
madores son utilizados para compatibilizar los niveles de impedancia del compensador
y del sistema de potencia AC. En la mayora de los casos los transformadores son uti-
lizados para minimizar el contenido de armnicos de tensin del inversor.
As mismo, el STATCOM est compuesto por cuatro partes bsicas: inversores, trans-
formadores, capacitores en el lado de corriente continua y un sistema de control como
se muestra en la gura 1.1.
Aunque existen variaciones de las topologas utilizadas, para el esquema discutido an-
teriormente, el STATCOM contina con el mismo funcionamiento, sin alteraciones.
Las partes integrantes del STATCOM se discuten a continuacin, identicando sus fun-
ciones y variaciones dentro de las topologas investigadas.
2.2.1. Inversores
2.2.1.1. Funcin
Los inversores son un conjunto de vlvulas con capacidad de corte y conduccin con-
troladas, y su funcin en el STATCOM es la generacin de la tensin alterna a partir de
la tensin de corriente continua en el lado de corriente continua del capacitor.
2.2.1.2. Vlvulas utilizadas
Las vlvulas que componen actualmente a los inversores son GTOs (Gate Turn Off
Thyristor), IGBTs (Insulated Gate Bipolar Transistor) e IGCTs (Integrated Gate
Commutated Thyristor). En la actualidad, con mayores capacidades de tensin de blo-
queo (6,5kV) y corriente de corte (6kA), estan siendo fabricados los HiGT (High-
Conductivity IGBT) e IEGTs (Injection Enhancement IGBT).
A pesar de que la capacidad de bloqueo y corte de los IGBTs son un poco menores que
los IGCTs, IEGTs y HiGT pueden soportar frecuencias de conmutacion ms altas. Es
por esta razn que an son ms utilizados.
2.2.1.3. Topologas
Las disposiciones ms comunes de los inversores del STATCOM estn compuestas
por inversores trifsicos de seis pulsos, como en la gura 2.3, inversores monofsicos
agrupados en una estructura trifsica, gura 2.4, y por los inversores de tres niveles tipo
neutral point clamped, o NPC, gura 2.5.
Figura 2.3: Inversor de seis pulsos
A
B
C
Figura 2.4: Inversor trifsico compuesto por 3 inversores puente H monofsicos
El inversor trifsico de seis pulsos, mostrado en la gura 2.3, tiene la ventaja de utilizar
menos vlvulas que los otros ejemplos mostrados. Idealmente se utilizan dos vlvulas
por cada fase, totalizando seis vlvulas.
A
B
C
C1
C2
N
Figura 2.5: Inversor de tres niveles NPC
Los armnicos presentes en la tension fase-fase del inversor de seis pulsos estn en el
orden de 6n 1, donde n= 0,1,2..., eliminando los armnicos
1
multiplos de 3 existentes
en la tensin fase neutro [Narain G. Hingorani, 2000].
Los inversores monofsicos se pueden conectar formando una estructura trifsica como
el mostrado en la gura 2.4.
Comparndolo con el inversor de seis pulsos, la estructura compuesta por inversores
monofsicos utiliza mucho ms vlvulas(12 vlvulas), es decir, el doble de lo que ne-
cesita el inversor de 6 pulsos, el cual, es una limitante importante cuando se desean
implementar inversores de ms pulsos como ocurre con la implemenrtacin del STAT-
COM.
Los inversores trifsicos de tres niveles, como el mostrado en la gura 2.5, han sido es-
tudiados y recomendados, como lo ms reciente en el desarrollo del STATCOM[Narain
G. Hingorani, 2000].
La topologa del neutral point clamped, solo utiliza 12 vlvulas, siendo 4 por cada
fase, adems de 6 diodos. El resultado obtenido de la tensin fase fase es parecido al
obtenido por la topologa trifsica compuesta por inversores monofsicos.
Las desventajas de estas topologas, est en la gran cantidad de vlvulas y diodos a uti-
lizar, requiriendo una lgica de accionamiento ms compleja. Tambin se ha vericado
que existen problemas con la equalizacin de las tensiones de los capacitores en el lado
de corriente continua [S. Kincic, ].
1
El anlisis armnico ser presentado en ms detalle en el siguiente captulo, para topologas del
inversor del STATCOM de seis y doce pulsos
2.2.1.4. Frecuencia de conmutacin
Las tcnicas de accionamiento utilizadas para los inversores se pueden dividir en dos
categoras cuando se aborda la frecuencia de conmutacin: tcnicas de accionamiento
en bajas y altas frecuencias.
Los accionamientos dados en bajas frecuencias son aquellos que utilizan frecuencias
prximas a la frecuencia de la red y se caracterizan por las formas de onda cuadradas,
por las tcnicas PAM, pulse amplitud modulation [Rashid, 2004]. Cuando se utiliza
bajo nmero de pulsos, el contenido armnico en ordenes bajos puede ser alto. Una
ventaja de esta tcnica est ligada a menores prdidas de energa en los elementos de
conmutacin.
Los accionamientos a altas frecuencias, caracterizadas por las tecnicas PWM, utilizan-
do frecuencias de conmutacin mucho mayores que la de la red del orden de algunos
kHz. Los armnicos producidos por estas tcnicas tienen ordenes mas altos que puden
ser ltrados ms facilmente . Tambin se ha vericado que estas tcnicas permiten res-
puestas dinmicas mas rpidas. No obstante con el aumento de los equipamientos, las
prdidas de los dispositivos de conmutacin aumentan, inviabilizando en muchos casos
este tipo de modulacin.
La tcnica de modulacin empleada en este trabajo para la lgica de disparo de los
dispositivos de conmutacin es la llamada SVPWM modulacin por ancho de pulso
del vector espacial, de la cual se hablar en detalle ms adelante.
2.2.2. Capacitor del lado de corriente continua
2.2.2.1. Funcin
La principal funcin del lado de corriente continua del STATCOM, es servir de fuente
de tensin DC, posibilitando la actuacin del inversor. Tambin un condensador del
lado de corriente continua sirve como acumulador temporal de energa, permitiendo su
intercambio entre el sistema elctrico y el STATCOM.
2.2.2.2. Topologas
El arreglo de capacitores ms empleado en el STATCOM, se muestra en la gura 2.6,
donde solo un capacitor se conecta al lado de corriente continua de los inversores.
Tambin, este tipo de topologas se puede extender en conguraciones de hasta 48
pulsos. Existen otros ejemplos, como el mostrado en la gura 2.7, donde para cada
conversor hay un capacitor independiente. En dicha gura se aprecia una sola fase, la
cual, pertenece a una conguracin trifsica .
Figura 2.6: Topologas de capacitores del lado de corriente continua
2.2.3. Transformadores
En el STATCOM son utilizados dos grupos de transformadores para sus dos funciones
principales. El transformador de acoplamiento, cuya funcin es la conexin del STAT-
COM con el sistema elctrico de potencia, a n de atender las limitaciones de tensin
de los dispositivos de conmutacin de los inversores. Y los transformadores de estruc-
tura magntica de reduccin de armnicos, los cuales a travs de formas especiales de
conexin, eliminan algunos de los armnicos generados por los inversores reduciendo
su contenido de inyeccin a la red. Existe la posibilidad de hacer transformadores que
posean ambas funciones, pero ste tipo de arreglo no es comn.
En cuanto al tranformador de acoplamiento, posee relaciones simples de conexin entre
primario y secundario, con conexiones de tipo estrella-estrella estrella-delta.
Cuando los transformadores se utilizan para reducir el contenido de armnicos, existen
dos tipos de arreglos bsicos: las topologas zig-zag y estrella-delta.
Figura 2.7: Topologas con capacitores individuales por cada inversor
2.2.3.1. Transformadores en zig-zag
El arreglo de transformadores llamados zig-zag est compuesto por transformadores
especiales montados en forma tal de poder generar desfasajes en las tensiones. En la
gura 2.8 se muestra un ejemplo de utilizacin de estos transformadores en un STAT-
COM de 12 pulsos.
0,0
30
C
0
30
Transformador
de Acople
Inversores
Desfasaje
de disparos
Figura 2.8: STATCOM de 12 pulsos utilizando transformadores zig-zag
Esquema para atraso de fase
Esquema para avance de fase
V
A11
V
A11 V
A21
V
A21
V
B11
V
B11
V
B21 V
B21
V
C11
V
C11 V
C21
V
C21
V
A12 V
A12
V
A22 V
A22
V
B12 V
B12
V
B22 V
B22
V
C12
V
C12
V
C22 V
C22
xV
a
xV
a
yV
a
yV
a
xV
b
xV
b
yV
b
yV
b
xV
c
xV
c
yV
c
yV
c
Figura 2.9: Esquemas de conexin de los transformadores zig-zag
En este trabajo se dene el primario como el lado de alta, conectado al sistema de
potencia denido por los ndices V
A11
,V
A12
, y el secundario se reere al lado de baja, el
cual, est conectado hacia el lado de los inversores indicados por los ndices V
A21
,V
A22
,
como se muestra en la gura 2.9.
Las conexiones se pueden describir segn las relaciones matemticas :
V
A(adelanto)
= x.V
A
y.V
B
(2.3)
V
B(adelanto)
= x.V
B
y.V
C
(2.4)
V
C(adelanto)
= x.V
C
y.V
A
(2.5)
V
A(atraso)
= x.V
A
y.V
C
(2.6)
V
B(atraso)
= x.V
B
y.V
A
(2.7)
V
C(atraso)
= x.V
C
y.V
B
(2.8)
En estas ecuaciones los valores de x e y son coecientes de proporcin para lograr un
valor de desfasaje , con amplitud de un valor unitario, relacionados de la siguiente
manera:
x = cos
3
3
sen (2.9)
y = 2
3
3
sen (2.10)
Valores de normalmente utilizados en la composicin multipulsos del STATCOM, se
muestran en la tabla 2.1.
Tabla 2.1: Valores de x y de y en funcin de
x y
0
o
1 0
7,5
o
0,916 0,151
15
o
0,816 0,229
22,5
o
0,703 0,442
30
o
0,577 0,577
37,5
o
0,442 0,703
45
o
0,229 0,816
52,5
o
0,151 0,916
El lado de los transformadores que corresponde al sistema AC, los devanados de los
transformadores se conectan en serie, y el lado que corresponde al inveror se conectan
en zig-zag, un ejemplo de stas conguraciones se pueden encontrar en [S.Mori, 1993].
ste arreglo de los transformadores, idealmente tiene la capacidad de hacer que se can-
celen armnicos, permitiendo apenas armnicos del orden de m6n1, donde n=1,2,3,...
y m=1,2,4,8,.... es el nmero de inversores de seis pulsos utilizados.
Por todo lo dicho anteriormente, existe grn desventaja en la utilizacin de este tipo de
transformadores zig-zag.
Debido a la complejidad de este transformador, su costo es mas elevado que el de otros
transformadores ms comunes. El montaje de las relaciones de fases son complejas y la
relacin de espiras necesita de alguna precisin, para obtener los desfasajes deseados y
por consecuencia los armnicos que se desean cancelar. Tambin es una desventaja de
este arreglo la relacin ja que deben tener cada uno de los tranformadores, a menos
que sean disponibles varios taps para lograr la relacin deseada.
0,0
15
C
0
45
Transformador
de Acople
Inversores
Desfasaje
de disparos
30
45
30
15
0,0
15
C
0
45
Transformador
de Acople
Inversores
Desfasaje
de disparos
30
45
30
15
7,5
37,5
22,5
52,5
7,5
22,5
37,5
52,5
Para cada uno de los arreglos multipulsos, mostrados en las guras 2.8, 2.10 y 2.11,
existe un valor de desfasaje especco entre los grupos de transformadores e inversores.
El desfasamiento se realiza a travs del control de disparos de los inversores y por
las relaciones de los transformadores zig-zag. Para obtener una forma de onda de 12
pulsos , es necesario obtener un desfasaje de 30
o
entre cada grupo de transformadores
y un desfasaje de 30
o
entre los pulsos de disparo de los dos inversores. De esta misma
forma se realizan los esquemas del STATCOM de 24 pulsos con desfasajes de 15
o
entre transformadores e inversores y el STATCOM 48 pulsos con desfasajes de 7,5
o
entre transformadores e inversores.
Con el aumento del nmero de pulsos las tensiones generadas por el STATCOM, han
de tener una mayor aproximacin a la forma de onda de la tensin senoidal.
2.2.3.2. Transformadores en estrella - delta
Otra forma de arreglo de los transformadores, es la llamada composicin estrella -
delta. Es una forma aproximada pero su implementacin es ms simple. Como lo dice
su nombre, este arreglo utiliza transformadores estrella-delta y estrella-estrella.
Unas de la ventajas existentes de la utilizacin de este tipo de transformadores es que
sern iguales de 2 en 2 para el STATCOM de 24 pulsos, y de 4 en 4 para el STATCOM
de 48 pulsos (mientras que para el esquema zig-zag se necesitan respectivamente 4 y 8
transformadores diferentes), y no existen complejas relaciones de los arrollamientos de
sus devanados.
Al igual que en la composicin zig-zag las seales de disparo tendrn los mismos des-
fasajes, los devanados del lado del sistema AC estataran en serie y los devanados del
lado del inversor estarn en estrella o en delta, mientras que los desfases entre los trans-
formadores seran de 30
o
(dependiendo de la secuencia de fases). Por estas conside-
raciones se dice de que no existe un cancelamiento total de las componentes armnicas
6n1, al contrario como ocurre con los esquemas zig-zag. Por eso muchos autores ha-
blan para referirse al STATCOM, que utilizan estos esquemas de quasi-24 y quasi-48
pulsos [H. Fujita, 1996] [C.A.C. Cavaliere, 2000].
A manera de ejemplo se presenta un esquema de un STATCOM de quasi- 24 pulsos en
la gura 2.12 .
En los esquemas estrella-delta, tambin es necesario observar que los transformadores
4 INVERSORES
DE SEIS PULSOS
SEALES
DE CONTROL
SISTEMA AC
TRIFASICO
YY
Y
YY
Y
Figura 2.12: Diagrama de un STATCOM quasi-24 pulsos (esquema de conexiones es-
trella - delta)
estrella-delta tienen una relacin de transformacin de
3 entre primario y secundario

y esto debe ser tenido en cuenta en las relaciones de tensin.
Captulo 3
Anlisis detallado del STATCOM de 6
y 12 pulsos
En este captulo se hace unlisis del STATCOM de seis y doce pulsos, se hallan las
expresiones matemticas de las tensiones y corrientes del STATCOM, para tener una
mejor comprensin acerca del intercambio de potencia rectiva entre el sistema AC y el
compensador. Adems se hace un anlis de Fourier, para saber que tipo de armnicos
aparecen en corriente y en tensin en la red de potencia.
26
3.1 Introduccin 27
3.1. Introduccin
En este captulo se presenta un anlisis detallado de un STATCOM de 6 y 12 pulsos. El
anlisis conlleva a tener una idea ms clara de como ocurre el intercambio de potencia
con el inversor, utilizando el capacitor como una fuente DC y de los modos de operacin
del STATCOM. Adems se ofrece una importante explicacin acerca de la disposicin
de un STATCOM multipulso.
3.2. Fuente inversora de Tensin (FIT)
La fuente inversora de tensin tal como se muestra en la gura 3.1, consta de una
conguracin trifsica, comunmente llamada inversor de seis pulsos, colocados parale-
lamente en tres ramas de a dos dispositivos de conmutacin (GTOs,IGBTs), por cada
rama, acompaados de sus respectivos diodos antiparalelos.
Las seales de control de disparos g
i
, aplicadas a los transistores se controlan tal que
cada una conduce sobre 180
o
, cuando el inversor se conecta a una carga resistiva. La
gura 3.2 muestra las seales de disparo de los transistores. El inversor se puede ver
como la combinacin de tres inversores monofsicos donde cada brazo produce una
seal desfasada 120
o
con respecto a las otras. Las seales de disparo estn desfasadas
60
o
entre s.
Vdc
Q1 Q3
Q4
Q6
Q2
D1 D3
D2
D6
D4
g1
g2
g3
g4
g6
R
R R
ia ib
ic
Q5
D5
g5
Figura 3.1: Fuente inversora de tensin
3.2 Fuente inversora de Tensin (FIT) 28
g1
g2
g3
g4
g5
g6
0
3
2
3

2 3 4 5
t
Figura 3.2: Seales de disparo
0
Vab
Vbc
Vca
-Vdc
-Vdc
-Vdc
+Vdc
+Vdc
+Vdc
3
2
3

2 3 4 5
t
Figura 3.3: Formas de onda de las tensiones lnea a lnea
La secuencia de conmutacin mostrada en la gura 3.2 genera las tensiones de lnea
Van
0
1/3 Vdc
2/3 Vdc
-1/3 Vdc
-2/3 Vdc
0
Vbn
Vcn
3
2
3

2 3 4 5
t
Figura 3.4: Formas de onda de las tensiones lnea a neutro
v
ab
(t),v
bc
(t) y v
ca
(t) ilustradas en la gura 3.3. En esta gura se nota que las tensiones
de lnea tienen un ancho de pulso de 120
o
con una magnitud pico de tensin, V
DC
. Por
cada intervalo de 60
o
existen tres diferentes modos de operacin [N. Mohan, 1995] ;
en el primer medio ciclo la secuencia de operacin 1-5-6 (g 3.5), 1-2-6 (g 3.6) y
1-2-3 (g 3.7). En orden para obtener las tensiones lnea a neutro, sern analizados los
diferentes modos de operacin correspondientes.
Secuencia 1-5-6; 0 t /3
v
an
= v
cn
=
1
3
V
DC
v
bn
=
2
3
V
DC
Secuencia 1-2-6; /3 t 2/3
v
an
=
2
3
V
DC
v
bn
= v
cn
=
1
3
V
DC
Secuencia 1-2-3; 2/3 t
v
an
= v
bn
=
1
3
V
DC
v
cn
=
2
3
V
DC
a
c
b
+V
an
+V
cn
+V
bn
n
V
dc
Figura 3.5: Circuito equivalente para la secuencia 1-5-6
a
c
b
+V
an
+V
cn
+V
bn
n
V
dc
a
c
b
+V
an
+V
cn
+V
bn
n
V
dc
3.3 Anlisis armnico 31
Los siguientes tres modos de operacin se obtienen en forma similar, donde ellos son
los negativos de los tres modos anteriores. La gura 3.4 muestra las tensiones lnea a
neutro v
an
(t), v
bn
(t) y v
cn
(t). La sincronizacin y la frecuencia de las tensiones genera-
das dependen directamente de la sincronizacin y frecuencia de las seales de disparo
y no del tipo de carga. La magnitud pico de tensin depende de la tensin DC.
3.3. Anlisis armnico
El contenido armnico de las tensiones lnea a lnea y lnea a neutro, se pueden obtener
aplicando un anlisis de Fourier para las formas de onda mostradas en las guras 3.3
y 3.4. La simetra impar de la tensin v
ab
(t) se obtiene atrasandola 30
o
, adems los
valores instantneos de v
ab
(t) basados en el anlisis de Fourier es dado por:
v
ab
(t) =
n=1
V
ab
n
sin(nt +

6
n) (3.1)
Tomando las indicaciones de [William H. Hayt Jr., 1998] acerca de la simetra impar
de la seal se puede hallar:
V
ab
n
(t) =
4
T
Z 5
6
6
V
DC
sin(nt)dt (3.2)
con T=2
as que:
V
ab
n
=
2
n
V
DC
(cos(
6
n) cos(
5
6
n)) (3.3)
De la ecuacin (3.3) se nota que slo los terminos 6r1 (n=6r1), donde r, son to-
dos los enteros positivos, tal que n= 1,5,7,11,13..., adems la ecuacin (3.3) se puede
reducir a:
V
ab
n
=
4
n
V
DC
cos(
6
n) n = 6r 1, r = 0, 1, 2, ... (3.4)
Finalmente la tensin v
ab
(t), se expresa por la ecuacin (3.5). Similarmente las tensio-
nes v
bc
(t) y v
ca
(t) se obtienen de la misma manera:
3.3 Anlisis armnico 32
v
ab
(t) =
4
V
DC
n
1
n
cos(
6
n)sin(nt +

6
n) (3.5)
v
bc
(t) =
4
V
DC
n
1
n
cos(
6
n)sin(nt
2
n) (3.6)
v
ca
(t) =
4
V
DC
n
1
n
cos(
6
n)sin(nt
7
6
n) (3.7)
n = 6r 1 r = 0, 1, 2, 3, ...
De una manera similar se obtiene el contenido armnico de las tensiones lnea a neutro
v
an
(t), v
bn
(t) y v
cn
(t). Analizando v
an
(t) se nota que tiene simetra impar, as que
v
an
(t) =
n=1
V
an
n
sin(nt)
donde:
V
an
n
=
4
3T
V
DC
_
Z
/3
0
sin(nt)dt +2
Z
2/3
/3
sin(nt)dt +
Z

2/3
sin(nt)dt
_
(3.8)
y T=2.
Por consiguiente:
V
an
n
=
2
3n
V
DC
_
cos
_
3
n
_
cos
_
2
3
n
_
+1(1)
n
_
(3.9)
Aunque las formas de onda de las tensiones lnea a lnea y lnea a neutro, son di-
ferentes, ellas presentan un contenido armnico similar. Las tensiones lnea a neutro
tambin poseen armnicos del orden n=6r 1, donde r son los enteros positivos tal que
n=1,5,7,11,13...., con relacin a esto, la ecuacin (3.9) se puede expresar por:
V
an
n
=
4
3n
V
DC
_
cos
_
3
n
_
+1
_
n=6r 1, r = 0, 1, 2, ... (3.10)
Las siguientes ecuaciones describen las tensiones lnea a neutro basado en el anlisis
3.4 Compensador esttico sncrono (STATCOM) basado en un inversor de seis
pulsos 33
de Fourier
v
an
(t) =
4
3
V
DC
n
1
n
_
cos
_
3
n
_
+1
_
sin(nt) (3.11)
v
bn
(t) =
4
3
V
DC
n
1
n
_
cos
_
3
n
_
+1
_
sin
_
nt
2
3
_
(3.12)
v
cn
(t) =
4
3
V
DC
n
1
n
_
cos
_
3
n
_
+1
_
sin
_
nt +
2
3
_
(3.13)
n = 6r 1 r = 0, 1, 2, 3, ...
De estas ecuaciones se puede notar que las tensiones de fase y de lnea estn desfasadas
por 30
o
, y que las tensiones de lnea son
3 veces las tensiones de fase.

3.4. Compensador esttico sncrono (STATCOM) basa-
do en un inversor de seis pulsos
Los compensadores estticos son dispositivos con la habilidad de generar y absorber
potencia reactiva y activa, pero la aplicacin mas comn es el intercambio de potencia
reactiva entre el sistema AC y el compensador.
El compensador esttico sncrono (STATCOM) basado en un inversor de seis pulsos
como el que se muestra en la gura 3.8, es un bloque de construccin bsica para
los compensadores de gran potencia. Aunque en aplicaciones de gran potencia esta
conguracin no es muy utilizada debido al gran nmero de armnicos que genera.
Consecuentemente se le hace un anlisis en corriente, para tener una base y comprender
mejor las conguraciones de inversores de ms pulsos.
3.4.1. Intercambio de potencia reactiva
Como se explic anteriormente el intercambio de potencia reactiva entre el sistema
AC y el compensador se controla por la variacin de la magnitud de la componente
fundamental de la tensin del inversor, por encima o por debajo de la tensin del sistema
AC.
El control del compensador se lleva a cabo por pequeas variaciones del ngulo de dis-
paro de los dispositivos semiconductores, a n de que la componente fundamental de la
pulsos 34
Vdc
Q1 Q3
Q4
Q6
Q2
D1 D3
D2
D6
D4
g1
g2
g3
g4
g6
L L
ia ib
ic
Q5
D5
g5
L
Figura 3.8: STATCOM con inversor de seis pulsos
tensin producida por el inversor se obligue a adelantar o atrasar la tensin del sistema
AC por unos pocos grados. Esto causa un ujo de potencia activa entrando o salien-
do del inversor modicando el valor de tensin DC del capacitor, y consecuentemente
la magnitud de la tensin en terminales del inversor, resultando un ujo de potencia
reactiva. Si el compensador solo suministra potencia reactiva, la potencia reactiva su-
ministrada por el capacitor ser cero y adems su tensin no cambiar.
3.4.1.1. Anlisis de las seales de corrientes AC
El ujo de corriente entre el compensador y el sistema AC se determina por la tensin
que aparece a travs del inductor. Tomando la tensin del sistema AC puramente se-
noidal e
an
(t) =V
m
sen(t), entonces la magnitud de la componente fundamental de la
tensin del inversor dada por la ecuacin (3.13) v
an
(t) = 0,6366V
DC
sen(t), entonces,
i
a
(t)
1
=
V
m
0,6366V
DC
L
cos(t); (componente fundamental)
i
a
(t)
n
=
0,6366V
DC
n
2
L
cos(nt); (componente del armnico n)
la componente fundamental de la corriente tambin se llama componente fundamental
I
q
.
I
a1
= I
q
=
V
m
0,6366V
DC
L
(3.14)
pulsos 35
I
an
= I
qn
=
0,6366V
DC
n
2
L
n > 1 (3.15)
Respecto a lo anterior la componente fundamental de la corriente estar en adelanto
cuando V
m
<0,6366V
DC
. Esto es, si la amplitud de la tensin del inversor se incrementa
por encima de la tensin del sistema AC (la corriente uye del conversor al sistema
AC); en este caso el compensador se ve como un capacitor desde el sistema AC. La
componente fundamental de la corriente estar atrasada cuando V
m
> 0,6366V
DC
, esto
es, si la amplitud de la tensin del inversor est por debajo de la del sistema AC (la
corriente uye del sistema AC al conversor). En este caso el compensador se ve como
un inductor desde el sistema AC.
De la ecuacin (3.15) se puede deducir que los armnicos solo uyen del compensador
al sistema AC.
Las ecuaciones de i
a
(t), se derivan teniendo en cuenta la tensin a travs de los ter-
minales del inductor sobre cada intervalo de conduccin de 60
o
(g 3.9). Las otras dos
corrientes de fase i
b
(t) e i
c
(t) sern idnticas, excepto desfasadas por 120
o
y 240
o
res-
pectivamente.
Vm
0
1/3 Vdc
2/3 Vdc
-1/3 Vdc
-2/3 Vdc
0
3
2
3
2 t
Figura 3.9: Tensiones de fase del compensador y del sistema AC
En todo tiempo la tensin v
L
(t) = L
d
dt
i
a
(t), donde v
L
(t) es la tensin instantnea del
inductor, la cual , es la diferencia de las tensiones instantneas del sistema AC, e
an
(t), y
pulsos 36
el STATCOM , v
an
(t). Las siguientes ecuaciones pueden ser escritas para cada intervalo
de conduccin de 60
o
.
Intervalo: 0 t /3
v
L
(t) =V
m
sen(t)
1
3
V
DC
= L
d
dt
i
a
(t)
integrando a ambos lados se tiene:
i
a
(t) =
V
m
L
(cos(t) 1)
1
3L
V
DC
t +I
0
(3.16)
donde, I
0
es la condicin inicial en t = 0; i
a
(0) = I
0
.
Intervalo: /3 t 2/3
v
L
(t) =V
m
sen(t)
2
3
V
DC
= L
d
dt
i
a
(t)
i
a
(t) =
V
m
L
_
cos(t)
1
2
_
_
2
3L
t
2
9L
_
V
DC
+I
1
(3.17)
donde, I
1
= i
a
(

3
)
intervalo 2/3 t
v
L
(t) =V
m
sen(t)
1
3
V
DC
= L
d
dt
i
a
(t)
integrando a ambos lados en ese intervalo:
i
a
(t) =
V
m
L
_
cos(t)
1
2
_
_
1
3L
t
2
9L
_
V
DC
+I
2
(3.18)
donde, I
2
= i
a
(
2
3
)
y es conocido que en estado estable la forma de onda de la corriente es
simtrica, adems:
i
a
(

2
) = 0
pulsos 37
Teniendo en cuenta la propiedad mostrada arriba, la condicin en estado estable I
0
, se
puede calcular como sigue, remplazando en la ecuacin (3.17):
i
a
(

2
) =
V
m
L
_
cos
_
2
_
1
2
_
_
2
3L
.

2
2
9L
_
V
DC
+I
1
= 0
I
1
=
V
m
2L
+

9L
V
DC
(3.19)
Ahora se reemplaza en la ecuacin (3.16),
i
a
(

3
) =
V
m
L
_
cos
_
3
_
1
_
1
3L
V
DC

3
+I
0
= 0
I
1
=
V
m
2L

9L
V
DC
+I
0
(3.20)
de las ecuaciones (3.19) y (3.20), se puede estimar
I
0
=
V
m
L
+
2
9L
V
DC
(3.21)
substituyendo la ecuacin (3.21) en (3.16), se derivan las ecuaciones en estado estable,
i
a
(t) =
V
m
L
cos(t)
_
1
3L
t
2
9L
_
V
DC
(3.22)
0 t /3
i
a
(t) =
V
m
L
cos(t)
_
2
3L
t

3L
_
V
DC
(3.23)
/3 t 2/3
i
a
(t) =
V
m
L
cos(t)
_
1
3L
t

9L
_
V
DC
(3.24)
2/3 t
En el intervalo t 2, las formas de la corriente son negativas respecto a las
corrientes denidas anteriormente en cada intervalo. La gura 3.10 describe las forma
de onda de las corrientes para la fase a. Para el caso de adelanto V
DC
=6V, V
m
=2,5V
y L = 3mH, para el caso de atraso en corriente V
DC
= 6V, V
m
= 4,5V y L = 3mH.
pulsos 38
0 0.01 0.02 0.03 0.04 0.05 0.06
1
0.5
0
0.5
1
Adelanto
Atraso
Figura 3.10: Corrientes entre el sistema AC y el compensador
3.4.1.2. Corriente en el capacitor
La corriente en el capacitor se construye a segmentos de las corrientes del sistema AC
y depende de los dispositivos semiconductores que estn conduciendo en los intervalos
de 60
o
. La corriente en el capacitor se puede deducir de la siguiente forma:
i
D1
(t) +i
D3
(t) +i
D5
(t) +i
DC
(t) = i
Q1
(t) +i
Q3
(t) +i
Q5
(t)
as que:
i
DC
(t) = i
Q1
(t) +i
Q3
(t) +i
Q5
(t) (i
D1
(t) +i
D3
(t) +i
D5
(t)) (3.25)
Intervalo; 0 t /3
i
DC
(t) = i
a
(t) +i
c
(t)
donde:
i
c
(t) =
V
m
L
cos
_
t +
2
3
_
_
1
3L
t +

9L
_
V
DC
asi que:
i
DC
(t) =
V
m
L
sen(t
6
)
_
2
3L
t

9L
_
V
DC
(3.26)
Intervalo; /3 t 2/3
i
DC
(t) = i
a
(t)
i
a
(t) =
V
m
L
cos(t)
_
2
3L
t

3L
_
V
DC
(3.27)
pulsos 39
Intervalo; 2/3 t
i
DC
(t) = i
a
(t) +i
b
(t)
donde:
i
b
(t) =
V
m
L
cos
_
t
2
3
_
_
1
3L
t
4
9L
_
V
DC
asi que:
i
DC
(t) =
V
m
L
sen(t
5
6
)
_
2
3L
t
5
9L
_
V
DC
(3.28)
Estas expresiones son equivalentes excepto cambiadas de fase por 60
o
entre s. Adems,
la forma de onda de la corriente en el capacitor en los tres periodos de conduccin es
idntica a la descrita por la ecuacin (3.26).
3.4.1.3. Tensin DC en el capacitor
El anlisis previo asume una tensin DC constante. Esto es equivalente a considerar
una capacitancia innita y consecuentemente un riple de tensin igual a cero. Si se
tiene una capacitancia nita, entonces se considera que existe un riple de tensin DC,
que depende del valor del capacitor y de la corriente que circule por dicho capacitor.
Asumiendo la corriente dada por la ecuacin (3.26), la tensin en el capacitor sobre los
primeros 60
o
de conduccin est dado por:
v
cap
(t) =
1
C
Z
t
0
i
DC
(t)dt +V
0
(3.29)
donde V
0
es la condicin inicial en t = 0, V
0
= v
cap
(0). Sustituyendo la ecuacin (3.26)
en la ecuacin (3.29) se tiene que:
v
cap
(t) =
V
m
2
LC
cos
_
t
6
_
1
3LC
V
DC
t
2
+

9LC
V
DC
t +
3V
m
2
2
LC
+V
0
(3.30)
El valor de V
0
se calcula de la componente promedio de la ecuacin (3.30) con periodo
T = /3. Con el mismo periodo de tiempo se determina el nivel de tensin V
DC
.
pulsos 40
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018
0.85
0.9
0.95
1
1.05
1.1
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018
0.98
0.99
1
1.01
1.02
1.03
1.04
t(seg)
t(seg)
a)
b)
Figura 3.11: Tensin del capacitor a) generando potencia reactiva; b) absorbiendo po-
tencia reactiva
V
DC
=
1
C
Z
3
0
v
cap
(t)dt (3.31)
V
DC
=
3
V
m
3
LC

3
243
3
LC
V
DC
+

3
162
3
LC
+
3
6
3
LC
V
m
+

3
V
0
_
(3.32)
Simplicando estas expresiones resulta que :
V
0
= 0,0889
_
V
m
2
LC
_
0,0609
_
1
2
LC
V
DC
_
+V
DC
(3.33)
Las formas de onda de corriente y tensin del capacitor obtenidas mediante las ecua-
ciones (3.26) y (3.30) se muestran en las guras 3.11 y 3.12. La gura 3.11 ilustra que
la tensin pico mximo V
pmax
ocurre cuando el compensador est generando potencia
reactiva (corriente en adelanto), t = /3.
V
pmax
= v
cap
_

6
_
=
_
V
m
2
LC
_
_

6
_
2
_
1
3LC
V
DC
_
+
_

9LC
V
DC
__

6
_
pulsos 41
+
3V
m
2
2
LC
+V
0
Simplicando se obtiene:
V
pmax
=0,0451
V
m
2
LC
+0,0305
_
V
DC
2
LC
_
+V
DC
(3.34)
El valor pico de tensin obtenido arriba es importante porque esta tensin pico es apli-
cada directamente sobre los dispositivos semiconductores, los cuales se deben especi-
car para soportar como mnimo este valor de tensin en estado estable.
Figura 3.12: Corrientes del capacitor a) generando potencia reactiva; b) absorbiendo
potencia reactiva
3.4.2. Intercambio de potencia activa y reactiva
Considerando una diferencia de fase entre el sistema AC y la componente funda-
mental de la tensin del compensador, como se muestra en la gura 3.13, permite el
intercambio de potencia activa y reactiva.
De manera similar para obtener las ecuaciones de corriente como en la seccin 3.4.1.1,
se procede de la misma forma:
Intervalo: 0 t /3
pulsos 42
Vm
0
1/3 Vdc
2/3 Vdc
-1/3 Vdc
-2/3 Vdc
0
3
2
3

2 t
Figura 3.13: Tensiones de fase del compensador y del sistema AC
v
L
(t) =V
m
sen(t )
1
3
V
DC
= L
d
dt
i
a
(t)
i
a
(t) =
V
m
L
(cos(t ) cos())
1
3L
V
DC
t +I
0
(3.35)
donde: I
0
es la condicin inicial en t = 0; i
a
(0) = I
0
.
Intervalo: /3 t 2/3
v
L
(t) =V
m
sen(t)
2
3
V
DC
= L
d
dt
i
a
(t)
i
a
(t) =
V
m
L
_
cos(t ) cos
_
__
_
2
3L
t
2
9L
_
V
DC
+I
1
(3.36)
donde: I
1
= i
a
(

3
)
intervalo 2/3 t
pulsos 43
v
L
(t) =V
m
sen(t )
1
3
V
DC
= L
d
dt
i
a
(t)
integrando a ambos lados en ese intervalo:
i
a
(t) =
V
m
L
_
cos(t ) cos
_
2
3

__
_
1
3L
t
2
9L
_
V
DC
+I
2
(3.37)
donde: I
2
= i
a
(
2
3
)
y adems teniendo en cuenta que i
a
(0) =i
a
(
), en estado estable la con-

dicin inicial se puede calcular como sigue:
De la ecuacin (3.37) se obtiene
I
0
=i
a
(
) =
V
m
L
_
cos() cos
_
2
3

__
+
_
1
3L
.
2
9L
_
V
DC
I
2
as que:
I
0
=
V
m
L
_
cos() cos
_
2
3

__
+

9L
V
DC
I
2
(3.38)
Ahora aplicando las condiciones iniciales para sus respectivas ecuaciones:
I
2
= i
a
(
2
3
) =
V
m
L
_
cos
_
2
3

_
cos
_
2
3

__
2
9L
V
DC
+I
1
(3.39)
I
1
= i
a
(

3
) =
V
m
L
_
cos
_
_
cos()
_

9L
V
DC
+I
0
(3.40)
Sustituyendo las ecuaciones (3.39) y (3.40) en la ecuacin (3.38) se puede estimar que:
I
0
=
V
m
L
cos() +
2
9L
V
DC
(3.41)
Por consiguiente de la ecuacin (3.41) se derivan las ecuaciones en estado estable,
i
a
(t) =
V
m
L
cos(t )
_
1
3L
t
2
9L
_
V
DC
(3.42)
0 t /3
3.5 Conversor de 12 pulsos 44
i
a
(t) =
V
m
L
cos(t )
_
2
3L
t

L
_
V
DC
(3.43)
/3 t 2/3
i
a
(t) =
V
m
L
cos(t )
_
1
3L
t

9L
_
V
DC
(3.44)
2/3 t
3.5. Conversor de 12 pulsos
El conversor de seis pulsos en muchas aplicaciones de gran potencia no es muy em-
pleado debido al gran contenido de armnicos que se generan. La combinacin de dos
conversores de seis pulsos ofrece mejor desempeo. Esta nueva conguracin se llama
STATCOM de doce pulsos.
Las tensin lnea a lnea, v
ab
(t), del conversor de seis pulsos, dada por la ecuacin (3.1),
se expresa a continuacin:
v
ab
(t) =V
ab
1
sen(t +30
o
) +V
ab
5
sen(5t +150
o
) +V
ab
7
sen(7t +210
o
)
+V
ab
11
sen(11t +330
o
) +V
ab
13
sen(13t +30
o
) +V
ab
17
sen(17t +150
o
)
+V
ab
19
sen(19t +210
o
) +V
ab
23
sen(23t +330
o
) +.... (3.45)
Si el compensador se conecta a un transformador Y Y, con relacin de transformacin
1 : 1, la tensin lnea a neutro estar dada por:
v
an
(t) =
1
3
(V
ab
1
sen(t) V
ab
5
sen(5t) V
ab
7
sen(7t) +V
ab
11
sen(11t)
+V
ab
13
sen(13t) V
ab
17
sen(17t) V
ab
19
sen(19t) +V
ab
23
sen(23t) +....) (3.46)
v
an
(t) =
1
3
n=1
V
ab
n
(1)
r
sen(nt) (3.47)
n = 6r 1, r = 0, 1, 2, ...
Ahora supngase que la tensin lnea a lnea producida por el segundo conversor est
atrasada 30
o
con respecto al otro conversor y con la misma magnitud. As que:
v
ab
(t)
2
=V
ab
1
sen(t) +V
ab
5
sen(5t) +V
ab
7
sen(7t) +V
ab
11
sen(11t)
+V
ab
13
sen(13t) +V
ab
17
sen(17t) +V
ab
19
sen(19t) +V
ab
23
sen(23t) +.... (3.48)
v
ab
(t)
2
=
n=1
V
ab
n
sen(nt) (3.49)
Si el segundo conversor se conecta a un transformador Y, con relacin de transfor-
macin 1 :
1
3
, la tensin lnea a neutro en el devanado secundario en Y debe ser:
v
anY
(t)
2
=
1
3
(V
ab
1
sen(t) +V
ab
5
sen(5t) +V
ab
7
sen(7t) +V
ab
11
sen(11t)
+V
ab
13
sen(13t) +V
ab
17
sen(17t) +V
ab
19
sen(19t) +V
ab
23
sen(23t) +....) (3.50)
v
anY
(t)
2
=
1
n=1
V
ab
n
sen(nt) (3.51)
n = 6r 1, r = 0, 1, 2, ...
Entonces, la tensin lnea a lnea en el devanado conectado en estrella es:
v
abY
(t)
2
=V
ab
1
sen(t +30
o
) V
ab
5
sen(5t +150
o
) V
ab
7
sen(7t +210
o
)
+V
ab
11
sen(11t +330
o
) +V
ab
13
sen(13t +30
o
) V
ab
17
sen(17t +150
o
)
V
ab
19
sen(19t +210
o
) +V
ab
23
sen(23t +330
o
) +....... (3.52)
v
abY
(t)
2
=
n=1
V
ab
n
(1)
r
sen(nt +

6
n) (3.53)
n = 6r 1, r = 0, 1, 2, ...
Las dos formas de onda dadas por las ecuaciones (3.1) y (3.53), se suman para dar
origen a la tensin de 12 pulsos,
v
ab
(t)
12
= v
ab
(t) +v
abY
(t)
2
v
ab
(t)
12
= 2(V
ab
1
sen(t +30
o
) +V
ab
11
sen(11t +330
o
) +V
ab
13
sen(13t +30
o
)
+V
ab
23
sen(23t +330
o
) +....) (3.54)
Estas formas de onda se muestran en las gura 3.14 y el arreglo de dos conversores
de seis pulsos conectados en paralelo a un mismo condensador se muestra en la gura
3.15.
Vab
0
0.5774 Vdc
1.1547Vdc
-0.5774Vdc
-1.1547Vdc
0
a)
Vdc
-Vdc
b)
f1 f2
f3
f4 f5 f6 f7
f8
Figura 3.14: a) v
ab
(t) y v
abY
(t)
2
; b) Tensin del conversor de 12 pulsos
ia(t)
ib(t)
ic(t)
C
+
-
Vdc
idc(t)
Y-Y
Y-
Conversor de
seis pulsos
Conversor de
seis pulsos
Figura 3.15: STATCOM de 12 pulsos
La tensin de lnea dada por la ecuacin (3.54) se expresa en series de Fourier como:
v
ab
(t)
12
=
n=1
V
ab12
n
sen(nt +

6
n)
n = 12r 1, r = 0, 1, 2, ...
3.5.1. Seales de corriente AC
Aplicando un procedimiento similar al presentado en secciones anteriores, se llega al
anlisis de corrientes AC. De antemano se puede ver de la ecuacin (3.54) que solo apa-
recen las componetes armnicas n= 12r1, (r=0,1,2,...) de tensiones de lnea multipli-
cadas por dos. Entonces se puede decir que la magnitud de la componente fundamental
para las tensiones de lnea del STATCOM de 12 pulsos es el doble de la magnitud de
la tensin de lnea del STATCOM de 6 pulsos (lo mismo ocurre para las tensiones de
fase).
Suponiendo una tensin para el sistema AC puramente senoidal e
an
(t) = V
m
sen(t),
y de lo dicho anteriormente, la magnitud de la componente fundamental de la tensin
lnea-neutro del inversor viene dada por la ecuacin v
an
(t) = 1,2732V
DC
sen(t), en-
tonces:
i
a
(t)
1
=
V
m
1,2732V
DC
L
cos(t); (componente fundamental)
i
a
(t)
n
=
1,2732V
DC
n
2
L
cos(nt); (componente del armnico n)
la componente fundamental de la corriente tambin se llama componente fundamental
I
q
.
I
a1
= I
q
=
V
m
1,2732V
DC
L
(3.55)
I
an
= I
qn
=
1,2732V
DC
n
2
L
n > 1 (3.56)
Respecto a lo anterior la componente fundamental de la corriente estar en adelanto
cuando V
m
< 1,2732V
DC
. Esto es, si la amplitud de la tensin del inversor se incre-
menta por encima de la tensin del sistema AC (la corriente uye del conversor al
sistema AC). En este caso el compensador se como un capacitor desde el sistema AC.
La componente fundamental de la corriente estar atrasada cuando V
m
> 1,2732V
DC
.
Esto es, si la amplitud de la tensin del inversor est por debajo de la del sistema AC (la
corriente uye del sistema AC al conversor). En este caso, el compensador se ve como
un inductor desde el sistema AC.
Captulo 4
Control vectorial de inversores
trifsicos
En ste captulo se estudia la tcnica de modulacin por ancho de pulso del vector
espacial (SVPWM), generalmente llamada control vectorial. El control vectorial SV-
PWM de inversores trifsicos alimentados por fuente de tensin (VSI o voltage source
inverter) constituye una poderosa herramienta para el desarrollo y control de los con-
vertidores de continua-alterna.
49
4.1 Introduccin 50
4.1. Introduccin
Existe una gran variedad de tcnicas de modulacin para convertidores de tensin conti-
nua a alterna (DC/AC), comunmente denominados inversores. La modulacin de 120
o
,
la modulacin de 180
o
, la modulacin por desplazamiento de fase y las modulaciones
de la anchura del pulso, son algunas de las tecnicas de modulacin ms conocidas. Ca-
da una de ellas con un concepto y funcionamiento especco y determinado [Rashid,
2004]. La eleccin de una tcnica de modulacin u otra depende de los factores que
caracterizan la aplicacin, como son entre otros:
Nivel de potencia a controlar.
Dispositivos semiconductores empleados.
Requisitos de la carga.
Caractersticas de la forma de onda de salida (Distorsin, amplitud, frecuen-
cia,...).
Una de las tcnicas ms utilizadas es la modulacin por ancho de pulso (PWM) con
sus multiples variantes: un nico pulso por semiperiodo, varios pulsos por semiperiodo,
varios pulsos con referencia senoidal, etc [Rashid, 2004] . Tiene como principal ventaja
la sencillez en su implementacin, pero como contrapartida, no se adaptan bien a la
respuesta dinmica en cargas no lineales y poseen elevadas prdidas por conmutacin
que se producen cuando se requiere una distorsin armnica total (THD) de bajo valor,
debido a que el nmero de conmutaciones debe ser elevado.
En este captulo se estudia la tcnica de modulacin por ancho de pulso del vector
espacial (SVPWM), generalmente llamada control vectorial. Su principal caracterstica
es que se sustituye todo el sistema trifsico por un slo vector en el que la frecuencia
queda reejada en su velocidad de giro con el paso del tiempo. Esto permite emplear
dicho vector para estudiar tanto los regmenes estacionarios como dinmicos en dichos
sistemas.
Son varias e importantes las ventajas que presenta el control vectorial aunque depende
de la aplicacin. Por ejemplo, cuando se trata de controlar un motor trifsico de induc-
cin, el control vectorial sobre el inversor trifsico permite desacoplar las variables del
4.2 Representacin de seales elctricas en el espacio vectorial 51
motor de induccin de forma que se logra un control independiente de la velocidad y
del par, equiparable al motor de continua [Vas, 1994]. O bien cuando lo que se pretende
es generar un sistema trifsico de tensiones senoidales (por ejemplo en sistemas de ge-
neracin de energa elctrica a partir de paneles solares y bateras), la implementacin
digital del control vectorial resulta mucho ms sencilla que la de un PWM senoidal
trifsico.
El Control Vectorial SVPWM de inversores trifsicos alimentados por fuente de ten-
sin (VSI voltage source inverter ) constituye una poderosa herramienta para el desa-
rrollo y control de los convertidores de continua-alterna.
4.2. Representacin de seales elctricas en el espacio
vectorial
La representacin vectorial con vectores espaciales se basa en la forma en que se pue-
de representar todo sistema compuesto por tres magnitudes x
1
, x
2
y x
3
, dependientes
del tiempo, y que en un marco de referencia bidimensional verican la relacin de la
ecuacin (4.1).
x
1
(t) +x
2
(t) +x
3
(t) = 0 (4.1)
En nuestro caso, en los sistemas elctricos trifsicos equilibrados, estas tres magnitudes
representan corrientes, tensiones, etc.
4.2.1. Sistemas trifsicos a cuatro hilos
En un sistema trifsico a cuatro hilos, las tensiones de las tres fases respecto al neutro
verican la relacin de la ecuacion (4.1), pudiendose expresar de la forma siguiente:
v
R
(t) +v
S
(t) +v
T
(t) = 0 (4.2)
Las tres tensiones de un sistema trifsico equilibrado de cuatro hilos (tres fases simples)
se puede representar en un marco de referencia bidimencional mediante tres
vectores espaciales desfasados entre si 120
o
tal y como se muestra en la gura 4.1.
Vr
Vt
Vs
120
Vr Vt + = -Vs
Figura 4.1: Tensiones de fase de un sistema trifsico
En la gura 4.1 se observa que sobre el plano bidimensional se han distribuido
las tres componentes de un sistema trifsico
V
R
,
V
S
,
V
T
,
desfasadas entre si 2/3 (120
o
), haciendo coincidir la primera de las componentes,
V
R
con el eje de las abscisas para tener un desfasaje de 0
o
con el eje de referencia y, de esta
manera, simplicar las ecuaciones que describen a cada una de las tensiones simples:
v
R
(t) =V
m
sen(t)
v
S
(t) =V
m
sen(t
2
3
)
v
T
(t) =V
m
sen(t +
2
3
)
donde V
m
representa el valor mximo de la seal de tensin en este caso, vericando
que:
V
rms
=
V
m
2
La evolucin temporal de las tres tensiones de un sistema trifsico se detallan en la
gura 4.2 .
0 2 4 6 8 10 12
0.5
0
0.5
Vr Vs Vt
t0
Angulo(rad)
t1 t2
Figura 4.2: Evolucin temporal de las tensiones de un sistema trifsico
t=t0
120 120
t=t1
120
t=t2
Vr
Vr
Vr
Vs
Vt
Vt
Vt
Vs
Vs
W
W
W
V(t)
V(t)
V(t)
Figura 4.3: Evolucin temporal del vector suma
Sobre la graca de la gura 4.2 se han marcado tres instantes determinados (t
0
, t
1
y t
2
)en
los que se tiene en cada uno de ellos, una terna de valores diferentes de las tensiones
intantneas. La gura 4.3 muestra en estos tres instantes las tres tensiones simples en
el sistema de referencia y la evolucin temporal del vector suma.
La composicin vectorial muestra como un sistema trifsico de tensiones v
R
(t), v
S
(t) y
v
T
(t), que son los valores instantneos de las tensiones equilibradas por fase, se reduce
a un nico vector
V(t), mediante una suma vectorial en cada instante de tiempo, con-

sistente en la suma del valor en ese instante de v
R
(t), del valor de v
S
(t) girado 120
o
y
del valor de v
T
(t) girado 240
o
. Este vector espacial complejo resultante, gira con una
frecuencia al rededor del origen en el sentido de las agujas del reloj.
La expresin matemtica de este vector resultante es la siguiente:
V (t) =
2
3
_
V
R
(t) +
V
S
(t) e
j120
o
+
V
T
(t) e
j 240
o
_
(4.3)
Donde
2
3
representa el factor de escala aplicado al vector resultante para que su mdulo
coincida con el mdulo de los vectores de cada fase.
4.2.2. Sistema trifsico de tres hilos
De forma equivalente el desarrollo anterior puede hacerse igualmente para sistemas
trifsicos equilibrados de tres hilos de tensiones compuestas, donde cada tensin es la
diferencia entre dos tensiones simples, y que verican la ecuacin (4.1)
V
RS
(t) +
V
ST
(t) +
V
TR
(t) = 0
V
RS
(t) =
V
R
(t)
V
S
(t) v
RS
(t) =V
LL
sen(t +

6
)
V
ST
(t) =
V
S
(t)
V
T
(t) v
ST
(t) =V
LL
sen(t
2
)
V
TR
(t) =
V
T
(t)
V
R
(t) v
TR
(t) =V
LL
sen(t +
5
6
)
donde se verica la siguiente relacin:
V
LL
=
3V
m
En la gura 4.4 se muestra una representacin espacial de los vectores correspondientes
a las tensiones compuestas y a las tensiones simples de un sistema trifsico de tensiones
equilibradas. Se puede apreciar como todo el sistema de tensiones compuestas es el
mismo que el de tensiones simples pero desfasado 30
o
y con distinto mdulo.
Vt
Vs
Vrs
Vr
Vtr
Vst
30
Figura 4.4: Relacin entre las tensiones de fase y tensiones de lnea
El vector espacial resultante que representa al sistema trifsico de tensiones compuestas
para cualquier instante de tiempo se expresa en la ecuacin (4.4).
V (t) =
2
3
_
v
RS
(t) +
v
ST
(t) e
j120
o
+
v
TR
(t) e
j 240
o
_
(4.4)
Este vector puede descomponerse en dos componentes correspondientes a las partes
real e imaginaria:
V (t) =
(t) + j
(t)
Recordando la formula de Euler:
e
j
= cos+ j sen
y desarrollando la ecuacin (4.4), se obtiene que las componentes real e imaginaria del
vector resultante del sistema trifsico son:
4.3 Estados de conmutacin en inversores VSI 56
(t) =
2
3
_
v
RS
(t)
1
2
v
ST
(t)
1
2
v
TR
(t)
_
(4.5)
Como es sabido que:
v
RS
(t) =v
ST
(t) v
TR
(t)
se obtiene que:
(t) = v
RS
(t) =V
LL
sen
_
+

6
_
(4.6)
para la parte imaginaria se tiene que:
(t) =
2
3
_
3
2
v
ST
(t)
3
2
v
TR
(t)
_
=
3
3
V
LL
_
sen
_
t
2
_
sen
_
t +
5
6
__
(4.7)
(t) =V
LL
cos
_
t +
5
6
_
(4.8)
De donde se deduce que el vector espacial resultante viene expresado por la ecuacin
(4.9).

V (t) =
(t) + j
(t) =V
LL
e
j(t
3
)
(4.9)
Con esta expresin, se consigue una composicin vectorial en el plano complejo de un
sistema trifsico, obteniendo un nico vector giratorio representativo de dicho sistema.
4.3. Estados de conmutacin en inversores VSI
Un inversor est formado por tres semipuentes y cada uno de ellos dispone de dos in-
terruptores, unidireccionales en tensin y bidireccionales en corrientes; uno en la parte
superior y otro en la parte inferior. Las tres semiramas se conectan por los extremos a
la barra de continua (DC) y desde los puntos medios de los ramas, mediante inductores
en serie se conecta a un sistema trifsico de tensiones, tal y como se aprecia en la gura
4.5.
v
i
S
ap
S
an
S
bp
S
bn
S
cp
S
cn
V
a
V
b
V
c
V
DC
Figura 4.5: Inversor trifsico formado por tres semipuentes
Hay que tener en cuentas dos reglas:
1. Nunca se puede cortocircuitar la fuente de continua (condensador en el lado DC),
lo que se traduce en que en una misma no pueden dispararse los dos semicon-
ductores a la vez. As por ejemplo, para la fase a, si S
ap
=1 , entonces
S
an
=0 , que representan respectivamente las funciones de conmutacin en
una sola rama para los terminales positivo y negativo.
2. Nunca se deben dejar en circuito abierto las inductancias dispuestas en el lado de
alterna, lo que se traduce que siempre debe haber algn semiconductor condu-
ciendo en cada rama.
Atendiendo a las dos reglas anteriores, se tienen ocho combinaciones posibles del in-
versor trifsico. En la gura 4.6 se muestran los ocho estados de conmutacin, donde
los tres estados de los interruptores conectados al positivo del barraje de DC sirven pa-
ra identicar el estado de conmutacin y se representan con vectores. De esta manera
cuando un interruptor est cerrado se representa con un 1 y cuando est abierto
con un 0 .
Estas ocho conbinaciones son denominadas vectores de tensin de conmutacin (SVV
o Switching Voltaje Vectors) o estados de conmutacin.
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
Vr
Vs
Vt
Vdc
V1(100) V2(110)
V3(010)
V4(011)
V5(001)
V6(101)
V7(111) V8(000)
Figura 4.6: Estados de los interruptores
Los vectores V
o
y V
7
son los llamados vectores nulos. Son vectores excepcionales por-
que tienen una tensin nula pero que sin embargo cumplen con las dos reglas mencio-
nadas anteriormente [Vargas, 1990].
Si consideramos el vector de tensin de conmutacin V
1
que corresponde a la combi-
nacin (100), se deduce de la gura 4.6 que la tensiones entre los terminales a-b-c en
la salida del inversor VSI son:
V
ab
= +V
DC
V
bc
= 0
V
ca
=V
DC
Sustituyendo dichos valores en al ecuacin (4.4) se obtiene el valor del estado de con-
mutacin:
V
1
=
2
3
Vdc30
o
(4.10)
Este ngulo est medido desde la referencia V
ab
por lo que queda en fase con el eje .
Vector Estados V
ab
V
bc
V
ca
Expresin SVV
V0 (000) 0 0 0

V
o
= 0 V
V1 (100) +Vdc 0 -Vdc

V
1
=
2
3
Vdc0
o
V
V2 (110) 0 +Vdc -Vdc

V
2
=
2
3
Vdc60
o
V
V3 (010) -Vdc +Vdc 0

V
3
=
2
3
Vdc120
o
V
V4 (011) -Vdc 0 +Vdc

V
4
=
2
3
Vdc180
o
V
V5 (001) 0 -Vdc +Vdc

V
5
=
2
3
Vdc240
o
V
V6 (101) +Vdc -Vdc 0

V
6
=
2
3
Vdc300
o
V
V7 (111) 0 0 0

V
7
= 0 V
Tabla 4.1: Estados de conmutacin
Procediendo de la misma forma para los dems vectores, se obtienen las expresiones de
los vectores de tensin de conmutacin o SVV, donde la expresin de los ocho vectores
se da referida sobre el eje (ver tabla 4.1 ).
La ecuacin general que representa a los vectores de tensin de conmutacin es:
V
r
=
_
2
3
V
dc
e
j
(r1)
3
r = 1...6
0 r = 0, 7
(4.11)
A partir de la tabla 4.1 se pueden representar los distintos vectores SVV sobre los ejes
y , tal como se muestra en la gura 4.7.
V1
V2 V3
V4
V5
V6
(100)
(110) (010)
(011)
(011) (101)
(000) (111)
V7 V0
1
2
3
4
5
6
Figure 4.7: Posicin de los diferentes vectores
Como se puede apreciar, los extremos de los vectores forman un hexgono regular. El
rea encerrada entre dos vectores adyacentes dene un sector y en total se dispone de
seis sectores. Los vectores V
0
y V
7
corresponden con el valor cero y se han representado
en el centro del plano.
Mediante la tcnica de modulacin por ancho de pulso del vector espacial SVPWM,
se aplican diferentes vectores directores durante unos tiempos determinados, que de-
penden del vector de referencia que representa al sistema trifsico que desea obtener
en la salida, consiguiendo de esa forma una seal de salida con una menor distorsin
armnica. Este vector de referencia acta como consigna y gira dentro del plano ,
pasando por todos los sectores: cuantos ms vectores de conmutacin intermedios se
generen entre los vectores directores, la onda de salida se aproximar ms a una senoide
y por tanto ser menor la distorsin armnica total (THD)
4.4 Modulacin de vectores espaciales 61
4.4. Modulacin de vectores espaciales
La modulacin de los vectores espaciales consiste en construir el vector espacial

V
re f
aplicando alternativamente los dos vectores espaciales adyacentes que denen el sector
en el que se sita el vector
V
re f
.
Para obtener las ecuaciones generales se detallan como se obtienen los tiempos de
aplicacin de los vectores adyacentes cuando el vector

V
re f
que se encuentra en el
sector 1 (uno), situado entre los vectores directores adyacentes
V
1
y

V
2
, se descompone
como una suma del vector
d
1
mas
d
2
tal y como se observa en la gura 4.8.
d
1

V
1
V
re f
d
2
V
2
Figure 4.8: Descomposicin del vector intermedio
V
re f
situado en el sector 1
En el transcurso de un periodo de conmutacin denominado T
S
, se activarn los inter-
ruptores correspondientes al vector
V
1
durante un tiempo T
1
que corresponde al ciclo
de trabajo 1 o d
1
, y tambin se activarn los interruptores correspondientes al vector
V
2
durante un tiempo T
2
que corresponde al ciclo de trabajo 2 o d
2
. Para obtener estos
tiempos se dice que:
Z
T
S
0
V
re f
dt =
Z
T
1
0
V
1
dt +
Z
T
1
+T
2
T
1
V
2
dt +
Z
T
S
T
1
+T
2
V
0
dt (4.12)
donde:
4.4 Modulacin de vectores espaciales 62
T
s
V
re f
= T
1
V
1
+T
2
V
2
(4.13)
Ahora bien, aplicando la ecuacin general (4.11) para los vectores
V
1
y
V
2
en la ecuacin
(4.13) se tiene que:
T
s
V
re f
_
cos
sen
_
= T
1
2
3
V
DC
_
1
0
_
+T
2
2
3
V
DC
_
cos(/3)
sen(/3)
_
por tanto:
T
1
= T
s
V
re f
V
dc
sen
_
_
T
2
= T
s
V
re f
V
dc
sen()
El resto del tiempo del periodo de conmutacin estar ocupado por los vectores nulos,
este tiempo ser igual a:
T
0
= T
s
T
1
T
2
Una vez se conoce como calcular los ciclos de trabajo de los vectores directores que
representan los diferentes estados de conmutacin, es necesario analizar las opciones de
distribucin de los ciclos de trabajo a lo largo del periodo de conmutacin. Como se ha
visto, se tienen los dos ciclos de trabajo correspondientes a los dos vectores directores
que determinan el sector en el que se sita el vector espacial

V
re f
y el tiempo muerto
que puede ser asignado a los vectores nulos, por lo que las posibles combinaciones
para repartir estos ciclos de trabajo a lo largo del periodo de conmutacin son muy
numerosas [Rashid, 2004].
Captulo 5
Modelado del STATCOM
En este captulo se expone un modelo matemtico del STATCOM y el sistema AC.
Se describen las tcnicas de control del STATCOM, pasando primero por la tcnica
SVPWM (tcnica de modulacin por ancho de pulso del vector espacial), y despues
la tcnica de control de fase. sta ltima muy utilizada para el control de inversores
multipulso.
63
5.1 Introduccin 64
5.1. Introduccin
A continuacin se presentar el modelo matemtico del STATCOM, del cual, se deriva-
r el sistema de control, para su posterior implementacin de modelos ms detallados
en la herramienta de simulacin ATP (Alternating Transients Program). Se presentarn
dos tcnicas de control: la tcnica de control con SVPWM y la tcnica de control de
fase del STATCOM.
El STATCOM se emplea principalmente en sistemas de potencia para mantener las
tensiones en las barras mediante la compensacin dinmica de reactivos capacitivos o
inductivos. Tambin, se puede centrar su control en la regulacin del factor de potencia,
en amortiguar oscilaciones del sistema o una combinacin de estos, dependiendo de la
aplicacin y del nivel de tensin en el cual se desea implementar.
5.2. Representacinvectorial de cantidades trifsicas ins-
tantneas
Tal y como se hizo la transformacion del vector espacial en el captulo 4, se toman
estos resultados como base para representar las tensiones de un sistema trifsico me-
diante un nico vector en el marco de referencia . Esta representacin es compa-
rable con las muy estudiadas teoras generalizadas de potencia instantneas analizadas
por [E.H. Watanabe, 1993]. De la representacin del vector espacial, puede extenderse
hasta la introduccin en un sistema de coordenadas ortogonal en el cual cada vector se
describe por medio de sus componentes . La transformacin de variables de fase
(sean corrientes o tensiones) a coordenadas est dado por la siguente matriz de
transformacin:
[C] =
2
3
_
_
1
1
2

1
2
0
3
2

3
2
1
2
1
2
1
2
_
_
; [C]
1
=
3
2
[C]
t
_
_
i
0
_
_
= [C]
_
_
i
a
i
b
i
c
_
_
;
_
_
v
0
_
_
= [C]
_
_
v
a
v
b
v
c
_
_
(5.1)
5.2 Representacin vectorial de cantidades trifsicas instantneas 65
Ahora volviendo a referenciar las deniciones de potencia dadas por [E.H. Watanabe,
1993], en el cual dene una potencia activa instantanea p igual a :
p(t) = v
a
(t)i
a
(t) +v
b
(t)i
b
(t) +v
b
(t)i
b
(t) (5.2)
El diagrama de la gura 5.1 muestra dos vectores, uno es la representacin de la tensin
de la lnea de transmisin en el punto de conexin del STATCOM y el otro vector
representa la corriente existente entre el sistema AC y el STATCOM. Utilizando la
ecuacion (5.1) y remplazandola en (5.2) se lleva a la siguiente expresin de la potencia
instantnea
eje-Vc
V
i
i
V
V i
eje-Va
Figura 5.1: Marco de referencia
p =
3
2
(v
i
a
+v
) (5.3)
p =
3
2
[v[ [i[ cos
Donde es el angulo existente entre los vectores de tensin y corriente . Con esto se
muestra claramente que las componentes de corrientes que estan en fase con la tensin
son las que contribuyen a la potencia activa. De la misma forma se muestra la denicin
de potencia reactiva (demostracin mostrada en [E.H. Watanabe, 1993]) como:
q =
3
2
(v
) (5.4)
5.2 Representacin vectorial de cantidades trifsicas instantneas 66
q =
3
2
[v[ [i[ sen
Donde la constante de 3/2 es puesta para que coincida con la clasica denicin de
fasores bajo condiciones de sistemas de potencia en estado estable y balanceado .
La gura 5.2 muestra una manipulacin del sistema de coordenadas, utilizado ms ade-
lante para la separacin de variables en el sistema de control. Se dene un nuevo sistema
de coordenadas por un eje d, que est siempre coincidiendo con el vector de tensin
instantnea, y un eje q en cuadratura con el eje d. Los ejes d q no son estacionarios
en el plano, ellos siguen la trayectoria del vector de tensin y las coordenadas d q son
dadas por la siguiente transformacin variante en el tiempo denida en muchos libros
como la transformacin de Park, utilizada primeramente para el anlisis de maquinas
AC [Grainger y Stevenson, 1996] :
eje-Vc
Vd
i
eje-Va
q
d
= |V|
Figura 5.2: Marco de referencia d q
[T] =
2
3
_
cos() cos(
2
3
) cos(+
2
3
)
sen() sen(
2
3
) sen(+
2
3
)
1
2
1
2
1
2
_
_
; [T]
1
=
3
2
[T]
t
; =tan
1
(
v
)
_
_
i
a
i
b
i
c
_
_
= [T]
1
_
_
i
d
i
q
0
_
_
;
_
_
v
a
v
b
v
c
_
_
= [T]
1
_
_
[v[
0
0
_
_
(5.5)
5.3 Ecuaciones del circuito equivalente 67
y sustituyendo en la ecuacin (5.2) se puede obtener que:
p =
3
2
[v[ i
d
; q =
3
2
[v[ i
q
(5.6)
Bajo condiciones balanceadas y en estado estable, las coordenadas de los vectores de
tensin y corriente en el marco de referencia sncrono, son cantidades constantes. Esta
caracterstica ser muy utilizada para el posterior anlisis y el control desacoplado de
las componentes de tensin y corriente.
5.3. Ecuaciones del circuito equivalente
La gura 5.3, muestra una representacin simplicada de un STATCOM, incluyendo
un capacitor en el lado DC, un inversor y unas inductancias en serie en las tres fases
conectadas al sistema AC. Estas inductancias contabilizan la reactancia de dispersin de
los transformadores. El circuito tambien incluye una resistencia shunt con el capacitor
para representar las prdidas por conmutacion del inversor y la resistencia serie con la
lnea para representar las prdidas del transformador.
INVERSOR
R
R
R
L
L
L
C
Rp
Va
Vb
Vc
a e
b
e
c e
Ia
Ib
Ic
Figura 5.3: Esquema general del STATCOM
De acuerdo con la gura 5.3 , en trminos instantneos, las ecuaciones del circuito en
el lado AC se pueden explicar de la siguiente manera:
_
_
di
a
dt
di
b
dt
di
c
dt
_
_
=
_
R
L
0 0
0
R
L
0
0 0
R
L
_
_
_
_
i
a
i
b
i
c
_
_
+
1
L
_
_
v
a
e
a
v
b
e
b
v
c
e
c
_
_
(5.7)
Aplicandole la ecuacin (5.5) a la parte izquierda de la ecuacin (5.7) se tiene:
5.4 Sistemas de control del STATCOM 68
d
dt
_
_
i
a
i
b
i
c
_
_
=
d
_
T
1
_
dt
=
d
_
T
1
_
dt
_
_
i
d
i
q
0
_
_
+T
1
_
_
di
d
dt
di
q
dt
0
_
_
Ahora remplazando la ecuacin (5.5) en la ecuacin (5.7) y premultiplicando por [T]
se tiene que:
T
d
_
T
1
_
dt
_
_
i
d
i
q
0
_
_
+TT
1
_
_
di
d
dt
di
q
dt
0
_
_
=T
_
R
L
0 0
0
R
L
0
0 0
R
L
_
_
T
1
_
_
i
d
i
q
0
_
_
+
1
L
TT
1
_
_
[v[ e
d
e
q
0
_
_
(5.8)
Sabiendo que:
T
d
_
T
1
_
dt
=
_
_
0 0
0 0
0 0 0
_
_
; =
d
dt
se puede reducir la ecuacin (5.8) como sigue:
_
di
d
dt
di
q
dt
_
=
_

R
L

R
L
__
i
d
i
q
_
1
L
_
e
d
[v[
e
q
_
(5.9)
5.4. Sistemas de control del STATCOM
Hasta hace relativamente poco, los esfuerzos de investigaciones de los esquemas de
control del STATCOM se enfocaban principalmente en inversores o fuentes de tensin
reguladas tipo multi-pulso (control de fase) [C.Shauder y A.Edris, 1995]. Sin embargo
los recientes avances tecnolgicos, han permitido el aumento de la eciencia de los se-
miconductores de potencia ante frecuencias de conmutacin ms elevadas, permitiendo
a su vez el empleo de tcnicas SPWM [P. Garcia Gonzalez, 2000] y SVPWM como la
que se estudi en el captulo 4.
5.4.1. Tcnica de control por intermedio de la modulacin SVPWM
Para establecer este tipo de control se escogi el esquema desarrollado por [SCHAU-
DER y MEHTA, 1993], el cual est basado en la transformacin de Park. La ventaja
principal de esta teora tiene que ver con la posibilidad de desacoplar el control de P y
Q del STATCOM, adems de transformar ambas seales inherentemente sinusoidales
en variables DC que pueden ser manejadas por un sistema de control lineal.
Para explicar este sistema de control se emplear como referencia el sistema mostrado
en la gura 5.3. En este sistema R y L representan la resistencia y reactancia del trans-
formador respectivamente; e
a
, e
b
y e
c
las tensiones alternas a la salida del STATCOM.
Por ltimo las tensiones y corrientes del sistema AC corresponden a V
a
, V
b
, V
c
, i
a
, i
b
e
i
c
. La inspeccin de la ecuacin (5.9) conlleva a decir que se tendr un control desaco-
plado de i
d
e i
q
. De acuerdo con esto, la ecuacin de estado del sistema se puede denir
por la ecuacion (5.10):
_
di
d
dt
di
q
dt
_
=
_

R
L
0
0
R
L
__
i
d
i
q
_
1
L
_
x
1
x
2
_
(5.10)
Siendo x
1
y x
2
unas variables auxiliares de control utilizadas para eliminar la depen-
dencia entre e
d
e i
q
y entre e
q
e i
d
respectivamente. Se denen entonces los siguientes
controladores :
x
1
= k
r1
(1+
1
r1
)
_
i
dre f
i
d
_
(5.11)
x
2
= k
r2
(1+
1
r2
)
_
i
qre f
i
q
_
(5.12)
Finalmente las siguientes expresiones muestran la relacion entre la tensin en el capa-
citor y la tensin de salida del inversor
e
d
= kv
dc
cos (5.13)
e
q
= kv
dc
sen (5.14)
donde es el ngulo de desfasaje entre la componente fundamental de tensin de salida
del inversor y la tensin del sistema AC. Este control pretende mantener el ngulo
lo mas pequeo posible, manteniendo unos valores de magnitud aceptables para que
ocurra el intercambio de reactivos entre el sistema AC y el STATCOM.
Va
PARK
Vb Vc
ia
ib
ic
id iq
PLL
Va Vb Vc
CONTROLADOR
DE CORRIENTE
idref
iqref
Vtd
Vtq
PARK
Generador
SVPWM
INVERSOR
Seal
de disparo
Vd
id iq
Figura 5.4: Esquema general del sistema de control del STATCOM

5.4.1.1. Implementacin digital del modelo del STATCOM
La implementacin digital del modelo del STATCOM se hizo a travs de la herramienta
de simulacin ATP/EMTP y su interfaz grca ATPDRAW. La implementacin del
sistema de control se dividi en tres bloques: el de transformacin/medicin, el de
control y el de generacin de seales de disparo, tal y como se muestra en la gura 5.4
En primer lugar se sensan las tensiones del sistema V
a
, V
b
, V
c
, es decir, del punto de
conexin; y las corrientes i
a
, i
b
e i
c
, que uyen entre el sistema AC y el STATCOM,
para luego aplicarle la tansformacin de Park y as obtener sus correspondientes com-
ponentes en el sistema de referencia sncrono d q.
En el esquema de control de la gura 5.4 hay dos controladores separados que son in-
dependientemente capaces de controlar las componentes de eje directo y en cuadratura
de la corriente de salida del inversor, con un mnimo acople entre ellos. Estos controla-
dores generan los valores deseados de v
td
y v
tq
, que son las componentes de eje directo
y en cuadratura de la tensin de salida del conversor. A dichos valores se le hace una
transformacin de coordenadas d q a las coordenadas , la cual, ingresa como
vector de referencia en el bloque generador SVPWM, para que a su vez genere las
seales de disparo de los dispositivos de conmutacin del inversor.
PI
+
-
Vdc*
Vdc
idref PI
-1
ilq
iqref
PI
+
-
Vrms*
Vrms
iqref
a)
b)
c)
Figura 5.5: Controladores adicionales para las seales de referencia del sistema de con-
trol
Para generar las seales de referencia del controlador de corriente como i
dre f
e i
qre f
,
bien pueden ser valores de referencia introducidos por el operador en otra forma estas
seales de referencia son generadas por unos controladores adicionales dependiendo
del tipo de control que se tenga en el momento, tal y como lo muestran la gura 5.5.
El esquema de control mostrado por la gura 5.5a) compara los valores de ten-
siones del lado DC del compensador para luego introducir la seal de error en
un controlador proporcional-integral (PI) que dar origen a la seal de referencia
i
dre f
Si lo que se quiere, es corregir el factor de potencia en una carga conectada en la
barra de conexion del STATCOM, entonces se utiliza como seal de referencia
el valor negativo de la componente en cuadratura de la corriente consumida por
la carga i
lq
, tal y como se muestra en la gura 5.5b).
Adems se presenta un controlador adicional mostrado en la gura 5.5c), si lo que
se quiere es regular el valor de tensin en el punto de conexin del STATCOM,
y as, generar la seal de referencia i
qre f
.
El bloque mostrado en la gura 5.4 llamado controlador de corriente, se relata a conti-
nuacin:
Como se dijo anteriormente el controlador desarrollado se basa en el sistema de control
analizado por [SCHAUDER y MEHTA, 1993]. Las ecuaciones de estado que represen-
tan la dinmica de las corrientes del inversor estn dadas por la ecuacin (5.9). Las dos
componentes son desacopladas para obtener dos modelos independientes de planta de
primer orden con entradas de control x
1
y x
2
, como se muestra en la gura 5.6, donde:
1
sL + R
1
sL + R
x1
x2
id
iq
Figura 5.6: Modelos del controlador desacoplados
Vd
PI
+
-
id
wL*id
PI
+
-
iqref
iq
wL*iq
-
+
+
-
-
idref
x1
x2
Vtd
Vtq
Figura 5.7: Controlador de corrientes
x
1
=[v[ v
td
+Li
q
x
2
=v
tq
Li
d
Basado en los modelos desacoplados mostrados en la gura 5.6, un par de controladores
idnticos PI se pueden disear fcilmente . El esquema de control completo de corriente
es el que se muestra en la gura 5.7. Este tiene dos lazos de control separados que son
capaces de controlar las corrientes i
d
e i
q
independientemente. En cada lazo de control la
correspondiente seal de retroalimentacin se resta a la seal de corriente de referencia,
resultando una seal de error que alimenta a un controlador PI. Estos controladores no
tienen salidas limitadas, tienen ganancias idnticas y generan las seales de control x
1
y x
2
. En ltimas, tal como se muestra en la gura 5.4, se trata de una fuente inversora
de tensin controlada por corriente.
INVERSOR
Logica de
Disparo
Controlador
Calculo de
Potencia Reactiva
+
Control de
Sincronismo
Seal de
sincronismo
Seales de disparo
Seal de control
Medicion
de tensiones
Medicion
de corrientes
Referencia de
potencia reactiva
Figura 5.8: Esquema de control del STATCOM para el control de fase
5.4.2. Tcnica de control de fase
En muchos casos, el control STATCOM es funcin de la necesidades del sistema: con-
trol de potencia reactiva, estabilidad transitoria, regulacin de tensin, etc. Tambin las
tcnicas para estos nes son diversas, con sus ventajas y desventajas.
En esta seccin se presentarn algunos componentes del STATCOM de forma que fun-
cione como un controlador de potencia reactiva.
5.4.2.1. Componentes bsicos del sistema de control
Los componentes bsicos del control del STATCOM, mostrados en la gura 5.8, son
un controlador de sincronismo y un controlador de reactivos.
El sincronismo, entre la tensiones del sistema y las tensiones generadas por el conjunto
de inversores y transformadores se lleva a cabo por un oscilador de bloqueo de fase,
PLL, Phase Locked Loop.
Conforme se mmuestra, la potencia reactiva est relacionada con la amplitud de las
tensiones generadas por el STATCOM y estas tensiones estn relacionadas con la ten-
sin DC del capacitor. As mismo un control que actue sobre la tensin del capacitor es
capaz de controlar la potencia reactiva generada por el STATCOM.
Para completar, un bloque es necesario para la lgica de disparo de los dispositivos
de conmutacin. Este bloque recibe una seal de sincronismo y una seal de control
de reactivos y a partir de estas dos seales se construye la lgica de disparo de los
dispositivos de conmutacin de los inversores del STATCOM.
5.4.2.2. PLL, phase locked loop
Los osciladores con bloqueo de fase, PLL, son circuitos utilizados para sincronizar
seales.
Existen varios tipos de PLLs y la gura 5.9 muestra un ejemplo que fu utilizado por
[C.A.C. Cavaliere, 2000]. El PLL de la gura es un circuito que reliza la deteccin
de frecuencia y fase de las seales de secuecia positiva de las tensiones del sistema
elctrico.
En el esquema de la gura 5.9 las tensiones medidas se colocan en funcin de dos
valores fase-fase v
ab
y v
bc
. Los valores de tesin: v
a
, v
b
, v
c
, se normalizan.
Las tensiones fase-fase normalizadas se multiplican por corrientes cticias creadas a
travs de la realimentacin de la seal frecuencia . El resultado de esta operacin son
las potencias P1 y P2 , las cuales son sumadas resultando una nica potencia P.
En rgimen permanente las corrientes cticias tienen un ngulo de 90
o
adelantado con
respecto a las tensiones. Cuando esto ocurre, la potecia P es cero, el circuito se estabi-
liza, y las frecuencias y fases de las tensiones quedan detectadas.
Para una implementacin real, un PLL necesita de una seal de reinicializacin de su
integrador. Esto ocurre porque un integrador no puede seguir integrando indenida-
mente sin que haya saturacin. En el esquema mostrado en la gura 5.9 , el integrador
es inicializado cuando el valor de t + alcanza los + radianes, en este mismo
instante, este integrador recibe una condicin inicial de radianes.
Cond.
inicial
W
wt+teta
Ia
Ic
Va
Vb
Vc
Vab
Vcb
+
+
_
_
+
+
P1
P
P2
w
-pi
sen(wt+teta)
Controlador
PI
Controlador
PI
reset
wt +teta=pi
sen(wt+teta+2pi/3)
Figura 5.9: Esquema del PLL, phase locked loop.
5.4.2.3. Control de reactivos
En el caso de que las tensiones del sistema AC y a la salida del inversor del STATCOM,
enten en fase, no existir ujo de potencia activa entre el sistema AC y el STATCOM.
As mismo:
1. En el caso de que la amplitud de la red sea mayor que la del STATCOM, ste fun-
ciona como un banco de inductores equivalente con potencia reactiva inductiva.
2. En el caso de que la amplitud de la red sea menor que la del STATCOM, ste
funciona como un banco de capacitores equivalente con potencia reactiva capa-
citiva.
La amplitud de tensin del STATCOM depende de la amplitud de tensin del capacitor.
As mismo, controlando la tensin del capacitor, es posible controlar la amplitud de
tensin del capacitor.
En el caso en que el desfasaje entre las tensiones del sistema AC y el STATCOM sea
diferente de cero, habr un ujo de potencia activa entre el sistema AC y el STATCOM.
De este modo:
sen(wt +teta)
sen(wt +teta+2*pi/3)
Integrador
reset
wt+teta=pi
Cond.
inicial
Contrlador
PI
-pi
W
wt+teta
Ia
Ic
Va
Vb
Vc
Vab
Vcb
+
+
_
_
+
+
P1
P
P2
Desfasaje
Logica de
disparo
Seales de
disparo
Figura 5.10: Esquema del PLL+ seal de desfasaje.
1. Si las tensines del sistema AC estuvieran adelantadas, > 0, con relacin al
STATCOM, habr un ujo de potencia activa, entrando hacia el STATCOM, car-
gando el capacitor y consecuentemente aumentando su tensin.
2. Si las tensines del sistema AC estuvieran atrasadas, <0, con relacin al STAT-
COM, habr un ujo de potencia activa, saliendo desde el STATCOM, descargan-
do el capacitor y consecuentemente disminuyendo su tensin.
Como se muestra en la gura 5.10, la seal de sincronismo sumado a una seal de
desfasaje, es la entrada de referencia para la lgica de disparo de los dispositivos de
conmutacin del STATCOM.
En resumen, lograr que haya un pequeo desfasaje, , entre la seal de sincronismo
de las tensiones del sistema AC y la tensin del STATCOM, alterar la tensin del
capacitor y, consecuentemente, se alterar la tensin AC del propio STATCOM, de este
modo permitiendo el control de potencia reactiva.
5.4.2.4. Generacin de la seal de desfasaje
La seal de desfasaje, , sumada a la seal de sincronismo es la responsable del con-
trol de potencia reactiva del STATCOM. El objetivo del control determina la forma de
clculo de esta seal.
Utilzando las ecuaciones de potencia imaginaria instantanea denidas por [SCHAU-
DER y MEHTA, 1993], las cuales, estn detalladas por la ecuacin (5.4).
PI
+
-
q*
q
Figura 5.11: Controlador P-I de potencia reactiva
Un control proporcional-integral, como el mostrado en la gura 5.11, utilizando el
clculo de la potencia imaginaria instantnea q, ecuacin (5.4), comparada con una
referencia q. Entonces la salida del controlador proporcional-integral, es el ngulo de
desfasaje, necesario para el control de potencia reactiva.
Captulo 6
Modelo reducido del STATCOM
A continuacin se presenta un modelo reducido del STATCOM, despreciando las fun-
ciones de los elementos de conmutacin. ste modelo dinmico se incluye en un sof-
tware de simulacin de sistemas de potencia ante grandes perturbaciones (HAGP),
para ver las mejoras que se pueden obtener con la inclusin de un compensador de
esta naturaleza en un sistema de potencia elctrico.
78
6.1 Introduccin 79
6.1. Introduccin
El desarrollo y uso de los controladores FACTS, se ha dirigido en muchas aplicaciones
a mejorar la estabilidad de los sistemas de potencia elctricos.
Distintos modelos han sido propuestos para representar los controladores FACTS en
anlisis estticos y dinmicos. El presente captulo describe en detalle el modelo pro-
puesto para estos tipos de estudios, en los cuales se incluyen el compensador esttico
sncrono STATCOM.
stos modelos se basan en los estudios de estabilidad de tensin y de ngulo en los
reportes hechos por Claudio Caizares [Canizares., 2000] [Canizares., 1999].
Aunque no estn dentro del alcance de ste trabajo, hacer referencia al anlisis din-
mico de estabilidad de un sistema elctrico de potencia, es importante a la hora de
comprender cmo el sistema se modela por un conjunto de ecuaciones algebraico-
diferenciales y de sus mtodos de solucin.
6.2. Modelos del STATCOM
La estructura bsica de un STATCOM basado en un control PWM se muestra en la
gura 6.1.
El modelo del STATCOMpropuesto aqu, se basa en la ecuacin de balance de potencia
P = P
dc
+P
perdidas
(6.1)
la cual, representa basicamente el balance de potencia entre la potencia AC del contro-
lador P, y la potencia P
dc
, bajo la operacin balanceada y a frecuencia fundamental.
6.2.1. Modelo de estabilidad transitoria
Asumiendo un sistema de potencia balanceado, a frecuancia fundamental, el controla-
dor
1
se puede representar en estudios de estabilidad transitoria empleando el modelo
1
Cuando se hable de controlador en ste captulo tambin se hace referencia al STATCOM
6.2 Modelos del STATCOM 80
bsico mostrado en la gura 6.1. Las ecuaciones algebraico-diferenciales correspon-
dientes a este modelo pueden ser escritas como:
_
_
x
c

m
_
_
= f
c
_
x
c
, , m, V, V
dc
, V
re f
, V
dcre f
_
(6.2)
V
dc
=
VI
CV
dc
cos()
G
c
C
V
dc
R
C
I
2
V
dc
0 =
_
_
PVI cos()
QVIsen()
PV
2
G+kV
dc
VGcos() +kV
dc
VBsen()
Q+V
2
B+kV
dc
VGsen() kV
dc
VBcos()
_
_
. .
g(, m, V, V
dc
, V
re f
, V
dcre f
)
Magnitud
Controlador
Magnitud
+
-
+
-
+
-
+
-
V
I
P+ jQ a : 1
R+ jX
kV
dc
k(PWM)
V
dc
C
R
c
V
re f
V
dcre f
Figura 6.1: Modelo reducido del STATCOM con control de tensin PWM
Todas las variables se muestran en la gura 6.1. La admitancia G+ jB = (R+ jX)
1
,
se emplea para representar la impedancia del transformador y ltros AC series, G
c
=
1
R
c
se emplea para modelar las prdidas por conmutacin de los inversores y de sus
correspondientes circuitos snubbers, los cuales, tienen un efecto directo en la dinmica
de tensin del capacitor. La constante k =
_
3/8 m es directamente proporcional al
ndice de modulacin m.
+
-
+
-
+
+
+
-
K
v
(1+sT
1
)
(1+sT
2
)
V
re f
V
m
o
m
K
p2
+
K
I2
s
V
dcre f
V
dc

o
Figura 6.2: Tcnica de control PWM para control de tensin

Las variables x
c
y funciones f
c
() en la ecuacin (6.2) representan y relacionan las va-
riables del sistema de control del STATCOM, donde las variables x
c
, varan y dependen
de la tcnica PWM o de control de fase empleada para el sistema de control.
En la gura 6.2 se muestra el sistema de control PWM aplicado al STATCOM. Se
observa que en el controlador PWM, la magnitud de la tensin AC se controla a travs
del ndice de modulacin m, que tiene un efecto directo en la magnitud de la tensin del
conversor. El ngulo de fase , tiene la funcin de controlar el ujo de potencia activa
P dentro del controlador y por consiguiente el cargue y descargue del condensador, y
es empleado directamente para el control de la magnitud de tensin del condensador.
6.2.2. Modelo en estado estable
Las ecuaciones para el modelo en estado estable del STATCOM
2
se pueden obtener a
partir del modelo de estabilidad propuesto por las ecuacion 6.2, donde las ecuaciones
diferenciales propias del STATCOM y de su sistema de control se igualan a cero. stas
2
ver apndice
ecuaciones diferenciales en ltimas son las variables de estado que el modelo aporta al
sistema, es decir:
0 =
_
_
f
c
_
x
c
, , m, V, V
dc
, V
re f
, V
dcre f
_
PG
c
V
2
dc
RI
2
g(, m, V, V
dc
, V
re f
, V
dcre f
)
_
_
(6.3)
Para el sistema de control PWM m
o
=
_
8/3
V
re f
V
dcre f
. .
k
o
.
El modelo representado aqu, conlleva a la ms adecuada representacin para anlisis
en estado estable del STATCOM, en comparacin con los modelos tpicos de ujo de
cargas basados en una representacin del STATCOM como una fuente de potencia re-
activa. En estos tipos de modelos, los lmites usualmente se representan como lmites
de potencia reactiva, es decir, el STATCOM se representa tpicamente como un con-
densador sncrono, utilizando un modelo de barra estandar tipo PV.
A continuacin se presenta un caso que sirve como ejemplo de las representaciones
tpicas del STATCOM como barra tipo PV.
6.2.2.1. Sistema de 3 mquinas y 9 nodos
Como ejemplo de aplicacin se reproducen los resultados de la Western System Coor-
dinating Council (WSSC) [P.M. Anderson, 1994] que consiste de tres generadores y
nueve nodos. Fig 6.3. Los datos necesarios para la solucin del problema de ujo de
cargas para el sistema son: la tabla, 6.1,la cual, muestra las condiciones iniciales para
cada uno de los nodos del sistema. El nodo 1 (uno), se toma como referencia para to-
dos los clculos, y la tabla 6.3 indica los parametros de las lneas de transmisin. La
tensin de transmisin es 230kV y todos los datos son proporcionados en una base de
100 MVA.
Los valores nales de tensiones y ngulos en la solucin del ujo de cargas
3
se mues-
tran en la tabla 6.5 y los ujos de potencia por cada una de las lneas del sistema de
transmisin se encuentran en la tabla 6.7, En este estudio se obtienen prdidas de po-
tencia activa de 4,654 MW. Se debe considerar que stos resultados son sin la inclusin
del STATCOM.
3
Flujo de cargas solucionado por la herramienta de simulacin HAGP [Rodrguez, 2003].
5
1
2 3
4
6
7
8
9
STATCOM
Figura 6.3: Sistema 9 barras y 3 mquinas de la WSSC
Cuadro 6.1: Datos iniciales de los nodos
Bus Tensi
p.u
ngulo
de fa-
se
Pg
MW
Qg
MVAr
Pd
MW
Qd
MVAr
Tipo
bus
1 1,04 0,0 0,0 0,0 0,0 0,0 Ref
2 1,025 0,0 163 0,0 0,0 0,0 PV
3 1,025 0,0 85 0,0 0,0 0,0 PV
4 1,0 0,0 0,0 0,0 0,0 0,0 PQ
5 1,0 0,0 0.0 0.0 90 30 PQ
6 1,0 0,0 0,0 0,0 0,0 0,0 PQ
7 1,0 0,0 0,0 0,0 100 35 PQ
8 1,0 0,0 0,0 0,0 0,0 0,0 PQ
9 1,0 0,0 0,0 0,0 125 50 PQ
Cuadro 6.3: Datos de las lneas
Nodo envio Nodo recep R(p.u.) X(p.u.) B(p.u.)
1 4 0,0 0,0576 0,0
4 5 0,017 0,092 0,158
5 6 0,039 0,17 0,358
3 6 0,0 0,0586 0,0
6 7 0,0119 0,1008 0,209
7 8 0,0085 0,072 0,149
8 2 0,0 0,0625 0,0
8 9 0,032 0,161 0,306
9 4 0,01 0,085 0,176
Cuadro 6.5: Tensiones nodales sin la inclusin del STATCOM
Bus Tensi
p.u
ng
(Grad.)
Pg
MW
Qg
MVAr
Pd
MW
Qd
MVAr
Tipo
bus
1 1,04 0,0 71,65 26,76 0,0 0,0 Ref
2 1,025 9,274 163 6,74 0,0 0,0 PV
3 1,025 4,661 85 -10,74 0,0 0,0 PV
4 1,026 -2,213 0,0 0,0 0,0 0,0 PQ
5 1,013 -3,685 0,0 0,0 90 30 PQ
6 1,033 1,965 0,0 0,0 0,0 0,0 PQ
7 1,016 0,727 0,0 0,0 100 35 PQ
8 1,026 3,717 0,0 0,0 0,0 0,0 PQ
9 0,996 -3,986 0,0 0,0 125 50 PQ
Cuadro 6.7: Flujos de potencia por las lneas sin la inclusin del STATCOM
Nodo envo Nodo recep Flujo envo Flujo retorno
1 4 71, 65+ j26, 76 71, 64 j23, 64
4 5 30, 7+ j0, 89 30, 53 j16, 41
5 6 59, 47 j13, 59 60, 82 j17, 96
3 6 85 j10, 74 85+ j14, 84
6 7 24, 18+ j3, 12 24, 09 j24, 31
7 8 75, 91 j10, 69 76, 38 j0, 82
8 2 163+ j9, 08 163+ j6, 74
8 9 86, 62 j8, 26 84, 32 j11, 45
9 4 40, 68 j38, 55 40, 94+ j22, 45
6.3 Anlisis dinmico de estabilidad 85
Teniendo en consideracin los resultados obtenidos se lleva a cabo ahora un estudio de
ujo de cargas del mismo sistema incluyendo un STATCOM en el nodo 5. Para ste
anlisis el STATCOM se modela como un nodo cticio conectado a la barra 5 como
un nodo tipo PV. Para este problema se asigna a este nodo cticio el numero 10. Los
valores adicionales son la reactancia del transformador de acoplamiento, x
t
= 0,01 pu;
la magnitud de la tensin deseada Ve=1.025 p.u. (V
e
es la tensin a la salida del inversor,
es decir, en el nodo cticio) y la potencia activa que aporta el STATCOM al sistema es
P=0.
Los resultados que se obtienen al resolver el problema de ujos de cargas se resume en
las tablas 6.9 y 6.11. En ste caso las prdidas de potencia activa fueron de 4,63 MW.
De los resultados mostrados por las tablas 6.9 y 6.11 se puede concluir que el perl de
tensin se ve mejorado con la inclusin del STATCOM, en otras palabras el STATCOM
proporciona la potencia reactiva necesaria de forma tal que el nodo cticio tenga una
magnitud de tensin constante que hace posible los valores nodales del sistema. Es
importante destacar que al utilizar el STATCOM se disminuyen las prdidas de potencia
activa en la red y adems como se aprecia en las tablas 6.7 y 6.11, los ujos de potencia
reactiva cercanos al nodo donde se encuentra instalado el STATCOM sufren una ligera
reduccin. Se observa tambin que el ujo de potencia reactiva de los nodos 4-5 cambia
de direccin al emplear el dispositivo de compensacin. De la tabla 6.9 se nota que
el STATCOM se encuentra aportando 12,31 MVAR de potencia reactiva al sistema y
adems respetando la restriccin que se le impuso de que el ujo de potencia activa
entre el STATCOM y el sistema de potencia sea cero.
Se puede agregar que el modelo para ujos de cargas, colocando al compensador como
una fuente de potencia reactiva, no es capaz de incluir las prdidas atribuidas a los
dispositivos de conmutacin y a la resistencia de los transformadores. Adems no se
pueden representar bien los lmites del ndice de modulacin o los del ngulo de fase del
controlador que son de particular importancia en los estudios de estabilidad de tensin.
6.3. Anlisis dinmico de estabilidad
Como se dijo anteriormente, el hablar del fenmeno de estabilidad no concierne tanto
dentro de los objetivos propuestos por ste trabajo, pero sirve como base fundamental
Cuadro 6.9: Valores nodales con la inclusin del STATCOM en el nodo 5
Bus Tensin
p.u
ng
(Grad.)
Pg
MW
Qg
MVAr
Pd
MW
Qd
MVAr
Tipo
bus
1 1,04 0,0 71,63 19,61 0,0 0,0 Ref
2 1,025 9,275 163 4,65 0,0 0,0 PV
3 1,025 4,652 85 -15.18 0,0 0,0 PV
4 1,03 -2,205 0,0 0,0 0,0 0,0 PQ
5 1,024 -3,725 0,0 0,0 90 30 PQ
6 1,035 1.963 0,0 0,0 0,0 0,0 PQ
7 1,018 0,737 0,0 0,0 100 35 PQ
8 1,027 3,725 0,0 0,0 0,0 0,0 PQ
9 0,999 -3,961 0,0 0,0 125 50 PQ
10 1,025 -3,725 0,0 12,31 0,0 0,0 PV
Cuadro 6.11: Flujos de potencia por las lneas con la inclusin del STATCOM
Nodo envo Nodo recep Flujo envo Flujo retorno
1 4 71,63+ j19,61 71,62 j16,67
4 5 30,64 j7,03 30,48 j8,81
5 6 59,52 j8,89 60,87 j23,15
3 6 85 j15,18 85+ j19,34
6 7 24,13+ j3,81 24,04 j25,08
7 8 75,96 j9,92 76,43 j1,65
8 2 163+ j11,16 163+ j4,65
8 9 86,57 j9,51 84,28 j10,41
9 4 40,72 j39,59 40,98+ j23,7
10 5 0,0 j12,30 0,0+ j12,31
en el sentido de que al crear un modelo dinmico del STATCOM con las ecuaciones
algebraico-diferenciales que lo describen, se puede hablar de las mejoras que se pueden
obtener al insertar un dispositivo de compensacin como el STATCOM en un sistema
de potencia elctrico. Sabiendo adems que dentro de las principales estrategias del
control de estabilidad de los sistemas elctricos se encuentran los FACTS, en cuanto
amortiguamiento de oscilaciones, control de tensin, etc [Narain G. Hingorani, 2000].
Dentro de los fenmenos de estabilidad se encuentra el fenmeno de estabilidad de
gran perturbacin, el cual, estudia la respuesta de los elementos dinmicos del sistema
de potencia (generadores, cargas, elementos de compensacin dinmicos, enlaces de
corriente continua, etc) y la variacin de las diferentes magnitudes de la red (tensiones,
ujos por las lneas, frecuencia, etc) ante la ocurrencia de perturbaciones severas y
repentinas tales como cortocircuitos, prdidas de generacin o prdidas de carga.
Un sistema elctrico de potencia, bajo condiciones normales de operacin, se encuentra
casi en equilibrio y slo experimenta desviaciones menores causadas por variaciones
pequeas, casi continuas, en la demanda. Cuando una perturbacin severa y repentina
tiene lugar, hay cambios casi instantneos en las condiciones iniciales. La respuesta del
sistema a este tipo de eventos es altamente no lineal.
Un sistema elctrico de potencia es un sistema dinmico cuya respuesta ante la ocurren-
cia de una perturbacin grande se determina por modelos adecuados de los dispositivos
que lo componen. El sistema se modela por un conjunto de ecuaciones algbricas y
diferenciales que describen su comportamiento. Tales ecuaciones se deben resolver ite-
rativamente a partir de condiciones iniciales, mediante mtodos numricos para poder
determinar con exactitud la respuesta del sistema. Slo mediante este procedimiento se
puede hallar la evolucin en el tiempo de magnitudes del sistema tales como tensiones,
ujos, corrientes, etc.
El comportamiento de los elementos dinmicos se modela en el espacio de estados por
medio de ecuaciones diferenciales [Kundur, 1994] :
x = f (x, v) (6.4)
Las ecuaciones algbricas no lineales estn asociadas a la red:
i(x, v) = Yv (6.5)
donde
x = vector de estado del sistema
v = vector de componentes real e imaginaria de tensiones de barra del sistema
i = vector de componentes real e imaginaria de corrientes de barra del sistema
Y =matriz de la red del sistema, similar en estructura a la matriz admitancia de barras
Los dos conjuntos de ecuaciones forman un sistema algbrico-diferencial de primer
orden. Este debe ser resuelto a partir de unas condiciones iniciales (x
0
,v
0
) para obtener
la respuesta en el tiempo del sistema.
El anlisis de gran perturbacin es complejo desde el punto de vista computacional
debido a los grandes sistemas de ecuaciones que se deben resolver utilizando pasos de
tiempo pequeos para evitar errores y problemas de estabilidad numrica.
El modelo del STATCOM se incluy en la Herramienta Software para el Anlisis de
Estabilidad de Gran Perturbacin de Sistemas Elctricos de Potencia (HAGP), elabo-
rada en MATLAB y perteneciente a la Escuela de Ingenieras Elctrica, Electrnica y
Telecomunicaciones [Rodrguez, 2003].
Un diagrama de ujo simplicado de la herramienta software HAGP se muestra en la
gura 6.4, donde se muestran los pasos bsicos para la solucin del sistema elctrico
de potencia, incluyendo sus respectivas perturbaciones, dado el caso que sucedan.
Los detalles acerca de la forma de solucin y de la descripcin de los pasos para la
inclusin de modelos en dicha herramienta de simulacin, se describen en detalle en
dos tesis de pregrado de trabajos anteriores pertenecientes a la Escuela de Ingenieras
Elctrica, Electrnica y Telecomunicaciones [Rodrguez, 2003] [Meja, 2005].
DATOS DEL SISTEMA
FLUJO DE CARGAS
VALORES INICIALES
VARIABLES DE ESTADO
SIMULACION
DINAMICA
FIN
INICIO SIMULACION
t=0
PERTURBACION
MODIFICAR
EL SISTEMA
RESOLVER
ECUACIONES
DE RED
1
SI
NO
SI
NO
t = t +dt
1
PERTURBACION
MODIFICAR
EL SISTEMA
SI
VISUALIZAR
RESULTADOS
t >= Tiempo de simulacion
NO
FIN
SOLUCION DEL SISTEMA DE
ECUACIONES ALGEBRAICO-
DIFERENCIALES, VALORES
VARIABLES DE ESTADO Y
TENSIONES DE BARRA
RESOLVER
ECUACIONES DE RED
OBTENER NUEVOS
VALORES DE TENSION
SI
NO
Figura 6.4: Diagrama de ujo simplicado del HAGP
Captulo 7
Resultados de simulaciones
En este captulo se describirn y mostrarn los resultados de los diferentes casos de si-
mulacin del dispositivo FACTS (STATCOM), con sus respectivos modelos y sistema de
control. Se detallan algunos ejemplos para lograr mostrar algunas de las aplicaciones
ms frecuentes que tiene ste dispositivo cuando se inserta en un sistema de potencia.
90
7.1 Introduccin 91
7.1. Introduccin
En el presente captulo se mostrar el comportamiento de las simulaciones de un mo-
delo digital detallado del STATCOM mediante la herramienta de simulacin de transi-
torios electromagnticos ATP/EMTP.
Se exponen algunos casos de simulacin, implementandose su control digital, los cua-
les, se van a diferenciar dependiendo de el caso simulado y de las variables a controlar ,
como son: factor de potencia, regulacin de tensin en el punto de conexin del STAT-
COM y control de potencia reactiva.
7.2. Ejemplos implementados en ATP/EMTP y HAGP
A continuacin se presentan los datos y resultados de los ejemplos de sistemas elctri-
cos de potencia implementados en la herramienta de simulacin ATP/EMTP y HAGP
([Rodrguez, 2003]) .
7.2.1. Caso1
En este ejemplo se presenta un sistema de referencia para simulacin digital en el do-
minio del tiempo de un D-STATCOM
1
, controlado por la tcnica de modulacin SV-
PWM, a continuacin se describir el caso:
Una conguracin tpica de un STATCOM se muestra en la gura 7.1. El sistema con-
siste de una fuente con tensin Vs con su impedancia equivalente, al cual se le conecta
una carga tipo inductiva. El lado AC del STATCOM se conecta al sistema de potencia
a travs de ltros inductivos.
La variable a compensar para este ejemplo es el factor de potencia de la carga. Donde
se le aplica la tcnica de control por intermedio de la modulacin SVPWM. El sistema
de control se muestra en la gura 7.2.
En este caso de estudio el STATCOMse opera en un modo de operacin para correccin
del factor de potencia, donde lo que se busca es compenasar toda la potencia reactiva
consumida por la carga. A continuacin se describe el sistema AC:
1
Se llama al STATCOM de esta forma porque su aplicacin es sobre un sistema de potencia de
distribucin
7.2 Ejemplos implementados en ATP/EMTP y HAGP 92
Vdc
R
s
L
s
R
L
X
L
R
f
L
f
V
s
Figura 7.1: Sistema de potencia para el caso de simulacin 1
Parametro valor unidad
Frecuencia fundamental, f 60 Hz
V
s
112 V
L
s
145 H
R
s
25 m
L
f
0,637 mH
R
f
0,15
C 4820 F
V
DC
230 V
X
L
1,543
R
L
4,16
Cuadro 7.1: Parmetros del sistema de referencia
Va
PARK
Vb Vc
ia
ib
ic
id iq
PLL
Va Vb Vc
CONTROLADOR
DE CORRIENTE
idref
iqref
Vtd
Vtq
PARK
Generador
SVPWM
INVERSOR
Seal
de disparo
Vd
id iq
Figura 7.2: Sistema de control para el caso de simulacin 1

PI
+
-
Vdc*
Vdc
idref PI
+
-1
ilq
iqref
a)
b)
Figura 7.3: Seales de referencia del sistema de control para el caso de simulacin 1
La carga se model mediante su equivalente RL, la fuente de tensin se model como
fuente de tensin AC trifsica Type 14 en ATP. El STATCOM se conecta a travs de
un interruptor tripolar a los 0,3ms de haber arrancado la fuente de tensin Vs, la carga
trifsica inicialmente est fuera, al cabo de 3ms se conecta por intermedio de un inte-
rruptor tripolar, donde los tres polos se cierran simultneamente. Cabe anotar que las
seeales de disparo empiezan a actuar a los 3ms.
Dado que se desea compensar toda la potencia reactiva por la carga, las entradas de
referencia al bloque de control de corriente en la gura 7.2, para las corrientes de
referencia i
dre f
e i
qre f
viene dado por los controladores de la gura 7.3.
0.08 0.085 0.09 0.095 0.1 0.105 0.11 0.115 0.12 0.125 0.13
150
100
50
0
50
100
150
t[s]
T
e
n
s
i
o
n
[
V
]
Figura 7.4: Tensines del sistema AC y a la salida del conversor
De los resultados de la simulacon, en la gura 7.4 se muestran las tensiones del sistema
AC y la seal de tensin a la salida del conversor, en la cual se puede notar que entre
ellas no existe desfase, el cual, es requisito para el intercambio de potencia reactiva
entre el sistema AC y el STATCOM. Tambin se muestra el espectro armnico de la
tensin de salida del STATCOM en la gura 7.5, luego en la gura 7.6 se muestran las
componentes en cuadratura de las corrientes de la carga y la corriente de compensacin
que uye entre el STATCOM y el sistema AC.
Puede notarse que al rededor de los 80 ms el compensador ya muestra una respuesta
aceptable. Tambin se puede decir que en este caso el compensador est actuando en
modo capacitivo y est entregando reactivos capacitivos al sistema AC.
En la gura 7.7 se muestra la componente en cuadratura de la corriente suministrada
por la fuente de tensin V
s
, despues del transitorio se puede notar que esta corriente
Figura 7.5: Armnicos en la tensin de salida del compensador
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
10
5
0
5
10
t[s]
C
o
r
r
i
e
n
t
e
s
[
A
]
IqL
Iqcom
Figura 7.6: Componentes en cuadratura de las corrientes de la carga y del compensador
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
5
0
5
10
15
t(s)
(
A
)
Iqs
Figura 7.7: Componente en cuadratura de la corriente suministrada por la fuente de
tensin
tiende a cero, con pequeas oscilaciones debido a los armnicos generados. Con esto
se cumple el objetivo de compensar toda la potencia reactiva cosumida por la carga.
0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5
15
10
5
0
5
10
15
t[s]
C
o
r
r
i
e
n
t
e
s
[
s
]
Iqcom
IqL
Figura 7.8: Compensacin dinmica del factor de potencia de la carga
0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25
140
150
160
170
180
190
200
210
220
t(s)
(
V
)
Figura 7.9: Dinmica de la tensin del capacitor en el lado DC del compensador
Con este mismo ejemplo de simulacin se conmutaron dos cargas, una de tipo inductivo
y otra de tipo capacitivo, el evento simulado fu que primero la simulacin empezara
con compensacin del factor de potencia para la carga inductiva, luego al cabo de 1
(un) segundo se desconecta la carga inductiva y se conecta al mismo tiempo una carga
capacitiva tipo RC que entregue la misma cantidad de reactivos, equivalente al de la
carga inductiva. En la gura 7.8 se muestran las componentes en cuadratura de las
corrientes de compensacin y carga. En esta gura se puede notar la compensacin
dinmica del factor de potencia de la carga conectada al sistema.
En la gura 7.9 se puede notar como el capacitor se carga y se descarga al pasar de
una compensacin del factor de potencia inductivo a una compensacin del factor de
potencia capacitivo, permitiendo el cambio del ujo de potencia reactiva.
7.2.2. Caso 2
El ejemplo que se relata a continuacin es un ejemplo muy empleado en [Canizares.,
2000].
El sistema de potencia utilizado para la simulacin se muestra en la gura 7.11. El
objetivo de la inclusin del STATCOM en la barra ocho del sistema de potencia es la
regulacin de tensin, para mejorar el perl de tensin en estado estable y en estado
dinmico. En las siguientes tablas se presentan los parmetros del sistema de potencia,
de la mquina sncrona y del STATCOM.
El STATCOM utilizado es una conguracin de 12 pulsos mostrado por la gura 7.10,
donde el transformador de acoplamiento posee una reactancia de dispersin X =7, 5 %
con una potencia de 150 MVA, los transformadores que conforman el circuito mag-
ntico de eliminacin de armnicos, poseen una rectncia de dispersin cada uno de
3,75 % en una potencia de 75 MVA/cu.
Inversor de 6
pulsos
Inversor de 6
pulsos
+30
o
Figura 7.10: Conguracin de un STATCOM de 12 pulsos
El Regulador Automtico de Tensin(AVR) del generador utilizado es el que se muestra
en la gura 7.12.
Variable Valor [p.u] Variable Valor [p.u]
Polo 2 H 2,7113
Ra 0,001096 D 0
X
l
0,15 T
/
d
0
6,19488
X
d
1,7 T
/
q
0
0
X
q
1,64 T
//
d
0
0,028716
X
/
d
0,238324 T
//
q
0
0,07496
X
//
d
0,18469 X
//
q
0,185151
Cuadro 7.2: Parmetros de la maquina sncrona, con bases de 200MVA y 13.8kV
Elemento R X B/2
1-2 0,0003 0,0684 0
2-3 0,0159 0,2275 0,0754
3-8 0,022 0,316 0,1047
4-3 0 0,08 0
5-6 0,0026 0,0379 0,0126
6-7 0 0,12 0
8-9 0,0026 0,0379 0,0126
9-10 0 0,12 0
Filtro 0,0087 -4,3 0
Cuadro 7.3: Parmetros del sistema de potencia en p.u., con potencia base de 240 MVA
Variable Valor p.u.
R 0
X 0,0145
Gc 2,16
C 0,0432
Cuadro 7.4: Parmetros del STATCOM en p.u., con potencia base de 150 MVA y 12
kV
B-1 B-2
B-3
B-4 B-5
B-6
B-8 B-9
B-10
Z Thevenin
Barra infinita
230 kV
13,8 kV
Generador
245,5 MW
952 MW
267 MVAr
321 MW
43 MVAr
STATCOM
Figura 7.11: Sistema de potencia para el caso de simulacin 2
750
(0,11+0,01s)(0,17+0,95s)
0,04s
1+s
1
1+0,001s
1
+
-
+
+
1
+
V
4
V
f
Figura 7.12: Regulador Automtico de Tensin, AVR
3,962(1.5+0,0115s)
(1+0,001s)(1+0,03s)
1
-
+
+
-
V
8
Figura 7.13: Sistema de control de fase

7.2.2.1. Modelo del STATCOM en el software HAGP
El modelo del sistema de control empleado se explico en el captulo anterior. Los pa-
rmetros elctricos del STATCOM y del sistema de potencia se relatan en la tablas 7.2,
7.3 y 7.4.
Tiempo Evento
0,00 s Inicio de simulacin
4,5 s Conexin de una gran carga en barra 6
4,65 s Desconexin de carga de la barra 6
10 s Fin de simulacin
Cuadro 7.5: Secuencia de eventos de la simulacin caso 2
La descripcin de los eventos para el caso simulado bajo la herramienta de simulacin
HAGP se detallan en la tabla 7.5.
0 1 2 3 4 5 6 7 8 9 10
0.9965
0.997
0.9975
TENSIN EN BARRA BUS
8

Tiempo(s)
M
a
g
n
i
t
u
d

d
e

t
e
n
s
i
n
(
p
u
)
Figura 7.14: Magnitud de tensin en la barra 8
Como puede observarse, las respuestas dinmicas dadas por el software HAGP, en
cuento a el ngulo del rotor del generador(gura 7.15) y la tensin en el barra 8 del
sistema de potencia(gura 7.14) ante la ocurrencia de la perturbacin descrita anterior-
mente dieren un poco de las respuestas dinmicas dadas con los modelos creados en
ATP.
Tambin se debe tener en cuenta que los sistemas de excitacin del generador son
diferentes en sta simulacin y en la siguiente, stos se cambiaron para mejorar la
dinmica de la tensin a la salida del generador.
0 1 2 3 4 5 6 7 8 9 10
55
55.5
56
56.5
57
NGULO DEL ROTOR DEL GENERATOR
Tiempo(s)
n
g
u
l
o

d
e
l

r
o
t
o
r
(
g
r
a
d
o
s
)
Figura 7.15: ngulo del rotor del generador
Es importante resaltar que estos tipos de modelos desprecian los armnicos creados por
las conmutaciones de los interruptores de potencia.
7.2.2.2. Tcnica de control de fase
El sistema de control implementado para el STATCOM, es por medio de la tcnica
de control de fase explicada en el captulo anterior donde el controlador utilizado se
muestra en la gura 7.13.
En la gura 7.11 se conecta una carga de valor elevado a la barra 6 a los 4,5s y despus
el interruptor colocado entre barras 3-5 se abre en un tiempo de 4,65s para desconectar
la carga (952 MW 267MVAR, con su equivalente RL en la herramienta ATP/EMTP).
Ahora se observar la perturbacin del sistema y de como actua el STATCOM como
regulador de tensin.
Con la inclusin del STATCOM en el sistema, el perl de tensin de la barra 8, se ve
mejorado tal y como lo muestra la gura 7.16, el cual pasa de 0,82 p.u. a 1 en p.u.
0 0.2 0.4 0.6 0.8 1 1.2
0
0.2
0.4
0.6
0.8
1
t(s)
(
p
u
)
Figura 7.16: Inclusin del STATCOM en barra 8 a los 0,5 seg
1.15 1.155 1.16 1.165 1.17 1.175 1.18 1.185 1.19 1.195 1.2
1
0.5
0
0.5
1
x 10
4
t(s)
(
V
)
Figura 7.17: Forma de onda de la tensin a la salida del inversor
Figura 7.18: Espectro armnico de tensin en la salida del inversor
4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8 6
0.2
0.4
0.6
0.8
0.9
1
t(s)
(
r
a
d
)
Figura 7.19: ngulo del Generador
4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8 6
0.1
0.05
0
0.05
0.1
0.15
t(s)
(
r
a
d
)
Figura 7.20: ngulo de desfase del control
La forma de onda de tensin de la barra ocho muestra un risado debido a los armnicos
incluidos por las conmutaciones de los interruptores del circuito inversor (gura 7.17).
Su espectro armnico se muestra en la gura 7.18, donde solo aparecen los armni-
cos 11 y 13, por la actuacin del ltro RC implementado en la salida del inversor del
STATCOM.
La gura 7.19 muestra el transitorio del ngulo del generador, dado en radianes. sta
se estabiliza al rededor del segundo despues de la perturbacin, apreciandose la contri-
bucin a la estabilidad de ngulo del sistema de potencia.
Tambin la gura 7.20, muestra el transitorio del ngulo de desfase existente entre la
tensin del sistema AC y la tensin en lado de baja del transformador de acople. Cuando
la perturbacin se aplica al sistema de potencia, la tensin del conversor sufre una caida
de tensin, por consiguiente el sistema de control, acta al principio del transitorio para
aumentar el angulo de desfase, y con sto proporcionar los reactivos sucientes para
mantener la tensin en el punto de conexin del STATCOM. Por tanto, sin ste desfase
el cargue y descargue del capacitor seria imposible de controlar, que en ltimas es la
tensin en el capacitor, la que se controla para producir el intercambio de potencia
reactiva y activa entre el sistema de potencia y el STATCOM.
7.2.2.3. Esfuerzos en los elementos de conmutacin
En la graca mostrada por la gura 7.21 se muestra la corriente que aparece en uno de
los elementos de conmutacin con un valor de corriente al rededor de 15kA.
Figura 7.21: Corriente en un dispositivo de conmutacin
Como el STATCOM, mejor el perl de tensin de la barra 8, entonces quiere decir
que las exigencias de corriente y de tensin son altas, para una conguracin de 12
pulsos, por tanto, es recomendable ampliar el nmero de pulsos para que la potencia se
distribuya en todos los conversores de 6 pulsos y as, atenuar la demanda de corriente.
7.2.3. Caso 3
El ejemplo correspondiente al caso 3 se implement con la herramienta de simulacin
ATP/EMTP. Se trata del modelado y simulacin de un STATCOM con una congu-
racin de 48 pulsos, con la tcnica de control de fase, como se explic en el captulo
anterior. Este ejemplo es referenciado en [C.A.C. Cavaliere, 2000].
El sistema de potencia se simula un STATCOM conectado directamente a una barra in-
nita, centrando su modelado en el STATCOM, del cual, sus componentes se modelan
por intermedio de libreras, las cuales, dan origen a la creacin de modelos ms comple-
jos, en base de los propios modelos de ATP. stas libreras son hechas bajo la subrutina
DATA BASE MODULE que posee la herramienta de simulacin ATP. La creacin de
libreras para representar los componentes del STATCOM, es muy til porque ayuda
a minimizar la cantidad de modelos de ATP montados sobre su interfaz graca (ATP-
Draw).
La barra innita est representada como una fuente de tensin con su impedancia equi-
valente de Thevenin. Su potencia de cortocircuito S
cc
= 2000MVA, con una tensin
nominal de 138 kV . La impedancia de Thevenin se calcula como:
X
th
=
V
2
nominal
S
cc
7.2.3.1. Inversor
La capacidad de interrupcin de corriente y soporte de tensin de los dispositivos de
conmutacin determinan la potencia de los inversores. Los dispositivos semiconduc-
tores se modelan por medio interruptores ideales en ATP Type 13 (interrptores con-
trolados por TACS), donde a cada dispositivo de conmutacin se le adiciona una red
snubber tipo RC. Para el modelo digital implementado se consider que los inversores
tendran una potencia de 12,5 MVA. Estos inversores se basan en dispositivos semicon-
ductores con caractersticas de bloqueo de 6kA/6kV (para una potencia de 100 MVA
son necesarios 8 inversores).
7.2.3.2. Tranformadores
Los transformadores se modelaron a travs del acople de modelos de transformadores
monofsicos saturables que posee ATP. Que necesitan como parmetros de entrada, las
resistencias y las inductancias de dispersin de cada devanado de los transformadores.
Se despreciaron las caractersticas de saturacin de los transformadores.
Para la conguracin estrella-delta de los transformadores para la estructura de re-
duccin de armnicos, donde se necesitan 4 transformadores estrella-estrella y cuatro
transformadores estrella-delta, cumpliendo respectivamente con las relaciones de trans-
formacin de 1:1 para los transformadores estrella-estrella y 1 :
3, para tranformado-
res estrella-delta.
Se consider una salida de los inversores de 4,3kV, donde la disposicin en serie de 8
transformadores dan una suma de 34,3 kV. Deniendose una base en tensin de estos
transformadores de 4,3 kV con una potencia base de 12,5 MVA , una reactancia del 8 %
y una resistencia de 0.1 %.
El transformador de acoplamiento tiene una relacin de transformacin 138kV Y
34,5kV (4 :
3) respectivamente. Con una reactancia de 5 % y una resistencia de

0.1 % , en una base de 138kV y una potencia de 100MVA.
El capacitor que se utilizo tiene un valor de 4820 F.
En la tabla 7.6 se dene la secuencia de los eventos de la simulacin.
Tiempo [s] Evento
0,00 Inicio de simulacin
0,05 Partida del STATCOM(sin el PLL activado)
0,25 Partida de los inversores y control de potecia reactiva
0,5 Modo inductivo- referencia de potencia reactiva +1 p.u.
0,75 Modo capacitivo- referencia de potencia reactiva -1 p.u.
1,00 Referencia nula-referencia de potencia reactiva cero
1,25 n de simulacin
Cuadro 7.6: Secuencia de eventos de la simulacin
A continuacin se muestran las grcas que denen la simulacin.
En los primeros instantes de simulacin se muestra la dinmica de partida de la seal
de frecuencia (gura 7.22), , obtenida del PLL, se muestra un transitorio de aproxi-
madamente de 150 ms para estabilizarse en una frecuencia de 376 rad/seg, tambin se
observa de este resultado una pequea variacin de la seal creada apartir de la conmu-
tacin de los interruptores a los 0,25 s de simulacin.
0 0.05 0.1 0.15 0.2 0.25 0.3
200
100
0
100
200
300
377
500
t(s)
r
a
d
/
s
Figura 7.22: Frecuencia en radianes por segundo
0 0.05 0.1 0.15 0.2 0.25 0.3
4
pi
2
0
2
pi
4
t(s)
r
a
d
/
s
Figura 7.23: Seal de sincronismo del PLL
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
0.4
0.2
0
0.2
0.4
0.6
0.8
1
t(s)
(
p
u
)
Figura 7.24: Potencia imaginaria instantnea inicial
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
2
1.5
1
0.5
0
0.5
1
1.5
2
t(s)
(
p
u
)
Figura 7.25: Potencia activa instantnea inicial
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
0
0.5
1
1.5
2
t(s)
(
p
u
)
Figura 7.26: Tensin del capacitor
Adems, la gura 7.23 muestra la seal de sincronismo, t, generada por el PLL. La
interferencia del accionamiento de los inversores no se nota mucho porque el PLL tiene
un integrador que acta como ltro.
La gura 7.24 muestra la seal de potencia imaginaria instantnea, q, durante el inter-
valo de partida del STATCOM.
La gura 7.25 muestra la seal de potencia activa instantnea, p, durante el intervalo
de partida del STATCOM.
La gura 7.26 muestra el transitorio que ocurre con la tensin del capacitor, donde se
producen unas oscilaciones para luego en aproximadamente 90ms se estabilice a su
valor nominal.
Los picos observados en las potencias imaginarias y activas instantneas son conse-
cuentes del cargado inicial del capacitor. En el instante de partida del STATCOM, el
capacitor est descargado y con tensin cero. Cuando los inversores se accionan, per-
miten el paso de energa del sistema al capacitor, para cargarlo a su tensin nominal,
que corresponde a un valor nulo de potencia imaginaria instantnea.
0.22 0.24 0.26 0.28 0.3 0.32 0.34
2
1.5
1
0.5
0
0.5
1
1.5
t(s)
(
p
u
)
Icomp
V48p
Vs
Figura 7.27: Tensin del sistema V
s
, tensin 48 pulsos y corriente de compensacin Ic
En la gura 7.27 se muestra la tensin del sistema AC, Vs, la tensin del conversor de
48 pulsos, y la corriente de compensacin de la fase a, normalizadas. El transitorio
de partida del STATCOM, tambin interere con los valores de tensin y corrientes.
La tensin de 48 pulsos, est relacionada con la tensin del capacitor, por este motivo
los transitorios de la tensin del capacitor, afectan a la tensin de 48 pulsos, pero a su
vez la tensin de 48 pulsos interactua con la tensin del sistema y afecta la corriente de
compensacin.
7.2.3.3. Resultados de la respuesta dinmica del STATCOM
La gura 7.28 muestra la potencia imaginaria del STATCOM y la potencia imaginaria
instantnea de referencia, en el intervalo de actuacin del control de potencia reactiva.
Con este resultado fu probada la capacidad del STATCOM para operar tanto en modo
inductivo, como en el modo capacitivo con la potencia deseada. La potencia imaginaria
sigue perfectamente la referencia del control.
0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2
1
0.5
0
0.5
1
t(s)
(
p
u
)
Figura 7.28: Potencia imaginaria instantnea
0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2
1
0.5
0
0.5
1
t(s)
(
p
u
)
Figura 7.29: Potencia activa instantnea
0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2
0
0.2
0.4
0.6
0.8
1
1.2
t(s)
(
p
u
)
Figura 7.30: Tensin del capacitor
La gura 7.29 muestra la potencia activa instantnea p, en este resultado se observan
las variaciones de potencia activa necesarias para el funcionamiento del STATCOM.
0.7 0.72 0.74 0.76 0.78 0.8 0.82 0.84 0.86 0.88
1.5
1
0.5
0
0.5
1
1.5
t(s)
(
p
u
)
Icomp V48p
Figura 7.31: Tensin del sistema V
s
, tensin 48 pulsos y corriente de compensacin Ic
En la gura 7.30 se muestra la tensin en el capacitor, respecto a los modos de com-
pensacin, cero, inductivo y capacitivo.
En al gura 7.31 se muestra la tensin del sistema AC, Vs, la tensin de 48 pulsos, V
48
,
y la corriente de compensacin en el instante que pasa del modo inductivo al modo
capacitivo.
Observar en la gura 7.31 que en el modo inductivo, la tensin del sistema AC, est
por encima de la tensin de 48 pulsos y la corriente de compensacin est atrasada 90
o
de las tensiones. Cuando pasa al modo capacitivo, la tensin de 48 pulsos, sobrepasa a
la tensin del sistema AC y la corriente de compensacin pasa a estar 90
o
adelantadas
de las tensiones.
7.2.3.4. Esfuerzos en los elementos de conmutacin
En las grcas de las guras 7.32, 7.33 y 7.34 se muestran las corrientes que cir-
culan por unos de los dispositivos de conmutacin modelados en este ejemplo con
ATP/EMTP, cuando el STATCOM se encuentra operando en modo capacitivo, induc-
tivo y nulo, respectivamente. La mxima corriente demandada para tal propsito es
cercana 2500 A .
Figura 7.32: Corriente que circula por uno de los elementos de conmutacin, en modo
inductivo
Figura 7.33: Corriente que circula por uno de los elementos de conmutacin, en modo
capacitivo
Figura 7.34: Corriente que circula por uno de los elementos de conmutacin, con refe-
rencia de potencia reactiva cero
Figura 7.35: Tensiones de los elementos de conmutacin, en modo inductivo
En las grcas de las guras 7.35, 7.36 y 7.37 se muestran las tensiones que aparecen
en los extremos de los dispositivos de conmutacin, cuando el STATCOM se encuentra
operando en modo capacitivo, inductivo y nulo, respectivamente.
Figura 7.36: Tensiones de los elementos de conmutacin, en modo capacitivo
Figura 7.37: Tensiones de los elementos de conmutacin, con referencia de potencia
reactiva cero
Cuando el STATCOM est operando en modo capacitivo, aparece la mxima tensin
entre los extremos de los dispositivos de conmutacin. sta tensin corresponde a unos
6600 V.
Es importante resaltar que en cuanto a las especicaciones dadas en la seccin 7.2.3.1
para los dispositivos de conmutacin (6kV/6kA), los valores de tensin de bloqueo
estan muy ajustados, an sabiendo que no se ha simulado el ejemplo para una pertur-
bacin del sistema considerable. Solo los lmites de corrientes no se superaron para los
propsitos de esta simulacin.
Captulo 8
Conclusiones y sugerencias
La apropiacin de este tipo de tecnologas por medio de modelos digitales es muy
importante, en el sentido de que pueden predecir el comportamiento del sistema y su
interaccin con otro tipos de dispositivos.
El STATCOM tiene la gran capacidad de producir potencia reactiva bien sea inductiva
o capacitiva, a partir de la misma energa que le entrega la red de energa elctrica.
Las tcnicas de control de ste tipo de dispositivo, dependen de la aplicacin que se
quiera en el momento, ya que las tcnicas de modulacion PWM ofrecen mayor ca-
pacidad de respuesta en comparacin con otras, como las tcnicas de control de fase,
sta ltima, muy utilizada en los ltimos aos para conguraciones de STATCOM-
multipulsos.
Se comprob la versatilidad que puede tener el STATCOM, en el sentido de la estrategia
y de las variables a compensar, factor de potencia, regulacin de tensin en un punto
de conexin, control de potencia reactiva.
A nivel de Distribucin es mas aconcejable montar equipamientos controlados por tc-
nicas PWM debido a su rapidez de respuesta, teniendo en cuenta que los armnicos
generados pueden ser atenuados por intermedio de ltros adecuados.
La creacin de los modelos digitales fu muy decisiva para el buen desempeo de las
simulaciones. El modelo del STATCOM y del sistema de potencia se llev a cabo por
la creacin de librerias, a travs de la subrutina DATA BASE MODULE que contiene
ATP para la creacin de modelos mas acordes con los requerimientos de la simulacin.
El porqu de tantos modulos creados para la simulacin del STATCOM de 48 pulsos
es que antes de simular un equipamento tan detallado, primero hay que garantizar la
115
8.1 Aportes originales del trabajo 116
estabilidad numrica. Las oscilaciones que se presentan se deben a la cantidad de ele-
mentos no lineales en la simulacin. Dichas oscilaciones son propias y se acentan an
ms en los programas de simulacon que utilizan en su algoritmo de solucin, la regla
trapezoidal.
En cuanto a las simulacines en ATP, es muy impotante encontrar un buen punto de
operacin de los sistemas de potencia (especialmente para sistemas de potencia media-
nos y grandes), para poder hacer comparacines mas precisas con otros programas de
simulacin. Estos se pueden mejorar utilizando rutinas propias de ujos de cargas de
ATP como lo son FIX SOURCE y CAO LOAD FLOW, para obtener unas condiciones
iniciales en estado estable ms acordes entre un software y otro.
8.1. Aportes originales del trabajo
Creacin de modelos detallados de los componentes bsicos del STATCOM en
ATPDraw y del sistema de potencia, ya que ste no cuenta con modelos base para
llevar a cabo este tipo de simulaciones.
Estructuracin de uno de los modelos mas representativos de los dispositivos
FACTS, en el caso especico de este trabajo: EL STATCOM.
Creacin de un modelo reducido del STATCOM para su posterior inclusin en
el software de simulacin de sistemas de potencia para el anlisis de estabilidad
ante grandes perturbaciones (HAGP).
8.2. Sugerencias para posteriores desarrollos
A nivel acadmico se sugiere la creacin de un banco de modulos de libreras en
ATP, para hacer las simulaciones posteriores mas rpidas y fciles, por ejemplo:
sistemas de protecciones, modulos estandar para el sistema de excitacin de la
maquina sncrona, entre otros.
Crear una estrategia para la linealizacin del sistema de potencia con sus dis-
positivos de control para encontrar de una manera mas fcil los parmetros de
los diferentes sistemas de control, que al nal darn una respuesta dinmica ms
acorde con la teora.
8.2 Sugerencias para posteriores desarrollos 117
Vericar nuevas topologas para el STATCOM, utilizando inversores trifsicos
de tres niveles, transformadores en zig-zag.
Investigar el desempeo dinmico del STATCOM utilizando adems del conden-
sador, elementos de almacenamiento de energa como:
Materiales superconductores
Bateras, etc.
Disear estrategias de control para sistemas de potencias desbalanceados
Bibliografa
C.A.C. CAVALIERE E.H. WATANABE M. A. ANLISE DE OPERAO DE
STATCOM EM SISTEMAS DESBALANCEADOS. CBA2000, FLORIANOPO-
LIS,BRASIL 2000.
CANIZARES. C. A. MODELING AND IMPLEMENTATION OF TCR AND VSI BASED
FACTS CONTROLLERS. Internal Report, ENEL and Politecnico di Milano,
Milan, Italy, OCT. 1999.
CANIZARES. C. A. POWER FLOW AND TRANSIENT STABILITY MODELS OF
FACTS CONTROLLERS FOR VOLTAGE AND ANGLE STABILITY STUDIES.
Proc. of IEEE/PES Winter Meeting,Singapure, JANUARY 2000.
C.SCHAUDER E.STACEY M. L. L. A. A. AEP UPFC PROJECT: INSTALLATION,
COMMISSIONING AND OPERATION OF THE +/- 160 MVA STACOM (PHASE
I). IEEE Trans. on Power Delivery, VOL. 13, NO. 4, PP. 15301535, OCTO-
BER 1998.
C.SHAUDER M.GERNHARDT E. T. L. T. Y A.EDRIS. DEVELOPMENT OF A +/-
100 MVAR STATIC CONDENSER FOR VOLTAGE CONTROL OF TRANSMISSION
SYSTEMS. IEEE Trans. on Power Delivery, VOL. 10, NO. 3, PP. 371 379,
JULY 1995.
EDWARDS C. ADVANCED STATIC VAR GENERATOR EMPLOYING GTO THYRIS-
TORS. IEEE PES Winter Power Meeting, VOL. PWRS-4, NO. 38WM109-1,
PP. 177179, SEPTEMBER 1988.
EDWARDS C. ADVANCED STATIC VAR GENERATOR EMPLOYING GTO THYRIS-
TORS. IEEE Transactions on Power Delivery, VOL. 3, NO. 4, PP. 16221627,
OCTOBER 1988.
118
BIBLIOGRAFA 119
E.H. WATANABE R.M. STEPHAN M. A. NEW CONCEPTS OF INSTANTANEO-
US ACTIVE Y REACTIVE POWERS IN ELECTRICAL SYSTEMS WITH GENERIC
LOADS. IEEE Trans. on Power Delivery, VOL. 8, NO. 2, PP. 697703, APRIL
1993.
GRAINGER J. J. Y STEVENSON W. D. Anlisis de Sistemas de Potencia. MCGRAW-
HILL, 1996, 1DA. EDITION.
GYUGYI L. DYNAMIC COMPENSATION OF AC TRANSMISSION LINES BY SOLID
STATE SYNCHRONOS VOLTAGE SOURCE. IEEE Transactions on Power Deli-
very, VOL. 9, NO. 2, PP. 904911, APRIL 1994.
H. FUJITA S. TOMINAGA H. A. ANALYSIS AND DESIGN OF A DC VOLTAGE
CONTROLLED STATIC VAR COMPENSATOR USING QUAD SERIES VOLTAGE
SOURCE INVERTERS. IEEE Trans. on Industry Applications, VOL. 32, NO. 4,
JULY/AUGUST 1996.
KUNDUR P. Power System Stability and Control. MCGRAW-HILL, INC., 1994, 1ST.
EDITION.
MEJA L. Herramienta software para el modelado y simulacin de sistemas elctricos
de potencia ante grandes perturbaciones. Diseo y elaboracin. TESIS PREGRA-
DO, ESCUELA DE INGENIERA ELCTRICA, ELECTRNICA Y TELECOMUNI-
CACIONES, UNIVERSIDAD INDUSTRIAL DE SANTANDER, ABRIL 2005.
N. MOHAN T.M. UNDELAND W. R. Power Electronics. JOHN WILEY AND SONS
INC., NEW YORK, 1995.
NARAIN G. HINGORANI L. G. UNDERSTANDING FACTS Concepts And Tecnology
Of Flexible AC Transmission Systems. IEEE PRESS, 2000, 1ST. EDITION.
N.G.HINGORANI. FLEXIBLE AC TRANSMISSION. IEEE Spectrum, VOL. 30, NO.
4, PP. 4045, APRIL 1988.
N.G.HINGORANI. FACTS-FLEXIBLE AC TRANSMISSION SYSTEMS. Works-
hop on the Future in High-Voltage Transmission Systems, Cincinnati, Ohio, NOV.
1990.
BIBLIOGRAFA 120
P. GARCIA GONZALEZ A. G. C. CONTROL SYSTEM FOR A PWM BASED STAT-
COM. IEEE Trans. on Power Delivery, VOL. 15, NO. 4, PP. 11401145,
OCTOBER 2000.
P.M. ANDERSON A. F. Power System Control and Stability. IEEE PRESS, 1994, 1A.
EDITION.
RASHID M. Electronica de Potencia: Circuitos, Dispositivos y Aplicaciones. PEAR-
SON PRENTICE HALL, 2004, 3RA. EDITION.
RODRGUEZ C. Herramienta software para el modelado y simulacin de sistemas
elctricos de potencia ante grandes perturbaciones. Diseo y elaboracin. TESIS
PREGRADO, ESCUELA DE INGENIERA ELCTRICA, ELECTRNICA Y TELE-
COMUNICACIONES, UNIVERSIDAD INDUSTRIAL DE SANTANDER, DICIEMBRE
2003.
S. KINCIC A. CHANDRA S. B. MULTILEVER INVERTER AND ITS LIMITATIONS
WHEN APLIED AS STATCOM.
SCHAUDER C. Y MEHTA H. VECTOR ANALYSIS AND CONTROL OF ADVAN-
CED STATIC VAR COMPENSATORS. Proc. Inst. Elect. Eng., VOL. 140, NO. 4,
PP. 299306, JULY 1993.
S.MORI K.MATSUNO M. M. DEVELOPMENT OF A LARGE VAR GENERATOR
USING SELF-COMMUTATED INVERTERS FOR IMPROVING POWER SYSTEM STA-
BILITY. IEEE Trans. Power Delivery, VOL. 8, NO. 1, PP. 371 377, FEB.
1993.
VARGAS H. R. Control Digital de un Inversor Trifsico. TESIS DE MAESTRA, ES-
CUELA DE INGENIERA ELCTRICA, ELECTRNICA Y TELECOMUNICACIO-
NES, UNIVERSIDAD INDUSTRIAL DE SANTANDER, DICIEMBRE 1990.
VAS P. Vector Control of AC Machines. OXFORD:CLAREDON PRESS, 1994,
1ST. EDITION.
WILLIAM H. HAYT JR. J. K. Anlisis de Circuitos en Ingeniera. MC GRAW HILL,
1998, 5TA. EDITION.
BIBLIOGRAFA 121
Y. SUMI Y. HARUMOTO T. H. M. Y. K. I. NEW STATIC VAR CONTROL USING
FORCE COMMUTATED INVERTERS. IEEE Transactions on Power Apparatus
and Systems, VOL. PAS 100, NO. 9, PP. 42164224, SEPTEMBER 1981.
Apndice A
Ecuaciones diferenciales para incluir
en el HAGP
A.1. Solucin en estado estable
A continuacin se presentan las ecuaciones para la solucin del ujo de cargas de un
sistema de potencia con la inclusin del STATCOM.
Para la solucin de un sistema de ecuaciones no-lineales, tradicionalmente se emplea
el mtodo iterativo de Newton-Raphson, que a continuacin se presenta en forma com-
pacta:
x(i +1) = x(i) +J
1
(i)y f [x(i)] (A.1)
donde el error viene dado por la relacin:
=[x(i +1) x(i)[
De acuerdo con el circuito equivalente mostrado en la gura A.1, se deducen las si-
guientes expresiones:
P
stat
=V
2
GkV
dc
VGcos() +kV
dc
VBsen() (A.2)
Q
stat
=V
2
BkV
dc
VGsen() +kV
dc
VBcos() (A.3)
122
A.1 Solucin en estado estable 123
+
-
V
P
L
+ jQ
L
P
stat
+ jQ
stat
Z = R+ jX
kV
dc
C
R
c
P
dc
Figura A.1: Circuito equivalente del STATCOM
+
-
+
+
K
v
(1+sT
2
)
V
re f
V
m
o
m
sT
1
v
m1 v
m2
Figura A.2: Tcnica de control PWM para control de tensin
De acuerdo a la gura A.2, se pueden deducir las ecuaciones que denen la dinmica
del controlador:
v
m1
t
=
1
T
2
[K
v
(V
re f
V) v
m1
] (A.4)
v
m1
t
==
1
T
1
(v
m1
v
m2
) (A.5)
Ahora bien, como se dijo anteriormente que igualando a cero las ecuaciones dinmi-
cas del controlador, se pueden encontrar las ecuaciones para estado estable de dicho
controlador, por tanto se tiene que:
f
1
=m+K
v
(V
re f
V) +m
0
= 0 (A.6)
A.2 Simulacin dinmica 124
Adems, para que se cumpla el balance de potencia, se tiene que:
f
2
= P
stat
P
perdidas
= P
dc
; P
perdidas=
G
c
V
2
dc
+RI
2
(A.7)
El solucin de ujos de potencias de Newton-Raphson se determina a continuacin:
_
_
J
11
J
12
0 0
J
21
J
22
0 0
0
f
1
V
f
1
k
0
0 0
f
2
k
f
2
_
_
[V[
k
_
=
_
_
P
Q
f
1
f
2
_
_
(A.8)
P
j
j
=
P
j
j
+
P
stat
j
Q
j
j
=
Q
j
j
+
Q
stat
j
P
dc
j
=
P
dc
j
+
P
dc
j
P
j
V
j
=
P
j
V
j
+
P
stat
V
j
Q
j
V
j
=
Q
j
V
j
+
Q
stat
V
j
P
dc
V
j
=
P
dc
V
j
+
P
dc
V
j
f
1
k
=
f
1
k
+
f
1
k
P
j
k
=
P
j
k
+
P
stat
k
Q
j
k
=
Q
j
k
+
Q
stat
k
P
dc
k
=
P
dc
k
+
P
dc
k
f
1
V
j
=
f
1
V
j
+
f
1
V
j
P
j
=
P
j
+
P
stat
Q
j
=
Q
j
+
Q
stat
P
dc
=
P
dc
+
P
dc
Los elementos del jacobiano que se encuentran en cero, algunos son reemplazados por
las derivadas anteriores.
Con la inclusin del STATCOM en un sistema de potencia el orden de la matriz jaco-
biana J, se aumenta en dos, tal y como se muestra en la ecuacin (A.8).
A.2. Simulacin dinmica
Un algoritmo de integracin numrica de un conjunto de ecuaciones diferenciales cal-
cula la evolucin en el tiempo de las variables de estado en instantes de tiempo discreto
(x
o
, ..., x
k
...), separados un intervalo de tiempo t.
La herramienta software utilizada (HAGP), utiliza un mtodo de integracin implcito
llamado regla trapezoidal [Kundur, 1994]. La regla trapezoidal calcula las variables de
estado en el paso n+1 por medio de la expresin:
x
n+1
= x
n
+
t
2
[ x
n
+ x
n+1
] (A.9)
Donde x
n
y v
n
son los valores del vector de estado y del vector de tensiones en t = t
n
,
y x
n+1
y v
n+1
son los valores del vector de estado y del vector de tensiones en t =t
n+1
.
Para t = t
n+1
tambin se debe satisfacer la ecuacin de red:
i(x
n+1
, v
n+1
) =Yv
n+1
(A.10)
El clculo de las variables de estado en el paso n +1 por medio de la ecuacin (A.9)
implica la resolucin de un sistema de ecuaciones no lineales. Ello se puede apreciar al
sustituir en la ecuacin (A.9) x por f(x):
x
n+1
= x
n
+
t
2
[f(x
n+1
, v
n+1
) +f(x
n
, v
n
)] (A.11)
Si se denen los vectores F y G como:
F(x
n+1
, v
n+1
) = x
n+1
x
n
t
2
[f(x
n+1
, v
n+1
) +f(x
n
, v
n
)] (A.12)
G(x
n+1
, v
n+1
) = Yv
n+1
i(x
n+1
, v
n+1
) (A.13)
los valores desconocidos x
n+1
y v
n+1
deben satisfacer la condicin:
F(x
n+1
, v
n+1
) = 0 (A.14)
G(x
n+1
, v
n+1
) = 0 (A.15)
Las ecuaciones (A.14) y (A.15) son ecuaciones algbricas lineales. Aplicando el mto-
do de Newton para solucionar estas ecuaciones:
_
x
k+1
n+1
v
k+1
n+1
_
=
_
x
k
n+1
v
k
n+1
_
+
_
x
k
n+1
v
k
n+1
_
(A.16)
Los valores de x
k
n+1
y v
k
n+1
son hallados a partir de:
_
F(x
k
n+1
, v
k
n+1
)
G(x
k
n+1
, v
k
n+1
)
_
=
_
F
x
F
v
G
x
G
v
__
x
k
n+1
v
k
n+1
_
(A.17)
Los valores del jacobiano J en la ecuacin anterior se hallan para x
k
n+1
y v
k
n+1
, y este
tiene la siguiente estructura:
J =
_
F
x
F
v
G
x
G
v
_
=
_
A
D
B
D
C
D
(Y+Y
D
)
_
(A.18)
Las matrices A
D
, B
D
, C
D
y Y
D
dependen de los modelos de los elementos dinmi-
cos del sistema y las cargas estticas no lineales. Estas matrices tienen las siguientes
composiciones:
A
D
=
_
_
A
d1
0 0
0 A
d2
0
.
.
.
.
.
.
.
.
.
.
.
.
0 0 A
dm
_
_
B
D
=
_
_
B
d1
B
d2
.
.
.
B
dm
_
_
Y
D
=
_
_
Y
d1
0 0
0 Y
d2
0
.
.
.
.
.
.
.
.
.
.
.
.
0 0 Y
dN
_
_
C
D
=
_
C
d1
C
d2
C
dm
_
(A.19)
donde m es el nmero de elementos dinmicos del sistema y N el nmero de barras.
Las ecuaciones (A.14) y (A.15) pueden expresarse en trminos de las matrices como:
A
D
x
k
n+1
+B
D
v
k
n+1
= F(x
k
n+1
, v
k
n+1
) (A.20)
C
D
x
k
n+1
+(Y+Y
D
)v
k
n+1
= G(x
k
n+1
, v
k
n+1
) (A.21)
A partir de estas ecuaciones es posible hallar x
k
n+1
y v
k
n+1
a partir de:
x
k
n+1
= A
1
D
_
F
k
n+1
+B
D
v
k
n+1
_
(A.22)
(Y+Y
D
C
D
A
1
D
B
D
)v
k
n+1
= G
k
n+1
+C
D
A
1
D
F
k
n+1
(A.23)
Finalmente teniendo estos valores se hallan x
k+1
n+1
y v
k+1
n+1
a partir de la ecuacin A.16.
Hay que tener en cuenta que el procedimiento anterior es vlido cuando las funciones
A.3 Ecuaciones dinmicas del modelo del STATCOM 127
F y G son diferenciables. Cuando ocurre una discontinuidad tal como una falla u ope-
racin de interruptores, solamente se resuelven las ecuaciones de red para corregir los
valores del vector v en el instante en que ocurre la perturbacin, el vector de estado x
no se cambia. Luego se contina con el proceso de integracin.
A.3. Ecuaciones dinmicas del modelo del STATCOM
Las ecuaciones (A.4) y (A.5) son las ecuaciones dinmicas para el control del ndice de
modulacin. Por consiguiente, las ecuaciones que denen este controlador son:
v
m1
t
=
1
T
2
[K
v
(V
re f
V) v
m1
] (A.24)
v
m2
= K
v
T
1
T
2
(V
re f
V) +(1
T
1
T
2
)v
m2
(A.25)
Las ecuaciones que denen el controlador del ngulo , se relatan a continuacin:
v
d1
= K
p
(V
dcre f
V
dc
) (A.26)
v
d2
t
= K
I
(V
dcre f
V
dc
) (A.27)
La ecuacin dinmica de STATCOM es:
V
dc
=
VI
CV
dc
cos()
G
c
C
V
dc
R
C
I
2
V
dc
(A.28)
Las ecuaciones diferenciales (A.24), (A.27) y (A.28), representan las variables de esta-
do que el modelo aporta para la solucin del sistema .
Los vectores F se denen a continuacin:
F
1
= v
dc1
v
dc2
dt
2C
__
P
ac1
v
dc1
+
P
ac2
v
dc2
_
G
c
(v
dc1
+v
dc2
)
R
_
I
2
11
v
dc1
I
2
12
v
dc2
__
(A.29)
A.3 Ecuaciones dinmicas del modelo del STATCOM 128
F
2
= v
m11
v
m12
dt
2
1
T
2
_
2K
v
_
V
re f
(V
1
+V
2
)
+(v
m11
v
m12
)
_
(A.30)
F
3
= v
d21
v
d22
dtK
I
2
_
v
dc1
+v
dc2
2V
dcre f
(A.31)
Las componentes real e imaginaria de la corriente inyectada por el STATCOM al siste-
ma AC estn dadas por:
I
R
= GV
R
BV
I
kV
DC
(Gcos() Bsen()) (A.32)
I
I
= BV
R
+GV
I
kV
DC
(Bcos() +Gsen()) (A.33)
Las componentes real e imaginaria de la corriente inyectada por el STATCOM al siste-
ma se utilizan para hallar el vector G, el cual se expresa como:
G(x
n+1
, v
n+1
) = Yv
n+1
i(x
n+1
, v
n+1
)
Los elementos del jacobiano J de la solucin dinmica, los cuales corresponden a las
matrices A
D
, B
D
, C
D
y Y
D
, se encuentran calculando las ecuaciones diferenciales de
los vectores F y G, pertinentes para cada matriz, tal y como se explic con la ecuacin
(A.18).
Apndice B
Gua de usuario del HAGP
B.1. Introduccin
Este apndice contiene una gua de usuario del paquete HAGP (Herramienta Software
para el Anlisis de Estabilidad de Gran Perturbacin de Sistemas Elctricos de Poten-
cia). HAGP es un conjunto de funciones MATLAB para el anlisis de estabilidad de
ngulo y de tensin de sistemas elctricos de potencia. Este software corresponde a un
grn trabajo realizado por el Ing. Carlos Rodrguez S [Rodrguez, 2003] bajo la direc-
cin del Dr. Hermann Vargas, donde en primera instancia se analizaba la estabilidad de
ngulo de gran perturbacin de los sistemas de potencia elctricos, luego con la ayuda
de la ingeniera Lizette Meja [Meja, 2005], se incluyeron los modelos necesarios para
los estudios de estabilidad de tensin de los sistema elctricos de potencia.
Esta gua contiene exclusivamente los parmetros de los datos de entrada del STAT-
COM. Para obtener informacin de los generadores, su sistema de excitacin y control,
adems de otras instrucciones de simulacin y datos del sistema, consultar [Rodrguez,
2003].
B.2. Instalacin y Ejecucon
HAGP puede ser instalado en PC con versiones de MATLAB 5.3 o superiores, el pro-
ceso de instalacin comprende los siguientes pasos:
1. Descomprimir la carpeta de archivos HAGP 2.3.0 utilizando WinZip.
129
B.3 Datos de entrada 130
2. Copiar la carpeta con el conjunto de funciones HAGP en un directorio determi-
nado.
3. Adicionar al path de MATLAB la carpeta HAGP con todas las subcarpetas.
Para ejecutar HAGP, se escribe (en la ventana de comandos de MATLAB) la siguiente
instruccin:
> >hagp(nombre)
Donde nombre corresponde a un archivo de datos (sin la extensin .m).
Para nuestro caso, escribir > >hagp(facts).
B.3. Datos de entrada
La entrada de datos se realiza por medio de archivos de MATLAB que contienen datos
de los elementos dinmicos como se muestra a continuacin:
statcom = matriz de datos del STATCOM. Cada la de la matriz corresponde a
un conjunto de datos de un STATCOM determinado, cuyos parmetros corres-
ponden a cada columna.
1. Nmero de barra a la cual se conecta el STATCOM.
2. Status inicial del STATCOM. STATCOM en lnea = 1, STATCOM desconectado
= 0.
3. R = resistencia equivalente del transformador de acople, pu.
4. X = reactancia equivalente del transformador de acople, pu.
5. R
C
= resistencia equivalente en paralelo con el condensador, para contabilizar
perdidas de las redes snubber del inversor, pu.
6. C = Valor de capacitancia del lado DC del inversor, pu.
7. V
re f ac
= Tensin de referencia AC, pu.
B.3 Datos de entrada 131
8. V
dcre f
= Tensin de referencia en el lado DC del STATCOM, pu.
9. K
v
= ganancia del regulador de tensin, pu.
10. T
1
= constante de tiempo del controlador del ndice de modulacin, s.
11. T
2
= constante de tiempo del controlador del ndice de modulacin, s.
12. K
p
= ganancia del controlador proporcional para el control del ngulo .
13. K
I
= ganancia del controlador integral para el control del ngulo .
14. k
0
= valor inicial de la constante k, asociada al ndice de modulacin, m.
15.
0
= valor inicial del ngulo de la tensin a la salida del inversor
Notas:
Todas las constantes de tiempo deben ser diferentes de cero.
Para una mejor comprensin del archivo de entrada, ver el archivo de entrada
correspondiente al caso ejemplo dc1a.m, o mejor escribir el comando > >help
dc1a.
Los datos de entrada en pu de los modelos implementados deben estar en la base
del sistema.
Apndice C
Modelos y libreras en ATP
A continuacin se detallarn los diferentes modelos utilizados en las simulacines rea-
lizadas en ATP a traves de su interfaz grca ATPDRAW.
C.1. Caso de simulacin 1
Figura C.1: Caso de simulacin 1
132
C.1 Caso de simulacin 1 133
inv_vs.mod
Modelo que recibe la seal de referencia en
coordenadas y la magnitud de la tensin del
condensador. Genera los tiempos de disparo de la
tcnica de modulacin SVPWM .
MODEL INV_VS
comment ****************************************************************
* *
* ESTE PROGRAMA PERMITE REALIZAR EL CONTROL DE UN INVERSOR *
* POR MEDIO DE LA TCNICA DE MODULACIN DEL VECTOR ESPACIAL SVPWM *
************************************************************* endcomment
-- DEFINICIN DE LAS VARIABLES DE ENTRADA
INPUT VDC -- TENSIN EN EL CAPACITOR DEL INVERSOR (SE RECOMIENDA EMPLEAR UN INTEGRADOR)
INVALF -- COMPONENTE EN EL EJE ALFA DEL VECTOR ESPACIAL DE TENSIN
INVBET -- COMPONENTE EN EL EJE BETA DEL VECTOR ESPACIAL DE TENSIN
-- DEFINICION DE LOS DATOS QUE SON INTRODUCIDOS POR EL USUARIO
-- (AL HACER DOCLE CLICK SOBRE EL MODELO EN ATPDRAW)
DATA FREC -- FRECUENCIA DE LA SEAL DE SALIDA EN HERTZ
POSVS -- NMERO DE POSCICIONES DEL VECTOR ESPACIAL
-- EMPLEADO PARA DETERMINAR EL NMERO DE PULSOS DEL INVERSOR, V
-- ALORES (MLTIPLOS DE 6)
-- DEFINICIN DE LAS VARIABLES A EMPLEAR POR EL MODELO (SE INCLUYEN LAS SALIDAS)
VAR INVTPA -- SEAL DE CONTROL CONECTADA AL TERMINAL POSITIVO DEL INVERSOR DE LA FASE A
INVTNA -- SEAL DE CONTROL CONECTADA AL TERMINAL NEGATIVO DEL INVERSOR DE LA FASE A
INVTPB -- SEAL DE CONTROL CONECTADA AL TERMINAL POSITIVO DEL INVERSOR DE LA FASE B
INVTNB -- SEAL DE CONTROL CONECTADA AL TERMINAL NEGATIVO DEL INVERSOR DE LA FASE B
INVTPC -- SEAL DE CONTROL CONECTADA AL TERMINAL POSITIVO DEL INVERSOR DE LA FASE C
INVTNC -- SEAL DE CONTROL CONECTADA AL TERMINAL NEGATIVO DEL INVERSOR DE LA FASE C
FCMAG -- FACTOR DE CORRECCIN DE LA MAGNITUD (INVERSIN DE TRANSFORMACIN DEL VECTOR ESPACIAL) ((3^0.5)/2)
FCFASE -- FACTOR DE CORRECCIN DE LA FASE (INVERSIN DE TRANSFORMACIN DEL VECTOR ESPACIAL) 30
o
TEMP -- VARIABLE EMPLEADA PARA DETERCAR EL FIN DE LOS TIEMPOS (TPOCO) POR MEDIO DE LA
-- DIFERENCIA CON EL TIEMPO ACTUAL t
POBVS -- POSICIN BSICA DEL VECTOR ESPACIAL VALORES DE 1 A 6
POIVS -- POSICIN INTERMEDIA DEL VECTOR ESPACIAL, ACORDE A LOS CINCO VALORES DE TEMORIZACIN
-- TPOCO QUE INDICA CUANTO TIEMPO DEBE DURARA EN LAS POSICIONES BSICAS
J -- VARIABLE EMPLEADA PARA CICLOS FOR
-- VARIABLES QUE SON EMPLEADAS PARA EL USO DEL PROGRAMA TIEMPO
-- ENTRADAS
ANGVS -- DEL VECTOR ESPACIAL PARA GENERAR LAS TENSIONES DEL SISTEMA TRIFSICO
MAGVS -- MAGNITUD DEL VECTOR ESPACIAL PARA GENERAR LAS TENSIONES DEL SISTEMA TRIFSICO
VDCINV -- TENSIN EN EL CAPACITOR DEL INVERSOR
TTVS -- TIEMPO TOTAL EN QUE DEBE UBICARSE EL VECTOR ESPCIAL EN UNA POSICIN (T1+T2+T0=1/(FRECUENCIA*POSVS))
MAGVS1
-- SALIDAS
TPOCO[1..5] -- TIEMPO DE LAS POSICIONES COMPUESTAS (DOS POSICIONES A TRAVES DEL TIEMPO PARA
-- OBTENER UNA EQUIVALENTE) DEL VECTOR ESPACIAL CALCULADA POR MEDIO DEL SUBPROGRAMA TIEMPO
OUTPUT INVTPA -- SEAL DE CONTROL CONECTADA AL TERMINAL POSITIVO DEL INVERSOR DE LA FASE A
INVTNA -- SEAL DE CONTROL CONECTADA AL TERMINAL NEGATIVO DEL INVERSOR DE LA FASE A
INVTPB -- SEAL DE CONTROL CONECTADA AL TERMINAL POSITIVO DEL INVERSOR DE LA FASE B
INVTNB -- SEAL DE CONTROL CONECTADA AL TERMINAL NEGATIVO DEL INVERSOR DE LA FASE B
INVTPC -- SEAL DE CONTROL CONECTADA AL TERMINAL POSITIVO DEL INVERSOR DE LA FASE C
INVTNC -- SEAL DE CONTROL CONECTADA AL TERMINAL NEGATIVO DEL INVERSOR DE LA FASE C
------------------------------------------------------------------
-- MODELO EMPLEADO PARA CALCULAR LOS TIMEPOS DEL VECTOR ESPACIAL--
------------------------------------------------------------------
MODEL TIEMPO
INPUT ANGVS -- ANGULO DEL VECTOR ESPACIAL PARA GENERAR LAS TENSIONES DEL SISTEMA TRIFSICO
MAGVS -- MAGNITUD DEL VECTOR ESPACIAL PARA GENERAR LAS TENSIONES DEL SISTEMA TRIFSICO
VDCINV -- TENSIN EN EL CAPACITOR DEL INVERSOR
TTVS -- TIEMPO TOTAL EN QUE DEBE UBICARSE EL VECTOR ESPCIAL EN UNA POSICIN (T1+T2+T0=1/(FRECUENCIA*POSVS))
OUTPUT TPOCO[1..5] -- TIEMPO DE LAS POSICIONES COMPUESTAS (DOS POSICIONES A TRAVES DEL TIEMPO PARA
-- OBTENER UNA EQUIVALENTE) DEL VECTOR ESPACIAL CALCULADA POR MEDIO DEL SUBPROGRAMA TIEMPO
VAR V1 -- PROYECCIN DEL VECTOR ESPCIAL SOBRE LA POSICIN BSICA 1 (0 GRADOS)
V2 -- PROYECCIN DEL VECTOR ESPCIAL SOBRE LA POSICIN BSICA 2 (60 GRADOS)
T1 -- TIEMPO EN QUE SE POSICIONA EL VECTOR ESPACIAL EN LA POSICIN BSICA 1
TPOCO[1..5] -- VARIABLE DE SALIDA
ANG1 -- NGULO COMPRENDIDO ENTRE 0
o
Y 60
o
K
EXEC
ANG1:=ANGVS -- SE AJUSTA EL NGULO A VALORES POSITIVOS (ENTRE 0
o
Y 360
o
)
-- SUBRUTINA EMPLEADA PARA MODIFICAR EL NGULO A UN VALOR ENTRE 0 Y 60
o
FOR K:=6 TO 1 BY -1 DO
IF ANG1>=(60*PI*K/180) THEN
ANG1:=ANG1-(60*PI*K/180)
ENDIF
ENDFOR
V1:=2*MAGVS*SIN(60*PI/180-ANG1)/(3**0.5)
V2:=2*MAGVS*SIN(ANG1)/(3**0.5)
T1:=V1/VDCINV*TTVS
T2:=V2/VDCINV*TTVS
T0:=TTVS-T1-T2
TPOCO[1]:=T1/2
TPOCO[5]:=T1/2
TPOCO[2]:=T2/2
TPOCO[4]:=T2/2
TPOCO[3]:=T0
ENDEXEC
ENDMODEL
-- A CONTINUACIN SE MUESTRAN LOS PULSOS (CONEXIN A LOS TERMINALES POSITIVO 1 O NEGATIVO 0
-- DE LAS FASES DE ACUERDO A LA TCNICA DE MODULACIN SVPWM POR MEDIO DE UNA TABLA DE VALORES
-----------------------------------------------------------------
--FUNCIN PARA DETERMINAR LOS DISPAROS DE LOS INVERSORES FASEA
-----------------------------------------------------------------
FUNCTION DISPAROA POINTLIST (1,1), (2,1), (3,1), (4,1), (5,1),
(6,1), (7,0), (8,0), (9,0), (10,1),
(11,0),(12,0),(13,1),(14,0),(15,0),
(16,0),(17,0),(18,0),(19,0),(20,0),
(21,0),(22,1),(23,1),(24,1),(25,0),
(26,1),(27,1),(28,0),(29,1),(30,1)
-----------------------------------------------------------------
--FUNCIN PARA DETERMINAR LOS DISPAROS DE LOS INVERSORES FASE B
-----------------------------------------------------------------
FUNCTION DISPAROB POINTLIST (1,0), (2,1), (3,1), (4,1), (5,0),
(6,1), (7,1), (8,0), (9,1), (10,1),
(11,1),(12,1),(13,1),(14,1),(15,1),
(16,1),(17,0),(18,0),(19,0),(20,1),
(21,0),(22,0),(23,1),(24,0),(25,0),
(26,0),(27,0),(28,0),(29,0),(30,0)
------------------------------------------------------------------
--FUNCIN PARA DETERMINAR LOS DISPAROS DE LOS INVERSORES FASE C
-----------------------------------------------------------------
FUNCTION DISPAROC POINTLIST (1,0), (2,0), (3,1), (4,0), (5,0),
(6,0), (7,0), (8,0), (9,0), (10,0),
(11,0),(12,1),(13,1),(14,1),(15,0),
(16,1),(17,1),(18,0),(19,1),(20,1),
(21,1),(22,1),(23,1),(24,1),(25,1),
(26,1),(27,0),(28,0),(29,0),(30,1)
-- EN PRIMER LUGAR SE INICIALIZAN LAS VARIABLES DEL MODELO
INIT
TTVS:=1/(FREC*POSVS)
TEMP:=0
TPOCO[1..5]:=0
FCMAG:=(3**0.5)/2
FCFASE:=-30*PI/180
POBVS:=1
POIVS:=5
INVTPA:=1
INVTNA:=0
INVTPB:=0
INVTNB:=1
INVTPC:=1
INVTNC:=0
ENDINIT
EXEC
IF (((t-TEMP)>=TPOCO[POIVS]) AND (t>=0.003))
THEN
-- CLCULO DEL NGULO DEL VECTOR ESPACIAL (valores entre 0
o
y 360
o
)
ANGVS:=ATAN(INVBET/(INVALF))
IF ((INVALF>=0) AND (INVBET<=0))
THEN
ANGVS:=ANGVS+2*PI
ELSE
IF INVALF<=0 THEN
ANGVS:=ANGVS+PI
ENDIF
ENDIF
--LUEGO SE DETERMINA LA MAGNITUD DEL VECTOR ESPCIAL
MAGVS:=((INVALF**2)+(INVBET**2))**0.5
-- LUEGO SE EMPLEAN LOS FACTORES DE CORRECCIN PARA LLAMAR AL PROGRAMA DEL
-- CLCULO DE LOS TIEMPOS
MAGVS:=((3**0.5)/2)*MAGVS
ANGVS:=ANGVS-PI/6
IF (ANGVS>2*PI)
THEN
ANGVS:=ANGVS-2*PI
ENDIF
IF (ANGVS<(0))
THEN
ANGVS:=ANGVS+2*PI
ENDIF
-- SE DETECTA EL FIN DE CADA TEMPORIZACIN Y SE PROCEDE A REALIZAR LA
-- CONMUTACIN AL SIGUIENTE ESTADO
IF ((POIVS=5) OR ((TPOCO[POIVS+1]<=TIMESTEP) AND (POIVS=4)))
THEN
POIVS:=1
-- SE DETECTA LA POSICIN BSICA EN LA QUE SE ENCUENTRA EL VECTOR ESPACIAL 1 - 6
FOR J:=1 TO 6 BY 1
DO
IF (((J-1)*60*PI/180<ANGVS) AND (ANGVS<=(J*60*PI)))
THEN
POBVS:=J
ENDIF
ENDFOR
VDCINV:=VDC
-----------------------------------------
-- LLAMADO DE LA RUTINA PARA CLCULO NORMAL DE TIEMPOS
-----------------------------------------
USE TIEMPO AS TIEMPO ITERATE
INPUT VDCINV:=VDCINV
MAGVS:=MAGVS
ANGVS:=ANGVS
TTVS:=TTVS
OUTPUT TPOCO[1..5]:=TPOCO[1..5]
ENDUSE
-----------------------------------------
ELSE
-- SI UN TIEMPO DE CONMUTACIN ES INFERIOR AL PASO DE SIMULACIN SE
-- ELIMINA LA CONSMUTACIN SE SE PASA A LA SIGUIENTE POSICIN INTERMEDIA
IF ((TPOCO[POIVS+1]<=TIMESTEP) AND (POIVS<=3))
THEN
POIVS:=POIVS+2
ELSE
POIVS:=POIVS+1
ENDIF
ENDIF
TEMP:=t
INVTPA:=DISPAROA(5*(POBVS-1)+POIVS)
INVTNA:=1-DISPAROA(5*(POBVS-1)+POIVS)
INVTPB:=DISPAROB(5*(POBVS-1)+POIVS)
INVTNB:=1-DISPAROB(5*(POBVS-1)+POIVS)
INVTPC:=DISPAROC(5*(POBVS-1)+POIVS)
INVTNC:=1-DISPAROC(5*(POBVS-1)+POIVS)
ENDIF
-- write( VCAPACITOR :, VDC)
ENDEXEC
ENDMODEL
iqcal.mod
Este modelo hace la transformacin de Park de
tensiones o corrientes.
MODEL IQCAL
-- EN ESTA SECCIN SE DEFINEN LAS ENTRADAS DEL MODELO,DEFINIDAS EN EL GRFICO *.MOD
INPUT IFASEA -- CORRIENTE DE LA FASE A
IFASEB -- CORRIENTE DE LA FASE B
IFASEC -- CORRIENTE DE LA FASE C
TETA
-- EN ESTA SECCIN SE DEFINAN LAS VARIABLES DEL MODELO
-- QUE CORRESPONDEN EN NUESTRO CASO A LAS SALIDAS Y LAS TENSIONES DE LNEA
VAR CORD
CORQ
-- EN ESTA SECCIN SE DEFINEN LAS SALIDAS DEL MODELO DEFINIDAS EN EL GRFICO *.MOD
OUTPUT CORD
CORQ
-- EN PRIMER LUGAR SE INICIALIZAN LAS SALIDAS, PARA QUE EL MODELO NO PRESENTE ERRORES
-- DURANTE LA SIMULACIN PORQUE NO SE CONOCEN LOS VALORES DE LA VARIABLE
INIT
CORQ:=0
CORD:=0
ENDINIT
EXEC
CORD:=(2/3)*(IFASEA*cos(TETA)+IFASEB*cos(TETA-2*PI/3)
+IFASEC*cos(TETA+2*PI/3))
CORQ:=(2/3)*(IFASEA*sin(TETA)+IFASEB*sin(TETA-2*PI/3)
+IFASEC*sin(TETA+2*PI/3))
--WRITE( ANGULO, TETA*180/pi)
--WRITE( IQ, CORQ)
--WRITE( IFAA, IFASEA)
ENDEXEC
ENDMODEL
Figura C.2: Caso de simulacin 2
invcav.mod
Inversor trifsico con interruptores ideales
controlados por TACS. Por cada interruptor, existe
un circuito snubber tipo RC con diodo antiparalelo
BEGIN NEW DATA CASE --NOSORT--
DATA BASE MODULE
$ERASE
ARG,FASEXA,FASEXB,FASEXC,VDCPOS,VDCNEG
ARG,GATE11,GATE12,GATE13,GATE14,GATE15,GATE16
ARG,RSNUBB,CSNUBB,VDIODE
NUM,RSNUBB,CSNUBB,VDIODE
DUM,DWFS11,DWFS12,DWFS13,DWFS14,DWFS15,DWFS16
DUM,SWTX11,SWTX13,SWTX15,SWTX12,SWTX14,SWTX16
DUM,SWPS11,SWPS13,SWPS15,SNUB11,SNUB13,SNUB15
DUM,SNUB12,SNUB14,SNUB16,DIOD11,DIOD12,DIOD13
DUM,DIOD14,DIOD15,DIOD16,SWPS12,SWPS14,SWPS16
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
C < n 1>< n 2><ref1><ref2>< R >< A >< B ><Leng><><>0
DWFS14FASEXA .001 0
DWFS11SNUB11 CSNUBB 0
SWPS14SNUB14 CSNUBB 0
SWPS11SNUB11 RSNUBB 0
DWFS14SNUB14 RSNUBB 0
FASEXADWFS11 .001 0
FASEXBDWFS13 .001 0
FASEXCDWFS15 .001 0
DWFS12DIOD12 .001 0
DWFS14DIOD14 .001 0
DWFS16DIOD16 .001 0
SWTX11SWPS11 .001 0
SWTX13SWPS13 .001 0
SWTX15SWPS15 .001 0
SWPS12SWTX12 .001 0
SWPS14SWTX14 .001 0
SWPS16SWTX16 .001 0
VDCPOSSWPS11 .001 0
VDCPOSSWPS13 .001 0
VDCPOSSWPS15 .001 0
VDCNEGSWPS12 .001 0
VDCNEGSWPS14 .001 0
VDCNEGSWPS16 .001 0
DWFS11DIOD11 .001 0
DWFS13DIOD13 .001 0
DWFS15DIOD15 .001 0
FASEXBDWFS16 .001 0
FASEXCDWFS12 .001 0
/SWITCH
C < n 1>< n 2>< Tclose ><Top/Tde >< Ie ><Vf/CLOP >< type >
13DWFS11SWTX11 GATE11 0
13SWTX16DWFS16 GATE16 0
11DIOD11SWPS11 VDIODE 0
11SWPS14DIOD14 VDIODE 0
/SOURCE
BEGIN NEW DATA CASE
C
$PUNCH,INVCAV.LIB
BEGIN NEW DATA CASE
BLANK
logic1.atp
Lgica de disparo de los dispositivos de
conmutacin. Recibe la seal de sincronismo, un
desfasaje y monta los disparos de los interruptores
del inversor.
DATA BASE MODULE
$ERASE
ARG,WTXXXX,DELTAX,STTSWT,GATEUA,GATEUB,GATEUC,GATEDA,GATEDB,GATEDC
ARG,VALORX
NUM,VALORX
DUM,VAGATE,VBGATE,VCGATE,GATE11,GATE13,GATE15
/TACS
98VAGATE =SIN(WTXXXX-PI/2+DELTAX)*STTSWT
98VBGATE =SIN(WTXXXX-PI/2+DELTAX-2*PI/3)*STTSWT
98VCGATE =SIN(WTXXXX-PI/2+DELTAX+2*PI/3)*STTSWT
98GATE1160+ZERO +ZERO +UNITY VAGATEZERO
98GATE1360+ZERO +ZERO +UNITY VBGATEZERO
98GATE1560+ZERO +ZERO +UNITY VCGATEZERO
98GATEUA =GATE11*STTSWT
98GATEDA =( .NOT. (GATE11))*STTSWT
98GATEUB =GATE13*STTSWT
98GATEDB =( .NOT. (GATE13))*STTSWT
98GATEUC =GATE15*STTSWT
98GATEDC =( .NOT. (GATE15))*STTSWT
C 1 2 3 4 5 6 7 8
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
/SWITCH
/SOURCE
BEGIN NEW DATA CASE
C
$PUNCH,LOGIC1.LIB
BEGIN NEW DATA CASE
BLANK
ruptor.atp
Mdulo que representa un interruptor tripolar que
cierra los tres polos al recibir una seal de control.
DATA BASE MODULE
$ERASE
ARG,VA0000,VB0000,VC0000,VA0003,VB0003,VC0003,CLOSFA
DUM,VA0001,VB0001,VC0001,VA0002,VB0002,VC0002
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
VA0000VA0001 .001 0
VA0001 1.E6 0
VB0000VB0001 .001 0
VB0001 1.E6 0
VC0000VC0001 .001 0
VC0001 1.E6 0
VA0002 1.E6 0
VB0002 1.E6 0
VC0002 1.E6 0
VA0002VA0003 .001 0
VB0002VB0003 .001 0
VC0002VC0003 .001 0
VC0001VC0002 1.E9 0
VB0001VB0002 1.E9 0
VA0001VA0002 1.E9 0
/SWITCH
13VA0001VA0002 CLOSFA 0
13VB0001VB0002 CLOSFA 0
13VC0001VC0002 CLOSFA 0
BEGIN NEW DATA CASE
C
$PUNCH
BEGIN NEW DATA CASE
BLANK
pll.atp
Mdulo que detecta la frecuencia de las tensiones de
fase y genera la seal de sincronismo.
DATA BASE MODULE
C Diode_C. Diode circuit with series RD and LD, and snubber RS and CS.
C Diode points from AN to CA. Series RD and LD between CA and TP
$ERASE
ARG,VAXXXX,VBXXXX,VCXXXX,WXXXXX,WTTETA
ARG,VBASEX,GPIPPP,GPIPPI
NUM,VBASEX,GPIPPP,GPIPPI
DUM,P10000,P20000,IAREFX,ICREFX,P00000,RESETX
DUM,PPI000,OUTPPP,OUTPPI,CONDTN,C1XXXX,ZOH000
/TACS
98CONDTN =2*pi
1OUTPPI +PPI000 1.
1.
1.
0ZOH000 +RESETX 1.
98WTTETA58 +WXXXXX 1. 1.ZOH000C1XXXX
98P10000 =(VAXXXX-VBXXXX)*IAREFX/((VBASEX)*SQRT(2/3))
98P20000 =(VCXXXX-VBXXXX)*ICREFX/((VBASEX)*SQRT(2/3))
98OUTPPP =GPIPPP*P00000
98P00000 =P10000+P20000
98PPI000 =GPIPPI*P00000
98IAREFX =SIN(WTTETA)
98ICREFX =SIN(WTTETA+2*PI/3)
98C1XXXX =0
98RESETX60+UNITY +ZERO +ZERO WTTETACONDTN
0WXXXXX +OUTPPI +OUTPPP 1.
C 1 2 3 4 5 6 7 8
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
/SWITCH
/SOURCE
BEGIN NEW DATA CASE
C
$PUNCH,PLLCV1.LIB
BEGIN NEW DATA CASE
BLANK
tranyy.atp
Transformador trifsico compuesto por tres
transformadores monofsicos, con una conguracin
estrella-estrella. Por el lado secundario se encuentra
abierto para permitir conexiones.
DATA BASE MODULE
$ERASE
ARG,VA0001,VB0001,VC0001,VA0002,VB0002,VC0002,VASCRO,VBSCRO,VCSCRO
ARG,TNPRIM,TNSCND,LLPRIM,LLSCND,RRPRIM,RRSCND,RSTMAG
NUM,TNPRIM,TNSCND,LLPRIM,LLSCND,RRPRIM,RRSCND,RSTMAG
DUM,TX0001,TX0002,TX0003,GNDPRM,GNDSCD
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
TRANSFORMER TX0001RSTMAG 0
9999
1VA0001GNDPRM RRPRIMLLPRIMTNPRIM
2VA0002VASCRO RRSCNDLLSCNDTNSCND
9999
1VB0001GNDPRM RRPRIMLLPRIMTNPRIM
2VB0002VBSCRO RRSCNDLLSCNDTNSCND
9999
1VC0001GNDPRM RRPRIMLLPRIMTNPRIM
2VC0002VCSCRO RRSCNDLLSCNDTNSCND
GNDSCD 1.E6 0
GNDPRM 1.E6 0
TX0001 1.E6 0
TX0002 1.E6 0
TX0003 1.E6 0
VA0002TX0001 1.E6 0
VB0002TX0002 1.E6 0
VC0002TX0003 1.E6 0
TX0001VA0001 1.E6 0
TX0002VB0001 1.E6 0
TX0003VC0001 1.E6 0
GNDPRMVA0001 1.E6 0
GNDPRMVB0001 1.E6 0
GNDSCDVC0002 1.E6 0
GNDSCDVA0002 1.E6 0
GNDSCDVB0002 1.E6 0
GNDPRMVC0001 1.E6 0
/SWITCH
/SOURCE
C < n 1><>< Ampl. >< Freq. ><Phase/T0>< A1 >< T1 >< TSTART >< TSTOP >
/INITIAL
/OUTPUT
BEGIN NEW DATA CASE
C
$PUNCH,TRANYY.LIB
BEGIN NEW DATA CASE
BLANK
tranyy.atp
Transformador trifsico compuesto por tres
transformadores monofsicos, con una conguracin
delta-estrella. Por el lado secundario se encuentra
abierto para permitir conexiones.
DATA BASE MODULE
$ERASE
ARG,VA0001,VB0001,VC0001,VA0002,VB0002,VC0002,VASCRO,VBSCRO,VCSCRO
ARG,TNPRIM,TNSCND,LLPRIM,LLSCND,RRPRIM,RRSCND,RSTMAG
NUM,TNPRIM,TNSCND,LLPRIM,LLSCND,RRPRIM,RRSCND,RSTMAG
DUM,TX0001,TX0002,TX0003,GNDPRM
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
9999
1VA0001VB0001 RRPRIMLLPRIMTNPRIM
2VA0002VASCRO RRSCNDLLSCNDTNSCND
9999
1VB0001VC0001 RRPRIMLLPRIMTNPRIM
2VB0002VBSCRO RRSCNDLLSCNDTNSCND
9999
1VC0001VA0001 RRPRIMLLPRIMTNPRIM
2VC0002VCSCRO RRSCNDLLSCNDTNSCND
VA0002TX0001 1.E6 0
VB0002TX0002 1.E6 0
VC0002TX0003 1.E6 0
TX0001VA0001 1.E6 0
TX0002VB0001 1.E6 0
TX0003VC0001 1.E6 0
VB0001VA0001 1.E6 0
VC0003VB0001 1.E6 0
VA0001VC0001 1.E6 0
GNDPRMVA0002 1.E6 0
GNDPRMVB0002 1.E6 0
GNDPRMVC0002 1.E6 0
TX0001 1.E6 0
TX0002 1.E6 0
TX0003 1.E6 0
/SWITCH
/SOURCE
/INITIAL
/OUTPUT
BEGIN NEW DATA CASE
C
$PUNCH,TRANDY.LIB
BEGIN NEW DATA CASE
BLANK
gndstr.atp
Este mdulo representa un snubber numrico cuya
funcin es estabilizar oscilaciones conectando las
fases en un puente neutro a travs de resistencias de
1 M.
DATA BASE MODULE
$ERASE
ARG,FASEA1,FASEB1,FASEC1
ARG,FASEA2,FASEB2,FASEC2
ARG,GNDYXX
DUM,DELTAA,DELTAB,DELTAC
/BRANCH
DELTAAGNDYXX 1.E6 0
DELTABGNDYXX 1.E6 0
DELTACGNDYXX 1.E6 0
DELTAAFASEA2 .0001 0
DELTABFASEB2 .0001 0
DELTACFASEC2 .0001 0
FASEA1DELTAA .0001 0
FASEB1DELTAB .0001 0
FASEC1DELTAC .0001 0
BEGIN NEW DATA CASE
C
$PUNCH
BEGIN NEW DATA CASE
BLANK
tr_ac_vs.mod
Este mdulo hace la transformacin de las seales
monofsicas de tensin o corrientes, al sistema de
coordenadas .
MODEL TR_AC_VS
comment ****************************************************************
* *
* ESTE PROGRAMA REALIZA LA TRANSFORMACIN DEL VECTOR ESPACIAL *
* LAS ENTRADAS CORRESPONDEN A LAS TENSIONES DE FASE DE UN SISTEMA *
* TRIFSICO, Y LAS SALIDAS CORRESPONDEN A LA TRANSFORMACON DEL *
* VECTOR ESPACIAL, CON LAS COMPONENTES ALFA, BETA Y CERO *
************************************************************* endcomment
INPUT VFASEA -- TENSIN DE LA FASE A
VFASEB -- TENSIN DE LA FASE B
VFASEC -- TENSIN DE LA FASE C
VAR VALFA -- PROYECCIN DEL VECTOR ESPACIAL SOBRE EL EJE ALFA
VBETA -- PROYECCIN DEL VECTOR ESPACIAL SOBRE EL EJE BETA
VCERO -- COMPONENTE DE SECUENCIA CERO DEL VECTOR ESPACIAL
VAB -- VALOR INSTANTNEO DE LA TENSIN DE FASE AB
VBC -- VALOR INSTANTNEO DE LA TENSIN DE FASE AB
VCA -- VALOR INSTANTNEO DE LA TENSIN DE FASE AB
OUTPUT VALFA -- PROYECCIN DEL VECTOR ESPACIAL SOBRE EL EJE ALFA
INIT
VALFA:=0
VBETA:=0
VCERO:=0
ENDINIT
-- EN ESTA SECCIN SE DEFINE EL PROGRAMA, EL CUAL SER EJECUTADO
-- CADA INTERVALO DE TIEMPO DEFINIDO COMO PASO DE LA SIMULACIN
EXEC
-- SE CALCULAN LAS TENSIONES DE LNEA
VAB:=VFASEA-VFASEB
VBC:=VFASEB-VFASEC
VCA:=VFASEC-VFASEA
-- SE CALCULAN LAS COMPONENTES DEL VECTOR ESPACIAL, SEGUN LIBRO PETER BASS
VALFA:=(2/3)*(VFASEA-0.5*VFASEB-0.5*VFASEC)
VBETA:=(2/3)*(SQRT(3)/2)*(VFASEB-VFASEC)
VCERO:=(2/3)*(VAB+VBC+VCA)
ENDEXEC
ENDMODEL
A continuacin se presenta el listado que corresponde al archivo *.ATP del caso de
simulacin 2.
BEGIN NEW DATA CASE
C --------------------------------------------------------
C Generated by ATPDRAW marzo, mircoles 22, 2006
C A Bonneville Power Administration program
C Programmed by H. K. Hidalen at SEfAS - NORWAY 1994-2002
C --------------------------------------------------------
$DUMMY, XYZ000
C dT >< Tmax >< Xopt >< Copt >
5.E-6 10.
500 1 1 1 1 0 0 1 0
TACS HYBRID
/MODELS
MODELS
/TACS
92TQ -1. 10.
92TITA -1. 10.
90BUS4C -1. 10.
90BUS4A -1. 10.
92OMEGA -1. 10.
92IF1 -1. 10.
2XX0336 +XX0334 750.
1.
.0187 .1062 .0095
1XX0308 +VT4 1.
1.
1. 0.001
90BUS4B -1. 10.
90BUS10B -1. 10.
90BUS10A -1. 10.
90BUS10C -1. 10.
98DELTA =(TITA-OMEGAR*TIMEX-PI/2)*180/PI
99PG =TQ*OMEGA
1XX0340 +XX0336 0.04
1.
1. 1.
98XX0334 = XX0332 - XX0340
98XX0332 = XX0310 - XX0308
98XX0310 =1
98XX0112 = XX0336 + XX0346
98XX0346 =1
90BUS3C -1. 10.
90BUS3A -1. 10.
90BUS3B -1. 10.
98VT10 =(SQRT(BUS10A**2+BUS10B**2+BUS10C**2)/230000)*SQRT(2)
90BUS8B -1. 10.
90BUS8C -1. 10.
90BUS8A -1. 10.
98VF = 1. * XX0112
98XX0496 =TETA3
11STTSWT 1. .505 10.
11XX0502 1. .5 10.
2XX0532 +XX0516 2.1
1.5 .0115
11. .041 3.E-5
98XX0516 = XX0520 - XX0462
98XX0520 =1
98VDC =VDCPOS-VDCNEG
90VDCPOS -1. 10.
90VDCNEG -1. 10.
98XX0518 =-DELT+DEFASE*PI/180
98XX0524 =-DELT-30*PI/180+DEFASE*PI/180
98XX0530 =BUS8A
98VRMS8 66 +XX0530 60.
0DELT +XX0532 1.-0.174 .1745
98DEFASE =90+90-30-10
90TETA3 -1. 10.
98XX0462 =VRMS8*sqrt(6)/230000
33PG
33DELTA
33DELTA
33VT8
33WT1
33DELT
33VDC
33XX0141
33XX0462
33VT4
/MODELS
INPUT
IX0001 {tacs(BUS8A )}
IX0002 {tacs(BUS8B )}
IX0003 {tacs(BUS8C )}
OUTPUT
VA1
VB1
XX0281
TETA2
MODEL TR_AC_VS
comment ****************************************************************
* *
************************************************************* endcomment
INIT
VALFA:=0
VBETA:=0
VCERO:=0
ENDINIT
EXEC
VAB:=VFASEA-VFASEB
VBC:=VFASEB-VFASEC
VCA:=VFASEC-VFASEA
ENDEXEC
ENDMODEL
MODEL ANGUL
comment ****************************************************************
* *
************************************************************* endcomment
INPUT VALFA -- TENSIN DE LA FASE A
VBETA -- TENSIN DE LA FASE B
VAR ANGVS
OUTPUT ANGVS
INIT
ANGVS:=0
ENDINIT
EXEC
-- CLCULO DEL NGULO DEL VECTOR ESPACIAL (valores entre 0
o
y 360
o
)
ANGVS:=ATAN(VBETA/VALFA)
IF ((VALFA>=0) AND (VBETA<=0))
THEN
ANGVS:=ANGVS+2*PI
ELSE
IF VALFA<=0 THEN
ANGVS:=ANGVS+PI
ENDIF
ENDIF
-- WRITE( ANGULO, ANGVS*180/pi)
ENDEXEC
ENDMODEL
USE TR_AC_VS AS TR_AC_VS
INPUT
VFASEA:= IX0001
VFASEB:= IX0002
VFASEC:= IX0003
OUTPUT
VA1 :=VALFA
VB1 :=VBETA
XX0281:=VCERO
ENDUSE
USE ANGUL AS ANGUL
INPUT
VALFA:= VA1
VBETA:= VB1
OUTPUT
TETA2 :=ANGVS
ENDUSE
ENDMODELS
C 1 2 3 4 5 6 7 8
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
TRANSFORMER X0034A 1.E5 0
9999
1X0298AX0298B .005.16838 13.8
2BUS3A .0146.773132.79
TRANSFORMER X0034A X0034B
1X0298BX0298C
2BUS3B
TRANSFORMER X0034A X0034C
1X0298CX0298A
2BUS3C
X0306AX0068A .06613 40. 0
X0306BX0068B .06613 40. 0
X0306CX0068C .06613 40. 0
X0060AX0061A .57322.158 0
X0060BX0061B .57322.158 0
X0060CX0061C .57322.158 0
X0061A .15163 0
X0061B .15163 0
X0061C .15163 0
X0060A .15163 0
X0060B .15163 0
X0060C .15163 0
X0066A 51.5138.307 0
X0066B 51.5138.307 0
X0066C 51.5138.307 0
X0068ABUS3A 3.5046133.01 0
X0068BBUS3B 3.5046133.01 0
X0068CBUS3C 3.5046133.01 0
BUS3A .9073 0
BUS3B .9073 0
BUS3C .9073 0
BUS3A BUS8A 4.8492184.75 0
BUS3B BUS8B 4.8492184.75 0
BUS3C BUS8C 4.8492184.75 0
BUS8A 1.26 0
BUS8B 1.26 0
BUS8C 1.26 0
BUS3A 1.26 0
BUS3B 1.26 0
BUS3C 1.26 0
X0068A .9073 0
X0068B .9073 0
X0068C .9073 0
BUS8A X0083A .57322.158 0
BUS8B X0083B .57322.158 0
BUS8C X0083C .57322.158 0
X0083A .15163 0
X0083B .15163 0
X0083C .15163 0
BUS8A .15163 0
BUS8B .15163 0
BUS8C .15163 0
BUS10A 161.8957.524 0
BUS10B 161.8957.524 0
BUS10C 161.8957.524 0
X0320AX0066A 46.772 0
X0320BX0066B 46.772 0
X0320CX0066C 46.772 0
X0123A .00522 700. 0
X0123B .00522 700. 0
X0123C .00522 700. 0
X0083ABUS10A 46.772 0
X0083BBUS10B 46.772 0
X0083CBUS10C 46.772 0
VDCNEGVDCPOS 4.5E4 0
XX0490 .001 0
TRANSFORMER X0225A 1.E5 0
9999
1VOA VOB .369 12.
2X0486A 136. 132.8
TRANSFORMER X0225A X0225B
1VOB VOC
2X0486B
TRANSFORMER X0225A X0225C
1VOC VOA
2X0486C
NETRY1 1.E8 0
VDCNEGVDCPOS 2.25 0
X0294A 1.9176 2. 0
X0294B 1.9176 2. 0
X0294C 1.9176 2. 0
VOA .00522 1028. 0
VOB .00522 1028. 0
VOC .00522 1028. 0
$INCLUDE, C:\ATP\project\usp\logic1.lib, XX0496, XX0518, STTSWT, XX0348 $$
, XX0354, XX0360, XX0366, XX0372, XX0378, 100.
$INCLUDE, C:\ATP\project\usp\pll.lib, BUS8A#, BUS8B#, BUS8C#, XX0141, WT1### $$
, 1.38E5, 100., 5.E3
$INCLUDE, C:\ATP\project\usp\invcav.lib, XX0420, XX0422, XX0424, VDCPOS $$
, VDCNEG, XX0348, XX0378, XX0354, XX0366, XX0360, XX0372, 500., 0.5 $$
, 0.0
$INCLUDE, C:\ATP\project\usp\logic1.lib, XX0496, XX0524, STTSWT, XX0384 $$
, XX0390, XX0396, XX0402, XX0408, XX0414, 100.
$INCLUDE, C:\ATP\project\usp\invcav.lib, XX0426, XX0428, XX0430, VDCPOS $$
, VDCNEG, XX0384, XX0414, XX0390, XX0402, XX0396, XX0408, 500., 0.5 $$
, 0.0
$INCLUDE, C:\ATP\project\usp\tranyy.lib, XX0450, XX0452, XX0454, VOA### $$
, VOB###, VOC###, XX0444, XX0446, XX0448, 3500., 3500.,0.0314,0.0314 $$
, 0.0, 0.0, 1.E5
$INCLUDE, C:\ATP\project\usp\trandy.lib, XX0438, XX0440, XX0442, XX0444 $$
, XX0446, XX0448, XX0490, XX0490, XX0490, 6063., 3500.,0.0314,0.0314 $$
, 0.0, 0.0, 1.E5
$INCLUDE, C:\ATP\project\usp\gndelt.lib, XX0432, XX0434, XX0436, XX0438 $$
, XX0440, XX0442
$INCLUDE, C:\ATP\project\usp\gndstr.lib, XX0460, XX0458, XX0456, XX0450 $$
, XX0452, XX0454, NETRY1
$INCLUDE, C:\ATP\project\usp\ruptor.lib, XX0420, XX0422, XX0424, XX0460 $$
, XX0458, XX0456, XX0502
, XX0434, XX0436, XX0502
/SWITCH
BUS4A X0298A -1. 10. 0
BUS4B X0298B -1. 10. 0
BUS4C X0298C -1. 10. 0
BUS3A X0060A -1. 4.65 0
BUS3B X0060B -1. 4.65 0
BUS3C X0060C -1. 4.65 0
X0061AX0320A 4.5 10. 0
X0061BX0320B 4.5 10. 0
X0061CX0320C 4.5 10. 0
BUS8A X0486A .5 10. 0
BUS8B X0486B .5 10. 0
BUS8C X0486C .5 10. 0
TETA2 TETA3 MEASURING 0
/SOURCE
14X0306A 0 1.38E5 60. -1. 10.
14X0306B 0 1.38E5 60. -120. -1. 10.
14X0306C 0 1.38E5 60. 120. -1. 10.
60TETA2 0 -1. 10.
59BUS4A 8.E3 60. -5.
BUS4B
BUS4C
PARAMETER FITTING 2.
1 1 1 2 1. 1. 200. 13.8 300.
BLANK
.001096 .15 1.7 1.64 .238324 .378 .18469 .185151
6.19488 .028716 .07496 .1053 .1
1 1.0 2.7113
BLANK
11
21
31
41
51
BLANK
74OMEGA 2
73TQ 14
74TITA 1
71VF
73IF1 4
FINISH
/INITIAL
/OUTPUT
VOA VOB VOC
BLANK TACS
BLANK MODELS
BLANK BRANCH
BLANK SWITCH
BLANK SOURCE
BLANK INITIAL
BLANK OUTPUT
BLANK PLOT
BEGIN NEW DATA CASE
BLANK
La gura C.3 muestra parte del sitema de potencia de un STATCOM con una barra
innita, la cual est representada por una fuente trifasica tipo TACS(60).
La gura C.4 muestra parte del montaje del STATCOM de 48 pulsos, en la cual se
presentan solo dos inversoren de seis pulsos. El mismo montaje es equivalente para los
siguientes seis inversores de seis pulsos.
La gura C.5 muestra como se genera la seal de desfasaje a partir del controlador PI.
Figura C.3: Caso de simulacin 3 - parte 1
medida.atp
Este mdulo sensa las seales de tensin y corriente
por las lneas
DATA BASE MODULE
$ERASE
ARG,VA0001,VB0001,VC0001,VA0003,VB0003,VC0003
ARG,VAXXXX,VBXXXX,VCXXXX,IAXXXX,IBXXXX,ICXXXX
DUM,VA0002,VB0002,VC0002
/TACS
90VA0001 100.
90VB0001 100.
90VC0001 100.
91VA0002 100.
91VB0002 100.
91VC0002 100.
98VAXXXX =VA0001
98IAXXXX =VA0002
98VBXXXX =VB0001
98IBXXXX =VB0002
98VCXXXX =VC0001
98ICXXXX =VC0002
C 1 2 3 4 5 6 7 8
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
VA0001VA0002 .001 0
VB0001VB0002 .001 0
VC0001VC0002 .001 0
VA0001 1.E9 0
VB0001 1.E9 0
VC0001 1.E9 0
VA0002 1.E9 0
VB0002 1.E9 0
VC0002 1.E9 0
/SWITCH
VB0002VB0003 MEASURING 0
VC0002VC0003 MEASURING 0
VA0002VA0003 MEASURING 0
BEGIN NEW DATA CASE
C
$PUNCH,MEDIDA.LIB
BEGIN NEW DATA CASE
BLANK
A continuacin se presenta el listado que corresponde al archivo *.ATP del caso de
simulacin 3.
BEGIN NEW DATA CASE
C --------------------------------------------------------
C Generated by ATPDRAW abril, lunes 17, 2006
C A Bonneville Power Administration program
C Programmed by H. K. Hidalen at SEfAS - NORWAY 1994-2002
C --------------------------------------------------------
$DUMMY, XYZ000
C dT >< Tmax >< Xopt >< Copt >
2.5E-6 1.5
500 1 1 1 1 0 0 1 0
TACS HYBRID
/MODELS
MODELS
/TACS
98WTX =WT1
11STTSWT 1. .25 10.
11CNTINV 1. .25 10.
90VDC -1. 10.
90NETRY3 -1. 10.
98FASEX =DELTAX*STTQRF+CORETO*PI/180
11STTQRF 1. .25 10.
90NETRY5 -1. 10.
23QS1 2. .5 .25 .5 10.
23QS2 1. 1. .5 .5 10.
98XX1378 =FASEX+30*PI/180
98QREFX =(QS1-1)*QS2*STTQRF
98ERROQ1 =(QREF-QPU)*STTQRF
98POSFSA =POSSTT*POSAMP*1000*SIN(WTXXXX+(POSFSE)*PI/180)
11POSSTT 1. .05 10.
98POSFSB =POSSTT*POSAMP*1000*SIN(WTXXXX+(POSFSE)*PI/180-2*PI/3)
98POSFSC =POSSTT*POSAMP*1000*SIN(WTXXXX+(POSFSE)*PI/180+2*PI/3)
98POSFSE =0
1XX1534 +DELTAX 1.
200.
200. 1.
98XX1540 =2.5*PI/180
98VDCPU =VDCPOS/5500
98VDCFIL =VDCFLT/5500
98DELTAG =DELTAX*180/PI
98DELTGF =DX12*180/PI
98XX1510 =FASEX+0
98XX1512 =FASEX+37.5*PI/180
98XX1514 =FASEX+7.5*PI/180
98XX1516 =FASEX+45*PI/180
98XX1518 =FASEX+15*PI/180
98XX1520 =FASEX+52.5*PI/180
98POSAMP =112.67
98WTXXXX =2*PI*60*TIMEX
98XX1538 =-2.5*PI/180
98CORETO =-70.885
98P1 =1.5*(VALFAX*IALFAX+VBETAX*IBETAX)
98Q1 =1.5*(VBETAX*IALFAX-VALFAX*IBETAX)
98QPU =Q1/100E6
90VALFAX -1. 10.
90VBETAX -1. 10.
90IALFAX -1. 10.
90IBETAX -1. 10.
90NETRY1 -1. 10.
98VA12P1 =(VA61-NETRY1)+(VA62-VB62)/SQRT(3)
98VB12P1 = (VB61-NETRY1)+(VB62-VC62)/SQRT(3)
98V24P1 =(VA12P1-VB12P1)+(VA12P2-VB12P2)
98V24P2 =(VA12P3-VB12P3)+(VA12P4-VB12P4)
98VB12P2 =(VB63-NETRY3)+(VB64-VC64)/SQRT(3)
98V48P =V24P1+V24P2
98VDCPOS =VDC
90NETRY7 -1. 10.
11XX1368 1. .05 10.
98V48PU =V48P/48.78E3
98VXXXPU =VXXXXA/112.675E3
98ITXXPU =IXXXXA/418.369
1XX1374 +VDC 1.
500.
500. 1.
1DX12 +XX1534 1.
200.
200. 1.
1XX1536 +QREFX 1.
50.
50. 1.
1QREF +XX1536 1.
50.
50. 1.
1VDCFLT +XX1374 1.
500.
500. 1.
98XX1526 =FASEX+22.5*PI/180
98VPU =SQRT(VXXXXA**2+VXXXXB**2+VXXXXC**2)/1.38E5
33WT1
33XX0044
33IXXXXA
33VDCFIL
33DELTAG
33DELTGF
33VDCPU
33VXXXXA
33V48P
33QREF
33QPU
33VDCPOS
33DELTAX
33P1
33V48PU
33VXXXPU
33ITXXPU
33Q1
33VPU
/MODELS
INPUT
IX0001 {tacs(IXXXXA)}
IX0002 {tacs(IXXXXB)}
IX0003 {tacs(IXXXXC)}
IX0004 {tacs(VXXXXA)}
IX0005 {tacs(VXXXXB)}
IX0006 {tacs(VXXXXC)}
OUTPUT
XX0473
XX0477
XX0153
XX0482
XX0483
XX0159
MODEL ITRAC_VE
comment ****************************************************************
* *
************************************************************* endcomment
INIT
VALFA:=0
VBETA:=0
VCERO:=0
ENDINIT
EXEC
VCERO:=(2/3)*(VFASEA+VFASEB+VFASEC)
ENDEXEC
ENDMODEL
MODEL TR_AC_VS
comment ****************************************************************
* *
************************************************************* endcomment
INIT
VALFA:=0
VBETA:=0
VCERO:=0
ENDINIT
EXEC
VAB:=VFASEA-VFASEB
VBC:=VFASEB-VFASEC
VCA:=VFASEC-VFASEA
ENDEXEC
ENDMODEL
USE ITRAC_VE AS ITRAC_VE
INPUT
VFASEA:= IX0001
VFASEB:= IX0002
VFASEC:= IX0003
OUTPUT
XX0473:=VALFA
XX0477:=VBETA
XX0153:=VCERO
ENDUSE
USE TR_AC_VS AS TR_AC_VS
INPUT
VFASEA:= IX0004
VFASEB:= IX0005
VFASEC:= IX0006
OUTPUT
XX0482:=VALFA
XX0483:=VBETA
XX0159:=VCERO
ENDUSE
ENDMODELS
C 1 2 3 4 5 6 7 8
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
VDC 4591. 1
NETRY1 1.E6 0
NETRY3 1.E6 0
X0001A 1.E6 1. 0
X0001B 1.E6 1. 0
X0001C 1.E6 1. 0
NETRY5 1.E6 0
X0005A 1.E6 1. 0
X0005B 1.E6 1. 0
X0005C 1.E6 1. 0
POSFSAXX0395 .001 0
POSFSBXX0397 .001 0
POSFSCXX0399 .001 0
XX0395XX0401 25.256 0
XX0397XX0403 25.256 0
XX0399XX0405 25.256 0
XX1338 1.E6 0
XX0401 1.E6 0
XX0403 1.E6 0
XX0405 1.E6 0
X0009A 1.E6 1. 0
X0009B 1.E6 1. 0
X0009C 1.E6 1. 0
XX0473IALFAX .001 0
XX0477IBETAX .001 0
XX0482VALFAX .001 0
XX0483VBETAX .001 0
X0017A 1.E6 1. 0
X0017B 1.E6 1. 0
X0017C 1.E6 1. 0
X0025A 1.E6 1. 0
X0025B 1.E6 1. 0
X0025C 1.E6 1. 0
NETRY7 1.E6 0
X0021A 1.E6 1. 0
X0021B 1.E6 1. 0
X0021C 1.E6 1. 0
X0029A 1.E6 1. 0
X0029B 1.E6 1. 0
X0029C 1.E6 1. 0
VDC XX0858 .1 0
VDC XX0854 .1 0
VDC XX1010 .1 0
VDC XX1006 .1 0
X0033A 1.E6 1. 0
X0033B 1.E6 1. 0
X0033C 1.E6 1. 0
X0013A 1.E6 1. 0
X0013B 1.E6 1. 0
X0013C 1.E6 1. 0
X0037A 1.E6 1. 0
X0037B 1.E6 1. 0
X0037C 1.E6 1. 0
VDC XX1160 .1 0
VDC XX1156 .1 0
VDC XX1346 .1 0
VDC XX1330 .1 0
$INCLUDE, C:\ATP\project\usp\pllcav.lib, VXXXXA, VXXXXB, VXXXXC, XX0044 $$
, WT1###, 1.38E5, 100., 5.E3
$INCLUDE, C:\ATP\project\usp\logic1.lib, WTX###, XX1378, STTSWT, XX0748 $$
, XX0754, XX0760, XX0744, XX0746, XX0750, 100.
$INCLUDE, C:\ATP\project\usp\invcav.lib, XX0816, XX0818, XX0820, XX0858 $$
, ######, XX0748, XX0750, XX0754, XX0744, XX0760, XX0746, 500., 0.5 $$
, 1.
, XX0790, XX0796, XX0780, XX0782, XX0786, 100.
, ######, XX0784, XX0786, XX0790, XX0780, XX0796, XX0782, 500., 0.5 $$
, 1.
$INCLUDE, C:\ATP\project\usp\tranyy.lib, X0017A, X0017B, X0017C, XX0864 $$
, XX0866, XX1194, XX0836, XX0838, XX0840, 1., 1., 0.157, 0.157 $$
,7.4E-4,7.4E-4, 0.0
$INCLUDE, C:\ATP\project\usp\trandy.lib, X0029A, X0029B, X0029C, XX0836 $$
, XX0838, XX0840, XX1210, XX1208, XX1206,1.7321, 1., 0.474, 0.157 $$
,0.0022,7.4E-4, 0.0
, XX1496, XX1498
, XX1492, XX1494, NETRY1
, XX0850, XX0848, CNTINV
, XX0908, XX0914, XX0822, XX0876, XX0896, 100.
, ######, XX0880, XX0896, XX0908, XX0822, XX0914, XX0876, 500., 0.5 $$
, 1.
, XX0944, XX0950, XX0934, XX0936, XX0940, 100.
, ######, XX0938, XX0940, XX0944, XX0934, XX0950, XX0936, 500., 0.5 $$
, 1.
, XX1208, XX1206, XX0988, XX0990, XX0992, 1., 1., 0.157, 0.157 $$
,7.4E-4,7.4E-4, 0.0
, XX0990, XX0992, XX1204, XX1202, XX1200,1.7321, 1., 0.474, 0.157 $$
,0.0022,7.4E-4, 0.0
, XX1042, XX1044
, XX1038, XX1040, NETRY3
, XX1056, XX1062, XX1046, XX1048, XX1052, 100.
, ######, XX1050, XX1052, XX1056, XX1046, XX1062, XX1048, 500., 0.5 $$
, 1.
, XX1094, XX1100, XX1084, XX1086, XX1090, 100.
, ######, XX1088, XX1090, XX1094, XX1084, XX1100, XX1086, 500., 0.5 $$
, 1.
, XX1202, XX1200, XX1138, XX1140, XX1142, 1., 1., 0.157, 0.157 $$
,7.4E-4,7.4E-4, 0.0
, XX1140, XX1142, XX1404, XX1402, XX1400,1.7321, 1., 0.474, 0.157 $$
,0.0022,7.4E-4, 0.0
, XX1196, XX1198
, XX1190, XX1192, NETRY5
, XX1224, XX1232, XX1212, XX1214, XX1218, 100.
$INCLUDE, C:\ATP\project\usp\ruptor.lib, X0013A, X0013B, X0013C, XX1452 $$
, XX1454, XX1456, XX1368
, ######, XX1216, XX1218, XX1224, XX1212, XX1232, XX1214, 500., 0.5 $$
, 1.
, XX1262, XX1268, XX1252, XX1254, XX1258, 100.
, ######, XX1256, XX1258, XX1262, XX1252, XX1268, XX1254, 500., 0.5 $$
, 1.
$INCLUDE, C:\ATP\project\usp\medida.lib, X0009A, X0009B, X0009C, X0013A $$
, X0013B, X0013C, VXXXXA, VXXXXB, VXXXXC, IXXXXA, IXXXXB, IXXXXC
, XX1402, XX1400, XX1312, XX1314, XX1316, 1., 1., 0.157, 0.157 $$
,7.4E-4,7.4E-4, 0.0
$INCLUDE, C:\ATP\project\usp\trandy.lib, XX1300, XX1308, XX1310, XX1312 $$
, XX1314, XX1316, XX1338, XX1338, XX1338,1.7321, 1., 0.474, 0.157 $$
,0.0022,7.4E-4, 0.0
, XX1396, XX1398
, XX1392, XX1394, NETRY7
$INCLUDE, C:\ATP\project\usp\ruptor.lib, XX0401, XX0403, XX0405, X0009A $$
, X0009B, X0009C, XX1368
$INCLUDE, C:\ATP\project\usp\gndelt.lib, X0037A, X0037B, X0037C, XX0864 $$
, XX0866, XX1194
$INCLUDE, C:\ATP\project\usp\trfaco.lib, XX1452, XX1454, XX1456, X0037A $$
, X0037B, X0037C, 4.,1.7321, 12.63, 0.789,0.0952, 0.006, 0.0
$INCLUDE, C:\ATP\project\usp\medida.lib, XX1490, XX1492, XX1494, X0017A $$
, X0017B, X0017C, VA61##, VB61##, VC61##, XX0602, XX0603, XX0604
$INCLUDE, C:\ATP\project\usp\medida.lib, XX1388, XX1396, XX1398, XX1300 $$
, XX1308, XX1310, VA68##, VB68##, VC68##, XX0686, XX0687, XX0688
$INCLUDE, C:\ATP\project\usp\pilimi.lib, ERROQ1, DELTAX, XX1540, XX1538 $$
, 1., 50.
/SWITCH
/SOURCE
60POSFSA 0 -1. 10.
60POSFSB 0 -1. 10.
60POSFSC 0 -1. 10.
60XX0473 0 -1. 10.
60XX0477 0 -1. 10.
60XX0482 0 -1. 10.
60XX0483 0 -1. 10.
/INITIAL
/OUTPUT
BLANK TACS
BLANK MODELS
BLANK BRANCH
BLANK SWITCH
BLANK SOURCE
BLANK INITIAL
BLANK OUTPUT
BLANK PLOT
BEGIN NEW DATA CASE
BLANK

PDF

Cargado por

Copyright:

Formatos disponibles

PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

PDF

Cargado por

Copyright:

Formatos disponibles

COMPENSADOR ESTTICO SNCRONO ( STATCOM ):

V Fasor de tensin de barra

I Fasor de corriente de barra

3 entre primario y secundario

3 veces las tensiones de fase.

), en estado estable la con-

V(t), mediante una suma vectorial en cada instante de tiempo, con-

Figure 4.8: Descomposicin del vector intermedio

Figura 5.4: Esquema general del sistema de control del STATCOM

Figura 6.2: Tcnica de control PWM para control de tensin

Figura 7.2: Sistema de control para el caso de simulacin 1

Figura 7.13: Sistema de control de fase

3) respectivamente. Con una reactancia de 5 % y una resistencia de

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.