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.

Two significant trends of recent power system evolution are: (1) increasing installa-tion of dynamic loads and distributed generation resources in distribution systems; (2) large-scale renewable energy integration at the transmission system level. A majority of these devices interface with power systems through power electronic converters. However, existing transient stability (TS) simulators are inadequate to represent the dynamic behavior of these devices accurately. On the other hand, simulating a large system using an electromagnetic transient (EMT) simulator is computationally impractical. EMT-TS hybrid simulation approach is an alternative to address these challenges. Furthermore, to thoroughly analyze the increased interactions among the transmission and distribution systems, an integrated modeling and simulation approach is essential.

The thesis is divided into three parts. The first part focuses on an improved hybrid simulation approach and software development. Compared to the previous work, the pro-posed approach has three salient features: three-sequence TS simulation algorithm, three-phase/three-sequence network equivalencing and flexible switching of the serial and par-allel interaction protocols.

The second part of the thesis concentrates on the applications of the hybrid simula-tion tool. The developed platform is first applied to conduct a detailed fault-induced de-layed voltage recovery (FIDVR) study on the Western Electricity Coordinating Council (WECC) system. This study uncovers that after a normally cleared single line to ground fault at the transmission system could cause air conditioner motors to stall in the distribu-tion systems, and the motor stalling could further propagate to an unfaulted phase under certain conditions. The developed tool is also applied to simulate power systems inter-faced with HVDC systems, including classical HVDC and the new generation voltage source converter (VSC)-HVDC system.

The third part centers on the development of integrated transmission and distribution system simulation and an advanced hybrid simulation algorithm with a capability of switching from hybrid simulation mode to TS simulation. Firstly, a modeling framework suitable for integrated transmission and distribution systems is proposed. Secondly, a power flow algorithm and a diakoptics based dynamic simulation algorithm for the integrated transmission and distribution system are developed. Lastly, the EMT-TS hybrid simulation algorithm is combined with the diakoptics based dynamic simulation algorithm to realize flexible simulation mode switching to increase the simulation efficiency.

The standard optimal power flow (OPF) problem is an economic dispatch (ED) problem combined with transmission constraints, which are based on a static topology. However, topology control (TC) has been proposed in the past as a corrective mechanism to relieve overloads and voltage violations. Even though the benefits of TC are presented by several research works in the past, the computational complexity associated with TC has been a major deterrent to its implementation. The proposed work develops heuristics for TC and investigates its potential to improve the computational time for TC for various applications. The objective is to develop computationally light methods to harness the flexibility of the grid to derive maximum benefits to the system in terms of reliability. One of the goals of this research is to develop a tool that will be capable of providing TC actions in a minimal time-frame, which can be readily adopted by the industry for real-time corrective applications.

A DC based heuristic, i.e., a greedy algorithm, is developed and applied to improve the computational time for the TC problem while still maintaining the ability to find quality solutions. In the greedy algorithm, an expression is derived, which indicates the impact on the objective for a marginal change in the state of a transmission line. This expression is used to generate a priority list with potential candidate lines for switching, which may provide huge improvements to the system. The advantage of this method is that it is a fast heuristic as compared to using mixed integer programming (MIP) approach.

Alternatively, AC based heuristics are developed for TC problem and tested on actual data from PJM, ERCOT and TVA. AC based N-1 contingency analysis is performed to identify the contingencies that cause network violations. Simple proximity based heuristics are developed and the fast decoupled power flow is solved iteratively to identify the top five TC actions, which provide reduction in violations. Time domain simulations are performed to ensure that the TC actions do not cause system instability. Simulation results show significant reductions in violations in the system by the application of the TC heuristics.

Present distribution infrastructure is designed mainly for uni-directional power flow with well-controlled generation. An increase in the inverter-interfaced photovoltaic (PV) systems requires a thorough re-examination of the design, operation, protection and control of distribution systems. In order to understand the impact of high penetration of PV generation, this work conducts an automated and detailed modeling of a power distribution system. The simulation results of the modeled distribution feeder have been verified with the field measurements.

Based on the feeder model, this work studies the impact of the PV systems on voltage profiles under various scenarios, including reallocation of the PV systems, reactive power support from the PV inverters, and settings of the load-tap changing transformers in coordination with the PV penetration. Design recommendations have been made based on the simulation results to improve the voltage profiles in the feeder studied.

To carry out dynamic studies related to high penetration of PV systems, this work proposes a differential algebraic equation (DAE) based dynamic modeling and analysis method. Different controllers including inverter current controllers, anti-islanding controllers and droop controllers, are designed and tested in large systems. The method extends the capability of the distribution system analysis tools, to help conduct dynamic analyses in large unbalanced distribution systems.

Another main contribution of this work is related to the investigation of the PV impacts on the feeder protection coordination. Various protection coordination types, including fuse-fuse, recloser-fuse, relay-fuse and relay-recloser have been studied. The analyses provide a better understanding of the relay and recloser settings under different configurations of the PV interconnection transformers, PV penetration levels, and fault types.

A decision tree and fuzzy logic based fault location identification process has also been proposed in this work. The process is composed of the off-line training of the decision tree, and the on-line analysis of the fault events. Fault current contribution from the PV systems, as well as the variation of the fault resistance have been taken into consideration. Two actual fault cases with the event data recorded were used to examine the effectiveness of the fault identification process.

The grounding system in a substation is used to protect personnel and equipment. When there is fault current injected into the ground, a well-designed grounding system should disperse the fault current into the ground in order to limit the touch potential and the step potential to an acceptable level defined by the IEEE Std 80. On the other hand, from the point of view of economy, it is desirable to design a ground grid that minimizes the cost of labor and material. To design such an optimal ground grid that meets the safety metrics and has the minimum cost, an optimal ground grid application was developed in MATLAB, the OptimaL Ground Grid Application (OLGGA).

In the process of ground grid optimization, the touch potential and the step potential are introduced as nonlinear constraints in a two layer soil model whose parameters are set by the user. To obtain an accurate expression for these nonlinear constraints, the ground grid is discretized by using a ground-conductor (and ground-rod) segmentation method that breaks each conductor into reasonable-size segments. The leakage current on each segment and the ground potential rise (GPR) are calculated by solving a matrix equation involving the mutual resistance matrix. After the leakage current on each segment is obtained, the touch potential and the step potential can be calculated using the superposition principle.

A genetic algorithm is used in the optimization of the ground grid and a pattern search algorithm is used to accelerate the convergence. To verify the accuracy of the application, the touch potential and the step potential calculated by the MATLAB application are compared with those calculated by the commercialized grounding system analysis software, WinIGS.

The user's manual of the optimal ground grid application is also presented in this work.

Constraint relaxation by definition means that certain security, operational, or financial constraints are allowed to be violated in the energy market model for a predetermined penalty price. System operators utilize this mechanism in an effort to impose a price-cap on shadow prices throughout the market. In addition, constraint relaxations can serve as corrective approximations that help in reducing the occurrence of infeasible or extreme solutions in the day-ahead markets. This work aims to capture the impact constraint relaxations have on system operational security. Moreover, this analysis also provides a better understanding of the correlation between DC market models and AC real-time systems and analyzes how relaxations in market models propagate to real-time systems. This information can be used not only to assess the criticality of constraint relaxations, but also as a basis for determining penalty prices more accurately.

Constraint relaxations practice was replicated in this work using a test case and a real-life large-scale system, while capturing both energy market aspects and AC real-time system performance. System performance investigation included static and dynamic security analysis for base-case and post-contingency operating conditions. PJM peak hour loads were dynamically modeled in order to capture delayed voltage recovery and sustained depressed voltage profiles as a result of reactive power deficiency caused by constraint relaxations. Moreover, impacts of constraint relaxations on operational system security were investigated when risk based penalty prices are used. Transmission lines in the PJM system were categorized according to their risk index and each category was as-signed a different penalty price accordingly in order to avoid real-time overloads on high risk lines.

This work also extends the investigation of constraint relaxations to post-contingency relaxations, where emergency limits are allowed to be relaxed in energy market models. Various scenarios were investigated to capture and compare between the impacts of base-case and post-contingency relaxations on real-time system performance, including the presence of both relaxations simultaneously. The effect of penalty prices on the number and magnitude of relaxations was investigated as well.

This work presents research on practices in the day-ahead electric energy market, including replication practices and reliability coordinators used by some market operators to demonstrate the impact these practices have on market outcomes. The practice of constraint relaxations similar to those an Independent System Operator (ISO) might perform in day-ahead market models is implemented. The benefits of these practices are well understood by the industry; however, the implications these practices have on market outcomes and system security have not been thoroughly investigated. By solving a day-ahead market model with and without select constraint relaxations and comparing the resulting market outcomes and possible effects on system security, the effect of these constraint relaxation practices is demonstrated.

Proposed market solutions are often infeasible because constraint relaxation practices and approximations that are incorporated into market models. Therefore, the dispatch solution must be corrected to ensure its feasibility. The practice of correcting the proposed dispatch solution after the market is solved is known as out-of-market corrections (OMCs), defined as any action an operator takes that modifies a proposed day-ahead dispatch solution to ensure operating and reliability requirements. The way in which OMCs affect market outcomes is illustrated through the use of different corrective procedures. The objective of the work presented is to demonstrate the implications of these industry practices and assess the impact these practices have on market outcomes.

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.

ABSTRACT

This dissertation introduces a real-time topology monitoring scheme for power systems intended to provide enhanced situational awareness during major system disturbances. The topology monitoring scheme requires accurate real-time topology information to be effective. This scheme is supported by advances in transmission line outage detection based on data-mining phasor measurement unit (PMU) measurements.

A network flow analysis scheme is proposed to track changes in user defined minimal cut sets within the system. This work introduces a new algorithm used to update a previous network flow solution after the loss of a single system branch. The proposed new algorithm provides a significantly decreased solution time that is desired in a real- time environment. This method of topology monitoring can provide system operators with visual indications of potential problems in the system caused by changes in topology.

This work also presents a method of determining all singleton cut sets within a given network topology called the one line remaining (OLR) algorithm. During operation, if a singleton cut set exists, then the system cannot withstand the loss of any one line and still remain connected. The OLR algorithm activates after the loss of a transmission line and determines if any singleton cut sets were created. These cut sets are found using properties of power transfer distribution factors and minimal cut sets.

The topology analysis algorithms proposed in this work are supported by line outage detection using PMU measurements aimed at providing accurate real-time topology information. This process uses a decision tree (DT) based data-mining approach to characterize a lost tie line in simulation. The trained DT is then used to analyze PMU measurements to detect line outages. The trained decision tree was applied to real PMU measurements to detect the loss of a 500 kV line and had no misclassifications.

The work presented has the objective of enhancing situational awareness during significant system disturbances in real time. This dissertation presents all parts of the proposed topology monitoring scheme and justifies and validates the methodology using a real system event.

Understanding the graphical structure of the electric power system is important

in assessing reliability, robustness, and the risk of failure of operations of this criti-

cal infrastructure network. Statistical graph models of complex networks yield much

insight into the underlying processes that are supported by the network. Such gen-

erative graph models are also capable of generating synthetic graphs representative

of the real network. This is particularly important since the smaller number of tradi-

tionally available test systems, such as the IEEE systems, have been largely deemed

to be insucient for supporting large-scale simulation studies and commercial-grade

algorithm development. Thus, there is a need for statistical generative models of

electric power network that capture both topological and electrical properties of the

network and are scalable.

Generating synthetic network graphs that capture key topological and electrical

characteristics of real-world electric power systems is important in aiding widespread

and accurate analysis of these systems. Classical statistical models of graphs, such as

small-world networks or Erd}os-Renyi graphs, are unable to generate synthetic graphs

that accurately represent the topology of real electric power networks { networks

characterized by highly dense local connectivity and clustering and sparse long-haul

links.

This thesis presents a parametrized model that captures the above-mentioned

unique topological properties of electric power networks. Specically, a new Cluster-

and-Connect model is introduced to generate synthetic graphs using these parameters.

Using a uniform set of metrics proposed in the literature, the accuracy of the proposed

model is evaluated by comparing the synthetic models generated for specic real

electric network graphs. In addition to topological properties, the electrical properties

are captured via line impedances that have been shown to be modeled reliably by well-studied heavy tailed distributions. The details of the research, results obtained and

conclusions drawn are presented in this document.