CN106059838B  A kind of reliability of relay protection calculation method and device  Google Patents
A kind of reliability of relay protection calculation method and device Download PDFInfo
 Publication number
 CN106059838B CN106059838B CN201610617812.XA CN201610617812A CN106059838B CN 106059838 B CN106059838 B CN 106059838B CN 201610617812 A CN201610617812 A CN 201610617812A CN 106059838 B CN106059838 B CN 106059838B
 Authority
 CN
 China
 Prior art keywords
 state
 current search
 network state
 node
 space tree
 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.)
 Expired  Fee Related
Links
Classifications

 H—ELECTRICITY
 H04—ELECTRIC COMMUNICATION TECHNIQUE
 H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
 H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
 H04L41/14—Network analysis or design
 H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network

 H—ELECTRICITY
 H04—ELECTRIC COMMUNICATION TECHNIQUE
 H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
 H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
 H04L41/08—Configuration management of networks or network elements
 H04L41/0803—Configuration setting
 H04L41/0823—Configuration setting characterised by the purposes of a change of settings, e.g. optimising configuration for enhancing reliability
 H04L41/0836—Configuration setting characterised by the purposes of a change of settings, e.g. optimising configuration for enhancing reliability to enhance reliability, e.g. reduce downtime
Abstract
The present invention discloses a kind of reliability of relay protection calculation method and device, obtains the Making by Probability Sets of multimodal capacityconstrained matrix and each default element manipulation in kind under each mode according to the historical traffic of default element in kind each in intelligent substation；Building node can keep to the side failure stochasticflow networks model；Corresponding effective status set is obtained based on stochasticflow networks model；According to effective status set, the reliability of relay protection under stochasticflow networks model is calculated.Based on method disclosed above, using stochasticflow networks model to the Calculation of Reliability of relay protection, to achieve the purpose that carry out the reliability of complicated electric power system quantitative calculating.
Description
Technical field
The present invention relates to reliability assessment technical field more particularly to a kind of reliability of relay protection calculation method and dresses
It sets.
Background technique
Protective relaying device as one of the critical elements in secondary equipment in power system, can and alarm by direct shadow
Ring the safety and reliability of Operation of Electric Systems.Longtime running practice have shown that, protective relaying device hardware and software itself lose
The probability very little of effect and builtin selfdiagnostic function often can detecte this kind of failure event, therefore believe under smart grid background
Breath stream becomes the key factor for restricting protective relaying device reliability.
Currently, existing reliability of relay protection appraisal procedure is based primarily upon instantaneous state probability, half markoff process
Or monte carlo simulation methodology carries out quantitative analysis to relay protection and carried out using reliability block diagram and minimal path sets algorithm can
It is calculated by property.But the appraisal procedure only considered connectivity of the information flow from source point to meeting point, the operating status of acquisition is only
Entirely ineffective and two kinds of normal operation, in view of the factor of influence element manipulation, therefore potentially network failure and information
Reliability decrease caused by current mass declines does not embody.
In view of this, existing technical solution, when calculating reliability of relay protection, there are asking for calculated result inaccuracy
Topic.
Summary of the invention
In view of this, the present invention provides a kind of reliability of relay protection calculation method and device, to solve existing technology
The problem for the reliability of relay protection result inaccuracy that scheme calculates.Technical solution is as follows:
A kind of reliability of relay protection calculation method is applied to reliability of relay protection computing device, comprising:
According to the multimodal capacityconstrained matrix of historical traffic acquisition of default element in kind each in intelligent substation and respectively
Making by Probability Sets of a default element manipulation in kind under each mode；
Building node can keep to the side failure stochasticflow networks model G=(V, E, C, P)；
Wherein, node V is the physical connection cluster tool of each default element in kind, and side E is the default material object
The set of element, C are the multimodal capacityconstrained matrix, and P is each default element manipulation in kind under each mode
Making by Probability Sets；
Corresponding effective status set W is obtained based on the stochasticflow networks model G；
According to the effective status set W, the reliability of relay protection under the stochasticflow networks model G is calculatedWherein, l is lower boundary point Y in the effective status set W_{d}Number.
It is preferably, described that corresponding effective status set W is obtained based on the stochasticflow networks model G, comprising:
Each network state is obtained based on the stochasticflow networks model G, and generates state space tree, wherein Ge Gesuo
State the node that network state is the state space tree；
Zero network state is generated based on the network state, and using zero network state as the state space tree
First father node；
Determine a default element in kind of network state described in the state space tree as current search member
Part；
Determine that next network state is current search node in the state space tree, and according to default minimal path sets meter
Calculation method obtains the information source node of the state space tree and the minimal path sets of information meeting point；
Judge whether the information source node and the network topology of the information meeting point are connected to according to the minimal path sets；
If be not connected to, back tracking operation is executed, and returns and executes next network in the determination state space tree
State is current search node, the step for；
If connection, maximum stream calculation is started according to default maximumflow algorithm, obtains the network of the current search node
State capacity simultaneously judges whether the network state capacity is less than preset need capacity；
If the network state capacity is not less than preset need capacity, the network state capacity is determined as the lower bound
Point Y_{d}It is stored in state set L, executes the back tracking operation, and return next in the execution determination state space tree
A network state is current search node, the step for；
If the network state capacity be less than preset need capacity, judge the current search node mode serial number whether
Less than preset mode serial number maximum value；
If the mode serial number of the current search node is less than preset mode serial number maximum value, returns and execute the determining institute
Stating next network state in state space tree is current search node, the step for；
If the mode serial number of the current search node is not less than preset mode serial number maximum value, the backtracking behaviour is executed
Make, and returning to next network state in the execution determination state space tree is current search node, the step for；
Wherein, the back tracking operation includes:
Judge the current search element whether be the state space tree final search element；
If the current search element is not the final search element of the state space tree, by the state space tree
In a upper network state be the current search node, and determine next default member in kind in the state space tree
Part is as the current search element；
If the current search element is the final search element of the state space tree, will be in the state space tree
A upper network state is the current search node, and judges whether the current search node is zero network state；
If the current search node is not zero network state, determine next described pre in the state space tree
If element in kind is as the current search element；
If the current search node is zero network state, the state set L is determined as the effective status
Set W, and terminate maximum stream calculation.
Preferably, the value range of Probability p of each default element manipulation in kind under each mode is 0~1.
Preferably, the default minimal path sets calculation method includes: connection matrix method.
Preferably, the default maximumflow algorithm includes: extensions path FordFulkerson algorithm.
A kind of reliability of relay protection computing device, comprising: obtain module, model construction module, effective status set and obtain
Modulus block and Calculation of Reliability module；
The acquisition module, it is multimodal for being obtained according to the historical traffic of default element in kind each in intelligent substation
The Making by Probability Sets of capacityconstrained matrix and each default element manipulation in kind under each mode；
The model construction module, for construct node can keep to the side failure stochasticflow networks model G=(V, E, C, P)；
Wherein, node V is the physical connection cluster tool of each default element in kind, and side E is the collection of the default element in kind
It closes, C is the multimodal capacityconstrained matrix, and P is Making by Probability Sets of each default element manipulation in kind under each mode；
The effective status set obtains module, for obtaining corresponding effective shape based on the stochasticflow networks model G
State set W；
The Calculation of Reliability module, for calculating the stochasticflow networks model G according to the effective status set W
Under reliability of relay protectionWherein, l is lower boundary point Y in the effective status set W_{d}
Number.
Preferably, it includes: state space tree generation unit, the generation of zero network state that the effective status set, which obtains module,
Unit searches for element determination unit, determines computing unit, the first judging unit, the first execution unit, backtracking execution unit, the
Two judging units, the second execution unit, third judging unit and third execution unit；
The state space tree generation unit, for obtaining each network state based on the stochasticflow networks model G, and
Generate state space tree, wherein each network state is the node of the state space tree；
The zero network state generation unit, for generating zero network state based on the network state, and by described zero
First father node of the network state as the state space tree；
Described search element determination unit, one for determining network state described in the state space tree are described pre
If element in kind is as current search element；
The determining computing unit, for determining, next network state is current search section in the state space tree
Point, and the information source node of the state space tree and the minimal path of information meeting point are obtained according to default minimal path sets calculation method
Collection；
First judging unit, for judging the information source node and the information meeting point according to the minimal path sets
Network topology whether be connected to；If so, triggering the second judgment unit；If not, triggering first execution unit；
First execution unit for triggering backtracking execution unit, and enters the determining computing unit；
The second judgment unit is used for according to the maximum stream calculation of maximumflow algorithm starting is preset, and acquisition is described currently to search
The network state capacity of socket point simultaneously judges whether the network state capacity is less than preset need capacity；If it is not, described in triggering
Second execution unit；If so, triggering the third judging unit；
Second execution unit, for the network state capacity to be determined as the lower boundary point Y_{d}It is stored in state set
It closes in L, triggers the backtracking execution unit, and trigger the determining computing unit；
The third judging unit, for judging whether the mode serial number of the current search node is less than preset mode sequence
Number maximum value；It is to trigger the determining computing unit；It is no, trigger the third execution unit；
The third execution unit for triggering the backtracking execution unit, and enters the determining computing unit；
The backtracking execution unit, comprising: the 4th judging unit, the 4th execution unit, the 5th judging unit, the 5th execute
Unit and the 6th execution unit；
4th judging unit, for judging whether the current search element is finally searching for the state space tree
Rope element；If not, triggering the 4th execution unit；If so, triggering the 5th judging unit；
4th execution unit, for being the current search section by a network state upper in the state space tree
Point, and determine that next default element in kind is as the current search element in the state space tree；
5th judging unit, for being the current search section by a network state upper in the state space tree
Point, and judge whether the current search node is zero network state；If it is not, triggering the 5th execution unit；If so,
Trigger the 6th execution unit；
5th execution unit, for determining that next default element in kind is as institute in the state space tree
State current search element；
6th execution unit for the state set L to be determined as the effective status set W, and terminates most
Big stream calculation.
It compares and the prior art, what the present invention realized has the beneficial effect that
A kind of reliability of relay protection calculation method provided by the present invention and device above, according to each in intelligent substation
The historical traffic of a default element in kind obtains multimodal capacityconstrained matrix and each default element manipulation in kind in each mode
Under Making by Probability Sets；Building node can keep to the side failure stochasticflow networks model；It is obtained based on stochasticflow networks model corresponding
Effective status set；According to effective status set, the reliability of relay protection under stochasticflow networks model is calculated.Based on abovementioned public affairs
The method opened, using stochasticflow networks model to the Calculation of Reliability of relay protection, to reach reliable to complicated electric power system
Property carry out the purpose of quantitative calculating.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The embodiment of invention for those of ordinary skill in the art without creative efforts, can also basis
The attached drawing of offer obtains other attached drawings.
Fig. 1 is a kind of reliability of relay protection calculation method flow chart disclosed in the embodiment of the present invention one；
Fig. 2 is a kind of reliability of relay protection calculation method partial process view disclosed in the embodiment of the present invention two；
Fig. 3 is a kind of reliability of relay protection calculation method partial process view disclosed in the embodiment of the present invention two；
Fig. 4 is a kind of reliability of relay protection computing device structure schematic diagram disclosed in the embodiment of the present invention three；
Fig. 5 is a kind of reliability of relay protection computing device partial structure diagram disclosed in the embodiment of the present invention four；
Fig. 6 is a kind of reliability of relay protection computing device partial structure diagram disclosed in the embodiment of the present invention four.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Embodiment one
A kind of reliability of relay protection calculation method disclosed in the embodiment of the present invention one is applied to reliability of relay protection meter
Device is calculated, flow chart is as shown in Figure 1, reliability of relay protection calculation method includes:
S101 obtains multimodal capacityconstrained matrix according to the historical traffic of default element in kind each in intelligent substation
With Making by Probability Sets of each default element manipulation in kind under each mode；
S102, building node can keep to the side failure stochasticflow networks model G=(V, E, C, P)；
Wherein, node V is the physical connection cluster tool of each default element in kind, and side E is default element in kind
Set, C is multimodal capacityconstrained matrix, and P is Making by Probability Sets of each default element manipulation in kind under each mode；
During executing step S102, during building node can keep to the side the stochasticflow networks model of failure,
The physical connection cluster tool V={ v of each default element in kind_{i} 1≤i≤n } node as stochasticflow networks model, in advance
If the set E={ e of element in kind_{i} 1≤i≤a } side as stochasticflow networks model, C={ c_{ij}1≤i≤n+a,0≤j≤
h_{i}It is used as multimodal capacityconstrained matrix, P={ p_{ij}1≤i≤n+a,0≤j≤h_{i}It is that each default element manipulation in kind exists
Making by Probability Sets under each mode, wherein h_{i}For the corresponding mode serial number of ith of element in element；
S103 obtains corresponding effective status set W based on stochasticflow networks model G；
S104 calculates the reliability of relay protection under stochasticflow networks model G according to effective status set WWherein, l is lower boundary point Y in effective status set W_{d}Number.
It should be noted that the value range of Probability p of each default element manipulation in kind under each mode is 0~1, institute
There is Probability p of the default in kind element manipulation under each mode and is 1.
Reliability of relay protection calculation method disclosed by the embodiments of the present invention, according to default material object each in intelligent substation
The historical traffic of element obtains the probability set of multimodal capacityconstrained matrix and each default element manipulation in kind under each mode
It closes；Building node can keep to the side failure stochasticflow networks model；Corresponding effective status collection is obtained based on stochasticflow networks model
It closes；According to effective status set, the reliability of relay protection under stochasticflow networks model is calculated.Based on method disclosed above,
Using stochasticflow networks model to the Calculation of Reliability of relay protection, complicated electric power system reliability is quantified to reach
The purpose of calculating.
Embodiment two
Learning method is indicated based on a kind of multimodal data disclosed in the embodiments of the present invention one, as illustrated in FIG. 1
In step S103, the specific implementation procedure of corresponding effective status set W is obtained based on stochasticflow networks model G, such as Fig. 2 institute
Show, comprising:
S201 obtains each network state based on stochasticflow networks model G, and generates state space tree, wherein each net
Network state is the node of the state space tree；
During executing step S201, corresponding state space Ω=(Y is obtained based on stochasticflow networks model G_{1},
Y_{2},…,Y_{i},Y_{m}), and state space tree is generated according to each network state, wherein each network state Y_{i}It is empty for the state
Between the node set.
S202 generates zero network state based on network state, and using zero network state as the first father of state space tree
Node；
During executing step S202, based on obtained each network state Y_{i}Generate zero network state, and by zero
First father node, that is, initiating searches point of the network state as state space tree.
S203 determines a default element in kind of network state described in state space tree as current search element；
S204 determines that next network state is current search node in state space tree, and according to default minimal path sets
Calculation method obtains the information source node of state space tree and the minimal path sets of information meeting point；
S205 judges whether information source node and the network topology of the information meeting point are connected to according to minimal path sets；
S206 executes back tracking operation if be not connected to, and returns to execution and determine next network state in state space tree
For current search node, the step for；
S207 starts maximum stream calculation according to default maximumflow algorithm, obtains the network of current search node if connection
State capacity simultaneously judges whether network state capacity is less than preset need capacity；
Network state capacity is determined as lower boundary point Y if network state capacity is not less than preset need capacity by S208_{d}It deposits
It is stored in state set L, executes back tracking operation, and return to execution and determine that next network state is currently to search in state space tree
Socket point, the step for；
S209 judges whether the mode serial number of current search node is small if network state capacity is less than preset need capacity
In preset mode serial number maximum value；
S210 is returned described in executing if the mode serial number of the current search node is less than preset mode serial number maximum value
Determine that next network state is current search node in the state space tree, the step for；
S211 executes back tracking operation if the mode serial number of current search node is not less than preset mode serial number maximum value, and
It returns to execute and determines that next network state is current search node in the state space tree, the step for；
Wherein, the specific implementation procedure of back tracking operation is as shown in Figure 3, comprising:
S301, judge current search element whether be state space tree final search element；
S302, if current search element is not the final search element of state space tree, by state space tree upper one
A network state is current search node, and determines that next default element in kind is as current search member in state space tree
Part；
S303, if current search element is the final search element of state space tree, by upper one in state space tree
Network state is current search node, and judges whether current search node is zero network state；
S304 determines next default element in kind in state space tree if current search node is not zero network state
As current search element；
State set L is determined as effective status set W, and terminate if current search node is zero network state by S305
Maximum stream calculation.
It should be noted that default minimal path sets calculation method includes but is not limited to connection matrix method, it can be according to reality
It needs to select suitable minimal path sets calculation method；Default maximumflow algorithm includes but is not limited to extensions path Ford
Fulkerson algorithm, the method that suitable calculating maxflow can be selected according to actual needs.
Reliability of relay protection calculation method disclosed by the embodiments of the present invention, according to default material object each in intelligent substation
The historical traffic of element obtains the probability set of multimodal capacityconstrained matrix and each default element manipulation in kind under each mode
It closes；Building node can keep to the side failure stochasticflow networks model；Corresponding effective status collection is obtained based on stochasticflow networks model
It closes；According to effective status set, the reliability of relay protection under stochasticflow networks model is calculated.Based on method disclosed above,
Using stochasticflow networks model to the Calculation of Reliability of relay protection, complicated electric power system reliability is quantified to reach
The purpose of calculating.
Embodiment three
Based on the reliability of relay protection calculation method that each embodiment of aforementioned present invention provides, the present embodiment three is then corresponding public
The reliability of relay protection computing device of the abovementioned reliability of relay protection calculation method of execution, structural schematic diagram such as Fig. 4 are opened
Shown, reliability of relay protection computing device includes 400: obtaining module 401, model construction module 402, effective status set and obtains
Modulus block 403 and Calculation of Reliability module 404；Wherein,
Module 401 is obtained, it is multimodal for being obtained according to the historical traffic of default element in kind each in intelligent substation
The Making by Probability Sets of capacityconstrained matrix and each default element manipulation in kind under each mode；
Model construction module 402, for construct node can keep to the side failure stochasticflow networks model G=(V, E, C, P)；Its
In, node V is the physical connection cluster tool of each default element in kind, and side E is the set of the default element in kind,
C is the multimodal capacityconstrained matrix, and P is Making by Probability Sets of each default element manipulation in kind under each mode；
Effective status set obtains module 403, for obtaining corresponding effective status set based on stochasticflow networks model G
W；
Calculation of Reliability module 404, for according to effective status set W, the relay calculated under stochasticflow networks model G to be protected
Protect reliabilityWherein, l is lower boundary point Y in effective status set W_{d}Number.
Reliability of relay protection computing device disclosed by the embodiments of the present invention obtains module according to each in intelligent substation
The historical traffic of default material object element obtains multimodal capacityconstrained matrix and each default element manipulation in kind under each mode
Making by Probability Sets；Model construction module building node can keep to the side failure stochasticflow networks model；Effective status set obtains mould
Block is based on stochasticflow networks model and obtains corresponding effective status set；Calculation of Reliability module obtains module according to effective status
Set calculates the reliability of relay protection under stochasticflow networks model.Based on device disclosed above, using stochasticflow networks mould
Type is to the Calculation of Reliability of relay protection, to achieve the purpose that carry out complicated electric power system reliability quantitative calculating.
Example IV
Reliability of relay protection computing device and attached drawing 4 in conjunction with disclosed in abovedescribed embodiment three, the present embodiment four are also disclosed
A kind of reliability of relay protection computing device, wherein effective status set obtain module 403 structural schematic diagram as shown in figure 5,
It includes: state space tree generation unit 501, zero network state generation unit that effective status set, which obtains module 403,
502, it searches for element determination unit 503, determine computing unit 504, the first judging unit 505, the first execution unit 506, backtracking
Execution unit 507, second judgment unit 508, the second execution unit 509, third judging unit 510 and third execution unit 511；
State space tree generation unit 501 for obtaining each network state based on stochasticflow networks model G, and generates
State space tree, wherein each network state is the node of the state space tree；
Zero network state generation unit 502 for generating zero network state based on network state, and zero network state is made
For the first father node of state space tree；
Element determination unit 503 is searched for, the default element in kind of one for determining network state in state space tree is made
For current search element；
Determine computing unit 504, for determining in state space tree that next network state is current search node, and root
The information source node of state space tree and the minimal path sets of information meeting point are obtained according to default minimal path sets calculation method；
First judging unit 505, for judging that the network of the information source node and information meeting point is opened up according to minimal path sets
It flutters and whether is connected to；If so, triggering second judgment unit 508；If not, the first execution unit 506 of triggering；
First execution unit 506 for triggering backtracking execution unit 507, and enters and determines computing unit 504；
Second judgment unit 508, for obtaining current search node according to the maximum stream calculation of maximumflow algorithm starting is preset
Network state capacity and judge whether network state capacity is less than preset need capacity；If it is not, the second execution unit of triggering
509；If so, triggering third judging unit 510；
Second execution unit 509, for network state capacity to be determined as lower boundary point Y_{d}It is stored in state set L, touches
The execution unit 507 that traces back is sent back to, and enters and determines computing unit 504；
Third judging unit 510, for judging whether the mode serial number of current search node is less than preset mode serial number most
Big value；It is to trigger and determine computing unit 504；It is no, trigger third execution unit 511；
Third execution unit 511 for triggering the backtracking execution unit, and enters and determines computing unit 504；
Wherein, the structural schematic diagram for recalling execution unit 507 is as shown in Figure 6, comprising: the 4th judging unit the 601, the 4th is held
Row unit 602, the 5th judging unit 603, the 5th execution unit 604 and the 6th execution unit 605；
4th judging unit 601, for judge current search element whether be state space tree final search element；Such as
Fruit is no, triggers the 4th execution unit 602；If so, the 5th judging unit 603 of triggering；
4th execution unit 602 for being current search node by a network state upper in state space tree, and determines
Next default element in kind is as the current search element in state space tree；
5th judging unit 603 for being current search node by a network state upper in state space tree, and judges
Whether current search node is zero network state；If it is not, the 5th execution unit 604 of triggering；If so, the 6th execution unit of triggering
605；
5th execution unit 604, for determining that next default element in kind is as current search member in state space tree
Part；
6th execution unit 605 for state set L to be determined as effective status set W, and terminates maximum stream calculation.
Reliability of relay protection computing device disclosed by the embodiments of the present invention obtains module according to each in intelligent substation
The historical traffic of default material object element obtains multimodal capacityconstrained matrix and each default element manipulation in kind under each mode
Making by Probability Sets；Model construction module building node can keep to the side failure stochasticflow networks model；Effective status set obtains mould
Each unit in block is based on stochasticflow networks model and obtains corresponding effective status set；Calculation of Reliability module obtains module root
According to effective status set, the reliability of relay protection under stochasticflow networks model is calculated.Based on device disclosed above, using with
Machine flow network model carries out quantitative calculating to complicated electric power system reliability to reach to the Calculation of Reliability of relay protection
Purpose.
A kind of reliability of relay protection calculation method provided by the present invention and device are described in detail above, this
Apply that a specific example illustrates the principle and implementation of the invention in text, the explanation of above example is only intended to
It facilitates the understanding of the method and its core concept of the invention；At the same time, for those skilled in the art, think of according to the present invention
Think, there will be changes in the specific implementation manner and application range, in conclusion the content of the present specification should not be construed as pair
Limitation of the invention.
It should be noted that all the embodiments in this specification are described in a progressive manner, each embodiment weight
Point explanation is the difference from other embodiments, and the same or similar parts between the embodiments can be referred to each other.
For the device disclosed in the embodiment, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, phase
Place is closed referring to method part illustration.
It should also be noted that, herein, relational terms such as first and second and the like are used merely to one
Entity or operation are distinguished with another entity or operation, without necessarily requiring or implying between these entities or operation
There are any actual relationship or orders.Moreover, the terms "include", "comprise" or its any other variant are intended to contain
Lid nonexclusive inclusion, so that the element that the process, method, article or equipment including a series of elements is intrinsic,
It further include either the element intrinsic for these process, method, article or equipments.In the absence of more restrictions,
The element limited by sentence "including a ...", it is not excluded that in the process, method, article or equipment including the element
In there is also other identical elements.
The foregoing description of the disclosed embodiments enables those skilled in the art to implement or use the present invention.
Various modifications to these embodiments will be readily apparent to those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, of the invention
It is not intended to be limited to the embodiments shown herein, and is to fit to and the principles and novel features disclosed herein phase one
The widest scope of cause.
Claims (5)
1. a kind of reliability of relay protection calculation method, which is characterized in that be applied to reliability of relay protection computing device, packet
It includes:
Multimodal capacityconstrained matrix and each institute are obtained according to the historical traffic of default element in kind each in intelligent substation
State Making by Probability Sets of the default element manipulation in kind under each mode；
Building node can keep to the side failure stochasticflow networks model G=(V, E, C, P)；
Wherein, node V is the physical connection cluster tool of each default element in kind, and side E is the default element in kind
Set, C is the multimodal capacityconstrained matrix, and P is probability of each default element manipulation in kind under each mode
Set；
Corresponding effective status set W is obtained based on the stochasticflow networks model G；
According to the effective status set W, the reliability of relay protection under the stochasticflow networks model G is calculatedWherein, l is lower boundary point Y in the effective status set W_{d}Number；Wherein,
It is described that corresponding effective status set W is obtained based on the stochasticflow networks model G, comprising:
Each network state is obtained based on the stochasticflow networks model G, and generates state space tree, wherein each net
Network state is the node of the state space tree；
Zero network state is generated based on the network state, and using zero network state as the first of the state space tree
Father node；
Determine a default element in kind of network state described in the state space tree as current search element；
Determine that next network state is current search node in the state space tree, and according to default minimal path sets calculating side
Method obtains the information source node of the state space tree and the minimal path sets of information meeting point；
Judge whether the information source node and the network topology of the information meeting point are connected to according to the minimal path sets；
If be not connected to, back tracking operation is executed, and returns and executes next network state in the determination state space tree
For current search node, the step for；
If connection, maximum stream calculation is started according to default maximumflow algorithm, obtains the network state of the current search node
Capacity simultaneously judges whether the network state capacity is less than preset need capacity；
If the network state capacity is not less than preset need capacity, the network state capacity is determined as the lower boundary point Y_{d}
It is stored in state set L, executes the back tracking operation, and return and execute next net in the determination state space tree
Network state is current search node, the step for；
If the network state capacity is less than preset need capacity, judge whether the mode serial number of the current search node is less than
Preset mode serial number maximum value；
If the mode serial number of the current search node is less than preset mode serial number maximum value, returns and execute the determination shape
Next network state is current search node in state space tree, the step for；
If the mode serial number of the current search node is not less than preset mode serial number maximum value, the back tracking operation is executed, and
Returning and executing next network state in the determination state space tree is current search node, the step for；
Wherein, the back tracking operation includes:
Judge the current search element whether be the state space tree final search element；
If the current search element is not the final search element of the state space tree, in the state space tree
One network state is the current search node, and determines that next default element in kind is made in the state space tree
For the current search element；
If the current search element is the final search element of the state space tree, by the state space tree upper one
A network state is the current search node, and judges whether the current search node is zero network state；
If the current search node is not zero network state, next default reality in the state space tree is determined
Construction element is as the current search element；
If the current search node is zero network state, the state set L is determined as the effective status set
W, and terminate maximum stream calculation.
2. the method according to claim 1, wherein each default element manipulation in kind is under each mode
The value range of Probability p is 0~1.
3. the method according to claim 1, wherein the default minimal path sets calculation method includes: contact square
The tactical deployment of troops.
4. the method according to claim 1, wherein the default maximumflow algorithm includes: extensions path Ford
Fulkerson algorithm.
5. a kind of reliability of relay protection computing device characterized by comprising obtain module, model construction module, effective shape
State set obtains module and Calculation of Reliability module；
The acquisition module, for obtaining multimodal capacity according to the historical traffic of default element in kind each in intelligent substation
The Making by Probability Sets of constraint matrix and each default element manipulation in kind under each mode；
The model construction module, for construct node can keep to the side failure stochasticflow networks model G=(V, E, C, P)；Wherein,
Node V is the physical connection cluster tool of each default element in kind, and side E is the set of the default element in kind, and C is
The multimodal capacityconstrained matrix, P are Making by Probability Sets of each default element manipulation in kind under each mode；
The effective status set obtains module, for obtaining corresponding effective status collection based on the stochasticflow networks model G
Close W；
The Calculation of Reliability module, for calculating under the stochasticflow networks model G according to the effective status set W
Reliability of relay protectionWherein, l is lower boundary point Y in the effective status set W_{d}?
Number；Wherein,
It includes: state space tree generation unit, zero network state generation unit, search member that the effective status set, which obtains module,
Part determination unit, determine computing unit, the first judging unit, the first execution unit, backtracking execution unit, second judgment unit,
Second execution unit, third judging unit and third execution unit；
The state space tree generation unit for obtaining each network state based on the stochasticflow networks model G, and generates
State space tree, wherein each network state is the node of the state space tree；
The zero network state generation unit, for generating zero network state based on the network state, and by zero network
First father node of the state as the state space tree；
Described search element determination unit, a default reality for determining network state described in the state space tree
Construction element is as current search element；
The determining computing unit, for determining, next network state is current search node in the state space tree, and
The information source node of the state space tree and the minimal path sets of information meeting point are obtained according to default minimal path sets calculation method；
First judging unit, for judging the net of the information source node and the information meeting point according to the minimal path sets
Whether network topology is connected to；If so, triggering the second judgment unit；If not, triggering first execution unit；
First execution unit for triggering backtracking execution unit, and enters the determining computing unit；
The second judgment unit, for obtaining the current search section according to the maximum stream calculation of maximumflow algorithm starting is preset
The network state capacity of point simultaneously judges whether the network state capacity is less than preset need capacity；If it is not, triggering described second
Execution unit；If so, triggering the third judging unit；
Second execution unit, for the network state capacity to be determined as the lower boundary point Y_{d}It is stored in state set L
In, the backtracking execution unit is triggered, and trigger the determining computing unit；
The third judging unit, for judging whether the mode serial number of the current search node is less than preset mode serial number most
Big value；It is to trigger the determining computing unit；It is no, trigger the third execution unit；
The third execution unit for triggering the backtracking execution unit, and enters the determining computing unit；
The backtracking execution unit, comprising: the 4th judging unit, the 4th execution unit, the 5th judging unit, the 5th execution unit
With the 6th execution unit；
4th judging unit, for judge the current search element whether be the state space tree final search member
Part；If not, triggering the 4th execution unit；If so, triggering the 5th judging unit；
4th execution unit, for being the current search node by a network state upper in the state space tree,
And determine that next default element in kind is as the current search element in the state space tree；
5th judging unit, for being the current search node by a network state upper in the state space tree,
And judge whether the current search node is zero network state；If it is not, triggering the 5th execution unit；If so, touching
Send out the 6th execution unit described；
5th execution unit is worked as described in next default element conduct in kind for determining in the state space tree
Preceding search element；
6th execution unit for the state set L to be determined as the effective status set W, and terminates maxflow
It calculates.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201610617812.XA CN106059838B (en)  20160729  20160729  A kind of reliability of relay protection calculation method and device 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201610617812.XA CN106059838B (en)  20160729  20160729  A kind of reliability of relay protection calculation method and device 
Publications (2)
Publication Number  Publication Date 

CN106059838A CN106059838A (en)  20161026 
CN106059838B true CN106059838B (en)  20190312 
Family
ID=57196868
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201610617812.XA Expired  Fee Related CN106059838B (en)  20160729  20160729  A kind of reliability of relay protection calculation method and device 
Country Status (1)
Country  Link 

CN (1)  CN106059838B (en) 
Families Citing this family (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN109861395B (en) *  20190306  20220617  中国南方电网有限责任公司  Method for dynamically controlling state machine of intelligent substation relay protection equipment 
CN110210531B (en) *  20190515  20211102  北京航空航天大学  Fuzzy multistate manufacturing system task reliability evaluation method based on extended random flow network 
Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN102882273A (en) *  20120906  20130116  广东省电力调度中心  Quantitative calculation method and system for reliability of relay protection system of intelligent substation 
CN102945317A (en) *  20121026  20130227  华北电力大学  Reliability assessment method for relay protection device in consideration of software and human factors 
CN104616212A (en) *  20150206  20150513  广东电网有限责任公司电力调度控制中心  Relay protection system reliability analysis method and system 
CN104680020A (en) *  20150311  20150603  上海毅昊自动化有限公司  SCDbased relay protection system reliability online evaluation system 

2016
 20160729 CN CN201610617812.XA patent/CN106059838B/en not_active Expired  Fee Related
Patent Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN102882273A (en) *  20120906  20130116  广东省电力调度中心  Quantitative calculation method and system for reliability of relay protection system of intelligent substation 
CN102945317A (en) *  20121026  20130227  华北电力大学  Reliability assessment method for relay protection device in consideration of software and human factors 
CN104616212A (en) *  20150206  20150513  广东电网有限责任公司电力调度控制中心  Relay protection system reliability analysis method and system 
CN104680020A (en) *  20150311  20150603  上海毅昊自动化有限公司  SCDbased relay protection system reliability online evaluation system 
NonPatent Citations (1)
Title 

电力系统继电保护课程设计模式的探索与实践;何瑞文 等;《电力系统及其自动化学报》;20090630;第21卷(第3期);第125128页 
Also Published As
Publication number  Publication date 

CN106059838A (en)  20161026 
Similar Documents
Publication  Publication Date  Title 

Huang et al.  Bad data injection in smart grid: attack and defense mechanisms  
Huang et al.  Modeling cascading failures in smart power grid using interdependent complex networks and percolation theory  
Zhao et al.  Power system fault diagnosis based on history driven differential evolution and stochastic time domain simulation  
Wang et al.  System islanding using minimal cutsets with minimum net flow  
Cavalieri et al.  Models for seismic vulnerability analysis of power networks: comparative assessment  
CN103308824B (en)  Power system fault diagnostic method based on probability Petri net  
Lin et al.  Division algorithm and interconnection strategy of restoration subsystems based on complex network theory  
Chopade et al.  Critical infrastructure interdependency modeling: Using graph models to assess the vulnerability of smart power grid and SCADA networks  
CN109217296A (en)  A kind of power grid key node discrimination method based on weighting trend entropy and voltage stabilization  
CN104734870B (en)  A kind of software fault propagation law discovery method based on cellular automata  
CN109767352A (en)  A kind of power information physics emerging system safety situation evaluation method  
CN106056466B (en)  Bulk power grid critical circuits recognition methods based on FPgrowth algorithm  
CN106059838B (en)  A kind of reliability of relay protection calculation method and device  
CN107247450A (en)  Circuit breaker failure diagnostic method based on Bayesian network  
CN110350522A (en)  A kind of electric system vulnerable line identifying method based on Weighted H index  
Deka et al.  Data attack on strategic buses in the power grid: Design and protection  
Zhu et al.  Riskaware vulnerability analysis of electric grids from attacker's perspective  
CN108183897A (en)  A kind of information physical emerging system safety risk estimating method  
CN110022293A (en)  A kind of electric network information physics emerging system methods of risk assessment  
CN108090677A (en)  A kind of critical infrastructures reliability assessment method  
CN107450517A (en)  A kind of safe diagnosability determination method of Stochastic discrete event systems and system  
CN110048901A (en)  A kind of Fault Locating Method of power telecom network, device and equipment  
CN110380903A (en)  A kind of power telecom network fault detection method, device and equipment  
Zhou et al.  Robustness analysis of power system dynamic process and repair strategy  
Atkins et al.  The structure of electrical networks: a graph theory based analysis 
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  
CF01  Termination of patent right due to nonpayment of annual fee  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20190312 Termination date: 20200729 