基于二维投射的特征指标对煤矿设备状态识别的敏感性和聚类性研究方法

摘要

本发明公开了一种基于二维投射的特征指标对煤矿设备状态识别的聚类性和敏感性研究方法，涉及煤矿主扇风机和提升机。本发明的研究方法是采集不同设备状态下的振动信号并加工处理得到24种特征指标，所述24种特征指标包含13个时域特征指标和11个频域特征指标，然后利用二维投射方法将所述24个特征指标投射到二维空间，提出最佳投射方向矩阵有效的避免了单一投射方向下在寻优过程中对新空间的探索能力有限、容易收敛到局部最优解、结果稳定性可靠性差的缺点，研究在最佳投影方向矩阵下所述24种特征指标的投影分布特点；提出用极差系数、平均差系数、离散系数、主轴线相对系数和均值系数来研究24种特征指标的投影值，从而实现所述24特征指标对不同故障的敏感性和聚类性的衡量和评价。

主权项

一种基于二维投射的特征指标对煤矿设备状态识别的聚类性和敏感性研究方法，其特征在于，该方法包括以下步骤：(1)在煤矿设备关键部位——轴承的水平、垂直、轴向三个方向上安装加速度振动传感器，不间断地监测设备运行过程中的振动数据，然后将所述振动数据发送至后端接收和处理系统，并存入相关数据文件；(2)对加速度振动传感器中提取的振动数据{a_i,i＝1,2,...,n}进行处理，具体步骤为：①对所述振动数据{a_i,i＝1,2,...,n}进行傅里叶变换得到{x_i,i＝1,2,...,n}；②按下列公式计算出设备某运行状态的24种特征指标的样本：{x_avi}，{x_pi}，{x_rmsi}，{x_ri}，{D_xi}，{x_p‑pi}，{α_i}，{β_i}，{C_fi}，{S_fi}，{I_fi}，{CL_fi}，{K_vi}，{f_avgi}，{f_bi}，{f_bbi}，{f_vi}，{f_rvi}，{S_i}，{E_r1i}，{E_r2i}，{E_r3i}，{E_r4i}，{E_r5i}，绝对均值：$x_{av}=\frac{1}{n}\underset{i=1}{\overset{n}{\Sigma }}\left|x_i\right|;$峰值：x_p＝max|x_i|；有效值：$x_{rms}=\sqrt{\frac{1}{n}\underset{i=1}{\overset{n}{\Sigma }}{x}_i^2};$方根幅值：$x_r={\left(\frac{1}{n}\underset{i=1}{\overset{n}{\Sigma }}\sqrt{\left|x_i\right|}\right)}^2;$方差：$D_x=\frac{1}{n}\underset{i=1}{\overset{n}{\Sigma }}{(x_i-x_{av})}^2;$峰‑峰值：x_p‑p＝max(x_i)‑min(x_i)；偏态指标：$\alpha =\sqrt{\frac{1}{6n}}\underset{i=1}{\overset{n}{\Sigma }}{\left(\frac{x_i-x_{av}}{\sqrt{D_x}}\right)}^3;$峭度指标：$\beta =\sqrt{\frac{n}{24}}\left[\underset{i=1}{\overset{n}{\Sigma }}{\left(\frac{x_i-x_{av}}{\sqrt{D_x}}\right)}^4-3\right];$峰值指标：$C_f=\frac{x_p}{x_{rms}};$波形指标：$S_f=\frac{x_{rms}}{x_{av}};$脉冲指标：$I_f=\frac{x_p}{x_{rms}};$裕度指标：${\mathrm{CL}}_f=\frac{x_p}{x_r};$变异系数：$K_v=\frac{\sqrt{D_x}}{x_{av}};$重心频率：$f_{avg}=\frac{{\int }_0^{\infty }fp\left(f\right)df}{{\int }_0^{\infty }p\left(f\right)df};$均方频率：$f_b=\frac{{\int }_0^{\infty }{f}^{2}p\left(f\right)df}{{\int }_0^{\infty }p\left(f\right)df};$均方根频率：$f_{bb}={\left(\frac{{\int }_0^{\infty }{f}^{2}p\left(f\right)df}{{\int }_0^{\infty }p\left(f\right)df}\right)}^{1/2};$频率方差：$f_v=\frac{{\int }_0^{\infty }{\left(f-f_{avg}\right)}^{2}p\left(f\right)df}{{\int }_0^{\infty }p\left(f\right)df};$频率标准差：$f_{rv}={\left[\frac{{\int }_0^{\infty }{(f-f_{avg})}^{2}p\left(f\right)df}{{\int }_0^{\infty }p\left(f\right)df}\right]}^2;$谱峰稳定指数：$S=\sqrt{\frac{\underset{i=1}{\overset{n/2}{\Sigma }}\left\{{f_i}^2·p\left(f_i\right)\right\}}{\underset{i=1}{\overset{n/2}{\Sigma }}p\left(f_i\right)}}/\sqrt{\frac{\underset{i=1}{\overset{n/2}{\Sigma }}\left\{{f_i}^4·p\left(f_i\right)\right\}}{\underset{i=1}{\overset{n/2}{\Sigma }}{f_i}^2·p\left(f_i\right)}};$第一频带相对能量：$E_{r1}={\int }_0^{B_f}p\left(f\right)df/{\int }_0^{F_s}p\left(f\right)df;$第二频带相对能量：$E_{r2}={\int }_{B_f}^{2B_f}p\left(f\right)df/{\int }_0^{F_s}p\left(f\right)df;$第三频带相对能量：$E_{r3}={\int }_{2B_f}^{3B_f}p\left(f\right)df/{\int }_0^{F_s}p\left(f\right)df;$第四频带相对能量：$E_{r4}={\int }_{3B_f}^{4B_f}p\left(f\right)df/{\int }_0^{F_s}p\left(f\right)df;$第五频带相对能量：$E_{r5}={\int }_{4B_f}^{5B_f}p\left(f\right)df/{\int }_0^{F_s}p\left(f\right)df;$所述24种特征指标计算公式中f表示信号的频率，p(f)表示信号的功率谱，B_f表示1/5频段值，F_s表示最高频率值；(3)寻求24种特征指标最佳投射方向矩阵的具体步骤为：①24种特征指标向量X^*的训练样本的预处理：设备处于某一工作状态，所述状态可为正常状态或某一典型故障状态，步骤(2)②所述24种特征指标的样本构成训练样本空间可以用p×n的矩阵X^*＝{x^*(i,j)|i＝1,2,...,p；j＝1,2,...,n}来表示，p表示训练样本的容量，n表示特征指标的数量；为消除各特征指标量纲的影响并统一其值的波动范围，需要对样本数据进行归一化处理：其中，x_max(j)，x_min(j)分别表示原始数据x^*(i,j)第j个指标的最大值和最小值，X＝{x(i,j)i＝1,2,...,p；j＝1,2,...,n}表示X^*＝{x^*(i,j)i＝1,2,...,p；j＝1,2,...,n}归一化处理后的序列；②将24种特征指标向量X^*进行投射，分析其投射值分布特点：所述投射就是把p维数据{x(i,j)|i＝1,2,...,p；j＝1,2,...,n}转化为以单位长度向量a＝{a(1),a(2),a(3),...,a(p)}为投射方向的一维投射值z(j)：$z\left(j\right)=\underset{i=1}{\overset{p}{\Sigma }}a\left(i\right)x(i,j),j=1,2,...,n,$其中a表示单位长度向量，然后按照投射值{z(j)|j＝1,2,...,n}的分布特点进行归类；③为了将所述步骤(3)②中各特征指标的投射值区域分开来，构造投射目标函数以寻求最佳的投射方向矩阵A＝[a₁,a₂,...,a_m]，投射目标函数指标可表达成：Q(a)＝S_zD_z，其中，S_z用投射值z(i)标准差来描述类间距离；D_z用投射值z(i)局部密度来描述类内密度，即$S_z=\sqrt{\frac{\underset{i=1}{\overset{n}{\Sigma }}{\left(z\right(i)-E(z\left)\right)}^2}{n-1}},$$D_z=\underset{i=1}{\overset{n}{\Sigma }}\underset{j=1}{\overset{n}{\Sigma }}(R-r(i,j\left)\right)·u(R-r(i,j\left)\right),$其中，E(z)为序列z＝{z(i)|i＝1,2,...,n}的平均值；R为局部密度的窗口半径，它的取值范围为：r_max＝max(r(i,j))，p为特征指标的数量；r(i,j)表示样本之间的距离，r(i,j)＝|z(i)‑z(j)|；单位阶跃函数$u\left(t\right)=\left\{\begin{array}{cc}1,& t\ge 0\\ 0,& t<0\end{array}\right.,$最大化目标函数：Max:Q(a)＝S_z·D_z，约束条件：$\underset{j=1}{\overset{p}{\Sigma }}{a}^2\left(j\right)=1,$在最佳投射方向矩阵A＝[a₁,a₂,...,a_m]下，每个特征指标的投射值w_i(i＝1,2,...,s)都相对集中成带状分布，借助所述带状分布特点衡量特征指标对设备状态识别的聚类性和敏感性；(4)利用所述投影值带状分布实现24种特征指标对煤矿设备状态识别的敏感度和聚类性的衡量和评价的具体步骤：①计算衡量各特征指标敏感度和聚类性的五种离散性指标：极差系数：${\mathrm{XD}}_{JCi}=\frac{{\mathrm{JC}}_i-JC}{JC}$平均差系数：${\mathrm{XD}}_{PJCi}=\frac{{\mathrm{PJC}}_i-PJC}{PJC}$离散系数：${\mathrm{XD}}_{LSi}=\frac{{\mathrm{LS}}_i-LS}{LS}$主轴线相对系数：${\mathrm{XD}}_{zi}=\frac{Z_{zi}-Z_z}{Z_z}$均值系数：${\mathrm{XD}}_{mi}=\frac{Z_{mi}-Z_m}{Z_m}$其中，极差JC＝max(z(i))‑min(z(i))；z(i)表示各特征指标在最佳投射方向矩阵上的投影值；平均差Z_z表示特征指标在最佳投射方向矩阵下的投影值主轴线；离散度投影值平均值JC、PJC、LS、Z_z、Z_m分别表示正常状态下的极差、平均差、离散度、主轴线和平均值；而JC_i、PJC_i、LS_i、Z_zi、Z_mi分别表示某一状态下的极差、平均差、离散度、主轴线和平均值；②判断特征指标的敏感性和聚类性如果步骤(4)①中计算得到的极差系数XD_JCi＞0.2、平均差系数XD_PJCi＞0.4、离散度系数XD_LSi＞0.3、主轴线相对系数XD_zi＞2和平均值系数XD_mi＞2，说明所述特征指标对煤矿设备的状态改变敏感性强，否则说明所述特征指标对煤矿设备的状态改变敏感性不强；当离散度LS＞1.5时，说明所述特征指标对煤矿设备的状态改变聚类性强，说明所述特征指标对煤矿设备的状态改变聚类性不强；在工程中可以针对不同设备的不同故障状态，分别计算得到相应的极差系数、平均差系数、离散度系数、主轴线相对系数和平均值系数，判断特征指标是否对所述故障状态具有较强的敏感性和聚类性。