CN108196165A - 基于样本协方差矩阵最大特征值的电网异常状态检测方法 - Google Patents
基于样本协方差矩阵最大特征值的电网异常状态检测方法 Download PDFInfo
- Publication number
- CN108196165A CN108196165A CN201810018458.8A CN201810018458A CN108196165A CN 108196165 A CN108196165 A CN 108196165A CN 201810018458 A CN201810018458 A CN 201810018458A CN 108196165 A CN108196165 A CN 108196165A
- Authority
- CN
- China
- Prior art keywords
- matrix
- sliding window
- power grid
- maximum eigenvalue
- covariance matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
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
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- 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
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于样本协方差矩阵最大特征值的电网异常状态检测方法,包括以下步骤:步骤1:构造数据源矩阵Xs;步骤2:获得滑动窗口矩阵X;步骤3:滑动窗口矩阵X的标准化;步骤4:获得样本协方差矩阵S;步骤5:求取样本协方差矩阵最大特征值;步骤6:电网状态异常越限判别:判断是否该最大特征值大于阈值,若成立,则判定电网发生状态异常,发出警告;否则,当前无异常状态,返回步骤2继续执行状态异常检测流程。本发明解决了传统平均谱半径方法在低信噪比场景下进行电网异常状态潜在的失效问题,同时通过简化计算环节,节省了传统基于随机矩阵理论的电网异常状态检测方法的计算耗时,显著提升了计算效率。
Description
技术领域
本发明属于电网异常检测技术领域,具体涉及一种基于样本协方差矩阵最大特征值(Maximum Eigenvalue of Sample Covariance Matrix,MESCM)的电网异常状态检测方法。
背景技术
以同步相量测量装置(Synchronized Phasor Measurement Units,PMU)为基础的广域测量系统(Wide Area Measurement System,WAMS)日趋成熟,所产生的数据量正呈指数级增长。将大数据技术引入传统电力系统分析中,深入开展基于数据驱动的运行状态数据挖掘、提取、分析与融合的研究,实现电网运行状态的“大数据思维”分析与评价,对于我国“互联网+”智能电网的发展具有重要的理论意义。
大数据本质上是一种方法论或认知论,它认为数据是一种主体的内部机理特征或运行演变规律的外在表象。数据驱动是大数据的核心思想,即主要行为决策的制定主要取决于数据分析,而并非物理建模分析或直觉经验,具体来说,是认知者通过分析所收集到的整体/局部数据的统计特性来认知主体。
当前,主流的数据驱动方法有神经网络、支持向量机(Support Vector Machine,SVM)、聚类法、主成份分析法(Principal Component Analysis,PCA)及随机矩阵理论(Random Matrix Theory,RMT)等。神经网络一个最主要的特点在于其具有大量可调的自由参数,这使得其构建起的模型具有较高的灵活性。但另一方面却缺乏有力的理论指导和支撑,大多数情况下仍过分依赖于经验,带有一定程度的随机性;SVM具有严谨的理论基础,仅仅需要少量的样本,对样本的维数不敏感等优势,但过度依赖核函数等问题成为发展的阻碍;聚类是一个无监督学习过程,因此确定最佳聚类数是一项困难的工作。PCA可以发现和辨别一些事件,且具有一定的冗余和容错性,但结果与训练过程强相关,训练过程不合适(如主元空间选择不当)或者系统中发生未经过训练的事件时,PCA的效果将不尽人意。
随机矩阵理论是一种具有普适性的方法,无需详细物理模型,可以从高维角度认识复杂系统的行为特征。一方面,从基于RMT的电力系统分析应用研究进展角度来看,有文献提出了一种基于RMT的电力系统大数据应用架构,并给出了输电网运行状态异常检测方法。在此基础上,有文献利用增广矩阵,采用平均谱半径(Mean Spectral Radius,MSR)指标,进一步提出了一种配电网运行状态相关性分析方法。有文献采用了相同评价指标,从整体上分析了不同扰动对于电力系统暂态稳定性影响程度和影响范围的分析。继而有文献进一步利用历史数据和实时数据建立了随机矩阵模型,通过数据融合,借助平均谱半径指标,基于IEEE 39母线算例,实现了静态稳定态势评估。但广域测量系统中同步相量测量数据涉及现场测量点多、电磁环境复杂、通信距离远的特点,有较高的信噪比降低风险。而现有基于RMT的电力系统分析应用研究中,鲜有涉及上述方法在低信噪比场景下的适应性。此外,传统平均谱半径方法在低信噪比场景下具有电网异常状态潜在的失效问题,而且传统基于随机矩阵理论的电网异常状态检测方法的计算耗时,计算效率低。
发明内容
本发明的目的在于克服上述缺点,提出一种能解决了传统平均谱半径方法在低信噪比场景下进行电网异常状态潜在的失效问题,同时通过减少计算环节,节省了传统基于随机矩阵理论的电网异常状态检测方法的计算耗时,显著提升了计算效率的基于样本协方差矩阵最大特征值的电网异常状态检测方法。
本发明的一种基于样本协方差矩阵最大特征值(Maximum Eigenvalue of SampleCovariance Matrix,MESCM)的电网异常状态检测方法,包括以下步骤:
步骤1:数据源矩阵Xs的构造。假设一个电网有N≥1个同步相量测量单元接收端,在任意采样时刻ti,所接受的信号可构成一个列向量,如式(1)所示,
xs(ti)=(x1,x2,...,xN)T (1)
将各个采样时刻的数据按照时间顺序排列,可以构成一个二维矩阵,即数据源矩阵Xs,如式(2)所示,
步骤2:滑动窗口矩阵X的获得。采用滑动窗口技术分析量测数据,从数据源矩阵Xs中取得当前采样时刻ti的N×T维的滑动窗口矩阵X,其中N是采样数据的维度,单位:个;T是窗口宽度,单位:个;矩阵行列比c,如式(3)所示,
且满足以下条件,即c∈(0,1],为比值,无单位。
步骤3:滑动窗口矩阵X的标准化。对滑动窗口矩阵X进行归一化处理,得到标准非厄米特矩阵(Non-Hermitian Matrix)如式(4)所示。
式中1≤i≤N,1≤j≤T,xi=(xi,1,xi,2,…,xi,T);σ(xi)分别为xi的均值和标准差;分别为的均值和标准差,且
步骤4:样本协方差矩阵S的获得。求矩阵的样本协方差矩阵S,如式(5)所示,
式中上标H表示复共轭转置。
步骤5:样本协方差矩阵最大特征值λmax的求取。计算样本协方差矩阵S的特征值,并从中筛选出最大特征值作为电网异常状态检测指标λmax。
步骤6:电网状态异常越限判别。判断是否该最大特征值λmax大于阈值γ,若λmax≥γ成立,则判定电网发生状态异常,发出警告。否则,当前无异常状态,令i=i+1,返回步骤2继续执行状态异常检测流程。
其中,步骤6所述阈值的设定,如式(6)所示,
式中c为前述滑动窗口矩阵X的矩阵行列比,kγ为结合异常事件分级告警的思路,设置的阈值裕度,一般为1.2。
其中,步骤1所述采样数据类型的选择为母线电压幅值V、母线电压相角θ以及发电机功角信号δ。
其中,步骤1所述采样数据维数N的选择:各采样数据维数应结合应用需求,以满足所监测区域电力系统的系统可观性为前提,并考虑一定冗余度。从中国大区联合电力系统应用需求来看,采样数据涉及的对象一般宜包括其范围内所有的1000kV变电站1000kV和500kV母线、所有的500kV变电站500kV和220kV母线、所有的接入500kV电压等级电网的大型区域性发电厂500kV母线和机端电压母线、重要的接入220kV电压等级电网的中型局域性发电厂220kV母线和机端电压母线、以及重要的220kV变电站220kV母线。从中国省级电力系统应用需求来看,采样数据涉及的对象一般宜包括其范围内所有的1000kV变电站1000kV和500kV母线、所有的500kV变电站500kV和220kV母线、所有的接入500kV电压等级电网的大型区域性发电厂500kV母线和机端电压母线、所有的接入220kV电压等级电网的中型局域性发电厂220kV母线和机端电压母线、以及所有的220kV变电站220kV母线、重要的接入110kV电压等级电网的中小型局域性发电厂110kV母线和机端电压母线、以及重要的110kV变电站110kV母线。
其中,步骤2所述滑动窗口宽度T的选择,该变量表示滑动窗口内离散采样点的数量。就暂态事件异常检测应用来说,一般所述滑动窗口时间宽度为10~30秒,采样频率100Hz,即T=1000~3000。就稳态事件异常检测应用来说,一般所述滑动窗口时间宽度为300~900秒,采样频率10Hz,即T=3000~9000。
本发明与现有技术相比,具有以下优点:
1、相较于传统的平均谱半径分析法,本发明基于样本协方差矩阵最大特征值的方法抗噪性能更好、计算耗时更少。
2.本发明方法在非完整信息情形下能够一定程度上实现有效检测,具有一定的鲁棒性。
附图说明
图1是本发明的流程图;
图2是实施例中10机39母线电网接线结构示意图;
图3是案例一完整信息下的MESCM和MSR方法的结果比较;
图4是案例一的谱分布示意图;
图5是案例一非完整信息下的MESCM和MSR方法的结果比较;
图6是案例二的MESCM和MSR方法的结果比较;
图7是案例二的谱分布示意图。
具体实施方式
以下结合附图和实施例详细描述本发明的具体实施方式,但本发明不受所述具体实施例所限。
考虑传输过程中信道噪声的影响,定义异常状态检测模型如式(1)所示。
Xs=XP+m×η (1)
其中XP为PMU信号矩阵,η为噪声矩阵,本发明均采用高斯噪声,m为噪声幅值。噪声大小会影响电网异常状态检测的效果。为此,定义数据源矩阵的信噪比(Signal-to-NoiseRatio,SNR)如式(2)所示。
其中Tr(·)为矩阵的迹。
依据本发明步骤,如图1所示,
步骤1:数据源矩阵Xs的构造。假设一个电网有N≥1个同步相量测量单元接收端,在任意采样时刻ti,所接受的信号可构成一个列向量,如式(3)所示,
xs(ti)=(x1,x2,...,xN)T (3)
将各个采样时刻的数据按照时间顺序排列,可以构成一个二维矩阵,即数据源矩阵Xs,如式(4)所示,
步骤2:滑动窗口矩阵X的获得。采用滑动窗口技术分析量测数据,从数据源矩阵Xs中取得当前采样时刻ti的N×T维的滑动窗口矩阵X,其中N是采样数据的维度,单位:个;T是窗口宽度,单位:个;矩阵行列比c,如式(5)所示,
且满足以下条件,即c∈(0,1],为比值,无单位。
步骤3:滑动窗口矩阵X的标准化。对滑动窗口矩阵X进行归一化处理,得到标准非厄米特矩阵(Non-Hermitian Matrix)如式(6)所示。
式中1≤i≤N,1≤j≤T,xi=(xi,1,xi,2,…,xi,T);σ(xi)分别为xi的均值和标准差;分别为的均值和标准差,且
步骤4:样本协方差矩阵S的获得。求矩阵的样本协方差矩阵S,如式(7)所示。
式中上标H表示复共轭转置;
值得注意的是,样本协方差矩阵S的经验谱分布函数(Empirical SpectralDistribution,ESD)服从M-P律,如式(8)和(9):
其中
式中a和b分别表示谱密度函数中的最小特征值和最大特征值,也就是说S的特征值分布具有一定范围,是有界的。当有异常事件发生时系统的随机性被破坏,将不满足统计规律,最大特征值会超过其边界范围,本发明方法正是利用了这一特点进行异常状态的检测。
步骤5:样本协方差矩阵最大特征值λmax的求取。计算样本协方差矩阵S的特征值,并从中筛选出最大特征值作为电网异常状态检测指标λmax;
步骤6:电网状态异常越限判别。判断是否该最大特征值λmax大于阈值γ,若λmax≥γ成立,则判定电网发生状态异常,发出警告。否则,当前无异常状态,令i=i+1,返回步骤2继续执行状态异常检测流程。
其中,步骤6所述阈值的设定,如式(10)所示,
式中c为前述滑动窗口矩阵X的矩阵行列比,kγ为结合异常事件分级告警的思路,设置的阈值裕度,一般为1.2。
其中,步骤1所述采样数据类型的选择为母线电压幅值V、母线电压相角θ以及发电机功角信号δ。
其中,步骤1所述采样数据维数N的选择:各采样数据维数应结合应用需求,以满足所监测区域电力系统的系统可观性为前提,并考虑一定冗余度。从中国大区联合电力系统应用需求来看,采样数据涉及的对象一般宜包括其范围内所有的1000kV变电站1000kV和500kV母线、所有的500kV变电站500kV和220kV母线、所有的接入500kV电压等级电网的大型区域性发电厂500kV母线和机端电压母线、重要的接入220kV电压等级电网的中型局域性发电厂220kV母线和机端电压母线、以及重要的220kV变电站220kV母线。从中国省级电力系统应用需求来看,采样数据涉及的对象一般宜包括其范围内所有的1000kV变电站1000kV和500kV母线、所有的500kV变电站500kV和220kV母线、所有的接入500kV电压等级电网的大型区域性发电厂500kV母线和机端电压母线、所有的接入220kV电压等级电网的中型局域性发电厂220kV母线和机端电压母线、以及所有的220kV变电站220kV母线、重要的接入110kV电压等级电网的中小型局域性发电厂110kV母线和机端电压母线、以及重要的110kV变电站110kV母线。。
其中,步骤2所述滑动窗口宽度T的选择,该变量表示滑动窗口内离散采样点的数量。滑动窗口的宽度与采样频率及去噪性能强相关。较大的滑动窗口宽度,将增加获取的信号数据,去噪能力将增强,但会掩盖较为微弱的异常信号;较小的滑动窗口宽度,将对异常信号更加敏感,但在采样频率高的情况下,原始信号中的信息将被破坏。故应结合采样频率和状态变量的数量来确定窗口的宽度。就暂态事件异常检测应用来说,一般所述滑动窗口时间宽度为10~30秒,采样频率100Hz,即T=1000~3000。就稳态事件异常检测应用来说,一般所述滑动窗口时间宽度为300~900秒,采样频率10Hz,即T=3000~9000。。
为了验证本文所提方法在异常负荷跃变以及三相接地故障下的识别能力,算例分析在一个IEEE39节点标准系统和一个南方电网规划系统中开展。具体为:在IEEE39节点系统中分别设置了节点15的负荷异常跃变及节点15至节点16线路的三相接地短路故障;在南方电网规划系统中设置贵州区域某220kV联络输电线路三相短路接地故障。借助PSS/E34.2软件开展时域仿真,以获取模拟的PMU数据,进而构造数据源矩阵,仿真步长为0.01s。依据图1所示步骤,利用MATLAB2014a软件编制算法程序,通过与传统的平均谱半径分析方法计算结果的比较,验证其有效性以及低信噪比场景下的适应性。
案例一:一个IEEE 10机39母线系统
1、完整信息下异常负荷跃变测试
分别在ρ=4.2×103的高信噪比场景及ρ=260的低信噪比场景下进行研究,共2500个采样点,其中采样时刻t 0至t500中无异常事件发生,设置从采样时刻t501起母线15有功功率由320MW跃变至360MW,采样时刻t1500后恢复正常。
IEEE 10机39母线接线结构图,如图2所示。选取除平衡母线39外其他38个母线中负荷母线的电压相角和发电机母线的功角数据,构成38维数据源矩阵进行分析。设滑动窗口T=120,即一个窗口内需要120组采样数据(包括1组当前时刻和119组历史数据),因此本发明方法获得的MESCM指标和平均谱半径分析法获得的MSR指标变化曲线从采样时刻t120开始。该测试中阈值设定为γ=2.9。
依次对每个滑动时间窗口构成的矩阵,按照图1所述步骤进行计算,可得到ρ=4.2×103的高信噪比场景及ρ=260的低信噪比场景下的MESCM指标变化曲线,如图3(a)所示,同时,由传统平均谱半径分析方法,可获得MSR指标变化曲线,如图3(b)所示。图3(a)和(b)中,粗虚线分别代表MESCM和MSR指标的阈值或内径,实线分别代表为高信噪比环境下MESCM和MSR值,细虚线分别代表低信噪比环境下MESCM和MSR值。
由该预想扰动事件描述,可知在采样时刻t501及t1500,系统分别有母线负荷阶跃上升和阶跃下降行为发生。对比图3(a)和(b)可以发现,一方面,在高信噪比场景下,系统随机性被打破,MESCM及MSR均明显越过了各自阈值,说明了两种方法能够有效检测电网状态异常。另一方面,在低信噪比场景下,由于噪声的干扰,MSR指标并未越过其内径,不能实现有效检测。而此时MESCM的最大值小于高信噪比下的最大值,在采样时刻t526及t1527突破其阈值,表征MESCM法能够检测出状态异常事件的发生。
此外,从上述算例的谱分布角度来看,可绘制两种信噪比场景下特征值分布、平均谱半径以及内外径情况如图4所示。图中深色长虚线为MSR外径,深色短虚线为MSR内径,上三角和下三角分别代表高信噪比和低信噪比环境下的特征值,浅色短虚线和浅色长虚线分别代表高信噪比和低信噪比环境下的MSR。
从图4可以看出,高信噪比环境下,采样时刻t 1476,MSR跌落至最低值0.776,说明MSR随着事件的发生而偏离了内径限值,或者说特征值分布出现了“坍缩”现象,证明了该方法对于状态异常检测的有效性。而低信噪比环境下,采样时刻t1472,MSR跌落至最低值0.8816。此时,特征值分布仍符合单环定律,但并未越限,表明该方法在低信噪比环境下难以有效检测异常负荷跃变。
2、非完整信息下三相短路故障测试
仍然采用与前述相同的2种信噪比环境设定,设置采样时刻t1000母线15至母线16线路首端发生三相接地短路故障,5个采样时刻后切除故障线路,共2000个采样时刻。
为验证非完整信息下的适应性,剔除15、16母线及与之相连的14、18、17、21、24等这几个母线的信息。选取除平衡母线39外的其他31个母线中负荷母线的电压,电压相角;发电机母线的电压,功角等数据构成62维数据源矩阵。设滑动窗口T=200,阈值γ=3.01,开展所提MESCM方法分析,结果如图5(a)所示,类似地给出MSR法计算结果如图5(b)所示。
由图5中(a)与(b)的实线可以发现,高信噪比环境中非完整信息情形下,两种方法均能够有效的识别出三相短路故障呈现的异常状态。然而,对比图5中(a)与(b)的粗虚线可以发现,低信噪比环境中非完整信息情形下,MSR在采样时刻t 1021达到最小值0.8746,但该值仍然大于其内径值,说明此时该方法不能有效检测异常状态,而MESCM方法则能够有效识别。
案例二:一个实际电网系统
1、大规模电网的适用性分析
以一个南方电网1592机7877母线规划系统为对象,考虑ρ=4.4×10-3的低信噪比场景,采样时刻t1000设置贵州区域内某220kV输电线路三相短路故障,5个采样时刻后切除故障线路,共2000个采样点。取各省级区域共200个主要500kV母线的电压相角构造200维数据源矩阵,设滑动窗口T=300,阈值γ=3.96,由所提MESCM方法可得结果如图6(a)所示,类似地给出MSR法计算结果如图6(b)所示。此外,从谱分布角度,给出相应结果如图6所示。
观察图6(a)可以发现,受噪声影响,MESCM从采样时刻t1078开始出现较为明显的上升趋势,延迟于该局部区域220kV三相短路故障事件发生后78个采样时刻,在采样时刻t1098越过阈值并发出异常状态报警。而从图6(b)可以看出,MSR在采样时刻t1228跌落至最小值0.79,仍大于其内径,难以实现有效检测。图7谱分布情况呈现出200个特征值点仍收敛于双环之内,表明了MSR分析法易受噪声影响而失效。
2、计算耗时分析
利用MATLABR2014a软件,在主频3.2GHz CPU,8G RAM的计算机上进行了MESCM法和MSR法的计算效率比较分析。将上述三个算例涉及的不同维数规模的数据源矩阵所需计算耗时列于表1中。
表1两种方法的计算耗时比较
由表1可知,两种算法的计算耗时均随着数据源矩阵维数增加而上升。但MESCM法的耗时相对较少,大约仅需传统MSR法的1/5。
综上所述,实施例中2个不同规模电网算例验证了本发明提出的一种基于MESCM的电网异常状态检测方法的有效性,相较于以传统的平均谱半径分析法为代表的现有技术,本发明基于样本协方差矩阵最大特征值的方法抗噪性能更好、计算耗时更少,在电力系统态势感知技术领域,具有明显的新颖性和创造性。同时本发明方法在非完整信息情形下能够实现检测,具有一定的鲁棒性,具有较好的实用性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,任何未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (8)
1.一种基于样本协方差矩阵最大特征值的电网异常状态检测方法,包括以下步骤:
步骤1:构造数据源矩阵Xs:利用电网多个同步相量测量单元接收端,将各个采样时刻的数据按照时间顺序排列,构成一个多维矩阵,即数据源矩阵Xs;
步骤2:获得滑动窗口矩阵X:采用滑动窗口技术分析量测数据,从数据源矩阵Xs中取得当前采样时刻的滑动窗口矩阵X;
步骤3:滑动窗口矩阵X的标准化:对滑动窗口矩阵X进行归一化处理;
步骤4:获得样本协方差矩阵S:对归一化处理的滑动窗口矩阵,求样本协方差矩阵S;
步骤5:求取样本协方差矩阵最大特征值:计算样本协方差矩阵S的特征值,并从中筛选出最大特征值作为电网异常状态检测指标;
步骤6:电网状态异常越限判别:判断是否该最大特征值大于阈值,若成立,则判定电网发生状态异常,发出警告;否则,当前无异常状态,返回步骤2继续执行状态异常检测流程。
2.如权利要求1所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述步骤1中的构造数据源矩阵Xs:利用电网多个同步相量测量单元接收端,在任意采样时刻t i,所接受的信号可构成一个列向量,如式(1)所示,
将各个采样时刻的数据按照时间顺序排列,可以构成一个二维矩阵,即数据源矩阵Xs,如式(2)所示,
。
3.如权利要求1所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述步骤2中的获得滑动窗口矩阵X:采用滑动窗口技术分析量测数据,从数据源矩阵Xs中取得当前采样时刻t i的N×T维的滑动窗口矩阵X,其中N是采样数据的维数,单位:个;T是滑动窗口宽度,单位:个;矩阵行列比c,如式(3)所示,
且满足以下条件,即c∈(0,1],为比值,无单位。
4.如权利要求1所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述步骤3中的滑动窗口矩阵X的标准化:对滑动窗口矩阵X进行归一化处理,得到标准非厄米特矩阵(Non- Hermitian Matrix) ,如式(4)所示,
式中1,;分别为的均值和标准差;、分别为的均值和标准差,且;。
5.如权利要求1所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述步骤4中的求滑动窗口矩阵X归一化处理得到的标准非厄米特矩阵的样本协方差矩阵S,如式(5)所示,
(5)
式中上标H表示复共轭转置。
6.根据权利要求3所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述步骤6中的阈值的设定,如式(6)所示,
式中c为前述滑动窗口矩阵X的矩阵行列比,kγ为结合异常事件分级告警的思路,设置的阈值裕度,一般为1.2。
7.如权利要求1至6中任一项所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述步骤1中采样数据的类型为母线电压幅值V、母线电压相角θ以及发电机功角信号δ。
8.如权利要求6所述的基于样本协方差矩阵最大特征值的电网异常状态检测方法,其特征在于:所述滑动窗口宽度,就暂态事件异常检测应用来说,滑动窗口时间宽度为10~30秒,采样频率100Hz,即T=1000~3000;就稳态事件异常检测应用来说,滑动窗口时间宽度为300~900秒,采样频率10Hz,即T=3000~9000。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810018458.8A CN108196165A (zh) | 2018-01-09 | 2018-01-09 | 基于样本协方差矩阵最大特征值的电网异常状态检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810018458.8A CN108196165A (zh) | 2018-01-09 | 2018-01-09 | 基于样本协方差矩阵最大特征值的电网异常状态检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108196165A true CN108196165A (zh) | 2018-06-22 |
Family
ID=62588841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810018458.8A Pending CN108196165A (zh) | 2018-01-09 | 2018-01-09 | 基于样本协方差矩阵最大特征值的电网异常状态检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108196165A (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108647838A (zh) * | 2018-06-27 | 2018-10-12 | 南京师范大学 | 一种基于随机矩阵理论和在线序列极限学习机的电网态势感知方法 |
CN109061387A (zh) * | 2018-08-23 | 2018-12-21 | 国网内蒙古东部电力有限公司电力科学研究院 | 一种基于圆环定理的电网异常状态判别方法 |
CN109242321A (zh) * | 2018-09-17 | 2019-01-18 | 国网河北省电力有限公司电力科学研究院 | 用户电力负荷在线分析方法及终端设备 |
CN109361484A (zh) * | 2018-11-13 | 2019-02-19 | 上海电器科学研究所(集团)有限公司 | 一种电力系统时间同步数据的传输方法 |
CN109507526A (zh) * | 2018-10-29 | 2019-03-22 | 山东大学 | 基于同步相量测量和协方差矩阵理论的配电网故障诊断方法 |
CN109657613A (zh) * | 2018-12-19 | 2019-04-19 | 贵州大学 | 基于幂法和并行计算技术的大规模电网异常负荷识别方法 |
CN109657943A (zh) * | 2018-12-06 | 2019-04-19 | 中国科学院深圳先进技术研究院 | 风电场机组运行状态的动态评估方法、装置及电子设备 |
CN109800191A (zh) * | 2019-01-25 | 2019-05-24 | 中科驭数(北京)科技有限公司 | 用于计算序列数据的协方差的方法及装置 |
CN109884469A (zh) * | 2019-03-06 | 2019-06-14 | 山东理工大学 | 配电网故障区段与故障时刻的判定方法 |
CN109901022A (zh) * | 2019-04-08 | 2019-06-18 | 山东理工大学 | 基于同步量测数据的配电网区域定位方法 |
CN110045227A (zh) * | 2019-03-23 | 2019-07-23 | 广西电网有限责任公司电力科学研究院 | 一种基于随机矩阵与深度学习的配电网故障诊断方法 |
CN110580492A (zh) * | 2019-04-30 | 2019-12-17 | 上海铁大电信科技股份有限公司 | 一种基于小幅波动检测的轨道电路故障前兆发现方法 |
CN110632455A (zh) * | 2019-09-17 | 2019-12-31 | 武汉大学 | 一种基于配网同步量测大数据的故障检测定位方法 |
CN110988587A (zh) * | 2019-11-13 | 2020-04-10 | 上海恒能泰企业管理有限公司 | 一种基于最大最小特征值法的配网异常检测方法 |
CN111007351A (zh) * | 2019-11-13 | 2020-04-14 | 上海恒能泰企业管理有限公司 | 一种基于高维随机矩阵特征根检测的配网异常检测方法 |
CN112308299A (zh) * | 2020-10-19 | 2021-02-02 | 新奥数能科技有限公司 | 用于电力系统负荷预测模型的样本数据提取方法和装置 |
CN113899963A (zh) * | 2021-08-30 | 2022-01-07 | 科华数据股份有限公司 | 交流系统晃电检测方法、装置及终端设备 |
CN114062850A (zh) * | 2021-11-17 | 2022-02-18 | 江南大学 | 一种双阈值电网早期故障检测方法 |
CN114784973A (zh) * | 2022-05-18 | 2022-07-22 | 苏州纬讯光电科技有限公司 | 一种单芯高压电缆终端故障预警方法 |
WO2023241327A1 (zh) * | 2022-06-14 | 2023-12-21 | 无锡隆玛科技股份有限公司 | 一种最大特征向量电网异常定位方法 |
WO2023241326A1 (zh) * | 2022-06-14 | 2023-12-21 | 无锡隆玛科技股份有限公司 | 基于样本协方差矩阵最大特征值变化率的电网异常检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150346286A1 (en) * | 2012-06-26 | 2015-12-03 | Siemens Aktiengesellschaft | Fault detection in energy supply networks |
CN105699804A (zh) * | 2016-01-22 | 2016-06-22 | 吉林大学 | 一种配电网大数据故障检测与定位方法 |
CN106932686A (zh) * | 2017-03-21 | 2017-07-07 | 国网上海市电力公司 | 基于最大最小特征值法的电网异常快速检测方法 |
CN107015079A (zh) * | 2017-03-28 | 2017-08-04 | 国网上海市电力公司 | 一种基于广域测量的配电系统网络异常检测方法 |
CN107132454A (zh) * | 2017-05-04 | 2017-09-05 | 国网上海市电力公司 | 基于随机矩阵谱半径法的电网异常快速检测方法 |
-
2018
- 2018-01-09 CN CN201810018458.8A patent/CN108196165A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150346286A1 (en) * | 2012-06-26 | 2015-12-03 | Siemens Aktiengesellschaft | Fault detection in energy supply networks |
CN105699804A (zh) * | 2016-01-22 | 2016-06-22 | 吉林大学 | 一种配电网大数据故障检测与定位方法 |
CN106932686A (zh) * | 2017-03-21 | 2017-07-07 | 国网上海市电力公司 | 基于最大最小特征值法的电网异常快速检测方法 |
CN107015079A (zh) * | 2017-03-28 | 2017-08-04 | 国网上海市电力公司 | 一种基于广域测量的配电系统网络异常检测方法 |
CN107132454A (zh) * | 2017-05-04 | 2017-09-05 | 国网上海市电力公司 | 基于随机矩阵谱半径法的电网异常快速检测方法 |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108647838A (zh) * | 2018-06-27 | 2018-10-12 | 南京师范大学 | 一种基于随机矩阵理论和在线序列极限学习机的电网态势感知方法 |
CN108647838B (zh) * | 2018-06-27 | 2022-04-01 | 南京师范大学 | 一种基于随机矩阵理论和在线序列极限学习机的电网态势感知方法 |
CN109061387A (zh) * | 2018-08-23 | 2018-12-21 | 国网内蒙古东部电力有限公司电力科学研究院 | 一种基于圆环定理的电网异常状态判别方法 |
CN109242321A (zh) * | 2018-09-17 | 2019-01-18 | 国网河北省电力有限公司电力科学研究院 | 用户电力负荷在线分析方法及终端设备 |
CN109242321B (zh) * | 2018-09-17 | 2020-10-09 | 国网河北省电力有限公司电力科学研究院 | 用户电力负荷在线分析方法及终端设备 |
CN109507526B (zh) * | 2018-10-29 | 2020-01-07 | 山东大学 | 基于同步相量测量和协方差矩阵理论配电网故障诊断方法 |
CN109507526A (zh) * | 2018-10-29 | 2019-03-22 | 山东大学 | 基于同步相量测量和协方差矩阵理论的配电网故障诊断方法 |
CN109361484A (zh) * | 2018-11-13 | 2019-02-19 | 上海电器科学研究所(集团)有限公司 | 一种电力系统时间同步数据的传输方法 |
CN109657943A (zh) * | 2018-12-06 | 2019-04-19 | 中国科学院深圳先进技术研究院 | 风电场机组运行状态的动态评估方法、装置及电子设备 |
CN109657613A (zh) * | 2018-12-19 | 2019-04-19 | 贵州大学 | 基于幂法和并行计算技术的大规模电网异常负荷识别方法 |
CN109657613B (zh) * | 2018-12-19 | 2023-05-02 | 贵州大学 | 基于幂法和并行计算技术的大规模电网异常负荷识别方法 |
CN109800191A (zh) * | 2019-01-25 | 2019-05-24 | 中科驭数(北京)科技有限公司 | 用于计算序列数据的协方差的方法及装置 |
CN109800191B (zh) * | 2019-01-25 | 2020-04-24 | 中科驭数(北京)科技有限公司 | 用于计算序列数据的协方差的方法及装置 |
CN109884469A (zh) * | 2019-03-06 | 2019-06-14 | 山东理工大学 | 配电网故障区段与故障时刻的判定方法 |
WO2020177706A1 (zh) * | 2019-03-06 | 2020-09-10 | 山东理工大学 | 配电网故障区段与故障时刻的判定方法 |
CN110045227B (zh) * | 2019-03-23 | 2019-12-17 | 广西电网有限责任公司电力科学研究院 | 一种基于随机矩阵与深度学习的配电网故障诊断方法 |
CN110045227A (zh) * | 2019-03-23 | 2019-07-23 | 广西电网有限责任公司电力科学研究院 | 一种基于随机矩阵与深度学习的配电网故障诊断方法 |
CN109901022B (zh) * | 2019-04-08 | 2021-12-28 | 山东理工大学 | 基于同步量测数据的配电网区域定位方法 |
CN109901022A (zh) * | 2019-04-08 | 2019-06-18 | 山东理工大学 | 基于同步量测数据的配电网区域定位方法 |
CN110580492A (zh) * | 2019-04-30 | 2019-12-17 | 上海铁大电信科技股份有限公司 | 一种基于小幅波动检测的轨道电路故障前兆发现方法 |
CN110632455A (zh) * | 2019-09-17 | 2019-12-31 | 武汉大学 | 一种基于配网同步量测大数据的故障检测定位方法 |
CN111007351A (zh) * | 2019-11-13 | 2020-04-14 | 上海恒能泰企业管理有限公司 | 一种基于高维随机矩阵特征根检测的配网异常检测方法 |
CN110988587A (zh) * | 2019-11-13 | 2020-04-10 | 上海恒能泰企业管理有限公司 | 一种基于最大最小特征值法的配网异常检测方法 |
CN112308299A (zh) * | 2020-10-19 | 2021-02-02 | 新奥数能科技有限公司 | 用于电力系统负荷预测模型的样本数据提取方法和装置 |
CN112308299B (zh) * | 2020-10-19 | 2024-04-19 | 新奥数能科技有限公司 | 用于电力系统负荷预测模型的样本数据提取方法和装置 |
CN113899963A (zh) * | 2021-08-30 | 2022-01-07 | 科华数据股份有限公司 | 交流系统晃电检测方法、装置及终端设备 |
CN114062850A (zh) * | 2021-11-17 | 2022-02-18 | 江南大学 | 一种双阈值电网早期故障检测方法 |
CN114062850B (zh) * | 2021-11-17 | 2022-09-23 | 江南大学 | 一种双阈值电网早期故障检测方法 |
CN114784973A (zh) * | 2022-05-18 | 2022-07-22 | 苏州纬讯光电科技有限公司 | 一种单芯高压电缆终端故障预警方法 |
WO2023241327A1 (zh) * | 2022-06-14 | 2023-12-21 | 无锡隆玛科技股份有限公司 | 一种最大特征向量电网异常定位方法 |
WO2023241326A1 (zh) * | 2022-06-14 | 2023-12-21 | 无锡隆玛科技股份有限公司 | 基于样本协方差矩阵最大特征值变化率的电网异常检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108196165A (zh) | 基于样本协方差矩阵最大特征值的电网异常状态检测方法 | |
Tong et al. | Detection and classification of transmission line transient faults based on graph convolutional neural network | |
CN105403811B (zh) | 风电场电网故障诊断方法和装置 | |
Lazzaretti et al. | Novelty detection and multi-class classification in power distribution voltage waveforms | |
CN106443310B (zh) | 一种基于som神经网络的变压器故障检测方法 | |
CN109738776A (zh) | 基于lstm的风机变流器开路故障识别方法 | |
CN103245881B (zh) | 一种基于潮流分布特征的配电网故障分析方法及装置 | |
CN101464964B (zh) | 一种设备故障诊断的支持向量机模式识别方法 | |
Zhang et al. | Fault localization in electrical power systems: A pattern recognition approach | |
CN106249101A (zh) | 一种智能配电网故障辨识方法 | |
Yang et al. | Bad data detection algorithm for PMU based on spectral clustering | |
CN108268905A (zh) | 一种基于支持向量机的变压器故障诊断方法及系统 | |
CN109193650A (zh) | 一种基于高维随机矩阵理论的电网薄弱点评估方法 | |
Hassani et al. | Unsupervised concrete feature selection based on mutual information for diagnosing faults and cyber-attacks in power systems | |
Xu et al. | Real-time event detection and feature extraction using PMU measurement data | |
CN105406461B (zh) | 一种配电网停电事件的负荷自适应动态监测方法 | |
CN110533331A (zh) | 一种基于输电线路数据挖掘的故障预警方法和系统 | |
Kazerooni et al. | Literature review on the applications of data mining in power systems | |
Jia et al. | Defect prediction of relay protection systems based on LSSVM-BNDT | |
Yu et al. | Fault location in distribution system using convolutional neural network based on domain transformation | |
Zhang et al. | A weighted random forest approach to improve predictive performance for power system transient stability assessment | |
Sarlak et al. | SVM‐based method for high‐impedance faults detection in distribution networks | |
Uddin et al. | On the protection of power system: Transmission line fault analysis based on an optimal machine learning approach | |
CN107462810A (zh) | 一种适用于有源配电网的故障区段定位方法 | |
CN110705831A (zh) | 电力系统故障后功角失稳模式预判模型构建方法及其应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180622 |