<?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">1074916</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2022.1074916</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>An optimal energy storage system sizing determination for improving the utilization and forecasting accuracy of photovoltaic (PV) power stations</article-title>
<alt-title alt-title-type="left-running-head">Li et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2022.1074916">10.3389/fenrg.2022.1074916</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Bin</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2058058/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Mingzhe</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yan</surname>
<given-names>Shiye</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Yifan</given-names>
</name>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shi</surname>
<given-names>Bowen</given-names>
</name>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Ye</surname>
<given-names>Jilei</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1925789/overview"/>
</contrib>
</contrib-group>
<aff>
<institution>School of Energy Science and Engineering</institution>, <institution>Nanjing Tech University</institution>, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1259467/overview">Liansong Xiong</ext-link>, Xi&#x2019;an Jiaotong University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2063666/overview">Yang Liu</ext-link>, ENGIE, France</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1463159/overview">Huimin Wang</ext-link>, University of Electronic Science and Technology of China, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2067707/overview">Kai Zhou</ext-link>, Electric Power Group LLC, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2063651/overview">Tao Shi</ext-link>, Nanjing University of Posts and Telecommunications, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jilei Ye, <email>yejilei@njtech.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>1074916</elocation-id>
<history>
<date date-type="received">
<day>20</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>02</day>
<month>11</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Li, Li, Yan, Zhang, Shi and Ye.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Li, Li, Yan, Zhang, Shi and Ye</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>As a new type of flexible regulation resource, energy storage systems not only smooth out the fluctuation of new energy generation but also track the generation scheduling combined with new energy power to enhance the reliability of new energy system operations. In recent years, installing energy storage for new on-grid energy power stations has become a basic requirement in China, but there is still a lack of relevant assessment strategies and techno-economic evaluation of the size determination of energy storage systems from the perspective of new energy power stations. Therefore, this paper starts from summarizing the role and configuration method of energy storage in new energy power stations and then proposes multidimensional evaluation indicators, including the solar curtailment rate, forecasting accuracy, and economics, which are taken as the optimization targets for configuring energy storage systems in PV power stations. Lastly, taking the operational data of a 4000&#xa0;MWPV plant in Belgium, for example, we develop six scenarios with different ratios of energy storage capacity and further explore the impact of energy storage size on the solar curtailment rate, PV curtailment power, and economics. The method proposed in this paper is effective for the performance evaluation of large PV power stations with annual operating data, realizes the automatic analysis on the optimal size determination of energy storage system for PV power stations, and verifies the rationality of the principle for configuring energy storage for PV power stations in some regions of China.</p>
</abstract>
<kwd-group>
<kwd>PV power station</kwd>
<kwd>solar curtailment rate</kwd>
<kwd>forecasting accuracy</kwd>
<kwd>economic analysis</kwd>
<kwd>energy storage system sizing</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Highlights</title>
<p>
<list list-type="simple">
<list-item>
<p>1) This paper starts by summarizing the role and configuration method of energy storage in new energy power station and then proposes a new evaluation index system, including the solar curtailment rate, forecasting accuracy, and economics, which are taken as the optimization targets for configuring energy storage system in PV power stations.</p>
</list-item>
<list-item>
<p>2) It calculates the forecasting accuracy of PV power and the penalty cost for each month in terms of penalty rules in a certain region, provides the specific algorithm process of energy storage sizing, and derives the simulation results.</p>
</list-item>
<list-item>
<p>3) The method proposed in this paper is applicable to the analysis of the operating characteristics of large PV power stations within the whole year and realizes the automatic analysis of the optimal scheme for the configuration with energy storage on PV power stations.</p>
</list-item>
</list>
</p>
</sec>
<sec sec-type="intro" id="s2">
<title>1 Introduction</title>
<p>In recent years, solar energy has been widely regarded as one of the largest green energy and plays a major role in modern power system (<xref ref-type="bibr" rid="B25">Ndebele, 2020</xref>). Combined with the &#x201c;carbon peak and carbon neutral&#x201d; target proposed by China (<xref ref-type="bibr" rid="B28">Qian et al., 2022</xref>), it is expected that by 2030, the installed capacity of domestic PV power generation will exceed 1000&#xa0;GW (<xref ref-type="bibr" rid="B15">Kut et al., 2021</xref>., Meanwhile, by the end of 2021, the cumulative installed solar power capacity in Germany reached 66.5 GW and is expected to exceed 100&#xa0;GW by 2030. In addition, Australia&#x2019;s installed capacity has also reached 25.3&#xa0;GW (<xref ref-type="bibr" rid="B1">Ali et al., 2019</xref>; <xref ref-type="bibr" rid="B3">Apeh et al., 2022</xref>).</p>
<p>However, there remain two major issues when integrating it into the power grid due to the consequent intermittent and fluctuating features (<xref ref-type="bibr" rid="B6">Chao, 2011</xref>; <xref ref-type="bibr" rid="B14">Kim et al., 2020</xref>). In addition, the grid dispatching center uses the forecasting power of the PV power station as reference for scheduling plan, but there is still a large deviation between forecasting power and actual generation (<xref ref-type="bibr" rid="B13">Keeratimahat et al., 2021</xref>), which may not only cause the occurrence of PV curtailment (<xref ref-type="bibr" rid="B5">Bird et al., 2016</xref>; <xref ref-type="bibr" rid="B33">Zhang et al., 2022a</xref>) but also result in low system reliability and security due to inadequate backup resources. Furthermore, it is possible to face high penalty for PV power station (<xref ref-type="bibr" rid="B4">Azimov and Avezova, 2022</xref>).</p>
<p>Energy storage is one of the most effective solutions to smooth out new energy power fluctuations (<xref ref-type="bibr" rid="B7">Chen et al., 2021</xref>; <xref ref-type="bibr" rid="B31">Yang et al., 2022</xref>), promote high penetration of grid-connected green energy, and reduce the forecasting deviations of PV power, because of its flexible dual characteristics of charging and discharging (<xref ref-type="bibr" rid="B23">Lutsenko and Fetsov, 2020</xref>). Therefore, the development of energy storage system will promote energy conservation and emission reduction (<xref ref-type="bibr" rid="B20">Lu et al., 2022</xref>). However, the initial investment cost of energy storage is high (<xref ref-type="bibr" rid="B9">Du et al., 2022</xref>), which makes it necessary to figure out a reasonable configuring capacity. J. <xref ref-type="bibr" rid="B17">Li et al. (2020</xref>) propose a capacity optimization method for combined PV and storage systems, which considers the power allocation for PV and storage systems with the objective of economic optimality; P. D. <xref ref-type="bibr" rid="B21">Lund (2018</xref>) considered the PV self-consumption, as well as the sensitivity of the storage system size of weather, and finally obtained the economically optimal PV energy storage system configuration; R. <xref ref-type="bibr" rid="B22">Luthander et al. (2019</xref>) consider the matching of PV generation with demand and the transfer of solar energy through energy storage systems, which can achieve zero-energy buildings that meet the relevant requirements of EU regulations; and R. <xref ref-type="bibr" rid="B8">Dom&#xed;nguez et al. (2020</xref>) analyzed the operation and construction of the European power system in 2050 and proposed a stochastic model to study the future renewable energy construction planning. The results showed that new renewable energy generation capacity in Europe should reach 881&#xa0;GW by 2050, taking into account the installation of energy storage system, and until then, CO<sub>2</sub> emissions will be reduced by 77%. F. <xref ref-type="bibr" rid="B32">Zhang et al. (2020a</xref>) proposed a predictive control strategy of energy storage system to improve AGC regulation capacity for high PV penetration system; meanwhile, a protective charging and discharging method is developed to extend energy storage system life. O. <xref ref-type="bibr" rid="B2">Alrawi et al. (2022</xref>) utilized empirical evidence and an economic model to evaluate rooftop PV systems in Qatar with different economic tools and collected real data, and S. <xref ref-type="bibr" rid="B29">Ud-Din Khan et al. (2022</xref>) developed theoretical models to evaluate the economy of integrated PV systems. Various types of lithium-ion and liquid-flow batteries are then evaluated technically and economically to determine appropriate type of storage. E. <xref ref-type="bibr" rid="B26">Oh and Son (2020</xref>), based on the uncertainty of wind power forecast, proposed a method of determining the theoretical energy storage capacity considering the stochastic characteristics of uncertainty. Finally, the effectiveness of the model is verified by the actual wind power generation data and the forecast data; V. Jani and H. Abdi (<xref ref-type="bibr" rid="B12">Jani and Abdi, 2018</xref>) established three different and incompatible objective functions considering operation cost, voltage deviation, and air emission to determine the optimal configuration of energy storage capacity. Two multiobjective hybrid algorithms are used to solve the problem. Simulation results show the effectiveness of the proposed method; S. Garip and S. Ozdemir (<xref ref-type="bibr" rid="B11">Garip and Ozdemir, 2022</xref>) chose energy cost minimization as the objective function. The proposed algorithm can obtain the optimal PV and energy storage size with minimum cost using the new energy management method and the PSO algorithm.</p>
<p>In summary, there have been many studies on energy storage sizing in PV power systems, but there are few sizing models with consideration of assessment indicators in terms of local new energy policy and economic factors. From the side of new energy generation, installing energy storage systems not only can improve the operating characteristics of PV power station but can also indirectly improve the system reliability and environmental protection. Therefore, this paper builds a new energy storage sizing model for large PV power stations, with the aim to improve the solar curtailment rate and forecasting accuracy. First, with processing original data of PV power generation, the annual curtailment and forecasting accuracy of PV power are quantified. Then, we calculate the forecasting accuracy of PV power and the penalty cost for each month in terms of penalty rules in a certain region, provide the specific algorithm process of energy storage sizing, and derive the simulation results. A comprehensive energy storage system size determination strategy is obtained with the trade-off among the solar curtailment rate, the forecasting accuracy, and financial factors, which provides a practical reference to determine energy storage size for PV power station and further verifies the feasibility of energy storage system in the high PV penetration power system.</p>
</sec>
<sec id="s3">
<title>2 Role of energy storage in PV power stations and deployment rules in China</title>
<sec id="s3-1">
<title>2.1 Roles of energy storage systems in PV power stations</title>
<p>Chinese renewable energy enters a new stage of high-quality leap; in the first half of 2022, nonfossil energy power generation accounted for 83% of new power generation installed capacity, while renewable energy generation exceeded 1.1 billion kW (<xref ref-type="bibr" rid="B27">Pablo-Romero et al., 2021</xref>). The generating capacity of hydropower, wind power, and solar power has increased by 20.3%, 7.8%, and 13.5% from the last year, respectively. Large-scale wind power and PV on-grid has posed a huge challenge to the safety and reliability of power system (<xref ref-type="bibr" rid="B16">Lei et al., 2022</xref>; <xref ref-type="bibr" rid="B19">Liu et al., 2022</xref>). From the perspective of system operation, installing energy storage not only mitigates PV power fluctuation but also realizes real-time scheduling of tracking power and enhance the grid&#x2019;s ability to consume new energy (<xref ref-type="bibr" rid="B30">Wang et al., 2019</xref>); in addition, the joint operation of PV power generation and energy storage can play the role of peak shaving and valley filling and also participate in auxiliary services (<xref ref-type="bibr" rid="B24">Motalleb et al., 2016</xref>). On the other side, from the perspective of new energy power plant, installing reasonably energy storage can improve the accuracy of PV forecasting (<xref ref-type="bibr" rid="B18">Liu et al., 2020</xref>), which may not only reduce the penalty cost but also further reduce the solar curtailment rate and improve net revenue of PV power station. In general, most optimization objectives of combing with PV power generation and energy storage focus on system operation, including meeting the technical requirements of grid power quality, minimizing power fluctuations, providing stable output power, and improving the utilization rate of new energy. With different optimization objectives, sizing procedures of energy storage system are developed, but the size determination strategy of energy storage systems from the perspective of renewable power plant is not enough. The sizing of energy storage systems in PV power plants is closely related to the operation mode, market rules, and financial factors. Installing energy storage system with reasonable capacity is necessary for power plant operation; therefore, an optimal sizing strategy of energy storage system in PV power plants is very important and meaningful.</p>
</sec>
<sec id="s3-2">
<title>2.2 Deployment rules of energy storage in PV power stations in China</title>
<p>So far in 2021, the deployment rules of energy storage for new energy plant have been put forward in 24 provinces of China, of which governments have made clear requirements for energy storage supporting distributed PV. In all configuring rules of energy storage, the highest proportion of energy storage capacity requirements in Shandong Zaozhuang is 15%&#x2013;30% of the installed PV rated capacity, and the duration time can be 2&#x2013;4&#xa0;h, while in some regions of Henan and Shanxi province the capacity of energy storage reaches 20%. Mostly the duration time is 2&#xa0;h, and 3&#xa0;h is required for the market-based on-grid project in Hebei province. We investigate the relevant rules in Jiangsu province in terms of energy storage capacity requirements, the solar curtailment rate, and the minimum standards for forecasting assessment. Specifically, Jiangsu Development and Reform Commission issued the &#x201c;Provincial Development and Reform Commission on the province&#x2019;s 2021 photovoltaic power generation project market notice&#x201d;, in which the energy storage system with 8% of PV rated power generation and 2&#xa0;h duration is needed principally for projects located in south of the Yangtze River; for those projects located in north of the Yangtze River, 10% of PV rated power and 2&#xa0;h duration are required. In terms of the solar curtailment rate, according to the &#x201c;Clean Energy Consumption Action Plan (2018&#x2013;2020)&#x201d; developed by the National Development and Reform Commission and the National Energy Administration, the utilization rate of PV power generation must be higher than 95% and the solar curtailment rate must be less than 5% (<xref ref-type="bibr" rid="B35">Zhang et al., 2022b</xref>). For assessment indicators, according to the specific requirements of the &#x201c;Jiangsu Province Electricity Grid Operation Management Rules&#x201d;, the time interval is set at 15&#xa0;min, and the PV forecasting qualification rate is also defined in the rules. The criteria for qualification assessment are calculated monthly, and the penalty expense for every month is also standardized.</p>
<p>Based on aforementioned regulations and standards, this paper will figure out the economics of PV power plants with different sizes of energy storage systems. In order to calculate the relevant evaluation indicators from the massive dataset, this paper uses statistics to carry out the economic analysis of energy storage configuration scheme.</p>
</sec>
</sec>
<sec id="s4">
<title>3 Statistical analysis and extraction of feature indicators with PV power station data processing</title>
<sec id="s4-1">
<title>3.1 Solar curtailment rate</title>
<p>Through the statistics and analysis on the sample dataset of PV plant, the PV curtailment power at each sample point is measured, as shown in <xref ref-type="disp-formula" rid="e1">Eq. 1</xref>.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>max</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">I</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the PV curtailment power at the sampling poinI i, and <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the actual and forecasting power at the sampling point i, respectively. When <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, all of the PV power can be dispatched; when <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, there is a surplus of PV power, and the phenomenon of solar curtailment occurs without energy storage system (<xref ref-type="bibr" rid="B34">Zhang et al., 2020b</xref>). By calculating the monthly curtailment power and solar curtailment rate (<xref ref-type="bibr" rid="B10">Fang et al., 2017</xref>), actual on-grid power of PV station in each month can be indirectly deduced. The total monthly solar curtailment power and solar curtailment rate is shown in <xref ref-type="disp-formula" rid="e2">Eqs 2</xref>, <xref ref-type="disp-formula" rid="e3">3</xref>.<disp-formula id="e2">
<mml:math id="m7">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the PV curtailment power dataset for each month, and <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> denotes the number of days at the month&#xa0;m.<disp-formula id="e3">
<mml:math id="m10">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>where &#x3bb;(m) denotes the dataset of solar curtailment rate in each month.</p>
</sec>
<sec id="s4-2">
<title>3.2 Forecasting qualification rate of PV power generation</title>
<sec id="s4-2-1">
<title>3.2.1 Monthly forecasting qualification rate</title>
<p>According to &#x201c;Jiangsu Province Electricity Grid Operation Management Rules,&#x201d; the short- and medium-term PV forecasting power are computed by 96 points, and we obtain the qualification rate from <xref ref-type="disp-formula" rid="e4">Eq. 4</xref>. The next day qualification rate is required not less than 90%. The qualification rate of PV forecasting power is determined by the ratio of qualified points to total points for every month. When the unqualified rate exceeds 2%, PV power station will be punished by &#xa5;10 per 10,000&#xa0;kW of rated capacity.<disp-formula id="e4">
<mml:math id="m11">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mo>%</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>100</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>In <xref ref-type="disp-formula" rid="e4">Eq. 4</xref>, <inline-formula id="inf8">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the actual power at the sample point i, <inline-formula id="inf9">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the forecasted power at the sample point i, and <inline-formula id="inf10">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is PV rated capacity.</p>
</sec>
<sec id="s4-2-2">
<title>3.2.2 Statistical analysis process</title>
<p>In order to statistically calculate the number of unqualified PV forecasting points from a huge dataset and quantified the penalty cost for each month, this paper proposes a method to count the monthly unqualified PV forecasting points and specify the penalty cost; the process is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Statistical process of unqualified forecasting points for PV power station. The steps are as follows.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g001.tif"/>
</fig>
<p>
<statement content-type="step" id="Step_1">
<label>Step 1</label>
<p>Input data, calculate the forecasting power deviation and qualification rate for each point using <xref ref-type="disp-formula" rid="e4">Eq. 4</xref>.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_2">
<label>Step 2</label>
<p>The number of sampling points for each month is calculated to determine the specific data range from each month&#x2019;s dataset, as shown in <xref ref-type="disp-formula" rid="e5">Eq. 5</xref>.<disp-formula id="e5">
<mml:math id="m15">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>where S(m) represents the specific data range corresponding to the data of month&#xa0;m in the annual dataset, and d(m) represents the number of days in month&#xa0;m. When and only when m &#x3d; 1, <inline-formula id="inf11">
<mml:math id="m16">
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>&#x3002;</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_3">
<label>Step 3</label>
<p>Count the number of unqualified PV forecasting points in each month as shown in <xref ref-type="disp-formula" rid="e6">equation 6</xref>, and calculate the percentage of unqualified points in each month as shown in <xref ref-type="disp-formula" rid="e7">Eq. 7</xref>.<disp-formula id="e6">
<mml:math id="m17">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>90</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m18">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>where N(m) is the number of unqualified points in month&#xa0;m, and &#x3b1;[S(m)] indicates the qualification rate in month&#xa0;m. In <xref ref-type="disp-formula" rid="e7">Eq. 7</xref>, &#x3c6;(m) is the percentage of monthly unqualified points.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_4">
<label>Step 4</label>
<p>Determine whether the percentage of unqualified points is greater than 2% per month. If it is greater than 2%, calculate the monthly penalty, as shown in <xref ref-type="disp-formula" rid="e8">Eq. 8</xref>, and if less than or equal to 2%, the penalty is zero in the month.</p>
<p>Let <inline-formula id="inf12">
<mml:math id="m19">
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, then<disp-formula id="e8">
<mml:math id="m20">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>d</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:munderover>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>90</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x222a;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x3c6;</mml:mi>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>%</mml:mo>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <inline-formula id="inf13">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the absolute value of the PV forecasting deviation at the sample point i, <inline-formula id="inf14">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> denotes the penalty in month&#xa0;m, and <inline-formula id="inf15">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents penalty cost of per unit unqualified forecasting PV power.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_5">
<label>Step 5</label>
<p>Determine whether all of the monthly data have been analyzed, otherwise return to step 2.</p>
</statement>
</p>
<p>
<statement content-type="step" id="Step_6">
<label>Step 6</label>
<p>Output: calculation results including penalty and the number of unqualified points for each month.</p>
</statement>
</p>
</sec>
</sec>
</sec>
<sec id="s5">
<title>4 Sizing determination strategy of energy storage considering assessment indicators and financial factor</title>
<p>The sizing process of energy storage capacity is programmed with MATLAB, as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. By extracting and processing the raw dataset after performing multiple cycles, the feature indicators mentioned in <xref ref-type="sec" rid="s2">section 2</xref> are developed.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Implementation process of energy storage capacity size determination.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g002.tif"/>
</fig>
<p>The steps are as follows:</p>
<p>
<statement>
<p>Step 1Input basic dataset and initial parameters including energy storage capacity <inline-formula id="inf16">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and rated power <inline-formula id="inf17">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> for different scenarios, cost per unit capacity <inline-formula id="inf18">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, cost per unit power <inline-formula id="inf19">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and PV on-grid tariff <inline-formula id="inf20">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mi mathvariant="normal">v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</statement>
</p>
<p>
<statement>
<p>Step 2Input power constraints, that is, when the first cycle is performed, the PV curtailment power of each sample point is calculated by <xref ref-type="disp-formula" rid="e1">Eq. 1</xref> to form curtailment dataset that corresponds to each time point, and the curtailment power dataset is formed by <xref ref-type="disp-formula" rid="e9">Eq. 9</xref> after filtering the adsorption by energy storage. If the curtailment power at a certain moment is greater than the rated power of energy storage, the power cannot be absorbed and set as zero; meanwhile, curtailment dataset is updated.<disp-formula id="e9">
<mml:math id="m29">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>&#x222a;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mn>1,35040</mml:mn>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where <inline-formula id="inf21">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the PV curtailment dataset composed by consumption power of energy storage.<disp-formula id="e10">
<mml:math id="m31">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.25</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>95</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:mrow>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">q</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>95</mml:mn>
<mml:mo>&#x223c;</mml:mo>
<mml:mn>96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mn>1,365</mml:mn>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">j</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mn>96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>95,96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where <inline-formula id="inf22">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mo>,</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the total amount of PV curtailment power that has been consumed by the energy storage until sampling point j in day&#xa0;T.</p>
<p>Step 3Energy constraints: the energy storage system is set one cycle per day; therefore, the dataset is converted from a scale of 96&#x2a;365 to a scale of 1&#x2a;365. Then the second cycle is carried out on the basis of the end of the first cycle, which takes into account whether the daily energy storage capacity reaches 90% of the rated capacity, and also tracks the real-time curtailment power being consumed at every 15&#xa0;min sample point. The energy is calculated per day as shown in Eq. 10.</p>
<p>In order to extend the life of battery energy storage, the charging and discharging capacity limit is stipulated as 90% of the rated capacity. Therefore, when the storage charging energy is greater than 90% of the rated capacity, the storage capacity needs to be defaulted on the day.<disp-formula id="e11">
<mml:math id="m33">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">T</mml:mi>
<mml:mo>,</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.9</mml:mn>
<mml:mo>&#x2a;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m34">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">q</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="normal">j</mml:mi>
<mml:mo>&#x223c;</mml:mo>
<mml:mn>96</mml:mn>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>At this point, the energy storage battery has reached saturation after being completely charged during the day, and no longer absorbs the PV curtailment power; therefore, the real-time dataset of the solar curtailment power that can be absorbed by the battery is updated again.</p>
</statement>
</p>
<p>
<statement>
<p>Step 4When all data are cycled through for 365 days in a year, the net revenue of energy storage C is obtained as shown in <xref ref-type="disp-formula" rid="e13">Eq. 13</xref>. The results of the solar curtailment rate, C, and the loss cost are outputted.<disp-formula id="e13">
<mml:math id="m35">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>365</mml:mn>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
</p>
</statement>
</p>
</sec>
<sec id="s6">
<title>5 Case studies</title>
<sec id="s6-1">
<title>5.1 Basic data</title>
<p>In this paper, the original data of 35,040 sample points with an interval of 15&#xa0;min for a 4000&#xa0;MW&#xa0;PV plant in Belgium are used to study the size determination of energy storage. From the rules mentioned in <xref ref-type="sec" rid="s2">section 2.2</xref>, the energy storage capacity ratio is set as shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Set value of energy storage power and capacity.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">The rated power ratio of energy storage with PV (%)</th>
<th align="left">Power/MW</th>
<th align="left">Duration/h</th>
<th align="left">The energy storage system capacity/MWh</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">8</td>
<td align="char" char=".">320</td>
<td align="char" char=".">2</td>
<td align="char" char=".">640</td>
</tr>
<tr>
<td align="left">10</td>
<td align="char" char=".">400</td>
<td align="char" char=".">2</td>
<td align="char" char=".">800</td>
</tr>
<tr>
<td align="left">11</td>
<td align="char" char=".">440</td>
<td align="char" char=".">2</td>
<td align="char" char=".">880</td>
</tr>
<tr>
<td align="left">12</td>
<td align="char" char=".">480</td>
<td align="char" char=".">2</td>
<td align="char" char=".">960</td>
</tr>
<tr>
<td align="left">13</td>
<td align="char" char=".">520</td>
<td align="char" char=".">2</td>
<td align="char" char=".">1,040</td>
</tr>
<tr>
<td align="left">15</td>
<td align="char" char=".">600</td>
<td align="char" char=".">2</td>
<td align="char" char=".">1,200</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>From <xref ref-type="table" rid="T1">Table 1</xref>, six scenarios with energy storage power ratios of 8%, 10%, 11%, 12%, 13%, and 15% are built. To discuss the techno-economic indicators proposed in <xref ref-type="sec" rid="s2">section 2</xref>, the initial parameters are described in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Main parameter settings.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Parameter</th>
<th align="left">Setting</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">The forecasting data and raw data of PV power</td>
<td align="left">96&#x2a;365,96&#x2a;365</td>
</tr>
<tr>
<td align="left">Rated power <inline-formula id="inf23">
<mml:math id="m36">
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>,Capacity <inline-formula id="inf24">
<mml:math id="m37">
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">Configured at 8%&#x2013;15% of installed PV capacity</td>
</tr>
<tr>
<td align="left">Cost per unit capacity <inline-formula id="inf25">
<mml:math id="m38">
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">1.3&#xa5;/Wh</td>
</tr>
<tr>
<td align="left">Cost per unit power <inline-formula id="inf26">
<mml:math id="m39">
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">S</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">0.08&#xa5;/W</td>
</tr>
<tr>
<td align="left">PV on-grid tariff <inline-formula id="inf27">
<mml:math id="m40">
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mi mathvariant="normal">v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">0.5&#xa5;/kWh</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s6-2">
<title>5.2 Analysis of feature indicators</title>
<p>The PV curtailment power of each sample point can be obtained from dataset by <xref ref-type="disp-formula" rid="e1">Eq. 1</xref>, and then the actual power and curtailment power of PV plant for the whole year is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Actual power and curtailment power of PV plant in a whole year.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g003.tif"/>
</fig>
<p>From <xref ref-type="fig" rid="F3">Figure 3</xref>, it can be seen that the solar curtailment situation of the PV plant is quite serious. To further quantify the solar curtailment, the monthly solar curtailment power and solar curtailment rate are performed by <xref ref-type="disp-formula" rid="e2">Eqs 2</xref>, <xref ref-type="disp-formula" rid="e3">3</xref>, respectively. Then, actual PV on-grid power for each month can be figured out, as shown in <xref ref-type="fig" rid="F4">Figure 4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Monthly solar curtailment power and solar curtailment rate.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g004.tif"/>
</fig>
<p>In <xref ref-type="fig" rid="F4">Figure 4</xref>, it is clear that the solar curtailment rate in January is as high as 19.5%, the lowest in June is only 5.5%, and the average solar curtailment rate is as high as 8.07%. The solar curtailment occurred with low forecasting accuracy at some time points, which can be alleviated by installing storage devices for PV power plant. To further investigate the impact of PV forecasting accuracy on solar curtailment rate, the actual power generated at each point is compared with the forecasting power of the previous day, as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Actual power vs forecasting power of PV power station.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g005.tif"/>
</fig>
<p>As can be seen in <xref ref-type="fig" rid="F5">Figure 5</xref>, there are still many sample points with significant forecasting deviation. When dispatching power is too high, PV power cannot meet the load demand, and when dispatching power is too low, PV power generation is surplus. According to the regulation of &#x201c;Jiangsu Province Electricity Grid Operation Management Rules,&#x201d; the unqualified points of PV forecasting power, the forecasting deviation, and the qualification rate of each sample point can be provided. The results are shown in <xref ref-type="fig" rid="F6">Figure 6</xref>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Forecasting deviation and qualification rate of PV power station.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g006.tif"/>
</fig>
<p>From <xref ref-type="fig" rid="F6">Figure 6</xref>, it can be seen that the deviation between the forecasting and the actual power of some points is large, and there also exist many unqualified forecasting points.</p>
</sec>
<sec id="s6-3">
<title>5.3 Results and discussion</title>
<sec id="s6-3-1">
<title>5.3.1 Simulation results</title>
<p>To get the total number of unqualified points for each month from the huge dataset and thus calculating penalty cost, the simulation results are outputted by completing steps one to six in <xref ref-type="sec" rid="s3">section 3</xref>, as shown in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Statistics and simulation results.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Month no.</th>
<th align="left">Day no.</th>
<th align="left">Point no.</th>
<th align="left">Monthly unqualified points number</th>
<th align="left">Percentage of monthly unqualified point (%)</th>
<th align="left">Whether &#x2265;2%</th>
<th align="left">The sum of unqualified predicted power(MW)</th>
<th align="left">Penalty cost (&#xa5;10,000 unit)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">42</td>
<td align="char" char=".">1.4</td>
<td align="left">N</td>
<td align="left">&#x2014;</td>
<td align="char" char=".">0</td>
</tr>
<tr>
<td align="left">2</td>
<td align="char" char=".">28</td>
<td align="char" char=".">2,688</td>
<td align="char" char=".">166</td>
<td align="char" char=".">6.18</td>
<td align="left">Y</td>
<td align="left">95,280</td>
<td align="char" char=".">9.53</td>
</tr>
<tr>
<td align="left">3</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">219</td>
<td align="char" char=".">7.34</td>
<td align="left">Y</td>
<td align="left">130,702</td>
<td align="char" char=".">13.07</td>
</tr>
<tr>
<td align="left">4</td>
<td align="char" char=".">30</td>
<td align="char" char=".">2,880</td>
<td align="char" char=".">348</td>
<td align="char" char=".">12.08</td>
<td align="left">Y</td>
<td align="left">200,558</td>
<td align="char" char=".">20.06</td>
</tr>
<tr>
<td align="left">5</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">226</td>
<td align="char" char=".">7.6</td>
<td align="left">Y</td>
<td align="left">127,955</td>
<td align="char" char=".">12.8</td>
</tr>
<tr>
<td align="left">6</td>
<td align="char" char=".">30</td>
<td align="char" char=".">2,880</td>
<td align="char" char=".">208</td>
<td align="char" char=".">7.2</td>
<td align="left">Y</td>
<td align="left">121,992</td>
<td align="char" char=".">12.2</td>
</tr>
<tr>
<td align="left">7</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">216</td>
<td align="char" char=".">7.3</td>
<td align="left">Y</td>
<td align="left">127,626</td>
<td align="char" char=".">12.76</td>
</tr>
<tr>
<td align="left">8</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">178</td>
<td align="char" char=".">6.0</td>
<td align="left">Y</td>
<td align="left">91,644</td>
<td align="char" char=".">9.16</td>
</tr>
<tr>
<td align="left">9</td>
<td align="char" char=".">30</td>
<td align="char" char=".">2,880</td>
<td align="char" char=".">181</td>
<td align="char" char=".">6.28</td>
<td align="left">Y</td>
<td align="left">110,362</td>
<td align="char" char=".">11.04</td>
</tr>
<tr>
<td align="left">10</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">48</td>
<td align="char" char=".">1.6</td>
<td align="left">N</td>
<td align="left">&#x2014;</td>
<td align="char" char=".">0</td>
</tr>
<tr>
<td align="left">11</td>
<td align="char" char=".">30</td>
<td align="char" char=".">2,880</td>
<td align="char" char=".">29</td>
<td align="char" char=".">1</td>
<td align="left">N</td>
<td align="left">&#x2014;</td>
<td align="char" char=".">0</td>
</tr>
<tr>
<td align="left">12</td>
<td align="char" char=".">31</td>
<td align="char" char=".">2,976</td>
<td align="char" char=".">31</td>
<td align="char" char=".">1.04</td>
<td align="left">N</td>
<td align="left">&#x2014;</td>
<td align="char" char=".">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>From the results in <xref ref-type="table" rid="T3">Table 3</xref>, months of the forecasting power assessment in 2021 is unqualified, and the total penalty to be paid is about &#xa5;1,006,200, which seriously affects the efficiency of system dispatching and causes additional expenses for the PV power plant. To solve this problem, energy storage system should be installed for PV power stations. The simulation results for six scenarios with different energy storage sizes are shown in <xref ref-type="table" rid="T4">Table 4</xref>.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Simulation results with different energy storage sizes.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="left">The power ratio of energy storage with PV (%)</th>
<th rowspan="2" align="left">Power/MW</th>
<th rowspan="2" align="left">Capacity/MWh</th>
<th rowspan="2" align="left">Annual investment cost/&#xa5;10,000 unit</th>
<th rowspan="2" align="left">Annual solar curtailment reduction/MWh</th>
<th rowspan="2" align="left">Net income/&#xa5;10,000 unit</th>
<th colspan="2" align="left">Solar curtailment rate</th>
<th colspan="2" align="left">Annual unqualified point</th>
<th colspan="2" align="left">Penalty/&#xa5;10,000 unit</th>
</tr>
<tr>
<th align="left">Without energy storage</th>
<th align="left">With energy storage</th>
<th align="left">Without energy storage</th>
<th align="left">With energy storage</th>
<th align="left">Without energy storage</th>
<th align="left">With energy storage</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">8</td>
<td align="char" char=".">320</td>
<td align="char" char=".">640</td>
<td align="char" char=".">8,576</td>
<td align="char" char=".">175,503</td>
<td align="char" char=".">23.67</td>
<td align="char" char=".">8.07%</td>
<td align="char" char=".">4.32%</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">100.62</td>
<td align="char" char=".">100.62</td>
</tr>
<tr>
<td align="left">10</td>
<td align="char" char=".">400</td>
<td align="char" char=".">800</td>
<td align="char" char=".">10,720</td>
<td align="char" char=".">220,935</td>
<td align="char" char=".">105.81</td>
<td align="char" char=".">8.07%</td>
<td align="char" char=".">3.34%</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">100.62</td>
<td align="char" char=".">100.62</td>
</tr>
<tr>
<td align="left">11</td>
<td align="char" char=".">440</td>
<td align="char" char=".">880</td>
<td align="char" char=".">11,792</td>
<td align="char" char=".">242,183</td>
<td align="char" char=".">317.17</td>
<td align="char" char=".">8.07%</td>
<td align="char" char=".">2.89%</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">1,689</td>
<td align="char" char=".">100.62</td>
<td align="char" char=".">92.52</td>
</tr>
<tr>
<td align="left">12</td>
<td align="char" char=".">480</td>
<td align="char" char=".">960</td>
<td align="char" char=".">12,864</td>
<td align="char" char=".">260,905</td>
<td align="char" char=".">181.25</td>
<td align="char" char=".">8.07%</td>
<td align="char" char=".">2.5%</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">1,526</td>
<td align="char" char=".">100.62</td>
<td align="char" char=".">85.55</td>
</tr>
<tr>
<td align="left">13</td>
<td align="char" char=".">520</td>
<td align="char" char=".">1,040</td>
<td align="char" char=".">13,936</td>
<td align="char" char=".">276,752</td>
<td align="char" char=".">&#x2212;98.37</td>
<td align="char" char=".">8.07%</td>
<td align="char" char=".">2.15%</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">1,399</td>
<td align="char" char=".">100.62</td>
<td align="char" char=".">79.76</td>
</tr>
<tr>
<td align="left">15</td>
<td align="char" char=".">600</td>
<td align="char" char=".">1,200</td>
<td align="char" char=".">16,080</td>
<td align="char" char=".">308,086</td>
<td align="char" char=".">&#x2212;675.7</td>
<td align="char" char=".">8.07%</td>
<td align="char" char=".">1.5%</td>
<td align="char" char=".">1892</td>
<td align="char" char=".">1,175</td>
<td align="char" char=".">100.62</td>
<td align="char" char=".">68.06</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>To describe the effect of energy storage size on each assessment indicator, more discussions including the influence on the net income, the solar curtailment rate, and the forecasting qualification rate are developed.</p>
</sec>
<sec id="s6-3-2">
<title>5.3.2 Impact of energy storage capacity on net income</title>
<p>From the perspective of the net income of energy storage, with the increase of the storage capacity, the net income of energy storage shows a trend of first increase and then decrease, as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. As we know, energy storage is charged when the PV power is in surplus and discharged at night, an increasing income present with the increase of power ratio, and the maximum annual net income reaches 3.17 million. However, when the power ratio exceeds 11%, the income begins to decline, and when the ratio exceeds by 12%, the income is even negative. With the increase of the storage capacity, the investment cost of energy storage rises and the income obtained from the PV consumption cannot offset the annual investment cost of energy storage, as shown in <xref ref-type="fig" rid="F8">Figure 8</xref>.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Net income of energy storage with different rated power ratios.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g007.tif"/>
</fig>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Net income and annual investment cost of energy storage.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g008.tif"/>
</fig>
</sec>
<sec id="s6-3-3">
<title>5.3.3 Impact of energy storage capacity on the solar curtailment rate</title>
<p>When the power ratio is 8%, the annual solar curtailment rate is already lower than 5%. With the increase of the storage capacity, the solar curtailment rate decreased continuously from 4.32% to 1.5%, and the PV utilization rate is greatly improved, as shown in <xref ref-type="fig" rid="F9">Figure 9</xref>. However, with the large size of energy storage, the net income of energy storage is negative when the power ratio exceeds 12% even the improvement of the PV consumption rate. Therefore, the size determination of energy storage should involve the solar curtailment rate and economics.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Solar curtailment rate and net income with different power ratios.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g009.tif"/>
</fig>
<p>When the capacity ratio is 11% (440MW/880&#xa0;MWh), the maximum net revenue of energy storage is obtained. Here, the solar curtailment rate is 2.89%, which is far below the specified value. The comparison of the solar curtailment power at each point without and with energy storage is shown in <xref ref-type="fig" rid="F10">Figure 10</xref>. It is obvious for reducing the curtailment for PV power stations.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Solar curtailment power without and with energy storage.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g010.tif"/>
</fig>
</sec>
<sec id="s6-3-4">
<title>5.3.4 Effect of energy storage capacity on the forecasting qualification rate</title>
<p>As the rated power ratio increases, the penalty and annual unqualified forecasting points decreases, as shown in <xref ref-type="fig" rid="F11">Figure 11</xref>. The penalty remains the same level from 8% to 10%, and the number of annual unqualified forecasting points does not change, which means that the energy storage size is not enough.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Annual unqualified forecasting power points and penalty with different rated power ratios.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g011.tif"/>
</fig>
<p>The net income of energy storage does not increase continuously while improving the forecasting qualification rate and reducing the penalty, as seen in <xref ref-type="fig" rid="F7">Figure 7</xref>. The net income is negative with the ratio of 13%. The highest income appears with the ratio of 11%; meanwhile, the penalty cost decreases by &#xa5;81,000 and the unqualified forecasting point number drops by 203. When the ratio is 12%, the penalty cost and unqualified points reduction are &#xa5;150,700and 330, respectively. However, the net income is 1.36 million lower than that of 11%. Therefore, the reasonable capacity will be confirmed as 440MW/880&#xa0;MWh. The deviation of forecasting power and actual power is clearly improved, as shown in <xref ref-type="fig" rid="F12">Figure 12</xref>.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Deviation between forecasting and actual power.</p>
</caption>
<graphic xlink:href="fenrg-10-1074916-g012.tif"/>
</fig>
<p>The aforementioned results show that the maximum annual net income and lower solar curtailment rate of 2.89% can be realized with 440MW/880&#xa0;MWh energy storage. Also, in this scenario, the deviation of forecasting and actual power improved largely and the forecasting qualification rate enhanced clearly.</p>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s7">
<title>6 Conclusion</title>
<p>In view of the lack of relevant policies and comprehensive techno-economic indicators when configuring energy storage at the PV power station side, an optimal size determination strategy of energy storage is proposed by considering assessment indicators and the economics from the perspective of PV plant according to the &#x201c;Assessment Rules for New Energy Grid Connection in Jiangsu Province of China&#x201d; in this paper. The feature indicators are automatically extracted from the annual operation data of PV plant, and different scenarios are established for evaluating the reduction of solar curtailment rate and the improvement of forecasting accuracy. The methods proposed in this paper are verified step by step based on the dataset of PV power station, and the following conclusions are summarized from the simulation results.<list list-type="simple">
<list-item>
<p>1) Taking a 4400&#xa0;MW&#xa0;PV plant in Belgium, as example, the results show that when the rated power ratio is 11%, the highest annual net income of energy storage can be achieved. The solar curtailment rate is reduced to 2.89%, and total number of annual unqualified points of PV forecasting power is reduced by 203, which can effectively reduce the penalty cost of PV power plants.</p>
</list-item>
<list-item>
<p>2) The method proposed in this paper is applicable to the extraction of the operating characteristics of large PV power stations within the whole year and realizes the automatic analysis of the optimal sizing of energy storage from the side of PV power station.</p>
</list-item>
<list-item>
<p>3) For specific PV power plant, the size of energy storage should be determined by multidimensional optimization combined with the annual operating characteristics of PV power plants and local assessment rules, in favor of improving the techno-economic indicators of the joint operation of PV power stations and energy storage.</p>
</list-item>
</list>
</p>
<p>In summary, the method proposed in this paper is reasonable for the performance evaluation of large PV power stations with annual operating data and realizes the automatic analysis of the optimal size determination of energy storage systems for PV power stations, which will provide a generalized sizing method of energy storage for PV power stations in different regions.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s8">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s9">
<title>Author contributions</title>
<p>BL: writing and editing, establishing the model, data collation and analysis, and first author. ML: data collation and reviewing. SY, YZ, and BS: survey. JY: writing and methodology, reviewing and editing, and corresponding author.</p>
</sec>
<sec id="s10">
<title>Funding</title>
<p>This work was supported by Technical Standard of Shanghai 2020 &#x201c;Technology Innovation Action Plan&#x201d; (20DZ2205400).</p>
</sec>
<sec sec-type="COI-statement" id="s11">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ali</surname>
<given-names>S. M. H.</given-names>
</name>
<name>
<surname>Lenzen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Tyedmers</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Optimizing 100% renewable grids through shifting residential water&#x2010;heater load</article-title>. <source>Int. J. Energy Res.</source> <volume>43</volume> (<issue>4</issue>), <fpage>1479</fpage>&#x2013;<lpage>1493</lpage>. <pub-id pub-id-type="doi">10.1002/er.4416</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alrawi</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Bayram</surname>
<given-names>I. S.</given-names>
</name>
<name>
<surname>Koc</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Al-Ghamdi</surname>
<given-names>S. G.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Economic viability of rooftop photovoltaic systems and energy storage systems in Qatar</article-title>. <source>Energies</source> <volume>15</volume> (<issue>9</issue>), <fpage>3040</fpage>. <pub-id pub-id-type="doi">10.3390/en15093040</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Apeh</surname>
<given-names>O. O.</given-names>
</name>
<name>
<surname>Meyer</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Overen</surname>
<given-names>O. K.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Contributions of solar photovoltaic systems to environmental and socioeconomic aspects of national development&#x2014;A review</article-title>. <source>Energies</source> <volume>15</volume> (<issue>16</issue>), <fpage>5963</fpage>. <pub-id pub-id-type="doi">10.3390/en15165963</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Azimov</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>Avezova</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Sustainable small-scale hydropower solutions in Central Asian countries for local and cross-border energy/water supply</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>167</volume>, <fpage>112726</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2022.112726</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bird</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lew</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Milligan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Carlini</surname>
<given-names>E. M.</given-names>
</name>
<name>
<surname>Estanqueiro</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Flynn</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>Wind and solar energy curtailment: A review of international experience</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>65</volume>, <fpage>577</fpage>&#x2013;<lpage>586</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2016.06.082</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chao</surname>
<given-names>H.-p.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Efficient pricing and investment in electricity markets with intermittent resources</article-title>. <source>Energy Policy</source> <volume>39</volume> (<issue>7</issue>), <fpage>3945</fpage>&#x2013;<lpage>3953</lpage>. <pub-id pub-id-type="doi">10.1016/j.enpol.2011.01.010</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Power configuration scheme for battery energy storage systems considering the renewable energy penetration level</article-title>. <source>Front. Energy Res.</source> <volume>9</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2021.718019</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dom&#xed;nguez</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Carri&#xf3;n</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Oggioni</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Planning and operating a renewable-dominated European power system under uncertainty</article-title>. <source>Appl. Energy</source> <volume>258</volume>, <fpage>113989</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.113989</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Du</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Jing</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Blockchain based peer-to-peer energy trading between wind power producer and prosumers in short-term market</article-title>. <source>Front. Energy Res.</source> <volume>10</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2022.923292</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Optimal sizing of utility-scale photovoltaic power generation complementarily operating with hydropower: A case study of the world&#x2019;s largest hydro-photovoltaic plant</article-title>. <source>Energy Convers. Manag.</source> <volume>136</volume>, <fpage>161</fpage>&#x2013;<lpage>172</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2017.01.012</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garip</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ozdemir</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimization of PV and battery energy storage size in grid-connected microgrid</article-title>. <source>Appl. Sci.</source> <volume>12</volume> (<issue>16</issue>), <fpage>8247</fpage>. <pub-id pub-id-type="doi">10.3390/app12168247</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jani</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Abdi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Optimal allocation of energy storage systems considering wind power uncertainty</article-title>. <source>J. Energy Storage</source> <volume>20</volume>, <fpage>244</fpage>&#x2013;<lpage>253</lpage>. <pub-id pub-id-type="doi">10.1016/j.est.2018.09.017</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Keeratimahat</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Bruce</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>MacGill</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Analysis of short-term operational forecast deviations and controllability of utility-scale photovoltaic plants</article-title>. <source>Renew. Energy</source> <volume>167</volume>, <fpage>343</fpage>&#x2013;<lpage>358</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2020.11.090</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Vudata</surname>
<given-names>S. P.</given-names>
</name>
<name>
<surname>Bhattacharyya</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lima</surname>
<given-names>F. V.</given-names>
</name>
<name>
<surname>Turton</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Dynamic optimal dispatch of energy systems with intermittent renewables and damage model</article-title>. <source>Mathematics</source> <volume>8</volume> (<issue>6</issue>), <fpage>868</fpage>. <pub-id pub-id-type="doi">10.3390/math8060868</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kut</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Pietrucha-Urbanik</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Tch&#xf3;rzewska-Cie&#x15b;lak</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Reliability-oriented design of a solar-PV deployments</article-title>. <source>Energies</source> <volume>14</volume> (<issue>20</issue>), <fpage>6535</fpage>. <pub-id pub-id-type="doi">10.3390/en14206535</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lei</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Economic dispatch of distribution network with dispersed wind power considering network reconfiguration</article-title>. <source>Front. Energy Res.</source> <volume>10</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2022.942350</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>The capacity allocation method of photovoltaic and energy storage hybrid system considering the whole life cycle</article-title>. <source>J. Clean. Prod.</source> <volume>275</volume>, <fpage>122902</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2020.122902</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Grassi</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Spadacini</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Pignari</surname>
<given-names>S. A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Physically based modeling of hand-assembled wire bundles for accurate EMC prediction</article-title>. <source>IEEE Trans. Electromagn. Compat.</source> <volume>62</volume> (<issue>3</issue>), <fpage>914</fpage>&#x2013;<lpage>922</lpage>. <pub-id pub-id-type="doi">10.1109/temc.2019.2922455</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xiong</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Qian</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Phase locked-loop with decaying DC transient removal for three-phase grids</article-title>. <source>Int. J. Electr. Power &#x26; Energy Syst.</source> <volume>143</volume>, <fpage>108508</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2022.108508</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Sheng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Economic benefit evaluation system of green building energy saving building Technology based on entropy weight method</article-title>. <source>Processes</source> <volume>10</volume> (<issue>2</issue>), <fpage>382</fpage>. <pub-id pub-id-type="doi">10.3390/pr10020382</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lund</surname>
<given-names>P. D.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Capacity matching of storage to PV in a global frame with different loads profiles</article-title>. <source>J. Energy Storage</source> <volume>18</volume>, <fpage>218</fpage>&#x2013;<lpage>228</lpage>. <pub-id pub-id-type="doi">10.1016/j.est.2018.04.030</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luthander</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Nilsson</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Wid&#xe9;n</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Aberg</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Graphical analysis of photovoltaic generation and load matching in buildings: A novel way of studying self-consumption and self-sufficiency</article-title>. <source>Appl. Energy</source> <volume>250</volume>, <fpage>748</fpage>&#x2013;<lpage>759</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.05.058</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lutsenko</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Fetsov</surname>
<given-names>S. S.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Effect of side walls shape on charging and discharging performance of thermal energy storages based on granular phase change materials</article-title>. <source>Renew. Energy</source> <volume>162</volume>, <fpage>466</fpage>&#x2013;<lpage>477</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2020.08.029</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Motalleb</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Reihani</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Ghorbani</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Optimal placement and sizing of the storage supporting transmission and distribution networks</article-title>. <source>Renew. Energy</source> <volume>94</volume>, <fpage>651</fpage>&#x2013;<lpage>659</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2016.03.101</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ndebele</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Assessing the potential for consumer-driven renewable energy development in deregulated electricity markets dominated by renewables</article-title>. <source>Energy Policy</source> <volume>136</volume>, <fpage>111057</fpage>. <pub-id pub-id-type="doi">10.1016/j.enpol.2019.111057</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oh</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Son</surname>
<given-names>S.-Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Theoretical energy storage system sizing method and performance analysis for wind power forecast uncertainty management</article-title>. <source>Renew. Energy</source> <volume>155</volume>, <fpage>1060</fpage>&#x2013;<lpage>1069</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2020.03.170</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pablo-Romero</surname>
<given-names>M. P.</given-names>
</name>
<name>
<surname>S&#xe1;nchez-Braza</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Galyan</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Renewable energy use for electricity generation in transition economies: Evolution, targets and promotion policies</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>138</volume>, <fpage>110481</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2020.110481</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qian</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Spatio-temporal pattern of coupling coordination between urban development and ecological environment under the &#x201c;double carbon&#x201d; goal: A case study in anhui, China</article-title>. <source>Sustainability</source> <volume>14</volume> (<issue>18</issue>), <fpage>11277</fpage>. <pub-id pub-id-type="doi">10.3390/su141811277</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ud-Din Khan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wazeer</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Almutairi</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Alanazi</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Techno-economic analysis of solar photovoltaic powered electrical energy storage (EES) system</article-title>. <source>Alex. Eng. J.</source> <volume>61</volume> (<issue>9</issue>), <fpage>6739</fpage>&#x2013;<lpage>6753</lpage>. <pub-id pub-id-type="doi">10.1016/j.aej.2021.12.025</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Meng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Hydro-thermal-wind-photovoltaic coordinated operation considering the comprehensive utilization of reservoirs</article-title>. <source>Energy Convers. Manag.</source> <volume>198</volume>, <fpage>111824</fpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2019.111824</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimal coordinated bidding strategy of wind and solar system with energy storage in day-ahead market</article-title>. <source>J. Mod. Power Syst. Clean. Energy</source> <volume>10</volume> (<issue>1</issue>), <fpage>192</fpage>&#x2013;<lpage>203</lpage>. <pub-id pub-id-type="doi">10.35833/mpce.2020.000037</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>MPC based control strategy for battery energy storage station in a grid with high photovoltaic power penetration</article-title>. <source>Int. J. Electr. Power &#x26; Energy Syst.</source> <volume>115</volume>, <fpage>105448</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2019.105448</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Qiao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>X.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>IASA-based capacity allocation optimization study for a hydro&#x2013;pumped storage&#x2013;photovoltaic&#x2013;wind complementary clean energy base</article-title>. <source>Front. Energy Res.</source> <volume>10</volume>. <pub-id pub-id-type="doi">10.3389/fenrg.2022.891225</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Dou</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Optimal operation of photothermal power station in regional power grid with high permeability of new energy</article-title>. <source>Sustain. Comput. Inf. Syst.</source> <volume>27</volume>, <fpage>100388</fpage>. <pub-id pub-id-type="doi">10.1016/j.suscom.2020.100388</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Elia Campana</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Bi</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Egusquiza</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Capacity configuration of a hydro-wind-solar-storage bundling system with transmission constraints of the receiving-end power grid and its techno-economic evaluation</article-title>. <source>Energy Convers. Manag.</source> <volume>270</volume>, <fpage>116177</fpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2022.116177</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>