Power systems are transforming into more complex and stressed systems each day. These stressed conditions could lead to a slow decline in the power grid's voltage profile and sometimes lead to a partial or total blackout. This phenomenon can be identified by either solving a power flow problem or using measurement-based real-time monitoring algorithms. The first part of this thesis focuses on proposing a robust power flow algorithm for ill-conditioned systems. While preserving the stable nature of the fixed point (FP) method, a novel distributed FP equation is proposed to calculate the voltage at each bus. The proposed algorithm's performance is compared with existing methods, showing that the proposed method can correctly find the solutions when other methods cannot work due to high condition number matrices. It is also empirically shown that the FP algorithm is more robust to bad initialization points. The second part of this thesis focuses on identifying the voltage instability phenomenon using real-time monitoring algorithms. This work proposes a novel distributed measurement-based monitoring technique called voltage stability index (VSI). With the help of PMUs and communication of voltage phasors between neighboring buses, the processors embedded at each bus in the smart grid perform simultaneous online computations of VSI. VSI enables real-time identification of the system's critical bus with minimal communication infrastructure. Its benefits include interpretability, fast computation, and low sensitivity to noisy measurements. Furthermore, this work proposes the ``local static-voltage stability index" (LS-VSI) that removes the minimal communication requirement in VSI by requiring only one PMU at the bus of interest. LS-VSI also solves the issue of Thevenin equivalent parameter estimation in the presence of noisy measurements. Unlike VSI, LS-VSI incorporates the ZIP load models and load tap changers (LTCs) and successfully identifies the bifurcation point considering ZIP loads' impact on voltage stability. Both VSI and LS-VSI are useful to monitor the voltage stability margins in real-time using the PMU measurements from the field. However, they cannot indicate the onset of voltage emergency situations. The proposed LD-VSI uses the dynamic measurements of the power system to identify the onset of a voltage emergency situation with an alarm. Compared to existing methods, it is shown that it is more robust to PMU measurement noise and can also identify the voltage collapse point while the existing methods have issues with the same.

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.

Low frequency oscillations (LFOs) are recognized as one of the most challenging problems in electric grids as they limit power transfer capability and can result in system instability. In recent years, the deployment of phasor measurement units (PMUs) has increased the accessibility to time-synchronized wide-area measurements, which has, in turn, enabledthe effective detection and control of the oscillatory modes of the power system. This work assesses the stability improvements that can be achieved through the coordinated wide-area control of power system stabilizers (PSSs), static VAr compensators (SVCs), and supplementary damping controllers (SDCs) of high voltage DC (HVDC) lines, for damping electromechanical oscillations in a modern power system. The improved damping is achieved by designing different types of coordinated wide-area damping controllers (CWADC) that employ partial state-feedback. The first design methodology uses a linear matrix inequality (LMI)-based mixed H2/Hinfty control that is robust for multiple operating scenarios. To counteract the negative impact of communication failure or missing PMU measurements on the designed control, a scheme to identify the alternate set of feedback signals is proposed. Additionally, the impact of delays on the performance of the control design is investigated. The second approach is motivated by the increasing popularity of artificial intelligence (AI) in enhancing the performance of interconnected power systems. Two different wide-area coordinated control schemes are developed using deep neural networks (DNNs) and deep reinforcement learning (DRL), while accounting for the uncertainties present in the power system. The DNN-CWADC learns to make control decisions using supervised learning; the training dataset consisting of polytopic controllers designed with the help of LMI-based mixed H2/Hinfty optimization. The DRL-CWADC learns to adapt to the system uncertainties based on its continuous interaction with the power system environment by employing an advanced version of the state-of-the-art deep deterministic policy gradient (DDPG) algorithm referred to as bounded exploratory control-based DDPG (BEC-DDPG). The studies performed on a 29 machine, 127 bus equivalent model of theWestern Electricity Coordinating Council (WECC) system-embedded with different types of damping controls have demonstrated the effectiveness and robustness of the proposed CWADCs.

The availability of data for monitoring and controlling the electrical grid has increased exponentially over the years in both resolution and quantity leaving a large data footprint. This dissertation is motivated by the need for equivalent representations of grid data in lower-dimensional feature spaces so that machine learning algorithms can be employed for a variety of purposes. To achieve that, without sacrificing the interpretation of the results, the dissertation leverages the physics behind power systems, well-known laws that underlie this man-made infrastructure, and the nature of the underlying stochastic phenomena that define the system operating conditions as the backbone for modeling data from the grid.

The first part of the dissertation introduces a new framework of graph signal processing (GSP) for the power grid, Grid-GSP, and applies it to voltage phasor measurements that characterize the overall system state of the power grid. Concepts from GSP are used in conjunction with known power system models in order to highlight the low-dimensional structure in data and present generative models for voltage phasors measurements. Applications such as identification of graphical communities, network inference, interpolation of missing data, detection of false data injection attacks and data compression are explored wherein Grid-GSP based generative models are used.

The second part of the dissertation develops a model for a joint statistical description of solar photo-voltaic (PV) power and the outdoor temperature which can lead to better management of power generation resources so that electricity demand such as air conditioning and supply from solar power are always matched in the face of stochasticity. The low-rank structure inherent in solar PV power data is used for forecasting and to detect partial-shading type of faults in solar panels.

The electric power system (EPS) is an extremely complex system that has operational interdependencies with the water delivery and treatment system (WDTS). The term water-energy nexus is commonly used to describe the critical interdependencies that naturally exist between the EPS and water distribution systems (WDS). Presented in this work is a framework for simulating interactions between these two critical infrastructure systems in short-term and long-term time-scales. This includes appropriate mathematical models for system modeling and for optimizing control of power system operation with consideration of conditions in the WDS. Also presented is a complete methodology for quantifying the resilience of the two interdependent systems.

The key interdependencies between the two systems are the requirements of water for the cooling cycle of traditional thermal power plants as well as electricity for pumping and/or treatment in the WDS. While previous work has considered the dependency of thermoelectric generation on cooling water requirements at a high-level, this work considers the impact from limitations of cooling water into network simulations in both a short-term operational framework as well as in the long-term planning domain.

The work completed to set-up simulations in operational length time-scales was the development of a simulator that adequately models both systems. This simulation engine also facilitates the implementation of control schemes in both systems that take advantage of the knowledge of operating conditions in the other system. Initial steps for including the influence of anticipated water availability and water rights attainability within the combined generation and transmission expansion planning problem is also presented. Lastly, the framework for determining the infrastructural-operational resilience (IOR) of the interdependent systems is formulated.

Adequately modeling and studying the two systems and their interactions is becoming critically important. This importance is illustrated by the possibility of unforeseen natural or man-made events or by the likelihood of load increase in the systems, either of which has the risk of putting extreme stress on the systems beyond that experienced in normal operating conditions. Therefore, this work addresses these concerns with novel modeling and control/policy strategies designed to mitigate the severity of extreme conditions in either system.

Power systems are undergoing a significant transformation as a result of the retirements of conventional coal-fired generation units and the increasing integration of converter interfaced renewable resources. The instantaneous renewable generation penetration as a percentage of the load served in megawatt (MW), in some areas of the United States (U.S.) sometimes approaches over 50 percent. These changes have introduced new challenges for reliability studies considering the two functional reliability aspects, i.e., adequacy and the dynamic security or operating reliability.

Adequacy assessment becomes more complex due to the variability introduced by renewable energy generation. The traditionally used reserve margin only considers projected peak demand and would be inadequate since it does not consider an evaluation of off-peak conditions that could also be critical due to the variable renewable generation. Therefore, in order to address the impact of variable renewable generation, a probabilistic evaluation that studies all hours of a year based on statistical characteristics is a necessity to identify the adequacy risks. On the other hand, the system dynamic behavior is also changing. Converter interfaced generation resources have different dynamic characteristics from the conventional synchronous units and inherently do not participate in grid regulation functions such as frequency control and voltage control that are vital to maintaining operating reliability. In order to evaluate these evolving grid characteristics, comprehensive reliability evaluation approaches that consider system stochasticity and evaluate both adequacy and dynamic security are important to identify potential system risks in this transforming environment.

The grand transition of electric grids from conventional fossil fuel resources to intermittent bulk renewable resources and distributed energy resources (DERs) has initiated a paradigm shift in power system operation. Distributed energy resources (i.e. rooftop solar photovoltaic, battery storage, electric vehicles, and demand response), communication infrastructures, and smart measurement devices provide the opportunity for electric utility customers to play an active role in power system operation and even benefit financially from this opportunity. However, new operational challenges have been introduced due to the intrinsic characteristics of DERs such as intermittency of renewable resources, distributed nature of these resources, variety of DERs technologies and human-in-the-loop effect. Demand response (DR) is one of DERs and is highly influenced by human-in-the-loop effect. A data-driven based analysis is implemented to analyze and reveal the customers price responsiveness, and human-in-the-loop effect. The results confirm the critical impact of demographic characteristics of customers on their interaction with smart grid and their quality of service (QoS). The proposed framework is also applicable to other types of DERs. A chance-constraint based second-order-cone programming AC optimal power flow (SOCP-ACOPF) is utilized to dispatch DERs in distribution grid with knowing customers price responsiveness and energy output distribution. The simulation shows that the reliability of distribution gird can be improved by using chance-constraint.

