Bulk Power Systems Emergency Control Based on Machine Learning Algorithms and Phasor Measurement Units Data: A State-of-the-Art Review

Abstract

Modern electrical power systems are characterized by a high rate of transient processes, the use of digital monitoring and control systems, and the accumulation of a large amount of technological information. The active integration of renewable energy sources contributes to reducing the inertia of power systems and changing the nature of transient processes. As a result, the effectiveness of emergency control systems decreases. Traditional emergency control systems operate based on the numerical analysis of power system dynamic models. This allows for finding the optimal set of preventive commands (solutions) in the form of disconnections of generating units, consumers, transmission lines, and other primary grid equipment. Thus, the steady-state or transient stability of a power system is provided. After the active integration of renewable sources into power systems, traditional emergency control algorithms became ineffective due to the time delay in finding the optimal set of control actions. Currently, machine learning algorithms are being developed that provide high performance and adaptability. This paper contains a meta-analysis of modern emergency control algorithms for power systems based on machine learning and synchronized phasor measurement data. It describes algorithms for determining disturbances in the power system, selecting control actions to maintain transient and steady-state stability, stability in voltage level, and limiting frequency. This study examines 53 studies piled on the development of a methodology for analyzing the stability of power systems based on ML algorithms. The analysis of the research is carried out in terms of accuracy, computational latency, and data used in training and testing. The most frequently used textual mathematical models of power systems are determined, and the most suitable ML algorithms for use in the operational control circuit of power systems in real time are determined. This paper also provides an analysis of the advantages and disadvantages of existing algorithms, as well as identifies areas for further research.

1. Introduction

Emergency control of electrical power systems (EC EPS) is a crucial element in ensuring reliable power supply to consumers and promptly preventing the consequences of emerging situations [1]. EC is classified into two types: local and centralized. Local EC is employed for protecting a small energy district or individual equipment from overcurrent, unacceptable voltage oscillations, loss of transient (TS) or static small signal (SSS) in the generating equipment of a selected power plant [2], and EPS separation [3] in case of stability loss, with frequency-based balancing of the allocated power district [4]. Local EC algorithms are widely used in practice and have demonstrated high efficiency throughout their operational period. However, they often overlook the impact of control actions (CA) on the EPS. EPS is a connected structure characterized by dynamic properties with a high degree of uncertainty during transients [5]. In such conditions, ensuring stability in a separate EPS fragment can lead to a loss of stability over remote connections. To prevent such phenomena, centralized emergency control is utilized, providing stability for the EPS [6]. Modern EPS is undergoing changes related to the basic principles of building EC systems, driven by several factors:

• Significant digitalization of the production, transmission, and distribution processes of electricity [7];
• Implementation of digital monitoring and management systems of EPS [8];
• Integration of a significant number of renewable energy sources (RES) [9,10,11];
• Accumulation of a significant amount of data describing the transient processes in EPS [12];
• Installation of flexible AC transmission devices (FACTS) [13];
• Using power storage devices [14];
• Tightening of the rules for the operation of the electricity market, which leads to an increase in active power flows [15].

Such changes directly affect the accuracy and reliability of the specification as follows [16,17]:

• The integration of RES in EPS leads to an increase in the rate of transients due to a decrease in the total inertia of EPS;
• An increase in the probability of loss of ESP stability due to an increase in active power flows through the elements of the electrical network;
• The high level of digitalization makes it necessary to consider the digital security of critical EPS sustainability infrastructures.

Currently, there are instances of inefficient operation of EC EPS due to insufficient speed and adaptability, as well as vulnerability to cyber-attacks [18,19,20,21]. On the other hand, the availability of large amounts of data and a developed infrastructure capable of receiving telemetry information about the status of EPS, transmitting commands for the implementation of CAs from a centralized control center, along with numerous systems for the mathematical analysis of processes occurring in EPS, enable the implementation of EC algorithms based on machine learning (ML) [22,23,24]. These algorithms are characterized by significant speed [25], adaptability, and self-learning capabilities. Unlike deterministic approaches to analyzing the mathematical model of EPS based on numerical methods, often used to determine the optimal volume of CA, ML algorithms operate with minimal time delay.

ML algorithms are used to solve the following tasks:

• Assessment of the technical condition of equipment [26,27];
• Determination of the disturbance type [28,29];
• EPS stability Analysis [30,31,32,33];
• Forecasting EPS load schedules [34,35,36,37,38];
• Forecasting RES power generation [39,40,41,42,43].

However, due to the high complexity of implementation, the use of ML for the EC is quite limited. From the point of view of ensuring EPS stability, ML algorithms allow for a certain set of features to calculate the optimal set of effects aimed at changing the state of the electrical network elements or changing the settings of control devices to preserve TS, SSS, and voltage stability. The main feature of using ML algorithms for EPS EC is related to the preparation and processing of data for training and test samples. The accuracy and adaptability of the trained ML algorithm directly depend on the quality, completeness, and representativeness of the data. Therefore, special attention is paid to the issue of the data used in the examined studies.

The objective of this paper is a meta-analysis of research focused on developing Emergency Control of Electrical Power Systems (EC EPS) algorithms using ML methods and data from synchronized phasor measurement units (PMUs). It examines the use of ML algorithms for selecting CAs to maintain TS, SSS, and voltage stability, and to ensure required frequency levels. Additionally, infrastructure solutions that enable the implementation of ML algorithms in real EPS are also considered. The scientific novelty of this work lies in identifying the advantages and disadvantages of existing EC EPS algorithms based on ML methods, as well as pinpointing directions for future research. The motivation of this article is to fill the gap in the systematization and analysis of research aimed at developing EPS EC methods to maintain TS, SSS, and voltage stability, and to ensure the required frequency level based on ML algorithms.

2. Existing Algorithms for Emergency Control of Power Systems

Figure 1 presents a general diagram of the Emergency Control of Electrical Power Systems (EC EPS) algorithm with settings defined in the offline mode [44]. To gather information on the flow of active and reactive powers, and the current loads of electrical network elements, as well as frequency and voltage values, current and voltage transformers are utilized. Measured continuous values of EPS performance parameters are converted into discrete, time-synchronized vector representations in Phasor Measurement Unit (PMU) devices [45]. Subsequently, these discrete synchrophasor values are fed into the disturbance detector block, which, based on the type of incident, sends a signal to the settings block, the structure of which is detailed in Figure 2.

The CAs settings block is a matrix structure that establishes the correlation between the type of disturbance and the necessary volume of CA to ensure the stability of the protected electrical power systems (EPS). Figure 2 shows the existing EPS EC implementation model. In this model, signals of the beginning of an accident (Fault 1, Fault 2, Fault 3) are transmitted via communication channels from substations or power plants to a single control center. The signal of the beginning of an accident is understood as a signal of disconnection of the power line, SG, and operation of the relay protection device. According to a predetermined logic, signals for the implementation of CAs are issued for each considered signal of the beginning of an accident. The signal data are also transmitted via communication channels at a power plant or substation. In Figure 2, potential disturbances are listed vertically, which may include emergency shutdowns of power lines and short circuits. The available EC EPS options are displayed horizontally: shutdown of synchronous generators, load shedding, unloading of synchronous generators, etc. To configure the EC EPS algorithm, it is essential to determine the linkage between the detected disturbance and the required control actions. In Figure 2, this correlation is indicated by color. The figure also exemplifies how Fault 1 is addressed by transmitting signals to implement CA 1 and CA 2. In Figure 2, the non-zero cells of the matrix are highlighted, defining the correspondence between Fault and CA.

