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
- calculation
- distributed
- fault set
- node
- static security
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 44
- 230000003068 static effect Effects 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000012216 screening Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000005283 ground state Effects 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 abstract 1
- 238000005192 partition Methods 0.000 description 10
- 238000000638 solvent extraction Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000013097 stability assessment Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
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 field of safety analysis of a power system, and particularly relates to a distributed static safety analysis method.
Background
Static safety analysis is to study whether overload exists in other elements in a system and whether bus voltage exceeds the limit after a certain operation element quits operation due to faults according to the principle of N-1 for a network in operation or a network in a certain research state. In the traditional static security analysis, because each regional scheduling center is independent from each other and only has all models of power grids of the jurisdiction region, each regional scheduling center firstly needs to perform equivalence processing on an external power grid and then can analyze a local power grid. At present, an external grid is often equivalent to a static load or a generator in a dispatching platform, under the mode, joint static safety analysis of interconnected regional power grids cannot be performed, and the influence of a fault in one regional power grid on another regional power grid cannot be reflected.
In order to improve the accuracy of the analysis, the documents z.haibo, z.boming, s.hongbin, and a.ran, "a new Distributed Power flow algorithm based on asynchronous operation," in 2006International Conference on Power System Technology,2006, pp.1-7, Z.w Liu and M.b Liu, "Distributed Reactive Power Optimization Computing in Multi-area Power systems Using equation evaluation," in 2010International Conference on electric Equivalent Engineering (ice), 2010, 3659-3663, respectively, maintain the external network and correct the values in the south Equivalent process to reflect the accuracy of the external network. Documents z.li, j.wang, h.sun, and q.guo, "Transmission control mapping artifacts of distribution Grids," IEEE trans.power sys, vol.31, No.2, pp.1659-1660, and mar.2016. a master-slave algorithm is used, and a full power flow is solved by an alternating iteration method. The alternative iteration method is that a plurality of areas are calculated in sequence, after each area is calculated, boundary information is sent to another area, the other area can be calculated, and therefore the operation is repeated, and only one area is calculated at each moment. The document Y.Chen and C.Shen, "A Jacobian-free Newton-GMRES (m) method with adaptive preconditioners and bits application for power flow calculations," IEEE trans. Power Syst., vol.21, No.3, pp.1096-1103,2006. The documents Z.ren et al, "Distributed power flow configuration network loss allocation and load factor of subareas," in 201635 th Chinese Control reference (CCC),2016, pp.2820-2824. A Distributed power flow model is established that takes into account the distribution of network losses among the regions.
In addition, a fault set (CCS) needs to be screened first when performing N-1 static safety analysis in the power system. The 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), "2015 IEEE innovation, 2015, pp.1-6" and The document P.A.Kaplonvich and K.S.Tutstsingsyn, "faFaseelection of N-2 connectivity for online security assessment," in 2013IEEEPower Energy General Meeting,2013, pp.1-5. based on DC trend, The system is rapidly scanned using a distributed factor method to screen out The fault set and The lines in The fault set are calculated with full trend, which is a currently common method. The document X.LIU, "Power System Dynamic virtualization under extreme Transmission Line contacts," Master's Thesis, McGill University,2007, uses the feature root sensitivity method to screen the fault set, but only for faults that produce large perturbations. Documents c.a. baine, n.acharya, s.veda, and n.r. chaudhuri, "Fast contract entrance screening and ranking for small signal stability assessment," in 2014IEEE PES General Meeting | Conference exception, 2014, pp.1-5.
The existing dispatching platform only has real-time data of the power grid in the jurisdiction area, and the equivalence processing is uniformly carried out on the external power grid. The most common method at present is to equate the boundary node to PQ node type, and the injected power of the node is always constant. In an actual power grid, the open end of an internal line will inevitably affect the power on a tie line, and then the power fluctuation is transmitted to an external network at the opposite end. Therefore, the existing analysis method has the defects of accuracy and is not matched with the actual operation condition of the power grid. The invention uses a distributed power flow method to achieve the same result of combined power flow calculation of a plurality of regional power grids, and realizes the combined static safety analysis of the interconnected power system in the existing dispatching mode.
In addition, the existing fault set screening mode is based on direct current flow, and a distributed factor method is used for line disconnection scanning. The direct current power flow only considers the active power of the line, so the analysis accuracy is not high. On the basis of the existing regional dispatching operation mode, a distributed static safety analysis system is established, all regional dispatching centers perform simultaneous calculation, and a small amount of boundary information is exchanged through a Wide Area Network (WAN) under the unified management of a coordination layer, so that the mutual influence among regional power networks is considered, and the same result as the joint calculation is achieved. Joint calculation refers to unified calculation in which complete data and models are integrated by each regional power grid in an interconnected power system and a centralized method is used.
Disclosure of Invention
The invention aims to provide a distributed static safety analysis method, which uses distributed fault set screening, carries out sequencing based on the electrical distance of line nodes, carries out distributed load flow calculation on lines in a fault set and improves the calculation speed by multithreading on the premise of ensuring the analysis accuracy. In the method, each power grid simultaneously participates in analysis and calculation, so that the method has more accurate input and output response relation compared with an equivalence method, and can accurately track various changes of an external power grid. The JFNG algorithm (named Jacobi-Free Newton-GMRES in full, which is an inaccurate Newton method with a preprocessing mechanism and based on generalized minimum residual error, can be used for solving a nonlinear equation system and does not need to generate a Jacobian matrix explicitly) used by the coordination side has higher convergence compared with an alternative iteration method. In addition, the distributed fault set screening based on the node electrical distance and the calculation mode of multithreading/process can effectively improve the calculation efficiency of static safety analysis.
In order to achieve the purpose, the invention provides the following scheme:
a distributed static security analysis method, comprising: step 1: the load flow calculation and the related instruction initialization refer to that the coordination side initializes the N-1 static safety analysis and sends a message to the calculation side. Step 2: the distributed fault set screening means that each regional power grid generates a node admittance matrix in ground state load flow calculation, the inverse of the node admittance matrix is a node impedance matrix which reflects the magnitude of the electrical distance between each node, the electrical distances of each node to boundary nodes are sorted from small to large, and lines are grouped according to the connected nodes and are sequentially placed in a fault set. And step 3: and (3) recovering the load flow section after each previous calculation is finished at the calculation side, reading the information of the next line from the fault set, modifying the network topological structure, and performing the next calculation, wherein when all the lines in the fault set are calculated, or the calculation is terminated when the shutdown criterion is met.
Optionally, the coordination side in step 1 refers to the assignment of the message headers of multiple threads and the setting of the parameters related to the JFNG algorithm to be completed in the initialization process.
Optionally, the calculation side in step 1 refers to that topology analysis and ground state power flow calculation are required in the initialization process.
Optionally, the shutdown criterion of step 3 is to calculate the lines in the fault set, classify and compare the out-of-limit conditions according to the grouping of the lines, and terminate the calculation if no new out-of-limit information is generated when the next group is compared with the previous group as the electrical distance increases.
Optionally, the coordination side in step 3 may start multiple threads to perform calculation at the same time, and the corresponding calculation side needs to start multiple processes.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the distributed static safety analysis is provided on the basis of the existing dispatching platform, and the mutual influence among regional power grids can be analyzed and evaluated, so that the accuracy of the N-1 static safety analysis of the interconnected power system is improved.
A distributed fault set screening method is provided. And each regional power grid sorts and groups the lines according to the electrical distance of the boundary nodes, puts the lines into a fault set, and improves the calculation efficiency by using a multithreading/progress mode to perform parallel analysis and calculation on the lines in the fault set.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a distributed static security analysis method according to the present invention
FIG. 2 is a flowchart of a distributed static security analysis method of the present invention at step 101;
FIG. 3 is a flowchart of a step 102 of a distributed static security analysis method of the present invention;
FIG. 4 is a flow chart of step 103 of a distributed static security analysis method of the present invention;
fig. 5 is a partition diagram of a partitioned interconnected power system according to a distributed static security analysis method of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a distributed static safety analysis method, which uses distributed fault set screening, carries out sequencing based on the electrical distance of line nodes, carries out distributed load flow calculation on lines in a fault set and improves the calculation speed by multithreading on the premise of ensuring the analysis accuracy. In the method, each power grid simultaneously participates in analysis and calculation, so that the method has more accurate input and output response relation compared with an equivalence method, and can accurately track various changes of an external power grid. The JFNG algorithm (named Jacobi-Free Newton-GMRES in full, which is an inaccurate Newton method with a preprocessing mechanism and based on generalized minimum residual error, can be used for solving a nonlinear equation system and does not need to generate a Jacobian matrix explicitly) used by the coordination side has higher convergence compared with an alternative iteration method. In addition, the distributed fault set screening based on the node electrical distance and the calculation mode of multithreading/process can effectively improve the calculation efficiency of static safety analysis.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
Example 1:
as shown in fig. 1 to 4, which are flowcharts of an embodiment of the present invention, a distributed static security analysis method includes the steps of:
101: the load flow calculation and the related instruction initialization refer to that the coordination side initializes the N-1 static safety analysis and sends a message to the calculation side.
102: the distributed fault set screening means that each regional power grid generates a node admittance matrix in ground state load flow calculation, the inverse of the node admittance matrix is a node impedance matrix which reflects the magnitude of the electrical distance between each node, the electrical distances of each node to boundary nodes are sorted from small to large, and lines are grouped according to the connected nodes and are sequentially placed in a fault set.
103: and (3) recovering the load flow section after each previous calculation is finished at the calculation side, reading the information of the next line from the fault set, modifying the network topological structure, and performing the next calculation, wherein when all the lines in the fault set are calculated, or the calculation is terminated when the shutdown criterion is met.
The coordination side of step 101 refers to the assignment of the multi-threaded message header and the setting of the parameters related to the JFNG algorithm to be completed in the initialization process.
The calculation side of step 101 means that topology analysis and ground state power flow calculation are required in the initialization process.
The shutdown criteria of step 103 is to calculate the lines in the fault set and classify and compare the out-of-limit conditions according to the grouping of the lines, and with the increase of the electrical distance, if no new out-of-limit information is generated in the next group compared with the previous group, the calculation is terminated.
The coordinating side of step 103 may start multiple threads to perform computation at the same time, and the corresponding computing side needs to start multiple processes.
Example 2:
taking a two-partition interconnected network as an example, as shown in fig. 5, it is divided into three sub-partitions, two of which are regional networks and one of which is a tie-line partition.
Wherein,anddenotes S1The boundary nodes of the partitions inject power,anddenotes S1Voltage and phase angle of boundary nodes of the partition.Andthen is S of the representation2And (5) partitioning. The diagram uses 0 to represent the balance node. V10And theta10Denotes S1Voltage and phase angle values, V, of the zone balancing nodes20And theta20Then corresponds to S2And (5) partitioning. P10And Q10Denotes S1Power value of partition balancing node, P20And Q20Represents then S2And (5) partitioning. After partitioning, it is ensured that there is a balancing node in each partition, and typically a PV node may be selected as the balancing node, e.g. partition S1Of balancing node P10And V10Is known in advance, and theta10And Q10Is unknown; there should also be a lead balancer in each section, whose voltage and phase angle are known, to provide voltage and phase angle references for the entire interconnected power system. P for output power of generatorGAnd QGIndicating that the load is using PLAnd QLAnd (4) showing.
The type of boundary node is generally treated as PQ nodeMemory for recordingThen for S1And S2Comprises the following steps:
for the tie line division, then haveAndthus:
for the non-dominant balancing machine, take partition 1 as example P10Known as θ10Is unknown. Thus for different theta10Different active powers can be calculated
Andboth represent power flow equations. The coordination equation is therefore:
when the distributed power flow calculation result and the combined power flow calculation result are the same,and Δ P10Should all be equal to 0.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (3)
1. A distributed static security analysis method, the method comprising:
step 1: load flow calculation and related instruction initialization, namely, initializing N-1 static security analysis and sending a message to a calculation side by a coordination side; the coordination side is used for completing the allocation of a multi-thread message header and the setting of related parameters of the JFNG algorithm in the initialization process; the calculation side refers to that topology analysis and ground state load flow calculation are required in the initialization process;
step 2: distributed fault set screening means that each regional power grid generates a node admittance matrix in ground state load flow calculation, the inverse of the node admittance matrix is a node impedance matrix which reflects the magnitude of the electrical distance between each node, each node is ordered from small to large in terms of the electrical distance between boundary nodes, and circuits are grouped according to the connected nodes and are sequentially placed in a fault set;
and step 3: and (3) recovering the load flow section after each previous calculation is finished at the calculation side, reading the information of the next line from the fault set, modifying the network topological structure, and performing the next calculation, wherein when all the lines in the fault set are calculated, or the calculation is terminated when the shutdown criterion is met.
2. A distributed static safety analysis method according to claim 1, wherein the shutdown criteria of step 3 is to calculate the lines in the fault set and classify and compare the out-of-limit conditions according to the grouping of the lines, and as the electrical distance increases, if no new out-of-limit information is generated in the next group compared with the previous group, the calculation is terminated.
3. A distributed static security analysis method according to claim 1, wherein the coordination side in step 3 can start multiple threads to perform computation at the same time, and the corresponding computation side needs to start 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 | |
| CN103279638B (en) | Large power grid overall situation on-line integrated quantitative evaluation method based on response | |
| Lai et al. | Distributed power control for DERs based on networked multiagent systems with communication delays | |
| CN104794541B (en) | A kind of power distribution network running optimizatin method based on simulated annealing and cone optimization | |
| CN113098711A (en) | Power distribution Internet of things CPS (control system) management and control method and system based on cloud edge cooperation | |
| Zhao et al. | Graph-based Preconditioning Conjugate Gradient Algorithm for" N-1" Contingency Analysis | |
| CN105656039A (en) | Recognition method for vulnerable lines of power system | |
| CN105680442B (en) | Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure | |
| Zhang et al. | Adaptive and dynamic adjustment of fault detection cycles in cloud computing | |
| CN103324858A (en) | Three-phase load flow state estimation method of power distribution network | |
| CN106877309B (en) | Stability judgment method for direct-current micro-grid in master-slave control mode | |
| CN116613746A (en) | OS 2-based distributed power grid modeling model | |
| He et al. | Power system frequency situation prediction method based on transfer learning | |
| CN104935249B (en) | Photovoltaic power generation system stability verification method and apparatus | |
| CN111834996B (en) | Power grid line loss calculation method and device | |
| Hogan et al. | Towards effective clustering techniques for the analysis of electric power grids | |
| CN110675276A (en) | Method and system for inversion droop control of direct current power transmission system | |
| Sadeghi et al. | Fully decentralized multi-agent coordination scheme in smart distribution restoration: Multilevel consensus | |
| KR102220716B1 (en) | Apparatus for reconfiguring distribution line and method thereof | |
| Yesil et al. | Preliminary Studies on Dynamic Reduction of the Turkish Transmission Network | |
| CN112688332B (en) | Hybrid alternating current-direct current microgrid power flow analysis method and system based on sequence component coordinates | |
| Xue et al. | Effective and robust case screening for transient stability assessment | |
| CN109390943A (en) | Based on the online equivalent the whole network Analysis of Short-Circuit Current method and system of external network | |
| Li et al. | Parallel power flow calculation based on multi-port inversed matrix method | |
| Gao et al. | Fast emergency control strategy calculation based on dynamic equivalence and integral sensitivity |
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 |