Incerteza total de AEP em projetos eólicos: decomposição e propagação até o P90

← Notas Técnicas

O P90 de AEP é o número que o lender usa para calcular se o projeto cobre a dívida. Não é o P50 — a estimativa central de energia — mas a estimativa que tem 90% de probabilidade de ser superada ao longo da vida do projeto. A diferença entre P50 e P90 é inteiramente determinada pela incerteza total da análise: quanto maior a incerteza, maior o gap, pior o DSCR, mais difícil o financiamento.

Entender como a incerteza total é construída — quais componentes a compõem, como se propagam e quais dominam o resultado — é o que separa um laudo de recurso que fecha financiamento de um que gera condições precedentes ou é devolvido para revisão. Esta nota apresenta a metodologia de decomposição e propagação de incertezas aplicada em laudos eólicos da Renewene, com os valores típicos observados em projetos onshore no Brasil.

Incerteza vs variabilidade: a distinção que importa

O primeiro equívoco comum na interpretação de P50/P90 é confundir incerteza (epistêmica) com variabilidade natural (aleatória). São origens distintas com tratamento estatístico diferente:

• Variabilidade interanual é a flutuação natural do recurso eólico ano a ano — causada por variações nos padrões climáticos, ENSO, NAO e teleconexões de larga escala. É real, existe no mundo físico e não pode ser eliminada com mais dados ou melhor metodologia. Contribui para a incerteza de longo prazo porque não sabemos em quais anos o projeto operará dentro de sua vida útil.
• Incerteza de medição é o desconhecimento sobre o valor verdadeiro de uma grandeza — causada por limitações do instrumento, método de calibração, ou representatividade da campanha. Pode ser reduzida com melhor instrumentação, campanhas mais longas e metodologia mais rigorosa.
• Incerteza de modelo é a diferença entre o modelo e a realidade — causada por simplificações nos modelos de wake, curva de potência, cisalhamento. Pode ser parcialmente reduzida com modelos mais sofisticados, mas nunca eliminada.

Na prática, as três contribuem para o P90 de forma combinada via propagação de incertezas (RSS). A variabilidade interanual é invariavelmente o componente dominante — o que significa que aumentar a precisão de medição além de certo ponto tem retorno decrescente se o P50-P90 gap for principalmente determinado pelo clima.

Propagação por soma quadrática (RSS)

O método padrão da indústria para combinar incertezas independentes é a soma quadrática das raízes (Root Sum of Squares — RSS), conforme descrito no Guia para a Expressão da Incerteza de Medição (GUM, BIPM/JCGM 100:2008). Para incertezas expressas como desvio padrão (1σ) de cada componente i:

σ_total = √(σ₁² + σ₂² + σ₃² + … + σₙ²)

A hipótese implícita é que as incertezas são independentes e não correlacionadas. Na prática, algumas correlações existem — por exemplo, um ano de medição meteorologicamente atípico afeta tanto a estimativa de recurso quanto a MCP — e o RSS pode tanto subestimar quanto superestimar a incerteza total dependendo do sinal da correlação. Laudos conservadores tratam as correlações de forma explícita; laudos básicos assumem independência total.

Uma vez calculada σ_total, a relação entre P-values é imediata pela distribuição normal:

P90 = P50 × (1 − 1.282 × σ_total) P75 = P50 × (1 − 0.674 × σ_total) P99 = P50 × (1 − 2.326 × σ_total)

O fator 1.282 é o percentil 90 da distribuição normal padrão (z₀.₉₀). A validade da distribuição normal para AEP — em vez de uma distribuição assimétrica — é justificada pelo teorema do limite central quando muitos componentes independentes se combinam, e é a convenção adotada universalmente pela indústria.

Os dez componentes de incerteza

A decomposição de incerteza de AEP pode ser organizada em quatro grupos. Os valores abaixo representam intervalos típicos observados em projetos onshore no Brasil (σ, 1 desvio padrão, em % da AEP estimada). Cada projeto terá valores específicos dependendo da qualidade da campanha e da complexidade do sítio:

Componente	Grupo	Típico (σ)	Faixa
Variabilidade interanual	Recurso	5.5%	4.0–7.0%
Wake losses (esteira)	Modelo	2.8%	1.5–4.0%
Correlação MCP / reanálise	Recurso	2.3%	1.5–3.5%
Modelo de fluxo (MMD)	Modelo	2.1%	1.0–3.5%
Curva de potência	Modelo	1.4%	0.8–2.0%
Disponibilidade de turbinas	Operação	1.2%	0.5–2.0%
Instrumento / LiDAR	Medição	1.1%	0.5–1.5%
Extrapolação vertical	Recurso	0.9%	0.3–1.5%
QA/QC · disponibilidade dados	Medição	0.7%	0.3–1.5%
Perdas elétricas	Operação	0.5%	0.2–1.0%
Total RSS		7.35%	5.0–10.0%

A variabilidade interanual domina amplamente — ela representa 56% da variância total (σ²). Wake losses e correlação MCP são os próximos maiores contribuintes. Instrumentação e QA/QC, quando bem executados, têm contribuição modesta — o que significa que investir em melhor medição reduz pouco o P90 se a variabilidade interanual não for endereçada com campanha longa.

Do P50 ao P90: a distribuição completa

Com σ_total = 7.35%, a distribuição de AEP em torno do P50 estimado determina todos os P-values relevantes para o financiamento. O P90 resulta em 90.6% do P50 — uma diferença de 9.4 pontos percentuais que, para um projeto de 300 GWh/ano de P50, representa 28.2 GWh de gap entre a estimativa central e o valor para fins de DSCR.

A área sombreada em vermelho representa 10% da probabilidade acumulada — a chance de a AEP real ficar abaixo do P90. Para o lender, o P90 é o nível de energia que o projeto cobre em 9 de cada 10 anos. Com σ = 7.35%, o P90 está 9.4 pontos percentuais abaixo do P50 — gap que define a margem do DSCR.

Sensibilidade: o que mais importa para reduzir o gap

Como cada componente contribui para o gap P50−P90 de forma proporcional à sua variância (σᵢ²), não ao seu desvio padrão (σᵢ), o componente com maior σ tem contribuição desproporcionalmente maior ao quadrado. Para um projeto hipotético com P50 = 300 GWh/ano:

P50 − P90 = 300 × 1.282 × 0.0735 = 28.2 GWh/ano

A contribuição de cada componente para esses 28.2 GWh é proporcional a σᵢ²/σ_total². O gráfico de tornado abaixo mostra essas contribuições em GWh — o formato que torna a conversa sobre incerteza concreta para investidores e lenders:

A variabilidade interanual consome 56% do budget de incerteza (em variância) — mais do que todos os outros componentes somados. Dobrar a precisão de instrumentação reduziria a contribuição da medição de ~3.1% para ~0.8% do gap, economizando ~0.7 GWh de P50-P90 num projeto de 300 GWh. A mesma energia e recursos investidos em estender a campanha de 1 para 3 anos reduziriam a variabilidade interanual de ~6.5% para ~4.0%, economizando ~5-8 GWh.

Como reduzir o gap P50−P90

A decomposição em FIG. 3 tem implicação direta para onde investir recursos na campanha:

1. Campanha longa (≥ 3 anos) é o investimento de maior retorno. A variabilidade interanual reduz na proporção de 1/√n anos de campanha. Uma campanha de 3 anos versus 1 ano reduz esta componente de ~6.5% para ~3.8%, economizando ~8 GWh de P90 num projeto de 300 GWh/ano. Nenhum outro investimento em qualidade tem retorno comparável.
2. Modelagem de wake com incerteza explícita é o segundo maior retorno. Wake losses é o segundo componente (14.5% do total). Modelos de alta fidelidade (LES, RANS) reduzem a incerteza do modelo de 2.8% para ~1.5-2.0%, economizando ~2-4 GWh de P90.
3. Escolha e validação da reanálise MCP (NT.02). A incerteza de MCP de 2.3% pode ser reduzida a ~1.5% com ensemble de duas reanálises e período concorrente de 3+ anos. Veja NT.02 para metodologia detalhada.
4. Instrumentação e QA/QC (NT.01, NT.03) têm retorno marginal baixo a partir de certo patamar — esses componentes somados respondem por menos de 4% do gap total. Melhorar QA/QC além do mínimo normativo não move o P90 de forma significativa.

Documentação exigida no laudo de recurso

Para que um P90 seja bankável, o laudo deve apresentar a decomposição de incerteza de forma explícita e rastreável — não apenas o número final de P90. Os lenders e ITAs verificam:

• Tabela de incertezas com todos os componentes, valores de σ em % da AEP e justificativa para cada valor adotado;
• Referência metodológica para cada componente (ex.: variabilidade interanual estimada pelo método de Carta et al. 2002 ou equivalente);
• Teste de sensibilidade mostrando como o P90 mudaria se cada componente variasse ±30% em torno do valor adotado;
• Demonstração de que os componentes são suficientemente independentes para aplicação de RSS — ou, onde correlações existem, correção explícita;
• Cálculo explícito de P50, P75, P90 e P99, com os fatores z correspondentes declarados.

Laudos que apresentam apenas "P90 = X GWh" sem a decomposição são invariavelmente devolvidos pelo ITA com solicitação de dados de suporte. A ausência da decomposição equivale, para fins de due diligence, a uma caixa preta: o número pode ser razoável, mas não pode ser verificado.