<?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">1260251</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2023.1260251</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>Research on multi-market strategies for virtual power plants with hydrogen energy storage</article-title>
<alt-title alt-title-type="left-running-head">Zhang 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.2023.1260251">10.3389/fenrg.2023.1260251</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Wenyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shen</surname>
<given-names>Yu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Xuanyuan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Ming</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ren</surname>
<given-names>Weixi</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xu</surname>
<given-names>Xiaochuan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Yuyuan</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2380408/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>State Grid Jibei Zhangjiakou Wind and Solar Energy Storage and Transportation New Energy Co., Ltd.</institution>, <addr-line>Zhangjiakou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>State Grid Jibei Electric Power Co.</institution>, <addr-line>Hebei</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>School of Economics and Management</institution>, <institution>North China Electric Power University</institution>, <addr-line>Beijing</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/2339597/overview">Shenbo Yang</ext-link>, Beijing University of Technology, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/787075/overview">Qingkun Tan</ext-link>, State Grid Corporation of China, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2388862/overview">Liu Dunnan</ext-link>, North China Electric Power University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Yuyuan Zhang, <email>120202106034@ncepu.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>03</day>
<month>10</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1260251</elocation-id>
<history>
<date date-type="received">
<day>17</day>
<month>07</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>17</day>
<month>08</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Zhang, Shen, Wang, Li, Ren, Xu and Zhang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Zhang, Shen, Wang, Li, Ren, Xu and Zhang</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 the main body of resource aggregation, Virtual Power Plant (VPP) not only needs to participate in the external energy market but also needs to optimize the management of internal resources. Different from other energy storage, hydrogen energy storage systems can participate in the hydrogen market in addition to assuming the backup supplementary function of electric energy. For the Virtual Power Plant Operator (VPPO), it needs to optimize the scheduling of internal resources and formulate bidding strategies for the electric-hydrogen market based on external market information. In this study, a two-stage model is constructed considering the internal and external interaction mechanism. The first stage model optimizes the operation of renewable energy, flexible load, extraction storage, and hydrogen energy storage system based on the complementary characteristics of internal resources; the second stage model optimizes the bidding strategy to maximize the total revenue of the electricity energy market, auxiliary service market and hydrogen market. Finally, a typical scenario is constructed and the rationality and effectiveness of the strategy are verified. The results show that the hybrid VPP with hydrogen storage has better economic benefits, resource benefits and reliability.</p>
</abstract>
<kwd-group>
<kwd>virtual power plant</kwd>
<kwd>electricity market</kwd>
<kwd>hydrogen market</kwd>
<kwd>hydrogen storage system</kwd>
<kwd>bidding strategy</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Sustainable Energy Systems</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<sec id="s1-1">
<title>1.1 Background and motivation</title>
<p>With the implementation of China&#x2019;s &#x201c;double carbon&#x201d; strategy, new energy sources such as wind power and photovoltaic will see more rapid development, and the penetration rate of new energy sources will continue to increase, which will increase the impact of new energy power fluctuations on the safety and stability of microgrid and its access system and the difficulty of operation and scheduling (<xref ref-type="bibr" rid="B3">Ashish et al., 2021</xref>). On the one hand, the uncertainty of renewable energy power generation puts higher demands on the grid&#x2019;s regulation capacity, and smoothing the volatility and indirectness of renewable energy power generation is an urgent problem to be solved. On the other hand, with the increasing installed capacity of renewable energy, the phenomenon of abandoning wind, light and water is getting more and more serious, and the full consumption of renewable energy and the improvement of energy utilization efficiency is also one of the challenges being faced (<xref ref-type="bibr" rid="B9">Huang et al., 2021</xref>; <xref ref-type="bibr" rid="B13">Liu et al., 2021</xref>). As a smart energy system, VPP has the characteristics of bidirectional trend, which can be used as a controllable power source to strengthen the power supply capacity to the grid and carry out peak regulation, and as a controllable load to increase the power consumption with the system to achieve the valley filling, smooth the system output and demand, and provide a guarantee for the stability of the power system (<xref ref-type="bibr" rid="B18">Nosratabadi et al., 2017</xref>; <xref ref-type="bibr" rid="B21">Saleh Sadeghi et al., 2022</xref>; <xref ref-type="bibr" rid="B4">Chen et al., 2023</xref>). In addition to providing security for the operation of the power system to obtain compensation revenue, VPPs can also participate in the power market at all levels as a flexible resource for capacity, power, and ancillary services to obtain economic benefits through market transactions (<xref ref-type="bibr" rid="B5">Dai et al., 2022</xref>). Energy storage can make up for the inherent defects of new energy in terms of random volatility and fundamentally solve the problem of a high percentage of new energy consumption. The charging and discharging characteristics of energy storage can smooth out the system power fluctuations, improve the new energy consumption capacity, reduce the frequency of power fluctuation impact of a microgrid on its access system, and realize the friendly grid connection of microgrid (<xref ref-type="bibr" rid="B12">Li et al., 2021</xref>). Hydrogen storage enables the smooth operation of power systems through the conversion of hydrogen energy to electrical energy. When electricity is sufficient, the technology of hydrogen production by electrolysis of water is used to make full use of electricity by storing hydrogen; when electricity output is insufficient, the stored hydrogen energy can be used by hydrogen fuel cells or hydrogen combustion turbines to re-generate electricity and return to the system. As a clean and efficient energy source, hydrogen energy storage can play an important role in VPPs (<xref ref-type="bibr" rid="B7">Furat et al., 2022</xref>; <xref ref-type="bibr" rid="B19">Qiu et al., 2022</xref>).</p>
</sec>
<sec id="s1-2">
<title>1.2 Literature review</title>
<p>The research on energy storage in VPPs mainly includes market participation strategy, capacity allocation, optimal scheduling, and benefit allocation. This study focuses on the research from the perspective of market strategy for VPP.</p>
<p>Many scholars have researched the strategy of VPPs participating in the power market. VPPs as flexible resource-rich subjects can play an important role in the electricity market (<xref ref-type="bibr" rid="B22">Shafiekhani et al., 2019</xref>). In addition, with the development of integrated energy systems, multi-energy complementarity has gradually become a trend for resource optimization within VPPs, which can also participate in multi-energy markets as independent subjects (<xref ref-type="bibr" rid="B17">Naughton et al., 2020</xref>; <xref ref-type="bibr" rid="B10">Ju et al., 2022</xref>). The current research on the participation of VPPs in the electricity market mainly includes the design of market mechanisms (<xref ref-type="bibr" rid="B1">Ahmad, 2022</xref>; <xref ref-type="bibr" rid="B16">Morteza et al., 2022</xref>) and the study of trading strategies (<xref ref-type="bibr" rid="B23">Tang and Yang, 2019</xref>; <xref ref-type="bibr" rid="B5">Dai et al., 2022</xref>; <xref ref-type="bibr" rid="B25">Zheng et al., 2022</xref>). For example, Rahimi Mahdi et al. constructed a VPP with wind turbines, PV, conventional generators, energy storage systems, and controllable loads and proposed a strategy for VPP participation in day-ahead and real-time electricity markets considering demand response (<xref ref-type="bibr" rid="B15">Mahdi et al., 2022</xref>). Bo Li et al. proposed a market participation strategy and compared the benefits of VPPs with and without energy storage to participate in electricity energy markets and ancillary services markets (<xref ref-type="bibr" rid="B11">Li and Ghiasi, 2021</xref>). Alahyari, Arman, et al. constructed a VPP consisting of wind power, energy storage, and flexible load, and proposed an optimization strategy considering the stochasticity of renewable energy output and the uncertainty of electricity market price uncertainty (<xref ref-type="bibr" rid="B2">Alahyari et al., 2020</xref>). Henao, MM et al. proposed a bidding strategy for VPP participation in the market while determining the optimal size of the energy storage system (<xref ref-type="bibr" rid="B8">Henao Michelle and Oviedo Jairo Jos&#xe9;, 2022</xref>). Appino, Riccardo R et al. considered the uncertainty of renewable energy output and the volatility of energy prices to optimize the strategy of VPPs with hydrogen energy storage participating in the real-time electricity market (<xref ref-type="bibr" rid="B20">Han et al., 2021</xref>). The above studies take VPP bidding as the research object and put forward a two-tier strategy of internal resource optimization and energy management considering uncertainty. Specifically, these studies only consider VPPS &#x2018;strategies for participating in the electricity market, ignoring the multiple energy attributes of VPPS.</p>
<p>Some scholars have also considered how a multi-energy coupled VPP can participate in a multi-energy market. For example, Zhang, Tao et al. introduced a VPP consisting of a natural gas network, power-to-gas equipment, flexible loads, and energy storage, and studied an optimization model considering dual energy markets by developing different scenarios (<xref ref-type="bibr" rid="B24">Zhang and Hu, 2022</xref>). Ju Liwei et al. constructed a VPP including cogeneration units, wind turbines, power and thermal storage systems, and controllable loads, and considered the impact of different energy market price fluctuations on the profit risk of VPP (<xref ref-type="bibr" rid="B10">Ju et al., 2022</xref>). Liu Xiaoou constructed a VPP operation model including wind turbines, electric vehicles, gas turbines, and controllable loads, and constructed an optimal carbon-electricity integration bidding strategy for VPPs by further analyzing the carbon-electricity integration market characteristics (<xref ref-type="bibr" rid="B14">Liu, 2022</xref>).</p>
<p>Although many studies have considered the strategy of VPP participation in multi-energy markets, the role of electric energy storage in a single electricity market is still only considered for energy storage resources in VPP (mainly electric energy storage resources, excluding thermal and gas storage devices). Unlike electrochemical energy storage, hydrogen energy storage can participate not only in the electric market but also in the hydrogen market. However, few studies have been conducted to evaluate the participation of hydrogen energy storage systems in VPPs. Unlike other energy storage systems, HSSs can not only fully consume the abandoned wind, light, and water resources as energy storage to ensure power supply, but also serve as a source of hydrogen feedstock to enhance the added value within the system and improve the ability of VPPs to participate in external markets (<xref ref-type="bibr" rid="B26">Zheng et al., 2020</xref>; <xref ref-type="bibr" rid="B14">Liu, 2022</xref>).</p>
</sec>
<sec id="s1-3">
<title>1.3 Contribution and research structure</title>
<p>Most of the above studies have considered how VPPs without hydrogen energy storage participate in the power market or multi-energy market, and few studies have explored the role of hydrogen energy storage systems in VPPs and how VPPs with hydrogen energy storage participate in the multi-energy market. In this study, the resource complementary characteristics of renewable energy, flexible load, pumped storage, and hydrogen storage are considered, a two-layer optimization model is constructed, and an external multi-energy market bidding and internal resource optimization strategy are proposed. The structure of the VPP is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Schematic diagram of VPP structure.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g001.tif"/>
</fig>
<p>The VPP consists of renewable energy generating units, pumped storage plants, flexible loads, and HSSs. The VPP operator collects market information and gives operating instructions to the VPP. The VPP optimizes its internal resources through the instructions given by the VPP operator.</p>
<p>Traditionally, the VPP participates in the power market mainly through peak-hour sales on the power side, peak-to-valley arbitrage on the energy storage side, and demand response on the load side to provide ancillary services.</p>
<p>In this study, the VPP participates in the day-ahead energy market and the peaking auxiliary service market and participates in the day-ahead power market as a price receiver. After obtaining the market information (purchase and sale prices, peak and valley hours, <italic>etc.</italic>) of the energy market and the peaking market (peak hours and peak compensation prices for peak filling and peak shaving) released by the power dispatching agency, the VPP coordinates its internal resources according to the bidding strategy. After coordinating internal bidding resources according to the bidding strategy, the company further adjusts the bidding plan in the energy market and the peaking market through optimization and reports the information on the electricity bidding in the energy market and the peaking market to the power dispatching agency. Both HSSs and pumped storage plants can participate in the peaking auxiliary service market, but for HSSs, the total amount of resources for their participation in either the electric auxiliary service market or the hydrogen market is fixed, so decisions need to be made based on prices.</p>
<p>The possible innovations and contributions of this study are as follows.<list list-type="simple">
<list-item>
<p>(1) The application of hydrogen energy storage systems in VPPs is explored.</p>
</list-item>
<list-item>
<p>(2) The participation of VPPs in external markets and the optimization of internal resources are both considered.</p>
</list-item>
<list-item>
<p>(3) The role of a combined clean energy storage strategy for participation in multiple markets is considered for VPP.</p>
</list-item>
<list-item>
<p>(4) The impact of energy market price changes on the market participation strategy of VPPs is evaluated.</p>
</list-item>
</list>
</p>
<p>The rest of this paper is structured as follows: <xref ref-type="sec" rid="s2">Section 2</xref> constructs a two-layer model of VPP participation in the electricity-hydrogen market. Case studies are performed and analyzed in <xref ref-type="sec" rid="s3">Section 3</xref>. <xref ref-type="sec" rid="s4">Section 4</xref> highlights the conclusions.</p>
</sec>
</sec>
<sec id="s2">
<title>2 A model for VPP participation in multiple markets</title>
<sec id="s2-1">
<title>2.1 Bidding strategy for the upper layer of the VPP</title>
<sec id="s2-1-1">
<title>2.1.1 Objective function</title>
<p>The goal for the VPP operator is to maximize total system revenue, which includes revenue from the sale of electricity in the energy market, revenue from the provision of ancillary services, and revenue from the sale of hydrogen in the hydrogen market.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>24</mml:mn>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>Where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the time-of-day price information obtained by the VPP operator in the ancillary services market and the electricity market, respectively. <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the market price of hydrogen, which is a fixed price. <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the quoted quantities of the VPP operator for the ancillary services market, the electricity market, and the hydrogen market, respectively.</p>
<p>Among them, regulating peak auxiliary services are further divided into peak-shaving auxiliary services and valley-filling auxiliary services, and the service prices of peak-shaving auxiliary services and valley-filling auxiliary services are different, therefore, the revenue that VPP operators can obtain by providing auxiliary services can be refined into the following model.<disp-formula id="e2">
<mml:math id="m8">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the valley-filling price and peak-shaving price respectively, <inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the valley-filling capacity of the VPP operator, and <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the peak-shaving capacity of the VPP operator.<disp-formula id="e3">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<inline-formula id="inf11">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf12">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the charging power of the pumped storage plant and HSS involved in valley filling and the load increase of the flexible load involved in valley filling, respectively; <inline-formula id="inf13">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf14">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the discharging power of the pumped storage plant and HSS involved in peak shaving and the load reduction of the flexible load involved in valley filling, respectively.</p>
</sec>
<sec id="s2-1-2">
<title>2.1.2 Constraints</title>
<p>VPPO&#x2019;s bids in the energy and ancillary services markets affect each other, and the limited and fixed flexibility resources available to VPPs require VPP operators to choose the allocation of resources between the energy and ancillary services markets.<disp-formula id="e5">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where <inline-formula id="inf15">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>P</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the system flexibility of the VPP at moment <inline-formula id="inf16">
<mml:math id="m21">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, provided by the lower layer.</p>
</sec>
</sec>
<sec id="s2-2">
<title>2.2 Optimization model for the lower layer of the VPP</title>
<sec id="s2-2-1">
<title>2.2.1 Objective function</title>
<p>For the lower-layer model of the VPP, the objective is to minimize the total system operating cost. Its system operating costs <inline-formula id="inf17">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> include the penalty cost of wind and light abandonment <inline-formula id="inf18">
<mml:math id="m23">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, the penalty cost of load loss <inline-formula id="inf19">
<mml:math id="m24">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, the cost of power purchase <inline-formula id="inf20">
<mml:math id="m25">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>E</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, the cost of equipment start-up and shutdown <inline-formula id="inf21">
<mml:math id="m26">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, and the cost of equipment operation <inline-formula id="inf22">
<mml:math id="m27">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>.<disp-formula id="e6">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mn>24</mml:mn>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>E</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m29">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
<disp-formula id="e8">
<mml:math id="m30">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
<disp-formula id="e9">
<mml:math id="m31">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>E</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
<inline-formula id="inf23">
<mml:math id="m32">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf24">
<mml:math id="m33">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf25">
<mml:math id="m34">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the abandonment penalty tariff, the loss of load penalty tariff, and the purchase tariff, respectively, and <inline-formula id="inf26">
<mml:math id="m35">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf27">
<mml:math id="m36">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf28">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the abandoned power, the lost power, and the purchasing power, respectively.<disp-formula id="e10">
<mml:math id="m38">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
<inline-formula id="inf29">
<mml:math id="m39">
<mml:mrow>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf30">
<mml:math id="m40">
<mml:mrow>
<mml:msubsup>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the positive and negative standby cost coefficients of the system equipment, and <inline-formula id="inf31">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b1;</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the start-stop status of the equipment in the period of <inline-formula id="inf32">
<mml:math id="m42">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. A factor of 0 indicates a shutdown status, and a factor of 1 indicates a start-up status.<disp-formula id="e11">
<mml:math id="m43">
<mml:mrow>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<inline-formula id="inf33">
<mml:math id="m44">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi>p</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> indicates the unit operating cost of the equipment and <inline-formula id="inf34">
<mml:math id="m45">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>O</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> indicates the output of the equipment.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Constraints</title>
<sec id="s2-2-2-1">
<title>2.2.2.1 Power balance constraint</title>
<p>
<disp-formula id="e12">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>g</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>Among them, <inline-formula id="inf35">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf36">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf37">
<mml:math id="m49">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>g</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf38">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf39">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> indicate wind power, photovoltaic power, pumped storage power, fuel cell power, and purchasing power, respectively. <inline-formula id="inf40">
<mml:math id="m52">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf41">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf42">
<mml:math id="m54">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf43">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf44">
<mml:math id="m56">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> indicate the adjusted flexible load, electrolyzer power, pumped storage power, sold power, and abandoned power, respectively.</p>
</sec>
<sec id="s2-2-2-2">
<title>2.2.2.2 Standby capacity constraint</title>
<p>To realize the full consumption of wind power and PV, VPP keeps the corresponding positive and negative standby to cope with the deviation of wind power and PV output through a pumped storage power station and flexible load.<disp-formula id="e13">
<mml:math id="m57">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m58">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>Where <inline-formula id="inf45">
<mml:math id="m59">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf46">
<mml:math id="m60">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the positive deviations of wind power and PV output. <inline-formula id="inf47">
<mml:math id="m61">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf48">
<mml:math id="m62">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the negative deviations between wind power and PV output. <inline-formula id="inf49">
<mml:math id="m63">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf50">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the positive and negative standby of pumped storage plants, <inline-formula id="inf51">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>u</mml:mi>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf52">
<mml:math id="m66">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the positive and negative standby of flexible loads.</p>
</sec>
<sec id="s2-2-2-3">
<title>2.2.2.3 Equipment output constraints</title>
<p>The wind power output constraints are as follows.<disp-formula id="e15">
<mml:math id="m67">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>w</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<inline-formula id="inf53">
<mml:math id="m68">
<mml:mrow>
<mml:msubsup>
<mml:mi>p</mml:mi>
<mml:mi>w</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the rated output power of wind power, <inline-formula id="inf54">
<mml:math id="m69">
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf55">
<mml:math id="m70">
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf56">
<mml:math id="m71">
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf57">
<mml:math id="m72">
<mml:mrow>
<mml:msub>
<mml:mi>v</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, are the access wind speed, rated wind speed and cut-off wind speed and actual wind speed of the system, respectively.</p>
<p>The photovoltaic output constraint is as follows.<disp-formula id="e16">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>I</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
<inline-formula id="inf58">
<mml:math id="m74">
<mml:mrow>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the amount of solar radiation per unit area, <inline-formula id="inf59">
<mml:math id="m75">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the rated installed capacity of PV, and <inline-formula id="inf60">
<mml:math id="m76">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the overall system efficiency, generally taken as 0.8.</p>
<p>The power output constraint of the pumped storage power plant is as follows.<disp-formula id="e17">
<mml:math id="m77">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>g</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
<inline-formula id="inf61">
<mml:math id="m78">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf62">
<mml:math id="m79">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the minimum and maximum technical output of the pumped storage plant under power generation conditions, respectively.<disp-formula id="e18">
<mml:math id="m80">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
<inline-formula id="inf63">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the operating state variable of the pumped storage plant at <inline-formula id="inf64">
<mml:math id="m82">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> time of pumping conditions, and <inline-formula id="inf65">
<mml:math id="m83">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the rated pumping power of the pumped storage plant.<disp-formula id="e19">
<mml:math id="m84">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>g</mml:mi>
</mml:msubsup>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
</p>
<p>Among them, pumped storage power plants are not pumped storage and discharge at the same time.</p>
<p>The hydrogen energy storage system output constraints are as follows.<disp-formula id="e20">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(20)</label>
</disp-formula>
<disp-formula id="e21">
<mml:math id="m86">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(21)</label>
</disp-formula>
<inline-formula id="inf66">
<mml:math id="m87">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the amount of hydrogen produced, <inline-formula id="inf67">
<mml:math id="m88">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the coefficient of hydrogen conversion, and <inline-formula id="inf68">
<mml:math id="m89">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the efficiency of hydrogen conversion. <inline-formula id="inf69">
<mml:math id="m90">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the coefficient of hydrogen to electricity, <inline-formula id="inf70">
<mml:math id="m91">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the efficiency of hydrogen to electricity, and <inline-formula id="inf71">
<mml:math id="m92">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the amount of hydrogen consumed.<disp-formula id="e22">
<mml:math id="m93">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(22)</label>
</disp-formula>
<inline-formula id="inf72">
<mml:math id="m94">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the amount of hydrogen sold per hour, and <inline-formula id="inf73">
<mml:math id="m95">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the amount of hydrogen stored in real-time.</p>
<p>The cumulative hydrogen storage capacity of the hydrogen storage facility satisfies the following constraints.<disp-formula id="e23">
<mml:math id="m96">
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(23)</label>
</disp-formula>
<disp-formula id="e24">
<mml:math id="m97">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(24)</label>
</disp-formula>
<disp-formula id="e25">
<mml:math id="m98">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(25)</label>
</disp-formula>
<disp-formula id="e26">
<mml:math id="m99">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(26)</label>
</disp-formula>
<inline-formula id="inf74">
<mml:math id="m100">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf75">
<mml:math id="m101">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf76">
<mml:math id="m102">
<mml:mrow>
<mml:msubsup>
<mml:mi>H</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the rated power of the electrolyzer, the maximum discharge power of the fuel cell, and the maximum hydrogen storage capacity of the hydrogen storage tank, respectively.<disp-formula id="e27">
<mml:math id="m103">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>e</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
<label>(27)</label>
</disp-formula>
</p>
<p>The electrolysis tank and the fuel cell do not work simultaneously.</p>
</sec>
<sec id="s2-2-2-4">
<title>2.2.2.4 Equipment start/stop constraint</title>
<p>For pumped storage plants and HSSs, frequent start-ups and shutdowns not only increase the cost but also affect the normal use of the equipment and increase the wear and tear of the equipment. Therefore, the number of starts and stops of energy storage systems in a typical operating day needs to be kept within a reasonable range.<disp-formula id="e28">
<mml:math id="m104">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2264;</mml:mo>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(28)</label>
</disp-formula>
<inline-formula id="inf77">
<mml:math id="m105">
<mml:mrow>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> indicates the number of equipment starts and stops, and <inline-formula id="inf78">
<mml:math id="m106">
<mml:mrow>
<mml:msup>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> indicates the maximum number of equipment starts and stops in a typical day.</p>
</sec>
<sec id="s2-2-2-5">
<title>2.2.2.5 Peak shaving constraints</title>
<p>For flexible loads, the following peaking constraints are satisfied.<disp-formula id="e29">
<mml:math id="m107">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(29)</label>
</disp-formula>
<inline-formula id="inf79">
<mml:math id="m108">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the adjustable amount of flexible load, <inline-formula id="inf80">
<mml:math id="m109">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the maximum adjustable amount of flexible load. The flexible load needs to participate in the peaking market according to the period of the peaking auxiliary service market for load adjustment, and the amount of peaking at other time is 0.<disp-formula id="e30">
<mml:math id="m110">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(30)</label>
</disp-formula>
<disp-formula id="e31">
<mml:math id="m111">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x394;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(31)</label>
</disp-formula>
<inline-formula id="inf81">
<mml:math id="m112">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf82">
<mml:math id="m113">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the flexible load filling and peak shaving power respectively, and <inline-formula id="inf83">
<mml:math id="m114">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>v</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf84">
<mml:math id="m115">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the filling and peak shaving flags respectively, which are 0&#x2013;1 variables and are not 1&#xa0;at the same time.</p>
</sec>
<sec id="s2-2-2-6">
<title>2.2.2.6 Loss of load constraint</title>
<p>
<disp-formula id="e32">
<mml:math id="m116">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(32)</label>
</disp-formula>
<inline-formula id="inf85">
<mml:math id="m117">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>P</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> indicates the predicted value of the flexibility load at <inline-formula id="inf86">
<mml:math id="m118">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> time.</p>
</sec>
</sec>
</sec>
<sec id="s2-3">
<title>2.3 Method</title>
<p>In this study, a combination of particle swarm optimization (PSO) algorithm and Cplex solver is used to solve the problem as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>, which is implemented in the MATLAB 2019a platform. The particle swarm algorithm has been widely used as one of the heuristic algorithms for solving optimization problems. Compared with other heuristics such as genetic algorithm (GA) and ant colony optimization (ACO), the PSO algorithm enhances global and local exploration capabilities through a flexible and balanced mechanism. The specific process is as follows.<list list-type="simple">
<list-item>
<p>(1) Initialize the particles, given the initialized velocity and position for each particle in the upper model for the decision variables of the electric energy market bid volume, the auxiliary service market bid volume, and the hydrogen market bid volume.</p>
</list-item>
<list-item>
<p>(2) The lower layer Cplex solver optimizes the internal resources of the VPP according to the lower layer objective function to obtain the optimal power purchase and sale strategy and updates the number of flexibility resources within the VPP to feedback to the upper layer.</p>
</list-item>
<list-item>
<p>(3) The upper-layer particle swarm algorithm updates the bid quantity based on the flexibility resource quantity fed back from the lower layer and calculates the total revenue to evaluate the change in total revenue.</p>
</list-item>
<list-item>
<p>(4) If the model has reached the termination condition, the optimization search process ends and the algorithm terminates; otherwise, the velocity and position of the particles are updated and return to step (2).</p>
</list-item>
</list>
</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Method flow chart.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g002.tif"/>
</fig>
</sec>
</sec>
<sec id="s3">
<title>3 Case study</title>
<sec id="s3-1">
<title>3.1 Parameters and scenario design</title>
<p>Since it is difficult for VPPs to influence market prices in both the electricity and hydrogen markets, this paper sets VPPs as recipients of market prices and makes decisions on quoted quantities based on price information. In this paper, the VPP operator rationalizes the allocation of flexibility resources within the VPP by aggregating them according to the demand of the power dispatching agency and the price information obtained from the power market and the hydrogen market to achieve the maximum total system revenue. Renewable energy generation, flexible loads, pumped storage plants, and HSSs can all participate in the peaking demand of the power market, and HSSs can also participate in the hydrogen market. The peaking periods and the prices of power purchase and sale and auxiliary services are shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. VPPs can participate in the electricity energy market in all 24 time periods.</p>
<p>The price of the electricity market is shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. According to the output characteristics of wind power and PV and load demand characteristics, four typical days in spring, summer, autumn, and winter are selected in this study, as shown in <xref ref-type="fig" rid="F4">Figure 4</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Price parameters.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g003.tif"/>
</fig>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Typical daily output and load.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g004.tif"/>
</fig>
<p>In this paper, the maximum adjustment of the flexible load is set to 20% of the load in that period. The access condition of the auxiliary service market is that the single-day bidding power is not less than 5,000&#xa0;kW. During the peak-shaving period, the maximum power purchase of the VPP is 250&#xa0;kWh. The other parameters of this paper are 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>Parameter table.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Type</th>
<th align="left">Parameters</th>
<th align="left">Value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="5" align="left">Cost parameters</td>
<td align="left">Hydrogen plant construction costs</td>
<td align="left">2,010 CNY/k W</td>
</tr>
<tr>
<td align="left">Hydrogen production equipment operation and maintenance costs</td>
<td align="left">80.4 CNY/k W</td>
</tr>
<tr>
<td align="left">Fuel cell construction costs</td>
<td align="left">4550 CNY/kw</td>
</tr>
<tr>
<td align="left">Fuel cell operation and maintenance costs</td>
<td align="left">182 CNY/kw</td>
</tr>
<tr>
<td align="left">Hydrogen storage equipment costs</td>
<td align="left">65 CNY/kg</td>
</tr>
<tr>
<td rowspan="2" align="left">Equipment operating parameters</td>
<td align="left">P-H efficiency</td>
<td align="left">65%</td>
</tr>
<tr>
<td align="left">P-H-P efficiency</td>
<td align="left">44%</td>
</tr>
<tr>
<td align="left">Price</td>
<td align="left">Hydrogen price</td>
<td align="left">4.5 CNY/Nm<sup>3</sup>
</td>
</tr>
<tr>
<td rowspan="4" align="left">Equipment capacity parameters</td>
<td align="left">Pumped storage plant capacity</td>
<td align="left">300&#xa0;kW</td>
</tr>
<tr>
<td align="left">Electrolytic water hydrogen production plant capacity</td>
<td align="left">100&#xa0;kW</td>
</tr>
<tr>
<td align="left">Hydrogen storage plant capacity</td>
<td align="left">250&#xa0;m<sup>3</sup>
</td>
</tr>
<tr>
<td align="left">Fuel cell capacity</td>
<td align="left">45&#xa0;kW</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Based on the participation of subjects within the VPP and the combination strategy of clean energy storage, four typical scenarios are constructed in this study, as shown in <xref ref-type="table" rid="T2">Table 2</xref> below.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Scenarios design.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th align="center">Wind and PV</th>
<th align="center">Flexible load</th>
<th align="center">Pumped storage</th>
<th align="center">HSSs</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Scenario 1</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">Scenario 2</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
<td align="left"/>
</tr>
<tr>
<td align="left">Scenario 3</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
<td align="left"/>
<td align="left">&#x2713;</td>
</tr>
<tr>
<td align="left">Scenario 4</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
<td align="left">&#x2713;</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2">
<title>3.2 Results</title>
<sec id="s3-2-1">
<title>3.2.1 Method validity verification</title>
<p>In order to quantitatively analyze the effectiveness and applicability of the PSO algorithm to solve the model built in this paper, the GA and ACO algorithms are selected as a comparison, and their running time, iteration times, and result efficiency are compared. Since there are many scenarios involved in this paper, Case 4 (spring) is taken to verify the effectiveness of the proposed solution algorithm as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Comparison results of different algorithms.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g005.tif"/>
</fig>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Analysis of VPP bid results</title>
<sec id="s3-2-2-1">
<title>3.2.2.1 Electric energy market bidding results</title>
<p>
<xref ref-type="fig" rid="F6">Figure 6</xref> shows the bidding results of the electric energy market. The bidding results of scenario 1 and scenario 2 in the electricity energy market in different seasons are consistent. There are two reasons for this: First, pumped storage does not directly participate in the upper grid electric energy market, but provides ancillary services by prioritizing the optimization of resources within the VPP and reducing the flexible load, because the price of purchasing power is higher than the price of sold power. Reducing power purchases is the strategy that achieves the greatest total net benefits. Second, in both Scenario 1 and Scenario 2 renewable energy generation provides electric energy services during normal and peak periods. For pumped storage, the price difference between peak and flat hours does not make it profitable, so renewable energy generation participates in the electricity energy market during normal hours. The bidding results for Scenario 3 and Scenario 4 are the same and lower than Scenario 1 and Scenario 2 in different seasons. After aggregating hydrogen storage, the total revenue obtained from hydrogen production through wind power and PV in the ordinary period is higher than the revenue from electricity sales by participating in the electricity energy market. Therefore, for scenarios 3 and 4, renewable energy generation in the normal period does not participate in the electricity market.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Electric energy market bidding results.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g006.tif"/>
</fig>
</sec>
<sec id="s3-2-2-2">
<title>3.2.2.2 Ancillary services market bidding results</title>
<p>The bidding results of typical days in different seasons under different scenarios are shown in <xref ref-type="fig" rid="F7">Figure 7</xref>. It can be seen that in scenario 1, with only renewable energy generation and flexible load participating in the ancillary services, the benchmark conditions for the ancillary services market can barely be met only in summer. In the other three seasons, VPP cannot participate in the ancillary services market, but only in the electric energy market and the hydrogen market. In the case of aggregated single hydrogen storage (scenario 3), VPP can participate in the ancillary services market only in summer and autumn and barely meets the benchmark conditions in autumn. In spring, it comes close to meeting the auxiliary service benchmark conditions. In contrast, for scenarios 2 and 4, VPP can participate in the ancillary services market in all four seasons, mainly because of the larger capacity of the aggregated pumped storage plant compared to just the aggregated HSS. The regulation capacity is stronger when both pumped storage power plants and hydrogen storage participate in the auxiliary service market.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Ancillary services market bidding results.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g007.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F8">Figure 8</xref> shows the auxiliary service market bidding by specific periods. The positive value indicates participation in peak-shaving auxiliary services, and the negative value indicates participation in valley-filling auxiliary services. In Scenario 1, since the energy storage system is not aggregated, the valley filling capacity is only realized through flexible load demand response, so the auxiliary service capacity is low. The peak-shaving auxiliary service is mainly satisfied by wind power and photovoltaic power output. Since VPP in Scenario 1 only reaches the peak-shaving benchmark in summer, it can only participate in the auxiliary service bidding in summer. Unlike Scenario 1, VPP can participate in auxiliary services in all seasons in Scenario 2. In addition, it can be seen that both the valley filling capacity and peak shaving capacity of Scenario 2 are better than that of Scenario 1. This is because by aggregating pumped storage plants in Scenario 2, VPP can provide valley filling service through pumped storage to purchase electricity from the distribution network in the valley filling auxiliary service phase and reduce the demand for the flexible load to the upper grid through pumped storage in the peak shaving auxiliary service phase. In addition, pumped storage can achieve peak-to-valley arbitrage through peak-to-valley price difference. Compared to Scenario 2, the auxiliary service capability of the VPP alone aggregated HSS in Scenario 3 is inferior. This is because the capacity of the pumped storage plant and the HSS aggregated by the VPP are different, and the capacity of the HSS is only one-third of the pumped storage. However, it can be seen that compared to Scenario 1, after aggregating the HSS, the VPP can participate in the auxiliary service market in both summer and autumn, and its peak-shaving and valley-filling capacity is slightly improved. In scenario 4 where VPP aggregates HSS and pumped storage power plant at the same time, the auxiliary service capacity is significantly increased.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Ancillary services market bidding results by period.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g008.tif"/>
</fig>
</sec>
<sec id="s3-2-2-3">
<title>3.2.2.3 Hydrogen market bidding results</title>
<p>The results of the hydrogen market bidding are shown in <xref ref-type="fig" rid="F9">Figure 9</xref>. It can be seen that VPP does not participate in the hydrogen market bidding in scenario 3 in both summer and fall, which is because the VPP reaches the threshold for bidding in the ancillary services market after aggregating hydrogen storage resources. Since the overall benefits of energy storage participation in the electricity energy market and ancillary services market are higher than those in the hydrogen market, all of the hydrogen storage in this scenario is used to participate in the electricity market by aggregating resources within the VPP. Unlike summer and fall, it is difficult for VPPs with only aggregated HSSs in spring and winter to meet the entry threshold for the ancillary services market. When it is impossible to participate in the auxiliary service market, the benefit of just participating in the electric energy market is lower than the benefit of participating in the hydrogen market, so all the hydrogen is used to participate in the hydrogen market.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Hydrogen market bidding results.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g009.tif"/>
</fig>
<p>Compared to scenario 3, in scenario 4 when VPP aggregates both pumped storage and HSSs, the threshold for ancillary services can be reached in all four seasons. The pumped storage system is prioritized to meet the electricity demand, and the HSS can be used to trade in the hydrogen market after meeting the output of the auxiliary service threshold. Therefore the hydrogen market bidding volume for Scenario 4 is higher than that of Scenario 3.</p>
</sec>
</sec>
<sec id="s3-2-3">
<title>3.2.3 VPP internal resource optimization results</title>
<sec id="s3-2-3-1">
<title>3.2.3.1 Scenario 1</title>
<p>
<xref ref-type="fig" rid="F10">Figure 10</xref> shows the demand response of the flexible load. The higher demand response during peak hours in summer and winter is attributed to the maximum purchasable power designed in this paper. If the flexible load is still larger than the maximum purchasable power after the demand response, the loss of load phenomenon will occur and the energy storage output is needed to compensate.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Demand response results for flexible loads.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g010.tif"/>
</fig>
</sec>
<sec id="s3-2-3-2">
<title>3.2.3.2 Scenario 2</title>
<p>
<xref ref-type="fig" rid="F11">Figure 11</xref> shows the charging and discharging of the pumped storage power plant on a typical day. The pumped storage power plant gives priority to consuming abandoned electricity from renewable energy sources, and if the capacity limit of the pumped storage power plant has not yet been reached, it can buy electricity in the distribution network to achieve peak-valley arbitrage. It can be seen that pumped storage power plants store energy in the valley hours and discharge it in the peak hours. Since pumped storage power plants operate consistently in all seasons, the summer season is used as an example.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Charging and discharging of the pumped storage power station.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g011.tif"/>
</fig>
</sec>
<sec id="s3-2-3-3">
<title>3.2.3.3 Scenario 3</title>
<p>As shown in <xref ref-type="fig" rid="F12">Figure 12</xref>, the operation of the HSS is consistent in spring and winter, and summer and fall. In spring and winter, since the VPP cannot participate in the ancillary services market, the fuel cells of the HSS do not work and all the hydrogen produced is used for storage and traded in the hydrogen market on the second day. In summer and fall, VPP can participate in the ancillary services market, and the hydrogen produced by the HSS is first used by the fuel cell, and the excess hydrogen is sold in the market.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Hydrogen production, storage, and utilization.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g012.tif"/>
</fig>
<p>According to <xref ref-type="fig" rid="F13">Figure 13</xref>, it can be seen that the HSS does not generate electricity in spring and winter. In the summer and autumn peak hours, electricity is generated in the valley, and during flat hours electricity is consumed to produce hydrogen. In addition, the HSS purchases electricity to produce hydrogen during the flat hours in both autumn and winter because the maximum capacity limit of the hydrogen production unit has not yet been reached after consuming the abandoned electricity.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Charging and discharging of hydrogen energy storage system.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g013.tif"/>
</fig>
</sec>
<sec id="s3-2-3-4">
<title>3.2.3.4 Scenario 4</title>
<p>According to <xref ref-type="fig" rid="F14">Figure 14</xref>, during the valley hours, the abandoned electricity from renewable energy sources is used for hydrogen production in priority, and the excess abandoned electricity is used for pumped storage power plants. Due to the capacity of the pumped storage plant, it gives priority to consuming the abandoned power, and the part that has not yet reached the capacity limit can be satisfied by purchasing power in the grid. The renewable energy output in flat hours is used for the HSS and pumped storage power station because the sum of the reduced power purchase cost and auxiliary service revenue is higher than the revenue from power sales in the electricity energy market. During peak hours, the VPP prioritizes the consumption of pumped storage energy, and when the pumped storage energy does not meet the demand, the fuel cell starts to make up for the shortfall. This is because pumped storage power plants can only gain revenue through the electricity market, while HSSs can gain revenue from both auxiliary services and the hydrogen market. The revenue of HSS participating in the electric energy market and auxiliary service market is lower than its participation in the hydrogen market, but the sum of the reduction of the power purchase cost and the increase of auxiliary service revenue due to the generation of HSS in peak hours is higher than the revenue in the hydrogen market, so the HSS chooses to participate in the hydrogen market based on meeting the load demand in peak hours.</p>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Charging and discharging of combined energy storage.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g014.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F15">Figure 15</xref>, the situation tends to be consistent in spring and autumn, with hydrogen consumption for fuel cell generation in the second peak hour, due to the priority use of pumped storage generation and fuel cell start-up to make up for the shortfall when pumped storage cannot meet the load demand. It can be seen that in summer, the fuel cell does not work because the pumped storage power plant output is sufficient to match the flexible load during peak hours. In winter, the fuel cell starts working in the first peak hour because the pumped storage power plant output is difficult to meet the load demand. Compared to the other three seasons, the amount of hydrogen stored in winter is close to zero at the end of a typical day. This is due to the high load demand in winter when hydrogen production is almost entirely used for fuel cells and no excess hydrogen is traded in the hydrogen market.</p>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Hydrogen production, storage, and utilization.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g015.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec id="s3-3">
<title>3.3 Discussion</title>
<sec id="s3-3-1">
<title>3.3.1 Benefits analysis</title>
<sec id="s3-3-1-1">
<title>3.3.1.1 Economic benefits</title>
<p>
<xref ref-type="table" rid="T3">Table 3</xref> shows the total revenue of the VPP for a typical day in spring. 1) In terms of revenue in the electricity market, scenario 3 has the lowest revenue because both scenario 1 and scenario 3 do not reach the peaking benchmark of the ancillary services market and cannot participate in the ancillary services market but only in the electricity energy market. For scenario 3, the revenue of renewable energy generation for hydrogen production in flat hours is higher than the revenue from participating in the electricity energy market, so the revenue of the electricity energy market for scenario 3 is lower than that of scenario 1. 2) The hydrogen revenue of scenario 3 is greater than that of scenario 4. In scenario 3, VPP has not reached the threshold of the electricity auxiliary service market, and the hydrogen revenue from HSS is greater than the revenue from electricity sales of fuel cells based on meeting the load demand. As for scenario 4, because VPP has reached the auxiliary service benchmark condition, the revenue from auxiliary service is greater than the revenue from hydrogen sales, so the excess hydrogen from HSS is only used for trading in the hydrogen market on the premise of maximizing the demand for auxiliary service that can be participated. 3) The power purchase cost of scenario 3 is higher than that of scenario 2 because the capacity of HSS is lower than that of pumped storage. 4) Compared to Scenario 2 and Scenario 4, Scenario 3 still has the abandonment penalty due to the capacity of the HSS. 5) After aggregating the energy storage system, there is no load loss in the VPP. The revenue for a typical winter day is shown in <xref ref-type="table" rid="T4">Table 4</xref> and is similar to that of spring.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Total bidding revenue for a typical spring day/CNY.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center"/>
<th align="center">Electricity market revenue</th>
<th align="center">Hydrogen market revenue</th>
<th align="center">Power purchase cost</th>
<th align="center">Abandonment penalty</th>
<th align="center">Loss of load penalty</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Scenario 1</td>
<td align="center">2,291.00</td>
<td align="center">-</td>
<td align="center">4,203.36</td>
<td align="center">227.57</td>
<td align="center">90.85</td>
</tr>
<tr>
<td align="center">Scenario 2</td>
<td align="center">5,017.65</td>
<td align="center">-</td>
<td align="center">2,303.97</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 3</td>
<td align="center">2028.35</td>
<td align="center">834.20</td>
<td align="center">3,467.58</td>
<td align="center">67.62</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 4</td>
<td align="center">5,479.00</td>
<td align="center">657.75</td>
<td align="center">2082.45</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Total bidding revenue for a typical winter day/CNY.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center"/>
<th align="center">Electricity market revenue</th>
<th align="center">Hydrogen market revenue</th>
<th align="center">Power purchase cost</th>
<th align="center">Abandonment penalty</th>
<th align="center">Loss of load penalty</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Scenario 1</td>
<td align="center">2,113.84</td>
<td align="center">-</td>
<td align="center">6,204.31</td>
<td align="center">160.85</td>
<td align="center">480.68</td>
</tr>
<tr>
<td align="center">Scenario 2</td>
<td align="center">4,765.59</td>
<td align="center">-</td>
<td align="center">2,716.25</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 3</td>
<td align="center">1861.18</td>
<td align="center">194.54</td>
<td align="center">3,743.61</td>
<td align="center">20.63</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 4</td>
<td align="center">5,275.20</td>
<td align="center">39.84</td>
<td align="center">2,436.50</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>According to <xref ref-type="table" rid="T5">Tables 5</xref>, <xref ref-type="table" rid="T6">6</xref>, the situation is similar in summer and autumn. 1) In the summer of Scenario 1, VPP can participate in the auxiliary service market so its revenue in the electricity market is higher than that in the autumn; 2) Compared with the spring, in the summer and autumn of Scenario 3, VPP can participate in the auxiliary service market so its revenue is higher than that in Scenario 1, but its revenue is lower than that in Scenario 2 due to the capacity of the HSS; 3) There is no revenue in the hydrogen market in both the summer and fall, which is because the sum of the power purchase cost to compensate for the reduction of flexible load demand and the revenue from ancillary services is higher than the income in the hydrogen market, so all the hydrogen is used for fuel cells.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Total bidding revenue for a typical summer day/CNY.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center"/>
<th align="center">Electricity market revenue</th>
<th align="center">Hydrogen market revenue</th>
<th align="center">Power purchase cost</th>
<th align="center">Abandonment penalty</th>
<th align="center">Loss of load penalty</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Scenario 1</td>
<td align="center">3,455.54</td>
<td align="center">-</td>
<td align="center">6,129.52</td>
<td align="center">236.29</td>
<td align="center">377.62</td>
</tr>
<tr>
<td align="center">Scenario 2</td>
<td align="center">5,489.15</td>
<td align="center">-</td>
<td align="center">2,413.66</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 3</td>
<td align="center">4,210.09</td>
<td align="center">0</td>
<td align="center">3,642.17</td>
<td align="center">75.40</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 4</td>
<td align="center">5,996.65</td>
<td align="center">953.97</td>
<td align="center">2,181.61</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Total bidding revenue for a typical autumn day/CNY.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center"/>
<th align="center">Electricity market revenue</th>
<th align="center">Hydrogen market revenue</th>
<th align="center">Power purchase cost</th>
<th align="center">Abandonment penalty</th>
<th align="center">Loss of load penalty</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Scenario 1</td>
<td align="center">2,363.97</td>
<td align="center">-</td>
<td align="center">5,345.60</td>
<td align="center">227.57</td>
<td align="center">56.83</td>
</tr>
<tr>
<td align="center">Scenario 2</td>
<td align="center">5,081.84</td>
<td align="center">-</td>
<td align="center">4,683.79</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 3</td>
<td align="center">3,741.48</td>
<td align="center">0</td>
<td align="center">3,467.58</td>
<td align="center">67.62</td>
<td align="center">0</td>
</tr>
<tr>
<td align="center">Scenario 4</td>
<td align="center">5,549.17</td>
<td align="center">641.40</td>
<td align="center">2082.44</td>
<td align="center">0</td>
<td align="center">0</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-3-1-2">
<title>3.3.1.2 Resource benefits</title>
<p>Restricted by the electricity market, wind power, and PV can only provide auxiliary services during peak-shaving service hours. During valley hours and flat hours, there will be power abandonment because the superior grid has a strong supply capacity itself and does not need additional power output. According to <xref ref-type="fig" rid="F16">Figure 16</xref>, there is no power abandonment in Scenario 2 and Scenario 4, because compared with Scenario 1, VPP aggregates energy storage resources and can fully utilize the resources of wind and PV. In scenario 3, the VPP still has the abandonment phenomenon, which is limited by the capacity of the aggregated hydrogen storage, but the abandonment phenomenon is significantly alleviated compared to scenario 1.</p>
<fig id="F16" position="float">
<label>FIGURE 16</label>
<caption>
<p>Abandoned power.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g016.tif"/>
</fig>
</sec>
<sec id="s3-3-1-3">
<title>3.3.1.3 Reliability benefits</title>
<p>As shown in <xref ref-type="fig" rid="F17">Figure 17</xref>, the VPP will have a loss of load only in Scenario 1. Due to the limitation of the power available to be purchased from the grid, there will be load loss during peak hours because the VPP load demand is higher than the power available to be purchased. In other scenarios, there is no load loss because the VPP aggregates energy storage resources and can cover its power shortfall.</p>
<fig id="F17" position="float">
<label>FIGURE 17</label>
<caption>
<p>Loss of load.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g017.tif"/>
</fig>
</sec>
</sec>
<sec id="s3-3-2">
<title>3.3.2 Sensitivity analysis</title>
<p>
<xref ref-type="fig" rid="F18">Figure 18</xref> shows the sensitivity analysis of the changes in electricity and hydrogen prices. The market participation strategy of the VPP changes as a result of price changes in the electricity and ancillary services markets and the hydrogen market. In Scenario 1 and Scenario 2, a single change in electricity market price causes a linear change in total revenue because no hydrogen market is involved. In scenarios 3 and 4, total revenue does not show a linear change due to changes in the electricity market price as well as the hydrogen market price. The VPP chooses the optimal strategy that maximizes the total revenue due to the changes in hydrogen and electricity prices, which implies a change in the proportion of hydrogen resources participating in the electricity and hydrogen markets based on the existing prices. In addition, as the prices in the electricity and hydrogen markets increase simultaneously, the total revenue does not always increase, but peaks at a certain price. This is because the amount of hydrogen resources allocated in the electricity market and the hydrogen market does not change when the electricity price and the hydrogen price increase in the same proportion. But when the price changes in the electricity market and the hydrogen market in different proportions, hydrogen resources do not necessarily prioritize the demand for ancillary services, but rather tilt more resources to the higher-priced hydrogen market.</p>
<fig id="F18" position="float">
<label>FIGURE 18</label>
<caption>
<p>Sensitivity analysis of the prices.</p>
</caption>
<graphic xlink:href="fenrg-11-1260251-g018.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>4 Conclusion</title>
<p>This paper investigates the bidding strategy for a VPP in the electricity and hydrogen markets and constructs a two-layer model that considers the bidding strategy and internal resource operation optimization strategy based on the complementary characteristics of renewable power, flexible load, pumped storage, and hydrogen storage resources. The internal resources are optimally dispatched and the optimized flexibility capability is fed back to the VPP operator to adjust the bidding strategy until the optimal bidding strategy is achieved, taking into account the objectives of maximizing the total revenue of the VPP to participate in the external market and minimizing the internal abandonment penalty and load loss penalty. In this paper, the particle swarm algorithm combined with the Cplex solver is used for solving the problem, and the results show that.<list list-type="simple">
<list-item>
<p>(1) The inside-outside two-layer optimization model constructed in this paper can guarantee that the VPP makes optimal bidding decisions based on information from the external power and hydrogen markets and optimally adjusts internal resources based on the bidding situation to achieve the goal of maximizing total net benefits.</p>
</list-item>
<list-item>
<p>(2) In the current market environment, pumped storage power plants can achieve peak-to-valley arbitrage by interacting with the grid. However, for HSSs, it is uneconomical to compensate for peak hour load demand by purchasing power in the valley hours. However, in the case of the hydrogen market, HSS can earn hydrogen sales by purchasing power from the grid to produce hydrogen in the valley and the flat hours.</p>
</list-item>
<list-item>
<p>(3) In the absence of a benchmark threshold for ancillary services, the benefits of HSSs to reduce the cost of power purchase by reducing the demand for flexible loads are smaller than the benefits of participating in the hydrogen market but larger than the benefits of participating in the electricity market. However, due to the threshold of auxiliary services, HSS will sacrifice part of the benefits of direct participation in the hydrogen market and give priority to participating in the electricity auxiliary services market.</p>
</list-item>
<list-item>
<p>(4) For a VPP with multiple complementary energy storage, it is more profitable to fully consider the electricity-hydrogen market and the multi-functional properties of hydrogen storage than to merely equate the HSS with other energy storage as a backup power source. (<xref ref-type="bibr" rid="B6">Emmanouil et al., 2022</xref>).</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<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="s6">
<title>Author contributions</title>
<p>WZ: Conceptualization, Writing&#x2013;review and editing. YS: Writing&#x2013;original draft. XW: Writing&#x2013;original draft. ML: Investigation, Methodology, Writing&#x2013;original draft. WR: Data curation, Software, Writing&#x2013;review and editing. XX: Software, Data curation, Writing&#x2013;original draft. YZ: Visualization, Supervision, Writing&#x2013;original draft.</p>
</sec>
<sec id="s7">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of interest</title>
<p>WZ, YS, ML, WR, and XX were employed by State Grid Jibei Zhangjiakou Wind and Solar Energy Storage and Transportation New Energy Co., Ltd. XW was employed by State Grid Jibei Electric Power Co.</p>
<p>The remaining author declares 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>
<p>The Reviewer LD declared a shared affiliation with the author YZ at the time of the review.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<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>Ahmad</surname>
<given-names>R. J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A stochastic model for participation of virtual power plants in futures markets, pool markets and contracts with withdrawal penalty</article-title>. <source>J. Energy Storage</source>, <volume>50</volume>, <fpage>104334</fpage>. <pub-id pub-id-type="doi">10.1016/j.est.2022.104334</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alahyari</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ehsan</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Moghimi</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Managing distributed energy resources DERs through virtual power plant technology VPP: A stochastic information-gap decision theory igdt approach</article-title>. <source>Iran. J. Sci. Technol. Trans. Electr. Eng.</source> <volume>44</volume> (<issue>1</issue>), <fpage>279</fpage>&#x2013;<lpage>291</lpage>. <pub-id pub-id-type="doi">10.1007/s40998-019-00248-w</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ashish</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Myron</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Dmitrii</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Eugene</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Joey</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Christian</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Transition pathway towards 100% renewable energy across the sectors of power, heat, transport, and desalination for the Philippines</article-title>. <source>Renew. Sustain. Energy Rev.</source>, <volume>144</volume>. <fpage>110934</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2021.110934</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Optimal bidding strategy of renewable-based virtual power plant in the day-ahead market</article-title>. <source>Int. J. Electr. Power Energy Syst.</source>, <volume>239</volume> <fpage>122379</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.122379</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>An equilibrium model of the electricity market considering the participation of virtual power plants</article-title>. <source>Energy Rep.</source> <volume>8</volume> (<issue>13</issue>), <fpage>506</fpage>&#x2013;<lpage>514</lpage>. <pub-id pub-id-type="doi">10.1016/j.egyr.2022.08.096</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Emmanouil</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Nikolaos</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Theocharis</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Challenges of integrating hydrogen energy storage systems into nearly zero-energy ports</article-title>. <source>Energy</source>, <volume>241</volume>. <fpage>122878</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.122878</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Furat</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Martin</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Shafiullah</surname>
<given-names>G. M.</given-names>
</name>
</person-group> <source>Hydrogen production for energy: An overview</source>. <publisher-name>Elsevier BV</publisher-name>, <publisher-loc>Amsterdam, Netherlands</publisher-loc>, <year>2022</year>, <fpage>3847</fpage>&#x2013;<lpage>3869</lpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Henao Michelle</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Oviedo Jairo Jos&#xe9;</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Bidding strategy for VPP and economic feasibility study of the optimal sizing of storage systems to face the uncertainty of solar generation modelled with IGDT</article-title>. <source>Energies</source> <volume>15</volume> (<issue>3</issue>), <fpage>953</fpage>. <pub-id pub-id-type="doi">10.3390/en15030953</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A two-stage planning and optimization model for water - hydrogen integrated energy system with isolated grid</article-title>. <source>J. Clean. Prod.</source> <volume>313</volume>, <fpage>127889</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2021.127889</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ju</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Nearly-zero carbon optimal operation model and benefit allocation strategy for a novel virtual power plant using carbon capture, power-to-gas, and waste incineration power in rural areas</article-title>. <source>Appl. Energy</source> <volume>310</volume>, <fpage>118618</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2022.118618</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Ghiasi</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A new strategy for economic virtual power plant utilization in electricity market considering energy storage effects and ancillary services</article-title>. <source>J. Electr. Eng. Technol.</source> <volume>16</volume> (<issue>6</issue>), <fpage>2863</fpage>&#x2013;<lpage>2874</lpage>. <pub-id pub-id-type="doi">10.1007/s42835-021-00811-8</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Energy management for hybrid energy storage system in electric vehicle: A cyber-physical system perspective</article-title>. <source>Energy</source> <volume>230</volume>, <fpage>120890</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.120890</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Miller</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Cholette Michael</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Gerard</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Glenn</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A multi-dimension clustering-based method for renewable energy investment planning</article-title>. <source>Renew. Energy</source> <volume>172</volume>, <fpage>651</fpage>&#x2013;<lpage>666</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2021.03.056</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Research on bidding strategy of virtual power plant considering carbon-electricity integrated market mechanism</article-title>. <source>Int. J. Electr. Power Energy Syst.</source> <volume>137</volume>, <fpage>107891</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2021.107891</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mahdi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Fatemeh Jahanbani</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Oluwaseyi</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Ali Jahanbani</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Two-stage interval scheduling of virtual power plant in day-ahead and real-time markets considering compressed air energy storage wind turbine</article-title>. <source>J. Energy Storage</source>, <volume>45</volume>. <fpage>103599</fpage>. <pub-id pub-id-type="doi">10.1016/j.est.2021.103599</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Morteza</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Abdollah</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Omid</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Miadreza</surname>
<given-names>S.-K.</given-names>
</name>
<name>
<surname>Catal&#xe3;o Jo&#xe3;o</surname>
<given-names>P. S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimal bidding strategy of a renewable-based virtual power plant including wind and solar units and dispatchable loads</article-title>. <source>Energy</source> <volume>239</volume> <fpage>122379</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2021.122379</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Naughton</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Mancarella</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Optimization of multi-energy virtual power plants for providing multiple market and local network services</article-title>. <source>Electr. Power Syst. Res.</source> <volume>189</volume> (<issue>189</issue>), <fpage>106775</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2020.106775</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nosratabadi</surname>
<given-names>S. M.</given-names>
</name>
<name>
<surname>Hooshmand</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Gholipour</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>A comprehensive review on microgrid and virtual power plant concepts employed for distributed energy resources scheduling in power systems</article-title>. <source>Renew. Sustain. ENERGY Rev.</source> <volume>67</volume>, <fpage>341</fpage>&#x2013;<lpage>363</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2016.09.025</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qiu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimal planning of Cross-regional hydrogen energy storage systems considering the uncertainty</article-title>. <source>Appl. Energy</source> <volume>326</volume>, <fpage>119973</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2022.119973</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>&#xc1;ngel</surname>
<given-names>J. O. G.</given-names>
</name>
</person-group>, (<year>2021</year>). <article-title>Energy-based stochastic MPC for integrated electricity-hydrogen VPP in real-time markets</article-title>. <source>Electr. Power Syst. Res.</source> <volume>195</volume> <fpage>106738</fpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2020.106738</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Saleh Sadeghi</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Majid</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Amirhossein</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Hamidreza</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Masoud Aliakbar</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Economic assessment of multi-operator virtual power plants in electricity market: A game theory-based approach</article-title>. <source>Sustain. Energy Technol. Assessments</source> <volume>53</volume>, <fpage>102733</fpage>. <pub-id pub-id-type="doi">10.1016/j.seta.2022.102733</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shafiekhani</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Badri</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Catal&#xe3;o</surname>
<given-names>J. P. S.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Strategic bidding of virtual power plant in energy markets: A Bi-level multi-objective approach</article-title>. <source>Int. J. Electr. Power &#x26; Energy Syst.</source> <volume>113</volume> <fpage>208</fpage>-<lpage>219</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijepes.2019.05.023</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>W.-J.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>H.-T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Optimal operation and bidding strategy of a virtual power plant integrated with energy storage systems and elasticity demand response</article-title>. <source>IEEE Access</source> <volume>7</volume>, <fpage>79798</fpage>&#x2013;<lpage>79809</lpage>. <pub-id pub-id-type="doi">10.1109/access.2019.2922700</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimal scheduling strategy of virtual power plant with power-to-gas in dual energy markets</article-title>. <source>IEEE Trans. Industry Appl.</source> <volume>58</volume> (<issue>2</issue>), <fpage>2921</fpage>&#x2013;<lpage>2929</lpage>. <pub-id pub-id-type="doi">10.1109/TIA.2021.3112641</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Mei</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A peer-to-peer energy trading market embedded with residential shared energy storage units</article-title>. <source>Appl. Energy</source>, <volume>308</volume>. <fpage>118400</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2021.118400</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Development of renewable energy multi-energy complementary hydrogen energy system (A case study in China): A review</article-title>. <source>Energy Explor. Exploitation</source> <volume>38</volume> (<issue>6</issue>), <fpage>2099</fpage>&#x2013;<lpage>2127</lpage>. <pub-id pub-id-type="doi">10.1177/0144598720953512</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>