Estudo do reaproveitamento energético de ar quente numa Olaria

Como o sector cerâmico é um consumidor intensivo de energia, este trabalho teve como objectivo principal a elaboração de um plano de optimização do desempenho energético da olaria número três da Fábrica Cerâmica de Valadares. Para o efeito, efectuou-se o levantamento energético desta fracção autónoma. O valor total obtido para os ganhos térmicos foi de 8,7x107 kJ/dia, sendo 82% desta energia obtida na combustão do gás natural. Por outro lado, as perdas energéticas rondam os 8,2x107 kJ/dia, sendo o ar de exaustão e a envolvente os principais responsáveis, com um peso de 42 % e 38%, respectivamente. Tendo em conta estes valores, estudaram-se várias medidas de isolamento da cobertura, pavimento, paredes e saída de ar através de fendas do edifício. No caso do isolamento da cobertura sugeriu-se a substituição das telhas de fibrocimento e do isolamento actualmente existentes por painéis sandwich de cobertura. Esta acção permite uma poupança de 64.796€/ano, com um investimento de 57.029€ e o seu período de retorno de 0,9 anos. O Valor Actualizado Líquido (VAL) no 5º ano foi de 184.069€, com uma Taxa Interna de Rentabilidade (TIR) de 92%. Para isolar o pavimento, sugeriu-se a utilização de placas de poliuretano expandido (PU) de 20mm de espessura. Assim, consegue-se uma poupança de 7.442 €/ano, com um investimento de 21.708€, e um tempo de retorno 2,9 anos. No final do 5º ano de vida útil do projecto, o VAL é de 4.070€ e a TIR 7%. Relativamente ao isolamento das paredes e pilares, sugeriu-se a utilização de placas de PU (30mm), recobertas com chapa de ferro galvanizado. O tempo de retorno do investimento é de 1,5 anos, uma vez que, o investimento é de 13.670€ e a poupança anual será de 9.183€. Esta solução apresenta no último ano um VAL de 12.835€ e uma TIR de 22%. No isolamento das fendas do edifício, sugeriu-se a redução de 20% da sua área livre. Esta medida de optimização implica um investimento de 8.000€, revelando-se suficientemente eficaz, pois apresenta um tempo de retorno de 0,67 anos. O VAL e a TIR da solução no último ano de vida útil do projecto de investimento são de 36.835€ e 35%, respectivamente. Por fim, sugeriu-se ainda a instalação de um sistema de controlo que visa o aproveitamento de ar quente proveniente do forno, instalado no piso inferior à olaria, para pré-aquecer o ar alimentado aos geradores de calor. Esta medida implicaria um investimento de 4.000€, com um tempo de retorno de 2,4 anos e uma poupança anual é de 1.686€. O investimento é aconselhável, já que, no 5º ano, o VAL é de 1.956€ e a TIR é de 17%.

Optimização energética da olaria de uma cerâmica

Na indústria cerâmica o consumo de energia é elevado, fazendo com que este custo represente uma parte significativa dos custos totais de produção das peças. De forma a diminuir esta dependência, a energia deve ser gerida de forma contínua e eficazmente. O presente trabalho consistiu na análise da situação energética e na elaboração de propostas de optimização da etapa de conformação que ocorre na Olaria número quatro da Fábrica Cerâmica de Valadares, S.A. Determinou-se o rendimento efectivo da Olaria, tendo-se obtido um valor de 24,7%. As perdas térmicas ocorrem na Olaria, a nível da envolvente, da ventilação, da exaustão de gases e da inércia térmica, representando, respectivamente, 18122 MJ, 50222 MJ, 39228 MJ e 4338 MJ por semana de trabalho. Numa última fase sugeriram-se algumas medidas de optimização energética. A primeira medida visa uma melhoria na manutenção dos geradores, um aumento na gama de temperaturas de funcionamento dos geradores e uma minimização dos tempos de abertura dos portões. Na segunda medida propõe-se a diminuição da percentagem de excesso de ar para 10%, equivalendo a uma poupança de 8839 €/ano. Na terceira medida avaliou-se a possibilidade da aplicação de um permutador de calor de modo a aproveitar os gases de combustão. Esta permitiria uma poupança de 119 €/ano, no entanto, devido ao elevado tempo de retorno do investimento (12,6 anos) considerou-se que esta medida não era viável. A quarta proposta relaciona-se com a optimização da ventilação da Olaria por aumento do ciclo de renovação de ar para 5 h, promovendo uma poupança de 8583 € anuais. Como última sugestão de optimização, aconselhou-se a diminuição do volume da olaria em 6935 m3. Com esta proposta é possível obter uma poupança de 4993 €/ano. Esta medida envolve um investimento de 12000 €, sendo o tempo de retorno do investimento de 2,4 anos. Das cinco propostas estudadas concluiu-se que quatro são viáveis permitindo uma melhoria do funcionamento da Olaria e uma poupança significativa na factura energética.

Estudo da valorização de um resíduo de ácidos gordos livres

Com o desenvolvimento económico a gestão de recursos é um desafio que a sociedade moderna enfrenta. Um desenvolvimento sustentável, que procure satisfazer as necessidades da geração actual e futuras, é fundamental. É pois necessário encontrar novas fontes de energia e reduzir/valorizar resíduos. O objectivo deste trabalho prende-se com a valorização de um resíduo de ácidos gordos livres como fonte alternativa de energia. Estes ácidos gordos são obtidos na purificação primária da glicerina proveniente da produção de biodiesel por via alcalina, cedida pela empresa SOCIPOLE SA. Inicialmente separaram-se os ácidos gordos livres e ésteres presentes na glicerina, e caracterizou-se o produto obtido. Efectuaram-se dois ensaios, obtendo-se 42,9% e 53,8% em massa de ácidos gordos livres e 41,7% e 34,0% de ésteres no 1º e 2º ensaio respectivamente. Além dos ésteres e ácidos gordos livres, no 1º produto existem 15,4% de matéria não esterificável e no 2º 12,2%. Na primeira parte do trabalho, estudou-se a possibilidade de utilização directa deste produto como combustível. Determinou-se o PCI, obtendo-se 29,8 e 29,4 MJ/kg para o 1º e 2º ensaio respectivamente. Estes valores são superiores aos que se encontram no despacho 17313/2008 de 26 de Junho de 2008, para biogasolina e biodiesel (27,2 MJ/kg), significando que os ácidos gordos livres separados têm forte potencial para serem usados como combustível. Efectuou-se também um estudo de corrosão de aço do corpo e tubos de uma caldeira através de ensaios de polarização, usando um Potenciostato/Galvanostato, e por variação da massa de amostras mergulhadas nos ácidos gordos livres a 105ºC e à temperatura ambiente. Nas condições estudadas, os ácidos gordos não promoveram corrosão nas amostras. Por fim, fez-se a combustão de uma amostra dos ácidos gordos livres num bico de Meeker e verificou-se a formação de depósitos equivalentes a 1,93% da massa inicial. Na segunda parte estudou-se a possibilidade de produzir biodiesel por esterificação, dos ácidos gordos livres. Fez-se um estudo prévio para verificar qual a melhor razão molar de amostra/metanol a utilizar, optando-se por 1:6. Após cinco etapas de esterificação, obteve-se um produto final com 85,6% de ésteres. Não foram então atingidos os 96,5% exigidos pela norma europeia EN 14214, devido à presença de matéria não esterificável na matéria-prima, que não foi possível separar totalmente do produto durante o processo (esterificação, lavagem e secagem). Pela análise de todos os parâmetros, o produto obtido não pode ser considerado biodiesel, não sendo possível ser comercializado como tal. A análise dos resultados aponta para a opção da valorização do resíduo de ácidos gordos na queima directa em geradores e caldeiras que deverá no entanto, ser precedida de um estudo mais detalhado que inclua a análise da constituição da matéria não esterificável, para garantir uma combustão segura.

Um novo gerador central de padrões de locomoção para geração de trajectórias em tempo real de robôs quadrúpedes

Neste trabalho estuda-se a geração de trajectórias em tempo real de um robô quadrúpede. As trajectórias podem dividir-se em duas componentes: rítmica e discreta. A componente rítmica das trajectórias é modelada por uma rede de oito osciladores acoplados, com simetria 4 2 Z  Z . Cada oscilador é modelado matematicamente por um sistema de Equações Diferenciais Ordinárias. A referida rede foi proposta por Golubitsky, Stewart, Buono e Collins (1999, 2000), para gerar os passos locomotores de animais quadrúpedes. O trabalho constitui a primeira aplicação desta rede à geração de trajectórias de robôs quadrúpedes. A derivação deste modelo baseia-se na biologia, onde se crê que Geradores Centrais de Padrões de locomoção (CPGs), constituídos por redes neuronais, geram os ritmos associados aos passos locomotores dos animais. O modelo proposto gera soluções periódicas identificadas com os padrões locomotores quadrúpedes, como o andar, o saltar, o galopar, entre outros. A componente discreta das trajectórias dos robôs usa-se para ajustar a parte rítmica das trajectórias. Este tipo de abordagem é útil no controlo da locomoção em terrenos irregulares, em locomoção guiada (por exemplo, mover as pernas enquanto desempenha tarefas discretas para colocar as pernas em localizações específicas) e em percussão. Simulou-se numericamente o modelo de CPG usando o oscilador de Hopf para modelar a parte rítmica do movimento e um modelo inspirado no modelo VITE para modelar a parte discreta do movimento. Variou-se o parâmetro g e mediram-se a amplitude e a frequência das soluções periódicas identificadas com o passo locomotor quadrúpede Trot, para variação deste parâmetro. A parte discreta foi inserida na parte rítmica de duas formas distintas: (a) como um offset, (b) somada às equações que geram a parte rítmica. Os resultados obtidos para o caso (a), revelam que a amplitude e a frequência se mantêm constantes em função de g. Os resultados obtidos para o caso (b) revelam que a amplitude e a frequência aumentam até um determinado valor de g e depois diminuem à medida que o g aumenta, numa curva quase sinusoidal. A variação da amplitude das soluções periódicas traduz-se numa variação directamente proporcional na extensão do movimento do robô. A velocidade da locomoção do robô varia com a frequência das soluções periódicas, que são identificadas com passos locomotores quadrúpedes.

Robôs quadrúpedes: geração de trajectórias em tempo real usando sistemas dinâmicos não-lineares

A geração de trajectórias de robôs em tempo real é uma tarefa muito complexa, não existindo ainda um algoritmo que a permita resolver de forma eficaz. De facto, há controladores eficientes para trajectórias previamente definidas, todavia, a adaptação a variações imprevisíveis, como sendo terrenos irregulares ou obstáculos, constitui ainda um problema em aberto na geração de trajectórias em tempo real de robôs. Neste trabalho apresentam-se modelos de geradores centrais de padrões de locomoção (CPGs), inspirados na biologia, que geram os ritmos locomotores num robô quadrúpede. Os CPGs são modelados matematicamente por sistemas acoplados de células (ou neurónios), sendo a dinâmica de cada célula dada por um sistema de equações diferenciais ordinárias não lineares. Assume-se que as trajectórias dos robôs são constituídas por esta parte rítmica e por uma parte discreta. A parte discreta pode ser embebida na parte rítmica, (a.1) como um offset ou (a.2) adicionada às expressões rítmicas, ou (b) pode ser calculada independentemente e adicionada exactamente antes do envio dos sinais para as articulações do robô. A parte discreta permite inserir no passo locomotor uma perturbação, que poderá estar associada à locomoção em terrenos irregulares ou à existência de obstáculos na trajectória do robô. Para se proceder á análise do sistema com parte discreta, será variado o parâmetro g. O parâmetro g, presente nas equações da parte discreta, representa o offset do sinal após a inclusão da parte discreta. Revê-se a teoria de bifurcação e simetria que permite a classificação das soluções periódicas produzidas pelos modelos de CPGs com passos locomotores quadrúpedes. Nas simulações numéricas, usam-se as equações de Morris-Lecar e o oscilador de Hopf como modelos da dinâmica interna de cada célula para a parte rítmica. A parte discreta é modelada por um sistema inspirado no modelo VITE. Medem-se a amplitude e a frequência de dois passos locomotores para variação do parâmetro g, no intervalo [-5;5]. Consideram-se duas formas distintas de incluir a parte discreta na parte rítmica: (a) como um (a.1) offset ou (a.2) somada nas expressões que modelam a parte rítmica, e (b) somada ao sinal da parte rítmica antes de ser enviado às articulações do robô. No caso (a.1), considerando o oscilador de Hopf como dinâmica interna das células, verifica-se que a amplitude e frequência se mantêm constantes para -50.2. A extensão do movimento varia de forma directamente proporcional à amplitude. No caso das equações de Morris-Lecar, quando a componente discreta é embebida (a.2), a amplitude e a frequência aumentam e depois diminuem para - 0.170.5 Pode concluir-se que: (1) a melhor forma de inserção da parte discreta que menos perturbação insere no robô é a inserção como offset; (2) a inserção da parte discreta parece ser independente do sistema de equações diferenciais ordinárias que modelam a dinâmica interna de cada célula. Como trabalho futuro, é importante prosseguir o estudo das diferentes formas de inserção da parte discreta na parte rítmica do movimento, para que se possa gerar uma locomoção quadrúpede, robusta, flexível, com objectivos, em terrenos irregulares, modelada por correcções discretas aos padrões rítmicos.

Desenvolvimento de uma metodologia baseada em redes neuronais para o cálculo da energia produzida por aerogeradores

O uso da energia eólica para a produção de eletricidade apresenta na última década um crescimento apreciável. Monitorizar o desempenho dos aerogeradores torna-se um processo incontornável, quer por motivos financeiros, quer por questões operacionais. Os investimentos despendidos na construção de parques eólicos são muito consideráveis, pelo que é essencial a análise constante dos aspetos preponderantes no retorno do investimento. A maximização da energia produzida por cada aerogerador é o objetivo principal da monitorização dos parques eólicos. Os sistemas Supervisory Control and Data Acquisition (SCADAs) instalados nos parques eólicos permitem uma supervisão em tempo real relativamente ao estado e funcionamento dos aerogeradores, adquirindo uma elevada importância na avaliação dos rendimentos energéticos e anomalias de funcionamento, garantido desta forma melhorias de produtividade. O objetivo deste trabalho é estimar a energia produzida pelos aerogeradores quando ocorrem falhas de comunicação com o seu contador interno ou avaria do mesmo. A ocorrência destas situações não permite a monitorização da energia produzida durante esse período. Foram analisados dados operacionais dos aerogeradores relativos a um parque eólico localizado na zona Norte de Portugal, sendo usados os dados recolhidos pelo sistema SCADA sobre a forma de médias de 10 min referentes ao período de janeiro de 2011 a agosto 2011. O desempenho da rede neuronal depende da qualidade e quantidade do conjunto de dados usados para o treino da rede. Os dados usados devem representar de forma fiel o estado que se pretende para o equipamento. Para a obtenção do objetivo proposto foi fundamental a identificação das grandezas disponíveis a utilizar no método de cálculo da energia produzida. Os resultados obtidos com aplicação das redes neuronais no método de cálculo da energia produzida por aerogeradores demonstram que independentemente do período de indisponibilidade da informação referente à energia produzida é possível estimar o valor da mesma.

Complex order biped rhythms

Animal locomotion is a complex process, involving the central pattern generators (neural networks, located in the spinal cord, that produce rhythmic patterns), the brainstem command systems, the steering and posture control systems and the top layer structures that decide which motor primitive is activated at a given time. Pinto and Golubitsky studied an integer CPG model for legs rhythms in bipeds. It is a four-coupled identical oscillators' network with dihedral symmetry. This paper considers a new complex order central pattern generator (CPG) model for locomotion in bipeds. A complex derivative Dα±jβ, with α, β ∈ ℜ+, j = √-1, is a generalization of the concept of an integer derivative, where α = 1, β = 0. Parameter regions where periodic solutions, identified with legs' rhythms in bipeds, occur, are analyzed. Also observed is the variation of the amplitude and period of periodic solutions with the complex order derivative.

Escalonamento de geradores térmicos com integração de produção eólica: uma abordagem baseada em programação inteira mista

A alta e crescente participação da energia eólica na matriz da produção traz grandes desafios aos operadores do sistema na gestão da rede e planeamento da produção. A incerteza associada à produção eólica condiciona os processos de escalonamento e despacho económico dos geradores térmicos, uma vez que a produção eólica efetiva pode ser muito diferente da produção prevista. O presente trabalho propõe duas metodologias de otimização do escalonamento de geradores térmicos baseadas em Programação Inteira Mista. Pretende-se encontrar soluções de escalonamento que minimizem as influências negativas da integração de energia eólica no sistema elétrico. Inicialmente o problema de escalonamento de geradores é formulado sem considerar a integração da energia eólica. Posteriormente foi considerada a penetração da energia eólica no sistema elétrico. No primeiro modelo proposto, o problema é formulado como um problema de otimização estocástico. Nesta formulação todos os cenários de produção eólica são levados em consideração no processo de otimização. No segundo modelo, o problema é formulado como um problema de otimização determinística. Nesta formulação, o escalonamento é feito para cada cenário de produção eólica e no fim determina-se a melhor solução por meio de indicadores de avaliação. Foram feitas simulações para diferentes níveis de reserva girante e os resultados obtidos mostraram que a alta participação da energia eólica na matriz da produção põe em causa a segurança e garantia de produção devido às características volátil e intermitente da produção eólica e para manter os mesmos níveis de segurança é preciso dispor no sistema de capacidade reserva girante suficiente capaz de compensar os erros de previsão.

Maximização da valorização energética em aterros sanitários recorrendo a equipamentos baseados no Ciclo Orgânico de Rankine (ORC)

Atualmente, os aterros sanitários representam uma solução para a gestão e tratamento dos resíduos sólidos urbanos. Da deposição, ocorrem duas formas de emissões ao longo do tempo, a produção de biogás e de lixiviados, que resultam sobretudo da decomposição da matéria orgânica. Um dos principais constituintes do biogás é o metano, o qual tem elevado poder calorífico. O presente trabalho aborda, a maximização da valorização energética em aterros sanitários, recorrendo a equipamentos baseados no Ciclo Orgânico de Rankine (ORC) para a produção de eletricidade. É apresentado como caso de estudo a central de valorização energética da Suldouro, em Sermonde, que produz eletricidade a partir do biogás resultante da decomposição da matéria orgânica depositada em aterro. O biogás é utilizado como combustível para os motogeradores utilizados para o seu aproveitamento energético, sendo que apenas cerca de 40% do potencial energético contido no biogás é transformado em eletricidade, registando-se perdas sobretudo nas emissões dos gases de exaustão e na água de arrefecimento dos motores. Para avaliação do potencial da recuperação energética dos gases de escape é avaliado o desempenho termodinâmico do ciclo ORC. Para tal foi desenvolvida uma ferramenta em MATLAB utilizando como modelo a configuração do ORC com recuperador de calor. O cálculo das propriedades termodinâmicas dos fluidos foi obtido através da criação de uma sub-rotina que chama o programa CoolProp. Este programa restitui propriedades como a entalpia, entropia, pressões e temperaturas em cada ponto do ciclo, permitindo assim ao utilizador otimizar o tempo na obtenção de resultados. A avaliação económica é fundamental na tomada de decisões por parte do investidor e dos financiadores do projeto. É então apresentada a análise económica e efetuada uma análise de sensibilidade, onde foram efetuadas variações nos vetores mais importantes de forma a poder avaliar-se o impacto em termos da sua rentabilidade. A ferramenta desenvolvida permite obter de forma prática, os três indicadores económicos extremamente influentes no que se refere à tomada de decisão. A utilização dos sistemas ORC e os seus benefícios não se esgotam na maximização dos aproveitamentos da valorização energética em aterros sanitários. Também a recuperação de calor para a produção de energia elétrica pode ter um impacto importante em muitos setores intensivos de energia, contribuindo significativamente para a redução do consumo e aumentando a eficiência de todo o processo de produção.

Modified Particle Swarm Optimization Applied to Integrated Demand Response and DG Resources Scheduling

The elastic behavior of the demand consumption jointly used with other available resources such as distributed generation (DG) can play a crucial role for the success of smart grids. The intensive use of Distributed Energy Resources (DER) and the technical and contractual constraints result in large-scale non linear optimization problems that require computational intelligence methods to be solved. This paper proposes a Particle Swarm Optimization (PSO) based methodology to support the minimization of the operation costs of a virtual power player that manages the resources in a distribution network and the network itself. Resources include the DER available in the considered time period and the energy that can be bought from external energy suppliers. Network constraints are considered. The proposed approach uses Gaussian mutation of the strategic parameters and contextual self-parameterization of the maximum and minimum particle velocities. The case study considers a real 937 bus distribution network, with 20310 consumers and 548 distributed generators. The obtained solutions are compared with a deterministic approach and with PSO without mutation and Evolutionary PSO, both using self-parameterization.

A Hybrid Simulated Annealing approach to handle Energy Resource Management considering an intensive use of Electric Vehicles

The massification of electric vehicles (EVs) can have a significant impact on the power system, requiring a new approach for the energy resource management. The energy resource management has the objective to obtain the optimal scheduling of the available resources considering distributed generators, storage units, demand response and EVs. The large number of resources causes more complexity in the energy resource management, taking several hours to reach the optimal solution which requires a quick solution for the next day. Therefore, it is necessary to use adequate optimization techniques to determine the best solution in a reasonable amount of time. This paper presents a hybrid artificial intelligence technique to solve a complex energy resource management problem with a large number of resources, including EVs, connected to the electric network. The hybrid approach combines simulated annealing (SA) and ant colony optimization (ACO) techniques. The case study concerns different EVs penetration levels. Comparisons with a previous SA approach and a deterministic technique are also presented. For 2000 EVs scenario, the proposed hybrid approach found a solution better than the previous SA version, resulting in a cost reduction of 1.94%. For this scenario, the proposed approach is approximately 94 times faster than the deterministic approach.

Simulação da Central de Cogeração da Refinaria de Matosinhos (Petrogal)

Mestrado em Engenharia Química – Ramo Optimização Energética na Indústria Química

Virtual Power Player Using Demand Response to Deal with Unexpected Low Wind Power Generation

Demand response is assumed as an essential resource to fully achieve the smart grids operating benefits, namely in the context of competitive markets and of the increasing use of renewable-based energy sources. Some advantages of Demand Response (DR) programs and of smart grids can only be achieved through the implementation of Real Time Pricing (RTP). The integration of the expected increasing amounts of distributed energy resources, as well as new players, requires new approaches for the changing operation of power systems. The methodology proposed in this paper aims the minimization of the operation costs in a distribution network operated by a virtual power player that manages the available energy resources focusing on hour ahead re-scheduling. When facing lower wind power generation than expected from day ahead forecast, demand response is used in order to minimize the impacts of such wind availability change. In this way, consumers actively participate in regulation up and spinning reserve ancillary services through demand response programs. Real time pricing is also applied. The proposed model is especially useful when actual and day ahead wind forecast differ significantly. Its application is illustrated in this paper implementing the characteristics of a real resources conditions scenario in a 33 bus distribution network with 32 consumers and 66 distributed generators.

Maximizing the Social Welfare of Virtual Power Players Operation in Case of Excessive Wind Power

The integration of growing amounts of distributed generation in power systems, namely at distribution networks level, has been fostered by energy policies in several countries around the world, including in Europe. This intensive integration of distributed, non-dispatchable, and natural sources based generation (including wind power) has caused several changes in the operation and planning of power systems and of electricity markets. Sometimes the available non-dispatchable generation is higher than the demand. This generation must be used; otherwise it is wasted if not stored or used to supply additional demand. New policies and market rules, as well as new players, are needed in order to competitively integrate all the resources. The methodology proposed in this paper aims at the maximization of the social welfare in a distribution network operated by a virtual power player that aggregates and manages the available energy resources. When facing a situation of excessive non-dispatchable generation, including wind power, real time pricing is applied in order to induce the increase of consumption so that wind curtailment is minimized. This method is especially useful when actual and day-ahead resources forecast differ significantly. The distribution network characteristics and concerns are addressed by including the network constraints in the optimization model. The proposed methodology has been implemented in GAMS optimization tool and its application is illustrated in this paper using a real 937-bus distribution network with 20.310 consumers and 548 distributed generators, some of them non-dispatchable and with must take contracts. The implemented scenario corresponds to a real day in Portuguese power system.

Distributed Energy Resources Management using Plug-in Hybrid Electric Vehicles as a Fuel-Shifting Demand Response Resource

In the smart grids context, distributed energy resources management plays an important role in the power systems’ operation. Battery electric vehicles and plug-in hybrid electric vehicles should be important resources in the future distribution networks operation. Therefore, it is important to develop adequate methodologies to schedule the electric vehicles’ charge and discharge processes, avoiding network congestions and providing ancillary services. This paper proposes the participation of plug-in hybrid electric vehicles in fuel shifting demand response programs. Two services are proposed, namely the fuel shifting and the fuel discharging. The fuel shifting program consists in replacing the electric energy by fossil fuels in plug-in hybrid electric vehicles daily trips, and the fuel discharge program consists in use of their internal combustion engine to generate electricity injecting into the network. These programs are included in an energy resources management algorithm which integrates the management of other resources. The paper presents a case study considering a 37-bus distribution network with 25 distributed generators, 1908 consumers, and 2430 plug-in vehicles. Two scenarios are tested, namely a scenario with high photovoltaic generation, and a scenario without photovoltaic generation. A sensitivity analyses is performed in order to evaluate when each energy resource is required.