多参数融合的球磨机振动特征提取方法 | [PDF全文] |
球磨机是物料粉磨的关键设备,由于矿石内部结构的不均匀性及外部环境的变化,球磨机难以保持在最佳工作状态[1, 2].在传统的破碎过程中,现场技术人员凭借经验来判断处理球磨机的内部运行状态,此基于手动操作的主观判断具有很大的局限性,且误差较大[3].目前,常用的方法是收集球磨机轴承,筒体,底座等的振动信号[4].根据振动信号与球磨机负荷之间的关系,确定球磨机的负荷状态[5, 6].由于球磨机磨矿过程是非线性、大延迟、强耦合的复杂过程,单一特征无法完整的反映球磨机负荷状态的变化.因此,将多个特征联合表征球磨机负荷状态的研究成为趋势[7].
由于振动信号的时域和频域特性在不同负载条件下不同.文中利用了时域特性和频域特性的优点,并结合D-S证据理论,提出了一种振动信号特征提取的多参数融合方法[8].通过对实验数据的分析,结果表明所提出的多参数融合的振动信号特征向量能够有效的反映球磨机的负荷状态,具有一定的优越性.
1 振动信号时、频域特征提取采用时域统计分析法,利用振动信号在时域上的波形、振幅、能量等信息,提取球磨机振动信号的有量纲特征参数(峰峰值、均值、标准差)和无量纲特征参数(峭度值、偏度值)[9-12].
振动信号的频谱特性可以有效地反映不同负载条件下振动信号的结构特征.通过统计方法描述振动信号的频域特性,以找到反映球磨机负载状态的频域特性[13-15].具体实现步骤分述如下.
1)线性差分方程用于表示观测数据的离散序列,如式(1)所示:
(1) |
式(1)中,w(n-i)表示白噪声序列.
2)对式(1)中的线性方程进行Z变换,并输入白噪声的谱密度
(2) |
式(2)中,
3)将bi=0、z=ejw代入式(2)可得AR模型的功率谱估计值为:
(3) |
式(3)中,P为AR模型的阶数;AR模型参数aj可以从观察到的数据序列的自相关序列中的计算导出.
2 多参数融合特征向量模型的构造分别测量球磨机在不同负荷条件下运行过程中轴承的振动信号[16-18],并分析各振动信号的时域和频域特征参数,构成N种球磨机负荷的振动信号的特征参数的一维矢量是:
(4) |
假设在最佳负荷状态下球磨机振动信号的特征参数指标矢量为:
(5) |
综合向量F和O可得mass函数向量:
(6) |
式(6)中,
同样,利用建立球磨机负荷特性预测模型的方法,可以建立n种负荷条件下l个振动特性参数的矩阵形式:
(7) |
基于D-S证据不确定性推理,对M进行多特征参数融合以获得得到新的mass函数向量MN[19-20]:
(8) |
式(8)中,第j个元素为:
对实验球磨机进行转速控制和能耗监测,并将传感器安装在球磨机两侧的支撑轴承座上进行振动信号采集,实验后的矿石进行筛分并分级处理[21].振动测试实验平台如图 1所示.
填充率分别为10%、20%、30%、40%和50%,料球比为0.6进行磨矿实验.对测取的振动信号进行时域、频域特征分析,如表 1所列.
点击放大 |
从表 1可知,不同负载振动信号的时域指标不同,可以很好地表征负载参数.时域指标峰峰值和均值在各个时段内发生变化,由于物料不断被磨细,振动信号的振幅减小,峰峰值和平均值更能反映振动信号的大小.时域指标方差、偏度和峭度更多的是描述振动信号的分布,随着磨矿的进行,其幅值也发生变化.不同负载的功率谱能量分布主要在频段一:1710~2848 Hz、频段二:2188~3855 Hz、频段三:3855~5830 Hz和频段四:5830~8686 Hz中发生变化.
根据对大量实验磨矿后矿石筛分与功率之比可得,填充率30%为最佳负荷状态下球磨机振动信号的特征参数指标矢量.由表 1中的数据可计算出将各个振动特征参数的mass函数值,结果如表 2所列.
点击放大 |
从表 2可以看出,在时域特征内,偏度呈现出先减小后增加的变化,其他时域特征参数主要反映在数值的微小波动中,虽然在不同负载条件下信号时域特性不同,但特征差异不是特别明显;在频域特征内,不同负荷状态功率谱总能量存在不同,负荷参数为(30,0.6)时总能量最小,负荷参数为(50,0.6)时总能量最大.但是, 负荷参数为(10,0.6)和负荷参数为(20,0.6)的总能量之间的差异却非常小,所以单独采用频域特征反映球磨机的负荷状态也可能存在误判.
为了解决单一特征值存在的缺陷,提高特征值反映球磨机负荷状态的准确性,利用D-S证据不确定性推理将振动信号的时域特征和频域特征进行数据融合,建立球磨机在不同负荷状态和振动特征参数之间的联系.融合计算出多振动特性参数的mass函数值,如表 3所列.
点击放大 |
由表 3可知,负荷参数为(30,0.6)时多特征参数融合后的特征向量的能量值最小,负荷为(40,0.6)时融合后的特征向量的能量最大,且不同负荷参数下融合后的特征向量的能量值存在明显的差异.对比表 2和表 3可得,不同负荷参数下信号特征存在差异,但是在单一特征的情况相邻负荷参数间振动信号的特征差异并不显著,而多参数融合后的特征在不同负荷参数下差异非常明显.
4 结论1)当填充率小于40%,在振动特征参数中,总能量、峰峰值、均值的mass函数值最大,对负荷预测贡献率较大.当填充率达到50%,有量纲时域特征的贡献率最小;随着填充率的增大,方差、频段一和频段二的mass函数值呈减小趋势,贡献率减小,偏差、频段三和频段四呈增大趋势,贡献率增大.
2)多参数融合后的各填充率下最大mass值的振动特征参数都不同,更加明显的表征不同负荷状态,避免了依靠单个特征建立特征向量模型的不准确性,为球磨机负荷的准确预测提供了新的方法.
[1] |
石立, 张国旺, 肖骁. 金属矿山选矿厂磨矿分级自动控制研究现状[J].
金属材料与冶金工程, 2014, 42(1): 43–48.
|
[2] |
梁朝霞.球磨机负荷检测软测量建模方法的研究[D].沈阳: 东北大学, 2008: 4-6.
http://cdmd.cnki.com.cn/Article/CDMD-10145-2010055713.htm |
[3] |
罗小燕, 陈慧明, 卢小江, 等. 基于网格搜索与交叉验证的SVM磨机负荷预测[J].
中国测试, 2017, 43(1): 132–135.
|
[4] |
胡智勇, 胡杰鑫, 谢里阳, 等. 滚动轴承振动信号处理方法综述[J].
中国工程机械学报, 2016, 14(6): 525–531.
|
[5] |
汤健, 赵立杰, 岳恒, 等. 基于多源数据特征融合的球磨机负荷软测量[J].
浙江大学学报(工学版), 2010, 44(7): 1406–1413.
DOI: 10.3785/j.issn.1008-973X.2010.07.031.
|
[6] |
周越, 司刚全, 曹晖, 等. 功率谱分析在筒式钢球磨煤机内存煤量测量中的应用研究[J].
工业仪表和自动化装置, 2006, 21(6): 26–30.
|
[7] |
田亮, 曾德良, 刘鑫屏, 等. 基于数据融合的球磨机最佳负荷工作点判断[J].
热能动力工程, 2004(2): 198–201.
DOI: 10.3969/j.issn.1001-2060.2004.02.024.
|
[8] |
阮竹恩, 李翠平, 李仲学. 地下突水隐患识别与风险预测的改进D-S算法[J].
金属矿山, 2014(9): 172–176.
|
[9] |
杨国亮, 康乐乐, 朱松伟, 等. 基于TLBO算法优化的球磨机FBEL控制方案研究[J].
江西理工大学学报, 2018, 39(1): 80–86.
|
[10] |
GUGEL K S, MOON R M. Automated mill control using vibration signal processing[C]//Charleston, USA: Cement Industry Technical Conference Record, IEEE, 2007.
https://www.researchgate.net/publication/224703321_Automated_mill_control_using_vibration_signal_processing |
[11] |
柴天佑. 复杂工业过程运行优化与反馈控制[J].
自动化学报, 2013, 39(11): 1744–1757.
|
[12] |
WEI D H, CRAIG I K. Grinding mill circuits-a survey of control and economic concerns[J].
International Journal of Mineral Process, 2009, 90(14): 56–66. |
[13] |
TANG J, CHAI T Y, YU W, et al. Feature extraction and selection based on vibration spectrum with application to estimating the load parameters of ball mill in grinding process[J].
Control Engineering Practice, 2012, 20(10): 991–1004. DOI: 10.1016/j.conengprac.2012.03.020. |
[14] |
任海.转速率与填充率对球磨机磨矿效果的影响研究[D].昆明: 昆明理工大学, 2011, 22-24.
http://cdmd.cnki.com.cn/Article/CDMD-10674-1012262179.htm |
[15] |
JAYASUNDARA C T, YANG R Y, YU A B, et al. Discrete particle simulation of particle in isa mill-effect of grinding medium properties[J].
Chemical Engineering Journal, 2008(135): 103–112. |
[16] |
畅晓亮, 董为民, 曹雪丽. 料球比对球磨机磨矿效果影响的试验研究[J].
新技术新工艺, 2012(2): 44–46.
DOI: 10.3969/j.issn.1003-5311.2012.02.015.
|
[17] |
吕权息, 汪思源, 张翔. 振动信号在球磨机料位监测系统的运用研究[J].
控制系统及其应用, 2002, 21(2): 36–39.
|
[18] |
GUGEL K, PALACIOS G. Improving ball mill control with modern tools based on digital signal processing(DSP) technology[C]// Dallas Cement Industry Technical Conference, IEEE-IAS/PCA, 2003.
|
[19] |
李时辉, 刘山葆. 一种新型数据融合方法在火灾监测中的应用研究[J].
江西理工大学学报, 2015, 36(3): 85–88.
|
[20] |
汤健, 柴天佑, 赵立杰, 等. 融合时/频信息的磨矿过程磨机负荷软测量[J].
控制理论与应用, 2012, 29(5): 564–570.
|
[21] |
代少升, 张辛, 杜江, 等. 磨音信号检测与处理方法研究[J].
计算机测量与控制, 2018, 43(11): 80–84.
|