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 PDF

Info

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
Application number
CN201610617812.XA
Other languages
Chinese (zh)
Other versions
CN106059838A (en
Inventor
徐庆焯
贺小旺
何瑞文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201610617812.XA priority Critical patent/CN106059838B/en
Publication of CN106059838A publication Critical patent/CN106059838A/en
Application granted granted Critical
Publication of CN106059838B publication Critical patent/CN106059838B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L41/00Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
    • H04L41/14Network analysis or design
    • H04L41/145Network analysis or design involving simulating, designing, planning or modelling of a network
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L41/00Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
    • H04L41/08Configuration management of networks or network elements
    • H04L41/0803Configuration setting
    • H04L41/0823Configuration setting characterised by the purposes of a change of settings, e.g. optimising configuration for enhancing reliability
    • H04L41/0836Configuration 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 multi-modal capacity-constrained 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 stochastic-flow networks model;Corresponding effective status set is obtained based on stochastic-flow networks model;According to effective status set, the reliability of relay protection under stochastic-flow networks model is calculated.Based on method disclosed above, using stochastic-flow 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

A kind of reliability of relay protection calculation method and device
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 built-in self-diagnostic 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 multi-modal capacity-constrained 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 stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow networks model G;
According to the effective status set W, the reliability of relay protection under the stochastic-flow networks model G is calculatedWherein, l is lower boundary point Y in the effective status set WdNumber.
It is preferably, described that corresponding effective status set W is obtained based on the stochastic-flow networks model G, comprising:
Each network state is obtained based on the stochastic-flow 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 maximum-flow 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 YdIt 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 maximum-flow algorithm includes: extensions path Ford-Fulkerson 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 multi-modal for being obtained according to the historical traffic of default element in kind each in intelligent substation The Making by Probability Sets of capacity-constrained matrix and each default element manipulation in kind under each mode;
The model construction module, for construct node can keep to the side failure stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow networks model G State set W;
The Calculation of Reliability module, for calculating the stochastic-flow 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 Wd 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 stochastic-flow 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 maximum-flow 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 YdIt 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 multi-modal capacity-constrained matrix and each default element manipulation in kind in each mode Under Making by Probability Sets;Building node can keep to the side failure stochastic-flow networks model;It is obtained based on stochastic-flow networks model corresponding Effective status set;According to effective status set, the reliability of relay protection under stochastic-flow networks model is calculated.Based on above-mentioned public affairs The method opened, using stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow networks model of failure, The physical connection cluster tool V={ v of each default element in kindi| 1≤i≤n } node as stochastic-flow networks model, in advance If the set E={ e of element in kindi| 1≤i≤a } side as stochastic-flow networks model, C={ cij|1≤i≤n+a,0≤j≤ hiIt is used as multi-modal capacity-constrained matrix, P={ pij|1≤i≤n+a,0≤j≤hiIt is that each default element manipulation in kind exists Making by Probability Sets under each mode, wherein hiFor the corresponding mode serial number of i-th of element in element;
S103 obtains corresponding effective status set W based on stochastic-flow networks model G;
S104 calculates the reliability of relay protection under stochastic-flow networks model G according to effective status set WWherein, l is lower boundary point Y in effective status set WdNumber.
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 multi-modal capacity-constrained matrix and each default element manipulation in kind under each mode It closes;Building node can keep to the side failure stochastic-flow networks model;Corresponding effective status collection is obtained based on stochastic-flow networks model It closes;According to effective status set, the reliability of relay protection under stochastic-flow networks model is calculated.Based on method disclosed above, Using stochastic-flow 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 multi-modal 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 stochastic-flow networks model G, such as Fig. 2 institute Show, comprising:
S201 obtains each network state based on stochastic-flow 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 stochastic-flow networks model G1, Y2,…,Yi,Ym), and state space tree is generated according to each network state, wherein each network state YiIt 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 YiGenerate 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 maximum-flow 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 S208dIt 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 maximum-flow algorithm includes but is not limited to extensions path Ford- Fulkerson algorithm, the method that suitable calculating max-flow 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 multi-modal capacity-constrained matrix and each default element manipulation in kind under each mode It closes;Building node can keep to the side failure stochastic-flow networks model;Corresponding effective status collection is obtained based on stochastic-flow networks model It closes;According to effective status set, the reliability of relay protection under stochastic-flow networks model is calculated.Based on method disclosed above, Using stochastic-flow 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 above-mentioned 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 multi-modal for being obtained according to the historical traffic of default element in kind each in intelligent substation The Making by Probability Sets of capacity-constrained matrix and each default element manipulation in kind under each mode;
Model construction module 402, for construct node can keep to the side failure stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow networks model G W;
Calculation of Reliability module 404, for according to effective status set W, the relay calculated under stochastic-flow networks model G to be protected Protect reliabilityWherein, l is lower boundary point Y in effective status set WdNumber.
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 multi-modal capacity-constrained 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 stochastic-flow networks model;Effective status set obtains mould Block is based on stochastic-flow 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 stochastic-flow networks model.Based on device disclosed above, using stochastic-flow 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 above-described 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 stochastic-flow 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 maximum-flow 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 YdIt 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 multi-modal capacity-constrained 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 stochastic-flow networks model;Effective status set obtains mould Each unit in block is based on stochastic-flow 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 stochastic-flow 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 non-exclusive 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:
Multi-modal capacity-constrained 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 stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow networks model G;
According to the effective status set W, the reliability of relay protection under the stochastic-flow networks model G is calculatedWherein, l is lower boundary point Y in the effective status set WdNumber;Wherein,
It is described that corresponding effective status set W is obtained based on the stochastic-flow networks model G, comprising:
Each network state is obtained based on the stochastic-flow 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 maximum-flow 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 Yd 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 maximum-flow 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 multi-modal 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 stochastic-flow 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 multi-modal capacity-constrained 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 stochastic-flow networks model G Close W;
The Calculation of Reliability module, for calculating under the stochastic-flow 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 Wd? 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 stochastic-flow 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 maximum-flow 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 YdIt 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 max-flow It calculates.
CN201610617812.XA 2016-07-29 2016-07-29 A kind of reliability of relay protection calculation method and device Expired - Fee Related CN106059838B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610617812.XA CN106059838B (en) 2016-07-29 2016-07-29 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) 2016-07-29 2016-07-29 A kind of reliability of relay protection calculation method and device

Publications (2)

Publication Number Publication Date
CN106059838A CN106059838A (en) 2016-10-26
CN106059838B true CN106059838B (en) 2019-03-12

Family

ID=57196868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610617812.XA Expired - Fee Related CN106059838B (en) 2016-07-29 2016-07-29 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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109861395B (en) * 2019-03-06 2022-06-17 中国南方电网有限责任公司 Method for dynamically controlling state machine of intelligent substation relay protection equipment
CN110210531B (en) * 2019-05-15 2021-11-02 北京航空航天大学 Fuzzy multi-state manufacturing system task reliability evaluation method based on extended random flow network

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102882273A (en) * 2012-09-06 2013-01-16 广东省电力调度中心 Quantitative calculation method and system for reliability of relay protection system of intelligent substation
CN102945317A (en) * 2012-10-26 2013-02-27 华北电力大学 Reliability assessment method for relay protection device in consideration of software and human factors
CN104616212A (en) * 2015-02-06 2015-05-13 广东电网有限责任公司电力调度控制中心 Relay protection system reliability analysis method and system
CN104680020A (en) * 2015-03-11 2015-06-03 上海毅昊自动化有限公司 SCD-based relay protection system reliability online evaluation system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102882273A (en) * 2012-09-06 2013-01-16 广东省电力调度中心 Quantitative calculation method and system for reliability of relay protection system of intelligent substation
CN102945317A (en) * 2012-10-26 2013-02-27 华北电力大学 Reliability assessment method for relay protection device in consideration of software and human factors
CN104616212A (en) * 2015-02-06 2015-05-13 广东电网有限责任公司电力调度控制中心 Relay protection system reliability analysis method and system
CN104680020A (en) * 2015-03-11 2015-06-03 上海毅昊自动化有限公司 SCD-based relay protection system reliability online evaluation system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
电力系统继电保护课程设计模式的探索与实践;何瑞文 等;《电力系统及其自动化学报》;20090630;第21卷(第3期);第125-128页

Also Published As

Publication number Publication date
CN106059838A (en) 2016-10-26

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 FP-growth 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. Risk-aware 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 non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190312

Termination date: 20200729