Estudo das limitações impostas ao amortecimento de oscilações eletromecânicas pelos zeros da FTMA de controladores suplementares

Study of the limitations imposed by the zeros of the OLTF of supplementary controllers on the damping of eletromechanic oscillations

SISTEMAS ELÉTRICOS DE POTÊNCIA

Estudo das limitações impostas ao amortecimento de oscilações eletromecânicas pelos zeros da FTMA de controladores suplementares

Study of the limitations imposed by the zeros of the OLTF of supplementary controllers on the damping of eletromechanic oscillations

Ricardo Frangiosi de Moura^I; Marcos Amorielle Furini^II; Percival Bueno de Araujo^III

^IInstituto Federal do Tocantins - Campus Palmas - IFTO AE 310 Sul, Avenida LO 05, s/n, 77021-090, Palmas - TO. [e-mail]

^IIInstituto Federal de São Paulo - Campus Votuporanga - IFSP Av. Jerônimo Figueira da Costa, n. 3014, 15503-110, Votuporanga - SP. [e-mail]

IIIDEE - FEIS - UNESP Av. José Carlos Rossi 1370, 15385-000, Ilha Solteira - SP. [e-mail]

RESUMO

O principal objetivo deste trabalho é analisar a influência de diferentes sinais de entrada para controladores suplementares de amortecimento (POD - Power Oscillation Damping) na localização e controle de dispositivos FACTS UPFC (Unified Power Flow Controller), visando o amortecimento de oscilações de baixa frequência no Sistema Elétrico de Potência (SEP). Serão analisados três sinais de entrada, sendo dois sinais localmente disponíveis (corrente elétrica e fluxo de potência ativa na linha de instalação do UPFC) e um sinal remoto (velocidade angular). Para este fim, o Modelo de Sensibilidade de Potência (MSP) é utilizado para a representação do SEP. De maneira usual, uma análise de sensibilidade, através do método dos resíduos, fornece a localização para instalação do conjunto UPFC/POD e o procedimento para ajuste dos parâmetros do controlador para o amortecimento de um polo de interesse. Entretanto, será mostrado neste trabalho que devido à proximidade de algum zero ao polo de interesse, 0 maior resíduo não necessariamente fornecerá a melhor localização para o conjunto UPFC/POD, fato verificado com maior incidência quando da utilização de sinais locais para a entrada dos controladores. Análises comparativas da eficiência dos sinais de entrada para o POD são apresentadas para o sistema New England composto por 10 geradores, 39 barras e 46 linhas de transmissão.

Palavras-chave: estabilidade a pequenas perturbações, estabilizadores suplementares, FACTS, sistemas elétricos de potência, UPFC, zeros.

ABSTRACT

The main objective of this work is to analyze the influence of different input signals for supplementary damping controller (POD - Power Oscillation Damping) on the location and control of FACTS devices UPFC (Unified Power Flow Controller), in order to damp low frequency oscillations in Electric Power System (EPS). Will analyze three input signals, two signals available locally (electric current and power flow in UPFC line installation) and a remote signal (angular velocity). To this purpose, the Power Sensitivity Model (PSM) is used to represent the EPS. In the usual way, a sensitivity analysis, using the residue method, provides the location for installation of the UPFC/POD and the procedure for setting the parameters of the controller to damp a pole of interest. However, it will be shown in this work due to the proximity of a zero to an eigenvalue of interest, the largest residue does not necessarily provide the best location for all UPFC / POD, as observed more frequently in local signals. Comparative analyses of efficiency of input signals to the POD are presented to the New England system with 10 generators, 39 buses and 46 transmission lines.

Keywords: small-signal stability, supplementary stabilizers, FACTS, electric power system, UPFC, zeros.

INTRODUÇÃO

A análise do comportamento das oscilações eletromecânicas (baixa frequência) de um SEP pode ser realizada através do estudo da estabilidade a pequenas perturbações. A utilização de Estabilizadores do Sistema de Potência (ESP) é realizada de maneira eficaz, econômica e confiável há várias décadas (?). Sua localização (malha de controle da tensão do gerador) e filosofia de aplicação evidenciam que sua principal atuação é no amortecimento das oscilações de modo local, ou seja, associadas aos geradores eletricamente próximos. Quando ajustados para amortecer oscilações do tipo interárea (grupo de geradores de uma área, oscilando contra geradores de outra área) seu desempenho não é tão bom, podendo vir a prejudicar modos locais antes estáveis.

Entretanto, com o avanço da tecnologia empregada na eletrônica de potência, dispositivos automáticos de compensação foram desenvolvidos, os quais receberam o nome de FACTS (Flexible Alternating Current Transmission Systems) que podem introduzir amortecimento adicional às oscilações eletromecânicas (?). Para que ocorra a adição de amortecimento de maneira efetiva, é necessária a introdução de um controlador adicional (por exemplo, do tipo POD) à malha de controle dos dispositivos FACTS (Noroozian and Andersson, 1995). Além disso, para que o controlador POD introduza amortecimento aos modos eletromecânicos problemáticos, deve-se instalá-lo em locais onde possa fornecer o maior amortecimento para o modo oscilatório considerado sem que haja prejuízo aos outros modos estáveis. Nesse contexto, a escolha do sinal de entrada para o conjunto FACTS/POD é de suma importância para a eficiência de sua atuação.

A escolha do sinal de entrada deve ser realizada de forma que este tenha grande observabilidade do modo oscilatório a ser amortecido e que seja localmente disponível (Watanabe et al., 1998). Neste trabalho serão analisados três sinais: corrente elétrica, fluxo de potência ativa e velocidade angular.

A melhor localização para instalação de dispositivos FACTS pode ser obtida da análise dos resíduos da função de transferência em malha aberta (FTMA) do FACTS a ser inserido no sistema (Martins and Lima, 1990). O resíduo obtido para cada localização possível do controlador fornece informações sobre a controlabilidade e observabilidade do conjunto entrada-saída para o controlador em um determinado autovalor. Portanto, o maior resíduo fornecerá maior amortecimento ao modo de oscilação selecionado (Chang and Xu, 2007). Entretanto, a utilização do resíduo como único índice para localização de controladores em sistemas elétricos pode não ser eficaz.

Desta forma, destaca-se que além da análise dos polos da FTMA do controlador é necessária a análise dos zeros desta mesma FTMA, os quais estão intimamente relacionados ao sinal de entrada utilizado e consequentemente, à possibilidade de inclusão de amortecimento (Jones and Andersson, 1998; Mhaskar and Kulkarni, 2006). Portanto, uma má localização dos zeros da FTMA (i.e., proximidade ao polo de interesse) limita a possibilidade de controle do polo de interesse em malha fechada, mesmo que o resíduo obtido para este controlador seja o de maior magnitude (Martins and Lima, 1990).

O estudo sobre a influência dos zeros na localização e controle do dispositivo FACTS UPFC equipado com controlador POD para diferentes sinais de entrada é o foco deste trabalho, onde o SEP foi representado pelo MSP (Deckmann and da Costa, 1994) e todas as simulações e resultados foram obtidos através de implementações utilizando-se o software MATLAB^®.

2 MODELO DE SENSIBILIDADE DE POTÊNCIA

O MSP é uma ferramenta de análise linear do comportamento dinâmico de sistemas elétricos de potência e possui uma característica importante que é a preservação de todas as barras do sistema elétrico no modelo, não ocorrendo a redução às barras internas dos geradores. Como consequência, há a manutenção da integridade da estrutura do sistema elétrico de potência (Deckmann and da Costa, 1994).

Neste modelo, o comportamento dinâmico do sistema elétrico de potência pode ser descrito em função das variáveis de estado (Δx), algébricas (Δz) e de entradas (Δµ), de acordo com (1) onde Δy é um conjunto de saídas desejadas.

Na equação (1), J1, J2, J3 e J4 são submatrizes que relacionam as variáveis de estado (velocidade angular do rotor (Δω) o ângulo interno do rotor (Δδ), a tensão interna de eixo em quadratura (ΔE'_q) e a tensão de campo (ΔE_fd)) com as variáveis algébricas (módulo (ΔV) e o ângulo (Δθ) da tensão terminal das barras do sistema). As entradas do sistema linearizado são a potência mecânica da turbina (ΔP_m), a tensão de referência (ΔV_ref) do regulador automático de tensão de cada gerador e as cargas ativa (ΔP_L) e reativa (ΔQ_L), relacionadas através das submatrizes B1 e B2. As submatrizes C1 e C2 fornecem as saídas desejadas. Nota-se de (1) que o MSP realiza uma representação completa das variáveis envolvidas na modelagem do SEP. Desta forma, a inclusão de novos dispositivos tais como os FACTS ou POD é facilitada.

3 CONTROLADOR FACTS UPFC

O UPFC pode ser considerado um controlador de representação genérica das funções dos controladores FACTS, pois é capaz de atuar de forma simultânea ou seletiva sobre todos os parâmetros que determinam o fluxo de potência em uma linha de transmissão (ângulo, tensão e impedância) (Gyugyi, 1992).

3.1 Inclusão do UPFC no SEP

A inclusão do UPFC no sistema elétrico de potência pode ser feita através de um circuito equivalente composto por uma fonte de tensão (

_s) ideal e uma fonte de corrente (Ī_sh) s ideal, representando, respectivamente, o conversor em série e o conversor em derivação (Meng and So, 2000).

O controle do fluxo de potência na linha de transmissão pode ser feito através do controle da magnitude e do ângulo da tensão série injetada, através dos parâmetros r e y, respectivamente. Assim, a tensão _s pode ser expressa conforme descrito em (2).

A tensão série injetada pelo UPFC pode ser decomposta em duas parcelas, uma em fase (v_p) e outra em quadratura (v_q) com a corrente da linha de transmissão. A atuação do conversor em derivação é representada pelo controle do nível da tensão pela injeção de uma corrente elétrica em quadratura com a tensão da barra inicial (i) de instalação (representada por i, também mostrada na Fig. 2).

Portanto, a inclusão deste dispositivo pode ser feita considerando três variáveis de controle v e i), as quais são adicionadas às variáveis de estado e, através da linearização das equações que as definem, gera o modelo dinâmico linearizado do UPFC conforme mostrado em (3) (Meng and So, 2000).

As constantes de tempo T_vp,T_vq e T_iq representam o atraso inerente do dispositivo de controle.

De acordo com o conjunto de equações (3) nota-se que o UPFC irá operar no modo de compensação reativa, modulando as variações da tensão v_q, através do POD (Moura et al., 2010).

Os controladores ESP e POD utilizados neste trabalho possuem estrutura idêntica e na Fig. 3 é mostrada sua estrutura básica.

Estes controladores são compostos por um ganho estático K (ajustável para fornecer o amortecimento desejado); um bloco washout (definido pela constante de tempo T) e as constantes de tempo T₁, T₂, T₃ e T₄ (responsáveis pela compensação de fase necessária do sinal de entrada escolhido para gerar torque de amortecimento) (Kundur, 1994).

3.2 Localização e Projeto dos Controladores ESP E POD

O correto funcionamento dos controladores é extremamente dependente da localização e do ajuste dos parâmetros para uma determinada condição de operação. Desta forma, o valor dos resíduos das funções de transferência entre a entrada e a saída do controlador (para um determinado modo oscilatório de interesse), é utilizado para se definir a sua localização no sistema elétrico de modo a proporcionar uma contribuição efetiva de amortecimento a modos oscilatórios instáveis ou fracamente amortecidos.

O resíduo Rijk (equação (4)) fornece informações sobre a controlabilidade e observabilidade do conjunto entrada-saída (Cj -Bk) para o controlador em um determinado autovalor (Xi) de acordo com os autovetores direito ( &i) e esquerdo (pi) obtidos da matriz de estados. Portanto, o maior resíduo fornecerá maior amortecimento ao modo de oscilação selecionado e consequentemente dará a melhor localização para o controlador.

O projeto do controlador se resume a calcular as constantes de tempo T₁ = T₃, T₂ = T₄ e o ganho K de forma a introduzir a compensação de fase necessária ao deslocamento do autovalor de interesse. Esse processo é descrito em Aboul-Ela et al. (1996) para controladores ESP, sendo que este mesmo procedimento é aplicável para controladores de dispositivos FACTS (Chang and Xu, 2007). Para isso considere que o ângulo a ser compensado pelo controlador seja β,ω_i seja a frequência em rad/s do modo oscilatório de interesse e que λides é a posição desejada para o autovalor. Com estas definições, o conjunto de equações (5) fornece o procedimento para a obtenção dos parâmetros dos controladorores ESP e POD.

3.3 Influência dos Zeros no Projeto do Controlador

A localização dos zeros da FTMA de um controlador é determinante no comportamento do polo da função de transferência em malha fechada (FTMF) do sistema a ser controlado. Desta forma a localização dos zeros da FTMA de um controlador está intimamente relacionada à facilidade ou dificuldade de controle sobre um sistema (Martins et al., 1992).

Como os polos da FTMF tendem a se deslocar para os zeros da FTMA quando o ganho tende ao infinito, é importante que haja uma boa distância entre os polos críticos e os zeros a fim de se obter um movimento adequado do autovalor de interesse com o aumento do ganho (Martins and Lima, 1990).

Se para uma dada localização do controlador no sistema, o conjunto entrada-saída fornecer boas condições de controlabilidade e observabilidade para o polo de interesse, esta posição será eficiente no amortecimento deste polo de acordo com o método dos resíduos. Observe que neste método focase o projeto no polo de interesse e nenhuma informação é obtida em relação aos zeros da FTMA. Em vista disso, em alguns casos, uma posição satisfatória obtida pelo método dos resíduos pode não fornecer amortecimento adequado devido à má localização dos zeros da FTMA.

Portanto, a localização de zeros em uma determinada região do plano complexo pode limitar a eficácia do controlador projetado pelo método dos resíduos, ou seja, o amortecimento desejado poderá não ser obtido. Assim, a escolha do sinal de entrada do controlador é extremamente importante pois pode fornecer zeros mal localizados.

Neste trabalho, os zeros são obtidos através das funções disponíveis na biblioteca do software MATLAB^®, sendo que o algoritmo utilizado é apresentado de maneira completa em (Emami-Naeini and Dooren, 1982).

3.4 Sinais de Entrada para o Controlador POD

Neste trabalho são analisados três sinais de entrada para o controlador para se verificar qual fornece a melhor atuação com menor valor de ganho e boa separação entre polos e zeros. A corrente elétrica e o fluxo de potência ativa na linha entre as barras i e j (de impedância Z_ij = R_ij + jX_ij) (linha onde está instalado o UPFC) serão considerados (equações (6) e (7), respectivamente). Note que estes sinais estão disponíveis localmente.

Também será utilizada a velocidade angular como sinal de entrada para o POD, sendo que pelo fato desta grandeza ser uma variável de estado, este sinal já está disponível na formulação apresentada na equação (1) e pode ser prontamente utilizado, ao contrário dos sinais anteriores que são dependentes da combinação das variáveis algébricas do modelo.

4 SISTEMA MULTIMÁQUINAS UTILIZADO

O sistema teste utilizado é conhecido como New England e seus dados completos podem ser encontrados em (Araujo and Zanetta Jr., 2001). A partir do diagrama unifilar apresentado na Fig. 4 pode-se afirmar que neste sistema existem duas áreas, sendo que o sistema New York (área 1) é representado de maneira compacta pelo gerador 10.

A partir de um ponto de operação considerado caso base são calculados os autovalores da matriz de estados e os de maior interesse (dominantes) são mostrados na Tabela 1, bem como a frequência de oscilação e o coeficiente de amortecimento (ξ) associados. Dos valores apresentados pode-se afirmar que este sistema possui 4 modos eletromecânicos instáveis, sendo 3 do tipo local e 1 do tipo interárea (conforme se pode concluir a partir das frequências típicas associadas a estes autovalores).

Desta maneira, medidas devem ser tomadas para se tentar estabilizar o sistema de potência. Neste trabalho, a solução para se fornecer amortecimento suplementar ao modo interárea será através da inclusão do conjunto UPFC/POD no sistema. Para os modos locais será a inclusão dos dispositivos ESPs.

5 SIMULAÇÕES E RESULTADOS

5.1 Projeto dos ESPs

A alocação dos controladores ESPs pode ser feita a partirdos resíduos da FTMA. Na Tabela 2 são mostrados os módulos dos resíduos associados às possíveis localizações dos controladores ESPs e aos modos instáveis descritos na Tabela 1, excetuando a máquina geradora 10, que representa a área vizinha.

Os geradores 9 e 5 possuem grande participação na formação de todos os modos locais instáveis, conforme mostrado na Tabela 2, sendo portanto as melhores localizações para instalação dos ESPs.

Na Tabela 3 tem-se os parâmetros dos ESPs instalados nas máquinas geradoras 9 e 5, projetados de acordo com a técnica apresentada na seção 3.2.

Após a inclusão dos ESPs, os 3 modos locais antes instáveis são amortecidos (conforme apresentado na Tabela 4), bem como o modo interárea. Mas, ainda assim, o modo local 1 e o modo interárea são fracamente amortecidos. Para melhorar o desempenho do sistema será incluído o conjunto UPFC/POD para se obter maior amortecimento.

5.2.1 Sinal de Entrada -Velocidade Angular

O amortecimento de oscilações eletromecânicas do tipo interárea pode ser realizado através de sinais remotos, como a velocidade angular. Porém, para que esse amortecimento seja efetivo é necessário que o sinal carregue informações sobre ambas as áreas e, portanto a utilização da diferença entre a velocidade angular de dois geradores (um em cada área) é prática comum (Klein et al., 1991). Para o sistema New England foi considerada a diferença entre a velocidade angular do gerador 10 (área 1) e o gerador 6 (área 2) que possui maior participação no modo interárea, obtido a partir do cálculo dos fatores de participação, conforme a Fig. 5.

5.2 Projeto do Controlador POD

Nesta seção serão apresentados os ajustes do controlador POD para cada um dos sinais de entrada considerados e o objetivo é que o autovalor de interesse (interárea) possua coeficiente de amortecimento igual a 0,1 pu. É importante destacar que o projeto do controlador POD e as análises sobre os efeitos dos zeros no comportamento dinâmico do SEP foram realizadas considerando-se os ESPs em atuação no sistema elétrico de potência.

5.2.1 Sinal de Entrada - Velocidade Angular

O amortecimento de oscilações eletromecânicas do tipo interárea pode ser realizado através de sinais remotos, como a velocidade angular. Porém, para que esse amortecimento seja efetivo é necessário que o sinal carregue informações sobre ambas as áreas e, portanto a utilização da diferença entre a velocidade angular de dois geradores (um em cada área) é prática comum (Klein et al., 1991). Para o sistema New England foi considerada a diferença entre a velocidade angular do gerador 10 (área 1) e o gerador 6 (área 2) que possui maior participação no modo interárea, obtido a partir do cálculo dos fatores de participação, conforme a Fig. 5.

Na Tabela 5 é mostrado o módulo dos resíduos associados às possíveis, e mais significativas localizações do conjunto UPFC/POD no sistema de transmissão.

Dos dados mostrados na Tabela 5 conclui-se que a linha de interligação entre as barras 30 e 10 é a melhor localização para o conjunto UPFC/POD. Em vista disso realiza-se o projeto do controlador POD conforme apresentado na seção 3.2, onde os parâmetros obtidos são mostrados na Tabela 6.

Como esperado a compensação de fase exigida para o sinal de velocidade é pequena, pois o amortecimento de oscilações eletromecânicas é obtido por torques elétricos em fase com esse sinal. Observe que o ganho necessário ao amortecimento aumenta com a diminuição do módulo do resíduo.

A utilização deste sinal remoto é eficiente para todas as localizações analisadas, pois em todos os casos o amortecimento obtido é bem próximo do desejado, de acordo com o apresentado na Tabela 7.

O alto desempenho deste sinal é explicado pela boa separação entre o polo (λ) de interesse (interárea) e os zeros (z) associados à FTMA do controlador conforme a Tabela 8, favorecendo seu deslocamento para o semiplano esquerdo do plano complexo.

Na Fig. 6 é mostrado o comportamento dos polos e do zero para o UPFC/POD instalado na linha 30 - 10 para o ganho do controlador em 0, 0,5K_POD, K,2 K_POD e 3K_POD(K_POD o ganho determinado pelo projeto).

Através do gráfico é possível visualizar que o zero está inicialmente bem distante do polo obtido para ganho nulo, ou seja, o polo da FTMA do controlador. Para os polos obtidos em malha fechada, nota-se que para cada aumento de ganho o polo se desloca ainda mais para o semi-plano esquerdo do plano complexo, inserindo assim maior amortecimento ao sistema. O comportamento para as demais linhas é análogo ao da linha 30 - 10.

Simulações extras permitiram concluir que o amortecimento inserido ao modo interárea fica limitado pela ocorrência de uma bifurcação do tipo sela do modo oscilatório relacionado ao controlador, causando instabilidade do tipo monotônica.

5.2.2 Sinal de Entrada - Fluxo de Potência Ativa

Realizando um procedimento semelhante ao adotado para a velocidade angular, são calculados os módulos dos resíduos associados às possíveis e mais significantes localizações do conjunto UPFC/POD no sistema, considerando o fluxo de potência ativa na linha de instalação do controlador como sinal de entrada. Os resultados obtidos são apresentados na Tabela 9.

Comparando os dados das Tabelas 5 e 9 conclui-se que para ambos os sinais o ranking de classificação das localizações para o conjunto UPFC/POD é o mesmo. Entretanto, o módulo dos resíduos referentes ao sinal fluxo de potência ativa é muito maior quando comparado ao sinal velocidade angular. Esse fato é explicado pela alta observabilidade do modo interárea nesse sinal de entrada, o que ocasiona um menor valor de ganho para os controladores que utilizam o fluxo de potência ativa como entrada, conforme apresentado na Tabela 10.

Outro fato importante é que a compensação de fase exigida com este sinal é muito maior quando comparada à exigida pela velocidade angular. Desta forma, com o conjunto UPFC/POD atuando, foram obtidos os coeficientes de amortecimento nas respectivas linhas do SEP para um coeficiente desejado de 0,1 pu. Observou-se que os melhores coeficientes de amortecimento ocorrem nas linhas de transmissão 30 31 e 37-38, conforme a Tabela 11.

Observa-se dos dados da Tabela 11 que apesar da linha 30-10 possuir o maior resíduo não se obteve o amortecimento desejado, ficando abaixo do coeficiente obtido pela linha 37 38 que foi classificada como a quarta opção pelo método dos resíduos. Este fato é explicado pela existência de um zero mal posicionado da FTMA do controlador associado ao polo de interesse, atraindo-o e dificultando o seu deslocamento.

Na Tabela 12 são apresentadas as distâncias entre o polo (λ) de interesse e o zero (z) da FTMA e conclui-se que as linhas 30 -31 e 37 -38 possuem as duas maiores distâncias (|λ -z|).

Na Fig. 7 é mostrado o comportamento dos polos e do zero na linha de transmissão 30 -10. A distância entre o polo e o zero para o ganho nulo é de praticamente a metade em relação às duas maiores distâncias obtidas para o sistema. Parao ganho de projeto K_PODo polo sofre um pequeno deslocamento no semiplano esquerdo do plano complexo, porém, esse deslocamento não é suficiente para inserir a quantidade de amortecimento desejada ao sistema.

Analisando o sistema apenas pela teoria dos resíduos (Tabela 9), a linha 37 -38 não deveria fornecer o coeficiente de amortecimento mostrado na Tabela 11 já que o seu resíduo é bem menor que o da linha 30 -10. Para se compreender este fato considere a Fig. 8 onde é mostrado o deslocamento do polo em função da variação de K_POD para a linha 37 -38.

Observa-se que a distância do zero ao polo de interesse é bem maior quando comparada com a da linha 30 -10. A maior distância e o fato de o zero estar localizado no semiplano esquerdo do plano complexo é que faz com que sua força de atração traga o polo ainda mais para o semiplano esquerdo, inserindo assim uma quantidade de amortecimento maior ao sistema. Daí, conclui-se que a localização dos zeros influencia diretamente na maior ou menor facilidade que o controlador tem de inserir amortecimento às oscilações do SEP.

5.2.3 Sinal de Entrada -Corrente Elétrica

Na Tabela 13 é apresentado o módulo dos resíduos associados às melhores localizações do UPFC/POD quando da utilização do módulo da corrente elétrica na linha de instalação como sinal de entrada

Comparando-se os dados da Tabela 13 com os da Tabela 9 nota-se que a classificação é alterada para a linha 30-31 mas a magnitude dos resíduos para este sinal local é muito menor quando comparada ao obtido com o fluxo de potência ativa.

De maneira análoga realiza-se o projeto dos controladores através dos resíduos mostrados na Tabela 13 e os respectivos parâmetros são apresentados na Tabela 14.

Os coeficientes de amortecimento obtidos para o modo interárea com a atuação dos controladores cujos parâmetros são apresentados na Tabela 14 para um coeficiente desejado de 0,1 pu são mostrados na Tabela 15, de onde se conclui que a linha 30 -31 forneceu o melhor amortecimento, fato já evidenciado para o sinal fluxo de potência ativa. Entretanto, neste caso o resíduo forneceu a localização corretamente. Observa-se também que a localização nas demais linhas não foi eficiente na obtenção do amortecimento desejado. Isto se deve a uma má separação entre o polo de interesse e o zero associado, conforme mostrado na Tabela 16.

O sistema de energia elétrica é continuamente submetido a pequenas variações de cargas, com os consequentes ajustes na geração, através da variação da abertura do dispositivo de captação de água (no caso de usinas hidrelétricas), ocasionando um aumento ou diminuição da potência mecânica de entrada dos geradores (Kundur, 1994). Em vista disso, para a análise do desempenho dinâmico do sistema elétrico frente a uma pequena perturbação, considerou-se que para um pequeno aumento na carga do sistema, corresponderá um pequeno ajuste na geração. Este ajuste será aqui representado por um degrau de 0,05 pu na potência mecânica de entrada da unidade geradora 1. Além disso, os geradores 9 e 5 estão equipados com ESPs ajustados com os parâmetros mostrados na Tabela 3.

Dos dados da Tabela 16 conclui-se que a separação entre o polo e o zero da FTMA para a linha 30 -31 é praticamente o dobro da separação para a linha 30 -10, fato já observado para o sinal fluxo de potência. Portanto, o descolamento do polo de interesse (modo interárea) é favorecido na linha 30-31.

Na Fig. 9 é apresentado o comportamento dos polos e do zero para a linha 30-31 . Verifica-se claramente através do deslocamento dos polos que, aumentando-se o ganho do controlador o polo de interesse caminha em direção ao zero.

5.3 Efeitos Sobre o Comportamento Dinâmico do SEP

O sistema de energia elétrica é continuamente submetido a pequenas variações de cargas, com os consequentes ajustes na geração, através da variação da abertura do dispositivo de captação de água (no caso de usinas hidrelétricas), ocasionando um aumento ou diminuição da potência mecânica de entrada dos geradores (Kundur, 1994). Em vista disso, para a análise do desempenho dinâmico do sistema elétrico frente a uma pequena perturbação, considerou-se que para um pequeno aumento na carga do sistema, corresponderá um pequeno ajuste na geração. Este ajuste será aqui representado por um degrau de 0,05 pu na potência mecânica de entrada da unidade geradora 1. Além disso, os geradores 9 e 5 estão equipados com ESPs ajustados com os parâmetros mostrados na Tabela 3.

Na Fig. 10 é mostrado o comportamento do ângulo interno da máquina geradora 6 (pertencente à área 2), em relação ao ângulo interno da unidade geradora 10 (pertencente à área 1), Δδ₆ - Δδ₁₀. As diferentes curvas se referem ao conjunto UPFC/POD instalado nas linhas 30 -10 e 37 -38, considerando-se o fluxo de potência ativa como entrada para o controlador onde ele encontra-se instalado, pois se enquadra na filosofia de utilização de sinais locais (Aboul-Ela et al., 1996). As linhas em questão possuem o maior e o menor resíduo, respectivamente, conforme dados apresentados na Tabela 9. Analisando as curvas mostradas conclui-se que os comportamentos são semelhantes e que o sistema para os dois casos obtém um bom amortecimento. Através da análise apenas dos resíduos não seria possível chegar a esta conclusão já que a linha 37 - 38 possui um resíduo praticamente 3 vezes menor que o resíduo da linha 30 - 10. Assim, não basta apenas que para um dado sinal de entrada do controlador, este seja instalado no local de maior resíduo. É desejável que o zero possua uma boa separação do polo obtido a partir da FTMA e que este zero esteja bem localizado no plano complexo.

Nos gráficos da Fig. 11 é mostrado o comportamento dinâmico do ângulo interno do gerador 6 (Δδ₆ - Δδ₁₀) considerando-se o conjunto UPFC/POD instalado nas melhores posições, as quais foram obtidas pelas análises dos resíduos e zeros da FTMA do controlador POD (para os sinais fluxo de potência ativa e corrente elétrica instalação na 30 - 31; para o sinal velocidade angular instalação na 30 - 10). Assim, é possível analisar para qual sinal de entrada do controlador o sistema apresentou um amortecimento global mais significativo.

A corrente elétrica e o fluxo de potência ativa (sinais disponíveis no local de instalação do conjunto UPFC/POD) apresentam curvas muito semelhantes, demonstrando que para estas condições, o amortecimento fornecido foi suficiente para que se obtivesse um bom desempenho do sistema elétrico de potência quando submetido a uma pequena perturbação. No caso do uso da velocidade angular como sinal de entrada do controlador, esperava-se um desempenho melhor do SEP quando comparado com os obtidos para os outros dois sinais de entrada. Isto porque o coeficiente de amortecimento fornecido ao modo interárea com o sinal velocidade angular foi o maior obtido conforme apresentado na Tabela 17. Porém, a curva relativa ao sinal da velocidade angular apresentado na Fig. 11 tem o pior comportamento para os 3 sinais de entrada utilizados. Isto pode ser explicado analisando-se os dados mostrados nas Tabelas 17, 18 e 19, nas quais são apresentados os quatro modos oscilatórios de interesse considerandose a inclusão dos ESPs e do conjunto UPFC/POD para cada um dos sinais de entrada sob análise, ajustados com os parâmetros mostrados nas Tabelas 6, 10 e 14.

Analisando-se os autovalores mostrados na Tabela 17 podese afirmar que o conjunto UPFC/POD forneceu um ótimo coeficiente de amortecimento ao modo interárea mas prejudicou o amortecimento do modo local 1, que passou a ter um coeficiente de amortecimento associado menor que o mostrado na Tabela 4 (atuação apenas dos ESPs). Isto não ocorre para os outros sinais utilizados pois o maior amortecimento inserido ao modo interárea não prejudicou os modos locais, conforme atestam os valores mostrados nas Tabelas 18 e 19.

Assim, conclui-se que analisando-se apenas o amortecimento do modo interárea, a velocidade angular é um ótimo sinal de entrada para o controlador mas para o comportamento global do sistema ele não é eficaz pois dificulta o amortecimento dos modos oscilatórios restantes e consequentemente afeta a estabilidade do sistema elétrico. Desta maneira os sinais de entrada potência ativa e corrente elétrica são escolhas melhores pois além de inserir amortecimento ao modo interárea não deteriora em demasia os demais modos de oscilação.

A dificuldade observada para o amortecimento do modo local 1 devido a aplicação do UPFC/POD com sinal de entrada velocidade angular, é melhor compreendida analisando-se a localização dos zeros e dos polos relacionados a este modo, obtidos através da FTMA do controlador, considerando os 3 sinais de entrada (Fig. 12).

Observa-se que o zero relacionado ao modo local 1 obtido para a FTMA do UPFC/POD utilizando a velocidade angular como sinal de entrada está localizado no semiplano direito do plano complexo. Esta má localização faz com que para este sinal de entrada, o polo relacionado ao modo local 1 será atraído pelo zero e consequentemente ocorrerá a instabilidade na medida em que se aumenta o ganho do controlador. Assim, para o ganho de projeto obtido com a velocidade angular como entrada do controlador, obtém-se um ótimo coeficiente de amortecimento em relação ao modo interárea, porém, o modo local 1 é bastante afetado, prejudicando o comportamento dinâmico do sistema como um todo.

Na Fig. 13 pode ser melhor observada a atração que o zero exerce sobre o polo relacionado ao modo local 1, com o aumento do ganho do controlador que utiliza como sinal de entrada a velocidade angular.

Observe que variando-se o ganho do controlador de 0 a 1000 vezes o valor do ganho de projeto (K_POd), é possível visualizar o deslocamento do polo em direção ao zero. Portanto, para que um determinado modo de oscilação não seja negativamente afetado pelo ganho do controlador, é desejável que todos os zeros das FTMA estejam localizados no semiplano esquerdo do plano complexo ou o mais próximo deste.

6 CONCLUSÕES

O principal objetivo deste trabalho foi analisar a eficiência de três sinais de entrada para controladores suplementares de amortecimento, acoplados ao dispositivo FACTS UPFC para o amortecimento das oscilações eletromecânicas de modo interárea.

Concluiu-se das simulações realizadas e resultados obtidos que índices como os resíduos das funções de transferência fornecem informações extremamente valiosas, tanto para se averiguar as ações a ser tomadas para estabilização do sistema, quanto para se verificar a atuação das ações adotadas em cada modo oscilatório. Os resíduos podem determinar a melhor localização do dispositivo FACTS no sistema elétrico de potência, inserindo alto nível de amortecimento às oscilações eletromecânicas do sistema elétrico de potência, porém não é sensível à presença de zeros.

Constatou-se também que quanto mais à esquerda do semiplano esquerdo o zero estiver localizado, melhor para a inserção de amortecimento, pois este zero vai atrair o polo como foi verificado através dos resultados apresentados. Assim, a introdução do amortecimento desejado depende da localização dos zeros da função de transferência em malha aberta do conjunto UPFC/POD. Quanto maior a distância entre o polo de interesse e o zero relacionado a ele, melhor a eficiência do controlador.

Dos sinais analisados, a velocidade angular (sinal remoto) apresentou resultados satisfatórios quanto ao amortecimento inserido ao modo interárea mas retirou amortecimento do modo local 1, prejudicando o comportamento dinâmico global do sistema. A corrente elétrica e o fluxo de potência ativa são sinais locais, mas ocorrem zeros próximos ao polo de interesse, limitando a inserção de amortecimento e, neste caso, uma localização fornecida pelo método dos resíduos pode não ser eficaz.

Portanto, conclui-se que a localização, o sinal de entrada e o ajuste dos parâmetros de controladores POD no sistema elétrico de potência são tarefas complexas que exigem alto grau de conhecimento do sistema a ser controlado e que o estudo do posicionamento e deslocamento dos polos e zeros da FTMA e da FTMF dos controladores fornecem informações valiosas para o controle do sistema elétrico de potência.

Dever ser ressaltado que os estudos aqui apresentados utilizando-se o dispositivo UPFC também são válidos para controladores ESPs e os demais FACTS existentes.

AGRADECIMENTOS

Agradecemos o suporte financeiro da CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior) no desenvolvimento desta pesquisa.

Artigo submetido em 23/09/2010 (Id.: 1195)

Revisado em 25/11/2010, 14/04/2011, 09/06/2011

Aceito sob recomendação do Editor Associado Prof. Luís Fernando C. Alberto

• Aboul-Ela, M., Sallam, A., McCalley, J. and Fouad, A. (1996). Damping controller design for power system oscillations using global signals, IEEE Transactions on Power Systems 11(2):767-773.
• Araujo, P. and Zanetta Jr., L. (2001). Pole placement method using the system matrix transfer function and sparsity, International Journal of Electrical Power EnergySystems 23(3):173-178.
• Chang, Y. and Xu, Z. (2007). A novel svc supplementary controller based on wide area signals, Electric Power Systems Research 77(12):1569-1574.
• Deckmann, S. and da Costa, V. (1994). A power sensitivity model for electromechanical oscillation studies, IEEE Transactions on Power Systems 9(2):965-971.
• Emami-Naeini, A. and Dooren, P. V. (1982). Computation of zeros of linear multivariable systems, Automatica 18(4):4 1 5-430.
• Gyugyi, L. (1992). Unified power-flow control concept for flexible ac transmission systems, IEE Proceedings Generation, Transmission and Distribution 139(4):323-331.
• Jones, L. E. and Andersson, G. (1998). Application of modal analysis of zeros to power systems control and stability, Electric Power Systems Research 46(3):205-211.
• Klein, M., Rogers, G. and Kundur, P. (1991). A fundamental study of inter-area oscillations in power systems, IEEE Transactions on Power Systems 6(3):914-921.
• Kundur, P. (1994). Power System Stability and Control, New York: McGraw-Hill.
• Martins, N. and Lima, L. (1990). Determination of suitable locations for power system stabilizers and static var compensators for damping electromechanical oscillations in large scale power systems, IEEE Transactions on Power Systems 5(4):1455-1469.
• Martins, N., Pinto, H. and Lima, L. (1992). Efficient methods for finding transfer function zeros of power systems, IEEE Transactions on Power Systems 7(3):1350-1361.
• Meng, Z. and So, P. (2000). A current injection upfc model for enhancing power system dynamic performance, IEEE Power Engineering Society Winter Meeting, 2000., Vol. 2, pp. 1544-1549.
• Mhaskar, U. and Kulkarni, A. (2006). Power oscillation damping using facts devices: modal controllability, observability in local signals, and location of transfer function zeros, IEEE Transactions on Power Systems 21(1):285-294.
• Moura, R. F., Furini, M. A. and Araujo, P. B. (2010). Análise de resíduos e zeros da função de transferência de controladores suplementares de dispositivos facts upfc para diferentes sinais de entrada, XVIII Congresso Brasileiro de Automática, Bonito/MS, pp. 3369-3375.
• Noroozian, N. and Andersson, G. (1995). Damping of Watanabe, E. H., Barbosa, P. G., Almeida, K. C. and Taranto, G. N. (1998). Tecnologia facts- tutorial, SBA Controle & Automação 9(1):39-55.
• Watanabe, E. H., Barbosa, P. G., Almeida, K. C. and Taranto, G. N. (1998). Tecnologia facts- tutorial, SBA Controle & Automação 9(1):39-55.

Datas de Publicação

• Publicação nesta coleção
  28 Maio 2012
• Data do Fascículo
  Abr 2012

Histórico

• Revisado
  09 Jun 2011