<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Archiving and Interchange DTD v2.3 20070202//EN" "archivearticle.dtd">
<article article-type="methods-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">1342770</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1342770</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Methods</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Optimized operational approach for multi-type reactive power compensation to enhance the grid integration strength of new energy clusters</article-title>
<alt-title alt-title-type="left-running-head">Li 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.2023.1342770">10.3389/fenrg.2023.1342770</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Li</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Qin</surname>
<given-names>Jianru</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2584584/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Hu</surname>
<given-names>Heng</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Qinglei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2586576/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Luo</surname>
<given-names>Lei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>State Grid Shaanxi Electric Power Company Ltd.</institution>, <addr-line>Xian</addr-line>, <addr-line>Shanxi</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Sichuan Energy Internet Research Institute</institution>, <institution>Tsinghua University</institution>, <addr-line>Chengdu</addr-line>, <addr-line>Sichuan</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/901136/overview">Donghui Zhang</ext-link>, Tsinghua University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2586184/overview">Shu Zhang</ext-link>, Sichuan University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2586272/overview">Tianwen Zheng</ext-link>, Semantic Web Company, Austria</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jianru Qin, <email>Qinjr_LK@126.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>03</day>
<month>01</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1342770</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>11</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>12</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Li, Qin, Hu, Zhang and Luo.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Li, Qin, Hu, Zhang and Luo</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>The insufficient system strength in the high-proportion new energy access area has gradually emerged as a crucial factor contributing to the transient overvoltage issue. Therefore, it is imperative to propose a reactive power optimization operation mode that takes into consideration both the power grid strength and system operating voltage of the new energy cluster system. Firstly, the relationship between the evaluation index of power grid strength and the performance of system voltage response is elucidated, while analyzing the influence mechanism of various reactive power compensation devices on the power grid strength of new energy cluster systems. Then, a reactive power operation optimization model is proposed to maximize the strength of the system grid and minimize the voltage deviation. To solve this problem, a hybrid approach combining genetic algorithm and CPLEX solver is employed. Finally, the effectiveness of the proposed method is validated through a typical simulation example.</p>
</abstract>
<kwd-group>
<kwd>power grid strength</kwd>
<kwd>reactive power optimization</kwd>
<kwd>voltage control</kwd>
<kwd>multiple target</kwd>
<kwd>reactive power compensation device</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Smart Grids</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>The issue of voltage stability has emerged as a prominent concern in ensuring the safety and reliability of large-scale power grids in recent years, necessitating an immediate enhancement of dynamic reactive power reserves. On one hand, the issue of &#x201c;Strong DC and weak AC&#x201d; in power grids becomes prominent with the advancement and implementation of DC transmission technology (<xref ref-type="bibr" rid="B17">RAO et al., 2015</xref>; <xref ref-type="bibr" rid="B13">LI et al., 2018</xref>; <xref ref-type="bibr" rid="B16">RAO et al., 2019</xref>). Moreover, due to the dense DC drop point and relatively close electrical distance of the load center at the receiving end, a single AC failure often leads to simultaneous multiple DC commutation failures (<xref ref-type="bibr" rid="B12">LIU et al., 2015</xref>; <xref ref-type="bibr" rid="B26">XU et al., 2019</xref>). The construction of large-scale new energy bases, primarily wind power and photovoltaic, in remote areas has led to prominent issues such as a weak DC power grid and insufficient short-circuit capacity. Additionally, the decline in dynamic reactive power reserve and voltage stability has become more significant (<xref ref-type="bibr" rid="B28">ZHANG et al., 2011</xref>; <xref ref-type="bibr" rid="B5">CHEN et al., 2016</xref>). On the contrary, the rapid increase in load on the asynchronous motor at the load center will result in a decline in bus voltage and further depletion of reactive power within the system, ultimately leading to voltage collapse and large-scale power outage incidents. It is imperative to ensure an adequate supply of dynamic reactive power on a large scale to address this insufficiency in voltage support.</p>
<p>The integration of large-scale new energy into the power grid has prompted experts and scholars, both domestic and foreign, to propose the concept of power grid strength as a measure of the interaction between traditional power grids and connected new energy power electronic equipment. This concept is characterized by two dimensions: frequency support strength and voltage support strength (<xref ref-type="bibr" rid="B3">Australian Energy Market Operator, 2020b</xref>; <xref ref-type="bibr" rid="B2">Australian Energy Market Operator, 2020a</xref>). The voltage support strength refers to the system&#x2019;s ability to maintain voltage stability during significant disturbances, and it is commonly assessed using the short circuit ratio index (<xref ref-type="bibr" rid="B25">XIN et al., 2020</xref>; <xref ref-type="bibr" rid="B30">ZHU et al., 2021</xref>; <xref ref-type="bibr" rid="B23">WU et al., 2022</xref>). In order to assess the strength of new energy grid-connected systems, scholars have proposed the concept of critical short-circuit ratio (<xref ref-type="bibr" rid="B20">WANG et al., 2019</xref>; <xref ref-type="bibr" rid="B4">CAO et al., 2023</xref>) as an index for evaluating the short-circuit ratio. The voltage support capability of the system is assessed by comparing the real-time short-circuit ratio with the critical short-circuit ratio. Due to the existence of various short-circuit ratio index construction methods for different Scale of new energy grid connection and new energy access forms (<xref ref-type="bibr" rid="B9">Lin et al., 2022a</xref>), there is a challenge in ensuring consistency between the calculation of critical short-circuit ratio index and the boundary value used for system strength assessment (<xref ref-type="bibr" rid="B22">Wu et al., 2018</xref>; <xref ref-type="bibr" rid="B10">Lin et al., 2022b</xref>; <xref ref-type="bibr" rid="B27">YANG et al., 2023</xref>). The short-circuit ratio and critical short-circuit ratio calculation methods proposed in literature (<xref ref-type="bibr" rid="B10">Lin et al., 2022b</xref>) can be employed for the evaluation of each new energy station&#x2019;s short-circuit ratio within the new energy multi-feed system. These methods enable identification of weak links in the new energy cluster system and contribute to formulating strategies for strengthening the weak power grid system. Therefore, this paper primarily adopts these literature-based approaches to assess the strength of the system power grid.</p>
<p>The existing literature on voltage control through reactive power compensation can be categorized into two distinct groups. The first type refers to <italic>in situ</italic> reactive power compensation (<xref ref-type="bibr" rid="B8">Kirn and Lee, 2005</xref>; <xref ref-type="bibr" rid="B1">Abdel-Rahman et al., 2006</xref>; <xref ref-type="bibr" rid="B15">Park and Baek, 2012</xref>; <xref ref-type="bibr" rid="B11">LIU et al., 2019</xref>; <xref ref-type="bibr" rid="B6">Fan et al., 2020</xref>; <xref ref-type="bibr" rid="B21">Wen and ZHOU, 2020</xref>). The literature (<xref ref-type="bibr" rid="B6">Fan et al., 2020</xref>) suggests connecting the fast reactive power compensation device to the heavy-load distribution station and coordinating it with the capacitor bank and on-load regulating transformer. The literature (<xref ref-type="bibr" rid="B21">Wen and ZHOU, 2020</xref>) proposed a voltage control technology for distribution stations based on the coordination between SVC and on-load regulator transformers. The literature (<xref ref-type="bibr" rid="B1">Abdel-Rahman et al., 2006</xref>) presents a SVC access strategy incorporating two-stage slope regulation and two-stage voltage regulation control. The literature (<xref ref-type="bibr" rid="B15">Park and Baek, 2012</xref>) proposes a voltage coordination control method based on local reactive power compensation, which is developed through extensive research on shunt capacitors, reactors, and on-load regulating transformers. The literature (<xref ref-type="bibr" rid="B8">Kirn and Lee, 2005</xref>) presents a reactive power compensation approach utilizing artificial neural network. The literature (<xref ref-type="bibr" rid="B11">LIU et al., 2019</xref>) presents a centralized/local adaptive reactive voltage coordination control approach for active distribution networks. The limitation of the above method is obvious, the local compensation does not take into account the tidal flow of the power system, and it has little improvement on the economic and stable operation of the power system. The second approach involves coordinating reactive power compensation devices in the power system to achieve system voltage control (<xref ref-type="bibr" rid="B7">FU et al., 2018</xref>; <xref ref-type="bibr" rid="B14">NI and LI, 2020</xref>; <xref ref-type="bibr" rid="B19">WANG et al., 2020</xref>; <xref ref-type="bibr" rid="B24">XIAO et al., 2020</xref>; <xref ref-type="bibr" rid="B18">TANG et al., 2021</xref>; <xref ref-type="bibr" rid="B29">ZHU and ZHAO, 2021</xref>). The literature (<xref ref-type="bibr" rid="B24">XIAO et al., 2020</xref>) presents a reactive power optimization model that is based on sensitivity analysis. The characteristics of voltage control in load center SVC were analyzed in reference (<xref ref-type="bibr" rid="B19">WANG et al., 2020</xref>). The literature (<xref ref-type="bibr" rid="B29">ZHU and ZHAO, 2021</xref>) proposes a voltage control method that utilizes static synchronous compensation and coordinated optimization of on-load regulator transformers for enhanced performance. The literature (<xref ref-type="bibr" rid="B18">TANG et al., 2021</xref>) presents a regional voltage control approach that can be applied to the penetration of FACTS devices (<xref ref-type="bibr" rid="B14">NI and LI, 2020</xref>). A voltage adaptive coordination control method is proposed in reference (<xref ref-type="bibr" rid="B7">FU et al., 2018</xref>). To summarize, the existing studies solely focus on voltage stability as the single objective during stable system operation, neglecting the requirements for system strength in weak grid scenarios of new energy clusters. Hence, this paper aims to enhance the power grid strength of new energy clusters by proposing a multi-type reactive power compensation optimization method that considers constraints related to short-circuit ratios. This approach enables simultaneous consideration of both voltage stability objectives and system power grid strength requirements.</p>
<p>Taking into account the aforementioned considerations, the relationship between the evaluation index of power grid strength and the voltage response performance of the system is examined, and the influence mechanism of various reactive power compensation devices on the power grid strength in a new energy cluster system is analyzed. Subsequently, a model for optimizing reactive power operation is proposed to maximize system grid strength while minimizing voltage deviation. A hybrid strategy combining genetic algorithm and CPLEX solver is employed to solve this problem. Finally, a typical simulation example is presented to validate the effectiveness of the proposed method.</p>
</sec>
<sec id="s2">
<title>2 Study on the relationship between power grid strength and system voltage response performance</title>
<p>In the power system with high penetration rate of new energy, system strength is a common concern, the system &#x201c;strong&#x201d; and &#x201c;weak&#x201d; is a relative concept, not simply based on the capacity of the new energy equipment connected to the AC grid to evaluate, that is, can not be simply judged as a new energy grid connected to the small capacity of the system, its grid strength is higher than the new energy grid connected to the large capacity of the power system. The short-circuit ratio is currently widely employed to characterize the system strength of new energy clusters (<xref ref-type="bibr" rid="B10">Lin et al., 2022b</xref>), as shown in <xref ref-type="disp-formula" rid="e1">Eq. 1</xref>:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mtext>SCR</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mtext>ac</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x22c5;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mi>Z</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mtext>pu</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>Where, <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mtext>ac</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the short-circuit capacity of the AC system; <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the rated capacity of new energy equipment; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the rated voltage of grid-connected bus; <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:mi>Z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the system impedance; <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mtext>pu</mml:mtext>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the per unit value of system impedance.</p>
<p>The above traditional short-circuit ratio indicators have two limitations for the new energy multi-feed system. Firstly, it is challenging to establish a standardized benchmark capacity due to the varying installed capacities of new energy in each network. Secondly, when considering the impact of junction voltage, it is necessary to take into account not only the directly connected new energy but also other new energy feeds. The proposed short-circuit ratio index <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="normal">S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> based on capacity calculation utilizes the ratio between the AC system&#x2019;s short-circuit capacity and the equivalent grid-connected capacity of new energy. As shown in Eq. <xref ref-type="disp-formula" rid="e2">2</xref>.<disp-formula id="e2">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>ac</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>eq</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>ac</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">U</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>Where, <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the voltage support strength at the <italic>i</italic> node of the system; <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>eq</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the equivalent grid-connected capacity of new energy at the system node <italic>i</italic>; <inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> respectively are new energy installed capacity of system nodes <italic>i</italic>, <italic>j</italic>; <inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">U</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf12">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> respectively are the system <italic>i</italic>, <italic>j</italic> node and voltage; <inline-formula id="inf13">
<mml:math id="m15">
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf14">
<mml:math id="m16">
<mml:mrow>
<mml:msubsup>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>&#x2a;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> respectively are self-impedance and mutual impedance of the system; &#x2a; represents the conjugation operation.</p>
<p>The expression of the relationship between short circuit ratio and junction voltage can be further derived as per the literature derivation (<xref ref-type="bibr" rid="B10">Lin et al., 2022b</xref>).<disp-formula id="e3">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="normal">U</mml:mi>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi mathvariant="normal">N</mml:mi>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>E</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mrow>
<mml:mtext>eq</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>Where, <inline-formula id="inf15">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="normal">U</mml:mi>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the voltage support strength at grid-connected point <italic>i</italic>; <inline-formula id="inf16">
<mml:math id="m19">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>U</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the voltage change at grid-connected point <italic>i</italic>.</p>
<p>The short-circuit ratio index, as indicated by <xref ref-type="disp-formula" rid="e3">formula (3)</xref>, exhibits an inverse relationship with the voltage change at the junction point. A higher short-circuit ratio corresponds to a smaller voltage change at the junction point, indicating enhanced system anti-disturbance capability and increased system strength.</p>
<p>In order to facilitate the analysis of the coupling mechanism between the short-circuit ratio index and the reactive power at the system node, the short-circuit ratio shown in Eq. <xref ref-type="disp-formula" rid="e2">2</xref> is calculated and simplified as follows:<disp-formula id="e4">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>ac</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>j</mml:mi>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>r</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>j</mml:mi>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>r</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>Z</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>Z</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>In Eq. <xref ref-type="disp-formula" rid="e4">4</xref>, <inline-formula id="inf17">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mtext>SCR</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the voltage support strength at the junction point <italic>i</italic>; <italic>P</italic>
<sub>
<italic>i</italic>
</sub> and <italic>P</italic>
<sub>
<italic>j</italic>
</sub> are the active power output values of new energy at connecting points <italic>i</italic> and <italic>j</italic> respectively; <italic>Q</italic>
<sub>
<italic>i</italic>
</sub> and <italic>Q</italic>
<sub>
<italic>j</italic>
</sub> are the reactive power output values of merging points i and j respectively. <italic>R</italic>
<sub>
<italic>ji</italic>
</sub> is the influence factor of voltage interaction between <italic>j</italic> and <italic>i</italic>.</p>
<p>The focus of this paper is primarily on evaluating the operational strength of the short-circuit ratio system, whereby the actual output of new energy at parallel points is utilized for calculating the equivalent short-circuit ratio, rather than relying on the installed capacity of new energy.</p>
</sec>
<sec id="s3">
<title>3 Influence analysis of multi-type reactive power compensation devices on power grid strength of new energy cluster system</title>
<p>In order to enhance the voltage performance of new energy grid-connected systems and strengthen weak grid systems of new energy clusters, common methods for reactive power compensation include the installation of static var compensators (SVC) and static var generators (SVG), as well as utilizing the reactive power emitted by the new energy equipment itself for reactive power compensation. The various methods of reactive power compensation, however, employ distinct mechanisms to enhance the strength of the system grid. The subsequent section will elucidate the mechanism employed by each method during the calculation process for short-circuit ratio.</p>
<p>The calculation process of the short-circuit ratio index is illustrated in <xref ref-type="fig" rid="F1">Figure 1</xref>. Initially, the admittance matrix of the system is established based on its structure (including lines, transformers, <italic>etc.</italic>), incorporating the subtransient reactance of conventional units into the admittance matrix while disregarding the impedance of the inverter interface for new energy equipment, then the system impedance matrix is derived by taking the inverse of the admittance matrix. According to the principle of symmetrical short-circuit calculation, the per unit value of node short-circuit current is obtained by taking the reciprocal of the self-impedance of each system node. Consequently, the short-circuit capacity for each node in the system is determined. Finally, <xref ref-type="disp-formula" rid="e4">formula (4)</xref> and <xref ref-type="disp-formula" rid="e5">5</xref> are utilized to measure the system strength of each new energy node and obtain their respective short-circuit ratios.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>The process of short-circuit ratio calculation.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g001.tif"/>
</fig>
<p>The short-circuit ratio calculation process reveals that the inclusion of various reactive power compensation methods impacts the calculation of each node&#x2019;s short-circuit ratio, subsequently influencing the determination of new energy output limits for each node and ultimately affecting the absorption capacity of each new energy station. Among them, the input of SVC and synchronous compensator primarily impact the short-circuit capacity of each node in the system by modifying the admittance matrix. However, direct access to tuning the synchronous compensator affects the self-admittance matrix of the access node location. SVC is considered as an electrical admittance matrix parallel to the bus, which explains why synchronous compensator has a better compensation effect compared to SVC. The reactive power of Static Var Generator (SVG) and renewable energy sources mainly increases the reactive power of the system node power supply. <xref ref-type="table" rid="T1">Table 1</xref> summarizes the impact mechanism of each reactive power compensation mode on short-circuit ratio calculation.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Influence mechanism of reactive power compensation on new energy consumption.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Modes of reactive power compensation</th>
<th align="center">Influence mechanism</th>
<th align="center">Impact effect</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">SVC</td>
<td align="center">Input capacitance/reactance values are incorporated into the admittance matrix</td>
<td align="center">Self-impedance and mutual impedance increase</td>
</tr>
<tr>
<td align="center">SVG</td>
<td align="center">Participate in short-circuit ratio calculation as a separate power supply</td>
<td align="center">The number of current sources involved in short-circuit ratio calculation increases</td>
</tr>
<tr>
<td align="center">New energy itself reactive power compensation</td>
<td align="center">Does not change the number of current sources involved in the short-circuit ratio calculation</td>
<td align="center">Increases the virtual reactive power of the source output</td>
</tr>
<tr>
<td align="center">Synchronous compensator</td>
<td align="center">Equivalent is a voltage source with inductive internal impedance</td>
<td align="center">Increase short-circuit capacity and reduce self-impedance and mutual impedance</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In order to quantitatively analyze the impact of the aforementioned reactive power compensation methods on the system&#x2019;s power grid strength, a comparative analysis is conducted using a new energy gathering and transmission terminal system as an example to assess the effectiveness of reactive power compensation methods on <italic>K</italic>
<sub>SCR</sub>. The calculation results are presented in <xref ref-type="table" rid="T2">Table 2</xref>, where WT/PV indicates that the reactive power of the new energy station is used to compensate the reactive power. The calculation results demonstrate that the implementation of any of the four reactive power compensation methods can effectively enhance the short-circuit ratio at the grid-connected node of new energy. The increase in <italic>K</italic>
<sub>SCR</sub> for SVC is the smallest among them, with a growth rate ranging from 1% to 8%, and the majority of which are below 5%. The two modes of SVG compensation and new energy reactive power compensation exhibit a similar magnitude of improvement for <italic>K</italic>
<sub>SCR</sub>, ranging from 6% to 11%. The addition of the synchronous compensator has a significantly noticeable impact on <italic>K</italic>
<sub>SCR</sub>, resulting in an evident improvement. The installation of the distributed synchronous compensator at new energy power stations under the gathering station leads to an enhancement level exceeding 40% for <italic>K</italic>
<sub>SCR</sub>, with the highest improvement reaching 83%. Therefore, the utilization of synchronous compensator for reactive power compensation and enhancement of power grid strength has become a prominent focus within the industry.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Comparison of <italic>K</italic>
<sub>SCR</sub> calculation results of some new energy stations under different reactive power compensation methods.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Station number</th>
<th colspan="5" align="center">
<italic>K</italic>
<sub>SCR</sub>
</th>
<th colspan="4" align="center">Upgrade rate/%</th>
</tr>
<tr>
<th align="center">No measures</th>
<th align="center">SVC</th>
<th align="center">SVG</th>
<th align="center">WT/PV</th>
<th align="center">Synchronous compensator</th>
<th align="center">SVC</th>
<th align="center">SVG</th>
<th align="center">WT/PV</th>
<th align="center">Synchronous compensator</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">1</td>
<td align="center">1.7479</td>
<td align="center">1.8838</td>
<td align="center">1.9307</td>
<td align="left">1.9312</td>
<td align="center">3.1985</td>
<td align="center">7.8</td>
<td align="center">10.5</td>
<td align="center">10.5</td>
<td align="center">83.0</td>
</tr>
<tr>
<td align="center">2</td>
<td align="center">1.1196</td>
<td align="center">1.1615</td>
<td align="center">1.1959</td>
<td align="center">1.1962</td>
<td align="center">1.8214</td>
<td align="center">3.7</td>
<td align="center">6.8</td>
<td align="center">6.8</td>
<td align="center">62.7</td>
</tr>
<tr>
<td align="center">3</td>
<td align="center">1.1717</td>
<td align="center">1.1836</td>
<td align="center">1.2566</td>
<td align="center">1.2549</td>
<td align="center">1.9702</td>
<td align="center">1.0</td>
<td align="center">7.2</td>
<td align="center">7.1</td>
<td align="center">68.1</td>
</tr>
<tr>
<td align="center">4</td>
<td align="center">1.3321</td>
<td align="center">1.3989</td>
<td align="center">1.4430</td>
<td align="center">1.4431</td>
<td align="center">2.1384</td>
<td align="center">5.0</td>
<td align="center">8.3</td>
<td align="center">8.3</td>
<td align="center">60.5</td>
</tr>
<tr>
<td align="center">5</td>
<td align="center">1.2786</td>
<td align="center">1.3373</td>
<td align="center">1.3772</td>
<td align="center">1.3774</td>
<td align="center">1.9742</td>
<td align="center">4.6</td>
<td align="center">7.7</td>
<td align="center">7.7</td>
<td align="center">54.4</td>
</tr>
<tr>
<td align="center">6</td>
<td align="center">1.1904</td>
<td align="center">1.2400</td>
<td align="center">1.2735</td>
<td align="center">1.2738</td>
<td align="center">1.7684</td>
<td align="center">4.2</td>
<td align="center">7.0</td>
<td align="center">7.0</td>
<td align="center">48.5</td>
</tr>
<tr>
<td align="center">7</td>
<td align="center">1.3917</td>
<td align="center">1.4247</td>
<td align="center">1.5130</td>
<td align="center">1.5106</td>
<td align="center">2.1633</td>
<td align="center">2.4</td>
<td align="center">8.7</td>
<td align="center">8.5</td>
<td align="center">55.4</td>
</tr>
<tr>
<td align="center">8</td>
<td align="center">1.3156</td>
<td align="center">1.3807</td>
<td align="center">1.4188</td>
<td align="center">1.4192</td>
<td align="center">1.9546</td>
<td align="center">4.9</td>
<td align="center">7.8</td>
<td align="center">7.9</td>
<td align="center">48.6</td>
</tr>
<tr>
<td align="center">9</td>
<td align="center">1.2239</td>
<td align="center">1.2828</td>
<td align="center">1.3173</td>
<td align="center">1.3175</td>
<td align="center">1.7952</td>
<td align="center">4.8</td>
<td align="center">7.6</td>
<td align="center">7.6</td>
<td align="center">46.7</td>
</tr>
<tr>
<td align="center">10</td>
<td align="center">1.3128</td>
<td align="center">1.3795</td>
<td align="center">1.4186</td>
<td align="center">1.4187</td>
<td align="center">1.9699</td>
<td align="center">5.1</td>
<td align="center">8.1</td>
<td align="center">8.1</td>
<td align="center">50.1</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s4">
<title>4 A mathematical model for optimal reactive power operation considering power grid strength and voltage fluctuations</title>
<p>During conventional operations for optimizing reactive power in a power grid, once the placement of each reactive power compensation device is determined, adjustments are made to their outputs (e.g., SVGs and synchronous compensators) in order to comply with system constraints. This ensures economic and secure operation of the grid while achieving minimal voltage deviations at every node. The analysis of the impact of various types of reactive power compensation devices on the power grid strength in the new energy cluster system, as described in Part 2, reveals that adjusting the output of these devices can effectively enhance the short-circuit ratio and overall system strength. Therefore, the paper comprehensively examines the regulatory impact of reactive power compensation devices on system voltage and power grid strength It proposes an optimized operational mode for reactive power optimization that considers both power grid strength requirements and voltage fluctuations.</p>
<sec id="s4-1">
<title>4.1 Objective function</title>
<p>The objective function is defined as the minimum sum of voltage offsets across all nodes in the system, aiming to address the issue of voltage fluctuation. The specific expression is shown as Eq. <xref ref-type="disp-formula" rid="e6">6</xref>: <disp-formula id="e6">
<mml:math id="m23">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x3d;</mml:mo>
</mml:mrow>
<mml:mi>min</mml:mi>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mn>24</mml:mn>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>Where, <italic>M</italic> represents the total number of nodes in the system; <inline-formula id="inf18">
<mml:math id="m24">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> Indicates the voltage deviation of system node <italic>i</italic>. <inline-formula id="inf19">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the Per unit voltage at node <italic>i</italic> of the system.</p>
<p>According to the short-circuit ratio calculation formula and the traditional reactive power-voltage adjustment process, the reactive power required by the short-circuit ratio adjustment process of the new energy node can only be matched with the reactive power required by the traditional reactive power-voltage adjustment process when there is a forward shift in system node voltage (i.e., when the voltage of the new energy node exceeds 1pu). The goal of minimizing voltage deviation and maximizing the power grid strength at new energy junction points necessitates the establishment of a deficiency penalty index &#x394;K when the voltage per unit value exceeds 1 and the short-circuit ratio index is below 2. To achieve this, we establish the following objective function to minimize the deficiency penalty in power grid strength. As shown in Eqs. <xref ref-type="disp-formula" rid="e7">7</xref>, <xref ref-type="disp-formula" rid="e8">8</xref>.<disp-formula id="e7">
<mml:math id="m26">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x3d;</mml:mo>
</mml:mrow>
<mml:mi>min</mml:mi>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mn>24</mml:mn>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math id="m27">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>1000</mml:mn>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>,</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x26;</mml:mo>
</mml:mrow>
<mml:msubsup>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>ESCR</mml:mtext>
</mml:msubsup>
<mml:mrow>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mo>,</mml:mo>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>Where, <italic>N</italic> represents the number of grid-connected new energy sources; <inline-formula id="inf20">
<mml:math id="m28">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the penalty index value of insufficient power grid strength at the new energy connection point <italic>i</italic>; <inline-formula id="inf21">
<mml:math id="m29">
<mml:mrow>
<mml:msubsup>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>ESCR</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the short-circuit ratio index value at the new energy connection point <italic>i</italic>.</p>
<p>In summary, the multi-objective function is shown as <xref ref-type="disp-formula" rid="e9">Eq. 9</xref>:<disp-formula id="e9">
<mml:math id="m30">
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>F</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>Where, <italic>a</italic> and <italic>&#xdf;</italic> are the weight coefficients.</p>
</sec>
<sec id="s4-2">
<title>4.2 Constraint function</title>
<p>The control variables for reactive power optimization include the synchronous compensator and SVG reactive power output, while the state variable is represented by the voltage of each node in the system. The establishment of the following operational constraint equation is necessary to ensure the economical and safe operation of the system by maintaining control variables and state variables within a limited range.<list list-type="simple">
<list-item>
<p>(1) Power flow equilibrium equation. The equation constraint in this paper is based on the active and reactive power balance of each node in the system. As shown in Eq. <xref ref-type="disp-formula" rid="e10">10</xref>.</p>
</list-item>
</list>
<disp-formula id="e10">
<mml:math id="m31">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>G</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>new</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SVG</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SC</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>B</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>Where, <inline-formula id="inf22">
<mml:math id="m32">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>G</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the active power output of conventional power supply of node <italic>i</italic> at time <italic>t</italic>; <inline-formula id="inf23">
<mml:math id="m33">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>new</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the active power output of new energy at node <italic>i</italic> at time <italic>t</italic>; <italic>Q</italic>
<sub>
<italic>i,t</italic>
</sub> represents the reactive power output of node i power supply at time t; <inline-formula id="inf24">
<mml:math id="m34">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the active power output of the node <italic>i</italic> load at time <italic>t</italic>; <inline-formula id="inf25">
<mml:math id="m35">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the load reactive power of node <italic>i</italic> at time <italic>t</italic>; <inline-formula id="inf26">
<mml:math id="m36">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SVG</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> means that the SVG of node <italic>i</italic> emits reactive power at time <italic>t</italic>; <inline-formula id="inf27">
<mml:math id="m37">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SC</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the reactive power emitted by the synchronous compensator of node <italic>i</italic> at time <italic>t</italic>; <italic>U</italic>
<sub>
<italic>j,t</italic>
</sub> represents the voltage amplitude of node j at time t; <italic>G</italic>
<sub>
<italic>ij,t</italic>
</sub>, <italic>B</italic>
<sub>
<italic>ij,t</italic>
</sub> represent the real and imaginary parts of the node admittance matrix at time t, respectively; <italic>&#x3b8;</italic>
<sub>
<italic>ij</italic>
</sub> represents the phase Angle difference between nodes <italic>i</italic> and <italic>j</italic>.<list list-type="simple">
<list-item>
<p>(2) Control variable inequality constraints. The reactive power optimization problem addressed in this paper encompasses the output reactive power control variables of SVG and the output reactive power control variables of the synchronous compensator, which must satisfy the following inequality constraints. As shown in Eqs. <xref ref-type="disp-formula" rid="e11">11</xref>, <xref ref-type="disp-formula" rid="e12">12</xref>.</p>
</list-item>
</list>
<disp-formula id="e11">
<mml:math id="m38">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SVG</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>min</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SVG</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SVG</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m39">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SC</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>min</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SC</mml:mtext>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SC</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>Where, <inline-formula id="inf28">
<mml:math id="m40">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SVG</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>min</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf29">
<mml:math id="m41">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SVG</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represent the lower limit and upper limit of SVG reactive power installed on node <italic>i</italic> respectively; <inline-formula id="inf30">
<mml:math id="m42">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SC</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>min</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf31">
<mml:math id="m43">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>SC</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> Indicates the lower limit and upper limit of reactive power output of synchronous compensator at node <italic>i</italic> respectively.<list list-type="simple">
<list-item>
<p>(3) State variable constraints. The state variable in the optimization problem is the node voltage, which must satisfy both an upper bound and a lower bound. As shown in <xref ref-type="disp-formula" rid="e13">Eq. 13</xref>:</p>
</list-item>
</list>
<disp-formula id="e13">
<mml:math id="m44">
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>min</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>Where, <inline-formula id="inf32">
<mml:math id="m45">
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>min</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf33">
<mml:math id="m46">
<mml:mrow>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the minimum and maximum voltage of node <italic>i</italic> respectively.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Model solving method</title>
<p>The genetic algorithm (GA) is an adaptive global optimization probability search algorithm that simulates the genetic and evolutionary processes of organisms in their natural environment. The variables that need to be solved are encoded into a chromosome for genetic modification, while the objective function corresponding to a set of variable values is transformed into an individual fitness. The initial value is randomly assigned, and after the process of natural selection, individuals that meet the criteria are selected. Genetic operations such as replication, crossover, and mutation are then performed on these individuals to obtain the succeeding generation. The process is iterated until an optimal solution is obtained or the genetic algorithm reaches a predetermined threshold.</p>
<p>When applying genetic algorithms to power system optimization, it is important to recognize that the power flow solution within the system is subject to various constraints. The effectiveness of these solutions is evaluated based on an objective function, whereby those with lower evaluation values are discarded while only those with higher evaluation values have the opportunity to iterate their characteristics into subsequent rounds of solutions. Ultimately, this iterative process aims towards achieving optimization.</p>
<sec id="s5-1">
<title>5.1 Coding</title>
<p>The genetic algorithm is well-suited for mixed integer optimization, making the coding process between the reactive power optimization problem and chromosome highly convenient. The optimization scheme can be represented by each chromosome, thus this paper employs the real number coding method. The synchronous compensator and SVG output reactive power as control variables. Therefore, the chromosome of the control variable can be represented as <xref ref-type="disp-formula" rid="e14">Eq. 14</xref>:<disp-formula id="e14">
<mml:math id="m47">
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SVG</mml:mtext>
</mml:msubsup>
<mml:mrow>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
</mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi mathvariant="normal">S</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>V</mml:mi>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mrow>
<mml:mfenced open="|" close="" separators="|">
<mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SC</mml:mtext>
</mml:msubsup>
<mml:mrow>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
</mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>C</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SC</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>Where, <italic>N</italic>
<sub>
<italic>s</italic>
</sub> and <italic>N</italic>
<sub>
<italic>c</italic>
</sub> represent the total number of SVG and synchronous compensator respectively; <inline-formula id="inf34">
<mml:math id="m48">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SVG</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf35">
<mml:math id="m49">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mtext>SC</mml:mtext>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the reactive power output of the <italic>i</italic> SVG and the synchronous compensator respectively at time <italic>t</italic>.</p>
</sec>
<sec id="s5-2">
<title>5.2 Solution step</title>
<p>
<statement content-type="step" id="Step_1">
<label>Step 1</label>
<p>Input system network parameters, including transmission lines, transformers, shunt capacitors, sub-transient reactance of synchronous generators, load constant impedance. Then calculate the system admittance matrix <italic>Y</italic>
</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_2">
<label>Step 2</label>
<p>Utilize the reactive power output of synchronous synchronous compensator and the SVG as control variables, carry out genetic coding, and generate the initial population;</p>
</statement>
</p>
<p>
<statement content-type="step" id="step_3">
<label>Step 3</label>
<p>Invert the admittance matrix Y to derive the system impedance matrix Z, and extract the self-impedance and mutual impedance of the new energy connection point, SVG installation point, and synchronous synchronous compensator installation point to form a new power supply impedance matrix <bold>
<italic>Z</italic>
</bold>
<sup>new</sup>. The short circuit capacity of the new energy connection point, SVG installation point, and synchronous synchronous compensator installation point is determined based on the power impedance matrix <bold>
<italic>Z</italic>
</bold>
<sup>new</sup>. Then, for each individual of the initial population in step 2, the short-circuit ratio index value under different SVG and synchronous synchronous compensator reactive power output is calculated.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_4">
<label>Step 4</label>
<p>Calculate the voltage value of each node for every individual in the initial population from Step 2, considering different reactive power outputs of SVG and synchronous modulation cameras.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_5">
<label>Step 5</label>
<p>The objective function value of the multi-type reactive power compensation optimization operation method is calculated, and genetic operations are performed on the population, including selection, crossover, and mutation, to obtain a new population.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_6">
<label>Step 6</label>
<p>Determine if the maximum number of iterations has been reached. If so, terminate the calculation and output the result; otherwise, increment the iteration count by one (<inline-formula id="inf36">
<mml:math id="m50">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>) and proceed to step 2.</p>
<p>The algorithm flow is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
</statement>
</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Algrithm flow chart.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g002.tif"/>
</fig>
</sec>
</sec>
<sec id="s6">
<title>6 Case study</title>
<sec id="s6-1">
<title>6.1 Introduction to the case study system</title>
<p>The effectiveness of the proposed method is verified by using the improved IEEE9-node system as an example in this paper, as depicted in <xref ref-type="fig" rid="F3">Figure 3</xref>. On the basis of the original IEEE9-node system, thermal power units with capacities of 100MW, 150&#xa0;MW and 150&#xa0;MW are respectively placed on the original generator nodes, and 250&#xa0;MW wind power, 250&#xa0;MW wind power and 250&#xa0;MW PV are respectively placed at nodes 4, 6 and 8. At the same time, 50Mvar synchronous compensator &#x2b;250Mvar SVG, 200Mvar SVG and 200Mvar SVG are placed at nodes 4, 6 and 8 respectively. This configuration establishes a weak power grid system predominantly driven by new energy sources.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Improved IEEE9 node example system.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g003.tif"/>
</fig>
</sec>
<sec id="s6-2">
<title>6.2 The scenario setting of the comparative example</title>
<p>The present study aims to investigate the impact of various reactive power compensation devices on the strength and voltage of the system grid. To achieve this, two different schemes have been devised in this paper:</p>
<p>Scheme 1: The conventional approach to reactive power optimization solely focuses on the objective of minimizing voltage deviation.</p>
<p>Scheme 2: Reactive power optimization (reactive power optimization model proposed in this paper) takes into account both the maximization of new energy grid-connected strength and the minimization of voltage deviation.</p>
</sec>
<sec id="s6-3">
<title>6.3 Comparative analysis of daily operating voltage</title>
<p>As depicted in <xref ref-type="fig" rid="F4">Figure 4</xref>, node 4 represents a new energy node that encounters the issue of overvoltage. Considering the limitations imposed by the system power grid&#x2019;s strength, the voltage deviation between node 4 and node 1 is comparatively smaller when compared to conventional reactive power optimization strategies.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The comparison of daily voltage for node 4.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g004.tif"/>
</fig>
<p>Taking the operating voltage of nodes 7 as an example, it can be observed from <xref ref-type="fig" rid="F5">Figure 5</xref> that in the absence of a power supply at node 7 within the system, no overvoltage issue arises. However, this observation fails to satisfy the essential conditions of the proposed reactive power optimization model in this paper, which aims to maximize power grid strength. Consequently, compared to the reactive power optimization strategy presented in this study, traditional approaches exhibit smaller voltage deviations and fluctuations.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>The comparison of daily voltage for node 7.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g005.tif"/>
</fig>
</sec>
<sec id="s6-4">
<title>6.4 Comparative analysis of voltage offset</title>
<p>Taking node 8&#x2019;s voltage deviation as an example, <xref ref-type="fig" rid="F6">Figure 6</xref> illustrates that the traditional reactive power optimization strategy results in a maximum voltage offset of 0.072 and a minimum voltage offset of 0.0026. However, by incorporating the goal of maximizing power grid strength, the maximum voltage offset is reduced to 0.0639 while achieving a minimum voltage offset of zero. Hence, the proposed reactive power optimization method outperforms the traditional approach.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>The voltage offset of node 8.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g006.tif"/>
</fig>
</sec>
<sec id="s6-5">
<title>6.5 Comparative analysis of short circuit ratio in daily operation</title>
<p>The analysis of <xref ref-type="fig" rid="F7">Figure 7</xref> reveals that, in consideration of the strength constraints imposed by the system power grid, node 4 exhibits a higher daily operating short-circuit ratio compared to the critical short-circuit ratio 2. However, under the traditional reactive power optimization strategy, there are certain time periods where the operating short-circuit ratio of node 4 falls below the critical value of 2. This discrepancy fails to fully satisfy the strength requirements of the system power grid and hampers large-scale integration of new energy sources.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>The comparison of short-circuit ratio in daily operations for node 4.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g007.tif"/>
</fig>
</sec>
<sec id="s6-6">
<title>6.6 Comparative analysis of daily running reactive power output</title>
<p>The analysis in <xref ref-type="fig" rid="F8">Figure 8</xref> reveals that nodes 4 and 8 necessitate lesser daily reactive power compensation, considering the strength constraints of the system power grid, compared to the conventional reactive power optimization strategy. This approach effectively reduces the demand for reactive power compensation capability within the system and enhances its dynamic reactive power reserve against significant disturbances.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Comparison of reactive power output of nodes 4, 6 and 8.</p>
</caption>
<graphic xlink:href="fenrg-11-1342770-g008.tif"/>
</fig>
<p>In summary, the proposed optimal operation model of reactive power compensation, considering the constraint of new energy limit grid-connected strength, can simultaneously achieve the objectives of minimizing voltage deviation and maximizing power grid strength under overvoltage conditions.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s7">
<title>7 Conclusion</title>
<p>The paper analyzes the influence mechanism of multi-type reactive power compensation devices on the power grid strength of new energy cluster systems, and proposes a reactive power operation optimization model that aims to maximize the power grid strength while minimizing voltage deviation. The specific conclusions are as follows.<list list-type="simple">
<list-item>
<p>(1) The influence of various reactive power compensation methods on short-circuit ratio calculation is analyzed based on the calculation process, and the impact of multi-type reactive power compensation methods on enhancing the power grid strength is quantified.</p>
</list-item>
<list-item>
<p>(2) The calculation results demonstrate that both the reactive power compensation method of the new energy equipment itself and the SVG reactive power compensation method exhibit similar effects in enhancing system power grid strength, surpassing the SVC compensation method. Notably, the installation of synchronous condenser yields the most pronounced effect on improving new energy absorption capacity.</p>
</list-item>
<list-item>
<p>(3) After considering the objective of maximizing the system grid strength, the reactive power compensation required by the new energy interconnection points is reduced compared to that of traditional reactive power optimization strategies. This effectively enhances the dynamic reserve capacity of reactive power in dealing with significant disturbances.</p>
</list-item>
</list>
</p>
<p>The next research direction of this paper is to establish a comprehensive reactive power resource optimization planning model, encompassing various types of reactive power compensation measures, in order to enhance the new energy consumption capacity of weak power grids while considering economic factors.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s8">
<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 author.</p>
</sec>
<sec id="s9">
<title>Author contributions</title>
<p>LiL: Validation, Writing&#x2013;review and editing. JQ: Conceptualization, Methodology, Writing&#x2013;original draft, Writing&#x2013;review and editing. HH: Data curation, Writing&#x2013;review and editing. QZ: Investigation, Writing&#x2013;review and editing. LeL: Investigation, Writing&#x2013;review and editing.</p>
</sec>
<sec id="s10">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study received funding from State Grid Shaanxi Electric Power Co., Ltd. Science and Technology Project, SGSN0000DKJS2313584.</p>
</sec>
<sec sec-type="COI-statement" id="s11">
<title>Conflict of interest</title>
<p>Authors LiL, QZ, and LeL were employed by State Grid Shaanxi Electric Power Company Ltd.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
<p>The handling editor DZ declared a past collaboration with the author JQ.</p>
<p>The authors declare that this study received funding from State Grid Shaanxi Electric Power Co., Ltd. The funder had the following involvement in the study: collection of data.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<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>Abdel-Rahman</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Youssef</surname>
<given-names>F. M. H.</given-names>
</name>
<name>
<surname>Saber</surname>
<given-names>A. A.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>New static var compensator control strategy and coordination with under-load tap changer</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>21</volume> (<issue>3</issue>), <fpage>1630</fpage>&#x2013;<lpage>1635</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2005.858814</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="book">
<collab>Australian Energy Market Operator</collab> (<year>2020a</year>). <source>Renewable integration study: stage 1 report</source>.</citation>
</ref>
<ref id="B3">
<citation citation-type="book">
<collab>Australian Energy Market Operator</collab> (<year>2020b</year>). <source>System strength in the NEM explained</source>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhuang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Orienta-tion and application prospect of virtual synchronous generator in new power system</article-title>. <source>Automation Electr. Power Syst.</source> <volume>47</volume> (<issue>4</issue>), <fpage>190</fpage>&#x2013;<lpage>207</lpage>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Min</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Study on the mechanism of transient voltage stability of new energy source with power electronic interface</article-title>. <source>Power Syst. Prot. Control</source> <volume>44</volume> (<issue>9</issue>), <fpage>15</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.7667/PSPC150481</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ai</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Static voltage stability margin evolution based on gradient boosting regression tree</article-title>. <source>Electr. Meas. &#x26; Instrum.</source> <volume>57</volume> (<issue>20</issue>), <fpage>39</fpage>&#x2013;<lpage>45</lpage>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Adaptive and coordinated volt/var optimization for unbalanced active distribution networks of multiple voltage levels</article-title>. <source>Power Syst. Technol.</source> <volume>42</volume> (<issue>7</issue>), <fpage>2136</fpage>&#x2013;<lpage>2147</lpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kirn</surname>
<given-names>G. W.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>K. Y.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Coordination control of ULTC transformer and STATCOM based on an artificial neural network</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>20</volume> (<issue>2</issue>), <fpage>580</fpage>&#x2013;<lpage>586</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2005.846205</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>Y. U.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2022a</year>). <article-title>Overview of strength quantification indexes of power system with power electronic equipment</article-title>. <source>Proc. CSEE</source> <volume>42</volume> (<issue>2</issue>), <fpage>499</fpage>&#x2013;<lpage>515</lpage>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>Y. U.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2022b</year>). <article-title>Short circuit ratio index analysis and critical short circuit ratio calculation of renewable energy grid-connected system</article-title>. <source>Proc. CSEE</source> <volume>42</volume> (<issue>03</issue>), <fpage>919</fpage>&#x2013;<lpage>929</lpage>.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ye</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
</person-group>, (<year>2019</year>). <article-title>Adaptive centralized/local Volt-VAR control in active distribution network</article-title>. <source>Power Syst. Clean Energy</source> <volume>35</volume> (<issue>11</issue>), <fpage>20</fpage>&#x2013;<lpage>29</lpage>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Research on reactive compensation strategies for improving stability level of sending-end of 750 kV grid in Northwest China</article-title>. <source>Proc. CSEE</source> <volume>35</volume> (<issue>5</issue>), <fpage>1015</fpage>&#x2013;<lpage>1022</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>C. A. O.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Emergency control of synchronous condenser to suppress DC continuous commuta tion failure</article-title>. <source>Automation Electr. Power Syst.</source> <volume>42</volume> (<issue>22</issue>), <fpage>91</fpage>&#x2013;<lpage>97</lpage>. <pub-id pub-id-type="doi">10.7500/AEPS20170726005</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ni</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Overview of unified power quality conditioner research development</article-title>. <source>Power Syst. Prot. Control</source> <volume>48</volume> (<issue>20</issue>), <fpage>177</fpage>&#x2013;<lpage>187</lpage>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Park</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Baek</surname>
<given-names>Y. S.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Coordination control of voltage between STATCOM and reactive power compensation devices in steady-state</article-title>. <source>J. Electr. Eng. Technol.</source> <volume>7</volume> (<issue>5</issue>), <fpage>689</fpage>&#x2013;<lpage>697</lpage>. <pub-id pub-id-type="doi">10.5370/jeet.2012.7.5.689</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rao</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Leng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tu</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Analysis of the global HVDC power transmission development and the sug gestion of the HVDC export</article-title>. <source>South. Power Syst. Tech. nology</source> <volume>13</volume> (<issue>10</issue>), <fpage>1</fpage>&#x2013;<lpage>7</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rao</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Practice and analysis of UHVDC transmission</article-title>. <source>High. Volt. Eng.</source> <volume>41</volume> (<issue>8</issue>), <fpage>2481</fpage>&#x2013;<lpage>2488</lpage>. <pub-id pub-id-type="doi">10.13336/j.1003-6520.hve.2015.08.002</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Ming</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Ouyang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Reactive power and voltage coordinationcontrol strategy of dfig-based wind farm under wind speed variations</article-title>. <source>Guang Dong Electr. Power</source> <volume>34</volume> (<issue>5</issue>), <fpage>19</fpage>&#x2013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1007-290X.2021.005.003</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Coordinated control of steady-state voltage for UHVDC/AC power systems with large-scale synchronous condenser integration</article-title>. <source>Power Syst. Prot. Control</source> <volume>48</volume> (<issue>24</issue>), <fpage>120</fpage>&#x2013;<lpage>127</lpage>. <pub-id pub-id-type="doi">10.19783/j.cnki.pspc.200060</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A quantita-tive evaluation approach for critical operating SCR of MMC system based on small-signal stability and operating con-straints</article-title>. <source>Proc. CSEE</source> <volume>39</volume> (<issue>10</issue>), <fpage>2853</fpage>&#x2013;<lpage>2863</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.180793</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wen</surname>
<given-names>Q. U.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Implementation of a new reactive power compensation mode for 10 kV power users</article-title>. <source>Power Capacitor React. Power Compens.</source> <volume>41</volume> (<issue>4</issue>), <fpage>130</fpage>&#x2013;<lpage>132</lpage>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Di</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Javadi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Malyscheff</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Hong</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>J. N.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Assessing impact of renewable energy integration on system strength using site-dependent short circuit ratio</article-title>. <source>IEEE Trans. Sustain. Energy</source> <volume>9</volume> (<issue>3</issue>), <fpage>1072</fpage>&#x2013;<lpage>1080</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2017.2764871</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Stability analysis of dense wind power area based on short circuit ratio index</article-title>. <source>Electr. Power Autom. Equip.</source> <volume>42</volume> (<issue>08</issue>), <fpage>72</fpage>&#x2013;<lpage>78</lpage>.</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Voltage stability of weak sending-end system with large-scale grid-connected photovoltaic power plants</article-title>. <source>Electr. Power</source> <volume>53</volume> (<issue>11</issue>), <fpage>31</fpage>&#x2013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.11930/j.issn.1004-9649.202007045</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xin</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Gan</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Ju</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Generalized short circuit ratio of power systems with multiple power elec-tronic devices: analysis for various renewable power generations</article-title>. <source>Proc. CSEE</source> <volume>40</volume> (<issue>17</issue>), <fpage>5516</fpage>&#x2013;<lpage>5526</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.200827</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Operation characteristic of synchronous condenser and SVC for dynamic re active power compensation of DC converter station</article-title>. <source>J. North China Electr. Power Univ. Nat. Sci. Edition)</source> <volume>46</volume> (<issue>1</issue>), <fpage>47</fpage>&#x2013;<lpage>53</lpage>. <pub-id pub-id-type="doi">10.3969/j.ISSN.1007-2691.2019.01.07</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Tran-sient Angle instability mode and mechanism of distributed synchronous condensers in renewable energy station</article-title>. <source>Automation Electr. Power Syst.</source> <volume>47</volume> (<issue>3</issue>), <fpage>12</fpage>&#x2013;<lpage>18</lpage>.</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Transient voltage stability study about a regional grid integrated with wind power</article-title>. <source>Electr. Power</source> <volume>44</volume> (<issue>9</issue>), <fpage>17</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1004-9649.2011.09.004</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Voltage control strategy of a UHV AC power grid considering the influence of line power flow fluctuation on bus voltage</article-title>. <source>Power Syst. Prot. Control</source> <volume>49</volume> (<issue>6</issue>), <fpage>76</fpage>&#x2013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.19783/j.cnki.pspc.200567</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Qu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Chun</surname>
<given-names>L. I. U.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Improved cal-culation method of equivalent short-circuit ratio for power generation cluster of renewable energy</article-title>. <source>Automation Electr. Power Syst.</source> <volume>45</volume> (<issue>22</issue>), <fpage>74</fpage>&#x2013;<lpage>82</lpage>.</citation>
</ref>
</ref-list>
</back>
</article>