A Gray-Box Approach for Decentralized Grid-Equivalent Model Identification

Resumo Técnico: Uma Abordagem de Caixa-Cinza para Identificação de Modelo Equivalente de Rede Descentralizado

Enunciação do Problema
A integração de fontes de energia renovável e eletrônica de potência em redes elétricas tornou necessária uma modelagem baseada em dados robusta e análise de estabilidade. Um desafio crítico neste domínio é a identificação de modelos de impedância da rede. Abordagens centralizadas tradicionais lutam com a sincronização de escala de tempo entre renováveis de baixa inércia e geradores síncronos de alta inércia. Além disso, métodos descentralizados existentes frequentemente dependem de suposições restritivas que falham em redes dinâmicas multiagente. Especificamente, muitos algoritmos atuais exigem identificação de sistema coordenada (onde os conversores excitam sequencialmente), múltiplos pontos de operação ou a suposição irrealista de que a rede se comporta como um barramento infinito "rígido" com uma tensão equivalente perfeitamente constante.

Em uma configuração multi-conversor, sinais locais em um Ponto de Conexão Comum (PCC) são influenciados por entradas de tensão de outros conversores atuando como a "rede". Tentar estimar uma única impedância aparente sem desacoplar essas influências ativas leva a resultados enviesados e imprecisos. O problema central abordado é a necessidade de um algoritmo descentralizado minimamente invasivo, capaz de estimar simultaneamente a impedância equivalente passiva da rede e a tensão equivalente ativa da rede, sem interromper as operações da rede ou exigir comunicação entre conversores.

Metodologia
Os autores propõem um algoritmo de identificação de caixa-cinza operando no domínio da frequência. A abordagem modela a rede a partir da perspectiva de um Conversor de Fonte de Tensão (VSC) local usando um circuito equivalente de Thévenin, onde a rede é representada por uma admitância equivalente passiva ($Y_i(s)$) e uma tensão equivalente ativa ($\Delta \tilde{v}_i(s)$).

A metodologia prossegue através das seguintes etapas:

1. Modelagem do Sistema: A rede é modelada como VSCs interconectados por linhas resistivo-indutivas. As dinâmicas de pequeno sinal são expressas em coordenadas complexas (derivadas de referenciais $dq$). A relação entre a corrente local ($\Delta i_i$) e a tensão ($\Delta v_i$) é formulada como um mapa bilinear envolvendo o vetor de parâmetros desconhecido $\theta_i$ (representando as razões de resistência e indutância da linha) e o fator de acoplamento dinâmico desconhecido $\tilde{h}_i$ (representando a influência da tensão da rede).
2. Estrutura de Caixa-Cinza: O algoritmo distingue entre a impedância da rede "passiva", que é modelada parametricamente (estrutura de baixa ordem, resistivo-indutiva), e a tensão da rede "ativa", que é modelada não parametricamente. Essa distinção reconhece que a tensão ativa da rede surge das reações de malha fechada de outros conversores e cargas, que são desconhecidas e potencialmente complexas.
3. Variáveis Instrumentais (VI) e Discriminação de Frequência: Para abordar a correlação entre a tensão local e a tensão equivalente da rede (que causa viés na estimação por mínimos quadrados padrão), os autores empregam Variáveis Instrumentais. Eles utilizam o sinal de excitação de banda larga local (uma sequência binária pseudoaleatória) como instrumento.
   • Discriminação de Frequência (DF): Um processo de filtragem identifica e descarta pontos de frequência não confiáveis. Isso inclui remover frequências com baixa coerência entre excitação e resposta, frequências dominadas por dinâmicas fundamentais da rede (baixa frequência) ou aliasing/harmônicos (alta frequência), e frequências onde a admitância estimada viola restrições de passividade.
4. Estimação em Duas Etapas:
   • Etapa 1 (Estimação de Parâmetros): Um algoritmo de mínimos quadrados com restrições estima a impedância equivalente paramétrica ($\theta_i$) usando os dados de frequência filtrados e variáveis instrumentais.
   • Etapa 2 (Estimação de Tensão): Com os parâmetros fixos, um filtro de Kalman (formulado como um observador de entrada desconhecida) estima a tensão equivalente da rede não paramétrica ($\Delta \tilde{v}_i$) através do espectro de frequências.

Principais Contribuições
O artigo destaca três contribuições primárias:

1. Execução Totalmente Descentralizada: O algoritmo opera independentemente em cada VSC usando estritamente medições locais do PCC. Requer zero comunicação ou temporização coordenada com outros conversores, tornando-o escalável para redes com fontes e cargas arbitrárias.
2. Separação de Entidades Ativas e Passivas: Ao estimar a impedância paramétrica e a tensão não paramétrica separadamente no domínio da frequência, o método desacopla as dinâmicas físicas da linha de transmissão das reações de controle ativo de outros VSCs. Isso resolve erros graves de acoplamento cruzado que tipicamente corrompem a identificação no domínio do tempo.
3. Robustez: O método é validado sob condições não ideais, incluindo redes de baixa inércia (barramentos não rígidos), excitações simultâneas multiagente, desvios na amplitude e frequência nominais, parâmetros de linha heterogêneos e ruído harmônico.

Resultados
O algoritmo foi testado via simulações numéricas em um sistema interconectado de 5 conversores operando em modo de formação de rede.

• Configuração: 20 experimentos de Monte-Carlo foram conduzidos com parâmetros de controle de droop e impedâncias de linha variáveis. Todos os conversores injetaram excitações mínimas de banda larga (0,001 p.u.) aderindo aos padrões harmônicos IEEE 519.
• Estimação de Impedância: A admitância equivalente estimada ($Y_1(s)$) mostrou alta precisão, com um erro médio de magnitude de 0,26 dB (3,13%) e um erro de fase de 0,07 graus (0,08%).
• Estimação de Tensão: O algoritmo reconstruiu com sucesso a magnitude e a fase da tensão equivalente de pequeno sinal. Os resultados no domínio do tempo (via transformada inversa de Fourier) demonstraram que a tensão estimada correspondia estreitamente aos valores reais em ambas as coordenadas $d$ e $q$.
• Comparação: Os resultados foram comparados contra Estimativas Empíricas de Função de Transferência (ETFE) brutas, que mostraram viés significativo em baixas frequências devido a interações de malha fechada, validando a necessidade do pré-processamento no domínio da frequência e da abordagem de VI proposta.

Significância e Alegações
Os autores afirmam que este trabalho permite a identificação robusta de modelos equivalentes de rede sob condições de rede ativa onde métodos descentralizados anteriores falham. Ao remover a exigência de excitação coordenada ou a suposição de um barramento infinito rígido, a abordagem de caixa-cinza proposta permite análise de estabilidade precisa e projeto de controle adaptativo em sistemas de energia modernos e dinâmicos. O artigo conclui que o modelo identificado pode servir como base para futuros algoritmos de controle adaptativo descentralizado e análise de estabilidade, embora note que a aplicação do algoritmo a hardware do mundo real e sua extensão para incluir geradores síncronos e esquemas de controle diversos são direções para trabalhos futuros.