Programas Computacionais para Estudos de Harmônicos em Sistemas Elétricos
Programas Computacionais para Estudos de Harmônicos em Sistemas Elétricos
Neste artigo são abordados vários programas de tensão, mas na Fig. 2b é possível observar harmônicos de
computacionais para SEP que possuem ferramentas para corrente de 3a e 5a ordem.
análise da qualidade de energia, incluindo harmônicos. Neste
Espectro corrente da fase A
estudo foram consultados manuais dos programas e algumas Espectro tensão da fase A
4
versões de demonstração. Devido à constante atualização dos 200
I de pico (A)
3
V de pico (V)
150
programas, é possível que mudanças já tenham ocorrido em 100 2
0 0
manuais também não apresentam detalhes sobre modelos de 1 2 3 4 5 6 7 8 9 1 2 3 4 5 6 7 8 9
Ordem
componentes utilizados nas simulações recomendando-se Ordem
100
Tensão (V)
1,0 Tensão
0
-1,0 Corrente equivalentes harmônicas da rede para cada frequência de
-100 -3,0 interesse (varredura em frequência), possibilitando ainda o
-200 -5,0
cálculo das tensões e correntes para cada frequência nas
Tempo (ms)
diversas barras e linhas do sistema.
Fig. 1 – Tensão e corrente no domínio do tempo.
Para a análise harmônica, os métodos utilizados são
Neste domínio é possível observar que a forma de onda da varredura em frequência, injeção de corrente e fluxo de
corrente não é totalmente senoidal devido à existência de potência harmônico. A maioria dos programas possui
componentes harmônicas, porém não é possível identificar a ferramentas para dimensionamento de filtros para harmônicos.
frequência e a amplitude dos harmônicos presentes neste sinal.
(c) Cálculo da Distorção Harmônica Total
Esta informação somente é visualizada com a aplicação da
Finalmente, utilizando o princípio da superposição, é
Transformada Rápida de Fourier (FFT) para a obtenção dos
possível calcular a distorção harmônica total. Como exemplo
componentes harmônicos de tensão e corrente do sistema,
para um sistema elétrico, na Fig. 3 é apresentado o resultado
conforme mostrado na Fig. 2. Na Fig. 2a não há harmônicos
da varredura em frequência em uma barra de 4,16 kV, onde é
3
mostrada a impedância harmônica para cada freqüência, no realizado utilizando os modelos com correções típicas das
intervalo de 60 a 1200 Hz. Na fig. 4 são mostradas as tensões concessionárias NGC (Reino Unido) e EDF (França).
harmônicas resultantes no domínio da freqüência, quando um
conversor de 6 pulsos é ligado nesta barra como fonte de (3) Transformadores
harmônicos, calculados com o uso do programa PowerFactory. Os programas também permitem a representação de
transformadores de dois e três enrolamentos por suas
300
Amplitude (ohms)
20
15 (5) Bibliotecas de Modelos
10
5 Geralmente, os programas possuem bibliotecas de modelos
0
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 típicos de componentes do sistema elétrico. Entretanto, em
Ordem (x 60Hz)
muitos casos, os modelos disponíveis possuem valores típicos
Fig. 4 – Tensões harmônicas na presença de um conversor de 6 pulsos.
diferentes daqueles utilizados no Brasil, tais como tensões
B. Modelos de Componentes Elétricos para Análise de nominais diferentes, frequência e potências nominais dos
Harmônicos transformadores. Nestes casos, o usuário deve criar sua
A modelagem dos componentes dos sistemas elétricos para própria biblioteca.
estudos harmônicos é um tema bastante discutido na literatura C. Técnicas de Simulação do Sistema Elétrico
[6], [7]. Existem diversos modelos para a maioria dos
Outro aspecto importante é a escolha do tipo de sistema a
componentes. Os modelos dependem do tipo de rede em
ser simulado: monofásico ou trifásico.
estudo e dos níveis de precisão desejados. Alguns programas,
Os modelos monofásicos, representando somente a
com o objetivo de conseguir simplicidade e rapidez, utilizam
sequência positiva, geralmente são suficientes para sistemas
modelos simplificados [2]. Estes modelos podem levar a
de transmissão que geralmente são equilibrados.
resultados não satisfatórios em algumas situações.
Já nos demais sistemas, quando há a presença de cargas ou
(1) Modelos de Cargas bancos de capacitores monofásicos, ou harmônicos múltiplos
A maioria dos programas permite aos usuários escolherem de três, são necessários modelos trifásicos que trabalham com
os modelos de cargas a serem utilizados. Assim nos programas as sequências zero, positiva e negativa [4]. Nesta situação
PTW e no PSS/E, as cargas podem ser consideradas na também é possível realizar simulações de cargas harmônicas
simulação como modelo RL série ou paralelo. Os programas gerando harmônicos não característicos. O programa Power
Cymharmo, Superharm, HarmZS permitem também a Factory permite a simulações de sistemas desequilibrados,
modelagem de cargas utilizando modelos conforme o padrão calculando as distorções harmônicas de tensão e corrente por
CIGRÉ ou arranjos mistos. fase ou por sequência.
D. Facilidades de Uso e Aspectos Comerciais
(2) Modelos de Linhas de Transmissão e de
Distribuição (1) Interface Gráfica
As linhas podem ser representadas pelos modelos RL, Nem todos os programas existentes possuem interface
modelo π ou π com parâmetros distribuídos. gráfica. Em alguns, a interface somente está disponível para
saída de resultados. Isto dificulta, bem como aumentam as
2.a. Efeito Pelicular chances de erro na montagem e visualização do sistema e no
Alguns programas consideram o efeito pelicular sobre a preenchimento das características dos componentes.
resistência das linhas como no Cymharmo, Power Factory,
(2) Importação e Exportação de Resultados
PSS/E e Superharm, enquanto que em outros, como o PTW,
Outra característica importante é a função de exportação e
este efeito é desprezado.
importação de dados para outros aplicativos, o que pode
A modelagem do efeito pelicular é variada nos diversos
facilitar o tratamento posterior dos resultados e a elaboração
programas. Este efeito é modelado mantendo a relação X/R
de relatórios. Vários programas geram gráficos que não são de
constante como no programa PSS/E, como também pode ser
fácil compreensão ou que possuem tamanhos de letras e cores
informado pelo usuário, através de dados tabelados como no
que dificultam imprimi-los diretamente. Além disto, em
PSS/E, Superharm e Power Factory. No Power Factory, o
muitas situações, várias simulações podem ser necessárias e os
efeito pelicular também pode ser calculado utilizando uma
resultados devem ser comparados entre si para uma melhor
equação polinomial. Já no Cymharmo, este cálculo pode ser
avaliação.
4
(3) Versões de Programas para Avaliação o usuário necessita fazer um cadastro para aquisição posterior
A maioria dos programas analisados disponibiliza versões ou para solicitação aos representantes. Apesar de a maioria
para avaliação pelo usuário, incluindo versões com limitações dos fabricantes serem de origem estrangeira, vários possuem
do número de barras, versões com funções reduzidas ou representações no Brasil, fato que pode facilitar a troca de
versões que funcionam durante períodos determinados como experiências, obtenção de maiores informações ou de versões
15 ou 30 dias. Em alguns casos, estas versões podem ser de testes.
obtidas diretamente nos site dos fabricantes. Porém, em geral,
TABELA I
CARACTERÍSTICA DOS PROGRAMAS PARA ESTUDOS DE HARMÔNICOS
Módulo de
1 Interface Domínio p/ Versões Citações no
Programa País Ano análise de 2 3
gráfica harmônicos de teste IEEE
harmônicos
PSPICE USA 1984 Sim Principal Tempo Sim 202
ATP-EMTP 1984 Sim Principal Tempo (4) 75
PSCAD/EMTDC Canadá 1993 Sim Principal Tempo/ Frequência Sim 237
EMTP-RV Canadá 1987 Sim Principal Tempo Sim 52
MICROTRAN Canadá 1987 Parcial Principal Tempo/ Frequência Sim 2
PSIM USA 1997 Sim Principal Tempo/ Frequência Sim 61
MATLAB/SIMULINK USA 1984 Sim Principal Tempo Sim 909
PSAF / CYMFLOW Canadá 1986 Sim CYMharmo Frequência Sim 5
PTW / SKM USA 1972 Sim Hi-Wave Frequência Sim 3
ETAP USA 1986 Sim Principal Frequência Sim 6
POWER FACTORY Alemanha 1986 Sim Principal Frequência Sim 25
EDSA USA 1986 Sim Principal Frequência - 3
PSS/E SIEMENS USA/Alemanha 2004/1976 Sim PSS SINCAL’s Frequência Sim 5
SUPERHARM USA 1984 Parcial Principal Frequência Sim 3
HARM / CEPEL Brasil - Parcial Harmzs Frequência Sim 1
ASPEN USA 1986 Sim DistriView Frequência Sim -
SPARDO Colômbia 1992 Sim mp Power Frequência - -
1
Início de comercialização do programa quando disponível ou criação/aquisição da empresa de software
2
Incluem versões de teste, para universidades ou de estudante. Alguns demos necessitam de cadastro.
3
Busca no portal IEEE XPLORE = Nome do programa + harmonic em março 2011.
4
Este programa pode ser obtido gratuitamente por usuários que não participam da venda de outros programas de transitórios eletromagnéticos
semelhantes ao ATP.
indutâncias, resistências, capacitâncias e linhas sem perdas,
IV. PRINCIPAIS PROGRAMAS AVALIADOS incluindo uma chave e uma fonte de excitação.
Nesta seção são apresentados os programas avaliados por Divergências entre os desenvolvedores do EMTP levaram à
nome, endereço eletrônico e uma breve descrição das criação do ATP que atualmente possui um programa de
características. Na tabela II são mostrados a relação dos apoio chamado ATPDRAW, que proporcionou uma
programas, algumas de suas características e o número de interface gráfica para entrada de dados via Windows e
citações dos mesmos para análise harmônica através de buscas também um módulo de saída dos dados. O ATP pode
no portal IEEE explore. trabalhar no domínio do tempo utilizando valores discretos
ou no domínio da frequência. O software possui os
A. Programas no Domínio do Tempo principais elementos para modelagem de sistemas de
1. SPICE - www.orcad.com - O software SPICE foi transmissão, distribuição e industriais, além de possuir
desenvolvido para a simulação de circuitos eletrônicos em ferramentas necessárias para a modelagem de diversos
1972, na Universidade da Califórnia. Porém é possível elementos como dispositivos de eletrônica de potência [9].
simular e modelar circuitos e cargas harmônicas de
sistemas elétricos de potência. Elementos da rede elétrica 3. EMTP-RV - www.emtp.com - É também um programa
podem ser representados utilizando valores em pu e desenvolvido a partir do EMTP, em 1987 no Canadá. É
modelados utilizando componentes discretos, ativos ou utilizado para análise de transientes nos sistemas elétricos;
passivos. Desta forma podem-se criar modelos apropriados realiza estudos em regime permanente, no domínio do
de linha, cargas e outros elementos da rede e realizar a tempo, cálculo de fluxo de carga, análise de harmônicos e
simulação no domínio do tempo. Atualmente, versões varredura em frequência.
modificadas semelhantes ao SPICE, são comercializadas
por várias companhias usando vários nomes [12]. 4. PSCAD-EMTDC - https://pscad.com/ - foi desenvolvido
por Manitoba HVDC Research Center no Canadá. Possui
2. ATP-EMTP - www.emtp.org - O ATP é originado a partir interface gráfica e possibilita simular respostas transitórias
do EMTP, desenvolvido inicialmente na década de 1960 na em sistemas elétricos. Atualmente o desenvolvedor do
Alemanha. Quando foi criado, o programa só permitia PSCAD lançou o módulo independente FDHAP específico
modelagem de circuitos monofásicos através de para análise harmônica utilizando linhas de comando. Este
5
realiza o cálculo da amplitude e fase do DHT-V nos nós monitoramento em tempo real de sistemas elétricos de
para diversos tipos de fontes, varredura em frequência e transmissão, distribuição e industrial. Composto por
projeto de filtros [8]. módulo principal e outros quatro módulos adicionais:
módulo de simulação da monitoração, módulo de
5. MICRO TRAN - www.microtran.com - Também derivado gerenciamento de sistema de energia, módulo de controle
do EMTP, é desenvolvido pela empresa Microtran Power de demanda e módulo de automação de subestações. Possui
Systems Analysis Corporation, desde 1987 no Canadá. entrada e a saída de resultados por interfaces gráficas
Realiza a simulação de transitórios em sistemas elétricos de [18],[19].
potência causados por raios ou surtos de chaveamento.
Realiza também a análise de partida de motores, fenômenos 4. POWER FACTORY - www.digsilent.de – Programa da
de ressonância e análise da qualidade da energia, incluindo empresa Alemã Digsilent, comercializado desde 1986, que
harmônicos. Conta atualmente com interface gráfica [10], realiza a simulação e análise de sistemas elétricos de
[11]. potência de transmissão, distribuição e industriais. Executa
diversas funções como a análise de fluxo de carga, curto-
6. PSIM - www.powersimtech.com/ – Programa da empresa circuito, estabilidade do sistema, coordenação da proteção e
Powersim. Realiza a simulação de sistemas de eletrônica de análise harmônica. Realiza análise harmônica no módulo
potência, como conversores, acionamento de motores e principal com entrada de dados e saída de resultados
dinâmica de sistemas elétricos. Utiliza interface gráfica. A através de interface gráfica. Permite a simulação detalhada
partir da versão 9, foi incluída a simulação de fontes para redes equilibradas e desequilibradas exigindo desta
alternativas de energia com modelos de turbinas eólicas e forma um bom conhecimento do software para as
células solares. Permite a simulação de componentes do simulações [20], [21], [22], [23], [24], [25].
sistema elétrico e a criação de modelos para análise
específicas como a modelagem de cargas harmônicas [13]. 5. EDSA - www.poweranalytics.com – O software é
direcionado para projeto, simulação e análise de sistemas
7. MATLAB/SIMULINK - www.mathworks.com – O elétricos de potência. Comercializado desde 1986, possui
MATLAB é um software interativo de alto desempenho também módulos para análise em tempo real e para
voltado para a computação técnica, da empresa The controle de SEP. Possui módulo de qualidade de energia
MathWorks. Possui diversas ferramentas de cálculo que calcula a distorções harmônicas. Possui entrada de
numérico, simulação de sistemas e ferramentas para o dados e saída de resultados através de interface gráfica [26].
desenvolvimento de aplicações. Na plataforma do Matlab
tem-se o Simulink, que é a ferramenta para modelar, 6. PSS/E–Siemens - www.energy.siemens.com/hq/en/
simular e analisar sistemas dinâmicos. Dentre as diversas services/power-transmission-distribution/power-
bibliotecas do Simulink, tem-se o SimPowerSystem que technologies-international/ - Software desenvolvido pela
modela e simula sistemas elétricos de potência [14], com empresa Shaw Power Technologies, parte do grupo
ferramentas de cálculo das distorções harmônicas de Siemens desde 2004. Composto de vários módulos, é um
corrente, tensão e cálculo de valores eficazes. programa integrado de SEP que realiza simulação, análise e
otimização, principalmente na transmissão. O módulo PSS
B. Programas no Domínio da Frequência
SINCAL’s realiza análise de harmônicos [27], [28].
1. CYMFLOW - www.cyme.com - Módulo principal do 7. SUPERHARM - www.electrotek.com – É um dos módulos
software PSAF, atualmente parte do grupo Cooper autônomos do pacote de softwares PQSoft da empresa
Industries. Realiza a simulação de sistemas elétricos de Electrotek Concepts para análise de qualidade de energia e
transmissão, distribuição e industrial e possui sub-módulo eficiência energética. Este módulo é específico para
específico para harmônicos. Possui interface gráfica para modelagem e simulação de harmônicos em sistemas
entrada e saída de dados [15]. elétricos. A entrada de dados é através de linhas de código
e a saída de resultados por interface gráfica [29].
2. PTW DAPPER - www.skm.com – É um software da
empresa SKM Systems Analysis presente no mercado 8. HARM - www.cepel.br – Desenvolvido pelo CEPEL
desde 1972. Tem módulos integrados de simulação e especificamente para estudo de harmônicos nos sistemas de
análise de redes elétricas em sistemas de transmissão, potência. A versão mais atual do software é o HARMZS,
distribuição e sistemas industriais. Possui interface gráfica que substituiu o HARMZW. A entrada de dados é através
e permite diversas análises como correntes de falta, co- de linhas de código e a saída de resultados por interface
geração e análise de transitórios. O módulo HI-HAVE é o gráfica [30].
módulo responsável pela análise harmônica. Realizam 9. SPARDO - www.energyco.com/ – Pacote de programas da
varredura em frequência, cálculos da THD e projeto de Energy Computer Systems. O módulo SPARDO mp
filtros [16],[17]. Distribution realiza análise operacional, otimização e
planejamento de sistemas do distribuição como cálculo de
3. ETAP - www.etap.com – Programa para simulação e
6
fluxo de carga, análise de curto-circuito e gerenciamento de [14] Matlab. “User Guides.” Disponível on line <www.mathworks.com/help
/techdoc/index.html>. Acesso em 10/03/2011.
carregamento de transformadores. O módulo SPARDO mp [15] Cymharmo. “User’s Guide and Reference Manual.” March 2006, 104p.
Power realiza cálculos, modelagem e análise de sistemas de [16] Power Tools for Windows – “Tutorial.” 2008. 360p.
transmissão, distribuição e industriais, cálculos para ajuste [17] Power Tools for Windows – “Hi-Ware reference Manual.” 2006. 96p.
[18] ETAP – “User guide – Chapter 21- Harmonic analysis.” 48p.
da proteção e análise de harmônicos [31].
[19] ETAP 7.0.0 Demo – “Getting Started.” 2009. 144p.
[20] PowerFactory – “User’s Manual, Version 14.”Gomaringen, 2008.
10. ESA - www.easypower.com – Pacote de programas da 1100p.
empresa Easy Power. composto por vários módulos que [21] DIgSILENT “Technical Documentation - Induction Machine - TechRef
realizam modelagem, análise e monitoramento de redes. ElmAsm V2.1”, 2010, 21p.
[22] DIgSILENT “Technical Documentation - General Load Model -
Possui módulos para cálculo do fluxo de carga, análise de TechRef ElmLod V4”, 2008, 26p.
curto circuito, estabilidade do sistema, partida de motores e [23] DIgSILENT “Technical Documentation - Two-Winding Transformer (3-
o módulo Spectrum para análise de harmônicos [32]. Phase) - TechRef ElmTr2 V3”, 2007, 28p.
[24] DIgSILENT “Technical Documentation - Overhead Line Models - build
511”, 2009, 15p.
V. CONCLUSÕES [25] Wasilewski, J.; Wiechowski, W. and Bak, C. L. “Harmonic domain
modeling of a distribution system using the DIgSILENT PowerFactory
Os programas computacionais são uma ferramenta software.” In IEEE conference Future Power Systems. Amsterdam.
indispensável na análise de harmônicos nos sistemas elétricos. 2005.
Há diversos programas que fazem análise no domínio do [26] Paladin DesignBase Programs. Disponível on line
<www.poweranalytics.com> Acesso em 20/03/2011.
tempo e no domínio da frequência, sendo que os programas no [27] PSS/E 30 – Volume II - Program Application Guide, 2004. 598p.
domínio da frequência mostram-se mais adequados com maior [28] PSS/E Brochure. Disponível on line < www.energy.siemens.com>
rapidez e facilidades na modelagem para os grandes sistemas. Acesso em 20/03/2011.
Alguns programas não explicitam os modelos utilizados nas [29] SUPERHARM - User´s Guide – Version 4.3 – November 2002. 182p.
[30] CEPEL - HarmZs - Manual de Utilização do Programa. V. 1.7,2007,
simulações, não permitindo mesmo aos usuários mais 176p.
experientes avaliar o nível de precisão dos resultados [31] SPARD mp Power. Productos. Disponível on line
fornecidos. <http://www.energyco.com/SPP_e.html> Acesso em 10/08/2011.
[32] ESA-Spectrum Harmonic Analysis Software. Disponível on line
<www.easypower.com/products/EasyPower/asyPower_family.php>.
VI. AGRADECIMENTOS
Os autores agradecem à FAPEMIG e ao CNPq pelo suporte VIII. BIOGRAFIA
financeiro aos projetos e às atividades de pesquisa. Marcio Melquíades Silva é professor efetivo do Centro Federal de Educação
Tecnológica de Minas Gerais – CEFET-MG Possui graduação em engenharia
VII. REFERÊNCIAS elétrica, é Mestre pelo CEFET-MG e doutorando no programa de Pós-
Graduação em Engenharia Elétrica na Universidade Federal de Minas. Suas
[1] Dias G. A. “Harmônicas em Sistemas Industriais.” Porto Alegre. áreas de interesse são: Eficiência energética, qualidade da energia e
EDIPUCRS, 2º Ed., 2002. pp-159-168. instalações elétricas.
[2] Bam, L.; Jewell, W. “Review: Power System Analysis Software Tools.”
In IEEE Power Engineering Society General Meeting. 2005.
[3] Bengiamin, N.N; Holcomb, F.H. “PC-Basead Power systems software: Manuel Losada y Gonzalez é professor associado do Departamento de
comparing functions and features.” In IEEE Computer Applications in Engenharia da Universidade Federal de Minas Gerais (UFMG). Recebeu os
Power. Pp. 35-39; 1992. graus de barachel, mestre e doutor em Engenharia Elétrica pela UFMG nos
[4] Variz, Abílio Manuel. “Cálculo do Fluxo de Harmônicas em Sistemas anos de 1978, 1995 e 2000, respectivamente. Suas áreas de pesquisa incluem
Trifásicos Utilizando o Método de Injeção de Correntes.” Tese sistemas elétricos de potência, eficiência energética e qualidade de energia.
(Doutorado em engenharia elétrica) - COPPE/UFRJ, Rio de Janeiro,
2006. Eduardo Nohme é Professor associado do Departamento de Engenharia
[5] Chatterjea, A.; Mahmoud, A.; Harris, D.L. “The use of PSpice as a Elétrica da Universidade Federal de Minas Gerais (UFMG). Mestre pela PUC-
power system harmonics.” In IEEE Power Symposium. Pp.40–49. 1989 RJ em 1982 e Doutor pela UFSC em 1994 em Sistemas de Energia Elétrica,
[6] Task Force On Harmonic Modeling And Simulation, IEEE Power Eng. área de concentração estabilidade e segurança, incorporando mais
Soc. T&D Committee. “Impact of Aggregate Linear Load Modeling on recentemente em sua atuação na análise e operação de sistemas, os estudos
Harmonic Analysis: A Comparison of Common Practice and Analytical relativos à integração à rede de fontes renováveis/alternativas de energia.
Models.” In IEEE Transactions on Power Delivery, Vol. 18, No. 2, April
2003, pp. 625-630. Selênio Rocha Silva é professor titular do Departamento de Engenharia
[7] Xu, W. “Status and future directions of power system harmonic
Elétrica da Universidade Federal de Minas Gerais (UFMG). Recebeu os graus
analysis.” In Proceedings of the IEEE PES General Meeting, Colorado,
de bacharel e mestre em Engenharia Elétrica pela UFMG respectivamente nos
Vol. 1, pp.756-761, 2004.
anos de 1980 e 1984, e o título de doutor em Engenharia Elétrica pela
[8] Manitoba Hvdc Research Centre Inc. “Applications of PSCAD /
EMTDC.” 2008. 114p. Universidade Federal da Paraíba em Campina Grande (atual UFCG) em 1988.
[9] Velasco, Juan A. Martinez. “El paquete de simulación ATP. Nuevas Seus interesses incluem máquinas elétricas, acionamentos elétricos, qualidade
prestaciones.” Disponível em <http://www.iitree-unlp.org.ar/ da energia e sistemas de geração de fontes alternativas de energia.
caue/PapersRevista/2000-N3-El%20paquete%20ATP.pdf> Acesso em
03/03/2011.
[10] Micro Tran. “Reference Manual.” 2002. 21p.
[11] Mtdiagram.”Reference Manual - Graphical user interface and network
drafting tool for Microtran”. 190p.
[12] Veronese, P. R. “Simuladores de circuitos eletrônicos – mini-curso.”
Disponível em <http://www.sel.eesc.usp.br/sasel>. Acesso em
03/03/2011.
[13] PSIM. “User’s Guide. Version 9.0.” January 2010. 216P.