<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Energy Res.</journal-id>
<journal-title>Frontiers in Energy Research</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Energy Res.</abbrev-journal-title>
<issn pub-type="epub">2296-598X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1352385</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2024.1352385</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Energy Research</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Virtual inertia control to active support of the variable-speed wind turbine in variable frequency limit time</article-title>
<alt-title alt-title-type="left-running-head">Tingyun et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2024.1352385">10.3389/fenrg.2024.1352385</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Tingyun</surname>
<given-names>Gu</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2598149/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Houyi</surname>
<given-names>Zhang</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bowen</surname>
<given-names>Li</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Junyi</surname>
<given-names>Mao</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Qiang</surname>
<given-names>Fan</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yutao</surname>
<given-names>Xu</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Qihui</surname>
<given-names>Feng</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
</contrib-group>
<aff>
<institution>Electric Power Research Institute of Guizhou Power Grid Co., Ltd.</institution>, <addr-line>Guiyang</addr-line>, <addr-line>Guizhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1822992/overview">Yaser Qudaih</ext-link>, Higher Colleges of Technology, United Arab Emirates</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2602504/overview">Mansoor Janjua</ext-link>, Higher Colleges of Technology, United Arab Emirates</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/881527/overview">Gaber Magdy</ext-link>, Aswan University, Egypt</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Gu Tingyun, <email>gutingyun@126.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>01</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1352385</elocation-id>
<history>
<date date-type="received">
<day>08</day>
<month>12</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>03</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Tingyun, Houyi, Bowen, Junyi, Qiang, Yutao and Qihui.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Tingyun, Houyi, Bowen, Junyi, Qiang, Yutao and Qihui</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Constraining the virtual inertia of variable-speed wind turbines within the variable frequency limit time will be the key to improving the wind power&#x2019;s reliability in actively supporting the system frequency. First, this paper analyzes the controllable characteristics and evaluation methods of the virtual inertia of the wind turbine. Second, it quantifies the variable frequency limit time of the grid with high-proportion wind power and considers the wind turbine speed and power margin to restrict the virtual inertia of the wind turbine during the period of frequency drop or increase to provide reliable inertia support to ensure the safety of the system frequency. On this basis, the inertia demand under the frequency safety of the system was analyzed, and the virtual inertia control strategy of the variable-speed wind turbine in the variable frequency limit time was proposed using the speed tracking performance of the wind turbine. Finally, the grid with a high-proportion wind power simulation system is built to verify that the variable-speed wind turbine can reliably satisfy the inertia demand of system frequency modulation under the proposed control strategy and significantly improve its active support performance to the power grid.</p>
</abstract>
<kwd-group>
<kwd>wind power generation</kwd>
<kwd>virtual inertia</kwd>
<kwd>frequency support</kwd>
<kwd>rotor speed tracking</kwd>
<kwd>power grid</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Smart Grids</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Improving the grid-connected support capability of wind turbines has become a key technology for enabling high-proportion wind power consumption and achieving the &#x201c;30&#xb7;60 carbon target.&#x201d; Under maximum power tracking control, large-scale random and intermittent wind power is injected into the power grid, and the synchronous generator is unable to ensure the safe operation of the system due to the reduction of the installed proportion (<xref ref-type="bibr" rid="B3">Fern&#xe1;ndez-Guillam&#xf3;n et al., 2019</xref>; <xref ref-type="bibr" rid="B11">Mensou et al., 2020</xref>; <xref ref-type="bibr" rid="B32">Zhu et al., 2021</xref>). At present, the frequency active support control of variable-speed wind turbines has attracted much attention.</p>
<p>Variable-speed wind turbines can actively participate in the frequency adjustment of the power grid by adding virtual inertia and primary frequency controllers based on power tracking control (<xref ref-type="bibr" rid="B26">Yuan et al., 2017</xref>; <xref ref-type="bibr" rid="B13">Ochoa and Martinez, 2018</xref>; <xref ref-type="bibr" rid="B27">Yuan et al., 2018</xref>). In the early stage of frequency change, if the wind turbine does not have an inertial response, it will lead to a significant increase in the frequency change rate, and a large frequency drop or increase will seriously threaten the safety of system operation (<xref ref-type="bibr" rid="B17">Sun et al., 2010</xref>; <xref ref-type="bibr" rid="B18">Tielens and Hertem, 2016</xref>; <xref ref-type="bibr" rid="B2">Das et al., 2020</xref>). In order to solve the problem of low inertia in wind turbines, <xref ref-type="bibr" rid="B4">Fu et al. (2017)</xref> and <xref ref-type="bibr" rid="B14">Qi et al. (2022)</xref> defined the virtual inertia of wind turbines. Different from the constant inherent inertia of the synchronous generator, the virtual moment of inertia through additional control not only depends on the inherent inertia, initial speed, and control parameters of the wind turbine but is also closely related to the change in the system frequency, resulting in its complex and changeable characteristics while also renders it more flexible and controllable. <xref ref-type="bibr" rid="B10">Mauricio et al. (2009)</xref>; <xref ref-type="bibr" rid="B5">Gautam et al. (2011)</xref>; and <xref ref-type="bibr" rid="B6">Li et al. (2016)</xref> proposed a virtual inertia control strategy for wind turbines and used the differential link df/dt to simulate the rotor motion characteristics of synchronous generators. On this basis, <xref ref-type="bibr" rid="B19">ULLAH et al. (2008)</xref>; <xref ref-type="bibr" rid="B9">Makrini et al. (2017)</xref> further added a speed recovery module to ensure the speed recovery of the wind turbine after the inertial response, avoiding the secondary frequency drop caused by excessive inertial support during system recovery. Effective inertia support can avoid large and rapid changes in frequency, and allowing for sufficient adjustment time for subsequent primary frequency modulation. In order to reduce the static deviation of the system frequency after primary frequency modulation, the same active frequency droop control as the synchronous generator was introduced into the active power control of the wind turbine. <xref ref-type="bibr" rid="B15">Ramtharan et al. (2007)</xref> and <xref ref-type="bibr" rid="B28">Zertek et al. (2012)</xref> proposed a deloading control method for wind turbines and simulated the primary frequency modulation characteristics of synchronous generators through variable pitch control. <xref ref-type="bibr" rid="B1">Almeida and Lopes (2007)</xref> and <xref ref-type="bibr" rid="B29">Zhang et al. (2012)</xref> used the overspeed control of the wind turbine to complete the deloading operation of the wind turbine and cooperated with the virtual inertia to improve the frequency modulation performance of the wind turbine. In theory, the power electronic converter gives the wind turbine flexible power control performance. However, the characteristics of virtual frequency modulation are affected by many factors, and its engineering practicability still requires to be explored.</p>
<p>The performance evaluation of the virtual inertia and primary frequency modulation of the wind turbine is generally difficult to solve (<xref ref-type="bibr" rid="B23">Wu and Infield, 2014</xref>; <xref ref-type="bibr" rid="B31">Zhu et al., 2018</xref>). Among them, the additional primary frequency control is relatively simple. When the load-shedding reserve capacity is sufficient, the wind turbine can have the same power distribution capability as the synchronous generator using the active power frequency droop characteristics (<xref ref-type="bibr" rid="B20">Vidyanandan and Senroy, 2013</xref>; <xref ref-type="bibr" rid="B24">Xiao et al., 2017</xref>). However, it is difficult to evaluate the virtual inertia of the wind turbine, which requires discussion at present. Because the initial wind speed determines the kinetic energy reserve of the wind turbine and the wind turbine speed will be affected by power tracking control, additional inertia control, and system frequency characteristics at the same time, the virtual inertia is always in dynamic change, which is not only complicated in quantification but also has operating risk. <xref ref-type="bibr" rid="B8">Ma et al. (2017)</xref>; <xref ref-type="bibr" rid="B30">Zhu et al. (2017)</xref>; and <xref ref-type="bibr" rid="B25">Ying et al. (2018)</xref> analyzed that although the introduction of virtual inertia can reduce the frequency change rate, it has the risk of aggravating the system power oscillation. <xref ref-type="bibr" rid="B16">Shi et al. (2011)</xref> and <xref ref-type="bibr" rid="B7">Lucas et al. (2019)</xref> used the small disturbance analysis method to prove that the virtual inertia of the wind turbine would lead to the right shift of the characteristic root and weaken the stability of the system. Reasonable quantification of virtual inertia can not only provide a basis for accurately predicting the effect of frequency active support but also be the key to evaluating the operation risk of additional control. <xref ref-type="bibr" rid="B22">Wen and Lin (2021)</xref> and <xref ref-type="bibr" rid="B21">Wang et al. (2022)</xref> proposed quantitative requirements for virtual inertia to ensure the frequency security of the grid through constraints such as frequency change rate and minimum value. However, the expected virtual inertia support is not only a prerequisite but also an urgent problem to be solved. Although the source of virtual inertia has been extended to various renewed energy sources and loads, the nonlinear time-varying characteristics make it difficult to quantify the inertia, and a reasonable time scale is urgently needed to constrain the multi-source inertia to complete the evaluation (<xref ref-type="bibr" rid="B22">Wen and Lin, 2021</xref>). In summary, the wind power frequency active support technology is not mature at present, and the promotion and application still lack the quantification and control methods of virtual inertia.</p>
<p>In order to improve the reliability of wind turbine frequency active support control according to the system frequency safety requirements, this paper proposes a quantitative method and control strategy of virtual inertia using the inertia support time and wind turbine operation characteristics constraints. First, based on the definition of virtual inertia of the wind turbine, the frequency characteristics of the wind power grid-connected system and the problems of inertia evaluation are analyzed. Second, the frequency conversion extreme time is calculated, and the wind turbine speed and power margin are considered to constrain and quantify the wind turbine virtual inertia. Third, according to the virtual inertia demand and quantitative index, the virtual inertia control strategy of the variable-speed wind turbine in the extreme time of frequency conversion is proposed by speed tracking control. Finally, a simulation system with high-proportion wind power is built to verify the improvement effect of the proposed evaluation method and control strategy on the reliability of wind turbine frequency active support.</p>
<p>The remainder of this paper is organized as follows: <xref ref-type="sec" rid="s2">Section 2</xref> presents the virtual inertia evaluation method for wind turbines. <xref ref-type="sec" rid="s3">Section 3</xref> presents the calculation method of frequency conversion extreme time for the system with high-proportion wind power. <xref ref-type="sec" rid="s4">Section 4</xref> presents the virtual inertia constraint and evaluation of wind turbines. <xref ref-type="sec" rid="s5">Section 5</xref> presents the virtual inertia demand and control of wind turbines. <xref ref-type="sec" rid="s6">Section 6</xref> presents the numerical results and demonstrates the effectiveness of the proposed method. <xref ref-type="sec" rid="s7">Section 7</xref> concludes the paper and presents future work.</p>
</sec>
<sec id="s2">
<title>2 Virtual inertia evaluation of the wind turbine</title>
<sec id="s2-1">
<title>2.1 Wind turbine frequency active support characteristics</title>
<p>After the system with a high proportion of wind power encounters an active power disturbance, the system frequency support should be completed by the synchronous generator set and wind turbine. In order to improve the utilization rate of wind energy, wind farms usually do not use the deloading operation mode to reserve for primary frequency regulation. In this mode, the virtual inertial response will be the frequency active support function that wind power urgently needs. In the wind power grid-connected system, the synchronous generator and the wind turbine should jointly complete the inertial support, and the primary frequency modulation function is borne by the synchronous generator and the energy storage.</p>
<p>Taking a short-term frequency drop as an example, the dynamic response of wind turbines and synchronous generators participating in frequency adjustment is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. In this figure, <italic>t</italic>
<sub>0</sub> is the frequency drop time, <italic>t</italic>
<sub>b</sub> is the primary frequency modulation action time, <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is the time when the frequency drops to the lowest value, and <italic>t</italic>
<sub>p</sub> is the end time of the primary frequency modulation.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Power support of generators and dynamic response of system frequency.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g001.tif"/>
</fig>
<p>At time <italic>t</italic>
<sub>0</sub>, the load suddenly increases &#x394;<italic>P</italic>
<sub>d</sub>, and the system frequency drops sharply. The synchronous generator set and wind turbine need to respond quickly to the frequency change, and the system power demand is compensated by the fast active power support, that is, &#x394;<italic>P</italic>
<sub>e</sub> &#x3d; &#x394;<italic>P</italic>
<sub>d</sub>.</p>
<p>At the <italic>t</italic>
<sub>0</sub>&#x223c;<italic>t</italic>
<sub>
<italic>f</italic>m</sub> stage, after the wind turbine starts the inertia support control, it will share the unbalanced power borne by the synchronous generator, thereby slowing down the system frequency drop speed. The common virtual inertia control of wind turbines adopts a differential link, and the power response is defined as Eq. <xref ref-type="disp-formula" rid="e1">(1)</xref>:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where &#x394;<italic>P</italic>
<sub>w</sub> is the wind turbine inertial support power and <italic>K</italic>
<sub>I</sub> is the differential control coefficient.</p>
<p>As shown in <xref ref-type="fig" rid="F1">Figure 1</xref>, during a frequency drop (d<italic>f</italic>/d<italic>t</italic> &#x3c;0), the power response of the virtual inertia controller is &#x394;<italic>P</italic>
<sub>w</sub> &#x3e;0. At this stage, the synchronous generator and the wind turbine jointly maintain the inertia support power and satisfy the load demand &#x394;<italic>P</italic>
<sub>d</sub>. The wind turbine reduces the kinetic energy demand of the synchronous generator by releasing the kinetic energy.</p>
<p>At the <italic>t</italic>
<sub>
<italic>f</italic>m</sub> &#x223c; <italic>t</italic>
<sub>p</sub> stage, during frequency recovery (d<italic>f</italic>/d<italic>t</italic> &#x3e;0), the power response of the virtual inertia controller is &#x394;<italic>P</italic>
<sub>w</sub> &#x3c;0. At this stage, the wind turbine absorbs power from the system under virtual inertia control, and the wind turbine changes from releasing kinetic energy to absorbing kinetic energy, which increases the primary frequency regulation burden of the synchronous generator and slows down the frequency recovery speed.</p>
<p>Considering that the virtual inertia has the problem of grid-connected security, at present, GB/T19963.1&#x2013;2021 &#x201c;wind farm access power system technical regulations&#x201d; have stipulated that the wind turbine is required to have the inertia support function after being connected to the system, but the additional controller must satisfy the following conditions: &#x394;<italic>f</italic>&#xd7;d<italic>f</italic>/d<italic>t</italic> &#x3e;0. The wind turbine should set the time <italic>t</italic>
<sub>
<italic>f</italic>m</sub> as the end time of the virtual inertia control to satisfy the grid-connected requirements.</p>
<p>It is worth noting that it is necessary to detect the frequency signal at the grid-connected point of the wind turbine when determining the end <italic>t</italic>
<sub>
<italic>f</italic>m</sub> of the inertial support. The differential signal has high-frequency noise in the field test. Combined with the actual test situation, it is difficult to determine whether &#x394;<italic>f</italic> &#xd7; d<italic>f</italic>/d<italic>t</italic> is greater than zero. In addition, at the end of the moment of inertia support, the wind turbine has released kinetic energy, resulting in a decrease in speed. After restoring to the maximum power tracking control, the output power of the wind turbine will still be lower than the initial level before the frequency drop. This is because the wind turbine speed has dropped. Therefore, even if the inertial control can be terminated in time, the wind turbine will still absorb the power, accelerate the rotor, store kinetic energy, and complete the maximum power tracking target; however, this will not be conducive to frequency recovery during primary frequency modulation.</p>
<p>It can be seen from <xref ref-type="fig" rid="F1">Figure 1</xref> that during the active frequency support period, the virtual inertia control of the wind turbine cannot continuously provide effective active power support for the system. Combined with the inertia start-up conditions specified in the grid-connected standard, this paper defines the effective time of the inertia response as the frequency conversion extreme time <italic>t</italic>
<sub>
<italic>f</italic>m</sub>; that is, after the system is disturbed, the frequency drops or increases from the initial value to the maximum frequency deviation.</p>
<p>Before adopting virtual inertia control, if the wind turbine can predict the extreme time of frequency conversion, it not only ensures the safe removal of the additional control and effectively avoids the misoperation caused by the large monitoring error of the frequency differential signal but also helps solve the problem of quantitative evaluation of virtual inertia, and then it provides a more reliable calculation method for the inertia demand under the system frequency safety.</p>
</sec>
<sec id="s2-2">
<title>2.2 Virtual inertia evaluation based on frequency conversion extremum time</title>
<p>The inertia of the wind power grid-connected system is composed of the inherent inertia of the synchronous generator and the virtual inertia of the wind turbine. The inertia time constant of the synchronous generator depends on the mechanical characteristics of the rotor, which can be expressed as Eq. <xref ref-type="disp-formula" rid="e2">(2)</xref> (<xref ref-type="bibr" rid="B32">Zhu et al., 2021</xref>)<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mtext>kg</mml:mtext>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>J</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <italic>J</italic>
<sub>g</sub> is the mechanical moment of inertia of the synchronous generator; <italic>&#x3c9;</italic>
<sub>e</sub> is the synchronous speed; <italic>p</italic>
<sub>g</sub> is the number of rotor poles; <italic>E</italic>
<sub>kg</sub> is the rotor kinetic energy stored at the rated speed of the synchronous generator; and <italic>S</italic>
<sub>g</sub> is the rated capacity of the synchronous generator.</p>
<p>Referring to the definition of the inertia time constant of a synchronous generator, the virtual inertia time constant <italic>H</italic>
<sub>vir</sub> of a variable-speed wind turbine can be expressed as (<xref ref-type="bibr" rid="B14">Qi et al., 2022</xref>)<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>J</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>where <italic>H</italic>
<sub>w</sub> is the inherent inertia time constant of the wind turbine; &#x394;<italic>&#x3c9;</italic>
<sub>r</sub> is the change in the wind turbine speed; <italic>&#x3c9;</italic>
<sub>r0</sub> is the initial speed of the wind turbine before the inertial response; &#x394;<italic>&#x3c9;</italic>
<sub>e</sub> is the synchronous speed variation; <italic>p</italic>
<sub>w</sub> is the number of wind turbine poles; <italic>S</italic>
<sub>w</sub> is the rated capacity of the wind turbine; and <italic>J</italic>
<sub>vir</sub> is the virtual moment of inertia of the wind turbine, which can be expressed as<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>J</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>J</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where <italic>J</italic>
<sub>w</sub> is the inherent moment of inertia of the wind turbine.</p>
<p>According to the definition of the virtual inertia of the wind turbine in <xref ref-type="disp-formula" rid="e3">(3)</xref> and <xref ref-type="disp-formula" rid="e4">(4)</xref>, it can be seen that if &#x394;<italic>&#x3c9;</italic>
<sub>r</sub>/&#x394;<italic>&#x3c9;</italic>
<sub>e</sub> &#x3e;&#x3e;1, the virtual inertia will be much larger than the inherent inertia of the wind turbine and even has better inertial support ability than the synchronous generator set. However, in practical control, <italic>H</italic>
<sub>vir</sub> is determined by <italic>&#x3c9;</italic>
<sub>r0</sub>, <italic>H</italic>
<sub>w</sub>, and &#x394;<italic>&#x3c9;</italic>
<sub>r</sub>/&#x394;<italic>&#x3c9;</italic>
<sub>e</sub>, where <italic>&#x3c9;</italic>
<sub>r0</sub> depends on the real-time wind speed, which can be expressed as Eq. <xref ref-type="disp-formula" rid="e5">(5)</xref>
<disp-formula id="e5">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mi>R</mml:mi>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where &#x3bb;<sub>opt</sub> is the optimal tip speed ratio; <italic>v</italic> is the real-time wind speed; <italic>R</italic> is the radius of the wind wheel; <italic>&#x3c9;</italic>
<sub>max</sub> is the maximum speed allowed for the stable operation of the wind turbine; <italic>v</italic>
<sub>min</sub> is the cut-in wind speed; and <italic>v</italic>
<sub>max</sub> is the maximum wind speed before the pitch control starts.</p>
<p>In the control process, it is necessary to comprehensively consider the dynamic changes in the inherent inertia <italic>H</italic>
<sub>w</sub>, wind speed, additional power, and system frequency of the wind turbine to obtain the effect of &#x394;<italic>&#x3c9;</italic>
<sub>r</sub>&#x3e;&#x3e;&#x394;<italic>&#x3c9;</italic>
<sub>e</sub>. Due to the coupling problem of multiple variables, the virtual inertia of the wind turbine can be much higher than the inherent inertia in theory, but the dynamic changes in the control process will make it difficult to predict its support performance and seriously weaken its application value.</p>
<p>In the current virtual inertia evaluation method, due to the real-time fluctuation of wind turbine speed and system angular frequency, &#x394;<italic>&#x3c9;</italic>
<sub>r</sub>/&#x394;<italic>&#x3c9;</italic>
<sub>e</sub> needs to be monitored in real-time. In the absence of time-scale constraints, the virtual inertia of the wind turbine will be an uncertain parameter. If reliable inertia reserves cannot be obtained, it will be impossible to predict the drop or increase in the amplitude of frequency after disturbance, which makes the key problem of frequency security early warning difficult to solve.</p>
<p>To solve this problem, a time scale must be introduced to constrain the virtual inertia. In this paper, the extreme time of frequency conversion is taken as the effective time of inertia response of the wind turbine. In the range of <italic>t</italic>
<sub>0</sub>&#x223c;<italic>t</italic>
<sub>
<italic>f</italic>m</sub>, if &#x394;<italic>&#x3c9;</italic>
<sub>e</sub> is equal to the unit value of &#x394;<italic>f</italic>
<sub>max</sub>, then <italic>H</italic>
<sub>vir</sub> can be expressed as<disp-formula id="e6">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#xd7;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mi>R</mml:mi>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>According to (<xref ref-type="disp-formula" rid="e6">6</xref>), &#x394;<italic>&#x3c9;</italic>
<sub>r</sub> determines the size of <italic>H</italic>
<sub>vir</sub> when other parameters are known. At this time, the virtual inertia of the wind turbine can be quantitatively evaluated according to the speed variation in the frequency conversion extreme time. Therefore, the key to evaluating the virtual inertia of the wind turbine is to calculate the extreme value time of the system frequency conversion and obtain the speed change after the inertia response of the wind turbine.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Frequency conversion extreme time of system with high-proportion wind power</title>
<sec id="s3-1">
<title>3.1 Frequency conversion extreme time calculation</title>
<p>Considering the limitation of the inherent inertia of the wind turbine on the speed, as well as the real-time changes in the wind turbine speed and system frequency, not only the potential of virtual inertia control remains unclear, but there are also still many problems in quantitative evaluation. In order to solve this problem, this paper will use the frequency conversion extreme time and the wind turbine operating state to constrain the virtual inertia for quantitative evaluation. Among them, the calculation process of the system frequency conversion extreme time as a key parameter is as follows.</p>
<p>In the system with high-proportion wind power, the aggregation model structure of the synchronous generator set and wind turbine is shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. The model includes the synchronous generator, wind turbine, and load frequency response module. The parameter design of each module is shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>System frequency response aggregation model.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g002.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>System frequency response model parameters and value range.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Model parameter</th>
<th align="center">Value range</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">System inertia time constant <italic>H</italic>
</td>
<td align="center">4&#x223C;10&#xa0;s</td>
</tr>
<tr>
<td align="center">System damping coefficient <italic>D</italic>
</td>
<td align="center">0&#x223C;1</td>
</tr>
<tr>
<td align="center">System load regulation coefficient <italic>K</italic>
<sub>L</sub>
</td>
<td align="center">0&#x223C;2</td>
</tr>
<tr>
<td align="center">Synchronous machine adjustment coefficient <italic>&#x3c3;</italic>
</td>
<td align="center">0.03&#x2013;0.05</td>
</tr>
<tr>
<td align="center">Turbine equivalent inertia time constant <italic>T</italic>
</td>
<td align="center">6&#x223C;10&#xa0;s</td>
</tr>
<tr>
<td align="center">Turbine characteristic coefficient <italic>a</italic>
</td>
<td align="center">0&#x3c;<italic>a</italic>&#x3c;1</td>
</tr>
<tr>
<td align="center">Wind power penetration <italic>k</italic>
</td>
<td align="center">\</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to <xref ref-type="fig" rid="F2">Figure 2</xref>, the system frequency response equation can be expressed as<disp-formula id="e7">
<mml:math id="m7">
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>where &#x394;<italic>P</italic>
<sub>m</sub> is the synchronous generator power response signal; &#x394;<italic>P</italic>
<sub>w</sub> is the wind turbine power response signal; &#x394;<italic>P</italic>
<sub>L</sub> is the load power response signal; and <italic>H</italic> is the system inertia time constant. Among them, &#x394;<italic>P</italic>
<sub>m</sub> depends on the synchronous generator adjustment coefficient and turbine parameters, which can be expressed as<disp-formula id="e8">
<mml:math id="m8">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mfrac>
<mml:mo>&#xd7;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>T</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>where &#x394;<italic>P</italic>
<sub>L</sub> depends on the load regulation coefficient, equal to <italic>K</italic>
<sub>L</sub>&#xd7;&#x394;<italic>f</italic>, and <italic>H</italic> depends on the system capacity and wind power penetration, which can be expressed as<disp-formula id="e9">
<mml:math id="m9">
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where <italic>k</italic> is the wind power penetration and <italic>S</italic>
<sub>B</sub> is the rated capacity of the system.</p>
<p>The system frequency variation can be calculated using <xref ref-type="disp-formula" rid="e7">(7)</xref> and expressed as<disp-formula id="e10">
<mml:math id="m10">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mi>&#x3be;</mml:mi>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mi>b</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="&#x7c;">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>The inverse Laplace transform is applied to (<xref ref-type="disp-formula" rid="e10">10</xref>), and the system frequency response &#x394;<italic>f</italic>(<italic>t</italic>) is<disp-formula id="e11">
<mml:math id="m11">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="&#x7c;">
<mml:mrow>
<mml:mfrac>
<mml:mi>&#x3c3;</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>&#x3be;</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3c6;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e11">(11)</xref>, each parameter is expressed as<disp-formula id="e12">
<mml:math id="m12">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:msqrt>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mi>&#x3be;</mml:mi>
<mml:mi>T</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mi>b</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>T</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>&#x3be;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:msqrt>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>&#x3be;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>T</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="&#x7c;">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mi>w</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>b</mml:mi>
<mml:msqrt>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>&#x3be;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mo>,</mml:mo>
<mml:mtext>and&#x2009;</mml:mtext>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>arctan</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mi>&#x3be;</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>arctan</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mfrac>
<mml:msqrt>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>&#x3be;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3be;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>Without considering the power response of the wind turbine, &#x394;<italic>P</italic>
<sub>w</sub> &#x3d; 0 is taken, and &#x394;<italic>P</italic>
<sub>d</sub> is regarded as a step disturbance. At this time, the derivative of <xref ref-type="disp-formula" rid="e11">(11)</xref> is obtained, and the extreme time of frequency conversion of the system can be expressed as<disp-formula id="e13">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>arctan</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>&#x3be;</mml:mi>
<mml:mi>T</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
<p>According to (<xref ref-type="disp-formula" rid="e11">11</xref>)<xref ref-type="disp-formula" rid="e12">&#x2013;</xref>(<xref ref-type="disp-formula" rid="e13">13</xref>), the extreme time of frequency conversion is not only related to the inertia of conventional synchronous generators but also affected by wind power penetration and wind turbine power response. Among them, conventional synchronous generators can take typical parameter settings. The following section further discusses the influence of wind turbine-related parameters on the frequency conversion extreme time.</p>
</sec>
<sec id="s3-2">
<title>3.2 Influencing factors of frequency conversion extreme time</title>
<sec id="s3-2-1">
<title>3.2.1 Wind power penetration</title>
<p>From <xref ref-type="disp-formula" rid="e8">(8)</xref> and <xref ref-type="disp-formula" rid="e9">(9)</xref>, it can be seen that under the influence of wind power penetration, the system inertia time constant becomes <italic>H</italic>
<sub>g</sub> (1-<italic>k</italic>), and the adjustment coefficient becomes &#x3c3;/(1-<italic>k</italic>); that is, wind power penetration will reduce the system inertia time constant. However, at the same time, the system adjustment coefficient is equivalently increased. When the two parameters change at the same time, the influence of the two on the frequency conversion extreme time has an offset effect. Therefore, the influence of wind power penetration on the frequency conversion extreme time is not significant. In order to verify this conclusion, the typical data in <xref ref-type="table" rid="T1">Table 1</xref> are substituted into (<xref ref-type="disp-formula" rid="e12">12</xref>) and (<xref ref-type="disp-formula" rid="e13">13</xref>) to calculate the intermediate variables and the frequency conversion extreme time, as shown in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Parameter values and calculation results of frequency expression.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th colspan="3" align="center">Parameter value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">
<italic>H</italic> &#x3d; 5&#xa0;s</td>
<td align="center">
<italic>&#x3c3;</italic> &#x3d; 0.04</td>
<td align="center">
<italic>D</italic> &#x3d; 0</td>
</tr>
<tr>
<td align="center">
<italic>T</italic> &#x3d; 8&#xa0;s</td>
<td align="center">
<italic>a</italic> &#x3d; 0.3</td>
<td align="center">
<italic>K</italic>
<sub>L</sub> &#x3d; 1</td>
</tr>
<tr>
<td colspan="3" align="center">
<italic>k</italic> &#x3d; 0% computation</td>
</tr>
<tr>
<td align="center">
<italic>b</italic> &#x3d; 0.570</td>
<td align="center">
<italic>w</italic> &#x3d; 0.299</td>
<td align="center">
<italic>&#x3be;</italic> &#x3d; 0.855</td>
</tr>
<tr>
<td align="center">
<italic>&#x3b1;</italic> &#x3d; 7.218</td>
<td align="center">
<italic>&#x3c6;</italic> &#x3d; &#x2212;8.34<sup>&#x3002;</sup>
</td>
<td align="center">
<italic>t</italic>
<sub>
<italic>f</italic>m</sub> &#x3d; 2.305&#xa0;s</td>
</tr>
<tr>
<td colspan="3" align="center">
<italic>k</italic> &#x3d; 20% computation</td>
</tr>
<tr>
<td align="center">
<italic>b</italic> &#x3d; 0.573</td>
<td align="center">
<italic>w</italic> &#x3d; 0.279</td>
<td align="center">
<italic>&#x3be;</italic> &#x3d; 0.873</td>
</tr>
<tr>
<td align="center">
<italic>&#x3b1;</italic> &#x3d; 7.668</td>
<td align="center">
<italic>&#x3c6;</italic> &#x3d; &#x2212;7.492<sup>&#x3002;</sup>
</td>
<td align="center">
<italic>t</italic>
<sub>
<italic>f</italic>m</sub> &#x3d; 2.298&#xa0;s</td>
</tr>
<tr>
<td colspan="3" align="center">
<italic>k</italic> &#x3d; 40% computation</td>
</tr>
<tr>
<td align="center">
<italic>b</italic> &#x3d; 0.577</td>
<td align="center">
<italic>w</italic> &#x3d; 0.249</td>
<td align="center">
<italic>&#x3be;</italic> &#x3d; 0.902</td>
</tr>
<tr>
<td align="center">
<italic>&#x3b1;</italic> &#x3d; 8.671</td>
<td align="center">
<italic>&#x3c6;</italic> &#x3d; &#x2212;6.619<sup>&#x3002;</sup>
</td>
<td align="center">
<italic>t</italic>
<sub>
<italic>f</italic>m</sub> &#x3d; 2.295&#xa0;s</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The calculation results in <xref ref-type="table" rid="T2">Table 2</xref> are brought into (<xref ref-type="disp-formula" rid="e11">11</xref>), and the power disturbance &#x394;<italic>P</italic>
<sub>d</sub> &#x3d; 0.05&#xa0;pu is taken to obtain the frequency response under different wind power penetration rates, as shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>System frequency response curves under different wind power penetration.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g003.tif"/>
</fig>
<p>It can be seen from <xref ref-type="fig" rid="F3">Figure 3</xref> that the wind power penetration rate will affect the system frequency change rate, the maximum frequency deviation, and the steady-state frequency deviation after primary frequency regulation, but it has little effect on the frequency conversion extreme time. According to the calculation results in <xref ref-type="table" rid="T2">Table 2</xref>, 2.3&#xa0;s is used as the typical value for system frequency conversion extreme time in this paper.</p>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Additional inertial response</title>
<p>In the extreme time of frequency conversion, when the inertia support potential of the wind turbine is stimulated to the maximum extent, the power response of the wind turbine can be regarded as a step change, and the power reference value will jump from the maximum power tracking control signal to the rated power or the minimum power output in a short time. Under this condition, the limit power support signal of the wind turbine can be expressed as<disp-formula id="e14">
<mml:math id="m14">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>where <italic>P</italic>
<sub>n</sub> is the rated power of the wind turbine; <italic>P</italic>
<sub>min</sub> is the power corresponding to the minimum speed of the stable operation of the wind turbine; and <italic>P</italic>
<sub>opt</sub> is the active output power of the wind turbine under MPPT control, which can be expressed as<disp-formula id="e15">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>where <italic>k</italic>
<sub>opt</sub> is the maximum power tracking coefficient and <italic>&#x3c9;</italic>
<sub>0</sub> is the cut-in speed.</p>
<p>When &#x394;<italic>P</italic>
<sub>w</sub> is equivalent to a step signal, its effect is equivalent to the equivalent reduction of &#x394;<italic>P</italic>
<sub>d</sub>, and it does not affect the parameters in <italic>G</italic>(<italic>t</italic>). Therefore, the power response of the wind turbine only affects the increase and decrease in the speed of frequency and the frequency deviation extremum and does not affect the frequency conversion extremum time.</p>
<p>In summary, when evaluating the inertia of the wind turbine in the extreme time of frequency conversion, the calculation results of <xref ref-type="disp-formula" rid="e13">(13)</xref> can be regarded as the action time of the virtual inertia of the wind turbine. At this time, the wind turbine speed variation can be calculated by the wind turbine running state constraint, and then the virtual inertia of the wind turbine can be estimated.</p>
</sec>
</sec>
</sec>
<sec id="s4">
<title>4 Virtual inertia constraint and evaluation of wind turbines</title>
<sec id="s4-1">
<title>4.1 Wind turbine operating state constraints</title>
<p>The virtual inertia of the wind turbine depends on its kinetic energy reserve and is closely related to its initial speed <italic>&#x3c9;</italic>
<sub>r0</sub>. In addition, the initial active power <italic>P</italic>
<sub>we0</sub> of the wind turbine determines the adjustment range of the power support &#x394;<italic>P</italic>
<sub>w</sub>, which in turn affects the kinetic energy release or absorption capacity of the wind turbine.</p>
<sec id="s4-1-1">
<title>4.1.1 Kinetic energy reserve constraint</title>
<p>The wind turbine energy reserve depends on the allowable value of the speed change during the frequency conversion extreme period. When the disturbance power &#x394;<italic>P</italic>
<sub>d</sub> &#x3e;0, the system has a power shortage, and the wind turbine needs to release the rotor kinetic energy; when &#x394;<italic>P</italic>
<sub>d</sub> &#x3c;0, the system has a power surplus, and the wind turbine needs to absorb energy and store it as rotor kinetic energy. Therefore, the kinetic energy reserve constraint of wind turbines can be expressed as Eq. <xref ref-type="disp-formula" rid="e16">(16)</xref>
<disp-formula id="e16">
<mml:math id="m16">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mtext>kw</mml:mtext>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mtext>kwmax</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>J</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>min</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mi>J</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>max</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>where &#x394;<italic>E</italic>
<sub>kwmax</sub> is the wind turbine energy reserve; <italic>&#x3c9;</italic>
<sub>r0</sub> is the wind turbine speed at the beginning of the inertial response; <italic>&#x3c9;</italic>
<sub>min</sub> and <italic>&#x3c9;</italic>
<sub>max</sub> are the allowable minimum and maximum speeds for the stable operation of the wind turbine, respectively, usually 0.7&#xa0;pu and 1.2&#xa0;pu (<xref ref-type="bibr" rid="B4">Fu et al., 2017</xref>). Under the constraint of kinetic energy reserve, the variation range of &#x394;<italic>&#x3c9;</italic>
<sub>r</sub> in the frequency conversion extreme time is expressed as<disp-formula id="e17">
<mml:math id="m17">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>where &#x394;<italic>&#x3c9;</italic>
<sub>rmax1</sub> is the maximum speed variation under the constraint of kinetic energy reserve.</p>
</sec>
<sec id="s4-1-2">
<title>4.1.2 Wind turbine power constraint</title>
<p>In order to ensure the safe operation of the wind turbine, the support power provided by the virtual inertia control is not allowed to exceed the limit. According to (<xref ref-type="disp-formula" rid="e14">14</xref>) and (<xref ref-type="disp-formula" rid="e15">15</xref>), the wind turbine power constraint is<disp-formula id="e18">
<mml:math id="m18">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>wmax</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>min</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>In the process of the inertial response, the rotor motion equation of the wind turbine can be expressed as<disp-formula id="e19">
<mml:math id="m19">
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>wm</mml:mtext>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mtext>we</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>where <italic>P</italic>
<sub>we</sub> is the electromagnetic power output by the wind turbine and <italic>P</italic>
<sub>wm</sub> is the mechanical power captured by the wind turbine, <italic>P</italic>
<sub>wm</sub> &#x3d; <italic>k</italic>
<sub>opt</sub>
<italic>&#x3c9;</italic>
<sub>r0</sub>
<sup>3</sup>.</p>
<p>Combining (<xref ref-type="disp-formula" rid="e18">18</xref>) and (<xref ref-type="disp-formula" rid="e19">19</xref>), under the constraint of the rated power of the wind turbine, the variation range of &#x394;<italic>&#x3c9;</italic>
<sub>r</sub> in the extreme time of frequency conversion is<disp-formula id="e20">
<mml:math id="m20">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mtext>&#x2003;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mtext>opt</mml:mtext>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi>min</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(20)</label>
</disp-formula>where &#x394;<italic>&#x3c9;</italic>
<sub>rmax2</sub> is the maximum speed variation under the wind turbine power constraint.</p>
<p>It can be seen from <xref ref-type="disp-formula" rid="e20">(20)</xref> that the size of &#x394;<italic>&#x3c9;</italic>
<sub>rmax2</sub> is determined using <italic>H</italic>
<sub>w</sub>
<italic>, k</italic>
<sub>opt</sub>, <italic>&#x3c9;</italic>
<sub>r0</sub>, and other parameters. When the capacity and type of the wind turbine are determined, <italic>H</italic>
<sub>w</sub> and <italic>k</italic>
<sub>opt</sub> can be obtained, and &#x394;<italic>&#x3c9;</italic>
<sub>rmax2</sub> only depends on &#x394;<italic>&#x3c9;</italic>
<sub>rmax2</sub>.</p>
<p>In order to satisfy the kinetic energy reserve constraint and the wind turbine power constraint at the same time, according to (<xref ref-type="disp-formula" rid="e17">17</xref>) and (<xref ref-type="disp-formula" rid="e20">20</xref>), the maximum speed variation of the wind turbine &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> can be expressed as Eq. <xref ref-type="disp-formula" rid="e21">(21)</xref>
<disp-formula id="e21">
<mml:math id="m21">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="&#x7c;">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mi>max</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(21)</label>
</disp-formula>
</p>
<p>Taking the typical parameters of a 2&#xa0;MW doubly fed wind turbine as an example, take <italic>H</italic>
<sub>w</sub> &#x3d; 4&#xa0;s, <italic>k</italic>
<sub>opt</sub> &#x3d; 1/1.23, <italic>&#x3c9;</italic>
<sub>min</sub> &#x3d; 0.7&#xa0;pu, <italic>&#x3c9;</italic>
<sub>max</sub> &#x3d; 1.2&#xa0;pu, and <italic>t</italic>
<sub>
<italic>f</italic>m</sub> &#x3d; 2.3&#xa0;s. According to (17), (20), and (21), &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> corresponding to different initial speeds of the wind turbine under kinetic energy reserve constraints and power constraints is obtained as follows.</p>
<p>According to <xref ref-type="fig" rid="F4">Figure 4</xref>, when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e;0, the wind turbine needs to release the kinetic energy of the rotor. When the initial speed of the wind turbine is small, the kinetic energy reserve of the wind turbine is small, but the inertial support power available is large. At this time, &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> is determined using the kinetic energy reserve constraint; when the initial speed of the wind turbine is large, its kinetic energy reserve is large, but the available inertial support power is small. At this time, &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> is determined using the wind turbine power constraint. When the initial speed of the wind turbine is 0.905&#xa0;pu, &#x394;<italic>&#x3c9;</italic>
<sub>rmax1</sub> &#x3d; &#x394;<italic>&#x3c9;</italic>
<sub>rmax2</sub>, and &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> reaches the maximum of 0.205&#xa0;pu at this initial speed. When &#x394;<italic>P</italic>
<sub>d</sub> &#x3c;0, when the wind turbine needs to absorb energy, the relationship between &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> and speed constraint is completely opposite to that when &#x394;<italic>P</italic>
<sub>d</sub> &#x3c;0. When the initial speed of the wind turbine is small, &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> depends on the power constraint of the wind turbine. When the initial speed is large, &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> depends on the kinetic energy reserve constraint. When the initial speed is 1.07&#xa0;pu, &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub> reaches the maximum of 0.13&#xa0;pu.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Evaluation results of the maximum speed variation of the wind turbine.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g004.tif"/>
</fig>
<p>In summary, in the extreme time of frequency conversion, under the constraint of the wind turbine&#x2019;s operating state, when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0, the range of rotational speed variation provided is 0&#x2013;0.205&#xa0;pu; when &#x394;<italic>P</italic>
<sub>d</sub> &#x3c;0, the range of rotational speed variation is 0&#x2013;0.13&#xa0;pu.</p>
</sec>
</sec>
<sec id="s4-2">
<title>4.2 Virtual inertia evaluation of the wind turbine</title>
<p>According to (<xref ref-type="disp-formula" rid="e6">6</xref>), the virtual inertia time constant of the wind turbine also depends on &#x394;<italic>f</italic>
<sub>max</sub>. Since the lower limit of power system frequency security is 48&#xa0;Hz, the frequency deviation extremum of the system should be limited within this range after the inertial response of the wind turbine. Since the wind turbine has a fast power response capability, it can be set within the frequency safety range to complete the predetermined virtual inertia support. Therefore, &#x394;<italic>f</italic>
<sub>max</sub> &#x3d; 2&#xa0;Hz is used to conservatively estimate the virtual inertia of the wind turbine.</p>
<p>According to the evaluation results of <xref ref-type="sec" rid="s4-1">Section 4.1</xref> &#x394;<italic>&#x3c9;</italic>
<sub>rmax</sub>, the range of the virtual inertia time constant of the wind turbine during the frequency conversion extreme time is as follows.</p>
<p>In the diagram, <italic>H</italic>
<sub>virmax1</sub> and <italic>H</italic>
<sub>virmax2</sub> are the maximum inertia time constants of the wind turbine under rotor kinetic energy constraint and power constraint, respectively. When the above two constraints are satisfied simultaneously, <italic>H</italic>
<sub>virmax</sub> can be expressed as Eq. <xref ref-type="disp-formula" rid="e22">(22)</xref>
<disp-formula id="e22">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(22)</label>
</disp-formula>
</p>
<p>It can be seen from <xref ref-type="fig" rid="F5">Figure 5</xref> that when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0 and the speed of the rotor kinetic energy released by the wind turbine decreases, <italic>H</italic>
<sub>virmax</sub> depends on the rotor kinetic energy constraint at low speeds, and <italic>H</italic>
<sub>virmax</sub> depends on the wind turbine power constraint at high speeds. The <italic>H</italic>
<sub>virmax</sub> range is 0&#x2013;18.61&#xa0;s; when &#x394;<italic>P</italic>
<sub>d</sub>&#x3c;0, when the wind turbine absorbs energy and the speed increases, <italic>H</italic>
<sub>virmax</sub> depends on the wind turbine power constraint at low speed, and <italic>H</italic>
<sub>virmax</sub> depends on the rotor kinetic energy constraint at high speed. The <italic>H</italic>
<sub>virmax</sub> range is 0&#x2013;13.85&#xa0;s.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Evaluation results of the maximum virtual inertia of the wind turbine.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g005.tif"/>
</fig>
<p>However, &#x394;<italic>f</italic>
<sub>max</sub> in the actual response process of the system is related to the disturbance power. When the disturbance power is small, &#x394;<italic>f</italic>
<sub>max</sub> will be less than 2&#xa0;Hz, and <italic>H</italic>
<sub>virmax</sub> of the wind turbine will be greater than the evaluation result, which can ensure that the wind turbine has a credible inertial support capability.</p>
<p>According to the above analysis, the wind turbine can show strong inertial support performance in a short time by changing the rotational speed through additional control. In theory, the wind turbine can cope with large power disturbances and has frequency support potential. In the actual control process, a large change in the speed in a short time will bring about problems such as power overshoot, increased mechanical load, and difficulty in speed recovery. Therefore, the wind turbine should perform an inertial response under the premise of considering the system&#x2019;s inertia demand.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Virtual inertia demand and control of wind turbines</title>
<sec id="s5-1">
<title>5.1 Virtual inertia requirements</title>
<p>The inertia of the wind power high-proportion system is significantly reduced, and the virtual inertia of the wind turbine should be set according to the system frequency safety requirements. At present, foreign research has put forward relevant standards for microgrid island operation, requiring that the system frequency change rate (d<italic>f</italic>/d<italic>t</italic>) is not higher than &#xb1;0.5&#xa0;Hz/s (<xref ref-type="bibr" rid="B21">Wang et al., 2022</xref>). The maximum frequency change rate of the system occurs at the initial stage of the power disturbance. The primary frequency modulation has not yet acted, and the unbalanced power of the rotor side of the synchronous machine is the largest. At this time, the system frequency change rate can be expressed as Eq. <xref ref-type="disp-formula" rid="e23">(23)</xref>
<disp-formula id="e23">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(23)</label>
</disp-formula>
</p>
<p>According to the above formula, when the frequency change rate constraint is determined, the system inertia requirement can be expressed as<disp-formula id="e24">
<mml:math id="m24">
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mo>&#x2265;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi>min</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mfenced open="" close="|" separators="&#x7c;">
<mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(24)</label>
</disp-formula>
</p>
<p>where <italic>H</italic>
<sub>min</sub> is the minimum inertia requirement of the system.</p>
<p>When the wind turbine provides virtual inertia, the inertia time constant of the system with a high proportion of wind power can be further expressed using <xref ref-type="disp-formula" rid="e9">(9)</xref> as<disp-formula id="e25">
<mml:math id="m25">
<mml:mrow>
<mml:mi>H</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mi>k</mml:mi>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(25)</label>
</disp-formula>
</p>
<p>Combining (<xref ref-type="disp-formula" rid="e24">24</xref>) and (<xref ref-type="disp-formula" rid="e25">25</xref>), the inertia demand of the wind turbine is expressed as<disp-formula id="e26">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>k</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mo>&#x7c;</mml:mo>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(26)</label>
</disp-formula>
</p>
<p>Taking <italic>H</italic>
<sub>g</sub> &#x3d; 5s and (d<italic>f</italic>/d<italic>t</italic>)<sub>max</sub> &#x3d; &#xb1;0.5&#xa0;Hz/s, the virtual inertia requirement of the wind turbine is calculated using the brought-in <xref ref-type="disp-formula" rid="e26">(26),</xref> as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Wind turbine virtual inertia demand assessment results.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g006.tif"/>
</fig>
<p>It can be seen from <xref ref-type="fig" rid="F6">Figure 6</xref> that the demand of system frequency security for the virtual inertia for the wind turbine mainly depends on disturbance power and wind power penetration. When the disturbance power is small and the wind power penetration rate is low, only the synchronous machine inertia <italic>H</italic>
<sub>g</sub> can meet the system inertia demand. When the permeability is lower than 20%, the inertia of the synchronous machine can cope with more than 8% of the disturbance. However, with the increase in wind power penetration, the ability of the system to cope with power disturbances is gradually weakened. When the wind power penetration reaches 50%, the system can cope with disturbances of less than 5%, which is not enough to cope with typical faults. At this time, the wind turbine should provide the necessary inertial support according to the system requirements. The inertia demand of the system increases with the increase in power disturbance. Under the same power disturbance, the virtual inertia demand of the wind turbine can be calculated according to the wind power penetration rate. When &#x394;<italic>P</italic>
<sub>d</sub> &#x3d; 0.15&#xa0;pu, the system inertia requirement is 7.5&#xa0;s. When <italic>k</italic> &#x3d; 20%, the wind turbine needs to provide 17.5&#xa0;s of virtual inertia to make the total inertia of the system reach 7.5&#xa0;s. When <italic>k</italic> &#x3d; 50%, the wind turbine needs to provide 10&#xa0;s of virtual inertia to make the total inertia of the system reach 7.5&#xa0;s.</p>
</sec>
<sec id="s5-2">
<title>5.2 Virtual inertia control strategy</title>
<p>Combined with the system frequency safety requirements, the wind turbine can use its speed tracking performance to complete the virtual inertial support within the frequency conversion extreme time. The structure of the wind turbine virtual inertia controller proposed in this paper is shown in <xref ref-type="fig" rid="F7">Figure 7</xref>, which consists of two modules: virtual inertia control and virtual inertia evaluation. Among them, the virtual inertia evaluation module has two evaluation functions: the maximum inertia of the wind turbine and the system inertia demand.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Wind turbine virtual inertia controller structure.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g007.tif"/>
</fig>
<p>In the wind turbine virtual inertia evaluation module, the maximum speed variation and maximum inertia are evaluated according to the operating state and frequency conversion extreme time constraint, which are used as the limiting conditions of the wind turbine inertia demand and speed variation demand evaluation results. For the evaluation of the wind turbine&#x2019;s virtual inertia demand, according to (<xref ref-type="disp-formula" rid="e26">26</xref>), disturbance power and wind power permeability can be introduced to calculate the wind turbine&#x2019;s inertia demand. According to the evaluation results of the wind turbine&#x2019;s inertia demand, the wind turbine&#x2019;s speed demand can be calculated. The specific analysis is as follows:</p>
<p>The real-time equivalent inertia time constant <italic>H</italic>
<sub>vir</sub> in the process of wind turbine inertia response can be expressed as<disp-formula id="e27">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(27)</label>
</disp-formula>
</p>
<p>In the extreme time of frequency conversion, both sides of <xref ref-type="disp-formula" rid="e27">(27)</xref> are integrated at the same time to obtain<disp-formula id="e28">
<mml:math id="m28">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn>0</mml:mn>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:msubsup>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(28)</label>
</disp-formula>where &#x394;<italic>f</italic>
<sub>max</sub> is the frequency deviation value corresponding to <italic>t</italic>
<sub>
<italic>f</italic>m</sub>, which can be calculated according to <xref ref-type="disp-formula" rid="e11">(11)</xref> under the assumption that the system inertia meets the demand. According to the calculation results, the simultaneous formulas <xref ref-type="disp-formula" rid="e19">(19)</xref> and <xref ref-type="disp-formula" rid="e28">(28)</xref> are solved, and the rotational speed variation that meets the inertia demand of the wind turbine can be expressed as<disp-formula id="e29">
<mml:math id="m29">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mtext>vir</mml:mtext>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mtext>&#x2002;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(29)</label>
</disp-formula>
</p>
<p>According to the above formula, the required speed variation &#x394;<italic>&#x3c9;</italic>
<sub>r</sub> of the wind turbine at different initial speeds can be obtained by introducing the operating state parameters and inertia requirements of the wind turbine. Then, within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> time, wind turbine corresponding speed can be calculated as <italic>&#x3c9;</italic>
<sub>r1</sub> &#x3d; &#x394;<italic>&#x3c9;</italic>
<sub>r</sub>&#x2b;<italic>&#x3c9;</italic>
<sub>r0</sub>.</p>
<p>The wind turbine virtual inertia control module generates additional power signals through speed tracking control to make the wind turbine track the reference speed in real-time, and the reference speed <italic>&#x3c9;</italic>
<sub>r_ref</sub> is set according to <italic>&#x3c9;</italic>
<sub>r1</sub>. First, in order to ensure that the wind turbine speed can change from <italic>&#x3c9;</italic>
<sub>r0</sub> to <italic>&#x3c9;</italic>
<sub>r1</sub> within <italic>t</italic>
<sub>
<italic>f</italic>m</sub>, the speed-tracking control principle based on the wind turbine rotor motion equation is designed as follows:<disp-formula id="e30">
<mml:math id="m30">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mi>t</mml:mi>
</mml:msubsup>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:msubsup>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(30)</label>
</disp-formula>
</p>
<p>where <italic>K</italic>
<sub>p</sub> is the speed-tracking controller parameter.</p>
<p>According to (<xref ref-type="disp-formula" rid="e30">30</xref>), the value of <italic>K</italic>
<sub>p</sub> in <italic>t</italic>
<sub>
<italic>f</italic>m</sub> can be calculated as<disp-formula id="e31">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(31)</label>
</disp-formula>where <italic>&#x3c9;</italic>
<sub>r_ref</sub> &#x3d; <italic>&#x3c9;</italic>
<sub>r1</sub>&#xb1;&#x394;<italic>&#x3c9;</italic>
<sub>r1</sub> and &#x394;<italic>&#x3c9;</italic>
<sub>r1</sub> is the speed change margin.</p>
<p>According to (<xref ref-type="disp-formula" rid="e31">31</xref>), the value of <italic>K</italic>
<sub>p</sub> can be calculated to realize the change in the wind turbine speed from <italic>&#x3c9;</italic>
<sub>r0</sub> to <italic>&#x3c9;</italic>
<sub>r1</sub> in <italic>t</italic>
<sub>
<italic>f</italic>m</sub> to realize the quantitative control of the wind turbine speed.</p>
<p>Second, in order to ensure the smooth exit of the inertial support power of the wind turbine after the <italic>t</italic>
<sub>
<italic>f</italic>m</sub> moment, the first-order inertial link is introduced, as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. When the wind turbine speed changes to <italic>&#x3c9;</italic>
<sub>r1</sub> at <italic>t</italic>
<sub>0</sub>&#x2b;<italic>t</italic>
<sub>
<italic>f</italic>m</sub>, the switches S<sub>1</sub> and S<sub>2</sub> switch from 0 to 1, and the speed tracking control exits through the first-order inertia link. Under the above control, the wind turbine inertial support power signal can be expressed as<disp-formula id="e32">
<mml:math id="m32">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="&#x7c;">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>_</mml:mo>
<mml:mtext>ref</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(32)</label>
</disp-formula>where <italic>T</italic> is the first-order inertial link time parameter, taking 5&#xa0;s.</p>
<p>Under the control of <xref ref-type="disp-formula" rid="e32">(32)</xref>, the inertial support power of the wind turbine is the largest at the initial moment of disturbance, which can effectively suppress the frequency change rate. After the frequency reaches the extreme point, the rotational speed changes to the demand evaluation value &#x3c9;r1, and the smooth exit of the inertial support power is realized under the action of the first-order inertial link.</p>
</sec>
</sec>
<sec id="s6">
<title>6 Simulation verification</title>
<sec id="s6-1">
<title>6.1 Simulation system</title>
<p>In order to verify the effectiveness of the virtual inertia evaluation and control method of the wind turbine proposed in this paper, the IEEE-9 bus system (<xref ref-type="bibr" rid="B12">None, 2015</xref>) with high-proportion wind power shown in <xref ref-type="fig" rid="F8">Figure 8</xref> is built using DIgSILENT/PowerFactory simulation software. The test system includes three thermal power plants (G1, G2, and G3) and a wind farm (DFIG), which can change the wind power penetration by adjusting the capacity of thermal power units and wind turbines. Assuming that the wind speed remains unchanged at 8&#xa0;m/s, the main parameters of the system are shown in <xref ref-type="table" rid="T3">Table 3</xref>. The load mutation occurs when <italic>t</italic>
<sub>0</sub> &#x3d; 2.0&#xa0;s.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Simulation model of the wind power high permeability system.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g008.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>System parameters.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameters of the wind turbine</th>
<th align="left">Numerical value</th>
<th align="left">Parameters of the wind turbine</th>
<th align="left">Numerical value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Rated voltage/kV</td>
<td align="left">0.69</td>
<td align="left">Frequency bias coefficient</td>
<td align="left">0.04</td>
</tr>
<tr>
<td align="left">Rated power/kVA</td>
<td align="left">2,222</td>
<td align="left">Rated frequency/Hz</td>
<td align="left">50</td>
</tr>
<tr>
<td align="left">Inherent inertia/s</td>
<td align="left">4</td>
<td align="left">Inherent inertia/s</td>
<td align="left">5</td>
</tr>
<tr>
<td align="left">Wind speed/m/s</td>
<td align="left">8</td>
<td align="left"/>
<td align="left"/>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s6-2">
<title>6.2 Frequency conversion extreme time verification</title>
<p>In order to verify the correctness of the frequency conversion extreme time theory calculation, the test system adopts the following parameter settings to analyze the frequency response characteristics after disturbance.</p>
<p>As shown in <xref ref-type="table" rid="T4">Table 4</xref>, in example 1, the wind power penetration rate in the test system is <italic>k</italic> &#x3d; 20%. At the period of 2.0&#xa0;s, the power disturbance is encountered, and the variation is &#x394;<italic>P</italic>
<sub>d</sub> &#x3d; 0.06&#xa0;pu. According to <xref ref-type="fig" rid="F6">Figure 6</xref>, the inherent inertia of the synchronous machine can meet the system frequency safety requirements, and the wind turbine is not required to provide virtual inertia. When the wind turbine is not controlled by virtual inertia, the dynamic response of system frequency and synchronous machine power is shown in <xref ref-type="fig" rid="F9">Figure 9</xref>.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Parameter setting of simulation example 1</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">System parameter</th>
<th align="center">Numerical value</th>
<th align="center">Calculation</th>
<th align="center">Numerical value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Installed gross capacity/MW</td>
<td align="center">1,500</td>
<td align="center">Wind power penetration <italic>k</italic>
</td>
<td align="center">20%</td>
</tr>
<tr>
<td align="center">Total load/MW</td>
<td align="center">900</td>
<td align="center">Power disturbance &#x394;<italic>P</italic>
<sub>d</sub>
</td>
<td align="center">0.06&#xa0;pu</td>
</tr>
<tr>
<td align="center">Wind turbine capacity/MW</td>
<td align="center">300</td>
<td align="center">System inertia requirements <italic>H</italic>
</td>
<td align="center">3&#xa0;s</td>
</tr>
<tr>
<td align="center">Load mutation amount/MW</td>
<td align="center">90</td>
<td align="center">Wind turbine inertia demand <italic>H</italic>
<sub>vir</sub>
</td>
<td align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>
<bold>(a1)</bold>: The change curve of frequency with time &#x394;P<sub>d</sub> &#x3e; 0. <bold>(a2)</bold>: The change curve of PG with time &#x394;P<sub>d</sub> &#x3e; 0. <bold>(b1)</bold>:The change curve of frequency with time &#x394;P<sub>d</sub> &#x3c; 0. <bold>(b2)</bold>:The change curve of PG with time &#x394;P<sub>d</sub> &#x3c; 0.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g009.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F9">Figure 9</xref>, the synchronous machine responds to frequency changes by increasing/decreasing power when the load changes suddenly. The simulation results of <xref ref-type="fig" rid="F9">Figure 9</xref> are compared with the theoretical calculation results, as shown in <xref ref-type="table" rid="T5">Table 5</xref>.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Calculation results of simulation examples.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Parameter</th>
<th align="center">Simulation result</th>
<th align="center">Calculation result</th>
<th align="center">Deviation (%)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Frequency conversion extreme time/s</td>
<td align="center">2.34</td>
<td align="center">2.3</td>
<td align="center">1.74</td>
</tr>
<tr>
<td align="center">Maximum frequency change rate/Hz/s</td>
<td align="center">0.357</td>
<td align="center">0.375</td>
<td align="center">4.7</td>
</tr>
<tr>
<td align="center">Maximum frequency deviation/Hz</td>
<td align="center">0.328</td>
<td align="center">0.318</td>
<td align="center">3.14</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>According to <xref ref-type="table" rid="T5">Table 5</xref>, when the system generator governor is set to a typical parameter, the error between the simulation and calculation results of the system frequency conversion extreme time, the maximum frequency change rate, and the maximum frequency deviation is within 5%. The simulation results verify the accuracy of the calculation and analysis results of the system frequency response model.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s6-3">
<title>6.3 Frequency modulation effect comparison of virtual inertia control</title>
<p>In order to verify the frequency modulation effect of the virtual inertia control strategy proposed in this paper, the power disturbance is increased and compared with the commonly used differential control. The system parameters in example 2 are set as <xref ref-type="table" rid="T6">Table 6</xref>:</p>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Parameter setting of simulation example 2.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">System parameter</th>
<th align="center">Numerical value</th>
<th align="center">Calculation parameter</th>
<th align="center">Numerical value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Installed gross capacity/MW</td>
<td align="center">1,500</td>
<td align="center">Wind power penetration <italic>k</italic>
</td>
<td align="center">20%</td>
</tr>
<tr>
<td align="center">Total load/MW</td>
<td align="center">900</td>
<td align="center">Power disturbance &#x394;<italic>P</italic>
<sub>d</sub>
</td>
<td align="center">0.12&#xa0;pu</td>
</tr>
<tr>
<td align="center">Wind turbine capacity/MW</td>
<td align="center">300</td>
<td align="center">System inertia requirements <italic>H</italic>
</td>
<td align="center">6s</td>
</tr>
<tr>
<td align="center">Load mutation amount/MW</td>
<td align="center">180</td>
<td align="center">Wind inertia demand <italic>H</italic>
<sub>vir</sub>
</td>
<td align="center">10&#xa0;s</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In example 2, the load disturbance variation increases to 180&#xa0;MW. The test system uses formula <xref ref-type="disp-formula" rid="e26">(26)</xref> to analyze the frequency safety, and the demand for the virtual inertia of the wind turbine is 10&#xa0;s. In order to meet the needs of inertial support, the wind turbine adopts the following two control strategies for comparison: 1) when differential control is adopted, <italic>K</italic>
<sub>I</sub> is 20 to meet the inertia demand; 2) when the proposed control is adopted, speed tracking is used. According to the formula <xref ref-type="disp-formula" rid="e27">(27)</xref>, when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0 and &#x394;<italic>P</italic>
<sub>d</sub> &#x3c; 0, the required speed changes in the wind turbine are 0.03&#xa0;pu and 0.029&#xa0;pu, respectively, and the <italic>K</italic>
<sub>p</sub> values obtained using formula <xref ref-type="disp-formula" rid="e30">(30)</xref> are 3.2 and 3.3, respectively. Under the two control strategies, the dynamic response of the system is shown in <xref ref-type="fig" rid="F10">Figures 10</xref>, <xref ref-type="fig" rid="F11">11</xref>.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>
<bold>(A1)</bold>:The first case using differential control. <bold>(A2)</bold>:The second case using differential control. <bold>(A3)</bold>:The third case using differential control. <bold>(B1)</bold>:The first case using speed tracking control. <bold>(B2)</bold>:The second case using speed tracking control. <bold>(B3)</bold>:The third case using speed tracking control.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>
<bold>(A1)</bold>:The first case using differential control. <bold>(A2)</bold>:The second case using differential control. <bold>(A3)</bold>:The third case using differential control. <bold>(B1)</bold>:The first case using speed tracking control. <bold>(B2)</bold>:The second case using speed tracking control. <bold>(B3)</bold>:The third case using speed tracking control.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g011.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F10">Figure 10</xref>, without additional control, the wind turbine hardly responds to frequency changes. The maximum frequency drop amplitude of the system &#x394;<italic>f</italic>
<sub>max</sub> is &#x2212;0.61&#xa0;Hz, and the maximum frequency change rate (d<italic>f</italic>/d<italic>t</italic>)<sub>max</sub> is &#x2212;0.67&#xa0;Hz/s, which has exceeded the allowable value of frequency safety.</p>
<p>When the differential control is adopted, the speed change in the wind turbine is 0.025&#xa0;pu within the frequency conversion extreme time <italic>t</italic>
<sub>
<italic>f</italic>m</sub>, and the maximum frequency change rate (d<italic>f</italic>/d<italic>t</italic>)<sub>max</sub> and the drop amplitude &#x394;<italic>f</italic>
<sub>max</sub> of the system are reduced to &#x2212;0.54&#xa0;Hz/s and &#x2212;0.58&#xa0;Hz, respectively, but still do not meet the maximum frequency change rate constraint. The <italic>H</italic>
<sub>vir</sub> value of the wind turbine in <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is 8.79&#xa0;s, which does not meet the inertia demand of the wind turbine under this example. In addition, since the differential control cannot exit at <italic>t</italic>
<sub>
<italic>f</italic>m</sub>, the overshoot occurs during the frequency recovery process, which is not conducive to frequency security and stability.</p>
<p>When the control strategy proposed in this paper is adopted, the wind turbine speed variation is 0.028&#xa0;pu, the system frequency change rate (d<italic>f</italic>/d<italic>t</italic>)<sub>max</sub> and the drop amplitude &#x394;<italic>f</italic>
<sub>max</sub> are reduced to &#x2212;0.49&#xa0;Hz/s and &#x2212;0.52&#xa0;Hz, respectively, and <italic>H</italic>
<sub>vir</sub> of the wind turbine within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is 10.98&#xa0;s, which meets the system&#x2018;s demand for virtual inertia. According to <xref ref-type="fig" rid="F10">Figure 10</xref>(b1)&#x223c;(b3), the smooth exit of the inertial support power of the wind turbine can be realized in the speed tracking control, which reduces the overshoot in the system frequency recovery process. In addition, under the action of the first-order inertial link, the inertial support power of the wind turbine is reduced to 0 after experiencing delay, and the speed begins to recover later, which avoids the wind turbine absorbing power in the frequency recovery stage and is conducive to system frequency recovery. The frequency modulation effect is better than the traditional differential control.</p>
<p>Similarly, when &#x394;<italic>P</italic>d &#x3c;0, the wind turbine speed variation under differential control is 0.026&#xa0;pu, the system frequency change rate (d<italic>f</italic>/d<italic>t</italic>)<sub>max</sub> and the drop amplitude &#x394;<italic>f</italic>
<sub>max</sub> are reduced to 0.54&#xa0;Hz/s and 0.58&#xa0;Hz, respectively. The <italic>H</italic>
<sub>vir</sub> value of the wind turbine within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is 9.14&#xa0;s, which does not meet the inertia demand of the wind turbine under this example; under the control strategy proposed in this paper, the wind turbine speed variation is 0.029&#xa0;pu, and the system frequency change rate (d<italic>f</italic>/d<italic>t</italic>)<sub>max</sub> and the drop amplitude &#x394;<italic>f</italic>
<sub>max</sub> are reduced to 0.48&#xa0;Hz/s and 0.52&#xa0;Hz, respectively. The <italic>H</italic>
<sub>vir</sub> value of the wind turbine within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is 11.37&#xa0;s, which meets the system&#x2019;s demand for virtual inertia.</p>
<p>In summary, in both cases of &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0 and &#x394;<italic>P</italic>
<sub>d</sub> &#x3c; 0, the proposed control can make the wind turbine show a virtual inertia that meets the system requirements within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> by tracking the speed and can achieve a smooth exit of the inertial support power after <italic>t</italic>
<sub>
<italic>f</italic>m</sub>, avoiding frequency overshoot. The test results show that compared with the current differential control, the inertial support proposed in this paper aims to meet the inertia requirements of the system and combines the virtual inertia constraint of the wind turbine to complete the design of control parameters, which is more conducive to the reliability of the wind turbine to improve frequency support.</p>
</sec>
<sec id="s6-4">
<title>6.4 Evaluation and verification of the maximum inertia of the wind turbine</title>
<p>In order to further test the maximum inertia performance of the wind turbine under the operating state constraints, we will further increase the power disturbance and improve the wind power penetration rate; the system parameters in example 3 are set as <xref ref-type="table" rid="T7">Table 7</xref>:</p>
<table-wrap id="T7" position="float">
<label>TABLE 7</label>
<caption>
<p>Parameter setting of simulation example 3.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">System parameter</th>
<th align="center">Numerical value</th>
<th align="center">Calculation parameter</th>
<th align="center">Numerical value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Installed gross capacity/MW</td>
<td align="center">1,500</td>
<td align="center">Wind power penetration <italic>k</italic>
</td>
<td align="center">30%</td>
</tr>
<tr>
<td align="center">Total load/MW</td>
<td align="center">900</td>
<td align="center">Power disturbance &#x394;<italic>P</italic>
<sub>d</sub>
</td>
<td align="center">0.3&#xa0;pu</td>
</tr>
<tr>
<td align="center">Wind turbine capacity/MW</td>
<td align="center">450</td>
<td align="center">System inertia requirements <italic>H</italic>
</td>
<td align="center">15&#xa0;s</td>
</tr>
<tr>
<td align="center">Load mutation amount/MW</td>
<td align="center">450</td>
<td align="center">Wind inertia demand <italic>H</italic>
<sub>vir</sub>
</td>
<td align="center">50&#xa0;s</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In example 3, the wind power penetration rate is increased to 30%, and the disturbance variation reaches 450&#xa0;MW. In order to ensure that the system frequency change rate meets the constraints, the test system&#x2019;s inertia requirement is set to 15&#xa0;s, and the virtual inertia requirement of the wind turbine is set to 50&#xa0;s. In this case, the inertia demand of the system is far greater than the maximum virtual inertia that the wind turbine can provide, so the wind turbine should provide the maximum inertia support according to the initial operating state during the frequency conversion extreme time. By changing the wind speed, the initial speed of the wind turbine is set to 0.8&#xa0;pu and 1.1&#xa0;pu, and the dynamic response of the system frequency, wind turbine speed, and power output is obtained when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0 and &#x394;<italic>P</italic>
<sub>d</sub> &#x3c; 0, as shown in <xref ref-type="fig" rid="F12">Figures 12</xref>&#x2013;<xref ref-type="fig" rid="F15">15</xref>.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Frequency response when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g012.tif"/>
</fig>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>
<bold>(a1)</bold>:The first case when it is w<sub>ro</sub>&#x3d;1.1. <bold>(a2)</bold>:The second case when it is w<sub>ro</sub>&#x3d;1.1. <bold>(b1)</bold>:The first case when it is w<sub>ro</sub>&#x3d;0.8. <bold>(b2)</bold>:The second case when it is w<sub>ro</sub>&#x3d;0.8.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g013.tif"/>
</fig>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Frequency response when &#x394;P<sub>d</sub> &#x3c; 0.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g014.tif"/>
</fig>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>
<bold>(a1)</bold>:The first case when it is w<sub>ro</sub>&#x3d;1.1. <bold>(a2)</bold>:The second case when it is w<sub>ro</sub>&#x3d;1.1. <bold>(b1)</bold>:The first case when it is w<sub>ro</sub>&#x3d;0.8. <bold>(b2)</bold>:The second case when it is w<sub>ro</sub>&#x3d;0.8.</p>
</caption>
<graphic xlink:href="fenrg-12-1352385-g015.tif"/>
</fig>
<p>When the initial speed of the wind turbine is 1.1&#xa0;pu, the theoretical analysis of <xref ref-type="fig" rid="F4">Figures 4</xref>, <xref ref-type="fig" rid="F5">5</xref> shows that the maximum speed change in the wind turbine in <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is 0.062&#xa0;pu and the maximum inertia is 7.44&#xa0;s. As shown in <xref ref-type="fig" rid="F13">Figure 13</xref>, when &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0, although the kinetic energy of the rotor that the wind turbine can release is large, due to the power constraint, the speed change in the wind turbine during the frequency conversion extreme time is only 0.06&#xa0;pu, which can increase the maximum frequency deviation by 0.3&#xa0;Hz. The maximum inertia <italic>H</italic>
<sub>vir</sub> &#x3d; 7.56&#xa0;s, allows for the smooth exit of wind turbine inertial support after the frequency conversion extreme moment, under the action of the first-order inertial link.</p>
<p>When the initial speed of the wind turbine is <italic>&#x3c9;</italic>
<sub>r0</sub> &#x3d; 0.8&#xa0;pu, the maximum speed change in the wind turbine in <italic>t</italic>
<sub>
<italic>f</italic>m</sub> is 0.1&#xa0;pu and the maximum inertia is 8&#xa0;s. As shown in <xref ref-type="fig" rid="F13">Figure 13</xref>, due to the constraint of the rotor kinetic energy reserve, the inertia response of the wind turbine immediately exits after the speed of the wind turbine drops to the minimum value of 0.7&#xa0;pu at the extreme time of frequency conversion, which will cause a secondary drop in frequency. However, due to the large inertial support power available at this speed, the maximum frequency deviation can be increased by 0.4&#xa0;Hz, and the maximum inertia of the wind turbine during the extreme time of frequency conversion is <italic>H</italic>
<sub>vir</sub> &#x3d; 10&#xa0;s.</p>
<p>Similarly, when &#x394;<italic>P</italic>
<sub>d</sub> &#x3c;0, it can be seen from <xref ref-type="fig" rid="F4">Figures 4</xref>, <xref ref-type="fig" rid="F5">5</xref> that the theoretical values of the maximum speed change in the wind turbine within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> are 0.1&#xa0;pu and 0.034&#xa0;pu at the initial speeds of 1.1&#xa0;pu and 0.8&#xa0;pu, respectively. The maximum inertia is 12&#xa0;s and 2.75&#xa0;s, respectively.</p>
<p>According to the test results in <xref ref-type="fig" rid="F14">Figures 14</xref>, <xref ref-type="fig" rid="F15">15</xref>, when <italic>&#x3c9;</italic>
<sub>r0</sub> &#x3d; 1.1&#xa0;pu, due to the constraint of the rotor kinetic energy reserve, the inertia response of the wind turbine increases to 1.2&#xa0;pu and the frequency continues to increase. However, due to the large inertial support power provided by the wind turbine, the maximum frequency deviation can be reduced by 0.4&#xa0;Hz, and the maximum inertia <italic>H</italic>
<sub>vir</sub> &#x3d; 13.75&#xa0;s. When <italic>&#x3c9;</italic>
<sub>r0</sub> &#x3d; 0.8&#xa0;pu, due to the small inertial support power provided by the wind turbine, the speed change in the wind turbine is 0.035&#xa0;pu in the extreme time of frequency conversion, which can reduce the maximum frequency deviation by 0.2&#xa0;Hz, and the maximum inertia <italic>H</italic>
<sub>vir</sub> &#x3d; 3.5&#xa0;s.</p>
<p>In summary, when the wind turbine needs to provide the maximum inertia support in the face of large disturbances, for the two cases of &#x394;<italic>P</italic>
<sub>d</sub> &#x3e; 0 and &#x394;<italic>P</italic>
<sub>d</sub>&#x3c; 0, the maximum speed variation of the wind turbine within <italic>t</italic>
<sub>
<italic>f</italic>m</sub> under the constraint of the operating state is basically consistent with the theoretical analysis results in this paper. Since &#x394;<italic>f</italic>
<sub>max</sub> is lower than 2&#xa0;Hz under this disturbance, the maximum inertia of the wind turbine is slightly larger than the conservative evaluation result at &#x394;<italic>f</italic>
<sub>max</sub> &#x3d; 2&#xa0;Hz, as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. The simulation test results verify the accuracy of the maximum inertia evaluation analysis method of the wind turbine proposed in this paper.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s7">
<title>7 Conclusion</title>
<p>In order to solve the problem of evaluation and control of the virtual inertia of the wind turbine, this paper proposes the constraint and control method of the virtual inertia in the extreme time of frequency conversion, which is beneficial to improving the frequency active support performance of the wind turbine. The following conclusions are obtained through theoretical analysis and simulation verification:<list list-type="simple">
<list-item>
<p>1) The virtual inertia obtained by the rapid adjustment of the active power of the wind turbine is closely related to many factors, such as frequency change and its own operating state. The dynamic change in virtual inertia in the process of frequency support will greatly increase the difficulty of system frequency safety warning, and it is impossible to accurately predict the control effect.</p>
</list-item>
<list-item>
<p>2) In this paper, the frequency response model of the wind power grid-connected system is established, and the calculated frequency conversion extreme time is introduced into the virtual inertia evaluation. Combined with the wind turbine speed and power output state, the virtual inertia constraint is completed, which can provide more reliable inertia support data for the wind turbine frequency response.</p>
</list-item>
<list-item>
<p>3) Within the allowable range of the system frequency change rate, this paper evaluates the inertia requirements of the system and proposes a virtual inertia control method within the frequency conversion extreme time using the wind turbine speed tracking control performance. The test results show that the proposed control meets the support requirements of the system frequency more reliably and avoids the secondary fluctuation of frequency caused by inertia overshoot.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s8">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s9">
<title>Author contributions</title>
<p>GT: conceptualization, validation, and writing&#x2013;original draft. ZH: data curation and writing&#x2013;review and editing. LB: methodology and writing&#x2013;original draft. MJ: project administration and writing&#x2013;review and editing. FaQ: formal analysis and writing&#x2013;review and editing. XY: resources and writing&#x2013;review and editing. FeQ: investigation and writing&#x2013;original draft.</p>
</sec>
<sec sec-type="funding-information" id="s10">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<sec sec-type="COI-statement" id="s11">
<title>Conflict of interest</title>
<p>Authors GT, ZH, LB, MJ, FaQ, XY, and FeQ were employed by Electric Power Research Institute of Guizhou Power Grid Co., Ltd.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Almeida</surname>
<given-names>R. G.</given-names>
</name>
<name>
<surname>Lopes</surname>
<given-names>J. A. P.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Participation of doubly fed induction wind generators in system frequency regulation</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>22</volume> (<issue>3</issue>), <fpage>944</fpage>&#x2013;<lpage>950</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2007.901096</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Das</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Nu&#xf1;o</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>A. Cutululis</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Frequency stability of power system with large share of wind power under storm conditions</article-title>. <source>J. Mod. Power Syst. Clean Energy</source> <volume>8</volume> (<issue>2</issue>), <fpage>219</fpage>&#x2013;<lpage>228</lpage>. <pub-id pub-id-type="doi">10.35833/mpce.2018.000433</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fern&#xe1;ndez-Guillam&#xf3;n</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>G&#xf3;mez-L&#xe1;zaro</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Muljadi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Molina-Garc&#xed;a</surname>
<given-names>&#xc1;.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Power systems with high renewable energy sources: a review of inertia and frequency control strategies over time</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>115</volume>, <fpage>109369</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2019.109369</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Hei</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Active participation of variable speed wind turbine in inertial and primary frequency regulations</article-title>. <source>Electr. Power Syst. Res.</source> <volume>147</volume>, <fpage>174</fpage>&#x2013;<lpage>184</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2017.03.001</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gautam</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Goel</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Ayyanar</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Vittal</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Harbour</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Control strategy to mitigate the impact of reduced inertia due to doubly fed induction generators on large power systems</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>26</volume> (<issue>1</issue>), <fpage>214</fpage>&#x2013;<lpage>224</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2010.2051690</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Rui</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2016</year>). &#x201c;<article-title>The integrated control strategy for primary frequency control of DFIGs based on virtual inertia and pitch control[C]</article-title>,&#x201d; in <conf-name>Innovative Smart Grid Technologies-ASIA, IEEE</conf-name>, <conf-loc>Melbourne, Australia</conf-loc>, <conf-date>28 November - 1 December 2016</conf-date>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lucas</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Campos-Gaona</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Anaya-Lara</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Assessing the impact of DFIG synthetic inertia provision on power system small-signal stability</article-title>. <source>Energies</source> <volume>12</volume> (<issue>18</issue>), <fpage>3440</fpage>. <pub-id pub-id-type="doi">10.3390/en12183440</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Qiu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Thorp</surname>
<given-names>J. S.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Research on the impact of DFIG virtual inertia control on power system small-signal stability considering the phase-locked loop:</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>32</volume> (<issue>3</issue>), <fpage>2094</fpage>&#x2013;<lpage>2105</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2016.2594781</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Makrini</surname>
<given-names>A. E.</given-names>
</name>
<name>
<surname>Karkri</surname>
<given-names>Y. E.</given-names>
</name>
<name>
<surname>Boukhriss</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>LVRT control strategy of DFIG based wind turbines combining passive and active protection</article-title>. <source>Int. J. Renew. Energy Res.</source> <volume>7</volume> (<issue>3</issue>), <fpage>1258</fpage>&#x2013;<lpage>1269</lpage>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mauricio</surname>
<given-names>J. M.</given-names>
</name>
<name>
<surname>Marano</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Exp&#xf3;sito</surname>
<given-names>A. G.</given-names>
</name>
<name>
<surname>Martinez Ramos</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Frequency regulation contribution through variable-speed wind energy conversion systems</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>24</volume> (<issue>1</issue>), <fpage>173</fpage>&#x2013;<lpage>180</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2008.2009398</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mensou</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Essadki</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Nasser</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bououlid Idrissi</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>A direct power control of a DFIG based WECS during symmetrical voltage dips</article-title>. <source>Prot. Control Mod. Power Syst.</source> <volume>5</volume> (<issue>1</issue>), <fpage>5</fpage>&#x2013;<lpage>47</lpage>. <pub-id pub-id-type="doi">10.1186/s41601-019-0148-y</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<collab>None</collab> (<year>2015</year>). <article-title>Definition and classification of power system stability IEEE/CIGRE joint task force on stability terms and definitions</article-title>. <source>Gold. Res. Thoughts</source> <volume>4</volume> (<issue>3</issue>), <fpage>1387</fpage>&#x2013;<lpage>1401</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ochoa</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Martinez</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Frequency dependent strategy for mitigating wind power fluctuations of a doubly-fed induction generator wind turbine based on virtual inertia control and blade pitch angle regulation</article-title>. <source>Renew. Energy</source> <volume>128</volume> (<issue>PT.A</issue>), <fpage>108</fpage>&#x2013;<lpage>124</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2018.05.047</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qi</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Madonski</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ke</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Tracking-differentiator-based dynamic virtual inertial control of offshore wind power plant for frequency regulation</article-title>. <source>Int. J. Electr. power energy Syst.</source> <volume>141</volume>, <fpage>108150</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2022.108150</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ramtharan</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ekanayake</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Jenkins</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Frequency support from doubly fed induction generator wind turbines</article-title>. <source>IET Renew Power Gener</source> <volume>1</volume> (<issue>1</issue>), <fpage>3</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg:20060019</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>L. B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>L. Z.</given-names>
</name>
<name>
<surname>Ni</surname>
<given-names>Y. X.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Analysis of impact of grid-connected wind power on small signal stability</article-title>. <source>Wind Energy</source> <volume>14</volume> (<issue>4</issue>), <fpage>517</fpage>&#x2013;<lpage>537</lpage>. <pub-id pub-id-type="doi">10.1002/we.440</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2010</year>). &#x201c;<article-title>A review on analysis and control of small signal stability of power systems with large scale integration of wind power</article-title>,&#x201d; in <conf-name>International Conference on Power System Technology, IEEE</conf-name>, <conf-loc>Zhejiang, Hangzhou, China</conf-loc>, <conf-date>24-28 October 2010</conf-date>.</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tielens</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Hertem</surname>
<given-names>D. V.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The relevance of inertia in power systems</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>55</volume>, <fpage>999</fpage>&#x2013;<lpage>1009</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2015.11.016</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ullah</surname>
<given-names>N. R.</given-names>
</name>
<name>
<surname>Thiringer</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Karlsson</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Temporary primary frequency control support by variable speed wind turbines&#x2014; potential and applications</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>23</volume> (<issue>2</issue>), <fpage>601</fpage>&#x2013;<lpage>612</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2008.920076</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vidyanandan</surname>
<given-names>K. V.</given-names>
</name>
<name>
<surname>Senroy</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Primary frequency regulation by deloaded wind turbines using variable droop</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>28</volume> (<issue>2</issue>), <fpage>837</fpage>&#x2013;<lpage>846</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2012.2208233</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wenfeng</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Minimum inertia estimation of power system considering dynamic frequency constraints</article-title>. <source>Proc. CSEE</source> <volume>42</volume> (<issue>1</issue>), <fpage>114</fpage>&#x2013;<lpage>127</lpage>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Minimum inertia requirement assessment of microgrids in islanded and grid-connected modes</article-title>. <source>Proc. CSEE</source> <volume>41</volume> (<issue>6</issue>), <fpage>2040</fpage>&#x2013;<lpage>2053</lpage>.</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Infield</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Power system frequency management challenges-a new approach to assessing the potential of wind capacity to aid system frequency stability</article-title>. <source>IET Renew. Power Gener.</source> <volume>8</volume> (<issue>7</issue>), <fpage>733</fpage>&#x2013;<lpage>739</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2013.0424</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Scholbrock</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Muljadi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Gevorgian</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Evaluation of different inertial control methods for variable-speed wind turbines simulated by fatigue, aerodynamic, structures and turbulence (FAST)</article-title>. <source>IET Renew. Power Gener.</source> <volume>11</volume> (<issue>12</issue>), <fpage>1534</fpage>&#x2013;<lpage>1544</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2017.0123</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ying</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>X. M.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Impact of inertia control of DFIG-based WT on electromechanical oscillation damping of SG</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>33</volume> (<issue>3</issue>), <fpage>3450</fpage>&#x2013;<lpage>3459</lpage>. <pub-id pub-id-type="doi">10.1109/tpwrs.2018.2801283</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Integrated wind turbine controller with virtual inertia and primary frequency responses for grid dynamic frequency support</article-title>. <source>IET Renew. Power Gener.</source> <volume>11</volume> (<issue>8</issue>), <fpage>1129</fpage>&#x2013;<lpage>1137</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2016.0465</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Guan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Che</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Virtual inertia adaptive control of a doubly fed induction generator (DFIG) wind power system with hydrogen energy storage</article-title>. <source>Energies</source> <volume>11</volume> (<issue>4</issue>), <fpage>904</fpage>. <pub-id pub-id-type="doi">10.3390/en11040904</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zertek</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Verbic</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Pantos</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Optimised control approach for frequency-control contribution of variable speed wind turbines</article-title>. <source>IET Renew. Power Gener.</source> <volume>6</volume> (<issue>1</issue>), <fpage>17</fpage>&#x2013;<lpage>23</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2010.0233</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Z. S.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y. Z.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G. J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Coordinated frequency regulation by doubly fed induction generator-based wind power plants</article-title>. <source>IET Renew Power Generation</source> <volume>6</volume> (<issue>1</issue>), <fpage>38</fpage>&#x2013;<lpage>47</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2010.0208</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<etal/>
</person-group> (<year>2017</year>). <article-title>Impact of wind turbine on power oscillation characteristic in region network and integrated control strategy</article-title>. <source>Electr. Meas. Instrum.</source> <volume>54</volume> (<issue>1</issue>), <fpage>33</fpage>&#x2013;<lpage>38</lpage>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>X. H.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Review on power system frequency regulation with high wind power permeability</article-title>. <source>Energy Power Eng.</source> <volume>10</volume> (<issue>8</issue>), <fpage>366</fpage>&#x2013;<lpage>382</lpage>. <pub-id pub-id-type="doi">10.4236/epe.2018.108023</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Comprehensive coordinated control strategy of PMSG-based wind turbine for system inertia support</article-title>. <source>IET Renew. Power Gener.</source> <volume>15</volume> (<issue>9</issue>), <fpage>1915</fpage>&#x2013;<lpage>1926</lpage>. <pub-id pub-id-type="doi">10.1049/rpg2.12115</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>