<?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">739439</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2021.739439</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>Design of Control Strategy and Effect Evaluation for Primary Frequency Regulation of Wind Storage System</article-title>
<alt-title alt-title-type="left-running-head">Li et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Primary Frequency Regulation of Wind/ESS</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Cuiping</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Zheshen</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Junhui</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1394665/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ma</surname>
<given-names>Yunbao</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zou</surname>
<given-names>Jiajun</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<label>
<sup>1</sup>
</label>Key Laboratory of Modern Power System Simulation and Control and Renewable Energy Technology, Ministry of Education (Northeast Electric Power University), <addr-line>Jilin</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<label>
<sup>2</sup>
</label>Wenzhou Power Supply Company, State Grid Zhejiang Electric Power Company, <addr-line>Wenzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<label>
<sup>3</sup>
</label>Zaozhuang Power Supply Company, State Grid Shandong Electric Power Company, <addr-line>Zaozhuang</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<label>
<sup>4</sup>
</label>Jilin Province Electric Power Company of National Grid, <addr-line>Changchun</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/1222560/overview">Bo Yang</ext-link>, Kunming University of Science and Technology, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1409606/overview">Xiaomeng Ai</ext-link>, Huazhong University of Science and Technology, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/957164/overview">Binyu Xiong</ext-link>, Wuhan University of Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Junhui Li, <email>lijunhui@neepu.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Smart Grids, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>09</month>
<year>2021</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>9</volume>
<elocation-id>739439</elocation-id>
<history>
<date date-type="received">
<day>10</day>
<month>07</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>02</day>
<month>08</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2021 Li, Zhang, Li, Ma and Zou.</copyright-statement>
<copyright-year>2021</copyright-year>
<copyright-holder>Li, Zhang, Li, Ma and Zou</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>Because of wind sources&#x2019; volatility and intermittent news, wind power has made a more and more heavy burden on the power system. This paper analyzes in detail the traditional control method, parallel control strategy and serial control strategy of the wind storage system, and combines the advantages of the two to propose an optimal control strategy. A comprehensive performance evaluation method for the primary frequency regulation of the ESS participating in the power grid is proposed based on the power system operation requirements. In the example, the frequency modulation performance of the optimal control strategy is verified by the evaluation method described in this paper in the Chinese frequency adjustment market.</p>
</abstract>
<kwd-group>
<kwd>wind power</kwd>
<kwd>energy storage system</kwd>
<kwd>primary frequency regulation</kwd>
<kwd>control strategy</kwd>
<kwd>effect evaluation</kwd>
</kwd-group>
<contract-num rid="cn001">U1766204</contract-num>
<contract-sponsor id="cn001">Natural Science Foundation of Jilin Province<named-content content-type="fundref-id">10.13039/100007847</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The total installed capacity of wind power in China is rising. As of August 2019, the total installed capacity of wind power in China has reached 198 million kilowatts (<xref ref-type="bibr" rid="B19">Shen et&#x20;al., 2018</xref>; <xref ref-type="bibr" rid="B30">Yan and Zhao, 2018</xref>; <xref ref-type="bibr" rid="B25">Wind power grid-connected, 2019</xref>). Wind power itself has the unfavorable characteristics of large fluctuation and intermittent power output (<xref ref-type="bibr" rid="B18">Ramtharan et&#x20;al., 2007</xref>; <xref ref-type="bibr" rid="B24">Teninge et&#x20;al., 2009</xref>; <xref ref-type="bibr" rid="B15">Muyeen et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B31">Zertek et&#x20;al., 2012</xref>; <xref ref-type="bibr" rid="B4">Dang et&#x20;al., 2014</xref>; <xref ref-type="bibr" rid="B13">Luo et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B17">Peng et&#x20;al., 2015</xref>), After wind power is connected to the grid, it will cause a large burden of frequency modulation on the grid. This burden may cause the power system to collapse (<xref ref-type="bibr" rid="B12">Li et&#x20;al., 2018</xref>; <xref ref-type="bibr" rid="B32">Zhang et&#x20;al., 2018</xref>; <xref ref-type="bibr" rid="B1">Annamraju and Nandiraju, 2019</xref>; <xref ref-type="bibr" rid="B2">Ayyarao and Ayyarao, 2019</xref>; <xref ref-type="bibr" rid="B9">Li et&#x20;al., 2021</xref>). Many studies have been carried out the researches on the influence of wind power on power system frequency.</p>
<p>The following articles mainly solve the problem by improving the frequency regulation capability of wind turbines: In (<xref ref-type="bibr" rid="B21">Sun et&#x20;al., 2018</xref>; <xref ref-type="bibr" rid="B28">Xing et&#x20;al., 2018</xref>), since the wind turbine has no inertia, a control method for reducing the load of the wind turbine is proposed. To increase the reserve capacity of the frequency modulation system, the load capacity of the wind turbines is reduced. In (<xref ref-type="bibr" rid="B27">Wu et&#x20;al., 2013</xref>; <xref ref-type="bibr" rid="B22">Tan et&#x20;al., 2016</xref>), To solve the problem of insufficient capacity of wind turbine frequency modulation reserve, a method for tracking the power curve of sub-optimal wind turbines is proposed. This method guarantees the ability of wind power to participate in frequency modulation by withdrawing from the MPPT mode of the wind turbine. But this method reduces the economy of wind turbines. In summary, the wind turbine&#x2019;s control method for frequency modulation is mainly at the reduced economy or load capacity of wind turbines. Due to its characteristics, the response speed of wind turbines is difficult to meet the needs of power system primary frequency modulation.</p>
<p>The capacity of ESS can be configured as needed, and at the same time, it responds extremely fast to the demand for frequency modulation (<xref ref-type="bibr" rid="B14">Miao et&#x20;al., 2015</xref>; <xref ref-type="bibr" rid="B5">Hu et&#x20;al., 2016</xref>; <xref ref-type="bibr" rid="B11">Li et&#x20;al., 2017</xref>; <xref ref-type="bibr" rid="B23">Tariq et&#x20;al., 2018</xref>; <xref ref-type="bibr" rid="B26">Wu et&#x20;al., 2018</xref>; <xref ref-type="bibr" rid="B16">Pandzic and Bobanac, 2019</xref>; <xref ref-type="bibr" rid="B33">Zhao et&#x20;al., 2019</xref>). ESS can assist wind turbines to solve the problem of frequency modulation.</p>
<p>In the related literature of ESS assisting wind turbine to solve the problem of frequency modulation, in (<xref ref-type="bibr" rid="B7">Jing et&#x20;al., 2005</xref>; <xref ref-type="bibr" rid="B20">Sun et&#x20;al., 2019</xref>), the method of adjusting the frequency in stages solves the problem of economic decline caused by using too much energy storage, and uses the ESS to solve the part of the power system where the frequency changes rapidly. In (<xref ref-type="bibr" rid="B6">Hu et&#x20;al., 2014</xref>), the frequency section control method of wind turbine and the ESS is proposed. This method mainly uses Fourier transform to distribute the high frequency components to the ESS and the remaining components to the wind turbine. In (<xref ref-type="bibr" rid="B29">Yan et&#x20;al., 2016</xref>; <xref ref-type="bibr" rid="B34">Zheng et&#x20;al., 2017</xref>), use the ESS to assist the wind turbine to improve the inertial response capability of the wind turbine, thereby improving the overall frequency regulation capability of the wind storage system.</p>
<p>In order to improve the frequency modulation performance, the above literature is mainly to improve the capability of wind turbine&#x2019;s frequency modulation or use energy storage to improve the frequency modulation performance of the wind turbine. But does not consider how to allocate to the wind turbine and ESS when the frequency modulation instruction is issued. Therefore, this paper proposes an optimal control method based on series control and parallel control based on the related issues of ESS auxiliary wind turbines participating in the primary frequency regulation of the power system. This method can take into account both economical and frequency regulation performance. The effectiveness of the above strategy is verified through simulation cases. Finally, a method is proposed for evaluating the performance of the ESS participating in the primary frequency modulation. The results prove that the optimized control method has the best comprehensive performance in the context of the frequency modulation market in China.</p>
</sec>
<sec id="s2">
<title>Establishment of Mathematical Model Related to Frequency Modulation</title>
<p>In this chapter, the models of traditional generating units participating in the power system, such as hydropower units and thermal power units, use commonly used traditional models (<xref ref-type="bibr" rid="B8">Kundur, 2001</xref>).</p>
<p>Because the selection of relevant parameters of power system frequency modulation is more complicated, and many parameters are empirical parameters, this paper chooses the parameters in documents (<xref ref-type="bibr" rid="B14">Miao et&#x20;al., 2015</xref>) and (<xref ref-type="bibr" rid="B8">Kundur, 2001</xref>) as the main source of this&#x20;paper.</p>
<p>The specific mathematical model is as follows.</p>
<sec id="s2-1">
<title>Mathematical Model of Wind Turbine</title>
<p>At present, the mainstream wind turbines have zero inertia characteristics, so it is necessary to adopt certain technical means to make the wind turbines have the ability of conventional primary frequency modulation. And the main control methods of wind turbine are &#x201c;virtual inertia control&#x201d; and &#x201c;pitch control.&#x201d; Both of them simulate the frequency response of traditional units in power system by changing the active output of wind turbines. The transfer function to build the virtual inertia is shown below: <disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <italic>k</italic>
<sub>
<italic>vd</italic>
</sub> is the virtual inertia constant of the wind turbine, which is usually 8; <italic>T</italic>
<sub>
<italic>wind1</italic>
</sub> is the inertial response time constant of the wind turbine, and the value is usually 0.1&#xa0;s.</p>
<p>The Transfer Function Formula for Pitch Control is as Follows<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <italic>k</italic>
<sub>
<italic>change</italic>
</sub> Is the frequency modulation coefficient for pitch control, usually 20; <italic>T</italic>
<sub>
<italic>wind2</italic>
</sub> is the pitch control time constant of the wind turbine,&#x20;3 s.</p>
<p>In summary, the overall frequency response transfer function formula of the wind turbine is as follows:<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
</sec>
<sec id="s2-2">
<title>Mathematical Model of the ESS</title>
<p>The advantage of the ESS is that its capacity can be configured according to demand, and at the same time can respond to the fast-changing frequency modulation needs. Configuring the ESSs in wind farms can alleviate the pressure of frequency modulation of power systems. Since the ESS is deployed in the wind farm, the inertial response coefficient and the primary frequency modulation coefficient of the ESS are the same as the wind generator. The mathematical model of the ESS is as follows:<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>s</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>where T<sub>Energy</sub> is the frequency response time constant of ESS, generally 0.3&#xa0;s.</p>
</sec>
</sec>
<sec id="s3">
<title>Control Strategy Design</title>
<p>After the power system sends down frequency modulation instructions, the traditional unit and wind storage system will allocate frequency modulation tasks according to how many frequency modulation problems are caused. The specific distribution method is as follows:<disp-formula id="equ1">
<mml:math id="m5">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2217;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m6">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>p</mml:mi>
<mml:mo>&#x2217;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m7">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is primary frequency regulation command responded by the wind storage system, unit MW; <inline-formula id="inf2">
<mml:math id="m8">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is frequency modulation command that the traditional unit responds to, unit MW; p is the proportion of wind power installed capacity in all units participating in frequency modulation, 0.2; P<sub>all</sub> is the total frequency modulation command issued by the power system.</p>
<p>In traditional frequency modulation reserve resources, thermal power units and hydropower units are allocated in the following ways.<disp-formula id="e7">
<mml:math id="m9">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>&#x2217;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math id="m10">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2217;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <inline-formula id="inf3">
<mml:math id="m11">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mtext>P</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>
<sub>wat</sub> and <inline-formula id="inf4">
<mml:math id="m12">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mtext>P</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>
<sub>f</sub> are the regulators of hydroelectric units and thermal power units, respectively. k is the proportion of thermal power units,&#x20;0.8.</p>
<sec id="s3-1">
<title>Control Strategy of Wind-Storage System</title>
<p>The wind turbine and the ESS can be divided into three control modes according to the task assignment when receiving the frequency modulation instruction: serial mode, parallel mode, and optimal operation&#x20;mode.</p>
<sec id="s3-1-1">
<title>Serial Control Mode of Wind-Storage System</title>
<p>Serial control is a classic control strategy of wind/storage system participating in frequency modulation. The schematic diagram is shown in <xref ref-type="fig" rid="F1">Figure&#x20;1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Schematic diagram of wind storage system adopting serial control.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g001.tif"/>
</fig>
<p>In the serial control mode, the frequency modulation work distribution of the wind turbine and the ESS should be carried out according to the following principles,<disp-formula id="e9">
<mml:math id="m13">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
<disp-formula id="e10">
<mml:math id="m14">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>m</mml:mi>
<mml:msup>
<mml:mi>v</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
<disp-formula id="e11">
<mml:math id="m15">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:msubsup>
<mml:mi>V</mml:mi>
<mml:mi>m</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m17">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>where <italic>m</italic> is the quality of wind turbine, kg; v is the speed of wind turbine, m/s; <inline-formula id="inf5">
<mml:math id="m18">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the response time for virtual inertia, s; <inline-formula id="inf6">
<mml:math id="m19">
<mml:mi>&#x3c1;</mml:mi>
</mml:math>
</inline-formula> is the air density, the standard value is 1.29&#xa0;kg/m<sup>3</sup>; <italic>R</italic>
<sub>
<italic>1</italic>
</sub> is the radius of sweep area, m; <italic>V</italic>
<sub>
<italic>m</italic>
</sub> is wind speed, m/s; <italic>C</italic>
<sub>
<italic>p</italic>
</sub> is the percentage of wind energy utilization; <inline-formula id="inf7">
<mml:math id="m20">
<mml:mi>&#x3b2;</mml:mi>
</mml:math>
</inline-formula> is windward pitch angle of wind turbine; <inline-formula id="inf8">
<mml:math id="m21">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the change of upwind pitch angle of wind turbine under serial control; <inline-formula id="inf9">
<mml:math id="m22">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> Is the overall frequency modulation task volume of the wind storage system, MW; <inline-formula id="inf10">
<mml:math id="m23">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the virtual inertial power regulation for wind turbines, MW; <inline-formula id="inf11">
<mml:math id="m24">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the wind turbine pitch control to regulate the active power, MW; <inline-formula id="inf12">
<mml:math id="m25">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf13">
<mml:math id="m26">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the adjustment of the ESS and the wind turbine,&#x20;MW.</p>
<p>The frequency modulation flow chart using serial control is shown in <xref ref-type="fig" rid="F2">Figure&#x20;2</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Flow chart of serial frequency modulation for wind storage combined system.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g002.tif"/>
</fig>
<p>When the wind storage combined system responds to the system frequency modulation instruction by serial control mode, the wind turbine first responds to the frequency modulation instruction of the power system. At this time, the wind turbine controls its output power through virtual inertia and pitch control to suppress the frequency modulation of the power system. The unfinished frequency modulation components of wind turbines are allocated to the&#x20;ESS.</p>
</sec>
<sec id="s3-1-2">
<title>Parallel Control Mode of Wind-Storage System</title>
<p>When a parallel control strategy is adopted to participate in primary frequency modulation of power system, the control strategy sketch is shown in <xref ref-type="fig" rid="F3">Figure&#x20;3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Parallel frequency regulation schematic diagram of wind storage combined system.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g003.tif"/>
</fig>
<p>Where <italic>K</italic>
<sub>1</sub> is the control factor that the wind-storage system adopts a parallel control strategy to allocate tasks when receiving a frequency modulation instruction.</p>
<p>In the parallel control mode, the frequency modulation work distribution of wind turbine and ESS should be carried out according to the following principles,<disp-formula id="e14">
<mml:math id="m27">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">bin</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mtext>&#xa0;</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2217;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="e15">
<mml:math id="m28">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mtext>bin</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>K</mml:mi>
</mml:mrow>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2217;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="italic">P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m29">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">change-</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
<disp-formula id="e17">
<mml:math id="m30">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>m</mml:mi>
<mml:msubsup>
<mml:mi>v</mml:mi>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
<disp-formula id="e18">
<mml:math id="m31">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">change-</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi mathvariant="normal">&#x3d;</mml:mi>
<mml:mfrac>
<mml:mi mathvariant="normal">1</mml:mi>
<mml:mi mathvariant="normal">2</mml:mi>
</mml:mfrac>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>22</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msubsup>
<mml:mi>V</mml:mi>
<mml:mi>m</mml:mi>
<mml:mn>3</mml:mn>
</mml:msubsup>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>p</mml:mi>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>When adopting parallel control strategy, <inline-formula id="inf14">
<mml:math id="m32">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the energy output of the ESS, MW; <italic>v</italic>
<sub>
<italic>1</italic>
</sub> is the speed of wind turbine, m/s; <inline-formula id="inf15">
<mml:math id="m33">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the virtual inertial response time of wind turbines, s; <inline-formula id="inf16">
<mml:math id="m34">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>&#x3001;<inline-formula id="inf17">
<mml:math id="m35">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the changes of pitch angle and pitch angle of wind turbine;<inline-formula id="inf18">
<mml:math id="m36">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>&#x3001;<inline-formula id="inf19">
<mml:math id="m37">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mi mathvariant="normal">&#x2019;</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are wind turbine pitch control and virtual inertia output,&#x20;MW.</p>
<p>The schematic diagram of the process of wind storage system adopting parallel control strategy to participate in grid frequency modulation is shown in <xref ref-type="fig" rid="F4">Figure&#x20;4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Parallel frequency regulation flow chart of wind storage joint System.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g004.tif"/>
</fig>
<p>When wind turbine adopts parallel control, the flow chart of <italic>K</italic>
<sub>
<italic>1</italic>
</sub> value determination is shown in <xref ref-type="fig" rid="F5">Figure&#x20;5</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Flow chart of parallel frequency regulation parameters calculation for wind storage combined system.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g005.tif"/>
</fig>
<p>Considering the actual operation requirements of the power system, the <italic>K</italic>
<sub>
<italic>1</italic>
</sub> value in this paper is taken as&#x20;0.5.</p>
<p>When adopting a parallel control strategy, the wind turbine and the ESS will respond to the primary frequency regulation commands of the power system according to their respective frequency modulation task coefficients <italic>K</italic>
<sub>
<italic>1</italic>
</sub>. Wind turbines change their output through virtual inertia and pitch control. The ESS will respond to frequency modulation instructions quickly within rated power. Wind turbines and ESS s work together to complete the primary frequency modulation task of the power system.</p>
</sec>
<sec id="s3-1-3">
<title>Optimal Control Mode of Wind-Storage System</title>
<p>The optimized control strategy proposed in this paper mainly considers the capacity allocation of the ESS and satisfies the constraints of serial control and parallel control. At each stage of wind power generation participating in power system frequency regulation, the control method is changed according to the actual situation and the operation status of the ESS to meet the power system frequency regulation needs.<list list-type="simple">
<list-item>
<p>1) Target function construction</p>
</list-item>
</list>
</p>
<p>Taking into account the actual operation requirements of the power system, the technical indicators of primary frequency modulation mainly include the maximum frequency deviation and the steady-state frequency deviation of the power system, and the economic indicators mainly include the capacity demand and power demand of the ESS. The objective function of the optimized control strategy of the wind storage system is constructed.<disp-formula id="e19">
<mml:math id="m38">
<mml:mrow>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>where <italic>C</italic>
<sub>
<italic>1</italic>
</sub> and <italic>C</italic>
<sub>
<italic>3</italic>
</sub> are the frequency modulation performance index penalty factor and energy storage cost penalty factor, both 0.5; <inline-formula id="inf20">
<mml:math id="m39">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf21">
<mml:math id="m40">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the maximum frequency deviation and steady-state frequency deviation of the power system, respectively, HZ. <italic>S</italic> is the cost coefficient that changes according to different control strategies.</p>
<p>The frequency modulation flow chart of the optimized control strategy is shown in <xref ref-type="fig" rid="F6">Figure&#x20;6</xref>. Where <italic>K</italic> and <italic>1-K</italic> are the frequency allocation ratio of the ESS using serial control and parallel control, respectively; The system changes the proportion of series control and parallel control at time&#x20;<italic>t</italic>
<sub>
<italic>0</italic>
</sub>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Schematic diagram of optimal control strategy for wind storage system.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g006.tif"/>
</fig>
<p>When the wind storage system adopts the optimal control strategy, the total energy output of the wind storage system is as follows.<disp-formula id="e20">
<mml:math id="m41">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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>&#x2264;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(20)</label>
</disp-formula>where <inline-formula id="inf22">
<mml:math id="m42">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf23">
<mml:math id="m43">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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>&#x2264;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the energy output of the wind-storage combined system at time 0&#x20;&#x3c; t&#x20;&#x3c; t<sub>0</sub> and time t<sub>0</sub> &#x3c;&#x20;t &#x3c; T under the optimized control strategy, respectively, MW; T is the maintenance time of primary frequency modulation, 30&#xa0;s.<disp-formula id="e21">
<mml:math id="m44">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(21)</label>
</disp-formula>
</p>
<p>In the period of 0&#x20;&#x3c; t &#x3c; t<sub>0</sub>, <inline-formula id="inf24">
<mml:math id="m45">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf25">
<mml:math id="m46">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the total output of the serial control and parallel control of the wind storage system under the optimized control strategy, respectively, MW. <inline-formula id="inf26">
<mml:math id="m47">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf27">
<mml:math id="m48">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the output of the serial control and parallel control of the wind turbine under the optimized control strategy, respectively, MW. <inline-formula id="inf28">
<mml:math id="m49">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf29">
<mml:math id="m50">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the output of serial control and parallel control of ESS under optimized control strategy, respectively, MW.<disp-formula id="e22">
<mml:math id="m51">
<mml:mrow>
<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>&#x2264;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<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>&#x2264;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>u</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>e</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(22)</label>
</disp-formula>
</p>
<p>In the formula, the number 2 indicates that the time period is t<sub>0</sub> &#x3c; t &#x3c; T.</p>
<p>The optimized control strategy operation flowchart is shown in <xref ref-type="fig" rid="F7">Figure&#x20;7</xref>.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Flow chart of optimal control strategy for wind storage combined system.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g007.tif"/>
</fig>
<p>The purpose of the optimal control strategy is to reduce the frequency fluctuation of the power system, and reduce the cost of energy storage, and improve the economy.<list list-type="simple">
<list-item>
<p>2) Setting control parameters</p>
</list-item>
</list>
</p>
<p>Calculate <italic>K</italic> and <italic>t</italic>
<sub>
<italic>0</italic>
</sub> with the control variable method. The calculation method is shown in <xref ref-type="fig" rid="F8">Figure&#x20;8</xref>. Where <italic>K</italic>
<sub>
<italic>0</italic>
</sub> is the initial set value; <italic>f</italic>
<sub>min</sub> is the lowest frequency drop point in the simulation process; <italic>f</italic>
<sub>s</sub> is the steady-state frequency value. When the value of <italic>K</italic> is 0.85, the overall frequency modulation effect of the system is the&#x20;best.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Flow chart for the parameters of the optimizing control strategy.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g008.tif"/>
</fig>
<p>The calculation method of <italic>t</italic>
<sub>
<italic>0</italic>
</sub> is shown in <xref ref-type="fig" rid="F9">Figure&#x20;9</xref>. Where <italic>T</italic>
<sub>
<italic>n</italic>
</sub> is the initial time set point. Through simulation verification, when <italic>t</italic>
<sub>
<italic>0</italic>
</sub> takes 0.45&#xa0;s, the frequency modulation effect of the system is the&#x20;best.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Flow chart of time acquisition in optimal control strategy.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g009.tif"/>
</fig>
</sec>
</sec>
<sec id="s3-2">
<title>Demand Analysis of ESS Capacity</title>
<p>This section mainly solves the problem of capacity allocation of the ESSs, including the determination of the constraints of capacity allocation of the ESSs and the establishment of configuration processes.</p>
<sec id="s3-2-1">
<title>Constraints of ESS Configuration</title>
<p>Energy storage systems need to meet the constraints of state of charge (SOC). At a certain time t, SOC is calculated by the following formula:<disp-formula id="e23">
<mml:math id="m52">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x222b;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mi>E</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(23)</label>
</disp-formula>where <italic>t</italic>
<sub>
<italic>1</italic>
</sub> is the time when the ESS starts to work, s; <italic>E</italic> is the ESS rated capacity,&#x20;MWh.</p>
<p>The SOC should meet the constraints of the upper and lower limits:<disp-formula id="e24">
<mml:math id="m53">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:mi>C</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>O</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(24)</label>
</disp-formula>
</p>
<p>
<italic>SOC</italic>
<sub>min</sub> and <italic>SOC</italic>
<sub>max</sub> are the lower and upper limits of charge-discharge of the ESS, respectively, 0.05&#x2013;0.95.</p>
<p>The ESS should also satisfy the output power <italic>P</italic>
<sub>
<italic>e</italic>
</sub> constraint, as follows:<disp-formula id="e25">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(25)</label>
</disp-formula>
</p>
<p>Under normal operation<italic>, P</italic>
<sub>min</sub> and <italic>P</italic>
<sub>max</sub> are the lower limit and upper limit of the output power of the energy storage system, respectively,&#x20;MW.</p>
</sec>
<sec id="s3-2-2">
<title>ESS Configuration Method</title>
<p>The configuration method of the capacity and power of the ESS is as follows:<list list-type="simple">
<list-item>
<p>1) N &#x3d; 1, import the first frequency data;</p>
</list-item>
<list-item>
<p>2) The output data P<sub>1</sub> of the ESS is derived and the rated power P of the ESS is obtained;</p>
</list-item>
<list-item>
<p>3) The capacity E<sub>1</sub> of the ESS is calculated according to the energy output during the frequency modulation period;</p>
</list-item>
<list-item>
<p>4) Considering the charging and discharging efficiency of&#x20;the&#x20;ESS, the rated capacity E of the ESS can be obtained;</p>
</list-item>
<list-item>
<p>5) N &#x3d; N &#x2b;&#x20;1.</p>
</list-item>
</list>
</p>
<p>In summary, the following <xref ref-type="fig" rid="F10">Figure&#x20;10</xref> configuration flowchart can be obtained.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Flow chart of capacity allocation for ESS.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g010.tif"/>
</fig>
<p>The energy storage system configuration results of the three control strategies are shown in <xref ref-type="table" rid="T1">Table&#x20;1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Configuration results of three control strategies</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Control strategy</th>
<th align="center">Rated power (MW)</th>
<th align="center">Rated capacity (MWh)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Parallel control</td>
<td align="char" char=".">6.54</td>
<td align="char" char=".">1.062</td>
</tr>
<tr>
<td align="left">Serial control</td>
<td align="char" char=".">4.84</td>
<td align="char" char=".">0.524</td>
</tr>
<tr>
<td align="left">Optimize control</td>
<td align="char" char=".">5.84</td>
<td align="char" char=".">0.636</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
</sec>
<sec id="s4">
<title>Evaluation Method of Frequency Modulation Effect</title>
<p>The frequency modulation effect evaluation method proposed in this paper comprehensively analyzes the frequency modulation effect and economy and can be used to compare the comprehensive frequency modulation performance of different control strategies.</p>
<sec id="s4-1">
<title>Evaluation System</title>
<p>The calculation formula of comprehensive frequency modulation performance evaluation is as follows:<disp-formula id="e26">
<mml:math id="m55">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2217;</mml:mo>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(26)</label>
</disp-formula>
</p>
<p>
<italic>G</italic> is the score of the final frequency modulation effect; <italic>g</italic>
<sub>
<italic>i</italic>
</sub> is the score of the <italic>i-th</italic> factor affecting the effect of frequency modulation; <italic>D</italic>
<sub>
<italic>i</italic>
</sub> is the weighting factor of the <italic>i-th</italic> factor that affects the frequency modulation effect; <italic>n</italic> is the number of factors that affect the frequency modulation effect.</p>
<p>The calculation method of the maximum frequency deviation <italic>G</italic>
<sub>
<italic>dev</italic>
</sub> of the power grid is as follows:<disp-formula id="e27">
<mml:math id="m56">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2217;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(27)</label>
</disp-formula>
<inline-formula id="inf30">
<mml:math id="m57">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf31">
<mml:math id="m58">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the lower limit frequency and upper limit frequency of the frequency change range <inline-formula id="inf32">
<mml:math id="m59">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, respectively; <inline-formula id="inf33">
<mml:math id="m60">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf34">
<mml:math id="m61">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the lower limit score and upper limit score in the frequency change range <inline-formula id="inf35">
<mml:math id="m62">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, respectively.</p>
<p>The calculation method of power system steady-state frequency deviation <italic>G</italic>
<sub>
<italic>sta</italic>
</sub> is as follows:<disp-formula id="e28">
<mml:math id="m63">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2217;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(28)</label>
</disp-formula>where <inline-formula id="inf36">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf37">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the lower limit frequency and upper limit frequency of the frequency change range <inline-formula id="inf38">
<mml:math id="m66">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, respectively; <inline-formula id="inf39">
<mml:math id="m67">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf40">
<mml:math id="m68">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the lower limit score and upper limit score in the frequency change range <inline-formula id="inf41">
<mml:math id="m69">
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, respectively.</p>
<p>The calculation formula of the power system frequency modulation economy <italic>G</italic>
<sub>
<italic>eco</italic>
</sub> score is as follows:<disp-formula id="e29">
<mml:math id="m70">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>c</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2217;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(29)</label>
</disp-formula>where <inline-formula id="inf42">
<mml:math id="m71">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>c</mml:mi>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf43">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the cost and revenue of energy storage system participating in frequency modulation, respectively; <inline-formula id="inf44">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf45">
<mml:math id="m74">
<mml:mrow>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the lower limit and upper limit of the score in the range, respectively.</p>
</sec>
<sec id="s4-2">
<title>Evaluation Criterion</title>
<p>
<xref ref-type="table" rid="T2">Table&#x20;2</xref> gives the scoring principles for evaluating the comprehensive performance of frequency modulation:</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Score table of frequency regulation effect.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Maximum deviation</th>
<th align="char" char=".">&#x2264;0.2</th>
<th align="char" char=".">0.2&#x223c;0.3</th>
<th align="char" char=".">0.3&#x223c;0.4</th>
<th align="char" char=".">0.4&#x223c;0.5</th>
<th colspan="2" align="char" char=".">&#x2265;0.5</th>
</tr>
<tr>
<th align="left">Score</th>
<th align="char" char=".">10</th>
<th align="char" char=".">8&#x223c;10</th>
<th align="char" char=".">6&#x223c;8</th>
<th align="char" char=".">4&#x223c;6</th>
<th colspan="2" align="char" char=".">&#x2264;3</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Steady deviation</td>
<td align="char" char=".">&#x2264;0.2</td>
<td align="char" char=".">0.2&#x223c;0.25</td>
<td align="char" char=".">0.25&#x223c;0.3</td>
<td align="char" char=".">0.3&#x223c;0.35</td>
<td align="char" char=".">0.35&#x223c;0.4</td>
<td align="char" char=".">0.4</td>
</tr>
<tr>
<td align="left">Score</td>
<td align="center">10</td>
<td align="center">8&#x223c;10</td>
<td align="center">6&#x223c;8</td>
<td align="center">4&#x223c;6</td>
<td align="center">2&#x223c;4</td>
<td align="center">&#x3c;2</td>
</tr>
<tr>
<td align="left">Cost</td>
<td align="center">&#x2264;P</td>
<td align="center">P&#x223c;1.2P</td>
<td align="center">1.2P&#x223c;1.4P</td>
<td align="center">1.4P&#x223c;1.6P</td>
<td align="center">1.6P&#x223c;1.8P</td>
<td align="center">&#x2265;1.8P</td>
</tr>
<tr>
<td align="left">Score</td>
<td align="center">10</td>
<td align="center">8&#x223c;10</td>
<td align="center">6&#x223c;8</td>
<td align="center">4&#x223c;6</td>
<td align="center">2&#x223c;4</td>
<td align="center">&#x2264;2</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Considering that the limit value of power system under-frequency load shedding is 0.5Hz, the lowest frequency limit in the table is set to 0.5&#xa0;Hz; In the case of limit wind power permeability (0.5), the steady frequency of primary frequency modulation is 0.6&#xa0;Hz under the same disturbance, so 0.6&#xa0;Hz is set as the upper limit of steady frequency value. Because there is no compensation in the primary frequency modulation market at present in China when the total benefit-cost ratio of primary frequency modulation in power system is greater than or equal to 1, the score is the highest, and vice versa, the score decreases.</p>
</sec>
</sec>
<sec id="s5">
<title>Simulation Analysis</title>
<sec id="s5-1">
<title>Simulation Parameters</title>
<p>The method described in this paper is based on the windows10 64bit operating system, and the relevant programs are written using matlab 2018b software, Visio2013 is used for flowchart drawing, and Origin2018 software is used for waveform drawing. with the load set to 1000 MW, the wind farm rated power of 200&#xa0;MW, and the load fluctuation of 100&#xa0;MW (0.1p.u).</p>
</sec>
<sec id="s5-2">
<title>Simulation Results With Different Control Strategies</title>
<p>First of all, it is necessary to verify whether the configuration of the wind-storage system has a positive effect on power system frequency modulation&#x3002;When adopting the optimal control strategy, perform the following <xref ref-type="fig" rid="F11">Figure&#x20;11</xref> simulation analysis:</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Frequency regulation results with/without wind power and ESS.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g011.tif"/>
</fig>
<p>The results prove that the configuration of wind/storage system is conducive to power system frequency modulation. <xref ref-type="table" rid="T3">Table&#x20;3</xref> shows the specific values of the impact of the configuration of the wind storage system on frequency, including steady-state frequency deviation and maximum frequency deviation.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Frequency deviation with/without wind power and ESS.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Data type</th>
<th align="center">Maximum deviation (Hz)</th>
<th align="center">Steady deviation (Hz)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">With wind-storage system</td>
<td align="char" char=".">&#x2212;0.425</td>
<td align="char" char=".">&#x2212;0.205</td>
</tr>
<tr>
<td align="left">Without wind-storage system</td>
<td align="char" char=".">&#x2212;0.647</td>
<td align="char" char=".">&#x2212;0.249</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>With the addition of wind storage system, the lowest frequency is increased by 0.22&#xa0;Hz and the steady frequency is increased by 0.044&#xa0;Hz.</p>
<sec id="s5-2-1">
<title>Contribution of Frequency Modulation Resources Under Different Control Strategies</title>
<p>The following is an analysis of two types of frequency modulation resources, including energy storage systems and wind farm.<list list-type="simple">
<list-item>
<p>1) Energy storage system output</p>
</list-item>
</list>
</p>
<p>The primary frequency modulation output of the energy storage system under the same disturbance is simulated using parallel, serial and optimal control strategies. The simulation results are shown in the following <xref ref-type="fig" rid="F12">Figure&#x20;12</xref>.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Comparison of energy storage capacity under different control strategies.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g012.tif"/>
</fig>
<p>In the process of traditional frequency modulation resource startup, the energy storage system responds quickly to frequency changes after disturbances occur in the power system. The energy storage output rises rapidly from 0 to the maximum in 0.3&#xa0;s. The maximum energy storage output is 0.03 p.u under parallel control, 0.023 p.u under serial control and 0.03 p.u under optimal control. Therefore, parallel control and optimal control of the energy storage system can provide greater energy support at the early stage of disturbance occurrence and raise the frequency to the lowest&#x20;point.</p>
<p>The use of energy storage resources in parallel control is stable at 0.015 p.u., the series control is almost 0, and the optimal control is between the two, stable at 0.0047 p.u. The optimized control of the energy storage output power is small, that is, the use of energy storage resources is small, and it can provide backup capacity for the next frequency modulation work.<list list-type="simple">
<list-item>
<p>2) Wind farm output</p>
</list-item>
</list>
</p>
<p>Wind farms are an important part of the frequency modulation resources of power systems. The control strategy should consider improving the frequency modulation performance of wind farms. The following <xref ref-type="fig" rid="F13">Figure&#x20;13</xref> is an analysis of the wind farm frequency modulation capabilities of three different control strategies.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Comparison of wind power under different control strategies.</p>
</caption>
<graphic xlink:href="fenrg-09-739439-g013.tif"/>
</fig>
<p>In the start-up stage of wind turbine, the virtual inertia response of wind farm under each control strategy releases the kinetic energy of wind turbine rotor within 0.1&#xa0;s. This will delay the fast frequency drop of power system, and each control strategy has little influence on this process.</p>
<p>In the process of pitch control, because wind power takes on different tasks in different control strategies, the control effect of the control strategy is reflected in the process of pitch output. Serial control mainly uses wind power as the backup energy for ESS, so the maximum output power of wind power is 0.024 p. u. In parallel control, wind power and energy storage work in parallel mode, so the output power of wind power is smaller than 0.017 p. u., and the optimal control strategy is 0.021p.u. between the above two. This proves that the optimized control strategy can reasonably use the capacity of wind power frequency modulation and reduce the amount of energy storage action.</p>
</sec>
<sec id="s5-2-2">
<title>Frequency Modulation Effect Comparison</title>
<p>The following <xref ref-type="sec" rid="s13">Supplementary Figure S1</xref> is a simulation analysis of the frequency modulation effects of the three control strategies using Simulink.</p>
<p>
<xref ref-type="table" rid="T4">Table&#x20;4</xref> is the frequency modulation effect data of three control strategies, including maximum frequency deviation and steady-state frequency deviation.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Frequency deviation comparison.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Control strategy</th>
<th align="center">Maximum deviation (Hz)</th>
<th align="center">Steady deviation (Hz)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Parallel control</td>
<td align="char" char=".">&#x2212;0.387</td>
<td align="char" char=".">&#x2212;0.187</td>
</tr>
<tr>
<td align="left">Serial control</td>
<td align="char" char=".">&#x2212;0.442</td>
<td align="char" char=".">&#x2212;0.241</td>
</tr>
<tr>
<td align="left">Optimize control</td>
<td align="char" char=".">&#x2212;0.401</td>
<td align="char" char=".">&#x2212;0.205</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>It can be drawn from the table above that the maximum frequency deviation and the steady-state frequency deviation of the optimized control strategy are both between serial control and parallel control, and are superior to serial control.</p>
<p>Generally speaking, the optimal control strategy can provide effective support for the power system by utilizing the fast response characteristics of the energy storage system in the start-up stage of traditional units and wind power. After the start-up of the unit, on the premise of maximizing the frequency modulation capability of the traditional unit, the energy storage capacity is reduced and the economy of frequency modulation is guaranteed.</p>
</sec>
<sec id="s5-2-3">
<title>Cost Analysis Under Different Control Strategies</title>
<p>Taking the lithium batteries which are widely used in frequency modulation of power system as the research object, the economy of lithium batteries under three control strategies is analyzed. The simulation parameters are shown in the following&#x20;<xref ref-type="table" rid="T5">Table&#x20;5</xref>.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Economy parameters of frequency regulation.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameter name</th>
<th align="center">Recovery cost</th>
<th align="center">Abandoned wind income</th>
<th align="center">Unit power price</th>
<th align="center">Unit capacity price</th>
<th align="center">Time (s)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Parameter values</td>
<td align="center">4,046&#xa0;$/t</td>
<td align="center">80&#xa0;$/MWh</td>
<td align="center">738,462&#xa0;$/MW</td>
<td align="center">215,385&#xa0;$/MWh</td>
<td align="center">30</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>There are no compensation measures for battery energy storage participating in primary frequency modulation of power grid (<xref ref-type="bibr" rid="B3">Chen et&#x20;al., 2016</xref>). But battery energy storage can be recycled. For example, lithium batteries can recycle non-ferrous metals for secondary use (<xref ref-type="bibr" rid="B10">Li et&#x20;al</xref>). At the same time, the use of energy storage systems can also reduce wind curtailment. <xref ref-type="table" rid="T6">Table&#x20;6</xref> compares the costs and profits of energy storage systems.</p>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Cost comparison under different control strategies.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Control strategy</th>
<th align="center">Profit ($)</th>
<th align="center">Cost ($)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Parallel Control</td>
<td align="center">1,428,154</td>
<td align="center">5,058,277</td>
</tr>
<tr>
<td align="left">Serial Control</td>
<td align="center">1,026,262</td>
<td align="center">3,687,015</td>
</tr>
<tr>
<td align="left">Optimize Control</td>
<td align="center">1,177,108</td>
<td align="center">4,449,600</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>From the above data, it can be concluded that the cost of the parallel control strategy under the same interference is the highest. Optimal control reduces the energy storage workload after the traditional frequency modulation resources are started, so the economy is better. Serial control has lower revenue due to less initial power and capacity.</p>
<p>According to the evaluation method proposed in Chapter 4, the three control strategies are scored and compared. The scores are shown in <xref ref-type="table" rid="T7">Table&#x20;7</xref>.</p>
<table-wrap id="T7" position="float">
<label>TABLE 7</label>
<caption>
<p>Scoring table for different control strategies.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Control strategy</th>
<th align="center">Lowest frequency score</th>
<th align="center">Frequency steady score</th>
<th align="center">Cost score</th>
<th align="center">Total score</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Parallel control</td>
<td align="char" char=".">7.2</td>
<td align="center">10</td>
<td align="char" char=".">4.255</td>
<td align="char" char=".">6.427</td>
</tr>
<tr>
<td align="left">Serial control</td>
<td align="char" char=".">5.95</td>
<td align="char" char=".">9.35</td>
<td align="char" char=".">5.59</td>
<td align="char" char=".">6.72</td>
</tr>
<tr>
<td align="left">Optimize control</td>
<td align="char" char=".">7.05</td>
<td align="char" char=".">9.75</td>
<td align="char" char=".">5.92</td>
<td align="char" char=".">7.122</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to the results in the above table, the overall score of optimized control is higher than that of serial control and parallel control, which proves that the overall performance of optimized control is higher than that of serial control and parallel control after a comprehensive analysis of the economy and frequency modulation effect.</p>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s6">
<title>Conclusion</title>
<p>Aiming at the control strategy of large-scale wind storage system participating in power system primary frequency modulation, this paper made the following work:</p>
<p>In this paper, the serial and parallel control modes of ESS and wind turbine are studied based on the traditional frequency modulation mode of power system. Considering the economy and frequency modulation effect, an optimal control strategy based on serial and parallel of wind-storage system is proposed, and the objective function of the optimal control strategy is established. And the control variable method is used to optimize the relevant variables in the optimization strategy.</p>
<p>The comprehensive evaluation method of wind storage combined system participating in primary frequency regulation of power grid is studied. Considering that the current primary frequency modulation evaluation method cannot meet the requirement of paying equal attention to both the effect of frequency modulation and the economy of frequency modulation in the power grid. Considering the effect of frequency modulation under different control strategies and the economic factors of wind storage combined system participating in primary frequency modulation, a comprehensive evaluation method of primary frequency modulation are designed.</p>
<p>In the analysis part of the calculation example, the same interference is added to the three different control strategies, and the frequency modulation effect is analyzed using the simulation system. The results show that the optimized control reduces the maximum frequency deviation by 0.041&#xa0;Hz and the steady-state frequency deviation by 0.036&#xa0;Hz than the serial control. Compared with the parallel control strategy, the rated power is reduced by 0.7&#xa0;MW and the rated capacity is reduced by 0.462&#xa0;MWh.</p>
<p>According to the proposed evaluation method, the comprehensive frequency modulation effect scores of the three control strategies are evaluated. The results show that the optimal control strategy proposed in this paper scores 7.122 points, which is higher than 6.72 points for serial control and 6.427 points for parallel control. This proves that the overall performance of optimized control is higher than that of serial control and parallel control after a comprehensive analysis of the economy and frequency modulation effect.</p>
</sec>
<sec id="s7">
<title>Outlook</title>
<p>This paper optimizes the parallel and serial control in the traditional frequency modulation method of wind storage system, and puts forward the optimal control strategy of this paper, and finally achieved certain results. At the same time, this paper still has some areas that can be improved. For example, the effect of the optimal control strategy proposed in this paper may not be better than the serial control strategy and the parallel control strategy under certain conditions. The optimal control strategy is the optimal method proposed after comprehensively considering the cost of energy storage and the effect of frequency modulation. It may not be applicable in some occasions where the accuracy of frequency modulation needs to be improved (while the cost of frequency modulation is allowed to be ignored). In fact, the cost of the current energy storage system is still high, so its economy is still a very important indicator in the power system. If the cost of the energy storage system can be significantly reduced in the future, then this strategy still has room for continuous improvement.</p>
</sec>
</body>
<back>
<sec id="s8">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s13">Supplementary Material</xref>, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s9">
<title>Author Contributions</title>
<p>CL: Conceptualization, Methodology, Writing- Original draft preparation. ZZ: Writing- Original draft preparation. JL: Methodology. YM: Data curation. JZ: Translation and Validation.</p>
</sec>
<sec id="s10">
<title>Funding</title>
<p>This work was supported in part by Joint Foundation of Natural Science Foundation of Jilin Province (No. 2020122352JC) and the National Natural Science Foundation of China (No. U1766204).</p>
</sec>
<sec sec-type="COI-statement" id="s11">
<title>Conflict of Interest</title>
<p>Author ZZ was employed by the company Wenzhou Power Supply Company, State Grid Zhejiang Electric Power Company. Author YM was employed by the company Zaozhuang Power Supply Company, State Grid Shandong Electric Power Company. Author JZ was employed by the company Jilin Province Electric Power Company of National&#x20;Grid.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s12" sec-type="disclaimer">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s13">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2021.739439/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2021.739439/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material>
<label>SUPPLEMENTARY FIGURE S1</label>
<caption>
<p>Comparison chart of frequency modulation effect.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image1.png" id="SM1" mimetype="application/png" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Annamraju</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Nandiraju</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Robust Frequency Control in a Renewable Penetrated Power System: an Adaptive Fractional Order-Fuzzy Approach</article-title>. <source>Prot. Control. Mod. Power Syst.</source> <volume>4</volume> (<issue>4</issue>), <fpage>181</fpage>&#x2013;<lpage>195</lpage>. <pub-id pub-id-type="doi">10.1186/s41601-019-0130-8</pub-id> </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ayyarao</surname>
<given-names>T. S. L. V.</given-names>
</name>
<name>
<surname>Ayyarao</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Modified Vector Controlled DFIG Wind Energy System Based on Barrier Function Adaptive Sliding Mode Control</article-title>. <source>Prot. Control. Mod. Power Syst.</source> <volume>4</volume> (<issue>4</issue>), <fpage>34</fpage>&#x2013;<lpage>41</lpage>. <pub-id pub-id-type="doi">10.1186/s41601-019-0119-3</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Dayu.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Lizi.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Liguo.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Control Strategy of Energy System for Frequency Regulation and Evaluation of Investment Income</article-title>. <source>Mod. Electric Power</source> <volume>33</volume> (<issue>1</issue>), <fpage>80</fpage>&#x2013;<lpage>86</lpage>. <pub-id pub-id-type="doi">10.19725/j.cnki.1007-2322.2016.01.013</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Seuss</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Suneja</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Harley</surname>
<given-names>R. G.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>SoC Feedback Control for Wind and ESS Hybrid Power System Frequency Regulation</article-title>. <source>IEEE J.&#x20;Emerg. Sel. Top. Power Electron.</source> <volume>2</volume> (<issue>1</issue>), <fpage>79</fpage>&#x2013;<lpage>86</lpage>. <pub-id pub-id-type="doi">10.1109/jestpe.2013.2289991</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>Zechun.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>Rui.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Linlin.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Joint Operation Optimization of Wind-Storage Union with Energy Storage Participating Frequency Regulation</article-title>. <source>Power Syst. Tech.</source> <volume>40</volume> (<issue>08</issue>), <fpage>2251</fpage>&#x2013;<lpage>2257</lpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2016.08.001</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>Zechun.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X. I. E.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Fang.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Research on Automatic Generation Control Strategy Incorporating Energy Storage Resources</article-title>. <source>Proc. CSEE</source> <volume>34</volume> (<issue>29</issue>), <fpage>5080</fpage>&#x2013;<lpage>5087</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.2014.29.009</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jing</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>Yong.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>Jiahua.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Research on Subsection and Layer Control Strategy of Doublyfed Variable Speed Wind Turbine</article-title>. <source>Power Syst. Tech.</source> <volume>29</volume> (<issue>9</issue>), <fpage>15</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2005.09.004</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kundur</surname>
<given-names>Prabha.</given-names>
</name>
</person-group> (<year>2001</year>). <source>Power System Stability and Control</source>. <publisher-loc>Beijing</publisher-loc>: <publisher-name>China Electric Power Press</publisher-name>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Xing</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Improving Wind Power Integration by Regenerative Electric Boiler and Battery Energy Storage Device</article-title>. <source>Int. J.&#x20;Electr. Power Energ. Syst.</source> <volume>131</volume> (<issue>131</issue>), <fpage>107039</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2021.107039</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="other">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Jianbo.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>X. U.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>Zhongguang.</given-names>
</name>
</person-group> <article-title>An Overview on the Current Technologies of Recycling Spent Lithium-Ion Bat-Teries</article-title>. <source>Chin. J.&#x20;Rare Met.</source> <volume>43</volume> (<issue>02</issue>), <fpage>201</fpage>&#x2013;<lpage>212</lpage>. <pub-id pub-id-type="doi">10.13373/j.cnki.cjrm.xy17060021</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Jianlin.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Huimeng.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>Xiaodong.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Overview on Key Applied Technologies of Large-Scale Distributed Energy Storage</article-title>. <source>Power Syst. Tech.</source> <volume>41</volume> (<issue>10</issue>), <fpage>3365</fpage>&#x2013;<lpage>3375</lpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2017.1101</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Junhui.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Xichao.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Gangui.</given-names>
</name>
<name>
<surname>Ge</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Survey on Frequency Regulation Technology in High Wind Penetration Power System</article-title>. <source>Power Syst. Prot. Control.</source> <volume>46</volume>, <fpage>163</fpage>&#x2013;<lpage>170</lpage>. <pub-id pub-id-type="doi">10.7667/PSPC170002</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>K. P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Coordinated Operational Planning for Wind Farm with Battery Energy Storage System</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>6</volume> (<issue>1</issue>), <fpage>253</fpage>&#x2013;<lpage>262</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2014.2367550</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Miao</surname>
<given-names>Fufeng.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>Xisheng.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>Zhiping.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Analysis of Frequency Characteristics of Power System Based on Wind Farm-Energy Storage Combined Frequency Regulation</article-title>. <source>High Voltage Eng.</source> <volume>41</volume> (<issue>07</issue>), <fpage>2209</fpage>&#x2013;<lpage>2216</lpage>. </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Muyeen</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Hasanien</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Tamura</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Reduction of Frequency Fluctuation for Wind Farm Connected Power Systems by an Adaptive Artificial Neural Network Controlled Energy Capacitor System</article-title>. <source>IET Renew. Power Gener.</source> <volume>6</volume> (<issue>4</issue>), <fpage>226</fpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg.2010.0126</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pandzic</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Bobanac</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>An Accurate Charging Model of Battery Energy Storage</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>34</volume> (<issue>2</issue>), <fpage>1416</fpage>&#x2013;<lpage>1426</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2018.2876466</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Dou</surname>
<given-names>Pengchong.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Yuwei.</given-names>
</name>
<name>
<surname>Fengquan</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Application of Microgrid Technology in Active Distribution Network</article-title>. <source>Electric Power Automation Equipment</source> <volume>35</volume> (<issue>4</issue>), <fpage>8</fpage>&#x2013;<lpage>16</lpage>. <pub-id pub-id-type="doi">10.16081/j.issn.1006-6047.2015.04.002</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ramtharan</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Jenkins</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Ekanayake</surname>
<given-names>J.&#x20;B.</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>0</lpage>. <pub-id pub-id-type="doi">10.1049/iet-rpg:20060019</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shen</surname>
<given-names>Guanye.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Bingliang.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Economic Allocation for Energy Storage System Considering Wind Power</article-title>. <source>J.&#x20;Northeast Electric Power Univ.</source> <volume>38</volume> (<issue>04</issue>), <fpage>27</fpage>&#x2013;<lpage>34</lpage>. </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>Bohao.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>Yong.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>Y. E.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Integrated Frequency Control Strategy for Wind Power Cluster with Mul-Tiple Temporal-Spatial Scale Coordination Based on H-DMPC</article-title>. <source>Proc. CSEE</source> <volume>39</volume> (<issue>1</issue>), <fpage>155</fpage>&#x2013;<lpage>167</lpage>. <pub-id pub-id-type="doi">10.1109/tec.2015.2476827</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>Yi.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Zekun.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Di.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>An Optimization Control Strategy for Intelligent UFLS Using Dynamically Aggregated Controllable Loads</article-title>. <source>Proc. CSEE</source> <volume>38</volume> (<issue>7</issue>), <fpage>1913</fpage>&#x2013;<lpage>1921&#x2b;2207</lpage>. <pub-id pub-id-type="doi">10.13334/j.0258-8013.pcsee.1707407</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Meegahapola</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Muttaqi</surname>
<given-names>K. M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A Suboptimal Power-Point-Tracking-Based Primary Frequency Response Strategy for DFIGs in Hybrid Remote Area Power Supply Systems</article-title>. <source>IEEE Trans. Energ. Convers.</source> <volume>31</volume> (<issue>1</issue>), <fpage>93</fpage>&#x2013;<lpage>105</lpage>. <pub-id pub-id-type="doi">10.1109/tec.2015.2476827</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tariq</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Maswood</surname>
<given-names>A. I.</given-names>
</name>
<name>
<surname>Gajanayake</surname>
<given-names>C. J.</given-names>
</name>
<name>
<surname>Gupta</surname>
<given-names>A. K.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Modeling and Integration of a Lithium-Ion Battery Energy Storage System with the More Electric Aircraft 270&#x20;V DC Power Distribution Architecture</article-title>. <source>IEEE Access</source> <volume>6</volume>, <fpage>41785</fpage>&#x2013;<lpage>41802</lpage>. <pub-id pub-id-type="doi">10.1109/ACCESS.2018.2860679</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Teninge</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Jecu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Roye</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Bacha</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Duval</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Belhomme</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Contribution to Frequency Control through Wind Turbine Inertial Energy Storage</article-title>. <source>IET Renew. Power Gener.</source> <volume>3</volume> (<issue>3</issue>), <fpage>358</fpage> <pub-id pub-id-type="doi">10.1049/iet-rpg.2008.0078</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="book">
<collab>Wind power grid-connected operation in the first three quarters</collab> (<year>2019</year>). <source>Wind Power Grid-Connected Operation in the First Three Quarters of 2019</source>. <publisher-loc>Beijing</publisher-loc>: <publisher-name>National Energy Administration</publisher-name>.</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Jun.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Zongxiang.</given-names>
</name>
<name>
<surname>Qiao</surname>
<given-names>Ying.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Optimal Operation of Wind Farm with Hybrid Storage Devices Considering Efficiency Characteristics of Dynamic Charging and Discharging</article-title>. <source>Automation Electric Power Syst.</source> <volume>42</volume> (<issue>11</issue>), <fpage>41</fpage>&#x2013;<lpage>47</lpage>. <pub-id pub-id-type="doi">10.7500/AEPS20180326002</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Zishuang.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>Jilai.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>Xiyun.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>DFIG&#x2019;s Frequency Regulation Method Only for High Wind Speed with Subop-Timal Power Tracking</article-title>. <source>Trans. China Electrotechnical Soc.</source> <volume>28</volume> (<issue>05</issue>), <fpage>112</fpage>&#x2013;<lpage>119</lpage>. </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xing</surname>
<given-names>Pengxiang.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>Qiaoming.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Gang.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Response Characteristics and Mechanism Analysis about Virtual Inertia Control of Wind Generators</article-title>. <source>High Voltage Eng.</source> <volume>44</volume> (<issue>04</issue>), <fpage>1302</fpage>&#x2013;<lpage>1310</lpage>. <pub-id pub-id-type="doi">10.13336/j.1003-6520.hve.20180329033</pub-id> </citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>Gangui.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Yibo.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>Z. H. O. N. G.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Frequency Control Strategy for Wind Storage Combined System</article-title>. <source>Electric Power Construction</source> <volume>37</volume> (<issue>12</issue>), <fpage>55</fpage>&#x2013;<lpage>60</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-7229.2016.12.007</pub-id> </citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>Gangui.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Weizhe.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>ZHANG Liyu Research on Integrated Control of Deloading Frequency Regulation for Variable Speed and Variable Pitch Angle Wind Turbines</article-title>. <source>J.&#x20;Northeast Electric Power Univ.</source> <volume>38</volume> (<issue>05</issue>), <fpage>1</fpage>&#x2013;<lpage>8</lpage>. </citation>
</ref>
<ref id="B31">
<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>A Novel Strategy for Variable-Speed Wind Turbines&#x27; Participation in Primary Frequency Control</article-title>. <source>IEEE Trans. Sustain. Energ.</source> <volume>3</volume> (<issue>4</issue>), <fpage>791</fpage>&#x2013;<lpage>799</lpage>. <pub-id pub-id-type="doi">10.1109/tste.2012.2199773</pub-id> </citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Xu.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Yunlong.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>Shuai.</given-names>
</name>
<name>
<surname>Zha</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Retrospect and Prospect of Research on Frequency Regulation Technol-Ogy of Power System by Wind Power</article-title>. <source>Power Syst. Tech.</source> <volume>42</volume> (<issue>06</issue>), <fpage>1793</fpage>&#x2013;<lpage>1803</lpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2018.0359</pub-id> </citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>Jingjing.</given-names>
</name>
<name>
<surname>Min</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Xinqin.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Coordinated Control Strategy of Wind Power and Energy Storage in Frequency Regulation Based on Torque Limit Control</article-title>. <source>Trans. China Electrotechnical Soc.</source> <volume>34</volume> (<issue>23</issue>), <fpage>4982</fpage>&#x2013;<lpage>4990</lpage>. </citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>Zhong.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Zhenyong.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Weihua.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Analysis and Comparison of Primary Frequency Control Technology for Wind Power and thermal Power Unit</article-title>. <source>Electric Power Automation Equipment</source> <volume>37</volume> (<issue>12</issue>), <fpage>92</fpage>&#x2013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.16081/j.issn.1006-6047.2017.12.013</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>