Virtual power plant low-carbon economic dispatching method and system

Abstract

The present application relates to the technical field of power system dispatching, in particular to a virtual power plant low-carbon economic dispatching method and system, the method first acquires historical output data of each virtual power plant, obtains output prediction data considering uncertainty according to time series prediction; then constructs a carbon emission reduction responsibility allocation space, determines the carbon emission reduction responsibility weight of each virtual power plant through topological homeomorphism mapping, solves the target function with the lowest carbon emission cost combining the carbon emission constraint index, and obtains the optimal output planning sequence; finally, a discrete-continuous hybrid dynamics model is constructed to realize the collaborative scheduling of multiple energy storage devices and handle the non-integer multiple time difference value of different scheduling periods, the present application constructs a dynamic carbon emission reduction responsibility allocation mechanism through topological mapping theory, quantifies the output prediction uncertainty using probability measure theory, and optimizes the scheduling of each virtual power plant.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of power system dispatching, in particular to a virtual power plant low-carbon economic dispatching method and system for realizing the economic optimal operation of a virtual power plant group under the condition of considering carbon emission constraints. BACKGROUND

[0002] As a new type of power organization form, a virtual power plant can effectively improve the consumption capacity of renewable energy and reduce the overall carbon emission of a system by integrating various distributed energy resources and controllable loads.

[0003] The existing virtual power plant dispatching method mainly has the following problems: firstly, the carbon emission reduction responsibility allocation mechanism is rigid, usually adopting a simple linear allocation method, which cannot adapt to the actual situation of different virtual power plants and the dynamically changing operation environment; secondly, the output prediction uncertainty is not handled well, mostly adopting a deterministic prediction or a simple interval estimation, which cannot accurately quantify and propagate the uncertainty; thirdly, the multi-energy storage device collaborative dispatching is not optimized, especially there is a fault problem in the connection processing of different dispatching periods; and fourthly, the large-scale virtual power plant optimization solving efficiency is low, which cannot meet the real-time dispatching demand.

[0004] Therefore, there is an urgent need for a virtual power plant low-carbon economic dispatching method and system capable of realizing dynamic fair allocation of carbon emission reduction responsibility, accurately handling output prediction uncertainty, and optimizing multi-energy storage device collaborative dispatching. SUMMARY

[0005] The purpose of the present application is to provide a virtual power plant low-carbon economic dispatching method and system, aiming to solve the technical problems of unreasonable carbon emission reduction responsibility allocation, insufficient handling of output prediction uncertainty, and insufficient optimization of multi-energy storage device collaborative dispatching in the prior art.

[0006] The present application provides a virtual power plant low-carbon economic dispatching method, which comprises the following steps:

[0007] Obtaining the historical output data of each virtual power plant in the power grid within a set dispatching period, predicting the output data of each virtual power plant within the set dispatching period according to time series, and obtaining the output prediction data considering uncertainty;

[0008] Based on the output prediction data considering uncertainty, a carbon emission reduction responsibility allocation space is constructed, each virtual power plant is mapped as a vector point in the space, the carbon emission reduction responsibility weight of each virtual power plant is determined through topological homeomorphism mapping, the determined carbon emission reduction responsibility weight is substituted into the carbon emission constraint index, an optimization problem with the target of minimizing the carbon emission cost is constructed, and the output optimal planning sequence of each virtual power plant within the set dispatching period is obtained by solving the optimization problem;

[0009] According to the output prediction data considering uncertainty and the output optimal planning sequence of each virtual power plant, a discrete-continuous hybrid dynamics model is constructed to realize the collaborative scheduling of multiple energy storage devices, handle the time difference non-integer multiple situation of different scheduling periods, optimize the scheduling of each virtual power plant, and optimize the power grid.

[0010] As preferred, the step of constructing a carbon emission reduction responsibility allocation space based on the output prediction data considering uncertainty specifically includes:

[0011] An n-dimensional responsibility vector space is constructed, where n is the number of virtual power plants, and each virtual power plant is represented as a vector point in the space;

[0012] Total responsibility constraints, non-negative constraints and upper limit constraints are introduced to form a convex polyhedron in the responsibility space, representing all feasible responsibility allocation schemes;

[0013] A time parameter is introduced into the responsibility space to dynamically adjust the responsibility allocation over time, and the smooth transition of the responsibility allocation is maintained through continuous deformation.

[0014] As preferred, the step of determining the carbon emission reduction responsibility weight of each virtual power plant through topological homeomorphism mapping specifically includes:

[0015] A carbon quota space is defined, where points represent quota allocation schemes;

[0016] A mapping function from the carbon quota space to the responsibility space is constructed to ensure that the mapping satisfies the continuity condition and reversibility;

[0017] The mapping function is designed to follow the weight guiding principle, quota adaptation principle, fairness principle and efficiency principle;

[0018] Real-time monitoring of carbon dioxide intensity and power grid operating state variables, and distributing environmental change information to each virtual power plant through a parameter propagation network;

[0019] Adjustment step and adjustment direction vector are defined to realize the iterative adjustment of responsibility allocation, and constraint projection is used to ensure that the adjusted responsibility allocation is still within the feasible region.

[0020] As preferred, the step of obtaining the output prediction data considering uncertainty specifically includes:

[0021] The output of each virtual power plant at time t is modeled as a random variable, and a stochastic process is constructed to describe the time evolution of the output;

[0022] An output state space and a probability measure are defined to form a complete probability space;

[0023] The statistical characteristics of historical output time series are analyzed to identify periodicity, trend and randomness;

[0024] Constructing the output distribution model, including the marginal distribution and conditional distribution;

[0025] Calculating the point prediction value based on the conditional expectation, and quantifying the uncertainty of the prediction to construct the prediction interval.

[0026] As a preferred, the carbon emission constraint index specifically includes:

[0027] The carbon emission and responsibility matching constraint ensures that the difference between the total carbon emission of each virtual power plant and the carbon emission reduction responsibility value does not exceed the set threshold;

[0028] The total emission constraint ensures that the total emission of the system does not exceed the upper limit;

[0029] The equipment operation constraint includes the upper and lower limits of the output, the climbing rate constraint;

[0030] The energy storage device specific constraint includes the charging and discharging rate, the state of charge range, etc.;

[0031] By introducing a disturbance set to represent all possible prediction errors, a robust constraint is constructed to ensure that the constraint still satisfies under uncertainty conditions.

[0032] As a preferred, the step of solving the objective function with the lowest carbon emission cost specifically includes:

[0033] Constructing the objective function f, considering factors such as the discounted price of the energy storage battery, the dispatching period, the energy storage efficiency, the energy consumption efficiency, the discount rate and the carbon price, etc.

[0034] The overall optimization problem is decomposed into a network of sub-problems with self-similarity, forming a fractal structure;

[0035] Designing a hierarchical solving strategy, including solving the basic sub-problems, intermediate sub-problems and global problems from bottom to top;

[0036] Implementing an interlayer information transmission mechanism, and designing a convergence judgment standard to control the iteration process;

[0037] The global coordinator supervises the overall optimization process to solve the conflict between local optimal and global optimal.

[0038] As a preferred, the step of constructing a discrete-continuous hybrid dynamics model according to the output prediction data considering uncertainty and the output optimal planning sequence of each virtual power plant specifically includes: inputting the output prediction data considering uncertainty and the output optimal planning sequence of each virtual power plant into the discrete-continuous hybrid dynamics model to obtain the output optimal planning sequence of each virtual power plant. As the initial condition of the system state, the output optimal planning sequence is inputted into the discrete-continuous hybrid dynamics model to obtain the output optimal planning sequence of each virtual power plant. As a system control target; design the state transition of the discrete part to describe the scheduling decision point; design the continuous part to describe the dynamic evolution of the system between decision points; construct a state variable vector to represent the system state; calculate the Lyapunov index spectrum of the system to analyze the structural characteristics of the chaotic attractor; define the fast time scale and the slow time scale, optimize the overall scheduling strategy on the slow time scale, and perform real-time adjustment and control on the fast time scale.

[0039] As preferred, the step of realizing the collaborative scheduling of multiple energy storage devices specifically includes:

[0040] Defining an energy storage state vector, including state of charge, charging and discharging power, and working mode;

[0041] Defining a state transition function to describe the change of state under action;

[0042] Construct a set of state constraints to ensure that the state is within a safe range;

[0043] Design a hierarchical control structure, including formulating an overall scheduling strategy, allocating scheduling tasks for each energy storage device, and executing specific charging and discharging control;

[0044] Implement a load sharing mechanism and fault tolerance strategy between devices to ensure that a single device failure does not affect the overall scheduling.

[0045] As preferred, the step of processing the non-integer multiple of the time difference between different scheduling periods specifically includes:

[0046] Defining the time deviation as the difference between the integer multiple of the start time of the second scheduling period and the start time of the first scheduling period;

[0047] Design a time interpolation function to generate a transition state at a non-integer multiple of time;

[0048] Implement a state smoothing transition strategy to avoid scheduling faults;

[0049] Construct a state reconstruction method to maintain the continuity of scheduling;

[0050] When the difference between the start time of the second scheduling period and the start time of the first scheduling period is not an integer multiple, calculate the charging and discharging scheduling operation of the second energy storage battery, and execute the operation.

[0051] A virtual power plant low-carbon economic scheduling system, comprising:

[0052] An output prediction module for obtaining historical output data of each virtual power plant in the power grid within a set scheduling period, predicting the output data of each virtual power plant within the set scheduling period according to time series, and obtaining output prediction data considering uncertainty;

[0053] An optimal planning solving module is configured to construct a carbon emission reduction responsibility distribution space based on the output prediction data considering uncertainty, map each virtual power plant as a vector point in the space, determine the carbon emission reduction responsibility weight of each virtual power plant through topological homeomorphism mapping, solve a target function with the lowest carbon emission cost in combination with a carbon emission constraint index, and obtain an optimal output planning sequence of each virtual power plant in a set scheduling period.

[0054] An optimal scheduling module is configured to construct a discrete-continuous hybrid dynamic model based on the output prediction data considering uncertainty and the optimal output planning sequence of each virtual power plant, realize collaborative scheduling of multiple energy storage devices, and handle a time difference non-integer multiple situation of different scheduling periods, to perform optimal scheduling on each virtual power plant and optimize the power grid.

[0055] The present application has the following beneficial effects:

[0056] 1. A carbon emission reduction responsibility dynamic distribution mechanism based on topological mapping is constructed, which can realize fair distribution and smooth adjustment of responsibility according to the actual situation of the virtual power plant and the dynamically changing operating environment, and improve the rationality and adaptability of the distribution compared with the traditional static distribution method.

[0057] 2. An uncertainty output prediction and constraint construction method based on probability measure is introduced, which realizes accurate quantification and propagation of uncertainty, improves prediction accuracy and system robustness, and improves the ability of the system to adapt to large-scale disturbances by 30% while reducing prediction deviation by 15%-20%.

[0058] 3. A discrete-continuous hybrid dynamic model and a multi-dimensional collaborative scheduling mechanism are adopted, which solves the fault problem in the connection processing of different scheduling periods, improves the utilization efficiency of the energy storage device by about 25%, prolongs the service life of the device and reduces the maintenance cost.

[0059] 4. A multi-level optimization solving method with fractal structure is applied, which reduces the calculation complexity from O(n³) to O(n·logn), supports simultaneous scheduling of hundreds of virtual power plants, and reduces the operating cost by 7%-12% under the condition of meeting the same carbon emission reduction target. BRIEF DESCRIPTION OF DRAWINGS

[0060] Figure 1 A flowchart of the virtual power plant low-carbon economic scheduling method of the present application;

[0061] Figure 2 A schematic diagram of the uncertainty output prediction and constraint construction method based on probability measure of the present application;

[0062] Figure 3 A schematic diagram of the multi-level optimization solving method with fractal structure of the present application;

[0063] Figure 4The structure block diagram of the virtual power plant low-carbon economic dispatching system is shown in the figure. DETAILED DESCRIPTION

[0064] Please refer to Figure 1 - Figure 4 The application will be further described below in conjunction with the drawings and examples.

[0065] Referring to Figure 1 The application provides a virtual power plant low-carbon economic dispatching method, comprising the following steps:

[0066] Firstly, the historical output data of each virtual power plant in the power grid within a set dispatching period is obtained, the output data of each virtual power plant within the set dispatching period is predicted according to the time sequence, and the output prediction data considering uncertainty is obtained.

[0067] Then, based on the output prediction data considering uncertainty, a carbon emission reduction responsibility allocation space is constructed, each virtual power plant is mapped as a vector point in the space, the carbon emission reduction responsibility weight of each virtual power plant is determined through topological homeomorphism mapping, the objective function with the lowest carbon emission cost is solved in combination with the carbon emission constraint index, and the optimal planning sequence of the output of each virtual power plant within the set dispatching period is obtained.

[0068] Finally, according to the output prediction data considering uncertainty and the optimal planning sequence of the output of each virtual power plant, a discrete-continuous hybrid dynamics model is constructed, the collaborative dispatching of multiple energy storage devices is realized, the time difference non-integer multiple situation of different dispatching periods is handled, each virtual power plant is optimally dispatched, and the power grid is optimized.

[0069] In an embodiment of the application, the step of constructing the carbon emission reduction responsibility allocation space based on the output prediction data considering uncertainty specifically comprises: constructing an n-dimensional responsibility vector space, wherein n is the number of virtual power plants, each virtual power plant is represented as a vector point in the space; introducing total responsibility constraints, non-negative constraints and upper limit constraints to form a convex polyhedron in the responsibility space, representing all feasible responsibility allocation schemes; introducing a time parameter into the responsibility space to dynamically adjust the responsibility allocation with time, and maintaining the smooth transition of the responsibility allocation through continuous deformation.

[0070] Specifically, the n-dimensional responsibility vector space can be expressed as , wherein the basis vector corresponds to the n virtual power plants. Any point represents a possible responsibility allocation scheme, wherein represents the carbon emission reduction responsibility value of the i-th virtual power plant. In the space, a metric function is defined to measure the difference between different responsibility allocation schemes, for example, the Euclidean distance can be used:

[0071] ,

[0072] wherein, is the distance between two responsibility allocation schemes and , is the carbon reduction responsibility value of the i-th virtual power plant in scheme , is the carbon reduction responsibility value of the i-th virtual power plant in scheme , n is the total number of virtual power plants, denotes the summation operation over all virtual power plants.

[0073] The constrained boundary of the responsibility space is jointly defined by the following three constraints:

[0074] 1. Total responsibility constraint: wherein denotes the summation of carbon reduction responsibility values of all virtual power plants, is the total amount of carbon reduction responsibility, this constraint ensures that the total sum of responsibilities of all virtual power plants is equal to the total amount of responsibility set by the system;

[0075] 2. Non-negative constraint: wherein is the carbon reduction responsibility value of the i-th virtual power plant, this constraint ensures that the responsibility value of each virtual power plant is non-negative;

[0076] 3. Upper limit constraint: wherein is the carbon reduction responsibility value of the i-th virtual power plant, is the maximum bearable responsibility of the i-th virtual power plant, this constraint prevents individual virtual power plants from bearing too much responsibility.

[0077] These constraints jointly constitute a convex polyhedron D in the responsibility space, representing all feasible responsibility allocation schemes. To achieve dynamic adjustment of responsibility allocation, a time parameter t is introduced, making the point P(t) in the responsibility space change over time, forming a responsibility trajectory γ(t). In practical applications, the adjustment frequency of responsibility allocation can be set according to system requirements, such as once an hour, once a day, or once a week. Preferably, the adjustment frequency can be set to once every 4 hours, which can timely respond to changes in system state, and also avoid system instability caused by too frequent adjustments.

[0078] In another embodiment of the present application, the step of determining the carbon emission reduction responsibility weight of each virtual power plant by topological homeomorphism mapping specifically comprises: defining a carbon quota space, wherein a point represents a quota allocation scheme; constructing a mapping function from the carbon quota space to the responsibility space, so that the mapping satisfies the continuity condition and the reversibility; designing the mapping function to follow the weight guiding principle, the quota adaptation principle, the fairness principle and the efficiency principle; monitoring the carbon dioxide intensity and the power grid operation state variables in real time, distributing environmental change information to each virtual power plant through a parameter propagation network; defining an adjustment step and an adjustment direction vector, realizing iterative adjustment of responsibility allocation, and ensuring that the adjusted responsibility allocation is still within the feasible region through constraint projection.

[0079] Specifically, a carbon quota space Q is defined, wherein a point represents a quota allocation scheme, is the carbon quota of the i-th virtual power plant. A mapping function is constructed, which maps the quota allocation to the responsibility allocation. The design of the mapping function needs to satisfy the following conditions:

[0080] 1. Continuity condition: a small change in the quota leads to a small change in the responsibility, that is, for any , there exists such that when , there is ; wherein represents the distance between two points and in the quota space, represents the distance between the corresponding points and in the responsibility space, and are positive real numbers;

[0081] 2. Reversibility: different quota allocations correspond to different responsibility allocations, that is, for any , there is ; wherein and are different points in the quota space, and are the corresponding points in the responsibility space.

[0082] The specific form of the mapping function can be designed as:

[0083] ,

[0084] wherein is the carbon emission reduction responsibility value of the i-th virtual power plant, is its carbon quota, is its carbon emission reduction responsibility weight, is the total amount of carbon emission reduction responsibility, and is the weight coefficient, satisfying , , is the total carbon quota of all virtual power plants, is the proportion of the carbon quota of the i-th virtual power plant in the total quota.

[0085] In practical applications, and can be flexibly adjusted according to system requirements. For example, when it is desired to give more consideration to the weight factor, the value of can be increased; when it is desired to give more consideration to the quota factor, the value of can be increased. Preferably, , so that both the weight factor and the quota factor are considered, achieving a relatively balanced responsibility allocation.

[0086] In order to realize dynamic adjustment of responsibility allocation, it is necessary to monitor the carbon dioxide intensity and the grid operation state variables (including load level, renewable energy output, etc.) in real time, and distribute environmental change information to each virtual power plant through a parameter propagation network. The iterative adjustment process of responsibility allocation can be represented as:

[0087] ,

[0088] wherein, is the responsibility allocation vector at time t, containing the responsibility values of all virtual power plants, is the responsibility allocation vector at time , is the adjustment step, which is a non-negative scalar, is the adjustment direction vector, indicating the direction of responsibility adjustment, is the time step.

[0089] The adjustment step can be dynamically set according to the system state, for example, when the system state changes dramatically, the step can be appropriately reduced to ensure the stability of the adjustment; when the system state changes gently, the step can be appropriately increased to speed up the convergence. Preferably, the adjustment step can be set between 0.01 and 0.1, and the specific value can be dynamically adjusted according to actual conditions.

[0090] In order to ensure that the adjusted responsibility allocation is still within the feasible domain D, a constraint projection is needed:

[0091] ,

[0092] wherein, is the projected responsibility allocation vector, This represents the projection operation onto a convex polyhedron D, ensuring that the projected points lie within the feasible region and satisfy all constraints. Constrained projection ensures that responsibility assignment always satisfies the overall responsibility constraint, non-negativity constraint, and upper bound constraint.

[0093] In another embodiment of the present invention, the steps of obtaining power output prediction data considering uncertainty specifically include: modeling the power output of each virtual power plant at time t as a random variable, constructing a stochastic process to describe the temporal evolution of the power output; defining the power output state space and probability measure to form a complete probability space; analyzing the statistical characteristics of historical power output time series to identify periodic, trend and random components; constructing a power output distribution model, including marginal distribution and conditional distribution; calculating point prediction values ​​based on conditional expectation, quantifying the uncertainty of the prediction, and constructing a prediction interval.

[0094] Specifically, the output of the i-th virtual power plant at time t is modeled as a random variable. Constructing a stochastic process Describe the time evolution of the output force. Define the output state space S and the probability measure P, and introduce... - Algebra F forms a complete probability space Where S is the set of all possible output states, and F is the power output state on S. - Algebra (containing all measurable events), where P is a probability measure function that maps events in F to probability values ​​in the interval [0,1].

[0095] Statistical analysis of historical output time series includes calculating the mean, variance, and autocorrelation function. Taking the autocorrelation function as an example, its calculation formula is:

[0096] ,

[0097] in, The autocorrelation function represents the interval. The correlation between the two output values ​​over time This represents the expectation operation. Let i be the output of the i-th virtual power plant at time t. For its in Constant effort The average value of the output. The variance of output, This is due to the time delay. By analyzing the autocorrelation function, the periodicity of the force can be identified. For example, for photovoltaic power generation, the autocorrelation function typically peaks at integer multiples of 24 hours, reflecting the diurnal periodicity of photovoltaic power generation.

[0098] The power output distribution model can be constructed using either parametric or non-parametric methods. For parametric methods, a suitable distribution function, such as the normal distribution or Weibull distribution, can be selected based on the statistical characteristics of the power output. For non-parametric methods, techniques such as kernel density estimation can be used to directly estimate the distribution from the data. In practical applications, different types of virtual power plants may have different power output distribution characteristics, requiring the selection of an appropriate distribution model based on specific circumstances.

[0099] The point prediction value calculated based on conditional expectation can be expressed as:

[0100] ,

[0101] in, For the i-th virtual power plant in Predicting output at any given moment Indicates that in the known Under the condition of constant exertion, The expected value of output at any given time is given by k, which is the size of the observation window for historical data. In practical applications, the choice of k needs to balance prediction accuracy and computational complexity, and is generally taken as a value between 12 and 24, that is, considering data from 12 to 24 historical time points.

[0102] The uncertainty of a forecast can be quantified by the forecast interval:

[0103] ,

[0104] in, This is the lower limit of the prediction interval. The upper limit of the prediction interval can be determined by setting the confidence level β. Generally, β = 0.95 is taken, which means that there is a 95% probability that the actual output will fall within the prediction interval.

[0105] In another embodiment of the present invention, the carbon emission constraint indicators specifically include: carbon emission and responsibility matching constraints to ensure that the difference between the total carbon emissions of each virtual power plant and the carbon emission reduction responsibility value does not exceed a set threshold; total emission constraints to ensure that the total system emissions do not exceed the upper limit; equipment operation constraints, including upper and lower limits of output and ramp rate constraints; energy storage equipment specific constraints, including charge and discharge rates and state of charge range; and robust constraints are constructed by introducing a perturbation set to represent all possible prediction errors to ensure that the constraints are still satisfied under uncertainty conditions.

[0106] Specifically, the carbon emission-responsibility matching constraint can be expressed as:

[0107] ,

[0108] in, This represents the total carbon emissions of the i-th virtual power plant. Its carbon emission reduction responsibility value The absolute value of the difference To set a threshold, representing the maximum allowable deviation, This means that this constraint should be satisfied for all virtual power plants. In practical applications, The value can be set according to the system tolerance, and is generally 100%. That is, a deviation of 5% is allowed.

[0109] Total emissions constraints can be expressed as:

[0110] ,

[0111] in, This represents the sum of carbon emissions from all virtual power plants. This represents the maximum allowable total carbon emissions allowed by the system.

[0112] Equipment operating constraints include upper and lower limits for output and gradient rate constraints:

[0113] ,

[0114] ,

[0115] in, Let i be the output of the i-th virtual power plant at time t. For its minimum output limit, Limiting its maximum output, This represents the absolute value of the change in output between two adjacent time points. Let i be the maximum ramp rate of the i-th virtual power plant. This means that these constraints should be met for all virtual power plants and at all points in time.

[0116] Energy storage device-specific constraints include charge / discharge rate constraints and state of charge range constraints:

[0117] ,

[0118] ,

[0119] in, Let t represent the charging and discharging power of the i-th energy storage device at time t (positive value indicates discharging, negative value indicates charging). It is the negative of its maximum charging power. Its maximum discharge power, Let t be its state of charge at time t, representing the percentage of battery charge. This is the lower limit of its state of charge. This is the upper limit of its state of charge. It is indicated that these constraints should be satisfied for all energy storage devices and all time points.

[0120] To handle the uncertainty of the output prediction, a set of disturbances U is introduced to represent all possible prediction errors:

[0121]

[0122] where U is the set of disturbances, which is an ellipsoidal set, u is an n-dimensional prediction error vector, each component corresponds to the prediction error of a virtual power plant, denotes the transpose of u, is the covariance matrix of n x n-dimensional prediction error, which describes the correlation between the prediction errors of different virtual power plants, is the inverse matrix of is the inverse matrix of is a robustness level parameter that controls the size of the ellipsoidal set, is its square value. By adjusting the value of , the degree of conservatism of the robust constraint can be controlled. Generally speaking, the value of can be in the range of 1 to 3, with smaller values corresponding to lower conservatism and larger values corresponding to higher conservatism. Preferably, , the robust constraint can cover about 95% of the prediction error cases, achieving a good balance between conservatism and economy.

[0123] Based on the set of disturbances, a robust constraint is constructed:

[0124]

[0125] where is the constraint function, x is the decision variable vector containing all variables to be optimized, and u is the prediction error vector, indicates that the constraint should be satisfied for any prediction error u in the set of disturbances U. This robust constraint ensures that the constraint condition is satisfied under any possible prediction error, thereby enhancing the robustness of the system.

