Análisis y simulación de VSC’s

con modulación PWM para uso

en sistemas HVDC
Rafael Antonio RESUMEN 1. INTRODUCCIÓN
Peña Suesca1
Un sistema de Transmisión DC a Alta Ten- La importancia de la transmisión DC, se
César Leonardo sión (HVDC) requiere de conversores electró- debe entre varias razones, a los problemas de
Trujillo Rodríguez2 nicos por su característica de convertir energía inestabilidad que presentan los sistemas AC
eléctrica AC en DC y viceversa. Existen bási- cuando se transporta energía en distancias lar-
camente dos tipos de configuraciones posibles gas, que obligan a centros de investigación y a
de conversores trifásicos para este proceso de empresas multinacionales a buscar mecanismos
conversión: el Conversor de Fuente de Co- que permitan solucionarlos. Una solución es la
rriente (CSC) y el Conversor de Fuente de transmisión DC a alto voltaje o HVDC (High
Voltaje (VSC). Voltage Direct Current) y a nivel mundial las
En el presente artículo se analizarán los interconexiones HVDC han permitido solven-
conversores trifásicos VSC con modulación tar problemas de inestabilidad y adicionalmente
PWM, los balances de potencia activa y reactiva, interconectar de manera sencilla sistemas don-
la interconexión asincrónica de sistemas trifásicos de la tensión y la frecuencia no son compati-
con enlaces VSC y el dimensionamiento de al- bles (caso Brasil – Argentina). Los sistemas
gunos elementos de la parte de potencia, utili- HVDC se han difundido y son cada vez más
zando el SimPowerSystems de MatLab®. comunes por ejemplo, el proyecto Cabora
Bassa inaugurado en 1978 que une a Apollo
Palabras clave: Conversor de fuente de con Mozambique Sur África en extensión de
voltaje VSC, PWM, sistemas HVDC. 1360Km a través de un circuito bipolar de
533kV y potencia de 1920 MW [1], así como
Analysis and simulation of VSC with PWM el proyecto comisionado para el 2010-2020 que
modulation for use in HVDC systems une a Jinshajiang con Shanghai China en exten-
sión de 2000Km a una tensión de 800kV y
ABSTRACT potencia de 24.8 GW [2].
A High Voltage Direct Current Transmission
system (HVDC) requires of electronic Hoy por hoy en Colombia el Plan de Ex-
converters by its characteristic to turn electrical pansión 2004 - 2018 de la UPME [3], tiene
AC energy in DC energy or vice versa . Basically previsto la interconexión eléctrica entre Colom-
exist two types of possible configurations of bia y Panamá, con posibilidades de interco-
three-phase converters for this process of nexión HVDC, y es por este motivo que se
conversion: Current Source Converter (CSC) requiere investigar en todas las etapas relacio-
and Voltaje Source Converter (VSC). nadas con el proceso de transmisión DC, sien-
do una etapa de gran importancia, la etapa de
In the present paper will be analyzed three- conversión de energía.
phase converters VSC with modulation PWM,
the balance of active power and reactive En este artículo se analiza la topología
power, the asynchronous interconnection of conversora de fuente de voltaje o VSC (Vol-
1
Profesor de la Facultad de
three-phase systems with VSC converters and taje Source Converter) a PWM para siste-
Ingeniería de la Universidad
Distrital.
the sizing of some elements of the part of mas de transmisión y conversión HVDC y
2
power, using the SimPowerSystems of se dan algunos elementos para diseñar la sa-
Profesor de la Facultad de
Ingeniería de la Universidad MatLab®. lida de la etapa de potencia de una unidad
Distrital.

Key words: HVDC systems, PWM, voltage conversora VSC, para subestaciones de sis
Artículo recibido en Junio de
2007, aprobado para publicación source converter VSC. temas de potencia. 30
en Agosto de 2007

30 Vol.11 No.2
REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS
Ingeniería
 2. MODELO BÁSICO DE riando en el tiempo para tener una función
CONVERSOR TRIFÁSICO DC/AC Van(t) variable en el tiempo y controlada por
δ(t). La frecuencia de δ(t) se conoce como fre-
2.1. Modo Inversor DC/AC cuencia de modulación .
El modo inversor permite convertir una Despejando δ(t) de (1) se tiene:
tensión DC en una tensión AC de amplitud y
frecuencia variable. En la figura 1 se muestra 1 1 V an (t )
δ (t ) = + (2)
la etapa básica de potencia de un conversor 2 2 E
PWM trifásico en modo inversor con carga
resistiva y en la figura 2 el modelo simulink Para obtener el valor deseado V an (t) =
correspondiente. Emsen(ωt+ϕ), donde Em es el valor pico de
tensión, ω es la frecuencia angular y ϕ es el
ángulo de desfasamiento, se requiere,

1 1 Em 1 1
δ (t ) =
+ sen(ωt + ϕ ) = + m ∗ sen(ω mt + ϕ ) (3)
2 2 E 2 2
Em
∴ m= y Em ≤ E
E

Para el caso trifásico es indispensable que los

voltajes Van(t), Vbn(t) y Vcn(t), estén desfasados
2π/3, por lo tanto:
In Mean

+
i
- iE1(t)
Val or Medio IE1

Puente Trifásico
IE1
δA(t) = ½ + ½ m*sen(ωt),
δB(t) = ½ + ½ m*sen(ωt-2π/3) y
125 kV Inversor SPWM
g Salida
+
i signal rms
A +
-

125 kV1
-
B

C
LB
LA
ILA RMS
ILA

δC(t) = ½ + ½ m*sen(ωt+2π/3),
i
+ iE2(t)
- LC

100MWa +
que establecen los voltajes de fase:
100MWb v signal rms
-

100MWc Van
Van RM S

IE1 Van(t) = E*m*sen(ωt),

Mediciones
Vbn(t) = E*m*sen(ωt-2π/3) y
Figura 2. Modelo simulink inversor
i signal rms
+
-

LB
trifásico DC/AC básico.
LA
ILA RMS
ILA
Vcn(t) = E*m*sen(ωt+2π/3),
Para las figuras 1 y 2, se tienen los siguientes
100MWc
Para obtener una tensión sinusoidal de una
100MWb
100MWa +
-
v signal rms

Van valores, E1 = E2 = 125kV, S = 300MW, fm =

Van RM S

frecuencia determinada, con un reducido nú- 60Hz, fc = 1800Hz y m = 0.9.

mero de armónicos, se emplea la técnica cono-
cida como modulación por ancho de pulso Los datos de 150kVDC por polo (300kVDC
PWM (Pulse Width Modulation) [4]. Para ello polo a polo) y potencias de 550MW corres-
se compara una señal de control sinusoidal, de ponden a valores máximos aceptados en la ac-
la misma frecuencia que la tensión de salida de- tualidad para enlaces DC con conversores VSC,
seada, con una forma de onda triangular de de acuerdo a la bibliografía consultada [6, 7, 8,
mayor frecuencia. La frecuencia de la onda trian- 9]. Los valores de 125kVDC por polo, 300MW
gular, establece la frecuencia de conmutación del y longitud total de 571km utilizados en la si-
inversor (interruptores del inversor) y se man- mulación, corresponden a valores tomados del
tiene constante. Plan de Expansión 2004 - 2018 de la UPME,
el cual tiene previsto la interconexión eléctrica
La expresión que gobierna el voltaje de sali- entre Colombia y Panamá, con posibilidades
da para la fase a del inversor está dada por [5] de interconexión HVDC.

Van (t ) = E[2δ (t ) − 1] (1) 2.1.1. Criterio para selección de las

inductancias de conmutación LA, LB y LC.
Siendo E: el voltaje DC y δ: la relación del
Considerando los valores anteriores como va-
ancho del pulso T1 (tiempo de encendido de
lores base, para el lado AC se puede establecer:
los interruptores) con respecto al periodo T,
ZBASE = (V BASE) 2/S BASE = (125kV/√2) 2/
31 δ(t) significa que el ancho del pulso está va-
100MW = 78Ω.

REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS Vol.11 No.2 31
Ingeniería
 Teniendo en cuenta que, las inductancias LA, Para los datos presentados anteriormente, los
LB y LC deben ser lo suficientemente grandes valores fueron:
para disminuir los ΔI y bloquear las altas fre- Tabla I. Valores de voltaje y corriente
cuencias de conmutación, pero a su vez que a soportado por los dispositivos
S3φ =
60Hz sus reactancias sean tan pequeñas que IE1DC VcnRMS ICRMS
P E1&E2 =
3 VcnRMS*
QL =
(E1+E2)*IE1DC ILCRMS2*XLC
ILCRMS
las caídas de tensión a esta frecuencia no sean 1159A 78kV 1235A 289.75MW 288.99MW 23.80MVAr
considerables, entonces se escogen valores
XLA, XLB y XLC que estén entre 0.15 a 0.20 en El contenido armónico de los voltajes VAn(t)
por unidad (p.u.). y Van(t) se presenta en la figura 4. La compara-
Para XLA = XLB = XLC = 0.2 p.u. en valores ción de los dos espectros de frecuencia mues-
reales X LA = 0.2*78 = 15.625Ω y L A = tra que el THD pasa de 108.25% a 13.34%
(15.625)/(377) ≈ 40mH. La potencia reactiva por efecto de la inductancia de conmutación.
Q nominal a 60Hz de la reactancia de la fase Sin embargo se hace necesario seguir reducien-
A es QX LA = I A2*X LA = (S A/V A) 2* X LA = do el conjunto de armónicos más significati-
20MVAr. Teniendo en cuenta los reactivos vos presentes alrededor de 1800Hz, 3600Hz y
de armónicos superiores, entonces QXLA = 5400Hz, para disminuir el THD y mejorar la
QX LB = QX LC > 20MVAr. calidad de las formas de onda de los voltajes
de salida. Dicha reducción de armónicos se
Las formas de onda de VAB(t), VAn(t), Van(t), puede lograr mediante la implementación de
Vbn(t) y Vcn(t), correspondientes a la simula- filtros sintonizados.
ción, se muestran en la figura 3.
2.1.2. Criterios para el Cálculo
de Filtros en el lado AC
• Colocar en paralelo a la salida filtros serie
sintonizados [10] que sean corto circuito a
las frecuencias de los armónicos a eliminar.
• Dimensionar los capacitores de los filtros
para que en la frecuencia de 60Hz com-
pensen los reactivos de la inductancia de
conmutación.

Figura 3. Formas de onda de V AB(t), V An(t), Van(t), Vbn(t) y Vcn(t).

En la figura 3 se pueden apreciar los trenes

de pulsos de los voltajes de salida antes de las
inductancias de conmutación, VAB(t) y VAn(t), y
los voltajes de salida Van(t), Vbn(t) y Vcn(t) que
presentan rizado apreciable debido al conteni-
do armónico de orden superior que no alcan-
zó a ser rechazado en su totalidad por las
inductancias de conmutación.
En la tabla I se presentan los valores de co-
rriente y voltaje obtenidos de la simulación
IE1DC, VcnRMS e ICRMS y los valores calculados
de potencia P E1&2, S3ϕy QL. Los resultados son
inferiores a los valores nominales porque la
potencia nominal de 300MW se cumple para
voltaje nominal 88.39kV que se obtiene con m
32
= 1 y caída despreciable en XLc. Figura 4. Espectro de frecuencia de (a) VAn(t) y (b) Van(t).

