CN106021682B - Sub-synchronous oscillation analysis method and device based on impedance network model - Google Patents
Sub-synchronous oscillation analysis method and device based on impedance network model Download PDFInfo
- Publication number
- CN106021682B CN106021682B CN201610320261.0A CN201610320261A CN106021682B CN 106021682 B CN106021682 B CN 106021682B CN 201610320261 A CN201610320261 A CN 201610320261A CN 106021682 B CN106021682 B CN 106021682B
- Authority
- CN
- China
- Prior art keywords
- impedance
- sub
- synchronous oscillation
- polymerization
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000010355 oscillation Effects 0.000 title claims abstract description 106
- 238000004458 analytical method Methods 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 claims abstract description 76
- 238000006116 polymerization reaction Methods 0.000 claims abstract description 75
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000013016 damping Methods 0.000 claims description 39
- 238000013459 approach Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000002776 aggregation Effects 0.000 claims description 12
- 238000004220 aggregation Methods 0.000 claims description 12
- 230000004044 response Effects 0.000 claims description 12
- 239000013598 vector Substances 0.000 claims description 11
- 230000005540 biological transmission Effects 0.000 claims description 9
- 230000006870 function Effects 0.000 description 14
- 230000008569 process Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The invention discloses a kind of sub-synchronous oscillation analysis methods and device based on impedance network model, wherein method includes: acquisition impedance model；Impedance network model is established by element interconnected method according to impedance model；Star triangular transformation is combined to be polymerize to obtain the polymerization impedance of impedance network model to impedance network model using series-parallel；According to the sub-synchronous oscillation characteristic of polymerization impedance analysis system.This method provides the new methods of sub-synchronous oscillation modeling and stability analysis, can be simple and efficient the stability and oscillating characteristic of decision-making system.
Description
Technical field
The present invention relates to Power System Stability Analysis technical fields more particularly to a kind of time same based on impedance network model
Walk analysis of the oscillation method and apparatus.
Background technique
Under the trend that modern power systems power electronics degree is increasingly deepened, electric electronic current change technology not only changes
The form of power grid primary equipment, and its distinctive multiple dimensioned sequential controlling mechanism and low inertia broadband response characteristic, it is just aobvious
The behavioral characteristics for changing electric system are write, new stability problem is brought.In recent years, especially with Flow varying power supplies such as wind-powered electricity generation, photovoltaics
Novel sub-synchronous oscillation (Subsynchronous Oscillation, SSO) problem caused by extensive access power grid is especially prominent
Out.In addition, the negative resistance property that current transformer formula constant power load embodies in subsynchronous frequency, the phaselocked loop of adjacent current transformer
(PLL) return circuit coupling and current transformer control participate in the series/parallel resonant tank that filter circuit and power grid are constituted, and may all swash
Send out sub-synchronous oscillation risk.This magnanimity converter device (wind-powered electricity generation, photovoltaic, inverter etc.)-power grid interaction causes new
Type sub-synchronous oscillation seriously threatens the equipment safety of modern power network, system stabilization and power quality, has become restriction wind in China
The bottleneck factor that the new energy such as light dissolve on a large scale.
For the analysis method of electric system power electronics bring stability problem, to continue to use traditional electromagnetism mostly temporary
Impedance analysis method under state emulation, Eigenvalues analysis or simplified condition.The electromagnetic transient simulation software being widely used at present is (such as
PSCAD, EMTP) only support the modeling and analysis of single step-length (fixed rate), it is difficult to it is changeable to realize to take into account precision and efficiency
Flow the time-domain-simulation under device, complex AC/DC network, magnanimity operating condition；Traditional Eigenvalues analysis method faces dimension calamity problem,
It is difficult to be suitable for the practical bulk power grid including thousands of current transformers, tens of thousands of a nodes.In recent years, it is widely used in power electronics to set
Standby and electric system repercussion study impedance model analytic approach provides new approaches.In practical applications, impedance model has
There is following advantage: 1) being derived the impedance model and total system impedance model of available each subsystem, and physics is anticipated
Justice is relatively unambiguous；2) when changing system parameter, the impedance model of one or several subsystem is only influenced, to overall impedance model shadow
Sound is small；3) the Nyquist stability criterion based on impedance model can be used and judge system stability, it is visual in image.Under simplified condition
Impedance analysis method in order to derive conveniently, corresponding simplification is had also been made to the control strategy of controller in system, has ignored portion
The dynamic characteristic of sub-controller.Although these simplify the foundation that operation is conducive to system impedance model, bringing be can not be ignored
Analytical error.
Summary of the invention
The purpose of the present invention is intended to solve above-mentioned one of technical problem at least to a certain extent.
For this purpose, the first purpose of this invention is to propose a kind of sub-synchronous oscillation analysis side based on impedance network model
Method, this method provides the new methods of sub-synchronous oscillation modeling and stability analysis, can be simple and efficient the stability of decision-making system
And oscillating characteristic.
Second object of the present invention is to propose a kind of sub-synchronous oscillation analytical equipment based on impedance network model.
In order to achieve the above object, according to a first aspect of the present invention embodiment propose it is a kind of time same based on impedance network model
Walk analysis of the oscillation method, comprising the following steps: obtain impedance model；It is established according to the impedance model by element interconnected method
Impedance network model；Using in series and parallel in conjunction with the Equivalent Networks such as star triangular transformation method to impedance network model progress abbreviation
Obtain the polymerization impedance of the impedance network model；According to the sub-synchronous oscillation characteristic of the polymerization impedance analysis system.
The sub-synchronous oscillation analysis method based on impedance network model of the embodiment of the present invention, first acquisition impedance model,
Impedance network model is established by element interconnected method then according to impedance model, then using in series and parallel in conjunction with star triangular transformation etc.
Equivalent Network method polymerize impedance network model to obtain the polymerization impedance of impedance network model, finally according to polymerization impedance
The sub-synchronous oscillation characteristic of analysis system.This method provides the new method of sub-synchronous oscillation modeling and stability analysis, Ke Yijian
The stability and oscillating characteristic of single efficiently decision-making system.
In addition, the sub-synchronous oscillation analysis method according to the above embodiment of the present invention based on impedance network model can be with
With following additional technical characteristic:
In some instances, the acquisition impedance model specifically includes: using Small Signal Model Analysis Method, initially sets up each
The condition of small signal equation model of component, using component i exit voltage and current as interface variables, the shape of each component linearisation
State equation model can arrange as canonical form:Wherein, Δ indicates incremental computations；ΔxiExpression removes
Other state variable increment column vectors outside port voltage, electric current；Δui、ΔiiRespectively indicate port voltage, current increment arrange to
Amount；aijIndicate the coefficient matrix of respective dimensions；I, j ∈ I={ 1,2,3 }, subscript i indicate i-th of component；Become by Laplce
It changes, algebraic equation model of the state equation model in the domain s is expressed as:
Wherein, s indicates Laplace operator；The exit component i voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween pass
System, such as following formula: Δ ui(s)=Zi(s)·Δii(s), wherein Zi(s) the equivalent impedance model of component i is indicated.
In some instances, when component i is one port network component, the ZiIt (s) is single transmission function；As component i
For multiport network component, the ZiIt (s) is transfer function matrix.
It in some instances, include: frequency characteristic according to the method for the sub-synchronous oscillation of the polymerization impedance analysis system
Analytic approach, polymerization RLC series circuit analytic approach, zero pole point stability methods.
It in some instances, is frequency-response analysis method according to the method for the polymerization impedance analysis system sub-synchronous oscillation
When, analyze the real part imaginary part of the polymerization impedance with frequency variation curve, when the polymerization impedance imaginary part there are zero crossing and
When zero crossing is located in subsynchronous frequency range, if corresponding real part is positive value, sub-synchronous oscillation has positive damping, secondary same
Step oscillation convergence, if corresponding real part is zero or negative value, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation is held
Continuous or diverging；It, will when method according to the polymerization impedance analysis system sub-synchronous oscillation is polymerization RLC series circuit analytic approach
The model aggregation for polymerizeing impedance is an equivalent Second-Order RLC Filter Circuit series circuit, and when R is greater than 0, sub-synchronous oscillation has positive damping,
Sub-synchronous oscillation is stablized, and when R is less than or equal to 0, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or hair
It dissipates；When method according to the sub-synchronous oscillation of the polymerization impedance analysis system is zero pole point stability methods, described in calculating
It polymerize the zero pole point of impedance, when the zero pole point real part of the polymerization impedance is all negative, system is stablized；When described poly-
When closing the zero pole point of impedance has positive real part, system is unstable.
In order to achieve the above object, a kind of time based on impedance network model that embodiment proposes according to the second aspect of the invention
Synchronized oscillation analytical equipment, comprising: module is obtained, for obtaining impedance model；Module is established, for according to the impedance model
Impedance network model is established by element interconnected method；Aggregation module, for using in series and parallel in conjunction with networks such as star triangular transformations
Equivalence method polymerize the impedance network model to obtain the polymerization impedance of the impedance network model；Analysis module is used
In the sub-synchronous oscillation characteristic according to the polymerization impedance analysis system.
The sub-synchronous oscillation analytical equipment based on impedance network model of the embodiment of the present invention, first acquisition module obtain resistance
Anti- model then sets up module according to impedance model and establishes impedance network model by element interconnected method, then aggregation module is adopted
The Equivalent Networks methods such as star triangular transformation are combined to be polymerize to obtain impedance network model to impedance network model with series-parallel
It polymerize impedance, ultimate analysis module is according to the sub-synchronous oscillation characteristic for polymerizeing impedance analysis system.This method provides subsynchronous
The new method of oscillation modeling and stability analysis, can be simple and efficient the stability and oscillating characteristic of decision-making system.
In addition, the sub-synchronous oscillation analytical equipment according to the above embodiment of the present invention based on impedance network model can be with
With following additional technical characteristic:
In some instances, the acquisition module is specifically used for: using Small Signal Model Analysis Method, initially sets up each component
Condition of small signal equation model, using component i exit voltage and current as interface variables, the state side of each component linearisation
Journey model can arrange as canonical form:Wherein, Δ indicates incremental computations；ΔxiIt indicates to remove port
Other state variable increment column vectors outside voltage, electric current；Δui、ΔiiRespectively indicate port voltage, current increment column vector；
aijIndicate the coefficient matrix of respective dimensions；I, j ∈ I={ 1,2,3 }, subscript i indicate i-th of component；By Laplace transform,
Algebraic equation model of the state equation model in the domain s is expressed as:Its
In, s indicates Laplace operator；The exit component i voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relationship,
Such as following formula: Δ ui(s)=Zi(s)·Δii(s), wherein Zi(s) the equivalent impedance model of component i is indicated.
In some instances, when component i is one port network component, the ZiIt (s) is single transmission function；As component i
For multiport network component, the ZiIt (s) is transfer function matrix.
In some instances, the analysis method of the analysis module includes: frequency-response analysis method, polymerization RLC series electrical
Road analytic approach, zero pole point stability methods.
In some instances, when the analysis method of the analysis module is frequency-response analysis method, the polymerization resistance is analyzed
Anti- real part imaginary part with frequency variation curve, when the polymerization impedance imaginary part is there are zero crossing and zero crossing is located at subsynchronous frequency
When within the scope of rate, if corresponding real part is positive value, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stablized, if right
Answering real part is zero or negative value, then sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging；The analysis
It is an equivalent Second-Order RLC Filter Circuit by the model aggregation for polymerizeing impedance when the analysis method of module is polymerization RLC series circuit analytic approach
Series circuit, when R is greater than 0, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation convergence is subsynchronous when R is less than or equal to 0
Oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging；The analysis method of the analysis module is that zero pole point is steady
When way of qualitative analysis, the zero pole point of the polymerization impedance is calculated, the zero pole point real part and if only if the polymerization impedance is all
When negative, system is stablized；When the zero pole point of the polymerization impedance has positive real part, system is unstable.
The additional aspect of the present invention and advantage will be set forth in part in the description, and will partially become from the following description
Obviously, or practice through the invention is recognized.
Detailed description of the invention
Above-mentioned and/or additional aspect of the invention and advantage will become from the description of the embodiment in conjunction with the following figures
Obviously and it is readily appreciated that, in which:
Fig. 1 is the process of the sub-synchronous oscillation analysis method according to an embodiment of the invention based on impedance network model
Figure；
Fig. 2 is the schematic diagram of power electronics electric system according to an embodiment of the invention；
Fig. 3 is impedance network model schematic according to an embodiment of the invention；
Fig. 4 is the equivalent aggregation process schematic of impedance network model according to an embodiment of the invention；
Fig. 5 is impedance (real part, void of power electronics electric system polymerization impedance according to an embodiment of the invention
Portion)-frequency characteristic；
Fig. 6 is the structure of the sub-synchronous oscillation analytical equipment according to an embodiment of the invention based on impedance network model
Schematic diagram.
Specific embodiment
The embodiment of the present invention is described below in detail, examples of the embodiments are shown in the accompanying drawings, wherein from beginning to end
Same or similar label indicates same or similar element or element with the same or similar functions.Below with reference to attached
The embodiment of figure description is exemplary, it is intended to is used to explain the present invention, and is not considered as limiting the invention.
The invention proposes the sub-synchronous oscillation analysis methods based on impedance network model of power electronics electric system
And device, i.e., converter device is equivalent to input-output impedance model, is later impedance network with system simultaneous, and will resistance
Anti- network is gradually polymerized to an impedance, by analyzing the frequency characteristic of impedance, realizes the judgement of system SSO stability characteristic (quality).
Below with reference to the accompanying drawings describe the embodiment of the present invention sub-synchronous oscillation analysis method based on impedance network model and
Device.
Fig. 1 is the process of the sub-synchronous oscillation analysis method according to an embodiment of the invention based on impedance network model
Figure.
As shown in Figure 1, the sub-synchronous oscillation analysis method based on impedance network model of the embodiment of the present invention includes following
Step:
Step 101, impedance model is obtained.
Specifically, power electronic equipment and the Controlling model of other power system component transient Models and full size, analysis
The Inner Constitution and connection relationship of each component in power electronics electric system, for various groups in power electronics electric system
Part, the input voltage of N-port component and the relationship of input current are as follows under the domain s:Switch to Vector-Matrix Form are as follows: V (s)=Z (s) I (s).
Wherein, Vi(s)、IiIt (s) is i-th of the port voltage and port current of N-terminal component；ZijIt (s) is impedance matrix pair
The impedance between port is answered, there is following form:
Wherein, in which: Z (s) be impedance matrix, Z (s) depend on component type, the connection relationship of internal component and its
Control system.For one port network component, Z (s) is single transmission function form；For multiport network component, Z (s) is
The matrix of transmission function.
It should be noted that in order to obtain the impedance model of each component of power electronics electric system, can use but not
It is limited to Small Signal Model Analysis Method, Model Distinguish method.
For example, using Small Signal Model Analysis Method, the condition of small signal equation model of each component is initially set up, with group
Part i exit voltage and current is interface variables, and the state equation model of each component linearisation can arrange as canonical form:Wherein, Δ indicates incremental computations；ΔxiIndicate that other states in addition to port voltage, electric current become
Measure increment column vector；Δui、ΔiiRespectively indicate port voltage, current increment column vector；aijIndicate the coefficient square of respective dimensions
Battle array；I, j ∈ I={ 1,2,3 }, subscript i indicate i-th of component.
By Laplace transform, algebraic equation model of the state equation model in the domain s is expressed as:Wherein, s indicates Laplace operator.
The exit component i voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relationship, such as following formula: Δ ui(s)
=Zi(s)·Δii(s), wherein Zi(s) the equivalent impedance model of component i is indicated.
Step 102, impedance network model is established by element interconnected method according to impedance model.
Specifically, it is described in detail in conjunction with Fig. 2, as shown in Fig. 2, with the large-scale AC and DC power of certain converter device containing magnanimity
For system, Load flow calculation is carried out under specific operation, is obtained various components port electrical parameter, is obtained according to aforementioned modeling method
After taking the impedance of various components, system is constructed using element interconnected method (Component Connection Method, CCM)
Whole impedance network model.I.e. on the basis of original electric system topological structure, various components are used into its impedance model
It replaces, is connected between the component impedance with electric connecting relation using ideal line, so that each impedance simultaneous be got up i.e.
Impedance network model as shown in Figure 3 can be obtained.
In modeling process, equivalent, such as wind power plant's blower type and operating condition is consistent can be carried out according to certain principle
Unit can be used single large capacity equivalence unit come it is equivalent, to simplify operation, but it is general not to its blower interior control system into
Row simplifies.
Step 103, the Equivalent Networks methods such as star triangular transformation are combined to polymerize impedance network model using series-parallel
Obtain the polymerization impedance of impedance network model.
Specifically, after obtaining power electronics power system impedance network model, abbreviation is carried out to impedance network.Abbreviation
When, abbreviation can be carried out using the method for series-parallel knot and the transformation of star triangular network according to the connection relationship between impedor；
Abbreviation can also be carried out by the transformation calculations between impedance matrix.Finally obtain the polymerization impedance of whole network, expression formula one
As be more complicated transmission function form.
For example: Fig. 3 show the large-scale AC and DC power system impedance network model of Fig. 2 device of device containing multiple-variable flow,
The polymerization impedance Z (s) of whole system is obtained after the abbreviation process shown in Fig. 4.First by 1 impedance of bulk power grid and HVDC impedance string
Connection, photovoltaic power plant impedance is connected with 1 impedance of route, wind field impedance is connected with 2 impedance of route, current transformer formula load impedance and route
6 impedances series connection, i.e.,
ImpedanceAnd impedanceIt connects after parallel connection with 3 impedance of route, i.e.,
ImpedanceImpedanceWith impedanceParallel connection, i.e.,
4 impedance of route, 5 impedance of route and tie line impedance carry out delta-Y transformation, i.e.,
ImpedanceWith impedance ZY1(s) series connection, impedance ZY2(s) with 2 impedance Z of bulk power gridS2(s) series connection, impedance ZS2(s)
With 3 impedance Z of bulk power gridS3(s) it connects, i.e.,
ImpedanceWith impedanceParallel connection, i.e.,
Impedance ZY(s) and impedanceSystematic polymerization impedance Z (s) is obtained after series connection, i.e.,
Step 104, according to the sub-synchronous oscillation characteristic of polymerization impedance analysis system.
Specifically, it after the equivalent impedance for obtaining power electronics power system network, can be carried out using a variety of methods
Stability and oscillating characteristic analysis, it is typical as frequency-response analysis method, polymerization RLC series circuit analytic approach, zero pole point are stablized
Property analytic approach, (broad sense) Nyquist stability methods etc..
For example, zero pole point method for analyzing stability.After the transmission function for obtaining systematic polymerization impedance Z (s), zero can be taken
Pole stability methods judge the stability of system, the i.e. zero pole point of computing system polymerization impedance Z (s), are polymerize by judgement
The real and imaginary parts feature of impedance Z (s) zero pole point judges the oscillating characteristic of system.
When the zero pole point real part for polymerizeing impedance Z (s) is all negative, system is stablized；When polymerization impedance Z (s)
When zero pole point has positive real part, system is unstable.Particularly, when the imaginary part pair of polymerization impedance Z (s) some or certain several zero pole points
When resonance frequency being answered to be located in subsynchronous frequency range, there is a possibility that SSO occurs in system, if corresponding to zero pole point real part at this time
It is negative, SSO convergence；If corresponding zero pole point real part is non-negative, exists and continue or dissipate SSO risk.
For example, frequency-response analysis method.Frequency characteristic is carried out to power electronics electric system frequency domain polymerization impedance Z (s)
Analysis can be analyzed using polymerization impedance Z (s) real part, imaginary part with frequency variation characteristics.
Analysis polymerization impedance Z (s) real part imaginary part with frequency variation curve, when polymerization impedance Z (s) imaginary part there are zero crossing and
When zero crossing is located in subsynchronous frequency range, if its corresponding real part is positive value at this time, show that SSO has positive damping, SSO is steady
It is fixed, if corresponding real part is zero or negative value, show that SSO has zero damping or negative damping, SSO persistent oscillation or diverging.
It is illustrated in figure 5 under certain power electronics electric system operating condition and polymerize impedance real part imaginary part with frequency variation song
Line, it can be seen that imaginary part curve zero crossing respective frequencies are located in subsynchronous frequency range, and real part is negative under respective frequencies
Value, so SSO has negative damping under operating condition, SSO is unstable.
For example, polymerization RLC series circuit analytic approach.It is bent to polymerize impedance Z (s) real part imaginary part frequency characteristic by analysis system
After line, system resonance frequencies ω can be obtainedr.Near series resonance frequency ω | 0≤| ω-ωr| < h } (h is the normal of very little
Number), equivalent impedance model Z (s) is polymerized to an equivalent Second-Order RLC Filter Circuit series circuit.Wherein, equivalent resistance R takes Z (j ωr)
Real part, and the parameter value of equivalent inductance L and equivalent capacity C can be by quasi- by nonlinear parameter in resonance frequency vertex neighborhood
It closes optimization method (as used least square method) to be calculated, polymerization impedance model is further equivalent to polymerization rlc circuit.
Based on obtained polymerization second-order circuit parameter, the damping and frequency of oscillation (i.e. SSO damping of second-order circuit can be calculated
And frequency), and then carry out quantization SSO analysis, calculation formula is as follows:
As a result, as R > 0, SSO has positive damping, and SSO stablizes；Conversely, SSO has zero damping or negative damping, SSO is lasting
Oscillation or diverging.
It is corresponding with the sub-synchronous oscillation analysis method provided by the above embodiment based on impedance network model, it is of the invention
A kind of embodiment also provides a kind of sub-synchronous oscillation analytical equipment based on impedance network model, since the embodiment of the present invention provides
Sub-synchronous oscillation analysis method based on impedance network model and above-mentioned several embodiments provide based on impedance network model
Sub-synchronous oscillation analysis method there is the same or similar technical characteristic, therefore aforementioned time same based on impedance network model
The embodiment of step analysis of the oscillation method is also applied for the sub-synchronous oscillation provided in this embodiment based on impedance network model point
Analysis apparatus is not described in detail in the present embodiment.As shown in fig. 6, should be analyzed based on the sub-synchronous oscillation of impedance network model
Device can include: obtain module 10, establish module 20, aggregation module 30 and analysis module 40.
Wherein, module 10 is obtained for obtaining impedance model.
Module 20 is established for establishing impedance network model by element interconnected method according to impedance model.
Aggregation module 30 be used for using combine in series and parallel the Equivalent Networks method such as star triangular transformation to impedance network model into
Row polymerization obtains the polymerization impedance of impedance network model.
Analysis module 40 is used for the sub-synchronous oscillation characteristic according to polymerization impedance analysis system.
In some instances, it obtains module 10 to be specifically used for: using Small Signal Model Analysis Method, initially setting up each component
Condition of small signal equation model, using component i exit voltage and current as interface variables, the state equation of each component linearisation
Model can arrange as canonical form:Wherein, Δ indicates incremental computations；ΔxiIt indicates except port electricity
Other state variable increment column vectors outside pressure, electric current；Δui、ΔiiRespectively indicate port voltage, current increment column vector；aij
Indicate the coefficient matrix of respective dimensions；I, j ∈ I={ 1,2,3 }, subscript i indicate i-th of component.
By Laplace transform, algebraic equation model of the state equation model in the domain s is expressed as:Wherein, s indicates Laplace operator.
The exit component i voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relationship, such as following formula: Δ ui(s)
=Zi(s)·Δii(s), wherein Zi(s) the equivalent impedance model of component i is indicated.
In some instances, when component i is one port network component, the ZiIt (s) is single transmission function.As component i
For multiport network component, the ZiIt (s) is transfer function matrix.
In some instances, the analysis method of analysis module 40 includes: frequency-response analysis method, polymerization RLC series circuit
Analytic approach, zero pole point stability methods.
In some instances, when the analysis method of analysis module 40 is frequency-response analysis method, the reality of analysis polymerization impedance
Portion's imaginary part with frequency variation curve, when polymerization impedance imaginary part is there are zero crossing and zero crossing is located in subsynchronous frequency range
When, if corresponding real part is positive value, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stablized, if corresponding real part is
Zero or negative value, then sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging.
When the analysis method of analysis module 40 is polymerization RLC series circuit analytic approach, it is by the model aggregation for polymerizeing impedance
One equivalent Second-Order RLC Filter Circuit series circuit, when R is greater than 0, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stablized, when R is small
When being equal to 0, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging.
The analysis method of analysis module be zero pole point stability methods when, calculate polymerization impedance zero pole point, when and only
When the zero pole point real part for polymerizeing impedance is all negative, system is stablized；When the zero pole point for polymerizeing impedance has positive real part, system
It is unstable.
The present invention is applicable not only to traditional subsynchronous oscillation of electrical power system analysis, and is suitable for inclusion in a large amount of electric power electricity
The electric system of son equipment provides a kind of sub-synchronous oscillation and builds especially for wind, the high electric system of light power generation permeability
The new method of mould and stability analysis.Impedance model is a kind of external characteristics model, for because of various reasons, (such as manufacturing firm's secrecy is needed
Want) it can not provide for internal control system model and the electrical circuitry equipment of parameter, one kind is provided accurately without regard in equipment
The expression method of portion's private information, and the order and parameter of model can be determined by discrimination method.Since impedance is convenient for
Using simple circuit transformations (such as series-parallel, star-triangular transformation) Lai Shixian abbreviation and polymerization, to overcome the small letter of tradition
The dimension disaster topic that the modeling of number state equation faces.After impedance network is aggregated into single impedance, simple impedance zero can be used
Pole, impedance frequency characteristics etc. quantify the stability and oscillating characteristic of decision-making system, clear with physical meaning, succinct high
The advantages of effect.
The sub-synchronous oscillation analytical equipment based on impedance network model of the embodiment of the present invention, first acquisition module obtain resistance
Anti- model then sets up module according to impedance model and establishes impedance network model by element interconnected method, then aggregation module is adopted
The Equivalent Networks methods such as star triangular transformation are combined to be polymerize to obtain impedance network model to impedance network model with series-parallel
It polymerize impedance, ultimate analysis module is according to the sub-synchronous oscillation for polymerizeing impedance analysis system.This method provides sub-synchronous oscillations
The new method of modeling and stability analysis, can be simple and efficient the stability and oscillating characteristic of decision-making system.
In the description of the present invention, it is to be understood that, term " first ", " second " are used for description purposes only, and cannot
It is interpreted as indication or suggestion relative importance or implicitly indicates the quantity of indicated technical characteristic.Define as a result, " the
One ", the feature of " second " can explicitly or implicitly include at least one of the features.In the description of the present invention, " multiple "
It is meant that at least two, such as two, three etc., unless otherwise specifically defined.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is included at least one embodiment or example of the invention.In the present specification, schematic expression of the above terms are not
It must be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be in office
It can be combined in any suitable manner in one or more embodiment or examples.In addition, without conflicting with each other, the skill of this field
Art personnel can tie the feature of different embodiments or examples described in this specification and different embodiments or examples
It closes and combines.
Any process described otherwise above or method description are construed as in flow chart or herein, and expression includes
It is one or more for realizing specific logical function or process the step of executable instruction code module, segment or portion
Point, and the range of the preferred embodiment of the present invention includes other realization, wherein can not press shown or discussed suitable
Sequence, including according to related function by it is basic simultaneously in the way of or in the opposite order, Lai Zhihang function, this should be of the invention
Embodiment person of ordinary skill in the field understood.
Those skilled in the art are understood that realize all or part of step that above-described embodiment method carries
It suddenly is that relevant hardware can be instructed to complete by program, the program can store in a kind of computer-readable storage medium
In matter, which when being executed, includes the steps that one or a combination set of embodiment of the method.
Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example
Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, modifies, replacement and variant.
Claims (4)
1. a kind of sub-synchronous oscillation analysis method based on impedance network model, which comprises the following steps:
Obtain impedance model；
Impedance network model is established by element interconnected method according to the impedance model；
Using being polymerize to obtain the impedance network model in series and parallel to the impedance network model in conjunction with star triangular transformation
It polymerize impedance；
According to the sub-synchronous oscillation characteristic of the polymerization impedance analysis system；
Wherein, the acquisition impedance model specifically includes:
Using Small Signal Model Analysis Method, the condition of small signal equation model of each component is initially set up, with the exit component i voltage
It is interface variables with electric current, the state equation model of each component linearisation can arrange as canonical form:Wherein, Δ indicates incremental computations；ΔxiIndicate that other states in addition to port voltage, electric current become
Measure increment column vector；Δui、ΔiiRespectively indicate port voltage, current increment column vector；aijIndicate the coefficient square of respective dimensions
Battle array；I, j ∈ I={ 1,2,3 }, subscript i indicate i-th of component；
By Laplace transform, algebraic equation model of the state equation model in the domain s is expressed as:Wherein, s indicates Laplace operator；
The exit component i incremental voltage Δ u is obtained by mathematical operationiWith increment current Δ iiBetween relationship are as follows: Δ ui(s)
=Zi(s)·Δii(s), wherein Zi(s) the equivalent impedance model of component i is indicated；
Method according to the sub-synchronous oscillation of the polymerization impedance analysis system includes: frequency-response analysis method, polymerization RLC series connection
Circuit analysis method, zero pole point stability methods；
When method according to the polymerization impedance analysis system sub-synchronous oscillation is frequency-response analysis method, the polymerization resistance is analyzed
Anti- real part imaginary part is with frequency variation curve；When it is described polymerization impedance imaginary part is there are zero crossing and zero crossing is located at subsynchronous frequency
When within the scope of rate, if corresponding real part is positive value, sub-synchronous oscillation has positive damping, sub-synchronous oscillation convergence, if right
Answering real part is zero or negative value, then sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging；
When method according to the polymerization impedance analysis system sub-synchronous oscillation is polymerization RLC series circuit analytic approach, it will polymerize
The model aggregation of impedance is an equivalent Second-Order RLC Filter Circuit series circuit, and when R is greater than 0, sub-synchronous oscillation has positive damping, secondary same
Step oscillation convergence, when R is less than or equal to 0, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging；
When method according to the polymerization impedance analysis system sub-synchronous oscillation is zero pole point stability methods, calculate described poly-
The zero pole point for closing impedance, when the zero pole point real part of the polymerization impedance is all negative, system is stablized；When the polymerization
When the zero pole point of impedance has positive real part, system is unstable.
2. sub-synchronous oscillation analysis method as described in claim 1, which is characterized in that when component i be one port network component,
The ZiIt (s) is single transmission function；
When component i is multiport network component, the ZiIt (s) is transfer function matrix.
3. a kind of sub-synchronous oscillation analytical equipment based on impedance network model characterized by comprising
Module is obtained, for obtaining impedance model；
Module is established, for establishing impedance network model by element interconnected method according to the impedance model；
Aggregation module is polymerize to obtain the resistance in conjunction with star triangular transformation to the impedance network model in series and parallel for using
The polymerization impedance of anti-network model；
Analysis module, for the sub-synchronous oscillation characteristic according to the polymerization impedance analysis system；
Wherein, the acquisition module is specifically used for:
Using Small Signal Model Analysis Method, the condition of small signal equation model of each component is initially set up, with the exit component i voltage
It is interface variables with electric current, the state equation model of each component linearisation can arrange as canonical form:Wherein, Δ indicates incremental computations；ΔxiIndicate that other states in addition to port voltage, electric current become
Measure increment column vector；Δui、ΔiiRespectively indicate port voltage, current increment column vector；aijIndicate the coefficient square of respective dimensions
Battle array；I, j ∈ I={ 1,2,3 }, subscript i indicate i-th of component；
By Laplace transform, algebraic equation model of the state equation model in the domain s is expressed as:Wherein, s indicates Laplace operator；
The exit component i incremental voltage Δ u is obtained by mathematical operationiWith increment current Δ iiBetween relationship are as follows: Δ ui(s)
=Zi(s)·Δii(s), wherein Zi(s) the equivalent impedance model of component i is indicated；
Method according to the sub-synchronous oscillation of the polymerization impedance analysis system includes: frequency-response analysis method, polymerization RLC series connection
Circuit analysis method, zero pole point stability methods；
When method according to the polymerization impedance analysis system sub-synchronous oscillation is frequency-response analysis method, the polymerization resistance is analyzed
Anti- real part imaginary part is with frequency variation curve；When it is described polymerization impedance imaginary part is there are zero crossing and zero crossing is located at subsynchronous frequency
When within the scope of rate, if corresponding real part is positive value, sub-synchronous oscillation has positive damping, sub-synchronous oscillation convergence, if right
Answering real part is zero or negative value, then sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging；
When method according to the polymerization impedance analysis system sub-synchronous oscillation is polymerization RLC series circuit analytic approach, it will polymerize
The model aggregation of impedance is an equivalent Second-Order RLC Filter Circuit series circuit, and when R is greater than 0, sub-synchronous oscillation has positive damping, secondary same
Step oscillation convergence, when R is less than or equal to 0, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or diverging；
When method according to the polymerization impedance analysis system sub-synchronous oscillation is zero pole point stability methods, calculate described poly-
The zero pole point for closing impedance, when the zero pole point real part of the polymerization impedance is all negative, system is stablized；When the polymerization
When the zero pole point of impedance has positive real part, system is unstable.
4. sub-synchronous oscillation analytical equipment as claimed in claim 3, which is characterized in that when component i be one port network component,
The ZiIt (s) is single transmission function；
When component i is multiport network component, the ZiIt (s) is transfer function matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610320261.0A CN106021682B (en) | 2016-05-13 | 2016-05-13 | Sub-synchronous oscillation analysis method and device based on impedance network model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610320261.0A CN106021682B (en) | 2016-05-13 | 2016-05-13 | Sub-synchronous oscillation analysis method and device based on impedance network model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106021682A CN106021682A (en) | 2016-10-12 |
| CN106021682B true CN106021682B (en) | 2019-08-09 |
Family
ID=57096888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610320261.0A Active CN106021682B (en) | 2016-05-13 | 2016-05-13 | Sub-synchronous oscillation analysis method and device based on impedance network model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106021682B (en) |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106599526B (en) * | 2017-02-17 | 2020-02-11 | 华北电力大学（保定） | Impedance analysis method for sub-synchronous oscillation of thermal power generating unit |
| CN107069811B (en) * | 2017-04-12 | 2019-07-26 | 清华大学 | Impedance network modeling and method for analyzing stability based on synchronous reference coordinate system |
| CN107221962B (en) * | 2017-08-02 | 2019-05-28 | 华北电力大学 | The subsynchronous stable region control method and controller of photovoltaic DC-to-AC converter control parameter |
| CN107769172B (en) | 2017-11-09 | 2019-03-01 | 清华大学 | The system-level protection system and method for secondary/supersynchronous resonance/oscillation |
| CN108154315B (en) * | 2018-01-19 | 2021-11-05 | 国电南瑞科技股份有限公司 | Grid-connected converter subsynchronous oscillation risk analysis method considering influence of phase-locked loop |
| CN108667008B (en) * | 2018-04-20 | 2020-07-14 | 上海交通大学 | Universal small signal modeling method for subsynchronous oscillation analysis of power system |
| CN108462192B (en) * | 2018-04-28 | 2020-05-05 | 清华大学 | Mode analysis method for broadband oscillation of power system |
| CN108448595B (en) * | 2018-05-15 | 2022-11-25 | 中国电力科学研究院有限公司 | Photovoltaic-series compensation system time domain simulation small disturbance analysis method and system |
| CN108667048B (en) * | 2018-05-31 | 2021-04-27 | 清华大学 | Frequency domain stability judging method and device for oscillation stability of new energy grid-connected system |
| CN109672217B (en) * | 2018-12-13 | 2022-03-01 | 华北电力大学 | Subsynchronous oscillation stability quantitative analysis method for wind turbine generator grid-connected system |
| CN109802406B (en) * | 2019-04-01 | 2022-05-17 | 云南电网有限责任公司 | Method for analyzing resonance stability of flexible direct-current transmission access system |
| CN110176762B (en) * | 2019-05-16 | 2022-08-09 | 国电南瑞科技股份有限公司 | Frequency domain impedance-based subsynchronous oscillation risk online evaluation method and device |
| CN110912125B (en) * | 2019-12-10 | 2022-02-08 | 东北电力大学 | High-proportion power electronic system resonance risk identification method |
| CN111146804A (en) * | 2020-01-02 | 2020-05-12 | 南方电网科学研究院有限责任公司 | Method and device for judging oscillation stability of wind power-flexible direct current transmission system |
| CN111327061B (en) * | 2020-02-12 | 2021-08-03 | 南方电网科学研究院有限责任公司 | Method and device for judging oscillation stability based on apparent impedance |
| CN111222297B (en) * | 2020-02-27 | 2022-11-08 | 华北电力大学 | Method for analyzing influence factors of subsynchronous oscillation of doubly-fed wind generator |
| CN112163532A (en) * | 2020-09-30 | 2021-01-01 | 国网冀北电力有限公司电力科学研究院 | Oscillation mode identification method and device of power system |
| CN113725910B (en) * | 2021-10-08 | 2022-05-03 | 南通大学 | Stability analysis and quantitative evaluation method for wind power plant grid-connected system |
| CN116822436B (en) * | 2023-06-30 | 2024-02-27 | 四川大学 | Oscillation risk sensitivity analysis method for direct-current transmission end alternating-current system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104978625A (en) * | 2015-06-10 | 2015-10-14 | 清华大学 | Convergence RLC circuit model-based method for analyzing subsynchronous resonance in power system |
| CN105281324A (en) * | 2015-10-10 | 2016-01-27 | 江苏省电力公司电力经济技术研究院 | Subsynchronous oscillation assessment method of electric power system |
-
2016
- 2016-05-13 CN CN201610320261.0A patent/CN106021682B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104978625A (en) * | 2015-06-10 | 2015-10-14 | 清华大学 | Convergence RLC circuit model-based method for analyzing subsynchronous resonance in power system |
| CN105281324A (en) * | 2015-10-10 | 2016-01-27 | 江苏省电力公司电力经济技术研究院 | Subsynchronous oscillation assessment method of electric power system |
Non-Patent Citations (2)
| Title |
|---|
| A Small-Signal Impedance Method for analyzing the SSR of Series-Compensated DFIG-Based Wind Farms;Huakun Liu 等;《2015 IEEE Power&Energy Society General Meeting》;20150730;第1-5页 * |
| 直驱风机风电场与交流电网相互作用引发次同步振荡的机理与特性分析;谢小荣 等;《中国电机工程学报》;20160505;第36卷(第9期);第2366-2372页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106021682A (en) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106021682B (en) | Sub-synchronous oscillation analysis method and device based on impedance network model | |
| CN107069811B (en) | Impedance network modeling and method for analyzing stability based on synchronous reference coordinate system | |
| Ambriz-Perez et al. | Advanced SVC models for Newton-Raphson load flow and Newton optimal power flow studies | |
| CN107664954A (en) | A kind of electromagnetic transient in power system visual simulating model automatic forming method | |
| CN102819641B (en) | Large-scale power distribution network integral model simplification method applicable to electromagnetic transient simulation | |
| CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
| CN104866665A (en) | Hybrid simulation method including power electronic equipment based on interface equivalence and interaction | |
| CN103049617B (en) | Method for simplifying passive electromagnetic transient-state simulation model of large-scale power distribution network | |
| CN105938578A (en) | Large-scale photovoltaic power station equivalent modeling method based on clustering analysis | |
| CN112149280B (en) | Method for obtaining power grid strength of multi-feed-in system of new energy containing SVG | |
| CN108054757A (en) | A kind of embedded idle and voltage N-1 Close loop security check methods | |
| CN103455948B (en) | A kind of distribution system multi-dimensional multi-resolution Modeling and the method for analysis | |
| CN104375876B (en) | Electromagnetical transient emulation method is immunized in a kind of 0+ errors under input quantity catastrophe | |
| CN101527455A (en) | Interconnected electric network distributed current calculating method on the basis of alternation and iteration of current module | |
| CN103077268A (en) | State space automatic modeling method orienting electromagnetic transient simulation of power system | |
| CN109638931A (en) | The power converter system model predictive control method and system of more DC-DC parallel connections | |
| CN106208099A (en) | A kind of Method for Reactive Power Optimization in Power based on bi-level programming and application thereof | |
| Lai et al. | Robust distributed cooperative control for DC mircogrids with time delays, noise disturbances, and switching topologies | |
| Wan et al. | Modelling, analysis and virtual parallel resistor damping control of VSC‐based DC grid using master–slave control mode | |
| CN109657276A (en) | Electromagnetical transient emulation method and system | |
| CN104700205B (en) | A kind of method for changing electricity grid network topological structure and selecting paralleling compensating device | |
| CN108347057B (en) | LCC-MMC mixed direct-current power grid alternating iteration load flow calculation method | |
| CN204229391U (en) | Hybrid simulation interface and there is the electric system hybrid real-time simulation equipment of this interface | |
| CN108808702A (en) | Time-lag power system electromechanic oscillation mode computational methods based on low order IGD-LMS algorithms | |
| CN110474329A (en) | A kind of continuous tide variable step size method and system based on polynomial regression technology |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |