The present application claims priority to Chinese Patent Application No. 202111361656.2, entitled “MODEL PREDICTIVE VIRTUAL SYNCHRONOUS GENERATOR INVERTER CONTROL METHOD CONSIDERING FREQUENCY-MOVEMENT DIRECTION”, filed on Nov. 17, 2021 to the China National Intellectual Property Administration, which is incorporated herein by reference in its entirety.
The present invention belongs to the technical field of power electronics, and in particular, relates to a model predictive virtual synchronous generator inverter control method considering a frequency-movement direction.
The description in this section merely provides background information related to the present invention and does not necessarily constitute the prior art.
As an effective utilization form of distributed energy, micro-grid can transmit electric energy and provide grid support, which is an important way for effective utilization of renewable energy and an important supplement to a large power grid. The micro-grid has two modes of grid-connected operation and island operation. In an island mode, the micro-grid loses the support of the large power grid, and thus can support the common bus voltage and frequency only by relying on its internal distributed energy, micro-gas turbine and energy storage. Because the micro-grid with distributed energy lacks rotational kinetic energy of synchronous generators in the large power grid, especially when the proportion of distributed energy increases gradually, the inertia of micro-grid decreases and the stability of micro-grid weakens. A virtual synchronous generator technology can simulate the electromechanical transient characteristics of a synchronous generator, and make an inverter in the island mode have the external characteristics of the synchronous generator, such as inertia, damping, primary frequency modulation, reactive voltage regulation, etc. Therefore, the virtual synchronous generator technology can effectively increase the system inertia, which facilitates improving the stability of micro-grid in the island mode.
The inventors found that the traditional virtual synchronous generator technology has a fixed inertia parameter, which is difficult to adapt to the complex micro-grid environment. When the inertia of virtual synchronous generator is large, it is difficult for the system frequency to regress to the fundamental frequency from which it once deviates, which is disadvantageous for the system power allocation adjustment. On the contrary, when the inertia of the virtual synchronous generator is small, the anti-interference ability of the micro-grid is weak, the system stability is poor, and the load switching can easily lead to system oscillation or instability, and the control effect is poor.
In order to overcome the above-mentioned deficiencies of the prior art, the present invention provides a model predictive virtual synchronous generator inverter control method considering a frequency-movement direction for improving system elasticity. Frequency-movement direction is determined by the absolute value change between the frequency value and the fundamental frequency value.
To achieve the above objective, one or more embodiments of the present invention provide the following technical solutions:
in a first aspect, the present invention discloses the model predictive virtual synchronous generator inverter control method considering a frequency-movement direction, including:
obtaining a current frequency absolute value based on fundamental frequency of an inverter and a steady-state absolute value based on the fundamental frequency of the inverter;
determining a current frequency-movement state by comparing the above two absolute values: in a state of deviating from the fundamental frequency or a state of regressing to the fundamental frequency;
setting a corresponding prediction output horizon according to the frequency-movement state, and constructing a cost function considering the frequency-movement state; and
calculating an optimal virtual power increment value based on the cost function, then calculating an optimal virtual power, inputting the optimal virtual power to a swing equation of a virtual synchronous generator, and performing frequency adjustment according to an output value.
As a further embodiment, the obtaining a current frequency absolute value based on fundamental frequency of an inverter and a steady-state absolute value based on the fundamental frequency of the inverter includes:
detecting a system frequency of output of an inverter in a micro-grid island mode, calculating a system steady-state frequency from a current power of the inverter, and determining the current frequency absolute value based on the fundamental frequency and the steady-state absolute value based on the fundamental frequency.
As a further embodiment, when determining the frequency-movement state, if the current frequency absolute value based on the fundamental frequency is less than the steady-state absolute value based on the fundamental frequency, it is in a state of deviating from the fundamental frequency; and
if the current frequency absolute value based on the fundamental frequency is greater than the steady-state absolute value based on the fundamental frequency, it is in a state of regressing to the fundamental frequency.
As a further embodiment, the cost function considering the frequency-movement state includes a prediction horizon minimization target for a frequency-movement fluctuation value and includes a prediction horizon minimization target for a frequency regression fluctuation value.
As a further embodiment, according to different frequency-movement states, assignments of the deviation state weight factor matrix, the regression state weight factor matrix and the power fluctuation penalty factor matrix are different,
where in a frequency-deviation state, priority is given to realizing the optimization of the prediction horizon output for the frequency-deviation fluctuation value through assignment; and
in a frequency regression state, priority is given to realizing the optimization of the prediction horizon output for the frequency regression fluctuation value through assignment.
As a further embodiment, the optimal virtual power is calculated from the optimal virtual power increment value and the virtual power.
As a further embodiment, after calculating the optimal virtual power, a step of correcting the optimal virtual power is further included, whether the calculated optimal virtual power exceeds a virtual power selectable range is determined, and if so, the current frequency absolute value based on the fundamental frequency of the inverter and the steady-state absolute value based on the fundamental frequency of the inverter are re-obtained, and the optimal virtual power is calculated again;
otherwise, the calculated optimal virtual power is directly input into the swing equation of the virtual synchronous generator, and the calculated virtual angular velocity is output to an integral link to achieve optimal virtual power frequency adjustment.
In a second aspect, the present invention discloses a model predictive virtual synchronous generator inverter control system considering a frequency-movement direction, including:
a frequency-movement state determination module configured to obtain a current frequency absolute value based on fundamental frequency of an inverter and a steady-state absolute value based on the fundamental frequency of the inverter,
determining a frequency-movement state by comparing the above two absolute values: in a state of deviating from the fundamental frequency or a state of regressing to the fundamental frequency;
a frequency-adjustment module configured to set a corresponding prediction output horizon according to the frequency-movement state, and construct a cost function considering the frequency-movement state; and
calculate an optimal virtual power increment value based on the cost function, then calculate an optimal virtual power, input the optimal virtual power to a swing equation of a virtual synchronous generator, and perform frequency adjustment according to an output value.
The above one or more technical solutions have the following beneficial effects:
The advantages of additional aspects of the present invention will be given in part in the description which follows and, in part, will be obvious from the following description, or may be learned by practice of the present invention.
The accompanying drawings constituting a part of the present invention are used for providing a further understanding of the present invention. The exemplary examples of the present invention and descriptions thereof are used for explaining the present invention, and do not constitute an improper limitation of the present invention.
It should be noted that, the following detailed descriptions are all exemplary, and are intended to provide further descriptions for the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those usually understood by a person of ordinary skill in the art to which the present invention belongs.
It should be noted that the terms used herein are merely used for describing specific implementations, and are not intended to limit exemplary implementations of the present invention.
The embodiments of the present invention and the features of the embodiments may be combined with each other without conflict.
As described in the background, according to the existing virtual synchronous generator technology for inverters, it is difficult for the system frequency to regress to the fundamental frequency from which it once deviates under large inertia conditions. Under the condition of small inertia, the micro-grid system has weak anti-interference ability and poor system stability, which leads to poor control effect. In order to solve this problem, this embodiment provides a model predictive virtual synchronous generator inverter control method considering a frequency-movement direction.
Specifically, as shown in
Step 1: a system frequency ωm of output of an inverter in a micro-grid island mode is detected, a system steady-state frequency ωref is calculated from the current output power of the inverter, and the frequency absolute value |ωm0(k)| based on the fundamental frequency is an absolute value of the difference between the system frequency ωm and the fundamental frequency ω0; and the steady-state frequency absolute value |ωref0(k)| based on the fundamental frequency is an absolute value of the difference between the steady-state frequency ωref and the fundamental frequency ω0.
Specifically, the frequency absolute value |ωm0(k)| at the kth instant based on the fundamental frequency ω0 and the steady-state absolute value |ωref0(k)| based on the fundamental frequency ω0 are determined, with the calculation formula thereof as follows:
In this step, the system frequency ωm is an output frequency output by the inverter, and the steady-state frequency ωref is derived from the curve of droop control P−ω of the inverter.
Step 2: in order to determine the system frequency-movement state, an absolute value |ωm0(k)| of the current frequency value ωm based on the fundamental frequency is compared with an absolute value |ωref0(k)| of the steady-state frequency value ωref based on the fundamental frequency, with the determination formula thereof as follows:
In order to realize the frequency control considering the frequency-movement state, it is necessary to determine the current system frequency state.
when |ωm0(k)|−|ωref0(k)|<0, it is determined that the system frequency ωm is in a state of deviating from the fundamental frequency; and
when |ωm0(k)|−|ωref0(k)|>0, it is determined that the system frequency ωm is in a state of regressing to the fundamental frequency.
Step 3: a corresponding prediction output horizon is set according to the frequency-movement state, and a cost function considering the frequency-movement state is constructed, including a prediction horizon minimization target for a frequency-movement fluctuation value and a prediction horizon minimization target for a frequency regression fluctuation value.
When the system frequency om is in a state of deviating from the fundamental frequency, a prediction horizon output Wpd (j+1|j) of the frequency-movement fluctuation value is constructed, i.e., the value of j+1 instant is predicted from the j instant (j=k, k+1 . . . ), and an expected value R(k+1)=[0 0 . . . 0]T in the deviation state (an expected value at k+1 instant) is set; and when the system frequency ωm is in a state of regressing to the fundamental frequency, the prediction horizon output Wpb(j+1|j) of the frequency regression fluctuation value is constructed, and an expected value Rb(k+1)=[0 0 . . . 0]T in the regression state is set.
Further, by optimizing the virtual power increment value in an incremental model, the system frequency ωm can be adjusted effectively in different states.
In order to adapt to different optimization objectives in a frequency-deviation state and a frequency regression state, the cost function Fp considering the frequency-movement state is constructed.
The derivation process is described in detail below:
The incremental model is derived from a discrete model of a swing equation of the virtual synchronous generator technology, and its calculation formula is as follows:
where
is a control period, J is a moment of inertia, Pm is virtual power, Pout is output power and D is a damping coefficient.
ωd(j)=ωm(j)−ωm(k)(j=k,k+1, . . . ,p)
W
pd(j+1|j)=SxΔωd(j)+Iωd(j)+SdΔPout(j)+SuΔP(j)
ωb(j)=ωm(j)−ωrefj=k,k+1, . . . ,p
Where ωref is a steady-state frequency expected value.
Further, the prediction horizon output Wpb(j+1|j) for the frequency regression fluctuation value ωb is
W
pb(j+1|j)=SxΔωb(j)+Iωb(j)+SdΔPout(j)+SuΔP(j)
Therefore, when the frequency deviates from the fundamental frequency, in order to achieve slow frequency movement, the prediction horizon output Wpd(j+1|j) of the frequency-deviation fluctuation value ωd should be 0, i.e., Rd(k+1)=[0 0 . . . 0]T; and when the frequency regresses from the deviation value to the fundamental frequency, in order to achieve frequency-accelerated regression, the prediction horizon output Wpb(j+1|j) of the frequency regression fluctuation value ωb is 0, i.e., Rb(k+1)=[0 0 . . . 0]T.
Further, the prediction horizon output Wpd(j+1|j) for the frequency-deviation fluctuation value ωd and the prediction horizon output Wpb(j+1|j) for the frequency regression fluctuation value wb are integrated, and the cost function considering the frequency-movement state is constructed as follows:
F
p=min{∥ΓdWpd(k+1|k)−Rd(k+1)∥2+∥ΓbWpd(k+1|k)−Rb(k+1)∥2+∥ΓpΔP(k)∥2}
where Γd is a deviation state weight factor matrix, Γb is a regression state weight factor matrix, and Γ is a power fluctuation penalty factor matrix.
Step 4; in order to achieve the frequency optimization control targets under different movement states, the assignments of matrixes Γd, Γb and Γ are different according to different frequency-movement states.
When in the state of the system frequency deviating from the fundamental frequency, a weighting factor and a penalty factor are assigned to optimize the frequency change rate to be the minimum, and an optimal virtual power increment value under the frequency-movement state is calculated to slow down the deviation of the system frequency ωm and effectively reduce the frequency change rate.
Specifically, in the frequency-deviation state, priority is given to realizing the optimization of the prediction horizon output Wpd(j+1|j) for the frequency-deviation fluctuation value ωd, and its assignment is as follows:
frequency-deviation state:
When in the state of the frequency regressing to the fundamental frequency state, the weight factor and the penalty factor are assigned to optimize the difference between a steady-state value and a frequency to be the minimum, and the optimal virtual power increment value is calculated.
Specifically, when in the state of the system frequency regressing to the fundamental frequency state, the weighting factor and the penalty factor are assigned to calculate the optimal virtual power increment value in the frequency regression state to accelerate the system frequency ωm to regress to the fundamental frequency, and effectively and rapidly increase the frequency change rate, priority is given to realizing the dominance of optimization of the prediction horizon output Wpd(j+1|j) for the frequency regression fluctuation value ωb, and the assignment thereof is as follows:
According to the extreme value theory, the optimal virtual power increment ΔP*(k) under different frequency-movement states is calculated.
Step 5: the optimal virtual power is calculated according to the optimal virtual power increment value and the virtual power, and the expression thereof is
P*(k)=ΔP*(k)+P(k)
In order to ensure the reliability of the system, it needs to determine whether the optimal virtual power P*(k) exceeds an adjustable range thereof, and when correcting the correctness of the virtual power, whether the optimal virtual power exceeds a selectable range of the virtual power, and if so, step 1 is performed; otherwise, the calculated virtual power is input into a swing equation of a virtual synchronous generator, i.e.,
the calculated virtual angular velocity ωm is input to an integral link to achieve optimal virtual power frequency adjustment.
According to the present invention, the system frequency is ultimately adjusted by injecting the optimal virtual power into the virtual synchronous power generation swing equation considering the frequency-movement state.
It should be noted that, the technical solutions of the present invention include virtual synchronous generator increment model establishment, system frequency-movement state determination, cost function considering the frequency-movement state, optimization of virtual power increment, and virtual power adjustment of frequency. The present invention determines a system frequency-movement state by using an absolute value of a frequency-movement based on the fundamental frequency, and further sets a frequency optimization target in different frequency-movement states.
The present invention considers the frequency-movement direction, and is based on a model predictive method of the virtual synchronous generator, realizing realizes the adjustment of the system frequency by optimizing the virtual power in different movement states. When the frequency deviates from the fundamental frequency, the frequency changes slowly to suppress the frequency-deviation; and when the frequency regresses to the fundamental frequency, the frequency change is accelerated to cause the frequency to quickly regress to the fundamental frequency. In addition, the characteristics of real-time optimization of model prediction improve the control effect. The present invention can solve the effect on the frequency caused micro-grid load switching, and the implementation method is simple and reliable. Therefore, it is of great significance for the application of the inverter.
The present invention is a real-time optimization process to achieve the optimal regulation of the system frequency in a frequency-deviation state and a frequency regression state to improve the performance of the inverter.
In addition, the present invention is not limited to specific micro-grid practical requirements, is not limited to inverter topologies, and is applicable to micro-grid single-phase bus and three-phase bus forms. In addition, the present invention is not limited to a DC side power supply form of the inverter, and is applicable to different situations such as low voltage, medium voltage, and high voltage, and has strong scalability and practicality.
Simulation Case
It can be seen from the above simulation results that the model predictive virtual synchronous generator inverter control method considering a frequency-movement direction provided by the present invention can effectively improve the dynamic characteristic capability of the inverter to regulate frequency.
In another embodiment, a virtual synchronous generator inverter is further included, which is controlled by using the above-mentioned method to achieve different frequency adjustment targets, i.e.: when the frequency deviates from the fundamental frequency, the frequency change rate is effectively reduced to achieve slow frequency-movement; and when the frequency regresses to the fundamental frequency from the non-fundamental frequency value, the frequency change rate is increased effectively and quickly, so that the frequency regresses quickly.
An objective of this embodiment is to provide a computing device including a memory, a processor and a computer program stored on the memory and executable on the processor, and when the processor executes the program, the steps of the method described above are implemented.
An objective of this embodiment is to provide a computer-readable storage medium.
The computer-readable storage medium has a computer program stored thereon, and the program when executed by a processor, performs the steps of the above method.
An objective of this embodiment is to provide a model predictive virtual synchronous generator inverter control system considering a frequency-movement direction, which includes:
a frequency-movement state determination module configured to obtain a current frequency absolute value based on fundamental frequency of an inverter and a steady-state absolute value based on the fundamental frequency of the inverter;
determining a frequency-movement state by comparing the above two absolute values: in a state of deviating from the fundamental frequency or a state of regressing to the fundamental frequency;
a frequency-adjustment module configured to set a corresponding prediction output horizon according to the frequency-movement state, and construct a cost function considering the frequency-movement state; and
calculate an optimal virtual power increment value based on the cost function, then calculate an optimal virtual power, input the optimal virtual power to a swing equation of a virtual synchronous generator, and perform frequency adjustment according to an output value.
The steps involved in the devices of Embodiments II, III and IV above correspond to those of method Embodiment I, and for specific implementations, reference can be made to the relevant descriptions of Embodiment I. The term “computer-readable storage medium” shall mean a single medium or multiple media including one or more sets of instructions; It should also be understood to include any medium capable of storing, encoding, or carrying a set of instructions for execution by the processor and causing the processor to perform any of the methods of the present invention.
It will be appreciated by a person skilled in the art that the various modules or steps of the present invention described above may be implemented in a general-purpose computing device. Optionally, they may be implemented in program code executable by the computing device, such that they may be stored in a storage device and executed by the computing device, or they may be fabricated as individual integrated circuit modules separately, or multiple modules or steps thereof may be fabricated as a single integrated circuit module. The present invention is not limited to any specific combination of hardware and software.
The specific implementations of the present invention are described above with reference to the accompanying drawings, but are not intended to limit the protection scope of the present invention. A person skilled in the art should understand that various modifications or deformations may be made without creative efforts based on the technical solutions of the present invention, and such modifications or deformations shall fall within the protection scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
2021113616562 | Nov 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2022/130367 | 11/7/2022 | WO |