CN100583136C - 一种高炉炼铁专家系统决策合成方法 - Google Patents
一种高炉炼铁专家系统决策合成方法 Download PDFInfo
- Publication number
- CN100583136C CN100583136C CN200710164605A CN200710164605A CN100583136C CN 100583136 C CN100583136 C CN 100583136C CN 200710164605 A CN200710164605 A CN 200710164605A CN 200710164605 A CN200710164605 A CN 200710164605A CN 100583136 C CN100583136 C CN 100583136C
- Authority
- CN
- China
- Prior art keywords
- decision
- making
- blast furnace
- composite operator
- model
- 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.)
- Expired - Fee Related
Links
- 238000003723 Smelting Methods 0.000 title claims abstract description 12
- 238000001308 synthesis method Methods 0.000 title claims description 9
- 239000013598 vector Substances 0.000 claims abstract description 51
- 238000012549 training Methods 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 35
- 239000002131 composite material Substances 0.000 claims description 131
- 238000012417 linear regression Methods 0.000 claims description 39
- 238000012706 support-vector machine Methods 0.000 claims description 17
- 229910052717 sulfur Inorganic materials 0.000 claims description 15
- 230000002194 synthesizing effect Effects 0.000 claims description 15
- 229910052799 carbon Inorganic materials 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 6
- 238000002790 cross-validation Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 abstract description 10
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 5
- 229910052742 iron Inorganic materials 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000010276 construction Methods 0.000 abstract description 2
- 230000006399 behavior Effects 0.000 abstract 1
- XWHPIFXRKKHEKR-UHFFFAOYSA-N iron silicon Chemical compound [Si].[Fe] XWHPIFXRKKHEKR-UHFFFAOYSA-N 0.000 description 13
- 238000010189 synthetic method Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 239000003245 coal Substances 0.000 description 4
- 238000003786 synthesis reaction Methods 0.000 description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 3
- 230000004927 fusion Effects 0.000 description 3
- 229920001296 polysiloxane Polymers 0.000 description 3
- 229910052710 silicon Inorganic materials 0.000 description 3
- 239000010703 silicon Substances 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000003034 coal gas Substances 0.000 description 2
- 239000000571 coke Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000003062 neural network model Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000004886 process control Methods 0.000 description 2
- 208000011580 syndromic disease Diseases 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003137 locomotive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000006276 transfer reaction Methods 0.000 description 1
Landscapes
- Manufacture Of Iron (AREA)
Abstract
本发明公开了一种高炉炼铁专家系统决策合成方法。以多种类型高炉专家系统的单独决策为输入,构造了四种决策合成算子模型,采用线性和非线性回归技术对决策合成算子模型进行参数辨识,以高炉炉况特征向量及各类型高炉专家系统的历史表现为选择依据,把最优化决策合成算子选择问题转化为一个分类器构造问题,采用非线性多分类支持向量机模型进行分类器建模,通过快速分类器训练算法得到最优分类器,进而利用该分类器对决策合成算子模型参数进行递推式的二次校正,由最优分类器动态选择出适合于当前炉况的最优决策合成算子形成高炉操作优化决策。对高炉冶炼过程的智能控制与决策的合成具有普遍的通用性,能提高高炉专家系统决策准确度与系统鲁棒性。
Description
技术领域
本发明涉及高炉炼铁专家系统领域的知识融合及决策合成方法,尤其涉及高炉炼铁专家系统连续型决策的合成方法。
背景技术
高炉炼铁过程是一个高度复杂的过程,其运动机制往往具有非线性、大时滞、大噪声、分布参数等特征。同时,高炉本身是一个集传热与化学反应耦合的开放体。以高炉专家系统的方式实现高炉冶炼过程的智能控制与决策,是冶金科技发展的前沿课题[1]。按专家系统以何种机制形成决策的不同,传统高炉专家系统可以划分为以两种:1)基于知识推理的专家系统(例如文献[2]):以高炉炼铁过程的实际生产经验为主要知识来源,形成知识库进而通过(模糊)逻辑推理等方法来产生目标决策;2)基于理论模型的高炉专家系统(例如文献[3]):从理论角度分析高炉炼铁过程中的物理化学过程,建立相应的数学模型并以理论优化计算方式产生目标决策。大量实践证明,上述两种专家系统各具优势,例如,基于知识推理的专家系统在炉况异常的情况下,其决策的准确度要优于基于理论模型的专家系统,而在炉况正常的时候,基于理论模型的专家系统往往能给出更准确的结果。因此,利用知识融合及决策合成技术巧妙地把各种类型高炉专家系统的知识决策融合在一起,能在系统鲁棒性和决策精度等方面大大提高传统高炉专家系统的质量。
参考文献:
[1].刘祥官,刘芳.高炉炼铁过程优化与智能控制系统.北京:冶金工业出版社,2003.
[2].易粟,马竹梧.高炉异常炉况诊断专家系统.冶金自动化,2003,26(1).
[3].郜传厚,周志敏,邵之江.高炉冶炼过程的混沌性解析.物理学报,2005,54(4).
发明内容
本发明的目的是提供一种高炉炼铁专家系统决策合成方法。
包括如下步骤:
1)以多种类型高炉专家系统的单独决策为输入,构造四种应用于决策值为连续实数的决策合成算子模型,即朴素选择决策合成算子、简单平均决策合成算子、加权平均决策合成算子和非线性回归决策合成算子;
2)采用线性回归和非线性回归技术分别辨识出上述四种决策合成算子的模型参数;
3)将最优决策合成算子选择问题转化为多分类器问题;
4)以多分类支持向量机为分类器模型,采用快速SMO算法为训练方法求解出最优分类器;
5)利用最优分类器对决策合成算子的模型参数进行递推式的二次校正,得到最终的决策合成算子模型;
6)依次利用最优分类器和最终决策合成算子进行决策合成,获得高炉操作优化决策。
所述的朴素选择决策合成算子:选择其中一个专家系统决策为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,那么共有M种具体的决策合成算子,记为S1{h=d1},S2{h=d2},K,SM{h=dM},其中Si{h=di},i=1,2,K M,表示决策合成
算子Si按h=di的方式产生合成后的决策。
所述的简单平均决策合成算子:以多种类型高炉专家系统决策的平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合算子为
表示决策合成算子SM+1按 的方式产生合成后的决策。
所述的加权平均决策合成算子:以多种类型高炉专家系统决策的加权平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,M个高炉专家系统单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合成算子为
表示决策合成算子SM+2按 的方式产生合成后的决策,其中λi>0,i=1,2,K M是各专家系统决策的权重分配,且满足
所述的非线性回归决策合成算子:以多种类型高炉专家系统决策为回归项构造非线性回归模型,并以回归模型输出值作为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合成算子为SM+3{h=f(d1,d2,K,dM)},其中f是非线性函数,采用非线性回归支持向量机模型来表示f,记该决策合成算子为
表示决策合成算子SM+3按 的方式产生合成后的决策,其中V是支持向量的个数,βi是核函数系数,Ki是与第i个支持向量对应的非线性核函数。
所述的采用线性回归和非线性回归技术辨识决策合成算子模型参数:以多种类型高炉专家系统单独决策值的历史数据和真实决策值构造训练样本集Z={Z1,Z2,K ZN},其中N>0是历史数据的总数,即样本总数,Zi={di1,di2,K,diM,yi}是第i组样本,dik,k=1,2,K M,表示专家k在i时的单独决策,yi是真实的决策。以Z为训练样本集,采用最小二乘算法辨识加权平均决策合成算子的模型参数,即式 中的λi,i=1,2,K M,采用快速SMO支持向量机训练算法辨识非线性回归决策合成算子模型的参数,即式 中的βi和支持向量。
所述的将最优决策合成算子选择问题转化为多分类器问题:将步骤2)中辨识后的四种决策合成算子模型分别应用于训练样本Z={Z1,Z2,K ZN},对于每组样本Zi={di1,di2,K,diM,yi},得出各种决策合成算子产生的决策合成值{hi1,hi2,K,hiM,hi,M+1,hi,M+2,hi,M+3},以真实决策值yi为参考标准,可确认出针对样本Zi的最优决策合成算子其数学表达式为
采用交叉验证方法选择出高炉炉况特征向量为输入向量,记训练样本Z={Z1,Z2,K ZN}对应的炉况特征向量数据为X={X1,X2,K XN},Xi=[xi1,xi2,K,xiQ],其中Q是高炉特征向量的维数。根据式 以[xi1,xi2,K,xiQ]为输入向量,以oi为正确的类别标签,形成用于分类问题的训练样本C={C1,C2,K CN},Ci={xi1,xi2,K xiQ,oi},至此,将最优决策合成算子选择问题转化为多分类器问题。
所述的采用多分类支持向量机为分类器模型,以快速SMO算法为训练方法求解出最优分类器:以非线性多分类支持向量机模型来对多分类器问题建模,以C={C1,C2,K CN},Ci={xi1,xi2,K xiQ,oi}为训练样本数据,调用快速SMO支持向量机训练算法,记训练好后的分类器为F,其数学表达式为,
F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3},
其中X是炉况特征向量空间,即分类器F是由炉况特征向量空间到决策合成算子集合空间的一个映射。
所述的利用最优分类器对决策合成算子的模型参数进行递推式的二次校正:将步骤4)的最优分类器,F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3},应用于炉况数据X={X1,X2,K XN},根据映射F将X划分了M+3个子集,记为X1,X2,K,XM,XM+1,XM+2,XM+3,由X={X1,X2,K XN}与决策合成算子模型原始训练样本Z={Z1,Z2,K ZN}的一一对应关系,F亦将Z划分成M+3个子集,记为Z1,Z2,K,ZM,ZM+1,ZM+2,ZM+3,与加权平均决策合成算子SM+2对应的样本子集是ZM+2,由于ZM+2对SM+2模型参数辨识更具指导意义,因此增强样本子集ZM+2的信度,对SM+2模型进行二次校正,以SM+2为初始模型,以ZM+2为需要强调的更新样本,采用带遗忘因子的递推最小二乘辨识算法完成对SM+2的二次校正,修正后的模型记为S′M+2,同理,与非线性回归决策合成算子SM+3对应的样本子集是ZM+3,对SM+3模型进行二次校正,以SM+3为初始模型,以ZM+3为需要更新样本,采用带遗忘因子的递推支非线性回归持向量机训练算法完成对SM+3的二次校正,修正后的模型记为S′M+3。
所述的依次利用最优分类器和最终决策合成算子形成高炉操作优化决策:以当前炉况数据特征向量X′=[x′1,x′2,K,x′Q]为最优分类器的输入向量,其中Q是炉况数据特征向量的维数,利用式F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3}最优分类选出最适合于当前炉况的决策合成算子So∈{S1,S2,K,SM,SM+1,S′M+2,S′M+3},被选出的最终决策合成算子So以各类型高炉专家系统单独决策值[d′1,d′2,K,d′M]为模型的输入,产生合成后的高炉操作优化决策h′=So(d′1,d′2,K,d′M)。
本发明对高炉冶炼过程的智能控制与决策的合成具有普遍的通用性,能提高高炉专家系统决策准确度与系统鲁棒性,提供的决策合成算子模型即考虑了线性模型又考虑了非线性模型,对于高炉炼铁过程中各种连续型决策具有通用性;采用的线性回归和非线性回归算法求解时间对样本数目大小不敏感,且有很好的收敛性和鲁棒性,可靠性高,保证了决策合成问题求解的实时性;将最优决策合成算子选择问题转化为多分类器问题有巧妙性,使得决策合成策略可以根据当前炉况动态选择最佳的合成方案;利用最优分类器对决策合成算子的模型参数进行递推式的二次校正能够很好地提高决策合成算子的模型准确度。将本发明提出的高炉炼铁专家决策合成方法,应用于某大型钢铁企业高炉铁水硅含量预报的多高炉专家系统决策合成问题中,取得了非常好的效果,合成后的高炉铁水硅含量预报值的准确度显著优于其中任一高炉专家系统单独给出的预报值,给高炉平稳可靠运行提供技术上的保障。
具体实施方式
本发明旨在建立一种相对简单但行之有效的高炉炼铁专家系统决策合成方法,分析了高炉炼铁工艺过程的特点,考虑了连续型决策的合成方法。连续型决策(决策值为连续实数)在高炉炼铁过程控制中扮演重要角色,例如冶炼过程控制中核心和基础的炉温控制问题、硅含量预报问题等都属于连续型决策问题。本发明针对连续型决策的特点,以各类型高炉专家系统的单独决策为输入,构造了几种行之有效的决策合成算子模型,采用线性回归技术和非线性回归技术对决策合成算子模型进行参数辨识,并以高炉炉况特征向量和各类型高炉专家系统的历史表现为选择最优决策合成算子的依据,并巧妙地把该最优化决策合成算子选择问题转化为一个分类器构造问题,采用专门为小样本问题设计的非线性多分类支持向量机模型进行分类器建模,通过快速分类器训练算法得到最优分类器,再利用最优分类器对决策合成算子模型参数进行递推式的二次校正,应用时,由该最优分类器动态选择出适合于当前炉况的最优决策合成算子,最后由选出的决策合成算子形成高炉操作优化决策。本发明对高炉冶炼过程的智能控制与决策具有普遍的通用性,能大大提高高炉专家系统决策准确度与系统鲁棒性。
下面以高炉铁水硅含量预报专家系统的决策合成为应用实例,阐述本发明提出的高炉炼铁专家决策合成方法的具体实施步骤。对高炉复杂系统的建模和控制是当今冶金科技发展的前沿课题,其中对高炉炉温的预测和控制是难点所在。由于高炉过程的高度复杂性和封闭性,导致了直接测量的困难。鉴于铁水硅含量间接反应炉内温度的变化,长期以来,冶金专家们建立了一系列高炉铁水含硅量预报的专家系统:基于多元回归模型的专家系统,基于时间序列模型的专家系统,基于神经网络模型的专家系统,基于知识的推理的专家系统等。上述专家系统各具优势,利用决策合成方法取长补短能大大提高高炉铁水含硅量预报专家系统的决策准确度与系统鲁棒性。
高炉炼铁专家决策合成方法包括如下步骤:
1)针对高炉炼铁过程中连续型决策的特性,以各类型高炉专家系统的单独决策为输入,构造四种应用于决策值为连续实数的决策合成算子模型,即朴素选择决策合成算子、简单平均决策合成算子、加权平均决策合成算子和非线性回归决策合成算子。上述决策合成算子模型即考虑了线性模型又考虑了非线性模型,对于高炉炼铁过程中各种连续型决策(如高炉铁水硅含量预报)具有通用性,在高炉铁水硅含量预报的应用实例中,采用四种类型专家系统,它们分别是:基于多元回归模型的专家系统,基于时间序列模型的专家系统,基于神经网络模型的专家系统和基于知识的推理的专家系统。
2)采用线性回归和非线性回归技术分别辨识出上述四种决策合成算子的模型参数。在高炉铁水硅含量预报的应用实例中,线性回归采用了最小二乘方法,非线性回归采用了非线性回归支持向量机。
3)将最优决策合成算子选择问题转化为多分类器问题;
4)以多分类支持向量机为分类器模型,采用快速SMO算法为训练方法求解出最优分类器;
5)利用最优分类器对决策合成算子的模型参数进行递推式的二次校正,得到最终的决策合成算子模型;
6)依次利用最优分类器和最终决策合成算子进行决策合成,获得高炉操作优化决策,在高炉铁水硅含量预报的应用实例中,决策合成结果是硅含量预报值。
所述的朴素选择决策合成算子:其合成策略是选择其中一个专家决策为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,那么共有M种具体的决策合成算子,记为S1{h=d1},S2{h=d2},K,SM{h=dM},其中Si{h=di},i=1,2,K M,表示决策合成算子Si按h=di的方式产生合成后的决策。在高炉铁水硅含量预报的应用实例中,M=4。
所述的简单平均决策合成算子:其合成策略是以各专家决策的平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合算子为
表示决策合成算子SM+1按 的方式产生合成后的决策,简单平均决策合成算子对噪声不敏感具有好的鲁棒性。
所述的加权平均决策合成算子:其合成策略是以各专家决策的加权平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,M个高炉专家系统单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合成算子为
表示决策合成算子SM+2按 的方式产生合成后的决策,其中λi>0,i=1,2,K M是各专家决策的权重分配,且满足 加权平均决策合成算子通过权重分配能较准确得产生合成后的决策。
所述的非线性回归决策合成算子:以各专家决策为回归项构造非线性回归模型,并以回归模型输出值作为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合成算子为SM+3{h=f(d1,d2,K,dM)},其中f是非线性函数,采用非线性回归支持向量机模型来表示f,记该决策合成算子为
表示决策合成算子SM+3按 的方式产生合成后的决策,其中V是支持向量的个数,βi是核函数系数,Ki是与第i个支持向量对应的非线性核函数,在高炉铁水硅含量预报的应用实例中,Ki采用的是径向基核函数[1]。非线性回归决策合成算子能够对有效地描述问题的非线性特征。
所述的采用线性回归和非线性回归技术辨识决策合成算子模型参数:以各类型高炉专家系统单独决策值的历史数据和真实决策值构造训练样本集Z={Z1,Z2,K ZN},其中N>0是历史数据的总数,即样本总数,Zi={di1,di2,K,diM,yi}是第i组样本,dik,k=1,2,K M,表示专家k在i时的单独决策,yi是真实的决策。以在高炉铁水硅含量预报为例,dik表示的是表示专家k在i时的单给出的硅含量百分比,yi是测量出的实际硅含量百分比。以Z为训练样本集,采用最小二乘算法辨识加权平均决策合成算子的模型参数,即式 中的λi,i=1,2,K M,采用快速SMO支持向量机训练算法辨识非线性回归决策合成算子模型的参数,即式 中的βi和支持向量。
所述的将最优决策合成算子选择问题转化为多分类器问题:将步骤2)中辨识后的四种决策合成算子模型分别应用于训练样本Z={Z1,Z2,K ZN},对于每组样本Zi={di1,di2,K,diM,yi},得出各种决策合成算子产生的决策合成值{hi1,hi2,K,hiM,hi,M+1,hi,M+2,hi,M+3},以真实决策值yi为参考标准,可确认出针对样本Zi的最优决策合成算子其数学表达式为
采用交叉验证方法选择出高炉炉况特征向量为输入向量,记训练样本Z={Z1,Z2,K ZN}对应的炉况特征向量数据为X={X1,X2,K XN},Xi=[xi1,xi2,K,xiQ],其中Q是高炉特征向量的维数。以在高炉铁水硅含量预报为例,高炉炉况特征向量包含了如下物理量:铁量差、透气性、喷煤量、风温、料批、风量、富氧量、热风压力、炉顶压力、喷煤量、热风温度、炉顶温度、煤气CO、CO2的含量、矿焦比。根据式 以[xi1,xi2,K,xiQ]为输入向量,以oi为正确的类别标签,形成用于分类问题的训练样本C={C1,C2,K CN},Ci={xi1,xi2,K xiQ,oi},至此,将最优决策合成算子选择问题转化为多分类器问题。
所述的采用多分类支持向量机为分类器模型,以快速SMO算法为训练方法求解出最优分类器:以非线性多分类支持向量机模型来对多分类器问题建模,采用sigmoid核函数,以C={C1,C2,K CN},C1={xi1,xi2,K xiQ,oi}为训练样本数据,调用快速SMO支持向量机训练算法,记训练好后的分类器为F,其数学表达式为,
F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3},
其中X是炉况特征向量空间,即分类器F是由炉况特征向量空间到决策合成算子集合空间的一个映射。
所述的利用最优分类器对决策合成算子的模型参数进行递推式的二次校正:将步骤4)的最优分类器,F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3},应用于炉况数据X={X1,X2,K XN},根据映射F将X划分了M+3个子集,记为X1,X2,K,XM,XM+1,XM+2,XM+3,由X={X1,X2,K XN}与决策合成算子模型原始训练样本Z={Z1,Z2,K ZN}的一一对应关系,F亦将Z划分成M+3个子集,记为Z1,Z2,K,ZM,ZM+1,ZM+2,ZM+3,与加权平均决策合成算子SM+2对应的样本子集是ZM+2,由于ZM+2对SM+2模型参数辨识更具指导意义,因此增强样本子集ZM+2的信度,对SM+2模型进行二次校正,以SM+2为初始模型,以ZM+2为需要强调的更新样本,采用带遗忘因子的递推最小二乘辨识算法完成对SM+2的二次校正,修正后的模型记为S′M+2,其数学表达式为,
S′M+2=El(SM+2,ZM+2)
其中Rl表示带遗忘因子的递推最小二乘辨识算法。同理,与非线性回归决策合成算子SM+3对应的样本子集是ZM+3,由于ZM+3对SM+3模型参数辨识更具指导意义,因此增强样本子集ZM+3的信度,对SM+3模型进行二次校正,以SM+3为初始模型,以ZM+3为需要强调的更新样本,采用带遗忘因子的递推支非线性回归持向量机训练算法完成对SM+3的二次校正,修正后的模型记为S′M+3,其数学表达式为,
S′M+3=Rs(SM+3,ZM+3)
其中Rs表示带遗忘因子的递推支非线性回归持向量机训练算法。
所述的依次利用最优分类器和最终决策合成算子形成高炉操作优化决策:以当前炉况数据特征向量X′=[x′1,x′2,K,x′Q]为最优分类器的输入向量,其中Q是炉况数据特征向量的维数,在高炉铁水硅含量预报的应用实例中,高炉炉况特征向量包含了如下物理量:铁量差、透气性、喷煤量、风温、料批、风量、富氧量、热风压力、炉顶压力、喷煤量、热风温度、炉顶温度、煤气CO、CO2的含量、矿焦比。利用式F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3}最优分类选出最适合于当前炉况的决策合成算子So∈{S1,S2,K,SM,SM+1,S′M+2,S′M+3},被选出的最终决策合成算子So以各类型高炉专家系统单独决策值[d′1,d′2,K,d′M]为模型的输入,产生合成后的高炉操作优化决策h′=So(d′1,d′2,K,d′M)。
Claims (3)
1.一种高炉炼铁专家系统决策合成方法,其特征在于包括如下步骤:
1)以多种类型高炉专家系统的单独决策为输入,构造四种应用于决策值为连续实数的决策合成算子模型,即朴素选择决策合成算子、简单平均决策合成算子、加权平均决策合成算子和非线性回归决策合成算子;
2)采用线性回归和非线性回归技术分别辨识出上述四种决策合成算子的模型参数;
3)将最优决策合成算子选择问题转化为多分类器问题;
4)以多分类支持向量机为分类器模型,采用快速SMO算法为训练方法求解出最优分类器;
5)利用最优分类器对决策合成算子的模型参数进行递推式的二次校正,得到最终的决策合成算子模型;
6)依次利用最优分类器和最终决策合成算子进行决策合成,获得高炉操作优化决策;
所述的朴素选择决策合成算子:选择其中一个专家系统决策为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,那么共有M种具体的决策合成算子,记为S1{h=d1},S2{h=d2},K,SM{h=dM},其中Si{h=di},i=1,2,K M,表示决策合成算子Si按h=di的方式产生合成后的决策
所述的简单平均决策合成算子:以多种类型高炉专家系统决策的平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合算子为
表示决策合成算子SM+1按 的方式产生合成后的决策;
所述的加权平均决策合成算子:以多种类型高炉专家系统决策的加权平均值为合成后的决策,假设有M个高炉专家系统参与决策合成,M个高炉专家系统单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合成算子为:
表示决策合成算子SM+2按 的方式产生合成后的决策,其中λi>0,i=1,2,K M是各专家系统决策的权重分配,且满足
所述的非线性回归决策合成算子:以多种类型高炉专家系统决策为回归项构造非线性回归模型,并以回归模型输出值作为合成后的决策,假设有M个高炉专家系统参与决策合成,它们单独给出的决策值分别为d1,d2,K,dM,合成后的决策记为h,记该决策合成算子为SM+3{h=f(d1,d2,K,dM)},其中f是非线性函数,采用非线性回归支持向量机模型来表示f,记该决策合成算子为
表示决策合成算子SM+3按 的方式产生合成后的决策,其中V是支持向量的个数,βi是核函数系数,Ki是与第i个支持向量对应的非线性核函数;
所述的将最优决策合成算子选择问题转化为多分类器问题:将步骤2)中辨识后的四种决策合成算子模型分别应用于训练样本Z={Z1,Z2,K ZN},对于每组样本Zi={di1,di2,K,diM,yi},得出各种决策合成算子产生的决策合成值{hi1,hi2,K,hiM,hi,M+1,hi,M+2,hi,M+3},以真实决策值yi为参考标准,可确认出针对样本Zi的最优决策合成算子其数学表达式为
采用交叉验证方法选择出高炉炉况特征向量为输入向量,记训练样本Z={Z1,Z2,K ZN}对应的炉况特征向量数据为X={X1,X2,K XN},Xi=[xi1,xi2,K,xiQ],其中Q是高炉特征向量的维数据式
所述的利用最优分类器对决策合成算子的模型参数进行递推式的二次校正:将步骤4)的最优分类器,F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3},应用于炉况数据X={X1,X2,K XN},根据映射F将X划分了M+3个子集,记为X1,X2,K,XM,XM+1,XM+2,XM+3,由X={X1,X2,K XN}与决策合成算子模型原始训练样本Z={Z1,Z2,K ZN}的一一对应关系,F亦将Z划分成M+3个子集,记为Z1,Z2,K,ZM,ZM+1,ZM+2,ZM+3,与加权平均决策合成算子SM+2对应的样本子集是ZM+2,由于ZM+2对SM+2模型参数辨识更具指导意义,因此增强样本子集ZM+2的信度,对SM+2模型进行二次校正,以SM+2为初始模型,以ZM+2为需要强调的更新样本,采用带遗忘因子的递推最小二乘辨识算法完成对SM+2的二次校正,修正后的模型记为S′M+2,同理,与非线性回归决策合成算子SM+3对应的样本子集是ZM+3,对SM+3模型进行二次校正,以SM+3为初始模型,以ZM+3为需要更新样本,采用带遗忘因子的递推支非线性回归持向量机训练算法完成对SM+3的二次校正,修正后的模型记为S′M+3;
所述的依次利用最优分类器和最终决策合成算子形成高炉操作优化决策:以当前炉况数据特征向量X′=[x′1,x′2,K,x′Q]为最优分类器的输入向量,其中Q是炉况数据特征向量的维数,利用式F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3}最优分类选出最适合于当前炉况的决策合成算子So∈{S1,S2,K,SM,SM+1,S′M+2,S′M+3},被选出的最终决策合成算子So以各类型高炉专家系统单独决策值[d′1,d′2,K,d′M]为模型的输入,产生合成后的高炉操作优化决策h′=So(d′1,d′2,K,d′M)。
2.根据权利要求1所述的一种高炉炼铁专家系统决策合成方法,其特征在于所述的采用线性回归和非线性回归技术辨识决策合成算子模型参数:以多种类型高炉专家系统单独决策值的历史数据和真实决策值构造训练样本集Z={Z1,Z2,K ZN},其中N>0是历史数据的总数,即样本总数,Zi={di1,di2,K,diM,yi}是第i组样本,dik,k=1,2,K M,表示专家k在i时的单独决策,yi是真实的决策;以Z为训练样本集,采用最小二乘算法辨识加权平均决策合成算子的模型参数,即式 中的λi,i=1,2,K M,采用快速SMO支持向量机训练算法辨识非线性回归决策合成算子模型的参数,即式 中的βi和支持向量。
3.根据权利要求1所述的一种高炉炼铁专家系统决策合成方法,其特征在于所述的采用多分类支持向量机为分类器模型,以快速SMO算法为训练方法求解出最优分类器:以非线性多分类支持向量机模型来对多分类器问题建模,以C={C1,C2,K CN},Ci={xi1,xi2,KiQ,oi}为训练样本数据,调用快速SMO支持向量机训练算法,记训练好后的分类器为F,其数学表达式为,
F:X→{S1,S2,K,SM,SM+1,SM+2,SM+3},
其中X是炉况特征向量空间,即分类器F是由炉况特征向量空间到决策合成算子集合空间的一个映射。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710164605A CN100583136C (zh) | 2007-12-21 | 2007-12-21 | 一种高炉炼铁专家系统决策合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710164605A CN100583136C (zh) | 2007-12-21 | 2007-12-21 | 一种高炉炼铁专家系统决策合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101211425A CN101211425A (zh) | 2008-07-02 |
CN100583136C true CN100583136C (zh) | 2010-01-20 |
Family
ID=39611436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710164605A Expired - Fee Related CN100583136C (zh) | 2007-12-21 | 2007-12-21 | 一种高炉炼铁专家系统决策合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100583136C (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102031319B (zh) * | 2009-09-30 | 2013-05-29 | 鞍钢股份有限公司 | 一种高炉铁水含硅量的预报方法 |
CN102925602B (zh) * | 2012-11-22 | 2014-11-26 | 武汉钢铁(集团)公司 | 高炉操作炉型维护方法 |
CN103729571B (zh) * | 2014-01-23 | 2017-02-15 | 中南大学 | 一种高炉炼铁过程一氧化碳利用率的建模方法 |
CN106909972B (zh) * | 2017-01-09 | 2020-04-10 | 浙江大学 | 一种传感器数据校准模型的学习方法 |
CN111383723B (zh) * | 2020-04-16 | 2023-06-27 | 山东钢铁股份有限公司 | 一种高炉炼铁成分预控方法 |
-
2007
- 2007-12-21 CN CN200710164605A patent/CN100583136C/zh not_active Expired - Fee Related
Non-Patent Citations (6)
Title |
---|
基于支持向量机的高炉铁水硅含量多类别分类. 渐令,龚淑华,王义康.浙江大学学报(理学版),第34卷第3期. 2007 |
基于支持向量机的高炉铁水硅含量多类别分类. 渐令,龚淑华,王义康.浙江大学学报(理学版),第34卷第3期. 2007 * |
高炉冶炼过程的混沌性解析. 郜传厚,周志敏,邵之江.物理学报,第54卷第4期. 2005 |
高炉冶炼过程的混沌性解析. 郜传厚,周志敏,邵之江.物理学报,第54卷第4期. 2005 * |
高炉异常炉况诊断专家系统. 易粟,徐用懋,马竹梧.冶金自动化. 2002 |
高炉异常炉况诊断专家系统. 易粟,徐用懋,马竹梧.冶金自动化. 2002 * |
Also Published As
Publication number | Publication date |
---|---|
CN101211425A (zh) | 2008-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109447346B (zh) | 基于灰色预测与神经网络组合模型转炉氧耗量预测方法 | |
CN102031319B (zh) | 一种高炉铁水含硅量的预报方法 | |
WO2021189739A1 (zh) | 一种基于生产计划的钢铁企业氧气负荷预测方法 | |
CN100583136C (zh) | 一种高炉炼铁专家系统决策合成方法 | |
CN107918368B (zh) | 钢铁企业煤气产生量与消耗量的动态预测方法及设备 | |
Xie et al. | Robust stochastic configuration network multi-output modeling of molten iron quality in blast furnace ironmaking | |
CN102163261A (zh) | 一种基于案例推理的钢水温度预测方法 | |
JPWO2020152993A1 (ja) | 金属材料の設計支援方法、予測モデルの生成方法、金属材料の製造方法、及び設計支援装置 | |
CN111178621A (zh) | 一种电采暖负荷预测支持向量回归机模型的参数优化方法 | |
Jiang et al. | Real-time moisture control in sintering process using offline–online NARX neural networks | |
CN105574297B (zh) | 自适应高炉铁水硅含量趋势预报方法 | |
CN1753010A (zh) | 用于钢铁企业能源优化调度的分类建模和滚动求解方法 | |
CN108845501A (zh) | 一种基于懒惰学习的高炉铁水质量自适应优化控制方法 | |
CN114066069A (zh) | 一种组合权重的副产煤气发生量预测方法 | |
CN114036827B (zh) | 基于分解的高炉炼铁多目标碳减排方法 | |
Tian et al. | A new AdaBoost. IR soft sensor method for robust operation optimization of ladle furnace refining | |
Ji et al. | Application of the improved the ELM algorithm for prediction of blast furnace gas utilization rate | |
CN116722605B (zh) | 一种基于物联网的配电网调度优化方法 | |
Ding et al. | Forecast of pv power generation based on residual correction of markov chain | |
CN111798023A (zh) | 一种炼钢烧结生产中综合焦比预测的方法 | |
CN116504329A (zh) | 基于煤数据图像化的卷积神经网络焦炭热态质量预测方法 | |
Zhang et al. | An ensemble model based on weighted support vector regression and its application in annealing heating process | |
CN103135444A (zh) | 一种钢铁生产能耗免疫预测控制模型 | |
CN116127345B (zh) | 基于深度聚类生成对抗网络的转炉炼钢工艺模式设计方法 | |
Liu et al. | Using multivariate grey model and principal component analysis to modeling the blast furnace |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100120 Termination date: 20111221 |