CN107579516B - Method for improving state estimation calculation speed of power system - Google Patents
Method for improving state estimation calculation speed of power system Download PDFInfo
- Publication number
- CN107579516B CN107579516B CN201710823171.8A CN201710823171A CN107579516B CN 107579516 B CN107579516 B CN 107579516B CN 201710823171 A CN201710823171 A CN 201710823171A CN 107579516 B CN107579516 B CN 107579516B
- Authority
- CN
- China
- Prior art keywords
- matrix
- column
- zero
- reactive
- state estimation
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000011159 matrix material Substances 0.000 claims abstract description 119
- 238000004458 analytical method Methods 0.000 claims abstract description 42
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 31
- 238000005516 engineering process Methods 0.000 claims abstract description 18
- 238000005259 measurement Methods 0.000 claims description 34
- 230000011664 signaling Effects 0.000 claims description 11
- 238000012986 modification Methods 0.000 claims description 8
- 230000004048 modification Effects 0.000 claims description 8
- 238000010845 search algorithm Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims 1
- 239000002994 raw material Substances 0.000 claims 1
- 238000007781 pre-processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a method for improving the state estimation calculation speed of a power system, which realizes the rapid calculation of an active and reactive Jacobian matrix and an information matrix by using a multithread parallel algorithm based on an OpenMP shared memory programming mode, reduces the calculated amount by using a sparse technology through zero-removing operation in the matrix multiplication process, and accelerates the decomposition speed of the information matrix based on non-zero-element symbol analysis and a numerical decomposition method in the factor decomposition process, thereby improving the overall calculation speed of large-scale power grid state estimation.
Description
Technical Field
The invention relates to a method for improving the state estimation calculation speed of a power system, and belongs to the technical field of power automation.
Background
With the construction of dispatching systems and model data centers of all levels of national provinces of intelligent dispatching technical support systems (D5000) and the development of power grid scale, the state estimation calculation scale of a power system is increased day by day, the state estimation of actual production operation in the current intelligent dispatching technical support system, particularly the state estimation of the operation of the model data center, the modeling range of the state estimation calculation scale is from 1000kV to 10kV feeder outlet lines of an extra-high voltage feeder, and the state estimation calculation scale is very large, so that the calculation efficiency of the state estimation calculation scale cannot meet the requirements of real-time calculation analysis. Therefore, improving the calculation efficiency of large-scale grid state estimation becomes a focus of attention of many researchers.
At present, in the state estimation calculation of the power system, the calculation of an information matrix and the LU factorization of a coefficient matrix of a linear equation set occupy most of the state estimation time, wherein for a large-scale sparse linear system, the proportion of the factorization calculation time to the total calculation time of solving the equation set is large, and the speed of improving the information matrix and the LU factorization plays a significant role in accelerating the state of the whole power system.
Disclosure of Invention
In order to solve the technical problem, the invention provides a method for improving the calculation speed of the state estimation of the power system.
In order to achieve the purpose, the invention adopts the technical scheme that:
a method for improving the calculation speed of power system state estimation comprises the following steps,
converting a relational power grid model in a dispatching control system into a hierarchical power grid model;
reading the measurement of the SCADA in the dispatching system, and associating the measurement with an element in the hierarchical power grid model;
determining a network topology, namely a node-branch model, according to the remote signaling and the connection relation of each element;
node sorting is carried out based on the measurement number associated with the node;
forming a functional Jacobian matrix H based on OpenMP technology according to a node-branch modelpAnd a reactive Jacobian matrix Hq;
Active information matrix G calculated based on sparse matrix technology and OpenMP technologypAnd a reactive information matrix Gq,GpIs a (n-1) × (n-1) -dimensional matrix, GqThe method is characterized in that the method is an n × n-dimensional matrix, wherein n is the number of nodes calculated by a power grid;
to active information matrix GpThe analysis of the row number sets of the non-zero elements of each column of the lower triangular matrix L and the upper triangular matrix U is completed in sequence from the 1 st column to the n-1 st column according to a non-zero element symbol analysis method;
for reactive information matrix GqThe analysis of the row number sets of the non-zero elements of each column of the lower triangular matrix L and the upper triangular matrix U is completed in sequence from the 1 st column to the nth column according to a non-zero element symbol analysis method;
to active information matrix GpAccording to the Left Looking LU numerical decomposition calculation analysis method, combining the obtained non-zero row number set to sequentially perform Left Looking LU numerical decomposition from the 1 st row to the n-1 st row;
for reactive information matrix GqAccording to the Left Looking LU numerical decomposition calculation analysis method, combining the obtained non-zero row number set to sequentially carry out Left Looking LU numerical decomposition from the 1 st row to the n th row;
iteratively solving the system of equations G Δ x repeatedly during the state estimation computation processk=HTR-1·[z-h]And correcting the state quantity until DeltaxkSatisfying a specified convergence criterion; wherein G is an active information matrix GpOr reactive information matrix Gq,xkFor the calculation of the kth iteration the node voltage amplitude or phase angle, i.e. the state quantity, xk+1=xk+ΔxkH is a nonlinear measurement vector function, z is a measurement vector, and when G is an active information matrix GpR is thenH is HpWhen G is a reactive information matrix GqR is thenH is Hq,Andrespectively an active measurement error variance matrix and a reactive measurement error variance matrix;
when Δ xkAnd when the iteration threshold is smaller than the iteration threshold, stopping iterative computation.
The converted hierarchical network model is stored in a hierarchical database, and measurements are also stored in the hierarchical database and associated with elements in the hierarchical network model.
Remote signaling in measurement comprises a breaker state and a disconnecting link state; when the state estimation is calculated based on a new power grid model, network topology analysis is carried out based on the whole network elements; when the state of the circuit breaker and the disconnecting link exceeds the threshold value compared with the last state change number, network topology analysis is carried out based on the whole network elements; and when the current state of the circuit breaker and the disconnecting link does not exceed the threshold value compared with the last state change number, only carrying out local network topology analysis on the plant station with the changed remote signaling state according to the voltage grade.
Decoupling an active modification equation and a reactive modification equation by adopting a rapid decomposition state estimation algorithm based on the least square principle, and calculating an active Jacobian matrix H in parallelpAnd a reactive Jacobian matrix Hq。
According to the formed active Jacobian matrix HpJacobian matrix H without powerqActive power measurement error variance matrixSum reactive power measurement error variance matrixActive information matrix calculation based on sparse matrix technology and OpenMP technologyAnd a reactive information matrix
Setting a certain column vector of the coefficient matrix A ═ LU as b, and obtaining a non-zero element structure of a column vector x by solving Lx ═ b;
let β be { i | b }i≠0},χ={j|xjNot equal to 0 represents a collection of nodes of non-zero elements in b and x, respectively, where biIs the ith row element, x, in column vector bjIs the jth row element in the column vector x, where i ∈ [1, n'],j∈[1,n′]When the coefficient matrix A in the state estimation calculation is an active information matrix GpN is n-1, and the coefficient matrix A is a reactive information matrix G in the state estimation calculationq,n′=n;
Suppose that the k' -1 column of L has been computed, and its corresponding directed graph is G (L)k′-1);
L, U column k' is a set of non-zero nodes derived by the following criteria,
wherein lijNot equal to 0 denotes G (L)k′-1) Where there is an edge, x, from node j to node iiIs the ith row element in the column vector x;
when solving for the non-zero structure of x, the zero values produced by the numerical cancellation are ignored, since lij*xjResult in no matter biWhether or not it is 0, xiIs non-zero;
a and L are stored sparsely according to columns, the row number of a certain column of non-zero elements of A is retrieved, the non-zero elements of the column are retrieved from the corresponding column of L through the row number, and a non-zero set of x is obtained through a depth search algorithm.
The invention achieves the following beneficial effects: 1. the fast calculation of an active jacobian matrix, a reactive jacobian matrix and an information matrix is realized by using a multithreading parallel algorithm based on an OpenMP shared memory programming mode, the calculated amount is reduced by using a sparse technology through zero elimination operation in a matrix multiplication process, and the decomposition speed of the information matrix is accelerated based on non-zero element symbol analysis and a numerical decomposition method in a factor decomposition process, so that the overall calculation speed of large-scale power grid state estimation is improved; 2. based on the local topology technology, the time consumed by topology analysis during remote signaling preprocessing in large-scale power grid state estimation can be greatly reduced, and the occurrence probability of time consumed by topology analysis based on a whole network model before and after remote signaling preprocessing is effectively avoided; 3. in the LU decomposition process, a method of symbol analysis and numerical decomposition separation is adopted, so that the calculation of zero elements in the decomposition process can be effectively avoided.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a flow of parallel computation of an information matrix;
FIG. 3 is a schematic diagram of non-zero structure analysis;
FIG. 4 is a schematic diagram of solving L, U column j.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
As shown in fig. 1, a method for increasing the calculation speed of power system state estimation includes the following steps:
Because the power grid structure of the power system has a hierarchical relationship, in order to facilitate analysis and calculation, an analysis and calculation method based on a hierarchical database is adopted in state estimation and calculation, firstly, a relational power grid model in a regulation and control system is converted into a hierarchical power grid model, the power grid model is stored in the hierarchical database, the hierarchical database stores the power grid model on a state estimation application server hard disk in a binary file mode in a disk file mapping mode, elements in the hierarchical database correspond to elements in the relational database in a scheduling and control system one by one, and calculation and analysis are carried out based on the hierarchical database during calculation and calculation.
And 2, reading the measurement of the SCADA in the dispatching system, and associating the measurement with an element in the hierarchical power grid model.
The state estimation obtains the measurement of the SCADA from the dispatching control system, wherein the remote signaling comprises a breaker state and a disconnecting link state, the remote measurement comprises line active, reactive and current, transformer active, reactive and current, gear, load active, reactive and current, generator active, reactive and current and bus voltage, and as elements of the hierarchical database and elements of the relational database in the dispatching control system are in one-to-one correspondence through keywords, the measurement acquired by the SCADA is easily related to corresponding elements in the hierarchical database.
And 3, determining a network topology, namely a node-branch model, according to the state of the circuit breaker, the state of the disconnecting link and the connection relation of each element.
When the state estimation is calculated based on a new power grid model, network topology analysis is carried out based on the whole network elements; and when the current state of the circuit breaker and the disconnecting link exceeds the threshold value compared with the last state change number, performing network topology analysis based on the whole network elements. In an actual power grid, the states of a plurality of circuit breakers and disconnecting links rarely change during normal operation, so that when the number of the state changes of the circuit breakers and the disconnecting links at this time is not more than a threshold value compared with the number of the state changes at the last time, only the plant stations with the changed remote signaling states are subjected to local network topology analysis according to the voltage levels. Based on the local topology technology, the time consumed by topology analysis during remote signaling preprocessing in large-scale power grid state estimation can be greatly reduced, and the occurrence probability of time consumed by topology analysis based on a whole network model before and after remote signaling preprocessing is effectively avoided.
Step 4, node sequencing is carried out based on the measurement number of node association, and a functional Jacobian matrix H is formed based on OpenMP (open multiprocessing) technology according to a node-branch modelpAnd a reactive Jacobian matrix Hq。
Decoupling an active modification equation and a reactive modification equation by adopting a rapid decomposition state estimation algorithm based on the least square principle, and calculating an active Jacobian matrix H in parallelpAnd a reactive Jacobian matrix Hq。
Step 5, according to the formed active Jacobian matrix HpJacobian matrix H without powerqTherein are provided withPower measurement error variance matrixSum reactive power measurement error variance matrixActive information matrix calculation based on sparse matrix technology and OpenMP technologyAnd a reactive information matrixGpIs a (n-1) × (n-1) -dimensional matrix, GqThe method is characterized in that the method is an n × n-dimensional matrix, wherein n is the number of nodes calculated by a power grid;
implementation of G using OpenMP techniquesp、GqBased on sparse techniques, computing Gp、GqThe specific flow of a certain column of values is shown in fig. 2.
Step 6, for the active information matrix GpThe analysis of the row number sets of the non-zero elements of each column of the lower triangular matrix L and the upper triangular matrix U is completed in sequence from the 1 st column to the n-1 st column according to a non-zero element symbol analysis method; for reactive information matrix GqAnd the analysis of the row number sets of the non-zero elements of the columns of the lower triangular matrix L and the upper triangular matrix U is completed in sequence from the 1 st column to the nth column according to a non-zero element symbol analysis method.
Setting a coefficient matrix A as LU column vector as b, and solving Lx as b to obtain a non-zero element structure of a column vector x;
let β be { i | b }i≠0},χ={j|xjNot equal to 0 represents a collection of nodes of non-zero elements in b and x, respectively, where biIs the ith row element, x, in column vector bjIs the jth row element in the column vector x, where i ∈ [1, n'],j∈[1,n′]When the coefficient matrix A in the state estimation calculation is an active information matrix GpN is n-1, and the coefficient matrix A is a reactive information matrix G in the state estimation calculationq,n′=n;
Suppose that the k' -1 column of L has been computed and the corresponding directed graph isG(Lk′-1);
L, U column k' is a set of non-zero nodes derived by the following criteria,
wherein lijNot equal to 0 denotes G (L)k′-1) Where there is an edge, x, from node j to node iiIs the ith row element in the column vector x;
when solving for the non-zero structure of x, the zero values produced by the numerical cancellation are ignored, since lij*xjResult in no matter biWhether or not it is 0, xiA non-zero, non-zero element structural analysis diagram is shown in fig. 3;
and B, because A and L are stored sparsely according to columns, the row number of a certain column of non-zero elements of A is retrieved, the non-zero elements of the column are retrieved from the column corresponding to L through the row number, and the non-zero set of x is obtained through a depth search algorithm.
Step 7, for the active information matrix GpAccording to a Left Looking LU numerical decomposition calculation analysis method, combining the obtained non-zero row number set to sequentially perform Left Looking LU numerical decomposition from the 1 st column to the n-1 st column; for reactive information matrix GqAnd according to a Left Looking LU numerical decomposition calculation analysis method, combining the obtained non-zero row number set to sequentially perform Left Looking LU numerical decomposition from the 1 st column to the n th column.
The Left Looking LU numerical decomposition is sequentially decomposed from the 1 st column to the last column, namely, one column vector is calculated each time, and numerical calculation is carried out according to the symbol analysis result of the column, so that the numerical values of all non-zero elements of the column are obtained.
Solving L, U column j is schematically shown in FIG. 4, where matrix block L in FIG. 4 isj、L′j、UjAnd vector aj、a′j、lj、uj、x′jIs defined as follows;
aj=(a1j,…,a(j-1)j)T
a′j=(ajj,…,a′nj)T
lj=(ljj,…,l′cj)T
uj=(u1j,…,u(j-1)j)T
x′j=(xjj,…,x′nj)T
wherein, x'jJ column j to n' th element set of x, ajjElement of j row and j column of A, ajIs the 1 st to j-1 st element set, a 'of the j column of A'jJ to n' th element set of column j of A, ujThe element sets from 1 st element to j-1 st element in the jth column of the U;
when calculating L, U column j element, first solve the lower trigonometric equation set Ljuj=aj(Ljuj=ajFor concrete solving expressions, and Lx ═ b above is a general expression) to obtain ujThen solve for x'j=a′j-L′jujX'jSelecting principal element to obtain principal element ujj,ujjFor the jth element of selected jth column of U, and calculate Uj=x′j,lj=x′j/ujjSo that the following trigonometric equation set L is solved based on column 1, …, j-1 of Ljuj=ajGiving L, U a value for column j.
In the LU decomposition process, a method of symbol analysis and numerical decomposition separation is adopted, so that the calculation of zero elements in the decomposition process can be effectively avoided.
Step 8, repeatedly and iteratively solving equation set G delta x in the state estimation calculation processk=HTR-1·[z-h]And correcting the state quantity until DeltaxkSatisfying a specified convergence criterion; wherein G is an active information matrix GpOr reactive information matrix Gq,xkFor the calculation of the kth iteration the node voltage amplitude or phase angle, i.e. the state quantity, xk+1=xk+ΔxkH is a nonlinear measurement vector function, z is a measurement vector, and when G is an active information matrix GpR is thenH is HpWhen G is a reactive information matrix GqR is thenH is Hq,Andthe error variance matrix of active measurement and the error variance matrix of reactive measurement are respectively.
Step 9, when Δ xkAnd when the iteration threshold is smaller than the iteration threshold, stopping iterative computation. And calculating the load flow values of all equipment in the power grid according to the node voltage amplitude and the phase angle obtained by state estimation calculation, wherein the load flow values comprise active and reactive values of equipment such as a line, a transformer, a load and a generator, and the reactive values and the assessment indexes of the reactive value of the capacitive reactance device.
The method realizes the rapid calculation of the active and reactive Jacobian matrixes and the information matrix by using a multithread parallel algorithm based on an OpenMP shared memory programming mode, reduces the calculated amount by using a sparse technology through zero elimination operation in the matrix multiplication process, and accelerates the decomposition speed of the information matrix based on non-zero-element symbol analysis and a numerical decomposition method in the factor decomposition process, thereby improving the overall calculation speed of large-scale power grid state estimation.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (5)
1. A method for improving the calculation speed of state estimation of a power system is characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
converting a relational power grid model in a dispatching control system into a hierarchical power grid model;
reading the measurement of the SCADA in the dispatching system, and associating the measurement with an element in the hierarchical power grid model;
determining a network topology, namely a node-branch model, according to the state of the circuit breaker, the state of the disconnecting link and the connection relation of each element;
remote signaling in measurement comprises a breaker state and a disconnecting link state; when the state estimation is calculated based on a new power grid model, network topology analysis is carried out based on the whole network elements; when the state of the circuit breaker and the disconnecting link exceeds the threshold value compared with the last state change number, network topology analysis is carried out based on the whole network elements; when the current state change number of the circuit breaker and the disconnecting link does not exceed a threshold value compared with the last state change number, only carrying out local network topology analysis on the plant station with the changed remote signaling state according to the voltage grade;
node sorting is carried out based on the measurement number associated with the node;
forming a functional Jacobian matrix H based on OpenMP technology according to a node-branch modelpAnd a reactive Jacobian matrix Hq;
Active information matrix G calculated based on sparse matrix technology and OpenMP technologypAnd a reactive information matrix Gq,GpIs a (n-1) × (n-1) -dimensional matrix, GqThe method is characterized in that the method is an n × n-dimensional matrix, wherein n is the number of nodes calculated by a power grid;
to active information matrix GpThe analysis of the row number sets of the non-zero elements of each column of the lower triangular matrix L and the upper triangular matrix U is completed in sequence from the 1 st column to the n-1 st column according to a non-zero element symbol analysis method;
for reactive information matrix GqAccording to non-zero symbolsThe number analysis method completes the analysis of the row number sets of the non-zero elements of each column of the lower triangular matrix L and the upper triangular matrix U from the 1 st column to the nth column in sequence;
to active information matrix GpAccording to the Left Looking LU numerical decomposition calculation analysis method, combining the obtained non-zero row number set to sequentially perform Left Looking LU numerical decomposition from the 1 st row to the n-1 st row;
for reactive information matrix GqAccording to the Left Looking LU numerical decomposition calculation analysis method, combining the obtained non-zero row number set to sequentially carry out Left Looking LU numerical decomposition from the 1 st row to the n th row;
iteratively solving the system of equations G Δ x repeatedly during the state estimation computation processk=HTR-1·[z-h]And correcting the state quantity until DeltaxkSatisfying a specified convergence criterion; wherein G is an active information matrix GpOr reactive information matrix Gq,xkFor the kth iterative calculation the node voltage amplitude or phase angle, xk+1=xk+ΔxkH is a nonlinear measurement vector function, z is a measurement vector, and when G is an active information matrix GpR is thenH is HpWhen G is a reactive information matrix GqR is thenH is Hq,Andrespectively an active measurement error variance matrix and a reactive measurement error variance matrix;
when Δ xkAnd when the iteration threshold is smaller than the iteration threshold, stopping iterative computation.
2. The method for improving the calculation speed of the state estimation of the power system according to claim 1, wherein: the converted hierarchical network model is stored in a hierarchical database, and measurements are also stored in the hierarchical database and associated with elements in the hierarchical network model.
3. The method for improving the calculation speed of the state estimation of the power system according to claim 1, wherein: decoupling an active modification equation and a reactive modification equation by adopting a rapid decomposition state estimation algorithm based on the least square principle, and calculating an active Jacobian matrix H in parallelpAnd a reactive Jacobian matrix Hq。
4. The method for improving the calculation speed of the state estimation of the power system according to claim 1, wherein: according to the formed active Jacobian matrix HpJacobian matrix H without powerqActive power measurement error variance matrixSum reactive power measurement error variance matrixActive information matrix calculation based on sparse matrix technology and OpenMP technologyAnd a reactive information matrix
5. The method for improving the calculation speed of the state estimation of the power system according to claim 1, wherein:
setting a certain column vector of the coefficient matrix A ═ LU as b, and obtaining a non-zero element structure of a column vector x by solving Lx ═ b;
let β be { i | b }i≠0},χ={j|xjNot equal to 0, representing non-zero elements of b and x, respectivelySet of nodes, wherein biIs the ith row element, x, in column vector bjIs the jth row element in the column vector x, where i ∈ [1, n'],j∈[1,n′]When the coefficient matrix A in the state estimation calculation is an active information matrix GpN is n-1, and the coefficient matrix A is a reactive information matrix G in the state estimation calculationq,n′=n;
Suppose that the k' -1 column of L has been computed, and its corresponding directed graph is G (L)k′-1);
L, U column k' is a set of non-zero nodes derived by the following criteria,
wherein lijNot equal to 0 denotes G (L)k′-1) Where there is an edge, x, from node j to node iiIs the ith row element in the column vector x;
when solving for the non-zero structure of x, the zero values produced by the numerical cancellation are ignored, since lij*xjResult in no matter biWhether or not it is 0, xiIs non-zero;
a and L are stored sparsely according to columns, the row number of a certain column of non-zero elements of A is retrieved, the non-zero elements of the column are retrieved from the corresponding column of L through the row number, and a non-zero set of x is obtained through a depth search algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710823171.8A CN107579516B (en) | 2017-09-13 | 2017-09-13 | Method for improving state estimation calculation speed of power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710823171.8A CN107579516B (en) | 2017-09-13 | 2017-09-13 | Method for improving state estimation calculation speed of power system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107579516A CN107579516A (en) | 2018-01-12 |
CN107579516B true CN107579516B (en) | 2020-08-25 |
Family
ID=61033126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710823171.8A Active CN107579516B (en) | 2017-09-13 | 2017-09-13 | Method for improving state estimation calculation speed of power system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107579516B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108985622B (en) * | 2018-07-13 | 2020-08-04 | 清华大学 | Power system sparse matrix parallel solving method and system based on DAG |
CN109274091B (en) * | 2018-10-15 | 2021-03-26 | 同济大学 | Transmission and distribution integrated parallel state estimation method |
CN111062610B (en) * | 2019-12-16 | 2022-07-29 | 国电南瑞科技股份有限公司 | Power system state estimation method and system based on information matrix sparse solution |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1069517A1 (en) * | 1999-07-16 | 2001-01-17 | ABB Power Automation AG | Visual description of voltage states in an electrical switching graph |
CN102427227B (en) * | 2011-10-18 | 2013-08-21 | 清华大学 | Quick correction decoupling power system state estimating method considering zero injection constraint |
-
2017
- 2017-09-13 CN CN201710823171.8A patent/CN107579516B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN107579516A (en) | 2018-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112217202B (en) | Distributed new energy, energy storage and power distribution network planning method considering flexibility investment | |
Lin et al. | Short-term power prediction for photovoltaic power plants using a hybrid improved Kmeans-GRA-Elman model based on multivariate meteorological factors and historical power datasets | |
Duan et al. | A combined short-term wind speed forecasting model based on CNN–RNN and linear regression optimization considering error | |
CN106447098B (en) | Photovoltaic ultra-short-term power prediction method and device | |
CN107579516B (en) | Method for improving state estimation calculation speed of power system | |
CN114004139B (en) | Photovoltaic power generation power prediction method | |
CN107609141B (en) | Method for performing rapid probabilistic modeling on large-scale renewable energy data | |
CN111062610B (en) | Power system state estimation method and system based on information matrix sparse solution | |
CN107436971B (en) | Improved Latin hypercube sampling method suitable for non-positive stereotype correlation control | |
CN110543720B (en) | State estimation method based on SDAE-ELM pseudo-measurement model | |
CN113988477A (en) | Photovoltaic power short-term prediction method and device based on machine learning and storage medium | |
CN105512808A (en) | Power system transient stability assessment method based on big data | |
CN114792156A (en) | Photovoltaic output power prediction method and system based on curve characteristic index clustering | |
CN110866633A (en) | Micro-grid ultra-short term load prediction method based on SVR support vector regression | |
CN111695724A (en) | Wind speed prediction method based on hybrid neural network model | |
CN111898729A (en) | Power transmission line fault reason identification method and system | |
Huang et al. | Short-term load forecasting based on a hybrid neural network and phase space reconstruction | |
Luo et al. | A multi-step prediction method for wind power based on improved TCN to correct cumulative error | |
CN108694475B (en) | Short-time-scale photovoltaic cell power generation capacity prediction method based on hybrid model | |
CN105207255B (en) | A kind of power system peak regulation computational methods suitable for wind power output | |
CN105095659B (en) | Coordinate distributed state estimation method to province based on cloud computing | |
CN103577890A (en) | Clustering point wind power forecasting method based on tuple matching | |
CN116128211A (en) | Wind-light-water combined short-term optimization scheduling method based on wind-light uncertainty prediction scene | |
Chen et al. | Short-term photovoltaic power prediction based on LGBM-XGBoost | |
CN115173421A (en) | Probability optimal power flow calculation method based on progressive Latin hypercube sampling |
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 |