In recent years, with the increasing penetration of solar generation, the uncertainty and variability of the power system generation also have increased. Power systems always require a balance between generation and load. The generation of the conventional generators must be scheduled to meet the total net load of the system with the variability and uncertainty of the solar resources integrated. The ability to match generation to load requires certain flexibility of the conventional generation units as well as a flexible transmission network to deliver the power. In this work, given the generation flexibility primarily reflected in the ramping rates, as well as the minimum and maximum output of the generation units, the transmission network flexibility is assessed using the metric developed in this work.

The main topic of this thesis is the examination of the transmission system flexibility using time series power flows (TSPFs). First, a TSPFs program is developed considering the economic dispatch of all the generating stations, as well as the available ramping rate of each generating unit. The time series power flow spans a period of 24 hours with 5-minute time interval and hence includes 288 power flow snapshots. Every power flow snapshot is created based on the power system topology and the previous system state. These power flow snapshots are referred to as the base case power flow below.

Sensitivity analysis is then conducted by using the TSPFs program as a primary tool, by fixing all but one of the system changes which include: solar penetration, wires to wires interconnection, expected retirements of coal units and expected participation in the energy

imbalance market. The impact of each individual change can be evaluated by the metric developed in the following chapters.

The Java-Bali power system is the biggest power system in Indonesia. On September 5th, 2018 at 11:26 AM, a region in the East Java-Bali subsystem suffered a blackout due to a single line to ground fault that disrupted the stability of the interconnected system and caused cascaded tripping.

This thesis presents the results of an evaluation of the dynamic performance of the East Java-Bali subsystem. It involves the static and dynamic simulations of the sequence of events that led to the East Java Bali subsystem blackout, especially the impact of the loss of a set of 500 kV transmission lines, which in reality was suspected to be the main cause.

The basic calculations related to power system state and familiarization with PSS/E (a commercial power system analysis software package) are first demonstrated. A simple 3-bus system test is taken as an example. The steady state characteristics of the active and reactive power injection, voltage and phase angle are calculated manually and compared to the PSS/E simulation results. As for the dynamic characteristics, short circuit current, electrical and mechanical power, rotor angle, and fault clearing time are determined by observing the plots of the simulation results. Based on understanding of the PSS/E modeling and simulation, the configuration, generation, and loading of the simplified East Java-Bali subsystem is evaluated. The generators (including the excitation system and governor) and transmission lines parameters are updated, as the reference model for the study. The model is validated by the actual data (active power flow) before the fault. Single line to ground fault and loss of generation disturbances were simulated to observe the stability of the system.

The analysis of the blackout is conducted through the simulation results based on all relevant documentation (such as fault report and sequence of events). With respect to the sequence of events (a single line to ground fault on the 500kV transmission lines, overload on 150kV transmission lines and tripping of power plants), several simulations of the East Java-Bali subsystem operations provided in the official blackout report are evaluated. Finally, the undervoltage load shedding strategy is evaluated and proposed as a solution to mitigate the blackout in the East Java-Bali subsystem.

The simulations reveal some interesting results regarding the operational characteristics of the East Java-Bali subsystem before the disturbances and during the cascaded tripping.

Reliable and secure operation of bulk power transmission system components is an important aspect of electric power engineering. Component failures in a transmission network can lead to serious consequences and impact system reliability. The operational health of the transmission assets plays a crucial role in determining the reliability of an electric grid. To achieve this goal, scheduled maintenance of bulk power system components is an important activity to secure the transmission system against unanticipated events. This thesis identifies critical transmission elements in a 500 kV transmission network utilizing a ranking strategy.

The impact of the failure of transmission assets operated by a major utility company in the Southwest United States on its power system network is studied. A methodology is used to quantify the impact and subsequently rank transmission assets in decreasing order of their criticality. The analysis is carried out on the power system network using a node breaker model and steady state analysis. The light load case of spring 2019, peak load case of summer 2023 and two intermediate load cases have been considered for the ranking. The contingency simulations and power flow studies have been carried out using a commercial power flow study software package, Positive Sequence Load Flow (PSLF). The results obtained from PSLF are analyzed using Matlab to obtain the desired ranking. The ranked list of transmission assets will enable asset managers to identify the assets that have the most significant impact on the overall power system network performance. Therefore, investment and maintenance decisions can be made effectively. A conclusion along with a recommendation for future work is also provided in the thesis.