Summary of the invention
The object of the present invention is to provide a kind of non-linear voltage control method and systems, can be to virtual synchronous generator
End voltage is controlled, to reach the oscillation for slowing down system in transient process, improves the purpose of dynamic response process.
To achieve the above object, the present invention provides following schemes:
A kind of non-linear voltage control method, which comprises
The calculation formula about quantity of state is determined according to the control system of virtual synchronous generator, obtains the control system
Nonlinear system model, the quantity of state includes the angular frequency and generator rotor angle of the virtual synchronous generator;
The nonlinear system model is converted into Linear system model by changes in coordinates;
The optimum control rate of the end voltage of the virtual synchronous generator in the Linear system model is solved,
Obtain the optimum control rate of linear system middle-end voltage;
In the optimum control rate generation of the linear system middle-end voltage, is returned in the nonlinear system, nonlinear system is obtained
Described in virtual synchronous generator end voltage optimum control rate;
Regulated and controled using the optimum control rate of virtual synchronous generator voltage described in the nonlinear system described virtual
The amplitude of synchronous motor end voltage, to slow down the oscillation of system in transient process.
Optionally, the calculation formula determined according to the control system of virtual synchronous generator about quantity of state, obtains
The nonlinear system model of the control system, specifically includes:
Turned according to the input of each controlling unit, the virtual synchronous generator in the control circuit of the control system
Electromagnetic torque in square, the control system determines the calculation formula of control process;
Obtain the electromagnetic torque, the electromagnetic torqueWherein, E is the virtual synchronous engine
End voltage, U is network voltage, XLFor line impedance, δ is the generator rotor angle of the virtual synchronous generator, and ω is virtual synchronous hair
Motor angular frequency;
To be brought into the calculation formula of the control process, obtain about the quantity of state formula wherein, the end voltage E of the virtual synchronous generator is control variable, and X is the quantity of state,
For the derivative of the quantity of state, X=[δ ω]T,
Optionally, described that the nonlinear system model is converted into Linear system model by changes in coordinates, it is specific to wrap
It includes:
It will be in nonlinear systemIt is transformed into linear system, is obtained by coordinate transformWherein, Z is the state variable under linear system,For the derivative of the state variable under linear system, A is shape
State variable coefficient matrix coefficient matrix, B are the coefficient matrix of corresponding control amount, and v is the optimum control of linear system end voltage
Rate,
Optionally, the optimum control of the end voltage of the virtual synchronous generator in the Linear system model
Rate is solved, and is obtained the optimum control rate of linear system middle-end voltage, is specifically included:
Obtain the performance indicator formula that linear system seeks optimum control rateWherein, Q
It is symmetrical positive definite weight matrix with R;
According to the formulaLinear optimal control is sought using linear optimal active control algorithm
Rate v processed.
A kind of non-linear voltage control system, the system comprises:
Quantity of state calculation formula determination unit, for being determined according to the control system of virtual synchronous generator about quantity of state
Calculation formula, obtain the nonlinear system model of the control system, the quantity of state includes the virtual synchronous generator
Angular frequency and generator rotor angle;
Linear system model conversion unit, for the nonlinear system model to be converted to linear system by changes in coordinates
System model;
Linear system optimum control rate computing unit, for carrying out the virtual synchronous power generation to the Linear system model
The solution of set end voltage optimum control rate obtains the optimum control rate of linear system middle-end voltage;
Nonlinear system optimal control rate computing unit, for by the optimum control rate generation of the linear system middle-end voltage
It returns in the nonlinear system, obtains the optimum control rate of virtual synchronous generator voltage described in nonlinear system;
Voltage control unit is held, for utilizing the optimal of virtual synchronous generator voltage described in the nonlinear system
Control rate regulates and controls the amplitude of the virtual synchronous motor terminal voltage, to slow down the oscillation of system in transient process.
Optionally, the quantity of state calculation formula determination unit, specifically includes:
First state amount calculation formula determines subelement, for each control in the control circuit according to the control system
Link, the input torque of the virtual synchronous generator, the electromagnetic torque in the control system determine the calculating of control process
Formula;
Electromagnetic torque obtains subelement, for obtaining the electromagnetic torque, the electromagnetic torqueIts
In, E is the end voltage of the virtual synchronous engine, and U is network voltage, XLFor line impedance, δ is virtual synchronous power generation
The generator rotor angle of machine, ω are virtual synchronous generator angular frequency;
Second quantity of state calculation formula determines subelement, and being used for willIt is brought into the meter of the control process
It calculates in formula, obtains the formula about the quantity of stateWherein, the end of the virtual synchronous generator
Voltage E is control variable, and X is the quantity of state,For the derivative of the quantity of state, X=[δ ω]T,
Optionally, the Linear system model conversion unit, specifically includes:
Linear system model transforming subunit, for by being become in nonlinear system by coordinate
It changes and is transformed into linear system, obtain wherein, Z is the state variable under linear system, is under linear system
State variable derivative, A be state variable coefficient matrix coefficient matrix, B be corresponding control amount coefficient matrix, v is line
The optimum control rate of property system end voltage,
Optionally, the linear system optimum control rate computing unit, specifically includes:
Performance indicator formula obtains subelement, and the performance indicator formula of optimum control rate is sought for obtaining linear systemWherein, Q and R is symmetrical positive definite weight matrix;
Linear system optimum control rate computation subunit, for according to the formulaBenefit
Linear optimal contro8 rate v is sought with linear optimal active control algorithm.
The specific embodiment provided according to the present invention, the invention discloses following technical effects: invention introduces frequency,
The parameters such as generator rotor angle, line impedance and rotary inertia regulate and control the end voltage of virtual synchronous generator, make the control for holding voltage
Not only and voltage deviation and idle deviation it is related, also accurately embody the influence of the physical quantity in other transient processes.Pass through
Nonlinear system is linearized, linear optimal contro8 rate is sought, its generation is returned in original system, obtain under nonlinear system most
Excellent control rate.To slow down the oscillation of system in transient process, dynamic performance is improved, the stabilization of system is improved
Property.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of non-linear voltage control method and systems, can be to virtual synchronous generator
End voltage is controlled, to reach the oscillation for slowing down system in transient process, improves the purpose of dynamic response process.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Fig. 1 is the Basic Topological of virtual synchronous generator, as shown in Figure 1, virtual synchronous generator includes three contraries
Become device, filter, the calculating of grid-connected port power, VSG control algolithm and SVPWM and modulates 5 submodules.Virtual synchronous generator
It (VSG) is there is inverter similar with synchronous generator by the mechanical property and electromagnetic property of imitation synchronous generator
Characteristic provides the purpose that inertia is supported and damping is supported to be reached for power grid.
As shown in Figure 1, e=[ea eb ec]T, u=[ua ub uc]T, i=[ia ib ic]T, respectively virtual synchronous power generation
Machine three-phase induction electromotive force, output end voltage and grid-connected current;RsAnd LsRespectively refer to virtual stator armature resistance and synchronous electricity
Sense;PeWith QeIt is the active power and reactive power of VSG output respectively.
As shown in Figure 1, virtual synchronous generator mainly includes main circuit and control system.Wherein, main circuit is conventional
Gird-connected inverter topology, including renewable energy source (can be considered prime mover), DC/AC converter and filter circuit etc. are (corresponding same
Walk the energy converting between mechanical process of generator);Control system is to realize the core of virtual synchronous generator, mainly includes virtual
Synchronous generator ontology model and control algolithm, the former mainly intends the electromagnetic relationship and machinery of synchronous generator from mechanism upper mold
Movement, the latter then mainly simulate the active frequency modulation and the features such as idle pressure regulation of synchronous generator from external characteristics.
Fig. 2 is traditional virtual synchronous generator voltage control schematic diagram, as shown in Fig. 2, in traditional VSG control,
Voltage is held to be codetermined by the idle and idle reality output deviation of the deviation and setting of setting voltage and practical end voltage.
As shown in Fig. 2, the end voltage expression of VSG isIt follows that
The end voltage of VSG is only influenced by idle deviation and voltage deviation.And in actual synchronous generator, end voltage is by excitation electricity
Stream determines, the size of exciting current in addition to being influenced by active reactive, also with frequency, generator rotor angle and the multiple objects of output electric current etc.
The influence of reason amount, therefore the amplitude of its voltage real-time change in transient process can improve dynamic response and improve system
Stability.However, current VSG implements this control strategy not yet, cause its response in transient process unsatisfactory.
Renewable energy mainly passes through that VSG is grid-connected, and current VSG is mainly using sagging control or inertia control, active frequency
Rate control and the decoupling of reactive power/voltage control two parts, this control mode can guarantee VSG safe and stable operation in the steady state, but
When failure occurs, system oscillation is violent, and dynamic response deteriorates, and system stabilization is on the hazard.
Fig. 3 is the control flow schematic diagram of the end voltage of virtual synchronous of embodiment of the present invention generator, as shown in figure 3, empty
The rate-determining steps of the end voltage of quasi- synchronous generator are as follows:
Step 301: the calculation formula about quantity of state being determined according to the control system of virtual synchronous generator, obtains institute
The nonlinear system model of control system is stated, the quantity of state includes the angular frequency and generator rotor angle of the virtual synchronous generator;
Step 302: the nonlinear system model is converted into Linear system model by changes in coordinates;
Step 303: the virtual synchronous generator voltage optimum control rate being carried out to the Linear system model and is asked
Solution, obtains the optimum control rate of linear system middle-end voltage;
Step 304: the optimum control rate generation of the linear system middle-end voltage being returned in the nonlinear system, is obtained non-
The optimum control rate of virtual synchronous generator voltage described in linear system;
Step 305: being regulated and controled using the optimum control rate of virtual synchronous generator voltage described in the nonlinear system
The amplitude of the virtual synchronous motor terminal voltage, to slow down the oscillation of system in transient process.
In step 301, it specifically includes: according to each controlling unit, the void in the control circuit of the control system
Electromagnetic torque in the input torque of quasi- synchronous generator, the control system determines the calculation formula of control process.
Obtain the electromagnetic torque, the electromagnetic torqueWherein, E is the virtual synchronous engine
End voltage, U is network voltage, XLFor line impedance, δ is the generator rotor angle of the virtual synchronous generator, and ω is virtual synchronous hair
Motor angular frequency;
It willIt is brought into the calculation formula of the control process, obtains the calculating about the quantity of state
Formula.
Fig. 4 is virtual synchronous generated power-frequency control system schematic diagram, as shown in figure 4, TmAnd TeRespectively indicate VSG
Input torque and electromagnetic torque.ω is virtual synchronous generator angular frequency, ω0For the specified angular frequency of power grid, D and J are respectively VSG
Damped coefficient and virtual inertia.The method provided in the control system and step 301 provided according to Fig. 4, available needle
Wherein to the calculation formula about the quantity of state of the control system, the virtual synchronous power generation
The end voltage E of machine is control variable, and it is the derivative of the quantity of state, X=[δ ω] that X, which is the quantity of state,T,
It in step 302, will be in nonlinear systemLinear system is transformed by coordinate transform
In system, obtainWherein, Z is the state variable under linear system,For leading for the state variable under linear system
Number, A are state variable coefficient matrix coefficient matrix, and B is the coefficient matrix of corresponding control amount, and v is linear system end voltage
Optimum control rate,
In step 303, the performance indicator formula that linear system seeks optimum control rate is obtainedWherein, Q and R is symmetrical positive definite weight matrix;
According to the formulaLinear optimal control is sought using linear optimal active control algorithm
Rate v processed.
Before being coordinately transformed, whether can be tested first to above system with exact linearization method.One affine
Nonlinear system is capable of the necessary and sufficient condition of exact linearization method are as follows:
1) near equalization point, matrixOrder do not change, and be equal to system rank
Number;Equalization point, that is, system incipient stability operation point, referred to ω=ω0, δ=δ0It is corresponding,
WhereinRefer to g along the n rank Lie bracket in the direction f.
2) vector fieldIt is pairing at equalization point.
Since system (3) is second-order system, its single order Lie derivatives need to be only calculated, i.e.,
By calculating it is found that system can exact linearization method in region Ω={ δ, ω | δ ≠ 0, π }.For actual motion
System, generator rotor angle meet 0 < δ < π, therefore are believed that the system can be by exact linearization method.
Inspection obtain above system can be with exact linearization method after, carry out converting linear system for nonlinear system
Specific steps:
1) D is enabled1={ g }, D2={ g, adfG }, select the vector field of two Line independentsMake its satisfactionA satisfactory group selection is:
2) mapping X=F (W) is calculated, corresponding equation isX0=[δ0,ω0]T
Wherein, w1And w2For the state variable under new mapping,It indicates with X0For initial value, using w2 as independent variable,
WithFor the integral of derivative, i.e.,
Indicate withFor initial value, with w1For independent variable, withFor the integral of derivative, know accordingly
Inverse mapping is
3) by calculating in F-1Under f (X) export mappingTo obtain f(0)(w), f(0)(w) for followed by
R1Vector field needed for transformation.
It inverts first and maps F-1Jacobian matrix, i.e.,
It can thus be concluded that
Original system is turned into linear system
Calculate transformation R1
Wherein,WithIt is R1State variable under transformation.
Definition converts T
Thus it acquires
Wherein,WithFor the intermediate variable of transform.
Definition
z1=w2=δ-δ0
Transformed linear system, which can be obtained, is
Wherein, Z=[z1,z2]TFor the state variable under linear system, A is state variable coefficient matrix coefficient matrix, and B is
The coefficient matrix of corresponding control amount, value are respectively
The solution that optimum control rate is carried out to the system after linearisation, initially sets up its performance indicator
Wherein, Q and R is symmetrical positive definite weight matrix, usually be can be taken as
Linear optimal contro8 rate v* useable linear optimal control algorithm acquires, the optimum control for being optimal performance indicator
Rate is
V=-R-1BTP*Z
Wherein, P* is symmetrical constant matrices, is Riccati non trivial solution, i.e. P* meets
ATP+PA-PBR-1BTP+Q=0
It can thus be appreciated that
v*=-z1-1.73z2
In this result generation, is back to original system, can obtain the corresponding optimum control rate of former nonlinear system is
The control method of the end voltage of virtual synchronous generator provided by the invention introduces frequency, generator rotor angle, line impedance
The end voltage of virtual synchronous generator is regulated and controled with parameters such as rotary inertias, make hold voltage control not only with voltage deviation
It is related with idle deviation, also accurately embody the influence of the physical quantity in other transient processes.By to nonlinear system line
Property, linear optimal contro8 rate is sought, its generation is returned in original system, obtains the optimum control rate under nonlinear system.To
The oscillation for slowing down system in transient process, improves dynamic performance, improves the stability of system.
In order to achieve the above objectives, the present invention also provides a kind of non-linear voltage control system, Fig. 5 is the embodiment of the present invention
The structural schematic diagram of the control system of the end voltage of virtual synchronous generator, as shown in figure 5, the system comprises:
Quantity of state calculation formula determination unit, for being determined according to the control system of virtual synchronous generator about state
The calculation formula of amount, obtains the nonlinear system model of the control system, and the quantity of state includes the virtual synchronous power generation
The angular frequency and generator rotor angle of machine;
Linear system model conversion unit, for the nonlinear system model to be converted to linear system by changes in coordinates
System model;
Linear system optimum control rate computing unit, for carrying out the virtual synchronous power generation to the Linear system model
The solution of set end voltage optimum control rate obtains the optimum control rate of linear system middle-end voltage;
Nonlinear system optimal control rate computing unit, for by the optimum control rate generation of the linear system middle-end voltage
It returns in the nonlinear system, obtains the optimum control rate of virtual synchronous generator voltage described in nonlinear system;
Voltage control unit is held, for utilizing the optimal of virtual synchronous generator voltage described in the nonlinear system
Control rate regulates and controls the amplitude of the virtual synchronous motor terminal voltage, to slow down the oscillation of system in transient process.
Wherein, the quantity of state calculation formula determination unit, specifically includes:
First state amount calculation formula determines subelement, for each control in the control circuit according to the control system
Link, the input torque of the virtual synchronous generator, the electromagnetic torque in the control system determine the calculating of control process
Formula;
Electromagnetic torque obtains subelement, for obtaining the electromagnetic torque, the electromagnetic torqueIts
In, E is the end voltage of the virtual synchronous engine, and U is network voltage, XLFor line impedance, δ is virtual synchronous power generation
The generator rotor angle of machine, ω are virtual synchronous generator angular frequency;
Second quantity of state calculation formula determines subelement, by that will be brought into based on the control process
Calculate formula in, obtain about the quantity of state formula wherein, the end of the virtual synchronous generator
Voltage E is control variable, and it is the derivative of the quantity of state, X=[δ ω] that X, which is the quantity of state,T,
The Linear system model conversion unit, specifically includes:
Linear system model transforming subunit, for by being become in nonlinear system by coordinate
It changes and is transformed into linear system, obtain wherein, Z is the state variable under linear system, is under linear system
State variable derivative, A be state variable coefficient matrix coefficient matrix, B be corresponding control amount coefficient matrix, v is line
The optimum control rate of property system end voltage,
The linear system optimum control rate computing unit, specifically includes:
Performance indicator formula obtains subelement, and the performance indicator formula of optimum control rate is sought for obtaining linear systemWherein, Q and R is symmetrical positive definite weight matrix;
Linear system optimum control rate computation subunit, for according to the formulaBenefit
Linear optimal contro8 rate v is sought with linear optimal active control algorithm.
Considered in end voltage (end voltage) control system of virtual synchronous generator provided by the invention frequency, generator rotor angle,
Influence of the parameters such as line impedance and rotary inertia to the end voltage of virtual synchronous generator, slows down system in transient process
Oscillation, improves dynamic performance, improves the stability of system.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other
The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment
For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part
It is bright.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said
It is bright to be merely used to help understand method and its core concept of the invention;At the same time, for those skilled in the art, foundation
Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.