A method and system for equivalent aggregation of output power of distributed power supply clusters

Abstract

This invention relates to the field of power system operation and control technology, and particularly to an equivalent aggregation method and system for the output power of distributed generation clusters. The purpose of this invention is to construct an equivalent aggregation method and system for cluster output power. First, output models of the power of generators, photovoltaics, wind turbines, and energy storage devices in the cluster are constructed. Then, a cluster network power flow technical constraint model is constructed, and node voltage and branch power flow constraints are established based on linear branch power flow equations. Next, an inscribed ellipsoid model of the equivalent domain for cluster output power aggregation is constructed. The parameter values ​​of the inscribed ellipsoid model are solved by constructing and solving a specific optimization problem. Compared with existing methods, the proposed equivalent aggregation method for cluster output power significantly accelerates the calculation speed and greatly improves the accuracy of the solution results, enabling rapid equivalent aggregation of cluster output power and possessing significant value in practical applications.

Description

Technical Field

[0001] This invention relates to the field of power system operation and control technology, and in particular to an equivalent aggregation method and system for the output power of a distributed power generation cluster. Background Technology

[0002] For the scheduling of distributed generation clusters, it is often necessary to know the adjustable range of the cluster's power before scheduling, so that the upper-level power grid can formulate a reasonable and feasible scheduling plan. Therefore, it is essential to solve the equivalent aggregation model of distributed generation clusters. In existing technologies, a cluster with multiple types of power generation equipment has a large number of devices with different characteristics. At the same time, due to the characteristics of the power grid, there is a certain degree of coupling between the output power of the devices. In addition, since the devices themselves may have constraints such as energy constraints and ramping constraints, the coupling relationships between the decision variables of the cluster are relatively tight. Solving the equivalent aggregation model involves high dimensionality, making it quite difficult. Summary of the Invention

[0003] To address the aforementioned issues, this invention proposes an equivalent aggregation method and system for cluster output power. It employs a high-dimensional polyhedral inscribed ellipsoid method to construct an inscribed ellipsoid model of the equivalent domain for cluster output power aggregation, thereby achieving equivalent aggregation of the cluster.

[0004] An equivalent aggregation method for the output power of a distributed power supply cluster, the method comprising:

[0005] Construct output models of the power of each device in the distributed power cluster to obtain the active power and reactive power of each device;

[0006] Based on the active and reactive power of each device, a power flow technical constraint model for the cluster network is constructed.

[0007] Based on the power flow constraint model of the cluster network, an inscribed ellipsoid model of the equivalent domain of cluster output power aggregation is constructed, and the parameter values ​​of the inscribed ellipsoid model are solved. The inscribed ellipsoid model serves as a constraint condition to characterize the feasible domain of equivalent aggregation of the output power of the distributed power cluster.

[0008] Furthermore, the construction of the power output model for each device in the cluster includes constructing a synchronous generator device model, a photovoltaic device model, a wind turbine device model, and an energy storage device model.

[0009] Furthermore, the construction of the synchronous generator equipment model specifically includes:

[0010] The following model is constructed based on the parameters of the synchronous generator:

[0011]

[0012]

[0013]

[0014] Equation (1) represents the constraint on the active power output of generator equipment i at time t, Equation (2) represents the constraint on the reactive power output of generator equipment i at time t, and Equation (3) represents the ramp constraint of generator equipment i from time t-1 to time t. These represent the active and reactive power output at time t of the generator device labeled i; These represent the minimum and maximum active power output of the generator device labeled i, respectively. These represent the minimum and maximum reactive power output of the generator device labeled i, respectively. Δt represents the downhill and uphill climbing parameters of the generator device labeled i, respectively; Δt represents the time interval between two adjacent moments.

[0015] Furthermore, the construction of the photovoltaic equipment model specifically includes:

[0016] The following model is constructed based on the parameters of the photovoltaic equipment:

[0017]

[0018]

[0019]

[0020] Equation (4) represents the constraint on the active power output of the photovoltaic device labeled j at time t, Equation (5) represents the constraint on the reactive power output of the photovoltaic device labeled j at time t, and Equation (6) represents the power capacity constraint of the photovoltaic device labeled j at time t. These represent the active and reactive power output at time t of the photovoltaic device labeled j; These represent the minimum and maximum active power output of the photovoltaic device labeled j, respectively. These represent the minimum and maximum reactive power output of the photovoltaic device labeled j, respectively. This indicates the maximum capacity of the photovoltaic device labeled j.

[0021] Furthermore, the construction of the wind turbine equipment model specifically includes:

[0022] The following model is constructed based on the parameters of the wind turbine equipment:

[0023]

[0024]

[0025]

