The objective of this thesis is to detect certain cyber attacks in a power distribution ener-gy management system in a Smart Grid infrastructure. In the Smart Grid, signals are sent be-tween the distribution operator and the customer on a real-time basis. Signals are used for auto-mated energy management, protection and energy metering. This thesis aims at making use of various signals in the system to detect cyber attacks. The focus of the thesis is on a cyber attack that changes the parameters of the energy management system. The attacks considered change the set points, thresholds for energy management decisions, signal multipliers, and other digitally stored parameters that ultimately determine the transfer functions of the components. Since the distribution energy management system is assumed to be in a Smart Grid infrastructure, customer demand is elastic to the price of energy. The energy pricing is represented by a distribution loca-tional marginal price. A closed loop control system is utilized as representative of the energy management system. Each element of the system is represented by a linear transfer function. Studies are done via simulations and these simulations are performed in Matlab SimuLink. The analytical calculations are done using Matlab.

Signals from the system are used to obtain the frequency response of the component transfer functions. The magnitude and phase angle of the transfer functions are obtained using the fast Fourier transform. The transfer function phase angles of base cases (no attack) are stored and are compared with the phase angles calculated at regular time intervals. If the difference in the phase characteristics is greater than a set threshold, an alarm is issued indicating the detection of a cyber attack.

The developed algorithm is designed for use in the envisioned Future Renewable Electric Energy Delivery and Management (FREEDM) system. Examples are shown for the noise free and noisy cases.

This thesis addresses the issue of making an economic case for bulk energy storage in the Arizona bulk power system. Pumped hydro energy storage (PHES) is used in this study. Bulk energy storage has often been suggested for large scale electric power systems in order to levelize load (store energy when it is inexpensive [energy demand is low] and discharge energy when it is expensive [energy demand is high]). It also has the potential to provide opportunities to avoid transmission and generation expansion, and provide for generation reserve margins. As the level of renewable energy resources increases, the uncertainty and variability of wind and solar resources may be improved by bulk energy storage technologies.

For this study, the MATLab software platform is used, a mathematical based modeling language, optimization solvers (specifically Gurobi), and a power flow solver (PowerWorld) are used to simulate an economic dispatch problem that includes energy storage and transmission losses. A program is created which utilizes quadratic programming to analyze various cases using a 2010 summer peak load from the Arizona portion of the Western Electricity Coordinating Council (WECC) system. Actual data from industry are used in this test bed. In this thesis, the full capabilities of Gurobi are not utilized (e.g., integer variables, binary variables). However, the formulation shown here does create a platform such that future, more sophisticated modeling may readily be incorporated.

The developed software is used to assess the Arizona test bed with a low level of energy storage to study how the storage power limit effects several optimization outputs such as the system wide operating costs. Large levels of energy storage are then added to see how high level energy storage affects peak shaving, load factor, and other system applications. Finally, various constraint relaxations are made to analyze why the applications tested eventually approach a constant value. This research illustrates the use of energy storage which helps minimize the system wide generator operating cost by "shaving" energy off of the peak demand.

The thesis builds on the work of another recent researcher with the objectives of strengthening the assumptions used, checking the solutions obtained, utilizing higher level simulation languages to affirm results, and expanding the results and conclusions.

One important point not fully discussed in the present thesis is the impact of efficiency in the pumped hydro cycle. The efficiency of the cycle for modern units is estimated at higher than 90%. Inclusion of pumped hydro losses is relegated to future work.

The safety issue in an electrical power distribution system is of critical importance. In some circumstances, even the continuity of service has to be compromised for a situation that can cause a hazard to the public. A downed conductor that creates an electrical path between a current carrying conductor and ground pose a potential lethal hazard to anyone in the near proximity. Electric utilities have yet to find a fully accepted and reliable method for detecting downed conductors even with decades of research.

With the entry of more automation and a smarter grid in the different layers of distribution power system supply, new doors are being opened and new feasible solutions are waiting to be explored. The 'big data' and the infrastructures that are readily accessible through the smart metering system is the base of the work and analysis performed in this thesis. In effect, the new technologies and new solutions are an artifact of the Smart Grid effort which has now reached worldwide dimensions. A solution to problems of overhead distribution conductor failures / faults that use simple methods and that are easy to implement using existing and future distribution management systems is presented.

A European type distribution system using three phase supply is utilized as the test bed for the concepts presented. Fault analysis is performed on the primary and the secondary distribution system using the free downloadable software OpenDSS. The outcome is a set of rules that can be implemented either locally or central using a voltage based method. Utilized in the distribution management systems the operators will be given a powerful tool to make the correct action when a situation occurs. The test bed itself is taken from an actual system in Norway.

