Distributed Load Control in Multiphase Radial Networks

The current power grid is on the cusp of modernization due to the emergence of distributed generation and controllable loads, as well as renewable energy. On one hand, distributed and renewable generation is volatile and difficult to dispatch. On the other hand, controllable loads provide significant potential for compensating for the uncertainties. In a future grid where there are thousands or millions of controllable loads and a large portion of the generation comes from volatile sources like wind and solar, distributed control that shifts or reduces the power consumption of electric loads in a reliable and economic way would be highly valuable.

Load control needs to be conducted with network awareness. Otherwise, voltage violations and overloading of circuit devices are likely. To model these effects, network power flows and voltages have to be considered explicitly. However, the physical laws that determine power flows and voltages are nonlinear. Furthermore, while distributed generation and controllable loads are mostly located in distribution networks that are multiphase and radial, most of the power flow studies focus on single-phase networks.

This thesis focuses on distributed load control in multiphase radial distribution networks. In particular, we first study distributed load control without considering network constraints, and then consider network-aware distributed load control.

Distributed implementation of load control is the main challenge if network constraints can be ignored. In this case, we first ignore the uncertainties in renewable generation and load arrivals, and propose a distributed load control algorithm, Algorithm 1, that optimally schedules the deferrable loads to shape the net electricity demand. Deferrable loads refer to loads whose total energy consumption is fixed, but energy usage can be shifted over time in response to network conditions. Algorithm 1 is a distributed gradient decent algorithm, and empirically converges to optimal deferrable load schedules within 15 iterations.

We then extend Algorithm 1 to a real-time setup where deferrable loads arrive over time, and only imprecise predictions about future renewable generation and load are available at the time of decision making. The real-time algorithm Algorithm 2 is based on model-predictive control: Algorithm 2 uses updated predictions on renewable generation as the true values, and computes a pseudo load to simulate future deferrable load. The pseudo load consumes 0 power at the current time step, and its total energy consumption equals the expectation of future deferrable load total energy request.

Network constraints, e.g., transformer loading constraints and voltage regulation constraints, bring significant challenge to the load control problem since power flows and voltages are governed by nonlinear physical laws. Remarkably, distribution networks are usually multiphase and radial. Two approaches are explored to overcome this challenge: one based on convex relaxation and the other that seeks a locally optimal load schedule.

To explore the convex relaxation approach, a novel but equivalent power flow model, the branch flow model, is developed, and a semidefinite programming relaxation, called BFM-SDP, is obtained using the branch flow model. BFM-SDP is mathematically equivalent to a standard convex relaxation proposed in the literature, but numerically is much more stable. Empirical studies show that BFM-SDP is numerically exact for the IEEE 13-, 34-, 37-, 123-bus networks and a real-world 2065-bus network, while the standard convex relaxation is numerically exact for only two of these networks.

Theoretical guarantees on the exactness of convex relaxations are provided for two types of networks: single-phase radial alternative-current (AC) networks, and single-phase mesh direct-current (DC) networks. In particular, for single-phase radial AC networks, we prove that a second-order cone program (SOCP) relaxation is exact if voltage upper bounds are not binding; we also modify the optimal load control problem so that its SOCP relaxation is always exact. For single-phase mesh DC networks, we prove that an SOCP relaxation is exact if 1) voltage upper bounds are not binding, or 2) voltage upper bounds are uniform and power injection lower bounds are strictly negative; we also modify the optimal load control problem so that its SOCP relaxation is always exact.

To seek a locally optimal load schedule, a distributed gradient-decent algorithm, Algorithm 9, is proposed. The suboptimality gap of the algorithm is rigorously characterized and close to 0 for practical networks. Furthermore, unlike the convex relaxation approach, Algorithm 9 ensures a feasible solution. The gradients used in Algorithm 9 are estimated based on a linear approximation of the power flow, which is derived with the following assumptions: 1) line losses are negligible; and 2) voltages are reasonably balanced. Both assumptions are satisfied in practical distribution networks. Empirical results show that Algorithm 9 obtains 70+ times speed up over the convex relaxation approach, at the cost of a suboptimality within numerical precision.

Análisis de emplazamiento de una central undimotriz en la Costa Vasca.

[ES]El objetivo de este proyecto es intentar conseguir, mediante la energía que ofrece el mar, un sistema de energía eléctrica autónoma para un pueblo de la Costa Vasca. Para ello, es necesario un estudio que permita obtener el flujo de potencia y la dirección óptima donde situar la central undimotriz, así como el mejor sistema para su aprovechamiento.

Linking energy policy, electricity generation and transmission using strong sustainability and co-optimization

