A gestão do risco associado à negociação de energia eléctrica através de contratos bilaterais e mercado spot

Num mercado de electricidade competitivo onde existe um ambiente de incerteza, as empresas de geração adoptam estratégias que visam a maximização do lucro, e a minimização do risco. Neste contexto, é de extrema importância para desenvolver uma estratégia adequada de gestão de risco ter em conta as diferentes opções de negociação de energia num mercado liberalizado, de forma a suportar a tomada de decisões na gestão de risco. O presente trabalho apresenta um modelo que avalia a melhor estratégia de um produtor de energia eléctrica que comercializa num mercado competitivo, onde existem dois mercados possíveis para a transacção de energia: o mercado organizado (bolsa) e o mercado de contratos bilaterais. O produtor tenta maximizar seus lucros e minimizar os riscos correspondentes, seleccionando o melhor equilíbrio entre os dois mercados possíveis (bolsa e bilateral). O mercado de contratos bilaterais visa gerir adequadamente os riscos inerentes à operação de mercados no curto prazo (mercado organizado) e dar o vendedor / comprador uma capacidade real de escolher o fornecedor com que quer negociar. O modelo apresentado neste trabalho faz uma caracterização explícita do risco no que diz respeito ao agente de mercado na questão da sua atitude face ao risco, medido pelo Value at Risk (VaR), descrito neste trabalho por Lucro-em-Risco (PAR). O preço e os factores de risco de volume são caracterizados por um valor médio e um desvio padrão, e são modelizados por distribuições normais. Os resultados numéricos são obtidos utilizando a simulação de Monte Carlo implementado em Matlab, e que é aplicado a um produtor que mantém uma carteira diversificada de tecnologias de geração, para um horizonte temporal de um ano. Esta dissertação está organizada da seguinte forma: o capítulo 1, 2 e 3 descrevem o estado-da-arte relacionado com a gestão de risco na comercialização de energia eléctrica. O capítulo 4 descreve o modelo desenvolvido e implementado, onde é também apresentado um estudo de caso com uma aplicação do modelo para avaliar o risco de negociação de um produtor. No capítulo 5 são apresentadas as principais conclusões.

Planeamento e optimização em redes celulares móveis auto-geridas assentes em LTE

O objectivo deste trabalho consiste em avaliar os benefícios das Self Organizing Networks (SON), no que concerne ao planeamento e optimização de redes Long Term Evolution (LTE), não só através do seu estudo, como também através do desenvolvimento e teste de algoritmos, que permitem avaliar o funcionamento de algumas das suas principais funções. O estudo efectuado sobre as SON permitiu identificar um conjunto de funções, tais como a atribuição automática de Physical Cell Id (PCI), o Automatic Neighbour Relation (ANR) e a optimização automática de parâmetros de handover, que permitem facilitar ou mesmo substituir algumas das tarefas mais comuns em planeamento e optimização de redes móveis celulares, em particular, redes LTE. Recorrendo a um simulador LTE destinado à investigação académica, em código aberto e desenvolvido em Matlab®, foi desenvolvido um conjunto de algoritmos que permitiram a implementação das funções em questão. Para além das funções implementadas, foram também introduzidas alterações que conferem a este simulador a capacidade de representar e simular redes reais, permitindo uma análise mais coerente dos algoritmos desenvolvidos. Os resultados obtidos, para além de evidenciarem claramente o benefício dos algoritmos desenvolvidos, foram ainda comparados com os obtidos pela ferramenta profissional de planeamento e optimização Atoll®, tendo-se verificado a franca proximidade de desempenho em algumas das funções. Finalmente, foi desenvolvida uma interface gráfica que permite o desenho, configuração e simulação de cenários, bem como a análise de resultados.

Conversor multinível NPC de cinco níveis como ondulador de tensão ligado à rede

O crescimento da utilização de accionamentos electromecânicos de velocidade variável entre outros dispositivos que necessitam de tensões elevadas, na ordem dos kV e com elevados níveis de qualidade, despertou o interesse pelos conversores multinível. Este tipo de conversor consegue alcançar elevadas tensões de funcionamento e simultaneamente melhorar a qualidade das formas de onda de tensão e corrente nas respectivas fases. Esta dissertação de mestrado tem por objectivo apresentar um estudo sobre o conversor multinível com díodos de ligação ao neutro (NPC – neutral point clamped), de cinco níveis utilizado como ondulador de tensão ligado à rede. O trabalho começa por desenvolver o modelo matemático do conversor multinível com díodos de ligação ao neutro de cinco níveis e a respectiva interligação com a rede eléctrica. Com base no modelo do conversor são realizadas simulações numéricas desenvolvidas em Matlab-Simulink. Para controlo do trânsito de energia no conversor é utilizando controlo por modo de deslizamento aplicado às correntes nas fases. As simulações efectuadas são comparadas com resultados de simulação obtidos para um ondulador clássico de dois níveis. Resultados de simulação do conversor multinível são posteriormente comparados com resultados experimentais para diferentes valores de potências activa e reactiva. Foi desenvolvido um protótipo experimental de um conversor multinível com díodos de ligação ao neutro de cinco níveis e a respectiva electrónica associada para comando e disparo dos semicondutores de potência.

Caracterização de regiões encefálicas no cérebro normal por tensor de difusão e aplicação à epilepsia pós-traumática

Mestrado em Radiações Aplicadas às Tecnologias da Saúde.

Orientação de campo da máquina assíncrona: controlo preditivo

Nesta dissertação de Mestrado foi realizado o estudo teórico e a respectiva comprovação experimental da aplicação do Controlo por Orientação de Campo para a máquina assíncrona trifásica permitindo a imposição das correntes com recurso ao Controlo Preditivo. Foi efectuado inicialmente a modelização do sistema, ondulador de tensão trifásico e máquina assíncrona trifásica, e o estudo do Controlo Preditivo mais adequado a aplicar ao sistema. Depois de realizado o Controlo por Orientação de Campo foi projectado e dimensionado um anel exterior para controlo de velocidade. Os ensaios com recurso inicialmente à simulação numérica em ambiente Matlab/Simulink e posterior experimentação laboratorial, permitem a validação dos modelos teóricos e da técnica de controlo utilizada. Foram realizados testes em regime estacionário e regime dinâmico para diferentes referências de velocidade. Os ensaios efectuados em simulação e experimentais foram comparados com a solução de imposição das correntes através de controladores histeréticos, nesta dissertação designada por técnica clássica. A implementação laboratorial do protótipo, foi realizada através do projecto, dimensionamento e construção de vários módulos para o efeito. O rectificador trifásico tem entrada compatível com as tensões nominais da rede, inclui também uma resistência de pré-carga dos condensadores, e uma resistência de dissipação para limitar a tensão no barramento de tensão DC. O módulo integrado de potência para funcionamento como ondulador de tensão trifásico tem uma estrutura compacta, onde estão incluídos drivers de disparo dos seus semicondutores, e algumas protecções inerentes ao seu bom funcionamento. Para implementação do Controlo por Orientação de Campo e o anel exterior para controlo de velocidade, foi usado um controlador digital de sinal do fabricante dSPACE.

Modelação do comportamento dinâmico de sistemas barragem-fundação-albufeira: formulação no espaço de estados com amortecimento generalizado

Com este trabalho mostra-se a importância da utilização integrada de modelos numéricos e de resultados da observação do comportamento dinâmico de barragens com vista ao controlo da segurança de sistemas barragem-fundação-albufeira. Neste trabalho utilizou-se o MATLAB para desenvolver um programa de elementos finitos 3D (DySSA 1.0, “Dynamic State Space Analysis”) para análise do comportamento dinâmico de sistemas barragem-fundação-albufeira considerando uma formulação em deslocamentos com a albufeira discretizada em elementos finitos com módulo de distorção nulo. Utiliza-se uma formulação no espaço de estados o que permite considerar amortecimento não proporcional à massa e à rigidez e conduz a modos de vibração complexos. Apresentam-se os fundamentos do método dos elementos finitos na perspetiva da sua implementação computacional para análise dinâmica de estruturas, e apresenta-se sumariamente o programa DySSA 1.0, o qual permite a análise estática e dinâmica de estruturas utilizando elementos finitos tridimensionais isoparamétricos do segundo grau (tipo cubo com 20 pontos nodais), e correspondentes elementos de junta (com 16 pontos nodais, 8 por cada face). O programa foi testado com base na análise dinâmica de uma parede em consola submetida à pressão hidrodinâmica. Os resultados numéricos foram comparados com soluções analíticas e com resultados experimentais obtidos no laboratório de estruturas do ISEL. Descrevem-se os fundamentos da análise dinâmica de estruturas no domínio do tempo, e referem-se os princípios em que se baseiam as metodologias de identificação modal no domínio da frequência. Analisa-se o comportamento dinâmico da barragem do Cabril em termos de frequências naturais e configurações modais, utilizando um modelo 3D representativo do sistema. Os resultados deste modelo numérico são comparados com resultados experimentais obtidos com base no sistema de monitorização dinâmica em contínuo, que foi instalado pelo LNEC em 2008.

