Nesta pesquisa, é desenvolvido um sistema autônomo de recarga de VE alimentado por um sistema fotovoltaico bifacial e armazenamento de energia em bateria (BESS), conforme mostrado na Figura 1, utilizando um método de MPPT controlado por RNA. Um circuito equivalente de um único diodo é usado para modelar o arranjo fotovoltaico bifacial, levando em conta a temperatura, resistência série, resistência de derivação, irradiância solar e irradiância no lado traseiro. A saída do sistema fotovoltaico é fornecida à estação de recarga de VE por meio de um conversor CC–CC elevador, enquanto o BESS auxilia na regulação da tensão do barramento CC. O Algoritmo de Otimização do Secretário (SBOA) é utilizado para otimizar os ganhos do controlador PI e os parâmetros do filtro. A modelagem dos componentes do sistema e a metodologia proposta são descritas nas seguintes subseções.
Modelagem de Componentes
1. Sistema fotovoltaico bifacial
Os módulos fotovoltaicos ajudam a alcançar a corrente e tensão de saída exigidas. No presente estudo, dois módulos fotovoltaicos monofaciais foram utilizados para emular as características de geração de energia de matrizes bifaciais sob diferentes condições de albedo da superfície9. O ganho bifacial da irradiância BGg é determinado pela Equação (1)
(1)
Aqui, a irradiação frontal é Gf, e a irradiância traseira é Gr. A corrente máxima fotovoltaica bifacial é Imp, e a tensão é Vmp. Finalmente, a saída produzida pelo sistema bifacial com o número escolhido de módulos em série Ns e em paralelo Np é dada pela Equação (2). O controle do sistema fotovoltaico é apresentado na Figura 2
(2)
2. Sistema de armazenamento de bateria (BSS)
As baterias são um meio crucial para armazenar energia elétrica18. Este trabalho utiliza baterias de íons de lítio como sistema de armazenamento de energia. Entre as opções disponíveis de baterias recarregáveis, este trabalho concentra-se nas baterias de íons de lítio, que são consideradas ideais para veículos elétricos devido à sua eficiência e densidade de potência. As equações (3) e (4) descrevem a bateria de íons de lítio com restrições.
(3)
As especificações do PV, da bateria de armazenamento e do conversor são resumidas na Tabela 2, enquanto a Tabela 3 apresenta a alocação do fluxo de potência entre o conjunto fotovoltaico, a unidade de armazenamento de energia e as cargas conectadas.
(4)
3. Conversor CC –CA
O inversor com ponto neutro clampado (NPC) é amplamente adotado em sistemas modernos de conversão de energia para fontes renováveis. Em comparação com a estrutura convencional de dois níveis, a configuração NPC oferece menor distorção harmônica total (THD), menor interferência eletromagnética e menor tensão sobre os dispositivos semicondutores. O elo CC dessa topologia é dividido em três níveis de potencial: +Vdc/2, 0 e -Vdc/2, por meio de capacitores emparelhados e diodos clampadores que formam a junção neutra. Controlando adequadamente a comutação, o inversor gera três estados de tensão discretos nos terminais de saída, melhorando assim a qualidade da forma de onda e a eficiência da conversão de energia. O ponto neutro do barramento CC em um inversor NPC desempenha um papel fundamental no equilíbrio das tensões dos capacitores e na garantia do funcionamento adequado do inversor. O esquema de controle do inversor é mostrado na Figura 3.
4. Conversor Boost CC-CC
Conversores CC-CC elevam a tensão de entrada contínua para uma tensão de saída mais alta controlando o ciclo de trabalho de um dispositivo de comutação de alta frequência, mantendo a eficiência. Ele possui dois modos de operação. Durante o estado LIGADO, o diodo está polarizado inversamente e a fonte de entrada fornece energia ao indutor, que a armazena como um campo magnético, enquanto o capacitor de saída fornece energia à carga. Durante o estado DESLIGADO, o diodo fica polarizado diretamente. Nesse caso, a energia armazenada no indutor é liberada através do diodo para o capacitor conectado na saída e à carga, combinando-se com a fonte de entrada para produzir uma tensão de saída maior que a tensão de entrada. Variando o ciclo de trabalho do PWM, a tensão de saída do conversor varia de acordo com a relação
, onde D é o ciclo de trabalho. É utilizado em sistemas solares, de células de combustível, veículos elétricos e de armazenamento.
Esquema de controle de RNA para PMP
Neste trabalho, um MPPT baseado em RNA foi selecionado para coletar a potência máxima gerada pelo sistema solar17. Os conjuntos de dados de irradiação e temperatura utilizados para o treinamento da RNA foram gerados no MATLAB/Simulink sob condições operacionais ambientais variáveis, correspondentes a níveis de irradiância solar de 800–1000 W/m2 e faixas de temperatura de 20–25 °C, com um ciclo de trabalho (D). A RNA é treinada para minimizar o erro quadrático médio (EQM) entre as saídas obtidas e as desejadas (Op,
) para extrair a potência máxima do sistema fotovoltaico bifacial. Aqui, a irradiação solar e a temperatura são utilizadas como entradas para a RNA, e o ciclo de trabalho é a saída para controlar o conversor elevador. A estrutura do modelo de RNA desenvolvido para MPPT, conforme apresentada na Figura 4, é considerada neste trabalho. Tabela Suplementar 1 apresenta as vantagens e desvantagens do método proposto em comparação com outros métodos padrão.
Uma RNA consiste em três camadas principais: uma camada de saída (OPL), uma camada de entrada (IPL) e uma camada oculta (HIL). Isso ajuda a transferir os dados entre a IPL e a HIL. Posteriormente, ela é produzida pelos pesos nas conexões que ligam a IPL e a HIL. Nesse contexto, os cálculos são executados com um certo viés aplicado à variável HL, e os resultados obtidos são acumulados na variável OPL. Aqui, é selecionada uma RNA do tipo LMBP [17]. Os pesos das conexões são ajustados durante o treinamento mediante a medição do erro para alcançar a saída desejada. Neste caso, o LMBP é utilizado para treinar RNAs com o MSE como função de desempenho. O algoritmo LMBP utiliza as derivadas obtidas para atualizar os pesos, os quais exibem as propriedades de aprendizado eficaz e convergência acelerada.
Cada neurônio em uma rede de perceptron multicamadas possui uma função de soma e uma função de ativação. No entanto, existem certos pesos numéricos (wpk) que conectam esses neurônios entre os níveis. Quando as entradas são multiplicadas pelos pesos conforme especificado na Equação (5). A função sigmoide não linear é considerada com a minimização do erro médio quadrático dada na Equação (6).
(5)
(6)
SBOA otimizou os parâmetros do filtro e de controle do sistema proposto
Atualmente, os algoritmos meta-heurísticos desempenham um papel fundamental na resolução de problemas de engenharia. Figura 5 apresenta uma classificação dos algoritmos. O algoritmo SBOA23 imita o comportamento do SB para sobreviver em seu habitat natural. As fases de exploração e exploração do comportamento de caça do Secretário-Ave são modeladas no SBOA. A fase de exploração do algoritmo imita o comportamento do SB ao capturar cobras, enquanto sua fase de exploração imita o comportamento de desviar de predadores como águias. O processo de solução começa gerando valores aleatórios dentro de seus respectivos limites para cada SB na população e calculando o valor da função objetivo executando o modelo Simulink, tratando o valor de cada SB como um parâmetro de projeto. Com base no valor objetivo, cada SB altera sua posição para representar uma solução melhor. A modelagem matemática do SBOA é discutida abaixo:
Fase inicial de preparação
Primeiro, é necessário identificar as soluções iniciais que foram utilizadas para iniciar a busca por uma minimização típica da função objetivo F(Y). Nesse caso, a população aleatória inicial Y = [Y1, Y2, ......YN] de SB's é formada por meio de N soluções iniciais. A Equação (7) inicializa a população e Yi representa uma solução.
(7)
Onde os limites inferiores e superiores das variáveis de escolha são denotados por lb e ub. Em [0, 1], r é um número aleatório. N é a dimensão do problema. Além disso, o valor de aptidão da solução Fi = F(Yi) é usado para quantificar a qualidade de Yi.
Estratégia de caça dos SB
Encontrar, comer e atacar a presa são as três fases principais da caça do SB. Todo o processo de caça foi dividido em três intervalos de tempo iguais, t < 1/3T, 1/3T < t < 2/3T e 2/3T < t < T. Aqui, t é a iteração atual e T é o número máximo de iterações. Esses intervalos correspondem às três fases da predação do SB: busca, consumo e ataque à presa. Essas divisões são baseadas nas estatísticas biológicas das fases de caça da ave e nas durações temporais de cada fase. Consequentemente, cada fase do SBOA é modelada da seguinte forma:
Busca por presas
O SB deve caçar presas a uma distância segura durante esta fase. Para coletar dados suficientes em toda a região de busca, o primeiro passo nos algoritmos de otimização exige uma exploração mais intensa. Utilizando as posições dos outros dois SBs como orientação, o SB pode investigar regiões adicionais possíveis. Para aumentar ainda mais a diversidade do algoritmo, introduzem-se assim processos de mutação diferencial. Quando t é menor que 1/3 de T, a Equação (8) é usada para atualizar a posição de cada indivíduo. Ele procurará, consumirá e atacará presas. Consequentemente, cada fase do SBOA é modelada da seguinte forma:
(8)
Onde, Yi(t) é a posição atual do iº SB, Ynewi(t) fornece a solução atualizada, dois indivíduos Yi(t), Yr2(t) foram escolhidos aleatoriamente da população existente. Um vetor aleatório, para a estratégia de seleção de parâmetros, R1 possui 1 × N escolhido aleatoriamente de [0, 1] para um equilíbrio eficaz entre as fases de exploração e exploração. A solução atualizada é considerada somente se resultar em um valor objetivo melhorado em comparação com a solução previamente obtida.
Consumo de presas
O SB utiliza agilidade ágil e se move para pairar ao redor da cobra após detectar uma presa potencial. A paciência da presa será usada para diminuir suas defesas, observando e atraindo os oponentes enquanto circulam. Outros SB ajustam sua postura próximo à presa. O sucesso na caça será significativamente aumentado com essa abordagem. A Equação (9) ilustra como o movimento browniano (B) é usado para simular o movimento aleatório do SB quando 1/3 T < t < 2/3 T.
(9)
Em que o B é implementado no MATLAB usando a distribuição normal padrão por meio do vetor gerado aleatoriamente r(1,M). Em seguida, os SB utilizam a Equação (10) para atualizar suas posições.
(10)
Em que Ybest(t) é a melhor solução e o fator de escala exponencial exp((t/T)4) é escolhido para aumentar a capacidade de exploração a cada iteração.
Atacando a presa
Após consumir continuamente, a vítima esgotará sua energia. As aves-secretárias devem lançar o ataque agora. Aqui, diferentes métodos de ataque, como passos incessantes/saltos longos esporádicos em um curto período de tempo, são simulados usando a técnica de voo de Lévy. As características são descritas pelas Equações (11) e (12). As opções candidatas estão próximas da resposta ótima no momento, pois as SB se aproximam rapidamente da presa. Este plano será implementado quando t > 2/3T.
(11)
(12)
Onde, LF representa a estratégia de voo de Levy e
, s= 0,01 e φ = 1,5 são valores fixos. μ e υ são valores arbitrários que estão no intervalo [0, 1]. Onde, τ denota gama e η = 0,5, o fator de ponderação 0,5 foi adotado, conforme proposto na formulação original do SBOA, para regular o tamanho do passo e manter um equilíbrio entre grandes saltos exploratórios e refinamento local. Tabela 4 apresenta os valores utilizados para a RNA e o algoritmo SB no sistema desenvolvido.
Estratégia de fuga para SB
Ao caçar outras presas na natureza, as AVs correm o risco de serem caçadas. Águias, falcões, raposas e chacais são os principais adversários com os quais precisam lidar. Elas devem usar uma variedade de técnicas de evasão para se protegerem ou protegerem suas presas ao perceberem perigo. Este algoritmo simula os métodos de fuga ao modelar modos de corrida (D2) e camuflagem (D1). Ao enfrentar adversários, as avoantes inicialmente se misturam ao ambiente para permanecerem seguras. Aqui, as avoantes procuram evitar ótimos locais no algoritmo ao atualizarem suas posições ao redor da presa (o melhor indivíduo). Elas usarão táticas de voo ou corrida rápida para permanecerem seguras caso não consigam evadir o inimigo. Para referência, um indivíduo aleatório Yrand é escolhido como líder para evitar ficar restrito a um ótimo local. As Equações (13) e (14) são usadas pelas avoantes para atualizar suas posições. O fluxograma do SBOA é mostrado na Figura 6.
(13)
(14)
Representação das variáveis de projeto
As variáveis do problema neste estudo incluem os parâmetros de projeto, como os ganhos do controlador PI de conversores CC-CC, bem como as resistências e indutâncias de filtros e conversores. A Equação (15) apresenta a representação das variáveis de controle. Os limites dos parâmetros de controle estão listados na Tabela Suplementar 2.
(15)
Função de aptidão (FF)
A minimização da THD é selecionada como função objetivo (Obj) neste estudo. O valor máximo de FF corresponde ao valor mínimo de Obj, dado nas Equações (16) e (17).
(16)
(17)