12.07.2015 Views

análise de observabilidade em sistemas elétricos ... - NUPET - UTFPR

análise de observabilidade em sistemas elétricos ... - NUPET - UTFPR

análise de observabilidade em sistemas elétricos ... - NUPET - UTFPR

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

UNIVERSIDADE TECNOLÓGICA FEDERAL DO PARANÁDEPARTAMENTO ACADÊMICO DE ELETROTÉCNICACURSO DE ENGENHARIA INDUSTRIAL ELÉTRICA / AUTOMAÇÃOALLISON DINIZ NOCERAGIANCARLO GOMESANÁLISE DE OBSERVABILIDADE EM SISTEMAS ELÉTRICOS DEPOTÊNCIATRABALHO DE CONCLUSÃO DE CURSOCURITIBA2013


ALLISON DINIZ NOCERAGIANCARLO GOMESANÁLISE DE OBSERVABILIDADE EM SISTEMAS ELÉTRICOS DEPOTÊNCIAProposta <strong>de</strong> Trabalho <strong>de</strong> Conclusão <strong>de</strong> Curso <strong>de</strong>Graduação, apresentado à disciplina <strong>de</strong>Metodologia Aplicada ao TCC, do curso <strong>de</strong>Engenharia Industrial Elétrica com ênfase <strong>em</strong>Automação, do Departamento Acadêmico <strong>de</strong>Eletrotécnica (DAELT) da Universida<strong>de</strong>Tecnológica Fe<strong>de</strong>ral do Paraná (<strong>UTFPR</strong>), comorequisito parcial para a obtenção do título <strong>de</strong>Engenheiro Eletricista com ênfase <strong>em</strong>Automação.Orientador: Prof. Dr. Raphael Augusto <strong>de</strong> SouzaBeneditoCURITIBA2013


LISTA DE ILUSTRAÇÕESFigura 1 – Fluxograma do Processo ......................................................................... 12


LISTA DE ABREVIATURAS E SIGLASEE – Estimação <strong>de</strong> EstadoIEEE – Institute of Electrical and Electronics EngineersSEP – Sist<strong>em</strong>a Elétrico <strong>de</strong> PotênciaUTR – Unida<strong>de</strong> Terminal R<strong>em</strong>ota


SUMÁRIO1 INTRODUÇÃO ......................................................................................................... 61.1 TEMA ................................................................................................................. 71.1.1 Delimitação do T<strong>em</strong>a .................................................................................. 71.2 PROBLEMA E PREMISSAS ................................................................................. 91.3 OBJETIVOS ........................................................................................................ 101.3.1 Objetivo Geral ............................................................................................... 101.3.2 Objetivos Específicos.................................................................................... 101.4 JUSTIFICATIVA .................................................................................................. 101.5 PROCEDIMENTOS METODOLÓGICOS ............................................................ 111.6 ESTRUTURA DO TRABALHO ............................................................................ 121.7 CRONOGRAMA .................................................................................................. 131.8 REFERÊNCIAS ................................................................................................... 13


7estimar as variáveis <strong>de</strong> estado. Essa dinâmica é o foco principal <strong>de</strong>ssa proposta <strong>de</strong>estudo, i<strong>de</strong>ntificar as ilhas e procurar a restauração da observabilida<strong>de</strong>, além <strong>de</strong>analisar métodos entre os muitos existentes e conhecer as suas diferenças evantagens particulares.Os sist<strong>em</strong>as <strong>de</strong> medição instalados são constituídos usualmente <strong>de</strong>medidores do tipo fluxo <strong>de</strong> potência ativa e reativa nos ramos ou linhas da re<strong>de</strong>,relativo a duas barras; injeção <strong>de</strong> potência ativa e reativa e magnitu<strong>de</strong>s das tensõesnos barramentos, relativos a uma barra da re<strong>de</strong> [1] .Essas medidas comumente são processadas <strong>de</strong> forma a haver<strong>em</strong>redundâncias, o que aumenta a confiabilida<strong>de</strong> da estimação <strong>de</strong> estado, além <strong>de</strong>suprir falhas <strong>em</strong> medidores e possibilitar a <strong>de</strong>tecção <strong>de</strong> erros grosseiros. Essasmedições são aferidas <strong>em</strong> sub-estações que po<strong>de</strong>m ser localizadas a gran<strong>de</strong>sdistâncias umas das outras, consi<strong>de</strong>rando as dimensões <strong>de</strong> alguns SEP, causando,<strong>de</strong>sta maneira possíveis falhas que po<strong>de</strong>m ser t<strong>em</strong>porárias ou não.1.1 TEMAAnálise <strong>de</strong> Observabilida<strong>de</strong> <strong>em</strong> Sist<strong>em</strong>as Elétricos <strong>de</strong> Potência.1.1.1 Delimitação do T<strong>em</strong>aA partir <strong>de</strong>ssa introdução po<strong>de</strong>-se inferir a importância do t<strong>em</strong>a escolhido, queespecificamente trata da análise da observabilida<strong>de</strong> <strong>em</strong> sist<strong>em</strong>as elétricos <strong>de</strong>potência; a avaliação da possibilida<strong>de</strong> <strong>de</strong> restaurar a observabilida<strong>de</strong>, seja porpseudo-medidas ou inserção <strong>de</strong> novas medidas; a i<strong>de</strong>ntificação <strong>de</strong> ilhas observáveise <strong>de</strong> medidas críticas, que são as medidas essenciais à observabilida<strong>de</strong>.Estudos com esse enfoque já são realizados <strong>de</strong>s<strong>de</strong> a década <strong>de</strong> 60, com oex<strong>em</strong>plo <strong>de</strong> Schweppe (SCHWEPPE; WILDES; RON, 1968; SCHWEPPE; WILDES,1970; SCHWEPPE; RON, 1970; e SCHWEPPE, 1970) que traçou tópicos qu<strong>em</strong>ostravam a complexida<strong>de</strong> e necessida<strong>de</strong> <strong>de</strong> se conhecer mais sobre o assunto [2] .Depois, a partir da década <strong>de</strong> 70, vieram estudos que já buscavam <strong>de</strong> formamais específica o foco no t<strong>em</strong>a, porém ainda s<strong>em</strong> apresentar efetiva aplicaçãoprática. Ex<strong>em</strong>plificam esse caso os métodos <strong>de</strong> Handschin e Bongers (1972) eFetzer e An<strong>de</strong>rson (1975), entre outros. Ainda nessa década, outro método que


8merece <strong>de</strong>staque é o <strong>de</strong> Cl<strong>em</strong>ents e Wollenberg (1975), que relataram umprocedimento on<strong>de</strong> a análise das medidas <strong>de</strong> fluxo é feita antes e então a análisedas medidas <strong>de</strong> injeção. Aprimorando esse algoritmo <strong>de</strong> Cl<strong>em</strong>ets e Wollenberg(1975), All<strong>em</strong>ong et al (1980) <strong>de</strong>u o primeiro passo da nova década nesses estudos[2] .Nessa década <strong>de</strong> 80, foram introduzidos os conceitos <strong>de</strong> observabilida<strong>de</strong>algébrica, numérica e topológica e com o conceito <strong>de</strong> teoria <strong>de</strong> grafos <strong>de</strong>senvolveuseum algoritmo para essa análise. Pouco mais tar<strong>de</strong> <strong>em</strong> 1985, Monticelli e Wu<strong>de</strong>senvolveram um método numérico a partir da fatoração triangular da matriz ganhoque faz a verificação da observabilida<strong>de</strong> da re<strong>de</strong> e, caso não seja observável comoum todo, permite i<strong>de</strong>ntificar as ilhas [2] .Já durante a década <strong>de</strong> 90, Bretas (1996) apresentou o método <strong>de</strong>observabilida<strong>de</strong> que utiliza o conceito <strong>de</strong> caminhos <strong>de</strong> fatoração, associado àfatoração triangular da matriz Ganho, para <strong>de</strong>terminação das ilhas observáveis erestauração da observabilida<strong>de</strong> [1] .Mais recente, o trabalho <strong>de</strong> Benedito, London Jr. e Bretas (2009),<strong>de</strong>senvolveu um método que será estudado nesse trabalho, o método parai<strong>de</strong>ntificação <strong>de</strong> ilhas observáveis, através do conceito <strong>de</strong> caminhos <strong>de</strong> fatoração,associado à fatoração da matriz Jacobiana. Esse método possibilita a i<strong>de</strong>ntificação<strong>de</strong> ilhas observáveis <strong>de</strong> uma forma bastante direta, s<strong>em</strong> exigir a solução <strong>de</strong>equações algébricas e apresenta baixo esforço computacional [1] .Com tantos estudos realizados e a constante necessida<strong>de</strong> <strong>de</strong> melhoria <strong>de</strong>conhecimentos e técnicas no assunto, que precisavam acompanhar a evolução ecrescimento dos sist<strong>em</strong>as elétricos, atualmente a análise <strong>de</strong> observabilida<strong>de</strong> atingiuo ponto <strong>em</strong> que po<strong>de</strong> ser feita para sist<strong>em</strong>as <strong>de</strong> alto grau <strong>de</strong> complexida<strong>de</strong> e queenvolv<strong>em</strong> inúmeras barras.Dentre os métodos atuais, os métodos numéricos costumam ser os maissimples por visar rotinas já disponíveis nos programas <strong>de</strong> estimadores <strong>de</strong> estado,exig<strong>em</strong> menos esforço computacional que os métodos topológicos. Esses, por suavez, fundamentam-se na teoria <strong>de</strong> grafos e apresentam a necessida<strong>de</strong> da criação <strong>de</strong>rotinas específicas, exig<strong>em</strong> o maior esforço computacional por ter<strong>em</strong> naturezacombinatória complexa com um grau diretamente proporcional ao tamanho enúmero <strong>de</strong> barras e ramos do SEP (Benedito, 2009).


10A partir da <strong>de</strong>finição <strong>de</strong> forma específica dos probl<strong>em</strong>as alvo da proposta e <strong>de</strong>seus efeitos <strong>de</strong> modo geral, po<strong>de</strong>-se elaborar as pr<strong>em</strong>issas e suposições relativasao resultado esperado do estudo proposto. Nessa fase inicial do trabalho, projetamseresultados relativos primeiramente ao aprendizado e manuseio <strong>de</strong> métodos jáconhecidos para a verificação da observabilida<strong>de</strong> e i<strong>de</strong>ntificação <strong>de</strong> ilhasobserváveis e medidas críticas. Com isso espera-se obter o conhecimento dosmétodos e <strong>de</strong> ferramentas comuns entre eles. A partir disso, a expectativa volta-separa o <strong>de</strong>senvolvimento <strong>de</strong> um software que possa <strong>de</strong>senvolver com o métodonumérico da fatoração da matriz jacobiana, um novo método para a mesmautilização. Desse, a hipótese que se levanta é fazê-lo capaz <strong>de</strong> simular sist<strong>em</strong>as <strong>de</strong>potência <strong>de</strong> testes <strong>de</strong> 14 a 30 barras do IEEE <strong>de</strong> modo a obter a confiabilida<strong>de</strong> daanálise, e a restauração da observabilida<strong>de</strong>, caso necessário, que permaneçamesmo com a perda <strong>de</strong> 1 ou 2 medidas quaisquer ou, até mesmo, <strong>de</strong> uma UTRqualquer, que é o critério técnico mais aceito [2] .1.3 OBJETIVOS1.3.1 Objetivo GeralDesenvolver um software para análise <strong>de</strong> observabilida<strong>de</strong> <strong>em</strong> sist<strong>em</strong>as elétricos<strong>de</strong> potência.1.3.2 Objetivos EspecíficosVerificar se um sist<strong>em</strong>a é observável ou não;I<strong>de</strong>ntificar ilhas observáveis;I<strong>de</strong>ntificar medidas críticas do sist<strong>em</strong>a <strong>de</strong> medidas utilizado;Restaurar a observabilida<strong>de</strong> do sist<strong>em</strong>a <strong>de</strong> potência.1.4 JUSTIFICATIVAOs sist<strong>em</strong>as <strong>de</strong> potência mo<strong>de</strong>rnos representam os maiores e maiscomplexos sist<strong>em</strong>as já construídos pelo hom<strong>em</strong>, o que exige técnicascomputacionais cada vez mais avançadas e precisas para construir, manter e operaresses sist<strong>em</strong>as (STEVENSON, 1994).Um sist<strong>em</strong>a elétrico <strong>de</strong> potência t<strong>em</strong> como finalida<strong>de</strong> suprir as exigências dosseus consumidores, e para alcançar esse objetivo é indispensável continuida<strong>de</strong>,


11qualida<strong>de</strong> e economia no fornecimento da energia elétrica. Desta maneira, uma dasmetodologias utilizadas para realizar essa operação e controle é o processo <strong>de</strong>Estimação <strong>de</strong> Estado, no qual se <strong>de</strong>termina as variáveis <strong>de</strong> estado do SEP a partirum conjunto <strong>de</strong> medidas redundantes.Entretanto, é comum <strong>em</strong> sist<strong>em</strong>as elétricos não se obter medições suficientes,das variáveis <strong>de</strong> estado, para o controle do sist<strong>em</strong>a observado, <strong>de</strong>vido a probl<strong>em</strong>asnos barramentos ou nas linhas <strong>de</strong> transmissão. Muitas vezes são utilizadas pseudomedidas,baseadas <strong>em</strong> dados históricos, na tentativa <strong>de</strong> restaurar o funcionamentodo sist<strong>em</strong>a. Todavia, <strong>de</strong>pen<strong>de</strong>ndo da complexida<strong>de</strong> do SEP essas pseudo-medidasnão são suficientes para restauração do sist<strong>em</strong>a, sendo necessária uma anális<strong>em</strong>ais profunda do probl<strong>em</strong>a.Nesse contexto, po<strong>de</strong>-se notar o quão importante é <strong>de</strong>senvolver umaferramenta capaz <strong>de</strong> realizar esse controle <strong>de</strong> medidas, atuar na segurançapreventiva, mostrando medidas críticas e ainda restaurar a observabilida<strong>de</strong> <strong>de</strong> umsist<strong>em</strong>a <strong>de</strong> potência.1.5 PROCEDIMENTOS METODOLÓGICOSA primeira ação a ser realizada é a fundamentação teórica, esta adquiridacom a revisão bibliográfica <strong>de</strong> livros, artigos, dissertações e, eventualmente,consultas na literatura disponibilizada pela biblioteca digital do IEEE. Então, seráiniciada a fase <strong>de</strong> estudos mat<strong>em</strong>áticos, utilizando o software MATLAB para realizarcálculos necessários e, principalmente, a fatoração triangular (eliminação <strong>de</strong> Gauss)<strong>de</strong> matrizes. De modo geral, po<strong>de</strong>-se <strong>de</strong>stacar esta etapa <strong>de</strong>scrita anteriormente,pois o projeto é baseado, <strong>em</strong> sua essência, na análise da fatoração da matriz e seusfatores, i<strong>de</strong>ntificando a <strong>de</strong>pendência linear existente entre as diversas barras dosist<strong>em</strong>a.O fluxograma mostrado na Figura 1 foi elaborado para facilitar o entendimentodos procedimentos seguidos no programa para realizar a análise <strong>de</strong> observabilida<strong>de</strong>.


12Figura 1 – Fluxograma do Processo.Fonte: Própria.Depois <strong>de</strong> finalizado o programa, começará a etapa das simulações, as quaisserão feitas com sist<strong>em</strong>as-testes <strong>de</strong> 14 a 30 barras do IEEE. Concluindo essa fase,<strong>de</strong> maneira eficiente e satisfatória, é possível partir para a etapa final do projeto,relatar os resultados e evi<strong>de</strong>nciar as conclusões.1.6 ESTRUTURA DO TRABALHOO trabalho será organizado da seguinte maneira:Capítulo 1 – Introdução, contendo a <strong>de</strong>limitação do t<strong>em</strong>a, probl<strong>em</strong>as epr<strong>em</strong>issas, objetivos geral, objetivos específicos, justificativa e procedimentosmetodológicos.Capítulo 2 – Embasamento teórico sobre Estimação <strong>de</strong> Estado <strong>em</strong> sist<strong>em</strong>aselétricos <strong>de</strong> potência, evi<strong>de</strong>nciando conceitos específicos da análise <strong>de</strong>observabilida<strong>de</strong>.Capítulo 3 – Desenvolvimento do programa para análise <strong>de</strong> observabilida<strong>de</strong>.Utilizando o software MATLAB será <strong>de</strong>senvolvido o algoritmo que verifique se osist<strong>em</strong>a é observável, além disso, i<strong>de</strong>ntifique as ilhas observáveis e possíveismedidas críticas. As simulações serão feitas com sist<strong>em</strong>as-testes <strong>de</strong> 14 a 30 barrasdo IEEE.


13Capítulo 4 – Apresentação dos resultados das simulações. Serãoapresentadas às telas das simulações que mostram se o sist<strong>em</strong>a analisado éobservável como um todo ou não. Se existir<strong>em</strong> ilhas observáveis no sist<strong>em</strong>a essastelas mostrarão as ilhas e as barras que estão contidas nelas, b<strong>em</strong> como as barrasque possu<strong>em</strong> medidas críticas.Capítulo 5 – Conclusão.1.7 CRONOGRAMATabela 1 – Cronograma para realização do TCC.Ativida<strong>de</strong>s2012 2013 2014Dez Jan Fev Mar Abr Mai Jun Jul Ago Set Out Nov Dez Jan Fev MarFundamentaçãoTeóricaReuniõesDesenvolvimento daPropostaEntrega daPropostaRedação do TCC 1EmbasamentoMat<strong>em</strong>áticoDesenvolvimento doAlgoritmoImpl<strong>em</strong>entação doAlgoritmoSimulaçõesRedação do TCC 2Conclusão dasAtivida<strong>de</strong>sFonte: Própria.1.8 REFERÊNCIAS[1] BENEDITO, Raphael A. S.; London Jr., João B.A.; Bretas, Newton G. Análise<strong>de</strong> Observabilida<strong>de</strong> <strong>em</strong> Sist<strong>em</strong>as Elétricos <strong>de</strong> Potência: I<strong>de</strong>ntificação <strong>de</strong> IlhasObserváveis Através <strong>de</strong> Caminhos <strong>de</strong> Fatoração Associados à MatrizJacobiana. XI Simpósio <strong>de</strong> Especialistas <strong>em</strong> Planejamento da Operação eExpansão Elétrica. 2009.


14[2] VIGLIASSI, Marcos Paulo. Projeto <strong>de</strong> Sist<strong>em</strong>as <strong>de</strong> Medição Confiáveispara Efeito <strong>de</strong> Estimação <strong>de</strong> Estado Via Algoritmos Evolutivos e Matriz H t Δ .2009. Dissertação (Mestrado <strong>em</strong> Engenharia Elétrica) – Escola <strong>de</strong> Engenharia <strong>de</strong>São Carlos, Universida<strong>de</strong> <strong>de</strong> São Paulo, São Carlos, 2009.[3] BENEDITO, Raphael Augusto <strong>de</strong> Souza. Avaliação da Qualida<strong>de</strong> eRedundância das Medidas para Estimação <strong>de</strong> Estado <strong>em</strong> Sist<strong>em</strong>as Elétricos <strong>de</strong>Potência. 2007. 207 f. Dissertação (Mestrado <strong>em</strong> Engenharia Elétrica) – Escola <strong>de</strong>Engenharia <strong>de</strong> São Carlos, Universida<strong>de</strong> <strong>de</strong> São Paulo, São Carlos, 2007.BRETAS, N.G. (1996). Network Observability: A theory based on triangularfactorization and path graph concepts - part I. International conference onAutomatic Control, cida<strong>de</strong> do Porto, Portugal, Set<strong>em</strong>bro.STEVENSON, William D. Jr.; GRAINGER, John J. Power Syst<strong>em</strong> Analysis. 1 st ed.Singapore: Mc Graw-Hill. 1994.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!