Ferroressonância em Transformadores Alimentados por Cabos Longos

¿Por qué energizar un transformador a través de un cable subterráneo largo genera sobretensiones destructivas de hasta 4.5 p.u.? Desliza este dossier técnico pa

Ing. Francisco Ramírez

Fenomenologia da Ressonância Harmônica e Ferrorressonância em Sistemas Cabo-Transformador

No projeto e na operação de sistemas de distribuição e subtransmissão em média e alta tensão (MT/AT), a substituição de linhas aéreas por cabos subterrâneos de isolação extrudada — principalmente polietileno reticulado (XLPE) — introduz modificações substanciais nos parâmetros distribuídos da rede. A reatância capacitiva transversal de um cabo subterrâneo é ordens de magnitude menor do que a de uma linha aérea equivalente, devido à elevada permissividade relativa do dielétrico (εr2.33.0\varepsilon_r \approx 2.3 - 3.0) e à reduzida distância entre condutor e blindagem. Quando um trecho extenso de cabo subterrâneo alimenta um transformador de potência operando em vazio ou com carga muito reduzida, estabelece-se um circuito oscilante de alta indutância e capacitância capaz de exibir fenômenos de ressonância linear e ferrorressonância.

A ressonância linear se manifesta quando a reatância indutiva equivalente da rede (XLX_L) se iguala exatamente à reatância capacitiva do cabo (XCX_C) a uma determinada frequência, seja ela a fundamental ou uma frequência harmônica específica originada por cargas não lineares. No entanto, o fenômeno mais destrutivo e complexo nessas topologias é a ferrorressonância: um fenômeno dinâmico não linear superarmônico ou subarmônico que surge da interação entre a capacitância distribuída do cabo subterrâneo (CcC_c) e a indutância magnetizante não linear e saturável (Lm(i)L_m(i)) do núcleo ferromagnético do transformador.

Diferente da ressonância linear RLC pura, a ferrorressonância não apresenta uma única frequência de ressonância definida pela equação clássica f0=12πLCf_0 = \frac{1}{2\pi \sqrt{LC}}. Sendo o núcleo do transformador um elemento altamente não linear sujeito a histerese e saturação magnética, a indutância magnetizante variará drasticamente em função do fluxo instantâneo ϕ(t)\phi(t) ou da corrente de magnetização im(t)i_m(t):

Lm(i)=dλdi=NdΦdiL_m(i) = \frac{d\lambda}{di} = \frac{N \cdot d\Phi}{di}

Onde λ\lambda é o fluxo concatenado total, NN é o número de espiras do enrolamento e Φ\Phi é o fluxo magnético no núcleo. Em estado não saturado, a indutância magnetizante é extremamente elevada (Lm10100HL_m \approx 10 - 100 H), mas ao entrar profundamente no regime de saturação, a permeabilidade relativa do material ferromagnético (μr\mu_r) cai drasticamente para valores próximos à permeabilidade do váCuO (μ0\mu_0), reduzindo la indutância incremental a valores equivalentes à indutância de dispersão do enrolamento (LkmHL_k \approx mH).

Esta variação abrupta da indutância permite que o circuito satisfaça a condição de ressonância capacitiva-indutiva para múltiplos valores de tensão e frequência, dando origem a saltos de estado (jump resonance), bistabilidade, sobretensões harmônicas sustentadas com pico distorcido e transientes térmico-dielétricos severos que ameaçam a integridade operacional dos ativos elétricos.

Mecanismos de Excitação e Configurações Críticas

A ferrorressonância em cabos subterrâneos alimentando transformadores é desencadeada comumente por meio de manobras não simultâneas de equipamentos de manobra, ruptura monofásica de condutores ou queima assimétrica de fusíveis de média tensão. As três principais configurações de excitação são detalhadas a seguir:

  • Operação Monofásica de Equipamentos de Manobra (Abertura ou Fechamento Assíncrono): Durante o fechamento ou abertura polo a polo de seccionadores ou disjuntores monofásicos, ou em caso de falha de um polo em um disjuntor tripolar, uma ou duas fases del cabo ficam desenergizadas do lado da fonte, mas permanecem acopladas capacitivamente às fases energizadas através das capacitâncias entre fases (CmCm) e fase-terra (CgCg) do cabo subterrâneo. Esse acoplamento injeta uma corrente capacitiva no enrolamento energizado do transformador, deslocando o ponto de operação para a zona de saturação.
  • Condutor Rompido com Cabo em Vazio: Se um condutor de fase se romper no trecho de alimentação por cabo antes do transformador, a seção isolada do cabo permanece alimentada galvanicamente através da impedância do enrolamento do transformador e da capacitância das outras fases, formando um circuito de retorno ressonante com o neutro flutuante ou aterrado.
  • Acoplamento Capacitivo entre Circuitos Paralelos: Cabos subterrâneos de múltiplos circuitos dispostos no mesmo banco de dutos ou canaleta podem transferir energia capacitiva significativa a um transformador fora de serviço, mas não aterrado, induzindo ferrorressonância por acoplamento capacitivo parasita.

Formulação Matemática Avançada e Bifurcação de Estados

O comportamento dinâmico da ferrorressonância é modelado analiticamente por meio de equações diferenciais ordinárias não lineares. Simplificando a topologia para um circuito equivalente monofásico alimentado por uma fonte senoidal v(t)=Vmcos(ωt+θ)v(t) = V_m \cos(\omega t + \theta) através de uma capacitância de cabo CcC_c e com resistência de perdas série/paralelo RR, a relação entre o fluxo magnético concatenado ϕ(t)\phi(t) e a corrente magnetizante i(t)i(t) do transformador é aproximada por um polinômio ímpar de grau nn (onde habitualmente n=7,9o11n = 7, 9 o 11):

i(ϕ)=aϕ(t)+bϕ(t)ni(\phi) = a \cdot \phi(t) + b \cdot \phi(t)^n

Aplicando a Segunda Lei de Kirchhoff à malha série formada pela fonte de tensão, pela capacitância do cabo e pela indutância magnetizante saturável do transformador, a Equação Diferencial do Sistema responde à equação do tipo Duffing forzada:

d2ϕdt2+1RsCcdϕdt+1Cc(aϕ+bϕn)=ωVmsin(ωt+θ)\frac{d^2\phi}{dt^2} + \frac{1}{R_s C_c} \frac{d\phi}{dt} + \frac{1}{C_c} \left( a \cdot \phi + b \cdot \phi^n \right) = \omega V_m \sin(\omega t + \theta)

Onde RsR_s representa as perdas equivalentes por efeito Joule e perdas no ferro (corrientes parasitas de Foucault e perdas por histerese no núcleo). A não linearidade introduzida pelo termo bϕnb \cdot \phi^n invalida a aplicação do princípio da superposição e gera um mapa de fase multifurcado.

Análise no Espaço de Estados e Diagramas de Bifurcação

Para uma análise trifásica rigorosa, o sistema é formulado no espaço de estados vetorial:

x˙(t)=f(x(t),u(t),μ)\mathbf{\dot{x}}(t) = \mathbf{f}(\mathbf{x}(t), \mathbf{u}(t), \mathbf{\mu})

Onde o vetor de estado x(t)=[ϕa,ϕb,ϕc,vca,vcb,vcc]T\mathbf{x}(t) = [\phi_a, \phi_b, \phi_c, vca, vcb, vcc]^T contém os fluxos nos núcleos de cada fase e as tensões nos capacitores do cabo; u(t)\mathbf{u}(t) é a excitação externa de tensão; e μ\mathbf{\mu} representa o vetor de parâmetros físicos (Lc,Cc,Rcore,longituddecableL_c, C_c, Rcore, longitud de cable).

À medida que o parâmetro de controle (por exemplo, o comprimento do cabo ou a tensão do sistema VrmsVrms) é incrementado continuamente, as soluções estacionárias do sistema sofrem bifurcações topológicas. Os tipos de ferrorressonância resultantes, categorizados segundo a teoria de sistemas dinâmicos, são:

  • Modo Fundamental (Síncrono 1/1): A resposta em corrente e tensão é periódica com o mesmo período T=2π/ωT = 2\pi/\omega que a rede elétrica. Caracteriza-se por correntes de pico muito elevadas e hipersaturação repetida em cada semiciclo.
  • Modo Subarmônico (1/n): A resposta é periódica com um período que é um múltiplo inteiro do período da fonte (Ts=nTT_s = n \cdot T, tipicamente n=3,5n=3, 5). O fluxo magnético oscila a frequências de 16.6Hz16.6 Hz ou 20Hz20 Hz (para redes de 50/60 Hz), gerando um zumbido característico de baixíssima frequência e grave saturação térmica.
  • Modo Quase-Periódico (Quase-periódico): As frequências da resposta não são comensuráveis com a frequência da rede. A trajetória no espaço de fases forma um toroide e não se repete no tempo.
  • Modo Caótico: A trajetória no espaço de estados é determinista, mas extremamente sensível às condições iniciais (Efeito Borboleta). O espectro de frequência é contínuo e plano (ruído de banda larga), e as sobretensões exibem picos imprevisíveis e estocásticos de magnitude devastadora.
ExponentedeLyapunov:λL=limt1tln(δx(t)δx(0))>0ComportamientoCaoˊticoExponente de Lyapunov: \lambda_L = \lim_{t \to \infty} \frac{1}{t} \ln \left( \frac{\|\delta \mathbf{x}(t)\|}{\|\delta \mathbf{x}(0)\|} \right) > 0 \quad \Rightarrow \quad Comportamiento Caótico

Análise Forense de Modos de Falha Dielétrica e Térmica

A ocorrência de eventos ferrorressonantes não detectados ou de longa duração em subsistemas cabo-transformador resulta indefectivelmente em falhas catastróficas. A análise forense de engenharia elétrica classifica as destruições em quatro mecanismos principais de degradação e colapso:

Colapso Térmico por Sobrefluxo Magnético (V/fV/f)

Durante a ferrorressonância em modo fundamental ou subarmônico, o fluxo magnético no núcleo ultrapassa amplamente o fluxo nominal de saturação (Bs1.72.0TeslaB_s \approx 1.7 - 2.0 Tesla), atingindo valores de excursão de até 2.3Tesla2.3 Tesla. Este sobrefluxo extremo satisfaz a condição:

(Vf)real>1.251.50p.u.\left( \frac{V}{f} \right)_{ real } > 1.25 - 1.50 p.u.

Dado que o núcleo ferromagnético está completamente saturado, as linhas de fluxo magnético não podem ser confinadas dentro das chapas de aço-silício de grão orientado (GO) e se dispersam para os meios circundantes de menor relutância: o óleo dielétrico, as ferragens mecânicas de fixação (culatras e tirantes), as blindagens magnéticas e as paredes de chapa de aço do tanque do transformador. O fluxo de dispersão induz correntes parasitas de Foucault massivas na estrutura metálica externa e nos parafusos de fixação. A densidade de potência dissipada por unidade de volume escala-se segundo:

PFoucault=Kf2Bpeak2tmdP_{ Foucault } = K \cdot f^2 \cdot B_{ peak }^2 \cdot t^d_m

Onde tdt_d representa a espessura do material condutor e KK uma constante estrutural. Este fenômeno produz pontos quentes localizados (hot-spots) com temperaturas que superam rapidamente de 300C300^\circ C a 500C500^\circ C em minutos, pirolisando o óleo isolante (geração acelerada de gases dissolvidos como etileno C2H4C_2H_4, hidrogênio H2H_2 e metano CH4CH_4) e degradando irreversivelmente a isolação sólida de celulose (papel Kraft) por meio de despolimerização térmica brusca.

Perfuração Dielétrica por Sobretensões Deslocadas e Harmônicas

As sobretensões de origem ferrorressonante diferem das sobretensões de impulso atmosférico (BIL) ou de manobra (SIL). Elas têm um caráter sustentado no tempo (desde segundos até horas) com formas de onda distorcidas que apresentam elevadas taxas de variação de tensão com relação ao tempo (dV/dtdV/dt). Em transformadores com neutro isolado ou aterrado por meio de alta impedância, o deslocamento do ponto neutro durante uma manobra monofásica eleva a tensão das fases sadias a valores fase-fase (3p.u.\sqrt{3} p.u.) ou superiores (até 3.54.5p.u.3.5 - 4.5 p.u.).

A tensão aplicada sobre o dielétrico do cabo XLPE e sobre a isolação entre camadas/bobinas do enrolamento do transformador supera o gradiente de rigidez dielétrica crítica (EcritEcrit). Isso desencadeia:

  • Árvores de Água (Water Treeing) e Árvores Elétricas (Electrical Treeing) em XLPE: As sobretensões harmônicas aceleram o envelhecimento por descargas parciais (DP) em cavidades microscópicas do dielétrico de polietileno, erosionando o material por meio de descargas dendríticas contínuas até atingir a perfuração total para a terra.
  • Falha da Isolação Entre Espiras em Transformadores: As elevadas harmônicas de alta frequência impõem uma distribuição não uniforme de tensão ao longo do enrolamento primário. O primeiro percentual de espiras próximo ao terminal de entrada suporta até 70%70\% do pico da sobretensión imposta, provocando a perfuração do verniz ou papel do condutor e um curto-circuito entre espiras adjacentes.

Efeitos Térmico-Destrutivos em Para-raios de Óxido Metálico (MOV)

Os para-raios de óxido de zinco (MOV - Metal Oxide Varistors) instalados na transição cabo-transformador são dimensionados para absorver transientes de energia impulsiva de curta duração (microssegundos a milissegundos). Durante um evento de ferrorressonância sustentada, a tensão nos terminais excede continuamente a Tensão de Operação Contínua (MCOV - Maximum Continuous Operating Voltage).

O para-raios entra em zona de condução ativa contínua de alta corrente. A energia dissipada dentro das pastilhas de ZnO é calculada segundo:

Eabsorbed=0tferrovmov(t)imov(t)dtE_{ absorbed } = \int0^{t_{ ferro }} v_{ mov }(t) \cdot i_{ mov }(t) \, dt

Como a duração tferrot_{ ferro } é contínua, la energia acumulada ultrapassa a Capacidade Nominal de Absorção Térmica de Energia do varistor (usualmente expressa em kJ/kVkJ/kV de MCOV, típica de 4a12kJ/kV4 a 12 kJ/kV). Ao superar o limite térmico, as pastilhas de ZnO sofrem uma disparada térmica (thermal runaway), reduzindo sua resistência intrínseca a zero e provocando um curto-circuito fase-terra violento com explosão mecânica e projeção de fragmentos de porcelana ou polímero.

Normatização e Marco de Critérios Internacionais (IEEE / IEC)

O projeto de subsistemas integrados por cabos e transformadores requer um cumprimento estrito das prescrições normativas para garantir a imunidade contra fenômenos ferrorressonantes e a invariância operacional.

As principais normas internacionais aplicáveis são:

  • IEEE C57.105: IEEE Guide for Application of Transformer Connections in Three-Phase Distribution Systems. Fornece diretrizes detalhadas sobre as conexões de transformadores propensas à ferrorressonância e define limites práticos para o comprimento crítico de cabos subterrâneos.
  • IEC 60071-1 / IEEE Std 1313.1: Insulation Coordination - Part 1: Definitions, principles and rules. Define os níveis de isolamento padronizados, os fatores de sobretensão temporária (TOV) e os critérios de avaliação de sobretensões sustentadas não lineares.
  • IEC 60076-1 / IEEE C57.12.00: Power Transformers - Part 1: General. Especifica as tolerâncias de sobrefluxo magnético contínuo (V/f1.10p.u.V/f \le 1.10 p.u. a plena carga, V/f1.05p.u.V/f \le 1.05 p.u. em vazio).
  • IEEE 519: IEEE Standard for Harmonic Control in Electric Power Systems. Estabelece os limites admissíveis de distorção harmônica total de tensão (THDv5%THD_v \le 5\% em MT, 1.52.5%\le 1.5 - 2.5\% em AT) e injeção de corrente harmônica.
Parâmetro Operational / Dielétrico Norma de Referência Limite de Critério Crítico Mecanismo de Falha Incipiente Consequência Térmica/Dielétrica no Ativo
Fator de Sobrefluxo (V/fV/f) IEC 60076-1 / IEEE C57.12.00 >1.10p.u.> 1.10 p.u. (Contínuo)
>1.25p.u.> 1.25 p.u. (t>10st > 10 s)
Dispersão do fluxo magnético para partes estruturais não laminadas. Pirólise do óleo, ponto quente nas culatras (>250C>250^\circ C), degeneração do papel isolante.
Sobretensão Temporária (TOV) IEC 60071-1 >1.50p.u.> 1.50 p.u. (t>1st > 1 s)
>2.50p.u.> 2.50 p.u. (t>0.1st > 0.1 s)
Superação da rigidez dielétrica à frequência industrial. Descargas parciais em XLPE, perfuração entre camadas/entre espiras nos enrolamentos.
Absorção Energética em MOV IEC 60099-4 / IEEE C62.11 Eabs>EthE_{ abs } > E_{ th }
(410kJ/kV4 - 10 kJ/kV)
Excesso de dissipação de potência em regime de condução contínua MCOV. Disparada térmica (Thermal runaway), curto-circuito interno e explosão do para-raios.
Comprimento Relativo do Cabo (LcaboL_{ cabo }) IEEE C57.105 Capacitância do cabo equivalente supera a capacitância crítica CcritC_{ crit }. Formação da malha ressonante série entre a capacitância do cabo e a indutância saturada. Bifurcação caótica ou modo subarmônico autossustentado ao operar equipamentos de manobra monofásicos.
Distorção Harmônica de Tensão (THDvTHD_v) IEEE 519 / IEC 61000-3-6 THDv>5.0%THD_v > 5.0\% (MT)
THDv>2.5%THD_v > 2.5\% (AT)
Injeção harmônica coincidente com frequências naturais de ressonância paralela. Perdas adicionais por correntes de Foucault no cobre e ressonância paralela severa.

Estratégias de Projeto, Mitigação e Topologia de Rede

A prevenção e supressão de fenômenos ferrorressonantes em sistemas alimentados por cabos subterrâneos longos exige uma combinação de topologias de rede robustas, critérios de manobra estritos e elementos de amortecimento passivo/ativo.

Eliminação da Operação Monofásica (Sistemas de Interrupção Tripolar)

A medida preventiva mais efetiva em nível de projeto primário é proibir a instalação de fusíveis monofásicos ou seccionadores de manobra polo a polo na cabeceira da linha subterrânea de alimentação de transformadores de distribuição ou subtransmissão. Em seu lugar, requer-se o uso exclusivo de disjuntores ou seccionadores sob carga de acionamento mecânico tripolar simultâneo. Ao garantir que as três fases se conectem e desconectem com uma dispersão temporal entre polos inferior a 2ms2 ms, elimina-se a injeção de tensão capacitiva através de fases abertas.

Seleção da Topologia e Construção do Núcleo do Transformador

A construção física do núcleo magnético influencia drasticamente a suscetibilidade à ferrorressonância:

  • Núcleos de 3 Colunas (3-Leg Core): Em um núcleo trifásico de três colunas coplanares, o fluxo de sequência zero não dispõe de um caminho de material ferromagnético fechado para retornar. O fluxo de sequência zero deve sair do núcleo e retornar através do óleo e do tanque metálico, o que implica uma indutância de sequência zero (L0L_0) extremamente baixa (tipicamente 0.10.2p.u.0.1 - 0.2 p.u.). Essa baixa indutância fornece um amortecimento natural drástico contra modos ferrorressonantes que envolvam o neutro, tornando esses transformadores extraordinariamente resistentes ao fenômeno.
  • Núcleos de 5 Colunas ou Banco de Transformadores Monofásicos: Em transformadores com núcleo de 5 colunas (ou três unidades monofásicas independentes), existe um caminho de baixíssima relutância para o fluxo de sequência zero através das colunas laterais de retorno. Consequentemente, L0LdL_0 \approx L_d (alta indutância magnetizante de sequência zero), aumentando exponencialmente o risco de ferrorressonância diante de aberturas monofásicas. Deve-se evitar essa topologia se a alimentação for realizada por meio de cabo subterrâneo longo sem disjuntor tripolar.

Resistências de Amortecimento (Damping Resistors) em Terciários em Delta Aberto

Quando o transformador possui uma conexão secundária ou terciária em delta aberto (Open Delta / Broken Delta), pode-se instalar uma resistência de amortecimento (RdR_d) conectada entre os terminais do enrolamento em delta aberto. Em condições normais trifásicas equilibradas, a tensão vetorial somada através do delta é zero (VΔ=0VV_{\Delta} = 0 V), de modo que a resistência não dissipa potência em regime permanente.

Ao ocorrer um desequilíbrio ferrorressonante ou deslocamento de neutro, surge uma tensão de sequência zero (3V03V_0) nos terminais do delta aberto. A resistência dissipa a energia oscilante capacitiva introduzindo um fator de amortecimento (ζ>1/2\zeta > 1/\sqrt{2}) que colapsa as oscilações ressonantes. O valor crítico da resistência de amortecimento referida ao secundário é calculado mediante a relação:

Rd13ωCc(NtercNprim)2R_d \le \frac{1}{3 \cdot \omega \cdot C_c' \cdot \left(\frac{Nterc}{Nprim}\right)^2}

Onde CcC_c' é la capacitância equivalente do cabo vista do primario, e Nterc/NprimNterc/Nprim é a relação de transformação dos enrolamentos.

Manobra Controlada (Point-on-Wave Switching) e Resistências de Pré-inserção (PIR)

Em redes de transmissão e subtransmissão, o emprego de disjuntores dotados de relés de manobra sincronizada ou de resistências de pré-inserção (PIR) permite mitigar o pico de sobretensão e a corrente de magnetização (inrush current) ao energizar cabos subterrâneos e transformadores. Ao energizar o disjuntor no instante exato do valor de pico da onda de tensão para cargas capacitivas, ou na passagem por zero da tensão para a indutância magnetizante do transformador, minimiza-se a componente contínua de fluxo (Φdc=0\Phi_{dc} = 0), prevenindo a entrada do núcleo na região de saturação profunda.

Metodologia de Análise e Aplicação Práctica com o Vexten Suite

Para ilustrar o fluxo de trabalho de análise rigorosa de ferrorressonância e sobretensões em um projeto industrial de alta engenharia, utiliza-se a plataforma integrada Vexten Suite, executando uma sequência metodológica coordenada por meio dos módulos de cálculo de curto-circuito (IEC 60909), dimensionamento de cabos (IEC 60287) e simulação de transientes e resposta em frequência.

Fluxo de Trabalho Metodológico no Vexten Suite

  1. Passo: Caracterização da Rede e Potência de Curto-Circuito (Módulo Vexten Short-Circuit / IEC 60909):

    Modela-se a rede equivalente de alta tensão calculando a potência de curto-circuito inicial (SscSsc'') e a impedância complexa equivalente do sistema (Zk=Rk+jXkZ_k = R_k + j X_k) no ponto de acoplamento comum (PCC).

  2. Passo: Modelagem do Cabo Subterrâneo em Alta Frequência (Módulo Vexten Cable Manager / IEC 60287):

    Calculam-se os parâmetros elétricos distribuídos do circuito trifásico de cabos subterrâneos de XLPE à temperatura máxima de serviço (90C90^\circ C), obtendo a resistência CA (RcRc), a indutância distribuída (LcLc) e a capacitância de sequência positiva e zero (C1,C0C1, C0).

  3. Passo: Cálculo da Impedância de Saída e Varredura de Frequência (Vexten Frequency Scan):

    Calcula-se a impedância de entrada vista a partir dos terminais do transformador em função da ordem harmônica (hh), determinando os pontos de ressonância paralela (antirressonância) e ressonância série:

    Zin(h)=Rc+jhω1Lc+jhω1Lm(i)1h2ω12Lm(i)Cc+jhω1RcCcZ_{ in }(h) = \frac{R_c + j h \omega_1 L_c + j h \omega_1 L_m(i)}{1 - h^2 \omega_1^2 L_m(i) C_c + j h \omega_1 R_c C_c}
  4. Passo: Determinação da Capacitância Crítica e Avaliação do Risco de Ferrorressonância (Vexten Transient Engine):

    Compara-se a capacitância total instalada do cabo (Ctotal=CcLongitudC_{ total } = C_c \cdot Longitud) com o valor limite de capacitância crítica (CcritC_{ crit }) definido segundo a relação empírico-analítica da IEEE C57.105:

    Ccrit=Sn2πf1VLL2(%Im100)KcoreC_{ crit } = \frac{S_n}{2\pi f_1 \cdot V_{L-L}^2} \cdot \left( \frac{\%I_m}{100} \right) \cdot K_{ core }

    Onde SnS_n é a potência nominal do transformador, VLLV_{L-L} é a tensão nominal de linha, %Im\%I_m é a corrente magnetizante em porcentagem, e KcoreK_{ core } é um fator topológico do núcleo (Kcore0.01K_{ core } \approx 0.01 para 3 colunas, Kcore1.0K_{ core } \approx 1.0 para 5 colunas ou monofásicos).

Caso de Estudo Numérico Desenvolvido

Considere um sistema de distribuição industrial com os seguintes parâmetros de campo:

  • Tensão nominal do sistema: VLL=33kVV_{L-L} = 33 kV (50Hz50 Hz).
  • Potência de curto-circuito no nó da rede: Ssc=500MVASsc'' = 500 MVA (X/R=10X/R = 10).
  • Cabo subterrâneo de alimentação: Circuito trifásico unipolar XLPE 1×240mm21 \times 240 mm ^2 Cobre, de comprimento L=6.5kmL = 6.5 km.
    • Capacitância do cabo por unidade de comprimento: Cc=0.26 μF/kmC_c = 0.26 \ \mu F/km.
    • Capacitância total do enlace: Ctotal=6.5×0.26=1.69 μFC_{ total } = 6.5 \times 0.26 = 1.69 \ \mu F.
  • Transformador de potência receptor: Sn=5.0MVAS_n = 5.0 MVA, 33/11kV33/11 kV, conexão Dy11.
    • Núcleo tipo 5 colunas (alta indutância de sequência zero).
    • Corrente magnetizante em vazio: %Im=0.35%\%I_m = 0.35\%.
    • Reatância de dispersão: ek=7.5%e_k = 7.5\%.
  • Manobra avaliada: Abertura monofásica acidental de um fusível ou polo na cabeceira do cabo subterrâneo.

Resolução Analítica e Simulação no Vexten Suite

Cálculo da Capacitância Total do Cabo:

Ctotal=1.69×106FC_{ total } = 1.69 \times 10^{-6} F

A reatância capacitiva à frequência fundamental (50Hz50 Hz) resulta em:

XC=12π501.69×106=15.309×104=1883.5 ΩX_C = \frac{1}{2\pi \cdot 50 \cdot 1.69 \times 10^{-6}} = \frac{1}{5.309 \times 10^{-4}} = 1883.5 \ \Omega

Cálculo da Indutância Magnetizante Não Saturada do Transformador:

A corrente magnetizante nominal por fase (referida ao primário de 33kV33 kV) é:

In=Sn3VLL=5×106333×103=87.48AI_n = \frac{S_n}{\sqrt{3} \cdot V_{L-L}} = \frac{5 \times 10^6}{\sqrt{3} \cdot 33 \times 10^3} = 87.48 A
Im=In(0.35100)=87.48×0.0035=0.3062AI_m = I_n \cdot \left(\frac{0.35}{100}\right) = 87.48 \times 0.0035 = 0.3062 A

A reatância magnetizante equivalente não saturada por fase (XmX_m) é aproximada por:

Xm=VLNIm=33000/30.3062=19052.550.3062=62222.5 ΩX_m = \frac{V_{L-N}}{I_m} = \frac{33000 / \sqrt{3}}{0.3062} = \frac{19052.55}{0.3062} = 62222.5 \ \Omega
Lm=Xm2π50=62222.5314.16198.06HL_m = \frac{X_m}{2\pi \cdot 50} = \frac{62222.5}{314.16} \approx 198.06 H

Determinação da Capacitância Crítica (CcritC_{ crit }) para Núcleo de 5 Colunas:

Aplicando o critério da norma IEEE C57.105 para uma rede de 33kV33 kV com núcleo de 5 colunas (Kcore=1.0K_{ core } = 1.0):

Ccrit=5.0×1062π50(33000)2(0.35100)1.0=5.0×1063.421×10110.0035=5.115×108F=0.05115 μFC_{ crit } = \frac{5.0 \times 10^6}{2\pi \cdot 50 \cdot (33000)^2} \cdot \left( \frac{0.35}{100} \right) \cdot 1.0 = \frac{5.0 \times 10^6}{3.421 \times 10^{11}} \cdot 0.0035 = 5.115 \times 10^{-8} F = 0.05115 \ \mu F

Avaliação Forense do Risco:

Ctotal (1.69 μF)Ccrit (0.05115 μF)C_{ total } \ (1.69 \ \mu F ) \gg C_{ crit } \ (0.05115 \ \mu F )

A capacitância instalada do cabo subterrâneo supera em mais de 33 vezes o limite crítico de ferrorressonância. A relação entre a reatância capacitiva do cabo e a reatância magnetizante na zona de saturação parcial (XmXsat200800 ΩX_m \to X_{ sat } \approx 200 - 800 \ \Omega) intercepta perfeitamente a curva de histerese do transformador.

Ao simular no Vexten Suite a abertura monofásica da fase A, a tensão da fase aberta escala para uma resposta ferrorressonante de Modo Fundamental com picos de sobretensão em estado estacionário de:

Vpico,ferrorresonante=3.12p.u.(Vfasetierra=59.4kVpico)V_{ pico, ferrorresonante } = 3.12 p.u. \quad (V_{ fase-tierra } = 59.4 kV _{ pico })

Junto com uma corrente de pico hipersaturada no enrolamento primário de Ipico=14.2AI_{ pico } = 14.2 A (46veces46 veces a corrente magnetizante nominal), com uma distorção harmônica total de tensão THDv=42.8%THD_v = 42.8\% dominada pela 3ª e 5ª harmônicas.

Ações Corretivas e Redesenho Implementado no Vexten Suite

Para mitigar de forma definitiva o risco iminente de falha dielétrica e colapso térmico identificado na análise do Vexten Suite, prescrevem-se e integram-se as duas seguintes soluções técnicas:

  1. Substituição dos Equipamentos de Manobra de Cabeceira: Eliminação completa de seccionadores/fusíveis monofásicos e instalação de um disjuntor tripolar motorizado com disparo simultâneo dos três polos (tempo de dispersão entre polos Δt<1.5ms\Delta t < 1.5 ms).
  2. Instalação de Resistência de Amortecimento no Terciário em Delta Aberto: Adota-se um transformador de potencial/medição auxiliar com enrolamento secundário em delta aberto de 110V110 V. O valor da resistência de amortecimento RdR_d requerida, dimensionada a partir do motor de simulação do Vexten Suite para dissipar a energia de oscilação capacitiva, é:
    Rd=13ωCtotal(11033000/3)2=13314.161.69×106(0.005774)2=15.309×101=18.83 ΩR_d = \frac{1}{3 \cdot \omega \cdot C_{ total } \cdot \left(\frac{110}{33000/\sqrt{3}}\right)^2} = \frac{1}{3 \cdot 314.16 \cdot 1.69 \times 10^{-6} \cdot (0.005774)^2} = \frac{1}{5.309 \times 10^{-1}} = 18.83 \ \Omega

    Seleciona-se comercialmente um banco de resistências de potência de níquel-cromo de 18 Ω18 \ \Omega, capaz de dissipar uma potência contínua de 5kW5 kW durante transientes de desequilíbrio temporário.

Ao executar novamente a análise dinâmica no Vexten Suite com as modificações topológicas aplicadas, a resposta transitória diante de perturbações assimétricas decae exponencialmente para zero em menos de fdamp<40msf_{ damp } < 40 ms (2 ciclos de frequência industrial), retornando a tensão e o fluxo magnético aos parâmetros nominais estáveis, sem incorrer em zonas de saturação ou sobretensões dielétricas perigosas.