第36卷第4期 北京科技大学学报 Vol.36 No.4 2014年4月 Journal of University of Science and Technology Beijing Apr.2014 基于主客观证据融合的高炉悬料预测方法 曹卫华2》,杜楠2”,安剑奇12)四,吴敏2) 1)中南大学信息科学与工程学院,长沙4100832)先进控制与智能自动化湖南省工程实验室,长沙410083 ☒通信作者,E-mail:anjianqi@csu.cd.cn 摘要针对高炉关键异常炉况悬料难以预测的问题,基于DS证据理论,提出一种综合模糊专家推理和后验概率最小二乘 支持向量机的悬料预测方法.首先,结合高炉生产过程和悬料现象,分析悬料形成的内在机理:其次,通过模糊专家推理提取 基于专家规则的主观证据,再通过建立后验概率最小二乘支持向量机模型提取基于数据内在客观规律的客观证据:最后,基 于DS证据理论完成主客观证据融合,实现悬料预测.该方法充分利用专家经验和最小二乘支持向量机的自学习能力,能够 提高预测精度.仿真结果表明本文提出的方法有效、准确。 关键词高炉;预测:模糊逻辑;最小二乘逼近:支持向量机:信总融合 分类号T℉543.1 Prediction method of blast furnace hanging based on fusion of subjective and objective evidences CAO Wei-hua,DU Nan2),AN Jian-,WU Min' 1)School of Information Science and Engineering,Central South University,Changsha 410083,China 2)Hunan Engineering Laboratory for Advanced Control and Intelligent Automation,Changsha 410083.China Corresponding author,E-mail:anjianqi@csu.edu.cn ABSTRACT Aiming at the difficulty of predicting blast furnace hanging,a prediction method was proposed for the hanging based on the D-S evidence theory and in combination with fuzzy expert inference and a posterior probability least squares support vector machine. Firstly,the causes of hanging are obtained by mechanism analysis in consideration of blast furnace operations and hanging phenomena. Secondly,subjective evidences are extracted by fuzzy expert reasoning,while a posterior probability least squares support vector machine model is developed to extract objective evidences.Finally,in order to predict the hanging precisely,the subjective and objec- tive evidences are fused based on the D-S evidence theory,which makes full use of the expertise and the self-learning ability of the least squares support vector machine.Simulation results illustrate that the proposed method can make accurate prediction of the hanging. KEY WORDS blast furnaces;prediction:fuzzy logic:least squares approximations:support vector machines:information fusion 高炉炼铁是钢铁治金生产过程中的关键环主要依靠专家和操作人员的经验,从而不可避免的 节四,维持高炉的稳定高效运行,是整个钢铁生产 受到人为主观因素的影响,出现预测错误的情况 企业正常生产的重要保障.高炉悬料是发生频率较 如何有效地预测高炉悬料,提示操作人员及时采取 高且危害极大的一种异常炉况,会严重破坏高炉的 适当的调控措施,避免悬料的发生或炉况的进一步 正常治炼行程.由于高炉炼铁是在密闭容器内进行 恶化,在实际生产中至关重要 的,其运行状态无法直接检测回,致使悬料的预测 高炉悬料预测方法大致分为两类:一类是基于 收稿日期:20130909 基金项目:湖南省自然科学省市(衡阳)联合基金资助项目(12刀8012):国家自然科学基金资助项目(61203017):国家高技术研究发展计划资 助项目(2012A4040307):国家自然科学基金重点资助项目(61333002) DOI:10.13374/j.issn1001-053x.2014.04.013;http://journals.ustb.edu.en
第 36 卷 第 4 期 2014 年 4 月 北京科技大学学报 Journal of University of Science and Technology Beijing Vol. 36 No. 4 Apr. 2014 基于主客观证据融合的高炉悬料预测方法 曹卫华1,2) ,杜 楠1,2) ,安剑奇1,2) ,吴 敏1,2) 1) 中南大学信息科学与工程学院,长沙 410083 2) 先进控制与智能自动化湖南省工程实验室,长沙 410083 通信作者,E-mail: anjianqi@ csu. edu. cn 摘 要 针对高炉关键异常炉况悬料难以预测的问题,基于 D--S 证据理论,提出一种综合模糊专家推理和后验概率最小二乘 支持向量机的悬料预测方法. 首先,结合高炉生产过程和悬料现象,分析悬料形成的内在机理; 其次,通过模糊专家推理提取 基于专家规则的主观证据,再通过建立后验概率最小二乘支持向量机模型提取基于数据内在客观规律的客观证据; 最后,基 于 D--S 证据理论完成主客观证据融合,实现悬料预测. 该方法充分利用专家经验和最小二乘支持向量机的自学习能力,能够 提高预测精度. 仿真结果表明本文提出的方法有效、准确. 关键词 高炉; 预测; 模糊逻辑; 最小二乘逼近; 支持向量机; 信息融合 分类号 TF543. 1 Prediction method of blast furnace hanging based on fusion of subjective and objective evidences CAO Wei-hua1,2) ,DU Nan1,2) ,AN Jian-qi 1,2) ,WU Min1,2) 1) School of Information Science and Engineering,Central South University,Changsha 410083,China 2) Hunan Engineering Laboratory for Advanced Control and Intelligent Automation,Changsha 410083,China Corresponding author,E-mail: anjianqi@ csu. edu. cn ABSTRACT Aiming at the difficulty of predicting blast furnace hanging,a prediction method was proposed for the hanging based on the D-S evidence theory and in combination with fuzzy expert inference and a posterior probability least squares support vector machine. Firstly,the causes of hanging are obtained by mechanism analysis in consideration of blast furnace operations and hanging phenomena. Secondly,subjective evidences are extracted by fuzzy expert reasoning,while a posterior probability least squares support vector machine model is developed to extract objective evidences. Finally,in order to predict the hanging precisely,the subjective and objective evidences are fused based on the D-S evidence theory,which makes full use of the expertise and the self-learning ability of the least squares support vector machine. Simulation results illustrate that the proposed method can make accurate prediction of the hanging. KEY WORDS blast furnaces; prediction; fuzzy logic; least squares approximations; support vector machines; information fusion 收稿日期: 2013--09--09 基金项目: 湖南省自然科学省市( 衡阳) 联合基金资助项目( 12JJ8012) ; 国家自然科学基金资助项目( 61203017) ; 国家高技术研究发展计划资 助项目( 2012AA040307) ; 国家自然科学基金重点资助项目( 61333002) DOI: 10. 13374 /j. issn1001--053x. 2014. 04. 013; http: / /journals. ustb. edu. cn 高炉炼铁是钢铁冶金生产过程中的关键环 节[1],维持高炉的稳定高效运行,是整个钢铁生产 企业正常生产的重要保障. 高炉悬料是发生频率较 高且危害极大的一种异常炉况,会严重破坏高炉的 正常冶炼行程. 由于高炉炼铁是在密闭容器内进行 的,其运行状态无法直接检测[2],致使悬料的预测 主要依靠专家和操作人员的经验,从而不可避免的 受到人为主观因素的影响,出现预测错误的情况. 如何有效地预测高炉悬料,提示操作人员及时采取 适当的调控措施,避免悬料的发生或炉况的进一步 恶化,在实际生产中至关重要. 高炉悬料预测方法大致分为两类: 一类是基于
第4期 曹卫华等:基于主客观证据融合的高炉悬料预测方法 ·507· 模糊专家推理的主观专家系统方法:另一类是 应,目的是将铁矿石治炼成合格的生铁回.即炉料 基于统计学习的客观智能预测方法,包括神经网 在下降的过程中,与上升的煤气流发生氧化还原反 络B-、支持向量机等方法.其中,专家系统利 应,生成液态生铁.悬料是指由于炉内料柱透气性 用产生式规则表示专家知识,采用基于规则匹配的 与煤气流运动极不适应而导致炉料停止下降的炉况 模糊专家推理方法进行推理,实现悬料预测.专家 失常现象0,是高炉治炼生产过程中时常发生的异 系统中的规则是从炼铁知识和操作规程中提炼得到 常炉况,严重影响高炉冶炼生产的正常进行 的,是专家在长期生产实践中总结出来的知识,主观 1.1高炉生产过程 上反映了悬料的本质特征.但是高炉治炼环境时刻 高炉冶炼是将铁矿石还原成生铁的连续生产过 变化,不同时刻的悬料模式差异较大,而专家系统缺 程.烧结矿、焦炭、石灰石等固体原料按规定配比通 乏有效的知识更新手段,难以及时获取高炉悬料的 过高炉炉顶布料装置投入炉内形成炉料,炉料在重 新模式,自适应性不足.基于统计学习的智能预测 力作用下向下运动.热风炉输出的热风经下部风口 方法,利用悬料样本数据进行训练和学习,通过学习 鼓入,并在风口回旋区与喷吹燃料和风口焦炭燃烧 算法建立预测模型,客观地抽取出每一个输入参数 生成煤气:煤气被炉缸与炉顶压力差推动经炉料间 对悬料影响的权值因子,实现对悬料的预测.该类 隙向上运动,进入高炉项部.炉料在逐步下降的过 方法应用灵活,具有较强的自学习能力和适应性;但 程中,受到逆流向上运动的高温煤气的还原作用,不 其仅从统计学角度出发,挖掘数据之间存在的客观 断地被加热和还原而成为熔融状态,并以液态形式 规律,缺乏与高炉己知经验规律的结合,模型精度严 通过焦炭层的空隙向下滴落,最终在炉缸形成炉渣 重依赖于悬料样本数据,由于悬料样本数据属于小 和铁水. 样本集,难以保证悬料样本的泛化性,致使预测模型 1.2悬料成因分析 难以客观的完全反映悬料本质规律,预测效果不佳 悬料是高炉炉热、炉凉和管道等异常炉况发展 高炉治炼过程复杂,存在众多干扰,造成检测信 到一定程度的必然结果.从炼铁工艺角度分析,炉 息的客观不确定性,而且悬料和参数之间的异常状 内料柱透气性因炉温过热而变差,导致风压增大,致 态之间没有确定的对应关系,也存在很大的客观不 使热风造成的压差增大而托起料柱:同时在炉墙上 确定性.同时,专家系统由于自适应不足,而且模糊 黏结的软熔炉料,阻止了边缘煤气流的上升和料柱 隶属函数的设定通常基于专家经验,具有很大主观 的下降,进而形成悬料:从流体力学角度分析,炉内 性,致使推理结果存在主观不确定性:智能预测方法 煤气流向上运动的进程受阻,在高炉下部区域内形 由于样本泛化性不足导致预测结果存在客观不确 成湍流并产生巨大的托力,从而阻碍炉料的下行,最 定性. 终形成悬料 针对悬料与参数之间,以及专家系统和智能预 通过上述两个角度的分析,总结悬料形成的内 测方法存在主、客观不确定性的问题,本文结合DS 在机理如下:炉内煤气流的通道突然被阻塞,造成料 证据理论,提出一种基于主客观证据融合的高炉悬 柱某一层面的煤气压力快速升高,加大了对该层面 料预测方法.首先,结合高炉生产过程分析悬料成 以上料柱的托力,随着托力增加导致该层面以上的 因,并综合现场生产经验和数据分析,确定悬料预测 料柱停止运动,从而形成悬料,该层面以下的料柱正 输入参数.在此基础上,一方面基于专家经验,以专 常下降 家规则的形式描述征兆参数与悬料的关系,通过模 糊专家推理,提取主观证据对悬料的可信度;另一方 2基于主客观证据融合的悬料预测方案 面,基于后验概率最小二乘支持向量机(least 正常炉况时,高炉各项过程检测参数在正常范 squares support vector machine,LSSVM)模型,利用 围内小幅波动,变化趋势相适应(如风量和风压的 LSSVM的自学习能力,从统计学角度挖掘征兆参数 变化趋势相一致等).当即将发生异常炉况时,部分 与悬料间存在的客观联系,提取客观证据对悬料的 检测参数会有明显的异常波动,经过现场专家多年 可信度.最后,基于D一S证据理论对主、客观证据 的总结分析,悬料发生前的主要征兆有:(1)风压逐 可信度进行融合推理,实现对悬料的预测 渐上升,风量缓慢减少;(2)料速持续减小,料尺越 来越宽;(3)料柱透气性变差,炉顶压力降低,全压 ,高炉生产过程与悬料成因分析 差升高 高炉生产过程伴随着极其复杂的物理、化学反 结合专家经验,通常选取风量、风压、透气性指
第 4 期 曹卫华等: 基于主客观证据融合的高炉悬料预测方法 模糊专家推理的主观专家系统方法[3--4]; 另一类是 基于统计学习的客观智能预测方法,包括神经网 络[5--6]、支持向量机[7--8]等方法. 其中,专家系统利 用产生式规则表示专家知识,采用基于规则匹配的 模糊专家推理方法进行推理,实现悬料预测. 专家 系统中的规则是从炼铁知识和操作规程中提炼得到 的,是专家在长期生产实践中总结出来的知识,主观 上反映了悬料的本质特征. 但是高炉冶炼环境时刻 变化,不同时刻的悬料模式差异较大,而专家系统缺 乏有效的知识更新手段,难以及时获取高炉悬料的 新模式,自适应性不足. 基于统计学习的智能预测 方法,利用悬料样本数据进行训练和学习,通过学习 算法建立预测模型,客观地抽取出每一个输入参数 对悬料影响的权值因子,实现对悬料的预测. 该类 方法应用灵活,具有较强的自学习能力和适应性; 但 其仅从统计学角度出发,挖掘数据之间存在的客观 规律,缺乏与高炉已知经验规律的结合,模型精度严 重依赖于悬料样本数据,由于悬料样本数据属于小 样本集,难以保证悬料样本的泛化性,致使预测模型 难以客观的完全反映悬料本质规律,预测效果不佳. 高炉冶炼过程复杂,存在众多干扰,造成检测信 息的客观不确定性,而且悬料和参数之间的异常状 态之间没有确定的对应关系,也存在很大的客观不 确定性. 同时,专家系统由于自适应不足,而且模糊 隶属函数的设定通常基于专家经验,具有很大主观 性,致使推理结果存在主观不确定性; 智能预测方法 由于样本泛化性不足导致预测结果存在客观不确 定性. 针对悬料与参数之间,以及专家系统和智能预 测方法存在主、客观不确定性的问题,本文结合 D--S 证据理论,提出一种基于主客观证据融合的高炉悬 料预测方法. 首先,结合高炉生产过程分析悬料成 因,并综合现场生产经验和数据分析,确定悬料预测 输入参数. 在此基础上,一方面基于专家经验,以专 家规则的形式描述征兆参数与悬料的关系,通过模 糊专家推理,提取主观证据对悬料的可信度; 另一方 面,基于后验概率最小二乘支持向量机 ( least squares support vector machine,LSSVM) 模型,利用 LSSVM 的自学习能力,从统计学角度挖掘征兆参数 与悬料间存在的客观联系,提取客观证据对悬料的 可信度. 最后,基于 D--S 证据理论对主、客观证据 可信度进行融合推理,实现对悬料的预测. 1 高炉生产过程与悬料成因分析 高炉生产过程伴随着极其复杂的物理、化学反 应,目的是将铁矿石冶炼成合格的生铁[9]. 即炉料 在下降的过程中,与上升的煤气流发生氧化还原反 应,生成液态生铁. 悬料是指由于炉内料柱透气性 与煤气流运动极不适应而导致炉料停止下降的炉况 失常现象[10],是高炉冶炼生产过程中时常发生的异 常炉况,严重影响高炉冶炼生产的正常进行. 1. 1 高炉生产过程 高炉冶炼是将铁矿石还原成生铁的连续生产过 程. 烧结矿、焦炭、石灰石等固体原料按规定配比通 过高炉炉顶布料装置投入炉内形成炉料,炉料在重 力作用下向下运动. 热风炉输出的热风经下部风口 鼓入,并在风口回旋区与喷吹燃料和风口焦炭燃烧 生成煤气; 煤气被炉缸与炉顶压力差推动经炉料间 隙向上运动,进入高炉顶部. 炉料在逐步下降的过 程中,受到逆流向上运动的高温煤气的还原作用,不 断地被加热和还原而成为熔融状态,并以液态形式 通过焦炭层的空隙向下滴落,最终在炉缸形成炉渣 和铁水. 1. 2 悬料成因分析 悬料是高炉炉热、炉凉和管道等异常炉况发展 到一定程度的必然结果. 从炼铁工艺角度分析,炉 内料柱透气性因炉温过热而变差,导致风压增大,致 使热风造成的压差增大而托起料柱; 同时在炉墙上 黏结的软熔炉料,阻止了边缘煤气流的上升和料柱 的下降,进而形成悬料; 从流体力学角度分析,炉内 煤气流向上运动的进程受阻,在高炉下部区域内形 成湍流并产生巨大的托力,从而阻碍炉料的下行,最 终形成悬料. 通过上述两个角度的分析,总结悬料形成的内 在机理如下: 炉内煤气流的通道突然被阻塞,造成料 柱某一层面的煤气压力快速升高,加大了对该层面 以上料柱的托力,随着托力增加导致该层面以上的 料柱停止运动,从而形成悬料,该层面以下的料柱正 常下降. 2 基于主客观证据融合的悬料预测方案 正常炉况时,高炉各项过程检测参数在正常范 围内小幅波动,变化趋势相适应( 如风量和风压的 变化趋势相一致等) . 当即将发生异常炉况时,部分 检测参数会有明显的异常波动,经过现场专家多年 的总结分析,悬料发生前的主要征兆有: ( 1) 风压逐 渐上升,风量缓慢减少; ( 2) 料速持续减小,料尺越 来越宽; ( 3) 料柱透气性变差,炉顶压力降低,全压 差升高. 结合专家经验,通常选取风量、风压、透气性指 ·507·
·508 北京科技大学学报 第36卷 数、顶压、压差和料速等参数作为悬料的征兆参 LSSVM模型,提取客观证据可信度. 数叫.为了综合利用参数的变化趋势及快慢程度 3.1基于模糊专家推理的主观证据提取 信息,增加对参数变化的敏感性,本文选取风量偏 高炉炉况预测专家系统中,采用模糊理论的隶 差、风压梯度、透气性指数偏差、压差梯度、顶压梯度 属度函数,能够有效表示征兆参数和炉况变化的模 和料速均值(即平均料速)六个参数,作为悬料预测 糊状态,实现炉况预测的不确定性推理.首先, 的输入参数. 利用征兆参数的特征值构造隶属度函数,实现参数 在DS证据理论中,确定辨识框架中命题的基 的模糊化处理:然后,基于产生式规则进行模糊专家 本可信度是关键☒,本文拓展了D-S证据理论中 推理,计算预测规则对炉况的支持度,提取主观证据 基本可信度的内涵,不局限于传感器的可信度,也包 可信度 括专家规则的推理结果可信度和后验概率LSSVM 通过分析悬料征兆参数的变化特点,其征兆参 模型的预测结果可信度.同时,针对专家系统推理 数特征值的隶属度分别采用升半梯形(x)和降半 结果存在主观不确定性的特点,本文称其为主观证 梯形2(x)隶属度函数计算. 据:针对智能预测方法的预测结果存在客观不确定 0 x:≤aai 性的特点,本文称其为客观证据.基于主客观证据 xi-da 融合的高炉悬料预测方案包括以下几个方面:首先, u1(x)= d2-da aaa2 结合输入参数的模糊化处理,进行模糊推理,提取主 式中:x,和x2分别表示风压梯度和压差梯度;aa和 观证据可信度:其次,通过悬料训练样本确定LSSVM a2为常数,表示对应参数的函数转折点. 炉况分类器和对应后验概率模型的相关参数,进而 x:≤b1; 对测试样本进行处理,提取客观证据可信度;最后, 基于DS证据融合理论,对主、客观证据可信度进 2(x: b2 -ba bab2 判断,实现悬料预测,从而综合利用专家知识和 (2) LSSVM的自学习能力,降低不确定性,提高预测精 式中:x1、x2、x和x分别表示风量偏差、透气性指数 度.该方案的总体结构框图如图1所示 偏差、顶压梯度和料速均值;b和b2为常数,表示对 客观诊断 应参数的函数转折点. 后验概率 客观证据 输人参数 LSSVM模型可信度D-S 根据正常炉况征兆参数的变化特点,其征兆参 炉过程检测参 证 悬料预测结果 数特征值的隶属度采用梯形山3(x)隶属度函数 主观诊断 融合 计算. 输人参数 主观证据 模糊专家推理 x,≤ci 、专家知识 可信度 xi-ca 图1悬料预测总体结构框图 CaC4 主客观证据提取是本文提出的悬料预测方法的 (3) 关键环节,也是主客观证据融合的基础.首先,根据 式中:x,(i=1,2,…,6)分别表示风量偏差、风压梯 悬料相关的关键专家规则,基于模糊专家推理提取 度、透气性指数偏差、压差梯度、顶压梯度和料速均 主观证据可信度:其次,针对神经网络存在泛化能力 值,cac2c和ca为常数,表示对应参数的函数转 差、易陷入局部最优解等缺陷,以及标准支持向量机 折点 (support vector machine,SVM)存在只能给出确定的 本文采用基于产生式的不精确知识表示法,结 样本类别标号,无法给出对每个类别隶属度(即样 合可信度表示专家知识.规则的一般形式如下 本的后验概率)的问题,本文通过建立后验概率 所示:
北 京 科 技 大 学 学 报 第 36 卷 数、顶 压、压差和料速等参数作为悬料的征兆参 数[11]. 为了综合利用参数的变化趋势及快慢程度 信息,增加对参数变化的敏感性,本文选取风量偏 差、风压梯度、透气性指数偏差、压差梯度、顶压梯度 和料速均值( 即平均料速) 六个参数,作为悬料预测 的输入参数. 在 D--S 证据理论中,确定辨识框架中命题的基 本可信度是关键[12],本文拓展了 D--S 证据理论中 基本可信度的内涵,不局限于传感器的可信度,也包 括专家规则的推理结果可信度和后验概率 LSSVM 模型的预测结果可信度. 同时,针对专家系统推理 结果存在主观不确定性的特点,本文称其为主观证 据; 针对智能预测方法的预测结果存在客观不确定 性的特点,本文称其为客观证据. 基于主客观证据 融合的高炉悬料预测方案包括以下几个方面: 首先, 利用产生式规则描述专家知识,总结悬料预测规则, 结合输入参数的模糊化处理,进行模糊推理,提取主 观证据可信度; 其次,通过悬料训练样本确定LSSVM 炉况分类器和对应后验概率模型的相关参数,进而 对测试样本进行处理,提取客观证据可信度; 最后, 基于 D--S 证据融合理论,对主、客观证据可信度进 行融合推理,并根据判决规则对预测结果进行辅助 判 断,实 现 悬 料 预 测,从而综合利用专家知识和 LSSVM 的自学习能力,降低不确定性,提高预测精 度. 该方案的总体结构框图如图 1 所示. 图 1 悬料预测总体结构框图 Fig. 1 Overall structure diagram of hanging prediction 3 基于模糊专家推理和后验概率 LSSVM 模型的主客观证据提取 主客观证据提取是本文提出的悬料预测方法的 关键环节,也是主客观证据融合的基础. 首先,根据 悬料相关的关键专家规则,基于模糊专家推理提取 主观证据可信度; 其次,针对神经网络存在泛化能力 差、易陷入局部最优解等缺陷,以及标准支持向量机 ( support vector machine,SVM) 存在只能给出确定的 样本类别标号,无法给出对每个类别隶属度( 即样 本的后验 概 率) 的 问 题,本文通过建立后验概率 LSSVM 模型,提取客观证据可信度. 3. 1 基于模糊专家推理的主观证据提取 高炉炉况预测专家系统中,采用模糊理论的隶 属度函数,能够有效表示征兆参数和炉况变化的模 糊状态,实现炉况预测的不确定性推理[13]. 首先, 利用征兆参数的特征值构造隶属度函数,实现参数 的模糊化处理; 然后,基于产生式规则进行模糊专家 推理,计算预测规则对炉况的支持度,提取主观证据 可信度. 通过分析悬料征兆参数的变化特点,其征兆参 数特征值的隶属度分别采用升半梯形 μ1 ( x) 和降半 梯形 μ2 ( x) 隶属度函数计算. μ1 ( xi ) = 0 xi≤ai1 ; xi - ai1 ai2 - ai1 ai1 < xi≤ai2 ; 1 xi > ai2 . i = 1,2 ( 1) 式中: x1和 x2 分别表示风压梯度和压差梯度; ai1 和 ai2为常数,表示对应参数的函数转折点. μ2 ( xi ) = 1 xi≤bi1 ; bi2 - xi bi2 - bi1 bi1 < xi≤bi2 ; 0 xi > bi2 . i = 1,2,3,4. ( 2) 式中: x1、x2、x3和 x4分别表示风量偏差、透气性指数 偏差、顶压梯度和料速均值; bi1和 bi2为常数,表示对 应参数的函数转折点. 根据正常炉况征兆参数的变化特点,其征兆参 数特征值的隶属度采用梯形 μ3 ( x ) 隶 属 度 函 数 计算. μ3 ( xi ) = 0 xi≤ci1 ; xi - ci1 ci2 - ci1 ci1 < xi≤ci2 ; 1 ci2 < xi≤ci3 ; ci4 - xi ci4 - ci3 ci3 < xi≤ci4 ; 0 xi > ci4 . i = 1,2,…,6. ( 3) 式中: xi ( i = 1,2,…,6) 分别表示风量偏差、风压梯 度、透气性指数偏差、压差梯度、顶压梯度和料速均 值,ci1、ci2、ci3和 ci4 为常数,表示对应参数的函数转 折点. 本文采用基于产生式的不精确知识表示法,结 合可信度表示专家知识. 规则的一般形式如下 所示: ·508·
第4期 曹卫华等:基于主客观证据融合的高炉悬料预测方法 ·509· IF E THEN H (CF (H,E)). f(x)=w(x)+b (6) 其中,E为规则前提,H为规则结论;CF(H,E)为规 为LSSVM炉况预测分类器中样本x的标准输出值. 则的可信度,反映了前提条件与结论的联系强度,其 在构造最优线性决策函数时,根据结构风险最 作用域为D,1],由专家经验确定.本文依据主观上 小化原则,同时引入松弛因子:来解决线性不可分 反映悬料本质特征的关键专家知识,定义如下悬 问题,建立如下优化目标函数: 料相关预报规则: (7) F风量降低AND风压升高THEN风量情况 mi(w,=子w2+c∑贫, (CF=0.9): s.t.y:w'中(x)+b]=1-,i=1,2,…,N. F透气性指数降低OR压差升高OR顶压降低 (8) THEN透气性情况(CF=0.9): 式中,C为惩罚因子,实现对错分样本的惩罚. F平均料速降低THEN料速情况(CF=I). 采用拉格朗日法将上述优化问题转换成求解如 对照上述悬料相关预报规则,本文定义正常炉 下线性方程组: 况相关的预报规则如下所示: 「0 F风量波动正常AND风压波动正常THEN风 yK(x1x)+1/C… yyxK(xx) 量情况(CF=0.8): F透气性指数波动正常AND压差波动正常 yxYK(N,) yNYAK(XNXN)+1/G AND顶压波动正常THEN透气性情况(CF=O.9): F平均料速正常THEN料速情况(CF=I). 基于上述预测规则,进行模糊专家推理,提取风 (9) 量、透气性和料速情况对炉况的支持度,即为主观证 据可信度. 采用最小二乘法求出回归系数&:和偏差b,则 3.2基于后验概率LSSVM模型的客观证据提取 N LSSVM是标准支持向量机的一种扩展的,而 f(x)=∑ay:K(xx)+b. (10) 后验概率LSSVM模型,是建立在标准LSSVM基础 其中,K(x,x;)=中(x)中(x)是核函数,根据高炉炉 上的,通过后验概率模型将其输出值映射为概率值, 况分类特点,选取高斯径向基函数作为LSSVM炉况 提高了LSSVM的适应性.在本文中,基于后验概率 预测分类器的核函数,即 LSSVM模型的客观证据提取包括两个步骤:首先建 K(x,x)=exp x-x; (11) 立面向悬料和正常炉况二类分类问题的LSSVM炉 况预测分类器:其次通过LSSVM的后验概率模型将 式中,σ为核宽度. 分类器的标准输出映射为概率值,提取客观证据可 3.2.2后验概率LSSVM模型 信度. 含两个参数M和T的Sigmoid函数形式灵活, 3.2.1 LSSVM炉况预测分类器 而且以其作为输出概率模型的后验概率SVM模型, 将风量偏差、风压梯度、透气性指数偏差、压差 在实际应用中效果良好可,因此本文采用该函数作 梯度、顶压梯度和料速均值六个参数作为LSSVM炉 为LSSVM炉况分类器的后验概率模型.其概率输 况预测分类器的输入向量,记为向量x·将对应高 出形式为 炉炉况预测结果作为炉况预测分类器的输出,记为 (12) y:y:∈{-1,1}为类别标签,1代表悬料,-1代表 p(y=1f)=1+p(M听x)+T' 正常炉况.由输入、输出构成炉况分类器的样本集 p(=-1lf(x))=1-p(y=1f(x).(13) S={(xy)li=1,2,…,N,x:eR(n=6), 式中,f(x)为LSSVM炉况预测分类器的标准输出, y∈{-1,1}}, (4) p(y=1f(x))和p(y=-1f(x))分别表示在输出 则LSSVM炉况预测分类器通过非线性映射函数 值为f(x)的条件下样本x属于1类(悬料)和-1 中(x)将输入空间样本映射到高维特征空间,并在此 类(正常炉况)的概率,参数M和T控制Sigmoid函 数的形态,其优化策略可通过如下最大似然问题 空间构造最优线性决策函数,即分类函数 y(x)=sgn [f(x)] (5) 解决: min F(z)= 其中 =(M.T)T
第 4 期 曹卫华等: 基于主客观证据融合的高炉悬料预测方法 IF E THEN H ( CF( H,E) ) . 其中,E 为规则前提,H 为规则结论; CF( H,E) 为规 则的可信度,反映了前提条件与结论的联系强度,其 作用域为[0,1],由专家经验确定. 本文依据主观上 反映悬料本质特征的关键专家知识[14],定义如下悬 料相关预报规则: IF 风量降低 AND 风压升高 THEN 风量情况 ( CF = 0. 9) ; IF 透气性指数降低 OR 压差升高 OR 顶压降低 THEN 透气性情况( CF = 0. 9) ; IF 平均料速降低 THEN 料速情况( CF = 1) . 对照上述悬料相关预报规则,本文定义正常炉 况相关的预报规则如下所示: IF 风量波动正常 AND 风压波动正常 THEN 风 量情况( CF = 0. 8) ; IF 透气性指数波动正常 AND 压差波动正常 AND 顶压波动正常 THEN 透气性情况( CF = 0. 9) ; IF 平均料速正常 THEN 料速情况( CF = 1) . 基于上述预测规则,进行模糊专家推理,提取风 量、透气性和料速情况对炉况的支持度,即为主观证 据可信度. 3. 2 基于后验概率 LSSVM 模型的客观证据提取 LSSVM 是标准支持向量机的一种扩展[15],而 后验概率 LSSVM 模型,是建立在标准 LSSVM 基础 上的,通过后验概率模型将其输出值映射为概率值, 提高了 LSSVM 的适应性. 在本文中,基于后验概率 LSSVM 模型的客观证据提取包括两个步骤: 首先建 立面向悬料和正常炉况二类分类问题的 LSSVM 炉 况预测分类器; 其次通过 LSSVM 的后验概率模型将 分类器的标准输出映射为概率值,提取客观证据可 信度. 3. 2. 1 LSSVM 炉况预测分类器 将风量偏差、风压梯度、透气性指数偏差、压差 梯度、顶压梯度和料速均值六个参数作为 LSSVM 炉 况预测分类器的输入向量,记为向量 xi . 将对应高 炉炉况预测结果作为炉况预测分类器的输出,记为 yi,yi∈{ - 1,1} 为类别标签,1 代表悬料,- 1 代表 正常炉况. 由输入、输出构成炉况分类器的样本集 S = { ( xi,yi ) | i = 1,2,…,N,xi∈Rn ( n = 6) , yi∈{ - 1,1} } , ( 4) 则 LSSVM 炉况预测分类器通过非线性映射函数 ( x) 将输入空间样本映射到高维特征空间,并在此 空间构造最优线性决策函数,即分类函数 y( x) = sgn[f( x) ], ( 5) 其中 f( x) = wT ( x) + b ( 6) 为 LSSVM 炉况预测分类器中样本 x 的标准输出值. 在构造最优线性决策函数时,根据结构风险最 小化原则,同时引入松弛因子 ξi 来解决线性不可分 问题,建立如下优化目标函数: min ω,b,ξ J( w,ξ) = 1 2 ‖w‖2 + C ∑ N i = 1 ξ 2 i , ( 7) s. t. yi [wT ·( x) + b]= 1 - ξi,i = 1,2,…,N. ( 8) 式中,C 为惩罚因子,实现对错分样本的惩罚. 采用拉格朗日法将上述优化问题转换成求解如 下线性方程组[16]: 0 y1 … yN y1 y1 y1K( x1,x1 ) + 1 /C … y1 yN K( x1,xN ) yN yN y1K( xN,x1 ) … yN yN K( xN,xN ) + 1 /C · b α1 α N = 0 1 1 . ( 9) 采用最小二乘法求出回归系数 αi 和偏差 b,则 f( x) = ∑ N i = 1 αiyiK( x,xi ) + b. ( 10) 其中,K( x,xi ) = ( x) ( xi ) 是核函数,根据高炉炉 况分类特点,选取高斯径向基函数作为 LSSVM 炉况 预测分类器的核函数,即 K( x,xi ) = ( exp - ‖x - xi‖2 σ2 ) . ( 11) 式中,σ 为核宽度. 3. 2. 2 后验概率 LSSVM 模型 含两个参数 M 和 T 的 Sigmoid 函数形式灵活, 而且以其作为输出概率模型的后验概率 SVM 模型, 在实际应用中效果良好[17],因此本文采用该函数作 为 LSSVM 炉况分类器的后验概率模型. 其概率输 出形式为 p( y = 1 | f( x) ) = 1 1 + exp( Mf( x) + T) , ( 12) p( y = - 1 | f( x) ) = 1 - p( y = 1 | f( x) ) . ( 13) 式中,f( x) 为 LSSVM 炉况预测分类器的标准输出, p( y = 1 | f( x) ) 和 p( y = - 1 | f( x) ) 分别表示在输出 值为 f( x) 的条件下样本 x 属于 1 类( 悬料) 和 - 1 类( 正常炉况) 的概率,参数 M 和 T 控制 Sigmoid 函 数的形态,其优化策略可通过如下最大似然问题 解决: min z = ( M,T) T F( z) = ·509·
·510 北京科技大学学报 第36卷 ,log(p:)+(1-t,)log(1-p)].(14) 式中,K= ∑Πm,(,)E1, (18) 对式(12)求解参数M和T,确定LSSVM后验概率 模型,获得输出概率,实现客观证据可信度的提取 lBel(A)-Bel(A2)>e2 则A为炉况预测结果,即当某种炉况具有最大的基 4基于D一S证据理论的主客观证据融合 本概率分配值,而且基本概率分配值大于61,同时 方法 其信任函数值与另一种炉况的信任函数值相差大于 基于主客观证据,如何进行有效的融合,是本文 82时,该种炉况即为预测结果.参照文献21]关于 提出的悬料预测方法的核心问题.DS证据理论作 DS证据理论决策准则中相关阈值的取值范围:&1 为一种理论完善的不确定推理方法,在多源信息融 不小于0.6,2不小于0.15.本文结合现场工业实 合阿和故障诊断四等方面得到了广泛应用。因此, 际情况,选取61=0.6,62=0.2. 基于DS证据理论,完成主客观证据融合,实现悬 5仿真结果及分析 料预测,是一个良好的途径 在DS证据理论中,由互不相容的基本命题组 针对某钢铁企业2200m3高炉的实际运行情况, 成的完备集合称为辨识框架,常用⊙表示.它是对 采用五点三次平滑算法对现场采集的数据进行平滑 某一问题的所有可能答案的一个集合,但只有一个 处理,消除噪声等干扰,然后分别选取不同生产条件 答案是正确的.本文定义辨识框架⊙={A,B},A和 下专家认定的炉况状态和相应的征兆数据,包括正 B分别为悬料和正常炉况.通过模糊专家推理提取 常炉况300组,悬料炉况150组.对征兆参数进行 风量情况、透气性情况和料速情况三条规则对炉况 特征提取,获得风量偏差、风压梯度、透气性指数偏 的支持度,即主观据可信度P(A)和P,(B),其中 差、压差梯度、顶压梯度及料速均值六个特征参数, i=1,2,3;通过后验概率LSSVM模型提取客观证据 并定义炉况状态的类别标签(即1代表悬料,-1代 可信度P,(A)和P,(B).则基于DS证据理论的主 表正常炉况),共同构成炉况样本数据.随机选取悬 客观证据融合,主要包括以下步骤: 料和正常炉况各60%样本数据作为训练样本,建立 步骤1对主、客观证据可信度进行归一化处 后验概率LSSVM炉况预测模型;余下的40%作为 理,获得炉况的基本概率分配函数m:(),又称为炉 测试样本,对本文提出的预测方法进行仿真验证 况的基本可信度,表示主、客观证据对相应炉况的信 首先,在训练LSSVM炉况预测分类器的过程 任程度, 中,采用交叉验证法确定LSSVM炉况预测分类器的 P:(a) 惩罚因子和核参数,C=8,σ=0.4:其次,将训练样 m,(90)=pW+P(®0c8,i=1,2,…4 本分为3份,其中任意2份用于训练C=8,σ=0.4 (15) 的LSSVM炉况预测分类器,余下1份用于拟合后验 步骤2根据DS证据理论中基于正交和的合 概率模型.重复进行三次训练和拟合,将M和T参 成法则,计算融合后的炉况基本概率分配函数 数三次求解结果的均值作为参数M和T终值,M= m(),表示融合后的证据对相应炉况的信任程度, -4.572,T=0.315 通过上述后验概率LSSVM模型,利用测试样 ∑Πm,(,) n9=i=1 本,获得单一基于D-S证据理论的主观证据融合炉 m(a)= -,H0C⊙,0≠⑦,(16) 1-K 况预测结果(即表1中方法I),以及单一后验概率
北 京 科 技 大 学 学 报 第 36 卷 - ∑ N i = 1 [ti log( pi ) + ( 1 - ti ) log( 1 - pi ) ]. ( 14) 式中, pi = 1 1 + exp( Mf( x) + T) , ti = N + + 1 N + + 2 yi = 1, 1 N - + 2 y { i = - 1, i = 1,2,…,N, N + 和 N - 分别是 y = 1 类和 y = - 1 类的样本数. 采用具有逆向线性搜索特点的牛顿迭代法[18], 对式( 12) 求解参数 M 和 T,确定 LSSVM 后验概率 模型,获得输出概率,实现客观证据可信度的提取. 4 基于 D--S 证据理论的主客观证据融合 方法 基于主客观证据,如何进行有效的融合,是本文 提出的悬料预测方法的核心问题. D--S 证据理论作 为一种理论完善的不确定推理方法,在多源信息融 合[19]和故障诊断[20]等方面得到了广泛应用. 因此, 基于 D--S 证据理论,完成主客观证据融合,实现悬 料预测,是一个良好的途径. 在 D--S 证据理论中,由互不相容的基本命题组 成的完备集合称为辨识框架,常用 Θ 表示. 它是对 某一问题的所有可能答案的一个集合,但只有一个 答案是正确的. 本文定义辨识框架 Θ = { A,B} ,A 和 B 分别为悬料和正常炉况. 通过模糊专家推理提取 风量情况、透气性情况和料速情况三条规则对炉况 的支持度,即主观据可信度 Pi ( A) 和 Pi ( B) ,其中 i = 1,2,3; 通过后验概率 LSSVM 模型提取客观证据 可信度 P4 ( A) 和 P4 ( B) . 则基于 D--S 证据理论的主 客观证据融合,主要包括以下步骤: 步骤 1 对主、客观证据可信度进行归一化处 理,获得炉况的基本概率分配函数 mi ( θ) ,又称为炉 况的基本可信度,表示主、客观证据对相应炉况的信 任程度, mi ( θ) = Pi ( θ) Pi ( A) + Pi ( B) ,θΘ,i = 1,2,…,4 ( 15) 步骤 2 根据 D--S 证据理论中基于正交和的合 成法 则,计算融合后的炉况基本概率分配函数 m( θ) ,表示融合后的证据对相应炉况的信任程度, m( θ) = ∩ ∑θj = θ ∏ 4 i = 1 mi ( θj ) 1 - K ,θΘ,θ≠,( 16) 式中,K = ∩ ∑θj = ∏ 4 i = 1 mi ( θj ) < 1; 步骤 3 根据融合后的基本概率分配函数,计 算炉况信任函数 Bel( θ) ,表示证据对相应炉况为真 的信任程度, Bel( θ) = m( θ) ,θΘ; ( 17) 步骤 4 根据判决规则,确定炉况预测结果,并 输出炉况的信任函数值,即为相应炉况的发生概率, 其判决规则如下为 m( A1 ) = max{ m( Ai ) ,AiΘ,i = 1,2} , m( A1 ) > ε1, Bel( A1 ) - Bel( A2 ) > ε2 { . ( 18) 则 A1为炉况预测结果,即当某种炉况具有最大的基 本概率分配值,而且基本概率分配值大于 ε1,同时 其信任函数值与另一种炉况的信任函数值相差大于 ε2 时,该种炉况即为预测结果. 参照文献[21]关于 D--S 证据理论决策准则中相关阈值的取值范围: ε1 不小于 0. 6,ε2 不小于 0. 15. 本文结合现场工业实 际情况,选取 ε1 = 0. 6,ε2 = 0. 2. 5 仿真结果及分析 针对某钢铁企业 2200 m3 高炉的实际运行情况, 采用五点三次平滑算法对现场采集的数据进行平滑 处理,消除噪声等干扰,然后分别选取不同生产条件 下专家认定的炉况状态和相应的征兆数据,包括正 常炉况 300 组,悬料炉况 150 组. 对征兆参数进行 特征提取,获得风量偏差、风压梯度、透气性指数偏 差、压差梯度、顶压梯度及料速均值六个特征参数, 并定义炉况状态的类别标签( 即 1 代表悬料,- 1 代 表正常炉况) ,共同构成炉况样本数据. 随机选取悬 料和正常炉况各 60% 样本数据作为训练样本,建立 后验概率 LSSVM 炉况预测模型; 余下的 40% 作为 测试样本,对本文提出的预测方法进行仿真验证. 首先,在训练 LSSVM 炉况预测分类器的过程 中,采用交叉验证法确定 LSSVM 炉况预测分类器的 惩罚因子和核参数,C = 8,σ = 0. 4; 其次,将训练样 本分为 3 份,其中任意 2 份用于训练 C = 8,σ = 0. 4 的 LSSVM 炉况预测分类器,余下 1 份用于拟合后验 概率模型. 重复进行三次训练和拟合,将 M 和 T 参 数三次求解结果的均值作为参数 M 和 T 终值,M = - 4. 572,T = 0. 315. 通过上述后验概率 LSSVM 模型,利用测试样 本,获得单一基于 D--S 证据理论的主观证据融合炉 况预测结果( 即表 1 中方法Ⅰ) ,以及单一后验概率 ·510·
第4期 曹卫华等:基于主客观证据融合的高炉悬料预测方法 ·511· LSSVM模型(即表中方法Ⅱ)和本文所提预测方法 方法同样具有很大的优势,而三种预测方法在不同 的炉况预测结果,炉况预测概率曲线如图2和图3 治炼状态下的预测准确率有所不同,是由于不同治 所示,炉况预测结果准确率如表1所示.其中,单一 炼状态下炉况样本参数幅值及其波动大小的不同所 主观证据的融合结果和本文所提的预测方法,根据 引起的 式(18)来确定炉况预测结果:单一后验概率LSSVM 表2治炼状态划分准则 模型的输出以0.6为决策值,判断炉况预测结果. Table 2 Smelting state division criterion 表1炉况预测结果准确率对比 治炼状态 治炼强度,1/(m3·d1) Table 1 Comparison of prediction accuracy for furnace condition 强化 1>1.0 炉况测试 炉况预测 炉况预测 常规 0.9≤1≤1.0 预测 样本数 结果数 准确率/% 弱化 1<0.9 方法 悬料 正常 悬料 正常 悬料正常 12 方法I 60 120 51 100 85.00 83.33 1.0 方法Ⅱ 60 120 53 104 88.33 86.67 0.8 本文方法60 120 56 111 93.3392.50 0.6 1.2 实际发生概率 0.4 实际发生概率 。-本文所提方法须测概岸 1.0 0.2 。后验概率LSSVM模型预测概率 0.8 。基于D-S的主观证据融合预测概率 10 19 0.6 样本序列 0.4 ©本文所提方法预测概率 图4强化冶炼状态下悬料预测概率曲线 0.2 。后验概率1SSVM模型预测概率 →基于D-S的主观证据险合演测概率 Fig.4 Predicted probubility curves of hanging under intensity smel- 05101520253035404550560 ting 样本序列 图2悬料侦测概率曲线 实际发生概率 Fig.2 Predicted probability curves of hanging 1.2 →基于D-S的主观证据融合预测概率 1.0 ·本文所提方法预测概 0.8 0.2 。一后验概率LSSVM模型顶测做率 基于D-S的主观证据融合预测概率 0.6 10 152025303537 0.4 样本序列 一实际发生概率 0.2 ·本文所提方法预测概率 图5强化治炼状态下正常炉况预测概率曲线 一后验概率LSSVM模型预测概率 Fig.5 Predicted probability curves of normal furnace condition under 20 40 60 80 100 120 intensity smelting 样本序列 图3正常炉况预测概率曲线 1.2 Fig.3 Predicted probability curves of normal fumace condition 1.0 根据该钢铁企业高炉生产现场的操作制度,依 据冶炼强度I(t·(m3·d)-1)将高炉治炼过程划分 0.6 实际发生概率 为不同的治炼状态,具体如表2所示.按表2所述 4 。本文所提方法面测概率 治炼状态划分标准,分析各炉况测试样本所处的治 02 一后验概率LSSVM模型预测概幸 基于D-S的主观证据融合预测概率 炼状态,并统计各冶炼状态所含炉况测试样本数及 10 15 20 23 相应的炉况预测结果.相应的统计结果如图4~图 样本序列 9和表3~表5所示. 图6常规冶炼状态下悬料预测概率曲线 从表1及表3~表5可以看出,在整体和不同 Fig.6 Predicted probability curves of hanging under conventional 治炼状态下,本文所提预测方法相比两种单一预测 smelting
第 4 期 曹卫华等: 基于主客观证据融合的高炉悬料预测方法 LSSVM 模型( 即表中方法Ⅱ) 和本文所提预测方法 的炉况预测结果,炉况预测概率曲线如图 2 和图 3 所示,炉况预测结果准确率如表 1 所示. 其中,单一 主观证据的融合结果和本文所提的预测方法,根据 式( 18) 来确定炉况预测结果; 单一后验概率 LSSVM 模型的输出以 0. 6 为决策值,判断炉况预测结果. 表 1 炉况预测结果准确率对比 Table 1 Comparison of prediction accuracy for furnace condition 预测 方法 炉况测试 样本数 炉况预测 结果数 炉况预测 准确率/% 悬料 正常 悬料 正常 悬料 正常 方法Ⅰ 60 120 51 100 85. 00 83. 33 方法Ⅱ 60 120 53 104 88. 33 86. 67 本文方法 60 120 56 111 93. 33 92. 50 图 2 悬料预测概率曲线 Fig. 2 Predicted probability curves of hanging 图 3 正常炉况预测概率曲线 Fig. 3 Predicted probability curves of normal furnace condition 根据该钢铁企业高炉生产现场的操作制度,依 据冶炼强度 I ( t·( m3 ·d) - 1 ) 将高炉冶炼过程划分 为不同的冶炼状态,具体如表 2 所示. 按表 2 所述 冶炼状态划分标准,分析各炉况测试样本所处的冶 炼状态,并统计各冶炼状态所含炉况测试样本数及 相应的炉况预测结果. 相应的统计结果如图 4 ~ 图 9 和表 3 ~ 表 5 所示. 从表 1 及表 3 ~ 表 5 可以看出,在整体和不同 冶炼状态下,本文所提预测方法相比两种单一预测 方法同样具有很大的优势,而三种预测方法在不同 冶炼状态下的预测准确率有所不同,是由于不同冶 炼状态下炉况样本参数幅值及其波动大小的不同所 引起的. 表 2 冶炼状态划分准则 Table 2 Smelting state division criterion 冶炼状态 冶炼强度,I/( t·m - 3 ·d - 1 ) 强化 I > 1. 0 常规 0. 9≤ I ≤ 1. 0 弱化 I < 0. 9 图 4 强化冶炼状态下悬料预测概率曲线 Fig. 4 Predicted probability curves of hanging under intensity smelting 图 5 强化冶炼状态下正常炉况预测概率曲线 Fig. 5 Predicted probability curves of normal furnace condition under intensity smelting 图 6 常规冶炼状态下悬料预测概率曲线 Fig. 6 Predicted probability curves of hanging under conventional smelting ·511·
·512 北京科技大学学报 第36卷 1.2 表4常规治炼状态下炉况预测结果准确率对比 Table 4 Comparison of prediction accuracy for fumace condition under 1.0 conventional smelting 炉况测试 炉况预测 炉况预测 预测 样本数 结果数 准确率/% 一实际发生概率 0.4 方法 。本文所提方法预测概率 悬料正常 悬料正常 悬料正常 02 。-后验概率LSSVM模型预测概率 基于D-S的主观证据融合预测概率 方法I 25 52 22 44 88.0084.61 101520253035404550 方法Ⅱ 25 52 22 45 88.00 86.54 样本序列 本文方法 25 52 24 49 96.0094.23 图7常规治炼状态下正常炉况预测概率曲线 Fig.7 Predieted probability curves of normal furnace condition under 表5弱化治炼状态下炉况预测结果准确率对比 conventional smelting Table 5 Comparison of prediction accuracy for fumace condition under smelting weakening 1.2 炉况测试 炉况预测 炉况预测 1.0 预测 样本数 结果数 准确率/% 方法 悬料正常 悬料正常 悬料正常 0.6 方法I 16 31 13 26 81.2583.87 实际发生概率 0.4 本文所提方法测概率 方法Ⅱ 16 31 14 27 87.5087.10 0.2 。-后验概率LSSVM模型预测概率 基于D-S的主观证据融合预测概率 本文方法 16 31 15 29 93.75 93.55 10 1516 样本序列 在强化治炼状态下,常采用高压操作,喷吹燃料 图8弱化治炼状态下悬料预测概率曲线 和高风温等操作,使风量、风压、透气性指数、压差、 Fig.8 Predicted probability curves of hanging under smelting weake 顶压的幅值及其波动较大,同时缩短了治炼周期,使 ning 料速均值偏大.因此,在该状态下的本文炉况样本 参数幅值普遍偏大 在常规冶炼状态下,高炉顶压,风温和风量等参 1.0 数处于常规状态,使风量、风压、透气性指数、压差、 顶压的幅值及其波动处于中等程度,同时治炼周期 3 0.6 也处于常规状态,使料速均值处于中等程度.因此, 实际发生概率 。本文所提方法预测概率 在该状态下的本文炉况样本参数幅值普遍处于中等 0.2 后验概率LSSVM模型预测概率 程度. →基于D-S的主观证据陆合预测概率 10 152025 在弱化治炼状态下,常采用慢风和降低高炉顶 30 样木序列 压等操作,使风量、风压、透气性指数、压差、顶压的 图9弱化治炼状态下正常炉况预测概率曲线 幅值及其波动较小,同时延长了治炼周期,使料速均 Fig.9 Predicted probability curves of normal furnace condition under 值偏小.因此,在该状态下的本文炉况样本参数幅 smelting weakening 值普遍偏小 单一基于D-S证据理论的主观证据融合炉况 表3强化治炼状态下炉况预测结果准确率对比 预测方法中,炉况预报规则和隶属度函数的构造主 Table 3 Comparison of prediction accuracy for furnace condition under 要是基于常规治炼状态下的专家知识,因此该方法 intensity smelting 整体上预测准确率偏低,悬料和正常炉况的预测准 炉况测试 炉况预测 炉况预测 预测 确率分别为85.00%和83.33%,如表1所示:而在 样本数 结果数 准确率/% 方法 常规治炼状态下预测准确率较高,悬料和正常炉况 悬料正常 悬料 正常 悬料正常 的预测准确率分别为88.00%和84.61%,如表4所 方法I 19 37 16 30 84.2181.08 示.强化冶炼状态下,其悬料样本数据幅值相对于 方法Ⅱ 19 37 17 32 89.4786.49 常规治炼状态的悬料样本数据,除料速均值偏大外, 本文方法19 37 17 33 89.4789.19 其他参数与预报规则和隶属度函数符合程度较高
北 京 科 技 大 学 学 报 第 36 卷 图 7 常规冶炼状态下正常炉况预测概率曲线 Fig. 7 Predicted probability curves of normal furnace condition under conventional smelting 图 8 弱化冶炼状态下悬料预测概率曲线 Fig. 8 Predicted probability curves of hanging under smelting weakening 图 9 弱化冶炼状态下正常炉况预测概率曲线 Fig. 9 Predicted probability curves of normal furnace condition under smelting weakening 表 3 强化冶炼状态下炉况预测结果准确率对比 Table 3 Comparison of prediction accuracy for furnace condition under intensity smelting 预测 方法 炉况测试 样本数 炉况预测 结果数 炉况预测 准确率/% 悬料 正常 悬料 正常 悬料 正常 方法Ⅰ 19 37 16 30 84. 21 81. 08 方法Ⅱ 19 37 17 32 89. 47 86. 49 本文方法 19 37 17 33 89. 47 89. 19 表 4 常规冶炼状态下炉况预测结果准确率对比 Table 4 Comparison of prediction accuracy for furnace condition under conventional smelting 预测 方法 炉况测试 样本数 炉况预测 结果数 炉况预测 准确率/% 悬料 正常 悬料 正常 悬料 正常 方法Ⅰ 25 52 22 44 88. 00 84. 61 方法Ⅱ 25 52 22 45 88. 00 86. 54 本文方法 25 52 24 49 96. 00 94. 23 表 5 弱化冶炼状态下炉况预测结果准确率对比 Table 5 Comparison of prediction accuracy for furnace condition under smelting weakening 预测 方法 炉况测试 样本数 炉况预测 结果数 炉况预测 准确率/% 悬料 正常 悬料 正常 悬料 正常 方法Ⅰ 16 31 13 26 81. 25 83. 87 方法Ⅱ 16 31 14 27 87. 50 87. 10 本文方法 16 31 15 29 93. 75 93. 55 在强化冶炼状态下,常采用高压操作,喷吹燃料 和高风温等操作,使风量、风压、透气性指数、压差、 顶压的幅值及其波动较大,同时缩短了冶炼周期,使 料速均值偏大. 因此,在该状态下的本文炉况样本 参数幅值普遍偏大. 在常规冶炼状态下,高炉顶压,风温和风量等参 数处于常规状态,使风量、风压、透气性指数、压差、 顶压的幅值及其波动处于中等程度,同时冶炼周期 也处于常规状态,使料速均值处于中等程度. 因此, 在该状态下的本文炉况样本参数幅值普遍处于中等 程度. 在弱化冶炼状态下,常采用慢风和降低高炉顶 压等操作,使风量、风压、透气性指数、压差、顶压的 幅值及其波动较小,同时延长了冶炼周期,使料速均 值偏小. 因此,在该状态下的本文炉况样本参数幅 值普遍偏小. 单一基于 D--S 证据理论的主观证据融合炉况 预测方法中,炉况预报规则和隶属度函数的构造主 要是基于常规冶炼状态下的专家知识,因此该方法 整体上预测准确率偏低,悬料和正常炉况的预测准 确率分别为 85. 00% 和 83. 33% ,如表 1 所示; 而在 常规冶炼状态下预测准确率较高,悬料和正常炉况 的预测准确率分别为 88. 00% 和 84. 61% ,如表 4 所 示. 强化冶炼状态下,其悬料样本数据幅值相对于 常规冶炼状态的悬料样本数据,除料速均值偏大外, 其他参数与预报规则和隶属度函数符合程度较高, ·512·
第4期 曹卫华等:基于主客观证据融合的高炉悬料预测方法 ·513· 因此该方法的悬料预测准确率居中,为84.21%,如 提高了证据对炉况的支持度,降低了不确定性,保证 表3所示;而正常炉况样本数据幅值相对于常规治 了悬料预测的可靠性和可信度,有利于把握炉况的 炼状态的正常炉况样本数据普遍偏大,因此该方法 发展趋势. 的正常炉况预测准确率偏小,为81.08%,如表3所 表6主客观证据可信度 示.弱化治炼状态下,其悬料样本数据幅值相对于 Table 6 Credibility of subjective and objective evidences 常规治炼状态的悬料样本数据,除料速均值与预报 证据类型 悬料(A) 正常炉况(B) 规则和隶属度函数符合程度较高外,其他参数普遍 P,(A)=0.51 P1(B)=0.44 偏小,因此该方法的悬料预测准确率偏小,为 主观证据 P2(M)=0.56 P2(B)=0.42 81.25%,如表5所示:而正常炉况样本数据幅值相 P3(A)=0.66 P3(B)=0.39 对于常规冶炼状态的正常炉况样本,除料速均值偏 客观证据 P4(A)=0.77 P4(B)=0.23 小外,其他参数与预报规则和隶属度函数符合程度 较高,因此该方法的正常炉况预测准确率居中,为 表7证据融合结果 83.87%,如表5所示. Table 7 Fusion results of evidences 单一后验概率LSSVM模型预测方法中,炉况预 证据融合类型 悬料(A) 正常炉况(B) 测准确率与样本数据息息相关,不同治炼状态下样 主观证据融合 m,(A)=0.724 m,(B)=0.276 本数据数目从高到低依次为常规、强化和弱化治炼 主客观证据融合 m(A)=0.898 m(B)=0.102 状态.该方法相应的悬料预测准确率从高到低按治 炼状态依次为强化、常规和弱化冶炼状态,分别为 6 结论 89.47%、88.00%和87.50%,如表3~表5所示:该 (1)针对高炉关键异常炉况悬料难以预测的问 方法相应的正常炉况预测准确率从高到低依次按治 题,本文提出一种基于主客观证据融合的高炉悬料 炼状态为弱化、常规和强化治炼状态,分别为 87.10%、86.54%和86.49%,如表3~表5所示.总 预测方法.该方法基于D-S证据理论,恰当且有效 体上该方法在不同治炼状态下的预测准确率相差不 地综合利用了专家知识和LSSVM的自学习能力,为 高炉悬料预测提供一种新的思路 大,具有较强的自学习能力和适应性,并和样本数目 具有一致性.此外,由于炉况样本在不同治炼状态 (2)将本文所提预测方法与单一基于D-S证 下分布广泛,该方法的炉况预测准确率在整体及不 据理论的主观证据融合方法及单一后验概率LSSVM 同治炼状态中,普遍高于单一基于D-S证据理论的 模型进行仿真对比研究,仿真结果表明本文所提预 主观证据融合炉况预测方法. 测方法在悬料预测准确率和不确定性两方面的综合 本文所提基于主客观证据融合的高炉悬料预测 表现均优于两种单一预测方法. 方法,整体上对悬料和正常炉况的预报准确率分别 (3)对本文所提预测方法在不同治炼状态下的 为93.33%和92.50%,均高于两种单一预测方法, 悬料预测效果进行对比分析,结果表明本文所提预 测方法在不同治炼状态均具有较好的预测精度 如表1所示.同时,如表3表5所示,在不同冶炼 状态下,炉况样本数据存在显著差异的情况下,本文 参考文献 所提预测方法基于主、客观证据对炉况的支持度,综 [1]Wu S L,Tuo B Y,Zhang L.H,et al.Influence of coke reactivity 合了专家经验和LSSVM的自学习能力,均能够有效 on the ferric burden reduction of the lumpy zone in a blast furnace 预测炉况的发展趋势,并提高了悬料预测精度 J Univ Sci Technol Beijing,2013,35(3):282 针对某次常规治炼状态下悬料采集的样本数 (吴胜利,庹必阳,张丽华,等.焦炭反应性对高炉块状带含 据,本文所提预测方法提取的主、客观证据可信度如 铁炉料还原的影响.北京科技大学学报,2013,35(3):282) 2] Liu D X,Li X L,Zhou X,et al.Application of the fuzzy C-means 表6所示,基于DS证据理论的主观和主客观证据 clustering algorithm in blast furnace burden surface identification. 融合结果如表7所示. J Univ Sci Technal Beijing,2012,34(6):683 根据表7和式(18)可知,该次炉况预测结论为 (刘德馨,李晓理,周翔,等.模糊C均值聚类算法在高炉料 悬料,符合实际情况.从图2~图9,以及表6和表7 面分类中的应用.北京科技大学学报,2012,34(6):683) B] 中可以看出,本文所提预测方法的单次炉况预报概 Yi S,Xu Y M,MaZ W.An expert system for abnormal status di- agnosis on blast furnace.Metall Ind Autom,2002,26(1):15 率,相比两种单一预测方法更接近于实际发生概率。 (易粟,徐用懋,马竹捂.高炉异常炉况诊断专家系统.治金 同时,通过基于D-S证据理论的主客观证据融合, 自动化,2002,26(1):15)
第 4 期 曹卫华等: 基于主客观证据融合的高炉悬料预测方法 因此该方法的悬料预测准确率居中,为 84. 21% ,如 表 3 所示; 而正常炉况样本数据幅值相对于常规冶 炼状态的正常炉况样本数据普遍偏大,因此该方法 的正常炉况预测准确率偏小,为 81. 08% ,如表 3 所 示. 弱化冶炼状态下,其悬料样本数据幅值相对于 常规冶炼状态的悬料样本数据,除料速均值与预报 规则和隶属度函数符合程度较高外,其他参数普遍 偏小,因此该方法的悬料预测准确率偏小,为 81. 25% ,如表 5 所示; 而正常炉况样本数据幅值相 对于常规冶炼状态的正常炉况样本,除料速均值偏 小外,其他参数与预报规则和隶属度函数符合程度 较高,因此该方法的正常炉况预测准确率居中,为 83. 87% ,如表 5 所示. 单一后验概率 LSSVM 模型预测方法中,炉况预 测准确率与样本数据息息相关,不同冶炼状态下样 本数据数目从高到低依次为常规、强化和弱化冶炼 状态. 该方法相应的悬料预测准确率从高到低按冶 炼状态依次为强化、常规和弱化冶炼状态,分别为 89. 47% 、88. 00% 和 87. 50% ,如表 3 ~ 表 5 所示; 该 方法相应的正常炉况预测准确率从高到低依次按冶 炼状 态 为 弱 化、常规和强化冶炼状态,分 别 为 87. 10% 、86. 54% 和 86. 49% ,如表 3 ~ 表 5 所示. 总 体上该方法在不同冶炼状态下的预测准确率相差不 大,具有较强的自学习能力和适应性,并和样本数目 具有一致性. 此外,由于炉况样本在不同冶炼状态 下分布广泛,该方法的炉况预测准确率在整体及不 同冶炼状态中,普遍高于单一基于 D--S 证据理论的 主观证据融合炉况预测方法. 本文所提基于主客观证据融合的高炉悬料预测 方法,整体上对悬料和正常炉况的预报准确率分别 为 93. 33% 和 92. 50% ,均高于两种单一预测方法, 如表 1 所示. 同时,如表 3 ~ 表 5 所示,在不同冶炼 状态下,炉况样本数据存在显著差异的情况下,本文 所提预测方法基于主、客观证据对炉况的支持度,综 合了专家经验和 LSSVM 的自学习能力,均能够有效 预测炉况的发展趋势,并提高了悬料预测精度. 针对某次常规冶炼状态下悬料采集的样本数 据,本文所提预测方法提取的主、客观证据可信度如 表 6 所示,基于 D--S 证据理论的主观和主客观证据 融合结果如表 7 所示. 根据表 7 和式( 18) 可知,该次炉况预测结论为 悬料,符合实际情况. 从图 2 ~ 图 9,以及表 6 和表 7 中可以看出,本文所提预测方法的单次炉况预报概 率,相比两种单一预测方法更接近于实际发生概率. 同时,通过基于 D--S 证据理论的主客观证据融合, 提高了证据对炉况的支持度,降低了不确定性,保证 了悬料预测的可靠性和可信度,有利于把握炉况的 发展趋势. 表 6 主客观证据可信度 Table 6 Credibility of subjective and objective evidences 证据类型 悬料( A) 正常炉况( B) P1 ( A) = 0. 51 P1 ( B) = 0. 44 主观证据 P2 ( A) = 0. 56 P2 ( B) = 0. 42 P3 ( A) = 0. 66 P3 ( B) = 0. 39 客观证据 P4 ( A) = 0. 77 P4 ( B) = 0. 23 表 7 证据融合结果 Table 7 Fusion results of evidences 证据融合类型 悬料( A) 正常炉况( B) 主观证据融合 ms ( A) = 0. 724 ms ( B) = 0. 276 主客观证据融合 m( A) = 0. 898 m( B) = 0. 102 6 结论 ( 1) 针对高炉关键异常炉况悬料难以预测的问 题,本文提出一种基于主客观证据融合的高炉悬料 预测方法. 该方法基于 D--S 证据理论,恰当且有效 地综合利用了专家知识和 LSSVM 的自学习能力,为 高炉悬料预测提供一种新的思路. ( 2) 将本文所提预测方法与单一基于 D--S 证 据理论的主观证据融合方法及单一后验概率LSSVM 模型进行仿真对比研究,仿真结果表明本文所提预 测方法在悬料预测准确率和不确定性两方面的综合 表现均优于两种单一预测方法. ( 3) 对本文所提预测方法在不同冶炼状态下的 悬料预测效果进行对比分析,结果表明本文所提预 测方法在不同冶炼状态均具有较好的预测精度. 参 考 文 献 [1] Wu S L,Tuo B Y,Zhang L H,et al. Influence of coke reactivity on the ferric burden reduction of the lumpy zone in a blast furnace. J Univ Sci Technol Beijing,2013,35( 3) : 282 ( 吴胜利,庹必阳,张丽华,等. 焦炭反应性对高炉块状带含 铁炉料还原的影响. 北京科技大学学报,2013,35( 3) : 282) [2] Liu D X,Li X L,Zhou X,et al. Application of the fuzzy C-means clustering algorithm in blast furnace burden surface identification. J Univ Sci Technol Beijing,2012,34( 6) : 683 ( 刘德馨,李晓理,周翔,等. 模糊 C 均值聚类算法在高炉料 面分类中的应用. 北京科技大学学报,2012,34( 6) : 683) [3] Yi S,Xu Y M,Ma Z W. An expert system for abnormal status diagnosis on blast furnace. Metall Ind Autom,2002,26( 1) : 15 ( 易粟,徐用懋,马竹捂. 高炉异常炉况诊断专家系统. 冶金 自动化,2002,26( 1) : 15) ·513·
·514 北京科技大学学报 第36卷 4]Bi X G,Yang X P,Li H Y,et al.Study on the prediction expert math.Iron Steel,2001,36(3):12 system for abnormal furnace conditions.Henan Metall,2011, (吴胜利,刘茂林.基于模糊数学的高炉炉况预测模型.钢 19(4):5 铁,2001,36(3):12) (毕学功,杨绪平,李宏玉,等。高炉异常炉况预报专家系统 [14]Pan L,Yuan Q,Tang L,et al.Fault diagnosis for blast fumace 河南治金,2011,19(4):5) production process based on intelligent fuzzy system.Cent South 5]Liang D,Bai C G,Qiu G B,et al.Research on intelligence diag- Unig Sci Technol,2007,38 (Suppl 1)129 nose method for blast fumace operation.fron Steel Res,2006, (潘炼,袁琦,唐磊,等.基于智能模糊系统的高炉炉况故障 18(9):56 诊断.中南大学学报:自然科学版,2007,38(增刊1):129) (梁栋,白晨光,邱贵宝,等。高炉炉况的智能化判断方法. [15]Guo H,Wang L,Liu H P.Integrating kernel principal compo- 钢铁研究学报,2006,18(9):56) nent analysis with least squares support vector machines for time [6]Lu HS,Gao B,Zhao L G,et al.Neural network expert system of series forecasting problems.J Unir Sci Technol Beijing,2006, forecasting blast fumace operational conditions.IUni Sci Technol 28(3):303 Beng,2002,24(3):276 (郭辉,王玲,刘贺平.基于核主成分分析与最小二乘支持 (卢虎生,高斌,赵利国,等.高炉炉况判断神经网络专家系 向量机结合处理时间序列预测问题.北京科技大学学报, 统.北京科技大学学报,2002,24(3):276) 2006,28(3):303) Qu F,Wu M,Cao W H,et al.BF status diagnosis approach [16]Lu HJ,Gao Y T,Lu X N,et al.Range interval optimization of based on SVM.Iron Steel,2007,42(10)17 prestressed anchors based on the least squares support vector ma- (曲飞,吴敏,曹卫华,等.基于支持向量机的高炉炉况诊断 chine and artificial fish swarm algorithm.J Unir Sci Technol Bei- 方法.钢铁,2007,42(10):17) mg,2010,32(1):133 [8]An J Q,Wu M,He Y,et al.Blast furnace status diagnosis based (卢宏建,高永涛,卢小娜,等.基于最小二乘支持向量机和 on burden surface temperature field feature.J Zhejiang Unir Eng 人工鱼群算法的预应力锚杆布置间距优化.北京科技大学学 Sci,2010,44(7):1276 报,2010,32(1):133) (安剑奇,吴敏,何勇,等.基于料面温度场特征的高炉炉况 [17]He MG,Yin G F,Lin LJ,et al.Recognition model of ultrason- 诊断方法.浙江大学学报:工学版,2010,44(7):1276) ic flaws based on PSVM.J Sichuan Univ Eng Sci Ed,2010, 9]Gao C H,Jian L,Chen J M,et al.Data-driven modeling and pre- 42(6):232 dictive algorithm for complex blast fumace ironmaking process.Ac- (何明格,殷国官,林丽君,等.基于概率支持向量机原理的 ta Autom Sin,2009,35(6):725 超声缺陷识别模型研究.四川大学学报:工程科学版,2010, (皓传厚,渐令,陈积明,等.复杂高炉炼铁过程的数据驱动 42(6):232) 建模及预测算法.自动化学报,2009,35(6):725) [18]Lin H T,Lin C J,Weng R C.A note on Platt's probabilistic [0]Xu F F,Cheng S S,Zhao H B,et al.Model of pressure drop of outputs for support vector machines.Mach Learn,2007,68(3): lump zone in blast fumace.J /ron Steel Res,2013,25(5):8 267 (续飞飞,程树森,赵宏博,等.高炉块状带压差模型.钢铁 [19]Paksoy A,Goktiirk M.Information fusion with Dempster-Shafer 研究学报,2013,25(5):8) evidence theory for software defect prediction.Procedia Comput [11]Liang J B,Wu M.Hu B.Application of diagnosis model of fur- Sa,2011,3:600 nace working status in LianGang No.6 Blast Furnace.Equip 20]Oukhellou L,Debiolles A,Denoeux T,et al.Fault diagnosis in Manuf Technol,2009(6):167 railway track circuits using Dempster-Shafer classifier fusion.Eng (梁剑波,吴敏,胡波.炉况诊断系统在涟钢6“高炉的运用 Appl Artif Intell,2010,23 (1)117 装备制造技术,2009(6):167) 21]Xu X B,Wang Y C,Wen C L.Information-fusion method for [12]Baraldi P,Compare M,Zio E.Maintenance policy performance fault diagnosis based on reliability evaluation of evidence.Control assessment in presence of imprecision based on Dempster-Shafer Theory Appl,2011,28(4):504 theory of evidence.Inf Sci,2013,245(1)112 (徐晓滨,王玉成,文成林.评估诊断证据可靠性的信息融 [13]Wu S L,Liu M L.Forecasting model of BF status based on fuzzy 合故障诊断方法.控制理论与应用,2011,28(4):504)
北 京 科 技 大 学 学 报 第 36 卷 [4] Bi X G,Yang X P,Li H Y,et al. Study on the prediction expert system for abnormal furnace conditions. Henan Metall,2011, 19( 4) : 5 ( 毕学功,杨绪平,李宏玉,等. 高炉异常炉况预报专家系统. 河南冶金,2011,19( 4) : 5) [5] Liang D,Bai C G,Qiu G B,et al. Research on intelligence diagnose method for blast furnace operation. J Iron Steel Res,2006, 18( 9) : 56 ( 梁栋,白晨光,邱贵宝,等. 高炉炉况的智能化判断方法. 钢铁研究学报,2006,18( 9) : 56) [6] Lu H S,Gao B,Zhao L G,et al. Neural network expert system of forecasting blast furnace operational conditions. J Univ Sci Technol Beijing,2002,24( 3) : 276 ( 卢虎生,高斌,赵利国,等. 高炉炉况判断神经网络专家系 统. 北京科技大学学报,2002,24( 3) : 276) [7] Qu F,Wu M,Cao W H,et al. BF status diagnosis approach based on SVM. Iron Steel,2007,42( 10) : 17 ( 曲飞,吴敏,曹卫华,等. 基于支持向量机的高炉炉况诊断 方法. 钢铁,2007,42( 10) : 17) [8] An J Q,Wu M,He Y,et al. Blast furnace status diagnosis based on burden surface temperature field feature. J Zhejiang Univ Eng Sci,2010,44( 7) : 1276 ( 安剑奇,吴敏,何勇,等. 基于料面温度场特征的高炉炉况 诊断方法. 浙江大学学报: 工学版,2010,44( 7) : 1276) [9] Gao C H,Jian L,Chen J M,et al. Data-driven modeling and predictive algorithm for complex blast furnace ironmaking process. Acta Autom Sin,2009,35( 6) : 725 ( 皓传厚,渐令,陈积明,等. 复杂高炉炼铁过程的数据驱动 建模及预测算法. 自动化学报,2009,35( 6) : 725) [10] Xu F F,Cheng S S,Zhao H B,et al. Model of pressure drop of lump zone in blast furnace. J Iron Steel Res,2013,25( 5) : 8 ( 续飞飞,程树森,赵宏博,等. 高炉块状带压差模型. 钢铁 研究学报,2013,25( 5) : 8) [11] Liang J B,Wu M,Hu B. Application of diagnosis model of furnace working status in LianGang No. 6 Blast Furnace. Equip Manuf Technol,2009( 6) : 167 ( 梁剑波,吴敏,胡波. 炉况诊断系统在涟钢 6# 高炉的运用. 装备制造技术,2009( 6) : 167) [12] Baraldi P,Compare M,Zio E. Maintenance policy performance assessment in presence of imprecision based on Dempster-Shafer theory of evidence. Inf Sci,2013,245( 1) : 112 [13] Wu S L,Liu M L. Forecasting model of BF status based on fuzzy math. Iron Steel,2001,36( 3) : 12 ( 吴胜利,刘茂林. 基于模糊数学的高炉炉况预测模型. 钢 铁,2001,36( 3) : 12) [14] Pan L,Yuan Q,Tang L,et al. Fault diagnosis for blast furnace production process based on intelligent fuzzy system. J Cent South Univ Sci Technol,2007,38( Suppl 1) : 129 ( 潘炼,袁琦,唐磊,等. 基于智能模糊系统的高炉炉况故障 诊断. 中南大学学报: 自然科学版,2007,38( 增刊1) : 129) [15] Guo H,Wang L,Liu H P. Integrating kernel principal component analysis with least squares support vector machines for time series forecasting problems. J Univ Sci Technol Beijing,2006, 28( 3) : 303 ( 郭辉,王玲,刘贺平. 基于核主成分分析与最小二乘支持 向量机结合处理时间序列预测问题. 北京科技大学学报, 2006,28( 3) : 303) [16] Lu H J,Gao Y T,Lu X N,et al. Range interval optimization of prestressed anchors based on the least squares support vector machine and artificial fish swarm algorithm. J Univ Sci Technol Beijing,2010,32( 1) : 133 ( 卢宏建,高永涛,卢小娜,等. 基于最小二乘支持向量机和 人工鱼群算法的预应力锚杆布置间距优化. 北京科技大学学 报,2010,32( 1) : 133) [17] He M G,Yin G F,Lin L J,et al. Recognition model of ultrasonic flaws based on PSVM. J Sichuan Univ Eng Sci Ed,2010, 42( 6) : 232 ( 何明格,殷国富,林丽君,等. 基于概率支持向量机原理的 超声缺陷识别模型研究. 四川大学学报: 工程科学版,2010, 42( 6) : 232) [18] Lin H T,Lin C J,Weng R C. A note on Platt’s probabilistic outputs for support vector machines. Mach Learn,2007,68( 3) : 267 [19] Paksoy A,Gktürk M. Information fusion with Dempster-Shafer evidence theory for software defect prediction. Procedia Comput Sci,2011,3: 600 [20] Oukhellou L,Debiolles A,Denux T,et al. Fault diagnosis in railway track circuits using Dempster-Shafer classifier fusion. Eng Appl Artif Intell,2010,23( 1) : 117 [21] Xu X B,Wang Y C,Wen C L. Information-fusion method for fault diagnosis based on reliability evaluation of evidence. Control Theory Appl,2011,28( 4) : 504 ( 徐晓滨,王玉成,文成林. 评估诊断证据可靠性的信息融 合故障诊断方法. 控制理论与应用,2011,28( 4) : 504) ·514·