[0126] In still another embodiment of the present application, the step of solving the objective function with the lowest carbon emission cost specifically comprises: constructing an objective function f, considering factors such as the discounted price of energy storage batteries, scheduling time period, energy storage efficiency, energy consumption efficiency, discount rate and carbon price; decomposing the overall optimization problem into a network of sub-problems with self-similarity to form a fractal structure; designing a hierarchical solving strategy, including solving the basic sub-problems, intermediate sub-problems and global problems from bottom to top; implementing an inter-layer information transmission mechanism and designing a convergence judgment standard to control the iteration process; supervising the overall optimization process through a global coordinator to solve the conflict between local optimum and global optimum.

[0127] ​Specifically, the carbon emission reduction responsibility weight of each virtual power plant is determined by a topological homeomorphism mapping will be substituted into the carbon emission and responsibility matching constraint In combination with the total emission constraint , the equipment operation constraint and the energy storage device specific constraint, an optimization problem is constructed to minimize the carbon emission cost.

[0128] Specifically, the objective function f can be expressed as:

[0129] ,

[0130] where f is the carbon emission cost objective function, represents the summation of all n virtual power plants, represents the summation of all T scheduling periods, is the discounted electricity price obtained by the energy storage battery in the i-th virtual power plant (yuan / kWh), is the charging and discharging power of the energy storage battery in the i-th virtual power plant at period t (kW), is the length of a scheduling period (hours), is the energy consumption efficiency of the energy storage battery (dimensionless, taking a value between 0 and 1), is the energy storage efficiency of the energy storage battery in the i-th virtual power plant at period t (dimensionless, taking a value between 0 and 1), is the total energy consumption of the energy storage battery in the i-th virtual power plant at period t (kWh), and r is the discount rate (dimensionless, representing the time value of money), is the discount factor at period t, and p is the carbon price (yuan / ton), is the total carbon emission of the i-th virtual power plant (tons).

[0131] In practical applications, the values of each parameter need to be determined according to specific circumstances. For example, the discounted electricity price can be set according to the electricity market price and the energy storage characteristics, generally between 0.3 yuan / kWh and 0.8 yuan / kWh; the scheduling period length can be set according to system demand, generally 15 minutes, 30 minutes or 1 hour; the energy storage efficiency is influenced by factors such as battery type and operating state, generally between 0.85 and 0.95; the discount rate can be set according to the time value of money, generally between 0.05 and 0.1; the carbon price is influenced by market factors, generally between 50 yuan / ton and 200 yuan / ton.

[0132] The overall optimization problem is decomposed into a network of self-similar sub-problems, forming a fractal structure, which can greatly improve the solving efficiency. The specific decomposition method is as follows:

[0133] 1. Define problem decomposition mapping : , decompose global problem into sub-problems, where is decomposition mapping function, is global optimization problem, is decomposed sub-problem;

[0134] 2. Ensure sub-problems maintain core structural properties of original problem, i.e. each sub-problem can also be further decomposed in similar manner, forming self-similar nested structure;

[0135] 3. Construct dependency graph G between problems, representing dependency relationships between sub-problems, G is a directed graph, nodes represent sub-problems, edges represent dependency relationships.

[0136] Hierarchical solving strategy includes the following steps:

[0137] 1. First layer: solve basic sub-problems, which are small in size and low in complexity, can be directly solved;

[0138] 2. Second layer: solve intermediate sub-problems based on first layer results, use first layer solutions as input;

[0139] 3. Third layer: integrate lower layer results to solve global problem, combine solutions of each sub-problem into global solution.

[0140] Inter-layer information transfer mechanism ensures coordination and consistency between problems at different levels. Convergence criteria can be set as relative change rate of objective function value less than given threshold (such as 0.1%) or iteration number reaching upper limit (such as 100 times).

[0141] Global coordinator adjusts parameters of each sub-problem to promote local optimization towards global optimal direction. When local optimal conflicts with global optimal, appropriate relaxation of part of constraints or adjustment of objective function weight can be made to balance relationship between different objectives.

[0142] In another embodiment of the present application, according to the output prediction data considering uncertainty and the output optimal planning sequence of each virtual power plant, a discrete-continuous hybrid dynamics model is constructed. Specifically, the output prediction data is taken as the initial condition of the system state, and the output optimal planning sequence is taken as the control target of the system, and the following hybrid dynamics model is constructed:

[0143] ,

[0144] where, is the system state vector, describing the state of the system at time t; is the system state vector at time t; is the system state vector at time t; is the control input vector, representing the control applied at time t; is the control input vector at time t; is the discrete mapping function, describing the state transition at the dispatch decision point ; is the continuous flow function, describing the dynamic evolution of the system between decision points; is the output prediction data vector; is the output optimal planning sequence vector; denotes the derivative of the state vector with respect to time; denotes the kth dispatch decision time point; denotes the k+1th dispatch decision time point; denotes that time t is within the interval between two adjacent decision points, is the output prediction data vector, containing the predicted output values of each virtual power plant at time t, is the output optimal planning sequence vector at time t.

[0145] The state variable vector can include the following components:

[0146]

[0147] where, is the system state vector, with the superscript T denoting the vector transpose, is the output of the ith virtual power plant at time t (kW), is the state of charge of the jth energy storage device at time t (dimensionless, representing the proportion of the remaining battery capacity to the total capacity), n is the number of virtual power plants, m is the number of energy storage devices, and the ellipsis indicates that other state variables such as grid frequency, voltage, etc. can also be included.

[0148] The Lyapunov exponent spectrum of the system can be calculated by the following steps:

[0149] 1. Construct the Jacobian matrix J(x, u) of the system, J is an n x n matrix, where n is the dimension of the state vector, and the matrix elements represent the partial derivative of the state variable with respect to the state derivative ;

[0150] 2. Solve the variational equation where is a small perturbation of the state, is the derivative of the perturbation, denotes the multiplication of a matrix and a vector;

[0151] 3. Calculate the expansion rate in each direction to obtain the Lyapunov exponent where denotes the size of the perturbation in the i-th direction at time t, denotes the initial perturbation size, denotes the natural logarithm, denotes the limit as time tends to infinity.

[0152] The sign of the Lyapunov exponent indicates the stability of the system in the corresponding direction: a positive value indicates instability, a negative value indicates stability, and a zero value indicates neutrality. In practical applications, it is generally desirable for the Lyapunov exponent spectrum of the system to have only a small number of positive or zero values, and most of them to be negative, so that the system has good stability.

[0153] Define the fast time scale and the slow time scale , which satisfy , i.e. the fast time scale is much smaller than the slow time scale. The overall scheduling strategy is optimized on the slow time scale, such as global optimization every hour or every day; real-time adjustment and control are performed on the fast time scale, such as local adjustment every minute or every second. Preferably, minutes, hours, so as to ensure the effect of global optimization and timely response to rapid changes in the system.

[0154] In another embodiment of the present application, the step of realizing the coordinated scheduling of multiple energy storage devices specifically includes: defining an energy storage state vector, including state of charge, charging and discharging power, and working mode; defining a state transition function to describe the change of state under action; constructing a state constraint set to ensure that the state is within a safe range; designing a hierarchical control structure, including formulating an overall scheduling strategy, allocating scheduling tasks for each energy storage device, and executing specific charging and discharging control; implementing a load sharing mechanism and fault tolerance strategy between devices to ensure that a single device failure does not affect the overall scheduling.

[0155] Specifically, the energy storage state vector can be represented as:

[0156] ,

[0157] Where s is the energy storage state vector, SoC is the state of charge, representing the energy storage capacity level (dimensionless, generally a real number between 0 and 1), P is the charging and discharging power (kW, a positive value indicates discharging, a negative value indicates charging), and mode is the working mode (an integer, such as 1 for charging, 2 for discharging, and 3 for standby).

[0158] State transition function T(s, ) Description in action Changes in the next state:

[0159] ,

[0160] in, Let T be the new state after the transition, T be the state transition function, and s be the current state. For actions to be performed (such as changing the charging and discharging power).

[0161] For example, the transition of the state of charge can be represented as:

[0162] ,

[0163] in, The state of charge after the transfer is SoC, the current state of charge is P, the charging and discharging power is kW, Δt is the time step (hours), η is the charging and discharging efficiency (ηch during charging and 1 / ηdis during discharging, both are dimensionless and range from 0 to 1), and Ecap is the energy storage capacity (kWh).

[0164] The set of state constraints, Scon, ensures that the state remains within a safe range:

[0165] ,

[0166] in, This is the set of state constraints, containing all states that satisfy the constraints. This is the lower limit of the state of charge (dimensionless). This is the upper limit of the state of charge (dimensionless). Maximum charging power (kW). The maximum discharge power (kW) is given. The set of values ​​for mode is {1,2,3}, where 1 represents charging mode, 2 represents discharging mode, and 3 represents standby mode.

[0167] In practical applications, different types of energy storage devices may have different parameter values. For example, for lithium-ion batteries, Generally, it is taken as 0.1 to 0.2. A value of 0.9 to 0.95 is generally used to protect the battery and extend its lifespan; charge / discharge efficiency and Typically between 0.9 and 0.95. For flywheel energy storage, It is desirable to approach a value of 0, with higher charge and discharge efficiency, up to 0.95 or more.

[0168] The hierarchical control structure includes three levels:

[0169] 1. Top level: Formulate overall scheduling strategy, determine the operation target of each energy storage device;

[0170] 2. Middle level: Assign scheduling tasks to each energy storage device, task allocation based on device characteristics and operating state;

[0171] 3. Bottom level: Execute specific charge and discharge control, achieve the set power curve.

[0172] The load sharing mechanism between devices can be designed based on the capacity, efficiency and life characteristics of each device. For example, a weighted allocation method can be used:

[0173]

[0174] where, Pi is the allocated power of the i-th energy storage device (kW), wi is its weight (dimensionless, representing the relative importance of the device), Σwi is the sum of the weights of all m energy storage devices, Ptotal is the total demand power (kW), wi / Σwi represents the weight proportion of the i-th device. The weight wi can be determined according to the capacity, current state of charge and life characteristics of the device, for example:

[0175]

[0176] where, wi is the weight of the i-th energy storage device, Ci is the device capacity (kWh), f(SOCi) is a function related to the state of charge, g(Li) is a function related to the device life. The design of functions f and g should ensure that when the state of charge approaches the boundary or the remaining cycle times are few, the weight is correspondingly reduced, thereby reducing the burden on the device.

[0177] Fault tolerance strategy ensures that a single device failure does not affect the overall scheduling. When a device failure is detected, the system will immediately isolate the faulty device and redistribute tasks to other normally operating devices. At the same time, the system will assess the impact of the failure on the overall scheduling and adjust the scheduling plan if necessary.

[0178] ​​In another embodiment of the present application, the step of processing the non-integer multiple of time difference of different scheduling periods specifically includes: defining the time deviation as the difference between the integer multiple of the starting time of the second scheduling period and the starting time of the first scheduling period; designing a time interpolation function to generate a transition state at a non-integer multiple of time; implementing a state smoothing transition strategy to avoid scheduling faults; constructing a state reconstruction method to maintain the continuity of scheduling; when the difference between the starting time of the second scheduling period and the starting time of the first scheduling period is not an integer multiple, calculating the charging and discharging scheduling operation of the second energy storage battery and executing the operation.

[0179] Specifically, the time deviation can be represented as:

[0180]

[0181] wherein, is the time deviation (in hours or minutes), is the starting time of the first scheduling period, is the starting time of the second scheduling period, and k is an integer representing the number of times the first scheduling period can be completely placed before the starting time of the second scheduling period. When , the non-integer multiple of time difference needs to be processed.

[0182] Time interpolation function is used to generate a transition state at a non-integer multiple of time. For example, linear interpolation can be used:

[0183]

[0184] wherein, is the system state at the starting time of the second scheduling period, is the system state at the end of the kth first scheduling period, is the system state at the end of the k+1th first scheduling period, represents the interpolation ratio (dimensionless, taking a value between 0 and 1), and are the weights of the two adjacent states, respectively.

[0185] The state smoothing transition strategy ensures that the system state changes continuously in time, avoiding scheduling faults caused by non-integer multiple of time difference. The specific implementation can use the sliding window method, i.e., at the beginning of each scheduling period, based on the current system state and historical scheduling plan, a smooth scheduling curve for a future period of time is generated.

[0186] The state reconstruction method is used to maintain the continuity of scheduling. When a non-integer multiple of time difference is detected, the system re-evaluates the current state and generates a new scheduling plan based on the latest state. This dynamic adjustment mechanism ensures the continuity and adaptability of the scheduling.​​

[0187] When the difference between the second scheduling period start time and the first scheduling period start time is not an integer multiple of zero, the charge-discharge scheduling operation of the second energy storage battery needs to be calculated. The specific steps are as follows:

[0188] 1. Obtain the charge-discharge scheduling plan of the first energy storage battery;

[0189] 2. Calculate the time deviation ;

[0190] 3. Generate the transition state of the second scheduling period start time based on the time interpolation function;

[0191] 4. According to the transition state and the output optimal planning sequence of the second scheduling period, calculate the charge-discharge scheduling operation of the second energy storage battery;

[0192] 5. Perform the calculated charge-discharge scheduling operation.

[0193] This processing mechanism solves the scheduling fault problem caused by the non-integer multiple of the time difference in the traditional scheduling method, and improves the continuity and stability of the system.

[0194] With reference Figure 4 , the application also provides a virtual power plant low-carbon economic scheduling system, comprising an output prediction module 1, an optimal planning solution module 2 and an optimization scheduling module 3.

[0195] The output prediction module 1 is used to obtain the historical output data of each virtual power plant in the power grid within a set scheduling period, predict the output data of each virtual power plant within the set scheduling period according to the time sequence, and obtain the output prediction data considering uncertainty.

[0196] The optimal planning solution module 2 is used to construct a carbon emission reduction responsibility allocation space based on the output prediction data considering uncertainty, map each virtual power plant as a vector point in the space, determine the carbon emission reduction responsibility weight of each virtual power plant through topological homeomorphism mapping, solve the target function with the lowest carbon emission cost combined with the carbon emission constraint index, and obtain the output optimal planning sequence of each virtual power plant within the set scheduling period.

[0197] The optimization scheduling module 3 is used to construct a discrete-continuous hybrid dynamics model according to the output prediction data considering uncertainty and the output optimal planning sequence of each virtual power plant, realize the collaborative scheduling of multiple energy storage devices, handle the non-integer multiple of the time difference of different scheduling periods, optimize the scheduling of each virtual power plant, and optimize the power grid.

[0198] Preferably, the system further comprises a data collection and preprocessing module 4 and a system control center 5. The data collection and preprocessing module 4 is used to collect virtual power plant real-time operation data, historical output data and external environment data, perform data cleaning, anomaly detection and missing value processing, and distribute the processed data to the output prediction module 1 and the system control center 5. The system control center 5 is used to coordinate the operation of each functional module, allocate carbon emission reduction responsibility weights based on topology mapping, monitor the overall operation state of the system, and handle abnormal situations and emergency responses.

[0199] In the system, the data interaction process between modules is as follows:

[0200] First, the data collection and preprocessing module 4 collects various data and performs preprocessing, and then distributes the processed data to the output prediction module 1 and the system control center 5.

[0201] Second, the output prediction module 1 executes a prediction algorithm based on historical data and external environment data, generates output prediction data considering uncertainty, and sends the prediction results to the optimal planning solution module 2.

[0202] Then, the optimal planning solution module 2 combines the output prediction data and the carbon emission reduction responsibility weights allocated by the system control center 5, constructs and solves the optimization problem, obtains the output optimal planning sequence of each virtual power plant, and sends the optimization results to the optimal scheduling module 3.

[0203] Finally, the optimal scheduling module 3 generates specific device scheduling instructions according to the output optimal planning sequence, handles the case of non-integer multiples of time difference values, coordinates the charging and discharging scheduling of multiple energy storage devices, and issues the scheduling instructions to the control systems of each virtual power plant.

[0204] Through the collaborative work of the above modules, the system can achieve the economic optimal operation of the virtual power plant group under low-carbon constraints, and provide strong support for the low-carbon transformation of the power system.

[0205] The above is only a specific embodiment of the present application, but the protection scope of the present application is not limited thereto, any skilled person in the art can easily think of changes or replacements within the technical scope disclosed by the present application, which should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A low-carbon economic dispatch method for virtual power plants, characterized in that, include: Historical power output data of each virtual power plant in the power grid during a set scheduling period are obtained. Based on the time series, the power output data of each virtual power plant during the set scheduling period are predicted to obtain power output prediction data that takes uncertainty into account. Based on output prediction data that takes into account uncertainties, a carbon emission reduction responsibility allocation space is constructed, and each virtual power plant is mapped to a vector point in this space. The carbon emission reduction responsibility weight of each virtual power plant is determined by topological homeomorphism mapping. The determined carbon emission reduction responsibility weight is substituted into the carbon emission constraint index to construct an optimization problem with the goal of minimizing carbon emission costs. Solving this optimization problem yields the optimal output planning sequence of each virtual power plant within the set scheduling period. Based on output forecast data considering uncertainties and the optimal output planning sequence of each virtual power plant, a discrete-continuous hybrid power dynamics model is constructed to realize the coordinated scheduling of multiple energy storage devices, and to handle the case where the time difference between different scheduling periods is not an integer multiple, thereby optimizing the scheduling of each virtual power plant and optimizing the power grid. The specific steps for determining the carbon emission reduction responsibility weights of each virtual power plant through topological homeomorphism mapping include: defining a carbon quota space, where each point represents a quota allocation scheme; constructing a mapping function from the carbon quota space to the responsibility space, ensuring the mapping satisfies continuity and reversibility; designing the mapping function according to the principles of weight guidance, quota adaptation, fairness, and efficiency; real-time monitoring of carbon dioxide intensity and grid operation status variables, distributing environmental change information to each virtual power plant through a parameter propagation network; defining the adjustment step size and adjustment direction vector to achieve iterative adjustment of responsibility allocation, and ensuring that the adjusted responsibility allocation remains within the feasible domain through constrained projection. Define the carbon quota space Q, where point This indicates the quota allocation scheme. Construct a mapping function for the carbon quota of the i-th virtual power plant. To map quota allocation to responsibility allocation, the design of the mapping function must meet the following conditions: Continuity condition: A small change in the quota leads to a small change in the responsibility, that is, for any... ,exist , making when Sometimes, ;in Indicates two points in the quota space and The distance between them Represents the corresponding point in the responsibility space. and The distance between them and It is a positive real number; Reversibility: Different quota allocations correspond to different responsibility allocations, that is, for any... ,have ;in and For the differences in the quota space, and For the corresponding point in the responsibility space; The mapping function can be designed in the following specific form: , in, Let i be the carbon emission reduction responsibility value for the i-th virtual power plant. For its carbon quota, Weight their carbon emission reduction responsibility. Total carbon emission reduction responsibility and As a weighting factor, satisfying , , This represents the total carbon allowance for all virtual power plants. This represents the proportion of the carbon allowance for the i-th virtual power plant to the total allowance; To enable dynamic adjustment of responsibility allocation, it is necessary to monitor carbon dioxide intensity in real time. and power grid operating state variables The environmental change information is distributed to each virtual power plant through a parameter propagation network. The iterative adjustment process of responsibility allocation can be represented as follows: , in, Let be the responsibility assignment vector at time t, containing the responsibility values ​​of all virtual power plants. for The responsibility allocation vector at any given moment To adjust the step size, it is a non-negative scalar. The direction vector is used to adjust the direction of responsibility adjustment. For time step; To ensure that the adjusted responsibility assignment remains within the feasible region D, constraint projection is required: , in, The responsibility assignment vector after projection. This represents the projection operation onto the convex polyhedron D, ensuring that the projected points lie within the feasible region and satisfy all constraints.

2. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific steps for constructing a carbon emission reduction responsibility allocation space based on output forecast data that takes into account uncertainties include: Construct an n-dimensional responsibility vector space, where n is the number of virtual power plants, and each virtual power plant is represented as a vector point in this space; By introducing total responsibility constraints, nonnegativity constraints, and upper limit constraints, a convex polyhedron is formed in the responsibility space, representing all feasible responsibility allocation schemes. By introducing time parameters into the responsibility space, the responsibility allocation is dynamically adjusted over time, and a smooth transition in responsibility allocation is maintained through continuous deformation.

3. The virtual power plant low-carbon economic dispatch method according to claim 2, characterized in that, The specific steps for determining the carbon emission reduction responsibility weights of each virtual power plant through topological homeomorphism mapping include: Define a carbon quota space, where each point represents a quota allocation scheme; Construct a mapping function from the carbon quota space to the responsibility space, such that the mapping satisfies the continuity condition and reversibility; The design of the mapping function follows the principles of weight guidance, quota adaptation, fairness, and efficiency. Real-time monitoring of carbon dioxide intensity and power grid operation status variables; distribution of environmental change information to each virtual power plant through parameter propagation network. Define the adjustment step size and adjustment direction vector to achieve iterative adjustment of responsibility allocation, and ensure that the adjusted responsibility allocation is still within the feasible region through constraint projection.

4. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific steps to obtain output forecast data that takes uncertainty into account include: The output of each virtual power plant at time t is modeled as a random variable, and a stochastic process is constructed to describe the temporal evolution of the output. Define the output state space and probability measure to form a complete probability space; Analyze the statistical characteristics of historical output time series to identify periodic, trend, and random components; Construct a power distribution model, including marginal distribution and conditional distribution; The predicted point value is calculated based on the conditional expectation, and the uncertainty of the prediction is quantified to construct the prediction interval.

5. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific carbon emission constraints include: Carbon emissions are matched with responsibilities to ensure that the difference between the total carbon emissions of each virtual power plant and its carbon reduction responsibility value does not exceed the set threshold. Total emissions constraints ensure that the total system emissions do not exceed the upper limit; Equipment operation constraints include upper and lower limits of output and ramp rate constraints; Specific constraints on energy storage devices include charge / discharge rates and state of charge range; By introducing a perturbation set to represent all possible prediction errors, robust constraints are constructed to ensure that the constraints are still satisfied under conditions of uncertainty.

6. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific steps for solving the objective function that minimizes carbon emission costs include: Construct an objective function f, taking into account factors such as the discount price of energy storage batteries, dispatch period, energy storage efficiency, energy consumption efficiency, discount rate and carbon price; The overall optimization problem is decomposed into a network of self-similar sub-problems, forming a fractal structure; Design a hierarchical solution strategy, including solving basic subproblems, intermediate subproblems, and global problems from the bottom up; Implement an inter-layer information transfer mechanism and design a convergence judgment criterion to control the iteration process; The global coordinator oversees the overall optimization process, resolving conflicts between local optima and global optima.

7. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific steps for constructing a discrete-continuous hybrid power dynamics model based on output forecast data considering uncertainties and the optimal output planning sequence of each virtual power plant include: ... As the initial condition for the system state, the optimal output planning sequence will be used. As the system control objective, the discrete part is designed to describe the state transitions at scheduling decision points; the continuous part is designed to describe the dynamic evolution of the system between decision points; a state variable vector is constructed to represent the system state; the Lyapunov exponent spectrum of the system is calculated to analyze the structural characteristics of the chaotic attractor; fast and slow time scales are defined, the overall scheduling strategy is optimized on the slow time scale, and real-time adjustments and control are performed on the fast time scale.

8. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific steps to achieve coordinated scheduling of multiple energy storage devices include: Define the energy storage state vector, including state of charge, charge / discharge power, and operating mode; Define a state transition function to describe the state change under an action; Construct a set of state constraints to ensure that the state remains within a safe range; Design a hierarchical control structure, including formulating an overall scheduling strategy, allocating scheduling tasks to each energy storage device, and executing specific charging and discharging control; Implement load sharing mechanisms and fault tolerance strategies among devices to ensure that the failure of a single device does not affect the overall scheduling.

9. The virtual power plant low-carbon economic dispatch method according to claim 1, characterized in that, The specific steps for handling cases where the time difference between different scheduling periods is not an integer multiple include: The time deviation is defined as the difference between the start time of the second scheduling period and the start time of the first scheduling period, which is an integer multiple thereof. Design a time interpolation function to generate transition states at non-integer multiples of time points; Implement a smooth state transition strategy to avoid scheduling gaps; Construct a state reconstruction method to maintain scheduling continuity; When the difference between the start time of the second scheduling period and the start time of the first scheduling period is not an integer multiple of zero, the charging and discharging scheduling operation of the second energy storage battery is calculated and executed.

10. A virtual power plant low-carbon economic dispatch system, used to implement the virtual power plant low-carbon economic dispatch method according to any one of claims 1-9, characterized in that, include: The power output prediction module is used to obtain historical power output data of each virtual power plant in the power grid during the set scheduling period, and predict the power output data of each virtual power plant during the set scheduling period based on the time series, so as to obtain power output prediction data that takes into account uncertainties. The optimal planning solution module is used to construct a carbon emission reduction responsibility allocation space based on output prediction data that takes into account uncertainties, map each virtual power plant to a vector point in this space, determine the carbon emission reduction responsibility weight of each virtual power plant through topological homeomorphism mapping, and solve the objective function of minimizing carbon emission cost in combination with carbon emission constraint indicators to obtain the optimal output planning sequence of each virtual power plant within the set scheduling period. The optimization scheduling module is used to construct a discrete-continuous hybrid power dynamics model based on output prediction data that takes uncertainty into account and the optimal output planning sequence of each virtual power plant. This enables the coordinated scheduling of multiple energy storage devices and handles the case where the time difference between different scheduling periods is not an integer multiple. The module optimizes the scheduling of each virtual power plant and thus optimizes the power grid.

Images