Ferrorressonância em Reatores de Derivação por Manobra Monofásica

Introdução e Marco Teórico dos Fenômenos Não Lineares em Sistemas de Alta Tensão

A integração de sistemas de transmissão em corrente alternada de ultra alta tensão (UHVAC) e extra alta tensão (EHVAC) introduz complexidades operativas significativas associadas à gestão da potência reativa capacitiva gerada pelas linhas de transmissão em vazio ou com carga leve. Para mitigar as sobretensões a frequência industrial e controlar o perfil de tensão sazonal, o uso de reatores de derivação (shunt reactors), seja de núcleo de ferro com entreferros distribuídos ou de núcleo de ar blindados, é uma prática padrão indispensável. Não obstante, a interação entre a indutância não linear do núcleo magnético do reator e a capacitância de sequência positiva e zero da linha de transmissão cria um circuito LC ressonante exposto a fenômenos de ferro-ressonância e saturação magnética transitória severa durante manobras operativas desbalanceadas, particularmente a comutação monofásica.

Desde a perspectiva da teoria de circuitos eletromagnéticos avançados, a ferro-ressonância não é uma condição de ressonância linear clássica, mas um modo de oscilação não linear caracterizado pela coexistência de múltiplos estados estáveis (multi-estabilidade) para os mesmos parâmetros de rede. O sistema pode saltar abruptamente de um estado de baixa tensão e corrente quase sinusoidal para um estado ferro-ressonante de alta tensão, caracterizado por sobretensões temportemporary (TOV) sustentadas, distorção harmônica massiva, correntes de excitação extremamente elevadas e solicitações térmicas e dielétricas extremas sobre os enrolamentos do reator e os equipamentos associados da subestação.

O análise formal desse fenômeno requer a formulação matemática rigorosa do fluxo magnético no núcleo do reator acoplado com as equações diferenciais estocásticas e determinísticas do sistema elétrico trifásico com desequilíbrios introduzidos pela abertura ou fechamento sequencial de polos em disjuntores de alta tensão (GIS ou AIS).

Equações Governantes e Modelagem Não Linear do Núcleo Magnético

Para modelar com precisão a saturação magnética e a ferro-ressonância, é fundamental representar a característica de magnetização não linear do núcleo de ferro do reator de derivação. A relação entre o fluxo magnético encadeado \lambda(t) e a corrente de excitação i(t) é expressa mediante funções analíticas avançadas que incorporam os efeitos de histerese e as correntes parasitas (perdas no núcleo). Uma aproximação clássica utilizada em estudos de transitórios eletromagnéticos (EMTP) baseia-se na função de Fröhlich modificada ou em expansões em série de potências ímpares para evitar descontinuidades:

i(t) = a \lambda(t) + b \lambda(t)^{2n+1}

Onde $a$ representa a inversa da indutância diferencial na zona linear ( $L_0^{-1}$ ), $b$ é o coeficiente de saturação do núcleo, e $n$ é um inteiro positivo que determina a agudeza do joelho de saturação. A equação diferencial que rege o circuito série equivalente composto pela resistência da linha $R$ , a indutância parasita ou a do próprio reator com sua característica não linear, e a capacitância efetiva do sistema $C$ vista a partir dos terminais do reator é formulada como segue:

\frac{d^2\lambda}{dt^2} + \frac{R}{L(\lambda)} \frac{d\lambda}{dt} + \frac{1}{C L(\lambda)} \lambda = \frac{1}{C} v_s(t)

Onde v_s(t) = V_m \sin(\omega t + \phi) representa a fonte de tensão sinusoidal do sistema, e a indutância aparente ou diferencial é definida como L(\lambda) = \frac{d\lambda}{di}. Quando o núcleo entra em saturação profunda, L(\lambda) diminui drasticamente para valores próximos à indutância de ar do enrolamento ( $Lair$ ), o que altera de maneira abrupta a frequência de ressonância natural do circuito:

\omega_0(\lambda) = \frac{1}{\sqrt{C L(\lambda)}}

Durante uma manobra de comutação monofásica (por exemplo, a abertura de uma única fase de um disjuntor devido à operação unipolar de religamento ou um defeito monofásico eliminado), as fases sãs continuam acopladas capacitivamente à fase aberta. O reator conectado à fase aberta fica flutuando através de uma rede capacitiva acoplada às outras duas fases energizadas e à terra. Isto estabelece um circuito sintonizado onde a capacitância equivalente $Ceq$ e a indutância não linear L(\lambda) satisfazem as condições para a excitação de modos sub-harmônicos (\omega_0 = \omega/2, \omega/3, \dots), harmônicos (3\omega, 5\omega, etc.) ou quase-periódicos.

Análise de Manobras de Comutação Monofásica e Origem do Desequilíbrio

As manobras de abertura e fechamento unipolar em sistemas EHVAC e UHVAC são operativas ou de proteção por natureza (ex. religamento automático unipolar após falhas monofásicas temporárias). No intervalo de tempo em que uma fase está aberta no disjuntor mas o reator permanece conectado no lado da linha, a topologia da rede muda drasticamente. As capacitâncias mútuas entre fases ( $Cab, Cbc, Cca$ ) e as capacitâncias para a terra ( $Can, Cbn, Ccn$ ) desempenham um papel crítico.

Considerando uma linha de transmissão transposta de comprimento $l$ , a matriz de capacitâncias por unidade de comprimento transforma-se, ao abrir-se a fase $A$ , em uma rede Thévenin equivalente vista a partir dos terminais do reator da fase $A$ . A tensão induzida eletrostática e eletromagneticamente a partir das fases sãs $B$ e $C$ energiza o circuito ressonante. Se a reatância indutiva do reator à frequência industrial se aproxima ou cruza o valor da reatância capacitiva equivalente da rede desequilibrada, o sistema entra em uma zona de instabilidade severa.

Parâmetro de Operação / Rede	Condição Normal (Simétrica)	Condição de Comutação Monofásica	Limite Normativo (IEC / IEEE)
Tensão de Sequência Zero ( $V_0$ )	Praticamente 0% a 1% de $V_1$	Pode alcançar até 30% - 50% de $V_1$	IEEE Std 1036: Máx 5% sob condições assimétricas
Corrente de Excitação do Reator	Corrente nominal de projeto (I_n \le 1.0 p.u. )	Picos de até 8 - 15 p.u. por saturação profunda	IEC 60289: Tolerância térmica de sobrecorrente transitória
Distorção Harmônica Total de Tensão (THDv)	< 1.5\% em barras EHV	15\% - 40\% com predomínio de harmônicos 2°, 3° e 5°	IEEE Std 519: Limite estrito de 1.5% a 2.0% em EHV
Fator de Sobretensão Temporal (TOV)	$1.0 p.u.$	$1.8 p.u. - 2.6 p.u.$ sustentados por vários ciclos	IEC 60071-2: Curvas de suportabilidade temporal a 10s - 30s

Análise Forense de Falhas e Consequências em Equipamentos de Subestação

Quando um reator de derivação experimenta ferro-ressonância sustentada induzida por comutação monofásica, as consequências físicas e dielétricas sobre os ativos da subestação são catastróficas se as proteções não atuarem com a velocidade adequada. A seguir, detalham-se os mecanismos de falha forense por cada tipo de equipamento:

Falhas em Enrolamentos de Reatores de Derivação

As correntes de saturação altamente distorcidas geram forças eletrodinâmicas radiais e axiais massivas sobre os condutores do enrolamento. Estas forças superam os limites elásticos do cobre e dos suportes mecânicos, provocando o deslocamento físico das espiras, a deformação dos canais de refrigeração e a abrasão do isolamento de papel celulósico. Do mesmo modo, as perdas suplementares por correntes parasitas nos enrolamentos e nos elementos estruturais metálicos do tanque devido ao conteúdo harmônico elevado (h \ge 3) produzem pontos quentes locais (hotspots) com temperaturas que superam os 250 °C, degradando rapidamente o óleo isolante e gerando gases combustíveis característicos ( $CO, CO₂, CH_4, C_2H_2$ ).

Deterioração Dielétrica em Cabos de Alta Tensão e Buchas (Bushing)

As sobretensões temporárias (TOV) ferro-ressonantes com componentes de alta frequência submetem o isolamento polimérico (XLPE) ou impregnado em óleo dos cabos de conexão e das buchas a esforços dielétricos repetitivos. Esses esforços aceleram os fenômenos de descarga parcial (PD) interna, erosionando a estrutura molecular do isolamento até provocar a perfuração dielétrica sob a forma de árvore elétrica (electrical treeing) e a conseguinte falha para a terra catastrófica com explosão da bucha.

Comportamento Anômalo em Disjuntores (Circuit Breakers)

Durante a interrupção de correntes com alto conteúdo harmônico e componente de corrente contínua (DC) decrescente geradas pelo regime ferro-ressonante, os disjuntores a SF6 enfrentam severas dificuldades na passagem por zero da corrente. Isto incrementa a probabilidade de re-ignições (reignitions) ou re-arco (restrikes), injetando transitórios de alta frequência adicionais que agravam o colapso dielétrico geral da rede.

Estratégias de Projeto Prático e Mitigação Avançada

A mitigação efetiva da instabilidade por ferro-ressonância e saturação magnética em reatores de derivação exige a implementação de critérios de projeto multidisciplinares tanto na fase de engenharia da rede quanto na especificação dos equipamentos. As fórmulas e fatores de projeto críticos detalham-se a seguir:

Seleção da Característica de Linearidade do Núcleo

Para evitar a saturação prematura ante desequilíbrios monofásicos, o projeto do reator deve garantir que o joelho da curva de saturação magnética se situe pelo menos entre 140% e 160% da tensão nominal de operação contínua:

V_{ knee } \ge 1.45 \cdot V_{n\_phase}

Adicionalmente, o uso de reatores com núcleo de ar blindados elimina por completo a não linearidade introduzida pelos núcleos ferromagnéticos, eliminando o risco de ferro-ressonância. No entanto, seu maior tamanho, peso e campo magnético disperso requerem estudos detalhados de compatibilidade eletromagnética (EMC).

Implementação de Resistores de Amortecimento e Pára-raios

Quando são empregados reatores de núcleo de ferro, a instalação de pára-raios de óxido metálico (MO Surge Arresters) de alta capacidade energética diretamente nos terminais do reator ajuda a limitar as cristas de tensão ferro-ressonante. A energia absorvida pelo pára-raios durante um evento transitório prolongado calcula-se mediante a integral de Joule:

W_{ surge } = \int0^{t_{ clear }} v(t) i(t) dt

Da mesma forma, a incorporação de um esquema de aterramento através de um resistor de amortecimento no neutro do reator (desde que a configuração do sistema o permita) dissipa a energia oscilatória de sequência zero, quebrando a condição de ressonância.

Aplicação Prática e Análise através da Vexten Suite

No ambiente de engenharia avançada da Vexten Academy, a validação dos estudos de ferro-ressonância e saturação é executada utilizando módulos integrados que combinam os padrões internacionais IEC e IEEE. O fluxo de trabalho computacional na plataforma Vexten Suite abrange as seguintes etapas normativas:

Cálculo de Curto-circuito e Níveis de Contribuição segundo IEC 60909 / IEEE 141

O primeiro passo consiste em determinar a potência de curto-circuito da rede no ponto de interconexão do reator para avaliar a relação de rigidez do sistema ( $SCR$ - Short Circuit Ratio):

SCR = \frac{S_{ sc\_3ph }}{Q_{ reactor }}

Se o $SCR$ é baixo ( $SCR < 3$ ), a rede é considerada fraca, o que incrementa exponencialmente a suscetibilidade a sobretensões por ferro-ressonância. O motor de cálculo da Vexten Suite processa matrizes de impedância nodal $Zbus$ para extrair as componentes simétricas exatas durante a abertura unipolar do disjuntor.

Derating Harmônico e Dimensionamento de Cabos segundo IEC 60287 / NEC 310

A presença de correntes harmônicas elevadas induzidas pela saturação magnética do reator altera as perdas por efeito Joule e por correntes parasitas nos cabos de potência associados à subestação. O fator de derating por harmônicos ( $FHD$ ) calcula-se na Vexten Suite aplicando os coeficientes especificados na IEC 60287:

I_{ max\_cable } = I_{ base } \cdot \sqrt{\frac{1}{\sum_{h=1}^{k} I_h^2 \cdot (1 + Y_h)}}

Onde $I_h$ é a magnitude relativa do harmônico $h$ e $Y_h$ representa o fator de perdas suplementares dependente da frequência para os condutores e blindagens metálicas do cabo.

Mitigação de Fator de Potência e Ressonância Harmônica

Para evitar a coincidência entre a frequência de ressonância paralela do sistema de compensação reativa e os harmônicos característicos gerados pela saturação, a Vexten Suite realiza uma varredura de frequência (Frequency Scan) calculando a impedância de entrada Zin(j\omega) na faixa de $0 Hz$ a $3 kHz$ :

Zin(j\omega) = \left( Rsys(\omega) + j\omega Lsys(\omega) \right) \parallel \left( Rreact(\omega) + j\omega L_{\lambda}(\omega) \right) \parallel \left( \frac{1}{j\omega Ceq} \right)

Os resultados desta varredura permitem aos engenheiros de projeto ajustar os parâmetros de sintonia dos filtros de harmônicos ou modificar a indutância nominal dos reatores de derivação para garantir que os polos de impedância fiquem estritamente afastados das frequências harmônicas dominantes (2^\circ, 3^\circ, 5^\circ, 7^\circ), assegurando assim a estabilidade operativa total do sistema de transmissão EHVAC/UHVAC sob qualquer contingência de comutação monofásica.

Conclusões e Recomendações de Engenharia

A instabilidade por ferro-ressonância e a saturação magnética em reatores de derivação durante manobras de comutação monofásica representam um desafio crítico na engenharia de sistemas elétricos de alta potência. A análise rigorosa mediante modelos não lineares de indutância dependente do fluxo magnético, combinada com a aplicação estrita de normativas internacionais (IEC 60289, IEC 60071, IEEE Std 1036 e IEEE Std 519), é indispensável para prevenir falhas catastróficas nos ativos da subestação. A utilização de ferramentas de simulação avançada como a plataforma Vexten Suite garante que as especificações de projeto, os fatores de derating térmico e os estudos de varredura de frequência mitiguem eficazmente o risco de sobretensões temporárias e distorção harmônica extrema, consolidando a confiabilidade e segurança operativa das redes de transmissão do futuro.