发明内容
针对上述问题,本发明的目的是提供一种风机传动系统故障特征提取方法,该方法能准确将各种融合在一起的故障特征信息分离、提取、分析,具有较高的故障特征识别精度,解决了故障趋势预测随机不确定性、非线性等关键技术问题。
为实现上述目的,本发明采取以下技术方案:一种风机传动系统故障特征提取方法,其包括以下步骤:(1)利用现有的旋转机械远程在线样本采集及传输实验平台获取现场风机传动系统的实测数据,通过设置在风电机组旋转部件上的多个传感器采集到多个通道的振动信号xi(t),i=1、2、…n,其中,n为传感器个数;(2)根据独立成分分析方法对振动信号xi(t)进行盲源分离,得到原始独立振动信号源sj(t)的近似信号源yj(t),j=1、2、…m;(3)根据近似信号源yj(t)构造信号矩阵Zb×a,并压缩该信号矩阵Zb×a获取新的信号特征矩阵Z′q×a;(4)将新的信号特征矩阵Z′q×a放大并方阵化,构造特征向量T:①取常数矩阵Ba×q为:
其中,e取[2,5]之间的整数;②令Z″=Z′q×aBa×q,获得元素个体被放大的最终特征方阵Z″,令行列式|λhE-Z″|=0,则得到特征方阵Z″的特征值λ1、λ2…λq,h=1,2,…,q,并构造特征向量T=(λ1,λ2…λq);(5)对近似信号源yj(t)进行小波包分解,获得分解系数,并构造系数向量P;(6)根据步骤(4)、步骤(5)中获得的特征向量T=(λ1,λ2…λq)和系数向量P计算各向量元素的方差DT与DP;(7)根据各向量元素的方差DT与DP计算得到故障角θ的正切值tanθ,利用正切值tanθ对风机传动系统故障特征进行评估与预测。
所述步骤(3)中,获取所述新的信号特征矩阵Z′q×a包括以下步骤:①假设近似信号源yj(t)个数为a,向量Y为b维,构造信号矩阵Zb×a为:
②取常数矩阵Aq×b为:
其中,q≤b,且③令Z′q×a=Aq×bZb×a,则信号矩阵Zb×a元素个数被压缩,进而获得新的信号特征矩阵Z′q×a。
所述步骤(2)中,基于所述独立成分分析方法的盲源分离方法如下:①假设有n个振动信号xi(t),且n个振动信号xi(t)的向量形式用X′表示;i=1、2、…n;②对n个振动信号xi(t)的向量X′进行中心化处理,即令X′-E[X′]=X″,使向量X′均值为0;③根据独立成分分析算法中的白化公式对向量X″进行白化处理,进而得到向量X;④根据独立成分分析算法中的分离矩阵W和向量X构造出原始独立振动信号源sj(t)的近似信号源yj(t)的向量Y=WX。
所述步骤(5)中,构造所述系数向量P包括以下步骤:①假设尺度数为l,则原始空间被分解成2l个子空间,将子空间统一标记为Wl 0、Wl 1、Wl 2…Wl 2l-1,其中,Wl 2p表示分解末节点的低频部分,Wl 2p+1表示分解末节点的高频部分,其中p=0、1、2…2l-1-1;令每一个子空间的正交归一基为则有:
其中,h0(k)和h1(k)分别是与尺度函数和小波函数对应的正交镜像滤波器(QMF)的系数,和分别是和空间中的一组正交归一基,且有②在小波包的分解中,在节点(l+1,p)处的小波包系数由下式给出:
式中的③在节点(l,p)处的小波包系数由下式重建:式中,和分别是和每两个点插入一个零后所得到的序列;④由末层节点的小波包分解系数dl(k)构造系数向量P:
所述步骤(7)中,利用所述正切值tanθ对所述风机传动系统故障特征进行评估与预测步骤如下:①计算风机传动系统在正常工作条件下的各向量元素的方差DT0与DP0,则故障角θ的正切值tanθ为:
其中,各向量元素的方差DT与DP是在随机状态下获得的,利用tanθ在值域(-∞,∞)上单调递增这一性质将区间(-∞,∞)再细分成许多个小区间:(-∞,α]、[α,β)…[τ,∞),α、β和τ分别表示各个区间的节点;
②根据已有故障信息再将这些小区间确定为反映不同故障特征故障程度的故障区间,当再次获得新tanθ值时,根据其所在区间判断风机传动系统的故障类型、故障程度及预估故障发展趋势。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于采用了独立成分分析方法(ICA),能有效地从信号中分离出统计独立的非高斯信号源的线性或非线性组合,从而实现分离提纯混叠成采集信号的原始的独立信号源,并判断该源信号是否在向故障的方向发展,达到提前预防故障的目的。2、本发明由于通过最佳小波包分解获得分解系数和对信号矩阵元素个数的压缩与元素本身的放大以及最大信息量抽取等方法,提高了故障信息分离提取效果与运行速度,增强了对原始信号中微弱故障信号的识别能力,解决了故障趋势预示随机不确定性、非线性等问题。3、本发明由于采用将信号特征矩阵压缩成新的信号特征矩阵,并将新的信号特征矩阵放大并方阵化,对新的信号特征矩阵元素个体的放大,目的是使得能够表征被检测部件状态的每一个矩阵元素的特质充分暴露出来,同时获得了信号方阵,最终可以获得更为简洁的表征被检测部件状态的特征向量,提高故障特征识别精度。综上所述,本发明可以广泛在旋转机械故障趋势预测中应用。
具体实施方式
本发明提供一种风机传动系统故障特征提取方法,该方法以故障趋势预测为目的,通过独立成分分析方法(ICA),从信号中分离出统计独立的非高斯信号源的线性或非线性组合,从而实现分离提纯故障相关特征量,为风电机组在线监测预警提供有效依据,保障风电机组安全可靠运行。下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明的风机传动系统故障特征提取方法包括以下步骤:
(1)利用现有的旋转机械远程在线样本采集及传输实验平台获取现场风机传动系统的实测数据,通过设置在风电机组旋转部件(如齿轮箱等)上的多个传感器采集到多个通道的振动信号xi(t),i=1、2、…n,其中,n为通道个数(即传感器个数)。
(2)根据独立成分分析方法(ICA),对振动信号xi(t)进行盲源分离,得到原始独立振动信号源sj(t)的近似信号源yj(t),j=1、2、…m。其中,基于ICA算法的盲源分离方法如下:
①假设有n个振动信号xi(t)(i=1、2、…n),且n个振动信号xi(t)的向量形式用X′表示;
②对n个振动信号xi(t)的向量X′进行中心化处理,即令X′-E[X′]=X″,使其均值为0,其中,E表示数学期望;
③根据ICA算法中的白化公式对向量X″进行白化处理,白化处理如下:
a)首先求解向量X″的协方差矩阵Cx,Cx=E[X″(X″)T];
b)根据协方差矩阵Cx,得到以协方差矩阵Cx的单位范数特征向量为列的矩阵F=(e1…en),其中,ei(i=1,2,…,n)为协方差矩阵Cx的单位范数特征向量;
c)根据协方差矩阵Cx,得到以协方差矩阵Cx的特征值为对角元素的对角矩阵D=diag(d1…dn),其中,di(i=1,2,…,n)为协方差矩阵Cx的特征值;
d)将步骤b)和步骤c)中的矩阵F和对角矩阵D代入ICA算法中的白化公式得到向量X。
④根据ICA算法中的分离矩阵W和向量X构造出原始独立振动信号源sj(t)的近似信号源yj(t)的向量Y=WX;
(3)根据近似信号源yj(t)构造信号矩阵Zb×a,并压缩该信号矩阵Zb×a获取新的信号特征矩阵Z′q×a;其步骤为:
①假设近似信号源yj(t)个数为a,向量Y为b维,构造信号矩阵Zb×a为:
②取常数矩阵Aq×b为:
其中,q≤b,且
③令Z′q×a=Aq×bZb×a,则信号矩阵Zb×a元素个数被压缩,进而获得新的信号特征矩阵Z′q×a。
(4)将新的信号特征矩阵Z′q×a放大并方阵化,构造特征向量T:
①取常数矩阵Ba×q为:
其中,e取[2,5]之间的整数,且根据需要可自由选取e值,但是,一旦选定某一数值,就需在每次针对同一或同类监测目标应用本方法时持续使用该数值,否则最终会产生错误结果。
②令Z″=Z′q×aBa×q,则可获得元素个体被放大的最终特征方阵Z″,令行列式|λhE-Z″|=0,则得到特征方阵Z″的特征值λ1、λ2…λq,h=1,2,…,q,并构造特征向量T为:
T=(λ1,λ2…λq)。
(5)对近似信号源yj(t)进行小波包分解,获得分解系数,并构造系数向量P;
近似信号源yj(t)的小波包分解是将信号投影到小波包基上,获得一系列系数,之后用这一系列系数刻画近似信号源yj(t)的特征,其步骤如下:
①假设尺度数(亦即分解层数)为l,则由完全二叉树关系可知原始空间V0被分解成2l个子空间,将子空间统一标记为Wl 0、Wl 1、Wl 2…Wl 2l-1,其中,Wl 2p表示分解末节点的低频部分,Wl 2p+1表示分解末节点的高频部分,其中p=0、1、2…2l-1-1。令每一个子空间的正交归一基为则有:
其中,h0(k)和h1(k)分别是与尺度函数和小波函数对应的正交镜像滤波器(QMF)的系数,和分别是和空间中的一组正交归一基,且有当p=0,l=0时,函数J0(t-k)即是原始空间V0中的正交归一基O(t-k),即尺度函数;Wl 1中的正交归一基J1(2-lt-k)即是Wl中的正交归一基Jl,k(t),即小波函数。
②在小波包的分解中,在节点(l+1,p)处的小波包系数由下式给出:
式中的
③在节点(l,p)处的小波包系数可由下式重建:
式中,和分别是和每两个点插入一个零后所得到的序列。
④由末层(即第l层)节点的小波包分解系数dl(k)构造系数向量P:
(6)根据步骤(4)、步骤(5)中获得的特征向量T=(λ1,λ2…λq)和系数向量计算各向量元素的方差DT与DP:
其中,与分别表示特征值λ1,λ2…λq与分解系数的平均值,即:
(7)根据各向量元素的方差DT与DP计算得到故障角θ的正切值tanθ,利用故障角θ的正切值tanθ对风机传动系统故障特征进行评估与预测;
①按照步骤(1)至步骤(5),计算风机传动系统在正常工作条件下获得的特征向量T0和系数向量P0,根据步骤(6)计算各向量元素的方差DT0与DP0,则故障角θ的正切值tanθ为:
其中,各向量元素的方差DT与DP是在随机状态下获得的,利用tanθ在值域(-∞,∞)上单调递增这一性质将区间(-∞,∞)再细分成许多个小区间:(-∞,α]、[α,β)…[τ,∞)。其中,α、β和τ分别表示各个区间的节点,可以通过多次实验获得。
②根据已有故障信息再将这些小区间确定为反映不同故障特征故障程度的故障区间,当再次获得新tanθ值时,即可根据其所在区间判断风机传动系统的故障类型、故障程度及预估故障发展趋势。
上述各实施例仅用于说明本发明,各个步骤都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别步骤进行的改进和等同变换,均不应排除在本发明的保护范围之外。