CN106021682A - Subsynchronous oscillation analysis method and device based on impedance network model - Google Patents
Subsynchronous oscillation analysis method and device based on impedance network model Download PDFInfo
- Publication number
- CN106021682A CN106021682A CN201610320261.0A CN201610320261A CN106021682A CN 106021682 A CN106021682 A CN 106021682A CN 201610320261 A CN201610320261 A CN 201610320261A CN 106021682 A CN106021682 A CN 106021682A
- Authority
- CN
- China
- Prior art keywords
- impedance
- sub
- synchronous oscillation
- model
- polymerization
- 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
-
- 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
Abstract
The invention discloses a subsynchronous oscillation analysis method and device based on an impedance network model, wherein the method comprises the following steps of obtaining an impedance model; building the impedance network model by an element interconnection method according to the impedance network; aggregating the impedance network model by network equivalence methods such as series-parallel connection combined with star triangular transformation to obtain aggregation impedance of the impedance network model; and analyzing the subsynchronous oscillation features of a system according to the aggregation impedance. The method provides a new method for subsynchronous oscillation modeling and stable analysis; and the stability and the oscillation features of the system can be simply and efficiently judged.
Description
Technical field
The present invention relates to Power System Stability Analysis technical field, particularly relate to a kind of based on impedance network model secondary same
Step analysis of the oscillation method and apparatus.
Background technology
Under the trend that modern power systems power electronics degree is deepened day by day, electric electronic current change technology not only changes
The form of electrical network primary equipment, and its distinctive multiple dimensioned sequential control mechanism and low inertia broadband response characteristic, the most aobvious
Write the behavioral characteristics changing power system, bring new stability problem.In recent years, especially with the Flow varying power supply such as wind-powered electricity generation, photovoltaic
Novel sub-synchronous oscillation (Subsynchronous Oscillation, the SSO) problem that extensive access electrical network causes is the most prominent
Go out.Additionally, 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 the series/parallel resonant tank that current transformer control participation filter circuit and electrical network are constituted, all may swash
Send out sub-synchronous oscillation risk.It is new that the interaction of this magnanimity converter device (wind-powered electricity generation, photovoltaic, inverter etc.)-electrical network causes
Equipment safety, system stability and the power quality of type sub-synchronous oscillation serious threat modern power network, has become restriction wind in China
The bottleneck factor that the new forms of energy such as light are dissolved on a large scale.
It is temporary that the analysis method of the stability problem brought for power system power electronics continues to use traditional electromagnetism mostly
Impedance analysis method under state emulation, Eigenvalues analysis or reduced condition.The electromagnetic transient simulation software that is widely used at present (as
PSCAD, EMTP) only support modeling and the analysis of single step-length (fixed rate), it is difficult to take into account precision changeable to realize with efficiency
Time-domain-simulation under stream device, complex AC/DC network, magnanimity operating mode;Traditional Eigenvalues analysis method faces dimension calamity problem,
It is difficult to the actual bulk power grid being applicable to include thousands of current transformers, tens thousand of nodes.In recent years, it is widely used in power electronics to set
Standby and power system repercussion study impedance model analytic process provides new approaches.In actual applications, impedance model tool
Have a following advantage: 1) through impedance model and the total system impedance model that can obtain each subsystem of deriving, and physics meaning
Justice is relatively unambiguous;2), when changing systematic parameter, the impedance model of one or several subsystem is only affected, to overall impedance model shadow
Ring little;3) Nyquist stability criterion based on impedance model can be used to judge system stability, visual in image.Under reduced condition
Impedance analysis method in order to derive conveniently, the control strategy of system middle controller be have also been made corresponding simplification, have ignored portion
The dynamic characteristic of sub-controller.Although these simplify operation and are conducive to the foundation of system impedance model, bring and can not be ignored
Analytical error.
Summary of the invention
The purpose of the present invention is intended to solve one of above-mentioned technical problem the most to a certain extent.
To this end, the first of the present invention purpose is to propose a kind of sub-synchronous oscillation analysis side based on impedance network model
Method, this method provides sub-synchronous oscillation modeling and the new method of stability analysis, can the stability of simply efficient 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.
For reaching above-mentioned purpose, embodiment proposition is a kind of based on impedance network model secondary same according to a first aspect of the present invention
Step analysis of the oscillation method, comprises the following steps: obtain impedance model;Set up by element interconnected method according to described impedance model
Impedance network model;Use connection in series-parallel to combine the Equivalent Network methods such as star triangular transformation and described impedance network model is carried out abbreviation
Obtain the polymerization impedance of described impedance network model;Sub-synchronous oscillation characteristic according to described polymerization impedance analysis system.
The sub-synchronous oscillation based on impedance network model of the embodiment of the present invention analyzes method, first obtains impedance model,
Set up impedance network model then according to impedance model by element interconnected method, then use connection in series-parallel to combine star triangular transformation etc.
Impedance network model is polymerized the polymerization impedance obtaining impedance network model by Equivalent Network method, finally according to polymerization impedance
The sub-synchronous oscillation characteristic of analysis system.This method provide the new method of sub-synchronous oscillation modeling and stability analysis, Ke Yijian
The stability of single efficiently decision-making system and oscillating characteristic.
It addition, sub-synchronous oscillation based on impedance network model according to the above embodiment of the present invention analysis method is all right
There is following additional technical characteristic:
In some instances, described acquisition impedance model specifically includes: use Small Signal Model Analysis Method, initially sets up each
The condition of small signal equation model of assembly, with assembly i exit voltage and current as interface variables, the linearizing shape of each assembly
State equation model can arrange as canonical form:Wherein, Δ represents incremental computations;ΔxiExpression removes
Other state variable increment column vectors outside port voltage, electric current;Δui、ΔiiRepresent respectively port voltage, current increment arrange to
Amount;aijRepresent the coefficient matrix of respective dimensions;I, j ∈ I={1,2,3}, subscript i represents i-th assembly;Become by Laplce
Changing, described state equation model algebraic equation model representation in s territory is:
Wherein, s represents Laplace operator;Assembly i exit voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween pass
System, such as following formula: Δ ui(s)=Zi(s)·Δii(s), wherein, ZiS () represents the equiva lent impedance model of assembly i.
In some instances, when assembly i is one port network assembly, described ZiS () is single transmission function;As assembly i
For multiport network assembly, described ZiS () is transfer function matrix.
In some instances, include according to the method for the sub-synchronous oscillation of described polymerization impedance analysis system: frequency characteristic
Analytic process, polymerization RLC series circuit analytic process, zero pole point stability methods.
In some instances, the method according to described polymerization impedance analysis system sub-synchronous oscillation is frequency-response analysis method
Time, analyze the real part imaginary part of described polymerization impedance with frequency variation curve, when the imaginary part of described polymerization impedance exist zero crossing and
When zero crossing is positioned at subsynchronous frequency range, if the real part of correspondence be on the occasion of, then sub-synchronous oscillation has a positive damping, secondary with
Step vibration convergence, if corresponding real part is zero or negative value, then sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation is held
Continue or dissipate;When method according to described polymerization impedance analysis system sub-synchronous oscillation is for polymerization RLC series circuit analytic process, will
The model aggregation of polymerization impedance is an equivalent Second-Order RLC Filter Circuit series circuit, and when R is more than 0, sub-synchronous oscillation has positive damping,
Sub-synchronous oscillation is stable, 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 sends out
Dissipate;When the method for the sub-synchronous oscillation according to described polymerization impedance analysis system is zero pole point stability methods, calculate described
The zero pole point of polymerization impedance, the zero pole point real part of described polymerization impedance that and if only if all negative time, system stability;When described poly-
When the zero pole point of conjunction impedance has positive real part, system is unstable.
For reaching above-mentioned purpose, a kind of based on impedance network model time of embodiment proposition according to the second aspect of the invention
Synchronized oscillation analytical equipment, including: acquisition module, it is used for obtaining impedance model;Set up module, for according to described impedance model
Impedance network model is set up by element interconnected method;Aggregation module, is used for using connection in series-parallel to combine the networks such as star triangular transformation
Described impedance network model is polymerized the polymerization impedance obtaining described impedance network model by equivalence method;Analyze module, use
In the sub-synchronous oscillation characteristic according to described 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 and sets up impedance network model according to impedance model by element interconnected method, and repolymerization module is adopted
Combine the Equivalent Network methods such as star triangular transformation impedance network model is polymerized to obtain impedance network model with connection in series-parallel
Polymerization impedance, ultimate analysis module is according to the sub-synchronous oscillation characteristic of polymerization impedance analysis system.This method provide subsynchronous
Vibration modeling and the new method of stability analysis, can the simply efficiently stability of decision-making system and oscillating characteristic.
It addition, sub-synchronous oscillation analytical equipment based on impedance network model according to the above embodiment of the present invention is all right
There is following additional technical characteristic:
In some instances, described acquisition module specifically for: use Small Signal Model Analysis Method, initially set up each assembly
Condition of small signal equation model, with assembly i exit voltage and current as interface variables, the linearizing state side of each assembly
Journey model can arrange as canonical form:Wherein, Δ represents incremental computations;ΔxiRepresent except port
Other state variable increment column vectors outside voltage, electric current;Δui、ΔiiRepresent port voltage, current increment column vector respectively;
aijRepresent the coefficient matrix of respective dimensions;I, j ∈ I={1,2,3}, subscript i represents i-th assembly;By Laplace transform,
Described state equation model algebraic equation model representation in s territory is:Its
In, s represents Laplace operator;Assembly i exit voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relation,
Such as following formula: Δ ui(s)=Zi(s)·Δii(s), wherein, ZiS () represents the equiva lent impedance model of assembly i.
In some instances, when assembly i is one port network assembly, described ZiS () is single transmission function;As assembly i
For multiport network assembly, described ZiS () is transfer function matrix.
In some instances, the analysis method of described analysis module includes: frequency-response analysis method, polymerization RLC series electrical
Road analytic process, zero pole point stability methods.
In some instances, when the analysis method of described analysis module is frequency-response analysis method, analyze described polymerization resistance
Anti-real part imaginary part is with frequency variation curve, when the imaginary part of described polymerization impedance exists zero crossing and zero crossing is positioned at subsynchronous frequency
Time in the range of rate, if the real part of correspondence be on the occasion of, then sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stable, 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 dissipates;Described analysis
When the analysis method of module is for polymerization RLC series circuit analytic process, it is an equivalent Second-Order RLC Filter Circuit by the model aggregation of polymerization impedance
Series circuit, when R is more than 0, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is restrained, when R is less than or equal to 0, subsynchronous
Vibration has zero damping or negative damping, and sub-synchronous oscillation continues or dissipates;The analysis method of described analysis module is that zero pole point is steady
During way of qualitative analysis, calculating the zero pole point of described polymerization impedance, the zero pole point real part of described polymerization impedance that and if only if is all
Time negative, system stability;When the zero pole point of described polymerization impedance has positive real part, system is unstable.
Aspect and advantage that the present invention adds will part be given in the following description, and part will become from the following description
Obtain substantially, or recognized by the practice of the present invention.
Accompanying drawing explanation
Above-mentioned and/or the additional aspect of the present invention and advantage are from combining the accompanying drawings below description to embodiment and will become
Substantially with easy to understand, wherein:
Fig. 1 is the flow process that based on impedance network model according to an embodiment of the invention sub-synchronous oscillation analyzes method
Figure;
Fig. 2 is the schematic diagram of power electronics power 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, the void of the polymerization impedance of power electronics power system according to an embodiment of the invention
Portion)-frequency characteristic;
Fig. 6 is the structure of based on impedance network model according to an embodiment of the invention sub-synchronous oscillation analytical equipment
Schematic diagram.
Detailed description of the invention
Embodiments of the invention are described below in detail, and the example of described embodiment is shown in the drawings, the most from start to finish
Same or similar label represents same or similar element or has the element of same or like function.Below with reference to attached
The embodiment that figure describes is exemplary, it is intended to is used for explaining the present invention, and is not considered as limiting the invention.
The present invention proposes the sub-synchronous oscillation based on impedance network model of power electronics power system and analyzes method
And device, input-output impedance model will be equivalent to by converter device, be impedance network with system simultaneous afterwards, and will resistance
Anti-network is progressively polymerized to an impedance, by analyzing the frequency characteristic of impedance, it is achieved 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 based on impedance network model analyze method and
Device.
Fig. 1 is the flow process that based on impedance network model according to an embodiment of the invention sub-synchronous oscillation analyzes method
Figure.
Include following as it is shown in figure 1, the sub-synchronous oscillation based on impedance network model of the embodiment of the present invention analyzes method
Step:
Step 101, obtains impedance model.
Specifically, power electronic equipment and other power system component transient Model and the Controlling model of full size, analyze electricity
The Inner Constitution of each assembly and annexation in power electronization power system, for assemblies various in power electronics power system, s
Under territory, the input voltage of N-port assembly and the relation of input current are as follows:
Transferring Vector-Matrix Form to is: V (s)=Z (s) I (s).
Wherein, Vi(s)、IiS () is i-th port voltage and the port current of N end assembly;ZijS () is impedance matrix pair
Answer the impedance between port, there is following form:
Wherein, wherein: Z (s) is impedance matrix, Z (s) depend on component type, the annexation of internal component and
Control system.For one port network assembly, Z (s) is single transmission functional form;For multiport network assembly, Z (s) is
The matrix of transmission function.
It should be noted that for the impedance model obtaining each assembly of power electronics power system, can use but not
It is limited to Small Signal Model Analysis Method, Model Distinguish method.
For example, use Small Signal Model Analysis Method, initially set up the condition of small signal equation model of each assembly, with group
Part i exit voltage and current is interface variables, and the linearizing state equation model of each assembly can arrange as canonical form:Wherein, Δ represents incremental computations;ΔxiRepresent that other states in addition to port voltage, electric current become
Amount increment column vector;Δui、ΔiiRepresent port voltage, current increment column vector respectively;aijRepresent the coefficient square of respective dimensions
Battle array;I, j ∈ I={1,2,3}, subscript i represents i-th assembly.
By Laplace transform, described state equation model algebraic equation model representation in s territory is:Wherein, s represents Laplace operator.
Assembly i exit voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relation, such as following formula: Δ ui(s)
=Zi(s)·Δii(s), wherein, ZiS () represents the equiva lent impedance model of assembly i.
Step 102, sets up impedance network model according to impedance model by element interconnected method.
Specifically, describe in detail in conjunction with Fig. 2, as in figure 2 it is shown, with certain large-scale AC and DC power containing magnanimity converter device
As a example by system, under specific operation, carry out Load flow calculation, obtain each component ports electric parameter, obtain according to aforementioned modeling method
After taking the impedance of each assembly, use element interconnected method (Component Connection Method, CCM) constructing system
Overall impedance network model.I.e. on the basis of original power system topological structure, each assembly is used its impedance model
Replace, have between the component impedance of electric connecting relation and use ideal line to connect, thus each impedance simultaneous is got up i.e.
Available impedance network model as shown in Figure 3.
In modeling process, equivalence can be carried out according to certain principle, as consistent in wind power plant's blower fan type and operating mode
It is equivalent that unit can use single Large Copacity equivalence unit, and with simplified operation, but its blower interior control system of general tree is entered
Row simplifies.
Step 103, uses connection in series-parallel to combine the Equivalent Network methods such as star triangular transformation and is polymerized impedance network model
Obtain the polymerization impedance of impedance network model.
Specifically, after obtaining power electronics power system impedance network model, impedance network is carried out abbreviation.Abbreviation
Time, the method for connection in series-parallel knot and the conversion of star triangular network according to the annexation between impedor, can be used to carry out abbreviation;
Abbreviation can also be carried out by the transformation calculations between impedance matrix.Finally obtain the polymerization impedance of whole network, its expression formula one
As be more complicated transmission functional form.
Illustrate: Fig. 3 show the Fig. 2 large-scale AC and DC power system impedance network model containing multiple-variable flow device device,
Polymerization impedance Z (s) of whole system is obtained after abbreviation process shown in Fig. 4.First by bulk power grid 1 impedance and HVDC impedance string
Connection, the impedance of photovoltaic power plant is connected with circuit 1 impedance, wind field impedance is connected with circuit 2 impedance, current transformer formula load impedance and circuit
6 impedance series connection, i.e.
ImpedanceAnd impedanceConnect, i.e. with circuit 3 impedance after parallel connection
ImpedanceImpedanceWith impedanceParallel connection, i.e.
Circuit 4 impedance, circuit 5 impedance carry out delta-Y transformation, i.e. with tie line impedance
ImpedanceWith impedance ZY1(s) series connection, impedance ZY2(s) and bulk power grid 2 impedance ZS2(s) series connection, impedance ZS2(s)
With bulk power grid 3 impedance ZS3S () is connected, i.e.
ImpedanceWith impedanceParallel connection, i.e.
Impedance ZY(s) and impedanceSystematic polymerization impedance Z (s) is obtained, i.e. after series connection
Step 104, according to the sub-synchronous oscillation characteristic of polymerization impedance analysis system.
Specifically, after the equiva lent impedance obtaining power electronics power system network, multiple method can be used to carry out
Stability and oscillating characteristic analysis, typical such as frequency-response analysis method, to be polymerized RLC series circuit analytic process, zero pole point stable
Property analytic process, (broad sense) Nyquist stability methods etc..
Such as, zero pole point method for analyzing stability.After obtaining the transmission function of systematic polymerization impedance Z (s), zero can be taked
Limit stability methods judges the stability of system, i.e. calculates the zero pole point of systematic polymerization impedance Z (s), is polymerized by judgement
The real part of impedance Z (s) zero pole point and imaginary part feature judge the oscillating characteristic of system.
And if only if polymerization impedance Z (s) zero pole point real part all negative time, system stability;When polymerization impedance Z (s)
When zero pole point has positive real part, system is unstable.Especially, when the imaginary part pair of polymerization impedance Z (s) certain or certain several zero pole point
When answering resonant frequency to be positioned at subsynchronous frequency range, there is the probability that SSO occurs in system, if now corresponding zero pole point real part
Being negative, SSO restrains;If corresponding zero pole point real part non-negative, exist and continue or dissipate SSO risk.
Such as, frequency-response analysis method.Power electronics power system frequency domain polymerization impedance Z (s) is carried out frequency characteristic
Analyze, polymerization impedance Z (s) real part, imaginary part can be used with frequency variation characteristics analysis.
Analyze polymerization impedance Z (s) real part imaginary part with frequency variation curve, when polymerization impedance Z (s) imaginary part exist zero crossing and
When zero crossing is positioned at subsynchronous frequency range, if now the real part of its correspondence be on the occasion of, show that SSO has positive damping, SSO is steady
Fixed, if corresponding real part is zero or negative value, show that SSO has zero damping or a negative damping, SSO persistent oscillation or dissipate.
It is illustrated in figure 5 under certain power electronics power system operating mode polymerization impedance real part imaginary part with frequency change song
Line, it can be seen that imaginary part curve zero crossing respective frequencies is positioned at subsynchronous frequency range, and under respective frequencies, real part is negative
Value, so SSO has negative damping under operating mode, SSO is unstable.
Such as, polymerization RLC series circuit analytic process.Bent by analyzing systematic polymerization impedance Z (s) real part imaginary part frequency characteristic
After line, available system resonance frequencies ωr.Near series resonance frequency ω | 0≤| ω-ωr| < (h is the least normal to h}
Number), equiva lent 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 intending by nonlinear parameter in resonant frequency vertex neighborhood
Close optimization method (as used method of least square) to be calculated, be equivalent to polymerization impedance model further be polymerized rlc circuit.
Based on the polymerization second-order circuit parameter obtained, damping and frequency of oscillation (the i.e. SSO damping of second-order circuit can be calculated
And frequency), and then carry out quantization SSO analysis, its computing formula is as follows:
Thus, as R > 0 time, SSO has positive damping, and SSO is stable;Otherwise, SSO has zero damping or negative damping, and SSO is lasting
Vibrate or dissipate.
Method is analyzed corresponding with the sub-synchronous oscillation based on impedance network model that above-described embodiment provides, the present invention's
A kind of embodiment also provides for a kind of sub-synchronous oscillation analytical equipment based on impedance network model, owing to the embodiment of the present invention provides
Sub-synchronous oscillation based on impedance network model analyze that method and above-mentioned several embodiments provide based on impedance network model
Sub-synchronous oscillation analyze method there is same or analogous technical characteristic, therefore aforementioned based on impedance network model time with
The sub-synchronous oscillation based on impedance network model that the embodiment of step analysis of the oscillation method is also applied for the present embodiment and provides divides
Analysis apparatus, is not described in detail in the present embodiment.As shown in Figure 6, should sub-synchronous oscillation analysis based on impedance network model
Device comprises the steps that acquisition module 10, sets up module 20, aggregation module 30 and analyze module 40.
Wherein, acquisition module 10 is used for obtaining impedance model.
Set up module 20 for setting up impedance network model according to impedance model by element interconnected method.
Aggregation module 30 is used for using connection in series-parallel to combine the Equivalent Network methods such as star triangular transformation and enters impedance network model
Row polymerization obtains the polymerization impedance of impedance network model.
Analyze module 40 for the sub-synchronous oscillation characteristic according to polymerization impedance analysis system.
In some instances, acquisition module 10 specifically for: use Small Signal Model Analysis Method, initially set up each assembly
Condition of small signal equation model, with assembly i exit voltage and current as interface variables, the linearizing state equation of each assembly
Model can arrange as canonical form:Wherein, Δ represents incremental computations;ΔxiRepresent except port electricity
Other state variable increment column vectors outside pressure, electric current;Δui、ΔiiRepresent port voltage, current increment column vector respectively;aij
Represent the coefficient matrix of respective dimensions;I, j ∈ I={1,2,3}, subscript i represents i-th assembly.
By Laplace transform, state equation model algebraic equation model representation in s territory is:Wherein, s represents Laplace operator.
Assembly i exit voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relation, such as following formula: Δ ui(s)
=Zi(s)·Δii(s), wherein, ZiS () represents the equiva lent impedance model of assembly i.
In some instances, when assembly i is one port network assembly, described ZiS () is single transmission function.As assembly i
For multiport network assembly, described ZiS () is transfer function matrix.
In some instances, the analysis method analyzing module 40 includes: frequency-response analysis method, polymerization RLC series circuit
Analytic process, zero pole point stability methods.
In some instances, when the analysis method of analysis module 40 is frequency-response analysis method, the reality of polymerization impedance is analyzed
Portion's imaginary part is with frequency variation curve, when the imaginary part of polymerization impedance exists zero crossing and zero crossing is positioned at subsynchronous frequency range
Time, if the real part of correspondence be on the occasion of, then sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stable, 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 dissipates.
When the analysis method of analysis module 40 is for polymerization RLC series circuit analytic process, by the model aggregation of polymerization impedance it is
One equivalent Second-Order RLC Filter Circuit series circuit, when R is more than 0, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stable, when R is little
When equal to 0, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or dissipates.
Analyze the analysis method of module when being zero pole point stability methods, calculate the zero pole point of polymerization impedance, when and only
When be polymerized impedance zero pole point real part all negative time, system stability;When the zero pole point being polymerized impedance has positive real part, system
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 power system of son equipment, especially for wind, the light generating high power system of permeability, it is provided that a kind of sub-synchronous oscillation is built
Mould and the new method of stability analysis.Impedance model is a kind of external characteristics model, for because various reasons are (as manufacturing firm's secrecy needs
Want) electrical circuitry equipment of internal control system model and parameter cannot be provided for, it is provided that a kind of accurately without regard in equipment
The expression method of portion's private information, and exponent number and the parameter of model can be determined by discrimination method.Owing to impedance is easy to
Use simple circuit transformations (such as connection in series-parallel, star-triangular transformation) to realize abbreviation and polymerization, thus overcome traditional little letter
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
Limit, impedance frequency characteristics etc. quantify stability and the oscillating characteristic of decision-making system, have physical meaning clear and definite, the highest
The advantage 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 and sets up impedance network model according to impedance model by element interconnected method, and repolymerization module is adopted
Combine the Equivalent Network methods such as star triangular transformation impedance network model is polymerized to obtain impedance network model with connection in series-parallel
Polymerization impedance, ultimate analysis module is according to the sub-synchronous oscillation of polymerization impedance analysis system.This method provide sub-synchronous oscillation
Modeling and the new method of stability analysis, can the simply efficiently stability of decision-making system and oscillating characteristic.
In describing the invention, it is to be understood that term " first ", " second " are only used for describing purpose, and can not
It is interpreted as instruction or hint relative importance or the implicit quantity indicating indicated technical characteristic.Thus, define " the
One ", the feature of " second " can express or implicitly include at least one this feature.In describing the invention, " multiple "
It is meant that at least two, such as two, three etc., unless otherwise expressly limited specifically.
In the description of this specification, reference term " embodiment ", " some embodiments ", " example ", " specifically show
Example " or the description of " some examples " etc. means to combine this embodiment or example describes specific features, structure, material or spy
Point is contained at least one embodiment or the example of the present invention.In this manual, to the schematic representation of above-mentioned term not
Identical embodiment or example must be directed to.And, the specific features of description, structure, material or feature can be in office
One or more embodiments or example combine in an appropriate manner.Additionally, in the case of the most conflicting, the skill of this area
The feature of the different embodiments described in this specification or example and different embodiment or example can be tied by art personnel
Close and combination.
In flow chart or at this, any process described otherwise above or method description are construed as, and expression includes
One or more is for realizing the module of code, fragment or the portion of the executable instruction of the step of specific logical function or process
Point, and the scope of the preferred embodiment of the present invention includes other realization, wherein can not by shown or discuss suitable
Sequence, including according to involved function by basic mode simultaneously or in the opposite order, performs function, and this should be by the present invention
Embodiment person of ordinary skill in the field understood.
Those skilled in the art are appreciated that and realize all or part of step that above-described embodiment method is carried
Suddenly the program that can be by completes to instruct relevant hardware, and described program can be stored in a kind of computer-readable storage medium
In matter, this program upon execution, including one or a combination set of the step of embodiment of the method.
Although above it has been shown and described that embodiments of the invention, it is to be understood that above-described embodiment is example
Property, it is impossible to being interpreted as limitation of the present invention, those of ordinary skill in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, revises, replaces and modification.
Claims (10)
1. a sub-synchronous oscillation based on impedance network model analyzes method, it is characterised in that comprise the following steps:
Obtain impedance model;
Impedance network model is set up by element interconnected method according to described impedance model;
Use connection in series-parallel to combine the Equivalent Network methods such as star triangular transformation described impedance network model carries out polymerization to obtain described
The polymerization impedance of impedance network model;
Sub-synchronous oscillation characteristic according to described polymerization impedance analysis system.
2. sub-synchronous oscillation as claimed in claim 1 analyzes method, it is characterised in that described acquisition impedance model specifically wraps
Include:
Use Small Signal Model Analysis Method, initially set up the condition of small signal equation model of each assembly, with assembly i exit voltage
Being interface variables with electric current, the linearizing state equation model of each assembly can arrange as canonical form:Wherein, Δ represents incremental computations;ΔxiRepresent that other states in addition to port voltage, electric current become
Amount increment column vector;Δui、ΔiiRepresent port voltage, current increment column vector respectively;aijRepresent the coefficient square of respective dimensions
Battle array;I, j ∈ I={1,2,3}, subscript i represents i-th assembly;
By Laplace transform, described state equation model algebraic equation model representation in s territory is:Wherein, s represents Laplace operator;
Assembly i exit voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relation, such as following formula: Δ ui(s)=Zi
(s)·Δii(s), wherein, ZiS () represents the equiva lent impedance model of assembly i.
3. sub-synchronous oscillation as claimed in claim 2 analyzes method, it is characterised in that when assembly i is one port network assembly,
Described ZiS () is single transmission function;
When assembly i is multiport network assembly, described ZiS () is transfer function matrix.
4. sub-synchronous oscillation as claimed in claim 1 analyzes method, it is characterised in that according to described polymerization impedance analysis system
The method of sub-synchronous oscillation include: frequency-response analysis method, polymerization RLC series circuit analytic process, zero pole point stability analysis
Method.
5. sub-synchronous oscillation as claimed in claim 4 analyzes method, it is characterised in that
When method according to described polymerization impedance analysis system sub-synchronous oscillation is frequency-response analysis method, analyze described polymerization resistance
Anti-real part imaginary part is with frequency variation curve, when the imaginary part of described polymerization impedance exists zero crossing and zero crossing is positioned at subsynchronous frequency
Time in the range of rate, if the real part of correspondence be on the occasion of, then sub-synchronous oscillation has a positive damping, and sub-synchronous oscillation is restrained, 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 dissipates;
When method according to described polymerization impedance analysis system sub-synchronous oscillation is for polymerization RLC series circuit analytic process, will polymerization
The model aggregation of impedance is an equivalent Second-Order RLC Filter Circuit series circuit, and when R is more than 0, sub-synchronous oscillation has positive damping, secondary same
Step vibration 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 dissipates;
When method according to described polymerization impedance analysis system sub-synchronous oscillation is zero pole point stability methods, calculate described poly-
Close the zero pole point of impedance, the zero pole point real part of described polymerization impedance that and if only if all negative time, system stability;When described polymerization
When the zero pole point of impedance has positive real part, system is unstable.
6. a sub-synchronous oscillation analytical equipment based on impedance network model, it is characterised in that including:
Acquisition module, is used for obtaining impedance model;
Set up module, for setting up impedance network model according to described impedance model by element interconnected method;
Aggregation module, is used for using connection in series-parallel to combine the Equivalent Network methods such as star triangular transformation and carries out described impedance network model
Polymerization obtains the polymerization impedance of described impedance network model;
Analyze module, for the sub-synchronous oscillation characteristic according to described polymerization impedance analysis system.
7. sub-synchronous oscillation analytical equipment as claimed in claim 6, it is characterised in that described acquisition module specifically for:
Use Small Signal Model Analysis Method, initially set up the condition of small signal equation model of each assembly, with assembly i exit voltage
Being interface variables with electric current, the linearizing state equation model of each assembly can arrange as canonical form:Wherein, Δ represents incremental computations;ΔxiRepresent that other states in addition to port voltage, electric current become
Amount increment column vector;Δui、ΔiiRepresent port voltage, current increment column vector respectively;aijRepresent the coefficient square of respective dimensions
Battle array;I, j ∈ I={1,2,3}, subscript i represents i-th assembly;
By Laplace transform, described state equation model algebraic equation model representation in s territory is:Wherein, s represents Laplace operator;
Assembly i exit voltage Δ u is obtained by mathematical operationiWith electric current Δ iiBetween relation, such as following formula: Δ ui(s)=Zi
(s)·Δii(s), wherein, ZiS () represents the equiva lent impedance model of assembly i.
8. sub-synchronous oscillation analytical equipment as claimed in claim 7, it is characterised in that when assembly i is one port network assembly,
Described ZiS () is single transmission function;
When assembly i is multiport network assembly, described ZiS () is transfer function matrix.
9. sub-synchronous oscillation analytical equipment as claimed in claim 6, it is characterised in that the analysis method bag of described analysis module
Include: frequency-response analysis method, polymerization RLC series circuit analytic process, zero pole point stability methods.
10. sub-synchronous oscillation analytical equipment as claimed in claim 9, it is characterised in that
When the analysis method of described analysis module is frequency-response analysis method, analyze the real part imaginary part of described polymerization impedance with frequency
Change curve, when there is zero crossing in the imaginary part of described polymerization impedance and zero crossing is positioned at subsynchronous frequency range, if right
The real part answered on the occasion of, then sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stable, 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 dissipates;
When the analysis method of described analysis module is for polymerization RLC series circuit analytic process, it is one by the model aggregation of polymerization impedance
Individual equivalence Second-Order RLC Filter Circuit series circuit, when R is more than 0, sub-synchronous oscillation has positive damping, and sub-synchronous oscillation is stable, when R is less than
During equal to 0, sub-synchronous oscillation has zero damping or negative damping, and sub-synchronous oscillation continues or dissipates;
When the analysis method of described analysis module is zero pole point stability methods, calculate the zero pole point of described polymerization impedance, when
And if only if the zero pole point real part of described polymerization impedance all negative time, system stability;When the zero pole point of described polymerization impedance has
During positive real part, system is unstable.
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 true CN106021682A (en) | 2016-10-12 |
CN106021682B 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) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599526A (en) * | 2017-02-17 | 2017-04-26 | 华北电力大学(保定) | Impedance analyzing method for sub-synchronous oscillation of thermal power unit |
CN107069811A (en) * | 2017-04-12 | 2017-08-18 | 清华大学 | Impedance network modeling and method for analyzing stability based on reference synchronization coordinate system |
CN107221962A (en) * | 2017-08-02 | 2017-09-29 | 华北电力大学 | The subsynchronous stable region control method and controller of photovoltaic DC-to-AC converter control parameter |
CN107769172A (en) * | 2017-11-09 | 2018-03-06 | 清华大学 | The system-level protection system and method for secondary/supersynchronous resonance/vibration |
CN108154315A (en) * | 2018-01-19 | 2018-06-12 | 国电南瑞科技股份有限公司 | A kind of grid-connected converter sub-synchronous oscillation risk analysis method for considering phaselocked loop and influencing |
CN108448595A (en) * | 2018-05-15 | 2018-08-24 | 中国电力科学研究院有限公司 | A kind of photovoltaic-string complement system time-domain-simulation Small Perturbation Analysis and system |
CN108462192A (en) * | 2018-04-28 | 2018-08-28 | 清华大学 | The pattern analysis method of electric system broadband oscillation |
CN108667008A (en) * | 2018-04-20 | 2018-10-16 | 上海交通大学 | A kind of general Approach for Modeling of Small-Signal for subsynchronous oscillation of electrical power system analysis |
CN108667048A (en) * | 2018-05-31 | 2018-10-16 | 清华大学 | The frequency domain of new energy grid connection system oscillatory stability sentences steady method and device |
CN109672217A (en) * | 2018-12-13 | 2019-04-23 | 华北电力大学 | Wind turbines grid-connected system sub-synchronous oscillation Convenience stable criterion method |
CN109802406A (en) * | 2019-04-01 | 2019-05-24 | 云南电网有限责任公司 | A method of analysis flexible DC transmission access system resonance stability |
CN110176762A (en) * | 2019-05-16 | 2019-08-27 | 国电南瑞科技股份有限公司 | A kind of sub-synchronous oscillation risk online evaluation method and device based on frequency domain impedance |
CN110912125A (en) * | 2019-12-10 | 2020-03-24 | 东北电力大学 | 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 |
CN111222297A (en) * | 2020-02-27 | 2020-06-02 | 华北电力大学 | Method for analyzing influence factors of subsynchronous oscillation of doubly-fed wind generator |
CN111327061A (en) * | 2020-02-12 | 2020-06-23 | 南方电网科学研究院有限责任公司 | Method and device for judging oscillation stability based on apparent impedance |
CN112163532A (en) * | 2020-09-30 | 2021-01-01 | 国网冀北电力有限公司电力科学研究院 | Oscillation mode identification method and device of power system |
CN113725910A (en) * | 2021-10-08 | 2021-11-30 | 南通大学 | Stability analysis and quantitative evaluation method for wind power plant grid-connected system |
CN116822436A (en) * | 2023-06-30 | 2023-09-29 | 四川大学 | 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 |
---|
HUAKUN LIU 等: "A Small-Signal Impedance Method for analyzing the SSR of Series-Compensated DFIG-Based Wind Farms", 《2015 IEEE POWER&ENERGY SOCIETY GENERAL MEETING》 * |
谢小荣 等: "直驱风机风电场与交流电网相互作用引发次同步振荡的机理与特性分析", 《中国电机工程学报》 * |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599526A (en) * | 2017-02-17 | 2017-04-26 | 华北电力大学(保定) | Impedance analyzing method for sub-synchronous oscillation of thermal power unit |
CN107069811A (en) * | 2017-04-12 | 2017-08-18 | 清华大学 | Impedance network modeling and method for analyzing stability based on reference synchronization 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 |
CN107221962A (en) * | 2017-08-02 | 2017-09-29 | 华北电力大学 | The subsynchronous stable region control method and controller of photovoltaic DC-to-AC converter control parameter |
CN107769172A (en) * | 2017-11-09 | 2018-03-06 | 清华大学 | The system-level protection system and method for secondary/supersynchronous resonance/vibration |
US11476670B2 (en) | 2017-11-09 | 2022-10-18 | Tsinghua University | System-level protection system and method for sub/super-synchronous resonance/oscillation |
CN108154315A (en) * | 2018-01-19 | 2018-06-12 | 国电南瑞科技股份有限公司 | A kind of grid-connected converter sub-synchronous oscillation risk analysis method for considering phaselocked loop and influencing |
CN108154315B (en) * | 2018-01-19 | 2021-11-05 | 国电南瑞科技股份有限公司 | Grid-connected converter subsynchronous oscillation risk analysis method considering influence of phase-locked loop |
CN108667008A (en) * | 2018-04-20 | 2018-10-16 | 上海交通大学 | A kind of general Approach for Modeling of Small-Signal for subsynchronous oscillation of electrical power system analysis |
CN108667008B (en) * | 2018-04-20 | 2020-07-14 | 上海交通大学 | Universal small signal modeling method for subsynchronous oscillation analysis of power system |
CN108462192A (en) * | 2018-04-28 | 2018-08-28 | 清华大学 | The pattern analysis method of electric system broadband oscillation |
CN108448595B (en) * | 2018-05-15 | 2022-11-25 | 中国电力科学研究院有限公司 | Photovoltaic-series compensation system time domain simulation small disturbance analysis method and system |
CN108448595A (en) * | 2018-05-15 | 2018-08-24 | 中国电力科学研究院有限公司 | A kind of photovoltaic-string complement system time-domain-simulation Small Perturbation Analysis and system |
CN108667048A (en) * | 2018-05-31 | 2018-10-16 | 清华大学 | The frequency domain of new energy grid connection system oscillatory stability sentences steady method and device |
CN108667048B (en) * | 2018-05-31 | 2021-04-27 | 清华大学 | Frequency domain stability judging method and device for oscillation stability of new energy grid-connected system |
CN109672217A (en) * | 2018-12-13 | 2019-04-23 | 华北电力大学 | Wind turbines grid-connected system sub-synchronous oscillation Convenience stable criterion method |
CN109802406B (en) * | 2019-04-01 | 2022-05-17 | 云南电网有限责任公司 | Method for analyzing resonance stability of flexible direct-current transmission access system |
CN109802406A (en) * | 2019-04-01 | 2019-05-24 | 云南电网有限责任公司 | A method of analysis flexible DC transmission access system resonance stability |
CN110176762B (en) * | 2019-05-16 | 2022-08-09 | 国电南瑞科技股份有限公司 | Frequency domain impedance-based subsynchronous oscillation risk online evaluation method and device |
CN110176762A (en) * | 2019-05-16 | 2019-08-27 | 国电南瑞科技股份有限公司 | A kind of sub-synchronous oscillation risk online evaluation method and device based on frequency domain impedance |
CN110912125B (en) * | 2019-12-10 | 2022-02-08 | 东北电力大学 | High-proportion power electronic system resonance risk identification method |
CN110912125A (en) * | 2019-12-10 | 2020-03-24 | 东北电力大学 | 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 |
CN111327061A (en) * | 2020-02-12 | 2020-06-23 | 南方电网科学研究院有限责任公司 | Method and device for judging oscillation stability based on apparent impedance |
CN111327061B (en) * | 2020-02-12 | 2021-08-03 | 南方电网科学研究院有限责任公司 | Method and device for judging oscillation stability based on apparent impedance |
CN111222297A (en) * | 2020-02-27 | 2020-06-02 | 华北电力大学 | Method for analyzing influence factors of subsynchronous oscillation of doubly-fed wind generator |
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 |
CN113725910A (en) * | 2021-10-08 | 2021-11-30 | 南通大学 | Stability analysis and quantitative evaluation method for wind power plant grid-connected system |
CN116822436A (en) * | 2023-06-30 | 2023-09-29 | 四川大学 | Oscillation risk sensitivity analysis method for direct-current transmission end alternating-current system |
CN116822436B (en) * | 2023-06-30 | 2024-02-27 | 四川大学 | Oscillation risk sensitivity analysis method for direct-current transmission end alternating-current system |
Also Published As
Publication number | Publication date |
---|---|
CN106021682B (en) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021682B (en) | Sub-synchronous oscillation analysis method and device based on impedance network model | |
CN106372339B (en) | The multi-rate simulating method and device of power electronics electric system | |
Ambriz-Perez et al. | Advanced SVC models for Newton-Raphson load flow and Newton optimal power flow studies | |
CN111179121B (en) | Power grid emergency control method based on expert system and deep reverse reinforcement learning | |
Kataoka et al. | Voltage stability limit of electric power systems with generator reactive power constraints considered | |
CN102222919B (en) | Power system reactive power optimization method based on improved differential evolution algorithm | |
CN107664954A (en) | A kind of electromagnetic transient in power system visual simulating model automatic forming method | |
CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
CN103049617B (en) | Method for simplifying passive electromagnetic transient-state simulation model of large-scale power distribution network | |
CN104104081B (en) | A kind of uncertain tidal current analysis method of non-iterative based on optimization method | |
CN103077268B (en) | Towards the state space method for automatic modeling of electromagnetic transient in power system emulation | |
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 | |
CN108828981A (en) | Rail traffic vehicles system Hardware In The Loop Simulation Method, apparatus and system | |
CN106208099A (en) | A kind of Method for Reactive Power Optimization in Power based on bi-level programming and application thereof | |
CN111181164B (en) | Improved master-slave split transmission and distribution cooperative power flow calculation method and system | |
Nigam et al. | A controller hardware-in-the-loop testbed: Verification and validation of microgrid control architectures | |
CN109657276A (en) | Electromagnetical transient emulation method and system | |
CN105281326A (en) | Slow coherency partition method considering damping influence of stator | |
WO2022052146A1 (en) | Heavy overload check method for load transfer decision of open-loop power grid | |
Tu et al. | Smart transformer modelling and hardware in-the-loop validation | |
Liao et al. | Machine learning methods for power flow control of multi-active-bridge converters | |
Li et al. | Construction of power grid digital twin model based on GAN | |
CN110233476A (en) | Voltage stability assessment method and relevant apparatus during a kind of black starting-up |
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 |