Geometric Time-Domain Identification of Three-Phase Load Equivalents from Terminal Measurements

Resumo Técnico: Identificação Geométrica no Domínio do Tempo de Equivalentes de Carga Trifásica

Definição do Problema
Os sistemas de potência modernos exigem parâmetros de carga em tempo real para proteção, monitoramento e controle adaptativo. Embora os métodos existentes estimem parâmetros estáticos ou dinâmicos agregados a partir de eventos ou quantidades no domínio dos fasores, eles frequentemente falham em recuperar equivalentes energéticos de nível de elemento diretamente de formas de onda terminais instantâneas. Um desafio específico surge em sistemas de três fios onde nenhuma tensão de neutro é medida. Nesses casos, as variáveis de ramo estão ocultas atrás de tensões linha-linha e correntes de linha restritas ($i_a + i_b + i_c = 0$), tornando difícil distinguir entre equivalentes delta e wye ou determinar se os dados medidos contêm informação suficiente para suportar um modelo físico escolhido.

Metodologia
O artigo propõe um método geométrico no domínio do tempo que interpreta as formas de onda de tensão e corrente medidas como trajetórias em espaços de sinais euclidianos. A metodologia central estende uma formulação de identificação geométrica de fase única publicada anteriormente para sistemas de três e quatro fios.

1. Formulação Geométrica: O método constrói um "spacor" (vetor espacial) a partir das variáveis medidas e de suas derivadas ou integrais necessárias (ex: $v, \dot{v}, \int v$). As Leis de Kirchhoff para Corrente e Tensão (KCL/KVL) confinam essas trajetórias a subespaços de menor dimensão dentro de um espaço euclidiano de maior dimensão.
2. Recuperação de Parâmetros: Os parâmetros equivalentes ($R, L, C$) são recuperados comparando a geometria da trajetória medida (o subespaço osculante) com o subespaço do modelo. No domínio simbólico, isso resulta em quocientes de forma fechada de coeficientes de lâmina do hipervolume osculante.
3. Implementação de Dados Amostrados: Para lidar com dados do mundo real, a teoria é implementada como um estimador de matriz por janela.
   • Janelamento: Os dados são processados em janelas locais finitas (tipicamente 0,5–2 ciclos).
   • Reconstrução de Derivadas: Derivadas e primitivas são computadas usando projeção harmônica, filtros FIR ou diferenças finitas regularizadas para evitar o ruído de diferenciação pontual de alta ordem.
   • Normalização: As equações matriciais são normalizadas para remover efeitos de escala dependentes da unidade.
4. Camadas de Diagnóstico: O método incorpora verificações explícitas antes de aceitar uma solução:
   • Identificabilidade: Testa o posto geométrico e o número de condição da trajetória local para garantir que a janela contenha informação suficiente (harmônicos ou transientes suficientes) para resolver a ordem do modelo escolhido.
   • Admissibilidade de Energia: Um portão de validação física que verifica se o equivalente identificado fecha o balanço de energia terminal instantâneo ($p_{PCC} = p_{dissipado} + \dot{W}_{armazenada}$).

Principais Contribuições
Em relação ao trabalho anterior de fase única, este artigo introduz:

• Formulação Trifásica Unificada: Um framework geométrico para conexões wye e delta em sistemas de três e quatro fios, abordando especificamente as restrições acopladas de sistemas de três fios sem um neutro medido.
• Análise de Observabilidade: Uma análise dependente da topologia identificando quais formas de ramo são diretamente recuperáveis. Prova-se que configurações delta em série e wye em paralelo são estruturalmente singulares em sistemas de três fios devido à falta de potenciais de fase absolutos ou correntes de ramo internas.
• Estimador de Matriz por Janela: Uma realização numérica da construção geométrica usando janelas finitas, derivadas por projeção harmônica e diagnósticos explícitos de posto/condicionamento.
• Portão de Admissibilidade de Energia: Um critério que rejeita equivalentes ajustados numericamente se eles falharem em satisfazer a lei de conservação de energia instantânea da topologia assumida.
• Triagem de Baixa Informação: Um mecanismo para reportar janelas com posto geométrico insuficiente (ex: excitação de frequência única pura para modelos de alta ordem) como evidência de baixa confiança em vez de forçar um ajuste.

Resultos
O método foi validado através de simulações numéricas usando dados limpos, ruído de medição (até 40 dB SNR), quantização de 12 bits, desajustes de ganho e atrasos de sensores.

• Precisão: Em casos controlados e sem ruído, o estimador de janela reproduz constantes de ramo analíticas com precisão de máquina. Sob condições realistas (ruído de 60 dB, janelas de 0,5–2 ciclos), os valores recuperados de $R, L$ e $C$ permanecem dentro de uma banda de engenharia de 5%.
• Robustez: O método identifica com sucesso trajetórias informativas enquanto expõe casos estruturalmente degenerados.
• Modos de Falha: O portão de admissibilidade de energia sinaliza efetivamente vieses de tempo coerentes. Por exemplo, um atraso de 50 µs no sensor de corrente fez com que os erros de parâmetro subissem para ~37%, o que o método identificou corretamente como um problema de sincronização em vez de um equivalente confiável.
• Validação: Quinze das vinte e uma configurações de janela prática no estudo de estresse atenderam aos critérios de alta confiança (cobertura $\ge$ 80%, erro relativo $\le$ 5%, resíduo de energia $\le$ 5%).

Significância e Alegações
O artigo alega que seu valor primário reside em determinar um circuito equivalente que, para uma topologia admissível assumida, fecha o balanço de energia terminal instantâneo da carga medida. Os elementos $R, L$ e $C$ identificados formam um "equivalente fisicamente estruturado" em vez de um ajuste de caixa preta.

Os autores enfatizam que o método não pretende determinar unicamente a localização interna da energia de uma carga (pois múltiplas topologias podem satisfazer os dados terminais), mas sim fornecer um equivalente fisicamente consistente para uma topologia escolhida. A significância do trabalho é a introdução de um "portão de confiança": o estimador retorna um equivalente auditável e energeticamente consistente ou declara a janela como não informativa, evitando assim o emprego de parâmetros não confiáveis em sistemas de proteção ou controle. O trabalho é apresentado como um passo em direção à contabilidade de energia em nível de elemento, distinto da estimação de estado do sistema.