CN105699804A - Big data fault detection and positioning method for power distribution network - Google Patents

Big data fault detection and positioning method for power distribution network Download PDF

Info

Publication number
CN105699804A
CN105699804A CN201610045643.7A CN201610045643A CN105699804A CN 105699804 A CN105699804 A CN 105699804A CN 201610045643 A CN201610045643 A CN 201610045643A CN 105699804 A CN105699804 A CN 105699804A
Authority
CN
China
Prior art keywords
matrix
data
distribution network
covariance matrix
pmu
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
Application number
CN201610045643.7A
Other languages
Chinese (zh)
Other versions
CN105699804B (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.)
Jilin University
Original Assignee
Jilin 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 Jilin University filed Critical Jilin University
Priority to CN201610045643.7A priority Critical patent/CN105699804B/en
Publication of CN105699804A publication Critical patent/CN105699804A/en
Application granted granted Critical
Publication of CN105699804B publication Critical patent/CN105699804B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/22Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • Y04S10/52Outage or fault management, e.g. fault detection or location

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供一种配电网大数据故障检测与定位方法,属于智能电网领域。该方法利用同步相量测量单元PMU接收配电网电压采样数据,对各路PMU接收到的电压进行预处理并构造接收矩阵,应用随机矩阵理论获得接收数据协方差矩阵的无偏估计值,特征值分解该协方差矩阵获得相应的主成分,应用主成分分析计算线性回归系数,应用该系数计算相对近似误差,对比相对近似误差与门限值的大小,实现配电网故障检测与定位。本发明所提方法克服了在多个PMU测量下传统协方差矩阵估计有偏的问题,可实现配电网故障的实时检测与定位。

The invention provides a method for detecting and locating big data faults in a distribution network, which belongs to the field of smart grids. This method uses the synchronized phasor measurement unit PMU to receive the voltage sampling data of the distribution network, preprocesses the voltage received by each PMU and constructs the receiving matrix, and uses the random matrix theory to obtain the unbiased estimate of the covariance matrix of the received data. Decompose the covariance matrix to obtain the corresponding principal components, apply the principal component analysis to calculate the linear regression coefficient, apply the coefficient to calculate the relative approximation error, compare the relative approximation error and the threshold value, and realize the fault detection and location of the distribution network. The method proposed by the invention overcomes the problem of biased traditional covariance matrix estimation under multiple PMU measurements, and can realize real-time detection and location of distribution network faults.

Description

一种配电网大数据故障检测与定位方法A big data fault detection and location method for distribution network

技术领域technical field

本发明属于智能电网领域,具体涉及一种配电网大数据故障检测与定位方法。The invention belongs to the field of smart grids, and in particular relates to a method for detecting and locating big data faults in distribution networks.

背景技术Background technique

配电网是向最终用户供电的网络,是电力系统发电、输电、配电、用电中向用户供电的主要环节。配电网处于电力系统的末端,地域分布广、电网规模大、设备种类多、网络连接多样、运行方式多变等使其运行状态分析具有一定的困难和复杂性。在智能电网系统中,随着同步测量装置(PhasorMeasurementUnit,PMU)广泛使用,电网运行数据的实时处理问题引起了广泛关注。The distribution network is a network that supplies power to end users, and is the main link in power generation, transmission, distribution, and consumption of power systems to supply power to users. The distribution network is at the end of the power system, and its wide geographical distribution, large-scale power grid, various types of equipment, various network connections, and variable operation modes make the analysis of its operating status difficult and complex. In the smart grid system, with the widespread use of the synchronous measurement unit (PhasorMeasurementUnit, PMU), the real-time processing of grid operation data has attracted widespread attention.

应用同步测量装置数据分析电网的监测、监测和控制已经取得显著进展。Lyapunov指数的电压相量被利用来监视短期电压稳定的稳定性,利用离散傅里叶变换可以实现PMU的自适应技术传输线路故障检测和定位,此外,应用相量角度测量和系统拓扑结构可实现配电网故障检测。Significant progress has been made in the monitoring, monitoring and control of power grids using data analysis from synchronized measuring devices. The voltage phasor of the Lyapunov exponent is exploited to monitor the stability of the short-term voltage stability, the adaptive technique of the PMU can be realized by using the discrete Fourier transform to detect and locate the transmission line fault, moreover, the application of the phasor angle measurement and the system topology can realize Distribution network fault detection.

从研究的角度来看,同步测量相量的大型部署带来的一个主要问题是接收数据矩阵维数的增加,即PMU数目和采样点数以不变的比例系数同时趋于无穷大。在现有的基于空间相关性的配电网故障检测与定位方法中,其关键步骤是计算接收数据的协方差矩阵的无偏估计。传统的计算协方差矩阵的方法基于最大似然这一思想,将样本协方差矩阵作为真实协方差矩阵的无偏估计值,然而,这一计算方法的适用条件是PMU数目较小,且样本数远远大于PMU数目。当将上述传统的基于最大似然估计的协方差矩阵计算方法直接应用到大规模部署的PMU测量装置时,将出现较大的估计偏差,导致相应的算法性能明显下降。From a research point of view, a major problem brought about by the large-scale deployment of synchrometric phasors is the increase in the dimension of the received data matrix, that is, the number of PMUs and the number of sampling points tend to infinity at the same time with a constant scaling factor. In the existing distribution network fault detection and location methods based on spatial correlation, the key step is to calculate the unbiased estimate of the covariance matrix of the received data. The traditional method of calculating the covariance matrix is based on the idea of maximum likelihood, and the sample covariance matrix is used as an unbiased estimate of the real covariance matrix. However, the applicable condition of this calculation method is that the number of PMUs is small and the number of samples Much larger than the number of PMUs. When the above-mentioned traditional method of covariance matrix calculation based on maximum likelihood estimation is directly applied to a large-scale deployment of PMU measurement devices, large estimation deviations will occur, resulting in a significant decrease in the performance of the corresponding algorithm.

此外,同步相量测量单元的大规模使用还将明显导致接收数据量的增加,增加了分析数据空间相关性的难度,从数学角度分析,当配电网中出现故障时,必然导致数据空间相关性发生变化,因此,在PMU数目较多的条件下,探索降低接收数据的维度,提升分析数据空间相关性的有效性,是有待解决的另一个问题。In addition, the large-scale use of synchrophasor measurement units will obviously lead to an increase in the amount of received data, which increases the difficulty of analyzing the spatial correlation of data. From a mathematical point of view, when a fault occurs in the distribution network, it will inevitably lead to spatial correlation of data. Therefore, under the condition of a large number of PMUs, it is another problem to be solved to explore and reduce the dimension of the received data and improve the effectiveness of analyzing the spatial correlation of the data.

随机矩阵理论是研究在矩阵的维数和样本数同时趋于无穷时,协方差矩阵的精确估计方法、特征值与特征向量的分布特性等的有效数学工具。本发明所提出的方法在充分分析随机矩阵理论的基础上,研究出适用于多PMU测量下下的电压接收数据的协方差矩阵的无偏估计方法,病讲主成分分析应用到配电网电压故障检测与定位中,以提升算法的实时性和有效性。Random matrix theory is an effective mathematical tool to study the precise estimation method of covariance matrix, the distribution characteristics of eigenvalues and eigenvectors when the dimension of the matrix and the number of samples tend to infinity at the same time. The method proposed in the present invention, on the basis of fully analyzing the random matrix theory, develops an unbiased estimation method for the covariance matrix of the voltage receiving data under multi-PMU measurement, and applies principal component analysis to the distribution network voltage Fault detection and location to improve the real-time and effectiveness of the algorithm.

发明内容Contents of the invention

本发明提供了一种配电网大数据故障检测与定位方法。用于解决在PMU数目较多的情况下,接收数据协方差矩阵估计不准确,数据维数增加限制算法实用性、难以对配电网故障进行实时诊断和定位的问题。The invention provides a big data fault detection and location method of distribution network. It is used to solve the problems that in the case of a large number of PMUs, the estimation of the covariance matrix of the received data is inaccurate, the increase of the data dimension limits the practicability of the algorithm, and it is difficult to diagnose and locate the faults of the distribution network in real time.

本发明采取的方案是,包含如下步骤:The scheme that the present invention takes is, comprises the following steps:

(1)根据配电网拓扑结构布放PMU,形成配电网大数据测量装置;(1) Arrange PMUs according to the topological structure of the distribution network to form a big data measurement device for the distribution network;

(2)应用步骤(1)中的测量装置接收配电网电压数据;(2) The measuring device in the application step (1) receives the distribution network voltage data;

(3)对接收的电压数据进行预处理;(3) Preprocessing the received voltage data;

(4)应用预处理后的数据构造数据接收矩阵;(4) Applying the preprocessed data to construct a data receiving matrix;

(5)应用随机矩阵理论计算数据接收矩阵的协方差矩阵的无偏估计;(5) Apply random matrix theory to calculate the unbiased estimate of the covariance matrix of the data receiving matrix;

(6)对协方差矩阵的无偏估计进行特征值分解,提取相应的主成分;(6) Decompose the eigenvalues of the unbiased estimate of the covariance matrix, and extract the corresponding principal components;

(7)应用主成分分析计算线性回归系数;(7) Apply principal component analysis to calculate the linear regression coefficient;

(8)应用线性回归系数计算相对近似误差;(8) Applying the linear regression coefficient to calculate the relative approximation error;

(9)对比相对近似误差与门限值的大小,实现配电网故障检测与定位。(9) Comparing the relative approximation error and the threshold value to realize the detection and location of distribution network faults.

上述方法中描述的接收数据预处理过程主要为计算接收的电压数据的均值,以此为基础从接收数据中减去均值;The received data preprocessing process described in the above method is mainly to calculate the average value of the received voltage data, and subtract the average value from the received data on this basis;

本发明所述步骤(4)构造数据接收矩阵,其具体构造方法为每一个PMU接收的数据构成接收矩阵的列,每一时刻的采样数据构成接收矩阵的行,假设PMU个数为N,对N路PMU数据进行M次采样,第i(i=0,...,N)个PMU接收到的电压采样数据表示为其中上标T为矩阵的转置,则接收数据矩阵为YM×N:=[y(1),.y(2),..,y(N-1),y(N)]。Step (4) of the present invention constructs the data reception matrix, and its specific construction method is that the data that each PMU receives constitutes the row of the reception matrix, and the sampling data at each moment constitutes the row of the reception matrix, assuming that the number of PMUs is N, for N channels of PMU data are sampled M times, and the voltage sampling data received by the i (i=0,...,N)th PMU is expressed as Wherein the superscript T is the transposition of the matrix, then the received data matrix is Y M×N :=[y (1) ,.y (2) ,..,y (N-1) ,y (N) ].

本发明所述步骤(5)应用随机矩阵理论计算接收电压数据矩阵的协方差矩阵的无偏估计,其具体过程是在OAS估计中,通过权衡低偏差和低方差得到估计函数定义如下:Step (5) of the present invention applies random matrix theory to calculate the unbiased estimate of the covariance matrix of the received voltage data matrix, and its specific process is to obtain the estimated function by weighing low deviation and low variance in OAS estimation It is defined as follows:

minmin pp EE. {{ || || ΣΣ ~~ -- RR YY || || Ff 22 }}

其中:E为数学期望,以及Where: E is the mathematical expectation, and

ΣΣ ~~ == (( 11 -- pp )) RR ~~ YY ++ pp Uu ~~

其中的样本协方差矩阵,p为收缩因子,用于减小均方误差,通常在0与1之间取值,收缩目标U定义如下:in The sample covariance matrix of , p is the shrinkage factor, used to reduce the mean square error, usually takes a value between 0 and 1, and the shrinkage target U is defined as follows:

Uu ~~ == TT rr (( RR ~~ YY )) NN II

式中Tr为矩阵的迹,I为N维单位矩阵,假设各个样本之间为独立同分布,则In the formula, Tr is the trace of the matrix, I is the N-dimensional unit matrix, assuming that the samples are independent and identically distributed, then

pp pp == (( 11 -- 22 NN )) TT rr (( RR YY 22 )) ++ TrTr 22 (( RR YY )) (( Mm ++ 11 -- 22 NN )) TT rr (( RR YY 22 )) ++ (( 11 -- Mm NN )) TrTr 22 (( RR YY ))

RY为真实协方差矩阵,在实际应用中,直接求解真实协方差矩阵RY是不可行的,通过OAS估计对上式进行迭代,得到近似协方差矩阵首次迭代时使用RY的一个初始假设值,p0可以取0到1之间的任意值,之后通过不断迭代对求得的协方差矩阵进行修正,直到迭代收敛:R Y is the real covariance matrix. In practical applications, it is not feasible to directly solve the real covariance matrix R Y. The above formula is iterated through OAS estimation to obtain an approximate covariance matrix An initial hypothetical value of RY is used in the first iteration, p 0 can take any value between 0 and 1, and then modify the obtained covariance matrix through continuous iteration until the iteration converges:

pp ii ++ 11 == (( 11 -- 22 NN )) TT rr (( ΣΣ ~~ ii RR ~~ YY )) ++ TrTr 22 (( ΣΣ ~~ ii )) (( Mm ++ 11 -- 22 NN )) TT rr (( ΣΣ ~~ ii RR ~~ YY )) ++ (( 11 -- Mm NN )) TrTr 22 (( ΣΣ ~~ ii ))

ΣΣ ~~ ii ++ 11 == (( 11 -- pp ii ++ 11 )) RR ~~ YY ++ pp ii ++ 11 Uu ~~

当上式收敛时,可以得到pOAS When the above formula converges, we can get p OAS

pp Oo AA SS == mm ii nno {{ (( 11 -- 22 NN )) TT rr (( RR ~~ YY 22 )) ++ TrTr 22 (( RR ~~ YY )) (( Mm ++ 11 -- 22 NN )) [[ TT rr (( RR ~~ YY 22 )) -- TrTr 22 (( RR ~~ YY )) NN ]] ,, 11 }}

利用得到的pOAS得到最终估计出的协方差矩阵的无偏估计Use the resulting p OAS to get an unbiased estimate of the final estimated covariance matrix for

RR ^^ YY == (( 11 -- pp Oo AA SS )) RR ~~ YY ++ pp Oo AA SS Uu ~~ ..

本发明所述步骤(6)中提取主成分,其过程是对协方差矩阵的无偏估计进行特征值分解Extract principal component in step (6) of the present invention, its process is the unbiased estimate to covariance matrix Do eigenvalue decomposition

11 Mm RR ^^ YY αα jj == ll jj αα jj ,, jj == 11 ,, 22 ,, ...... ,, NN

其中为协方差矩阵无偏估计的N个特征值,为与特征值对应的N个特征向量;in unbiased estimate for the covariance matrix The N eigenvalues of is the N eigenvectors corresponding to the eigenvalues;

的N个特征值进行降序排列,其对应的特征向量随特征值排序;将τ设置为方差阈值,利用累计方差贡献率选择出刚好满足的前L个特征向量,构造L维主成分PCm。Will The N eigenvalues of the eigenvalues are sorted in descending order, and the corresponding eigenvectors are sorted with the eigenvalues; set τ as the variance threshold, and use the cumulative variance contribution rate to select the one that just satisfies The first L eigenvectors of , construct the L-dimensional principal component PCm.

本发明所述步骤(7)中计算线性回归系数,其方法是将矩阵YM×N在L维主成分子空间进行投影,得到矩阵YYCalculate the linear regression coefficient in the step (7) of the present invention, its method is that matrix Y M * N is carried out projection in L dimension principal component subspace, obtains matrix YY

YY=Y·PCmYY=Y·PCm

YY BB :: == [[ ythe y bb 11 ,, ythe y bb 22 ,, ...... ,, ythe y bb SS ]] ∈∈ RR Mm ×× SS

其中S<L,yb满足Where S<L, y b satisfies

cc oo sthe s &theta;&theta; == (( ythe y bb ii &CenterDot;&Center Dot; ythe y bb jj )) // (( || ythe y bb ii || &CenterDot;&CenterDot; || ythe y bb jj || )) &ap;&ap; 00 ,, ii ,, jj == 11 ,, 22 ,, ...... ,, SS

利用YB表示yi,其中yi∈YM×N则线性回归系数vi满足Use Y B to denote y i , where y i ∈ Y M×N and Then the linear regression coefficient v i satisfies

ythe y ii &ap;&ap; &Sigma;&Sigma; jj == 11 mm vv jj ii &CenterDot;&CenterDot; ythe y bb jj == YY BB &CenterDot;&CenterDot; vv ii

进一步可得线性回归系数vi Further linear regression coefficient v i can be obtained

vv ii == (( YY BB TT YY BB )) -- 11 YY BB TT ythe y ii

其中上标-1为矩阵的逆。where the superscript -1 is the inverse of the matrix.

本发明所述步骤(8)中计算相对近似误差,其方法是假设检测阶段近似的第n(n=1,2,…,N-S)个PMU在t时刻的数据为 Calculating the relative approximation error in the step (8) of the present invention, its method is to assume that the data of the nth (n=1, 2,..., NS) PMU that detects stage approximation at time t is

xx ~~ (( tt )) nno :: == YY BB ~~ (( tt )) &CenterDot;&CenterDot; vv ii

式中为在t时刻计算得到的YB,将近似得到的与t时刻的实际电压采样数据z(t)n做差,得到绝对近似误差 In the formula is the Y B calculated at time t, and the approximate Make a difference with the actual voltage sampling data z(t) n at time t to get the absolute approximate error

ee ^^ (( tt )) nno == xx ~~ (( tt )) nno -- zz (( tt )) nno

由于在电力系统中,故障发生时变量的变化幅度不是很大,所以是非常小的,这导致不能准确地判断出故障是否发生,为避免这一问题,因此将除以训练阶段对应PMU的的平均值ers n,可以将故障时刻的近似误差数据放大,做到相对近似误差r(t)n Because in the power system, the variable range is not very large when a fault occurs, so is very small, which makes it impossible to accurately determine whether a fault has occurred. To avoid this problem, the Divided by the training phase corresponding to the PMU The average value of er s n can amplify the approximate error data at the time of the fault, so that the relative approximate error r(t) n

rr (( tt )) nno == ee ^^ (( tt )) nno erer sthe s nno &times;&times; 100100 %%

本发明所述步骤(9)中实现配电网故障检测与定位,其方法是通过历史事故的PMU数据,确定一个阈值η,当In step (9) of the present invention, realize distribution network fault detection and location, its method is to determine a threshold value n by the PMU data of historical accident, when

|r(t)n|≥η|r(t) n |≥η

时,可判定故障发生,t为故障时间,即表征对故障的检测,n为故障PMU,表征对故障的定位。When , it can be determined that a fault occurs, t is the fault time, which represents the detection of the fault, and n is the faulty PMU, which represents the location of the fault.

本发明的优点是:The advantages of the present invention are:

(1)本发明所提方法研究了随机矩阵理论中多种协方差矩阵无偏估计的计算方法,通过对比分析将OAS重构方法引入到配电网大数据处理中,克服了传统样本协方差矩阵计算方法估计偏差较大的问题,提升了相应算法的适用性;(1) The method proposed in the present invention studies the calculation methods of various covariance matrix unbiased estimates in the random matrix theory, and introduces the OAS reconstruction method into the distribution network big data processing through comparative analysis, which overcomes the traditional sample covariance The problem of large estimation deviation of the matrix calculation method improves the applicability of the corresponding algorithm;

(2)本发明所提方法在获得接收数据协方差矩阵无偏估计值的基础上,应用主成分分析降低接收数据的维数,提升了配电网电压空间相关性分析的有效性以及相应算法的计算有效性;(2) On the basis of obtaining the unbiased estimated value of the covariance matrix of the received data, the proposed method of the present invention applies principal component analysis to reduce the dimension of the received data, and improves the validity of the distribution network voltage spatial correlation analysis and the corresponding algorithm The computational validity of

(3)本发明所提方法将相对近似误差作为评价配电网电压空间相关性变化的标准,通过比较相对近似误差与门限值的关系确定配电网故障是否发生以及发生的时刻和位置,可实现电压故障的实时检测与定位。(3) The method proposed in the present invention regards the relative approximation error as the standard for evaluating the spatial correlation change of the distribution network voltage, and determines whether the distribution network fault occurs and the time and position of the occurrence by comparing the relationship between the relative approximation error and the threshold value, Real-time detection and location of voltage faults can be realized.

附图说明Description of drawings

图1是本发明所采用的配电网电压数据测量装置示意图;Fig. 1 is a schematic diagram of a distribution network voltage data measuring device adopted in the present invention;

图2是本发明的流程图;Fig. 2 is a flow chart of the present invention;

图3是IEEE39节点系统图;Fig. 3 is a diagram of IEEE39 node system;

图4是特征值贡献比选择结果图;Fig. 4 is the selection result diagram of eigenvalue contribution ratio;

图5是配电网电压故障检测结果图;Fig. 5 is a diagram of the detection result of the distribution network voltage fault;

图6是配电网电压故障定位结果图。Figure 6 is a diagram of the distribution network voltage fault location results.

具体实施方式detailed description

包括下列步骤:Include the following steps:

(1):根据配电网拓扑结构布放PMU,形成配电网大数据测量装置;(1): Arrange PMUs according to the topological structure of the distribution network to form a big data measurement device for the distribution network;

依据配电网拓扑结构布放PMU,所形成的测量装置如图1所示,其特点是PMU均匀布放在整个配电网;The PMUs are laid out according to the topological structure of the distribution network, and the resulting measurement device is shown in Figure 1, and its characteristic is that the PMUs are evenly placed in the entire distribution network;

(2):应用步骤(1)中的测量装置接收配电网电压数据;(2): The measuring device in the application step (1) receives the distribution network voltage data;

在所需检测区域安放PMU装置,用来分别接收各个区域的PMU数据,将得到的原始数据输送到本地相位数据集中器中,之后把各个本地数据集中器的数据汇总,输送到公司数据集中器,并进行存储;Place the PMU device in the required detection area to receive the PMU data of each area respectively, and send the obtained raw data to the local phase data concentrator, and then aggregate the data of each local data concentrator and send it to the company data concentrator , and store it;

(3):对接收数据进行预处理;(3): Preprocessing the received data;

计算接收的电压数据u的均值v,以此为基础从接收数据中减去均值,获得预处理后的电压数据y=u-v;Calculate the average value v of the received voltage data u, and subtract the average value from the received data on this basis to obtain the preprocessed voltage data y=u-v;

(4):应用预处理后的数据构造数据接收矩阵;(4): apply the preprocessed data to construct a data receiving matrix;

每一个PMU接收的数据构成接收矩阵的列,每一时刻的采样数据构成接收矩阵的行,假设PMU个数为N,对N路PMU数据进行M次采样,第i(i=0,...,N)个PMU接收到的电压采样数据表示为其中上标T为矩阵的转置,则接收数据矩阵为YM×N:=[y(1),.y(2),..,y(N-1),y(N)];The data received by each PMU constitutes the columns of the receiving matrix, and the sampled data at each moment constitutes the rows of the receiving matrix. Assuming that the number of PMUs is N, the data of N channels of PMUs are sampled M times, and the i (i=0,... ., N) The voltage sampling data received by the PMUs is expressed as Wherein the superscript T is the transposition of the matrix, then the received data matrix is Y M×N :=[y (1) ,.y (2) ,..,y (N-1) ,y (N) ];

(5):应用随机矩阵理论计算数据接收矩阵的协方差矩阵的无偏估计;(5): Apply random matrix theory to calculate the unbiased estimate of the covariance matrix of the data receiving matrix;

在OAS估计中,通过权衡低偏差和低方差得到估计函数定义如下:In OAS estimation, the estimated function is obtained by balancing low bias and low variance It is defined as follows:

minmin pp EE. {{ || || &Sigma;&Sigma; ~~ -- RR YY || || Ff 22 }}

其中:E为数学期望,以及Where: E is the mathematical expectation, and

&Sigma;&Sigma; ~~ == (( 11 -- pp )) RR ~~ YY ++ pp Uu ~~

其中为YM×N的样本协方差矩阵,p为收缩因子,用于减小均方误差,通常在0与1之间取值。收缩目标U定义如下:in is the sample covariance matrix of Y M×N , and p is the shrinkage factor, which is used to reduce the mean square error, and usually takes a value between 0 and 1. The contraction target U is defined as follows:

Uu ~~ == TT rr (( RR ~~ YY )) NN II

式中Tr为矩阵的迹,I为N维单位矩阵,假设各个样本之间为独立同分布,则In the formula, Tr is the trace of the matrix, I is the N-dimensional unit matrix, assuming that the samples are independent and identically distributed, then

pp pp == (( 11 -- 22 NN )) TT rr (( RR YY 22 )) ++ TrTr 22 (( RR YY )) (( Mm ++ 11 -- 22 NN )) TT rr (( RR YY 22 )) ++ (( 11 -- Mm NN )) TrTr 22 (( RR YY ))

RY为真实协方差矩阵,在实际应用中,直接求解真实协方差矩阵RY是不可行的,通过OAS估计对上式进行迭代,得到近似协方差矩阵首次迭代时使用RY的一个初始假设值,p0可以取0到1之间的任意值,之后通过不断迭代对求得的协方差矩阵进行修正,直到迭代收敛:R Y is the real covariance matrix. In practical applications, it is not feasible to directly solve the real covariance matrix R Y. The above formula is iterated through OAS estimation to obtain an approximate covariance matrix An initial hypothetical value of RY is used in the first iteration, p 0 can take any value between 0 and 1, and then modify the obtained covariance matrix through continuous iteration until the iteration converges:

pp ii ++ 11 == (( 11 -- 22 NN )) TT rr (( &Sigma;&Sigma; ~~ ii RR ~~ YY )) ++ TrTr 22 (( &Sigma;&Sigma; ~~ ii )) (( Mm ++ 11 -- 22 NN )) TT rr (( &Sigma;&Sigma; ~~ ii RR ~~ YY )) ++ (( 11 -- Mm NN )) TrTr 22 (( &Sigma;&Sigma; ~~ ii ))

&Sigma;&Sigma; ~~ ii ++ 11 == (( 11 -- pp ii ++ 11 )) RR ~~ YY ++ pp ii ++ 11 Uu ~~

当上式收敛时,可以得到pOAS When the above formula converges, we can get p OAS

pp Oo AA SS == mm ii nno {{ (( 11 -- 22 NN )) TT rr (( RR ~~ YY 22 )) ++ TrTr 22 (( RR ~~ YY )) (( Mm ++ 11 -- 22 NN )) &lsqb;&lsqb; TT rr (( RR ~~ YY 22 )) -- TrTr 22 (( RR ~~ YY )) NN &rsqb;&rsqb; ,, 11 }}

利用得到的pOAS得到最终估计出的协方差矩阵的无偏估计Use the resulting p OAS to get an unbiased estimate of the final estimated covariance matrix for

RR ^^ YY == (( 11 -- pp Oo AA SS )) RR ~~ YY ++ pp Oo AA SS Uu ~~

(6):对协方差矩阵进行特征值分解,提取相应的主成分;(6): Decompose the eigenvalues of the covariance matrix and extract the corresponding principal components;

对协方差矩阵的无偏估计进行特征值分解Unbiased Estimation of the Covariance Matrix Do eigenvalue decomposition

11 Mm RR ^^ YY &alpha;&alpha; jj == ll jj &alpha;&alpha; jj ,, jj == 11 ,, 22 ,, ...... ,, NN

其中为协方差矩阵无偏估计的N个特征值,为与特征值对应的N个特征向量;in unbiased estimate for the covariance matrix The N eigenvalues of is the N eigenvectors corresponding to the eigenvalues;

的N个特征值进行降序排列,其对应的特征向量随特征值排序;将τ设置为方差阈值,利用累计方差贡献率选择出刚好满足的前L个特征向量,构造L维主成分PCm;Will The N eigenvalues of the eigenvalues are sorted in descending order, and the corresponding eigenvectors are sorted with the eigenvalues; set τ as the variance threshold, and use the cumulative variance contribution rate to select the one that just satisfies The first L eigenvectors of , construct the L-dimensional principal component PCm;

(7):应用主成分分析计算线性回归系数;(7): Applying principal component analysis to calculate linear regression coefficients;

将矩阵YM×N在L维主成分子空间进行投影,得到矩阵YYProject the matrix Y M×N in the L-dimensional principal component subspace to obtain the matrix YY

YY=Y·PCmYY=Y·PCm

YY BB :: == &lsqb;&lsqb; ythe y bb 11 ,, ythe y bb 22 ,, ...... ,, ythe y bb SS &rsqb;&rsqb; &Element;&Element; RR Mm &times;&times; SS

其中S<L,yb满足Where S<L, y b satisfies

cc oo sthe s &theta;&theta; == (( ythe y bb ii &CenterDot;&Center Dot; ythe y bb jj )) // (( || ythe y bb ii || &CenterDot;&CenterDot; || ythe y bb jj || )) &ap;&ap; 00 ,, ii ,, jj == 11 ,, 22 ,, ...... ,, SS

利用YB表示yi,其中yi∈YM×N则线性回归系数vi满足Use Y B to denote y i , where y i ∈ Y M×N and Then the linear regression coefficient v i satisfies

ythe y ii &ap;&ap; &Sigma;&Sigma; jj == 11 mm vv jj ii &CenterDot;&Center Dot; ythe y bb jj == YY BB &CenterDot;&Center Dot; vv ii

进一步可得线性回归系数vi Further linear regression coefficient v i can be obtained

vv ii == (( YY BB TT YY BB )) -- 11 YY BB TT ythe y ii

其中上标-1为矩阵的逆;where the superscript -1 is the inverse of the matrix;

(8):应用线性回归系数计算相对近似误差;(8): Applying the linear regression coefficient to calculate the relative approximation error;

假设检测阶段近似的第n(n=1,2,…,N-S)个PMU在t时刻的数据为 Assume that the data of the nth (n=1,2,...,NS) PMU approximated in the detection stage at time t is

xx ~~ (( tt )) nno :: == YY BB ~~ (( tt )) &CenterDot;&CenterDot; vv ii

式中为在t时刻计算得到的YB,将近似得到的与t时刻的实际电压采样数据z(t)n做差,得到绝对近似误差 In the formula is the Y B calculated at time t, and the approximate Make a difference with the actual voltage sampling data z(t) n at time t to get the absolute approximate error

ee ^^ (( tt )) nno == xx ~~ (( tt )) nno -- zz (( tt )) nno

由于在电力系统中,故障发生时变量的变化幅度不是很大,所以是非常小的,这导致不能准确地判断出故障是否发生。为避免这一问题,因此将除以训练阶段对应PMU的的平均值ers n,可以将故障时刻的近似误差数据放大,做到相对近似误差r(t)n Because in the power system, the variable range is not very large when a fault occurs, so is very small, which makes it impossible to accurately determine whether a fault has occurred. To avoid this problem, the Divided by the training phase corresponding to the PMU The average value of er s n can amplify the approximate error data at the time of the fault, so that the relative approximate error r(t) n

rr (( tt )) nno == ee ^^ (( tt )) nno erer sthe s nno &times;&times; 100100 %%

(9):对比相对近似误差与门限值的大小,实现配电网故障检测与定位,(9): Comparing the size of the relative approximation error and the threshold value, to realize the fault detection and location of the distribution network,

通过历史事故的PMU数据,确定一个阈值η,当Through the PMU data of historical accidents, determine a threshold η, when

|r(t)n|≥η|r(t) n |≥η

时,可判定故障发生,t为故障时间,即表征对故障的检测,n为故障PMU,表征对故障的定位。When , it can be determined that a fault occurs, t is the fault time, which represents the detection of the fault, and n is the faulty PMU, which represents the location of the fault.

下面通过一些实验数据分析本发明所提出的配电网大数据故障诊断与定位方法的性能。仿真实验1和仿真实验2所采用的仿真软件为MATLAB软件。The performance of the distribution network big data fault diagnosis and location method proposed by the present invention is analyzed below through some experimental data. The simulation software used in simulation experiment 1 and simulation experiment 2 is MATLAB software.

仿真实验:本实验用以分析特征值的贡献比以及所提方法的有效性,实验过程中采用IEEE39节点模型产生的电压数据形成接收数据矩阵,起节点系统图如图3所示,该模型中假设第9个PMU和第30个PMU为故障位置,故障发生时间为分别为第235个采样点和第274个采样点,特征值贡献币比的仿真结果如图4所示,配电网电压故障检测与定位结果分别如图5和图6所示,分析图4可知,应用39维的数据矩阵进行故障检测与故障定位时,可选择10个主特征值所对应的特征向量作为主成分,分析图5可知,所提出的方法在采样点标号分别为235和274点的相对近似误差明显大于其他采样点,即实现了电压故障的实时检测;分析图6可知,在第9个PMU和第30个PMU时,所提方法计算的相对近似误差明显高于其他PMU,即可实现对电压故障的准确定位。Simulation experiment: This experiment is used to analyze the contribution ratio of eigenvalues and the effectiveness of the proposed method. During the experiment, the voltage data generated by the IEEE39 node model is used to form the receiving data matrix. The starting node system diagram is shown in Figure 3. In this model Assuming that the 9th PMU and the 30th PMU are the fault locations, and the fault occurrence time is the 235th sampling point and the 274th sampling point respectively, the simulation results of the eigenvalue contribution ratio are shown in Figure 4, the distribution network voltage The results of fault detection and location are shown in Fig. 5 and Fig. 6 respectively. Analyzing Fig. 4, it can be seen that when using a 39-dimensional data matrix for fault detection and fault location, the eigenvectors corresponding to the 10 main eigenvalues can be selected as the principal components. Analyzing Figure 5, we can see that the relative approximation errors of the proposed method at the sampling points labeled 235 and 274 are significantly larger than other sampling points, that is, real-time detection of voltage faults is realized; analyzing Figure 6, we can see that at the ninth PMU and the When there are 30 PMUs, the relative approximation error calculated by the proposed method is significantly higher than that of other PMUs, and the accurate location of voltage faults can be realized.

Claims (7)

1.一种配电网大数据故障检测与定位方法,其特征在于包含如下步骤:1. A distribution network big data fault detection and location method, is characterized in that comprising the following steps: (1)根据配电网拓扑结构布放PMU,形成配电网大数据测量装置;(1) Arrange PMUs according to the topological structure of the distribution network to form a big data measurement device for the distribution network; (2)应用步骤(1)中的测量装置接收配电网电压数据;(2) The measuring device in the application step (1) receives the distribution network voltage data; (3)对接收的电压数据进行预处理;(3) Preprocessing the received voltage data; (4)应用预处理后的数据构造数据接收矩阵;(4) Applying the preprocessed data to construct a data receiving matrix; (5)应用随机矩阵理论计算数据接收矩阵的协方差矩阵的无偏估计;(5) Apply random matrix theory to calculate the unbiased estimate of the covariance matrix of the data receiving matrix; (6)对协方差矩阵的无偏估计进行特征值分解,提取相应的主成分;(6) Decompose the eigenvalues of the unbiased estimate of the covariance matrix, and extract the corresponding principal components; (7)应用主成分分析计算线性回归系数;(7) Apply principal component analysis to calculate the linear regression coefficient; (8)应用线性回归系数计算相对近似误差;(8) Applying the linear regression coefficient to calculate the relative approximation error; (9)对比相对近似误差与门限值的大小,实现配电网故障检测与定位。(9) Comparing the relative approximation error and the threshold value to realize the detection and location of distribution network faults. 2.根据权利要求1所述一种配电网大数据故障检测与定位方法,其特征在于:步骤(4)构造数据接收矩阵的具体构造方法为每一个PMU接收的数据构成接收矩阵的列,每一时刻的采样数据构成接收矩阵的行,假设PMU个数为N,对N路PMU数据进行M次采样,第i(i=0,...,N)个PMU接收到的电压采样数据表示为其中上标T为矩阵的转置,则接收数据矩阵为 Y M &times; N : = &lsqb; y ( 1 ) , . y ( 2 ) , .. , y ( N - 1 ) , y ( N ) &rsqb; . 2. according to claim 1, a kind of distribution network big data fault detection and location method is characterized in that: the concrete construction method of step (4) constructing data reception matrix is the row that the data that each PMU receives constitutes reception matrix, The sampling data at each moment constitutes the row of the receiving matrix. Assuming that the number of PMUs is N, the data of N PMUs is sampled M times, and the voltage sampling data received by the i (i=0,...,N) PMU Expressed as Where the superscript T is the transpose of the matrix, then the received data matrix is Y m &times; N : = &lsqb; the y ( 1 ) , . the y ( 2 ) , .. , the y ( N - 1 ) , the y ( N ) &rsqb; . 3.根据权利要求1所述一种配电网大数据故障检测与定位方法,其特征在于:步骤(5)的具体过程是在OAS估计中,通过权衡低偏差和低方差得到估计函数定义如下:3. A kind of distribution network big data fault detection and localization method according to claim 1, is characterized in that: the specific process of step (5) is in OAS estimation, obtains estimation function by weighing low deviation and low variance It is defined as follows: mm ii nno pp EE. {{ || || &Sigma;&Sigma; ~~ -- RR YY || || Ff 22 }} 其中:E为数学期望,以及Where: E is the mathematical expectation, and &Sigma;&Sigma; ~~ == (( 11 -- pp )) RR ~~ YY ++ pp Uu ~~ 其中为YM×N的样本协方差矩阵,p为收缩因子,用于减小均方误差,通常在0与1之间取值。收缩目标U定义如下:in is the sample covariance matrix of Y M×N , and p is the shrinkage factor, which is used to reduce the mean square error, and usually takes a value between 0 and 1. The contraction target U is defined as follows: Uu ~~ == TT rr (( RR ~~ YY )) NN II 式中Tr为矩阵的迹,I为N维单位矩阵,假设各个样本之间为独立同分布,则In the formula, Tr is the trace of the matrix, I is the N-dimensional unit matrix, assuming that the samples are independent and identically distributed, then pp pp == (( 11 -- 22 NN )) TT rr (( RR YY 22 )) ++ TrTr 22 (( RR YY )) (( Mm ++ 11 -- 22 NN )) TT rr (( RR YY 22 )) ++ (( 11 -- Mm NN )) TrTr 22 (( RR YY )) RY为真实协方差矩阵,在实际应用中,直接求解真实协方差矩阵RY是不可行的,通过OAS估计对上式进行迭代,得到近似协方差矩阵首次迭代时使用RY的一个初始假设值,p0可以取0到1之间的任意值,之后通过不断迭代对求得的协方差矩阵进行修正,直到迭代收敛:R Y is the real covariance matrix. In practical applications, it is not feasible to directly solve the real covariance matrix R Y. The above formula is iterated through OAS estimation to obtain an approximate covariance matrix An initial hypothetical value of RY is used in the first iteration, p 0 can take any value between 0 and 1, and then modify the obtained covariance matrix through continuous iteration until the iteration converges: pp ii ++ 11 == (( 11 -- 22 NN )) TT rr (( &Sigma;&Sigma; ~~ ii RR ~~ YY )) ++ TrTr 22 (( &Sigma;&Sigma; ~~ ii )) (( Mm ++ 11 -- 22 NN )) TT rr (( &Sigma;&Sigma; ~~ ii RR ~~ YY )) ++ (( 11 -- Mm NN )) TrTr 22 (( &Sigma;&Sigma; ~~ ii )) &Sigma;&Sigma; ~~ ii ++ 11 == (( 11 -- pp ii ++ 11 )) RR ~~ YY ++ pp ii ++ 11 Uu ~~ 当上式收敛时,可以得到pOAS When the above formula converges, we can get p OAS pp Oo AA SS == mm ii nno {{ (( 11 -- 22 NN )) TT rr (( RR ~~ YY 22 )) ++ TrTr 22 (( RR ~~ YY )) (( Mm ++ 11 -- 22 NN )) &lsqb;&lsqb; TT rr (( RR ~~ YY 22 )) -- TrTr 22 (( RR ~~ YY )) NN &rsqb;&rsqb; ,, 11 }} 利用得到的pOAS得到最终估计出的协方差矩阵的无偏估计为:Use the resulting p OAS to get an unbiased estimate of the final estimated covariance matrix for: RR ^^ YY == (( 11 -- pp Oo AA SS )) RR ~~ YY ++ pp Oo AA SS Uu ~~ .. 4.根据权利要求1一种配电网大数据故障检测与定位方法,其特征在于:所述步骤(6)中提取主成分的过程是对协方差矩阵的无偏估计进行特征值分解:4. according to claim 1 a kind of distribution network big data fault detection and localization method, it is characterized in that: the process of extracting principal components in the described step (6) is the unbiased estimation to covariance matrix Perform an eigenvalue decomposition: 11 Mm RR ^^ YY &alpha;&alpha; jj == ll jj &alpha;&alpha; jj ,, jj == 11 ,, 22 ,, ...... ,, NN 其中为协方差矩阵无偏估计的N个特征值,为与特征值对应的N个特征向量;in unbiased estimate for the covariance matrix The N eigenvalues of is the N eigenvectors corresponding to the eigenvalues; 的N个特征值进行降序排列,其对应的特征向量随特征值排序;将τ设置为方差阈值,利用累计方差贡献率选择出刚好满足的前L个特征向量,构造L维主成分PCm;。Will The N eigenvalues of the eigenvalues are sorted in descending order, and the corresponding eigenvectors are sorted with the eigenvalues; set τ as the variance threshold, and use the cumulative variance contribution rate to select the one that just satisfies The first L eigenvectors of , construct the L-dimensional principal component PCm;. 5.根据权利要求1所述一种配电网大数据故障检测与定位方法,其特征在于:所述步骤(7)中计算线性回归系数的方法是将矩阵YM×N在L维主成分子空间进行投影,得到矩阵YY:5. a kind of distribution network big data fault detection and location method according to claim 1, is characterized in that: the method for calculating linear regression coefficient in the described step (7) is to matrix Y M * N in L dimension main component The molecular space is projected to obtain the matrix YY: YY=Y·PCmYY=Y·PCm YY BB :: == &lsqb;&lsqb; ythe y bb 11 ,, ythe y bb 22 ,, ...... ,, ythe y bb SS &rsqb;&rsqb; &Element;&Element; RR Mm &times;&times; SS 其中S<L,yb满足:Where S<L, y b satisfies: cc oo sthe s &theta;&theta; == (( ythe y bb ii &CenterDot;&CenterDot; ythe y bb jj )) // (( || ythe y bb ii || &CenterDot;&CenterDot; || ythe y bb jj || )) &ap;&ap; 00 ,, ii ,, jj == 11 ,, 22 ,, ...... ,, SS 利用YB表示yi,其中yi∈YM×N则线性回归系数vi满足:Use Y B to denote y i , where y i ∈ Y M×N and Then the linear regression coefficient v i satisfies: ythe y ii &ap;&ap; &Sigma;&Sigma; jj == 11 mm vv jj ii &CenterDot;&Center Dot; ythe y bb jj == YY BB &CenterDot;&CenterDot; vv ii 进一步可得线性回归系数viFurther linear regression coefficient v i can be obtained: vv ii == (( YY BB TT YY BB )) -- 11 YY BB TT ythe y ii 其中上标-1为矩阵的逆。where the superscript -1 is the inverse of the matrix. 6.根据权利要求1所述一种配电网大数据故障检测与定位方法,其特征在于:所述步骤(8)中计算相对近似误差的方法是假设检测阶段近似的第n(n=1,2,…,N-S)个PMU在t时刻的数据为 6. According to claim 1, a method for detecting and locating big data faults in a distribution network is characterized in that: the method for calculating the relative approximation error in the step (8) is to assume that the nth (n=1 ,2,…,NS) PMU data at time t is xx ~~ (( tt )) nno :: == YY BB ~~ (( tt )) &CenterDot;&CenterDot; vv ii 式中为在t时刻计算得到的YB,将近似得到的与t时刻的实际电压采样数据z(t)n做差,得到绝对近似误差 In the formula is the Y B calculated at time t, and the approximate Make a difference with the actual voltage sampling data z(t) n at time t to get the absolute approximate error ee ^^ (( tt )) nno == xx ~~ (( tt )) nno -- zz (( tt )) nno 除以训练阶段对应PMU的的平均值可以将故障时刻的近似误差数据放大,做到相对近似误差r(t)n Will Divided by the training phase corresponding to the PMU average of The approximation error data at the fault moment can be amplified to achieve a relative approximation error r(t) n rr (( tt )) nno == ee ^^ (( tt )) nno erer sthe s nno &times;&times; 100100 %% .. 7.根据权利要求1所述一种配电网大数据故障检测与定位方法,其特征在于:所述步骤(9)中实现配电网故障检测与定位的方法是通过历史事故的PMU数据,确定一个阈值η,当7. according to claim 1, a kind of distribution network big data fault detection and location method is characterized in that: in the described step (9), the method for realizing distribution network fault detection and location is by the PMU data of historical accidents, Determine a threshold η, when |r(t)n|≥η|r(t) n |≥η 时,可判定故障发生,t为故障时间,即表征对故障的检测,n为故障PMU,表征对故障的定位。When , it can be determined that a fault occurs, t is the fault time, which represents the detection of the fault, and n is the faulty PMU, which represents the location of the fault.
CN201610045643.7A 2016-01-22 2016-01-22 A kind of power distribution network big data fault detection and location method Expired - Fee Related CN105699804B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610045643.7A CN105699804B (en) 2016-01-22 2016-01-22 A kind of power distribution network big data fault detection and location method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610045643.7A CN105699804B (en) 2016-01-22 2016-01-22 A kind of power distribution network big data fault detection and location method

Publications (2)

Publication Number Publication Date
CN105699804A true CN105699804A (en) 2016-06-22
CN105699804B CN105699804B (en) 2018-05-15

Family

ID=56228572

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610045643.7A Expired - Fee Related CN105699804B (en) 2016-01-22 2016-01-22 A kind of power distribution network big data fault detection and location method

Country Status (1)

Country Link
CN (1) CN105699804B (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106338666A (en) * 2016-08-29 2017-01-18 上海交通大学 Power grid stability detection method based on free probability theory
CN106356994A (en) * 2016-08-29 2017-01-25 上海交通大学 Grid stability judging method based on grid PMU big data
CN106548021A (en) * 2016-11-03 2017-03-29 华中科技大学 A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix
CN106932686A (en) * 2017-03-21 2017-07-07 国网上海市电力公司 Power network exception method for quick based on minimax method of characteristic
CN106950446A (en) * 2017-03-22 2017-07-14 国网上海市电力公司 The abnormal quick detection of power network based on PCA and localization method
CN107132454A (en) * 2017-05-04 2017-09-05 国网上海市电力公司 The abnormal quick determination method of power network based on random matrix spectral radius method
CN108152675A (en) * 2017-12-21 2018-06-12 华中科技大学 It is determined and fault zone localization method based on the fault moment of Random Matrices Theory
CN108196165A (en) * 2018-01-09 2018-06-22 贵州大学 Power grid abnormal state detection method based on sample covariance matrix maximum eigenvalue
CN108734602A (en) * 2018-04-09 2018-11-02 广东电网有限责任公司 A kind of user's electricity stealing identification system based on higher-dimension random matrix
CN109490712A (en) * 2018-12-11 2019-03-19 吉林大学 A kind of Power System Faults Detection method
CN109507526A (en) * 2018-10-29 2019-03-22 山东大学 Fault Diagnosis Method for Distribution Networks based on synchronous phasor measurement and covariance matrix theory
CN109901022A (en) * 2019-04-08 2019-06-18 山东理工大学 Distribution network regional location method based on synchronous measurement data
CN109975594A (en) * 2019-02-28 2019-07-05 北京交通大学 A kind of phasor principal component analytical method for data compression in synchronized measurement system
CN110137944A (en) * 2019-04-24 2019-08-16 国网山东省电力公司莱芜供电公司 A kind of voltage stability disturbance source locating method based on Random Matrices Theory
CN110632455A (en) * 2019-09-17 2019-12-31 武汉大学 A Fault Detection and Location Method Based on Synchronous Measurement Big Data of Distribution Network
CN110673060A (en) * 2019-09-25 2020-01-10 山东大学 Fault diagnosis method for distribution network based on synchrophasor measurement and random matrix theory
CN111273126A (en) * 2020-01-20 2020-06-12 湖南大学 Power distribution network topology rapid sensing method
CN114167214A (en) * 2021-12-02 2022-03-11 浙江大学 Edge-cloud collaborative fault detection method for low-voltage distribution network based on random matrix theory
CN114827951A (en) * 2022-04-24 2022-07-29 重庆长安汽车股份有限公司 Vehicle network quality analysis method and system based on vehicle terminal and storage medium

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5303112A (en) * 1990-10-26 1994-04-12 S & C Electric Company Fault detection method and apparatus
CN102832617B (en) * 2012-09-03 2014-08-13 武汉大学 Large power grid transient state stabilization analyzing method based on precision pattern discrimination
CN103901317B (en) * 2012-12-24 2017-02-08 天津市电力公司 Power distribution fault studying and judging method based on multisystem collaboration
CN103838959A (en) * 2013-12-18 2014-06-04 国网上海市电力公司 Method for applying partial least squares regression to power distribution network harmonic source positioning and detecting
CN104931848B (en) * 2015-03-01 2018-05-08 华北电力大学(保定) A kind of line fault of electrical power system method of discrimination based on wide area measurement information

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106338666B (en) * 2016-08-29 2018-12-25 上海交通大学 A kind of grid stability detection method based on free probability theory
CN106356994A (en) * 2016-08-29 2017-01-25 上海交通大学 Grid stability judging method based on grid PMU big data
CN106356994B (en) * 2016-08-29 2018-09-11 上海交通大学 A kind of grid stability method of discrimination based on power grid PMU big datas
CN106338666A (en) * 2016-08-29 2017-01-18 上海交通大学 Power grid stability detection method based on free probability theory
CN106548021A (en) * 2016-11-03 2017-03-29 华中科技大学 A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix
CN106548021B (en) * 2016-11-03 2018-01-05 华中科技大学 A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix
CN106932686A (en) * 2017-03-21 2017-07-07 国网上海市电力公司 Power network exception method for quick based on minimax method of characteristic
CN106932686B (en) * 2017-03-21 2019-07-09 国网上海市电力公司 Power grid exception rapid detection method based on minimax method of characteristic
CN106950446A (en) * 2017-03-22 2017-07-14 国网上海市电力公司 The abnormal quick detection of power network based on PCA and localization method
CN107132454A (en) * 2017-05-04 2017-09-05 国网上海市电力公司 The abnormal quick determination method of power network based on random matrix spectral radius method
CN108152675A (en) * 2017-12-21 2018-06-12 华中科技大学 It is determined and fault zone localization method based on the fault moment of Random Matrices Theory
CN108196165A (en) * 2018-01-09 2018-06-22 贵州大学 Power grid abnormal state detection method based on sample covariance matrix maximum eigenvalue
CN108734602A (en) * 2018-04-09 2018-11-02 广东电网有限责任公司 A kind of user's electricity stealing identification system based on higher-dimension random matrix
CN108734602B (en) * 2018-04-09 2021-08-03 广东电网有限责任公司 User electricity stealing behavior identification system based on high-dimensional random matrix
CN109507526A (en) * 2018-10-29 2019-03-22 山东大学 Fault Diagnosis Method for Distribution Networks based on synchronous phasor measurement and covariance matrix theory
CN109507526B (en) * 2018-10-29 2020-01-07 山东大学 Fault diagnosis method for distribution network based on synchrophasor measurement and covariance matrix theory
CN109490712A (en) * 2018-12-11 2019-03-19 吉林大学 A kind of Power System Faults Detection method
CN109975594A (en) * 2019-02-28 2019-07-05 北京交通大学 A kind of phasor principal component analytical method for data compression in synchronized measurement system
CN109975594B (en) * 2019-02-28 2021-11-30 北京交通大学 Phasor principal component analysis method for data compression in synchronous measurement system
CN109901022A (en) * 2019-04-08 2019-06-18 山东理工大学 Distribution network regional location method based on synchronous measurement data
CN109901022B (en) * 2019-04-08 2021-12-28 山东理工大学 Power distribution network area positioning method based on synchronous measurement data
CN110137944A (en) * 2019-04-24 2019-08-16 国网山东省电力公司莱芜供电公司 A kind of voltage stability disturbance source locating method based on Random Matrices Theory
CN110632455A (en) * 2019-09-17 2019-12-31 武汉大学 A Fault Detection and Location Method Based on Synchronous Measurement Big Data of Distribution Network
CN110673060A (en) * 2019-09-25 2020-01-10 山东大学 Fault diagnosis method for distribution network based on synchrophasor measurement and random matrix theory
CN111273126A (en) * 2020-01-20 2020-06-12 湖南大学 Power distribution network topology rapid sensing method
CN111273126B (en) * 2020-01-20 2021-01-29 湖南大学 Power distribution network topology rapid sensing method
CN114167214A (en) * 2021-12-02 2022-03-11 浙江大学 Edge-cloud collaborative fault detection method for low-voltage distribution network based on random matrix theory
CN114167214B (en) * 2021-12-02 2022-09-06 浙江大学 Edge-cloud collaborative fault detection method for low-voltage distribution network based on random matrix theory
CN114827951A (en) * 2022-04-24 2022-07-29 重庆长安汽车股份有限公司 Vehicle network quality analysis method and system based on vehicle terminal and storage medium

Also Published As

Publication number Publication date
CN105699804B (en) 2018-05-15

Similar Documents

Publication Publication Date Title
CN105699804B (en) A kind of power distribution network big data fault detection and location method
CN111900731B (en) A PMU-based Power System State Estimation Performance Evaluation Method
WO2015176565A1 (en) Method for predicting faults in electrical equipment based on multi-dimension time series
CN104330471B (en) Lamb wave time-varying probability model monitoring method for aviation structure damage
CN111679158A (en) Distribution network fault identification method based on similarity of synchronous measurement data
CN104519571A (en) Indoor positioning method based on RSS (Received Signal Strength)
CN101726688A (en) Method for diagnosing multi-data-source information fusion-based power system fault
CN110190600A (en) A three-phase distribution network topology identification method based on AMI measurement neighbor regression
CN104836223A (en) Power grid parameter error and bad data coordinated identification and estimation method
CN107681685A (en) A Probabilistic Power Flow Calculation Method for Power System Considering Nonlinear Correlation of Photovoltaics
CN105548718A (en) System harmonic impedance calculation method based on mixed overall least square method
CN105375484A (en) PMU-based electric power system distributed dynamic-state estimation method
CN109524972B (en) Low-frequency oscillation parameter estimation method based on GSO and SVM algorithms
CN106300345A (en) Based on the low-frequency oscillation parameter identification method improving Prony algorithm
Chandler A spectral method for estimating parameters in rainfall models
CN117039890B (en) Distribution network prediction-assisted interval state estimation method for network attack detection
Huang et al. False phasor data detection under time synchronization attacks: A neural network approach
CN114116832A (en) A data-driven abnormality identification method for distribution network
CN105071388B (en) A kind of State Estimation for Distribution Network based on Maximum-likelihood estimation
Xia et al. Outlier reconstruction based distribution system state estimation using equivalent model of long short-term memory and Metropolis-Hastings sampling
CN107462810A (en) A kind of fault section location method suitable for active power distribution network
CN107025500B (en) Method for identifying key deployment position of active power distribution network measurement equipment
CN102280877A (en) Method for identifying parameter of poor branch of power system through a plurality of measured sections
CN103825270B (en) A kind of power distribution network three-phase state estimates the processing method of Jacobian matrix constant
Wu et al. Real time anomaly detection in wide area monitoring of smart grids

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180515

CF01 Termination of patent right due to non-payment of annual fee