Ingeniería Eléctrica

Tensões de Passo e Toque em Subestações conforme IEEE 80

El diseño de mallas de puesta a tierra en subestaciones eléctricas no se limita a lograr un valor arbitrario de resistencia inferior a 1 ohmio. Bajo una falla a

Ing. Francisco Ramírez

Fundamentos Eletrofisiológicos e Critérios de Fibrilação Ventricular

O projeto de sistemas de malha de aterramento em subestações elétricas de alta e extra-alta tensão tem como objetivo primário a preservação da vida humana e a integridade dos ativos eletromecânicos. Durante um evento de curto-circuito franco à terra, a injeção massiva de corrente no solo gera gradientes de potencial superficiais que podem submeter o corpo de um operador ou transeunte a diferenças de potencial críticas.

A resposta fisiológica do organismo diante da passagem de corrente alternada em frequência industrial (50/60 Hz) é dominada pela indução de fibrilação ventricular, fenômeno no qual os miócitos ventriculares perdem a sincronia contrátil, anulando o débito cardíaco. A probabilidade de fibrilação depende da magnitude da corrente, da duração da exposição, da trajetória através do corpo e da fase do ciclo cardíaco em que ocorre o choque elétrico (especialmente a fase vulnerável correspondente à onda T no eletrocardiograma).

IB=ktsI_B = \frac{k}{\sqrt{t_s}}

No âmbito da norma IEEE Std 80, Charles Dalziel estabeleceu empiricamente que a energia tolerável pelo corpo humano antes de atingir o limiar de 0,5% de probabilidade de fibrilação ventricular é regida por uma constante de energia de choque SbS_b. Para uma população com peso corporal de aproximadamente 50 kg, a constante k50=0,116k50 = 0,116, enquanto para 70 kg, k70=0,157k70 = 0,157. Portanto, a corrente máxima admissível pelo corpo IBI_B para uma duração de choque tst_s (onde 0,03 \le t_s \le 3,0 s ) é definida analytically por:

IB,50kg=0,116ts[A],IB,70kg=0,157ts[A]I_{B,50 kg } = \frac{0,116}{\sqrt{t_s}} \quad [ A ], \qquad I_{B,70 kg } = \frac{0,157}{\sqrt{t_s}} \quad [ A ]

Por outro lado, a Comissão Eletrotécnica Internacional em seu relatório técnico IEC/TS 60479-1 adota uma abordagem probabilística não linear baseada nos estudos de Biegelmeier. Em vez de assumir uma relação puramente inversamente proporcional à raiz quadrada do tempo para toda a janela temporal, a IEC 60479-1 define curvas de tempo/corrente (zonas AC-1 a AC-4.3). A impedância do corpo humano ZBZ_B na IEC não é uma resistência fixa, mas uma rede resistivo-capacitiva dependente da tensão de toque (UTU_T), da trajetória da corrente e do estado de umidade da pele:

ZB(UT)=Rint+Rp1(UT)1+jωRp1Cp1+Rp2(UT)1+jωRp2Cp2Z_B(U_T) = Rint + \frac{Rp1(U_T)}{1 + j\omega Rp1 Cp1} + \frac{Rp2(U_T)}{1 + j\omega Rp2 Cp2}

Enquanto a IEEE Std 80 simplifica conservadoramente a resistência do corpo humano a um valor constante puramente resistivo R_B = 1000\,\Omega (representativo do trajeto mão-pé ou mão-mão com contato firme), as normas IEC 61936-1 e IEC 60479-1 aplicam curvas dinâmicas onde a impedância efetiva decai drasticamente à medida que a tensão de toque supera 200 V, alcançando valores assintóticos próximos à resistência interna do corpo (Rint \approx 500 a 750\,\Omega).

Arcabouço Matemático e Formulação Rigorosa segundo a IEEE Std 80

Circuitos Equivalentes de Thevenin e Limites de Tolerabilidade

