CN115065077A - Control method and device for energy storage voltage type current converter - Google Patents

Control method and device for energy storage voltage type current converter Download PDF

Info

Publication number
CN115065077A
CN115065077A CN202210944503.9A CN202210944503A CN115065077A CN 115065077 A CN115065077 A CN 115065077A CN 202210944503 A CN202210944503 A CN 202210944503A CN 115065077 A CN115065077 A CN 115065077A
Authority
CN
China
Prior art keywords
mpc
mathematical model
energy storage
power
output
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.)
Withdrawn
Application number
CN202210944503.9A
Other languages
Chinese (zh)
Inventor
何敏强
杨沛豪
孙钢虎
兀鹏越
王小辉
柴琦
寇水潮
高峰
关国璋
刘鹏峰
薛国龙
胡世铭
詹新民
林开俊
刘文武
杨振
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Thermal Power Research Institute Co Ltd
Huaneng Luoyuan Power Generation Co Ltd
Original Assignee
Xian Thermal Power Research Institute Co Ltd
Huaneng Luoyuan Power Generation Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Thermal Power Research Institute Co Ltd, Huaneng Luoyuan Power Generation Co Ltd filed Critical Xian Thermal Power Research Institute Co Ltd
Priority to CN202210944503.9A priority Critical patent/CN115065077A/en
Publication of CN115065077A publication Critical patent/CN115065077A/en
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241The oscillation concerning frequency
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/50Controlling the sharing of the out-of-phase component
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Power Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Ac-Ac Conversion (AREA)

Abstract

The application provides a control method of an energy storage voltage type converter, which comprises the following steps: acquiring a VSG rotor motion equation of the energy storage voltage type converter, and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model; converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation at a second preset moment; establishing a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation to minimize an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, and carrying out VSG power control according to the target active power reference value to realize MPC control of the energy storage voltage type converter. The method and the device adopt a two-cycle delay compensation control strategy to carry out advanced control on system variables, and improve control precision.

Description

Control method and device for energy storage voltage type current converter
Technical Field
The present disclosure relates to the field of energy storage converter control technologies, and in particular, to a method and an apparatus for controlling an energy storage voltage type converter.
Background
As an energy storage technology, one of the key technologies of energy revolution, the energy storage technology has received much attention in the industry in recent years because it can provide various auxiliary services such as peak shaving, frequency modulation, emergency response and the like for the power grid. In order to realize friendly grid connection of an energy storage system and provide stable voltage and frequency support for a power grid, research on a control strategy of an energy storage converter needs to be developed.
At present, in the field of energy storage converter control, double closed-loop control and dead-beat control are mostly adopted to realize dynamic response of voltage and frequency. But the conventional control strategy cannot maintain the stability of the asynchronous energy storage converter control system under the high permeability of the distributed power supply.
Disclosure of Invention
The present application is directed to solving, at least in part, one of the technical problems in the related art.
Therefore, a first objective of the present application is to provide a Control method for an energy storage voltage type converter, which solves the technical problem that the existing method cannot maintain the stability of an asynchronous energy storage converter Control system under high permeability of a distributed power supply, and adopts a two-period delay compensation Control strategy, i.e., a two-step MPC (Model Predictive Control) Control method, to perform advanced Control on system variables, accurately sample and offset delay effects, and adopts a two-step MPC frequency deviation power constraint function to improve Control accuracy, implement steady-state Control of the energy storage voltage type converter, and avoid frequency oscillation.
A second objective of the present application is to provide a control device for a storage voltage type inverter.
In order to achieve the above object, an embodiment of the first aspect of the present application provides a control method for a converter of a storage voltage type, including: acquiring a VSG (Virtual Synchronous Generation) rotor motion equation of the energy storage voltage type converter, and converting the VSG rotor motion equation into a Virtual angular frequency change rate mathematical model; converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation at a second preset moment; establishing a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation to minimize an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, and carrying out VSG power control according to the target active power reference value to realize MPC control of the energy storage voltage type converter.
Optionally, in an embodiment of the present application, obtaining a VSG rotor equation of motion of the energy storage voltage type converter includes:
acquiring an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type current converter;
and simulating the energy storage voltage type converter into a synchronous generator model according to the output active power mathematical model and the output reactive power mathematical model to obtain a VSG rotor motion equation of the energy storage voltage type converter.
Optionally, in an embodiment of the present application, before converting the virtual angular frequency change rate mathematical model into a virtual angular frequency adjustment MPC discrete domain mathematical equation according to the power MPC discrete domain matrix equation at the second preset time, the method further includes:
obtaining an active power change rate mathematical model and a reactive power change rate mathematical model of the energy storage voltage type current converter;
discretizing the active power change rate mathematical model and the reactive power change rate mathematical model to obtain an output active power MPC mathematical model, an output reactive power MPC mathematical model and an active and reactive power discrete domain matrix equation of the energy storage voltage type converter at a first preset moment;
and establishing a power MPC discrete domain matrix equation at a second preset moment according to the active power discrete domain matrix equation and the reactive power discrete domain matrix equation.
Optionally, in an embodiment of the present application, obtaining an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type converter includes:
constructing a current change rate equation of the energy storage voltage type current converter, and performing Clark conversion on the current change rate equation to obtain
Figure DEST_PATH_IMAGE001
A current change rate mathematical model under a coordinate system;
the energy storage voltage type current converter is obtained through a current change rate mathematical model to be static in two phases
Figure 312861DEST_PATH_IMAGE001
A voltage change rate mathematical model under a coordinate system;
and obtaining an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type current converter according to the current change rate mathematical model and the voltage change rate mathematical model.
Optionally, in an embodiment of the present application, obtaining an active power change rate mathematical model and a reactive power change rate mathematical model of the energy storage voltage type converter includes:
obtaining an instantaneous change rate mathematical model of the output power of the energy storage voltage type current converter according to the output active power mathematical model and the output reactive power mathematical model;
substituting the current change rate mathematical model and the voltage change rate mathematical model into the output power instantaneous change rate mathematical model to obtain the two-phase static state of the energy storage voltage type current converter
Figure 432126DEST_PATH_IMAGE001
An active power change rate mathematical model and a reactive power change rate mathematical model under a coordinate system.
Optionally, in an embodiment of the present application, the VSG rotor equation of motion of the energy storage voltage type converter is expressed as:
Figure 887379DEST_PATH_IMAGE002
wherein the content of the first and second substances,Jin order to be a virtual moment of inertia,
Figure DEST_PATH_IMAGE003
respectively VSG mechanical torque, electromagnetic torque and damping torque,
Figure 521622DEST_PATH_IMAGE004
in order to be the angular frequency of the frequency,
Figure DEST_PATH_IMAGE005
is the active power reference value and is,
Figure 87733DEST_PATH_IMAGE006
active power is output for the VSG, D is the damping coefficient,
Figure DEST_PATH_IMAGE007
in order to be the nominal angular frequency,
Figure 327084DEST_PATH_IMAGE008
virtual electrical angles for the VSG.
Optionally, in an embodiment of the present application, the discrete domain matrix equation of the power MPC at the second preset time is expressed as:
Figure DEST_PATH_IMAGE009
wherein the content of the first and second substances,
Figure 953238DEST_PATH_IMAGE010
is a matrix of coefficients, and is,
Figure DEST_PATH_IMAGE011
the active power prediction value is controlled for the two-step model,
Figure 74777DEST_PATH_IMAGE012
the reactive power prediction value is controlled for the two-step model,
Figure DEST_PATH_IMAGE013
in order to output the mathematical model of the active power MPC,
Figure 382262DEST_PATH_IMAGE014
in order to output a mathematical model of reactive power MPC,
Figure DEST_PATH_IMAGE015
is two phases at rest
Figure 538437DEST_PATH_IMAGE001
Under the coordinate system
Figure 335492DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure DEST_PATH_IMAGE017
is two phases static
Figure 944327DEST_PATH_IMAGE001
Under the coordinate system
Figure 789924DEST_PATH_IMAGE018
Shaft ac mains side voltage.
Optionally, in an embodiment of the present application, the virtual angular frequency rate of change mathematical model is represented as:
Figure DEST_PATH_IMAGE019
wherein the content of the first and second substances,
Figure 800605DEST_PATH_IMAGE020
in order to be the angular frequency of the frequency,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 768561DEST_PATH_IMAGE022
is the VSG output power variation;
the virtual angular frequency adjustment quantity MPC discrete domain mathematical equation is expressed as:
Figure DEST_PATH_IMAGE023
wherein, the first and the second end of the pipe are connected with each other,
Figure 864693DEST_PATH_IMAGE024
the predicted angular frequency change amount at time k +2,
Figure DEST_PATH_IMAGE025
for the predicted angular frequency change at time k +1, A, B is a coefficient matrix,
Figure 513980DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 379168DEST_PATH_IMAGE022
is the VSG output power variation.
Optionally, in an embodiment of the present application, the two-step MPC frequency deviation power constraint function is expressed as:
Figure 518025DEST_PATH_IMAGE026
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE027
representing a two-step MPC frequency deviation power constraint function,
Figure 101453DEST_PATH_IMAGE028
representing the system frequency deviation weighting function at time k +2,
Figure DEST_PATH_IMAGE029
and representing the weight function of the active power output by the energy storage VSG at the k +2 moment.
Optionally, in an embodiment of the present application, the output active power MPC mathematical model and the output reactive power MPC mathematical model of the energy storage voltage type converter at the first preset time are expressed as:
Figure 554431DEST_PATH_IMAGE030
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE031
in order to output the mathematical model of the active power MPC,
Figure 274126DEST_PATH_IMAGE032
in order to output a mathematical model of reactive power MPC,
Figure DEST_PATH_IMAGE033
in order to sample the control period of the system,
Figure 583884DEST_PATH_IMAGE034
is a filter inductor of an LC filter circuit,
Figure DEST_PATH_IMAGE035
is two phases static
Figure 592292DEST_PATH_IMAGE001
Under the coordinate system
Figure 645698DEST_PATH_IMAGE036
The shaft energy storage system outputs a voltage which,
Figure DEST_PATH_IMAGE037
is two phases at rest
Figure 485478DEST_PATH_IMAGE001
Under the coordinate system
Figure 966138DEST_PATH_IMAGE018
The shaft energy storage system outputs a voltage which,
Figure 461842DEST_PATH_IMAGE038
is time k
Figure 53360DEST_PATH_IMAGE016
The square of the voltage on the ac mains side of the shaft,
Figure DEST_PATH_IMAGE039
is time k
Figure 747647DEST_PATH_IMAGE018
The square of the voltage on the ac mains side of the shaft,
Figure 664787DEST_PATH_IMAGE040
is a filter inductor of an LC filter circuit,
Figure DEST_PATH_IMAGE041
the active power is output for the energy storage at the moment k,
Figure 444524DEST_PATH_IMAGE042
the reactive power is output for the energy storage at the moment k,
Figure DEST_PATH_IMAGE043
Figure 777416DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure DEST_PATH_IMAGE045
In that
Figure 591789DEST_PATH_IMAGE016
A shaft,
Figure 414251DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 681284DEST_PATH_IMAGE046
is the angular frequency;
the active and reactive power discrete domain matrix equation of the energy storage voltage type converter at the first preset moment is expressed as follows:
Figure DEST_PATH_IMAGE047
wherein the content of the first and second substances,
Figure 817868DEST_PATH_IMAGE010
in the form of a matrix of coefficients,
Figure 486746DEST_PATH_IMAGE048
in order to output the mathematical model of the active power MPC,
Figure DEST_PATH_IMAGE049
in order to output a mathematical model of reactive power MPC,
Figure 480110DEST_PATH_IMAGE050
the active power is output for the energy storage at the moment k,
Figure DEST_PATH_IMAGE051
the reactive power is output for the energy storage at the moment k,
Figure 234440DEST_PATH_IMAGE052
is time k
Figure 909135DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure DEST_PATH_IMAGE053
is time k
Figure 432520DEST_PATH_IMAGE018
Shaft ac mains side voltage.
In order to achieve the above object, a second aspect of the present application provides a control device for an inverter of a storage voltage type, including:
the acquisition module is used for acquiring a VSG rotor motion equation of the energy storage voltage type converter and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model;
the conversion module is used for converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to the power MPC discrete domain matrix equation at the second preset moment;
and the control module is used for establishing a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation, minimizing an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, carrying out VSG power control according to the target active power reference value, and realizing MPC control of the energy storage voltage type converter.
The control method and the device for the energy storage voltage type current converter solve the technical problem that the existing method cannot maintain the stability of an asynchronous energy storage current converter control system under the high permeability of a distributed power supply, a two-period delay compensation control strategy, namely a two-step MPC control method, is adopted to carry out advanced control on system variables, accurately sample and offset delay influence, and a two-step MPC frequency deviation power constraint function is adopted to improve the control precision, realize the steady-state control of the energy storage voltage type current converter and avoid frequency oscillation.
Additional aspects and advantages of the present application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the present application.
Drawings
The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
fig. 1 is a flowchart of a control method of an energy storage voltage type converter according to an embodiment of the present application;
fig. 2 is a topology diagram of a storage PCS circuit of a control method of a storage voltage type converter according to an embodiment of the present application;
fig. 3 is a schematic structural diagram of a control device of an energy storage voltage type inverter according to a second embodiment of the present application.
Detailed Description
Reference will now be made in detail to embodiments of the present application, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are exemplary and intended to be used for explaining the present application and should not be construed as limiting the present application.
The following describes a control method and apparatus for a storage voltage type converter according to an embodiment of the present application with reference to the drawings.
Fig. 1 is a flowchart of a control method of an energy storage voltage type converter according to an embodiment of the present disclosure.
As shown in fig. 1, the control method of the energy storage voltage type converter includes the following steps:
step 101, obtaining a VSG rotor motion equation of the energy storage voltage type converter, and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model.
In the embodiment of the application, an energy storage PCS current change rate equation is constructed according to kirchhoff's voltage law, wherein the energy storage PCS current change rate equation is expressed as:
Figure 862364DEST_PATH_IMAGE054
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE055
and L represents the equivalent inductance of the line,
Figure 103990DEST_PATH_IMAGE056
for energy storage of the PCS alternating current three-phase current,
Figure DEST_PATH_IMAGE057
for storing energy PCS alternating current three-phase voltage,
Figure 316796DEST_PATH_IMAGE058
Rthe equivalent resistance of the line is shown,
Figure 960267DEST_PATH_IMAGE045
is the three-phase voltage of an alternating current power grid,
Figure 295434DEST_PATH_IMAGE040
Figure 24355DEST_PATH_IMAGE034
Figure DEST_PATH_IMAGE059
the LC filter circuit is formed by the following steps,
Figure 34993DEST_PATH_IMAGE060
Figure DEST_PATH_IMAGE061
is an equivalent load.
Clark conversion is carried out on the current change rate equation of the energy storage PCS, and the conversion is carried out
Figure 532971DEST_PATH_IMAGE001
A mathematical model of rate of change of current in a coordinate system, wherein the mathematical model of rate of change of current is expressed as:
Figure 304618DEST_PATH_IMAGE062
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE063
Figure 255256DEST_PATH_IMAGE064
outputting current for energy storage system
Figure DEST_PATH_IMAGE065
In that
Figure 75445DEST_PATH_IMAGE016
A shaft,
Figure 427929DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 370477DEST_PATH_IMAGE066
Figure DEST_PATH_IMAGE067
for outputting voltage to energy storage system
Figure 808411DEST_PATH_IMAGE057
In that
Figure 229028DEST_PATH_IMAGE016
A shaft,
Figure 639281DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 752731DEST_PATH_IMAGE043
Figure 677961DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 902269DEST_PATH_IMAGE045
In that
Figure 963766DEST_PATH_IMAGE016
A shaft,
Figure 248117DEST_PATH_IMAGE018
An axial component.
Current change rate through energy storage PCS
Figure 598327DEST_PATH_IMAGE001
Energy storage PCS (Power System) obtained by current change rate mathematical model under coordinate system is static in two phases
Figure 626326DEST_PATH_IMAGE001
A mathematical model of rate of change of voltage under a coordinate system, wherein the mathematical model of rate of change of voltage is represented as:
Figure 807908DEST_PATH_IMAGE068
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE069
Figure 997581DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 897404DEST_PATH_IMAGE070
In that
Figure 463515DEST_PATH_IMAGE016
A shaft,
Figure 702866DEST_PATH_IMAGE018
The axial component, E is the net side voltage amplitude,
Figure 63440DEST_PATH_IMAGE071
is the angular frequency.
According to the stored energy PCS in
Figure 450559DEST_PATH_IMAGE001
Current change rate mathematical model and energy storage PCS (process control system) in coordinate system are static in two phases
Figure 820361DEST_PATH_IMAGE001
The voltage change rate mathematical model under the coordinate system obtains an energy storage PCS output active power and reactive power mathematical model, wherein the output active power and reactive power mathematical model is expressed as follows:
Figure 710957DEST_PATH_IMAGE072
wherein, P is an active power mathematical model, Q is a reactive power mathematical model,
Figure 508011DEST_PATH_IMAGE073
Figure 54530DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 228023DEST_PATH_IMAGE070
In that
Figure 238704DEST_PATH_IMAGE016
A shaft,
Figure 206660DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 302792DEST_PATH_IMAGE063
Figure 14396DEST_PATH_IMAGE074
outputting current for energy storage system
Figure 817267DEST_PATH_IMAGE065
In that
Figure 956124DEST_PATH_IMAGE016
A shaft,
Figure 539552DEST_PATH_IMAGE018
An axial component.
Outputting an active power and reactive power mathematical model according to the energy storage PCS, and simulating the energy storage PCS into a synchronous generator model to obtain a VSG rotor motion equation and a reactive power regulation equation, wherein the VSG rotor motion equation is expressed as:
Figure 789268DEST_PATH_IMAGE075
wherein the content of the first and second substances,Jin order to be a virtual moment of inertia,
Figure 774542DEST_PATH_IMAGE003
respectively VSG mechanical torque, electromagnetic torque and damping torque,
Figure 84300DEST_PATH_IMAGE004
in order to be the angular frequency of the frequency,
Figure 92707DEST_PATH_IMAGE076
is the active power reference value and is,
Figure 146114DEST_PATH_IMAGE077
active power is output for the VSG, D is the damping coefficient,
Figure 985894DEST_PATH_IMAGE007
in order to be the nominal angular frequency,
Figure 466554DEST_PATH_IMAGE078
virtual electrical angles for the VSG.
In the VSG control system, the virtual rotational inertia J enables the energy storage PCS to have inertia and damping coefficient in the power and frequency adjusting processDSo that the energy storage PCS has the capability of suppressing grid power oscillations. The VSG control also has an excitation regulation inertia, and the reactive regulation equation is expressed as:
Figure 758995DEST_PATH_IMAGE079
wherein the content of the first and second substances,uis a virtual internal potential of the VSG,
Figure 616093DEST_PATH_IMAGE080
is the effective value of the rated voltage,
Figure 310379DEST_PATH_IMAGE081
for the deviation between the virtual internal potential and the rated voltage,
Figure 899624DEST_PATH_IMAGE082
in order to obtain the reactive power regulation coefficient,
Figure 944940DEST_PATH_IMAGE083
the reactive power is output for the VSG,
Figure 340149DEST_PATH_IMAGE084
is a reactive power reference value.
And 102, converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to the power MPC discrete domain matrix equation at the second preset time.
In the embodiment of the application, energy storage PCS outputs active power and a reactive power mathematical model to be derived with respect to time, and an energy storage PCS output power instantaneous change rate mathematical model is obtained, wherein the output power instantaneous change rate mathematical model is as follows:
Figure 888942DEST_PATH_IMAGE085
wherein P is an active power mathematical model, Q is a reactive power mathematical model,
Figure 976984DEST_PATH_IMAGE043
Figure 244017DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 380600DEST_PATH_IMAGE045
In that
Figure 783900DEST_PATH_IMAGE016
A shaft,
Figure 42843DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 797172DEST_PATH_IMAGE063
Figure 534184DEST_PATH_IMAGE064
outputting current for energy storage system
Figure 57569DEST_PATH_IMAGE065
In that
Figure 487414DEST_PATH_IMAGE016
A shaft,
Figure 401143DEST_PATH_IMAGE018
An axial component.
The stored energy is PCS in
Figure 941846DEST_PATH_IMAGE001
Current change rate mathematical model and energy storage PCS (process control system) in coordinate system are static in two phases
Figure 585317DEST_PATH_IMAGE001
Substituting the voltage change rate mathematical model under the coordinate system into the energy storage PCS output power instantaneous change rate mathematical model to obtain the energy storage PCS which is static in two phases
Figure 920483DEST_PATH_IMAGE001
The mathematical models of the change rates of the active power and the reactive power under the coordinate system are expressed as follows:
Figure 649405DEST_PATH_IMAGE086
wherein P is an active power mathematical model, Q is a reactive power mathematical model,
Figure 728219DEST_PATH_IMAGE035
is two phases at rest
Figure 163880DEST_PATH_IMAGE001
Under the coordinate system
Figure 669948DEST_PATH_IMAGE036
The shaft energy storage system outputs a voltage which,
Figure 886165DEST_PATH_IMAGE037
is two phases at rest
Figure 768671DEST_PATH_IMAGE001
Under the coordinate system
Figure 121154DEST_PATH_IMAGE018
The shaft energy storage system outputs a voltage which,
Figure 63703DEST_PATH_IMAGE040
is a filter inductor of an LC filter circuit,
Figure 439320DEST_PATH_IMAGE034
is a filter inductor of an LC filter circuit,
Figure 859937DEST_PATH_IMAGE043
Figure 332507DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 180377DEST_PATH_IMAGE045
In that
Figure 105608DEST_PATH_IMAGE016
A shaft,
Figure 329916DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 656992DEST_PATH_IMAGE046
is the angular frequency.
Make energy storage PCS at two-phase standstill
Figure 879026DEST_PATH_IMAGE001
Discretizing the mathematical model of the change rate of the active power and the reactive power under the coordinate system to obtain the energy storage PCS (process control system) output active power and reactive power MPC mathematical model and an active power and reactive power discrete domain matrix equation at a first preset moment, wherein the first preset moment can bekThe +1 time, k represents an arbitrary time, the first preset time is the next time of the k time,
kthe mathematical model of the energy storage PCS at the +1 moment for outputting the active power and the reactive power MPC is expressed as follows:
Figure 291553DEST_PATH_IMAGE030
wherein the content of the first and second substances,
Figure 53972DEST_PATH_IMAGE031
in order to output the mathematical model of the active power MPC,
Figure 235555DEST_PATH_IMAGE032
in order to output a mathematical model of reactive power MPC,
Figure 690807DEST_PATH_IMAGE033
in order to sample the control period of the device,
Figure 590630DEST_PATH_IMAGE034
is a filter inductor of an LC filter circuit,
Figure 94424DEST_PATH_IMAGE035
is two phases at rest
Figure 130513DEST_PATH_IMAGE001
Under the coordinate system
Figure 756666DEST_PATH_IMAGE036
The shaft energy storage system outputs a voltage which,
Figure 143785DEST_PATH_IMAGE037
is two phases at rest
Figure 248008DEST_PATH_IMAGE001
Under the coordinate system
Figure 404182DEST_PATH_IMAGE018
The shaft energy storage system outputs a voltage which,
Figure 138920DEST_PATH_IMAGE038
is time k
Figure 747756DEST_PATH_IMAGE016
The square of the voltage on the ac mains side of the shaft,
Figure 921248DEST_PATH_IMAGE039
is time k
Figure 931930DEST_PATH_IMAGE018
The square of the voltage on the ac mains side of the shaft,
Figure 899886DEST_PATH_IMAGE040
is a filter inductor of an LC filter circuit,
Figure 996018DEST_PATH_IMAGE041
the active power is output for the energy storage at the moment k,
Figure 442043DEST_PATH_IMAGE042
the reactive power is output for the energy storage at the moment k,
Figure 504633DEST_PATH_IMAGE043
Figure 643491DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 226919DEST_PATH_IMAGE045
In that
Figure 476635DEST_PATH_IMAGE016
A shaft,
Figure 461908DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 771667DEST_PATH_IMAGE046
in order to be the angular frequency of the frequency,
according to a discrete domain mathematical model standard form equation:
Figure 780074DEST_PATH_IMAGE087
kthe energy storage PCS active power and reactive power discrete domain matrix equation at the +1 moment is expressed as follows:
Figure 833481DEST_PATH_IMAGE047
wherein the content of the first and second substances,
Figure 673261DEST_PATH_IMAGE010
in the form of a matrix of coefficients,
Figure 153921DEST_PATH_IMAGE048
in order to output the mathematical model of the active power MPC,
Figure 446362DEST_PATH_IMAGE049
in order to output a mathematical model of reactive power MPC,
Figure 303459DEST_PATH_IMAGE050
the active power is output for the energy storage at the moment k,
Figure 935429DEST_PATH_IMAGE051
the reactive power is output for the energy storage at the moment k,
Figure 586990DEST_PATH_IMAGE052
is time k
Figure 632306DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure 27516DEST_PATH_IMAGE053
is time k
Figure 576309DEST_PATH_IMAGE018
The shaft is supplied with a voltage on the ac mains side,
Figure 664350DEST_PATH_IMAGE010
expressed as:
Figure 869067DEST_PATH_IMAGE088
Figure 802388DEST_PATH_IMAGE089
wherein the content of the first and second substances,
Figure 471266DEST_PATH_IMAGE033
in order to sample the control period of the device,
Figure 730209DEST_PATH_IMAGE035
is two phases at rest
Figure 218960DEST_PATH_IMAGE001
Under the coordinate system
Figure 221551DEST_PATH_IMAGE036
The shaft energy storage system outputs a voltage which,
Figure 744936DEST_PATH_IMAGE037
is two phases at rest
Figure 846884DEST_PATH_IMAGE001
Under the coordinate system
Figure 88510DEST_PATH_IMAGE018
The shaft energy storage system outputs a voltage which,
Figure 629212DEST_PATH_IMAGE040
is a filter inductor of an LC filter circuit,
Figure 7104DEST_PATH_IMAGE034
is a filter inductor of an LC filter circuit,
Figure 607850DEST_PATH_IMAGE052
is time k
Figure 71192DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure 353269DEST_PATH_IMAGE053
is time k
Figure 851246DEST_PATH_IMAGE018
The shaft is supplied with a voltage on the ac mains side,
Figure 357314DEST_PATH_IMAGE046
is the angular frequency.
Due to the inherent cycle delay phenomenon of the energy storage PCS control system in the sampling and calculating links, the three-phase voltage of the alternating current power grid at the k moment in the MPC control link is caused
Figure 573532DEST_PATH_IMAGE090
The sampling value cannot be applied to the sampling period, and as errors are accumulated continuously, larger deviation of the control system is caused. In order to suppress control deviation caused by periodic delay, a two-period delay compensation control strategy, namely a two-step MPC control method, is adopted to carry out advanced control on system variables, accurately sample and offset delay influence, and improve control precision.
According to the energy storage PCS at the first preset moment, outputting an active power discrete domain matrix equation and a reactive power discrete domain matrix equation, establishing an energy storage PCS power MPC discrete domain matrix equation at the second preset moment, wherein the first preset moment can bekThe +1 time, k represents any time, the first preset time is the next time of the k time, and the second preset time may bek+2, k represents an arbitrary time, the second predetermined time beingkThe next time instant of the time instant +1,
kthe +2 moment energy storage PCS power MPC discrete domain matrix equation is expressed as:
Figure 456037DEST_PATH_IMAGE009
wherein the content of the first and second substances,
Figure 808521DEST_PATH_IMAGE010
in the form of a matrix of coefficients,
Figure 423173DEST_PATH_IMAGE011
the active power prediction value is controlled for the two-step model,
Figure 861108DEST_PATH_IMAGE012
the reactive power prediction value is controlled for the two-step model,
Figure 547304DEST_PATH_IMAGE091
in order to output the mathematical model of the active power MPC,
Figure 754294DEST_PATH_IMAGE014
in order to output a mathematical model of reactive power MPC,
Figure 867744DEST_PATH_IMAGE015
is two phases at rest
Figure 792975DEST_PATH_IMAGE001
Under the coordinate system
Figure 954966DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure 282042DEST_PATH_IMAGE017
is two phases at rest
Figure 566393DEST_PATH_IMAGE001
Under the coordinate system
Figure 713340DEST_PATH_IMAGE018
Shaft ac mains side voltage.
According to the method, the MPC control is used as an upper-level control strategy of an energy storage PCS control link, the MPC control output value is used as an energy storage PCS to output active and reactive power reference values, VSG power is corrected in real time, the stability of a control system is improved, and frequency oscillation is avoided. To suppress frequency oscillations, the present application employs a two-step MPC frequency deviation power constraint function.
The MPC is used as a front-end control strategy of the VSG, the output of the MPC is the input of the VSG, and the MPC is used for correcting the active power reference value of the VSG, so that the effect of improving the control stability of the VSG is achieved.
Converting a VSG rotor motion equation of the energy storage voltage type converter into a virtual angular frequency change rate mathematical model, wherein the virtual angular frequency change rate mathematical model is expressed as:
Figure 741339DEST_PATH_IMAGE019
wherein the content of the first and second substances,
Figure 922922DEST_PATH_IMAGE020
in order to be the angular frequency of the frequency,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure 378174DEST_PATH_IMAGE092
Figure 215680DEST_PATH_IMAGE093
as a virtual amount of angular frequency adjustment,
Figure 781790DEST_PATH_IMAGE094
in order to be at the nominal angular frequency,
Figure 817879DEST_PATH_IMAGE095
Figure 444033DEST_PATH_IMAGE022
for the amount of VSG output power variation,
Figure 565573DEST_PATH_IMAGE096
is the active power reference value and is,
Figure 935374DEST_PATH_IMAGE077
active power is output for the VSG.
Converting the virtual angular frequency change rate mathematical model into a virtual angular frequency adjustment MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation of the energy storage PCS at a second preset time, wherein the second preset time may bekAt the time +2, the time point,
the virtual angular frequency adjustment MPC discrete domain mathematical equation is expressed as:
Figure 29232DEST_PATH_IMAGE023
wherein the content of the first and second substances,
Figure 826287DEST_PATH_IMAGE024
the predicted angular frequency change amount at time k +2,
Figure 435123DEST_PATH_IMAGE025
for the predicted angular frequency change at time k +1, A, B is a coefficient matrix,
Figure 343036DEST_PATH_IMAGE097
Figure 353717DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 321673DEST_PATH_IMAGE098
in order to be the nominal angular frequency,
Figure 417805DEST_PATH_IMAGE099
Figure 67092DEST_PATH_IMAGE022
for the amount of VSG output power variation,
Figure 197859DEST_PATH_IMAGE100
is the active power reference value and is,
Figure 336717DEST_PATH_IMAGE077
the active power is output for the VSG,
wherein A, B is represented as:
Figure 654565DEST_PATH_IMAGE101
Figure 169860DEST_PATH_IMAGE102
wherein the content of the first and second substances,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure 155134DEST_PATH_IMAGE033
in order to sample the time for the system,
Figure 402576DEST_PATH_IMAGE103
is a time constant.
103, establishing a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation to minimize an MPC frequency deviation power constraint function value, performing vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, and performing VSG power control according to the target active power reference value to realize MPC control of the energy storage voltage type converter.
In order to correct the VSG power in real time, improve the control system stability, and avoid frequency oscillation, it is necessary to minimize the MPC frequency deviation power constraint function value. In the energy storage PCS control system, a two-step power MPC control strategy is combined with VSG control to form a closed-loop control system. Wherein the output of the MPC control is compared with the VSG active power reference value
Figure 207721DEST_PATH_IMAGE096
Vector addition is carried out to obtain a new active power reference value, the new active power reference value participates in VSG power control, and VSG outputs active power
Figure 261127DEST_PATH_IMAGE077
Reactive power
Figure 100907DEST_PATH_IMAGE104
Sum net side angular frequency
Figure 581567DEST_PATH_IMAGE071
Is an input to the MPC. The VSG power reference value is continuously corrected through the two-step MPC frequency deviation power constraint function, when the network side frequency rises, the MPC active output is negative, and the VSG active power reference value is enabled to be
Figure 139587DEST_PATH_IMAGE096
Reduce, thereby reducing VSG output power
Figure 668789DEST_PATH_IMAGE077
And further, the rise of the network side frequency is suppressed. When the frequency of the network side is reduced, the active output of the MPC is positive, so that the VSG active power reference value is enabled
Figure 628655DEST_PATH_IMAGE105
Increase to increase VSG output power
Figure 280216DEST_PATH_IMAGE077
And further, the network side frequency drop is suppressed.
According to the control method of the energy storage voltage type converter, the VSG rotor motion equation of the energy storage voltage type converter is obtained, and the VSG rotor motion equation is converted into a virtual angular frequency change rate mathematical model; converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation at a second preset moment; establishing a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation to minimize an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, and carrying out VSG power control according to the target active power reference value to realize MPC control of the energy storage voltage type converter. Therefore, the technical problem that the existing method cannot maintain the stability of the control system of the asynchronous energy storage converter under the high permeability of the distributed power supply can be solved, a two-period delay compensation control strategy, namely a two-step MPC control method, is adopted to carry out advanced control on system variables, accurately sample and offset delay influence, and a two-step MPC frequency deviation power constraint function is adopted to improve the control precision, realize the steady-state control of the energy storage voltage type converter and avoid frequency oscillation.
Further, in the embodiment of the present application, obtaining the VSG rotor motion equation of the energy storage voltage type converter includes:
acquiring an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type current converter;
and simulating the energy storage voltage type converter into a synchronous generator model according to the output active power mathematical model and the output reactive power mathematical model to obtain a VSG rotor motion equation of the energy storage voltage type converter.
Further, in this embodiment of the present application, before converting the virtual angular frequency change rate mathematical model into the virtual angular frequency adjustment MPC discrete domain mathematical equation according to the power MPC discrete domain matrix equation at the second preset time, the method further includes:
obtaining an active power change rate mathematical model and a reactive power change rate mathematical model of the energy storage voltage type current converter;
discretizing the active power change rate mathematical model and the reactive power change rate mathematical model to obtain an output active power MPC mathematical model, an output reactive power MPC mathematical model and an active and reactive power discrete domain matrix equation of the energy storage voltage type current converter at the first preset moment;
and establishing a power MPC discrete domain matrix equation at a second preset moment according to the active power discrete domain matrix equation and the reactive power discrete domain matrix equation.
Further, in this embodiment of the present application, obtaining an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type converter includes:
constructing a current change rate equation of the energy storage voltage type current converter, and performing Clark conversion on the current change rate equation to obtain
Figure 59953DEST_PATH_IMAGE001
A current change rate mathematical model under a coordinate system;
the energy storage voltage type current converter is obtained by a mathematical model of the current change rate in a two-phase static state
Figure 455162DEST_PATH_IMAGE001
A voltage change rate mathematical model under a coordinate system;
and obtaining an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type current converter according to the current change rate mathematical model and the voltage change rate mathematical model.
Further, in this embodiment of the present application, obtaining an active power change rate mathematical model and a reactive power change rate mathematical model of the energy storage voltage type converter includes:
obtaining an output power instantaneous change rate mathematical model of the energy storage voltage type current converter according to the output active power mathematical model and the output reactive power mathematical model;
substituting the current change rate mathematical model and the voltage change rate mathematical model into the output power instantaneous change rate mathematical model to obtain the energy storage voltage type current converter between twoPhase stationary
Figure 269535DEST_PATH_IMAGE001
An active power change rate mathematical model and a reactive power change rate mathematical model under a coordinate system.
Further, in the embodiment of the present application, the VSG rotor motion equation of the energy storage voltage type converter is expressed as:
Figure 91997DEST_PATH_IMAGE002
wherein the content of the first and second substances,Jin order to be a virtual moment of inertia,
Figure 296713DEST_PATH_IMAGE003
respectively VSG mechanical torque, electromagnetic torque and damping torque,
Figure 495614DEST_PATH_IMAGE004
in order to be the angular frequency of the frequency,
Figure 164492DEST_PATH_IMAGE005
is the active power reference value and is,
Figure 157856DEST_PATH_IMAGE006
active power is output for the VSG, D is the damping coefficient,
Figure 912185DEST_PATH_IMAGE007
in order to be the nominal angular frequency,
Figure 914777DEST_PATH_IMAGE008
is the VSG virtual electrical angle.
Further, in the embodiment of the present application, the discrete domain matrix equation of the power MPC at the second preset time is expressed as:
Figure 375845DEST_PATH_IMAGE009
wherein the content of the first and second substances,
Figure 540110DEST_PATH_IMAGE010
in the form of a matrix of coefficients,
Figure 781735DEST_PATH_IMAGE011
the active power prediction value is controlled for the two-step model,
Figure 322438DEST_PATH_IMAGE012
the reactive power predicted value is controlled for the two-step model,
Figure 700330DEST_PATH_IMAGE013
in order to output the mathematical model of the active power MPC,
Figure 301076DEST_PATH_IMAGE014
in order to output a mathematical model of reactive power MPC,
Figure 764418DEST_PATH_IMAGE015
is two phases static
Figure 46495DEST_PATH_IMAGE001
Under the coordinate system
Figure 278893DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure 50540DEST_PATH_IMAGE017
is two phases at rest
Figure 1178DEST_PATH_IMAGE001
Under the coordinate system
Figure 883684DEST_PATH_IMAGE018
Shaft ac mains side voltage.
Further, in the embodiment of the present application, the virtual angular frequency change rate mathematical model is expressed as:
Figure 236168DEST_PATH_IMAGE019
wherein the content of the first and second substances,
Figure 110540DEST_PATH_IMAGE004
in order to be the angular frequency of the frequency,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure 548474DEST_PATH_IMAGE106
Figure 234670DEST_PATH_IMAGE093
as a virtual amount of angular frequency adjustment,
Figure 441661DEST_PATH_IMAGE094
in order to be the nominal angular frequency,
Figure 555110DEST_PATH_IMAGE095
Figure 480341DEST_PATH_IMAGE022
for the amount of change in the VSG output power,
Figure 704649DEST_PATH_IMAGE096
is the active power reference value and is,
Figure 969408DEST_PATH_IMAGE077
outputting active power for the VSG;
the virtual angular frequency regulating quantity MPC discrete domain mathematical equation is expressed as follows:
Figure 988180DEST_PATH_IMAGE023
wherein, the first and the second end of the pipe are connected with each other,
Figure 400707DEST_PATH_IMAGE024
the predicted angular frequency change amount at time k +2,
Figure 428705DEST_PATH_IMAGE025
for the predicted angular frequency change at time k +1, A, B is a coefficient matrix,
Figure 610288DEST_PATH_IMAGE097
Figure 65540DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 637467DEST_PATH_IMAGE107
in order to be the nominal angular frequency,
Figure 203578DEST_PATH_IMAGE108
Figure 505246DEST_PATH_IMAGE022
for the amount of VSG output power variation,
Figure 131399DEST_PATH_IMAGE109
is the active power reference value and is,
Figure 252939DEST_PATH_IMAGE077
the active power is output for the VSG,
wherein A, B is represented as:
Figure 622741DEST_PATH_IMAGE101
Figure 716598DEST_PATH_IMAGE110
wherein the content of the first and second substances,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure 513653DEST_PATH_IMAGE033
in order to sample the time for the system,
Figure 122489DEST_PATH_IMAGE111
is a time constant.
The method comprises the steps of taking an energy storage PCS power MPC discrete domain matrix equation as an upper-layer control system of a VSG rotor motion equation and a reactive power regulation equation of an energy storage converter, taking an MPC control output value as a VSG active power reference value and a VSG reactive power reference value, correcting VSG power in real time, improving the stability of the control system, and avoiding frequency oscillation so as to inhibit the frequency oscillation.
Further, in the embodiment of the present application, the two-step MPC frequency deviation power constraint function is expressed as:
Figure 30402DEST_PATH_IMAGE026
wherein the content of the first and second substances,
Figure 41084DEST_PATH_IMAGE027
representing a two-step MPC frequency deviation power constraint function,
Figure 9040DEST_PATH_IMAGE028
representing the system frequency deviation weighting function at time k +2,
Figure 105172DEST_PATH_IMAGE029
and representing the weight function of the active power output by the energy storage VSG at the k +2 moment.
In order to obtain a frequency deviation weight function, an energy storage PCS power simulation synchronous machine control equation is established:
Figure 754459DEST_PATH_IMAGE112
wherein the content of the first and second substances,Jin order to be a virtual moment of inertia,
Figure 619646DEST_PATH_IMAGE100
is the active power reference value and is,
Figure 758504DEST_PATH_IMAGE077
active power is output for the energy storage PCS,Das damping coefficient, virtual moment of inertiaJThe energy storage PCS has inertia and damping coefficient in the process of power and frequency adjustmentDSo that the energy storage PCS has the capability of suppressing grid power oscillations,
according to the energy storage PCS power simulation synchronous machine control equation, establishing a frequency deviation weight function as follows:
Figure 341932DEST_PATH_IMAGE019
wherein the content of the first and second substances,
Figure 857227DEST_PATH_IMAGE020
in order to be the angular frequency of the frequency,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure 576921DEST_PATH_IMAGE092
Figure 824363DEST_PATH_IMAGE093
as a virtual amount of angular frequency adjustment,
Figure 895087DEST_PATH_IMAGE094
in order to be at the nominal angular frequency,
Figure 948494DEST_PATH_IMAGE095
Figure 788274DEST_PATH_IMAGE022
for the amount of change in the VSG output power,
Figure 268934DEST_PATH_IMAGE096
is the active power reference value and is,
Figure 826954DEST_PATH_IMAGE077
the active power is output for the VSG,
converting the frequency deviation weight function into a virtual angular frequency adjustment quantity MPC discrete domain mathematical equation according to a second preset moment energy storage PCS power MPC discrete domain matrix equation:
Figure 356155DEST_PATH_IMAGE023
wherein the content of the first and second substances,
Figure 50442DEST_PATH_IMAGE024
the predicted angular frequency change amount at time k +2,
Figure 967582DEST_PATH_IMAGE025
for the predicted angular frequency change at time k +1, A, B is a coefficient matrix,
Figure 747320DEST_PATH_IMAGE097
Figure 142529DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 956901DEST_PATH_IMAGE107
in order to be the nominal angular frequency,
Figure 779364DEST_PATH_IMAGE108
Figure 984080DEST_PATH_IMAGE022
for the amount of VSG output power variation,
Figure 182980DEST_PATH_IMAGE109
is the active power reference value and is,
Figure 851859DEST_PATH_IMAGE077
the active power is output for the VSG,
wherein A, B is represented as:
Figure 845223DEST_PATH_IMAGE101
Figure 599552DEST_PATH_IMAGE113
wherein the content of the first and second substances,Din order to be a damping coefficient of the damping,Jin order to be a virtual moment of inertia,
Figure 336564DEST_PATH_IMAGE114
is a system ofThe time of sampling is such that,
Figure 797632DEST_PATH_IMAGE103
is a time constant.
According to a virtual angular frequency adjustment quantity MPC discrete domain mathematical equation, a two-step MPC frequency deviation power constraint function can be expressed as follows:
Figure 227476DEST_PATH_IMAGE115
wherein the content of the first and second substances,
Figure 469102DEST_PATH_IMAGE116
is composed ofkThe time virtual angular frequency regulating quantity is the controlled output quantity,
Figure 744226DEST_PATH_IMAGE117
in order to control the amount of input in the domain,
Figure 387696DEST_PATH_IMAGE016
for the minimum prediction horizon for a two-step MPC,
Figure 722863DEST_PATH_IMAGE018
for the maximum prediction horizon of a two-step MPC,pin order to control the range of the control,
Figure 389468DEST_PATH_IMAGE118
are weighting factors used to reduce the amount of MPC control adjustment.
The problem of MPC optimization for virtual angular frequency adjustments with frequency deviation power constraint functions can be expressed as:
Figure 468282DEST_PATH_IMAGE119
wherein the content of the first and second substances,
Figure 966259DEST_PATH_IMAGE120
for the purpose of the frequency deviation power constraint function,
Figure 737906DEST_PATH_IMAGE121
the predicted angular frequency variation at time k +2+ i,
Figure 688545DEST_PATH_IMAGE122
the predicted angular frequency change at time k +1+ i, A, B is a coefficient matrix,
Figure 571050DEST_PATH_IMAGE092
Figure 923534DEST_PATH_IMAGE093
as a virtual amount of angular frequency adjustment,
Figure 803765DEST_PATH_IMAGE123
in order to be the nominal angular frequency,
Figure 241700DEST_PATH_IMAGE095
Figure 662317DEST_PATH_IMAGE022
for the amount of VSG output power variation,
Figure 134887DEST_PATH_IMAGE096
is the active power reference value and is,
Figure 248336DEST_PATH_IMAGE077
the active power is output for the VSG,
Figure 173567DEST_PATH_IMAGE124
is the angular frequency at time k + i,
Figure 335558DEST_PATH_IMAGE125
the angular frequency at time k + i-1,
Figure 397055DEST_PATH_IMAGE126
the predicted angular frequency variation at time k + i-1,
Figure 681406DEST_PATH_IMAGE127
predicted angular frequency variation at time k +2+ iThe minimum value of (a) is determined,
Figure 93932DEST_PATH_IMAGE128
the maximum value of the predicted angular frequency change amount at the time k +2+ i.
Further, in this embodiment of the present application, the mathematical model of the output active power MPC and the mathematical model of the output reactive power MPC of the energy storage voltage type converter at the first preset time are expressed as follows:
Figure 121931DEST_PATH_IMAGE030
wherein the content of the first and second substances,
Figure 303514DEST_PATH_IMAGE031
in order to output the mathematical model of the active power MPC,
Figure 430870DEST_PATH_IMAGE032
in order to output the mathematical model of the reactive power MPC,
Figure 330693DEST_PATH_IMAGE033
in order to sample the control period of the device,
Figure 896803DEST_PATH_IMAGE034
is a filter inductor of an LC filter circuit,
Figure 198472DEST_PATH_IMAGE035
is two phases at rest
Figure 559046DEST_PATH_IMAGE001
Under the coordinate system
Figure 946165DEST_PATH_IMAGE036
The shaft energy storage system outputs a voltage which,
Figure 315966DEST_PATH_IMAGE037
is two phases at rest
Figure 144245DEST_PATH_IMAGE001
Coordinate systemLower part
Figure 941300DEST_PATH_IMAGE018
The shaft energy storage system outputs a voltage which,
Figure 550136DEST_PATH_IMAGE038
is time k
Figure 723628DEST_PATH_IMAGE016
The square of the voltage on the ac mains side of the shaft,
Figure 734309DEST_PATH_IMAGE039
is time k
Figure 702265DEST_PATH_IMAGE018
The square of the voltage on the ac mains side of the shaft,
Figure 736080DEST_PATH_IMAGE040
is a filter inductor of an LC filter circuit,
Figure 447685DEST_PATH_IMAGE041
the active power is output for the energy storage at the moment k,
Figure 312872DEST_PATH_IMAGE042
the reactive power is output for the energy storage at the moment k,
Figure 451730DEST_PATH_IMAGE043
Figure 35158DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 284873DEST_PATH_IMAGE045
In that
Figure 207830DEST_PATH_IMAGE016
A shaft,
Figure 517589DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 588313DEST_PATH_IMAGE046
is the angular frequency;
the active and reactive power discrete domain matrix equation of the energy storage voltage type converter at the first preset moment is expressed as follows:
Figure 641720DEST_PATH_IMAGE047
wherein the content of the first and second substances,
Figure 481500DEST_PATH_IMAGE010
in the form of a matrix of coefficients,
Figure 962160DEST_PATH_IMAGE048
in order to output the mathematical model of the active power MPC,
Figure 254601DEST_PATH_IMAGE049
in order to output a mathematical model of reactive power MPC,
Figure 879459DEST_PATH_IMAGE050
the active power is output for the energy storage at the moment k,
Figure 573746DEST_PATH_IMAGE051
the reactive power is output for the energy storage at the moment k,
Figure 225307DEST_PATH_IMAGE052
is time k
Figure 270624DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure 665833DEST_PATH_IMAGE053
is time k
Figure 214626DEST_PATH_IMAGE018
Shaft ac mains side voltage.
Fig. 2 is a topology diagram of a storage PCS circuit of a method for controlling a storage voltage type converter according to an embodiment of the present application.
As shown in figure 2 of the drawings, in which,
Figure 240351DEST_PATH_IMAGE129
in order to store the dc side bus voltage of the PCS,
Figure 507384DEST_PATH_IMAGE130
Figure 706284DEST_PATH_IMAGE065
for energy storage PCS alternating current A, B, C three-phase voltage and current,
Figure 109584DEST_PATH_IMAGE090
is the three-phase voltage of an alternating current power grid,
Figure 368527DEST_PATH_IMAGE040
Figure 122856DEST_PATH_IMAGE034
Figure 797551DEST_PATH_IMAGE059
the LC filter circuit is formed by the following steps,
Figure 320936DEST_PATH_IMAGE060
Figure 750780DEST_PATH_IMAGE061
is an equivalent load.
Fig. 3 is a schematic structural diagram of a control device of an energy storage voltage type inverter according to a second embodiment of the present application.
As shown in fig. 3, the control apparatus for a converter of a storage voltage type includes:
the acquiring module 10 is used for acquiring a VSG rotor motion equation of the energy storage voltage type converter and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model;
the conversion module 20 is configured to convert the virtual angular frequency change rate mathematical model into a virtual angular frequency adjustment MPC discrete domain mathematical equation according to the power MPC discrete domain matrix equation at the second preset time;
and the control module 30 is configured to establish a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation, minimize an MPC frequency deviation power constraint function value, perform vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, perform VSG power control according to the target active power reference value, and implement MPC control of the energy storage voltage type converter.
The control device of the energy storage voltage type converter comprises an obtaining module, a calculating module and a control module, wherein the obtaining module is used for obtaining a VSG rotor motion equation of the energy storage voltage type converter and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model; the conversion module is used for converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to the power MPC discrete domain matrix equation at the second preset moment; and the control module is used for establishing a two-step MPC frequency deviation power constraint function according to a virtual angular frequency adjustment MPC discrete domain mathematical equation, minimizing an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, carrying out VSG power control according to the target active power reference value, and realizing MPC control of the energy storage voltage type converter. Therefore, the technical problem that the existing method cannot maintain the stability of the control system of the asynchronous energy storage converter under the high permeability of the distributed power supply can be solved, a two-period delay compensation control strategy, namely a two-step MPC control method, is adopted to carry out advanced control on system variables, accurately sample and offset delay influence, and a two-step MPC frequency deviation power constraint function is adopted to improve the control precision, realize the steady-state control of the energy storage voltage type converter and avoid frequency oscillation.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example" or "some examples" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the application. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Moreover, various embodiments or examples and features of various embodiments or examples described in this specification can be combined and combined by one skilled in the art without being mutually inconsistent.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present application, "plurality" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing steps of a custom logic function or process, and alternate implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present application.
The logic and/or steps represented in the flowcharts or otherwise described herein, e.g., an ordered listing of executable instructions that can be considered to implement logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium could even be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, via for instance optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.
It should be understood that portions of the present application may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the various steps or methods may be implemented in software or firmware stored in memory and executed by a suitable instruction execution system. If implemented in hardware, as in another embodiment, any one or combination of the following techniques, which are known in the art, may be used: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application specific integrated circuit having an appropriate combinational logic gate circuit, a Programmable Gate Array (PGA), a Field Programmable Gate Array (FPGA), or the like.
It will be understood by those skilled in the art that all or part of the steps carried by the method for implementing the above embodiments may be implemented by hardware related to instructions of a program, which may be stored in a computer readable storage medium, and when the program is executed, the program includes one or a combination of the steps of the method embodiments.
In addition, functional units in the embodiments of the present application may be integrated into one processing module, or each unit may exist alone physically, or two or more units are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode. The integrated module, if implemented in the form of a software functional module and sold or used as a stand-alone product, may also be stored in a computer readable storage medium.
The storage medium mentioned above may be a read-only memory, a magnetic or optical disk, etc. Although embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application, and that variations, modifications, substitutions and alterations may be made to the above embodiments by those of ordinary skill in the art within the scope of the present application.

Claims (11)

1. A method for controlling a converter of a stored energy voltage type, comprising:
acquiring a VSG rotor motion equation of an energy storage voltage type converter, and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model;
converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation at a second preset moment;
and establishing a two-step MPC frequency deviation power constraint function according to the virtual angular frequency adjustment value MPC discrete domain mathematical equation to minimize an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, and carrying out VSG power control according to the target active power reference value to realize MPC control of the energy storage voltage type converter.
2. The method for controlling the energy storage voltage type converter according to claim 1, wherein the obtaining of the VSG rotor motion equation of the energy storage voltage type converter comprises:
acquiring an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type current converter;
and simulating the energy storage voltage type converter into a synchronous generator model according to the output active power mathematical model and the output reactive power mathematical model to obtain a VSG rotor motion equation of the energy storage voltage type converter.
3. A method for controlling a converter according to claim 2, wherein before said converting said virtual angular frequency change rate mathematical model into a virtual angular frequency adjustment MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation at a second predetermined time, further comprising:
obtaining an active power change rate mathematical model and a reactive power change rate mathematical model of the energy storage voltage type current converter;
discretizing the active power change rate mathematical model and the reactive power change rate mathematical model to obtain an output active power MPC mathematical model, an output reactive power MPC mathematical model and an active and reactive power discrete domain matrix equation of the energy storage voltage type converter at the first preset moment;
and establishing a power MPC discrete domain matrix equation at a second preset moment according to the active power discrete domain matrix equation and the reactive power discrete domain matrix equation.
4. The method for controlling a storage voltage type converter according to claim 3, wherein the obtaining of the mathematical model of the output active power and the mathematical model of the output reactive power of the storage voltage type converter comprises:
constructing a current change rate equation of the energy storage voltage type current converter, and performing Clark conversion on the current change rate equation to obtain
Figure 126165DEST_PATH_IMAGE001
A current change rate mathematical model under a coordinate system;
obtaining the static state of the energy storage voltage type current converter in two phases through the current change rate mathematical model
Figure 288157DEST_PATH_IMAGE001
A voltage change rate mathematical model under a coordinate system;
and obtaining an output active power mathematical model and an output reactive power mathematical model of the energy storage voltage type current converter according to the current change rate mathematical model and the voltage change rate mathematical model.
5. The method for controlling the energy storage voltage type converter according to claim 4, wherein the step of obtaining the mathematical model of the active power change rate and the mathematical model of the reactive power change rate of the energy storage voltage type converter comprises:
obtaining an output power instantaneous change rate mathematical model of the energy storage voltage type current converter according to the output active power mathematical model and the output reactive power mathematical model;
substituting the current change rate mathematical model and the voltage change rate mathematical model into the output power instantaneous change rate mathematical model to obtain the condition that the energy storage voltage type current converter is static at two phases
Figure 615233DEST_PATH_IMAGE001
An active power change rate mathematical model and a reactive power change rate mathematical model under a coordinate system.
6. The method for controlling an energy storage voltage type converter according to claim 1, wherein the VSG rotor motion equation of the energy storage voltage type converter is expressed as:
Figure 899583DEST_PATH_IMAGE002
wherein J is a virtual moment of inertia,
Figure 46531DEST_PATH_IMAGE003
respectively VSG mechanical torque, electromagnetic torque and damping torque,
Figure 74530DEST_PATH_IMAGE004
in order to be the angular frequency of the frequency,
Figure 193796DEST_PATH_IMAGE005
is the active power reference value and is,
Figure 649048DEST_PATH_IMAGE006
active power is output for the VSG, D is the damping coefficient,
Figure 548871DEST_PATH_IMAGE007
in order to be at the nominal angular frequency,
Figure 114981DEST_PATH_IMAGE008
is the VSG virtual electrical angle.
7. A method for controlling a converter according to claim 1, wherein the discrete domain matrix equation of the power MPC at the second predetermined time is expressed as:
Figure 151070DEST_PATH_IMAGE009
wherein the content of the first and second substances,
Figure 777224DEST_PATH_IMAGE010
is a matrix of coefficients, and is,
Figure 898763DEST_PATH_IMAGE011
the active power prediction value is controlled for the two-step model,
Figure 206248DEST_PATH_IMAGE012
the reactive power prediction value is controlled for the two-step model,
Figure 362423DEST_PATH_IMAGE013
in order to output the mathematical model of the active power MPC,
Figure 159478DEST_PATH_IMAGE014
in order to output a mathematical model of reactive power MPC,
Figure 768313DEST_PATH_IMAGE015
is two phases at rest
Figure 676227DEST_PATH_IMAGE001
Under the coordinate system
Figure 686908DEST_PATH_IMAGE016
The voltage on the ac mains side of the shaft,
Figure 592547DEST_PATH_IMAGE017
is two phases static
Figure 688679DEST_PATH_IMAGE001
Under the coordinate system
Figure 400283DEST_PATH_IMAGE018
Shaft ac mains side voltage.
8. A method of controlling a converter according to claim 1 wherein the mathematical model of the rate of change of virtual angular frequency is expressed as:
Figure 531050DEST_PATH_IMAGE019
wherein the content of the first and second substances,
Figure 669907DEST_PATH_IMAGE020
angular frequency, D damping coefficient, J virtual moment of inertia,
Figure 987756DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 440734DEST_PATH_IMAGE022
is the VSG output power variation;
the virtual angular frequency regulating quantity MPC discrete domain mathematical equation is expressed as follows:
Figure 426008DEST_PATH_IMAGE023
wherein the content of the first and second substances,
Figure 735767DEST_PATH_IMAGE024
the predicted angular frequency change amount at time k +2,
Figure 540911DEST_PATH_IMAGE025
the predicted angular frequency change at time k +1, A, B is a coefficient matrix,
Figure 594318DEST_PATH_IMAGE021
as a virtual amount of angular frequency adjustment,
Figure 434098DEST_PATH_IMAGE022
is the VSG output power variation.
9. A method of controlling a converter according to claim 1 wherein the two-step MPC frequency deviation power constraint function is expressed as:
Figure 852441DEST_PATH_IMAGE026
wherein the content of the first and second substances,
Figure 410461DEST_PATH_IMAGE027
representing a two-step MPC frequency deviation power constraint function,
Figure 1980DEST_PATH_IMAGE028
representing the system frequency deviation weighting function at time k +2,
Figure 961846DEST_PATH_IMAGE029
to representAnd the energy storage VSG outputs an active power weight function at the k +2 moment.
10. A method for controlling a storage voltage type converter according to claim 3, characterized in that the mathematical models of the output active power MPC and the output reactive power MPC of the storage voltage type converter at the first predetermined moment are expressed as:
Figure 613407DEST_PATH_IMAGE030
wherein the content of the first and second substances,
Figure 393144DEST_PATH_IMAGE031
in order to output the mathematical model of the active power MPC,
Figure 788353DEST_PATH_IMAGE032
in order to output a mathematical model of reactive power MPC,
Figure 540409DEST_PATH_IMAGE033
in order to sample the control period of the device,
Figure 362871DEST_PATH_IMAGE034
is a filter inductor of an LC filter circuit,
Figure 629904DEST_PATH_IMAGE035
is two phases at rest
Figure 828804DEST_PATH_IMAGE001
Under the coordinate system
Figure 497683DEST_PATH_IMAGE036
The shaft energy storage system outputs a voltage which,
Figure 491047DEST_PATH_IMAGE037
is two phases at rest
Figure 183060DEST_PATH_IMAGE001
Under the coordinate system
Figure 185651DEST_PATH_IMAGE018
The shaft energy storage system outputs a voltage which,
Figure 709036DEST_PATH_IMAGE038
is time k
Figure 873301DEST_PATH_IMAGE016
The square of the voltage on the ac mains side of the shaft,
Figure 114926DEST_PATH_IMAGE039
is time k
Figure 655629DEST_PATH_IMAGE018
The square of the voltage on the ac mains side of the shaft,
Figure 33521DEST_PATH_IMAGE040
is a filter inductor of an LC filter circuit,
Figure 589528DEST_PATH_IMAGE041
the active power is output for the energy storage at the moment k,
Figure 52870DEST_PATH_IMAGE042
the reactive power is output for the energy storage at the moment k,
Figure 397264DEST_PATH_IMAGE043
Figure 629662DEST_PATH_IMAGE044
for three-phase voltage of AC mains
Figure 401309DEST_PATH_IMAGE045
In that
Figure 351947DEST_PATH_IMAGE016
A shaft,
Figure 172136DEST_PATH_IMAGE018
The axial component of the magnetic flux is,
Figure 524620DEST_PATH_IMAGE046
is the angular frequency;
the active and reactive power discrete domain matrix equation of the energy storage voltage type converter at the first preset moment is expressed as follows:
Figure 467168DEST_PATH_IMAGE047
wherein the content of the first and second substances,
Figure 905103DEST_PATH_IMAGE010
in the form of a matrix of coefficients,
Figure 591299DEST_PATH_IMAGE048
in order to output the mathematical model of the active power MPC,
Figure 798289DEST_PATH_IMAGE049
in order to output a mathematical model of reactive power MPC,
Figure 849422DEST_PATH_IMAGE050
the active power is output for the energy storage at the moment k,
Figure 774653DEST_PATH_IMAGE051
the reactive power is output for the energy storage at the moment k,
Figure 998960DEST_PATH_IMAGE052
is time k
Figure 326037DEST_PATH_IMAGE016
The shaft is supplied with a voltage on the ac mains side,
Figure 344808DEST_PATH_IMAGE053
is time k
Figure 757335DEST_PATH_IMAGE018
Shaft ac mains side voltage.
11. A control apparatus for a converter of the energy storage voltage type, comprising:
the acquisition module is used for acquiring a VSG rotor motion equation of the energy storage voltage type converter and converting the VSG rotor motion equation into a virtual angular frequency change rate mathematical model;
the conversion module is used for converting the virtual angular frequency change rate mathematical model into a virtual angular frequency regulating quantity MPC discrete domain mathematical equation according to a power MPC discrete domain matrix equation at a second preset time;
and the control module is used for establishing a two-step MPC frequency deviation power constraint function according to the virtual angular frequency adjustment value MPC discrete domain mathematical equation, minimizing an MPC frequency deviation power constraint function value, carrying out vector addition on an output value controlled by the MPC and a VSG active power reference value to obtain a target active power reference value, carrying out VSG power control according to the target active power reference value, and realizing MPC control of the energy storage voltage type converter.
CN202210944503.9A 2022-08-08 2022-08-08 Control method and device for energy storage voltage type current converter Withdrawn CN115065077A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210944503.9A CN115065077A (en) 2022-08-08 2022-08-08 Control method and device for energy storage voltage type current converter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210944503.9A CN115065077A (en) 2022-08-08 2022-08-08 Control method and device for energy storage voltage type current converter

Publications (1)

Publication Number Publication Date
CN115065077A true CN115065077A (en) 2022-09-16

Family

ID=83207771

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210944503.9A Withdrawn CN115065077A (en) 2022-08-08 2022-08-08 Control method and device for energy storage voltage type current converter

Country Status (1)

Country Link
CN (1) CN115065077A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116223900A (en) * 2022-12-30 2023-06-06 广州视骁科技有限公司 Method, device, electronic equipment and storage medium for determining alternating current circuit power

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113595140A (en) * 2021-07-30 2021-11-02 西安热工研究院有限公司 Method for establishing MPC weight value function of energy storage converter
CN113595116A (en) * 2021-07-30 2021-11-02 西安热工研究院有限公司 Method for establishing MPC discrete mathematical model of energy storage converter
CN113783459A (en) * 2021-09-15 2021-12-10 华能洋浦热电有限公司 Energy storage commutation device model prediction control method and system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113595140A (en) * 2021-07-30 2021-11-02 西安热工研究院有限公司 Method for establishing MPC weight value function of energy storage converter
CN113595116A (en) * 2021-07-30 2021-11-02 西安热工研究院有限公司 Method for establishing MPC discrete mathematical model of energy storage converter
CN113783459A (en) * 2021-09-15 2021-12-10 华能洋浦热电有限公司 Energy storage commutation device model prediction control method and system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
侯龙等: "采用新型MPC提高储能PCS虚拟同步机控制稳定性研究", 《国外电子测量技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116223900A (en) * 2022-12-30 2023-06-06 广州视骁科技有限公司 Method, device, electronic equipment and storage medium for determining alternating current circuit power
CN116223900B (en) * 2022-12-30 2024-01-09 广州视骁科技有限公司 Method, device, electronic equipment and storage medium for determining alternating current circuit power

Similar Documents

Publication Publication Date Title
CN115313524B (en) Photovoltaic power generation grid-connected control method and system based on network-building type converter
CN103683319A (en) Grid connected inverter control method based on hysteresis modulation under unbalanced grid voltage condition
CN110266056B (en) Virtual synchronous power generation control method based on light storage independent grid connection
CN114172183A (en) Multi-target cooperative control method and device for energy storage in power system
CN115065077A (en) Control method and device for energy storage voltage type current converter
CN109286200B (en) Control method and control system of variable-speed constant-frequency wind turbine generator
CN116845886B (en) Multi-port autonomous photovoltaic system network construction control method based on model prediction
CN110266044A (en) A kind of micro-grid connection control system and method based on energy accumulation current converter
Lajouad et al. Robust adaptive nonlinear controller of wind energy conversion system based on permanent magnet synchronous generator
CN114759575A (en) Virtual synchronous double-fed fan subsynchronous oscillation suppression method and system
CN116470581A (en) Method and system for controlling operation of voltage source type virtual synchronous machine
Belfedal et al. Comparison of PI and direct power control with SVM of doubly fed induction generator
CN115102241A (en) Control method and device for network-building type double-fed fan and computer readable storage medium
CN114977249A (en) Stability analysis method for energy storage converter control system
Shehata Direct power control of wind‐turbine‐driven DFIG during transient grid voltage unbalance
CN115065076A (en) Energy storage voltage type converter MPC and VSG coordination control method and device
Jiang et al. Secondary reactive power balancing and voltage stability in microgrid using prioritized centralized controller
CN112271759A (en) Method, system and device for optimizing wind power frequency modulation parameters under multi-frequency modulation resource framework
CN113013926A (en) Distributed grid-connected power generation system sequence impedance aggregation method and system
Du et al. A virtual synchronous generator adopting dynamic damping without frequency detection
Pourjafari et al. A new approach in sensor-less vector control of stand-alone doubly fed induction generator
Revelles et al. Comparison of current control structures for three-phase four-wire systems in natural frame
CN112968452B (en) Control method and device of converter
Khan Development of a power conditioner for a PMSG-based wind energy system integrated into a weak grid
Amran et al. Improvement performances of Doubly Fed Induction Generator via MPPT Strategy using Model Reference Adaptive Control based Direct Power Control with Space Vector Modulation

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
WW01 Invention patent application withdrawn after publication

Application publication date: 20220916

WW01 Invention patent application withdrawn after publication