Control para eliminación completa de tensión de secuencia negativa en una microred eléctrica basada en varios convertidores de potencia

TRABAJO FINAL DE MÁSTER
TÍTULO: Con ol pa a eliminación comple a de ensión de secuencia
nega i a en una mic o ed eléc ica basada en a ios con e ido es
de po encia
AUTOR: BREÑA DIAZ, LUIS ALBERTO
FECHA DE PRESENTACIÓN: 27/05/2025
APELLIDOS: BREÑA DIAZ NOMBRES: LUIS ALBERTO
TITULACIÓN: Con ol pa a eliminación comple a de ensión de secuencia
nega i a en una mic o ed eléc ica basada en a ios con e ido es de
po encia
PLAN: 2024
DIRECTOR: MIGUEL CASTILLA FERNANDEZ
DEPARTAMENTO: INGENIERÍA ELECTRÓNICA
CALIFICACIÓN DEL TFM
DATA DE LECTURA: 27/05/2025
Es e p oyec o iene en cuen a los aspec os medioambien ales:

Sí

No
TRIBUNAL
PRESIDENTE
RAFAEL RAMON
RAMOS
LARA
SECRETARIO
JAUME
MIRET
TOMAS
VOCAL
RAFAEL
VIDAL
FERRÉ
RESUMEN
En las mic o edes eléc icas, la p esencia de ensiones de secuencia
nega i a debido a desbalances de ca ga o pe u baciones puede a ec a el
endimien o de los con e ido es, educi la ida ú il de los equipos y gene a
ine iciencias en la dis ibución de ene gía. Es e abajo p esen a un esquema
de con ol pa a la eliminación comple a de la ensión de secuencia nega i a
en una mic o ed con múl iples con e ido es de po encia.
La es a egia p opues a se basa en un con ol de gene ación de co ien e
de e e encia a pa i de la ensión de la ed i ual, c eada de mane a
independien e en cada con e ido . Pa a la de ección y sepa ación de las
componen es de secuencia nega i a, se emplea el Second-O de
Gene alized In eg a o (SOGI), un mé odo e icien e que pe mi e ex ae la
secuencia posi i a y nega i a sin necesidad de il os adicionales. Su
implemen ación mejo a la p ecisión del cálculo de e e encia y acili a la
compensación ac i a.
El con ol se implemen a en el ma co de e e encia es á ico al a-be a, lo
que pe mi e una mi igación e ec i a de los desequilib ios. A a és de
simulaciones en MATLAB/Simulink, se alida la e ec i idad del esquema,
demos ando la eliminación casi o al de la ensión de secuencia nega i a en
cada con e ido . Es o mejo a la es abilidad de la mic o ed y educe las
pé didas en los con e ido es.
Es e abajo o ece una solución e icien e pa a mejo a la calidad de la
ene gía en mic o edes con al a pene ación de ene gías eno ables y
ca gas desbalanceadas, asegu ando una ope ación más es able y con iable
del sis ema.
Palab as cla e (máximo 10):
Mic o ed
Desbalance
Secuencia posi i a
Secuencia nega i a
Con e ido
Red Vi ual
Ma lab
Simulink
SOGI
Calidad de ene gía
ABSTRACT
In elec ical mic og ids, he p esence o nega i e sequence ol ages due
o unbalanced loads o ne wo k dis u bances can a ec con e e
pe o mance, educe equipmen li espan, and c ea e ine iciencies in powe
dis ibu ion. This wo k p esen s a con ol scheme o he comple e elimina ion
o nega i e sequence ol age in a mic og id wi h mul iple powe con e e s.
The p oposed s a egy is based on a e e ence cu en gene a ion con ol
de i ed om he i ual g id ol age, c ea ed independen ly in each
con e e . To de ec and sepa a e he nega i e sequence componen s, he
Second O de Gene alized In eg a o (SOGI) is employed like an e icien
me hod ha ex ac s posi i e and nega i e sequences wi hou equi ing
addi ional il e s. I s implemen a ion enhances e e ence calcula ion
accu acy and acili a es ac i e compensa ion.
The con ol is implemen ed in he s a iona y e e ence ame alpha-be a,
allowing e ec i e mi iga ion o imbalances. Simula ions in MATLAB/Simulink
alida e he e ec i eness o he scheme, demons a ing an almos comple e
elimina ion o nega i e sequence ol age in each con e e . This imp o es
he mic og id s abili y and educes losses in he con e e s.
This wo k p o ides an e icien solu ion o enhance powe quali y in
mic og ids wi h high pene a ion o enewable ene gy and unbalanced loads,
ensu ing a mo e s able and eliable sys em ope a ion.
Keywo ds (10 maximum):
Mic og id
Unbalanced Load
Posi i e Sequence
Nega i e Sequence
Con e e
Vi ual G id
Ma lab
Simulink
SOGI
Powe Quali y
SUMARIO
1. INTRODUCCIÓN .......................................................................................................................... 8
1.1 CONTEXTO Y MOTIVACIÓN ....................................................................................................... 8
1.2 PLANTEAMIENTO DEL PROBLEMA Y JUSTIFICACIÓN ................................................................. 9
1.3 OBJETIVOS DE LA INVESTIGACIÓN ............................................................................................ 9
2. MARCO TEÓRICO ....................................................................................................................... 10
2.1 MICRORREDES ELÉCTRICAS .................................................................................................... 10
2.2 CONVERTIDORES DE POTENCIA EN MICRORREDES ................................................................ 11
2.3 COMPONENTES DE SECUENCIA POSITIVA Y NEGATIVA .......................................................... 12
2.4 CONTROL DE LA TENSIÓN DE SECUENCIA NEGATIVA EN MICRORREDES ................................ 14
3. DESARROLLO METODOLÓGICO .................................................................................................. 15
3.1 ARQUITECTURA DE LA MICRORRED ESTUDIADA .................................................................... 15
3.2 MODELADO DE MICRORRED SEPIC ......................................................................................... 16
3.3 ESTRATEGIA DE CONTROL PROPUESTA .................................................................................. 21
3.3.1 EXTRACCIÓN DE SECUENCIAS DE TENSIÓN POSITIVA Y NEGATIVA ................................... 21
3.3.1.1 TRANSFORMADA DE CLARKE ................................................................................... 21
3.3.1.2 EXTRACTOR DE SECUENCIAS DE DOBLE CANAL (SOGI + PSNC) ................................ 22
3.3.2 ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA A TRAVÉS DE UNA RED VIRTUAL .. 25
3.4 IMPLEMENTACIÓN DE SISTEMA EN ENTORNO MATLAB/SIMULINK ................................................ 29
3.4.1 MODELADO DE CADA CONVERTIDOR DE POTENCIA ......................................................... 29
3.4.2 MODELADO DE LÍNEA INTEGRADA A MICRORRED PARA CADA CONVERTIDOR DE POTENCIA
……………………………………………………………………………………………………………………………………….30
3.4.2.1 INVERSOR TRIFÁSICO ............................................................................................... 31
3.4.2.2 INDUCTANCIA DE FILTRO ......................................................................................... 31
3.4.2.3 CAPACITOR DE FILTRO .............................................................................................. 32
3.4.2.4 INDUCTANCIA LCL DE ACOPLAMIENTO .................................................................... 32
3.4.2.5 TRANSFORMADOR DE ACOPLAMIENTO ................................................................... 33
3.4.3 MODELADO DE FUENTE DE TENSIÓN PROGRAMABLE ...................................................... 33
3.4.4 MODELADO TOTAL DEL SISTEMA ...................................................................................... 35
4. SIMULACIÓN DEL SISTEMA ........................................................................................................ 36
4.1 CONDICIONES DE SIMULACIÓN .............................................................................................. 36
4.2 ESCENARIOS CONSIDERADOS ................................................................................................. 36
4.2.1 CONTROLADOR ACTIVADO ÚNICAMENTE CON INVERSOR G1 .......................................... 38
4.2.2 CONTROLADORES ACTIVADOS EN INVERSORES G1 Y G2 .................................................. 39
4.2.3 CONTROLADORES ACTIVADOS EN INVERSORES G1, G2 Y G3 ............................................ 40
4.2.4 CONTROLADORES ACTIVADOS EN INVERSORES G1 A G4 .................................................. 41
4.2.5 CONTROLADORES ACTIVADOS EN INVERSORES G1 A G5 .................................................. 43
4.2.6 CONTROLADORES ACTIVADOS EN LOS SEIS INVERSORES (G1 A G6) ................................. 44
5. ANÁLISIS DE RESULTADOS ......................................................................................................... 49
5.1. ESCENARIO 1: CONTROLADOR ACTIVADO ÚNICAMENTE CON INVERSOR G1 ......................... 49
5.2 ESCENARIO 2: CONTROLADORES ACTIVADOS EN INVERSORES G1 Y G2 .................................. 50
5.3 ESCENARIO 3: CONTROLADORES ACTIVADOS EN INVERSORES G1, G2 Y G3 ............................ 51
5.4 ESCENARIO 4: CONTROLADORES ACTIVADOS EN INVERSORES G1 A G4 .................................. 52
5.5 ESCENARIO 5: CONTROLADORES ACTIVADOS EN INVERSORES G1 A G5 .................................. 53
5.6 ESCENARIO 6: CONTROLADORES ACTIVADOS EN LOS SEIS INVERSORES (G1 A G6)...…….……… 54
6. ANÁLISIS DE IMPACTO AMBIENTAL, SOCIOECONÓMICO Y ÉTICO ............................................... 55
6.1 ANÁLISIS DE IMPACTO AMBIENTAL ........................................................................................ 55
6.2 ANÁLISIS DE IMPACTO SOCIOECONÓMICO ............................................................................ 55
6.3 ANÁLISIS ÉTICO ....................................................................................................................... 56

7. CONCLUSIONES Y TRABAJO FUTURO ......................................................................................... 57
7.1 CONCLUSIONES ...................................................................................................................... 57
7.2 TRABAJO FUTURO ................................................................................................................... 57
8. AGRADECIMIENTOS ................................................................................................................... 59
9. BIBLIOGRAFÍA ............................................................................................................................ 60
10. ANEXOS ................................................................................................................................ 62
10.1 ANEXO 1 (CODIGO GENERAL DE CONTROL) ............................................................................ 62
10.2 ANEXO 2 (CÓDIGO PARA OBTENCIÓN DE GRÁFICOS) ............................................................. 69
SUMARIO DE FIGURAS
FIGURA 1. EJEMPLO DE MICRORRED ELÉCTRICA ...................................................................................................... 8
FIGURA 2. COMPONENTES DE TENSIÓN DE SECUENCIA POSITIVA Y NEGATIVA EN UNA TENSIÓN DESEQUILIBRADA ................ 12
FIGURA 3. DIAGRAMAS DE BLOQUES DE CONTROL DE CONVERTIDORES GRID-FEEDING ................................................. 16
FIGURA 4. ESQUEMÁTICO DETALLADO DE MICRORRED SEPIC ................................................................................. 17
FIGURA 5. ESQUEMÁTICO ELÉCTRICO DE MICRORRED SEPIC ................................................................................... 17
FIGURA 6. ESQUEMÁTICO EN ENTORNO DE SIMULACIÓN MATLAB/SIMULINK DE MICRORRED SEPIC .............................. 18
FIGURA 7. EXTRACTOR DE SECUENCIAS DE TENSIÓN POSITIVA Y NEGATIVA .................................................................. 21
FIGURA 8. REPRESENTACIÓN DE TRANSFORMADA DE CLARKE EN BLOQUE DE CONTROL MATLAB/SIMULINK ..................... 22
FIGURA 9. REPRESENTACIÓN EN BLOQUES FUNCIONALES PARA EXTRACTOR DE SECUENCIA DOBLE CANAL (SOGI + PNSC) . 23
FIGURA 10. REPRESENTACIÓN DE ECUACIONES SOGI + PNSC EN BLOQUE DE CONTROL MATLAB/SIMULINK .................... 24
FIGURA 11. ESQUEMÁTICO ELÉCTRICO DE SECUENCIA POSITIVA POR CONVERTIDOR DE POTENCIA ................................... 25
FIGURA 12. ESQUEMÁTICO ELÉCTRICO DE SECUENCIA NEGATIVA POR CONVERTIDOR DE POTENCIA .................................. 25
FIGURA 13. REPRESENTACIÓN ELÉCTRICA DE RED VIRTUAL ..................................................................................... 26
FIGURA 14. REPRESENTACIÓN EN BLOQUES FUNCIONALES INTRODUCIENDO ESTRATEGIA DE CONTROL DE RED VIRTUAL ...... 26
FIGURA 15. CODIFICACIÓN DE VOLTAJES VIRTUALES EN MATLAB/SIMULINK ............................................................... 27
FIGURA 16. CODIFICACIÓN DE CORRIENTES ACTIVAS Y REACTIVAS DE SECUENCIAS POSITIVAS Y NEGATIVAS DE REFERENCIA
PROCEDENTES DE RED VIRTUAL EN MATLAB/SIMULINK ................................................................................ 28
FIGURA 17. CODIFICACIÓN DE CORRIENTES DE REFERENCIA EN MARCO ΑΒ PROCEDENTE DE RED VIRTUAL EN
MATLAB/SIMULINK ............................................................................................................................... 28
FIGURA 18. MODELO EN SIMULINK DE INVERSOR TRIFÁSICO TIPO GRID-FEEDING CON MEDICIÓN DE TENSIONES DE SALIDA .. 30
FIGURA 19. MODELO DE LÍNEA DE ACOPLAMIENTO ENTRE CONVERTIDOR DE POTENCIA Y MICRORRED EN SIMULINK ........... 30
FIGURA 20. PARAMETRIZACIÓN DE CADA IGBT CON DIODO ANTIRRETORNO EN CADA CONVERTIDOR DE POTENCIA ........... 31
FIGURA 21. PARAMETRIZACIÓN DE INDUCTANCIA DE FILTRO EN CADA CONVERTIDOR DE POTENCIA ................................ 31
FIGURA 22. PARAMETRIZACIÓN DE CAPACITANCIA DE FILTRO EN CADA CONVERTIDOR DE POTENCIA ............................... 32
FIGURA 23. PARAMETRIZACIÓN DE INDUCTANCIA LCL DE ACOPLAMIENTO EN CADA CONVERTIDOR DE POTENCIA ............. 32
FIGURA 24. PARAMETRIZACIÓN DE TRANSFORMADOR EN CADA CONVERTIDOR DE POTENCIA ........................................ 33
FIGURA 25. MODELADO DE FUENTE DE TENSIÓN PROGRAMABLE EN MATLAB/SIMULINK ............................................. 34
FIGURA 26. CODIFICACIÓN DE FUENTE DE TENSIÓN PROGRAMABLE EN MATLAB/SIMULINK .......................................... 34
FIGURA 27. MODELADO TOTAL DEL SISTEMA BASADO EN MICRORRED SEPIC EN ENTORNO MATLAB/SIMULINK ............... 35
FIGURA 28. RESPUESTA DE CONTROLADOR EN ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA EN EL INVERSOR G1 ...... 38
FIGURA 29. RESPUESTA DE CONTROLADOR EN ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA EN LOS INVERSORES G1 Y
G2 ..................................................................................................................................................... 39
FIGURA 30. RESPUESTA DE CONTROLADOR EN ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA EN LOS INVERSORES G1, G2
Y G3 ................................................................................................................................................... 40
FIGURA 31. RESPUESTA DE CONTROLADOR EN ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA EN LOS INVERSORES G1,
G2, G3 Y G4 ....................................................................................................................................... 42
FIGURA 32. RESPUESTA DE CONTROLADOR EN ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA EN LOS INVERSORES G1,
G2, G3, G4 Y G5 ................................................................................................................................. 43
FIGURA 33. RESPUESTA DE CONTROLADOR EN EVALUACIÓN INICIAL (ESTUDIO 1) DE ELIMINACIÓN DE TENSIÓN DE SECUENCIA
NEGATIVA EN LOS SEIS INVERSORES DE LA MICRORRED G1, G2, G3, G4, G5 Y G6 ............................................. 45
FIGURA 34. RESPUESTA DE CONTROLADOR SEGUNDA EVALUACIÓN (ESTUDIO 2) DE ELIMINACIÓN DE TENSIÓN DE SECUENCIA
NEGATIVA EN LOS SEIS INVERSORES DE LA MICRORRED G1, G2, G3, G4, G5 Y G6 ............................................. 46
FIGURA 35. RESPUESTA DE CONTROLADOR EN TERCERA EVALUACIÓN (ESTUDIO 3) DE ELIMINACIÓN DE TENSIÓN DE SECUENCIA
NEGATIVA EN LOS SEIS INVERSORES DE LA MICRORRED G1, G2, G3, G4, G5 Y G6 ............................................. 46
FIGURA 36. RESPUESTA DE CONTROLADOR EN CUARTA EVALUACIÓN (ESTUDIO 4) DE ELIMINACIÓN DE TENSIÓN DE SECUENCIA
NEGATIVA EN LOS SEIS INVERSORES DE LA MICRORRED G1, G2, G3, G4, G5 Y G6 ............................................. 47
FIGURA 37. RESPUESTA DE CONTROLADOR EN QUINTA EVALUACIÓN (ESTUDIO 5) DE ELIMINACIÓN DE TENSIÓN DE SECUENCIA
NEGATIVA EN LOS SEIS INVERSORES DE LA MICRORRED G1, G2, G3, G4, G5 Y G6 ............................................. 47
SUMARIO DE TABLAS
TABLA 1. DESCRIPCIÓN, REPRESENTACIÓN Y VALORES DE PARÁMETROS CORRESPONDIENTE A MICRORRED SEPIC ............. 37
TABLA 2. EVALUACIÓN DE LA ESTABILIDAD SEGÚN LA RELACIÓN ENTRE INDUCTANCIA VIRTUAL Y DE LÍNEA EN EL INVERSOR G1
......................................................................................................................................................... 38
TABLA 3. EVALUACIÓN DE LA ESTABILIDAD SEGÚN LA RELACIÓN ENTRE INDUCTANCIA VIRTUAL Y DE LÍNEA EN LOS INVERSORES G1
Y G2 ................................................................................................................................................... 39
TABLA 4. EVALUACIÓN DE LA ESTABILIDAD SEGÚN LA RELACIÓN ENTRE INDUCTANCIA VIRTUAL Y DE LÍNEA EN LOS INVERSORES
G1, G2 Y G3 ....................................................................................................................................... 41
TABLA 5. EVALUACIÓN DE LA ESTABILIDAD SEGÚN LA RELACIÓN ENTRE INDUCTANCIA VIRTUAL Y DE LÍNEA EN LOS INVERSORES
G1, G2, G3 Y G4 ................................................................................................................................. 42
TABLA 6. EVALUACIÓN DE LA ESTABILIDAD SEGÚN LA RELACIÓN ENTRE INDUCTANCIA VIRTUAL Y DE LÍNEA EN LOS INVERSORES
G1, G2, G3, G4 Y G5 ........................................................................................................................... 44
TABLA 7. ESTUDIOS DE ESTABILIDAD DE ELIMINACIÓN DE TENSIÓN DE SECUENCIA NEGATIVA SEGÚN LA RELACIÓN ENTRE
INDUCTANCIA VIRTUAL Y DE LÍNEA EN LOS SEIS INVERSORES DE LA MICRORRED G1, G2, G3, G4, G5 Y G6 ............. 48
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
8
1. INTRODUCCIÓN
Las mic o edes eléc icas, con o madas po múl iples con e ido es de po encia que
in eg an uen es de gene ación dis ibuida, equie en es a egias de con ol a anzadas
pa a ga an iza una ope ación es able y de al a calidad. Uno de los p oblemas c í icos
en es os sis emas es la p esencia de componen es de ensión de secuencia nega i a,
causados po desequilib ios en la ed o pe u baciones ex e nas, los cuales gene an
pé didas adicionales, deg adación de equipos y a ec ación en la calidad de la ene gía.
Es e abajo p opone un esquema de con ol o ien ado a la eliminación comple a de
la ensión de secuencia nega i a en una mic o ed basada en múl iples con e ido es.
Se desa olla un en oque de con ol dis ibuido y coo dinado que ac úa di ec amen e
sob e los con e ido es pa a compensa ac i amen e las componen es de secuencia
nega i a. La me odología incluye el modelado dinámico de la mic o ed, el diseño del
sis ema de con ol y su alidación median e simulaciones y análisis de desempeño bajo
di e en es condiciones ope a i as.
Figu a 1. Ejemplo de Mic o ed Eléc ica
1.1 CONTEXTO Y MOTIVACIÓN
La c ecien e pene ación de uen es de gene ación dis ibuida en los sis emas
eléc icos ha impulsado el desa ollo de mic o edes, las cuales pe mi en una ope ación
au ónoma y coo dinada de ecu sos ene gé icos locales. Es as mic o edes se
ca ac e izan po inclui con e ido es elec ónicos de po encia, que desempeñan un
papel undamen al en la conexión y con ol de uen es como paneles o o ol aicos,
u binas eólicas, bancos de ba e ías, en e o os.
Sin emba go, el en o no de ope ación de una mic o ed es dinámico y es á suje o a
condiciones como desequilib ios de ca ga, a iaciones en la gene ación, asime ías en
la ed de dis ibución y ansi o ios eléc icos. Es as condiciones pueden gene a
componen es de ensión de secuencia nega i a, las cuales de e io an la calidad del
suminis o eléc ico, causan ib aciones mecánicas en mo o es, inc emen an pé didas
po calen amien o y disminuyen la e iciencia de los con e ido es.
En es e con ex o, el con ol a anzado de los con e ido es es cla e pa a asegu a
una ope ación segu a, e icien e y es able de la mic o ed. La posibilidad de elimina
comple amen e la ensión de secuencia nega i a ep esen a un a ance signi ica i o en
el diseño de sis emas de con ol dis ibuidos, capaces de mejo a la calidad de ene gía
y ex ende la ida ú il de los equipos conec ados a la ed.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
15
3. DESARROLLO METODOLÓGICO
Es e capí ulo p esen a el desa ollo me odológico de la es a egia de con ol
p opues a pa a la eliminación ac i a de la ensión de secuencia nega i a en una
mic o ed i ásica compues a exclusi amen e po con e ido es g id- eeding. A
di e encia de con igu aciones adicionales que emplean con e ido es g id- o ming pa a
es abiliza la ensión, el sis ema conside ado en es a in es igación depende
comple amen e de la ed ex e na pa a su e e encia de ensión y ecuencia, lo que
impone e os adicionales en é minos de es abilidad y calidad de suminis o bajo
condiciones de desbalance.
El mé odo p opues o se undamen a en la implemen ación dis ibuida de un
Gene ado de Co ien e de Re e encia (GCR) basado en la ensión de una ed i ual
balanceada, es imada localmen e po cada con e ido . Es a ed i ual ac úa como
ma co de compa ación pa a de ec a y cuan i ica la p esencia de secuencia nega i a
en el sis ema. Median e écnicas de p ocesamien o en el dominio es á ico al a-be a y la
gene ación de co ien es de compensación especí icas, cada con e ido con ibuye
ac i amen e a mi iga la dis o sión causada po el desbalance, sin necesidad de una
a qui ec u a de con ol je á quico ni una uen e dominan e de ensión.
3.1 ARQUITECTURA DE LA MICRORRED ESTUDIADA
La mic o ed empleada como caso de es udio en es e abajo, denominada SEPIC
Mic oG id, es á con o mada po seis con e ido es elec ónicos de po encia del ipo g id-
eeding, conec ados en pa alelo a un pun o común de acoplamien o (PCC). Es a
con igu ación ep esen a un en o no con olado que pe mi e e alua con al a p ecisión
el compo amien o dinámico [10] del sis ema bajo dis in os escena ios ope a i os,
incluyendo desequilib ios de ensión, pe u baciones i ásicas y allos en la ed.
Cada con e ido de po encia es á alimen ado po una uen e de co ien e con inua
(DC) y conec ado a la ed a a és de una e apa de con e sión ull-b idge y un il o de
salida. Es os il os se diseñan con el obje i o de a enua las componen es a mónicas y
ga an iza una adecuada calidad de la co ien e inyec ada. Además, cada con e ido
incluye un lazo de con ol digi al, implemen ado median e un con olado independien e
en cada con e ido de po encia, que pe mi e ejecu a en iempo eal la es a egia de
compensación de la componen e de secuencia nega i a. La igu a 3 mues a un
diag ama de bloques de es e lazo de con ol.
La es a egia de con ol p opues a se basa en un en oque comple amen e dis ibuido,
donde cada con e ido ope a de o ma au ónoma e independien e, u ilizando una ed
i ual local pa a es ima la componen e de secuencia nega i a. A pa i de es a
es imación, se gene a una co ien e de e e encia cuya inyección iene como obje i o
cancela la ensión de secuencia nega i a en el PCC. Es a a qui ec u a e i a la
necesidad de comunicación en e con e ido es, lo que inc emen a la obus ez del
sis ema an e allos de enlace y acili a la escalabilidad de la mic o ed.
El sis ema cuen a además con elemen os de p ueba que pe mi en induci
condiciones ano males de ope ación. Es os incluyen impedancias desequilib adas,
disposi i os de desconexión de ase y mecanismos pa a simula allos con olados.
Tales elemen os son undamen ales pa a alida la e icacia de la es a egia de con ol
p opues a en e a pe u baciones eales, sin comp ome e la in eg idad del sis ema.

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
16
Finalmen e, la mic o ed inco po a una in aes uc u a de supe isión [11] que,
aunque no in e iene di ec amen e en el con ol p ima io, pe mi e egis a las a iables
eléc icas ele an es y ealiza ajus es expe imen ales. Es a pla a o ma esul a
especialmen e ú il pa a la e i icación de los esul ados median e simulaciones y
p uebas de labo a o io, apo ando e idencia cuan i a i a sob e la eliminación e ec i a de
la componen e de secuencia nega i a en condiciones dinámicas eales.
Figu a 3. Diag amas de bloques de con ol de con e ido es G id-Feeding
3.2 MODELADO DE MICRORRED SEPIC
El modelado de los con e ido es de po encia implemen ados en la mic o ed SEPIC,
mos ada en las igu as 4, 5 y 6, se ha desa ollado u ilizando el en o no de simulación
MATLAB/Simulink, con un en oque cen ado en la ep esen ación p ecisa del
compo amien o dinámico de cada uno de los seis con e ido es ipo g id- eeding. Cada
con e ido se conec a a la ed a a és de un il o LCL, lo que pe mi e mi iga el
con enido a mónico de las co ien es inyec adas y mejo a la calidad de ene gía en el
pun o de acoplamien o común.
En el modelo implemen ado, cada con e ido es á cons i uido po una e apa de po encia
basada en un in e so i ásico con olado median e modulación po ancho de pulso
(PWM), el cual es gobe nado po una es uc u a de con ol en lazo ce ado [12]. Es a
es uc u a incluye bloques especí icos pa a el seguimien o de e e encias de co ien e
gene adas a pa i de una ed i ual independien e pa a cada con e ido , lo que pe mi e
implemen a la es a egia de con ol p opues a pa a la eliminación de la ensión de
secuencia nega i a. Cada con olado gene a las señales de dispa o que de e minan el
compo amien o del in e so , basándose en el e o en e las co ien es medidas y las
e e encias deseadas.
Los con e ido es se conec an a la ed median e sus espec i os il os LCL,
compues os po dos induc ancias sepa adas po un capaci o de acoplamien o, jun o
con esis encias pa ási as que ep esen an las pé didas en los componen es. En el
modelo ambién se ha incluido la impedancia de acoplamien o en e con e ido es, la
cual es esencial pa a ep esen a la in e acción en e unidades en la mic o ed.
Adicionalmen e, en el modelo de Simulink se ep esen an los elemen os de medición
y los lazos de ealimen ación necesa ios pa a la implemen ación del con ol en iempo
eal [13], incluyendo las ans o maciones a ma cos de e e encia es á icos
(abc/al abe a) y las ans o maciones in e sas, así como los bloques de con ol
p opo cional- esonan e (PRES) que egulan las componen es de co ien e ac i a y
eac i a.
Es e en o no de modelado pe mi e ealiza p uebas p ecisas sob e la e ec i idad del
con ol dis ibuido diseñado, e aluando de mane a indi idual el compo amien o de cada
con e ido y su con ibución al echazo de la componen e de secuencia nega i a.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
17
Asimismo, acili a el análisis de la es abilidad del sis ema en e a pe u baciones,
condiciones de ope ación desequilib adas y a iaciones en los pa áme os de la ed.
Figu a 4. Esquemá ico de allado de Mic o ed SEPIC
Figu a 5. Esquemá ico eléc ico de Mic o ed SEPIC
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
18
Figu a 6. Esquemá ico en en o no de simulación Ma lab/Simulink de Mic o ed SEPIC
El modelamien o ma emá ico indi idual de cada uno de los seis con e ido es en
modo g id- eeding (G1 al G6) pe mi e desc ibi con p ecisión cómo in e ac úan con la
ed, cómo se p opagan las pe u baciones de secuencia nega i a y cómo las
ca ac e ís icas eléc icas locales in luyen en la es abilidad global del sis ema [14].
Dado que es os con e ido es ope an como uen es de co ien e con oladas, su
compo amien o depende di ec amen e de la e e encia de co ien e impues a po el
sis ema de con ol, y de las condiciones de ensión en el pun o de acoplamien o. A a és
del modelado, es posible ep esen a dinámicamen e la inyección de co ien e ac i a y
eac i a, así como las in e acciones en e nodos a a és de las impedancias de línea,
il os de salida y caídas de ensión asociadas.
Es e modelado no solo pe mi e e alua el e ec o de las asime ías y desbalances de
ensión sob e el desempeño del sis ema, sino que ambién posibili a el diseño e
implemen ación de es a egias a anzadas de con ol, o ien adas a la sup esión de
componen es de secuencia nega i a, la mejo a de la calidad de ene gía y la op imización
de la dis ibución de po encia en e los seis con e ido es de po encia.
Es e en oque es esencial pa a el diseño de es a egias de con ol que ga an icen una
ope ación equilib ada, obus a y e icien e de oda la mic o ed, especialmen e en
escena ios con múl iples gene ado es dis ibuidos, ca gas no lineales o allas
desequilib adas. La mic o ed mos ada en la igu a 5 se modela median e el conjun o
de ecuaciones disc e as que se esc iben seguidamen e:
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
19
𝑉0
,[𝑘]=#−𝑤1(𝐿20 +𝐿0&)𝐼30
,[𝑘]−#𝑤1𝐿&(C𝐼30
,[𝑘]+𝐼3&
,[𝑘]D−#𝑤1𝐿(4C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+
#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+𝐼35
,[𝑘]+𝐼36
,[𝑘]D#−𝑤1𝐿7C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+
𝐼35
,[𝑘]+𝐼36
,[𝑘]D+#𝑉7
,[𝑘] (3.1)
𝑉&
,[𝑘]=#−𝑤1𝐿2&C𝐼30
,[𝑘]+𝐼3&
,[𝑘]D−#𝑤1𝐿&(C𝐼30
,[𝑘]+𝐼3&
,[𝑘]D−#𝑤1𝐿(4C𝐼30
,[𝑘]+
𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+𝐼35
,[𝑘]+𝐼36
,[𝑘]D#−𝑤1𝐿7C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+
𝐼34
,[𝑘]+𝐼35
,[𝑘]+𝐼36
,[𝑘]D+#𝑉7
,[𝑘] (3.2)
𝑉(
,[𝑘]=#−𝑤1𝐿2(C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+#𝐼35
,[𝑘]+𝐼36
,[𝑘]D−#𝑤1𝐿(4C𝐼30
,+𝐼3&
,+
#𝐼3(
,+#𝐼35
,+𝐼36
,[𝑘]D[𝑘]−𝑤1𝐿7C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+𝐼35
,[𝑘]+𝐼36
,[𝑘]D+
#𝑉7
,[𝑘]# (3.3)
𝑉4
,[𝑘]=#−𝑤1𝐿24C𝐼30
,+𝐼3&
,+#𝐼3(
,+#𝐼35
,+𝐼36
,[𝑘]D[𝑘]−#𝑤1𝐿7C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+
#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+𝐼35
,[𝑘]+𝐼36
,[𝑘]D+#𝑉7
,[𝑘] (3.4)
𝑉5
,[𝑘]=#−𝑤1(𝐿25 +𝐿(5)(𝐼35
,[𝑘]+𝐼36
,[𝑘])−𝑤1𝐿(4C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+
#𝐼35
,[𝑘]+𝐼36
,[𝑘]D−#𝑤1𝐿7C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+𝐼35
,[𝑘]+𝐼36
,[𝑘]D+#𝑉7
,[𝑘]
(3.5)
𝑉6
,[𝑘]=#−𝑤1(𝐿26 +𝐿56)(𝐼36
,[𝑘])−𝑤1𝐿(5C#𝐼35
,[𝑘]+𝐼36
,[𝑘]D+𝑤1𝐿(4C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+
#𝐼3(
,[𝑘]+#𝐼35
,[𝑘]+𝐼36
,[𝑘]D−#𝑤1𝐿7C𝐼30
,[𝑘]+𝐼3&
,[𝑘]+#𝐼3(
,[𝑘]+𝐼34
,[𝑘]+𝐼35
,[𝑘]+
𝐼36
,[𝑘]D+#𝑉7
,[𝑘] (3.6)
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
20
donde:
• 𝑉8)0
,[𝑘]#: Secuencia de ensión nega i a pa a los seis con e ido es de
po encia (G1 al G6).
• 𝑤𝑜: F ecuencia angula nominal del sis ema.
• 𝐿29: Induc ancia conec ada a cada con e ido de po encia donde i = 1, 2, 3,
4, 5, 6.
• 𝐿&(,𝐿(4,𝐿(5,𝐿56,𝐿7: Induc ancias de líneas en e con e ido es de po encia.
• 𝐼39
,#: Co ien es de secuencia nega i a gene adas po cada in e so donde i =
1, 2, 3, 4, 5, 6.
• 𝑉7
,[𝑘]: Tensión de secuencia nega i a en el nodo común a ed.

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
21
3.3 ESTRATEGIA DE CONTROL PROPUESTA
La es a egia de con ol p opues a en es e abajo iene como obje i o p incipal la
eliminación comple a de la ensión de secuencia nega i a (𝑉,) en el pun o común de
acoplamien o de la mic o ed SEPIC. Es a es a egia se implemen a de o ma
descen alizada [15], median e un con ol local pa a cada uno de los seis con e ido es
de po encia g id- eeding [16]. Todos los con e ido es ienen un con ol basado en un
esquema de gene ado de e e encia de co ien e con ed i ual independien e eniendo
como obje i o p incipal la sup esión comple a de la componen e de ensión de secuencia
nega i a. Además, es e con ol asegu a el man enimien o del suminis o de po encia
ac i a desde cada con e ido hacia la ed, así como asegu a la compa ibilidad con
ope ación en desequilib io sin a ec a la calidad del suminis o i ásico de modo que en
el dominio abc las ampli udes de ensión sean iguales.
3.3.1 EXTRACCIÓN DE SECUENCIAS DE TENSIÓN POSITIVA Y
NEGATIVA
La sepa ación de las componen es de secuencia posi i a y nega i a de la ensión
i ásica es un paso undamen al en la es a egia de con ol p opues a, ya que pe mi e
iden i ica y compensa ac i amen e el desequilib io p esen e en la ed a a és de la
sepa ación de las componen es de secuencia posi i a y nega i a. La co ec a ex acción
de es as secuencias, al como se mues a en la Figu a 7, pe mi e no solo e alua la
calidad de la ene gía, sino ambién diseña la es a egia de con ol que mi igue los
e ec os ad e sos sob e los con e ido es de po encia y la es abilidad global de la
mic o ed.
Figu a 7. Ex ac o de secuencias de ensión posi i a y nega i a
3.3.1.1 TRANSFORMADA DE CLARKE
La T ans o mada de Cla ke es una he amien a ma emá ica undamen al u ilizada en
el análisis y con ol de sis emas eléc icos i ásicos, especialmen e cuando se equie e
sepa a o analiza enómenos elacionados con la calidad de ene gía [17], como
secuencias de ensión posi i a y nega i a pe mi iéndonos con e i señales i ásicas
desbalanceadas ( !, ",# %)#en un sis ema de coo denadas o ogonales de dos
componen es ( :,# ;).
La ep esen ación ma icial gene al pa a ensiones y co ien es de la T ans o mada
de Cla ke se mues a a con inuación:
/ :
;1=#1/3#H2 −1 −1
0√3 −√3KL !
"
%M (3.7)
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
22
Hi:
i;K=#1/3#H2 −1 −1
0√3 −√3KLi!
i"
i%M (3.8)
La igu a 8 mues a el código pa a es as dos ans o madas incluido en una Ma lab
unc ion.
Figu a 8. Rep esen ación de T ans o mada de Cla ke en bloque de con ol Ma lab/Simulink
3.3.1.2 EXTRACTOR DE SECUENCIAS DE DOBLE CANAL
(SOGI + PSNC)
El Ex ac o de Secuencia de Doble Canal nos pe mi e log a una ex acción p ecisa
de las componen es de secuencia posi i a y nega i a lo que es undamen al pa a el
diseño de la es a egia de con ol a pa i de la composición de dos bloques uncionales
p incipales: El In eg ado Gene alizado de Segundo O den (SOGI) y el Cálculo de
Secuencias Posi i a y Nega i a (PNSC).
El bloque SOGI se enca ga de p ocesa las ensiones ans o madas al ma co αβ, es
deci ( :#𝑦# ;), ob enidas median e la ans o mada de Cla ke. Cada componen e es
p ocesada de mane a independien e a a és de un in eg ado de segundo o den
sin onizado a la ecuencia undamen al de la ed ω, en donde 𝜔=#𝜔<=2𝜋60,
conside ando una ganancia, denominada k, con un alo √2 a a és de las siguien es
ecuaciones:
###########################𝑆𝑂𝐺𝐼=#=!
>?"(𝑠)=# @!A
A#)@!# (3.9)
########################𝐷(𝑠)=#=!
=(𝑠)=# >B<A
A#)>B!A)B# (3.10)
#############################################𝑄(𝑠)=#3=!
=(𝑠)=# >B!#
A#)>B!A)B# (3.11)
Es e p ocedimien o pe mi e ob ene señales de la componen e undamen al,
minimizando la in luencia de a mónicos y uidos de al a ecuencia. La es uc u a del
SOGI emplea una ealimen ación c uzada pa a mejo a la exac i ud en la gene ación de
la cuad a u a y aumen a la selec i idad de la banda de paso al como se mues a en la
Figu a 9.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
23
Figu a 9. Rep esen ación en bloques uncionales pa a Ex ac o de Secuencia Doble Canal (SOGI +
PNSC)
Pa e del p ocedimien o inal, las señales il adas y en cuad a u a son
pos e io men e u ilizadas en el bloque PNSC. Es e módulo ealiza combinaciones
especí icas en e las señales il adas y sus cuad a u as pa a calcula las componen es
de secuencia posi i a y nega i a:
• La secuencia posi i a (𝑣C
),#𝑣D
)) se ob iene combinando 𝑣C
< y 𝑣D
< de o ma cons uc i a
pa a cap u a la o ación de la ed en sen ido posi i o.
• La secuencia nega i a (𝑣C
,,#𝑣D
,) se calcula con combinaciones que e lejan la
in e sión de la di ección de o ación de las pe u baciones de secuencia nega i a.
Finalmen e p ocedemos a ep esen a ma emá icamen e la ampli ud y ase a pa i
de las secuencias posi i as y nega i as ob enidas:
a. Tensión de secuencia posi i a.
𝑉)=<(𝑣C
))&+C𝑣D
)D& (3.11)
b. Tensión de secuencia nega i a.
𝑉,=<(𝑣C
,)&+C𝑣D
,D& (3.12)
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
24
c. Fase a pa i de secuencias.
cos𝜑=#=$
%$=$
&,$='
%$='
&
E%E&# (3.13)
sin𝜑=#=$
%$='
&,$=$
&$='
&
E%E& (3.14)
𝜑=𝑎𝑡𝑎𝑛2(sin𝜑,cos𝜑) (3.15)
La igu a 10 mues a el código que pe mi e calcula las ensiones de secuencia, sus
ampli udes y ase en una Ma lab unc ion.
Figu a 10. Rep esen ación de ecuaciones SOGI + PNSC en bloque de con ol Ma lab/Simulink
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
31
3.4.2.1 INVERSOR TRIFÁSICO
El bloque del in e so gene a ensiones al e nas i ásicas median e una es uc u a
puen e de IGBTs con olados po señales PWM. El con ol de conmu ación es
p opo cionado po el bloque "In e e Con ol", que calcula las señales de e e encia
pa a log a la inyección adecuada de co ien e ac i a y/o eac i a, de acue do con el
con ol g id- eeding implemen ado. En la Figu a 20 se mues a la co ec a
pa ame ización de cada IGBT con diodo an i e o no.
Figu a 20. Pa ame ización de cada IGBT con diodo an i e o no en cada Con e ido de Po encia
3.4.2.2 INDUCTANCIA DE FILTRO
Después del in e so , las señales i ásicas pasan po una induc ancia de salida. Es a
se modela como una combinación se ie de una esis encia (R ) y una induc ancia (L 1)
los cuales son pa ame izados en la Figu a 21. Su p opósi o p incipal es a enua los
a mónicos gene ados po el p oceso de conmu ación, mejo a la o ma de onda de
co ien e y limi a los ansi o ios.
Figu a 21. Pa ame ización de Induc ancia de Fil o en cada Con e ido de Po encia

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
32
3.4.2.3 CAPACITOR DE FILTRO
Pos e io men e, se incluye una e apa capaci i a conec ada en es ella. Es a es á
compues a po una esis encia de pé didas Rc y una capaci ancia C 1, los cuales son
pa ame izados en la Figu a 22. Jun os, o man un il o LCL en combinación con la
siguien e induc ancia. Es a e apa mejo a el desempeño del sis ema en la mi igación de
dis o sión a mónica, especialmen e en ecuencias de conmu ación.
Figu a 22. Pa ame ización de Capaci ancia de Fil o en cada Con e ido de Po encia
3.4.2.4 INDUCTANCIA LCL DE ACOPLAMIENTO
La segunda induc ancia del il o LCL conec a la e apa capaci i a con el
ans o mado . Es e bloque ep esen a la induc ancia del lado de ed y es impo an e la
co ec a pa ame ización como se mues a en la igu a 23, que o ma pa e del diseño
del il o pasi o. Su alo es seleccionado pa a log a un comp omiso en e la a enuación
de a mónicos y la es abilidad del sis ema.
Figu a 23. Pa ame ización de Induc ancia LCL de Acoplamien o en cada Con e ido de Po encia
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
33
3.4.2.5 TRANSFORMADOR DE ACOPLAMIENTO
Finalmen e, se u iliza un ans o mado i ásico con impedancias se ie
ep esen adas po R 1 ( esis encia) y L 1 (induc ancia de uga) donde son
pa ame izados en la Figu a 24. Es e componen e ealiza el aislamien o gal ánico en e
el in e so y la mic o ed, y ambién adap a los ni eles de ensión eque idos pa a la
conexión al bus común.
Figu a 24. Pa ame ización de T ans o mado en cada Con e ido de Po encia
3.4.3 MODELADO DE FUENTE DE TENSIÓN PROGRAMABLE
En la p esen e sección se desc ibe el modelado ma emá ico y la implemen ación de
una uen e de ensión p og amable que pe mi e gene a componen es de secuencia
posi i a y nega i a, con oladas en magni ud y ase. Es e ipo de modelado nos pe mi e
emula una ed p incipal que puede p esen a desequilib ios con olados en su ensión,
pe mi iendo es udia el compo amien o de los con e ido es g id- eeding en e a
di e en es condiciones de secuencia nega i a.
La Figu a 25 mues a el bloque de simulación co espondien e. Den o del bloque,
se cons uyen las señales de ensión de ase !, "$y %#como suma de las componen es
de secuencia posi i a y nega i a. También se calcula la magni ud e ec i a de cada ase,
e ique ada como V!, V", V% que pe mi e analiza el impac o del desequilib io gene ado.
El código implemen ado en la Figu a 26 nos pe mi e ealiza el cálculo de las
magni udes de ase, componen es balanceadas de secuencia posi i a y nega i a,
además de la composición de las ensiones de ase log ando ga an iza la condición
i ásica en suma ( !+# "+## %=0), incluso cuando hay desequilib ios in encionados
en e las componen es.
Es e modelo pe mi e gene a un en o no de simulación con olado con desequilib ios
p og amables, necesa ios pa a alida la e icacia del con olado p opues o en la
eliminación de la secuencia nega i a. Al pode ajus a Vp, Vn y phi, se simulan dis in os
g ados de desequilib io y se obse a cómo el con e ido (g id- eeding) esponde bajo
el esquema de con ol desa ollado.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
34
Figu a 25. Modelado de Fuen e de Tensión P og amable en Ma lab/Simulink
Figu a 26. Codi icación de Fuen e de Tensión P og amable en Ma lab/Simulink
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
35
3.4.4 MODELADO TOTAL DEL SISTEMA
El modelo comple o desa ollado en Simulink ep esen a una mic o ed i ásica
compues a po seis con e ido es de po encia llamados In e so es en el en o no
Simulink (In e e 1 a In e e 6), una uen e de ensión p og amable, una ed de
in e conexión con impedancias dis ibuidas y con olado es independien es de
eliminación de secuencia nega i a pa a cada In e so . La Figu a 27 mues a la
es uc u a gene al del sis ema implemen ado.
Es e en o no de simulación en Simulink ha sido diseñado pa a analiza el
compo amien o colec i o de múl iples con e ido es conec ados a una ed
desequilib ada, con el in de alida el con ol p opues o pa a elimina la componen e de
secuencia nega i a en las co ien es inyec adas.
El modelo pe mi e ealiza dis in as p uebas en un en o no con olado, ales como:
• Ac i a o desac i a con e ido es indi idualmen e.
• Va ia el ni el de desequilib io gene ado po la uen e.
• E alua la calidad de la ensión y co ien e en el nodo común y en cada pun o de
conexión.
• Compa a el compo amien o del sis ema con y sin la es a egia de
compensación.
Además, se han inco po ado bloques de medición que pe mi en egis a las
a iables eléc icas ele an es ( ensiones y co ien es i ásicas), an o a la salida de
cada con e ido como en el nodo cen al.
Es e modelo p opo ciona una pla a o ma lexible y p ecisa pa a alida el con ol
p opues o y analiza su e ec i idad en la eliminación de la secuencia nega i a,
especialmen e en la mic o ed SEPIC con seis con e ido es de po encia del ipo g id-
eeding llamados In e so es.
Figu a 27. Modelado o al del sis ema basado en Mic o ed SEPIC en en o no Ma lab/Simulink
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
36
4. SIMULACIÓN DEL SISTEMA
En es e capí ulo se p esen a la simulación de la mic o ed implemen ada en Simulink.
El obje i o p incipal es alida el compo amien o del sis ema an e la p esencia de
ensión de secuencia nega i a y comp oba la e icacia del con olado p opues o pa a
elimina la componen e de ensión de secuencia nega i a en un en o no con múl iples
escena ios donde se hace la p ueba del con ol p opues o.
4.1 CONDICIONES DE SIMULACIÓN
Las simulaciones se ealiza on u ilizando un paso de mues eo de 𝑡A=1𝑥10,6#𝑠𝑒𝑔,
lo cual pe mi e cap u a con su icien e p ecisión las dinámicas ápidas del sis ema y una
ecuencia undamen al de la mic o ed de 60Hz.
La mic o ed incluye seis con e ido es de po encia llamados In e so es den o del
modelado en Simulink (In e e 1 a In e e 6), cada uno conec ado a a és de un
ans o mado , una línea con impedancia y un il o LCL, donde los pa áme os eléc icos
usados en el modelo se mues an en la Tabla 1.
4.2 ESCENARIOS CONSIDERADOS
Con el in de e i ica el co ec o uncionamien o del con olado diseñado pa a la
eliminación de la componen e de secuencia nega i a, se plan ea on di e en es
escena ios de p ueba en los que se a inc emen ando g adualmen e el núme o de
con e ido es (in e so es) que aplican el con ol.
El obje i o p incipal es e alua la espues a del sis ema cuando el con olado se
ac i a de o ma p og esi a: p ime o en un solo in e so , luego haciendo uso de dos
in e so es, luego es y así sucesi amen e has a llega a los seis in e so es uncionando
simul áneamen e con la es a egia de con ol p opues a. Es o no solo pe mi e obse a
cómo se compo a el sis ema en cada caso, sino ambién ajus a de mane a adecuada
los pa áme os de la ed i ual en cada con olado , es deci , la induc ancia y la
esis encia i uales.
Es os pa áme os son de e minan es pa a ga an iza la es abilidad del sis ema y una
dis ibución e ec i a de la co ien e de secuencia nega i a. Si no se eligen
co ec amen e, puede p oduci se una sob ecompensación o incluso ines abilidad en e
con e ido es.
Es e en oque p og esi o pe mi e iden i ica con p ecisión cómo se dis ibuye la
compensación en e los in e so es y qué alo es de pa áme os i uales son más
adecuados cuando más unidades es án in oluc adas en el p oceso de eliminación de
secuencia nega i a. De es e modo, se log a op imiza el diseño del con olado y
ga an iza un compo amien o obus o, incluso en con igu aciones con múl iples
con e ido es ope ando al mismo iempo.

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
37
ITEM
DESCRIPCIÓN DE
PARÁMETRO
REPRESENTACIÓN
VALOR
UNIDADES
1
Vol aje Co ien e Con inua
Vcc
400
V
2
Vol aje Nominal
Vo
110
V ms
3
Vol aje de Secuencia Posi i a
Vp_g id
155
V
4
Vol aje de Secuencia Nega i a
Vn_g id
35
V
5
F ecuencia de ed
𝜔
2
𝜋
60
ad/s
6
Induc ancia de Fil os Induc i os
L
3
mH
7
Resis encia de Fil os Induc i os
R
0.01
Ω
8
Capaci ancia de Fil os
Capaci i os
C
20
uF
9
Resis encia de Fil os
Capaci i os
Rc
24
Ω
10
Induc ancia de Fil os LCL
L_LCL
1.75
mH
11
Resis encias de Fil os LCL
R_LCL
0.01
Ω
12
Induc ancia de Red
Lg
2
mH
13
Resis encia de Red
Rg
0.01
Ω
14
Induc ancia de Línea Z12
L_z12
2
mH
15
Induc ancia de Línea Z23
L_z23
0.8
mH
16
Induc ancia de Línea Z34
L_z34
0.8
mH
17
Induc ancia de Línea Z35
L_z35
2
mH
18
Induc ancia de Línea Z56
L_z56
0.8
mH
19
Induc ancia de T ans o mado 1
L_T1
1
mH
20
Resis encia de T ans o mado 1
R_T1
0.5
Ω
21
Induc ancia de T ans o mado 2
L_T2
1
mH
22
Resis encia de T ans o mado 2
R_T2
0.5
Ω
23
Induc ancia de T ans o mado 3
L_T3
0.6
mH
24
Resis encia de T ans o mado 3
R_T3
1.13
Ω
25
Induc ancia de T ans o mado 4
L_T4
0.6
mH
26
Resis encia de T ans o mado 4
R_T4
1.13
Ω
27
Induc ancia de T ans o mado 5
L_T5
0.6
mH
28
Resis encia de T ans o mado 5
R_T5
1.13
Ω
29
Induc ancia de T ans o mado 6
L_T6
1.8
mH
30
Resis encia de T ans o mado 6
R_T6
1.95
Ω
Tabla 1. Desc ipción, ep esen ación y alo es de pa áme os co espondien e a Mic o ed
SEPIC
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
38
4.2.1 CONTROLADOR ACTIVADO ÚNICAMENTE CON
INVERSOR G1
Es e p ime escena io se ac i a el con olado de compensación de secuencia
nega i a únicamen e en el in e so G1. Es e con olado es á basado en un gene ado
de co ien e de e e encia que u iliza la ensión de una ed i ual, la cual inco po a una
induc ancia y una esis encia i ual ajus ables pa a modi ica el compo amien o
dinámico del sis ema de o ma local.
Los in e so es G2 a G6 pe manecen desac i ados, lo que signi ica que no gene an
co ien e du an e es a p ueba. De es a mane a, se e alúa exclusi amen e la espues a
del in e so G1 an e la p esencia de desequilib io en la ensión de ed y ob eniendo una
eliminación de secuencia nega i a a o able al como se mues a en la Figu a 28.
Además de e i ica la e ec i idad del con olado , en es a e apa se buscó iden i ica
cuál es la elación adecuada en e la induc ancia i ual L y la induc ancia ísica de
línea, pa a asegu a una espues a es able y e ec i a. En la Tabla 2 se mues an los
esul ados ob enidos al a ia es a elación.
Figu a 28. Respues a de con olado en eliminación de ensión de secuencia nega i a en el
in e so G1
In e so
Re e encia
L_ /L_i
L_
(mH)
Resul ado
G1
0.5
3.3
Ines able
1.0
6.6
Ok (Mejo esul ado)
2
13.2
Ok
Tabla 2. E aluación de la es abilidad según la elación en e induc ancia i ual y de línea en el
in e so G1
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
39
4.2.2 CONTROLADORES ACTIVADOS EN INVERSORES G1 Y G2
En es e segundo escena io se ac i an los con olado es de compensación de
secuencia nega i a en los in e so es G1 y G2, mien as que los in e so es es an es
(G3 a G6) pe manecen desac i ados. Al igual que en el caso an e io , cada con olado
gene a su co ien e de e e encia a pa i de la ensión de la ed i ual local, sin
necesidad de comunicación en e unidades. Es a con igu ación pe mi e e alua cómo
dos con e ido es con con ol ac i o compa en la a ea de elimina la componen e
nega i a y cómo a ec a la selección de pa áme os i uales en cada uno.
El compo amien o del sis ema mos ó mejo as signi ica i as en compa ación con el
caso an e io . Las co ien es inyec adas po G1 y G2 esul a on más simé icas y, en
conjun o, con ibuye on a una mayo educción de la componen e de secuencia nega i a
en el nodo común, al como se obse a en la igu a 29. Sin emba go, se obse ó que la
coo dinación implíci a en e los con olado es depende ue emen e de la co ec a
selección de los pa áme os i uales de cada in e so . La Tabla 3 esume los esul ados
ob enidos al a ia la elación en e la induc ancia i ual y la induc ancia ísica de línea
pa a cada in e so .
Figu a 29. Respues a de con olado en eliminación de ensión de secuencia nega i a en los
in e so es G1 y G2
In e so
Re e encia
L_ /L_i
L_
(mH)
Resul ado
G1
0,5
4,6
Ines able
1
9,2
Ok (Mejo esul ado)
2
18,4
Ok
G2
0,5
2,3
Ines able
1
4,6
Ok (Mejo esul ado)
2
9,2
Ok
Tabla 3. E aluación de la es abilidad según la elación en e induc ancia i ual y de línea en los
in e so es G1 y G2
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
40
4.2.3 CONTROLADORES ACTIVADOS EN INVERSORES G1, G2
Y G3
En es e e ce escena io se ac i a el con olado de compensación de secuencia
nega i a en es in e so es: G1, G2 y G3 donde cada con olado u iliza su p opio
gene ado de co ien e de e e encia basado en la ensión de una ed i ual local
mien as los in e so es G4 a G6 pe manecen desac i ados. El obje i o de es e caso es
obse a cómo se compo a el sis ema cuando la compensación se dis ibuye en e es
unidades ac i as, y e alua el impac o de los pa áme os i uales sob e la es abilidad y
e ec i idad de la es a egia.
Du an e las p uebas se iden i icó, una ez más, que la selección adecuada de la
induc ancia i ual L en cada in e so es undamen al pa a asegu a un epa o
adecuado del es ue zo de compensación sin gene a ines abilidades y log ando una
co ec a eliminación de la secuencia nega i a al como se mues a en la Figu a 30. La
Tabla 4 esume los esul ados ob enidos al a ia la elación en e la induc ancia i ual
y la induc ancia ísica de línea pa a cada in e so .
Figu a 30. Respues a de con olado en eliminación de ensión de secuencia nega i a en los
in e so es G1, G2 y G3
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
47
Figu a 36. Respues a de con olado en cua a e aluación (Es udio 4) de eliminación de ensión de
secuencia nega i a en los seis in e so es de la mic o ed G1, G2, G3, G4, G5 y G6
Figu a 37. Respues a de con olado en quin a e aluación (Es udio 5) de eliminación de ensión de
secuencia nega i a en los seis in e so es de la mic o ed G1, G2, G3, G4, G5 y G6

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
48
In e so
Re e encia
L_ /L_i
L_
(mH)
Resul ado
G1
0,5
26,1
Ines able
1
52,2
Es udio 1 (Acep able con oscilaciones)
1.05
54.6
Es udio 2 (Acep able con oscilaciones)
1.02
53
Es udio 3 (Acep able con oscilaciones)
0.96
50.2
Es udio 4 (Resul ado con oscilaciones
educidas)
1.12
58.2
Es udio 5 (Acep able con oscilaciones)
2
104,4
Ines able
G2
0,5
13,05
Ines able
1
26,1
Ok (Mejo esul ado)
2
52,2
Ines able
G3
0,5
4,35
Ines able
1
8,7
Ok (Mejo esul ado)
2
17,4
Ines able
G4
0,5
2,18
Ines able
1
4,35
Ok (Mejo esul ado)
2
8,7
Ines able
G5
0,5
13,05
Ines able
1
26,1
Ok (Mejo esul ado)
2
52,2
Ines able
G6
0,5
3,3
Ines able
1
52,2
Es udio 1 (Acep able con oscilaciones)
1.06
54.88
Es udio 2 (Acep able con oscilaciones)
1.02
53
Es udio 3 (Acep able con oscilaciones)
1.01
52.88
Es udio 4 (Resul ado con oscilaciones
educidas)
1.17
60.88
Es udio 5 (Acep able con oscilaciones)
2
13,2
Ines able
Tabla 7. Es udios de es abilidad de eliminación de ensión de secuencia nega i a según la elación
en e induc ancia i ual y de línea en los seis in e so es de la mic o ed G1, G2, G3, G4, G5 y G6
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
49
5. ANÁLISIS DE RESULTADOS
5.1 ESCENARIO 1: CONTROLADOR ACTIVADO ÚNICAMENTE
CON INVERSOR G1
Cuando solo el in e so G1 iene ac i ado el con olado de compensación de
secuencia nega i a, se obse a una educción pa cial de la componen e nega i a en la
co ien e que inyec a, en compa ación con el caso sin ningún con olado ac i o. Sin
emba go, como el es o de los in e so es (G2 a G6) pe manecen desac i ados, el e ec o
global sob e la ed es limi ado, ya que la compensación se ealiza solo desde un pun o
del sis ema.
Las o mas de onda de co ien e de G1 mues an una mejo a cla a en la sime ía
en e ases, y el análisis de secuencia indica una educción impo an e de la componen e
nega i a local. No obs an e, en el nodo común de la mic o ed oda ía se egis a una
co ien e de secuencia nega i a signi ica i a, lo que con i ma que la acción de un único
in e so no es su icien e pa a anula comple amen e el desequilib io cuando hay
múl iples uen es de dis o sión o con e ido es desa endidos.
Respec o a la es abilidad, los esul ados expe imen ales indican que el sis ema es
muy sensible al alo seleccionado pa a la induc ancia i ual L . Como se mues a en
la Tabla 4.1, cuando la elación L /Li es in e io a 1 (po ejemplo, 0.5), se p esen a
ines abilidad. Es o se mani ies a en oscilaciones en la co ien e de salida y en la ensión
i ual, indicando una espues a demasiado ápida o ag esi a del con olado .
Po el con a io, al usa una elación L /Li = 1.0, se log a un buen comp omiso en e
elocidad de espues a y es abilidad, siendo es e el mejo esul ado egis ado en la
abla. Un alo más ele ado ( elación 2.0) ambién o ece es abilidad, pe o con una
espues a algo más len a.
En esumen, los p incipales hallazgos del escena io 1 son:
• La es a egia de con ol basada en la ed i ual unciona co ec amen e cuando
se aplica de o ma indi idual.
• Un solo in e so no es su icien e pa a compensa odo el desequilib io del
sis ema.
• La selección de L es c í ica. Una elación L /Li = 1.0 se p esen a como una
opción óp ima.
• Es e escena io si e como base pa a escala la es a egia a más unidades,
dis ibuyendo la a ea de compensación de mane a coope a i a.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
50
5.2 ESCENARIO 2: CONTROLADORES ACTIVADOS EN
INVERSORES G1 Y G2
En es e escena io se ac i an los con olado es de compensación de secuencia
nega i a en los in e so es G1 y G2, mien as que los demás in e so es pe manecen
desac i ados. Es e caso pe mi e obse a cómo se compo a el sis ema cuando la a ea
de compensación comienza a dis ibui se en e dos unidades. También se e alúa si la
coope ación en e ambos con olado es mejo a la calidad de la co ien e en el nodo
común y cómo in luye la co ec a selección de los pa áme os i uales.
Los esul ados mues an una mejo a no able espec o al escena io an e io . La
componen e nega i a de co ien e en el nodo común disminuye de o ma más
signi ica i a y las o mas de onda de co ien e en G1 y G2 p esen an mayo sime ía
en e ases. Es o indica que la inclusión de un segundo con olado mejo a la e icacia
global de la compensación.
Sin emba go, es a mejo a depende en g an medida de la co ec a elección de la
induc ancia i ual en cada in e so . Al igual que en el escena io an e io , se obse a
que si la elación L /Li es demasiado baja, el sis ema se uel e ines able. Es o se
mani es ó median e oscilaciones no deseadas en la co ien e de salida, especialmen e
cuando ambos in e so es u ilizaban una elación de 0.5.
Po el con a io, cuando se selecciona una elación L /Li igual a 1.0 en ambos
con olado es, el sis ema ope a de o ma es able y e icien e. Los con olado es no
in e ie en en e sí, y la compensación se dis ibuye de o ma adecuada. Una elación
mayo , como 2.0, ambién man iene la es abilidad, aunque con una espues a más len a.
Es e escena io demues a que es posible log a una compensación más e ec i a
cuando más de un in e so pa icipa ac i amen e, siemp e que los pa áme os i uales
es én bien ajus ados. También se con i ma que no se equie e comunicación en e
con olado es, ya que cada uno esponde únicamen e a su medición local de ensión.
En esumen, los esul ados del escena io 2 indican que:
• La pa icipación de un segundo in e so mejo a signi ica i amen e la eliminación
de la componen e de secuencia nega i a.
• La es abilidad del sis ema sigue siendo sensible a los alo es seleccionados de
induc ancia i ual.
• La elación L /Li igual a 1.0 esul a nue amen e la más adecuada.
• La es a egia es escalable, siemp e que los con olado es se ajus en de o ma
cohe en e.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
51
5.3 ESCENARIO 3: CONTROLADORES ACTIVADOS EN
INVERSORES G1, G2 Y G3
En es e escena io se ac i an los con olado es de compensación de secuencia
nega i a en los in e so es G1, G2 y G3. El es o de los in e so es (G4 a G6)
pe manecen desac i ados. Es a con igu ación pe mi e e alua cómo se compo a el
sis ema cuando la compensación se dis ibuye en e es unidades, lo cual se ap oxima
más a un en o no ope a i o ípico en una mic o ed con múl iples uen es.
Los esul ados mues an una mejo a cla a en la sime ía de las co ien es de ase
inyec adas po los in e so es ac i os. Al mismo iempo, la co ien e o al en el nodo
común p esen a una educción signi ica i a en la componen e de secuencia nega i a, lo
que indica que la compensación conjun a es más e ec i a que en los escena ios
an e io es. Es o con i ma que dis ibui el es ue zo en e a ios in e so es mejo a el
desempeño global sin necesidad de comunicación en e ellos.
Al analiza la in luencia de los pa áme os i uales, se obse ó que una ez más la
elación en e la induc ancia i ual y la induc ancia de línea es de e minan e pa a
man ene la es abilidad del sis ema. En los es in e so es, el uso de una elación L /Li
= 0.5 esul ó en ines abilidad. En es os casos se de ec a on oscilaciones y
compo amien os no deseados en las a iables de co ien e y ensión i ual.
Po o o lado, al es ablece una elación L /Li = 1.0 en G1, G2 y G3, el sis ema se
es abiliza y se alcanza el mejo esul ado en é minos de desempeño dinámico y calidad
de co ien e. Es e alo pe mi e una espues a su icien emen e ápida sin comp ome e
la es abilidad. Cuando se u ilizan elaciones mayo es, como 2.0, el sis ema se man iene
es able, aunque con una espues a más len a, especialmen e du an e ansi o ios.
Es e caso ambién pe mi ió e i ica que los es in e so es pueden ope a con sus
con olado es ac i os sin in e e encias, siemp e que los pa áme os es én bien
ajus ados. Cada uno compensa localmen e la componen e nega i a de acue do con la
ensión que pe cibe, sin necesidad de sinc onización con el es o.
En esumen, los esul ados del escena io 3 mues an que:
• El aumen o del núme o de in e so es ac i os mejo a signi ica i amen e la
compensación de secuencia nega i a.
• La es abilidad y e iciencia del sis ema depende de que odos los con olado es
engan una con igu ación cohe en e de induc ancia i ual.
• La elación L /Li = 1.0 con inúa siendo el pun o de mejo endimien o.
• El sis ema es escalable y sigue uncionando co ec amen e con con olado es
dis ibuidos y au ónomos.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
52
5.4 ESCENARIO 4: CONTROLADORES ACTIVADOS EN
INVERSORES G1 A G4
En es e escena io se ac i an los con olado es de compensación de secuencia
nega i a en los in e so es G1, G2, G3 y G4, mien as que los in e so es G5 y G6
pe manecen desac i ados. El obje i o de es e caso ue e alua el compo amien o del
sis ema cuando la mayo ía de los con e ido es de po encia pa icipan en la
compensación, y obse a si el epa o de la a ea en e cua o unidades pe mi e
alcanza un buen desempeño sin comp ome e la es abilidad.
Los esul ados ob enidos mos a on una mejo a conside able en la calidad de la
co ien e inyec ada al nodo común. Las o mas de onda de co ien e p esen an sime ía
en e ases en los cua o in e so es ac i os, y la componen e de secuencia nega i a
disminuye no ablemen e en compa ación con los escena ios an e io es. Es o indica que
la acción conjun a de los con olado es es e ec i a pa a con a es a el desequilib io de
ensión impues o po la ed.
En é minos de es abilidad, el sis ema se mos ó muy sensible a la elección de la
induc ancia i ual. En los cua o in e so es, solo la elación L /Li igual a 1.0 pe mi ió un
compo amien o es able. Cuando se u iliza on alo es po debajo ( elación 0.5) o po
encima ( elación 2.0), el sis ema p esen ó ines abilidades e iden es, ales como
oscilaciones en la co ien e y des iaciones en las ensiones i uales.
Es os esul ados e ue zan la impo ancia de man ene una cohe encia en la
con igu ación de los pa áme os i uales cuando se ac i a el con ol en múl iples
con e ido es. Aunque los con olado es ope an de mane a au ónoma, sus e ec os se
combinan en el nodo común, po lo que un desequilib io en la con igu ación puede
a ec a el compo amien o global del sis ema.
En esumen, los p incipales hallazgos del escena io 4 son los siguien es:
• Cua o in e so es ac i os pe mi en una compensación e ec i a de la
componen e nega i a de la ed.
• La selección de la induc ancia i ual sigue siendo un ac o c í ico pa a asegu a
la es abilidad.
• Se con i ma que la elación L /Li igual a 1.0 es la más adecuada pa a odos los
in e so es.
• El sis ema puede man ene se es able y e icien e si se ga an iza una
con igu ación cohe en e en e odos los con olado es ac i os.
Es e análisis alida la escalabilidad del en oque p opues o y sien a las bases pa a
ac i a la compensación en los cinco con e ido es en el escena io siguien e.

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
53
5.5 ESCENARIO 5: CONTROLADORES ACTIVADOS EN
INVERSORES G1 A G5
En es e escena io se ac i an los con olado es de compensación de secuencia
nega i a en los in e so es G1, G2, G3, G4 y G5. Solo el in e so G6 pe manece
desac i ado. La inalidad de es e caso es e alua la e ec i idad del esquema de con ol
cuando p ác icamen e oda la mic o ed con ibuye ac i amen e a la eliminación de la
componen e de secuencia nega i a, y examina cómo se compo a el sis ema con cinco
unidades abajando de o ma au ónoma pe o simul ánea.
Los esul ados ob enidos mues an que el sis ema log a un al o g ado de sime ía en
las co ien es de ase inyec adas po los in e so es ac i os. Además, la componen e de
secuencia nega i a medida en el nodo común se educe de o ma signi ica i a,
alcanzando alo es mínimos en compa ación con los escena ios an e io es. Es o
con i ma que, al aumen a el núme o de con olado es ac i os, el es ue zo de
compensación se dis ibuye de o ma más e icien e, educiendo la ca ga indi idual de
cada in e so .
Sin emba go, al igual que en escena ios an e io es, el sis ema demos ó se muy
sensible a la selección de la induc ancia i ual. Se obse ó que únicamen e la elación
L /Li igual a 1.0 ga an izó la es abilidad del sis ema en odos los in e so es. Cuando se
u iliza on elaciones más bajas o al as, se p esen a on compo amien os ines ables
como oscilaciones o desajus es en e la co ien e gene ada y la e e encia calculada.
Es e compo amien o con i ma que, aunque los con olado es ope an de mane a
descen alizada y sin comunicación, es esencial que los pa áme os de la ed i ual
es én bien ajus ados y sean cohe en es en e odos los con e ido es. Una mala
con igu ación en uno solo de ellos puede al e a la es abilidad gene al del sis ema.
En esumen, los p incipales hallazgos del escena io 5 son los siguien es:
• Cinco in e so es ac i os pe mi en una compensación casi comple a de la
componen e de secuencia nega i a en el sis ema.
• El sis ema es al amen e e icien e, pe o exige p ecisión en la con igu ación de los
pa áme os i uales.
• La elación L /Li igual a 1.0 sigue siendo el único alo que ga an iza es abilidad
en odos los casos.
• La cohe encia en los pa áme os en e con olado es es cla e pa a man ene un
compo amien o equilib ado y es able.
Es e escena io demues a que la es a egia p opues a se adap a bien a
con igu aciones con un g an núme o de con e ido es y an icipa un uncionamien o
óp imo cuando se ac i e el con ol en los seis in e so es.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
54
5.6 ESCENARIO 6: CONTROLADORES ACTIVADOS EN LOS
SEIS INVERSORES (G1 A G6)
En es e escena io se ac i a la es a egia de con ol de eliminación de ensión de
secuencia nega i a en los seis in e so es que con o man la mic o ed. Es a
con igu ación ep esen a la máxima exp esión del esquema p opues o, ya que odos los
con e ido es de po encia con ibuyen de o ma simul ánea a la compensación del
desequilib io de ed median e el gene ado de co ien e de e e encia basado en ensión
i ual.
Los esul ados gene ales indican que el sis ema log a una compensación e ec i a,
con una no able educción de la componen e de secuencia nega i a en el pun o común
de acoplamien o. Sin emba go, ambién se obse a que la es abilidad global del sis ema
depende ue emen e de la co ec a elección de la induc ancia i ual en cada in e so .
Tal como se de alla en la Tabla 7, en los in e so es G2, G3, G4 y G5 la mejo
espues a se ob u o con una elación L /Li = 1.0, mos ando un compo amien o es able
y una sup esión e icien e de la secuencia nega i a. En con as e, los in e so es G1 y
G6 p esen a on mayo sensibilidad, po lo que se ealiza on es udios adicionales con
alo es lige amen e ajus ados de L , que pe mi ie on obse a una mejo a p og esi a
en la es abilidad.
En pa icula , el Es udio 4, co espondien e a una elación L /Li = 0.96 en G1 y L /Li
= 1.01 en G6, a ojó el esul ado más a o able, combinando es abilidad mejo ada con
una compensación e ec i a. Es e hallazgo indica que, en cie os nodos, un ajus e ino
de los pa áme os i uales puede op imiza el compo amien o dinámico del sis ema
comple o.
Po an o, el análisis de es e escena io con i ma que:
• La ac i ación simul ánea del con ol en odos los in e so es pe mi e una
sup esión p ác icamen e o al de la ensión de secuencia nega i a.
• El sis ema man iene su es uc u a dis ibuida sin necesidad de comunicación
en e nodos, pe o equie e una con igu ación cohe en e y, en algunos casos,
indi idualizada de los pa áme os i uales.
• Es posible mejo a la es abilidad del sis ema con lige os ajus es en in e so es
cla e, lo que sugie e que es a egias de sin onización ina pod ían complemen a
el diseño gene al.
Es e escena io demues a que la es a egia p opues a es o almen e escalable y
adap able a sis emas más complejos, siemp e que se ga an ice una adecuada
con igu ación de cada unidad.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
55
6. ANÁLISIS DE IMPACTO AMBIENTAL,
SOCIOECONÓMICO Y ÉTICO
6.1 ANÁLISIS DE IMPACTO AMBIENTAL
El sis ema de con ol desa ollado iene un e ec o di ec o y ele an e sob e el en o no,
ya que mejo a el uncionamien o de mic o edes eléc icas al educi el desequilib io en
la ensión de ed. Es a mejo a se aduce en una educción de pé didas eléc icas y en
una mayo e iciencia en la ope ación de los equipos conec ados, lo que con ibuye al
uso más acional de los ecu sos ene gé icos.
Además, al acili a una ope ación más es able en edes que in eg an gene ación
eno able, como la sola o la eólica, se a o ece la in eg ación de uen es limpias. Es o
apoya la ansición ene gé ica y la disminución del impac o ambien al asociado al uso
de combus ibles ósiles. Si bien el sis ema no in e iene di ec amen e sob e uen es de
gene ación, su capacidad pa a mejo a la es abilidad y calidad del suminis o iene un
e ec o posi i o sob e el uso e icien e de la ene gía.
6.2 ANÁLISIS DE IMPACTO SOCIOECONÓMICO
Desde el pun o de is a social y económico, una solución como la p esen ada puede
se de u ilidad en comunidades que ope an con edes aisladas o con in aes uc u a
limi ada. La posibilidad de mejo a la calidad de la ene gía con con olado es dis ibuidos
y sin necesidad de comunicación en e nodos ep esen a una al e na i a accesible,
especialmen e en zonas u ales o egiones donde la expansión de edes eléc icas
con encionales es cos osa o in iable.
La mejo a en la es abilidad del sis ema eléc ico iene un impac o di ec o en la calidad
de ida de los usua ios, ya que educe in e upciones, luc uaciones y allas en los
equipos eléc icos. Es o ambién puede aduci se en una disminución de cos os
ope a i os, an o pa a usua ios esidenciales como pa a pequeñas indus ias, lo que
puede ene un e ec o posi i o en la economía local.
Po o o lado, la implemen ación de es e ipo de soluciones pod ía omen a la
o mación de p o esionales y écnicos especializados en con ol, elec ónica de po encia
y ges ión de edes mode nas, impulsando el desa ollo de capacidades locales.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
56
6.3 ANÁLISIS ÉTICO
Desde el pun o de is a é ico, el p oyec o esponde a p incipios undamen ales del
eje cicio de la ingenie ía, como el comp omiso con el bienes a social, la sos enibilidad
y el uso esponsable del conocimien o. Desa olla ecnologías que hagan más
e icien es los sis emas eléc icos sin depende de g andes in e siones, y que puedan
se eplicadas en dis in as escalas, ep esen a una o ma de con ibui a una dis ibución
más jus a de los ecu sos ecnológicos.
El en oque descen alizado y au ónomo de los con olado es ambién pe mi e que las
comunidades o usua ios engan mayo con ol sob e sus sis emas eléc icos, lo que
p omue e la au onomía ene gé ica y el espe o po las necesidades locales. En es e
sen ido, la solución plan eada no busca impone un modelo cen alizado, sino adap a se
a di e en es con ex os con c i e ios écnicos y humanos.
En conjun o, el abajo ealizado no solo busca esol e un p oblema écnico, sino
ambién p opone una he amien a que pueda se pa e de sis emas eléc icos más
e icien es, accesibles y sos enibles.
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
63
i isemp y(q p im_be a_be o e)
q p im_be a_be o e = 0;
end
%Replacing alues
p im_be a= p im_be a_be o e;
p im_alpha= p im_alpha_be o e;
q p im_alpha=q p im_alpha_be o e;
q p im_be a=q p im_be a_be o e;
%% Designing SOGI
%SOGI o alpha
p im_alpha= p im_alpha_be o e+ s*wo*(k*( alpha- p im_alpha)-
q p im_alpha);
q p im_alpha=q p im_alpha_be o e+ s*(wo* p im_alpha);
%SOGI o be ha
p im_be a= p im_be a_be o e+ s*wo*(k*( be a- p im_be a)-q p im_be a);
q p im_be a=q p im_be a_be o e+ s*(wo* p im_be a);
%Re eshing memo ies comple ing he in eg a e equa ions
p im_be a_be o e= p im_be a;
p im_alpha_be o e= p im_alpha;
q p im_alpha_be o e=q p im_alpha;
q p im_be a_be o e=q p im_be a;
%% Designing he PNSC
_alpha_pos=0.5*( p im_alpha-q p im_be a);
_be a_pos=0.5*(q p im_alpha+ p im_be a);
_be a_neg=0.5*( p im_be a-q p im_alpha);
_alpha_neg=0.5*(q p im_be a+ p im_alpha);
% Ampli udes and phase
Vp = sq ( _alpha_pos^2 + _be a_pos^2); %posi i e sequence ol age
Vn = sq ( _alpha_neg^2 + _be a_neg^2); %nega i e sequence ol age
sin_phi=( _alpha_pos* _be a_neg + _alpha_neg* _be a_pos)/(Vp*Vn);
cos_phi=( _alpha_pos* _alpha_neg - _be a_pos* _be a_neg)/(Vp*Vn);
% Building he phase (phi)
phi=a an2(sin_phi,cos_phi);
%Designing algho i mic o adding 2*pi alue in phi alue
i phi<0
phi=phi+2*pi;
end
%% Designing condi ionals i he Vol age is nega i e mus be ze o
% Va Vb Vc measu emen s
Vaa = (Vp^2+Vn^2+2*Vp*Vn*cos(phi));
i (Vaa<0)
Vaa=0;
end
Va=sq (Vaa);

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
64
Vbb = (Vp^2+Vn^2+2*Vp*Vn*cos(phi+2*pi/3));
i (Vbb < 0)
Vbb = 0;
end
Vb=sq (Vbb);
Vcc = (Vp^2+Vn^2+2*Vp*Vn*cos(phi-2*pi/3));
i (Vcc<0)
Vcc = 0;
end
Vc=sq (Vcc);
%% VIRTUAL GRID VOLTAGE
%ialpha de i a e
pe sis en ialpha _be o e
i isemp y(ialpha _be o e)
ialpha _be o e = ialpha;
end
de _ialpha _d = (ialpha - ialpha _be o e )/ s;
ialpha _be o e = ialpha;
%ibe a de i a e
pe sis en ibe a _be o e
i isemp y(ibe a _be o e)
ibe a _be o e = ibe a;
end
de _ibe a _d = (ibe a - ibe a _be o e )/ s;
ibe a _be o e = ibe a;
%Ob aining he i ual ol ages
alpha_ = alpha - R * ialpha - L *(de _ialpha _d );
be a_ = be a - R * ibe a - L *(de _ibe a _d );
%% SOGI FOR VIRTUAL GRID VOLTAGE DEVELOPED
pe sis en p im_be a _be o e p im_alpha _be o e
i isemp y( p im_be a _be o e)
p im_be a _be o e = 0;
end
i isemp y( p im_alpha _be o e)
p im_alpha _be o e = 0;
end
pe sis en q p im_alpha _be o e q p im_be a _be o e
i isemp y(q p im_alpha _be o e)
q p im_alpha _be o e = 0;
end
i isemp y(q p im_be a _be o e)
q p im_be a _be o e = 0;
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
65
end
%Replacing alues
p im_be a_ = p im_be a _be o e;
p im_alpha_ = p im_alpha _be o e;
q p im_alpha_ =q p im_alpha _be o e;
q p im_be a_ =q p im_be a _be o e;
%SOGI o alpha i ual
p im_alpha_ = p im_alpha _be o e+ s*wo*(k*( alpha_ - p im_alpha_ )-
q p im_alpha_ );
q p im_alpha_ =q p im_alpha _be o e+ s*(wo* p im_alpha_ );
%SOGI o be a i ual
p im_be a_ = p im_be a _be o e+ s*wo*(k*( be a_ - p im_be a_ )-
q p im_be a_ );
q p im_be a_ =q p im_be a _be o e+ s*(wo* p im_be a_ );
%Re eshing memo ies comple ing he in eg a e equa ions
p im_be a _be o e= p im_be a_ ;
p im_alpha _be o e= p im_alpha_ ;
q p im_alpha _be o e=q p im_alpha_ ;
q p im_be a _be o e=q p im_be a_ ;
%% Designing he PNSC o VIRTUAL GRID
_alpha_pos_ =0.5*( p im_alpha_ -q p im_be a_ );
_be a_pos_ =0.5*(q p im_alpha_ + p im_be a_ );
_be a_neg_ =0.5*( p im_be a_ -q p im_alpha_ );
_alpha_neg_ =0.5*(q p im_be a_ + p im_alpha_ );
% Ampli udes and phase
Vp_ = sq ( _alpha_pos_ ^2 + _be a_pos_ ^2); %posi i e sequence
ol age
Vn_ = sq ( _alpha_neg_ ^2 + _be a_neg_ ^2); %nega i e sequence
ol age
sin_phi_ =( _alpha_pos_ * _be a_neg_ +
_alpha_neg_ * _be a_pos_ )/(Vp_ *Vn_ );
cos_phi_ =( _alpha_pos_ * _alpha_neg_ -
_be a_pos_ * _be a_neg_ )/(Vp_ *Vn_ );
% Building he phase (phi)
phi_ =a an2(sin_phi_ ,cos_phi_ );
%Designing algho i mic o adding 2*pi alue in phi alue
i phi_ <0
phi_ =phi_ +2*pi;
end
%% Designing i ual condi ions i he i ual ol age is nega i e
mus be ze o
% Va Vb Vc measu emen s
Vaa_ = (Vp_ ^2+Vn_ ^2+2*Vp_ *Vn_ *cos(phi_ ));
i (Vaa_ <0)
Vaa_ =0;
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
66
end
Va_ =sq (Vaa_ );
Vbb_ = (Vp_ ^2+Vn_ ^2+2*Vp_ *Vn_ *cos(phi_ +2*pi/3));
i (Vbb_ < 0)
Vbb_ = 0;
end
Vb_ =sq (Vbb_ );
Vcc_ = (Vp_ ^2+Vn_ ^2+2*Vp_ *Vn_ *cos(phi_ -2*pi/3));
i (Vcc_ <0)
Vcc_ = 0;
end
Vc_ =sq (Vcc_ );
%% Powe con ol loop (Re e ence Gene a o om he i ual g id)
% Injec ing ac i e and eac i e cu en h ough posi i e sequence
IP_pos_ e = 2/3*P_ e /(Vp_ +0.0001); % Ac i e posi i e cu en
IQ_pos_ e = 0; % Reac i e posi i e cu en
IP_neg_ e = 0; % Ac i e nega i e cu en
IQ_neg_ e = (Vn_ - Vn_ e ) / (wo*L ); % Reac i e nega i e cu en
% Ac i a I_Q_neg_ e solo si el iempo es mayo a 0.25 s
i ime < 0.2
IQ_neg_ e = 0;
end
i ime >= 0.2
IQ_neg_ e = (Vn_ - Vn_ e ) / (wo * L );
end
I_pos_ e = sq (IP_pos_ e ^2 + IQ_pos_ e ^2);
i I_pos_ e > I_ a ed
IP_pos_ e = (IP_pos_ e * I_ a ed)/I_pos_ e ;
IQ_neg_ e = (IQ_neg_ e * I_ a ed)/I_pos_ e ;
I_pos_ e = sq (IP_pos_ e ^2 + IQ_pos_ e ^2);
end
%% Using alues om Vi ual GRID
%Building Cu en Loop om ac i e and eac i e posi i e cu en s
injec ed
ialpha_ e = (( _alpha_pos_ *IP_pos_ e +
_be a_pos_ *IQ_pos_ e )/(Vp_ +0.0001))+(( _alpha_neg_ *IP_neg_ e +
_be a_neg_ *IQ_neg_ e )/(Vn_ +0.0001));
ibe a_ e = (( _be a_pos_ *IP_pos_ e -
_alpha_pos_ *IQ_pos_ e )/(Vp_ +0.0001))+(( _be a_neg_ *IP_neg_ e -
_alpha_neg_ *IQ_neg_ e )/(Vn_ +0.0001));
%% cu en con ol loop
pe sis en yialpha_1 zialpha_1 yibe a_1 zibe a_1
i isemp y(yialpha_1)
yialpha_1=0;
end
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
67
i isemp y(zialpha_1)
zialpha_1=0;
end
i isemp y(yibe a_1)
yibe a_1 = 0;
end
i isemp y(zibe a_1)
zibe a_1=0;
end
yialpha = yialpha_1 + s*wo*(2*damp_ ac o *(ki*(ialpha_ e -ialpha)-
yialpha_1)-(wo*zialpha_1));
zialpha = zialpha_1 + s*yialpha;
yialpha_1 = yialpha;
zialpha_1 = zialpha;
yibe a = yibe a_1 + s*wo*(2*damp_ ac o *(ki*(ibe a_ e -ibe a)-
yibe a_1)-(wo*zibe a_1));
zibe a = zibe a_1 + s*yibe a;
yibe a_1 = yibe a;
zibe a_1 = zibe a;
dalpha=(kp*(ialpha_ e -ialpha) + yialpha + alpha)*2/Vdc;
dbe a =(kp*(ibe a_ e -ibe a) + yibe a + ibe a)*2/Vdc;
%% Ac i e and eac i es powe s
p = (3/2) * ( alpha * ialpha + be a * ibe a);
q = (3/2) * ( be a * ialpha - alpha * ibe a);
%% Lowpass Fil e o P y Q
pe sis en P_1 Q_1
i isemp y(P_1)
P_1 = 0;
Q_1 = 0;
end
% LPF pa ame e s
alpha = s * (wo/2);
P = P_1 + alpha * (p - P_1);
Q = Q_1 + alpha * (q - Q_1);
P_1 = P;
Q_1 = Q;
%% Ob aining Ia Ib Ic
Ia = sq (IP_pos_ e ^2 + IQ_neg_ e ^2 +
2*IP_pos_ e *IQ_neg_ e *cos(phi-pi/2));
Ib = sq (IP_pos_ e ^2 + IQ_neg_ e ^2 +
2*IP_pos_ e *IQ_neg_ e *cos(phi+2*pi/3-pi/2));
Ic = sq (IP_pos_ e ^2 + IQ_neg_ e ^2 +
2*IP_pos_ e *IQ_neg_ e *cos(phi-2*pi/3-pi/2));
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
68
Imax = max([Ia,Ib,Ic]);
%% Ac i e and eac i es powe s
% il e o he cu en
pe sis en Imax_1
i isemp y(Imax_1)
Imax_1 = 0;
end
Imax = Imax_1 + alpha*(Imax-Imax_1);
Imax_1 = Imax;
i Imax > I_ a ed
IP_pos_ e = (IP_pos_ e * I_ a ed )/Imax;
Ia = Ia*I_ a ed/Imax;
Ib = Ib*I_ a ed/Imax;
Ic = Ic*I_ a ed/Imax;
Imax = max([Ia,Ib,Ic]);
end
end

Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
69
10.2 ANEXO 2 (CÓDIGO PARA OBTENCIÓN DE GRÁFICOS)
clc; clea ; close all;
%% In his sec ion we a e decla ing he
s = 1e-6; % Sample Time
wo = 2 * pi * 60;
R = 0.5; % T ans o me Resis ance
L = 1e-3; % T ans o me Induc ance
Ri = 0.01; % Line Resis ance
Li = 1.6e-3; % Line Induc ance
Rg = 0.01; % G id Resis ance
Lg = 2e-3; % G id Induc ance
RL = 24; % Load Resis ance
V_neg_ e =0; % Nega i e ol age e e ence
sim("Exe cise_TFM_LBD.slx");
%% Vol ages and Cu en s plo ed
igu e (1)
subplo (3,2,1)
plo ( , Va/155, 'b', 'DisplayName', 'Va ( )','LineWid h', 1.2);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Vb/155, ' ', 'DisplayName', 'Vb ( )','LineWid h', 1.2);
hold on;
plo ( , Vc/155, 'm', 'DisplayName', 'Vc ( )','LineWid h', 1.2);
ylabel('Va Vb Vc (V) (P.u)');
xlabel('Time (s)');
legend('Va','Vb','Vc')
g id on;
subplo (3,2,2)
plo ( , P, 'b','DisplayName', 'Pc (KW)','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Q, ' ','DisplayName', 'Qc (KW)','LineWid h', 1.5);
ylabel('Pc - Qc [Kw]')
xlabel('Time (s)');
legend('P','Q')
g id on;
subplo (3,2,3)
plo ( ,IP_pos_ e , 'b', 'DisplayName', 'Ipp (A)','LineWid h', 1.5);
hold on;
plo ( ,IQ_neg_ e , ' ', 'DisplayName', 'Iqn (A)','LineWid h', 1.5);
ylabel('Ipp - Iqn [A]')
xlabel('Time (s)');
legend('Ipp', 'Iqn');
g id on;
subplo (3,2,4)
plo ( , Vp_g id, 'b', 'DisplayName', 'Vps (V)', 'LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 10);
hold on;
plo ( ,Vp, ' ', 'DisplayName', 'Vpc (V)','LineWid h', 1.5);
plo ( ,Vp_ , ' ', 'DisplayName', 'Vpc (V)','LineWid h', 1.5);
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
70
ylabel('Vps - Vpc [A]')
xlabel('Time (s)');
legend('Vpg', 'Vp', 'Vp_ ')
g id on;
subplo (3,2,5)
plo ( , Ia_i, 'b', 'DisplayName', 'I_a','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Ib_i, ' ', 'DisplayName', 'I_b','LineWid h', 1.5);
plo ( , Ic_i,'m', 'DisplayName', 'I_c','LineWid h', 1.5);
hold on;
plo ( , Imax, 'k','LineWid h', 1.5)
ylabel('Iabc_i - Imax (A)');
xlabel('Time (s)');
legend('Ia', 'Ib', 'Ic','Imax')
g id on;
subplo (3,2,6)
plo ( , Vn_g id, 'b', 'DisplayName', 'Vn_s (V)','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( ,Vn, ' ', 'DisplayName', 'Vn (V)','LineWid h', 1.5);
plo ( ,Vn_ , 'm', 'DisplayName', 'Vn (V)','LineWid h', 1.5);
ylabel('Vns - Vnc [V]')
xlabel('Time (s)');
legend('Vng', 'Vn', 'Vn_ ')
g id on;
%
%% INVERTER 2
igu e (2)
subplo (3,2,1)
plo ( , Va_2/155, 'b', 'DisplayName', 'Va ( )','LineWid h', 1.2);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Vb_2/155, ' ', 'DisplayName', 'Vb ( )','LineWid h', 1.2);
hold on;
plo ( , Vc_2/155, 'm', 'DisplayName', 'Vc ( )','LineWid h', 1.2);
ylabel('Va Vb Vc (V) (P.u)');
xlabel('Time (s)');
legend('Va','Vb','Vc')
g id on;
subplo (3,2,2)
plo ( , P_2, 'b','DisplayName', 'Pc (KW)','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Q_2, ' ','DisplayName', 'Qc (KW)','LineWid h', 1.5);
ylabel('Pc - Qc [Kw]')
xlabel('Time (s)');
legend('P','Q')
g id on;
subplo (3,2,3)
plo ( ,IP_pos_ e _2, 'b', 'DisplayName', 'Ipp (A)','LineWid h', 1.5);
hold on;
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
71
plo ( ,IQ_neg_ e _2, ' ', 'DisplayName', 'Iqn (A)','LineWid h', 1.5);
ylabel('Ipp - Iqn [A]')
xlabel('Time (s)');
legend('Ipp', 'Iqn');
g id on;
subplo (3,2,4)
plo ( , Vp_g id, 'b', 'DisplayName', 'Vps (V)', 'LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 10);
hold on;
plo ( ,Vp_2, ' ', 'DisplayName', 'Vpc (V)','LineWid h', 1.5);
plo ( ,Vp_ 2, ' ', 'DisplayName', 'Vpc (V)','LineWid h', 1.5);
ylabel('Vps - Vpc [A]')
xlabel('Time (s)');
legend('Vpg', 'Vp', 'Vp_ ')
g id on;
subplo (3,2,5)
plo ( , Ia_i2, 'b', 'DisplayName', 'I_a','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Ib_i2, ' ', 'DisplayName', 'I_b','LineWid h', 1.5);
plo ( , Ic_i2,'m', 'DisplayName', 'I_c','LineWid h', 1.5);
hold on;
plo ( , Imax_2, 'k','LineWid h', 1.5)
ylabel('Iabc_i - Imax (A)');
xlabel('Time (s)');
legend('Ia', 'Ib', 'Ic','Imax')
g id on;
subplo (3,2,6)
plo ( , Vn_g id, 'b', 'DisplayName', 'Vn_s (V)','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( ,Vn_2, ' ', 'DisplayName', 'Vn (V)','LineWid h', 1.5);
plo ( ,Vn_ 2, 'm', 'DisplayName', 'Vn (V)','LineWid h', 1.5);
ylabel('Vns - Vnc [V]')
xlabel('Time (s)');
legend('Vng', 'Vn', 'Vn_ ')
g id on;
%% INVERTER 3
igu e (3)
subplo (3,2,1)
plo ( , Va_3/155, 'b', 'DisplayName', 'Va ( )','LineWid h', 1.2);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Vb_3/155, ' ', 'DisplayName', 'Vb ( )','LineWid h', 1.2);
hold on;
plo ( , Vc_3/155, 'm', 'DisplayName', 'Vc ( )','LineWid h', 1.2);
ylabel('Va Vb Vc (V) (P.u)');
xlabel('Time (s)');
legend('Va','Vb','Vc')
g id on;
subplo (3,2,2)
plo ( , P_3, 'b','DisplayName', 'Pc (KW)','LineWid h', 1.5);
Con ol pa a eliminación comple a de ensión de secuencia nega i a en una mic o ed
eléc ica basada en a ios con e ido es de po encia
LUIS BREÑA DIAZ
72
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Q_3, ' ','DisplayName', 'Qc (KW)','LineWid h', 1.5);
ylabel('Pc - Qc [Kw]')
xlabel('Time (s)');
legend('P','Q')
g id on;
subplo (3,2,3)
plo ( ,IP_pos_ e _3, 'b', 'DisplayName', 'Ipp (A)','LineWid h', 1.5);
hold on;
plo ( ,IQ_neg_ e _3, ' ', 'DisplayName', 'Iqn (A)','LineWid h', 1.5);
ylabel('Ipp - Iqn [A]')
xlabel('Time (s)');
legend('Ipp', 'Iqn');
g id on;
subplo (3,2,4)
plo ( , Vp_g id, 'b', 'DisplayName', 'Vps (V)', 'LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 10);
hold on;
plo ( ,Vp_3, ' ', 'DisplayName', 'Vpc (V)','LineWid h', 1.5);
plo ( ,Vp_ 3, ' ', 'DisplayName', 'Vpc (V)','LineWid h', 1.5);
ylabel('Vps - Vpc [A]')
xlabel('Time (s)');
legend('Vpg', 'Vp', 'Vp_ ')
g id on;
subplo (3,2,5)
plo ( , Ia_i3, 'b', 'DisplayName', 'I_a','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( , Ib_i3, ' ', 'DisplayName', 'I_b','LineWid h', 1.5);
plo ( , Ic_i3,'m', 'DisplayName', 'I_c','LineWid h', 1.5);
hold on;
plo ( , Imax_3, 'k','LineWid h', 1.5)
ylabel('Iabc_i - Imax (A)');
xlabel('Time (s)');
legend('Ia', 'Ib', 'Ic','Imax')
g id on;
subplo (3,2,6)
plo ( , Vn_g id, 'b', 'DisplayName', 'Vn_s (V)','LineWid h', 1.5);
se (gca, "linewid h", 2, " on size", 8);
hold on;
plo ( ,Vn_3, ' ', 'DisplayName', 'Vn (V)','LineWid h', 1.5);
plo ( ,Vn_ 3, 'm', 'DisplayName', 'Vn (V)','LineWid h', 1.5);
ylabel('Vns - Vnc [V]')
xlabel('Time (s)');
legend('Vng', 'Vn', 'Vn_ ')
g id on;
%% INVERTER 4
igu e (4)
subplo (3,2,1)
plo ( , Va_4/155, 'b', 'DisplayName', 'Va ( )','LineWid h', 1.2);
se (gca, "linewid h", 2, " on size", 8);