Next Article in Journal
The Quest for Renewable Energy—Effects of Different Asphalt Mixes and Laboratory Loading on Piezoelectric Energy Harvesters
Previous Article in Journal
Retrofit Proposals for Energy Efficiency and Thermal Comfort in Historic Public Buildings: The Case of the Engineering Faculty’s Seat of Sapienza University
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment of User Aggregators in Demand Bidding Markets

1
Department of Electrical Engineering, Cheng-Shiu University, Kaohsiung 833, Taiwan
2
School of Mechanical and Electrical Engineering, Tan Kah Kee College, Xiamen University, Zhangzhou 363105, China
*
Author to whom correspondence should be addressed.
Energies 2023, 16(1), 156; https://doi.org/10.3390/en16010156
Submission received: 30 November 2022 / Revised: 17 December 2022 / Accepted: 19 December 2022 / Published: 23 December 2022

Abstract

:
This paper mainly discusses the demand bidding and risk management of user aggregators by considering profit and risk. The covariance matrix of demand price was used to analyze the risk model under an uncertain demand price. By considering revenue and cost, the demand bidding strategy of user aggregators was derived to obtain the maximum profit. By using a risk-tolerance parameter β, a new demand bidding model for the user aggregators that takes both risk and profit into consideration was formulated. We simulated the risk posed by fluctuating demand prices for user aggregators using this model. Finally, this paper proposes Feasible Particle Swarm Optimization (FPSO) to solve the demand bidding model of user aggregators. Through the dynamic adjustment of control factor parameters in the FPSO, we changed the behavioral characteristics of various types of particles, improved the search efficiency and stability of particles in high-dimensional space, and sought the optimal solution for the system as a whole. This paper provides a parameter adjustment mechanism, improves the capability of algorithm implementation, and increases the probability of finding the optimal solution. The simulation results suggest that a tradeoff between profit and risk needs to be considered in the search process. By doing so, enterprises’ abilities in terms of operation and management control can be enhanced, and effective demand management can be achieved.

1. Introduction

In recent years, the overall investment environment in Taiwan has been booming, and electricity consumption has reached a record high every year. The gap between peak load and off-peak load continues to increase every year. The ratio between the summer peak load during the day and the off-peak load at night is about 1:0.61, and the maximum summer load is about 1.4 times that of the maximum non-summer load [1]. This shows that there is still a great difference in load between different seasons and between the peak load and off-peak load of the power system. In order to reflect the power supply cost difference between the peak and off-peak power consumption in summer and non-summer months, the Taiwan Power Company (TPC) implements an electricity price system based on seasons and time to control the peak load and improve the power supply [2]. Electricity users can cooperate with the power company to implement load management and reduce power consumption in peak periods without affecting their electricity demand or normal operation. Therefore, it could improve their willingness to participate the implementation performance of this system. In addition, electricity users could also obtain some revenue due to favorable feedback and compensation from the power company. All of these could be more beneficial through the utilization rate of power equipment and enhancing the quality of electricity utilization. If operators could effectively implement these to reduce the problem of excessive peak load through shifting the electricity consumption during peak periods to off-peak periods, achieve power balance, and reduce the construction demand of new power plants, it would effectively improve the electricity supply and ensure that power companies and enterprise users achieve a mutually beneficial situation. Therefore, for load management, the research and development of a demand management strategy is one of the most important policy directions [3].
The demand response (DR) strategy is one of the measures taken by power companies to reduce electricity consumption, and this strategy could provide incentives to reduce electric charges. During periods when the power supply is low or the cost is high, electricity users could cooperate to reduce the agreed electricity load and consumption, which could have a very significant effect on system load reduction in peak periods. This strategy could reduce the burden of electricity charges on users for electricity consumption, and it has gradually started to attract the attention of enterprise users. In the demand response mechanism, such as the planned electricity reduction measures or demand bidding strategy [4], the electricity load and consumption reduced by electricity users could obtain high benefits through participating in the demand bidding strategy. Under the framework of a demand bidding strategy, power companies could provide electricity price discounts or other market incentive mechanisms, and the aggregator representative of electricity users could reduce the electricity load and consumption during periods of low peak electricity supply for many users according to the needs of the power company by means of demand bidding. The concept of Virtual Peaking Capacity (VPC) [5] could generate electric resources not only from the central power system but also from the distributed power generation system at the power load end under the “real-time” and “two-way” approaches. Therefore, in comparison with traditional systems, the electricity supply and demand are multi-variant and more flexible. In addition, through the implementation of energy saving plans, it could reduce the electricity consumption demand of power users.
In recent years, demand management strategies for the electricity market have mainly focused on the effects of various operation strategies on profit. Some studies have effectively demonstrated how to solve these problems in the electricity market. Reference [6] incorporated price-based demand response to perform the optimal dispatch schedule in the day-ahead electricity market. A new demand bidding strategy for the smart building aggregator was proposed to perform demand response in the electricity market [7]. References [8,9] proposed a bi-level model for the demand bidding strategy of load agents with incentive-based demand response in day-ahead electricity markets. Based on game theory, [10] proposed an optimal bidding strategy for demand response in the electricity market. A dynamic demand bidding strategy for an aggregator was proposed to participate in the frequency regulation market [11]. A demand response model between the integrated energy production base, load aggregator, and user is established, which is solved by a mixed-integer quadratic programming–multi-verse optimizer distributed algorithm [12]. Reference [13] proposed a novel approach for incorporating incentive-based and price-based demand response programs in long-term generational investment planning. Reference [14] established an energy management platform by integrating distributed power and energy storage to perform the demand bidding strategy. The common disadvantage of the methods described is the lack of guarantee that the profit may be at risk due to the uncertainties of the electricity price. While some risk was expected in the formation and operation of demand response, addressing uncertainties would be a key issue in studying this topic.
Risk management has become an essential condition for the sustainable operation of the power generation industry [15,16]. The demand bidding strategy of user aggregators is particularly important for the degree of risk of profit [17]. The risk assessment would analyze the possible risks faced by each system during power dispatching and pricing and discuss the maximum profit of demand bidding for user aggregators and the corresponding risk after considering the risk limit. However, in the demand bidding strategy, representatives of users could integrate many different electricity users and sign contracts with them under the participation conditions of demand response [18,19,20]. This study analyzed the demand bidding of user aggregators and considered the maximum profit after taking profit and risk into consideration. The demand price variability is used to generate demand price variation to establish a covariance matrix of demand price [21], which is used to analyze the risk level of demand bidding for user aggregators. The effect of risk is explicitly introduced in the bidding strategy problem, considering the variance in demand prices. The tradeoff model between profit and risk is properly addressed within a given risk tolerance. Finally, this paper proposed a Feasible Particle Swarm Optimization (FPSO) to solve the demand bidding model of user aggregators. In the FPSO procedure, the dynamic control parameters are embedded in the particle swarm of the FPSO to improve the behavior patterns of each particle swarm and increase its search efficiency and accuracy in high dimensions. Different modifications in the moving patterns of FPSO are proposed to search the feasible space more effectively. The results can help decision-makers optimize the tradeoffs between maximum profit and minimum risk. In addition to helping decision-makers improve system operation safety, energy efficiency, and risk management, this study can also explore the economic benefits of different decisions in the architecture and planning of energy management systems.

2. Risk Model for Demand Bidding

In Taiwan, the model of planning electricity consumption reduction allows users to evaluate their characteristics for obtaining their maximum profit. Three models were developed by the Taiwan Power Company, which included an 8 day monthly reduction model, a 6 h daily reduction model, and a 2 h daily reduction model [22]. The implementation period was from June 1 to September 30 each year [22]. In the 8 day monthly reduction model, the suppression of electricity consumption occurs from Monday to Friday every month, and the suppression of electricity consumption time is 10 a.m. to 5 p.m. In the 6 h daily reduction model, the suppression of electricity consumption occurs from Monday to Friday every month, and the suppression of electricity consumption times are 10–12 a.m. and 1–5 p.m. In the 2 h daily reduction model, the suppression of electricity consumption occurs from Monday to Friday every month, and the suppression of electricity consumption time is 1–3 p.m. This paper adopted the 6 h daily reduction model to analyze the risk of demand bidding by user aggregators. The suppression of electricity consumption occurs from 10 a.m. to 12 p.m. and from 1 p.m. to 5 p.m. daily during non-holidays. This study utilized the uncertain demand price in the electricity market to generate the demand price variation and establish the covariance matrix of the demand price to analyze the risk of demand bidding by user aggregators. The procedure is described as follows:
  • In the day-ahead demand-bidding market, the revenue of the user aggregators is defined as follows:
    R e v n u e = t = 1 T λ t p t
    in the day-ahead demand bidding market, the average revenue of the user aggregator can be obtained by adding the power demand that can be sold at the estimated power expectation for each period (the total demand amount of user aggregators). The expected revenue operator in total periods can be changed as follows:
    R exp = t = 1 T E λ t { λ t } p t
    where p t is the power demand (MW) bid by the user aggregator during the period t , and R e s t is the expected total revenue of the user aggregator. λ t is the demand price of the user aggregator during the period t . T is the total period of demand bidding (the total peak period was 6 h which included from 10 a.m. to 12 p.m. and from 1 p.m. to 5 p.m. E λ T is the expected value operator of the random variable λ T ;
  • During the 6 h of the peak period, the profits of each period affect each other. The total variation value of revenue during the peak period of 6 h is calculated as follows:
    V a r λ 1 , , λ T { t = 1 T λ t p t } = t = 1 T j = 1 T p i V i j p j
    V a r λ 1 , , V a r λ T is the variance operator of the random variable.
    V = ( V i j ) is the T × T covariance matrix of the demand price λ 1 , , λ T .
    p i and p j are the power demand (MW) bid by the user aggregators, and they are converted into rows and columns. Because the demand price is the only random variable, the variation of total profit can be expressed by taking the demand price as the covariance matrix. The covariance matrix of the d day is:
    V = E λ 1 , , λ T { ( d true d est ) ( d true d est ) T }
    Λ d = [ λ 1 , , λ T ] T ;
  • If the bidding history data of the demand trading market are collected up to d 1 day, the covariance matrix formula of d day can be expressed by the actual and predicted values as follows:
    V est = 1 D i = 1 D ( Λ i true Λ i est ) ( Λ i true Λ i est ) T
    where ture represents the superscript of the actual value, and D is the total number of days that is the greatest and contains d 1 day. When Equation (5) is directly used, due to the nature of the price characteristics of the demand trading market, they may have multiple seasonal characteristics and high variability, as well as unusual purchase prices affected by high load;
  • In order to obtain an accurate prediction, it could be modified with the exponentially weighted moving average equation [23]:
    V est = ( 1 α ) i = 1 D α i 1 ( Λ D i + 1 true Λ D i + 1 est ) ( Λ D i + 1 true Λ D i + 1 est ) T
    where Λ is the data on the demand price of past user aggregators, which is multiplied by the weight value α ( 0 < α < 1 ) . The closer to the estimated d day, the greater the weight value; the further from the estimated d day, the more exponentially the weight value decays. Therefore, old data have less influence on the variance and covariance because they generate outliers due to excessive load. Equations (5) and (6) are both modified formulas representing the covariance matrix V . However, the larger the D value, the smaller the unreasonable estimation offset. In addition, the smaller the estimation offset, the more accurate the estimation result;
  • Regarding the demand bidding planning of the user aggregators, the demand planning strategy for maximum profit can be formulated as follows:
    maximize p 1 , p T t = 1 T ( λ t est p t n = 1 N U t , n F c t t , n ) , p 1 , , p T Π
    F c t n = a n P c t n 2 + b n P c t n + c n
    where U t , n is the purchase status of the user aggregators for downstream users during the period t . F c t t , n and P c t n are the cost function and electricity amount of the user aggregators in the demand bidding, respectively. a , b , and c are the demand bidding curves of user aggregators. Π is a feasible solution. N is the number of user aggregators who took part in the demand bidding during the time period t . The demand bidding must face the tradeoff between the maximum profit and the minimum risk. If the user aggregators seek to minimize risk regardless of profit, the risk minimization procedure can be formulated as follows:
    minimize p 1 , p T i = 1 T j = 1 T p i V i j est p j , p 1 , , p T Π ;
  • V est is the risk variation of the demand price of the utility during the peak hours, while p i and p j are the demand bids by the user aggregators. The user aggregators are most interested in the best demand bidding to make profits, and the demand bidding of these user aggregators have the maximum profit with the minimum risk. In order to compromise these two conflicting goals, the best choice is complemented by a risk-tolerance parameter β . Therefore, the demand bidding for the user aggregators, taking both risk and maximum profit into consideration, can be formulated as Equation (10):
    maximize p 1 , p T t = 1 T ( λ t est p t i = 1 N U t , i F c t t , i ) β i = 1 T j = 1 T p i V i j est p j , p 1 , , p T Π
    in addition, it should meet the restriction conditions of period t as in Equation (11):
    n = 1 N U t , n F c t t , n = P t , i t = 1 , 2 , , T
    profit and risk are two different objectives. Hence, a compromise solution is required to solve the demand bidding problem. Equation (10) is an objective function in this paper. The objective function adopted the concepts of “profit” and “risk” to obtain the maximum profit by controlling the “risk”. When the risk-tolerance parameter β is higher, the profit is lower.

3. Feasible Particle Swarm Optimization

In 1995, PSO was developed by [24] to simulate numerical analysis. The PSO algorithm is based on the mechanism of birds’ swarming behavior. When birds forage, each bird is like a particle. Each particle keeps in mind its current best position ( p b e s t ), and the best position of the group (gbst) in the population. The disadvantage of PSO is that the convergence speed is higher in the initial search for solving the optimization problem. In the later stage, the particle swarm gradually moves toward the optimal solution of the swarm. Thus, the diversity of the whole swarm is lost, and the particle can easily fall into the local optimal solution.
In PSO, the position and velocity of particles are defined in Equations (12) and (13):
v i t + 1 = w × v i t + c 1 × r a n d × ( p b e s t i t x i t ) + c 2 × r a n d × ( g b e s t t x i t )
x i t + 1 = x i t + v i t + 1
where x i t is the position of particle i at iteration t , and v i t is the velocity of particle i during iteration t . p b e s t i t is the best position of a particle i at iteration t , and g b e s t t is the best position of all particles at iteration t . c 1 and c 2 are learning constants that influence the forward speed of the particle. In this paper, c 1 and c 2 are 2.05. w is the speed weight of the current generation, and w means a small variance when the particle changes its position; otherwise, the variance is large.
PSO has the characteristics of movement, evolution, elimination, and multiple variability. The movement in the PSO process is only dependent on the movement distance generated, and the movement information is transmitted between the particles. It is easy to converge on the local optimum. FPSO introduces an auto-tuning scheme that allocates “infeasible solution ( g b e s t inf )” and “feasible solution ( g b e s t f )” to the solution space by using random parameters to make the PSO search for the global optimum more efficient, as shown in Figure 1. In Figure 1, the “infeasible solution ( g b e s t inf )” space may be an attractive solution in the next generation, and improvement of the auto-tuning solution starts with finding the best direction to take advantage of more opportunity to obtain the global optimum. The velocity of particles in the FPSO system is defined in Equation (14):
i f r < 0.3 t h e n v i , j t + 1 = w × v i , j t + c 1 × r a n d × ( p b e s t i , j x i , j t ) + c 2 × r a n d × ( g b e s t inf , j x i , j t ) e l s e i f 0.3 r 0.7 t h e n v i , j t + 1 = w × v i , j t + c 1 × r a n d × ( p b e s t i , j x i , j t ) + c 2 × r a n d × ( g b e s t j x i , j t ) e l s e r > 0.7 t h e n v i , j t + 1 = w × v i , j t + c 1 × r a n d × ( p b e s t i , j x i , j t ) + c 2 × r a n d × ( g b e s t f , j x i , j t )
Figure 2 presents a flowchart of the solution algorithm.

4. Case Study and Analysis

In this paper, the utility firm purchased demand from user aggregators for a total of 6 h. The demand price announced by the utility firm is shown in Table 1, and the demand to be provided after the bidding is shown in Table 2. The model of the covariance matrix constructed by the user aggregators with the use of the bidding history data of the demand trading market is shown in Table 3 [25]. In addition, the risk-tolerance parameter was set as β = 0.05; the higher β has the lower the revenue of the user aggregators. After predicting the demand price of the utility firm, the covariance matrix is constructed, and the risk-tolerance parameter is set so that the profit maximization of the user aggregators with risk assessment can be calculated.
Two cases with a demand bidding curve for user aggregators are analyzed to assess the feasibility of the proposed algorithm. The two cases are expressed as:
Case 1.
Profit maximization of the user aggregators with the quadratic bidding curve;
Case 2.
Profit maximization of the user aggregators with different bidding curves.
All data were obtained from TPC Power Development Planning [26].
Case 1.
Profit maximization of the user aggregators with the quadratic bidding curve.
In this case, 13 user aggregators were taken as samples using the quadratic bidding curve. Table 4 shows the parameter coefficients of the bidding curve for the user aggregators.
FPSO is used to optimize the amount and price of bidding demand as shown in Table 5. Because the demand prices for users 4 and 5 are relatively high, they have no winning bids. In these 6 h, the demand price provided by the utility firm to the user aggregators is USD 170,192.1, and the risk cost calculated by the Vest covariance matrix and the risk tolerance parameter β = 0.05 is USD 13,078.08. The demand purchase cost paid by the user aggregators to the downstream users is USD 120,741, and the average purchase price per unit is USD 52.78 MW. Therefore, the maximum profit of the user aggregators with risk assessment is USD 34,793.28.
When the demand purchase is a quadratic curve, the demand purchase amount of the user aggregators also affects the profit directly. When the sum of demands in Table 5 is 100%, the cases of 50% and 150% demand are calculated. The profit of the user aggregators is shown in Figure 3. When the risk tolerance parameter β = 0.05, the profit at 50% demand is USD 20,314.10, and the profit at 150% demand is USD 18,363.59. Both of them are less than the profit of 100% demand, which is USD 34,793.28. Therefore, it could be found that the user aggregators must effectively evaluate the purchase cost of downstream users and select the appropriate demand to provide while participating in the purchase from the utility firm.
Case 2.
Profit maximization of the user aggregators with different bidding curves.
The user aggregators consisted of 13 users. The bidding curves are linear, quadratic, and segment linear, as shown in Figure 4, and their detailed data parameters are shown in Table 6. Users 3, 6, 9, and 12 are segmented linearly with forbidden zones, and purchase demand could not be purchased from them between the minimum and maximum non-operational zones.
Table 7 shows the amount and price of bidding demand. In Table 7, the bidding curves of users 4 and 5 were relatively high, while the cost of user 12 became expensive in the later period. In these 6 h, the demand purchase price provided by the utility firm to the user aggregators is USD 170,192.1, and the risk cost calculated by the Vest covariance matrix and the risk tolerance parameter β = 0.05 is USD 13,078.07. The demand purchase cost paid by the user aggregators to the downstream users is USD 93,795.15. Therefore, the maximum profit of the user aggregators with risk assessment is USD 74,673.67. Because Case 1 has a quadratic purchase curve, the profit of Case 2 is higher than that of Case 1.
When the demand purchase curve is different, the demand purchase amount of the user aggregators will also affect the profit directly. When the sum of demand shown in Table 7 was 100%, the cases of 50% and 150% demand were calculated. The profit of the user aggregators is shown in Figure 5. The profit of 50% demand is USD 45,700.71, the profit of 100% demand is USD 74,673.67, and the profit of 150% demand is USD 77,030.91. It was found that the profit of the DR user aggregators at 100% and 150% became gradually flat and could not generate more profit.
Based on the measurement formulas of Case 1 and Case 2, it can be seen that different user demand purchase curves had a great influence on the profit of the user aggregators. Therefore, the proportion of the user aggregators that participate in the demand purchase by the utility firm also needs to be effectively evaluated to calculate the maximum profit value for the user aggregators.

5. Conclusions

This study mainly discussed the demand bidding model of user aggregators and considered the maximum profit after taking profit and risk into consideration. The covariance matrix of demand prices was used to establish the demand price variability and derive the risk model of demand bidding. By using a risk-tolerance parameter β, a new demand bidding model for the user aggregators that takes both risk and profit into consideration was formulated. All constraints, including the load balance and bidding capacity constraints, were considered in this paper. This paper used Feasible Particle Swarm Optimization (FPSO) to perform simulations and analyses. Two cases with different bidding curves of user aggregators were used to simulate the efficiency of the proposed algorithm. We found that the user aggregators can effectively select the appropriate risk and demand during the bidding procedure. It is believed that the proposed algorithms could provide relevant practitioners with a simple and rapid tool to cope with possible uncertain electric price changes in the future and improve the practitioners’ abilities in managing dispatching risks and making cost evaluations, as well as generating capability constraints.

Author Contributions

First author, project idea, related experiences, system model, manuscript—English revision, C.-J.T.; experimental and conducted the simulations, C.-S.T.; project assistance and manuscript—preparation as the corresponding author, M.-T.T. All authors discussed the simulation results and approved the publication. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Science and Technology Council, Taiwan (Grant Nos. MOST 111-2221-E-230-002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Taiwan Power Company. Available online: http://www.taipower.com.tw/content/new_info/ (accessed on 1 April 2022).
  2. Taiwan Power Company. Time-of-Use Rate for Industry. Available online: http://www.taipower.com.tw/electricity (accessed on 1 April 2022).
  3. Glazunova, A. Development of a Day-Ahead Demand Side Management Strategy to Improve the Microgrid Efficiency. IFAC-Papers OnLine 2022, 55, 256–261. [Google Scholar] [CrossRef]
  4. Chung, C.Y. The study of the electricity analysis and demand response strategy of shopping malls. Master’s Thesis, Cheng-Shing University, Kaohsiung, Taiwan, 2021. [Google Scholar]
  5. Pasetti, M.; Rinaldi, S.; Manerba, D. A virtual plant architecture for the demand-side management of smart prosumers. Appl. Sci. 2018, 8, 432. [Google Scholar] [CrossRef] [Green Version]
  6. Duan, Q. A price-based demand response scheduling model in day-ahead electricity market. In Proceedings of the 2016 IEEE Power and Energy Society General Meeting (PESGM), Boston, MA, USA, 17–21 July 2016. [Google Scholar]
  7. Safar, A.; Tahmasebi, M.; Pasupuleti, J. Smart buildings aggregator bidding strategy as a negawatt demand response resources in the spinning reserve electricity market. In Proceedings of the 2019 IEEE PES Innovative Smart Grid Technologies Europe (ISGT-Europe), Bucharest, Romania, 29 September–2 October 2019. [Google Scholar]
  8. Jia, Y.; Mi, Z.; Yu, Y.; Song, Z.; Liu, L.; Sun, X. Purchase bidding strategy for load agent with the incentive-based demand response. IEEE Access 2019, 7, 58626–58637. [Google Scholar] [CrossRef]
  9. Dadashi, M.; Haghifam, S.; Zare, K.; Haghifam, M.R.; Abapour, M. Short-term scheduling of electricity retailers in the presence of Demand Response Aggregators: A two-stage stochastic bi-Level programming approach. Energy 2020, 205, 117926. [Google Scholar] [CrossRef]
  10. Abapour, S.; Mohammadi, B.; Hagh, M.T. Robust bidding strategy for demand response aggregators in electricity market based on game theory. J. Clean. Prod. 2020, 243, 118393. [Google Scholar] [CrossRef]
  11. Liu, X.; Li, Y.; Lin, X.; Guo, J.; Shi, Y.; Shen, Y. Dynamic bidding strategy for a demand response aggregator in the frequency regulation market. Appl. Energy 2022, 314, 118998. [Google Scholar] [CrossRef]
  12. Li, S.; Zhang, L.; Nie, L.; Wang, J. Trading strategy and benefit optimization of load aggregators in integrated energy systems considering integrated demand response: A hierarchical Stackelberg game. Energy 2022, 249, 123678. [Google Scholar] [CrossRef]
  13. Pourramezan, A.; Samadi, M. A novel approach for incorporating incentive-based and price-based demand response programs in long-term generation investment planning. Int. J. Electr. Power Energy Syst. 2022, 142, 108315. [Google Scholar] [CrossRef]
  14. Tang, W.; Yang, H.T. Optimal operation and bidding strategy of a virtual power plant integrated with energy storage systems and elasticity demand response. IEEE Access 2019, 7, 79798–79809. [Google Scholar] [CrossRef]
  15. Dahraie, M.V.; Kermani, H.R.; Moghadda, A.A. Risk-averse probabilistic framework for scheduling of virtual power plants considering demand response and uncertainties. Int. J. Electr. Power Energy Syst. 2020, 121, 106126. [Google Scholar] [CrossRef]
  16. Liang, Z.; Alsafasfeh, Q.; Jin, T.; Pourbabak, H.; Su, W. Risk-constraints optimal energy management for virtual power plants considering correlated demand response. IEEE Trans. Smart Grid 2019, 10, 1577–1587. [Google Scholar] [CrossRef]
  17. Zare, K.; Nojavan, S. Risk-based optimal bidding strategy of generation company in day-ahead electricity market using information gap decision theory. Int. J. Electr. Power Energy Syst. 2013, 147, 83–92. [Google Scholar]
  18. Kim, H.J.; Kim, M.K. Risk-based hybrid energy management with developing bidding strategy and advanced demand response of grid-connected microgrid based on stochastic/information gap decision theory. Int. J. Electr. Power Energy Syst. 2021, 131, 107046. [Google Scholar] [CrossRef]
  19. Tian, M.T.; Yan, S.R.; Tian, X.T.; Nojavan, S.; Jermsittiparsrt, K. Risk and profit-based bidding and offering strategies for pumped hydro storage in the energy market. J. Clean. Prod. 2020, 256, 120715. [Google Scholar] [CrossRef]
  20. Fazlalipour, P.; Ehsan, M.; Mohammadi-Ivatloo, B. Risk-aware stochastic bidding strategy of renewable micro-grids in day-ahead and real-time markets. Energy 2019, 171, 689–700. [Google Scholar] [CrossRef]
  21. Alexander, C. Risk Management and Analysis- Volume 1 Measuring and Modeling Financial Risk; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2000. [Google Scholar]
  22. Taiwan Power Company. The Demand Response Load Management Strategy. Available online: https://www.taipower.com.tw/tc/news2.aspx (accessed on 15 December 2021).
  23. Lucas, A.; Zhang, X. Score-driven exponentially weighted moving averages and Value-at-Risk forecasting. Int. J. Forecast. 2016, 32, 293–302. [Google Scholar] [CrossRef] [Green Version]
  24. Kennedy, J.; Eberhart, R.C. Particle Swarm Optimization. In Proceedings of the IEEE International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  25. Lin, W.M.; Gow, H.J.; Tsai, M.T. Electricity Price Forecasting Using Enhanced Probability Neural Network. Energy Convers. Manag. 2010, 51, 2707–2714. [Google Scholar] [CrossRef]
  26. Taiwan Power Company. The Sustainable Operation of White Paper for Taiwan Power Company; Taiwan Power Company: Taipei, Taiwan, 2021. [Google Scholar]
Figure 1. The auto-tuning scheme of FPSO.
Figure 1. The auto-tuning scheme of FPSO.
Energies 16 00156 g001
Figure 2. Flowchart of the solution algorithm.
Figure 2. Flowchart of the solution algorithm.
Energies 16 00156 g002
Figure 3. Profit diagram of bidding demand in Case 1.
Figure 3. Profit diagram of bidding demand in Case 1.
Energies 16 00156 g003
Figure 4. Bidding curves of user aggregators in Case 2.
Figure 4. Bidding curves of user aggregators in Case 2.
Energies 16 00156 g004
Figure 5. Profit diagram of bidding demand in Case 2.
Figure 5. Profit diagram of bidding demand in Case 2.
Energies 16 00156 g005
Table 1. Demand price (USD/MW).
Table 1. Demand price (USD/MW).
Hour101113141516
Price66.60978.37285.69878.27372.07263.342
Table 2. Demand provided of user aggregators.
Table 2. Demand provided of user aggregators.
Hour101113141516
DR(MW)352.81364.18417.69396.41387.3795370.56
Table 3. Covariance matrix of electricity prices.
Table 3. Covariance matrix of electricity prices.
Hour101113141516
100.0036−0.00180.0054−0.0036−0.00180.0018
11−0.00180.0009−0.00270.00180.0009−0.0009
130.0054−0.00270.0081−0.0054−0.00270.0027
14−0.00360.0018−0.00540.00360.0018−0.0018
15−0.00180.0009−0.00270.00180.0009−0.0009
160.0018−0.00090.0027−0.0018−0.00090.0009
Table 4. Bidding curve of user aggregators.
Table 4. Bidding curve of user aggregators.
AggregatorMax. (MW)Min. (MW)abC
117.130.6933.659.4705
228.550.94240.936.903
34550.35740.1528.771
44550.60564.50
575100.42162.591.34
675100.70845.75172.83
782.5150.31339.8564.783
882.5300.29833.1578.596
982.5300.27735.580.132
1022.540.5212416.65105.51
1128.550.1632.1522.292
123050.0144.7510.787
1316.531.6129.430.745
Table 5. Bidding demand amount and price in Case 1.
Table 5. Bidding demand amount and price in Case 1.
UnitDemand Amount (MW/hour)Total (MW)Total Purchase Price (USD)Purchase Price Per Unit (USD/MW)
101113141516
117.0917.1017.1017.0517.1017.07102.524715.0445.99
216.3415.6318.0417.1118.2214.83100.175902.0458.92
342.3745.0045.0044.9045.0045.00267.2715,155.6756.71
4000000000
5000000000
613.6519.8023.1921.3222.7116.94117.618097.6368.85
742.5742.9759.5054.1546.9347.21293.3316,635.6356.71
863.9671.1682.5068.9772.7165.05424.3523,549.4055.49
962.9259.4375.1776.3169.7969.21412.8223,064.6555.87
1022.5022.5022.5022.5022.3922.50134.894459.8333.06
1128.5028.5028.5028.5028.5028.34170.846404.5737.49
1230.0029.9430.0029.9230.0030.00179.868167.3445.41
1312.1011.9716.5015.2814.6513.8784.354599.6054.53
Table 6. Bidding curves of the different user aggregators.
Table 6. Bidding curves of the different user aggregators.
Aggregator1234567
Max.17.128.54545757582.5
Min.3555101015
A00.942000.42100
b133.6540.940.1564.562.545.7539.85
C9.470536.90328.77172.28291.34172.8364.783
b20048.180054.90
Not Operate Min. 23 40
Not Operate Max. 27 45
Aggregator 8910111213
Max. 82.582.522.528.53016.5
Min. 30304553
A 0.298000.1600
b1 33.1535.516.6532.1544.7529.4
C 78.59680.132105.5122.29210.78730.745
b2 042.60053.70
Not Operate Min. 54 16
Not Operate Max. 59 19
The operating curve of aggregators 1, 4, 7, 10, and 13 is a linear function. The operating curve of aggregators 2, 5, 8, and 11 is a quadratic function. The operating curve of aggregators 3, 6, 9, and 12 is a segment linear function.
Table 7. Bidding demand amount and price in Case 2.
Table 7. Bidding demand amount and price in Case 2.
UnitThe Demand Amount (MW/Hour)Total
(MW)
Total Purchase
Price (USD)
Purchase Price
Per Unit (USD/MW)
101113141516
117.1015.6517.1017.0817.0917.10101.123459.4334.21
20.000.0014.770.000.000.0014.77846.8057.32
345.0044.4045.0045.0045.0044.97269.3710,934.0240.59
4000000000
5000000000
638.8375.0075.0075.0075.0064.84403.6718,090.7044.82
781.1382.5082.5082.5082.5082.50493.6320,059.6640.64
854.7343.3652.0846.8540.8634.65272.5213,276.1148.72
947.7235.7953.8753.9253.3552.94297.5912,120.8940.73
1022.5022.5022.5022.5022.5022.50135.002880.8121.34
1128.5028.2928.5028.5028.5028.50170.796402.6537.49
120010.188.156.695.5430.572260.6573.96
1316.5016.5016.5016.5016.5016.4798.973094.2431.26
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

Tien, C.-J.; Tu, C.-S.; Tsai, M.-T. Risk Assessment of User Aggregators in Demand Bidding Markets. Energies 2023, 16, 156. https://doi.org/10.3390/en16010156

AMA Style

Tien C-J, Tu C-S, Tsai M-T. Risk Assessment of User Aggregators in Demand Bidding Markets. Energies. 2023; 16(1):156. https://doi.org/10.3390/en16010156

Chicago/Turabian Style

Tien, Ching-Jui, Chia-Sheng Tu, and Ming-Tang Tsai. 2023. "Risk Assessment of User Aggregators in Demand Bidding Markets" Energies 16, no. 1: 156. https://doi.org/10.3390/en16010156

APA Style

Tien, C. -J., Tu, C. -S., & Tsai, M. -T. (2023). Risk Assessment of User Aggregators in Demand Bidding Markets. Energies, 16(1), 156. https://doi.org/10.3390/en16010156

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