CN107480109B - Incomplete parallel state estimation realization method based on network structure dynamic partitioning - Google Patents
Incomplete parallel state estimation realization method based on network structure dynamic partitioning Download PDFInfo
- Publication number
- CN107480109B CN107480109B CN201710696694.0A CN201710696694A CN107480109B CN 107480109 B CN107480109 B CN 107480109B CN 201710696694 A CN201710696694 A CN 201710696694A CN 107480109 B CN107480109 B CN 107480109B
- Authority
- CN
- China
- Prior art keywords
- network
- state estimation
- voltage
- calculation
- voltage network
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
The invention relates to a state estimation incomplete parallel implementation method based on network structure dynamic partitioning, which comprises the following steps: high voltage network equivalence; medium and low voltage network equivalence; performing serial state estimation calculation on the high-voltage network; carrying out parallel state estimation and parallel computation on the medium-low voltage network; and correcting a calculation result. The state estimation incomplete parallel implementation method based on the network structure dynamic partition is provided based on the evolution of the power grid structure in the current extra-high voltage power grid construction process, is based on the equivalence of an original network, is not based on the coefficient matrix block formed after the topological analysis, and is simple to implement; because the scale of the high-voltage network is small, the serial calculation consumes little time and can be almost ignored; the calculation efficiency can be greatly improved without carrying out parallelization transformation on the original serial state estimation algorithm.
Description
Technical Field
The invention belongs to the field of electric power system dispatching automation, and provides a state estimation incomplete parallel implementation method based on network structure dynamic partitioning.
Background
The state estimation of the power system is an important component of an energy management system in a modern power grid, the scale of the power grid is larger and larger along with the expansion of the scale of the power grid and the change of a management mode, particularly along with the gradual construction of a national (branch) model, a provincial model and a data center, the model center of a provincial power grid comprises the power grid of more than 10kV, the number of nodes is tens of thousands, and the state estimation is used as a basic application of all power grid analysis to provide mature data for analysis and calculation for other applications. Since real-time control applications have high real-time requirements for mature data, high-performance state estimation methods need to be researched. At present, two technologies are mainly used for network partitioning, one technology is a branch cutting method, a current source is mainly added among sub-network connecting lines to replace the mutual influence of coupling among sub-network state quantities, the other technology is a node tearing method, a voltage source is mainly added after a network is divided into sub-networks to replace the coupling influence among the sub-network state quantities, a coefficient matrix reflecting the network characteristics of a power system is divided into a plurality of sub-matrixes, the sub-matrixes are coordinated and interconnected, and the sub-matrixes can be independently solved, so that the parallel connection is achieved. The branch cutting method and the node tearing method are based on coefficient matrix partitioning formed after topological analysis, are complex to realize, need to perform parallelization transformation on an original state estimation algorithm, and are low in calculation efficiency.
Disclosure of Invention
Aiming at the problems, the invention provides the incomplete parallel realization method for the state estimation based on the dynamic partition of the network structure, which can carry out dynamic partition according to the structural characteristics of a power grid, does not need to carry out parallel transformation on the original state estimation algorithm, and greatly improves the calculation efficiency by using a mode of combining serial calculation and parallel calculation.
In order to solve the problems, the invention adopts the technical scheme that: the incomplete parallel state estimation realizing method based on dynamic partitioning of the network structure comprises the following steps:
step 1, high-voltage network equivalence: the high-voltage network of 500kV and above is processed by taking a 500kV bus as a boundary, and the high-voltage side of a 500kV transformer is processed into an equivalent load, so that an independent high-voltage network isolated from a 220kV network is formed.
Step 2, medium and low voltage network equivalence: equating 500kV incoming lines and outgoing lines of a 500kV transformer substation into a generator; thereby forming a medium and low voltage network which is separated from the ring network of 500 kV.
And 3, performing serial state estimation calculation on the high-voltage network: and (3) calculating the high-voltage network formed based on the step (1) by utilizing an existing serial state estimation algorithm.
The scale of the high-voltage network formed in the step 1 is very small, and generally, more than 100 transformer substations of a provincial power grid of 500kV are not more than 100, so that the number of calculated nodes is not more than 100, the calculation time is very small and can be almost ignored, and the state estimation is performed on the high-voltage network.
Step 4, performing parallel state estimation and parallel computation on the medium-low voltage network: and (3) based on the plurality of medium and low voltage networks formed in the step (2), performing parallel calculation by using a multithreading technology and utilizing the existing serial state estimation algorithm, wherein the calculated voltage and the phase angle of the 500kV bus calculated in the step (3) are used as a reference voltage and a reference phase angle during calculation.
Step 5, correcting a calculation result: and (3) adding the voltage phase angle of each bus of each medium and low voltage network calculated in the step (4) with the voltage phase angle of the 500kV bus calculated in the step (3), and calculating the voltage reference nodes of the whole network.
Since the voltage reference point of each medium and low voltage network in the step 4 is the result calculated in the step 3, the voltage in the whole calculation result obtained in the step 4 is consistent with the result obtained by the conventional whole network calculation, and since the reference point of each node voltage phase angle in the step 4 is the 500kV bus, the voltage reference point of the whole network needs to be calculated according to the calculation result in the step 3.
The method mainly solves the problem of low state estimation calculation efficiency of the large-scale provincial power grid, provides timely and reliable calculation data for other high-level application software, firstly realizes the decoupling of calculation tasks by dynamically partitioning the provincial power grid, and then completes the calculation of the whole power grid by using a parallel calculation technology.
The state estimation incomplete parallel implementation method based on the network structure dynamic partition is provided based on the evolution of the power grid structure in the current extra-high voltage power grid construction process, is based on the equivalence of an original network, is not based on the coefficient matrix block formed after the topological analysis, and is simple to implement; because the scale of the high-voltage network is small, the serial calculation consumes little time and can be almost ignored; the calculation efficiency can be greatly improved without carrying out parallelization transformation on the original serial state estimation algorithm.
Drawings
Fig. 1 shows the grid structure characteristics of a provincial power grid to which the present invention is applied.
Fig. 2 shows the equivalent of the high voltage network of the present invention.
Fig. 3 is the equivalent of the medium and low voltage network of the present invention.
Detailed Description
Along with the construction of ultrahigh voltage power grids in China, high-voltage grid frames are increasingly strong, 500kV and 220kV power grids are gradually looped off, and the characteristics shown in figure 1 are presented, a 500kV network is operated as a looped network, a plurality of power supply areas are formed by the 220kV network and the 500kV network, the looped network is formed by 220kV substations in the same power supply area, but the looped network is not closed between the 220kV substations in each power supply area.
The incomplete parallel state estimation realizing method based on dynamic partitioning of the network structure comprises the following steps:
step 1, high-voltage network equivalence: as shown in fig. 2, the high voltage network of 500kV or more is processed with the 500kV bus as a boundary, and the high voltage side of the 500kV transformer is processed into an equivalent load, thereby forming an independent high voltage network isolated from the 220kV network.
Step 2, medium and low voltage network equivalence: as shown in fig. 3, the 500kV incoming line and outgoing line of the 500kV substation are equivalent to a generator; thereby forming a medium and low voltage network which is separated from the ring network of 500 kV.
And 3, performing serial state estimation calculation on the high-voltage network: and (3) calculating the high-voltage network formed based on the step (1) by utilizing an existing serial state estimation algorithm.
The scale of the high-voltage network formed in the step 1 is very small, and generally, more than 100 transformer substations of a provincial power grid of 500kV are not more than 100, so that the number of calculated nodes is not more than 100, the calculation time is very small and can be almost ignored, and the state estimation is performed on the high-voltage network.
Step 4, performing parallel state estimation and parallel computation on the medium-low voltage network: and (3) based on the plurality of medium and low voltage networks formed in the step (2), performing parallel calculation by using a multithreading technology and utilizing the existing serial state estimation algorithm, wherein the calculated voltage and the phase angle of the 500kV bus calculated in the step (3) are used as a reference voltage and a reference phase angle during calculation.
Step 5, correcting a calculation result: and (3) adding the voltage phase angle of each bus of each medium and low voltage network calculated in the step (4) with the voltage phase angle of the 500kV bus calculated in the step (3), and calculating the voltage reference nodes of the whole network.
Since the voltage reference point of each medium and low voltage network in the step 4 is the result calculated in the step 3, the voltage in the whole calculation result obtained in the step 4 is consistent with the result obtained by the conventional whole network calculation, and since the reference point of each node voltage phase angle in the step 4 is the 500kV bus, the voltage reference point of the whole network needs to be calculated according to the calculation result in the step 3.
Claims (1)
1. The method for realizing incomplete parallel state estimation based on dynamic partitioning of a network structure is characterized by comprising the following steps of based on the equivalence of an original network, not performing parallel transformation on an original state estimation algorithm, and using a mode of combining serial calculation and parallel calculation:
step 1, high-voltage network equivalence;
step 2, medium and low voltage network equivalence;
step 3, carrying out serial state estimation calculation on the high-voltage network;
step 4, performing parallel state estimation and parallel computation on the medium-low voltage network;
step 5, correcting a calculation result;
the high-voltage network equivalence method in the step 1 comprises the following steps: the 500kV bus is used as a boundary, a high-voltage network of 500kV and above is processed, and the high-voltage side of a 500kV transformer is processed into an equivalent load, so that an independent high-voltage network isolated from a 220kV network is formed;
the specific method of the medium-low voltage network equivalence in the step 2 is as follows: equating 500kV incoming lines and outgoing lines of a 500kV transformer substation into a generator; thereby forming a medium and low voltage network which is separated from the ring network of 500 kV;
calculating by using the existing serial state estimation algorithm based on the high-voltage network formed in the step 1;
based on the plurality of medium and low voltage networks formed in the step (2), performing parallel calculation by using a multithreading technology and utilizing an existing serial state estimation algorithm, wherein during calculation, the calculated voltage and phase angle of the 500kV bus calculated in the step (3) are used as a reference voltage and a reference phase angle;
the calculation result correction in the step 5 specifically comprises: and (3) adding the voltage phase angle of each bus of each medium and low voltage network calculated in the step (4) with the voltage phase angle of the 500kV bus calculated in the step (3), and calculating the voltage reference nodes of the whole network.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710696694.0A CN107480109B (en) | 2017-10-23 | 2017-10-23 | Incomplete parallel state estimation realization method based on network structure dynamic partitioning |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710696694.0A CN107480109B (en) | 2017-10-23 | 2017-10-23 | Incomplete parallel state estimation realization method based on network structure dynamic partitioning |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107480109A CN107480109A (en) | 2017-12-15 |
| CN107480109B true CN107480109B (en) | 2021-06-01 |
Family
ID=60599573
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710696694.0A Active CN107480109B (en) | 2017-10-23 | 2017-10-23 | Incomplete parallel state estimation realization method based on network structure dynamic partitioning |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107480109B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109390943B (en) * | 2018-12-27 | 2021-09-28 | 国电南瑞科技股份有限公司 | Whole-network short-circuit current analysis method and system based on external network online equivalence |
| CN112072668B (en) * | 2020-08-31 | 2022-07-26 | 华南理工大学 | Convex distributed optimal power flow solving method for power system |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101958580B (en) * | 2010-10-15 | 2015-01-28 | 上海申瑞继保电气有限公司 | Division calculation method for real time online trend flow of large power grid |
| US9281686B2 (en) * | 2012-11-06 | 2016-03-08 | General Electric Technology Gmbh | State estimation for cooperative electrical grids |
| CN103036230B (en) * | 2012-12-10 | 2016-03-23 | 上海市电力公司 | A kind of Dynamic Equivalence of the alternating current-direct current series-parallel connection bulk power grid based on engineer applied |
| CN103353921B (en) * | 2013-06-18 | 2016-08-31 | 国家电网公司 | One determines method based on heuristic automatic network subregion power system parallel flow |
| CN103595039B (en) * | 2013-10-17 | 2015-08-05 | 广东电网公司茂名供电局 | The distributed cyclization risk analysis method of network system |
| EP3127204A1 (en) * | 2014-04-04 | 2017-02-08 | Danmarks Tekniske Universitet | Thévenin equivalent based static contingency assessment |
| CN106250640B (en) * | 2016-08-04 | 2019-05-14 | 山东大学 | A kind of layering Dynamic Equivalence suitable for area power grid |
-
2017
- 2017-10-23 CN CN201710696694.0A patent/CN107480109B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107480109A (en) | 2017-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111525688A (en) | Electric power Internet of things power distribution network system based on 5G network and interconnection method | |
| CN104393617B (en) | A kind of topological structure of multiterminal flexible direct current power distribution network dynamic simulation experiment system | |
| CN107480109B (en) | Incomplete parallel state estimation realization method based on network structure dynamic partitioning | |
| CN102750416A (en) | Topological subnetting method of electromagnetic transient simulation containing switching characteristic circuit | |
| CN104979845A (en) | Multi-band high voltage direct current power transmission supplementary damping control method based on multi-stage linear optimal theory | |
| CN103474989A (en) | Network reconstruction method based on sensitivity analysis | |
| CN110120673B (en) | Distributed transmission and distribution cooperative reactive power optimization method and system based on Thevenin equivalent parameter identification | |
| CN112684290B (en) | Efficient calculation method for interelectrode short-circuit fault current of flexible direct-current power grid | |
| CN104200276B (en) | A kind of intelligent distribution network reconstructing method of feature based load injection | |
| CN106849154B (en) | Photovoltaic power station node harmonic voltage and harmonic source sensitivity relation analysis method | |
| Geth et al. | The PlanGridEV distribution grid simulation tool with EV models | |
| CN108471142B (en) | Distributed power grid frequency synchronization and active power distribution control method | |
| Meng et al. | Equivalent modeling and simulation for PV system on dynamic clustering equivalent strategy | |
| Yi et al. | An unified power flow method of hybrid AC/DC power system with VSC-MTDC | |
| CN103825453B (en) | A kind of structure changes multi input high-gain DC converter | |
| CN110380418B (en) | Newton-Raphson power flow calculation optimization method based on graph decomposition | |
| CN103368950B (en) | The method of communication between the coupled equipment of a kind of Electric Power Communication Data processor | |
| Tang et al. | Fpga-based real-time simulation for multiple energy storage systems | |
| CN111082522A (en) | Distribution network state monitoring and fault processing system based on ubiquitous power Internet of things | |
| Duan et al. | Research on flexible power distribution unit and its key technologies for energy internet | |
| Zhang et al. | Research on the construction mode of Ningxia distribution network Internet of things | |
| CN104917386B (en) | It is applicable the mode of resonance multi-input direct current transformer of new energy | |
| Sun et al. | Topology identification method based on multi-agent for outage areas in multi-contact distribution network | |
| Zhou et al. | Analysis of primary frequency regulation characteristics of PV power plant considering communication delay | |
| Di et al. | Research and application of low-voltage flexible interconnection technology of distribution 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |