TFG Jesus Lopez Merino 2016
TFG Jesus Lopez Merino 2016
TFG Jesus Lopez Merino 2016
Estudio y Simulación de un
I
Abstract
The aim of this project is the analysis of the so called, STATic synchronous COMpensator
(STATCOM), which is one of the most important power electronics applications in the
field of flexible alternating current transmission system (FACTS).
In recent times, several researches about STATCOM have been developed resulting in
numerous publications in journals and technical reports. Therefore to be able to acquire
a deep knowledge on this subject is an arduous and complicated task. This work presents
a useful collection of the most important pieces of information found about this topic, to
be used by future students and researchers interested in the subject.
The control of the system for the reactive power compensation has been implemented
as a cascade tuned PI-controller, which offers the DC bus voltage control, power factor
correction and the AC voltage control in an efficient manner.
II
Índice general
Página
Resumen I
Abstract II
Índice de Figuras V
Índice de Tablas IX
Lista de Acrónimos X
1 Introducción 1
1.1 Situación actual . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Problemática . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 Objetivo del trabajo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4 Estructura del documento . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
III
Índice general
3.6.2 Lazo de potencia reactiva . . . . . . . . . . . . . . . . . . . . . . . 50
3.6.3 Lazo de tensión eficaz del PCC . . . . . . . . . . . . . . . . . . . . 53
3.6.4 Lazo de tensión en la etapa de continua . . . . . . . . . . . . . . . . 55
3.7 Anti-windup para reguladores PI . . . . . . . . . . . . . . . . . . . . . . . 59
3.8 Implementación del SPWM . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4 Resultados de la Simulación 68
4.1 STATCOM como regulador de tensión . . . . . . . . . . . . . . . . . . . . 70
4.1.1 Primer caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.1.2 Segundo caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
4.1.3 Tercer caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.1.4 Cuarto caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
4.2 STATCOM como corrector del factor de potencia . . . . . . . . . . . . . . 83
4.2.1 Primer caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
4.2.2 Segundo caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.2.3 Tercer caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
4.2.4 Cuarto caso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
5 Conclusión 97
5.1 Objetivos cumplidos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.2 Posibles mejoras del modelo presentado . . . . . . . . . . . . . . . . . . . . 97
5.3 Trabajos futuros . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
6 Presupuesto 100
6.1 Presupuesto personal. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
6.2 Análisis económico de la inversión . . . . . . . . . . . . . . . . . . . . . . . 103
Bibliografı́a 111
Apéndices 112
V
Índice de figuras
Tabla 2.1 Modo de actuación de los FACTS, ante problemas del sistema eléctrico 9
Tabla 4.1 Caso 1, como controlador de la tensión eficaz del punto de conexión
a la red . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
Tabla 4.2 Caso 2, como controlador de la tensión eficaz del punto de conexión
a la red . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Tabla 4.3 Caso 3, como controlador de la tensión eficaz del punto de conexión
a la red . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
Tabla 4.4 Caso 4, como controlador de la tensión eficaz del punto de conexión
a la red . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
Tabla 4.5 Caso 1, como corrector del factor de potencia . . . . . . . . . . . . . 83
Tabla 4.6 Caso 2, como corrector del factor de potencia . . . . . . . . . . . . . 88
Tabla 4.7 Caso 3, como corrector del factor de potencia . . . . . . . . . . . . . 91
Tabla 4.8 Caso 4, como corrector del factor de potencia . . . . . . . . . . . . . 94
IX
Lista de Acrónimos
X
Lista de Acrónimos
Introducción
La industria del suministro eléctrico está sufriendo una profunda transformación mundial
debido a que el mercado se orienta hacia fuentes de energı́a renovables y la demanda
energética está en constante aumento. En dicho contexto, la expansión de las instalaciones
está siendo frustrada en muchas ocasiones por las exigencias medioambientales, el uso del
suelo y las presiones reguladoras, todas estas barreras evitan la concesión de licencias de
construcción de nuevas lı́neas de transmisión y plantas de generación. Una de las opciones
disponibles para gestionar y maximizar los recursos de transmisión ya existentes, con
elevados niveles de seguridad y estabilidad, es la utilización de los sistemas flexibles de
transmisión de alterna (FACTS) [1].
En los últimos tiempos los dispositivos FACTS son ampliamente conocidos, por su elevada
capacidad de control del sistema eléctrico. La significativa rapidez de respuesta de los
FACTS es uno de sus progresos más relevantes, la cual permite una mejora notable en
las ondas de tensión e intensidad, reaccionando en periodos de tiempo muy cortos con
respuestas idóneas en estado estacionario o incluso en condiciones de transitorio. Además
de la ventajosa dinámica de estos dispositivos, si se analiza factores estáticos tales como
las pérdidas mecánicas, también existen beneficios [2].
Uno de los dispositivos FACTS de nueva generación que emplea algunos de los desarrollos
electrónicos más recientes es el compensador estático sı́ncrono de potencia reactiva
(STATCOM), representado en la Figura 1.1. Dicho dispositivo presenta un funcionamiento
más complejo y un control más amplio, con la finalidad de actuar como un compensador
de potencia reactiva de una serie de cargas eléctricas. El resultado generado por fenómenos
circunstanciales que provocan variaciones de la tensión, las cuales llegan a ser muy
pronunciadas e inadmisibles, es un problema que un compensador de potencia puede
resolver eficazmente [2].
El objeto de estudio de este proyecto es la configuración y modo de operación de dicho
dispositivo (STATCOM) en los sistemas de distribución de media tensión. El análisis
general se focaliza en el control de dos de sus funciones principales [3]:
Regular la tensión del nudo de la red al que se conecta (PCC).
Corregir el factor de potencia del PCC.
1
1.1. Situación actual
STATCOM 20 3000
TCSC 10 2000
Tabla 1.1. Estimación del número de dispositivos FACTS en el mundo y su potencia total
instalada [2]
1.2. Problemática
Para lograr un funcionamiento idóneo del sistema eléctrico se debe lograr evitar la
presencia de armónicos, que impiden que se cumpla las condiciones necesarias para el
funcionamiento correcto del sistema, tales como mantener la tensión y la frecuencia
fundamental en valores constantes. Otro aspecto relevante es el factor de potencia
que será un ı́ndice de la cantidad de potencia activa y en definitiva, potencia útil
que se está logrando aprovechar. Dichas exigencias deben lograrse para cualquier carga
representativa del consumidor y ante cualquier perturbación que provoque variaciones en
la tensión o la intensidad del sistema [7].
la reducción de la vida útil de los elementos que forman parte de dicho sistema. En
definitiva, se desea un factor de potencia unitario para aumentar la potencia real o útil
y evitar la inversión en dispositivos capaces de soportar intensidades superiores. Estas
circunstancias ponen en valor la capacidad del STATCOM de actuar como corrector del
factor de potencia [3].
Otra alternativa es la de regular la tensión del nudo al que se conecta, para mantener
un valor constante frente a las perturbaciones producidas en el sistema. La complejidad
de ofrecer una solución adecuada se eleva en circunstancias donde las variaciones que se
deben hacer frente son enormes y rápidas. Un ejemplo idóneo de condiciones de carga
con constantes cambios de extremada rapidez es el horno de arco eléctrico, donde el
inicio de fundición es muy irregular y desequilibrado en un pequeño periodo de tiempo,
lo que exige una respuesta rápida y eficaz del compensador dinámico. Aunque no existe
una fórmula empı́rica representativa del concepto denominado como calidad de energı́a,
se encuentran parámetros reveladores, tales como la distorsión armónica total (THD) o
la variación máxima de la tensión eficaz del suministro en un determinado periodo de
tiempo. Cabe destacar que la tensión es un elemento muy sensible a las fluctuaciones,
al que el STATCOM ofrece una solución. Para mantener la seguridad y estabilidad del
sistema eléctrico, las instalaciones de suministro deben situarse en unos lı́mites del orden
de ± 5 % de la tensión nominal, medido sobre un determinado periodo de tiempo (de
pocos minutos a horas) [3].
Este problema tradicionalmente habı́a sido solucionado por las compañı́as mediante la
interconexión de diferentes sistemas eléctricos. Sin embargo, esto implica hacer más
tendidos de lı́neas de transportes que encarecen el presupuesto. Dispositivos como el
STATCOM son una solución más flexible, que ofrece un resultado satisfactorio tanto
de manera local como en un área más amplia del sistema eléctrico.
7
2.1. Sistemas flexibles de transmisión de corriente alterna, FACTS
Figura 2.1. SVC para estabilizar un sistema de transmisión de 735 kV en Chénier, Canada [9]
Lı́mites de tensión:
Lı́mites térmicos:
Flujos circulares:
Tabla 2.1. Modo de actuación de los FACTS, ante problemas del sistema eléctrico [1]
Figura 2.2. Dos TCSC instalados por PGCIL (Power Grid Corporation of India Ltd) [10]
(a) Tensión e intensidad caracterı́stica del (b) Tensión e intensidad caracterı́stica del
STATCOM SVC
Otros factores beneficiosos que suponen un añadido a todas las ventajas operativas
comentadas del STATCOM, son los bajos costes de operación y mantenimiento. Todas
estas ventajas se fundamentan en dispositivos que a diferencia de los tiristores ofrecen
la posibilidad de apagado, aumentando considerablemente su margen de libertad de
actuación [3].
En la actualidad, el empleo de inversores en fuente de tensión para la interconexión en el
sector de la distribución energética es una práctica ampliamente extendida. El próximo
paso para la evolución del STATCOM es la combinación con fuentes de energı́a en la
etapa de continua (DC), lo que posibilitarı́a la combinación de potencia activa y reactiva,
logrando un equilibrado de operación de la red de nivel superior [2].
Figura 2.5. Diagrama fasorial con una carga puramente resistiva [3]
Figura 2.6. Diagrama fasorial con una carga puramente reactiva [3]
Figura 2.7. Diagrama fasorial con una carga puramente capacitiva [3]
Tanto las capacitancias como las reactancias solo producen o absorben potencia reactiva
y solo suponen las pérdidas propias del devanado. Como se observa en la Figura 2.8, la
tensión sufre variaciones en función de la naturaleza de la carga como se ha analizado.
Para mantener la tensión del nudo al que se conecta en el valor deseado, se requiere una
potencia reactiva determinada (QvR ), que suministra el STATCOM como compensador
de reactiva. Dicha potencia provoca variaciones en el punto donde se conecta la carga
(Qs = QvR + Qc ), además es dependiente de la variación del propio fasor ∆V , que
representa la caı́da de la tensión en la impedancia interna expresada en (2.1) [3].
P
∆V = (Rs + jXs ) (2.1)
V
2
Xs P − Rs Qs 2
2
Rs P + Xs Qs
|E| = V +
+
(2.2)
V V
Figura 2.9. Diagrama fasorial de las intensidades, en el modo de corrección del factor de
potencia [3]
Vk VvR
P = sin δvR (2.3)
Xl
Vk2 Vk VvR
Q= − cos δvR (2.4)
Xl Xl
Dónde:
Ik : Vector de la intensidad.
Vk : Vector de la tensión en el nudo de acoplamiento a la red.
VvR : Vector de la tensión de salida del convertidor.
YvR : Matriz de admitancia.
Los parámetros del control tales como la magnitud y el carácter de la potencia reactiva y
activa necesaria para la meta fijada, son establecidos a raı́z de un flujo de carga. Existe
una amplia aceptación de la validez y seguridad que proporciona métodos iterativos como
el Newton-Raphson, a la hora de ofrecer una solución al flujo de carga existente en las
lı́neas de transmisión de la red eléctrica en régimen permanente.
El intercambio de potencia de reactiva y activa entre ambos puntos queda definido en las
siguientes expresiones [1]:
∗ ∗ ∗
SvR = VvR IvR = VvR YvR (VvR − Vk∗ ) (2.7)
2
PvR = VvR GvR + VvR Vk [GvR cos(δvR − θk ) + BvR sin(δvR − θk )] (2.8)
2
QvR = −VvR BvR + VvR Vk [GvR sin(δvR − θk ) − BvR cos(δvR − θk )] (2.9)
PvR = Vk2 GvR + Vk VvR [GvR cos(θk − δvR ) + BvR sin(θk − δvR )] (2.10)
QvR = −Vk2 BvR + Vk VvR [GvR sin(θk − δvR ) − BvR cos(θk − δvR )] (2.11)
Cabe destacar que en este caso se ha optado por convertidores convencionales de fuente
en tensión de dos niveles en lugar de los multinivel. Sin embargo, estos últimos presentan
ventajas, especialmente en aplicaciones de media y gran potencia, dado que permiten
operar en un mayor número de niveles de la tensión de DC, logrando una representación
mas fidedigna de la onda sinusoidal de la tensión de salida del convertidor [15].
Los convertidores además de clasificarse por su configuración, pueden distribuirse en dos
categorı́as en función de la fuente de entrada utilizada en el lado de continua (DC) [3]:
1. Los inversores en fuente de tensión (VSC) en el lado de continua, generalmente un
condensador. Esta opción permite un flujo bidireccional, dependiendo del signo de
la corriente que puede ser negativo o positivo.
2. Los inversores en fuente de corriente (CSCs). Esta opción presenta como ventajas
principales que las inductancias que se emplean son más económicas que los
condensadores, porque presentarı́a menos pérdidas de Joule en el sistema de
almacenamiento, y además el exceso de tensión tienen mayor probabilidad de
ser perjudicial para el condensador que el exceso de las intensidades para las
inductancias [16]. No obstante, los STATCOM basados en el empleo de los CSC
no están desarrollados y no se comercializan, principalmente por sus problemas
potenciales. Dichos defectos pueden llegar a ser subsanados por nuevas tecnologı́as
desarrolladas en la superconducción a elevadas temperaturas o temperaturas medias
(como el recientemente descubrimiento del boruro de magnesio con una temperatura
de transición de 39 K [16] [17]).
Aunque la opción de emplear una fuente de corriente serı́a más idónea por la elevada
eficiencia de sus dispositivos, en general en aplicaciones de compensación de reactiva, se
hace uso de una fuente de tensión, como es el caso del STATCOM.
Figura 2.15. Clase de semiconductor totalmente controlado idóneo, en función del nivel de
potencia y frecuencia del convertidor [3]
50 50 50
0 0 0
1
Señal portadora
Tensión (p.u.)
0.5 Señal referencia
−0.5
−1
100 200 300 400 500 600
Tensión (p.u.)
−1
f portadora
mf =
f moduladora
Vm,moduladora
ma =
Vf,portadora
Para este tipo de modulación mientras su valor sea inferior a la unidad (ma ≤ 1),
la modulación es lineal y el valor pico de la moduladora no sobrepasa al de la
portadora, siendo el modo de funcionamiento correcto. Pero en el caso de que dicho
ı́ndice supere la unidad (ma > 1), el valor pico de la señal moduladora superarı́a al
de la portadora, produciéndose el fenómeno denominado sobremodulación. Como se
observa en la Figura 2.18, en la sobremodulación se produce la pérdida de pulsos y
se pierde la relación lineal entre el ı́ndice de modulación y la amplitud de la tensión
del convertidor, imposibilitando el control de la corriente inyectada a la red [19].
1 Señal portadora
Tensión (p.u.)
Señal referencia
0.5
−0.5
−1
Tensión (p.u.)
−1
En los próximos capı́tulos, siempre que se emplee el término del ı́ndice de modulación, se
corresponde al de la amplitud. Dado que es el factor relevante a la hora de determinar la
magnitud de la tensión de salida del convertidor.
Figura 2.19. Circuito equivalente monofásico para la tensión de salida del inversor [26]
Dónde:
VvR : Valor eficaz de la tensión de fase de la salida del convertidor.
δvR : Ángulo de fase de la tensión de salida del convertidor.
Ik : Valor eficaz de la intensidad de salida del convertidor.
γk : Ángulo de fase de la intensidad de salida del convertidor.
Vk : Valor eficaz de la tensión del PCC.
3 · Vk2 2 3 · Vk · VvRmax 2
P 2 + (Q + ) =( ) (2.12)
Xf Xf
pico
2 VvR f ase
ma = (2.13)
VDC
P 2 + Q2 = (3 · Vk · Imax
2
) (2.14)
Los puntos de operación seguros vienen definidos por la intersección de los cı́rculos que
describen dichos lı́mites fı́sicos, como se muestra en la Figura 2.20, donde se establece la
relación entre las tensiones del sistema y las potencias proporcionadas por el convertidor
[13]. Algunas consideraciones relevantes de esta representación son:
La proyección del eje directo representa la potencia reactiva entregada a la red,
cuando actua como inversor y el eje de cuadratura, representa la potencia activa.
Se emplea una constante para lograr una representación del diagrama fasorial del
esquema de la Figura 2.20, donde se observa que el lı́mite establecido de la potencia
está directamente relacionado con el módulo de la caı́da de tensión que se produce
en la reactancia Xf Ik . Dicha constante es:
Xf
k=
3 · Vk
2.5
2
Qmax(P=0)(p.u.)
1.5
0.5
−0.5
−1
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8
VDC(p.u.)
1.2
1.15
1.1
(p.u.)
1.05
max(P=0)
0.95
Q
0.9
0.85
0.8
0.15 0.16 0.17 0.18 0.19 0.2
X (p.u.)
f
34
Figura 3.1. Esquema del caso de estudio [12]
3 3
P = (vds ids + vqs iqs ) = |v| |i| cos φ (3.2)
2 2
3 3
Q = (vqs ids − vds iqs ) = |v| |i| sin φ (3.3)
2 2
Figura 3.3. Representación de los fasores del sistema en el plano α - β y en el plano d-q [29]
En la Figura 3.3 se muestra los fasores del sistema referido a los ejes α - β y a los ejes
d-q . Donde θ es el ángulo del fasor de la tensión de la red con el eje real del sistema
de referencia estacionario. Como se observa en dicha representación, el vector tensión del
PCC (V~K ) coincide con el eje directo (d), mientras que la componente de dicho vector
en el eje de cuadratura (q) es nula [29]. Ası́ pues los vectores de la tensión e intensidad
quedan definidos en los ejes de referencia d-q como:
ia id
ib = [c]−1 iq (3.4)
ic 0
va |v|
vb = [c]−1 0 (3.5)
vc 0
La potencia activa y reactiva se simplificarı́an sobre los ejes de referencia d-q del siguiente
modo:
3
P = |v| id (3.6)
2
3
Q = − |v| iq (3.7)
2
La teorı́a general del lazo de seguimiento de fase (PLL) y de la opción implementada para
las caracterı́sticas de este sistema, el SRF-PLL queda analizado en el Apéndice D. Esta
técnica se apoya en un sistema de referencia sı́ncrono, sobre el cual se proyecta el vector
de la tensión de red.
A raı́z de la medida de la señal del vector tensión del punto de conexión a la red
trifásica (PCC) y el bloque de transformación abc/dq del esquema general de control
del STATCOM que se muestra en la Figura 3.4 (a), se obtiene las componentes en el
sistema de referencia d-q. La posición angular del SRF-PLL se modifica coincidiendo
con el ángulo de fase del vector de tensión del PCC, logrando la sincronización de las
frecuencias de las variables fundamentales. La implementación del bloque PLL mostrado
en la Figura 3.4 (b) consiste en un regulador PI (proporcional-integral), como se observa
en la Figura 3.4 (c) [30].
Figura 3.4. Bloques del esquema general del control del STATCOM
Las variables de estado de dicho circuito serán las intensidades y las tensiones. Las
ecuaciones del modelo matemático que describen el funcionamiento del sistema [26] [29],
quedan expresadas en el sistema de referencia d-q-o y son válidas indistintamente de si el
convertidor opera como inversor o rectificador, cambiando exclusivamente el signo de las
intensidades.
vvR a ia ia vk a
vvR b = Rf ib + Lf d ib + vk b (3.8)
dt
vvR c ic ic vk c
2π 2π
~vkαβ = vk a + vk b e j 3 + vk c e−j 3 (3.9)
d~ikαβ
~vkαβ = ~vvR
αβ
− Rf ~ikαβ − Lf (3.10)
dt
d~ikdq
~vkdq = dq
~vvR ~ dq
− Rf ik − Lf − jLf w~ikdq (3.11)
dt
dikd
d
vkd = vvR − Rf ikd − Lf + jLf wikq (3.12)
dt
dikq
vkq = q
vvR − Rf ikq − Lf − jLf wikd (3.13)
dt
El fasor de la tensión de la red (~vk ) coincide con el eje directo (d), mientras que las
componentes jLw~ikq y −jLw~ikd son los denominados términos de acoplamiento. Dichas
componentes suponen que cualquier cambio de consigna en el eje directo tiene un efecto
indeseado en el eje de cuadratura y viceversa, impidiendo la ventaja del control vectorial
de desacoplar las magnitudes en cada eje. Por esta razón, como se observa en la Figura
3.6 (a) del esquema general de control, son eliminados para lograr el desacople de ambos
ejes deseado.
(a) Términos del acoplamiento en el esquema del (b) Implementación de los terminos de
circuito general acoplamiento
dq dq 3 d d
PDC = PAC = v i (3.14)
2 k k
dVDC
PDC = −VDC IDC = −VDC (CDC ) (3.15)
dt
La implementación en Simulink presentada en la Figura 3.7 (b) se corresponde a la
obtención de la tensión de la etapa de continua mediante la ecuación de la energı́a (3.16).
Z 2
1 dVDC
WDC = − PDC = CDC (3.16)
2 dt
Sin embargo, existe otros parámetros que requieren una definición constante y dificultan
establecer los lı́mites exactos de las magnitudes comentadas, destacando:
1. El ángulo de fase de la tensión de salida del STATCOM (δvR ).
2. El ángulo de fase de la tensión del punto de unión a la red (θk ).
Dichos factores producen complicaciones debido a que la relación θk − δvR , determina
el flujo de energı́a entre el convertidor y el PCC [32]. Generalmente se suele emplear
la aproximación θko = δvRo , donde θko es el ángulo de fase inicial de la tensión del PCC
cuando el STATCOM no entrega potencia activa. Sin embargo, los valores obtenidos serı́an
una aproximación cercana a la realidad, puesto que realmente existe una potencia activa
suministrada por la red para solventar las pérdidas del filtro, del transformador y del
propio convertidor, ası́ que en la práctica es necesario un valor aproximado de δvR .
A raı́z de la expresión (2.13) en el punto 2.3.4 y considerando que se ha establecido
un ı́ndice de modulación inferior a la unidad, para evitar la sobremodulación ante
perturbaciones o cambios de consignas severas, los valores empleados son:
ma = 0.95
VDC = 230 [V]
El máximo valor pico de la tensión de fase del inversor posible (VvRmax ) en dichas
condiciones serı́a 109.25 V, mientras que en su defecto el término mı́nimo se establecerá en
un valor nulo. Es necesario resaltar que la eliminación de los términos de acoplamiento que
se definió en el punto 3.3, ejercen un efecto en los reguladores de intensidad que deberá ser
considerado y analizado adecuadamente para definir los lı́mites de la tensión.
A la hora de establecer los lı́mites operacionales de la intensidad, primero se debe
establecer una relación directa entre sus componentes directa y de cuadratura. En primera
instancia atendiendo a la ecuación (2.14) en el punto 2.3.4.
P 2 + Q2 = (3 · Imax · Vk )2 (3.17)
3 3
( · Id · Vkd )2 + ( · Iq · Vkd )2 = (3 · Imax · Vkd )2 (3.18)
2 2
Por otro lado, empleando la ecuación de los lı́mites de funcionamiento en los ejes P/Q
mostrada en (2.12) de la sección 2.3.4 y considerando de nuevo que la componente de
cuadratura de la tensión del PCC es nula, se obtiene la ecuación (3.20).
2
2 · Vk2d 2 2 · Vkd · VvRmax 2
(Id · Vkd ) + (Iq · Vkd + ) =( ) (3.21)
Xf Xf
Idmin = −Idmax
Iqmin = −Iqmax
LÍMITES ESTABLECIDOS
VvRmax [V ] 109.25
VvRmin [V ] 0
3.5. Reguladores PI
Un regulador es un elemento capaz de proporciona respuestas y soluciones de control
válidas y eficaces frente a perturbaciones o problemas del entorno real [34] [35], en este
caso el sistema eléctrico. Dichos reguladores pueden constar de tres términos [36]:
El término proporcional (P) que otorga fuerza y rapidez en la acción de control,
pasando la señal de error por una ganancia.
El término integral (I) que reduce el error estacionario hasta que es nulo, a través
de compensaciones de baja frecuencia.
El término derivativo (D) que mejora la respuesta transitoria a través de la
compensación de alta frecuencia.
El uso generalizado del regulador PID, el cual contiene los tres elementos, ha incentivado
el número de investigaciones que analizan las ventajas y desventajas de su actuación
[37] [36], ofreciendo soluciones a los problemas relacionados con los términos derivativo e
integral.
Us 1
GP ID (s) = = kp + + Td s (3.22)
Es Ti s
Donde:
U (s): Señal de control.
E (s): Señal del error.
kp : Ganancia proporcional.
Ti : Constante de tiempo de acción integral.
Td : Constante de tiempo de acción derivativa
ki Tn s + 1
GP I = kp + = (3.23)
s Ti s
kp 1
Tn = Ti =
ki ki
1
Tn s + 1 1 k
Gla = GP I · Gint · Gplanta = · · (3.24)
Ti s Tlcint s + 1 τplanta s + 1
Donde [39]:
Gla (s) – Planta en lazo abierto.
Glc (s) – Planta en lazo cerrado.
GP I (s) - Regulador PI.
Gint (s) – Lazo interno (dinámica rápida).
Gplanta (s) - Sistema a regular (dinámica lenta).
El método empleado para los parámetros del controlador del lazo de intensidad es el de
la cancelación de la mayor constante de tiempo [13] [39]. La respuesta más rápida la
proporciona el actuador y por lo tanto, la constante de tiempo del lazo más interno es el
periodo de conmutación del convertidor (τconvertidor ).
La función de la planta del lazo de intensidad definida en (3.26), se basa en la relación
de la tensión de salida del convertidor y la intensidad, expresada en (3.25), a partir de
las ecuaciones (3.12) y (3.13). Esta es idéntica tanto para la componente directa como de
cuadratura de la corriente.
Id (s) Iq (s) 1
= = (3.25)
Vd (s) Vq (s) Lf s + Rf
1 1
k Rf
Gplanta = = L
(3.26)
1 + τplanta s 1 + Rff s
Lf
TnI s + 1 = s+1 (3.27)
Rf
Lf
TnI = [s] (3.28)
Rf
A raiz de dicha cancelación se halla la expresión del lazo abierto (3.29), a partir de la cual
se define la ecuación de lazo cerrado (3.30).
1
Rf
Gla = (3.29)
TiI s(1 + τconv s)
1 1
Gla Rf ∼ Rf 1
Glc = = 1 = 1 = (3.30)
1 + Gla Rf
+ TiI s · (1 + τconv s) Rf
+ TiI s 1 + Rf TiI s
El hecho de que se opere a bajas frecuencias (en este caso 60 Hz), posibilita la aproximación
de despreciar el término TiI τconv s2 y obtener la constante de tiempo del lazo cerrado de
la intensidad simplificada en la expresión (3.31).
1
Rf 1
Glc = 1 = (3.32)
Rf
+ TiI s · (1 + τconv s) 1 + Rf TiI s + Rf TiI τconv s2
k
Glc = s2
(3.33)
1+ 2 wεo s + wo2
ε
wo = 2 (3.34)
Rf TiI
1
Rf TiI τconv = (3.35)
wo2
2τconv s
TiI = [ ] (3.36)
Rf Ω
Utilizando los términos (3.28) y (3.36), se hallan las ganancias que representan el término
integral y proporcional del regulador de intensidad (3.37) y (3.38).
1 Ω
kiI = [ ] (3.37)
TiI s
Lf [mH] 0.637
Rf [Ω] 0.15
TiI [ Ωs ] 0.0082
kiI [ Ωs ] 121.5
0.8
Referencia
0.6 Control
0.4
iq(p.u)
0.2
−0.2
−0.4
0.95 0.96 0.97 0.98 0.99 1 1.01 1.02 1.03 1.04 1.05
Tiempo (s)
En las Figuras 3.10 y 3.11, se produce una variación significativa de las consignas iniciales,
para representar el efecto que se produce ante perturbaciones o fallas. La respuesta de
los reguladores con los parámetros asignados es idónea, rápida y precisa, validando los
planteamientos teóricos formulados.
1
Referencia
0.8 Control
0.6
id(p.u.)
0.4
0.2
0.95 0.96 0.97 0.98 0.99 1 1.01 1.02 1.03 1.04 1.05
Tiempo (s)
La función del lazo abierto de reactiva (3.41) presenta una estructura semejante al de
corriente.
TnQ s + 1 1 1 −3Ud
GlaQ = GP IQ · Glc · Gplanta = · · · (3.41)
TiQ s 1 + Tlcint s 1 + τconv s 2
−3Ud
Kq =
2
La mayor constante de tiempo en este caso se trata del lazo interno de regulación de la
intensidad, que se cancela con el cero del regulador TnQ , como se expresa en (3.42) [39]
[13]. El lazo abierto de reactiva definido en (3.41), queda simplificado a la expresión (3.43).
Kq
GlaQ = (3.43)
TiQ s · (1 + τconv s)
GlaQ 1
GlcQ = = TiQ TiQ
(3.44)
1 + GlaQ 1+ s+ τ s2
Kq Kq conv
1
GlcQ = TiQ
(3.45)
1+ Kq
s
TiQ
TlcQ = (3.46)
Kq
ε · Kq
wo = 2 (3.47)
TiQ
TiQ 1
τconv = 2 (3.48)
Kq wo
Considerando que: ε = √1
2
1 1
kiQ = [ ] (3.50)
TiQ V s
1
kpQ = kiQ · TnQ [ ] (3.51)
V
El análisis de las unidades mostrado en (3.52), sirve para comprobar la veracidad de los
cálculos estimados en este punto. En esta ocasión la variable de entrada es la potencia
reactiva inyectada o recibida y la variable de salida es la intensidad, de este modo quedarı́a:
Z
1 1
U [A] = kp [ ]E[V A] + kii [ ] Edt[V As] (3.52)
V Vs
La Figura 3.12 muestra el bloque del regulador de potencia reactiva en el control (a) y la
implementación del regulador PI (b), cuyos parámetros se definen en la Tabla 3.3.
Ud [V ] 93.9
TiQ [V s] 0.174
kpQ [ V1 ] 0.0071
A fin de validar el análisis del lazo de potencia reactiva, se somete a dicho regulador a
un cambio brusco en su consigna inicial, para evaluar su respuesta ante perturbaciones o
fallas. Como se muestra en la Figura 3.13, el regulador opera de forma eficiente y eficaz.
Dicho análisis se ofrece en magnitudes por unidad.
2.4
Referencia
2.2
Control
2
1.8
1.6
Qk(p.u.)
1.4
1.2
0.8
0.6
0.4
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Tiempo (s)
El lazo de tensión eficaz del PCC se emplea para mantener dicha tensión en los valores
nominales. El análisis y planteamiento de este lazo es semejante al de reactiva, dado
que ambas funciones son excluyentes y utilizan un regulador PI que responde ante las
variaciones de potencia reactiva producidas por la carga, mediante respuestas controladas
de la componente de cuadratura de la intensidad (iq ).
La ecuación que representa el lazo abierto del valor eficaz de la tensión (3.53), presenta
una estructura muy parecida al de la potencia reactiva (3.41). La relación entre la variable
de entrada (la tensión) y la de salida (la intensidad), atiende a la ley de Ohm y queda
establecida por la propia impedancia del sistema (ZBASE ).
TnVrms s + 1 1 1
GlaVrms = GP IVrms · Glc · Gplanta = · · · ZBASE (3.53)
TiVrms s 1 + Tlcint s 1 + τconv s
ZBASE 1
GlaVrms = · (3.55)
TiVrms s 1 + τconv s
GlaVrms ZBASE
GlcVrms = = (3.56)
1 + GlaVrms ZBASE + TiVrms s + TiVrms τconv s2
TiVrms
TlcVrms = (3.57)
ZBASE
Atendiendo a la ecuación (3.56) y en esta ocasión, sin realizar la simplificación anterior,
mediante la comparativa con la expresión (3.33) como en los casos anteriores, obtenemos
la constante TiVrms , mediante la fórmula (3.58).
1 1
kiVrms = [ ] (3.59)
TiVrms Ωs
1
kpVrms = kiVrms · TnVrms [ ] (3.60)
Ω
El análisis de las unidades expresado en (3.61), demuestra la coherencia de los términos
obtenidos. Donde la variable de entrada es la tensión y la de salida es la intensidad.
Z
1 1
U [A] = kp [ ]E[V ] + kii [ ] Edt[V s] (3.61)
Ω Ωs
1
kiVrms [ Ωs ] 612.48
kpVrms [ Ω1 ] 0.756
Tabla 3.4. Parámetros del regulador PI de la tensión eficaz del punto de conexión a la red
La Figura 3.14 presenta el bloque del regulador PI del valor eficaz de tensión del PCC en
el esquema general de control (a) y su implementación (b). Los parámetros del regulador
se resumen en la Tabla 3.4.
(a) Bloque del regulador PI del valor eficaz (b) Implementación del regulador del valor
de la tensión del PCC, en el esquema general eficaz de la tensión del PCC
Figura 3.14. Regulador PI del valor eficaz de la tensión del punto de acople a la red
1.07
Referencia
1.06 Control
1.05
1.04
Vk(p.u.)
1.03
1.02
1.01
0.99
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Tiempo (s)
Figura 3.15. Respuesta del regulador PI del valor eficaz de la tensión del PCC
TnVDC s + 1 1 1
GlaVDC = GP IVDC · Glc · Gplanta = · · · Gplanta (3.62)
TiVDC s 1 + Tlcint s 1 + τconv s
dVDC
Ic = C (3.63)
dt
Aplicando la transformada de Laplace a la ecuación (3.63), puesto que es un estudio en
el modelo del tiempo continuo, se obtiene (3.64).
VDC (s) 1
Gplanta = = (3.65)
Ic (s) Cs
El método de regulación aplicado en este lazo será el del óptimo simétrico como
consecuencia de la presencia de dos polos en el origen [41]. El empleo de dicha técnica
permite alcanzar respuestas eficientes frente a perturbaciones y no frente a variaciones en
la variable controlada, como si era el objetivo primordial en el método de la cancelación
de la mayor constante de tiempo. Otra caracterı́stica que favorece su elección es que para
el caso del condensador de un convertidor conectado a la red, la constante de planta es
elevada [39]. De la lógica propia de dicho método y atendiendo al diagrama de Bode de
su función de transferencia, se alcanza la expresión aproximada (3.66).
TnVDC s + 1 1
GlaVDC = · (3.66)
TiVDC s Cs
T jw + 1
nVDC c 1 TnVDC
GlaVDC (s = jwc ) = 1 =
· ∼
= (3.67)
TiVDC jwc
|Cjwc | TiVDC Cwc
TnVDC
wc = (3.68)
TiVDC C
2 1
= wc = (3.69)
TnVDC 2Tlcint
Del planteamiento de las expresiones (3.69) y (3.70), se deducen las constantes de tiempo
(3.71) y (3.72).
2 TnVDC
= wc = (3.70)
TnVDC TiVDC C
Tn2V 2
32τconv
TiVDC = DC
= [s · Ω] (3.72)
2C C
A partir de las constantes de tiempo, se logran las componentes directa e integral del
regulador PI de continua (3.73) y (3.74), debido a la relación directa existente entre
ambas variables.
1 1
kiVDC = [ ] (3.73)
TiVDC s·Ω
1
kpVDC = kiVDC · TnVDC [ ] (3.74)
Ω
Para concluir se formula el análisis de las unidades en (3.75), donde vuelve a comprobarse
que los valores obtenidos son razonables. En esta ocasión la variable de entrada es VDC y
la de salida se trata de id .
Z
1 1
U [A] = kp [ ]E[V ] + kii [ ] Edt[V s] (3.75)
Ω Ωs
1
kiVDC [ Ωs ] 395.3
kpVDC [ Ω1 ] 1.952
1.44
Referencia
1.43
Control
1.42
1.41
VDC(p.u.)
1.4
1.39
1.38
1.37
1.36
1.35
1.34
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Tiempo (s)
0.2
0.15
0.1
0.05
PvR(p.u.)
−0.05
−0.1
−0.15
Figura 3.18. Potencia activa del convertidor, generada por la respuesta del regulador PI de la
tensión de la etapa de continua
umax ; u > umax
us = u; umin < u < umax (3.76)
umin ; u < umin
Figura 3.20. Esquema general del regulador PI, donde se estudia el efecto del windup
En la Figura 3.21 (a) se observa que el anti-windup evita que el error siga aumentando,
debido a que produce que la salida del término integral sea nula al alcanza un valor
cercano al lı́mite fı́sico. De este modo, una vez desconectada la carga el regulador retorna
rápidamente al punto de operación idóneo. Sin embargo, en la opción en la que no se ha
saturado el término integral, se genera un efecto indeseado enorme en el regulador que no
deberı́a existir.
700
Salida (con anti−windup) (2)
600
Entrada (con anti−windup) (1)
500 Salida (sin anti−windup) (2)
Entrada (sin anti−windup) (1)
400
Error(e)(p.u.)
300
200
100
−100
−200
−300
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
(a) Comparación de la saturación del término integral en ambos casos
80
Salida (con anti−windup) (4)
60
Entrada (con anti−windup) (3)
40 Salida (sin anti−windup) (4)
Señal saturada(us)(p.u.)
−20
−40
−60
−80
−100
−120
−140
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
(b) Comparación de la saturación en ambos casos
Como se observa en la Figura 3.21 (b), en este último caso el controlador permanece
saturado hasta un segundo después de la desconexión de la carga, a pesar de que el
sistema se situa en un punto de operación donde no se alcanza los lı́mites fı́sicos. El
windup acumulado provoca una respuesta del controlador con un retraso significativo en
su tiempo de estabilización y con perturbaciones considerables, que incluso provocan la
pérdida del control [44], como ocurre en las potencias del caso analizado, como se muestra
en la Figura 3.22.
3 P (sin anti−windup)
Q (sin anti−windup)
2.5
P (con anti−windup)
2
Q (con anti−windup)
Potencias(p.u.)
1.5
0.5
−0.5
−1
Figura 3.22. Efecto del windup en las potencias suministradas por el convertidor
En este estudio se ha optado por frenar la integración cuando los reguladores se saturan,
la implementación del anti-windup ha sido idéntica en todos los reguladores, como se
observa en la Figura 3.23. La ecuación del regulador PI se define en (3.77) [42].
un = KP e + KI η (3.77)
Considerando que:
dη
η̇ ≡ =e
dt
Figura 3.24. Representación esquemática del sistema general con el convertidor real [45]
1
Modulada
0.9 Portadora
0.8 Modular
0.7
0.6
Señales
0.5
0.4
0.3
0.2
0.1
0
0.2 0.2005 0.201 0.2015 0.202 0.2025 0.203
Tiempo (s)
El VSC simula a partir del tren de pulsos la señal sinusoidal, la cual presenta un elevado
contenido de armónicos a frecuencias elevadas producto de la propia conmutación. En la
Figura 3.24 se remarca el filtro pasivo a la salida del propio STATCOM, el cual elimina
o atenúa los armónicos de orden superior de la modulación, como se observa en la Figura
3.27. En dicha imagen se establece la comparación entre la señal de tensión de la red (Vk ),
tras la actuación del filtro paso bajo y la señal de tensión generada por el convertidor
(VvR ), para demostrar la influencia de dicho elemento.
A pesar de que en el punto 2.3.3 se incidió en que un valor elevado de la frecuencia de
conmutación favorece la eficiencia del filtrado, dicho término queda restringido por las
pérdidas de encendido de los interruptores, lo que limita a operar a frecuencias menores
[13] [16].
Los filtros pasivos se clasifican en tres clases principales, los L, los LC y los LCL, donde
todas las componentes reactivas están conectadas a una pequeña resistencia para aumentar
la estabilidad en comparación con los filtros puramente reactivos [13]. En este estudio se
ha optado por el modelo más sencillo, que presenta la ventaja de un bajo coste inicial.
En general, los parámetros principales que se atienden de un filtro son la frecuencia de
resonancia, la distorsión armónica eliminada y la capacidad de compensación reactiva que
permite, este último aspecto se analizó en el punto 2.3.4 [46].
2 Va Vb Vc Va Vb Vc
vR vR vR k k k
1.5
0.5
V(p.u.)
−0.5
−1
−1.5
−2
0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22
Tiempo (s)
1.5
Va Vb Vc Va Vb Vc
ideal ideal ideal real real real
0.5
Vk(p.u.)
−0.5
−1
−1.5
0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22
Tiempo (s)
1.425
1.42
(p.u.)
1.415
DC
V
1.41
1.405
1.4
0.2 0.202 0.204 0.206 0.208 0.21 0.212 0.214 0.216 0.218 0.22
Tiempo (s)
(b) Tensión de la etapa de continua
Resultados de la Simulación
Es este capı́tulo se ofrece un análisis de la respuesta del STATCOM frente a las distintas
perturbaciones planteadas, a raı́z de la implementación analizada en el Capı́tulo 3. En
las secciones 4.1 y 4.2, se muestran los resultados obtenidos en el dominio del tiempo
continuo manteniendo la consideración de que la actuación del convertidor se asemeja
al de una fuente de tensión alterna promediada, dirigida por el control. En los últimos
casos se expone un análisis complementario en el dominio del tiempo discreto donde el
STATCOM opera como un convertidor de dos niveles con modulación sinusoidal por ancho
de pulso (SPWM), con una frecuencia de conmutación de 1620 Hz.
El punto 4.1 muestra un análisis de la actuación del STATCOM como regulador de la
tensión eficaz del PCC, mientras que la sección 4.2 se centra en su función como corrector
del factor de potencia.
En términos generales las cargas industriales presentan generalmente un carácter
inductivo. El STATCOM responde a los cambios de dichas cargas producidos por
maniobras o conmutaciones. Dichas variaciones producen huecos de tensión, lo que se
traduce en gastos y costes en el mantenimiento. Por esta razón, existe un amplio abanico
de investigaciones sobre los motivos o causas que generan los huecos de tensión, sus
consecuencias y efectos [47] [48] [49].
La implementación del modelo de control del STATCOM analizada en el Capı́tulo 3 se
someterá a caı́das breves de tensión en un tiempo inferior a 5 segundos del orden del 95 %
de la tensión nominal, como representación del efecto producido por la conexión de cargas
industriales lineales de carácter inductivo. En su defecto, se planteará también el caso de
desconexiones, ejemplificando la situación que se produce en horarios nocturnos cuando
la industria permanece inactiva provocando un aumento de la tensión indeseado, que en
las simulaciones se situará en torno al 105 % de la nominal.
Las caracterı́sticas del circuito de estudio de la Figura 4.1, se definen en la sección 2.2.3 y
los valores de los elementos empleados se presentan en el Apéndice A. La “carga” puede
ser medida por su intensidad, pero en términos del sistema eléctrico se trata de la potencia
reactiva (Q) de la carga, la cual es directamente responsable de la caı́da de tensión o del
factor de potencia. La carga inicial del sistema será inductiva, dado que es la situación
habitual y su definición se simplifica en su potencia aparente equivalente, compuesta de:
68
PCARGA = 3 [kW ]
La carga auxiliar, que se muestra en la Figura 4.1, sirve como ejemplificación de las
variaciones producidas en el sistema y se define por sus potencias activa y reactiva. Se
ha establecido distintos casos de estudio en los que se remarca detalles relevantes de la
actuación del STATCOM.
En todos los casos realizados, se presenta un resumen del valor de la carga inicial, la carga
auxiliar y por último, el resultado de la interacción de ambas, mediante la nomenclatura:
CI : Carga inicial.
CA : Carga auxiliar.
Tiempo(s)
Caso 1
0-1 1-2 2-3
CI CA CT CI CA CT CI CA CT
Tabla 4.1. Caso 1, como controlador de la tensión eficaz del punto de conexión a la red
El primer caso de estudio consiste en la conexión de una carga inductiva elevada de una
potencia reactiva de 5 kVAr, la cual será desconectada posteriormente para demostrar
la eficiencia del control implementado. Las caracterı́sticas de la carga equivalente quedan
resumidas en la Tabla 4.1.
En la Figura 4.2, se observa que inicialmente la carga hubiera provocado un descenso
del 2 % de la tensión nominal, sin la actuación del STATCOM. No obstante, empleando
un STATCOM sin ningún tipo de restricción se logra mantener la tensión en el valor
nominal. Transcurrido un segundo, la incorporación de la carga auxiliar supone un
aumento significativo de la potencia reactiva inductiva, que produce una caı́da de la
tensión adicional del 2 % de la nominal, acumulando un descenso total del 4 % del nominal,
que el STATCOM ideal es capaz de solventar. Para alcanzar dicho resultado el dispositivo
inyecta potencia reactiva, como se muestra en la Figura 4.4. En dicha representación, se
observa la rápida respuesta del convertidor cuando se retorna al estado inicial, tras la
desconexión de la carga auxiliar.
Otro detalle relevante es que el STATCOM suministra potencia reactiva al PCC por el
carácter inductivo de la carga conectada inicialmente, lo cual se demuestra en el hecho de
que la intensidad está atrasada respecto a la tensión del PCC como se analizó en el punto
2.2.1, dicho desfase se muestra en la Figura 4.3.
1.03
Sin STATCOM
1.02 Con STATCOM
1.01
1
Vk(p.u.)
0.99
0.98
0.97
0.96
0.95
0,5 1 1,5 2 2,5
Tiempo (s)
1
Vk
Variables fundamentales en el PCC(p.u.)
0.8
I
0.6 k
0.4
0.2
−0.2
−0.4
−0.6
−0.8
−1
0 0.1 0.2
Tiempo (s)
1.4
Sin límites
1.3
Límite Ik
1.2
1.1
Q (p.u.)
0.9
vR
0.8
0.7
0.6
0.5
0.4
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
Figura 4.4. Comparación de la potencia reactiva suministrada por el STATCOM, sin lı́mites
fı́sicos o considerando exclusivamente el lı́mite de la intensidad
1.03
Sin límites
1.025
Límite Ik
1.02
1.015
1.01
V (p.u.)
1.005
k
0.995
0.99
0.985
0.98
0,5 1 1,5 2 2,5 3
Tiempo (s)
Figura 4.5. Comparación de la regulación de la tensión del PCC por el STATCOM, sin lı́mites
fı́sicos o considerando exclusivamente el lı́mite de la intensidad
1.2
i
d
1 iq
0.8
Ik(p.u.)
0.6
0.4
0.2
−0.2
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
El control impide sobrepasar el lı́mite de la intensidad máxima soportada por los IGBTs,
mostrado en la Tabla 3.1 del punto 3.4, considerando los valores base del Apéndice
B, la componente de cuadratura de la intensidad (iq ) no deberı́a superar la unidad en
magnitudes por unidad, como se observa en la Figura 4.6. Mientras que como se demuestra
en la Figura 4.7, la potencia activa en régimen permanente es nula, porque el regulador de
la tensión de continua actúa eficazmente manteniéndola constante ante las perturbaciones
del sistema.
0.04
Límite Ik
0.03 Sin límites
0.02
PvR(p.u.)
0.01
−0.01
−0.02
−0.03
0 0,5 1 1,5 2 2,5 3
Tiempo (s)
Figura 4.7. Comparación de la potencia activa suministrada por el STATCOM, sin lı́mites
fı́sicos o considerando exclusivamente el lı́mite de la intensidad
Tiempo(s)
Caso 2
0-1 1-2 2-3
CI CA CT CI CA CT CI CA CT
Tabla 4.2. Caso 2, como controlador de la tensión eficaz del punto de conexión a la red
El segundo caso se resume en la Tabla 4.2, en el que se decide introducir una carga
inductiva auxiliar menos severa que en el caso anterior, de una potencia reactiva de 4
kVAr. En este ejemplo se compara dos modelos del STATCOM, en uno se ha implementado
exclusivamente el lı́mite de la intensidad máxima soportada por los IGBTs (Ik ), mientras
que en el otro se ha considerado también la limitación de la máxima tensión de salida del
convertidor (VvR ) para evitar superar el ı́ndice de modulación de 0.95.
1.3
Límite Ik
1.2
Límite VvR
1.1
1
Q (p.u.)
0.9
vR
0.8
0.7
0.6
0.5
0.4
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
1.03
Límite I
1.025 k
Límite VvR
1.02
1.015
1.01
V (p.u.)
1.005
k
0.995
0.99
0.985
0.98
0 0.5 1 1,5 2 2,5 3
Tiempo (s)
1
Límite I
k
0.9
Límite VvR
0.8
I (p.u.)
0.7
k
0.6
0.5
0.4
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
0.96
0.95
Índice de modulación (ma)
0.94
0.93
0.92
0.91
0.9
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2
Tiempo (s)
Figura 4.11. Índice de modulación del STATCOM, considerando el lı́mite de la tensión máxima
de salida del convertidor
En los dos últimos casos de estudio se muestra los resultados obtenidos a raı́z de la
implementación del convertidor real que opera mediante la técnica de modulación SPWM,
a partir de los pulsos obtenidos por el control implementado.
Tiempo(s)
Caso 3
0-1 1-2 2-3
CI CA CT CI CA CT CI CA CT
Tabla 4.3. Caso 3, como controlador de la tensión eficaz del punto de conexión a la red
1.025
Vk
1.02 real
Vk
1.015 promediada
1.01
1.005
V (p.u.)
1
k
0.995
0.99
0.985
0.98
0.975
0 0.5 1 1,5 2 2,5 3
Tiempo (s)
Figura 4.12. Regulación de la tensión del punto de conexión a la red por el STATCOM,
considerando sus lı́mites fı́sicos
1.6
Q
vR
1.4 real
Q
vR
1.2 promediada
1
QvR(p.u.)
0.8
0.6
0.4
0.2
−0.2
0 0.5 1 1,5 2 2,5 3
Tiempo (s)
Figura 4.13. Potencia reactiva suministrada por el STATCOM, considerando sus lı́mites fı́sicos
Ante una carga auxiliar inductiva que supone una potencia reactiva total de 5.5 kVAr,
el STATCOM operarı́a eficazmente inyectando la potencia reactiva necesaria, como se
observa en la Figura 4.13. Cabe destacar el hecho de que para cumplir dicha finalidad el
STATCOM debe inyectar 9 kVAr de potencia reactiva, sin embargo al PCC solo llega 8
kVAr, como se demuestra en la Figura 4.14.
1.1
Qk
1 real
Qk
promediada
0.9
0.8
Q (p.u.)
0.7
k
0.6
0.5
0.4
Figura 4.14. Potencia reactiva suministrada por el STATCOM tras la actuación del filtro,
considerando sus lı́mites fı́sicos
Otro dato a señalar son las fluctuaciones o variaciones que se produce en la corriente
real, producto de los armónicos con frecuencia inferiores generados por la conmutación
y que no han sido eliminadas por el filtro paso bajo, como se analizó en 3.8. Dichos
armónicos generan errores en las mediciones de las variables fundamentales, utilizadas para
la retroalimentación del control, lo que conlleva una dificultad añadida en la definición de
los lı́mites y en poder impedir de este modo la sobremodulación.
Dicho análisis demuestra la relevancia de optimizar el filtro paso bajo. Además de que
como se ha remarcado en el Capı́tulo 2, es importante aprovechar detalles técnicos
como la topologı́a del convertidor y la técnica de modulación empleada para favorecer
la eliminación interna de dichas fluctuaciones.
14
13
12
Qfiltro(%)
11
10
7
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
Figura 4.15. Porcentaje de la potencia reactiva suministrada por el STATCOM que se pierde
en el filtro
1
I
k
real
0.9 I
k
promediada
0.8
I (p.u.)
0.7
k
0.6
0.5
0.4
0.018 0.5 1 1.5 2 2.5 3
Tiempo (s)
Tiempo(s)
Caso 4
0-1 1-2 2-3
CI CA CT CI CA CT CI CA CT
Tabla 4.4. Caso 4, como controlador de la tensión eficaz del punto de conexión a la red
El último caso resumido en la Tabla 4.4, se trata de la conexión de una carga capacitiva
que provoca que el STATCOM que suministraba potencia reactiva, cambie de carácter y
comience absorberla en el intervalo de tiempo de 1 a 2 segundos. En la Figura 4.17 se
comprueba que el STATCOM al no alcanzar los lı́mites fı́sicos responde adecuadamente.
En dicho periodo de tiempo, la carga total equivale a una potencia capacitiva de 11 kVAr,
donde la potencia absorbida por el STATCOM no alcanza los 8 KVAr, como se demuestra
en la Figura 4.17, esto se debe a las pérdidas de potencia reactiva que se producen en los
elementos empleados en la conexión a la red, como el filtro.
2
Q
vR
real
1.5
Q
vR
promediada
1
QvR(p.u.)
0.5
−0.5
−1
−1.5
0 0.5 1 1,5 2 2,5 3
Tiempo (s)
Figura 4.17. Potencias absorbida por el STATCOM, considerando sus lı́mites fı́sicos
1
Vk Ik
0.8 real real
0.6
0.4
0.2
−0.2
−0.4
−0.6
−0.8
−1
1.4 1.41 1.42 1.43 1.44 1.45 1.46 1.47 1.48 1.49 1.5
Tiempo (s)
1.4
Vk VvR Vk VvR
real real promediada promediada
1.3
1.2
1.1
V (p.u.)
0.9
0.8
0.7
Figura 4.19. Magnitud de la tensión del PCC por la actuación del STATCOM, considerando
sus lı́mites fı́sicos
id iq id iq
1 real real promediada promediada
0.8
0.6
0.4
Ik(p.u.)
0.2
−0.2
−0.4
−0.6
−0.8
−1
0 0.5 1 1,5 2 2,5 3
Tiempo (s)
Figura 4.20. Intensidad que circula entre el PCC y el STATCOM, considerando sus lı́mites
fı́sicos
Este hecho unido al carácter inductivo de la red y las lı́neas de transmisión, provoca que si
se compara el caso 3 y 4, se comprueba que cuando la carga total equivalente es capacitiva,
el STATCOM es capaz de compensar una mayor cantidad de energı́a reactiva.
Por último, en la Figura 4.20 se observa que las fluctuaciones de las componentes de
la intensidad real no son significativas. En dicha representación se muestra como la
componente de cuadratura (iq ) cambia de sentido, puesto que el STATCOM deja de
suministrar potencia reactiva para comenzar a absorberla, mientras que la componente
directa (id ) no varı́a su signo, porque la red sigue suministrando la potencia activa
consumida por los dispositivos utilizados en la conexión del STATCOM al PCC, tales
como el filtro.
Tiempo(s)
Caso 1
0-1 1-2 2-3
CI CA CT CI CA CT CI CA CT
El primer caso de estudio se resume en la Tabla 4.5, en el que se conecta una carga
inductiva auxiliar que consume una potencia reactiva de 9 kVAr y una potencia activa de
2 kW. La carga auxiliar se desconectará posteriormente transcurrido un intervalo de un
segundo, para demostrar la eficiencia y dinámica del control implementado.
Con STATCOM
1.2
Sin STATCOM
Factor de potencia
0.8
0.6
0.4
0.2
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
Figura 4.21. Demostración de la actuación eficiente del STATCOM, como corrector del factor
de potencia
En la Figura 4.21, se analiza como el empleo del STATCOM permite obtener un factor
de potencia unitario, a pesar de que inicialmente si no se empleara dicho dispositivo
este presentarı́a un valor inferior a 0.8. Incluso tras la incorporación de la carga auxiliar
inductiva en el intervalo transcurrido entre los segundos 1 y 2, dicho término serı́a inferior
a 0.4.
Cabe resaltar el detalle de que al mantenerse la componente de la potencia activa de la
carga auxiliar que consume 2 kW, el factor de potencia es mayor en el periodo final donde
la potencia activa equivalente es de 5 kW, mientras que en el instante inicial era de 3
kW. De este modo, se demuestra que en tres franjas de tiempo con una carga equivalente
de caracterı́sticas distintas, el STATCOM cumple el objetivo de un factor de potencia
unitario.
Dicho factor de potencia supone que la potencia reactiva del generador es nula, operando
a plena capacidad para suministrar exclusivamente la potencia activa o útil exigida por
el sistema, como se observa en la Figura 4.22 (a) y como queda expuesto en la sección
2.2.1. Sin embargo, en dicha representación se muestra un pequeño error despreciable que
provoca que la potencia reactiva del generador no sea del todo nula, esto es consecuencia
de las pérdidas reactivas de los elementos empleados en la conexión al PCC, las cuales
dependen de las variaciones de las tensiones y las intensidades, estas últimas especialmente
significativas en el intervalo de 1 a 2 segundos.
A pesar de que la carga equivalente consume una potencia activa de 5 kW, el generador
suministra casi 7 kW como consecuencia de tener que subsanar las pérdidas de los
elementos conectados al PCC, especialmente el filtro. Este elemento supone casi 2 kW
adicionales que el generador debe subsanar, cuando la corriente es elevada, como se observa
en la Figura 4.22 (b).
P
generador
0.8
Qgenerador
Potencias (p.u.)
0.6
0.4
0.2
−0.2
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
(a) Potencias suministrada por la red tras la actuación del STATCOM ideal, sin lı́mites fı́sicos
0.8
0.6
0.4
Pfiltro(p.u.)
0.2
−0.2
−0.4
−0.6
−0.8
−1
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
(b) Potencia activa consumida en el filtro
Un dato significativo es que para corregir el factor de potencia de una carga inductiva
de una potencia reactiva equivalente de 12 KVAr, el STATCOM inyecta una potencia
reactiva algo superior a 14 KVAr, como se observa en la Figura 4.23, mientras que en el
primer caso del anterior punto 4.1, para una carga inductiva bastante inferior de 8 KVar,
el STATCOM necesitaba inyectar más de 13 KVAr, como se presentaba en la Figura
4.4. El STATCOM necesita inyectar menos potencia reactiva para corregir el factor de
potencia que para regular la tensión eficaz del PCC. La explicación teórica de que cumplir
esta función requiera menos energı́a, se debe a que para regular el factor de potencia en
su valor unitario, solo se necesita compensar la carga inductiva y lograr una potencia
reactiva nula en el PCC, mientras que mantener la tensión del PCC en el valor nominal
exige un esfuerzo adicional de potencia reactiva. Otro motivo menos relevante es que la
intensidad necesaria entre el STATCOM y la red para corregir el factor de potencia es
menor y por lo tanto, se produce una menor cantidad de pérdidas en la conexión al PCC.
1.6
Límite Ik
1.4 Sin límites
1.2
1
Q (p.u.)
0.8
vR
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
Figura 4.23. Comparación de la potencia reactiva suministrada por el STATCOM, sin lı́mites
fı́sicos o considerando exclusivamente el lı́mite de la intensidad
Límite I
k
0.6 Sin límites
0.4
PvR(p.u.)
0.2
−0.2
−0.4
Figura 4.24. Comparación de la potencia activa suministrada por el STATCOM, sin lı́mites
fı́sicos o considerando exclusivamente el lı́mite de la intensidad
La implementación del lı́mite de la máxima intensidad soportada por los IGBTs se cumple,
puesto que si se emplea los valores base del Apéndice B en dicho valor mostrado en la
Tabla 3.1 del punto 3.4, la componente de cuadratura de la intensidad (iq ) no deberı́a
superar la unidad en magnitudes por unidad, como se muestra en la Figura 4.25.
1
id
iq
0.8
0.6
I (p.u.)
k
0.4
0.2
Tiempo(s)
Caso 2
0-1 1-2 2-3
CI CA CT CI CA CT CI CA CT
El segundo caso se resume en la Tabla 4.6, en el que se introduce una carga inductiva
auxiliar menor de una potencia reactiva de 7 kVAr. Este punto sirve para establecer
una comparación entre dos modelos, en el primero se ha implementado el lı́mite de la
máxima intensidad soportada por los IGBTs y en el segundo se incluye la restricción de
la máxima tensión de salida del STATCOM permitida para evitar la sobremodulación.
Como se muestra en la Figura 4.26, si el STATCOM presentará solo la restricción de
la intensidad, cumplirı́a el objetivo de corregir el factor de potencia, no obstante cabe
señalar que el modelo donde se implementa el lı́mite de la máxima tensión de salida
del STATCOM aunque no cumple dicho objetivo, logra un factor de potencia aceptable
próximo a la unidad.
1.05
Límite Ik
1.04
Límite VvR
1.03
Factor de potencia
1.02
1.01
0.99
0.98
0.97
0.96
0.95
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
1.4
Límite V
vR
1.2 Límite Ik
0.8
QvR
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3
Tiempo (s)
El motivo por el que al considerar dicha restricción no cumpla su finalidad, se debe a que
el STATCOM no tiene la capacidad de suministrar la potencia reactiva necesaria para
compensar totalmente la carga, como se observa en la Figura 4.27. No obstante, el control
mantiene la tensión de continua de ambos modelos constante, logrando que la potencia
activa del convertidor sea nula, como se muestra en la Figura 4.28.
0.2
P (p.u.)
0
vR
−0.2
−0.4
−0.6
0.98
Índice de modulación (ma)
0.96
0.94
0.92
0.9
0.88
1 1.2 1.4 1.6 1.8 2
Tiempo (s)
Figura 4.29. Índice de modulación del STATCOM, considerando ambos lı́mites fı́sicos
Los dos últimos ejemplos presentados en 4.2 se han simulado en un periodo de dos
segundos, donde se añade la actuación convertidor real que emplea la técnica de
modulación SPWM. El tercer caso se resume en la Tabla 4.7. En la Figura 4.30, se
demuestra que el STATCOM cumple la finalidad de mantener el factor de potencia en su
valor unitario, sin sobrepasar los lı́mites fı́sicos ante una carga auxiliar inductiva de una
potencia reactiva total de 9 kVAr. En el caso real, la propia modulación genera armónicos
en la potencia proporcionada por el convertidor provocando que el factor de potencia no
sea compensado de manera uniforme y constante, como si ocurre en el caso promediado.
Tiempo (s)
Caso 3
0-1 1-2
CI CA CT CI CA CT
Real
1.15
Promediada
1.1
Factor de potencia
1.05
0.95
0.9
0.85
0.8
0.75
0.7
0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8
Tiempo (s)
Figura 4.30. Factor de potencia proporcionado por el STATCOM, considerando sus lı́mites
fı́sicos
2
Q
vR
real
Q
vR
promediada
1.5
Q (p.u.)
1
vR
0.5
Figura 4.31. Potencias suministradas por el STATCOM, considerando sus lı́mites fı́sicos
1.5
Qk
real
Qk
promediada
1
Qk(p.u.)
0.5
Figura 4.32. Potencias suministradas por el STATCOM tras la actuación del filtro,
considerando sus lı́mites fı́sicos
15
14
13
12
11
Qfiltro(%)
10
5
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Tiempo (s)
Figura 4.33. Porcentaje de la potencia reactiva suministrada por el STATCOM que se pierde
en el filtro
Ik
1.2 real
Ik
promediada
1
I (p.u.)
0.8
k
0.6
0.4
0.2
0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8
Tiempo (s)
Figura 4.34. Magnitud de la intensidad que circula entre el punto de conexión a la red y el
STATCOM
Para una carga inductiva equivalente pero con el STATCOM actuando como regulador de
la tensión del PCC, el valor era ligeramente inferior al 8 % como se presenta en la Figura
4.16, porque la corriente era de 0.45 en magnitudes por unidad. Esta diferencia recae en la
necesidad de una potencia reactiva e intensidad menor para cumplir la función de corregir
el factor de potencia, respecto a la requerida para regular la tensión del PCC en su valor
nominal.
Tiempo (s)
Caso 4
0-1 1-2
CI CA CT CI CA CT
El último caso resumido en la Tabla 4.8, refleja la actuación del STATCOM ante la
aparición de una carga capacitiva, que para lograr su compensación deja de operar como
inversor y comienza actuar como rectificador absorbiendo energı́a en el intervalo de tiempo
de 1 a 2 segundos, como se ejemplifica en la Figura 4.35.
0.5
Q
k
real
Qk
promediada
0
Qk(p.u.)
−0.5
−1
0 0,5 1 1,5 2
Tiempo (s)
Figura 4.35. Potencias suministradas por el STATCOM, considerando sus lı́mites fı́sicos
1.3
Real
1.2 Promediada
1.1
Factor de potencia
0.9
0.8
0.7
0.6
0.5
0 0.5 1 1.5 2
Tiempo (s)
(a) Corrector del factor de potencia por el STATCOM, considerando los lı́mites fı́sicos
3
Pgen
real
2.5
Qgen
real
2
P
gen
Potencias(p.u.)
promediada
1.5 Qgen
promediada
0.5
−0.5
−1
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Tiempo (s)
(b) Potencias suministrada por el generador tras la actuación del STATCOM, considerando los
lı́mites fı́sicos
Conclusión
Este proyecto sirve como representación del control de las dos principales funciones de
un STATCOM conectado a una red. El siguiente capı́tulo es un resumen de los objetivos
cumplidos, posibles mejoras en el control del modelo implementado y por último, algunos
de los trabajos futuros que podrı́an realizarse para ampliar la información ofrecida en este
estudio.
97
5.2. Posibles mejoras del modelo presentado
Presupuesto
En este capı́tulo se estima una aproximación de los costes que supondrı́a el proyecto
realizado desde el punto de vista empresarial. En el punto 6.1 se elabora el presupuesto
personal de la realización del estudio “Estudio y Simulación de un STATCOM conectado
a una red”. Mientras que en la sección 6.2 se ofrece un análisis de los factores a considerar
a la hora de estudiar la inversión de dicho dispositivo.
100
6.1. Presupuesto personal.
50 · 6 + 133 · 6
Coste mensual Autonomo = = 91,5 [e/mes] (6.1)
12
2. Costes operacionales: Este punto incluye los recursos informáticos necesarios para
la elaboración de este estudio, desde los dispositivos de trabajo como el ordenador
portátil hasta el pago de las licencias de los programas informáticos empleados.
Atendiendo a la referencia [54] se obtiene los precios de las licencias de Matlab
y Simulink sin profundizar en los paquetes complementarios, se estima un coste
mensual aproximado de ambos mediante las expresiones (6.2) y (6.3).
2000
Coste mensual M atlab = = 166,67 [e/mes] (6.2)
12
3000
Coste mensual Simulink = = 250 [e/mes] (6.3)
12
Otro objeto de análisis son los costes indirectos tales como el pago de la hipoteca del
crédito pedido para la compra de un local, que se utiliza como oficina y los gastos
derivados de su utilización, como la comunidad, el agua y la luz.
3. Subcontratación: Muchas empresas requieren la necesidad de contratar servicios
secundarios, que aunque están directamente relacionados con su actividad, no son
la especialidad a la que se dedican y les resulta más rentable optar por la vı́a de
la subcontratación. En el caso planteado ante el montaje de una pequeña empresa
nueva, contratar una gestorı́a para encargarse del tedioso trabajo administrativo,
es una situación lógica y común, por lo que se ha considerado en el presupuesto
planteado.
Una vez calculado los coste mostrados en la Tabla 6.1, se elabora una factura como la
mostrada en la Tabla 6.2, donde se incluye el 21 % del IVA mediante la expresión (6.4),
obteniendo el importe final a cobrar por el servicio.
1. Personal
2. Costes operacionales
2.1 Hardware
2.2 Software
3. Subcontratación
1. Personal 9823.5
3. Subcontratación 270
El STATCOM ofrece entres sus principales aplicaciones regular la tensión del PCC
y corregir el factor de potencia. Sin embargo, aunque la calidad de la energı́a queda
ampliamente mejorada por su actuación, la generalización del STATCOM queda lastrada
por la inversión de capital necesaria. Se ofrece un resumen de las principales premisas
económicas que se atienden a la hora de elaborar un análisis de la inversión [50]. El
objetivo final es extraer una conclusión fidedigna mediante el empleo de una herramienta
de análisis económico ampliamente conocida, VAN (Valor actual neto). Dicho parámetro
es el producto de los beneficios y costes relacionados con la implantación del propio
STATCOM.
Un factor clave es la potencia del dispositivo, que determina su capacidad para cumplir
la función fijada por el cliente y por consiguiente, su precio. Los principales beneficios
de su operación han quedado remarcados en la sección 2.2.1, pero cabe destacar otros
efectos secundarios, como reducir los costes laborales de errores producidos por el trabajo
en condiciones desfavorables (intensidades elevadas) y producen pérdidas del valor de la
producción, el gasto de reconstruir o reparar los materiales perdidos o el equipamiento
dañado, además de los ingresos (oportunidades) perdidos por no cumplir el contrato.
Los principales factores a considerar en los beneficios son el cumplimiento de las siguientes
funciones:
Regular la tensión del punto de conexión, analizado en 4.1. En este apartado se
incide en considerar la eficacia del STATCOM con los parámetros Want y Wpos . El
beneficio del empleo del STATCOM se determina con la expresión (6.5) [50].
Donde:
• Φsag : Beneficio debido a la mitigación del sag de tensión [e].
• ksag : Coste unitario por pérdidas [e/kWh].
• Ccons : Capacidad del consumidor [kW]
• tint : Duración de la interrupción [h].
• Want , Wpos : Factores adimensionales empleados para estimar la variación de
magnitud antes y después de la mitigación.
Corregir el factor de potencia en el PCC, analizado en 4.2. El beneficio de la
reducción del factor de potencia queda expresado mediante las ecuaciones (6.6)
y (6.7) [50].
ΦP F = kP F · Qcom (6.6)
s s
2 2
1 − P Fant 1 − P Fval
Qcom = PL ( 2
− 2
) (6.7)
P Fant P Fval
Donde:
• ΦP F : Beneficio de corregir el factor de potencia [e].
• Qcom : Cantidad de potencia reactiva compensada [kVAr]
• PL : Potencia activa de la carga [kW].
• P Fant : Factor de potencia de la instalación anterior.
• P Fval : Mı́nimo factor de potencia aceptable.
• kP F : Coste de un bajo factor de potencia [e/kWh].
Otro factor son las pérdidas energéticas evitadas por el STATCOM, con la suma de dichos
parámetros se establecerı́a los beneficios de su instalación. Desde el punto de vista de los
gastos, se debe analizar la contratación del personal necesario, dicho gasto se analiza en
un presupuesto similar al presentado en el punto 6.1, el importe inicial de la compra y en
menor medida los gastos de mantenimiento y operación.
Los costes fijos de la inversión inicial del STATCOM pueden relacionarse exclusivamente
con su capacidad de potencia reactiva [50], sin embargo como se ha estudiado en este
proyecto dicha inversión se fundamenta en las dimensiones del condensador y los IGBTs,
además de la configuración del convertidor que influye en el número de componentes
empleados. Para determinar dicho precio se debe realizar un estudio del mercado, tras la
definición de las condiciones de operación requeridas. Además se debe analizar los gastos
de mantenimiento definidos en (6.8) [55].
tX
vida
100 + λ
Cop = Cman · ( ) (6.8)
t=1
100
Donde:
Cop : Costes operacional total [e].
Cman : Costes anual de mantenimiento [e].
λ : Porcentaje anual de incremento del coste de mantenimiento [ %].
tvida : Ciclo de vida estimado para el STATCOM en condiciones estipuladas [años].
La teorı́a del análisis de los costes es mucho más compleja que el resumen ofrecido en
este punto, deberı́a ser considerado el valor tiempo del dinero, la necesidad de requerir
prestamos para la inversión con sus correspondientes intereses, la legislación vigente que
se debe cumplir y los costes de no llegar a los plazos. El sumatorio de todos las cuestiones
comentadas acrecientan la cuantı́a final de la instalación. El VAN es el resultado de
comparar costes y beneficios, siendo un parámetro fiable para aceptar o desechar la idea
de un proyecto [56], el cual se define en la expresión (6.9).
n
X Bi − Ci i
V AN = −Cf + ( r ) (6.9)
t=1
1 + 100
Donde:
Cf : Costes fijos [e].
Bi : Beneficios totales del año i [e].
Ci : Costes totales del año i [e].
r : Índice de rentabilidad.
n : Periodo de análisis [años].
Del análisis presentado se puede extraer la conclusión de que una inversión de esta ı́ndole
será rentable, en primera instancia si el punto de la distribución energética donde se
acople el dispositivo presenta unas condiciones deplorables desde el punto de vista del
consumidor, por una gran cantidad de incidencias que perjudican la calidad energética.
Otro factor determinante, es elaborar un estudio del dimensionamiento del STATCOM
necesario atendiendo a su capacidad, eficacia y vida útil.
2015 / 2016
Dic Ene Feb Mar Abr May Jun Jul Ago Sep Oct
1. Documentación
1.1. Búsqueda
1.2. Comprensión
2. Implementación
2.1. Estudio PLL
2.2. Planta real
2.3. Reguladores
2.4. Flujo de carga
2.5. SPWM
2.6. Lı́mites
3. Memoria
3.1. Introducción
3.2. Estado del arte
3.3. Control (1)
3.4. Control (2)
3.5. Simulaciones
3.6. Conclusión
3.7. Presupuesto
4. Corrección
4.1. Revisión
4.2. Corrección
5. Presentación
5.1. Preparación
5.2. Exposición
106
Bibliografı́a
[2] Xiao-Ping Zhang, Christian Rehtanz y Bikash Pal. Flexible AC transmission systems:
modelling and control. Springer Science & Business Media, 2006.
[3] Enrique Acha, Vassilios Agelidis, Olimpo Anaya y Timothy John Eastham Miller.
Power Electronic Control in Electrical Systems. Elsevier, 2001.
[4] Codra. Alstom Grid uses a SCADA system for its power compen-
sators. http://www.energy.siemens.com/mx/en/power-transmission/facts/
static-var-compensator-plus/#content=References.
[7] Timothy John Eastham Miller y Charles Concordia. Reactive power control in electric
systems, volume 2. Wiley New York, 1982.
[8] N. G. Hingorani. FACTS Technology - State of the Art, Current Challenges and
the Future Prospects. In Power Engineering Society General Meeting, 2007. IEEE,
pages 1–4, June 2007.
[10] ABB. TCSC for stable transmission of surplus power from Eastern to Western In-
dia. https://library.e.abb.com/public/1e36fc1bf524064ec125787800258583/
A02-0185%20E_Raipur_LR.pdf.
[11] Aysen Arsoy, Yilu Liu, Shen Chen, Zhiping Yang, ML Crow y PF Ribeiro. Dynamic
performance of a static synchronous compensator with energy storage. In Power
Engineering Society Winter Meeting, 2001. IEEE, volume 2, pages 605–610. IEEE,
2001.
107
[12] C. K. Sao, P. W. Lehn, M. R. Iravani y J. A. Martinez. A benchmark system for
digital time-domain simulation of a pulse-width-modulated D-STATCOM. IEEE
Transactions on Power Delivery, 17(4):1113–1120, Oct 2002.
[13] Farhad Shahnia, Sumedha Rajakaruna y Arindam Ghosh. Static Compensators
(STATCOMs) in Power Systems. Springer, 2015.
[14] Marı́a Ángeles Martı́n Prats. Nuevas técnicas de modulación vectorial para
convertidores electrónicos de potencia multinivel. Tesis Doctoral, Universidad de
Sevilla, 2003.
[15] Jih-Sheng Lai y Fang Zheng Peng. Multilevel converters-a new breed of power
converters. IEEE Transactions on industry applications, 32(3):509–517, 1996.
[16] Bingsen Wang y Jimmie J Cathey. DSP-controlled, space-vector PWM, current
source converter for STATCOM application. Electric Power Systems Research,
67(2):123–131, 2003.
[17] Saswato R Das. The sensible superconductor [magnesium diboride]. IEEE Spectrum,
39(7):34–37, 2002.
[18] P Sanchis, Ó Alonso, A Guerrero y L Marroyo. IGBT, caracterı́sticas y utilización.
MUNDO ELECTRONICO, pages 80–80, 1999.
[19] López Mesa Diana Jimena, Camacho Muñoz Guillermo Alberto, Dı́az Chávez, Jaime
Oscar y Gaviria López Carlos Alberto. Modulación PWM aplicada a inversores
trifásicos dentro del esquema de accionamientos eléctricos AC., 2007.
[20] Carlos Mancebo Durán. Control de motores eléctricos con baterı́as electroquı́micas.
Trabajo Fin de Grado, Universidad Carlos III de Madrid., 2015.
[21] A Schonung y H Stemmler. Static frequency changers with subharmonic control in
conjunction with reversible variable speed ac drives. Brown Boveri Rev, 51(8/9):555–
577, 1964.
[22] Pedro Rodrı́guez Cortés. Aportaciones a los acondicionadores activos de corriente
en derivación para redes trifásicas de cuatro hilos. Tesis Doctoral, Universitat
Politècnica de Catalunya, 2005.
[23] Joachim Holtz. Pulsewidth modulation for electronic power conversion. Proceedings
of the IEEE, 82(8):1194–1214, 1994.
[24] SL Capitaneanu, B de Fornel, M Fadel, J Faucher y A Almeida. Control-Graphical
and algebraic synthesis for PWM methods. EPE-European Power Electronics and
Drives Journal, 11(3):16–28, 2001.
[25] Alfredo De la Barrera González, AA Arellano Sánchez, Indrid Reséndiz Maldonado
y Marco A Javier Sánchez. Análisis de PWM sinusoidal (SPWM) utilizando Matlab
para generar el código de conmutación digital. IEEE ROC&C, 2015.
[26] Mónica Chinchilla Sánchez. Control de un sistema de generación eólica de velocidad
variable con generador sı́ncrono multipolar de imanes permanentes acoplado a red.
Tesis Doctoral, Universidad Carlos III de Madrid, 2001.
108
[27] Robert H Park. Two-reaction theory of synchronous machines generalized method
of analysis-part i. Transactions of the American Institute of Electrical Engineers,
48(3):716–727, 1929.
[28] Eugenie Clark. A method for artificial insemination in viviparous fishes. Science
(New York, NY), 112(2920):722, 1950.
[29] Joaquı́n Eloy-Garcı́a Carrasco. Control directo de potencia de convertidores
electrónicos conectados a la red. Tesis Doctoral, Universidad Carlos III de Madrid.
Departamento de Ingenierı́a Eléctrica,, 2007.
[30] D. Dong, B. Wen, D. Boroyevich, P. Mattavelli y Y. Xue. Analysis of Phase-Locked
Loop Low-Frequency Stability in Three-Phase Grid-Connected Power Converters
Considering Impedance Interactions. IEEE Transactions on Industrial Electronics,
62(1):310–321, Jan 2015.
[31] Yankui Zhang, Yan Zhang, Bei Wu y Jian Zhou. Power injection model of STATCOM
with control and operating limit for power flow and voltage stability analysis. Electric
Power Systems Research, 76(12):1003–1010, 2006.
[32] Claudio A Canizares, Massimo Pozzi, Sandro Corsi y Edvina Uzunovic. STATCOM
modeling for voltage and angle stability studies. International Journal of Electrical
Power & Energy Systems, 25(6):431–441, 2003.
[33] Shigeo Morimoto, Yoji Takeda, Takao Hirasa y Katsunori Taniguchi. Expansion of
operating limits for permanent magnet motor by current vector control considering
inverter capacity. IEEE Transactions on Industry Applications, 26(5):866–871, 1990.
[34] John G Ziegler y Nathaniel B Nichols. Optimum settings for automatic controllers.
trans. ASME, 64(11), 1942.
[35] Joseba Quevedo y Teresa Escobet. Digital control: Past, present and future of PID
control. Elsevier Science Inc., 2000.
[36] Yun Li, Kiam Heong Ang y G. C. Y. Chong. PID control system analysis and design.
IEEE Control Systems, 26(1):32–41, Feb 2006.
[37] P Cominos y N Munro. PID controllers: recent tuning methods and design to
specification. IEE Proceedings-Control Theory and Applications, 149(1):46–53, 2002.
[38] Y Li, W Feng, XK Zhu, KC Tan, X Guan y KH Ang. PIDeasy and automated
generation of optimal PID controllers. Asia-Pacific Conference on Control &
Measurement, 3 rd, Dunhuang, China, pages 29–33, 1998.
[39] Universidad Carlos III de Madrid. Notas de la materia Control de Máquinas
Eléctricas, 2016.
[40] Athanasios Papoulis. Sistemas y circuitos: digitales y analógicos. Marcombo, 1989.
[41] Jeffrey W Umland y Mohammed Safiuddin. Magnitude and symmetric optimum
criterion for the design of linear control systems: what is it and how does it compare
with the others? IEEE Transactions on Industry Applications, 26(3):489–497, 1990.
109
[42] A. Scottedward Hodel y C. E. Hall. Variable-structure PID control to prevent
integrator windup. IEEE Transactions on Industrial Electronics, 48(2):442–451, Apr
2001.
[43] D Vrancic. Design of anti-windup and bumpless transfer protection. Part II:
PID controller tuning by multiple integration method. Tesis Doctoral, PhD thesis,
University of Ljubljana, J. Stefan Institute, Ljubljana, Slovenia, 1996.
[44] Antonio Visioli. Modified anti-windup scheme for PID controllers. IEE Proceedings-
Control Theory and Applications, 150(1):49–54, 2003.
[45] Guopeng Zhao, Jinjun Liu y Minxiao Han. Design of input filters considering the
stability of STATCOM systems. Journal of Power Electronics, 11(6):904–913, 2011.
[46] Cunping Wang, Xianggen Yin, Minghao Wen, Jian Liu, Qing Xiong y Bin Zhang.
Structure and parameters design of output LC filter in D-STATCOM. In Power
System Technology (POWERCON), 2010 International Conference on, pages 1–6.
IEEE, 2010.
[47] Math HJ Bollen, Irene YH Gu, Peter GV Axelberg y Emmanouil Styvaktakis.
Classification of underlying causes of power quality disturbances: deterministic versus
statistical methods. EURASIP Journal on Applied Signal Processing, 2007(1):172–
172, 2007.
[48] Emmanouil Styvaktakis, Irene YH Gu y Math HJ Bollen. Classification of power
system events: voltage dips. In Harmonics and Quality of Power, 2000. Proceedings.
Ninth International Conference on, volume 2, pages 745–750. IEEE, 2000.
[49] William E Brumsickle, Robert S Schneider, Glen A Luckjiff, Deepak M Divan y
Mark F McGranaghan. Dynamic sag correctors: cost-effective industrial power line
conditioning. IEEE Transactions on Industry Applications, 37(1):212–217, 2001.
[50] Kittaya Somsai y Thanatchai Kulworawanichpong. Cost estimation for reactive
power compensation in distribution power system by using D-STATCOM, 2013.
[51] Zhiping Yang, Mariesa L Crow, Chen Shen y Lingli Zhang. The steady state
characteristics of a StatCom with energy storage. In Power Engineering Society
Summer Meeting, 2000. IEEE, volume 2, pages 669–674. IEEE, 2000.
[52] Arturo Zuniga. Manual de presupuesto para proyectos. http://es.slideshare.
net/cpczuniga/manual-de-presupuesto-para-proyectos.
[53] Ministerio de Empleo y Seguridad Social. Base de cotización. http:
//www.seg-social.es/Internet_1/Trabajadores/CotizacionRecaudaci10777/
Regimenes/RegimenEspecialTrab10724/index.htm.
[54] MathWorks. Precio de licencias individuales. http://es.mathworks.com/
pricing-licensing/index.html?prodcode=ML.
[55] TS Chung y Ge Shaoyun. A recursive lp-based approach for optimal capacitor
allocation with cost-benefit consideration. Electric Power Systems Research,
39(2):129–136, 1996.
110
[56] Richard Robinson, Uno Danielson y Martin Snaith. Road Maintenance Management-
Concepts and Systems, 1998.
111
Apéndices
112
Apéndice A
Frecuencia fundamental f= 60 Hz
113
Apéndice B
Valores base
2
I AC base = 3
Sbase /Vbase [A]
114
Apéndice C
Transformación al sistema de
referencia d-q-o
115
Un aspecto fundamental para el desarrollo teórico del tema en cuestión son los métodos de
descomposición de la corriente y la teorı́a de potencia. Los enfoques dominantes fueron los
de Fryze [7] y Budeanu [8], a finales de los años 20. La caracterización de la potencia en el
dominio de la frecuencia fue planteada por Budenau, quien definió a raı́z de la interacción
de la tensión e intensidad, tres componentes, la activa P, la reactiva Q y por último, una
potencia de distorsión D. Fryze empleo dicha representación para el dominio temporal,
como queda ampliamente expuesto en la referencia [1].
Respecto al enfoque vectorial del sistema trifásico fue Fortescue [9] quien en 1918 introdujo
el concepto de las componentes simétricas. Dicho método posibilitó estudiar sistemas
trifásicos incluso en un estado de desequilibrio, con un análisis matemático y fı́sico de los
fenómenos producidos en dichas condiciones.
Lyon [10] planteó la transformación de las variables, considerando el ejemplo de la tensión,
obteniendo la siguiente expresión:
+
~v va
~v − = [T+−0 ] vb (C.2)
~v 0 vc
Donde:
2π
α = ej 3
Se considera [1]:
va ,vb y vc representan las variables instantáneas de las tres fases.
α es el operador de Fortescue.
La norma de la matriz de transformación [T+−0 ] es igual a uno.
La variable instantánea v 0 es siempre real.
Las variables instantáneas ~v + y ~v − y pueden representarse como vectores espaciales
cuyas trayectorias en el plano complejo dependen de funciones temporales va ,vb y
vc .
Las variables ~v + ,~v − y v 0 no deben confundirse con los vectores instantáneos de
secuencia positiva, negativa y homopolar.
Clarke [5] empleo elementos reales en lugar de complejos en su matriz de transformación,
expreso dicho sistema referencia con la nomenclatura α − β − γ, como queda evidenciado
en las fórmula (C.3)-(C.6).
116
V~αβγ = [Tαβγ ] V~abc (C.3)
√ +
vα √2<(~v )
vβ = 2=(~v + ) (C.4)
vγ v0
1 1
1 −
√ 2
−√2
2
[Tαβγ ] = · 0 23 − 23 (C.5)
3 1 1 1
2 2 2
1 0 1
√
[Tαβγ ]−1 = [Tαβγ ]T = − 21 2√
3
1 (C.6)
1 3
−2 − 2 1
Park [6] extendió la transformación de Clarke, para expresar las variables de fase sobre
unos ejes de referencia rotativos (d-q) que giraban a velocidad w, y ocupaban la posición
θ = wt a lo largo del tiempo. La transformación propuesta por Park se resume en las
expresiones (C.7)-(C.10):
vd va
vq = [Tdqz ] [Tαβγ ] vb (C.8)
vz vc
cos(θ) sin(θ) 0
[Tdqz ] = − sin(θ) cos(θ) 0 (C.9)
0 0 1
cos(θ) cos(θ − 2π/3) cos(θ + 2π/3)
2
[Tθ ] = − sin θ − sin(θ − 2π/3) − sin(θ + 2π/3) (C.10)
3 1 1 1
2 2 2
En la Figura C.1 se observan las componentes del vector V~ sobre el sistema de referencia
rotativo d-q-z. Esta concepción del sistema trifásico resulta de gran utilidad en la definición
de potencias instantáneas. Las teorı́as posteriores más significativas que contribuyeron al
desarrollo del enfoque vectorial en la caracterización de intensidades y potencias, no serán
remarcadas en este proyecto, pero son analizadas en la referencia [1].
117
Figura C.1. Sistema de coordenadas α -β -γ y d-q-z [1]
118
Bibliografı́a
119
Apéndice D
El PLL es una técnica empleada para lograr la sincronización del control con la frecuencia
fundamental de la planta. Para cumplir dicha finalidad se emplea el vector de la tensión,
cuyo ángulo de fase contiene información imprescindible, destacando los sistemas que
emplean electrónica de potencia. No obstante, existe otros factores determinantes tales
como las dificultades y problemas originados en el acoplamiento entre la actuación del
PLL y la impedancia de la carga [1].
En términos generales la detección de la secuencia positiva de la tensión es una cuestión
crucial, por ejemplo en sistemas flexibles de trasmisión en corriente alterna (FACTS) [2].
Dicha función en definitiva posibilita sincronizar variables de salida de los convertidores
[3] y calcular el flujo de potencia reactiva y activa [4].
Existen diferentes enfoques para la detección del módulo y la fase de la componente de la
secuencia positiva en la tensión de red. La elección del método empleado dependerá de la
capacidad y precisión a la hora de determinar caracterı́sticas tales como el ángulo de fase,
la frecuencia y la amplitud, en situaciones donde la red sufre desequilibrios por cargas,
fallas o desbalanceos.
Soluciones tales como obtener el módulo y fase de las mismas mediante detectores de valor
de pico y de paso por cero [5], son insuficientes debido a que solo se detecta la tensión de la
red una vez cada medio ciclo, provocando un retardo de respuesta del sistema, además que
la presencia de armónicos o los desequilibrios de las tensiones de red generan mediciones
erróneas. Dicho método funcionarı́a idealmente si fuera un sistema con un conjunto de
sinusoides perfectas y equilibradas.
Estas circunstancias convierten como solución ideal el empleo de un lazo de enganche de
fase (PLL- Phase Locked Loop), el cual podrı́a definirse como un sistema digital compuesto
por tres elementos principales [6], como se evidencia en la representación esquemática de
la Figura D.1:
Detector de fase (PD). Se establece la diferencia entre la referencia y la salida.
Bucle de filtrado (LF).
Oscilador controlado digitalmente con una señal ya sincronizada con la entrada.
120
Figura D.1. Modelo esquemático de un PLL [6]
Generalmente, las distintas estructuras disponibles del PLL se diferencian mayormente por
la implementación del PD. Los factores más relevantes para alcanzar el comportamiento
u operación correcta de la estructura PLL son [7]:
La estabilidad del PLL, buscando una estructura robusta capaz de ofrecer
respuestas idóneas incluso en las situaciones más adversas. Una forma de lograrlo
es estableciendo la distancia correcta entre los polos del bucle cerrado.
El comportamiento del PLL ante componentes con distorsiones, logrando la
capacidad de rechazo adecuada ante las perturbaciones. Es importante el destino
de la aplicación, puesto que sirve como referencia para conocer aproximadamente la
distorsión presente en vectores como el de la tensión.
El SRF-PLL conduce a buenos resultados en términos generales, sin embargo para sistemas
donde se debe hacer frente armónicos relevantes en la tensión de red, dicha técnica puede
presentar un error permanente. Una alternativa empleada para solucionar dicha situación
es la de separar las secuencias positiva y negativa de las tensiones y retroalimentar solo
la positiva [8]. Otra técnica utilizadas para hacer frente a armónicos de tensión de orden
elevado es la de reducir el ancho de banda del SRF-PLL para lograr atenuar dichos
armónicos [2], pero dicha solución presenta ciertos inconvenientes que se traducen en una
actuación deficiente con un error constante, cuando se produce desequilibrios o transitorios
en las tensiones de frecuencia fundamental [8] [1].
Además de que el comportamiento de los PLLs lineales con estructuras más sencillas en
sistemas de baja frecuencia como redes débiles o microredes, donde existen impedancias de
distinta ı́ndole o cargas no lineales, no el más efectivo o adecuado, debido a que el propio
PLL adquiere un carácter no lineal. El control de los sistemas bajo estas circunstancias
ofrece una respuesta errónea, como se demuestra en diversos estudios [8] [9].
En el caso de necesitarse un compartimiento óptimo bajo dichas circunstancias, donde el
SRF-PLL presenta un error constante, se puede emplear una topologı́a que aumente la
exactitud de detección de fase, un ejemplo es el DSRF-PLL [2].
121
Bibliografı́a
122