CN115276038A - Distributed alternating current micro-grid frequency recovery method and system - Google Patents

Distributed alternating current micro-grid frequency recovery method and system Download PDF

Info

Publication number
CN115276038A
CN115276038A CN202210944156.XA CN202210944156A CN115276038A CN 115276038 A CN115276038 A CN 115276038A CN 202210944156 A CN202210944156 A CN 202210944156A CN 115276038 A CN115276038 A CN 115276038A
Authority
CN
China
Prior art keywords
distributed
local
information
power generation
generation unit
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
CN202210944156.XA
Other languages
Chinese (zh)
Other versions
CN115276038B (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.)
Zhejiang University ZJU
Hainan Institute of Zhejiang University
Original Assignee
Zhejiang University ZJU
Hainan Institute of Zhejiang University
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 Zhejiang University ZJU, Hainan Institute of Zhejiang University filed Critical Zhejiang University ZJU
Priority to CN202210944156.XA priority Critical patent/CN115276038B/en
Publication of CN115276038A publication Critical patent/CN115276038A/en
Application granted granted Critical
Publication of CN115276038B publication Critical patent/CN115276038B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241The oscillation concerning frequency
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00002Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by monitoring
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00004Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by the power network being locally controlled
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00006Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

The invention discloses a distributed alternating current micro-grid frequency recovery method and a distributed alternating current micro-grid frequency recovery system, which relate to the technical field of power generation control and specifically comprise the following steps: collecting three-phase voltage information, current information, droop controller output angular frequency information and voltage, current and angular frequency information of a neighbor distributed power generation unit connected with a directional communication topology; calculating a distributed error index according to the acquired information; judging whether the local distributed power generation unit is attacked or not; if yes, each distributed power generation unit solves local distributed optimal auxiliary control input based on a system equation and an integral performance index equation of the distributed power generation unit; calculating an optimal angular frequency reference value; the invention adopts a complete distribution type, each DG has a distributed performance index belonging to the DG, global information is not needed, information is exchanged with adjacent DGs only based on a sparse communication network, consensus is achieved in a cooperation mode, and the requirement on a centralized performance index is eliminated.

Description

Distributed alternating current micro-grid frequency recovery method and system
Technical Field
The invention relates to the technical field of power generation control, in particular to a distributed alternating current micro-grid frequency recovery method and system.
Background
With the massive access of distributed energy and the deep application of information technology in the smart grid, the information flow and the energy flow in the active power distribution network are deeply merged, so that the information flow becomes a typical information physical system (CPS). Due to the fact that CPS marketing behaviors of the active power distribution network are gradually widened, interaction between users and the power distribution network is increasingly frequent, and the active power distribution network faces a severe challenge in terms of dealing with information safety risks. The ac microgrid is a component of an active power distribution network CPS, which is connected to a main network in the form of an autonomous system via a Point of Common Coupling (PCC). When the micro-grid operates in a grid-connected mode, the voltage support and the frequency synchronization of the micro-grid are strictly controlled by the main grid. When the grid fails, the microgrid will quickly switch to island mode operation. Due to the loss of the support of the voltage and the frequency of the power grid, the power grid must be autonomous by a control system to realize the synchronization of the voltage and the frequency.
Conventional microgrid systems are typically controlled using a hierarchical control architecture. The main control is realized locally in the power supply unit, the operation is carried out under a rapid time scale, only difference adjustment can be carried out, the difference adjustment is used for balancing the generated energy and the load demand of the microgrid, and droop control is generally adopted. The quadratic control operates on a medium time scale to compensate for voltage and frequency deviations generated by the main control in case of mode switching or failure. And the third control is operated under a longer time scale, and the micro-grid is controlled to operate under a better economic condition according to the flow between the micro-grid and the main grid. The secondary control ensures that the micro-grid operates under a nominal condition, and plays a vital role in improving the electric energy quality of a user side or a load side.
Centralized control relies on a central controller (MGCC) and a highly interconnected communication network, and the control system is robust. However, the scalability, flexibility and reconfigurability of the centralized control system are poor, and plug and play cannot be realized. When the MGCC is attacked by malicious information, the overall paralysis of the micro-grid can be caused. The distributed control does not need a communication system, can greatly reduce the complexity of the system and has higher convergence speed. However, the control method cannot optimize the output of each unit generally, and it is difficult to realize consistent coordination of distributed generation units (DG). Therefore, how to guarantee consistent cooperation of DGs by using a sparse communication network is an urgent problem to be solved by those skilled in the art.
Disclosure of Invention
In view of this, the present invention provides a distributed ac microgrid frequency recovery method and system, which overcome the above-mentioned drawbacks.
In order to achieve the above purpose, the invention provides the following technical scheme:
a frequency recovery method based on a distributed alternating current micro-grid comprises the following specific steps:
collecting three-phase voltage information v of local alternating current microgrid bus oi And three-phase current information i oi Local droop controller outputting angular frequency information omega i And the voltage v of the neighboring distributed generation units connected to the communication topology oj Current i oj Angular frequency information omega j
According to the three-phase voltage information v of the local AC micro-grid bus oi Three-phase current information i oi Local droop controller outputting angular frequency information omega i And the voltage v of the neighboring distributed generation units connected to the communication topology oj Current i oj Angular frequency information omega j Calculating a distributed error index e i
According to a distributed error index e i Judging the attack condition of the local distributed power generation unit;
if any distributed power generation unit is attacked, each distributed power generation unit collaboratively solves local distributed optimal auxiliary control input based on the system equation and integral performance index of each distributed power generation unit
Figure BDA0003783795010000021
Based on locally distributed optimal auxiliary control input
Figure BDA0003783795010000022
Calculating an optimal angular frequency reference value
Figure BDA0003783795010000023
Optionally, the distributed error index e of any distributed power generation unit i The calculation formula of (2) is as follows:
Figure BDA0003783795010000024
in the formula, a i,j Is an element of the adjacency matrix; omega i Outputting angular frequency information for a local droop controller; omega j Angular frequency information for adjacent distributed power generation units; y is i Outputting a value for the local system; omega ref Is a reference voltage; b i Connecting elements of the gain matrix for the reference values; j is a neighbor node.
Optionally, the system equation of any distributed power generation unit is:
Figure BDA0003783795010000031
in the formula (I), the compound is shown in the specification,
Figure BDA0003783795010000032
is the first derivative of the system state; x is the number of i Is the local system state; y is i Outputting a value for the local system; u. of i A local auxiliary control input; i is a local node; A. b and C are both system matrixes.
Optionally, a local auxiliary control input u i The expression of (a) is:
Figure BDA0003783795010000033
in the formula, R is a positive definite constant matrix; r is i Is a positive scalar quantity; i is a local node; b is a system matrix; x is the number of n Is a state variable of the nth distributed generation unit; p i,n Lambda is a multiplier of Langerday i A matrix of coefficients related to the state variables; xi i Is a lambd of a lambertian multiplier i A matrix of coefficients independent of the state variables.
Optional, locally distributed optimal auxiliary control input
Figure BDA0003783795010000034
The acquisition steps are as follows:
constructing a system equation and an integral performance index equation of any distributed power generation unit, and constructing an integral performance index J i The expression of (c) is:
Figure BDA0003783795010000035
in the formula, e i Is a local error index; u. of i A local auxiliary control input; q. q.s i 、r i Are all positive scalars; q and R are both positive definite constant matrixes;
introducing lagrange multiplier lambda i Obtaining local auxiliary control input u based on system equation and integral performance index equation of each distributed power generation unit i The expression of (2);
for lagrange multiplier lambda i System equation and integral performance index equation of each distributed power generation unit and auxiliary control input u i The expression of (2) is subjected to multiple combined operations to obtain an optimal state track;
computing local distributed optimal auxiliary control inputs based on optimal state trajectories
Figure BDA0003783795010000041
Optionally, any optimum angular frequency reference value
Figure BDA0003783795010000042
The calculation formula of (2) is as follows:
Figure BDA0003783795010000043
in the formula u i A local auxiliary control input; d P,i Is a droop control coefficient; p is i The output active power for the local VSI.
A frequency recovery system based on a distributed alternating current micro-grid comprises a communication network, a state collector, an attack detection module, a distributed optimal cooperative controller, a power controller, a voltage controller, a current controller and a voltage source type inverter;
a communication network for transmission of data;
a state collector for collecting three-phase voltage information v of local AC microgrid bus oi Three-phase current information i oi Local droop controller outputting angular frequency information omega i And the voltage v of the adjacent distributed generation unit of any one of the distributed generation units oj Current i oj Angular frequency information omega j
An attack detection step, which is used for judging whether each distributed power generation unit is attacked or not and sending a control signal to the distributed optimal cooperative controller;
distributed optimal cooperative controller for computing local distributed optimal auxiliary control input
Figure BDA0003783795010000044
And an optimum angular frequency reference value
Figure BDA0003783795010000045
A power controller for controlling the power according to the optimal angular frequency reference value
Figure BDA0003783795010000046
Completing the compensation of the frequency deviation;
voltage controller ofAccording to a reference voltage
Figure BDA0003783795010000047
Deriving a reference current for an input current controller
Figure BDA0003783795010000048
A current controller for controlling the current according to the reference current
Figure BDA0003783795010000049
Obtaining a reference voltage for controlling a PWM generator to generate a PWM signal according to the output current of the voltage source type inverter;
voltage source type inverter for controlling angular frequency omega according to PWM signal and power controller i And outputting the voltage.
Optionally, the system is divided into a physical layer and an information layer;
the physical layer comprises n distributed power generation units, and the n distributed power generation units are connected to a local microgrid bus through a voltage source type inverter;
the information layer is a directed connection communication topology containing a spanning tree.
Optionally, the directional connection communication topology is represented by G (V, E, a, B), where V is a node set; e is a directed edge set; a is an adjacency matrix of the topology; b is the connection gain matrix.
Optionally, the microgrid comprises n distributed power generation units and m integrated loads.
According to the technical scheme, the invention discloses a distributed AC micro-grid frequency recovery method and a distributed AC micro-grid frequency recovery system, and compared with the prior art, the distributed AC micro-grid frequency recovery method and the distributed AC micro-grid frequency recovery system have the following advantages:
1. the invention adopts complete distribution, each DG has the own distributed performance index, global information is not needed, information is exchanged with the adjacent DGs only on the basis of a sparse communication network, consensus is achieved in a cooperation mode, and the requirement on the centralized performance index is eliminated;
2. the invention has universality to different communication topologies;
3. the method provided by the invention can realize the optimal frequency consistency recovery of each distributed generator set under the attack of the microgrid by using the minimum control energy, and the convergence speed is accelerated by considering the synchronism among the generator sets.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the prior art descriptions will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 is a schematic diagram of the structure of the information layer and the physical layer of the present invention;
FIG. 2 is a schematic diagram of a four-unit two-integrated-load microgrid topology structure according to the present invention;
FIG. 3 is a schematic diagram of the communication network topology of FIG. 2;
FIG. 4 is a schematic diagram of a single distributed generation unit local control system configuration of the present invention;
FIG. 5 is a schematic flow chart of the method of the present invention.
Detailed Description
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 only a part of the embodiments of the present invention, and not all of the embodiments. 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 discloses a distributed alternating current microgrid frequency recovery system, which is divided into a physical layer and an information layer as shown in figure 1;
and in the physical layer, the n distributed generation units are connected to the microgrid bus through Voltage Source Inverters (VSIs). And circuit equivalent impedance is formed between the buses.
And the information layer comprises a sparse and efficient directed connection communication topology of the spanning tree. The directed connectivity topology can be represented by G (V, E, A, B); wherein V = { V = 1 ,v 2 ,...,v n The is a node set; e is a directed edge set; n is a radical of hydrogen i Representing a set of neighbor nodes, N i ={v j |(i,j)∈E};
Figure BDA0003783795010000061
Is a contiguous matrix of the topology, A = [ a = [) i,j ]Wherein a is i,j When the directional communication is carried out, the value is 1, otherwise, the value is 0;
Figure BDA0003783795010000062
to connect the gain matrix, wherein b if node i can receive the information of the leader node i The value is 1, otherwise 0. The sparse and efficient communication network involved in the embodiment is a directed connectivity topology including a spanning tree, and meets the condition of inducing multi-MAS cooperation consensus.
The graph Laplace matrix of the communication topology is
Figure BDA0003783795010000063
The microgrid in this embodiment is composed of n distributed power generation units and m integrated loads, as shown in fig. 2, the microgrid is a microgrid with four units and two integrated loads, and the communication topology of the microgrid may be six or more than six in fig. 3, as shown in fig. 3.
For the inside of each distributed power generation unit, a block diagram of frequency recovery optimal distributed cooperative control under malicious attack is shown in fig. 4, and comprises a communication network, a state collector, an attack detection module, a distributed optimal cooperative controller, a power controller, a voltage controller, a current controller and a voltage source type inverter;
a communication network for transmission of data; a state collector for collecting three-phase voltage information v of local AC microgrid bus oi Three-phase current information i oi Local droop controller output angular frequency information omega i And the voltage v of the adjacent distributed generation unit of any one of the distributed generation units oj Current i oj Angular frequency information omega j (ii) a An attack detection step, which is used for judging whether each distributed power generation unit is attacked or not and sending a control signal to the distributed optimal cooperative controller; distributed optimal cooperative controller for computing local distributed optimal auxiliary control input
Figure BDA0003783795010000071
And an optimum angular frequency reference value
Figure BDA0003783795010000072
A power controller for controlling the power according to the optimal angular frequency reference value
Figure BDA0003783795010000073
Completing the compensation of the frequency deviation; a voltage controller for controlling the voltage according to the reference voltage
Figure BDA0003783795010000074
Deriving a reference current for an input current controller
Figure BDA0003783795010000075
A current controller for controlling the current according to the reference current
Figure BDA0003783795010000076
Obtaining a reference voltage for controlling a PWM generator to generate a PWM signal according to the output current of the voltage source type inverter; voltage source type inverter for controlling angular frequency omega of power controller according to PWM signal i And outputting the voltage.
In another embodiment, a state estimator is further included for estimating states of neighboring nodes.
The embodiment further includes a frequency recovery method based on the distributed ac microgrid, which is characterized in that, as shown in fig. 5, the specific steps are as follows:
step 1, collecting three-phase voltage information v of local alternating current micro-grid bus oi Three-phase current information i oi Local droop controller outputting angular frequency information omega i Andvoltage v of neighbor distributed generation units connected in directed communication topology oj Current i oj Angular frequency information omega j
Step 2, according to the three-phase voltage information v of the local alternating current micro-grid bus oi Three-phase current information i oi Local droop controller output angular frequency information omega i And voltage v of a neighboring distributed generation unit connected to a directed communication topology oj Current i oj Angular frequency information omega j Calculating a distributed error index e i (ii) a The method comprises the following specific steps:
the secondary frequency recovery uses an optimal distributed cooperation strategy. A state collector for collecting the three-phase voltage and current information v of the local AC microgrid bus oi 、i oi And angular frequency information omega output by the local droop controller i (ii) a And according to the voltage voj, the current ioj and the angular frequency information omega of the adjacent distributed power generation units j And calculating any distributed error index e by the following formula i
Figure BDA0003783795010000081
In the formula, omega ref For reference voltage, 50Hz was taken. Each DG utilizes its own information and neighbor information to complete frequency recovery with a high degree of consistency.
Step 3, according to the distributed error index e i Judging the state of any distributed power generation unit under attack; the method specifically comprises the following steps:
distributed error index e i And simultaneously transmitting an attack detection link and a distributed optimal cooperative controller. And in the attack detection link, the condition that the frequency fluctuation exceeds +/-0.2 Hz is considered as being attacked or failed, when the bus frequency of the local distributed power generation unit exceeds the limit, the distributed optimal cooperative controller is automatically triggered, and the frequency cooperative recovery is started together with the n distributed power generation units in the microgrid.
Step 4, if the local distributed power generation units are attacked, each distributed power generation unit is based on the system equation and the integral performance index of the distributed power generation unitCollaborative solving of locally distributed optimal auxiliary control inputs
Figure BDA0003783795010000082
The method specifically comprises the following steps:
the system equation and the integral performance index equation of any distributed power generation unit are as follows:
Figure BDA0003783795010000083
wherein
Figure BDA0003783795010000084
Are all positive definite parameter matrices, q i And r i Are all positive scalar parameters. The performance index is divided into two parts, the first part is the cost of the control process, and the second part is the control input energy. Obviously, to achieve the control objective, it is desirable to minimize both process cost and control energy. e.g. of the type i Is a local error index; u. of i Is a local auxiliary control input.
Carrying out input and output feedback linearization on a local control system of any distributed power generation unit, and constructing the following system equation:
Figure BDA0003783795010000091
the system comprises n DGs, i =1, 2.., n,
Figure BDA0003783795010000092
and
Figure BDA0003783795010000093
respectively, system state, system output, and control input.
Figure BDA0003783795010000094
Is a system matrix.
The idea of solving the functional extreme value with constraint condition according to the variational method is introducedMultiplier lambda i And converting the obtained solution into an unconstrained lagrangian function extremum problem for solving, wherein the unconstrained lagrangian function extremum problem specifically comprises the following steps:
Figure BDA0003783795010000095
Figure BDA0003783795010000096
λ i (∞)=0 (6);
Figure BDA0003783795010000097
wherein, the formula (4) and the formula (5) are adjoint equations, the formula (6) is an edge condition, and the formula (7) is a control equation. N is satisfied for all i =1,2 i Hamiltonian for ith DG
Figure BDA0003783795010000098
Will be lambda i Substitution of formula (7) gives:
Figure BDA0003783795010000099
for lambda i Derivation of
Figure BDA00037837950100000910
Further, the following equations (4) and (8) are substituted to obtain:
Figure BDA0003783795010000101
after substituting formula (9) into formula (5), obtaining:
Figure BDA0003783795010000102
wherein L is a Laplacian matrix of the system communication topology,
Figure BDA0003783795010000103
i is a local node; j is a neighbor node.
This formula is satisfied for all i =1,2. Wherein
Figure BDA0003783795010000104
Is a matrix with 1 for the ith entry and 0 for the remaining entries. Since P and ξ are time-invariant matrices, so
Figure BDA0003783795010000105
And
Figure BDA0003783795010000106
are all 0. The equation of formula (10) is split into two parts on both sides: and [ x ] 1 x 2 …x n ] T Relevant moiety, [ x ] 1 x 2 …x n ] T An extraneous portion. The following two equations are constructed:
Figure BDA0003783795010000107
Figure BDA0003783795010000108
obtaining P from each of the equations (11) and (12) i,1 ,P i,2 ,...,P i,n And xi i Further carrying out the step (8) to obtain u i . Then u is removed i (i =1, 2.., n) the belt-in (4) yields the following differential equation about the state of the system:
Figure BDA0003783795010000111
solving the differential equation to obtain the optimal state track
Figure BDA0003783795010000112
Will P i,1 ,P i,2 ,...,P i,n 、ξ i And
Figure BDA0003783795010000113
surrogate formula (8) to obtain local distributed optimal auxiliary control input
Figure BDA0003783795010000114
Step 5, inputting according to local distributed optimal auxiliary control
Figure BDA0003783795010000115
Calculating an optimal angular frequency reference value
Figure BDA0003783795010000116
Reference value of optimal angular frequency
Figure BDA0003783795010000117
And inputting the droop controller so as to complete the compensation of the frequency deviation of the primary controller.
Figure BDA0003783795010000118
The acquisition steps are as follows:
the primary control uses a local controller of a single distributed power generation unit, including droop control and voltage-current double closed-loop control. Three-phase voltage and current information v acquired from alternating current micro-grid bus oi 、i oi Then, the power is input to a power controller after being subjected to abc/dq orthogonal decomposition. The power controller comprises a power calculating device and a droop controller, wherein the power calculating device calculates and obtains active power and reactive power P of the bus end of the microgrid according to input voltage and current information i 、Q i . The droop controller is according to the droop expression:
Figure BDA0003783795010000119
obtaining a reference voltage input to a voltage controller
Figure BDA00037837950100001110
And angular frequency omega inputted to the voltage source type inverter i . The voltage controller is based on the reference voltage
Figure BDA00037837950100001111
Deriving a reference current for an input current controller
Figure BDA00037837950100001112
The current controller is based on the reference current
Figure BDA00037837950100001113
And the output current of the voltage source type inverter obtains a reference voltage for controlling the PWM generator to generate the PWM signal. Further based on the PWM signal and the angular frequency omega given by the droop controller i Controlling the inverter output.
The frequency deviation caused by the droop control needs to be compensated by designing a secondary control. To obtain the auxiliary control input for the distributed optimal cooperative controller, equation (14) is differentiated to obtain the following differential equation:
Figure BDA0003783795010000121
Figure BDA0003783795010000122
and
Figure BDA0003783795010000123
are respectively omega i
Figure BDA0003783795010000124
And P i A differential of (c). u. u i Is the auxiliary control input to the ith distributed generation unit.
Based on the equations (14) and (15), it is possible to obtain:
Figure BDA0003783795010000125
the method provided by the invention is completely distributed, each DG has a distributed performance index, global information is not needed, information is exchanged with adjacent DGs only on the basis of a sparse communication network, consensus is achieved in a cooperative mode, and the requirement on a centralized performance index is eliminated; the method has universality to different communication topologies; the method provided by the invention can realize the optimal frequency recovery of each distributed generator set under the attack of the microgrid by using the minimum control energy. And the convergence speed is accelerated by considering the synchronism among the units.
In the present specification, the embodiments are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (10)

1. A frequency recovery method based on a distributed alternating current microgrid is characterized by comprising the following specific steps:
collecting three-phase voltage information v of local AC micro-grid bus oi Three-phase current information i oi Local droop controller output angular frequency information omega i And voltage v of a neighboring distributed generation unit connected to a directed communication topology oj Current i oj Angular frequency information omega j
According to local exchange microThree-phase voltage information v of power grid bus oi And three-phase current information i oi Local droop controller output angular frequency information omega i And the voltage v of the neighboring distributed generation units connected to the communication topology oj Current i oj Angular frequency information omega j Calculating a distributed error indicator e i
According to a distributed error index e i Judging the attack condition of the local distributed power generation unit;
if any distributed power generation unit is attacked, each distributed power generation unit collaboratively solves local distributed optimal auxiliary control input based on the system equation and integral performance index of each distributed power generation unit
Figure FDA0003783789000000011
Based on locally distributed optimal auxiliary control input
Figure FDA0003783789000000012
Calculating an optimal angular frequency reference value
Figure FDA0003783789000000013
2. The frequency recovery method based on the distributed alternating current micro-grid according to claim 1, wherein the distributed error index e of any distributed power generation unit i The calculation formula of (2) is as follows:
Figure FDA0003783789000000014
in the formula, a i,j Is an element of the adjacency matrix; omega i Outputting angular frequency information for a local droop controller; omega j Angular frequency information for adjacent distributed power generation units; y is i Outputting a value for the local system; omega ref Is a reference voltage; b is a mixture of i Connecting elements of the gain matrix for the reference values; j is a neighbor node.
3. The frequency recovery method based on the distributed alternating current micro-grid according to claim 1, wherein a system equation of any distributed power generation unit is as follows:
Figure FDA0003783789000000015
in the formula (I), the compound is shown in the specification,
Figure FDA0003783789000000021
is the first derivative of the system state; x is the number of i Is the local system state; y is i Outputting a value for the local system; u. u i A local auxiliary control input; i is a local node; A. b and C are both system matrixes.
4. The method of claim 3, wherein the local auxiliary control input u is a frequency recovery input i The expression of (a) is:
Figure FDA0003783789000000022
in the formula, R is a positive definite constant matrix; r is i Is a positive scalar quantity; i is a local node; b is a system matrix; x is the number of n Is a state variable of the nth distributed generation unit; p is i,n Is a lambd of a lambertian multiplier i A matrix of coefficients relating to the state variables; xi shape i Lambda is a multiplier of Langerday i A matrix of coefficients independent of the state variables.
5. The method of claim 1, wherein the local distributed optimal auxiliary control input is used as a frequency recovery input for the distributed ac microgrid
Figure FDA0003783789000000023
The step of obtainingThe method comprises the following steps:
constructing a system equation and an integral performance index equation of any distributed power generation unit, and constructing an integral performance index J i The expression of (c) is:
Figure FDA0003783789000000024
in the formula, e i Is a local error index; u. u i A local auxiliary control input; q. q of i 、r i Are all positive scalars; q and R are both positive definite constant matrixes;
introducing lagrange multiplier lambda i Obtaining local auxiliary control input u based on system equation and integral performance index equation of each distributed power generation unit i The expression of (1);
for lagrange multiplier lambda i System equation and integral performance index equation of each distributed power generation unit and auxiliary control input u i The expression of (2) is subjected to multiple combined operations to obtain an optimal state track;
computing local distributed optimal auxiliary control input based on optimal state trajectory
Figure FDA0003783789000000025
6. The frequency recovery method based on the distributed alternating current micro-grid according to claim 1, wherein any optimal angular frequency reference value
Figure FDA0003783789000000026
The calculation formula of (2) is as follows:
Figure FDA0003783789000000031
in the formula u i A local auxiliary control input; d P,i Is a droop control coefficient; p i Active for output of local VSIAnd (3) power.
7. A frequency recovery system based on a distributed alternating current micro-grid is characterized by comprising a communication network, a state collector, an attack detection module, a distributed optimal cooperative controller, a power controller, a voltage controller, a current controller and a voltage source type inverter;
a communication network for transmission of data;
the state collector is used for collecting three-phase voltage information v of a local alternating current micro-grid bus oi Three-phase current information i oi Local droop controller outputting angular frequency information omega i And the voltage v of the adjacent distributed generation unit of any one of the distributed generation units oj Current i oj Angular frequency information omega j
An attack detection step, which is used for judging whether each distributed power generation unit is attacked or not and sending a control signal to the distributed optimal cooperative controller;
distributed optimal cooperative controller for computing local distributed optimal auxiliary control input
Figure FDA0003783789000000032
And an optimum angular frequency reference value
Figure FDA0003783789000000033
A power controller for controlling the power according to the optimal angular frequency reference value
Figure FDA0003783789000000034
Completing the compensation of the frequency deviation;
a voltage controller for controlling the voltage according to the reference voltage
Figure FDA0003783789000000035
Deriving a reference current for an input current controller
Figure FDA0003783789000000036
A current controller for controlling the current according to the reference current
Figure FDA0003783789000000037
Obtaining a reference voltage for controlling a PWM generator to generate a PWM signal according to the output current of the voltage source type inverter;
voltage source type inverter for controlling angular frequency omega of power controller according to PWM signal i And outputting the voltage.
8. The frequency recovery system based on the distributed alternating current micro-grid according to claim 7, wherein the system is divided into a physical layer and an information layer;
the physical layer comprises n distributed power generation units, and the n distributed power generation units are connected to a local microgrid bus through a voltage source type inverter;
the information layer is a directed connection communication topology containing a spanning tree.
9. The frequency recovery system based on the distributed ac microgrid of claim 8, wherein the directional connection communication topology is represented by G (V, E, a, B), wherein V is a node set; e is a directed edge set; a is an adjacency matrix of the topology; b is the connection gain matrix.
10. The frequency recovery system based on the distributed alternating current microgrid of claim 7, characterized in that the microgrid is composed of n distributed power generation units and m integrated loads.
CN202210944156.XA 2022-08-05 2022-08-05 Alternating-current micro-grid frequency recovery method and system based on distribution Active CN115276038B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210944156.XA CN115276038B (en) 2022-08-05 2022-08-05 Alternating-current micro-grid frequency recovery method and system based on distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210944156.XA CN115276038B (en) 2022-08-05 2022-08-05 Alternating-current micro-grid frequency recovery method and system based on distribution

Publications (2)

Publication Number Publication Date
CN115276038A true CN115276038A (en) 2022-11-01
CN115276038B CN115276038B (en) 2023-09-12

Family

ID=83748977

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210944156.XA Active CN115276038B (en) 2022-08-05 2022-08-05 Alternating-current micro-grid frequency recovery method and system based on distribution

Country Status (1)

Country Link
CN (1) CN115276038B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116886360A (en) * 2023-07-10 2023-10-13 浙江大学 Fully distributed island micro-grid attack elastic control system and method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2532888A1 (en) * 2011-06-08 2012-12-12 Siemens Aktiengesellschaft Arrangement for generating a control signal for controlling a power output of a power generation system
US20160276830A1 (en) * 2015-03-16 2016-09-22 Board Of Regents, The University Of Texas System System and Method for Distributed Control of an Electrical Network
CN112736975A (en) * 2020-12-30 2021-04-30 上海交通大学 Alternating current micro-grid droop coefficient optimization method
CN112769127A (en) * 2020-12-16 2021-05-07 浙江工业大学 Alternating current micro-grid frequency attack detection and recovery method based on distributed intermediate observer
CN114512983A (en) * 2022-03-02 2022-05-17 国网浙江省电力有限公司信息通信分公司 Distributed power supply elastic control method for network attack

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2532888A1 (en) * 2011-06-08 2012-12-12 Siemens Aktiengesellschaft Arrangement for generating a control signal for controlling a power output of a power generation system
US20160276830A1 (en) * 2015-03-16 2016-09-22 Board Of Regents, The University Of Texas System System and Method for Distributed Control of an Electrical Network
CN112769127A (en) * 2020-12-16 2021-05-07 浙江工业大学 Alternating current micro-grid frequency attack detection and recovery method based on distributed intermediate observer
CN112736975A (en) * 2020-12-30 2021-04-30 上海交通大学 Alternating current micro-grid droop coefficient optimization method
CN114512983A (en) * 2022-03-02 2022-05-17 国网浙江省电力有限公司信息通信分公司 Distributed power supply elastic control method for network attack

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Y. CHEN: "《Distributed Event-Triggered Secondary Control for Islanded Microgrids With Proper Trigger Condition Checking Period》", 《 IEEE TRANSACTIONS ON SMART GRID》, vol. 13, no. 2, pages 837 - 848 *
陈郁林 等: "《虚假数据注入攻击下的微电网分布式协同控制》", 《电力系统自动化》, vol. 45, no. 5, pages 97 - 103 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116886360A (en) * 2023-07-10 2023-10-13 浙江大学 Fully distributed island micro-grid attack elastic control system and method
CN116886360B (en) * 2023-07-10 2024-02-09 浙江大学 Fully distributed island micro-grid attack elastic control system and method

Also Published As

Publication number Publication date
CN115276038B (en) 2023-09-12

Similar Documents

Publication Publication Date Title
CN108363306B (en) Micro-grid distributed controller parameter determination method based on linear quadratic optimization
Bidram et al. A multiobjective distributed control framework for islanded AC microgrids
Lai et al. Distributed multiagent-oriented average control for voltage restoration and reactive power sharing of autonomous microgrids
Lu et al. Hierarchical distributed control approach for multiple on-site DERs coordinated operation in microgrid
CN108879690B (en) AC/DC hybrid microgrid data drive control method
CN105391071A (en) Multifunctional grid-connected parallel inverter group intelligent control method used in microgrid
CN110071514B (en) Consistency droop control method for power distribution and voltage frequency recovery
Raeispour et al. Robust Distributed Disturbance-Resilient $ H_\infty $-Based Control of Off-Grid Microgrids With Uncertain Communications
Nazari et al. Communication-failure-resilient distributed frequency control in smart grids: Part I: Architecture and distributed algorithms
Liu et al. Distributed reinforcement learning to coordinate current sharing and voltage restoration for islanded DC microgrid
CN115276038A (en) Distributed alternating current micro-grid frequency recovery method and system
Wu et al. Distributed optimal operation of hierarchically controlled microgrids
CN104779607A (en) Distributed coordination control method and system for direct-current micro-grid
Peng et al. Fully distributed discrete-time control of DC microgrids with zip loads
WO2023143284A1 (en) Distributed adaptive control method for microgrid energy storage system under unknown disturbance
Shotorbani et al. A distributed non-Lipschitz control framework for self-organizing microgrids with uncooperative and renewable generations
CN111987712B (en) Energy storage droop coefficient distributed control method for multi-voltage-level direct-current power distribution network
Zholbaryssov et al. Distributed enforcement of phase-cohesiveness for frequency control of islanded inverter-based microgrids
Doostinia et al. Distributed event-triggered average consensus control strategy with fractional-order local controllers for DC microgrids
Zuo et al. Enhanced power and energy coordination for batteries under the real-time closed-loop, distributed microgrid control
He et al. A novel control method based on consensus algorithm for microgrids
Yu et al. Optimal Cooperative Secondary Control for Islanded DC Microgrids via a Fully Actuated Approach
Lin et al. Distributed optimized consensus-based frequency synchronization in isolated ac microgrids
Zhao et al. Optimal distributed coordinated reinforcement learning for secondary voltage control in time-delayed microgrid
Aalipour et al. Proportional power-sharing control of distributed generators in microgrids

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
CB02 Change of applicant information

Address after: 310030 866 Yuhangtang Road, Xihu District, Hangzhou City, Zhejiang Province

Applicant after: ZHEJIANG University

Applicant after: Hainan Research Institute of Zhejiang University

Address before: 310030 866 Yuhangtang Road, Xihu District, Hangzhou City, Zhejiang Province

Applicant before: ZHEJIANG University

Applicant before: Research Institute of Hainan Zhejiang University

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant