CN106532685A - Generalized impedance criterion calculation method for stability analysis of grid-connected inverter and application - Google Patents
Generalized impedance criterion calculation method for stability analysis of grid-connected inverter and application Download PDFInfo
- Publication number
- CN106532685A CN106532685A CN201610946122.9A CN201610946122A CN106532685A CN 106532685 A CN106532685 A CN 106532685A CN 201610946122 A CN201610946122 A CN 201610946122A CN 106532685 A CN106532685 A CN 106532685A
- Authority
- CN
- China
- Prior art keywords
- network port
- inverter
- side network
- generalized
- impedance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
The invention discloses a generalized impedance criterion calculation method for stability analysis of a grid-connected inverter and an application. The method comprises the following steps of building a small signal impedance model of an inverter side network port and a small signal impedance model of a power grid side network port under polar coordinates; and calculating generalized impedances of the inverter side network port and the power grid side network port according to the built small signal impedance model of the inverter side network port and small signal impedance model of the power grid side network port under the polar coordinates, and dividing the two generalized impedances obtained in the second step to obtain the ratio of the generalized impedances of the inverter side network port and the power grid side network port as a generalized impedance criterion. The obtained generalized impedance criterion can be used for small signal stability analysis of a grid-connected inverter system, the problem that the analysis accuracy is affected by a traditional impedance method due to the presence of coupling is solved, a foundation is laid for clearing power grid oscillation generated when new energy is connected to a power grid on a large scale, and the method disclosed by the invention is an important method for designing an inverter controller to improve the stability.
Description
Technical field
The present invention relates to a kind of stabilization of power grids criterion computational methods and application, more particularly, to a kind of for grid-connected inverse
Become generalized impedance criterion computational methods and the application of device stability analysis.
Background technology
With the development of power system, increasing power electronic equipment starts to be applied to power system, including for
The application related to flexible AC transmitting system of the combining inverter of new forms of energy access, HVDC transmission system.Although electric power
Electronic equipment has the advantages that flexibly controllable, but the switching device of power electronic equipment can produce harmonic wave, harm electricity in electrical network
Source quality and system stability.When new forms of energy access electrical network, in order to improve power delivery capabilities and save line cost, usually passing
Defeated circuit adopts series capacitor compensation, but the interaction between series compensation capacitance and the control system of power electronic equipment
The sub-synchronous oscillation of electrical network may be caused.Between the multiple dynamic devices of power system, the multiple stage fan converter of such as wind field
Controller between interaction be also possible to cause system oscillation unstability.
Nowadays the method for analysis system stability mainly has：Time-domain-simulation method, characteristic value calculating method, complex torque coefficients,
Impedance method etc..Wherein, time-domain simulation method can directly observe the response curve of each state variable, it can be difficult to differentiating each
The damping characteristic of oscillation mode, and required amount of calculation is more huge, is not suitable for studying large scale system；Characteristic value is calculated
The details with regard to system oscillation can be drawn, but when system complex, the solution of characteristic value becomes difficult；Multiple moment coefficient
Method is a simple and effective analysis method, but the method is only applicable to one machine infinity bus system, it is difficult to tackle multi-computer system
Complex oscillation；Impedance method is another kind of easy analysis method, and the method is the combining inverter table accessed for new forms of energy
The output impedance being shown as on being connected to voltage source or current source.And power supply output impedance and network are analyzed using Nyquist criterion
The ratio of input impedance, so as to judge the stability of system.
The problem of Conventional impedance method is, when impedance model is set up under rest frame, to deposit between positive-negative sequence impedance
In coupling, ignoring coupling terms can make the judgement of the stability of a system obtain the conclusion of mistake.When output impedance is under synchronous coordinate system
During modeling, equally exist between the d axles and q axles of synchronous coordinate system and couple.Therefore common Nyquist criterion cannot still make
With.Although can so analyze complex with analysis system stability using broad sense Nyquist criterion, lose using resistance
The advantage of anti-method Simplified analysis.
The content of the invention
The invention discloses a kind of generalized impedance criterion computational methods and application for combining inverter stability analysis, profit
Grid-connected inverter system stability is judged and to analyze with Nyquist criterion, referred to as generalized impedance method.
The technical solution used in the present invention is that have three steps：
The first step, sets up the small signal impedance model and the grid side network port of the inverter side network port under polar coordinates
Small signal impedance model；
Second step, according to the small signal impedance model and the network port of the inverter side network port set up under polar coordinates
Small signal impedance model, calculating tries to achieve the generalized impedance of the inverter side network port and the grid side network port；
3rd step, two generalized impedances that second step is obtained is divided by and obtains the inverter side network port and the network port
Generalized impedance ratio, as generalized impedance criterion.
Described polar coordinates are the polar coordinates with amplitude and phase angle as representation.
For inverter and the one machine infinity bus system of electrical network composition, the wave filter in inverter adopts LCL filter, inverse
The control mode for becoming device is added bicyclic vector controlled using lock, and inverter meets following condition：The inner ring control of inverter has
The electric voltage feed forward and decoupling item of full remuneration, and consider that electric voltage feed forward is fully compensated the dynamic with inverter inner ring and phaselocked loop,
The line resistance of the outer shroud dynamic, the equivalent time delay of PWM links and filter inductance resistance and grid side network of inverter is ignored
Disregard.
Described inverter side network includes the output inductor of inverter and inverter, the inverter side network-side
Mouthful small signal impedance model expression be：
Wherein, Δ I is disturbed for network port current amplitude, and I is network port electric current stable state amplitude,For network port electricity
Stream phase angle disturbance, Δ U are network port voltage disturbance, and U is network port voltage stable state amplitude, and Δ δ is network port voltage phase
Angle disturbs,For inverter side network port admittance matrix, YgRepresent the non-zero of inverter side network port admittance matrix
Element, YgExpression formula is：
Wherein, s is Laplace operator, GpllThe transmission function of (s) for phaselocked loop, Gpll(s)=(Kppll+Kipll/ s)/s,
KppllFor the proportionality coefficient of PI links in phaselocked loop, KipllFor the integral coefficient of PI links in phaselocked loop, G2S () is current inner loop
The transmission function of PI controllers, G2(s)=Kp2+Ki2/ s, Kp2For the proportionality coefficient of inner ring PI link, Ki2For inner ring PI link
Integral coefficient, LfFor the output inductor of inverter.
The grid side network port includes the filter capacitor C of inverter LCL filter1, comprising LCL filter grid side
The substitutional connection inductance L of inductance1With circuit series compensation capacitance.
The small signal impedance model of the described grid side network port is first represented with the loop equation of below equation：
Wherein, L represents the component type of inductance, and C represents the component type of electric capacity, and k represents the element number of grid side, j
The node serial number of grid side is represented, j=1 represents the node serial number of the grid side network port, and j=2 represents substitutional connection inductance L1
With circuit series compensation capacitance C2Between node, φjFor the power-factor angle of node j, YL(1,1)Represent inductance L1With 1 phase of node
The admittance matrix of pass, YL(1,2)Represent the inductance L1 admittance matrixs related to node 2, YC(1,1)Represent electric capacity C1It is related to node 1
Admittance matrix, YC(2,2)Represent electric capacity C2The admittance matrix related to node 2；
YL(1,1)And YL(1,2)Calculated using following formula：
YC(1,1)And YC(2,2)Calculated using following formula：
Wherein, YL(k,j)Represent inductance LkThe admittance matrix related to node j, YC(k,j)Represent electric capacity CkIt is related to node j
Admittance matrix, s is Laplace operator, and ω is system angular frequency, and L represents inductance class component, and C represents electric capacity class component, and k is represented
The element number of grid side network, j represent the node serial number of grid side network, and j=1 represents that the node of the grid side network port is compiled
Number, j=2 represents substitutional connection inductance L1With circuit series compensation capacitance C2Between node, φjFor the power factor of node j
Angle, L1The substitutional connection inductance of the grid side network port is represented, its substitutional connection is grid side filter inductance and line inductance
Inductance and；
Using Y=(YL(1,1)+YC(1,1))-YL(1,2)(YL(2,2)+YC(2,2))-1YL(1,1)Simplify loop equation to eliminate internal section
Point, it is below equation to eliminate the grid side network port equation after internal node, and as the small-signal of the grid side network port
Impedance model：
Wherein, Y represents the small-signal admittance matrix of the grid side network port.
The small signal impedance model of the inverter side network port and the grid side network port adopts below equation table
Show：
Wherein, a represents that Δ I is affected by Δ U, and b represents that Δ I is affected by U Δ δ, and c is representedAffected by Δ U, d tables
ShowAffected by U Δ δ；
Then the generalized impedance of two kinds of network ports is tried to achieve using below equation：
ZG=(ad-bc)/a
Wherein, ZGRepresent the generalized impedance of the network port.
Specifically, the generalized impedance computing formula of the inverter side network port is as follows：
Specifically, the generalized impedance computing formula of the grid side network port is as follows：
ZG_grid=(Y (1,1) Y (2,2)-Y (1,2) Y (2,1))/Y (1,1)
Wherein, Y (1, be 1) element of Y matrix the first row first rows, Y (1, be 2) element of Y matrix the first row secondary series,
Y (2, it is 1) element of Y matrixes the second row first row, Y (2, it is 2) element of Y matrixes the second row secondary series.
Described generalized impedance ZGInverse be generalized admittance YG, generalized impedance ZGWith generalized admittance YGRelation similar to
Conventional impedance and the relation of admittance, generalized admittance YGDefinition in the following ways：
When network port current amplitude disturbance Δ I is 0, the generalized admittance of the network port is current phase angle correlative
With the ratio of network port voltage phase angle correlative U Δ δ, that is, meet：
In formula, I is stable state port current, and U is stable state port voltage, and Δ I is network port current disturbing,For network
Port current phase angle is disturbed, and Δ δ is disturbed for network port voltage phase angle.
Described generalized impedance criterion can be used for the stability for judging grid-connected inverter system.
The inventive method is finally calculated generalized impedance Z of the inverter side network portG_invWith the grid side network port
Generalized impedance ZG_gridRatio ZG_inv/ZG_gridAs the foundation or controller design of grid-connected inverter system judgement of stability
With reference to.It is embodied as drawing both nyquist plots, by judging that nyquist plot with the intersection point of negative real axis is
It is no beyond (- 1, j0) point or within obtain the whether stable result of grid-connected inverter system.
The invention has the beneficial effects as follows：
The present invention calculates the generalized impedance criterion for obtaining can be used for judging the stability of grid-connected inverter system, for analyzing
Inverter model there is the characteristic of decoupling, it is to avoid cause the coupling terms of Conventional impedance method in ignore inverter model
Error, improves the accuracy of analysis.
The present invention lays a good foundation to clear the electrical network vibration produced when extensive new forms of energy access electrical network, is to improve inversion
The important evidence of device controller stability.
Description of the drawings
Fig. 1 is that the generalized impedance method of the present invention is calculated and analysis process figure.
Inverter one machine infinity bus system structured flowcharts of the Fig. 2 for embodiment.
Fig. 3 is system by phasor diagram during little interference.
Fig. 4 be the ratio of the generalized impedance of the inverter side network port and the grid side network port under different line lengths how
Qwest's curve.
Fig. 5 is the system voltage current waveform under dq coordinate systems before and after the change of line inductance size.
Specific embodiment
Below in conjunction with the accompanying drawings and specific embodiment is described in further detail to the present invention.
The embodiment and its flow process that stability of a system judgement is carried out using the inventive method is as shown in Figure 1.
The first step, sets up the small signal impedance model and the grid side network port of the inverter side network port under polar coordinates
Small signal impedance model；
The instantiation with reference to shown in Fig. 2, it is considered to inverter one machine infinity bus system, its median filter adopt typical LCL
Filtering, the control strategy of inverter are added bicyclic vector controlled using modal lock, and inverter meets hypothesis below condition：
Assume 1：The inner ring control of inverter is with electric voltage feed forward and decoupling item, and considers that electric voltage feed forward is fully compensated.
Assume 2：Consider the dynamic of inverter inner ring and phaselocked loop, the equivalent time delay of the dynamic and PWM links of outer shroud is neglected
Slightly disregard.
Assume 3：Ignore filter inductance resistance and line resistance.
The dynamic characteristic of the phaselocked loop of described inverter is represented by：
The state equation of inverter filtering inductance is：
After considering decoupling and the feedforward factor of inner ring control, the dynamical equation of interior ring controller is：
Wherein, s is Laplace operator；UsdAnd UsqThe respectively d of the node that inverter filtering inductance is connected with inverter
Shaft voltage and q shaft voltages；UsdrAnd UsqrThe d shaft voltage reference values of the node that respectively inverter filtering inductance is connected with inverter
With q shaft voltage reference values；UdAnd UqThe d shaft voltages and q shaft voltages of the node that respectively inverter filtering inductance is connected with network；
IdAnd IqRespectively the d axles component of inverter output current and q axle components；IdrAnd IqrRespectively d axles of inverter output current
Component reference value and q axle component reference values；ω is the system angular frequency obtained by phaselocked loop；Gpll(s)=(Kppll+Kipll/s)/s
For the transmission function of phaselocked loop；ω0For system fundametal compoment angular frequency；LfFor the inductance value of inverter side filter inductance；G2(s)
For the transmission function of current inner loop PI controller, G2(s)=Kp2+Ki2/s。
Fig. 3 show system by phasor diagram during little interference.As seen from the figure, the synchronous seats of xy and dq two are had in the figure
Mark system, wherein xy coordinate systems are the synchronous angular frequency with grid sidesThe coordinate system of rotation, dq coordinate systems are then with phaselocked loop frequency
Rotatable coordinate axis on the basis of rate ω；θpllFor the angle between dq coordinate systems and xy coordinate systems, its size is by phaselocked loop according to inverse
Become the q shaft voltage U of the node that device filter inductance is connected with networkqIt is determined that；θIAnd θURespectively electric current and voltage are in dq coordinate systems
Phase angle,The phase angle of electric current and voltage in xy coordinate systems is respectively with δ.
Formula (1)-(3) are linearized, the dynamical equation under inverter microvariations can be obtained：
(4)-(6) can be rewritten as dynamic under the inverter microvariations represented under polar coordinates by the phasor relation according to Fig. 3
State equation：
According to Fig. 2, the angular relationship in above formula meets：
ΔθU=Δ δ-Δ θpll (10)
Simultaneous (10) and (7), can be by Δ θpllIt is expressed as with Δ δ：
Time delay due to not considering PWM links, therefore have：
Again by assuming 2, the vibration of present invention concern is time/vibration of supersynchronous frequency range and inner ring bandwidth far above outer shroud,
Therefore it is believed that outer shroud exports constant, i.e. Δ IdrWith Δ IqrIt is 0.Therefore, bring formula (8) and (9) into formula (13) to obtain：
(10)-angular relationship shown in (12) is substituted into formula (14) again to obtain：
In formula,
Formula (15) is the small signal impedance model of the inverter side network port under polar coordinates.
In the grid-connected typical circuit of single inverter shown in Fig. 2, network port circuit includes the output filtering of inverter
Electric capacity C1, it is considered to the effective transmission circuit inductance L after the net side inductance of wave filter1, circuit series compensation capacitance C2, Yi Jiwu
Poor bulk power grid.
With i, j represents element two ends node serial number respectively, represents element number with k, can obtain, linearizing inductance unit
Part state equation is：
Note φiFor the power-factor angle of i-th node, then formula (16) is represented by under polar coordinates：
For ease of writing, above formula can be denoted as：
ΔIk=YL(k,i)ΔUi-YL(k,j)ΔUj (18)
Wherein,
Similarly, linearizing capacitance state equation is：
It is expressed as under polar coordinates：
For ease of writing, (21) can be denoted as：
ΔIk=YC(k,i)ΔUi-YC(k,j)ΔUj (22)
Wherein,
Can obtain further according to the topological row network loop equation shown in Fig. 2：
Especially, when inverter is filtered using pure inductance, Y in formulaC(1,1)For 0, when circuit is without series compensation capacitance,
YC(2,2)For 0.
Can be calculated, the small-signal admittance matrix of the grid side network port is：
Y=(YL(1,1)+YC(1,1))-YL(1,2)(YL(2,2)+YC(2,2))-1YL(1,1) (25)
Under polar coordinates, the small signal impedance model of the grid side network port is：
Second step, according to the small signal impedance model and the network port of the inverter side network port set up under polar coordinates
Small signal impedance model, calculating tries to achieve the generalized impedance of the inverter side network port and the grid side network port；
Define (generalized admittance YG)：When port current amplitude disturbance Δ I is 0, the generalized admittance of port is current phase angle
CorrelativeWith the ratio of port voltage phase angle correlative U Δ δ, that is, meet：
Similar to Conventional impedance and the relation of admittance, the inverse of generalized admittance is generalized impedance ZG。
In formula, I is stable state port current, and U is stable state port voltage, and Δ I is that port current is disturbed,For port current
Phase angle is disturbed, and Δ δ is disturbed for port voltage phase angle.
If the small signal impedance model of polar coordinates lower port is：
Then the generalized impedance of the port can be tried to achieve according to following formula：
ZG=(ad-bc)/a (29)
Can according to the small signal impedance model of the inverter side network port and the grid side network port under (29) and polar coordinates
In the hope of generalized impedance Z of the inverter side network portG_invAnd generalized impedance Z of the grid side network portG_grid。
3rd step, two generalized impedances that second step is obtained is divided by and obtains the inverter side network port and the network port
Generalized impedance ratio, as generalized impedance criterion.
With reference to concrete example, application of the ratio of generalized impedance in terms of system stability analysis is illustrated.In order that asking
Topic simplifies, and the inverter considered in example assumes that DC voltage is constant, and assumes that inverter injects watt current I to electrical networkdref
With reactive current Iqref；Inverter is filtered using pure inductance, and inverters parameter is as shown in table 1.In simulations, it is considered to circuit
Inductance is changed into 0.7pu from 0.35pu, to simulate the synchronised grids for accessing varying strength.
Table 1 emulates the parameter of inverter model used
Using the inverter side network port and the generalized impedance of the grid side network port, whether system can stably be made
Judge.For the broad sense of the machine net reciprocal effect characteristic equation of the system, the inverter side network port and the grid side network port
The ratio of impedance is ZG_grid(s)/ZG_invS (), draws respectively when line inductance is 0.35pu and 0.7pu the Nai Kui corresponding to it
This Tequ line, as a result as shown in Figure 4.
As seen from the figure, the nyquist plot of open-loop transfer function Mid Frequency be approximately two turn clockwise it is ellipse
Circle, and negative real axis is passed through clockwise.With the increase of line length, the intersection point positioned at nyquist plot and negative real axis can be to the left
It is mobile.When line inductance is 0.35pu, intersection point does not surround (- 1, j0) point positioned at (- 1, j0) point right side, nyquist plot,
Therefore system stability.When line inductance is 0.7pu, intersection point is surrounded clockwise positioned at (- 1, j0) point left side, nyquist plot
(- 1, j0) point 2 is enclosed.It means that former closed-loop system characteristic equation has two limits for being located at RHP, therefore work as circuit
Series connection additional inductor and after dying down, system is unstable.
What the Digital Simulation carried out in MATLAB softwares demonstrated generalized impedance is used for having for stability of a system criterion
Effect property.Simulation model is built in MATLAB/SIMULINK softwares, the sampling of control system and controls frequency and is in emulation
The frequency of 4kHz, SVPWM is 4kHz, and simulation step length is 5us.
In emulation, in t=5s, line inductance is changed into 0.7pu from 0.35pu, represents the increase of contact line length.Sit in dq
Under mark system, the system voltage current waveform before and after line inductance size changes is as shown in Figure 5.Found out by figure, in series connection additional electrical
Before sense, the vibration of system gradually decays, system stability.After extra inductance of connecting, systems stay vibration will not
Decay, system there occurs unstability.Contrast MATLAB simulation results and theoretical analysis result, the two is consistent.
By above-mentioned simulated example, it can be seen that generalized impedance method proposed by the present invention can accurate analysis system it is steady
It is qualitative.Analysis process avoids Conventional impedance method because ignoring coupling terms and the problem of impact analysis accuracy.In the present invention
Generalized impedance not only can judge the stability of system than criterion, be also used as in inverter control system design, improve inverse
Become the foundation of device stability.The method achieve system depression of order analysis, for the stability analysis of complication system have it is important
Meaning.
Above-mentioned specific embodiment is used for illustrating the present invention, rather than limits the invention, the present invention's
In spirit and scope of the claims, any modifications and changes made to the present invention both fall within the protection model of the present invention
Enclose.
Claims (9)
1. a kind of generalized impedance criterion computational methods for combining inverter stability analysis, it is characterised in that mainly including following
Three steps：
The first step, set up under polar coordinates the small signal impedance model and the grid side network port of the inverter side network port it is little
Signal impedance model；
Second step, according to the little of the small signal impedance model and the network port of the inverter side network port set up under polar coordinates
Signal impedance model, calculating try to achieve the generalized impedance of the inverter side network port and the grid side network port；
3rd step, two generalized impedances that second step is obtained is divided by and obtains the wide of the inverter side network port and the network port
The ratio of adopted impedance, as generalized impedance criterion.
2. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 1, its
It is characterised by：Described polar coordinates are the polar coordinates with amplitude and phase angle as representation.
3. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 1, its
It is characterised by：For inverter and the one machine infinity bus system of electrical network composition, the wave filter in inverter adopts LCL filter,
The control mode of inverter is added bicyclic vector controlled using lock, and inverter meets following condition：The inner ring control tool of inverter
Have an electric voltage feed forward and decoupling item of full remuneration, and consider electric voltage feed forward be fully compensated and inverter inner ring and phaselocked loop it is dynamic
State, the line resistance of the outer shroud dynamic, the equivalent time delay of PWM links and filter inductance resistance and grid side network of inverter are neglected
Slightly disregard.
4. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 1, its
It is characterised by：Described inverter side network includes the output inductor of inverter and inverter, the inverter side network
The small signal impedance model expression of port is：
Wherein, Δ I is disturbed for network port current amplitude, and I is network port electric current stable state amplitude,For network port electric current phase
Angle disturbs, and Δ U is network port voltage disturbance, and U is network port voltage stable state amplitude, and Δ δ disturbed for network port voltage phase angle
It is dynamic,For inverter side network port admittance matrix, YgThe nonzero element of inverter side network port admittance matrix is represented,
YgExpression formula is：
Wherein, s is Laplace operator, GpllThe transmission function of (s) for phaselocked loop, Gpll(s)=(Kppll+Kipll/ s)/s, Kppll
For the proportionality coefficient of PI links in phaselocked loop, KipllFor the integral coefficient of PI links in phaselocked loop, G2S () is controlled for current inner loop PI
The transmission function of device processed, G2(s)=Kp2+Ki2/ s, Kp2For the proportionality coefficient of inner ring PI link, Ki2For the integration of inner ring PI link
Coefficient, LfFor the output inductor of inverter.
5. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 1, its
It is characterised by：The grid side network port includes the filter capacitor C of inverter LCL filter1, comprising LCL filter electrical network
The substitutional connection inductance L of side inductance1With circuit series compensation capacitance C2；
The small signal impedance model of the described grid side network port is first represented with the loop equation of below equation：
Wherein, L represents the component type of inductance, and C represents the component type of electric capacity, and k represents the element number of grid side, and j is represented
The node serial number of grid side, j=1 represent the node serial number of the grid side network port, and j=2 represents substitutional connection inductance L1With line
Road series compensation capacitance C2Between node, φjFor the power-factor angle of node j, YL(1,1)Represent inductance L1With the node of j=1
Related admittance matrix, YL(1,2)Represent the inductance L1 admittance matrixs related to the node of j=2, YC(1,1)Represent electric capacity C1And j=
The related admittance matrix of 1 node, YC(2,2)Represent electric capacity C2The admittance matrix related to the node of j=2；
YL(1,1)And YL(1,2)Calculated using following formula：
YC(1,1)And YC(2,2)Calculated using following formula：
Wherein, YL(k,j)Represent inductance LkThe admittance matrix related to node j, YC(k,j)Represent electric capacity CkThe admittance related to node j
Matrix, s is Laplace operator, and ω is system angular frequency, and L represents inductance class component, and C represents electric capacity class component, and k represents electrical network
The element number of side network, j represent the node serial number of grid side network, and j=1 represents the node serial number of the grid side network port, j
=2 represent substitutional connection inductance L1With circuit series compensation capacitance C2Between node, φjFor the power-factor angle of node j, L1
The substitutional connection inductance of the grid side network port is represented, its substitutional connection is the inductance of grid side filter inductance and line inductance
With；
Using Y=(YL(1,1)+YC(1,1))-YL(1,2)(YL(2,2)+YC(2,2))-1YL(1,1)Simplify loop equation to eliminate internal node,
It is below equation to eliminate the grid side network port equation after internal node, and as the small signal impedance of the grid side network port
Model：
Wherein, Y represents the small-signal admittance matrix of the grid side network port.
6. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 1, its
It is characterised by：The small signal impedance model of the inverter side network port and the grid side network port adopts below equation table
Show：
Wherein, a represents affecting parameters of the Δ I by Δ U, and b represents affecting parameters of the Δ I by U Δ δ, and c is representedAffected by Δ U
Parameter, d are representedBy the affecting parameters of U Δ δ；
Then the generalized impedance of two kinds of network ports is tried to achieve using below equation：
ZG=(ad-bc)/a
Wherein, ZGRepresent the generalized impedance of the network port.
7. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 6, its
It is characterised by：
Specifically, the generalized impedance computing formula of the described inverter side network port is as follows：
Specifically, the generalized impedance computing formula of the described grid side network port is as follows：
ZG_grid=(Y (1,1) Y (2,2)-Y (1,2) Y (2,1))/Y (1,1)
Wherein, Y (1, be 1) element of Y matrix the first row first rows, Y (1, be 2) element of Y matrix the first row secondary series, Y (2,
1) it is the element of Y matrixes the second row first row, and Y (2, it is 2) element of Y matrixes the second row secondary series.
8. a kind of generalized impedance criterion computational methods for combining inverter stability analysis according to claim 6, its
It is characterised by：Described generalized impedance ZGInverse be generalized admittance YG, generalized admittance YGDefinition in the following ways：
When network port current amplitude disturbance Δ I is 0, the generalized admittance of the network port is current phase angle correlativeWith net
The ratio of network port voltage phase angle correlative U Δ δ, that is, meet：
In formula, I is stable state port current, and U is stable state port voltage, and Δ I is network port current disturbing,For network port electricity
Stream phase angle disturbance, Δ δ are disturbed for network port voltage phase angle.
9. the application of the generalized impedance criterion for being obtained according to the arbitrary methods described of claim 1-8, it is characterised in that：Described is wide
Adopted impedance criterion can be used for the stability for judging grid-connected inverter system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610946122.9A CN106532685B (en) | 2016-10-26 | 2016-10-26 | For the generalized impedance criterion calculation method of gird-connected inverter stability analysis and application |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610946122.9A CN106532685B (en) | 2016-10-26 | 2016-10-26 | For the generalized impedance criterion calculation method of gird-connected inverter stability analysis and application |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106532685A true CN106532685A (en) | 2017-03-22 |
| CN106532685B CN106532685B (en) | 2019-02-01 |
Family
ID=58292819
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610946122.9A Active CN106532685B (en) | 2016-10-26 | 2016-10-26 | For the generalized impedance criterion calculation method of gird-connected inverter stability analysis and application |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106532685B (en) |
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526328A (en) * | 2016-12-08 | 2017-03-22 | 浙江大学 | Generalized impedance measurement and calculation method applicable to power grid and networking equipment |
| CN106992548A (en) * | 2017-05-22 | 2017-07-28 | 哈尔滨工业大学 | A kind of control method for improving grid-connection converter stability |
| CN107317357A (en) * | 2017-07-17 | 2017-11-03 | 湖南大学 | The small-signal admittance modeling of LCL type combining inverter and method for analyzing stability |
| CN107437818A (en) * | 2017-07-07 | 2017-12-05 | 华中科技大学 | Suppress the control method of light current three-phase LCL type grid-connection converter phaselocked loop small-signal disturbance off the net |
| CN107565545A (en) * | 2017-07-28 | 2018-01-09 | 西南交通大学 | A kind of tractor net coupled system stability criteria computational methods based on impedance return ratio matrix |
| CN107608933A (en) * | 2017-08-11 | 2018-01-19 | 浙江大学 | A kind of inverter generalized impedance measuring method for applying disturbance based on secondary side |
| CN107942172A (en) * | 2017-12-01 | 2018-04-20 | 华北电力科学研究院有限责任公司 | The damping characteristic test system and method for wind turbine and flexible DC power transmission association system |
| CN109103903A (en) * | 2018-09-13 | 2018-12-28 | 华北电力大学 | A kind of judgment method causing sub-synchronous oscillation for straight drive blower |
| CN109449958A (en) * | 2018-11-29 | 2019-03-08 | 云南电网有限责任公司 | A kind of double-fed blower grid-connected system method for analyzing stability |
| CN109450284A (en) * | 2018-11-23 | 2019-03-08 | 广东雅达电子股份有限公司 | A kind of three-phase inverter method for analyzing stability considering impedance coupler |
| CN109638880A (en) * | 2018-12-21 | 2019-04-16 | 浙江大学 | Current transformer grid-connected system stability improvement method based on admittance shaping range |
| CN109687511A (en) * | 2018-12-21 | 2019-04-26 | 浙江大学 | Current transformer grid-connected system vibrational stabilization judgment method based on compound circuit mode power |
| CN110011364A (en) * | 2019-04-29 | 2019-07-12 | 华北电力大学（保定） | A kind of active power load fluctuation that reduces is to the control method of system stabilizing influence |
| WO2019136944A1 (en) * | 2018-01-12 | 2019-07-18 | 清华大学 | Grid-side sub-synchronous damping controller all-condition optimization method and device |
| CN110032074A (en) * | 2019-05-22 | 2019-07-19 | 中国科学院光电技术研究所 | A kind of double compensation device design method of two-way feedforward disturbance observer |
| CN110286275A (en) * | 2019-05-22 | 2019-09-27 | 清华大学 | The method and device of system stability when a kind of determining equipment is grid-connected |
| CN110676880A (en) * | 2019-10-18 | 2020-01-10 | 湖南中能合创电气科技有限公司 | Three-phase inverter stability analysis method based on SISO system theory |
| CN111146807A (en) * | 2020-01-03 | 2020-05-12 | 浙江大学 | Method for judging small interference stability of converter grid-connected system under polar coordinate |
| CN111211573A (en) * | 2020-01-09 | 2020-05-29 | 中国科学院电工研究所 | Operation stability analysis method for AC/DC power distribution and utilization system |
| CN112003271A (en) * | 2020-08-10 | 2020-11-27 | 浙江大学 | Converter access alternating current micro-grid stability analysis method based on distributed impedance criterion |
| CN112217191A (en) * | 2019-07-11 | 2021-01-12 | 华北电力大学（保定） | Direct-current power distribution network stability analysis method based on node impedance matrix |
| CN112670992A (en) * | 2021-01-22 | 2021-04-16 | 上海交通大学 | Method and system for analyzing stability and correcting instability of power distribution network containing energy router |
| CN113419588A (en) * | 2021-07-01 | 2021-09-21 | 湘潭大学 | Two-stage cascade converter system stability analysis method based on impedance characteristics |
| CN113746344A (en) * | 2021-09-09 | 2021-12-03 | 上海海事大学 | Impedance model modeling method of high-frequency isolation two-stage battery energy storage converter |
| CN115688344A (en) * | 2022-12-29 | 2023-02-03 | 国网江西省电力有限公司电力科学研究院 | Broadband oscillation analysis method and system of multi-converter grid-connected system |
| WO2023010659A1 (en) * | 2021-08-05 | 2023-02-09 | 合肥工业大学 | Stability criterion for grid-connected converter control loop under weak grid condition |
| CN116305805A (en) * | 2023-01-31 | 2023-06-23 | 四川大学 | Model construction method for eliminating stability reduced-order analysis error of large-scale converter |
| CN117220330A (en) * | 2023-05-31 | 2023-12-12 | 山东大学 | Inverter impedance stability domain quantification method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060116857A1 (en) * | 2004-11-30 | 2006-06-01 | Sevic John F | Method and apparatus for model extraction |
| CN101692109A (en) * | 2009-10-13 | 2010-04-07 | 清华大学 | Impedance measurement method based on system discrimination |
| CN102495290A (en) * | 2011-12-28 | 2012-06-13 | 上海海事大学 | Obtaining apparatus for port impedance characteristic of alternating current power electronic module and obtaining method thereof |
| CN104638678A (en) * | 2015-03-16 | 2015-05-20 | 国网山西省电力公司电力科学研究院 | Multi-fan grid-connected resonance analysis method |
-
2016
- 2016-10-26 CN CN201610946122.9A patent/CN106532685B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060116857A1 (en) * | 2004-11-30 | 2006-06-01 | Sevic John F | Method and apparatus for model extraction |
| CN101692109A (en) * | 2009-10-13 | 2010-04-07 | 清华大学 | Impedance measurement method based on system discrimination |
| CN102495290A (en) * | 2011-12-28 | 2012-06-13 | 上海海事大学 | Obtaining apparatus for port impedance characteristic of alternating current power electronic module and obtaining method thereof |
| CN104638678A (en) * | 2015-03-16 | 2015-05-20 | 国网山西省电力公司电力科学研究院 | Multi-fan grid-connected resonance analysis method |
Cited By (43)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106526328A (en) * | 2016-12-08 | 2017-03-22 | 浙江大学 | Generalized impedance measurement and calculation method applicable to power grid and networking equipment |
| CN106526328B (en) * | 2016-12-08 | 2019-02-05 | 浙江大学 | It is a kind of suitable for the generalized impedance of power grid and networked devices measurement and calculation method |
| CN106992548A (en) * | 2017-05-22 | 2017-07-28 | 哈尔滨工业大学 | A kind of control method for improving grid-connection converter stability |
| CN107437818A (en) * | 2017-07-07 | 2017-12-05 | 华中科技大学 | Suppress the control method of light current three-phase LCL type grid-connection converter phaselocked loop small-signal disturbance off the net |
| CN107437818B (en) * | 2017-07-07 | 2019-08-30 | 华中科技大学 | Inhibit the control method of the small signal disturbance of three-phase LCL type grid-connection converter phaselocked loop under weak grid |
| CN107317357A (en) * | 2017-07-17 | 2017-11-03 | 湖南大学 | The small-signal admittance modeling of LCL type combining inverter and method for analyzing stability |
| CN107317357B (en) * | 2017-07-17 | 2019-08-30 | 湖南大学 | The small signal admittance modeling of LCL type gird-connected inverter and method for analyzing stability |
| CN107565545A (en) * | 2017-07-28 | 2018-01-09 | 西南交通大学 | A kind of tractor net coupled system stability criteria computational methods based on impedance return ratio matrix |
| CN107608933A (en) * | 2017-08-11 | 2018-01-19 | 浙江大学 | A kind of inverter generalized impedance measuring method for applying disturbance based on secondary side |
| CN107608933B (en) * | 2017-08-11 | 2020-05-19 | 浙江大学 | Inverter generalized impedance measurement method based on secondary side disturbance |
| CN107942172A (en) * | 2017-12-01 | 2018-04-20 | 华北电力科学研究院有限责任公司 | The damping characteristic test system and method for wind turbine and flexible DC power transmission association system |
| WO2019136944A1 (en) * | 2018-01-12 | 2019-07-18 | 清华大学 | Grid-side sub-synchronous damping controller all-condition optimization method and device |
| CN109103903A (en) * | 2018-09-13 | 2018-12-28 | 华北电力大学 | A kind of judgment method causing sub-synchronous oscillation for straight drive blower |
| CN109450284A (en) * | 2018-11-23 | 2019-03-08 | 广东雅达电子股份有限公司 | A kind of three-phase inverter method for analyzing stability considering impedance coupler |
| CN109449958A (en) * | 2018-11-29 | 2019-03-08 | 云南电网有限责任公司 | A kind of double-fed blower grid-connected system method for analyzing stability |
| CN109638880B (en) * | 2018-12-21 | 2022-04-05 | 浙江大学 | Converter grid-connected system stability improving method based on admittance shaping range |
| CN109638880A (en) * | 2018-12-21 | 2019-04-16 | 浙江大学 | Current transformer grid-connected system stability improvement method based on admittance shaping range |
| CN109687511A (en) * | 2018-12-21 | 2019-04-26 | 浙江大学 | Current transformer grid-connected system vibrational stabilization judgment method based on compound circuit mode power |
| CN109687511B (en) * | 2018-12-21 | 2020-08-11 | 浙江大学 | Converter grid-connected system oscillation stability judgment method based on complex circuit modal power |
| CN110011364B (en) * | 2019-04-29 | 2023-10-13 | 华北电力大学（保定） | Control method for reducing influence of active power load fluctuation on system stability |
| CN110011364A (en) * | 2019-04-29 | 2019-07-12 | 华北电力大学（保定） | A kind of active power load fluctuation that reduces is to the control method of system stabilizing influence |
| CN110286275A (en) * | 2019-05-22 | 2019-09-27 | 清华大学 | The method and device of system stability when a kind of determining equipment is grid-connected |
| CN110032074B (en) * | 2019-05-22 | 2022-04-19 | 中国科学院光电技术研究所 | Double compensator design method of double-path feedforward disturbance observer |
| CN110032074A (en) * | 2019-05-22 | 2019-07-19 | 中国科学院光电技术研究所 | A kind of double compensation device design method of two-way feedforward disturbance observer |
| CN112217191A (en) * | 2019-07-11 | 2021-01-12 | 华北电力大学（保定） | Direct-current power distribution network stability analysis method based on node impedance matrix |
| CN112217191B (en) * | 2019-07-11 | 2023-07-28 | 华北电力大学（保定） | DC power distribution network stability analysis method based on node impedance matrix |
| CN110676880B (en) * | 2019-10-18 | 2023-03-31 | 湖南中能合创电气科技有限公司 | Three-phase inverter stability analysis method based on SISO system theory |
| CN110676880A (en) * | 2019-10-18 | 2020-01-10 | 湖南中能合创电气科技有限公司 | Three-phase inverter stability analysis method based on SISO system theory |
| CN111146807A (en) * | 2020-01-03 | 2020-05-12 | 浙江大学 | Method for judging small interference stability of converter grid-connected system under polar coordinate |
| CN111211573A (en) * | 2020-01-09 | 2020-05-29 | 中国科学院电工研究所 | Operation stability analysis method for AC/DC power distribution and utilization system |
| CN112003271A (en) * | 2020-08-10 | 2020-11-27 | 浙江大学 | Converter access alternating current micro-grid stability analysis method based on distributed impedance criterion |
| CN112003271B (en) * | 2020-08-10 | 2021-12-07 | 浙江大学 | Converter access alternating current micro-grid stability analysis method based on distributed impedance criterion |
| CN112670992B (en) * | 2021-01-22 | 2023-11-07 | 上海交通大学 | Power distribution network stability analysis and instability correction method and system containing energy router |
| CN112670992A (en) * | 2021-01-22 | 2021-04-16 | 上海交通大学 | Method and system for analyzing stability and correcting instability of power distribution network containing energy router |
| CN113419588A (en) * | 2021-07-01 | 2021-09-21 | 湘潭大学 | Two-stage cascade converter system stability analysis method based on impedance characteristics |
| WO2023010659A1 (en) * | 2021-08-05 | 2023-02-09 | 合肥工业大学 | Stability criterion for grid-connected converter control loop under weak grid condition |
| CN113746344B (en) * | 2021-09-09 | 2023-08-01 | 上海海事大学 | Impedance model modeling method of high-frequency isolation two-stage battery energy storage converter |
| CN113746344A (en) * | 2021-09-09 | 2021-12-03 | 上海海事大学 | Impedance model modeling method of high-frequency isolation two-stage battery energy storage converter |
| CN115688344A (en) * | 2022-12-29 | 2023-02-03 | 国网江西省电力有限公司电力科学研究院 | Broadband oscillation analysis method and system of multi-converter grid-connected system |
| CN116305805A (en) * | 2023-01-31 | 2023-06-23 | 四川大学 | Model construction method for eliminating stability reduced-order analysis error of large-scale converter |
| CN116305805B (en) * | 2023-01-31 | 2023-12-19 | 四川大学 | Model construction method for eliminating stability reduced-order analysis error of large-scale converter |
| CN117220330A (en) * | 2023-05-31 | 2023-12-12 | 山东大学 | Inverter impedance stability domain quantification method and system |
| CN117220330B (en) * | 2023-05-31 | 2024-03-26 | 山东大学 | Inverter impedance stability domain quantification method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106532685B (en) | 2019-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106532685A (en) | Generalized impedance criterion calculation method for stability analysis of grid-connected inverter and application | |
| CN106786776B (en) | It is a kind of to utilize the method for correcting generalized impedance method analysis grid-connected inverter system stability | |
| CN109449958B (en) | A kind of double-fed blower grid-connected system method for analyzing stability | |
| CN104638678B (en) | A kind of grid-connected harmonic analysis method of multi fan | |
| CN105006825B (en) | A kind of electric power electric transformer and its control method of high quality of power supply output | |
| CN107093901B (en) | A kind of machine-electricity transient model and emulation mode of Distributed Power Flow controller | |
| CN108418253A (en) | The impedance modeling of current-control type virtual synchronous generator and method for analyzing stability | |
| CN110165709A (en) | Consider the virtual synchronous machine grid-connected inverting system stability method for improving of sampling time delay | |
| CN105119322B (en) | A kind of control method for improving grid-connected voltage source type inverter system stability | |
| CN108880300B (en) | Double-fed fan rectifier impedance calculation method based on double closed-loop control | |
| CN108923463A (en) | Consider the frequency coupling modeling method of the single-phase LCL type gird-connected inverter of phaselocked loop | |
| CN104113059A (en) | Node current injection method-based modeling method of unified power flow controller | |
| CN106410810A (en) | UPFC control unit with additional damping control function | |
| CN105914783A (en) | Island type microgrid small interference stability analysis method | |
| CN106549399A (en) | A kind of APF DC side voltage control methods in parallel based on sliding formwork PI complex control algorithms | |
| CN112467786B (en) | Small signal model analysis method for virtual synchronous machine of hybrid microgrid converter | |
| CN108847670A (en) | A kind of harmonic instability analysis method of double-fed blower grid side converter | |
| CN111654052A (en) | Flexible direct current converter modeling device and method based on dynamic phasor method | |
| CN108347060A (en) | A kind of power electronics interface power Reduced Modeling Methods and system | |
| CN111146807A (en) | Method for judging small interference stability of converter grid-connected system under polar coordinate | |
| CN110048447B (en) | H between flexible direct current transmission system station ∞ Decoupling controller | |
| CN113708399B (en) | DC voltage control time scale dynamic stability analysis method and device | |
| CN107425544A (en) | The electro-magnetic transient stabilising arrangement and method of a kind of distributed generation system | |
| Le et al. | Stability analysis of grid-connected inverter system containing virtual synchronous generator under time delay and parameter uncertainty | |
| CN112000018A (en) | Robust fault-tolerant control module, method and system based on residual generator |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |