CN112686144A - 一种基于磨音信号的矿石球磨过程负荷识别方法 - Google Patents
一种基于磨音信号的矿石球磨过程负荷识别方法 Download PDFInfo
- Publication number
- CN112686144A CN112686144A CN202011595567.XA CN202011595567A CN112686144A CN 112686144 A CN112686144 A CN 112686144A CN 202011595567 A CN202011595567 A CN 202011595567A CN 112686144 A CN112686144 A CN 112686144A
- Authority
- CN
- China
- Prior art keywords
- signal
- grinding
- noise
- load
- grinding sound
- 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.)
- Granted
Links
- 238000000227 grinding Methods 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000005236 sound signal Effects 0.000 title claims abstract description 42
- 230000008569 process Effects 0.000 title claims abstract description 16
- 238000000498 ball milling Methods 0.000 title claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims abstract description 33
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 24
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 238000007781 pre-processing Methods 0.000 claims abstract description 3
- 230000006870 function Effects 0.000 claims description 24
- 230000009467 reduction Effects 0.000 claims description 15
- 238000012549 training Methods 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 11
- 210000002569 neuron Anatomy 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 9
- 230000007613 environmental effect Effects 0.000 claims description 8
- 229920006395 saturated elastomer Polymers 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000036961 partial effect Effects 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 230000002829 reductive effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012300 Sequence Analysis Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 239000002184 metal Substances 0.000 abstract 1
- 230000001629 suppression Effects 0.000 abstract 1
- 238000009826 distribution Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000002401 inhibitory effect Effects 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000011410 subtraction method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种基于磨音信号的矿石球磨过程负荷识别方法,首先对采集到的原始信号进行预处理,包括消除直流分量、滤波,然后提出基于自回归(AR)谱估计的几何谱减法抑制噪声信号,提高信噪比;其次基于集合经验模态分解(EEMD)方法将磨音信号分解为K个本征模态分量IMF,利用排列熵选取可靠性高的模态分量进行磨音信号重构;接着计算每种负荷类型下重构信号的盒分形维数,并作为最终负荷分类依据;最后建立基于模糊C均值聚类的Bagging和极限学习机的负荷识别模型,实现负荷识别。本发明对国内某多金属选厂实际生产负荷进行识别,结果表明该识别模型能够准确地识别不同磨机负荷状态,具有良好的泛化性能和识别精度。
Description
技术领域
本发明涉及磨矿过程检测技术领域,具体涉及一种基于磨音信号的矿石球磨过程负荷识别方法。
背景技术
磨机负荷是指磨机内瞬时的全部装载量,包括新给矿量、循环负荷、水量及介质装载量等。磨机研磨物料过程的机理比较复杂,具有惯性大和时滞时间长等特性,使得磨机负荷最佳工作点具有时变性。但磨机内部状态黑箱,机理复杂难以针对该过程建立精确的模型。磨机内部负荷状态常出现三种状态:欠负荷、正常、过负荷,其中磨机在欠负荷状态时,不能充分发挥磨机的效率,而磨机处于过负荷状态时,会出现“涨肚”现象,使得磨机无法正常工作,甚至损坏磨机。因此,对磨机这种能量消耗巨大而能量利用率却极低的设备,准确检测磨机负荷状态是实现磨矿过程优化控制的一个关键。
为了实现磨矿过程的自动控制及安全运行,需要对磨机内部负荷状态进行检测。但因为磨机为全封闭重型旋转设备,其内部状态无法直接检测,因此只能通过其他参数对其进行间接的检测。磨机研磨声音是一个比较好能反应磨机内部状态的量。市场上也有一定的磨音分析产品,但由于各个选厂的应用场景不同,导致效果差强人意,国外产品存在价格昂贵不便维护等缺点,而国内难以提供有效准确的监测产品。现阶段,大多数厂家仍旧依靠经验丰富的领域专家依靠“听音”对磨机负荷状态进行判别,该方法实时性强,可以无损检测磨机负荷,但人工听音主观性强。因此,为磨机安全运转等信息监测提供一种无接触、无应激的自动化监测方法,对提高选矿厂生产效益具有重要的现实意义。
发明内容
本发明的目的是提供一种基于磨音信号的矿石球磨过程负荷识别方法。为实现上述目的,本发明的技术方案如下:
步骤1,基于多录音设备进行磨音信号检测。
步骤2,先对采集到的原始信号进行预处理,包括消除直流分量,滤波处理,再利用基于自回归谱估计地几何谱减法对两个录音设备的信号进行做差,达到抑制噪声信号的效果。
步骤3,基于集合经验模态分解(EEMD)方法将磨音信号分解为K个本征模态分量IMF,然后选择排列熵大的模态分量进行磨音信号重构,最后计算每种负荷类型下重构信号的分形维数,作为负荷识别模型的输入。
步骤4,进行基于模糊C均值聚类的Bagging和极限学习机(Bag-ELM)模型的搭建与训练,将步骤3得到的负荷分类特征输入到相应的Bag-ELM中进行识别。
进一步的,所述步骤1中,将录音设备一固定在离球磨机20cm-60cm的位置,主要采集球磨机研磨声音和环境噪声,采集到的时域信号表示为Y1(n);录音设备二背向磨机,固定在距离球磨机1.8m-2.2m的位置,主要采集环境噪声,采集到的时域信号表示为Y2(n),因而有:
其中,Y1(n)为磨音信号s(n)和环境噪声信号δ(n)的线性叠加,Y2(n)是环境噪声信号的近似值,n=1,2,...,N,N为时域信号采样点数量;
进一步的,所述步骤2的具体操作为:由式(1)可知现场采集的信号中包含干扰噪声,必须先进行适当的预处理,才能进行特征提取;为消除***直流漂移带来的影响,对原始信号做以下处理:
利用窗函数方法对磨音信号进行FIR滤波器的设计,设计方法如下:
确定窗函数Wd的类型;
根据待求滤波器的理想频率响应Hd(ejω)求出理想单位脉冲响应hd(n);
计算滤波器的单位脉冲响应h(n)
h(n)=hd(n)Wd(n) (4)其中,Wd(n)为窗函数序列。
FIR滤波器的输出为
滤波后从带噪磨音信号的功率谱中减去噪声功率谱,从而得到较为纯净的磨音功率谱,方法如下:
根据带噪磨音信号、纯净磨音信号和噪声信号之间在复平面上的几何关系,计算通用增益函数,其极坐标形式如下:
其中,{aY,aX,aD}和{θY,θX,θD}分别是带噪磨音信号、纯净磨音信号和噪声信号的幅值和相位;
得到定义增益函数为噪声信号与带噪磨音信号幅值比:
其中,定义CYD=cos(θY-θD),CXD=cos(θX-θD);
最终计算纯净磨音功率谱:
根据不相关随机序列的统计均值为零的原理,对M次分解得到的IMF分量做均值处理,以消除多次加入高斯白噪声对真实IMF的影响,得到EEMD分解的最终结果:
计算IMF分量与磨音信号之间的排列熵,对各个分量的排列熵进行由大到小排序,前三分之二的IMF分量保留,并作为有效的IMF分量;剩余IMF分量则当作干扰分量去除,根据选出的有效IMF分量进行重构,得到重构磨音信号s′(n);
把网格边长α逐步放大到nα,其中n属于正整数;
其中,i=1,2,…,Ns/n,n=1,2,…,N,N<Ns,Ns为重构信号的采样点数,d(nα)为中间计算算子;
令Nnα为S的网格计数,可表示为
Nnα=d(nα)/(nα)+1 (14)
其中Nnα≥1;
在lgnα-lgNnα中选择线性较好的一段为信号无标度区间[n1,n2],则此区间内应满足lgNnα=a lg nα+b,n1≤n≤n2;然后根据最小二乘法进行直线拟合,得出的斜率即为分形盒维数DB;
进一步的,步骤4中所述的基于模糊C均值聚类的Bag-ELM算法如下:
对数据集D={(xz,yz)},z=1,2,…,m,其中训练集为Dtrain={(xq,yq)},q=1,2,…,a,测试集为Dtest={(xp,yp)},p=1,2,…,b,对Dtrain进行boostrap自助采样,得到子训练集Dt={(xe,ye)},e=1,2,…,s,t=1,2,…,T,一共可采样出T个含s个样本的训练集,其中xq,xp和xe为重构磨音信号的分形盒维数,yp,yq和ye为负荷标签值;对磨音信号进行多尺度时间序列分析,将采样点数为H的时间序列,分割成τ个新时间序列,对每个新时间序列进行盒分形维数计算,即样本特征维数为τ;
采用极限学习机(ELM)作为基学习器的学习算法,则弱分类器ht为
其中,H为隐含层节点数,f为隐含层激活函数,βi为第i个隐层神经元与输出层之间的权重,ωi为第i个隐层神经元与输入层之间的权重,bi为第i个隐层神经元的阈值,ωi和bi可以随机生成;
按照投票结果对T个ELM分类器进行模糊C均值聚类分析,其优化目标函数为
其中,c为聚类中心个数,设为T的平方根,并向上取整,λ为模糊指数,通常设为2,uit和||ht-vi||分别为第t个分类器对第i个聚类中心的隶属度和欧氏距离,且必须满足约束条件
根据式(18)更新隶属度函数,根据式(19)更新聚类中心
最终得到c个聚类中心,从c个类别中分别选择一个ELM分类器组成新的集合,
所选分类器最接近各组的聚类中心,按照投票机制得出强分类器g:
所得式(20)为最终的负荷识别模型;在实际应用中根据步骤一到步骤三,得到尺度为τ的分形盒维数作为模型输入特征,将其通过c个ELM子模型得到识别结果,最后投票判断当前负荷状态。本发明提供的识别模型能够准确地识别不同磨机负荷状态,具有良好的泛化性能和识别精度。正常工况下的正确识别率为84%,识别率标准差为0.0952;偏空磨工况下的正确识别率为78%,识别率标准差为0.0831;空磨工况下的正确识别率为96%,识别率标准差为0.0121;偏饱磨磨工况下的正确识别率为92%,识别率标准差为0.0218;饱磨磨工况下的正确识别率为100%,识别率标准差为0.0031,表明该识别模型能够准确地识别磨机负荷状态,具有良好的泛化性能和识别精度。
附图说明
图1为本发明的总体流程图;
图2为几何谱减法分解示意图;
图3为降噪后空磨工况下的能量分布图;
图4为降噪后空磨工况下的声谱图;
图5为降噪后偏空磨工况下的能量分布图;
图6为降噪后偏空磨工况下的声谱图;
图7为降噪后正常工况下的能量分布图;
图8为降噪后正常工况下的声谱图;
图9为降噪后偏饱磨工况下的能量分布图;
图10为降噪后偏饱磨工况下的声谱图;
图11为降噪后饱磨工况下的能量分布图;
图12为降噪后饱磨工况下的声谱图。
具体实施方式
下面将结合附图说明本发明的具体实施方式,应当理解附图中示出和描述的实施方式仅是示例性的,意在阐述本发明的原理和方法,而并非限制本发明的范围。
如图1所示,本发明所述的是一种基于磨音信号的球磨过程负荷识别方法。首先,在磨机周围设置两个录音设备,采集多源信号,减少噪声的采集;再对采集到的原始信号进行基于AR谱估计的几何谱减法降噪处理,达到抑制噪声信号的效果;然后基于集合经验模态分解(EEMD)方法分解磨音信号,利用排列熵选取可靠性高的子信号进行磨音信号重构,并计算每种负荷类型下重构信号的分形维数,作为最终负荷分类指标;最后建立基于C均值聚类的Bagging和极限学习机(Bag-ELM)模型,将得到的负荷分类指标输入到相应的Bag-ELM中进行识别。此处采用国内某选厂的一段球磨流程的磨音数据,并对其进行样本切割,共得到2400组工业历史数据进行分析,其中将1900组数据作为训练数据,剩余500组作为测试数据。
步骤1,基于多录音设备进行信号检测,将录音设备一固定在离球磨机20cm-60cm的位置,主要采集球磨机研磨声音和环境噪声,采集到的时域信号表示为Y1(n);录音设备二背向磨机,固定在距离球磨机1.8m-2.2m的位置,主要采集环境噪声,采集到的时域信号表示为Y2(n),因而有:
其中,Y1(n)为磨音信号s(n)和环境噪声信号δ(n)的线性叠加,Y2(n)是环境噪声信号的近似值,n=1,2,...,N,N为时域信号采样点数量;
步骤2,由式(1)可知现场采集的信号中包含干扰噪声,必须先进行适当的处理,才能进行特征提取,主要包括消除直流分量,滤波处理和基于AR谱估计的几何谱减法降噪;
为消除***直流漂移带来的影响,对原始信号做以下处理:
由于磨音和故障信号的频率一般都不低于500Hz,所以500Hz以下的信号可以是认为磨机运转噪声引起,设计一个截止频率为500Hz的FIR高通滤波器将其直接滤去,利用窗函数方法对磨音信号进行FIR滤波器的设计,其性能指标为:通带截止频率ft为550Hz,阻带截止频率fz为500Hz,过渡带fg为50Hz,最小阻带衰减As为50dB,滤波器的设计方法如下:
选择汉宁窗作为窗函数Wd;
根据待求滤波器的理想频率响应Hd(ejω)求出理想单位脉冲响应hd(n);
计算滤波器的单位脉冲响应h(n)
h(n)=hd(n)Wd(n) (4)
滤波器的输出为
几何谱减法将复平面中的带噪信号频谱表示为纯净信号和噪声信号的矢量和,可以为谱分解方法提供有价值的见解,可以简单理解为从带噪磨音信号的功率谱中减去噪声功率谱,从而得到较为纯净的磨音功率谱,方法如下:
根据带噪磨音信号、纯净磨音信号和噪声信号之间在复平面上的几何关系,如图2所示,计算通用增益函数,其极坐标形式如下:
其中,{aY,aX,aD}和{θY,θX,θD}分别是带噪磨音信号、纯净磨音信号和噪声信号的幅值和相位;
得到定义增益函数为噪声信号与带噪磨音信号幅值比:
其中,定义CYD=cos(θY-θD),CXD=cos(θX-θD);
最终计算纯净磨音功率谱:
降噪后对各个工况下的磨音进行可视化,如图3至图12所示,对应不同的负荷参数,灰度越重代表该频率段能量越集中,磨音信号的频谱主要集中在500Hz-6000Hz之间。空磨工况下能量分布整体较高;正常工况下能量分布相对均匀且整体偏低,饱磨工况下低频段能量分布较集中,高频段能量较低;偏空磨和偏饱磨则介于三者之间。现实情况中无法确定初始的信噪比,为了体现降噪效果,人为产生一定信噪比的带噪磨音。因此按照设定的信噪比产生高斯白噪声叠加到原始磨音信号中。表1为原始磨音信号与降噪后磨音信号的信噪比,可以观察到利用基于AR谱的几何谱减法降噪后磨音信号有了明显的改善;
步骤3,基于集合经验模态分解(EEMD)方法将磨音信号分解为若干个本征模态分量IMF,每个本征模态分量代表具有物理含义的单一振动模态,然后选择排列熵大的模态分量进行磨音信号重构,最后计算每种负荷类型下重构信号的多尺度分形维数,并作为最终负荷分类指标;
EEMD分解过程如下:
根据不相关随机序列的统计均值为零的原理,对10次分解得到的IMF分量做均值处理,以消除多次加入高斯白噪声对真实IMF的影响,得到EEMD分解的最终结果:
计算IMF分量与磨音信号之间的排列熵,对各个分量的排列熵进行由大到小排序,前三分之二的IMF分量保留,并作为有效的IMF分量;剩余IMF分量则当作干扰分量去除,根据选出的有效IMF分量进行重构,得到重构磨音信号s′(n);
将重构信号序列作无量纲图形化处理可得到一条自仿射曲线,具有不规则性和近似的自相似性,可以用分形维数来描述其复杂程度,其中盒维数计算简单,能从全局反映信号的复杂程度和非线性,适用于周期性的磨音信号的描述,盒维数是指覆盖一个信号的最小盒子数,盒子的尺度不同,信号的盒维数不同,一维离散信号的分形盒维数的大小在1和2之间,信号的复杂程度越高,分形盒维数越大;
把网格边长α逐步放大到nα,其中n属于正整数;
其中,i=1,2,…,Ns/n,n=1,2,…,N,N<Ns,Ns为采样点数,d(nα)为中间计算算子;
令Nnα为S的网格计数,可表示为
Nnα=d(nα)/(nα)+1 (14)
其中Nnα≥1;
在lg nα-lg Nnα中选择线性较好的一段为信号无标度区间[n1,n2],则此区间内应满足lg Nnα=a lg nα+b,n1≤n≤n2;然后根据最小二乘法进行直线拟合,得出的斜率即为分形盒维数DB;
步骤4,进行基于模糊C均值聚类的Bagging和极限学习机(Bag-ELM)模型的搭建与训练,将步骤3得到的负荷分类特征输入到相应的Bag-ELM中进行识别;
对数据集D={(xz,yz)},z=1,2,…,2400进行划分,其中训练集为Dtrain={(xq,yq)},q=1,2,…,1900,测试集为Dtest={(xp,yp)},p=1,2,…,500,对Dtrain进行boostrap自助采样,得到子训练集Dt={(xe,ye)},e=1,2,…,190,t=1,2,…,10,一共可采样出10个含190个样本的训练集,其中xq,xp和xe为重构磨音信号的特征,为增强数据的可靠性,将10000个采样点作为一个样本进行分割,对每1024个采样点的序列计算盒分形维数,共得到9个分形盒维数,因此样本的输入特征为9维,yp,yq和ye为负荷标签值;
采用极限学习机(ELM)作为基学习器的学习算法,则弱分类器ht为
其中,H为隐含层节点数,f为隐含层激活函数,βi为第i个隐层神经元与输出层之间的权重,ωi为第i个隐层神经元与输入层之间的权重,bi为第i个隐层神经元的阈值,ωi和bi可以随生成;
按照投票结果对T个ELM分类器进行C均值聚类分析,其优化目标函数为
其中,c为聚类中心个数,取值为5,λ为模糊指数,通常设为2,uit和||ht-vi||分别为第t个分类器对第i个聚类中心的隶属度和欧氏距离,且必须满足约束条件
根据式(18)更新隶属度函数,根据式(19)更新聚类中心
最终得到5个聚类中心,分别从5个类别中分别选择的一个ELM分类器组成新的集合,所选分类器最接近各组的聚类中心,按照投票机制得出强分类器g:
所得式(20)为最终的负荷识别模型;在实际应用中,根据步骤一到步骤三,得到尺度为9的多尺度分形盒维数作为模型输入特征,将其通过5个ELM子模型得到识别结果,最后投票判断当前负荷状态。
用本发明方法对五种工况进行识别,模型性能统计结果见表2,总体的正确识别率为90%,Kappa系数为0.8441;正常工况下的正确识别率为84%;偏空磨工况下的正确识别率为78%;空磨工况下的正确识别率为96%;偏饱磨磨工况下的正确识别率为92%;饱磨磨工况下的正确识别率为100%,表明该识别模型能够准确地识别磨机负荷状态,具有良好的泛化性能和识别精度。
表1
表2
Claims (2)
1.一种基于磨音信号的矿石球磨过程负荷识别方法,其特征在于包括以下步骤:
步骤1,基于多录音设备进行音频信号检测,将录音设备一固定在离球磨机一段距离的位置,采集球磨机研磨声音和环境噪声,采集到的时域信号表示为Y1(n);录音设备二背向磨机,固定在距离球磨机一段距离的位置,采集环境噪声,采集到的时域信号表示为Y2(n),有:
其中,Y1(n)为磨音信号s(n)和环境噪声信号δ(n)的线性叠加,Y2(n)是环境噪声信号的近似值,n=1,2,…,N,N为时域信号采样点数量;
步骤2,对采集到的时域信号进行预处理以抑制噪声信号,首先消除直流分量:
然后设计一个FIR高通滤波器对时域信号进行磨机固有频率滤波降噪,滤波器的设计方法如下:
确定窗函数Wd的类型;
根据待求滤波器的理想频率响应Hd(ejω)求出理想单位脉冲响应hd(n);
计算滤波器的单位脉冲响应h(n)
h(n)=hd(n)Wd(n) (4)
其中,Wd(n)为窗函数序列。
滤波器的输出为
滤波后从带噪磨音信号的功率谱中减去噪声功率谱,从而得到纯净的磨音功率谱,方法如下:
根据带噪磨音信号、纯净磨音信号和噪声信号之间在复平面上的几何关系,计算通用增益函数,其极坐标形式如下:
其中,{aY,aX,aD}和{θY,θX,θD}分别是带噪磨音信号、纯净磨音信号和噪声信号的幅值和相位;
得到定义增益函数为噪声信号与带噪磨音信号幅值比:
其中,定义CYD=cos(θY-θD),CXD=cos(θX-θD);
最终计算纯净磨音功率谱:
步骤3,基于集合经验模态分解方法,即EEMD方法,将磨音信号分解为K个本征模态分量IMF,然后选择排列熵大的模态分量进行磨音信号重构,最后计算每种负荷类型下重构信号的多尺度分形维数,并作为最终负荷分类指标;
EEMD分解过程如下:
计算IMF分量与磨音信号之间的排列熵,对各个分量的排列熵进行由大到小排序,前三分之二的IMF分量保留,并作为有效的IMF分量;剩余IMF分量则当作干扰分量去除,根据选出的有效IMF分量进行重构,得到重构磨音信号s′(n);
分形维数计算方法如下:
把网格边长α逐步放大到nα,其中n属于正整数;
其中,i=1,2,…,Ns/n,n=1,2,…,N,N<Ns,Ns为重构信号的采样点数,d(nα)为中间计算算子;
令Nnα为S的网格计数,可表示为
Nnα=d(nα)/(nα)+1 (14)
其中Nnα≥1;
在lgnα-lgNnα中选择线性较好的一段为信号无标度区间[n1,n2],则此区间内满足lgNnα=algnα+b,n1≤n≤n2;然后根据最小二乘法进行直线拟合,得出的斜率即为分形盒维数DB;
步骤4,进行基于模糊C均值聚类的Bagging和极限学习机ELM模型的搭建与训练,将步骤3得到的负荷分类特征输入到相应的Bag-ELM中进行识别,基于模糊C均值聚类的Bagging-ELM算法如下:
对数据集D={(xz,yz)},z=1,2,…,m,其中训练集为Dtrain={(xq,yq)},q=1,2,…,a,测试集为Dtest={(xp,yp)},p=1,2,…,b,对Dtrain进行boostrap自助采样,得到子训练集Dt={(xe,ye)},e=1,2,…,s,t=1,2,…,T,一共采样出T个含s个样本的训练集,其中xq、xp和xe为重构磨音信号的分形盒维数,yp、yq和ye为负荷标签值;对磨音信号进行多尺度时间序列分析,将采样点数为H的时间序列,分割成τ个新时间序列,对每个新时间序列进行盒分形维数计算,即样本特征维数为τ;
采用极限学习机,即ELM,作为基学习器的学习算法,则弱分类器ht为
其中,H为隐含层节点数,f为隐含层激活函数,βi为第i个隐层神经元与输出层之间的权重,ωi为第i个隐层神经元与输入层之间的权重,bi为第i个隐层神经元的阈值,ωi和bi随机生成;
按照投票结果对T个ELM分类器进行模糊C均值聚类分析,其优化目标函数为
其中,c为聚类中心个数,设为T的平方根,并向上取整,λ为模糊指数,设为2,uit和||ht-vi||分别为第t个分类器对第i个聚类中心的隶属度和欧氏距离,且必须满足约束条件
根据式(18)更新隶属度函数,根据式(19)更新聚类中心
最终得到c个聚类中心。从c个类别中分别选择一个ELM分类器组成新的集合,所选分类器最接近各组的聚类中心,按照投票机制得出强分类器g:
式(20)为最终的负荷识别模型;在实际应用中根据步骤一到步骤三,得到尺度为τ的多尺度分形盒维数作为模型输入特征,将其通过c个ELM子模型得到识别结果,最后投票判断当前负荷状态。
2.如权利要求1所述的基于磨音信号的矿石球磨过程负荷识别方法,其特征在于:根据专家经验划分五种工况为空磨、偏空磨、正常、偏饱磨、饱磨,录音设备一固定在离球磨机的距离为20cm-60cm,录音设备二背向磨机,固定在距离球磨机1.8m-2.2m的位置,分别采集五种工况下的磨音信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011595567.XA CN112686144B (zh) | 2020-12-29 | 2020-12-29 | 一种基于磨音信号的矿石球磨过程负荷识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011595567.XA CN112686144B (zh) | 2020-12-29 | 2020-12-29 | 一种基于磨音信号的矿石球磨过程负荷识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112686144A true CN112686144A (zh) | 2021-04-20 |
CN112686144B CN112686144B (zh) | 2022-03-08 |
Family
ID=75455283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011595567.XA Active CN112686144B (zh) | 2020-12-29 | 2020-12-29 | 一种基于磨音信号的矿石球磨过程负荷识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112686144B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113838054A (zh) * | 2021-11-29 | 2021-12-24 | 山东精良海纬机械有限公司 | 基于人工智能的机械零件表面损伤检测方法 |
CN114077851A (zh) * | 2021-11-22 | 2022-02-22 | 河北工业大学 | 基于fsvc的球磨机工况识别方法 |
CN115602191A (zh) * | 2022-12-12 | 2023-01-13 | 杭州兆华电子股份有限公司(Cn) | 一种变压器声纹检测***的噪声消除方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110251807A1 (en) * | 2009-01-26 | 2011-10-13 | Geneva Cleantech Inc. | Automatic detection of appliances |
CN105138782A (zh) * | 2015-09-02 | 2015-12-09 | 上海大学 | 基于eemd-elm的非平稳脉动风速高精度预测方法 |
CN105528636A (zh) * | 2015-12-04 | 2016-04-27 | 中国人民解放军61599部队计算所 | 一种基于模糊推理的磨机负荷参数软测量方法 |
CN110111012A (zh) * | 2019-05-13 | 2019-08-09 | 中南大学 | 一种基于稳态图像特征的接触网负荷识别方法 |
CN110470477A (zh) * | 2019-09-19 | 2019-11-19 | 福州大学 | 一种基于ssae和ba-elm的滚动轴承故障诊断方法 |
CN111259869A (zh) * | 2020-03-11 | 2020-06-09 | 北京慧飒科技有限责任公司 | 一种非侵入式电气集群负载故障识别方法 |
EP3685137A1 (en) * | 2017-09-23 | 2020-07-29 | Nanoprecise SCI Corp. | System and method for automated fault diagnosis and prognosis for rotating equipment |
-
2020
- 2020-12-29 CN CN202011595567.XA patent/CN112686144B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110251807A1 (en) * | 2009-01-26 | 2011-10-13 | Geneva Cleantech Inc. | Automatic detection of appliances |
CN105138782A (zh) * | 2015-09-02 | 2015-12-09 | 上海大学 | 基于eemd-elm的非平稳脉动风速高精度预测方法 |
CN105528636A (zh) * | 2015-12-04 | 2016-04-27 | 中国人民解放军61599部队计算所 | 一种基于模糊推理的磨机负荷参数软测量方法 |
EP3685137A1 (en) * | 2017-09-23 | 2020-07-29 | Nanoprecise SCI Corp. | System and method for automated fault diagnosis and prognosis for rotating equipment |
CN110111012A (zh) * | 2019-05-13 | 2019-08-09 | 中南大学 | 一种基于稳态图像特征的接触网负荷识别方法 |
CN110470477A (zh) * | 2019-09-19 | 2019-11-19 | 福州大学 | 一种基于ssae和ba-elm的滚动轴承故障诊断方法 |
CN111259869A (zh) * | 2020-03-11 | 2020-06-09 | 北京慧飒科技有限责任公司 | 一种非侵入式电气集群负载故障识别方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114077851A (zh) * | 2021-11-22 | 2022-02-22 | 河北工业大学 | 基于fsvc的球磨机工况识别方法 |
CN114077851B (zh) * | 2021-11-22 | 2024-04-23 | 河北工业大学 | 基于fsvc的球磨机工况识别方法 |
CN113838054A (zh) * | 2021-11-29 | 2021-12-24 | 山东精良海纬机械有限公司 | 基于人工智能的机械零件表面损伤检测方法 |
CN113838054B (zh) * | 2021-11-29 | 2022-02-25 | 山东精良海纬机械有限公司 | 基于人工智能的机械零件表面损伤检测方法 |
CN115602191A (zh) * | 2022-12-12 | 2023-01-13 | 杭州兆华电子股份有限公司(Cn) | 一种变压器声纹检测***的噪声消除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112686144B (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112686144B (zh) | 一种基于磨音信号的矿石球磨过程负荷识别方法 | |
Zhu et al. | Wavelet analysis of sensor signals for tool condition monitoring: A review and some new results | |
Wu et al. | Induction machine fault detection using SOM-based RBF neural networks | |
CN102750543B (zh) | 一种基于bud谱峭度的暂态电能质量扰动分类识别方法 | |
CN113780055B (zh) | 一种momeda与压缩感知的滚动轴承故障诊断方法 | |
CN103398769B (zh) | 基于小波ggd特征和不均衡k-均值下采样集成svm的变压器在线故障检测方法 | |
CN111879397B (zh) | 高压断路器储能机构故障诊断方法 | |
Yan et al. | A feature selection framework-based multiscale morphological analysis algorithm for fault diagnosis of rolling element bearing | |
CN112326210A (zh) | 一种声振信号联合1d-cnn的大型电机故障诊断方法 | |
CN107728018A (zh) | 一种电力电缆现场局部放电信号的降噪方法 | |
CN111414893B (zh) | 基于vmd精细复合多尺度散布熵的转子故障特征提取方法 | |
CN106153179A (zh) | 中速磨煤机振动故障诊断方法 | |
CN109605128B (zh) | 一种基于功率谱熵差的铣削颤振在线检测方法 | |
Zhao et al. | Probabilistic principal component analysis assisted new optimal scale morphological top-hat filter for the fault diagnosis of rolling bearing | |
CN116304570B (zh) | 一种基于EEMD联合Chebyshev滤波的水轮机故障信号去噪方法 | |
Zheng et al. | An adaptive group sparse feature decomposition method in frequency domain for rolling bearing fault diagnosis | |
Mubaraali et al. | Intelligent fault diagnosis in microprocessor systems for vibration analysis in roller bearings in whirlpool turbine generators real time processor applications | |
CN114139820A (zh) | 一种用于非侵入式电能负荷预测的改进模态分解方法 | |
Ding et al. | A fault feature extraction method of motor bearing using improved LCD | |
CN113988125A (zh) | 一种基于改进同步压缩变换的扭振信号瞬时频率提取方法 | |
Zhang et al. | A novel hybrid compound fault pattern identification method for gearbox based on NIC, MFDFA and WOASVM | |
Cui et al. | Egram based SVD method for gear fault diagnosis | |
CN114897016A (zh) | 一种基于多源频谱特征的风机轴承故障智能诊断方法 | |
Chen et al. | A novel two-stage feature extraction based on multi-fractal spectrum with CI-VMD for early fault diagnosis of bearing | |
Li | Recognition method of non-stationary mechanical vibration signal based on convolution neural network |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20220217 Address after: Yuelu District City, Hunan province 410083 Changsha Lushan Road No. 932 Applicant after: CENTRAL SOUTH University Applicant after: Changsha Mining and Metallurgy Research Institute Co., Ltd Address before: Yuelu District City, Hunan province 410083 Changsha Lushan Road No. 932 Applicant before: CENTRAL SOUTH University |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |