With the increasing penetration levels of distributed energy resources along distribution feeders, the importance of load modeling has grown significantly and therefore it is important to have an accurate representation of the distribution system in the planning and operation studies. Although, currently, most of the power system studies are being done using positive sequence commercial software packages for computational convenience purposes, it comes at the cost of reduced accuracy when compared to the more accurate electromagnetic transient (EMT) simulators (but more computationally intensive). However, it is expected, that in the next several years, the use of EMT simulators for large-scale system studies would become a necessity to implement the ambitious renewable energy targets adopted by many countries across the world. Currently, the issue of developing more accurate EMT feeder and load models has yet to be addressed. Therefore, in the first phase of this work, an optimization algorithm to synthesize an EMT distribution feeder and load model has been developed by capturing the current transients when three-phase voltage measurements (obtained from a local utility) are played-in as input, from events such as sub-transmission faults, to the synthesized model. Using the developed algorithm, for the proposed feeder model, both the load composition and the load parameters have been estimated. The synthesized load model has a load composition which includes impedance loads, single-phase induction motor (SPHIM) loads and three-phase induction motor loads. In the second phase of this work, an analytical formulation of a 24 V EMT contactor is developed to trip the air conditioner EMT SPHIM load, in the feeder and load model developed in Phase 1 of this work, under low voltage conditions. Additionally, a new methodology is developed, to estimate and incorporate the trip and reconnection settings of the proposed EMT contactor model to trip, reconnect and stall the SPHIMs in a positive sequence simulator (PSLF) for single-line to ground faults. Also, the proposed methodology has been tested on a modified three-segment three-phase feeder model using a local utility’s practical feeder topological and loading information. Finally, the developed methodology is modified to accommodate three-phase faults in the system.

Ensuring reliable operation of large power systems subjected to multiple outages is a challenging task because of the combinatorial nature of the problem. Traditional methods of steady-state security assessment in power systems involve contingency analysis based on AC or DC power flows. However, power flow based contingency analysis is not fast enough to evaluate all contingencies for real-time operations. Therefore, real-time contingency analysis (RTCA) only evaluates a subset of the contingencies (called the contingency list), and hence might miss critical contingencies that lead to cascading failures.This dissertation proposes a new graph-theoretic approach, called the feasibility test (FT) algorithm, for analyzing whether a contingency will create a saturated or over-loaded cut-set in a meshed power network; a cut-set denotes a set of lines which if tripped separates the network into two disjoint islands. A novel feature of the proposed approach is that it lowers the solution time significantly making the approach viable for an exhaustive real-time evaluation of the system. Detecting saturated cut-sets in the power system is important because they represent the vulnerable bottlenecks in the network. The robustness of the FT algorithm is demonstrated on a 17,000+ bus model of the Western Interconnection (WI). Following the detection of post-contingency cut-set saturation, a two-component methodology is proposed to enhance the reliability of large power systems during a series of outages. The first component combines the proposed FT algorithm with RTCA to create an integrated corrective action (iCA), whose goal is to secure the power system against post-contingency cut-set saturation as well as critical branch overloads. The second component only employs the results of the FT to create a relaxed corrective action (rCA) that quickly secures the system against saturated cut-sets. The first component is more comprehensive than the second, but the latter is computationally more efficient. The effectiveness of the two components is evaluated based upon the number of cascade triggering contingencies alleviated, and the computation time. Analysis of different case-studies on the IEEE 118-bus and 2000-bus synthetic Texas systems indicate that the proposed two-component methodology enhances the scope and speed of power system security assessment during multiple outages.

The holomorphic embedding method (HEM) applied to the power-flow problem (HEPF) has been used in the past to obtain the voltages and flows for power systems. The incentives for using this method over the traditional Newton-Raphson based nu-merical methods lie in the claim that the method is theoretically guaranteed to converge to the operable solution, if one exists.

In this report, HEPF will be used for two power system analysis purposes:

a. Estimating the saddle-node bifurcation point (SNBP) of a system

b. Developing reduced-order network equivalents for distribution systems.

Typically, the continuation power flow (CPF) is used to estimate the SNBP of a system, which involves solving multiple power-flow problems. One of the advantages of HEPF is that the solution is obtained as an analytical expression of the embedding parameter, and using this property, three of the proposed HEPF-based methods can es-timate the SNBP of a given power system without solving multiple power-flow prob-lems (if generator VAr limits are ignored). If VAr limits are considered, the mathemat-ical representation of the power-flow problem changes and thus an iterative process would have to be performed in order to estimate the SNBP of the system. This would typically still require fewer power-flow problems to be solved than CPF in order to estimate the SNBP.

Another proposed application is to develop reduced order network equivalents for radial distribution networks that retain the nonlinearities of the eliminated portion of the network and hence remain more accurate than traditional Ward-type reductions (which linearize about the given operating point) when the operating condition changes.

Different ways of accelerating the convergence of the power series obtained as a part of HEPF, are explored and it is shown that the eta method is the most efficient of all methods tested.

The local-measurement-based methods of estimating the SNBP are studied. Non-linear Thévenin-like networks as well as multi-bus networks are built using model data to estimate the SNBP and it is shown that the structure of these networks can be made arbitrary by appropriately modifying the nonlinear current injections, which can sim-plify the process of building such networks from measurements.

Electric power system security assessment is one of the most important requirements for operational and resource planning of the bulk power system ensuring safe operation of the power system for all credible contingencies. This deterministic approach usually provides a conservative criterion and can result in expensive bulk system expansion plans or conservative operating limits. Furthermore, with increased penetration of converter-based renewable generation in the electric grid, the dynamics of the grid are changing. In addition, the variability and intermittency associated with the renewable energy sources introduce uncertainty in the electricity grid. Since security margins have direct economic impact on the utilities; more clarity is required regarding the basis on which security decisions are made. The main objective of this work is to provide an approach for risk-based security assessment (RBSA) to define dynamic reliability standards in future electricity grids. RBSA provides a measure of the security of the power system that combines both the likelihood and the consequence of an event.

A novel approach to estimate the impact of transient stability is presented by modeling several important protection systems within the transient stability analysis. A robust operational metric to quantify the impact of transient instability event is proposed that incorporates the effort required to stabilize any transiently unstable event. The effect of converter-interfaced renewable energy injection on system reliability is investigated us-ing RBSA. A robust RBSA diagnostics tool is developed which provides an interactive user interface where the RBSA results and contingency ranking reports can be explored and compared based on specific user inputs without executing time domain simulations or risk calculations, hence providing a fast and robust approach for handling large time domain simulation and risk assessment data. The results show that RBSA can be used effectively in system planning to select security limits. Comparison of RBSA with deterministic methods show that RBSA not only provides less conservative results, it also illustrates the bases on which such security decisions are made. RBSA helps in identifying critical aspects of system reliability that is not possible using the deterministic reliability techniques.

Alternate sources of energy such as wind, solar photovoltaic and fuel cells are coupled to the power grid with the help of solid state converters. Continued deregulation of the power sector coupled with favorable government incentives has resulted in the rapid growth of renewable energy sources connected to the distribution system at a voltage level of 34.5kV or below. Of late, many utilities are also investing in these alternate sources of energy with the point of interconnection with the power grid being at the transmission level. These converter interfaced generation along with their associated control have the ability to provide the advantage of fast control of frequency, voltage, active, and reactive power. However, their ability to provide stability in a large system is yet to be investigated in detail. This is the primary objective of this research.

In the future, along with an increase in the percentage of converter interfaced renewable energy sources connected to the transmission network, there exists a possibility of even connecting synchronous machines to the grid through converters. Thus, all sources of energy can be expected to be coupled to the grid through converters. The control and operation of such a grid will be unlike anything that has been encountered till now. In this dissertation, the operation and behavior of such a grid will be investigated. The first step in such an analysis will be to build an accurate and simple mathematical model to represent the corresponding components in commercial software. Once this bridge has been crossed, conventional machines will be replaced with their solid state interfaced counterparts in a phased manner. At each stage, attention will be devoted to the control of these sources and also on the stability performance of the large power system.

This dissertation addresses various concerns regarding the control and operation of a futuristic power grid. In addition, this dissertation also aims to address the issue of whether a requirement may arise to redefine operational reliability criteria based on the results obtained.

A nonlinear dynamic model for a passively cooled small modular reactor (SMR) is developed. The nuclear steam supply system (NSSS) model includes representations for reactor core, steam generator, pressurizer, hot leg riser and downcomer. The reactor core is modeled with the combination of: (1) neutronics, using point kinetics equations for reactor power and a single combined neutron group, and (2) thermal-hydraulics, describing the heat transfer from fuel to coolant by an overall heat transfer resistance and single-phase natural circulation. For the helical-coil once-through steam generator, a single tube depiction with time-varying boundaries and three regions, i.e., subcooled, boiling, and superheated, is adopted. The pressurizer model is developed based upon the conservation of fluid mass, volume, and energy. Hot leg riser and downcomer are treated as first-order lags. The NSSS model is incorporated with a turbine model which permits observing the power with given steam flow, pressure, and enthalpy as input. The overall nonlinear system is implemented in the Simulink dynamic environment. Simulations for typical perturbations, e.g., control rod withdrawal and increase in steam demand, are run. A detailed analysis of the results show that the steady-state values for full power are in good agreement with design data and the model is capable of predicting the dynamics of the SMR. Finally, steady-state control programs for reactor power and pressurizer pressure are also implemented and their effect on the important system variables are discussed.

As the world embraces a sustainable energy future, alternative energy resources, such as wind power, are increasingly being seen as an integral part of the future electric energy grid. Ultimately, integrating such a dynamic and variable mix of generation requires a better understanding of renewable generation output, in addition to power grid systems that improve power system operational performance in the presence of anticipated events such as wind power ramps. Because of the stochastic, uncontrollable nature of renewable resources, a thorough and accurate characterization of wind activity is necessary to maintain grid stability and reliability. Wind power ramps from an existing wind farm are studied to characterize persistence forecasting errors using extreme value analysis techniques. In addition, a novel metric that quantifies the amount of non-stationarity in time series wind power data was proposed and used in a real-time algorithm to provide a rigorous method that adaptively determines training data for forecasts. Lastly, large swings in generation or load can cause system frequency and tie-line flows to deviate from nominal, so an anticipatory MPC-based secondary control scheme was designed and integrated into an automatic generation control loop to improve the ability of an interconnection to respond to anticipated large events and fluctuations in the power system.

This dissertation presents innovative techniques to develop performance-based models and complete transient models of induction motor drive systems with vector controls in electro-magnetic transient (EMT) and positive sequence transient stability (PSTS) simulation programs. The performance-based model is implemented by obtaining the characteristic transfer functions of perturbed active and reactive power consumptions with respect to frequency and voltage perturbations. This level of linearized performance-based model is suitable for the investigation of the damping of small-magnitude low-frequency oscillations. The complete transient model is proposed by decomposing the motor, converter and control models into d-q axes components and developing a compatible electrical interface to the positive-sequence network in the PSTS simulators. The complete transient drive model is primarily used to examine the system response subject to transient voltage depression considering increasing penetration of converter-driven motor loads.

For developing the performance-based model, modulations are performed on the supply side of the full drive system to procure magnitude and phase responses of active and reactive powers with respect to the supply voltage and frequency for a range of discrete frequency points. The prediction error minimization (PEM) technique is utilized to generate the curve-fitted transfer functions and corresponding bode plots. For developing the complete drive model in the PSTS simulation program, a positive-sequence voltage source is defined properly as the interface of the model to the external system. The dc-link of the drive converter is implemented by employing the average model of the PWM converter, and is utilized to integrate the line-side rectifier and machine-side inverter.

Numerical simulation is then conducted on sample test systems, synthesized with suitable characteristics to examine performance of the developed models. The simulation results reveal that with growing amount of drive loads being distributed in the system, the small-signal stability of the system is improved in terms of the desirable damping effects on the low-frequency system oscillations of voltage and frequency. The transient stability of the system is also enhanced with regard to the stable active power and reactive power controls of the loads, and the appropriate VAr support capability provided by the drive loads during a contingency.

Power flow calculation plays a significant role in power system studies and operation. To ensure the reliable prediction of system states during planning studies and in the operating environment, a reliable power flow algorithm is desired. However, the traditional power flow methods (such as the Gauss Seidel method and the Newton-Raphson method) are not guaranteed to obtain a converged solution when the system is heavily loaded.

This thesis describes a novel non-iterative holomorphic embedding (HE) method to solve the power flow problem that eliminates the convergence issues and the uncertainty of the existence of the solution. It is guaranteed to find a converged solution if the solution exists, and will signal by an oscillation of the result if there is no solution exists. Furthermore, it does not require a guess of the initial voltage solution.

By embedding the complex-valued parameter α into the voltage function, the power balance equations become holomorphic functions. Then the embedded voltage functions are expanded as a Maclaurin power series, V(α). The diagonal Padé approximant calculated from V(α) gives the maximal analytic continuation of V(α), and produces a reliable solution of voltages. The connection between mathematical theory and its application to power flow calculation is described in detail.

With the existing bus-type-switching routine, the models of phase shifters and three-winding transformers are proposed to enable the HE algorithm to solve practical large-scale systems. Additionally, sparsity techniques are used to store the sparse bus admittance matrix. The modified HE algorithm is programmed in MATLAB. A study parameter β is introduced in the embedding formula βα + (1- β)α^2. By varying the value of β, numerical tests of different embedding formulae are conducted on the three-bus, IEEE 14-bus, 118-bus, 300-bus, and the ERCOT systems, and the numerical performance as a function of β is analyzed to determine the “best” embedding formula. The obtained power-flow solutions are validated using MATPOWER.

For a (N+1)-bus power system, possibly 2N solutions exists. One of these solutions

is known as the high-voltage (HV) solution or operable solution. The rest of the solutions

are the low-voltage (LV), or large-angle, solutions.

In this report, a recently developed non-iterative algorithm for solving the power-

flow (PF) problem using the holomorphic embedding (HE) method is shown as

being capable of finding the HV solution, while avoiding converging to LV solutions

nearby which is a drawback to all other iterative solutions. The HE method provides a

novel non-iterative procedure to solve the PF problems by eliminating the

non-convergence and initial-estimate dependency issues appeared in the traditional

iterative methods. The detailed implementation of the HE method is discussed in the

report.

While published work focuses mainly on finding the HV PF solution, modified

holomorphically embedded formulations are proposed in this report to find the

LV/large-angle solutions of the PF problem. It is theoretically proven that the proposed

method is guaranteed to find a total number of 2N solutions to the PF problem

and if no solution exists, the algorithm is guaranteed to indicate such by the oscillations

in the maximal analytic continuation of the coefficients of the voltage power series

obtained.

After presenting the derivation of the LV/large-angle formulations for both PQ

and PV buses, numerical tests on the five-, seven- and 14-bus systems are conducted

to find all the solutions of the system of nonlinear PF equations for those systems using

the proposed HE method.

After completing the derivation to find all the PF solutions using the HE method, it

is shown that the proposed HE method can be used to find only the of interest PF solutions

(i.e. type-1 PF solutions with one positive real-part eigenvalue in the Jacobian

matrix), with a proper algorithm developed. The closet unstable equilibrium point

(UEP), one of the type-1 UEP’s, can be obtained by the proposed HE method with

limited dynamic models included.

The numerical performance as well as the robustness of the proposed HE method is

investigated and presented by implementing the algorithm on the problematic cases and

large-scale power system.