<?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">852509</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2022.852509</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>Hydrogen-Electric Coupling Coordinated Control Strategy of Multi-Station Integrated System Based on the Honeycomb Topology</article-title>
<alt-title alt-title-type="left-running-head">Wang et al.</alt-title>
<alt-title alt-title-type="right-running-head">Honeycomb Multi-Station</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Wenyong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1247744/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Huo</surname>
<given-names>Qunhai</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1039977/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Ningyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1515775/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ni</surname>
<given-names>Jianfu</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1643861/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhu</surname>
<given-names>Jinda</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yu</surname>
<given-names>Zuofei</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wei</surname>
<given-names>Tongzhen</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Institute of Electrical Engineering</institution>, <institution>Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>University of Chinese Academy of Sciences</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>NARI Technology Co., Ltd.</institution>, <addr-line>Nanjing</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/993078/overview">Haoran Ji</ext-link>, Tianjin University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1188255/overview">Shenxi Zhang</ext-link>, Shanghai Jiao Tong University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1664569/overview">Jie Jian</ext-link>, Tianjin University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Qunhai Huo, <email>huoqunhai@mail.iee.ac.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>31</day>
<month>03</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>852509</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>01</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>07</day>
<month>03</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Wang, Huo, Zhang, Ni, Zhu, Yu and Wei.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Wang, Huo, Zhang, Ni, Zhu, Yu and Wei</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>With the high-proportion accession of renewable energy and randomness of the load side in the new energy power system, unbalanced feeder power and heavy overload of the transformer caused by massive access of highly uncertain source loads become more and more serious. In order to solve the aforementioned problems, a honeycomb topology of the multi-station integrated system is proposed. The soft open point (SOP) is used as the key integrated equipment of the internal unit of a multi-station integrated system. The honeycomb grid structure is composed of flexible nodes, and the multi-station integrated system is composed of multi-network flexible interconnection. Based on the characteristics of the regional resource endowment, hydrogen energy flow is deeply coupled in parts of honeycomb grids. In order to improve the reliability and flexibility of the multi-station integrated unit, the structure of the new multi-station integrated unit, the power balance constraints on the unit, and the switching process of SOP control mode are studied. At the same time, the hydrogen electricity coupling structure and the coordinated control strategy of hydrogen electricity conversion are proposed to solve the problem of deep application of hydrogen energy. Finally, the effectiveness of the proposed multi-station integrated system is verified by using three simulation models.</p>
</abstract>
<kwd-group>
<kwd>honeycomb topology</kwd>
<kwd>soft open point</kwd>
<kwd>multi-station integration</kwd>
<kwd>hydrogen electric coupling</kwd>
<kwd>coordination control</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>With the development of energy reform, the development of high-proportion renewable energy accession and high-proportion power electronic equipment application, and supply side randomness and demand side randomness in power system is accelerated (<xref ref-type="bibr" rid="B2">Ding et al., 2020</xref>). With the increasing plug and play loads connected to the power grid and the increasing proportion of distributed power generation and electric vehicles, the operation of the traditional power grid is faced with a series of problems, such as feeder power imbalance, heavy overload of transformer, difficulty in capacity expansion of transmission, and distribution network (<xref ref-type="bibr" rid="B7">Liu et al., 2020</xref>; <xref ref-type="bibr" rid="B10">Sarantakos et al., 2021</xref>). Meanwhile, the conflict between the randomness of large-scale clean energy and the demand for power quality of new electrical equipment has intensified (<xref ref-type="bibr" rid="B9">Mudaliyar and Mishra 2021</xref>; <xref ref-type="bibr" rid="B12">Sun et al., 2021</xref>; <xref ref-type="bibr" rid="B17">Yang et al., 2021</xref>). Therefore, it is urgent to upgrade the network structure of the existing power grid to solve the aforementioned problems.</p>
<p>In recent years, researchers have proposed many solutions, among which multi-station integration is an important direction (<xref ref-type="bibr" rid="B18">Zhang et al., 2019</xref>; <xref ref-type="bibr" rid="B16">Xincong et al., 2020</xref>; <xref ref-type="bibr" rid="B19">Zhang et al., 2020</xref>; <xref ref-type="bibr" rid="B5">Li et al., 2021</xref>). Multi-station integration is &#x201c;a system that integrates multiple functional stations related to energy and information communication adjacent to geographical locations such as distributed renewable energy power station, electric vehicle charging station, energy storage power station, data center, and 5G base station.&#x201d; At present, the research on the multi-station integrated system is still in its infancy. The problem of solving the integration of four functional stations through mixed integer linear programing (MILP), namely, fast charging stations, plug-in electric vehicles, renewable energy, and energy storage systems is analyzed (<xref ref-type="bibr" rid="B8">Moradzadeh and Abdelaziz 2020</xref>); the low-cost operation problem of the integration of three functional stations of electric bus charging stations, photovoltaic, and battery integrated energy storage system is studied (<xref ref-type="bibr" rid="B22">Zhuang and Liang 2021</xref>); a multi-index comprehensive evaluation method based on the analytic hierarchy process and risk entropy weight is proposed (<xref ref-type="bibr" rid="B21">Zhu et al., 2022</xref>), which verifies that the proposed comprehensive evaluation method can objectively and comprehensively evaluate the operation of the multi-station integrated energy system and assist it to optimize the operation mode and strategy; A day-ahead optimal operation model of the multi-station integration system is proposed (<xref ref-type="bibr" rid="B20">Zhu et al., 2021</xref>), which can improve the operation economy of integration station, while reducing carbon emissions and ensuring power supply reliability. A multi-station integrated duopoly profit distribution mechanism based on the Stackelberg model is constructed (<xref ref-type="bibr" rid="B15">Xie et al., 2021</xref>). Most of the existing studies on the multi-station integration are simple and nearly a combination of multiple functional stations, focusing on research scenarios, coordination and optimization, but the research depth of integration and integration between the functional stations is insufficient. The topology design does not fully consider the high-quality power supply and consumption demand of the functional stations. There are still deficiencies in the power flow mutual assistance of multi-functional stations in the system, and the complementary characteristics between the multiple functional stations are not fully brought into play. Compared with other energy storage methods, hydrogen energy storage has obvious advantages in large-scale and long-term energy storage. However, the overall conversion efficiency of &#x201c;electricity to hydrogen to electricity&#x201d; is low, and the application scenario and efficient conversion mode of distributed generation hydrogen production are not clear.</p>
<p>The optimal operation of SOP in the distribution networks can not only solve feeder imbalance but also improve the penetration of distributed generation and realize the flexible operation of the distribution network (<xref ref-type="bibr" rid="B3">Ji et al., 2019</xref>; <xref ref-type="bibr" rid="B6">Li et al., 2019</xref>; <xref ref-type="bibr" rid="B1">Chen et al., 2020</xref>; <xref ref-type="bibr" rid="B11">Sun et al., 2020</xref>; <xref ref-type="bibr" rid="B14">Wang et al., 2020</xref>). To solve the aforementioned problems, a new honeycomb multi-station integrated system, which integrates various functional stations through SOP, is proposed, and a mesh-flexible interconnected multi-station integrated system centered on flexible nodes is established. Taking the micro-energy network as the basic unit of the multi-station integrated system, the integration of electric energy flow and hydrogen energy flow is realized in the unit. The new multi-station integration unit structure, power to balance constraint, and control mode switching of SOP are presented. At the same time, the hydrogen electric coupling structure and coordinated control strategy are proposed to realize the deep application of hydrogen electric coupling. Finally, three simulation examples verify the effectiveness and superiority in the honeycomb topology based on SOP. The major contributions of this study are as follows:<list list-type="simple">
<list-item>
<p>(1) Based on the honeycomb distribution network, a honeycomb multi-station integrated topology is proposed, which integrates various functional stations through SOP. The system stability and networking flexibility are improved, and the consumption of renewable energy is promoted.</p>
</list-item>
<list-item>
<p>(2) Considering hydrogen energy as a flexible resource, the flow and conversion of hydrogen-electric energy in the honeycomb distribution network are analyzed. Then, the coordinated control strategy of hydrogen electric coupling power is proposed, which provides a solution for the deep coupling problem of &#x201c;electricity&#x2013;energy flow&#x201d; and &#x201c;hydrogen&#x2013;energy flow.&#x201d;</p>
</list-item>
<list-item>
<p>(3) Three simulation examples verify that the proposed honeycomb multi-station integrated system can effectively improve the stability of the system and promote the consumption of renewable energy.</p>
</list-item>
</list>
</p>
</sec>
<sec id="s2">
<title>2 Distribution System of Multi-Station Integration Based on Honeycomb Topology</title>
<sec id="s2-1">
<title>2.1 New Honeycomb Topology</title>
<p>This article continues the idea of the honeycomb power grid (<xref ref-type="bibr" rid="B13">Wang et al., 2018</xref>; <xref ref-type="bibr" rid="B4">Jiang et al., 2019</xref>) and proposes a flexible interconnected hydrogen electric coupling distribution grid structure, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. In the honeycomb topology, a multi-station integrated unit is established with the flexible node as the central area, and multiple units form a ring network structure through the interconnection of the honeycomb power flow. The multi-station integrated unit can adopt AC, DC, or AC/DC hybrid power distributed mode, which has a series of advantages, such as high power supply reliability, flexible power distributed mode, convenient access to new energy, and so on. As an important core energy station, hydrogen energy is deeply integrated with the honeycomb power grid. At the same time, energy storage power station and hydrogen energy station participate in the coordinated operation of the whole energy grid system as flexible resources. The reliability analysis of this topology is shown in <xref ref-type="sec" rid="s11">Supplementary Appendix SA1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Schematic diagram of honeycomb grid topology.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g001.tif"/>
</fig>
<p>The honeycomb multi-station integrated system based on SOP has the following characteristics:<list list-type="simple">
<list-item>
<p>(1) It has high compatibility and can access various types of distributed power and AC/DC loads, which is conducive to the upgrading of the existing power grid.</p>
</list-item>
<list-item>
<p>(2) Each multi-station integrated unit is a relatively independent autonomous system. The basic self-balance of supply and demand in the unit to the greatest extent can be realized, and the power can be freely exchanged with adjacent units through flexible equipment.</p>
</list-item>
<list-item>
<p>(3) The networking mode is flexible and can adapt to the characteristics of various load types of power grid in the future, which can be connected separately by AC or DC or by mixed AC and DC.</p>
</list-item>
<list-item>
<p>(4) The power supply reliability is high. When the system has small disturbance, there are more power sources that can provide power support for the fault side than the traditional power grid. Large disturbance in the traditional power grid will lead to large-area power failure, while the honeycomb power grid only needs to isolate the fault area.</p>
</list-item>
<list-item>
<p>(5) The ability to absorb renewable energy is strong, which can avoid the voltage fluctuation caused by the randomness and indirectness of renewable energy power generation, and the local absorption and grid connection collection of renewable energy reduce the long-distance transmission of electric energy.</p>
</list-item>
</list>
</p>
</sec>
<sec id="s2-2">
<title>2.2 Structure of Hydrogen Electric Coupling</title>
<p>With the rapid development of new energy, the current power system energy flow has gradually expanded on one-way flow to two-way flow, and from single energy access to the power grid to multi-state access to the power grid. The power system is changing from single energy carrier to multi energy carrier, which brings great challenges to the operation of the distribution system. The renewable utilization rate of the electrolytic water hydrogen production device at a larger time scale should be improved. The &#x201c;green hydrogen&#x201d; technology of hydrogen production of renewable energy to deal with the volatility and consumption of wind power and photovoltaic power generation has become an important path for high-efficiency consumption of renewable energy. It has become an important research topic to maximize the local utilization of renewable energy through hydrogen electric coupling technology.</p>
<p>Based on SOP, wind power station, photovoltaic power station, energy storage power station, hydrogen production station, hydrogen gas turbine, and electric vehicle, this study couples the electric energy flow and hydrogen energy flow through the multi-station integrated system. The schematic diagram of the internal hydrogen electric energy flow and conversion is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. The honeycomb multi-station integrated system converts the new energy that cannot be absorbed by the power grid or cannot be easily used by the power grid into the hydrogen energy storage through &#x201c;hydrogen energy flow&#x201d; so as to provide hydrogen for fuel cell electric vehicles and gas engines and realize 100% local absorption of new energy.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Schematic diagram of internal hydrogen/electric energy flowed and conversion.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g002.tif"/>
</fig>
<p>As an important part of the honeycomb multi-station integrated hydrogen electric coupling system, wind power station and photovoltaic power station provide local energy for the load. When the electric energy generated by the wind power station and photovoltaic power station system is surplus, the surplus electric energy is first obtained by electrolyzing water in the hydrogen production station and stored in the hydrogen storage tank, and then stored in the energy storage power station. If the electric energy generated by the power generation system is insufficient to meet the power demand of users, energy will be supplied by the auxiliary power generation system of gas turbine and energy storage power station. If the user load cannot be met, the load demand of the system can be met through the mutual assistance of power flow of adjacent multi-station integrated units.</p>
</sec>
<sec id="s2-3">
<title>2.3 Structure of Multi-Station Integrated Unit</title>
<p>As the key equipment of the multi-station integrated system, the three-port SOP highly integrates multiple functional stations (distributed renewable energy power station, electric vehicle charging station, energy storage power station, hydrogen energy station, etc.) in the multi-station integrated unit, as shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Multi-station integrated subsystem based on the three-port SOP.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g003.tif"/>
</fig>
<p>In practical engineering application, hydrogen energy can be reasonably allocated according to local resources, and electric energy can form various distribution structures such as &#x201c;three-station integration&#x201d; and &#x201c;five-station integration&#x201d; according to the actual characteristics of the existing power grid. The multi-station integrated unit has the structural characteristics of diversification, replication, and expansion.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Control Strategy Analysis</title>
<sec id="s3-1">
<title>3.1 Soft Open Point Control</title>
<p>The SOP has different control modes according to different operating states of the system. The switching of SOPs&#x2019; control mode in this study is shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Three-port SOP control mode.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Mode</th>
<th align="center">VSC1</th>
<th align="center">VSC2</th>
<th align="center">VSC3</th>
<th align="center">State</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Mode 1</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>dc</italic>
</sub>
<italic>&#x2212;Q</italic>
</td>
<td align="left">
<italic>P&#x2212;Q</italic>
</td>
<td align="left">
<italic>P&#x2212;Q</italic>
</td>
<td align="left">Normal</td>
</tr>
<tr>
<td align="left">Mode 2</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>ac</italic>
</sub>
<italic>&#x2212;f</italic>
</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>dc</italic>
</sub>
<italic>&#x2212;Q/P&#x2212;Q</italic>
</td>
<td align="left">
<italic>P&#x2212;Q/U</italic>
<sub>
<italic>dc</italic>
</sub>
<italic>&#x2212;Q</italic>
</td>
<td align="left">bus 1 fault</td>
</tr>
<tr>
<td align="left">Mode 3</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>dc</italic>
</sub>
<italic>&#x2212;Q</italic>
</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>ac</italic>
</sub>
<italic>&#x2212;f</italic>
</td>
<td align="left">
<italic>P&#x2212;Q</italic>
</td>
<td align="left">bus 2 fault</td>
</tr>
<tr>
<td align="left">Mode 4</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>dc</italic>
</sub>
<italic>&#x2212;Q</italic>
</td>
<td align="left">
<italic>P&#x2212;Q</italic>
</td>
<td align="left">
<italic>U</italic>
<sub>
<italic>ac</italic>
</sub>
<italic>&#x2212;f</italic>
</td>
<td align="left">bus 3 fault</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In case of the feeder fault, the port on the fault side needs to switch to <italic>U</italic>
<sub>
<italic>ac</italic>
</sub>&#x2212;<italic>f</italic> control mode to ensure the continuous power operation of the important loads. If the port on the fault side is in the <italic>U</italic>
<sub>
<italic>dc</italic>
</sub>&#x2212;<italic>Q</italic> control mode, the non-fault port needs to be switched to the <italic>U</italic>
<sub>
<italic>dc</italic>
</sub>&#x2212;<italic>Q</italic> control mode. After the feeder fault is eliminated, the control mode is switched to the normal operation state. When switching to the control mode, it is necessary to consider the control mode switching strategy to avoid the transient impact of the converter output current, resulting in oscillation and impact on the operation stated switching of the multi-station integrated unit.</p>
<p>The switching control method of controller state tracking and introducing inertia link is adopted. The controller block diagram is shown in <xref ref-type="fig" rid="F4">Figure 4</xref>. (The control strategy in the figure takes VSC1 as an example.)</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>SOP control mode.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g004.tif"/>
</fig>
<p>The controller state tracking is to take <italic>u</italic>
<sub>
<italic>d</italic>
</sub> and <italic>u</italic>
<sub>
<italic>q</italic>
</sub> of the current control mode as the input signals of the system, and <italic>u</italic>
<sub>
<italic>d</italic>
</sub> and <italic>u</italic>
<sub>
<italic>q</italic>
</sub> of the other control modes track <italic>u</italic>
<sub>
<italic>d&#x2a;_ref</italic>
</sub> and <italic>u</italic>
<sub>
<italic>q&#x2a;_ref</italic>
</sub> in real time, respectively, so as to ensure that the output voltage can match when the converter is switched from one working mode to another, and reduce the oscillation or overshoot caused by the system operation mode switching. Finally, the smooth switching of the three port SOP working mode can be realized. At the same time, the inertia link is introduced into the controller, which can increase or decrease the hysteresis of <italic>i</italic>
<sub>
<italic>ref</italic>
</sub> during system switching, so as to suppress the impact caused by the sudden change of <italic>i</italic>
<sub>
<italic>ref</italic>
</sub>, and make the reference current of the inner loop gradually adjust to the required value of the constant voltage and constant frequency. In addition, the power flowed to transfer command and logic switch switching command of the SOP controller are provided by the control center of the multi-station integration system, so as to coordinate with other functional stations.</p>
</sec>
<sec id="s3-2">
<title>3.2 Power Balanced Constraint Model of Multi-Station Integrated Unit</title>
<p>In order to realize the stable operation of the multi-station fusion unit, this section gives the power model of each functional station in the multi-station fusion unit. The model contains many integer variables to indicate the operational state of resources and many continuous variables, such as the power flow variables.</p>
<sec id="s3-2-1">
<title>3.2.1 Output Model of Photovoltaic</title>
<p>The output power function of the photovoltaic is as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>T</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>In the formula, P<sub>STC</sub> is the output power of photovoltaic panels under standard conditions; <italic>k</italic>
<sub>
<italic>c</italic>
</sub> is the temperature coefficient; <italic>T</italic>
<sub>
<italic>PV</italic>
</sub> is the actual temperature of the photovoltaic panels; <italic>T</italic>
<sub>
<italic>ref</italic>
</sub> is the standard ambient temperature; <italic>A</italic>
<sub>
<italic>r</italic>
</sub> is the area of the photovoltaic panel; and <italic>&#x3b7;</italic>
<sub>
<italic>PV</italic>
</sub> is the efficiency of the photovoltaic panels.</p>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Output Model of Wind Turbine</title>
<p>The output power function of the wind turbine is as follows:<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:msup>
<mml:mi>v</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf1">
<mml:math id="m3">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mi>r</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf2">
<mml:math id="m4">
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mi>r</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>; <italic>P</italic>
<sub>
<italic>r</italic>
</sub>, <italic>v</italic>
<sub>
<italic>ci</italic>
</sub>, <italic>v</italic>
<sub>
<italic>r</italic>
</sub>, and <italic>v</italic>
<sub>
<italic>co</italic>
</sub> are the rated power, cut-in wind speed, rated wind speed, and cut-out wind speed of the fan, respectively.</p>
</sec>
<sec id="s3-2-3">
<title>3.2.3 Model of Energy Storage Equipment</title>
<p>Assuming that the maximum charging power of the energy storage device is equal to the maximum discharge power, it is satisfied in the charging and discharging process:<disp-formula id="e3">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mtext>,</mml:mtext>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m6">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mtext>,&#xa0;</mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>In the formula, <italic>P</italic>
<sub>
<italic>ESS,max</italic>
</sub> and <italic>P</italic>
<sub>
<italic>ESS,min</italic>
</sub> represent the upper and lower limits of charge and discharge power for the energy storage devices, respectively.</p>
<p>The charge state at the end of discharge is equal to the initial charge state.<disp-formula id="e5">
<mml:math id="m7">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<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>S</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>d</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<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>T</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:msubsup>
<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>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
</mml:msubsup>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>In the formula, SOC<sub>
<italic>i</italic>
</sub> and SOC<sub>(i-1)</sub> are the charge state of battery energy storage at time <italic>i</italic> and <italic>i</italic>-1, respectively; <italic>&#x3b7;</italic>
<sub>
<italic>c</italic>
</sub> and <italic>&#x3b7;</italic>
<sub>
<italic>d</italic>
</sub> are the charging efficiency and discharging efficiency of the battery energy storage, respectively.</p>
</sec>
<sec id="s3-2-4">
<title>3.2.4 Models of IDC and 5G</title>
<p>To simplify the analysis, IDC and 5G are treated as the constant power models.<disp-formula id="e6">
<mml:math id="m8">
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mn>5</mml:mn>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mn>5</mml:mn>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
</sec>
<sec id="s3-2-5">
<title>3.2.5 Model of Electrolytic Cell</title>
<p>The hydrogen&#x2013;electric conversion device consists of three parts: electrolytic cell, hydrogen storage tank, and hydrogen gas engine.</p>
<p>The output power of the electrolyzer can be expressed as follows:<disp-formula id="e7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>In the formula, P<sub>GE</sub>(t) is the input power of the cell and &#x3b7;<sub>GE</sub> is the working efficiency of the electrolyzer.</p>
</sec>
<sec id="s3-2-6">
<title>3.2.6 Model of Hydrogen Gas Engine</title>
<p>The output power of the hydrogen gas engine can be expressed as follows:<disp-formula id="e8">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</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:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>H</mml:mi>
<mml:mi>E</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:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <italic>P</italic>
<sub>
<italic>tank-HES</italic>
</sub> (t) is the input power of the hydrogen gas engine and <italic>&#x3b7;</italic>
<sub>
<italic>HES</italic>
</sub> is the efficiency of the hydrogen gas engine.</p>
</sec>
<sec id="s3-2-7">
<title>3.2.7 Model of Hydrogen Storage Tank</title>
<p>A hydrogen storage tank is used to store hydrogen produced by the electrolytic water. It can also provide hydrogen for the hydrogen gas engine and improve system flexibility. The mathematical model of the energy storage can be expressed as follows:<disp-formula id="e9">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>where <italic>E</italic>
<sub>
<italic>tank</italic>
</sub>(<italic>t</italic>) is the energy stored in the hydrogen storage tank at time t, E<sub>tank</sub> (<italic>t</italic>-1) is the energy stored in the hydrogen storage tank at time t-1, <italic>&#x3b7;</italic>
<sub>
<italic>tank</italic>
</sub> is the working efficiency of the hydrogen storage tank, and &#x3b7;<sub>GE_DC/DC</sub> and &#x3b7;<sub>HES_DC/DC</sub> are the working efficiency of DC/DC at the port side of the electrolyzer and hydrogen gas engine, respectively.</p>
<p>The maximum output power of the electrolytic cell and the maximum output power of the hydrogen gas engine are limited by their capacity and the remaining storage energy of the hydrogen storage tank, which are expressed as follows:<disp-formula id="e10">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</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:mi mathvariant="normal">&#x394;</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:msub>
<mml:mtext>&#xa0;and</mml:mtext>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
<disp-formula id="e11">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where C<sub>GE</sub> and <italic>C</italic>
<sub>
<italic>HES</italic>
</sub> are the capacities of the electrolyzer and hydrogen gas engine, respectively; <italic>E</italic>
<sub>
<italic>tank_min</italic>
</sub> and <italic>E</italic>
<sub>
<italic>tank_max</italic>
</sub> are the lower and upper limits of the storage energy of the hydrogen storage tank, respectively.</p>
<p>In the multi-station integrated unit shown in <xref ref-type="fig" rid="F3">Figure 3</xref>, <italic>P</italic>
<sub>MSI</sub>(<italic>t</italic>), <italic>P</italic>
<sub>IDC</sub>(<italic>t</italic>), <italic>P</italic>
<sub>VSC</sub>(<italic>t</italic>), <italic>P</italic>
<sub>load</sub>(<italic>t</italic>), <italic>P</italic>
<sub>5G</sub>(<italic>t</italic>), <italic>P</italic>
<sub>ESS</sub>(<italic>t</italic>), <italic>P</italic>
<sub>HES</sub>(<italic>t</italic>), <italic>P</italic>
<sub>GE</sub>(<italic>t</italic>), <italic>P</italic>
<sub>PV</sub>(<italic>t</italic>), <italic>P</italic>
<sub>PEV</sub>(<italic>t</italic>), and <italic>P</italic>
<sub>HFEV</sub>(<italic>t</italic>), respectively, represent the exchange power between the multi-station integrated system and the adjacent multi-station integrated system at t time, the load power of the data center, the input/output power of each port of SOP, the power of the conventional switchable load, the load power of the 5G base station, the charge/discharge power of the energy storage power station, the power consumed by hydrogen production in the hydrogen energy station, the output power of the gas engine, the output power of the photovoltaic power station, charge/discharge power of the pure electric vehicle, and the absorbed power of the hydrogen fuel electric vehicle; <italic>&#x3be;</italic>
<sub>load1/2/3</sub> is the load-shedding rate of conventional load; <italic>&#x3bb;</italic>
<sub>MSI</sub>&#x3001;<italic>&#x3bb;</italic>
<sub>VSC</sub>&#x3001;<italic>&#x3bb;</italic>
<sub>ESS</sub>&#x3001;<italic>&#x3bb;</italic>
<sub>HES</sub>&#x3001;<italic>&#x3bb;</italic>
<sub>GE</sub>&#x3001;<italic>&#x3bb;</italic>
<sub>PV</sub>&#x3001;<italic>&#x3bb;</italic>
<sub>PEV</sub>&#x3001;and <italic>&#x3bb;</italic>
<sub>HFEV</sub> are 0/1 variable, which can be expressed as follows:<disp-formula id="e12">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>i</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>Exchange</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>power</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>with</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>not</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>exchange</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>power</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>with</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>,</mml:mtext>
<mml:mn>2</mml:mn>
<mml:mtext>,</mml:mtext>
<mml:mn>3</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>i</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mn>1</mml:mn>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>VSC</mml:mtext>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>port</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>normal</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>VSC</mml:mtext>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>port</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>fault</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>,</mml:mtext>
<mml:mn>2</mml:mn>
<mml:mtext>,</mml:mtext>
<mml:mn>3</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>i</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mn>1</mml:mn>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>EES</mml:mtext>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>charge</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>or</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>discharge</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>EES</mml:mtext>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>,</mml:mtext>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="e15">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>i</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>HES</mml:mtext>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>production</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>HES</mml:mtext>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</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:mo>(</mml:mo>
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mtext>,</mml:mtext>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>GE</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>normal</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>GE</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</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>
<disp-formula id="e17">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>PV</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>normal</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>PV</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
<disp-formula id="e18">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>WT</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>normal</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>WT</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
<disp-formula id="e19">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>PEV</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>charge</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>or</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>discharge</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>PEV</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>is</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
<disp-formula id="e20">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</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:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>1</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>HFEV</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>hydrogen</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>consumption</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mtext>HFEV</mml:mtext>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mtext>standby</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(20)</label>
</disp-formula>
</p>
<p>The power models on the three feeders of SOP can be expressed as follows:<disp-formula id="e21">
<mml:math id="m23">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(21)</label>
</disp-formula>
<disp-formula id="e22">
<mml:math id="m24">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mn>5</mml:mn>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(22)</label>
</disp-formula>
<disp-formula id="e23">
<mml:math id="m25">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x3d;</mml:mo>
<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:mrow>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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>&#x2b;</mml:mo>
<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:mrow>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(23)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e23">eq. 23</xref>, M<sub>1</sub> represents the number of pure electric vehicles and M<sub>2</sub> represents the number of hydrogen-fueled electric vehicles. In this article, the feeder line loss and power electronic converter loss are ignored.</p>
<p>The multi-station integrated system should also meet the hydrogen energy power conservation:<disp-formula id="e24">
<mml:math id="m26">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<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:mrow>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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:mtable columnalign="left">
<mml:mtr columnalign="left">
<mml:mtd columnalign="left">
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(24)</label>
</disp-formula>
</p>
<p>Among them, &#x3b1; represents the energy conversion rate of electric hydrogen production; &#x3b2; represents the utilization rate of hydrogen energy, that is, the remaining hydrogen energy in the hydrogen pipeline can be expressed as (1-<italic>&#x3b2;</italic>)[<italic>&#x3bb;</italic>
<sub>HES_1</sub>
<italic>P</italic>
<sub>HES_1</sub>(<italic>t</italic>)&#x2b;<italic>&#x3bb;</italic>
<sub>HES_2</sub>
<italic>P</italic>
<sub>HES_2</sub>(<italic>t</italic>)].</p>
<p>SOP should meet the active power balance constraint. Compared with the system loss, the active power loss of VSC<italic>i</italic> is very small and can be ignored in the optimization process. The active power emitted or absorbed meets the following equation:<disp-formula id="e25">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">VSC_1</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="italic">&#x2b;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">VSC_2</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="italic">&#x2b;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">VSC_3</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="italic">t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="italic">&#x3d;0</mml:mi>
<mml:mi mathvariant="italic">.</mml:mi>
</mml:mrow>
</mml:math>
<label>(25)</label>
</disp-formula>
</p>
<p>Substituting <xref ref-type="disp-formula" rid="e21">Eqs 21</xref>&#x2013;<xref ref-type="disp-formula" rid="e24">24</xref> into <xref ref-type="disp-formula" rid="e25">Eq. 25</xref>, the power model of the multi-station integrated system can be obtained as follows:<disp-formula id="e26">
<mml:math id="m28">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mn>5</mml:mn>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<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:mn>2</mml:mn>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>_</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<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:mrow>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3b1;</mml:mi>
</mml:mfrac>
<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:mrow>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<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:mn>3</mml:mn>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3be;</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>V</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>W</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<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:mn>3</mml:mn>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>j</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:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>&#x3b1;</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(26)</label>
</disp-formula>
</p>
<p>At the same time, each functional station involved in the multi-station integrated system also needs to meet its corresponding operation constraints.</p>
</sec>
</sec>
<sec id="s3-3">
<title>3.3 Coordinated Control Strategy of Hydrogen Electric Coupling System</title>
<p>Considering the actual dynamic response delay characteristics of the electrolytic cell and gas engine, when the pressure <italic>p</italic>
<sub>H2</sub> of the hydrogen storage tank is within the normal range, both the electrolytic cell and gas engine can realize power transfer according to their respective power reference values and aiming at balancing the power difference between the supply and demand. At this time, the main function of the energy storage power station is to stabilize the slow dynamic response process when the power of the electrolytic cell and gas engine changes by its good transient performance. The power coordination control strategy flow of the honeycomb multi-station integrated hydrogen electric coupling system is shown in <xref ref-type="fig" rid="F5">Figure 5</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Flowchart of the system power coordination control strategy.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g005.tif"/>
</fig>
<p>When the PV/wind output is less than the load demand (<italic>&#x394;P</italic> &#x3c; 0) and the pressure <italic>p</italic>
<sub>H2</sub> of the hydrogen storage tank is within the normal range, the hydrogen production station is shut down (<italic>P</italic>
<sub>
<italic>HES</italic>
</sub> &#x3d; 0), and the gas turbine is used as the auxiliary power supply (<italic>P</italic>
<sub>
<italic>GE_ref</italic>
</sub> &#x3d; &#x7c;<italic>&#x394;P&#x7c;</italic>). If the electric field energy reserve of the energy storage power station is sufficient at this time (<italic>U</italic>
<sub>
<italic>ESS</italic>
</sub> &#x3e; <italic>U</italic>
<sub>
<italic>ESS_min</italic>
</sub>), the energy storage power station will release energy to make up for the power shortage of the dynamic process of the gas engine (mode 1). If the energy storage power station is not full at this time (<italic>U</italic>
<sub>
<italic>ESS</italic>
</sub> &#x3c; <italic>U</italic>
<sub>
<italic>ESS_max</italic>
</sub>), the energy storage power station will absorb the residual energy of the system to stabilize the residual power in the dynamic process of the gas engine (mode 2).</p>
<p>When the PV/wind output is greater than the load demand (<italic>&#x394;P</italic> &#x3e; 0), the pressure <italic>p</italic>
<sub>H2</sub> of the hydrogen storage tank is within the normal range. At this time, the gas turbine is shut down (<italic>P</italic>
<sub>
<italic>HES</italic>
</sub> &#x3d; 0), and the hydrogen production station converts the residual energy of the system into hydrogen and stores it in the hydrogen storage tank (<italic>P</italic>
<sub>
<italic>HES_ref</italic>
</sub> &#x3d; <italic>&#x394;P</italic>). If the electric field energy reserve of the energy storage power station is sufficient at this time (<italic>U</italic>
<sub>
<italic>ESS</italic>
</sub> &#x3e; <italic>U</italic>
<sub>
<italic>ESS_min</italic>
</sub>), the energy storage power station will release energy to make up for the power shortage of the dynamic process of the gas engine (mode 3). If the energy storage power station is not full at this time (<italic>U</italic>
<sub>
<italic>ESS</italic>
</sub> &#x3c; <italic>U</italic>
<sub>
<italic>ESS_max</italic>
</sub>), the energy storage power station will absorb the system residual energy to balance the residual power in the dynamic process of the hydrogen production station (mode 4).</p>
<p>When the PV/wind output is less than the load demand (<italic>&#x394;P</italic> &#x3c; 0), the pressure <italic>p</italic>
<sub>H2</sub> of the hydrogen storage tank is within the abnormal range. At this time, the hydrogen production station is shut down (<italic>P</italic>
<sub>
<italic>HES</italic>
</sub> &#x3d; 0), and the gas engine is shut down (P<sub>GE</sub> &#x3d; 0). If the electric field energy reserve of the energy storage power station is sufficient (<italic>U</italic>
<sub>
<italic>ESS</italic>
</sub> &#x3e; <italic>U</italic>
<sub>
<italic>ESS_min</italic>
</sub>), the energy storage power station will release energy to balance the power shortage of the system (mode 5). If the energy storage power station is not full at this time (<italic>U</italic>
<sub>
<italic>ESS</italic>
</sub> &#x3c; <italic>U</italic>
<sub>
<italic>ESS_max</italic>
</sub>), the energy storage power station will absorb the residual power of the system (mode 6). If the aforementioned mode conditions cannot be met, the energy storage plant will also shut down (<italic>P</italic>
<sub>
<italic>ESS_ref</italic>
</sub> &#x3d; 0) (mode 7).</p>
</sec>
</sec>
<sec id="s4">
<title>4 Example Simulation Verification and Analysis</title>
<p>In order to verify the feasibility of the proposed multi-station integrated system topology and control strategy, a simulation model of the hydrogen electric coupling multi-station integrated system based on three-port SOP is built in Matlab/Simulink environment, and its specific dynamic response characteristics are tested and analyzed.</p>
<sec id="s4-1">
<title>4.1 Example 1: Operation Verification of Three-Port Soft Open Point</title>
<p>VSC<italic>i</italic> (<italic>i</italic> &#x3d; 1,2,3) under the normal operation of the three-port SOP is controlled by U<sub>dc</sub>&#x2212;Q, P&#x2212;Q, and P&#x2212;Q, respectively. Setting 1.0 s to start the three-port SOP, the control center provides control instructions. The power flow control response curve of starting SOP is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Power flow control response curve.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g006.tif"/>
</fig>
<p>In <xref ref-type="fig" rid="F6">Figure 6</xref>, after a short period of power flow adjustment in SOP, the power of bus 1, bus 2, and bus 3 will soon reach equilibrium. Through the application of SOP to the multi-station integrated system, the power flow regulation between the distributed networks can be easily realized, which can provide important support for the further coordinated operation of the multi-station integrated system.</p>
</sec>
<sec id="s4-2">
<title>4.2 Example 2: Power Flow Transferring and Power Supply Restoration of Multi-Station Integrated Unit</title>
<p>It is assumed that the active loads of the three ports after ignoring the internal self-balancing part of the micro-grid are 100, 100, and 160&#xa0;kW, respectively. The feeder load balancing is the control target, and the power flowing into the DC side of SOP is positive. During the normal operation, the active power of the three ports is 20, 20, and 40&#xa0;kW, respectively.<list list-type="simple">
<list-item>
<p>(1) At 0.2&#xa0;s, 90-kw load is input to bus 2. After the regulation of SOP, the active power of the three ports is 50, &#x2212;40, and &#x2212;10&#xa0;kW, respectively. At 0.4&#xa0;s, 60-kW load is cut off in bus 3. After the regulation of SOP, the active power of the three ports is 30, &#x2212;60, and 30&#xa0;kW, respectively. <xref ref-type="fig" rid="F7">Figure 7</xref> shows the system load on/off simulation waveform.</p>
</list-item>
</list>
</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Simulation waveform of the load on/off. <bold>(A)</bold> SOP Bus voltage at dc side. <bold>(B)</bold> Three phase current of <italic>I<sub>a</sub>, I<sub>b</sub>, I<sub>c</sub>
</italic> (Bus1 as an example). <bold>(C)</bold> Active power of VSC1, VSC2, VSC3.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g007.tif"/>
</fig>
<p>The load on/off simulation waveform in <xref ref-type="fig" rid="F7">Figure 7</xref> shows that when the load is input or removed, the SOP can realize the rapid power flow transfer in the system, and the multi-station integration unit can reach a stable operation state within 20&#xa0;ms.<list list-type="simple">
<list-item>
<p>(2) When it is set at 0.3&#xa0;s, the upper power grid connected to bus 3 fails and the 100-kW load is cut off. At this time, VSC3 is switched from P&#x2212;Q control to U<sub>ac</sub>&#x2212;<italic>f</italic> control. <xref ref-type="fig" rid="F8">Figure 8</xref> shows the simulation waveform of load voltage on the fault side.</p>
</list-item>
</list>
</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Simulation waveform of the load voltage at the fault side.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g008.tif"/>
</fig>
<p>This example verifies the effectiveness of the SOP control mode switching strategy. It can be seen from the simulation waveform in <xref ref-type="fig" rid="F8">Figure 8</xref> that under the action of SOP smooth switching control, when the power grid fails, the feeder voltage at the fault side can quickly recover and stabilize, and the uninterrupted and stable power supply of power loss load in the multi-station integration unit can be realized. At the same time, the feeder voltage on the fault side also quickly reaches a stable state. From this, we can conclude that the power flowed transfer and rapid power supply recovery of the multi-station integration unit can be realized through the power regulation of SOP, so as to ensure the flexibility and stability of the system operation. In addition, regulating the multi-station integration unit through SOP can provide basic guarantee of the coordination and interaction among the multiple functional stations.</p>
</sec>
<sec id="s4-3">
<title>4.3 Example 3: Energy Flow and Realization Process of Hydrogen Electric Coupling System</title>
<p>This example analysis adopts the typical scenario of the multi-station integration system, as shown in <xref ref-type="fig" rid="F9">Figure 9</xref>, which includes the three-port SOP, renewable energy power stations A1-A3, sensitive load B0, conventional loads B1-B3, bus 1-4, energy storage power station E1, hydrogen production station H1, hydrogen gas turbine G1, and hydrogen storage station.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Case of scenario.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g009.tif"/>
</fig>
<p>During the normal operation of bus 1-3, the functional stations are in independent operation or partially coupled operation, and the energy utilization efficiency is not high. In case of power grid failure, the existing equipment on the line is relied on to provide a short-term power support. After bus 1-3 are interconnected with SOP and hydrogen electric coupling equipment is added, one-way flow of energy to multi-directional flow can be changed, and the single station operation switches to the multi-station coordinated operation. In this example, 100% new energy generation is considered, and the power loss of power equipment and lines is ignored. Based on the one operation day, the operation conditions of three bus 1-3 are given, respectively, as shown in <xref ref-type="fig" rid="F10">Figures 10</xref>&#x2013;<xref ref-type="fig" rid="F12">12</xref>.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Operation of bus 1.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Operation of bus 2.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g011.tif"/>
</fig>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Operation of bus 3.</p>
</caption>
<graphic xlink:href="fenrg-10-852509-g012.tif"/>
</fig>
<p>The daily power consumption loaded power curve of the load B0-B3 is shown in <xref ref-type="sec" rid="s11">Supplementary Appendix SA2</xref>. The daily output curve of A1-A3 renewable energy powered station is shown in <xref ref-type="sec" rid="s11">Supplementary Appendix SA3</xref>. The minimum transferable power of SOP port is set to 10&#xa0;kW. When the power is less than 10&#xa0;kW and cannot be transferred, it is absorbed through H1 of the hydrogen production station. Bus 1 grid becomes faulty from 17:00 to 19:00. The example scenario is analyzed as follows:</p>
<p>During the period from 01:00 to 13:00, the power generation of renewable energy in the system is greater than the load. The maximum consumption of renewable energy is realized through the power flow transfer of SOP, and the excess power is stored in the energy storage power station E1 or the energy that cannot be consumed is made into hydrogen through the hydrogen production station H1 for storage. During the period from 14:00 to 16:00, the energy storage power station E1 reaches its charged state, and the excess energy is completely consumed by the hydrogen production station H1. During the period from 17:00 to 19:00, due to the power grid failure of bus 1, the power support is first provided by the discharge of energy storage power station E1. At 19:00, the power of the energy storage power station is insufficient, and the hydrogen combustion gassed engine G1 burns hydrogen to supplement the power. During the period from 20:00 to 24:00, the renewable energy power generation of the system is less than the load. Bus 1-3 provide power through the power grid and makes up for the insufficient power of the system through the power flow transfer of SOP.</p>
<p>Through the example analysis, the integration of renewable energy power stations, energy storage power stations, flexible switching stations, hydrogen production stations and other functional stations, and the deep coupling of electric energy flow and hydrogen energy flow through the multi-station integrated system can achieve 100% consumption of new energy. By participating in the coordinated operation of the whole energy network system as a flexible resource together with the hydrogen energy station, the power supply timed and system power supply reliability of sensitive loads can be improved, and the comprehensive utilization rate of system resources can also be improved.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Conclusion and Prospect</title>
<p>Aiming at the operation flexibility and reliability requirements of the renewable energy power system with double high and double random characteristics, a honeycomb multi-station integrated system composed of SOP is proposed. The multi-station integrated unit structure, power to balance constraint, and SOP control mode switching are studied. At the same time, the hydrogen electric coupling structure and coordinated control strategy are also presented. The flexibility and reliability of the multi-station integrated system is verified through the three simulation examples. The proposed multi-station integrated system can not only adapt to the development trend of the power grid in the future but also has a series of advantages, such as simple structure, easy replication and expansion, reliable power supply, economic operation, flexible control, and so on. Through the in-depth application of the honeycomb topology and SOP, the multi-station integrated system can improve the comprehensive performance of the distribution network, the consumption capacity of distributed generation, and the reliability of the system power supplies.</p>
<p>The research content of this article is expected to promote the practical application process of multi-station integrated technology in the future power grid and accelerate the popularization and application of multi-station integrated technology in the future intelligent distributed network.</p>
</sec>
</body>
<back>
<sec id="s6">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s11">Supplementary Material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author Contributions</title>
<p>WW: conceptualization, methodology, and writing&#x2014;original draft. QH: software and methodology; NZ: methodology; JN, JZ, and ZY: validation; and TW: supervision.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>This work is supported by the National Key R&#x26;D Program of China (No.2021YFB2601604) and the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA21050304).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of Interest</title>
<p>Author JN, JZ, ZY were employed by NARI Technology Co.,Ltd. </p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2022.852509/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2022.852509/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet1.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zha</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>A Novel Node Flexibility Evaluation Method of Active Distribution Network for SNOP Integration</article-title>. <source>IEEE J. Emerging Selected Top. Circuits Syst.</source> <volume>11</volume> (<issue>1</issue>), <fpage>188</fpage>&#x2013;<lpage>198</lpage>. <pub-id pub-id-type="doi">10.1109/JETCAS.2020.3039535</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ding</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Shahidehpour</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Multiperiod Distribution System Restoration with Routing Repair Crews, Mobile Electric Vehicles, and Soft-Open-Point Networked Microgrids</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>11</volume> (<issue>6</issue>), <fpage>4795</fpage>&#x2013;<lpage>4808</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2020.3001952</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ji</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Robust Operation of Soft Open Points in Active Distribution Networks with High Penetration of Photovoltaic Integration</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>10</volume> (<issue>1</issue>), <fpage>280</fpage>&#x2013;<lpage>289</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2018.2833545</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Conception, Key Technology and Prospect of Honeycomb-Shape Active Distribution Network</article-title>. <source>Automation Electric Power Syst.</source> <volume>43</volume> (<issue>17</issue>), <fpage>1</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.7500/AEPS20180730017</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). &#x201c;<article-title>Construction Method of &#x201c;Multi-Station Integration&#x201d; Operation System Based on Ubiquitous Power Internet of Things</article-title>,&#x201d; in <conf-name>2021 International Conference on Artificial Intelligence and Smart Systems (ICAIS)</conf-name>, <fpage>1343</fpage>&#x2013;<lpage>1346</lpage>. <pub-id pub-id-type="doi">10.1109/ICAIS50930.2021.9395800</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Optimal Operation of Soft Open Points in Active Distribution Networks under Three-Phase Unbalanced Conditions</article-title>. <source>IEEE Trans. Smart Grid</source> <volume>10</volume> (<issue>1</issue>), <fpage>380</fpage>&#x2013;<lpage>391</lpage>. <pub-id pub-id-type="doi">10.1109/tsg.2017.2739999</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>A Bi-Level Interval Robust Optimization Model for Service Restoration in Flexible Distribution Networks</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>36</volume> (<issue>3</issue>), <fpage>1843</fpage>&#x2013;<lpage>1855</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2020.3038400</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moradzadeh</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Abdelaziz</surname>
<given-names>M. M. A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>A New MILP Formulation for Renewables and Energy Storage Integration in Fast Charging Stations</article-title>. <source>IEEE Trans. Transp. Electrific.</source> <volume>6</volume> (<issue>1</issue>), <fpage>181</fpage>&#x2013;<lpage>198</lpage>. <pub-id pub-id-type="doi">10.1109/tte.2020.2974179</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mudaliyar</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Mishra</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Real-Time Coordinated Control of Low-Voltage DC Distribution Network with Soft Opening Point</article-title>. <source>IEEE Trans. Power Electron.</source> <volume>36</volume> (<issue>6</issue>), <fpage>7123</fpage>&#x2013;<lpage>7137</lpage>. <pub-id pub-id-type="doi">10.1109/tpel.2020.3037020</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sarantakos</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Zografou-Barredo</surname>
<given-names>N.-M.</given-names>
</name>
<name>
<surname>Huo</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Greenwood</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A Reliability-Based Method to Quantify the Capacity Value of Soft Open Points in Distribution Networks</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>36</volume> (<issue>6</issue>), <fpage>5032</fpage>&#x2013;<lpage>5043</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2021.3071088</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Optimized Two-Time Scale Robust Dispatching Method for the Multi-Terminal Soft Open Point in Unbalanced Active Distribution Networks</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>12</volume> (<issue>1</issue>), <fpage>587</fpage>&#x2013;<lpage>598</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2020.3013386</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Optimized Two-Time Scale Robust Dispatching Method for the Multi-Terminal Soft Open Point in Unbalanced Active Distribution Networks</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>12</volume> (<issue>1</issue>), <fpage>587</fpage>&#x2013;<lpage>598</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2020.3013386</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Development and Characteristic Analysis of Flexibility in Smart Distribution Network</article-title>. <source>Automation Electric Power Syst.</source> <volume>42</volume> (<issue>10</issue>), <fpage>13</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.7500/AEPS20171012002</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chung</surname>
<given-names>C. Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Coordinated Planning of Converter-Based DG Units and Soft Open Points Incorporating Active Management in Unbalanced Distribution Networks</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>11</volume> (<issue>3</issue>), <fpage>2015</fpage>&#x2013;<lpage>2027</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2019.2950168</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Duopoly Benefit Distribution Mechanism Under Multi-Station Integration Based on Stackelberg Model</article-title>. <source>Power Syst. Technol.</source> <volume>45</volume> (<issue>10</issue>), <fpage>4009</fpage>&#x2013;<lpage>4015</lpage>. <pub-id pub-id-type="doi">10.1109/APPEEC48164.2020.9220407</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Xincong</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Mingze</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2020</year>). &#x201c;<article-title>Location Evaluation of Multi-Station Integration Based on AHP-Fuzzy Comprehensive Evaluation Method</article-title>,&#x201d; in <conf-name>2020 12th IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC)</conf-name>, <fpage>1</fpage>&#x2013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1109/APPEEC48164.2020.9220407</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Real-Time Coordinated Scheduling for ADNs with Soft Open Points and Charging Stations</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>36</volume> (<issue>6</issue>), <fpage>5486</fpage>&#x2013;<lpage>5499</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2021.3070036</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). &#x201c;<article-title>Comprehensive Evaluation Method of Safety and Benefit for Multi-Station Integration</article-title>,&#x201d; in <conf-name>2019 IEEE Sustainable Power and Energy Conference (iSPEC)</conf-name>, <fpage>2618</fpage>&#x2013;<lpage>2623</lpage>. <pub-id pub-id-type="doi">10.1109/iSPEC48194.2019.8975335</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). &#x201c;<article-title>Planning Strategy of Multi-Station Integration for Substations in Power Systems</article-title>,&#x201d; in <conf-name>2020 IEEE International Conference on Advances in Electrical Engineering and Computer Applications( AEECA)</conf-name>, <fpage>576</fpage>&#x2013;<lpage>581</lpage>. <pub-id pub-id-type="doi">10.1109/AEECA49918.2020.9213535</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Day-ahead Optimal Operation of Multi-Station Fusion System Considering Carbon Emission</article-title>. <source>Power Demand Side Manage.</source> <volume>23</volume> (<issue>05</issue>), <fpage>40</fpage>&#x2013;<lpage>45</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1009-1831.2021.05.009</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Mu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Dai</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Qian</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Multi-Index Comprehensive Evaluation of Multi-Station Integrated Energy System Based on Analytic Hierarchy Process and Risk Entropy Weigh</article-title>. <comment>Electrical Measurement &#x26; Instrumentation</comment>. <comment>Available at: <ext-link ext-link-type="uri" xlink:href="http://kns.cnki.net/kcms/detail/23.1202.th.%2020220114.1628.005.html">http://kns.cnki.net/kcms/detail/23.1202.th. 20220114.1628.005.html</ext-link> (Accessed 02 06, 2022)</comment>. </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhuang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Stochastic Energy Management of Electric Bus Charging Stations with Renewable Energy Integration and B2G Capabilities</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>12</volume> (<issue>2</issue>), <fpage>1206</fpage>&#x2013;<lpage>1216</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2020.3039758</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>