Variação Sazonal da Resistência de Aterramento por Seca ou Congelamento e Tratamento com Géis Condutivos

Termodinâmica e Transporte Iônico do Solo: Mecanismos Físicos da Variação Sazonal

A resistividade elétrica do solo ($\rho$) não é um parâmetro estático de projeto; trata-se de uma variável termodinâmica e eletroquímica altamente dependente do estado higrométrico, da concentração de sais dissolvidos, da porosidade da matriz geológica e da temperatura absoluta do meio. Em sistemas de potência, subestações elétricas e usinas de geração, desconsiderar as oscilações sazonais introduz discrepâncias críticas entre os cálculos nominais baseados em normas como a IEEE Std 80 ou a IEC 60364-5-54 e o comportamento real do sistema de aterramento (SPAT) diante de eventos de curto-circuito ou descargas atmosféricas.

Condução eletrolítica, porosidade e saturação hídrica (Equação de Archie generalizada)

A condução de corrente elétrica em meios porosos não saturados ocorre predominantemente através da fase líquida contida no espaço poroso intersticial mediante transporte iônico migratório. A matriz mineral de quartzo, feldspatos ou silicatos atua usualmente como um isolante dielétrico com resistividades intrínsecas superiores a $10^6\ \Omega\cdot m$. A condutividade global do terreno ($\sigma = 1/\rho$) é modelada por meio de uma extensão da Lei de Archie para meios não saturados e multicomponentes:

Onde:

• $\rho_w$: Resistividade do eletrólito poroso ($\Omega\cdot m$), governada pela atividade e mobilidade dos íons ($Na ^+$, $Ca ^{2+}$, $Mg ^{2+}$, $Cl ^-$, $SO _4^{2-}$, $HCO _3^-$).
• $\phi$: Porosidade total do solo (fração volumétrica de vazios, $0 < \phi < 1$).
• $S_w$: Grau de saturação de água nos poros ($V_{\text{água}} / Vvazios$, $0 \le S_w \le 1$).
• $a$: Fator de tortuosidade lítica (tipicamente entre 0,5 e 1,5).
• $m$: Fator de cimentação da matriz porosa ($1,3 \le m \le 2,5$).
• $n$: Expoente de saturação hídrica (geralmente $n \approx 2,0$).

Sob condições de estiagem severa, a sucção mátrica ($\psi_m$) ascende exponencialmente conforme a curva de retenção de água no solo (SWRC, na sigla em inglês), expelindo a água capilar livre. À medida que $S_w \to 0$, a continuidade percolativa das trajetórias de condução se rompe abruptamente (limiar de percolação), forçando a resistividade aparente do terreno ($\rho_{app}$) a aumentar em até 2 a 3 ordens de grandeza em relação ao seu estado de saturação de projeto.

Termocriogenia do terreno: Transição de fase da água e efeito do permafrost sazonal

O declínio térmico abaixo do ponto de congelamento ($T < 0\ ^\circ C$) desencadeia uma transição de fase termodinâmica de primeira ordem na água livre dos macroporos. O gelo possui uma estrutura cristalina hexagonal fechada na qual a mobilidade protônica pelo mecanismo de Grotthuss é extremamente baixa, exibindo uma resistividade intrínseca na ordem de $10^5$ a $10^7\ \Omega\cdot m$.

Onde $\theta_u(T)$ representa o conteúdo volumétrico de água líquida não congelada remanescente nas microporosidades à temperatura $T$, $\theta_0$ é o conteúdo de água total prévio ao congelamento, $T_f$ é a temperatura de congelamento rebaixada pela salinidade porosa, e $\beta, \gamma$ são coeficientes empíricos texturais. A presença de camadas de gelo superficiais não apenas anula o transporte por difusão iônica, mas também induz o fenômeno de criossucção, atraindo umidade de horizontes inferiores em direção à frente de congelamento ativo, dessecando e enrijecendo dieletricamente as camadas subjacentes imediatas ao eletrodo.

Dependência térmica da mobilidade iônica e Lei de Arrhenius

Acima do ponto de congelamento ($T > 0\ ^\circ C$), a resistividade do eletrólito é controlada pela viscosidade dinâmica ($\eta$) do solvente aquoso e pela difusividade dos íons segundo a relação de Nernst-Einstein. A variação térmica da resistividade pode ser expressa formalmente por meio do modelo semilinearizado ou pelo formalismo de Arrhenius:

