CN106294022A - A kind of Jacobian matrix redundancy storage method for static security analysis - Google Patents
A kind of Jacobian matrix redundancy storage method for static security analysis Download PDFInfo
- Publication number
- CN106294022A CN106294022A CN201610662283.5A CN201610662283A CN106294022A CN 106294022 A CN106294022 A CN 106294022A CN 201610662283 A CN201610662283 A CN 201610662283A CN 106294022 A CN106294022 A CN 106294022A
- Authority
- CN
- China
- Prior art keywords
- jacobian matrix
- node
- electrical network
- redundancy
- matrix
- 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
Links
- 239000011159 matrix material Substances 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000003860 storage Methods 0.000 title claims abstract description 21
- 238000004458 analytical method Methods 0.000 title claims abstract description 15
- 230000003068 static effect Effects 0.000 title claims abstract description 15
- 230000005283 ground state Effects 0.000 claims abstract description 14
- 230000005611 electricity Effects 0.000 claims description 2
- 239000013589 supplement Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 7
- 239000013598 vector Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000005520 electrodynamics Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F11/00—Error detection; Error correction; Monitoring
- G06F11/07—Responding to the occurrence of a fault, e.g. fault tolerance
- G06F11/14—Error detection or correction of the data by redundancy in operation
- G06F11/1479—Generic software techniques for error detection or fault masking
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of Jacobian matrix redundancy storage method for static security analysis, the method comprises the steps: that (1), according to the admittance matrix Y of ground state electrical network before electrical network basic data generation electric network element fault, is stored by sparse format；(2) according to the sparsity structure of the admittance matrix Y of ground state electrical network, generate the Jacobian matrix J of sparse format, nodes all in electrical network are pressed PQ node processing, reserves storage area and voltage magnitude increment Delta V of non-zero entry for PV node redundancy in Jacobian matrix J2Storage area, i.e. form redundancy reactive power imbalance equation corresponding to PV node；(3) when generating the update equation group of any one concrete Fault load flow of ground state electrical network, Jacobian matrix J carries out decoupling process according to the actual node type of faulty grids.The present invention solves the problem that batch processing Load Flow Jacobian Matrix structure is inconsistent, provides data structure basis for batch processing power flow parallel computing.
Description
Technical field
The invention belongs to High performance computing in power system application, particularly relate to a kind of refined for static security analysis
Than matrix redundancy storage method.
Background technology
Load flow calculation is most widely used general, the most basic and most important electric computing of one in power system.In power train
In the research of the system method of operation and programme, it is required for carrying out Load flow calculation to compare the method for operation or to plan power supply plan
Feasibility, reliability and economy.And in actual production process, off-line trend and online power flow calculate all meters to trend
Calculate speed and have the requirement that this is higher.
The sparse format of Jacobian matrix is relevant to admittance battle array, and the most relevant to node type, node type does not changes
Becoming, Jacobian matrix dimension is constant, and sparse format is constant；If site type definition is different, even if being generated by same electrical network
The exponent number of Jacobian matrix and sparse format also differ.Fault load flow in static security analysis is to send out on ground state electrical network
Raw some broken strings produce, and any one N-1 disconnection fault at most affects 4 elements of ground state grid nodes admittance battle array, wherein
Nondiagonal element (transadmittance) may be become null element by non-zero entry, say, that the sparse format of any one Fault load flow is all base
One subset (refering in particular to non-zero entry quantity and form) of state trend sparse format, N-2 and other more complicated static security analysis
Also it is similar to.This guarantees admittance battle array during Load flow calculation, when using batch processing trend method, quote the one of admittance battle array data
Cause property.
Batch processing band wagon technique has strict requirements to data form, needs the data of the process of each Load flow calculation
Keep structure consistent.Newton-Laphson method needs to calculate Jacobian matrix, and Jacobian matrix is sparse vectors simultaneously
Coefficient matrix, needs also exist for ensureing the sparse format concordance of Jacobian matrix.
Summary of the invention
Goal of the invention: for the deficiencies in the prior art, the invention provides a kind of Jacobi for static security analysis
Matrix redundancy storage method, this method solves structure inconsistence problems in static security analysis Fault load flow Jacobian matrix,
The merging that can realize GPU memory access accesses, and memory access efficiency is greatly improved, and accelerates the calculating speed of the static security analysis that GPU accelerates.
Load flow calculation: electrodynamic noun, refers at given power system network topology, component parameters and generating, load parameter
Under the conditions of, calculate active power, reactive power and voltage distribution in power network.
Admittance matrix: that set up based on the Equivalent admittance of system element, describe each node voltage of electric power networks and
The matrix of relation between injection current.
A kind of Jacobian matrix redundancy storage method for static security analysis of the present invention, described method includes walking as follows
Rapid:
(1) according to the admittance matrix Y of ground state electrical network before electrical network basic data generation electric network element fault, deposit by sparse format
Storage；
(2) according to the sparsity structure of the admittance matrix Y of ground state electrical network, generate the Jacobian matrix J of sparse format, generating
During, nodes all in electrical network are pressed PQ node processing, redundancy opens up Jacobian matrix J, for PV node at Jacobian matrix J
Middle redundancy reserves storage area and voltage magnitude increment Delta V of non-zero entry2Storage area, i.e. form redundancy corresponding to PV node
Reactive power imbalance equation；
(3) when generating the update equation group of any one concrete Fault load flow of ground state electrical network, Jacobian matrix J according to
The actual node type of faulty grids carries out decoupling process: the idle amount of unbalance Δ Q supplemented by PV node redundancy2It is set to 0, will
All PV node redundancies supplement Jacobian matrix non-zero entry in element be set to 0, by redundancy corresponding for each PV node without
The element in Jacobian matrix corresponding with voltage magnitude Delta in merit unbalanced power equation is set to 1.
Wherein, described in described step (1), electrical network master data includes node serial number, number of nodes, branch number, branch road
Quantity, branch road reactance, the connected node serial number of branch road and parallel branch reactance.
First beneficial effect: compared with prior art, present invention have the advantage that, in the present invention, memory space is by all
Node is all that PQ node redundancy is reserved, and has a unified sparse format, thus the data structure during each Load flow calculation
Keep consistent, be suitable for the SIMT framework hardware such as GPU and data are merged the demand visited, realize the merging of GPU memory access on this basis
Access, memory access efficiency is greatly improved, the calculating speed of the static security analysis that GPU accelerates can be accelerated；Secondly concrete event is being generated
Process according to electrical network actual node type when hindering the Jacobian matrix of trend, it is ensured that the redundant data of interpolation does not interferes with former
Solving of update equation group；The last present invention solves the problem that batch processing Load Flow Jacobian Matrix structure is inconsistent, for static state
In safety analysis, batch Fault load flow parallel computation provides data structure basis.
Accompanying drawing explanation
Fig. 1 is the schematic diagram of the redundancy storage method of Jacobian matrix of the present invention.
Detailed description of the invention
Below in conjunction with the accompanying drawings technical scheme is described further.
As it is shown in figure 1, a kind of Jacobian matrix redundancy storage method for static security analysis of the present invention, described method
Comprise the steps:
(1) inputting the basic data of electrical network, before generating electric network element fault according to electrical network basic data, ground state electrical network leads
Receive matrix Y, store by sparse format；
(2) according to the sparsity structure of the admittance matrix Y of ground state electrical network, generate the Jacobian matrix J of sparse format, generating
During, for forming unified data form, and consider the node type conversion requirements in Load flow calculation, by joints all in electrical network
Pressing PQ node processing, redundancy opens up Jacobian matrix J, reserves the storage of non-zero entry for PV node redundancy in Jacobian matrix J
Deposit space and voltage magnitude increment Delta V2Storage area, i.e. form redundancy reactive power imbalance equation corresponding to PV node；
In the present embodiment, redundancy reactive power imbalance equation by all PV nodes is placed on the lower end of update equation group；According to admittance
All grid nodes are pressed PQ node processing by matrix, and the storage reserving non-zero entry for PV node redundancy in Jacobian matrix is empty
Between, in figure, the position of two null matrix 0 is i.e. the storage area of the non-zero entry reserved for PV node redundancy；
(3) when generating the update equation group of any one concrete Fault load flow of ground state electrical network, Jacobian matrix J according to
The actual node type of faulty grids carries out decoupling process；
The concrete step that processes is: as it is shown in figure 1, the PQ node of numbered i, arrange with merit imbalance Δ P (i) and idle not
Balance Δ Q (i) two equations, add voltage magnitude increment Delta V1(i) and two variablees of phase angle increment Delta θ (i), Δ V1I () is electricity
Pressure amplitude value increment Delta V1Middle i-th element, Δ θ (i) is i-th element in phase angle incremental vector Δ θ；The PV node of numbered j
Row write Δ P (j) equation, add Δ θ (j) variable, in the equation of redundancy Δ Q (j) of PV node, and unit in the non-zero entry that redundancy is reserved
Element is set to zero, i.e. the element that the respective redundant of two null matrix 0 positions in Fig. 1 is reserved is set to 0, with voltage magnitude Delta
ΔV2Element in j Jacobian matrix that () is corresponding is set to 1, i.e. forms unit diagonal matrix I as shown in Figure 1, corresponding nothing
Merit unbalanced power amount Δ Q2It is predisposed to 0.In Fig. 1, what H, N, J, L stored is the real effective information of Jacobian matrix, passes through
Above-mentioned process, the redundancy voltage magnitude Δ V obtained2=0, meet the requirement that PV node voltage amplitude is constant, it is ensured that PV node is superfluous
The remaining idle amount of unbalance Δ Q supplemented2Voltage magnitude increment Delta V that corresponding equation group and redundancy are supplemented2Do not affect correction side
Correctly solving of journey group.
Claims (2)
1. the Jacobian matrix redundancy storage method for static security analysis, it is characterised in that: described method include as
Lower step:
(1) according to the admittance matrix Y of ground state electrical network before electrical network basic data generation electric network element fault, store by sparse format；
(2) according to the sparsity structure of the admittance matrix Y of ground state electrical network, the Jacobian matrix J of sparse format is generated, in the process of generation
In, nodes all in electrical network are pressed PQ node processing, redundancy opens up Jacobian matrix J, for PV node in Jacobian matrix J superfluous
The storage area of remaining reserved non-zero entry and voltage magnitude increment Delta V2Storage area, i.e. form redundancy corresponding to PV node idle
Unbalanced power equation；
(3) when generating the update equation group of any one concrete Fault load flow of ground state electrical network, Jacobian matrix J is according to fault
The actual node type of electrical network carries out decoupling process: the idle amount of unbalance Δ Q supplemented by PV node redundancy2It is set to 0, will be all
In the non-zero entry of the Jacobian matrix that PV node redundancy supplements, element is set to 0, by the idle merit of redundancy corresponding for each PV node
With voltage magnitude Delta Δ V in rate imbalance equation2The corresponding element in Jacobian matrix is set to 1.
Jacobian matrix redundancy storage method for static security analysis the most according to claim 1, it is characterised in that:
Electrical network master data described in described step (1) includes node serial number, number of nodes, branch number, branch road quantity, branch road electricity
Anti-, the connected node serial number of branch road and parallel branch reactance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610662283.5A CN106294022B (en) | 2016-08-12 | 2016-08-12 | A kind of Jacobian matrix redundancy storage method for static security analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610662283.5A CN106294022B (en) | 2016-08-12 | 2016-08-12 | A kind of Jacobian matrix redundancy storage method for static security analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106294022A true CN106294022A (en) | 2017-01-04 |
| CN106294022B CN106294022B (en) | 2019-03-29 |
Family
ID=57669514
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610662283.5A Active CN106294022B (en) | 2016-08-12 | 2016-08-12 | A kind of Jacobian matrix redundancy storage method for static security analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106294022B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109344361A (en) * | 2018-08-27 | 2019-02-15 | 南昌大学 | A kind of electric power system tide calculate in the method that is quickly formed of Jacobian matrix |
| CN109447839A (en) * | 2018-09-28 | 2019-03-08 | 南昌大学 | The random method generated and form electric system Jacobian matrix of nonzero element |
| CN109599881A (en) * | 2018-12-19 | 2019-04-09 | 杭州电子科技大学 | A kind of power grid frequency modulation pressure regulation method based on lithium manganate battery energy-storage system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070027642A1 (en) * | 2005-07-15 | 2007-02-01 | Chang Gung University | Method for Calculating Power Flow Solution of a Power Transmission Network that Includes Interline Power Flow Controller (IPFC) |
| US20090182518A1 (en) * | 2008-01-14 | 2009-07-16 | Chia-Chi Chu | Method of Calculating Power Flow Solution of a Power Grid that Includes Generalized Power Flow Controllers |
| CN101621200A (en) * | 2009-08-04 | 2010-01-06 | 大连海事大学 | Rectangular coordinate Newton load flow calculation method |
| CN101957891A (en) * | 2010-09-14 | 2011-01-26 | 天津理工大学 | Static safety analysis system of modularization power grid on basis of communication system |
| CN103617150A (en) * | 2013-11-19 | 2014-03-05 | 国家电网公司 | GPU (graphic processing unit) based parallel power flow calculation system and method for large-scale power system |
-
2016
- 2016-08-12 CN CN201610662283.5A patent/CN106294022B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070027642A1 (en) * | 2005-07-15 | 2007-02-01 | Chang Gung University | Method for Calculating Power Flow Solution of a Power Transmission Network that Includes Interline Power Flow Controller (IPFC) |
| US20090182518A1 (en) * | 2008-01-14 | 2009-07-16 | Chia-Chi Chu | Method of Calculating Power Flow Solution of a Power Grid that Includes Generalized Power Flow Controllers |
| CN101621200A (en) * | 2009-08-04 | 2010-01-06 | 大连海事大学 | Rectangular coordinate Newton load flow calculation method |
| CN101957891A (en) * | 2010-09-14 | 2011-01-26 | 天津理工大学 | Static safety analysis system of modularization power grid on basis of communication system |
| CN103617150A (en) * | 2013-11-19 | 2014-03-05 | 国家电网公司 | GPU (graphic processing unit) based parallel power flow calculation system and method for large-scale power system |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109344361A (en) * | 2018-08-27 | 2019-02-15 | 南昌大学 | A kind of electric power system tide calculate in the method that is quickly formed of Jacobian matrix |
| CN109344361B (en) * | 2018-08-27 | 2022-05-20 | 南昌大学 | Method for quickly forming Jacobian matrix in power system load flow calculation |
| CN109447839A (en) * | 2018-09-28 | 2019-03-08 | 南昌大学 | The random method generated and form electric system Jacobian matrix of nonzero element |
| CN109447839B (en) * | 2018-09-28 | 2021-08-06 | 南昌大学 | Method for randomly generating non-zero elements and forming Jacobian matrix of power system |
| CN109599881A (en) * | 2018-12-19 | 2019-04-09 | 杭州电子科技大学 | A kind of power grid frequency modulation pressure regulation method based on lithium manganate battery energy-storage system |
| CN109599881B (en) * | 2018-12-19 | 2022-01-11 | 杭州电子科技大学 | Power grid frequency and voltage modulation method based on lithium manganate battery energy storage system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106294022B (en) | 2019-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106532711B (en) | Change the Newton load flow calculation method of Jacobian matrix with iteration and node type | |
| CN100466412C (en) | United accident-prevention rehearse method for provincial and city power system | |
| CN101572409B (en) | Self-adaptive device and method for analyzing global power flow of generation, transmission and distribution | |
| Eltamaly et al. | Load flow analysis by gauss-seidel method; a survey | |
| CN101895116B (en) | Online available transmission capacity computing method based on distribution factor | |
| CN111626633B (en) | Self-energy-storage flexible interconnected power distribution network expansion planning method based on comprehensive evaluation index | |
| CN112910012B (en) | Power distribution system elasticity improvement evaluation method, storage medium and computing equipment | |
| Biserica et al. | Network partition for coordinated control in active distribution networks | |
| CN106294022B (en) | A kind of Jacobian matrix redundancy storage method for static security analysis | |
| Wang et al. | ADMM-based distributed active and reactive power control for regional AC power grid with wind farms | |
| CN107808200B (en) | Short-time maintenance plan optimization method and device for power distribution network | |
| CN106410811B (en) | Iteration small impedance branches endpoint changes the tidal current computing method of Jacobian matrix for the first time | |
| CN105305429A (en) | Three-phase reclosing sequence resetting method for improving transient frequency stability of two-unit electric power system | |
| CN107230982A (en) | A kind of micro-capacitance sensor linearizes tidal current computing method | |
| CN106712029B (en) | The Newton load flow calculation method of small impedance branches PQ endpoint change Jacobian matrix | |
| Dezaki et al. | A new method based on sensitivity analysis to optimize the placement of SSSCs | |
| CN114665481A (en) | Urban power grid partitioning method and device based on power grid node vectorization | |
| CN113629787A (en) | Generator reactive power coordination mutual-aid compensation method considering sensitivity of center point | |
| CN108808683A (en) | Alternating current-direct current power grid load flow calculation method containing multi-balancing machine and multimachine joint pressure regulation | |
| Alkano et al. | Distributed MPC for Power-to-Gas facilities embedded in the energy grids | |
| CN105005831B (en) | A kind of calculation method of the quasi- direct current Dynamic Optimal Power Flow Problem based on electric power system tide coupled relation | |
| CN110518603A (en) | The complete distribution power system load flow calculation method of the mesolow of three-phase imbalance containing looped network | |
| Yang et al. | Distribution feeder parameter estimation without synchronized phasor measurement by using radial basis function neural networks and multi-run optimization method | |
| CN103904649A (en) | Generator polymerization method suitable for models with different orders | |
| Liu et al. | Distributed Optimal Operation of AC/DC Hybrid Distribution Network Considering Wind Power Uncertainty |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |