CN113537579A - 一种集成csa和h-elm的水泥生产质量预测方法 - Google Patents
一种集成csa和h-elm的水泥生产质量预测方法 Download PDFInfo
- Publication number
- CN113537579A CN113537579A CN202110721214.8A CN202110721214A CN113537579A CN 113537579 A CN113537579 A CN 113537579A CN 202110721214 A CN202110721214 A CN 202110721214A CN 113537579 A CN113537579 A CN 113537579A
- Authority
- CN
- China
- Prior art keywords
- learning machine
- extreme learning
- crow
- elm
- hidden layer
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/06—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
- G06N3/061—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons using biological neurons, e.g. biological neurons connected to an integrated circuit
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/088—Non-supervised learning, e.g. competitive learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Software Systems (AREA)
- Strategic Management (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Economics (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Neurology (AREA)
- Game Theory and Decision Science (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种集成CSA和H‑ELM的水泥生产质量预测方法,包括以下步骤:S1:采集与水泥生产质量预测相关的数据;S2:使用灰色关联分析法对步骤S1收集到的数据进行初步筛选;S3:初始化分层极限学习机的参数,利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型;S4:利用所述预测模型对水泥生产质量进行预测。本方面通过CSA算法优化了H‑ELM隐含层神经元个数,舍去人工调参的繁琐操作,提高了模型的实用性和预测精度,同时优化后的模型可以实时的,智能和高精度的预测水泥烧成***回转窑中产品的生产质量。
Description
技术领域
本发明涉及水泥生产技术领域,更具体地,涉及一种集成CSA和H-ELM的水泥生产质量预测方法。
背景技术
水泥与塑料、金属、木材俗称国民经济四大材料。水泥工业是国民经济和社会发展的重要支柱产业,是我国经济持续增长的重要支撑力量,其自动化程度直接彰显一个国家的整体基础工业实力。随着国家大力推动城市化建设,加快工业发展进度,截止到2016年,我国水泥年平均产量在25亿吨左右,约占全球水泥产量五成以上。但我国的水泥工业相对自动化程度较低,单位水泥熟料能耗比发达国家高50%以上,并且在水泥生产过程会产生大量的气态污染物以及粉尘颗粒,与此同时,水泥的生产质量还无法得到保证。因此,在保证水泥生产线稳定工作的前提下,如何提升自动化程度,水泥生产质量,降低水泥生产能耗,减少污染物排放成为了一个社会焦点问题。
我国水泥工业主要采用新型干法水泥生产技术。主要由生料制备***,烧成***和水泥制备***组成。烧成***是水泥生产过程中最为重要的工艺环节,在烧成***中,水泥生料会发生复杂的化学物理反应,其过程直接影响水泥生产质量高低以及能源利用率。游离氧化钙(下文简写为f-CaO)含量是评价烧成***生产熟料质量的关键性指标,其含量高低对水泥的生产效率以及能耗有着直接性影响。在水泥生产企业,工作人员需要及时了解f-CaO含量,以便调整工况。现阶段f-CaO含量检测,主要是使用化学分析法以及荧光分析法。然而,化学分析法需要先对熟料进行采样,然后实验室检测,不能及时的反映水泥生产质量。荧光分析法使用的设备价格较高,实际水泥生产企业出于成本的考虑并没有大规模装备,因此使用率极低。水泥企业主要采用化学分析法,当工作人员发现f-CaO含量测量值不符合预期时,再根据自身经验手动调整水泥生产设备的某些重要的工艺参数已达到对f-CaO含量的控制,但此方法会导致调整参数操作存在滞后性,难以满足实际生产工况实时性要求,同时工作人员依照自身经验进行调整,具有一定的盲目性,若调整不佳,反而影响水泥生产设备能耗,污染物排放,质量等。因此,有必要提出一种全新的f-CaO含量预测方法,在保障预测准确的前提下,降低预测时滞,提升产品质量。不仅可以对水泥生产过程变量具有更多的掌控,而且对生产进行指导,具有重要意义。
人工智能技术的发展,为水泥生产质量预测提供了新的研究方向,通过建立工艺参数的预测模型,达到对f-CaO含量实时检测的目的,还可以有效的弥补工业参数检测方法的不足。因此,研究f-CaO含量预测模型以及优化方法对水泥生产应用技术具有重要的价值。
现有技术中主要通过针对建立工艺参数的预测模型实现对f-CaO含量预测,然而,随着水泥生产要求的逐年提升,工艺越来越复杂,传统的人工智能算法已经无法满足使用要求,逐渐暴露出参数选取不合理,泛化能力差等缺点。因此,采用单一的预测模型,由于自身的局限性,往往会出现预测效果不佳的情况。针对单一模型的局限性,许多的学者提出智能算法优化神经网络参数来提升f-CaO含量预测性能,然而,由于目前群体优化算法(如上述GA,GWO,PSO算法)在全局搜索,局部搜索上依然采用随机的方式,影响精度的参数过多,从而使神经网络的预测性能难以满足水泥生产技术的高要求。
公开日为2019年02月15日,公开号为CN109342703A的中国专利公开了一种水泥熟料游离钙含量测量方法及***。所述测量方法包括:获取与所述水泥熟料游离钙FCaO含量有关的相关变量参数;根据所述相关变量参数,确定每个所述相关变量参数的时间序列;根据所述时间序列建立基于多变量时间序列的卷积残差网络的熟料FCaO软测量模型;根据所述熟料FCaO软测量模型确定水泥熟料游离钙含量。该专利使用单一的预测模型,同样会出现预测效果不佳的情况。
发明内容
本发明提供一种集成CSA和H-ELM的水泥生产质量预测方法,利用乌鸦算法对预测模型进行优化,达到提升预测精度的目的,解决水泥烧成***回转窑水泥质量的实时预测问题。
为解决上述技术问题,本发明的技术方案如下:
一种集成CSA和H-ELM的水泥生产质量预测方法,包括以下步骤:
S1:采集与水泥生产质量预测相关的数据;
S2:使用灰色关联分析法对步骤S1收集到的数据进行初步筛选;
S3:初始化分层极限学习机(简称H-ELM)的参数,利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法(简称CSA)优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型;
S4:利用所述预测模型对水泥生产质量进行预测。
优选地,步骤S1中所述与水泥生产质量预测相关的数据包括与烧成***回转窑相关的喂料量、煤粉流量、分解炉出口温度、物料温度、二次风温、压力、回转窑电机电流、转速、高温风机转速、回转窑尾温、回转窑头负压、回转窑尾压力、回转窑烟室氧气含量、烟室氮氧化物含量、三次风温和三次风压力。
优选地,步骤S2中利用灰色关联分析法对步骤S1收集到的数据进行初步筛选,具体为:
设自变量为X=(xi(l)|i=1,2,...,h,l=1,2,...,M),因变量为T={t(l)|l=1,2,...,M},其中h表示自变量的数目,M表示样本的数目,令Δi(l)=|t(l)-xi(l)|,则灰色关联度分析法计算得到的灰关联度为:
式中,ρ为分辨系数,0<ρ<1;
因变量和自变量的灰色关联度γi如下:
设置灰色关联度阈值,选择灰色关联度大于灰色关联度阈值的数据作为数据样本,完成初步筛选。
优选地,所述步骤S3中初始化分层极限学习机的参数具体为:
随机初始化隐含层个数、神经元个数、隐含层连接权值和阈值。
优选地,步骤S3中利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型,具体为:
a、将分层极限学习机的隐层神经元个数组成乌鸦算法的初始种群集合,初始化乌鸦算法的参数;
b、以分层极限学习机的训练结果误差作为适应度函数,评估乌鸦的位置;
c、择优选择乌鸦种群集合中乌鸦的新位置,评估乌鸦的新位置,如果新位置比当前位置适应度值高,更新乌鸦位置至新位置;否则,乌鸦的新位置随机更新;
d、重复步骤c直到迭代次数达到阈值,迭代终止,输出当前最佳位置,得到分层极限学习机的最优隐层神经元个数;
e、分层极限学习机利用FISTA算法计算每一隐含层网络权重,计算该层的实际输出并且作为下一层输入;
f、运用最小二乘法求得分层极限学习机最后一层隐层和输出层之间的权值,输出预测结果。
优选地,步骤c中乌鸦的新位置更新,具体如下:
择优选择乌鸦种群集合中乌鸦的新位置:
Xi,iter+1=Xi,iter+ri×fl×(Mj,iter-Xi,iter)
式中,向量Xi,iter表示第i只乌鸦Crowi在第i次迭代中的位置,Xi,iter+1表示第i只乌鸦Crowi在第i+1次迭代中的位置,Mj,iter表示第j只乌鸦Crowi在第i次的记忆,ri为0-1之间的随机数,fl为乌鸦移动的步长;
乌鸦的新位置随机更新:
Xi,iter=a random position
Mi,iter更新为:
式中,f()表示乌鸦算法的适应度函数。
优选地,步骤e中所述分层极限学习机,具体为:
所述分层极限学习机是一种具有多层特性的结构,其分层体现在整个模型的训练过程分为无监督特征提取和有监督特征分类。
优选地,在无监督特征提取阶段,基于极限学习机自动编码器ELM-AE思想,并使用基于l1-范数范数最小化的快速阈值收缩算法优化隐含层权值,分类部分采用基于岭回归的ELM模型作为分类模型;
在无监督学习阶段,分层极限学习机利用基于极限学习机自动编码器ELM-AE作为特征提取工具,ELM-AE中,输出等于输入,隐藏节点通过随机正交权重和偏差将输入数据投影到不同或相等的维度,实现原始数据在不同维度的特征表达;
基于极限学习机自动编码器ELM-AE训练过程可以用以下优化公式表述为:
式中,H表示隐含层输出矩阵,X表示每一层的输入矩阵,β表示输出权值;
采用l1-范数对β进行约束,此时基于极限学习机自动编码器ELM-AE可以表示为如下约束问题:
为了明确l1优化问题,可以将上式从写成公式:
Oβ=p(β)+q(β)
其中,p(β)=||Hβ-X||2,q(β)=||β||l1,再使用基于l1-范数范数最小化的快速阈值收缩算法,即FISTA算法优化隐含层权值;
经过n层的无监督学习,最终可获得高层稀疏特征,每个隐含层的输出可以表示为:
Hi=g(Hi-1·β)
式中,Hi代表第i隐含层输出矩阵,Hi-1代表第i-1隐含层输出矩阵,g(·)代表隐藏层的激活函数。
优选地,所述FISTA算法优化隐含层权值:
F1:计算光滑凸函数P的利普希兹常量L;
F2:初始化参数yl=β0,tl=1,y为隐层输出权重,t为迭代次数,然后对于k(k≥1)进行如下迭代:
βk=pL(yk)
其中pL可以由下式给出:
其中β代表隐含层输出权重,βk-1代表上一次迭代的隐含层输出权重;
F3:通过上述迭代,计算出输出权值β。
优选地,在有监督特征分类阶段,利用最小二乘法计算最后一层的隐含层输出矩阵。
与现有技术相比,本发明技术方案的有益效果是:
本发明通过CSA算法优化了H-ELM隐含层神经元个数,舍去人工调参的繁琐操作,提高了模型的实用性和预测精度,同时优化后的模型可以实时的,智能和高精度的预测水泥烧成***回转窑中产品的生产质量,克服了传统单一预测模型参数选取不合理,泛化能力差,预测精度不高的缺点。
附图说明
图1为本发明的方法流程图。
图2为本发明提出的模型训练集与真实值对比示意图。
图3为本发明提出的模型测试集与真实值对比示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
本实施例提供一种集成CSA和H-ELM的水泥生产质量预测方法,如图1所示,包括以下步骤:
S1:采集与水泥生产质量预测相关的数据;
S2:使用灰色关联分析法对步骤S1收集到的数据进行初步筛选;
S3:初始化分层极限学习机的参数,利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型CSA-H-ELM;
S4:利用所述预测模型对水泥生产质量进行预测。
通过了解整个水泥的生产流程可以知道,烧成***回转窑末端熟料的f-CaO含量是衡量水泥质量的关键性指标,其含量主要由生原料的配比以及烧成***加工情况决定。生料的配比严格按照要求进行,所以生料的配比相对而言是稳定的。而整个烧成***回转窑内,工况复杂,能够直接或间接影响烧成质量的工艺参数众多,所以在采集数据阶段,重点选取烧成***回转窑的工艺参数。步骤S1中所述与水泥生产质量预测相关的数据包括与烧成***回转窑相关的喂料量、煤粉流量、分解炉出口温度、物料温度、二次风温、压力、回转窑电机电流、转速、高温风机转速、回转窑尾温、回转窑头负压、回转窑尾压力、回转窑烟室氧气含量、烟室氮氧化物含量、三次风温和三次风压力。
烧成***回转窑的煅烧工艺具有非线性,强耦合,干扰因素多的特点,这对构建一个合适的f-CaO含量预测模型会产生影响。因此需要初步筛选工艺参数,将对f-CaO含量影响不大的工艺参数进行剔除。步骤S2中利用灰色关联分析法对步骤S1收集到的数据进行初步筛选,具体为:
设自变量为X=(xi(l)|i=1,2,...,h,l=1,2,...,M),因变量为T={t(l)|l=1,2,...,M},其中h表示自变量的数目,M表示样本的数目,令Δi(l)=|t(l)-xi(l)|,则灰色关联度分析法计算得到的灰关联度为:
式中,ρ为分辨系数,0<ρ<1;
因变量和自变量的灰色关联度γi如下:
设置灰色关联度阈值,选择灰色关联度大于灰色关联度阈值的数据作为数据样本,完成初步筛选。
采用灰色关联分析法各个工艺参数与f-CaO含量之间的相关性。得到的各工艺参数与f-CaO含量之间的灰关联度如表1。
通过采集的200组数据,利用灰色关联分析法计算工艺参数与f-CaO含量灰关联值。通过表1可以发现,烧成***回转窑内不同的工艺参数与f-CaO含量的关联程度存在明显的差异。为了更好的建立工艺参数与f-CaO含量之间的函数关系,现设定灰关联度阈值为0.750,则初步选取喂料量,煤粉流量,分解炉出口温度,物料温度,二次风温,回转窑电机电流,回转窑尾温,回转窑烟室氧气含量,烟室氮氧化物含量9个工艺参数作为f-CaO含量预测模型的输入集合,具体表示为:
X'={x1(l),x2(l),...,x9(l)|l=1,2,...,200}
f-CaO含量变量集合表示为下式:
T'={t(l)|l=1,2,...,200}
将X',T'进行归一化处理得到X,T。
表1工艺参数与水泥f-CaO含量的灰关联度
所述步骤S3中初始化分层极限学习机的参数具体为:
随机初始化隐含层个数、神经元个数、隐含层连接权值和阈值。
步骤S3中利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型,具体为:
a、将分层极限学习机的隐层神经元个数组成乌鸦算法的初始种群集合,初始化乌鸦算法的参数;
b、以分层极限学习机的训练结果误差作为适应度函数,评估乌鸦的位置;
c、择优选择乌鸦种群集合中乌鸦的新位置,评估乌鸦的新位置,如果新位置比当前位置适应度值高,更新乌鸦位置至新位置;否则,乌鸦的新位置随机更新;
d、重复步骤c直到迭代次数达到阈值,迭代终止,输出当前最佳位置,得到分层极限学习机的最优隐层神经元个数;
e、分层极限学习机利用FISTA算法计算每一隐含层网络权重,计算该层的实际输出并且作为下一层输入;
f、运用最小二乘法求得分层极限学习机最后一层隐层和输出层之间的权值,输出预测结果。
步骤c中乌鸦的新位置更新,具体如下:
择优选择乌鸦种群集合中乌鸦的新位置:
Xi,iter+1=Xi,iter+ri×fl×(Mj,iter-Xi,iter)
式中,向量Xi,iter表示第i只乌鸦Crowi在第i次迭代中的位置,Xi,iter+1表示第i只乌鸦Crowi在第i+1次迭代中的位置,Mj,iter表示第j只乌鸦Crowi在第i次的记忆,ri为0-1之间的随机数,fl为乌鸦移动的步长;
乌鸦的新位置随机更新:
Xi,iter=a random position
Mi,iter更新为:
式中,f()表示乌鸦算法的适应度函数。
步骤e中所述分层极限学习机,具体为:
所述分层极限学习机是一种具有多层特性的结构,其分层体现在整个模型的训练过程分为无监督特征提取和有监督特征分类。
在无监督特征提取阶段,基于极限学习机自动编码器ELM-AE思想,并使用基于l1-范数范数最小化的快速阈值收缩算法优化隐含层权值,分类部分采用基于岭回归的ELM模型作为分类模型;
在无监督学习阶段,分层极限学习机利用基于极限学习机自动编码器ELM-AE作为特征提取工具,ELM-AE中,输出等于输入,隐藏节点通过随机正交权重和偏差将输入数据投影到不同或相等的维度,实现原始数据在不同维度的特征表达;
基于极限学习机自动编码器ELM-AE训练过程可以用以下优化公式表述为:
式中,H表示隐含层输出矩阵,X表示每一层的输入矩阵,β表示输出权值;
采用l1-范数对β进行约束,此时基于极限学习机自动编码器ELM-AE可以表示为如下约束问题:
为了明确l1优化问题,可以将上式从写成公式:
Oβ=p(β)+q(β)
其中,p(β)=||Hβ-X||2,q(β)=||β||l1,再使用基于l1-范数范数最小化的快速阈值收缩算法,即FISTA算法优化隐含层权值;
经过n层的无监督学习,最终可获得高层稀疏特征,每个隐含层的输出可以表示为:
Hi=g(Hi-1·β)
式中,Hi代表第i隐含层输出矩阵,Hi-1代表第i-1隐含层输出矩阵,g(·)代表隐藏层的激活函数。
所述FISTA算法优化隐含层权值:
F1:计算光滑凸函数P的利普希兹常量L;
F2:初始化参数yl=β0,tl=1,y为隐层输出权重,t为迭代次数,然后对于k(k≥1)进行如下迭代:
βk=pL(yk)
其中pL可以由下式给出:
其中β代表隐含层输出权重,βk-1代表上一次迭代的隐含层输出权重;
F3:通过上述迭代,计算出输出权值β。
在有监督特征分类阶段,利用最小二乘法计算最后一层的隐含层输出矩阵。
本次共采集历史数据200组,随机选取150组作为CSA-H-ELM模型的训练样本,输入集合为X1,输出集合为T1。剩下50组数据作为CSA-H-ELM模型的测试样本,输入集合为X2,输出集合为T2。通过训练样本训练CSA-H-ELM模型,完善模型的网络结构以及参数,建立f-CaO含量预测模型。其中H-ELM模型共有9个输入变量,1个输出变量,隐含层层数经过测试选取3层,隐含层神经元个数由CSA模型确定;CSA模型共有3个输入变量,3个输出变量,设置种群规模为30,步长为fl=2,感知概率AP=0.5,搜索空间Q为[10,400].
通过对本实施例提出的CSA-H-ELM模型进行训练,训练值与真实值对比结果如图2所示。
通过对训练完成的Crows-H-ELM模型进行验证测试,测试值与真实值对比结果如图3所示。
此外,本实施例还将该方法与其他预测方法进行了比较,包括ELM模型,基于GA的BPNN模型(GA-BPNN),基于GWO的SVM模型(GWO-SVM),基于QPSO的SVM模型(QPSO-SVM),H-ELM模型,条件相同。表2显示了RMSE和消耗时间性能比较。表2显示,本发明所提出的CSA-H-ELM算法无论是在训练过程,还是测试过程,RMSE均小于上述5个方法,这是因为H-ELM采用-范数的快速迭代阈值法可以实现特征的自动提取并且映射在高维空间,实现特征的稀疏表示,大大挖掘了输入数据中的有用信息。对比表2中H-ELM和所提出的CSA-H-ELM算法可以发现,CSA-H-ELM算法在训练测试时间上略长于H-ELM,这是因为使用改进的乌鸦算法搜寻最为合适的隐含层神经元个数需要花费额外的时间,而CSA-H-ELM算法的RMSE更小,说明所提出的改进乌鸦算法确实能有效的降低预测误差,提高预测精度,证明了所提出算法的有效性。
表2
相同或相似的标号对应相同或相似的部件;
附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制;
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.一种集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,包括以下步骤:
S1:采集与水泥生产质量预测相关的数据;
S2:使用灰色关联分析法对步骤S1收集到的数据进行初步筛选;
S3:初始化分层极限学习机的参数,利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型;
S4:利用所述预测模型对水泥生产质量进行预测。
2.根据权利要求1所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,步骤S1中所述与水泥生产质量预测相关的数据包括与烧成***回转窑相关的喂料量、煤粉流量、分解炉出口温度、物料温度、二次风温、压力、回转窑电机电流、转速、高温风机转速、回转窑尾温、回转窑头负压、回转窑尾压力、回转窑烟室氧气含量、烟室氮氧化物含量、三次风温和三次风压力。
4.根据权利要求1所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,所述步骤S3中初始化分层极限学习机的参数具体为:
随机初始化隐含层个数、神经元个数、隐含层连接权值和阈值。
5.根据权利要求4所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,步骤S3中利用步骤S2中筛选得到的数据对分层极限学习机进行训练时,利用乌鸦算法优化分层极限学习机的隐层神经元个数,获取分层极限学习机的最优隐层神经元个数,得到预测模型,具体为:
a、将分层极限学习机的隐层神经元个数组成乌鸦算法的初始种群集合,初始化乌鸦算法的参数;
b、以分层极限学习机的训练结果误差作为适应度函数,评估乌鸦的位置;
c、择优选择乌鸦种群集合中乌鸦的新位置,评估乌鸦的新位置,如果新位置比当前位置适应度值高,更新乌鸦位置至新位置;否则,乌鸦的新位置随机更新;
d、重复步骤c直到迭代次数达到阈值,迭代终止,输出当前最佳位置,得到分层极限学习机的最优隐层神经元个数;
e、分层极限学习机利用FISTA算法计算每一隐含层网络权重,计算该层的实际输出并且作为下一层输入;
f、运用最小二乘法求得分层极限学习机最后一层隐层和输出层之间的权值,输出预测结果。
6.根据权利要求5所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,步骤c中乌鸦的新位置更新,具体如下:
择优选择乌鸦种群集合中乌鸦的新位置:
Xi,iter+1=Xi,iter+ri×fl×(Mj,iter-Xi,iter)
式中,向量Xi,iter表示第i只乌鸦Crowi在第i次迭代中的位置,Xi,iter+1表示第i只乌鸦Crowi在第i+1次迭代中的位置,Mj,iter表示第j只乌鸦Crowi在第i次的记忆,ri为0-1之间的随机数,fl为乌鸦移动的步长;
乌鸦的新位置随机更新:
Xi,iter=a random position
Mi,iter更新为:
式中,f()表示乌鸦算法的适应度函数。
7.根据权利要求6所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,步骤e中所述分层极限学习机,具体为:
所述分层极限学习机是一种具有多层特性的结构,其分层体现在整个模型的训练过程分为无监督特征提取和有监督特征分类。
8.根据权利要求7所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,在无监督特征提取阶段,基于极限学习机自动编码器ELM-AE思想,并使用基于l1-范数范数最小化的快速阈值收缩算法优化隐含层权值,分类部分采用基于岭回归的ELM模型作为分类模型;
在无监督学习阶段,分层极限学习机利用基于极限学习机自动编码器ELM-AE作为特征提取工具,ELM-AE中,输出等于输入,隐藏节点通过随机正交权重和偏差将输入数据投影到不同或相等的维度,实现原始数据在不同维度的特征表达;
基于极限学习机自动编码器ELM-AE训练过程可以用以下优化公式表述为:
式中,H表示隐含层输出矩阵,X表示每一层的输入矩阵,β表示输出权值;
采用l1-范数对β进行约束,此时基于极限学习机自动编码器ELM-AE可以表示为如下约束问题:
为了明确l1优化问题,可以将上式从写成公式:
Oβ=p(β)+q(β)
其中,p(β)=||Hβ-X||2,q(β)=||β||l1,再使用基于l1-范数范数最小化的快速阈值收缩算法,即FISTA算法优化隐含层权值;
经过n层的无监督学习,最终可获得高层稀疏特征,每个隐含层的输出可以表示为:
Hi=g(Hi-1·β)
式中,Hi代表第i隐含层输出矩阵,Hi-1代表第i-1隐含层输出矩阵,g(·)代表隐藏层的激活函数。
10.根据权利要求9所述的集成CSA和H-ELM的水泥生产质量预测方法,其特征在于,在有监督特征分类阶段,利用最小二乘法计算最后一层的隐含层输出矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110721214.8A CN113537579A (zh) | 2021-06-28 | 2021-06-28 | 一种集成csa和h-elm的水泥生产质量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110721214.8A CN113537579A (zh) | 2021-06-28 | 2021-06-28 | 一种集成csa和h-elm的水泥生产质量预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113537579A true CN113537579A (zh) | 2021-10-22 |
Family
ID=78126054
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110721214.8A Pending CN113537579A (zh) | 2021-06-28 | 2021-06-28 | 一种集成csa和h-elm的水泥生产质量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113537579A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117743620A (zh) * | 2024-02-21 | 2024-03-22 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种大岩土数智*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202946A (zh) * | 2016-07-18 | 2016-12-07 | 燕山大学 | 基于深度信念网络模型的水泥熟料游离钙含量预测方法 |
CN109933881A (zh) * | 2019-03-06 | 2019-06-25 | 武汉大学 | 一种基于优化深度信念网络的电力电子电路故障诊断方法 |
CN112329370A (zh) * | 2020-11-02 | 2021-02-05 | 南京航空航天大学 | 一种基于极限学习机的SiC功率管在线故障诊断方法 |
-
2021
- 2021-06-28 CN CN202110721214.8A patent/CN113537579A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202946A (zh) * | 2016-07-18 | 2016-12-07 | 燕山大学 | 基于深度信念网络模型的水泥熟料游离钙含量预测方法 |
CN109933881A (zh) * | 2019-03-06 | 2019-06-25 | 武汉大学 | 一种基于优化深度信念网络的电力电子电路故障诊断方法 |
CN112329370A (zh) * | 2020-11-02 | 2021-02-05 | 南京航空航天大学 | 一种基于极限学习机的SiC功率管在线故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
FANG DU等: ""An Effective Hierarchical Extreme Learning Machine based Multimodal Fusion Framework"", 《NEUROCOMPUTING》, vol. 322, pages 1 - 21 * |
史永胜等: ""基于分层极限学习机的锂离子电池剩余使用寿命预测"", 《电子器件》, vol. 44, no. 3, pages 652 - 667 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117743620A (zh) * | 2024-02-21 | 2024-03-22 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种大岩土数智*** |
CN117743620B (zh) * | 2024-02-21 | 2024-05-03 | 中国有色金属工业昆明勘察设计研究院有限公司 | 一种大岩土数智*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109508818B (zh) | 一种基于LSSVM的在线NOx预测方法 | |
CN113095550B (zh) | 基于变分递归网络和自注意力机制的空气质量预测方法 | |
CN106202946A (zh) | 基于深度信念网络模型的水泥熟料游离钙含量预测方法 | |
CN110610061B (zh) | 一种融合多源信息的混凝土塌落度高精度预测方法 | |
CN113177358B (zh) | 一种基于模糊细粒度特征提取的水泥质量软测量方法 | |
CN111079978B (zh) | 基于逻辑回归和增强学习的煤与瓦斯突出预测的方法 | |
CN113589693B (zh) | 基于邻域优化的水泥工业分解炉温度模型预测控制方法 | |
CN114721263B (zh) | 基于机器学习和智能优化算法的水泥分解炉智慧调控方法 | |
CN106327004A (zh) | 一种基于熟料质量指标的水泥烧成过程优化方法 | |
CN113537579A (zh) | 一种集成csa和h-elm的水泥生产质量预测方法 | |
CN117289652A (zh) | 一种基于多元宇宙优化的数控机床主轴热误差建模方法 | |
CN114239400A (zh) | 基于局部双加权概率隐变量回归模型的多工况过程自适应软测量建模方法 | |
CN109992844A (zh) | 一种基于adqpso-svr模型的锅炉飞灰含碳量预测方法 | |
CN114334025A (zh) | 一种水泥熟料煅烧环境变量的构建及验证方法 | |
CN112836902B (zh) | 一种基于改进型bp神经网络的煤燃烧发热量预测方法 | |
CN107844679B (zh) | 一种分段式高炉铁水硅含量预测方法及装置 | |
CN112598050A (zh) | 生态环境数据质量控制方法 | |
CN109101683B (zh) | 燃煤分质利用与清洁预处理***热解釜的模型更新方法 | |
CN112541296A (zh) | 一种基于pso-lssvm的so2预测方法 | |
CN111639821A (zh) | 一种水泥炉窑生产能耗预测方法及*** | |
CN117290752A (zh) | 一种磨煤机煤粉细度的预测方法 | |
CN116631534A (zh) | 一种基于数据增强变压器网络的水泥熟料游离钙预测方法 | |
CN110648023A (zh) | 基于二次指数平滑改进gm(1,1)的数据预测模型的建立方法 | |
CN115271189A (zh) | 一种数据驱动的水泥熟料游离氧化钙的预测方法 | |
CN115392125A (zh) | 水泥回转窑的温度预测方法 |
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 |