Deslocamento de Neutro e TOV em Mineração a Grande Altitude: Fator K conforme NFPA 70B

Contextualização Eletrofísica de Operações Mineradoras em Grandes Altitudes

A operação de sistemas elétricos de potência em instalações mineradoras situadas em altitudes superiores a 3.000 metros acima do nível do mar (m.a.n.m.) impõe restrições severas à rigidez dielétrica do meio isolante e às capacidades de dissipação térmica dos equipamentos. Em altitudes elevadas, a densidade relativa do ar ($\delta$) diminui drasticamente de acordo com a lei barométrica fundamental:

Onde $P$ e $T$ representam a pressão e a temperatura absoluta na altitude $h$ (m.a.n.m.), $P_0$ e $T_0$ são as condições padrão ao nível do mar (101,325 kPa e 288,15 K), $g$ é a aceleração da gravidade, $M$ é a massa molar do ar seco e $R$ é a constante universal dos gases. Esta redução de densidade altera o livre caminho médio dos elétrons livres, incrementando a probabilidade de ionização por colisão sob a aplicação de um campo elétrico, segundo a Lei de Paschen:

Como resultado direto, as distâncias de isolamento no ar (fase-terra e fase-fase), bem como as distâncias de escoamento (creepage distance), experimentam uma degradação não linear.

De acordo com a norma IEC 60071-2 e IEEE C37.100.1, o fator de correção por altitude para o isolamento externo ($K_a$) é definido segundo a formulação matemática:

Onde $m$ é um parâmetro empírico que depende do tipo de tensão aplicada (frequência industrial, impulso de manobra ou impulso atmosférico/raio) e do valor da tensão suportável. Para instalações mineradoras situadas a 4.500 m.a.n.m., o fator $K_a$ tipicamente oscila entre 1,45 e 1,60, o que exige sobredimensionar a tensão suportável nominal à frequência industrial ($V_w$) de disjuntores, isoladores e painéis (switchgears) de média tensão segundo:

Simultaneamente, a dissipação térmica por convecção natural é drasticamente atenuada devido à massa específica reduzida do fluido refrigerante (ar). O fator de derating térmico ($kalt$) aplicado à corrente nominal de equipamentos como transformadores, barramentos condutores e disjuntores é aproximado mediante a norma IEEE C57.96:

Onde $\alpha \approx 0,005$ a $0,008$ segundo a geometria do elemento radiante e a relação entre perdas por radiação e convecção. A interação simultânea da degradação dielétrica e do derating térmico estabelece um cenário crítico para o surgimento de sobretensões temporárias (TOV) e fenômenos de deslocamento de neutro em redes com neutro isolado ou aterrado por alta impedância (HRG).

Física do Deslocamento de Neutro e Sobretensões Temporárias (TOV)

O deslocamento de neutro é um fenômeno eletrodinâmico que ocorre em sistemas trifásicos não solidamente aterrados quando a simetria das admitâncias de fase para a terra é rompida ou ocorre uma falta monofásica à terra (SLG - Single Line-to-Ground). Considere-se um sistema de distribuição minerador em média tensão alimentando cargas não lineares onde o neutro do sistema $N$ está conectado à terra através de uma admitância $Y_n = 1/RNGR + j\omega C_n$.

Aplicando a lei das correntes de Kirchhoff no nó de neutro $N$ em relação ao potencial de terra de referência $G$, a tensão de deslocamento de neutro ($VNG$) diante de um desequilíbrio de admitâncias fasoriais por fase ($Y_A, Y_B, Y_C$) e tensões geradas ($E_A, E_B, E_C$) é expressa rigorosamente como:

Onde a admitância de cada fase $i \in \{A, B, C\}$ inclui a capacitância distribuída do sistema de cabos subterrâneos/blindados ($C_{0,i}$), as admitâncias de isolamento e os componentes de filtro passivo conectados a essa fase:

Sob uma condição de falta monofásica à terra de alta impedância ou franca na fase $A$ ($R_f \to 0$), a admitância $Y_A \to \infty$. O limite da equação de deslocamento de neutro conduz a $VNG \to -E_A$. Portanto, o vetor de tensão neutro-terra iguala-se à magnitude da tensão fase-neutro pré-falta, mas com sinal oposto. As tensões fasoriais das fases não faltosas ($VBG$ e $VCG$) experimentam uma elevação súbita atingindo a tensão fase-fase pré-falta ($VLL$):

O fator de sobretensão temporária ($kTOV$) para as fases não faltosas em um sistema HRG é definido mediante a razão entre a tensão de pico fase-terra durante a falta e a tensão de pico nominal fase-neutro:

No domínio das componentes simétricas, a sobretensão temporária por deslocamento de neutro é analisada avaliando-se as redes de sequência positiva ($Z_1$), sequência negativa ($Z_2$) e sequência zero ($Z_0$). A relação de tensões durante uma falta monofásica à terra é dada por:

Quando a relação $X_0 / X_1 \gg 3$ e $R_0 / X_1 \gg 1$ (característica de sistemas não aterrados ou aterrados por alta resistência), o sistema é classificado como não efetivamente aterrado segundo a norma IEEE 142. Sob estas condições, o coeficiente $kTOV$ não se limita a $1,38$ p.u. (limite de redes efetivamente aterradas), mas atinge níveis superiores a $1,732$ p.u. em regime permanente, com transientes de oscilação LC de alta frequência que podem alcançar picos de:

A presença de capacitâncias distribuídas elevadas em alimentadores de mineração subterrânea ou cavas a céu aberto (cabos de arrasto SHD-GC de alta capacitância linear $C_0 \approx 0,3 - 0,6 \, \mu F/km$) combinada com a indutância não linear de transformadores de potencial (TPs) ou transformadores de força em vazio produz um acoplamento ferroressonante. A ferroressonância paralelo é desencadeada quando a reatância capacitiva de sequência zero ($XC0 = 1 / (\omega C_0)$) se iguala à reatância magnetizante não linear do transformador saturado ($X_m(i)$):

Sob ferroressonância, o neutro oscila caoticamente, gerando sobretensões extremas sustentadas ($TOV > 3,0 p.u.$) com frequências subarmônicas ($1/3 f_1$, $1/2 f_1$) ou harmônicas que superam frequentemente a tensão suportável dielétrica do ar corrigida pela altitude ($V_b / K_a$).

Harmônicos de Cargas Mineradoras e Fator K segundo a NFPA 70B e IEEE C57.110

As plantas de processamento mineral operam com cargas não lineares de alta potência, tais como inversores de frequência de média tensão (VFDs tipo VSI com topologia neutral point clamped NPC de 3 níveis, 6/12/18 pulsos ou AFE), acionamentos de moinhos SAG e de bolas acionados por cicloconversores ou motores síncronos alimentados por inversores LCI. Estas cargas injetam correntes harmônicas de ordem espectral:

Onde $p$ é o número de pulsos do conversor e $k \in \{1, 2, 3, \dots\}$. Para retificadores de 6 pulsos, $h \in \{5, 7, 11, 13, 17, 19, \dots\}$. Em sistemas com desequilíbrios ou trifásicos a 4 fios, emergem harmônicos triplens ($h = 3, 9, 15, 21, \dots$), os quais possuem natureza de sequência zero ($Z_0$).

As correntes harmônicas triplens não se cancelam no nó de neutro dos transformadores em conexão estrela com neutro acessível; pelo contrário, somam-se aritmeticamente no condutor de neutro:

A circulação de correntes harmônicas pelos enrolamentos dos transformadores causa um aquecimento severo devido ao incremento não linear das perdas por correntes de Foucault no cobre (eddy current losses, $PEC$) e perdas suplementaires em partes estruturais (stray losses, $POS$). O Fator K (K-Factor), definido na norma IEEE C57.110 e adotado na norma NFPA 70B para a avaliação de integridade do isolamento em manutenção preventiva/preditiva, quantifica o efeito de aquecimento de correntes harmônicas em relação à corrente fundamental:

Onde $I_h$ é a corrente RMS da harmônica de ordem $h$, e $I_1$ é a corrente RMS fundamental. O fator de perdas por correntes de Foucault ($FHL$) utilizado para realizar o derating da capacidade nominal do transformador é expresso como:

Onde $THD_I = \frac{\sqrt{\sum_{h=2}^{\infty} I_h^2}}{I_1}$ é a Distorção Harmônica Total de corrente. A capacidade de carga reduzida do transformador em pu (per-unit, $Pmax(pu)$) sob ambiente harmônico e de grande altitude é recalculada mediante a combinação da IEEE C57.110 e do fator de altitude:

Onde $PLL(pu)$ são as perdas totais a plena carga calculadas em pu, e $P_{EC-R}(pu)$ são as perdas por correntes de Foucault no enrolamento calculadas à frequência fundamental em condições nominais.

Segundo a norma NFPA 70B (Recommended Practice for Electrical Equipment Maintenance), a acumulação de distorção harmônica elevada combinada com elevações de temperatura acelera o processo de degradação do isolamento térmico mediante a equação de Arrhenius para a vida útil esperada do isolamento ($\mathcal{L}$):

Onde $E_a$ é a energia de ativação do dielétrico, $R$ a constante dos gases e $Thotspot$ a temperatura do ponto mais quente do enrolamento (hotspot), a qual se eleva dramaticamente segundo:

A NFPA 70B exige para instalações industriais em grandes altitudes a medição obrigatória mediante termografia infravermelha ajustada pela transmitância atmosférica corrigida pela altitude, análise de gases dissolvidos no óleo (DGA - Dissolved Gas Analysis) com limiares reduzidos de formação de etileno e acetileno devido a descargas parciais induzidas por TOV, e a verificação periódica do Fator K do sistema de distribuição.

Análise Forense de Falhas Eletromecânicas Integradas

A combinação do deslocamento de neutro, das sobretensões temporárias (TOV), da atenuação dielétrica por baixa pressão atmosférica e do estresse térmico harmônico desencadeia mecanismos de falha multifásicos nos equipamentos da mina. A seguir, realiza-se uma análise forense detalhada por tipo de ativo crítico:

Isolamento em Cabos de Potência (XLPE / EPR)

Em cabos de alta e média tensão (por exemplo, 13,8 kV ou 34,5 kV isolados em XLPE), um deslocamento de neutro sustenta uma tensão fase-terra de $\sqrt{3} VLN$ nos condutores das fases sãs durante horas caso o sistema HRG não elimine a falta rapidamente. A tensão de pico aplicada ao dielétrico interno é:

Esta sobretensão reduz drasticamente a Tensão de Incepção de Descargas Parciais (CIV - Partial Discharge Inception Voltage). A presença de microcavidades (voids) dentro do isolante XLPE gera um gradiente de campo elétrico local $Evoid = \epsilon_{r} \cdot E_0$. Quando $Evoid$ supera a tensão de ruptura do gás contido dentro da cavidade (cuja rigidez cai segundo a lei de Paschen se houver migração de pressão ou gradiente de altitude em terminações), iniciam-se descargas parciais sustentadas. Estas descargas erosionam a matriz polimérica formando arborescências elétricas (electrical treeing), acelerando a ruptura dielétrica catastrófica do cabo.

Painéis de Média Tensión (Switchgear) e Isolamento no Ar

O ar dentro do invólucro do switchgear atua como meio isolante primário entre os barramentos de fase e a estrutura aterrada. A 4.200 m.a.n.m., a rigidez dielétrica do ar reduz-se em um fator aproximado de $\delta \approx 0,62$. Durante um evento de TOV associado a uma falta à terra em um sistema HRG com elevada presença de 3.º e 9.º harmônicos:

• A distância de isolamento no ar (clearance) requerida para suportar impulsos de manobra ou transientes de reabertura de arco deve ser incrementada em $45\% - 60\%$ segundo a norma IEC 60664-1 Tabela A.2.
• Se o painel foi projetado para operação ao nível do mar ($V_{\text{disrupción}} = 95 kV BIL$), a 4.200 m.a.n.m. a tensão suportável BIL real desce para $95 \times 0,62 = 58,9 kV BIL$.
• Um transiente de TOV com pico de 48 kV fase-terra sobreposto a harmônicos de alta frequência ultrapassa o limite dielétrico do ar degradado, provocando a abertura de um arco fase-terra que evolui rapidamente para um curto-circuito trifásico destrutivo por ionização massiva do volume do painel (evento de Arc Flash).

Transformadores de Força e Fator K

Em transformadores, o aumento da tensão de neutro eleva o ponto de operação da curva de magnetização do núcleo ($B-H$). Se a TOV elevar a densidade de fluxo $B$ acima do ponto de saturação ($Bsat \approx 1,7 - 1,8 Tesla$):

A corrente de magnetização converte-se em impulsos pontiagudos ricos em harmônicos ímpares e triplens ($Imag \propto B^{n}, n > 9$). Isto incrementa exponencialmente as perdas no ferro ($P_Fe \propto f \cdot B^2 + f^2 \cdot B^2$) e induz maiores perdas por correntes de Foucault ($PEC$). O calor não dissipado devido ao $kalt$ degrada termicamente o papel isolante celulósico, liberando monóxido e dióxido de carbono ($CO, CO₂$), e gerando pontos quentes (hotspots) que induzem a formação de compostos furânicos, reduzindo o Grau de Polimerização (DP) do papel de 1000 para valores críticos (< 200), o que assinala o fim da vida útil mecânica e a falha iminente entre espiras do enrolamento.

Parâmetro Eletromecânico / Isolamento	Condição Padrão (0 m.a.n.m., $THD_I < 5\%$)	Condição Severa Mineradora (4.500 m.a.n.m., $K \ge 13, TOV = 1,73 \, pu$)	Consequência Térmica e Dielétrica Imposta
Rigidez Dielétrica do Ar ($E_b$)	$3,0 kV/mm$	$\sim 1,85 kV/mm$ ($\delta \approx 0,61$)	Ruptura dielétrica fase-terra em painéis de média tensão e isoladores.
Tensão Fase-Terra em Fases Sãs ($VLG$)	$1,0 p.u.$ ($VLL / \sqrt{3}$)	$1,732 p.u.$ sustentada (ou $> 2,5 p.u.$ transiente)	Envelhecimento acelerado do XLPE por descargas parciais (treeing).
Fator K de Perdas Harmônicas	$K = 1,0$	$K = 13,0 - 20,0$	Perdas $PEC$ incrementadas em $1300\% - 2000\%$ nos enrolamentos.
Temperatura do Ponto Quente ($Thotspot$)	$98^\circ C - 110^\circ C$	$> 145^\circ C$ (sem derating específico)	Destruição da cadeia polimérica do papel Kraft (DP < 200).
Tensão de Incepção de Descargas Parciais (CIV)	$100\%$ do valor nominal	$\sim 60\% - 65\%$ do valor nominal	Geração sustentada de efeito coroa e degradação do isolamento externo.
Capacidade de Dissipação dos Radiadores ($h_c$)	$100\%$ de eficiência convectiva	$\sim 70\%$ de eficiência convectiva	Superaquecimento generalizado do fluido isolante/refrigerante.

Estratégias Avançadas de Mitigação e Critérios de Projeto

Para garantir a continuidade operacional e a segurança integral da rede elétrica mineradora sob condições combinadas de grandes altitudes, distorção harmônica e risco de TOV por deslocamento de neutro, devem ser implementadas as seguintes metodologias de projeto e engenharia de detalhe:

Dimensionamento do Resistor de Aterramento de Neutro (NGR) para Controle de TOV e Altitude

Em sistemas HRG em grandes altitudes, a corrente escolhida para o resistor de aterramento ($INGR$) deve ser estritamente maior do que a corrente total de carga capacitiva da rede em sequência zero ($IC0$):

Onde $C_{0, total}$ é a capacitância acumulada de todos os cabos, bancos de capacitores e capacitâncias parasitas da rede. Ao garantir $RNGR \le \frac{1}{3 \omega C_{0, total}}$, dissipa-se a energia capacitiva armazenada durante a falta monofásica, amortecendo completamente os picos de sobretensão por reabertura de arco (arcing grounds) e limitando $kTOV \le 1,732 p.u.$.

Adicionalmente, o elemento físico R (banco de resistores de liga de aço inoxidável ou cromo-níquel) deve ser dimensionado considerando o derating térmico por altitude para a dissipação do calor gerado durante uma falta sustentada ($tfalla = 10 s$ ou contínuo segundo a NFPA 70B):

Seleção e Coordenação de Isolamento de Para-Raios (Surge Arresters)

A seleção de para-raios de óxido de zinco (ZnO) em instalações mineradoras em grandes altitudes sujeitas a TOV requer uma avaliação rigorosa da margem entre a Tensão de Operação Contínua Máxima (MCOV ou $V_c$) e a curva de capacidade TOV do para-raios ($VTOV(t)$). A MCOV deve ser superior à tensão máxima fase-terra na presença de deslocamento de neutro:

O para-raios deve ser capaz de absorber a energia de sobretensão sem sofrer disparo térmico (thermal runaway). A capacidade de absorção de energia específica ($E_K$ em kJ/kV de MCOV) é desclassificada por altitude devido à refrigeração reduzida do invólucro polimérico:

Da mesma forma, a distância de escoamento externa do para-raios é dimensionada aplicando-se um fator mínimo de $31 mm/kV$ de tensão máxima do sistema, corrigido pelo fator $K_a$ para ambientes com contaminação mineradora severa (poeira condutora metálica ou salina).

Filtragem Harmônica e Desassintonia para Evitar Ressonância Paralelo

Para mitigar a injeção do espectro harmônico e prevenir o superaquecimento de transformadores por Fator K, são projetados filtros passivos desassintonizados (detuned filters) ou filtros tipo C (C-type damped filters). A frequência de sintonia ($fst$) deve ser posicionada abaixo da quinta harmônica ($h=5$), tipicamente em $h_r = 4,2 - 4,7$, para evitar ressonância paralelo entre os bancos de capacitores e a indutância da rede ($Lsys + Ltrafo$):

Em presença de inversores de frequência de alta potência, a implementação de Filtros Ativos de Potência (APF) de topologia Multinível conectados em paralelo injeta em oposição de fase a corrente harmônica exata requerida pela carga ($IAPF(t) = -Ih(t)$), garantindo que $THD_I < 5\%$ no Ponto de Acoplamento Comum (PCC) segundo a norma IEEE 519, reduzindo efetivamente o Fator K visto pelo transformador de alimentação para $K \approx 1,0$.

Aplicação e Metodologia de Cálculo com Vexten Suite

A plataforma avançada de engenharia elétrica Vexten Suite integra módulos de simulação numérica e analítica projetados para resolver simultaneamente as equações complexas de curto-circuito, derating de isolantes por altitude, fluxo harmônico e capacidade de condução de cabos em ambientes críticos como a mineração em grandes altitudes. A seguir, detalha-se a metodologia de cálculo e o fluxo de projeto aplicado mediante os módulos do Vexten Suite.

Módulo de Curto-Circuito e Impedâncias Vexten (IEC 60909 / IEEE 141)

O módulo Vexten Short-Circuit & Grounding Analyzer calcula a matriz de impedâncias de sequência ($Z_1, Z_2, Z_0$) ajustada pela temperatura operacional dos condutores e pela resistividade do solo em grande altitude. Durante uma falta monofásica à terra em uma rede HRG mineradora, o software resolve a corrente inicial de curto-circuito simétrica ($Ik1''$) e ignora os efeitos da impedância de carga convencional segundo a IEC 60909-0:

Dado que em uma rede HRG o termo de impedância de sequência zero é dominado pela resistência do neutro ($Z_0 \approx 3 RNGR$), a corrente de falta à terra é aproximada no Vexten por:

Simultaneamente, o motor do Vexten calcula o perfil fasorial de tensões nos barramentos sãos, computando diretamente o fator $k_{TOV, bar}$ e avaliando o deslocamento do neutro $VNG$ para cada cenário de operacionalidade da mina.

Módulo de Dimensionamento de Cabos e Incremento por Harmônicos Vexten (IEC 60287 / NEC 310)

O módulo Vexten Cable Sizing & Harmonic Derating Engine calcula a capacidade de condução de corrente efetiva (ampacidade, $Iamp$) de cabos unipolares ou tripolares de média e baixa tensão (por exemplo, tipo SHD-GC ou XLPE) sob a ação simultânea da altitude e da corrente de neutro harmônica mediante a equação multidimensional de derating:

Onde o fator de degradação por correntes harmônicas no neutro ($kharm$) é determinado internamente no Vexten avaliando-se o aquecimento do condutor neutro e dos condutores de fase de acordo com a norma NEC 310.15(E) e IEC 60287-1-1. Se a porcentagem de harmônicos de terceira ordem ($h=3$) superar $33\%$, a corrente no neutro supera a corrente de fase, e o Vexten ajusta automaticamente a seção do condutor de neutro para $150\% - 200\%$ da seção da fase, recalculando o gradiente térmico camada por camada:

Onde $T_1, T_2, T_3, T_4$ representam as resistências térmicas das camadas do cabo (isolamento, cobertura, leito e solo/ar) modificadas no código do Vexten segundo a densidade volumétrica do ar na altitude inserida pelo usuário.

Módulo de Fator K e Análise de Ressonância Vexten (IEEE C57.110 / IEEE 519)

O módulo Vexten Harmonic & Resonance Analyzer realiza a varredura de impedância espectral do sistema ($Z(f)$) desde $50 Hz$ até $2500 Hz$ em passos de $1 Hz$. O algoritmo detecta os pontos de inflexão onde a derivada da impedância muda de sinal ($\frac{d|Z|}{df} = 0$), identificando frequências de ressonância paralelo ($f_p$) e série ($f_s$).

A partir do espectro harmônico de corrente inserido ou exportado de medidores de qualidade da energia elétrica (arquivos PQDIF/COMTRADE), o Vexten calcula automaticamente:

1. O Fator K do sistema segundo a IEEE C57.110.
2. O fator de perdas por correntes de Foucault $FHL$.
3. A capacidade de carga com derating do transformador em kVA ($kVA_{derated} = kVA_{nominal} \cdot Pmax(pu)$).
4. O perfil de elevação de temperatura do ponto mais quente (hotspot, $Thotspot$) e o índice de perda de vida útil segundo o modelo de Arrhenius da NFPA 70B.

Algoritmo de Cálculo Integrado no Vexten Suite

O procedimento a seguir descreve o fluxo de cálculo computacional executado pelo Vexten Suite para resolver a interação entre o isolamento por altitude, o deslocamento de neutro e a distorção por Fator K:

Mediante esta arquitetura computacional, o Vexten Suite provê uma ferramenta integral de alta precisão que elimina os riscos de falhas por sobretensões temporárias, arcos elétricos dielétricos e colapsos térmicos de isolamento em sistemas de potência mineradores de alta complexidade em grandes altitudes.