CN108539798A - Energy-storage system Secondary Control strategy based on Model Predictive Control - Google Patents
Energy-storage system Secondary Control strategy based on Model Predictive Control Download PDFInfo
- Publication number
- CN108539798A CN108539798A CN201810349876.5A CN201810349876A CN108539798A CN 108539798 A CN108539798 A CN 108539798A CN 201810349876 A CN201810349876 A CN 201810349876A CN 108539798 A CN108539798 A CN 108539798A
- Authority
- CN
- China
- Prior art keywords
- voltage
- control
- angular frequency
- variable
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004146 energy storage Methods 0.000 title claims abstract description 17
- 238000011217 control strategy Methods 0.000 title abstract 2
- 238000005457 optimization Methods 0.000 claims abstract description 25
- 238000005096 rolling process Methods 0.000 claims abstract description 4
- 230000001105 regulatory effect Effects 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 14
- 230000010354 integration Effects 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000003990 capacitor Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000000470 constituent Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 238000009828 non-uniform distribution Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 229960001948 caffeine Drugs 0.000 claims 1
- 230000001276 controlling effect Effects 0.000 claims 1
- 230000005611 electricity Effects 0.000 claims 1
- 239000012467 final product Substances 0.000 claims 1
- RYYVLZVUVIJVGH-UHFFFAOYSA-N trimethylxanthine Natural products CN1C(=O)N(C)C(=O)C2=C1N=CN2C RYYVLZVUVIJVGH-UHFFFAOYSA-N 0.000 claims 1
- 230000001052 transient effect Effects 0.000 abstract description 5
- 238000000034 method Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
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/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Electrical Variables (AREA)
- Inverter Devices (AREA)
Abstract
The invention discloses a kind of energy-storage system Secondary Control strategy based on Model Predictive Control.Existing frequency, the PI controls of voltage Secondary Control scheme generally use, in load fluctuation, frequency can not be stablized in reference frequency value, and frequency and the transient response of voltage are slower.The technical solution adopted by the present invention is:First, the inverter of energy-storage system and its state-space model of controlling unit are derived, secondly, establishes the prediction model of angular frequency and voltage;Then through rolling optimization and feedback compensation, the control variable of angular frequency and voltage, is finally respectively added on active power/angular frequency and reactive power/voltage droop control by the control variable for obtaining angular frequency and voltage, realize that frequency retrieval to the precise proportions of a reference value and reactive power distribute.There is the present invention faster frequency, voltage transient to respond, and control accuracy is high, highly practical.
Description
Technical Field
The invention belongs to the technical field of distributed power generation and micro-grids, and particularly relates to a secondary voltage regulation and frequency modulation strategy of an energy storage system based on model predictive control.
Background
The micro-grid is composed of devices such as a Distributed Generator (DG) including an energy storage system and a diesel generator, a power converter and a load, and impact caused when the micro-grid is connected to a power grid can be relieved through interactive supplementation with a large power grid. When a microgrid is in an island mode, a plurality of DGs are usually operated in parallel to provide power for loads, and droop control with a plug and play characteristic is the most common control mode. The traditional P-f droop characteristic cannot stabilize the frequency at a rated value when the load changes, and the line impedance in the low-voltage microgrid is usually resistive or resistive, so that the traditional Q-V droop control cannot realize proportional distribution of reactive power, and each output voltage deviates from a voltage reference value.
In order to improve the distribution accuracy of reactive power, some scholars introduce virtual impedance at the output end of an inverter, change the characteristics of output impedance and improve the power distribution accuracy of the inverter, but the virtual impedance cannot solve the influence of different output lines on the current sharing of a system, so that some scholars adopt a secondary regulation strategy, and add secondary regulation optimization quantities on reference values of frequency and voltage respectively to realize the frequency recovery to a rated value and the proportional distribution of the reactive power according to the rated capacity; however, the traditional secondary regulation strategy adopts PI control, which affects the transient characteristics of frequency and voltage when the load fluctuates.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a secondary regulation strategy of an energy storage system based on model predictive control so as to realize the optimal control of the voltage and the frequency of the energy storage system.
Therefore, the invention adopts the following technical scheme: the model predictive control-based energy storage system secondary regulation strategy comprises the following steps:
deducing a state space model of the inverter and a control link of the inverter, and establishing a prediction model of angular frequency and voltage;
respectively obtaining an optimized objective function of angular frequency and voltage according to an expected value, obtaining a plurality of control increments when the function takes a minimum value through rolling optimization of the optimized objective function, obtaining an output predicted value by combining an initial predicted value at the last sampling moment, and obtaining a corrected output predicted value after feedback correction with an actual measured value to be used as an initial predicted value at the next moment; and simultaneously, taking the instant control increment to form a control variable as a secondary regulating variable for regulating the angular frequency and the voltage.
Because the PI controller can only utilize past and current measured values, and adopts a secondary regulation strategy based on PI control, transient characteristics are easily influenced when the frequency and voltage of an energy storage device which is connected into a microgrid and serves as a distributed power supply are regulated. Research and analysis have found that Model Predictive Control (MPC) predicts future information using not only current and past measurements but also a predictive model, adjusts control variables through roll optimization, and feedback corrects the predicted output to maintain the output at a desired value, taking into account the constraints associated with the controlled object.
The method utilizes the characteristic that Model Predictive Control (MPC) can predict future process information and accurately track an expected value, establishes a prediction model of angular frequency and voltage by deducing a state space model of an inverter of an energy storage system and a control link thereof, obtains control variables of the angular frequency and the voltage through rolling optimization and feedback correction, and realizes frequency recovery to a reference value and accurate proportional distribution of reactive power as secondary adjustment optimization quantities of active power/angular frequency and reactive power/voltage droop characteristics. Compared with a secondary regulation strategy of PI control, the secondary regulation strategy based on MPC has faster frequency and voltage transient response characteristics when the load fluctuates, the frequency is always maintained at a frequency reference value when the load fluctuates, the deviation of each output voltage and the reference voltage is greatly reduced, and meanwhile, the reactive power is accurately distributed according to the rated capacity proportion.
As a supplement to the above technical solution, the step of constructing the prediction model of angular frequency and voltage includes:
step 1): sampling to obtain the voltage u at the output end of the inverteriAnd current iiThe voltage u at the output end of the filteroAnd current ioVoltage of common connection point is Vs;uoConversion into dq-axis variable u by parkod、uoqIn the same way, get io、ui、iiAnd VsDq axis variable of (1); the angular frequency of the output voltage is set to omega when the inverter operates stably0Derivation of ii、uoAnd ioThe equation of state of (a) is:
and
wherein L isf、CfAnd RfRespectively a filter inductor, a capacitor and a resistor; l, R are respectively line inductance and resistance; omegan=2πfnAt a nominal angular frequency, wherein fn=50Hz;
Step 2): calculating active power and reactive power, p ═ u respectivelyodiod+uoqioq,q=uoqiod-uodioqWith a cut-off angular frequency of ωcThe low-pass filter has an average power of P, Q, i.e.s is a Laplace transform differential operator;
using P, Q to obtain droop control parameters of angular frequency and voltage, and adding secondary regulation quantity to obtainWhere ω is 2 π f, ωn=2πfn,fn、UnReference values of frequency and voltage are respectively; pn、QnRespectively corresponding active power and reactive power when the frequency and the voltage are at reference values; m and n are P-omega and Q-V droop coefficients respectively; omegas、uMRespectively an added angular frequency variable and a voltage variable; using the voltage value generated by droop control as the reference of the voltage control outer loopSignals, i.e.
u* od=Un+n(Qn-Q)+uM、u* oq=0;
Step 3): defining a voltage difference state variableNamely, pair (u)* od-uod) Integration is performed to obtain a variableSame reason pair (u)* oq-uoq) Integral derived variableDesign voltage outer loop control equationWherein k isvp、kviRespectively are the proportional and integral parameters of the voltage outer ring; h is the regulated output current ioThe current signal i generated by the voltage outer loop control equation* id、i* iqAs a current inner loop reference signal;
defining a current difference state variableNamely pair (i)* id-iid) Integration to obtain the variable gammadThe same thing is to (i)* iq-iiq) Integration to obtain the variable gammaq. Designing the control equation of the current inner loop asWherein k isip、kiiRespectively are the proportion and integral parameters of the inner ring; neglecting the dynamic influence of the switching part at higher frequencies, the voltage signal u generated by the current inner loop control equation is considered* id、u* iqAre respectively equal to inversionOutput voltage u of the deviceid、uiq;
Step 4): from ω to ωn+m(Pn-P)+ωsObtaining an angular frequency, and obtaining an angle δ of a dq axis coordinate defined by the inverter through integration, namely δ ═ ω dt, wherein all variables are defined under the dq axis coordinate system, and a common bus voltage is (V)d,Vq) Will (V)d,Vq) Transforming to dq axis coordinate system defined by the inverter to obtain
According to delta, P, Q,γdq(γd、γq)、iidq(iid、iiq)、uodq(uod、uoq)、iodq(iod、ioq) Establishing small signal state space model of inverter and control link thereofWhereinIs a state variable, Δ u ═ Δ ωsΔuMΔVdΔVq]TIn order to control the variables of the plant,
matrix ADG、BDGIn (I)id、Iiq、Iod、Ioq、Vod、VoqAre respectively state space variablesiid、iiq、iod、ioq、uod、uoqCorresponding state values when the inverter operates stably;
step 5): obtaining the frequency and voltage U output by each distributed power supply by using a central controller, and calculating the average angular frequencyAnd average voltageEstablishing an angular frequency prediction model with an output prediction value as an average angular frequencyWherein, Cω=[0-m 0 0 0 0 0 0 0 0 0 0 0],Dω=[1 0 0],Δuω=[ΔωsΔVdΔVq]T,
Matrix BDGωIn (II)id、Iiq、Iod、Ioq、Vod、VoqAre respectively a state space variable iid、iiq、iod、ioq、uod、uoqCorresponding state values when the inverter operates stably;
establishing a voltage prediction model with an output prediction value as an average voltage
Wherein, CU=[0 0–n 0 0 0 0 0 0 0 0 0 0],DU=[1 0 0],ΔuU=[ΔuMΔVdΔVq]T,
With a sampling period TsDiscretizing the prediction models of the angular frequency and the voltage to respectively obtain
Wherein
In the formula, the triangular symbol Δ represents a small signal change amount of a certain variable.
In addition to the above-described solutions, the desired values ω are respectively specifiedn=2πfn,fn50Hz and Un311V, the target optimization function for angular frequency is Wω、WsRespectively representing the weight coefficients for suppressing the tracking error and the control increment change, LωRepresenting the predicted step size, M, of the angular frequency objective optimization functionωRepresenting the control step with the constraint ofωsmax=2π×0.5rad/s、ωsminEach of-2 π × 0.5rad/s is a control variable ωs(k) Upper and lower limits of (d); control increment Δ ωs(k) Respectively is Δ ωsmax、Δωsmin(ii) a The upper limit and the lower limit of the predicted output quantity are respectively Determination of a control increment [ Delta omega ] from an objective optimization function of angular frequencys(k),…,Δωs(k+Mω-1)]TTaking the real-time control increment delta omegas(k) Angular frequency control variable omega with last sampling moments(k-1) constitutes a k-time controlled variable ωs(k) And adding the angular frequency control circuit into active power/angular frequency droop control to realize secondary regulation of the angular frequency.
In addition to the above technical solution, the objective optimization function of the voltage is
WU、WuRespectively representing the weight coefficients for suppressing the tracking error and the control increment change, LURepresenting the predicted step size, M, of the voltage objective optimization functionURepresents a control step size; the constraint condition is
The maximum deviation allowed by the output voltage is 5 percent of the reference value, and u is takenMmin=-15V,uMmax15V, thenDetermination of a control increment [ Delta u ] from a target optimization function of the voltageM(k),…,ΔuM(k+MU-1)]TTaking the real-time control increment delta uM(k) With the voltage control variable u at the last sampling instantM(k-1) constituent k time controlled variable uM(k) Added to the reactive power/voltage droop control to obtain a resultant voltage, i.e. Un+n(Qn-Q)+uMAnd the influence of load fluctuation on the deviation of the output voltage from the reference voltage is reduced.
In addition to the above-described solution, the deviation of the combined voltage from the actually measured average voltage is PI-regulated, i.e. the deviation isIn the formula, kQpFor proportional coefficient of PI regulation, kQiFor integral coefficient of PI regulation, s is Laplace transform differential operator, and the obtained regulating quantity delta UQAnd adding the voltage to the synthesized voltage to serve as a final voltage reference value so as to inhibit the influence of line impedance on the non-uniform distribution of the reactive power and realize the proportional distribution of the reactive power according to the rated capacity.
The method utilizes the characteristic that a model predictive control algorithm has accurate tracking on an expected value to obtain the secondary regulating quantity of angular frequency and voltage, and respectively adds the secondary regulating quantity into active power/angular frequency droop control and reactive power/voltage droop control, so that the system frequency is always maintained at a frequency reference value while the reactive power is distributed according to the rated capacity proportion, and the deviation of each output voltage and the reference voltage is greatly reduced.
Detailed Description
The present invention will be further described with reference to the following embodiments.
The embodiment provides an energy storage system secondary regulation strategy based on model predictive control, which comprises the following steps:
step 1): sampling to obtain the voltage u at the output end of the inverteriAnd current iiThe voltage u at the output end of the filteroAnd current ioVoltage of common connection point is Vs;uoConversion into dq-axis variable u by parkod、uoqIn the same way, get io、ui、iiAnd VsDq axis variable of (1); the angular frequency of the output voltage is set to omega when the inverter operates stably0Derivation of ii、uoAnd ioThe equation of state of (a) is:
and
wherein L isf、CfAnd RfRespectively a filter inductor, a capacitor and a resistor; l, R are respectively line inductance and resistance; omegan=2πfnAt a nominal angular frequency, wherein fn=50Hz。
Step 2): calculating active power and reactive power, p ═ u respectivelyodiod+uoqioq,q=uoqiod-uodioqWith a cut-off angular frequency of ωcThe low-pass filter has an average power of P, Q, i.e.
s is a Laplace transform differential operator;
using P, Q to obtain droop control parameters of angular frequency and voltage, and adding secondary regulation quantity to obtainWhere ω is 2 π f, ωn=2πfn,fn、UnReference values of frequency and voltage are respectively; pn、QnRespectively corresponding active power and reactive power when the frequency and the voltage are at reference values; m and n are P-omega and Q-V droop coefficients respectively; omegas、uMRespectively an added angular frequency variable and a voltage variable; the voltage value generated by droop control is used as a reference signal for the outer loop of the voltage control, i.e.
u* od=Un+n(Qn-Q)+uM、u* oq=0。
Step 3): defining a voltage difference state variableNamely, pair (u)* od-uod) Integration is performed to obtain a variableSame reason pair (u)* oq-uoq) Integral derived variableDesign voltage outer loop control equationWherein k isvp、kviRespectively are the proportional and integral parameters of the voltage outer ring; h is the regulated output current ioThe current signal i generated by the voltage outer loop control equation* id、i* iqAs a current inner loop reference signal;
defining a current difference state variableNamely pair (i)* id-iid) Integration to obtain the variable gammadThe same thing is to (i)* iq-iiq) Integration to obtain the variable gammaq. Designing the control equation of the current inner loop as
Wherein k isip、kiiRespectively are the proportion and integral parameters of the inner ring; neglecting the dynamic influence of the switching part at higher frequencies, the voltage signal u generated by the current inner loop control equation is considered* id、u* iqAre respectively equal to the output voltage u of the inverterid、uiq。
Step 4): from ω to ωn+m(Pn-P)+ωsObtaining an angular frequency, and obtaining an angle δ of a dq axis coordinate defined by the inverter through integration, namely δ ═ ω dt, wherein all variables are defined under the dq axis coordinate system, and a common bus voltage is set to be ^ ω dt(Vd,Vq) Will (V)d,Vq) Transforming to dq axis coordinate system defined by the inverter to obtain
According to delta, P, Q,γdq(γd、γq)、iidq(iid、iiq)、uodq(uod、uoq)、iodq(iod、ioq) Establishing small signal state space model of inverter and control link thereofWhereinIs a state variable, Δ u ═ Δ ωsΔuMΔVdΔVq]TIn order to control the variables of the plant,
matrix ADG、BDGIn (I)id、Iiq、Iod、Ioq、Vod、VoqAre respectively a state space variable iid、iiq、iod、ioq、uod、uoqAnd corresponding state values when the inverter operates stably.
Step 5): obtaining the frequency and voltage U output by each distributed power supply by using a central controller, and calculating the average angular frequencyAnd average voltageEstablishing an angular frequency prediction model with an output prediction value as an average angular frequencyWherein, Cω=[0-m 0 0 0 0 0 0 0 0 0 0 0],Dω=[1 0 0],Δuω=[ΔωsΔVdΔVq]T,
Matrix BDGωIn (II)id、Iiq、Iod、Ioq、Vod、VoqAre respectively a state space variable iid、iiq、iod、ioq、uod、uoqAnd corresponding state values when the inverter operates stably.
Establishing a voltage prediction model with an output prediction value as an average voltage
Wherein, CU=[0 0–n 0 0 0 0 0 0 0 0 0 0],DU=[1 0 0],ΔuU=[ΔuMΔVdΔVq]T,
With a sampling period TsDiscretizing the prediction models of the angular frequency and the voltage to respectively obtain,
wherein
In the formula, the triangular symbol Δ represents a small signal change amount of a certain variable.
Respectively given desired values ωn=2πfn,fn50Hz and Un311V, the target optimization function for angular frequency isWω、WsRespectively representing the weight coefficients for suppressing the tracking error and the control increment change, LωRepresenting the predicted step size, M, of the angular frequency objective optimization functionωIndicating the control step size.
The constraint condition isωsmax=2π×0.5rad/s、ωsminEach of-2 π × 0.5rad/s is a control variable ωs(k) Upper and lower limits of (d); control increment Δ ωs(k) Respectively is Δ ωsmax、Δωsmin(ii) a The upper limit and the lower limit of the predicted output quantity are respectively Determination of a control increment [ Delta omega ] from an objective optimization function of angular frequencys(k),…,Δωs(k+Mω-1)]TTaking the real-time control increment delta omegas(k) Angular frequency control variable omega with last sampling moments(k-1) constitutes a k-time controlled variable ωs(k) And adding the angular frequency control circuit into active power/angular frequency droop control to realize secondary regulation of the angular frequency.
The objective optimization function of the voltage is
WU、WuRespectively representing the weight coefficients for suppressing the tracking error and the control increment change, LURepresenting the predicted step size, M, of the voltage objective optimization functionUIndicating the control step size. The constraint condition is
The maximum deviation allowed by the output voltage is 5 percent of the reference value, and u is takenMmin=-15V,uMmax15V, thenDetermination of a control increment [ Delta u ] from a target optimization function of the voltageM(k),…,ΔuM(k+MU-1)]TTaking the real-time control increment delta uM(k) With the voltage control variable u at the last sampling instantM(k-1) constituent k time controlled variable uM(k) Added to the reactive power/voltage droop control to obtain a resultant voltage, i.e. Un+n(Qn-Q)+uMAnd the influence of load fluctuation on the deviation of the output voltage from the reference voltage is reduced.
The deviation of the resulting voltage from the actually measured average voltage is subjected to PI regulation (proportional-integral regulation), i.e.In the formula, kQpFor proportional coefficient of PI regulation, kQiFor the integral coefficient of PI regulation, s is the Laplace transform differential operator. The adjustment amount DeltaU to be obtainedQAnd adding the voltage to the synthesized voltage to serve as a final voltage reference value so as to inhibit the influence of line impedance on the non-uniform distribution of the reactive power and realize the proportional distribution of the reactive power according to the rated capacity.
The above-mentioned embodiments only express the embodiments of the present invention, and therefore, should not be interpreted as limiting the scope of the present invention, and should not be interpreted as limiting the structure of the present invention in any way. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention.
Claims (5)
1. The energy storage system secondary regulation strategy based on model predictive control is characterized by comprising the following steps:
deducing a state space model of the inverter and a control link of the inverter, and establishing a prediction model of angular frequency and voltage;
respectively obtaining an optimized objective function of angular frequency and voltage according to an expected value, obtaining a plurality of control increments when the function takes a minimum value through rolling optimization of the optimized objective function, obtaining an output predicted value by combining an initial predicted value at the last sampling moment, and obtaining a corrected output predicted value after feedback correction with an actual measured value to be used as an initial predicted value at the next moment; and simultaneously, taking the instant control increment to form a control variable as a secondary regulating variable for regulating the angular frequency and the voltage.
2. The model predictive control-based energy storage system quadratic regulation strategy of claim 1, wherein the step of constructing the predictive model of angular frequency and voltage comprises:
step 1): sampling to obtain the voltage u at the output end of the inverteriAnd current iiThe voltage u at the output end of the filteroAnd current ioVoltage of common connection point is Vs;uoConversion into dq-axis variable u by parkod、uoqIn the same way, get io、ui、iiAnd VsDq axis variable of (1); the angular frequency of the output voltage is set to omega when the inverter operates stably0Derivation of ii、uoAnd ioThe equation of state of (a) is:
and
wherein L isf、CfAnd RfRespectively a filter inductor, a capacitor and a resistor; l, R are respectively line inductance and resistance; omegan=2πfnAt a nominal angular frequency, wherein fn=50Hz;
Step 2): calculating active power and reactive power, p ═ u respectivelyodiod+uoqioq,q=uoqiod-uodioqWith a cut-off angular frequency of ωcThe low-pass filter has an average power of P, Q, i.e.s is a Laplace transform differential operator;
by usingP, Q obtaining droop control parameters of angular frequency and voltage, and adding secondary adjustment amount to obtain final productWhere ω is 2 π f, ωn=2πfn,fn、UnReference values of frequency and voltage are respectively; pn、QnRespectively corresponding active power and reactive power when the frequency and the voltage are at reference values; m and n are P-omega and Q-V droop coefficients respectively; omegas、uMRespectively an added angular frequency variable and a voltage variable; the voltage value generated by droop control is used as a reference signal of a voltage control outer ring, i.e. u* od=Un+n(Qn-Q)+uM、u* oq=0;
Step 3): defining a voltage difference state variableNamely, pair (u)* od-uod) Integration is performed to obtain a variableSame reason pair (u)* oq-uoq) Integral derived variableDesign voltage outer loop control equationWherein k isvp、kviRespectively are the proportional and integral parameters of the voltage outer ring; h is the regulated output current ioThe current signal i generated by the voltage outer loop control equation* id、i* iqAs a current inner loop reference signal;
defining a current difference state variableNamely pair (i)* id-iid) Integration to obtain the variable gammadThe same thing is to (i)* iq-iiq) Integration to obtain the variable gammaqDesigning the control equation of the current inner loop asWherein k isip、kiiRespectively are the proportion and integral parameters of the inner ring; neglecting the dynamic influence of the switching part at higher frequencies, the voltage signal u generated by the current inner loop control equation is considered* id、u* iqAre respectively equal to the output voltage u of the inverterid、uiq;
Step 4): from ω to ωn+m(Pn-P)+ωsObtaining an angular frequency, and obtaining an angle δ of a dq axis coordinate defined by the inverter through integration, namely δ ═ ω dt, wherein all variables are defined in the dq axis coordinate system, and a common bus voltage is (V)d,Vq) Will (V)d,Vq) Transforming to dq axis coordinate system defined by the inverter to obtainAccording to delta, P, Q,γd、γq、iid、iiq、uod、uoq、iod、ioqEstablishing small signal state space model of inverter and control link thereofWhereinIn order to be a state variable, the state variable,
Δu=[ΔωsΔuMΔVdΔVq]Tfor controlling variables, matrix ADG、BDGIn (I)id、Iiq、Iod、Ioq、Vod、VoqAre respectively a state space variable iid、iiq、iod、ioq、uod、uoqCorresponding state values when the inverter operates stably;
step 5): obtaining the frequency and voltage U output by each distributed power supply by using a central controller, and calculating the average angular frequencyAnd average voltageEstablishing an angular frequency prediction model with an output prediction value as an average angular frequency
Wherein, Cω=[0 -m 0 0 0 0 0 0 0 0 0 0 0],Dω=[1 0 0],Δuω=[ΔωsΔVdΔVq]T,
Matrix BDGωIn (II)id、Iiq、Iod、Ioq、Vod、VoqAre respectively a state space variable iid、iiq、iod、ioq、uod、uoqCorresponding state values when the inverter operates stably;
electricity with output prediction value as average voltagePressure prediction model
Wherein, CU=[0 0 –n 0 0 0 0 0 0 0 0 0 0],DU=[1 0 0],ΔuU=[ΔuMΔVdΔVq]T,
With a sampling period TsDiscretizing the prediction models of the angular frequency and the voltage to respectively obtain
Wherein,
in the formula, the triangular symbol Δ represents a small signal change amount of a certain variable.
3. The model predictive control-based energy storage system quadratic regulation strategy of claim 2, characterized in that a desired value ω is respectively givenn=2πfn,fn50Hz and Un311V, the target optimization function for angular frequency isWω、WsRespectively representing the weight coefficients for suppressing the tracking error and the control increment change, LωRepresenting the predicted step size, M, of the angular frequency objective optimization functionωRepresenting the control step with the constraint ofωsmax=2π×0.5rad/s、ωsminControl is carried out for-2 pi x 0.5rad/s respectivelySystem variable omegas(k) Upper and lower limits of (d); control increment Δ ωs(k) Respectively is Δ ωsmax、Δωsmin(ii) a The upper limit and the lower limit of the predicted output quantity are respectively Determination of a control increment [ Delta omega ] from an objective optimization function of angular frequencys(k),…,Δωs(k+Mω-1)]TTaking the real-time control increment delta omegas(k) Angular frequency control variable omega with last sampling moments(k-1) constitutes a k-time controlled variable ωs(k) And adding the angular frequency control circuit into active power/angular frequency droop control to realize secondary regulation of the angular frequency.
4. The energy storage system secondary regulation strategy of claim 3, wherein the objective optimization function of voltage isWU、WuRespectively representing the weight coefficients for suppressing the tracking error and the control increment change, LURepresenting the predicted step size, M, of the voltage objective optimization functionURepresenting the control step with the constraint ofThe maximum deviation allowed by the output voltage is 5 percent of the reference value, and u is takenMmin=-15V,uMmax15V, thenDetermination of a control increment [ Delta u ] from a target optimization function of the voltageM(k),…,ΔuM(k+MU-1)]TTaking the real-time control increment delta uM(k) With the voltage control variable u at the last sampling instantM(k-1) constituent k time controlled variable uM(k) Added to the reactive power/voltage droop control to obtain a resultant voltage, i.e. Un+n(Qn-Q)+uMAnd the influence of load fluctuation on the deviation of the output voltage from the reference voltage is reduced.
5. The energy storage system secondary regulation strategy of claim 4, characterized in that the deviation of the resulting voltage from the actual measured average voltage is PI regulated, i.e. theIn the formula, kQpFor proportional coefficient of PI regulation, kQiFor integral coefficient of PI regulation, s is Laplace transform differential operator, and obtained regulating quantity delta UQAnd adding the voltage to the synthesized voltage to serve as a final voltage reference value so as to inhibit the influence of line impedance on the non-uniform distribution of the reactive power and realize the proportional distribution of the reactive power according to the rated capacity.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810349876.5A CN108539798B (en) | 2018-04-18 | 2018-04-18 | Secondary regulation strategy of energy storage system based on model predictive control |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810349876.5A CN108539798B (en) | 2018-04-18 | 2018-04-18 | Secondary regulation strategy of energy storage system based on model predictive control |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108539798A true CN108539798A (en) | 2018-09-14 |
CN108539798B CN108539798B (en) | 2020-12-18 |
Family
ID=63477874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810349876.5A Active CN108539798B (en) | 2018-04-18 | 2018-04-18 | Secondary regulation strategy of energy storage system based on model predictive control |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108539798B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109217709A (en) * | 2018-10-15 | 2019-01-15 | 深圳市安和威电力科技股份有限公司 | Bi-directional power conversion AC-DC control system and method based on IGBT |
CN109598088A (en) * | 2018-12-20 | 2019-04-09 | 中国矿业大学 | A kind of belt conveyor belt speed setting control method |
CN110148956A (en) * | 2019-05-07 | 2019-08-20 | 万克能源科技有限公司 | A kind of battery energy storage system auxiliary AGC control method based on MPC |
CN110855155A (en) * | 2019-12-04 | 2020-02-28 | 兰州交通大学 | Screen grid power supply control method based on model predictive control |
CN110932336A (en) * | 2019-11-25 | 2020-03-27 | 云南电网有限责任公司临沧供电局 | Low-voltage distribution station voltage control method and system based on model predictive control |
CN112350352A (en) * | 2020-11-20 | 2021-02-09 | 西安热工研究院有限公司 | Method for increasing energy storage reactive power regulation rate |
CN113054690A (en) * | 2021-03-17 | 2021-06-29 | 华翔翔能科技股份有限公司 | Voltage control method and system based on event triggering and electronic equipment |
CN113629748A (en) * | 2021-10-11 | 2021-11-09 | 国网江西省电力有限公司电力科学研究院 | Five-level energy storage converter cascade model prediction control method and device |
CN116466287A (en) * | 2023-06-20 | 2023-07-21 | 贵州海纳储能技术有限公司 | Automatic calibration method for on-line inverter parallel system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5612859A (en) * | 1995-05-03 | 1997-03-18 | Lucent Technologies Inc. | Battery polarity switch for generating low frequency voltage waveforms |
CN104377826A (en) * | 2014-11-19 | 2015-02-25 | 云南电网公司电力科学研究院 | Active power distribution network control strategy and method |
-
2018
- 2018-04-18 CN CN201810349876.5A patent/CN108539798B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5612859A (en) * | 1995-05-03 | 1997-03-18 | Lucent Technologies Inc. | Battery polarity switch for generating low frequency voltage waveforms |
CN104377826A (en) * | 2014-11-19 | 2015-02-25 | 云南电网公司电力科学研究院 | Active power distribution network control strategy and method |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109217709A (en) * | 2018-10-15 | 2019-01-15 | 深圳市安和威电力科技股份有限公司 | Bi-directional power conversion AC-DC control system and method based on IGBT |
CN109598088B (en) * | 2018-12-20 | 2022-02-11 | 中国矿业大学 | Belt speed setting control method for belt conveyor |
CN109598088A (en) * | 2018-12-20 | 2019-04-09 | 中国矿业大学 | A kind of belt conveyor belt speed setting control method |
CN110148956A (en) * | 2019-05-07 | 2019-08-20 | 万克能源科技有限公司 | A kind of battery energy storage system auxiliary AGC control method based on MPC |
CN110932336A (en) * | 2019-11-25 | 2020-03-27 | 云南电网有限责任公司临沧供电局 | Low-voltage distribution station voltage control method and system based on model predictive control |
CN110932336B (en) * | 2019-11-25 | 2020-11-13 | 云南电网有限责任公司临沧供电局 | Low-voltage distribution station voltage control method and system based on model predictive control |
CN110855155A (en) * | 2019-12-04 | 2020-02-28 | 兰州交通大学 | Screen grid power supply control method based on model predictive control |
CN110855155B (en) * | 2019-12-04 | 2021-06-18 | 兰州交通大学 | Screen grid power supply control method based on model predictive control |
CN112350352A (en) * | 2020-11-20 | 2021-02-09 | 西安热工研究院有限公司 | Method for increasing energy storage reactive power regulation rate |
CN113054690A (en) * | 2021-03-17 | 2021-06-29 | 华翔翔能科技股份有限公司 | Voltage control method and system based on event triggering and electronic equipment |
CN113629748A (en) * | 2021-10-11 | 2021-11-09 | 国网江西省电力有限公司电力科学研究院 | Five-level energy storage converter cascade model prediction control method and device |
CN116466287A (en) * | 2023-06-20 | 2023-07-21 | 贵州海纳储能技术有限公司 | Automatic calibration method for on-line inverter parallel system |
CN116466287B (en) * | 2023-06-20 | 2023-09-22 | 贵州海纳储能技术有限公司 | Automatic calibration method for on-line inverter parallel system |
Also Published As
Publication number | Publication date |
---|---|
CN108539798B (en) | 2020-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108539798B (en) | Secondary regulation strategy of energy storage system based on model predictive control | |
KR102236082B1 (en) | Microgrid dynamic stability control system and microgrid dynamic stability control method | |
CN109687526B (en) | Island micro-grid layered distributed control strategy based on consistency theory | |
WO2020252813A1 (en) | Double-layer adaptive inertia control method and device for inverter interfaced distributed generator | |
US9588557B2 (en) | Reactive following for distributed generation and loads of other reactive controller(s) | |
CN111541274A (en) | Island microgrid control strategy based on virtual synchronous generator characteristics | |
CN108199396B (en) | Virtual excitation closed-loop control system of energy storage inverter and design method thereof | |
US20180076622A1 (en) | Expanded Reactive Following for Distributed Generation and Loads of Other Reactive Controller(s) | |
CN110460112B (en) | Micro-grid sag control method based on power deviation value consistency control | |
CN103259268B (en) | Microgrid reactive-voltage control device and control method thereof | |
CN105024390B (en) | Micro-grid battery energy storage system frequency modulation control method based on BP neural network | |
KR101520248B1 (en) | Method and Apparatus for Controlling Doubly-fed Induction Generator using Adaptive Backstepping Control Scheme | |
CN107612025A (en) | Current-control type inverter improves control method in microgrid | |
CN112260287A (en) | Novel photovoltaic microgrid droop reactive compensation control method | |
US12044214B2 (en) | Relating to voltage control in renewable energy power plants | |
CN110350538B (en) | Micro-grid coordination control method based on active demand side response | |
CN115589036A (en) | Method for adjusting parallel active power equalization of multiple inverters based on self-adaptive virtual impedance | |
CN109586305B (en) | Power distribution network operation control strategy based on flexible multi-state switch | |
CN113328429B (en) | Parallel-off-grid smooth switching and off-grid rolling optimized scheduling method, equipment and medium | |
CN112636392B (en) | Single-stage multi-terminal hybrid micro-grid structure suitable for low-voltage house and control method thereof | |
CN112310972B (en) | Coordinated control method for AVC system of master station and static reactive compensation device of station end | |
CN112350382A (en) | Energy storage active power adjusting method with angular frequency compensation module | |
CN112350352A (en) | Method for increasing energy storage reactive power regulation rate | |
CN115549504A (en) | Control method of three-level energy storage converter | |
CN113765160A (en) | Distributed energy and grid-connected reactive power control method and device thereof, and storage medium |
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 |