Simulação e optimização de tempo de recuperação de um modulador de impulso de alta tensão

presente dissertação tem como objectivo a análise e simulação de um Modulador de Impulso de Alta Tensão (este circuito tem como finalidade substituir o modulador actual do ISOLDE/CERN), utilizando programas de simulação e de implementação de modelos matemáticos respectivamente no Orcad-Pspice e o Matlab-Simulink. Estudou-se o modulador de impulso de alta tensão projectado para o ISOLDE/CERN, e analisou se o seu funcionamento em função do tempo de recuperação e sobreelevação da tensão na carga, os dois parâmetros mais importantes do sistema. Descreveu-se a operação do modulador de impulso de alta tensão baseado no gerador de MARX, desenvolvido no ISEL em pareceria com o CERN, estudou-se o funcionamento deste tipo de circuito baseado na tecnologia de semicondutores. Desenvolveu-se um modelo matemático que descreve a operação do modulador e implementou-se em MATLAB. Fez-se uma análise comparativa dos resultados obtidos experimentalmente com os de simulações obtidos em MATLAB e em PSPICE, e estudaram-se as limitações destes programas na análise do circuito.

Fotodetector e dispositivo WDM integrados

Ao longo deste trabalho é apresentada a caracterização optoelectrónica de uma estrutura semicondutora empilhada de fotodíodos PIN (Positive-Intrinsic-Negative), baseados em silício amorfo hidrogenado (a-Si:H - Hydrogenated Amorphous Silicon) e siliceto de carbono amorfo hi-drogenado (a-SiC:H - Hydrogenated Amorphous Silicon Carbide), em que ambos funcionam como filtros ópticos na zona visível do espectro electromagnético e cuja sensibilidade espectral na região do visível é modulada pelo sinal de tensão eléctrico aplicado e pela presença de polarização óptica adicional (radiação de fundo). Pretende-se utilizar a característica de sensor de cor destes dispositivos semicondutores para realizar a demultiplexagem de sinais ópticos e desenvolver um algoritmo que permita fazer o reco-nhecimento autónomo do sinal transmitido em cada canal, tendo em vista a utilização de vários ca-nais para a transmissão de sinais a curta distância. A transmissão destes sinais deverá ser suportada no meio de transmissão fibra óptica, que constituirá uma importante mais-valia na optimização do sistema WDM (Wavelength Division Mul-tiplexing), permitindo optimizar a transmissão de sinais. Pelas suas capacidades intrínsecas, as fi-bras ópticas de plástico (POF - Plastic Optical Fibers) são uma solução adequada para a transmis-são de sinais no domínio visível do espectro electromagnético a curtas distâncias. Foi realizada uma sucinta caracterização optoelectrónica da estrutura semicondutora sob diferentes condições de iluminação, variando o comprimento de onda e a iluminação de fundo que influencia a resposta espectral do dispositivo semicondutor, variando as cores dos fundos inciden-tes, variando o lado incidente do fundo sobre a estrutura semicondutora, variando a intensidade des-ses mesmos fundos incidentes e também variando a frequência do sinal de dados. Para a transmissão dos sinais de dados foram utilizados três dispositivos LED (Light-Emitting Diode) com as cores vermelho (626nm), verde (525nm) e azul (470nm) a emitir os respec-tivos sinais de dados sobre a estrutura semicondutora e onde foram aplicadas diversas configurações de radiação de fundo incidente, variando as cores dos fundos incidentes, variando o lado incidente do fundo sobre a estrutura semicondutora e variando também a intensidade desses mesmos fundos incidentes. Com base nos resultados obtidos ao longo deste trabalho, foi possível aferir sobre a influên-cia da presença da radiação de fundo aplicada ao dispositivo, usando diferentes sequências de dados transmitidos nos vários canais. Sob polarização inversa, e com a aplicação de um fundo incidente no lado frontal da estrutura semicondutora os valores de fotocorrente gerada são amplificados face aos valores no escuro, sendo que os valores mais altos foram encontrados com a aplicação do fundo de cor violeta, contribuindo para tal, o facto do sinal do canal vermelho e canal verde serem bastan-te amplificados com a aplicação deste fundo. Por outro lado, com a aplicação dos fundos incidentes no lado posterior da estrutura semi-condutora, o sinal gerado não é amplificado com nenhuma cor, no entanto, a aplicação do fundo de cor azul proporciona a distinção do sinal proveniente do canal azul e do canal vermelho, sendo que quando está presente um sinal do canal vermelho, o sinal é fortemente atenuado e com a presença do sinal do canal azul o sinal gerado aproxima-se mais do valor de fotocorrente gerada com a estru-tura no escuro. O algoritmo implementado ao longo deste trabalho, permite efectuar o reconhecimento au-tónomo da informação transmitida por cada canal através da leitura do sinal da fotocorrente forne-cida pelo dispositivo quando sujeito a uma radiação de fundo incidente violeta no lado frontal e uma radiação de fundo incidente azul no lado posterior. Este algoritmo para a descodificação dos sinais WDM utiliza uma aplicação gráfica desenvolvida em Matlab que com base em cálculos e compara-ções de sinal permite determinar a sequência de sinal dos três canais ópticos incidentes. O trabalho proposto nesta tese é um módulo que se enquadra no desenvolvimento de um sistema integrado de comunicação óptica a curta distância, que tem sido alvo de estudo e que resulta das conclusões de trabalhos anteriores, em que este dispositivo e outros de configuração idêntica foram analisados, de forma a explorar a sua utilização na implementação da tecnologia WDM den-tro do domínio do espectro visível e utilizando as POF como meio de transmissão.

Sistema de medida de impedâncias em RF

O objectivo deste trabalho consiste na construção de um protótipo de um aparelho para medir a componente real e imaginária da impedância de dispositivos em RF. Inicialmente é efectuado um estudo do significado de impedância e de outros conceitos a ela associados. Assim como dos vários métodos existentes para a sua medição. Com base nesse estudo é adoptado um tipo de medidor para implementação. Após a implementação do protótipo são efectuados testes de modo a caracterizar o desempenho do mesmo. O protótipo possui um microcontrolador que permite gerir os vários módulos que o compõem e é também o responsável pela comunicação com o utilizador a nível local e com um PC ou terminal remoto. Os dados das medições são enviados para o PC através de uma ligação série para posterior tratamento dos mesmos através de software desenvolvido em ambiente Matlab. Pretende-se também demonstrar que é possível actualmente construir um aparelho de custo moderado mas cujas prestações, em certas aplicações, são adequadas, evitando assim a necessidade de utilizar outros aparelhos de laboratório, normalmente utilizados para este fim, cujo custo é cem vezes superior.

Prensa electromagnética com recuperação de energia para corte e moldagem de chapas e tubos de AL

Foi desenvolvido um conversor de potência e atuador mecânico para a moldagem e corte, por ação de pressão magnética, de chapas e tubos de Al, com uma energia máxima de descarga de 10kJ. O conversor é composto por duas malhas de descarga em paralelo e duas malhas de recuperação de energia. O circuito é capaz de gerar uma corrente de pico de 50kA com uma taxa de variação máxima de 2kA/s através de um atuador, recuperar até 32% da sua energia inicial e diminuir o tempo de carga dos bancos de condensadores no mesmo valor, reduzindo assim a potência da fonte de alimentação primária Foram construídos vários atuadores de forma a otimizar o processo, considerando as várias funções pretendidas, como a deformação ou corte de chapas e compressão de tubos. O circuito elétrico aproximado desenvolvido em Matlab/Simulink foi validado, considerando apenas o funcionamento da malha primária sem o atuador e a dinâmica associada, sendo capaz de simular as respostas do sistema para várias situações específicas, tornando-se numa ferramenta para otimização do mesmo. Experimentalmente, os resultados obtidos provam ser possível cortar chapas Al de 0,5 e 0,8mm, com apenas 13% da energia total do circuito, e comprimir tubos de Al com 2mm de espessura e 50mm de diâmetro com apenas 2,4kJ. A topologia do circuito e a construção da máquina tiveram em conta vários aspetos que melhoram a proteção de pessoas e equipamentos e devida à sua configuração este é capaz de suportar variações de capacidade nos bancos de condensadores e variações de indutância nas bobinas de recuperação de energia sem se danificar.

Reactive power compensation by EPSO technique

This paper presents a methodology to address reactive power compensation using Evolutionary Particle Swarm Optimization (EPSO) technique programmed in the MATLAB environment. The main objective is to find the best operation point minimizing power losses with reactive power compensation, subjected to all operational constraints, namely full AC power flow equations, active and reactive power generation constraints. The methodology has been tested with the IEEE 14 bus test system demonstrating the ability and effectiveness of the proposed approach to handle the reactive power compensation problem.

Optimal dispatch with reactive power compensation by Genetic Algorithm

This paper presents a Unit Commitment model with reactive power compensation that has been solved by Genetic Algorithm (GA) optimization techniques. The GA has been developed a computational tools programmed/coded in MATLAB. The main objective is to find the best generations scheduling whose active power losses are minimal and the reactive power to be compensated, subjected to the power system technical constraints. Those are: full AC power flow equations, active and reactive power generation constraints. All constraints that have been represented in the objective function are weighted with a penalty factors. The IEEE 14-bus system has been used as test case to demonstrate the effectiveness of the proposed algorithm. Results and conclusions are dully drawn.

A probabilistic methodology for distributed generation location in isolated electrical service area

Distributed generation unlike centralized electrical generation aims to generate electrical energy on small scale as near as possible to load centers, interchanging electric power with the network. This work presents a probabilistic methodology conceived to assist the electric system planning engineers in the selection of the distributed generation location, taking into account the hourly load changes or the daily load cycle. The hourly load centers, for each of the different hourly load scenarios, are calculated deterministically. These location points, properly weighted according to their load magnitude, are used to calculate the best fit probability distribution. This distribution is used to determine the maximum likelihood perimeter of the area where each source distributed generation point should preferably be located by the planning engineers. This takes into account, for example, the availability and the cost of the land lots, which are factors of special relevance in urban areas, as well as several obstacles important for the final selection of the candidates of the distributed generation points. The proposed methodology has been applied to a real case, assuming three different bivariate probability distributions: the Gaussian distribution, a bivariate version of Freund’s exponential distribution and the Weibull probability distribution. The methodology algorithm has been programmed in MATLAB. Results are presented and discussed for the application of the methodology to a realistic case and demonstrate the ability of the proposed methodology for efficiently handling the determination of the best location of the distributed generation and their corresponding distribution networks.

Solving a signalized traffic intersection problem with an hyperbolic penalty function

Mathematical Program with Complementarity Constraints (MPCC) finds many applications in fields such as engineering design, economic equilibrium and mathematical programming theory itself. A queueing system model resulting from a single signalized intersection regulated by pre-timed control in traffic network is considered. The model is formulated as an MPCC problem. A MATLAB implementation based on an hyperbolic penalty function is used to solve this practical problem, computing the total average waiting time of the vehicles in all queues and the green split allocation. The problem was codified in AMPL.

Inexact restoration approaches to solve mathematical program with complementarity constraints

Mathematical Program with Complementarity Constraints (MPCC) finds applica- tion in many fields. As the complementarity constraints fail the standard Linear In- dependence Constraint Qualification (LICQ) or the Mangasarian-Fromovitz constraint qualification (MFCQ), at any feasible point, the nonlinear programming theory may not be directly applied to MPCC. However, the MPCC can be reformulated as NLP problem and solved by nonlinear programming techniques. One of them, the Inexact Restoration (IR) approach, performs two independent phases in each iteration - the feasibility and the optimality phases. This work presents two versions of an IR algorithm to solve MPCC. In the feasibility phase two strategies were implemented, depending on the constraints features. One gives more importance to the complementarity constraints, while the other considers the priority of equality and inequality constraints neglecting the complementarity ones. The optimality phase uses the same approach for both algorithm versions. The algorithms were implemented in MATLAB and the test problems are from MACMPEC collection.