CN105787815A - Method and system for screening N-2 fault - Google Patents
Method and system for screening N-2 fault Download PDFInfo
- Publication number
- CN105787815A CN105787815A CN201610168704.9A CN201610168704A CN105787815A CN 105787815 A CN105787815 A CN 105787815A CN 201610168704 A CN201610168704 A CN 201610168704A CN 105787815 A CN105787815 A CN 105787815A
- Authority
- CN
- China
- Prior art keywords
- branch road
- residue
- cut
- overload
- contribution rate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012216 screening Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000009826 distribution Methods 0.000 claims description 61
- 238000012546 transfer Methods 0.000 claims description 43
- 239000011159 matrix material Substances 0.000 claims description 22
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 claims description 8
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 claims description 8
- 239000004744 fabric Substances 0.000 claims description 3
- 230000004888 barrier function Effects 0.000 claims description 2
- 238000012502 risk assessment Methods 0.000 abstract description 4
- 238000004880 explosion Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 240000004272 Eragrostis cilianensis Species 0.000 description 1
- 108010074506 Transfer Factor Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Economics (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Computing Systems (AREA)
Abstract
The invention relates to a method and a system for screening an N-2 fault. An active power limit value for each branch in the power grid, active power variations before and after breaking of each remaining branch for a target branch and an active power value before breaking of each remaining branch for the target branch are acquired firstly; according to the above three values, an overload contribution rate for the target branch by breaking of each remaining branch is acquired respectively; the remaining branches are grouped according to the overload contribution rates, two branches whose overload contribution rate sum is no smaller than an overload contribution rate standard value are screened according to the grouped condition, and an N-2 fault is obtained. By adopting the scheme of the invention, each possible combination of two branches does not need to be screened, grouped screening is carried out according to the grouped condition, the huge amount of computation which is difficult to calculate caused by combination explosion can be avoided, and the computation amount for power grid risk analysis is greatly reduced.
Description
Technical field
The present invention relates to power grid security risk field, particularly relate to screening technique and the system of N-2 fault.
Background technology
Along with new forms of energy and the introducing of electricity market, and the load sharply increased, the safe operation face of electrical network
Face bigger challenge.Owing to multiple failure has considerable influence to the safe operation of electrical network, relevant departments are
Through having formulated corresponding rules and regulation, it is desirable to when the more much higher heavy fault of electrical network probability of happening, should ensure that
, the most there is not the phenomenons such as branch road overload, cutting load in electrical network properly functioning.But, despite accordingly
Code, but it is to be resolved to still suffer from many problems on technology realizes.One of them important problem is
When multiple failure is analyzed, the multiple forecast failure enormous amount caused due to multiple shot array, formed and be difficult to process
Data computation burden, considerably limit the power grid security risk analysis of real-time online.Therefore, modern electric
System is badly in need of a multiple failure screening technique fast and efficiently, establishes base for online power grid risk analysis
Plinth.
At present, the screening of excessive risk fault is studied by existing many scholars.Generally speaking, substantially may be used
To be divided into two classes, a class is to concentrate from the suitable forecast accident of the probability that breaks down that to filter out consequence degree higher
Fault, two classes are to concentrate from the forecast accident that consequence degree is suitable to filter out the fault that probability of happening is higher.According to
Characterize the difference of index of damage sequence degree, first kind method can be divided into again dynamic security index method and quiet
State safety index method.Mostly these methods are to investigate each forecast accident respectively, then general by simplify
Rate computational methods or simply spend index, filter out fault.Its common problem is due to multiple shot array, treats
Investigation forecast accident enormous amount, the method for calculating probability simplified despite the use of or degree index, but its efficiency
It is still the lowest, it is impossible to meet the online evaluation of power grid security risk.
Summary of the invention
Based on this, it is necessary to excessive for data computation burden in existing Contingency screening method, inefficient
Problem, it is provided that the screening technique of a kind of N-2 fault and system.
The screening technique of a kind of N-2 fault, comprises the following steps:
Obtain the active power limiting value of Targeted Tributary in electrical network;
Obtain Targeted Tributary each active power variable quantity remained before and after branch road cut-offs respectively and target in electrical network
Branch road active power value before each residue branch road cut-offs respectively;
According to active power limiting value, each meritorious power variation and the active power value of each Targeted Tributary, point
Do not obtain each residue branch breaking overload contribution rate to Targeted Tributary;
According to default overload contribution rate is interval, each residue branch road is grouped, it is thus achieved that residue branch road group;
In the residue branch road group more than or equal to the first overload contribution rate boundary value preset, filter out mesh
The overload contribution rate sum of mark branch road is more than or equal to two first residues of the overload contribution rate standard value preset
Two first residue branch roads are cut-off simultaneously are considered as N-2 fault by branch road;
One article the is filtered out in the residue branch road group more than or equal to the second overload contribution rate boundary value preset
Two residue branch roads, screen in the residue branch road group less than or equal to the 3rd overload contribution rate boundary value preset
Go out one article of the 3rd residue branch road, make the second residue branch road that overload contribution rate and the 3rd residue of Targeted Tributary to be propped up
Targeted Tributary is transshipped contribution rate sum more than or equal to the overload contribution rate standard value preset by road;By second
Residue branch road and the 3rd residue branch road cut-off simultaneously and are considered as N-2 fault.
A kind of screening system of N-2 fault, including with lower unit:
First acquiring unit, for obtaining the active power limiting value of Targeted Tributary in electrical network;
Second acquisition unit, respectively remains having before and after branch road cut-offs respectively for obtaining Targeted Tributary in electrical network
Merit power variation and Targeted Tributary active power value before each residue branch road cut-offs respectively;
3rd acquiring unit, for propping up according to active power limiting value, each meritorious power variation and each target
The active power value on road, obtains each residue branch breaking overload contribution rate to Targeted Tributary respectively;
Packet screening unit, for being grouped each residue branch road according to the overload contribution rate preset is interval, it is thus achieved that
Residue branch road group;In the residue branch road group more than or equal to the first overload contribution rate boundary value preset, sieve
Select and Targeted Tributary is transshipped contribution rate sum more than or equal to two of overload contribution rate standard value preset
Residue branch road, two residue branch roads cut-off simultaneously and are considered as N-2 fault;More than or equal to the second mistake preset
Carry in the residue branch road group of contribution rate boundary value and filter out the first residue branch road, less than or equal to the preset
The residue branch road group of three overload contribution rate boundary values filters out the second residue branch road, makes the first residue branch road pair
The overload contribution rate of Targeted Tributary and the second residue branch road to the overload contribution rate sum of Targeted Tributary more than or etc.
In default overload contribution rate standard value;First residue branch road and the second residue branch road cut-off simultaneously and are considered as N-2
Fault.
According to the scheme of the invention described above, it is first to obtain the active power limiting value of each bar branch road in electrical network,
Targeted Tributary active power variable quantity before and after each residue branch road cut-offs respectively and Targeted Tributary prop up in each residue
Road cut-off respectively before active power value, obtain each residue branch breaking respectively to target according to these three data
The overload contribution rate of branch road;Further according to overload contribution rate, each residue branch road is grouped, according to packet situation
Filter out overload contribution rate sum and be not less than two branch roads of overload contribution rate standard value, obtain N-2 fault.
Due to the build-in attribute of electrical network, the overload contribution rate of major part branch road is the least, the only mistake of fraction branch road
Carry contribution rate bigger, therefore use the solution of the present invention to enter without two branch combinations possible to each
Row filter, but screen according to packet situation, thus what being difficult to of avoiding that multiple shot array causes calculated
Huge amount of calculation, is greatly reduced the amount of calculation that power grid risk is analyzed.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the screening technique of N-2 fault in one of them embodiment;
Fig. 2 is the overload contribution rate distribution map of a network system in one of them embodiment;
Fig. 3 is the overload contribution rate distribution map of a network system in one of them embodiment;
Fig. 4 is the structural representation of the screening system of N-2 fault in one of them embodiment.
Detailed description of the invention
For making the purpose of the present invention, technical scheme and advantage clearer, below in conjunction with accompanying drawing and enforcement
Example, is described in further detail the present invention.Should be appreciated that detailed description of the invention described herein
Only in order to explain the present invention, do not limit protection scope of the present invention.
Shown in Figure 1, for the embodiment of screening technique of the N-2 fault of the present invention.In this embodiment
The screening technique of N-2 fault, comprises the following steps:
Step S101: obtain the active power limiting value of Targeted Tributary in electrical network;
Step S102: obtain Targeted Tributary each active power remained before and after branch road cut-offs respectively in electrical network and become
Change amount and Targeted Tributary active power value before each residue branch road cut-offs respectively;
Step S103: meritorious according to active power limiting value, each meritorious power variation and each Targeted Tributary
Performance number, obtains each residue branch breaking overload contribution rate to Targeted Tributary respectively;
Step S104: according to default overload contribution rate is interval, each residue branch road is grouped, it is thus achieved that residue branch road
Group;In the residue branch road group more than or equal to the first overload contribution rate boundary value preset, filter out mesh
The overload contribution rate sum of mark branch road is more than or equal to two first residues of the overload contribution rate standard value preset
Two first residue branch roads are cut-off simultaneously are considered as N-2 fault by branch road;
Step S105: sieve in the residue branch road group more than or equal to the second overload contribution rate boundary value preset
Select one second residue branch road, propping up less than or equal to the residue of the 3rd overload contribution rate boundary value preset
Lu Zuzhong filters out one article of the 3rd residue branch road, make the second residue branch road to the overload contribution rate of Targeted Tributary and
Targeted Tributary is transshipped contribution rate sum more than or equal to the overload contribution rate standard preset by the 3rd residue branch road
Value;Second residue branch road and the 3rd residue branch road are cut-off simultaneously and be considered as N-2 fault.
In the present embodiment, first obtaining the active power limiting value of each bar branch road in electrical network, Targeted Tributary is respectively
Remain the active power variable quantity before and after branch road cut-offs respectively and Targeted Tributary before each residue branch road cut-offs respectively
Active power value, obtain the overload tribute to Targeted Tributary of each residue branch breaking respectively according to these three data
Offer rate;Further according to overload contribution rate, each residue branch road is grouped, filters out two according to packet situation
The binary combination on road, obtains N-2 fault.Due to the build-in attribute of electrical network, the overload contribution of major part branch road
Rate is the least, and the only overload contribution rate of fraction branch road is bigger, therefore use the solution of the present invention without
The binary combination that each is possible is screened, but screens according to packet situation, thus avoid
The huge amount of calculation that what multiple shot array caused be difficult to calculates, is greatly reduced the amount of calculation that power grid risk is analyzed.
Preferably, overload contribution rate interval, the first overload contribution rate boundary value, the second overload contribution rate boundary
Value, the 3rd overload contribution rate boundary value and overload contribution rate standard value are default, can be according to actual needs
It is adjusted.
In a specific embodiment, according toObtain each residue branch breaking to target
The overload contribution rate of branch road, in formula, Cα, βRepresent that branch road β cut-offs the overload contribution rate to branch road α, Δ fα,βTable
Show that branch road β cut-offs the active power variable quantity of before and after branch road α, fαRepresent that branch road β cut-offs the meritorious of front branch road α
Performance number, fαlimitRepresent the active power limiting value of branch road α.
When power system is properly functioning, branch breaking can cause system load flow to redistribute, consequently, it is possible to draw
Play some branch road overload.In order to another branch road overload is made the size contributed by a certain branch breaking of quantization signifying,
Overload contribution rate is used to characterize.The definition of overload contribution rate is that branch road β cut-offs before and after's branch road α active power
Variable quantity cut-off the ratio of active power margin of front branch road α with branch road β, active power margin refers to that branch road α's is meritorious
Power threshold and branch road β cut-off the difference of the active power of front branch road α.
Wherein in an embodiment, before and after obtaining Targeted Tributary each residue branch road cut-offfing respectively in electrical network
The step of active power variable quantity comprises the following steps:
Obtain each residue branch road respectively and Targeted Tributary is cut-off distribution factor;
Obtain the active power value before cut-offfing of each residue branch road respectively;
According to respectively cut-offfing distribution factor and each active power value before cut-offfing remaining branch road, obtain target
Branch road active power variable quantity before and after each residue branch road cut-offs respectively.
In the present embodiment, it is not necessary to special detection Targeted Tributary gaining merit before and after each residue branch road cut-offs respectively
Power variation, can cut-off distribution factor by each residue branch road and respectively remain branch road Targeted Tributary
The active power value before cut-offfing to calculate acquisition, convenient.
In a specific embodiment, according toObtain Targeted Tributary at each residue branch breaking
Active power variable quantity front and back, in formula, dα, βRepresent that branch road α is cut-off distribution factor, f by branch road ββTable
Show branch road β cut-off before active power value, Δ fα,βRepresent that branch road β cut-offs the active power change of before and after branch road α
Change amount.The sensitivity to other branch breaking of branch breaking distribution factor (LODF) Branch Power Flow, its energy
Change and characterize the impact on other Branch Power Flow of a certain branch breaking.
Wherein in an embodiment, obtain each residue branch road and Targeted Tributary is cut-off the step of distribution factor
Comprise the following steps:
Obtain the power transfer distribution factor of the Targeted Tributary when each residue branch road occurs power to shift respectively, with
And the power transfer distribution factor of self when respectively residue branch road occurs power to shift respectively;
Power transfer distribution factor according to each Targeted Tributary and the power transfer distribution factor of each residue branch road,
Obtain each residue branch road and Targeted Tributary is cut-off distribution factor.
In the present embodiment, Targeted Tributary is cut-off distribution factor by each residue branch road respectively by each residue branch road
When there is power transfer, the power transfer distribution factor of each Targeted Tributary and the power transfer point of each residue branch road
The cloth factor determines.In electrical network, it is easily obtained power transfer distribution factor, is also easy for obtaining each residue and props up
Targeted Tributary is cut-off distribution factor by road.
In a specific embodiment, according toObtain each residue branch road to Targeted Tributary
Cut-off distribution factor, in formula, dα, βRepresent that branch road α is cut-off distribution factor, p by branch road βα, βRepresent and propping up
There is the power transfer distribution factor of branch road α, p during power transfer in road ββ, βRepresenting occurs power to turn at branch road β
The power transfer distribution factor of self during shifting.Generally, pβ, β≠ 1, work as pβ, βWhen=1, dα, βNothing
Definition, now represents that electrical network splits into two systems because branch road β cut-offs.
In another specific embodiment, according toObtain the power of each Targeted Tributary
Transfer distribution factor;
In formula, pα, βRepresent the power transfer distribution factor of the branch road α when branch road β occurs power transfer, x α
Represent the reactance of branch road α, eαRepresenting row vector, its element is at the start node of branch road α and terminal node pair
That answers is classified as 1 and-1, and remaining is 0, and B represents the susceptance matrix of all branch roads,Represent column vector, its yuan
The plain behavior 1 and-1 corresponding with terminal node at the start node of branch road β, remaining is 0;
According toObtain the power transfer distribution factor of each residue branch road;
In formula, pβ, βRepresent the power transfer distribution factor of self, x when branch road β occurs power transferβRepresent
The reactance of branch road β, eβRepresenting row vector, its element is corresponding with terminal node at the start node of branch road β
Being classified as 1 and-1, remaining is 0.
Distribution factor is the electric system state quantity linear approximation to the sensitivity that node injecting power changes,
In the present invention, except branch breaking distribution factor, further relate to power transfer distribution factor.Power transfer distribution
The factor (PTDF) is the sensitivity that power between bus nodes is shifted by Branch Power Flow, and it can quantization signifying mother
Between line node power transfer amount on a certain branch road trend affect size.Its definition can use following mathematics public
Formula represents:
Wherein, alpha represents that branch road, the Latin alphabet i and j represent bus nodes, T(i,j)Represent from joint
Point i injects active power T, extracts power T, Δ f out from node jαRepresent that the active power on branch road α is at power
Transfer T(i,j)Variable quantity before and after generation.It addition, when node i is connected by branch road β with node j, pα,(i,j)Can
It is written as pα,β。
Power transfer factor can be calculated by equation below:
Wherein, xαRepresent the reactance of branch road α, eαBeing row vector, its element is at the start node of branch road α and end
What only node was corresponding is classified as 1 and-1, and remaining is 0, and B represents the susceptance matrix of all branch roads, e(i,j)Represent row
Vector, its element is respectively 1 and-1 at the i-th row and jth row, and remaining is 0.
In a preferred embodiment, for an electrical network with N bar branch road, its overload contribution rate can
To form the Matrix C of a N × N, its Elements C (α, β)=cα,β.Matrix C can obtain as follows.
(1) PTDF matrix is calculated
PTDF=X-1·φ·B-1·φT
(2) LODF matrix is calculated
LODF=PTDF (E diag (PTDF))-1
(3) branch power nargin matrix is calculated
MARGIN=LIMIT-PF
(4) computation overload contribution rate matrix
C=LODF.* (ONES PFT)./(MARGIN·ONEST)
Wherein, symbol " " representing matrix multiplication, the product of respective element in " .* " representing matrix, even X=Y.*Z,
Then (((i, j), the division of respective element in " ./" representing matrix, " diag " represents nondiagonal element X for i, j) × Z for i, j)=Y
Element zero setting obtains diagonal matrix, and PTDF is that power shifts distribution factor matrix, PTDF (α, β)=pα,β, X is
One diagonal matrix, its element is the reactance of branch road, and φ is node branch road incidence matrix, and B is all branch roads
Susceptance matrix, LODF is branch breaking distribution factor matrix, LODF (α, β)=dα,β, E is unit matrix,
LIMIT and PF is column vector, and its element represents branch power limiting value and branch road active power value respectively,
ONES be element be the column vector of 1.
As a example by two actual electric power saving systems, its overload contribution rate is calculated, and statistical analysis, knot
Fruit is as shown in Figure 2 and Figure 3.In figure, abscissa is subinterval numbering, and subinterval is by transshipping contribution rate
Interval (0,1) be divided into what 200 subintervals number consecutively 1 to 200 were formed, ordinate is each
The frequency of branch road in subinterval.In figure, scatter diagram is formed with above-mentioned horizontal stroke, ordinate, and curve is to scatterplot
The nonlinear fitting of figure.It will be seen that branch road overload contribution rate has long-tail characteristic distributions, i.e. major part
The overload contribution rate on road is the least, and only fraction is bigger.This feature is to be dilute by electrical network admittance matrix
Dredge what this electrical network build-in attribute of matrix caused.It practice, overload contribution rate is that long-tail is distributed this most just
Feature, result in the N-2 Contingency screening strategy that the present invention announced and has the highest efficiency.
Packet based on overload contribution rate, refers to for a certain bar branch road, according to other branch road mistake to this branch road
Carry the size of contribution rate, other branch road is grouped.Illustrate as a example by branch road i below.
Other branch road is C to the overload contribution rate of branch road ii=C (i :)=(ci,1,ci,2..., ci,k,…,ci,n-1), wherein
C (i :) represents the i-th all row of row of overload contribution rate matrix.According to following rule, other branch road is grouped:
O={k | ci,k>=1}, A={k | 0.8≤ci,k< 1}, B={k | 0.5≤ci,k< 0.8}, C={k | 0.1≤ci,k< 0.5},
D={k | 0.01≤ci,k< 0.1}, E={k | 0.005≤ci,k<0.01}.Group O is by transshipping contribution rate more than or equal to 1
Branch road forms, and according to the definition of overload contribution rate, in them, any one cut-offs and all can cause branch road i mistake
Carrying, these branch roads form serious N-1 fault.When power system meets N-1 check, group O is empty set,
Do not comprise any branch road.Branch road in group A to the overload contribution rate of branch road i between 0.8 and 1, other
The implication of group can be analogized.
For convenience, following two definition are given:
Definition 1: Class1 mistake, refers to that catastrophe failure being treated as in screening process is noncritical failure.
Definition 2: type 2 mistake, refers to that noncritical failure being treated as in screening process is catastrophe failure.
When certain two branch road cut-offs and causes power flow transfer simultaneously, if causing other branch road to transship, then this two
Bar cut-offs branch road one serious N-2 fault of composition;When certain two branch road cut-offs and causes power flow transfer simultaneously,
If not resulting in other branch road overload, then these two are cut-off branch road one non-serious N-2 fault of composition.Should refer to
Going out, in serious N-2 fault, any one branch breaking will not cause other branch road to transship.If in system
Be implicitly present in the N-1 that other branch road can be caused to transship, then these N-1 branch roads should individually take out and process,
N-2 is not formed, because any one the branch combinations in they and remaining all can cause other to pass by with other branch road
Carrying, the overload contribution rate of these N-1 branch roads is more than or equal to 1.
According to the definition of overload contribution rate, for branch road i, if ignore two cut-off between branch road influence each other,
And there are two branch road j and branch road k so that ci,j+ci,k>=1, wherein ci,j< 1, ci,k< 1, then and if only if props up
Road j and branch road k cut-offs simultaneously and branch road i can be caused to transship, and branch road j and branch road k forms a serious N-2 event
Barrier.But, ignore completely two cut-off between branch road influence each other and can increase the frequency of Class1 mistake,
And this is irrational, because the frequency reducing Class1 mistake should be the primary and foremost purpose of Contingency screening.Separately
On the one hand, due to two cut-off between branch road influence each other more complicated, it is difficult to by simple mode will
It quantifies.Therefore, considering the factor in terms of above-mentioned two, the present invention is to influencing each other between two branch roads
Give qualitative consideration, i.e. work as ci,j+ci,k>=0.85, wherein ci,j< 1, ci,k< when 1, just by branch road j and branch road
K cut-offs simultaneously and is considered as serious N-2 fault.
Based on above-mentioned thinking, the serious N-2 source of failure that branch road i transships can be caused in following several ways: (1)
The N-2 that in group ABC=A ∪ B ∪ C, any two branch roads are formed;Article (2) one, branch road is from organizing AB=A ∪
B, another branch road is from organizing D, the N-2 of formation;Article (3) one, branch road is from group A, branch road from
Group E, the N-2 of formation.
When specifically applying, can arrange overload contribution rate standard value is 0.85, the first overload contribution rate boundary value
It is 0.1, is being group ABC, sieve more than or equal to the residue branch road group of the first overload contribution rate boundary value preset
Select Targeted Tributary is transshipped contribution rate sum more than or equal to 0.85 two residue branch roads, two residues
Branch road cut-offs simultaneously and is considered as serious N-2 fault;Can also arrange the second overload contribution rate boundary value is 0.5, the
Three overload contribution rate boundary values are 0.1, and corresponding residue branch road group is group ABDE, according to Targeted Tributary
The overload contribution rate sum condition more than or equal to 0.85, is divided into two kinds of situations: one is a residue branch road
Screening from group AB, another residue branch road screens from group D, and two residue branch roads cut-off simultaneously and are considered as
Serious N-2 fault;Another kind is that a residue branch road screens from group A, and another residue branch road is from group E
Middle screening, two residue branch roads cut-off simultaneously and are considered as serious N-2 fault.
For each branch road, find the serious N-2 fault that this branch road can be caused to transship, and add serious N-2
Failure collection.Finally remove duplicate keys, obtain final serious N-2 failure collection.
In network system, not all multiple failure all can cause the more serious consequence of ratio, therefore,
All of anticipation multiple failure need not be carried out detailed analysis.On the other hand, in Contingency screening method
Under help, it can be carried out detailed analysis by calculating the fault being used for screening in resource set and formulates
Corresponding counter-measure, thus it is greatly reduced the amount of calculation that power grid risk is analyzed.
It is the inherent characteristics that long-tail is distributed this electrical network that the present invention takes full advantage of overload contribution rate between branch road, root
According to the size of overload contribution rate, respective branch is grouped, then filters out overload contribution rate sum the least
In two branch roads of overload contribution rate standard value, obtain serious N-2 fault.This method avoid and investigate respectively
Each forecast failure, has the highest screening efficiency, lays a good foundation for online power grid risk assessment.
According to the screening technique of above-mentioned N-2 fault, the present invention also provides for the screening system of a kind of N-2 fault,
Embodiment with regard to the screening system of the N-2 fault of the present invention is described in detail below.
Shown in Figure 4, for the embodiment of screening system of the N-2 fault of the present invention.In this embodiment
The screening system of N-2 fault includes the first acquiring unit 201, second acquisition unit 202, the 3rd acquiring unit
203, it is grouped screening unit 204;
First acquiring unit 201, for obtaining the active power limiting value of Targeted Tributary in electrical network;
Second acquisition unit 202, before and after obtaining Targeted Tributary each residue branch road cut-off respectively in electrical network
Active power variable quantity and Targeted Tributary active power value before each residue branch road cut-offs respectively;
3rd acquiring unit 203, for according to active power limiting value, each meritorious power variation and each target
The active power value of branch road, obtains each residue branch breaking overload contribution rate to Targeted Tributary respectively;
Packet screening unit 204, for being grouped each residue branch road according to the overload contribution rate preset is interval, obtains
Branch road group must be remained;In the residue branch road group more than or equal to the first overload contribution rate boundary value preset,
Filter out and Targeted Tributary is transshipped contribution rate sum more than or equal to the two of the overload contribution rate standard value preset
Two residue branch roads are cut-off simultaneously are considered as N-2 fault by bar residue branch road;More than or equal to the preset
The residue branch road group of two overload contribution rate boundary values filters out the first residue branch road, less than or equal to presetting
The 3rd overload contribution rate boundary value residue branch road group in filter out the second residue branch road, make the first residue prop up
Road Targeted Tributary is transshipped contribution rate and second residue branch road the overload contribution rate sum of Targeted Tributary is more than
Or equal to the overload contribution rate standard value preset;First residue branch road and the second residue branch road are cut-off simultaneously and regard
For N-2 fault.
Wherein in an embodiment, second acquisition unit 202 obtains each residue branch road respectively to Targeted Tributary
Cut-off distribution factor, do not obtain the active power value before cut-offfing of each residue branch road, according to respectively cut-offfing point
Active power value before the cloth factor and each residue branch breaking, obtains Targeted Tributary at each residue branch road respectively
Active power variable quantity before and after cut-offfing.
Wherein in an embodiment, second acquisition unit 202 obtains, at each residue branch road, power occurs respectively
The power transfer distribution factor of Targeted Tributary during transfer, and when each residue branch road occurs power to shift respectively certainly
The power transfer distribution factor of body, shifts distribution factor and each residue branch road according to the power of each Targeted Tributary
Power transfer distribution factor, obtains each residue branch road and Targeted Tributary is cut-off distribution factor.
The screening system of the N-2 fault of the present invention and the screening technique one_to_one corresponding of the N-2 fault of the present invention,
Above-mentioned N-2 fault screening technique embodiment illustrate technical characteristic and beneficial effect be all applicable to
In the embodiment of the screening system of N-2 fault.
Each technical characteristic of embodiment described above can combine arbitrarily, for making description succinct, the most right
The all possible combination of each technical characteristic in above-described embodiment is all described, but, if these skills
There is not contradiction in the combination of art feature, is all considered to be the scope that this specification is recorded.
Embodiment described above only have expressed the several embodiments of the present invention, and it describes more concrete and detailed,
But can not therefore be construed as limiting the scope of the patent.It should be pointed out that, for this area
For those of ordinary skill, without departing from the inventive concept of the premise, it is also possible to make some deformation and change
Entering, these broadly fall into protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be with appended power
Profit requires to be as the criterion.
Claims (10)
1. the screening technique of a N-2 fault, it is characterised in that comprise the following steps:
Obtain the active power limiting value of Targeted Tributary in electrical network;
Obtain each active power change remained before and after branch road cut-offs respectively in described electrical network of described Targeted Tributary
Amount and described Targeted Tributary active power value before each described residue branch road cut-offs respectively;
According to having of described active power limiting value, each described active power variable quantity and each described Targeted Tributary
Merit performance number, obtains each described residue branch breaking overload contribution rate to described Targeted Tributary respectively;
According to default overload contribution rate is interval, each described residue branch road is grouped, it is thus achieved that residue branch road group;
In the residue branch road group more than or equal to the first overload contribution rate boundary value preset, filter out institute
The overload contribution rate sum stating Targeted Tributary is more than or equal to two first of the overload contribution rate standard value preset
Described two described first residue branch roads are cut-off simultaneously are considered as N-2 fault by residue branch road;
One article the is filtered out in the residue branch road group more than or equal to the second overload contribution rate boundary value preset
Two residue branch roads, screen in the residue branch road group less than or equal to the 3rd overload contribution rate boundary value preset
Go out one article of the 3rd residue branch road, make described second residue branch road to the overload contribution rate of described Targeted Tributary and institute
State the 3rd residue branch road and described Targeted Tributary is transshipped contribution rate sum more than or equal to described default overload
Contribution rate standard value;Described second residue branch road and described 3rd residue branch road are cut-off simultaneously and be considered as N-2 event
Barrier.
The screening technique of N-2 fault the most according to claim 1, it is characterised in that obtain described mesh
The step of mark branch road active power variable quantity before and after each residue branch road cut-offs respectively in described electrical network include with
Lower step:
Obtain each described residue branch road respectively and described Targeted Tributary is cut-off distribution factor;
Obtain the active power value before cut-offfing of each described residue branch road respectively;
Distribution factor and the active power value before cut-offfing of each described residue branch road is cut-off according to each,
Obtain described Targeted Tributary active power variable quantity before and after each described residue branch road cut-offs respectively.
The screening technique of N-2 fault the most according to claim 2, it is characterised in that obtain each described
The step that described Targeted Tributary is cut-off distribution factor by residue branch road comprises the following steps:
Obtain the power transfer distribution that power described Targeted Tributary when shifting occurs respectively at each described residue branch road
The factor, and the power transfer distribution factor of self when each described residue branch road occurs power to shift respectively;
Power transfer distribution factor according to each described Targeted Tributary and the power transfer point of each described residue branch road
The cloth factor, obtains each described residue branch road and described Targeted Tributary is cut-off distribution factor.
The screening technique of N-2 fault the most according to claim 3, it is characterised in that:
According toObtain the power transfer distribution factor of each described Targeted Tributary;
In formula, pα, βRepresent the power transfer distribution factor of the branch road α when branch road β occurs power transfer, xα
Represent the reactance of branch road α, eαRepresenting row vector, its element is at the start node of branch road α and terminal node pair
That answers is classified as 1 and-1, and remaining is 0, and B represents the susceptance matrix of all branch roads,Represent column vector, its yuan
The plain behavior 1 and-1 corresponding with terminal node at the start node of branch road β, remaining is 0;
According toObtain the power transfer distribution factor of each described residue branch road;
In formula, pβ, βRepresent the power transfer distribution factor of self, x when branch road β occurs power transferβRepresent
The reactance of branch road β, eβRepresenting row vector, its element is corresponding with terminal node at the start node of branch road β
Being classified as 1 and-1, remaining is 0.
The screening technique of N-2 fault the most according to claim 3, it is characterised in that:
According toObtain each described residue branch road and described Targeted Tributary cut-off distribution factor,
In formula, dα, βRepresent that branch road α is cut-off distribution factor, p by branch road βα, βRepresenting occurs power to turn at branch road β
The power transfer distribution factor of branch road α, p during shiftingβ, βRepresent the power of self when branch road β occurs power transfer
Transfer distribution factor.
The screening technique of N-2 fault the most according to claim 2, it is characterised in that:
According toObtain described Targeted Tributary active power before and after each described residue branch breaking
Variable quantity, in formula, dα, βRepresent that branch road α is cut-off distribution factor, f by branch road ββBefore representing that branch road β cut-offs
Active power value, Δ fα,βRepresent that branch road β cut-offs the active power variable quantity of before and after branch road α.
The screening technique of N-2 fault the most according to claim 1, it is characterised in that:
According toObtain each described residue branch breaking the overload of described Targeted Tributary is contributed
Rate, in formula, cα, βRepresent that branch road β cut-offs the overload contribution rate to branch road α, Δ fα,βBefore representing that branch road β cut-offs
The active power variable quantity of rear branch road α, fαRepresent that branch road β cut-offs the active power value of front branch road α, fαlimit
Represent the active power limiting value of branch road α.
8. the screening system of a N-2 fault, it is characterised in that include with lower unit:
First acquiring unit, for obtaining the active power limiting value of Targeted Tributary in electrical network;
Second acquisition unit, is used for obtaining described Targeted Tributary each residue branch road in described electrical network and cut-offs respectively
Active power variable quantity front and back and described Targeted Tributary wattful power before each described residue branch road cut-offs respectively
Rate value;
3rd acquiring unit, for according to described active power limiting value, each described active power variable quantity and
The active power value of each described Targeted Tributary, obtains each described residue branch breaking respectively to described Targeted Tributary
Overload contribution rate;
Packet screening unit, for each described residue branch road being grouped according to the overload contribution rate preset is interval,
Obtain residue branch road group;In the residue branch road group more than or equal to the first overload contribution rate boundary value preset,
Filter out and described Targeted Tributary is transshipped contribution rate sum more than or equal to the overload contribution rate standard value preset
Two described residue branch roads, described two described residue branch roads are cut-off simultaneously and are considered as N-2 fault;Greatly
In or equal to preset second overload contribution rate boundary value residue branch road group in filter out the first residue branch road,
In the residue branch road group less than or equal to the 3rd overload contribution rate boundary value preset, filter out the second residue prop up
Road, makes described first residue branch road that the overload contribution rate and described second of described Targeted Tributary is remained branch road pair
The overload contribution rate sum of described Targeted Tributary is more than or equal to described default overload contribution rate standard value;Will
Described first residue branch road and described second residue branch road cut-off simultaneously and are considered as N-2 fault.
The screening system of N-2 fault the most according to claim 8, it is characterised in that described second obtains
Take unit to obtain each described residue branch road respectively described Targeted Tributary is cut-off distribution factor, do not obtain each institute
State the active power value before cut-offfing of residue branch road, according to cut-offfing distribution factor and each described surplus described in each
Active power value before remaining branch breaking, obtains described Targeted Tributary before each described residue branch road cut-offs respectively
After active power variable quantity.
The screening system of N-2 fault the most according to claim 9, it is characterised in that described second
Acquiring unit obtains and occurs the power of power described Targeted Tributary when shifting to shift respectively at each described residue branch road
Distribution factor, and the power transfer distribution factor of self when each described residue branch road occurs power to shift respectively,
According to each described Targeted Tributary power transfer distribution factor and each described residue branch road power transfer distribution because of
Son, obtains each described residue branch road and described Targeted Tributary is cut-off distribution factor.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610168704.9A CN105787815B (en) | 2016-03-22 | 2016-03-22 | N-2 fault screening method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610168704.9A CN105787815B (en) | 2016-03-22 | 2016-03-22 | N-2 fault screening method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105787815A true CN105787815A (en) | 2016-07-20 |
CN105787815B CN105787815B (en) | 2020-01-07 |
Family
ID=56391688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610168704.9A Active CN105787815B (en) | 2016-03-22 | 2016-03-22 | N-2 fault screening method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105787815B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106451463A (en) * | 2016-11-09 | 2017-02-22 | 广东电网有限责任公司珠海供电局 | A method and apparatus for analyzing a thermal-stable section in a ring network system |
CN109167356A (en) * | 2018-10-16 | 2019-01-08 | 云南电网有限责任公司 | A kind of N-2 Contingency screening method |
CN109541361A (en) * | 2018-12-28 | 2019-03-29 | 深圳供电局有限公司 | A kind of public affairs low pressure side branch's load heavy duty prior-warning device and method |
CN110794263A (en) * | 2019-12-17 | 2020-02-14 | 国网新疆电力有限公司昌吉供电公司 | Method for positioning fault section of power distribution network line with distributed power supply |
CN112688314A (en) * | 2020-12-16 | 2021-04-20 | 国网河南省电力公司漯河供电公司 | Method for calculating static safety analysis of N-1 of power system based on PMU (phasor measurement Unit) measured data |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013812A (en) * | 2007-02-07 | 2007-08-08 | 浙江大学 | Method for configuring PMU taking regard of high risk cascading failure path |
-
2016
- 2016-03-22 CN CN201610168704.9A patent/CN105787815B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013812A (en) * | 2007-02-07 | 2007-08-08 | 浙江大学 | Method for configuring PMU taking regard of high risk cascading failure path |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106451463A (en) * | 2016-11-09 | 2017-02-22 | 广东电网有限责任公司珠海供电局 | A method and apparatus for analyzing a thermal-stable section in a ring network system |
CN106451463B (en) * | 2016-11-09 | 2019-03-29 | 广东电网有限责任公司珠海供电局 | The method and apparatus of the steady section of heat in a kind of analysis loop network system |
CN109167356A (en) * | 2018-10-16 | 2019-01-08 | 云南电网有限责任公司 | A kind of N-2 Contingency screening method |
CN109167356B (en) * | 2018-10-16 | 2022-03-18 | 云南电网有限责任公司 | N-2 fault screening method |
CN109541361A (en) * | 2018-12-28 | 2019-03-29 | 深圳供电局有限公司 | A kind of public affairs low pressure side branch's load heavy duty prior-warning device and method |
CN110794263A (en) * | 2019-12-17 | 2020-02-14 | 国网新疆电力有限公司昌吉供电公司 | Method for positioning fault section of power distribution network line with distributed power supply |
CN112688314A (en) * | 2020-12-16 | 2021-04-20 | 国网河南省电力公司漯河供电公司 | Method for calculating static safety analysis of N-1 of power system based on PMU (phasor measurement Unit) measured data |
CN112688314B (en) * | 2020-12-16 | 2023-10-27 | 国网河南省电力公司漯河供电公司 | PMU-based method for calculating N-1 static security analysis |
Also Published As
Publication number | Publication date |
---|---|
CN105787815B (en) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105787815A (en) | Method and system for screening N-2 fault | |
Dabachi et al. | Energy consumption, energy price, energy intensity environmental degradation, and economic growth nexus in African OPEC countries: evidence from simultaneous equations models | |
CN101291062B (en) | Analysis method of DC power flow in dynamic topology electric network | |
CN103985066A (en) | Method for evaluating power system static risk based on hybrid power flow | |
CN107871216A (en) | A kind of recognition methods of power distribution network fragility node | |
CN103985065A (en) | Method for evaluating electric power system risk based on fault pre-scanning | |
CN105528645A (en) | Frangibility prediction method for large power grid | |
CN104992266A (en) | Method of determining power grid node importance degree and system thereof | |
CN105354675A (en) | Key transmission section identification-based cascading failure analysis method for alternating current/direct current power network | |
CN104092210A (en) | Online static security analysis method with N-2 opening rapid scanning function | |
US20220083655A1 (en) | Chip security analysis method based on petri net | |
CN105205571A (en) | Risk-considered urban power network operation security assessment method | |
CN112329376B (en) | Monte Carlo simulation-based substation system shock resistance toughness quantitative evaluation method | |
CN107947151A (en) | A kind of recognition methods of electric system key node | |
CN103823837A (en) | Key transmission section searching method based on fuzzy clustering and shortest path | |
CN106780134A (en) | The determination method and system of the wind power generating algorithm apoplexy electricity condition number upper limit | |
CN102570450A (en) | Static reliability assessment method for complex power grid | |
CN106655201A (en) | Security domain-based safe optimization and control method for electric power thermal stability | |
CN108182511A (en) | It is a kind of based on Demand Side Response reserve value assessment method of the sum of ranks than method | |
CN106684917A (en) | Positional dynamic optimization and allowed capacity calculating method and system for distributed power supplies | |
CN103400213A (en) | Backbone net rack survivability assessment method based on LDA (Linear Discriminant Analysis) and PCA (Principal Component Analysis) | |
CN105024645A (en) | Matrix evolution-based photovoltaic array fault location method | |
CN104950205A (en) | Quantitative analysis method for effect of AC (alternating current) system failure on DC (direct current) line | |
CN104867063A (en) | Power grid multiple major failure recognition method based on a linear weighting method | |
CN105427200A (en) | Power plant fault diagnosis method based on fuzzy theory and data fusion |
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 |