CN105281329B - A kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities - Google Patents
A kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities Download PDFInfo
- Publication number
- CN105281329B CN105281329B CN201510740422.7A CN201510740422A CN105281329B CN 105281329 B CN105281329 B CN 105281329B CN 201510740422 A CN201510740422 A CN 201510740422A CN 105281329 B CN105281329 B CN 105281329B
- Authority
- CN
- China
- Prior art keywords
- sequence
- generator
- battle array
- impedance
- formula
- 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
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to a kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities, belong to power system stability and control technical field.The admittance battle array of each sequence is inverted and obtains the Impedance Matrix of each sequence by the present invention;Port Impedance battle array is tried to achieve using dual-port theory first, and then tries to achieve comprehensive impedance battle array, it is inverted and obtains complex admittance battle array;According to complex admittance battle array modification positive sequence network, so that the positive sequence admittance battle array that is expanded;The equivalent transadmittance between 2 generator nodes is calculated using Ward equivalent methods, thus the generator frequency variation track under two kinds of reclosing time sequences is calculated, realize the quantitative evaluation to frequency variation track to obtain transient frequency stability margin, choose coincidence scheme of the corresponding time sequence of coincidence of transient frequency stability margin the greater for optimization.The present invention makes the analysis of complex fault calculate greatly simplified;Transient frequency stability when 2 machine systems carry out single-phase lock can be effectively improved, is worked well.
Description
Technical field
The present invention relates to a kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities, category
In power system stability and control technical field.
Background technology
Frequency is an important operating index in power system, is also the pith of power system stability operation.It
Reflecting active power that generator sent and load, (including burden with power, the loss in network that power plant is consumed have
The burden with power that workload and various electrical equipments are consumed) between dynamic equilibrium relation.In order to meet user's request,
The general fluctuation range no more than ± 0.2 Η z of the skew of frequency in power system, when system is not to be strict with, system
Frequency can not exceed ± 0.5Hz fluctuation range, so the fluctuation of frequency is to surpass in power system normal work
Cross allowed scope.If the fluctuation of frequency is beyond the scope allowed it is necessary to making frequency by corresponding control measure
Maintain in allowed limits.
After disturbing, the dynamic process of frequency of each department not only has the distribution of time, and has a distribution in space,
Just says different at the time of, the dynamic process of frequency difference of same node or areal;In synchronization, different nodes or
The dynamic process of the different regional frequency of person is also different.The amplitude of different regions hunting of frequency is also different.
Transient frequency is generally divided into system-wide transient frequency, regional transient frequency, generator transient frequency and general mother
The transient frequency of line.System-wide transient frequency refers to be acted on the total equivalent rotary inertia of system by system imbalance power summation
Produced incremental speed changes with time.Regional transient frequency refer to this area's imbalance power summation act on this area be used to
Property center total equivalent rotary inertia on produced by incremental speed change with time.Generator node transient frequency is to generate electricity
Machine incremental speed changes with time, namely with generator speed directly corresponding frequency.General bus transient frequency be secondly
The incremental speed that the space vector of voltage of phase voltage instantaneous value synthesis rotates in phase plane changes with time.
Single-pole reclosing sequential is adjusted, it is therefore an objective to the Secondary Shocks again reduced to system of being tried one's best when coinciding with permanent fault
With uneven energy, weaken the adverse effect to power network, improve the transmittability of network.It is single-phase under the stable visual angle of transient frequency
Adjusting for reclosing time sequence needs an index, and this index is exactly transient frequency stability margin.When system transient modelling frequency is steady
When determining nargin more than 0, system is stable, and numerical value is more big more is conducive to system stable;When system transient modelling frequency stabilization nargin
During less than 0, system is unstable, and the smaller stabilization for being more unfavorable for system, passes through the size of transient frequency stability margin value
To carry out the best match of single-pole reclosing sequential.
Dual-port is theoretical, Ward equivalent methods are a kind of simplifying methods of network, and it is by the electric power networks of complicated, high-order
Letter is simple low order network, carries out such simplification and not only can guarantee that computational accuracy, while computational efficiency can be greatly improved.
Patent of the present invention calculates generator frequency versus time curve under two kinds of time sequence of coincidence, and computing system transient frequency is steady
Determine nargin and realize the quantitative evaluation that track is changed over time to frequency, selecting system transient frequency margin of safety the greater is corresponding
Single-phase time sequence of coincidence is used as optimal coincidence setting program.
The content of the invention
The invention provides a kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities,
Placed an order the optimization problem of timing scheme of coinciding for solving the stable visual angle of transient frequency in 2 electromechanical Force systems.
The technical scheme is that:A kind of single-phase time sequence of coincidence for improving 2 electro-mechanical force system transient modelling frequency stabilities is whole
Determine method, when singlephase earth fault occurs for the transmission line of alternation current of 2 electromechanical Force systems, there are two kinds of failures in single-phase coincidence
I.e. singlephase earth fault and single-phase wire break failure, form the fault network of dual-port;Symmetrical component method is utilized according to fault type
Three sequence sequence nets are obtained, the admittance battle array of each sequence is write according to positive sequence, negative phase-sequence, zero-sequence network row, the admittance battle array of each sequence is inverted and obtains each sequence
Impedance Matrix;Port Impedance battle array is tried to achieve using dual-port theory first, and then tries to achieve comprehensive impedance battle array, it is inverted and integrated
Admittance battle array;According to complex admittance battle array modification positive sequence network, so that the positive sequence admittance battle array that is expanded;2 are calculated using Ward equivalent methods
Equivalent transadmittance between platform generator node, thus calculates the generator frequency variation track under two kinds of reclosing time sequences,
Realize the quantitative evaluation to frequency variation track to obtain transient frequency stability margin, choose transient frequency stability margin larger
The corresponding time sequence of coincidence of person is the coincidence scheme of optimization.
Methods described is comprised the following steps that:
(1) analog line head end input single-pole reclosing
1) positive sequence, negative phase-sequence, zero sequence admittance battle array are formed
Analog line head end puts into single-pole reclosing, equivalent to the multiple event that there occurs single-phase earthing and single-phase end broken string
Barrier, is initially formed positive sequence, negative phase-sequence, zero sequence admittance battle array under the multiple failure, is designated as Y1、Y2、Y0;Wherein, the method for admittance battle array is formed such as
Under:The number of nonzero element is equal to the connected earth-free circuitry number of corresponding node in each row off-diagonal element of admittance battle array;Admittance battle array
Each diagonal element, i.e., the self-admittance Y of each nodeiiEqual to the admittance sum on respective nodes institute's chord road:Admittance
The transadmittance Yij of each off-diagonal element of battle array is equal to connect between 2 nodes the negative value of admittance:Yij=-yij;
In formula, yijFor the admittance on institute's chord road between node i and node j;
2) positive sequence, negative phase-sequence, zero sequence impedance battle array are asked for
To positive sequence, negative phase-sequence, zero sequence admittance battle array Y1、Y2、Y0Invert, as shown in formula (1), obtain positive sequence, negative phase-sequence, zero sequence impedance
Battle array, is designated as Z1、Z2、Z0;
3) port Impedance is asked for using dual-port theory
1. positive sequence, negative phase-sequence, zero sequence self-impedance and the mutual impedance of port 1 and port 2 are asked for;
If corresponding node serial number is m, n at generation singlephase earth fault, n is 0 node of the earth;Single-phase end disconnection fault
Locate corresponding node serial number for p, q;It is that the corresponding port of m, n is denoted as port 1 by node serial number;Node serial number is that p, q are corresponding
Port is denoted as port 2;
By step 2) in the three sequence Impedance Matrixes that calculate corresponding element substitute into formula (2), try to achieve port 1 and port 2
Positive sequence, negative phase-sequence, the self-impedance of zero sequence and mutual impedance;
Z11(1)=Z11(2)=Zmm(1)+Znn(1)-2Zmn(1)
Z12(1)=Z12(2)=Zmp(1)+Znq(1)-Zmq(1)-Znp(1)
Z21(1)=Z21(2)=Zpm(1)+Zqn(1)-Zpn(1)-Zqm(1)
Z22(1)=Z22(2)=Zpp(1)+Zqq(1)-2Zpq(1)
Z11(0)=Zmm(0)+Znn(0)-2Zmn(0) (2)
Z12(0)=Zmp(0)+Znq(0)-Zmq(0)-Znp(0)
Z21(0)=Zpm(0)+Zqn(0)-Zpn(0)-Zqm(0)
Z22(0)=Z22(0)=Zpp(0)+Zqq(0)-2Zpq(0)
In formula, Z11(1)、Z22(1)The positive sequence self-impedance of port 1 and port 2 is represented respectively;Z12(1)、Z21(1)Port is represented respectively
Positive sequence mutual impedance between 1 and port 2;Z11(2)、Z22(2)The negative phase-sequence self-impedance of port 1 and port 2 is represented respectively;Z12(2)、Z21(2)
The negative phase-sequence mutual impedance between port 1 and port 2 is represented respectively;Z11(0)、Z22(0)Represent that the zero sequence of port 1 and port 2 hinders certainly respectively
It is anti-;Z12(0)、Z21(0)The zero sequence mutual impedance between port 1 and port 2 is represented respectively;Zuv(1)(u=m, n, p, q;V=m, n, p, q)
Represent Impedance Matrix Z1In u row v column elements, Zuv(2)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z2In u rows v arrange
Element, Zuv(0)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z0In u row v column elements;
2. the port Impedance battle array of the positive sequence, negative phase-sequence, zero sequence of port 1 and port 2 is sought
Obtained each sequence impedance is calculated by formula (2), positive sequence, negative phase-sequence, zero sequence port Impedance battle array is formed, Z is denoted as respectively(1)、
Z(2)、Z(0), as shown in formula (3);
4) comprehensive impedance battle array Z is asked forF
1. obtained each sequence impedance element is calculated according to by formula (2), row write negative phase-sequence, zero sequence port comprehensive impedance battle array Z, such as formula
(4) shown in;
2. negative phase-sequence shown in formula (4), zero sequence port comprehensive impedance battle array Z are substituted into formula (5) counter circuit Impedance Matrix ZL;
In formula,CTFor C transposition;Z′st(s=a, b, c;L=a, b, c), during s=l, represent back
Road s self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;
Formula (5) is 3. eliminated to the impedance element of the corresponding impedance elements of closed-loop path c, retention loop a and loop b by formula (6)
Element, obtains comprehensive impedance battle array ZF, as shown in formula (7);
In formula, Zsl(s=a, b;L=a, b) during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the mutual impedance behind closed-loop path between loop s and loop l;Z′sl(s=a, b;L=a, b) during s=l, represents loop s
Self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;Z′sc(s=a, b) between expression loop s and loop c
Mutual impedance;Z′cl(l=a b) represents the mutual impedance between loop c and loop l;Z′ccRepresent loop c self-impedance;
In formula, Zsl(s=a, b;L=a, b), during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the mutual impedance behind closed-loop path between loop s and loop l;
5) formula (8) is inverted and obtains complex admittance battle array, be denoted as YF, as shown in formula (8)
In formula, ysl(s=a, b;L=a, b), during s=l, represents to eliminate the self-admittance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the transadmittance behind closed-loop path between loop s and loop l;
6) complex admittance battle array Y is utilizedFIn element modification positive sequence admittance battle array Y1Middle corresponding element, amending method is shown in formula
(9) positive sequence that, is expanded admittance battle array Y1E
In formula, y 'uv(1) (u=m, n;V=m, n) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u
=m, n;V=m, n) represents admittance battle array Y1In u row v column elements;
y′uv(1)(u=p, q;V=p, q) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u=p, q;v
=p, q) represent admittance battle array Y1In u row v column elements;
7) according to extension positive sequence admittance battle array Y1E, the equivalent mutual conductance between 2 generator nodes is calculated using ward equivalent methods
Receive, be designated as Y12;Specific method is as follows:
1. the set constituted generator 1 and the corresponding node of generator 2 is denoted as remaining node group in boundary node B, network
Into set be denoted as external node E, will extension positive sequence admittance battle array Y1EMatrix in block form form is write, matrix in block form can be obtained and represented
Network equation such as formula (10) shown in;
In formula, YBBRepresent the admittance battle array corresponding to all boundary nodes;YBERepresent all boundary nodes and all outside segments
The corresponding admittance battle array of point, and YBEEqual to YEBTransposition;YEERepresent the admittance battle array corresponding to all external nodes;Represent institute
There is the column voltage vector of boundary point;Represent the column voltage vector of all external nodes;Represent the injection of all boundary nodes
Electric current;Represent the Injection Current of all external nodes;
2. the matrix in block form in formula (10) is substituted into (11) formula and calculated, obtain boundary admittance battle array
In formula,It is the boundary admittance battle array after equivalence, it includes the tribute of the equivalent branch road produced after external network abbreviation
Offer;Therefore after the equivalent transadmittance between the node of generator 1 that retains in 2 machine systems and the node of generator 2 is exactly equivalence
Boundary admittance, so there is formula (12) establishment:
In formula, G12For equivalent transadmittance Y12Real part, referred to as conductance, and G12> 0;B12For equivalent transadmittance Y12Void
Portion, referred to as susceptance, and B12> 0;
8) by step 7) the middle equivalent transadmittance Y calculated between obtained generator 1 and generator 212In substitution formula (13),
Local derviation K of the active power to respective generator rotor angle of generator 1 and generator 2 is calculated respectively11、K22;
In formula, PG1、PG2The respectively active power of generator 1 and generator 2;δ1、δ2Respectively generator 1 and generator
2 generator rotor angle; Represent that the active power of generator 1 and generator 2 exists to the partial derivative of generator's power and angle respectively
The numerical value of initial time;Eq1、Eq2The respectively electromotive force of generator 1 and generator 2;δ0Phase at expression system initial work location
Angle;|Y12| represent equivalent transadmittance Y12Amplitude,γ represents equivalent transadmittance Y12Phase angle, γ=
arctg(B12/G12), and γ ∈ (0, pi/2);
9) by step 8) in calculate obtained K11、K222 machine systems is temporary when the single-phase coincidence of head end is calculated in substitution formula (14)
The transient state angular frequency of state angular frequency, generator 1 and generator 2 is denoted as Δ ω respectively1(t)、Δω2(t);
In formula, D1、D2The damped coefficient of generator 1 and generator 2 is represented respectively;ΔPm10、ΔPm20Generator is represented respectively
1 and generator 2 occur chugging;M1、M2Represent the rotary inertia of generator 1,2; ω0For
The rated angular velocity of generator, t is the time;
10) because the relation between angular frequency and frequency meets formula (15), by step 9) in calculate obtained generator 1,2
Angular frequency Δ ω1(t)、Δω2(t) formula (16) is substituted into, the transient frequency Δ f of generator 1,2 when calculating head end single-phase coincidence1
(t)、Δf2(t);
The π (15) of f=ω/2
In formula, f is frequency;ω is angular frequency;
(2) analog line end input single-pole reclosing
1) positive sequence, negative phase-sequence, zero sequence admittance battle array are formed
Analog line end puts into single-pole reclosing, equivalent to the multiple event that there occurs single-phase earthing and single-phase head end broken string
Barrier, is initially formed positive sequence, negative phase-sequence, zero sequence admittance battle array under the multiple failure, is designated as Y '1、Y′2、Y′0;
2) by Y '1、Y′2、Y′0Substitution formula (1), positive sequence when obtaining the single-phase coincidence of line end of inverting, negative phase-sequence, zero sequence resistance
Anti- battle array, is designated as Z '1、Z′2、Z′0;
3) 3 in repeat step (1)) port Impedance is asked for using dual-port theory, it is designated as Z '(1)、Z′(2)、Z′(0);
4) 4 in repeat step (1)) comprehensive impedance battle array is asked for, it is designated as Z 'F;
5) 5 in repeat step (1)) complex admittance battle array is obtained, it is denoted as Y 'F;
6) 6 in repeat step (1)) utilize complex admittance battle array Y 'FMiddle element modification positive sequence admittance battle array Y '1Middle corresponding element,
The positive sequence that is expanded admittance battle array Y '1E;
7) 7 in repeat step (1)) according to extension positive sequence admittance battle array Y '1E, 2 generators are calculated using ward equivalent methods
Equivalent transadmittance between node, is designated as Y '12;
8) 8 in repeat step (1)) local derviation of the active power to respective generator rotor angle of generator 1 and generator 2 is calculated respectively
K′11、K′22;
9) 9 in repeat step (1)) the transient state angular frequency difference of generator 1 and generator 2 when calculating the single-phase coincidence in end
It is denoted as Δ ω '1(t)、Δω′2(t);
10) 10 in repeat step (1)) the transient frequency Δ f ' of generator 1,2 when calculating the single-phase coincidence in end1
(t)、Δf′2(t);
(3) frequency shift (FS) security quantitative evaluation index
Obtained generator transient frequency is calculated to step (1), (2) using formula (17) and carries out quantitative evaluation, circuit is obtained
The transient frequency stability margin η of generator 1,2 during head end single-phase coincidenceG1、ηG2, generator 1,2 during line end single-phase coincidence
Transient frequency stability margin η 'G1、η′G2;
η=[fext-(fcr-kTcr)] × 100% (17)
In formula, fcrAnd TcrIt is the frequency shift (FS) threshold value of generator and the duration allowed, f respectivelyextRefer to transient state mistake
The limiting value of generator frequency in journey, k is critical frequency shift duration to be converted into the commutation factor of frequency, η for just or
Negative value represents that transient frequency is stable or unstable respectively;
The transient frequency stability margin of each generator in 2 machine systems is obtained, each generator transient frequency stability margin is taken most
The small transient frequency stability margin being worth as system;
System transient modelling frequency stabilization nargin during the single-phase coincidence of head end is calculated by formula (18), is denoted as η1;
η1=min { ηG1,ηG2} (18)
System transient modelling frequency stabilization nargin during the single-phase coincidence in end is calculated by formula (19), is denoted as η2;
η2=min { η 'G1,η′G2} (19)
(4) draw a conclusion
System transient modelling frequency stabilization nargin η during the single-phase coincidence of alternative route first and end1、η2, choose the greater corresponding
Single-phase time sequence of coincidence as output optimization time sequence of coincidence scheme;If η1> η2, then single-phase coincidence is put into by circuit head end first;
Conversely, putting into single-phase coincidence first by line end.
The present invention operation principle be:
1. dual-port is theoretical
(1) port Impedance is asked for using dual-port theory
When at F1 and at F2 unbalanced fault occurs for power system simultaneously, accessed simultaneously equivalent at F1 and at F2
The impedance of asymmetrical three-phase, is made up of at failure two parts circuit, and a part is the unsymmetric circuit of faults situation, separately
A part is the reflection symmetrical original system of triphase parameter.Due to there is failure at two, one is shown as at each failure in sequence net
Individual port, therefore the sequence net of double failure is two-port network, i.e. two-port network.
1. positive sequence, negative phase-sequence, zero sequence self-impedance and the mutual impedance of port 1 and port 2 are asked for.
If corresponding node serial number is m, n at generation singlephase earth fault (n is 0 node of the earth);Line end is single-phase disconnected
Corresponding node serial number is p, q at line.It is that the corresponding port of m, n is denoted as port 1 by node serial number;Node serial number is p, q correspondence
Port be denoted as port 2.
Corresponding element in positive sequence, negative phase-sequence, zero sequence impedance battle array is substituted into formula (20), try to achieve port 1 and port 2 positive sequence,
Negative phase-sequence, the self-impedance of zero sequence and mutual impedance.
Z11(1)=Z11(2)=Zmm(1)+Znn(1)-2Zmn(1)
Z12(1)=Z12(2)=Zmp(1)+Znq(1)-Zmq(1)-Znp(1)
Z21(1)=Z21(2)=Zpm(1)+Zqn(1)-Zpn(1)-Zqm(1)
Z22(1)=Z22(2)=Zpp(1)+Zqq(1)-2Zpq(1)
Z11(0)=Zmm(0)+Znn(0)-2Zmn(0) (20)
Z12(0)=Zmp(0)+Znq(0)-Zmq(0)-Znp(0)
Z21(0)=Zpm(0)+Zqn(0)-Zpn(0)-Zqm(0)
Z22(0)=Z22(0)=Zpp(0)+Zqq(0)-2Zpq(0)
In formula, Z11(1)、Z22(1)The positive sequence self-impedance of port 1 and port 2 is represented respectively;Z12(1)、Z21(1)Port is represented respectively
Positive sequence mutual impedance between 1 and port 2;Z11(2)、Z22(2)The negative phase-sequence self-impedance of port 1 and port 2 is represented respectively;Z12(2)、Z21(2)
The negative phase-sequence mutual impedance between port 1 and port 2 is represented respectively;Z11(0)、Z22(0)Represent that the zero sequence of port 1 and port 2 hinders certainly respectively
It is anti-;Z12(0)、Z21(0)The zero sequence mutual impedance between port 1 and port 2 is represented respectively;Zuv(1)(u=m, n, p, q;V=m, n, p, q)
Represent Impedance Matrix Z1In u row v column elements, Zuv(2)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z2In u rows v arrange
Element, Zuv(0)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z0In u row v column elements;
2. the port Impedance battle array of the positive sequence, negative phase-sequence, zero sequence of port 1 and port 2 is sought.
Obtained each sequence impedance is calculated by formula (20), positive sequence, negative phase-sequence, zero sequence port Impedance battle array is formed, Z is denoted as respectively(1)、
Z(2)、Z(0), as shown in formula (21).
(2) comprehensive impedance battle array is asked for
1. negative phase-sequence, zero sequence port comprehensive impedance battle array Z are write according to two-port network row, as shown in formula (22).
2. negative phase-sequence shown in formula (22), zero sequence port comprehensive impedance battle array Z are substituted into formula (23) counter circuit Impedance Matrix ZL。
In formula,CTFor C transposition;Z′st(s=a, b, c;L=a, b, c), during s=l, represent back
Road s self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;
3. formula (23) is eliminated into the corresponding impedance elements of closed-loop path c, retention loop a and loop b impedance by formula (24)
Element, obtains comprehensive impedance battle array ZF, as shown in formula (25).
In formula, Zsl(s=a, b;L=a, b) during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the mutual impedance behind closed-loop path between loop s and loop l;Z′sl(s=a, b;L=a, b) during s=l, represents loop s
Self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;Z′sc(s=a, b) between expression loop s and loop c
Mutual impedance;Z′cl(l=a b) represents the mutual impedance between loop c and loop l;Z′ccRepresent loop c self-impedance;
In formula, Zsl(s=a, b;L=a, b), during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the mutual impedance behind closed-loop path between loop s and loop l;
(3) formula (25) is inverted and obtains complex admittance battle array, be denoted as YF, as shown in formula (26).
In formula, ysl(s=a, b;L=a, b), during s=l, represents to eliminate the self-admittance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the transadmittance behind closed-loop path between loop s and loop l;
(4) complex admittance battle array Y is utilizedFIn element modification positive sequence admittance battle array Y1Middle corresponding element, amending method is shown in formula
(27) the admittance battle array that, is expanded Y1E。
In formula, y 'uv(1)(u=m, n;V=m, n) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u=
m,n;V=m, n) represents admittance battle array Y1In u row v column elements;
y′uv(1)(u=p, q;V=p, q) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u=p, q;v
=p, q) represent admittance battle array Y1In u row v column elements.
2.Ward equivalent methods
The node of former network is represented with set N.The part referred to as external network that abbreviation falls is intended to, its set of node is represented with E.Protect
The node of subnetwork is stayed to be represented with collection G is retained.Then there are G ∩ E=N.The associated section of external network node is neutralized retaining collection
Point composition boundary node set, is represented with B.The part of outside segments of getting along well point set association is internal node collection, is represented with I.If will lead
Receive short matrix representation network equation by I, B, E gather divide, then can write out the network equation represented with matrix in block form form such as
Under:
Eliminate the voltage quantities of external nodeHave
In formula,It is the boundary admittance matrix after equivalence, it includes the equivalent branch road that is produced after external network abbreviation
Contribution.
The order of equation time of formula (29) is low compared with the order of formula (28), therefore easily calculated.Due to equivalent boundary admittance square
Battle arrayAn addition Item is added relative to original YBB, and addition Item is full battle array under many circumstances, soIt is non-
Null element is more than YBB, i.e.,By the sparse part of some lost.WhenDimension is relatively low, i.e., solved when boundary node is less after simplifying
Network equation (29) is faster than solving former network equation (28).
3. single-phase time sequence of coincidence is to the Influencing Mechanism of system transient modelling frequency stability
The formation of 3.1 admittance battle array
2 machine systems are as shown in Figure 1 under static load model.2 machine systems are powered through double back transmission line of alternation current to load.When
When singlephase earth fault occurs for non-head end any point on circuit, when putting into single-pole reclosing by circuit head end and end respectively,
Calculate the change of bus frequency when first closing head end and first closing end.
The changing rule of 2 machine system each internal potential frequency of generator when original machine power is undergone mutation, is not examined in analysis
Consider the influence of regulating system and part throttle characteristics, focus on empty to disturbance place and frequency in the research of distribution.For ease of right
The comprehensive analysis of various phenomenons, carries out relevant theory deduction first.
If 1#With 2#Equivalent transadmittance between machine interior nodes is Y12=G12+jB12, the discussion of 2 machine system equivalent transadmittances,
In Fig. 1, if line impedance is ZL=rL+jxL, load impedance represents with constant impedance, i.e. ZDi=RDi+jXDiOr YDi=GDi-
jBDi, i=1,2.Wherein, ZLFor line impedance, rLFor line resistance, xLFor line reactance;ZDiFor load constant impedance, RDiFor
Load resistance, XDiFor load reactance;YDiFor load constant admittance, GDiFor load conductance, BDiFor susceptance.
Set againAnd(wherein, xdiFor the reactance of generator, yiFor leading for two generator nodes
Receive, yLFor the admittance of circuit.) can obtain 2 machine systems admittance battle array be:
Load bus is eliminated to obtain:In formula
Order
It can obtain:
R12=-rL(xd2BD2+xd1BD1+1)+rL(GD1GD2-BD1BD2)xd1xd2
+xL(xd1GD1+xd2GD2)+xd1xd2(GD1+GD2)+xd1xd2xL(GD1BD2+GD2BD1)
X12=-(xd1+xd2+xL)-rL(xd1xd2+xd2GD2+xd1xd2(GD1BD2+GD2BD1))
-xd1xd2(BD1+BD2)-xL(xd1BD1+xd2BD2)-xd1xd2xL(BD1BD2-GD1GD2)
General rLVery little, so there is R12> 0, X12< 0, therefore G12> 0, B12> 0.Y12=G12+jB12As madeγ=arctg (B12/G12), γ ∈ (0, pi/2).
The transient frequency of 3.2 systems
When system is uniform damping, i.e. D1/M1=D2/M2During=σ, the analytic solutions of frequency procedure can be obtained:
In formula,Y12=G12+jB12, wherein,
G12It is the real part of admittance, and G12> 0;B12It is the imaginary part of admittance, and B12> 0.|Y12| represent equivalent transadmittance Y12Amplitude, andγ represents admittance Y12Phase angle, and γ=arctg (B12/G12), then there are γ ∈ (0, pi/2);Generator
The 1 and active power local derviation K respectively to respective generator rotor angle of generator 211、K22;PG1、PG2Generator 1 and generator 2 are represented respectively
Have power;δ1、δ2The generator rotor angle of generator 1 and generator 2 is represented respectively;δ0Phase angle at expression system initial work location;Di、
(i=1,2) represent the damped coefficient of generator;ΔPmi0(i=1,2) represent the chugging that prime mover occurs;M1、M2Represent the rotary inertia of generator;ω0For generator
Rated angular velocity.
4. transient frequency stability margin
One group of two-element list [(f shown in Fig. 2cr1,Tcr1),(fcr2,Tcr2),...,(fcrl,Tcrl)] it is used to description one
The TFDA problems of bus.If for the threshold value f of each frequencycrlLess than power network rated frequency (fr) two-element list, bus
Actual frequency is less than threshold value fcrlDuration be both less than corresponding threshold value Tcrl;Simultaneously for each threshold value fcrl
Higher than power network rated frequency (fr) two-element list, bus actual frequency exceed threshold value fcrlDuration it is both less than corresponding
Threshold value Tcrl, then it is assumed that the transient frequency of the bus is stable.Otherwise, it is believed that the transient frequency of the bus is unstable.Its
In, fcr1,...,fcrlFor one group of given frequency threshold value;Tcr1,...,TcrlFor the given frequency threshold value correspondence duration
Threshold value.
The stability margin of transient frequency is calculated according to transient frequency versus time curve, it is typically that transient frequency is stable
Nargin is denoted as ηfb, formula (32) is its definition:
ηfb=[fext-(fcr-kTcr)] × 100% (32)
Wherein fcrAnd TcrIt is the frequency shift (FS) threshold value of generator and the duration allowed, f respectivelyextRefer to transient process
The limiting value of middle generator frequency, k is the commutation factor for critical frequency shift duration being converted into frequency, ηfbFor just (or
It is negative) value represents that transient frequency is stable (or unstable).
The beneficial effects of the invention are as follows:
1st, using dual-port theory solve it is single-phase coincide with permanent fault when port battle array, calculate the analysis of complex fault
It is greatly simplified.
2nd, in the equivalent transadmittance between asking for 2 generator nodes, using Ward equivalent methods to extension positive sequence admittance battle array
Abbreviation is carried out, network size is not only reduced, and improve computational efficiency.
3rd, the stabilization of general frequency is to be compared by the steady-state value of frequency response with given threshold value, so as to judge
Frequency is stable or unstable.The present invention obtains characterizing system transient modelling frequency stability using the quantitative evaluation of transient frequency
Quantizating index, so as to realize the preferred of single-phase time sequence of coincidence scheme.
4th, a large amount of simulation results show, patent of the present invention can effectively improve transient frequency when 2 machine systems carry out single-phase lock
Stability, works well.
Brief description of the drawings
Fig. 1 is 2 machine system wiring figures;Including 2 generator G1、G2With 2 transformer T1、T2By transmission line of alternation current
Load is powered;E in figureq1、Eq2Respectively G1、G2Electromotive force;PDA、QDARespectively bus A load active power and idle
Power;PDC、QDCRespectively bus C load active power and reactive power;PDD、QDDRespectively bus D load active power
And reactive power;
Fig. 2 is description TFDA two-element list;In figure, f is frequency;frFor power network rated frequency;fcr1、fcr2、fcr3Respectively
Different frequency threshold values;T1、T2、T3Respectively frequency threshold value fcr1、fcr2、fcr3Corresponding time gate threshold value;When t/s is
Between/second;
Fig. 3 is bus 1-bus D circuit head end single phase ground faults, generator G during head end first single-phase coincidence1's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Fig. 4 is bus 1-bus D circuit head end single phase ground faults, generator G during head end first single-phase coincidence2's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Fig. 5 is bus 1-bus D circuit head end single phase ground faults, generator G during the first single-phase coincidence of end1's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Fig. 6 is bus 1-bus D circuit head end single phase ground faults, generator G during the first single-phase coincidence of end2's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Fig. 7 is bus 2-bus D circuit head end single phase ground faults, generator G during head end first single-phase coincidence1's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Fig. 8 is bus 2-bus D circuit head end single phase ground faults, generator G during head end first single-phase coincidence2's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Fig. 9 is bus 2-bus D circuit head end single phase ground faults, generator G during the first single-phase coincidence of end1's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second);
Figure 10 is bus 2-bus D circuit head end single phase ground faults, generator G during the first single-phase coincidence of end2's
Frequency variation curve;In figure, generator frequency curve (Hz) is generator frequency (hertz);Time (s) is the time (second).
Embodiment
Embodiment 1:As Figure 1-10 shows, a kind of single-phase time sequence of coincidence for improving 2 electro-mechanical force system transient modelling frequency stabilities
, when singlephase earth fault occurs for the transmission line of alternation current of 2 electromechanical Force systems, there are two kinds of events in single-phase coincidence in setting method
Barrier is singlephase earth fault and single-phase wire break failure, forms the fault network of dual-port;Symmetrical components are utilized according to fault type
Method obtains three sequence sequence nets, and the admittance battle array of each sequence is write according to positive sequence, negative phase-sequence, zero-sequence network row, and the admittance battle array of each sequence is inverted and obtains each
The Impedance Matrix of sequence;Port Impedance battle array is tried to achieve using dual-port theory first, and then tries to achieve comprehensive impedance battle array, it is inverted and obtains comprehensive
Close admittance battle array;According to complex admittance battle array modification positive sequence network, so that the positive sequence admittance battle array that is expanded;Calculated using Ward equivalent methods
Equivalent transadmittance between 2 generator nodes, thus calculates the generator frequency variation track under two kinds of reclosing time sequences,
Realize the quantitative evaluation to frequency variation track to obtain transient frequency stability margin, choose transient frequency stability margin larger
The corresponding time sequence of coincidence of person is the coincidence scheme of optimization.
Methods described is comprised the following steps that:
(1) analog line head end input single-pole reclosing
1) positive sequence, negative phase-sequence, zero sequence admittance battle array are formed
Analog line head end puts into single-pole reclosing, equivalent to the multiple event that there occurs single-phase earthing and single-phase end broken string
Barrier, is initially formed positive sequence, negative phase-sequence, zero sequence admittance battle array under the multiple failure, is designated as Y1、Y2、Y0;Wherein, the method for admittance battle array is formed such as
Under:The number of nonzero element is equal to the connected earth-free circuitry number of corresponding node in each row off-diagonal element of admittance battle array;Admittance battle array
Each diagonal element, i.e., the self-admittance Y of each nodeiiEqual to the admittance sum on respective nodes institute's chord road:Admittance
The transadmittance Y of each off-diagonal element of battle arrayijIt is equal to connect between 2 nodes the negative value of admittance:Yij=-yij;
In formula, yijFor the admittance on institute's chord road between node i and node j;
2) positive sequence, negative phase-sequence, zero sequence impedance battle array are asked for
To positive sequence, negative phase-sequence, zero sequence admittance battle array Y1、Y2、Y0Invert, as shown in formula (1), obtain positive sequence, negative phase-sequence, zero sequence impedance
Battle array, is designated as Z1、Z2、Z0;
3) port Impedance is asked for using dual-port theory
1. positive sequence, negative phase-sequence, zero sequence self-impedance and the mutual impedance of port 1 and port 2 are asked for;
If corresponding node serial number is m, n at generation singlephase earth fault, n is 0 node of the earth;Single-phase end disconnection fault
Locate corresponding node serial number for p, q;It is that the corresponding port of m, n is denoted as port 1 by node serial number;Node serial number is that p, q are corresponding
Port is denoted as port 2;
By step 2) in the three sequence Impedance Matrixes that calculate corresponding element substitute into formula (2), try to achieve port 1 and port 2
Positive sequence, negative phase-sequence, the self-impedance of zero sequence and mutual impedance;
Z11(1)=Z11(2)=Zmm(1)+Znn(1)-2Zmn(1)
Z12(1)=Z12(2)=Zmp(1)+Znq(1)-Zmq(1)-Znp(1)
Z21(1)=Z21(2)=Zpm(1)+Zqn(1)-Zpn(1)-Zqm(1)
Z22(1)=Z22(2)=Zpp(1)+Zqq(1)-2Zpq(1)
Z11(0)=Zmm(0)+Znn(0)-2Zmn(0) (2)
Z12(0)=Zmp(0)+Znq(0)-Zmq(0)-Znp(0)
Z21(0)=Zpm(0)+Zqn(0)-Zpn(0)-Zqm(0)
Z22(0)=Z22(0)=Zpp(0)+Zqq(0)-2Zpq(0)
In formula, Z11(1)、Z22(1)The positive sequence self-impedance of port 1 and port 2 is represented respectively;Z12(1)、Z21(1)Port is represented respectively
Positive sequence mutual impedance between 1 and port 2;Z11(2)、Z22(2)The negative phase-sequence self-impedance of port 1 and port 2 is represented respectively;Z12(2)、Z21(2)
The negative phase-sequence mutual impedance between port 1 and port 2 is represented respectively;Z11(0)、Z22(0)Represent that the zero sequence of port 1 and port 2 hinders certainly respectively
It is anti-;Z12(0)、Z21(0)The zero sequence mutual impedance between port 1 and port 2 is represented respectively;Zuv(1)(u=m, n, p, q;V=m, n, p, q)
Represent Impedance Matrix Z1In u row v column elements, Zuv(2)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z2In u rows v arrange
Element, Zuv(0)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z0In u row v column elements;
2. the port Impedance battle array of the positive sequence, negative phase-sequence, zero sequence of port 1 and port 2 is sought
Obtained each sequence impedance is calculated by formula (2), positive sequence, negative phase-sequence, zero sequence port Impedance battle array is formed, Z is denoted as respectively(1)、
Z(2)、Z(0), as shown in formula (3);
4) comprehensive impedance battle array Z is asked forF
1. obtained each sequence impedance element is calculated according to by formula (2), row write negative phase-sequence, zero sequence port comprehensive impedance battle array Z, such as formula
(4) shown in;
2. negative phase-sequence shown in formula (4), zero sequence port comprehensive impedance battle array Z are substituted into formula (5) counter circuit Impedance Matrix ZL;
In formula,CTFor C transposition;Z′st(s=a, b, c;L=a, b, c), during s=l, represent back
Road s self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;
Formula (5) is 3. eliminated to the impedance element of the corresponding impedance elements of closed-loop path c, retention loop a and loop b by formula (6)
Element, obtains comprehensive impedance battle array ZF, as shown in formula (7);
In formula, Zsl(s=a, b;L=a, b) during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the mutual impedance behind closed-loop path between loop s and loop l;Z′sl(s=a, b;L=a, b) during s=l, represents loop s
Self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;Z′sc(s=a, b) between expression loop s and loop c
Mutual impedance;Z′cl(l=a b) represents the mutual impedance between loop c and loop l;Z′ccRepresent loop c self-impedance;
In formula, Zsl(s=a, b;L=a, b), during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the mutual impedance behind closed-loop path between loop s and loop l;
5) formula (8) is inverted and obtains complex admittance battle array, be denoted as YF, as shown in formula (8)
In formula, ysl(s=a, b;L=a, b), during s=l, represents to eliminate the self-admittance of loop s behind closed-loop path;During s ≠ l,
Represent to eliminate the transadmittance behind closed-loop path between loop s and loop l;
6) complex admittance battle array Y is utilizedFIn element modification positive sequence admittance battle array Y1Middle corresponding element, amending method is shown in formula
(9) positive sequence that, is expanded admittance battle array Y1E
In formula, y 'uv(1)(u=m, n;V=m, n) represents extension positive sequence admittance battle array Y1EIn u row v column elements;
yuv(1)(u=m, n;V=m, n) represents admittance battle array Y1In u row v column elements;
y′uv(1)(u=p, q;V=p, q) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u=p, q;v
=p, q) represent admittance battle array Y1In u row v column elements;
7) according to extension positive sequence admittance battle array Y1E, the equivalent mutual conductance between 2 generator nodes is calculated using ward equivalent methods
Receive, be designated as Y12;Specific method is as follows:
1. the set constituted generator 1 and the corresponding node of generator 2 is denoted as remaining node group in boundary node B, network
Into set be denoted as external node E, will extension positive sequence admittance battle array Y1EMatrix in block form form is write, matrix in block form can be obtained and represented
Network equation such as formula (10) shown in;
In formula, YBBRepresent the admittance battle array corresponding to all boundary nodes;YBERepresent all boundary nodes and all outside segments
The corresponding admittance battle array of point, and YBEEqual to YEBTransposition;YEERepresent the admittance battle array corresponding to all external nodes;Represent institute
There is the column voltage vector of boundary point;Represent the column voltage vector of all external nodes;Represent the injection of all boundary nodes
Electric current;Represent the Injection Current of all external nodes;
2. the matrix in block form in formula (10) is substituted into (11) formula and calculated, obtain boundary admittance battle array
In formula,It is the boundary admittance battle array after equivalence, it includes the tribute of the equivalent branch road produced after external network abbreviation
Offer;Therefore after the equivalent transadmittance between the node of generator 1 that retains in 2 machine systems and the node of generator 2 is exactly equivalence
Boundary admittance, so there is formula (12) establishment:
In formula, G12For equivalent transadmittance Y12 real part, referred to as conductance, and G12> 0;B12For equivalent transadmittance Y12Void
Portion, referred to as susceptance, and B12> 0;
8) by step 7) the middle equivalent transadmittance Y calculated between obtained generator 1 and generator 212In substitution formula (13),
Local derviation K of the active power to respective generator rotor angle of generator 1 and generator 2 is calculated respectively11、K22;
In formula, PG1、PG2The respectively active power of generator 1 and generator 2;δ1、δ2Respectively generator 1 and generator
2 generator rotor angle; Partial derivative of the active power to generator's power and angle of generator 1 and generator 2 is represented respectively
In the numerical value of initial time;Eq1、Eq2The respectively electromotive force of generator 1 and generator 2;δ0At expression system initial work location
Phase angle;|Y12| represent equivalent transadmittance Y12Amplitude,γ represents equivalent transadmittance Y12Phase angle, γ
=arctg (B12/G12), and γ ∈ (0, pi/2);
9) by step 8) in calculate obtained K11、K222 machine systems is temporary when the single-phase coincidence of head end is calculated in substitution formula (14)
The transient state angular frequency of state angular frequency, generator 1 and generator 2 is denoted as Δ ω respectively1(t)、Δω2(t);
In formula, D1、D2The damped coefficient of generator 1 and generator 2 is represented respectively;ΔPm10、ΔPm20Generator is represented respectively
1 and generator 2 occur chugging;M1、M2Represent the rotary inertia of generator 1,2; ω0For
The rated angular velocity of generator, t is the time;
10) because the relation between angular frequency and frequency meets formula (15), by step 9) in calculate obtained generator 1,2
Angular frequency Δ ω1(t)、Δω2(t) formula (16) is substituted into, the transient frequency Δ f of generator 1,2 when calculating head end single-phase coincidence1
(t)、Δf2(t);
The π (15) of f=ω/2
In formula, f is frequency;ω is angular frequency;
(2) analog line end input single-pole reclosing
1) positive sequence, negative phase-sequence, zero sequence admittance battle array are formed
Analog line end puts into single-pole reclosing, equivalent to the multiple event that there occurs single-phase earthing and single-phase head end broken string
Barrier, is initially formed positive sequence, negative phase-sequence, zero sequence admittance battle array under the multiple failure, is designated as Y '1、Y′2、Y′0;
2) by Y '1、Y′2、Y′0Substitution formula (1), positive sequence when obtaining the single-phase coincidence of line end of inverting, negative phase-sequence, zero sequence resistance
Anti- battle array, is designated as Z '1、Z′2、Z′0;
3) 3 in repeat step (1)) port Impedance is asked for using dual-port theory, it is designated as Z '(1)、Z′(2)、Z′(0);
4) 4 in repeat step (1)) comprehensive impedance battle array is asked for, it is designated as Z 'F;
5) 5 in repeat step (1)) complex admittance battle array is obtained, it is denoted as Y 'F;
6) 6 in repeat step (1)) utilize complex admittance battle array Y 'FMiddle element modification positive sequence admittance battle array Y '1Middle corresponding element,
The positive sequence that is expanded admittance battle array Y '1E;
7) 7 in repeat step (1)) according to extension positive sequence admittance battle array Y '1E, 2 generators are calculated using ward equivalent methods
Equivalent transadmittance between node, is designated as Y '12;
8) 8 in repeat step (1)) local derviation of the active power to respective generator rotor angle of generator 1 and generator 2 is calculated respectively
K′11、K′22;
9) 9 in repeat step (1)) the transient state angular frequency difference of generator 1 and generator 2 when calculating the single-phase coincidence in end
It is denoted as Δ ω '1(t)、Δω′2(t);
10) 10 in repeat step (1)) the transient frequency Δ f ' of generator 1,2 when calculating the single-phase coincidence in end1
(t)、Δf′2(t);
(3) frequency shift (FS) security quantitative evaluation index
Obtained generator transient frequency is calculated to step (1), (2) using formula (17) and carries out quantitative evaluation, circuit is obtained
The transient frequency stability margin η of generator 1,2 during head end single-phase coincidenceG1、ηG2, generator 1,2 during line end single-phase coincidence
Transient frequency stability margin η 'G1、η′G2;
η=[fext-(fcr-kTcr)] × 100% (17)
In formula, fcrAnd TcrIt is the frequency shift (FS) threshold value of generator and the duration allowed, f respectivelyextRefer to transient state mistake
The limiting value of generator frequency in journey, k is critical frequency shift duration to be converted into the commutation factor of frequency, η for just or
Negative value represents that transient frequency is stable or unstable respectively;
The transient frequency stability margin of each generator in 2 machine systems is obtained, each generator transient frequency stability margin is taken most
The small transient frequency stability margin being worth as system;
System transient modelling frequency stabilization nargin during the single-phase coincidence of head end is calculated by formula (18), is denoted as η1;
η1=min { ηG1,ηG2} (18)
System transient modelling frequency stabilization nargin during the single-phase coincidence in end is calculated by formula (19), is denoted as η2;
η2=min { η 'G1,η′G2} (19)
(4) draw a conclusion
System transient modelling frequency stabilization nargin η during the single-phase coincidence of alternative route first and end1、η2, choose the greater corresponding
Single-phase time sequence of coincidence as output optimization time sequence of coincidence scheme;If η1> η2, then single-phase coincidence is put into by circuit head end first;
Conversely, putting into single-phase coincidence first by line end.
Embodiment 2:As Figure 1-10 shows, a kind of single-phase time sequence of coincidence for improving 2 electro-mechanical force system transient modelling frequency stabilities
, when singlephase earth fault occurs for the transmission line of alternation current of 2 electromechanical Force systems, there are two kinds of events in single-phase coincidence in setting method
Barrier is singlephase earth fault and single-phase wire break failure, forms the fault network of dual-port;Symmetrical components are utilized according to fault type
Method obtains three sequence sequence nets, and the admittance battle array of each sequence is write according to positive sequence, negative phase-sequence, zero-sequence network row, and the admittance battle array of each sequence is inverted and obtains each
The Impedance Matrix of sequence;Port Impedance battle array is tried to achieve using dual-port theory first, and then tries to achieve comprehensive impedance battle array, it is inverted and obtains comprehensive
Close admittance battle array;According to complex admittance battle array modification positive sequence network, so that the positive sequence admittance battle array that is expanded;Calculated using Ward equivalent methods
Equivalent transadmittance between 2 generator nodes, thus calculates the generator frequency variation track under two kinds of reclosing time sequences,
Realize the quantitative evaluation to frequency variation track to obtain transient frequency stability margin, choose transient frequency stability margin larger
The corresponding time sequence of coincidence of person is the coincidence scheme of optimization.
Embodiment 3:As Figure 1-10 shows,
By taking 2 machine systems shown in Fig. 1 as an example, 8 buses, the node parameter such as table of each generator and bus are had in system
Shown in 1, the voltage class of generator 1 is 18.0kV, and the voltage class of generator 1 is 13.8kV, and busbar voltage grade is 230kV.Line
The circuit per unit value parameter such as resistance, reactance and the admittance on road is as shown in table 2.8 line voltage distribution grade is 230kV.Transformer
The parameter such as resistance, reactance and no-load voltage ratio as shown in table 3, the inertia time constant of generator, the reactance of d-axis transient state, quadrature axis transient state electricity
The parameter such as anti-is as shown in table 4.The power reference value for the parameter per unit value being related in table 2- tables 4 takes 100MVA.
A phase singlephase earth faults occur for 2 machine system busbar 1-bus D circuit head ends 0s shown in Fig. 1, and 0.1s circuit A phases are first
End tripping operation, 0.2s circuit A phases end tripping operation, 1s circuit A phases head end is overlapped, the tripping operation of 1.2s faulty line first and ends three-phase.Formed
Positive sequence, negative phase-sequence, zero sequence admittance battle array when faulty line A phases head end when 1s is overlapped, end is broken, lead to positive sequence, negative phase-sequence, zero sequence
Battle array of receiving is inverted, and obtains positive sequence, negative phase-sequence, zero sequence impedance battle array, and positive sequence, negative phase-sequence, the resistance of zero sequence port are asked for using two-port network theory
Anti- battle array, is eliminated behind closed-loop path, obtains comprehensive impedance battle array, and comprehensive impedance battle array is inverted, and is obtained complex admittance battle array, is utilized WARD etc.
Value method, eliminates after the external node except generator 1 and generator 2, obtains the equivalent transadmittance between generator 1 and generator 2,
The transient state angular frequency that formula (14) tries to achieve generator 1 and generator 2 when bus 1-bus D circuit A phases head end is overlapped is substituted into,
Transformation relation according to formula (15) between angular frequency and frequency, obtains sending out when bus 1-bus D circuit A phases head end is overlapped
The transient frequency of motor 1 and generator 2.
A phase singlephase earth faults, the tripping operation of 0.1s circuit A phases head end, 0.2s lines occur for bus 1-bus D circuit head ends 0s
Road A phases end is tripped, and 1s circuit A phases end is overlapped, the tripping operation of 1.2s faulty line first and ends three-phase.Faulty line is formed in 1s
When A phases end overlap, head end broken string when positive sequence, negative phase-sequence, zero sequence admittance battle array, to positive sequence, negative phase-sequence, zero sequence admittance battle array invert, obtain
To positive sequence, negative phase-sequence, zero sequence impedance battle array, positive sequence, negative phase-sequence, zero sequence port Impedance battle array are asked for using two-port network theory, cancellation is closed
Close behind loop, obtain comprehensive impedance battle array, comprehensive impedance battle array is inverted, obtain complex admittance battle array, using WARD equivalent methods, cancellation is removed
After the external node of generator 1 and generator 2, the equivalent transadmittance between generator 1 and generator 2 is obtained, formula is substituted into
(14) the transient state angular frequency of generator 1 and generator 2 when bus 1-bus D circuit A phases end is overlapped is tried to achieve, according to formula (15)
Transformation relation between shown angular frequency and frequency, obtains generator 1 and generating when bus 1-bus D circuit A phases end is overlapped
The transient frequency of machine 2.
The transient frequency change curve of generator 1 and generator 2 is carried out when being overlapped to bus 1-bus D circuit A phases head end
Quantitative evaluation (as shown in Figure 3,4), chooses generator 1 and the transient frequency stability margin smaller 74.12% of generator 2 as first
The transient frequency stability margin of system, is shown in table 5 when end is overlapped.Generator 1 when being overlapped to bus 1-bus D circuit A phases end
Quantitative evaluation (as shown in Figure 5,6) is carried out with the transient frequency change curve of generator 2, generator 1 and the transient state of generator 2 is chosen
The transient frequency stability margin of system, is shown in table 5 when frequency stabilization nargin smaller 78.44% overlaps as end.Bus 1-
Bus D circuit A phase head end failures, the transient frequency stability margin of system, chooses the greater when contrast first and end is overlapped
78.44% corresponding end overlaps the time sequence of coincidence scheme as optimization, is shown in table 6.
Embodiment 4:As Figure 1-10 shows,
By taking 2 machine systems shown in Fig. 1 as an example, 8 buses, the node parameter such as table of each generator and bus are had in system
Shown in 1, the voltage class of generator 1 is 18.0kV, and the voltage class of generator 1 is 13.8kV, and busbar voltage grade is 230kV.Line
The circuit per unit value parameter such as resistance, reactance and the admittance on road is as shown in table 2.8 line voltage distribution grade is 230kV.Transformer
The parameter such as resistance, reactance and no-load voltage ratio as shown in table 3, the inertia time constant of generator, the reactance of d-axis transient state, quadrature axis transient state electricity
The parameter such as anti-is as shown in table 4.The power reference value for the parameter per unit value being related in table 2- tables 4 takes 100MVA.
A phase singlephase earth faults occur for 2 machine system busbar 2-bus D circuit head ends 0s shown in Fig. 1, and 0.1s circuit A phases are first
End tripping operation, 0.2s circuit A phases end tripping operation, 1s circuit A phases head end is overlapped, the tripping operation of 1.2s faulty line first and ends three-phase.Formed
Positive sequence, negative phase-sequence, zero sequence admittance battle array when faulty line A phases head end when 1s is overlapped, end is broken, lead to positive sequence, negative phase-sequence, zero sequence
Battle array of receiving is inverted, and obtains positive sequence, negative phase-sequence, zero sequence impedance battle array, and positive sequence, negative phase-sequence, the resistance of zero sequence port are asked for using two-port network theory
Anti- battle array, is eliminated behind closed-loop path, obtains comprehensive impedance battle array, and comprehensive impedance battle array is inverted, and is obtained complex admittance battle array, is utilized WARD etc.
Value method, eliminates after the external node except generator 1 and generator 2, obtains the equivalent transadmittance between generator 1 and generator 2,
The transient state angular frequency that formula (14) tries to achieve generator 1 and generator 2 when bus 2-bus D circuit A phases head end is overlapped is substituted into,
Transformation relation according to formula (15) between angular frequency and frequency, obtains sending out when bus 2-bus D circuit A phases head end is overlapped
The transient frequency of motor 1 and generator 2.
A phase singlephase earth faults, the tripping operation of 0.1s circuit A phases head end, 0.2s lines occur for bus 2-bus D circuit head ends 0s
Road A phases end is tripped, and 1s circuit A phases end is overlapped, the tripping operation of 1.2s faulty line first and ends three-phase.Faulty line is formed in 1s
When A phases end overlap, head end broken string when positive sequence, negative phase-sequence, zero sequence admittance battle array, to positive sequence, negative phase-sequence, zero sequence admittance battle array invert, obtain
To positive sequence, negative phase-sequence, zero sequence impedance battle array, positive sequence, negative phase-sequence, zero sequence port Impedance battle array are asked for using two-port network theory, cancellation is closed
Close behind loop, obtain comprehensive impedance battle array, comprehensive impedance battle array is inverted, obtain complex admittance battle array, using WARD equivalent methods, cancellation is removed
After the external node of generator 1 and generator 2, the equivalent transadmittance between generator 1 and generator 2 is obtained, formula is substituted into
(14) the transient state angular frequency of generator 1 and generator 2 when bus 2-bus D circuit A phases end is overlapped is tried to achieve, according to formula (15)
Transformation relation between shown angular frequency and frequency, obtains generator 1 and generating when bus 2-bus D circuit A phases end is overlapped
The transient frequency of machine 2.
The transient frequency change curve of generator 1 and generator 2 is carried out when being overlapped to bus 2-bus D circuit A phases head end
Quantitative evaluation (as shown in Figure 7,8), chooses generator 1 and the transient frequency stability margin smaller 89.01% of generator 2 as first
The transient frequency stability margin of system, is shown in table 5 when end is overlapped.Generator 1 when being overlapped to bus 2-bus D circuit A phases end
Quantitative evaluation (as shown in Fig. 9,10) is carried out with the transient frequency change curve of generator 2, generator 1 is chosen and generator 2 is temporary
The transient frequency stability margin of system, is shown in table 5 when state frequency stabilization nargin smaller 69.73% overlaps as end.Bus
2-bus D circuit A phase head end failures, the transient frequency stability margin of system, chooses the greater when contrast first and end is overlapped
89.01% corresponding head end overlaps the time sequence of coincidence scheme as optimization, is shown in table 6.
For in above-described embodiment, table 1- tables 6 are as follows:
The machine system node parameter of table 12
The machine system line per unit value parameter of table 22
Note:Power reference value is 100MWA
The machine system transformer parameter of table 32
Note:Power reference value is 100MWA
The machine system generator parameter of table 42
Note:Power reference value is 100MWA
In table 4, TjFor the inertia time constant of generator, unit is s;x′dFor d-axis transient state reactance;x′qFor quadrature axis transient state
Reactance;xdFor the unsaturated synchronous reactance of d-axis;xqFor the unsaturated synchronous reactance of quadrature axis;T′d0For d-axis transient state open circuit time constant;
T′q0For quadrature axis transient state open circuit time constant;xLFor stator leakage reactance;x″dFor d axle subtranient reactances;x″qFor q axle subtranient reactances;
T″d0For d axle time time constants;T″q0For q axle time time constants.
When the transmission line of alternation current single phase ground fault failure of 2 machine systems, generator 1 and generator 2 are in difference
The transient frequency stability margin of system is shown in table 5 under single-phase time sequence of coincidence.Choose the corresponding single-phase time sequence of coincidence of nargin the greater
As the coincidence scheme of optimization, the single-phase time sequence of coincidence scheme of raising transient frequency stability as shown in table 6 is obtained.
The transient frequency stability margin of table 5
Table 6 optimizes time sequence of coincidence scheme
Above in conjunction with accompanying drawing to the present invention embodiment be explained in detail, but the present invention be not limited to it is above-mentioned
Embodiment, can also be before present inventive concept not be departed from the knowledge that those of ordinary skill in the art possess
Put and make a variety of changes.
Claims (2)
1. a kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities, it is characterised in that:When 2
In the presence of two kinds of failures it is single-phase earthing in single-phase coincidence when singlephase earth fault occurs for the transmission line of alternation current of electromechanical Force system
Failure and single-phase wire break failure, form the fault network of dual-port;Three sequence sequences are obtained using symmetrical component method according to fault type
Net, the admittance battle array of each sequence is write according to positive sequence, negative phase-sequence, zero-sequence network row, and the admittance battle array of each sequence is inverted and obtains the Impedance Matrix of each sequence;Adopt
Try to achieve port Impedance battle array first with dual-port theory, and then try to achieve comprehensive impedance battle array, it is inverted and obtains complex admittance battle array;According to
Complex admittance battle array modification positive sequence network, so that the positive sequence admittance battle array that is expanded;2 generator nodes are calculated using Ward equivalent methods
Between equivalent transadmittance, thus calculate the generator frequency variation track under two kinds of reclosing time sequences, realize and frequency is become
Change the quantitative evaluation of track to obtain transient frequency stability margin, choose the corresponding coincidence of transient frequency stability margin the greater
Sequential is the coincidence scheme of optimization.
2. the single-phase time sequence of coincidence setting method according to claim 1 for improving 2 electro-mechanical force system transient modelling frequency stabilities,
It is characterized in that:Methods described is comprised the following steps that:
(1) analog line head end input single-pole reclosing
1) positive sequence, negative phase-sequence, zero sequence admittance battle array are formed
Analog line head end puts into single-pole reclosing, first equivalent to the multiple failure that there occurs single-phase earthing and single-phase end broken string
Positive sequence, negative phase-sequence, zero sequence admittance battle array under the multiple failure are initially formed, Y is designated as1、Y2、Y0;Wherein, the method for forming admittance battle array is as follows:Lead
The number of nonzero element is equal to the connected earth-free circuitry number of corresponding node in each row off-diagonal element of battle array of receiving;Admittance battle array is each diagonal
Element, i.e., the self-admittance Y of each nodeiiEqual to the admittance sum on respective nodes institute's chord road:Admittance battle array is each non-
The transadmittance Y of diagonal elementijIt is equal to connect between 2 nodes the negative value of admittance:Yij=-yij;
In formula, yijFor the admittance on institute's chord road between node i and node j;
2) positive sequence, negative phase-sequence, zero sequence impedance battle array are asked for
To positive sequence, negative phase-sequence, zero sequence admittance battle array Y1、Y2、Y0Invert, as shown in formula (1), obtain positive sequence, negative phase-sequence, zero sequence impedance battle array, note
For Z1、Z2、Z0;
3) port Impedance is asked for using dual-port theory
1. positive sequence, negative phase-sequence, zero sequence self-impedance and the mutual impedance of port 1 and port 2 are asked for;
If corresponding node serial number is m, n at generation singlephase earth fault, n is 0 node of the earth;At single-phase end disconnection fault pair
The node serial number answered is p, q;It is that the corresponding port of m, n is denoted as port 1 by node serial number;Node serial number is the corresponding port of p, q
It is denoted as port 2;
By step 2) in the three sequence Impedance Matrixes that calculate corresponding element substitute into formula (2), try to achieve port 1 and port 2 positive sequence,
Negative phase-sequence, the self-impedance of zero sequence and mutual impedance;
In formula, Z11(1)、Z22(1)The positive sequence self-impedance of port 1 and port 2 is represented respectively;Z12(1)、Z21(1)The He of port 1 is represented respectively
Positive sequence mutual impedance between port 2;Z11(2)、Z22(2)The negative phase-sequence self-impedance of port 1 and port 2 is represented respectively;Z12(2)、Z21(2)Point
Negative phase-sequence mutual impedance that Biao Shi be between port 1 and port 2;Z11(0)、Z22(0)Represent that the zero sequence of port 1 and port 2 hinders certainly respectively
It is anti-;Z12(0)、Z21(0)The zero sequence mutual impedance between port 1 and port 2 is represented respectively;Zuv(1)(u=m, n, p, q;V=m, n, p, q)
Represent Impedance Matrix Z1In u row v column elements, Zuv(0)(u=m, n, p, q;V=m, n, p, q) represent Impedance Matrix Z0In u rows v arrange
Element;
2. the port Impedance battle array of the positive sequence, negative phase-sequence, zero sequence of port 1 and port 2 is sought
Obtained each sequence impedance is calculated by formula (2), positive sequence, negative phase-sequence, zero sequence port Impedance battle array is formed, Z is denoted as respectively(1)、Z(2)、
Z(0), as shown in formula (3);
4) comprehensive impedance battle array Z is asked forF
1. obtained each sequence impedance element is calculated according to by formula (2), row write negative phase-sequence, zero sequence port comprehensive impedance battle array Z, such as formula (4)
It is shown;
2. negative phase-sequence shown in formula (4), zero sequence port comprehensive impedance battle array Z are substituted into formula (5) counter circuit Impedance Matrix ZL;
In formula,CTFor C transposition;Z′sl(s=a, b, c;L=a, b, c), during s=l, represent loop s's
Self-impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;
Formula (5) is 3. eliminated to the impedance element of the corresponding impedance elements of closed-loop path c, retention loop a and loop b by formula (6), obtained
To comprehensive impedance battle array ZF, as shown in formula (7);
In formula, Zsl(s=a, b;L=a, b) during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l, represent
Eliminate the mutual impedance between loop s and loop l behind closed-loop path;Z′sl(s=a, b;L=a, b) during s=l, represent loop s from
Impedance;During s ≠ l, the mutual impedance between loop s and loop l is represented;Z′sc(s=a, it is b) mutual between expression loop s and loop c
Impedance;Z′cl(l=a b) represents the mutual impedance between loop c and loop l;Z′ccRepresent loop c self-impedance;
In formula, Zsl(s=a, b;L=a, b), during s=l, represents to eliminate the self-impedance of loop s behind closed-loop path;During s ≠ l, represent
Eliminate the mutual impedance between loop s and loop l behind closed-loop path;
5) formula (7) is inverted and obtains complex admittance battle array, be denoted as YF, as shown in formula (8)
In formula, ysl(s=a, b;L=a, b), during s=l, represents to eliminate the self-admittance of loop s behind closed-loop path;During s ≠ l, represent
Eliminate the transadmittance between loop s and loop l behind closed-loop path;
6) complex admittance battle array Y is utilizedFIn element modification positive sequence admittance battle array Y1Middle corresponding element, amending method is shown in formula (9), obtains
To extension positive sequence admittance battle array Y1E
In formula, y 'uv(1)(u=m, n;V=m, n) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u=m, n;
V=m, n) represents admittance battle array Y1In u row v column elements;
y′uv(1)(u=p, q;V=p, q) represents extension positive sequence admittance battle array Y1EIn u row v column elements;yuv(1)(u=p, q;V=p,
Q) admittance battle array Y is represented1In u row v column elements;
7) according to extension positive sequence admittance battle array Y1E, the equivalent transadmittance between 2 generator nodes, note are calculated using ward equivalent methods
For Y12;Specific method is as follows:
1. the set constituted generator 1 and the corresponding node of generator 2 is denoted as what remaining node in boundary node B, network was constituted
Set is denoted as external node E, by extension positive sequence admittance battle array Y1EMatrix in block form form is write, the net that matrix in block form is represented can be obtained
Shown in network equation such as formula (10);
In formula, YBBRepresent the admittance battle array corresponding to all boundary nodes;YBERepresent all boundary nodes and all external nodes institute
Corresponding admittance battle array, and YBEEqual to YEBTransposition;YEERepresent the admittance battle array corresponding to all external nodes;Represent all sides
The column voltage vector of boundary's point;Represent the column voltage vector of all external nodes;Represent the injection electricity of all boundary nodes
Stream;Represent the Injection Current of all external nodes;
2. the matrix in block form in formula (10) is substituted into (11) formula and calculated, obtain boundary admittance battle array
In formula,It is the boundary admittance battle array after equivalence, it includes the contribution of the equivalent branch road produced after external network abbreviation;
Therefore the border after the equivalent transadmittance between the node of generator 1 that retains in 2 machine systems and the node of generator 2 is exactly equivalence
Admittance, so there is formula (12) establishment:
In formula, G12For equivalent transadmittance Y12Real part, referred to as conductance, and G12>0;B12For equivalent transadmittance Y12Imaginary part, be referred to as
Susceptance, and B12>0;
8) by step 7) the middle equivalent transadmittance Y calculated between obtained generator 1 and generator 212In substitution formula (13), respectively
Local derviation K of the active power of calculating generator 1 and generator 2 to respective generator rotor angle11、K22;
In formula, PG1、PG2The respectively active power of generator 1 and generator 2;δ1、δ2Respectively generator 1 and generator 2
Generator rotor angle;Represent the active power of generator 1 and generator 2 to the partial derivative of generator's power and angle first respectively
The numerical value at moment beginning;Eq1、Eq2The respectively electromotive force of generator 1 and generator 2;δ0Phase at expression system initial work location
Angle;|Y12| represent equivalent transadmittance Y12Amplitude,γ represents equivalent transadmittance Y12Phase angle, γ=
arctg(B12/G12), and γ ∈ (0, pi/2);
9) by step 8) in calculate obtained K11、K22The transient state angle of 2 machine systems when the single-phase coincidence of head end is calculated in substitution formula (14)
The transient state angular frequency of frequency, generator 1 and generator 2 is denoted as Δ ω respectively1(t)、Δω2(t);
In formula, D1、D2The damped coefficient of generator 1 and generator 2 is represented respectively;ΔPm10、ΔPm20The He of generator 1 is represented respectively
The chugging that generator 2 occurs;M1、M2Represent the rotary inertia of generator 1,2;ω0For
The rated angular velocity of generator, t is the time;
10) because the relation between angular frequency and frequency meets formula (15), by step 9) the middle angle for calculating obtained generator 1,2
Frequency Δ ω1(t)、Δω2(t) formula (16) is substituted into, the transient frequency Δ f of generator 1,2 when calculating head end single-phase coincidence1(t)、
Δf2(t);
The π (15) of f=ω/2
In formula, f is frequency;ω is angular frequency;
(2) analog line end input single-pole reclosing
1) positive sequence, negative phase-sequence, zero sequence admittance battle array are formed
Analog line end puts into single-pole reclosing, first equivalent to the multiple failure that there occurs single-phase earthing and single-phase head end broken string
Positive sequence, negative phase-sequence, zero sequence admittance battle array under the multiple failure are initially formed, Y ' is designated as1、Y′2、Y′0;
2) by Y '1、Y′2、Y′0Substitution formula (1), positive sequence, negative phase-sequence, the zero sequence impedance battle array inverted when obtaining the single-phase coincidence of line end,
It is designated as Z '1、Z′2、Z′0;
3) 3 in repeat step (1)) port Impedance is asked for using dual-port theory, it is designated as Z '(1)、Z′(2)、Z′(0);
4) 4 in repeat step (1)) comprehensive impedance battle array is asked for, it is designated as Z 'F;
5) 5 in repeat step (1)) complex admittance battle array is obtained, it is denoted as Y 'F;
6) 6 in repeat step (1)) utilize complex admittance battle array Y 'FMiddle element modification positive sequence admittance battle array Y '1Middle corresponding element, is obtained
To extension positive sequence admittance battle array Y '1E;
7) 7 in repeat step (1)) according to extension positive sequence admittance battle array Y '1E, 2 generator nodes are calculated using ward equivalent methods
Between equivalent transadmittance, be designated as Y '12;
8) 8 in repeat step (1)) local derviation K ' of the active power to respective generator rotor angle of generator 1 and generator 2 is calculated respectively11、
K′22;
9) 9 in repeat step (1)) the transient state angular frequency of generator 1 and generator 2 is denoted as respectively when calculating the single-phase coincidence in end
Δω′1(t)、Δω′2(t);
10) 10 in repeat step (1)) the transient frequency Δ f ' of generator 1,2 when calculating the single-phase coincidence in end1(t)、Δf′2
(t);
(3) frequency shift (FS) security quantitative evaluation index
Obtained generator transient frequency is calculated to step (1), (2) using formula (17) and carries out quantitative evaluation, circuit head end is obtained
The transient frequency stability margin η of generator 1,2 during single-phase coincidenceG1、ηG2, the transient state of generator 1,2 during line end single-phase coincidence
Frequency stabilization nargin η 'G1、η′G2;
η=[fext-(fcr-kTcr)] × 100% (17)
In formula, fcrAnd TcrIt is the frequency shift (FS) threshold value of generator and the duration allowed, f respectivelyextRefer in transient process and send out
The limiting value of electric machine frequency, k is the commutation factor for critical frequency shift duration being converted into frequency, and η is positive or negative value point
Biao Shi not transient frequency stabilization or unstable;
The transient frequency stability margin of each generator in 2 machine systems is obtained, each generator transient frequency stability margin minimum value is taken
It is used as the transient frequency stability margin of system;
System transient modelling frequency stabilization nargin during the single-phase coincidence of head end is calculated by formula (18), is denoted as η1;
η1=min { ηG1, ηG2} (18)
System transient modelling frequency stabilization nargin during the single-phase coincidence in end is calculated by formula (19), is denoted as η2;
η2=min { η 'G1,η′G2} (19)
(4) draw a conclusion
System transient modelling frequency stabilization nargin η during the single-phase coincidence of alternative route first and end1、η2, choose the greater corresponding single-phase
Time sequence of coincidence as output optimization time sequence of coincidence scheme;If η1>η2, then single-phase coincidence is put into by circuit head end first;Conversely,
Put into single-phase coincidence first by line end.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510740422.7A CN105281329B (en) | 2015-11-04 | 2015-11-04 | A kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510740422.7A CN105281329B (en) | 2015-11-04 | 2015-11-04 | A kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105281329A CN105281329A (en) | 2016-01-27 |
CN105281329B true CN105281329B (en) | 2017-09-22 |
Family
ID=55149877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510740422.7A Active CN105281329B (en) | 2015-11-04 | 2015-11-04 | A kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105281329B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105811410B (en) * | 2016-05-03 | 2018-06-15 | 昆明理工大学 | A kind of single-phase time sequence of coincidence setting method for inhibiting multi-machine power system transient frequency rate of descent |
CN105958429B (en) * | 2016-05-03 | 2018-10-02 | 昆明理工大学 | The single-phase time sequence of coincidence setting method of power frequency fault component distance protection malfunction rate in a kind of reduction ac and dc systems |
CN110086206A (en) * | 2019-03-11 | 2019-08-02 | 昆明理工大学 | A kind of extension one machine infinity bus system Transient angle stability analysis method containing blower |
CN110739700B (en) * | 2019-11-11 | 2023-05-26 | 中国大唐集团科学技术研究院有限公司西北电力试验研究院 | Coincident time sequence optimization method for determining load installation positions of all levels of system |
CN111965478B (en) * | 2020-06-30 | 2023-10-20 | 中国电力科学研究院有限公司 | Method and system for determining short-circuit current by considering phase shifting effect of transformer |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007072074A (en) * | 2005-09-06 | 2007-03-22 | Chugoku Electric Power Co Inc:The | Power failure/instantaneous voltage-down area display system, power failure/instantaneous voltage-down area display method, and power failure/instantaneous voltage-down area display program |
CN100590947C (en) * | 2007-11-08 | 2010-02-17 | 国网南京自动化研究院 | Large electric network on-line preventing control method based on static state and transient safety steady mode |
CN101989766B (en) * | 2009-08-03 | 2012-12-12 | 南京南瑞继保电气有限公司 | Method for diagnosing faults of power system on the basis of protecting action chain |
-
2015
- 2015-11-04 CN CN201510740422.7A patent/CN105281329B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105281329A (en) | 2016-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Du et al. | Modeling of grid-forming and grid-following inverters for dynamic simulation of large-scale distribution systems | |
Kalcon et al. | Small-signal stability analysis of multi-terminal VSC-based DC transmission systems | |
CN105281329B (en) | A kind of single-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities | |
Ma et al. | Analysis on application of a current-source based DFIG wind generator model | |
Jiang et al. | Damping performance analysis of IPFC and UPFC controllers using validated small-signal models | |
Elizondo et al. | Three-phase unbalanced transient dynamics and powerflow for modeling distribution systems with synchronous machines | |
Ranković et al. | A three-phase state estimation in active distribution networks | |
Silva et al. | Impact of distributed generation on the IEEE 34 node radial test feeder with overcurrent protection | |
CN107121604A (en) | A kind of unsymmetrical short-circuit current dc component damping time constant acquisition methods | |
Haddadi et al. | Power system test cases for EMT-type simulation studies | |
Pan et al. | A practical short-circuit current calculation method for DFIG-based wind farm considering voltage distribution | |
CN106026083B (en) | A kind of Power Network Transient Stability analysis method based on MATLAB | |
Erlich et al. | Selective damping of inter area oscillations using phasor measurement unit (PMU) signals | |
Bahravar et al. | Modeling and coordinated control of IPFC and SVC to improve the low frequency oscillation damping | |
CN105305429B (en) | A kind of three-phase time sequence of coincidence setting method for improving 2 electro-mechanical force system transient modelling frequency stabilities | |
CN108344943A (en) | A kind of DFIG low voltage ride-through testing methods and device | |
Jain et al. | Studying the impact of solar PV on power system dynamics using integrated transmission and distribution network models | |
Godoy Enciso et al. | Frequency response analysis under faults in weak power systems | |
Su et al. | A new approach for integration of two distinct types of numerical simulator | |
CN105811410B (en) | A kind of single-phase time sequence of coincidence setting method for inhibiting multi-machine power system transient frequency rate of descent | |
Siddique et al. | Analysis of Transient Stability with SSSC and UPFC with Multi-Band PSS in Two Area Multi-Machine System | |
Mohammed et al. | Structure preserving model analysis for power system transient stability | |
Annamalai et al. | Development of AC network dynamic equivalents for large scale power system | |
Cañizares et al. | Angle, voltage, and frequency stability | |
Alves et al. | Comparative Study of Different Values of Inertia for Islanded Systems with Energy Storage Systems |
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 |