The design of a sustainable electricity generation and transmission system is based on the established science of anthropogenic climate change and the realization that depending on imported fossil-fuels is becoming a measure of energy insecurity of supply. A model is proposed which integrates generation fuel mix composition, assignment of plants and optimized power flow, using Portugal as a case study. The result of this co-optimized approach is an overall set of generator types/fuels which increases the diversity of Portuguese electricity supply, lowers its dependency on imported fuels by 14.62% and moves the country towards meeting its regional and international obligations of 31% energy from renewables by 2020 and a 27% reduction in greenhouse gas emissions by 2012, respectively. The quantity and composition of power generation at each bus is specified, with particular focus on quantifying the amount of distributed generation. Based on other works, the resultant, overall distributed capacity penetration of 19.02% of total installed generation is expected to yield positive network benefits. Thus, the model demonstrates that national energy policy and technical deployment can be linked through sustainability and, moreover, that the respective goals may be mutually achieved via holistic, integrated design. ©2009 IEEE.

Continuous-time subspace flows related to the symmetric eigenproblem

The classes of continuous-time flows on Rn×p that induce the same flow on the set of p- dimensional subspaces of Rn×p are described. The power flow is briefly reviewed in this framework, and a subspace generalization of the Rayleigh quotient flow [Linear Algebra Appl. 368C, 2003, pp. 343-357] is proposed and analyzed. This new flow displays a property akin to deflation in finite time. © 2008 Yokohama Publishers.

Thyristors

Thyristors are usually three-terminal devices that have four layers of alternating p-type and n-type material (i.e. three p-n junctions) comprising its main power handling section. In contrast to the linear relation which exists between load and control currents in a transistor, the thyristor is bistable. The control terminal of the thyristor, called the gate (G) electrode, may be connected to an integrated and complex structure as a part of the device. Thyristors are used to approximate ideal closed (no voltage drop between anode and cathode) or open (no anode current flow) switches for control of power flow in a circuit. This differs from low-level digital switching circuits that are designed to deliver two distinct small voltage levels while conducting small currents (ideally zero). Thyristor circuits must have the capability of delivering large currents and be able to withstand large externally applied voltages. All thyristor types are controllable in switching from a forward-lockingstate (positive potential applied to the anode with respect to the cathode, with correspondingly little anode current flow) into a forward-conduction state (large forward anode current flowing, with a small anode-cathode potential drop). Most thyristors have the characteristic that after switching from a forward-blocking state into the forward-conduction state, the gate signal can be removed and the thyristor will remain in its forward-conduction mode. This property is termed "latching" and is an important distinction between thyristors and other types of power electronic devices. © 2007 Elsevier Inc. All rights reserved.

Finite element beam propagation method for analysis of plasmonic waveguide

The basic idea of the finite element beam propagation method (FE-BPM) is described. It is applied to calculate the fundamental mode of a channel plasmonic polariton (CPP) waveguide to confirm its validity. Both the field distribution and the effective index of the, fundamental mode are given by the method. The convergence speed shows the advantage and stability of this method. Then a plasmonic waveguide with a dielectric strip deposited on a metal substrate is investigated, and the group velocity is negative for the fundamental mode of this kind of waveguide. The numerical result shows that the power flow direction is reverse to that of phase velocity.

Analysis of the optical confinement factor in semiconductor lasers

We derive formulas for the optical confinement factor Gamma from Maxwell's equations for TE and TM modes in the slab waveguide. The numerical results show that the formulas yield correct mode gain for the modes propagating in the waveguide. We also compare the formulas with the standard definition of Gamma as the ratio of power flow in the active region to the total power flow. The results show that the standard definition will underestimate the difference of optical confinement factors between TE and TM modes, and will underestimate the difference of material gains necessary for polarization insensitive semiconductor laser amplifiers. It is important to use correct optical confinement factors for designing polarization insensitive semiconductor laser amplifiers. For vertical cavity surface-emitting lasers, the numerical results show that Gamma can be defined as the proportion of the product of the refractive index and the squared electric field in the active region. (C) 1996 American Institute of physics.

Efficient computation algorithm for calculating load contributions to line flows and losses

In a deregulated power system, it is usually required to determine the shares of each load and generation in line flows, to permit fair allocation of transmission costs between the interested parties. The paper presents a new method of determining the contributions of each load to line flows and losses. The method is based on power-flow topology and has the advantage of being the least computationally demanding of similar methods.

Transmission loss allocation in a deregulated electrical energy market

This paper presents a new method for transmission loss allocation in a deregulated electrical power market. The proposed method is based on physical flow through transmission lines. The contributions of individual loads to the line flows are used as basis for allocating transmission losses to different loads. With minimum assumptions, that sound to be reasonable and cannot be rejected, a novel loss allocation formula is derived. The assumptions made are: a number of currents sharing a transmission line distribute themselves over the cross section in the same manner; that distribution causes the minimum possible power loss. Application of the proposed method is straightforward. It requires only a solved power flow and any simple algorithm for power flow tracing. Both active and reactive powers are considered in the loss allocation procedure. Results of application show the accuracy of the proposed method compared with the commonly used procedures.

A method for determining generators’ shares in loads, line flows and losses

This paper presents a new method for calculating the individual generators’ shares in line flows, line losses and loads. The method is described and illustrated on active power flows, but it can be applied in the same way to reactive power flows. Starting from a power flow solution, the line flow matrix is formed. This matrix is used for identifying node types, tracing the power flow from generators downstream to loads, and to determine generators’ participation factors to lines and loads. Neither exhaustive search nor matrix inversion is required. Hence, the method is claimed to be the least computationally demanding amongst all of the similar methods.

Surface plasmons in layered structures with semiconductor and metallic films

The complete spectrum of eigenwaves including surface plasmon polaritons (SPP), dynamic (bulk) and complex waves in the layered structures containing semiconductor and metallic films has been explored. The effects of loss, geometry and the parameters of dielectric layers on the eigenmode spectrum and, particularly, on the SPP modes have been analysed using both the asymptotic and rigorous numerical solutions of the full-wave dispersion equation. The field and Poynting vector distributions have been examined to identify the modes and elucidate their properties. It has been shown that losses and dispersion of permittivity qualitatively alter the spectral content and the eigenwave properties. The SPP counter-directional power fluxes in the film and surrounding dielectrics have been attributed to vortices of power flow, which are responsible for the distinctive features of SPP modes. It has been demonstrated for the first time that the maximal attainable slow-wave factor of the SPP modes guided by thin Au films at optical frequencies is capped not by losses but the frequency dispersion of the actual Au permittivity. © 2009 EDP Sciences.

A Tower-Shadow Model for Wind Turbines Using a Wavelet-Prony Method

The increasing penetration of wind generation on the Island of Ireland has been accompanied by close investigation of low-frequency pulsations contained within active power flow. A primary concern is excitation of low-frequency oscillation modes already present on the system, particularly the 0.75 Hz mode as a consequence of interconnection between the Northern and Southern power system networks. In order to determine whether the prevalence of wind generation has a negative effect (excites modes) or positive impact (damping of modes) on the power system, oscillations must be measured and characterised. Using time – frequency methods, this paper presents work that has been conducted to extract features from low-frequency active power pulsations to determine the composition of oscillatory modes which may impact on dynamic stability. The paper proposes a combined wavelet-Prony method to extract modal components and determine damping factors. The method is exemplified using real data obtained from wind farm measurements.

Wind Farm Induced Oscillations

The impact of power fluctuations arising from fixed-speed wind turbines on the magnitude and frequency of inter-area oscillations was investigated. The authors used data acquisition equipment to record the power flow on the interconnector between the Northern Ireland and Republic of Ireland systems. By monitoring the interconnector oscillation using a fast Fourier transform, it was possible to determine the magnitude and frequency of the inter-area oscillation between the Northern Ireland electricity system and that of the electricity supply board. Analysis was preformed to determine the relationship (if any) between the inter-area oscillation and the observed wind power generation at the corresponding time. Subsequently, regression analysis was introduced to model this relationship between the FFT output and the wind power generation. The effect of conventional generators on the magnitude and frequency of the inter-area oscillation was also considered.

Python PSS/E Simulation to Test Efficacy of Proposed PMU Based WAMS and Potential WAMPAC Applications

A PSS/E 32 model of a real section of the Northern Ireland electrical grid was dynamically controlled with Python 2.5. In this manner data from a proposed wide area monitoring system was simulated. The area is of interest as it is a weakly coupled distribution grid with significant distributed generation. The data was used to create an optimization and protection metric that reflected reactive power flow, voltage profile, thermal overload and voltage excursions. Step changes in the metric were introduced upon the operation of special protection systems and voltage excursions. A wide variety of grid conditions were simulated while tap changer positions and switched capacitor banks were iterated through; with the most desirable state returning the lowest optimization and protection metric. The optimized metric was compared against the metric generated from the standard system state returned by PSS/E. Various grid scenarios were explored involving an intact network and compromised networks (line loss) under summer maximum, summer minimum and winter maximum conditions. In each instance the output from the installed distributed generation is varied between 0 MW and 80 MW (120% of installed capacity). It is shown that in grid models the triggering of special protection systems is delayed by between 1 MW and 6 MW (1.5% to 9% of capacity), with 3.5 MW being the average. The optimization and protection metric gives a quantitative value for system health and demonstrates the potential efficacy of wide area monitoring for protection and control.

Droop setting design for multi-terminal HVDC grids considering voltage deviation impacts

In multi-terminal high voltage direct current (HVDC) grids, the widely deployed droop control strategies will cause a non-uniform voltage deviation on the power flow, which is determined by the network topology and droop settings. This voltage deviation results in an inconsistent power flow pattern when the dispatch references are changed, which could be detrimental to the operation and seamless integration of HVDC grids. In this paper, a novel droop setting design method is proposed to address this problem for a more precise power dispatch. The effects of voltage deviations on the power sharing accuracy and transmission loss are analysed. This paper shows that there is a trade-off between minimizing the voltage deviation, ensuring a proper power delivery and reducing the total transmission loss in the droop setting design. The efficacy of the proposed method is confirmed by simulation studies.