CN106294022A - A kind of Jacobian matrix redundancy storage method for static security analysis - Google Patents

Abstract

The invention discloses a Jacobian matrix redundant storage method for static security analysis, which comprises the following steps: (1) generating the admittance matrix Y of the ground state power grid before the fault of the power grid element according to the basic data of the power grid, and storing it in a sparse format ; (2) According to the sparse structure of the admittance matrix Y of the ground state power grid, generate a sparse format Jacobian matrix J, treat all nodes in the power grid as PQ nodes, and reserve non-zero redundancy for PV nodes in the Jacobian matrix J The storage space of the element and the storage space of the voltage amplitude increment ΔV ₂ form the redundant reactive power unbalance equation corresponding to the PV node; (3) when generating the correction equations for any specific fault power flow of the base state power grid, The Jacobian matrix J is decoupled according to the actual node type of the fault grid. The invention solves the problem of inconsistency of the structure of the Jacobian matrix of the batch processing power flow, and provides a data structure basis for the parallel calculation of the batch processing power flow.

Description

A Jacobian Matrix Redundancy Storage Method for Static Security Analysis

technical field

The invention belongs to the application field of high-performance computing in power systems, and in particular relates to a Jacobian matrix redundant storage method for static security analysis.

Background technique

Power flow calculation is the most widely used, basic and important electrical calculation in power system. In the study of power system operation mode and planning scheme, power flow calculation is required to compare the feasibility, reliability and economy of the operation mode or planning power supply scheme. In the actual production process, both offline power flow and online power flow calculations have relatively high requirements on the calculation speed of power flow.

The sparse format of the Jacobian matrix is related to the admittance matrix and the node type. The node type does not change, the Jacobian matrix dimension remains unchanged, and the sparse format remains unchanged; if the node type definitions are different, even the Jacobian generated by the same grid The order and sparse format of the ratio matrix are also different. The fault power flow in the static safety analysis is generated by several disconnections on the ground state power grid, and any N-1 disconnection fault affects at most 4 elements of the node admittance matrix of the base state power grid, among which the off-diagonal elements (mutual admittance) It may change from non-zero elements to zero elements, that is to say, the sparse format of any fault power flow is a subset of the sparse format of the ground state power flow (specifically refers to the number and format of non-zero elements), N-2 and other more complex static Security analysis is similar. This ensures the consistency of the admittance matrix in the power flow calculation and the admittance matrix data referenced when using the batch processing power flow method.

The batch power flow technology has strict requirements on the data format, and the data in each power flow calculation process needs to maintain a consistent structure. The Newton-Raphson method needs to calculate the Jacobian matrix, which is also the coefficient matrix of the sparse linear equation system, and it is also necessary to ensure the sparse format consistency of the Jacobian matrix.

Contents of the invention

Purpose of the invention: Aiming at the deficiencies of the prior art, the present invention provides a Jacobian matrix redundant storage method for static security analysis, which solves the problem of structural inconsistency in the fault flow Jacobian matrix for static security analysis, and can realize GPU Combined memory access greatly improves memory access efficiency and speeds up the calculation speed of GPU-accelerated static security analysis.

Power flow calculation: a term in electricity science, which refers to the calculation of the distribution of active power, reactive power and voltage in the power grid under the given power system network topology, component parameters, and power generation and load parameters.

Admittance matrix: A matrix established on the basis of the equivalent admittance of system components to describe the relationship between the voltage of each node of the power network and the injected current.

The present invention is a Jacobian matrix redundant storage method for static security analysis, said method comprising the following steps:

(1) Generate the admittance matrix Y of the ground state power grid before the power grid component failure according to the basic data of the power grid, and store it in a sparse format;

(2) According to the sparse structure of the admittance matrix Y of the ground-state power grid, generate the Jacobian matrix J in a sparse format. During the generation process, treat all nodes in the power grid as PQ nodes, and open up the Jacobian matrix J redundantly, which is a PV node In the Jacobian matrix J, the storage space of non-zero elements and the storage space of the voltage amplitude increment ΔV ₂ are redundantly reserved, that is, the redundant reactive power unbalance equation corresponding to the PV node is formed;

(3) When generating the correction equations for any specific fault power flow of the base state grid, the Jacobian matrix J is decoupled according to the actual node type of the fault grid: the reactive power unbalance ΔQ ₂ supplemented by PV node redundancy is set as 0, set the elements in the non-zero elements of the Jacobian matrix of all PV node redundancy supplements to 0, and set the Jacobian corresponding to the voltage amplitude variable increment in the redundant reactive power imbalance equation corresponding to each PV node The elements in the matrix are set to 1.

Wherein, the basic grid data in the step (1) includes node numbers, node numbers, branch numbers, branch numbers, branch reactances, branch connected node numbers and parallel branch reactances.

Beneficial effects: Compared with the prior art, the present invention has the following advantages: First, the storage space in the present invention is reserved for PQ node redundancy according to all nodes, and has a unified sparse format, so that the data in each power flow calculation process The structure is consistent, which is suitable for the data merge access requirements of GPU and other SIMT architecture hardware. On this basis, the merge access of GPU memory access is realized, which greatly improves the memory access efficiency and can speed up the calculation speed of GPU-accelerated static security analysis; When generating the Jacobian matrix of the specific fault power flow, it is processed according to the actual node type of the power grid, so as to ensure that the added redundant data will not affect the solution of the original correction equations; finally, the present invention solves the problem of inconsistent structure of the Jacobian matrix of the batch processing power flow, and provides Parallel calculation of batch fault power flow in static security analysis provides the basis of data structure.

Description of drawings

FIG. 1 is a schematic diagram of the redundant storage method of the Jacobian matrix of the present invention.

detailed description

The technical solution of the present invention will be further described below in conjunction with the accompanying drawings.

As shown in Figure 1, a kind of Jacobian matrix redundant storage method that the present invention is used for static safety analysis, described method comprises the following steps:

(1) Input the basic data of the power grid, generate the admittance matrix Y of the ground state power grid before the failure of the power grid components according to the basic data of the power grid, and store it in a sparse format;

(2) According to the sparse structure of the admittance matrix Y of the ground state power grid, the Jacobian matrix J in a sparse format is generated. During the generation process, in order to form a unified data format and consider the node type conversion requirements in the power flow calculation, All nodes are treated as PQ nodes, and the Jacobian matrix J is opened for redundancy, and storage space for non-zero elements and voltage amplitude increment ΔV ₂ is reserved for PV nodes in the Jacobian matrix J redundantly, that is, PV nodes are formed Corresponding redundant reactive power imbalance equation; in this embodiment, the redundant reactive power imbalance equations of all PV nodes are placed at the lower end of the correction equation group; according to the admittance matrix, all grid nodes are treated as PQ nodes, Reserve non-zero element storage space for PV node redundancy in the Jacobian matrix. The position of the two zero matrix 0 in the figure is the non-zero element storage space reserved for PV node redundancy;

(3) When generating the correction equations for any specific fault power flow of the base state grid, the Jacobian matrix J is decoupled according to the actual node type of the fault grid;

The specific processing steps are: as shown in Figure 1, for the PQ node numbered i, write two equations of active power unbalance ΔP(i) and reactive power unbalance ΔQ(i), and add voltage amplitude increment ΔV ₁ (i ) and the phase angle increment Δθ(i) two variables, ΔV ₁ (i) is the i-th element in the voltage amplitude increment ΔV ₁ , and Δθ(i) is the i-th element in the phase angle increment vector Δθ; The PV node numbered j writes the ΔP(j) equation and adds the Δθ(j) variable. In the redundant ΔQ(j) equation of the PV node, the elements in the redundant reserved non-zero elements are set to zero, that is, The corresponding redundant reserved elements in the two zero matrix 0 positions in 1 are set to 0, and the elements in the Jacobian matrix corresponding to the voltage amplitude variable increment ΔV ₂ (j) are set to 1, which is formed as shown in Figure 1 For the unit diagonal matrix I shown in , the corresponding reactive power unbalance ΔQ ₂ is preset to 0. In Figure 1, H, N, J, L store the real and effective information of the Jacobian matrix. Through the above-mentioned processing, the redundant voltage amplitude ΔV ₂ = 0 is obtained, which meets the requirement that the voltage amplitude of the PV node remains unchanged. It is ensured that the equations corresponding to the reactive power unbalance ΔQ ₂ supplemented by PV node redundancy and the voltage amplitude increment ΔV ₂ supplemented by redundancy will not affect the correct solution of the correction equations.

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 V₂Storage 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 redundancy₂It 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 equation₂The 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.