1. Introduction
With escalating energy demand and advancements in the energy sector, carbon dioxide emissions have surged annually. It is crucial to advance decarbonization in power generation, develop renewable energy source (RESs) technologies, foster the transition towards RESs, and build a new power system predominantly reliant on RESs for achieving China’s carbon peaking and carbon neutrality goals [
1].
The conventional power system depends on substantial capacity of thermal units and the secure, stable operation of backup system for wind and solar power. However, renewable resources like wind and solar power are highly susceptible to climatic, seasonal, and temperature variations, introducing significant uncertainty and volatility. As the integration of these resources into the grid increases, the system faces amplified uncertainty at both the generation and consumption ends, complicating operation and challenging its security, stability, and dynamic balance [
2]. Thus, it is of paramount importance to study the optimal scheduling of systems with a high proportion of RESs. Demand response (DR) can effectively integrate flexible loads, reducing the use of high-cost thermal units [
3]; thereby, reducing the load in peak periods can avoid unit capacity expansions and ease bottlenecks in the power system [
4]. DR strategies based on demand-side resource management have emerged as crucial for ensuring reliable power supply and the safe and stable operation of new power systems due to their flexible regulation, high potential for adjustment, minimal grid impact, and low costs [
5]. In the current study, researchers generally categorize DR into price-based DR (PDR) and incentive-based DR (IDR) [
6,
7]. PDR guides electricity consumers to optimize their own load profiles through the setting of price mechanisms like time-of-use (TOU) tariffs and categorized electricity price. IDR includes modes of operation such as interruptible load, incentive coupon issuance, and direct load control, which do not involve tariffs and customers will not bear direct benefit losses [
8]. In Ref. [
9], a solid estimation of consumer demand and its relationship with TOU tariffs is offered. Ref. [
10] proposes a dynamic tariff scheme based on linear regression, which guides residential customers to participate in DR and addresses the optimization challenge of DR using a particle swarm optimization algorithm. Ref. [
11] introduces an incentive-based DR (IDR) strategy, employing two-stage optimization approach to minimize the cost of flexibility management and improve the peak-to-average ratio of the community microgrid load profile. Ref. [
12] explores an optimization problem for operating energy hubs (EHs), aiming to enhance customer benefit and reduce the overall operational cost of the EH system, in the context of both PDR and IDR programs. However, the two DR mechanisms mentioned above currently suffer from the following shortcomings: PDR is an involuntary regulation method and is affected by differences in customers’ perceptions of the DR benefits, which is not customer-friendly, and its large-scale rollout can lead to the over-regulation or rebound of loads [
13,
14]. IDR widely adopts customer baseline load (CBL) to evaluate customer contributions, but in large-scale DR, the calculation of CBL for large number of customers faces enormous data calculation, storage and communication pressure, and the accuracy and authenticity of CBL cannot be guaranteed [
15].
For these reasons, some scholars proposed the strategy of customer directrix load (CDL), a desired load profile that offsets power system fluctuations effectively [
16,
17]. Customers can make their load profiles as close as possible to the CDL in order to reap the benefits. CDL-based DR is simple and easy to implement, has the characteristics of widespread promotion and normalized implementation, and can well support the construction of new power systems, so it has received extensive attention from scholars. Ref. [
18] proposes a calculation model of nodal customer directrix load (NCDL) with robust security, which can guarantee that the power flow across all lines remains within established limits for any real response outcomes within the designated DR range. The above literature takes the associated network constraints and the uncertainty of flexible loads into account, however has not yet considered the uncertainty of RES output, which is also a factor that should not be ignored in scheduling optimization. The scenario optimization method is one of the main approaches to study the uncertainty of RES, and uncertainty is considered by generating sets of scenarios through Monte Carlo simulation. This method can accurately characterize the random variables. While the current CDL-based DR optimal scheduling model does not consider hydro units, the modeling of joint optimal scheduling problems of cascade hydropower plants has emerged as a challenging and significant problem in the field of optimal scheduling. Ref. [
19] introduces a scheduling strategy for a thermal-hydro-wind-solar system that incorporates the impacts of water head and discharges within a single function. This strategy enables the transformation of the head-dependent hydropower conversion function into a piecewise linear function, which can be efficiently solved using a mixed-integer programming (MIP) model.
To address the above challenges, we propose a multi-sources stochastic coordinated scheduling model based on CDL-based DR. Compared with traditional methods, the model we propose integrates the uncertainty of RES through various scenarios generated using Monte Carlo simulation. By examining the complementary traits between wind and solar generation, we can establish the fundamental functions of both. The scenarios can be generated by using the Copula function, which is used to describe the complementary traits of wind and solar output. A fast scenario reduction technique based on improved k-means clustering and the SBR algorithm [
20] is employed to balance precision with computational efficiency. The dynamic characteristics of hydro units, which are not currently considered, are also included to improve energy efficiency and stable supply and reduce energy volatility through multi-energy complementarity. Finally, case simulations were carried out by means of a modified IEEE 6-bus power system arithmetic example, and the results of the case study indicate that the proposed model can effectively improve system flexibility, reduce the curtailment of RES in a high-proportion RES power system, and obtain the intermediate scheduling solutions in the base case with the uncertainty of RESs.
The main contributions are presented as follows:
- (1)
Compared with traditional stochastic optimization methods where the objective function is to minimize the cost of system operation in all scenarios, the proposed model aims to obtain a dispatch solution with a stable operating cost in the base case. This dispatch solution allows for the power system to operate under the forecasted base case and can safely redispatch all units in response to real-time fluctuations in RES output.
- (2)
We introduce the dynamic characteristics of cascade hydro units and CDL-based DR in day-ahead stochastic scheduling model with RES uncertainty for the first time. The complementary characteristics of wind and solar are considered and modeled using the t-Copula function. These all aim to improve energy efficiency, ensure stable supply, and reduce energy fluctuations through multi-sources complementary and coordinated scheduling.
We concisely encapsulate this paper’s structure: after this introduction,
Section 2 introduces the concept of CDL, defines the response effect metrics and incentives for CDL.
Section 3 introduces our methodology, focusing on the novel day-ahead stochastic scheduling model incorporating CDL-based DR. In
Section 4, we detail the formulation of our model, elucidating on the use of Copula theory and Monte Carlo simulations for addressing uncertainties in RESs.
Section 5 applies the proposed model to a modified IEEE-6 bus system, showcasing the practical application and effectiveness of our approach through case studies. The discussion in
Section 6 extends the analysis of our findings, contemplating their broader implications and suggesting avenues for future research.
Section 7 synthesizes our contributions, emphasizing the effectiveness of our novel stochastic scheduling method in enhancing RES consumption and robustness of system operation.
3. Operation Decision Model of Multi-Sources Power System with CDL-Based DR
In this section, we will introduce the calculation model of CDL and the scheduling model of a multi-sources power system. In the first stage, the DR center obtains the CDL based on the day-ahead forecasted RES output, load demand, and other relevant constraints. At this point, the flexible load
participates in the decision as a variable. Then, the decision-making process of customers who participate in DR is simulated by using the customers decision model introduced in
Section 2.2. The actual load profile
after DR can be obtained in this customers decision model. In the second stage, based on the simulation results of the first stage, the scheduling model of multi-sources power system is solved. At this point, the flexible load
participates in the decision as a constant with the same value as
.
3.1. Deterministic Scheduling Model
3.1.1. Objective Function
The purpose of the water-thermal-wind-solar scheduling model is to minimize the system costs (including fuel cost and unit start-stop cost) while considering the cost associated with curtailed RES output and loss of load and to promote the absorption of RES. It is noted that DR costs are reflected ex-post, meaning that the response targets are set decoupled from DR costs. Therefore, DR costs do not need to be considered in the objective function as shown:
where
is the start-up and shutdown cost of the thermal power units;
is the operating cost of the thermal power units;
is the cost of RES abandonment; and
is the loss of load cost;
where
is the number of flexible generators;
and
are the start-up and shutdown costs of thermal unit
at
time;
,
,
are the cost coefficient of thermal unit
;
is active power thermal unit
at
time;
is the penalty coefficient of RES abandonment;
is the maximum power of RES generation output at
time;
is the power of RES generation output at
time;
is the penalty coefficient of loss of load; and
is the power of loss of load at
time.
3.1.2. Constraints
The constraints of the model mainly consist of system constraints, thermal power unit constraints, cascade hydro unit constraints, wind power and solar output constraints, and load collinear constraints [
21]. The system constraints mainly include the load balancing constraint (13) and the direct current (DC) power flow constraint (14)
where
Pht is the output of hydro unit
h at
t time;
Pdt is the demand value of electric load
d at
t time;
SF is the transfer matrix;
PLmax is the maximum capacity matrix on each transfer line;
are the scheduling vectors of thermal power, hydro power, wind power, and solar power, and
are the correlation matrix of buses with thermal power, hydro power, wind power, and solar power;
is the dispatch vector of bus load; and
is the incidence matrix of bus load.
The thermal power unit constraints mainly include the unit capacity constraint (15), minimum start/stop time constraints (16) and (17), start/stop cost constraints (18) and (19), and upward and downward climb constraints (20) and (21):
where
is the commitment state of unit
i at
t time,
are the minimum on/off time of unit
i,
are the timers for the on/off of unit
i at
t time, which record the times of startup and shutdown of the unit,
sui,
sdi are the startup/shutdown costs of unit
i, and
are the ramp-up and ramp-down rates of unit
i.
The hydro unit constraints between upstream and downstream reservoirs need to be considered for the cascade hydro unit, and its main constraints are:
where
is the initial and final storage capacity of the hydro unit
h;
Qht is the reservoir discharge volume of the hydro unit
h at the moment of
t;
Vht is the storage capacity of the hydro unit
h at
t time;
is the minimal reservoir discharge volume of the hydro unit
h;
is the maximal reservoir discharge volume of the hydro unit
h;
is the minimum/maximum storage capacity of the hydro unit
h, and
rht is the natural incoming water volume of the hydro unit
h at
t time.
Equations (26) and (27) are the water head function of the cascade hydro unit, which represents the hydropower conversion relationship.
where
Hht is the water head level of the hydro unit
h at time
t, the value of which is related to the physical dimensions of the terraced hydroelectric power plant,
h0,h, α
h is the physical constant with respect to the hydro unit
h, and
ηh is the water-to-power conversion factor.
The RES constraints ensure that the dispatched output of the RESs will not exceed the predicted value at time:
3.1.3. Piecewise Linearization of the Hydropower Conversion Function
Ref. [
19] introduced extra integer variables to convert the nonlinear water-to-power conversion function into a piecewise linear function. Ref. [
21] utilized the heuristics to convert the transformation profile into piecewise linear functions. The previous function is divided into a lattice where each element is divided into two triangles. Therefore, the hydropower conversion function can be expressed in terms of piecewise linearization, then it can be incorporated into the MIP model.
As noted above, we used a
grid to divide
and
into subintervals
and
where
,
. Therefore, Equation (27) can be expressed as (30)–(34):
Every grid is divided into two triangles in the upper left and lower right corners. represents the location index of the upper left triangles. represents the location index of the lower right triangles.
3.1.4. Abstract Formulation
The model we proposed can be expressed in an abstract form to ease the introduction of the stochastic model, as shown in (35)–(38):
where
is a binary vector representing start-stop states, start-stop actions, and auxiliary variables of the linearized hydropower conversion function.
is a continuous vector representing the scheduling decisions for each energy source, and
are abstract matrices and vectors related to the cost and constraint coefficients. Constraint (36) ensures that
is a binary vector. Constraint (37) represents the constraints of binary variables, and Constraint (38) represents the conditions of a system with both binary and continuous variables. Specifically, Equation (35) represents the objective function of the model, i.e., (8)–(12).
represents the start-up and shutdown cost of thermal power units i.e., (9).
represents all costs related to thermal, hydro, wind and solar generation, including all penalty costs, i.e., (10)–(12). Equation (36) represents that
only takes values of 0 or 1. Equation (37) represents constraints related only to the start/stop state of units, i.e., (16)–(19). Equation (38) represents rest constraints related to outputs of various units.
3.2. Stochastic Scheduling Model
The deterministic scheduling model only considers accurate prediction information for scheduling. However, in practice, the DR center prefers to take into account as much uncertainty as possible a day ahead, based on the scenarios, to find a relatively stable operation cost. This study introduces a stochastic scheduling model for a thermal-hydro-wind-solar system that takes into account the uncertainties of RESs as shown in Equations (39)–(44). The design allows the system to function based on a base-case scenario using projected data, with the capability to securely re-dispatch thermal and hydro units to address the real-time variabilities of RESs.
Given the unpredictable nature of wind and photovoltaic generation due to their intermittent characteristics, operators would favor minimal fluctuations in operating costs across various realizations of renewable generation. The objective function (39) consists of the cost of base-case operations and the anticipated variance between other potential scenarios and the base case.
where
is the base-case scheduling decision;
is the scheduling decision in other potential scenarios;
represents the serial number of typical scenarios;
is the weight of the scenario
, and the value of
is set as the probability of occurrence of typical scenarios calculated at
Section 4.2. Scenario Reduction.
represents the base-case costs, i.e., deterministic scheduling model.
represents the absolute value of difference in operating costs excluding start/stop costs between all typical scenarios and the base case.
Equations (40)–(44) represent operational constraints. Specifically, flexible resources like adaptable generators with rapid ramping capabilities are reallocated to balance the electric load. The dispatches in these scenarios are further interconnected through constraint (44), which regulates the corrective ramping capacity of the generating units [
22,
23].
where
are abstract matrices and vectors of cost and constraint coefficient correlations. Specifically, Equation (40) represents that
only takes values of 0 or 1. Equation (41) represents constraints related only to the start/stop state of units, i.e., (16)–(19). Equation (42) represents rest constraints related to outputs of various units in base case. Equation (43) represents rest constraints related to outputs of various units in all typical scenarios. Equation (44) represents the constraints that can further couple the dispatch solution in the base case and typical scenarios.
The absolute value function in Equation (39) can be solved linearly by the big-M method of operational optimization. We have:
where
, a binary value vector, indicates the positive/negative components of
.
represents a sufficiently large value.
Based on the above models, the CDL that can consume as much RESs as possible at the system level can be obtained, and the DR center only needs to publish this profile to inform the customers of the load adjustment target.
Figure 2 shows the CDL profiles derived from both deterministic and stochastic scenarios. It is evident that the CDL profile generated under deterministic conditions exhibits a larger peak-to-valley difference compared to its stochastic counterpart, thus more effectively guiding customers towards achieving peak shaving objectives. However, the deterministic approach yields less robust outcomes. In contrast, the stochastic scenario accounts for uncertainties, necessitating additional committed units to maintain system security. This method results in more robust scheduling decisions that avoid excessive conservatism.
5. Case Studies
The effectiveness of the proposed day-ahead stochastic model was tested in a modified 6-bus power system [
19]. The system node wiring diagram is shown in
Figure 4. There are three conventional flexible thermal power units G1, G2, G3, two hydro units H1, H2, three loads L1, L2, L3, one wind power unit W, and one solar power unit PV.
The detailed parameters of these generators are shown in
Table 1,
Table 2,
Table 3 and
Table 4, with data for the water-lined generators. The historical data and forecasted values of load, wind, and solar power generation are derived from Ref. [
27]. The forecasted values of electrical load, wind, and solar generation are shown in
Figure 5.
is set to 100
, and
is set to 1000
. All case studies were solved using Gurobi 9.5.2 on a personal computer.
5.1. Deterministic Scheduling Case
5.1.1. Case 1
In this case, the dispatch model for the multi-source power system without DR can only dispatch for deterministic forecast values of electrical load and wind and solar generation.
Table 5 shows the impact of RES penetration on the associated dispatch costs in Case 1. The total costs consist of the loss of the load penalty, curtailed RES output costs, and power generation costs. With RES penetration increasing, total costs decrease and then increase mainly due to the loss of the load penalty and curtailed RES output cost. The load loss penalty decreases first but no longer decreases when the penetration rate is 0.4, and the curtailed RES output costs rise rapidly. The cost of thermal power generation has been decreasing, a trend attributed to the rising share of RES installations. This increase in the proportion of renewable sources reduces the generation costs for thermal power plants. When the RES penetration reaches 0.4, the limit of load balancing is reached. In the absence of DR, increasing the penetration of RES at this time will only increase curtailed RES output and reduce the output of thermal units without contributing to the source-load balance of the system and reducing loss of load.
Loss of load and curtailed RES under different RES penetration rates are shown in the
Figure 6. The rate of loss of load decreased from 7.42% to 2.97% and then did not change again. Curtailed RES output rate continued to increase from 46.35% to 78.15%. Obviously, at low permeabilities, like 0.1, the curtailed RES output rate reached 46.35%. This is a very high value. These results show that in a new power system, the reduction in installed capacity of conventional units leads to a huge challenge in the source-load balance of the system and the absorption of RESs.
5.1.2. Case 2
The rate of flexible load participation in DR was set at 0.3. The influence of RES penetration on scheduling outcomes in case 1 is presented in
Table 6. We have examined three scenarios where flexible loads can be seamlessly adjusted to the optimal load profile, i.e., CDL. The results compared to case 1 show that CDL-based DR can effectively mitigate loss of load and reduce the curtailed RES output. When RES penetration is 0.1, the rate of curtailed RES output reduces from 46.35% in case 1 to 0%. Additionally, when RES penetration is 0.6, the rate of curtailed RES output reduces from 78.15% in case 1 to 54.43% because here there is a large RES output, but the load is small in comparison. Although it increases the DR costs, the total costs and the thermal power generation costs are much lower.
5.2. Stochastic Scheduling Case
The complementary traits of wind-solar generations were modeled by using the Copula theory proposed in
Section 4.1 to obtain 10,000 scenarios of wind-solar generations. Subsequently, we employed the fast-forward method presented in
Section 4.2 to obtain reduced data sets of scenarios, providing 5 scenarios as a reasonable approximation of the 10,000 scenarios, as shown in
Figure 7.
Table 7 shows the value of the weighting coefficients
in the objective function (39) obtained in scenario reduction.
In this section, we consider the uncertainties associated with wind and solar power generation, and the following scenarios are explored.
5.2.1. Case 3
Case 3 considers stochastic scheduling based on case 2. The rate of flexible load participating in DR was set at 0.3 and the RES penetration was set at 0.3. The customers adjust their load profiles as much as possible to improve the similarity metric
in order to obtain more incentives. In case 3, the customer decision-making process was modeled by using the response evaluation model mentioned in
Section 2.2. As described in
Section 3, we obtained the value of
and
in first stage. Then we calculated the value of
. We can constrain the variable
by setting the value of
as described in (6). For ease of discussion, we incremented the value of
by 0.05. Thus, the value of the right-hand side of inequality (6) was increased by 0.05 from 0.75. Specially, because in practice it is impossible for
to take the value 1, the last band was set to 0.99 instead of 1. Therefore, the similarity metric
was divided into six bands from 0.75 to 0.99.
Figure 8 shows the shaping of the CDL profiles and actual customers’ load profiles at different response effects: during time periods 1 to 7, the units take on less load, and the CDL guides the customers to shift as much load as possible into that time period in the process of improving similarity
, boosting the load and thus filling in the valleys. Meanwhile, during time periods 9 to 16 and 17 to 23, the units take on more load, and the CDL guides the customers to shift as much load as possible out of that time period in the process of improving similarity metric
, reducing the load and thus achieving peak shaving.
The various costs at different values of
are shown the
Table 8. Building on the observation from
Table 8, the diminishing overall system costs with the increase in the similarity metric
underscores the efficiency of DR mechanisms. As
ascends, indicating higher customer participation in DR, it not only compensates for the DR costs incurred but also reduces operational and imbalance costs significantly. This trend demonstrates the value of aligning consumer behavior with system needs, leveraging flexible loads to enhance grid stability and integrate RES more effectively. This analysis suggests that a proactive DR strategy, supported by incentives and robust forecasting, can lead to a sustainable, cost-efficient energy system by mitigating the challenges posed by renewable intermittency and demand fluctuations. Therefore, the system can not only reduce the traditional generating units and use more RES units, which is conducive to the protection of the environment and reduce the use of fossil energy, but also obtain great economic benefits.
5.2.2. Case 4
In this case, with the aim to discuss the influence of forecasting errors in wind and solar generation, we reduced forecasting errors on the basis of case 3 by using the scenario modification method proposed in Ref. [
19]. The error rate was set within 0.05. The forecasted wind and solar generation values and scenarios are displayed in
Figure 9. In comparison with
Figure 7, it can be observed that larger forecasting errors result in greater deviations in wind and solar generation, as well as higher levels of uncertainties. Forecasted errors in wind and solar generation were set to less than 0.05, and we discuss stochastic dispatch when the similarity metric
was set to 0.9.
Figure 10 shows the results of unit commitment (UC) in case 3 and case 4. It should be noted that throughout the day hydro units H1, H2, along with thermal units G1, G3, remain in continuous operation. The UC solutions we obtained in case 3 differ from those in case 4. Due to the higher level of uncertainties in case 3, additional units must be activated to ensure the security of the system. In particular, there is a significant fluctuation in wind and solar generation in hour 12, during which the system shows insufficient ramping capacity. Therefore, thermal unit G2 provides enough ramping capabilities in case 3 by operating in additional hour 12, and the operation cost in case3 is USD 214,223.60, which is 11.37% more than that of case 4.
Table 9 shows the operation costs in case 4, and it is noted that the total costs, loss of load penalty, and thermal power generation costs are lower than those of case 3.
The base-case dispatches of H1 and H2 are increased from 3652.70 MWh in case 3 to 3751.45 MWh in case 4, and the dispatch of G3 was decreased from 1073.65 MWh in case 3 to 876.96 MWh in case 4. To ensure the system security under increased uncertainty and larger deviations in Case 3, an additional spinning reserve must be allocated to offset the fluctuations in wind and solar power generation from hydroelectric units. Consequently, less expensive hydro units operate at a suboptimal dispatch level in Case 3, but hydro units in case 4 are operating at a better dispatch level. As a result, the total costs are much lower in case 4.
6. Discussion
The simulation results in each case reveal the following insights: firstly, hydro units, as an effective regulatory resource, enhance the thermal-hydro-wind-solar system’s capacity to absorb RESs. However, hydro generation is influenced by hydrological conditions, seasonal variations, and other factors. Secondly, compared to scheduling models without DR, this study’s approach encourages consumers to adjust their electricity usage proactively, engaging more actively in DR. With increased participation in CDL-based DR, the system experiences significant peak shaving and valley filling effects, thereby enhancing the system’s reliability and economic efficiency. This, in turn, substantially facilitates the integration of a high proportion of new energy sources. Thirdly, the uncertainty of RES and the prediction errors in their generation affect operational costs. Larger prediction errors lead to increased levels of uncertainty, necessitating that hydro units provide more spinning reserve to compensate for the variability in wind and solar power generation. Additionally, more thermal units are required to ensure system reliability and stability. In summary, the proposed stochastic model offers a practical approach for identifying an optimal base-case dispatch strategy, thereby minimizing cost fluctuations amidst uncertainties.
As noted in Ref. [
28], the expansion of the number of scenarios considerably increases the computational load.
Table 10 shows the average computation time for all parameters in cases 1–4. It is observed that when not considering uncertainties in RESs, the computation time in case 1 is 4.09 s and that in case 2 is 4.78 s. When considering uncertainties in RESs, the computation time is increased significantly in case 3 and case 4. The computation time in case 3 with higher-level uncertainties is 458.59 s. It is about three times as long as in case 4. The reason is that more reserves are required to be turned on to balance real-time fluctuations in RES output by turning on additional units. Therefore, it needs longer time to obtain the solution in case 3.
7. Conclusions
We present a day-ahead stochastic scheduling method for day-ahead energy systems with CDL-based DR. The complementary traits between wind and solar power generation are modeled through Copula theory, and the Monte Carlo method is applied to simulate the wind and photovoltaic output errors to establish a set of stochastic scenarios for RES output. The objective function aims at minimizing the costs associated with RES integration, DR and load shedding, incorporating the uncertainties of cascaded hydroelectric units and RES to develop a day-ahead economic dispatch model that leverages multiple energy sources complementarily. Our model aims to obtain a dispatch solution with stable operating cost in the base case. This dispatch solution allows for the power system to operate under the forecasted base case and can safely redispatch all units in response to re-al-time fluctuations in RES output. Additionally, the model includes safety constraints to ensure more secure dispatch outcomes and effectively prevent flow overloading.
Our study demonstrates the potential of a day-ahead stochastic scheduling model with CDL to enhance the flexibility and efficiency of power systems. However, practical implementation faces several challenges: The deployment of this model requires advanced infrastructure capable of handling high volumes of data from diverse energy sources. Effective implementation relies on sophisticated communication protocols to ensure seamless interaction between the power system’s components. The protocols must ensure data integrity, security, and timeliness to facilitate the dynamic scheduling process. The regulatory environment plays a crucial role in the model’s implementation. Successful implementation requires coordination among all stakeholders, including power generators, grid operators, consumers, and regulatory bodies. Overcoming these challenges requires concerted efforts from all stakeholders and researchers.