Next Article in Journal
Bibliometric Analysis of Multi-Criteria Decision-Making (MCDM) Methods in Environmental and Energy Engineering Using CiteSpace Software: Identification of Key Research Trends and Patterns of International Cooperation
Previous Article in Journal
Impact of Aromatic Hydrocarbons on Emissions in a Custom-Built High-Pressure Combustor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of Two-Stage Economic Optimization Operation of Virtual Power Plants Considering Uncertainty

by
Hao Sun
1,2,
Yanmei Liu
3,
Penglong Qi
1,
Zhi Zhu
1,
Zuoxia Xing
1,* and
Weining Wu
1
1
School of Electrical Engineering, Shenyang University of Technology, Shenyang 110870, China
2
SPIC Northeast Electric Power Co., Ltd., Shenyang 110181, China
3
Material Branch, State Grid Liaoning Electric Power Co., Ltd., Shenyang 110004, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(16), 3940; https://doi.org/10.3390/en17163940
Submission received: 1 July 2024 / Revised: 31 July 2024 / Accepted: 7 August 2024 / Published: 8 August 2024
(This article belongs to the Section C: Energy Economics and Policy)

Abstract

:
In a highly competitive electricity spot market, virtual power plants (VPPs) that aggregate dispersed resources face various uncertainties during market transactions. These uncertainties directly impact the economic benefits of VPPs. To address the uncertainties in the economic optimization of VPPs, scenario analysis is employed to transform the uncertainties of wind turbines (WTs), photovoltaic (PV) system outputs, and electricity prices into deterministic problems. The objective is to maximize the VPP’s profits in day-ahead and intra-day markets (real-time balancing market) by constructing an economic optimization decision model based on two-stage stochastic programming. Gas turbines and electric vehicles (EVs) are scheduled and traded in the day-ahead market, while flexible energy storage systems (ESS) are deployed in the real-time balancing market. Based on simulation analysis, under the uncertainty of WTs and PV system outputs, as well as electricity prices, the proposed model demonstrates that orderly charging of EVs in the day-ahead stage can increase the revenue of the VPP by 6.1%. Additionally, since the ESS can adjust the deviations in day-ahead bid output during the intra-day stage, the day-ahead bidding strategy becomes more proactive, resulting in an additional 3.1% increase in the VPP revenue. Overall, this model can enhance the total revenue of the VPP by 9.2%.

1. Introduction

With the global transition in energy structures and the deepening development of power markets, the electricity spot market has become increasingly significant as a crucial component for the efficient allocation of power resources and market-based transactions [1]. In the electricity spot market, distribution networks serve as the critical link between the generation side and the consumption side, with their operational status and efficiency directly impacting overall market performance and user experience. However, the widespread integration of distributed energy resources, such as WTs and PV systems, along with the diversity and uncertainty of user loads, poses substantial challenges to the secure, stable operation and economic dispatch of distribution networks [2].
VPPs emerge as an innovative energy management technology that aggregates and manages various distributed resources within the distribution network, forming a unified, optimizable virtual entity. This approach provides a novel solution to address uncertainties in distribution networks [3,4]. In the context of the electricity spot market, VPPs not only participate as independent market entities but also leverage their flexible dispatch and resource optimization capabilities to effectively mitigate market uncertainties and achieve economically efficient distribution network operations.
Research on VPPs in electricity spot market transactions has primarily focused on bidding strategies, classifying VPPs as either a “price taker” or “price maker”. When a VPP acts as a “price taker”, it does not influence the price formation process and typically formulates its trading strategy based on price forecasts, making it highly susceptible to price fluctuations. This approach is particularly suitable for VPPs on the distribution network side. Reference [5] proposes a bidding strategy to maximize the overall profit of a VPP by considering day-ahead market clearing prices, power forecasts, and output constraints of distributed wind turbines (WTs) and photovoltaic (PV) systems, along with network structure constraints. Reference [6] introduces a multi-time scale economic dispatch strategy for VPPs, encompassing day-ahead bidding and intra-day operations to enhance economic benefits. Reference [7] presents a two-stage stochastic programming model with risk constraints, addressing the optimal bidding strategy for VPPs in day-ahead, intra-day, and spinning reserve markets. Although the “price taker” strategy facilitates VPPs’ participation in market transactions, developing flexible and effective bidding strategies necessitates comprehensive consideration of market bidding rules, operational costs, competitor bids, and operational risks to maximize profits. Reference [8] incorporates real-time balancing prices and competitor bidding strategies, using conditional value-at-risk to measure bidding risk, and establishes a bi-level optimal bidding model to enhance VPP operational profits. Reference [9] introduces a multi-market model for optimal VPP bidding strategies, converting fixed prices to dynamic response prices to improve demand response efficacy, achieving a win–win for operators and internal consumers. Additionally, VPPs can engage in hydrogen and carbon trading markets, efficiently utilizing renewable energy to reduce carbon emissions and contribute to environmental sustainability [10,11].
A critical research topic in VPPs is formulating economic optimization operation strategies in the electricity spot market while considering uncertainties. Distribution network uncertainties mainly arise from the stochastic nature of renewable energy output, unpredictable user loads, and dynamic market price variations. These uncertainties not only increase the complexity of VPP operations but also impose higher requirements for economic optimization. Currently, there are three primary methods for addressing uncertainties in VPPs: fuzzy theory, robust optimization, and stochastic programming. Reference [12] employ fuzzy theory to systematically represent imprecise and uncertain information, balancing economic and risk factors in VPPs optimization scheduling decisions to enhance operational efficiency and adaptability. Fuzzy theory also improves short-term load forecasting accuracy in VPPs [13]. However, the complexity of VPP system models and multiple uncertainties limit the practical application efficacy of fuzzy optimization [14]. References [15,16] adopt robust optimization, converting models into robust optimization models using strong duality and robust optimization theories, ensuring stable VPP profits under worst-case scenarios. Reference [17] proposes a data-driven robust optimization model to mitigate multiple uncertainties and enhance decision-making, while reference [18] discusses the importance of robust scheduling methods for adapting to dynamic operating conditions in VPPs. In VPP market bidding strategies, reference [19] emphasizes the importance of considering energy storage costs and transmission congestion in robust bidding strategies to improve market competitiveness and adaptability. Reference [20] incorporates demand response mechanisms into VPPs based on robust optimization, proven to enhance operational flexibility and resilience to renewable energy fluctuations. However, references [21,22] indicate that robust optimization often yields overly conservative results, sacrificing some economic benefits and requiring significant computational resources to address multiple uncertainties. Stochastic programming allows for the consideration of uncertainties in renewable energy generation and market prices, enabling VPPs to make informed, robust trading decisions [23]. Reference [24] uses probability density functions of wind speed and solar irradiance distributions, employing scenario generation and reduction techniques to create representative scenarios, describing WT and PV system output uncertainties using typical scenario sets. References [8,25] develop stochastic programming models to determine optimal VPP bidding and pricing strategies, maximizing profits and optimizing dispatch under renewable energy production and demand response uncertainties. The integration of EVs with VPP operations presents new challenges and opportunities. Reference [26] employs stochastic programming to optimize VPPs’ capacity allocation, efficiently meeting EV charging demands, managing EV charging schedules, minimizing costs, and enhancing overall benefits. Recently, advancements in information processing technologies facilitated more accurate acquisition of uncertainty variable probability distributions, reducing computational loads while efficiently utilizing data, thus broadening the applicability of stochastic programming in addressing VPPs uncertainties.
In summary, traditional VPP optimization strategies often overlook the real-time and flexibility requirements of the electricity spot market, particularly when facing challenges posed by uncertain renewable energy output and uncoordinated EV charging. To improve the economic efficiency and stability of VPPs in the spot market, it is essential to propose new optimization strategies to better manage intra-day deviation penalty costs and resource scheduling uncertainties. Therefore, this paper proposes a two-stage optimization strategy for VPPs’ day-ahead and intra-day operations, considering intra-day deviation penalty costs. Specifically, this strategy dispatches gas turbines and EVs in the day-ahead market while deploying flexible ESS in the intra-day real-time balancing market. The proposed bidding strategy and optimization algorithm aim to increase day-ahead market profits, reduce instability from uncoordinated EV charging, and achieve a win–win situation for VPPs and EV users. By using ESS to smooth output deviations, the strategy enhances VPPs’ flexibility and participation enthusiasm in the intra-day market. This comprehensive approach considers multiple distributed resources and their uncertainties, improving overall economic efficiency and operational stability.
The main contributions of this research are as follows:
(1)
Optimized Economic Benefits: Significant enhancement of VPPs market revenue through the rational scheduling of EVs and ESS;
(2)
Improved System Stability: Reduction of power supply instability from uncoordinated EV charging and reliable operation through flexible energy storage dispatch;
(3)
Increased Market Participation: Lowering intra-day market penalty costs and boosting VPPs day-ahead bidding enthusiasm and flexibility, promoting efficient renewable energy utilization;
(4)
Promoted Energy Transition: Providing new theoretical and practical foundations for VPPs applications in the electricity spot market, advancing towards more sustainable and intelligent power systems.
Through this research, we aim to offer new theoretical and practical insights for the optimal operation of VPPs, facilitating their widespread application in power markets and driving energy transition and sustainable development.
The remainder of this paper is structured as follows: Section 2 introduces the VPP structure and electricity market transaction framework. Section 3 details the scenario analysis method. Section 4 presents the two-stage stochastic programming model for the VPP economic optimization. Section 5 provides computational analysis. Finally, Section 6 concludes this paper.

2. The VPP Composition and Electricity Market Trading Framework

2.1. Composition of the VPP

VPPs efficiently manage and utilize diverse energy sources within power systems by integrating traditional power plants, distributed renewable energy, ESS, and flexible loads such as EVs. This integration optimizes energy production and consumption, balancing supply and demand in real time [27]. Such a setup allows VPPs to autonomously schedule activities in the electricity spot market, enhancing their operational capabilities and economic benefits [28,29]. As the cost of energy storage decreases, VPPs can more effectively mitigate the risks associated with renewable energy integration by utilizing ESS and participating in electricity market transactions. This not only improves the overall stability of the power system but also enhances its economic viability. The structure of the VPP constructed in this paper is illustrated in Figure 1. It comprises various components such as WTs, PV systems, ESS, gas turbines, flexible loads, and EVs. As an integrated “source-load-storage” system, the VPP leverages advanced communication, measurement, and control technologies, as well as software systems, to aggregate resources. The VPP scheduling control center coordinates these diverse resource elements, enabling unified participation in grid dispatch and market transactions. This centralized management and optimization control of distributed energy sources, storage systems, controllable loads, and EVs mitigate the impact of uncertainties on the grid. The profitability of a VPP is twofold: externally, it participates in grid dispatch operations and electricity market transactions; internally, it profits from the coordinated and reciprocal scheduling of “source-load-storage” elements.

2.2. WT Power Generation Units

The Weibull distribution model for wind speed [30] is expressed by its probability density function (PDF) as follows:
F ( v ) = 1 e v c k
f ( v ) = k c v c k 1 e v c k
where v represents wind speed, k represents the shape parameter, and c represents the scale parameter. The shape parameter and the scale parameter can be computed using the following expressions:
k = ( σ v ) 1.086 ( 1 k 10 )
c = v Γ ( 1 + 1 k )
where v represents the mean wind speed, σ represents the standard deviation of wind speed, and Γ ( x ) represents the Gamma function.
The expression relating the power output of a WT to wind speed is
P W = 0 v < v ci , v v co k 1 v + k 2 v ci < v v r P r v r < v v co
k 1 = P r v r v ci
k 2 = k 1 v ci
where P W represents the output power generated by the WTs; P r represents the rated power of the WTs; v represents the wind speed affecting the WTs; v c i and v c o , respectively, signify the cut-in and cut-out wind speeds of the WTs; and v r   represents the rated wind speed of the WTs.

2.3. PV System Power Generation Units

It is generally accepted that solar irradiance follows a Beta distribution [31], with the PDF expressed as
f ( γ ) = Γ ( α + β ) Γ ( α ) + Γ ( β ) × ( γ γ m a x ) a 1 × ( 1 γ γ m a x ) b 1
where γ and γ m a x represent the solar irradiance and maximum solar irradiance, respectively. The parameters α and β are the two parameters of the Beta distribution, which can be calculated using the following expressions:
α = w w ( 1 w ) σ 2 1
β = ( 1 w ) [ w ( 1 w ) σ 2 1 ]
The expression relating the output power of a PV system to solar irradiance is
P P V = P p v γ γ e 0 γ γ e P p v γ γ e
where γ represents the rated solar irradiance of the PV system, and P p v represents the rated power of the PV system.

2.4. Gas Turbine Units

The mathematical model for the gas turbine is expressed as
P M T , t = η MT L NG V Δ t
η MT = a × ( p MT ) 3 + b × ( p MT ) 2 + c × ( p MT ) + d
where P M T , t represents the output power of the gas turbine during the time period t; η MT is the efficiency of the gas turbine; L NG represents the lower heating value of natural gas, which is taken as 9.7 kWh/m3 in this paper; V represents the volume of natural gas consumed during time period t; and a, b, c, and d are the fitting coefficients for the gas turbine.

2.5. ESS

The state of charge expression for the charging and discharging of ESS is
S O C = E r a t e E d E r a t e × 100 %
where SOC represents the state of charge of the ESS, E r a t e is the rated capacity of the ESS, and E d is the capacity already released by the ESS.
The state of charge expression for the ESS during charging in time period t is
S t = ( 1 δ t ) S t 1 + P E S S c h a r , t Δ t η c E r a t e
The state of charge expression for the ESS during discharging in time period t is
S t = ( 1 δ Δ t ) S t 1 P E S S d i s , t Δ t E r a t e η d
where S t represents the state of charge of the ESS during time period t; δ is the self-discharge rate; P E S S c h a r , t and P E S S d i s , t represent the charging and discharging power of the ESS during time period t, respectively; and η c and η d represent the charging and discharging efficiencies, respectively.

