CN103957032B - Load redistribution method for electric power coupling network to resist cascade failure - Google Patents
Load redistribution method for electric power coupling network to resist cascade failure Download PDFInfo
- Publication number
- CN103957032B CN103957032B CN201410156420.9A CN201410156420A CN103957032B CN 103957032 B CN103957032 B CN 103957032B CN 201410156420 A CN201410156420 A CN 201410156420A CN 103957032 B CN103957032 B CN 103957032B
- Authority
- CN
- China
- Prior art keywords
- node
- communication
- load
- network
- net
- 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
Abstract
The invention discloses a load redistribution method for an electric power coupling network to resist cascade failure. The load redistribution method for the electric power coupling network to resist the cascade failure comprises the steps of abstracting and simplifying equipment in a communication network and a power network as a coupling network topology model; establishing a node capacity model according to the characteristics of the coupling network; establishing a load redistribution model according to a distribution method of loads after node failure; establishing a cascade failure model according to the characteristics of node failure of the coupling network; attacking communications network nodes, obtaining the maximum function set of power grid nodes after the coupling network is stabilized, and selecting the optimal load redistribution method. The method improves the capacity of the electric power coupling network for resisting the cascade failure, is beneficial to reducing maintenance cost, and improves the economic benefits and the social benefits of an electrical power system.
Description
Technical field
The invention belongs to electric power networks and communication network resist ruins technical field, resist particularly to a kind of electric power coupling network
Cascading failure load re-distribution method.
Background technology
In recent years, the problem in some reality, the congestion of such as transportation network, electric power networks large-scale blackout etc., attract
The interest of a lot of scholar's research Complex Networks Theory, is greatly promoted prosperity and the development of Network Science.When very long one section
Between people study limitation in the research of single network, do not consider the impact to other networks for this network and other networks to this net
The impact of network.Going deep into it has been found that the relation between some networks is very close with research, the characteristic of a network is subject to
The impact of another or several network characteristic, the interaction between network gradually embodies, and we influence each other this
Network be referred to as coupling network.Electric power networks and communication network are exactly typical coupling network, and electrical network changes to intelligent, electrical network
The more intelligent coupling with communication network is stronger.The normal operation of electric power networks needs communication network to transmit the letter such as Monitoring and Controlling
Cease, communication network device needs electric power networks to provide electric power support for it simultaneously.The node failure of communication network is possible to cause electricity
Net interdependent node lost efficacy, and electrical network failure node may lead to the node failure of communication network in turn again, if under such cascade
Go, a glitch may finally cause significant impact to whole system.It is impact that load after node failure reassigns problem
The key factor of coupling network, chaotic reassigns the help very little to cascading failure for the strategy, it could even be possible to promoting cascade to lose
The generation of effect, makes coupling network more fragile, and the present invention selects appropriate load distribution plan according to the self character of coupling network
Slightly, improve the ability that coupling network resists cascading failure.
Content of the invention
The purpose of the present invention is to propose to a kind of electric power coupling network resists cascading failure load re-distribution method, its feature
It is, comprise the following steps：
Step 1：Device abstract in communication network and electric power networks is reduced to coupling network topological model G；
Step 2：Node capacity model is set up according to the feature of coupling network；
Step 3：Distribution method according to its load after node failure is set up load and is reassigned model；
Step 4：Cascading failure model is set up according to coupling network node failure feature；
Step 5：Node at random or in calculated attack communication network, obtains electricity after cascading failure under different allocation rule terminates
Maximum function group GNZP of net node, chooses optimal load redistribution method.
Described step 1 specifically includes following steps：
Step 1.1：Obtain the generator of normal work and 35kV and above transformer station in electrical network and be abstracted into power node
Vp；
Step 1.2：Obtain the communications equipment room of powerline network and wireless base station is abstracted into communication node Vc；
Step 1.3：The high voltage transmission line obtaining between grid nodes is abstracted into the side E between nodep, do not consider power distribution network
With the main electrical scheme in power plant, transformer station, merge the power transmission line of parallel lines on same tower, eliminate from ring and multiple circuit, so far electric power networks
Topological model builds and finishes, and uses Gp=(Vp,Ep) represent；
Step 1.4：Obtain the communication line between communication node or wireless channel is abstracted into the side E between communication nodec, extremely
This communication network topology model construction finishes, and uses Gc=(Vc,Ec) represent；
Step 1.5：Obtain the grid nodes that electric power is provided to each communication node, be abstracted into grid nodes to communication network section
The unidirectional side E of pointpc；
Step 1.6：Obtain the communication node that communication support is provided to each power node, be abstracted into communication network and arrive electrical network
Unidirectional side Ecp, so far coupling network topological model build finish, with G=(Gp,Gc,Epc,Ecp) represent；
In described step 1.2, communication node includes：Node and the Business Stream of communication service are provided for abstract power node
The node of up process, provides the node of communication service as abstract communication node for the abstract following electrical equipment of power node
Load.
Described step 2 specifically includes following steps：
Step 2.1：Determine the capacity model of communication net node
, in electrical network, electrical network normal work node is more, and the load of communication network is bigger for the load dependency of communication network.Alternately analyze
The cascading failure process of coupling network, if iteration n times cascading failure stops, nth iteration is referred to as state n；Wherein, n=0,1,
2,3 ..., N；
The total traffic that communication network undertakes within the state n unit interval is：
Wherein：
Φ represents the node set of normal work in state n electrical network,
kpniExpression statenThe Connected degree of grid nodes i,
Represent that grid nodes produce the factor of influence of message volume,
The traffic list of communication node i is shown as：
Wherein：
Ι represents the node set of normal work in state n communication network,
kcniThe Connected degree of expression state n communication net node i,
EnRepresent the Connected degree sum of Ι interior joint,
Represent that communication node undertakes the factor of influence of portfolio；
Each communication node constitutes communication node load vector in the portfolio of state n
Wherein：mcRepresent the nodes of communication network；
The maximum traffic capacity of communication net node is directly proportional to initial load, is expressed as：
Wherein：
δcFor the tolerance factor of communication net node, span (0.5-1),
CL01,kc0iPortfolio for original state communication node i and Connected degree,
W0,E0For the total portfolio of original state communication network and total Connected degree；
Step 2.2：Determine the capacity model of grid nodes
With communication network relation less, the power consumption of communication network is far smaller than the load of whole electrical network for the load of electrical network, so
When considering network load, we ignore the impact of communication network；
The average of node i load in the recent period during grid nodes i initial load power taking net normal work：
PL0i=(L1+L2+…+Ln)/n
Wherein：
LnIt is the average load of time period n interior nodes；
The load of state n grid nodes i constitutes grid nodes load vector：
Wherein：mpRepresent the node number of electrical network.
The capacity of grid nodes i is：
PCi=(1+ δp)PLe
Wherein：
PLeIt is rated capacity or the rated load of power node, each grid nodes gain of parameter can be consulted；
δpFor the tolerance factor of grid nodes, span (0.3-0.5).
Described step 3 is set up load reassignment model and is specifically included following steps：
Step 3.1：Determine the load vector set of the set S of node, the set Z of normal work node and network that lost efficacy
L；
Step 3.2：After determining that node i lost efficacy, incremental loading on node j：
Wherein:I is the set S interior joint of the node losing efficacy, and j is the set Z interior joint of normal work node, LiFor load
The load of vectorial set L interior joint i, L (Dij,θ,kj, φ) and represent incremental loading defined function, DijFor between two nodes away from
From kjNode degree for j, θ, φ are that load reassigns factor of influence, control scope and the uniformity of reassignment respectively, ε is negative
Carry distribution coefficient, Ω represents the set of all normal work nodes；
Step 3.3：Calculate node j loads total increment：
Wherein：
Φ represents the set of all failure nodes.
Step 3.4：The load capacity of each node after computational load reassignment：
Lj=Lj+ΔLj,
And and capacity CiRelatively, CiMaximum capacity for this node device；If Lj>Cj, it is considered as this node failure, every
The load arranging failure node after individual state is 0, is no longer participate in load and reassigns.
Reassign in model in load, adjust the value reassigning range parameter θ, allocation rule is divided into and faces distribution (θ recently
=inf), intermediate state (0<θ<Inf) and three kinds of situations of global assignment (θ=0), inf represents positive infinity.
In described step 4, coupling network cascading failure model concretely comprises the following steps：
Step 4.1：Attack the node in communication network, cause the inefficacy on communication net node and side, until communication network reaches surely
Determine state, update connection matrix E between netcp, this process is called state 1；
Step 4.2：If connection matrix E in state 1cpChange, electrical network net interior nodes and the inefficacy on side may be caused, reach
To stable state, status number adds 1, judges to connect E between netpcWhether changing, if changing, updating Epc, execute next step；If no
Become then cascading failure to terminate, jump to step 4.4；
Step 4.3：If connection matrix E in state 1pcChange, cause communication network net interior nodes and the inefficacy on side, reach
Stable state, status number adds 1, judges to connect E between netcpWhether changing, if changing, updating Ecp, return to step 4.2；If constant
Then cascading failure terminates, and executes next step；
Step 4.4：Coupling network cascading failure stops, and obtains maximum function group GNZP of grid nodes, and GNZP represents level
The abstractdesription of connection failure procedure.
In described cascading failure model, when analyzing the failure procedure on state n net interior nodes and side it is assumed that another network
Keep the stable state of n-1 constant；
In state n, the inefficacy on net interior nodes and side is specially：
The node failure of state n includes：Lose and between net, connect (Ecp, Epc) node failure, lose net in (Ec, Ep) connect section
Point failure and overload lost efficacy, and its process uses P respectivelyO, PIAnd PLRepresent.Wherein PIAnd PLIt is an interaction, lose in net and connect
Connect node failure, load reassignment can be caused to lead to node overload to lose efficacy, overload lost efficacy can lose in net again and connects, and led to node
Lost efficacy, did so with until stable state, so the maximum function group after state n is stable is represented by GNZn=Xn-2POPIPL=
Xn-2POPIL, Xn-2Represent the maximum function group of stable state n-2 in net, n-2 state to n-state needs cascading failure in experience net
Process, PILCause the node overload failure procedure of load reassignment for net internal segment point failure.
Between described loss net, connecting node lost efficacy and was specially：
Communication net node is lost Joint failure finger joint point between net and is lost electric power support, is considered as this communication node in cascading failure
Lost efficacy in time period；Grid nodes are lost Joint failure finger joint point between net and are lost communication support, lead to equipment fault, are considered as this and lead to
Letter node lost efficacy within the cascading failure time period.
In described loss net, connecting node lost efficacy and was specially：
Communication network or grid nodes lose the connection with present networks other node, become isolated node, are considered as this node
Lost efficacy within the cascading failure time period.
The inefficacy of described overload is specially：
The packet that communication network interior joint overload inefficacy finger joint point receives leads to mass data bag much larger than its disposal ability
It is dropped, be considered as this communication node and lost efficacy within the cascading failure time period；The overload of grid nodes lost efficacy and referred to that transformer exceeded
Its maximum capacity or generator exceed its peak load and lead to equipment fault, are considered as this power node within the cascading failure time period
Lost efficacy.
Described step 5 specifically includes：
Node at random or in calculated attack communication network, obtains the maximum function group of stable rear electrical network under different allocation rule
GNZP, selects optimum load re-distribution method.
The invention has the beneficial effects as follows reassigning unreasonable, the accident such as cascading failure causing for coupling network load,
A kind of electric power coupling network proposing resists cascading failure load re-distribution method, and this method has deeply considered enforcement regional power grid
With the topological structure of communication net node, according to nodal distance and Connected degree, the load of failure node is reassigned, choose
Excellent distribution method, thus allow coupling network have higher survivability.
Brief description
Fig. 1 is the flow chart that electric power coupling network resists cascading failure load re-distribution method；
Fig. 2 is coupling network cascading failure illustraton of model；
Fig. 3 is the illustraton of model of state n cascading failure.
Specific embodiment
For further illustrating the effect of technological means that the present invention taken and acquirement, below in conjunction with the accompanying drawings, to this
Bright technical scheme, carries out clear and complete description.
Below in conjunction with the accompanying drawings, preferred embodiment is elaborated.It is emphasized that the description below is merely exemplary
, rather than in order to limit the scope of the present invention and its application.
Fig. 1 resists the schematic flow sheet of cascading failure load re-distribution method for electric power coupling network of the present invention.The present invention
Including several steps once：
S1：By abstract for coupling network for topological model；
In acquisition methods enforcement regional power grid, the generator of normal work and the transformer station of more than 35KkV are abstracted into electricity first
Power node Vp；Obtain the communications equipment room of powerline network and wireless base station is abstracted into communication node Vc；Obtain between grid nodes
High voltage transmission line be abstracted into the side E between nodep, do not consider power distribution network and the main electrical scheme in power plant, transformer station, merge same bar
And the power transmission line of frame, elimination builds and finishes from ring and multiple circuit, so far electric power networks topological model, uses Gp=(Vp,Ep) table
Show；Obtain the communication line between communication node or wireless channel is abstracted into the side E between communication nodec, so far communication network open up
Flutter model construction to finish, use Gc=(Vc,Ec) represent；Obtain the grid nodes that electric power is provided to each communication node, be abstracted into electricity
Net node is to the unidirectional side E of communication net nodepc；Obtain the communication node that communication support is provided to each power node, be abstracted into
Communication network is to the unidirectional side E of electrical networkcp, so far coupling network topological model build finish, with G=(Gp,Gc,Epc,Ecp) represent.One
Individual simple embodiment model construction finishes as shown in Fig. 2 state 0.
S2：Set up node capacity model；
Node capacity model includes communication net node capacity model and grid nodes capacity model.Due to communication network and electrical network
It is coupling network, offered load respectively has its feature.
Communication net node capacity model：
Load in electrical network, the more communication networks of electrical network normal work node for the load dependency of communication network is bigger, state n,
Total message volume in communication network in unit interval is：
Wherein：
Φ represents the node set of normal work in state n electrical network,
kpniThe Connected degree of expression state n grid nodes i,
Represent the factor of influence of the message volume that grid nodes produce, respectively value 10 and 1.
The traffic list of communication node i is shown as：
Wherein：
Ι represents the node set of normal work in state n communication network,
kcniThe Connected degree of expression state n communication net node i,
EnRepresent the Connected degree sum of Ι interior joint,
Represent that communication node undertakes the factor of influence of portfolio,
Each communication node constitutes communication node load vector in the portfolio of state n
Wherein：mcRepresent the nodes of communication network.
The maximum traffic capacity of communication net node is directly proportional to initial load, is expressed as：
Wherein：
δcFor the tolerance factor of communication net node, δc=0.8,
CL01,kc0iPortfolio for original state communication node i and Connected degree,
W0,Ec0For the total portfolio of original state communication network and total Connected degree.
Connect and between net in conjunction with connecting in the net of Fig. 2 state 1 figure each point, can be obtained by above formula：
W0=140, Ec0=12, CL0=(23.3,23.3,11.7,35,23.3,23.3),
CC=(42,42,21,63,42,42).
Step 2.2：Determine the capacity model of grid nodes
With communication network relation less, the power consumption of communication network is far smaller than the load of whole electrical network for the load of electrical network, so
When considering network load, we ignore the impact of communication network.
The average of node i load in the recent period during grid nodes i initial load power taking net normal work：
PL0i=(L1+L2+…+Ln)/n
Wherein：
LnIt is the mean power of time period n interior nodes.
The load of state n grid nodes i constitutes grid nodes load vector：
Wherein：mpRepresent the node number of electrical network.
The capacity of grid nodes i is：
PCi=(1+ δp)PLe
Wherein：
PLeIt is rated capacity or the rated load of power node, each grid nodes gain of parameter can be consulted.
δpFor the tolerance factor of grid nodes, δp=0.3.
According to the coupling model of Fig. 2 state 1, we are at assignment grid nodes rated capacity：
PLe=(3150,6300,8000,1800,6300,5000), assignment grid nodes initial load is：
PL0=(2205,4410,5600,1500,4410,3500), obtain grid nodes capacity：
PC=(4095,8190,10400,2340,8190,6500), unit is KVA.
S3：Set up load and reassign model；
The present invention is loaded using identical in communication network and electrical network and reassigns model.
After node i lost efficacy, the load L of node iiOther nodes can be distributed to, the incremental loading of node j is：
Wherein, L (Dij,θ,kj, φ) and it is incremental loading defined function, DijFor the Floyd distance between i, j two node, kj
Node degree for j, θ, φ are that load reassigns factor of influence, control scope and the uniformity of reassignment respectively, and ε is load distribution
Coefficient, Ω represents the set of all normal work nodes.
Calculate node j loads total increment：
Φ represents the set of all failure nodes.
The load capacity of each node after load reassignment：Lj=Lj+ΔLj, make LjWith capacity CjRelatively, CjFor this node device
Maximum capacity, if Lj>Cj, it is considered as this node failure, update load vector, finally that failure node is negative in each state
Load is set to 0.
Reassign in model in load, the value that the value of φ is set to 1, ε is only set to behind 0.6, that is, node failure
The load having 60% can be assigned to other nodes.Adjust the value reassigning range parameter θ, allocation rule is divided into and faces recently point
Join (θ=inf), intermediate state (0<θ<Inf) and three kinds of situations of global assignment (θ=0), inf represents positive infinity.In the middle of wherein
Situation can arrange several values, such as θ=(1,2,4,10,50,100).
S4：Set up cascading failure model；
The cascading failure process of coupling network is divided into several states, in state n of analysis it is assumed that another network keeps
The stable state of state n-1 is constant.
Fig. 2 is an example of coupling network cascading failure：
State 1：Communication net node Vc1And Vc2Lost efficacy, cascading failure starts, in all nets being connected with this two nodes
Connect Joint failure and between net, node Vc3Lose the Joint failure with maximum function group；
State 2：Connect between some nets because laststate loses, lead to grid nodes Vp1And Vp2Lost efficacy it is assumed that due to
Load reassignment leads to grid nodes Vp4Lost efficacy, all Joint failures with this three point attachment；
State 3：Grid nodes lost efficacy and do not cause communication net node to lose efficacy further, and cascading failure stops.
Fig. 3 is the failure procedure of the node of each state：The node failure of state n includes：Lose and between net, connect (Ecp,
Epc) node failure, lose net in (Ec, Ep) connecting node lost efficacy and overload lost efficacy, its process uses P respectivelyO, PIAnd PLRepresent.Its
Middle PIAnd PLIt is an interaction, loses connecting node in net and lost efficacy, load reassignment can be caused to lead to node overload to lose efficacy,
Overload lost efficacy can lose in net again and connected, and led to node failure, did so with until stable state, so after state n is stable
Maximum function group is represented by GNZn=Xn-2POPIPL=Xn-2POPIL.
S5：Attack communication net node, obtain the maximum function group of grid nodes, choose optimal load redistribution method.
Attack and can be attacked and calculated attack using random according to regional feature during communication net node.To with calculated attack,
Cascading failure obtains the maximum function group of electrical network under different allocation rule after stopping, choosing the most person of node as this coupling network
Optimum load re-distribution method.For random attack can attack p=(0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,
0.9) node of ratio, obtains maximum function group GNZP of each p value under each redistribution method, draws p-GNZP curve, chooses
Excellent load re-distribution method.
Only combine accompanying drawing 2 coupling network model with calculated attack below to be illustrated, load reassigns rule and takes faces recently
Distribution (θ=inf), intermediate state (θ=1) and three kinds of situations of global assignment (θ=0), calculated attack makes communication net node Vc1With
Vc2Lost efficacy.According to cascading failure model of the present invention, when communication network choosing faces the method for salary distribution recently, three kinds of distribution methods of electrical network are remaining
Node ratio is 0.67,0.5,0.5 respectively；When communication network chooses distribution condition, the remaining node of three kinds of distribution methods of electrical network
Ratio is 0.67,0.5,0.5 respectively；When communication network selects global assignment situation, the remaining node ratio of three kinds of distribution methods of electrical network
It is 0.67,0.5,0.5 respectively.As can be seen here in the coupling model of Fig. 2 state 0, the load of electrical network reassigns mode and chooses
Closely face distribution effects preferably, communication network can randomly select a kind of method of salary distribution.
The above, the only present invention preferably specific embodiment, but protection scope of the present invention is not limited thereto,
Any those familiar with the art the invention discloses technical scope in, the change or replacement that can readily occur in,
All should be included within the scope of the present invention.Therefore, protection scope of the present invention should be with scope of the claims
It is defined.
Claims (1)
1. a kind of electric power coupling network resists cascading failure load re-distribution method it is characterised in that the method includes following step
Suddenly：
Step 1：It is reduced to coupling network topological model G by abstract for the facility information obtaining in communication network and electric power networks；Tool
Body includes：
Step 1.1：Obtain the generator of normal work and 35kV and above substation information in electrical network, be abstracted into electric power
Node Vp；
Step 1.2：Obtain communications equipment room and the wireless base station information of powerline network, be abstracted into communication node Vc；Wherein
Communication node includes：There is provided the node of communication service and its node of the up process of Business Stream for abstract power node, for taking out
The following electrical equipment of power node of elephant provides the load as abstract communication node for the node of communication service；
Step 1.3：The high voltage transmission line obtaining between grid nodes is abstracted into the side E between nodep, do not consider power distribution network and send out
Power plant, the main electrical scheme of transformer station, merge the power transmission line of parallel lines on same tower, eliminate from ring and multiple circuit, so far electric power networks topology
Model construction finishes, and uses Gp=(Vp,Ep) represent；
Step 1.4：Obtain the communication line between communication node or wireless channel is abstracted into the side E between communication nodec, so far lead to
Communication network topological model builds and finishes, and uses Gc=(Vc,Ec) represent；
Step 1.5：Obtain the grid nodes that electric power is provided to each communication node, be abstracted into grid nodes and arrive communication net node
Unidirectional side Epc；
Step 1.6：Obtain the communication node that communication support is provided to each power node, be abstracted into communication network unidirectional to electrical network
Side Ecp, so far coupling network topological model build finish, with G=(Gp,Gc,Epc,Ecp) represent；
Step 2：Node capacity model is set up according to the feature of coupling network；Specifically include：
Step 2.1：Determine the capacity model of communication net node
, in electrical network, electrical network normal work node is more, and the load of communication network is bigger for the load dependency of communication network；Alternately analysis coupling
The cascading failure process of network, if iteration n times cascading failure stops, nth iteration is referred to as state n；Wherein, n=0,1,2,
3…,N；
The total traffic that communication network undertakes within the state n unit interval is：
Wherein：
Φ represents the node set of normal work in state n electrical network,
kpniThe Connected degree of expression state n grid nodes i,
Represent that grid nodes produce the factor of influence of message volume；
The traffic list of communication node i is shown as：
Wherein：
Ι represents the node set of normal work in state n communication network,
kcniThe Connected degree of expression state n communication net node i,
EnRepresent the Connected degree sum of Ι interior joint,
Represent that communication node undertakes the factor of influence of portfolio；
Each communication node constitutes communication node load vector in the portfolio of state n
Wherein：mcRepresent the nodes of communication network；
The maximum traffic capacity of communication net node is directly proportional to initial load, is expressed as：
Wherein：
δcFor the tolerance factor of communication net node, span (0.5-1),
CL0i,Kc0iPortfolio for original state communication node i and Connected degree,
W0,E0For the total portfolio of original state communication network and total Connected degree；
Step 2.2：Determine the capacity model of grid nodes
The average of node i load in the recent period during grid nodes i initial load power taking net normal work：
PL0i=(L1+L2+…+Ln)/n
Wherein：
LnIt is the average load of time period n interior nodes；
The load of state n grid nodes i constitutes grid nodes load vector：
Wherein：mpRepresent the node number of electrical network；
The capacity of grid nodes i is：
PCi=(1+ δp)PLe
Wherein：
PLeIt is rated capacity or the rated load of power node；
δpFor the tolerance factor of grid nodes, span (0.3-0.5)；
Step 3：Distribution method according to its load after node failure is set up load and is reassigned model；Specifically include following steps：
Step 3.1：Determine the load vector set L of the set S of node, the set Z of normal work node and network that lost efficacy；
Step 3.2：After determining that node i lost efficacy, incremental loading on node j：
Wherein:I is the set S interior joint of the node losing efficacy, and j is the set Z interior joint of normal work node, LiFor load vector
The load of set L interior joint i, L (Dij,θ,kj, φ) and represent incremental loading defined function, DijFor the distance between two nodes, kj
Node degree for j, θ, φ are that load reassigns factor of influence, control scope and the uniformity of reassignment respectively, and ε is load distribution
Coefficient, Ω represents the set of all normal work nodes；
Step 3.3：Calculate node j loads total increment：
Φ represents the set of all failure nodes；
Step 3.4：The load capacity of each node after computational load reassignment：
Lj=Lj+ΔLj,
And with its capacity CjRelatively, CjMaximum capacity for communication net node；If Lj>Cj, it is considered as this node failure, in each shape
The load arranging failure node after state is 0, is no longer participate in load and reassigns；Adjust the value reassigning range parameter θ, distribution rule
Then it is divided into three kinds of situations:Face distribution θ=inf recently, inf represents just infinite, intermediateness 0<θ<Inf and global assignment θ=0；
Step 4：Cascading failure model is set up according to coupling network node failure feature；Concretely comprise the following steps：
Step 4.1：Attack the node in communication network, cause the inefficacy on communication net node and side, until communication network reaches stable shape
State, updates connection matrix E between netcp, this process is called state 1；The inefficacy on net interior nodes and side is specially：
The node failure of state n includes：Lose and between net, connect (Ecp, Epc) node failure, lose net in (Ec, Ep) connecting node mistake
Effect and overload lost efficacy, and its process uses P respectivelyO, PIAnd PLRepresent；Wherein PIAnd PLIt is an interaction, lose and in net, connect section
Point failure, can cause load reassignment to lead to node overload to lose efficacy, and overload lost efficacy can lose in net again and connects, and led to net interior connection
Node failure, until stable state, so the maximum function group after state n terminates is expressed as GNZn=Xn-2POPIPL=Xn- 2POPIL, abstract representation cascading failure process；Wherein, Xn-2Represent the maximum function group of stable state n-2, n-2 state is to n shape
State needs cascading failure process in experience net, PILCause the node overload failure procedure of load reassignment for net internal segment point failure；
When analyzing the failure procedure of state n node and side it is assumed that another network keeps the stable state of n-1 constant；
Wherein, lose connecting node inefficacy to specifically include：
1) lose connecting node between net and lost efficacy and be：Communication net node is lost Joint failure finger joint point between net and is lost electric power support, depending on
Lost efficacy within the cascading failure time period for this communication node；Grid nodes are lost Joint failure finger joint point between net and are lost communication
Hold, lead to equipment fault, be considered as this grid nodes and lost efficacy within the cascading failure time period；
2) lose connecting node in net and lost efficacy and be：Communication network or grid nodes lose the connection with present networks other node, become
Isolated node, is considered as this node and lost efficacy within the cascading failure time period；
Wherein, transshipping inefficacy is：The packet that communication network interior joint overload inefficacy finger joint point receives is led much larger than its disposal ability
Cause mass data to be coated discarding, be considered as this communication node and lost efficacy within the cascading failure time period；The overload of grid nodes lost efficacy
Refer to transformer and exceed the load of its maximum capacity or generator exceed its peak load and lead to equipment fault, be considered as this power node
Lost efficacy within the cascading failure time period；
Step 4.2：Connection matrix E in laststatecpChange, cause electrical network net interior nodes and the inefficacy on side, until electrical network
Reach stable state, status number adds 1, judge to connect E between netpcWhether changing, if changing, updating Epc, execute next step；If
Constant, cascading failure terminates, and jumps to step 4.4；
Step 4.3：Connection matrix E in laststatepcChange, cause communication network net interior nodes and the inefficacy on side, Zhi Daotong
Letter Wondertek arrives stable state, and status number adds 1, judges connection E between netcpWhether changing, if changing, updating Ecp, return to step
4.2；If constant, cascading failure terminates, and executes next step；
Step 4.4：Coupling network cascading failure terminates, and obtains maximum function group GNZP of grid nodes；
Step 5：Node at random or in calculated attack communication network, obtains cascading failure under different allocation rule and terminates rear electrical network section
Maximum function group GNZP of point, chooses optimal load redistribution method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410156420.9A CN103957032B (en) | 2014-04-17 | 2014-04-17 | Load redistribution method for electric power coupling network to resist cascade failure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410156420.9A CN103957032B (en) | 2014-04-17 | 2014-04-17 | Load redistribution method for electric power coupling network to resist cascade failure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103957032A CN103957032A (en) | 2014-07-30 |
| CN103957032B true CN103957032B (en) | 2017-02-22 |
Family
ID=51334271
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410156420.9A Expired - Fee Related CN103957032B (en) | 2014-04-17 | 2014-04-17 | Load redistribution method for electric power coupling network to resist cascade failure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103957032B (en) |
Families Citing this family (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104298593B (en) * | 2014-09-23 | 2017-04-26 | 北京航空航天大学 | SOA system reliability evaluation method based on complex network theory |
| CN105469317B (en) * | 2015-11-18 | 2019-10-15 | 中国电力科学研究院 | A kind of reliability of power communication network analysis method |
| CN105591386A (en) * | 2015-12-29 | 2016-05-18 | 中国电力科学研究院 | Method and apparatus for mitigating cascading failures of coupling network of power grid and communication network |
| CN105656198B (en) * | 2015-12-29 | 2019-02-15 | 中国电力科学研究院 | A kind of acquisition methods of powerline network redundant path strategy |
| CN105976206A (en) * | 2016-05-10 | 2016-09-28 | 东南大学 | Capacity pricing method of user access project |
| CN106339803B (en) * | 2016-08-24 | 2019-11-01 | 天津大学 | A kind of power distribution network methods of risk assessment considering communication failure |
| CN106371422B (en) * | 2016-08-31 | 2018-08-31 | 北京航空航天大学 | A method of prediction critical infrastructures fault propagation |
| CN106789376B (en) * | 2017-03-24 | 2020-09-22 | 大连大学 | Construction method of cascade failure model of command control network with hierarchical structure |
| CN107092984B (en) * | 2017-04-12 | 2020-10-09 | 北京航空航天大学 | Network function end node propagation prediction method based on cascade failure |
| CN107453919B (en) * | 2017-08-31 | 2021-01-15 | 中国人民解放军国防科技大学 | Complex network node importance evaluation method and system |
| CN107769962A (en) * | 2017-09-19 | 2018-03-06 | 贵州电网有限责任公司 | A kind of communication network failure cascade venture influence analysis method of attack resistance |
| CN107919985B (en) * | 2017-11-08 | 2020-10-27 | 南京邮电大学 | Application of complex network cascade failure capacity load architecture |
| CN109146173B (en) * | 2018-08-20 | 2022-03-01 | 中国民航大学 | Method for modeling and optimizing cascade failure survivability of airspace sector network |
| CN108847996A (en) * | 2018-08-30 | 2018-11-20 | 浙江工业大学 | A kind of delay restorative procedure of dynamic network |
| CN109919801B (en) * | 2019-03-15 | 2022-09-30 | 中山大学 | Coupling method and device based on node importance of power system |
| CN110224930B (en) * | 2019-05-30 | 2021-10-29 | 国网福建省电力有限公司 | Power communication network service recovery method based on logic cascade failure model |
| CN112597614A (en) * | 2019-09-17 | 2021-04-02 | 华为技术有限公司 | Method and device for generating physical network topological graph |
| CN111193629B (en) * | 2020-01-14 | 2022-02-01 | 西安电子科技大学 | Fault propagation method for dynamic load cascade failure of multilayer information network |
| CN112615365B (en) * | 2020-12-08 | 2022-04-19 | 国网四川省电力公司经济技术研究院 | Smart power grid vulnerability key point identification method and device |
| CN115242429B (en) * | 2022-06-09 | 2023-07-04 | 北京航空航天大学 | Command information-dependent network collapse method based on coupling weighting degree deviation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1901321A (en) * | 2006-07-20 | 2007-01-24 | 王甘雨 | Information communication network for power grid production and managing |
| CN102214920A (en) * | 2011-06-10 | 2011-10-12 | 华北电力大学 | Circuit cluster-based method for analyzing grid cascading failure based on |
| WO2012035592A1 (en) * | 2010-09-15 | 2012-03-22 | 株式会社 東芝 | Information communication device and information communication method |
| CN103050971A (en) * | 2013-01-30 | 2013-04-17 | 武汉大学 | Power grid cascading failure inhibition method based on operation limit matching model |
-
2014
- 2014-04-17 CN CN201410156420.9A patent/CN103957032B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1901321A (en) * | 2006-07-20 | 2007-01-24 | 王甘雨 | Information communication network for power grid production and managing |
| WO2012035592A1 (en) * | 2010-09-15 | 2012-03-22 | 株式会社 東芝 | Information communication device and information communication method |
| CN102214920A (en) * | 2011-06-10 | 2011-10-12 | 华北电力大学 | Circuit cluster-based method for analyzing grid cascading failure based on |
| CN103050971A (en) * | 2013-01-30 | 2013-04-17 | 武汉大学 | Power grid cascading failure inhibition method based on operation limit matching model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103957032A (en) | 2014-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103957032B (en) | Load redistribution method for electric power coupling network to resist cascade failure | |
| Li et al. | Multicast routing for decentralized control of cyber physical systems with an application in smart grid | |
| Lévesque et al. | Communications and power distribution network co-simulation for multidisciplinary smart grid experimentations | |
| CN102522747B (en) | Calculation method of power supply capability of power distribution system | |
| CN105896672B (en) | The charging method of mobile robot in a kind of wireless charging sensor network system | |
| CN109167349B (en) | Power CPS survivability quantitative evaluation method considering load optimization reconfiguration | |
| Wang et al. | Cascading risk assessment in power-communication interdependent networks | |
| CN105940425B (en) | Power control system, method and information transmission capability control system, method | |
| CN103298009B (en) | The modeling of complicated wireless network control system and fault propagation method | |
| CN109039766B (en) | Power CPS network risk propagation threshold determination method based on seepage probability | |
| CN105186578B (en) | There is the distributed automatic scheduling method of power system accurately calculating network loss ability | |
| CN112398946A (en) | Source network load storage scheduling control system and configuration method | |
| CN102904815B (en) | Scale-free network based router-level topology modeling method | |
| CN106786977B (en) | Charging scheduling method of electric vehicle charging station | |
| CN107241181A (en) | Grid power transmission circuit wireless monitor system | |
| CN108681786B (en) | Distributed power generation site selection planning method based on power supply community structure | |
| CN102983570A (en) | Islanding method for power distribution network comprising distributed generation and considering interconnection switch | |
| Zhang et al. | Distributed optimal active power dispatch under constraints for smart grids | |
| CN105633955B (en) | A kind of multifunctional solid current limiter Optimal Configuration Method | |
| CN101989743B (en) | Energy-saving power generation dispatching optimization method based on direct current power flow | |
| CN106101010A (en) | VBRAS flow load balance method and device based on end office's cloud | |
| CN107835130A (en) | A kind of flow allocation method and device | |
| CN107559182A (en) | A kind of water pump control method and control system | |
| CN104821604A (en) | Multi-photovoltaic generator group coordination control method based on event trigger mechanism | |
| CN103944748B (en) | Network-key-node self-similar-traffic generation simplification method based on genetic algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211021 Address after: 102206 East Building A207, xiaolizhuang Road South (No. 9, Beier Village), Beier village, Shahe Town, Changping District, Beijing Patentee after: Beijing huazeyuntong Network Technology Co.,Ltd. Address before: 102206 No. 2 Nong Road, Zhu Xin Zhuang, Beijing, Changping District Patentee before: NORTH CHINA ELECTRIC POWER University |
|
| 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: 20170222 |