Onde $\alpha_T$ é o coeficiente de temperatura da resistividade do solo (tipicamente $\alpha_T \approx 0,02$ a $0,03\ K ^{-1}$ para eletrólitos diluídos), $E_a$ é a energia de ativação para o transporte iônico ($\approx 0,15 - 0,22\ eV$), e $k_B$ é a constante de Boltzmann. Uma redução de temperatura de $25\ ^\circ C$ para $1\ ^\circ C$ aumenta a resistividade do solo em mais de 70% exclusivamente pela redução da agitação térmica e da cinética de deriva eletroforética.

Modelagem Matemática da Resistividade Equivalente e Parâmetros da Malha de Terra

As oscilações climatológicas modificam o solo em profundidade, conformando estratigrafias dinâmicas. A modelagem não pode ser considerada homogênea ($\rho = constante$), mas sim formulada como um problema de estratificação horizontal bicamada ou multicamada dependente do tempo.

Estratificação sazonal do solo: Perfil bicamada dinâmico variável no tempo

Durante a temporada de estiagem ou de congelamento, a camada superficial de profundidade $h_1$ experimenta uma alteração extrema em sua resistividade ($\rho_1(t)$), enquanto os estratos profundos ($\rho_2$) permanecem relativamente protegidos das influências meteorológicas devido à inércia térmica e hidrológica da crosta. O coeficiente de reflexão eletromagnético interfacial $K$ é definido analiticamente por:

O potencial eletrostático $V(r)$ na superfície de um solo bicamada devido a uma injeção pontual de corrente $I$ é obtido resolvendo-se a equação de Laplace ($\nabla^2 V = 0$) em coordenadas cilíndricas mediante a integral de Hankel:

Quando ocorre um congelamento severo ou seca superficial extrema, $\rho_1 \gg \rho_2$, implicando que $K \to -1$. Por outro lado, em um leito rochoso subjacente com superfície umedecida, $K \to +1$. Em ambos os casos limite, as séries infinitas de coeficientes de reflexão condicionam a distribuição dos gradientes de potencial na superfície e a dispersão tridimensional de corrente para as camadas profundas da terra.

Fórmulas analíticas de resistência de aterramento sob condições limite (Schwarz, Sverak e IEEE 80)

Para uma malha combinada com condutores horizontais e hastes verticais embutidas em um solo não homogêneo, o método analítico rigoroso de C. W. Schwarz calcula a resistência total de aterramento ($R_g$) considerando o acoplamento mútuo entre a retícula horizontal ($R_1$) e as hastes verticais ($R_2$):

Onde as componentes são definidas a partir do comprimento total do condutor da malha ($L_c$), da área ocupada ($A$), da profundidade de enterramento ($h$), do raio do condutor ($r$), do número de hastes verticais ($n_R$) e do comprimento unitário da haste ($L_r$):

Onde $k_1, k_2$ são coeficientes geométricos de forma da malha, $b$ é o raio da haste vertical e $\rho_{app}$ é a resistividade aparente ponderada que deve sintetizar a degradação da camada superior. Se a camada superior de espessura $h_1 \ge h$ congela ou desseca, a contribuição da componente $R_1$ fica virtualmente isolada ($R_1 \to \infty$), forçando que a totalidade da corrente $I_G$ seja dissipada através do subsistema de hastes verticais $R_2$, superdimensionando a densidade de corrente superficial $J_s$ e elevando drasticamente a resistência equivalente global.

Cálculo de potenciais de passo, toque e GPR crítico durante eventos climáticos adversos

De acordo com a norma IEEE Std 80, a tensão de passo tolerável ($Estep$) e a tensão de toque tolerável ($Etouch$) para um ser humano de 50 kg ou 70 kg são governadas pelo fator de redução superficial $C_s$, o qual depende da resistividade da camada superficial do terreno ($\rho_s$) e da espessura da brita/camada de acabamento ($h_s$):

Durante épocas de seca, a desidratação da subcamada imediata ($\rho \gg \rho_s$) degrada o valor do potencial de malha ($E_m$) e de passo ($E_s$). O incremento sazonal da resistência global da malha $R_g(t)$ eleva a Elevação de Potencial de Terra (GPR, Ground Potential Rise):

Onde $I_G$ é a corrente simétrica de falta à terra de projeto, $D_f$ é o fator de decremento por assimetria da componente de corrente contínua (CC subtransitória), e $S_f$ é o fator de divisão de corrente que drena pelos cabos para-raios (cabos guarda). Um aumento sazonal de $R_g(t)$ por um fator de $\times 4$ a $\times 10$ faz com que o GPR exceda amplamente as tensões suportáveis dielétricas dos isolamentos secundários, cabos de controle e blindagens de telecomunicações acopladas ao SPAT.

Análise Forense de Falhas Elétricas Induzidas por Degradação Sazonal do Terreno

A falha de um sistema de aterramento devido a variações climáticas não se manifesta de forma linear ou previsível por meio de medições estáticas tradicionais na primavera ou no outono; sua manifestação é destrutiva e cataclísmica durante faltas assimétricas monofásicas ($1\Phi-T$) ou descargas atmosféricas diretas coincidentes com o pico de estiagem estival ou de congelamento invernal.

Perda de referência de neutro, sobretensões dinâmicas temporárias (TOV) e falha de transformadores

Em sistemas com transformadores de potência conectados em estrela solidamente aterrada (Yg), a impedância de sequência zero do sistema ($Z_0$) está intrinsecamente ligada à impedância do eletrodo de aterramento do neutro ($Z_g = R_g + j X_g$).

Sob uma falta à terra em uma das fases (fase A) em um regime de solo degradado no qual $R_g \gg X_1$, a relação de sequências se altera drasticamente: $R_0 / X_1 > 2$ e $X_0 / X_1 > 3$. O sistema deixa de operar como um sistema solidamente aterrado e transita dinamicamente para um regime de neutro isolado ou aterrado por alta impedância. A sobretensão temporária induzida nas fases sadias (fases B e C) é descrita vetorialmente como:

Esta condição de sobretensão sustentada na frequência industrial ($TOV = \sqrt{3}\ p.u.$) excede a máxima tensão de operação contínua ($U_c$ / MCOV) dos para-raios instalados e supera a rigidez dielétrica do isolamento menor de transformadores, enrolamentos terciários e terminais de cabos de MT/AT, provocando perfurações dielétricas destrutivas por efeito corona e avalanches de descargas parciais.

Inoperância ou recalcitrância de para-raios e descarregadores de sobretensão (Surge Arresters / DPS)

Os para-raios de óxido de zinco ($ZnO$) dependem de uma trajetória de impedância transitória ultrabaixa para a terra a fim de escoar a onda de corrente impulsiva ($8/20\ \mu s$ ou $10/350\ \mu s$). A impedância transitória de impulso $Z_p(t)$ de uma haste vertical sob congelamento severo ou solo seco não se comporta de forma puramente resistiva; ela apresenta uma forte indutância série parasita ($L_p$) combinada a uma resistência impulsiva $R_i$ agravada pela escassez de portadores iônicos livres:

Quando o solo superficial está congelado ou seco, o campo elétrico crítico de ionização do solo ($E_0 \approx 300\ kV/m$) necessário para formar canais de microfaíscas de plasma que reduzam a resistência impulsiva não é alcançado localmente. Consequentemente, o termo $i(t) \cdot R_i(t)$ assume valores da ordem de centenas de quilovolts, transferindo a totalidade da frente de onda da sobretensão para o isolamento de barramentos da subestação, cubículos blindados em $SF _6$ ou transformadores nominalmente protegidos.

Falha de coordenação de proteções de falta à terra (ANSI 50N/51N, 67N) por elevação da impedância de retorno

A corrente de curto-circuito monofásica à terra ($Ik1''$) calculada conforme a norma IEC 60909 é formulada analiticamente por meio de suas impedâncias de sequência direta ($Z_1$), inversa ($Z_2$) e zero ($Z_0$):

Se a resistência sazonal do SPAT da subestação ou das estruturas da linha de transmissão se eleva drasticamente por congelamento ou dessecação ($R_g(t) \gg Z_1$), a magnitude da corrente $Ik1''$ é reduzida severamente para valores abaixo do limiar de partida (pickup) dos relés de sobrecorrente temporizados ou instantâneos (ANSI 51N/50N). A falta converte-se em um defeito de alta impedância (Hi-Z), permanecendo sustentada no sistema, o que induz arcos elétricos intermitentes, ignição de vegetação circundante e a elevação prolongada de potenciais de passo e toque letais no perímetro externo das instalações.

Matriz Comparativa de Normalização e Limites Críticos de Operação

A seguinte matriz técnica condensa os limites operacionais, as implicações de segurança e as condições de falha segundo os principais padrões internacionais de engenharia elétrica frente à deterioração por condições ambientais extremas.

Parâmetro / Norma Técnica	Limite Padrão Nominal (Condição Base)	Comportamento Sob Estiagem Severa ($S_w < 0,1$)	Comportamento Sob Congelamento ($T < -10\ ^\circ C$)	Modo de Falha Forense / Consequência Operacional
Resistência de Aterramento ($R_g$) IEEE Std 81 / IEEE Std 142	Subestações AT: $R_g \le 1,0\ \Omega$ Distribuição/MT: $R_g \le 5,0\ \Omega$	Aumento de $300\%$ a $1200\%$ ($3,0\ \Omega \le R_g \le 15,0\ \Omega$)	Aumento de $800\%$ a $3500\%$ ($10,0\ \Omega \le R_g \le 80,0\ \Omega$)	Violação do GPR; inoperância de esquemas de disparo de neutro; disrupção inversa (backflashover) em cadeias de isoladores.
Tensão de Toque Tolerável ($Etouch$) IEEE Std 80 / IEC 61936-1	$Etouch \le \frac{1000 + 1,5 C_s \rho_s}{\sqrt{t_s}}$ Base: $\rho_s \approx 2500\ \Omega\cdot m$ (brita)	Aparente elevação superficial na brita seca, porém colapso do perfil subsuperficial: $E_m$ supera o limite de segurança.	A camada congelada satura o fator de redução $C_s \to 0$. Queda drástica da tensão admissível para o corpo humano.	Fibrilação ventricular e óbito de operadores durante manobras de disjuntores ou seccionadoras no pátio de AT.
Limite de Atuação ANSI 51N IEC 60255-151 / IEEE C37.112	Sensibilidade ajustada entre $10\% - 20\%$ da corrente nominal de carga ou disparo rápido.	$Ik1''$ cai abaixo do limiar de partida devido à alta impedância do laço de sequência zero.	A corrente de falta decresce a ponto de se confundir com desequilíbrios normais de fase por carga assimétrica.	Falta não eliminada; dano térmico irreversível em transformadores de força por TOV prolongada; colapso de cabos.
Resistência de Isolamento e MCOV de Para-raios IEC 60099-4 / IEEE C62.11	$V_{ residual } \le BIL _{ equipamento } / 1,4$ Tensão de operação contínua $U_c \ge 1,05 VLN$	A TOV dinâmica eleva a tensão nos terminais para $\sqrt{3} VLN$, ultrapassando a curva de suportabilidade de TOV (TOV capability curve).	A impedância de dispersão $R_i$ dispara; o para-raios não absorve energia adequadamente, descarregando sobre o enrolamento primário.	Explosão termomecânica dos invólucros poliméricos ou de porcelana dos para-raios de óxido de zinco ($ZnO$).
Gradiente de Tensão Perimétrico IEC 62305 / NFPA 780	Separação equipotencial com $Vstep < 5\ kV$ na borda da malha a $1\ m$.	Expansão geométrica do cone de potencial para o exterior da subestação.	A camada superficial congelada força a emersão de linhas de campo elétrico altamente concentradas lateralmente.	Descargas disruptivas em direção a cercas perimétricas metálicas externas, eletrocutando transeuntes ou animais.

Físico-Química dos Géis e Compostos de Tratamento Condutivo

Para mitigar a variabilidade climática sazonal, utilizam-se compostos de melhoria eletroquímica que modificam localmente a interface eletrodo-solo, estabilizando o teor de umidade, otimizando a concentração iônica e reduzindo a impedância de contato metal-solo.

Classificação de agentes: Bentonita sódica, géis à base de poliacrilamidas e matrizes de carbono liofilizadas

Existem diferenças estruturais e eletrodinâmicas críticas entre os diversos agentes condicionadores empregados na engenharia:

• Bentonita Sódica (Montmorilonita Hidratada): Argila mineral filossilicatada composta por camadas octaédricas de alumina intercaladas entre duas camadas tetraédricas de sílica ($2:1$). Possui a capacidade de expandir de 12 a 15 vezes seu volume seco original na presença de água livre. Sua condutividade é predominantemente iônica por troca de cátions $Na ^+$. Contudo, sob condições de seca térmica persistente, sofre contração volumétrica irreversível (fissuramento reticular), perdendo o contato físico-elétrico direto com o eletrodo metálico e elevando drasticamente a resistência de contato.
• Géis de Reticulação Polimérica (Poliacrilamidas e Hidrogéis Iônicos): Polímeros superabsorventes formados por cadeias hidrofílicas tridimensionais unidas por ligações cruzadas covalentes. Retêm soluções eletrolíticas concentradas de sais higroscópicos ($MgSO _4$, $CaCl _2$, silicatos de potássio) por meio de pontes de hidrogênio. Não apresentam retração mecânica volumétrica durante a dessorção hídrica e retêm a umidade sob tensões mátricas de até $1,5\ MPa$, impedindo o congelamento a temperaturas inferiores a $-18\ ^\circ C$ devido ao rebaixamento crioscópico propiciado pelo soluto aprisionado na matriz.
• Cimentos Condutivos e Matrizes Carbonosas à Base de Grafite: Compostos de cura hidráulica à base de cimento Portland dopados com micropartículas de grafite cristalino dessulfurado ($C \ge 99\%$). A condução é predominantemente eletrônica (tunelamento quântico / teoria da percolação quântica) e não iônica, com resistividades intrínsecas secas de $0,02$ a $0,1\ \Omega\cdot m$. Não dependem da água porosa para conduzir, tornando-se imunes à dessecação e ao congelamento estrutural, além de fornecerem uma barreira física permanente contra a corrosão galvânica.

Higroscopia, capacidade de troca catiônica (CTC) e cinética de dessorção hídrica

A Capacidade de Troca Catiônica (CTC) quantifica a densidade de cargas negativas na superfície da matriz coloidal capazes de reter cátions trocáveis em equilíbrio:

Enquanto um solo arenoso padrão exibe uma $CTC \le 5\ meq /100\ g$, um gel polimérico dopado ou uma bentonita tratada supera $80 - 120\ meq /100\ g$. A cinética de evaporação e retenção de umidade ao longo do tempo $S_w(t)$ sob estresse térmico é modelada pela equação de transporte difusivo não linear de Richards:

Os géis eletroquímicos de última geração modificam localmente a condutividade hidráulica não saturada $K(\theta)$, atuando como uma válvula osmótica reversa que impede a migração ascendente de água por termomigração em direção à superfície evaporativa durante períodos de estiagem extrema.

Corrosão galvânica induzida, pH e passivação de eletrodos de cobre e aço cobreado

Um condicionador químico não deve provocar a degradação acelerada do eletrodo metálico. O potencial eletroquímico de corrosão ($E_{ corr }$) e a densidade de corrente anódica de passivação ($J_{ corr }$) são calculados por meio das curvas de polarização de Tafel:

Para assegurar uma vida útil de projeto superior a 30 anos, em conformidade com a norma IEC 62561-7 (Requisitos para compostos de melhoria de aterramento):

• O pH do gel curado deve situar-se estritamente na região de imunidade / passivação termodinâmica do diagrama de Pourbaix para o cobre: $7,0 \le pH \le 10,5$.
• A concentração lixiviável de íons sulfato ($SO _4^{2-}$) e cloreto ($Cl ^-$) deve ser inferior a $0,05\%$ em massa seca, a fim de evitar a corrosão por pites (pitting corrosion) da camada de cobre em hastes de aço cobreado (copper-clad steel, espessura mínima nominal $\ge 254\ \mu m$).
• A resistividade própria do composto compactado deve ser menor que $0,2\ \Omega\cdot m$ ($20\ \Omega\cdot cm$) na condição saturada e menor que $1,0\ \Omega\cdot m$ sob condições desidratadas de ensaio laboratorial acelerado a $105\ ^\circ C$.

Metodologia de Engenharia, Dimensionamento e Simulação Avançada com o Vexten Suite

A integração da análise numérica avançada de malhas de aterramento com a simulação eletromecânica do sistema elétrico de potência é indispensável para garantir a conformidade normativa ao longo de todo o ciclo sazonal.

Inversão de sondagens Wenner/Schlumberger e parametrização de camadas no Vexten Grounding Engine

A primeira etapa metodológica consiste na execução rigorosa de ensaios de sondagem elétrica vertical (SEV) por meio do arranjo tetraeletródico de Wenner, em conformidade com a IEEE Std 81, utilizando espaçamentos progressivos $a = [1, 2, 4, 8, 16, 32]\ m$. A resistência medida $R(a)$ resulta na resistividade aparente:

No módulo de inversão geofísica não linear do Vexten Grounding Engine, implementa-se o algoritmo de otimização de Gauss-Newton com regularização de Tikhonov para solucionar o problema inverso e obter a estratigrafia multicamada de projeto:

Onde $\mathbf{p} = [\rho_1, \rho_2, \dots, \rho_n, h_1, H₂, \dots, h_{n-1}]^T$ representa o vetor de parâmetros geoelétricos, $\mathbf{L}$ é a matriz de suavização diferencial de primeira ordem, e $\lambda$ é o multiplicador de Lagrange do parâmetro de regularização.

Com o Vexten Suite, o engenheiro parametriza dois cenários envoltórios críticos (worst-case envelopes):

• Cenário de Estiagem / Seca Crítica: Modificação da camada 1 ($h_1 = 1,5\ m$) com a aplicação de um fator multiplicador hidrológico $\kappa_{ dry } \in [3,0;\ 8,0]$.
• Cenário de Inverno / Congelamento Extremo: Modificação da camada superficial ($h_1 = 0,8\ m$) com resistividade criogênica $\rho_1 \ge 10\,000\ \Omega\cdot m$.

Cossimulação de curto-circuito assimétrico IEC 60909 / IEEE 141 acoplada à impedância sazonal da malha

Diferentemente dos fluxos de trabalho desacoplados convencionais, o solucionador matricial de curto-circuito do Vexten Short-Circuit IEC 60909 executa uma cossimulação acoplada e iterativa integrada ao cálculo da malha de terra:

O valor de $R_g(t)$ calculado pelo módulo eletrostático retroalimenta dinamicamente a matriz nodal de admitâncias de sequência zero. Dessa maneira, o software recalcula com extrema precisão o fator de divisão de corrente real ($S_f$), computando a contribuição efetiva dos cabos para-raios de linhas aéreas, blindagens metálicas de cabos isolados dimensionados conforme a IEC 60287 / NBR 5410 / NBR 14039, e as correntes de retorno pelo solo profundo.

Otimização topológica de hastes profundas e valas tratadas com gel condutivo

Quando a simulação em terreno natural indica a violação dos potenciais máximos admissíveis de toque ($E_m > E_{ touch }$), o Vexten Grounding Engine permite formular e iterar a estratégia de mitigação por meio do condicionamento eletroquímico perimétrico e da instalação de poços profundos encapsulados.

O raio equivalente hidrodinâmico-elétrico ($req$) de uma haste vertical de raio metálico $r_0$ instalada em uma perfuração de raio $r_b$ preenchida com gel condutivo de resistividade $\rho_{ gel }$ imersa em um solo de resistividade local $\rho_{ solo }$ é modelado quantitativamente por:

Visto que $\rho_{ gel } \ll \rho_{ solo }$ (por exemplo, $0,1\ \Omega\cdot m$ contra $500\ \Omega\cdot m$), o termo exponencial tende a $e^0 = 1$, fazendo com que $req \to r_b$. De forma analítica, o poço tratado converte o raio de injeção de corrente do raio do núcleo metálico ($r_0 \approx 9,5\ mm$) para o raio geométrico total da perfuração ($r_b \approx 75 - 150\ mm$), reduzindo a resistência intrínseca de contato do eletrodo em mais de $40\%$ a $60\%$ e tornando-o substancialmente independente da dessecação superficial.

Ao forçar mecanicamente a dissipação de corrente em direção à camada inferior estável e protegida ($\rho_2$) através da haste encapsulada em matriz coloidal condutiva, assegura-se a estabilidade do sistema frente a flutuações climatológicas severas, garantindo a operação confiável, a coordenação seletiva das proteções de sobrecorrente e a segurança irrestrita das equipes de campo durante toda a vida útil da subestação.