2.6. EVs Units

EVs have dual attributes as both sources and loads due to their ability to charge and discharge, making them a resource with significant regulatory potential. Therefore, this paper considers EVs as controllable resources within a VPP. After the VPP signs agreements with EV owners, the control center can directly manage the charging and discharging of EVs when they are connected to the grid, according to the optimization strategy.

2.6.1. EVs Travel and Idle Periods

The primary function of EVs in a VPP is to meet the daily travel needs of residents. Given this characteristic, this paper optimizes the charging and discharging of EVs during the period from when owners return home and connect to the grid to when they leave for work and disconnect from the grid. This optimization achieves temporal and spatial energy transfer. Through refined management of EV charging and discharging, not only are the travel needs of owners met and the storage characteristics of EVs effectively utilized, but the economic benefits of the VPP are also enhanced. Based on statistical data, the PDF for the time periods when EV owners connect and disconnect from the grid is expressed as follows [32]:
f dep   ( x ) = 1 2 π σ dep   exp x μ dep   2 2 σ dep   2 , 0 < x μ dep   + 12 1 2 π σ dep   exp x 24 μ dep   2 2 σ dep   2 , μ dep   + 12 < x 24
f arr   ( x ) = 1 2 π σ arr   exp x + 24 μ arr   2 2 σ arr   2 , 0 < x μ arr   12 1 2 π σ arr   exp x μ arr   2 2 σ arr   2 , μ arr   12 < x 24
where x represents a particular time point within the 24 h period, μ dep   and   μ arr   are the expected values of the EVs disconnection and connection times, respectively; and σ dep   and σ arr   are the standard deviations of the EVs’ disconnection and connection times, respectively.

2.6.2. EVs’ Driving Distance

The PDF for the driving distance of EVs is expressed as [33]
L c ( d ) = 1 d σ c 2 π exp In d μ c 2 2 σ c 2
where d represents the driving distance of the EVs, μ c represents the expected value (mean) of the driving distance, and σ c represents the standard deviation of the driving distance.
The relationship between the initial State of Charge (SOC) and the driving distance for EVs is typically expressed as
S O C s t a r t = S O C e n d d P c E E V
where P c represents the EVs’ energy consumption per kilometer, E E V represents the battery capacity of a single EV, S O C e n d represents the state of charge (SOC) at the end of the previous charging session, while S O C s t a r t represents the SOC at the start of the current charging session.

2.6.3. EV Battery Pack Charging and Discharging Model

The state of charge of the EVs during charging in time period t is
S O C E V , t = S O C E V , t 1 + η E V c h a r P E V c h a r , t
The state of charge of the EVs during discharging in time period t is
S O C E V , t = S O C E V , t η E V d i s P E V d i s , t
where S O C E V , t   represents the state of charge of the EVs’ battery in time period t; P E V c h a r , t ,   and P E V d i s , t represent the charging or discharging power of the EVs during time period t; and η E V c h a r and η E V d i s represent the efficiency of charging and discharging for EVs.
After EVs are integrated into the VPP, the control center manages their charging and discharging schedules. However, due to the large number of EVs connected to the VPP, individually calculating the charging and discharging schedules for each vehicle would be computationally intensive. Therefore, this paper uses historical data on EVs to obtain PDFs for their departure times from and driving distances to the grid. Subsequently, employing Latin hypercube sampling and the Ng, Jordan, and Weiss (NJW) spectral clustering method, EVs with similar travel characteristics are clustered together. EVs within the same cluster, characterized by similar travel patterns, are treated as a single scenario set, enabling unified scheduling of each EV scenario set.

2.7. The Bidding Rules and Process for VPPs Participating in Spot Market Trading

The process of VPPs participating in electricity spot market trading is illustrated in Figure 2 as follows:
(1)
Prior to the conclusion of day D in the day-ahead electricity market trading, VPPs forecast the generation capacity of its internal distributed energy resources, internal load demands, and market electricity prices for day D + 1. Based on these forecasts, VPPs optimize the operation of its internal distributed resources with the goal of maximizing operational revenue. Subsequently, VPPs submit their bidding strategy for day D + 1 to the day-ahead electricity trading center (Independent System Operator, ISO). After the bidding process, the day-ahead electricity market operator announces the market clearing price and the awarded bid quantity for VPPs. VPPs then schedule the generation according to the awarded quantity on day D + 1 to ensure completion of the agreed-upon electricity transactions.
(2)
On day D + 1, VPPs control each distributed resource according to the day-ahead awarded bid quantity. Due to the difficulty in accurately predicting distributed WTs and PV system power outputs, there is often a deviation between the actual and planned outputs during operation. Consequently, there exists a discrepancy between the actual output of VPPs and the awarded bid quantity. If the actual output exceeds the awarded bid quantity, VPPs must sell the excess electricity in the intra-day market at a price lower than the market price. Conversely, if the actual operational output is lower than the awarded bid quantity, VPPs need to purchase additional electricity in the real-time balancing market at a price higher than the market price.

3. Scenario Analysis Method

3.1. Scenario Generation and Reduction for WTs and PV System Power Output

In this paper, we address the uncertainties in WTs and PV system power output using stochastic programming, where scenario generation and reduction techniques are employed to create a set of typical scenarios for WTs and PV system power output. This approach transforms the uncertainty problem into a deterministic one, allowing for model-based solutions. As illustrated in Figure 3, we calculate shape parameters based on historical wind speed and solar irradiance data, construct PDFs using Weibull and Beta functions, and employ Latin hypercube sampling and Kantorovich distance reduction methods to generate numerous wind speed and solar irradiance scenarios [34]. Corresponding power output functions are then used to generate typical WTs and PV system power output scenarios, with assigned weights representing the probability of each scenario occurring. This method accurately reflects the WTs and PV system power characteristics of the region, using the expected value of the generated scenarios as the actual output during real-time operation.

3.2. EVs Scenario Clustering

Upon integrating EVs into the VPPs, the VPPs’ control centers manage the charging and discharging schedules of the EVs. Given the large number of EVs, individually calculating and solving the charging and discharging plans for each vehicle would result in a significant computational burden. Therefore, we use historical data on EVs to obtain PDFs for their grid connection and disconnection times and travel distances. We then apply Latin hypercube sampling and the NJW spectral clustering method [35] to cluster EVs with similar travel characteristics. Each cluster, representing EVs with similar characteristics, is treated as a single scenario set, enabling unified scheduling for each EV scenario set.

3.3. Scenario Generation and Reduction for Market Prices

VPPs adopt a “price-taker” bidding strategy. We utilize the K-means clustering method to analyze historical day-ahead and real-time market price data, resulting in a set of typical joint price scenarios for day-ahead and real-time market prices [36].

4. Economic Optimization Decision Model for VPPs Based on Two-Stage Stochastic Programming

4.1. Two-Stage Optimization Strategy for VPPs

This paper assumes that virtual power plants (VPPs) adopt a “price taker” strategy when participating in the electricity spot market, with their bid quantities prioritized for clearing in the market. The participation of VPPs in the electricity spot market involves two stages. The first stage is the day-ahead market, where VPPs optimize the operation strategies of WTs, PV systems, gas turbines, and EVs and make decisions on bidding quantities before the clarity of WTs and PV system power outputs and electricity prices. The second stage is the real-time market, where considering the uncertainty in actual WTs and PV system outputs during the day, VPPs utilize ESS to pre-adjust the deviations in WTs and PV system outputs for each time period, thereby minimizing deviation penalty costs. The specific process is illustrated in Figure 4.

4.2. Objective Function

The profit of VPPs comes from revenue in the day-ahead market and the income (or expenditure) from balancing output deviations in the intra-day market. Due to the uncertainty of WTs and PV system outputs and day-ahead and real-time electricity prices, this paper employs scenario analysis methods to characterize WTs and PV system outputs. The decision model transforms the objective function into maximizing expected revenue across all scenarios, expressed as
R ^ v p p = m α x t = 1 T ω = 1 M π ω ( R t , ω D a + R t E V c h a r C t M T C t E V d i s C t L ) + m α x t = 1 T ω = 1 M π ω s = 1 N 1 γ s ν = 1 N 2 γ ν ( R t , s , ν R t C t E S S )
where w , s , and v represent the sets of scenarios for intra-day electricity prices, WTs, and PV system power outputs, respectively. π w , γ s , and γ v represent the probabilities associated with each scenario set for electricity prices, WT power, and PV system power. R t D a   represents the revenue from VPPs’ participation in day-ahead electricity market sales, R t E V c h a r represents revenue from EV charging within VPPs, and R t , s , ν R t refers to the revenue or costs incurred by VPPs in real-time balancing markets. C t M T , C t E V d i s , C t E S S , and C t L represent the operational costs of gas turbines, subsidies for EV discharging, costs associated with ESS charging and discharging, and adjustments for flexible loads, respectively.

4.2.1. Day-Ahead Electricity Market Revenue

After optimizing the internal power output schedule, VPPs submit the bid quantities for each time period of the next day. Upon successful bidding, each unit generates power according to the planned schedule in real-time. VPPs’ bid quantities and total sales revenue in the day-ahead market can be expressed as follows:
P t D a = P W , t D a + P P V , t D a + P M T , t + i = 1 N P E V d i s , i , t i = 1 N P E V c h a r , i , t P L , t
R t D a = λ t D a P t D a Δ t
where   P W , t D a ,   P P V , t D a , P M T , t , P E V c h a r , i , t , P E V d i s , i , t , and P L , t represent WT power generation forecasting, PV system power generation forecasting, power outputs of gas turbines, EV cluster charging power, EV cluster discharging power, and the load demand of VPPs during time period t, respectively. λ t D a represents the market clearing price in the day-ahead market.
The revenue from EV charging:
R t E V c h a r = λ t E V c h a r i = 1 N P E V c h a r , i , t Δ t
where λ t E V c h a r   represents the charging tariff for EVs.

4.2.2. Revenue or Penalty Costs from Real-Time Electricity Market Transactions

In the second stage, the intra-day market, deviations often exist between planned and actual WTs and PV system outputs (simulated scenarios generate actual outputs in this paper). To mitigate these output deviations, this paper employs ESS as balancing units:
P t , s , v R t = ( P W , t , s R t P W , t D a ) + ( P P V , t , v R t P P V , t D a ) + P E S S d i s , t , P E S S c h a r , t
where   P W , t , s R t and P P V , t , v R t represent the actual WTs and PV system output values composed of the scenario set.
R t , s , v R t = ( 1 θ ) λ t D a Δ P t , s , v + ( 1 + θ ) λ t D a Δ P t , s , v
Δ P t , s , v , + = max P t , s , v R t , 0
Δ P t , s , v = max P t , s , v R t , 0
The equation specifies that θ is the penalty factor, indicating that VPPs will sell below the market price and buy above it in the intra-day market; λ t D a represents the day-ahead market price.

4.2.3. Gas Turbine Fuel Costs

The gas turbines within VPPs consume natural gas during their operation. Therefore, the fuel cost of the gas turbine is calculated as follows:
C t M T = c M T P M T , t Δ t η MT L NG
where c M T represents the unit price of natural gas, η MT , L NG , and P MT are as mentioned earlier.

4.2.4. The Cost of Subsidies for EVs

Frequent control of EV charging and discharging by VPPs can shorten battery lifespan. Therefore, VPPs need to subsidize the discharging of EVs. The cost of subsidizing EVs is
C t E V = c E V d i s i = 1 N ( P E V d i s , i , t ) Δ t
where c E V d i s represents the unit compensation cost for VPPs to control the charging and discharging of EVs, and P E V d i s , i , t represents the discharge power of the EVs cluster at time t.

4.2.5. Operating Cost of ESS

ESS in VPPs shorten their lifespan with each charge and discharge cycle. Therefore, there exists an operating cost during operation:
C t E S S = c E S S ( P E S S c h a r , t + P E S S d i s , t ) Δ t
where c E S S represents the unit operating cost of ESS.

4.3. Constraints

4.3.1. Constraints on WTs and PV System Power Output

The constraints on WT and PV System power output are as shown in Equations (34) and (35).
0 P W , t D a P W , t m a x
0 P P V , t D a P P V , t m a x
where P W , t D a and P P V , t D a represent the actual output of WTs and PV system power during time period t, while P W , t m a x and P P V , t m a x represent the upper limits of WTs and PV system power output during time period t. The upper limits of WTs and PV system power output correspond to the forecasted output values for the day.

4.3.2. Constraints on ESS

The constraints on ESS are as shown in Equations (36)–(39).
α t c P E S S c h a r m i n P E S S c h a r , t α t c P E S S c h a r m a x
α t d P E S S d i s m i n P E S S d i s , t α t d P E S S d i s m a x
S O C E S S m i n S O C t S O C E S S m a x
α t c + α t d 1
where P E S S c h a r m a x and P E S S c h a r m i n represent the maximum and minimum charging power of the ESS, respectively; P E S S d i s m a x and P E S S d i s m i n represent the maximum and minimum discharging power of the ESS, respectively; S O C E S S m a x and S O C E S S m i n represent the maximum and minimum state of charge of the ESS, respectively; and α t   c and α t d represent the charging and discharging efficiencies of the ESS, respectively.

4.3.3. Constraints on Gas Turbine Operation

The constraints on gas turbine operation are as shown in Equations (40) and (41).
P M T min P M T , t P M T max
C M T P MT , t P MT , t 1 C M T +
where P M T max and P M T min represent the maximum and minimum power generation of the gas turbine, respectively; and C M T + and C M T represent the maximum power for the upward and downward ramping of the gas turbine, respectively.

4.3.4. Constraints on EVs

When VPPs control the charging and discharging of EVs according to scheduling strategies, EVs exhibit characteristics similar to ESS to external observers.
u t c P E V c h a r m i n P E V c h a r , i , t u t c P E V c h a r m a x
u t d P E V d i s m i n P E V d i s , i , t u t d P E V d i s m a x
S O C E V , i , t = S O C E V , i , t 1 + η E V c h a r P E V c h a r , i , t P E V d i s , i , t / η E V d i s
S O C E V , i m i n S O C E V , i , t S O C E V , i m a x
u t d + u t c 1
where S O C E V , i , t represents the state of charge of EVs cluster i at time t; S O C E V , i m a x and S O C E V , i m i n represent the maximum and minimum state of charge of EVs cluster i, respectively; P E V c h a r max , P E V c h a r m i n , P E V d i s max , and P E V d i s m i n represent the maximum and minimum charging and discharging power of EVs cluster i, respectively; u t c and u t d are binary variables representing the charging and discharging states of EVs cluster i at time t. These variables cannot both be 1 simultaneously in the same time period t.

5. Case Analysis

5.1. Basic Data

In this section, the VPP system aggregates distributed WTs and PV systems, gas turbines, EVs, ESS, and internal electrical loads. To verify the effectiveness of the two-stage economic optimization decision-making model for the VPP based on stochastic programming constructed in this chapter, the selected equipment parameters and price parameters for the case study are shown in Table 1 and Table 2.
The day-ahead forecast curves for the internal loads, PV power, and WT power of the VPP are shown in Figure 5.
In this paper, the models discussed in Section 3 and Section 4 were implemented using MATLAB 2020a. The Yalmip toolbox within MATLAB 2020a was employed, and the Gurobi solver was used to find solutions, with an average solving time of 125 s.

5.2. Scenario Generation and Clustering Results

5.2.1. WTs and PV System Scenario Generation and Reduction

In this paper, historical wind speed and solar irradiance data from a specific region were selected to construct 10 typical scenarios of WTs and PV system power output. These scenarios were utilized to simulate the actual intra-day output of WTs and PV system power. The generated WTs and PV system power output scenario set are illustrated in Figure 6.
Based on simulation results, the probabilities associated with each typical WT power output scenario are as follows: 0.09, 0.06, 0.10, 0.07, 0.18, 0.10, 0.17, 0.09, 0.05, and 0.09. Similarly, the probabilities for each typical PV system power output scenario are 0.08, 0.10, 0.08, 0.16, 0.06, 0.13, 0.12, 0.07, 0.13, and 0.07.

5.2.2. Clustering of Joint Electricity Price Scenarios

In this paper, K-means clustering was applied to historical electricity price data for day-ahead and intra-day markets. This clustering method yielded 10 typical joint electricity price scenarios, which are used to simulate the electricity prices when the VPP participates in market transactions. The joint electricity price scenario set is illustrated in Figure 7.
Based on simulation results, the probabilities associated with each typical electricity price scenario are as follows: 0.09, 0.06, 0.14, 0.06, 0.15, 0.05, 0.15, 0.10, 0.12, and 0.08.

5.2.3. Clustering Scenarios for Electric Vehicles

Considering the diverse range of EV models on the market, each with unique specifications, this paper simplifies the analysis by assuming all EVs belong to the same model type. The off-grid expected value for EVs is set to 7.45 with a standard deviation of 2.14, while the on-grid expected value is 17.30 with a standard deviation of 2.56. Using Latin hypercube sampling, numerical values for various characteristics of EVs were obtained. These values were then clustered into five clusters using the NJW spectral clustering method, with each cluster’s unified characteristics derived from the average values of all EVs within that cluster. The clustered characteristics of EVs are presented in Table 3.

5.3. Simulation Results Presentation and Analysis

5.3.1. First-Stage Day-Ahead Market Decision Results and Analysis

The internal dispatch strategy of the VPP in the day-ahead market, as depicted in Figure 8, leverages a “price-taker” bidding strategy, providing it with a competitive edge. Consequently, distributed WTs and PV system energy prioritize meeting internal load demands, with surplus electricity bid into the day-ahead market to generate revenue. During periods of low electricity prices, the gas turbine remains inactive due to higher generation costs compared to revenue from external sales, with WT power meeting internal demand. Conversely, during high-price periods, the gas turbine operates at full capacity as its generation costs are lower than revenue from external sales. Considering the gas turbine’s ramping constraints, it increases output continuously during 6–7 and 15–16 h and reduces output during 12–14 and 20–22 h. Additionally, EV discharge (negative) or charge (positive) during low-price periods from 1 to 5 and 22 to 24 h, and discharge during high-price periods from 17 to 21 h. This adjustment increases the VPP’s bid quantity during high-price periods, thereby enhancing its economic returns.
Figure 9 illustrates the correlation between the VPP’s bidding quantity in the day-ahead market, WTs and PV system power outputs, and clearing prices. During periods of higher WTs and PV system outputs, the VPP secures more successful bids in the day-ahead market. Similarly, higher clearing prices prompt increased bidding quantities. For instance, during 6–7, 13–14, and 20–21 h, despite lower clearing prices, the plant’s bids are higher due to greater WTs and PV system outputs. Conversely, during 8–11 and 17–19 h, when clearing prices are higher, the VPP also increases its bidding quantities. This reflects the plant’s ability to flexibly adjust bidding strategies based on market clearing prices and internal aggregation of distributed resources to maximize profits.
Different strategies for charging and discharging among the five clusters of EVs are depicted in Figure 10. The VPP optimally controls charging and discharging for each cluster to maximize benefits. Clusters 1, 2, 3, and 4 charge during periods of low market electricity prices while connected to the grid, meeting daily travel needs. During high-price grid-connected periods, they discharge to increase the VPP’s bid quantity in the day-ahead market, achieving a win–win situation for the plant operator and EV owners. Cluster 5, on the other hand, cannot exploit peak-demand pricing due to low electricity prices during grid-connected periods. Therefore, the VPP only regulates charging for this cluster and does not control discharging activities.

5.3.2. Second-Stage Intra-Day Market Decision Results and Analysis

The deviation in WTs and PV system power output is illustrated in Figure 11, depicting discrepancies between actual and planned outputs within the day. Variations in WTs and PV system output contribute to deviations between the VPP’s actual output during real-time operations and its forecasted bid quantities in the day-ahead market. Figure 12 shows the smoothed output deviations before and after adjustment. Positive deviations in the VPP’s output require selling surplus electricity in the real-time balancing market, while negative deviations necessitate purchasing electricity from the market. Before smoothing, actual WTs and PV system outputs during peak-price periods (8–11 and 17–19 h) are lower than planned outputs, resulting in the VPP needing to purchase additional electricity at higher prices in the real-time balancing market, incurring penalty costs due to deviations. To mitigate these costs, the VPP utilizes ESS to smooth out output deviations. After smoothing, the noticeable reduction in output deviations during 8–11 and 17–19 h indicates effective management of deviations through ESS.
The ESS within the VPP play a pivotal role in mitigating deviations in WTs and PV system output within the intra-day real-time balancing market. Beyond offsetting these deviations, they strategically plan charging and discharging activities based on the intra-day market prices and storage capacity. As depicted in Figure 12 and Figure 13, during low-priced periods, such as 1–2 h, the ESS procure electricity from the real-time balancing market at reduced rates for backup purposes. In periods with lower prices (3–5 h, 7 h, and 14–15 h) and higher actual WTs and PV system outputs compared to forecasts, the systems prioritize storing excess electricity for backup, later selling it when prices are higher. Conversely, during peak-price periods (8–11 h and 17–21 h), where penalties for shortfall in electricity supply are substantial, the ESS prioritize supplying the deficit and selling the remaining capacity at higher market prices, thereby maximizing revenue.

5.3.3. Economic Benefit Analysis

To validate the effectiveness of the proposed two-stage stochastic optimization model for the VPP operation, aiming to enhance its economic efficiency and stability, three scenarios are considered:
Scenario 1: The VPP aggregates WTs, PV systems, and gas turbines participating in both day-ahead and intra-day markets.
Scenario 2: The VPP aggregates WTs, PV systems, gas turbines, and EVs participating in both day-ahead and intra-day markets.
Scenario 3: The VPP aggregates WTs, PV systems, gas turbines, EVs participating in the day-ahead market, and ESS participating in the intra-day market.
Based on Table 4, Table 5 and Table 6, a comparison between Scenarios 1 and 2 reveals that incorporating EVs into the day-ahead market can significantly enhance the VPP earnings. By managing EV charging during off-peak electricity pricing periods to meet daily commuting needs and discharging surplus energy during peak pricing periods, the VPP’s revenue in the day-ahead electricity market increased from USD 7022.7401 to USD 7268.573, marking a gain of USD 245.8329. Although subsidizing EV charging imposes costs on vehicle owners, the VPP can strategically sell surplus energy during peak pricing, resulting in an overall revenue increase from USD 5696.3104 to USD 6041.7065—a rise of USD 345.3961, or approximately 6.1%. Moreover, implementing controlled EV charging and discharging schedules effectively mitigates the supply instability caused by unregulated EV charging.
Comparing Scenarios 2 and 3, where ESS are absent from the intra-day market, uncertainty in WTs and PV system outputs leads to penalty costs during periods of energy shortfall in the real-time balancing market. To minimize these penalties, Scenario 2 adopts a conservative strategy in the day-ahead market, reducing bid quantities. In contrast, Scenario 3 integrates ESS into the intra-day market, smoothing out output deviations and reducing penalty costs. Consequently, the VPP adopts a more proactive bidding strategy in the day-ahead market, increasing bid quantities and raising earnings from USD 7268.573 to USD 7435.6843—an increase of USD 167.1113. Additionally, by optimizing ESS charging and discharging strategies, earnings in the intra-day balancing market also improve, lifting the VPP’s total revenue from USD 6041.7065 to USD 6231.4033—a rise of USD 189.6967, or about 3.1%. Introducing ESS into the intra-day market enhances the VPP flexibility in day-ahead market decisions, thereby boosting system stability and overall operational profitability.

6. Conclusions

This paper introduces an innovative two-stage optimization strategy for the day-ahead and intra-day operation of VPPs, which incorporates intra-day deviation penalty costs. The strategy aims to achieve efficient VPP operation in the electricity market through precise scheduling of gas turbines, EVs, and flexible ESS. In the day-ahead stage, gas turbines and EVs are included in the scheduling framework. By leveraging the flexible charging and discharging capabilities of EVs, the strategy significantly enhances the VPP’s market revenue while mitigating the negative impact of uncoordinated EV charging on grid stability, thus creating a win–win situation for both VPP operators and EV users. In the intra-day market, ESS smooths out discrepancies between day-ahead forecasts and actual outputs. This approach successfully reduces the penalty costs faced by VPPs in the intra-day market. Simulation results indicate that through orderly management of EV charging in the day-ahead stage, the VPP’s revenue increased by 6.1%. In the intra-day stage, thanks to precise regulation by ESS, the strategy effectively mitigated deviations between day-ahead bids and actual outputs, resulting in an additional 3.1% revenue increase. Consequently, the overall revenue of the VPP improved by 9.2%. The proposed two-stage optimization strategy for VPPs provides precise scheduling of various resources, effectively addressing uncertainties in the electricity market. This significantly enhances the economic efficiency and stability of VPP operations, offering new insights and methods for the economic optimization of future VPPs.

Author Contributions

H.S. led the investigation, designed and implemented the methodology, and wrote the initial draft of the manuscript, providing critical insights and analysis. Y.L. played a key role in developing the software tools used in the study and contributed significantly to drafting the original manuscript. P.Q. was primarily responsible for refining the methodology and ensuring its robust application throughout the research. Z.Z. conducted the formal analysis, interpreting the results and contributing to the accuracy of the findings. Z.X. provided overall supervision of the research process and contributed to the review and editing of the manuscript to enhance its quality and coherence. W.W. offered valuable supervision and managed the necessary resources, ensuring the smooth execution of the study. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

Author Hao Sun was employed by the company SPIC Northeast Electric Power Co., Ltd., Author Yanmei Liu was employed by the company Material Branch, State Grid Liaoning Electric Power Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. 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.

References

  1. Liu, C.; Yang, R.; Wang, K.; Zhang, J. Community-Focused Renewable Energy Transition with Virtual Power Plant in an Australian City—A Case Study. Buildings 2023, 13, 844. [Google Scholar] [CrossRef]
  2. Cao, W.; Yu, J.; Xu, M. Optimization Scheduling of Virtual Power Plants Considering Source-Load Coordinated Operation and Wind–Solar Uncertainty. Processes 2024, 12, 11. [Google Scholar] [CrossRef]
  3. Ullah, Z.; Arshad; Hassanin, H. Modeling, Optimization, and Analysis of a Virtual Power Plant Demand Response Mechanism for the Internal Electricity Market Considering the Uncertainty of Renewable Energy Sources. Energies 2022, 15, 5296. [Google Scholar] [CrossRef]
  4. Cao, W.; Wang, S.; Xu, M. Optimal Scheduling of Virtual Power Plant Based on Latin Hypercube Sampling and Improved CLARA Clustering Algorithm. Processes 2022, 10, 2414. [Google Scholar] [CrossRef]
  5. Song, M.; Amelin, M. Purchase bidding strategy for a retailer with flexible demands in day-ahead electricity market. IEEE Trans. Power Syst. 2017, 32, 1839–1850. [Google Scholar] [CrossRef]
  6. Yi, Z.; Xu, Y.; Gu, W.; Wu, W. A Multi-time-scale Economic scheduling strategy for virtual power plant based on deferrable loads aggregation and disaggregation. IEEE Trans. Sustain. Energy 2020, 11, 1332–1346. [Google Scholar] [CrossRef]
  7. Vahedipour-Dahraie, M.; Rashidizade-Kermani, H.; Shafie-Khah, M.; Catalão, J.P. Risk-Averse Optimal Energy and Reserve Scheduling for Virtual Power Plants Incorporating Demand Response Programs. IEEE Trans. Smart Grid 2020, 12, 1405–1415. [Google Scholar] [CrossRef]
  8. Kardakos, E.G.; Simoglou, C.K.; Bakirtzis, A.G. Optimal Offering Strategy of a Virtual Power Plant: A Stochastic Bi-Level Approach. IEEE Trans. Smart Grid 2016, 7, 794–806. [Google Scholar] [CrossRef]
  9. Mei, S.; Tan, Q.; Liu, Y.; Trivedi, A.; Srinivasan, D. Optimal bidding strategy for virtual power plant participating in combined electricity and ancillary services market considering dynamic demand response price and integrated consumption satisfaction. Energy 2023, 284, 128592. [Google Scholar] [CrossRef]
  10. Wu, Q.; Li, C. A bi-level optimization framework for the power-side virtual power plant participating in day-ahead wholesale market as a price-maker considering uncertainty. Energy 2024, 304, 132050. [Google Scholar] [CrossRef]
  11. Tong, Y.; Meng, Z.; Qiu, P.; Cui, J.; Ma, J.; Tian, Y.; Wang, Z. The carbon trading operation optimization for virtual power plants of industrial parks considering wind power. J. Phys. Conf. Ser. 2023, 2474, 012032. [Google Scholar] [CrossRef]
  12. Zhang, L.; Zheng, H.; Cai, G.; Zhang, Z.; Wang, X.; Koh, L. Power-frequency oscillation suppression algorithm for ac microgrid with multiple virtual synchronous generators based on fuzzy inference system. IET Renew. Power Gener. 2022, 16, 1589–1601. [Google Scholar] [CrossRef]
  13. Islas, M.; Rubio, J.; Muñiz, S.; Ochoa, G.; Pacheco, J.; Meda-Campaña, J.; Mujica-Vargas, D.; Aguilar-Ibañez, C.; Gutierrez, G.J.; Zacarías, A. A fuzzy logic model for hourly electrical power demand modeling. Electronics 2021, 10, 448. [Google Scholar] [CrossRef]
  14. Abbasi, M.H.; Taki, M.; Rajabi, A.; Li, L.; Zhang, J. Coordinated operation of electric vehicle charging and wind power generation as a virtual power plant: A multi-stage risk constrained approach. Appl. Energy 2019, 239, 1294–1307. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Liu, F.; Wang, Z.; Su, Y.; Wang, W.; Feng, S. Robust scheduling of virtual power plant under exogenous and endogenous uncertainties. IEEE Trans. Power Syst. 2022, 37, 1311–1325. [Google Scholar] [CrossRef]
  16. Ren, L.; Peng, D.; Wang, D.; Li, J.; Zhao, H. Multi-objective optimal dispatching of virtual power plants considering source-load uncertainty in v2g mode. Front. Energy Res. 2023, 10, 983743. [Google Scholar] [CrossRef]
  17. Fang, F.; Yu, S.; Xin, X. Data-driven-based stochastic robust optimization for a virtual power plant with multiple uncertainties. IEEE Trans. Power Syst. 2022, 37, 456–466. [Google Scholar] [CrossRef]
  18. Naughton, J.; Han, W.; Cantoni, M.; Mancarella, P. Co-optimizing virtual power plant services under uncertainty: A robust scheduling and receding horizon dispatch approach. IEEE Trans. Power Syst. 2021, 36, 3960–3972. [Google Scholar] [CrossRef]
  19. Tan, Z.; Tan, Q.; Wang, Y. Bidding strategy of virtual power plant with energy storage power station and photovoltaic and wind power. J. Eng. 2018, 2018, 6139086. [Google Scholar] [CrossRef]
  20. Dong, J. Decision-making model of virtual power plant for participating in spot market transaction based on hybrid stochastic and robust approach. Am. J. Environ. Resour. Econ. 2019, 4, 32. [Google Scholar] [CrossRef]
  21. Alahyari, A.; Pozo, D. Performance-based virtual power plant offering strategy incorporating hybrid uncertainty modeling and risk viewpoint. Electr. Power Syst. Res. 2022, 203, 107632. [Google Scholar] [CrossRef]
  22. Alahyari, A.; Ehsan, M.; Pozo, D.; Farrokhifar, M. Hybrid uncertainty-based offering strategy for virtual power plants. IET Renew. Power Gener. 2020, 14, 2359–2366. [Google Scholar] [CrossRef]
  23. Ji, W.; Wang, Y.; Deng, X.; Zhang, M.; Ye, T. Distributionally robust optimal dispatch of virtual power plant based on moment of renewable energy resource. Energy Eng. 2022, 119, 1967–1983. [Google Scholar] [CrossRef]
  24. Ju, L.; Tan, Z.; Yuan, J.; Tan, Q.; Li, H.; Dong, F. A bi-level stochastic scheduling optimization model for a virtual power plant connected to a wind–photovoltaic–energy storage system considering the uncertainty and demand response. Appl. Energy 2016, 171, 184–199. [Google Scholar] [CrossRef]
  25. Geng, S.; Tan, C.; Niu, D.; Guo, X. Optimal allocation model of virtual power plant capacity considering electric vehicles. Math. Probl. Eng. 2021, 2021, 5552323. [Google Scholar] [CrossRef]
  26. Sun, W.; Liu, S.; Dong, H.; Huang, Q. Electric power dispatching of virtual power plant with electric vehicle. J. Phys. Conf. Ser. 2022, 2409, 012019. [Google Scholar] [CrossRef]
  27. Baringo, A.; Baringo, L.; Arroyo, J. Day-ahead self-scheduling of a virtual power plant in energy and reserve electricity markets under uncertainty. IEEE Trans. Power Syst. 2019, 34, 1881–1894. [Google Scholar] [CrossRef]
  28. Zhong, W.; Chen, J.; Liu, M.; Murad, M.; Milano, F. Coordinated control of virtual power plants to improve power system short-term dynamics. Energies 2021, 14, 1182. [Google Scholar] [CrossRef]
  29. Yu, X.; Zhao, W. Participation of digital twin virtual power plants in power market trading model. Int. J. Manag. Stud. Soc. Sci. Res. 2023, 5, 131–135. [Google Scholar] [CrossRef]
  30. Shukla, A.; Verma, K.; Kumar, R. Impact of EV fast charging station on distribution system embedded with wind generation. J. Eng. 2019, 2019, 4692–4697. [Google Scholar] [CrossRef]
  31. Karaki, S.H.; Chedid, R.B.; Ramadan, R. Probabilistic performance assessment of autonomous solar-wind energy conversion systems. IEEE Trans. Energy Convers. 1999, 14, 766–772. [Google Scholar] [CrossRef]
  32. Li, Y.; Zhang, S.; Xiao, X.; Wang, Y. Charging and discharging scheduling strategy of EVs considering demands of supply side and demand side under V2G mode. Electr. Power Autom. Equip. 2021, 41, 129–135. [Google Scholar]
  33. Qian, K.; Liu, Y.; Zhang, X. Scenario-based probabilistic power flow calculation of distribution systems with electric vehicle charging loads. Power Syst. Prot. Control 2020, 48, 62–70. [Google Scholar]
  34. Dupacova, J.; Growekuska, N.; Romisch, W. Scenario reduction in stochastic programming: An approach using probability metrics. Math. Program. 2000, 95, 493–511. [Google Scholar] [CrossRef]
  35. Wang, R.; Wu, J.; Cai, Z.; Liu, G.; Zhang, H.; Cai, J. Optimal Dispatch of Virtual Power Plant Containing Electric Vehicles in a Multi-coordinated Market. South. Power Syst. Technol. 2021, 15, 45–55. [Google Scholar]
  36. Cai, X.; Cui, Y.; Dong, K.; Meng, Z.; Pan, Y.; Yu, Z.; Wang, J.; Meng, X.; Yu, Y. Day-ahead optimal scheduling approach of wind-storage joint system based on improved K-means and MADDPG algorithm. Energy Storage Sci. Technol. 2021, 10, 2200–2208. [Google Scholar]
Figure 1. Structure of the VPP.
Figure 1. Structure of the VPP.
Energies 17 03940 g001
Figure 2. Process of VPPs participation in electricity market bidding and trading.
Figure 2. Process of VPPs participation in electricity market bidding and trading.
Energies 17 03940 g002
Figure 3. Flow chart of scenario analysis method.
Figure 3. Flow chart of scenario analysis method.
Energies 17 03940 g003
Figure 4. Two-stage economic optimization operation process for VPPs.
Figure 4. Two-stage economic optimization operation process for VPPs.
Energies 17 03940 g004
Figure 5. Day-ahead forecast curves for internal loads, WT power, and PV system power.
Figure 5. Day-ahead forecast curves for internal loads, WT power, and PV system power.
Energies 17 03940 g005
Figure 6. Typical scenario set of WTs’ power output.
Figure 6. Typical scenario set of WTs’ power output.
Energies 17 03940 g006
Figure 7. Day-ahead and intra-day market joint tariff scenario sets.
Figure 7. Day-ahead and intra-day market joint tariff scenario sets.
Energies 17 03940 g007
Figure 8. Internal scheduling strategy for VPP in the day-ahead electricity market.
Figure 8. Internal scheduling strategy for VPP in the day-ahead electricity market.
Energies 17 03940 g008
Figure 9. Day-ahead market best bid electricity.
Figure 9. Day-ahead market best bid electricity.
Energies 17 03940 g009
Figure 10. Charging and discharging strategies for different classes of EVs.
Figure 10. Charging and discharging strategies for different classes of EVs.
Energies 17 03940 g010
Figure 11. WTs and PV system output deviation.
Figure 11. WTs and PV system output deviation.
Energies 17 03940 g011
Figure 12. Flatten the deviation of front and rear output.
Figure 12. Flatten the deviation of front and rear output.
Energies 17 03940 g012
Figure 13. ESS charging and discharging strategies.
Figure 13. ESS charging and discharging strategies.
Energies 17 03940 g013
Table 1. Equipment parameters.
Table 1. Equipment parameters.
EquipmentParameterValue
Gas turbineMaximum/minimum power generation (MW)1.8
Power generation efficiency0.7
Up and down climbing rate (MW/h)0.5
EVsBattery capacity of a single EV (MWh)0.05
Power consumption per kilometer (MWh)0.00015
Maximum and minimum State of Charge0.9, 0.1
Maximum charging/discharging power (WM)0.01
Minimum charging/discharging power (MW)0
ESS
(composition of lithium batteries)
Rated capacity (MWh)4
Maximum and minimum state of charge0.9, 0.1
Maximum charging power (MW)0.4
Minimum charging power (MW)0
Table 2. Price parameters.
Table 2. Price parameters.
ParameterValue (USD)
EV charging price (USD/MWh)55.1436
EV discharge subsidy electricity price (USD/MWh)68.9295
Operating cost of ESS charging and discharging (USD/MWh)6.88
Natural gas (USD/m3)0.4136
Penalty factor0.2
Table 3. Characteristics of EVs after clustering.
Table 3. Characteristics of EVs after clustering.
ClusterOff-Grid Time (Hours/Minutes)On-Grid Time (Hours/Minutes)On-Grid SOC Status
13:5214:460.714
25:5516:260.645
37:3417:350.626
49:1220:240.668
511:1321:410.639
Table 4. Day-ahead market benefits.
Table 4. Day-ahead market benefits.
ScenarioDay-Ahead Market Revenue (USD)EV Charging Revenue (USD)Gas Turbine Fuel Cost (USD)EV Discharging Subsidy Cost (USD)
one7022.740101470.96930
two7268.573437.5591470.9693337.9958
three7435.6843437.5591470.9693337.9958
Table 5. Intra-day market benefits.
Table 5. Intra-day market benefits.
ScenarioReal-Time Balancing Market Revenue (USD)ESS Operating Cost (USD)
one144.53960
two144.53960
three219.99452.8689
Table 6. Total benefits of the VPP.
Table 6. Total benefits of the VPP.
ScenarioDay-Ahead Market Total Revenue (USD)Real-Time Balancing Market Total Revenue (USD)Total VPP Revenue (USD)
one5551.7708144.53965696.3104
two5897.1669144.53966041.7065
three6064.2782167.12516231.4033
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, H.; Liu, Y.; Qi, P.; Zhu, Z.; Xing, Z.; Wu, W. Study of Two-Stage Economic Optimization Operation of Virtual Power Plants Considering Uncertainty. Energies 2024, 17, 3940. https://doi.org/10.3390/en17163940

AMA Style

Sun H, Liu Y, Qi P, Zhu Z, Xing Z, Wu W. Study of Two-Stage Economic Optimization Operation of Virtual Power Plants Considering Uncertainty. Energies. 2024; 17(16):3940. https://doi.org/10.3390/en17163940

Chicago/Turabian Style

Sun, Hao, Yanmei Liu, Penglong Qi, Zhi Zhu, Zuoxia Xing, and Weining Wu. 2024. "Study of Two-Stage Economic Optimization Operation of Virtual Power Plants Considering Uncertainty" Energies 17, no. 16: 3940. https://doi.org/10.3390/en17163940

APA Style

Sun, H., Liu, Y., Qi, P., Zhu, Z., Xing, Z., & Wu, W. (2024). Study of Two-Stage Economic Optimization Operation of Virtual Power Plants Considering Uncertainty. Energies, 17(16), 3940. https://doi.org/10.3390/en17163940

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop