CN106786529B - A kind of distribution static security analysis method - Google Patents
A kind of distribution static security analysis method Download PDFInfo
- Publication number
- CN106786529B CN106786529B CN201611201901.2A CN201611201901A CN106786529B CN 106786529 B CN106786529 B CN 106786529B CN 201611201901 A CN201611201901 A CN 201611201901A CN 106786529 B CN106786529 B CN 106786529B
- Authority
- CN
- China
- Prior art keywords
- distributed
- node
- static security
- security analysis
- calculation
- 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.)
- Active
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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
Abstract
The present invention provides a kind of distributed static security analysis method, including Load flow calculation and dependent instruction initialization, the screening of distributed fault set and the distributed N-1 trend of multithreading/process and etc..On the basis of existing subdispatch operational mode, each subdispatch center calculates simultaneously, and a small amount of boundary information is exchanged by wide area network under the unified management of cooperation layer, to consider influencing each other between each regional power grid, and reach result identical with combined calculation.It is the same as a result, realizing the joint static security analysis of interconnected electric power system under existing scheduling method that the present invention reaches multiple regions power grid joint Load flow calculation using Distributed Power Flow method.It is screened in the present invention using distributed fault set, the electrical distance based on line node is ranked up, and carries out Distributed Power Flow calculating to the route in failure collection, improves calculating speed with multithreading/process under the premise of guaranteeing precision of analysis.
Description
Technical field
The invention belongs to the safety analysis field of electric system, especially a kind of distributed static security analysis method.
Background technique
Static security analysis is studied a certain by N-1 principle the network under running network or a certain research state
A operating element because failure it is out of service after, whether there is or not out-of-limit whether there is or not overload and busbar voltage for other elements in system.Traditional
In static security analysis, since each region control centre is mutually independent, and only whole moulds of compass of competency power grid
Type, therefore each subdispatch center carries out equivalent processing firstly the need of to external electrical network, could divide later local power grid
Analysis.At present in dispatching platform often by external network equivalent be static load or generator, in such a mode, interconnection region power grid
Joint static security analysis can not carry out, and influence of the failure in a regional power grid to another regional power grid can not be by
It embodies.
For relevant issues existing for external network equivalent method, in order to improve the accuracy of analysis, document Z.Haibo,
Z.Boming,S.Hongbin,and A.Ran,“A new distributed power flow algorithm between
multi-control-centers based on asynchronous iteration,”in 2006International
Conference on Power System Technology, 2006, pp.1-7. and document Z.w Liu and M.b Liu,
“Distributed Reactive Power Optimization Computing in Multi-area Power
Systems Using Ward Equivalent,”in 2010International Conference on Electrical
And Control Engineering (ICECE), 2010, pp.3659-3663. carries out Dai Weinan and Ward etc. to outer net respectively
Value, and Equivalent Model is modified in calculating process, more accurately to reflect external electrical network.Document Z.Li, J.Wang,
H.Sun,and Q.Guo,“Transmission Contingency Screening Considering Impacts of
Distribution Grids,”IEEE Trans.Power Syst.,vol.31,no.2,pp.1659–1660,Mar.2016.
Full trend is solved using principal and subordinate's algorithm, and by alternative iteration method, interior external power grid is involved in calculating in this method, and accuracy is more
Height, but principal and subordinate's convergence is low, the computational efficiency of alternating iteration is still to be improved.Alternative iteration method refers to multiple regions
Successively being calculated, after boundary information is issued another region after having been calculated by each region, another region can just be calculated,
Reciprocal with this, each moment can only have a region calculating.Document Y.Chen and C.Shen, " A Jacobian-
free Newton-GMRES(m)method with adaptive preconditioner and its application
for power flow calculations,”IEEE Trans.Power Syst.,vol.21,no.3,pp.1096–1103,
2006. improve the convergence of calculating using JFNG algorithm.Document Z.Ren et al., " Distributed power flow
considering network loss allocation and load factor of subareas,”in2016 35th
Chinese Control Conference (CCC), 2016, pp.2820-2824., which is established, considers via net loss interregional
The Distributed Power Flow model of distribution.
In addition, firstly the need of filtering out a failure collection (CCS) when carrying out N-1 static security analysis in electric system.
Document G.Zhou et al., " The static security analysis in power system based on
Spark Cloud Computing platform,”in Smart Grid Technologies-Asia(ISGT ASIA),
2015IEEE Innovative, 2015, pp.1-6. and document P.A.Kaplunovich and K.S.Turitsyn, " Fast
selection of N-2contingencies for online security assessment,”in 2013IEEE
Power Energy Society General Meeting, 2013, pp.1-5. is based on DC power flow, uses distribution factor method
Quickly scanning is carried out to filter out failure collection to system, and full Load flow calculation is carried out to the route in failure collection, this is mesh
Preceding common method.Document X.LIU, " Power System Dynamic Vulnerability under Extreme
Transmission Line Contingencies, " Master ' s Thesis, McGill University, make in 2007.
Failure collection is screened with eigenvalue sensitivity method, but is only used for generating the failure compared with large disturbances.Document C.A.Baone,
N.Acharya,S.Veda,and N.R.Chaudhuri,“Fast cont ingency screening and ranking
for small signal stability assessment,”in 2014IEEE PES General Meeting|
Conference Exposition, 2014, pp.1-5. pairs of eigenvalue sensitivity methods improve, sensitive using single order characteristic root
Degree method screens failure collection, but eigenvalue sensitivity method needs whole electric network models, therefore this method is current
It is difficult to realize in dispatching platform.
Existing dispatching platform only has the real time data of institute's compass of competency power grid, and carries out equivalence to external electrical network is unified
Processing.Most common method be PQ node type by boundary node equivalence, the injecting power of node remains permanent
It is fixed.And in actual electric network, the beginning of internal wiring will certainly impact the power on interconnection, and then power swing is passed
It is delivered to the external network of opposite end.Therefore there is deficiencies and power grid practical operation situation in accuracy for existing analysis method
And it mismatches.It is the same as a result, existing that the present invention reaches multiple regions power grid joint Load flow calculation using Distributed Power Flow method
The joint static security analysis of interconnected electric power system is realized under some scheduling methods.
Furthermore existing failure collection screening mode is to carry out line disconnection based on DC power flow, and using distribution factor method
Scanning.DC power flow only takes into account route active power, therefore the accuracy analyzed is not high.The present invention is in existing subdispatch
On the basis of operational mode, distributed static security analysis system is established, each subdispatch center calculates simultaneously, and is coordinating
Pass through wide area network (Wide Area Network, WAN) under the unified management of layer and exchange a small amount of boundary information, to consider each area
Influencing each other between the power grid of domain, and reach result identical with combined calculation.Combined calculation refers to each in interconnected electric power system
The unified calculation that regional power grid is carried out by complete data together with model integration and using the method for centralization.
Summary of the invention
The object of the present invention is to provide a kind of distributed static security analysis methods, are screened using distributed fault set,
Electrical distance based on line node is ranked up, and carries out Distributed Power Flow calculating to the route in failure collection, is being guaranteed
Calculating speed is improved with multithreading under the premise of precision of analysis.Each power grid simultaneously participates in analysis and calculates in this method,
Therefore there is more accurate input and output response relation compared to equivalent method, can accurately tracks the various change of external electrical network.Association
Adjust side use JFNG algorithm (full name Jacobi-Free Newton-GMRES, be it is a kind of with pretreatment mechanism based on wide
The Inexact Newton's method of adopted least residual can be used for solving Nonlinear System of Equations, not need explicitly to generate Jacobean matrix
Battle array.) alternative iteration method is compared with higher convergence.In addition, based on node electrical distance distributed fault collection screening with
And multithreading/process calculating mode can effectively improve the computational efficiency of static security analysis.
To achieve the above object, the present invention provides following schemes:
A kind of distribution static security analysis method, comprising: step 1: Load flow calculation and dependent instruction initialization refer to association
It adjusts side initialization N-1 static security analysis and sends a message to calculation side.Step 2: the screening of distributed fault set refers to each
Regional power grid can generate node admittance matrix in carrying out ground state Load flow calculation, and the inverse of node admittance matrix is node impedance square
Battle array, it reflects the size of electrical distance between each node, and each node sorts from small to large to boundary node electrical distance,
And route is grouped according to the node connected, it is sequentially placed into failure collection.Step 3: multithreading/process distribution
N-1 trend, calculation side is every carried out last round of calculating after, first trend section is restored, then from failure collection read under
One line information, and network topology structure is modified, next round calculating is carried out, when routes all in failure collection all have been calculated
Finish, or meets shutdown standard and then terminate calculating.
Optionally, the coordination side of step 1 refers to the distribution that the message header of multithreading is completed in initialization procedure and JFNG
The setting of algorithm relevant parameter.
Optionally, the calculation side of step 1 refers to and needs to carry out topological analysis and ground state Load flow calculation in initialization procedure.
Optionally, the shutdown standard of step 3 refers to and calculate and by out-of-limit situation by line to the route in failure collection
The grouping on road is sorted out and is compared, with the increase of electrical distance, if next group compared with upper one group it is not new out-of-limit
Information generates, then terminates calculating.
Optionally, the coordination side of step 3 can start multiple threads simultaneously and be calculated, and corresponding calculation side then needs out
Open multiple processes.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
Distributed static security analysis is proposed on the basis of existing dispatching platform, it can be between mutual regional power grid
Influence is analyzed and is assessed, to improve the accuracy of interconnected electric power system N-1 static security analysis.
Propose distributed failure collection screening technique.Each regional power grid is according to the electrical distance to boundary node to line
Road is ranked up and is grouped, and is put into failure collection, by using multithreading/process mode to the route in fault set simultaneously
Row analytical calculation is to improve computational efficiency.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is a kind of flow chart of distributed static security analysis method of the present invention
Fig. 2 is a kind of step 101 flow chart of distributed static security analysis method of the present invention；
Fig. 3 is a kind of step 102 flow chart of distributed static security analysis method of the present invention；
Fig. 4 is a kind of step 103 flow chart of distributed static security analysis method of the present invention；
Fig. 5 is a kind of subregion interconnected electric power system cutting figure of distributed static security analysis method of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of distributed static security analysis methods, are screened using distributed fault set,
Electrical distance based on line node is ranked up, and carries out Distributed Power Flow calculating to the route in failure collection, is being guaranteed
Calculating speed is improved with multithreading under the premise of precision of analysis.Each power grid simultaneously participates in analysis and calculates in this method,
Therefore there is more accurate input and output response relation compared to equivalent method, can accurately tracks the various change of external electrical network.Association
Adjust side use JFNG algorithm (full name Jacobi-Free Newton-GMRES, be it is a kind of with pretreatment mechanism based on wide
The Inexact Newton's method of adopted least residual can be used for solving Nonlinear System of Equations, not need explicitly to generate Jacobean matrix
Battle array.) alternative iteration method is compared with higher convergence.In addition, based on node electrical distance distributed fault collection screening with
And multithreading/process calculating mode can effectively improve the computational efficiency of static security analysis.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Embodiment 1:
It as shown in Figures 1 to 4, is a kind of flow chart of embodiment of the invention, a kind of distribution static security analysis method,
Method comprising steps of
101: Load flow calculation and dependent instruction initialization refer to and coordinate side initialization N-1 static security analysis and send to disappear
It ceases to calculation side.
102: the screening of distributed fault set refers to that each region power grid can generate node in carrying out ground state Load flow calculation
Admittance matrix, the inverse of node admittance matrix is nodal impedance matrix, it reflects the size of electrical distance between each node, will be each
A node sorts from small to large to boundary node electrical distance, and route is grouped according to the node connected, successively puts
Enter in failure collection.
103: the distributed N-1 trend of multithreading/process, calculation side is every carried out last round of calculating after, it is first disconnected to trend
Face is restored, then next line information is read from failure collection, and modify network topology structure, carries out next round meter
It calculates, is finished when routes all in failure collection all calculate, or meet shutdown standard and then terminate calculating.
The coordination side of step 101 refers to the distribution that the message header of multithreading is completed in initialization procedure and JFNG algorithm
The setting of relevant parameter.
The calculation side of step 101 refers to and needs to carry out topological analysis and ground state Load flow calculation in initialization procedure.
The shutdown standard of step 103 refers to and calculate and by out-of-limit situation by route to the route in failure collection
Grouping is sorted out and is compared, with the increase of electrical distance, if next group is compared not new more limit information with upper one group
It generates, then terminates calculating.
The coordination side of step 103 can start multiple threads simultaneously and be calculated, and corresponding calculation side then needs to open more
A process.
Embodiment 2:
By taking the interconnected network of two subregions as an example, as shown in figure 5, being classified as three child partitions, two of them are areas
Domain power grid, one is interconnection subregion.
Wherein,WithIndicate S1The boundary node injecting power of subregion,WithIndicate S1The boundary node of subregion
Voltage and phase angle.WithIt is then the S indicated2Subregion.Balance nodes are indicated using 0 in figure.V10And θ10It indicates
S1The voltage and angle values of partition balancing node, V20And θ20Then correspond to S2Subregion.P10And Q10Indicate S1The function of partition balancing node
Rate value, P20And Q20That indicate is then S2Subregion.It is ensured that there are a balance nodes in each subregion after subregion, can generally select
A PV node is selected as balance nodes, such as subregion S1Balance nodes, P10And V10It is known in advance, and θ10And Q10
It is then unknown；Should also there be a leading balancing machine in each subregion, voltage and phase angle, which are known, to be entire interconnection electricity
Force system provides voltage and phase angle benchmark.The output power P of generatorGAnd QGIt indicates, load uses PLAnd QLIt indicates.
The type of boundary node is generally handled as PQ node, noteIt is so right
In S1And S2Have:
To interconnection subregion, then haveAnd
Therefore:
For non-dominant balancing machine, the P by taking subregion 1 as an example10It is known that θ10It is unknown.Therefore for different θ10, can calculate
Obtain different active power
WithIndicate power flow equation.Therefore the equation of comptability are as follows:
When Distributed Power Flow calculating result is identical with joint calculation of tidal current,And Δ P10It should all wait
In 0.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said
It is bright to be merely used to help understand method and its core concept of the invention；At the same time, for those skilled in the art, foundation
Thought of the invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not
It is interpreted as limitation of the present invention.
Claims (3)
1. a kind of distribution static security analysis method, which is characterized in that the described method includes:
Step 1: Load flow calculation and dependent instruction initialization refer to and coordinate side initialization N-1 static security analysis and send message
To calculation side；The coordination side refers to that the distribution that the message header of multithreading is completed in initialization procedure is related to JFNG algorithm
The setting of parameter；The calculation side refers to and needs to carry out topological analysis and ground state Load flow calculation in initialization procedure；
Step 2: the screening of distributed fault set refers to that each region power grid can generate node in carrying out ground state Load flow calculation and lead
Receive matrix, the inverse of node admittance matrix is nodal impedance matrix, it reflects the size of electrical distance between each node, will be each
Node sorts from small to large to boundary node electrical distance, and route is grouped according to the node connected, is sequentially placed into
In failure collection；
Step 3: the distributed N-1 trend of multithreading/process, calculation side is every carried out last round of calculating after, first to trend section
Restored, then read next line information from failure collection, and modify network topology structure, carries out next round calculating,
It is finished when routes all in failure collection all calculate, or meets shutdown standard and then terminate calculating.
2. a kind of distributed static security analysis method according to claim 1, which is characterized in that the step 3 is stopped
Machine standard, refer to the route in failure collection carry out calculate and out-of-limit situation is sorted out and is compared by the grouping of route,
As the increase of electrical distance terminates calculating if next group is compared not new more limit information with upper one group and generated.
3. a kind of distributed static security analysis method according to claim 1, which is characterized in that the association of the step 3
It adjusts side that can start multiple threads simultaneously to be calculated, corresponding calculation side then needs to open multiple processes.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611201901.2A CN106786529B (en) | 2016-12-23 | 2016-12-23 | A kind of distribution static security analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611201901.2A CN106786529B (en) | 2016-12-23 | 2016-12-23 | A kind of distribution static security analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106786529A CN106786529A (en) | 2017-05-31 |
| CN106786529B true CN106786529B (en) | 2019-05-03 |
Family
ID=58897429
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611201901.2A Active CN106786529B (en) | 2016-12-23 | 2016-12-23 | A kind of distribution static security analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106786529B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109861214B (en) * | 2019-02-28 | 2022-10-04 | 中国电力科学研究院有限公司 | Method and system for judging weak line with stable transient power angle of regional power grid |
| CN111262249B (en) * | 2019-12-20 | 2024-01-30 | 国网新疆电力有限公司 | Power grid tide simulation method and system |
| CN110970939B (en) * | 2019-12-24 | 2021-07-30 | 华北电力大学 | Distributed energy cluster optimization method and system |
| CN111478333B (en) * | 2020-04-14 | 2021-11-30 | 广东电网有限责任公司广州供电局 | Parallel static security analysis method for improving power distribution network recovery after disaster |
| CN117154852B (en) * | 2023-10-27 | 2024-01-02 | 华南理工大学 | Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9099866B2 (en) * | 2009-09-01 | 2015-08-04 | Aden Seaman | Apparatus, methods and systems for parallel power flow calculation and power system simulation |
| CN105119268A (en) * | 2015-07-08 | 2015-12-02 | 国家电网公司 | Method for sorting transient severe faults of power system |
| CN106058856A (en) * | 2016-07-01 | 2016-10-26 | 中国南方电网有限责任公司 | Method for analyzing static security of grid rapidly |
-
2016
- 2016-12-23 CN CN201611201901.2A patent/CN106786529B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9099866B2 (en) * | 2009-09-01 | 2015-08-04 | Aden Seaman | Apparatus, methods and systems for parallel power flow calculation and power system simulation |
| CN105119268A (en) * | 2015-07-08 | 2015-12-02 | 国家电网公司 | Method for sorting transient severe faults of power system |
| CN106058856A (en) * | 2016-07-01 | 2016-10-26 | 中国南方电网有限责任公司 | Method for analyzing static security of grid rapidly |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106786529A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106786529B (en) | A kind of distribution static security analysis method | |
| Kroposki et al. | Autonomous energy grids: Controlling the future grid with large amounts of distributed energy resources | |
| Hu et al. | Improving frequency stability based on distributed control of multiple load aggregators | |
| CN106451556A (en) | Method and apparatus for determining capability of receiving distributed power supply by distributed network feeder line | |
| Lin et al. | Division algorithm and interconnection strategy of restoration subsystems based on complex network theory | |
| CN104716646B (en) | A kind of node Coupling Degrees method based on Injection Current | |
| CN106250640A (en) | A kind of layering Dynamic Equivalence being applicable to area power grid | |
| CN113098711A (en) | Power distribution Internet of things CPS (control system) management and control method and system based on cloud edge cooperation | |
| CN102253269A (en) | Integrated processing system for real-time data of electric power based on cloud computing and designing method | |
| CN103259289A (en) | Method for obtaining available transfer capability of wind power system with DFIG based on optimal power flow | |
| CN103838234A (en) | Generator wide area damping control system closed-loop testing method based on ADPSS | |
| CN102163845B (en) | Optimal configuration method of distributed generations (DG) based on power moment algorithm | |
| CN109193765A (en) | A kind of distributed photovoltaic cluster regulation method and device | |
| CN106877309B (en) | Stability judgment method for direct-current micro-grid in master-slave control mode | |
| Guzmán et al. | Microgrid power sharing framework for software defined networking and cybersecurity analysis | |
| CN109829599A (en) | The assemblage classification method and device of power distribution network based on high proportion renewable energy | |
| Cai et al. | A network model for the real-time communications of a smart grid prototype | |
| CN108649597A (en) | The analysis method that a kind of failure moving method and failure influence HVDC commutation failures | |
| CN104935249B (en) | Photovoltaic power generation system stability verification method and apparatus | |
| Chen et al. | Consensus control strategy with state predictor for virtual synchronous generators in isolated microgrid | |
| CN109818372A (en) | The control method and equipment of micro-capacitance sensor delay stability margin | |
| Lu et al. | Network simplification-based cluster coordinated optimization method for distributed PVs with inadequate measurement | |
| CN111834996B (en) | Power grid line loss calculation method and device | |
| CN109390943A (en) | Based on the online equivalent the whole network Analysis of Short-Circuit Current method and system of external network | |
| CN114243805A (en) | Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |