Power System Dispatch Based on Improved Scenario Division with Physical and Data-Driven Features

: In power systems with high penetration of renewable energy, traditional physical model-based optimal dispatch methods suffer from modeling difﬁculties and poor adaptability, while data-driven dispatch methods, represented by reinforcement learning, have the advantage of fast decision making and reﬂecting long-term beneﬁts. However, the performances of data-driven methods are much limited by the problem of distribution shift under insufﬁcient power system scenario samples in the training. To address this issue, this paper proposes an improved scenario division method by integrating the power system’s key physical features and the data-driven variational autoencoder (VAE)-generated features. Next, based on the scenario division results, a multi-scenario data-driven dispatch model is established. The effectiveness of the proposed method is veriﬁed by a simulation conducted on a real power system model in a province of China.


Introduction
Power systems are experiencing increasing penetration of renewable energy.The strong uncertainty of renewables makes the power system dispatching problem complex.So, the research of optimal dispatching methods for power systems with high renewable energy penetration is of great importance [1][2][3][4][5].
Traditional optimization-based dispatching methods suffer from modeling difficulties and calculating complexity, which make it difficult to meet the real-time decisions of complex power systems [6][7][8].In contrast, data-driven dispatching methods represented by reinforcement learning have the advantages of strong adaptability to environmental uncertainty, high real-time decision efficiency, and a balance between short-term and long-term benefits [9][10][11][12][13][14].
Data-driven dispatch methods can be further divided into single-scene and multi-scene dispatch.Single-scene data-driven decisions are relatively simple and widely used, but the problem of dispatching operating scenario data differences is usually ignored.Ref. [15] proposed microgrid dispatch strategies based on asynchronous dominant actor-critic (A3C) algorithms.Ref. [16] proposed a robust economic dispatch strategy for virtual power plants and ref. [17] proposed an energy storage system control strategy for a combined wind-solar storage system scenario based on deep reinforcement learning (DRL).Ref. [18] proposed an optimal dispatch method for integrated energy systems by using the proximal policy optimization (PPO) algorithm for considering security constraints.
The above data-driven dispatch methods in a single scenario usually suffer from the problem of a limited number of scenarios in the training samples, which leads to a large difference between the actual operating scenarios of the power system and the 1.
To address the limitations of data-driven dimensionality reduction methods, a datadriven dimensionality reduction method combined with power system physical indicators is proposed to construct system operational features.The extracted features are more interpretable and comprehensive than the traditional data-driven scenario division method.2.
Based on the above integrated physical and data-driven features, the Gaussian Mixture Model (GMM)-based power system scenario division method is proposed, which can effectively identify the differences in the features of each operation section.

3.
Based on scenario division, a multi-scenario data-driven dispatching model using a proximal strategy optimization algorithm is proposed for the rapid decision making of power system dispatch.
The remainder of the paper is organized as follows.Section 2 proposes the framework of multi-scenario data-driven dispatch.Section 3 constructs the power system's operational features.Section 4 proposes scenario division methods.Section 5 proposes a data-driven dispatch model of the power system.Section 6 validates the methods by simulations.Finally, Section 7 concludes the paper.

Framework of Power System Dispatch with Scenario Division
The framework of the proposed power system dispatch with scenario division is shown in Figure 1.
Firstly, the power system historical operation sections are obtained, which cover the operating data of thermal and renewable generations, loads, and the grid.These operation sections represent the operating states of the grid at a certain historical instant [24].Next, the data-driven features are extracted using a variational autoencoder (VAE) based on probability distribution, and the power system's key physical features are selected considering the generation-load features.Then, the GMM-based clustering method is used to divide the operation sections into scenarios, and the division is verified based on the Silhouette Coefficients (SC).Finally, in the offline stage, a reinforcement learning-based dispatch model is established for each divided scenario, and the scenario-specific agent is trained based on the PPO algorithm.In the online stage, the scenario-specific agent is used to realize fast online decision making for power system dispatch.Firstly, the power system historical operation sections are obtained, which cover the operating data of thermal and renewable generations, loads, and the grid.These operation sections represent the operating states of the grid at a certain historical instant [24].Next, the data-driven features are extracted using a variational autoencoder (VAE) based on probability distribution, and the power system's key physical features are selected considering the generation-load features.Then, the GMM-based clustering method is used to divide the operation sections into scenarios, and the division is verified based on the Silhouette Coefficients (SC).Finally, in the offline stage, a reinforcement learning-based dispatch model is established for each divided scenario, and the scenario-specific agent is trained based on the PPO algorithm.In the online stage, the scenario-specific agent is used to realize fast online decision making for power system dispatch.

Power System Operational Features
The power system operational data have a massive scale, large dimension, and multiple spatiotemporal scales.A feature extraction method integrating the data-driven and physical features is proposed and shown in Figure 2, which will provide a strong basis for the operational scenario division later.As shown in Figure 2, the system level features reflect the overall operating trend of the power system, and the nodal level features reflect the overall load and generation level of the power system.Therefore, at the system level and nodal level, physical feature selection and the construction of VAE-based feature extraction models are used for key feature mining in this paper.

Power System Operational Features
The power system operational data have a massive scale, large dimension, and multiple spatiotemporal scales.A feature extraction method integrating the data-driven and physical features is proposed and shown in Figure 2, which will provide a strong basis for the operational scenario division later.As shown in Figure 2, the system level features reflect the overall operating trend of the power system, and the nodal level features reflect the overall load and generation level of the power system.Therefore, at the system level and nodal level, physical feature selection and the construction of VAE-based feature extraction models are used for key feature mining in this paper.Firstly, the power system historical operation sections are obtained, which cover the operating data of thermal and renewable generations, loads, and the grid.These operation sections represent the operating states of the grid at a certain historical instant [24].Next, the data-driven features are extracted using a variational autoencoder (VAE) based on probability distribution, and the power system's key physical features are selected considering the generation-load features.Then, the GMM-based clustering method is used to divide the operation sections into scenarios, and the division is verified based on the Silhouette Coefficients (SC).Finally, in the offline stage, a reinforcement learning-based dispatch model is established for each divided scenario, and the scenario-specific agent is trained based on the PPO algorithm.In the online stage, the scenario-specific agent is used to realize fast online decision making for power system dispatch.

Power System Operational Features
The power system operational data have a massive scale, large dimension, and multiple spatiotemporal scales.A feature extraction method integrating the data-driven and physical features is proposed and shown in Figure 2, which will provide a strong basis for the operational scenario division later.As shown in Figure 2, the system level features reflect the overall operating trend of the power system, and the nodal level features reflect the overall load and generation level of the power system.Therefore, at the system level and nodal level, physical feature selection and the construction of VAE-based feature extraction models are used for key feature mining in this paper.Among the system level, peak load and average load present the trend of load, and renewable average, renewable peak-valley difference, and renewable peak present the trend of renewables.Among the nodal level, load levels reflect the operating condition of the power system, and high load levels indicate a high demand for electricity in the system.Renewable outputs reflect the system's renewable processing, and high renewable outputs indicate sufficient renewable energy resources.

Data-Driven Feature Extraction Using VAE
By maximizing information retention and minimizing reconstruction loss of the decoder, VAE can reduce the hidden layer feature space dimension and extract features.
Energies 2023, 16, 7520 4 of 14 Thus, the negative impact of generation-load uncertainty on subsequent power system dispatch operations can be effectively eliminated.The structural diagram of the feature extraction model based on VAE is shown in Figure 3.As shown in Figure 3, the mean and variance are extracted by the decoder to fit the data distribution.Then, the dataset is reconstructed through a decoder.During the network iteration process, the network is iteratively corrected by using the reconstruction error between the original dataset and the reconstructed dataset.
the power system, and high load levels indicate a high demand for electricity in the system.Renewable outputs reflect the system's renewable processing, and high renewable outputs indicate sufficient renewable energy resources.

Data-Driven Feature Extraction Using VAE
By maximizing information retention and minimizing reconstruction loss of the decoder, VAE can reduce the hidden layer feature space dimension and extract features.Thus, the negative impact of generation-load uncertainty on subsequent power system dispatch operations can be effectively eliminated.The structural diagram of the feature extraction model based on VAE is shown in Figure 3.As shown in Figure 3, the mean and variance are extracted by the decoder to fit the data distribution.Then, the dataset is reconstructed through a decoder.During the network iteration process, the network is iteratively corrected by using the reconstruction error between the original dataset and the reconstructed dataset.As shown in Figure 3, the VAE contains an encoder and a decoder.The encoder is the inferred network ( | ) q z x ρ with parameter ρ , which is used to generate the variable probability distribution of the hidden variable z.The decoder is the generative network ( | ) p x z θ with parameter θ , which reduces the approximate probability distribution of the original data z.The two probability distribution density models are constructed by multilayer neural networks [25][26][27][28], which can effectively characterize the high-dimensional data of the operation section.
The specific steps of VAE-based data-driven feature extraction proposed in this paper are as follows: (1) Pre-processing the original power system historical operation section data so that the data satisfy a normal distribution with a mean of 0 and a variance of 1. (2) Using the pre-processed data as the input to the VAE model, the encoder maps the input data to the mean vector and variance vector in the latent space.(3) Sampling from the latent space to obtain the hidden variables containing Gaussian noise.(4) The decoder maps the hidden variables in Step 3 back to the original high-dimensional data space and calculates the reconstruction loss.As shown in Figure 3, the VAE contains an encoder and a decoder.The encoder is the inferred network q ρ (z x) with parameter ρ, which is used to generate the variable probability distribution of the hidden variable z.The decoder is the generative network p θ (x|z) with parameter θ, which reduces the approximate probability distribution of the original data z.The two probability distribution density models are constructed by multilayer neural networks [25][26][27][28], which can effectively characterize the high-dimensional data of the operation section.
The specific steps of VAE-based data-driven feature extraction proposed in this paper are as follows: (1) Pre-processing the original power system historical operation section data so that the data satisfy a normal distribution with a mean of 0 and a variance of 1. (2) Using the pre-processed data as the input to the VAE model, the encoder maps the input data to the mean vector and variance vector in the latent space.(3) Sampling from the latent space to obtain the hidden variables containing Gaussian noise.(4) The decoder maps the hidden variables in Step 3 back to the original high-dimensional data space and calculates the reconstruction loss.(5) Minimizing the loss function using the backpropagation algorithm, updating the model parameters by stochastic gradient descent until the reconstruction loss is less than a threshold, and ensuring that the maximization evidence is under session.
The hidden variables obtained by the above steps are used as key data-driven features to characterize the high-dimensional data of the power system operation section data.

Physical Features Considering Generation-Load Features
Generation and load uncertainty is the main factor causing the difference in data distribution of the power system operation sections.The limitations of the neural networkbased VAE data dimensionality reduction method are more prominent when constructing features that divide multiple scenarios.Although its powerful nonlinear transformation capability can effectively extract the data-driven features of high-dimensional data, the data-driven features only consider the complex mapping relationship from the data level, which inevitably causes a certain degree of loss of the original information.Therefore, it is very necessary to consider the key physical features that can effectively portray the generation-load features, fully explore the information related to the generation-load uncertainty in the power system, and construct more refined scenario division features to fully reflect the variability of different operation sections.
During the Chinese Spring Festival in 2023, the maximum penetration rate of renewable energy in Guangdong Province reached 30%, which is a historical peak level.Under the influence of strong cold air, renewable output increases steeply from 2600 MW to 9190 MW, causing a lack of system adjustment capacity.This example demonstrates the complexity and fast changes of operational scenarios under increasing penetration of renewable energy.So, it is necessary to select key physical features to reflect the variability of power system operation sections considering the features of generation and load.According to refs.[29][30][31][32], the key physical features that consider the generation and load features are selected as shown in Table 1.In this paper, the above physical features are combined with VAE-based extraction of operation section features to construct features that divide multiple scenarios.This can improve the overall performance of the data-driven dispatch model based on scenario division.

GMM-Based Power System Scenario Division
The data-driven features extracted by VAE and the key physical features considering the generation-load features are used as the basis for scenario division.The GMM-based scenario division method is proposed in this section.The divided scenarios will support the construction of dispatch models later.
Firstly, based on the Akaike Information Criterion (AIC), the optimal number of scenario categories considering the fitting accuracy and the number of parameters is as follows: Then, the GMM is used to divide the features, and the screened features are input into the GMM for scenario division; each operational scenario category and operation section with obvious differences in generation-load features can be obtained, where the probability distribution of GMM is as follows: The number of sub-Gaussian models in the mixed model is K; P(•) is the probability density function; p(x|µ k , ∑ k ) , µ k , ∑ k , and α k are the first probability distribution model and its expectation, variance (or covariance), and model distribution weights, respectively, satisfying α k ≥ 0 and Finally, the SC [32] is used as a criterion for the effectiveness of scenario division to verify the rationality of scenario division.The contour coefficient reflects the degree of data density within the same scenario and the degree of data dispersion between different scenarios and can be used as an evaluation index to judge the effectiveness of scenario division: where a is the average distance between the current data and other data in the same category; b is the average distance between the current data and the closest data in another category.

5.
Power System Dispatch with Improved Scenario Division 5.1.Reinforcement Learning-Based System Dispatch Model The power system dispatch problem is a sequential decision-making problem which can be modeled as a Markovian decision process and solved using reinforcement learning.The Markov decision process is constructed as follows.

State Space
The state space at the time instant t includes the powers of thermal units, renewable units, loads, and branches: where P Gt and Q Gt are the active and reactive powers of unit G at the time t; P Lt and Q Lt are the active and reactive powers of the load at each node at the time t; V m t is the voltage amplitude of each node at the time t; Loss t is the grid loss (network loss value) at the time t; rho t indicates the actual branch current versus the rated current capacity at the time t.

Action Space
The agent action space A t at time instant t contains the active power P p Gt for thermal units and the active power P new Gt for renewable units:

Reward
To meet the multi-dispatch objectives of safety, reliability, economy, and environmental protection of the power system, the reward function is set as follows: 1.
Line-crossing limit situation (positive reward): where n line denotes the number of branches and rho lt denotes the current loading rate of the branch l at the time t.

2.
Consumption of renewable units (positive incentives): where N N denotes the number of renewable units; P Nkt and P max Nkt denote the active output and maximum output of k renewable units at the time t, respectively.

3.
Balance unit power over the limit (negative reward): where n balanced denotes the number of balancing units; p it denotes the actual active output and the maximum output of balancing units i at the time t.

4.
Unit operating costs (negative incentives): where N G denotes the total number of units, P Git denotes the actual active output of the unit i, and a i , b i , and c i denote the cost factor of the unit i, respectively.5.
Time steps (positive reward): where timestep is the step length of the agent's exploration.
The reward terms r 4t and r 5t are normalized as follows: In summary, the domain values of the reward terms r 1t , r 2t , and r 5t are [0, 1], and the domain values of the reward terms r 3t and r 4t are [-1, 0].The total reward formula is as follows: where r it is the award of the i term at time t, and a i is the coefficient.The setting coefficient a i appropriately can balance the weight of each reward function and avoid overlooking certain rewards.For certain specific goals, a i can also help the agent focus on a certain type of goal and focus on learning, which can enable the agent to learn the best strategy faster.In this paper, considering the dispatch needs of safety, green, and economy, a i is set as follows: where r 3t reflects the safety of dispatch decisions, so a 3 is set to 4; r 2t reflects the green environmental protection degree and economy of dispatch decisions, so a 2 is set to 2. In this paper, other goals are given secondary attention, so they are set to 1.

Reinforcement Learning-Based System Dispatch Model
The operational scenario-specific power system dispatch model based on the PPO algorithm is constructed in this subsection.

Off-Line Training
Each operational scenario agent is established and trained through continuous interaction with the power system simulation environment.The training steps are shown in Figure 4.

Online Decision
The execution of the dispatch takes the following steps: 1.
The dispatching center obtains the power system operation section state information. 2.
The power system operational scenario features are extracted based on the VAE, which fuses key physical features, and the GMM is used to determine the operational scenarios to which they belong.

3.
The scenario agent makes decisions based on the obtained power system operational state information, outputs joint actions, and issues power dispatching adjustment instructions.

4.
Each generating unit enters a new state after output regulation according to the schedule and at the same time, feeds the reward to the scenario agent.

Online Decision
The execution of the dispatch takes the following steps: 1.The dispatching center obtains the power system operation section state information.
2. The power system operational scenario features are extracted based on the VAE, which fuses key physical features, and the GMM is used to determine the operational scenarios to which they belong.3. The scenario agent makes decisions based on the obtained power system operational state information, outputs joint actions, and issues power dispatching adjustment instructions.4. Each generating unit enters a new state after output regulation according to the schedule and at the same time, feeds the reward to the scenario agent.

Simulation Setup
To verify the effectiveness of the proposed scenario division-based power system dis-

Simulation 6.1. Simulation Setup
To verify the effectiveness of the proposed scenario division-based power system dispatching method, a real power system model in a province of China is used for validation.The system contains 748 nodes, 55 renewable units, and 187 thermal units with a dispatch interval of 15 min.300 days of data are selected from the grid historical data, and the dataset is pre-processed and normalized before training.
The simulations are conducted in Pytorch framework on a PC with an i7-11700 CPU and 32 GB memory.The encoder and decoder of VAE use the same network structure with Relu activation function, a learning rate of 0.00001, and 9000 training rounds.

Simulation Result Analysis
(1) Feature extraction result based on VAE.The system operation section has 303 dimensions of information (containing 55 renewable units and 248 loads).The trend of the loss function of the VAE model is shown in Figure 5.As shown in Figure 5, at the initial stage of 0-1000 iterations, the loss function gradually increases.This is because in the initial stage, the model is learning the feature representation of the data and attempting to minimize reconstruction errors.After approximately 2000 iterations, the loss function begins to converge.This means that the VAE model has learned about the potential structure and feature expression of the data.In the simulation of the paper, 48 dimensional data-driven features were extracted from 303 dimensional original data through the VAE model.The simulations are conducted in Pytorch framework on a PC with an i7-11700 CPU and 32 GB memory.The encoder and decoder of VAE use the same network structure with Relu activation function, a learning rate of 0.00001, and 9000 training rounds.

Simulation Result Analysis
(1) Feature extraction result based on VAE.The system operation section has 303 dimensions of information (containing 55 renewable units and 248 loads).The trend of the loss function of the VAE model is shown in Figure 5.As shown in Figure 5, at the initial stage of 0-1000 iterations, the loss function gradually increases.This is because in the initial stage, the model is learning the feature representation of the data and attempting to minimize reconstruction errors.After approximately 2000 iterations, the loss function begins to converge.This means that the VAE model has learned about the potential structure and feature expression of the data.In the simulation of the paper, 48 dimensional data-driven features were extracted from 303 dimensional original data through the VAE model.(2) Scenario division results based on GMM.The optimal number of scenario division categories is determined to be three, based on the AIC criterion.The 300−day data are divided into three scenarios based on the GMM; the results of scenario division are shown in Figure 6.The basic idea of the GMM is to decompose data into a mixture of multiple Gaussian distribution components, where each Gaussian distribution component corresponds to a scenario.Therefore, the three types of scenarios divided by features based on the scenario division model represent the probability distribution of the three types of features.As shown in Figure 6, scenario B has relatively fewer days of scenario division compared with scenarios A and C.But the difference in score intervals among different scenarios is obvious, which effectively verifies the rationality of scenario division.(2) Scenario division results based on GMM.The optimal number of scenario division categories is determined to be three, based on the AIC criterion.The 300−day data are divided into three scenarios based on the GMM; the results of scenario division are shown in Figure 6.The basic idea of the GMM is to decompose data into a mixture of multiple Gaussian distribution components, where each Gaussian distribution component corresponds to a scenario.Therefore, the three types of scenarios divided by features based on the scenario division model represent the probability distribution of the three types of features.As shown in Figure 6, scenario B has relatively fewer days of scenario division compared with scenarios A and C.But the difference in score intervals among different scenarios is obvious, which effectively verifies the rationality of scenario division.Table 2 shows the comparison of the proposed method and the scenario division only based on VAE extraction features by using the SC values.From the table, it can be seen that the features constructed based on the method of this paper can effectively improve the division effect, and the correlation of each operation section in the operational scenario is higher.Table 2 shows the comparison of the proposed method and the scenario division only based on VAE extraction features by using the SC values.From the table, it can be seen that the features constructed based on the method of this paper can effectively improve the division effect, and the correlation of each operation section in the operational scenario is higher.(3) Dispatch results based on the proposed method.

VAE loss
To verify the effectiveness of the proposed method, the proposed method was compared with the convex optimization problem (OPT) and single scenario methods (original scenario), respectively.The comparison of results from different scenarios and methods are shown in Figure 7 and Table 3.As shown in Figure 7, from the comparison of the mean reward of scenario A and C based on the proposed method with the original scenario, it can be seen that scenario division can effectively improve the training effectiveness of agents and accelerate agents learning efficiency.The reason why the mean reward of scenario B is lower than the mean reward of the original scenario is that there are fewer samples in scenario B, making it difficult for reinforcement learning agents to achieve optimal decision making through training.In addition, the reason why the mean reward based on the OPT is higher than the mean reward based on the proposed method is that the OPT is usually based on mathematical models and objective functions for solving and have the ability to achieve global optimization.In contrast, reinforcement learning is a method of searching through trial and error and learning, which may make it difficult to find the global optimal solution in problems with large search space and complexity.But the average single-step solving time of the OPT is 19.82 s, and the decision time of the method proposed in this paper is 0.003 s.Compared to the OPT, the decision-making efficiency of the method proposed has been improved by 99.98%.In summary, the method proposed in this paper can effectively improve the efficiency of the reinforcement learning agent's training and decision making.To further validate the advantages of the proposed method, a comparison was made between the proposed method and the method based on VAE.The comparison results are shown in Figure 8.As shown in Figure 8, the mean reward values of the method proposed in this paper are significantly higher than the method based on VAE during the agent training process in both scenario A and scenario C. Compared to the method based on VAE, the mean reward value of the proposed method in scenario A is increased by 43.51% and that in scenario C is increased by 35.3%.Overall, the method proposed in this paper can prevent agents from making decisions directly on the dispatch scenarios with uncertainty distribution and reduce the learning difficulty.

Method
The Proposed Method Traditional Optimization Methods Average single-step decision time 0.003 s 19.82 s To further validate the advantages of the proposed method, a comparison was made between the proposed method and the method based on VAE.The comparison results are shown in Figure 8.As shown in Figure 8, the mean reward values of the method proposed in this paper are significantly higher than the method based on VAE during the agent training process in both scenario A and scenario C. Compared to the method based on VAE, the mean reward value of the proposed method in scenario A is increased by 43.51% and that in scenario C is increased by 35.3%.Overall, the method proposed in this paper can prevent agents from making decisions directly on the dispatch scenarios with uncertainty distribution and reduce the learning difficulty.Based on the analysis of power system dispatch results, the output power dispatch results of renewable energy with and without scenario division are shown Figure 9.The renewable energy consumption rate refers to the ratio of the actual output power of renewable energy used by the system to the total output power of renewable energy at the current time.As shown in Figure 9, the renewable energy consumption rate of nodes 13-22 with scenario division are significantly improved compared to without scenario division.Overall, the decision results of the proposed method are better in renewable energy consumption.
From the perspective of operating costs in the power system dispatch, unit operating cost with and without the scenario division are shown in Figure 10.As shown in Figure 10, the unit operating cost without the scenario division is USD 1,335,901, and the unit operating cost after scenario division is USD 1,266,477.The overall power system operating cost is reduced by 5.196%, and the peak of the unit operating cost is reduced.
renewable energy consumption rate refers to the ratio of the actual output power of renewable energy used by the system to the total output power of renewable energy at the current time.As shown in Figure 9, the renewable energy consumption rate of nodes 13-22 with scenario division are significantly improved compared to without scenario division.Overall, the decision results of the proposed method are better in renewable energy consumption.From the perspective of operating costs in the power system dispatch, unit operating cost with and without the scenario division are shown in Figure 10.As shown in Figure 10, the unit operating cost without the scenario division is USD 1,335,901, and the unit operating cost after scenario division is USD 1,266,477.The overall power system operating cost is reduced by 5.196%, and the peak of the unit operating cost is reduced.From the perspective of operating costs in the power system dispatch, unit operating cost with and without the scenario division are shown in Figure 10.As shown in Figure 10, the unit operating cost without the scenario division is USD 1,335,901, and the unit operating cost after scenario division is USD 1,266,477.The overall power system operating cost is reduced by 5.196%, and the peak of the unit operating cost is reduced.

Conclusions
A multi-scenario data-driven dispatch method based on scenario division integrating VAE-generated and physical features is proposed and verified by simulations.
To address the power system dispatch problem under multiple and fast-changing operational scenarios, firstly, VAE is used to extract the potential data-driven features, which can reduce the redundancy of high-dimensional features and accurately describe the operation features of generation-load data to the greatest extent.At the same time, the key physical features are selected, considering the features of generation-load, which are interpretable and can screen the operation features from the overall level and changing trend of generation-load.This provides a more comprehensive and interpretable feature.Based on these features, the GMM is used to divide the scenarios and distinguish the scenarios with obvious differences in generation-load operation features.Using SC to evaluate the effectiveness of scenario division.Finally, the dispatching model of each scenario is constructed and the precise dispatching decision under each scenario is given.Compared with a single scenario dispatching decision, this method is more accurate and reasonable for the actual dispatching requirements.

Figure 1 .
Figure 1.The framework of power system dispatch with scenario division.

Figure 1 .
Figure 1.The framework of power system dispatch with scenario division.

Figure 1 .
Figure 1.The framework of power system dispatch with scenario division.

Figure 3 .
Figure 3. Structural diagram of feature extraction model based on VAE.

Figure 3 .
Figure 3. Structural diagram of feature extraction model based on VAE.

Figure 6 .
Figure 6.Clustering results using the Gaussian mixture model.

Energies 2023 , 16 Figure 7 .Table 3 .Figure 7 .
Figure 7.Comparison of mean reward change curves for different scenarios and methods.Table3.Comparison of dispatch time between the traditional optimization method and the proposed method.MethodThe Proposed Method Traditional Optimization Methods Average single-step 0.003 s 19.82 s

Figure 7 .
Figure 7.Comparison of mean reward change curves for different scenarios and methods.

Figure 8 .
Figure 8.Comparison of rewards under different methods.

Figure 9 .
Figure 9. Output power results of renewable energy with and without scenario division.

Figure 10 .
Figure 10.Comparison of unit operating cost with and without the scenario division.

Figure 9 .
Figure 9. Output power results of renewable energy with and without scenario division.

Figure 9 .
Figure 9. Output power results of renewable energy with and without scenario division.

Figure 10 .
Figure 10.Comparison of unit operating cost with and without the scenario division.Figure 10.Comparison of unit operating cost with and without the scenario division.

Figure 10 .
Figure 10.Comparison of unit operating cost with and without the scenario division.Figure 10.Comparison of unit operating cost with and without the scenario division.

Table 1 .
Key physical features and implications.

Table 2 .
Comparison of clustering effect based on SC values.

Table 3 .
Comparison of dispatch time between the traditional optimization method and the proposed method.

Table 3 .
Comparison of dispatch time between the traditional optimization method and the proposed method.