CN117410990B - Distributed energy distributed control method and system based on local calculation - Google Patents
Distributed energy distributed control method and system based on local calculation Download PDFInfo
- Publication number
- CN117410990B CN117410990B CN202311718732.XA CN202311718732A CN117410990B CN 117410990 B CN117410990 B CN 117410990B CN 202311718732 A CN202311718732 A CN 202311718732A CN 117410990 B CN117410990 B CN 117410990B
- Authority
- CN
- China
- Prior art keywords
- representing
- distributed
- power
- control
- node
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 238000011217 control strategy Methods 0.000 claims abstract description 20
- 230000005855 radiation Effects 0.000 claims abstract description 10
- 239000000243 solution Substances 0.000 claims description 43
- 230000006870 function Effects 0.000 claims description 30
- 238000004146 energy storage Methods 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 9
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000010248 power generation Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000001364 causal effect Effects 0.000 claims description 3
- 230000002441 reversible effect Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims 1
- 238000004891 communication Methods 0.000 abstract description 5
- 230000007547 defect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 101100272279 Beauveria bassiana Beas gene Proteins 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Public Health (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to the technical field of power distribution network optimal scheduling, and discloses a distributed energy distributed control method and system for a power distribution network based on local calculation, wherein the distributed energy distributed control method comprises the following steps: acquiring resource information of a power grid, and constructing a distributed control objective function and constraint conditions of the power distribution network; constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy; constructing a polyhedron internal approximation of a feasible domain of the model by a radiation function simplified control method, and solving an optimal solution upper bound of the decentralized control model; constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a distributed control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound. Reactive voltage control is achieved without communication. The method can effectively solve the defects of high-efficiency, reliable and safe distribution of electric energy and overhigh communication cost, can improve the robustness and reliability of a power grid, and can reduce the expected running cost to the maximum extent.
Description
Technical Field
The invention relates to the technical field of power distribution network optimal scheduling, in particular to a distributed energy distributed control method and system based on local calculation.
Background
The distribution network refers to a network system that distributes electric energy transmitted from a transmission network to end users. Most of the traditional power distribution network systems are centralized control modes, namely, a central control center monitors and controls the whole power grid. With rapid development and large-scale access of renewable energy sources such as wind energy and solar energy, the traditional centralized control mode cannot meet the requirements of an electric power system under the conditions of new energy source access and load side participation. Meanwhile, the complexity of the power system is continuously improved along with the increase of the accessed distributed power supply, and the traditional centralized control mode cannot meet the requirements of the control system on response speed, control precision and reliability. For these problems, distributed control techniques have evolved. The distributed control technology distributes control tasks to each distributed power supply and load, so that local coordination and global coordination are realized. The method can improve the energy utilization rate, reliability and stability of the power grid, reduce the dependence of the power grid on central control, and improve the robustness and the anti-interference capability of the power grid.
Disclosure of Invention
The present invention has been made in view of the above-described problems.
Therefore, the technical problems solved by the invention are as follows: the method can effectively solve the defects that the traditional centralized control mode in the running and dispatching process of the power distribution network is difficult to meet the requirements on efficient, reliable and safe distribution of electric energy and the communication cost is too high, improve the robustness and reliability of the power distribution network, and simultaneously reduce the expected running cost to the maximum extent.
In order to solve the technical problems, the invention provides the following technical scheme: a distributed energy distributed control method based on local calculation comprises the following steps:
acquiring resource information of a power grid, and constructing a distributed control objective function and constraint conditions of the power distribution network;
constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy;
constructing a polyhedron internal approximation of a feasible domain of the model by a radiation function simplified control method, and solving an optimal solution upper bound of the decentralized control model;
constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a distributed control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound.
As a preferable scheme of the distributed energy source decentralized control method based on the local computing, the distributed energy source decentralized control method based on the local computing comprises the following steps: the resource information comprises line information, power information, energy storage information, current and voltage information of the power distribution network.
As a preferable scheme of the distributed energy source decentralized control method based on the local computing, the distributed energy source decentralized control method based on the local computing comprises the following steps: the distributed control objective function of the power distribution network comprises that the objective function of distributed control of the power distribution network is the sum of the minimum active network loss and the energy storage state;
;
wherein,representing the sum of the energy storage states +.>Representing active power loss of distribution networkConsumption (S)/(L)>Representing the state of energy storage->Is a coefficient vector +_>Indicating line->Active loss of->Represents the line resistance +.>Representing the active power of the line,/->Representing reactive power of the line, +.>Representing the voltage amplitude at bus 0, +.>Representing branchesConnected head-end node->Representing the total time period, +.>Indicating the total number of stored energy.
As a preferable scheme of the distributed energy source decentralized control method based on the local computing, the distributed energy source decentralized control method based on the local computing comprises the following steps: the constraint conditions comprise the following steps of network constraint of the power distribution network:
;
wherein,indicating line->Active power of transmission of +.>Indicating line->Is used for the transmission of reactive power,representation section->Voltage amplitude at>Representing node->Voltage amplitude of>Indicating line->Is (are) electric conduction>Indicating line->Susceptance of->Indicating line->Resistance of->Indicating line->Reactance of->Indicating line->Square of current amplitude>Indicating line->Middle node->And node->Phase angle difference between;
energy storage constraint:
;
;
wherein,is->State of charge of the energy storage device at the moment +.>Is the self-discharge coefficient of the energy storage device,representation->Charging power of the time-of-day energy storage device, < >>Representation->Discharge power of the time-of-day energy storage device, +.>Indicating the charging efficiency of the stored energy>Indicating the discharge efficiency of the stored energy>Indicating the rated capacity of the stored energy>Is indicated at +.>Energy storage->Is used for storing electric quantity;
;
wherein,representing the minimum value of the state of charge of the stored energy, < + >>Representing a maximum value of the stored state of charge;
distributed photovoltaic power generation output constraint:
;
wherein,representing the capacity of node i distributed photovoltaic, +.>Representing the photovoltaic utilization of node i for period T,representing the theoretical output of node i distributed photovoltaic at time t,/->Representing photovoltaic utilization limit,/->Representing the active output of the node i distributed photovoltaic; />Representing reactive power of the node i distributed photovoltaic;
load constraint;
;
wherein,representing node->Time->Active power demand of->Representing node->Time->Reactive power demand, < >>Representing node->、/>Active power injection->Representing node->、/>Reactive power injection is performed;
voltage amplitude constraint;
;
wherein,representing the upper limit of the node voltage amplitude,/->Representing a node voltage magnitude lower limit;
uncertainty model:
modeling active power injection and active and reactive power requirements of a photovoltaic inverter as discrete time random processes, perturbing the process with each busIn association, define perturbations as:
;
wherein,representing the disturbance function at busbar iGather (S)>Representing load active disturbance, < >>Representing load reactive disturbance and photovoltaic active disturbance; />Representing a disturbance value;
the full disturbance trajectory is defined as:
;
wherein R represents a natural number,representing the total number of disturbances.
As a preferable scheme of the distributed energy source decentralized control method based on the local computing, the distributed energy source decentralized control method based on the local computing comprises the following steps: constructing a discrete time state space model describing collective dynamics of the power distribution network on the basis of the independent resource model; from the following componentsSubsystem composition, each subsystem->Indicating connection to busbar->Storing state for each subsystem>Input->:
;
The state equation of the whole system represents:
;
wherein the matrixExpressed as follows>Represents the Kronecker product;
;
definition at timeEach subsystem i determines the local control input based on the available information, proposes a fully decentralized disturbance feedback control strategy, which is +_ in time for each subsystem i>Is limited to:
;
wherein,a causal measurable function representing the local disturbance history, the local control strategy defining subsystem i is defined as:
;
wherein,for decentralized control strategy->Representing a family of all allowable decentralized control resources;
defining an optimal decentralized control problem:
;
wherein,representing a set of feasible states; />Representing a set of possible control inputs; />Representing a model objective function representing a minimization of expected active power loss versus storage state sum expectations, the expression:
;
wherein,representing the desire; />The method comprises the steps of inputting an equivalent coefficient matrix; />In order to perturb the matrix of equivalent coefficients,representing the disturbance.
As a preferable scheme of the distributed energy source decentralized control method based on the local computing, the distributed energy source decentralized control method based on the local computing comprises the following steps: the radiation function comprises the steps of simplifying a control problem through the radiation function, and constructing a polyhedron internal approximation of a feasible domain of the optimal decentralized control problem of the model; replacing the convex quadratic constraint with the linear constraint to derive a convex setIs a polyhedral approximation of:
;
deterministic constantExpressed as:
;
according toApproximating the optimal decentralized control problem as:
;
wherein the coefficient matrixDesignated by the underlying problem data, m represents the number of this inequality constraint,N u the total number of inputs is indicated,N x representing the total number of states +.>Respectively represent state variables +>Input variable->And random variable->Is a lower coefficient matrix of (a);
a subspace of a feasible decentralized affine control strategy is defined based on the following:
;
wherein S represents a constraint set of distributed control information,requiring a linear operator +.>Satisfying the information structure embedded in the feasible decentralized control strategy set Γ;
;
wherein,representation->Unit vector of>Representing optimal solution->As an auxiliary variable, a control signal is provided,representing a positive cone->Representing a semi-positive fixation cone->Representing the input of the system, M represents the second moment matrix of the system disturbanceNot only reversible but also positive, and->Indicating the desire.
As a preferable scheme of the distributed energy source decentralized control method based on the local computing, the distributed energy source decentralized control method based on the local computing comprises the following steps: the optimal solution further comprises the steps of constructing an outer polyhedron approximation of a feasible domain of a completely dispersed disturbance feedback control strategy, and solving an optimal value lower bound of the outer approximation;
for nodesThe convex quadratic constraint is equivalent to a linear inequality:
;
according toApproximating optimal decentralized control as:
;
wherein,a representative trace; />Representing an optimal solution, which is an optimal solution lower bound of optimal decentralized control;respectively represent state variables +>Input variable->And random variable->Is a coefficient matrix of the upper coefficient;
solving:
;
will beAs a choice of an optimal solution.
The system adopting the distributed energy distributed control method based on the local calculation is characterized in that: the model construction module is used for acquiring the resource information of the power grid and constructing a distributed control objective function and constraint conditions of the power distribution network; constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy; the analysis module is used for constructing a polyhedron internal approximation of a feasible domain of the model through a radiation function simplified control method and solving an optimal solution upper bound of the decentralized control model; constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a distributed control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound.
A computer device, comprising: a memory and a processor; the memory stores a computer program characterized in that: the processor, when executing the computer program, implements the steps of the method of any of the present invention.
A computer-readable storage medium having stored thereon a computer program, characterized by: the computer program, when executed by a processor, implements the steps of the method of any of the present invention.
The invention has the beneficial effects that: the distributed energy distributed control method based on local calculation provided by the invention aims at the serious challenges brought to a power distribution system by increasingly increasing distributed resources: voltage amplitude at the bus rapidly fluctuates, substation tidal current foldback, and power quality problems due to the intermittence of renewable energy supplies. The control tasks are distributed to the distributed power supplies and the loads through a distributed control technology, and reactive voltage control is achieved without communication. The method can effectively solve the defects that the traditional centralized control mode in the running and dispatching process of the power distribution network is difficult to meet the requirements on efficient, reliable and safe distribution of electric energy and the communication cost is too high, improves the robustness and reliability of the power distribution network, and simultaneously reduces the expected running cost to the maximum extent.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for the description of the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is an overall flowchart of a distributed energy distributed control method for a power distribution network based on local computing according to a first embodiment of the present invention;
fig. 2 is a local voltage decentralized control diagram in a distributed energy decentralized control method of a power distribution network based on local computation according to a first embodiment of the present invention;
fig. 3 is a three-node actual engineering topology diagram of a distributed energy distributed control method of a power distribution network based on local computation according to a second embodiment of the present invention;
fig. 4 is a graph of reactive voltage decentralized control effect of scenario 1 according to a distributed energy distributed control method based on local computation according to a second embodiment of the present invention, where a represents voltage control effect and b represents reactive adjustment condition;
fig. 5 is a graph of reactive voltage decentralized control effect of scenario 2 according to a distributed energy source decentralized control method based on local computing according to a second embodiment of the present invention, where a represents voltage control effect and b represents reactive adjustment condition.
Detailed Description
So that the manner in which the above recited objects, features and advantages of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to the embodiments, some of which are illustrated in the appended drawings. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
Example 1
Referring to fig. 1-2, for one embodiment of the present invention, a distributed energy source decentralized control method for a power distribution network based on local computing is provided, including:
s1: and acquiring the resource information of the power grid, and constructing a distributed control objective function and constraint conditions of the power distribution network.
Further resource information includes line information, power information, energy storage information, and current voltage information of the power distribution network.
The decentralized control objective function of the power distribution network is as follows:
;
the objective function of distributed control of the power distribution network is to minimize the sum of active network losses and energy storage states. Wherein the method comprises the steps ofAndrespectively representing the sum of the energy storage states and the active power loss of the power distribution network. Constraint (2)/(>Representing the state of energy storage->Is a coefficient vector. Constraint (3)/(x)>Indicating line->The specific functional expression is shown in formula (4). Constraint (4)/(A)>Represents the line resistance +.>Representing the active and reactive power of the line, +.>Representing the voltage amplitude at bus 0.
The distribution network decentralized control constraint conditions are as follows:
a. network constraints of the distribution network:
;
;
equation (5) and equation (6) are the result of linearization of the power flow constraint, whereIndicating line->Is used for the transmission of active power,indicating line->Reactive power of transmission, +.>Respectively represent node->The magnitude of the voltage at that point,indicating line->Is/are>Representing nodes/>Phase angle deviation of (2); formula (7) represents a square constraint of voltage, +.>Indicating line->Resistance and reactance of (a); formula (8) represents a line current constraint, +.>Indicating line->Is the square of the current amplitude of (a).
b. Energy storage constraint
;
In the method, in the process of the invention,is->The state of charge of the energy storage device at the moment; />Is->The state of charge of the energy storage device at time +1; />Is the self-discharge coefficient of the energy storage device; />And->Respectively indicate->Charging and discharging power of the energy storage device at any time; />And->Respectively representing the charge and discharge efficiency of the stored energy; />Representing the rated capacity of the stored energy;is indicated at +.>Energy storage->Is used for storing electric quantity; />Indicating that the regulation step length is 1h.
In the method, in the process of the invention,and->Representing the minimum and maximum values of the stored state of charge.
c. Distributed photovoltaic power generation output constraint
;
In the method, in the process of the invention,the capacity of the distributed photovoltaic for node i; />Reactive power output at the moment of node i distributed photovoltaic t is provided; />Active output at the moment of node i distributed photovoltaic t; />Photovoltaic utilization for node i of period T, < >>The theoretical output force (the non-limit front output force) of the distributed photovoltaic of the node i at the moment t; />Is the minimum limit value of the photovoltaic utilization rate.
d. Load constraint
;
;/>
In the method, in the process of the invention,respectively represent node->Time->Active power and reactive power requirements of (a).Respectively represent node->Active reactive power injection.
e. Voltage amplitude constraint
;
In the method, in the process of the invention,respectively representing upper and lower limits of node voltage amplitude. According to the specification of GB/T40427-2021, a node voltage constraint range is set to be 0.93UN-1.07 UN, wherein UN is rated voltage (per unit value is 1).
f. Uncertainty model
Active power injection and active and reactive power requirements of a photovoltaic inverter are modeled as discrete time random processes. The disturbance process is associated with each busIn association, define perturbations as:
;
the full disturbance trajectory is defined as:
;
wherein,representing the number of disturbance tracks as 1+3nT。
S2: and constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy.
On the basis of the established independent resource model, a discrete time state space model describing the collective dynamics of the power distribution network is constructed. The system consists ofSubsystem composition, wherein each subsystem +.>Representing a linkIs connected to bus bar->Is a collection of resources. Storing local inputs for subsystems>The method comprises the steps of energy storage power generation condition and photovoltaic active power output:
;
the state equation of the whole system is expressed as follows:
;
wherein the matrixRepresenting disturbance factor, ++>Represents the Kronecker product, +.>Representing an n-order identity matrix:
;
definition at timeEach subsystem->Determining local control inputs based on available information, a fully decentralized disturbance feedback control strategy is proposed, i.e. for each subsystem +.>It is at time->Is limited to: />
;
In the method, in the process of the invention,causal measurable function representing the history of local disturbances, defining subsystem +.>Is defined as +.>. Wherein->For decentralized control strategy->Is a family of all allowable decentralized control resources. Fig. 2 is a local voltage decentralized control.
The optimal decentralized control problem is defined as follows:
;
wherein,representing a set of feasible states; />Representing a set of possible control inputs; />Representing a model objective function representing a desire to minimize the sum of the expected active power loss and the stored state, the specific expression is as follows:
;
wherein,representing the desire; />The method comprises the steps of inputting an equivalent coefficient matrix; />Is a disturbance equivalent coefficient matrix.
S3: by means of a radiation function simplified control method, a polyhedral internal approximation of a feasible domain of the model is built, and an optimal solution upper bound of the decentralized control model is solved.
Considering that the model (23) is an infinite-dimensional convex problem, the calculation solution is difficult to carry out. The control problem is simplified by the radial function, and a polyhedral internal approximation of the feasible region of the model (23) is constructed.
The convex quadratic constraint (12) is replaced by the following linear constraint to derive a convex setIs a polyhedral approximation of:
;
thus, the deterministic constantCan be expressed as:
;
according toThe model (23) can be approximated as: />
;
In the coefficient matrixMay be specified in terms of underlying issue data.
A subspace of a feasible decentralized affine control strategy is defined based on the following:
;
where S represents a constraint set of distributed control information,requiring a linear operator +.>The information structure embedded in the feasible set of decentralized control policies Γ is satisfied.
;
In the method, in the process of the invention,representation->The optimal value of the problem (29) is equal to the cost generated by the optimal affine control strategy of the problem (27). Adopts->Representing the optimal solution of problem (29), it is apparent that the optimal solution upper bound of problem (23), i.e.>Representing a positive cone->Representing a semi-positive fixation cone->Representing the input of the system, M represents the second moment matrix of the system disturbance, +.>Not only reversible but also positive, and->Indicating the desire.
S4: constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a distributed control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound.
Since the affine strategy calculated in S3 is suboptimal, the strategy induced suboptimal is constrained in this step by solving the lower bound of the optimal value of the problem (23). An outer polyhedral approximation of the feasible region of the problem (21) is constructed and the lower bound of the optimal value of the outer approximation is solved.
For nodesThe convex quadratic constraint (12) can be equivalent to the following linear inequality:
;
according toThe model (23) can be approximated as (31):
;
wherein,a representative trace; />Representing a suboptimal solution, which is an optimal solution lower bound for optimal decentralized control;respectively represent state variables +>Input variable->And random variable->Is a matrix of coefficients.
By usingRepresenting the optimal solution of problem (31), which is the lower boundary of the optimal solution of problem (23), the following solution can be obtained>Wherein->The solutions of (29) and (31), respectively, demonstrate that the affine strategy is a near-optimal solution choice.
A distributed energy distributed control system of a power distribution network based on local calculation is characterized in that: the model construction module is used for acquiring the resource information of the power grid and constructing a distributed control objective function and constraint conditions of the power distribution network; constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy; the analysis module is used for constructing a polyhedron internal approximation of a feasible domain of the model through a radiation function simplified control method and solving an optimal solution upper bound of the decentralized control model; constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a distributed control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, database, or other medium used in the various embodiments provided herein may include at least one of non-volatile and volatile memory. The nonvolatile memory may include read only memory, magnetic tape, floppy disk, flash memory, optical memory, high density embedded nonvolatile memory, resistive memory, magnetic memory, ferroelectric memory, phase change memory, graphene memory, and the like. Volatile memory can include random access memory, external cache memory, or the like. By way of illustration, and not limitation, RAM can take many forms, such as static random access memory or dynamic random access memory. The databases referred to in the various embodiments provided herein may include at least one of relational databases and non-relational databases. The non-relational database may include, but is not limited to, a blockchain-based distributed database, and the like.
The processors referred to in the embodiments provided herein may be general purpose processors, central processing units, graphics processors, digital signal processors, programmable logic units, quantum computing-based data processing logic units, etc., without being limited thereto.
Example 2
Referring to fig. 3-5, for one embodiment of the present invention, a distributed energy distributed control method for a power distribution network based on local computing is provided, and in order to verify the beneficial effects of the present invention, scientific demonstration is performed through economic benefit computing and simulation experiments.
Firstly, the section selects a specific power system comprising 3 nodes for testing, and the test result and the actual condition are verified through an algorithm. The voltage control is carried out by using the local voltage and reactive power information of the nodes, the topological structure of the system is shown in the following figure 3, 1, 2 and 3 represent buses, wherein the node 1 is a balance node, the node 2 is provided with energy storage equipment, the node 3 is provided with distributed photovoltaic, a certain reactive power regulation capability can be provided for the system, and the voltage deviation is limited within a range of +/-5 percent, so that simulation verification is carried out.
The voltage of the node 3 is higher in the case of moderate load during the photovoltaic large-power period, and as shown in fig. 4a, the reactive power is correspondingly adjusted according to the voltage so as to control the voltage within the limit value. In case of a higher voltage, the voltage is gradually reduced by reducing the reactive output as shown in fig. 4b, and after 4 iterations, it is substantially stabilized at the rated voltage 1p.u, at which time the reactive output is no longer reduced, keeping the voltage stable. Wherein fig. 4a is the corresponding voltage amplitude after reactive power adjustment and fig. 4b is the amplitude of reactive power adjustment.
Under the condition that no photovoltaic output is generated but the load is in large demand, the voltage of the node 3 is lower, as shown in fig. 5a, under the condition that the voltage is lower, the reactive output simulation effect is increased, as shown in fig. 5b, under the condition that the reactive output is increased, the voltage gradually rises, and after 4 iterations, the voltage is basically stabilized at the rated voltage 1p. Wherein fig. 5a is the corresponding voltage amplitude after reactive power adjustment, and fig. 5b is the amplitude of reactive power adjustment.
It should be noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that the technical solution of the present invention may be modified or substituted without departing from the spirit and scope of the technical solution of the present invention, which is intended to be covered in the scope of the claims of the present invention.
Claims (9)
1. The distributed energy distributed control method based on the local calculation is characterized by comprising the following steps of:
acquiring resource information of a power grid, and constructing a distributed control objective function and constraint conditions of the power distribution network;
constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy;
constructing a polyhedron internal approximation of a feasible domain of the model by a radiation function simplified control method, and solving an optimal solution upper bound of the decentralized control model;
constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a decentralized control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound;
wherein the uncertainty model in the constraint condition comprises:
active power injection and active and reactive power requirements of the photovoltaic inverter are modeled as discrete time random processes, a perturbation process is associated with each busbar i, and the perturbation is defined as:
wherein, xi i (t) represents the set of perturbation functions at bus i,representing load active disturbance, < >>Representing load reactive disturbance, +.>Representing an active disturbance of the photovoltaic; zeta type toy i (. Cndot.) represents a perturbation value;
on the basis of the independent resource model, the construction of the integral dynamic time state space model comprises the following steps:
consists of n subsystems, each subsystem i=1,..n represents a collection of resources connected to bus i, for each subsystem state x is stored i (t) input u i (t):
The state equation of the whole system represents:
x(t+1)=x(t)+Bu(t)
wherein, the matrix B is expressed as follows,represents the Kronecker product;
defining at time t, each subsystem i determines a local control input based on available information, proposing a fully decentralized disturbance feedback control strategy, for each subsystem i its control input at time t being limited to:
wherein, gamma i (. T) represents a causal measurable function of the local disturbance history, the local control strategy defining subsystem i is defined as:
γ i =(γ i (·,0),...,γ i (·,T-1))
wherein, gamma i For a decentralized control strategy Γ represents a family of all allowable decentralized control resources;
defining an optimal decentralized control problem:
minimize J(γ)
subject toγ∈Γ
x(t+1)=x(t)+Bu
u=γ(ξ)
wherein,representing a set of feasible states; />Representing a set of possible control inputs; j (γ) represents a model objective function representing a minimization of expected active power loss versus storage state summation expectations, expression:
J(γ)=E γ [c′x+(L u u+L ξ ξ)′Σ(L u u+L ξ ξ)]
wherein E is γ [·]Representing the desire; l (L) u The method comprises the steps of inputting an equivalent coefficient matrix; l (L) ξ For the disturbance equivalent coefficient matrix, ζ represents the disturbance.
2. The distributed energy distribution network distributed energy source decentralized control method based on local computing according to claim 1, wherein: the resource information comprises line information, power information, energy storage information, current and voltage information of the power distribution network.
3. The distributed energy distributed control method based on local computing as claimed in claim 2, wherein: the distributed control objective function of the power distribution network comprises that the objective function of distributed control of the power distribution network is the sum of the minimum active network loss and the energy storage state;
C toal =C E +C loss
δp ij (t)=r ij (p ij (t) 2 +q ij (t) 2 )/v 0 (t) 2
wherein C is E Representing the sum of the stored energy states, C loss Representing active power loss of a power distribution network, x represents an energy storage state, c' is a coefficient vector, and δp ij (t) represents the active loss of the line ij, r ij Representing the line resistance, p ij Representing the active power of the line, q ij Representing reactive power of the line, v 0 Representing the voltage amplitude at bus 0, i, j representing the first and last nodes connected by branch ij, T representing the total time period, n representing the total number of stored energy。
4. A distributed energy distribution network distributed energy source decentralized control method based on local computing according to claim 3, wherein: the constraints further include, distribution network constraints:
P ij =g ij (U i -U j )-b ij θ ij
Q ij =-b ij (U i -U j )-g ij θ ij
wherein P is ij Representing the transmitted active power, Q, of line ij ij Representing the reactive power, U, of the transmission of the line ij i Representing the voltage amplitude at node i, U j Representing the voltage amplitude, g, at node j ij Representing the conductance of line ij, b ij The susceptance, r, representing the line ij ij Representing the resistance, x, of the line ij ij Representing reactance of line ij, l ij Represents the square, θ of the current amplitude of line ij ij Representing the phase angle difference between node i and node j in line ij;
energy storage constraint:
wherein SOC (t) is the state of charge of the energy storage device at time t, sigma is the self-discharge coefficient of the energy storage device, and P ch (t) represents the charging power of the energy storage device at time t, P dis (t) represents the discharge power, eta of the energy storage device at the time t ch Represents the charging efficiency, eta of energy storage dis Indicating the discharge efficiency of the stored energy, E ess Indicating the rated capacity of the stored energy,representing the stored electricity quantity of the energy storage i at time t;
SOC min ≤SOC(t)≤SOC max
wherein SOC is min Representing the minimum value of the state of charge of the stored energy, SOC max Representing a maximum value of the stored state of charge;
distributed photovoltaic power generation output constraint:
η i ≥η min
wherein,representing the capacity, η, of a node i distributed photovoltaic i Represents the photovoltaic utilization of node i of period T, < >>Represents the theoretical output force, eta of the node i distributed photovoltaic at the moment t min Representing photovoltaic utilization limit,/->Representing the active output of the node i distributed photovoltaic; />Representing reactive power of the node i distributed photovoltaic;
load constraint;
wherein,active power demand, denoted node i time t, < >>Representing reactive power demand at time t of node i, p i (t) represents active power injection at nodes i, j, q i (t) represents reactive power injection at nodes i, j;
voltage amplitude constraint;
U max ≤U i,t ≤U min
wherein U is max Representing the upper limit of the node voltage amplitude, U min Representing a node voltage magnitude lower limit;
the full disturbance trajectory is defined as:
wherein R represents a natural number, N ξ Representing the total number of disturbances.
5. The distributed energy distributed control method based on local computing as claimed in claim 4, wherein: the method for simplifying the control method through the radiation function, constructing the polyhedral internal approximation of the feasible domain of the model, and solving the optimal solution upper bound of the decentralized control model comprises the following steps:
substituting the convex quadratic constraint with the linear constraint, deriving a polyhedral approximation of the convex set U (ζ):
deterministic constantExpressed as:
according to x Fx+ u Fu+ ξ FAnd the problem of optimal decentralized control is approximately expressed as that the xi is less than or equal to 0:
minimize E γ [c′x+(L u u+L ξ ξ)′Σ(L u u+L ξ ξ)]
subject toγ∈Γ
x Fx+ u Fu+ ξ Fξ≤0
x=x(0)+Bu
u=γ(ξ)
wherein the coefficient matrixSpecifying with underlying problem data, m representing the number of this inequality constraint, N u Representing the total number of inputs, N x The total number of states is indicated and, x F、 u F、 ξ Fa lower coefficient matrix respectively representing a state variable x, an input variable u and a random variable xi;
a subspace of a feasible decentralized affine control strategy is defined based on the following:
wherein S represents a distributed control information constraint set, and Q ε Γ requires a linear operator Q:satisfying the information structure embedded in the feasible decentralized control strategy set Γ;
minimize Tr((Q′L′ u ΣL u Q+(2L′ ξ ΣL u +e 1 c′B)Q+L′ ξ ΣL ξ )M)+c′x(0)
subject to Q∈S
Π∈R l×m ,/>
( u F+ x FB)Q+ x Fx(0)e′ 1 + ξ F+Z=0
Z=νe′ 1 +Π′W
wherein e 1 = (1, 0, …, 0) represents N ξ Unit vector of J in The optimal solution Z is represented as an auxiliary variable,representing a positive cone, ±representing a semi-positive fixed cone, v representing the input of the system, M representing a second moment matrix of the system disturbance, where m=e (ζζ') is both reversible and positive, where E (·) represents the desire.
6. The distributed energy distributed control method based on local computing as claimed in claim 5, wherein: the construction of the outer polyhedron approximation of the feasible domain of the problem, solving the optimal solution lower bound of the distributed control model, and acquiring the approximate optimal solution based on the upper bound and the lower bound comprises;
for node i=1..n convex quadratic constraint is equivalent to a linear inequality:
according toApproximating optimal decentralized control as:
minimize Tr((Q′L′ u ΣL u Q+(2L′ ξ ΣL u +e 1 c′B)Q+L′ ξ ΣL ξ )M)+c′x(0)
subject to Q∈S,
e′ 1 MZ′≥0
wherein Tr (·) represents the trace; j (J) out Representing an optimal solution, which is an optimal solution lower bound of optimal decentralized control;an upper coefficient matrix respectively representing a state variable x, an input variable u and a random variable xi;
solving to obtain:
J out ≤J opt ≤J in
will J opt As a choice of an optimal solution.
7. A system employing the distributed energy distributed control method for a power distribution network based on local computing as claimed in any one of claims 1 to 6, characterized in that:
the model construction module is used for acquiring the resource information of the power grid and constructing a distributed control objective function and constraint conditions of the power distribution network; constructing an overall dynamic time state space model, and establishing a completely dispersed disturbance feedback control strategy;
the analysis module is used for constructing a polyhedron internal approximation of a feasible domain of the model through a radiation function simplified control method and solving an optimal solution upper bound of the decentralized control model; constructing an outer polyhedron approximation of a problem feasible domain, solving an optimal solution lower bound of a distributed control model, and acquiring an approximate optimal solution based on the upper bound and the lower bound.
8. A computer device, comprising: a memory and a processor; the memory stores a computer program characterized in that: the processor, when executing the computer program, implements the steps of the method of any one of claims 1 to 6.
9. A computer-readable storage medium having stored thereon a computer program, characterized by: the computer program implementing the steps of the method of any of claims 1 to 6 when executed by a processor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311718732.XA CN117410990B (en) | 2023-12-14 | 2023-12-14 | Distributed energy distributed control method and system based on local calculation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311718732.XA CN117410990B (en) | 2023-12-14 | 2023-12-14 | Distributed energy distributed control method and system based on local calculation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117410990A CN117410990A (en) | 2024-01-16 |
CN117410990B true CN117410990B (en) | 2024-02-23 |
Family
ID=89494824
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311718732.XA Active CN117410990B (en) | 2023-12-14 | 2023-12-14 | Distributed energy distributed control method and system based on local calculation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117410990B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111507598A (en) * | 2020-04-08 | 2020-08-07 | 华北电力大学 | Power distribution system reliability calculation method considering demand side resource layering dispersion control |
CN113285461A (en) * | 2021-05-26 | 2021-08-20 | 合肥工业大学 | Power distribution network decentralized control method and system based on incomplete measurement |
CN116777234A (en) * | 2023-05-22 | 2023-09-19 | 中国电力科学研究院有限公司 | Method and system for constructing two-stage robust optimization model and formulating strategy |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112904731A (en) * | 2021-01-29 | 2021-06-04 | 清华大学 | Stability domain determination method for multi-time-lag distributed power information physical system |
-
2023
- 2023-12-14 CN CN202311718732.XA patent/CN117410990B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111507598A (en) * | 2020-04-08 | 2020-08-07 | 华北电力大学 | Power distribution system reliability calculation method considering demand side resource layering dispersion control |
CN113285461A (en) * | 2021-05-26 | 2021-08-20 | 合肥工业大学 | Power distribution network decentralized control method and system based on incomplete measurement |
CN116777234A (en) * | 2023-05-22 | 2023-09-19 | 中国电力科学研究院有限公司 | Method and system for constructing two-stage robust optimization model and formulating strategy |
Non-Patent Citations (2)
Title |
---|
Finite-Time Distributed Convex Optimization for Continuous-Time Multiagent Systems With Disturbance Rejection;Zhi Feng 等;《IEEE Transactions on Control of Network Systems》;20190905;第7卷(第2期);686 - 698 * |
分布式随机模型预测控制方法研究;戴荔;《中国博士学位论文全文数据库信息科技辑》;20160715(第07期);I140-1 * |
Also Published As
Publication number | Publication date |
---|---|
CN117410990A (en) | 2024-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Oshnoei et al. | Robust control scheme for distributed battery energy storage systems in load frequency control | |
Ruan et al. | Distributed voltage control in active distribution network considering renewable energy: A novel network partitioning method | |
Chen et al. | Consensus-based distributed control for photovoltaic-battery units in a DC microgrid | |
Rahman et al. | Optimization of virtual inertia considering system frequency protection scheme | |
CN112163304B (en) | Power transmission network redundancy constraint identification method, storage medium and computing equipment | |
Wang et al. | Two-stage full-data processing for microgrid planning with high penetrations of renewable energy sources | |
Chen et al. | A unified time scale intelligent control algorithm for microgrid based on extreme dynamic programming | |
CN112418496B (en) | Power distribution station energy storage configuration method based on deep learning | |
CN114430165A (en) | Micro-grid group intelligent coordination control method and device based on depth model prediction | |
CN115313519A (en) | Power distribution network energy storage optimal configuration method, device, equipment and storage medium | |
CN113872213B (en) | Autonomous optimization control method and device for power distribution network voltage | |
Hassanzadeh et al. | Hierarchical optimal allocation of BESS using APT-FPSO based on stochastic programming model considering voltage sensitivity and eigenvalues analyses | |
Zhang et al. | Distributed coordination regulation marginal cost control strategy for automatic generation control of interconnected power system | |
Akbari-Zadeh et al. | Dstatcom allocation in the distribution system considering load uncertainty | |
Yang et al. | Distributed optimal generation and load-side control for frequency regulation in power systems | |
CN109638837B (en) | New energy grid-connected-oriented dynamic load management method and system | |
CN117410990B (en) | Distributed energy distributed control method and system based on local calculation | |
CN112633699A (en) | Active power distribution network frame planning method and device, computer equipment and storage medium | |
CN117035335A (en) | Multi-stage energy storage and power transmission network collaborative planning method and system | |
Skiparev et al. | Reinforcement learning based MIMO controller for virtual inertia control in isolated microgrids | |
Hua et al. | Stochastic optimal control scheme for operation cost management in energy Internet | |
Yang et al. | Coordinated restoration of inverter‐based power sources and synchronous generators for the high renewable penetrated power system considering the dynamic frequency regulation capability | |
de Jongh et al. | Neural Predictive Control for the Optimization of Smart Grid Flexibility Schedules | |
Zhou | Simulation of photovoltaic absorption strategy for distribution network considering translational load | |
Yin et al. | Lazy deep Q networks for unified rotor angle stability framework with unified time-scale of power systems with mass distributed energy storage |
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 | ||
GR01 | Patent grant |