<?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">1373588</article-id>
<article-id pub-id-type="doi">10.3389/fenrg.2024.1373588</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 joint dispatch of wind, solar, hydro, and thermal power based on pumped storage power stations</article-title>
<alt-title alt-title-type="left-running-head">Jia et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenrg.2024.1373588">10.3389/fenrg.2024.1373588</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Jia</surname>
<given-names>Jun</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/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<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>Guangming</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhou</surname>
<given-names>Xiaoxiong</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shi</surname>
<given-names>Zhihan</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhu</surname>
<given-names>Mingxiang</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Lv</surname>
<given-names>Xiaodong</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2392877/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Transportation Engineering</institution>, <institution>Nanjing Tech University</institution>, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>College of Electrical Engineering and Control Science</institution>, <institution>Nanjing Tech University</institution>, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Taizhou College</institution>, <institution>Nanjing Normal University</institution>, <addr-line>Taizhou</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/2365354/overview">Minglei Bao</ext-link>, Zhejiang University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1902289/overview">Mojtaba Nedaei</ext-link>, University of Padua, Italy</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1472867/overview">Helena Mart&#xed;n</ext-link>, Universitat Politecnica de Catalunya, Spain</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Guangming Zhang, <email>zgm@njtech.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>03</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1373588</elocation-id>
<history>
<date date-type="received">
<day>20</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>07</day>
<month>03</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Jia, Zhang, Zhou, Shi, Zhu and Lv.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Jia, Zhang, Zhou, Shi, Zhu and Lv</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>In the context of energy conservation and emission reduction, the integration and consumption of large-scale wind and solar resources is an inevitable trend in future energy development. However, with the increase of wind and solar grid-connected capacity, the power system also requires more flexible resources to ensure safe operation. To enhance the economic efficiency of the complementary operation of wind, solar, hydro, and thermal sources, considering the peak regulation characteristics of different types of power sources, the study of the joint dispatch model of complementary utilization of various generation methods like wind, solar, hydro, thermal, and storage is of great significance for the economic dispatch of the power system. Existing studies mainly focus on traditional thermal power units or hydropower units, with few studies investigating the impact of pumped-storage power stations on the absorption of renewable energy. Firstly, this paper introduces the composition and function of each unit under the research framework and establishes a joint dispatch model for wind, solar, hydro, and thermal power. Secondly, the paper elaborates on the objective function within the model, mainly covering the operating costs of thermal power units, hydropower units, pumped storage, wind and solar units, the cost of discarding new energy, and the cost of load shedding. Subsequently, the paper presents the constraints of the system model, mainly the feasible boundaries for the operation of each unit within the system. Finally, The results of the calculations show that the proposed model reduces the total operating cost by 12% and the power abandonment rate by 82% compared to the conventional model. It is shown that the proposed model can not only significantly improve the economic efficiency of the system operation but also reduce the level of energy waste and load shedding, effectively enhancing the degree of energy utilization within the system.</p>
</abstract>
<kwd-group>
<kwd>multi-energy complementarity</kwd>
<kwd>hydropower unit</kwd>
<kwd>joint dispatch</kwd>
<kwd>pumped storage</kwd>
<kwd>renewable energy absorption</kwd>
<kwd>power system optimization</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Process and Energy Systems Engineering</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Developing a new power system adapted to the increasing proportion of new energy sources is a crucial measure for China to achieve its carbon peak and carbon neutrality goals on schedule and is essential for ensuring national energy security (<xref ref-type="bibr" rid="B14">Hou et al., 2023</xref>). Traditional thermal power units generate a large amount of carbon dioxide and other greenhouse gases when burning coal, significantly impacting global warming. Additionally, thermal power stations consume vast amounts of water resources during operation and produce air pollutants and solid waste, adversely affecting the environment and human health. The medium-to-long-term seasonal fluctuations, short-term randomness, and intermittent fluctuations of wind and solar power generation pose severe flexibility challenges for the new power system across daily, monthly, seasonal, and annual timescales, necessitating the exploration of new renewable flexibility resources (<xref ref-type="bibr" rid="B27">Nnamchi et al., 2023</xref>). Hydropower ranks second in installed capacity in China, with a 16.1% share in 2022, second only to thermal power units. Its strong regulation capability, combined with the random fluctuations of wind and solar power, forms a complementary system that outputs relatively smooth and stable high-quality power, effectively solving the challenges of wind and solar energy development (<xref ref-type="bibr" rid="B3">Bello et al., 2023</xref>). Prioritizing the use of existing power sources to achieve multi-energy complementarity and maximize the use of renewable energy is of significant research value (<xref ref-type="bibr" rid="B21">Lamb et al., 2023</xref>). Therefore, the functional positioning of hydropower should gradually shift from focusing on electricity generation to capacity support, thus promoting the integrated development of hydro-wind-solar complementarity. However, with the increasing capacity of wind and solar power, the issue of abandoning wind and solar energy is unavoidable, and conventional hydropower cannot effectively store the electricity generated from abandoned wind and solar power (<xref ref-type="bibr" rid="B17">Jin et al., 2023</xref>). Pumped storage, as a crucial technology for enhancing the absorption level of new energy, has developed rapidly in China (<xref ref-type="bibr" rid="B10">Garcia-Gonzalez et al., 2008</xref>; <xref ref-type="bibr" rid="B9">Feng et al., 2021</xref>; <xref ref-type="bibr" rid="B16">Huang et al., 2023</xref>; <xref ref-type="bibr" rid="B23">Liu et al., 2023</xref>; <xref ref-type="bibr" rid="B32">Wang et al., 2023</xref>). How to use pumped storage technology efficiently and economically to promote the absorption of new energy is an important direction of current research.</p>
<p>When the penetration rate of wind power increases to a certain extent, relying solely on thermal power to cope with the uncertainty of wind and solar output will lead to frequent starting and stopping of thermal power units, threatening the safety, stability, and economy of the power grid operation (<xref ref-type="bibr" rid="B35">Ye et al., 2023</xref>). In the analysis of wind and solar grid integration, research on the active output characteristics of the system mainly includes studies on the operating characteristics of wind-solar-thermal systems and the spatiotemporal complementarity of active power output from wind and solar energy (<xref ref-type="bibr" rid="B31">Sun and Harrison, 2019</xref>). Literature (<xref ref-type="bibr" rid="B30">Silva et al., 2016</xref>) used Pearson correlation coefficients to analyze the correlation between different power sources. Literature (<xref ref-type="bibr" rid="B13">Hong-Mei et al., 2013</xref>) assessed the natural characteristics of wind and solar output and the complementarity of wind, solar, and storage from the aspects of complementarity rate and smoothness. Literature (<xref ref-type="bibr" rid="B5">Cuiping et al., 2017</xref>) evaluated the operating characteristics of the photovoltaic-hydropower complementary system based on indicators such as the abandoned light ratio, the ratio of thermal power to load, and grid-connected revenue. Facing the numerous uncertainties of wind and photovoltaic grid integration, literature (<xref ref-type="bibr" rid="B4">Canizes et al., 2012</xref>) used the Monte Carlo method to simulate uncertain factors such as wind speed and natural inflow. Literature (<xref ref-type="bibr" rid="B12">Hinojosa and Velasquez, 2016</xref>) used scenario trees for uncertainty modeling, which includes a large amount of random information. Literature (<xref ref-type="bibr" rid="B2">Azizipanah-Abarghooee et al., 2012</xref>), based on the Monte Carlo concept, used the point estimation method to obtain the probabilistic statistical information of the variables in question.</p>
<p>Currently, hydropower, as the most widely used clean energy source, is also a hot topic in the research of multi-energy system optimization dispatch. For efficient use of the regulation capability of hydropower to satisfy the power balance, literature (<xref ref-type="bibr" rid="B8">Esmaeily et al., 2016</xref>) specifically conducted dispatch analysis for hydropower systems. Literature (<xref ref-type="bibr" rid="B1">Arce et al., 2002</xref>) considers the total power generation of hydropower stations as a constraint, taking into account the impact of hydropower units on system peak regulation, and optimizes the output of each unit in the system. Literature (<xref ref-type="bibr" rid="B11">Gromyko et al., 2023</xref>) included the water volume constraint of the hydropower station in the modeling of hydropower stations. Literature (<xref ref-type="bibr" rid="B36">Zanoli et al., 2023</xref>) used a piecewise linear function to fit the hydropower conversion curve. Literature (<xref ref-type="bibr" rid="B28">Rahman et al., 2022</xref>) analyzed the complementarity of wind and hydropower output on a time scale and established a linear programming model for the wind-hydro joint dispatch system. Literature (<xref ref-type="bibr" rid="B34">Yang et al., 2023</xref>) proposed a short-term stochastic optimization dispatch model for wind, water, and thermal multi-energy systems. Large-scale hydropower generation bases are mostly developed and utilized on the scale of river basin cascade hydropower station groups, with a high degree of coupling between cascading hydropower stations. To improve the overall efficiency of basin cascade hydropower, the operation of cascading hydropower stations needs to consider constraints such as water volume, storage capacity, and head in the upstream and downstream of the basin. However, traditional hydropower units, as one of the long-standing renewable energy technologies, hold an important position in the global energy structure. But this traditional method faces some challenges and shortcomings. Firstly, the energy output of traditional hydropower units is limited by the availability of water resources, especially in dry seasons or areas with scarce water resources, significantly affecting their power generation capacity. Additionally, the construction and operation of hydropower stations may have negative impacts on the local ecosystem and environment.</p>
<p>In this context, the development of pumped storage technology offers a new perspective. Pumped storage power stations, as an efficient method of energy storage, can store energy when electricity demand is low and release it during peak periods, thus optimizing energy allocation and utilization. This not only enhances the stability and reliability of the power grid but also provides an effective solution for the integration of renewable energy sources.</p>
<p>Pumped storage power stations are currently the world&#x2019;s most comprehensively evaluated large-scale electricity storage technology. Their basic principle involves using surplus electricity to pump water from lower elevations to higher reservoirs when electricity demand is low; then, when electricity demand is high, releasing this water to drive turbines and generate electricity. This technology is not only crucial for the stable operation of the power grid but also plays an important role in the integration of renewable energy sources. Early pumped storage power stations were mainly concentrated in Europe and North America, but in recent years, development in Asia, especially in China, has been rapid. Literature (<xref ref-type="bibr" rid="B18">Kanakasabapathy and Swarup, 2010</xref>) notes that China has made significant efforts in the development of pumped storage technology and has now become the world&#x2019;s largest pumped storage market. The efficiency of pumped storage power stations is affected by various factors, including hydraulic, mechanical, and electrical losses. Literature (<xref ref-type="bibr" rid="B15">Hu et al., 2012</xref>) points out that optimizing the design of turbines and pumps, as well as adopting advanced control strategies, can significantly improve the overall efficiency of pumped storage power stations. Although pumped storage power stations are a form of clean energy, their construction and operation have an impact on the environment that cannot be ignored. With the rapid development of renewable energy, the role of pumped storage power stations in the electricity market is becoming increasingly important. Literature (<xref ref-type="bibr" rid="B19">Katsuhiro et al., 2013</xref>) indicates that market mechanism reforms, such as implementing time-of-use electricity pricing policies, will help improve the economic viability of pumped storage power stations. At the same time, facing challenges from technological innovation and market competition, such as the rise of battery storage technology, pumped storage power stations need to continuously improve efficiency and reduce costs to maintain their competitiveness in the energy storage market. In the future, as the power system&#x2019;s demands for flexibility and reliability increase, pumped storage power stations will continue to play a key role in the power system. Literature (<xref ref-type="bibr" rid="B29">Sheng and Sun, 2014</xref>) predicts that, combined with digital technology and smart grids, pumped storage power stations will become more intelligent and efficient, providing stronger support for the stability of the grid and the integration of renewable energy. Literature (<xref ref-type="bibr" rid="B7">Ding et al., 2016</xref>) notes that the construction of pumped storage power stations can impact local hydrological conditions and ecosystems, necessitating detailed environmental impact assessments in site selection and design, hence, the construction of pumped storage power stations should not be carried out haphazardly. Maximizing the role of pumped storage power stations and adopting multi-energy joint dispatch based on pumped storage is a viable approach.</p>
<p>Joint dispatch refers to the collaborative work and optimized allocation of different types of energy sources, such as wind, solar, hydro, and thermal power. This concept is widely discussed in literature (<xref ref-type="bibr" rid="B39">Zhou et al., 2016</xref>; <xref ref-type="bibr" rid="B38">Zhou et al., 2017</xref>; <xref ref-type="bibr" rid="B37">Zhang et al., 2018</xref>). Through joint dispatch, it is possible to effectively balance and utilize the advantages and limitations of various energy sources, improve overall energy efficiency, reduce energy costs, and simultaneously decrease environmental pollution. Pumped storage power stations play a key role in joint dispatch systems. Literature (<xref ref-type="bibr" rid="B24">Menglin et al., 2018</xref>) points out that pumped storage can serve as a buffer energy storage facility, helping to balance the unstable output of wind and solar energy, and improve the stability and reliability of the power grid. Additionally, it can store energy during low electricity demand periods for use during peak times, thus enhancing the flexibility and economic efficiency of energy utilization. To achieve optimal joint dispatch, various factors need to be considered, including the prediction of energy production, demand-side response, and market price fluctuations. Literature (<xref ref-type="bibr" rid="B22">Lingamuthu and Mariappan, 2019</xref>; <xref ref-type="bibr" rid="B33">Xu et al., 2019</xref>; <xref ref-type="bibr" rid="B26">Nedaei and Walsh, 2022</xref>) indicates that multi-energy scheduling, pumped-storage power stations serve as energy storage systems to balance supply and demand, thereby enhancing the energy efficiency and stability of the system. Literature (<xref ref-type="bibr" rid="B25">Nedaei et al., 2023</xref>) is a commendable work that provides groundbreaking insights into inverter technology, offering substantial advancements in their design and efficiency through innovative modulation techniques.</p>
<p>In summary, this paper introduces pumped storage power stations and investigates the optimization dispatch problem of complementary systems including hydropower, wind power, solar power, thermal power, and pumped storage, fully exploring their potential for flexible regulation. The joint operation of wind, solar, water, and thermal power based on pumped storage power stations is not only a supplement and improvement to traditional energy systems but also a crucial step towards a cleaner, more efficient, and more sustainable energy future. This study aims to delve into the potential and implementation strategies of this model, providing practical cases and theoretical support for the global energy transition.</p>
</sec>
<sec id="s2">
<title>2 Basic system framework</title>
<p>The joint dispatch framework for the complementary utilization of multiple generation methods such as wind, solar, hydro, thermal, and storage established in this paper is shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. represents a system framework that integrates a variety of renewable and conventional energy sources into an electric power system. Within this system, wind power is depicted by symbols representing four wind turbines, signifying the conversion of wind energy into electrical power. Hydroelectric power is indicated by the illustration of a dam, highlighting the generation of electricity from the kinetic energy of flowing water. Solar photovoltaic power is symbolized by icons of solar panels, denoting the transformation of solar energy into electrical energy through the photovoltaic effect, and these three renewable sources are connected to inverters. Conventional thermal power generation is represented by an image of a factory, typically involving the combustion of fossil fuels to produce electricity. Lastly, the pumped storage plant represents an energy storage method that uses excess electricity to pump water into an elevated reservoir during low demand periods and releases the water flow to generate electricity when demand increases. These two conventional sources are linked to transformers. Collectively, this framework illustrates the unified management and delivery of electric power generated in various ways to meet the load requirements as needed, with an emphasis on the seamless integration of diverse energy generation methods. Within the considered framework, the objective function mainly focuses on the economic dispatch of the system, abandoned power, load shedding, etc., normalizing all comprehensive factors into economic dispatch costs, specifically the operating costs of each unit, the penalty costs for abandoning new energy power, and involuntary load shedding. The constraints include the output limits and ramp rate limits of each unit, as well as the pumping and generation power and capacity constraints of the pumped storage units.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>System framework diagram.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g001.tif"/>
</fig>
<p>In this system: Traditional thermal power generation faces issues such as ramp rate limitations and difficulties in coordinating high operational cost intervals, while the generation mode of hydropower units is more flexible compared to thermal power units, with their output and cost function being nearly linear, allowing them to complement thermal power generation and enhance the utilization rate of the system&#x2019;s generation resources. Pumped storage units, as a type of energy storage resource that can pump or release water resources, can work in conjunction with hydropower units to achieve flexible energy conversion, thereby deeply absorbing renewable energy sources (such as wind turbines and photovoltaics) within the system. Wind turbines and photovoltaics, as renewable energy sources within the system, are utilized as much as possible, in response to new energy policies.</p>
</sec>
<sec id="s3">
<title>3 Mathematical model</title>
<p>This model comprehensively considers the operating costs of thermal power units, hydropower units, pumped storage power stations, operating and maintenance costs of new energy generation units, penalty costs for abandoning new energy power, and involuntary load shedding costs. It establishes a dispatch model aimed at minimizing the total operating costs of generation, including wind power, photovoltaics, thermal power, hydropower, and pumped storage. The objective function is as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi mathvariant="italic">min</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>In Eq. <xref ref-type="disp-formula" rid="e1">1</xref>: where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the total operating cost of the system, <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the optimized dispatch cost of thermal power units, <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the optimized dispatch cost for renewable energy units (wind turbines, photovoltaics), <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the optimized dispatch cost for hydroelectric units, <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the optimized dispatch cost for pumped-storage, <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the penalty cost for curtailment of new energy sources, and <inline-formula id="inf7">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the penalty cost for involuntary load shedding.<list list-type="simple">
<list-item>
<p>(1) The operating costs of thermal power units</p>
</list-item>
</list>
</p>
<p>Flexible control of thermal power units is a crucial prerequisite for ensuring stable operation of power grids with high penetration of renewable energy. Renewable energy output is highly volatile, and to ensure their grid access space, sometimes even starting and stopping operations of thermal power units are required. Therefore, the operating costs of thermal power units primarily consist of coal consumption costs and start-stop costs, expressed by the formula:<disp-formula id="e2">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>In Eqs <xref ref-type="disp-formula" rid="e2">2</xref>&#x2013;<xref ref-type="disp-formula" rid="e4">4</xref>, <inline-formula id="inf8">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf9">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the coal consumption cost and the start-stop cost of the thermal power unit, respectively; <inline-formula id="inf10">
<mml:math id="m14">
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the scheduling time, <inline-formula id="inf11">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the total number of thermal power units; <inline-formula id="inf12">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf13">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf14">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the parameters for the coal consumption cost of electricity generation of the thermal power unit <inline-formula id="inf15">
<mml:math id="m19">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf16">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the start-stop cost of the unit <inline-formula id="inf17">
<mml:math id="m21">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> at the moment <inline-formula id="inf18">
<mml:math id="m22">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf19">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf20">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the power generation and operation status of the thermal power unit <inline-formula id="inf21">
<mml:math id="m25">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> in the time period <inline-formula id="inf22">
<mml:math id="m26">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, with <inline-formula id="inf23">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> representing operation and <inline-formula id="inf24">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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>
</inline-formula> representing stop (<xref ref-type="bibr" rid="B8">Esmaeily et al., 2016</xref>).<list list-type="simple">
<list-item>
<p>(2) The operating costs of hydroelectric units</p>
</list-item>
</list>
</p>
<p>Hydroelectric operating costs include the water consumption for electricity generation and the loss due to unit start-stop. The process of each unit start-stop can also be converted into equivalent water consumption. The formula is as follows:<disp-formula id="e5">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi mathvariant="normal">w</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>In the formula Eq. <xref ref-type="disp-formula" rid="e5">5</xref>: <inline-formula id="inf25">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the generation flow of hydroelectric unit <inline-formula id="inf26">
<mml:math id="m31">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> during period <inline-formula id="inf27">
<mml:math id="m32">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf28">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the equivalent water consumption flow for the start-stop of the hydroelectric unit <inline-formula id="inf29">
<mml:math id="m34">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf30">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the start-stop operation of unit <inline-formula id="inf31">
<mml:math id="m36">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> in period <inline-formula id="inf32">
<mml:math id="m37">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf33">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is a 0-1 variable (<xref ref-type="bibr" rid="B24">Menglin et al., 2018</xref>).<list list-type="simple">
<list-item>
<p>(3) Operating costs of pumped-storage power stations</p>
</list-item>
</list>
<disp-formula id="e6">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>c</mml:mi>
</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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<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>X</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>x</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>In the formula Eq. <xref ref-type="disp-formula" rid="e6">6</xref>: <inline-formula id="inf34">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf35">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represent the unit cost of electricity generation and the unit cost of pumping, respectively; <inline-formula id="inf36">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the total number of pumped-storage units; <inline-formula id="inf37">
<mml:math id="m43">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf38">
<mml:math id="m44">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>x</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the power generation and pumping power consumption of the pumped-storage power station at time <inline-formula id="inf39">
<mml:math id="m45">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, respectively (<xref ref-type="bibr" rid="B10">Garcia-Gonzalez et al., 2008</xref>).<list list-type="simple">
<list-item>
<p>(4) Operation and maintenance costs of new energy generation units</p>
</list-item>
</list>
</p>
<p>Wind and photovoltaic are both new energy generations, but due to the uncontrollable nature of wind speed and solar radiation, some operating costs are incurred during generation. It is assumed that the operation and maintenance costs of various types of units are directly proportional to the amount of electricity generated, the formula is:<disp-formula id="e7">
<mml:math id="m46">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>k</mml:mi>
</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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>In the formula Eq. <xref ref-type="disp-formula" rid="e7">7</xref>: <inline-formula id="inf40">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf41">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are respectively the unit operation and maintenance cost coefficients for wind farms and photovoltaic power stations; <inline-formula id="inf42">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the total number of units in wind farms and photovoltaic power stations, <inline-formula id="inf43">
<mml:math id="m50">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf44">
<mml:math id="m51">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the dispatch output values of wind farms and photovoltaic power stations at moment <inline-formula id="inf45">
<mml:math id="m52">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, respectively (<xref ref-type="bibr" rid="B31">Sun and Harrison, 2019</xref>).<list list-type="simple">
<list-item>
<p>(5) The penalty cost for curtailment of new energy sources</p>
</list-item>
</list>
</p>
<p>Considering the penalty for curtailment of wind and photovoltaic power, the formula is as follows:<disp-formula id="e8">
<mml:math id="m53">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>q</mml:mi>
</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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>In the formula Eq. <xref ref-type="disp-formula" rid="e8">8</xref>: <inline-formula id="inf46">
<mml:math id="m54">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf47">
<mml:math id="m55">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are respectively the penalty cost coefficients for curtailed wind and photovoltaic power; <inline-formula id="inf48">
<mml:math id="m56">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the total number of equipment for curtailment; <inline-formula id="inf49">
<mml:math id="m57">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf50">
<mml:math id="m58">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the amounts of curtailed wind and photovoltaic power at the corresponding moments.<list list-type="simple">
<list-item>
<p>(6) System load shedding penalty costs</p>
</list-item>
</list>
</p>
<p>To ensure the stability and quality of the power grid, when the actual output of the system deviates from the load demand, a load shedding penalty cost calculation function must be introduced, with penalties applied according to the size of the deviation, as follows:<disp-formula id="e9">
<mml:math id="m59">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>u</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<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:mi>T</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mi>&#x3b2;</mml:mi>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>In the formula Eq. <xref ref-type="disp-formula" rid="e9">9</xref>: <inline-formula id="inf51">
<mml:math id="m60">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the load shedding cost penalty factor, <inline-formula id="inf52">
<mml:math id="m61">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the power of load shedding at moment <inline-formula id="inf53">
<mml:math id="m62">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
<sec id="s4">
<title>4 Constraints</title>
<p>
<list list-type="simple">
<list-item>
<p>(1) Power balance constraints</p>
</list-item>
</list>
<disp-formula id="e10">
<mml:math id="m63">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>c</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>x</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<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:msubsup>
<mml:mi>P</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>u</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>q</mml:mi>
</mml:msub>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>In Eq. <xref ref-type="disp-formula" rid="e10">10</xref>: <inline-formula id="inf54">
<mml:math id="m64">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the power of thermal power units; <inline-formula id="inf55">
<mml:math id="m65">
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the power of hydroelectric units; <inline-formula id="inf56">
<mml:math id="m66">
<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> is the load of the power grid at moment <inline-formula id="inf57">
<mml:math id="m67">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.<list list-type="simple">
<list-item>
<p>(2) Power flow constraints</p>
</list-item>
</list>
<disp-formula id="e11">
<mml:math id="m68">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</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>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m69">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>The flow model represented by Eqs <xref ref-type="disp-formula" rid="e11">11</xref>, <xref ref-type="disp-formula" rid="e12">12</xref> is the Direct Current (DC) flow model. The rationale for establishing this model includes: 1) Given that the model scenario pertains to a regional system, where individual units are relatively close to each other, it is feasible to approximate multiple nodes within the system as a single node. 2) The model primarily focuses on the supply relationship between generation and load, under which consideration of active power takes precedence, thereby justifying the omission of reactive power. 3) As the distances between nodes within the considered scenario are relatively short, the voltage levels across these nodes are expected to be similar, allowing for the exclusion of voltage variations. 4) From a computational complexity perspective, disregarding phase angles and frequency variations within the flow model, and focusing solely on power flow, simplifies the computational complexity. Based on these considerations, the model under study can be effectively equated to a Direct Current flow model.</p>
<p>It should be noted that the DC flow model has significant advantages in certain situations, but it also has limitations, such as its inability to accurately describe phase angles and frequency variations found in Alternating Current (AC) flow models (<xref ref-type="bibr" rid="B20">Kile et al., 2015</xref>). However, within the context of this study, the DC flow model is adopted for calculation and analysis, given its applicability.</p>
<p>In the formula <inline-formula id="inf58">
<mml:math id="m70">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the reactance of the branch; <inline-formula id="inf59">
<mml:math id="m71">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the capacity between branches <inline-formula id="inf60">
<mml:math id="m72">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf61">
<mml:math id="m73">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the power transmission between nodes <inline-formula id="inf62">
<mml:math id="m74">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf63">
<mml:math id="m75">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf64">
<mml:math id="m76">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are respectively the voltage angles of <inline-formula id="inf65">
<mml:math id="m77">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> at moment <inline-formula id="inf66">
<mml:math id="m78">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.<list list-type="simple">
<list-item>
<p>(3) Thermal power unit operational constraints</p>
</list-item>
</list>
</p>
<p>Operational constraints of thermal power units mainly include: the maximum and minimum output of the units, ramping constraints of the units. The formula is as follows:<disp-formula id="e13">
<mml:math id="m79">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m80">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
<p>In the formula: <inline-formula id="inf67">
<mml:math id="m81">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf68">
<mml:math id="m82">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are respectively the minimum and maximum outputs of thermal power unit <inline-formula id="inf69">
<mml:math id="m83">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf70">
<mml:math id="m84">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</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="inf71">
<mml:math id="m85">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</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:math>
</inline-formula> are respectively the power generation of thermal power unit <inline-formula id="inf72">
<mml:math id="m86">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> at moments <inline-formula id="inf73">
<mml:math id="m87">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf74">
<mml:math id="m88">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf75">
<mml:math id="m89">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf76">
<mml:math id="m90">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are respectively the maximum and minimum ramp rates of unit <inline-formula id="inf77">
<mml:math id="m91">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. Meanwhile, the minimum operational and shutdown constraints of the unit are as follows:<disp-formula id="e15">
<mml:math id="m92">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x2200;</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m93">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x2200;</mml:mo>
<mml:mi>H</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>
</p>
<p>In Eqs <xref ref-type="disp-formula" rid="e15">15</xref>, <xref ref-type="disp-formula" rid="e16">16</xref>: wherein, <inline-formula id="inf78">
<mml:math id="m94">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the operating status of unit <inline-formula id="inf79">
<mml:math id="m95">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> at moment <inline-formula id="inf80">
<mml:math id="m96">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf81">
<mml:math id="m97">
<mml:mrow>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> refers to the hourly time period; <inline-formula id="inf82">
<mml:math id="m98">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
<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:math>
</inline-formula> is the operating status of unit <inline-formula id="inf83">
<mml:math id="m99">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> at moment <inline-formula id="inf84">
<mml:math id="m100">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. <inline-formula id="inf85">
<mml:math id="m101">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf86">
<mml:math id="m102">
<mml:mrow>
<mml:msub>
<mml:mi>T</mml:mi>
<mml:mrow>
<mml:mi>o</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> respectively represent the minimum operating and shutdown time of the unit, <inline-formula id="inf87">
<mml:math id="m103">
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the total dispatch time.<list list-type="simple">
<list-item>
<p>(4) Operational constraints of wind and photovoltaic power generation</p>
</list-item>
</list>
<disp-formula id="e17">
<mml:math id="m104">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(17)</label>
</disp-formula>
<disp-formula id="e18">
<mml:math id="m105">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>p</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>4</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>In Eqs <xref ref-type="disp-formula" rid="e17">17</xref>, <xref ref-type="disp-formula" rid="e18">18</xref>, <inline-formula id="inf88">
<mml:math id="m106">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf89">
<mml:math id="m107">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</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> represent the output power of wind and photovoltaic at time <inline-formula id="inf90">
<mml:math id="m108">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, respectively. <inline-formula id="inf91">
<mml:math id="m109">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf92">
<mml:math id="m110">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> denote the output power constraints of wind and photovoltaic at time <inline-formula id="inf93">
<mml:math id="m111">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, respectively.<list list-type="simple">
<list-item>
<p>(5) Hydroelectric operation constraints</p>
</list-item>
</list>
</p>
<p>Hydroelectric operation constraints include output constraints, hydroelectric energy conversion constraints, and daily flow constraints, as shown in the following formulas:<disp-formula id="e19">
<mml:math id="m112">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>h</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>s</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(19)</label>
</disp-formula>
<disp-formula id="e20">
<mml:math id="m113">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>A</mml:mi>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>h</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="m114">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(21)</label>
</disp-formula>
</p>
<p>In Eqs <xref ref-type="disp-formula" rid="e19">19</xref>&#x2013;<xref ref-type="disp-formula" rid="e21">21</xref>, <inline-formula id="inf94">
<mml:math id="m115">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf95">
<mml:math id="m116">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represent the maximum and minimum output of hydroelectric units <inline-formula id="inf96">
<mml:math id="m117">
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. <inline-formula id="inf97">
<mml:math id="m118">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>h</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the output of hydroelectric unit <inline-formula id="inf98">
<mml:math id="m119">
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> during time period <inline-formula id="inf99">
<mml:math id="m120">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. <inline-formula id="inf100">
<mml:math id="m121">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the hydroelectric conversion coefficient in m<sup>3</sup>/kWh, and <inline-formula id="inf101">
<mml:math id="m122">
<mml:mrow>
<mml:msub>
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the upstream head of hydroelectric unit during time period <inline-formula id="inf102">
<mml:math id="m123">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> in meters. <inline-formula id="inf103">
<mml:math id="m124">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf104">
<mml:math id="m125">
<mml:mrow>
<mml:msubsup>
<mml:mi>Q</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represent the upper and lower limits of daily flow for hydroelectric units.<list list-type="simple">
<list-item>
<p>(6) Pumped Storage Constraints</p>
</list-item>
</list>
</p>
<p>For simplicity, define the output power of pumped storage power plants as <inline-formula id="inf105">
<mml:math id="m126">
<mml:mrow>
<mml:msub>
<mml:mi>P</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>. If <inline-formula id="inf106">
<mml:math id="m127">
<mml:mrow>
<mml:msub>
<mml:mi>P</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 greater than or equal to zero, the system is in a discharging and generating state; otherwise, it is in a pumping and energy storage state. Pumped storage power constraints:<disp-formula id="e22">
<mml:math id="m128">
<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>min</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</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>P</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(22)</label>
</disp-formula>
</p>
<p>In the Eq. <xref ref-type="disp-formula" rid="e22">22</xref>, <inline-formula id="inf107">
<mml:math id="m129">
<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 mathvariant="italic">min</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf108">
<mml:math id="m130">
<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 mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represent the minimum and maximum pumping power of pumped storage power plants.</p>
<p>Hydroturbine power generation constraints:<disp-formula id="e23">
<mml:math id="m131">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>min</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>max</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(23)</label>
</disp-formula>
</p>
<p>In the Eq. <xref ref-type="disp-formula" rid="e23">23</xref>, <inline-formula id="inf109">
<mml:math id="m132">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the energy stored in the pumped storage power station at time <inline-formula id="inf110">
<mml:math id="m133">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf111">
<mml:math id="m134">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the hydroturbine efficiency; <inline-formula id="inf112">
<mml:math id="m135">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf113">
<mml:math id="m136">
<mml:mrow>
<mml:msubsup>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> are the minimum and maximum hydroturbine power generation capacities, respectively.</p>
<p>Operational constraints in Eq. <xref ref-type="disp-formula" rid="e24">24</xref>:<disp-formula id="e24">
<mml:math id="m137">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</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>(24)</label>
</disp-formula>
<inline-formula id="inf114">
<mml:math id="m138">
<mml:mrow>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the hydroturbine power generation at time <inline-formula id="inf115">
<mml:math id="m139">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>. Pumped storage power plants operate in either discharging and generating or pumping and energy storage modes, and these modes do not occur simultaneously.</p>
<p>Reservoir energy balance constraints in Eq. <xref ref-type="disp-formula" rid="e25">25</xref>:<disp-formula id="e25">
<mml:math id="m140">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>h</mml:mi>
<mml:mi>g</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(25)</label>
</disp-formula>
</p>
<p>In the equation, <inline-formula id="inf116">
<mml:math id="m141">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the energy stored in the pumped storage power station at time <inline-formula id="inf117">
<mml:math id="m142">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf118">
<mml:math id="m143">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the pump efficiency, and <inline-formula id="inf119">
<mml:math id="m144">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b7;</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the hydroturbine efficiency.</p>
<p>Reservoir capacity constraints in Eq. <xref ref-type="disp-formula" rid="e26">26</xref>:<disp-formula id="e26">
<mml:math id="m145">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>E</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(26)</label>
</disp-formula>
</p>
<p>In the equation, <inline-formula id="inf120">
<mml:math id="m146">
<mml:mrow>
<mml:msubsup>
<mml:mi>E</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> represents the maximum energy stored in the pumped storage power station at time <inline-formula id="inf121">
<mml:math id="m147">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
<sec id="s5">
<title>5 Case analysis</title>
<sec id="s5-1">
<title>5.1 Description of the algorithm</title>
<p>This paper focuses on the research of a county-level power grid in southern China, where the primary power sources include thermal power plants, hydroelectric power plants, and pumped storage power plants. Among them, there are 6 thermal power units in the thermal power plants, 4 hydroelectric units in the hydroelectric power plants, and 2 pumping units in the pumped storage power station. The hydroelectric conversion coefficient is assumed to be 90%. The upstream water head of the hydropower unit is 50&#xa0;m. In this region, the thermal power plants are the primary power source for supplying the load, while the hydroelectric power plants and pumped storage power plants participate in energy distribution as auxiliary forms to regulate it. The test cases were conducted on a laptop with AMD Ryzen 7 7735H with Radeon Graphics CPU, 3.2 GHz and 32 GB RAM. The MATLAB software with YALMIP toolbox and CPLEX solver were used to solve the optimization problems.</p>
<p>Detailed information on the parameters of each power plant is provided in the appendix. Wind and photovoltaic power forecasts as well as load forecasts within the regional grid are shown in <xref ref-type="fig" rid="F2">Figures 2</xref>, <xref ref-type="fig" rid="F3">3</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Predicted output of wind and photovoltaic power within the district-level grid.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g002.tif"/>
</fig>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Predicted load value within the district-level grid.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g003.tif"/>
</fig>
</sec>
<sec id="s5-2">
<title>5.2 Analysis of results</title>
<p>To validate the advantages of the proposed model, results are compared in the following four modes for further explanation. Mode 1: Initial mode, where the primary power generation is solely handled by thermal power plants, and there are no reservoirs or pumped storage in the system. Mode 2: Involves hydroelectric plants in scheduling optimization, working in coordination with thermal power plants for optimized operations, with reservoirs but no pumped storage. Mode 3: No hydroelectric plants participate in optimization, but there is auxiliary regulation by pumped storage plants, with no reservoirs but pumped storage available. Mode 4: The model established in this paper, which includes both reservoirs and pumped storage. Economic operating costs and curtailment rates under the four modes are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Economic operating costs and power abandonment in four models.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Model</th>
<th align="center">Total Running costs</th>
<th align="center">Unit operating costs</th>
<th align="center">Abandonment penalties/load-shedding costs</th>
<th align="center">Power abandonment rate</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Model 1</td>
<td align="center">11,184.93</td>
<td align="center">9,528.88</td>
<td align="center">1,656.05</td>
<td align="center">0.126</td>
</tr>
<tr>
<td align="center">Model 2</td>
<td align="center">11,076.28</td>
<td align="center">9,702.47</td>
<td align="center">1,373.81</td>
<td align="center">0.126</td>
</tr>
<tr>
<td align="center">Model 3</td>
<td align="center">10,572.54</td>
<td align="center">9,508.69</td>
<td align="center">1,063.84</td>
<td align="center">0.087</td>
</tr>
<tr>
<td align="center">Model 4</td>
<td align="center">9,923.078</td>
<td align="center">9,667.24</td>
<td align="center">255.836</td>
<td align="center">0.023</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>From <xref ref-type="table" rid="T1">Table 1</xref>, it is evident that in Mode 1, the total operating cost of the system is the highest. This is because only thermal power units generate electricity in the system, and there is no pumped storage and hydroelectric units for auxiliary regulation, leading to higher curtailment penalties/load shedding costs. In Mode 2, the participation of hydroelectric units significantly reduces curtailment penalties/load shedding costs, although the unit operating costs are relatively higher, resulting in a significant reduction in total operating costs. In Mode 3, the introduction of pumped storage reduces both unit operating costs and curtailment penalties/load shedding costs, emphasizing the importance of pumped storage. In Mode 4, although unit operating costs are slightly higher compared to Mode 3, other costs and curtailment rates are reduced, demonstrating that the coordination between pumped storage and hydroelectric units can optimize system scheduling performance.</p>
<p>Curtailment, load shedding, and new energy integration under the four modes are as follows in <xref ref-type="fig" rid="F4">Figures 4</xref>&#x2013;<xref ref-type="fig" rid="F15">15</xref>:</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Conditions of energy curtailment and load shedding in Mode 1.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Wind power consumption in Mode 1.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Photovoltaic power consumption in Mode 1.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Energy curtailment and load shedding in Mode 2.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g007.tif"/>
</fig>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Wind power consumption in Mode 2.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g008.tif"/>
</fig>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Photovoltaic power consumption in Mode 2.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g009.tif"/>
</fig>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Energy curtailment and load shedding in Mode 3.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Wind power consumption in Mode 3.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g011.tif"/>
</fig>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Photovoltaic power consumption in Mode 3.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g012.tif"/>
</fig>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Energy curtailment and load shedding in Mode 4.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g013.tif"/>
</fig>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>Wind power consumption in Mode 4.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g014.tif"/>
</fig>
<fig id="F15" position="float">
<label>FIGURE 15</label>
<caption>
<p>Photovoltaic power consumption in Mode 4.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g015.tif"/>
</fig>
<p>From the scenarios of renewable energy consumption, curtailment, and load shedding described above, it is evident that: In Modes 2 and 4, the system does not experience load shedding, as the output of the hydroelectric units can prevent load shedding incidents. In Mode 3, there is a small amount of load shedding due to the limited capacity or regulatory ability of pumped storage, which prevents the transfer of more energy to the load shedding gap, resulting in system load shedding. In Mode 4, the system&#x2019;s energy curtailment is minimized, and there is no load shedding, demonstrating the superiority of the model presented in this paper.</p>
<p>The operational results of each unit under Mode 4 are as follows:</p>
<p>
<xref ref-type="fig" rid="F16">Figure 16</xref> shows the power output of thermal power units, hydroelectric units, and pumped-storage power units over a 24-h period. The chart reveals that the power output of thermal power units (red bars) and hydroelectric units (black bars) is relatively stable, whereas the pumped-storage power units (purple line) exhibit significant negative and positive power outputs, reflecting their energy storage and release characteristics.</p>
<fig id="F16" position="float">
<label>FIGURE 16</label>
<caption>
<p>Operational results of each unit in Mode 4.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g016.tif"/>
</fig>
<p>The analysis of the joint operation of these three types of equipment demonstrates the critical role of pumped-storage in peak shaving and balancing the electricity grid load. During the night, when electricity demand is low and the production capacity of thermal and hydroelectric units might be excessive, pumped-storage units operate at negative power to absorb surplus electricity, thus avoiding resource waste. During peak electricity usage times, the pumped-storage units respond rapidly to increased load through positive power operation, supporting grid stability and reducing reliance on thermal power units.</p>
<p>The operation strategy of the pumped-storage power plant is as follows:</p>
<p>
<xref ref-type="fig" rid="F17">Figure 17</xref> shows the daily operational strategy of pumped storage. The graph shows the output variation of the pumped-storage units over a 24-h period: during low demand periods (such as nighttime), they exhibit negative output, meaning the units consume electricity to pump water to the upper reservoir; during high demand periods (usually daytime and evening), they switch to positive output, releasing water to generate electricity to meet high power demands.</p>
<fig id="F17" position="float">
<label>FIGURE 17</label>
<caption>
<p>Operational strategy of the pumped-storage power plant in Mode 4.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g017.tif"/>
</fig>
<p>The load adjustment before and after is as follows:</p>
<p>
<xref ref-type="fig" rid="F18">Figure 18</xref> shows the load data comparison before and after pumped storage participates in grid load regulation. By comparing the two curves - before regulation (unregulated) and after regulation (regulated), we can analyze the role of pumped storage in electricity system load management. The vertical axis represents the load (in KW), and the horizontal axis represents the time of day (0&#x2013;24&#xa0;h). The graph shows that in the unregulated scenario, the load curve exhibits significant peak-valley differences, meaning the power system will face high loads during certain periods (such as afternoon and evening) and lower loads at other times (such as late night and early morning).</p>
<fig id="F18" position="float">
<label>FIGURE 18</label>
<caption>
<p>Before and after load adjustment in Mode 4.</p>
</caption>
<graphic xlink:href="fenrg-12-1373588-g018.tif"/>
</fig>
<p>After the pumped-storage system is put into operation, the adjusted load curve becomes smoother, especially during the periods of original peak load. This indicates that pumped storage effectively reduces peak loads and lowers the operational pressure on the grid by storing energy (pumping water to the upper reservoir) during low load periods and releasing energy (generating electricity) during high load periods.</p>
<p>The operational strategy of pumped storage reflects its multiple advantages. Firstly, its peak-shaving capability, as shown in the positive output phase in <xref ref-type="fig" rid="F17">Figure 17</xref>, demonstrates the pumped storage system releasing stored water energy to generate electricity during high demand periods, helping the grid handle peak loads. This peak-shaving capability is a key component of the flexibility of the electricity system, particularly in the context of the increasing share of renewable energies like wind and solar power. Secondly, its energy storage capability is revealed during the negative output phase, highlighting the ability of pumped storage to store energy during low demand periods. In this way, it can convert surplus electrical energy into potential energy, increasing the storage capacity of the electricity system. Finally, its high efficiency and environmental friendliness. Pumped storage is one of the most mature large-scale electricity storage technologies, with a cycle efficiency of 70%&#x2013;80%. As a clean storage technology, pumped storage does not produce greenhouse gases or other pollutants during operation, helping to reduce reliance on fossil fuels and lower the carbon footprint of the electricity system. Through <xref ref-type="fig" rid="F18">Figure 18</xref>, we can visually observe how pumped storage adjusts its operational mode according to the daily cyclical changes in electricity demand, demonstrating its indispensable role and value in the electricity market. Furthermore, the involvement of pumped storage also helps to reduce the electricity system&#x2019;s dependence on traditional peak-shaving power plants (such as gas-fired power plants), thereby lowering operational costs and environmental impacts.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>This paper considers the coordinated dispatch of flexible resources such as pumped storage and hydropower units in traditional power systems and proposes a joint dispatch model for the complementary utilization of various generation methods like wind, solar, hydro, thermal, and storage. Conclusions are as follows:<list list-type="simple">
<list-item>
<p>(1) The joint dispatch model established in this paper for the complementary utilization of wind, solar, hydro, thermal, and storage generation methods has generally reduced the economic dispatch cost of the system, lowered the level of load shedding and the rate of abandoned power, and increased the absorption of wind and solar power;</p>
</list-item>
<list-item>
<p>(2) In the proposed model, hydropower units can buffer issues such as the high-cost operating intervals and startup/shutdown difficulties of traditional thermal power units, optimizing power generation in conjunction with thermal power units, thereby enhancing the economic efficiency of the system&#x2019;s operation;</p>
</list-item>
<list-item>
<p>(3) In the proposed model, pumped storage units can fully mobilize hydraulic resources, dispatching energy within the generation intervals where the original units bear high output costs or where there is surplus wind and solar power, thereby improving system load characteristics and reducing the burden of output on units within the system.</p>
</list-item>
</list>
</p>
<p>Although joint dispatch offers many advantages, it also faces numerous challenges in practical application, including technical complexity, cost investment, policies, and market mechanisms. Future research needs to focus more on technological innovation, cost-benefit analysis, and policy support mechanisms to promote the widespread application and development of this strategy.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<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="s8">
<title>Author contributions</title>
<p>JJ: Conceptualization, Data curation, Resources, Writing&#x2013;original draft. GZ: Funding acquisition, Project administration, Supervision, Writing&#x2013;review and editing. XZ: Conceptualization, Investigation, Writing&#x2013;review and editing. ZS: Methodology, Software, Validation, Visualization, Writing&#x2013;review and editing. MZ: Writing&#x2013;review and editing, Data curation, Methodology, Resources. XL: Writing&#x2013;review and editing, Formal Analysis, Validation.</p>
</sec>
<sec sec-type="funding-information" id="s9">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by Jiangsu Blue Project; Taizhou Science and Technology Support Plan (Social Development) Project (SSF20230014); Research Team Construction Project of Taizhou College, Nanjing Normal University: Machine Vision Deep Learning Research Team.</p>
</sec>
<ack>
<p>The authors would like to express their gratitude to all those who helped them during the writing of this paper.</p>
</ack>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenrg.2024.1373588/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenrg.2024.1373588/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet2.pdf" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arce</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ohishi</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Soares</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Optimal dispatch of generating units of the itaip &#xfa; hydroelectric plant</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>17</volume>, <fpage>154</fpage>&#x2013;<lpage>158</lpage>. <pub-id pub-id-type="doi">10.1109/59.982207</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Azizipanah-Abarghooee</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Niknam</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Roosta</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Malekpour</surname>
<given-names>A. R.</given-names>
</name>
<name>
<surname>Zare</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Probabilistic multiobjective wind-thermal economic emission dispatch based on point estimated method</article-title>. <source>Energy</source> <volume>37</volume> (<issue>1</issue>), <fpage>322</fpage>&#x2013;<lpage>335</lpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2011.11.023</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bello</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Suarez</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Drenkhan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Vega-J&#xe1;come</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Hydrological impacts of dam regulation for hydropower production: the case of Lake Sibinacocha, Southern Peru</article-title>. <source>J. Hydrology Regional Stud.</source> <volume>46</volume>, <fpage>101319</fpage>. <pub-id pub-id-type="doi">10.1016/j.ejrh.2023.101319</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Canizes</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Soares</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Vale</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Khodr</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Hybrid fuzzy Monte Carlo technique for reliability assessment in transmission power systems</article-title>. <source>Energy</source> <volume>45</volume> (<issue>1</issue>), <fpage>1007</fpage>&#x2013;<lpage>1017</lpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2012.06.049</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cuiping</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Junwu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Junhui</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Gangui</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Limin</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zehui</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Characteristic analysis of photovoltaic power generation and wind-photovoltaic joint power generation system</article-title>. <source>Power Syst. Clean Energy</source> <volume>33</volume> (<issue>01</issue>), <fpage>95</fpage>&#x2013;<lpage>102</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1674-3814.2017.01.015</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deane</surname>
<given-names>J. P.</given-names>
</name>
<name>
<surname>&#xd3; Gallach&#xf3;ir</surname>
<given-names>B. P.</given-names>
</name>
<name>
<surname>McKeogh</surname>
<given-names>E. J.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Techno-economic review of existing and new pumped hydro energy storage plant</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>14</volume>, <fpage>1293</fpage>&#x2013;<lpage>1302</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2009.11.015</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ding</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bie</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Robust Co-optimization to energy and ancillary service joint dispatch considering wind power uncertainties in real-time electricity markets</article-title>. <source>IEEE Trans. Sustain. Energy</source> <volume>7</volume>, <fpage>1547</fpage>&#x2013;<lpage>1557</lpage>. <pub-id pub-id-type="doi">10.1109/TSTE.2016.2561967</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Esmaeily</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Raeisi</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ahmadi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ahmadi</surname>
<given-names>M. R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>A note on short-term hydro-thermal scheduling</article-title>. <source>Energy Convers. Manag.</source> <volume>126</volume>, <fpage>1178</fpage>&#x2013;<lpage>1186</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2016.08.083</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Mai</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Nonlinear model predictive control for pumped storage plants based on online sequential extreme learning machine with forgetting factor</article-title>. <source>Complexity</source> <volume>2021</volume>, <fpage>1</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1155/2021/5692621</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garcia-Gonzalez</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>De</surname>
<given-names>l M.</given-names>
</name>
<name>
<surname>Santos</surname>
<given-names>L. M.</given-names>
</name>
<name>
<surname>Gonzalez</surname>
<given-names>A. M.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Stochastic joint optimization of wind generation and pumped-storage units in an electricity market</article-title>. <source>IEEE Trans. Power Syst.</source> <volume>23</volume> (<issue>2</issue>), <fpage>460</fpage>&#x2013;<lpage>468</lpage>. <pub-id pub-id-type="doi">10.1109/TPWRS.2008.919430</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gromyko</surname>
<given-names>P. V.</given-names>
</name>
<name>
<surname>Fedin</surname>
<given-names>K. V.</given-names>
</name>
<name>
<surname>Seleznev</surname>
<given-names>V. S.</given-names>
</name>
<name>
<surname>Kolesnikov</surname>
<given-names>Y. I.</given-names>
</name>
<name>
<surname>Ngomayezwe</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Monitoring of the hydraulic units operation of the Sayano&#x2010;Shushenskaya hydroelectric power plant using remote seismic observations</article-title>. <source>Earthq. Eng. Struct. Dyn.</source> <volume>52</volume>, <fpage>335</fpage>&#x2013;<lpage>349</lpage>. <pub-id pub-id-type="doi">10.1002/eqe.3762</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hinojosa</surname>
<given-names>V. H.</given-names>
</name>
<name>
<surname>Velasquez</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Improving the mathematical formulation of security-constrained generation capacity expansion planning using power transmission distribution factors and line outage distribution factors</article-title>. <source>Electr. Power Syst. Res.</source> <volume>140</volume>, <fpage>391</fpage>&#x2013;<lpage>400</lpage>. <pub-id pub-id-type="doi">10.1016/j.epsr.2016.06.002</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hong-Mei</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Jian-Yu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Qiu-Lan</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Manual complementariness optimization of electrical source involving intermittent power</article-title>. <source>Dianli Xit. Baohu yu Kongzhi/Power Syst. Prot. Control</source> <volume>41</volume> (<issue>9</issue>), <fpage>34</fpage>&#x2013;<lpage>39</lpage>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hou</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>F.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Study of the carbon neutral path in China: a literature review</article-title>. <source>Chin. J. Urban Environ. Stud.</source> <volume>11</volume>. <pub-id pub-id-type="doi">10.1142/S2345748123500082</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>A joint daily operational optimization model for wind power and pumped-storage plant</article-title>. <source>Automation Electr. Power Syst.</source> <volume>48</volume> (<issue>1</issue>) <fpage>571</fpage>&#x2013;<lpage>578</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-1026.2012.02.007</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Active power fluctuations and control in hydraulic disturbance of variable speed pumped storage plants</article-title>. <source>J. Energy Storage</source> <volume>60</volume>, <fpage>106666</fpage>. <pub-id pub-id-type="doi">10.1016/j.est.2023.106666</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jin</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Hong</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Gholinia</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>The evaluation distribution of runoff value on hydroelectric potential change-based RCPs scenarios and soft-computing: a case study</article-title>. <source>J. Electr. Eng. Technol.</source> <volume>18</volume>, <fpage>2551</fpage>&#x2013;<lpage>2565</lpage>. <pub-id pub-id-type="doi">10.1007/s42835-022-01343-5</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kanakasabapathy</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Swarup</surname>
<given-names>R. S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Bidding strategy for pumped-storage plant in pool-based electricity market</article-title>. <source>Energy Convers. Manag.</source> <volume>51</volume> (<issue>3</issue>), <fpage>572</fpage>&#x2013;<lpage>579</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2009.11.001</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Katsuhiro</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Juichiro</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Shiro</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Wada</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Imano</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Development of pump-turbine for seawater pumped storage power plant</article-title>. <source>Am. Soc. Civ. Eng.</source> <volume>47</volume> (<issue>5</issue>), <fpage>199</fpage>&#x2013;<lpage>202</lpage>. <pub-id pub-id-type="doi">10.1061/40440(1999)67</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Kile</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Uhlen</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Warland</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Kjolle</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A comparison of AC and DC power flow models for contingency and reliability analysis</article-title>. <conf-name>2014 Power Systems Computation Conference</conf-name>, <conf-loc>Wroclaw, Poland</conf-loc>, <conf-date>August 2014</conf-date>. <publisher-name>IEEE</publisher-name>. <pub-id pub-id-type="doi">10.1109/PSCC.2014.7038459</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lamb</surname>
<given-names>E. B.</given-names>
</name>
<name>
<surname>Devries</surname>
<given-names>D. R.</given-names>
</name>
<name>
<surname>Wright</surname>
<given-names>R. A.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Peaking hydropower and fish assemblages: an example from the Tallapoosa River, AL</article-title>. <source>River Res. Appl.</source> <volume>39</volume>, <fpage>1001</fpage>&#x2013;<lpage>1011</lpage>. <pub-id pub-id-type="doi">10.1002/rra.4127</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lingamuthu</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Mariappan</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Power flow control of grid connected hybrid renewable energy system using hybrid controller with pumped storage</article-title>. <source>Int. J. Hydrogen Energy</source> <volume>44</volume> (<issue>7</issue>), <fpage>3790</fpage>&#x2013;<lpage>3802</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijhydene.2018.12.092</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Short-term operation of cascade hydropower system sharing flexibility via high voltage direct current lines for multiple grids peak shaving</article-title>. <source>Renew. energy</source> <volume>213</volume>, <fpage>11</fpage>&#x2013;<lpage>29</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2023.05.095</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Menglin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xiaomeng</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Jiakun</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zuo</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>A systematic approach for the joint dispatch of energy and reserve incorporating demand response</article-title>. <source>Appl. Energy</source> <volume>230</volume>, <fpage>1279</fpage>&#x2013;<lpage>1291</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2018.09.044</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nedaei</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Faccio</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Walsh</surname>
<given-names>P. R.</given-names>
</name>
<name>
<surname>Rasul</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>Alain Bordas</surname>
<given-names>S. P.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Novel design, implementation, and performance optimization of inverters by considering the effect of modulation</article-title>. <source>Energy &#x26; Environ.</source>, <fpage>0958305X2311646</fpage>. <comment>(Early Access)</comment>. <pub-id pub-id-type="doi">10.1177/0958305x231164688</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nedaei</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Walsh</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Technical performance evaluation and optimization of a run-of-river hydropower facility</article-title>. <source>Renew. Energy</source> <volume>182</volume>, <fpage>343</fpage>&#x2013;<lpage>362</lpage>. <pub-id pub-id-type="doi">10.1016/j.renene.2021.10.021</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nnamchi</surname>
<given-names>S. N.</given-names>
</name>
<name>
<surname>Natukunda</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wanambwa</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Musiime</surname>
<given-names>E. B.</given-names>
</name>
<name>
<surname>Tukamuhebwa</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Wanazusi</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Effects of wind speed and tropospheric height on solar power generation: energy exploration above ground level</article-title>. <source>Energy Rep.</source> <volume>9</volume>, <fpage>5166</fpage>&#x2013;<lpage>5182</lpage>. <pub-id pub-id-type="doi">10.1016/j.egyr.2023.04.269</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rahman</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Farrok</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Haque</surname>
<given-names>M. M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Environmental impact of renewable energy source based electrical power plants: solar, wind, hydroelectric, biomass, geothermal, tidal, ocean, and osmotic</article-title>. <source>Renew. Sustain. energy Rev.</source> <volume>161</volume>, <fpage>112279</fpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2022.112279</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sheng</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>An economic dispatching strategy of peak load shifting by wind farm and pumped storage plant</article-title>. <source>Power Syst. Technol.</source> <volume>38</volume> (<issue>09</issue>), <fpage>2484</fpage>&#x2013;<lpage>2489</lpage>. <pub-id pub-id-type="doi">10.13335/j.1000-3673.pst.2014.09.027</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Silva</surname>
<given-names>A. R.</given-names>
</name>
<name>
<surname>Pimenta</surname>
<given-names>F. M.</given-names>
</name>
<name>
<surname>Assireu</surname>
<given-names>A. T.</given-names>
</name>
<name>
<surname>Spyrides</surname>
<given-names>M. H. C.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Complementarity of Brazils hydro and offshore wind power</article-title>. <source>Renew. Sustain. Energy Rev.</source> <volume>56</volume>, <fpage>413</fpage>&#x2013;<lpage>427</lpage>. <pub-id pub-id-type="doi">10.1016/j.rser.2015.11.045</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Harrison</surname>
<given-names>G. P.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Wind-solar complementarity and effective use of distribution network capacity</article-title>. <source>Appl. Energy</source> <volume>247</volume>, <fpage>89</fpage>&#x2013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.04.042</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Modeling and coordinated control for active power regulation of pumped storage&#x2010;battery integrated system under small&#x2010;disturbances</article-title>. <source>Energy Sci. Eng.</source> <volume>11</volume> (<issue>5</issue>), <fpage>1601</fpage>&#x2013;<lpage>1618</lpage>. <pub-id pub-id-type="doi">10.1002/ese3.1406</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Venkateshkumar</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Xing</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Modeling a pumped storage hydropower integrated to a hybrid power system with solar-wind power and its stability analysis</article-title>. <source>Appl. Energy</source> <volume>248</volume>, <fpage>446</fpage>&#x2013;<lpage>462</lpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2019.04.125</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Hongwei</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Jikai</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Hongpeng</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Optimal scheduling of thermoelectric interconnection virtual power plant considerign incentive demand response</article-title>. <source>Acta Energiae Solaris Sin.</source> <volume>44</volume> (<issue>4</issue>), <fpage>349</fpage>&#x2013;<lpage>358</lpage>. <pub-id pub-id-type="doi">10.19912/j.0254-0096.tynxb.2021-1524</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ye</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Real-time model predictive control study of run-of-river hydropower plants with data-driven and physics-based coupled model</article-title>. <source>J. Hydrology</source> <volume>617</volume>, <fpage>128942</fpage>. <pub-id pub-id-type="doi">10.1016/j.jhydrol.2022.128942</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zanoli</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Pepe</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Astolfi</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Luzi</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Reservoir advanced process control for hydroelectric power production</article-title>. <source>Processes</source> <volume>11</volume>, <fpage>300</fpage>. <pub-id pub-id-type="doi">10.3390/pr11020300</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Joint dispatch operation model of microgrid and distribution network considering user response to electricity price</article-title>. <source>Electr. Power Autom. Equip.</source> <volume>38</volume> (<issue>05</issue>), <fpage>177</fpage>&#x2013;<lpage>183</lpage>. <pub-id pub-id-type="doi">10.16081/j.issn.1006-6047.2018.05.026</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Multi-area generation-reserve joint dispatch approach considering wind power cross-regional accommodation</article-title>. <source>CSEE J. Power Energy Syst.</source> <volume>3</volume>, <fpage>74</fpage>&#x2013;<lpage>83</lpage>. <pub-id pub-id-type="doi">10.17775/CSEEJPES.2017.0010</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Chao</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>SC-PSO algorithm and IES-CCHP regional joint dispatch with TOU price</article-title>. <source>Electr. Power Autom. Equip.</source> <volume>36</volume> (<issue>12</issue>), <fpage>11</fpage>&#x2013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.16081/j.issn.1006-6047.2016.12.002</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>