<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">EE</journal-id>
<journal-id journal-id-type="nlm-ta">EE</journal-id>
<journal-id journal-id-type="publisher-id">EE</journal-id>
<journal-title-group>
<journal-title>Energy Engineering</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-0118</issn>
<issn pub-type="ppub">0199-8595</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">69610</article-id>
<article-id pub-id-type="doi">10.32604/ee.2025.069610</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Optimum Operation of Low-Voltage AC/DC Distribution Areas with Embedded DC Interconnections under Three-Phase Unbalanced Compensation Conditions</article-title>
<alt-title alt-title-type="left-running-head">Optimum Operation of Low-Voltage AC/DC Distribution Areas with Embedded DC Interconnections under Three-Phase Unbalanced Compensation Conditions</alt-title>
<alt-title alt-title-type="right-running-head">Optimum Operation of Low-Voltage AC/DC Distribution Areas with Embedded DC Interconnections under Three-Phase Unbalanced Compensation Conditions</alt-title>
</title-group>
<contrib-group>
<contrib id="author-1" contrib-type="author">
<name name-style="western"><surname>Tan</surname><given-names>Zhukui</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>Zhou</surname><given-names>Dacheng</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-3" contrib-type="author">
<name name-style="western"><surname>Deng</surname><given-names>Song</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-4" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Li</surname><given-names>Jikai</given-names></name><xref ref-type="aff" rid="aff-2">2</xref><email>lijikai@xj.cee-group.cn</email></contrib>
<contrib id="author-5" contrib-type="author">
<name name-style="western"><surname>Wu</surname><given-names>Zhuang</given-names></name><xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-6" contrib-type="author">
<name name-style="western"><surname>Feng</surname><given-names>Qihui</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-7" contrib-type="author">
<name name-style="western"><surname>Zhang</surname><given-names>Xuan</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<aff id="aff-1"><label>1</label><institution>Electric Power Research Institute of Guizhou Power Grid Co., Ltd.</institution>, <addr-line>Guiyang, 550002</addr-line>, <country>China</country></aff>
<aff id="aff-2"><label>2</label><institution>XJ Electric Co., Ltd.</institution>, <addr-line>Xuchang, 461000</addr-line>, <country>China</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Jikai Li. Email: <email>lijikai@xj.cee-group.cn</email></corresp>
</author-notes>
<pub-date date-type="collection" publication-format="electronic">
<year>2026</year>
</pub-date>
<pub-date date-type="pub" publication-format="electronic">
<day>27</day><month>2</month><year>2026</year>
</pub-date>
<volume>123</volume>
<issue>3</issue>
<elocation-id>4</elocation-id>
<history>
<date date-type="received">
<day>26</day>
<month>6</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>06</day>
<month>8</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2026 The Authors. Published by Tech Science Press.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>The Authors</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_EE_69610.pdf"></self-uri>
<abstract>
<p>This paper presents an optimal operation method for embedded DC interconnections based on low-voltage AC/DC distribution areas (EDC-LVDA) under three-phase unbalanced compensation conditions. It can optimally determine the transmission power of the DC and AC paths to simultaneously improve voltage quality and reduce losses. First, considering the embedded interconnected, unbalanced power structure of the distribution area, a power flow calculation method for EDC-LVDA that accounts for three-phase unbalanced compensation is introduced. This method accurately describes the power flow distribution characteristics under both AC and DC power allocation scenarios. Second, an optimization scheduling model for EDC-LVDA under three-phase unbalanced conditions is developed, incorporating network losses, voltage quality, DC link losses, and unbalance levels. The proposed model employs an improved particle swarm optimization (IPSO) two-layer algorithm to autonomously select different power allocation coefficients for the DC link and AC section under various operating conditions. This enables embedded economic optimization scheduling while maintaining compensation for unbalanced conditions. Finally, a case study based on the IEEE 13-node system for EDC-LVDA is conducted and tested. The results show that the proposed optimal operation method achieves a 100% voltage compliance rate and reduces network losses by 13.8%, while ensuring three-phase power balance compensation. This provides a practical solution for the modernization and upgrading of low-voltage power grids.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>Power loss optimization</kwd>
<kwd>low-voltage AC/DC distribution areas</kwd>
<kwd>embedded DC interconnections</kwd>
</kwd-group>
<funding-group>
<award-group id="awg1">
<funding-source>China Southern Power Grid Corporation</funding-source>
<award-id>GZKJXM20220041</award-id>
</award-group>
<award-group id="awg2">
<funding-source>National Key Research and Development Plan</funding-source>
<award-id>2022YFE0205300</award-id>
</award-group>
</funding-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label>
<title>Introduction</title>
<p>With the promotion of low-carbon power systems, a significant number of distributed photovoltaic systems are being integrated into rural distribution networks. The output of photovoltaic systems is characterized by strong randomness and intermittency, and a high level of integration on the low-voltage side worsens the imbalance of the three-phase distribution network, leading to issues such as unbalanced feeder loads, frequent voltage fluctuations that exceed limits, and reverse power flows [<xref ref-type="bibr" rid="ref-1">1</xref>&#x2013;<xref ref-type="bibr" rid="ref-4">4</xref>]. Embedded direct current interconnections (EDC) use flexible power electronics technology to create direct current links within traditional AC distribution areas, effectively addressing challenges related to dynamic capacity expansion and power quality in distribution networks, making it a promising direction for future flexible distribution networks [<xref ref-type="bibr" rid="ref-5">5</xref>,<xref ref-type="bibr" rid="ref-6">6</xref>].</p>
<p>The EDC-based low-voltage distribution areas (EDC-LVDA) consist of two types: interconnections between multiple areas and interconnections within a single area. These use AC-DC voltage source converters (VSC) to enhance power supply capabilities and improve power quality in distribution areas [<xref ref-type="bibr" rid="ref-7">7</xref>]. The interconnection between areas focuses on solving load balancing issues among different area transformers, while the interconnection within an area aims to manage low voltage at the end of the area, enable dynamic capacity expansion, and address three-phase imbalance issues [<xref ref-type="bibr" rid="ref-8">8</xref>,<xref ref-type="bibr" rid="ref-9">9</xref>]. However, commonly used three-phase three-leg (3P3L) VSCs have limited ability to handle asymmetric and non-linear loads, which prevents effective compensation of three-phase unbalanced power. Conversely, the three-phase four-leg (3P4L) VSC-based EDC-LVDC includes a separately controllable fourth leg that allows for independent control of each phase&#x2019;s power, providing an effective solution for three-phase unbalance compensation in low-voltage distribution areas [<xref ref-type="bibr" rid="ref-10">10</xref>,<xref ref-type="bibr" rid="ref-11">11</xref>].</p>
<p>The EDC-LVDC provides a DC path for power transmission, and the transmission power of the DC path will have a significant impact on the system&#x2019;s voltage distribution, loss characteristics, etc. Reasonable optimal regulation of the power in this DC path is crucial to ensuring the high power quality and highly reliable operation of EDC-LVDA. References [<xref ref-type="bibr" rid="ref-12">12</xref>,<xref ref-type="bibr" rid="ref-13">13</xref>] present a multi-time-scale optimization scheduling scheme for a flexible interconnected distribution network that integrates SOP and VSC technologies. Reference [<xref ref-type="bibr" rid="ref-14">14</xref>] proposes an optimization scheduling model for AC-DC low-voltage distribution networks based on multi-mode flexible interconnection, using convex optimization theory for model convexification and linearization. Reference [<xref ref-type="bibr" rid="ref-15">15</xref>] addresses the issue of source-load uncertainty in flexible interconnected distribution networks by proposing an optimization scheduling method based on improved model predictive control. Reference [<xref ref-type="bibr" rid="ref-16">16</xref>] introduces a comprehensive regulation strategy for three-phase imbalance and heavy-light load operation issues in flexible interconnected distribution areas, based on a four-bridge-arm intelligent soft switch, to meet the demands for comprehensive governance of these problems. However, previous research primarily focuses on load balancing power scheduling for interconnection between multiple distribution areas, with little attention given to interconnection scheduling within a single distribution area. Particularly under three-phase imbalance conditions, the imbalance scheduling commands for VSCs within the distribution area can significantly affect the distribution characteristics of the power flow, necessitating the establishment of a power flow calculation model and an optimization scheduling method that takes into account the management of three-phase imbalance conditions within the distribution area.</p>
<p>To address these issues, this paper proposes an optimum operation of EDC-LVDA under three-phase unbalanced compensation conditions, which can optimally determine the transmission power of the DC and AC path to simultaneously achieve optimization of voltage quality and losses. Firstly, a power flow calculation method for EDC-LVDA under three-phase imbalance compensation conditions is proposed, which considers the power flow distribution characteristics. Secondly, under three-phase unbalanced conditions, an optimization scheduling model for flexible low-voltage distribution areas is constructed, incorporating network losses, voltage quality, DC circuit losses, and imbalance levels. The Improved particle swarm optimization (IPSO) algorithm is employed to adaptively select the power distribution coefficients for DC and AC paths under different operating conditions, achieving embedded economic optimization scheduling while ensuring compensation for imbalance conditions. The potential applications of this method focus on urban and rural low-voltage distribution networks, renovation of aging power grids, and AC/DC hybrid microgrids. It can solve three-phase unbalance issues, improve voltage stability, reduce network losses, adapt to the integration of distributed energy resources, and is suitable for scenarios with large load fluctuations or in need of upgrading.</p>
<p>The paper is organized as follows: the second part of this paper introduces the operation of 3P4L-based EDC-LVDA under unbalanced compensation conditions. The third part introduces the power flow calculation method for EDC-LVDA under unbalanced compensation conditions. The fourth part introduces a scheduling optimization strategy for EDC-LVDA. The fifth part presents the case study analysis, and the sixth part is the conclusion of this paper.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Operation of EDC-LVDA under Unbalanced Compensation Conditions</title>
<p><xref ref-type="fig" rid="fig-1">Fig. 1</xref> shows the basic structure of EDC-LVDA. At both the beginning and the end of the distribution area, VSC converters are used to construct embedded DC links, enabling dynamic capacity expansion of the original AC distribution area and compensating for the three-phase imbalance of the distribution transformers in the area. The connections of different PVs and loads are unordered, and not all of them are three-phase symmetric connections, which is also the reason for the three-phase imbalance in the distribution network. To achieve the three-phase compensation of the AC area, the VSC uses the 3P4L-based VSC, its topology is plotted in <xref ref-type="fig" rid="fig-1">Fig. 1</xref>. The typical operation of the EDC-LVDA can be described as follows:</p>
<p><list list-type="order">
<list-item>
<p>The VSC1 at the beginning is responsible for controlling the DC bus voltage of the embedded DC link, allowing for flexible routing of power between the source and the load.</p></list-item>
<list-item>
<p>The VSC2 at the tail end outputs a precise phase-shifted signal based on the power signal received from the head end, which uses a split-phase active and reactive power control to achieve compensation for three-phase power imbalances. The details operation of active and reactive power control can be found in reference to [<xref ref-type="bibr" rid="ref-17">17</xref>].</p></list-item>
</list></p>
<fig id="fig-1">
<label>Figure 1</label>
<caption>
<title>The basic structure and operation of EDC-LVDA</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-1.tif"/>
</fig>
<p>The embedded DC link not only compensates for unbalanced power but also enables dynamic capacity expansion of the AC area, providing a new energy routing channel for the power loads in the area. The power processed by the DC-links will affect the voltage quality and power loss of the whole distribution area, thus the power analysis is presented.</p>
<p>The goal of the EDC-LVDA is to achieve three-phase power balance in transformers, thus can be obtained as
<disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>3</mml:mn></mml:mfrac><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>P</italic><sub><italic>Ti</italic></sub> (<italic>i</italic> &#x003D; <italic>a</italic>, <italic>b</italic>, <italic>c</italic>) is each phase power at the input transformer, <italic>P</italic><sub><italic>Li</italic></sub> (<italic>i</italic> &#x003D; <italic>a</italic>, <italic>b</italic>, <italic>c</italic>) is each phase load, which is unbalanced due to the strong randomness and intermittency of the load and PVs. <italic>P</italic><sub><italic>loss</italic></sub> is the power loss of the whole EDC-LVDA system.</p>
<p>Defining the power distribution ratio <italic>k</italic> (0 &#x003C; <italic>k</italic> &#x003C; 1) to stand for the ratio of the power processed by the EDC path that accounts for the total power, i.e.,
<disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula>where <italic>P</italic><sub>1<italic>i</italic></sub> (<italic>i</italic> &#x003D; <italic>a</italic>, <italic>b</italic>, <italic>c</italic>) is each phase power for VSC1. Thus, the power of the EDC <italic>P</italic><sub><italic>dc</italic></sub> can be obtained as
<disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>p</italic><sub><italic>total</italic></sub> is the total power of the loads in EDC-LVDA.</p>
<p>VSC2 is using the split-phase active and reactive power control to compensate for the load unbalance for each phase, thus the following stages can be obtained as
<disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>a</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:mi>k</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>T</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>P</italic><sub>2<italic>i</italic></sub> (<italic>i</italic> &#x003D; <italic>a</italic>, <italic>b</italic>, <italic>c</italic>) is each phase power for VSC2, which is unequal under the unbalanced compensation conditions.</p>
<p>Based on the above analysis, the power distribution ratio <italic>k</italic> directly determines the power processed by the DC path. If <italic>k</italic> &#x003D; 0, which means the power is fully transferred by the AC path, and <italic>k</italic> &#x003D; 1 means that the power is fully transferred by the EDC part. Meanwhile, due to the different characteristics of the AC and DC parts, a different power distribution ratio <italic>k</italic> will affect the voltage quality of each port and the power losses of the whole area. Therefore, how to allocate the power between the AC and DC routing paths, while ensuring three-phase balance in the embedded distribution transformer and optimizing the voltage quality of the AC/DC area and the overall network loss, is a critical issue that needs to be addressed urgently.</p>
</sec>
<sec id="s3">
<label>3</label>
<title>Power Flow Calculation Method for EDC-LVDA under Unbalanced Compensation Conditions</title>
<p>The foundation for optimizing the power distribution ratio <italic>k</italic>, which affects the losses and power quality in EDC-LVDA, is the accurate calculation of the power flow under imbalance conditions. This paper proposes a power flow calculation method for EDC-LVDA under three-phase imbalance compensation conditions, which considers the power flow distribution characteristics under AC-DC power distribution. The details are as follows.</p>
<p><xref ref-type="fig" rid="fig-2">Fig. 2</xref> shows the detailed calculation diagram of the power flow, which consists of five steps as follows.</p>
<fig id="fig-2">
<label>Figure 2</label>
<caption>
<title>The calculation diagram for EDC-LVDA under unbalanced compensation conditions</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-2.tif"/>
</fig>
<p>Step 1: Initialize and calculate the AC power flow to get the total power and load conditions.</p>
<p>Step 2: Calculate <italic>P</italic><sub><italic>dc</italic></sub> with <xref ref-type="disp-formula" rid="eqn-3">Eq. (3)</xref> and recalculate the AC power flow with the updated AC power [<xref ref-type="bibr" rid="ref-18">18</xref>].</p>
<p>Step 3: Calculate each phase power of the VSC1 and VSC2 with <xref ref-type="disp-formula" rid="eqn-2">Eqs. (2)</xref> and <xref ref-type="disp-formula" rid="eqn-4">(4)</xref>.</p>
<p>Step 4: Iterative calculation of AC power flow.</p>
<p>Step 5: If the calculations do not converge, re-run the iterative calculations.</p>
<p>Based on the above calculation diagram, the iterative AC/DC compensation process in the EDC-LVDA can be plotted in <xref ref-type="fig" rid="fig-3">Fig. 3</xref>, and the detailed descriptions can be obtained as follows.</p>
<fig id="fig-3">
<label>Figure 3</label>
<caption>
<title>The iterative compensation power process for AC-DC in an EDC-LVDA system</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-3.tif"/>
</fig>
<p><list list-type="simple">
<list-item><label>(1)</label><p>Set the initial values for each substation area and input the raw data.</p></list-item>
<list-item><label>(2)</label><p>First iteration: Determine the power transfer coefficient k and the parameters at the EDC. Set the control mode of the upstream VSC to dc-bus, while the downstream VSC is set to split-phase active and reactive power control. Use the following equation to balance the loads within the distribution area.</p></list-item>
</list>
<disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>c</mml:mi><mml:mi mathvariant="normal">&#x005F;</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>&#x03BB;</mml:mi><mml:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi><mml:mi mathvariant="normal">&#x005F;</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi><mml:mi mathvariant="normal">&#x005F;</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>D</mml:mi><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>n</mml:mi><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:mstyle><mml:mo>&#x00D7;</mml:mo><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mi>D</mml:mi><mml:mi>C</mml:mi><mml:mi>b</mml:mi><mml:mi>u</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x2212;</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi><mml:mi mathvariant="normal">&#x005F;</mml:mi><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>P</italic><sub><italic>d</italic>1</sub> and <italic>P</italic><sub><italic>d</italic>2</sub> are the input and output power of the EDC, and <italic>&#x03BB;</italic><sub><italic>VSCLoss</italic></sub> is the loss parameter of the VSC converter, <italic>Z</italic><sub><italic>DCBranch</italic></sub> is the resistance of the DC path, <italic>V</italic><sub><italic>DCBus</italic></sub> is the DC bus voltage, <italic>P</italic><sub><italic>dcloss</italic></sub> is the power loss of the DC distribution lines.
<list list-type="simple">
<list-item><label>(3)</label><p>Second iteration: During the first calculation of the AC-DC power flow, balance the loads within the substation area using three-phase averages.</p></list-item>
<list-item><label>(4)</label><p>Perform the second AC-DC iteration. If the three phases within the substation area are still unbalanced, introduce new active and reactive power balancing measures and inject power from the tail end to achieve balance.</p></list-item>
<list-item><label>(5)</label><p>Record the iterative values of the AC-DC power flow and compare the differences with the previous iteration to determine whether convergence has occurred.</p></list-item>
</list></p>
</sec>
<sec id="s4">
<label>4</label>
<title>Optimization Operation for EDC-LVDA under Unbalanced Compensation Conditions</title>
<p>Due to the different PV and load conditions at different times, the optimal <italic>k</italic> for the minimal power loss and voltage qualities is different. And the main point of the proposed optimum operation is to dynamically regulate the <italic>k</italic> under different load and PV conditions to achieve the minimum power loss and high-quality voltages. The optimal operation consists of two parts: the optimization model and the optimization solving algorithm. Among them, the optimization model includes two components: the objective function and the constraints. The details of each part are introduced as follows.</p>
<sec id="s4_1">
<label>4.1</label>
<title>Objective Function</title>
<p>The objective of optimization is to achieve a comprehensive optimization of voltage quality and system losses while ensuring full compensation of three-phase imbalances. Thus, the objective function can be represented as a weighted sum of the distribution network losses and voltage quality as follows.
<disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:mo movablelimits="true" form="prefix">min</mml:mo><mml:mi>F</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mfrac><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mfrac><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>s</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula>where <italic>f</italic><sub>1</sub> and <italic>f</italic><sub>2</sub> represent the objective function values for the distribution network loss and voltage quality, respectively. <italic>S</italic><sub>1</sub> and <italic>S</italic><sub>2</sub> are the reference values for <italic>f</italic><sub>1</sub> and <italic>f</italic><sub>2</sub>, which are based on the values of <italic>f</italic><sub>1</sub> and <italic>f</italic><sub>2</sub> when the distribution area has not inserted an EDC part. <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> and <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> are the weights assigned to <italic>f</italic><sub>1</sub> and <italic>f</italic><sub>2</sub>, respectively.</p>
<sec id="s4_1_1">
<label>4.1.1</label>
<title>Objective Function of Power Loss f<sub>1</sub> in EDC-LVDA</title>
<p>The power loss <italic>f</italic><sub><italic>1</italic></sub> in the EDC-LVDA system consists of the VSC converter loss <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, AC path loss <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>A</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, and DC path loss <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>, which can be obtained as
<disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></disp-formula></p>
<p>The loss of each part can be obtained as</p>
<p><disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mi>T</mml:mi></mml:mrow></mml:munder><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>&#x2208;</mml:mo><mml:msup><mml:mi>L</mml:mi><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:munder><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>&#x03D5;</mml:mi><mml:mo>;</mml:mo><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mo>,</mml:mo><mml:mi>b</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo>,</mml:mo><mml:mi>n</mml:mi></mml:mrow></mml:munder><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03D5;</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:msup><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mtext>&#x00A0;</mml:mtext><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mi>T</mml:mi></mml:mrow></mml:munder><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>&#x2208;</mml:mo><mml:msup><mml:mi>L</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:munder><mml:mn>2</mml:mn><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo stretchy="false">)</mml:mo><mml:mrow><mml:mn>2</mml:mn></mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mrow><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo><mml:msubsup><mml:mi>f</mml:mi><mml:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>s</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>a</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mi>T</mml:mi></mml:mrow></mml:munder><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>z</mml:mi><mml:mo>&#x2208;</mml:mo><mml:mi>L</mml:mi></mml:mrow></mml:munder><mml:munder><mml:mo movablelimits="false">&#x2211;</mml:mo><mml:mrow><mml:mi>&#x03D5;</mml:mi><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mo>,</mml:mo><mml:mi>b</mml:mi><mml:mo>,</mml:mo><mml:mi>c</mml:mi><mml:mo>,</mml:mo><mml:mi>n</mml:mi></mml:mrow></mml:munder><mml:mrow></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03D5;</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03D5;</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:mi>t</mml:mi></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where the symbol <italic>Z</italic> denotes the branch number connecting nodes <italic>i</italic> and <italic>j</italic>. <italic>I</italic><sub><italic>z</italic>,<italic>&#x03D5;,t</italic></sub> and <inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03D5;</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> represent the current flowing through a single line in the AC branch and the VSC equivalent loss branch at time <italic>t</italic>. <inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>&#x03D5;</mml:mi></mml:mrow><mml:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the resistance of the VSC equivalent loss branch. <italic>&#x03D5;</italic> indicates the phases <italic>a</italic>, <italic>b</italic>, <italic>c</italic>, and the neutral line n. <inline-formula id="ieqn-8"><mml:math id="mml-ieqn-8"><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the current flowing through the DC branch at time <italic>t</italic>. <italic>R</italic><sub><italic>z</italic>,<italic>&#x03C6;</italic></sub> and <inline-formula id="ieqn-9"><mml:math id="mml-ieqn-9"><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>z</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> are the line resistances of a single line in the AC and DC branches, respectively.</p>
</sec>
<sec id="s4_1_2">
<label>4.1.2</label>
<title>Objective Function of Power Quality f<sub>2</sub> in EDC-LVDA</title>
<p>The power quality for each node can be indicated by the voltage deviation. Thus, <italic>f</italic><sub>2</sub> is specified by the voltage <italic>u</italic><sub><italic>i</italic></sub> of phase <italic>i</italic> at the first end of the station and can be expressed as:
<disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mfrac><mml:mo>&#x00D7;</mml:mo><mml:mn>100</mml:mn><mml:mi mathvariant="normal">&#x0025;</mml:mi></mml:math></disp-formula></p>
<p>where <italic>u</italic><sub><italic>e</italic></sub> is the rated voltage.</p>
</sec>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Constraint Functions</title>
<sec id="s4_2_1">
<label>4.2.1</label>
<title>Operation Power Constraint of VSC</title>
<p>The processed power of VSC must satisfy the following equations as
<disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msup><mml:mi>P</mml:mi><mml:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msup><mml:mi>Q</mml:mi><mml:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>&#x2A7D;</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msup><mml:mi>S</mml:mi><mml:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></disp-formula></p>
<p>where <italic>Q</italic><sup><italic>VSC</italic></sup> and <italic>P</italic><sup><italic>VSC</italic></sup> are the processed active and reactive power of the VSC, respectively. And the <italic>S</italic><sup><italic>VSC</italic></sup> is the power capacity of the VSC.</p>
</sec>
<sec id="s4_2_2">
<label>4.2.2</label>
<title>Power Balance Constraint of the Nodes</title>
<p>The balance constraints of active power and reactive power at node <italic>i</italic> must satisfy the following equation as
<disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>+</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:mrow><mml:mi>L</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</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:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</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:mrow><mml:mi>D</mml:mi><mml:mi>T</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>Q</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>+</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:mrow><mml:mi>L</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</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:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</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:mrow><mml:mi>D</mml:mi><mml:mi>T</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>where <italic>P</italic><sub><italic>in,i,t</italic></sub> and <italic>Q</italic><sub><italic>in,int</italic></sub> are the injected active and reactive power of the node <italic>i</italic> at time <italic>t</italic>, respectively. <inline-formula id="ieqn-10"><mml:math id="mml-ieqn-10"><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:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-11"><mml:math id="mml-ieqn-11"><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:mrow><mml:mi>l</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> denote the active and reactive power of conventional loads at time <italic>t</italic>, respectively. <inline-formula id="ieqn-12"><mml:math id="mml-ieqn-12"><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:mrow><mml:mi>D</mml:mi><mml:mi>T</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-13"><mml:math id="mml-ieqn-13"><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:mrow><mml:mi>D</mml:mi><mml:mi>T</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> represent the active and reactive power injected into the node <italic>i</italic> by the distribution transformer at time <italic>t</italic>. <inline-formula id="ieqn-14"><mml:math id="mml-ieqn-14"><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:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-15"><mml:math id="mml-ieqn-15"><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:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> are the active and reactive power of the voltage source converter (VSC) at time <italic>t</italic>, respectively; <inline-formula id="ieqn-16"><mml:math id="mml-ieqn-16"><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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-17"><mml:math id="mml-ieqn-17"><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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> indicate the active and reactive power generated by PV systems at time <italic>t</italic>.</p>
</sec>
<sec id="s4_2_3">
<label>4.2.3</label>
<title>Current Constraints of the Branches</title>
<p>At time <italic>t</italic>, the current <italic>I</italic><sub><italic>z,t</italic></sub> flowing through branch <italic>z</italic> needs to satisfy.
<disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:mrow><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x2A7D;</mml:mo><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p>where <italic>I</italic><sub><italic>z,max</italic></sub> is the maximum value of the current in branch <italic>z</italic>.</p>
</sec>
<sec id="s4_2_4">
<label>4.2.4</label>
<title>Operation Constraints of the DC System</title>
<p>The power balance of the DC nodes can be expressed as
<disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</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:mrow><mml:mi>L</mml:mi><mml:mi>o</mml:mi><mml:mi>a</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msubsup><mml:mo>=</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:mrow><mml:mi>V</mml:mi><mml:mi>S</mml:mi><mml:mi>C</mml:mi><mml:mo>,</mml:mo><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>+</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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:math></disp-formula></p>
<p>where <inline-formula id="ieqn-18"><mml:math id="mml-ieqn-18"><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>n</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the injected power of node <italic>i</italic> at time <italic>t</italic>, and the <inline-formula id="ieqn-19"><mml:math id="mml-ieqn-19"><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:mrow><mml:mi>v</mml:mi><mml:mi>s</mml:mi><mml:mi>c</mml:mi><mml:mo>.</mml:mo><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the injected power by the VSC of node <italic>i</italic> at time <italic>t</italic>.</p>
<p>The DC voltage at all the time and nodes must be between the allowable range, which can be expressed as
<disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup><mml:mo>&#x2A7D;</mml:mo><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mtext>dc</mml:mtext></mml:mrow></mml:mrow></mml:msubsup><mml:mo>&#x2A7D;</mml:mo><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></disp-formula></p>
<p>where <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the dc voltage of the node <italic>i</italic> at time <italic>t</italic>. <inline-formula id="ieqn-21"><mml:math id="mml-ieqn-21"><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-22"><mml:math id="mml-ieqn-22"><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> are the allowable minimum and maximum values of the dc voltage.</p>
<p>Similarly, the current limitation of each branch can be expressed as
<disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:mo>&#x2212;</mml:mo><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow><mml:mrow><mml:mrow><mml:mtext>dc</mml:mtext></mml:mrow></mml:mrow></mml:msubsup><mml:mo>&#x2A7D;</mml:mo><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mrow><mml:mtext>dc</mml:mtext></mml:mrow></mml:mrow></mml:msubsup><mml:mo>&#x2A7D;</mml:mo><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow><mml:mrow><mml:mrow><mml:mtext>dc</mml:mtext></mml:mrow></mml:mrow></mml:msubsup></mml:math></disp-formula></p>
<p>where <inline-formula id="ieqn-23"><mml:math id="mml-ieqn-23"><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the dc current of the branch <italic>z</italic> at time <italic>t</italic>. <inline-formula id="ieqn-24"><mml:math id="mml-ieqn-24"><mml:msubsup><mml:mi>I</mml:mi><mml:mrow><mml:mi>z</mml:mi><mml:mo>,</mml:mo><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-25"><mml:math id="mml-ieqn-25"><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> are the allowable maximum value of the dc current.</p>
</sec>
<sec id="s4_2_5">
<label>4.2.5</label>
<title>Power Constraints of the PV Nodes</title>
<p>The output power of the PV system must satisfy the following equations as
<disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:msup><mml:mrow><mml:mo>(</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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msup><mml:mrow><mml:mo>(</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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>&#x2A7D;</mml:mo><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:msubsup><mml:mi>S</mml:mi><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></disp-formula></p>
<p>where <inline-formula id="ieqn-26"><mml:math id="mml-ieqn-26"><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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-27"><mml:math id="mml-ieqn-27"><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:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> are the active and reactive power of the PV in node <italic>i</italic> at time <italic>t</italic>, respectively. <inline-formula id="ieqn-28"><mml:math id="mml-ieqn-28"><mml:msubsup><mml:mi>S</mml:mi><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>P</mml:mi><mml:mi>V</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> is the power capacity of the PV inverters.</p>
</sec>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Solving Algorithm</title>
<p>In this paper, a two-layer solving algorithm based on the improved Particle Swarm Optimization (IPSO) algorithm is used for model solving. The upper layer is the power flow model, and the lower layer is the IPSO optimization model. IPSO gradually approaches the optimal solution based on the flying experience of the swarm and the learning process of individuals, achieving optimization of the problem, and is widely applied in the field of power system optimization scheduling [<xref ref-type="bibr" rid="ref-19">19</xref>&#x2013;<xref ref-type="bibr" rid="ref-21">21</xref>]. The detailed steps of IPSO are shown in <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, the details operation steps can be obtained as follows:</p>
<fig id="fig-4">
<label>Figure 4</label>
<caption>
<title>Two-layer solving algorithm process based on the PSO method</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-4.tif"/>
</fig>
<p>Step 1: Solve the upper model. Define decision variables, such as the power factor transferred between substations and the transfer relationships within the substations.</p>
<p>Step 2: Input data such as photovoltaic output and load to calculate the initial values for power flow.</p>
<p>Step 3: Solve the lower-level optimization. Initialize the particle swarm, where each particle represents a potential solution. The initial positions are random but constrained by the feasible region of the decision variables.</p>
<p>Step 4: Set the IPSO parameters, such as inertia weight and individual learning factors.</p>
<p>Step 5: Update using velocity and position updating formulas and calculate the fitness function value.</p>
<p>Step 6: Perform iterative searching, evaluate the objective function value at the new position, and update individual and global best positions.</p>
<p>Step 7: Determine whether the power flow has converged. If not, repeat the cycle. If converged, output the optimization results.</p>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Case Study Analysis</title>
<p>To verify the effectiveness of the proposed optimization scheduling method, an IEEE 13-node network-based EDC-LVDA case study is built and tested. The topology of the EDC-LVDA is shown in <xref ref-type="fig" rid="fig-5">Fig. 5</xref>, in which a DC is embedded between the first and last nodes of the AC system, and the capacity of the DC interconnection converter is 500 kVA. The resistance of the embedded DC line is 0.015 (per unit value, same below), and the converter parameters are R &#x003D; 0.001 and X &#x003D; 0.002. The maximum absolute value of the iterative convergence error is set to be less than &#x03B5; (1 &#x00D7; 10<sup>&#x2212;9</sup>). VSC1 uses DC-bus voltage control, while VSC2 employs split-phase active and reactive power control. The photovoltaic inverter has a capacity of 20 kVA, and the maximum output of the single-phase photovoltaic system is 10 kW. The output power curves of the photovoltaic nodes at different times are shown in <xref ref-type="fig" rid="fig-6">Fig. 6</xref>. <xref ref-type="table" rid="table-1">Table 1</xref> summarizes the connected rated power under different nodes, and the load utilization factor of each node at various times is shown in <xref ref-type="fig" rid="fig-7">Fig. 7</xref>.</p>
<fig id="fig-5">
<label>Figure 5</label>
<caption>
<title>The topology of the IEEE 13-node network based EDC-LVDA case study</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-5.tif"/>
</fig><fig id="fig-6">
<label>Figure 6</label>
<caption>
<title>Output power of each PV node</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-6.tif"/>
</fig><table-wrap id="table-1">
<label>Table 1</label>
<caption>
<title>The connected rated power loads under different nodes</title>
</caption>
<table>
<colgroup>
<col/>
<col/>
<col/>
<col/>
</colgroup>
<thead>
<tr>
<th>Nodes</th>
<th>Phase A (kW <bold>&#x002B;</bold> jkvar)</th>
<th>Phase B (kW <bold>&#x002B;</bold> jkvar)</th>
<th>Phase C (kW <bold>&#x002B;</bold> jkvar)</th>
</tr>
</thead>
<tbody>
<tr>
<td>2</td>
<td>10 &#x002B; j2</td>
<td>15 &#x002B; j2</td>
<td>5 &#x002B; j2</td>
</tr>
<tr>
<td>3</td>
<td>10 &#x002B; j4</td>
<td>16 &#x002B; j6</td>
<td>16 &#x002B; j5</td>
</tr>
<tr>
<td>4</td>
<td>10 &#x002B; j4</td>
<td>12 &#x002B; j4</td>
<td>12 &#x002B; j5</td>
</tr>
<tr>
<td>5</td>
<td>/</td>
<td>17 &#x002B; j6</td>
<td>20 &#x002B; j7</td>
</tr>
<tr>
<td>6</td>
<td>/</td>
<td>23 &#x002B; j10</td>
<td>/</td>
</tr>
<tr>
<td>7</td>
<td>10 &#x002B; j5</td>
<td>20 &#x002B; j10</td>
<td>15 &#x002B; j10</td>
</tr>
<tr>
<td>8</td>
<td>5 &#x002B; j2</td>
<td>/</td>
<td>/</td>
</tr>
<tr>
<td>9</td>
<td>/</td>
<td>6 &#x002B; j2</td>
<td>17 &#x002B; j8</td>
</tr>
<tr>
<td>10</td>
<td>12 &#x002B; j7</td>
<td>/</td>
<td>/</td>
</tr>
<tr>
<td>11</td>
<td>20 &#x002B; j11</td>
<td>10 &#x002B; j6</td>
<td>15 &#x002B; j8</td>
</tr>
<tr>
<td>12</td>
<td>10 &#x002B; j4</td>
<td>5 &#x002B; j3</td>
<td>15 &#x002B; j7</td>
</tr>
<tr>
<td>13</td>
<td>20 &#x002B; j8</td>
<td>15 &#x002B; j5</td>
<td>10 &#x002B; j6</td>
</tr>
</tbody>
</table>
</table-wrap><fig id="fig-7">
<label>Figure 7</label>
<caption>
<title>The load utilization factor at different times</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-7.tif"/>
</fig>
<p><xref ref-type="fig" rid="fig-8">Fig. 8</xref> shows the convergence curves of power flow and the optimization method. <xref ref-type="fig" rid="fig-8">Fig. 8a</xref> shows the convergence curves of maximum absolute difference <inline-formula id="ieqn-29"><mml:math id="mml-ieqn-29"><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi mathvariant="normal">&#x005F;</mml:mi><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>a</mml:mi><mml:mi mathvariant="normal">&#x005F;</mml:mi><mml:mrow><mml:mi>o</mml:mi><mml:mi>l</mml:mi><mml:mi>d</mml:mi><mml:mi>i</mml:mi></mml:mrow></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">|</mml:mo></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> in the power flow calculation process. It can be seen that the power flow can converge quickly and meet the convergence requirements. Based on the above parameters in the EDC-LVDA system, the conventional PSO, genetic algorithm (GA), and the proposed IPSO-based algorithm are respectively used to solve the objective and optimized model functions proposed in this paper, and the convergence curves are shown in <xref ref-type="fig" rid="fig-8">Fig. 8b</xref>. The IPSO-based algorithm can converge quickly and obtain the minimum optimal convergence value of 39.1932, showing a better optimization effect.</p>
<fig id="fig-8">
<label>Figure 8</label>
<caption>
<title>Convergence curves of power flow and optimization method. (<bold>a</bold>) Power Plow Calculation process; (<bold>b</bold>) optimization process under different algorithms</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-8.tif"/>
</fig>
<p>The optimal power distribution ratio <italic>k</italic> during an entire day has been obtained using the proposed optimization scheduling algorithm, as shown in <xref ref-type="fig" rid="fig-9">Fig. 9</xref>. The corresponding total system loss is depicted in <xref ref-type="fig" rid="fig-10">Fig. 10</xref>. It can be observed that compared to a fixed power distribution ratio, the proposed method can dynamically calculate the optimal power distribution ratio <italic>k</italic>, obtaining a significantly lower network loss, with a total loss reduction of 13.8% in a day.</p>
<fig id="fig-9">
<label>Figure 9</label>
<caption>
<title>The optimized power distribution ratio</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-9.tif"/>
</fig><fig id="fig-10">
<label>Figure 10</label>
<caption>
<title>Power loss comparison</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-10.tif"/>
</fig>
<p><xref ref-type="fig" rid="fig-11">Fig. 11</xref> shows the active and reactive power data of each phase at the transformer of the distribution area under the optimal power distribution ratio <italic>k</italic>. With the proposed optimization scheduling method, the power of each phase can maintain a three-phase power balance despite fluctuations in load and photovoltaic power. This sufficiently demonstrates that the proposed EDC-LVDA scheduling strategy can achieve balanced load distribution in the distribution area even under unbalanced load conditions.</p>
<fig id="fig-11">
<label>Figure 11</label>
<caption>
<title>The active and reactive power data of each phase at the transformer and the times</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-11.tif"/>
</fig>
<p><xref ref-type="fig" rid="fig-12">Fig. 12</xref> shows the voltage conditions of each node in the distribution area for each phase. The voltages at all nodes remain within the acceptable standard range of [0.95, 1.05], indicating that the proposed method can effectively address the low voltage issues in the distribution area and ensure the voltage quality of the system.</p>
<fig id="fig-12">
<label>Figure 12</label>
<caption>
<title>The voltage conditions of each node in the distribution area for each phase</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_69610-fig-12.tif"/>
</fig>
</sec>
<sec id="s6">
<label>6</label>
<title>Conclusions</title>
<p>Aiming to address the challenge of power optimization scheduling within an EDC-LVDA, this paper proposes an optimization scheduling strategy for interlinked flexible distribution networks that takes into account three-phase unbalance compensation. Firstly, a power flow calculation method for low-voltage flexible distribution areas that considers three-phase unbalance compensation is proposed. This method can effectively describe the power distribution characteristics under both AC and DC power allocation. Secondly, an optimization scheduling model for flexible low-voltage distribution areas has been constructed, which accounts for network losses, voltage quality, losses in the DC section, and three-phase imbalance. This model can adaptively select the power distribution coefficients for the DC and AC sections under different operating conditions, achieving embedded economic optimization while ensuring compensation for unbalanced conditions. Finally, a case study based on the IEEE 13 bus system is established to validate the proposed method. The results indicate that the proposed method can achieve 100% voltage compliance rate and a 13.8% reduction in network losses, which can support the integration of distributed energy, adapt to new loads such as electric vehicles, strengthen grid resilience, and provide a practical solution for the modernization of low-voltage power grids.</p>
<p>Besides, the proposed approach is designed for the optimal operation of EDC interconnections within a single area. While its current version does not apply to EDC interconnections between different areas, the extension of this approach to cross-area EDC interconnection scenarios will be elaborated in our future publications.</p>
</sec>
</body>
<back>
<ack>
<p>None.</p>
</ack>
<sec>
<title>Funding Statement</title>
<p>This work was supported by the key technology project of China Southern Power Grid Corporation (GZKJXM20220041), and partly by the National Key Research and Development Plan (2022YFE0205300).</p>
</sec>
<sec>
<title>Author Contributions</title>
<p>The authors confirm contribution to the paper as follows: Conceptualization, Zhukui Tan; methodology and software, Dacheng Zhou and Song Deng; data curation, Jikai Li and Zhuang Wu; writing&#x2014;original draft preparation, Dacheng Zhou; writing&#x2014;review and editing, Qihui Feng; supervision, Xuan Zhang and Zhukui Tan. All authors reviewed the results and approved the final version of the manuscript.</p>
</sec>
<sec sec-type="data-availability">
<title>Availability of Data and Materials</title>
<p>The original contributions presented in the study are included in the article; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec>
<title>Ethics Approval</title>
<p>Not applicable.</p>
</sec>
<sec sec-type="COI-statement">
<title>Conflicts of Interest</title>
<p>The authors declare no conflicts of interest to report regarding the present study.</p>
</sec>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>[1]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wei</surname> <given-names>YM</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>K</given-names></string-name>, <string-name><surname>Kang</surname> <given-names>JN</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>W</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>XY</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>X</given-names></string-name></person-group>. <article-title>Policy and management of carbon peaking and carbon neutrality: a literature review</article-title>. <source>Engineering</source>. <year>2022</year>;<volume>8</volume>(<issue>7</issue>):<fpage>52</fpage>&#x2013;<lpage>63</lpage> doi:<pub-id pub-id-type="doi">10.1016/j.eng.2021.12.018</pub-id>.</mixed-citation></ref>
<ref id="ref-2"><label>[2]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname> <given-names>M</given-names></string-name>, <string-name><surname>Long</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Guo</surname> <given-names>S</given-names></string-name>, <string-name><surname>Xiao</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Shi</surname> <given-names>T</given-names></string-name>, <string-name><surname>Xiang</surname> <given-names>X</given-names></string-name>, <etal>et al</etal></person-group>. <article-title>Stability boundary characterization and power quality improvement for distribution networks</article-title>. <source>Energies</source>. <year>2024</year>;<volume>17</volume>(<issue>24</issue>):<fpage>6215</fpage>. doi:<pub-id pub-id-type="doi">10.3390/en17246215</pub-id>.</mixed-citation></ref>
<ref id="ref-3"><label>[3]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Amrani</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Beladel</surname> <given-names>A</given-names></string-name>, <string-name><surname>Kouzou</surname> <given-names>A</given-names></string-name>, <string-name><surname>Rodriguez</surname> <given-names>J</given-names></string-name>, <string-name><surname>Abdelrahem</surname> <given-names>M</given-names></string-name></person-group>. <article-title>Four-wire three-level NPC shunt active power filter using model predictive control based on the grid-tied PV system for power quality enhancement</article-title>. <source>Energies</source>. <year>2024</year>;<volume>17</volume>(<issue>15</issue>):<fpage>3822</fpage>. doi:<pub-id pub-id-type="doi">10.3390/en17153822</pub-id>.</mixed-citation></ref>
<ref id="ref-4"><label>[4]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname> <given-names>G</given-names></string-name>, <string-name><surname>Li</surname> <given-names>J</given-names></string-name>, <string-name><surname>Yan</surname> <given-names>K</given-names></string-name>, <string-name><surname>Bian</surname> <given-names>J</given-names></string-name></person-group>. <article-title>Centralized-distributed scheduling strategy of distribution network based on multi-temporal hierarchical cooperative game</article-title>. <source>Energy Eng</source>. <year>2025</year>;<volume>122</volume>(<issue>3</issue>):<fpage>1113</fpage>&#x2013;<lpage>36</lpage>. doi:<pub-id pub-id-type="doi">10.32604/ee.2025.059558</pub-id>.</mixed-citation></ref>
<ref id="ref-5"><label>[5]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Guo</surname> <given-names>P</given-names></string-name>, <string-name><surname>Tian</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Yuan</surname> <given-names>Z</given-names></string-name>, <string-name><surname>Qu</surname> <given-names>L</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>XY</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>XP</given-names></string-name></person-group>. <article-title>Future-proofing city power grids: fid-based efficient interconnection strategies for major load-centred environments</article-title>. <source>IET Renew Power Gener</source>. <year>2024</year>;<volume>18</volume>(<issue>15</issue>):<fpage>3003</fpage>&#x2013;<lpage>19</lpage>. doi:<pub-id pub-id-type="doi">10.1049/rpg2.13027</pub-id>.</mixed-citation></ref>
<ref id="ref-6"><label>[6]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Mtolo</surname> <given-names>D</given-names></string-name>, <string-name><surname>Dorrell</surname> <given-names>D</given-names></string-name>, <string-name><surname>Pillay Carpanen</surname> <given-names>R</given-names></string-name></person-group>. <article-title>Balancing of low-voltage supply network with a smart utility controller leveraging distributed customer energy sources</article-title>. <source>Energies</source>. <year>2023</year>;<volume>16</volume>(<issue>23</issue>):<fpage>7707</fpage>. doi:<pub-id pub-id-type="doi">10.3390/en16237707</pub-id>.</mixed-citation></ref>
<ref id="ref-7"><label>[7]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Guo</surname> <given-names>X</given-names></string-name>, <string-name><surname>Gan</surname> <given-names>L</given-names></string-name>, <string-name><surname>Liu</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Hu</surname> <given-names>W</given-names></string-name></person-group>. <article-title>Coordinated dispatching of flexible AC/DC distribution areas considering uncertainty of source and load</article-title>. <source>Electr Power Syst Res</source>. <year>2023</year>;<volume>225</volume>(<issue>042&#x2013;011</issue>):<fpage>109805</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.epsr.2023.109805</pub-id>.</mixed-citation></ref>
<ref id="ref-8"><label>[8]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Gao</surname> <given-names>C</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>C</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>Y</given-names></string-name>, <string-name><surname>Qu</surname> <given-names>R</given-names></string-name>, <string-name><surname>Niu</surname> <given-names>K</given-names></string-name>, <string-name><surname>Chen</surname> <given-names>W</given-names></string-name></person-group>. <article-title>Novel low-carbon optimal operation method for flexible distribution network based on carbon emission flow</article-title>. <source>Energy Eng</source>. <year>2025</year>;<volume>122</volume>(<issue>2</issue>):<fpage>785</fpage>&#x2013;<lpage>803</lpage>. doi:<pub-id pub-id-type="doi">10.32604/ee.2024.058705</pub-id>.</mixed-citation></ref>
<ref id="ref-9"><label>[9]</label><mixed-citation publication-type="book"><person-group person-group-type="author"><string-name><surname>Li</surname> <given-names>P</given-names></string-name>, <string-name><surname>Ji</surname> <given-names>H</given-names></string-name>, <string-name><surname>Yu</surname> <given-names>H</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>C</given-names></string-name></person-group>. <source>Flexible distribution networks</source>. <publisher-loc>Amsterdam, The Netherlands</publisher-loc>: <publisher-name>Elsevier</publisher-name>; <year>2024</year>. p. <fpage>3</fpage>&#x2013;<lpage>10</lpage>. doi:<pub-id pub-id-type="doi">10.1016/B978-0-12-823890-5.00002-7</pub-id>.</mixed-citation></ref>
<ref id="ref-10"><label>[10]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tan</surname> <given-names>G</given-names></string-name>, <string-name><surname>Guo</surname> <given-names>X</given-names></string-name>, <string-name><surname>Zhao</surname> <given-names>W</given-names></string-name>, <string-name><surname>Qi</surname> <given-names>L</given-names></string-name>, <string-name><surname>Sun</surname> <given-names>X</given-names></string-name></person-group>. <article-title>Second harmonic suppression for DC output voltage of three-phase four-leg PWM rectifier under unbalanced grid voltage conditions</article-title>. <source>IEEE Trans Power Electron</source>. <year>2025</year>;<volume>40</volume>(<issue>8</issue>):<fpage>11088</fpage>&#x2013;<lpage>106</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TPEL.2025.3551801</pub-id>.</mixed-citation></ref>
<ref id="ref-11"><label>[11]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Pereira</surname> <given-names>O</given-names></string-name>, <string-name><surname>Quir&#x00F3;s-Tort&#x00F3;s</surname> <given-names>J</given-names></string-name>, <string-name><surname>Valverde</surname> <given-names>G</given-names></string-name></person-group>. <article-title>Phase rebalancing of distribution circuits dominated by single-phase loads</article-title>. <source>IEEE Trans Power Syst</source>. <year>2021</year>;<volume>36</volume>(<issue>6</issue>):<fpage>5333</fpage>&#x2013;<lpage>44</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TPWRS.2021.3076629</pub-id>.</mixed-citation></ref>
<ref id="ref-12"><label>[12]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Wang</surname> <given-names>J</given-names></string-name>, <string-name><surname>Zhou</surname> <given-names>N</given-names></string-name>, <string-name><surname>Chung</surname> <given-names>CY</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>Q</given-names></string-name></person-group>. <article-title>Coordinated planning of converter-based DG units and soft open points incorporating active management in unbalanced distribution networks</article-title>. <source>IEEE Trans Sustain Energy</source>. <year>2020</year>;<volume>11</volume>(<issue>3</issue>):<fpage>2015</fpage>&#x2013;<lpage>27</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TSTE.2019.2950168</pub-id>.</mixed-citation></ref>
<ref id="ref-13"><label>[13]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname> <given-names>B</given-names></string-name>, <string-name><surname>Tang</surname> <given-names>W</given-names></string-name>, <string-name><surname>Cong</surname> <given-names>PW</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>XH</given-names></string-name>, <string-name><surname>Lou</surname> <given-names>CW</given-names></string-name></person-group>. <article-title>Multi-time scale optimal control in hybrid AC/DC distribution networks based on SOP and VSC</article-title>. <source>Adv Technol Electr Eng Energy</source>. <year>2017</year>;<volume>36</volume>(<issue>9</issue>):<fpage>11</fpage>&#x2013;<lpage>9</lpage>. (In Chinese). doi:<pub-id pub-id-type="doi">10.1049/joe.2018.8527</pub-id>.</mixed-citation></ref>
<ref id="ref-14"><label>[14]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhang</surname> <given-names>B</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>L</given-names></string-name>, <string-name><surname>Tang</surname> <given-names>W</given-names></string-name>, <string-name><surname>Li</surname> <given-names>G</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>C</given-names></string-name></person-group>. <article-title>Optimal planning of hybrid AC/DC low-voltage distribution networks considering DC conversion of three-phase four-wire low-voltage AC systems</article-title>. <source>J Mod Power Syst Clean Energy</source>. <year>2024</year>;<volume>12</volume>(<issue>1</issue>):<fpage>141</fpage>&#x2013;<lpage>53</lpage>. doi:<pub-id pub-id-type="doi">10.35833/mpce.2022.000404</pub-id>.</mixed-citation></ref>
<ref id="ref-15"><label>[15]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ge</surname> <given-names>L</given-names></string-name>, <string-name><surname>Zhang</surname> <given-names>W</given-names></string-name>, <string-name><surname>Yan</surname> <given-names>F</given-names></string-name>, <string-name><surname>Yuan</surname> <given-names>XD</given-names></string-name>, <string-name><surname>Yu</surname> <given-names>YZ</given-names></string-name></person-group>. <article-title>Optimal scheduling of flexible interconnected distribution network based on adaptive model predictive control</article-title>. <source>Electr Power Autom Equip</source>. <year>2020</year>;<volume>40</volume>(<issue>6</issue>):<fpage>15</fpage>&#x2013;<lpage>23</lpage>. (In Chinese). doi:<pub-id pub-id-type="doi">10.16081/j.epae.202005009</pub-id>.</mixed-citation></ref>
<ref id="ref-16"><label>[16]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname> <given-names>T</given-names></string-name>, <string-name><surname>Guo</surname> <given-names>Q</given-names></string-name>, <string-name><surname>Wang</surname> <given-names>X</given-names></string-name>, <string-name><surname>Ji</surname> <given-names>H</given-names></string-name>, <string-name><surname>Tu</surname> <given-names>C</given-names></string-name>, <string-name><surname>Hou</surname> <given-names>Y</given-names></string-name>, <etal>et al.</etal></person-group> <article-title>A comprehensive control strategy for three phase unbalance and light and heavy load in flexible interconnected distribution stations based on F-SOP</article-title>. <source>Power Syst Technol</source>. <year>2024</year>;<volume>48</volume>(<issue>10</issue>):<fpage>4358</fpage>&#x2013;<lpage>66</lpage>. (In Chinese). doi:<pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2023.1268</pub-id>.</mixed-citation></ref>
<ref id="ref-17"><label>[17]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Feng</surname> <given-names>Q</given-names></string-name>, <string-name><surname>Li</surname> <given-names>J</given-names></string-name>, <string-name><surname>Zhou</surname> <given-names>J</given-names></string-name>, <string-name><surname>Xu</surname> <given-names>Y</given-names></string-name>, <string-name><surname>He</surname> <given-names>Q</given-names></string-name>, <string-name><surname>Gao</surname> <given-names>Y</given-names></string-name></person-group>. <article-title>A coordinated power quality improvement control strategy for AC/DC hybrid distribution networks based on three-phase four-leg flexible interconnection converter</article-title>. <source>Front Energy Res</source>. <year>2025</year>;<volume>13</volume>:<fpage>1602269</fpage>. doi:<pub-id pub-id-type="doi">10.3389/fenrg.2025.1602269</pub-id>.</mixed-citation></ref>
<ref id="ref-18"><label>[18]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Bak&#x0131;r</surname> <given-names>H</given-names></string-name>, <string-name><surname>Guvenc</surname> <given-names>U</given-names></string-name>, <string-name><surname>Duman</surname> <given-names>S</given-names></string-name>, <string-name><surname>Kahraman</surname> <given-names>HT</given-names></string-name></person-group>. <article-title>Optimal power flow for hybrid AC/DC electrical networks configured with VSC-MTDC transmission lines and renewable energy sources</article-title>. <source>IEEE Syst J</source>. <year>2023</year>;<volume>17</volume>(<issue>3</issue>):<fpage>3938</fpage>&#x2013;<lpage>49</lpage>. doi:<pub-id pub-id-type="doi">10.1109/JSYST.2023.3248658</pub-id>.</mixed-citation></ref>
<ref id="ref-19"><label>[19]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Bak&#x0131;r</surname> <given-names>H</given-names></string-name>, <string-name><surname>Kahraman</surname> <given-names>HT</given-names></string-name>, <string-name><surname>Y&#x0131;lmaz</surname> <given-names>S</given-names></string-name>, <string-name><surname>Duman</surname> <given-names>S</given-names></string-name>, <string-name><surname>Guvenc</surname> <given-names>U</given-names></string-name></person-group>. <article-title>Dynamic switched crowding-based multi-objective particle swarm optimization algorithm for solving multi-objective AC-DC optimal power flow problem</article-title>. <source>Appl Soft Comput</source>. <year>2024</year>;<volume>166</volume>(<issue>3</issue>):<fpage>112155</fpage>. doi:<pub-id pub-id-type="doi">10.1016/j.asoc.2024.112155</pub-id>.</mixed-citation></ref>
<ref id="ref-20"><label>[20]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhu</surname> <given-names>H</given-names></string-name>, <string-name><surname>Liu</surname> <given-names>T</given-names></string-name></person-group>. <article-title>Rotor displacement self-sensing modeling of six-pole radial hybrid magnetic bearing using improved particle swarm optimization support vector machine</article-title>. <source>IEEE Trans Power Electron</source>. <year>2020</year>;<volume>35</volume>(<issue>11</issue>):<fpage>12296</fpage>&#x2013;<lpage>306</lpage>. doi:<pub-id pub-id-type="doi">10.1109/TPEL.2020.2982746</pub-id>.</mixed-citation></ref>
<ref id="ref-21"><label>[21]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Belgana</surname> <given-names>S</given-names></string-name>, <string-name><surname>Fortin-Blanchette</surname> <given-names>H</given-names></string-name></person-group>. <article-title>A novel neural network-based droop control strategy for single-phase power converters</article-title>. <source>Energies</source>. <year>2024</year>;<volume>17</volume>(<issue>23</issue>):<fpage>5825</fpage>. doi:<pub-id pub-id-type="doi">10.3390/en17235825</pub-id>.</mixed-citation></ref>
</ref-list>
</back></article>