<?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">1666514</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2025.1666514</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>Positive sequence reactive current differential protection of transmission lines connected to energy storage power station</article-title>
<alt-title alt-title-type="left-running-head">Zhu 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.2025.1666514">10.3389/fenrg.2025.1666514</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhu</surname>
<given-names>Haoyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3135442/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Xiaoran</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Fan</surname>
<given-names>Chunju</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Hu</surname>
<given-names>Yan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tai</surname>
<given-names>Nengling</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1113089/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Smart Energy, University of Shanghai Jiao Tong</institution>, <addr-line>Shanghai</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Schoool of Electronic Information and Electrical engineering, University of Shanghai Jiao Tong</institution>, <addr-line>Shanghai</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/2286699/overview">Muhammad Waseem</ext-link>, Maynooth University, Ireland</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/1079788/overview">M. Asim Amin</ext-link>, University of Genoa, Italy</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3178817/overview">Muhammad Saqib Ali</ext-link>, Zhejiang University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Haoyu Zhu, <email>zhy119031910026@sjtu.edu.cn</email>; Chunju Fan, <email>fanchunju@sjtu.edu.cn</email>; Yan Hu, <email>yanhu@sjtu.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>30</day>
<month>10</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>13</volume>
<elocation-id>1666514</elocation-id>
<history>
<date date-type="received">
<day>15</day>
<month>07</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>09</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Zhu, Wang, Fan, Hu and Tai.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Zhu, Wang, Fan, Hu and Tai</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>When the energy storage power station encounters a fault on the transmission line during charging, active component of its short-circuit current still maintains an inverse relationship with the positive-sequence voltage at its grid connection point, influenced by the converter control strategy. This leads to a large phase difference between the short-circuit currents on both sides of the transmission line, posing a risk of no-trip failure in conventional current differential protection. To address the above issues, this paper proposes a differential protection scheme for transmission line connected to energy storage power stations based on positive-sequence reactive current, which can effectively avoid the influence of energy storage charging and discharging state on the differential current protection. The feasibility of the positive-sequence reactive current differential protection for transmission line connected to energy storage power station is analyzed through theoretical derivation. To address the issue of protection sensitivity being affected by line capacitive current when the fault voltage drop is relatively low, capacitive current compensation is added to the positive-sequence reactive current differential protection criterion. Finally, performance testing was conducted through PSCAD simulation. Results show that the proposed method can eliminate the impact of energy storage charge and discharge differences on the current differential protection performance and has good performance under different fault conditions.</p>
</abstract>
<kwd-group>
<kwd>energy storage</kwd>
<kwd>differential current protection</kwd>
<kwd>reactive current</kwd>
<kwd>charging and discharging state</kwd>
<kwd>transmission line</kwd>
</kwd-group>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Sustainable Energy Systems</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>In recent years, electrochemical energy storage in various application forms has developed rapidly (<xref ref-type="bibr" rid="B1">Amin et al., 2023</xref>; <xref ref-type="bibr" rid="B11">Mohan et al., 2024</xref>; <xref ref-type="bibr" rid="B16">Tang, 2024</xref>). Among them, for grid-connected energy storage power stations, the Chinese national standards stipulate its dynamic reactive power support capability fault crossing, but do not stipulate the magnitude and direction of the active component of the short-circuit current at this time (<xref ref-type="bibr" rid="B14">National Standards of People&#x2019;s Republic of China, 2023</xref>). Therefore, in current engineering practice, when energy storage power station encounters a fault on the transmission line during charging, its short-circuit current is characterized as injecting positive-sequence reactive current into the grid connection point and absorbing positive-sequence active current. Such short-circuit current characteristics will cause a large phase difference between the currents on both sides of the transmission line, resulting in a high risk of maloperation for conventional current differential protection (CDP).</p>
<p>Some studies have proposed protection improvement methods for the issue of the CDP no-trip failure caused by excessive phase difference. The improvement methods proposed in references (<xref ref-type="bibr" rid="B5">Guo et al., 2025</xref>; <xref ref-type="bibr" rid="B8">Lan et al., 2023</xref>; <xref ref-type="bibr" rid="B19">Zang et al., 2022</xref>) are based on the fault characteristics of photovoltaic and direct-drive wind turbines, and cannot adapt to the energy storage in charging state. Among them, reference (<xref ref-type="bibr" rid="B5">Guo et al., 2025</xref>) proposed a new principle of differential protection combining restraining current and restraining voltage, where the current restraining criterion can improve protection sensitivity compared to conventional criteria, the analysis in the paper was only conducted when the current phase difference was below 120&#xb0;, which does not cover the range of current phase difference variations during energy storage charging. Reference (<xref ref-type="bibr" rid="B8">Lan et al., 2023</xref>) constructed a new differential protection principle by combining different phase currents on both sides of the line, but this method is based on the premise that the short-circuit current on the inverter-based power side is three-phase symmetrical. Reference (<xref ref-type="bibr" rid="B19">Zang et al., 2022</xref>) proposes a d-axis-based current differential protection scheme, but this scheme cannot adapt to the characteristics that energy storage absorb active current during charging.</p>
<p>The improvement methods proposed in references (<xref ref-type="bibr" rid="B7">Jia, 2022</xref>; <xref ref-type="bibr" rid="B21">Zhang et al., 2024</xref>; <xref ref-type="bibr" rid="B9">Liang et al., 2023</xref>; <xref ref-type="bibr" rid="B18">Zang et al., 2021</xref>; <xref ref-type="bibr" rid="B10">Mishra et al., 2022</xref>) can enhance the performance of current differential protection of the transmission line connected to energy storage power station under charging conditions, but there are still shortcomings. Among them, reference (<xref ref-type="bibr" rid="B7">Jia, 2022</xref>) proposes a current differential protection based on amplitude comparison, by utilizing the significant difference in short-circuit current amplitude between inverter-based power and the system. However, this method still faces the problem of insufficient sensitivity when the capacity of new energy station is large. Reference (<xref ref-type="bibr" rid="B21">Zhang et al., 2024</xref>) proposes an adaptive current differential improvement criterion that compensates for both current amplitude and phase for energy storage current differential protection. However, the constant coefficients in the amplitude compensation function and phase compensation function mentioned in the paper are not explained in terms of their meaning and selection principles, and there is no analysis or simulation on whether this criterion will cause protection maloperation during external faults. Reference (<xref ref-type="bibr" rid="B9">Liang et al., 2023</xref>) proposes an improved method for energy storage current differential protection. This method compensates both the magnitude and phase of the short-circuit current to improve protection sensitivity. However, selection principles of constant coefficients in the magnitude compensation function and phase compensation function mentioned in the paper are not explained, and there is no analysis or simulation on whether this criterion will cause protection maloperation during external faults. Reference (<xref ref-type="bibr" rid="B18">Zang et al., 2021</xref>) proposes an enhanced current differential protection that modifies the amplitude ratio and phase difference of currents on both sides of the line. This method can effectively improve the performance of current differential protection, but due to the need to ensure reliability during external faults, it still lacks sufficient sensitivity for internal faults when the current amplitude ratio is close to 1. Reference (<xref ref-type="bibr" rid="B10">Mishra et al., 2022</xref>) proposes a distribution network current differential protection scheme based on Q-axis current, but it does not consider the effect of line capacitive current, which may lead to insufficient sensitivity under conditions of high transition resistance. Moreover, the photovoltaic power capacity in its example is very small, which differs significantly from scenarios with large-capacity energy storage power stations accessing the grid. Reference (<xref ref-type="bibr" rid="B4">Farshad, 2021</xref>) uses the second harmonic components of Q-axis current and voltage to achieve protection of distribution networks containing photovoltaic power stations, but the paper does not address the protection performance when the inverter has harmonic suppression control strategies.</p>
<p>In addition, there are also some methods using indicators such as cosine similarity (<xref ref-type="bibr" rid="B15">Sirisha and Pradhan, 2020</xref>; <xref ref-type="bibr" rid="B22">Zheng et al., 2021</xref>; <xref ref-type="bibr" rid="B20">Zhang et al., 2019</xref>) and Pearson correlation coefficient to identify faults (<xref ref-type="bibr" rid="B3">Chen et al., 2018</xref>)- (<xref ref-type="bibr" rid="B6">Jia et al., 2018</xref>). These methods use a data-driven approach to avoid conventional protection performance issues caused by renewable energy sources, but outliers may have a significant impact on them.</p>
<p>In this paper, a positive-sequence reactive current differential protection suitable for transmission line connected to energy storage power station is proposed based on the fault characteristics of energy storage injecting.</p>
<p>Positive-sequence reactive current into the grid connection point during transmission line fault. This method first calculates the positive-sequence reactive current component on the local side using the positive-sequence voltage and positive-sequence current at the protection installation. Subsequently, based on the line capacitance parameters and bus voltage, the positive-sequence reactive current is compensated. The compensated positive-sequence reactive current component is used for longitudinal differential protection. This method can eliminate the phase difference caused by the absorption of active current during energy storage charging, avoiding the risk of no-trip failure.</p>
</sec>
<sec id="s2">
<title>2 Effect of energy storage charging on the CDP</title>
<sec id="s2-1">
<title>2.1 Energy storage converter control strategy</title>
<p>For the current grid-connected energy storage, its converter control strategy usually adopts a dual closed-loop control method with a power outer loop and a current inner loop (<xref ref-type="bibr" rid="B17">Telukunta et al., 2017</xref>; <xref ref-type="bibr" rid="B2">Chen and Mei, 2015</xref>; <xref ref-type="bibr" rid="B13">National Standards of People&#x2019;s Republic of China, 2021</xref>; <xref ref-type="bibr" rid="B12">National Standards of People&#x2019;s Republi c of China, 2024</xref>). When different faults occur on the transmission line connected to energy storage power station, in order to meet the dynamic reactive power support capability specified by the current national standards and the current limiting requirements of the converter (<xref ref-type="bibr" rid="B14">National Standards of People&#x2019;s Republic of China, 2023</xref>), the control targets for the positive sequence reactive and active components of the energy storage short-circuit current are shown in <xref ref-type="disp-formula" rid="e1">Equations 1</xref>, <xref ref-type="disp-formula" rid="e2">2</xref>.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mi>i</mml:mi>
<mml:mi>q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0.9</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>0.9</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0.9</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0.9</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>max</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:msubsup>
<mml:mi>i</mml:mi>
<mml:mi>q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0.9</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf1">
<mml:math id="m3">
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>C</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the positive sequence voltage at the energy storage and network point; <inline-formula id="inf2">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the rated current of the energy storage; <inline-formula id="inf3">
<mml:math id="m5">
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mi>i</mml:mi>
<mml:mi>q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf4">
<mml:math id="m6">
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> are the positive-sequence reactive current and active current, respectively. <inline-formula id="inf5">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the dynamic reactive current coefficient, with a value of 1.5-3, specified by reference (<xref ref-type="bibr" rid="B14">National Standards of People&#x2019;s Republic of China, 2023</xref>). When K1 is smaller, the reactive current is smaller, making the reactive current differential protection least likely to operate; Therefore, this paper chooses 1.5, the most unfavorable case, to test the performance of the proposed scheme; <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the converter amplitude limiting requirement, this paper takes <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:mn>1.5</mml:mn>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
<sec id="s2-2">
<title>2.2 Effect of active component of short-circuit current on the CDP</title>
<p>Influenced by the control strategy, no matter in the charging or discharging state, the direction of the positive-sequence reactive power component of the short-circuit current is always from the energy storage to the grid connection point when a short-circuit occurs in the transmission line connected to energy storage power station. However, the direction of its positive-sequence active current is directly related to the charging and discharging state of the energy storage.</p>
<p>
<xref ref-type="fig" rid="F1">Figure 1</xref> shows the different relationship under charging and discharging conditions between short-circuit current of energy storage and voltage of the grid connection point. <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the positive-sequence active components of energy storage short-circuit current under charging and discharging conditions, respectively. The directions of <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are opposite during charging and discharging states, resulting in a significant phase difference between phase currents <inline-formula id="inf12">
<mml:math id="m14">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf13">
<mml:math id="m15">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Characteristics of short-circuit current of energy storage.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g001.tif">
<alt-text content-type="machine-generated">Diagram illustrating phase and current vectors in a coordinate system. Vectors labeled \(I^&#x2b;_{d\_cha}\), \(I^&#x2b;_{d\_disc}\), \(I^&#x2b;_{phase\_cha}\), \(I^&#x2b;_{phase\_disc}\), \(I^&#x2b;_{q\_disc}\), \(I^&#x2b;_{q\_cha}\), and \(U^&#x2b;_{PCC}\) are shown with directional arrows, intersecting at an origin point.</alt-text>
</graphic>
</fig>
<p>From <xref ref-type="fig" rid="F1">Figure 1</xref>, it can be seen that the phase difference between <inline-formula id="inf14">
<mml:math id="m16">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf15">
<mml:math id="m17">
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is an obtuse angle. This phase characteristic of short-circuit current under charging conditions will significantly increase the phase difference of short-circuit currents on both sides of the line, thereby leading to the risk of no-trip failure of the CDP. This risk will be more obvious when the capacity of the energy storage plant increases. The performance problem of the CDP caused by rectified state of converter has been analyzed in detail in the literature (<xref ref-type="bibr" rid="B9">Liang et al., 2023</xref>).</p>
</sec>
</sec>
<sec id="s3">
<title>3 Modified current differential protection based on adaptive phase compensation and its setting methods</title>
<p>According to the analysis in <xref ref-type="sec" rid="s2">Section 2</xref>, it can be seen that no-trip failure of the CDP under the energy storage charging state is mainly caused by the active component of the short-circuit current provided by the energy storage power station. Therefore, this paper proposes a differential protection based on the positive-sequence reactive component of the short-circuit current to avoid the influence of the active component.</p>
<sec id="s3-1">
<title>3.1 Energy storage converter control strategy</title>
<p>In the system shown in <xref ref-type="fig" rid="F2">Figure 2</xref>, the feasibility of applying positive-sequence reactive current differential protection to the transmission line PQ is analyzed. The positive-sequence reactive currents at points P and Q are both calculated using the local voltage and current through d-q decomposition. Since only the positive-sequence reactive current is used, performance of the positive-sequence reactive current differential protection applied to the line PQ is basically unaffected by the difference between charging and discharging states. A feasibility analysis is conducted using the charging state as an example.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Characteristics of short-circuit current of energy storage.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g002.tif">
<alt-text content-type="machine-generated">Diagram showing an energy storage system connected to a main grid. Points P and Q are along the line, with switches k1 and k2 between them. The flow directs towards the main grid on the right.</alt-text>
</graphic>
</fig>
<sec id="s3-1-1">
<title>3.1.1 Fault in the line PQ</title>
<p>When a fault occurs at a point k1 on the line PQ in <xref ref-type="fig" rid="F2">Figure 2</xref>, the short-circuit current at bus Q is provided by the system side. At this time, the direction of the positive-sequence reactive current calculated from the positive-sequence voltage and positive-sequence current must be from bus Q to the fault point k1; bus P is the energy storage grid connection point. According to the fault ride-through standard, regardless of whether the energy storage is charging or discharging at this time, the direction of its positive-sequence reactive current is from bus P to the fault point k1.</p>
</sec>
<sec id="s3-1-2">
<title>3.1.2 Fault out of the line PQ</title>
<p>When the fault occurs at a point k2 on line QN in <xref ref-type="fig" rid="F2">Figure 2</xref>, the direction of the positive sequence reactive current at bus P is still bus-pointing to the point of fault.</p>
<p>Assuming that the steady state of the control target can be reached after the fault, the voltage-current relationship at the energy storage grid connection point P can be written, resulting in <xref ref-type="disp-formula" rid="e3">Equation 3</xref>.<disp-formula id="e3">
<mml:math id="m18">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf16">
<mml:math id="m19">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf17">
<mml:math id="m20">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the positive-sequence current and positive-sequence voltage at point P, respectively; <inline-formula id="inf18">
<mml:math id="m21">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the positive-sequence voltage at fault point k2; <inline-formula id="inf19">
<mml:math id="m22">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the positive-sequence voltage drop from point P to fault point k2; <inline-formula id="inf20">
<mml:math id="m23">
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the positive-sequence impedance from point P to fault point k2. Using the phase of <inline-formula id="inf21">
<mml:math id="m24">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> as the reference direction, draw the phasor diagram as shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Phasor diagram of transmission line fault during energy storage charging.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g003.tif">
<alt-text content-type="machine-generated">Geometric diagram illustrating vectors and angles. Points labeled \(I_P^&#x2b; &#x3d; I_Q^&#x2b;\), \(U_Q^&#x2b;\), \(U_P^&#x2b;\), \(U_{k2}^&#x2b;\), and \(\Delta U_{P-k2}^&#x2b;\) show vector relationships. Angles labeled \(\theta_l\) and LVRT indicate measurements. Arrows depict vector direction.</alt-text>
</graphic>
</fig>
<p>In <xref ref-type="fig" rid="F3">Figure 3</xref>, &#x2220;LVRT is the angle between the positive-sequence voltage and positive-sequence current at the energy storage grid connection point, determined by the energy storage control strategy; <inline-formula id="inf22">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the line impedance angle; the positive-sequence voltage <inline-formula id="inf23">
<mml:math id="m26">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> at bus Q is located between <inline-formula id="inf24">
<mml:math id="m27">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf25">
<mml:math id="m28">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>When the fault point is near point Q, <inline-formula id="inf26">
<mml:math id="m29">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf27">
<mml:math id="m30">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are basically the same. At this time, the phase lead-lag relationship between <inline-formula id="inf28">
<mml:math id="m31">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf29">
<mml:math id="m32">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> is most likely to change, draw the critical condition of <inline-formula id="inf30">
<mml:math id="m33">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf31">
<mml:math id="m34">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> parallel, as shown in <xref ref-type="fig" rid="F4">Figure 4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Phasor diagram of the critical condition.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g004.tif">
<alt-text content-type="machine-generated">Diagram illustrating a phasor representation involving vectors \( I_{Q}^&#x2b; \), \( i_{q}^&#x2b; \), \( i_{d}^&#x2b; \), \( U_{P}^&#x2b; \), and \( U_{k2}^&#x2b; \). Angles noted include \( \angle LVRT \) and \( \theta_l \), with \( \Delta U_{P-k2}^&#x2b; \) and corresponding angle relationships \( \pi - \angle LVRT \) and \( \theta_1 \).</alt-text>
</graphic>
</fig>
<p>From the geometric relationship shown in <xref ref-type="fig" rid="F4">Figure 4</xref>, it can be obtained that when <xref ref-type="disp-formula" rid="e4">Equation 4</xref> holds, <inline-formula id="inf32">
<mml:math id="m35">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> always leads <inline-formula id="inf33">
<mml:math id="m36">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>.<disp-formula id="e4">
<mml:math id="m37">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2264;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mstyle>
</mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>&#x3c0;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2220;</mml:mo>
<mml:mtext>LVRT</mml:mtext>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>In(4), <inline-formula id="inf34">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the line impedance angle, and <inline-formula id="inf35">
<mml:math id="m39">
<mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is taken as 1 according to the most unfavorable condition.</p>
<p>When the fault type is a metallic three-phase short circuit, the positive sequence voltage at the energy storage grid connection point is too low to maintain a stable charging state.</p>
<p>If the fault type is other metallic short circuit, <inline-formula id="inf36">
<mml:math id="m40">
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> reaches the limit value <inline-formula id="inf37">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and the value of <inline-formula id="inf38">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is taken as <inline-formula id="inf39">
<mml:math id="m43">
<mml:mrow>
<mml:mn>1.5</mml:mn>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> in this paper. Combining with <xref ref-type="disp-formula" rid="e1">Equations 1</xref>, <xref ref-type="disp-formula" rid="e4">4</xref> can be further simplified to <xref ref-type="disp-formula" rid="e5">Equation 5</xref>.<disp-formula id="e5">
<mml:math id="m44">
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>0.9</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>2.25</mml:mn>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e5">Equation 5</xref>, <inline-formula id="inf40">
<mml:math id="m45">
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the system potential. <inline-formula id="inf41">
<mml:math id="m46">
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is basically the same as line voltage drop during normal operation, usually less than 10% E, significantly smaller than the value on the right side of <xref ref-type="disp-formula" rid="e5">Equation 5</xref>. Therefore, at this time, the direction of the reactive current at point Q is from bus Q to the fault point k2, and the differential protection of line PQ will not operate incorrectly.</p>
</sec>
</sec>
<sec id="s3-2">
<title>3.2 Positive-sequence reactive current differential protection criterion</title>
<p>Based on the analysis results in Section A, the positive-sequence reactive current differential protection criterion (RCDP) can be designed as shown in <xref ref-type="disp-formula" rid="e6">Equation 6</xref>.<disp-formula id="e6">
<mml:math id="m47">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>In(6), <inline-formula id="inf42">
<mml:math id="m48">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf43">
<mml:math id="m49">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> are the positive-sequence reactive currents measured at points P and Q in <xref ref-type="fig" rid="F2">Figure 2</xref>, with the positive direction defined as from the bus to the line. <inline-formula id="inf44">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the operating current, set according to <xref ref-type="disp-formula" rid="e7">Equation 7</xref>.<disp-formula id="e7">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.1</mml:mn>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">k</mml:mi>
<mml:mo>.</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e7">Equation 7</xref>, the coefficient 0.1 indicates the error of the current transformer; <inline-formula id="inf45">
<mml:math id="m52">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">k</mml:mi>
<mml:mo>.</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the maximum positive-sequence reactive short-circuit current during an external fault; <inline-formula id="inf46">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the coefficient to measure whether equipment is of the same type, taken as 0.5 in this paper; <inline-formula id="inf47">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the non-periodic component coefficient, used to reflect the effect of non-periodic components, taken as 1.5 in this paper.</p>
</sec>
<sec id="s3-3">
<title>3.3 Performance analysis of the RCDP</title>
<sec id="s3-3-1">
<title>3.3.1 The effect of line capacitance</title>
<p>Feasibility of the RCDP has been analyzed in Section A in Chapter <sc>III</sc> based on short-circuit current characteristics of energy storage, but the effect of line distributed capacitance on protection performance under special fault conditions was not considered in the analysis process. In fact, as the transition resistance at the short-circuit point increases, the positive-sequence voltage at the energy storage grid connection point gradually rises to above 0.9 p. u. During this process, the positive-sequence reactive current provided by the energy storage power station gradually decreases to zero. For the system side, if the transition resistance continues to increase, the positive-sequence reactive component of its short-circuit current will also be reduced to the same order of magnitude as the line capacitance current, and thus the effect of line capacitance current on the RCDP is not negligible.</p>
<p>Since point P is the energy storage grid connection point, the reactive current at point P should be determined solely by the control strategy. The positive-sequence reactive current component at Q in <xref ref-type="disp-formula" rid="e8">Equation 8</xref> should be the superposition of the positive-sequence reactive component of the line capacitive current and the reactive component of the short-circuit current when the line capacitive effect is neglected. The positive sequence component of the line capacitance current at Q is noted to be <inline-formula id="inf48">
<mml:math id="m55">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>, and its direction is perpendicular to <inline-formula id="inf49">
<mml:math id="m56">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>.The short-circuit currents at P and Q when the effect of line capacitance is neglected are <inline-formula id="inf50">
<mml:math id="m57">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>s</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf51">
<mml:math id="m58">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>s</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>, respectively. <xref ref-type="disp-formula" rid="e6">Equation 6</xref> can be rewritten as <xref ref-type="disp-formula" rid="e8">Equation 8</xref>.<disp-formula id="e8">
<mml:math id="m59">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>s</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>s</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mover accent="true">
<mml:mi>I</mml:mi>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>In expression <xref ref-type="disp-formula" rid="e8">Equation 8</xref>, the left-side expression is the differential current, and the ratio of the differential current to the operating current <inline-formula id="inf52">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the protection sensitivity <inline-formula id="inf53">
<mml:math id="m61">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.The positive sequence component of the line capacitor current flows from the energy storage to the system. <inline-formula id="inf54">
<mml:math id="m62">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf55">
<mml:math id="m63">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>s</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> are in opposite directions. The presence of <inline-formula id="inf56">
<mml:math id="m64">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> will make the differential current smaller thus leading to a reduction in the sensitivity of the protection. Therefore, capacitive current compensation needs to be added to the RCDP.</p>
</sec>
<sec id="s3-3-2">
<title>3.3.2 The effect of operating current setting methods</title>
<p>Operating current <inline-formula id="inf57">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> in <xref ref-type="disp-formula" rid="e6">Equation 6</xref> is set to avoid the maximum positive-sequence reactive power unbalance current in case of out-of-area fault. However, magnitude of the positive-sequence reactive component of the short-circuit current varies significantly with different fault types and transition resistances, and sensitivity of the RCDP will be low if it is set only on the basis of the imbalance current in the most severe case. Therefore, it is necessary to construct the positive-sequence reactive current differential protection criterion with reference to the conventional ratio-braking current differential protection.</p>
</sec>
</sec>
<sec id="s3-4">
<title>3.4 Positive-sequence reactive current differential protection criterion considering line capacitance current compensation</title>
<p>Based on the above analysis, the positive-sequence reactive current differential protection considering line capacitive current compensation (CRCDP) is constructed as shown in <xref ref-type="disp-formula" rid="e9">Equation 9</xref>.<disp-formula id="e9">
<mml:math id="m66">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3e;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>max</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>w</mml:mi>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>Q</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>In (9), <inline-formula id="inf58">
<mml:math id="m67">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msup>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the compensated positive sequence reactive current at Q. The positive direction of both <inline-formula id="inf59">
<mml:math id="m68">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf60">
<mml:math id="m69">
<mml:mrow>
<mml:msup>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is defined as from the bus to the line. The calculation is performed according to (<xref ref-type="disp-formula" rid="e10">Equation 10</xref>), where <inline-formula id="inf61">
<mml:math id="m70">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>Q</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the capacitance of line PQ, <inline-formula id="inf62">
<mml:math id="m71">
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the angular frequency, and <inline-formula id="inf63">
<mml:math id="m72">
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is the voltage magnitude at Q; <inline-formula id="inf64">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the restraining coefficient, which is the same as the slope in the conventional percentage differential protection.<disp-formula id="e10">
<mml:math id="m74">
<mml:mrow>
<mml:msup>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>w</mml:mi>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>Q</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>Noting that the left side of <xref ref-type="disp-formula" rid="e9">Equation 9</xref> is the differential amount, the right side is the restraining amount, and ratio of the restraining amount to the differential amount is protection sensitivity. <inline-formula id="inf65">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the restraining current set with reference to the conventional ratio-braking type current differential protection. Since the voltage at the Q bus is used for capacitive current compensation on the Q side, rather than an integral calculation utilizing the voltage along the line, this may result in an error in the restraining amount. To avoid the restraining amount being too low due to capacitor current compensation errors, the maximum compensation error <inline-formula id="inf66">
<mml:math id="m76">
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>Q</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mover accent="true">
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>Q</mml:mi>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mo>&#xb7;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> needs to be set as the lower limit of the restraining amount to prevent protection maloperation.</p>
<p>Protection action logic of the CRCDP is shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. The protection action signal needs to remain effective for more than 5 m to avoid maloperation caused by transient fluctuations.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Protection action logic of the CRCDP.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g005.tif">
<alt-text content-type="machine-generated">Flowchart illustrating a decision-making process. It starts with a condition \(U^&#x2b;_{P-pu} \geq 0.1\). If &#x22;No,&#x22; it leads to &#x22;Protection blocking.&#x22; If &#x22;Yes,&#x22; it calculates \(I_{Q^&#x2b;}'&#x3d;I_{Q^&#x2b;} - wc_Q &#x7c;U^&#x2b;_{\Theta}&#x7c;\). Satisfying Equation (9) moves to checking if &#x22;Signal persists over 5 ms.&#x22; If this is satisfied, it results in &#x22;Action&#x22;; if not, it results in &#x22;Inaction.&#x22; &#x22;Disatisfy&#x22; cases lead to &#x22;Protection blocking.&#x22;</alt-text>
</graphic>
</fig>
</sec>
</sec>
<sec id="s4">
<title>4 Simulation verification</title>
<p>Build the simulation model according to the system structure shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. The equivalent impedance of the external system and line parameters are shown in <xref ref-type="table" rid="T1">Table 1</xref>. The lengths of lines PQ and QN are both 40 km, and the capacity of energy storage power station is 200 MW. The transformer has a voltage rating of 35kV/220 kV and a rated capacity of 300MVA. The control strategy and filter of the energy storage converter is shown in <xref ref-type="fig" rid="F7">Figures 7a,b</xref>. In <xref ref-type="fig" rid="F7">Figure 7b</xref>, <inline-formula id="inf67">
<mml:math id="m77">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is 1400 <inline-formula id="inf68">
<mml:math id="m78">
<mml:mrow>
<mml:mi mathvariant="normal">&#x3bc;</mml:mi>
<mml:mi mathvariant="normal">F</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf69">
<mml:math id="m79">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is 4.8 <inline-formula id="inf70">
<mml:math id="m80">
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf71">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is 0.000675 <inline-formula id="inf72">
<mml:math id="m82">
<mml:mrow>
<mml:mi mathvariant="normal">H</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula id="inf73">
<mml:math id="m83">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is 157 <inline-formula id="inf74">
<mml:math id="m84">
<mml:mrow>
<mml:mi mathvariant="normal">&#x3bc;</mml:mi>
<mml:mi mathvariant="normal">F</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. In the simulation, the positive sequence components of voltage and current are extracted using the frequency scanner module integrated in PSCAD 4.6. This module employs a Fourier filter algorithm.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Structure of the grid-connected system of energy storage power station <bold>(a)</bold> Control strategy <bold>(b)</bold> Filter.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g006.tif">
<alt-text content-type="machine-generated">Diagram showing an energy storage system connected to a main grid. The layout includes connection points P, Q, and N with series devices labeled F1 to F6. A breaker symbol and energy flow indicator are shown.</alt-text>
</graphic>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>
<sc>L</sc>ine parameters and system equivalent impedance.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Line positive sequence impedance</th>
<th align="left">0.074 &#x2b; j0.42 (&#x3a9;/km)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Line zero sequence impedance</td>
<td align="left">0.222 &#x2b; j1.26 (&#x3a9;/km)</td>
</tr>
<tr>
<td align="left">System Equivalent Impedance</td>
<td align="left">5 &#x2b; j31.42 (&#x3a9;)</td>
</tr>
<tr>
<td align="left">Line capacitance</td>
<td align="left">0.0605 (&#x3bc;F/km)</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Control strategy and filter of the energy storage converter. <bold>(a)</bold> Control strategy. <bold>(b)</bold> Filter.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g007.tif">
<alt-text content-type="machine-generated">Diagram (a) shows a control system with multiple PI controllers, inputs labeled \( P \), \( Q \), \( P_{\text{ref}} \), and \( Q_{\text{ref}} \), and outputs labeled \( V_d \) and \( V_q \) connected to SPWM. Diagram (b) shows an electrical circuit with a resistor \( R_{\text{damp}} \), an inductor \( L_{\text{damp}} \), and capacitors \( C_{\text{damp}} \) and \( C_{\text{filter}} \).</alt-text>
</graphic>
</fig>
<p>Points F1-F6 in <xref ref-type="fig" rid="F6">Figure 6</xref> are the fault locations to be used in the subsequent simulation, and the specific location descriptions of F1-F6 are shown in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>F1-F6 position description.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Position</th>
<th align="left">Explanation</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">F1</td>
<td align="left">5% away from bus P on the line PQ</td>
</tr>
<tr>
<td align="left">F2</td>
<td align="left">20% away from bus P on the line PQ</td>
</tr>
<tr>
<td align="left">F3</td>
<td align="left">Middle of the line PQ</td>
</tr>
<tr>
<td align="left">F4</td>
<td align="left">80% away from bus P on the line PQ</td>
</tr>
<tr>
<td align="left">F5</td>
<td align="left">95% away from bus P on the line PQ</td>
</tr>
<tr>
<td align="left">F6</td>
<td align="left">5% away from bus Q on the line QN</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Points F1-F6 in <xref ref-type="fig" rid="F6">Figure 6</xref> <xref ref-type="table" rid="T2">Table 2</xref> are the fault locations to be used in the subsequent simulation, and the specific location descriptions of F1-F6 are shown in.</p>
<sec id="s4-1">
<title>4.1 Tripping time comparison</title>
<p>In the simulation, the energy storage is set to operate in charging mode, and a fault occurs at the midpoint of the line PQ (F3 in <xref ref-type="fig" rid="F6">Figure 6</xref>) at 1.5 s. The fault type is set as metallic phase A ground fault. Operating curves of CDP, RCDP, and CRCDP is shown in <xref ref-type="fig" rid="F8">Figure 8</xref>.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Operating curves of CDP, RCDP, and CRCDP. <bold>(a)</bold> Trip signal of CDP. <bold>(b)</bold> Trip signal of RCDP. <bold>(c)</bold> Trip signal of CRCDP.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g008.tif">
<alt-text content-type="machine-generated">Three graphs show trip signals over time. Graph (a) has a black line indicating 17.05 milliseconds for tripping. Graph (b) shows a blue line with 8.82 milliseconds. Graph (c) features a red line with 8.76 milliseconds. All x-axes range from 1.45 to 1.8 seconds.</alt-text>
</graphic>
</fig>
<p>As shown in <xref ref-type="fig" rid="F8">Figure 8</xref>, the CDP, RCDP, and CRCDP can all correctly identify faults during metallic phase A ground fault. Among them, CDP requires the longest tripping time of 17.05 m. The tripping time of RCDP and CRCDP is significantly shorter than that of CDP, which is 8.82 m and 8.76 m, respectively.</p>
</sec>
<sec id="s4-2">
<title>4.2 Protection performance comparison during internal and external faults</title>
<sec id="s4-2-1">
<title>4.2.1 Internal faults</title>
<p>In the simulation, the energy storage is set to operate in charging mode, and a fault occurs at the midpoint of the line PQ (F3 in <xref ref-type="fig" rid="F6">Figure 6</xref>) at 1.5 s. The fault types and their corresponding transition resistances are shown in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Fault types and transition resistance.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Case number</th>
<th align="left">Fault type</th>
<th align="left">Fault resistance</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Case1</td>
<td align="left">AG</td>
<td align="left">0&#x3a9;</td>
</tr>
<tr>
<td align="left">Case2</td>
<td align="left">AG</td>
<td align="left">100&#x3a9;</td>
</tr>
<tr>
<td align="left">Case3</td>
<td align="left">AB</td>
<td align="left">10&#x3a9;</td>
</tr>
<tr>
<td align="left">Case4</td>
<td align="left">ABG</td>
<td align="left">100&#x3a9;</td>
</tr>
<tr>
<td align="left">Case5</td>
<td align="left">ABCG</td>
<td align="left">100&#x3a9;</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The sensitivities of CDP, RCDP, and CRCDP corresponding to each fault are shown in <xref ref-type="fig" rid="F9">Figures 9a&#x2013;e</xref>. When the sensitivity exceeds 1 and remains above 5 m (as mentioned in <xref ref-type="fig" rid="F5">Figure 5</xref>), the protection trips.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Protection performance of different types of internal faults. <bold>(a)</bold> AG0ohm. <bold>(b)</bold> AG100ohm. <bold>(c)</bold> AB10ohm. <bold>(d)</bold> ABG100ohm. <bold>(e)</bold> ABCG100ohm.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g009.tif">
<alt-text content-type="machine-generated">Five line graphs labeled (a) to (e) compare Ksen over time for CDP, RCDP, and CRCDP. Each graph shows different peak values and times. CDP is blue, RCDP black, and CRCDP red, with varying behavior in each graph across the 1.45 to 1.8 seconds range.</alt-text>
</graphic>
</fig>
<p>As shown in <xref ref-type="fig" rid="F9">Figure 9</xref>, when the fault occurs during charging state, sensitivity of the CDP is seriously insufficient in faults from Case 2 to Case 5, posing a significant risk of protection no-trip failure. Sensitivity of the RCDP is very high in metallic phase A ground fault (Case 1) and AB phase-to-phase ground fault (Case 3), but when the transition resistance is high (Case 2 and Case 4), the sensitivity is only 0.63 and 0.97 respectively, also indicating insufficient sensitivity. Sensitivity of the CRCDP meets the operating conditions under various fault types.</p>
</sec>
<sec id="s4-2-2">
<title>4.2.2 External faults</title>
<p>In the simulation, the energy storage is set to operate in charging mode, and a fault occurs at 5% from bus Q on line QN (F6 in <xref ref-type="fig" rid="F6">Figure 6</xref>) at 1.5 s. The fault types and their corresponding transition resistances are shown in <xref ref-type="table" rid="T3">Table 3</xref>. The sensitivities of CDP, RCDP, and CRCDP corresponding to each fault are shown in <xref ref-type="fig" rid="F10">Figures 10a&#x2013;e</xref>.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Protection performance of different types of external faults. <bold>(a)</bold> AG0ohm. <bold>(b)</bold> AG100ohm. <bold>(c)</bold> AB10ohm. <bold>(d)</bold> ABG100ohm. <bold>(e)</bold> ABCG100ohm.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g010.tif">
<alt-text content-type="machine-generated">Five line graphs display the Kesn values over time, labeled (a) to (e). Each graph compares three methods: CDP (black line), RCDP (blue line), and CRCDP (red line). Significant fluctuations appear around 1.5 seconds, with CRCDP showing the most variation before stabilizing. The time range spans from 1.45 to 1.8 seconds on the x-axis, while Kesn values are on the y-axis.</alt-text>
</graphic>
</fig>
<p>As shown in <xref ref-type="fig" rid="F10">Figure 10</xref>, when external faults occur, sensitivity of the CRCDP may only exceed 1 during the transient process after the fault, but the duration does not exceed 5 m, which does not meet the signal duration requirement for protection trip. Thus, no maloperation occurs. The CRCDP has sufficient reliability for external faults.</p>
</sec>
</sec>
<sec id="s4-3">
<title>4.3 Effect of transition resistance on protection performance</title>
<p>In the simulation, the energy storage is set to operate in charging mode. A three-phase short circuit occurs at the midpoint of the line PQ (F3 in <xref ref-type="fig" rid="F6">Figure 6</xref>) at 1.5 s. The transition resistances are 10&#x3a9;, 50&#x3a9;, 100&#x3a9;, and 200&#x3a9; respectively. The sensitivities of CDP, RCDP, and CRCDP corresponding to each fault are shown in <xref ref-type="fig" rid="F11">Figures 11a&#x2013;e</xref>.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Comparison of protection sensitivity at different transition resistances. <bold>(a)</bold> Fault resistance: 10ohm. <bold>(b)</bold> Fault resistance: 50ohm. <bold>(c)</bold> Fault resistance: 100ohm. <bold>(d)</bold> Fault resistance: 300ohm</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g011.tif">
<alt-text content-type="machine-generated">Four graphs labeled (a) to (d) display Ksen versus time from 1.45 to 1.8 seconds. Each graph compares CDP, RCDP, and CRCDP lines in black, blue, and red, respectively. Key points marked: (a) 2.68 and 0.48, (b) 3.23, (c) 1.38, (d) 1.64 and 0.13.</alt-text>
</graphic>
</fig>
<p>As can be seen from <xref ref-type="fig" rid="F10">Figure 10</xref>, the sensitivity of the CDP is seriously insufficient, even when the transition resistance is low (10&#x3a9;), its sensitivity is only 0.48, which can not meet the action requirements, and there is a greater risk of protection no-trip failure; the RCDP has a high sensitivity when the transition resistance is low, but its ability to resist the transition resistance is poor. When the transition resistance reaches 300&#x3a9;, its sensitivity is only 0.13, there is also a greater risk of protection no-trip failure; CRCDP maintains high sensitivity throughout the changes in transition resistance. Even when the transition resistance increases to 300&#x3a9;, its sensitivity reaches 1.64, ensuring accurate protection operation.</p>
</sec>
<sec id="s4-4">
<title>4.4 Effect of fault location on protection performance</title>
<p>In the simulation, the energy storage is set to operate in the charging state, and a phase A grounded short circuit occurs at positions F1-F5 in <xref ref-type="fig" rid="F6">Figure 6</xref> at 1.5 s, with a transition resistance of 100&#x3a9;. The sensitivity and operation status of CRCDP corresponding to faults at F1-F5 are shown in <xref ref-type="table" rid="T4">Table 4</xref>.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Sensitivity and protection action status of the CRCDP for faults F1 through F5.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Position</th>
<th align="left">
<inline-formula id="inf75">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="left">Performance of CRCDP</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">F1</td>
<td align="left">1.95</td>
<td align="left">correct action</td>
</tr>
<tr>
<td align="left">F2</td>
<td align="left">1.86</td>
<td align="left">correct action</td>
</tr>
<tr>
<td align="left">F3</td>
<td align="left">1.68</td>
<td align="left">correct action</td>
</tr>
<tr>
<td align="left">F4</td>
<td align="left">1.51</td>
<td align="left">correct action</td>
</tr>
<tr>
<td align="left">F5</td>
<td align="left">1.43</td>
<td align="left">correct action</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>As can be seen from <xref ref-type="table" rid="T4">Table 4</xref>, as the fault point moves from point F1 to point F5, the sensitivity of CRCDP gradually decreases, but it can always ensure normal protection operation.</p>
</sec>
<sec id="s4-5">
<title>4.5 Effect of energy storage capacity on protection performance</title>
<p>In the simulation, the energy storage capacities are set to 100 MW, 150 MW, 200 MW, and 250 MW respectively, with the energy storage operating in charging mode. At 1.5 s, phase A grounded short circuit occurs at the midpoint of the line PQ (P3 in <xref ref-type="fig" rid="F6">Figure 6</xref>), with a transition resistance of 100&#x3a9;. The CRCDP sensitivity corresponding to each energy storage capacity is shown in <xref ref-type="fig" rid="F12">Figure 12</xref>.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>CRCDP sensitivity at different energy storage capacities.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g012.tif">
<alt-text content-type="machine-generated">Line graph showing Ksen over time from 1.45 to 1.8 seconds for power levels 100 megawatts, 150 megawatts, 200 megawatts, and 250 megawatts. The 250 megawatts line spikes and stabilizes around 3. The others stabilize near 1.</alt-text>
</graphic>
</fig>
<p>As shown in <xref ref-type="fig" rid="F12">Figure 12</xref>, the energy storage capacity has a significant impact on the sensitivity of CRCDP. The larger the energy storage capacity, the higher the sensitivity. When the energy storage capacity is 100 MW, the sensitivity of CRCDP is relatively low at about 1.21, which can still ensure normal protection operation.</p>
</sec>
<sec id="s4-6">
<title>4.6 Effect of charging and discharging states on protection performance</title>
<p>Take phase A grounded short circuit as an example, the sensitivity and action status of CRCDP under different charging and discharging states are explained. The transition resistance is set to 100&#x3a9;, and the fault start time is 1.5 s. The fault locations are set to P3 and P6 in <xref ref-type="fig" rid="F6">Figure 6</xref>, respectively. For each fault location, the energy storage is set to be in charging and discharging states, respectively. The specific fault conditions are shown in <xref ref-type="table" rid="T5">Table 5</xref>.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Fault location and storage charge/discharge status.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Case number</th>
<th align="left">Charging and discharging state</th>
<th align="left">Fault position</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Case1</td>
<td align="left">Charging</td>
<td align="left">F3</td>
</tr>
<tr>
<td align="left">Case2</td>
<td align="left">Discharging</td>
<td align="left">F3</td>
</tr>
<tr>
<td align="left">Case3</td>
<td align="left">Charging</td>
<td align="left">F6</td>
</tr>
<tr>
<td align="left">Case4</td>
<td align="left">Discharging</td>
<td align="left">F6</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The CRCDP sensitivity and action status corresponding to case 1- case 4 in <xref ref-type="table" rid="T5">Table 5</xref> are shown in <xref ref-type="table" rid="T6">Table 6</xref>.</p>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>CRCDP sensitivity and action status corresponding to case 1- case 4.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Case number</th>
<th align="left">
<inline-formula id="inf76">
<mml:math id="m86">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th align="left">Performance of CRCDP</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Case1</td>
<td align="left">1.69</td>
<td align="left">correct action</td>
</tr>
<tr>
<td align="left">Case2</td>
<td align="left">1.26</td>
<td align="left">correct action</td>
</tr>
<tr>
<td align="left">Case3</td>
<td align="left">0.51</td>
<td align="left">correct no-action</td>
</tr>
<tr>
<td align="left">Case4</td>
<td align="left">0.08</td>
<td align="left">correct no-action</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>As can be seen from <xref ref-type="table" rid="T6">Table 6</xref>, CRCDP can correctly identify faults under both charging and discharging conditions. Under the same fault conditions, the sensitivity of CRCDP in the energy storage charging state is higher than that in the discharging state because the voltage on line PN during energy storage charging is lower than that during discharging, resulting in a larger positive sequence reactive current. The positive-sequence voltage amplitude of bus P corresponding to Case 1 and Case 2 is shown in <xref ref-type="fig" rid="F13">Figure 13</xref>.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Positive-sequence voltage magnitude of bus P corresponding to case1 and case2.</p>
</caption>
<graphic xlink:href="fenrg-13-1666514-g013.tif">
<alt-text content-type="machine-generated">Line graph showing voltage changes over time during discharge and charge. The blue line represents discharge, starting high at 210 kV, dropping sharply around 1.5 seconds, and stabilizing near 190 kV. The red line shows charge, beginning at 190 kV, quickly dropping to about 160 kV at 1.5 seconds, and leveling off. Time is on the x-axis (1.45 to 1.8 seconds), and voltage (kV) is on the y-axis (150 to 210 kV).</alt-text>
</graphic>
</fig>
</sec>
</sec>
<sec id="s5">
<title>5 Conclusion</title>
<p>The positive-sequence reactive current differential protection proposed in this paper can effectively reduce the risk of current differential protection no-trip failure caused by excessive phase difference of short-circuit current on both sides of the transmission line when the energy storage power station is charging, and has the following conclusions:<list list-type="simple">
<list-item>
<p>1. The method proposed in this paper only uses the positive-sequence reactive component and can adapt to during both charging and discharging states of the energy storage power station. In addition, since the reactive current support capability provided by PV, direct-drive wind turbines, etc. is similar to that of grid-connected energy storage under the current standards, the protection proposed in this paper can also be applied to the transmission lines of photovoltaic power stations and direct-drive wind farms.</p>
</list-item>
<list-item>
<p>2. The method has sufficient sensitivity under different transition resistances, fault locations and fault types, and performs better when the capacity of the energy storage plant is larger, which can adapt to the development trend of gradual growth of the energy storage power station capacity.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>HZ: Conceptualization, Data curation, Formal Analysis, Methodology, Validation, Visualization, Writing &#x2013; original draft. XW: Investigation, Validation, Writing &#x2013; original draft. CF: Project administration, Supervision, Writing &#x2013; review and editing. YH: Supervision, Writing &#x2013; review and editing. NT: Supervision, Writing &#x2013; review and editing.</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 and/or publication of this article. This work was supported by Shanghai Outstanding Academic Leaders Program (No.22XD1401400) and National Natural Science Foundation of China (No. 52337006).</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>
</sec>
<sec sec-type="ai-statement" id="s10">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<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>Amin</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Suleman</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Waseem</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Iqbal</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Aziz</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Faiz</surname>
<given-names>M. T.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Renewable energy maximization for pelagic islands network of microgrids through battery swapping using deep reinforcement learning</article-title>. <source>IEEE Access</source> <volume>11</volume>, <fpage>86196</fpage>&#x2013;<lpage>86213</lpage>. <pub-id pub-id-type="doi">10.1109/access.2023.3302895</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>L. J.</given-names>
</name>
<name>
<surname>Mei</surname>
<given-names>S. W.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>An integrated control and protection system for photovoltaic microgrids</article-title>. <source>CSEE J. Power Energy Syst.</source> <volume>1</volume> (<issue>1</issue>), <fpage>36</fpage>&#x2013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.17775/cseejpes.2015.00005</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Bo</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Similarity comparison based high-speed pilot protection for transmission line</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>33</volume> (<issue>2</issue>), <fpage>938</fpage>&#x2013;<lpage>948</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2017.2731994</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Farshad</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A pilot protection scheme for transmission lines of half-bridge MMC-HVDC grids using cosine distance criterion</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>36</volume> (<issue>2</issue>), <fpage>1089</fpage>&#x2013;<lpage>1096</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2020.3001878</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>New principle of line differential protection based on combined current and voltage restraint principle</article-title>. <source>Proc. CSEE</source> <volume>1-14</volume>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.232465</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jia</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>C. J.</given-names>
</name>
<name>
<surname>Xuan</surname>
<given-names>Z. W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>Y. Q.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Fault characteristics analysis and line protection design within a large-scale photovoltaic power plant</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>9</volume> (<issue>5</issue>), <fpage>4099</fpage>&#x2013;<lpage>4108</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2017.2648879</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jia</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Bi</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Pilot protection based on amplitude comparison for renewable power teed lines</article-title>. <source>CSEE J. Power Energy Syst.</source> <volume>8</volume> (<issue>6</issue>), <fpage>1519</fpage>&#x2013;<lpage>1529</lpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Lan</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2023</year>). &#x201c;<article-title>A novel principle of current differential protection for renewable energy station</article-title>,&#x201d; in <conf-name>2023 IEEE International Conference on Advanced Power System Automation and Protection (APAP)</conf-name>. <conf-loc>Xuchang, China</conf-loc>, <conf-date>08-12 October 2023</conf-date>, <publisher-name>IEEE</publisher-name>, <fpage>519</fpage>&#x2013;<lpage>523</lpage>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>An enhanced current differential protection for AC transmission lines connecting MMC-HVDC stations</article-title>. <source>IEEE Syst. J.</source> <volume>17</volume> (<issue>1</issue>), <fpage>892</fpage>&#x2013;<lpage>903</lpage>. <pub-id pub-id-type="doi">10.1109/jsyst.2022.3155881</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mishra</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Singh</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Pradhan</surname>
<given-names>A. K.</given-names>
</name>
<name>
<surname>Bajpai</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Protecting distribution systems with inverter-interfaced PV plants using Q-Axis components</article-title>. <source>IEEE Syst. J.</source> <volume>16</volume> (<issue>2</issue>), <fpage>1763</fpage>&#x2013;<lpage>1773</lpage>. <pub-id pub-id-type="doi">10.1109/jsyst.2021.3071627</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mohan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chi</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Asim Amin</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>A market-based real-time algorithm for congestion alleviation incorporating EV demand response in active distribution networks</article-title>. <source>Appl. Energy</source> <volume>356</volume>, <fpage>122426</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2023.122426</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="book">
<collab>National Standards of People&#x27;s Republic of China</collab> (<year>2024</year>). <source>Technical requirements for connecting photovoltaic power station to power system, standardization administration of the People&#x27;s Republic of China GB/T</source>, <fpage>19964</fpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="book">
<collab>National Standards of People&#x27;s Republic of China</collab> (<year>2021</year>). <source>Technical specification for connecting wind farm to power system&#x2014;Part 1: on shore wind power</source>, <volume>19963</volume>. <publisher-loc>Beijing</publisher-loc>: <publisher-name>Standardization Administration of the People&#x27;s Republic of China</publisher-name>.</citation>
</ref>
<ref id="B14">
<citation citation-type="book">
<collab>National Standards of People&#x27;s Republic of China</collab>. (<year>2023</year>). <source>Technical requirements for power conversion system of electrochemical energy storage system, standardization administration of the People&#x27;s Republic of China GB/T 34120-2023</source>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sirisha</surname>
<given-names>A. N. R. L.</given-names>
</name>
<name>
<surname>Pradhan</surname>
<given-names>A. K.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Cosine similarity based directional comparison scheme for subcycle transmission line protection</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>35</volume> (<issue>5</issue>), <fpage>2159</fpage>&#x2013;<lpage>2167</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2019.2962275</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Research progress and prospect of energy storage planning method for new power system</article-title>. <source>Automation Electr. Power Syst.</source> <volume>48</volume> (<issue>9</issue>), <fpage>178</fpage>&#x2013;<lpage>191</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Telukunta</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Pradhan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Agrawal</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Singh</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Srivani</surname>
<given-names>S. G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Protection challenges under bulk penetration of renewable energy resources in power systems: a review</article-title>. <source>CSEE J. Power Energy Syst.</source> <volume>3</volume> (<issue>4</issue>), <fpage>365</fpage>&#x2013;<lpage>379</lpage>. <pub-id pub-id-type="doi">10.17775/cseejpes.2017.00030</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Zang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). &#x201c;<article-title>A current differential protection scheme based on Q-axis component for distribution networks</article-title>,&#x201d; in <conf-name>2021 6th Asia Conference on Power and Electrical Engineering</conf-name>. <conf-loc>Chongqing, China</conf-loc>, <conf-date>08-11 April 2021</conf-date>, <publisher-name>IEEE</publisher-name>, <fpage>24</fpage>&#x2013;<lpage>28</lpage>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Ad-axis based current differential protection scheme for an active distribution network</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>7</volume> (<issue>2</issue>), <fpage>23</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1186/s41601-022-00243-0</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Pearson correlation coefficient of current derivatives based pilot protection scheme for long-distance LCC-HVDC transmission lines</article-title>,&#x201d; in <conf-name>2019 IEEE 8th International Conference on Advanced Power System Automation and Protection (APAP)</conf-name>. <conf-loc>Xi&#x27;an, China</conf-loc>, <conf-date>21-24 October 2019</conf-date>, <publisher-name>IEEE</publisher-name>, <fpage>1367</fpage>&#x2013;<lpage>1371</lpage>.</citation>
</ref>
<ref id="B21">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2024</year>). &#x201c;<article-title>Adaptive current differential protection principle for transmission line connected to energy storage power station based on amplitude and phase compensation</article-title>,&#x201d; in <conf-name>2024 3rd International Conference on Energy and Electrical Power Systems (ICEEPS)</conf-name>. <conf-loc>Guangzhou, China</conf-loc>, <conf-date>14-16 July 2024</conf-date>, <publisher-name>IEEE</publisher-name>, <fpage>49</fpage>&#x2013;<lpage>53</lpage>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Bi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Cosine similarity based line protection for large-scale wind farms</article-title>. <source>IEEE Trans. Industrial Electron.</source> <volume>68</volume> (<issue>7</issue>), <fpage>5990</fpage>&#x2013;<lpage>5999</lpage>. <pub-id pub-id-type="doi">10.1109/tie.2020.2998756</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>