CN107769191A - key transmission cross-section automatic identifying method based on network topology and flow data - Google Patents
key transmission cross-section automatic identifying method based on network topology and flow data Download PDFInfo
- Publication number
- CN107769191A CN107769191A CN201610685173.0A CN201610685173A CN107769191A CN 107769191 A CN107769191 A CN 107769191A CN 201610685173 A CN201610685173 A CN 201610685173A CN 107769191 A CN107769191 A CN 107769191A
- Authority
- CN
- China
- Prior art keywords
- section
- transmission cross
- network
- cut
- power
- 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
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a kind of key transmission cross-section automatic identifying method based on network topology and flow data, it includes：Prepare electric network data, simplified model forms power network connected graph；Network contraction and the analysis of node branch road, and subnet division network are carried out to power network connected graph；Based on transmission cross-section key feature, initial transmission cross-section is generated by verification；Transmission cross-section importance index is calculated, and screens key transmission cross-section；Export important transmission cross-section composition and importance index.The problem of technical scheme provided by the invention is solved in traditional technical scheme, and artificial experience identification transmission cross-section is difficult in adapt to extensive and complicated and changeable practical power systems changes of operating modes.
Description
Technical field
The present invention relates to powernet simulation analysis field, and in particular to a kind of pass based on network topology and flow data
Key transmission cross-section automatic identifying method.
Background technology
With the expansion of interconnected network region and exchange capacity, especially big power outage is to society caused by interconnected network failure
Life almost causes disastrous effect, thus ensures the safety in operation of interconnected network, avoids great or even especially big power failure thing
Therefore become more and more important.There occurs several large area blackouts, loss are extremely serious in the world.Therefore, add
Strong grid safety monitoring, prevent that safe operation of the generation of massive blackout accident to ensureing power system is particularly important.Due to
Modern power network is in large scale and complicated, and it is unpractical to carry out monitor closely to each element, and urgent after failure
Constraints is too many during control, and carrying out load flow rectification to whole system often has little time, and power grid operation personnel need targetedly
Strengthen the key monitoring to system core link, wherein key transmission cross-section is one of key monitoring object.
Transmission cross-section is the key concept during practical power systems analysis calculates, and whole power network is divided into some areas accordingly
Domain manages and analysis.In power network actual motion, key transmission cross-section is carried out to monitor and control in real time, once break down, this
The principle for not expanding local fault influence, keeps the integrality and ability to transmit electricity of transmission cross-section as far as possible.Ensure transmission cross-section
Security can effectively suppress the development of power outage, spread and involve scope, be significant.
Traditional transmission cross-section is drawn by operation of power networks expert by off-line analysis, artificial selection.It is limited to specially
The limitation of family's experience and offline mode calculation scale, the section manually found generally only reflect power network under the extreme method of operation
Weak link, the adaptability changed to on-line operation mode is not strong.With the expansion of power network scale and the access of new energy, power network
The method of operation is increasingly changeable, notable in particular with the access of extensive intermittent new energy, the randomness of electric network swim mode
Increase, its security feature may continually change, such as：Non-key section originally is likely to become key sections, on the contrary
It is as the same.The method of traditional artificial discovery power network section can not adapt to the quick change of power system operating mode, may cause to close
The omission of key section, or even jeopardize power network safety operation.
Therefore, it is necessary to a kind of key transmission cross-section automatic identifying method based on network topology and flow data is developed,
The key transmission link of reflection complicated electric power system operation rule is picked out, is power grid security monitoring in real time and prevention cascading failure
Generation provides decision-making foundation.
The content of the invention
In order to overcome the above-mentioned deficiencies of the prior art, the present invention provides a kind of key based on network topology and flow data
Transmission cross-section automatic identifying method, it comprises the following steps：
Step 1：Electric network data is collected, builds power network connected graph；
Step 2：Network contraction and the analysis of node branch road, and subnet division network are carried out to the power network connected graph；
Step 3：Transmission cross-section feature is verified, generates initial transmission cross-section；
Step 4：Transmission cross-section importance index is calculated, and screens key transmission cross-section；
Step 5：Export key transmission cross-section composition and importance index.
Electric network data in step 1 includes：Network topology data, electric network swim data and power network geographical zone data.
Power network connected graph G (V, E) in step 1 is：Electric power networks containing summit and side；
Wherein：V={ v1,v2,…,vnIt is vertex set, represent the bus nodes in power system；
E={ e1,e2,…,enBe side set, represent the transmission line of electricity in power system；
Connected graph G (V, E) summit includes：Summit on branch summit and ring.
Step 2 comprises the following steps：
1) subregion where endpoint node i being retracted to coupled node j.
If 2) node j is located on branch, continue to be contracted in loop subregion.
If shown in node j node degree such as formula d (j), the individual connected node contractions of k (j) are shared to node j；
If d (j)=1+k (j), node j belongs to branch node, by node j and shrink its region node incorporate into
Subregion where non-constricted node connected node j；
The rest may be inferred, until be contracted to by whole node contractions on branch to should loop subregion in.
Subnet division network comprises the following steps in step 2：
1) using geographical zone as initial subnet network.
2) plant stand pre-processes；
3) each border plant stand is calculated respectively to the average electrical distance with all plant stands in its connected each sub-network
DijWith average electrical distance DijThe average electrical distance D of other plant stands into sub-network0, border plant stand is placed in its distance most
In near sub-network；
Average electrical distance
Wherein, B is the geographical zone that plant stand i is connected, nBIt is the plant stand number that geographical zone B is included, DijIt is factory
The electrical distance stood between i and plant stand j；
4) repeat step 2) and step 3), until the sub-network of all border plant stands in power network stops changing.
Electrical distance D between plant standij, be respectively from plant stand i and plant stand j to the equivalent impedance between two-port network,
Voltage after numerically equal to plant stand i input unitary currents between plant stand i and plant stand j；
According to principle of stacking, the electrical distance D to stand between i and plant stand jijIt is as follows：
In formula：ZiiFor plant stand i self-impedance, ZjjFor plant stand j self-impedance, ZijFor plant stand i and plant stand j mutual impedance；
If impedance matrix Z exponent number is n, then MijIt is a n dimensional vector, plant stand i and plant stand j take 1 and -1 respectively, remaining
Position takes 0.
Plant stand pretreatment includes；If the plant stand on branch is different from the plant stand subregion on its connected ring, make its with
The subregion of plant stand is identical on connected ring；If plant stand is all located at same subnet network on other connected rings of plant stand on the ring
In A, and itself it is not in the sub-network, then the plant stand is divided into sub-network A.
Step 3 includes：Essential characteristic based on transmission cross-section, the interconnection set of cut set between each sub-network is searched for, passed through
Verification generates initial transmission cross-section.
The verification of initial transmission cross-section includes：
According to the essential characteristic of transmission cross-section, whether the interconnection set verified between the sub-network of search meets trend side
Constrained to consistency constraint and electrical link, be initial transmission cross-section if meeting, it is on the contrary then be not belonging to initial transmission cross-section.
Direction of tide consistency constraint includes：
Make Pa=min { P→, P←, Pb=max { P→, P←, if Pa/Pb≤ ε, then it is judged to meeting trend consistency constraint,
And by PbDirection is as the section tidal current direction, then carries out electrical link constraint test；
Otherwise, trend consistency constraint is unsatisfactory for, the cut set section is not as transmission cross-section；
Wherein, P→And P←The active power value on 2 directions of cut set section is represented respectively, ε is previously given threshold value, its
Span is 0~0.15.
Electrical link constraint includes：
If min { Di,k}≥Ds min, then judge to meet that electrical link constrains；Otherwise, it is determined that to be unsatisfactory for electrical link about
Beam；
Wherein, Ds minFor previously given threshold value, Di,kTo cut-off branch road lkWith other branch roads of section liOpening between (i ≠ k)
The disconnected factor.
Search between each sub-network and form the interconnection set of cut set, comprise the following steps：
1) judge whether all branch roads form cut set in interconnection L, if so, then interconnection L is a transmission cross-section, search
Terminate；Conversely, disconnecting all branch roads in interconnection L, interconnection L is included into the interconnection set cut-off；
2) calculated using DC power flow algorithm：Cut-off in interconnection after contained branch road disconnection, branch road in remaining interconnection
Branch breaking distribution factor；
3) interconnection in branch breaking distribution factor belonging to maximum branch road is selected to cut-off interconnection in advance；
4) judge whether pre- to cut-off interconnection trend consistent with the line direction of tide that has cut-off contact：
If so, then go to step 5)；Conversely, judge to cut-off interconnection trend in advance and whether cut-off the ratio between interconnection trend
Less than given threshold value λ；
If so, then go to step 5)；Otherwise, the breadth coefficient that cut-offs of maximum branch road is set to 0, and the connection belonging to by the branch road
Winding thread is put back to do not cut-off the set of interconnection after, go to step 3).
5) judge that can the pre- interconnection that cut-offs form power network cut set with having cut-off interconnection；
If can, interconnection is cut-off in advance and has cut-off interconnection composition transmission cross-section, and section search terminates；
If it is not, the pre- interconnection that cut-offs is disconnected, the set for having cut-off interconnection, repeat step 2 are included into)~step
5), until pre- cut-off interconnection and cut-off interconnection composition power network cut set.
Branch breaking distribution factor DklAfter being cut-off for circuit l, power is transferred to the ratio on circuit k thereon, such as following formula institute
Show：
In formula：xkAnd xlFor circuit k and circuit l reactance；Xl-lAnd XklRespectively utilize the Impedance Moment of DC flow model
Mutual impedance of the self-impedance and circuit k both ends at the circuit l both ends that battle array obtains on circuit l both ends.
Step 4 includes：The importance index of transmission cross-section is calculated according to electric network swim data, is selected from initial transmission cross-section
Go out to reflect the key transmission cross-section of electrical network weak link.
Transmission cross-section s importance index Ds is shown below in power system：
Ys(i)For transmission cross-section s nargin,For the mean power nargin of all branch roads of transmission cross-section s, P 'sIt is disconnected to transmit electricity
The normalization of face s power.
As transmission cross-section nargin Ys(i)When larger, trend overload and system unstability do not occur for transmission cross-section；
As transmission cross-section nargin Ys(i)During close to zero, the importance for approaching limit transmitted power section drastically increases.
If the transmission capacity nargin of two transmission cross-sections is identical, the big transmission cross-section of transimission power is more important；
If the transimission power of two transmission sections is identical, the small transmission cross-section of transmission capacity nargin is more important.
Transmission cross-section s power PsNormalized calculation formula is as follows：
In formula, max { Pk}、min{PkBe respectively transmission cross-section power in power network maximum and minimum value.
Branch road i power margin is as follows in transmission cross-section s：
Ps(i)For the power of branch road i in transmission cross-section s, Ps(i)maxFor its maximum transmission power.
The mean power nargin of all branch roads of transmission cross-section sIt is as follows：
M is the quantity of transmission cross-section s branch roads.
Compared with immediate prior art, technical scheme provided by the invention has the advantages that：
1st, the present invention is solved in traditional technical scheme, and artificial experience identification transmission cross-section is difficult in adapt to extensive and multiple
The problem of miscellaneous changeable practical power systems changes of operating modes.
2nd, the present invention can improve the efficiency and reliability of key transmission cross-section automatic identification processing.
3rd, the present invention can carry out the quantitative evaluation of comprehensive and reasonable to transmission cross-section safety and stability significance level.4th, present invention knot
Close power network geographical zone and sub-network division searches out transmission cross-section, more meet the actual demand of operation of power networks control.
Brief description of the drawings
Fig. 1 is the calculation flow chart of the present invention.
Embodiment
Technical scheme is described in further details with reference to Figure of description.
To achieve the above object, technical scheme provided by the invention includes：
Step 1：Prepare electricity grid network topological data, electric network swim data, power network geographical zone data, by model simplification
Power network connected graph is formed afterwards；
Step 2：Network contraction is carried out to power network connected graph, calculates electrical distance between each plant stand, and combines geographical zone and believes
Whole power network is divided into several sub-networks by breath；
Step 3：Definition and essential characteristic based on transmission cross-section, search for the interconnection of composition cut set between each sub-network
Set, initial transmission cross-section is generated by verification；
Step 4：Transmission cross-section importance index is calculated according to electric network swim data, is filtered out from initial transmission cross-section anti-
Reflect the key transmission cross-section of electrical network weak link.
Prepare electricity grid network topological data, electric network swim data, power network geographical zone data in step 1, after simplification
Power network connected graph is formed, wherein, step 1 specifically includes following steps：
Step 1-1：Prepare electricity grid network topological data, electric network swim data, power network geographical zone data.
Grid topology data, cut-offs situation according to what is switched in electric power networks, and electric power is calculated by topological analysis algorithm
The real-time architecture topology of network, it is the basis of the advanced application such as tidal current analysis in electric network analysis, state estimation.
Electric network swim data, according to the operation bar of the elements such as given electric network composition, electrical network parameter, generator and load
Part, after Load flow calculation, determine power system each several part steady-state operating condition parameter.Generally, flow data includes：In system
Each power supply and the power of load point, pivot point voltage, the voltage of equalization point and phase angle, the voltage magnitude of each bus nodes and phase
Angle, the power distribution of each branch road, power attenuation of network etc..
Power network geographical zone, the characteristics of independent mutually and interconnection, is had according to the administrative power network of different Utilities Electric Co.s in power network,
The whole network hair power transmission and transforming equipment is divided into different regions according to geographical distribution and affiliated Utilities Electric Co., run for dispatching of power netwoks
Control and security and stability analysis.
Step 1-2：Power network connected graph is formed after model simplification.
Power network connected graph, i.e. electric power networks can represent with a connected graph G (V, E) containing summit and side, wherein：V
={ v1,v2,…,vn, vertex set is represented, represents the bus nodes in power system；E={ e1,e2,…,en, represent side
Set, represents the transmission line of electricity in power system.
Model simplification, i.e., main grid structure is not influenceed according to various kinds of equipment model and annexation in grid topology data, deletion
The equipment such as the capacitor of topological structure, reactor, small branch road, reduce the scale of connected graph.
In step 2, network contraction is carried out to power grid topological graph, calculates electrical distance between each plant stand, and combine geographical zone
Whole power network is divided into several sub-networks by information；
Step 2-1：Network contraction and the analysis of node branch road are carried out to power grid topological graph.
Network shrinks, and to the connected graph G (V, E) of electric power networks structure composition, will contact close plant stand and shrinks together,
The scale on electrical network figure G summits is reduced, improves the efficiency of sub-area division and section search.
Figure G summit, which can be divided into, forms 2 kinds of summit on branch summit and ring.For form branch summit, its with ring
Interconnection between summit is minimal cut set, can be retracted on ring and not considered further that.Node degree refers to and the node is associated
Side bar number, the also known as degree of association.
(1) node that node degree is 1 is branch road endpoint node, should be consistent with subregion where the node being uniquely connected, thus
Subregion where endpoint node i to be retracted to connected node j.
(2) if node j is still located on branch, need to continue to shrink.If node j node degree is d (j), the individual phases of k (j) are shared
Even node contraction is to node j, if d (j)=1+k (j), node j belongs to branch node, all is contracted to it by node j is related
The node in region incorporates subregion where the non-constricted node being connected with j into, and the rest may be inferred, until being contracted to d (j*) ＞ 1+k
(j*) the node j* of relation, i.e., by the subregion of loop where whole node contractions on branch to j.
Step 2-2：Electrical distance between each plant stand is calculated, and whole power network is divided into several with reference to geographical partition information
Sub-network；
Electrical distance between plant stand, it is the two-port network equivalent impedance entered in terms of plant stand i, j, is numerically equal to from factory
Voltage after the i that stands injection unitary currents between plant stand i and j.
The expression of impedance enabled node impedance matrix elements is waited according to principle of stacking, between plant stand, between plant stand i, j it is electric away from
From can be expressed as：
In formula：Zii, Zjj, ZijCorresponding element in respectively electric network impedance matrix Z, ZiiFor plant stand i self-impedance, ZjjFor
Plant stand j self-impedance, ZijFor plant stand i and plant stand j mutual impedance；If impedance matrix Z exponent number is n, then MijIt is a n dimension row
Vector, 1 and -1 is taken respectively in the i-th row and jth row, remaining position takes 0.
Sub-network division methods, refer to the definition according to transmission cross-section and feature, the composition branch road of same transmission cross-section should
With identical power supply area and loading zone, therefore, according to electrical distance between plant stand and geographical zone information, with power plant and power transformation
Stand and for unit into whole power network is divided to different sub-networks, each sub-network represents power supply area and load different in power network
Area.
Comprise the following steps that：
1) using geographical zone as initial sub-network.
2) plant stand is pre-processed, specifically included：For plant stand on branch, if its point with plant stand on connected ring
Area is different, then makes it identical with the subregion of plant stand on connected ring；For plant stand on ring, if plant stand on its other connected ring
It is all located in same subnet network A, and itself is not in the sub-network, then the plant stand is divided into sub-network A.
3) to each border plant stand, it is calculated respectively and arrives the electrical distance of all plant stands in connected each sub-network,
And calculate the average electrical distance D of its other plant stand into this sub-network0, the plant stand is placed in the sub-network closest with it
In.
In formula：B is the geographical zone that plant stand i is connected；nBIt is the plant stand number that geographical zone B is included；DijIt is factory
The electrical distance stood between i, j, SiFor the set of the sub-network of plant stand connection.
4) step 2 and step 3 are repeated, untill the sub-network of all border plant stands in power network no longer changes.
In step 3, definition and essential characteristic based on transmission cross-section, search between each sub-network and form the contact of cut set
Line set, initial transmission cross-section is generated by verification；
Transmission cross-section, refer to interrelated, one group of transmission line of electricity being closely connected on electricity net safety stable or electricity transaction
(can also be transformer or generator etc.).
Generally, transmission cross-section possesses 3 essential characteristics：First, section should be a minimal cut set of power network；Second, break
The effective power flow direction of branch road should be basically identical in face, for the active less relative to other branch roads of underloading branch road or transmission
Branch road allows trend reverse；3rd, the contact formed between a plurality of branch road of section is close, mutual branch breaking sensitivity
It is larger.
Cut set, refer to a power network connected graph G set of fingers, if these branch roads, which are removed, will make G be separated into two portions
Point, but if removing a branch road less, figure will be connection.
Initial transmission cross-section, refers to after power grid topological graph is divided into different sub-systems, searches between obtained each sub-network
Interconnector form transmission cross-section.
Initial transmission cross-section verification, i.e., according to the essential characteristic of transmission cross-section, verify and get in touch between the subsystem searched
Whether line set meets direction of tide consistency constraint and electrical link constraint, is initial transmission cross-section if verifying if, instead
Be then not belonging to initial transmission cross-section.
Constraint 1：Direction of tide consistency constraint.
Use P→And P←The active power value on 2 directions of cut set section is represented respectively.Make Pa=min { P→, P←, Pb=max
{P→, P←, if Pa/Pb≤ ε (ε is previously given threshold value, and its usual span is [0,0.15]), then judge to meet trend
Consistency constraint, and by PbDirection carries out electrical link constraint test again afterwards as the section tidal current direction；Otherwise, it is determined that
Trend consistency constraint is unsatisfactory for, the cut set section is not as transmission cross-section.
Constraint 2：Electrical link constrains.
Section to meeting direction of tide consistency constraint does further electrical link constraint checking.Specifically, choose
The most heavy branch road l of trend in cut set sectionk, as branch road is cut-off, calculate the branch road and other branch roads of section liBetween (i ≠ k)
Cut-off factor Di,kIf min { Di,k}≥Ds min(Ds minIt is electric 0.2), then to judge to meet for previously given threshold value, general value
Gas contact constraint；Otherwise, it is determined that being unsatisfactory for electrical link constraint.
The calculation formula of branch breaking distribution factor is as follows.
In formula：DklPower is transferred to the ratio on circuit k thereon after being cut-off for circuit l；xkAnd xlFor circuit k and l electricity
It is anti-；Xl-lAnd XklThe self-impedance at the circuit l both ends respectively obtained using the impedance matrix of DC flow model and circuit k both ends
Mutual impedance on circuit l both ends；DminFor threshold value.
Search between each sub-network and form the interconnection set of cut set, it is comprised the following steps that：
1) judge whether all branch roads form cut set in interconnection L, if so, then interconnection L is a transmission cross-section in itself,
Search terminates；If it is not, disconnecting all branch roads in interconnection L, interconnection L is included into and has cut-off interconnection set；
2) calculated and cut-off in interconnection after contained branch road disconnection using DC power flow algorithm, the branch of branch road in remaining interconnection
Cut-off distribution factor in road；
3) interconnection belonging to the maximum branch road of branch breaking distribution factor is selected to cut-off interconnection in advance；
4) it is whether consistent with having cut-off interconnection trend direction to cut-off interconnection trend in advance for judgement：
If so, then go to 5)；If it is not, judge to cut-off interconnection trend in advance and cut-off whether the ratio between interconnection trend is less than
Given threshold value λ；
If so, then go to 5)；If it is not, the breadth coefficient value of cut-offfing of the maximum branch road of branch breaking breadth coefficient is set to 0, will
Interconnection belonging to the branch road, which is put back to, does not cut-off interconnection set, goes to 3).
5) judge that can the pre- interconnection that cut-offs form power network cut set with having cut-off interconnection；
Deleted in G is schemed and cut-off interconnection in advance and cut-off the side of interconnection, and carry out topological analysis, judged after deleting
Whether figure G is divided into 2 independent subgraphs.If so, then form power network cut set；If it is not, power network cut set is not formed then.
If 6) power network cut set can be formed, interconnection is cut-off in advance and has cut-off interconnection composition transmission cross-section, section search
Terminate；If it is not, then disconnecting the pre- interconnection that cut-offs, it is included into and has cut-off interconnection set, repeat step 2)~step 5), until pre-
Cut-off interconnection and cut-off interconnection and may make up power network cut set.
In step 4, transmission cross-section importance index is calculated according to electric network swim data, is filtered out from initial transmission cross-section
Reflect the key transmission cross-section of electrical network weak link.
Key transmission cross-section, refer in power network real time execution the transmission of electricity that trend load is heavy, margin of safety is small, needs key monitoring
Section, this kind of section ratio in numerous transmission cross-sections is less, but occupies critical role.
Section importance, consider section conveying trend, the average margin of safety of section and section branch road minimum safe nargin, it is fixed
A kind of importance index of section safety and stability of justice, by calculating importance index and being ranked up so as to identify key transmission
Section.
In operation of power networks, key transmission cross-section is typically trend weight, the less transmission cross-section of margin of safety, it is necessary to carry out
Key monitoring and analysis.Therefore, section tidal current and margin of safety have been considered, defining section importance index is：
As section nargin Ys(i)When larger, section is not susceptible to trend overload and system unstability, nargin are important for section
The influence of degree is smaller, and now the importance of section is mainly by its transmission power PiDetermine；As section nargin Ys(i)Smaller closer zero
When, Ys(i)As section importance determinant, approaching the importance of the section of limit transmitted power will drastically increase.If two
The transmission capacity nargin of transmission cross-section is identical, then the big transmission cross-section of transimission power is more important.If the transmission of two transmission sections
Power is identical, then the small transmission cross-section of transmission capacity nargin is more important.The important of cross sections can be calculated by above-mentioned formula
Degree, sequence can obtain the key sections of system.
Each variable calculation formula of section importance is as follows：
If sharing n transmission cross-section in power system, wherein section s power is Ps, according to formula (5) by PsIt is normalized to：
In formula, max { Pk}、min{PkBe respectively grid power transmission section power maximum and minimum value.
If the power of i-th branch road is P in section ss(i), its maximum transmission power is Ps(i)max, the power margin of the branch road
For：
If the section shares m bar branch roads, then the mean power nargin of all branch roads of the section is
The margin of safety of section is not only relevant with the average load degree of section branch road, also maximum with loading level in section
Branch road it is closely related.Once certain branch road overloads in section, it is easy to causes the excision of whole section, and is born in section
The maximum branch road of load degree is the branch road most easily overloaded in section.Therefore, the margin of safety of section not only with each branch road of section
Mean power nargin is relevant, also the power margin minimum value min { Y with branch road in sections(i)Closely related.
According to the method described above, using certain actual electric network as Knowledge Verification Model, according to the powernet service data at a certain moment
And geographical zone, power network electrical partitioning is carried out, automatic search obtains initial transmission cross-section, and picks out weight according to importance index
Want transmission cross-section.The recognition result of preceding 3 important transmission cross-sections is listed in table 1.The importance index of each section is shown in Table 2.
1 important transmission cross-section of table
| Section title | Section power | Section power limit | Section importance |
| Section 1 | 2426 | 4446 | 1.789 |
| Section 2 | 1419 | 3388 | 1.123 |
| Section 3 | 692 | 1865 | 0.743 |
The section importance result of calculation of table 2
Finally it should be noted that：The above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, institute
The those of ordinary skill in category field with reference to above-described embodiment still can to the present invention embodiment modify or
Equivalent substitution, these are applying for this pending hair without departing from any modification of spirit and scope of the invention or equivalent substitution
Within bright claims.
Claims (22)
1. a kind of key transmission cross-section automatic identifying method based on network topology and flow data, it is characterised in that including such as
Lower step：
Step 1：Electric network data is collected, builds power network connected graph；
Step 2：Network contraction and the analysis of node branch road, and subnet division network are carried out to the power network connected graph；
Step 3：Transmission cross-section feature is verified, generates initial transmission cross-section；
Step 4：Transmission cross-section importance index is calculated, and screens key transmission cross-section；
Step 5：Export key transmission cross-section composition and importance index.
2. key transmission cross-section automatic identifying method according to claim 1, it is characterised in that the electricity in the step 1
Network data includes：Network topology data, electric network swim data and power network geographical zone data.
3. key transmission cross-section automatic identifying method according to claim 1, it is characterised in that the electricity in the step 1
Net connected graph G (V, E) is：Electric power networks containing summit and side；
Wherein：V={ v1,v2,…,vnIt is vertex set, represent the bus nodes in power system；
E={ e1,e2,…,enBe side set, represent the transmission line of electricity in power system.
4. key transmission cross-section automatic identifying method according to claim 3, it is characterised in that the connected graph G (V, E)
Summit include：Summit on branch summit and ring.
5. key transmission cross-section automatic identifying method according to claim 4, it is characterised in that the step 2 is included such as
Lower step：
1) subregion where endpoint node i being retracted to coupled node j.
If 2) the node j is located on branch, continue to be contracted in loop subregion.
6. key transmission cross-section automatic identifying method according to claim 5, it is characterised in that set the section of the node j
Shown in point degree such as formula d (j), the individual connected node contractions of k (j) are shared to the node j；
If d (j)=1+k (j), the node j belongs to branch node, and the node j and the node for being contracted to its region are drawn
Subregion where being grouped into the non-constricted node being connected with the node j；
The rest may be inferred, until be contracted to by whole node contractions on branch to should loop subregion in.
7. key transmission cross-section automatic identifying method according to claim 1, it is characterised in that divided in the step 2
Sub-network comprises the following steps：
1) using geographical zone as initial subnet network.
2) plant stand pre-processes；
3) each border plant stand is calculated respectively to the average electrical distance D with all plant stands in its connected each sub-networkijWith
The average electrical distance DijThe average electrical distance D of other plant stands into the sub-network0, by the border plant stand be placed in
In its closest sub-network；
The average electrical distance
Wherein, B is the geographical zone that plant stand i is connected, nBIt is the plant stand number that geographical zone B is included, DijBe plant stand i and
Electrical distance between plant stand j；
4) repeat step 2) and step 3), until the sub-network of all border plant stands in power network stops changing.
8. key transmission cross-section automatic identifying method according to claim 7, it is characterised in that electric between the plant stand
Distance Dij, it is respectively from plant stand i and plant stand j to the equivalent impedance between two-port network, it is single is numerically equal to plant stand i inputs
Voltage after the electric current of position between plant stand i and plant stand j；
According to principle of stacking, the electrical distance D to stand between i and plant stand jijIt is as follows：
In formula：ZiiFor plant stand i self-impedance, ZjjFor plant stand j self-impedance, ZijFor plant stand i and plant stand j mutual impedance；
If impedance matrix Z exponent number is n, then MijIt is a n dimensional vector, plant stand i and plant stand j take 1 and -1 respectively, remaining position
Take 0.
9. key transmission cross-section automatic identifying method according to claim 7, it is characterised in that the plant stand pretreatment bag
Include；
If the plant stand on branch is different from the plant stand subregion on its connected ring, its subregion with plant stand on connected ring is made
It is identical；If plant stand is all located in same subnet network A on other connected rings of plant stand on the ring, and itself is not on
In the sub-network, then the plant stand is divided into sub-network A.
10. key transmission cross-section automatic identifying method according to claim 1, it is characterised in that the step 3 includes：
Essential characteristic based on transmission cross-section, the interconnection set of cut set between each sub-network is searched for, it is disconnected to generate initial transmission of electricity through verification
Face.
11. key transmission cross-section automatic identifying method according to claim 10, it is characterised in that the initial transmission of electricity is disconnected
The verification in face includes：
According to the essential characteristic of the transmission cross-section, whether the interconnection set verified between the sub-network of search meets trend side
Constrained to consistency constraint and electrical link, be initial transmission cross-section if meeting, it is on the contrary then be not belonging to initial transmission cross-section.
12. key transmission cross-section automatic identifying method according to claim 11, it is characterised in that the direction of tide one
The constraint of cause property includes：
Make Pa=min { P→, P←, Pb=max { P→, P←, if Pa/Pb≤ ε, then it is judged to meeting trend consistency constraint, and will
PbDirection is as the section tidal current direction, then carries out electrical link constraint test；
Otherwise, trend consistency constraint is unsatisfactory for, the cut set section is not as transmission cross-section；
Wherein, P→And P←The active power value on 2 directions of cut set section, P are represented respectivelyaIt is the minimum of the active power value
Value, PbIt is the maximum of the active power value, ε is previously given threshold value, and its span is 0~0.15.
13. key transmission cross-section automatic identifying method according to claim 12, it is characterised in that the electrical link is about
Beam includes：
If min { Di,k}≥Ds min, then judge to meet that electrical link constrains；Otherwise, it is determined that to be unsatisfactory for electrical link constraint；
Wherein, Ds minFor previously given threshold value, Di,kTo cut-off branch road lkWith other branch roads of section liBetween (i ≠ k) cut-off because
Son.
14. key transmission cross-section automatic identifying method according to claim 11, it is characterised in that search for each sub-network
Between form cut set interconnection set, comprise the following steps：
1) judge whether all branch roads form cut set in interconnection L, if so, then the interconnection L is a transmission cross-section, search
Terminate；Conversely, disconnecting all branch roads in the interconnection L, the interconnection L is included into the interconnection set cut-off；
2) calculated using DC power flow algorithm：Cut-off in interconnection after contained branch road disconnection, the branch road of branch road in remaining interconnection
Cut-off distribution factor；
3) interconnection in the branch breaking distribution factor belonging to maximum branch road is selected to cut-off interconnection in advance；
4) judge described pre- to cut-off interconnection trend and whether the line direction of tide for having cut-off contact is consistent：
If so, then go to step 5)；Conversely, judge that the pre- interconnection trend of cut-offfing has cut-off the ratio between interconnection trend with described
Whether given threshold value λ is less than；
If so, then go to step 5)；Otherwise, the breadth coefficient that cut-offs of maximum branch road is set to 0, and the interconnection belonging to by the branch road
Put back to after not cut-offfing the set of interconnection, go to step 3).
5) judge that can the pre- interconnection and the interconnection that cut-off of cut-offfing form power network cut set；
If can, the pre- interconnection that cut-offs forms transmission cross-section with the interconnection that cut-off, and section search terminates；
If it is not, the pre- interconnection that cut-offs is disconnected, the set for having cut-off interconnection, repeat step 2 are included into)~step
5), until the pre- interconnection that cut-offs forms power network cut set with the interconnection that cut-off.
15. key transmission cross-section automatic identifying method according to claim 14, it is characterised in that the branch breaking point
Cloth factor DklAfter being cut-off for circuit l, power is transferred to the ratio on circuit k thereon, is shown below：
In formula：xkAnd xlFor circuit k and circuit l reactance；Xl-lAnd XklRespectively obtained using the impedance matrix of DC flow model
Mutual impedance of the self-impedance and circuit k both ends at the circuit l both ends arrived on circuit l both ends.
16. key transmission cross-section automatic identifying method according to claim 1, it is characterised in that the step 4 includes：
The importance index of transmission cross-section is calculated according to electric network swim data, reflection electrical network weak link is selected from initial transmission cross-section
Key transmission cross-section.
17. key transmission cross-section automatic identifying method according to claim 16, it is characterised in that transmitted electricity in power system
Section s importance index Ds is shown below：
Ys(i)For the nargin of the transmission cross-section s,For the mean power nargin of all branch roads of transmission cross-section s, P 'sFor institute
State the normalization of transmission cross-section s power.
18. key transmission cross-section automatic identifying method according to claim 17, it is characterised in that
As the transmission cross-section nargin Ys(i)When larger, trend overload and system unstability do not occur for the transmission cross-section；
As the transmission cross-section nargin Ys(i)During close to zero, the importance for approaching limit transmitted power section drastically increases.
19. key transmission cross-section automatic identifying method according to claim 17, it is characterised in that
If the transmission capacity nargin of two transmission cross-sections is identical, the big transmission cross-section of transimission power is more important；
If the transimission power of two transmission sections is identical, the small transmission cross-section of transmission capacity nargin is more important.
20. key transmission cross-section automatic identifying method according to claim 17, it is characterised in that
The power P of the transmission cross-section ssNormalized calculation formula is as follows：
In formula, max { Pk}、min{PkBe respectively transmission cross-section power in power network maximum and minimum value.
21. key transmission cross-section automatic identifying method according to claim 17, it is characterised in that the transmission cross-section s
Middle branch road i power margin is as follows：
Ps(i)For the power of branch road i in the transmission cross-section s, Ps(i)maxFor its maximum transmission power.
22. key transmission cross-section automatic identifying method according to claim 17, it is characterised in that the transmission cross-section s
The mean power nargin of all branch roadsIt is as follows：
M is the quantity of the transmission cross-section s branch roads.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610685173.0A CN107769191A (en) | 2016-08-18 | 2016-08-18 | key transmission cross-section automatic identifying method based on network topology and flow data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610685173.0A CN107769191A (en) | 2016-08-18 | 2016-08-18 | key transmission cross-section automatic identifying method based on network topology and flow data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107769191A true CN107769191A (en) | 2018-03-06 |
Family
ID=61261424
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610685173.0A Pending CN107769191A (en) | 2016-08-18 | 2016-08-18 | key transmission cross-section automatic identifying method based on network topology and flow data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107769191A (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108808661A (en) * | 2018-06-07 | 2018-11-13 | 国电南瑞科技股份有限公司 | Consider the multiple failure static security analysis method of scissors crossing and passway for transmitting electricity |
| CN108879693A (en) * | 2018-07-10 | 2018-11-23 | 福州大学 | A kind of passway for transmitting electricity method for fast searching based on principal component analysis |
| CN108923416A (en) * | 2018-07-07 | 2018-11-30 | 福州大学 | A kind of passway for transmitting electricity method for fast searching drawing rate dimensionality reduction to route based on node |
| CN108964069A (en) * | 2018-07-19 | 2018-12-07 | 中国电力科学研究院有限公司 | A kind of more round Load Shedding According To Frequency emulated computation methods of power grid and system |
| CN109638838A (en) * | 2019-01-21 | 2019-04-16 | 广东电网有限责任公司 | The recognition methods of power grid key sections, device and electronic equipment |
| CN110084513A (en) * | 2019-04-26 | 2019-08-02 | 中国电力科学研究院有限公司 | Method is determined based on the grid power transmission section of distribution factor |
| CN110175324A (en) * | 2019-04-11 | 2019-08-27 | 贵州电网有限责任公司 | A kind of operation of power networks operational order method of calibration and system based on data mining |
| WO2020078109A1 (en) * | 2018-10-17 | 2020-04-23 | 中国电力科学研究院有限公司 | Method, device, and storage medium for identifying weak section of electrical power grid |
| CN111130116A (en) * | 2019-12-14 | 2020-05-08 | 贵州电网有限责任公司 | Scheduling operation power flow checking method based on key topology change item identification |
| CN112186741A (en) * | 2020-09-10 | 2021-01-05 | 贵州电网有限责任公司 | Power grid operation weak point distinguishing method based on BPA (Business Process analysis) simulation calculation |
| CN112541100A (en) * | 2020-12-02 | 2021-03-23 | 国网安徽省电力有限公司 | Power grid weak link automatic searching method based on topology analysis |
| CN112989129A (en) * | 2021-03-18 | 2021-06-18 | 国家电网有限公司 | Power transmission network key perception section identification method and system based on graph theory |
| CN114665470A (en) * | 2022-03-21 | 2022-06-24 | 清华大学 | Method and device for identifying key power transmission section of power system |
| CN115882452A (en) * | 2023-01-04 | 2023-03-31 | 温州电力建设有限公司 | New energy consumption capability analysis and evaluation method considering source load uncertainty |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976842A (en) * | 2010-10-29 | 2011-02-16 | 清华大学 | Method for automatically acquiring key sections based on electrical partitioning |
| CN103500254A (en) * | 2013-10-08 | 2014-01-08 | 华北电力大学（保定） | Power grid key section automatic searching method with geographical zone boundary optimization taken into consideration |
-
2016
- 2016-08-18 CN CN201610685173.0A patent/CN107769191A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976842A (en) * | 2010-10-29 | 2011-02-16 | 清华大学 | Method for automatically acquiring key sections based on electrical partitioning |
| CN103500254A (en) * | 2013-10-08 | 2014-01-08 | 华北电力大学（保定） | Power grid key section automatic searching method with geographical zone boundary optimization taken into consideration |
Non-Patent Citations (1)
| Title |
|---|
| 徐岩等: "基于地理分区的关键输电断面识别", 《电气应用》 * |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108808661A (en) * | 2018-06-07 | 2018-11-13 | 国电南瑞科技股份有限公司 | Consider the multiple failure static security analysis method of scissors crossing and passway for transmitting electricity |
| CN108923416A (en) * | 2018-07-07 | 2018-11-30 | 福州大学 | A kind of passway for transmitting electricity method for fast searching drawing rate dimensionality reduction to route based on node |
| CN108879693A (en) * | 2018-07-10 | 2018-11-23 | 福州大学 | A kind of passway for transmitting electricity method for fast searching based on principal component analysis |
| CN108964069A (en) * | 2018-07-19 | 2018-12-07 | 中国电力科学研究院有限公司 | A kind of more round Load Shedding According To Frequency emulated computation methods of power grid and system |
| CN108964069B (en) * | 2018-07-19 | 2022-04-15 | 中国电力科学研究院有限公司 | Multi-turn low-frequency load reduction simulation calculation method and system for power grid |
| WO2020078109A1 (en) * | 2018-10-17 | 2020-04-23 | 中国电力科学研究院有限公司 | Method, device, and storage medium for identifying weak section of electrical power grid |
| CN109638838A (en) * | 2019-01-21 | 2019-04-16 | 广东电网有限责任公司 | The recognition methods of power grid key sections, device and electronic equipment |
| CN110175324A (en) * | 2019-04-11 | 2019-08-27 | 贵州电网有限责任公司 | A kind of operation of power networks operational order method of calibration and system based on data mining |
| CN110084513A (en) * | 2019-04-26 | 2019-08-02 | 中国电力科学研究院有限公司 | Method is determined based on the grid power transmission section of distribution factor |
| CN111130116A (en) * | 2019-12-14 | 2020-05-08 | 贵州电网有限责任公司 | Scheduling operation power flow checking method based on key topology change item identification |
| CN111130116B (en) * | 2019-12-14 | 2022-08-30 | 贵州电网有限责任公司 | Scheduling operation power flow checking method based on key topology change item identification |
| CN112186741A (en) * | 2020-09-10 | 2021-01-05 | 贵州电网有限责任公司 | Power grid operation weak point distinguishing method based on BPA (Business Process analysis) simulation calculation |
| CN112541100A (en) * | 2020-12-02 | 2021-03-23 | 国网安徽省电力有限公司 | Power grid weak link automatic searching method based on topology analysis |
| CN112541100B (en) * | 2020-12-02 | 2024-05-31 | 国网安徽省电力有限公司 | Automatic searching method for weak links of power grid based on topology analysis |
| CN112989129A (en) * | 2021-03-18 | 2021-06-18 | 国家电网有限公司 | Power transmission network key perception section identification method and system based on graph theory |
| CN112989129B (en) * | 2021-03-18 | 2024-05-24 | 国家电网有限公司 | Power transmission network key perception section identification method and system based on graph theory |
| CN114665470A (en) * | 2022-03-21 | 2022-06-24 | 清华大学 | Method and device for identifying key power transmission section of power system |
| CN115882452A (en) * | 2023-01-04 | 2023-03-31 | 温州电力建设有限公司 | New energy consumption capability analysis and evaluation method considering source load uncertainty |
| CN115882452B (en) * | 2023-01-04 | 2023-05-26 | 温州电力建设有限公司 | New energy consumption capability analysis and evaluation method considering uncertainty of source load |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107769191A (en) | key transmission cross-section automatic identifying method based on network topology and flow data | |
| Kashem et al. | Network reconfiguration for enhancement of voltage stability in distribution networks | |
| CN102868161B (en) | Optimization method of network variable structure with distributed type power supply distribution system | |
| Ganguly et al. | A novel multi-objective PSO for electrical distribution system planning incorporating distributed generation | |
| CN109274095A (en) | Low-voltage distribution network users topology estimation method and system based on mutual information | |
| CN105656040B (en) | A kind of power system restoration power computational methods for considering black starting-up process | |
| CN109086913A (en) | A kind of transient stability evaluation in power system method and system based on deep learning | |
| US20120022713A1 (en) | Power Flow Simulation System, Method and Device | |
| CN114124728B (en) | Vulnerability coupling link identification method between power grid control system and power primary equipment | |
| CN105656198B (en) | A kind of acquisition methods of powerline network redundant path strategy | |
| Akaber et al. | Automated post-failure service restoration in smart grid through network reconfiguration in the presence of energy storage systems | |
| CN113705085A (en) | Modeling and risk assessment method for multi-level structure of smart power grid | |
| Ramadan et al. | Optimal reconfiguration for vulnerable radial smart grids under uncertain operating conditions | |
| CN102290798A (en) | Quick power supply restoring rule based on distribution network wire connecting way | |
| CN105139081B (en) | A kind of distribution system contact point planing method of high reliability | |
| CN106253344B (en) | A kind of electric power networks and its control system and control method, network scheduling device | |
| CN117458494A (en) | Simplified modeling and analyzing method and system for busbar nodes of tide model | |
| CN108667027A (en) | Ac and dc systems power flow transfer search containing flexible direct current and quantitative analysis method | |
| CN114386222A (en) | Power distribution network cooperative fault recovery method and system based on information physical system | |
| CN107017635A (en) | A kind of adaptive islanding control method and system | |
| CN106899021B (en) | The load flow calculation system of electric system | |
| Aghamohammadi et al. | Controlled islanding based on the coherency of generators and minimum electrical distance | |
| Yu et al. | Optimization of an offshore oilfield multi-platform interconnected power system structure | |
| CN110768260A (en) | Power grid cascading failure model building method based on electrical betweenness | |
| CN115622041A (en) | Comprehensive performance evaluation method of power distribution network system based on complex network |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180306 |