CN106294022B - 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 PDF

Info

Publication number
CN106294022B
CN106294022B CN201610662283.5A CN201610662283A CN106294022B CN 106294022 B CN106294022 B CN 106294022B CN 201610662283 A CN201610662283 A CN 201610662283A CN 106294022 B CN106294022 B CN 106294022B
Authority
CN
China
Prior art keywords
node
jacobian matrix
redundancy
power grid
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.)
Active
Application number
CN201610662283.5A
Other languages
Chinese (zh)
Other versions
CN106294022A (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201610662283.5A priority Critical patent/CN106294022B/en
Publication of CN106294022A publication Critical patent/CN106294022A/en
Application granted granted Critical
Publication of CN106294022B publication Critical patent/CN106294022B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply

Abstract

The invention discloses a kind of Jacobian matrix redundancy storage methods for static security analysis, and this method comprises the following steps: (1) generating the admittance matrix Y of ground state power grid before electric network element failure according to electrical network basic data, store by sparse format；(2) according to the sparsity structure of the admittance matrix Y of ground state power grid, generate the Jacobian matrix J of sparse format, nodes all in power grid are pressed into PQ node processing, are storage space and voltage magnitude increment Delta V that PV node redundancy in Jacobian matrix J reserves non-zero entry2Storage space, that is, form the corresponding redundancy reactive power imbalance equation of PV node；(3) when generating the update equation group of any one specific Fault load flow of ground state power grid, Jacobian matrix J carries out decoupling processing according to the actual node type of faulty grids.The present invention solves the problems, such as that batch processing Load Flow Jacobian Matrix structure is inconsistent, provides data structure basis for batch processing power flow parallel computing.

Description

A kind of Jacobian matrix redundancy storage method for static security analysis
Technical field
The invention belongs to High performance computing in power system application fields more particularly to a kind of for the refined of static security analysis Than matrix redundancy storage method.
Background technique
Load flow calculation is most widely used, most basic and most important a kind of electrical operation in electric system.In power train In the research of the method for operation of uniting and programme, require to carry out Load flow calculation to compare the method for operation or plan power supply plan Feasibility, reliability and economy.And in actual production process, no matter offline trend and online power flow calculate all to the meter of trend Calculate the requirement that speed has this relatively high.
The sparse format of Jacobian matrix is related to admittance battle array, while related to node type, and node type does not change Become, Jacobian matrix dimension is constant, and sparse format is constant；If site type definition is different, even if being generated by the same power grid Jacobian matrix order and sparse format it is also not identical.Fault load flow in static security analysis is sent out on ground state power grid What raw several broken strings generated, any one N-1 disconnection fault at most influences 4 elements of ground state grid nodes admittance battle array, wherein Nondiagonal element (transadmittance) may become null element by non-zero entry, that is to say, that the sparse format of any one Fault load flow is all base The a subset (refering in particular to non-zero entry quantity and format) of state trend sparse format, N-2 and other more complicated static security analysis Also similar.This guarantees admittance battle arrays when Load flow calculation, use one of reference admittance battle array data when batch processing trend method Cause property.
Batch processing band wagon technique has strict requirements to data format, 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 simultaneously sparse vectors Coefficient matrix also needs the sparse format consistency for guaranteeing Jacobian matrix.
Summary of the invention
Goal of the invention: in view of the deficiencies of the prior art, the present invention provides a kind of Jacobi for static security analysis Matrix redundancy storage method, this method solve structure inconsistence problems in static security analysis Fault load flow Jacobian matrix, The merging access that GPU memory access can be achieved, greatly improves memory access efficiency, accelerates the calculating speed for the static security analysis that GPU accelerates.
Load flow calculation: electrodynamic noun refers in given power system network topology, component parameters and power generation, load parameter Under the conditions of, calculate the distribution of active power, reactive power and voltage in power network.
Admittance matrix: established based on the Equivalent admittance of system element, each node voltage of description electric power networks and The matrix of relationship between Injection Current.
A kind of Jacobian matrix redundancy storage method for static security analysis of the present invention, the method includes walking as follows It is rapid:
(1) the admittance matrix Y that ground state power grid before electric network element failure is generated according to electrical network basic data, is deposited by sparse format Storage；
(2) according to the sparsity structure of the admittance matrix Y of ground state power grid, the Jacobian matrix J of sparse format is generated, is being generated In the process, nodes all in power grid are pressed into PQ node processing, redundancy opens up Jacobian matrix J, is PV node in Jacobian matrix J Middle redundancy reserves the storage space and voltage magnitude increment Delta V of non-zero entry2Storage space, that is, form the corresponding redundancy of PV node Reactive power imbalance equation；
(3) generate ground state power grid any one specific Fault load flow update equation group when, Jacobian matrix J according to The actual node type of faulty grids carries out decoupling processing: the idle amount of unbalance Δ Q that PV node redundancy is supplemented2It is set to 0, it will Element is set to 0 in the non-zero entry of the Jacobian matrix of all PV node redundancies supplement, by the corresponding redundancy of each PV node without Element in function unbalanced power equation in Jacobian matrix corresponding with voltage magnitude Delta is set to 1.
Wherein, power grid master data described in the step (1) includes node serial number, number of nodes, branch number, branch Quantity, branch reactance, the connected node serial number of branch and parallel branch reactance.
The utility model has the advantages that compared with prior art, the present invention has the following advantages: firstly, memory space is by all in the present invention Node is all that PQ node redundancy is reserved, has unified sparse format, thus the data structure during each Load flow calculation It is consistent, is suitble to the SIMT framework hardware such as GPU that data are merged with the demand visited, realize the merging of GPU memory access on this basis Access, greatly improves memory access efficiency, can accelerate the calculating speed of the static security analysis of GPU acceleration；Secondly specific event is being generated It is handled when hindering the Jacobian matrix of trend according to power grid actual node type, guarantees that the redundant data of addition will not influence original The solution of update equation group；The last present invention solves the problems, such as that batch processing Load Flow Jacobian Matrix structure is inconsistent, for static state The parallel computation of batch Fault load flow provides data structure basis in safety analysis.
Detailed description of the invention
Fig. 1 is the schematic diagram of the redundancy storage method of Jacobian matrix of the present invention.
Specific embodiment
Technical solution of the present invention is described further with reference to the accompanying drawing.
As shown in Figure 1, a kind of Jacobian matrix redundancy storage method for static security analysis of the present invention, the method Include the following steps:
(1) basic data of power grid is inputted, ground state power grid is led before generating electric network element failure according to electrical network basic data It receives matrix Y, is stored by sparse format；
(2) according to the sparsity structure of the admittance matrix Y of ground state power grid, the Jacobian matrix J of sparse format is generated, is being generated In the process, to form unified data format, and the node type conversion requirements in Load flow calculation are considered, by sections all in power grid PQ node processing is pressed, redundancy opens up Jacobian matrix J, is the storage that PV node redundancy in Jacobian matrix J reserves non-zero entry Deposit space and voltage magnitude increment Delta V2Storage space, that is, form the corresponding redundancy reactive power imbalance equation of PV node； The redundancy reactive power imbalance equation of all PV nodes is placed on to the lower end of update equation group in the present embodiment；According to admittance All grid nodes are pressed PQ node processing by matrix, are the storage sky that PV node redundancy in Jacobian matrix reserves non-zero entry Between, the position of two null matrix 0 is the storage space for the non-zero entry reserved for PV node redundancy in figure；
(3) generate ground state power grid any one specific Fault load flow update equation group when, Jacobian matrix J according to The actual node type of faulty grids carries out decoupling processing；
Specific processing step are as follows: as shown in Figure 1, number be i PQ node, column with function imbalance Δ P (i) and it is idle not Δ Q (i) two equations are balanced, voltage magnitude increment Delta V is added1(i) and two variables of phase angle increment Delta θ (i), Δ V1It (i) is electricity Pressure amplitude value increment Delta V1In i-th of element, Δ θ (i) is i-th of element in phase angle incremental vector Δ θ；The PV node that number is j Column write Δ P (j) equation, add Δ θ (j) variable, in the equation of the redundancy Δ Q (j) of PV node, member in the reserved non-zero entry of redundancy Element is set to zero, i.e. the element that the respective redundant of two 0 positions of null matrix in Fig. 1 is reserved is set to 0, with voltage magnitude Delta ΔV2(j) element in corresponding Jacobian matrix is set to 1, that is, forms unit diagonal matrix I as shown in Figure 1, corresponding nothing Function unbalanced power amount Δ Q2It is predisposed to 0.In Fig. 1, what H, N, J, L were stored is the real effective information of Jacobian matrix, is passed through Above-mentioned processing, the redundancy voltage magnitude Δ V found out2=0, meet the constant requirement of PV node voltage magnitude, guarantees that PV node is superfluous The idle amount of unbalance Δ Q of remaining supplement2The voltage magnitude increment Delta V of corresponding equation group and redundancy supplement2Do not influence amendment side The correct solution of journey group.

Claims (2)

1. a kind of Jacobian matrix redundancy storage method for static security analysis, it is characterised in that: the method includes such as Lower step:
(1) the admittance matrix Y that ground state power grid before electric network element failure is generated according to electrical network basic data, is stored by sparse format；
(2) according to the sparsity structure of the admittance matrix Y of ground state power grid, the Jacobian matrix J of sparse format is generated, in generating process In, nodes all in power grid are pressed into PQ node processing, wherein the step of PQ node processing is that all PQ nodes are arranged with function injustice Weigh Δ P and idle uneven Δ two equations of Q, adds voltage magnitude increment Delta V1With two variables of phase angle increment Delta θ；Redundancy is opened Jacobian matrix J is warded off, is storage space and voltage magnitude increment Delta that PV node redundancy in Jacobian matrix J reserves non-zero entry V2Storage space, that is, form the corresponding redundancy reactive power imbalance equation of PV node；
(3) when generating the update equation group of any one specific Fault load flow of ground state power grid, Jacobian matrix J is according to failure The actual node type of power grid carries out decoupling processing: the idle amount of unbalance Δ Q that PV node redundancy is supplemented2It is set to 0, will be owned Element is set to 0 in the non-zero entry of the Jacobian matrix of PV node redundancy supplement, by the idle function of the corresponding redundancy of each PV node In rate imbalance equation with voltage magnitude Delta Δ V2Element in corresponding Jacobian matrix is set to 1.
2. the Jacobian matrix redundancy storage method according to claim 1 for static security analysis, it is characterised in that: Electrical network basic data described in the step (1) includes node serial number, number of nodes, branch number, number of branches, branch electricity Anti-, the connected node serial number of branch and parallel branch reactance.
CN201610662283.5A 2016-08-12 2016-08-12 A kind of Jacobian matrix redundancy storage method for static security analysis Active CN106294022B (en)

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 CN106294022A (en) 2017-01-04
CN106294022B true 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)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109344361B (en) * 2018-08-27 2022-05-20 南昌大学 Method for quickly forming Jacobian matrix in power system load flow calculation
CN109447839B (en) * 2018-09-28 2021-08-06 南昌大学 Method for randomly generating non-zero elements and forming Jacobian matrix of power system
CN109599881B (en) * 2018-12-19 2022-01-11 杭州电子科技大学 Power grid frequency and voltage modulation method based on lithium manganate battery energy storage system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7321834B2 (en) * 2005-07-15 2008-01-22 Chang Gung University Method for calculating power flow solution of a power transmission network that includes interline power flow controller (IPFC)
US7813884B2 (en) * 2008-01-14 2010-10-12 Chang Gung University Method of calculating power flow solution of a power grid that includes generalized power flow controllers

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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

Also Published As

Publication number Publication date
CN106294022A (en) 2017-01-04

Similar Documents

Publication Publication Date Title
Liu et al. A transactive energy framework for coordinated energy management of networked microgrids with distributionally robust optimization
CN100466412C (en) United accident-prevention rehearse method for provincial and city power system
CN106532711B (en) Change the Newton load flow calculation method of Jacobian matrix with iteration and node type
Eltamaly et al. Load flow analysis by gauss-seidel method; a survey
CN101527455B (en) Interconnected electric network distributed current calculating method on the basis of alternation and iteration of current module
CN106407158B (en) A kind of batch processing isomorphism sparse matrix that GPU accelerates multiplies the processing method of full vector
CN103077480B (en) Safety checking method for power system
CN106294022B (en) A kind of Jacobian matrix redundancy storage method for static security analysis
CN105514971A (en) Flow calculation method suitable for microgrids in various operation modes
CN106684857A (en) Linear optimal power flow model containing unified power flow controller
CN103023013A (en) Method for calculating three-phase load flow of active distribution network
CN107204617A (en) The Interval Power Flow computational methods of Cartesian form based on linear programming
CN103956735A (en) Harmonic power flow analysis method of distributed power generation system
CN110120673A (en) Distributed transmission & distribution collaboration idle work optimization method and system based on Thevenin's equivalence parameter identification
CN106253288A (en) A kind of optimal load flow algorithm containing THE UPFC based on automatic differential
CN103515964A (en) Reactive compensation control method and reactive compensation control device
CN106410811B (en) Iteration small impedance branches endpoint changes the tidal current computing method of Jacobian matrix for the first time
CN112531756A (en) Distributed control method, system and equipment for electric quantity balance of energy storage system
CN113097994A (en) Power grid operation mode adjusting method and device based on multiple reinforcement learning agents
CN105305429A (en) Three-phase reclosing sequence resetting method for improving transient frequency stability of two-unit electric power system
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
CN110649624B (en) Power system load flow parallel computing method
CN105470946A (en) Dynamo-electric transient simulation algorithm for power distribution network containing distributed power supply
CN108808683A (en) Alternating current-direct current power grid load flow calculation method containing multi-balancing machine and multimachine joint pressure regulation

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