一种融合最小二乘支持向量机回归与粒子群优化的挥发窑工况操作模式优化方法与流程

本发明涉及一种挥发窑工况操作模式优化方法，具体的说就是一种融合最小二乘支持向量机回归(lssvm)与粒子群优化的挥发窑工况操作模式优化方法。
背景技术：
挥发窑广泛用于冶金、建材、化工和造纸等领域，挥发窑处理物料的过程是一个典型的热工反应过程，由于窑的主要部分在回转，检测仪器通常安装在窑的两端，导致关键工艺参数不能准确及时测量，窑内热工状况难以定量描述。目前挥发窑生产过程多依赖操作人员经验根据窑头火焰的燃烧情况手动调整窑转速、鼓风量等操作参数，这使得挥发窑生产过程技术经济指标波动较大，运行状况不稳定，单位产品的能耗、原材料消耗和有害气体排放较高，严重制约企业节能减排水平的提高。目前挥发窑工况识别常通过表征窑特征的模型实现，模型分为机理模型与数据驱动模型两类。由于挥发窑生产过程反应机理复杂，很难建立精确的机理模型实现操作参数的在线识别，基于机理模型的识别方法一般需要做一些假设简化建模过程，机理模型简化带入的误差以及生产扰动带来的不确定性导致模型的精确性不高，因此传统的机理模型很少有直接应用于挥发窑实际生产过程的报道。基于数据驱动模型的识别方法能依据大量的实际生产数据实现挥发窑操作参数的优化，但是神经网络存在结构(神经网络类型、隐层层数等)、隐层节点数难以确定等问题，降低了它的应用效果。同样对于利用模糊模型识别操作参数而言，由于挥发窑的生产条件一直在缓慢变化，操作规律也随之变化，故如何获取模糊规则成为解决问题的瓶颈，常见的基于专家知识的方法受限于领域专家的主观意志、知识水平和操作经验，很难及时更新，可能导致规则冗余或不完整。最小二乘支持向量机(leastsquaresupportvectormachine，lssvm)是一种新兴的统计学方法，目前在挥发窑生产中lssvm回归理论还应用得很少，将lssvm方法用于挥发窑工况操作模式识别和性能预测，有利于克服上述智能算法的不足。在挥发窑工况操作模式优化过程中，通常按产品质量、能耗、有害气体排放工艺指标进行评价，这实际上是一个多目标优化决策问题，目前在挥发窑的文献中还很少有这方面的报道。技术实现要素：技术问题：利用基于最小二乘支持向量机回归模型(lssvm)预测挥发窑工艺指标，并利用粒子群优化算法寻找优化的挥发窑工况操作模式，使挥发窑长期运行在优化的操作模式区域。技术方案：挥发窑工况操作模式定义为：由一定的输入条件f，与之对应的操作参数o，状态参数j和性能工艺指标i组成的向量[foji]，其中参数可以是连续或离散的。挥发窑工况操作模式优化定义为：根据挥发窑的输入条件以及当前状态,寻找合适的工况操作参数,使得挥发窑的产品质量、能耗、有害气体排放工艺指标达到最优。本发明提出的挥发窑工况操作模式优化方法分为2部分：(1)建立基于最小二乘支持向量机的挥发窑工况评价模型如附图1所示建立基于最小二乘支持向量机的挥发窑工况评价模型,设(xi，yi)＝[foji]i且{(xi，yi)|xi∈rn，yi∈r，i＝1，2，…n}，根据lssvm理论，利用结构风险最小化原则，函数拟合问题可描述为最优化问题：式中是将输入数据映射到高维特征空间的函数；ω∈rn，ei为误差，ei，b∈r，c＞0为惩罚系数，用于控制解的光滑度，其值越大代表对误差的惩罚力度越强，t为倒置。根据式(1)将模型转换到对偶空间加以解决，得到如下lagrange函数：式中αi∈r是lagrange乘子，分别对ω，ei，b，αi求偏导，并令偏导数为0：消去ω，ei整理得线性方程组(4)：式中y＝[y1…yn]，α＝[α1…αn]，et＝[1…1]，k为满足mercer定理的核函数：求解式(4)得α，b，且获得的非线性回归函数为：核函数k(xi，x)的选取支持向量机是一种基于核的学习方法，核函数选取对支持向量机性能有着重要的影响。为此，本发明分别利用表一所示的3种核函数建立lssvm模型评判挥发窑工况，并运用k层交叉验证(kfoldcross-validation)方法对模型预测性能评估，选取预测误差(rmse)最小的lssvm模型评判挥发窑工况。表一采用的3种核函数核函数的名称核函数的表达式多项式核函数kp(xi，x)＝(xtxi+1)u(u∈n)rbf核函数kr(xi，x)＝exp(-‖x-xi‖2/2σ2)(σ≠0∈r)sigmoid核函数ks(xi，x)＝tanh(a(xtxi)+c)(a，c∈r)表一中kp，kr，ks分别表示多项式核函数，rbf核函数，sigmoid核函数。多项式核函数中u用来设置多项式核函数的最高项次数；rbf核函数中σ为函数的宽度参数，控制函数的径向作用范围；sigmoid核函数中a，c用来设置核函数中的参数。(2)基于粒子群算法的挥发窑工况操作模式优化挥发窑的工艺指标定义为：i＝μ1(1-z/zmax)2+μ2(1-m/mmax)2+μ3(1-q/qmax)2(6)其中z，m，q分别为窑渣含锌量，窑尾烟气温度，废气含尘浓度预测值，zmax，mmax，qmax分别为窑渣含锌量，窑尾烟气温度，废气含尘浓度最大值，μ1，μ2，μ3为权重，且0＜z≤zmaxmmin≤m≤mmax0＜q≤qmaxstep1：初始化，设定粒子群算法寻优参数的取值范围包括学习因子c1，c2；最大进化代数tmax，在d维空间中随机产生m个粒子，第d个粒子(d＝1，2，…m)的位置表示为qdh＝(qd1，qd2…qdh)，飞行速度表示为vdh＝(vd1，vd2…vdh)，(h＝1，2，…d)，设定飞行速度的上下限值，根据工艺要求输入影响工艺指标值idh变化的参数取值范围包括锌渣含量z，窑尾烟气温度m，废气含尘浓度q取值范围以及权重μ1，μ2，μ3取值；step2：计算每个粒子在当前位置处的适应度值以及对应最优适应度值的最好位置记为qbestdh，且qbestdh＝fitness(idh)，在这里本发明根据第一部分步骤，建立基于最小二乘支持向量机的挥发窑工况评价模型，预测挥发窑的工艺指标idh，从而获得适应度值fitness(idh)，在整个群体中，所有粒子经历过的最优适应度值及对应的最好位置记为gbestdh，且gbestdh＝max(fitness(i11)，fitness(i21)…fitness(idh))；step3：在每次迭代过程中，每个粒子根据下面公式更新自己的速度和位置：vdh(t+1)＝vdh(t)+c1r1(qbestdh-qdh(t))+c2r2(gbestdh-qdh(t))qdh(t+1)＝qdh(t)+vdh(t+1)其中c1，c2为学习因子；r1，r2是[0，1]之间的随机数；t为当前进化代数。更新qbestdh和gbestdh，如果粒子当前的适应值qbestdh更优，则置qbestdh为当前值，并设qdh位置为d维空间中的当前位置；如果粒子当前的适应值gbestdh更优，则置gbestdh为当前粒子的适应值，并设gbestdh位置为d维空间中的当前位置；step4：判断寻优是否达到最大进化代数tmax，或|vdh(t+1)-vdh(t)|＜ε，ε是粒子相邻进化代数之间速度变化的阈值，且i0为临界性能指标值，如果是的则输出当前最优挥发窑操作模式即当前窑身转速，鼓风风压，鼓风风量，窑内负压值，否则t＝t+1转step2。有益效果：本发明所提挥发窑工况操作模式优化方法可以优化挥发窑工况操作模式，使系统长期处于较优操作模式下，同时可以改变挥发窑生产过程依赖操作经验作为主要控制手段的状况,有效抑制挥发窑生产过程技术经济指标波动较大，运行状况不稳定现象，有利于企业节能减排水平的提高。附图说明图1为基于最小二乘支持向量机回归的挥发窑工况评价模型结构图。图2为融合最小二乘支持向量机回归与粒子群优化的挥发窑工况操作模式优化方法原理图。具体实施方式：本发明提出的融合最小二乘支持向量机回归与粒子群优化的挥发窑工况操作模式优化方法具体实施方案详述如下:本发明提出的挥发窑工况操作模式优化方法分为2部分：(1)建立基于最小二乘支持向量机的挥发窑工况评价模型如附图1所示建立基于最小二乘支持向量机的挥发窑工况评价模型，影响挥发窑工艺指标预测的参数有窑身转速，鼓风风压，鼓风风量，窑内负压，入窑物料配比，进料量，煤焦比，将这7个参数作为挥发窑工况评价模型的输入，预测窑渣含锌量，窑尾烟气温度，废气含尘浓度，带入公式(6)计算工艺指标值作为输出。针对400个样本数据，其中250个作为训练数据用来建立基于最小二乘支持向量机的挥发窑工况评价模型，100个作为测试数据，50个作为校验数据用来测试模型的泛化能力，计算预测误差。设(xi，yi)＝[foji]i且{(xi，yi)|xi∈rn，yi∈r，i＝1，2，…n}，根据lssvm理论，利用结构风险最小化原则，函数拟合问题可描述为最优化问题：式中是将输入数据映射到高维特征空间的函数；ω∈rn，ei为误差，ei，b∈r，c＞0为惩罚系数，用于控制解的光滑度，c＝1，t为倒置。根据式(1)将模型转换到对偶空间加以解决，得到如下lagrange函数：式中αi∈r是lagrange乘子，分别对ω，ei，b，αi求偏导，并令偏导数为0：消去ω，ei整理得线性方程组(4)：式中y＝[y1…yn]，α＝[α1…αn]，et＝[1…1]，k为满足mercer定理的核函数：求解式(4)得α，b，且获得的非线性回归函数为：核函数k(xi，x)的选取支持向量机是一种基于核的学习方法，核函数选取对支持向量机性能有着重要的影响。为此，本发明分别利用表一所示的3种核函数建立lssvm模型评判挥发窑工况，并运用10层交叉验证方法对模型预测性能评估，选取预测误差(rmse)最小的lssvm模型评判挥发窑工况。表一采用的3种核函数核函数的名称核函数的表达式多项式核函数kp(xi，x)＝(xtxi+1)u(u∈n)rbf核函数kr(xi，x)＝exp(-‖x-xi‖2/2σ2)(σ≠0∈r)sigmoid核函数ks(xi，x)＝tanh(a(xtxi)+c)(a，c∈r)表一中kp，kr，ks分别表示多项式核函数，rbf核函数，sigmoid核函数。多项式核函数中u用来设置多项式核函数的最高项次数，u＝3；rbf核函数中σ为函数的宽度参数，控制函数的径向作用范围，σ＝1；sigmoid核函数中a，c用来设置核函数中的参数，a＝c＝1。(2)基于粒子群算法的挥发窑工况操作模式优化挥发窑的工艺指标定义为：i＝μ1(1-z/zmax)2+μ2(1-m/mmax)2+μ3(1-q/qmax)2(6)其中z，m，q分别为窑渣含锌量，窑尾烟气温度，废气含尘浓度预测值，zmax，mmax，qmax分别为窑渣含锌量，窑尾烟气温度，废气含尘浓度最大值，μ1，μ2，μ3为权重，且0＜z≤zmaxmmin≤m≤mmax0＜q≤qmaxstep1：初始化，设定粒子群算法寻优参数的取值范围包括学习因子c1＝c2＝2，最大进化代数tmax＝1000，将当前进化代数置为t＝1，在d维空间，且d＝3中随机产生m个粒子，且m＝40，第d个粒子(d＝1，2，…40)的位置表示为qdh＝(qd1，qd2…qdh)，飞行速度表示为vdh＝(vd1，vd2…vdh)，(h＝1，2，3)，在整个搜索空间上随机初始化整个粒子群的位置和速度；设定飞行速度的上下限值，0.5≤v≤2.15，根据工艺要求输入0＜z≤3％，400℃≤m≤800℃，0＜q≤30％，μ1＝0.4，μ2＝μ3＝0.3；step2：计算每个粒子在当前位置处的适应度值以及对应最优适应度值的最好位置记为qbestdh，且qbestdh＝fitness(idh)，在这里本发明根据第一部分步骤，建立基于最小二乘支持向量机的挥发窑工况评价模型，预测挥发窑的工艺指标idh，从而获得适应度值fitness(idh)，在整个群体中，所有粒子经历过的最优适应度值及对应的最好位置记为gbestdh，且gbestdh＝max(fitness(i11)，fitness(i21)…fitness(idh))；step3：step3：在每次迭代过程中，每个粒子根据下面公式更新自己的速度和位置：vdh(t+1)＝vdh(t)+c1r1(qbestdh-qdh(t))+c2r2(gbestdh-qdh(t))qdh(t+1)＝qdh(t)+vdh(t+1)其中r1，r2是[0，1]之间的随机数，且r1，＝r2＝0.5。更新qbestdh和gbestdh，如果粒子当前的适应值qbestdh更优，则置qbestdh为当前值，并设qdh位置为d维空间中的当前位置；如果粒子当前的适应值gbestdh更优，则置gbestdh为当前粒子的适应值，并设gbestdh位置为d维空间中的当前位置；step4：判断寻优是否达到最大进化代数tmax，或|vdh(t+1)-vdh(t)|＜ε，ε是粒子相邻进化代数之间速度变化的阈值，且i0为临界性能指标值，i0＝60％，如果是的则输出当前最优挥发窑操作模式即当前窑身转速，鼓风风压，鼓风风量，窑内负压值，否则t＝t+1转step2。当前第1页12