The primary distinction between the Emergency Control of Electrical Power Systems (EC EPS) algorithm in offline and online modes lies in the approach to creating a correspondence matrix (implemented in the control action settings block) that links the type of disturbance with the required volume of CAs. In the offline mode, the creation of a correspondence matrix is typically conducted once a year for a pre-prepared mathematical model of the EPS, and the calculation of optimal CA volumes is manually performed by specialists managing the operational control of the protected EPS. A significant drawback of this approach in tuning the EC EPS algorithm is its limited adaptability and accuracy. Additionally, the parameters of the mathematical model of the EPS used may vary considerably from the actual parameters due to the impact of external factors on the elements of the electrical network [46,47]:

• Heating of conductors affects active resistance;
• Fogs and atmospheric pollution affect the level of corona on the surfaces of conductors;
• Insulation aging affects active and reactive resistance;
• Wind speed affects the carrying capacity of the conductor.

To increase the accuracy and adaptability of EC EPS, an algorithm is used in which the formation of a correspondence matrix between the type of disturbance and the volume of CA is performed in the “online” mode. The structure of the EC EPS algorithm in the online mode is shown in Figure 3. In Figure 3, the blocks involved in the cyclic calculation are highlighted with a fill.

The primary distinction between the Emergency Control of Electrical Power Systems (EC EPS) algorithm in the “online” and “offline” modes is the calculation of the optimal volume of CA with cyclicality. Figure 3 illustrates a segment of the cyclic updating of telemetry data and the updating of the correspondence matrix between the type of disturbance and the volume of CA. This method of implementing the EC EPS algorithm allows for the consideration of external factors, significantly enhancing accuracy and adaptability. A major challenge of the EC EPS algorithm in the “online” mode is the predetermined list of accidents for which EPS stability is ensured through the implementation of CA. To overcome this limitation, a CA selection algorithm can be employed at the rate of the transient process (post fault) [48]. This algorithm is fully adaptive, utilizing information about the current disturbance and current EPS parameters.

Table 1. Analysis of EC EPS algorithms.

Algorithm Type	Advantages	Disadvantages
Offline	High reliability of CA implement	Low adaptability and accuracy
Online	Increased adaptability caused by cyclically updating information about EPS parameters	There is a possibility of an error occurring at one of the stages of the cyclic calculation, only a predetermined list of considered disturbances is considered
Post fault	Increased adaptability by considering any disturbance and current EPS parameters	High-performance requirements for the algorithm for selecting the optimal CA set

provides an analysis of existing EC EPS algorithms.

All existing Emergency Control of Electrical Power Systems (EC EPS) algorithms exhibit specific advantages and disadvantages. The highest adaptability is found in algorithms based on the “post fault” principle, while the lowest is in those implemented on the “offline” principle. The most stringent requirements for algorithm speed in EC EPS modes are associated with the “post fault” principle. Currently, to enhance the speed of algorithms operating on the “online” principle, parallel computing systems are actively utilized. Figure 4 depicts a diagram of the implementation of parallel calculations in the CA algorithm of EC EPS following the “online” principle. For parallel computations, two types of servers are employed: an operational server and a calculation servers. The operational server generates a mathematical model of the protected EPS, links the objects of the calculation model to measurements, assesses the state, maintains a retrospective of calculations, and transfers the calculation results to remote emergency control controllers that execute the logic of the matrix correlation between disturbances and the calculated optimal set of CAs. Calculations are parallelized based on the disturbances considered. A management agent, which monitors the availability and load of the calculation servers, is used to distribute emergency disturbances among them. The emergency disturbance number serves as an input parameter for the calculation server; after the computation, the control agent relays the results of determining the optimal CA to the operational server [49,50].

Despite the advanced system for parallelizing calculations of CAs impacts, its efficiency is constrained by the performance of the algorithm used for determining the optimal EPS control law. Currently, these algorithms are predominantly implemented based on the deterministic principle of analyzing a system of algebraic differential equations that describe the dynamic model of the protected EPS. The efficacy of this approach is considerably limited. Moreover, the parallelization system itself is a significant bottleneck in the EC EPS chain. Although the operational server and management agent are configured as a failover cluster, there still exists a non-zero probability of system failure, which could result in a loss of EPS protection capabilities.

Consequently, implementing Emergency Control of Electrical Power Systems (EC EPS) algorithms based on the “post fault” principle using deterministic approaches proves challenging. Currently, the most appropriate method for implementing algorithms on the “post fault” principle is the utilization of ML algorithms. These algorithms offer high performance without the necessity of employing a parallel calculation system within the EPS operational control loop. ML algorithms are capable of determining the type of disturbance, selecting the optimal set of control actions (CAs) to ensure SSS, TS, voltage stability, and maintaining acceptable voltage levels and frequency [51].

Traditional EPS EC methods are based on the analysis of algebraic differential systems of equations describing the dynamic EPS model. This approach has several disadvantages associated with significant time delays necessary for the numerical solution of differential equations, considering EPS elements with unknown replacement circuit partners. On the other hand, ML algorithms have high performance, and the ability to implement the principle of “model-free control”, in which only measurements are used for the operation of the algorithm.

3. Determination of the Disturbance Type

The task of identifying the type of disturbance in EPS falls under the domain of digital signal processing (DSP) of instantaneous current and voltage values [52]. Common disturbances monitored in EPS Emergency Control (EC) algorithms include outages of electrical network elements (such as lines, transformers, and bus systems), short circuits (SC), and outages of loads and synchronous generators (SG). Figure 5 illustrates the classification of disturbances in EPS.

To determine the fault detection type (FD), methods that incorporate the following ML algorithms are used:

• Support vector machine (SVM) [53,54,55];
• Multi-class support vector machine s (MMC-SVM) [56];
• K-nearest neighbors (KNN) [57];
• Probabilistic neural network (PNN) [58];
• Artificial neural networks (ANN) [59];
• Convolutional adversarial neural network (CANN) [60,61];
• Chebyshev neural network (ChNN) [62];
• Decision trees (DT) [63];
• Rule-based decision tree (RBDT) [64];
• Fuzzy logic (FL) [65];
• Deep learning (DL) [66,67];
• Pattern recognition (PR) [68].

The following DSP methods are used in combination with ML algorithms:

• Discrete wavelet transform (DWT) [56];
• Discrete Fourier transform (DFT) [63];
• Hilbert–Huang transform (HHT) [66].

Figure 6 shows the main ML algorithms to determine the type of disturbance in the EPS.

In the study in [53], the SVM algorithm is employed to identify and classify disturbances in EPS. This algorithm analyzes the instantaneous values of currents and voltages over a sliding window with a duration of a quarter of a power frequency period. Initial data for the analysis are derived from measurements in both time and frequency domains at the end of the power line under examination. A dataset consisting of 25,168 disturbances was utilized to train and test the model. The signal sampling rate was set to 32 samples per industrial frequency cycle. The trained model achieved an average accuracy of 99.89%. The study also conducts a comparative analysis of SVM models with other models such as DT, random forest (RF), KNN, and ANN.

The authors in [54] introduced a method for identifying disturbance types for power quality analysis using the SVM combined with the fusion of time domain descriptors (FTDD) algorithm for feature selection. The classification criteria include the normal mode, transient process, harmonic distortion, voltage sag, voltage surge, and voltage interruption. The methodology was tested with real data from a power quality analyzer based on the Arduino controller. The dataset comprised 500 scenarios, equally divided between training and testing. The model’s accuracy was compared with standard DSP algorithms at various signal noise levels, achieving 99.11% accuracy with SVM + FTDD across all noise levels.

In [55], a methodology using SVM combined with DWT was developed to identify the type and location of disturbances on overhead and cable power lines. Tested on a 230 kV power transmission line model, the methodology considered 2448 scenarios of various disturbance types and locations. The average accuracy for disturbance identification was 98.8%, with a location determination error of less than 3%.

Research [56] explored the use of the MMC-SVM algorithm applied to root mean square voltage measurements at EPS nodes. The algorithm’s first stage identifies the occurrence of a disturbance, while the second stage recognizes its type. Data included transient processes with SC on power lines, load nodes, and substations without power take-offs. White noise ranging from 10 to 90 dB was added to the data for testing. Classification accuracy without noise was 97.16%, but adding noise reduced accuracy by an average of 7%. The study in [57] proposed a method for identifying disturbance types in distribution networks, notable for requiring minimal data in the training set. The method was compared to ANN and DT algorithms, achieving an accuracy of 89.15%, significantly higher than the 24.09% and 43.43% accuracies of the ANN and DT algorithms, respectively. The study also examined the impact of signal sampling frequency on classification accuracy, finding that increasing the frequency from 960 Hz to 3840 Hz improved accuracy by 16.86%.

The authors of [58] proposed a method for identifying the type of disturbance based on the PNN algorithm using DWT to the original signals. The technique proposed by the authors consists of converting a signal from the time domain to the frequency domain and then using the frequency spectrum as input data to the PNN. The structure of a PNN network is as follows: input layer, sample layer, summation layer, and output layer. To test the proposed methodology, the simplest EPS model is used, consisting of two catches, SG, load, and power line. Different types of SCs in different parts of the transmission line were considered disturbances. The authors claim that the accuracy of the proposed method reaches 100%.

An example of using ANN to determine the type of disturbance is given in [59]. Testing was performed on the simplest EPS model, consisting of two nodes and a power line. The ANN structure consists of three layers; a more detailed description of the model is not provided. The initial data sampling frequency was 10 kHz; the authors determined that reducing the sampling frequency to 3.4 kHz does not lead to a decrease in classification accuracy but helps to increase performance.

The study in [60] considers the analysis of disturbances that arose in the EPS Association of Southeast Asian Nations. To identify the type of disturbance, the authors used the CANN algorithm. The work does not provide an analysis of accuracy and time delays. The study in [61] also uses the CANN algorithm to identify the type of disturbance in distribution networks. Information provided by Duke Energy was used as input data. The following features of the initial dataset were used: power line number, weather conditions, season, hour of day, number of the damaged phase, and number of protection devices in operation. The accuracy of the CANN model was compared with the linear regression (LR) model. The average accuracy of the CANN algorithm was 85.94%. The authors consider the problems of lack of information and imbalance of data on disturbances in distribution networks. In [62], the authors considered the problem of identifying the type of disturbance on a power transmission line with longitudinal compensation based on the ChNN algorithm. The ChNN model has a flat single-layer structure, the Chebyshev polygon is used as a functional extension. ChNN is trained using the recursive least squares method. Instantaneous phase current signals with a sampling frequency of 4 kHz are used as input data for ChNN. The proposed model is compared with the SVM, FL, and DT algorithms. The average accuracy of the proposed model was 99.43%, which exceeds all considered models.

To identify the type of disturbance, the study in [63] uses the DT algorithm. Current and voltage signals from both ends of the protected power line are used as initial data. The DFT algorithm is used to extract features from the source data. Testing was performed considering varying degrees of noise in the source data. When testing the proposed algorithm for identifying the type of disturbance, 100% accuracy was achieved. In [64], the RBDT algorithm was used to identify the type of disturbance. To generate the initial data for identification, the HHT algorithm was used, which allows the use of 1st and 2nd modes of the signal. The authors claim that the proposed algorithm is effective in detecting and classifying various types of EPS faults. The work [65] proposed a combined approach of using FL and DWT to identify the type of damage. For testing, a two-node EPS model with a power line equipped with a longitudinal compensation device was used. The proposed approach has been tested for disturbances of various types.

To identify faults in electrical distribution networks with an isolated neutral, the work [66] uses a combined method that combines the DL and HHT algorithms. Current and voltage signals are used as initial data. To identify the type of disturbance, a time-frequency image is used, which is supplied as input to the DL algorithm. The ANN used consists of five layers with a sigmoid activation function in each layer. Testing was performed on digital and physical EPS models. The average accuracy of the proposed model was 99.99%. In [67], the DL algorithm is also used to identify the top of the disturbance. Testing was performed on the simplest EPS model, consisting of two nodes and one branch. To form a matrix of signal energy coefficients, the DWT algorithm is used. It is important to note that DL algorithms use a graphical representation of the signals under consideration as input information. The accuracy when testing the algorithm was 98.00%. The authors in [68] use the PR method to identify the type of disturbance. Testing was performed on IEEE9 and IEEE39 EPS models. The measurements obtained from the PMU are used as input data.

A total of 16 articles were reviewed on the topic of identifying the type of disturbance based on ML algorithms.

Table 2. Analysis of the reviewed articles devoted to identifying the type of disturbance in EPS.

Characteristic	Articles, %
An analysis of the algorithm’s performance is given	25.0
A description of the data sample is provided	62.5
A mathematical model of the EPS is used	75.0
The data of real EPS is used	25.0
An analysis of noise and inaccuracies that impact the data is performed	25.0

provides a numerical analysis of the reviewed works.

ML-based disturbance type identification in EPS is a popular research topic, as supported by a significant number of papers and approaches. A general promising direction for the development of the topic of disturbance type identification is real-time testing with determination of algorithm delays and possible optimization, allowing the use of disturbance type identification algorithms for the EPS EC problem.

4. Ensuring Small Signal Stability

One of the main tasks of ensuring the reliability of big EPS (BEPS) is monitoring and managing the magnitude of active power flows along long power transmission lines connecting remote energy regions. The determining phenomenon limiting exchange flows between EPSs is SSS. This problem is effectively solved by limiting active power flows and using EC devices. To provide ML-based SSS, the following algorithms are used:

• Extreme gradient boosting (XGBoost) [69];
• Long short-term memory networks (LSTM) [70];
• ANN [71];
• DL [72];
• KNN [73];
• Particle swarm optimization k-nearest neighbors (PSO-KNN) [74];
• Generalized regression neural network (GRNN) [75].

Figure 7 shows the considered ML algorithms to determine the CA for SSS preservation.

In [69], a method for selecting the optimal CA volume for preserving SSS based on the XGBoost algorithm and analysis of the total damping coefficient EPS is proposed. To select CA, the following system of equations are used:

$$\left\{\begin{array}{c}{\displaystyle \sum _{i=1}^n}∆P_i^2\to min\\ P_i^{min}\le P_{i0}+∆P_i\le P_i^{max}\\ k_0+∆k\ge k_{lim}\\ {\displaystyle \sum _{i=1}^n}{C}_i∆P_i=∆k\end{array},\right.$$

(1)

where ΔP_i—change in active power of SG with number i, $P_i^{min}$—minimum active power of SG with number i, $P_i^{max}$—maximum active power of SG with number i, P_i0—initial active power of SG with number i, k₀—damping coefficient of SG with number i, Δk—change in the damping coefficient EPS when implementing CA, k_lim—minimum value of the damping coefficient EPS, and C_i—influence coefficient of changes in power SG with the number i on the total damping coefficient EPS.

The proposed methodology was tested on IEEE9 and IEEE39 models. To form a data sample, numerical modeling was used with a change in the total load of the EPS model from 70% to 130% with a corresponding change in the active powers SG. Minimum values for total EPS damping coefficients were also determined. During testing, the accuracy of the proposed method was determined to be 99.2% and 97.1% for the IEEE9 and IEEE39 models, respectively.

The authors of the study in [70] proposed a universal technique for selecting CAs at a real-time pace for preserving SSS and TS based on the LSTM algorithm. The scientific novelty of the proposed method is the real-time prediction of the damping coefficient EPS, and local and intersystem oscillatory modes. The following EPS mathematical models were used for testing: IEEE39, IEEE68, IEEE145. The delay of the proposed method was less than 1.5 periods of power frequency.

The study in [71] aims to develop an SSS EPS provisioning system based on the ANN algorithm. The system considers single outages of power lines, load changes, and single outages of generators. A numerical example shows the high efficiency of the proposed system.

In [72], the authors propose an EPS stabilizer structure based on the DL algorithm. The proposed stabilizer is designed to provide SSS in the presence of wind generation in the EPS. The EPS mathematical model was used for testing. It consists of a wind generation source and a power transmission line.

The paper in [73] presents a technique for estimating SSS based on the KNN algorithm. A platform is proposed for SSS analysis, the block diagram of which is shown in Figure 8.

The authors of [73] indicate the CA selection as a direction for platform development to ensure the required frequency level in isolated EPS. For the numerical study, the IEEE59 model was used, in which 8760 scenarios were calculated with subsequence CAs selection. The work shows the effectiveness of the proposed method.

The study in [74] served as a continuation of [73], focusing on developing an accelerated method for analyzing SSS. For this purpose, the authors utilized a platform they developed, which encompasses scenario calculation, modeling of EPS operations within the electricity market, and SSS analysis. A modified version of the IEEE59 model was employed for numerical simulations. The PSO-KNN algorithm was applied for feature selection and clustering.

In [75], the analysis of SSS introduces the concept of a reference point (RP), which is used to identify the region of EPS operating modes within a given SSS margin. A specially developed clustering algorithm, based on the analysis of negative Euclidean distance, is employed to delineate this area. For real-time SSS analysis, calculations are conducted on a sliding window, involving operations such as determining the location of the operating point within the stability region. When the operating point exceeds the stability region, the GRNN algorithm is utilized to identify the optimal set of control actions (CAs). The method was tested on the IEEE9 and IEEE118 models, demonstrating an accuracy of over 94%.

In Section 4, articles on CA selection methods for providing SSS EPS based on ML algorithms were reviewed.

Table 3. Analysis of reviewed articles devoted to the definition of CA for maintaining SSS EPS.

Characteristic	Articles, %
An analysis of the algorithm’s performance is given	28.5
A description of the data sample is provided	57.1
A mathematical model of the EPS is used	100.0
The data of real EPS is used	0.0
An analysis of noise and inaccuracies that impact the data is performed	0.0

shows a numerical analysis of the considered works.

Testing of the considered algorithms aimed at determining CA for saving SSS EPS was performed on mathematical models IEEE9, IEEE39, IEEE59, IEEE68, IEEE118, and IEEE145. For these models, problems of outliers and noise in the data are not considered. These problems are typical for real EPS. Also, in the works under review, insufficient attention is paid to the issue of using the developed algorithms in real time.

5. Providing Transient Stability

The increase in RES number, which is observed in modern EPS, significantly complicates the issue of maintaining TS due to a significant increase in the speed of transient processes. Therefore, the issue of maintaining TS has received more discussion in research compared to SSS.

To provide TS based on ML methods, the following algorithms are used:

• ANN [76];
• FL [77];
• Mixed-integer linear programming (MILP) [78];
• Deep belief network (DBN) [79];
• Core vector machine (CVM) [80];
• Convolutional neural network (CNN) [81];
• Stacked denoising autoencoder (SDAE) [82];
• Twin convolutional support vector machine (TCSVM) [83];
• Extreme learning machine (ELM) [84];
• XGBoost [85];
• Mahalanobis kernel regression (MKR) [86].

Figure 9 shows the considered ML algorithms determining the CA for TS.

A notable characteristic of modern BEPS is the extensive use of converters and high-voltage converter equipment. In [76], the control of converters is suggested as a means to maintain EPS TS. For the numerical experiment, a mathematical model of an EPS comprising two nodes, a transmission line, and an inverter was utilized. The ANN algorithm was employed for selecting the CAs.

The authors in [77] proposed a CA selection model for TS EPS based on the FL algorithm, which is described by the following expression:

$$R_m:IF x_1 is A_1^m and x_2 is\dots and x_n is A_n^m THEN y is B^m,$$

(2)

where x_i—input variables, $A_1^m$ and $B_1^m$—fuzzy sets that constitute R_m, m = 1, 2 … M.

For the numerical experiment, a Real-Time Digital Simulator (RTDS) was used, in which a two-machine EPS model was implemented. The accuracy of the proposed algorithm was more than 92%.

To select a CA for TS EPS, the authors of [78] proposed an algorithm based on MILP. Testing was performed using mathematical models IEEE9 and the 74-bus Nordic test system. For the IEEE9 model, a data sample of 2000 scenarios was generated, and the average CA selection accuracy was 99.48%, with a time delay of 450 ms. For the 74-bus Nordic test system model, the average accuracy was 97.42%.

The authors of the study in [79] proposed a technique for selecting CAs for TS EPS based on DBN. The methodology consists of two parts: offline and online. For the offline part, synthetic transient modeling, feature selection, data preparation, and DBN model training are performed. In the online part, real data entering the input layer of the trained model is used for analysis. The paper compares the proposed method with CNN, KNN, RF, and multilayer perceptron (MLP) algorithms. As a result of the comparison, it was found that the DBN algorithm allows us to obtain the smallest CA selection error among those considered.

One of the SVM algorithm limitations is the assumption of the possibility of linear separation of classes in the data sample. To analyze classes with nonlinear separation, the work [80] uses the CVM algorithm, which allows the use of a nonlinear kernel function. To synthesize the data sample, the authors propose to use three sources: historical data, data from PMU, and data from traditional measurement sources. The IEEE39 model is used for testing. To form a data sample, 10 load levels were considered in each node: 80–125%. At the same time, the total data sample size was 5.310. The average accuracy for the IEEE39 model was 93.04%. Also, the proposed methodology was tested on real data. Data from the North–Central–East China (NCE) and Eastern Interconnection (EI) EPS of the United States were reviewed. The characteristics of the data samples and the obtained CA selection accuracies for real data are given in Table 8.

To select the optimal set of CAs to preserve TS, the authors in [81] used the CNN algorithm. The authors proposed a unique graphical technique for describing the transient process, which is resistant to noise in the data, which consists of applying the DFT algorithm to the signals of voltage, phase, frequency change, and active and reactive power of each SG of the EPS under consideration. As a result, the authors obtain a description of the transient process in the form of a graphical representation.

Further, the resulting graphical image is used as input to the CNN algorithm. Testing of the proposed methodology is performed on the IEEE39 model and real data obtained from the Guangdong Power Grid (GPG). Characteristics of data samples and the obtained CA selection accuracies for real data are given in Table 8.

In [82], the SDAE algorithm is used to analyze TS. The proposed TS analysis technique consists of the following steps:

• Creating an autoencoder with noise reduction;
• The use of adaptive synthetic sampling;
• The synthesized data are decoded into the original space;
• Using a classifier based on the SDAE algorithm.

Testing was performed on IEEE39 and South Carolina 500-Bus System models. For the IEEE39 model, the accuracy was 98.78%; for the South Carolina 500-Bus System model, an accuracy of 98.10% was achieved. The authors highlight the use of data from PMU and consider changes in EPS topology as directions for future research.

In [83], the TCSVM algorithm is used to analyze TS. To analyze TS, the following structure of individual modules is used:

• Data generation module;
• Feature selection module;
• Module for predicting the trajectory of the transition process.

For the numerical experiment, the following EPS models were used: the Brazilian 7-Bus equivalent model, IEEE68, the two-area system, SAVNW, and IEEE24.

The authors in [84] used the ELM algorithm in combination with an equivalent one-machine infinite bus (OMIB) to analyze TS. The proposed method is based on data application only, without using the EPS model. For testing, a modified IEEE39 model was used with the addition of wind generation at six nodes of the EPS model. The following features were used for the ELM algorithm: angular position and speed of SG rotors, output active power SG, angular position and speed of wind generator rotors, voltage amplitudes, and phases in the nodes of the EPS model.

In [85], the XGBoost algorithm is used in combination with Factorization Machine (FM) to analyze TS. The paper presents the concept of automatic feature selection. For TS analysis, a combination of XGBoost and FM algorithms is used as follows: the data obtained from the PMU is the input to the XGBoost algorithm, which transforms the time series into a space–time matrix, which is the input to the FM algorithm. IEEE39, IEEE68, and IEEE140 models are used for testing. The authors compare the developed methodology with the RF, SVM, DT, and MLP algorithms. An analysis was made of the influence of noise in the source data on the reliability of the algorithm.

The authors of the study in [86] present a technique to analyze TS based on the MKR algorithm. The proposed technique is intended for use in real time. Testing is performed on the IEEE39 model.

In this section, 11 articles devoted to methods for selecting CAs to provide TS EPS based on ML algorithms were reviewed.

Table 4. Analysis of reviewed articles devoted to determining CA for TS EPS.

Characteristic	Articles, %
An analysis of the algorithm’s performance is given	45.4
A description of the data sample is provided	63.6
A mathematical model of the EPS is used	63.6
The data of real EPS is used	36.3
An analysis of noise and inaccuracies that impact the data is performed	27.2

shows a numerical analysis of the reviewed works.

A common problem with the considered algorithms is the insufficient elaboration of the use of data from real EPS, as well as testing the algorithms applied by using RTDS.

6. Providing Voltage Stability

Maintaining required voltage levels is a challenge common to both BEPS and isolated EPS. Oscillatory instabilities in voltage levels can trigger avalanche-like processes, leading to widespread consumer outages. To address this issue, various algorithms based on ML methods are employed to ensure VS:

• Deep reinforcement learning (DRL) [87];
• Emotional deep neural network (EDNN) [88];
• Back propagation neural network (BPNN) [89];
• ANN [90];
• Spatio-Temporal Graph Convolutional Network—Double Deep Q-Network (STGCN-DDQN) [91];
• DT [92];
• SVM [93];
• Artificial emotional reinforcement learning (AERL) [94].

Figure 10 shows the considered ML algorithms to determine the CA for VS.

In [87], the DRL algorithm is used to provide VS. An important feature of the DRL-based VS analysis technique proposed by the authors is the consideration of noise in the source data. For testing, the IEEE8500 model is used. The volume of the training sample was 500 modes. The average accuracy of the proposed method is 94.92%.

The main challenge of ANN is the ability to classify features with nonlinear separation. To overcome this problem, the authors of the study in [88] used the EDNN algorithm, which contains a different number of neurons in each layer, which makes it possible to model the nonlinear separation of classes in a data sample. The numerical experiment demonstrated the high accuracy of the proposed method.

In [89], the authors used the BPNN algorithm to save VS. The application of this algorithm is used to analyze the voltage stability region and select preventive EPS control measures. The proposed methodology was tested using the IEEE118 mathematical model. The data sample size was 5000 scenarios. The accuracy of the algorithm for the simulated data was 99.99%.

The authors of [90] considered the problem of preserving VS in isolated EPS. For the numerical experiment, the simplest IEEE4 model was used. For the developed model, an accuracy of 98.73% was shown.

In the study in [91], a DL algorithm is used to preserve VS. To determine the optimal CA volume, the authors use the voltage recovery criterion:

$$\left\{\begin{array}{c}{V}_1\le V\left(t\right)\le V_{max}, { t}_0\le t\le t_1\\ V_2\le V\left(t\right)\le V_{max}, { t}_1\le t\le t_2\\ V_3\le V\left(t\right)\le V_{max}, { t}_2\le t\le T_f\end{array}\right.,$$

(3)

where V₁, V₂, V₃, and V_max—voltage requirements at different stages of the transient process and the upper voltage limit accordingly, t₁ and t₂—time moments separating different stages of the transient process, and T_f—observation period. The Nordic Benchmark System model was used for testing; the data sample size was 18,000 scenarios. The study provides an analysis of the possibility of using the developed methodology in real time.

The authors of [92] used the DT algorithm to analyze VS; to reduce the size of the data sample, they used the multi-objective biogeography-based optimization (MOBBO) algorithm. Real Iranian EPS data were used for testing. The data sample size was 1898, the average accuracy was 95.40%.

The study in [93] proposed a methodology for estimating VS based on the SVM algorithm and data coming from the PMU. Testing was performed on the IEEE39 model. The developed VS analysis technique is aimed at application in real time. The accuracy of the proposed method reaches 99.99%.

The study in [94] presented a technique for analyzing VS based on the AERL algorithm. Optimization of stress levels in EPS is performed based on the minimization of the following objective function:

$$F=min\left({\mu }_1{\mathrm{C}}_{ds}^{\prime }+{\mu }_2{P}_{loss}+{\mu }_3{V}_d\right),$$

(4)

where µ₁, µ₂, µ₃—weight components of electricity production costs (${\mathrm{C}}_{ds}^{\prime }$), active power losses (P_loss), and steady voltage values (V_d).

The proposed methodology was tested using mathematical models IEEE57, IEEE118, and IEEE300. During numerical experiments, the authors demonstrated the high efficiency of the proposed method.

The section reviewed eight articles devoted to methods for selecting CAs to provide VS EPS based on ML algorithms.

Table 5. Analysis of reviewed articles devoted to determining CA for VS EPS.

Characteristic	Articles, %
An analysis of the algorithm’s performance is given	42.8
A description of the data sample is provided	50.0
A mathematical model of the EPS is used	71.4
The data of real EPS is used	28.5
An analysis of noise and inaccuracies that impact the data is performed	42.8

shows a numerical analysis of the reviewed works.

A common problem of the above-mentioned algorithms is insufficient consideration of using data from real EPS.

7. Providing Acceptable Frequency Levels

To solve the frequency control (FC) problem based on ML methods, the following algorithms are used:

• RF [95];
• DL [96];
• Convolutional neural network—long short-term memory (CNN-LSTM) [97];
• FL [98];
• CNN [99];
• SVM [100];
• Reinforcement learning (RL) [101];
• Particle swarm optimization (PSO) [102];
• Big bang big crunch (BBBC) [103];
• Genetic algorithm (GA) [104];
• Dueling Deep Q-Learning (DQN) [105].

Figure 11 shows the considered ML algorithms to determine CA for FC.

The authors of [95] used the RF algorithm to solve the FC problem. A high-voltage direct current (HVDC) system was used to control the frequency. As a test model, two fragments of real EPS connected with an HVDC device were used. The total data sample size was 200 scenarios. The average accuracy of the method was 98.91%. The study in [96] presented an FC technique based on the DL algorithm. The authors provided a detailed mathematical description of the FC problem. Testing was performed on the IEEE39 model. The work [97] provided a large study of the FC problem based on the CNN-LSTM algorithm. As input data for the CNN-LSTM algorithm the active powers of each SG, the active power imbalances of each SG, the voltage amplitude in each node, the voltage phase in each node, and the active power of each load were used. To use the CNN-LSTM algorithm, the input data were transformed into a weighted matrix structure, which is formed using tensor calculation. To reduce the number of features, the authors consistently use the following algorithms: Principal Component Analysis, Linear Discriminant Analysis, and t-distributed stochastic neighbor embedding (t-SNE). Testing was performed using IEEE39 and ACTIVSg500 models. To control the frequency in EPS, the authors of [98] proposed a controller based on the FL algorithm. A simplified dynamic model and an IEEE39 model were used as a test model. The authors suggested using the developed method in real time. In the study in [99], a CNN algorithm was used to control the frequency, which consisted of the following layers: input layer (representation of the input data as a tensor), convolutional layer, subsampling layer, and fully connected layer. The IEEE39 model was used for testing. In [100], the SVM algorithm was used to solve the FC problem. Testing was performed on IEEE39 and NRPG 246 models. In works [101,102,103,104,105], the RL, PSO, BBBC, GA, and DQN algorithms were used to solve the FC problem. In this section, 11 articles devoted to methods for selecting CAs for FC EPS based on ML algorithms were considered.

Table 6. Analysis of reviewed articles devoted to determining CA for maintaining TS EPS.

Characteristic	Articles, %
An analysis of the algorithm’s performance is given	27.2
A description of the data sample is provided	27.2
A mathematical model of the EPS is used	54.5
The data of real EPS is used	45.4
An analysis of noise and inaccuracies that impact the data is performed	27.2

provides a numerical analysis of the reviewed works.

A common problem with the above-considered algorithms is the insufficient study of using data from real EPS.

8. Directions for Future Research

Table 7. Directions for the development of EC EPS algorithms.

Direction	Tendency of Development
FD	Testing and adaptation of algorithms for real-time operation with determination of acceptable delays and sampling rates of source signals
SSS	Testing algorithms on data obtained from real EPS. Development of algorithms for processing noise and outliers in source data. Real-time testing.
TS	Testing algorithms on data obtained from real EPS. Using RTDS to determine the actual algorithm delay.
VS	Testing algorithms on data obtained from real EPS.
FC	Testing algorithms on data obtained from real EPS. Analysis of the possibility of use in real time.

provides directions for future research.

General directions for the development of EC EPS algorithms include increased attention to the use of data from PMU [105,106,107,108,109], the development of universal analysis algorithms SSS, TS, VS, and FC [110], the use of real-time modeling systems, development of adaptive local automation devices [111].

Table 8. EC EPS algorithm analyses.

Algorithm	Ref.	Data Source	Parameters	Accuracy, %	Merits	Drawbacks
	FD
SVM	[53]	IEEE9	Delay = 5 ms DataSetSize = 25,168	99.80	The model achieves high speed and maximum accuracy.	The method is intended only for analyzing disturbances on transmission lines. The work does not provide the obtained hyperparameters of the SVM model.
SVM + FTDD	[54]	Real data	Delay = 7.82 s DataSetSize =100	99.11	The FTDD algorithm is used to select features. The proposed technique is highly resistant to noise in the data.	Disturbances typical for analyzing the quality of electricity in low-voltage networks are considered. The application of the technique for high-voltage networks and disturbances in the form of SC has not been studied.
SVM + DWT	[55]	230-kV 60-Hz transmission line model	Delay = 60 ms DataSetSize = 2448	98.80	The following kernel functions were considered: linear, polynomial, Gaussian. The technique allows us to determine the type and location of the disturbance.	The method is intended only for analyzing disturbances on overhead and cable transmission lines.
MMC-SVM	[56]	IEEE13	DataSetSize = 756	97.16	High reliability and resistance to noise in the source data.	A simple radial EPS model is considered, but testing of the method on more complex models with ring connections is not provided. The resulting hyperparameters of the SVM algorithm are not given.
KNN	[57]	Real data	DataSetSize = 334	89.15	High reliability, and minimal training sample size.	No estimate of disturbance type classification delay.
PNN + DWT	[58]	Two buses one transmission line model	—	100.00	The paper presents an original approach to using DWT results as input data for PNN.	Testing was performed on the simplest EPS model, there is no description of the data sample, and there is no estimate of the method delay.
ANN	[59]	Two buses one transmission line model	DataSetSize = 49,500	82.79	The acceptable sampling rate of the source data has been determined.	The method is intended only for analyzing disturbances on transmission lines, sampling parameters are not given, and the delay of the method is not estimated.
CANN	[60]	Association of Southeast Asian Nations EPS data	DataSetSize = 12,890	—	Records of disturbances in real EPS have been used.	The paper does not provide an analysis of the disturbance classification accuracy, nor does it present an analysis of the CANN algorithm’s response time.
CANN	[61]	Duke Energy fault data	DataSetSize = 8376	85.94	Recordings of disturbances in real EPS were used.	An analysis of the time delay of the disturbance type identification algorithm is not provided.
ChNN	[62]	Two buses one transmission line model	DataSetSize = 57,600	99.44	An accurate and fast ChNN-based disturbance identification model is presented. A ChNN learning algorithm is proposed.	An analysis of the time delay of the disturbance type identification algorithm is not provided.
DT	[63]	Two buses one transmission line model	DataSetSize = 2000	100.00	An analysis of the identification accuracy dependence of disturbance type on the degree of noise in the original data was performed.	An analysis of the time delay of the disturbance type identification algorithm is not provided. Testing was performed on a simple model of a power transmission line.
RBDT+ HHT	[64]	Two buses one transmission line model	Delay = 20 ms	—	High efficiency in disturbance identification and classification.	The parameters of the training set and the obtained parameters of the DT algorithm are not described.
FL	[65]	Two buses one transmission line model	—	—	High efficiency in analyzing disturbances on transmission lines with longitudinal compensation.	The parameters of the training sample are not described, the performance of the algorithm is not assessed, and the accuracy of the algorithm is not assessed.
DL + HHT	[66]	Distributed EPS model	DataSetSize = 1752	99.99	High-efficiency, testing was performed on EPS physical and mathematical models.	The performance of the algorithm has not been assessed.
DL	[67]	Two buses one transmission line model	—	98.00	The developed algorithm is characterized by a low volume of the required training sample, adaptation to changes in the EPS topology, and resistance to noise and errors in the data.	The performance of the algorithm has not been assessed. Parameters of the training sample are not provided.
PR	[68]	1. IEEE9, 2. IEEE39	1. Delay = 32 ms 2. Delay = 96 ms	97.44	The proposed method is highly reliable and accurate.	Parameters of the training sample are not provided.
SSS
XGBoost	[69]	1. IEEE9, 2. IEEE39	DataSetSize = 6300	1. 99.2, 2. 97.1	The obtained values of the hyperparameters of the XGBoost algorithm are presented.	The method for determining the total damping coefficient EPS and its minimum value is not described. An analysis of the numerical delay values for CA selection is not provided.
LSTM	[70]	1. IEEE39, 2. IEEE68, 3. IEEE145	Delay = 30 ms	1. 99.45, 2. 98.40, 3. 98.46	A universal method for analyzing SSS and TS in real time is presented.	The procedure for selecting features for data sampling is not provided.
ANN	[71]	—	—	—	A methodology for analyzing SSS considering single outages of power lines, generators, and load changes is presented.	The procedure for selecting features for the data sample is not provided.
DL	[72]	Two buses one transmission line model	—	—	The issue of EPS control using a stabilizer is considered.	The procedure for selecting features for the data sample is not provided. There is no analysis of the numerical delay and accuracy of the method.
KNN	[73]	IEEE59	DataSetSize = 8760	—	A flexible system has been developed that allows modeling transient processes considering the rules of functioning of the electricity market.	No analysis of method accuracy and time delay is provided.
PSO-KNN	[74]	IEEE59	DataSetSize = 8760	—	A flexible technique for SSS analysis based on the PSO-KNN algorithm is presented.	An analysis of the method’s accuracy and time delay is not provided.
GRNN	[75]	IEEE118	Delay = 140 ms DataSetSize = 960	95.84	The paper introduces the concept of RP, which is intended to determine the area of sustainability of EPS.	No study is provided to determine the optimal sliding window size to effectively use the GRNN algorithm.
TS
ANN	[76]	Two buses one transmission line model	Delay = 40 ms DataSetSize = 10,736	—	The study provides a detailed methodology for data generation, processing, and use for CA selection based on the ANN algorithm.	Testing was performed on a simple EPS model.
FL	[77]	2-machine power system	—	92.00	Testing was performed using RTDS.	A simple EPS model was chosen for testing.
MILP	[78]	1. IEEE9, 2 74-bus Nordic test system	Delay = 450 ms DataSetSize = 2000	1. 94.8, 2. 97.42	The proposed method for selecting CAs for preserving TS is shown to be highly accurate.	Testing on real EPS data is not provided.
DBN	[79]	IEEE39	DataSetSize = 42,630	99.01	The accuracy of the proposed model is compared with the CNN, KNN, RF, and MLP algorithms.	The numerical delay of the model is not given, and testing on real data is not considered.
CVM	[80]	1. IEEE39, 2. NCE, 3. EI	1. DataSetSize = 5310 2. DataSetSize = 50,000 3. DataSetSize = 56,000	1. 93.04 2. 95.81 3. 99.50	Testing of the methodology on mathematical and real data is presented.	No numerical delay analysis is provided for determining optimal CAs.
CNN	[81]	1. IEEE39, 2. GPG	1. Delay = 12 ms 1. DataSetSize = 7200 2. Delay = 16 ms 2. DataSetSize = 18,000	1. 98.76 2. 89.35	Testing of the methodology on mathematical and real data is given, delays of algorithms are indicated, and an original method of presenting transient data in graphical form is given.	Real-time testing is not provided.
SDAE	[82]	1. IEEE39, 2. South Carolina 500-Bus System	1. DataSetSize = 5775 2. DataSetSize = 34,725	1. 98.78 2. 98.10	The TS analysis method is highly accurate.	An analysis of the possibility of using the algorithm in real time is not provided. Data from PMU and EPS topology changes are not considered.
TCSVM	[83]	1. Brazilian 7-Bus equivalent model, 2. IEEE68, 3. the two-area system SAVNW, 4. IEEE24	—	86.27	The TS analysis method is highly accurate.	Data samples are not described.
ELM	[84]	IEEE39	Delay = 0.082 s	89.4	The proposed method is highly efficient and fast.	An analysis of the operation of the proposed method on real data is not provided.
XGBoost + FM	[85]	1. IEEE39, 2. IEEE68, 3. IEEE140	—	1. 97.38 2. 97.93 3. 99.29	The proposed TS analysis method makes it possible to automatically select features, is resistant to noise, and has high accuracy.	An analysis of the operation of the proposed method on real data is not provided.
MKR	[86]	IEEE39	Delay = 3.66 ms DataSetSize = 80,000	98.68	The proposed TS analysis algorithm has high speed and accuracy.	The influence of noise in the source data on the accuracy of the algorithm is not shown.
VS
DRL	[87]	IEEE8500	DataSetSize = 500	94.92	The proposed method is resistant to noise in the source data and has high accuracy.	An analysis of the operation of the proposed method on real data is not provided. There is no study of the possibility of using the technique in real time.
EDNN	[88]	—	—	—	An original method for voltage control in EPS is presented.	An analysis of the operation of the proposed method on real data is not provided.
BPNN	[89]	IEEE118	DataSetSize = 5000	99.99	High accuracy of the trained algorithm.	An analysis of the operation of the proposed method on real data is not provided.
ANN	[90]	IEEE4	Delay = 0.302 s	98.73	The developed model has been shown to have high accuracy and acceptable performance.	Testing was performed on the simplest EPS model. There is no study of the possibility of using the technique in real time.
STGCN-DDQN	[91]	Real-world Swedish and Nordic power grid	DataSetSize = 18,000 Delay = 0.0079 s	—	The developed model is highly efficient and acceptable.	—
DT	[92]	Iranian 66 bus	DataSetSize = 1898	95.40	High efficiency.	There is no study of the possibility of using the technique in real time.
SVM	[93]	IEEE39	—	99.99	Possibility of application in real time.	An analysis of the operation of the proposed method on real data is not provided.
AERL	[94]	1. IEEE57, 2. IEEE118, 3. IEEE300	1. Delay = 21.68 s, 2. Delay = 67.20 s, 3. Delay = 248.64 s	—	High efficiency.	An analysis of the operation of the proposed method on real data is not provided. It is difficult to use the technique in real time.
FC
RF	[95]	EI&WECC	DataSetSize = 200	98.91	For FC the HVDC is used.	There is no study of the possibility of using the technique in real time.
DL	[96]	IEEE39	—	—	A detailed description of the mathematical apparatus of FC analysis is given.	An analysis of the operation of the proposed method on real data is not provided. There is no study of the possibility of using the technique in real time.
CNN-LSTM	[97]	1. IEEE39, 2. ACTIVSg500	1. DataSetSize = 1800 1. Delay = 0.024 s 2. DataSetSize = 2100 2. Delay = 0.0105 s	1. 99.57 2. 99.99	High accuracy, high performance, large number of test examples.	There is no study of the possibility of using the technique in real time.
FL	[98]	IEEE39	—	—	The original FC technique is presented.	A detailed study of the possibility of using the technique in real time is not provided.
CNN	[99]	IEEE39	DataSetSize = 1400	99.76	High efficiency.	There is no study of the possibility of using the technique in real time.
SVM	[100]	1. IEEE39, 2. NRPG 246	Delay = 0.02 s	1. 96.80, 2. 99.00	High efficiency.	There is no study of the possibility of using the technique in real time.
RL	[101]	Kundur’s Two-Area System	—	—	The influence of wind generation and electric vehicles on the process of frequency change in EPS is considered.	An analysis of the operation of the proposed method on real data is not provided.
PSO	[102]	—	—	—	The method does not involve specifying the parameters of EPS mathematical models.	An analysis of the operation of the proposed method on real data is not provided. A detailed description of the data sample is not provided.
BBBC	[103]	75-bus real power system	—	—	High performance, less overshoot, and fast damping of the regulation process.	A detailed description of the data sample is not provided.
GA	[104]	75-bus real power system	—	—	High performance.	No detailed description of the data sample is provided.
DQN	[105]	IEEE25	Delay = 0.02 s	—	High performance.	No detailed description of the data sample is provided.

shows the characteristics of each method proposed in the reviewed works. In Table 8, the following notations are used: DataSetSize—data sampling volume; Delay—algorithm delay for identifying one class.

9. Conclusions

A modern EPS has several features directly impacting the accuracy and adaptability of EC for preserving SSS, TS, and VS. These features are due to a decrease in total inertia, increased stochasticity in electricity generation and consumption processes, and a rise in the number of digital EPS monitoring and control devices.

To enhance the accuracy and adaptability of EC in EPS, ML algorithms are extensively used. These algorithms cater to a variety of tasks: identifying the type of disturbance, and selecting CAs for maintaining SSS, TS, VS, and FC. This paper presents an overview of ML algorithms for application in EPS EC.

The following conclusions are drawn from the above review:

• The most common class of algorithms used to select optimal CAs is DL. These algorithms make it possible to identify hidden correlations in the source data, while they do not assume a linear separation of classes in the data sample.
• To test the proposed methods for using ML algorithms for EPS EC, mathematical models are most often used. The most commonly used model is IEEE39.
• Most works provide a detailed description of the data sample but do not provide recommendations regarding the minimum number of scenarios that must be considered when training the model. In the reviewed works, the data sample size varies from 100 to 56,000.
• In the reviewed works, there is no description of the requirements for the homogeneity and representativeness of the source data.
• The authors of most of the reviewed studies determine the magnitude of the classification delay; however, there is virtually no information on real-time testing.
• In a minority of cases, the authors use real data for testing, which is characterized by the presence of noise, emissions, and omissions.
• When using data from PMU, methods combining ML and DSP algorithms are widely used.

In the practical implementation of ML algorithms for EPS management, the following difficulties can be identified: the problem of forming a representative training sample, the ambiguity of the procedure for determining an acceptable ML algorithm, and possible problems of ensuring cybersecurity and interpretability of CAs selection results. Therefore, the following directions for future research can be identified (they are listed in order of priority):

• Development of universal methods for analyzing SSS, TS, and VS.
• Testing algorithms on data obtained from real EPS. Development of algorithms for processing noise and outliers in the source data.
• Testing and adaptation of algorithms for real-time operation with determination of acceptable delays and sampling rates of source signals.
• Active use of data obtained from PMU.
• Development of infrastructure solutions for integrating ML algorithms into existing EPS EC business processes.
• Study of cybersecurity problems of using ML algorithms in the operational control loop of EPS. This task is the subject of extensive research and involves the analysis of possible risk factors for measuring systems, EPS management, data storage, and processing.
• Definition of requirements for the composition of data sampling features for EPS EC, minimum volume, completeness, and representativeness of the dataset.
• Analysis of microgrid EC methods with a significant proportion of RES, electricity consumers, distributed generation, and consumers with price-dependent consumption [112,113,114,115].

From the above analysis, it can be concluded that the most flexible ML algorithms are algorithms that simulate the functioning of biological neural networks (ANN, CANN, ChNN, CNN) providing the required level of flexibility and accuracy in the conditions of the dynamic nature of EPS. These algorithms identify hidden correlations in the data, ensure reliable operation in conditions of noise and emissions in the data, and also have high performance, which is essential for EPS management. These algorithms can effectively solve all the problems considered in this paper.

Based on the conducted meta-analysis, there were no shortcomings in the considered ML algorithms. The main difficulty in developing the methods of EPS EC based on ML algorithms relates to the formation of a data sample, as well as to the testing of trained algorithms in real time and on real data.

The integration of ML algorithms into existing EPS management business processes is not difficult, because existing solutions can be used in terms of changes and data collection of established and transient processes, and existing communication channels can be used to transmit signals to the implementation of CAs. New tasks in the practical implementation of ML algorithms in the operational EPS control loop are the organization of ML algorithms training, testing, administration, cybersecurity, and reliability.