CN106786529A - A kind of distributed static security analysis method - Google Patents

A kind of distributed static security analysis method Download PDF

Info

Publication number
CN106786529A
CN106786529A CN201611201901.2A CN201611201901A CN106786529A CN 106786529 A CN106786529 A CN 106786529A CN 201611201901 A CN201611201901 A CN 201611201901A CN 106786529 A CN106786529 A CN 106786529A
Authority
CN
China
Prior art keywords
calculation
distributed
fault set
node
analysis method
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.)
Granted
Application number
CN201611201901.2A
Other languages
Chinese (zh)
Other versions
CN106786529B (en
Inventor
郑惠萍
陈颖
刘新元
任正伟
曲莹
黄少伟
程雪婷
薄利明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Electric Power Research Institute of State Grid Shanxi Electric Power Co Ltd
Original Assignee
Tsinghua University
Electric Power Research Institute of State Grid Shanxi Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University, Electric Power Research Institute of State Grid Shanxi Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN201611201901.2A priority Critical patent/CN106786529B/en
Publication of CN106786529A publication Critical patent/CN106786529A/en
Application granted granted Critical
Publication of CN106786529B publication Critical patent/CN106786529B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001Methods 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, the step such as including the trends of distributed N 1 of Load flow calculation and dependent instruction initialization, distributed fault set screening and multithreading/process.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 and combined calculation identical result.The present invention reaches the same result of multiple regional power grid joint Load flow calculations using Distributed Power Flow method, and the joint static security analysis of interconnected electric power system are realized under existing scheduling method.Screened using distributed fault set in the present invention, the electrical distance based on line node is ranked up, and Distributed Power Flow calculating is carried out to the circuit in failure collection, calculating speed is improved with multithreading/process on the premise of precision of analysis is ensured.

Description

Distributed static security analysis method
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 capability assessment," in 2014IEEE PES general meeting | Conference exposure, 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θ10and 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 node, noteThen 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 (5)

1. A distributed static security analysis method, the 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.
2. The distributed static security analysis method according to claim 1, wherein the coordination side in step 1 refers to allocation of message headers to complete multithreading and setting of parameters related to the JFNG algorithm in the initialization process.
3. The distributed static security analysis method according to claim 1, wherein the calculation side in step 1 refers to that topology analysis and ground state power flow calculation are required during initialization.
4. 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.
5. 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.
CN201611201901.2A 2016-12-23 2016-12-23 A kind of distribution static security analysis method Active CN106786529B (en)

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 true CN106786529A (en) 2017-05-31
CN106786529B 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)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109861214A (en) * 2019-02-28 2019-06-07 中国电力科学研究院有限公司 Judge method, the system of regional power grid transient rotor angle stability weakness route
CN110970939A (en) * 2019-12-24 2020-04-07 华北电力大学 Distributed energy cluster optimization method and system
CN111262249A (en) * 2019-12-20 2020-06-09 国网新疆电力有限公司 Power grid flow simulation method and system
CN111478333A (en) * 2020-04-14 2020-07-31 广东电网有限责任公司广州供电局 Parallel static security analysis method for improving power distribution network recovery after disaster
CN117154852A (en) * 2023-10-27 2023-12-01 华南理工大学 Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint

Citations (3)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109861214A (en) * 2019-02-28 2019-06-07 中国电力科学研究院有限公司 Judge method, the system of regional power grid transient rotor angle stability weakness route
CN111262249A (en) * 2019-12-20 2020-06-09 国网新疆电力有限公司 Power grid flow simulation method and system
CN111262249B (en) * 2019-12-20 2024-01-30 国网新疆电力有限公司 Power grid tide simulation method and system
CN110970939A (en) * 2019-12-24 2020-04-07 华北电力大学 Distributed energy cluster optimization method and system
CN111478333A (en) * 2020-04-14 2020-07-31 广东电网有限责任公司广州供电局 Parallel static security analysis method for improving power distribution network recovery after disaster
CN111478333B (en) * 2020-04-14 2021-11-30 广东电网有限责任公司广州供电局 Parallel static security analysis method for improving power distribution network recovery after disaster
CN117154852A (en) * 2023-10-27 2023-12-01 华南理工大学 Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint
CN117154852B (en) * 2023-10-27 2024-01-02 华南理工大学 Power grid economic dispatching method based on alternate iteration under consideration of N-1 tide constraint

Also Published As

Publication number Publication date
CN106786529B (en) 2019-05-03

Similar Documents

Publication Publication Date Title
CN106786529B (en) A kind of distribution static security analysis method
Li et al. Dynamic equivalent modeling of two-staged photovoltaic power station clusters based on dynamic affinity propagation clustering algorithm
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
CN104090985A (en) Active disconnection optimum fracture surface searching method based on electrical distance
CN110854884B (en) Wind power collection area subsynchronous oscillation risk online assessment and early warning method
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
CN109755937A (en) A kind of regional power grid inertia calculation method and apparatus based on measurement
CN103324858A (en) Three-phase load flow state estimation method of power distribution network
Wu et al. Equivalent modeling method for regional decentralized photovoltaic clusters based on cluster analysis
CN106877309A (en) Stability judgment method for direct-current micro-grid in master-slave control mode
CN108649597B (en) Fault migration method and analysis method for influence of fault on HVDC commutation failure
Mao et al. Community-detection-based approach to distribution network partition
CN110120673A (en) Distributed transmission & distribution collaboration idle work optimization method and system based on Thevenin's equivalence parameter identification
CN116613746A (en) OS 2-based distributed power grid modeling model
CN104935249B (en) Photovoltaic power generation system stability verification method and apparatus
CN111834996B (en) Power grid line loss calculation method and device
CN110675276A (en) Method and system for inversion droop control of direct current power transmission system
KR102220716B1 (en) Apparatus for reconfiguring distribution line and method thereof
Sadeghi et al. Fully decentralized multi-agent coordination scheme in smart distribution restoration: Multilevel consensus
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