CN113572199B - Smooth switching method of network-forming type current converter based on model prediction algorithm - Google Patents
Smooth switching method of network-forming type current converter based on model prediction algorithm Download PDFInfo
- Publication number
- CN113572199B CN113572199B CN202110917689.4A CN202110917689A CN113572199B CN 113572199 B CN113572199 B CN 113572199B CN 202110917689 A CN202110917689 A CN 202110917689A CN 113572199 B CN113572199 B CN 113572199B
- Authority
- CN
- China
- Prior art keywords
- time
- moment
- angular frequency
- converter
- grid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/041—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a variable is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/388—Islanding, i.e. disconnection of local power supply from the network
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/40—Synchronising a generator for connection to a network or to another generator
- H02J3/44—Synchronising a generator for connection to a network or to another generator with means for ensuring correct phase sequence
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a smooth switching method of a network-building type converter based on a model prediction algorithm, which comprises the following steps of 100, discretizing an active frequency model of a VSG (voltage-to-grid generator) system of a simulation synchronous generator to obtain a state space model; 200, describing the state space model in an incremental form, wherein the incremental form model predicts the state of the system in one or more future periods through the state space models at the current time and the last time; step 300, setting a target function for predicting smooth switching control; step 400, obtaining an optimal solution of the objective function under constraint conditions, comparing the optimal solution with a power increment required by the previous moment at the current moment, and adding the optimal solution with the active power input value at the current moment to obtain the active power input value at the current moment; step 500, when the alternating current micro grid system is switched from an island operation mode to a grid-connected operation mode, the predictive smooth switching control is completed through the steps 100 to 400, and the obtained optimal input is input into the alternating current micro grid system.
Description
Technical Field
The patent belongs to the technical field of control, and particularly relates to a smooth switching method of a network-forming type converter based on a model prediction algorithm.
Background
Distributed Generation (DG) is operated in the form of an intelligent alternating current micro-grid and becomes a new development direction of a future power grid, but most of distributed generation is connected into the power grid through a power electronic converter with rapid response, inertia and damping are lacked, inertia of the whole system can be reduced after the distributed generation is connected into a large power grid, and when disturbance occurs in the system, frequency is easy to fluctuate and even exceeds a safe operation range, so that stable operation and dynamic response of the whole power system are influenced. The virtual synchronous generator technology is characterized in that a converter control algorithm is designed, so that the converter simulates the frequency and voltage control characteristics of a Synchronous Generator (SG) from the external characteristics, and the dynamic regulation performance of a power grid is improved.
In addition, the alternating-current micro-grid can be connected to a large power grid for grid-connected operation and can be separated from the large power grid for off-grid operation according to actual requirements, but in the operation mode conversion process of the alternating-current micro-grid, the problem that the output power of a micro-power source is not matched with a load is easily caused in the alternating-current micro-grid, so that transient current or voltage impact is caused, and therefore smooth switching control needs to be carried out on the alternating-current micro-grid. In the traditional control, PI controllers are usually adopted to realize the presynchronization of frequency and phase angle, the control method is really classic, but the presynchronization of frequency and phase angle is mediated by two PI controllers, and the frequency and the phase angle are correlated with each other, so that large fluctuation is easy to occur in the actual control due to the lack of unified control, and the operation of a power system is damaged. In addition, the PID control has the disadvantages of slow response speed, low control accuracy for a nonlinear system, and the like.
Disclosure of Invention
The present invention is proposed based on the above requirements in the prior art, and the technical problem to be solved by the present invention is to provide a method for smoothly switching a network-forming converter based on a model prediction algorithm, so as to perform dynamic response quickly and efficiently.
In order to solve the above problems, the technical solution provided by this patent includes:
the utility model provides a smooth switching method of a network-forming type converter based on a model prediction algorithm, which comprises the following steps: step 100, discretizing an active frequency model of a simulated synchronous generator VSG system to obtain a state space model of the active frequency of the discrete-time VSG system, wherein the state space model obtains a difference value between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the current moment through the difference value between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the previous moment, the mechanical power of a converter and the electromagnetic power output by the converter; step 200, describing the state space model in an incremental mode, wherein the incremental mode model obtains the difference value between the actual angular frequency and the rated angular frequency of the alternating current micro-grid at the later moment compared with the current moment through the difference value between the actual angular frequency and the rated angular frequency of the alternating current micro-grid at the current moment and the last moment, the mechanical power difference value of the converter and the electromagnetic power difference output by the converter, namely the required power increment, and predicting the state of the system in one or more periods in the future through the mode; step 300, setting an objective function of a predictive smooth handover control, said predictive smooth handover controlThe method is characterized in that the frequency and the phase angle of an alternating current micro-grid are controlled to be equal to those of a large power grid, performance indexes of the method comprise the phase angle and the angular speed, and introduced angular acceleration variation is adaptively adjusted through a weighting coefficient; step 400, obtaining an optimal solution of the objective function under constraint conditions, wherein the constraint conditions comprise a state space model in an increment form at the k + i moment, and the optimal solution set of the objective function is obtained as delta P M ＝[ΔP m (k) ΔP m (k+1) ... ΔP m (k+p-1)] T In which Δ P m (k) Representing the difference, Δ P, between the mechanical power input to the converter at time k and at time k-1 m (k + 1) represents the difference between the mechanical power input to the converter at the time k +1 and at the time k, Δ P m (k + p-1) represents a mechanical power difference value of the converter input at the k + p-1 moment and the k + p-2 moment, a first item in an optimal solution set is an optimal solution which is a power increment required by the current moment compared with the previous moment and is added with an active power input value of the current moment to obtain an active power input value of the current moment, and p is a prediction time domain; step 500, when the alternating current micro grid system is switched from an island operation mode to a grid-connected operation mode, the predictive smooth switching control is completed through the steps 100 to 400, and the obtained optimal input is input into the alternating current micro grid system.
Preferably, the model simulating the active frequency of the synchronous generator VSG can be expressed as:
wherein J is moment of inertia, D P To be damping coefficient, P m 、P e The mechanical power and the electromagnetic power of the converter are respectively, the given value and the actual output value of the active power are respectively, omega represents the virtual angular speed output by the converter, and omega represents the virtual angular speed output by the converter 0 Is the nominal angular velocity.
Preferably, the discretized active frequency model expression of the analog synchronous generator VSG system is as follows:
wherein, ω' (k + 1) represents the sum of the actual angular frequencies of the AC micro-grid at the moment k +1A difference in nominal angular frequency; omega' (k) represents the difference between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the moment k; p m (k) The mechanical power input into the converter at the moment k is shown; p e (k) Representing the electromagnetic power output by the converter at the moment k; y is c (k) Representing the controlled output variable at time k, T s Is the sampling time of the system, and
preferably, the state space model is described in incremental form:
wherein, Δ ω' (k + 1) represents the difference between the actual angular frequency of the ac microgrid at the time k +1 and the actual angular frequency of the ac microgrid at the time k; Δ ω' (k) represents the difference between the actual angular frequency of the AC microgrid at time k and at time k-1, Δ P m (k) Representing the difference, Δ P, between the mechanical power input to the converter at time k and at time k-1 e (k) Representing the difference in electromagnetic power output by the converter at time k and at time k-1, y c (k) Representing the controlled output variable at time k, y c (k-1) represents the controlled output variable at time k-1, C c Representing the output state coefficient matrix.
Preferably, the states in p prediction time domains in the future of the prediction system can be described as follows:
Y p (k+1|k)＝S x Δω′(k)+S u ΔP M (k)+S d ΔP e (k)+γω′(k)
wherein Y is p (k +1 purple k) as the predicted future pThe system state of the step is output as a variable,
S x ＝[A A 2 +A ... A p +A p-1 +...+A] T
ΔP M ＝[ΔP m (k) ΔP m (k+1) ... ΔP m (k+p-1)] T
S d ＝[B d AB d +B d ... A p-1 B d +A p-2 B d +...+B d ] T
γ＝[1 1 ... 1] T 。
preferably, the objective function includes two parts of phase angle and angular velocity, and is represented as:
wherein, gamma is yi Is the angular frequency weight coefficient, θ g (k + i) is the large grid phase angle at time k + i, ω g ' (k + i) is the large grid angular frequency at time k + i, [ theta ] (k + i | k) is the phase angle at time k + i VSG, [ omega ] (k + i | k) is the angular frequency at time k + i VSG,
wherein gamma is 0 As initial value of the weighting coefficient, k Γ The weighting adjustment coefficients.
Preferably, the optimization problem of model predictive control is expressed as follows by using constraints including an incremental state space model and the like:
wherein J (ω '(k)) is an objective function, ω' min (k + i) is the minimum value of VSG angular frequency at time k + i, ω' max (k + i) is the maximum value of VSG angular frequency at the time k + i, and Δ P is taken m (k) Δ P of the first term m (k) As the most importantThe preferred input power increment.
Preferably, the optimal input of the current system is the power increment required by the current time compared with the previous time, and then added to the active power input value of the previous time, that is: p m (k)＝P m (k-1)+ΔP m (k)。
Preferably, the obtained optimal input of the current system is input into an alternating current micro-grid system, and the alternating current micro-grid system comprises a simulation synchronous generator VSG unit, an active load and a synchronous generator set SG of the distributed generation DG.
Compared with the prior art, the invention has the advantages of high dynamic response speed, strong control purpose, high precision, flexible nonlinear constraint condition processing and easy satisfaction of the requirement of high-performance processing.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the embodiments of the present specification, and other drawings can be obtained by those skilled in the art according to these drawings.
FIG. 1 is a flow chart of steps of a network-forming type converter smooth switching method based on a model prediction algorithm;
FIG. 2 is a diagram of an AC microgrid system in accordance with an embodiment of the present invention;
fig. 3 is a block diagram of active-frequency control of model predictive control in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
For the purpose of facilitating understanding of the embodiments of the present application, the following description will be made in terms of specific embodiments with reference to the accompanying drawings, which are not intended to limit the embodiments of the present application.
Example 1
The embodiment provides a smooth switching method of a network-forming type converter based on a model prediction algorithm.
The method is a method for uniformly controlling frequency presynchronization and phase angle presynchronization by utilizing a model prediction algorithm aiming at an alternating current micro-grid system consisting of VSG units and active loads and a large grid system consisting of a generator set when the alternating current micro-grid system is converted from an island mode to a grid-connected mode, and refer to fig. 1-3.
The active frequency rotor equation of the analog synchronous generator VSG is:
j is a rotational inertia, and provides an inertial support when the system frequency oscillates; d P Simulating the damping oscillation capacity of the synchronous generator for a damping coefficient; p m 、P e The power is respectively the mechanical power and the electromagnetic power of the converter, and is also respectively the given value and the actual output value of the active power; ω represents the virtual angular velocity of the converter output, ω 0 Is the nominal angular velocity.
Inertia of the synchronous generator SG is simulated through the rotational inertia J, although frequency fluctuation of a system after disturbance can be relieved to a certain extent, the fundamental reason is power unbalance, and therefore power shortage needs to be rapidly and accurately complemented by utilizing model prediction control to achieve a control target.
Discretizing the active frequency rotor equation to obtain a state space model of the active frequency of the discrete time system, wherein the state space model is expressed as follows:
wherein ω' (k + 1) represents the difference between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the moment of k + 1; omega' (k) represents the difference value between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the moment k; p is m (k) Representing the mechanical power input into the converter at the moment k; p e (k) Representing the electromagnetic power output by the converter at the moment k; y is c (k) Representing the controlled output variable at time k.
And is
To introduce integration to reduce or eliminate static errors, the state space model is changed to an incremental form:
wherein, Δ ω' (k + 1) represents the difference between the actual angular frequency of the ac microgrid at the time k +1 and the actual angular frequency of the ac microgrid at the time k; delta omega' (k) represents the difference value of the actual angular frequency of the alternating-current microgrid at the moment k and the moment k-1 and represents the increment of the state variable; delta P m (k) The difference value of mechanical power input into the converter at the moment k and the moment k-1 is represented, and represents the increment of a control input variable; delta P e (k) Representing the difference value of the electromagnetic power output by the converter at the moment k and the moment k-1, and representing the measurable external interference variable increment; y is c (k) Representing a controlled output variable at the time k; y is c (k-1) represents the controlled output variable at the time of k-1, C c Representing the output state coefficient matrix.
And predicting the state of the system in p periods in the future by using the incremental state space model. The p is a prediction time domain.
In this embodiment, if p =3, the output equation for the future 3-step prediction of the system is:
Y 3 (k+1|k)＝S x Δω′(k)+S u ΔP M (k)+S d ΔP e (k)+γω′(k)
wherein, Y 3 (k +1 dichotomy k) is the system state output variable for predicting 3 steps in the future.
S x ＝[A A 2 +A A 3 +A 2 +A] T
ΔP M ＝[ΔP m (k) ΔP m (k+1) ΔP m (k+2)] T
S d ＝[B d AB d +B d A 2 B d +AB d +B d ] T
γ＝[1 1 1] T
And step 300, setting an objective function of the predictive smooth switching control, wherein the predictive smooth switching control aims to control the frequency and the phase angle of the alternating current micro-grid to be equal to those of the large grid, and performance indexes of the predictive smooth switching control comprise the phase angle and the angular speed.
The target of the prediction smooth switching control is to control the frequency and the phase angle of the alternating current micro-grid to be equal to those of the large grid so as to achieve the target of pre-synchronization control.
In the presynchronization process, when the phase angle difference is 0, the frequency difference is also 0, so the main performance index in the objective function is the phase angle difference.
In a very short period T s The relationship between voltage phase angle and angular frequency can be expressed as:
θ (k + 1) represents the voltage phase angle at time k + 1, θ (k) represents the voltage phase angle at time k, ω (k + 1) represents the angular frequency at time k + 1, ω 0 Is the nominal angular velocity.
From the above relational expression of the voltage phase angle and the angular frequency, if only the phase angle difference is considered to be reduced to 0, the alternating-current microgrid phase angle theta can rapidly approach the large power grid angular frequency theta g But at the same time, the angular frequency ω will also fluctuate greatly, even exceeding the safe operation range, which is not favorable for the stable operation of the system, so that the angular frequency difference needs to be added as a performance index.
The objective function derived from the above analysis is:
in the formula of gamma yi The angular frequency weight coefficient is used for representing the importance degree of each target item in the performance index. It can be seen that the main purpose of the objective function is to pass the large grid angular frequency ω at times k +1 to k +3 g ' sum phase angle θ g Predicted value of (2), 1The optimal angular frequency dynamic response is obtained by considering the back-stepping calculation, the purposes of phase angle pre-synchronization and frequency pre-synchronization can be achieved, the system frequency can be maintained in a normal operation range in the process, the optimal input power of the VSG at the current k moment is obtained by back stepping, and then the optimal control effect is achieved by inputting the optimal input power into the system. Theta g (k + i) is the large grid phase angle at time k + i, ω g 'is the large grid angular frequency at time k + i, θ (k + i | k) is the VSG phase angle at time k + i, and ω' (k + i | k) is the VSG angular frequency at time k + i.
The objective function comprises two parts of a phase angle theta and an angular speed omega, and not only is the phase difference ensured to be fast and accurately reduced to 0, but also the angular speed is prevented from changing too much to influence the stability of the system frequency in the presynchronization process.
When the weighting coefficient is smaller, the phase difference is a main control target, the phase difference is reduced rapidly, but meanwhile, the angular speed is equivalent to the frequency, and larger-range fluctuation is generated, so that the normal and stable operation of the system is influenced; when the weighting coefficient is large, the angular speed is equivalent to the frequency as a main control target, the frequency fluctuation range is small at the moment, the system is stable, the time for reducing the phase difference to 0 is greatly prolonged, and the pre-synchronization time is long, so that the quick recovery of the grid-connected state of the alternating-current micro-grid is not facilitated.
Therefore, the amount of angular velocity change can be introduced into the weighting coefficient for adaptive adjustment, and the specific equation is as follows:
Γ 0 as initial value of the weighting coefficient, k Γ The weighting adjustment coefficients. When the angular velocity difference is large, the weighting coefficient is large, which is beneficial to the stability of frequency, and when the angular velocity difference is small, the weighting coefficient is small, which is beneficial to the rapid implementation of phase pre-synchronization.
According to the invention, the phase angle difference and the frequency difference between the micro-grid and the large-grid are taken as performance indexes of a target function, and the weights of the two parts are timely adjusted through variable weighting coefficients, so that the stable operation and the good dynamic response characteristic of the micro-grid are ensured while the presynchronous control is realized.
Using constraints including an incremental form state space model, and the like, an optimization problem of model predictive control is expressed as:
wherein J (ω '(k)) is an objective function, ω' min (k + i) is the minimum value of VSG angular frequency at time k + i, ω' max (k + i) is the maximum value of the VSG angular frequency at the time k + i,
solving the optimal solution of the model predictive control MPC problem with the constraint, and taking the delta P m (k) Δ P of the first term m (k) As an input power increment.
Obtaining the optimal input of the current system:
P m (k)＝P m (k-1)+ΔP m (k)。
The alternating-current micro-grid system in the embodiment comprises a VSG unit, an active load and a SG of a distributed power generation DG.
Under the initial operation state, the VSG unit is in an island operation mode and operates independently under load, and secondary frequency modulation control is adopted at the moment to ensure the frequency stability of the alternating current micro-grid; the synchronous generator set simulates the operation of a large power grid and is disconnected with an alternating current micro-power grid.
When the alternating-current microgrid is switched from an island operation mode to a grid-connected operation mode, the VSG unit is switched from the secondary frequency modulation control 0 to the model predictive control MPC, namely, the steps 100 to 400 are carried out, so that the input power instruction of the VSG unit is adjusted, pre-synchronization operation is carried out, and after the pre-synchronization is completed, grid-connected switch-on is carried out, and grid-connected operation is completed.
The control of the MPC algorithm directly influences the pre-synchronization effect and the frequency dynamic characteristics of the power grid in the pre-synchronization process.
The above-mentioned embodiments, objects, technical solutions and advantages of the present application are described in further detail, it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present application, and are not intended to limit the scope of the present application, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present application should be included in the scope of the present application.
Claims (8)
1. A smooth switching method of a network-forming type converter based on a model prediction algorithm is characterized by comprising the following steps:
step 100, discretizing an active frequency model of a simulation synchronous generator VSG system to obtain a state space model of the active frequency of the discrete-time VSG system, wherein the state space model obtains a difference value between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the current moment through the difference value between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the previous moment, the mechanical power of a converter and the electromagnetic power output by the converter;
step 200, describing the state space model in an incremental mode, wherein the incremental mode model obtains the difference value between the actual angular frequency and the rated angular frequency of the alternating current micro-grid at the later moment compared with the current moment through the difference value between the actual angular frequency and the rated angular frequency of the alternating current micro-grid at the current moment and the last moment, the mechanical power difference value of the converter and the electromagnetic power difference output by the converter, namely the required power increment, and predicting the state of the system in one or more periods in the future through the mode;
step 300, setting a target function of predictive smooth switching control, wherein the predictive smooth switching control takes the frequency and the phase angle of an alternating current micro-grid as the target and controls the frequency and the phase angle to be equal to those of a large grid, performance indexes of the predictive smooth switching control comprise the phase angle and the angular speed, and the introduced angular acceleration variation is subjected to adaptive adjustment through a weighting coefficient; the objective function includes two parts, phase angle and angular velocity, and is expressed as:
wherein, gamma is yi Is the angular frequency weight coefficient, θ g (k + i) is the large grid phase angle at time k + i, ω g ' (k + i) is the large grid angular frequency at time k + i, [ theta ] (k + i | k) is the VSG phase angle at time k + i, [ omega ] (k + i | k) is the VSG angular frequency at time k + i,
wherein gamma is 0 As initial value of the weighting coefficient, k Γ Is a weighting adjustment coefficient;
step 400, obtaining the optimal solution of the target function under the constraint condition, wherein the constraint condition comprises a state space model of an increment form at the moment of k + i, and the optimal solution set of the target function is obtained as delta P M ＝[ΔP m (k)ΔP m (k + 1)...ΔP m (k+p-1)] T In which Δ P m (k) Representing the difference, Δ P, between the mechanical power input to the converter at time k and at time k-1 m (k + 1) represents the difference between the mechanical power input to the converter at the time k +1 and at the time k, Δ P m (k + p-1) represents a mechanical power difference value of the converter input at the k + p-1 moment and the k + p-2 moment, a first item in the optimal solution set is an optimal solution which is a power increment required by the current moment compared with the previous moment and is added with an active power input value of the current moment to obtain an active power input value of the current moment, and p is a prediction time domain;
step 500, when the alternating current micro grid system is switched from an island operation mode to a grid-connected operation mode, the predictive smooth switching control is completed through the steps 100 to 400, and the obtained optimal input is input into the alternating current micro grid system.
2. The method for smoothly switching the network-forming type converter based on the model prediction algorithm is characterized in that the model for simulating the active frequency of the VSG can be expressed as follows:
wherein J is moment of inertia, D P To be damping coefficient, P m 、P e Respectively mechanical power and electromagnetic power of the converter, and respectively active power set value and actual output value, wherein omega represents virtual angular velocity output by the converter, and omega represents virtual angular velocity output by the converter 0 Is the nominal angular velocity.
3. The method for smoothly switching the network-forming type converter based on the model prediction algorithm as claimed in claim 2, wherein the discretized active frequency model expression of the analog synchronous generator VSG system is as follows:
wherein ω' (k + 1) represents the difference between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the moment of k + 1; omega' (k) represents the difference between the actual angular frequency and the rated angular frequency of the alternating-current microgrid at the moment k; p m (k) Representing the mechanical power input into the converter at the moment k; p e (k) Representing the electromagnetic power output by the converter at the moment k; y is c (k) Representing the controlled output variable at time k, T s Is the sampling time of the system, and
4. the method according to claim 3, wherein the state space model is described in an incremental form by:
wherein, Δ ω' (k + 1) represents the difference between the actual angular frequency of the ac microgrid at the time k +1 and the actual angular frequency of the ac microgrid at the time k; Δ ω' (k) represents the difference between the actual angular frequency of the AC microgrid at time k and at time k-1, Δ P m (k) Representing the difference, Δ P, between the mechanical power input to the converter at time k and at time k-1 e (k) Representing the difference in electromagnetic power output by the converter at time k and at time k-1, y c (k) Representing the controlled output variable at time k, y c (k-1) represents the controlled output variable at time k-1, C c Representing the output state coefficient matrix.
5. The method for smoothly switching the network-structured converter based on the model prediction algorithm according to claim 4, wherein the states in p future prediction time domains of the prediction system can be described as follows:
Y p (k+1|k)＝S x Δω′(k)+S u ΔP M (k)+S d ΔP e (k)+γω′(k)
wherein, Y p (k +1 k) is a system state output variable for predicting p steps in the future,
S x ＝[A A 2 +A...A p +A p-1 +...+A] T
ΔP M ＝[ΔP m (k) ΔP m (k+1)...ΔP m (k+p-1)] T
S d ＝[B d AB d +B d ...A p-1 B d +A p-2 B d +...+B d ] T
γ＝[1 1...1] T 。
6. the method for smoothly switching the grid-structured converter based on the model prediction algorithm according to claim 5, wherein the constraint condition comprising the state space model in an incremental form is used to represent the optimization problem of model prediction control as follows:
wherein, J m (ω '(k)) is an objective function, ω' min (k + i) is the minimum value of VSG angular frequency at time k + i, ω' max (k + i) is the maximum value of VSG angular frequency at the time of k + i, and is taken as delta P m (k) The first term in (1) is taken as the optimum input power increment.
7. The method according to claim 6, wherein the optimal input of the current system is the power increment required by the current time compared with the previous time, and the power increment is added to the active power input value of the previous time, that is:
P m (k)＝P m (k-1)+ΔP m (k)。
8. the method for smoothly switching the grid-structured converter based on the model prediction algorithm according to claim 7, wherein the obtained optimal input of the current system is input into an alternating current micro-grid system, and the alternating current micro-grid system comprises a VSG (synchronous generator) unit, an active load and a SG (synchronous generator set) of a distributed generation DG.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110917689.4A CN113572199B (en) | 2021-08-11 | 2021-08-11 | Smooth switching method of network-forming type current converter based on model prediction algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110917689.4A CN113572199B (en) | 2021-08-11 | 2021-08-11 | Smooth switching method of network-forming type current converter based on model prediction algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113572199A CN113572199A (en) | 2021-10-29 |
| CN113572199B true CN113572199B (en) | 2022-12-27 |
Family
ID=78171194
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110917689.4A Active CN113572199B (en) | 2021-08-11 | 2021-08-11 | Smooth switching method of network-forming type current converter based on model prediction algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113572199B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114123243B (en) * | 2021-11-17 | 2023-10-27 | 山东大学 | Model prediction virtual synchronous machine inverter control method considering frequency offset direction |
| CN115036967B (en) * | 2022-06-09 | 2023-06-06 | 龙源（北京）风电工程设计咨询有限公司 | Three-phase converter switching method and system, readable medium and electronic equipment |
| CN116845886B (en) * | 2023-08-29 | 2023-12-05 | 华能江苏综合能源服务有限公司 | Multi-port autonomous photovoltaic system network construction control method based on model prediction |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109950922A (en) * | 2019-01-31 | 2019-06-28 | 东南大学 | A kind of multistep model predictive control method suitable for VSC-HVDC |
| WO2021035269A1 (en) * | 2019-08-30 | 2021-03-04 | Avl List Gmbh | Method and controller for model predictive control of coupled converters |
| CN113179059A (en) * | 2021-05-21 | 2021-07-27 | 南京理工大学 | Improved virtual synchronous generator model prediction control method and system |
-
2021
- 2021-08-11 CN CN202110917689.4A patent/CN113572199B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109950922A (en) * | 2019-01-31 | 2019-06-28 | 东南大学 | A kind of multistep model predictive control method suitable for VSC-HVDC |
| WO2021035269A1 (en) * | 2019-08-30 | 2021-03-04 | Avl List Gmbh | Method and controller for model predictive control of coupled converters |
| CN113179059A (en) * | 2021-05-21 | 2021-07-27 | 南京理工大学 | Improved virtual synchronous generator model prediction control method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113572199A (en) | 2021-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113572199B (en) | Smooth switching method of network-forming type current converter based on model prediction algorithm | |
| Tungadio et al. | Load frequency controllers considering renewable energy integration in power system | |
| Zhang et al. | Data-driven control for interlinked AC/DC microgrids via model-free adaptive control and dual-droop control | |
| CN108923460B (en) | Parameter configuration method for multi-machine parallel dynamic consistent response of micro-grid virtual synchronous machine | |
| WO2023088124A1 (en) | Model prediction-based frequency self-adaptive control method for virtual synchronizer inverter | |
| Zhang et al. | Optimal frequency control for virtual synchronous generator based AC microgrids via adaptive dynamic programming | |
| He et al. | An adaptive VSG control strategy of battery energy storage system for power system frequency stability enhancement | |
| Duan et al. | Hierarchical power flow control in smart grids: Enhancing rotor angle and frequency stability with demand-side flexibility | |
| Fan et al. | Nonlinear model predictive control of HVDC for inter-area oscillation damping | |
| WO2020214913A1 (en) | Systems and methods for power system stabilization and oscillation damping control | |
| Lyu et al. | Fuzzy control based virtual synchronous generator for self-adaptative control in hybrid microgrid | |
| Beus et al. | A model predictive control approach to operation optimization of an ultracapacitor bank for frequency control | |
| Mc Namara et al. | Control strategies for automatic generation control over MTDC grids | |
| Li et al. | A linear quadratic regulator with integral action of wind turbine based on aerodynamics forecasting for variable power production | |
| CN112350352A (en) | Method for increasing energy storage reactive power regulation rate | |
| CN117526355A (en) | Distributed photovoltaic and adjustable load participation grid frequency modulation method and system | |
| Subramanian et al. | A cooperative rate-based model predictive framework for flexibility management of DERs | |
| Farahani et al. | Intelligent control of static synchronous series compensator via an adaptive self-tuning PID controller for suppression of torsional oscillations | |
| Del Rivero et al. | Control of aggregated virtual synchronous generators including communication delay compensation | |
| CN113162063B (en) | Design method of multi-direct-current coordination controller for inhibiting ultralow frequency oscillation | |
| Jiang et al. | Secondary reactive power balancing and voltage stability in microgrid using prioritized centralized controller | |
| Zhao et al. | A Smooth Grid-connected Strategy for Photovoltaic-storage Virtual Synchronous Generator Based on Adaptive Model Predictive Control | |
| Kammer et al. | Robust distributed averaging frequency control of inverter-based microgrids | |
| CN103855716A (en) | Intelligent FLC-PID mixed STATCOM control method | |
| Zhang et al. | Adaptive grid-forming strategy for a photovoltaic storage system based on edge transfer PSO algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |