CN112200038B - 一种基于cnn的电力系统振荡类型的快速辨识方法 - Google Patents
一种基于cnn的电力系统振荡类型的快速辨识方法 Download PDFInfo
- Publication number
- CN112200038B CN112200038B CN202011050933.3A CN202011050933A CN112200038B CN 112200038 B CN112200038 B CN 112200038B CN 202011050933 A CN202011050933 A CN 202011050933A CN 112200038 B CN112200038 B CN 112200038B
- Authority
- CN
- China
- Prior art keywords
- oscillation
- power system
- cnn
- layer
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/002—Flicker reduction, e.g. compensation of flicker introduced by non-linear load
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2103/00—Details of circuit arrangements for mains or AC distribution networks
- H02J2103/30—Simulating, planning, modelling, reliability check or computer assisted design [CAD] of electric power networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Nonlinear Science (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于CNN的电力系统振荡类型的快速辨识方法,本发明提供一种分析准确的电力系统振荡类型的快速辨识方法,包括:根据EDSs数学模型生成电力系统振荡样本数据,采用平铺以及零填充算法对振荡样本数据进行预处理操作;根据辨识要求确立分类准则,按照分类准则划分标记电力系统振荡样本数据以供后续网络训练与测试;搭建CNN模型,输入训练样本进行网络训练,通过测试样本分类准确率确定网络训练完成与否;将待测振荡信号经滑窗取样输入CNN,通过输出分析完成对电力系统振荡类型的辨识。本发明具有对电力系统振荡类型快速辨识的有效性和可行性等优点。
Description
技术领域
本发明涉及电力系统稳定与控制技术领域,具体涉及一种基于卷积神经网络(convolutional neural network,CNN)的电力系统振荡类型的快速辨识方法。
背景技术
自电力系统诞生以来,振荡研究就成为了关乎电力系统动态性能与稳定性的重要研究之一。经过多年的研究,传统电力系统中低频振荡(low frequency oscillation,LFO)和次同步振荡(sub-synchronous oscillation,SSO)的特性已被充分揭示。它们的共同特性是由较大物理惯量的旋转单元所主导,尤其是大型同步发电机组。然而,由于可再生能源的快速发展,如太阳能、风能、地热能、生物质能等,大量的新能源并网发电装备通过电力电子逆变装置接入电网,这给电力系统的规划以及运行控制带来了巨大的挑战。传统控制策略的电力电子逆变器侧重于发电和电能质量的控制,具有响应速度快、低惯量等特性,大量的采用传统电力电子技术,会导致系统呈现出惯量减弱的趋势,严重危及电力系统的稳定运行。例如,2009年10月,美国德克萨斯州的双馈风力发电系统和串联补偿电网之间发生了频率约为20Hz的次同步谐振事故,造成消弧电路损坏和大规模断路故障。考虑到电力系统振荡对电力系统稳定运行的严重威胁,快速辨识出电力系统振荡类型对电网后续动态稳定性的评估及调整具有重要意义。
目前,电力系统振荡类型的辨识主要是模态辨识,对应的方法分为基于模型的方法和基于实测信号的方法。由于无需系统精确的模型和参数,后者在电力系统中得到了广泛应用。常见算法有基于快速傅里叶变换(fast fourier transformation,FFT)的算法、小波分析、Prony、希尔伯特-黄变换(HHT)等以及近年来发展起来的借助旋转不变技术估计信号参数(ESPRIT:estimation ofsignal parameters via rotational invariancetechnique),矩阵束方法(MP:matrix pencil algorithm),随机子空间方法(SSI:stochastic subspace identification)等。经过多年的研究与改进,这些方法可以较为准确地获取系统的模态,在抗噪性上也有一定的效果,但在辨识的数学机理上,这些算法将被测信号视作平稳随机过程,并且往往需要采集数秒钟乃至数分钟的数据。随着可再生能源并网比例的提高以及大量的电力电子设备接入电网,所产生的振荡表现出不同于传统电力系统振荡的新特性。由于可再生能源的发电机单元和AC/DC电网之间的复杂动态交互,所产生的振荡一般具有从几赫兹到几千赫兹的宽频范围,同时还具有强时变性、非线性、非平稳性以及强噪声干扰等特性。在这种背景下,传统电力系统中振荡信号近似平稳的假设难以保证,亟需发展新的振荡类型辨识方法来适应电力系统的发展。
近年来,深度学习(deep learning,DL)发展迅猛,在众多领域已经得到了成功的应用。在电力系统振荡类型辨识方面,也有学者开展了探索性的研究。有研究提出使用深度学习算法辨识电力系统振荡模式是区域振荡还是局部振荡。还有研究提出指数型衰减正弦神经网络(exponentially damped sinusoids neural network,EDSNN)的辨识方法。但这些研究,均存在对被测振荡信号非平稳性特性考虑不足的缺陷。目前,同步相量测量单元(PMU)的广域测量系统(WAMS)在电力系统中获得了广泛应用,为电力系统振荡分析提供了数据来源,让深度学习在电力系统振荡类型辨识方面的应用成为可能。同时,由于事先通过大量的数据学习,深度学习的引入还使得电力系统振荡类型的辨识具有神经网络类算法的共同优点——辨识速度迅速。而这往往是系统状态强时变的高比例可再生能源并网电力系统中迫切需要的一个特点。因此,随着可再生能源渗透率的进一步提高,电力系统中的振荡模式愈发频繁且复杂,如何借鉴深度学习算法实现电力系统振荡类型的在线快速辨识,成为了一个亟待解决的技术问题。
发明内容
本发明的目的在于针对高比例可再生能源并网的电力电子化系统中易出现的复杂振荡现象,提供一种分析准确,设计合理的基于CNN的电力系统振荡类型的快速辨识方法。
本发明可通过下述技术方案实现:
一种基于CNN的电力系统振荡类型的快速辨识方法,由以下步骤构成:
步骤1:生成振荡样本数据,即根据电力系统振荡信号数学模型为CNN模型生成训练样本数据与测试样本数据。
所选振荡信号数学模型为指数型衰减正弦量(exponentially dampedsinusoids,EDSs),其公式为:
式(1)中,x(tj)为tj时刻的信号;Ai为幅值;σi为衰减因子;fi为频率;为相位;j=0,1,2,…,n-1;n为采样点数;i=1,2,…,m;m为模型实际阶数;为第i阶信号模态加入的时间节点;η(tj)代表噪声信号。
步骤2:采用预处理算法处理电力系统振荡样本数据以便后续训练过程中的特征提取;
预处理算法采用平铺的方式将序列数据处理为n阶方阵,其公式为:
式(2)中,X为离散的时序序列,长度为N,由于目标矩阵为n阶方阵,因此N=n2;
此外,为了实现对特征矩阵的边缘滤波,对上述n阶方阵实施零填充操作,即在输入矩阵的边缘使用零值进行填充,使原输入矩阵转换为m(m>n)阶方阵。
步骤3:确定电力系统振荡类型的分类准则,根据分类准则为样本数据添加类别标签。
分类准则为按存在性分类。一般而言,电力系统振荡信号按不同频率可大致分为四类,分别为:低频振荡(f≤2.5Hz)、次同步振荡(2.5Hz<f≤45Hz)、超同步振荡(50Hz<f<100Hz)、中高频振荡(f≥100Hz);对于某一高阶振荡信号,其可能包含任一类型振荡模态或同时包含多种类型振荡模态,按存在性分类即分别判断该四种振荡模态于某一振荡信号中是否存在,再按存在与不存在分为两类进行辨识,因此本发明需按照四种振荡类型训练四个不同的分类网络。
步骤4:根据振荡样本数据特征以及分类需求构建CNN模型。
本发明所搭建的CNN模型,具体结构包含输入层、卷积层、池化层、全连接层、Dropout层以及输出层;
其中,输入层用以输入序列数据;卷积层为网络核心部分,负责将原始数据映射到隐层的特征空间;池化层用以缩小数据尺寸,减少全连接层的参数并加快训练速度;全连接层用以完成分布式特征表示到样本标记空间的映射;Dropout层通过令数据随机失活来预防过拟合现象的发生;输出层则输出最终的分类结果。
步骤5:将训练样本输入CNN模型,训练模型的分类能力。
CNN模型训练过程为一个循环迭代过程,主要通过前向传播算法和反向传播算法更新网络各层级状态与参数,具体包括以下步骤:
S51:初始化各网络参数,包括各神经元的权重和偏置等;
S52:输入训练样本数据;
S53:计算CNN网络内部各层级输出;
S54:计算当前序列索引预测输出及其损失函数;
S55:计算所有参数基于损失函数的偏导数,通过梯度下降法更新网络所有的参数;
S56:当损失函数不再下降或达到设定训练次数,完成训练。
进一步的,在CNN模型训练过程中,所述CNN网络各层级输出的计算公式为:
1)卷积层输出:
式(3)中,i=1,2,…,I;j=1,2,…,J;I=P-M+1;J=Q-N+1;P、Q为输入矩阵X的尺寸,km,n是M×N的卷积核的第m行,第n列的元素,xi+m-1,j+n-1是输入矩阵X中的第i+m-1行,第j+n-1列的元素。每个卷积核都会产生一个I×J的输出矩阵Y;
2)池化层输出(平均池化):
式(4)中,SR为池化区域面积,即池化区域元素个数;Yi,j为卷积层输出Y第i行第j列的元素;
3)全连接层输出:
式(5)中,hn为输出神经元的值;vm为输入神经元的值;wmn、bn分别为链接权值与偏置,m为输入神经元的数目;n为输出神经元数目。
进一步的,在CNN模型训练过程中,所述预测输出值的计算公式及其损失函数的选取分别为:
式中,V和c表示预测输出的权值和偏置;tanh表示双曲正切函数;N为每次送入网络训练的批量大小;K为类型数目;yn,ij表示第i批次及第j类型样本的实际标签,则表示第i批次及第j类型样本的预测输出值。
步骤6:将测试样本输入CNN模型,测试各网络模型的辨识准确率,并通过调节网络参数不断提升辨识准确率以完成网络训练,当辨识准确率达到90%以上时,停止网络参数的调节,完成网络训练获得训练完成的CNN模型。
步骤7:通过滑动时窗获取电力系统振荡信号实测数据。对某段振荡信号采用滑动时窗取样,其过程如式(9)所示:
式(9)中,X为离散的时序序列,长度为N;L为滑窗长度;fs为信号采样频率,由于滑动间隔为0.5s,则滑动步长等于fs/2;L1、L2、…、Ln+1即为通过滑动时窗获取的不同时段振荡信号实测数据。
依此法获取的待测信号按时序排列,通过对这些待测信号顺序逐一辨识,可以有效分析出待测振荡样本数据各模态的时变特征,以适应电力系统中振荡频繁且复杂的环境。
步骤8:处理待测振荡信号并输入CNN模型,根据网络输出分析振荡信号的所属类型。根据存在性分类准则,所有网络输出结果均为包含对应模态分段或不包含对应模态分段这两种情况。对于某段待测振荡信号,将其输入各网络模型并综合分析输出结果,即可获得其所包含的振荡类型。
本发明所述振荡类型辨识方法采用深度学习算法,选取CNN作为振荡样本训练模型,能够从较短时的振荡信号中迅速分析出其所属的类型,实现对电力系统振荡类型的在线快速辨识。
本发明与现有技术相比,具有如下的优点和有益效果:
(1)本发明引入了人工智能领域的深度学习算法,利用深度神经网络的分类能力实现电力系统振荡类型的辨识,相较于传统的信号模型分析法,辨识更快更准确,更加适用于电力系统中振荡类型的实时辨识。
(2)本发明充分考虑到真实电网中常见的不稳定与时变振荡信号,采用深度学习算法并且选择定性辨识解决这一问题,使其更加适用于高比例可再生能源并网电力电子化系统中相对复杂的运行工况与振荡环境,相较于传统的信号模型分析法会将其视作稳定与时不变信号处理,得到的辨识结果既不是模态的真实数值也无法体现时变特征更有优势。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明的工作流程图。
图2为本发明的CNN模型结构图。
图3为本发明的CNN训练流程图。
图4为含噪声情况下理想振荡信号图像。
图5为电力系统实测振荡信号图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
本发明中,S51-S56表示步骤5的详细子步骤。
如图1所示,一种基于CNN的电力系统振荡类型的快速辨识方法,包括如下步骤:步骤1:生成样本数据,即根据电力系统振荡信号数学模型为CNN模型生成训练样本数据与测试样本数据。
所选振荡信号数学模型为EDSs,其公式为:
式(1)中,x(tj)为tj时刻的信号;Ai为幅值;σi为衰减因子;fi为频率;为相位;i=1,2,…,m;j=0,1,2,…,n-1;n为采样点数;m为模型实际阶数;为第i阶信号模态加入的时间节点;η(tj)代表噪声信号。
本实施例中,振荡信号采样时间取1s,采样频率取400Hz,其中采样时间即振荡信号模态辨识所需最短时间。生成信号过程中,振荡信号各模态参数在合理范围内随机取值,最高阶数取到4。此外,为验证本发明抗噪性能,噪声信号η(t)信噪比(signal-noiseratio,SNR)取10dB。
步骤2:采用预处理算法处理电力系统振荡样本数据以便后续训练过程中的特征提取。
预处理算法采用平铺的方式将序列数据处理为n阶方阵,其公式为:
式(2)中,X为离散的时序序列,长度为N=400,由于目标矩阵为n阶方阵,因此n=20;
此外,为了实现对特征矩阵的边缘滤波,对上述20×20矩阵实施零填充操作,即在输入矩阵的边缘使用零值进行填充,使原输入矩阵转换为28阶方阵。
步骤3:确定电力系统振荡类型的分类准则,根据分类准则为样本数据添加类别标签。
分类准则为按存在性分类。一般而言,电力系统振荡信号按不同频率可大致分为四类,分别为:低频振荡(f≤2.5Hz)、次同步振荡(2.5Hz<f≤45Hz)、超同步振荡(50Hz<f<100Hz)、中高频振荡(f≥100Hz);对于某一高阶振荡信号,其可能包含任一类型振荡模态或同时包含多种类型振荡模态,按存在性分类即分别判断该四种振荡模态于某一振荡信号中是否存在,再按存在与不存在分为两类进行辨识,因此本发明需按照四种振荡类型训练四个不同的分类网络。
该分类方法可以实现不同阶数振荡信号的同时辨识,无需依赖传统辨识方法中难以实现的精确定阶过程,因此辨识结果更加可信。
步骤4:根据振荡样本数据特征以及分类需求构建CNN模型。
如图2所示,本发明所搭建的CNN模型具体结构为:输入层-卷积层1-池化层1-卷积层2-池化层2-全局平均池化层-Dropout层1-全连接层1-Dropout层2-全连接层2-Dropout层3-全连接层3-Softmax层-输出层。
其中,输入层用以输入序列数据,根据输入数据格式,将输入层型号设定为28×28×1;卷积层为网络核心部分,负责将原始数据映射到隐层的特征空间,其中卷积层1与卷积层2的卷积核大小均设定为4×4;池化层用以缩小数据尺寸,减少全连接层的参数并加快训练速度,其中池化层1与池化层2均采用最大池化,其大小均设定为2×2;全局平均池化层用以替代全连接层,可以有效减少参数数量并防止过拟合现象发生;全连接层与Softmax层用以完成分布式特征表示到样本标记空间的映射以及结果归一化;Dropout层置于全连接层之后,同样为了防止过拟合现象发生,其随机失活概率为50%;输出层则输出最终的分类结果。
步骤5:将训练样本输入CNN模型,训练模型的分类能力。
CNN模型训练过程为一个循环迭代过程,如图3所示,主要通过前向传播算法和反向传播算法更新网络各层级状态与参数,具体包括以下步骤:
S51:初始化各网络参数,包括各神经元的权重和偏置等;
S52:输入训练样本数据;
S53:计算CNN网络内部各层级输出;
其中,CNN网络内部各层级输出的计算公式为:
1)卷积层输出:
式(3)中,i=1,2,…,I;j=1,2,…,J;I=P-M+1;J=Q-N+1;P、Q为输入矩阵X的尺寸,km,n是M×N的卷积核的第m行,第n列的元素,xi+m-1,j+n-1是输入矩阵X中的第i+m-1行,第j+n-1列的元素。每个卷积核都会产生一个I×J的输出矩阵Y;
2)池化层输出(平均池化):
式(4)中,SR为池化区域面积,即池化区域元素个数;Yi,j为卷积层输出Y第i行第j列的元素;
3)全连接层输出:
式(5)中,hn为输出神经元的值;vm为输入神经元的值;wmn、bn分别为链接权值与偏置,m为输入神经元的数目;n为输出神经元数目。
S54:计算当前序列索引预测输出及其损失函数;
其中,预测输出值的计算公式及其损失函数的选取分别为:
式中,V和c表示预测输出的权值和偏置;tanh表示双曲正切函数;N为每次送入网络训练的批量大小;K为类型数目;yn,ij表示第i批次及第j类型样本的实际标签,则表示第i批次及第j类型样本的预测输出值。
S55:计算所有参数基于损失函数的偏导数,通过梯度下降法更新网络所有的参数;
S56:当损失函数不再下降或达到设定训练次数,完成训练。
步骤6将测试样本输入CNN模型,测试各网络模型的辨识准确率,并通过调节网络参数使得辨识准确率能够达到95%左右,得到训练完成的神经网络模型。
根据上述步骤获得训练完成的神经网络模型,步骤7与步骤8将由两个具体的实施例进行阐释。
实施例1:
为验证该算法是否能够辨识出系统在振荡过程中叠加新的振荡模态,构造如下理想振荡测试信号:
式(10)中,ε(t)表示阶跃函数,η(t)表示噪声信号。
该振荡信号图像如图4所示,信号长度2s,信噪比SNR=10dB。1s前,信号包含两个模态,其中频率f1=0.84Hz,f2=1.21Hz,属于低频振荡。t=1s时,引入一个新的振荡模态,其频率f3=8.8Hz,属于次同步振荡,因此1s后的振荡信号中同时包含低频振荡与次同步振荡两种振荡类型。
步骤7通过滑动时窗获取待测振荡信号样本。滑窗长度为1s,滑动间隔为0.5s,采样频率为400Hz。为了辨识振荡信号中发生的模态改变,需在1s前后分别使用滑窗取样。本实施例中,选取四段信号作为辨识对象,分别为0~1s、0.5~1.5s与1~2s。其中0.5~1.5s包含了模态发生变化的时间节点,用以验证本发明针对复杂振荡模态辨识的能力。
步骤8处理待测振荡信号并输入CNN模型,根据网络输出分析振荡类型的辨识结果,如下表1所示:
表1.基于实施例1的振荡类型的辨识结果表
参考上表1,可以发现该方法的辨识结果与理想振荡信号的真实类型完全吻合。当所辨识的振荡信号涉及新振荡类型的加入,该方法依旧适用。此外,该方法所需辨识时间仅1s,远快于传统振荡信号辨识方法,因此从各方面考虑,该方法都更具优势。
实施例2:
为了验证本发明的实际辨识效果,从电力系统中获取一段振荡信号实测数据。如图5所示,该段振荡信号由一处小扰动所激发,时间位于4s处。为了辨识该扰动所激发的振荡类型,因此截取扰动结束后的数据作为本实施例中的待测振荡信号。
步骤7:通过滑动时窗获取待测振荡信号样本。滑窗长度为1s,滑动间隔为0.5s,采样频率为400Hz。本实施例中,于扰动结束后滑取两段信号作为辨识对象,分别为4~5s与6~7s,以验证振荡类型是否会发生改变。
步骤8:处理待测振荡信号并输入CNN模型,根据网络输出分析振荡类型的辨识结果,如下表2所示:
表2.基于实施例2的振荡类型的辨识结果表
参考上表2,可以判断该系统中由扰动所激发的振荡类型为低频振荡。由于采用滑窗采样,该方法可以伴随振荡信号的产生不断实时返回所产生的振荡类型,判断是否存在新型振荡的加入。此外,1s的识别时间使得该方法能迅速获取辨识结果,因此相较于传统振荡信号辨识方法,本发明更具实用价值。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,方法包括以下步骤:
步骤1:根据电力系统振荡信号数学模型为CNN模型生成振荡样本数据,所述振荡样本数据包括训练样本数据和测试样本数据;
步骤2:采用预处理算法处理振荡样本数据以便后续训练过程中的特征提取;
步骤3:确定电力系统振荡类型的分类准则,根据分类准则为振荡样本数据添加类别标签;
步骤4:根据振荡样本数据特征以及分类需求构建CNN模型;
步骤5:将训练样本输入CNN模型,训练模型的分类能力;
步骤6:将测试样本输入CNN模型,测试各CNN模型的辨识准确率,并通过调节网络参数不断提升辨识准确率以完成网络训练;
步骤7:通过滑动时窗获取待测振荡信号实测数据;
步骤8:处理待测振荡信号并输入CNN模型,根据CNN模型输出分析振荡信号的所属类型;其中:
步骤2中所述预处理算法采用平铺的方式将序列数据处理为n阶方阵,其公式为:
式(2)中,X为离散的时序序列,长度为N,由于目标矩阵为n阶方阵,因此N=n2;所述n阶方阵采用零填充操作,将原输入矩阵转换为m阶方阵,m>n;
步骤3中所述分类准则为按存在性分类,具体方法为:
电力系统振荡信号按不同频率可分为四类,分别为:低频振荡、次同步振荡、超同步振荡、中高频振荡;对于某一高阶振荡信号,其可能包含任一类型振荡模态或同时包含多种类型振荡模态,按存在性分类即分别判断该四种振荡模态于某一振荡信号中是否存在,再按存在与不存在分为两类进行辨识,因此需按照四种振荡类型训练四个不同的分类网络;所述低频振荡的频率f满足f≤2.5Hz,所述次同步振荡的频率f满足2.5Hz<f≤45Hz,所述超同步振荡的f满足50Hz<f<100Hz,所述中高频振荡f满足f≥100Hz。
2.根据权利要求1所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,步骤1中所述振荡信号数学模型为指数型衰减正弦量,其公式为:
式(1)中,x(tj)为tj时刻的信号;Ai为幅值;σi为衰减因子;fi为频率;为相位;i=1,2,…,m;j=0,1,2,...,n-1;n为采样点数;m为模型实际阶数;为第i阶信号模态加入的时间节点;η(tj)代表噪声信号。
3.根据权利要求1所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,步骤4中所述CNN模型,具体结构包含输入层、卷积层、池化层、全连接层、Dropout层以及输出层;
所述输入层用以输入序列数据;卷积层用以将原始数据映射到隐层的特征空间;池化层用以缩小数据尺寸,减少全连接层的参数并加快训练速度;全连接层用以完成分布式特征表示到样本标记空间的映射;Dropout层通过令数据随机失活来预防过拟合现象的发生;输出层则输出最终的分类结果。
4.根据权利要求1所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,步骤5中所述网络训练过程为一个循环迭代过程,通过前向传播算法和反向传播算法更新网络各层级状态与参数,具体包括以下步骤:
S51:初始化各网络参数,包括各神经元的权重和偏置;
S52:输入训练样本数据;
S53:计算CNN网络内部各层级输出;
S54:计算当前序列索引预测输出值及其损失函数;
S55:计算所有参数基于损失函数的偏导数,通过梯度下降法更新网络所有的参数;
S56:当损失函数不再下降或达到设定训练次数,完成训练。
5.根据权利要求4所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,步骤S53中所述CNN网络各层级输出的计算公式为:
1)卷积层输出:
式(3)中,i=1,2,…,I;j=1,2,…,J;I=P-M+1;J=Q-N+1;P、Q为输入矩阵X的尺寸,km,n是M×N的卷积核的第m行,第n列的元素,xi+m-1,j+n-1是输入矩阵X中的第i+m-1行,第j+n-1列的元素;每个卷积核都会产生一个I×J的输出矩阵Y;
2)池化层输出(平均池化):
式(4)中,SR为池化区域面积,即池化区域元素个数;Yi,j为卷积层输出Y第i行第j列的元素;
3)全连接层输出:
式(5)中,hn为输出神经元的值;vm为输入神经元的值;wmn、bn分别为链接权值与偏置,m为输入神经元的数目;n为输出神经元数目。
6.根据权利要求5所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,步骤S54中所述预测输出的计算公式以及损失函数的选取分别为:
式(6)、(7)和(8)中,V和c表示预测输出的权值和偏置;tanh表示双曲正切函数;N为每次送入网络训练的批量大小;K为类型数目;yn,ij表示第i批次及第j类型样本的实际标签,则表示第i批次及第j类型样本的预测输出值。
7.根据权利要求1所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,所述步骤6中,当辨识准确率达到90%以上时,停止网络参数的调节,完成网络训练获得训练完成的CNN模型。
8.根据权利要求1所述的一种基于CNN的电力系统振荡类型的快速辨识方法,其特征在于,步骤7中所述通过滑动时窗获取待测振荡信号实测数据,如式(9)所示:
式(9)中,X为离散的时序序列,长度为N;L为滑窗长度;fs为信号采样频率,由于滑动间隔为0.5s,则滑动步长等于fs/2;L1、L2、…、Ln+1即为通过滑动时窗获取的不同时段振荡信号实测数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011050933.3A CN112200038B (zh) | 2020-09-29 | 2020-09-29 | 一种基于cnn的电力系统振荡类型的快速辨识方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011050933.3A CN112200038B (zh) | 2020-09-29 | 2020-09-29 | 一种基于cnn的电力系统振荡类型的快速辨识方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112200038A CN112200038A (zh) | 2021-01-08 |
| CN112200038B true CN112200038B (zh) | 2023-12-05 |
Family
ID=74008507
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011050933.3A Active CN112200038B (zh) | 2020-09-29 | 2020-09-29 | 一种基于cnn的电力系统振荡类型的快速辨识方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112200038B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113486965A (zh) * | 2021-07-14 | 2021-10-08 | 西南交通大学 | 一种车网电气耦合数据的异常辨识模型的训练方法 |
| CN115133522A (zh) * | 2022-05-24 | 2022-09-30 | 沈阳工业大学 | 风机并网次同步振荡监测网络构建与分层重构抑制方法 |
| CN115187154B (zh) * | 2022-09-14 | 2022-12-16 | 华中科技大学 | 一种基于神经网络的区域电网振荡源风险预测方法及系统 |
| CN115377999B (zh) * | 2022-09-26 | 2024-07-26 | 华北电力大学 | 一种基于宽频测量数据的次同步振荡辨识的方法 |
| CN115659220B (zh) * | 2022-10-25 | 2025-09-30 | 国网河北省电力有限公司电力科学研究院 | Ssci的振荡类型的识别方法、模型训练方法及设备 |
| CN117350170B (zh) * | 2023-11-20 | 2024-02-09 | 华北电力大学(保定) | 基于koopman深度神经网络的非线性振荡分析方法 |
| CN117892130B (zh) * | 2024-01-16 | 2025-03-14 | 云南大学 | 基于轻量卷积神经网络的工业过程振荡检测方法及系统 |
| CN118194239B (zh) * | 2024-05-15 | 2024-09-03 | 合肥工业大学 | 一种基于多模态图网络的输电系统宽频振荡参数辨识方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1981453A (zh) * | 2004-06-17 | 2007-06-13 | W5网络公司 | 低功率无线通信系统和协议 |
| CN104578115A (zh) * | 2015-01-26 | 2015-04-29 | 国网四川省电力公司经济技术研究院 | 一种基于相关函数的电力系统低频振荡模式辨识方法 |
| CN105473065A (zh) * | 2014-02-11 | 2016-04-06 | 皇家飞利浦有限公司 | 确定在cpr期间的自主循环恢复 |
| CN109255394A (zh) * | 2018-10-18 | 2019-01-22 | 国网天津市电力公司电力科学研究院 | 一种基于模式相似度的强迫振荡识别方法 |
| CN110236536A (zh) * | 2019-06-04 | 2019-09-17 | 电子科技大学 | 一种基于卷积神经网络的脑电高频振荡信号检测系统 |
| CN110261746A (zh) * | 2019-07-08 | 2019-09-20 | 清华大学深圳研究生院 | 基于振荡波电压周期衰减特性的电缆缺陷检测方法 |
| CA3105412A1 (en) * | 2018-07-06 | 2020-01-09 | Wobben Properties Gmbh | Wind energy system and method for identifying low-frequency oscillations in an electrical supply network |
| CN111046327A (zh) * | 2019-12-18 | 2020-04-21 | 河海大学 | 适用于低频振荡与次同步振荡辨识的Prony分析方法 |
| CN111160167A (zh) * | 2019-12-18 | 2020-05-15 | 北京信息科技大学 | 基于s变换深度卷积神经网络的主轴故障分类识别方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7108192B2 (en) * | 1999-09-17 | 2006-09-19 | Silverbrook Research Pty Ltd | Rotationally symmetric tags |
-
2020
- 2020-09-29 CN CN202011050933.3A patent/CN112200038B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1981453A (zh) * | 2004-06-17 | 2007-06-13 | W5网络公司 | 低功率无线通信系统和协议 |
| CN105473065A (zh) * | 2014-02-11 | 2016-04-06 | 皇家飞利浦有限公司 | 确定在cpr期间的自主循环恢复 |
| CN104578115A (zh) * | 2015-01-26 | 2015-04-29 | 国网四川省电力公司经济技术研究院 | 一种基于相关函数的电力系统低频振荡模式辨识方法 |
| CA3105412A1 (en) * | 2018-07-06 | 2020-01-09 | Wobben Properties Gmbh | Wind energy system and method for identifying low-frequency oscillations in an electrical supply network |
| CN109255394A (zh) * | 2018-10-18 | 2019-01-22 | 国网天津市电力公司电力科学研究院 | 一种基于模式相似度的强迫振荡识别方法 |
| CN110236536A (zh) * | 2019-06-04 | 2019-09-17 | 电子科技大学 | 一种基于卷积神经网络的脑电高频振荡信号检测系统 |
| CN110261746A (zh) * | 2019-07-08 | 2019-09-20 | 清华大学深圳研究生院 | 基于振荡波电压周期衰减特性的电缆缺陷检测方法 |
| CN111046327A (zh) * | 2019-12-18 | 2020-04-21 | 河海大学 | 适用于低频振荡与次同步振荡辨识的Prony分析方法 |
| CN111160167A (zh) * | 2019-12-18 | 2020-05-15 | 北京信息科技大学 | 基于s变换深度卷积神经网络的主轴故障分类识别方法 |
Non-Patent Citations (7)
| Title |
|---|
| A Convolution Neural Network Method for Power System Oscillation Type Identification;Qianyu Li等;《2020 IEEE 4th Conference on Energy Internet and Energy System Integration (EI2)》;499-504 * |
| A Novel Neural Network Approach for Power System Low Frequency Oscillation Mode Identification;Zhongting Shen等;《2019 IEEE International Symposium on Circuits and Systems (ISCAS)》;1-5 * |
| A Robust Direct Parameter Identification of Exponentially Damped Low-Frequency Oscillation in Power Systems;Zhaobi Chu等;《Journal of Sensors》;1-11 * |
| 基于小波变换和神经网络的暂态电能质量扰动自动识别;刘晓芳等;《继电器》(第23期);46-50 * |
| 基于改进贝叶斯分类法的电能质量扰动分类方法;张文涛等;《电网技术》(第07期);22-25 * |
| 基于神经网络的电力负荷预测方法研究;罗宁等;《自动化与仪器仪表》(第01期);157-160 * |
| 多直驱风机经VSC-HVDC并网系统场内/场网次同步振荡特性分析;邵冰冰等;《中国电机工程学报》;第40卷(第12期);3835-3847 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112200038A (zh) | 2021-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112200038B (zh) | 一种基于cnn的电力系统振荡类型的快速辨识方法 | |
| CN110059357B (zh) | 一种基于自编码网络的智能电能表故障分类检测方法及系统 | |
| CN112183368B (zh) | 基于lstm的电力系统低频振荡模态特征的快速辨识方法 | |
| CN112751345B (zh) | 基于lstm与相轨迹的电力系统低频振荡模态辨识方法 | |
| CN112200694B (zh) | 基于图神经网络的主导失稳模式识别模型构建及应用方法 | |
| CN114021433B (zh) | 一种电力系统的主导失稳模式识别模型的构建方法及应用 | |
| CN108333468B (zh) | 一种有源配电网下不良数据的识别方法及装置 | |
| CN113659565B (zh) | 一种新能源电力系统频率态势的在线预测方法 | |
| CN112434743B (zh) | 基于gil金属微粒局部放电时域波形图像的故障识别方法 | |
| CN113612237A (zh) | 一种在海上风电场中定位谐振诱发的次同步振荡源的方法 | |
| CN105678343A (zh) | 基于自适应加权组稀疏表达的水电机组噪声异常诊断方法 | |
| CN111398798B (zh) | 基于振动信号区间特征提取的断路器储能状态辨识方法 | |
| CN118275818B (zh) | 一种配网输电线路故障检测方法及系统 | |
| Rezaeian Koochi et al. | New phasor‐based approach for online and fast prediction of generators grouping using decision tree | |
| CN115097354B (zh) | 配电网故障线路确定方法、系统、电子设备及介质 | |
| CN110378286A (zh) | 一种基于dbn-elm的电能质量扰动信号分类方法 | |
| CN114779002B (zh) | 一种输电线路故障点定位方法、装置、设备及存储介质 | |
| CN116662899A (zh) | 一种基于自适应策略的含噪数据异常检测方法 | |
| Shahriyari et al. | A deep Learning-Based approach for comprehensive rotor angle Stability Assessment | |
| CN114066214A (zh) | 一种基于多融合卷积神经网络的电能质量分析方法 | |
| CN116821771B (zh) | 基于istoa-cnn的hvdc系统故障诊断方法及系统 | |
| CN112329535B (zh) | 基于cnn的电力系统低频振荡模态特征的快速辨识方法 | |
| CN117077546A (zh) | 基于数据驱动的电力系统负荷建模方法及其系统 | |
| CN116131313A (zh) | 特征量与暂态功角稳定性关联关系的可解释性分析方法 | |
| CN113657149B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |