一种高炉炼铁过程铁水硅含量动态软测量系统及方法
【专利摘要】本发明提供一种高炉炼铁过程铁水硅含量动态软测量系统及方法。系统包括：实际数据采集单元、归一化预处理单元、动态软测量单元；方法包括：获取高炉铁水硅含量动态软测量所需参数：当前炉次的工况参数、上一炉次的工况参数和上一炉次的铁水硅含量；对获取的高炉铁水硅含量动态软测量所需参数进行归一化预处理；采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量动态软测量。本发明与现有的人工测量或化验铁水硅含量相比，减少了操作人员的工作量，降低了人工操作引入的测量的不确定性，提高了测量的时效性和准确性，可信度高。本发明方法对高炉冶炼过程的铁水硅含量预报具有通用性，有助于实现高炉铁水质量闭环一体化控制与优化运行。
【专利说明】
一种高炉炼铁过程铁水硅含量动态软测量系统及方法
技术领域
[0001] 本发明涉及高炉冶炼自动化控制技术领域，特别涉及一种高炉炼铁过程铁水硅含 量动态软测量系统及方法。
【背景技术】
[0002] 高炉生产是在高温、高压、多物理场共存、化学反应与传递现象同时发生的密闭条 件下进行的持续动态时变非线性的多变量耦合系统，保持合理的炉温是维持高炉稳定顺 行、高产、优质、低耗的关键因素。在冶炼过程中，炉温控制的好坏直接影响炉况的波动，而 炉况的波动又决定炉温的控制模式，炉温"过高"或"过低"都容易诱发炉况异常。实际生产 中，实现自动化水平较高的全流程一体化的闭环控制与操作优化的关键技术在于精确合理 的高炉炉温控制模型和时效的在线检测。作为最复杂的逆流反应器之一的炼铁高炉，受当 前现有传感器测量上的影响，一般通过高炉铁水硅含量(即化学热)来间接地反映炉内温度 的变化和高炉炉缸的热状态。高炉铁水的硅含量成为反映炉内物理化学反应情况、热状态 和铁水质量的重要指标，其变化规律直接反映和决定了高炉是否能够稳定顺行，也是表征 高炉热状态及变化的重要标志。为获取高炉内部热状态，进而有效地控制高炉炉况的稳定 顺行，建立高炉铁水硅含量预报模型非常必要。受高炉运行参数严重耦合、非线性较强、对 测量设备有较高的要求等一系列原因的影响导致高炉铁水硅含量的机理模型难以建立。为 此，基于丰富的高炉运行数据，借助于当前的智能建模算法与统计学习理论，建立以数据驱 动的铁水硅含量软测量模型应运而生。
[0003] 专利公开号CN104899425A公开了"一种高炉铁水硅含量的变量选择预报方法"，以 炉顶压力、炉顶温度、料速、炉顶煤气中的co、co 2和上一炉硅含量作为硅预报的输入变量， 基于高斯核函数的标准支持向量机(SVR)预报下一次的含娃量。
[0004] 专利公开号CN101211383A公开了"一种高炉铁水硅含量的特征分析预报方法"，以 高炉铁水硅含量预报模型的高炉工艺参数为输入变量，采用动态独立成分分析方法对输入 变量进行特征提取，使用最小二乘支持向量机算法建立硅含量预报的动态递推模型，引入 遗传算法优化模型参数。
[0005] 专利公开号CN102031319A公开了"一种高炉铁水含硅量的预报方法"，使用移动平 均法对高炉铁水中的硅(Si)含量进行预报。该方法采用铁水含硅量短期、中期、长期的均值 作为预报的模型的基本数据，然后结合神经网络对模型进行训练，得到模型中连接权值。 [0006]上述专利的方法以及其相关的文献大多利用了高炉采集的全部相关变量作为输 入变量，充分利用丰富数据特征的同时也引入了较多的噪声和运行耗时大的问题。另外，在 实际生产过程中，受检测仪表和变送器等装置的故障以及复杂的电磁干扰的影响，测量数 据存在众多的未知干扰，其实际的工业背景对算法本身的鲁棒性有较高的要求，而上述专 利报道的方法未考虑过鲁棒性的问题，针对冶炼工况的变化以及异常的抖动，铁水硅含量 预测模型的泛化能力会大幅下降。此外，上述方法没有考虑输入输出时序及过程时滞关系， 因而不能很好地捕获冶炼过程的固有动态特性。综上所述，建立具有稀疏鲁棒性的高炉铁 水硅含量动态软测量模型非常必要。
【发明内容】
[0007] 针对现有技术的不足，本发明提供一种高炉炼铁过程铁水硅含量动态软测量系统 及方法。
[0008] 本发明的技术方案如下：
[0009] -种高炉炼铁过程铁水硅含量动态软测量系统，包括：
[0010]实际数据采集单元:获取高炉铁水硅含量动态软测量所需参数，包括：当前炉次的 工况参数、上一炉次的工况参数和上一炉次的铁水硅含量；
[0011] 归一化预处理单元:对获取的高炉铁水娃含量动态软测量所需参数进行归一化预 处理；
[0012] 动态软测量单元:采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量 动态软测量:高炉炼铁过程铁水硅含量动态软测量模型采用输出自反馈，以高炉铁水硅含 量动态软测量所需参数为输入，以当前炉次硅含量为输出值，动态预测铁水硅含量。
[0013] 所述实际数据采集单元，包括：
[0014] 辅助变量确定模块:利用主成分分析法对若干炉次的高炉过程工况参数进行相关 性分析，确定对于高炉铁水硅含量动态软测量贡献率最大的工况参数作为辅助变量，具体 包括:炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、喷煤量；
[0015] 输入变量确定模块:基于辅助变量、结合不同时刻输入参数时滞、当前炉次辅助变 量的测量值、上一炉次辅助变量的测量值及上一炉次铁水硅含量的测量值，确定高炉铁水 硅含量动态软测量所需参数即高炉炼铁过程铁水硅含量动态软测量模型的输入变量，包 括:当前炉次的工况参数、上一炉次的工况参数和上一炉次的铁水硅含量;高炉铁水硅含量 动态软测量所需的工况参数包括炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、设定 喷煤量；
[0016] 输入变量采集模块:采集高炉炼铁过程铁水硅含量动态软测量模型的输入变量。 [0017] 所述动态软测量单元，包括：
[0018] 样本处理模块:采集历史若干炉次的高炉铁水硅含量动态软测量所需参数及相应 炉次铁水硅含量，并进行归一化预处理后作为训练数据集；
[0019] ARMA模型构建模块:将训练数据集中的当前炉次的工况参数、上一炉次的工况参 数、上一炉次的铁水硅含量作为输入，相应当前炉次铁水硅含量作为输出，建立具有输出自 反馈的动态ARMA模型；
[0020] 稀疏化模块:通过提取训练数据集的映射集的极大无关组实现训练数据集的稀疏 化，构造稀疏最小二乘支持向量回归机模型；
[0021 ]鲁棒化模块:在稀疏最小二乘支持向量回归机模型中引入IGGIII加权函数，基于 稀疏化的训练数据集构造具有稀疏鲁棒性的最小二乘支持向量机模型，对训练数据集进行 回归建模，得到高炉炼铁过程铁水硅含量动态软测量模型；
[0022] 多目标遗传优化模块：以评价指标为适应度函数优化高炉炼铁过程铁水硅含量动 态软测量模型；
[0023] 模型预测模块:利用高炉炼铁过程铁水硅含量动态软测量模型预测当前炉次铁水 娃含量。
[0024] 利用所述高炉炼铁过程铁水硅含量动态软测量系统进行高炉炼铁过程铁水硅含 量动态软测量的方法，包括：
[0025] 步骤1、获取高炉铁水硅含量动态软测量所需参数，包括:当前炉次的工况参数、上 一炉次的工况参数和上一炉次的铁水硅含量；
[0026] 步骤2、对获取的高炉铁水硅含量动态软测量所需参数进行归一化预处理；
[0027] 步骤3、采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量动态软测 量:高炉炼铁过程铁水硅含量动态软测量模型采用输出自反馈，以高炉铁水硅含量动态软 测量所需参数为输入，以当前炉次硅含量为输出值，动态预测铁水硅含量。
[0028]所述工况参数的确定方法如下：
[0029] 利用主成分分析法对若干炉次的高炉过程工况参数进行相关性分析，确定对于高 炉铁水硅含量动态软测量贡献率最大的6个工况参数作为辅助变量，包括:炉腹煤气量、热 风温度、热风压力、富氧率、鼓风湿度、喷煤量；
[0030] 基于上述6个辅助变量、结合不同时刻输入参数时滞、当前炉次辅助变量的测量 值、上一炉次辅助变量的测量值及上一炉次铁水硅含量的测量值，确定高炉铁水硅含量动 态软测量所需参数即高炉炼铁过程铁水硅含量动态软测量模型的输入变量，包括：当前炉 次的工况参数、上一炉次的工况参数和上一炉次的铁水硅含量;高炉铁水硅含量动态软测 量所需的工况参数包括炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、设定喷煤量。
[0031] 所述高炉炼铁过程铁水硅含量动态软测量模型的建立方法如下：
[0032] 步骤3-1、采集历史若干炉次的高炉铁水硅含量动态软测量所需参数及相应炉次 铁水硅含量，并进行归一化预处理；
[0033]步骤3-2、将归一化预处理后的高炉铁水硅含量动态软测量所需参数作为训练数 据集；
[0034] 步骤3-3、将训练数据集中的当前炉次的工况参数、上一炉次的工况参数、上一炉 次的铁水硅含量作为输入，相应当前炉次铁水硅含量作为输出，建立具有输出自反馈的动 态ARMA模型：
[0035] 步骤3-4、通过提取训练数据集的映射集的极大无关组实现训练数据集的稀疏化， 构造稀疏最小二乘支持向量回归机模型，作为初始的高炉炼铁过程铁水硅含量动态软测量 模型；
[0036]步骤3-5、在稀疏最小二乘支持向量回归机模型中引入IGGIII加权函数，基于稀疏 化的训练数据集构造具有稀疏鲁棒性的最小二乘支持向量机模型；
[0037] 步骤3-6、利用具有稀疏鲁棒性的最小二乘支持向量机模型对训练数据集进行回 归建模，得到高炉炼铁过程铁水硅含量动态软测量模型；
[0038] 步骤3-7、确定高炉炼铁过程铁水硅含量动态软测量模型的多目标评价指标，以评 价指标为适应度函数，优化高炉炼铁过程铁水硅含量动态软测量模型。
[0039] 所述构造稀疏最小二乘支持向量回归机模型的方法如下：
[0040] 步骤3-4-1、将训练数据集映射到高维希尔伯特空间，得到相应的训练数据集即映 射集；
[0041] 步骤3-4-2、求解映射集的极大无关组，得到稀疏化的训练数据集：
[0042]步骤3-4-3、利用稀疏化的训练数据集，建立稀疏最小二乘支持向量回归机模型。 [0043] 有益效果：
[0044] 为了解决以上高炉冶炼过程铁水硅含量测量鲁棒性的不足，本发明利用具有稀疏 鲁棒性的最小二乘支持向量回归机（Sparse and Robust Least Squares Support Vector Regression，S-R-LS_SVR)与多目标遗传（Multi-Objective Genetic 0ptimization，M0GA) 技术，构造一个具有较强鲁棒性和输出自反馈结构、并考虑不同时刻输入输出数据时滞的 铁水硅含量动态软测量模型，同时利用现场采集的高炉数据对铁水硅含量进行动态软测 量。本发明针对高炉炼铁过程的关键工艺指标一一铁水硅含量难以直接在线测量且离线检 验过程滞后的难题，通过提取动态软测量模型的输入样本数据在其映射集的极大无关组来 实现样本的稀疏化，并将IGGIII加权函数引入，构造具有稀疏鲁棒性的最小二乘支持向量 机模型。与现有的人工测量或化验铁水硅含量相比，减少了操作人员的工作量，降低了人工 操作引入的测量的不确定性，提高了测量的时效性和准确性，可信度高。另外，本发明方法 对高炉冶炼过程的铁水硅含量预报具有普遍的通用性，可获得较好的预报精度，有助于实 现高炉铁水质量的闭环一体化控制与优化运行。
【附图说明】
[0045] 图1是本发明【具体实施方式】中高炉炼铁过程的测量仪表配置图；
[0046] 图2是本发明【具体实施方式】中高炉炼铁过程铁水硅含量动态软测量系统结构框 图；
[0047] 图3是本发明【具体实施方式】中高炉炼铁过程铁水硅含量动态软测量方法流程图；
[0048] 图4是本发明【具体实施方式】中高炉炼铁过程铁水硅含量动态软测量模型的建立方 法流程图；
[0049] 图5是本发明【具体实施方式】中铁水硅含量动态软测量效果图；
[0050] 图1中：1-高炉，2-热风炉，3-流量计，4-温度计，5压力计，6-湿度计，7-炉腹煤气量 测量分析仪，8-富氧率测量分析仪，9-数据采集装置，10-计算机系统。
【具体实施方式】
[0051] 下面结合附图对本发明的【具体实施方式】做详细说明。
[0052] 本实施方式中，实施高炉炼铁过程铁水硅含量动态软测量方法采用测量系统以及 计算机系统10;如图1所示，所述的测量系统包括：
[0053]三个流量计3,分别用于在线测量高炉1煤粉喷吹系统煤粉喷吹量、富氧流量、冷风 流量；
[0054] -个温度计4,用于在线测量高炉1热风炉2的热风温度；
[0055] 一个压力计5,用于在线测量高炉1热风炉2的热风压力；
[0056] -个湿度计6,用于在线测量高炉1热风炉2的鼓风湿度。
[0057]另外，测量系统还包括如下两个测量分析仪：
[0058] -个炉腹煤气量测量分析仪7通过流量计3测量得到的冷风流量、富氧流量以及煤 粉喷吹量，以及湿度计6测量得到的鼓风湿度，分析计算出炉腹煤气量参数；
[0059] -个富氧率测量分析仪8通过流量计3测量得到的冷风流量、富氧流量，以及湿度 计6测量得到的鼓风湿度，分析计算出富氧率参数。
[0060] 流量计3、压力计5、温度计4、湿度计6这些常规测量仪表安装于高炉冶炼系统的各 个相应位置。数据采集器9的输入端连接这些常规测量仪表的输出端，并通过通信总线连接 计算机系统10。
[0061] 计算机系统10中设有高炉炼铁过程铁水硅含量动态软测量系统，如图2所示，包 括：
[0062] 实际数据采集单元:获取高炉铁水硅含量动态软测量所需参数，包括：当前炉次的 工况参数、上一炉次的工况参数和上一炉次的铁水硅含量；
[0063] 归一化预处理单元:对获取的高炉铁水娃含量动态软测量所需参数进行归一化预 处理；
[0064]动态软测量单元:采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量 动态软测量:高炉炼铁过程铁水硅含量动态软测量模型采用输出自反馈，以高炉铁水硅含 量动态软测量所需参数为输入，以当前炉次硅含量为输出值，动态预测铁水硅含量。
[0065]实际数据采集单元，包括：
[0066] 辅助变量确定模块:利用主成分分析法对若干炉次的高炉过程工况参数进行相关 性分析，确定对于高炉铁水硅含量动态软测量贡献率最大的工况参数作为辅助变量，具体 包括:炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、喷煤量；
[0067] 输入变量确定模块:基于辅助变量、结合不同时刻输入参数时滞、当前炉次辅助变 量的测量值、上一炉次辅助变量的测量值及上一炉次铁水硅含量的测量值，确定高炉铁水 硅含量动态软测量所需参数即高炉炼铁过程铁水硅含量动态软测量模型的输入变量，包 括:当前炉次的工况参数、上一炉次的工况参数和上一炉次的铁水硅含量;高炉铁水硅含量 动态软测量所需的工况参数包括炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、设定 喷煤量。
[0068] 输入变量采集模块:采集高炉炼铁过程铁水硅含量动态软测量模型的输入变量。 [0069] 动态软测量单元，包括：
[0070]样本处理模块:采集历史若干炉次的高炉铁水硅含量动态软测量所需参数及相应 炉次铁水硅含量，并进行归一化预处理后作为训练数据集；
[0071 ] ARMA模型构建模块:将训练数据集中的当前炉次的工况参数、上一炉次的工况参 数、上一炉次的铁水硅含量作为输入，相应当前炉次铁水硅含量作为输出，建立具有输出自 反馈的动态ARMA模型；
[0072]稀疏化模块:通过提取训练数据集的映射集的极大无关组实现训练数据集的稀疏 化，构造稀疏最小二乘支持向量回归机模型；
[0073]鲁棒化模块:在稀疏最小二乘支持向量回归机模型中引入IGGIII加权函数，基于 稀疏化的训练数据集构造具有稀疏鲁棒性的最小二乘支持向量机模型，对训练数据集进行 回归建模，得到高炉炼铁过程铁水硅含量动态软测量模型；
[0074] 多目标遗传优化模块：以评价指标为适应度函数优化高炉炼铁过程铁水硅含量动 态软测量模型；
[0075] 模型预测模块:利用高炉炼铁过程铁水硅含量动态软测量模型预测当前炉次铁水 娃含量。
[0076]本实施方式的高炉炼铁过程铁水硅含量动态软测量方法，如图3所示，包括以下步 骤：
[0077]步骤1、获取高炉铁水硅含量动态软测量所需参数，包括:当前炉次的工况参数、上 一炉次的工况参数和上一炉次的铁水硅含量；
[0078]所述工况参数的确定方法如下：
[0079] 步骤1-1、利用主成分分析法（Principal Component Analysis，PCA)对若干炉次 的高炉过程工况参数进行相关性分析，确定对于高炉铁水硅含量动态软测量贡献率最大的 6个工况参数作为辅助变量，具体包括：炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿 度、喷煤量；
[0080] 利用主成分分析法对若干炉次的高炉过程工况参数进行相关性分析，具体如下：
[0081] (1)、特征中心化，即对输入样本数据XnXm=[Xl，X2, . . .，Xm]每一列的数据都减去该 列的均值，进而得到矩阵六#?。
[0082] (2)、计算AnXm的协方差矩阵BmXm，并计算协方差矩阵B mXm的特征向量和特征值。
[0083] (3)、选取若干大的特征值所对应的特征向量作为新的数据。
[0084]其中炉腹煤气量m，m3、热风温度u2，°C、热风压力u3，KPa、富氧率U4，％、鼓风湿度 115，冊、喷煤量116，1113/11的特征值之和占所有特征值的98.723%>98%，因此选择该6类变量 构成新的样本集作为铁水硅含量动态软测量的辅助变量；
[0085] 步骤1-2、基于上述6个辅助变量、结合不同时刻输入参数时滞、当前炉次辅助变量 的测量值、上一炉次辅助变量的测量值及上一炉次铁水硅含量的测量值，确定高炉铁水硅 含量动态软测量所需参数即高炉炼铁过程铁水硅含量动态软测量模型的输入变量，包括： 当前炉次的工况参数、上一炉次的工况参数和上一炉次的铁水硅含量;高炉铁水硅含量动 态软测量所需的工况参数包括炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、设定喷 煤量。
[0086] 确定如下13个变量为动态软测量模型的输入变量：
[0088]步骤2、对获取的高炉铁水硅含量动态软测量所需参数进行归一化预处理；
[0089]归一化预处理的具体实现方法为：
[0091]其中Xl为归一化之前的高炉炼铁过程铁水硅含量动态软测量模型的输入变量，< 为归一化之后的高炉炼铁过程铁水硅含量动态软测量模型的输入变量，m为^的均值，〇1为 Xi的标准差。
[0092] 步骤3、采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量动态软测 量:高炉炼铁过程铁水硅含量动态软测量模型采用输出自反馈，以高炉铁水硅含量动态软 测量所需参数为输入，以当前炉次硅含量为输出值，动态预测铁水硅含量。
[0093] 如图4所示，高炉炼铁过程铁水硅含量动态软测量模型的建立方法如下：
[0094] 步骤3-1、采集历史若干炉次的高炉铁水硅含量动态软测量所需参数及相应炉次 铁水硅含量，并进行归一化预处理；
[0095]步骤3-2、将归一化预处理后的高炉铁水硅含量动态软测量所需参数作为训练数 据集；
[0096] 步骤3-3、将训练数据集中的当前炉次的工况参数、上一炉次的工况参数、上一炉 次的铁水硅含量作为输入，相应当前炉次铁水硅含量作为输出，建立具有输出自反馈的动 态ARMA(Auto Regressive Moving Average，ARMA)模型：
[0097] y(t) = fARMA{ui(t)，U2(t)，U3(t)，U4(t)，U5(t)，U6(t)，
[0098] Ui(t-l)，U2(t-l)，U3(t-l)，U4(t-l)，U5(t-l)，U6(t-l)，y(t_l)}
[0099] 其中y(t)为当前炉次铁水硅含量的真实值。
[0100] 步骤3-4、通过提取训练数据集的映射集的极大无关组实现训练数据集的稀疏化， 构造稀疏最小二乘支持向量回归机模型，作为初始的高炉炼铁过程铁水硅含量动态软测量 丰旲型；
[0101] 通过提取训练数据集在其映射集的极大无关组，来实现训练数据集的稀疏化，降 低高炉炼铁过程铁水硅含量动态软测量模型的复杂度，构造稀疏最小二乘支持向量回归机 模型（Sparse Least Squares Support Vector Regression，S-LS_SVR)，具体实现方法为：
[0102] 步骤3-4-1、将训练数据集映射到高维希尔伯特空间，得到相应的训练数据集即映 射集；
[0103] 训练数据集{.w,G通过非线性函数映射到高维希尔伯特空间，其映射集为 4 = ，N为样本数量，X e RNX13为由步骤3-3所确定的高炉炼铁过程铁水硅含量动态软测 量模型的13个输入变量组成的输入样本矩阵，yi为步骤3-3所确定的动态ARMA模型输出，即
铁水硅含量，Ph)为非线性函数。虽然不能被确切的表达 m为映射集中线性相关的个数，λι，λ」均为加权系数，取非线性函数4?)为高斯函数K(x，xi)， 定义如下：
[0105] 其中，Xi为第i组输入样本，〇为高斯(径向基)函数的伸缩量。
[0106] 步骤3-4-2、求解映射集Α的极大无关组，得到稀疏化的训练数据集：
[0107] (a)初始化极大无关组Βμ= Φ，在集合S= (1，2, . . .，N)选取j = l时的输入样本放 到Bm中；
[0108] (b)在S中依次选取j = j+Ι，计算
[0110] (c)若目标函数mjn C5(A)<s，ε为设定的阈值，则说明_χ,)可以由冰(Λ；;ψββ Λ/}线性 表示，摒弃第i组输入样本;若目标函数n￥n 则说明列4)不可由丨线性 表示，则MX )}线性无关，将第i组输入样本放到集合BM中；
[0111] (d)若迭代次数j<N，则转到(b);否则终止迭代；
[0112] (e)将极大无关组Bm中对应训练数据集的元素取出组成新训练数据集 ，凡为新训练数据集的样本数量。Ψ5通过高斯函数映射后的集合为 象=加(匕,)丨U，即稀疏化的训练数据集；
[0113] 因为Ψ5是映射集Α的极大无关组，则权值向量ω eRn可表述为
[0115]步骤3-4-3、利用稀疏化的训练数据集，建立稀疏最小二乘支持向量回归机模型 为：
[0117] 其中，C为正则化参数，b为偏置项，吟=?.-鳥_为误差项，yi，i分别为铁水娃含量的 实际测量值和动态软测量的预测值。
[0118] 步骤3-5、在稀疏最小二乘支持向量回归机模型中引入IGGIII加权函数，基于稀疏 化的训练数据集构造具有稀疏鲁棒性的最小二乘支持向量机模型（Sparse robust Least Squares Support Vector Regression，S-R_LS_SVR);
[0119] 为了改善初始的高炉炼铁过程铁水硅含量动态软测量模型的鲁棒性能，即实际情 况与高炉炼铁过程铁水硅含量动态软测量模型有较严重的偏离时也能保证估计精度不受 破坏性的影响。因此，对误差项 ei引入加权因子Vl，从而得到具有稀疏鲁棒性的最小二乘支 持向量机模型：
[0121 ]式中：Vi = diag(vi，V2, · · ·，Vr)，vi由IGGIII加权函数决定，即：
[0123] 其中，θ = 1.438median | e-median(e) |为误差的标准差估计，median( ·)为取中位 数运算，ki，k2为相关系数，根据经验值kie [1，3]，k2e [3.0,6]，本实施方式选取ki = l .3，k2 =3.2〇
[0124] 步骤3-6、利用具有稀疏鲁棒性的最小二乘支持向量机模型对训练数据集进行回 归建模，得到高炉炼铁过程铁水硅含量动态软测量模型。
[0125] 优化的高炉炼铁过程铁水硅含量动态软测量模型中引入拉格朗日算子a=(ai， a2，· · ·，ar)得：
[01 Z7]根据极小值的最优条件，即令L(P，b，e，a)对β，b，e，a的偏导数为零，并消去变量e， a得：
[0129] 进而获得高炉炼铁过程铁水硅含量动态软测量模型：
[0130] f(xi) = ΦΓ(χ?) ω+b
[0131 ]步骤3-7、确定高炉炼铁过程铁水硅含量动态软测量模型的多目标评价指标，以评 价指标为适应度函数，优化高炉炼铁过程铁水硅含量动态软测量模型；
[0132]传统的建模性能指标大多采用均方根误差(Root Mean Square Error，RMSE)，然 而这只能表明建模误差尽可能小，对模型预测曲线与真实曲线的接近程度并未考虑。在实 际应用中，准确的变化趋势对动态过程的建模至关重要。为此，提出综合考虑RMSE和估计曲 线与实际曲线相关性的模型精度多目标评价指标。
[0133]由数理统计理论可知，两个数据矢量X和Y分别表示高炉炼铁过程铁水硅含量实际 值和软测量值，E[(X-E(X))(Y-E(Y))]称之为X与Y的协方差或者相关矩，记作C0V(X，Y)，E (X)，E(Y)分别是X和Y的期望，而该两个数据矢量的相关系数为
[0135]
是X和Y的方差。相关系数Ρχγ是衡量数据矢量X与Y关系程度的 量：| Ρχγ卜1表不X与Y之间的线性关系密切;而| Ρχγ | -0表不X与Y的相关性很差;若| Ρχγ | = 1 表示X与Υ依概率1存在着线性关系，而I ρχγ I -0表示X与Υ不相关。
[0136] 综上，高炉炼铁过程铁水硅含量动态软测量模型的多目标评价指标如下：
[0138] Fcei2= 1~Ρχγ
[0139] 其中，只分别为铁水硅含量的预测值与实际值。上述两个评价指标既 可以保证建模过程的平稳性和限制输出曲线的横向偏移量，又可以保证建模过程的准确性 和限制输出曲线的纵向偏移量。
[0140]以上述多目标评价指标为适应度函数，利用多目标遗传算法对高炉炼铁过程铁水 硅含量动态软测量模型S-R-LS-SVR进行优化:利用NSGA-II对正则化参数C和高斯核函数的 伸缩量σ进行优化。以模型的多目标评价指标为适应度函数，采用实数编码，通过基于进行 非支配快速排序和拥挤距离计算的种群进行二进制锦标赛选择，利用模拟二进制交叉和多 项式变异的遗传因子增强种群多样性，具体如下所示：
[0141] (a)初始化。设置当前代num=l，最大迭代次数为maxnum，种群大小为Q，变异率为 0.1，交叉率为0.9。
[0142] (b)实数编码。对正则化参数C和高斯核函数的伸缩量〇的解空间进行浮点数编码， 每条染色体对应一组(：和〇，编码形式为染色体R=[0，C]，该条染色体的基因数量为S = 2,含 有Q条随机生成的染色体的初始种群为P=[Ri，R2, . . .，Rq]t。将优化的参数C和σ作为基因放 在染色体中。
[0143] (c)计算多目标评价指标。根据种群Ρ中的参数(：和〇计算Q个个体的多目标评价指 标。
[0144] (d)判断是否满足终止条件。如果num>最大迭代次数maxnum，则优化结束；否则， 进入步骤(e)。
[0145] (e)计算拥挤距离并进行非支配排序。
[0146] (f)锦标赛选择算子。根据非支配排序的结果，选择支配层较低的染色体，若同一 支配层的染色体有多个，选择拥挤距离较大的染色体以获得种群的多样性。
[0147] (g)模拟二进制交叉。由于采用实数编码，则交叉后代是父代的线性组合：
[0149] 式中：u为产生于(0，1)均匀分布的随机数。当u>0.5时
；当11 ^0.5
，n。为交叉分布指数，i = 1，2为优化目标函数的个数。
[0150] (h)多项式变异。变异后的染色体为，其中B'B1分别为优 化变量的上界与下界，为变异的参数。当rk> 0.5时，
；当^<0.5，
:，rk为来自（0，1)均匀分布的随机数，为变异分布指数。
[0151] (i)重组和优选。将经上述模拟二进制交叉和多项式变异产生的新染色体与原种 群中的染色体进行混合重组，并选择非支配排序前N个染色体组成新的种群。num = num+l， 返回步骤(c)。
[0152] 为验证本发明的有效性，以炼铁厂容积为2600m3的炼铁高炉为实施对象。该高炉 对象安装了如下的测量系统，包括：
[0153] 横河DPharp EJA系列压力变送器用于测量高炉热风系统的热风压力；
[0154] HH-WLB差压流量计用于测量冷风流量；
[0155] A+K平衡流量计用于测量富氧流量；
[0156] JWSK-6CWDA空气湿度传感器用于测量鼓风湿度；
[0157] YHIT红外测温仪用于测量热风温度；
[0158] HDLWG-06煤粉流量计用于测量煤粉喷吹量.
[0159] 另外，测量系统还包括如下两个测量分析仪：
[0160] -个测量分析仪通过流量计测量得到的冷风流量、富氧流量以及煤粉喷吹量，以 及湿度计测量得到的鼓风湿度，分析计算出炉腹煤气量参数；
[0161]炉腹煤气量测量分析仪设置如下：
[0163]富氧率测量分析仪设置如下：
[0165] 实验中共采集了270炉次的样本数据，其中前200炉次生产数据作为建模时的训练 样本，后70炉次的数据作为测试样本。在模型训练和模型预测过程中使用的所有变量的采 样数据，均采用以铁水出炉炉次为单位的测量平均值作为采样和预报周期。经PCA进行相关 性分析，确定炉腹煤气量m、热风温度11 2、热风压力U3、富氧率U4、鼓风湿度U5、喷煤量U6六个 辅助变量;采用输出自反馈的动态ARMA模型，确定13个输入变量分别为:本炉次炉腹煤气量 111(〇、本炉次热风温度112(〇、本炉次水热风压力11 3(〇、本炉次富氧率114(〖）、本炉次鼓风湿 度1!5(〇、本炉次设定喷煤量U6(t)、上炉次炉腹煤气量m(t-l)、上炉次热风温度1! 2(卜1)、上 炉次水热风压力U3(t-1)、上炉次富氧率U4(t-1)、上炉次鼓风湿度115(〖-1)、上炉次设定喷煤 量 U6(t-1)、上炉次Si含量估计值Κ?-ι);对上述数据进行一致的平滑处理和归一化处理。然 后通过输入样本数据在其映射集的极大无关组实现输入数据集的稀疏化，降低软测量模型 的复杂度，并将IGGIII加权函数引入建模过程构造具有稀疏鲁棒的软测量模型，最终利用 对训练样本进行回归建模。同时在使用带有精英策略的多目标遗传算法优化正则项参数C 和径向基核函数的伸缩量σ时，采用的种群大小为30,最大进化代数为50,交叉率为0.9,变 异率为〇 . 1，初始化范围分别为λ e (〇 . 1，20 )，σ e (〇 . 1，1 〇 )。优化结果为：λ = 9.8，σ = 0.8631〇
[0166] 图5为铁水硅含量预报值跟踪真实检验值的效果，可以看出铁水硅含量动态软测 量预测值与实际值基本一致，误差较小，且变化趋势基本一致。此外，本方法有较强的鲁棒 性和泛化能力，相比其他方法有较高的优越性与先进性。因此是一种具有较高实用价值的、 低成本的、可连续测量的高炉铁水硅含量动态软测量方法。
[0167] 可以理解的是，以上关于本发明的具体描述，仅用于说明本发明而并非受限于本 发明实施例所描述的技术方案，本领域的普通技术人员应当理解，仍然可以对本发明进行 修改或等同替换，以达到相同的技术效果；只要满足使用需要，都在本发明的保护范围之 内。
【主权项】
1. 一种高炉炼铁过程铁水硅含量动态软测量系统，其特征在于，包括： 实际数据采集单元:获取高炉铁水硅含量动态软测量所需参数，包括：当前炉次的工况 参数、上一炉次的工况参数和上一炉次的铁水硅含量； 归一化预处理单元：对获取的高炉铁水娃含量动态软测量所需参数进行归一化预处 理； 动态软测量单元:采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量动态 软测量:高炉炼铁过程铁水硅含量动态软测量模型采用输出自反馈，以高炉铁水硅含量动 态软测量所需参数为输入，以当前炉次硅含量为输出值，动态预测铁水硅含量。2. 根据权利要求1所述的高炉炼铁过程铁水硅含量动态软测量系统，其特征在于，所述 实际数据采集单元，包括： 辅助变量确定模块:利用主成分分析法对若干炉次的高炉过程工况参数进行相关性分 析，确定对于高炉铁水硅含量动态软测量贡献率最大的工况参数作为辅助变量，具体包括： 炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、喷煤量； 输入变量确定模块:基于辅助变量、结合不同时刻输入参数时滞、当前炉次辅助变量的 测量值、上一炉次辅助变量的测量值及上一炉次铁水硅含量的测量值，确定高炉铁水硅含 量动态软测量所需参数即高炉炼铁过程铁水硅含量动态软测量模型的输入变量，包括：当 前炉次的工况参数、上一炉次的工况参数和上一炉次的铁水硅含量;高炉铁水硅含量动态 软测量所需的工况参数包括炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、设定喷煤 量； 输入变量采集模块:采集高炉炼铁过程铁水硅含量动态软测量模型的输入变量。3. 根据权利要求1所述的高炉炼铁过程铁水硅含量动态软测量系统，其特征在于，所述 动态软测量单元，包括： 样本处理模块:采集历史若干炉次的高炉铁水硅含量动态软测量所需参数及相应炉次 铁水硅含量，并进行归一化预处理后作为训练数据集； ARMA模型构建模块:将训练数据集中的当前炉次的工况参数、上一炉次的工况参数、上 一炉次的铁水硅含量作为输入，相应当前炉次铁水硅含量作为输出，建立具有输出自反馈 的动态ARM模型； 稀疏化模块:通过提取训练数据集的映射集的极大无关组实现训练数据集的稀疏化， 构造稀疏最小二乘支持向量回归机模型； 鲁棒化模块:在稀疏最小二乘支持向量回归机模型中引入IGGIII加权函数，基于稀疏 化的训练数据集构造具有稀疏鲁棒性的最小二乘支持向量机模型，对训练数据集进行回归 建模，得到高炉炼铁过程铁水硅含量动态软测量模型； 多目标遗传优化模块：以评价指标为适应度函数优化高炉炼铁过程铁水硅含量动态软 测量模型； 模型预测模块:利用高炉炼铁过程铁水硅含量动态软测量模型预测当前炉次铁水硅含 量。4. 利用权利要求1所述高炉炼铁过程铁水硅含量动态软测量系统进行高炉炼铁过程铁 水硅含量动态软测量的方法，其特征在于，包括： 步骤1、获取高炉铁水硅含量动态软测量所需参数，包括：当前炉次的工况参数、上一炉 次的工况参数和上一炉次的铁水硅含量； 步骤2、对获取的高炉铁水硅含量动态软测量所需参数进行归一化预处理； 步骤3、采用高炉炼铁过程铁水硅含量动态软测量模型进行铁水硅含量动态软测量:高 炉炼铁过程铁水硅含量动态软测量模型采用输出自反馈，以高炉铁水硅含量动态软测量所 需参数为输入，以当前炉次硅含量为输出值，动态预测铁水硅含量。5. 根据权利要求4所述的高炉炼铁过程铁水硅含量动态软测量的方法，其特征在于，所 述工况参数的确定方法如下： 利用主成分分析法对若干炉次的高炉过程工况参数进行相关性分析，确定对于高炉铁 水娃含量动态软测量贡献率最大的6个工况参数作为辅助变量，包括:炉腹煤气量、热风温 度、热风压力、富氧率、鼓风湿度、喷煤量； 基于上述6个辅助变量、结合不同时刻输入参数时滞、当前炉次辅助变量的测量值、上 一炉次辅助变量的测量值及上一炉次铁水硅含量的测量值，确定高炉铁水硅含量动态软测 量所需参数即高炉炼铁过程铁水硅含量动态软测量模型的输入变量，包括：当前炉次的工 况参数、上一炉次的工况参数和上一炉次的铁水硅含量;高炉铁水硅含量动态软测量所需 的工况参数包括炉腹煤气量、热风温度、热风压力、富氧率、鼓风湿度、设定喷煤量。6. 根据权利要求4所述的高炉炼铁过程铁水硅含量动态软测量的方法，其特征在于，所 述高炉炼铁过程铁水硅含量动态软测量模型的建立方法如下： 步骤3-1、采集历史若干炉次的高炉铁水硅含量动态软测量所需参数及相应炉次铁水 娃含量，并进行归一化预处理； 步骤3-2、将归一化预处理后的高炉铁水硅含量动态软测量所需参数作为训练数据集； 步骤3-3、将训练数据集中的当前炉次的工况参数、上一炉次的工况参数、上一炉次的 铁水硅含量作为输入，相应当前炉次铁水硅含量作为输出，建立具有输出自反馈的动态 ARMA模型： 步骤3-4、通过提取训练数据集的映射集的极大无关组实现训练数据集的稀疏化，构造 稀疏最小二乘支持向量回归机模型，作为初始的高炉炼铁过程铁水硅含量动态软测量模 型； 步骤3-5、在稀疏最小二乘支持向量回归机模型中引入IGGIII加权函数，基于稀疏化的 训练数据集构造具有稀疏鲁棒性的最小二乘支持向量机模型； 步骤3-6、利用具有稀疏鲁棒性的最小二乘支持向量机模型对训练数据集进行回归建 模，得到高炉炼铁过程铁水硅含量动态软测量模型； 步骤3-7、确定高炉炼铁过程铁水硅含量动态软测量模型的多目标评价指标，以评价指 标为适应度函数，优化高炉炼铁过程铁水硅含量动态软测量模型。7. 根据权利要求6所述的高炉炼铁过程铁水硅含量动态软测量的方法，其特征在于，所 述构造稀疏最小二乘支持向量回归机模型的方法如下： 步骤3-4-1、将训练数据集映射到高维希尔伯特空间，得到相应的训练数据集即映射 集； 步骤3-4-2、求解映射集的极大无关组，得到稀疏化的训练数据集： 步骤3-4-3、利用稀疏化的训练数据集，建立稀疏最小二乘支持向量回归机模型。
【文档编号】C21B7/24GK105886680SQ201610308346
【公开日】2016年8月24日
【申请日】2016年5月11日
【发明人】周平, 郭东伟, 宋贺达, 王宏
【申请人】东北大学