[发明专利]基于振动分析法的变压器故障分类和识别方法有效
申请号: | 201611003127.4 | 申请日: | 2016-11-15 |
公开(公告)号: | CN106646096B | 公开(公告)日: | 2019-07-02 |
发明(设计)人: | 李敏;陈果;石同春;沈大千;秦少鹏;向天堂;邓权伦;罗宇昆;高翔;陈大浩 | 申请(专利权)人: | 国网四川省电力公司广安供电公司;国家电网公司 |
主分类号: | G01R31/06 | 分类号: | G01R31/06;G01H17/00 |
代理公司: | 成都行之专利代理事务所(普通合伙) 51220 | 代理人: | 田甜 |
地址: | 638500 *** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 变压器故障 样本数据 特征矢量 振动分析 分类 本征模函数分量 变压器运行状态 经验模式分解 主成分分析法 本征模函数 变压器试验 变压器振动 测试样本 二维图像 公式计算 模式识别 输出模式 原始样本 坐标投影 有效地 降维 集合 邻近 直观 采集 | ||
1.基于振动分析法的变压器故障分类和识别方法,其特征在于,包括以下步骤:
S1:选取变压器试验对象,采集不同状态下变压器振动信号为样本数据;
S2:利用希尔伯特-黄变换中集合经验模式分解方法计算样本数据得到本征模函数;
S3:提取本征模函数分量中特征矢量V;
S4:利用主成分分析法对特征矢量进行降维,坐标投影到二维图像中;
S5:利用K邻近法对样本数据进行分类;
S6:利用距离公式计算测试样本与原始样本的距离;
S7:进行模式识别;
S8:输出模式识别中相对应的变压器故障类型;
所述步骤S4中,主成分分析法对特征矢量进行降维步骤如下:
利用EEMD提取变压器振动信号在正常状态、绕组故障、铁芯故障3种工况下特征向量V=[v1,v2,...,vz];
同时针对3种工况,各选择m组作为经验样本,得到4m×z的训练矩阵R;
对训练矩阵R做中心化处理得到矩阵A=[aij],
其中,aij是R中元素中心化后的数据,vij是训练矩阵R中第i行第j列的样本值,vi是训练矩阵R中每行的均值;
利用计算中心化矩阵A的协方差矩阵S,
计算中协方差矩阵S的特征值和特征向量,进行最大值排列,取前两个非零特征值对应的特征向量[α1,α2]为投影方向;
对R进行所示的[α1,α2]方向的投影转换,得到训练样本的二维图;
式中,T表示矩阵转置,Y1、Y2分别为二维投影的横、纵坐标;
所述步骤S7中,模式识别方法如下:
选择K个样本作为测试样本的K个近邻,设w为变压器状态类型,K个组样本中,来自w1状态类型的样本有M1个,来自w2状态类型的样本有M2个,…,来自wc状态类型的样本有Mc个,若k1,k2,…,kc分别是k个近邻中属于w1,w2,…,wc类的样本数,则定义判别函为:
gi(Vt)=ki,i=1,2,…,c;
若gj(Vt)=max(ki),则测试样本Vi∈ωj,测试样本x对应的变压运行状态为wj变压器状态类型。
2.根据权利要求1所述的基于振动分析法的变压器故障分类和识别方法,其特征在于:在步骤S2中,利用希尔伯特-黄变换中集合经验模式分解方法计算样本数据得到本征模函数的步骤为:
向原始信号x(t)中加入M次的高斯白噪音序列ni(t),(i=1,2,…,M)得到Xi(t),即:
Xi(t)=x(t)+ni(t);
对Xi(t)分别进行经验模态分解,得到各个本征模函数的分量和余量,即:
其中hij为第i次加入高斯白噪声后,对Xi(t)进行分解得到的第j个本征模函数分量;rin为第i次加入高斯白噪声后,对Xi(t)进行分解后的余项;n为分解层数;
利用高斯白噪声频谱的零均值原理,将以上步骤对应的本征模函数分量进行总体平均运算,得到的EEMD后的本征模函数的分量为:
式中,hj(t)表示对原始信号进行EEMD分解得到的第j个本征模函数分量;M为加入白噪声的次数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于国网四川省电力公司广安供电公司;国家电网公司,未经国网四川省电力公司广安供电公司;国家电网公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201611003127.4/1.html,转载请声明来源钻瓜专利网。