There has been a considerable growth in distributed photovoltaic (PV) genera-tion and its integration in electric power distribution systems. This has led to a change in the distribution system infrastructure. Properly planned distributed gen-eration can offer a variety of benefits for system operations and enhance opera-tional performance of the distribution system. However, high penetration of PV resources can give rise to operating conditions which do not arise in traditional systems and one of the potential issues that needs to be addressed involves impact on power quality of the system with respect to the spectral distortion in voltages and currents.

The test bed feeder model representing a real operational distribution feeder is developed in OpenDSS and the feeder modeling takes into consideration the ob-jective of analysis and frequency of interest. Extensive metering infrastructure and measurements are utilized for validation of the model at harmonic frequencies. The harmonic study performed is divided into two sections: study of impact of non-linear loads on total harmonic voltage and current distortions and study of impact of PV resources on high frequency spectral distortion in voltages and cur-rents. The research work incorporates different harmonic study methodologies such as harmonic and high frequency power flow, and frequency scan study. The general conclusions are presented based on the simulation results and in addition, scope for future work is discussed.

With the power system being increasingly operated near its limits, there is an increasing need for a power-flow (PF) solution devoid of convergence issues. Traditional iterative methods are extremely initial-estimate dependent and not guaranteed to converge to the required solution. Holomorphic Embedding (HE) is a novel non-iterative procedure for solving the PF problem. While the theory behind a restricted version of the method is well rooted in complex analysis, holomorphic functions and algebraic curves, the practical implementation of the method requires going beyond the published details and involves numerical issues related to Taylor's series expansion, Padé approximants, convolution and solving linear matrix equations.

The HE power flow was developed by a non-electrical engineer with language that is foreign to most engineers. One purpose of this document to describe the approach using electric-power engineering parlance and provide an understanding rooted in electric power concepts. This understanding of the methodology is gained by applying the approach to a two-bus dc PF problem and then gradually from moving from this simple two-bus dc PF problem to the general ac PF case.

Software to implement the HE method was developed using MATLAB and numerical tests were carried out on small and medium sized systems to validate the approach. Implementation of different analytic continuation techniques is included and their relevance in applications such as evaluating the voltage solution and estimating the bifurcation point (BP) is discussed. The ability of the HE method to trace the PV curve of the system is identified.

This thesis addresses the issue of making an economic case for energy storage in power systems. Bulk energy storage has often been suggested for large scale electric power systems in order to levelize load; store energy when it is inexpensive and discharge energy when it is expensive; potentially defer transmission and generation expansion; and provide for generation reserve margins. As renewable energy resource penetration increases, the uncertainty and variability of wind and solar may be alleviated by bulk energy storage technologies. The quadratic programming function in MATLAB is used to simulate an economic dispatch that includes energy storage. A program is created that utilizes quadratic programming to analyze various cases using a 2010 summer peak load from the Arizona transmission system, part of the Western Electricity Coordinating Council (WECC). The MATLAB program is used first to test the Arizona test bed with a low level of energy storage to study how the storage power limit effects several optimization out-puts such as the system wide operating costs. Very high levels of energy storage are then added to see how high level energy storage affects peak shaving, load factor, and other system applications. Finally, various constraint relaxations are made to analyze why the applications tested eventually approach a constant value. This research illustrates the use of energy storage which helps minimize the system wide generator operating cost by "shaving" energy off of the peak demand.

The smart grid initiative is the impetus behind changes that are expected to culminate into an enhanced distribution system with the communication and control infrastructure to support advanced distribution system applications and resources such as distributed generation, energy storage systems, and price responsive loads. This research proposes a distribution-class analog of the transmission LMP (DLMP) as an enabler of the advanced applications of the enhanced distribution system. The DLMP is envisioned as a control signal that can incentivize distribution system resources to behave optimally in a manner that benefits economic efficiency and system reliability and that can optimally couple the transmission and the distribution systems. The DLMP is calculated from a two-stage optimization problem; a transmission system OPF and a distribution system OPF. An iterative framework that ensures accurate representation of the distribution system's price sensitive resources for the transmission system problem and vice versa is developed and its convergence problem is discussed. As part of the DLMP calculation framework, a DCOPF formulation that endogenously captures the effect of real power losses is discussed. The formulation uses piecewise linear functions to approximate losses. This thesis explores, with theoretical proofs, the breakdown of the loss approximation technique when non-positive DLMPs/LMPs occur and discusses a mixed integer linear programming formulation that corrects the breakdown. The DLMP is numerically illustrated in traditional and enhanced distribution systems and its superiority to contemporary pricing mechanisms is demonstrated using price responsive loads. Results show that the impact of the inaccuracy of contemporary pricing schemes becomes significant as flexible resources increase. At high elasticity, aggregate load consumption deviated from the optimal consumption by up to about 45 percent when using a flat or time-of-use rate. Individual load consumption deviated by up to 25 percent when using a real-time price. The superiority of the DLMP is more pronounced when important distribution network conditions are not reflected by contemporary prices. The individual load consumption incentivized by the real-time price deviated by up to 90 percent from the optimal consumption in a congested distribution network. While the DLMP internalizes congestion management, the consumption incentivized by the real-time price caused overloads.

Transmission expansion planning (TEP) is a complex decision making process that requires comprehensive analysis to determine the time, location, and number of electric power transmission facilities that are needed in the future power grid. This dissertation investigates the topic of solving TEP problems for large power systems. The dissertation can be divided into two parts. The first part of this dissertation focuses on developing a more accurate network model for TEP study. First, a mixed-integer linear programming (MILP) based TEP model is proposed for solving multi-stage TEP problems. Compared with previous work, the proposed approach reduces the number of variables and constraints needed and improves the computational efficiency significantly. Second, the AC power flow model is applied to TEP models. Relaxations and reformulations are proposed to make the AC model based TEP problem solvable. Third, a convexified AC network model is proposed for TEP studies with reactive power and off-nominal bus voltage magnitudes included in the model. A MILP-based loss model and its relaxations are also investigated. The second part of this dissertation investigates the uncertainty modeling issues in the TEP problem. A two-stage stochastic TEP model is proposed and decomposition algorithms based on the L-shaped method and progressive hedging (PH) are developed to solve the stochastic model. Results indicate that the stochastic TEP model can give a more accurate estimation of the annual operating cost as compared to the deterministic TEP model which focuses only on the peak load.

The subject of this thesis is distribution level load management using a pricing signal in a Smart Grid infrastructure. The Smart Grid implements advanced meters, sensory devices and near real time communication between the elements of the system, including the distribution operator and the customer. A stated objective of the Smart Grid is to use sensory information to operate the electrical power grid more efficiently and cost effectively. One potential function of the Smart Grid is energy management at the distribution level, namely at the individual customer. The Smart Grid allows control of distribution level devices, including distributed energy storage and distributed generation, in operational real time. One method of load control uses an electric energy price as a control signal. The control is achieved through customer preference as the customer allows loads to respond to a dynamic pricing signal. In this thesis, a pricing signal is used to control loads for energy management at the distribution level. The model for the energy management system is created and analyzed in the z-domain due to the envisioned discrete time implementation. Test cases are used to illustrate stability and performance by analytic calculations using Mathcad and by simulation using Matlab Simulink. The envisioned control strategy is applied to the Future Renewable Electric Energy Distribution Management (FREEDM) system. The FREEDM system implements electronic (semiconductor) controls and therefore makes the proposed energy management feasible. The pricing control strategy is demonstrated to be an effective method of performing energy management in a distribution system. It is also shown that stability and near optimal response can be achieved by controlling the parameters of the system. Addition-ally, the communication bandwidth requirements for a pricing control signal are evaluated.

As renewable energy becomes more prevalent in transmission and distribution systems, it is vital to understand the uncertainty and variability that accompany these resources. Microgrids have the potential to mitigate the effects of resource uncertainty. With the ability to exist in either an islanded mode or maintain connections with the main-grid, a microgrid can increase reliability, defer T&D; infrastructure and effectively utilize demand response. This study presents a co-optimization framework for a microgrid with solar photovoltaic generation, emergency generation, and transmission switching. Today unit commitment models ensure reliability with deterministic criteria, which are either insufficient to ensure reliability or can degrade economic efficiency for a microgrid that uses a large penetration of variable renewable resources. A stochastic mixed integer linear program for day-ahead unit commitment is proposed to account for uncertainty inherent in PV generation. The model incorporates the ability to trade energy and ancillary services with the main-grid, including the designation of firm and non-firm imports, which captures the ability to allow for reserve sharing between the two systems. In order to manage the computational complexities, a Benders' decomposition approach is utilized. The commitment schedule was validated with solar scenario analysis, i.e., Monte-Carlo simulations are conducted to test the proposed dispatch solution. For this test case, there were few deviations to power imports, 0.007% of solar was curtailed, no load shedding occurred in the main-grid, and 1.70% load shedding occurred in the microgrid.