Investigação dos Efeitos de Descargas Atmosféricas em Centrais Solares Conectadas a Redes de Transmissão

1. Introdução

A rápida integração de centrais solares fotovoltaicas (FV) de grande escala em redes de transmissão de alta tensão introduz novas vulnerabilidades a perturbações da rede, particularmente a descargas atmosféricas. Este artigo investiga a propagação de sobretensões induzidas por raios desde as linhas de transmissão até às centrais solares conectadas, uma questão crítica dada a sobreposição geográfica de regiões de alta irradiação solar e alta atividade de descargas atmosféricas. O estudo emprega simulações com o Programa de Transitórios Eletromagnéticos (EMTP) para modelar o sistema e avalia a eficácia dos pára-raios como uma medida de proteção primária.

Ideias-Chave

Descargas atmosféricas em linhas de transmissão podem induzir sobretensões severas no Ponto de Acoplamento Comum (PCC) das centrais solares.
A vulnerabilidade é agravada pelos longos percursos de cabos e pelos eletrónica de potência sensível (inversores) dentro das centrais FV.
As estratégias de proteção padrão concebidas para geração tradicional podem ser inadequadas para recursos distribuídos baseados em inversores, como a solar.

2. Metodologia & Modelagem do Sistema

A investigação baseia-se numa metodologia orientada por simulação, utilizando o software padrão da indústria EMTP-RV para uma modelagem precisa de transitórios eletromagnéticos.

2.1 Estrutura de Simulação EMTP

Todo o sistema — compreendendo a linha de transmissão, a rede de recolha da central solar, os transformadores e os dispositivos de proteção contra sobretensões — foi modelado no EMTP. Isto permite a análise no domínio do tempo de surtos de frente rápida com resolução de nanossegundos a microssegundos.

2.2 Modelo de Descarga Atmosférica & Central Solar

A descarga atmosférica é modelada utilizando a função de fonte de corrente de Heidler, um padrão para representar a corrente do canal: $i(t) = \frac{I_0}{\eta} \frac{(t/\tau_1)^n}{1+(t/\tau_1)^n} e^{-t/\tau_2}$. Os parâmetros $I_0$ (corrente de pico), $\tau_1$ (tempo de frente) e $\tau_2$ (tempo de cauda) foram variados. A central solar foi modelada como um circuito equivalente agregado, incluindo cabos DC, inversores e transformadores elevadores.

2.3 Configuração do Pára-Raios

Pára-raios de Varistor de Óxido Metálico (MOV) foram modelados em locais-chave: na torre da linha de transmissão perto do ponto de impacto e no ponto principal de ligação AC da central solar. A sua característica não linear V-I é dada por $i = k \cdot V^{\alpha}$, onde $k$ e $\alpha$ são constantes do dispositivo.

3. Cenários & Parâmetros de Simulação

3.1 Variação dos Parâmetros da Descarga

As simulações abrangeram uma gama de parâmetros realistas de descargas atmosféricas:

Corrente de Pico (I_p): 10 kA a 100 kA (representando tanto descargas negativas como positivas).
Tempo de Frente (t_f): 1 µs a 10 µs.
Tempo de Cauda (t_t): 20 µs a 200 µs.

Esta matriz permite avaliar o impacto tanto de descargas rápidas e de alta corrente como de eventos mais lentos e de maior duração.

3.2 Cenários de Distância do Impacto

Descargas atmosféricas foram simuladas a distâncias variáveis (ex.: 0,5 km, 1 km, 2 km) do ponto de ligação à rede da central solar ao longo da linha de transmissão. Foram considerados tanto impactos diretos no condutor de fase (falha de blindagem) como retroflashes devido a impactos na torre.

4. Resultados & Análise

4.1 Análise da Magnitude da Sobretensão

A métrica principal foi a magnitude da sobretensão transitória na barra AC da central solar. Sem pára-raios, as sobretensões frequentemente excediam 3,0 p.u. (por unidade) da tensão nominal do sistema para impactos dentro de 1 km, representando um risco severo para o isolamento dos inversores. A forma de onda da sobretensão é uma sobreposição do surto incidente e das reflexões dentro da rede interna de cablagem da central.

Descrição do Gráfico (Imaginado): Um gráfico de linhas mostraria a sobretensão (p.u.) no eixo Y versus a distância do impacto do raio (km) no eixo X. Duas linhas seriam traçadas: uma (vermelha, em declive acentuado) para o cenário sem pára-raios, mostrando tensões elevadas a curtas distâncias; e outra (azul, mais plana) para o cenário com pára-raios, mostrando tensões significativamente limitadas em todas as distâncias.

4.2 Análise Espectral de Fourier & Hilbert

Para além da magnitude no domínio do tempo, o estudo realizou análise espectral.

Transformada de Fourier: Revelou os componentes de frequência dominantes da sobretensão. Sem pára-raios, a energia concentrava-se em bandas de alta frequência (100 kHz - 1 MHz), que são particularmente danosas para dispositivos semicondutores. Com pára-raios, o espetro deslocou-se para frequências mais baixas.
Transformada de Hilbert-Huang (HHT) / Espetro Marginal: Esta análise tempo-frequência forneceu uma perceção de como a distribuição de energia evoluiu durante o evento transitório, mostrando a natureza não estacionária do surto e o efeito dinâmico de limitação do pára-raios.

4.3 Desempenho do Pára-Raios

Os pára-raios demonstraram elevada eficácia, limitando tipicamente as sobretensões para valores abaixo de 1,8 p.u., um nível geralmente dentro da capacidade de suporte dos inversores FV modernos (tipicamente classificados para 2,0-2,5 p.u. por curtos períodos). O requisito de absorção de energia para os pára-raios foi quantificado, o que é crítico para uma dimensionamento adequado.

Redução da Sobretensão de Pico

> 40%

Redução média com pára-raios instalados

Distância Crítica de Impacto

< 1 km

Impactos dentro desta gama causam o maior risco

5. Detalhes Técnicos & Formulação Matemática

O núcleo do modelo EMTP baseia-se na resolução das equações do telégrafo para a linha de transmissão, acopladas com modelos de componentes não lineares:

Linha de Transmissão (Modelo Dependente da Frequência): Resolvido usando o método das características: $\frac{\partial v}{\partial x} + L' \frac{\partial i}{\partial t} + R' i = 0$ e $\frac{\partial i}{\partial x} + C' \frac{\partial v}{\partial t} + G' v = 0$.
Modelo de Pára-Raios (MOV): A característica não linear por segmentos é frequentemente implementada usando o modelo $\alpha$-$k$ ou o modelo mais dinâmico de Pinceti-Giannettoni para rastreamento de energia.
Impedância do Inversor: A impedância de alta frequência do inversor FV, crucial para a divisão do surto, foi modelada como um circuito RLC paralelo baseado em projetos típicos de filtros.

6. Estrutura de Análise: Estudo de Caso

Cenário: Uma central solar de 100 MW conectada a uma linha de transmissão de 230 kV através de um transformador elevador 230/33 kV. Uma descarga atmosférica com I_p = 50 kA, t_f = 2 µs atinge uma torre a 0,8 km de distância, causando um retroflash.

Aplicação da Estrutura:

Configuração do Modelo: Construir o modelo EMTP com constantes de linha detalhadas, resistência de aterramento da torre (50 Ω) e impedância interna da central.
Execução de Base (Sem Proteção): Simular. Registar a sobretensão no PCC (~3,5 p.u., frequência dominante de 0,5 MHz).
Execução de Mitigação (Com Pára-Raios): Colocar pára-raios na torre atingida e no PCC. Re-simular. Registar a tensão limitada (~1,7 p.u., frequência dominante < 100 kHz).
Cálculo de Energia: Calcular a energia absorvida pelo pára-raios do PCC usando $W = \int v(t) \cdot i_{arrester}(t) dt$ para verificar se a sua classificação não é excedida.
Análise de Sensibilidade: Variar a resistência de aterramento e a impedância da central para ver o impacto na sobretensão.

Esta abordagem estruturada isola variáveis e quantifica os benefícios da proteção.

7. Perspetiva de Aplicação & Direções Futuras

Os resultados têm aplicações diretas no projeto e nos códigos de rede para grandes instalações solares:

Códigos de Rede Reforçados: Operadores de Sistemas de Transmissão (OST) como a PJM ou a ENTSO-E poderiam exigir estudos específicos de proteção contra sobretensões e especificações de pára-raios para centrais FV conectadas à rede em áreas de alta propensão a descargas atmosféricas (KERA).
Proteção Inteligente contra Surtos: Sistemas futuros poderiam integrar pára-raios habilitados para IoT que monitorizam o seu próprio estado e absorção de energia, comunicando com o SCADA da central para manutenção preditiva.
Esquemas de Proteção Híbridos: Combinar pára-raios MOV tradicionais com tecnologias emergentes como limitadores de corrente de falta em série (SFCL) ou grampos ativos baseados em semicondutores de banda larga poderia oferecer proteção superior com resposta mais rápida.
Integração de Gémeo Digital: Os modelos EMTP desenvolvidos nesta investigação podem formar a base de um gémeo digital para centrais solares operacionais, permitindo uma avaliação de risco em tempo real durante tempestades usando dados de redes de deteção de descargas atmosféricas (ex.: da Vaisala GLD360 ou Earth Networks).

8. Referências

Grebovic, S., Aksamovic, A., Filipovic, B., & Konjicija, S. (2025). Investigation of Lightning Effects on Solar Power Plants Connected to Transmission Networks. Artigo submetido ao IPST2025.
IEEE Std 1410-2010: IEEE Guide for Improving the Lightning Performance of Electric Power Overhead Distribution Lines.
CIGRE WG C4.408. (2013). Lightning Protection of Large Wind Turbine Blades. (Fornece metodologia relevante para estruturas de energia renovável).
Martinez, J. A., & Walling, R. A. (2013). EMTP Modeling of Inverter-Based Resources for Power System Dynamic Studies. IEEE Transactions on Power Delivery.
Vaisala. (2023). Annual Lightning Report 2022. [Online]. Disponível: https://www.vaisala.com
Isola, G., et al. (2020). Advanced Surge Arrester Models for Fast Transient Simulations in EMTP. Electric Power Systems Research.

9. Perspetiva do Analista: Ideia Central & Crítica

Ideia Central

Este artigo identifica corretamente uma linha de falha crítica, mas frequentemente subestimada, na transição energética: o conflito inerente entre a localização ótima das renováveis e a resiliência da rede. Os autores identificam que as próprias regiões com maior rendimento solar (zonas de cinturão solar) estão frequentemente co-localizadas com altos níveis isoceráunicos (dias de trovoada por ano). Isto não é uma mera coincidência; é um dilema fundamental de localização. A investigação muda efetivamente a narrativa de ver as centrais solares como cargas passivas e benignas para reconhecê-las como nós ativos e vulneráveis que importam e amplificam transitórios provenientes da rede, ameaçando a sua própria eletrónica de potência dispendiosa — sendo os inversores o calcanhar de Aquiles.

Fluxo Lógico

A lógica do artigo é robusta e segue um caminho clássico de avaliação de risco de engenharia: Identificação do Perigo → Modelagem do Sistema → Simulação da Consequência → Avaliação da Mitigação. Começa com o perigo plausível (descarga atmosférica no corredor de transmissão), modela a sua propagação através da complexa rede RLC de linhas e cablagem da central (usando a ferramenta EMTP validada pela indústria), quantifica a consequência danosa (sobretensão excedendo o Nível de Isolamento Básico do inversor) e, finalmente, testa uma ferramenta de mitigação padrão (pára-raios). A inclusão da análise tanto de Fourier como da Transformada de Hilbert-Huang adiciona uma camada valiosa, indo além do simples pico de tensão para compreender a assinatura no domínio da frequência da ameaça, que é mais relevante para a durabilidade dos semicondutores.

Pontos Fortes & Falhas

Pontos Fortes: O rigor metodológico é louvável. Usar o EMTP, o padrão de ouro para estudos de transitórios, confere credibilidade imediata. A variação de parâmetros (corrente, distância) fornece uma análise de sensibilidade útil. O foco na análise espectral é um passo à frente de muitos estudos puramente no domínio do tempo.

Falhas Críticas & Oportunidades Perdidas:

Ponto Cego Económico: O estudo para na eficácia técnica. Uma omissão flagrante é uma análise custo-benefício. Qual é o CAPEX/OPEX da proteção contra sobretensões recomendada versus o risco de falha do inversor (que pode custar milhões e incorrer em meses de inatividade)? Sem isto, as recomendações carecem de força acionável para os desenvolvedores de centrais.
Modelagem Estática: A central solar é modelada como um agregado passivo. Na realidade, os inversores controlam ativamente a tensão e a frequência. Sob um surto rápido, os seus laços de controlo podem interagir de forma imprevisível com o transitório, potencialmente agravando ou mitigando o evento. Esta resposta dinâmica do inversor é ignorada, uma simplificação que limita a precisão no mundo real, como notado em estudos dinâmicos por Martinez & Walling.
Mentalidade de Ponto Único de Falha: A solução é centralizada (pára-raios no PCC). Negligencia o potencial para uma estratégia de defesa em profundidade distribuída: pára-raios coordenados nas caixas combinadoras DC, terminais AC dos inversores e terminais dos transformadores, o que é uma prática comum no projeto moderno de centrais para proteger toda a cadeia de conversão de energia.

Ideias Acionáveis

Para utilities, desenvolvedores e OEMs:

Tornar Obrigatórios Estudos de Transitórios Específicos do Local: Os acordos de ligação à rede para centrais FV >20 MW em áreas propensas a descargas atmosféricas devem exigir um estudo EMTP detalhado como este, e não apenas uma lista de verificação de conformidade padrão. Isto deve ser defendido junto de organismos como o IEEE PES.
Desenvolver Especificações de Pára-Raios "Adaptadas às Renováveis": As normas para pára-raios MOV (IEEE C62.11) são genéricas. Os fabricantes de inversores e os produtores de pára-raios devem colaborar para definir características V-I otimizadas e classificações de energia para as formas de onda e ciclos de trabalho únicos vistos em aplicações FV.
Integrar Dados de Descargas Atmosféricas no SCADA da Central: Usar dados em tempo real de serviços como os da Vaisala para implementar um modo operacional de trovoada. Quando uma célula de tempestade estiver dentro de 10 km, a central poderia temporariamente reduzir a produção ou isolar-se se viável, reduzindo a exposição ao risco — uma forma de resiliência operacional inspirada em conceitos de inteligência na periferia da rede.
Financiar Investigação em Grampeamento Ativo: A indústria deve investir em I&D para proteção usando dispositivos SiC/GaN que podem grampear ativamente tensões em microssegundos, oferecendo proteção mais rápida e precisa do que os MOVs passivos, semelhante à forma como os drivers avançados revolucionaram a eletrónica de potência noutros campos.

Em conclusão, este artigo é um alerta vital que acerta na definição do problema mas apenas o resolve parcialmente. O seu verdadeiro valor reside em fornecer a evidência de simulação fundamental necessária para impulsionar padrões de proteção mais holísticos, economicamente fundamentados e tecnologicamente avançados para a rede dominada pela energia solar do amanhã.