CN109560547A - A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration - Google Patents

A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration Download PDF

Info

Publication number
CN109560547A
CN109560547A CN201910038858.XA CN201910038858A CN109560547A CN 109560547 A CN109560547 A CN 109560547A CN 201910038858 A CN201910038858 A CN 201910038858A CN 109560547 A CN109560547 A CN 109560547A
Authority
CN
China
Prior art keywords
distribution network
load
active
power
node
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.)
Pending
Application number
CN201910038858.XA
Other languages
Chinese (zh)
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 Power Grid Co Ltd
Zhuhai Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Guangdong Power Grid Co Ltd
Zhuhai Power Supply Bureau of Guangdong Power Grid Co Ltd
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 Power Grid Co Ltd, Zhuhai Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Guangdong Power Grid Co Ltd
Priority to CN201910038858.XA priority Critical patent/CN109560547A/en
Publication of CN109560547A publication Critical patent/CN109560547A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001Methods to deal with contingencies, e.g. abnormalities, faults or failures

Abstract

The invention discloses a kind of active distribution network N-1 safety evaluation method of consideration transmission & distribution collaboration, mainly comprises the steps that and establish active distribution network N-1 Dynamic Load Restoration Modeling, establish the transmission & distribution cooperated computing model based on principal and subordinate's splitting method and establish active distribution network N-1 index of security assessment system.The present invention is based on the methods of principal and subordinate's separation, it is contemplated that influence of the power transmission network for active distribution network, the index of security assessment after calculating active distribution network N-1 failure assess the safety of active distribution network from load level and voltage stabilization level is restored.This method can the safety accurately to power distribution network network structure and operation conditions evaluate, so as to avoid do not consider power transmission network influence and bring accuracy problem.

Description

A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration
Technical field
The present invention relates to distribution network technology fields, more particularly, to a kind of active distribution network N-1 of consideration transmission & distribution collaboration Safety evaluation method.
Background technique
Both at home and abroad is in for the main research with cooperation model the elementary step at present, mainly also focused at present therebetween The influence of operation and active distribution network have scholar to have studied active distribution network access the main function analysis for matching coupling network Voltage-controlled influence on power transmission network, but currently without method to active distribution network carry out security evaluation when consider transmission & distribution association Same influence, this will lead to missing of the security evaluation in accuracy.
The method of assessment security of distribution network does not consider after power distribution network N-1 failure fail-over policy to subtracting both at home and abroad at present The influence of few power loss loss, but after power distribution network N-1 failure, can generally exist recovery policy make to lose load area it is as small as possible with Guarantee maximum power supply, if not considering the load restoration strategy after power distribution network N-1 failure, does not meet actual conditions, will lead to peace Full assessment inaccuracy.
There are multi-voltage grades, a variety of methods of operation for main distribution net, huge for the main advanced computation amount with collaboration, Traditional tidal current computing method is unable to satisfy calculating demand.Consider main model and operating parameter, the data of real time data with fusion Amount, the multiplication of Load flow calculation operand;Meanwhile with the extensive access of distributed generation resource, cannot simply be modeled as PV node or PQ node, it is also necessary to consider the equation constraint that the control characteristic of distributed generation resource determines.Accordingly, it is considered to the calculating of main distribution collaboration The calculation for being more suitable for is needed to meet quick iteration demand.
Summary of the invention
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of active distribution network N- of consideration transmission & distribution collaboration 1 safety evaluation method.
In order to solve the above technical problems, the technical solution adopted by the present invention is that: a kind of active distribution considering transmission & distribution collaboration Net N-1 safety evaluation method, comprising the following steps:
S1: establishing active distribution network N-1 Dynamic Load Restoration Modeling, specifically includes step s11, s12 and step s13:
S11: establishing load restoration objective function, so that the minimum target the most of power loss load after power distribution network N-1 fault recovery Function:Wherein: LS,p,iFor active distribution network i-node burden with power, xS,ijFor the connection status of branch in model, 0 is disconnects, and 1 is closure, PS,ijFor the active power on power distribution network branch ij;
S12: the constraint condition of active distribution network N-1 Dynamic Load Restoration Modeling is established;
S13: Dynamic Load Restoration Modeling is improved.
S2: the transmission & distribution cooperated computing model based on principal and subordinate's splitting method is established, step s21 and s22 are specifically included:
S21: the global power flow algorithm of power transmission network and power distribution network is established;
S22: the global Load flow calculation process of transmission & distribution collaboration is determined.
S3: establishing active distribution network N-1 index of security assessment system, mainly includes step s31 and step s32;
S31: security evaluation coefficient is established;
S32: N-1 security evaluation is carried out to active distribution network.
Preferably, in the step s12, the constraint condition includes following constraint formula:
xS,ij∈{0,1},i,j∈φS,all (1)
Wherein, in formula: NS,node、NS,rootRespectively node total number, feeder line root node number;φS,outFor Fault Isolation The load bus set of power loss afterwards;φS,allFor available node sets all in network;J ∈ i is to have connection relationship with node i Node;Load at respectively failure front nodal point i is active and reactive;QS,ijFor on branch ij it is active with it is idle Trend numerical value;For the safe current limit value of branch ij;GS,ij、BS,ijFor the conductance and susceptance of branch ij;V S,iRespectively The voltage magnitude and its bound of node i;PS,iAnd QS,iIt is active and idle at i-node;θS,ijIt is between node i and j Phase angle difference,WithIt is the maximum value of the active and idle power output of distributed generation resource at i-node,WithIt is i section Distributed generation resource active and idle power output before failure at point, Δ PS,DG,iWith Δ QS,DG,iIt is i-node punishment cloth power failure The knots modification of the active and idle power output in front and back.
Constraint formula (1) represents the connection attribute of branch with 0-1 integer variable;Constraint formula (2) be closure branch number with Relationship between total node number;Constraint formula (3) (4) ensure that the non-power loss node in recovery front and back keeps original power constant;About The load that beam formula (5) (6) requires power loss node to restore is not less than 0 and is no more than failure preload original value;Constraint formula (2) and section Point power balance equation formula (3) to (6) the common guarantee connection of network and radial;Constraint formula (7) ensure that branch road Trend nonoverload;Constraint formula (8) to (11) is the relationship of voltage and power;Constraint formula (12) is node voltage amplitude constraint;About Beam formula (13) to (16) is node original loads and the relationship that distributed generation resource is contributed.
Preferably, in the step s13, due to constraint formula (8) to constraint formula (11) be it is non-convex, ask for convenience Constraint formula (8) to constraint formula (11) is converted to linear restriction (17) to linear restriction (20) by solution, as follows:
Pij=gij(Vi-Vj)-bijij) (19)
Qij=-bij(Vi-Vj)-gijij) (20)
Constraint formula (7) is converted into constraint formula (21),
For node voltage near per unit value, constraint formula (21) can be further modified to constraint formula (22)
Wherein, VS,0For the constant of voltage dimension.
Preferably, in the step s13, network-in-dialing when in order to guarantee to have zero injection node in a network is introduced Constraint formula (23) and constraint formula (24):
Improved Dynamic Load Restoration Modeling are as follows:
Constraint formula (1), constraint formula (2), constraint formula (12) to constraint formula (20), constraint formula (22) to constraint formula (24);
Preferably, in the step s21, power transmission network Load flow calculation equation is formula (29) are as follows:
Power transmission network Load flow calculation equation (30) are as follows:
Wherein, M,B,S represents power transmission network, marginal node, power distribution network, and k represents k-th of power distribution network;Equivalent nodal equation are as follows:
Wherein,It is power distribution network equivalence admittance matrix;Iteration auxiliary variableAre as follows:
Preferably, in the step s22, the global Load flow calculation process of transmission & distribution collaboration is mainly comprised the steps that
A: it obtainsIt is power distribution network equivalence admittance matrix;
B: initialization makes VB (k), k=0;
C: solve an equation (30) to equation (34), calculates SBS (k+1) and (SBS′)(k+1)
D: by (SBS′)(k+1)Solve an equation (29), obtains VB (k+1), k=k+1;
E: if max | VB (k+1)-VB (k)< ε, then terminate order;If not, then it is back to step C.
Preferably, in the step s31, three coefficient I are definedndex1、Index2And Index3, these three coefficients difference Indicate influence of influence, current operating condition of the network structure for load restoration for load restoration, current operating condition pair The influence of voltage level;
Wherein, the load total amount of power supply is lost after Fault Isolation, is usedIt indicates, referred to as total power loss load;It is total to lose The load total amount of power supply can not be connected in electric load by line loop operation, is usedIt indicates, it referred to as can not be extensive Multiple load;It may restore in load, since security constraint leads to the load total amount that cannot be sent a telegram in reply, useIt indicates, claims Not restore load;The load being finally recovered is denoted asIt is known as being resumed load;S-node indicates actively to match after restoring Whole number of nodes of power grid, VS,refIndicate reference voltage, αjIndicate the Risk rated ratio of node voltage deviation.
Preferably, the step s32 specifically includes the following steps:
A: initialization network structure enables i=1;
B: setting i-th line disconnects, and carries out N-1 scanning;
C: setting fault zone obtains that load may be restored;
D: load restoration optimization is carried out using improved Dynamic Load Restoration Modeling, is obtained
Lrestored, Loutage, LrestorableWith trend SBS (0), and calculation equation (31) is obtained to equation (34) (SBS′)(0)
E: with (SBS′)(0)For primary iteration value using transmission & distribution collaboration global Load flow calculation workflow management overall situation trend until Convergence obtains VS, j:
F: if i=N, the calculated result of each N-1 scanning is counted, calculates index of security assessment using formula (35);It is no Then, assignment i=i+1, and it is back to step b.
Compared with prior art, the beneficial effects of the present invention are:
The present invention is based on the methods of principal and subordinate's separation, it is contemplated that influence of the power transmission network for active distribution network calculates actively Index of security assessment after power distribution network N-1 failure, from load level and voltage stabilization level is restored, to the safety of active distribution network Property is assessed.N-1 fault scanning is carried out to active distribution network first, is then calculated using improved Dynamic Load Restoration Modeling extensive Multiple load cooperates with global tidal current computing method using transmission & distribution, calculates power distribution network after the active distribution network structure being restored Node voltage value, then Utilization assessment index system calculates index of security assessment, to complete entire security evaluation process.The party Method can the safety accurately to power distribution network network structure and operation conditions evaluate, so as to avoid power transmission network shadow is not considered It rings and bring accuracy problem.In addition, the present invention establishes active distribution network Dynamic Load Restoration Modeling, it may be considered that actively matching Fault recovery situation after power grid N-1 failure, the evaluation procedure more meet engineering reality.It is calculated to simplify, the present invention utilizes line Property method, linearizes AC power flow, avoids complicated calculations process, which can more preferably meet quickly Iteration demand.
Detailed description of the invention
Fig. 1 is flow chart of the invention;
Fig. 2 is the global Load flow calculation flow chart of transmission & distribution collaboration in the present invention;
Fig. 3 is the flow chart of active distribution network N-1 security evaluation in the present invention.
Specific embodiment
The present invention is further illustrated With reference to embodiment.Wherein, attached drawing only for illustration, What is indicated is only schematic diagram, rather than pictorial diagram, should not be understood as the limitation to this patent;Reality in order to better illustrate the present invention Example is applied, the certain components of attached drawing have omission, zoom in or out, and do not represent the size of actual product;To those skilled in the art For, the omitting of some known structures and their instructions in the attached drawings are understandable.
The same or similar label correspond to the same or similar components in the attached drawing of the embodiment of the present invention;It is retouched in of the invention In stating, it is to be understood that if the orientation or positional relationship for having the instructions such as term " on ", "lower", "left", "right" is based on attached drawing Shown in orientation or positional relationship, be merely for convenience of description of the present invention and simplification of the description, rather than indication or suggestion is signified Device or element must have a particular orientation, be constructed and operated in a specific orientation, therefore positional relationship is described in attached drawing Term only for illustration, should not be understood as the limitation to this patent, for the ordinary skill in the art, can To understand the concrete meaning of above-mentioned term as the case may be.
Embodiment
As shown in Figure 1, a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration, comprising the following steps:
S1: establishing active distribution network N-1 Dynamic Load Restoration Modeling, specifically includes step s11, s12 and step s13:s11: Load restoration objective function is established, so that the minimum objective function the most of power loss load after power distribution network N-1 fault recovery:Wherein: LS,p,iBurden with power for active distribution network in i-node, xS,ijFor The connection status of branch in model, 0 is disconnects, and 1 is closure, PS,ijFor the active power on power distribution network branch ij;
S12: the constraint condition of active distribution network N-1 Dynamic Load Restoration Modeling is established;
S13: Dynamic Load Restoration Modeling is improved.
S2: the transmission & distribution cooperated computing model based on principal and subordinate's splitting method is established, step s21 and s22 are specifically included:
S21: the global power flow algorithm of power transmission network and power distribution network is established;
S22: the global Load flow calculation process of transmission & distribution collaboration is determined.
S3: establishing active distribution network N-1 index of security assessment system, mainly includes step s31 and step s32;
S31: security evaluation coefficient is established;
S32: N-1 security evaluation is carried out to active distribution network.
Wherein, in step s12, constraint condition includes following constraint formula:
xS,ij∈{0,1},i,j∈φS,all (1)
Wherein, in formula: NS,node、NS,rootRespectively node total number, feeder line root node number;φS,outFor Fault Isolation The load bus set of power loss afterwards;φS,allFor available node sets all in network;J ∈ i is to have connection relationship with node i Node;Load at respectively failure front nodal point i is active and reactive;QS,ijFor on branch ij it is active with it is idle Trend numerical value;For the safe current limit value of branch ij;GS,ij、BS,ijFor the conductance and susceptance of branch ij;V S,iRespectively The voltage magnitude and its bound of node i;PS,iAnd QS,iIt is active and idle at i-node;θS,ijIt is between node i and j Phase angle difference,WithIt is the maximum value of the active and idle power output of distributed generation resource at i-node,WithIt is i section Distributed generation resource active and idle power output before failure at point, Δ PS,DG,iWith Δ QS,DG,iIt is i-node punishment cloth power failure The knots modification of the active and idle power output in front and back.
Constraint formula (1) represents the connection attribute of branch with 0-1 integer variable;Constraint formula (2) be closure branch number with Relationship between total node number;Constraint formula (3) (4) ensure that the non-power loss node in recovery front and back keeps original power constant;About The load that beam formula (5) (6) requires power loss node to restore is not less than 0 and is no more than failure preload original value;Constraint formula (2) and section Point power balance equation formula (3) to (6) the common guarantee connection of network and radial;Constraint formula (7) ensure that branch road Trend nonoverload;Constraint formula (8) to (11) is the relationship of voltage and power;Constraint formula (12) is node voltage amplitude constraint;About Beam formula (13) to (16) is node original loads and the relationship that distributed generation resource is contributed.
In addition, in step s13, it, will about since constraint formula (8) to constraint formula (11) is non-convex, solution for convenience Beam formula (8) to constraint formula (11) is converted to linear restriction (17) to linear restriction (20), as follows:
Pij=gij(Vi-Vj)-bijij) (19)
Qij=-bij(Vi-Vj)-gijij) (20)
Constraint formula (7) is converted into constraint formula (21),
For node voltage near per unit value, constraint formula (21) can be further modified to constraint formula (22)
Wherein, VS,0For the constant of voltage dimension.
Wherein, in step s13, network-in-dialing when in order to guarantee to have zero injection node in a network introduces constraint formula (23) and formula (24) are constrained:
Improved Dynamic Load Restoration Modeling are as follows:
Constraint formula (1), constraint formula (2), constraint formula (12) to constraint formula (20), constraint formula (22) to constraint formula (24);
Preferably, in the step s21, power transmission network Load flow calculation equation is formula (29) are as follows:
Power transmission network Load flow calculation equation (30) are as follows:
Wherein, M,B,S represents power transmission network, marginal node, power distribution network, and k represents k-th of power distribution network;Equivalent nodal equation are as follows:
Wherein,It is power distribution network equivalence admittance matrix;Iteration auxiliary variableAre as follows:
As shown in Fig. 2, the global Load flow calculation process of transmission & distribution collaboration mainly comprises the steps that A: obtaining in step s22 ?It is power distribution network equivalence admittance matrix;
B: initialization makes VB (k), k=0;
C: solve an equation (30) to equation (34), calculates SBS (k+1)And (SBS′)(k+1)
D: by (SBS′)(k+1)Solve an equation (29), obtains VB (k+1), k=k+1;
E: if max | VB (k+1)-VB (k)| < ε then terminates order;If not, then it is back to step C.
Wherein, in step s31, three coefficient I are definedndex1、Index2And Index3, these three coefficients respectively indicate net Influence of influence, current operating condition of the network structure for load restoration for load restoration, current operating condition are to voltage water Flat influence;
Wherein, the load total amount of power supply is lost after Fault Isolation, is usedIt indicates, referred to as total power loss load;It is total to lose The load total amount of power supply can not be connected in electric load by line loop operation, is usedIt indicates, it referred to as can not be extensive Multiple load;It may restore in load, since security constraint leads to the load total amount that cannot be sent a telegram in reply, useIt indicates, claims Not restore load;The load being finally recovered is denoted asIt is known as being resumed load;S-node indicates actively to match after restoring Whole number of nodes of power grid, VS,refIndicate reference voltage, αjIndicate the Risk rated ratio of node voltage deviation.
As shown in figure 3, step s32 specifically includes the following steps:
A: initialization network structure enables i=1;
B: setting i-th line disconnects, and carries out N-1 scanning;
C: setting fault zone obtains that load may be restored;
D: load restoration optimization is carried out using improved Dynamic Load Restoration Modeling, is obtained
Lrestored, Loutage, LrestorableWith trend SBS (0), and calculation equation (31) is obtained to equation (34) (SBS′)(0)
E: with (SBS′)(0)For primary iteration value using transmission & distribution collaboration global Load flow calculation workflow management overall situation trend until Convergence obtains VS, j:
F: if i=N, the calculated result of each N-1 scanning is counted, calculates index of security assessment using formula (35);It is no Then, assignment i=i+1, and it is back to step b.
In the present embodiment, N-1 fault scanning is carried out to active distribution network first, then using proposing in step s13 Improved Dynamic Load Restoration Modeling, which calculates, restores load, after the active distribution network structure being restored, using proposing in step s22 Transmission & distribution cooperate with global tidal current computing method, calculate the node voltage value of power distribution network at this time, then utilize the evaluation in step s31 Index system calculates index of security assessment, to complete entire security evaluation process.
Obviously, the above embodiment of the present invention be only to clearly illustrate example of the present invention, and not be pair The restriction of embodiments of the present invention.For those of ordinary skill in the art, may be used also on the basis of the above description To make other variations or changes in different ways.There is no necessity and possibility to exhaust all the enbodiments.It is all this Made any modifications, equivalent replacements, and improvements etc., should be included in the claims in the present invention within the spirit and principle of invention Protection scope within.

Claims (8)

1. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration, which comprises the following steps:
S1: establishing active distribution network N-1 Dynamic Load Restoration Modeling, specifically includes step s11, s12 and step s13:
S11: establishing load restoration objective function, so that the minimum objective function the most of power loss load after power distribution network N-1 fault recovery:Wherein: LS,p,iBurden with power for active distribution network in i-node, xS,ijFor The connection status of branch in model, 0 is disconnects, and 1 is closure, PS,ijFor the active power on power distribution network branch ij;
S12: the constraint condition of active distribution network N-1 Dynamic Load Restoration Modeling is established;
S13: Dynamic Load Restoration Modeling is improved.
S2: the transmission & distribution cooperated computing model based on principal and subordinate's splitting method is established, step s21 and s22 are specifically included:
S21: the global power flow algorithm of power transmission network and power distribution network is established;
S22: the global Load flow calculation process of transmission & distribution collaboration is determined.
S3: establishing active distribution network N-1 index of security assessment system, mainly includes step s31 and step s32;
S31: security evaluation coefficient is established;
S32: N-1 security evaluation is carried out to active distribution network.
2. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 1, feature exist In in the step s12, the constraint condition includes following constraint formula:
xS,ij∈{0,1},i,j∈φS,all (1)
Wherein, in formula: NS,node、NS,rootRespectively node total number, feeder line root node number;φS,outTo be lost after Fault Isolation The load bus set of electricity;φS,allFor available node sets all in network;J ∈ i is the section for having connection relationship with node i Point;Load at respectively failure front nodal point i is active and reactive;QS,ijFor the active and reactive power flow on branch ij Numerical value;For the safe current limit value of branch ij;GS,ij、BS,ijFor the conductance and susceptance of branch ij;V S,iRespectively node The voltage magnitude and its bound of i;PS,iAnd QS,iIt is active and idle at i-node;θS,ijIt is the phase angle between node i and j Difference,WithIt is the maximum value of the active and idle power output of distributed generation resource at i-node,WithIt is at i-node Distributed generation resource active and idle power output, Δ P before failureS,DG,iWith Δ QS,DG,iIt is i-node punishment cloth power failure front and back The knots modification of active and idle power output.
3. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 2, feature exist In, it is non-convex, solution for convenience due to constraining formula (8) to constraint formula (11) in the step s13, it will constraint formula (8) Linear restriction (17) is converted to linear restriction (20) to constraint formula (11), as follows:
Pij=gij(Vi-Vj)-bijij) (19)
Qij=-bij(Vi-Vj)-gijij) (20)
Constraint formula (7) is converted into constraint formula (21),
For node voltage near per unit value, constraint formula (21) can be further modified to constraint formula (22)
Wherein, VS,0For the constant of voltage dimension.
4. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 3, feature exist In, in the step s13, in order to guarantee in a network exist zero injection node when network-in-dialing, introduce constraint formula (23) with Constraint formula (24):
Improved Dynamic Load Restoration Modeling are as follows:
Constraint formula (1), constraint formula (2), constraint formula (12) to constraint formula (20), constraint formula (22) to constraint formula (24);
5. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 4, feature exist In in the step s21, power transmission network Load flow calculation equation is formula (29) are as follows:
Power transmission network Load flow calculation equation (30) are as follows:
Wherein, M,B,S represents power transmission network, marginal node, power distribution network, and k represents k-th of power distribution network;Equivalent nodal equation are as follows:
Wherein,It is power distribution network equivalence admittance matrix;Iteration auxiliary variableAre as follows:
6. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 5, feature exist In in the step s22, the global Load flow calculation process of transmission & distribution collaboration is mainly comprised the steps that
A: it obtainsIt is power distribution network equivalence admittance matrix;
B: initialization makes VB (k), k=0;
C: solve an equation (30) to equation (34), calculates SBS (k+1)And (SBS′)(k+1)
D: by (SBS′)(k+1)Solve an equation (29), obtains VB (k+1), k=k+1;
E: if max | VB (k+1)-VB (k)| < ε then terminates order;If not, then it is back to step C.
7. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 6, feature exist In defining three coefficient I in the step s31ndex1、Index2And Index3, these three coefficients respectively indicate network structure Influence for load restoration of influence, current operating condition for load restoration, current operating condition are to the shadow of voltage level It rings;
Wherein, the load total amount of power supply is lost after Fault Isolation, is usedIt indicates, referred to as total power loss load;Total power loss load In the load total amount of power supply can not be connected to by line loop operation, useIt indicates, can not referred to as restore negative Lotus;It may restore in load, since security constraint leads to the load total amount that cannot be sent a telegram in reply, useIt indicates, referred to as not Restore load;The load being finally recovered is denoted asIt is known as being resumed load;S-node indicates active distribution network after restoring Whole number of nodes, VS,refIndicate reference voltage, αjIndicate the Risk rated ratio of node voltage deviation.
8. a kind of active distribution network N-1 safety evaluation method for considering transmission & distribution collaboration according to claim 7, feature exist In, the step s32 specifically includes the following steps:
A: initialization network structure enables i=1;
B: setting i-th line disconnects, and carries out N-1 scanning;
C: setting fault zone obtains that load may be restored;
D: load restoration optimization is carried out using improved Dynamic Load Restoration Modeling, obtains Lrestorable, Loutage, Lestorable With trend SBS (0), and calculation equation (31) obtains (S to equation (34)BS′)(0)
E: with (SBS′)(0)For primary iteration value using the global Load flow calculation workflow management overall situation trend of transmission & distribution collaboration until receiving It holds back, obtains VS, j:
F: if i=N, the calculated result of each N-1 scanning is counted, calculates index of security assessment using formula (35);Otherwise, it assigns Value i=i+1, and it is back to step b.
CN201910038858.XA 2019-01-15 2019-01-15 A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration Pending CN109560547A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910038858.XA CN109560547A (en) 2019-01-15 2019-01-15 A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910038858.XA CN109560547A (en) 2019-01-15 2019-01-15 A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration

Publications (1)

Publication Number Publication Date
CN109560547A true CN109560547A (en) 2019-04-02

Family

ID=65873155

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910038858.XA Pending CN109560547A (en) 2019-01-15 2019-01-15 A kind of active distribution network N-1 safety evaluation method considering transmission & distribution collaboration

Country Status (1)

Country Link
CN (1) CN109560547A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111859585A (en) * 2020-07-06 2020-10-30 国电南瑞科技股份有限公司 Transmission and distribution cooperative equipment bearing capacity calculation method and system
CN112003277A (en) * 2020-08-21 2020-11-27 山东大学 Transmission and distribution cooperative load recovery optimization control method and system
CN113270871A (en) * 2020-02-17 2021-08-17 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 Flexible interconnection device capacity configuration optimization method based on power distribution network N-1 safety assessment
CN115347570A (en) * 2022-10-17 2022-11-15 国网浙江省电力有限公司宁波供电公司 Main-distribution cooperation-based regional power failure range analysis method
CN116307650A (en) * 2023-05-24 2023-06-23 东南大学溧阳研究院 Novel power distribution network source network load coordination random optimization operation method oriented to flexibility

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366315A (en) * 2013-07-24 2013-10-23 国家电网公司 Distribution network operating safety assessment method based on distribution network fault lost load recovery values
CN108599172A (en) * 2018-05-18 2018-09-28 广东电网有限责任公司佛山供电局 A kind of transmission & distribution net overall situation tidal current computing method based on artificial neural network

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103366315A (en) * 2013-07-24 2013-10-23 国家电网公司 Distribution network operating safety assessment method based on distribution network fault lost load recovery values
CN108599172A (en) * 2018-05-18 2018-09-28 广东电网有限责任公司佛山供电局 A kind of transmission & distribution net overall situation tidal current computing method based on artificial neural network

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孙宏斌等: "《发输配全局潮流计算——第一部分:数学模型和基本算法》", 《电网技术》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113270871A (en) * 2020-02-17 2021-08-17 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 Flexible interconnection device capacity configuration optimization method based on power distribution network N-1 safety assessment
CN113270871B (en) * 2020-02-17 2023-01-20 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 Flexible interconnection device capacity configuration optimization method based on power distribution network N-1 safety assessment
CN111859585A (en) * 2020-07-06 2020-10-30 国电南瑞科技股份有限公司 Transmission and distribution cooperative equipment bearing capacity calculation method and system
CN111859585B (en) * 2020-07-06 2022-07-29 国电南瑞科技股份有限公司 Transmission and distribution cooperative equipment bearing capacity calculation method and system
CN112003277A (en) * 2020-08-21 2020-11-27 山东大学 Transmission and distribution cooperative load recovery optimization control method and system
CN115347570A (en) * 2022-10-17 2022-11-15 国网浙江省电力有限公司宁波供电公司 Main-distribution cooperation-based regional power failure range analysis method
CN115347570B (en) * 2022-10-17 2023-01-24 国网浙江省电力有限公司宁波供电公司 Main-distribution cooperation-based regional power failure range analysis method
CN116307650A (en) * 2023-05-24 2023-06-23 东南大学溧阳研究院 Novel power distribution network source network load coordination random optimization operation method oriented to flexibility

Similar Documents

Publication Publication Date Title
Li et al. Double-mode energy management for multi-energy system via distributed dynamic event-triggered Newton-Raphson algorithm
Zhao et al. Coordinated restoration of transmission and distribution system using decentralized scheme
CN109560547A (en) A kind of active distribution network N-1 safety evaluation method considering transmission &amp; distribution collaboration
Roald et al. Chance-constrained AC optimal power flow: Reformulations and efficient algorithms
Hua et al. Voltage control for uncertain stochastic nonlinear system with application to energy Internet: Non-fragile robust H∞ approach
Khan et al. Distributed control algorithm for optimal reactive power control in power grids
Zhou et al. Distributed power management for networked AC–DC microgrids with unbalanced microgrids
Amini et al. Decomposition methods for distributed optimal power flow: panorama and case studies of the DC model
Mhanna et al. Component-based dual decomposition methods for the OPF problem
Zhao et al. Deep learning based model-free robust load restoration to enhance bulk system resilience with wind power penetration
Rezaei et al. Optimal IPFC signal selection and damping controller design using a novel current injection model in a multi-machine power system
Zamora-Cárdenas et al. State estimation of flexible AC transmission systems considering synchronized phasor measurements
Zhao et al. Optimal generator start-up sequence for bulk system restoration with active distribution networks
CN103887813A (en) Control method of wind power system operation based on wind power prediction uncertainty
Ding et al. Zonally robust decentralized optimization for global energy interconnection: Case study on northeast asian countries
Wang et al. A distributed proximal primal–dual algorithm for energy management with transmission losses in smart grid
Bie et al. Online multiperiod power dispatch with renewable uncertainty and storage: A two-parameter homotopy-enhanced methodology
Hu et al. Distributionally robust co-optimization of transmission network expansion planning and penetration level of renewable generation
Shchetinin et al. Decomposed algorithm for risk-constrained AC OPF with corrective control by series FACTS devices
Wu et al. Measurement-based online distributed optimization of networked distributed energy resources
Dong et al. A distributed power transfer limit calculation method for multi-area interconnection power networks
Han et al. Leader‐distributed follower‐decentralized control strategy for economic dispatch in cascaded‐parallel microgrids
Cavraro et al. Emergency voltage regulation in power systems via ripple-type control
Su et al. Interval TrendRank method for grid node importance assessment considering new energy
Mageshvaran et al. Steady state load shedding to prevent blackout in the power system using artificial bee colony algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20190402

RJ01 Rejection of invention patent application after publication