[0026] Equation (7) represents the constraint on the active power output of the wind turbine equipment with the label k at time t, Equation (8) represents the constraint on the reactive power output of the wind turbine equipment with the label k at time t, and Equation (9) represents the power capacity constraint of the wind turbine equipment with the label k at time t. These represent the active and reactive power output of the wind turbine device labeled k at time t, respectively. These represent the minimum and maximum active power output of the fan unit labeled k, respectively. These represent the minimum and maximum reactive power output of the fan equipment labeled k, respectively; This indicates the maximum capacity of the fan equipment labeled k.

[0027] Furthermore, the construction of the energy storage device model specifically includes:

[0028] The following model is constructed based on the parameters of the energy storage device:

[0029]

[0030]

[0031]

[0032]

[0033] Equation (10) represents the constraint on the active power output of the energy storage device with label l at time t, Equation (11) represents the constraint on the reactive power output of the energy storage device with label l at time t, and Equations (12)-(13) represent the energy constraints of the energy storage device with label l at time t. These represent the active and reactive power output at time t of the energy storage device labeled l; These represent the minimum and maximum active power output of the energy storage device labeled l, respectively. These represent the minimum and maximum reactive power output of the energy storage device labeled l, respectively. This represents the stored energy of the energy storage device labeled l at time t; These represent the minimum and maximum energy storage values ​​of the energy storage device labeled l.

[0034] Furthermore, the construction of the cluster network power flow technical constraint model includes: defining cluster nodes and injecting active and reactive power; calculating the linear branch power flow equations based on the active and reactive power injected into the nodes; and constructing node voltage and branch power flow constraint conditions by combining the linear branch power flow equations.

[0035] Furthermore, defining cluster nodes and injecting active and reactive power includes:

[0036] The cluster nodes include generators, photovoltaics, wind turbines, and energy storage devices in the cluster. Based on the active and reactive power of the generators, photovoltaics, wind turbines, and energy storage devices, the active and reactive power injected into the nodes are defined as follows:

[0037]

[0038]

[0039] Equation (14) represents the active power injected into node n at time t, and Equation (15) represents the reactive power injected into node n at time t. These represent the active and reactive power injected at node n at time t, respectively. These represent the active and reactive loads at node n at time t, respectively. Let each represent a set of numbers for the generator, photovoltaic, wind turbine, and energy storage devices at node n. These represent the active and reactive power output at time t of the generator device labeled i; These represent the active and reactive power output at time t of the photovoltaic device labeled j; These represent the active and reactive power output at time t of the wind turbine device labeled k; These represent the active and reactive power output of the energy storage device labeled l at time t, respectively.

[0040] Furthermore, the step of calculating the linear branch power flow equation based on the active and reactive power injected by the nodes includes:

[0041] Based on the pre-established decoupled linear power flow model, and according to the active and reactive power injected at each node, the following vectors are obtained: the vector composed of the voltage amplitudes of each node, the vector composed of the active power flow of each branch, and the vector composed of the reactive power flow:

[0042]

[0043]

[0044]

[0045] Wherein, equation (16) represents the voltage amplitude at time t, equation (17) represents the active power of each branch at time t, and equation (18) represents the reactive power of each branch at time t; A V A p A q B V B pB q It is the network parameter matrix, b V b p b q This is the network parameter vector, and the above network parameter values ​​are all calculated based on the decoupled linear power flow model; These represent the vectors consisting of the active and reactive power injected into each node at time t; v t This represents the vector composed of the voltage amplitudes of each node at time t; These represent the vectors composed of the active and reactive power flow of each branch at time t.

[0046] Furthermore, the construction of node voltage and branch power flow constraints by combining the linear branch power flow equations includes:

[0047] The constraints on node voltages and branch power flows are as follows:

[0048]

[0049]

[0050] Equation (19) represents the voltage amplitude constraints at time t, and Equation (20) represents the branch power capacity constraints from node i to node j at time t; v, These represent the vectors formed by the minimum and maximum voltage amplitudes of each node, respectively. This represents the maximum branch capacity from node i to node j;

[0051] In equation (20), the quadratic constraint condition represents a circular feasible region constraint condition. To simplify the quadratic constraint condition into a linear constraint condition, the feasible region of an inscribed regular polygon is used to approximate the constraint condition represented by the circular feasible region. The approximate inscribed regular polygon constraint is as follows:

[0052]

[0053] Where k = 1, 2, ..., n e n e This represents the number of sides of the equivalent approximation inscribed polygon, and its value is 8.

[0054] Furthermore, the construction of the inscribed ellipsoid model of the cluster output power aggregation equal domain includes:

[0055] Construct a decision variable vector x, by Composed of; among which These represent the active and reactive power outputs of the generator device labeled i at time t, respectively. These represent the active and reactive power outputs of the photovoltaic device labeled j at time t, respectively. These represent the active and reactive power output of the wind turbine device labeled k at time t, respectively. These represent the active and reactive power outputs of the energy storage device labeled l at time t, respectively.

[0056] The constraints described in equations (14)-(19) and (21) can be rearranged into the following matrix form:

[0057] Wx≤w (22)

[0058] Among them, matrix vector The constraint parameter is represented by n, whose value is obtained by rearranging the constraint conditions described in equations (14)-(19) and (21); x n represents the length of vector x; w This represents the length of vector w;

[0059] Based on the linear power flow model described in equations (17) and (18), the relationship between the output active power and the decision variable vector at each time point of the cluster grid connection is expressed as follows:

[0060] p0 = Dx (23)

[0061] Where p0 represents the vector of active power output and decision variables at each time point of the cluster grid connection; matrix The constraint parameter is represented, and its value can be obtained from the linear power flow model described in equations (17) and (18); n p This represents the length of vector p0;

[0062] Define the inscribed ellipsoid of the cluster output power aggregation isotropic region as E.

[0063]

[0064] Where z represents the independent variable of the aggregated domain; matrix sum vector Parameters representing a high-dimensional ellipsoid. The variable represents the possible values ​​of the ellipsoid; the operator ||·|| represents the norm of a vector; p0 represents a feasible point within the ellipsoid; x * This represents a vector of decision variables that satisfy the constraints.

[0065] Furthermore, the process of solving for the parameter values ​​of the inscribed ellipsoid model includes:

[0066] definition Let D be the matrix formed by the orthogonal basis extended from matrix D; define Let n be a matrix consisting of an orthogonal basis formed by extending the null space of matrix D, where n null =n x -np

[0067] Define matrix The values ​​are shown in equations (25)-(27);

[0068] W1:=WD1 (25)

[0069] W2:=WD2 (26)

[0070]

[0071] Solve the optimization problem shown below:

[0072]

[0073] stE≥0 (29)

[0074]

[0075]

[0076] Where i = 1, 2, ..., n w ; Represents the auxiliary variable matrix; Represents an auxiliary variable vector; α i Let represent the i-th element in vector α; (·) i The symbol represents the i-th row of a matrix or the i-th element of a vector; E≥0 indicates that the constraint matrix E is a positive semi-definite matrix;

[0077] By solving the optimization problems described in equations (28), (29), (30), and (31), the values ​​of the inscribed ellipsoid model matrix E and the vector e of the cluster output power aggregation equal domain can be obtained.

[0078] An equivalent aggregation system for the output power of a distributed power supply cluster includes: a power output model construction module, a cluster network power flow technical constraint model construction module, and an inscribed ellipsoid model construction module for the equivalent domain of cluster output power aggregation.

[0079] The power output model construction module is used to construct the power output model of each device in the distributed power cluster, and obtain the active power and reactive power of each device.

[0080] The cluster network power flow technical constraint model construction module is used to construct a cluster network power flow technical constraint model based on the active and reactive power of each device.

[0081] The module for constructing the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation is used to construct the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation based on the cluster network power flow technical constraint model, and solve the parameter values ​​of the inscribed ellipsoid model; wherein, the inscribed ellipsoid model serves as a constraint condition to characterize the feasible domain of equivalent aggregation of distributed power cluster output power.

[0082] An electronic device includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus;

[0083] Memory, used to store computer programs;

[0084] When the processor executes the program stored in the memory, it implements the above-described method for equal aggregation of the output power of a distributed power cluster.

[0085] A computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-described method for equal aggregation of output power of a distributed power cluster.

[0086] The present invention has at least the following beneficial effects:

[0087] Compared with existing technologies, the equivalent aggregation model of distributed power clusters solved by this invention has a simpler calculation process, a faster algorithm, and higher accuracy, achieving precise aggregation of the adjustable power relationships of distributed power clusters.

[0088] The solution results of this invention can be used as constraints for cluster scheduling, providing decision-making reference support for the optimized scheduling of distributed power supply clusters and laying the foundation for leveraging the auxiliary services of distributed power supply clusters.

[0089] The purpose of this invention is to construct an equivalent aggregation method and system for cluster output power. First, output models of the power of generators, photovoltaics, wind turbines, and energy storage devices in the cluster are constructed. Then, a technical constraint model of the cluster network power flow is constructed, establishing node voltage and branch power flow constraints based on linear branch power flow equations. Next, an inscribed ellipsoid model of the equivalent domain for cluster output power aggregation is constructed. By constructing and solving a specific optimization problem, the parameter values ​​of the inscribed ellipsoid model are obtained. Compared with existing methods, the proposed equivalent aggregation method for cluster output power significantly accelerates computation and greatly improves the accuracy of the solution results, enabling rapid equivalent aggregation of cluster output power and possessing significant value in practical applications.

[0090] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures pointed out in the description and the drawings. Attached Figure Description

[0091] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0092] Figure 1 This is a flowchart of the equivalent aggregation method for cluster output power according to an embodiment of the present invention;

[0093] Figure 2 This is a schematic diagram of the equivalent aggregation system for cluster output power in an embodiment of the present invention. Detailed Implementation

[0094] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0095] In a cluster of power generation equipment of various types, the equipment is numerous and has diverse characteristics. Due to the characteristics of the power grid, there is a certain degree of coupling between the output power of the equipment. At the same time, since the equipment itself may have constraints such as energy constraints and ramp-up constraints, the coupling between the decision variables of the cluster is relatively tight. Solving the equivalent aggregation model involves a high dimension, making it difficult to solve.

[0096] Therefore, in the first aspect, such as Figure 1 As shown, this invention proposes an equivalent aggregation method for the output power of a distributed power supply cluster, the method comprising:

[0097] Construct an output model of the power of each device in the cluster to obtain the active power and reactive power of each device;

[0098] Based on the active and reactive power of each device, a power flow technical constraint model for the cluster network is constructed.

[0099] Based on the power flow constraint model of the cluster network, an inscribed ellipsoid model of the equivalent value domain of cluster output power aggregation is constructed, and the parameter values ​​of the inscribed ellipsoid model are solved to obtain the inscribed ellipsoid model of the equivalent value domain of cluster output power aggregation. The inscribed ellipsoid model serves as a constraint condition to characterize the feasible domain of equivalent aggregation of distributed power cluster output power.

[0100] In this embodiment, the construction of the power output model of each device in the cluster includes constructing a synchronous generator device model, a photovoltaic device model, a wind turbine device model, and an energy storage device model.

[0101] In this embodiment, the construction of the synchronous generator equipment model specifically includes:

[0102] The following model is constructed based on the parameters of the synchronous generator:

[0103]

[0104]

[0105]

[0106] Equation (1) represents the constraint on the active power output of generator equipment i at time t, Equation (2) represents the constraint on the reactive power output of generator equipment i at time t, and Equation (3) represents the ramp constraint of generator equipment i from time t-1 to time t. These represent the active and reactive power output at time t of the generator device labeled i; These represent the minimum and maximum active power output of the generator device labeled i, respectively. These represent the minimum and maximum reactive power output of the generator device labeled i, respectively. Δt represents the downhill and uphill climbing parameters of the generator device labeled i, respectively; Δt represents the time interval between two adjacent moments.

[0107] In this embodiment, the construction of the photovoltaic device model specifically includes:

[0108] The following model is constructed based on the parameters of the photovoltaic equipment:

[0109]

[0110]

[0111]

[0112] Equation (4) represents the constraint on the active power output of the photovoltaic device labeled j at time t, Equation (5) represents the constraint on the reactive power output of the photovoltaic device labeled j at time t, and Equation (6) represents the power capacity constraint of the photovoltaic device labeled j at time t. These represent the active and reactive power output at time t of the photovoltaic device labeled j; These represent the minimum and maximum active power output of the photovoltaic device labeled j, respectively. These represent the minimum and maximum reactive power output of the photovoltaic device labeled j, respectively. This indicates the maximum capacity of the photovoltaic device labeled j.

[0113] In this embodiment, the construction of the wind turbine equipment model specifically includes:

[0114] The following model is constructed based on the parameters of the wind turbine equipment:

[0115]

[0116]

[0117]

[0118] Equation (7) represents the constraint on the active power output of the wind turbine equipment with the label k at time t, Equation (8) represents the constraint on the reactive power output of the wind turbine equipment with the label k at time t, and Equation (9) represents the power capacity constraint of the wind turbine equipment with the label k at time t. These represent the active and reactive power output of the wind turbine device labeled k at time t, respectively. These represent the minimum and maximum active power output of the fan unit labeled k, respectively. These represent the minimum and maximum reactive power output of the fan equipment labeled k, respectively; This indicates the maximum capacity of the fan equipment labeled k.

[0119] In this embodiment, constructing the energy storage device model specifically includes:

[0120] The following model is constructed based on the parameters of the energy storage device:

[0121]

[0122]

[0123]

[0124]

[0125] Equation (10) represents the constraint on the active power output of the energy storage device with label l at time t, Equation (11) represents the constraint on the reactive power output of the energy storage device with label l at time t, and Equations (12)-(13) represent the energy constraints of the energy storage device with label l at time t. These represent the active and reactive power output at time t of the energy storage device labeled l; These represent the minimum and maximum active power output of the energy storage device labeled l, respectively. These represent the minimum and maximum reactive power output of the energy storage device labeled l, respectively. This represents the stored energy of the energy storage device labeled l at time t; These represent the minimum and maximum energy storage values ​​of the energy storage device labeled l.

[0126] In this embodiment, the construction of the cluster network power flow technical constraint model includes: defining cluster nodes and injecting active and reactive power; calculating the linear branch power flow equations based on the active and reactive power injected into the nodes; and constructing node voltage and branch power flow constraint conditions by combining the linear branch power flow equations.

[0127] In this embodiment, defining cluster nodes and injecting active and reactive power includes:

[0128] The cluster nodes include generators, photovoltaics, wind turbines, and energy storage devices in the cluster. Based on the active and reactive power of the generators, photovoltaics, wind turbines, and energy storage devices, the active and reactive power injected into the nodes are defined as follows:

[0129]

[0130]

[0131] Equation (14) represents the active power injected into node n at time t, and Equation (15) represents the reactive power injected into node n at time t. These represent the active and reactive power injected at node n at time t, respectively. These represent the active and reactive loads at node n at time t, respectively. Let each represent a set of numbers for the generator, photovoltaic, wind turbine, and energy storage devices at node n. These represent the active and reactive power output at time t of the generator device labeled i; These represent the active and reactive power output at time t of the photovoltaic device labeled j; These represent the active and reactive power output at time t of the wind turbine device labeled k; These represent the active and reactive power output of the energy storage device labeled l at time t, respectively.

[0132] In this embodiment, the step of calculating the linear branch power flow equation based on the active and reactive power injected by the nodes includes:

[0133] Based on the pre-established decoupled linear power flow model, and according to the active and reactive power injected at each node, the following linear branch power flow equations are obtained:

[0134]

[0135]

[0136]

[0137] Wherein, equation (16) represents the voltage amplitude at time t, equation (17) represents the active power of each branch at time t, and equation (18) represents the reactive power of each branch at time t. A V A p A q B V B p B q It is the network parameter matrix, b V b p b q This is the network parameter vector, and the above network parameter values ​​are all calculated based on the decoupled linear power flow model; These represent the vectors consisting of the active and reactive power injected into each node at time t; v t This represents the vector composed of the voltage amplitudes of each node at time t; These represent the vectors composed of the active and reactive power flow of each branch at time t.

[0138] In this embodiment, the construction of node voltage and branch power flow constraints by combining the linear branch power flow equations includes:

[0139] The constraints on node voltages and branch power flows are as follows:

[0140]

[0141]

[0142] Equation (19) represents the voltage amplitude constraint at time t, and Equation (20) represents the branch power capacity constraint from node i to node j at time t. v , These represent the vectors formed by the minimum and maximum voltage amplitudes of each node, respectively. This represents the maximum branch capacity from node i to node j;

[0143] In equation (20), the quadratic constraint condition represents a circular feasible region constraint condition. To simplify the quadratic constraint condition into a linear constraint condition, the constraint condition represented by the circular feasible region is approximated by the feasible region of an inscribed regular polygon. The approximated inscribed regular polygon constraint is as follows:

[0144]

[0145] Where k = 1, 2, ..., n e n e This represents the number of sides of the equivalent approximation inscribed polygon, and its value is usually 8.

[0146] In this embodiment, constructing the inscribed ellipsoid model of the cluster output power aggregation isotropic domain includes:

[0147] Construct a decision variable vector x, by Composed of;

[0148] The constraints described in equations (1)-(19) and (21) can be rearranged into the following matrix form:

[0149] Wx≤w (22)

[0150] Among them, matrix vector The constraint parameter is represented, and its value can be obtained by rearranging the constraint conditions described in equations (1)-(19) and (21); n x n represents the length of vector x; w This represents the length of vector w;

[0151] Based on the linear power flow model described in equations (17)-(18), the relationship between the output active power and the decision variable vector at each time point of the cluster grid connection is expressed as follows:

[0152] p0 = Dx (23)

[0153] Where p0 represents the vector of active power output and decision variables at each time point of the cluster grid connection; matrix This represents the constraint parameter, whose value can be obtained from the linear power flow model described in equations (17)-(18); n p This represents the length of vector p0.

[0154] Define the inscribed ellipsoid of the cluster output power aggregation isotropic region as E.

[0155]

[0156] Where z represents the independent variable of the aggregated domain; matrix sum vector Parameters representing a high-dimensional ellipsoid. The variable represents the possible values ​​of the ellipsoid; the operator ||·|| represents the norm of a vector; p0 represents a feasible point within the ellipsoid; x * This represents a vector of decision variables that satisfy the constraints.

[0157] In this embodiment, solving for the parameter values ​​of the inscribed ellipsoid model includes:

[0158] definition Let D be the matrix formed by the orthogonal basis extended from matrix D; define Let n be a matrix consisting of an orthogonal basis formed by extending the null space of matrix D, where n null =n x -n p

[0159] Define matrix The values ​​are shown in equations (25)-(27);

[0160] W1:=WD1 (25)

[0161] W2:=WD2 (26)

[0162]

[0163] Solve the optimization problem shown below:

[0164]

[0165] stE≥0 (29)

[0166]

[0167]

[0168] Where i = 1, 2, ..., n w ; Represents the auxiliary variable matrix; Represents an auxiliary variable vector; α i Let represent the i-th element in vector α; (·) i The symbol represents the i-th row of a matrix or the i-th element of a vector; E≥0 indicates that the constraint matrix E is a positive semi-definite matrix;

[0169] By solving the optimization problems described in equations (28)-(31), the values ​​of the inscribed ellipsoid model matrix E and vector e of the cluster output power aggregation equal domain can be obtained.

[0170] Secondly, such as Figure 2As shown, this invention proposes an equivalent aggregation system for the output power of a distributed power supply cluster, comprising: a power output model construction module, a cluster network power flow technical constraint model construction module, and an inscribed ellipsoid model construction module for the equivalent domain of cluster output power aggregation.

[0171] The power output model building module is used to build the power output model of each device in the cluster, and obtain the active power and reactive power of each device.

[0172] The cluster network power flow technical constraint model construction module is used to construct a cluster network power flow technical constraint model based on the active and reactive power of each device.

[0173] The module for constructing the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation is used to construct the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation based on the power flow technical constraint model of the cluster network, and solve the parameter values ​​of the inscribed ellipsoid model to obtain the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation. The inscribed ellipsoid model serves as a constraint condition to characterize the feasible domain of equivalent aggregation of the output power of the distributed power cluster.

[0174] In practice, the implementation processes of the equivalent aggregation system for the output power of a distributed power cluster and the equivalent aggregation method for the output power of a distributed power cluster are one-to-one, and will not be elaborated here.

[0175] Thirdly, the present invention provides an electronic device, including a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus;

[0176] Memory, used to store computer programs;

[0177] When the processor executes the program stored in the memory, it implements the above-described method for equal aggregation of the output power of a distributed power cluster.

[0178] Fourthly, a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above-described method for equal aggregation of output power of a distributed power cluster.

[0179] This invention first constructs power constraint models for four types of equipment: synchronous generators, photovoltaics, wind turbines, and energy storage. Then, it constructs the network power flow equations and branch power flow constraints. Finally, it solves the inscribed ellipsoid model of the aggregated power domain of the cluster output.

[0180] The solution results of this invention can be used as constraints for cluster scheduling, providing decision-making reference support for the optimized scheduling of distributed power supply clusters and laying the foundation for leveraging the auxiliary services of distributed power supply clusters.

[0181] Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for equivalent aggregation of output power of a distributed power supply cluster, characterized in that, The method includes: Construct output models of the power of each device in the distributed power cluster to obtain the active power and reactive power of each device; Based on the active and reactive power of each device, a power flow technical constraint model for the cluster network is constructed. Based on the power flow constraint model of the cluster network, an inscribed ellipsoid model of the equivalent domain of cluster output power aggregation is constructed, and the parameter values ​​of the inscribed ellipsoid model are solved. The inscribed ellipsoid model serves as a constraint condition to characterize the feasible domain of equivalent aggregation of the output power of the distributed power cluster. The construction of output models for the power of each device in the distributed power cluster includes constructing output models for synchronous generator devices, photovoltaic devices, wind turbine devices, and energy storage devices. The construction of the cluster network power flow constraint model includes: defining cluster nodes and injecting active and reactive power; calculating the linear branch power flow equations based on the active and reactive power injected into the nodes; and constructing node voltage and branch power flow constraint conditions by combining the linear branch power flow equations. The parameter values ​​for solving the inscribed ellipsoid model include: definition For the matrix The matrix formed by the extended orthogonal basis; definition For the matrix The matrix formed by extending the null space into an orthogonal basis, where, ; Representing vectors Length; Represents a vector of decision variables; Representing vectors Length; This represents the output active power and decision variable vector at each time point of the cluster grid connection; Define matrix , , The values ​​are shown in equations (25)-(27); (25) (26) (27) Solve the optimization problem shown below: (28) (29) (30) (31) in, ; Represents the auxiliary variable matrix; , Represents an auxiliary variable vector; Representing vectors The first in One element; The symbol represents the first of the matrix. The row or the first row of a vector One element; Representing the constraint matrix It is a positive semi-definite matrix; Representing vectors The length of the vector; Indicates constraint parameters; By solving the optimization problems described in equations (28), (29), (30), and (31), the inscribed ellipsoid model matrix of the cluster output power aggregation equal domain is obtained. sum vector The value of .

2. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 1, characterized in that, The construction of the output model for the synchronous generator equipment specifically includes: The following model is constructed based on the parameters of the synchronous generator equipment: (1) (2) (3) In this context, equation (1) represents the constraint on the active power output of generator equipment i at time t, equation (2) represents the constraint on the reactive power output of generator equipment i at time t, and equation (3) represents the ramp-up constraint of generator equipment i from time t-1 to time t. , These represent the active and reactive power output at time t of the generator device labeled i; , These represent the minimum and maximum active power output of the generator device labeled i, respectively. , These represent the minimum and maximum reactive power output of the generator device labeled i, respectively. , These represent the downhill and uphill ramp parameters for the generator device labeled i, respectively. It represents the time interval between two adjacent moments.

3. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 1, characterized in that, The output model for constructing the photovoltaic device specifically includes: The following model is constructed based on the parameters of the photovoltaic equipment: (4) (5) (6) Equation (4) represents the constraint on the active power output of the photovoltaic device with label j at time t, Equation (5) represents the constraint on the reactive power output of the photovoltaic device with label j at time t, and Equation (6) represents the power capacity constraint of the photovoltaic device with label j at time t. , These represent the active and reactive power output at time t of the photovoltaic device labeled j; , These represent the minimum and maximum active power output of the photovoltaic device labeled j, respectively. , These represent the minimum and maximum reactive power output of the photovoltaic device labeled j, respectively. This indicates the maximum capacity of the photovoltaic device labeled j.

4. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 1, characterized in that, The output model for constructing the wind turbine equipment specifically includes: The following model is constructed based on the parameters of the wind turbine equipment: (7) (8) (9) Equation (7) represents the constraint on the active power output of the wind turbine equipment with the label k at time t, Equation (8) represents the constraint on the reactive power output of the wind turbine equipment with the label k at time t, and Equation (9) represents the power capacity constraint of the wind turbine equipment with the label k at time t. , These represent the active and reactive power output at time t of the wind turbine device labeled k; , These represent the minimum and maximum active power output of the fan unit labeled k, respectively. , These represent the minimum and maximum reactive power output of the fan unit labeled k, respectively. This indicates the maximum capacity of the fan equipment labeled k.

5. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 1, characterized in that, The output model for constructing the energy storage device specifically includes: The following model is constructed based on the parameters of the energy storage device: (10) (11) (12) (13) Equation (10) represents the constraint on the active power output of the energy storage device with label l at time t, Equation (11) represents the constraint on the reactive power output of the energy storage device with label l at time t, and Equations (12)-(13) represent the energy constraints of the energy storage device with label l at time t. , These represent the active and reactive power output at time t of the energy storage device labeled l; , These represent the minimum and maximum active power output of the energy storage device labeled l, respectively. , These represent the minimum and maximum reactive power output of the energy storage device labeled l, respectively. This represents the stored energy of the energy storage device labeled l at time t; , These represent the minimum and maximum energy storage values ​​of the energy storage device labeled 'l', respectively. It represents the time interval between two adjacent moments.

6. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 1, characterized in that, The process of defining cluster nodes and injecting active and reactive power includes: The cluster nodes include generators, photovoltaics, wind turbines, and energy storage devices in the cluster. Based on the active and reactive power of the generators, photovoltaics, wind turbines, and energy storage devices in the cluster, the active and reactive power injected into the nodes are defined as follows: (14) (15) Equation (14) represents the active power injected into node n at time t, and Equation (15) represents the reactive power injected into node n at time t. , These represent the active and reactive power injected at node n at time t, respectively. , These represent the active and reactive loads at node n at time t, respectively. , , , Let each represent a set of numbers for the generator, photovoltaic, wind turbine, and energy storage devices at node n. , These represent the active and reactive power output at time t of the generator device labeled i; , These represent the active and reactive power output at time t of the photovoltaic device labeled j; , These represent the active and reactive power output at time t of the wind turbine device labeled k; , These represent the active and reactive power output of the energy storage device labeled l at time t, respectively.

7. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 6, characterized in that, The calculation of the linear branch power flow equations based on the active and reactive power injected by the nodes includes: Based on the pre-established decoupled linear power flow model, and according to the active and reactive power injected at each node, the following vectors are obtained: the vector composed of the voltage magnitudes of each node, the vector composed of the active power flow of each branch, and the vector composed of the reactive power flow: (16) (17) (18) In this context, equation (16) represents the voltage amplitude at time t, equation (17) represents the active power of each branch at time t, and equation (18) represents the reactive power of each branch at time t. , , , , , It is the network parameter matrix. , , It is the network parameter vector. , , , , , , , , All were calculated based on the decoupled linear power flow model; , These represent the vectors composed of the active and reactive power injected into each node at time t, respectively. This represents the vector composed of the voltage amplitudes of each node at time t; , These represent the vectors composed of the active and reactive power flow of each branch at time t.

8. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 7, characterized in that, The combined linear branch power flow equations are used to construct node voltage and branch power flow constraints, including: The constraints on node voltages and branch power flows are as follows: (19) (20) Equation (19) represents the voltage amplitude constraint at time t, and Equation (20) represents the branch power capacity constraint from node i to node j at time t. , These represent the vectors formed by the minimum and maximum voltage amplitudes of each node, respectively. This represents the maximum branch capacity from node i to node j; Wherein, the quadratic constraint condition in equation (20) represents a circular feasible region constraint condition; the constraint condition represented by the circular feasible region is approximated by the feasible region of an inscribed regular polygon, and the approximated inscribed regular polygon constraint is as follows: (21) in, , This represents the number of sides of the equivalent approximation inscribed polygon, and its value is 8.

9. The method for equivalent aggregation of output power of a distributed power supply cluster according to claim 8, characterized in that, The construction of the inscribed ellipsoid model of the cluster output power aggregation equal domain includes: Construct decision variable vector ,Depend on , , , , , , , Composed of; among which , These represent the active and reactive power outputs of the generator device labeled i at time t, respectively. , These represent the active and reactive power outputs of the photovoltaic device labeled j at time t, respectively. , These represent the active and reactive power output of the wind turbine device labeled k at time t, respectively. , These represent the active and reactive power outputs of the energy storage device labeled l at time t, respectively. The constraints described in equations (14)-(19) and (21) can be rearranged into the following matrix form: (22) Among them, matrix ,vector The constraint parameter is represented by the value obtained by organizing the constraint conditions described in equations (14)-(19) and (21); Representing vectors Length; Representing vectors Length; Based on the linear power flow model described in equations (17) and (18), the relationship between the output active power and the decision variable vector at each time point of the cluster grid connection is expressed as follows: (23) in, This represents the vector of active power output and decision variables at each time point of the cluster grid connection; matrix The constraint parameter is represented by the value obtained according to the linear power flow model described in equations (17) and (18); Representing vectors Length; Define the inscribed ellipsoid of the cluster output power aggregation isotropic region as (24) in, The independent variable representing the domain of aggregation; matrix sum vector Parameters representing a high-dimensional ellipsoid. The variable representing the possible values ​​of the ellipsoid; operators The norm of a vector; This represents a feasible point within the ellipsoid region; This represents a vector of decision variables that satisfy the constraints.

10. An equivalent aggregation system for the output power of a distributed power supply cluster, characterized in that, include: Power output model building module, cluster network power flow technical constraint model building module, and cluster output power aggregation equal domain inscribed ellipsoid model building module; The power output model construction module is used to construct the power output model of each device in the distributed power cluster, and obtain the active power and reactive power of each device. The cluster network power flow technical constraint model construction module is used to construct a cluster network power flow technical constraint model based on the active and reactive power of each device. The module for constructing the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation is used to construct the inscribed ellipsoid model of the equivalent domain of cluster output power aggregation based on the cluster network power flow technical constraint model, and solve for the parameter values ​​of the inscribed ellipsoid model; wherein, the inscribed ellipsoid model serves as a constraint condition to characterize the feasible domain of equivalent aggregation of distributed power cluster output power. The construction of output models for the power of each device in the distributed power cluster includes constructing output models for synchronous generator devices, photovoltaic devices, wind turbine devices, and energy storage devices. The construction of the cluster network power flow constraint model includes: defining cluster nodes and injecting active and reactive power; calculating the linear branch power flow equations based on the active and reactive power injected into the nodes; and constructing node voltage and branch power flow constraint conditions by combining the linear branch power flow equations. The parameter values ​​for solving the inscribed ellipsoid model include: definition For the matrix The matrix formed by the extended orthogonal basis; definition For the matrix The matrix formed by extending the null space into an orthogonal basis, where, ; Representing vectors Length; Represents a vector of decision variables; Representing vectors Length; This represents the output active power and decision variable vector at each time point of the cluster grid connection; Define matrix , , The values ​​are shown in equations (25)-(27); (25) (26) (27) Solve the optimization problem shown below: (28) (29) (30) (31) in, ; Represents the auxiliary variable matrix; , Represents an auxiliary variable vector; Representing vectors The first in One element; The symbol represents the first of the matrix. The row or the first row of a vector One element; Representing the constraint matrix It is a positive semi-definite matrix; Representing vectors The length of the vector; Indicates constraint parameters; By solving the optimization problems described in equations (28), (29), (30), and (31), the inscribed ellipsoid model matrix of the cluster output power aggregation equal domain is obtained. sum vector The value of .

11. An electronic device, characterized in that, It includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; When a processor executes a program stored in memory, it implements an equivalent aggregation method for the output power of a distributed power cluster as described in any one of claims 1 to 9.

12. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements an equal-value aggregation method for the output power of a distributed power cluster as described in any one of claims 1 to 9.

Images