Disclosure of Invention
The invention aims to provide an information physical cooperative frequency modulation method and device under multi-distributed photovoltaic polymerization, which solve the problems of low frequency modulation speed, poor effect and low reliability in the prior art.
In order to achieve the above purpose, the invention is realized by adopting the following technical scheme:
In a first aspect, the present invention provides a method for information physical cooperative frequency modulation under multi-distributed photovoltaic polymerization, including:
Acquiring a total active power adjustment quantity;
Optimizing a communication network among the aggregates by improving a consistency algorithm and a preset route optimization scheduling strategy;
based on the total active power adjustment amount and the optimized communication network, calculating an active power adjustment amount of the polymer through a leader-follower consistency algorithm;
And calculating the active power adjustment quantity of each photovoltaic source in the polymer through a centralized proportional distribution algorithm based on the active power adjustment quantity of the polymer, and regulating the output of the photovoltaic sources through the active power adjustment quantity of the photovoltaic sources to finish frequency modulation.
Further, the preset route optimization scheduling strategy comprises an optimization target and constraint conditions;
Optimization target:
;
wherein, Representing an objective function,Representing the number of polymers,Represents theAggregation of adjacent polymers of individual polymers,The representation is from thePolymers to theTime delay in channels of individual aggregates,Representing a sign function,Represents theIndividual Polymer andCommunication connection relationship between polymers, if No.Individual Polymer andIndividual polymers are then1 OtherwiseIs 0;
the expression of the sign function is:
;
Where r is the argument of the sign function;
Constraint conditions:
;
;
wherein, The representation is from theConnection of the individual polymers to the starting position,From the start position to theConnection of individual polymers to each other,A set of adjacent relay routes representing a starting location;
;
;
wherein, The representation is from theConnection of the individual Polymer to the destination position,Representing from the destination location to theConnection relation of individual polymers,A set of adjacent relay routes representing destination locations;
;
;
wherein, The representation is from theConnection relation of individual aggregate to first relay route,The representation is from theThe relay routes to theConnection relation of individual polymers,Represents theA set of adjacent relay routes of the plurality of relay routes;
;
;
wherein, Represents theThe upper limit of the number of input channels for each relay route,Represents theThe upper limit of the number of output channels of each relay route.
Further, the improved consistency algorithm comprises:
Will simultaneously satisfy the following formulas AndSet as/>, in a communication networkAnd:
;
;
Wherein,Is the control gain parameter between the jth and ith polymers,Is the control gain parameter between the leader and the ith polymer,Represents theRatio of active power adjustment quantity to rated capacity at initial time of individual polymer,Represents theRatio of active power adjustment quantity to rated capacity at initial time of individual polymer,Representing the communication connection relationship between the leader and the ith polymer, if the leader and the ith polymer are connected1 OtherwiseFor 0, ΔP represents the total active power adjustment amount,Represents the rated capacity of the ith polymer, max represents the maximum overshoot,Representing the ratio of the active power adjustment amount to the rated capacity in the steady state of the leader,The ratio of the active power adjustment amount to the rated capacity in the steady state of the ith polymer is shown.
Further, the calculating, based on the total active power adjustment and the optimized communication network, the active power adjustment of the aggregate by a leader-follower consistency algorithm includes calculating the active power adjustment of the aggregate by the following formula:
;
;
;
wherein, Represents the active power adjustment amount of the ith polymer,Derivative representing the ratio of the active power adjustment quantity to the rated capacity of the ith polymer at time t,Representing the control input of the ith polymer,Represents the ratio of the active power adjustment amount to the rated capacity of the jth polymer at time t,Represents the ratio of the active power adjustment amount to the rated capacity of the ith polymer at time t,Represents the ratio of the active power adjustment amount to the rated capacity of the ith polymer,Representing the ratio of the active power adjustment quantity and the rated capacity of the leader at the time t;
wherein, The expression of (2) is:
;
wherein, AndThe value of (2) satisfies the formula:
。
Further, the active power adjustment quantity of each photovoltaic source in the polymer is calculated by a centralized proportional distribution algorithm, and the calculation is performed by the following formula:
;
wherein, Represents the active power adjustment quantity of the mth photovoltaic source in the ith polymer,Representing the capacity scaling factor of the mth photovoltaic source in the ith polymer,The active power adjustment amount of the ith polymer at time t is shown.
In a second aspect, the present invention further provides an information physical cooperative frequency modulation device under multi-distributed photovoltaic polymerization, including:
A data acquisition module configured to: acquiring a total active power adjustment quantity;
A communication network optimization module configured to: optimizing a communication network among the aggregates by improving a consistency algorithm and a preset route optimization scheduling strategy;
An aggregate adjustment quantity distribution module configured to: based on the total active power adjustment amount and the optimized communication network, calculating an active power adjustment amount of the polymer through a leader-follower consistency algorithm;
A photovoltaic source conditioning amount distribution module configured to: and calculating the active power adjustment quantity of each photovoltaic source in the polymer through a centralized proportional distribution algorithm based on the active power adjustment quantity of the polymer, and regulating the output of the photovoltaic sources through the active power adjustment quantity of the photovoltaic sources to finish frequency modulation.
Compared with the prior art, the invention has the following beneficial effects:
According to the information physical cooperative frequency modulation method and device under the multi-distributed photovoltaic polymerization, a multi-distributed photovoltaic polymerization cooperative framework based on a polymer and a photovoltaic source is designed, a plurality of photovoltaic power generation systems can be divided into a plurality of layers to carry out frequency modulation control, the complexity of communication and the calculated amount of an upper control device are reduced, and the response speed of the system is improved; the leader-following consistency algorithm can effectively reduce power fluctuation among polymers and improve the frequency stability of the whole photovoltaic power generation system by carrying out state synchronization and active power adjustment on the multi-polymer system; the multi-distributed photovoltaic polymerization cooperative architecture is suitable for multi-distributed photovoltaic polymerization systems with different scales, can adapt to the change and expansion of system capacity, and has stronger expandability; the influence of communication time delay on the frequency modulation effect in the communication network is effectively solved by improving a consistency algorithm and a preset route optimization scheduling strategy, and the conservation of the frequency modulation method on the problem of uncertainty of communication is reduced.
Detailed Description
The present invention will be further described with reference to the accompanying drawings, and the following examples are only for more clearly illustrating the technical aspects of the present invention, and are not to be construed as limiting the scope of the present invention.
Example 1
As shown in fig. 2, the invention designs a multi-distributed photovoltaic polymerization cooperative architecture, which comprises an upper polymer layer and a lower light Fu Qun, and a photovoltaic group comprises a plurality of distributed photovoltaic sources.
As shown in fig. 1, the invention provides an information physical cooperative frequency modulation method under multi-distributed photovoltaic polymerization, which comprises the following steps:
S1, acquiring the total active power adjustment quantity.
The total active power adjustment is derived from an automatic power generation control (automation generation control, AGC) command received by the aggregate.
S2, optimizing the communication network among the aggregates through improving a consistency algorithm and a preset route optimization scheduling strategy.
Aiming at the upper layer polymers, the physical distance between the polymers is longer, and the data transmission time is longer, so the time delay problem needs to be considered and solved in regulation and control.
Taking the ith polymer as an example, under the interference of communication delay, the controlled expression of the ratio of the active power adjustment quantity to the rated capacity is as follows:
(1);
(2);
wherein, Representing the ratio of the active power adjustment amount to the rated capacity of the jth polymer at time,Representing the delay of the transmission of data from the jth aggregate to the ith aggregate,RepresentationRatio of active power adjustment to rated capacity of time leader,Representing the delay of the leader transmitting data to the ith aggregator,The ratio of the active power adjustment amount to the rated capacity of the ith polymer at time t is shown.
In order to ensure that excessive overshoot and steady-state error do not exist in the regulation and control process under the influence of time delay, a control gain parameter is setAndThe/>, of the ith polymer is as followsThe tracking error expression of (2) is as follows:
(3);
Then, based on Laplace transformation final value theorem and integral theorem pair Processing is carried out to obtain the following expression:
(4);
Wherein s represents complex variables in Laplace transformation, which are eliminated during processing, so that only proper control gain parameters need to be selected similar to a non-delay regulation strategy AndAnd (3) satisfying the formula (5) ensures that the output regulation and control among the upper polymers under the influence of time delay still is distributed according to the capacity proportion.
(5);
Wherein,Is the control gain parameter between the jth and ith polymers,Is the control gain parameter between the leader and the ith polymer,Representing the ratio of the active power adjustment amount to the rated capacity at the initial time of the first polymer,Representing the ratio of the active power adjustment amount to the rated capacity at the initial time of the first polymer,Representing the communication connection relationship between the leader and the ith polymer, if the leader and the ith polymer are connected1 OtherwiseFor 0, ΔP represents the total active power adjustment amount,The rated capacity of the ith polymer is shown.
However, in the case of meeting the above requirements, it is also ensured that an excessive overshoot cannot exist in the control process, so that an overshoot calculation is required to select an appropriate control gain parameter, which is specifically as follows:
since the overshoot is maximum at the initial time of the leader-follower consistency algorithm execution, the calculated expression of the overshoot in the ith polymer regulation process is obtained as follows:
(6);
Wherein max represents the maximum overshoot allowed during operation, Representing the ratio of the active power adjustment amount to the rated capacity in the steady state of the leader,The ratio of the active power adjustment amount to the rated capacity in the steady state of the ith polymer is shown.
Therefore, only proper control gain parameters need to be selectedAndAnd the regulation and control tasks among the upper polymers can be ensured to be successfully executed by satisfying the formula (7).
(7);
However, as can be seen from the formula (6), once no proper control gain parameter ensures that the overshoot is smaller than the allowable maximum overshoot, the consistency control effect is lost, so that the invention considers the route optimization scheduling to improve the transmission speed and reduce the communication delay, thereby reducing the overshoot, and a specific preset route optimization scheduling strategy is as follows:
Optimization target:
(8);
wherein, Representing an objective function,Representing the number of polymers,Represents theAggregation of adjacent polymers of individual polymers,The representation is from thePolymers to theTime delay in channels of individual aggregates,Representing a sign function,Represents theIndividual Polymer andCommunication connection relationship between polymers, if No.Individual Polymer andIndividual polymers are then1 OtherwiseIs 0;
the expression of the sign function is:
(9);
Where r is the argument of the sign function;
Constraint conditions:
(1) There is only one output channel at the start position:
(10);
(11);
wherein, The representation is from theConnection of the individual polymers to the starting position,From the start position to theConnection of individual polymers to each other,A set of adjacent relay routes representing a starting location;
(2) There is only one input channel at the destination:
(12);
(13);
wherein, The representation is from theConnection of the individual Polymer to the destination position,Representing from the destination location to theConnection relation of individual polymers,A set of adjacent relay routes representing destination locations;
(3) Relay routing has both input and output channels:
(14);
(15);
wherein, The representation is from theConnection relation of individual aggregate to first relay route,The representation is from theConnection relation of individual relay routes to the first aggregate,Represents theA set of adjacent relay routes of the plurality of relay routes;
(4) The number of input channels and the number of output channels of the relay route do not exceed an allowable upper limit:
(16);
(17);
wherein, Represents theUpper limit of number of input channels of each relay route,Represents theThe upper limit of the number of output channels of each relay route.
The optimization model from the formula (8) to the formula (17) can be constructed to perform optimization updating on the communication network, and the adverse effect of longer communication delay on the frequency modulation effect is restrained by adopting a channel with higher transmission speed.
S3, calculating the active power adjustment quantity of the aggregate through a leader-follower consistency algorithm based on the total active power adjustment quantity and the optimized communication network.
As shown in fig. 3, in the present embodiment, for an upper-layer aggregate (fusion terminal), based on a leader-follower consistency algorithm, the equal capacity ratio is realized by means of a sparse communication network between aggregatesRegulating output, wherein the capacity ratio is the ratio of the active power regulating quantity to the rated capacity,Represents the rated capacity of the 1 st polymer,Represents the active power adjustment amount of the 1 st polymer,Represents the rated capacity of the Nth polymer,Indicating the active power adjustment amount of the nth polymer.
Taking the ith polymer as an example, the ratio of the active power adjustment amount to the rated capacity is defined asThe method comprises the following steps: /(I)(18);
The controlled expression of (2) is as follows:
(19);
(20);
wherein, Represents the active power adjustment amount of the ith polymer,Derivative representing the ratio of the active power adjustment quantity to the rated capacity of the ith polymer at time t,Representing the control input of the ith polymer,Represents the ratio of the active power adjustment amount to the rated capacity of the jth polymer at time t,Represents the ratio of the active power adjustment amount to the rated capacity of the ith polymer at time t,Representing the ratio of the active power adjustment amount to the rated capacity of the ith polymer,Representing the ratio of the active power adjustment quantity and the rated capacity of the leader at the time t;
wherein, The expression of (2) is:
(21)。
To ensure that excessive overshoot and steady state errors do not exist in the regulation and control process, a control gain parameter is set AndThe following are provided:
Tracking error expression of the ith polymer The following is shown:
(22);
Then, based on Laplace transformation final value theorem and integral theorem pair Processing is carried out to obtain the following expression:
(23);
therefore, only the proper selection is needed AndSo that equation (24) is satisfied, and then the/>, is found by equation (19) and equation (20)Thus, the active power adjustment quantity of each polymer is obtained, and the active power distribution among the polymers on the upper layer is completed.
(24)。
S4, calculating the active power adjustment quantity of each photovoltaic source in the polymer through a centralized proportional distribution algorithm based on the active power adjustment quantity of the polymer, and performing output regulation and control on the photovoltaic sources through the active power adjustment quantity of the photovoltaic sources to finish frequency modulation.
After each polymer determines the corresponding active power adjustment amount, the active power adjustment amount distribution of each photovoltaic source under the polymer is carried out according to the capacity proportionality coefficient of each photovoltaic source, and for the mth photovoltaic source in the ith polymer, the active power adjustment amount is distributedThe expression is as follows:
(25);
wherein, Represents the active power adjustment quantity of the mth photovoltaic source in the ith polymer,Representing the capacity scaling factor of the mth photovoltaic source in the ith polymer,The active power adjustment amount of the ith polymer at time t is shown.
Example 2
Based on the information physical cooperative frequency modulation method under the multi-distributed photovoltaic polymerization described in embodiment 1, the embodiment of the invention provides an information physical cooperative frequency modulation device under the multi-distributed photovoltaic polymerization, which comprises the following steps:
A data acquisition module configured to: acquiring a total active power adjustment quantity;
A communication network optimization module configured to: optimizing a communication network among the aggregates by improving a consistency algorithm and a preset route optimization scheduling strategy;
An aggregate adjustment quantity distribution module configured to: based on the total active power adjustment amount and the optimized communication network, calculating an active power adjustment amount of the polymer through a leader-follower consistency algorithm;
A photovoltaic source conditioning amount distribution module configured to: and calculating the active power adjustment quantity of each photovoltaic source in the polymer through a centralized proportional distribution algorithm based on the active power adjustment quantity of the polymer, and regulating the output of the photovoltaic sources through the active power adjustment quantity of the photovoltaic sources to finish frequency modulation.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing is merely a preferred embodiment of the present invention, and it should be noted that modifications and variations could be made by those skilled in the art without departing from the technical principles of the present invention, and such modifications and variations should also be regarded as being within the scope of the invention.