32 Vol.11 No.2
REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS
Ingeniería
 Para el diseño de los filtros serie sintoniza- Para V = 88kV en la tabla II se presentan los
dos se debe tener en cuenta que en la frecuen- valores de L y C de los tres filtros sintonizados.
cia ù1, la capacitancia C presenta reactancia
XC(ω1) = V2/Q y en ω2, presenta reactancia En la figura 6 se muestra que el rizado en las
XC(ω2) = XC(ω1)*ω1/ω2 y para que se dé la re- formas de onda de los voltajes de salida dismi-
sonancia en ω2, XL(ω2) = XC(ω2), es decir ω2L = nuyó considerablemente presentando THD de
XC(ω1)*ω1/ω2, que permiten obtener: 0.46% al colocar los tres filtros serie sintoniza-
dos, en comparación con el rizado y el THD
V 2 ω1 Q presentado en las figuras 3 y 4 respectivamente.
L = y C = (4 )
Q ω 22 ω 1V 2

donde Q son los reactivos a compensar, ω1

es la frecuencia angular a 60Hz y ω2 la frecuen-
cia del armónico que se quiere eliminar.
Para eliminar las frecuencias de 1800Hz,
3600Hz, 5400Hz y sus laterales correspondien-
tes en los voltajes de salida, se pueden diseñar
filtros LC de factor de calidad Q alto como los
que se muestran en las salidas de voltaje de la Figura 6. Espectro de frecuencia de Van(t)
figura 5 y en los que se han depreciado las com-
Teniendo en cuenta que los filtros son corto
ponentes resistivas. Es de anotar también que
circuito a las frecuencias de resonancia y par-
para atenuar los armónicos de orden superior
tiendo del valor del voltaje de los armónicos a
más alto se recomienda utilizar filtros con Q compensar calculados de la figura 6, 62.09kV,
bajo y mayor ancho de banda. Es importante 22.28kV y 14.33kV para 1800Hz, 3600Hz y
realizar estudios que establezcan los tipos de fil- 5400Hz respectivamente y de los valores de la
tros requeridos en una determinada instalación reactancia de conmutación a estas frecuencias,
dependiendo de la distorsión armónica y la in- las corrientes que circulan por los filtros son
terferencia generada permitidas. 137.2A, 24.62A y 10.55A, y las potencias en las
reactancias 606.7KVAr, 9.77KVAr y 1.19KVAr
In Mean

IE1
respectivamente a las frecuencias de resonancia.
i
Valor Medio IE1

125 kV
+
- iE1(t)
Puente Trifásico
Inversor SPWM
VAn(t) +
-
v A 60Hz la reactancia de los condensadores de
g Salida VAn

125 kV1
+
A

B
+
i
-
signal rms

ILA RMS
ILA
los filtros determinan la corriente obteniéndose
iE1(t) & i E2(t)

7.5KVAr, 1.82KVAr y 0.84KVAr respectivamen-

-
C LA
LB
i
+ iE2(t )
-
LC

Discrete,
Ts = 1e-006 s.
te en las reactancias de los filtros.
100MWa +
LCa v signal rms
100MWb -
LCb
Ts = 1e-6
100MWc LCc Van
Van RMS

Para el caso de cargas con generación es con-

Mediciones
veniente mencionar aspectos fundamentales del
flujo de potencia en líneas de transmisión de
Figura 5. Modelo simulink inversor alto voltaje, que por su carácter inductivo, tie-
trifásico DC/AC con filtros LC serie sintonizados. nen flujos de potencia activa que dependen de
las diferencias angulares y flujos de potencias
De la tabla I los reactivos Q a compensar son reactivas que dependen de las diferencias de
de 23.8MVAr, se toma Q ≈24MVAr que repar- tensión, dadas por las expresiones[11]:
tidos en los tres condensadores de los filtros
sintonizados a 1800Hz, 3600Hz y 5400Hz, co- VAVB V
P= sen(ϕA −ϕB) y Q = A (VA −VB cos(ϕA −ϕB)) (5)
rresponden a una compensación de 8MVAr a X X
60Hz por filtro para V = 125kV/√2.
En la figura 7 se presenta el modelo simulink
Tabla II. Modelo simulink inversor con carga activa. Las reactancias de conmuta-
trifásico DC/AC con filtros LC serie sintonizados.
ción en la etapa conversora PWM y las
Frecuencia (Hz) L (mH) C (µF) reactancias del equivalente Thevenin de la car-
1800 2.85 2.74 ga con generación o del sistema AC que recibe
3600 0.713 2.74 la potencia hacen que el flujo de potencia esté
33 5400 0.317 2.74 modelado por las expresiones anteriores.

REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS Vol.11 No.2 33
Ingeniería
 Considerando las reactancias del equivalente 2.2. Modo Rectificador AC/DC
Thevenin del orden del 0.15 % en p.u. para los Para valores de ϕ< 0o el conversor trabaja
valores trabajados a lo largo del artículo, las en modo rectificador. Para analizar este caso
reactancias resultan del orden de 11.7Ω y las en la figura 8 se muestra la etapa básica de
inductancias de 30mH aproximadamente. potencia de un conversor trifásico en modo
rectificador con carga resistiva de 100Ωy en la
Tomando el voltaje Thevenin del sistema AC
figura 9 el modelo simulink correspondiente.
como referencia para medir el desfasamiento
ϕ del conversor, se pueden expresar los valo-
res de δ del conversor como se indica a conti-
nuación:
δA(t) = ½ + ½ m*sen(ωmt+ϕ)
δB(t) = ½ + ½ m*sen(ωmt-2π/3+ϕ)
δC(t) = ½ + ½ m*sen(ωmt+2π/3+ϕ)
Figura 8. Modelo básico de rectificador trifásico AC/DC.
Para el circuito de la figura 7 fue necesario
agregar filtros LR paralelos en serie con cada +
i
-
In Mean

I Carga p DC IDC
SPWM

fase, para estabilizar el sistema. Con valor de Salida g

A
+
Carga p
+
-
v In Mean

a V Carga p DC VDC

voltaje en el lado DC de 125kVDC por polo y LTA

LA

b LB
B

C
-

valor de voltaje pico de 100kV en el genera- LTB

LTC
LC Puente Tri fásico
Rectificador
Carga n

dor AC trifásico, se hicieron las siguientes prue-

bas variando ϕ y m.
V1(t) V2(t) V3(t)

V
PQ
I Di screte,
PQ Ts = 1e-006 s.
P&Q
In Mean s ignal rms

i Valor Medio IE1

IE1 ILA T s = 1e-6
+ IE1(t) ILA RMS
-
+
+ v s ignal rms
Puente Tri fásico v -
SPWM -
Inversor Inversor
VAn Van
125 kV Van RMS
g Salida
+
i
A A A +
- a Medi ci ones
B LA B B
b
-
i E1(t) & iE2(t) LB

Figura 9. Modelo simulink rectificador

C C C c
125 kV1
LC

trifásico AC/DC básico.

LTC LT B LTA
i
+ IE2(t)
-

F3 F2 F1
+
v In Mean V3(t) V2(t) V1(t)
-

E2
E2 DC

Di screte,
Ts = 7.5e-006 s. Con valores pico de 125kV para V1(t), V2(t)
y V3(t) del sistema AC, m = 0.8 y ϕ = -10o para el
Medi ciones
Inversor

Figura 7. Modelo simulink inversor conversor, en la figura 10 se muestran las for-

trifásico DC/AC con carga activa
mas de onda de voltaje y corriente del lado
La tabla III presenta que en ϕ = 0o el flujo DC y en la figura 11 el espectro de frecuencia
de potencia activa a la carga es poco significa- de la señal de voltaje sobre las resistencias de
tivo y en ϕ = 20o existe flujo de potencia acti- carga en el lado DC.
va. Con m = 0.8 el voltaje teórico de salida del
conversor tiene valor pico Vpico = m*VDC =
0.8*125kV = 100kV y el flujo de potencia
reactiva es poco significativo. Con m < 0.8 el
flujo de potencia reactiva tiende a ser negativo
y con m > 0.8 el flujo de potencia reactiva
tiende a ser positivo.

Tabla III. Valores de P y Q para diferentes

valores de ϕ y m en modo inversor.
φ[grados] m P[MW] Q[MVAr]
0.7 -7.2 -34.2 Figura 10. Formas de onda de V carga(t) e I carga(t).
0 0.8 0.3 23.4
0.9 12.9 81.3 El espectro de frecuencias muestra presen-
0.7 124.2 -64.5 cia importante de armónicos en 180Hz, alre-
20 0.8 152.4 -13.74 dedor de 1800Hz, 3600Hz y otras frecuencias
0.9 180 44.1 de orden superior. 34

34 Vol.11 No.2
REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS
Ingeniería
 un rizado reducido. La figura 13 muestra dis-
minución del THD (7.96%) y se ve que es con-
veniente instalar un filtro adicional para elimi-
nar los armónicos cercanos a 1800Hz, el cual
se presenta en la figura 14.
Lp
+ i In Mean
-

I Carga p DC IDC
SPWM
Salida g Cp Cp 1 +
Carga p - v In Mean
+
A

Figura 11. Formas de onda de V carga(t) e I carga(t). a V Carga p DC VDC

LA B

b -
LT A LB C
Cn 1 Carga n
c LC Puente Tri fá si co Cn
LT B
Rectificador
Ln
LTC

2.2.1. Criterio para el Cálculo V1(t) V2(t) V3(t)

de Filtros en el lado DC
El filtro en el lado DC debe permitir el paso Discrete,
T s = 1e-006 s.

de la componente DC y cortocircuitar las com- T s = 1e -6

ponentes de alta frecuencia. Medicione s

Tomando por polo 125kVDC y 150MW Figura 14. Modelo simulink rectificador
trifásico AC/DC con filtros DC.
representados por R = V2/P = 100Ω, se pue-
den filtrar frecuencias cercanas y mayores a Considerando que se quieren reducir los ar-
1800Hz, instalando en cada polo de la salida mónicos de 1800Hz y superiores, con un filtro
DC del rectificador un condensador C que en L, se escoge XCp1 << XLp. Escogiendo el
cumpla RC >> TC. Entonces C >> (1/ condensador del filtro adicional de 20μF se tie-
(1800*100)) y se escoge C ≈ 20μF. ne XCp1 = 4.42Ω a 1800Hz y escogiendo XLp =
Con C incluido, en la figura 12 se muestra el 10*XCp1 se tiene Lp ≈ 4mH. En la figura 15 se
voltaje y la corriente del lado DC, 50kV y muestra que con la adición del nuevo filtro, el
500A respectivamente, y en la figura 13 su es- voltaje de salida en estado estable presenta una
pectro de frecuencia. distorsión armónica

Figura 15. Espectro de frecuencia de Vcarga(t)

con filtros adicionales en el lado DC.
Figura 12. Formas de onda de V carga(t) e Icarga(t) AC/DC.
Por otra parte, para mostrar el efecto de ϕ
sobre el flujo de potencia, en la figura 16 se
muestra la salida con ϕ = 0o y en la figura 17
se muestra la salida con ϕ = -20o

Figura 13. Espectro de frecuencia de Vcarga(t).

Se observa en las figura 12 que debido a la

35 presencia del condensador, el voltaje DC tiene Figura 16. Forma de onda de V carga(t) con ϕ = 0o.

REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS Vol.11 No.2 35
Ingeniería
 La figura 16 muestra que con ϕ = 0o el vol- SPWM Puente T rifásico
Lpfr

Lpr
571-km cabl e p

Lpi
Lpfi
+
i
- ip(t)
In Mean

Val or Medio Ip Ip

Puente Trifásico
+ v
-

VAn
Rectificador Rectificador SPWM

taje DC de salida es muy bajo y por tanto la

Inversor Inversor signal rms
Salida g Cpr Cpfi Cpi
Cpfr g Salida
+ +
A ILA1
i ILA1 RMS
a A + - a
LA B
b B LA1 b
-
LTa LB C ip(t) & in(t) -
c C LB1
LTb Cni c
LC Cnfr
Cnr Cnfi
LC1 +
LT c v signal rms

transferencia de potencia no es significativa. La

-
Lnfr 571-km cabl e n Lnfi
Carga_r i
+ - in(t) Van
Fa Fb Fc Van RMS
Lnr Lni
F3 F2 F1 Carga_i

V1(t) V2(t) V3(t)

+
- v In Mean

figura 17 muestra que con ϕ = -20o el voltaje

VCni DC VDC

Discrete,
T s = 1e-006 s.

Ts = 1e-6

DC de salida es de 97.1kV por polo y la po- Mediciones

Rectificador
Mediciones
Inversor

tencia de 188MW. Figura 18. Modelo simulink de sistema básico de

transmisión HVDC con conversores VSC.

En la figura19 se muestra el voltaje DC y la

corriente DC de salida de la primera estación.

Figura 17. Forma de onda de V carga(t) con ϕ = 20o.

En la tabla IV se muestran los resultados de la

simulación del conversor en modo rectificador
con voltajes pico de 125kV en los generadores
trifásicos, R = 200Ω en el lado DC y para dos Figura 19. Formas de onda de VDCr e IDCr
entregados por el rectificador.
valores de desfasamiento y de modulación.
Tabla IV. Valores de P y Q para diferentes
valores de ϕ y m en modo rectificador.
Los valores en estado estable de voltaje y
corriente que se muestran en la figura 19 son:
φ[grados] m P[MW] Q[MVAr] 120kV y 770A respectivamente por polo. La
0 0.95 66.3 310.8
0.95 268.5 -208.5
potencia transmitida es 184.8MW.
-20
0.75 225.3 28.5
La segunda estación funciona en modo in-
versor con los siguientes datos: Sistema DC de
De la tabla anterior se puede concluir que
alimentación de aproximadamente 120kV por polo.
para valores de ϕ negativos (-20o) el flujo de
Conversor Trifásico funcionando en Modo Inversor con
potencia hacia el conversor aumenta y con cam-
m = 0.9, ϕ = 0o (puede ser cualquier ángulo por no
bios del índice de modulación m de 0.95 a
tener generación a la salida) y fC = 1800Hz.
0.75 en ϕ = -20o la potencia reactiva cambia
en magnitud y pasa de adelanto a atraso. En la figura 20 se muestra el voltaje VABi(t) y
los voltajes Vani(t), Vbni(t) y Vcni(t) sobre la car-
3. MODELO BÁSICO DE TRANSMISIÓN ga en el inversor.
HVDC CON CONVERSORES VSC
En la figura 18 se muestran dos sistemas
trifásicos interconectados a través de un siste-
ma HVDC bipolar de 125kV, 300MW,
implementado con conversores VSC. Una de
las estaciones es alimentada desde un sistema
AC con generación, interconectada a través de
una línea de transmisión de 571km, con otra
estación con carga y sin generación. Los datos
de la línea son: 1.39*10-2Ω/km, 1.59*10-4H/
km y 2.31*10-7F/km. Figura 20. Formas de onda de
<VABi(t), Vani(t), Vbni(t) y Vcni(t) en el inversor.
La primera estación funciona en modo
rectificador con los siguientes datos: Sistema AC En el inversor los valores de voltaje son
de alimentación trifásico de 88kV y 100MW por fase 77.57kV por fase, entrega un potencia total e
a 60Hz.Conversor Trifásico funcionando en Modo 180.51MW, que son aproximadamente los
Rectificador con m = 0.9, ϕ = -20o y fC = 1800Hz. mismos recibidos desde el rectificador. 36

36 Vol.11 No.2
REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS
Ingeniería
 En la figura 21 se muestra un enlace VSC de + v
-
s ign al rm s

Va n
V

I
PQ

PQ
In Me an

Valor Medio Ipr

Ip r
In Mean
In Mean
V

I
PQ

PQ1
Van RMS P &Q + VAA(t)
Ipri P & Q1

dos sistemas AC con generación. Los sistemas

- v i Lpfr Valor Medio Ipri
+ i Lpfi Ipi
- ipr(t) + ipri (t) i Valor M edio Ipi + v
- + ip i(t) s ignal rm s
VAn - -
Lpr 571-km cable p Lpi
s ign al rm s + v
SPWM Puente T rifásico + Vaa RMS VAArm s
s ignal rm s
Rectificador Puente Trifásic o SPWM - v s ignal rm s -
Rectificador
Inve rsor Inversor
ILA Sa li da g Cp fi VAn1 ILA1 Van1
ILA RM S Cpr Cpi Van1 RM S
Cpfr g Sa lid a ILA1 RM S
+ i + +
- A A A i
A A A +

se representan con sus equivalentes Thevenin

-
B B LA B
B LA1 B B
-
C C LB C -
ip i(t) & i ni(t) L B1
C C C
LC Cnfr Cni
Cnr Cn fi LC1

Lnfr 571-km c able n Lnfi

+ i
- i ni(t)
LTa LT b LT c Fa Fb Fc

que para este caso son iguales.

Lnr Lni
F3
F2 F1
LT3 LT 2 LT1
+ v In Mean + v In Mean
- -

V1(t) V2(t) V3(t) VDCr VDCi

VCnr DC VCni DC
V3n(t) V2n(t) V1n(t)

Discrete,
Ts = 7.5e -006 s.

Debido a que no era parte del desarrollo del M edi ciones

Rec tificado r
Med iciones
In versor

documento incluir lazos de control que po- Figura 21. Modelo simulink de sistema de transmisión HVDC
drían haber sido para mantener el voltaje DC con conversores VSC y generación a ambos lados.

constante, la corriente DC constante o el vol- con ϕ = -20o y m = 0.98. Absorbe una poten-
taje AC constante, se incluyeron filtros serie de cia activa de 246.6MW que distribuye,
resistencia en paralelo con inductancia para 125.46MW a la estación inversora de 60Hz y
poder estabilizar el sistema. El criterio fue se- 95.88MW a la estación inversora de 50Hz. La
leccionar resistencias lo más bajas posibles para diferencia de potencia corresponde a las pér-
no tener pérdidas altas y en paralelo con didas presentes en los filtros serie, en la esta-
inductancias lo más bajas posibles para reducir ción rectificadora y en las líneas DC. También
aun más la corriente de las resistencias; pero se observa que maneja una potencia reactiva
los valores que estabilizaron el sistema no fue- en adelanto de 60.8MVAr.
ron tan bajos como se requería.
Las estaciones inversoras de 60Hz y 50Hz
Con valores de voltajes de generación de transforman la potencia DC recibida de la esta-
125kVpico y frecuencias de 60Hz se realizaron ción rectificadora y entregan potencia AC a sus
simulaciones que entregaron los datos que se correspondientes sistemas de 123.6MW y
presentan en la tabla 1 del Anexo 1. En esta 85.5MW respectivamente. La estación inversora
se observa: de 60Hz entrega potencia reactiva en atraso de
13.5MVAr y la de 50Hz entrega potencia reactiva
• El flujo de potencia reactiva en los
en adelanto de 49.5MVAr.
conversores es menor para índices de mo-
dulación cercanos a la unidad. Lo más importante de la simulación ante-
rior es evidenciar el hecho de que los siste-
• El incremento del desfasamiento en la es-
tación rectificadora aumentó el flujo de po- mas VSC pueden ser multiterminales es de-
tencia activa. cir a un mismo nodo DC se pueden conectar
diferentes estaciones con generación, con car-
• La disminución de dos grados en el ga o con generación y carga.
desfasamiento en la estación inversora au-
mentó un poco la potencia activa, contra- En este punto y considerando lo analizado a
riamente a lo que podría esperarse. La ra- lo largo del documento es importante obser-
zón fue el aumento del voltaje DC en la var que por su naturaleza los conversores VSC
línea que a su vez aumentó el voltaje AC a se comportan como fuentes de voltaje en el
la salida del inversor puesto que m perma- lado DC, mantienen la polaridad del voltaje
neció constante. Adicionalmente el flujo de DC independiente del sentido del flujo de
potencia activa se vio favorecido porque potencia y que por tanto varios conversores
se disminuyó el flujo de potencia reactiva. VSC del mismo nivel de voltaje DC se pue-
den conectar a un mismo barraje DC.
En la figura 22 se muestra el modelo
simulink de la interconexión de tres sistemas + v
-
s ign al rms

Van RMS
Van
V

P&Q
PQ

PQ
+ v
In Mean

Valor Medio Ipr

Ipr
In Mea n

Ipri
In Mean
VAA(t)
V

I
PQ

PQ1
Lpfr Valor M edio Ipri P & Q1

AC de 125kVpico con tres estaciones

- + i i Lpfi Ipi
- ipr(t) + i Valor Medio Ipi +
- ipri (t) + ipi(t) v signal rm s
- -
VAn Lpr 571-km c able p Lpi + v
s ign al rms SPWM Puente T rifásico + v Vaa RM S si gnal rms
Puente Trifásico VAArms sig nal rm s -
Rectificador Rectific ador SPWM -
Inversor Inversor
ILA Sa lid a g Cpfi Van1
ILA RMS Cpr Cpi VAn1 ILA1 Van1 RMS
Cpfr g Salida ILA1 RMS
i +
+ - A A A +
i
A A A + -

conversoras, una en modo rectificador a 60Hz

B B LA B
B LA1 B B
-
C C LB C -
ipi(t) & ini(t) LB1
C C C
LC Cnfr Cni
Cnr Cnfi LC1

Lnfr 571-km cable n Lnfi

i
+ - ini(t)
LT a LTb LT c Fa Fb Fc
Lnr Lni

y dos en modo inversor trabajando a 60Hz y

F3 F2 F1
LT 3 LT2 LT 1
+ +
- v In Mean - v In Mean

V1(t) V2(t) V3(t) VDCr VDCi

VCnr DC VCni DC
V3n(t) V2n(t) V1n(t)

Discrete,

50Hz respectivamente.
Ts = 7.5e-006 s.

V
PQ
In Mean I
In Mea n
VAA(t)1 PQ2
Ipri1 P &Q2
Mediciones M ediciones Valor Medio Ipri1
i Lpfi1 Ipi1
Rectificador Inversor + ipri(t) i Valor Medio Ipi1 +
- + ipi (t) v s ign al rms
- -
571-km cable p1 Lpi1 + v
+ v Vaa RMS1 signal rm s
Puente T rifásico VAArm s1 s ignal rms -
SPWM -
Inversor1 Inversor1
Cpfi1 Van2
Cpi1 VAn2 ILA2 Van1 RMS1
ILA1 RMS1

Se realizaron simulaciones sobre el modelo

g Sali da
+
i
A A A + -

B LA2 B B
ipi(t) & ini(t)1 -
C LB2 C C
Cni1
Cnfi1 LC2

simulink de la figura 22, y en la tabla 2 del

571-km cable n1 Lnfi1
+ i
- ini (t)
Lni1
F6 F5 F4
LT6 LT 5 LT4
+
- v In Mean

VCni DC1 VDCi1

V3n(t)1 V2n(t)1 V1n(t)1

Apéndice 1 se presentan los resultados de una

37
de las simulaciones. Allí se observan los resul- Figura 22. Modelo simulink de interconexión de
tados de la estación rectificadora trabajando tres sistemas AC a través de líneas DC.

REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS Vol.11 No.2 37
Ingeniería
 Adicionalmente con cada conversor VSC se tancia dada y una potencia transferida especifi-
puede producir voltaje AC controlando su ca, el nivel de tensión es el adecuado así como
amplitud, frecuencia y fase a través de δ. La si la opción DC es más viable que la AC se
anterior característica permite concluir que el establece a partir de análisis económicos, téc-
barraje DC puede ser el punto común de va- nicos y ambientales que no eran objeto del pre-
rios sistemas AC trabajando a diferentes fre- sente artículo y por eso no se abordaron. Es
cuencias e intercambiando potencia a través de de aclarar que una vez tomados los datos de
la conversión AC/DC y viceversa. voltaje, potencia y frecuencia de operación es
posible diseñar la etapa de potencia según lo
establecido a lo largo de este documento.
4. DATOS DE DISEÑO
DEL CONVERSOR VSC Los interruptores se implementan colocando
En las tablas V, VI y VII se condensan los varias válvulas en serie, cada una de ellas en arre-
datos que conciernen al diseño, obtenidos a lo glos de varios dispositivos (IGBT’s), dependien-
largo del artículo, con el fin de reunir los prin- do el número de los niveles nominales de volta-
cipales elementos de la etapa de potencia de je y corriente especificados por el fabricante.
un conversor VSC, de 300MW para operar en Los componentes requieren estar diseñados
una línea bipolar de ±125kVDC y alimentado para operar a las diferentes frecuencias de tra-
de un sistema trifásico AC de 125kVpico, 60Hz. bajo del conversor, por ejemplo las inductancias
La frecuencia de conmutación escogida es de de conmutación de 40mH de un conversor
1800Hz que facilita reconstruir adecuadamen- trifásico de 125kVpico por fase y potencia to-
te la señal sinusoidal y no es muy alta que pro- tal de 300MVA, deben poder conducir corrien-
duzca sobrepicos y dv/dt exagerados. tes de 1200A aproximadamente que a 60Hz
Tabla V. Valores de los filtros AC. producen potencia reactiva de 12002*377*0.04
Filtros AC H ≈ 22MVAr. El diseño de estas inductancias
Filtro Inductancia de Filtros sintonizados LC
serie conmutación debe tener en cuenta que los circuitos magnéti-
f 1800 3600 5400
R 2.74μF 2.74μF 2.74μF
cos operen apropiadamente en las frecuencias
13.79 (8MVAr (8MVAr (8MVAr de 1800Hz y armónicos superiores.
C
Ω a a a
40mH
LA 60Hz) 60Hz) 60Hz)
22MVAr
LB 0.713m Los condensadores de filtro en el lado DC
a 2.85mH 0.317mH
LC
60Hz H
L (606.7KV (1.19KVA
(9.77KV del conversor además de servir de filtros para
31.02 L Ar r
Ar
mH a a
1800Hz)
a
5400Hz) las frecuencias superiores están sometidos a
3600Hz)
componentes alternas producidas por la po-
Tabla VI. Valores de los filtros DC.
tencia reactiva que debe compensar el
Filtros DC
conversor.
Capacitivo Filtro L Inductancia de suavizamiento
Lpf 4mH Se requieren estudios más detallados sobre
Cp 20 µF Lp 8mH
Cpf 20 µF
filtros y en especial sobre los filtros serie de
inductancia en paralelo con resistencia en cada
Tabla VII. Dispositivos semiconductores.
Dispositivos semiconductores – IGBT’s
fase, que presentan consumo de potencia alto
Vmáximo Imáxima Vinterruptor No de dispositivos cuando se quiere estabilizar el sistema en valo-
Mayor a 360 por interruptor,
1700V 800A 3V considerando un 20% de margen res de tensión, corriente y potencia adecuados.
de sobretensión.

Cabe destacar que para este diseño se partió

de los datos de voltaje y potencia sugeridos
por la UPME para la interconexión entre Pa-
5. CONCLUSIONES
namá y Colombia en la alternativa DC. To- En el presente artículo se analizaron
mando estos datos se calcularon filtros AC y conversores trifásicos VSC con modulación
DC, así como número de dispositivos PWM funcionando en modo inversor y en
semiconductores relacionados con la elabora- modo rectificador y dichos resultados se va-
ción de la etapa electrónica de potencia. La lidaron a través del SimPowerSystems de
parte de diseño que contempla si para una dis- MatLab®.
38

38 Vol.11 No.2
REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS
Ingeniería
 Se evidenció como a través de m y ö en ä, es [4] J. M. Benavent, et al. «Electrónica de Potencia: Teoría y
Aplicaciones». Valencia: Alfaomega Marcombo, 2000. 235p.
posible manipular la amplitud y el desfasamiento
[5] D. Hart. «Electrónica de Potencia», Madrid: Prentice Hall,
de los voltajes con el fin de controlar de manera 2001, p. 472.
independiente la transferencia de potencia acti- [6] J. Gutiérrez, et al. «Transporte de Energía Eléctrica en
va y de potencia reactiva del sistema. Corriente Continua: estado actual y perspectivas», Auto-
mática e Instrumentación, España, Abril 2005, p. 72.
Otra característica adicional que se logro [7] B. Andersen, et al. «A new era in HVDC», Power Systems,
IEE Review, Marzo 2000.
constatar fue la posibilidad que ofrecen estos
[8] K. W. Kanngiesser, et al. «HVDC Systems and Their
conversores para realizar interconexiones Planning», Siemens. Publisher: EV HA 7T. Rev. 4, 618p.
asíncronas de sistemas AC, adicionalmente se [9] M. Bahrman, et al. «The ABCs of HVDC Transmission
mostró la capacidad que tienen para crecer Technologies», IEEE Power & Energy Magazine, Vol. 5,
modularmente con el fin de facilitar amplia- No. 2, Marzo 2007, p. 32.
ciones futuras. [10] J. Arrillaga. «High Voltage Direct Current Transmission»,
Institution of Engineering and Technology, Sep. 1998. 299p.
Se estableció una metodología que permitió [11] J. Grainger y W. Stevenson, «Análisis de sistemas de
potencia», México: McGraw-Hill, 1996, p. 740.
dimensionar dispositivos para implementar
prototipos de la etapa electrónica de potencia
de conversores VSC PWM. Rafael Antonio Peña Suesca
Nació en Bogotá, Colombia. Es Ingeniero Electrónico de la Uni-
Finalmente dicho trabajo permitió estable- versidad Distrital Francisco José de Caldas. Obtuvo su título de
Maestría en Ingeniería Eléctrica en la Universidad Nacional.
cer y consolidar las bases para seguir con la Actualmente se desempeña como Decano de la Facultad de
siguiente etapa de desarrollo del proyecto re- Ingeniería en la Universidad Distrital Francisco José de Caldas
y es director e investigador del grupo LIFAE donde realiza estu-
lacionada con el control de la potencia activa y dios sobre sistemas de transmisión DC y fuentes no convencio-
reactiva del sistema por métodos de compu- nales de energía. [email protected]
tación flexible y control lineal.
César Leonardo Trujillo Rodríguez
REFERENCIAS BIBLIOGRÁFICAS Nació en Bogotá, Colombia. Es Ingeniero Electrónico de la
Universidad Distrital Francisco José de Caldas. Obtuvo su
[1] V. Sood. «HVDC and FACTS Controllers», Boston, Kluwer título de Maestría en Ingeniería Eléctrica en la Universidad
Academic Publishers, 2004. 255p. Nacional de Bogotá, Colombia. Actualmente se desempeña
[2] M. Szechtman, et al. «800-kV HVDC on the Horizon», como profesor en Electrónica de Potencia en la Universidad
IEEE Power & Energy Magazine, Vol. 5, No. 2, Marzo Distrital Francisco José de Caldas y pertenece como investi-
2007, p. 61. gador al grupo LIFAE donde realiza estudios sobre calidad de
potencia eléctrica, sistemas de transmisión DC y fuentes no
[3] «Plan de Expansión Preliminar 2004-2018», Unidad de convencionales de energía. [email protected]
Planeación Minero Energética, Colombia, Ago. 2004.

ANEXO 1
Tabla I. Valores de P y Q para diferentes valores de ö y m en el enlace HVDC.
Estación rectificadora Línea DC Estación Inversora
φ[grados] m P[MW] Q[MVAr] VDC[kV] PDC[MW] φ[grados] m P[MW] Q[MVAr]
-17 0.8 166.74 214.5 87.15 139.96 20 0.95 117.9 -217.65
-20 0.98 232.5 -4.5 113.6 210.16 20 0.98 196.5 -93
-20 0.98 249.3 -73.5 126.4 223.98 18 0.98 209.4 -20.1

ANEXO 2
Tabla 2. Valores de P y Q en el enlace HVDC que interconecta tres sistemas AC.
Línea DC Línea DC
Estación rectificadora (60Hz) Estación Inversora (60Hz) Estación Inversora 1 (50Hz)
Superior Inferior
φ P Q VDC PDC φ P Q VDC PDC φ P Q
[kV [M [kV [M
[grad m [M [MV ] W] [grad m [M [MV ] W] [grad m [M [MV
os] W] Ar] os] W] Ar] os] W] Ar]
-20 0. 246. -60.8 127 125. 10 0. 123. 13.5 127 95. 8 0. 85.5 -49.5
39 98 6 .5 46 98 6 .5 88 9

REVISTA CIENTÍFICA Y TECNOLÓGICA DE LA FACULTAD DE INGENIERÍA, UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS Vol.11 No.2 39
Ingeniería