[发明专利]循环流化床生活垃圾焚烧锅炉炉膛出口烟气含氧量的预测系统及方法

摘要

本发明公开了一种循环流化床生活垃圾焚烧锅炉炉膛出口烟气含氧量的实时预测系统及方法，采用支持向量机算法和多种群遗传粒子群优化算法集成建模的方法，构建了一种快速经济且自适应更新的系统和方法对锅炉炉膛出口烟气含氧量进行实时预测，避开了繁琐复杂的机理建模工作。利用SVM算法的非线性动力学特性、泛化能力和实时预测能力来表征烟气含氧量的动态变化特性；利用粒子群优化算法对SVM算法惩罚参数C和核参数g进行优化，提高模型的泛化能力；引入遗传操作算子和多种群迁移机制，加快粒子群算法的收敛速度，提高粒子群优化算法解的多样性，减少粒子群算法寻优算陷入局部最优的可能性，提高算法的全局搜索能力和局部搜索能力。

主权项

一种循环流化床生活垃圾焚烧锅炉炉膛出口烟气含氧量的实时预测系统，其特征在于，该系统与循环流化床锅炉的集散控制系统以及生产管理系统相连，包括数据通讯接口和上位机，在上位机中对烟气含氧量预测模型进行训练和更新，然后将训练好的模型通过通讯接口送往集散控制系统、生产管理系统，所述上位机包括：信号采集模块。该模块用于采集CFB生活垃圾焚烧锅炉在焚烧指定生活垃圾时的运行工况状态参数和操作变量，并组成垃圾热值预测模型输入变量的训练样本矩阵X(m×n)，m表示样本个数，n表示变量的个数。数据预处理模块。对X(m×n)进行粗大误差处理和随机误差处理，剔除训练样本中的野值，排除锅炉停炉、压火、给料机堵塞等异常工况，训练样本输入变量经过归一化处理后映射到[0，1]区间内，得到标准化后的训练样本X*(m×n)。专家知识库模块。采用滚动时间窗口的方法不断更新训练样本，使得训练样本始终保持在最新的状态上，滚动时间窗口方法是指从当前时间开始，回溯L(单位秒)长度的时间尺寸。模型更新判定模块。检测当前预测模型的性能，当相对预测误差超过±5％时，则判定模型需要进行更新。智能建模模块。该模块先利用多种群遗传粒子群算法对SVM模型的惩罚参数C和核参数g进行寻优，然后将得到的最优参数组合赋给SVM模型，并以此为基础进行训练。具体步骤如下：2.1)初始化算法参数。多种群粒子群算法的最大寻优代数Tmax、最大惯性权重ωmax、最小惯性权重ωmin，速度更新系数R1、R2、R3、R4，种群数量pop、单个种群的粒子数量ind、移民率Pi、交叉概率Pc、变异概率Pm、遗传代沟Pe、遗传操作频度Pg。2.2)初始化种群。采用实数编码的方式，将SVM模型的惩罚参数C和核参数g有序编码在一个粒子当中，并将每个权值阈值随机生成为[0，200]之间的一个实数。2.3)初始化个体极值和群体极值。将每个粒子中包含的参数组合赋给SVM模型，并结合训练样本进行学习训练，利用训练得到的预测模型计算烟气含氧量预测值将预测值与实际测量值y*进行比较，并以误差平方和MSE作为粒子的适应度值fitness，适应度计算公式如下：fitness=MSE=Σi=1m(y^-y*)2---(2)]]>每个粒子的计算得到的适应度值作为该粒子本身的初始化极值，每个种群的当中MSE最低的值作为群体极值。2.4)更新粒子。根据最新的个体极值和全局极值，按照(3)式和(4)式更新粒子的速度vid(t)和位置xid(t)：vid(t+1)＝ωvid(t)+c1r1(pid‑xid(t))+c2r2(pgd‑xid(t))  (3)xid(t+1)＝xid(t)+vid(t+1)  (4)式中，t是粒子群优化算法的寻优代数，r1，r2是[0，1]之间的随机数，pid是指是指第i个粒子迄今为止寻搜索到最优位置，pgd指该种群迄今为止搜索最优位置。更进一步，为了改善基本粒子群算法容易陷入局部极值和收敛速度慢的缺陷，在PSO算法的基础上引进了动态加速常数c1、c2和惯性权重ω：ω=ωmax-ωmax-ωminTmaxt---(5)]]>c1=R1+R2×tTmax---(6)]]>c2=R3-R4×tTmax---(7)]]>其中，Tmax为最大寻优代数，ωmax为最大惯性权重，ωmin为最小惯性权重，R1、R2、R3、R4为常数。2.5)判断是否需要进行遗传操作。按照遗传操作频度Pg进行，通常每Pg代进行一次遗传操作。若要进行遗传操作则执行步骤2.6)，否则执行步骤2.7)2.6)执行遗传操作。遗传操作包括选择、交叉、变异、留优，具体如下。选择：像轮盘赌一样计算选择概率，等距离的选择个体，设N为需要选择的个体数目，选择指针的距离是1/N，第一个选择指针的位置由[0，1/N]区间内的随机数决定；交叉：采用算术重组的方法，即两个个体通过线性组合产生新的个体。假设在两个个体和之间进行算术交叉，则按式(8)和(9)得到两个新的个体和xA1t=λ1xA0t+λ2xB0t---(8)]]>xB1t=λ2xA0t+λ1xB0t---(9)]]>式中λ1和λ2为(0,1)之间的实数，且λ1+λ2＝1；变异：采取均匀性变异，假设s＝(v1,v2,...,vn)是父代，z＝(z1,z2,...,zn)是变异产生的后代。现在父代向量中随机地选择一个分量，假设是第k个，然后在其定义区间[ak,bk]中均匀随机地取一个数v'k代替vk以得到zk；留优：采用保留父代精英个体的策略，根据父代染色体的目标函数值，用过渡代染色体替换父代中目标函数值排名后Pe×100％的染色体，Pe为替换率(遗传代沟)，如此，父代中部分最优的个体得以保留，不会因为选择、交叉、变异等运算消失。2.7)粒子适应度值计算。按照式(2)计算更新之后粒子的适应度值。2.8)更新个体极值和群体极值。以适应度值为评价指标，比较当代粒子与上一代粒子之间的适应度值大小，如果当前粒子的适应度值优于上一代，则将当前粒子的位置设置为个体极值，否则个体极值保持不变。同时获取当代所有粒子适应度值最优的粒子，并与上一代最优粒子进行比较，如果当代最优粒子的适应度值优于上一代最优粒子的适应度值，则将当代粒子的最优适应度值设置为全局最优值，否则全局最优值保持不变。2.9)移民操作。采用多种群遗传粒子群算法，每一次寻优过程中，种群之间采用单向循环迁移的方式进行移民操作，第1个种群中的优秀个体迁往第2个，第2个迁往第3个，以此类推，直到最后一个迁往第一个。种群间的移民率Pi＝0.04，即表示用源种群中排名前Pi×100％的个体替换目标种群当中排名后Pi×100％的个体，以此完成种群之间最优知识的交流。2.10)算法停止条件算法判定。判断是否达到最大迭代次数或者到达预测精度的要求，如果没有达到则返回步骤2.4)，利用更新的粒子继续搜索，否则退出搜索，执行步骤2.11)。2.11)输出最优的参数组合粒子。2.12)将最优粒子当中的参数组合赋给SVM模型，并结合训练样本进行学习。2.13)验证模型的预测精度。将模型的预测值和实际值进行对比，计算相对预测误差。2.14)判断相对预测误差是否在±5％以内，如果满足要求则执行步骤2.15)，否则返回步骤2.12)，重新进行训练。2.15)输出满足要求的烟气含氧量预测模型。通讯模块。该模块将满足要求烟气含氧量预测模型传送给函数集散控制系统、生产管理系统。