Dissertations / Theses on the topic 'Dynamic stochastic optimal power flow'

High penetration of wind energy requires innovations in different areas of power engineering. Methods for improving wind energy and power system interconnection, control, and operation are proposed in this dissertation. A feed-forward transient compensation control scheme is proposed to enhance the low-voltage ride-through capability of wind turbines equipped with doubly fed induction generators. Stator-voltage transient compensation terms are introduced to suppress rotor-current overshoots and torque ripples during grid faults. A dynamic stochastic optimal power flow control scheme is proposed to optimally reroute real-time active and reactive power flow in the presence of high variability and uncertainty. The performance of the proposed power flow control scheme is demonstrated in test power systems with large wind plants. A combined energy-and-reserve wind market scheme is proposed to reduce wind production uncertainty. Variable wind reserve products are created to absorb part of the wind production variation. These fast wind reserve products can then be used to regulate system frequency and improve system security.

Nowadays, power systems are dealing with some new challenges raisedby the major changes that have been taken place since 80’s, e.g., deregu-lation in electricity markets, significant increase of electricity demands andmore recently large-scale integration of renewable energy resources such aswind power. Therefore, system operators must make some adjustments toaccommodate these changes into the future of power systems.One of the main challenges is maintaining the system stability since theextra stress caused by the above changes reduces the stability margin, andmay lead to rise of many undesirable phenomena. The other important chal-lenge is to cope with uncertainty and variability of renewable energy sourceswhich make power systems to become more stochastic in nature, and lesscontrollable.Flexible AC Transmission Systems (FACTS) have emerged as a solutionto help power systems with these new challenges. This thesis aims to ap-propriately utilize such devices in order to increase the transmission capacityand flexibility, improve the dynamic behavior of power systems and integratemore renewable energy into the system. To this end, the most appropriatelocations and settings of these controllable devices need to be determined.This thesis mainly looks at (i) rotor angle stability, i.e., small signal andtransient stability (ii) system operation under wind uncertainty. In the firstpart of this thesis, trajectory sensitivity analysis is used to determine themost suitable placement of FACTS devices for improving rotor angle sta-bility, while in the second part, optimal settings of such devices are foundto maximize the level of wind power integration. As a general conclusion,it was demonstrated that FACTS devices, installed in proper locations andtuned appropriately, are effective means to enhance the system stability andto handle wind uncertainty.The last objective of this thesis work is to propose an efficient solutionapproach based on Benders’ decomposition to solve a network-constrained acunit commitment problem in a wind-integrated power system. The numericalresults show validity, accuracy and efficiency of the proposed approach.

The Doctoral Degrees issued upon completion of the programme are issued by Comillas Pontifical University, Delft University of Technology and KTH Royal Institute of Technology. The invested degrees are official in Spain, the Netherlands and Sweden, respectively.QC 20141028

Massive infrastructure networks such as electric power, natural gas, or water systems play a pivotal role in everyday human lives. Development and operation of these networks is extremely capital-intensive. Moreover, security and reliability of these networks is critical. This work identifies and addresses a diverse class of computationally challenging and time-critical problems pertaining to these networks. This dissertation extends the state of the art on three fronts. First, general proofs of uniqueness for network flow problems are presented, thus addressing open problems. Efficient network flow solvers based on energy function minimizations, convex relaxations, and mixed-integer programming are proposed with performance guarantees. Second, a novel approach is developed for sample-efficient training of deep neural networks (DNN) aimed at solving optimal network dispatch problems. The novel feature here is that the DNNs are trained to match not only the minimizers, but also their sensitivities with respect to the optimization problem parameters. Third, control mechanisms are designed that ensure resilient and stable network operation. These novel solutions are bolstered by mathematical guarantees and extensive simulations on benchmark power, water, and natural gas networks.
Doctor of Philosophy
Massive infrastructure networks play a pivotal role in everyday human lives. A minor service disruption occurring locally in electric power, natural gas, or water networks is considered a significant loss. Uncertain demands, equipment failures, regulatory stipulations, and most importantly complicated physical laws render managing these networks an arduous task. Oftentimes, the first principle mathematical models for these networks are well known. Nevertheless, the computations needed in real-time to make spontaneous decisions frequently surpass the available resources. Explicitly identifying such problems, this dissertation extends the state of the art on three fronts: First, efficient models enabling the operators to tractably solve some routinely encountered problems are developed using fundamental and diverse mathematical tools; Second, quickly trainable machine learning based solutions are developed that enable spontaneous decision making while learning offline from sophisticated mathematical programs; and Third, control mechanisms are designed that ensure a safe and autonomous network operation without human intervention. These novel solutions are bolstered by mathematical guarantees and extensive simulations on benchmark power, water, and natural gas networks.

Recently, the environmental impact of aircraft and rising fuel prices have become an increasing concern in the aviation industry. To address these problems, organizations such as NASA have set demanding goals for reducing aircraft emissions, fuel burn, and noise. In an effort to reach the goals, a movement toward more-electric aircraft and electric propulsion has emerged. With this movement, the number of critical electrical loads on an aircraft is increasing causing power system reliability to be a point of concern. Currently, power system reliability is maintained through the use of back-up power supplies such as batteries and ram-air-turbines (RATs). However, the increasing power requirements for critical loads will quickly outgrow the capacity of the emergency devices. Therefore, reliability needs to be addressed when designing the primary power distribution system. Power system reliability is a function of component reliability and redundancy. Component reliability is often not determined until detailed component design has occurred; however, the amount of redundancy in the system is often set during the system architecting phase. In order to meet the capacity and reliability requirements of future power distribution systems, a method for redundancy allocation during the system architecting phase is needed. This thesis presents an aircraft power system design methodology that is based upon the engineering decision process. The methodology provides a redundancy allocation strategy and quantitative trade-off environment to compare architecture and technology combinations based upon system capacity, weight, and reliability criteria. The methodology is demonstrated by architecting the power distribution system of an aircraft using turboelectric propulsion. The first step in the process is determining the design criteria which includes a 40 MW capacity requirement, a 20 MW capacity requirement for the an engine-out scenario, and a maximum catastrophic failure rate of one failure per billion flight hours. The next step is determining gaps between the performance of current power distribution systems and the requirements of the turboelectric system. A baseline architecture is analyzed by sizing the system using the turboelectric system power requirements and by calculating reliability using a stochastic flow network. To overcome the deficiencies discovered, new technologies and architectures are considered. Global optimization methods are used to find technology and architecture combinations that meet the system objectives and requirements. Lastly, a dynamic modeling environment is constructed to study the performance and stability of the candidate architectures. The combination of the optimization process and dynamic modeling facilitates the selection of a power system architecture that meets the system requirements and objectives.

The contributions of the research are (a) an infrastructure of data acquisition systems that provides the necessary information for an automated EMS system enabling autonomous distributed state estimation, model validation, simplified protection, and seamless integration of other EMS applications, (b) an object-oriented, interoperable, and unified component model that can be seamlessly integrated with a variety of applications of the EMS, (c) a distributed dynamic state estimator (DDSE) based on the proposed data acquisition system and the object-oriented, interoperable, and unified component model, (d) a physically-based synchronous machine model, which is expressed in terms of the actual self and mutual inductances of the synchronous machine windings as a function of rotor position, for the purpose of synchronous machine parameters identification, and (e) a robust and highly efficient algorithm for the optimal power flow (OPF) problem, one of the most important applications of the EMS, based on the validated states and models of the power system provided by the proposed DDSE.

Orientador: Jose Roberto Sanches Mantovani
Resumo: Neste trabalho são propostos modelos matemáticos determinístico e estocástico de programação cônica de segunda ordem em coordenadas retangulares para o problema de fluxo de potência ótimo de geração e controle de potência reativa no sistemas elétricos de potência, considerando as minimização dos custos de geração de energia, perdas ativas da rede e emissão de poluentes no meio ambiente. Os modelos contemplam as principais características físicas e econômicas do problema estudado, assim como os limites operacionais do sistema elétrico. Os modelos são programados em linguagem AMPL e suas soluções são obtidas através do solver comercial CPLEX. Os sistemas testes IEEE30, IEEE118 e ACTIVSg200 são utilizados nas simulações computacionais dos modelos propostos. Os resultados obtidos pelo modelo determinístico desenvolvido são validados através de comparações com os resultados fornecidos pelo software MATPOWER , onde ambos consideram apenas a existência de gerações termoelétricas. No modelo estocástico utiliza-se a técnica de geração de cenários e considera-se um período de um ano (8760 horas), e geradores que utilizam fontes de geração renováveis e não renováveis.
Abstract: In this work we propose deterministic and stochastic mathematical models of second order conical programming in rectangular coordinates for the optimal power flow problem of reactive power generation and control in electric power systems, considering the minimization of energy generation costs, losses networks and emission of pollutants into the environment. The models contemplate the main physical and economic characteristics of the studied problem, as well as the operational limits of the electric system. The models are programmed in AMPL language and their solutions are obtained through the commercial solver CPLEX. The IEEE30, IEEE118 and ACTIVSg200 test systems are used in the computer simulations of the proposed models. The results obtained by the deterministic model developed are validated through comparisons with the results provided by the software MATPOWERR , where both consider only the existence of thermoelectric generations. The stochastic model uses the scenario generation technique and considers a period of one year (8760 hours), and generators using renewable and non-renewable generation sources.
Mestre

Le contexte de nos travaux de thèse est l'intégration de l'énergie éolienne sur les réseaux insulaires. Ces travaux sont soutenus par EDF SEI, l'opérateur électrique des îles françaises. Nous étudions un système éolien-stockage où un système de stockage d'énergie doit aider un producteur éolien à tenir, vis-à-vis du réseau, un engagement de production pris un jour à l'avance. Dans ce contexte, nous proposons une démarche pour l'optimisation du dimensionnement et du contrôle du système de stockage (gestion d'énergie). Comme les erreurs de prévision J+1 de production éolienne sont fortement incertaines, la gestion d'énergie du stockage est un problème d'optimisation stochastique (contrôle optimal stochastique). Pour le résoudre, nous étudions tout d'abord la modélisation des composants du système (modélisation énergétique du stockage par batterie Li-ion ou Sodium-Soufre) ainsi que des entrées (modélisation temporelle stochastique des entrées incertaines). Nous discutons également de la modélisation du vieillissement du stockage, sous une forme adaptée à l'optimisation de la gestion. Ces modèles nous permettent d'optimiser la gestion de l'énergie par la méthode de la programmation dynamique stochastique (SDP). Nous discutons à la fois de l'algorithme et de ses résultats, en particulier de l'effet de la forme des pénalisations sur la loi de gestion. Nous présentons également l'application de la SDP sur des problèmes complémentaires de gestion d'énergie (lissage de la production d'un houlogénérateur, limitation des rampes de production éolienne). Cette étude de l'optimisation de la gestion permet d'aborder l'optimisation du dimensionnement (choix de la capacité énergétique). Des simulations temporelles stochastiques mettent en évidence le fort impact de la structure temporelle (autocorrélation) des erreurs de prévision sur le besoin en capacité de stockage pour atteindre un niveau de performance donné. La prise en compte de paramètres de coût permet ensuite l'optimisation du dimensionnement d'un point de vue économique, en considérant les coûts de l'investissement, des pertes ainsi que du vieillissement. Nous étudions également le dimensionnement du stockage lorsque la pénalisation des écarts à l'engagement comporte un seuil de tolérance. Nous terminons ce manuscrit en abordant la question structurelle de l'interaction entre l'optimisation du dimensionnement et celle du contrôle d'un système de stockage, car ces deux problèmes d'optimisation sont couplés
The context of this PhD thesis is the integration of wind power into the electricity grid of small islands. This work is supported by EDF SEI, the system operator for French islands. We study a wind-storage system where an energy storage is meant to help a wind farm operator fulfill a day-ahead production commitment to the grid. Within this context, we propose an approach for the optimization of the sizing and the control of the energy storage system (energy management). Because day-ahead wind power forecast errors are a major source of uncertainty, the energy management of the storage is a stochastic optimization problem (stochastic optimal control). To solve this problem, we first study the modeling of the components of the system. This include energy-based models of the storage system, with a focus on Lithium-ion and Sodium-Sulfur battery technologies. We then model the system inputs and in particular the stochastic time series like day-ahead forecast errors. We also discuss the modeling of storage aging, using a formulation which is adapted to the control optimization. Assembling all these models enables us to optimize the energy management of the storage system using the stochastic dynamic programming (SDP) method. We introduce the SDP algorithms and present our optimization results, with a special interest for the effect of the shape of the penalty function on the energy control law. We also present additional energy management applications with SDP (mitigation of wind power ramps and smoothing of ocean wave power). Having optimized the storage energy management, we address the optimization of the storage sizing (choice of the rated energy). Stochastic time series simulations show that the temporal structure (autocorrelation) of wind power forecast errors have a major impact on the need for storage capacity to reach a given performance level. Then we combine simulation results with cost parameters, including investment, losses and aging costs, to build a economic cost function for sizing. We also study storage sizing when the penalization of commitment deviations includes a tolerance threshold. We finish this manuscript with a structural study of the interaction between the optimizations of the sizing and the control of an energy storage system, because these two optimization problems are coupled

Submitted by MARCO ALONSO MENESES ORIONDO (mameneseso@gmail.com) on 2016-04-29T17:33:57Z No. of bitstreams: 1 Meneses - Dissertação.pdf: 4087822 bytes, checksum: dc161c85a2ac9500af836619b1f9d6f3 (MD5)
Approved for entry into archive by Felipe Augusto Arakaki (arakaki@reitoria.unesp.br) on 2016-05-02T18:55:33Z (GMT) No. of bitstreams: 1 oriondo_mam_me_ilha.pdf: 4087822 bytes, checksum: dc161c85a2ac9500af836619b1f9d6f3 (MD5)
Made available in DSpace on 2016-05-02T18:55:33Z (GMT). No. of bitstreams: 1 oriondo_mam_me_ilha.pdf: 4087822 bytes, checksum: dc161c85a2ac9500af836619b1f9d6f3 (MD5) Previous issue date: 2016-03-10
Neste trabalho apresenta-se uma versão especializada da meta-heurística Algoritmo dos Fogos de Artifício de Busca Dinâmica com Mutação de Covariância (AFABDMC) para resolver o problema de Fluxo de Potência Ótimo (FPO) em sistemas de transmissão. No algoritmo proposto consideram-se como variáveis contínuas de controle a magnitude da tensão e a geração de potência ativa nas barras de geração e como variáveis de controle discretas o estado de operação dos shunts e a posição do comutador de taps em transformadores. Assim, o modelo para o problema é de Programação Não-Linear Inteira Mista (PNLIM). A estratégia de resolução adotada consiste em controlar, em cada iteração, os valores das variáveis de controle discretas utilizando-se a meta-heurística AFABDMC e a partir desses valores escolhidos pela meta-heurística, resolver um problema de Programação Não-Linear (PNL) que fornece os valores das variáveis de controle contínuas junto com o estado de operação do sistema. A meta-heurística AFABDMC foi escrita em linguagem MATLAB e o modelo do problema em AMPL. Os subproblemas de PNL foram resolvidos utilizando-se o solver KNITRO, sendo que a interface entre o MATLAB e o AMPL foi feita utilizando-se o AMPL API. Foram realizados testes com os sistemas IEEE de 14, 30, 57 e 118 barras e os resultados indicam que a metodologia proposta é capaz de encontrar soluções de muito boa qualidade para o problema.
In this work, a new specialized metaheuristic based on the Dynamic Search Fireworks Algorithm with Covariance Mutation (DSFWACM) is applied on the Optimal Power Flow (OPF) problem in transmission systems. In the proposed method, generator bus voltage magnitudes and active power generation are considered as continuous variables and the operating state of the shunts and transformer taps settings are considered as discrete variables. Thus, the model is a Mixed-Integer Nonlinear Programming (MINLP) problem. The adopted resolution strategy is to control, in each iteration, the value of the discrete control variables using the DSFWACM metaheuristic and from the metaheuristic’s chosen values, solve the Nonlinear Programming (NLP) problem that provides the values of the continuous control variables along with the system’s operation state. The DSFWACM metaheuristic was written in MATLAB and the problem model in AMPL. The NLP sub-problems were solved using the KNITRO solver, and the interface between MATLAB and AMPL was implemented using the AMPL API. Tests were conducted with the IEEE 14, 30, 57 and 118-bus systems and the results show that the proposed method is able to find high quality solutions to the problem.
CNPq: 132374/2011-0

Orientador: Rubén Romero Lázaro
Resumo: Neste trabalho apresenta-se uma versão especializada da meta-heurística Algoritmo dos Fogos de Artifício de Busca Dinâmica com Mutação de Covariância (AFABDMC) para resolver o problema de Fluxo de Potência Ótimo (FPO) em sistemas de transmissão. No algoritmo proposto consideram-se como variáveis contínuas de controle a magnitude da tensão e a geração de potência ativa nas barras de geração e como variáveis de controle discretas o estado de operação dos shunts e a posição do comutador de taps em transformadores. Assim, o modelo para o problema é de Programação Não-Linear Inteira Mista (PNLIM). A estratégia de resolução adotada consiste em controlar, em cada iteração, os valores das variáveis de controle discretas utilizando-se a meta-heurística AFABDMC e a partir desses valores escolhidos pela meta-heurística, resolver um problema de Programação Não-Linear (PNL) que fornece os valores das variáveis de controle contínuas junto com o estado de operação do sistema. A meta-heurística AFABDMC foi escrita em linguagem MATLAB e o modelo do problema em AMPL. Os subproblemas de PNL foram resolvidos utilizando-se o solver KNITRO, sendo que a interface entre o MATLAB e o AMPL foi feita utilizando-se o AMPL API. Foram realizados testes com os sistemas IEEE de 14, 30, 57 e 118 barras e os resultados indicam que a metodologia proposta é capaz de encontrar soluções de muito boa qualidade para o problema.
Abstract: In this work, a new specialized metaheuristic based on the Dynamic Search Fireworks Algorithm with Covariance Mutation (DSFWACM) is applied on the Optimal Power Flow (OPF) problem in transmission systems. In the proposed method, generator bus voltage magnitudes and active power generation are considered as continuous variables and the operating state of the shunts and transformer taps settings are considered as discrete variables. Thus, the model is a Mixed-Integer Nonlinear Programming (MINLP) problem. The adopted resolution strategy is to control, in each iteration, the value of the discrete control variables using the DSFWACM metaheuristic and from the metaheuristic’s chosen values, solve the Nonlinear Programming (NLP) problem that provides the values of the continuous control variables along with the system’s operation state. The DSFWACM metaheuristic was written in MATLAB and the problem model in AMPL. The NLP sub-problems were solved using the KNITRO solver, and the interface between MATLAB and AMPL was implemented using the AMPL API. Tests were conducted with the IEEE 14, 30, 57 and 118-bus systems and the results show that the proposed method is able to find high quality solutions to the problem.
Mestre

碩士
國立雲林科技大學
電機工程系碩士班
102
Optimal power flow (OPF) problem is to find the best solution in a complex and non-linear mathematical formulation. The OPF problem schedule related equipment of power system are regarded as the control variables which satisfies the constraints and minimize total generation cost, total emission and total real power loss. It can save further energy and carbon reduction and also reduce environmental pollution. This paper presents an Improved Artificial Bee Colony Algorithm for solving multi-objective optimal power flow (MOPF) problem. The main purpose is to find the best solution which includes active power output of generation, voltage magnitude of power generation, tap of transformer and shunt capacitor switching states. The three objectives of the OPF problem can be found effectively and converged from the proposed method that can execute update search as simulating the intelligent foraging behavior of a honey bee colony, and added chaotic queues. The load demanded between different intervals in daily is contained in solve OPF problem. The valve point effect of power plants, the generators output limits between different hour by usage of prohibited operating zones and ramp rate limit and other practical limitations are also considered. By considering above conditions, the dynamic multi-objective optimal power flow (DMOPF) problem will be larger and more complicated and it will get closer to the real power system. When dealing with multi-objective optimization problem, one of the objective functions is improved and other objective functions will become to be poor. They can not be improved at the same time because of the contradiction between each other. In this paper, the global Pareto optimal front which is composed by a set of non-dominated solution is found. The operators can select a set of non-dominated solutions appropriately according to different situation. To demonstrate the effectiveness of the proposed method, the MOPF and DMOPF problem are performed on the IEEE 30-BUS, IEEE 57-BUS, IEEE 118-BUS system. The results show the effectiveness of the proposed method to the OPF problem.

Yes
In this paper, a stochastic method is proposed to assess the amount of active and reactive power that can be injected/absorbed to/from grid within a distribution market environment. Also, the impact of wind power penetration on the reactive and active distribution-locational marginal prices is investigated. Market-based active and reactive optimal power flow is used to maximize the social welfare considering uncertainties related to wind speed and load demand. The uncertainties are modeled by Scenario-based approach. The proposed model is examined with 16-bus UK generic distribution system.
Supported by the Higher Education Ministry of Iraqi government.

Yes
In this paper, a stochastic method is proposed to assess the amount of active and reactive power that can be injected/absorbed to/from grid within a distribution market environment. Also, the impact of wind power penetration on the reactive and active distribution-locational marginal prices is investigated. Market-based active and reactive optimal power flow is used to maximize the social welfare considering uncertainties related to wind speed and load demand. The uncertainties are modeled by Scenario-based approach. The proposed model is examined with 16-bus UK generic distribution system.
Supported by the Higher Education Ministry of Iraqi government.

Electric power systems continue to increase in complexity because of the deployment of market mechanisms, the integration of renewable generation and distributed energy resources (DER) (e.g., wind and solar), the penetration of electric vehicles and other price sensitive loads. These revolutionary changes and the consequent increase in uncertainty and dynamicity call for significant modifications to power system operation models including unit commitment (UC), economic load dispatch (ELD) and optimal power flow (OPF). Planning and operation of these ???smart??? electric grids are expected to be impacted significantly, because of the intermittent nature of various supply and demand resources that have penetrated into the system with the recent advances. The main focus of this thesis is on the application of the Affine Arithmetic (AA) method to power system operational problems. The AA method is a very efficient and accurate tool to incorporate uncertainties, as it takes into account all the information amongst dependent variables, by considering their correlations, and hence provides less conservative bounds compared to the Interval Arithmetic (IA) method. Moreover, the AA method does not require assumptions to approximate the probability distribution function (pdf) of random variables. In order to take advantage of the AA method in power flow analysis problems, first a novel formulation of the power flow problem within an optimization framework that includes complementarity constraints is proposed. The power flow problem is formulated as a mixed complementarity problem (MCP), which can take advantage of robust and efficient state-of-the-art nonlinear programming (NLP) and complementarity problems solvers. Based on the proposed MCP formulation, it is formally demonstrated that the Newton-Raphson (NR) solution of the power flow problem is essentially a step of the traditional General Reduced Gradient (GRG) algorithm. The solution of the proposed MCP model is compared with the commonly used NR method using a variety of small-, medium-, and large-sized systems in order to examine the flexibility and robustness of this approach. The MCP-based approach is then used in a power flow problem under uncertainties, in order to obtain the operational ranges for the variables based on the AA method considering active and reactive power demand uncertainties. The proposed approach does not rely on the pdf of the uncertain variables and is therefore shown to be more efficient than the traditional solution methodologies, such as Monte Carlo Simulation (MCS). Also, because of the characteristics of the MCP-based method, the resulting bounds take into consideration the limits of real and reactive power generation. The thesis furthermore proposes a novel AA-based method to solve the OPF problem with uncertain generation sources and hence determine the operating margins of the thermal generators in systems under these conditions. In the AA-based OPF problem, all the state and control variables are treated in affine form, comprising a center value and the corresponding noise magnitudes, to represent forecast, model error, and other sources of uncertainty without the need to assume a pdf. The AA-based approach is benchmarked against the MCS-based intervals, and is shown to obtain bounds close to the ones obtained using the MCS method, although they are slightly more conservative. Furthermore, the proposed algorithm to solve the AA-based OPF problem is shown to be efficient as it does not need the pdf approximations of the random variables and does not rely on iterations to converge to a solution. The applicability of the suggested approach is tested on a large real European power system.

The deployment of Renewable Energy (RE)-based generation has experienced a sustained global growth in the recent decades, driven by many countries' interest in reducing greenhouse gas emissions and dependence on fossil fuel for electricity generation. This trend is also observed in remote off-grid systems (isolated microgrids), where local communities, in an attempt to reduce fossil fuel dependency and associated economic and environmental costs, and to increase availability of electricity, are favouring the installation of RE-based generation. This practice has posed several challenges to the operation of such systems, due to the intermittent and hard-to-predict nature of RE sources. In particular, this thesis addresses the problem of reliable and economic dispatch of isolated microgrids, also known as the energy management problem, considering the uncertain nature of those RE sources, as well as loads. Isolated microgrids feature characteristics similar to those of distribution systems, in terms of unbalanced power flows, significant voltage drops and high power losses. For this reason, detailed three-phase mathematical models of the microgrid system and components are presented here, in order to account for the impact of unbalanced system conditions on the optimal operation of the microgrid. Also, simplified three-phase models of Distributed Energy Resources (DERs) are developed to reduce the level of complexity in small units that have limited impact on the optimal operation of the system, thus reducing the number of equations and variables of the problem. The proposed mathematical models are then used to formulate a novel energy management problem for isolated microgrids, as a deterministic, multi-period, Mixed-Integer Nonlinear Programming (MINLP) problem. The multi-period formulation allows for a proper management of energy storage resources and multi-period constraints associated with the commitment decisions of DERs. In order to obtain solutions of the energy management problem in reasonable computational times for real-time, realistic applications, and to address the uncertainty issues, the proposed MINLP formulation is decomposed into a Mixed-Integer Linear Programming (MILP) problem, and a Nonlinear programming (NLP) problem, in the context of a Model Predictive Control (MPC) approach. The MILP formulation determines the unit commitment decisions of DERs using a simplified model of the network, whereas the NLP formulation calculates the detailed three-phase dispatch of the units, knowing the commitment status. A feedback signal is generated by the NLP if additional units are required to correct reactive power problems in the microgrid, triggering a new calculation MINLP problem. The proposed decomposition and calculation routines are used to design a new deterministic Energy Management System (EMS) based on the MPC approach to handle uncertainties; hence, the proposed deterministic EMS is able to handle multi-period constraints, and account for the impact of future system conditions in the current operation of the microgrid. In the proposed methodology, uncertainty associated with the load and RE-based generation is indirectly considered in the EMS by continuously updating the optimal dispatch solution (with a given time-step), based on the most updated information available from suitable forecasting systems. For a more direct modelling of uncertainty in the problem formulation, the MILP part of the energy management problem is re-formulated as a two-stage Stochastic Programming (SP) problem. The proposed novel SP formulation considers that uncertainty can be properly modelled using a finite set of scenarios, which are generated using both a statistical ensembles scenario generation technique and historical data. Using the proposed SP formulation of the MILP problem, the deterministic EMS design is adjusted to produce a novel stochastic EMS. The proposed EMS design is tested in a large, realistic, medium-voltage isolated microgrid test system. For the deterministic case, the results demonstrate the important connection between the microgrid's imbalance, reactive power requirements and optimal dispatch, justifying the need for detailed three-phase models for EMS applications in isolated microgrids. For the stochastic studies, the results show the advantages of using a stochastic MILP formulation to account for uncertainties associated with RE sources, and optimally accommodate system reserves. The computational times in all simulated cases show the feasibility of applying the proposed techniques to real-time, autonomous dispatch of isolated microgrids with variable RE sources.