CN114925880A - Virtual energy storage power plant distributed cooperation method based on non-precise alternative direction multiplier method - Google Patents

Virtual energy storage power plant distributed cooperation method based on non-precise alternative direction multiplier method Download PDF

Info

Publication number
CN114925880A
CN114925880A CN202210377568.XA CN202210377568A CN114925880A CN 114925880 A CN114925880 A CN 114925880A CN 202210377568 A CN202210377568 A CN 202210377568A CN 114925880 A CN114925880 A CN 114925880A
Authority
CN
China
Prior art keywords
energy storage
power
distributed
power plant
virtual
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210377568.XA
Other languages
Chinese (zh)
Other versions
CN114925880B (en
Inventor
蔡德福
王文娜
俞德华
吕莎
王枭
周鲲鹏
刘海光
陈汝斯
万黎
王涛
张良一
孙冠群
王尔玺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Wuhan University WHU
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Wuhan University WHU
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Wuhan University WHU, Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202210377568.XA priority Critical patent/CN114925880B/en
Priority claimed from CN202210377568.XA external-priority patent/CN114925880B/en
Publication of CN114925880A publication Critical patent/CN114925880A/en
Application granted granted Critical
Publication of CN114925880B publication Critical patent/CN114925880B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a distributed cooperation method of a virtual energy storage power plant based on a non-precise alternative direction multiplier method, wherein the virtual energy storage power plant is an effective means for realizing large-scale distributed energy storage grid-connected operation in a power distribution network, but large-scale equipment in the virtual energy storage power plant is optimized and cooperated, and the requirement of the power grid on the rapidity of response is often difficult to meet. The invention provides a distributed control framework of a virtual energy storage power plant for large-scale energy storage, which enables an energy storage polymer to follow a target load curve and simultaneously can avoid the out-of-limit of the voltage of a power grid. The method is based on the non-precise alternative direction multiplier method, the closed analytic solution of the energy storage local optimization problem can be obtained in the original variable updating, the complexity of optimization calculation is obviously reduced, and the cooperation efficiency of the energy storage equipment is improved.

Description

Virtual energy storage power plant distributed cooperation method based on non-precise alternative direction multiplier method
Technical Field
The invention relates to the field of virtual energy storage, in particular to a distributed cooperation method of a virtual energy storage power plant based on a non-precise alternative direction multiplier method.
Background
Under the background of renewable energy utilization, energy storage becomes an important composition and support technology of a smart grid and a renewable energy high-occupancy energy system, various auxiliary services such as frequency modulation and voltage regulation are provided for the operation of a power system, and the consumption level of renewable energy such as wind and light is effectively improved. The traditional centralized energy storage is limited by factors such as geographical conditions, the installed capacity acceleration is in a descending trend year by year, and the demand-side distributed energy storage units are rapidly developed. In 2019, the national institutes of development and improvement and the national energy resource bureau jointly release opinions about deepening the construction of the electric power market, and clearly show that third parties such as energy storage facilities are encouraged to participate in auxiliary services of the electric power system.
The distributed energy storage has the characteristics of small single machine capacity and large quantity and scale, and the active and reactive adjustment of the distributed energy storage is more flexible. How to optimize and coordinate the large-quantity and different-characteristic generalized energy storage devices in the distribution network provides a flexible and reliable power source for the power grid is a main target of distributed energy storage control strategy research. The virtual energy storage power plant is an extension of the concept of the virtual power plant and is the aggregation of a series of multi-element distributed energy storage systems, so that the virtual energy storage power plant has the capacity and the power grid supporting capacity similar to those of a centralized energy storage power station, can be scheduled by a system operator like independent water pumping energy storage, solves the problem of optimal utilization of energy storage resources covering a large area, and becomes an effective means for large-scale grid-connected utilization of energy storage in a distribution network. Considering a large number of generalized distributed energy storage units in a future power grid, the traditional centralized control is difficult to bear the corresponding communication cost, and the completely distributed control cannot achieve the effect of cooperative control through the cooperation among the energy storage units. The distributed control adopting the adjacent communication principle has the advantages of strong anti-interference performance, good expansibility, plug and play, privacy protection and the like, is more suitable for the aggregation control of the distributed energy storage system in the current background, and also conforms to the general trend of decentralization of the power system.
The distributed optimization algorithm allocates the original centralized optimization problem to each intelligent agent for cooperative solution, can realize quick optimization scheduling of resources in the system, and has the difficulty of obtaining the distributed solution method of the optimization problem with better convergence. In a large number of existing researches, a distributed optimization method is used for frequency and voltage regulation of a power system, and special optimization problem design can effectively utilize the sparsity characteristic of the structure of the power system by adopting an Alternating Direction Method of Multiplies (ADMM) and a dual decomposition method. Although the ADMM algorithm embodies better convergence and robustness than a general first-order optimization method, considering that each iteration of the algorithm needs to solve a new optimization problem, the speed and the like of solving the optimal control problem by using the ADMM are still to be improved.
Disclosure of Invention
The invention provides a distributed cooperation method of a virtual energy storage power plant based on a non-precise alternative direction multiplier method, which mainly solves the problem of cooperative optimization operation of large-scale energy storage equipment in the virtual energy storage power plant, reduces the calculation complexity of each energy storage intelligent body by adopting a non-precise dual optimization method, improves the cooperation efficiency among the energy storage equipment, and meets the requirements of a power grid on the response speed of an energy storage polymer at the moment of load follow.
A virtual energy storage power plant distributed cooperation method based on a non-precise alternative direction multiplier method comprises the following steps:
step 1: determining a distributed energy storage access position in a power distribution network system, and acquiring topological and structural parameters of the power distribution network system so as to acquire a resistance matrix and a reactance matrix of a system power distribution line;
and 2, step: according to the topology and structure parameters of the power distribution network system obtained in the step 1, determining the voltage safety constraint of the power grid and the operation constraint of the energy storage equipment, and establishing an optimized scheduling model of the virtual energy storage power plant, wherein the optimized scheduling model aims at realizing target load curve following and reducing equipment use cost;
and step 3: aiming at each distributed energy storage agent in the virtual energy storage power plant, converting the optimized scheduling model of the virtual energy storage power plant established in the step (2) into a standard optimization model, and highlighting coupling and non-coupling constraint conditions in the standard optimization model by adopting a vectorization characterization method;
and 4, step 4: converting the standard optimization model obtained in the step 3 into a dual problem form, and determining a distributed solving method of the standard optimization model in the dual problem: and when the algorithm reaches a given maximum cycle number, the iteration result of the original variable is used as the optimal power setting of the distributed energy storage equipment to guide the charging and discharging management of the actual equipment.
Further, for the power distribution network system in the step 1, topology and structure parameters of the power distribution network system and resistance and reactance matrixes of system power distribution lines are obtained, and on the basis of the topology and structure parameters, a DistFlow power flow model of the system is obtained:
V=R l P n +X l Q n +V 0 (1)
R l =M -T D r M -1 ,X l =M -T D x M -1 (2)
in formula (1), V represents a column vector consisting of all node voltages, P n ,Q n Representing a column vector formed by the injected power of all nodes, i.e.
Figure BDA0003590845520000031
Figure BDA0003590845520000041
All nodes in the system are aggregated
Figure BDA00035908455200000415
Is shown by V 0 =1 N v 0 ,v 0 For voltage at the point of common connection of the distribution network, 1 N Represents a column vector of length N with all elements 1; the formula (2) gives the resistance and reactance matrix of the power distribution network
Figure BDA0003590845520000042
M is a system diagram description
Figure BDA0003590845520000043
Is derived from the system topology, D r 、D x Respectively, a diagonal matrix formed by all the line resistances and reactances in epsilon, M -T Representing the inverse of M and transpose.
Further, the optimized scheduling model in step 2 is represented as:
Figure BDA0003590845520000044
Figure BDA0003590845520000045
Figure BDA0003590845520000046
Figure BDA0003590845520000047
Figure BDA0003590845520000048
Figure BDA0003590845520000049
in formula (3):
Figure BDA00035908455200000410
the active power and the reactive power of the ith energy storage equipment at the moment t; the running cost of stored energy is expressed by quadratic form with the parameter of
Figure BDA00035908455200000411
Wherein
Figure BDA00035908455200000412
A quadratic term coefficient representing the active use cost of the ith energy storage device at the moment t,
Figure BDA00035908455200000413
a first term coefficient representing the active use cost of the ith energy storage device at the moment t,
Figure BDA00035908455200000414
primary and secondary term coefficients respectively representing the reactive power use cost of the ith energy storage device at the moment t;
in the DistFlow power flow model according to the system in equation (4),
Figure BDA0003590845520000051
Figure BDA0003590845520000052
respectively representing system resistance and reactance matrix R l ,X l The (c) th column of (a),
Figure BDA0003590845520000053
the load power of the respective node i,
Figure BDA0003590845520000054
indicating that there is no set of pure load nodes storing energy;
in the formula (5), the reaction mixture is,
Figure BDA0003590845520000055
the target load curve represents that the energy storage polymer has the capacity of participating in power grid dispatching along with the instruction;
equation (6) represents the capacity limit of the energy storage device,
Figure BDA0003590845520000056
is the reactive power of the energy storage device,
Figure BDA0003590845520000057
apparent power capacity of the energy storage converter;
equation (7) represents the active power constraint of the energy storage device, where
Figure BDA0003590845520000058
Minimum and maximum active capacity for energy storage;
equation (8) represents the charge capacity constraint of the stored energy, where SoC min ,SoC max For the minimum and maximum charge capacity limits of energy storage,
Figure BDA0003590845520000059
the charging capacity at the ith energy storage time t is shown, and eta is the charging and discharging efficiency of the energy storage.
Further, the optimal scheduling model in step 2 converts the quadratic constraint of equation (6) into a linear constraint to facilitate subsequent solution:
Figure BDA00035908455200000510
where the parameter k takes 8.
Further, the step 3 is that
Figure BDA00035908455200000511
Wherein x is i,t Is the decision variable, Δ x, at the ith moment of energy storage t i,t And representing a decision variable of each energy storage unit in the form of the sum of the power output at the last moment and the power offset at the next moment, wherein the standard optimization model is represented as:
Figure BDA0003590845520000061
Figure BDA0003590845520000062
Figure BDA0003590845520000063
equation (11) is a coupling constraint condition, equation (12) is a non-coupling constraint condition, and Δ r in equation (12) i,t For the introduced relaxation variables, the conversion from inequality constraint to equality constraint is realized;
parameter matrices and vectors A of equations (10) - (12) i 、E、C、b 0 、e 0 、d i Are obtained according to the optimization models (3) - (9), i.e.
Figure BDA0003590845520000064
Figure BDA0003590845520000065
Figure BDA0003590845520000066
Figure BDA0003590845520000067
Wherein
Figure BDA0003590845520000071
Further, in the step 4, firstly, a standard optimization model of the virtual energy storage power plant is processed, and manual constraint is introduced to realize distributed solution:
Figure BDA0003590845520000072
Figure BDA0003590845520000073
Figure BDA0003590845520000074
in the formula (18), in the above,
Figure BDA0003590845520000075
g i,t (Δx i,t ) And h i,t (Δx i,t ) Respectively, smooth and non-smooth portions, y, of the objective function (10) i,t ,z i,t Respectively of formula (11) and formula (12)Lagrange multipliers of constraints; t is t ij Y representing adjacent nodes as auxiliary variables i,t Equal; l i The auxiliary variables are introduced to ensure the convexity of the subsequent model.
Further, the standard optimization model obtained by the equations (18) to (20) in the step 4 can be in a distributed solving form, and the augmented lagrangian function is constructed as follows:
Figure BDA0003590845520000081
wherein
Figure BDA0003590845520000082
Lagrange multipliers of the constraints (19) and (20);
aiming at the augmented Lagrangian function of the formula (21), the following non-precise alternative direction multiplier method is utilized to realize the complete distributed solution of the optimization model:
Figure BDA0003590845520000083
Figure BDA0003590845520000084
Figure BDA0003590845520000085
Figure BDA0003590845520000091
where k denotes the number of iterations, [ ■ ]] + Max (0, ■), given parameters of the algorithm are σ, τ, c,
Figure BDA0003590845520000092
is the introduced algorithm state variable.
Furthermore, aiming at the original variable update described by the formula (21), according to the non-precise minimization thought, an analytic solution of the local optimization problem of each energy storage intelligent agent is deduced, and the formula (21) is converted into:
Figure BDA0003590845520000093
wherein
Figure BDA0003590845520000094
To represent
Figure BDA0003590845520000095
To x and is at
Figure BDA0003590845520000096
A gradient value of (a) and
Figure BDA0003590845520000097
the optimal solution for equation (26) is:
Figure BDA0003590845520000101
Figure BDA0003590845520000102
thereby, according to what is obtained
Figure BDA0003590845520000103
And obtaining the optimal power setting of each energy storage device at the current moment.
The invention can avoid the out-of-limit of the power grid voltage while the energy storage polymer follows the target load curve, and can obtain a closed analytic solution of the energy storage local optimization problem in the original variable update based on the non-precise alternative direction multiplier method, thereby obviously reducing the complexity of optimization calculation and improving the cooperative efficiency of energy storage equipment.
Drawings
FIG. 1 is a schematic flow diagram of the non-exact alternative direction multiplier method of the present invention;
FIG. 2 is a diagram of a distributed control architecture for a virtual energy storage plant in an IEEE33 node power distribution network in accordance with the present invention;
FIG. 3 is a diagram of the effect of virtual energy storage plant control based on the imprecise alternative direction multiplier method in an IEEE33 node system;
FIG. 4 is a diagram of the effect of virtual energy storage plant control based on the inaccurate alternative direction multiplier method in an IEEE69 node system.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The embodiment of the invention provides a virtual energy storage power plant distributed cooperation method based on a non-precise alternative direction multiplier method, which comprises the following steps:
step 1: the structure of the power distribution network required to implement energy storage distributed voltage support is shown in fig. 2, topology and structure parameters of the power distribution network system are obtained, and a resistance matrix and a reactance matrix of the power distribution network system are obtained according to the formulas (1) - (2).
Step 2: determining voltage safety constraints of an actual power grid and operation constraints of energy storage system equipment, and constructing an optimal scheduling model of a virtual energy storage power plant according to formulas (3) - (8);
and step 3: and (3) defining coupling and non-coupling constraint conditions of each energy storage agent in the virtual energy storage power plant, converting the optimized scheduling model constructed in the step (2) into standard optimized models in the formulas (10) - (12), and determining a parameter matrix according to the formulas (14) - (16).
And 4, step 4: using non-exact alternative direction multipliersThe method (as shown in fig. 1) performs distributed solution on the standard optimization model, selects the allowed voltage range of the power distribution network as 0.95-1.05 p.u., the voltage at the public connection point of the power distribution network as 1p.u., and commonly accesses 20 energy storage devices and 4 photovoltaic devices in the power distribution network, and for the distributed optimization algorithm adopted by the invention, selects σ as 0.01, τ as 0.05, β as 2e3, c as 1,
Figure BDA0003590845520000111
Figure BDA0003590845520000112
the power rating of the stored energy is randomly selected in the range of 0.5MW to 1 MW.
Finally, the voltage support effect of distributed energy storage in this example is shown in fig. 3, and fig. 3(a) shows the aggregate power of the virtual energy storage power plants, so that the virtual energy storage power plants can accurately follow the target load curve under the distributed and centralized control framework; fig. 3(b) shows the operation cost of the virtual energy storage plant under the centralized and distributed frameworks, and it can be seen that the operation cost of the virtual energy storage plant under the distributed control is slightly higher than that under the centralized control; fig. 3(c) and fig. 3(d) compare the voltage support effect of centralized control and distributed control, and it can be seen that in the two control modes, the virtual energy storage power plant successfully avoids the voltage out-of-limit while following the target load curve. Fig. 4 shows that the virtual energy storage plant in the IEEE69 node system achieves similar control effect as the IEEE33 node system, wherein the centralized control strategy can give the optimal power setting of the energy storage device within 1 minute, and the proposed distributed control can achieve similar energy storage synergistic effect within 10 seconds.
The above description is only an embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (8)

1. A distributed cooperation method of a virtual energy storage power plant based on a non-precise alternative direction multiplier method is characterized by comprising the following steps:
step 1: determining a distributed energy storage access position in a power distribution network system, and acquiring topological and structural parameters of the power distribution network system so as to acquire a resistance matrix and a reactance matrix of a system power distribution line;
step 2: according to the topology and structure parameters of the power distribution network system obtained in the step 1, determining the voltage safety constraint of a power grid and the operation constraint of energy storage equipment, and establishing an optimized scheduling model of a virtual energy storage power plant, wherein the optimized scheduling model aims at realizing target load curve following and reducing equipment use cost;
and 3, step 3: aiming at each distributed energy storage agent in the virtual energy storage power plant, converting the optimized scheduling model of the virtual energy storage power plant established in the step (2) into a standard optimization model, and highlighting coupling and non-coupling constraint conditions in the standard optimization model by adopting a vectorization characterization method;
and 4, step 4: converting the standard optimization model obtained in the step 3 into a dual problem form, and determining a distributed solving method of the standard optimization model in the dual problem: and when the algorithm reaches a given maximum cycle number, the iteration result of the original variable is used as the optimal power setting of the distributed energy storage equipment to guide the charging and discharging management of the actual equipment.
2. The virtual energy storage power plant distributed cooperation method based on the imprecise alternative direction multiplier method as defined in claim 1, wherein topology and structure parameters of the power distribution network system in the step 1 and resistance and reactance matrixes of system power distribution lines are obtained, and based on the topology and structure parameters, a DistFlow power flow model of the system is obtained:
V=R l P n +X l Q n +V 0 (1)
R l =M -T D r M -1 ,X l =M -T D x M -1 (2)
in formula (1), V represents a column vector consisting of all node voltages, P n ,Q n Representing a column vector formed by the injected power of all nodes, i.e.
Figure FDA0003590845510000021
Figure FDA0003590845510000022
All nodes in the system are aggregated
Figure FDA0003590845510000023
Is represented by V 0 =1 N v 0 ,v 0 For voltages at common connection points of distribution networks, 1 N Represents a column vector of length N with all elements 1; the resistance and reactance matrix R of the power distribution network is shown in the formula (2),
Figure FDA0003590845510000024
M is a system diagram description
Figure FDA0003590845510000027
Is obtained from the system topology, D r 、D x Is a diagonal matrix formed by all the line resistances and reactances in epsilon, M -T Representing the inverse of M and transpose.
3. The distributed coordination method for virtual energy storage power plant based on non-precise alternative direction multiplier method according to claim 2, characterized in that said optimal scheduling model in step 2 is expressed as:
Figure FDA0003590845510000025
Figure FDA0003590845510000026
Figure FDA0003590845510000031
Figure FDA0003590845510000032
Figure FDA0003590845510000033
Figure FDA0003590845510000034
in formula (3):
Figure FDA0003590845510000035
the active power and the reactive power of the ith energy storage equipment at the moment t; the running cost of the stored energy is represented by a quadratic form with the parameters of
Figure FDA0003590845510000036
Wherein
Figure FDA0003590845510000037
A quadratic term coefficient representing the active use cost of the ith energy storage device at the moment t,
Figure FDA0003590845510000038
a first term coefficient representing the active use cost of the ith energy storage device at the moment t,
Figure FDA0003590845510000039
respectively representing the first time and the second time of the reactive use cost of the ith energy storage equipment at the moment tA term coefficient;
in the DistFlow power flow model according to the system in equation (4),
Figure FDA00035908455100000310
Figure FDA00035908455100000311
respectively represent a system resistance matrix and a reactance matrix R l ,X l The (c) th column of (a),
Figure FDA00035908455100000312
the load power of the respective node i,
Figure FDA00035908455100000313
representing the absence of a set of energy-storing pure load nodes;
in the formula (5), the reaction mixture is,
Figure FDA00035908455100000314
the target load curve represents that the energy storage polymer has the capacity of participating in power grid dispatching along with the instruction;
equation (6) represents the capacity limit of the energy storage device,
Figure FDA00035908455100000315
is the reactive power of the energy storage device,
Figure FDA00035908455100000316
apparent power capacity of the energy storage converter;
equation (7) represents the active power constraint of the energy storage device, where
Figure FDA00035908455100000317
Minimum and maximum active capacity for energy storage;
equation (8) represents the charge capacity constraint of the stored energy, where SoC min ,SoC max For the minimum and maximum charge capacity limits of energy storage,
Figure FDA00035908455100000318
the charging capacity at the ith energy storage time t is shown, and eta is the charging and discharging efficiency of the energy storage.
4. The distributed cooperation method of the virtual energy storage power plant based on the imprecise alternating direction multiplier method as defined in claim 3, wherein the optimal scheduling model in step 2 converts quadratic constraint of equation (6) into linear constraint to facilitate subsequent solution:
Figure FDA0003590845510000041
where the parameter k takes 8.
5. The distributed coordination method for virtual energy storage power plant based on non-precise alternative direction multiplier method according to claim 4, characterized in that said step 3 order
Figure FDA0003590845510000042
Figure FDA0003590845510000043
Wherein x i,t Is the decision variable, Δ x, at the ith moment of energy storage t i,t And representing a decision variable of each energy storage unit in a form of the sum of the power output at the last moment and the power offset at the next moment, wherein the standard optimization model is represented as follows:
Figure FDA0003590845510000044
Figure FDA0003590845510000045
Figure FDA0003590845510000046
equation (11) is a coupling constraint condition, equation (12) is a non-coupling constraint condition, and Δ r in equation (12) i,t For the introduced relaxation variables, the conversion from inequality constraint to equality constraint is realized;
parameter matrices and vectors A of equations (10) - (12) i 、E、C、b 0 、e 0 、d i Are obtained according to the optimization models (3) - (9), i.e.
Figure FDA0003590845510000051
Figure FDA0003590845510000052
Figure FDA0003590845510000053
Figure FDA0003590845510000054
Wherein
Figure FDA0003590845510000055
6. The distributed collaborative method of the virtual energy storage plant based on the imprecise alternative direction multiplier method of claim 5, wherein in step 4, a standard optimization model of the virtual energy storage plant is first processed, and artificial constraints are introduced to realize distributed solution:
Figure FDA0003590845510000061
Figure FDA0003590845510000062
Figure FDA0003590845510000063
in the formula (18), the reaction mixture is,
Figure FDA0003590845510000064
g i,t (Δx i,t ) And h i,t (Δx i,t ) Respectively smooth and non-smooth portions, y, in the objective function (10) i,t ,Z i,t Lagrange multipliers which are constraints of equations (11) and (12), respectively; t is t ij Y representing adjacent nodes as auxiliary variables i,t Equal; l. the i The auxiliary variables are introduced to ensure the convexity of the subsequent model.
7. The distributed cooperative method for the virtual energy storage power plant based on the imprecise alternative direction multiplier method as defined in claim 6, wherein the standard optimization model obtained in step 4 through equations (18) - (20) can be in a distributed solution form, and an augmented Lagrangian function of the standard optimization model is constructed as follows:
Figure FDA0003590845510000066
wherein
Figure FDA0003590845510000065
Lagrange multipliers with constraint conditions of (19) and (20);
aiming at the augmented Lagrangian function of the formula (21), the following non-precise alternative direction multiplier method is utilized to realize the complete distributed solution of the optimization model:
Figure FDA0003590845510000071
Figure FDA0003590845510000072
Figure FDA0003590845510000073
Figure FDA0003590845510000074
where k denotes the number of iterations, [ ■ ]] + Max (0, ■), given parameters of the algorithm are σ, τ, c,
Figure FDA0003590845510000075
is the introduced algorithm state variable.
8. The distributed coordination method for virtual energy storage power plant based on non-precise alternative direction multiplier method according to claim 7, characterized in that, for the original variable update described by equation (21), according to the non-precise minimization thought, the analytic solution of local optimization problem of each energy storage agent is derived, and equation (21) is converted into:
Figure FDA0003590845510000081
wherein
Figure FDA0003590845510000082
Represent
Figure FDA0003590845510000083
To x and in
Figure FDA0003590845510000084
A gradient value of (a) and
Figure FDA0003590845510000085
the optimal solution for equation (26) is:
Figure FDA0003590845510000086
Figure FDA0003590845510000087
thereby, according to what is obtained
Figure FDA0003590845510000091
And obtaining the optimal power setting of each energy storage device at the current moment.
CN202210377568.XA 2022-04-12 Virtual energy storage power plant distributed cooperation method based on inaccurate alternating direction multiplier method Active CN114925880B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210377568.XA CN114925880B (en) 2022-04-12 Virtual energy storage power plant distributed cooperation method based on inaccurate alternating direction multiplier method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210377568.XA CN114925880B (en) 2022-04-12 Virtual energy storage power plant distributed cooperation method based on inaccurate alternating direction multiplier method

Publications (2)

Publication Number Publication Date
CN114925880A true CN114925880A (en) 2022-08-19
CN114925880B CN114925880B (en) 2024-06-21

Family

ID=

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117236991A (en) * 2023-11-15 2023-12-15 广东工业大学 Distributed resource aggregation modeling method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110266038A (en) * 2019-05-28 2019-09-20 广东电网有限责任公司电力调度控制中心 A kind of more virtual plant distributed coordination regulation methods
CN111080022A (en) * 2019-12-23 2020-04-28 国网四川省电力公司经济技术研究院 Partition distributed coordination optimization method containing multiple benefit agents
US20210184494A1 (en) * 2019-10-29 2021-06-17 Tsinghua University Methods and apparatuses for identifying power feasible region of virtual power plant
CN112990596A (en) * 2021-03-31 2021-06-18 东南大学 Distributed optimization method for cooperative operation of active power distribution network and virtual power plant
CN113793029A (en) * 2021-09-14 2021-12-14 国网上海市电力公司 Virtual power plant optimal scheduling method and device
CN113890021A (en) * 2021-09-29 2022-01-04 国网综合能源服务集团有限公司 Multi-virtual power plant distributed transaction method considering constraint of power distribution network
CN114140022A (en) * 2021-12-10 2022-03-04 国网山西省电力公司电力科学研究院 Multi-virtual power plant distributed dynamic economic dispatching method and system

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110266038A (en) * 2019-05-28 2019-09-20 广东电网有限责任公司电力调度控制中心 A kind of more virtual plant distributed coordination regulation methods
US20210184494A1 (en) * 2019-10-29 2021-06-17 Tsinghua University Methods and apparatuses for identifying power feasible region of virtual power plant
CN111080022A (en) * 2019-12-23 2020-04-28 国网四川省电力公司经济技术研究院 Partition distributed coordination optimization method containing multiple benefit agents
CN111768054A (en) * 2019-12-23 2020-10-13 国网四川省电力公司经济技术研究院 Partition distributed coordination optimization method containing multiple benefit agents
CN112990596A (en) * 2021-03-31 2021-06-18 东南大学 Distributed optimization method for cooperative operation of active power distribution network and virtual power plant
CN113793029A (en) * 2021-09-14 2021-12-14 国网上海市电力公司 Virtual power plant optimal scheduling method and device
CN113890021A (en) * 2021-09-29 2022-01-04 国网综合能源服务集团有限公司 Multi-virtual power plant distributed transaction method considering constraint of power distribution network
CN114140022A (en) * 2021-12-10 2022-03-04 国网山西省电力公司电力科学研究院 Multi-virtual power plant distributed dynamic economic dispatching method and system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
赵晨;何宇俊;罗钢;龚超;赵越;张轩;陈启鑫;: "基于泛在互联的虚拟电厂参与实时市场模型", 电力建设, no. 06, 1 June 2020 (2020-06-01) *
陈厚合;王子璇;张儒峰;姜涛;李雪;李国庆;: "含虚拟电厂的风电并网系统分布式优化调度建模", 中国电机工程学报, no. 09, 5 May 2019 (2019-05-05) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117236991A (en) * 2023-11-15 2023-12-15 广东工业大学 Distributed resource aggregation modeling method

Similar Documents

Publication Publication Date Title
WO2021135332A1 (en) Edge computing-based multi-agent load regulation and control method
CN108363306B (en) Micro-grid distributed controller parameter determination method based on linear quadratic optimization
Mohsenian-Rad et al. Towards building an optimal demand response framework for DC distribution networks
CN112163700B (en) Electrochemical energy storage power station planning method considering cycle life of energy storage battery
CN105811407B (en) A kind of micro-capacitance sensor primary frequency modulation control method based on distributed Newton method
CN110265991B (en) Distributed coordination control method for direct-current micro-grid
CN107546743A (en) A kind of radial distribution networks distributed optimization trend method
Wang et al. Load curve smoothing strategy based on unified state model of different demand side resources
CN110826880A (en) Active power distribution network optimal scheduling method for large-scale electric vehicle access
CN109617147B (en) Power electronic transformer operation strategy optimization combination method
Rahman et al. Continuous-domain real-time distributed ADMM algorithm for aggregator scheduling and voltage stability in distribution network
CN107528321B (en) Random optimal distribution method for power of interconnected power grid alternating current-direct current parallel power transmission channel
CN115333110A (en) Power distribution network-microgrid group collaborative distributed optimization scheduling method and system based on ADMM
CN115545291A (en) Centralized sharing energy storage capacity configuration and pricing method based on master-slave game
CN116090739A (en) Virtual power plant optimal scheduling method comprehensively considering ordered charging of electric vehicles
CN104659812A (en) Multi-microgrid coordination control method based on predictive control
CN105426985B (en) Urban distribution network voltage distribution optimization method suitable for roof photovoltaic access
Dutta et al. Contribution of electric vehicles for frequency regulation in presence of diverse power sources and transmission links
CN107508318B (en) Active control method and system based on voltage sensitivity partition
CN114925880A (en) Virtual energy storage power plant distributed cooperation method based on non-precise alternative direction multiplier method
CN114709848A (en) Multiphase power distribution network energy storage optimization scheduling method considering power quality compensation
CN114925880B (en) Virtual energy storage power plant distributed cooperation method based on inaccurate alternating direction multiplier method
Xu et al. Asynchronous distributed dynamic economic dispatch in active distribution networks
CN109345030B (en) Multi-microgrid comprehensive energy system thermoelectric energy flow distribution type optimization method and device
CN111932057A (en) AC/DC micro-grid economic dispatching method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant