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 PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage 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
技术领域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:
其中:E为数学期望,以及Where: E is the mathematical expectation, and
其中的样本协方差矩阵,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:
式中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
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:
当上式收敛时,可以得到pOAS When the above formula converges, we can get p OAS
利用得到的pOAS得到最终估计出的协方差矩阵的无偏估计为Use the resulting p OAS to get an unbiased estimate of the final estimated covariance matrix for
本发明所述步骤(6)中提取主成分,其过程是对协方差矩阵的无偏估计进行特征值分解Extract principal component in step (6) of the present invention, its process is the unbiased estimate to covariance matrix Do eigenvalue decomposition
其中为协方差矩阵无偏估计的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
其中S<L,yb满足Where S<L, y b satisfies
利用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
进一步可得线性回归系数vi Further linear regression coefficient v i can be obtained
其中上标-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
式中为在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
由于在电力系统中,故障发生时变量的变化幅度不是很大,所以是非常小的,这导致不能准确地判断出故障是否发生,为避免这一问题,因此将除以训练阶段对应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
本发明所述步骤(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:
其中:E为数学期望,以及Where: E is the mathematical expectation, and
其中为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:
式中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
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:
当上式收敛时,可以得到pOAS When the above formula converges, we can get p OAS
利用得到的pOAS得到最终估计出的协方差矩阵的无偏估计为Use the resulting p OAS to get an unbiased estimate of the final estimated covariance matrix for
(6):对协方差矩阵进行特征值分解,提取相应的主成分;(6): Decompose the eigenvalues of the covariance matrix and extract the corresponding principal components;
对协方差矩阵的无偏估计进行特征值分解Unbiased Estimation of the Covariance Matrix Do eigenvalue decomposition
其中为协方差矩阵无偏估计的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
其中S<L,yb满足Where S<L, y b satisfies
利用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
进一步可得线性回归系数vi Further linear regression coefficient v i can be obtained
其中上标-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
式中为在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
由于在电力系统中,故障发生时变量的变化幅度不是很大,所以是非常小的,这导致不能准确地判断出故障是否发生。为避免这一问题,因此将除以训练阶段对应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
(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)
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)
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)
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 |
-
2016
- 2016-01-22 CN CN201610045643.7A patent/CN105699804B/en not_active Expired - Fee Related
Cited By (29)
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 |