Para determinar as tensões admissíveis de passo (VstepVstep) e de toque (VtouchVtouch), recorre-se a um circuito equivalente de Thevenin visto a partir dos pontos de contato do corpo com a superfície do terreno e as estruturas aterradas. A resistência equivalente do solo vista pelos pés do indivíduo é modelada considerando cada pé como um disco condutor plano de raio b=0,08mb = 0,08 m apoiado sobre a superfície de um semiespaço homogêneo de resistividade aparente \rho_s.

A resistência própria de um pé em contato com o solo é dada pela formulação de Maxwell:

Rfoot=ρs4bρs4(0,08)=3,125ρs3ρsRfoot = \frac{\rho_s}{4b} \approx \frac{\rho_s}{4(0,08)} = 3,125 \rho_s \approx 3\rho_s

Na condição de tensão de toque, ambos os pés encontram-se em paralelo a uma distância praticamente idêntica do ponto de falta, de modo que a resistência equivalente é R_{eq,touch} = \frac{Rfoot}{2} = 1,5\rho_s. Para a condição de tensão de passo, os pés encontram-se em série separados por uma distância normalizada de 1 metro, resultando em R_{eq,step} = 2 Rfoot = 6\rho_s.

Ao incorporar uma camada superficial de material de alta resistividade (tipicamente brita ou pedra britada de espessura hsh_s e resistividade \rho_s) sobre um terreno natural de resistividade \rho, deve-se introduzir um fator de redução ou derating C_s(h_s, K), o qual compensa o efeito de reflexão provocado pela descontinuidade dielétrica:

Cs=1+16bπn=1Kn1+(2nhsb)210,09(1ρρs)2hs+0,09C_s = 1 + \frac{16 b}{\pi} \sum_{n=1}^{\infty} \frac{K^n}{\sqrt{1 + \left(2 n \frac{h_s}{b}\right)^2}} \approx 1 - \frac{0,09 \left(1 - \frac{\rho}{\rho_s}\right)}{2 h_s + 0,09}

Onde o coeficiente de reflexão KK entre o terreno natural e a camada de brita superficial é expresso como:

K=ρρsρ+ρsK = \frac{\rho - \rho_s}{\rho + \rho_s}

Integrando essas deduções eletrodinâmicas, os potenciais máximos toleráveis calculados sob a IEEE Std 80 para um operador de 50 kg e 70 kg resultam nas seguintes equações fundamentais:

Vtouch,50=(1000+1,5Csρs)0,116ts,Vtouch,70=(1000+1,5Csρs)0,157tsV_{touch,50} = \left(1000 + 1,5 C_s \rho_s\right) \frac{0,116}{\sqrt{t_s}}, \qquad V_{touch,70} = \left(1000 + 1,5 C_s \rho_s\right) \frac{0,157}{\sqrt{t_s}}
Vstep,50=(1000+6,0Csρs)0,116ts,Vstep,70=(1000+6,0Csρs)0,157tsV_{step,50} = \left(1000 + 6,0 C_s \rho_s\right) \frac{0,116}{\sqrt{t_s}}, \qquad V_{step,70} = \left(1000 + 6,0 C_s \rho_s\right) \frac{0,157}{\sqrt{t_s}}

Elevação do Potencial de Terra (GPR) e Parâmetros de Projeto da Malha

A corrente de projeto que efetivamente se dissipa através da malha para o terreno infinito (IGI_G) não equivale à corrente de curto-circuito simétrica trifásica ou monofásica total nos terminais (IfI_f), mas é afetada pelo fator de divisão da corrente de falta (SfS_f) e pelo fator de decremento (DfD_f) devido à assimétria da componente de corrente contínua:

IG=SfDfIfI_G = S_f \cdot D_f \cdot I_f

O fator de decremento DfD_f quantifica o conteúdo aperiódico da corrente de curto-circuito em função da relação reatância/resistência (X/RX/R) da impedância de Thevenin no ponto de falta e da duração do tempo de eliminação tft_f:

Df=1+Tatf(1e2tf/Ta),Ta=XωR=X2πfRD_f = \sqrt{1 + \frac{T_a}{t_f} \left(1 - e^{-2 t_f / T_a}\right)}, \qquad T_a = \frac{X}{\omega R} = \frac{X}{2\pi f R}

A máxima elevação de potencial do sistema em relação a um ponto de referência remoto é denominada Ground Potential Rise (GPRGPR) e é calculada por:

GPR=IGRgGPR = I_G \cdot R_g

A resistência do sistema de aterramento (RgR_g) para malhas complexas interconectadas com hastes ou eletrodos verticais profundos é determinada com alta precisão analítica por meio da fórmula generalizada de Sverak:

Rg=ρ[1LT+120A(1+11+h20/A)]R_g = \rho \left[ \frac{1}{L_T} + \frac{1}{\sqrt{20 A}} \left( 1 + \frac{1}{1 + h \sqrt{20 / A}} \right) \right]

Onde LTL_T é o comprimento total dos condutores enterrados (malha horizontal mais eletrodos verticais), AA é a área ocupada pela malha em metros quadrados, e hh é a profundidade de sepultamento dos condutores da rede horizontal.

Tensão de Malha e Tensão de Passo Reais na Subestação

A tensão de malha real (EmE_m), que representa a pior condição de tensão de toque no centro do quadrante exterior mais crítico da malha, é formulada rigorosamente através do produto da resistividade aparente do terreno (\rho), do fator geométrico de espaçamento KmK_m, do fator de irregularidade corretiva KiK_i e da densidade linear de corrente injetada:

Em=ρIGKmKiLME_m = \frac{\rho \cdot I_G \cdot K_m \cdot K_i}{L_M}

Os fatores adimensionais KmK_m e KiK_i sintetizam a topologia da rede e são definidos analiticamente como:

Km=12π[ln(D216hd+(D+2h)28Ddh4d)+KiiKhln(8π(2n1))]K_m = \frac{1}{2\pi} \left[ \ln\left( \frac{D^2}{16 h d} + \frac{(D + 2h)^2}{8 D d} - \frac{h}{4d} \right) + \frac{Kii}{K_h} \ln\left( \frac{8}{\pi (2n - 1)} \right) \right]
Ki=0,644+0,148nK_i = 0,644 + 0,148 \cdot n

Onde DD é o espaçamento entre condutores paralelos, dd é o diâmetro do condutor da malha, hh é a profundidade de enterramento, nn é o fator geométrico equivalente derivado das dimensões e formato perimétrico (n = n_a \cdot n_b \cdot n_c \cdot n_d), K_h = \sqrt{1 + h/h_0} (com h0=1mh_0 = 1 m), e Kii=1Kii = 1 para malhas com eletrodos perimetrais ou Kii = 1/(2n)^{2/n} para malhas sem hastes.

O comprimento efetivo ponderado do condutor para o cálculo da tensão de malha (LML_M) incorpora a contribuição dos eletrodos verticais de aterramento (LrL_r):

LM=Lc+[1,55+1,22(LrLx2+Ly2)]LRL_M = L_c + \left[ 1,55 + 1,22 \left( \frac{L_r}{\sqrt{L_x^2 + L_y^2}} \right) \right] L_R

Analogamente, a tensão de passo calculada (EsE_s) no perímetro externo da subestação é avaliada com o coeficiente geométrico de passo KsK_s:

Es=ρIGKsKiLSE_s = \frac{\rho \cdot I_G \cdot K_s \cdot K_i}{L_S}
Ks=1π[12h+1D+h+1D(10,5n2)]K_s = \frac{1}{\pi} \left[ \frac{1}{2h} + \frac{1}{D + h} + \frac{1}{D} \left( 1 - 0,5^{n-2} \right) \right]
LS=0,75Lc+0,85LRL_S = 0,75 \cdot L_c + 0,85 \cdot L_R

Comportamento Transitório e Impacto de Descargas Atmosféricas segundo a IEC 62305

Impedância de Impulso versus Resistência em Frequência Industrial

A resposta de uma malha de aterramento ante a injeção de correntes do tipo raio reguladas pela série de normas IEC 62305 difere radicalmente do comportamento a 50/60 Hz. As descargas atmosféricas diretas ou induzidas apresentam frentes de onda ultrarrápidas normalizadas (forma de onda 10/350 \mu s para o primeiro impacto de raio segundo a IEC 62305-1, e 8/20 \mu s para impactos subsequentes), com conteúdos frequency estendendo-se na faixa de centenas de quilohertz a megahertz.

A estas frequências, a reatância indutiva distribuída longitudinal (LL') dos condutores de cobre ou aço galvanizado passa a ser o fator dominante sobre a condutância transversal do terreno (GG'). Como consequência, a corrente não tem tempo físico para se distribuir por toda a geometria da malha, concentrando-se nas imediações do ponto de impacto.

Zimp(t)=u(t)i(t)RgZimp(t) = \frac{u(t)}{i(t)} \neq R_g

O comprimento efetivo do condutor enterrado lel_e, além do qual qualquer extensão adicional do eletrodo não reduz a impedância de impulso aparente para um tempo de frente T1T_1, é expresso analiticamente segundo a IEC 62305-3 como:

le=ρT1μ01,6ρT1l_e = \sqrt{\frac{\rho \cdot T_1}{\mu_0}} \approx 1,6 \cdot \sqrt{\rho \cdot T_1}

Onde \mu_0 = 4\pi \times 10^{-7} H/m é a permeabilidade magnética do váCuO e T1T_1 é o tempo de frente do impulso em microssegundos (\mu s ).

Tensões de Passo e Toque Induzidas por Descargas Atmosféricas

Sob a ação de uma corrente impulsiva de raio i(t), o gradiente de potencial radial na superfície do solo gera tensões transitórias extremas que podem provocar descargas disruptivas no ar ao longo das extremidades inferiores do corpo humano (arco superficial ou flashover). A tensão de toque impulsiva transitória (Utouch,impU_{touch,imp}) e a tensão de passo transitória (Ustep,impU_{step,imp}) de acordo com o Anexo E da IEC 62305-3 acoplam-se por meio de componentes indutivas e resistivas puras:

Utouch,imp(t)=Rimpi(t)+Lgriddi(t)dtU_{touch,imp}(t) = Rimp \cdot i(t) + Lgrid \cdot \frac{di(t)}{dt}
U_{step,imp}(r) = \frac{\rhoeff \cdot i(t)}{2\pi} \left( \frac{1}{r} - \frac{1}{r + \Delta r} \right)

No caso de impactos de raio sobre a estrutura de blindagem (para-raios tipo Franklin, cabos para-raios ou pórticos), as elevadas derivadas temporais (di/dt > 100 kA/ \mu s ) geram potenciais indutivos parasitários nas descidas metálicas que excedem a rigidez dielétrica do ar, exigindo distâncias de separação seguras (ss) calculadas de acordo com a formulação:

skikckmls \ge k_i \cdot \frac{k_c}{k_m} \cdot l

Onde kik_i é o coeficiente que depende do nível de proteção contra raios (LPL I a IV), kck_c é o fator de partição geométrica da corrente nas descidas, kmk_m é o coeficiente do material isolante (ar ou sólido), e ll é a distância linear ao longo da descida até o ponto de equipotencialização mais próximo.

Matriz Comparativa Normativa: IEEE Std 80 vs. IEC 60479-1 / IEC 61936-1 / IEC 62305-3

Parâmetro / Critério de Projeto IEEE Std 80 (2013) IEC 60479-1 / IEC 61936-1 (2021) IEC 62305-3 (Proteção contra Raios)
Modelo de Peso Corporal Discreto: 50 kg (k=0,116k=0,116) ou 70 kg (k=0,157k=0,157) Populacional estatístico (Percentil 5%, 50%, 95%) Não aplicável (Risco avaliado por choque impulsivo)
Impedância do Corpo (RBR_B / ZBZ_B) Constante: R_B = 1000\,\Omega (resistivo puro) Não linear: Z_B = f(U_T), 500 a 1000\,\Omega R_B \approx 0\,\Omega (Arco elétrico em alta frequência)
Regime Temporal Crítico Frequência industrial (50/60 Hz), 0,03 \le t_s \le 3,0 s Frequência industrial (50/60 Hz), 0,01 \le t_s \le 10,0 s Transitório impulsivo (10/350\,\mu s , 8/20\,\mu s )
Resistência da Camada Superficial Derating analítico formal C_s(h_s, K) Fator de camada superficial similar ou empírico Isolamento de controle: asfalto 5cm5 cm (\ge 100 k \Omega\cdot m )
Limite de Tensão de Toque Vtouch = (1000 + 1,5 C_s \rho_s) \frac{k}{\sqrt{t_s}} Definido por curva limite U_{v,tol} = f(t_f, Z_B, Rfoot) Utouch \le 100 kV (com equipotencialização)
Fator de Divisão de Falta (SfS_f) Cálculo explícito com cabos para-raios e neutros Fator de redução rr (função de blindagens e cabos guarda) Fator de divisão kck_c baseado na topologia das descidas
Comportamento do Terreno Estratificação de resistividade multicamada a 50 Hz Modelos bicamada e semiespaço homogêneo Efeito de ionização do solo (E_0 \approx 300 kV/m )

Análise Forense de Falhas e Fenômenos de Transferência de Potencial

Mecanismos de Propagação de Potenciais Transferidos

Um dos modos de falha mais devastadores na engenharia forense de subestaciones corresponde à exportação ou importação inadvertida da Elevação do Potencial de Terra (GPRGPR) para fora dos limites perimétricos da instalação. Quando ocorre um curto-circuito à terra dentro da subestação, toda a malha de terra eleva-se a um potencial V=GPRV = GPR em relação ao potencial zero de referência em terra remota.

Vtrans=GPRΔVdropIGRgVtrans = GPR - \Delta Vdrop \approx I_G \cdot R_g

Os principais vetores físicos de acoplamento galvânico de potenciais transferidos compreendem:

  • Condutores Neutros de Redes de Média e Baixa Tensão: Se o neutro do enrolamento secundário do transformador de serviços auxiliares ou uma linha de distribuição rural estiver rigidamente aterrado na malha interna e se estender para fora da subestação com múltiplos aterramentos, o GPRGPR completo é aplicado diretamente sobre os isoladores das redes dos usuários finais, resultando em explosões de medidores residenciais, perfuração de eletrodomésticos e risco de eletrocussão em massa.
  • Blindagens Metálicas de Cabos de Potência e Telecomunicações: Blindagens de cabos subterrâneos de MT aterradas em ambas as extremidades transferem a corrente de retorno do solo para terminais remotos, provocando correntes circulantes de centenas de ampères que fundem a capa externa de polietileno (falha térmica por I2tI^2 t) e induzem sobretensões que perfuram as interfaces de instrumentação digital.
  • Tubulações de Condução Hidráulica ou Gás e Trilhos de Ferrovia: Elementos metálicos contínuos enterrados sem juntas isolantes dielétricas atuam como eletrodos de extensão infinita, expondo operadores em plantas de bombeamento ou passagens a nível distantes à totalidade do GPRGPR, onde a tensão de toque transferida V_{touch,trans} \approx GPR carece do efeito amortecedor da brita da subestação.

Comportamento em Subestações Isoladas a Gás (GIS)

Nas subestações isoladas a exofluoreto de enxofre (SF6SF _6), a operação de seccionadoras ou a ocorrência de descargas disruptivas gera Transitórios Muito Rápidos (VFT, Very Fast Transients) com tempos de subida da ordem de nanossegundos (tr<10nst_r < 10 ns). Essas frentes de onda não conseguem circular para a terra através das conexões convencionais de aterramento devido à reatância parasitária indutiva das barras de conexão planas (L \approx 1\,\mu H/m ).

Em consequência, origina-se o fenômeno de Elevação Transitória de Potencial do Invólucro (TEV, Transient Enclosure Voltage), onde a carcaça externa de alumínio da GIS experimenta oscilações de potencial de até várias dezenas de quilovolts em relação às estruturas metálicas de suporte e pisos circundantes, gerando microarcos secundários que destroem os cabos de controle do sistema de automação e proteção digital (SAS) e causam descargas por contato aos inspetores da planta.

Metodologia de Estratificação do Solo e Otimización de Malha

Inversão de Curvas de Sondagem Elétrica Vertical (SEV)

A precisão dos cálculos de VmV_m, VsV_s e RgR_g está estritamente vinculada à formulação do modelo geoelétrico do subsolo. A prospecção por meio do método tetrapontual de Wenner consiste em injetar corrente contínua ou de baixíssima frequência através de eletrodos externos (C1,C2C_1, C_2) e medir a queda de potencial nos eletrodos internos (P1,P2P_1, P_2), dispostos simetricamente com um espaçamento aa:

ρa(a)=2πaΔVI\rho_a(a) = 2\pi a \frac{\Delta V}{I}

Dado que o solo raramente é homogêneo, a resistividade aparente \rho_a varia em função da profundidade de penetração das linhas de corrente, que é proporcional ao espaçamento aa. Para modelar uma estrutura bicamada (camada superior de resistividade \rho_1 e espessura hh, sobre um semiespaço profundo de resistividade \rho_2), resolve-se analiticamente a equação integral de Stefanescu por meio da transformada de Hankel de ordem zero:

ρa(a)=ρ1[1+4n=1(Kn1+(2nha)2Kn4+(2nha)2)]\rho_a(a) = \rho_1 \left[ 1 + 4 \sum_{n=1}^{\infty} \left( \frac{K^n}{\sqrt{1 + \left(\frac{2 n h}{a}\right)^2}} - \frac{K^n}{\sqrt{4 + \left(\frac{2 n h}{a}\right)^2}} \right) \right]

Onde o fator de reflexão geoelétrico é:

K=ρ2ρ1ρ2+ρ1K = \frac{\rho_2 - \rho_1}{\rho_2 + \rho_1}

O algoritmo de otimização numérica (mínimos quadrados amortecidos não lineares do tipo Levenberg-Marquardt) busca minimizar a função objetivo de erro quadrático médio RMS:

varepsilonRMS=1Ni=1N(ρa,meas(ai)ρa,calc(ai)ρa,meas(ai))2×100%varepsilonRMS = \sqrt{\frac{1}{N} \sum_{i=1}^{N} \left( \frac{\rho_{a,meas}(a_i) - \rho_{a,calc}(a_i)}{\rho_{a,meas}(a_i)} \right)^2} \times 100\%

Topologias Avançadas de Mitigação de Gradientes Superficiais

Quando a resistividade superficial é desfavorável (\rho_1 \gg \rho_2) ou a corrente de falta IGI_G é massiva, o simples acréscimo de condutores uniformes na rede apresenta rendimentos decrescentes. As técnicas avançadas de otimização eletromagnética incluem:

  • Malha de Densificação Gradual (Compressão Perimétrica): Redução exponencial do espaçamento entre condutores à medida que se aproxima do limite externo da subestação. Isso uniformiza a densidade de corrente expelida por metro linear, aplainando os picos periféricos de potencial que disparam a tensão de toque externa e de passo.
  • Eletrodos Profissionais Verticais no Perímetro e Vértices: A inserção de hastes profundas (15 a 30 metros) perfuradas nos cantos da malha permite drenar corrente preferencialmente para os estratos inferiores mais condutivos (\rho_2), colapsando o GPRGPR geral e reduzindo substancialmente o gradiente superficial nas bordas críticas.
  • Camada Dielétrica de Alto Desempenho: Aplicação de camadas de brita lavada de granito (\rho_s \ge 3000\,\Omega\cdot m ) ou revestimento asfáltico (\rho_s \ge 10000\,\Omega\cdot m ) com espessura hs=0,10a0,15mh_s = 0,10 a 0,15 m. Deve-se verificar estritamente se o fator de redução CsC_s mantém sua eficácia ao longo da vida útil da planta, mitigando a contaminação por sedimentos finos e mato por meio de geotêxteis impermeáveis.

Implementação e Modelagem Computacional na Vexten Suite

Integração de Curtos-Circuitos e Avaliação Paramétrica

No fluxo de trabalho de projeto de engenharia dentro da plataforma integral Vexten Suite, o dimensionamento dos aterramentos é sincronizado nativamente com os módulos de cálculo de curto-circuito trifásico e assimétrico sob as normas IEC 60909 e IEEE 141 (Red Book). A determinação exata da componente de corrente de sequência zero (3I03I_0) e da impedância equivalente de Thevenin no nó da subestação alimenta diretamente o motor de otimização de terras.

A formulação matricial a seguir é implementada na Vexten Suite para a resolução do acoplamento indutivo e galvânico de múltiplas linhas aéreas com cabo para-raios (OPGW/aço) conectadas à subestação, permitindo determinar com total rigor o fator de divisão real SfS_f:

[Ig1Ig2Ign]=[Z11Z12Z1nZ21Z22Z2nZn1Zn2Znn]1[VGPREind,1VGPREind,2VGPREind,n]\begin{bmatrix} \mathbf{I}_{g1} \\ \mathbf{I}_{g2} \\ \vdots \\ \mathbf{I}_{gn} \end{bmatrix} = \begin{bmatrix} Z11 & Z12 & \cdots & Z1n \\ Z21 & Z22 & \cdots & Z2n \\ \vdots & \vdots & \ddots & \vdots \\ Zn1 & Zn2 & \cdots & Znn \end{bmatrix}^{-1} \begin{bmatrix} \mathbf{V}_{GPR} - \mathbf{E}_{ind,1} \\ \mathbf{V}_{GPR} - \mathbf{E}_{ind,2} \\ \vdots \\ \mathbf{V}_{GPR} - \mathbf{E}_{ind,n} \end{bmatrix}

Mapeamento 3D de Gradientes Superficiais e Verificación Normativa

O motor eletromagnético da Vexten Suite subdivide a malha de terra e os eletrodos verticais em segmentos cilíndricos discretos, aplicando o Método dos Momentos (MoM) para calcular a densidade de carga longitudinal qjq_j e a corrente de fuga radial em cada segmento. A distribuição de potencial em qualquer coordenada superficial (x, y, 0) é avaliada analiticamente por meio de integração numérica:

Φ(x,y,0)=14πσj=1MLjIj()(xx)2+(yy)2+(z)2d\Phi(x, y, 0) = \frac{1}{4\pi \sigma} \sum_{j=1}^{M} \int_{L_j} \frac{I_j(\ell')}{\sqrt{(x - x')^2 + (y - y')^2 + (z')^2}} \, d\ell'

Uma vez calculado o mapa tridimensional de potenciais, a Vexten Suite processa automaticamente a topologia da planta sobrepondo os equipamentos de manobra, transformadores de potência, cercas perimétricas e caminhos de operação. O sistema gera vetores espaciais de 1 metro em todas as direções a partir das massas condutoras acessíveis para contrastar ponto a ponto:

\max\left( \Phiestrutura - \Phi(x_{pé}, y_{pé}, 0) \right) \le V_{touch,adm}(t_s)
max(Φ(x+1,y,0)Φ(x,y,0))Vstep,adm(ts)\max\left( \Phi(x+1, y, 0) - \Phi(x, y, 0) \right) \le V_{step,adm}(t_s)

Qualquer transgressão aos limites normativos da IEEE Std 80 ou da IEC 61936-1 é destacada no ambiente gráfico por meio de gradientes de cor de alta resolução (heatmaps), permitindo que o engenheiro aplique de forma interativa hastes auxiliares, equalizadores de gradiente perimétricos ou modificações na espessura da camada de brita até certificar a conformidade total da instalação frente aos máximos requisitos de segurança humana e robustez operacional.