Facts HVDC
Facts HVDC
Facts HVDC
1
Escuela Politécnica Nacional (EPN), Quito, Ecuador
E-mail: mgonzalezc.ie@gmail.com;
franklin.quilumba@epn.edu.ec
impulsó el desarrollo de tecnologías capaces actualmente de
Resumen— En este artículo, se presenta una síntesis de la trabajar incluso con decenas de megavatios.
aplicación y ventajas de los dispositivos FACTS y de la transmisión La aplicación particular de la tecnología de la electrónica de
HVDC en los sistemas eléctricos de potencia. Posteriormente se potencia en sistemas de potencia, en conjunto con modernos
enfoca en el estudio del principio de funcionamiento y la sistemas de control, son los denominados Sistemas Flexibles de
modelación matemática para estado estable de los dispositivos
FACTS: SVC, STATCOM y TCSC y además enlaces HVDC de
Transmisión en Corriente Alterna: FACTS (Flexible AC
dos terminales basados en conversores tipo fuente de voltaje (VSC) Transmission Systems por sus siglas en inglés). Los dispositivos
específicamente para el cálculo de flujos de potencia en redes que FACTS se han convertido en las últimas décadas, en
contienen estos dispositivos. Para esto se identifica y presenta las herramientas indispensables a considerar en la planificación y
variables características de los diferentes tipos de conversores expansión de las redes eléctricas, no solo para reducir o incluso
estáticos que constituyen a estos dispositivos y que intervienen en eliminar sus limitantes, sino además para facilitar su operación
el algoritmo de Newton-Raphson para el cálculo de flujos de
potencia.
y control, y mejorar su confiablidad. Todo esto gracias a la
flexibilidad y rapidez de actuación de los sistemas de control y
Índices—Controladores FACTS, Electrónica de Potencia, Flujo a las ventajas que ofrecen los dispositivos electrónicos en
de Potencia, Transmisión HVDC. cuanto a controlabilidad.
109
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
Fig. 2. Efecto del control de las variables de la ecuación (1) para diferentes
valores operativos de las variables ܸ, ߠ y ݔ் expresadas como ܸԢ, ߠԢ y ݔ் Ԣ y
Fig. 1. Sistema de dos barras enlazadas a través de una línea de transmisión de ܸԢԢ, ߠԢԢ y ݔ் ԢԢ respectivamente.
impedancia xLT
IV. TIPOS DE CONTROLADORES FACTS
Este sistema representa la interconexión entre dos
subestaciones a través de una línea de transmisión corta de Los controladores FACTS se agrupan en 3 categorías [1]
impedancia ݔ் . El flujo de potencia desde la barra 1 hacia a la según el principio de compensación en la red y su forma de
barra 2, está determinado por la ecuación: conexión. Estas categorías son: controladores serie,
controladores paralelo y controladores combinados mostrados
ܸଵ ܸଶ (1)
en la Fig. 3.
ܲଵଶ ൌ ሺߠଵ െ ߠଶ ሻ
ݔ்
110
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
111
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
inglés es un controlador FACTS basado en tiristores que está ͻͲι ߙ ͳͺͲι (4)
conformado por una combinación de varias ramas: capacitores
fijos accionados mecánicamente (MSC o FC), reactores A continuación en la Fig. 8, se presentan las formas de onda
controlados por tiristores (TCR) y capacitores accionados por para dos ángulos de disparo de los tiristores del TCR.
tiristores (TSC); todas conectadas en paralelo con la red
eléctrica [3], como se muestra en la Fig. 6. El modelo en el que
se centra este documento es un SVC constituido por una rama
TCR y capacitores fijos.
112
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
al dimensionamiento de la capacidad de las impedancias análisis de Fourier, se obtiene que la impedancia a frecuencia
inductiva y capacitiva del SVC. fundamental del TCSC es igual a:
Fig. 11. Circuito Equivalente de un TCSC. Los tiristores son modelados como
un interruptor normal.
113
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
la ecuación (7) presenta polos, los cuales están determinados resonancia, tanto para la zona capacitiva como para la zona
por la siguiente expresión: inductiva.
A partir de estas consideraciones, se puede construir una
ሺʹ݊ െ ͳሻξ߱ߨܥܮ (13) curva característica V-I del TCSC de la misma manera que un
ߙ௦ ൌ ߨ െ ǡ݊ ൌ ͳǡʹǡ͵ ǥ SVC como se puede observar en la Fig. 14.
ʹ
5
VTCSC
x 10
2
kV
-2
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
5 VTIRISTORES
x 10
2
kV
0
-2
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
4 ILinea
x 10
1
kA
-1
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
4 IL
x 10
2
Fig. 14. Característica V-I (voltaje de compensación a corriente de línea) de
kA
0
-2
un TCSC en control de reactancia [4]
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
4 IC
x 10 La máxima compensación a la reactancia de la línea será
2
cuando la reactancia capacitiva equivalente alcance su valor
kA
0
-2
0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3
máximo, lo cual se obtiene en ߙ . Cabe destacar que para
Fig. 12. Formas de onda en un TCSC para un ángulo de disparo de 150 estado estable, trabajar en la zona inductiva no trae ningún
grados. beneficio en particular, sino todo lo contrario, debido a que la
reactancia del TCSC incrementaría la distancia eléctrica de la
línea, influyendo directamente sobre la estabilidad del sistema
y aumentando las pérdidas.
ZONA INDUCTIVA
C. STATCOM
D
El compensador estático sincrónico o STATCOM por sus
siglas en inglés (Static Compensator), mostrado
esquemáticamente en la Fig. 15, es un controlador FACTS de
conexión en derivación, cuyo conversor es un VSC controlado
D mediante PWM, ya sea basado en GTO o IGBT. Al ser un
D
dispositivo de conexión en paralelo a la red, el objetivo del
D
STATCOM es brindar compensación reactiva a la red mediante
ZONA CAPACITIVA
Fig. 13. Curva de la impedancia del TCSC frente al ángulo de disparo de los
tiristores.
114
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
potencia activa para compensar las pérdidas de conmutación y su máxima potencia reactiva capacitiva incluso en un voltaje
así mantener el voltaje DC constante [7]. tan bajo como 0.15 p.u. [7]. Esta característica le otorga una
ventaja muy amplia frente a dispositivos que tienen una
finalidad similar como el SVC, cuyo aporte de potencia reactiva
depende del voltaje de la barra. Gracias a esta característica, el
voltaje de conexión no es un limitante de este dispositivo, sobre
todo durante fallas o fenómenos transitorios.
115
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
A. Algoritmo para la Incorporación del SVC Los elementos del Jacobiano que deben recalcularse
Al insertar un SVC en la red, la información de la corresponden a aquellos que están relacionados con las dos
susceptancia equivalente de éste, no está incluida inicialmente barras de conexión:
en la matriz de admitancia de barra. Si ݇ es la barra de conexión
߲ܲ ߲ܲ ߲ܲ ߲ܲ ߲ܲ
del SVC, los elementos del Jacobiano relacionados ۍ ܸ ܸ ېοߠ ې
directamente con la inserción de este dispositivo y por ende que ߠ߲ ێ ߲ߠ ߲ܸ ߲ܸ ߲ߙ ۍ ۑ
ێ ۑ
necesitan ser modificados se pueden observar en la ecuación ߲ܲ ێ ߲ܲ ߲ܲ ߲ܲ ߲ܲ ۑ
ێοߠ ۑ
ߠ߲ ێ ܸ ܸ
߲ߠ ߲ܸ ߲ܸ ߲ߙ ێ ۑ ۑ
(18). El valor de la susceptancia del SVC ܤௌ debe ser añadida ߲ܳ ێ
߲ܳ ߲ܳ ߲ܳ ߲ܳ ێ ۑοܸ ۑ
al respectivo valor de ܤ para el cálculo de todas las derivadas. ێ ܸ ܸ ۑ
ߠ߲ ێ ߲ߠ ߲ܸ ߲ܸ ߲ߙ ܸ ێ ۑ ۑ
߲ܳ ێ ێ ۑ
߲ܳ ߲ܳ ߲ܳ ߲ܳ ۑοܸ
οܲ ߲ܲ ߲ܲ ߲ܲ ߠ߲ ێ ܸ ܸ ۑ ێۑ
ߠ߲ ۍ ې ܸ ېοߠ ې ߲ߠ ߲ܸ ߲ܸ ߲ߙ ۑ ܸ ێ
ۍ ߲ߙௌ ۍ ۑ ێ ۑ
߲ܸ
ێ ߲ܳ ێ ۑ ߲ܳ ߲ܳ ێ ۑοܸ ۑ
ێ߲ܲ ்ௌ ߲ܲ ்ௌ
ܸ
்߲ܲௌ
ܸ
்߲ܲௌ ்߲ܲௌ ێ ۑ ۑ
(26
ێοܳ ۑൌ ێێ
ܸ ێ ۑ (18) ߠ߲ ۏ ߲ߠ ߲ܸ ߲ܸ ߲ߙ ۏ ےοߙ ے
)
ێ ێ ۑ ߲ߠ ߲ܸ ߲ߙௌ ܸ ێ ۑ ۑ οܲ
ۑ ۍ ې
ێ ߲ܳێ ۑௌ ߲ܳௌ ߲ܳௌ ێ ۑ ۑ
ێ ۑ
ۏοܳௌ ߙ߲ ۏ ے ܸ ۏοߙௌ ے
ௌ ߲ܸ ߲ߙ ے
ௌ ێοܲ ۑ
ێ ۑ
Donde: ێ ۑ
ൌ ێοܳ ۑ
ܸଶ ܺ (19) ێ ۑ
ܳௌ ൌ െ ܺ െ ሺʹሺߨ െ ߙሻ ʹߙሻ൨
ܺ ܺ ߨ ێοܳ ۑ
ێ ۑ
ێ ۑ
Nótese que se añade el ángulo de disparo de los tiristores ۏο்ܲௌ ے
como variable de cálculo en el algoritmo. Las derivadas
correspondientes a esta variable son: Donde:
߲ܳௌ (20) (27)
ܸ ൌ ܳௌ െ ܸଶ ܤௌ ൌ ʹܳௌ ்ܲௌ ൌ ܲ ൌ െܸ ܸ ܤ ߠ ݊݁ݏ
߲ܸ
De igual manera que con el SVC, se añade al ángulo de
߲ܳ ߲ܳௌ ʹܸଶ (21) disparo de los tiristores del TCSC al algoritmo de cálculo del
ൌ ൌ ሾ
ሺʹߙሻ െ ͳሿ flujo de potencia. Las derivadas correspondientes a esta
߲ߙௌ ߲ߙௌ ߨܺ
variable son:
Cuando se controla el voltaje de la barra de conexión a un ்߲ܺௌ
ଶ ሺߪሻ
ൌ െʹܥଵ ሾͳ
ʹߪሿ ܥଶ ݇ ଶ
valor específico ܸ௦௧ , el vector de variables de estado se reduce ߲ߙ
ଶ ݇ߪ (28)
a: െ ܥଶ ݇ ʹߪ ݇ߪ
(22) ܥଶ ሾ ʹߪ ߪ െ ͳሿ
ࢄ ൌ ሾοߠ οߙ ሿ்
ௌ
116
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
barra de conexión ݇ del STATCOM. Los elementos que La dimensión del Jacobiano de un sistema HVDC de dos
relacionan al STATCOM con la red son: terminales con respecto a la de un STATCOM simple se duplica
ya que son añadidas las ecuaciones tanto del rectificador como
߲ܲ ߲ܲ ߲ܲ ߲ܲ ߲ܲ del inversor. El vector de variables para este caso se muestra a
ۍ ܸ௩ோ ݉ᇱ ېοߠ௩ோ ې
ێ ߲ߠ௩ೃ ߲ܸ௩ோ ߲߶ ߲݉ᇱ ߲ܤ ۍ ۑ continuación:
߲ܳ ێ ߲ܳ ߲ܳ ߲ܳ ߲ܳ ێ ۑοܸ௩ோ ۑ οߠ௩ೃ
ܸ௩ோ ݉ᇱ ێ ۑ ۍ ې
ߠ߲ ێ௩ ߲ܸ௩ோ ߲߶ ߲݉ᇱ ߲ܤ ܸ ێ ۑ௩ோ ۑ
ೃ
߲ܲ ێ ߲ܲ ߲ܲ ߲ܲ ߲ܲ ێ ۑ
ۑ ێοߠ ۑ
ێ ܸ௩ோ ݉ᇱ ߶ ۑ ێ௩ ۑ
ߠ߲ ێ௩ೃ ߲ܸ௩ோ ߲߶ ߲݉ᇱ ߲ܤ ێ ۑ ۑ
ێ ۑ ێοܸ௩ோ ۑ
߲ܳێ௩ோି ߲ܳ௩ோି ߲ܳ௩ோି ߲ܳ௩ோି ߲ܳ௩ோି ۑο݉ᇱ ܸ ێ௩ோ ۑ
ߠ߲ ێ ܸ௩ோ ݉ᇱ ۑ ێ ۑ
௩ೃ ߲ܸ௩ோ ߲߶ ߲݉ᇱ ߲ܤ ݉ ێ ۑᇱ
ێοܸ ۑ
ێ ۑ ௩ூ
߲ܳ ێ ܸ௩ோ
߲ܳ ߲ܳ
݉ᇱ
߲ܳ ߲ܳ ێ ۑ ۑ ێ
ܸ
ۑ
ߠ߲ ۏ௩ೃ ߲ܸ௩ோ ߲߶ ߲݉ᇱ ߲ܤ ۏ ے ܤ ے
(37) ێ ௩ூ ۑ
οܲ ێο߶௩ோ ۑ
ۍ ې
ێ ۑ ࢄ ൌ ێο݉ᇱ ۑ
ێோ ۑ
ێοܳ ۑ ᇱ
ێ ۑ ݉ ێோ ۑ
ێ ۑ ێοۑ ܤ
ൌێ οܲ ۑ ێೃ ۑ
ێ ۑ ێ ۑ
ێοܳ
௩ோି
ۑ ێο߶௩ூ ۑ
ێ ۑ ێο݉ூᇱ ۑ
ێ ۑ ݉ ێᇱ ۑ
ۏοܳ ے ێூ ۑ
ۏοܤ ے
Debido a que no existe transferencia de potencia reactiva
hacia la barra 0, además que la potencia activa total inyectada
en el nodo 0 es cero debido a que no existen ramales, para todas
las condiciones de control se cumple que:
117
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
VIII. APLICACIÓN
TABLA II
Los modelos presentados en la sección anterior, fueron PARÁMETROS DEL STATCOM
implementados en un programa computacional para el cálculo
de flujos de potencia escrito en lenguaje Python. Este programa Parámetro Valor Unidad
se denomina FACTS-PF y permite el cálculo de las variables de Resistencia del Conversor ܴ௩ோ 0.066 p.u,
estado que intervienen en los flujos de potencia en los sistemas Reactancia del Conversor ܺ௩ோ 0.833 p.u.
de potencia, además de las variables de los controladores Modo de Control Voltaje Local -
Consiga de Voltaje ܸ௦௧௧ 1.000 p.u
FACTS que puedan ser agregadas a ellos. Pérdidas de Conmutación 100 kW
Voltaje DC 1.4142 p.u
A. Aplicación del programa FACTS-PF a una red de prueba Potencia Reactiva Máxima (Ind.) 50 MVAr
con FACTS Potencia Reactiva Máxima (Cap.) 100 MVAr
Para validar los resultados del programa FACTS-PF, se
utilizó la red IEEE de 14 barras, a la cual se le aplicaron ciertas
TABLA III
modificaciones para observar el impacto de los controladores PARÁMETROS DELTCSC
FACTS en un sistema de potencia. En primer lugar la carga fue
incrementada en un 40% y posteriormente el compensador Parámetro Valor Unidad
sincrónico de la barra 6 fue sacado de servicio. Reactancia Inductiva Total ܺ 1.0325x10-3 p.u,
A esta red, como se puede ver en la Fig. 19, se agregan 3 Reactancia Capacitiva Total ܺ 1.5752x10-2 p.u.
dispositivos FACTS: Un SVC en la barra 14, un STATCOM en Modo de Control Flujo Controlado -
Consiga de P ்ܲௌ 100 MW
la barra 4 y un TCSC en serie con la línea 4-5 del lado de la Ángulo de Resonancia 145° -
barra 5. Los datos de líneas, transformadores, barras y cargas Rango de Operación 90° - 180° -
fueron tomados de [12].
A continuación se presentan los resultados de voltaje y
potencia generada en las barras de generación, obtenidas
cálculo de flujos de potencias por el programa FACTS-PF para
el sistema de IEEE de 14 barras con FACTS.
TABLA IV
RESULTADOS DEL FLUJO DE POTENCIA
118
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
119
XXVI JORNADAS EN INGENIERÍA ELÉCTRICA Y ELECTRÓNICA - EPN
[11] E. Acha, B. Kazemtabrizi, y L. M. Castro, «A New VSC-HVDC Model Franklin L. Quilumba (S’10-M’14)
for Power Flows Using the Newton-Raphson Method», IEEE Trans.
Power Syst., vol. 28, n.o 3, pp. 2602-2612, ago. 2013.
obtuvo el título de Ingeniero Eléctrico en
[12] «Power Systems Test Case Archive», University of Washington la Escuela Politécnica Nacional en Quito,
Electrical Engineering. [En línea]. Disponible en: Ecuador, en el 2008. Realizó sus estudios
https://www.ee.washington.edu/research/pstca/. [Accedido: 22-feb-
de posgrado en la Universidad de Texas
2016].
Arlington, en Arlington, Estados Unidos
de América, donde obtuvo el grado de
Master of Science in Electrical
Mario A. González nació en la ciudad Engineering y el título de Doctor of
de Quito, Ecuador en 1991. Finalizó sus Philosophy Ph.D. en el 2014.
estudios en la Carrera de Ingeniería Desde el 2014, el Dr. Quilumba es parte del cuerpo docente
Eléctrica en la Escuela Politécnica de la Escuela Politécnica Nacional, donde es Profesor
Nacional en Quito en 2015. Agregado. Sus áreas de interés son análisis, operación,
Actualmente se encuentra estabilidad y control de sistemas eléctricos de potencia,
concluyendo su proyecto de titulación modelación de carga, planeación de generación y transmisión;
para obtener el título de Ingeniero respuesta de la demanda; y predicción de carga.
Eléctrico en la misma universidad. Sus
áreas de interés son análisis, operación, y control de sistemas
eléctricos de potencia, y controladores FACTS.
120