<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">869844</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2022.869844</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>A Hierarchical Strategy for Multi-Objective Optimization of Distribution Network Considering DGs and V2G-Enabled EVs Integration</article-title>
<alt-title alt-title-type="left-running-head">Huang et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Multi-Objective Optimization of Distribution Network</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Huang</surname>
<given-names>Zhao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1550009/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ma</surname>
<given-names>Pengbo</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Mengmeng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Fang</surname>
<given-names>Baling</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1535278/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Ming</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/744885/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Railway Transportation</institution>, <institution>Hunan University of Technology</institution>, <addr-line>Zhuzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Key Laboratory for Micro-Nano Optoelectronic Devices of Ministry of Education</institution>, <institution>School of Physics and Electronics</institution>, <institution>Hunan University</institution>, <addr-line>Changsha</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/1259467/overview">Liansong Xiong</ext-link>, Nanjing Institute of Technology (NJIT), 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/1668472/overview">Long Zeng</ext-link>, Guangdong University of Technology, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1669713/overview">Pu Liu</ext-link>, Zhengzhou University of Light Industry, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Zhao Huang, <email>skyhuang1985@hotmail.com</email>; Baling Fang, <email>5911866@qq.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>10</day>
<month>03</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>869844</elocation-id>
<history>
<date date-type="received">
<day>05</day>
<month>02</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>18</day>
<month>02</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Huang, Ma, Wang, Fang and Zhang.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Huang, Ma, Wang, Fang and Zhang</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&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>As an important part of smart grid construction, the distribution network (DN) optimization problem has always attracted great attention, especially under the background that large-scale penetration of distributed generators (DGs) and electric vehicles (EVs) into building cluster poses both opportunities and challenges to the energy management. This research presents a hierarchical optimization strategy, for improving the safe and economical operation of DN considering the DGs and EVs integration. In Stage 1, the MPPT control model of DGs is designed to obtain the best energy conversion efficiency. In Stage 2, load models of EVs and battery energy storage system (BESS) under coordinate charging/discharging stimulated by a time-of-use incentive mechanism are established respectively, to achieve a load curve with a minimized peak-to-valley difference (PVD). In Stage 3, aiming for the best compromise between the active power loss and node voltage excursion, daily optimal scheduling of the static Var compensator (SVC) capacitors is dynamically worked out according to the varying power demand, as the solution for the defined multi-objective optimization problem. For enhancing the convergence speed, an advanced genetic algorithm with elite preservation strategy is employed. The proposed hierarchical strategy is demonstrated on an IEEE 33-node DN test case, and the simulation results show that first, the MPPT control ensures the maximum power outputs of DGs; next, power supply pressure could be relieved by the load shifting effects of the coordinated vehicle-to-grid (V2G) service and BESS configuration, reflected in the decreased load peak from 4,370.1 to 3,424.99&#xa0;kW, and the optimized PVD from 1763.8 to 703.8&#xa0;kW; meanwhile, <italic>via</italic> applicable power planning of the SVC components, optimized power loss and voltage quality can both be achieved, proving the feasibility of the optimization strategy, which promotes the economic and reliable operation of the DN system.</p>
</abstract>
<kwd-group>
<kwd>distribution network</kwd>
<kwd>distributed generators</kwd>
<kwd>electric vehicles</kwd>
<kwd>vehicle-to-grid</kwd>
<kwd>multi-objective optimization</kwd>
<kwd>coordinated dispatch</kwd>
</kwd-group>
<contract-num rid="cn002">20C0613</contract-num>
<contract-sponsor id="cn001">Natural Science Foundation of Hunan Province<named-content content-type="fundref-id">10.13039/501100004735</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Scientific Research Foundation of Hunan Provincial Education Department<named-content content-type="fundref-id">10.13039/100014472</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>With the sustainable development of the energy industry, the excessive exploitation of traditional fossil energy has brought increasingly severe environmental pollution problems such as energy resource shortage and climate change, making it extremely urgent to utilize renewable energy as an alternative for fossil fuels (<xref ref-type="bibr" rid="B2">Brockway et&#x20;al., 2019</xref>; <xref ref-type="bibr" rid="B22">Xiong et&#x20;al., 2020</xref>). In recent years, the large-scale penetration of renewable resources and new energy equipment alleviates the energy crisis, while also bringing rigorous challenges to the power system optimizing configuration and secondary energy rational utilization (<xref ref-type="bibr" rid="B16">Peidong et&#x20;al., 2009</xref>). For example, the randomness and intermittency of power generated from wind turbines (WTs) or photovoltaic cells (PVs) and the time-space dimensional decentralization of electric vehicle (EV) charging, induce more uncertainty for present distribution network (DN) operations than ever before, as well as problems involving harmonic pollution (<xref ref-type="bibr" rid="B19">Siahroodi et&#x20;al., 2021</xref>), three-phase voltage imbalance (<xref ref-type="bibr" rid="B10">Islam et&#x20;al., 2020</xref>), and transformers aging (<xref ref-type="bibr" rid="B5">Elbatawy and Morsi, 2022</xref>). Particularly, the enlarged peak-to-valley difference (PVD) of the load curve puts forward higher requirements for the power system. The safety and economy, expressed as node voltage excursion and active power loss respectively, are commonly regarded as crucial indices reflecting the operation state of the DN (<xref ref-type="bibr" rid="B4">Dong et&#x20;al., 2019</xref>). The serious peak load arisen from the aggregated charging behaviors of EV owners probably causes low-voltage even blackout, concurrently increasing system loss. Hence, guiding measures for coordinated regulation should be taken to minimize the negative influences of integrated renewable energy and EVs while satisfying the power and travel demands of customers.</p>
<p>Recently, the integrated energy system (IES) has been greatly supported by the Chinese government and turned into a research hotspot, since it benefits the integration of renewable energy and the coordinated development of the multi-energy system. Under the background of the IES, the management idea for the DN has changed from supply-oriented to demand-oriented (<xref ref-type="bibr" rid="B27">Zhang et&#x20;al., 2018</xref>), which integrates the functions of distributed generation control, real-time monitoring, information sharing, and market transactions. This makes for better accommodation and scheduling abilities to diversify distributed energy equipment.</p>
<p>To stabilize the output fluctuation of renewable power generation, the battery energy storage system (BESS) has become a key supporting unit to improve the compatibility level for the WT and PV of the DN. By power transferring through electrochemical charging/discharging, the BESS not only enhances the utilization of clean energies, bringing higher economic benefits to the regional grid but further improves the reliability of the power supply. A hierarchical coordinated control strategy with the fast-response BESS to suppress high fluctuations associated with DGs outputs has been developed based on the MPC framework (<xref ref-type="bibr" rid="B28">Zhang et&#x20;al., 2021</xref>). Moreover, as part of a vehicle-to-grid (V2G) system (<xref ref-type="bibr" rid="B14">Luo et&#x20;al., 2020</xref>), EV can also be regarded as a controllable power resource, to realize the bidirectional power flow between the building cluster and DN. According to previous studies (<xref ref-type="bibr" rid="B26">Zeng et&#x20;al., 2020</xref>; <xref ref-type="bibr" rid="B6">Fang et&#x20;al., 2021</xref>), the time-of-use (TOU) price mechanism is an effective method for integrated demand side management, which creates an economic incentive for users to adjust their charging/discharging time. For example, residential customers plug in their vehicles during off-peak hours and discharge in peak regions to produce V2G profits. This helps to reduce the risks of distribution transformer overloading and power outage and stabilizes the load fluctuation by valley filling and peak shaving, thus balancing the demand and supply. These findings could enlighten us about solving the optimizing configuration problem of the EV-integrated&#x20;DN.</p>
<p>Tremendous efforts have been devoted to DN optimization referring to multi-source cooperation. A voltage regulation model was proposed in a study by <xref ref-type="bibr" rid="B15">Mahmoudi et&#x20;al. (2021)</xref> to enhance the voltage quality of the power system with EVs penetration; nevertheless, the active loss is not included in the objective function. Contrarily, <xref ref-type="bibr" rid="B3">Chen G. et al. (2013)</xref> analyzed the model for DN reconfiguration, with the single optimization function of power loss. Zarei M et&#x20;al. designed a multi-objective optimization model for DN reconfiguration, considering the uncertainties of load variations and power production of DGs (<xref ref-type="bibr" rid="B25">Zarei and Zangeneh, 2017</xref>). For simultaneous optimization of power loss and voltage stability of the DN, reactive compensation measures, carried out through static Var compensator (SVC) capacitor generally in studies by (<xref ref-type="bibr" rid="B23">Xu et&#x20;al., 2021</xref>) and (<xref ref-type="bibr" rid="B7">GholamiFarkoush et&#x20;al., 2019</xref>), are always taken. Studies by (<xref ref-type="bibr" rid="B21">Wu et&#x20;al., 2021</xref>) and (<xref ref-type="bibr" rid="B20">Tan and Chen, 2020</xref>) mainly focus on power management under the predictive control strategy for charging EVs, while they are only regarded as stochastic charging loads rather than controllable units which could participate in energy dispatch, and the V2G functions are not included. In a study by <xref ref-type="bibr" rid="B11">Jiao et&#x20;al. (2021)</xref>, a multi-objective optimization model emphasizing the EV integration to a smart DN is described, with the EVs scheduled in response to TOU price. Since the charging or discharging plan is generated with fixed peak and valley periods, it may not be suitable for the time-varying&#x20;loads.</p>
<p>As analyzed above, earlier reports on the DN optimization operation concentrated on single or combinative objective functions of node voltage fluctuation and active power loss; however, these works entailed some deficiencies including the following aspects: (1) optimization scheduling did not involve wide-range types of energy resources, such as EVs or DGs, which were disregarded in some previous works; (2) utilization efficiency perfection of DGs was ignored; (3) EVs were only treated as stochastic charging loads; and (4) the dynamic dispatch of EVs was not considered.</p>
<p>To overcome the shortcomings of previous works, a hierarchical optimization strategy is proposed in this study, which combines the maximum power point tracking (MPPT) control of the DGs (Stage 1), coordinated V2G service, and BESS dispatch for minimized PVD <italic>via</italic> the positive guidance of the TOU mechanism (Stage 2), with dynamic reactive compensation for maintaining the reliable and economic operation of the DN (Stage 3). The tradeoff between node voltage excursion and active power loss is selected as the objective function to be optimized, and an improved genetic algorithm (GA) with elite preservation strategy (EPS) owning fast convergence is adopted to solve the best daily scheduling plan for the SVC component while meeting all constraints of the DN. An IEEE 33-bus system involving WT, PV, BESS, and EVs accessed to the building cluster is taken as the research object, to validate the effectiveness of the optimization model. In Stage 1, MPPT control improves the energy conversion efficiency of DGs and reduces the power supply pressure of the DN. In Stage 2, EVs participate in energy management as flexible loads, and functions including stabilization of the DGs&#x2019; output fluctuation and power load shifting are realized through coordinated dispatch of EVs and the BESS. In Stage 3, the multi-objective economic-secure optimization problem is dynamically worked out through SVC being put into operation, and the significant advantage of load shifting is highlighted by the statistical result comparison of relative voltage fluctuation and power loss, further confirming the effectiveness of the hierarchical scheduling approach.</p>
</sec>
<sec id="s2">
<title>2 A Hierarchical Optimization Framework</title>
<p>A smart DN is mostly configured as an energy cascade utilization system, enabling free power flow among multiple energy sources, and provides power to residents through building cluster. <xref ref-type="fig" rid="F1">Figure&#x20;1</xref> briefly illustrates the structure of the DN containing dispersed multi-source. Specifically, the demand side involves basic power load, BESS to be charged, and scattered EVs connected to the building cluster. The schedulable supply side includes diversified distributed energy sources (such as WT, PV, bidirectional EV chargers, and BESS), to satisfy the power demand of residents. EVs are not regarded as traditional loads since they can also act as power suppliers. For the V2G application, EV devices not only provide customers with charging services but also allow the owners to interact with the power grid and make V2G profits by selling excessive energy back to the grid at an appointed electricity price level. Overall, the optimization strategy of the multi-source DN is intended to maintain stability and improve the energy efficiency of the whole system through power control over the distributed energy resources based on their operation features. On the premise of that daily power demand, information is collected, and the impacts of EV traveling and distributed energy equipment working are confirmed; the efficient planning and operation of the supply side can be carried out according to users&#x2019; willingness, under the guidance of distributed cooperative control, to finally achieve the comprehensive optimization dispatch from the demand side response, energy interaction, and DN operation perspectives. In this study, a hierarchical optimization strategy for the smart DN with multi-source integration is proposed, and the framework is presented in <xref ref-type="fig" rid="F2">Figure&#x20;2</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Schematic of the multi-source DN.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g001.tif"/>
</fig>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Block diagram of the hierarchical optimization procedure.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g002.tif"/>
</fig>
<p>In Stage 1, the optimal operation control of distributed renewable energy generators is implemented. MPPT measures are taken to maximize the working efficiency of WT and PV, to make full use of clean energies to alleviate the power supply pressure of the grid. In Stage 2, from the economic and security views of the power system, minimizing the PVD concurrently suppressing the power fluctuation created by DGs is selected to be the optimization objective in this stage, and an incentive method for EVs coordinated scheduling is introduced, wherein TOU measure is utilized to provide an economic incentive to make EV owners charge during less congested regions and discharge in peak hours, combined with the variable power charge/discharge control of the BESS. Resultantly, the charge/discharge guidance of the connected EVs could be determined by figuring out the optimal peak and valley periods, and the real-time output plan of the BESS configuration associated with load shifting could also be solved. In Stage 3, based on the power load after peak shaking and valley filling function of EVs and BESS collaborative services, SVC is put into operation for reactive power compensation. The coordinated structure formulated as a comprehensive objective function aiming to concurrently minimize the active power loss and node voltage excursion of the DN is established, and the optimal daily output schedule of SVC capacitors can be dynamically obtained. The strategy not only maintains the balance of power supply and demand but favors the running costs and voltage quality, merging the high energy efficiency, economy, and security for system operation. Analysis of mathematical modeling of the hierarchical optimization problem is given&#x20;next.</p>
</sec>
<sec id="s3">
<title>3 Stage 1: MPPT Control of DGs</title>
<p>Renewable energy power generation such as WT and PV has the characteristics of non-pollution and sustainability but is intermittent and time-varying. The MPPT control method is employed to maximize the output efficiency and improve the response speed of the renewable power generation system, which is directly related to the effective utilization of wind and solar energy and the safe operation performances of the generation system. Specifically for the WT running characteristics, as provided in <xref ref-type="fig" rid="F3">Figure&#x20;3</xref>, under a constant blade pitch angle, the wind energy utilization coefficient changes with the tip speed ratio (<xref ref-type="bibr" rid="B13">Luo and Niu, 2016</xref>), and each wind velocity corresponds to an optimal mechanical speed of WT <italic>&#x3c9;</italic>
<sub>max</sub>, for example, <italic>S</italic> for wind velocity <italic>V</italic>
<sub>1</sub> and <italic>P</italic> for <italic>V</italic>
<sub>2</sub>, which makes the WT operate at MPP. By depicting all MPPs under different wind velocities, it forms a maximum power curve. Hence, it can be seen that satisfying the equation <italic>&#x3c9;</italic>
<sub>
<italic>m</italic>
</sub> &#x3d; <italic>&#x3c9;</italic>
<sub>max</sub> is the essential task for the MPPT control of WT. Similarly, for the PV generation system, which usually covers the PV panel, Boost circuit, inverter, and load, the duty cycle <italic>D</italic> of the switching device in the Boost circuit can be taken as the control parameter for MPPT to well reduce the complexity of the system. More precisely, the temperature changes mainly affect the PV output voltage, while the irradiation changes mainly affect the PV output current, and under a given temperature and sunlight intensity, the intersection point of the load-line with the PV characteristic determines the operating point (<xref ref-type="bibr" rid="B12">Koutroulis et&#x20;al., 2001</xref>). Consequently, the maximum power production is based on the load-line regulation by means of adjusting <italic>D</italic> to realize impedance matching, that is, the equivalent input resistance of the power converter equals the internal resistance of the PV cell (<xref ref-type="bibr" rid="B17">Podder et&#x20;al., 2019</xref>), to impel the working point to move toward the MPP correctly, and the process is declared in <xref ref-type="fig" rid="F4">Figure&#x20;4</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Characteristic curve of WT generation.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Characteristic curve of PV generation.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g004.tif"/>
</fig>
<p>Due to the connection through a gear box, the WT rotational speed <italic>&#x3c9;</italic>
<sub>
<italic>m</italic>
</sub> is directly proportional to the wind generator speed <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>; thus, the optimal <italic>&#x3c9;</italic>
<sub>max</sub> of WT corresponds to a reference <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>&#x2a; of the generator. According to <xref ref-type="fig" rid="F3">Figure&#x20;3</xref>, if the detected <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub> is equal to <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>&#x2a;, which indicates that the generation system operates at the MPP as expected, the excitation current value <italic>i</italic>
<sub>
<italic>c</italic>
</sub> of the wind power generator can be maintained. If the detected <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>&#x3c;<italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>&#x2a;, manifesting that the generation system works at Region <italic>&#x406;</italic>, and it is necessary to increase <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub> to approach the MPP; while <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>&#x3e;<italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>&#x2a; represents the operation state of the WT generation system at Region <italic>&#x2161;</italic>, and <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub> should be decreased to improve the conversion efficiency. Under such a non-equilibrium condition, regulation principles are established:<list list-type="simple">
<list-item>
<p>1) Region <italic>&#x406;</italic>: a positive speed variation &#x394;<italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub> illustrates that the working point is close to the MPP along the path 1 direction, and <italic>&#x3c9;</italic>
<sub>
<italic>m</italic>
</sub> increasing indicates that the output of the WT <italic>P</italic>
<sub>
<italic>m</italic>
</sub> may be greater than the generator <italic>P</italic>
<sub>
<italic>g</italic>
</sub>, and the excitation current <italic>i</italic>
<sub>
<italic>c</italic>
</sub> of the generator should be enhanced for power matching; contrastively, a negative &#x394;<italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub> signifies that the working point moves far away from the MPP along the path 3 direction, and the excitation current <italic>i</italic>
<sub>
<italic>c</italic>
</sub> needs to be rapidly reduced to satisfy <italic>P</italic>
<sub>
<italic>m</italic>
</sub> &#x3e; <italic>P</italic>
<sub>
<italic>g</italic>
</sub>, which drives the increase of <italic>&#x3c9;</italic>
<sub>
<italic>m</italic>
</sub> and <italic>&#x3c9;</italic>
<sub>
<italic>r</italic>
</sub>, and guarantee the correct moving direction.</p>
</list-item>
<list-item>
<p>2) Region <italic>&#x2161;</italic>: the WT generation system works in the high-speed area of the characteristic curves, and thus whether the working point moves along path 2 or 4, the generator output power <italic>P</italic>
<sub>
<italic>g</italic>
</sub> should be improved for optimal efficiency, that is, the excitation current <italic>i</italic>
<sub>
<italic>c</italic>
</sub> needs to be increased in this region.</p>
</list-item>
</list>
</p>
<p>Similar to WT analyzed above, the MPP tracking principle for PV generation also can be summarized as follows: if the PV generation system operates in Region <italic>&#x406;</italic>, the duty cycle <italic>D</italic> of Boost circuit should be increased, while <italic>D</italic> must be reduced as a response to the operating points in Region <italic>&#x2161;</italic>.</p>
</sec>
<sec id="s4">
<title>4 Stage 2: Coordinated Scheduling of Integrated EVs and BESS</title>
<p>The explosive growth of these DGs&#x2019; penetration brings power rush and uncertainty problems to the DN. Since the connected EVs and BESS can both function as power buffer for intermittent DGs and backup power supply, matching the bidirectional EVs&#x20;and BESS equipment for power smoothing is always considered as a simple and effective measure to solve the above problems.</p>
<sec id="s4-1">
<title>4.1 Coordinated Scheduling of Integrated EVs</title>
<p>In the presence of large-scale penetration of EVs into the building cluster, the load pressure, operating costs, and reliance on DN will be increased. Dispatch facilitates are coordinated to boost the utilization efficiency of EVs and improve the demand side flexibility. Due to the intrinsic advantages such as operational safety and energy density, the lithium-ion batteries of 25&#xa0;kWh are selected as the power battery for EVs (<xref ref-type="bibr" rid="B1">Barcellona et&#x20;al., 2019</xref>), and the charging power maintains at 2.5&#xa0;kW for a single vehicle. Concurrently, the safety threshold of the state of charge (SOC) is defined to be (10%, 90%); hence, it can be inferred that the charging/discharging could not last for more than 8&#xa0;h. This is supposing that the power consumption for an EV is 15&#xa0;kWh per 100&#xa0;km, and the theoretical endurance mileage can be calculated as 133&#xa0;km.</p>
<sec id="s4-1-1">
<title>4.1.1 Aggregated Charging of EVs</title>
<p>Aggregated charging of EVs could be explained as charging behavior only according to customers&#x2019; travel needs or living habits, without any guiding principle. The EV travel follows a probability density distribution represented in <xref ref-type="disp-formula" rid="e1">Eq. 1</xref> (<xref ref-type="bibr" rid="B18">Qian et&#x20;al., 2010</xref>) and depicted in <xref ref-type="fig" rid="F5">Figure&#x20;5</xref> as well, as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi mathvariant="normal">&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:mi>exp</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3bc;</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <italic>&#x3bc;</italic>
<sub>
<italic>d</italic>
</sub> &#x3d; 3.20, <italic>&#x3c3;</italic>
<sub>
<italic>d</italic>
</sub> &#x3d; 0.88, and <italic>x</italic> denotes the daily mileage of EV (mostly lies between 32 and 97&#xa0;km), so a charged 25-kWh battery can provide sufficient energy to meet daily driving requirements. Monte Carlo random sampling can be used to predict the charging power demand for a single EV during a day. In this study, a day is divided into 24 periods, with sampling of the states of EVs per hour. The charging expectation for a single vehicle is shown in <xref ref-type="fig" rid="F6">Figure&#x20;6</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Daily probability of driving distance.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Original charging requirements for a single EV.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g006.tif"/>
</fig>
</sec>
<sec id="s4-1-2">
<title>4.1.2 Coordinated Dispatch of EVs</title>
<p>The three-stage TOU price is designated to establish a positive guidance for the charging/discharging process of the EV owners. Depending on the load change of the power grid, the corresponding electricity price level of each period is confirmed, to achieve peak shaving and valley filling. In the modeling process for V2G service based on TOU, minimizing the PVD of the power load in the DN is chosen as the optimization objective, and some application background is specified as follows:<list list-type="simple">
<list-item>
<p>1) The battery power meets the driving demand of the EV owners, and there is no other power consumption behavior besides the normal travel.</p>
</list-item>
<list-item>
<p>2) User charge or discharge without exceeding the safe range of the SOC, and 80% of the total owners participate in the coordinated dispatch (<italic>d</italic>&#x20;&#x3d;&#x20;0.8).</p>
</list-item>
<list-item>
<p>3) <italic>m</italic>
<sub>
<italic>v</italic>
</sub>, <italic>m</italic>
<sub>
<italic>p</italic>
</sub>, and <italic>m</italic>
<sub>
<italic>f</italic>
</sub> are defined as the electricity price of the valley, peak, and normal period, respectively, and then the price model <italic>m</italic>
<sub>
<italic>t</italic>
</sub> can be described as follows:</p>
</list-item>
</list>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <italic>t</italic>
<sub>1</sub>&#x2013;<italic>t</italic>
<sub>4</sub> represents the start and end moment of the valley period and the peak period, respectively.<list list-type="simple">
<list-item>
<p>4) Before charging/discharging, the owners can query the current battery status of the EVs and independently select the charging/discharging time. The parameters <italic>t</italic>
<sub>
<italic>sc</italic>
</sub>, <italic>t</italic>
<sub>
<italic>c</italic>
</sub>, <italic>t</italic>
<sub>
<italic>sd</italic>
</sub>, and <italic>t</italic>
<sub>
<italic>d</italic>
</sub> are defined as follows:</p>
</list-item>
<list-item>
<p>
<italic>t</italic>
<sub>
<italic>sc</italic>
</sub>, <italic>t</italic>
<sub>
<italic>sd</italic>
</sub> -the start moment for users to charge/discharge; <italic>t</italic>
<sub>
<italic>c</italic>
</sub>, <italic>t</italic>
<sub>
<italic>d</italic>
</sub> -duration of the charging/discharge process.</p>
</list-item>
</list>
</p>
<p>EV users participating in the V2G dispatch can severally choose <italic>t</italic>
<sub>
<italic>sc</italic>
</sub> or <italic>t</italic>
<sub>
<italic>sd</italic>
</sub> according to <xref ref-type="disp-formula" rid="e3">Eqs 3</xref>, <xref ref-type="disp-formula" rid="e4">4</xref>:<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <italic>r</italic>
<sub>
<italic>c</italic>
</sub> and <italic>r</italic>
<sub>
<italic>d</italic>
</sub> are random numbers &#x2282; [0, 1] interval. It can be known that for the coordinated EVs dispatch model, <italic>t</italic>
<sub>1</sub>&#x2013;<italic>t</italic>
<sub>4</sub> determine the peak and valley areas, and hence are key factors for optimal V2G service quality, which lays an important foundation for the subsequent multi-objective DN scheduling.</p>
</sec>
</sec>
<sec id="s4-2">
<title>4.2 Optimal Planning of BESS</title>
<p>The BESS system is commonly connected to the power grid together with the DGs, for smoothing the power fluctuation of renewable power generation. Similar to EVs, BESS can also operate in both charge and discharge states, and it can run through all parts of the power system. Recent research on the application of BESS mainly focuses on two issues: a dispatch from the perspective of BESS cost and load shifting to directly achieve peak shaving and valley filling. Through BESS modeling, the best charge/discharge plan can be solved. BESS could work in both of fixed and variable power mode, considering that although the control process of the fixed power mode is relatively simple, the output of BESS may exceed the power demand; hence, the variable power control mode is selected in this study to guide the charging or discharging behaviors of the BESS system and satisfy the real-time power demand of residents dynamically and, finally, acquire a daily load curve with minimized PVD. Suppose that the initial SOC of the BESS is 10%, <italic>p</italic>
<sub>
<italic>c_bess</italic>
</sub>(<italic>t</italic>) and <italic>p</italic>
<sub>
<italic>d_bess</italic>
</sub>(<italic>t</italic>) denote the charge and discharge power demand of BESS at the <italic>t</italic>th moment, and <italic>&#x3b7;</italic>
<sub>
<italic>c</italic>
</sub> and <italic>&#x3b7;</italic>
<sub>
<italic>d</italic>
</sub> demonstrate the charge and discharge efficiency factors, respectively, Under the BESS power balance principle in <xref ref-type="disp-formula" rid="e5">Eq. 5</xref>, the variable power charge/discharge control is designed as the following steps and represented in <xref ref-type="fig" rid="F7">Figure&#x20;7</xref>:<disp-formula id="e5">
<mml:math id="m5">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>N</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>N</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<list list-type="simple">
<list-item>
<p>1) Collect temporary data of basic power demand, PV and WT power generation with MPPT control, and coordinated charging/discharging states of EVs in the regional&#x20;DN.</p>
</list-item>
<list-item>
<p>2) According to the gathered information in (1), determine the peak shaving line of the BESS, and calculate the capacity required for daily load peak shaving, which concurrently equals to the capacity to be absorbed for the BESS during off-peak&#x20;hours.</p>
</list-item>
<list-item>
<p>3) Set a horizontal line at the minimum load as the valley filling line and gradually move it upward with a certain step size until the output power balance is satisfied.</p>
</list-item>
<list-item>
<p>4) Confirm the optimal peak shaving and valley filling line, thus obtaining the solution of the daily output plan of the BESS and completing the BESS scheduling process.</p>
</list-item>
</list>
</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Variable power charge/discharge control flow of the BESS.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g007.tif"/>
</fig>
</sec>
</sec>
<sec id="s5">
<title>5 Stage 3: Multi-Objective Reactive Compensation of DN</title>
<p>The guiding ideology for the whole DN optimization can be concluded as follows: on the premise that the parameters including branches and loads are known and various constraints are satisfied, dynamically adjust the output of the SVC capacitor banks to optimize the comprehensive indices of active power loss and voltage excursion. Finally, the economic and safe operation of the multi-source DN can be achieved.</p>
<sec id="s5-1">
<title>5.1 Objective Function</title>
<p>The optimization objective function comprises two important issues: active power loss and node voltage excursion of the DN.<list list-type="simple">
<list-item>
<p>1) Active power&#x20;loss</p>
</list-item>
</list>
</p>
<p>The first objective is the minimization of the total active power loss of the DN, which can be mathematically modeled using <xref ref-type="disp-formula" rid="e6">Eq. 6</xref>:<disp-formula id="e6">
<mml:math id="m6">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where <italic>P</italic>
<sup>
<italic>&#x2032;</italic>
</sup>
<sub>
<italic>Lij</italic>
</sub> and <italic>P</italic>
<sub>
<italic>Lij</italic>
</sub>, denote the active power loss of the branch (<italic>i</italic>, <italic>j</italic>) before and after optimization, respectively.<list list-type="simple">
<list-item>
<p>2) Voltage excursion</p>
</list-item>
</list>
</p>
<p>Node voltage is an important indicator reflecting the security and service quality. To avoid all the voltages moving toward their maximum limits after optimization, the deviation of voltage from the rated value is chosen as an objective function as shown below:<disp-formula id="e7">
<mml:math id="m7">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mtext>N</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mtext>N</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>where <italic>U</italic>
<sup>
<italic>&#x2019;</italic>
</sup>
<sub>
<italic>j</italic>
</sub>, <italic>U</italic>
<sub>
<italic>j</italic>
</sub>, and <italic>U</italic>
<sub>
<italic>j</italic>N</sub>, respectively, represent the actual voltage of the node <italic>j</italic> before and after optimization, and the rated voltage of the node&#x20;<italic>j</italic>.</p>
<p>Combining the above two indicators, the objective function can be expressed as follows:<disp-formula id="e8">
<mml:math id="m8">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mi>F</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mtext>N</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msubsup>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mtext>N</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <italic>&#x3c9;</italic>
<sub>1</sub> and <italic>&#x3c9;</italic>
<sub>2</sub> are the weight coefficients for the two optimization objectives (<italic>&#x3c9;</italic>
<sub>1</sub> &#x3d; <italic>&#x3c9;</italic>
<sub>2</sub> &#x3d; 0.5). The penalty function is used to deal with the node voltage out-of-limit problem and <italic>U</italic>
<sub>
<italic>j</italic>min</sub> and <italic>U</italic>
<sub>
<italic>j</italic>max</sub> are the minimum and maximum voltage values of the node <italic>j</italic>, while <italic>&#x3bc;</italic> indicates the penalty factor (<italic>&#x3bc;</italic> &#x3d; 1000), and &#x394;<italic>U</italic>
<sub>
<italic>j</italic>
</sub> is defined as follows:<disp-formula id="e9">
<mml:math id="m9">
<mml:mrow>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
</sec>
<sec id="s5-2">
<title>5.2 Constraints</title>
<p>The basic constraints of DN optimization mainly include equality constraints (power flow) and inequality constraints as follows:<list list-type="simple">
<list-item>
<p>1) Equality constraints:</p>
</list-item>
</list>
<disp-formula id="e10">
<mml:math id="m10">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m11">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf2">
<mml:math id="m12">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, respectively, represent the output power of WT and PV at the <italic>t</italic>th moment; besides, <inline-formula id="inf3">
<mml:math id="m13">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf4">
<mml:math id="m14">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> denote the real-time discharge demands of dispersed EVs and BESS devices, which is positive in the discharging state and negative in the charging state. <inline-formula id="inf5">
<mml:math id="m15">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> expresses the actual purchased power from the grid if the power supply depending on the distributed energy equipment is less than user demand, and <inline-formula id="inf6">
<mml:math id="m16">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>B</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the basic load of the DN at the <italic>t</italic>th moment.<list list-type="simple">
<list-item>
<p>2) Inequality constraints:</p>
<list list-type="simple">
<list-item>
<p>1) SOC constraint of EVs and BESS:</p>
</list-item>
</list>
</list-item>
</list>
<disp-formula id="e11">
<mml:math id="m17">
<mml:mrow>
<mml:mn>10</mml:mn>
<mml:mo>%</mml:mo>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>90</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m18">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>p</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>W</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m19">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>p</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>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m20">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>k</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>where <inline-formula id="inf7">
<mml:math id="m21">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</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>W</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf8">
<mml:math id="m22">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</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>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> illustrate the maximum output capacities of the WT and PV generation units; <italic>C</italic>
<sub>
<italic>t</italic>
</sub> and <italic>C</italic>
<sub>max</sub> denote the output compensation capacity of the installed SVC capacitor bank at the <italic>t</italic>th moment, and its maximum capacity, respectively; <italic>C</italic>
<sub>
<italic>Nt</italic>
</sub> is the capacity of a single capacitor included in the capacitor bank at the <italic>t</italic>th moment; <italic>k</italic> is the number of single capacitors on each compensation&#x20;node.</p>
</sec>
<sec id="s5-3">
<title>5.3 Algorithm</title>
<p>The purpose of multi-objective optimization is to determine a set of Pareto solutions, which take each optimization objective into account. GA is a global probability searching tool, including selection, crossover, and mutation operators to retain the high-quality individuals in the population. Herein, in order to improve the convergence speed of the traditional GA, EPS policy is introduced, which is featured with that the excellent individuals owing the best fitness are gathered to form a sub population and directly copied to the next generation without participating in the crossover and mutation steps. It could protect the elite individuals from the disturbance of crossover and mutation operations in the traditional GA, thus improving the stability and convergence of the algorithm. Moreover, the crossover probability <italic>P</italic>
<sub>
<italic>c</italic>
</sub> and mutation probability <italic>P</italic>
<sub>
<italic>m</italic>
</sub> are adaptively adjusted following the fitness values, which are perfected as <xref ref-type="disp-formula" rid="e15">Eqs 15</xref>, <xref ref-type="disp-formula" rid="e16">16</xref> (<xref ref-type="bibr" rid="B9">Huang et&#x20;al., 2020</xref>):<disp-formula id="e15">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>F</mml:mi>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mtext>avg</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x2265;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mtext>avg</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mtext>avg</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>F</mml:mi>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mtext>avg</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x2265;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mtext>avg</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtable>
<mml:mtr>
<mml:mtd>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>F</mml:mi>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x2032;</mml:mo>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mtext>avg</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
</p>
<p>where <italic>F</italic>
<sub>avg</sub>, <italic>F</italic>
<sub>max</sub>, <italic>F</italic>
<sup>
<italic>&#x2032;</italic>
</sup>, and <italic>F</italic>
<sup>
<italic>&#x2033;</italic>
</sup>, demonstrate the average and maximum fitness value of individuals, the better fitness of the two crossover individuals, and the mutation individuals, respectively. <italic>k</italic>
<sub>1</sub>&#x2013;<italic>k</italic>
<sub>4</sub> are set to be 0.5, 0.9, 0.02, and 0.05. By improving <italic>P</italic>
<sub>
<italic>c</italic>
</sub> and <italic>P</italic>
<sub>
<italic>m</italic>
</sub>, individuals with higher fitness value than the average correspond to lower crossover and mutation probability, while individuals below the average correspond to higher <italic>P</italic>
<sub>
<italic>c</italic>
</sub> and <italic>P</italic>
<sub>
<italic>m</italic>
</sub>, which helps to duplicate the good individuals and eliminate the bad solutions. <xref ref-type="fig" rid="F8">Figure&#x20;8</xref> shows the flowchart of the algorithm. After the load flow calculation of the initialized population, 20% of the individuals with the best compromise effects between the two objective functions are reproduced directly to the next generation, while the rest are selected using the roulette wheel method and generated according to the crossover and mutation operations. The global optimal solution is acquired through an iterative calculation, with the group size 30 and iteration&#x20;500.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Flowchart of the optimizing process <italic>via</italic> the developed GA with EPS.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g008.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="results|discussion" id="s6">
<title>6 Results and Discussion</title>
<p>The standard IEEE 33-node power DN system shown in <xref ref-type="fig" rid="F9">Figure&#x20;9</xref> is taken for verification, and its load information is detailed in Reference (<xref ref-type="bibr" rid="B9">Huang et&#x20;al., 2020</xref>). The basic parameters of the network for the per-unit system are 10 MVA and 12.66&#xa0;kV. The proposed hierarchical optimization strategy is implemented on the MATLAB platform, and the research cycle is 24&#xa0;h. The IEEE 33-node DN mainly contains: (1) residential building load (original power load without DGs and EVs integration), which is assumed to evenly distribute on four nodes: 3, 10, 18, and 32. Considering the time-varying characteristics of the daily load, a series of residential building load data for a certain day are collected from the study by <xref ref-type="bibr" rid="B9">Huang et&#x20;al. (2020)</xref>, with the PVD of 1763.8&#xa0;kW, as shown in <xref ref-type="fig" rid="F10">Figure&#x20;10</xref>. (2) Renewable energy power generation. DG1 (WT) and DG2 (PV) are installed on nodes 2 and 5, respectively, with the optimal output power under MPPT control. For the day studied, meteorological data involving wind velocity, temperature, and sunlight intensity (<xref ref-type="bibr" rid="B8">Hosseinalizadeh et&#x20;al., 2016</xref>) associated with renewable energy power generation are depicted in <xref ref-type="fig" rid="F11">Figures 11</xref>, <xref ref-type="fig" rid="F12">12</xref>, with sampling per hour in view of the randomness and uncertainty of climate. (3) BESS, positioned at nodes 3 and 10, with individually rated power of 100&#xa0;kW and capacity of 500&#xa0;kWh. According to (<xref ref-type="bibr" rid="B24">Yao et&#x20;al., 2017</xref>), the charge/discharge efficiencies are selected as <italic>&#x3b7;</italic>
<sub>
<italic>c</italic>
</sub> &#x3d; 0.95, <italic>&#x3b7;</italic>
<sub>
<italic>d</italic>
</sub> &#x3d; 0.9, respectively. (4) SVC, placed on the weakest two nodes for compensation, and each capacitor bank is supposed to provide the maximum capacity <italic>C</italic>
<sub>
<italic>j</italic>max</sub> &#x3d; 1 Mvar. (5) EV charging loads, which can also be deemed as power suppliers under the V2G mode. It is assumed that 500&#xa0;EVs are accessed to the building cluster, with the charge/discharge power of 2.5&#xa0;kW. Furthermore, the convergence accuracy of power flow calculation is 10<sup>&#x2013;4</sup>.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>IEEE 33-node system.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Daily original power&#x20;load.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Daily wind velocity&#x20;data.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g011.tif"/>
</fig>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Daily temperature and irradiation intensity.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g012.tif"/>
</fig>
<sec id="s6-1">
<title>6.1 Discussion of Simulation Results in Stage 1</title>
<p>Based on the collected daily wind velocity, temperature, and irradiation intensity distribution data in <xref ref-type="fig" rid="F11">Figures 11</xref>, <xref ref-type="fig" rid="F12">12</xref>, the output of WT and PV under MPPT control can be obtained, as shown in <xref ref-type="fig" rid="F13">Figure&#x20;13</xref>. Despite the randomness of the climate, the DGs achieve the maximum output power within a short self-regulation time, concurrently with small oscillation and thus desired power equality. Concretely, from the original building load curve, it is obvious that the power system is most prone to blackout at about 19:00, with the peak load of 4,370.1&#xa0;kW. The MPPT control impels the WT and PV system efficiently operates with 631.82 and 23.28 kW, making the power demand optimized to 3,715&#xa0;kW under Scenario 1 (decreased by 14.99%), which is conducive to alleviating the power supply burden of the DN (<xref ref-type="fig" rid="F14">Figure&#x20;14</xref>). Besides, among all the MPPs, the PV system has the maximum daily output power of 610.35&#xa0;kW at 13:00, due to the high temperature and sunlight at noon, while the maximum power of WT (747.80&#xa0;kW) occurs at 2:00. These moments are not identified as peak power demand hours, and hence the high-efficiency operation of DGs may not be beneficial for valley filling of the load curve and affect the load rate as well as peak regulation management of the power system.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Output of DGs under MPPT control.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g013.tif"/>
</fig>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Load curve after MPPT control of DGs.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g014.tif"/>
</fig>
</sec>
<sec id="s6-2">
<title>6.2 Discussion of Simulation Results in Stage 2</title>
<p>According to the operation modes of EVs analyzed in <xref ref-type="sec" rid="s4">Section 4</xref>, the simulation research is carried out for three scenarios: Scenario 1 represents the power load with DGs compensation, whereas without integrated EVs or BESS; Scenario 2 is the aggregated charging mode of the 500&#xa0;EVs, which are connected to charge station without guiding measures; while 80% of the EVs are turned on to be coordinated V2G operation mode in Scenario 3. Based on Stage 1, fully considering the travel demands of the EV owners, the load curves with integrated EVs under Scenario 2 and 3 are obtained, respectively, as exhibited in <xref ref-type="fig" rid="F15">Figure&#x20;15</xref>.</p>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Load curves under different scenarios.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g015.tif"/>
</fig>
<p>It can be observed that the original power load between 17:00 and 20:00 is intensified by aggregated EVs charging behavior, which puts forward higher requirements for the power supply of the DN. Concretely, for the load curve under Scenario 1, the peak load of 3,715&#xa0;kW occurs at about 19:00, which would be further increased to 4,166.9&#xa0;kW under Scenario 2, meaning that the overloading even power outage is most likely to take place at 19:00. Contrastively, under Scenario 3, the introduction of TOU price fully makes use of the energy reserve capacity and bidirectional power flow of EVs and realizes the power load shifting by encouraging owners to charge in valley hours and discharge in peak hours. The optimal peak and valley periods are, solved as (15:00&#x2013;23:00) and (23:00&#x2013;7:00), respectively, and the peak load at 19:00 is optimized to 3,608.3&#xa0;kW (decreased by 13.41% compared to that of Scenario 2); concurrently, the PVD is evaluated to be 1,087.1&#xa0;kW, decreased by 41.75% in contrast to the PVD of 1866.3&#xa0;kW under Scenario 2. From the above analysis, it can be noted that residential power demands during peak hours can be effectively reduced, and the load curve is remarkably smoothed after peak shaving and valley filling, conducive to the subsequent active power loss and node voltage excursion optimization.</p>
<p>Through the power flow analysis, voltages for all of the nodes at 19:00 are also calculated (exhibited in <xref ref-type="fig" rid="F16">Figure&#x20;16</xref>). Supposing that a voltage dip lower than 0.95 p.u. or a voltage higher than 1.05 p.u. is deemed as voltage out-of-limit, then it is discovered that the integration of the building load increases the power supply burden of the DN system, since the voltage drop problems are more serious. Only 11 node voltages are better than the lower bound, and the worst two occur on nodes 17 and 32 of 0.8378 and 0.8294 p.u.. DG1 and DG2 separately connected to nodes 2 and 5 can further support the system voltage. After the MPPT control for DG input (Scenario 1), the lowest value is still on the nodes of 17 and 32, which are raised to 0.8535 and 0.8472 p.u., and the relative voltage excursion of each node adds up to 2.9347 p.u., decreased by 10.54% compared to 3.2806 p.u. of the original building load. Next, with the integration of the EVs, it is discovered that the system voltage quality problem is intensified under the EV aggregated charging (Scenario 2), which causes the voltages of the two weakest nodes to seriously drop to 0.8451 p.u. and 0.8373 p.u. Contrastively, the system voltage is overall heightened on the coordinated regulation mode (Scenario 3), and the voltages of nodes 17 and 32 are raised up to 0.8555 p.u. and 0.8495 p.u., slightly better than the DN system operates in Scenario 1. Moreover, <xref ref-type="table" rid="T1">Table&#x20;1</xref> illustrates the statistical results of the key indices involving the relative voltage excursion and power loss for the DN system under different scenarios. It is obvious that first, based on the flexible power supply of the DGs, MPPT control improves the utilization efficiency of clean energy and power capacity. Furthermore, due to the introduction of the TOU mechanism, the coordinated management of flexible loads inside the buildings not only boosts the utilization of EV but is conducive to the high efficiency and reliability of the DN operation. These follow the concept of energy conservation and emission reduction of the modern power system.</p>
<fig id="F16" position="float">
<label>FIGURE 16</label>
<caption>
<p>Voltage curves under different scenarios at 19:00.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g016.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Indices under different scenarios for the load peak at 19:00.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Regulation mode</th>
<th align="center">Relative voltage excursion (p.u.)</th>
<th align="center">Power loss (kW)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Without building load</td>
<td align="char" char=".">1.7007</td>
<td align="char" char=".">202.65</td>
</tr>
<tr>
<td align="left">Basic load</td>
<td align="char" char=".">3.2806</td>
<td align="char" char=".">797.63</td>
</tr>
<tr>
<td align="left">Scenario 1</td>
<td align="char" char=".">2.9347</td>
<td align="char" char=".">627.38</td>
</tr>
<tr>
<td align="left">Scenario 2</td>
<td align="char" char=".">3.1124</td>
<td align="char" char=".">712.45</td>
</tr>
<tr>
<td align="left">Scenario 3</td>
<td align="char" char=".">2.8934</td>
<td align="char" char=".">608.45</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Furthermore, to elaborate the regulation results for the next control stage, the other two scenarios are defined.</p>
<p>Scenario 4: variable power dispatch of BESS, based on Scenario&#x20;3.</p>
<p>Scenario 5: SVC, capacitor banks put into operation based on Scenario&#x20;4.</p>
<p>Considering the impact of the intermittent generation from the integrated DGs on the power grid, the BESS system is commonly necessary for power load smoothing. Since the BESS has the attributes of source and load, it could further flatten the load curve by reasonably transferring power demand from peak to valley. Based on the variable power charge/discharge dispatch control of the BESS in <xref ref-type="fig" rid="F7">Figure&#x20;7</xref>, it could be calculated that the peak shaving line (boundary 1) of 3,424.99&#xa0;kW and valley filling line (boundary 2) of 2,793.49&#xa0;kW make the PVD of the load curve reduce from 1,087.1&#xa0;kW (Scenario 3) to 703.80&#xa0;kW (Scenario 4, decreased about 35.26%). The peak shaving and valley filling effects and the real-time output of the BESS are shown in <xref ref-type="fig" rid="F17">Figures 17</xref>, <xref ref-type="fig" rid="F18">18</xref>, respectively. The optimal scheduling of the BESS further consolidates the load shifting function, making the load curve more flat and centralized and conducive to power saving, to ensure more economic and safe operation of the power grid on the basis of the MPPT control of DGs and coordinated dispatch of&#x20;EVs.</p>
<fig id="F17" position="float">
<label>FIGURE 17</label>
<caption>
<p>Load curves under different scenarios of EVs.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g017.tif"/>
</fig>
<fig id="F18" position="float">
<label>FIGURE 18</label>
<caption>
<p>Daily dispatch plan of the BESS.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g018.tif"/>
</fig>
</sec>
<sec id="s6-3">
<title>6.3 Discussion of Simulation Results in Stage 3</title>
<p>Ultimately, by using the improved GA with the EPS method analyzed in 5.3, the optimal compensation capacities of the SVC capacitor banks on nodes 17 and 32 are solved as 0.8513 and 0.9130 MVar, respectively. The risk of voltage out-of-limit is effectively reduced after reactive compensation, which can be testified by the final node voltage of 17 and 32 (0.9224 and 0.8961 p.u.), a significant improvement compared to values of Scenario 4 (0.8603 and 0.8520 p.u.). Based on the coordinated dispatch of EVs and BESS devices, <xref ref-type="fig" rid="F19">Figure&#x20;19</xref> declares the node voltages of the DN system at 19:00 before and after SVC capacitor compensation, emphasizing its excellent effect on voltage fluctuation suppression. After the proposed hierarchical optimization, the total voltage excursion is optimized to 2.0605 p.u. (decreased by 37.19% compared to that of the basic building load). More concretely, the daily output dispatch plan of the SVC can be dynamically solved out, as described in <xref ref-type="fig" rid="F20">Figure&#x20;20</xref>. It can be observed that with the daily power load changing, the output capacity of the SVC is adjusted dynamically, to concurrently reduce the node voltage fluctuation and active power&#x20;loss.</p>
<fig id="F19" position="float">
<label>FIGURE 19</label>
<caption>
<p>Voltage curves under different scenarios at 19:00.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g019.tif"/>
</fig>
<fig id="F20" position="float">
<label>FIGURE 20</label>
<caption>
<p>Daily dispatch plan of the SVC.</p>
</caption>
<graphic xlink:href="fenrg-10-869844-g020.tif"/>
</fig>
<p>In addition, as shown in <xref ref-type="table" rid="T2">Table&#x20;2</xref>, when the SVC devices are not installed (Scenario 4), the active power loss of the system is decreased to 575.99 kW; which will be further reduced to 514.90&#xa0;kW when the reactive compensation is available (Scenario 5). Compared with the power consumption of 797.63&#xa0;kW for the original building load, the reduction under Scenario 1&#x2013;5 achieves 21.34, 10.68, 23.72, 27.79, and 35.45%, respectively. Therefore, the MPPT control of the penetrated DGs, combined with the reasonable dispatch of EVs and BESS, could weaken the burden of the DN for power supply. Concurrently, the multi-objective reactive optimization approach also ensures the significant improvements of system power loss and node voltage quality, proving that the proposed hierarchical optimization strategy in this study could obtain outstanding performances. It is noteworthy that the response degree of total EV owners <italic>d</italic> has a crucial impact on the optimal effects. A higher <italic>d</italic> means more owners participating in the V2G management and the more obvious effect of minimizing the PVD resultantly, which makes for enhancing the initiative of users by economic interests and, concurrently, the economy and security operation of the DN system.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Indices before and after compensation at 19:00.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Compensation</th>
<th align="center">Relative voltage excursion (p.u.)</th>
<th align="center">Power loss (kW)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Scenario 4</td>
<td align="char" char=".">2.8133</td>
<td align="char" char=".">575.99</td>
</tr>
<tr>
<td align="left">Scenario 5</td>
<td align="char" char=".">2.0605</td>
<td align="char" char=".">514.90</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s7">
<title>7 Conclusion</title>
<p>As a demand-oriented energy supply system, a smart DN with multi-source coupling is a major strategic measure for domestic energy development and the trend of the global power system. However, owing to the hybrid access of user resources, natural resources, and time-space dispersed resources making the power supply diversified and unstable, the DN system faces great changes in architecture, which brings great challenges to the upper collaborative dispatch and optimization. Consequently, establishing on the fundamental principle of multi-source integration and DN optimization control, the hierarchical multi-source coordinated regulation strategy considering the penetration of DGs, EVs, and BESS is proposed in this study. Combining the efficient utilization of DGs, load shifting function of EVs and BESS, and the dynamic reactive compensation of SVC capacitor banks, the optimization strategy could effectively balance the tradeoff between active power loss and node voltage fluctuation, and its feasibility and effectiveness are illustrated through MATLAB simulation results. First, based on MPPT control, the self-optimizing processes of WT and PV output power are implemented for improving the power production efficiency of the DGs. Next, according to the dynamic load change, the TOU price mechanism is introduced to improve the initiative and utilization of bidirectional EVs and BESS by guiding the charging/discharging plans of the owners and achieving the peak shaving and valley filling function, which is significant for alleviating the power supply pressure of the DN. Finally, both total active power loss and node voltage fluctuation are selected as the optimization objectives, under the basic constraints, and the improved GA with EPS method is employed to determine the daily optimal SVC capacitor bank compensation plan for the multi-objective economic-secure scheduling problem. After the above hierarchical dispatch steps, the entire power loss and node voltage out-of-limit risk in the DN can be improved simultaneously, which promotes the economic and reliable operation of the regional DN. The research results in this study have good extensibility and may provide certain reference significance for exploring effective flexible load modeling and synergetic management technology of &#x201c;source-grid-storage-load&#x201d; system in the subsequent construction of smart DNs. Nevertheless, the application would be subject to some limitations. First, if other objective functions are also paid attention to, it is better to adopt efficient multi-objective solvers instead of a single-objective optimizer by regulating the weight coefficients &#x3c9;<sub>1</sub> and &#x3c9;<sub>2</sub>. Moreover, influences of parameter uncertainties, such as DGs and power load are not taken into account in this study. Finally, if different types of distributed energy equipment are to be surveyed, mathematical models should be modified according to the actual running characteristics.</p>
</sec>
</body>
<back>
<sec 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 authors.</p>
</sec>
<sec id="s9">
<title>Author Contributions</title>
<p>Conceptualization and methodology: ZH; software: PM and MW; modeling: ZH, PM, and MW; validation: ZH; writing&#x2014;original draft preparation: ZH; writing&#x2014;review and editing: BF and MZ; supervision: MZ; project administration:&#x20;BF.</p>
</sec>
<sec id="s10">
<title>Funding</title>
<p>This work was supported by the Natural Science Foundation of Hunan Province (2017JJ5044), the Scientific Research Foundation of Hunan Provincial Education Department (20C0613, 21C0428), and the Hunan High-Tech Industry Science and Technology Innovation Leading Foundation of China (2022GK4055).</p>
</sec>
<sec sec-type="COI-statement" id="s11">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="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>
<ack>
<p>Most sincere thanks are due to Miss Hong and Dr. Wang for the kind care during my injury and the great encouragements to me to finish this&#x20;work.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barcellona</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Piegari</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Villa</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Passive Hybrid Energy Storage System for Electric Vehicles at Very Low Temperatures</article-title>. <source>J.&#x20;Energ. Storage</source> <volume>25</volume>, <fpage>100833</fpage>. <pub-id pub-id-type="doi">10.1016/j.est.2019.100833</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brockway</surname>
<given-names>P. E.</given-names>
</name>
<name>
<surname>Owen</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Brand-Correa</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Hardt</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Estimation of Global Final-Stage Energy-Return-On-Investment for Fossil Fuels with Comparison to Renewable Energy Sources</article-title>. <source>Nat. Energ.</source> <volume>4</volume> (<issue>7</issue>), <fpage>612</fpage>&#x2013;<lpage>621</lpage>. <pub-id pub-id-type="doi">10.1038/s41560-019-0425-z</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Distribution System Reconfiguration Considering Distributed Generators and Plug-In Electric Vehicles [J]</article-title>. <source>Power Syst. Tech.</source> <volume>37</volume> (<issue>1</issue>), <fpage>82</fpage>&#x2013;<lpage>88</lpage>. </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Pu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Distributionally Robust Optimization Model of Active Distribution Network Considering Uncertainties of Source and Load</article-title>. <source>J.&#x20;Mod. Power Syst. Clean. Energ.</source> <volume>7</volume> (<issue>6</issue>), <fpage>1585</fpage>&#x2013;<lpage>1595</lpage>. <pub-id pub-id-type="doi">10.1007/s40565-019-0558-x</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elbatawy</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Morsi</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Integration of Prosumers with Battery Storage and Electric Vehicles via Transactive Energy</article-title>. <source>IEEE Trans. Power Deliv.</source> <volume>37</volume> (<issue>1</issue>), <fpage>383</fpage>&#x2013;<lpage>394</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrd.2021.3060922</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Multi-Objective Comprehensive Charging/Discharging Scheduling Strategy for Electric Vehicles Based on the Improved Particle Swarm Optimization Algorithm</article-title>. <source>Front. Energ. Res.</source> <volume>9</volume>, <fpage>811964</fpage>. <pub-id pub-id-type="doi">10.3389/fenrg.2021.811964</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gholami Farkoush</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wadood</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Khurshaid</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>C.-H.</given-names>
</name>
<name>
<surname>Rhee</surname>
<given-names>S.-B.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Minimizing Static VAR Compensator Capacitor Size by Using SMC and ASRFC Controllers in Smart Grid with Connected EV Charger</article-title>. <source>Int. J.&#x20;Electr. Power Energ. Syst.</source> <volume>107</volume>, <fpage>656</fpage>&#x2013;<lpage>667</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2018.12.029</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hosseinalizadeh</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Shakouri G</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Amalnick</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Taghipour</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Economic Sizing of a Hybrid (PV-WT-FC) Renewable Energy System (HRES) for Stand-Alone Usages by an Optimization-Simulation Model: Case Study of Iran</article-title>. <source>Renew. Sust. Energ. Rev.</source> <volume>54</volume>, <fpage>139</fpage>&#x2013;<lpage>150</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2015.09.046</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Multi-objective Optimization Strategy for Distribution Network Considering V2G-Enabled Electric Vehicles in Building Integrated Energy System [J]</article-title>. <source>Prot. Control. Mod. Power Syst.</source> <volume>5</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1186/s41601-020-0154-0</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Islam</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Hossain</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Islam</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Multiobjective Optimization Technique for Mitigating Unbalance and Improving Voltage Considering Higher Penetration of Electric Vehicles and Distributed Generation</article-title>. <source>IEEE Syst. J.</source> <volume>14</volume> (<issue>3</issue>), <fpage>3676</fpage>&#x2013;<lpage>3686</lpage>. <pub-id pub-id-type="doi">10.1109/jsyst.2020.2967752</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A Three-Stage Multi-Timescale Framework for Online Dispatch in a Microgrid with EVs and Renewable Energy [J]</article-title>. <source>IEEE Trans. Trans Electri</source>. <volume>1</volume> <fpage>1</fpage>. <pub-id pub-id-type="doi">10.1109/tte.2021.3098687</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Koutroulis</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Kalaitzakis</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Voulgaris</surname>
<given-names>N. C.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Development of a Microcontroller-Based, Photovoltaic Maximum Power point Tracking Control System</article-title>. <source>IEEE Trans. Power Electron.</source> <volume>16</volume> (<issue>1</issue>), <fpage>46</fpage>&#x2013;<lpage>54</lpage>. <pub-id pub-id-type="doi">10.1109/63.903988</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Niu</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A Novel Contra-rotating Power Split Transmission System for Wind Power Generation and its Dual MPPT Control Strategy [J]</article-title>. <source>IEEE Trans. Power Elect.</source> <volume>32</volume> (<issue>9</issue>), <fpage>6924</fpage>&#x2013;<lpage>6935</lpage>. </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Wan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Charging Scheduling Strategy for Different Electric Vehicles with Optimization for Convenience of Drivers, Performance of Transport System and Distribution Network</article-title>. <source>Energy</source> <volume>194</volume>, <fpage>116807</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2019.116807</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mahmoudi</surname>
<given-names>S. K.</given-names>
</name>
<name>
<surname>Golshan</surname>
<given-names>M. E. H.</given-names>
</name>
<name>
<surname>Zamani</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Coordinated Voltage Control Scheme for Transmission System Considering Objectives and Constraints of Network and Control Devices</article-title>. <source>Electric Power Syst. Res.</source> <volume>192</volume>, <fpage>106908</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2020.106908</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peidong</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yanli</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>jin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yonghong</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Lisheng</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xinrong</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Opportunities and Challenges for Renewable Energy Policy in China</article-title>. <source>Renew. Sust. Energ. Rev.</source> <volume>13</volume> (<issue>2</issue>), <fpage>439</fpage>&#x2013;<lpage>449</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2007.11.005</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Podder</surname>
<given-names>A. K.</given-names>
</name>
<name>
<surname>Roy</surname>
<given-names>N. K.</given-names>
</name>
<name>
<surname>Pota</surname>
<given-names>H. R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>MPPT Methods for Solar PV Systems: a Critical Review Based on Tracking Nature</article-title>. <source>IET Renew. Power Generation</source> <volume>13</volume> (<issue>10</issue>), <fpage>1615</fpage>&#x2013;<lpage>1632</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2018.5946</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qian</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Allan</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Modeling of Load Demand Due to EV Battery Charging in Distribution Systems [J]</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>26</volume> (<issue>2</issue>), <fpage>802</fpage>&#x2013;<lpage>810</lpage>. </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Siahroodi</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Mojallali</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Mohtavipour</surname>
<given-names>S. s.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Scenario-based Stochastic Framework for Harmonic Power Markets Using Plug-In Electric Vehicles</article-title>. <source>J.&#x20;Energ. Storage</source> <volume>35</volume>, <fpage>102290</fpage>. <pub-id pub-id-type="doi">10.1016/j.est.2021.102290</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tan</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Multi-objective Energy Management of Multiple Microgrids under Random Electric Vehicle Charging</article-title>. <source>Energy</source> <volume>208</volume>, <fpage>118360</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2020.118360</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>T. E.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A Model Predictive Control Approach in Microgrid Considering Multi-Uncertainty of Electric Vehicles</article-title>. <source>Renew. Energ.</source> <volume>163</volume>, <fpage>1385</fpage>&#x2013;<lpage>1396</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2020.08.137</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiong</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Modeling and Stability Issues of Voltage-Source Converter Dominated Power Systems: A Review</article-title>. <source>Csee Jpes</source>, <fpage>1</fpage>&#x2013;<lpage>18</lpage>. <pub-id pub-id-type="doi">10.17775/CSEEJPES.2020.03590</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Multi-objective Robust Optimization of Active Distribution Networks Considering Uncertainties of Photovoltaic</article-title>. <source>Int. J.&#x20;Electr. Power Energ. Syst.</source> <volume>133</volume>, <fpage>107197</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2021.107197</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hong</surname>
<given-names>Y. Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Novel Adaptive Multi-Clustering Algorithm-Based Optimal ESS Sizing in Ship Power System Considering Uncertainty [J]</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>33</volume> (<issue>1</issue>), <fpage>307</fpage>&#x2013;<lpage>316</lpage>. </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zarei</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zangeneh</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Multi-objective Optimization Model for Distribution Network Reconfiguration in the Presence of Distributed Generations [J]</article-title>. <source>Int. Trans. Electr. Energ. Syst.</source> <volume>27</volume> (<issue>12</issue>), <fpage>e2425</fpage>. <pub-id pub-id-type="doi">10.1002/etep.2425</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Hierarchical Bipartite Graph Matching Method for Transactive V2V Power Exchange in Distribution Power System [J]</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>12</volume> (<issue>1</issue>), <fpage>301</fpage>&#x2013;<lpage>311</lpage>. <pub-id pub-id-type="doi">10.1109/TSG.2020.3016597</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Targeted Poverty Alleviation Using Photovoltaic Power: Review of Chinese Policies</article-title>. <source>Energy Policy</source> <volume>120</volume>, <fpage>550</fpage>&#x2013;<lpage>558</lpage>. <pub-id pub-id-type="doi">10.1016/j.enpol.2018.06.004</pub-id> </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Or</surname>
<given-names>S. W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Chung</surname>
<given-names>C. Y.</given-names>
</name>
<name>
<surname>Voropai</surname>
<given-names>N. I.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Optimal Coordinated Control of Multi-Renewable-To-Hydrogen Production System for Hydrogen Fueling Stations</article-title>. <source>IEEE Trans. Ind. Applicat.</source>, <fpage>1</fpage>&#x2013;<lpage>1</lpage>. <pub-id pub-id-type="doi">10.1109/TIA.2021.3093841</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>