<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1491846</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2024.1491846</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Two-stage transient stability assessment using ensemble learning and cost sensitivity</article-title>
<alt-title alt-title-type="left-running-head">Lei et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2024.1491846">10.3389/fenrg.2024.1491846</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Lei</surname>
<given-names>Aoyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mei</surname>
<given-names>Yong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ma</surname>
<given-names>Dexin</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2833815/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Zining</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tao</surname>
<given-names>Wenwei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Huang</surname>
<given-names>Fangneng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Power Dispatch and Control Center of China Southern Power Grid</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>China Southern Power Grid Electric Power Research Institute</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1464268/overview">Cong Zhang</ext-link>, Hunan University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1518642/overview">Weike Mo</ext-link>, Jinan University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2634622/overview">Xiaodong Zheng</ext-link>, South China University of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Fangneng Huang, <email>leiay@csg.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>17</day>
<month>10</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1491846</elocation-id>
<history>
<date date-type="received">
<day>05</day>
<month>09</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>25</day>
<month>09</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Lei, Mei, Ma, Liu, Tao and Huang.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Lei, Mei, Ma, Liu, Tao and Huang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>This paper proposes a novel two-stage transient stability assessment (TSA) model that integrates ensemble learning with cost sensitivity to address the challenges posed by the integration of renewable energy and load fluctuations. The model employs CNNs as positive and negative classifiers to initially evaluate samples, with consistent results output directly. In cases of inconsistency, the sample is evaluated by a fair classifier, specifically an ELM, trained on critical samples. This approach significantly enhances the classification performance and credibility of the fair classifier, especially under imbalanced conditions, thereby improving the overall efficiency and accuracy of TSA. The proposed model demonstrates superior performance compared to single-stage models and other two-stage models, achieving high accuracy and robustness in transient stability assessment, particularly for critical samples.</p>
</abstract>
<kwd-group>
<kwd>transient stability</kwd>
<kwd>cost sensitivity</kwd>
<kwd>critical sample</kwd>
<kwd>ensemble learning</kwd>
<kwd>sample imbalance</kwd>
</kwd-group>
<contract-num rid="cn001">ZDKJXM20210063</contract-num>
<contract-sponsor id="cn001">China Southern Power Grid<named-content content-type="fundref-id">10.13039/501100005311</named-content>
</contract-sponsor>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Smart Grids</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The development of modern power systems presents new challenges to the safety of these systems. In recent years, frequent blackouts worldwide have had significant economic impacts (<xref ref-type="bibr" rid="B23">Wei Zhang and Zhang, 202</xref>). Transient stability is crucial to the safe operation of power systems. Once a disturbance occurs in the power system, early assessment of transient stability and timely implementation of emergency control measures can protect the system&#x2019;s stability (<xref ref-type="bibr" rid="B7">Guo et al., 2023</xref>; <xref ref-type="bibr" rid="B16">Meridji et al., 2023</xref>; <xref ref-type="bibr" rid="B19">Singh and Chauhan, 2023</xref>; <xref ref-type="bibr" rid="B22">Wang Y. et al., 2023</xref>; <xref ref-type="bibr" rid="B27">Zhu et al., 2023</xref>). Therefore, it is vital and necessary to design a TSA model that can quickly and accurately determine whether a power system will become unstable and infer the type of events that may cause instability.</p>
<p>Classical transient stability methods include time-domain simulation (TDS) and direct methods. TDS analyzes system stability through detailed modeling. As long as the component models and network structures are sufficiently accurate, the results are reliable. However, with the increasing complexity of power systems and rapid changes in electronic hardware, TDS-based models have become more complex, and their computational demands have increased, making TDS methods unsuitable for real-time transient stability assessment. The direct method, on the other hand, does not require numerical integration of dynamic processes and can quantitatively analyze system transient stability. However, direct method models are simple and may face applicability issues in complex power systems. Additionally, the evaluation results tend to be conservative. As a result, these methods currently cannot fully meet the practical needs of online transient stability assessment. To date, real-world power grids still require more advanced TSA solutions that can achieve high reliability and efficiency.</p>
<p>The rapid development of big data offers new ideas for transient stability assessment. TSA rules can be mined from large volumes of data using machine learning (ML) methods (<xref ref-type="bibr" rid="B15">Liu et al., 2023</xref>; <xref ref-type="bibr" rid="B13">Li et al., 2024</xref>; <xref ref-type="bibr" rid="B17">Shao et al., 2024</xref>). Once dynamic information is received, stability results can be output within milliseconds. After discovering underlying relationships offline through ML, they can be easily applied to online TSA. This data-driven approach almost eliminates online computation time and can effectively predict system stability.</p>
<p>Modern power systems usually remain stable after most disturbances due to their robustness. In practical applications, events causing power system instability are rare. These rare events lead to significant imbalances between stable (majority class) and unstable (minority class) transient data distributions. The imbalance in training data leads to high misclassification rates in TSA models. Moreover, misjudging the transient stability of power systems prevents timely responses and corresponding measures, leading to safety accidents and significant economic losses. Therefore, addressing the imbalance in training data is not only crucial for improving TSA accuracy but also has practical significance.</p>
<p>Existing methods for addressing class imbalance can be broadly categorized into data-driven and algorithm-driven approaches. Additionally, traditional data-driven methods attempt to oversample or undersample the data so that different classes appear in equal proportions in the training data. The main drawback of sampling-based methods is that they may lose important information in undersampling or overfit the training data in oversampling. Algorithm-driven methods primarily include cost-sensitive methods (<xref ref-type="bibr" rid="B20">Wang H. et al., 2023</xref>; <xref ref-type="bibr" rid="B3">Chen et al., 2022</xref>; <xref ref-type="bibr" rid="B14">Lin et al., 2022</xref>; <xref ref-type="bibr" rid="B21">Wang et al., 2020</xref>) and ensemble learning methods (<xref ref-type="bibr" rid="B18">Shen, 2023</xref>; <xref ref-type="bibr" rid="B2">Chen and Wang, 2021</xref>; <xref ref-type="bibr" rid="B1">Chen et al., 2021</xref>; <xref ref-type="bibr" rid="B26">Zhao et al., 2022</xref>; <xref ref-type="bibr" rid="B24">Wu et al., 2020</xref>). References (<xref ref-type="bibr" rid="B20">Wang H. et al., 2023</xref>) consider the imbalance in sample distribution in space by dividing samples into different levels of training sets and training the model using cost-sensitive approaches. Reference (<xref ref-type="bibr" rid="B18">Shen, 2023</xref>) proposes a comprehensive transient stability state assessment method for power systems based on machine learning, using multiple classifiers to alleviate the sample imbalance problem in transient stability samples. However, its drawbacks are also evident. The ensemble learning classifiers used are often weak classifiers with low resistance to interference and increased training time. Among these methods, the application of cost sensitivity in transient stability assessment is relatively limited. Although it helps mitigate sample imbalance, its blind use may lead to overfitting and increased misclassification rates. Applying ensemble learning alone to transient stability assessment can also mitigate sample imbalance, but due to the large number of ensemble classifiers, it may result in longer training times and weaker resistance to interference.</p>
<p>Recently, to mitigate the imbalance in transient stability samples, references (<xref ref-type="bibr" rid="B2">Chen and Wang, 2021</xref>; <xref ref-type="bibr" rid="B21">Wang et al., 2020</xref>) adopted a two-stage assessment approach combining cost sensitivity and ensemble learning. Specifically, three models are used, with two classifiers in the first stage, each biased towards one class of samples using cost-sensitive methods. A non-cost-sensitive classifier is then trained using the same training set as a fair classifier for second-stage evaluation. Although this approach combines the advantages of cost-sensitive learning and ensemble learning, it also overcomes some disadvantages. However, this method lacks credibility, mainly because the credibility of the fair classifier is low. This is because hard samples are evaluated by a non-cost-sensitive classifier, and the evaluation results are directly output. However, since the fair classifier is not selected based on the characteristics of the input samples, highly accurate and convincing evaluation results cannot be generated.</p>
<p>In transient stability assessment, what truly affects misclassification and missed judgments are the hard samples. In this study, we propose a two-stage transient stability assessment model based on the integration of ensemble learning and cost sensitivity. First, the initial sample set is input into positive and negative classifiers again, and samples with inconsistent evaluation results are collected to form a critical sample set. This critical sample set is used to train the ELM (fair classifier). The model combines the strengths of cost sensitivity and ensemble learning, greatly enhancing the credibility of the fair classifier while improving the overall reliability of the model.</p>
<p>The rest of this paper is structured as follows: <xref ref-type="sec" rid="s2">Section 2</xref> offers an overview of the transient stability assessment (TSA) problem and introduces the proposed two-stage TSA model. <xref ref-type="sec" rid="s3">Section 3</xref> details the implementation process of the model. <xref ref-type="sec" rid="s4">Section 4</xref> presents case studies and evaluates the model&#x2019;s performance. Finally, <xref ref-type="sec" rid="s5">Section 5</xref> summarizes the conclusions of the paper.</p>
</sec>
<sec id="s2">
<title>2 Integrated TSA model</title>
<sec id="s2-1">
<title>2.1 Transient stability assessment problem</title>
<p>Transient stability assessment (TSA) has become increasingly important in power systems, with the main challenge being how to overcome the inherent sample imbalance problem, which is determined by the robustness of power systems. To address the imbalance in transient stability samples, a two-stage assessment method combining cost sensitivity and ensemble learning has been widely adopted. However, due to the lack of optimization in the selection of the fair classifier, it is difficult to produce highly accurate and convincing evaluation results. The model proposed in this paper addresses the problem of transient stability assessment under imbalanced samples by considering the difficulty of classifying critical samples, combining cost sensitivity and ensemble learning, which greatly enhances the reliability of the fair classifier and the overall model reliability. Currently, the issue of sample imbalance in transient stability assessment cannot be ignored, as stable samples are the majority, unstable samples are fewer, and critical samples are even rarer. The training model must fully understand the characteristics of stable, unstable, and critical samples to classify them correctly. Therefore, solving the sample imbalance problem is crucial for transient stability assessment. The two-stage TSA model proposed in this paper aims to achieve accurate and reliable classification results even under imbalanced conditions.</p>
</sec>
<sec id="s2-2">
<title>2.2 Selection of positive and negative models and application of cost sensitivity</title>
<sec id="s2-2-1">
<title>2.2.1 Selection of positive and negative models</title>
<p>Unlike most previous evaluation model designs, this paper designs a two-stage transient stability assessment model based on the integration of ensemble learning and cost sensitivity. The advantage of this approach is that it can combine the strengths of various models and effectively overcome the drawbacks of low accuracy, low effectiveness, and high cost caused by sample imbalance.</p>
<p>The positive and negative models are the foundational models in the entire two-stage framework; their results are only valid when they are consistent. Otherwise, if the results are inconsistent, the assessment will enter the next stage without directly outputting the results. Therefore, to maximize the role of the positive and negative models, cost sensitivity is applied to these models. The positive model is inclined to evaluate samples as stable, while the negative model is inclined to evaluate samples as unstable. Under this design, when the evaluation results of the positive and negative models are consistent, the input samples must be far from the boundary, indicating that the samples are either extremely stable or extremely unstable. In such cases, the positive and negative models require strong feature extraction capabilities to fully recognize samples beyond the boundary.</p>
<p>Convolutional Neural Networks (CNNs) possess strong nonlinear expression and pattern recognition capabilities, making them suitable for handling complex dynamic processes and nonlinear problems in power systems (<xref ref-type="bibr" rid="B11">Lee et al., 2023</xref>). This allows CNN-based TSA models to more accurately reflect the transient stability status of the power system based on input operational parameters. Compared to traditional time-domain and direct methods, CNN-based TSA methods offer higher assessment accuracy and faster computation speed. They can meet the requirements of online real-time prediction, providing timely and accurate power system stability analysis results for dispatchers (<xref ref-type="bibr" rid="B6">Gu et al., 2024</xref>). Additionally, CNN models can adapt to different power system scenarios and demands by adjusting the network structure and parameters, offering higher flexibility and scalability. In summary, applying CNNs to power system transient stability assessment has advantages such as automatic feature extraction, strong nonlinear processing capabilities, high accuracy, real-time performance, and applicability to complex systems, making it a vital tool in power system operation and control. Therefore, we select CNN as the positive and negative models. With CNN&#x2019;s strong feature extraction capabilities, we can effectively handle the characteristics of stable and unstable sample regions, achieving initial simple classification. Furthermore, CNN-based methods can directly use low-level measurement data as input features, extracting multi-granularity information from measurement data through multi-size convolution kernels, thereby improving the accuracy of TSA in power systems (<xref ref-type="bibr" rid="B9">Jin et al., 2023</xref>).</p>
<p>For transient stability assessment, reasonable feature selection is needed to construct samples as input. This paper focuses on selecting steady-state quantities and fault characteristics for the forward-looking prediction of transient stability (<xref ref-type="bibr" rid="B8">Ji et al., 2022</xref>). Therefore, bus voltage, generator active and reactive power, and load active and reactive power are chosen as input features. Relevant data are obtained from measurements or time-domain simulations and input into the CNN model. The network structure of CNN is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>CNN network structure.</p>
</caption>
<graphic xlink:href="fenrg-12-1491846-g001.tif"/>
</fig>
<p>Based on the input electrical quantity matrix, the convolutional layer performs local learning through convolution kernels in the CNN model. The convolution operation for each convolution kernel is as <xref ref-type="disp-formula" rid="e1">Equation 1</xref>.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>&#x2a;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the output of the <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th convolutional surface, <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the weight matrix corresponding to the <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th convolution kernel, <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the input matrix, <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the bias term, and <inline-formula id="inf7">
<mml:math id="m8">
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the activation function (the ReLU function is selected here).</p>
<p>For the pooling layer, this paper adopts max pooling, as shown in <xref ref-type="disp-formula" rid="e2">Equation 2</xref>.<disp-formula id="e2">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>max</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>,</mml:mo>
<mml:mspace width="1em"/>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mspace width="1em"/>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the <inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th pooling surface, <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is a sub-block of the output matrix from the previous convolutional layer, and <inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the dimension of the sub-block of the output matrix from the previous convolutional layer.</p>
<p>The calculation formula for the fully connected layer is as <xref ref-type="disp-formula" rid="e3">Equation 3</xref>.<disp-formula id="e3">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>In the fully connected layer, <inline-formula id="inf12">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the output, <inline-formula id="inf13">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the input, <inline-formula id="inf14">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the weight matrix, and <inline-formula id="inf15">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the bias term. The ReLU function is selected for the activation function.</p>
<p>The fully connected layer passes the output values to the output layer, and the expression of the sigmoid function in the output layer is as <xref ref-type="disp-formula" rid="e4">Equation 4</xref>.<disp-formula id="e4">
<mml:math id="m19">
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <inline-formula id="inf16">
<mml:math id="m20">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the input value and <inline-formula id="inf17">
<mml:math id="m21">
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the function output.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Application of cost sensitivity</title>
<p>Typically, the approach to solving quantity imbalance is to assign higher weights to unstable samples. However, the misclassification of critical samples still occurs. Another approach is to assign higher weights to both imbalances simultaneously: first, higher weights are given to unstable samples, and then higher weights are assigned to misclassified samples. This approach can lead to overfitting (<xref ref-type="bibr" rid="B10">Kesici et al., 2023</xref>).</p>
<p>In this study, the model is trained by optimizing parameters to improve the fit to the training samples. The loss function is used to measure the difference between the predicted and actual values, and the impact of training samples on evaluation rules can be assessed through the loss function. This cost-sensitive approach makes the positive model inclined to evaluate samples as stable and the negative model inclined to evaluate samples as unstable.</p>
<p>For the transient stability assessment problem, the weight correction coefficient for model parameters can be represented by <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Weight coefficient matrix for model parameter modification.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Real label</th>
<th colspan="2" align="center">Predicted label</th>
</tr>
<tr>
<th align="center">Stable</th>
<th align="center">Unstable</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Stable</td>
<td align="center">
<inline-formula id="inf18">
<mml:math id="m22">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1,1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf19">
<mml:math id="m23">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">Unstable</td>
<td align="center">
<inline-formula id="inf20">
<mml:math id="m24">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf21">
<mml:math id="m25">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0,0</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Where <inline-formula id="inf22">
<mml:math id="m26">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1,1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0,0</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. <inline-formula id="inf23">
<mml:math id="m27">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf24">
<mml:math id="m28">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> represent the weight correction coefficients for model parameters when predicting stable conditions as unstable and unstable conditions as stable, respectively. Typically, <inline-formula id="inf25">
<mml:math id="m29">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>In different scenarios, the proportion of unstable samples to stable samples in the training set often differs, which causes the model trained under imbalanced samples to have a certain bias in its discriminative results. Therefore, this paper introduces a correction factor <inline-formula id="inf26">
<mml:math id="m30">
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> to adjust for the imbalance in sample numbers in different scenarios, as shown in <xref ref-type="disp-formula" rid="e5">Equations 5</xref>, <xref ref-type="disp-formula" rid="e6">6</xref>.<disp-formula id="e5">
<mml:math id="m31">
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m32">
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">uns</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where <inline-formula id="inf27">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the number of stable samples in the training set; <inline-formula id="inf28">
<mml:math id="m34">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">uns</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the number of unstable samples in the training set.</p>
<p>By adjusting the weight correction coefficients during the evaluation model parameter correction process, the model&#x2019;s loss function is as follows:<disp-formula id="e7">
<mml:math id="m35">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2032;</mml:mo>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:mi>&#x3b3;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x2212;</mml:mo>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
<mml:mi>ln</mml:mi>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
<mml:mi>&#x3b3;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>The model using <xref ref-type="disp-formula" rid="e7">Equation 7</xref> is the positive model. After applying cost sensitivity, this model tends to evaluate input samples as stable. Conversely, if the model needs to incline toward evaluating input samples as unstable, an additional correction factor <inline-formula id="inf29">
<mml:math id="m36">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is required. Since the original model tends to evaluate samples as stable due to sample imbalance, adding the coefficient <inline-formula id="inf30">
<mml:math id="m37">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> ensures that the model tends to evaluate input samples as unstable. Similarly, it is also necessary to adjust the imbalance of the number of samples in different scenarios here, as shown in <xref ref-type="disp-formula" rid="e8">Equations 8</xref>, <xref ref-type="disp-formula" rid="e9">9</xref>.<disp-formula id="e8">
<mml:math id="m38">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">uns</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
<disp-formula id="e9">
<mml:math id="m39">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>When <inline-formula id="inf31">
<mml:math id="m40">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, the evaluation model merely corrects for the bias introduced by the imbalanced training samples; when <inline-formula id="inf32">
<mml:math id="m41">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, the model increases its fit to unstable samples, reducing the misclassification probability for unstable samples.</p>
<p>The reorganized model&#x2019;s loss function is shown in <xref ref-type="disp-formula" rid="e10">Equation 10</xref>.<disp-formula id="e10">
<mml:math id="m42">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2032;</mml:mo>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1,0</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
<mml:mi>ln</mml:mi>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
<mml:mi>&#x3b8;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>0,1</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>The negative model applies the loss function from the above equation to tend toward evaluating input samples as unstable.</p>
</sec>
</sec>
<sec id="s2-3">
<title>2.3 Fair classifier</title>
<p>The fair classifier is used to classify hard samples, which are samples that the initial simple classification could not correctly categorize. Therefore, the sample set used to train the fair classifier is the critical sample set, and the model selection must consider this.</p>
<p>The fair classifier operates in the second stage of the assessment, requiring high accuracy and interpretability in its results. The sample set used to train the fair classifier is the critical sample set, which contains rich data information and high information entropy, requiring a suitable model to extract these features. For ELM training, the ideal sample should be data close to the stability boundary to obtain refined stability rules, allowing for effective learning of the critical region and a strong ability to correctly classify hard samples.</p>
<p>Given the above considerations, ELM is chosen as the fair classifier. Due to the distribution characteristics of the critical sample set, it matches the training requirements of the ELM model. If the training samples of ELM are from the critical region, the trained ELM will have a strong ability to correctly classify the critical sample set. ELM is used for training, but not with the original sample set; instead, it is trained using samples with conflicting evaluation results from two classifiers with opposing biases.</p>
<p>ELM can randomly select input layer parameters, and then use the Moore-Penrose generalized inverse to obtain the output layer weights with the smallest 2-norm. In ELM, only the number of hidden layer neuron nodes needs to be learned and adjusted, and the entire process does not require iteration. <xref ref-type="fig" rid="F2">Figure 2</xref> shows the ELM network structure.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Elm network structure.</p>
</caption>
<graphic xlink:href="fenrg-12-1491846-g002.tif"/>
</fig>
<p>Specifically, for a dataset <inline-formula id="inf33">
<mml:math id="m43">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> with <inline-formula id="inf34">
<mml:math id="m44">
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> samples, where <inline-formula id="inf35">
<mml:math id="m45">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="double-struck">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf36">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="double-struck">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf37">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the class label of the <inline-formula id="inf38">
<mml:math id="m48">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th sample, and the ELM output with <inline-formula id="inf39">
<mml:math id="m49">
<mml:mrow>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> hidden layer units is shown in <xref ref-type="disp-formula" rid="e11">Equation 11</xref>.<disp-formula id="e11">
<mml:math id="m50">
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>h</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where <inline-formula id="inf40">
<mml:math id="m51">
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf41">
<mml:math id="m52">
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the activation function; <inline-formula id="inf42">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="double-struck">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the weight vector of the <inline-formula id="inf43">
<mml:math id="m54">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th unit in the hidden layer; <inline-formula id="inf44">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="double-struck">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the multiplier of the <inline-formula id="inf45">
<mml:math id="m56">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th unit in the hidden layer; and <inline-formula id="inf46">
<mml:math id="m57">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="double-struck">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the bias of the <inline-formula id="inf47">
<mml:math id="m58">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-th unit in the hidden layer.</p>
<p>In ELM, the weight vector <inline-formula id="inf48">
<mml:math id="m59">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and the bias <inline-formula id="inf49">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are randomly assigned, and the multiplier <inline-formula id="inf50">
<mml:math id="m61">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> can be obtained based on all the data through the Moore-Penrose generalized inverse. Therefore, there is no need for the backpropagation (BP) neural network&#x2019;s necessary back-adjustment process (<xref ref-type="bibr" rid="B5">Fuqiang et al., 2023</xref>), making ELM&#x2019;s calculation speed significantly faster than classification algorithms based on optimization. This reduction in calculation time also reduces the online application time of the two-stage evaluation model (<xref ref-type="bibr" rid="B25">Zhang et al., 2017</xref>).</p>
</sec>
<sec id="s2-4">
<title>2.4 Overall evaluation process</title>
<p>After completing all the above designs, the framework for the proposed two-stage transient stability assessment model is established. Next, the overall evaluation process is designed.</p>
<p>The model proposed in this study is a two-stage evaluation model, where the second stage&#x2019;s evaluation somewhat depends on the results of the first stage. Therefore, it is necessary to design scoring for the first-stage evaluation results to complete the evaluation logic of the first and second stages.</p>
<p>Define <inline-formula id="inf51">
<mml:math id="m62">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf52">
<mml:math id="m63">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf53">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula id="inf54">
<mml:math id="m65">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the total score of the first stage, and <inline-formula id="inf55">
<mml:math id="m66">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf56">
<mml:math id="m67">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the scores of the positive and negative classifiers, respectively.</p>
<p>To improve efficiency, <inline-formula id="inf57">
<mml:math id="m68">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf58">
<mml:math id="m69">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are linked to the output results of the positive and negative classifiers. If the positive classifier output is stable (label 1), then <inline-formula id="inf59">
<mml:math id="m70">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is 1; otherwise, it is 0.</p>
<p>Define <inline-formula id="inf60">
<mml:math id="m71">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, the scoring logic is summarized in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Scoring table.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<inline-formula id="inf61">
<mml:math id="m72">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">
<inline-formula id="inf62">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">
<inline-formula id="inf63">
<mml:math id="m74">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<inline-formula id="inf64">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">0</td>
<td align="center">1</td>
</tr>
<tr>
<td align="center">
<inline-formula id="inf65">
<mml:math id="m76">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">1</td>
<td align="center">2</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to the model logic proposed in this paper, the first-stage evaluation results are credible, and the evaluation results are directly output only when the positive and negative classifier outputs are consistent. This occurs when <inline-formula id="inf66">
<mml:math id="m77">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is 0 or 2, meaning the first-stage evaluation results are reliable, and there is no need for the second-stage evaluation. If the positive and negative classifier outputs are inconsistent, the first-stage evaluation results are unreliable, and <inline-formula id="inf67">
<mml:math id="m78">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> equals 1, then the process proceeds to the second-stage evaluation.</p>
<p>The overall evaluation process can be simplified as shown in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Evaluation process.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<inline-formula id="inf68">
<mml:math id="m79">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="center">First stage</th>
<th align="center">Second stage</th>
<th align="center">Total score</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">0</td>
<td align="center">1</td>
<td align="center">0</td>
<td align="center">1</td>
</tr>
<tr>
<td align="center">1</td>
<td align="center">1</td>
<td align="center">1</td>
<td align="center">2</td>
</tr>
<tr>
<td align="center">2</td>
<td align="center">1</td>
<td align="center">0</td>
<td align="center">1</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>As shown in <xref ref-type="table" rid="T3">Table 3</xref>, when the total score is 1, the second-stage evaluation result is output as the model&#x2019;s evaluation result. When the total score is 2, that is, when <inline-formula id="inf73">
<mml:math id="m84">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> equals 1, the first-stage evaluation result is credible, and the first-stage evaluation result is output as the model&#x2019;s evaluation result.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Implementation process of the integrated TSA model</title>
<p>After selecting a specific example system, transient stability simulations are performed to collect steady-state data and create the initial sample set.</p>
<p>Based on the obtained initial dataset, the positive and negative classifiers are trained separately. After training, the initial dataset is re-input into the trained positive and negative classifiers for evaluation, and the critical dataset is organized based on the differing evaluation results from the positive and negative classifiers.</p>
<p>Since the fair classifier needs to correctly classify critical samples, ELM is selected as the fair classifier. By training ELM using the critical sample set, the ability of ELM to correctly classify critical samples is greatly enhanced, meeting the requirements of the fair classifier. The construction of the integrated TSA model is then completed.</p>
<p>When evaluating real-time data obtained from the power system, the measured data is input into the positive and negative classifiers for evaluation. If the evaluation results are consistent, the results are output, completing the assessment; if the evaluation results are inconsistent, the measured data is input into the fair classifier, and the evaluation results from the fair classifier are output, completing the assessment. The specific implementation process is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.<list list-type="simple">
<list-item>
<p>(1): Offline Training</p>
</list-item>
<list-item>
<p>(2): Online Application</p>
</list-item>
</list>
</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>(1) Offline training; (2) Online application.</p>
</caption>
<graphic xlink:href="fenrg-12-1491846-g003.tif"/>
</fig>
</sec>
<sec id="s4">
<title>4 Case analysis</title>
<p>Taking the New England 10-generator, 39-bus system as a case study, as shown in the <xref ref-type="fig" rid="F4">Figure 4</xref>, the entire bus system is divided into four regions according to adjacent buses. Four different load levels, 80%, 90%, 110%, and 120%, are considered, and each region is assigned a different load level, resulting in 256 different power flows,256 different operating modes. The generator output is adjusted accordingly to ensure the convergence of power flow calculations. The fault is set as a three-phase short-circuit fault, and an N-1 scan is performed on the system by sequentially setting faults on different lines. The fault duration is set to 0.2 s on lines without transformers, with the fault point located at 50% of the line&#x2019;s length. On lines with transformers, the fault point is set at the head section of the line.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Line diagram of IEEE 39-bus power system.</p>
</caption>
<graphic xlink:href="fenrg-12-1491846-g004.tif"/>
</fig>
<p>The selection of input features including bus voltage, generator active power and reactive power, as well as load active power and reactive power, is indeed widely used in data-driven transient stability assessment (<xref ref-type="bibr" rid="B12">Li et al., 2021</xref>; <xref ref-type="bibr" rid="B4">Du et al., 2021</xref>). Bus voltage reflects the system&#x2019;s voltage levels and can serve as an early warning for instability. Generator active power and reactive power respectively influence frequency stability and voltage regulation. Load active power affects the system&#x2019;s energy balance, while load reactive power is crucial for voltage stability. These features provide a solid foundation for predictive transient stability assessment.</p>
<p>A total of 11,776 samples are generated. By sampling from these 11,776 simulated samples, the training set, testing set, and validation set are obtained in a ratio of 8:1:1.</p>
<p>Transient stability assessment in power systems includes two categories: stable and unstable. The criterion for stability is the transient stability index (TSI). The details are given in <xref ref-type="disp-formula" rid="e12">Equation 12</xref>.<disp-formula id="e12">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">TSI</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>360</mml:mn>
<mml:mo>&#xb0;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo stretchy="false">&#x7c;</mml:mo>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">&#x7c;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mn>360</mml:mn>
<mml:mo>&#xb0;</mml:mo>
<mml:mo>&#x2b;</mml:mo>
<mml:mo stretchy="false">&#x7c;</mml:mo>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">&#x7c;</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>When the maximum rotor angle difference among generators exceeds <inline-formula id="inf74">
<mml:math id="m86">
<mml:mrow>
<mml:mn>360</mml:mn>
<mml:mo>&#xb0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>, TSI is less than 0, and the system is considered unstable; when the rotor angle difference is less than <inline-formula id="inf75">
<mml:math id="m87">
<mml:mrow>
<mml:mn>360</mml:mn>
<mml:mo>&#xb0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>, TSI is greater than 0, and the system is considered stable.</p>
<p>The following indicators are set as in <xref ref-type="disp-formula" rid="e13">Equations 13</xref>&#x2013;<xref ref-type="disp-formula" rid="e15">15</xref> to evaluate the performance of the model:<disp-formula id="e13">
<mml:math id="m88">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">ACC</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x22c5;</mml:mo>
<mml:mn>100</mml:mn>
<mml:mi>%</mml:mi>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m89">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x22c5;</mml:mo>
<mml:mn>100</mml:mn>
<mml:mi>%</mml:mi>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="e15">
<mml:math id="m90">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mi>A</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x22c5;</mml:mo>
<mml:mn>100</mml:mn>
<mml:mi>%</mml:mi>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
</p>
<p>After the training, the model is tested from the test set, and the performance test results of the model are obtained. At the same time, existing Two-Stage models (such as two-stage-LSTM, two-stage-GRU) are tested, which all use the same neural network (such as LSTM or GRU) for positive and negative sample classification and fair classification tasks. The weight coefficient matrix used for model parameter correction is shown in <xref ref-type="table" rid="T4">Table 4</xref>. Firstly, transient stability evaluation experiments were conducted on the single-stage model and the two-stage model, and the results were recorded and summarized as shown in <xref ref-type="table" rid="T5">Table 5</xref>.Vertical comparison shows that the two-stage method proposed in this paper outperforms the single-stage method in all indicators, with a significant reduction in misclassification rates and a marked improvement in classification accuracy. This validates the effectiveness of the two-stage model in complex tasks compared to a single model. <xref ref-type="fig" rid="F5">Figure 5</xref> shows the bar chart comparison of the performance between single-stage and two-stage models.Based on these findings, we further explored the effects of different neural network combinations within the two-stage framework. As such, this paper proposes the TCS-ELM model (Two-Stage Cost-Sensitive Ensemble Learning Model), using CNN for positive and negative sample classification in the first stage and ELM for fair classification in the second stage, fully leveraging CNN&#x2019;s feature extraction strengths and ELM&#x2019;s efficient classification capabilities. <xref ref-type="table" rid="T6">Table 6</xref> shows the performance testing results of various two-stage models and the proposed TCS-ELM model.Results indicate that the proposed TCS-ELM model, by using CNN in the first stage for stronger feature extraction capabilities and ELM in the second stage for improved classification accuracy, outperforms other two-stage models (such as Two-Stage-LSTM and Two-Stage-GRU) in all indicators, particularly exhibiting higher accuracy and robustness when handling critical samples. This demonstrates that the TCS-ELM model combines the advantages of CNN and ELM in their respective stages, significantly enhancing overall performance. <xref ref-type="fig" rid="F6">Figure 6</xref> presents the bar chart comparison of the performance between TCS-ELM and other two-stage models.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Weight coefficient matrix for model parameter modification.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Real label</th>
<th colspan="2" align="center">Predicted label</th>
</tr>
<tr>
<th align="center">Stable</th>
<th align="center">Unstable</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Stable</td>
<td align="center">
<inline-formula id="inf69">
<mml:math id="m80">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf70">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td align="center">Unstable</td>
<td align="center">
<inline-formula id="inf71">
<mml:math id="m82">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="center">
<inline-formula id="inf72">
<mml:math id="m83">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Comparison of single-stage and two-stage models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Model</th>
<th align="center">PACC</th>
<th align="center">PFD</th>
<th align="center">PFA</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">LSTM</td>
<td align="center">0.9202</td>
<td align="center">0.6535</td>
<td align="center">0.4444</td>
</tr>
<tr>
<td align="center">Two-Stage-LSTM</td>
<td align="center">0.9363</td>
<td align="center">0.5644</td>
<td align="center">0.2903</td>
</tr>
<tr>
<td align="center">GRU</td>
<td align="center">0.9253</td>
<td align="center">0.6238</td>
<td align="center">0.3968</td>
</tr>
<tr>
<td align="center">Two-Stage-GRU</td>
<td align="center">0.9465</td>
<td align="center">0.4257</td>
<td align="center">0.2564</td>
</tr>
<tr>
<td align="center">CNN</td>
<td align="center">0.9482</td>
<td align="center">0.3663</td>
<td align="center">0.2727</td>
</tr>
<tr>
<td align="center">Two-Stage-CNN</td>
<td align="center">0.9533</td>
<td align="center">0.2376</td>
<td align="center">0.2870</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Comparison of performance between single-stage and two-stage models. Model performance improvement: two-stage vs. baseline (LSTM, GRU, CNN).</p>
</caption>
<graphic xlink:href="fenrg-12-1491846-g005.tif"/>
</fig>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Comparison of TCS-ELM with other two-stage models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Model</th>
<th align="center">PACC</th>
<th align="center">PFD</th>
<th align="center">PFA</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">TCS-ELM</td>
<td align="center">0.9779</td>
<td align="center">0.1584</td>
<td align="center">0.1053</td>
</tr>
<tr>
<td align="center">Two-Stage-CNN</td>
<td align="center">0.9533</td>
<td align="center">0.2376</td>
<td align="center">0.2870</td>
</tr>
<tr>
<td align="center">Two-Stage-GRU</td>
<td align="center">0.9465</td>
<td align="center">0.4257</td>
<td align="center">0.2564</td>
</tr>
<tr>
<td align="center">Two-Stage-LSTM</td>
<td align="center">0.9363</td>
<td align="center">0.5644</td>
<td align="center">0.2903</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Comparison of performance between single-stage and two-stage models. Performance comparison of 4 deep learning models.</p>
</caption>
<graphic xlink:href="fenrg-12-1491846-g006.tif"/>
</fig>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>In conclusion, the two-stage transient stability assessment model proposed in this paper, which leverages a combination of cost sensitivity and ensemble learning, addresses the key challenge of sample imbalance while exhibiting robust classification capabilities, particularly for hard-to-classify samples. By integrating these advanced techniques, the model demonstrates exceptional performance, achieving high levels of accuracy without compromising the reliability and credibility of the evaluation results.</p>
<p>In contrast to widely used transient stability assessment methods, this model successfully overcomes several of their inherent limitations. Traditional models often struggle to balance speed, accuracy, and reliability, particularly in cases where the data distribution is uneven or where outliers skew the results. The proposed model not only mitigates these issues but also delivers rapid assessments, making it highly efficient for real-time or near-real-time applications. Its ability to provide convincing and trustworthy evaluations further strengthens its utility in practical scenarios where precise and timely decision-making is crucial, such as in power grid operations and other dynamic systems.</p>
<p>By offering a superior balance of speed, accuracy, and evaluation credibility, this model stands out as a significant advancement in the field of transient stability assessment, offering a practical solution that is both scalable and adaptable to diverse operating conditions.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>AL: Writing&#x2013;original draft. YM: Writing&#x2013;review and editing. DM: Writing&#x2013;review and editing. ZL: Writing&#x2013;review and editing. WT: Writing&#x2013;review and editing. FH: Writing&#x2013;original draft.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study received funding from the Science and Technology Project of China Southern Power Grid Co. Ltd. (ZDKJXM20210063).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
<p>The authors declare that this study received funding from the Science and Technology Project of China Southern Power Grid Co. Ltd. (ZDKJXM20210063). The funder had the following involvement with the study: the study design, the writing of this article.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Transient stability assessment model with parallel structure and data augmentation</article-title>. <source>Int. Trans. Electr. Energy Syst.</source> <volume>31</volume>, <fpage>12872</fpage>. <pub-id pub-id-type="doi">10.1002/2050-7038.12872</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Time-adaptive transient stability assessment based on gated recurrent unit</article-title>. <source>Int. J. Electr. Power and Energy Syst.</source> <volume>133</volume>, <fpage>107156</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2021.107156</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Imbalance correction method based on ratio of loss function values for transient stability assessment</article-title>. <source>CSEE J. Power Energy Syst</source>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Du</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A hierarchical power system transient stability assessment method considering sample imbalance</article-title>. <source>Energy Rep.</source> <volume>7</volume>, <fpage>224</fpage>&#x2013;<lpage>232</lpage>. <pub-id pub-id-type="doi">10.1016/j.egyr.2021.08.052</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fuqiang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Shihua</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xiaohui</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Pengcheng</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Application of ssa-bp neural network in uav point cloud hole repair</article-title>. <source>Bull. Surv. Mapp.</source>, <fpage>130</fpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Gu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). &#x201c;<article-title>Power system transient stability assessment based on svm-smote and one-dimensional convolutional neural network</article-title>,&#x201d; in <conf-name>2024 IEEE 2nd international conference on power science and Technology (ICPST) (IEEE)</conf-name>, <fpage>562</fpage>&#x2013;<lpage>568</lpage>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Qureshi</surname>
<given-names>N. M. F.</given-names>
</name>
<name>
<surname>Jarwar</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Shin</surname>
<given-names>D. R.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Ai-oriented smart power system transient stability: the rationality, applications, challenges and future opportunities</article-title>. <source>Sustain. Energy Technol. Assessments</source> <volume>56</volume>, <fpage>102990</fpage>. <pub-id pub-id-type="doi">10.1016/j.seta.2022.102990</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ji</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Power system transient voltage stability assessment based on deep residual network</article-title>. <source>Power Syst. Technol.</source> <volume>46</volume>, <fpage>2500</fpage>&#x2013;<lpage>2511</lpage>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jin</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Althubiti</surname>
<given-names>S. A.</given-names>
</name>
<name>
<surname>Alsenani</surname>
<given-names>T. R.</given-names>
</name>
<name>
<surname>Ghoneim</surname>
<given-names>M. E.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Transient stability assessment of power systems using support vector regressor and convolution neural network</article-title>. <source>Sustain. Comput. Inf. Syst.</source> <volume>37</volume>, <fpage>100826</fpage>. <pub-id pub-id-type="doi">10.1016/j.suscom.2022.100826</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kesici</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mahdi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yaslan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Seker</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Genc</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A novel integrated framework for real-time prediction of transient instabilities in power systems with their severity levels using deep learning</article-title>. <source>Sustain. Energy, Grids Netw.</source> <volume>36</volume>, <fpage>101221</fpage>. <pub-id pub-id-type="doi">10.1016/j.segan.2023.101221</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Chung</surname>
<given-names>S.-H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Power system transient stability assessment using convolutional neural network and saliency map</article-title>. <source>Energies</source> <volume>16</volume>, <fpage>7743</fpage>. <pub-id pub-id-type="doi">10.3390/en16237743</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Integrated assessment method for transient stability of power system under sample imbalance</article-title>. <source>Automation Electr. Power Syst.</source> <volume>45</volume>, <fpage>34</fpage>&#x2013;<lpage>41</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Z. Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Deep learning based on transformer architecture for power system short-term voltage stability assessment with class imbalance</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>189</volume>, <fpage>113913</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2023.113913</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A transient stability assessment model based on fault severity assignment</article-title>. <source>Front. Energy Res.</source> <volume>10</volume>, <fpage>822729</fpage>. <pub-id pub-id-type="doi">10.3389/fenrg.2022.822729</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>An automated and interpretable machine learning scheme for power system transient stability assessment</article-title>. <source>Energies</source> <volume>16</volume>, <fpage>1956</fpage>. <pub-id pub-id-type="doi">10.3390/en16041956</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Meridji</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Jo&#xf3;s</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Restrepo</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A power system stability assessment framework using machine-learning</article-title>. <source>Electr. Power Syst. Res.</source> <volume>216</volume>, <fpage>108981</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2022.108981</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>You</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>A novel data-driven lstm-saf model for power systems transient stability assessment</article-title>. <source>IEEE Trans. Industrial Inf.</source> <volume>20</volume>, <fpage>9083</fpage>&#x2013;<lpage>9097</lpage>. <pub-id pub-id-type="doi">10.1109/tii.2024.3379629</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shen</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Em-tsa: an ensemble machine learning-based transient stability assessment approach for operation of power systems</article-title>. <source>Math. Biosci. Eng.</source> <volume>20</volume>, <fpage>8226</fpage>&#x2013;<lpage>8240</lpage>. <pub-id pub-id-type="doi">10.3934/mbe.2023358</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Singh</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Chauhan</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A hybrid-extreme learning machine based ensemble method for online dynamic security assessment of power systems</article-title>. <source>Electr. Power Syst. Res.</source> <volume>214</volume>, <fpage>108923</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2022.108923</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023a</year>). <article-title>Svm based imbalanced correction method for power systems transient stability evaluation</article-title>. <source>ISA Trans.</source> <volume>136</volume>, <fpage>245</fpage>&#x2013;<lpage>253</lpage>. <pub-id pub-id-type="doi">10.1016/j.isatra.2022.10.039</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Transient stability assessment model with improved cost-sensitive method based on the fault severity</article-title>. <source>IET Generation, Transm. and Distribution</source> <volume>14</volume>, <fpage>4605</fpage>&#x2013;<lpage>4611</lpage>. <pub-id pub-id-type="doi">10.1049/iet-gtd.2020.0967</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2023b</year>). <article-title>Risk assessment of power imbalance for power systems with wind power integration considering governor ramp rate of conventional units</article-title>. <source>Electr. Power Syst. Res.</source> <volume>217</volume>, <fpage>109111</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2022.109111</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei Zhang</surname>
<given-names>Q. Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Spatial-temporal resilience assessment of distribution systems under typhoon coupled with rainstorm events</article-title>. <source>IEEE Trans. Industrial Inf</source>. <pub-id pub-id-type="doi">10.1109/TII.2024.3450079</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Two-stage transient stability prediction method of power system considering cost of misdetection and false alarm</article-title>. <source>Automation Electr. Power Syst.</source> <volume>44</volume>, <fpage>44</fpage>&#x2013;<lpage>56</lpage>.</citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Elm model for power system transient stability assessment</article-title>,&#x201d; in <source>2017 Chinese automation congress (CAC)</source> (<publisher-name>IEEE</publisher-name>), <fpage>5740</fpage>&#x2013;<lpage>5744</lpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>On-line transient stability assessment of a power system based on bagging ensemble learning</article-title>. <source>Power Syst. Prot. control</source> <volume>50</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Integrated data-driven power system transient stability monitoring and enhancement</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>39</volume>, <fpage>1797</fpage>&#x2013;<lpage>1809</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2023.3266387</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>