<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "http://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en" article-type="research-article" dtd-version="1.1">
<front>
<journal-meta>
<journal-id journal-id-type="pmc">EE</journal-id>
<journal-id journal-id-type="nlm-ta">EE</journal-id>
<journal-id journal-id-type="publisher-id">EE</journal-id>
<journal-title-group>
<journal-title>Energy Engineering</journal-title>
</journal-title-group>
<issn pub-type="epub">1546-0118</issn>
<issn pub-type="ppub">0199-8595</issn>
<publisher>
<publisher-name>Tech Science Press</publisher-name>
<publisher-loc>USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">27365</article-id>
<article-id pub-id-type="doi">10.32604/ee.2023.027365</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Virtual Synchronous Generator Adaptive Control of Energy Storage Power Station Based on Physical Constraints</article-title>
<alt-title alt-title-type="left-running-head">Virtual Synchronous Generator Adaptive Control of Energy Storage Power Station Based on Physical Constraints</alt-title>
<alt-title alt-title-type="right-running-head">Virtual Synchronous Generator Adaptive Control of Energy Storage Power Station Based on Physical Constraints</alt-title>
</title-group>
<contrib-group>
<contrib id="author-1" contrib-type="author">
<name name-style="western"><surname>Huang</surname><given-names>Yunfan</given-names></name><xref ref-type="aff" rid="aff-1">1</xref></contrib>
<contrib id="author-2" contrib-type="author">
<name name-style="western"><surname>Lv</surname><given-names>Qingquan</given-names></name><xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-3" contrib-type="author">
<name name-style="western"><surname>Zhang</surname><given-names>Zhenzhen</given-names></name><xref ref-type="aff" rid="aff-2">2</xref></contrib>
<contrib id="author-4" contrib-type="author" corresp="yes">
<name name-style="western"><surname>Dong</surname><given-names>Haiying</given-names></name><xref ref-type="aff" rid="aff-1">1</xref><email>hyDong@mail.lzjtu.cn</email></contrib>
<aff id="aff-1"><label>1</label><institution>School of New Energy and Power Engineering, Lanzhou Jiaotong University</institution>, <addr-line>Lanzhou, 730070</addr-line>, <country>China</country></aff>
<aff id="aff-2"><label>2</label><institution>State Grid Gansu Electric Power Research Institute</institution>, <addr-line>Lanzhou, 730070</addr-line>, <country>China</country></aff>
</contrib-group>
<author-notes>
<corresp id="cor1"><label>&#x002A;</label>Corresponding Author: Haiying Dong. Email: <email>hyDong@mail.lzjtu.cn</email></corresp>
</author-notes>
<pub-date date-type="collection" publication-format="electronic">
<year>2023</year></pub-date>
<pub-date date-type="pub" publication-format="electronic"><day>1</day><month>4</month><year>2023</year>
</pub-date>
<volume>120</volume>
<issue>6</issue>
<fpage>1401</fpage>
<lpage>1420</lpage>
<history>
<date date-type="received">
<day>26</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>12</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; 2023 Huang et al.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Huang et al.</copyright-holder>
<license xlink:href="https://creativecommons.org/licenses/by/4.0/">
<license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="TSP_EE_27365.pdf"></self-uri>
<abstract>
<p>The virtual synchronous generator (VSG) can simulate synchronous machine&#x2019;s operation mechanism in the control link of an energy storage converter, so that an electrochemical energy storage power station has the ability to actively support the power grid, from passive regulation to active support. Since energy storage is an important physical basis for realizing the inertia and damping characteristics in VSG control, energy storage constraints of the physical characteristics on the system control parameters are analyzed to provide a basis for the system parameter tuning. In a classic VSG control, its virtual inertia and damping coefficient remain unchanged. When the grid load changes greatly, the constant control strategy most likely result in the grid frequency deviation beyond the stable operation standard limitations. To solve this problem, a comprehensive control strategy considering electrified wire netting demand and energy storage unit state of charge (SOC) is proposed, and an adaptive optimization method of VSG parameters under different SOC is given. The energy storage battery can maintain a safe working state at any time and be smoothly disconnected, which can effectively improve the output frequency performance of energy storage system. Simulation results further demonstrated the effectiveness of the VSG control theoretical analysis.</p>
</abstract>
<kwd-group kwd-group-type="author">
<kwd>VSG</kwd>
<kwd>energy storage power station</kwd>
<kwd>physical constraints of energy storage</kwd>
<kwd>adaptive parameter</kwd>
<kwd>frequency performance</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<label>1</label>
<title>Introduction</title>
<p>Since the goals of carbon peaking and carbon neutrality and new-energy-oriented electric systems have been put forward, the power industry has established a new power system that treats renewable energy as the principal part of the reform goal in China. Since renewable energy&#x2019;s output is uncertain, the change of the power system&#x2019;s main source will further reduce its stability. The introduction of energy storage units into the power system shall not only improve the utilization rate of new energy but also actively carry out peak shaving and frequency modulation [<xref ref-type="bibr" rid="ref-1">1</xref>&#x2013;<xref ref-type="bibr" rid="ref-4">4</xref>].</p>
<p>With the innovation of battery technology, large-capacity centralized energy storage power stations continue to be used as power sources to provide energy support for the grid [<xref ref-type="bibr" rid="ref-5">5</xref>&#x2013;<xref ref-type="bibr" rid="ref-7">7</xref>], which are included in the grid-connected operation and auxiliary service management. Li et al. [<xref ref-type="bibr" rid="ref-8">8</xref>,<xref ref-type="bibr" rid="ref-9">9</xref>] concluded that the main functions of the energy storage power station are peak load regulation, long-term power supply, primary frequency regulation, stabilizing power fluctuation, standby power and tracking planned power generation.</p>
<p>According to Zhong [<xref ref-type="bibr" rid="ref-10">10</xref>], nowadays distributed generation system with the power electronic converter as the interface lacks the inertia and damping of the traditional synchronous generator (SG). In the wake of increasing new energy penetration rate and withdrawal of corresponding SGs, the system inertia level and primary frequency modulation resources gradually decrease. Under the disturbance of the same intensity, Delille et al. [<xref ref-type="bibr" rid="ref-11">11</xref>] had concluded that the reduction of the inertia level will worsen the maximum frequency offset and the maximum frequency change rate, so the power system is more vulnerable to power fluctuations and system failures.</p>
<p>Over the last decade, Zhong et al. [<xref ref-type="bibr" rid="ref-12">12</xref>,<xref ref-type="bibr" rid="ref-13">13</xref>] proposed a virtual synchronous generator (VSG), which gives power electronic converter of energy storage power station capacity to sustain inertia and damping of the electrified wire netting by imitating SG, and enhance its anti-interference ability, give a pledge to electrical grids&#x2019; safe and steady operation. In the distributed system described by Torres et al. [<xref ref-type="bibr" rid="ref-14">14</xref>], the converter takes an adaptive control of VSG, and relevant adaptive control rules are designed to optimize the inertia and damping coefficients to reduce the output power of energy storage devices. In addition, the anti-jamming capability of the system can also get improved by using a newly designed equation about frequency-dependent coefficient [<xref ref-type="bibr" rid="ref-15">15</xref>,<xref ref-type="bibr" rid="ref-16">16</xref>]. Such expositions are unsatisfactory because the above adaptive VSG control strategies need to accurately model the converter to get the best control.</p>
<p>On the other hand, fuzzy control uses a language method, which is easy to carry out human-computer interaction. The precise mathematical model does not need to be established in the control process. It has strong robustness and is suitable for solving the problems of nonlinearity, strong coupling, time-varying and delay in process control [<xref ref-type="bibr" rid="ref-17">17</xref>,<xref ref-type="bibr" rid="ref-18">18</xref>]. Ebrahimi et al. [<xref ref-type="bibr" rid="ref-19">19</xref>&#x2013;<xref ref-type="bibr" rid="ref-23">23</xref>] used fuzzy adaptive VSG control to reduce frequency fluctuation.</p>
<p>Hu et al. [<xref ref-type="bibr" rid="ref-24">24</xref>] expound the essence of inertia in an energy storage system with the theory that inertia is a form of energy, and an equivalent inertia calculation method for a given energy storage system is came up. On this basis, Yuan et al. [<xref ref-type="bibr" rid="ref-25">25</xref>] established the physical constraints of the energy storage side on the virtual synchronous machine plus its operating boundary.</p>
<p>Chen et al. [<xref ref-type="bibr" rid="ref-26">26</xref>] proposed an adaptive control strategy to achieve SOC balance between several VSGs based on paralleled energy storage systems in a decentralized manner. With the proposed method, no matter what the difference between VSGs&#x2019; initial values and interference is, the SOC balance can be easily achieved without a communication network. The above research fails to consider that when the energy storage system shuts down at the minimum SOC, the output power of VSG will be greatly reduced in a short time, which will seriously affect the power grid.</p>
<p>To sum up, a VSG adaptive control strategy considering SOC feedback of energy storage battery is proposed. Firstly, this paper lists physical constraints based on the energy storage side, thus given VSG&#x2019;s inertia, output power and change rate of the power have corresponding value ranges. Then divides the VSG parameter settings of inertia, damping, active and reactive reference values into three groups. Different control strategies are adopted under different SOC to realize parameter adaptive control, so as to smooth power fluctuations when VSG is off-grid while maintaining a certain frequency modulation performance. Finally, through the MATLAB platform, the validity of the above integrated control strategy is verified by simulation.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>VSG Control of Energy Storage Power Station</title>
<p>Based on the classical second-order model of SG with polar logarithm of 1, the second-order VSG model is established. In order to meet the seamless switching and smooth transition between isolated island operation and grid connected operation of energy storage power station, a voltage source VSG mode is adopted in this paper, where <xref ref-type="fig" rid="fig-1">Fig. 1</xref> is its topology diagram, where mechanical and electromagnetic equations can broadly be defined as follows:</p>
<p><disp-formula id="eqn-1"><label>(1)</label><mml:math id="mml-eqn-1" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>J</mml:mi><mml:mi>&#x03C9;</mml:mi><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mrow><mml:mtext>d</mml:mtext></mml:mrow><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>+</mml:mo><mml:mi>D</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>i</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msub><mml:mi>e</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>i</mml:mi><mml:mrow><mml:mi>a</mml:mi><mml:mi>b</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<fig id="fig-1">
<label>Figure 1</label>
<caption>
<title>Configuration of the VSG</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-1.tif"/>
</fig>
<p>In <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref>: <italic>P</italic><sub>m</sub> is the virtual machine power; <italic>P</italic><sub>e</sub> is the output active power; <italic>J</italic> is virtual inertia; <inline-formula id="ieqn-1"><mml:math id="mml-ieqn-1"><mml:mi mathvariant="bold-italic">&#x03C9;</mml:mi></mml:math></inline-formula> is the real-time electric angular velocity generated by VSG; <italic>D</italic> is the damping coefficient; <inline-formula id="ieqn-2"><mml:math id="mml-ieqn-2"><mml:msub><mml:mi mathvariant="bold-italic">&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> is the rated angular velocity of the power grid. Where in <xref ref-type="fig" rid="fig-1">Fig. 1</xref>, <italic>U</italic><sub>dc</sub> is the DC side voltage; <italic>e</italic><sub>abc</sub>, <italic>u</italic><sub>abc</sub> and <italic>i</italic><sub>abc</sub> are the potential, terminal voltage and output current of VSG, respectively; <italic>R</italic><sub>f</sub>, <italic>L</italic><sub>f</sub> and <italic>C</italic><sub>f</sub> are the filter resistance, filter inductance and filter capacitance of AC side, respectively. <italic>R</italic><sub>g</sub>, and <italic>L</italic><sub>g</sub> are the equivalent resistance and inductance of power grid.</p>
<p>In this paper, the VSG control is primarily composed of active frequency control and reactive voltage control. <italic>P</italic><sub>m</sub> is made up of active power given reference value <italic>P</italic><sub>ref</sub> and frequency deviation response command value:</p>
<p><disp-formula id="eqn-2"><label>(2)</label><mml:math id="mml-eqn-2" display="block"><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula>where <italic>K</italic><sub>f</sub> is the active frequency modulation coefficient of VSG.</p>
<p>VSG&#x2019;s transient potential <italic>E</italic> consists of three parts: <italic>E</italic><sub>0</sub>, the no-load potential, represents the machine terminal voltage when the converter runs off-grid without load; the machine terminal voltage regulation amount and the reactive power regulation amount. In summary, the given output voltage amplitude expression is:</p>
<p><disp-formula id="eqn-3"><label>(3)</label><mml:math id="mml-eqn-3" display="block"><mml:mi>E</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>v</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>U</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>q</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mrow><mml:mi>r</mml:mi><mml:mi>e</mml:mi><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>Q</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>In <xref ref-type="disp-formula" rid="eqn-3">Eq. (3)</xref>, <italic>K</italic><sub>v</sub> is the voltage regulation coefficient of the VSG; <italic>U</italic><sub>ref</sub> means the rated value of the machine terminal voltage; <italic>U</italic> is the converter output port line voltage; <italic>K</italic><sub>q</sub> represents the reactive power regulation coefficient of the VSG; <italic>Q</italic><sub>ref</sub> indicates system&#x2019;s reactive power command, <italic>Q</italic><sub>e</sub> indicates instantaneous reactive power output of machine terminal. From <xref ref-type="disp-formula" rid="eqn-1">Eqs. (1)</xref>&#x2013;<xref ref-type="disp-formula" rid="eqn-3">(3)</xref>, the VSG control scheme adopted in this paper is shown in <xref ref-type="fig" rid="fig-2">Fig. 2</xref>.</p>
<fig id="fig-2">
<label>Figure 2</label>
<caption>
<title>Configuration of VSG control</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-2.tif"/>
</fig>
<p>The output modulation voltage is controlled by an inner loop of inductor current and an outer loop of capacitor voltage, and the output 6-pulse PWM signal is modulated by Space Vector Pulse Width Modulation (SVPWM), which is shown in <xref ref-type="fig" rid="fig-3">Fig. 3</xref>. In the figure, <italic>L</italic><sub>v</sub> represents virtual impedance.</p>
<fig id="fig-3">
<label>Figure 3</label>
<caption>
<title>Voltage and current double closed-loop control</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-3.tif"/>
</fig>
</sec>
<sec id="s3">
<label>3</label>
<title>Physical Constraints for Safe and Stable Operation</title>
<p>In practical energy storage projects, the storage capacity of VSG is usually fixed, and its size is limited by factors such as environment and investment. Therefore, to ensure that the inertia and damping of the VSG are always subject to the energy storage during the setting process and the corresponding response characteristics of the system during the transient process, the corresponding physical constraints index of the VSG are designed, namely, SOC, the charge rate and discharge rate, and power change rate.</p>
<sec id="s3_1">
<label>3.1</label>
<title>Working Characteristics of Storage Batteries</title>
<p>Generally, the residual power of the energy storage unit is represented by SOC. The maximum capacity of a fully charged battery is <italic>S</italic><sub>CN</sub>, define <italic>SOC</italic><sub>0</sub> as the initial SOC when the energy storage converter starts to work, then the instantaneous charging and discharging current is defined as <italic>i</italic>(t), the expression of SOC is as follows:</p>
<p><disp-formula id="eqn-4"><label>(4)</label><mml:math id="mml-eqn-4" display="block"><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>=</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mfrac><mml:mrow><mml:mo>&#x222B;</mml:mo><mml:mi>i</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula></p>
<p><disp-formula id="eqn-5"><label>(5)</label><mml:math id="mml-eqn-5" display="block"><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub><mml:mo>&#x2264;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>&#x2264;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p><italic>V</italic><sub>0</sub> is defined as the initial working voltage of the battery, and the initial state of charge <italic>SOC</italic><sub>0</sub> is as follows:</p>
<p><disp-formula id="eqn-6"><label>(6)</label><mml:math id="mml-eqn-6" display="block"><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:msubsup><mml:mi>V</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mfrac></mml:math></disp-formula></p>
<p>In practical terms, a battery has a certain normal working range, as shown in <xref ref-type="disp-formula" rid="eqn-5">Eq. (5)</xref>. In order to avoid deep discharge and overcharge that accelerate the battery&#x2019;s life loss, when the SOC of the battery is less than <italic>SOC</italic><sub>min</sub>, the minimum value of SOC, the battery shall only be charged; when the SOC of the battery is greater than <italic>SOC</italic><sub>max</sub>, the maximum value of SOC, the battery should not be overcharged. As shown in <xref ref-type="fig" rid="fig-4">Fig. 4</xref>, when the lithium iron phosphate battery is discharged under the same temperature and constant current, the characteristic curve between the normal working regions can be approximately regarded as one straight line whose slope value is <italic>k</italic>, define the charge cut-off voltage as <italic>V</italic><sub>max</sub> when SOC &#x003D; <italic>SOC</italic><sub>max</sub>, and define the discharge cut-off voltage as <italic>V</italic><sub>min</sub> when SOC &#x003D; <italic>SOC</italic><sub>min</sub>.</p>
<fig id="fig-4">
<label>Figure 4</label>
<caption>
<title>Lithium iron phosphate battery (50 Ah) discharge curve</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-4.tif"/>
</fig>
<p>The evaluation method for the battery&#x2019;s charging or discharging rate is shown in the following <xref ref-type="disp-formula" rid="eqn-7">Eq. (7)</xref>, and its unit is <italic>C</italic>.</p>
<p><disp-formula id="eqn-7"><label>(7)</label><mml:math id="mml-eqn-7" display="block"><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>h</mml:mi><mml:mo>,</mml:mo><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>i</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula></p>
<p>The charging and discharging rate is a characteristic quantity to describe the energy storage battery&#x2019;s charging and discharging speed. During constant current discharge, the electric energy output by the battery with different discharge rates is different. The greater the current, the smaller the output capacity. The battery&#x2019;s maximum charge (discharge) rate is related to the physical properties of its own electrolyte. Excessive charge (discharge) rate will shorten the design working life of energy storage units. For better protection, the charging and discharging rates are set with maximum values <inline-formula id="ieqn-3"><mml:math id="mml-ieqn-3"><mml:msubsup><mml:mi mathvariant="bold-italic">&#x03B1;</mml:mi><mml:mrow><mml:mi mathvariant="bold-italic">c</mml:mi><mml:mi mathvariant="bold-italic">h</mml:mi></mml:mrow><mml:mrow><mml:mo mathvariant="bold" movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-4"><mml:math id="mml-ieqn-4"><mml:msubsup><mml:mi mathvariant="bold-italic">&#x03B1;</mml:mi><mml:mrow><mml:mi mathvariant="bold-italic">d</mml:mi><mml:mi mathvariant="bold-italic">i</mml:mi><mml:mi mathvariant="bold-italic">s</mml:mi></mml:mrow><mml:mrow><mml:mo mathvariant="bold" movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>.</p>
<p>In addition, the charging and discharging power <italic>P</italic><sub>ch</sub> and <italic>P</italic><sub>dis</sub> are respectively used to indicate the charging and discharging capacity that the battery can generate:</p>
<p><disp-formula id="eqn-8"><label>(8)</label><mml:math id="mml-eqn-8" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>k</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:msubsup><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:mtext>ch</mml:mtext></mml:mrow></mml:mrow><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup><mml:mo>=</mml:mo><mml:mo>&#x2212;</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>k</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:msubsup><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>It should be noted that <italic>V</italic><sub>min</sub> and <italic>V</italic><sub>max</sub> are the cut-off voltage of the battery when working at the maximum charge (discharge) rate.</p>
<p>The power change rate is an important factor that affects the dynamic response performance of VSG system or the life of the energy storage device. When the load power of the power grid changes, the VSG shall compensate for the fluctuation, therefore its output power changes accordingly. Especially when the power change rate is too big, it will lead to the chemical reaction inside the battery fail to keep up, leading to partial overcharge or partial overdischarge inside the battery.</p>
<p><disp-formula id="eqn-9"><label>(9)</label><mml:math id="mml-eqn-9" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mrow><mml:mo>|</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x2264;</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mi>&#x03B2;</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>h</mml:mi></mml:mrow></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>&#x2264;</mml:mo><mml:msub><mml:mi>&#x03B2;</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>&#x03B2;</italic><sub>ch</sub> and <italic>&#x03B2;</italic><sub>dis</sub> represent the limit values of the change rate of battery charge, discharge power. To avoid accelerated degradation of battery performance under severe dynamic loading conditions, the battery power change rate is usually limited to 2%&#x007E;20% of the maximum output power. As for the whole energy storage power station, its power limit and power change rate are limited to the installed capacity, and the maximum limit of 1-min active power change is 10% of the installed capacity.</p>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Physical Constraints of Energy Storage Side</title>
<p>As the energy supply end of the energy storage power station, the battery&#x2019;s performance finally determines the frequency response capability of VSG. The energy storage battery must always be kept in a normal working state so as to prolong its service life and avoid fire accidents caused by overheating. The theory limits the output power from SOC, discharge power, and change rate of discharge power, respectively.</p>
<p>Battery&#x2019;s SOC constraint: According to <xref ref-type="disp-formula" rid="eqn-5">Eq. (5)</xref>, when the battery is not within the normal operating range, the energy storage converter shall actively disconnect the corresponding batteries to ensure that those batteries will not discharge excessively.</p>
<p>Battery&#x2019;s discharge power constraint: Due to the limitation of the electrochemical reaction, usually battery&#x2019;s discharge power is supposed to be limited within a suitable range, which is shown in <xref ref-type="disp-formula" rid="eqn-10">Eq. (10)</xref>. When the discharge power of the energy storage converter exceeds this limit, the VSG shall make <italic>P</italic><sub>m</sub> equal to the rated discharge power of the battery.</p>
<p><disp-formula id="eqn-10"><label>(10)</label><mml:math id="mml-eqn-10" display="block"><mml:msub><mml:mrow><mml:mtext>P</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mtext>m</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>&#x2A7D;</mml:mo><mml:msubsup><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:mtext>dis</mml:mtext></mml:mrow></mml:mrow><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup></mml:math></disp-formula></p>
<p>Battery&#x2019;s change rate of discharge power constraint: as shown in <xref ref-type="disp-formula" rid="eqn-11">Eq. (11)</xref>. When the discharge power change rate exceeds this limit, the energy storage converter shall reduce the change amount of <italic>P</italic><sub>m</sub> until the constraints conditions are met.</p>
<p><disp-formula id="eqn-11"><label>(11)</label><mml:math id="mml-eqn-11" display="block"><mml:mrow><mml:mi mathvariant="normal">&#x0394;</mml:mi></mml:mrow><mml:msub><mml:mrow><mml:mtext>P</mml:mtext></mml:mrow><mml:mrow><mml:mrow><mml:mtext>m</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>&#x2A7D;</mml:mo><mml:msub><mml:mi>&#x03B2;</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:math></disp-formula></p>
<p>Virtually, the inertia and damping characteristics of VSG are controlled by energy storage units, so the setting of virtual inertia <italic>J</italic> and damping coefficient <italic>D</italic> should also take the constraints of the discharge characteristics of the battery into consideration. We need hardly say that the inertia of the SG is derived from the mechanical energy in the rotor. Likewise, the inertia of the VSG is actually derived from the potential energy in the battery. In <xref ref-type="disp-formula" rid="eqn-12">Eq. (12)</xref>, <italic>W</italic><sub>k</sub> means the rotor rotational kinetic energy of SG; <italic>E</italic><sub>vsg</sub> is the equivalent rotor kinetic energy of VSG; <italic>E</italic><sub>C</sub> indicates the electric potential energy of the battery.</p>
<p><disp-formula id="eqn-12"><label>(12)</label><mml:math id="mml-eqn-12" display="block"><mml:msub><mml:mi>W</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mi>J</mml:mi><mml:msup><mml:mrow><mml:mi mathvariant="normal">&#x03A9;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mrow><mml:mtext>VSG</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mi>J</mml:mi><mml:msup><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mrow><mml:mtext>C</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mi>C</mml:mi><mml:msup><mml:mi>u</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></disp-formula></p>
<p>While the battery is discharging, it can not be fully discharged. Most part of its potential energy is used to support VSG normal output and frequency modulation, and the rest part is used to maintain the SOC of the battery. Therefore, when the law of conservation of energy as <xref ref-type="disp-formula" rid="eqn-12">(12)</xref> shows and, hence, the energy used to maintain the lowest SOC needs to be deducted from the potential energy, as shown in <xref ref-type="disp-formula" rid="eqn-13">Eq. (13)</xref>.</p>
<p><disp-formula id="eqn-13"><label>(13)</label><mml:math id="mml-eqn-13" display="block"><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mi>J</mml:mi><mml:msup><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:mn>2</mml:mn></mml:mfrac><mml:mi>C</mml:mi><mml:msup><mml:mi>u</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math></disp-formula></p>
<p>Then, the next step is to obtain the battery&#x2019;s equivalent inertia <italic>J</italic><sub>eq</sub> considering constraints of SOC and maximum discharge rate:</p>
<p><disp-formula id="eqn-14"><label>(14)</label><mml:math id="mml-eqn-14" display="block"><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub><mml:mfrac><mml:mrow><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:msup><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mfrac><mml:msup><mml:mi>u</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>&#x2264;</mml:mo><mml:msubsup><mml:mi>&#x03B1;</mml:mi><mml:mrow><mml:msub><mml:mi></mml:mi><mml:mrow><mml:mi>d</mml:mi><mml:mi>i</mml:mi><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo movablelimits="true" form="prefix">max</mml:mo></mml:mrow></mml:msubsup><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>q</mml:mi></mml:mrow></mml:msub><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mi>C</mml:mi><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mfrac><mml:mrow><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:mi>C</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>S</mml:mi><mml:mi>O</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mo movablelimits="true" form="prefix">min</mml:mo></mml:mrow></mml:msub></mml:mrow><mml:msubsup><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mfrac></mml:math></disp-formula>where <italic>C</italic><sub>eq</sub> is defined as the equivalent capacitance value of the battery under the rated conditions according to <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref>. In <xref ref-type="disp-formula" rid="eqn-15">Eq. (15)</xref>, <italic>U</italic><sub>CN</sub> is the battery&#x2019;s rated voltage, <italic>I</italic><sub>CN</sub> is its rated current, and <italic>T</italic><sub>CN</sub> is its rated discharge time.</p>
<p><disp-formula id="eqn-15"><label>(15)</label><mml:math id="mml-eqn-15" display="block"><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mrow><mml:mtext>eq</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:msubsup><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>2</mml:mn><mml:msub><mml:mi>I</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mrow><mml:mtext>CN</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mfrac></mml:math></disp-formula></p>
</sec>
</sec>
<sec id="s4">
<label>4</label>
<title>Improved VSG Control Strategy Based on SOC Feedback</title>
<p>During the normal shutdown of the conventional generator set, the load shall be gradually removed, the load switch shall be disconnected, and the generator shall be shut down after running at idle speed for 3&#x2013;5 min. Since the capacity of the energy storage batteries is still very limited compared to the large grid, the energy storage power station will temporarily stop providing power grid support and frequency modulation services to prevent its energy storage units from discharging beyond the limit especially after one battery pack have been discharging to the power grid for a long period of time. Due to the discharge characteristics of the electric energy storage system, once the energy storage battery is disconnected from the energy storage converter, the energy storage converter AC output shall drop sharply, along with grid frequency becoming unstable, and VSG frequency modulation completely loses its energy support, which is a huge disturbance to the grid.</p>
<sec id="s4_1">
<label>4.1</label>
<title>Improved VSG Control Structure</title>
<p>When the SOC exceeds the energy storage battery&#x2019;s specified reasonable range, the battery must exit the frequency modulation process and use other battery packs to adjust the system frequency. It can be considered to appropriately reduce the battery output when the SOC is really low (still in the normal range) to form a variable VSG regulation strategy based on SOC, as shown in <xref ref-type="fig" rid="fig-5">Fig. 5a</xref>.</p>
<fig id="fig-5">
<label>Figure 5</label>
<caption>
<title>Improved VSG control structure</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-5.tif"/>
</fig>
<p><italic>P</italic><sub><italic>ref</italic>0</sub> and <italic>Q</italic><sub>ref0</sub> are the active and reactive output power commands received by the VSG during grid connection operation. These commands can be automatically generated according to the power shortage obtained by the computer or can be manually specified. <italic>P</italic><sub>ref</sub> and <italic>Q</italic><sub>ref</sub> are power commands given to VSG.</p>
<p>On the one hand, it can make the energy storage battery better maintain the SOC, slow down its performance degradation and prolong its service life. On the other hand, it can effectively avoid the adverse impact on the stability of the power grid when the SOC exceeds the limit. A schematic diagram of SOC partition of the energy storage battery is shown in <xref ref-type="fig" rid="fig-5">Fig. 5b</xref>.</p>
<p>In order to reduce the negative impact on the power grid when the energy storage power station is under power shortage and actively cuts off the network, it is necessary to imitate the normal shutdown working principle of the conventional generator set, adopt the buffer control strategy when the battery is at low SOC, actively reduce the power output of the energy storage system step by step, and select new inertia and damping parameters as shown in <xref ref-type="table" rid="table-1">Table 1</xref>, to ensure that the frequency fluctuations are still within the normal range during the gradual decline of the output. Meanwhile, this control strategy provides sufficient response time for other power supply frequency modulation systems.</p>
<table-wrap id="table-1">
<label>Table 1</label>
<caption>
<title>Adaptive optimization scheme of variable VSG control</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>State of charge</th>
<th align="center" colspan="2">Control parameters</th>
<th align="center" colspan="2">Reference value</th>
<th>Control method</th>
</tr>
</thead>
<tbody>
<tr>
<td>SOC<sub>b1</sub>-SOC<sub>max</sub></td>
<td><italic>J</italic> &#x003D; <italic>J</italic><sub>0</sub> &#x002B; <italic>J</italic>a</td>
<td><italic>D</italic> &#x003D; <italic>D</italic><sub>0</sub> &#x002B; <italic>D</italic>a</td>
<td><italic>P</italic><sub>ref</sub></td>
<td><italic>Q</italic><sub>ref</sub></td>
<td>Fuzzy control</td>
</tr>
<tr>
<td>SOC<sub>b2</sub>-SOC<sub>b1</sub></td>
<td><italic>J</italic> &#x003D; <italic>J</italic><sub>1</sub></td>
<td><italic>D</italic> &#x003D; <italic>D</italic><sub>1</sub></td>
<td>0.7<italic>P</italic><sub>ref</sub></td>
<td>0.7<italic>Q</italic><sub>ref</sub></td>
<td>Constant control</td>
</tr>
<tr>
<td>SOC<sub>min</sub>-SOC<sub>b2</sub></td>
<td><italic>J</italic> &#x003D; <italic>J</italic><sub>2</sub></td>
<td><italic>D</italic> &#x003D; <italic>D</italic><sub>2</sub></td>
<td>0.3<italic>P</italic><sub>ref</sub></td>
<td>0.3<italic>Q</italic><sub>ref</sub></td>
<td>Constant control</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Parameters J and D Influence on Frequency Characteristics</title>
<p>According to the derivation of the expression of virtual synchronous machine output active power and reactive power in reference [<xref ref-type="bibr" rid="ref-25">25</xref>], it can be concluded that:</p>
<p><disp-formula id="eqn-16"><label>(16)</label><mml:math id="mml-eqn-16" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:mtext>e</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mn>1</mml:mn><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:mstyle><mml:mrow><mml:mo>[</mml:mo><mml:mi>E</mml:mi><mml:mi>U</mml:mi><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B4;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msup><mml:mi>U</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>Q</mml:mi><mml:mrow><mml:mrow><mml:mtext>e</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mn>1</mml:mn><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:mstyle><mml:mrow><mml:mo>[</mml:mo><mml:mi>E</mml:mi><mml:mi>U</mml:mi><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>&#x03B4;</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>&#x2212;</mml:mo><mml:msup><mml:mi>U</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>Z</italic><sub>f</sub>, <italic>&#x03B8;</italic>, <italic>&#x03B4;</italic> are respectively the output filter impedance, impedance angle and output voltage phase of VSG, which can be expressed as:</p>
<p><disp-formula id="eqn-17"><label>(17)</label><mml:math id="mml-eqn-17" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:msup><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:msqrt></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mi>&#x03B8;</mml:mi><mml:mo>=</mml:mo><mml:mi>arctan</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>&#x03C9;</mml:mi><mml:msub><mml:mi>L</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x03B4;</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mo>&#x222B;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mstyle></mml:mtd></mml:mtr></mml:mtable></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>Calculate the voltage amplitude <italic>E</italic><sub>s</sub> and voltage phase <italic>&#x03B4;</italic><sub>s</sub> corresponding to the stable operating point:</p>
<p><disp-formula id="eqn-18"><label>(18)</label><mml:math id="mml-eqn-18" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mrow><mml:mrow><mml:mtext>ref</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msup><mml:mi>U</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B8;</mml:mi></mml:mrow><mml:mrow><mml:mi>U</mml:mi><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>arctan</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msub><mml:mi>Q</mml:mi><mml:mrow><mml:mrow><mml:mtext>ref</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi>U</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B8;</mml:mi></mml:mrow><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:mstyle></mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mrow><mml:mtext>ref</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msup><mml:mi>U</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>cos</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mi>&#x03B8;</mml:mi></mml:mrow><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mfrac></mml:mstyle></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>where <italic>P</italic><sub>e</sub> and <italic>&#x03B4;</italic> satisfy the following relationship:</p>
<p><disp-formula id="eqn-19"><label>(19)</label><mml:math id="mml-eqn-19" display="block"><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mrow><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi></mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mrow><mml:mi mathvariant="normal">&#x2202;</mml:mi></mml:mrow><mml:mi>&#x03B4;</mml:mi></mml:mrow></mml:mfrac><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:mi>&#x03B4;</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>E</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:mrow><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mi>U</mml:mi></mml:mrow><mml:msub><mml:mi>Z</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:mfrac><mml:mi>sin</mml:mi><mml:mo>&#x2061;</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>&#x03B8;</mml:mi><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>&#x03B4;</mml:mi><mml:mrow><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>According to <xref ref-type="disp-formula" rid="eqn-1">Eq. (1)</xref>, taking small signal model analysis method with reference to SG, and obtaining second-order transfer function between VSG&#x2019;s active power output and input, as shown in <xref ref-type="disp-formula" rid="eqn-20">Eq. (20)</xref>. And by analyzing the above second order system, this paper&#x2019;s VSG natural angular frequency and damping ratio are available in the following equations:</p>
<p><disp-formula id="eqn-20"><label>(20)</label><mml:math id="mml-eqn-20" display="block"><mml:mi>G</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo>(</mml:mo><mml:mi>s</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi>J</mml:mi><mml:msup><mml:mi>s</mml:mi><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:mrow><mml:mo>(</mml:mo><mml:mi>D</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:math></disp-formula></p>
<p><disp-formula id="eqn-21"><label>(21)</label><mml:math id="mml-eqn-21" display="block"><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mrow><mml:mtext>n</mml:mtext></mml:mrow></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:msqrt><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mrow><mml:mtext>E</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mrow><mml:mi>J</mml:mi></mml:mfrac></mml:mstyle></mml:msqrt></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x03B6;</mml:mi><mml:mo>=</mml:mo><mml:mn>0.5</mml:mn><mml:mrow><mml:mo>(</mml:mo><mml:mi>D</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mstyle displaystyle="true" scriptlevel="0"><mml:msqrt><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:msub><mml:mi>&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mrow><mml:mi>J</mml:mi><mml:msub><mml:mi>S</mml:mi><mml:mrow><mml:mi>E</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:msqrt></mml:mstyle></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p>The initial inertia <italic>J</italic><sub>0</sub> of the VSG controller is designed according to the natural oscillation frequency of SG (0.628&#x007E;15.700 rad/s). After determining the oscillation angular frequency of VSG, consider the concept of optimal second-order system and make &#x03B6; &#x003D; 0.707 to obtain faster response speed and lower overshoot. And so <italic>D</italic><sub>0</sub> can be obtained from <xref ref-type="disp-formula" rid="eqn-21">Eq. (21)</xref>.</p>
<p>When the VSG&#x2019;s SOC is at SOC<sub>b1</sub>-SOC<sub>max</sub>, under damping control (0 &#x003C; &#x03B6; &#x003C; 1) is adopted to ensure the energy storage system&#x2019;s dynamic and steady-state characteristics; When its SOC is at SOC<sub>min</sub>-SOC<sub>b1</sub>, to realize a smooth transition of output power, over damping control (&#x03B6; &#x003E; 1) is adopted.</p>
<p>The VSG&#x2019;s frequency response characteristics are analyzed as shown in <xref ref-type="fig" rid="fig-6">Fig. 6</xref>. The findings from the above studies suggest that when <italic>D</italic> is constant, the system&#x2019;s natural oscillation angular frequency &#x03C9;<italic>n</italic> and its damping ratio &#x03B6; gradually decreases with the increase of <italic>J</italic>. When <italic>J</italic> is constant and <italic>D</italic> grows, the &#x03C9;<italic>n</italic> remains unchanged, while &#x03B6; increases gradually. And as shown in <xref ref-type="fig" rid="fig-6">Fig. 6a</xref>, VSG&#x2019;s overshoot and rise time decrease when <italic>J</italic> increases, meanwhile its adjustment time increases. Therefore, the increase of <italic>J</italic> will make VSG&#x2019;s oscillation becomes more intense and slow down the system&#x2019;s response speed. What is shown in <xref ref-type="fig" rid="fig-6">Fig. 6b</xref>, VSG&#x2019;s overshoot, rise time and adjustment time decrease, but the peak time will move forward. So as <italic>D</italic> increases, the VSG frequency response curve is more stable. After summarizing and analyzing the various behaviors of frequency under different VSG parameters, the following fuzzy controller can be designed.</p>
<fig id="fig-6">
<label>Figure 6</label>
<caption>
<title>Influence of J and D with frequency modulation performance</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-6.tif"/>
</fig>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Fuzzy Controller Design</title>
<p>To achieve adaptive adjustment of virtual inertia and damping coefficient, this paper adopts a fuzzy controller.</p>
<p>Fuzzy adaptive control brings the system better robustness and adaptability. The nonlinear advantage of fuzzy control is that it does not need to build a complex mathematical model in the VSG control, and it has flexible adaptability. As <xref ref-type="fig" rid="fig-7">Fig. 7</xref> shows, this fuzzy controller consists of fuzzification, fuzzy reasoning and defuzzification. Set VSG&#x2019;s frequency deviation <inline-formula id="ieqn-5"><mml:math id="mml-ieqn-5"><mml:mi mathvariant="bold-italic">&#x03C9;</mml:mi></mml:math></inline-formula>-<inline-formula id="ieqn-6"><mml:math id="mml-ieqn-6"><mml:msub><mml:mi mathvariant="bold-italic">&#x03C9;</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> and its change rate d<inline-formula id="ieqn-7"><mml:math id="mml-ieqn-7"><mml:mi mathvariant="bold-italic">&#x03C9;</mml:mi></mml:math></inline-formula>/dt as controller&#x2019;s input, set inertia change value <italic>J</italic><sub>a</sub> and damping change value <italic>D</italic><sub>a</sub> as controller&#x2019;s output. Where <italic>k</italic><sub>e</sub> and <italic>k</italic><sub>ec</sub> are quantization factors, <italic>k</italic><sub>uJ</sub> and <italic>k</italic><sub>uD</sub> are scale coefficients. Set <italic>k</italic><sub>e</sub> &#x003D; 0.8, <italic>k</italic><sub>ec</sub> &#x003D; 0.015, <italic>k</italic><sub>uJ</sub> &#x003D; 0.2, <italic>k</italic><sub>uD</sub> &#x003D; 20.</p>
<fig id="fig-7">
<label>Figure 7</label>
<caption>
<title>Fuzzy controller structure diagram</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-7.tif"/>
</fig>
<p>The fuzzy sets that define the above four are {negative big, negative small, zero, positive small, positive big}, noted as {NB, NS, Z, PS, PB}. Set appropriate membership functions to limit the value range of <italic>J</italic> and <italic>D</italic> to avoid system power and change rate exceeding the limit, which is shown in <xref ref-type="fig" rid="fig-8">Fig. 8</xref>.</p>
<fig id="fig-8">
<label>Figure 8</label>
<caption>
<title>Membership functions of inputs and outputs</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-8.tif"/>
</fig>
<p>The fuzzy controller uses Mamdani fuzzy inference algorithm, the &#x201C;minimum-maximum&#x201D; synthesis rule, and the centroid defuzzification method. Specific fuzzy control rules adopted in this paper are shown in <xref ref-type="table" rid="table-2">Tables 2</xref> and <xref ref-type="table" rid="table-3">3</xref>.</p>
<table-wrap id="table-2">
<label>Table 2</label>
<caption>
<title><italic>J</italic><sub>a</sub> Fuzzy rules</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th rowspan="2">d&#x003C9;/dt</th>
<th align="center" colspan="5">&#x003C9;-&#x003C9;<sub>0</sub></th>
</tr>
<tr>
<th>NB</th>
<th>NS</th>
<th>Z</th>
<th>PS</th>
<th>PB</th>
</tr>
</thead>
<tbody>
<tr>
<td><bold>NB</bold></td>
<td>PB</td>
<td>PS</td>
<td>PS</td>
<td>NS</td>
<td>NB</td>
</tr>
<tr>
<td><bold>NS</bold></td>
<td>PS</td>
<td>PS</td>
<td>Z</td>
<td>NS</td>
<td>NS</td>
</tr>
<tr>
<td><bold>Z</bold></td>
<td>PS</td>
<td>Z</td>
<td>NS</td>
<td>Z</td>
<td>PS</td>
</tr>
<tr>
<td><bold>PS</bold></td>
<td>NS</td>
<td>NS</td>
<td>Z</td>
<td>PS</td>
<td>PS</td>
</tr>
<tr>
<td><bold>PB</bold></td>
<td>NB</td>
<td>NS</td>
<td>PS</td>
<td>PS</td>
<td>PB</td>
</tr>
</tbody>
</table>
</table-wrap><table-wrap id="table-3">
<label>Table 3</label>
<caption>
<title><italic>D</italic><sub>a</sub> Fuzzy rules</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th rowspan="2">d&#x003C9;/dt</th>
<th align="center" colspan="5">&#x003C9;-&#x003C9;<sub>0</sub></th>
</tr>
<tr>
<th>NB</th>
<th>NS</th>
<th>Z</th>
<th>PS</th>
<th>PB</th>
</tr>
</thead>
<tbody>
<tr>
<td><bold>NB</bold></td>
<td>PB</td>
<td>PB</td>
<td>PS</td>
<td>Z</td>
<td>PS</td>
</tr>
<tr>
<td><bold>NS</bold></td>
<td>PB</td>
<td>PS</td>
<td>Z</td>
<td>PS</td>
<td>PS</td>
</tr>
<tr>
<td><bold>Z</bold></td>
<td>PS</td>
<td>Z</td>
<td>Z</td>
<td>Z</td>
<td>PS</td>
</tr>
<tr>
<td><bold>PS</bold></td>
<td>PS</td>
<td>PS</td>
<td>Z</td>
<td>PS</td>
<td>PB</td>
</tr>
<tr>
<td><bold>PB</bold></td>
<td>PS</td>
<td>Z</td>
<td>PS</td>
<td>PB</td>
<td>PB</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Simulation Comparison and Analysis</title>
<p>To test the correctness and effectiveness of the theory and strategy mentioned above, three VSG control methods, Common VSG control (C-VSG), Conventional adaptive VSG control (S-VSG) and Improved VSG (I-VSG) control, are compared with each other. Note that the <italic>J</italic> and <italic>D</italic> of the C-VSG are constant. In addition, C-VSG does not consider the physical constraints on the energy storage side, and only power off and shutdown at <italic>SOC</italic><sub>min</sub>. The grid-connected simulation model of the energy storage power station as proposed in <xref ref-type="fig" rid="fig-1">Fig. 1</xref> has been constructed on the MATLAB platform.</p>
<p>The energy storage converter takes a VSG control method just the same way as <xref ref-type="fig" rid="fig-2">Fig. 2</xref>, and <xref ref-type="table" rid="table-4">Table 4</xref> shows the specific simulation parameters of VSG.</p>
<table-wrap id="table-4">
<label>Table 4</label>
<caption>
<title>System simulation parameters</title>
</caption>
<table frame="hsides">
<colgroup>
<col align="left"/>
<col align="left"/>
<col align="left"/>
</colgroup>
<thead>
<tr>
<th>Parameter</th>
<th>Symbol</th>
<th>Value</th>
</tr>
</thead>
<tbody>
<tr>
<td>Discharge voltage of energy storage battery</td>
<td><italic>U</italic><sub>dc</sub></td>
<td>800 V</td>
</tr>
<tr>
<td>Nominal voltage (phase to phase)</td>
<td><italic>u</italic><sub>AB</sub>, <italic>u</italic><sub>BC</sub>, <italic>u</italic><sub>CA</sub></td>
<td>380 V</td>
</tr>
<tr>
<td>Grid frequency</td>
<td><italic>f</italic><sub>N</sub></td>
<td>50 Hz</td>
</tr>
<tr>
<td>Rated active output power</td>
<td><italic>P</italic><sub>ref</sub></td>
<td>20 kW</td>
</tr>
<tr>
<td>Rated reactive output power</td>
<td><italic>Q</italic><sub>ref</sub></td>
<td>0 kW</td>
</tr>
<tr>
<td>Virtual inertia</td>
<td><italic>J</italic><sub>0</sub>, <italic>J</italic><sub>1</sub>, <italic>J</italic><sub>2</sub></td>
<td>(0.5, 0.3, 0.2) kg<bold>.</bold>m<sup>2</sup></td>
</tr>
<tr>
<td>Virtual damping</td>
<td><italic>D</italic><sub>0</sub>, <italic>D</italic><sub>1</sub>, <italic>D</italic><sub>2</sub></td>
<td>14, 25, 30</td>
</tr>
<tr>
<td>Switching frequency</td>
<td><italic>f</italic><sub>s</sub></td>
<td>5 kHz</td>
</tr>
<tr>
<td>Filter resistor</td>
<td><italic>R</italic><sub>f</sub></td>
<td>0.2 &#x03A9;</td>
</tr>
<tr>
<td>Filter inductance</td>
<td><italic>L</italic><sub>f</sub></td>
<td>5 &#x00D7; 10<sup>&#x2212;3</sup> H</td>
</tr>
<tr>
<td>Filter capacitor</td>
<td><italic>C</italic><sub>f</sub></td>
<td>2 &#x00D7; 10<sup>&#x2212;5</sup> F</td>
</tr>
<tr>
<td>Energy storage battery capacity</td>
<td><italic>S</italic><sub><italic>CN</italic></sub></td>
<td>1 A<bold>.</bold>h</td>
</tr>
<tr>
<td>State of charge of the battery</td>
<td><italic>SOC</italic><sub><italic>min</italic></sub>, <italic>SOC</italic><sub><italic>b1</italic></sub>,<break/><italic>SOC</italic><sub><italic>b2</italic></sub>, <italic>SOC</italic><sub><italic>max</italic></sub></td>
<td>5%, 6%, 7%, 95%</td>
</tr>
<tr>
<td>Virtual impedance</td>
<td><italic>L</italic><sub>v</sub></td>
<td>5 &#x00D7; 10<sup>&#x2212;3</sup> H</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>First, this paper verifies the influence of the discharge characteristics of the energy storage battery on the value range of <italic>J</italic> and <italic>D</italic>, as shown in <xref ref-type="fig" rid="fig-9">Fig. 9</xref>.</p>
<fig id="fig-9">
<label>Figure 9</label>
<caption>
<title>Simulations in case of different J and D</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-9.tif"/>
</fig>
<p>Set the system with a 5 kW, 1 kvar load. In this simulation experiment, the maximum discharge power limit of the energy storage battery is 12 kW. and the change threshold of discharge power is 80 kW/s. In the beginning, VSG operates with <italic>P</italic><sub>ref</sub> &#x003D; 4 kW and <italic>Q</italic><sub>ref</sub> &#x003D; 1 kvar. After 0.8 s, the active power command value of the system increases to 10 kW.</p>
<p>As shown in <xref ref-type="fig" rid="fig-9">Fig. 9</xref>, when the system power command changes, the system shows different response characteristics under different control parameters. The peak value of the active power change and its area surrounding the time axis is the required energy storage margin at this time.</p>

<p>In <xref ref-type="fig" rid="fig-9">Figs. 9a</xref>, <xref ref-type="fig" rid="fig-9">9c</xref>, when <italic>J</italic> &#x003D; 1.5, <italic>D</italic> &#x003D; 14 and <italic>J</italic> &#x003D; 0.5, <italic>D</italic> &#x003D; 3, the power fluctuation is large, the overshoot time is long, the maximum value exceeds the power limit, so the system response exceeds the capacity range of the energy storage unit.</p>

<p>As shown in <xref ref-type="fig" rid="fig-9">Fig. 9b</xref>, when <italic>J</italic> &#x003D; 0.1 and <italic>D</italic> &#x003D; 14, the maximum power change rate is80.164 kW/s, which is greater than the physical limit of the power change rate, exceeding its performance accessibility, so it cannot be achieved.</p>

<p>If <italic>J</italic> &#x003D; 0.5 and <italic>D</italic> &#x003D; 14, as shown in <xref ref-type="fig" rid="fig-9">Figs. 9a</xref>, <xref ref-type="fig" rid="fig-9">9b</xref>, when the active power command changes, its power change rate is up to 73.85 kW/s and is far from the limit value of the power change rate, so the system operates normally and safely. At this time, VSG is under damping control.</p>

<p>If <italic>J</italic> &#x003D; 0.5 and <italic>D</italic> &#x003D; 30, it can be seen that when VSG implements over-damping control, both power and power change rate meet operation constraints, and only the adjustment time is slightly extended. Therefore, the over-damping control can be adopted when the SOC is low. Since large system interference will not occur in a short time, the over-damping control with fixed parameters can maximize and weaken the adverse impact on the power grid when the output power drops.</p>
<p>Then, this paper verifies the function of the VSG fuzzy adaptive control module, as shown in <xref ref-type="fig" rid="fig-10">Fig. 10</xref>. In this comparative experiment, it is assumed that SOC is in the normal working zone, and compares three different VSG control modes.
<list list-type="order">
<list-item>
<p>Constant parameter VSG, C-VSG: <italic>J</italic> &#x003D; 0.5, <italic>D</italic> &#x003D; 14;</p>
</list-item>
<list-item>
<p>Conventional adaptive VSG, S-VSG: <italic>J</italic><sub>0</sub> &#x003D; 0.5, <italic>D</italic><sub>0</sub> &#x003D; 14;</p></list-item>
<list-item>
<p>Improved adaptive VSG (<xref ref-type="table" rid="table-1">Table 1</xref>), I-VSG: <italic>J</italic><sub>0</sub> &#x003D; 0.5, <italic>D</italic><sub>0</sub>&#x003D;14.</p>
</list-item>
</list></p>
<fig id="fig-10">
<label>Figure 10</label>
<caption>
<title>Load disturbance output comparison</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-10.tif"/>
</fig>
<p>The conventional adaptive control strategy as the control group is shown in <xref ref-type="disp-formula" rid="eqn-22">Eqs. (22)</xref> and <xref ref-type="disp-formula" rid="eqn-23">(23)</xref>:</p>
<p><disp-formula id="eqn-22"><label>(22)</label><mml:math id="mml-eqn-22" display="block"><mml:mi>J</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>&#x2264;</mml:mo><mml:mn>0</mml:mn></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mrow><mml:mo>|</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x2264;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>J</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>J</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>a</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>|</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi>&#x03C9;</mml:mi><mml:mo>&#x22C5;</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>&#x003E;</mml:mo><mml:mn>0</mml:mn><mml:mo>&#x2229;</mml:mo><mml:mrow><mml:mo>|</mml:mo><mml:mstyle displaystyle="true" scriptlevel="0"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi><mml:mi>&#x03C9;</mml:mi></mml:mrow><mml:mrow><mml:mi>d</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x003E;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mi>J</mml:mi></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula></p>
<p><disp-formula id="eqn-23"><label>(23)</label><mml:math id="mml-eqn-23" display="block"><mml:mi>D</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left left" rowspacing="4pt" columnspacing="1em"><mml:mtr><mml:mtd><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:mtd><mml:mtd><mml:mrow><mml:mo>|</mml:mo><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi>&#x03C9;</mml:mi><mml:mo>|</mml:mo></mml:mrow><mml:mo>&#x2264;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mrow><mml:mtext>D</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mi>D</mml:mi><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>b</mml:mi><mml:mrow><mml:mo>|</mml:mo><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi>&#x03C9;</mml:mi><mml:mo>|</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>|</mml:mo><mml:mi mathvariant="normal">&#x0394;</mml:mi><mml:mi>&#x03C9;</mml:mi><mml:mo>|</mml:mo></mml:mrow><mml:mo>&gt;</mml:mo><mml:msub><mml:mi>K</mml:mi><mml:mrow><mml:mrow><mml:mtext>D</mml:mtext></mml:mrow></mml:mrow></mml:msub></mml:mtd></mml:mtr></mml:mtable><mml:mo fence="true" stretchy="true" symmetric="true"></mml:mo></mml:mrow></mml:math></disp-formula>wherein, <italic>J</italic><sub>0</sub> and <italic>D</italic><sub>0</sub> are respectively the moments of inertia and damping coefficient of the VSG during stable operation; <italic>a</italic> and <italic>b</italic> are the adjustment coefficients of the moment of inertia and damping coefficient, respectively; <italic>K</italic><sub>J</sub> and <italic>K</italic><sub>D</sub> are change thresholds. The collaborative adaptive control flow chart is shown in <xref ref-type="fig" rid="fig-11">Fig. 11</xref>, where <italic>K</italic><sub>J</sub> &#x003D; 2.5, <italic>K</italic><sub>D</sub> &#x003D; 0.1, <italic>a</italic> &#x003D; 0.05, and <italic>b</italic> &#x003D; 5.</p>
<fig id="fig-11">
<label>Figure 11</label>
<caption>
<title>Flow chart of adaptive control for J and D</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-11.tif"/>
</fig>
<p><xref ref-type="fig" rid="fig-10">Fig. 10</xref> shows the output comparison of the VSG under different control strategies. The initial SOC of VSG under each control strategy is 50%. The simulation duration is 5 s, the load active power is 5 kW, and the reactive power is 1 kvar. The load is increased by 5 kW in 2 s, and the set value is restored in 4 s.</p>

<p>It can be seen from <xref ref-type="fig" rid="fig-10">Fig. 10a</xref> that VSG under parameter adaptive control consumes significantly less power. <xref ref-type="fig" rid="fig-10">Fig. 10b</xref> shows the response curves of the VSG output active power under different control strategies. When the load power suddenly increases by 5 kw at 2 s, the active power of C-VSG, S-VSG and I-VSG drops to 100,830 kW, 100,826 W and 101,077 W, respectively, and recovers to the rated output in 2.41, 2.43 and 2.39 s, respectively.</p>

<p><xref ref-type="fig" rid="fig-10">Fig. 10c</xref> shows that the adaptive change has a significant impact on the frequency change rate. When the frequency is subject to the initial load disturbance, the buffer effect of C-VSG control on the disturbance is general, and the maximum overshoot is 0.33 Hz; S-VSG control can better buffer frequency fluctuation, and the maximum overshoot is 0.25 Hz; The fuzzy adaptive VSG control has the best buffering effect on frequency fluctuation, and the maximum overshoot is 0.22 Hz.</p>

<p>When the load power suddenly increases, the frequency deviation of C-VSG exceeds 0.1 Hz, and the oscillation duration is about 0.6 s. As for S-VSG, its maximum frequency deviation is reduced to 0.093 Hz, and the adjustment time is reduced to 0.51 s. When I-VSG is used, the maximum frequency deviation is reduced to 0.086 Hz, and the adjustment time is reduced to 0.42 s. When the frequency deviation exceeds 0.1 Hz, the adaptive increase of D can significantly restrain the frequency deviation. Moreover, when the disturbance is eliminated, the fuzzy adaptive VSG control strategy can restore the frequency to stability as quickly as possible.</p>
<p>The range of RoCoF is also specified in <xref ref-type="fig" rid="fig-10">Fig.10d</xref> to ensure the system stability; the absolute maximum RoCoF with C-VSG method is 4.27 Hz/s. Using S-VSG and I-VSG control method, the maximum RoCoF becomes 3.50 and 4.00 Hz/s, respectively. Note that, using adaptive control, the maximum RoCoF has decreased by 6%&#x007E;18%, which improves the system frequency stability.</p>

<p><xref ref-type="fig" rid="fig-10">Figs. 10e</xref>, <xref ref-type="fig" rid="fig-10">10f</xref> show the changes of J and D under S-VSG control and I-VSG control. In <xref ref-type="fig" rid="fig-10">Fig. 10e</xref>, there are sharp spikes of the J waveform of S-VSG, which is caused by the uneven frequency change of the virtual rotor, sudden changes of d<italic>&#x03C9;</italic>/d<italic>t</italic> at the frequency inflection point.</p>

<p>Finally, this paper verifies the hierarchical power reduction control strategy when the SOC of the battery is in the buffer zone, as shown in <xref ref-type="fig" rid="fig-12">Fig. 12</xref>.</p>
<fig id="fig-12">
<label>Figure 12</label>
<caption>
<title>Comparisons for power and frequency with low SOC</title>
</caption>
<graphic mimetype="image" mime-subtype="tif" xlink:href="EE_27365-fig-12.tif"/>
</fig>
<p>The system is initially loaded with a 3 kW load. In the beginning, VSG operates with<italic>P</italic><sub>ref</sub> &#x003D; 16 kW and <italic>Q</italic><sub>ref</sub> &#x003D; 0 var. At 10 s, the system frequency suddenly increases by 0.2 Hz, and the frequency is recovered at 13 s.</p>
<p>To obtain a more obvious battery SOC curve, small capacity batteries are selected as the power supply in this test, and the graded power reduction operation strategy is started when the SOC reaches 7% and 6%, respectively.</p>
<p>By observing <xref ref-type="fig" rid="fig-12">Fig. 12a</xref>, it can be found that the SOC of C-VSG control reaches the minimum value at 27.37 s, then the battery system is in a limited discharge state and stops discharging and actively disconnects from the converter, and from <xref ref-type="fig" rid="fig-12">Fig. 12c</xref> we can see its output power drops sharply to a negative value, which means at that time, the energy storage converter did not supply power to the power grid, but consumed power from the power grid. The connection between the converter and the power grid should be immediately disconnected to reduce power loss. In practice, when the converter is going to stop power transmission, it should first disconnect the connection with the power grid side, then disconnect the DC input switch, and finally, the energy storage device can disconnect the interaction with the outside system.</p>

<p>More so, as shown in <xref ref-type="fig" rid="fig-12">Fig. 12b</xref>, from 27.37 s C-VSG&#x2019;s frequency management loses energy support, and the grid frequency is no longer maintained near the rated value.</p>

<p>Also, it can be seen from <xref ref-type="fig" rid="fig-12">Fig. 12c</xref> that the I-VSG control starts active power reduction operation at 5.326 s and reduces its output power again at 23.075 s. At 49.594 s, the SOC reaches the lowest value that can VSG be operated, and the battery is disconnected from the converter. Due to different control strategies, I-VSG markedly prolongs the discharge working time. Thus, the power shortage in unit time is reduced, and the impact of VSG shutdown on the power grid is reduced.</p>

<p>Due to the use of overdamping control, the system output power under the I-VSG control strategy does not overshoot, and its power change rate is smaller, as shown in <xref ref-type="fig" rid="fig-12">Fig. 12d</xref>. So the requirement for battery energy storage margin is smaller. The improved VSG control strategy effectively smoothes the short-term power fluctuation caused by the disconnection of the ionization network of the battery and retains sufficient frequency modulation capability during a sudden change in system frequency.</p>

<p>In the actual design of an energy storage power station, the capacity of the energy storage battery pack is far greater than 1Ah, and normally the installed capacity varies from 1 to 100 MW. Therefore, it is not only necessary to reasonably design the value of SOC<sub>b1</sub> and SOC<sub>b2</sub>, so that the operation time of each level of power reduction is about 30 s, but also to add multiple buffer zones when necessary, so as to limit the change rate of the active power of the energy storage plant to a reasonable range.</p>
</sec>
<sec id="s6">
<label>6</label>
<title>Conclusion</title>
<p>This article is based on fuzzy control theory, energy storage converter model and fuzzy adaptive control strategy of VSG, this paper studies the rationality of VSG control parameter selection when the energy storage unit is in a safe working state, and analyzes the relationship between two parameters, inertia and damping, and dynamic performance of the energy storage converter, so as to realize the adaptive VSG control strategy of the energy storage converter. The following conclusions are obtained through simulation analysis of a grid-connected model of an energy storage power station:
<list list-type="order">
<list-item>
<p>Considering the physical constraints of energy storage, the output power and its power change rate of VSG are always within the allowable range of battery performance. This control strategy sacrifices the dynamic regulation performance of the energy storage converter, thus ensuring the overall dynamic response capability of the whole energy storage system. Therefore, the energy storage battery will not have over-discharge and other adverse phenomena, which protects the safety and service life of the battery. In the actual construction of energy storage power stations, batteries with large capacity, high discharge rate and large threshold of power change rate should be selected as the power supply to cut down the negative impact on the output performance of the energy storage system.</p></list-item>
<list-item>
<p>When the load changes suddenly, fuzzy adaptive VSG control can better improve response characteristics of power grid frequency and output power, ultimately enhancing the robustness of energy storage systems.</p></list-item>
<list-item>
<p>The power grading and over damping operation described in this paper can also be applied to the shutdown preparation of the energy storage converter, and future research work will focus on the adaptive P-&#x03C9; Control coefficient and Q-E control coefficient based on SOC.</p></list-item>
</list></p>
</sec>
</body>
<back>
<sec><title>Funding Statement</title>
<p>This project is supported by the <funding-source>Science and Technology Project of State Grid Corporation of China</funding-source> (<award-id>W22KJ2722005</award-id>); <funding-source>Tianyou Innovation Team of Lanzhou Jiaotong University</funding-source> (<award-id>TY202009</award-id>).</p>
</sec>
<sec sec-type="COI-statement"><title>Conflicts of Interest</title>
<p>The authors declare that they have no conflicts of interest to report regarding the present study.</p>
</sec>
<ref-list content-type="authoryear">
<title>References</title>
<ref id="ref-1"><label>1.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Liu</surname>, <given-names>Z. C.</given-names></string-name>, <string-name><surname>Peng</surname>, <given-names>D. G.</given-names></string-name>, <string-name><surname>Zhao</surname>, <given-names>H. R.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>D. H.</given-names></string-name>, <string-name><surname>Liu</surname>, <given-names>Y. C.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Development prospects of energy storage participating in power system auxiliary services under dual carbon goals</article-title>. <source>Energy Storage Science and Technology</source><italic>,</italic> <fpage>1</fpage>&#x2013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.19799/j.cnki.2095-42392021.0431</pub-id></mixed-citation></ref>
<ref id="ref-2"><label>2.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Gao</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>H. Y.</given-names></string-name>, <string-name><surname>Wang</surname>, <given-names>Z. Q.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Planning for site selection and capacity determination of distributed energy storage in regional power grid</article-title>. <source>High Voltage Apparatus</source><italic>,</italic> <volume>56</volume><italic>(</italic><issue>8</issue><italic>),</italic> <fpage>52</fpage>&#x2013;<lpage>58</lpage>.</mixed-citation></ref>
<ref id="ref-3"><label>3.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Xi</surname>, <given-names>X. X.</given-names></string-name>, <string-name><surname>Xiong</surname>, <given-names>M. P.</given-names></string-name>, <string-name><surname>Yuan</surname>, <given-names>J. H.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Economy analysis of energy storage system in wind farm generation side</article-title>. <source>Smart Power</source><italic>,</italic> <volume>48</volume><italic>(</italic><issue>11</issue><italic>),</italic> <fpage>16</fpage>&#x2013;<lpage>21&#x0002B;47</lpage>.</mixed-citation></ref>
<ref id="ref-4"><label>4.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>He</surname>, <given-names>K. X.</given-names></string-name>, <string-name><surname>Ma</surname>, <given-names>S. L.</given-names></string-name>, <string-name><surname>Ma</surname>, <given-names>Z.</given-names></string-name>, <string-name><surname>Aoyu</surname>, <given-names>X. U. E.</given-names></string-name></person-group> (<year>2021</year>). <article-title>Energy storage technology development trend and policy environment analysis</article-title>. <source>Distributed Energy</source><italic>,</italic> <volume>6</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>45</fpage>&#x2013;<lpage>52</lpage>.</mixed-citation></ref>
<ref id="ref-5"><label>5.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Larcher</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Tarascon</surname>, <given-names>J. M.</given-names></string-name></person-group> (<year>2015</year>). <article-title>Towards greener and more sustainable batteries for electrical energy storage</article-title>. <source>Nature Chemistry</source><italic>,</italic> <volume>7</volume><italic>(</italic><issue>1</issue><italic>),</italic> <fpage>19</fpage>&#x2013;<lpage>29</lpage>. <pub-id pub-id-type="doi">10.1038/nchem.2085</pub-id>; <pub-id pub-id-type="pmid">25515886</pub-id></mixed-citation></ref>
<ref id="ref-6"><label>6.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Tang</surname>, <given-names>B.</given-names></string-name>, <string-name><surname>Shan</surname>, <given-names>L.</given-names></string-name>, <string-name><surname>Liang</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>J.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Issues and opportunities facing aqueous zinc-ion batteries</article-title>. <source>Energy &#x00026; Environmental Science</source><italic>,</italic> <volume>12</volume><italic>(</italic><issue>11</issue><italic>),</italic> <fpage>3288</fpage>&#x2013;<lpage>3304</lpage>. <pub-id pub-id-type="doi">10.1039/C9EE02526J</pub-id></mixed-citation></ref>
<ref id="ref-7"><label>7.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Blanc</surname>, <given-names>L. E.</given-names></string-name>, <string-name><surname>Kundu</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Nazar</surname>, <given-names>L. F.</given-names></string-name></person-group> (<year>2020</year>). <article-title>Scientific challenges for the implementation of Zn-ion batteries</article-title>. <source>Joule</source><italic>,</italic> <volume>4</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>771</fpage>&#x2013;<lpage>799</lpage>. <pub-id pub-id-type="doi">10.1016/j.joule.2020.03.002</pub-id></mixed-citation></ref>
<ref id="ref-8"><label>8.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname>, <given-names>J. L.</given-names></string-name>, <string-name><surname>Yuan</surname>, <given-names>X. D.</given-names></string-name>, <string-name><surname>Yu</surname>, <given-names>Z. G.</given-names></string-name>, <string-name><surname>Ge</surname>, <given-names>L.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Comments on power quality enhancement research for power grid by energy storage system</article-title>. <source>Automation of Electric Power Systems</source><italic>,</italic> <volume>43</volume><italic>(</italic><issue>8</issue><italic>),</italic> <fpage>15</fpage>&#x2013;<lpage>24</lpage>.</mixed-citation></ref>
<ref id="ref-9"><label>9.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Luo</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>J. L.</given-names></string-name></person-group> (<year>2012</year>). <article-title>Research progress of energy storage technology in power system</article-title>. <source>Power System and Clean Energy</source><italic>,</italic> <volume>28</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>71</fpage>&#x2013;<lpage>79</lpage>.</mixed-citation></ref>
<ref id="ref-10"><label>10.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhong</surname>, <given-names>Q. C.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Power-electronics-enabled autonomous power systems: Architecture and technical routes</article-title>. <source>IEEE Transactions on Industrial Electronics</source><italic>,</italic> <volume>64</volume><italic>(</italic><issue>7</issue><italic>),</italic> <fpage>5907</fpage>&#x2013;<lpage>5918</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2017.2677339</pub-id></mixed-citation></ref>
<ref id="ref-11"><label>11.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Delille</surname>, <given-names>G.</given-names></string-name>, <string-name><surname>Francois</surname>, <given-names>B.</given-names></string-name>, <string-name><surname>Malarange</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2012</year>). <article-title>Dynamic frequency control support by energy storage to reduce the impact of wind and solar generation on isolated power system&#x2019;s inertia</article-title>. <source>IEEE Transactions on Sustainable Energy</source><italic>,</italic> <volume>3</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>931</fpage>&#x2013;<lpage>939</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2012.2205025</pub-id></mixed-citation></ref>
<ref id="ref-12"><label>12.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhong</surname>, <given-names>Q. C.</given-names></string-name>, <string-name><surname>Weiss</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2011</year>). <article-title>Synchronverters: Inverters that mimic synchronous generators</article-title>. <source>IEEE Transactions on Industrial Electronics</source><italic>,</italic> <volume>58</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>1259</fpage>&#x2013;<lpage>1267</lpage>. <pub-id pub-id-type="doi">10.1109/TIE.2010.2048839</pub-id></mixed-citation></ref>
<ref id="ref-13"><label>13.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Zhong</surname>, <given-names>Q. C.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Virtual synchronous machines and autonomous power systems</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>37</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>336</fpage>&#x2013;<lpage>349</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.162325</pub-id></mixed-citation></ref>
<ref id="ref-14"><label>14.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Torres</surname>, <given-names>M. A. L.</given-names></string-name>, <string-name><surname>Lopes</surname>, <given-names>L. A. C.</given-names></string-name>, <string-name><surname>Moran</surname>, <given-names>L. A. T.</given-names></string-name>, <string-name><surname>Espinoza</surname>, <given-names>J. R. C.</given-names></string-name></person-group> (<year>2014</year>). <article-title>Self-tuning virtual synchronous machine: A control strategy for energy storage systems to support dynamic frequency control</article-title>. <source>IEEE Transactions on Energy Conversion</source><italic>,</italic> <volume>29</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>833</fpage>&#x2013;<lpage>840</lpage>. <pub-id pub-id-type="doi">10.1109/TEC.2014.2362577</pub-id></mixed-citation></ref>
<ref id="ref-15"><label>15.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Markovic</surname>, <given-names>U.</given-names></string-name>, <string-name><surname>Chu</surname>, <given-names>Z. D.</given-names></string-name>, <string-name><surname>Aristidou</surname>, <given-names>P.</given-names></string-name>, <string-name><surname>Hug</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2019</year>). <article-title>LQR-based adaptivevirtual synchronous machine for power systems with high inverter penetration</article-title>. <source> IEEE Transactions on Sustainable Energy</source><italic>,</italic> <volume>10</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>1501</fpage>&#x2013;<lpage>1512</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2018.2887147</pub-id></mixed-citation></ref>
<ref id="ref-16"><label>16.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Shi</surname>, <given-names>R. L.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Hu</surname>, <given-names>C.</given-names></string-name>, <string-name><surname>Xu</surname>, <given-names>H. Z.</given-names></string-name>, <string-name><surname>Gu</surname>, <given-names>J.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2018</year>). <article-title>Self-tuning virtual synchronous generator control for improving frequency stability in autonomous photovoltaic-diesel microgrids</article-title>. <source>Journal of Modern Power Systems and Clean Energy</source><italic>,</italic> <volume>6</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>482</fpage>&#x2013;<lpage>494</lpage>. <pub-id pub-id-type="doi">10.1007/s40565-017-0347-3</pub-id></mixed-citation></ref>
<ref id="ref-17"><label>17.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Li</surname>, <given-names>J.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>S.</given-names></string-name>, <string-name><surname>Xu</surname>, <given-names>S.</given-names></string-name></person-group> (<year>2008</year>). <article-title>Fuzzy control system design via fuzzy lyapunov functions</article-title>. <source>IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics)</source><italic>,</italic> <volume>38</volume><italic>(</italic><issue>6</issue><italic>),</italic> <fpage>1657</fpage>&#x2013;<lpage>1661</lpage>. <pub-id pub-id-type="doi">10.1109/TSMCB.2008.928224</pub-id>; <pub-id pub-id-type="pmid">19022736</pub-id></mixed-citation></ref>
<ref id="ref-18"><label>18.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Feng</surname>, <given-names>G.</given-names></string-name></person-group> (<year>2006</year>). <article-title>A survey on analysis and design of model-based fuzzy control systems</article-title>. <source>IEEE Transactions on Fuzzy Systems</source><italic>,</italic> <volume>14</volume><italic>(</italic><issue>5</issue><italic>),</italic> <fpage>676</fpage>&#x2013;<lpage>697</lpage>. <pub-id pub-id-type="doi">10.1109/TFUZZ.2006.883415</pub-id></mixed-citation></ref>
<ref id="ref-19"><label>19.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Andalib-Bin-Karim</surname>, <given-names>C.</given-names></string-name>, <string-name><surname>Liang</surname>, <given-names>X.</given-names></string-name>, <string-name><surname>Zhang</surname>, <given-names>H.</given-names></string-name></person-group> (<year>2018</year>). <article-title>Fuzzy-secondary-controller-based virtual synchronous generator control scheme for interfacing inverters of renewable distributed generation in microgrids</article-title>. <source> IEEE Transactions on Industry Applications</source><italic>,</italic> <volume>54</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>1047</fpage>&#x2013;<lpage>1061</lpage>. <pub-id pub-id-type="doi">10.1109/TIA.2017.2773432</pub-id></mixed-citation></ref>
<ref id="ref-20"><label>20.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Pradhan</surname>, <given-names>C. C.</given-names></string-name>, <string-name><surname>Bhende</surname>, <given-names>N.</given-names></string-name>, <string-name><surname>Samanta</surname>, <given-names>A. K.</given-names></string-name></person-group> (<year>2018</year>). <article-title>Adaptive virtual inertia-based frequency regulation in wind power systems</article-title>. <source>Renewable Energy</source><italic>,</italic> <volume>115</volume><italic>(</italic><issue>5</issue><italic>),</italic> <fpage>558</fpage>&#x2013;<lpage>574</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2017.08.065</pub-id></mixed-citation></ref>
<ref id="ref-21"><label>21.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Karimi</surname>, <given-names>A.</given-names></string-name>, <string-name><surname>Khayat</surname>, <given-names>Y.</given-names></string-name>, <string-name><surname>Naderi</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Dragi&#x010D;evi&#x0107;</surname>, <given-names>T.</given-names></string-name>, <string-name><surname>Mirzaei</surname>, <given-names>R.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2020</year>). <article-title>Inertia response improvement in AC microgrids: A fuzzy-based virtual synchronous generator control</article-title>. <source>IEEE Transactions on Power Electronics</source><italic>,</italic> <volume>35</volume><italic>(</italic><issue>4</issue><italic>),</italic> <fpage>4321</fpage>&#x2013;<lpage>4331</lpage>. <pub-id pub-id-type="doi">10.1109/TPEL.2019.2937397</pub-id></mixed-citation></ref>
<ref id="ref-22"><label>22.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Kerdphol</surname>, <given-names>T.</given-names></string-name>, <string-name><surname>Watanabe</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Hongesombut</surname>, <given-names>K.</given-names></string-name>, <string-name><surname>Mitani</surname>, <given-names>Y.</given-names></string-name></person-group> (<year>2019</year>). <article-title>Self-adaptive virtual inertia control-based fuzzy logic to improve frequency stability of microgrid with high renewable penetration</article-title>. <source>IEEE Access</source><italic>,</italic> <volume>7</volume><italic>,</italic> <fpage>76071</fpage>&#x2013;<lpage>76083</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2019.2920886</pub-id></mixed-citation></ref>
<ref id="ref-23"><label>23.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Ebrahimi</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Khajehoddin</surname>, <given-names>S. A.</given-names></string-name>, <string-name><surname>Karimi-Ghartemani</surname>, <given-names>M.</given-names></string-name></person-group> (<year>2019</year>). <article-title>An improved damping method for virtual synchronous machines</article-title>. <source>IEEE Transactions on Sustainable Energy</source><italic>,</italic> <volume>10</volume><italic>(</italic><issue>3</issue><italic>),</italic> <fpage>1491</fpage>&#x2013;<lpage>1500</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2019.2902033</pub-id></mixed-citation></ref>
<ref id="ref-24"><label>24.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Hu</surname>, <given-names>A. P.</given-names></string-name>, <string-name><surname>Yang</surname>, <given-names>B.</given-names></string-name>, <string-name><surname>Pan</surname>, <given-names>P. P.</given-names></string-name>, <string-name><surname>Li</surname>, <given-names>G. J.</given-names></string-name>, <string-name><surname>Tao</surname>, <given-names>Y. B.</given-names></string-name> <etal>et al.</etal></person-group> (<year>2018</year>). <article-title>Study on inertial characteristics of energy storage system with power electronic interface</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>38</volume><italic>(</italic><issue>17</issue><italic>),</italic> <fpage>4999</fpage>&#x2013;<lpage>5008&#x002B;5297</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.172143</pub-id></mixed-citation></ref>
<ref id="ref-25"><label>25.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Yuan</surname>, <given-names>C.</given-names></string-name>, <string-name><surname>Liu</surname>, <given-names>C.</given-names></string-name>, <string-name><surname>Zhao</surname>, <given-names>T. Y.</given-names></string-name>, <string-name><surname>Chen</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Xiao</surname>, <given-names>X. N.</given-names></string-name></person-group> (<year>2017</year>). <article-title>Research on operating boundary of virtual synchronous machine based on physical constraint of energy storage system</article-title>. <source>Proceedings of the CSEE</source><italic>,</italic> <volume>37</volume><italic>(</italic><issue>2</issue><italic>),</italic> <fpage>506</fpage>&#x2013;<lpage>516</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.161121</pub-id></mixed-citation></ref>
<ref id="ref-26"><label>26.</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name><surname>Chen</surname>, <given-names>M.</given-names></string-name>, <string-name><surname>Zhou</surname>, <given-names>D.</given-names></string-name>, <string-name><surname>Blaabjerg</surname>, <given-names>F.</given-names></string-name></person-group> (<year>2021</year>). <article-title>A decentralized adaptive soc balancing strategy in VSG-based islanded power system</article-title>. <source>IECON 2021&#x2013;47th Annual Conference of the IEEE Industrial Electronics Society</source>, vol. <comment>2021</comment>, pp. <fpage>1</fpage>&#x2013;<lpage>6</lpage>. <conf-loc>Toronto, ON, Canada</conf-loc>. <pub-id pub-id-type="doi">10.1109/IECON48115.2021.9589320</pub-id></mixed-citation></ref>
</ref-list>
</back>
</article>