CN103795058B - The air extract analysis of electric power system and system failure sort method - Google Patents
The air extract analysis of electric power system and system failure sort method Download PDFInfo
- Publication number
- CN103795058B CN103795058B CN201410049147.XA CN201410049147A CN103795058B CN 103795058 B CN103795058 B CN 103795058B CN 201410049147 A CN201410049147 A CN 201410049147A CN 103795058 B CN103795058 B CN 103795058B
- Authority
- CN
- China
- Prior art keywords
- node
- voltage
- fault
- stability margin
- theta
- 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
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 238000012544 monitoring process Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000012937 correction Methods 0.000 claims description 15
- 239000003990 capacitor Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 9
- 238000012804 iterative process Methods 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000006698 induction Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 claims description 3
- 230000007423 decrease Effects 0.000 abstract description 4
- 239000000243 solution Substances 0.000 description 11
- 230000006641 stabilisation Effects 0.000 description 6
- 238000011105 stabilization Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention proposes a kind of air extract analysis and system failure sort method of electric power system, comprising: based on the Newton iteration method determination quiescent voltage collapse point of Optimal Multiplier; The type of voltage collapse point is judged according to the characteristic of iteration convergence; Require checking system Failure risk situation according to stability margin and provide the sequence of stable fault; By fault parameter and use iterative method provide de-stabilise fault the order of severity sequence.The present invention can provide the voltage stability margin of electric power system online more quickly, carries out effective monitoring in real time to voltage stability; Stable fault and de-stabilise fault integrated ordered can when system jam the switching of the adjustment of online direction generator and reactive-load compensation equipment, also can instruct the configuration etc. of branch parameters adjustment, circuit increase and decrease, FACTS equipment by off-line, all have important meaning to the operation of electric power system and planning.
Description
Technical field
The present invention relates to the technical field of power system monitoring, analysis and control, the air extract analysis of specifically a kind of electric power system and system failure sort method.
Background technology
The Voltage-stabilizing Problems research starting of electric power system is later than the research of frequency stabilization problem, be started in the forties in last century, until after last century the seventies, Voltage-stabilizing Problems just starts to receive publicity as a special field, along with the development of modern power systems, the development of field of power transmission may lag behind the growth of system loading growth and power generation level, under electric power networks often operates in high load capacity level, face the critical restriction of voltage stabilization, cause a series of great voltage stabilization accident.
Static voltage stability problem has many analytical methods.Utilize system PV curve, be given in the voltage stability margin of load bus on certain load growing direction, it is the static voltage stability analysis method commonly used the most, the method has clear and definite physical background, analytical effect is directly perceived, has the analytical method of the maturations such as continuous tide, as shown in Figure 1, voltage stability margin λ is defined by following formula
In Fig. 1, P
maxcrefer to the active power level of system during voltage collapse point, P
initialrefer to system initial launch point active power level, P
maxIrefer to that system cloud gray model meets the maximum reactive power level of system under certain margin requirement, V
irefer to that system cloud gray model meets lowest voltage of a system level under certain margin requirement, V
crefer to voltage levvl during voltage collapse point.
The system failure such as the exiting of injection type equipment such as generator, parallel reactive compensator can produce direct impact with the Voltage Stability Level exited for system of the Branch Type such as system line, transformer equipment, the generation of generic failure all can cause the reduction of voltage stability margin λ, if when λ is greater than 0, can think that post-fault system still keeps stable, if λ is less than 0, then post-fault system will lose stable, be called unstability fault.
Voltage stability margin real-time assessment after Chinese invention patent (application number 201010140847.1) power transmission network catastrophe failure and method for optimally controlling, by judging that the solution point calculated in iterative trend step is true to separate or optimal solution carrys out the voltage stability of system after failure judgement, damped Newton method is utilized to solve, add an iterative process, recycling continuous tide solves stability margin, iterations is more, the time is longer, and does not have to realize the sequence to fault severity level.
According to simulation process information, Chinese invention patent (application number 201110368041.2) a kind of method for evaluating severity of power system fault, judges whether system there occurs merit angle, voltage, frequency unstability; The serious coefficient of this fault is had respectively according to merit angle, voltage, frequency three index calculate, aforementioned three the serious coefficients obtained are weighted combination and draw comprehensive serious coefficient, this fault severity level sort method relates to voltage stabilization and frequency stabilization, be not solve voltage stabilization contingency ranking name problem pointedly, and the foundation of contingency ranking is fail result, that fault is sorted on the impact that system produces, the complexity solving fault is not discussed, the traffic control method solving fault cannot be provided simultaneously.
Summary of the invention
The invention provides a kind of air extract analysis and system failure sort method of electric power system, the voltage stability margin of system can be provided online more quickly, effective monitoring is in real time carried out to the voltage stability of electric power system; Stable fault and de-stabilise fault integrated ordered can when system jam the switching of the adjustment of online direction generator output and reactive-load compensation equipment, also can instruct the configuration etc. of branch parameters adjustment, circuit increase and decrease, FACTS equipment by off-line, all have important meaning to the operation of electric power system and planning.
The technology used in the present invention means are as follows:
The air extract analysis of electric power system and system failure sort method, comprise the following steps:
(1) according to the monitoring result to electric power system, judge whether electric power system breaks down, if electric power system is not broken down, namely under normal operating mode, then carry out the Newton iteration method based on Optimal Multiplier from feasible zone outside, ask for the voltage stability margin of electric power system on predetermined load growing direction and determine voltage collapse critical point, then proceeding to step 2); If electric power system is broken down, then directly enter step 3);
(2) according to Newton iteration method convergence situation analysis voltage stability margin, and voltage collapse vertex type is distinguished;
(3) if electric power system is broken down, still according to the Newton iteration method of carrying out from feasible zone outside based on Optimal Multiplier, ask for voltage stability margin and determine voltage collapse critical point, in an iterative process, the voltage collapse critical point nargin asked under selecting described step (1) normal operating mode of iteration initial point;
(4) analyze the voltage stability margin that described step (3) is asked for, if voltage stability margin is not less than 0, be then stable fault, carry out order of severity sequence according to voltage stability margin size to stable fault, voltage stability margin value is less, and fault is more serious; If voltage stability margin is less than 0, be then unstability fault, first parametrization carried out to fault, then ask for voltage collapse critical point according to the Newton iteration method based on Optimal Multiplier from feasible zone outside, in an iterative process, voltage stability margin when iteration initial point selects fault parameter to be 0;
(5) calculate the voltage collapse critical point fault parameter that described step (4) is asked for, carry out order of severity sequence according to fault parameter size to unstability fault, fault parameter is larger, and fault is more serious;
(6) two kinds of sequences of combining step (4) and step (5), according to voltage stability margin value from big to small, then fault parameter order from small to large, to unified sequence of being out of order, instructs power system operation.
In aforesaid step (1), carry out the Newton iteration method based on Optimal Multiplier from feasible zone outside, ask for voltage stability margin and determine that the concrete grammar of voltage collapse critical point is:
The growing direction of load and generator is defined by following formula:
P
Li=P
Li0+λb
Pi
Q
Li=Q
Li0+λb
Qi
P
Gi=P
Gi0+λb
Gi
Wherein, λ is voltage stability margin, P
li0, Q
li0be respectively active power and reactive power that node i injects under normal condition, P
li, Q
libe respectively active power and reactive power that node i injects under current state, P
gi0, P
gifor the active power that node i generator injects under normal condition and under current state, b
pi, b
qi, b
gibe respectively that the load of node i is meritorious exerts oneself, idlely exert oneself and the change direction vector of generator output;
Power flow equation with parameter is expressed as:
f(x,λ)=f(x)-S=0
S=S
0+λb
Wherein, S
0, S is respectively under normal condition and current state lower node and generator injecting power vector, S
0=(P
li0, Q
li0, P
gi0), S=(P
li, Q
li, P
gi), b is node and generator injecting power change direction vector, b=(b
pi, b
qi, b
gi), x is state variable;
Adopt the Newton iteration method of Optimal Multiplier, iteration initial point is chosen and is met the voltage stability margin value of trend outside feasible zone, obtains state variable x in trend kth time iteration
(k)correction amount x
(k), Δ x
(k)=J
(k)-1f (x
(k)), wherein, J
(k)for the Jacobian matrix of kth time iteration, f (x
(k)) for kth time iteration is obtained state variable x
(k)substitute into power flow equation group, the concrete form of Jacobian matrix is:
Correction amount x is multiplied by with a scalar multiplier β
(k), then revise state variable x
(k), the sub-β of its Scalar Multiplication is tried to achieve by following target function:
F
i() represents i-th equation in equation group f (x, λ)=f (x)-S=0, and 2n is the number of equation, obtains scalar multiplier β by asking for F (β) extreme value, equation as shown in the formula:
The state variable x that scalar multiplier β is corresponding when being 0
*for the least square solution of power flow equation, voltage stability margin λ when β is 0
criticalbe the margin value that voltage collapse critical point is corresponding.
The specific implementation process of aforesaid step (2) is:
If in the iteration of step (1), occur that the situation of PV/PQ type conversion occurs certain node repeatedly, then the type of voltage collapse point is constraint induction type, and this node is voltage collapse point;
If in the iteration of step (1), occur that the situation of PV/PQ type conversion occurs certain several node repeatedly, then the type of voltage collapse point is constraint induction type, record these nodes, select 1 node i at every turn, calculate trend by described step (1), after trend convergence, according to final Jacobian matrix meter sensitivity, if meet the following conditions:
Then node i is voltage collapse point, wherein, and V
gifor node i generator voltage, Q
gifor the reactive power that node i generator injects,
If in the iteration of step (1), the correction of voltage stability margin is less than default precision, then the type of voltage collapse point is saddle junction type, after then obtaining least square solution, need to revise voltage stability margin λ further, search voltage collapse point, namely search the voltage stability margin λ of corresponding voltage collapse point
critical, be specially:
Definition ∑ is have between power flow equation to separate and without the border of separating between region, node clean injecting power vector S forms a space, S
λ, S ', S
mthree is vector power in space, and S ' is current power vector, definition S
mfor ∑ is at state variable x
*point nearest apart from current power vector S ' Euclidean distance on the section at place, definition S
λfor section and node injecting power change direction vector b intersection point, S
λwith S
criticaloverlap, S
criticalrepresent the vector power that voltage collapse critical point is corresponding, the correction amount λ of definition voltage stability margin is S ' and S
λbetween load parameter difference, then S
λcan following formula be expressed as:
S
λ=S′-Δλb
Wherein, θ
1represent vector power S
mwith the power flow equation f (x at least square solution place
*) between angle, θ
2represent vector power S
mand S
λbetween angle,
When ∑ is convex surface, the amendment type of voltage stability margin λ is:
λ
(k+1)=λ
(k)-Δλ
(k+1)
When ∑ is concave curved surface,
If Δ λ
(k)be greater than default precision, then the amendment type of voltage stability margin λ is:
If θ
1=90 ° or θ
2=90 ° and θ
1when ≠ 90 °, then according to the correction amount λ of following formula calculating voltage stability margin:
The amendment type of voltage stability margin λ is: λ
(k+1)=λ
(k)-Δ λ
(k+1).
Aforesaid step (5) is carried out parametrization to fault and is referred to, if reflection system failure parameters is μ, and parameter area is:
0<μ<1
Fault parameter μ represents that fault does not occur when being 0, represents that fault thoroughly occurs when being 1.
In aforesaid step (5), for different faults, the calculating of fault parameter μ is obtained by following parametrization power flow equation:
(1) the parametrization power flow equation that exits of single generator
(2) the parametrization power flow equation that exits of single shunt capacitor or reactor
(3) the parametrization power flow equation that exits of single load
(4) the parametrization power flow equation that exits of single branch road
If there is multiple faults, then parametrization power flow equation is the linear superposition of each single fault parameter power flow equation in this multiple faults;
Wherein, P
dithe active power of the load absorption of node i; P
gifor the active power that node i generator injects; Q
gifor the reactive power that node i generator injects; Q
difor the reactive power of node i load absorption; Q
sifor shunt capacitor capacity, Q
gmaxU, Q
gminUfor the upper and lower limit that generator reactive exports; Q
rifor the capacity of reactive-load compensation capacitor after fault; V
ifor the voltage magnitude of node i; I is all node set; θ
ijfor the phase angle difference between node i, j; B
ijfor the susceptance between admittance matrix interior joint i, j; G
ijfor the conductance between admittance matrix interior joint i, j; G
iifor the self-conductance of node i; B
iifor node i from susceptance; G
iinewfor branch road i-m break down after self-conductance in system admittance matrix; B
iinewfor branch road i-m break down after in system admittance matrix from susceptance.
By adopting above-mentioned technological means, the beneficial effect that the present invention has is:
1) operation is simple, and step is clear, and the data that on-line analysis can be obtained by system SCADA or PMU carry out Direct Analysis, and off-line analysis completes easily through emulation tool;
2) calculate simply, fast, use avoiding continuous tide based on Load Flow Method with Optimal Multiplier iterative analysis method and solve voltage collapse critical point close to situation about not easily restraining during collapse point from feasible zone outside, trend iteration initial point is clear and definite, decreases iterations, affects few by system scale;
3) parametrization is carried out to system jam, utilizing does not have the fault parameter between the 0-1 of physical significance to carry out fault severity level sequence, parameter for sorting is not only the sequence to the failure effect order of severity, mostly ordering scenario being fault generating process or on-line operation and solving difficulty;
4) by fault severity level, by after fault, whether unstability is divided into two large classes, and in two steps two class steps are sorted respectively, the ranking results carried out respectively has carried out unified Ordination, unified ranking results can the operation of guidance system and planning, by preconsolidation stress and on-the-spot operation, reduce the order of severity sequence that may break down.
Accompanying drawing explanation
Fig. 1 is Continuation power flow and voltage stability margin definition schematic diagram;
Fig. 2 is the schematic diagram carried out from feasible zone based on Optimal Multiplier Newton iteration;
Fig. 3 is the iteration geometrical model schematic diagram looking for voltage collapse point from area of feasible solutions;
Fig. 4 is the unified sequence schematic diagram describing stable fault and unstability fault;
Fig. 5 is air extract analysis and the system failure sort method flow chart of electric power system of the present invention.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail.
The present invention utilizes the Optimal Multiplier Newton iteration method of carrying out iteration from feasible zone outside, the method iterative target is for finding voltage collapse critical point, system voltage stabilizes nargin can be analyzed, this method avoid conventional continuous tide to calculate in the problem close to the convergence difficulties produced during collapse point, decrease iterations, improve solving speed, also can according to convergence characteristic identification voltage collapse type.
As shown in Figure 5, the inventive method comprises the following steps:
1, first according to the monitoring result to electric power system, judge whether electric power system breaks down, if electric power system is not broken down, namely under normal operating mode, then carry out the Newton iteration method based on Optimal Multiplier from feasible zone outside, ask for the voltage stability margin of electric power system on predetermined load growing direction and determine voltage collapse critical point, then proceeding to step 2; If electric power system is broken down, then directly enter step 3;
Under normal operating mode,
The growing direction of load and generator can be defined by following formula:
P
Li=P
Li0+λb
Pi
Q
Li=Q
Li0+λb
Qi
P
Gi=P
Gi0+λb
Gi
Wherein, λ is voltage stability margin, P
li0, Q
li0be respectively active power and reactive power that node i injects under normal condition, P
li, Q
libe respectively active power and reactive power that node i injects under current state, P
gi0, P
gifor the active power that node i generator injects under normal condition and under current state, b
pi, b
qi, b
gibe respectively that the load of node i is meritorious exerts oneself, idlely exert oneself and the change direction vector of generator output.Above node power can substitute into power flow equation,
Power flow equation with parameter can be expressed as:
f(x,λ)=f(x)-S=0
S=S
0+λb
Wherein, S
0, S is respectively under normal condition and current state lower node and generator injecting power vector, S
0=(P
li0, Q
li0, P
gi0), S=(P
li, Q
li, P
gi), b is node and generator injecting power change direction vector, b=(b
pi, b
qi, b
gi), x is state variable.
As shown in Figure 2, adopt the Newton iteration method of Optimal Multiplier, iteration initial point is chosen and is met the voltage stability margin value λ of trend outside feasible zone
0, along approaching direction as shown in the figure, find the nargin λ that voltage collapse critical point is corresponding
critical;
Load Flow Method with Optimal Multiplier utilizes the characteristic that power equation under rectangular coordinate is quadratic function, obtains the correction amount x of state variable in trend kth time iteration
(k), Δ x
(k)=J
(k)-1f (x
(k)), wherein, J
(k)for the Jacobian matrix of kth time iteration, f (x
(k)) for kth time iteration is obtained state variable x
(k)substitute into power flow equation group, the concrete form of Jacobian matrix is:
Correction amount x is multiplied by with a scalar multiplier β
(k), then revise state variable x
(k), the sub-β of its Scalar Multiplication is tried to achieve by following target function:
F
i() represents i-th equation in equation group f (x, λ)=f (x)-S=0, and 2n is the number of equation, obtains scalar multiplier β by asking for F (β) extreme value, equation as shown in the formula:
Before finding collapse critical point, carry out with iteration, β is more and more less, until be 0, now target function maintains one on the occasion of upper, corresponding state variable x
*for the least square solution of trend, corresponding Jacobian matrix J (x
*) unusual.
2, according to Newton iteration method convergence situation analysis voltage stability margin, and voltage collapse vertex type is distinguished;
If in the iteration of step 1, occur that the situation of PV/PQ type conversion occurs certain node repeatedly, this node is idle, and units limits causes system generation voltage collapse active constraint, and constraint induction type voltage collapse occurs, and this node is voltage collapse point;
If in the iteration of step 1, occur that the situation of PV/PQ type conversion occurs certain several node repeatedly, its mechanism and single node are changed similar repeatedly, then record these nodes, select 1 node i at every turn, calculate trend by described step 1, after trend convergence, according to final Jacobian matrix meter sensitivity, if meet the following conditions:
Then node i is the node causing constraint induction type voltage collapse, wherein, and V
gifor node i generator voltage, Q
gifor the reactive power that node i generator injects,
If in the iteration of step 1, the correction of voltage stability margin is less than default precision, default precision is depending on counting accuracy requirement, 0.001 is got in the present invention, then the type of voltage collapse point is saddle junction type, then, after obtaining least square solution, need to revise voltage stability margin λ further, search voltage collapse point, namely search the voltage stability margin λ of corresponding voltage collapse point
critical, be specially:
As shown in Figure 3, definition ∑ is have between power flow equation to separate and without the border of separating between region, the least square solution x of gained
*meet (1) f (x
*) be positioned on ∑, corresponding Jacobian matrix J (x
*) unusual; (2) J (x
*) left eigenvector ω corresponding to zero eigenvalue
*with ∑ at f (x
*) place is orthogonal.Node clean injecting power vector S forms a space, S
λ, S ', S
mthree is vector power in space, and definition S ' is current power vector, definition S
mfor ∑ is at state variable x
*point nearest apart from current power vector S ' Euclidean distance on the section at place, definition S
λfor section and node injecting power change direction vector b intersection point, S
λwith S
criticaloverlap, S
criticalrepresent the vector power that voltage collapse critical point is corresponding, the correction amount λ of definition voltage stability margin is S ' and S
λbetween load parameter difference, then S
λcan following formula be expressed as:
S
λ=S′-Δλb
Wherein, θ
1represent vector power S
mwith the power flow equation f (x at least square solution place
*) between angle, θ
2represent vector power S
mand S
λbetween angle,
When ∑ is convex surface, the amendment type of voltage stability margin λ is:
λ
(k+1)=λ
(k)-Δλ
(k+1)
When ∑ is concave curved surface,
If Δ λ
(k)be greater than default precision, then the amendment type of voltage stability margin λ is:
If θ
1=90 ° or θ
2=90 ° and θ
1when ≠ 90 °, then according to the correction amount λ of following formula calculating voltage stability margin:
The amendment type of voltage stability margin λ is: λ
(k+1)=λ
(k)-Δ λ
(k+1).
If 3 electric power systems are broken down, still according to the Newton iteration method of carrying out from feasible zone outside based on Optimal Multiplier, ask for voltage stability margin and determine voltage collapse critical point, in an iterative process, the voltage collapse critical point nargin that iteration initial point is asked under selecting described step 1 normal operating mode, i.e. λ
critical;
4, the voltage stability margin asked for of analytical procedure 3, if voltage stability margin is not less than 0, be then stable fault, carry out order of severity sequence according to voltage stability margin size to stable fault, voltage stability margin value is less, and fault is more serious; If voltage stability margin is less than 0, is then unstability fault, first carries out parametrization to fault, if reflection system failure parameters is μ, and parameter area is: 0 < μ < 1,
Fault parameter μ represents that fault does not occur when being 0, represent that fault thoroughly occurs when being 1, for a unstability fault, thoroughly there is rear power flow equation without solution in fault, therefore passes through from feasible zone based on the Newton iteration method of Optimal Multiplier, this load flow feasible region outside from μ=0, namely in an iterative process, voltage stability margin when iteration initial point selects fault parameter to be 0, along with μ increases gradually, find voltage collapse critical point, corresponding μ value is namely as the parameter of contingency ranking.μ value is larger, and fault is more serious, the most catastrophe failure of μ=1 correspondence.
5, calculation procedure 4) the voltage collapse critical point fault parameter asked for, carry out order of severity sequence according to fault parameter size to unstability fault, fault parameter is larger, and fault is more serious; The calculating of fault parameter μ is obtained by following parametrization power flow equation:
(1) the parametrization power flow equation that exits of single generator
(2) the parametrization power flow equation that exits of single shunt capacitor or reactor
(3) the parametrization power flow equation that exits of single load
(4) the parametrization power flow equation that exits of single branch road
When there is multiple faults, the system parameters power flow equation of multiple compound contingency is the linear superposition of above several situation, only adopts a parameter μ, and multiple faults occurs with certain parameter level, and weighs,
Wherein, P
dithe active power of the load absorption of node i; P
gifor the active power that node i generator injects; Q
gifor the reactive power that node i generator injects; Q
difor the reactive power of node i load absorption; Q
sifor shunt capacitor capacity, Q
gmaxU, Q
gminUfor the upper and lower limit that generator reactive exports; Q
rifor the capacity of reactive-load compensation capacitor after fault; V
ifor the voltage magnitude of node i; I is all node set; θ
ijfor the phase angle difference between node i, j; B
ijfor the susceptance between admittance matrix interior joint i, j; G
ijfor the conductance between admittance matrix interior joint i, j; G
iifor the self-conductance of node i; B
iifor node i from susceptance; G
iinewfor branch road i-m break down after self-conductance in system admittance matrix; B
iinewfor branch road i-m break down after in system admittance matrix from susceptance.
6, as shown in Figure 4, two kinds of sequences of combining step 4 and step 5, according to voltage stability margin value from big to small, then fault parameter order from small to large, to unified sequence of being out of order, instructs power system operation.
Claims (5)
1. the air extract analysis of electric power system and system failure sort method, is characterized in that, comprise the following steps:
(1) according to the monitoring result to electric power system, judge whether electric power system breaks down, if electric power system is not broken down, namely under normal operating mode, then carry out the Newton iteration method based on Optimal Multiplier from feasible zone outside, ask for the voltage stability margin of electric power system on predetermined load growing direction and determine voltage collapse critical point, then proceeding to step 2); If electric power system is broken down, then directly enter step 3);
(2) according to Newton iteration method convergence situation analysis voltage stability margin, and voltage collapse vertex type is distinguished;
(3) if electric power system is broken down, still according to the Newton iteration method of carrying out from feasible zone outside based on Optimal Multiplier, ask for voltage stability margin and determine voltage collapse critical point, in an iterative process, the voltage stability margin asked under selecting described step (1) normal operating mode of iteration initial point and voltage collapse critical point;
(4) voltage stability margin that described step (3) is asked for is analyzed, if voltage stability margin is not less than 0, is then stable fault, according to voltage stability margin size, order of severity sequence is carried out to stable fault, voltage stability margin value is less, and fault is more serious; If voltage stability margin is less than 0, be then unstability fault, first parametrization carried out to fault, then ask for voltage collapse critical point according to the Newton iteration method based on Optimal Multiplier from feasible zone outside, in an iterative process, voltage stability margin when iteration initial point selects fault parameter to be 0;
(5) calculate the voltage collapse critical point fault parameter that described step (4) is asked for, carry out order of severity sequence according to fault parameter size to unstability fault, fault parameter is larger, and fault is more serious;
(6) two kinds of sequences of combining step (4) and step (5), according to voltage stability margin value from big to small, then fault parameter order from small to large, to unified sequence of being out of order, instructs power system operation.
2. the air extract analysis of electric power system according to claim 1 and system failure sort method, it is characterized in that, in described step (1), carry out the Newton iteration method based on Optimal Multiplier from feasible zone outside, ask for voltage stability margin and determine that the concrete grammar of voltage collapse critical point is:
The growing direction of load and generator is defined by following formula:
P
Li=P
Li0+λb
Pi
Q
Li=Q
Li0+λb
Qi
P
Gi=P
Gi0+λb
Gi
Wherein, λ is voltage stability margin, P
li0, Q
li0be respectively active power and reactive power that node i injects under normal condition, P
li, Q
libe respectively active power and reactive power that node i injects under current state, P
gi0, P
gifor the active power that node i generator injects under normal condition and under current state, b
pi, b
qi, b
gibe respectively that the load of node i is meritorious exerts oneself, idlely exert oneself and the change direction vector of generator output;
Power flow equation with parameter is expressed as:
f(x,λ)=f(x)-S=0
S=S
0+λb
Wherein, S
0, S is respectively under normal condition and current state lower node and generator injecting power vector, S
0=(P
li0, Q
li0, P
gi0), S=(P
li, Q
li, P
gi), b is node and generator injecting power change direction vector, b=(b
pi, b
qi, b
gi), x is state variable;
Adopt the Newton iteration method of Optimal Multiplier, iteration initial point is chosen and is met the voltage stability margin value of trend outside feasible zone, obtains state variable x in trend kth time iteration
(k)correction amount x
(k),
Δ x
(k)=J
(k)-1f (x
(k)), wherein, J
(k)for the Jacobian matrix of kth time iteration, f (x
(k)) for kth time iteration is obtained state variable x
(k)substitute into power flow equation group, the concrete form of Jacobian matrix is:
Correction amount x is multiplied by with a scalar multiplier β
(k), then revise state variable x
(k), the sub-β of its Scalar Multiplication is tried to achieve by following target function:
F
i() represents i-th equation in equation group f (x, λ)=f (x)-S=0, and 2n is the number of equation, obtains scalar multiplier β by asking for F (β) extreme value, equation as shown in the formula:
The state variable x that scalar multiplier β is corresponding when being 0
*for the least square solution of power flow equation, voltage stability margin λ when β is 0
criticalbe the margin value that voltage collapse critical point is corresponding.
3. the air extract analysis of electric power system according to claim 1 and system failure sort method, it is characterized in that, the specific implementation process of described step (2) is:
If in the iteration of step (1), occur that the situation of PV/PQ type conversion occurs certain node repeatedly, then the type of voltage collapse point is constraint induction type, and this node is voltage collapse point;
If in the iteration of step (1), occur that the situation of PV/PQ type conversion occurs certain several node repeatedly, then the type of voltage collapse point is constraint induction type, record these nodes, select 1 node i at every turn, calculate trend by described step (1), after trend convergence, according to final Jacobian matrix meter sensitivity, if meet the following conditions:
Then node i is voltage collapse point, wherein, and V
gifor node i generator voltage, Q
gifor the reactive power that node i generator injects,
If in the iteration of step (1), the correction of voltage stability margin is less than default precision, then the type of voltage collapse point is saddle junction type, after then obtaining least square solution, need to revise voltage stability margin λ further, search voltage collapse point, namely search the voltage stability margin λ of corresponding voltage collapse point
critical, be specially:
Definition Σ has between power flow equation to separate and without the border of separating between region, node clean injecting power vector S forms a space, S
λ, S', S
mthree is vector power in space, and S' is current power vector, definition S
mfor Σ is at state variable x
*point nearest apart from current power vector S' Euclidean distance on the section at place, definition S
λfor section and node injecting power change direction vector b intersection point, S
λwith S
criticaloverlap, S
criticalrepresent the vector power that voltage collapse critical point is corresponding, the correction amount λ of definition voltage stability margin is S' and S
λbetween load parameter difference, then S
λcan following formula be expressed as:
S
λ=S'-Δλb
Wherein, θ
1represent vector power S
mwith the power flow equation f (x at least square solution place
*) between angle, θ
2represent vector power S
mand S
λbetween angle,
When Σ is convex surface, the amendment type of voltage stability margin λ is:
λ
(k+1)=λ
(k)-Δλ
(k+1)
When Σ is concave curved surface,
If Δ λ
(k)be greater than default precision, then the amendment type of voltage stability margin λ is:
If θ
1=90 ° or θ
2=90 ° and θ
1when ≠ 90 °, then according to the correction amount λ of following formula calculating voltage stability margin:
The amendment type of voltage stability margin λ is: λ
(k+1)=λ
(k)-Δ λ
(k+1).
4. the air extract analysis of electric power system according to claim 1 and system failure sort method, it is characterized in that, described step (5) is carried out parametrization to fault and is referred to, if reflection system failure parameters is μ, and parameter area is:
0<μ<1
Fault parameter μ represents that fault does not occur when being 0, represents that fault thoroughly occurs when being 1.
5. the air extract analysis of electric power system according to claim 1 and system failure sort method, it is characterized in that, in described step (5), for different faults, the calculating of fault parameter μ is obtained by following parametrization power flow equation:
(1) the parametrization power flow equation that exits of single generator
(2) the parametrization power flow equation that exits of single shunt capacitor or reactor
(3) the parametrization power flow equation that exits of single load
(4) the parametrization power flow equation that exits of single branch road
If there is multiple faults, then parametrization power flow equation is the linear superposition of each single fault parameter power flow equation in this multiple faults;
Wherein, P
dithe active power of the load absorption of node i; P
gifor the active power that node i generator injects; Q
gifor the reactive power that node i generator injects; Q
difor the reactive power of node i load absorption; Q
sifor shunt capacitor capacity, Q
gmaxU, Q
gminUfor the upper and lower limit that generator reactive exports; Q
rifor the capacity of reactive-load compensation capacitor after fault; V
ifor the voltage magnitude of node i; I is all node set; θ
ijfor the phase angle difference between node i, j; B
ijfor the susceptance between admittance matrix interior joint i, j; G
ijfor the conductance between admittance matrix interior joint i, j; G
iifor the self-conductance of node i; B
iifor node i from susceptance; G
iinewfor branch road i-m break down after self-conductance in system admittance matrix; B
iinewfor branch road i-m break down after in system admittance matrix from susceptance.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410049147.XA CN103795058B (en) | 2014-02-12 | 2014-02-12 | The air extract analysis of electric power system and system failure sort method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410049147.XA CN103795058B (en) | 2014-02-12 | 2014-02-12 | The air extract analysis of electric power system and system failure sort method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103795058A CN103795058A (en) | 2014-05-14 |
CN103795058B true CN103795058B (en) | 2015-10-28 |
Family
ID=50670497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410049147.XA Active CN103795058B (en) | 2014-02-12 | 2014-02-12 | The air extract analysis of electric power system and system failure sort method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103795058B (en) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104198854B (en) * | 2014-08-29 | 2016-09-28 | 国家电网公司 | A kind of online Branch Type Contingency screening and ranking method |
CN104617575B (en) * | 2015-01-28 | 2017-01-04 | 国电南瑞科技股份有限公司 | A kind of trend intangibility appraisal procedure under static security analysis forecast failure |
CN105068785B (en) * | 2015-04-22 | 2018-04-10 | 清华大学 | A kind of parallel calculating method and system |
CN104917174B (en) * | 2015-06-08 | 2017-05-03 | 北京交通大学 | Static voltage stability judgment method for active power distribution network |
CN105137244B (en) * | 2015-09-10 | 2018-06-19 | 南方电网科学研究院有限责任公司 | Method and system for realizing stable monitoring of static voltage in power system |
CN105406480A (en) * | 2015-12-31 | 2016-03-16 | 国网冀北电力有限公司电力科学研究院 | Voltage stability preventive control method and device |
CN106786548B (en) * | 2017-01-12 | 2019-05-17 | 天津大学 | A kind of optimal Reclosing Instant fast algorithm based on Dynamic Security Region |
CN106992516B (en) * | 2017-04-18 | 2020-02-14 | 国网上海市电力公司 | Probability static voltage stability margin optimization implementation method |
CN107147114B (en) * | 2017-06-06 | 2019-08-13 | 广西电网有限责任公司电力科学研究院 | A kind of method for visualizing for realizing power system steady state voltage stability domain |
CN107492888B (en) * | 2017-09-19 | 2019-08-20 | 河海大学 | A kind of voltage weakness zone appraisal procedure of direct current receiving end AC system |
CN107565571B (en) * | 2017-09-21 | 2019-10-08 | 中国农业大学 | A kind of method and device judging power system steady state voltage stability |
CN107843808B (en) * | 2017-10-17 | 2020-04-07 | 科华恒盛股份有限公司 | Power failure detection method and device for trailing of power grid |
CN109038572B (en) * | 2018-08-30 | 2021-05-07 | 国家电网公司华东分部 | Voltage safety range determination method considering faults of alternating current-direct current hybrid power grid |
CN109861214B (en) * | 2019-02-28 | 2022-10-04 | 中国电力科学研究院有限公司 | Method and system for judging weak line with stable transient power angle of regional power grid |
CN110224411B (en) * | 2019-07-12 | 2023-01-24 | 国网山东省电力公司临沂供电公司 | Point selection and volume fixing method for reactive voltage emergency coordination control system |
CN112557813B (en) * | 2020-11-26 | 2022-07-19 | 国网宁夏电力有限公司电力科学研究院 | Method for judging voltage stability of power grid under simultaneous fault of multiple loops of direct current |
CN112600192A (en) * | 2020-11-26 | 2021-04-02 | 国网北京市电力公司 | Fault analysis method and device, electronic equipment and computer readable storage medium |
CN112632453A (en) * | 2020-11-27 | 2021-04-09 | 上海电机学院 | Method, system, equipment and storage medium for calculating voltage breakdown point |
CN112990958A (en) * | 2021-01-19 | 2021-06-18 | 腾讯科技(深圳)有限公司 | Data processing method, data processing device, storage medium and computer equipment |
CN113937823B (en) * | 2021-09-30 | 2023-06-09 | 山东理工大学 | Static voltage stability enhancement control method adopting minimum number of generators to adjust |
CN116865272A (en) * | 2023-05-26 | 2023-10-10 | 华南理工大学 | Static voltage stability margin calculation method and device under N-1 fault of transmission and distribution system |
CN117067223B (en) * | 2023-10-16 | 2024-01-05 | 哈尔滨理工大学 | Six-foot robot free gait planning method based on motion stability estimation |
CN117215205B (en) * | 2023-11-09 | 2024-02-06 | 国网经济技术研究院有限公司 | DC system control parameter analysis method based on decision tree and ISS theory |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02280626A (en) * | 1989-04-21 | 1990-11-16 | Mitsubishi Electric Corp | Voltage and reactive power control for power system |
CN101247045A (en) * | 2008-03-20 | 2008-08-20 | 天津大学 | Electric voltage safety monitoring method based on voltage stabilization field in partition load space |
CN101841154A (en) * | 2010-04-07 | 2010-09-22 | 云南电力试验研究院(集团)有限公司 | Voltage stability margin real-time evaluation and optimum control method after grid major failure |
CN102403717A (en) * | 2011-11-18 | 2012-04-04 | 中国南方电网有限责任公司 | Method for evaluating severity of power system fault |
-
2014
- 2014-02-12 CN CN201410049147.XA patent/CN103795058B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02280626A (en) * | 1989-04-21 | 1990-11-16 | Mitsubishi Electric Corp | Voltage and reactive power control for power system |
CN101247045A (en) * | 2008-03-20 | 2008-08-20 | 天津大学 | Electric voltage safety monitoring method based on voltage stabilization field in partition load space |
CN101841154A (en) * | 2010-04-07 | 2010-09-22 | 云南电力试验研究院(集团)有限公司 | Voltage stability margin real-time evaluation and optimum control method after grid major failure |
CN102403717A (en) * | 2011-11-18 | 2012-04-04 | 中国南方电网有限责任公司 | Method for evaluating severity of power system fault |
Also Published As
Publication number | Publication date |
---|---|
CN103795058A (en) | 2014-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103795058B (en) | The air extract analysis of electric power system and system failure sort method | |
Zheng et al. | Regression tree for stability margin prediction using synchrophasor measurements | |
Haidar et al. | Transient stability evaluation of electrical power system using generalized regression neural networks | |
CN101635457B (en) | Electric network parameter estimation method based on parameter sensitivity of state estimation residual error | |
CN103972884B (en) | A kind of power system state estimation method | |
CN100438251C (en) | A large electric system vulnerable line identifying method | |
CN105976257A (en) | Power grid vulnerability evaluation method based on membership function fuzzy comprehensive evaluation method | |
CN110417011A (en) | A kind of online dynamic secure estimation method based on mutual information Yu iteration random forest | |
CN104917174B (en) | Static voltage stability judgment method for active power distribution network | |
CN106250640A (en) | A kind of layering Dynamic Equivalence being applicable to area power grid | |
CN102170131B (en) | Limit power calculating method of associated electricity transmission cross-section based on parallel calculating model | |
CN105762777A (en) | Pilot protection method containing multi-T-connection inverter interfaced distributed generation power distribution network | |
CN109066651A (en) | The calculation method of wind-powered electricity generation-load scenarios limit transmitted power | |
CN105303454A (en) | Vulnerability based power grid early-warning method | |
CN105243187A (en) | Transient stability measurement index based parameter influence degree analysis method | |
Cruz et al. | Algorithm development for power system contingency screening and ranking using voltage-reactive power performance index | |
CN102855382A (en) | Method for resolving three-phase short-circuit fault critical clearing time (CCT) of power system on line | |
CN105260521B (en) | A kind of load Fusion Modeling Method based on analytic hierarchy process (AHP) | |
Chouhan et al. | A literature review on optimal placement of PMU and voltage stability | |
CN104505827B (en) | Complicated electric power system closed loop control method based on response message | |
Augutis et al. | Application of Bayesian method for electrical power system transient stability assessment | |
CN103647284B (en) | A kind of voltage stabilization Forecasting Methodology of discontinuity surface problem when solving list | |
Seng et al. | Load flow analysis using improved Newton-Raphson method | |
Saunders et al. | Transient stability analysis using potential energy indices for determining critical generator sets | |
Althowibi et al. | On-line voltage collapse indicator for power 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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |