Resumos
Atualmente, a ANEEL (Agência Nacional de Energia Elétrica) aplica penalidades nas empresas de transmissão devido aos desligamentos programados e não-programados. O uso destas penalidades tem como objetivo assegurar a confiabilidade do serviço de transmissão em um ambiente competitivo. Este artigo propõe um método para minimizar as penalidades associadas com desligamentos de equipamentos de transmissão através da otimização dos intervalos das manutenções programadas. O modelo de otimização proposto neste artigo se baseia em novo método analítico que é capaz de estimar precisamente os valores esperados das penalidades.
Confiabilidade; Manutenção; Modelos Markovianos; Simulação Monte Carlo; Otimização
Currently, the National Agency for Electric Energy applies penalties to the transmission companies due to scheduled and non-scheduled outages. These penalties have as objective to assure the reliability of the transmission service in a competitive environment. This paper proposes a method to minimize the penalties associated to the transmission equipment outages through the optimization of the scheduled maintenance interval. The proposed optimization model is based on a new analytic method that estimates the expected penalties accurately.
Reliability; Maintenance; Markovian Models; Monte Carlo Simulation; Optimization
SISTEMAS ELÉTRICOS DE POTÊNCIA
Determinação do intervalo de manutenção programada da proteção de linhas de transmissão considerando-se penalidades associadas à indisponibilidade
Determination of scheduled maintenance interval in the protection of transmission lines considering the penalties associated to transmission equipament unavailabilities
Felipe E. L. PereiraI; Ricardo B. PradaI; Albert C. G. de MeloII; Anselmo Barbosa RodriguesIII; Maria da Guia da SilvaIII
IDepartamento de Engenharia Elétrica, PUC-Rio, Rua Marquês de São Vicente, 225, CEP 22453-900, Rio de Janeiro, RJ. felipel@puc-rio.br; prada@ele.puc-rio.br
IICEPEL, Av. Hum, S/N, Ilha da Cidade Universitária, CEP 2194-590 - Rio de Janeiro, RJ. albert@cepel.br; nebulok_99@yahoo.com
IIIDepartamento de Engenharia Elétrica, UFMA, Av. dos Portugueses, S/N, CEP 65080-040, São Luís, MA. guia@dee.ufma.br
RESUMO
Atualmente, a ANEEL (Agência Nacional de Energia Elétrica) aplica penalidades nas empresas de transmissão devido aos desligamentos programados e não-programados. O uso destas penalidades tem como objetivo assegurar a confiabilidade do serviço de transmissão em um ambiente competitivo. Este artigo propõe um método para minimizar as penalidades associadas com desligamentos de equipamentos de transmissão através da otimização dos intervalos das manutenções programadas. O modelo de otimização proposto neste artigo se baseia em novo método analítico que é capaz de estimar precisamente os valores esperados das penalidades.
Palavras chaves: Confiabilidade, Manutenção, Modelos Markovianos, Simulação Monte Carlo, Otimização.
RESUMO
Currently, the National Agency for Electric Energy applies penalties to the transmission companies due to scheduled and non-scheduled outages. These penalties have as objective to assure the reliability of the transmission service in a competitive environment. This paper proposes a method to minimize the penalties associated to the transmission equipment outages through the optimization of the scheduled maintenance interval. The proposed optimization model is based on a new analytic method that estimates the expected penalties accurately.
Keywords: Reliability, Maintenance, Markovian Models, Monte Carlo Simulation, Optimization.
1 INTRODUÇÃO
O modelo econômico do setor elétrico brasileiro baseia-se na desagregação entre os agentes de geração, transmissão, distribuição e comercialização de energia. A separação das atividades de geração e transporte de energia elétrica proporcionou a realização de leilões de concessões de linhas de transmissão e a participação de empresas transmissoras ou consórcios nestes leilões. As concessões das linhas de transmissão são dadas para as empresas de transmissão ou consórcios, que ofertarem as menores tarifas de transmissão (menores receitas anuais), desde que os valores ofertados pelos demais proponentes, sejam superiores a um determinado percentual estabelecido para o leilão. Caso seja inferior ao valor do percentual estipulado, o leilão continua sendo realizado através de lances sucessivos por meio de deságios. Devido às suas naturezas, os investimentos em linhas de transmissão são caracterizados por custos de implantação relativamente baixos e por uma receita operacional pré-estabelecida, de acordo com o lance vencedor do leilão, satisfazendo os limites máximos estipulados pela ANEEL. As empresas de transmissão ou consórcios vencedores firmam um Contrato de Prestação de Serviço de Transmissão (CPST) com o Operador Nacional do Sistema Elétrico (ONS), delegando-lhe o direito de administrar e coordenar o uso da rede básica e, em contrapartida, tem garantido o ressarcimento de seus custos e de seus investimentos remunerados. Os investimentos podem ser considerados de baixo risco, já que a receita associada é garantida através do CPST firmado. Consequentemente, o modelo econômico que norteia a prestação de serviços de transmissão de energia elétrica é fundamentado em uma estrutura reguladora. Ao contrário dos geradores e comercializadores, que têm liberdade para estabelecer seu preço no ambiente competitivo, a renda proveniente da transmissão é estabelecida pela ANEEL, como uma Receita Anual Permitida (RAP) pelo "aluguel" de seus ativos.
Os usuários do sistema (geradores, distribuidores e grandes consumidores) arcam com os investimentos necessários para se conectarem à rede básica e pagam ao ONS uma tarifa pelo uso do sistema. Esta tarifa é fixada de forma a cobrir o conjunto dos contratos do ONS com as empresas de transmissão ou consórcios, mais os custos da operação do sistema.
Uma das atribuições do agente regulador é determinar o desconto da RAP de uma empresa de transmissão em função da indisponibilidade de operação de seus equipamentos. Este desconto é denominado Parcela Variável Referente à Disponibilidade de Instalações da Rede Básica, ou simplesmente Parcela Variável (PV) e tem sido tema de grandes debates entre os agentes transmissores, a ANEEL, o ONS e instituições financeiras que têm investido em empreendimentos de expansão do setor elétrico nacional. As indisponibilidades dos ativos de transmissão podem ocorrer de duas maneiras: de forma voluntária ou programada (manutenções, conveniência operativa, manobras, e outras) e por saídas forçadas ou não programadas (falhas). Neste sentido, a fórmula estabelecida pela ANEEL (2007) nas licitações, para calcular a PV tem a forma geral:
onde:
PV = parcela variável por indisponibilidades de equipamentos;
PB = parcela equivalente ao duodécimo da RAP, associada à plena disponibilidade dos ativos que compõem uma instalação de uma empresa de transmissão de energia;
g (DDPi) = KpDDPi
h (DODi) =
Kp fator para desligamentos programados, sendo Kp = 10;
Ko fator para desligamentos não programados, sendo Ko = 150;
NP número de desligamentos programados da instalação durante o mês;
NO número de desligamentos não programados da instalação ao longo do mês;
D número de dias do mês;
DDPi duração em minutos de cada desligamento programado que ocorra no mês;
DODi duração em minutos de cada desligamento não programado que ocorra no mês.
A partir de (1) percebe-se que a PV é uma variável aleatória, ou seja, depende das indisponibilidades programadas (manutenção preventiva) e não programadas. Por sua vez, estas últimas, dependem das durações das falhas dos ativos de transmissão, a saber: falhas permanentes (queda de torres, falha de equipamento, falha humana, etc.), onde é necessária a intervenção da equipe de manutenção, e falhas transitórias (descarga atmosférica, queimadas), onde o ativo retorna à operação pela atuação de religamento manual ou automático. Neste contexto, torna-se evidente que as empresas de transmissão devem reduzir o custo decorrente da PV, de modo a maximizar o seu lucro anual. Uma das principais variáveis associadas com a minimização da PV é o intervalo de manutenção. Contudo, esta variável tem efeitos conflitantes sobre a PV. As atividades de manutenção melhoram a condição dos equipamentos. Consequentemente é possível minimizar as penalidades devido aos desligamentos não-programados reduzindo os intervalos de manutenção. Por outro lado, o aumento na frequência de manutenções causa uma elevação nas penalidades devido aos desligamentos programados. Devido a isto, as empresas de transmissão devem determinar os intervalos de manutenção visando obter um equilíbrio entre as penalidades causadas por desligamentos programados e não programados. Desta forma, o principal objetivo deste artigo é apresentar um método que minimize a PV através do ajuste do intervalo de manutenção das proteções das linhas de transmissão. A principal característica do modelo de otimização proposto neste artigo é o estabelecimento de um compromisso aceitável entre as parcelas da PV associadas com desligamentos programados e não-programados. Adicionalmente, o modelo de otimização da PV utiliza um método analítico que pode estimar precisamente o valor esperado da PV. Os resultados dos testes demonstraram que as soluções geradas pelo algoritmo de minimização da PV são melhores que aquelas obtidas via maximização da disponibilidade (Anders, 1990). Além disso, as estimativas da PV mostradas neste artigo são mais precisas que aquelas calculadas com base nos método analíticos de estimação do custo de interrupção (Billinton & Allan, 1996). O restante do artigo é organizado da seguinte forma:
i) na Seção 2 apresenta-se o modelo Markoviano usado para modelar a indisponibilidade das linhas de transmissão;
ii) o método analítico de estimação da PV é apresentado na Seção 3;
iii) na Seção 4 descreve-se o modelo de otimização da PV; e
iv) os resultados e as conclusões são apresentados nas Seções 5 e 6.
2 MODELO MARKOVIANO PARA EQUIPAMENTOS DE PROTEÇÃO DE LINHAS DE TRANSMISSÃO
Os modos de falha de equipamentos de proteção do sistema elétrico foram definidos por Siqueira (1999). A partir da análise dos modos de falhas, é possível identificar o conjunto dos estados possíveis de um sistema de proteção. Inicialmente, três estados são identificados, quando a proteção estiver disponível para a operação:
i) Normal: se a proteção estiver apta a executar sua missão, atuando conforme projetada. Em outras palavras, não apresentar defeito interno ou intrínseco, estando em condições de desligar o equipamento protegido na ocorrência de um defeito interno do mesmo e, simultaneamente, bloquear seu desligamento na ausência de defeito ou na ocorrência de defeitos externos ao processo.
ii) Falha causada por desligamento indevido do sistema de proteção: quando a proteção desligar intempestivamente o equipamento protegido, na ausência de uma falha no equipamento protegido, necessitando de reparo para voltar a operar (manutenção corretiva).
iii) Falha causada por recusa de atuação do sistema de proteção: é o caso onde a proteção apresenta uma falha que a impeça de atuar na ocorrência de uma falha no equipamento protegido. A proteção voltará ao estado normal após a realização de manutenção corretiva. Neste caso, o equipamento protegido será desligado pelas proteções de retaguarda de outros equipamentos, ou pelos sistemas de proteção contra falha de disjuntor, em consequência da recusa de atuação da proteção. Neste caso, o equipamento também exige a realização de uma manutenção corretiva.
iv) Teste: quando a proteção estiver sendo submetida a uma inspeção ou teste periódico, manutenção preventiva ou auto-teste.
O diagrama de Markov usado para representar as transições entre os quatro estados de operação definidos é mostrado na Figura 1 (Pereira, 2009), onde:
λ 12 taxa de transição entre os estados normal e teste;
µ21 taxa de transição entre os estados teste e normal;
λ 13 taxa de falha do estado normal para falha por desligamento indevido;
µ31 taxa de reparo do estado desligamento indevido para normal;
λ 14 taxa de falha do estado normal para falha causada por recusa;
µ41 taxa de reparo do estado falha causada por recusa para normal.
No sistema em estudo, agrupam-se os estados de falha devido à recusa de atuação do sistema de proteção e falha devido a desligamentos em um único estado não-conforme, doravante denominado de falha. Para manter a ocorrência lógica, o agrupamento só será viável se o novo modelo permanecer um sistema Markoviano. Para isto será suficiente que as taxas de transição do estado agregado para qualquer outro estado, permaneçam iguais às da cadeia original (Howard, 1971; Endrenyi, 1978). Os valores das probabilidades calculadas após agrupamento dos estados só são precisos em regime permanente, ou seja, o agrupamento de estados não pode ser usado quando se deseja calcular valores transitórios de probabilidade (Endrenyi, 1978). A Figura 2 (Pereira, 2009), ilustra o diagrama de transição entre estados do processo de Markov resultante, no qual o estado agregado é representado pelo número cinco, os demais estados permanecem com a numeração original.
Em geral, a taxa de transição do estado NORMAL para o estado falhado (λ 15), não é observável diretamente, uma vez que os defeitos internos à proteção são invisíveis ao operador (Siqueira, 1999). Em consequência, a taxa de falha de proteção de linhas de transmissão (λ 15 = λ f), principal variável da confiabilidade do sistema de proteção, terá de ser calculada através de outras variáveis observáveis, usando as equações do modelo. As variáveis observáveis são ocorrências e durações de eventos visíveis, tais como: frequências de recusas de atuação e atuações indevidas e os tempos médios de teste e de reparo. Estes parâmetros foram calculados em estudo realizado por Siqueira (1999) para um período de observação de onze anos. Portanto, a taxa de falha do sistema de proteção é determinada por:
onde:
Fr é a frequência de recusas de atuação;
Fe é a frequência de atuações indevidas;
m = 1/µ21 tempo médio de teste;
r = 1/µ51 tempo médio de reparo;
Ft é a periodicidade da manutenção.
Em (2), todos os elementos do lado direito podem ser avaliados por estatísticas dos dados históricos de desligamentos forçados. Fr e Fe representam, respectivamente, as frequências observadas de recusas e atuações indevidas no período, ambas observáveis, após um tempo longo de adoção da frequência de teste Ft. Com exceção destas frequências, os demais dados são característicos do sistema de proteção, e independentes da frequência de testes Ft. Determinam-se as probabilidades em regime permanente dos estados normal, teste e falha através da matriz de transição do diagrama de estados (Billinton & Allan, 1992) ilustrado na Figura 2:
A partir de (3) e (4) tem-se :
Finalmente, deve-se mencionar que neste artigo foi considerado que as frequências de atuações indevidas (Fe) e de recusas (Fr) de atuação são correlacionadas com os intervalos de manutenção (Ft-1) através de:
onde:
Feº(Frº) e δ e(δr) são o valor inicial e o fator de crescimento, respectivamente, associados com a frequência de atuações indevidas (recusas de atuação).
Neste artigo considera-se que a indisponibilidades das linhas de transmissão são causadas apenas pelas falhas e ações de manutenção nos seus equipamentos de proteção. Contudo, pode-se expandir o modelo Markoviano deste artigo para incluir falhas de outra natureza.
3 ESTIMAÇÃO DA PV USANDO O MODELO MARKOVIANO PROPOSTO
Nesta seção são apresentados os métodos analíticos usados na estimação do valor esperado da PV. Em 3.1 descreve-se uma técnica inspirada no cálculo de custos de interrupção (Billinton & Allan, 1996) que é denominada de Método Analítico Aproximado. Em 3.2 descreve-se uma técnica mais precisa que relaxa algumas aproximações usadas no modelo convencional. Esta técnica é denominada Método Analítico Exato.
3.1 MÉTODO ANALÍTICO APROXIMADO
A PV está associada com duas componentes: desligamentos programados e não programados. Ambas as componentes podem ser estimadas considerando-se que:
i) a duração de todos os desligamentos é igual a duração média dos desligamentos.
ii) O número de desligamentos é igual ao valor esperado da frequência dos desligamentos.
iii) Os desligamentos programados são associados com a duração e frequência do estado 2 na Figura 2.
iv) Os desligamentos não-programados estão associados com a duração e a frequência do estado 5 na Figura 2.
As considerações (i)-(iv) são válidas somente quando há pouca dispersão em torno dos valores esperados das durações dos desligamentos. Caso contrário, as penalidades devem ser estimadas usando-se o Método de Simulação Monte Carlo (Billinton & Allan, 1992; Billinton & Li, 1994; Billinton & Allan, 1996).
A partir das considerações (i), (ii) e (iii) é possível estimar o valor esperado da parcela da PV associada com desligamentos programados (E[PVp]) da forma:
onde D2 e F2 são a duração (em meses) e a frequência (em ocorrências por mês) do estado 2 na Figura 2.
Por outro lado, o valor esperado da PV referente aos desligamentos não-programados (E[PVo]) é calculado usandose as considerações (i), (ii) e (iv). Desta forma, tem-se que E[PVo] é:
a) para D5< 300 min.
E[PVº] PB ( Ko.F5.D5 .1440D)
b) para D5< 300 min.
onde:
D5 e F5 são a duração (em meses) e a frequência (em ocorrências por mês) do estado 5 na Figura 2,
O valor esperado da PV ( E[PV]. ) é:
Neste ponto, deve-se mencionar que os princípios usados para calcular o E[PV]. , através de (13), são os mesmos aplicados na estimação dos custos de interrupção do fornecimento de energia em redes de distribuição (Billinton & Allan, 1996). Os custos de interrupção também podem ser estimados usando-se os valores esperados das frequências de falhas e dos tempos de restauração nos pontos de carga. Devido a isto, as estimativas do E[PV]. e dos custos de interrupção estão sujeitos à mesma restrição: as variações em torno dos tempos médios de restauração e de desligamentos (programados e não-programados) devem ser pequenas.
3.2 MÉTODO ANALÍTICO EXATO
Os valores esperados da PV associados com desligamentos programados e não programados podem ser definidos como:
Os valores esperados mostrados em (14) e (15) são na verdade valores esperados de uma soma de variáveis aleatórias, ou seja (Miller & Childers, 2004):
onde:
µx é o valor esperado de Xi;
E[N] é o valor esperado da variável aleatória inteira N.
Xi com i = 1, N é conjunto de variáveis aleatórias com distribuições idênticas.
A partir de (16) tem-se que E[PVp] e E[PVo]
são:
Os números esperados de desligamentos programados e não programados podem ser obtidos calculando-se as frequências dos estados de teste e falha do modelo de Markov da Figura 2:
Por outro lado, os valores esperados das funções g(DDPi) e h(DODi) estão associados com o cálculo do valor esperado de funções de variáveis aleatórias:
onde: fX (x) é a função densidade de probabilidade da variável aleatória x e ρ (x) é a função da variável aleatória x e ρ(x) é a função da variável aleatória x.
Consequentemente, para estimar os valores esperados das funções g(DDPi) e h(DODi) deve-se conhecer as funções densidade de probabilidade associadas com as durações dos estados de falha e teste. Contudo, em uma cadeia de Markov homogênea, o tempo gasto em um estado tem distribuição exponencial com taxa de falha igual a soma de todas as taxas de transição que saem deste estado (Gubner, 2006). Devido a isto, tem-se que as funções densidade de probabilidade dos estados de falha (estado 5) e teste (estado 2) são dadas por:
onde f2 (t) e f5 (t) são as funções densidade de probabilidade associadas, respectivamente, com as durações dos estados 2 e 5 da Figura 2,.
Desta forma, os valores esperados das funções g(DDPi e h(DODi são:
Substituindo-se (19) e (21) em (17) encontra-se que a penalidade esperada associada com os desligamentos programados é:
Substituindo-se (20) e (22) em (18) tem-se que a penalidade esperada associada com os desligamentos não-programados é:
Finalmente, tem-se a expressão para a estimação do E[PV].
É importante observar que as expressões de E[PVp]
obtidas pelos métodos analíticos aproximado (10) e exato (23) são idênticas. As expressões de E[PVo]
geradas pelos métodos aproximado (11) e (12) e exato (24) são diferentes. Esta diferença é devido ao fato de que o método analítico aproximado estima E [PVo]
considerando que o valor esperado da função ρ (x) de uma variável aleatória x é igual a ρ (µX), onde µX é o valor esperado da variável aleatória. Esta aproximação não é válida para funções lineares por partes, tais como a parcela da PV devido aos desligamentos não programados (ver (22) ).
4 MODELOS DE OTIMIZAÇÃO USADOS PARA DETERMINAR O INTERVALO DE MANUTENÇÃO
Nesta seção são apresentados dois modelos de otimização usados para determinar os intervalos de manutenção das proteções das linhas de transmissão. O primeiro modelo se baseia na maximização da disponibilidade. O segundo modelo, que é uma das contribuições deste artigo, determina o intervalo de manutenção através da minimização do valor esperado da PV.
4.1 MODELO BASEADO NA MAXIMIZAÇÃO DA DISPONIBILIDADE
Geralmente os intervalos de manutenção são otimizados visando maximizar a disponibilidade de um equipamento (Anders, 1990). No caso do modelo Markoviano da Figura 2, a maximização da disponibilidade está associada com o problema de otimização:
maximizar
Existe apenas uma variável de decisão no problema de otimização: a frequência de manutenções programadas Ft (o inverso do intervalo de manutenção). Não é possível obter uma solução algébrica para o problema (26) a partir das condições de optimalidade, ou seja:= 0. Esta restrição é devido a presença dos termos exponenciais associados com λ f. Apesar disso, é possível resolver (26) usando-se métodos de otimização irrestrita, tais como: Busca Dicotômica, Busca de Fibonacci, Busca Áurea, Método de Interpolação Quadrática e Método de Interpolação Cúbica (Antoniou & Lu, 2007). Nestes métodos, a solução ótima é obtida baseandose somente no cálculo da função objetivo. Entretanto, não há garantia de que uma solução ótima de (26) seja um mínimo local do E[PV]. . Em outras palavras, não há garantia de que a maximização de P1 seja um problema de otimização equivalente ("proxy") a minimização do E[PV]. .
4.2 MODELO BASEADO NA MINIMIZAÇÃO DO VALOR ESPERADO DA PV
No panorama atual do setor elétrico nacional é mais adequado para as empresas de transmissão minimizar o E[PV]. do que maximizar a disponibilidade dos equipamentos. Este objetivo pode ser alcançado através da solução do problema de otimização:
onde:
O problema de otimização (27) tem as seguintes características de otimização em comum com (26): i) ambos são problemas de otimização não-linear irrestrita; ii) a única variável de decisão é a frequência de manutenções programadas Ft; iii) não é possível obter soluções algébricas para os dois problemas devido aos termos exponenciais associados com λ f. Há termos exponenciais na expressão de E [PVo] . iv) a solução de (27) pode ser obtida aplicando-se os mesmos métodos usados para resolver (26).
Ao contrário da solução de (26), é possível garantir que uma solução ótima de (27) seja um mínimo local da PV.
5 RESULTADOS
Nesta seção são apresentados os resultados obtidos com a aplicação dos modelos e técnicas propostos para estimar e otimizar o valor esperado da PV. Os resultados apresentados nas próximas subseções são organizados como:
i) Em 5.1 apresenta-se um estudo comparativo entre os métodos analíticos aproximado e exato para a estimação do valor esperado da PV; ii) Os intervalos de manutenção obtidos com os modelos de otimização descritos na Seção 4 são mostrados em 5.2. Além disso, são mostradas diversas análises de sensibilidade do valor esperado da PV com relação ao intervalo de manutenção.
5.1 ESTIMAÇÃO DO VALOR ESPERADO DA PV
Nesta seção é apresentada uma análise comparativa dos métodos analíticos usados para estimar o valor esperado da PV, isto é: os métodos aproximado e exato descritos nas Subseções 3.1 e 3.2, respectivamente. As estimativas do E[PV]. geradas por estes dois métodos serão comparadas com aquelas obtidas com a Simulação Monte Carlo (Billinton & Allan, 1992; Billinton & Li, 1994; Billinton & Allan, 1996). A Simulação Monte Carlo (SMC) é utilizada como parâmetro para avaliar a precisão das estimativas geradas pelos dois métodos analíticos, ou seja: o método analítico que produzir as estimativas mais próximas daquelas obtidas com a simulação será considerado o mais preciso. Esta comparação é realizada considerando-se que:
i) o tempo médio de duração do teste (m) é igual a 21,7490 horas (3, 0207 x 10-2 meses); ii) o fator de crescimento (δ e) e o valor inicial (Feº) da frequência de atuações indevidas são iguais a 0,01 e 0,03445, respectivamente; iii) o fator de crescimento (δ r) e o valor inicial () da frequência de recusas de atuação são iguais a 0,001 e 0,01867, respectivamente; iv) o intervalo de manutenção (Ft) é de 36 meses; v) o tempo médio de reparo (r) é variado de 1 até 20 horas; vi) a RAP é igual a R$ 25.800.000,00. Consequentemente, o pagamento base (PB) é igual R$ 2.150.000,00. vii) o tamanho da amostra usado para estimar a PV via simulação foi de 10.000 meses.
A estimação da PV está associada com os seguintes aspectos: i) frequência e a duração dos estados de teste e falha; ii) função não-linear da duração do estado de falha. Devido a isto, é mais adequado estimar a PV via SMC considerando o processo cronológico de transição de estados do equipamento. Neste artigo, o processo cronológico de operação do equipamento é simulado usando-se a amostragem de duração de estados (Billinton & Li, 1994). O sorteio das durações dos estados do equipamento foi realizado considerando-se que as durações têm distribuição exponencial. O algoritmo conceitual para estimar a PV via SMC é mostrado no Apêndice Apêndice .
Os E[PV]. estimados pelos dois métodos analíticos e pela simulação são mostrados na Figura 3. Pode-se concluir que as estimativas do E[PV]. geradas pelo método analítico exato são mais próximas daquelas obtidas via simulação Monte Carlo do que as estimativas calculadas pelo método analítico aproximado. Este fato é um indicativo de que o método analítico proposto neste artigo é mais robusto do que os métodos baseados no cálculo dos custos de interrupção para estimar o E[PV].. A precisão dos métodos analíticos com relação à simulação pode ser avaliada de forma mais precisa calculando-se os erros relativos das estimativas do E[PV]. obtidas com os métodos analíticos com relação aquelas geradas via simulação. O gráfico dos erros relativos entre os métodos analíticos aproximado e exato é mostrado na Figura 4.
O valor máximo do erro associado com o método analítico aproximado é igual a 50.7287%. Por outro lado, o valor máximo do erro associado com o método analítico exato é apenas 8.5239%. Além disso, o erro absoluto médio associado com os métodos analíticos aproximado e exato são iguais a 20.5925% e 3.9913%, respectivamente. Estes resultados demonstram quantitativamente que o método analítico exato proposto neste artigo é a técnica mais adequada para estimar o valor esperado da PV. É importante enfatizar que a simulação Monte Carlo é um método baseado na amostragem aleatória. Portanto, os resultados gerados pela simulação serão diferentes para cada sorteio de uma nova amostra. Contudo, os resultados ficam restritos a um intervalo definido por um nível de confiança especificado. Este intervalo contém a média populacional do índice estimado pela simulação através de uma média amostral. Consequentemente, espera-se que os valores estimados pelos métodos analíticos estejam contidos dentro do intervalo de confiança gerado pela simulação. Caso contrário, o valor esperado estimado por um modelo analítico não estará associado com a média amostral gerada pela simulação, isto é, o modelo analítico não é um bom estimador para o índice probabilístico desejado. Os intervalos de confiança associados com os valores das penalidades estimados via simulação, com o tempo de reparo variando de 1 até 20 horas, são mostrados na Figura 5. Também são mostradas nesta figura as penalidades estimadas pelos métodos aproximado e exato. Pode-se observar que os valores das penalidades estimados via método analítico exato estão contidos dentro dos intervalos de confiança gerados pela simulação. Em outras palavras, o método analítico exato é um bom estimador da média populacional das penalidades. A maior parte das penalidades calculadas pelo método analítico aproximado está fora do intervalo de confiança gerado pela simulação. Desta forma, a simulação não reconhece as estimativas geradas pelo método aproximado como uma média populacional das penalidades.
5.2 DETERMINAÇÃO DOS INTERVALOS DE MANUTENÇÃO VIA MODELOS DE OTIMIZAÇÃO
Nesta seção são apresentados os resultados obtidos com os dois modelos de otimização do intervalo de manutenção, ou seja: i) maximização da disponibilidade; ii) minimização do E[PV] .
Estes resultados foram obtidos para as mesmas condições definidas na Subseção 5.1. O tempo médio de reparo (r) foi fixado em 4,3861 horas (6,0918x10-3 meses) e a frequência de manutenções programadas (Ft) foi determinada pelos modelos de otimização citados.
Inicialmente, são analisados os efeitos das variações no intervalo de manutenção sobre a disponibilidade (probabilidade do estado normal). Estes efeitos são mostrados na Figura 6. Pode-se observar que a disponibilidade é uma função convexa do intervalo de manutenção no domínio [1,100]. Desta forma, é possível calcular o valor ótimo global desta probabilidade no intervalo [1,100]. O valor ótimo do intervalo de manutenção e da disponibilidade são iguais a 50,0497 meses e 0,9990, respectivamente. A maximização da disponibilidade foi realizada através da função fminbnd do toolbox de otimização do MATLAB (Branch, 1996). O método "default" para solução de problemas de otimização irrestrita na função fminbnd é Busca Áurea com interpolação polinomial (Antoniou & Lu, 2007).
Em seguida, foram analisadas as variações no E[PV] e nas suas componentes (E[ PVp] e E[ PVo] ) com relação ao intervalo de manutenção. As variações de E[ PVp] e E[ PVo] com relação ao intervalo de manutenção são mostradas na Figura 7. É importante citar que os valores de E[ PVp] e E [PVo] mostrados nesta figura foram obtidas usando-se o método analítico exato, pois este método foi identificado como o mais preciso para estimar o E[PV] nos testes realizados na Subseção 5.1.
A partir da Figura 7, Pode-se concluir que E[PVp] aumenta quando o intervalo de manutenção é diminuído. Por outro lado, pode-se observar que E[ PVo] é reduzido quando há um aumento na frequência de manutenções. Além disso, as variações no E[PV] (E [PVp] + E [PVo] ) com relação ao intervalo de manutenção são mostradas na Figura 8. Verificase que o E[PV] é uma função convexa do intervalo de manutenção no domínio [1,100]. Devido a isto, existe um valor ótimo global para o E[PV] no intervalo [1,100]. Os valores ótimos do E[PV] e do intervalo de manutenção obtidos com a função fminbndsão iguais a R$ 60.722,1707 e 24,0145 meses, respectivamente. Estes valores são identificados através das retas tracejadas na Figura 8. É importante observar que o intervalo de manutenção obtido com a maximização da probabilidade do estado normal (reta vertical pontilhada na Figura 8) é maior do que aquele obtido com a minimização do E[PV] . Em outras palavras, a maximização da disponibilidade resulta em um intervalo de manutenção subótimo com relação à minimização do E[PV]. Por exemplo, o E[PV] obtido através da solução do problema (26) (R$ 93.772,8275) é 35,2455% maior do que o E[PV]. associado com o problema (27) (R$ 60.722,1707).
Finalmente, é importante avaliar os percentuais associados com as parcelas do E[PV] obtidos com os modelos de maximização da disponibilidade e minimização do E[PV] . Estes percentuais são mostrados nas Figuras 9 e 10. Observa-se que a solução ótima do problema de minimização do E[PV] produz uma divisão mais equitativa da E[PV] , entre as suas componentes (E[PVp] e E[PVo]), do que a maximização da disponibilidade.
6 CONCLUSÕES
Neste artigo foi apresentado um modelo probabilístico para o sistema de proteção de linhas de transmissão que incorpora o efeito dos intervalos de manutenção nas taxas de falha do modelo. Os modelos probabilísticos do sistema de proteção foram usados para otimizar os intervalos de manutenção de acordo com dois paradigmas: maximização da disponibilidade do equipamento e minimização do valor esperado da Parcela Variável ( E[PV] ). Os testes realizados demonstraram que o último paradigma é mais adequado à realidade do setor elétrico nacional. Esta adequação é devido ao fato de que as penalidades estimadas com a minimização do E[PV] são menores do que aquelas calculadas com a maximização da disponibilidade dos equipamentos. Em ambos os paradigmas citados, foi possível otimizar os intervalos de manutenção usando-se apenas algoritmos para a solução de problemas de programação não-linear irrestrita unidimensional. Adicionalmente, foi proposto um método analítico exato para estimar o E[PV] . Os E[PV] calculados por este método são mais precisos do que aqueles obtidos com aproximações tradicionalmente usadas na estimação probabilística dos custos de interrupção.
Os custos associados com as manutenções programadas não são considerados no método de determinação do intervalo de manutenção proposto neste artigo. Isto é, a periodicidade da manutenção é determinada considerando somente as penalidades devido aos desligamentos programados e não programados. Contudo, os custos de manutenção (mão de obra, substituição de equipamentos, inspeção, etc.) podem ser adicionados diretamente à função objetivo (27) para minimizar simultaneamente as penalidades e os custos de manutenção. Também é possível incluir os custos de reparo do equipamento na função objetivo (Endrenyi, 1998). Esta estratégia de modelagem dos custos de manutenção tende a causar um aumento nos intervalos de manutenção, pois a componente associada com os custos de manutenção (frequência de manutenção vezes o custo de manutenção unitário) funciona como um termo de penalidade na função objetivo. Os aspectos discutidos acima serão abordados pelos autores em futuras publicações.
Artigo submetido em 14/05/2010 (Id.: 01149)
Revisado em 07/08/2010, 21/10/2010
Aceito sob recomendação do Editor Associado Prof. Julio Cesar Stacchini Souza
APÊNDICE
O algoritmo conceitual da SMC usada para estimar a PV é mostrado abaixo:
i) Faça PVest = 0, onde PVest é o valor estimado da PV.
ii) Repita os passos (iii)-(ix) para i = 1,...,NS, onde i é o contador de simulações e NS é tamanho da amostra de meses que serão sorteados (10.000).
iii) Faça tsim = 0 e PVsim = 0, onde tsim é o tempo de simulação PVsim e é o valor da PV para o i-ésimo mês simulado.
iv) Faça s = 0, onde s = 1, 2 e 5 indica que o equipamento está operando nos estados normal, teste ou falha, respectivamente.
v) Repita os passos (vi)-(viii) enquanto tsim < 1 mês.
vi) Sorteio do tempo de residência no estado atual e atualização da PV:
vi-a) Se s = 1, então t12 = - (λ 12) log (Z1) e t15 = - (λ 15) log (Z2) , onde Z1 e Z2 são dois números aleatórios com distribuição uniforme. Se min {t12 , t15 } = t12 então dsim=t12 e snew = 2, onde dsim é a duração do estado atual e s new é o próximo estado de operação do equipamento. Caso contrário (min{t12, t15} = t15), dsim = t15 e snew = 5).
vi-b) Se s = 2, então: dsim= - (1/µ21) log (Z1); snew1 e PVsim= PVsim + PB x Kp x dsim
vi-c) Se s = 5, então: dsim = - (1/µ51) log (Z1); e snew= 1 e atualize a PV de acordo com a seguinte equação:
vii)Atualize o tempo de simulação:
viii)Atualize o estado do equipamento: s = snew
ix)Atualize o valor estimado da PV: PVest = PVest + PVsim
x)Calcule o valor estimado da PV: PVest =/N S
Referências bibliográficas
- Agência Nacional de Energia Elétrica (2007) Resolução Normativa N 270 de 26 de junho de 2007 (?) , Disponível de: <http://www http://www.aneel.gov.br> [Acessado em 22 de março de 2010]
- Antoniou, A. & Lu, W., (2007) Practical Optimization: Algorithms and Engineering Applications. Springer.
- Anders, G.J. (1990) Probability Concepts in Electric Power Systems. John Wiley & Sons.
- Billinton, R. & Allan, R.N. (1992) Reliability Evaluation of Engineering Systems: Concepts and Techiniques. Plenum Press.
- Billinton, R. & Allan, R. N. (1996) Reliability Evaluation of Power Systems. 2nd ed. Plenum Press.
- Billinton, R. & Li, W. (1994) Reliability Assessment of Electric Power Systems Using Monte Carlo Methods. Plenum Press.
- Branch, M. A. (1996) MATLAB Optimization Toolbox: User's Guide. Math Works, Inc.
- Endrenyi, J., (1978) Reliability Modeling in Electric Power Systems, John Wiley & Sons.
- Endrenyi, J., Anders, G. J., & Leite da Silva, A. M. (1998), Probabilistic Evaluation of the Effect of Maintenance on Reliability - An Application. IEEE Trans. on Power Systems, 13(2) Maio, pp. 576-583.
- Gubner, J. A. (2006) Probability and Random Processes for Electrical and Computer Engineers. Cambridge University Press.
- Howard, R.A. (1971) Dynamic Probabilistic Systems -Volume I: Markov Models. John Wiley & Sons.
- Miller, S. L. & Childers, D. (2004) Probability and Random Processes: with Applications to Signal Processing and Communications. 2nd ed. Academic Press.
- Pereira, F. L. P. (2009), Determinação do Intervalo de Manutenção Programada da Proteção de Linhas de Transmissão Considerando-se Penalidades Associadas à Indisponibilidade. Tese de Doutorado, Pontifícia Universidade Católica do Rio de Janeiro.
- Siqueira, I. P. (1999) Processos de Decisão Markovianos em Sistemas de Segurança e Proteção. Dissertação de Mestrado, Universidade Federal de Pernambuco.
Apêndice
Datas de Publicação
-
Publicação nesta coleção
17 Jan 2012 -
Data do Fascículo
Out 2011
Histórico
-
Revisado
21 Out 2010 -
Recebido
14 Maio 2010