基于ukf神经网络的转炉炼钢工艺成本控制方法及系统的制作方法
【专利摘要】本发明提供一种基于UKF神经网络的转炉炼钢工艺成本控制方法及系统，其中的方法包括选择影响成本的控制参数；构建建模样本集；获得归一化样本集；构建前馈神经网络；采用UKF算法对前馈神经网络进行训练建模，获取神经网络参数；利用PSO粒子群算法对运用UKF算法所构建的模型进行优化，获取所构建模型的最值，并根据所构建模型的最值确定最优控制参数；根据最优控制参数成本值与所述建模样本集中的最小成本值的对比结果，确定转炉炼钢工艺的最小成本值。利用本发明，能够解决转炉炼钢成本高的问题。
【专利说明】
基于UKF神经网络的转妒炼钢工艺成本控制方法及系统
技术领域
[0001] 本发明设及炼钢技术领域，更为具体地，设及一种基于UK巧巾经网络的转炉炼钢工 艺成本控制方法及系统。
【背景技术】
[0002] 目前钢铁行业进入低谷、行业利润被无限压缩，只有降低自身成本才能寻求发展。 所W钢铁行业的降本增效是所有钢厂不懈的追求。而钢铁生产过程高溫、高危、高成本，无 法进行大规模现场。
[0003] 其中，实验碱性氧气转炉炼钢法是一种将铁水炼成钢水的炼钢过程。通过向烙池 供氧，发生氧化反应降低烙池中钢液含碳量，此炼钢法又称为转炉炼钢。通过虚拟炼钢模拟 实际冶炼过程，可为现场生产提供降本增效的可行性方案和指导性意见，具有重大意义和 经济效益。
[0004] 炉子的分类较多，较为普遍分类是顶吹转炉、底吹转炉和顶底复合吹转炉。在转炉 炼钢过程中，系统配料、操作过程等均会对炼钢的成本有着重要的作用，为进一步改进加入 原料配方、优化生产过程等生产参数，得到一个最为经济理想的冶炼过程，为企业提供优化 思路，节省成本。
[0005] 综上所述，为解决上述问题，基于虚拟炼钢模拟实际冶炼的思想，本发明提出了一 种基于UK巧巾经网络的转炉炼钢工艺成本控制方法。
【发明内容】
[0006] 鉴于上述问题，本发明的目的是提供一种基于UK巧申经网络的转炉炼钢工艺成本 控制方法及系统，能够解决转炉炼钢成本高的问题。
[0007] 本发明提供一种基于UK巧巾经网络的转炉炼钢工艺成本控制方法，包括：
[0008] 根据转炉炼钢的工艺选择影响成本的控制参数；
[0009] 利用模拟转炉炼钢平台采集不同控制参数的成本，构建建模样本集；
[0010] 将构建的建模样本集进行归一化处理，获得归一化样本集；
[0011] 根据所述归一化样本集构建前馈神经网络；
[0012] 采用UKF算法对所述前馈神经网络进行训练建模，获取神经网络参数；
[0013] 利用PSO粒子群算法对运用UKF算法所构建的模型进行优化，获取所构建模型的最 值，并根据所构建模型的最值确定最优控制参数；
[0014] 根据所述最优控制参数获取最优控制参数成本值；
[0015] 根据所述最优控制参数成本值与所述建模样本集中的最小成本值的对比结果，确 定转炉炼钢工艺的最小成本值。
[0016] 本发明还提供一种基于UK巧巾经网络的转炉炼钢工艺成本控制系统，包括
[0017] 控制参数选择单元，用于根据转炉炼钢的工艺选择影响成本的控制参数；
[0018] 建模样本集构建单元，用于利用模拟转炉炼钢平台采集不同控制参数的成本，构 建建模样本集；
[0019] 归一化样本集获取单元，用于将构建的建模样本集进行归一化处理，获得归一化 样本集；
[0020] 前馈神经网络构建单元，用于根据所述归一化样本集构建前馈神经网络；
[0021] 神经网络参数获取单元，用于采用UKF算法对所述前馈神经网络进行训练建模，获 取神经网络参数.
[0022] 最优控制参数获取单元，用于利用PSO粒子群算法对运用UKF算法所构建的模型进 行优化，获取所构建模型的最值，并根据所构建模型的最值确定最优控制参数；
[0023] 最优控制参数成本值获取单元，用于根据所述最优控制参数获取最优控制参数成 本值；
[0024] 最小成本值获取单元，用于根据所述最优控制参数成本值与所述建模样本集中的 最小成本值的对比结果，确定转炉炼钢工艺的最小成本值。
[0025] 从上面的技术方案可知，本发明提供的基于UK巧巾经网络的转炉炼钢工艺成本控 制方法及系统，在冶炼过程中的生产操作参数为信息载体，利用UK巧巾经网络挖掘原料配 方、操作参数与炼钢成本之间的关系;并通过智能优化算法利获取最低成本下的操作参数， 为实际生产最优生产提供指导，解决转炉炼钢成本较高的问题。
[0026] 为了实现上述W及相关目的，本发明的一个或多个方面包括后面将详细说明并在 权利要求中特别指出的特征。下面的说明W及附图详细说明了本发明的某些示例性方面。 然而，运些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外，本发明 旨在包括所有运些方面W及它们的等同物。
【附图说明】
[0027] 通过参考W下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面 理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：
[0028] 图1为根据本发明实施例的基于UK巧巾经网络的转炉炼钢工艺成本控制方法流程 示意图；
[0029] 图2为根据本发明实施例的基于UK巧巾经网络的转炉炼钢工艺成本控制系统逻辑 结构示意图；
[0030] 图3为根据本发明实施例的前馈神经网络结构示意图；
[0031] 图4为根据本发明实施例的所构建模型的训练样本效果图；
[0032] 图5为根据本发明实施例的所构建模型的测试样本预测精度效果图；
[0033] 图6为根据本发明实施例的PSO粒子群算法的函数获取极值流程示意图。
[0034] 在所有附图中相同的标号指示相似或相应的特征或功能。
【具体实施方式】
[0035] 在下面的描述中，出于说明的目的，为了提供对一个或多个实施例的全面理解，阐 述了许多具体细节。然而，很明显，也可W在没有运些具体细节的情况下实现运些实施例。
[0036] 针对前述提出的目前钢铁行业成本过高的问题，本发明提出了基于UK巧申经网络 的转炉炼钢工艺成本控制方法及系统，其中，本发明提出W冶炼过程中的生产操作参数为 信息载体，利用UK巧巾经网络方法挖掘原料配方、操作参数与炼钢成本之间的潜在规律;并 通过智能优化算法利用该规律获取最低成本下的操作参数，为企业的实际生产最优生产提 供指导。
[0037] 其中，需要说明的是，UKF化nscented Kalman Filter)，中文释义是无损卡尔曼滤 波、无迹卡尔曼滤波或者去芳香卡尔曼滤波，是无迹变换化T)和标准Kalman滤波体系的结 合，通过无损变换使非线性系统方程适用于线性假设下的标准Kalman滤波体系。
[0038] W下将结合附图对本发明的具体实施例进行详细描述。
[0039] 为了说明本发明提供的基于UK巧申经网络的转炉炼钢工艺成本控制方法，图1示出 了根据本发明实施例的基于UK巧巾经网络的转炉炼钢工艺成本控制方法流程。
[0040] 如图1所示，本发明提供的基于UK巧巾经网络的转炉炼钢工艺成本控制方法包括： Sl 10:根据转炉炼钢的工艺选择影响成本的控制参数；
[0041] S120:利用模拟转炉炼钢平台采集不同控制参数的成本，构建建模样本集；
[0042] S130:将构建的建模样本集进行归一化处理，获得归一化样本集；
[0043] S140:根据所述归一化样本集构建前馈神经网络；
[0044] S150:采UKF算法对前馈神经网络进行训练建模，获取神经网络参数；
[0045] S160:利用PSO粒子群算法对运用UKF算法所构建的模型进行优化，获取所构建模 型的最值，并根据所构建模型的最值确定最优控制参数；
[0046] S170:根据所述最优控制参数获取最优控制参数成本值；
[0047] S180:根据所述最优控制参数成本值与所述建模样本集中的最小成本值的对比结 果，确定转炉炼钢工艺的最小成本值。
[0048] 上述为本发明的基于UK巧申经网络的转炉炼钢工艺成本控制方法的流程，在步骤 SllO中，实际转炉炼钢工艺过程中，为了降低成本在保证热量足够的情况下，加入废钢、铁 矿石等提高出钢量;同时通过造渣材料的加入量、入炉铁水的溫度、出钢溫度等条件的控制 实现成本的降低。为此本发明采用铁水量、废钢量、造渣材料加入量、入炉铁水的溫度、出钢 溫度、白云石加入量、石灰石加入量、铁矿石加入量、氧气消耗量、氧枪位置等作为影响成本 的控制参数。
[0049] 具体地，确定转炉炼钢过程生产成本因素，成本影响因素构成的效率观测变量集 合知….?.}。其中，日1，日2，日3为决策变量，日4~日6操作变量，而~％为其他环境变量。选取 转炉炼钢过程的成本变量为性能观测变量：{yi}。
[0050] 本发明的实施例中，选取决策变量Ql为铁水量、决策变量02为废钢量、决策变量曰3 为造渣材料加入量，04~06为溫度控制，其余环境变量包括:除硫原料量、除憐原料量、除碳 原料量，选择炉炼钢过程性能变量yi为成本量
[0051 ]其中，影响成本的控制参数如表1所示：
[0化2]
[0053] 表1参数及符号表
[0054] 在步骤S120中，样本采集;利用模拟转炉炼钢平台采集不同的控制参数下的成本， 构建建模样本集[X;Y];
[0055] (1)针对第一步在炉炼钢过程中获取的生产成本因素、性能变量的样本数据。若采 集数据不匹配，多样本按较少数量样本采集时间间隔求取平均值，从而与较少数量样本进 行匹配，得到效率影响因素样本矩阵L如下：
[0化6]
[0057]性能样本矩阵：
[005引 Y=[yii yi2 …yiN]
[0化9]其中搞为生产因素个数;N为样本个数，Qik表示第i个变量的第k个观测值，i = l， 2,...,M;k -1,2,...,No
[0060] (2)构建输入变量集合{X1，X2，X3, . . .，xm}观测样本值：
[0061]
[0062] Y=[Yi] = [yii yi2 …yiN]
[006引此时Xl~X3为决策变量，X4~X6操作变量，%为其他环境变量。
[0064]其中，采集到数据如表2所示：
[00 化 I
[0066] 表2数据采集样本部分数据
[0067] 在步骤S130中，数据预处理。利用神经网络建模过程中，其隐含层节点函数为S型 函数，其值域为[-1,1];为提高建模过程精度，故而将所有的采集的样本进行归一化处理。 即：将样本集的参量值利用线性归一化方法映射到[-1,1]范围内，得到归一化的样本集
[006引具体地，利用线性映射方法，对得到的训练输入、输出样本X、Y进行去除量纲化处 理，得到新的训练输入、输出矩阵义、7,具体归一化算法如下：
[0069]
[0070]
[0071]
[0072]
[0073]
[0074] 其中：S胃、.*1。,。为设定输入变量归一化后数据范围的最大值、最小值；
[00巧]Xik为归一化前的第i个输入变量第k个样本值；
[0076] 相为归一化后第i个输入变量第k个样本值；
[0077] xi,min=min{x化 I l《k《N}
[007引 Xi,max=max{xik I l《k《N}
[0079] 兒《、沿>:。为设定输出变量归一化后数据范围的最大值、最小值；
[0080] yA为归一化前第j个输出变量的第k个采集样本值；
[0081] .?*为归一化后第j个输出变量的第k个值；
[00的]yj,max=max{yjk| l《k《N}
[008；3] yj,min=min{yjk| l《k《N}
[0084] 得到归一化后样本集：
[0085]
[0086]
[0087]在步骤S140中，构建3层前馈神经网络，输入变量为1 =巧，输出变量 y=[.nf。设置隐含层神经元个数为SI。图3示出了根据本发明实施例的前向神经网络结构， 如图3所示，前向神经网络输入层、隐含层、输出层通过权值、阔值进行连接。设置隐含层输 出函数为S巧巧数、输m民为线忡巧数:该3民神经网络函数表达式如下：
[008引 巧
[0089] 式中函数F(X)为S型函数，为隐含层激发函数，输出层激发函数为线性函数，Wlik, w2kj，blk，b2j分别代表输入层与隐含层的连接权值，隐含层和输出层的连接权值;隐含层阔 值;输出层阔值;巧n表不归一化的样本。
[0090] 在步骤S150中，本发明利用无迹卡尔曼滤波对构建的=层前向网络参数进行训 练，得到最终数学模型。：
[0091] (1)设3层前馈神经网络，记M为输入层神经元数，Si为隐层神经元数，1为输出层神 经元数.输入层至隐层神经元的连接权值，1';4 0 = 0，1，…，= 1，2,…邱），阔值为6! (* = U.…，A)， 隐层至输出层的连接权值<(* = 〇，U，VJ = 1，2，L，〇,阔值为別/ = 1，2.L .0则UKF順中所有权 值和阔值组成的狀杰巧量I为：
[0092]
[0093] 设I中的个数为n个值。
[0094] 设定非线性方程：
[0095] (4)
[0096] 其中@函数表达式参考公式(3 )，Xk为拙寸刻的神经网络输入样本。令《 k=0，Vk=0。 化为神经网络输出样本。
[0097] (2)设定UKF计算过程中控制采样点的分布状态参数a、待选参数K，W及非负权系 数e。
[0098] (3)计算化+1个〇点（即采样点，一个采样点即为一组I值似及O点(采样点）的相应 权重。其中n为状态矩阵的I维度。A = a2(n+K)-n
[0099] 化+1个采样点计算如下：
[0100]
[0101 ]每个采样点的权值如下：
[0102]
[0104]
[0103] (4)计算O点的一步状态预测4+i|t及状态变量协方差Pk+i|k;
[0105]
[0106] '(计算后为nXn的矩阵）（5)计 算输出的一步提前预测W及协方差;
[0111] (6)进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵；[0112]
[0107；
[0108；
[0109；
[0110；
[0113]
[0114]
[0115] (7)对获取的新样本数据(Xk+i，化+1)重新进行(2)~(6)步，直至所有样本对状态矩 阵、协方差矩阵、增益矩阵进行了更新。从而得到适应于所有样本状态矩阵。
[0116] (8)对最后一组样本得到状态矩阵X作为=层网络训练得到的权值和阔值。
[0117] (9)在得到网络参数各层权值、阔值后确定，利用UKF順构建的函数模型如下：
[0118] 具体地，在步骤S150中，采用UK巧巾经网络对通过转炉模拟实验所得数据进行建 模。得到神经网络参数w、b、e，如下所示。
[0119] 得到输入层到隐层的权值wl(20X10):
[0120]
[0121]
[0122] \i、厶 f
[0123] 隐层到输出层权值w2( I X 20):
[0124] w2=[-0.46 0.17 0.12 ... -0.075]
[0125] 输出层阔值b2(lXl):
[0126] b2 = -0.6319
[0127] 因此，图4示出了所构建模型的训练样本效果W及测试样本预测精度效果图，由模 型的相对误差可知，建模效果较好，随着样本的不断训练，模型精度越来越高，符合动态建 模的特性。
[0128] 在步骤S160中，运用PSO粒子群算法对决策参量在各自的上下限进行优化，得到最 佳决策参数。图6示出了根据本发明实施例的PSO粒子群算法的函数获取极值流程示意图。 如图6所示，Sl:粒子和速度初始化;S2粒子适应值计算；S3:寻找个体极值和群体极值;S4: 速度更新和位置更新；S5:粒子适应值计算；S6:个体极值和群体极值计算；S7:满足终止条 件;S8:结束。
[0129] 其中具体步骤如下：
[0130] (I)PSO粒子群算法是基于群体的，根据对环境的适应度将群体中的个体移动到好 的区域。然而它不对个体使用演化算子，而是将每个个体看作是D维捜索空间中，有n个粒子 组成的种群X=Ui X2…Xn),在捜索空间中W-定的速度飞行，运个速度根据它本身的飞 行经验和同伴的飞行经验来动态调整。
[0131] 其中，第i个微粒表示为Xi= [Xil Xi2…XiD]T，代表第i个粒子在D维捜索空间的位 置，亦代表问题的一个潜在解。根据目标函数即可算出每个粒子位置Xi对应的适应度值，第 i个粒子的速度为Vi = [Vii Vi2…Vid]t，它经历过的最好位置(有最好的适应值)记为Pi = [Pil Pi2…PiD]T，也称为Pbest(个体极值）。在群体所有微粒经历过的最好位置的索引号用 符号旨表不，目化=山1982'。98。]了，也称为旨6631;。
[0132] (2)初始化一群微粒(群体规模为m)，包括随机的位置和速度，随机赋值。
[0133] (3)评价每个微粒的适应度；
[0134] (4)对每个微粒，将它的适应值和它经历过的最好位置Pbest的作比较，如果较好， 则将其作为当前的最好位置Pbest;
[0135] (5)对每个微粒，将它的适应值和全局所经历最好位置gbest的作比较，如果较好， 贝幢新设置gbest的索引号；
[0136] (6)在每一次迭代的过程中粒子通过个体极值和全局极值更新自身的速度和位 置，更新公式如下：
[0137]
[013 引
[0139] 式中，W为惯性权重;d=l，2,…，D;i = l，2,…，n; k为当前迭代次数;Vid为粒子的速 度;Cl和C2为非负常数，成为加速因子;ri和K为分布在[0，1]之间的随机数。为防止粒子盲 目捜索，将其位置和速度限制在一定范围内[-Xmax, Xmax]、[-Vmax, Vmax]。第一部分为微粒先前 行为的惯性，第二部分为":认知（cognition)"部分，表示微粒本身的思考;第S部分为"社 会(social)"部分，表示微粒间的信息共享与相互合作。
[0140] (7)如未达到结束条件(通常为足够好的适应值或达到一个预设最大代数Gmax)，回 到算法参数。
[0141] 具体地，利用PSO粒子群算法优化步骤S150所得神经网络求取最值。
[0142] (1)构PSO粒子群算法优化的适应度函数，采用步骤S150所得神经网络作为适应度 函数，
[0143] 设隐含层神经元的激活函数天
，神经网络的第m组输入样本苯"的输出 息为：
[0144]
[0145] 其中，恥=[Wki，Wk2,…，WkM]，w、b、峽J步骤S150所求值。
[0146] (2)设晉决策巧量的巧化区间，即Xi.min《Xi《Xi.max;表3示出了决策巧量区间值。
[0147]
[0148] 表3决策变量区间表
[0149] (3)得到最佳控制参数组合;其中，表4示出了最优参数组合。
[0150]
[0151] 表4最优参数组合
[0152] 在步骤S170和步骤S180中，将所得最优控制参数组合带入转炉模型平台中进行测 试，得到实际的控制成本值，比较最优控制参数的成本值与采集样本的最小值成本值进行 比较，若计算的最优控制成本值小于采集样本的最小成本值，则说明计算结果有效，否则重 复上述所有过程;其中，表5示出了成本的最优值和模拟值。
[ni Wl
[0154] 表5成本数据比较
[0155] 由所得优化值进行模拟炼钢实验，在模拟过程中根据计算结果取符合实际操作值 反复实验，其最优操作得到最小成本为220.76 ($/t )，说明优化所得操作参数有效，吨钢成 本减少，系统效率得到了提高。达到了降低成本的目的。说明基于ELM神经网络的转炉炼钢 工艺成本优化控制方法有效。
[0156] 与上述方法相对应，本发明还提供一种基于UK巧巾经网络的转炉炼钢工艺成本控 制系统，图2示出了根据本发明实施例的基于UK巧申经网络的转炉炼钢工艺成本控制系统逻 辑结构。
[0157] 如图2所示，本发明提供的基于UK巧巾经网络的转炉炼钢工艺成本控制系统200包 括控制参数选择单元210、建模样本集构建单元220、归一化样本集获取单元230、前馈神经 网络构建单元240、神经网络参数获取单元250、最优控制参数获取单元260、最优控制参数 成本值获取单元270和最小成本值获取单元280。
[0158] 具体地，控制参数选择单元210,用于根据转炉炼钢的工艺选择影响成本的控制参 数；
[0159] 建模样本集构建单元220,用于利用模拟转炉炼钢平台采集不同控制参数的成本， 构建建模样本集；
[0160] 归一化样本集获取单元230,用于将构建的建模样本集进行归一化处理，获得归一 化样本集；
[0161] 前馈神经网络构建单元240,用于根据所述归一化样本集构建前馈神经网络；；
[0162] 神经网络参数获取单元250,用于采用UKF算法对所述前馈神经网络进行训练建 模，获取神经网络参数；
[0163] 最优控制参数获取单元260,用于利用PSO粒子群算法对运用UKF算法所构建的模 型进行优化，获取所构建模型的最值，并根据所构建模型的最值确定最优控制参数；
[0164] 最优控制参数成本值获取单元270,用于根据所述最优控制参数获取最优控制参 数成本值；
[0165] 最小成本值获取单元280,用于根据所述最优控制参数成本值与所述建模样本集 中的最小成本值的对比结果，确定转炉炼钢工艺的最小成本值。
[0166] 其中，控制参数选择单元210的控制参数包括铁水量、废钢量、造渣材料加入量、入 炉铁水的溫度、出钢溫度、白云石加入量、石灰石加入量、铁矿石加入量、氧气消耗量、氧枪 位置。
[0167] 其中，在本发明的实施例中，前馈神经网络构建单元240在根据所述归一化样本集 构建前馈神经网络的过程中，
[0168] 所述前馈神经网络包括输入层、隐含层和输出层，所述输入层与所述隐含层由权 值连接，所述隐含层与所述输出层由权值连接;所述前馈神经网络函数表达式如下：
[0169]
[0170] 其中，函数F(X)为S型函数，为隐含层激发函数，输出层激发函数为线性函数;Wlik, w2kj，blk，b2j分别代表输入层与隐含层的连接权值，隐含层和输出层的连接权值，隐含层阔 值，输出层阔值;蜗巧表示归一化的样本。
[0171] 其中，神经网络参数获取单元250在采用UKF算法对所述前馈神经网络进行训练建 模，获取神经网络参数的过程中，
[0172] 第一步：设前馈神经网络，记M为输入层神经元数，Si为隐层神经元数，1为输出层 神经元数;输入层至隐含层神经元的连接权值，非' =〇，1，''一= 一,），阔值为6;(* = 1，2，''，〇， 隐含层至输出层的连接权值= ，A;y = l，2，L，〇,阔值为刮/ = LU，〇,则UKF神经网 络中所有权值和阔值组成的状态变量I为：
[0173]
[0174] 设I中的个数为n个值；
[0175] 设定非线性方程：
[0176]
[0177] 其中，姑为即寸刻的神经网络输入样本;令O k = 0，Vk = 0，Yk为神经网络输出样本； [017引第二步:设定UKF计算过程中控制采样点的分布状态参数a、待选参数K，W及非负 权系数e;
[0179] 第；步:计算化+1个〇点^及O点的相应权重，其中，n为状态矩阵的I维度，A = a2(n +K) -n，其中，K为待选参数，a为分布状态参数；
[0180] 第四步:计算O点的一步状态预测及状态变量协方差Pk+i|k;
[0181] 第五步:计算输出的一步提前预测^及协方差/^< ;
[0182] 第六步:进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵；
[0183] 第屯步:对获取的新样本数据重新进行第二步至第六步步，直至所有样本对状态 矩阵、协方差矩阵、增益矩阵进行了更新；
[0184] 第八步:对最后一组样本得到状态矩阵X作为前馈神经网络训练得到的权值和阔 值；
[0185] 第九步:根据获得网络参数各层权值、阔值，利用UK巧巾经网络构建的函数模型如 下：
[0186]
[0187] 其中古堪反归一函数，)?../'从)分别代表网络函数和归一化函数。
[0188] 其中，最优控制参数获取单元260在利用PSO粒子群算法对运用UKF算法所构建的 模型进行优化，获取所构建模型的最值，并根据所构建模型的最值确定最优控制参数的过 程中，
[0189] 第一步:初始化一群微粒，其中，所述一群微粒的群体规模为m，包括随机的位置和 速度，随机赋值；
[0190] 第二步:评价每个微粒的适应度；
[0191] 第S步:对每个微粒，将它的适应值和它经历过的最好位置Pbest的作比较，如果较 好，则将其作为当前的最好位置Pbest;
[0192] 第四步:对每个微粒，将它的适应值和全局所经历最好位置gbest的作比较，如果较 好，则重新设置gbest的索引号；
[0193] 第五步:在每一次迭代的过程中粒子通过个体极值和全局极值更新自身的速度和 位置，更新公式如下：
[0194]
[0195]
[0196] 其中，W为惯性权重;d=l，2,…，D;i = l，2,…，n; k为当前迭代次数;Vid为粒子的速 度;C谢C2为非负常数，为加速因子;ri和K为分布在[(U]之间的随机数；
[0197] 第六步:重复第=步到第五步的操作，直到达到一个预设最大代数Gmax,将最后一 次迭代所得个体极值和全局极值作为优化所得最佳控制参数组合。
[0198] 通过上述实施方式可W看出，本发明提供的基于UK巧申经网络的转炉炼钢工艺成 本控制方法及系统，在冶炼过程中的生产操作参数为信息载体，利用UK巧巾经网络方法挖掘 原料配方、操作参数与炼钢成本之间的关系；并通过智能优化算法利获取最低成本下的操 作参数，为实际生产最优生产提供指导，解决转炉炼钢成本较高的问题。
[0199] 如上参照附图W示例的方式描述了根据本发明提出的基于UK巧申经网络的转炉炼 钢工艺成本控制方法及系统。但是，本领域技术人员应当理解，对于上述本发明所提出的基 于UK巧申经网络的转炉炼钢工艺成本控制方法及系统，还可W在不脱离本
【发明内容】
的基础 上做出各种改进。因此，本发明的保护范围应当由所附的权利要求书的内容确定。
【主权项】
1. 一种基于UKF神经网络的转炉炼钢工艺成本控制方法，包括:根据转炉炼钢的工艺选 择影响成本的控制参数； 利用模拟转炉炼钢平台采集不同控制参数的成本，构建建模样本集； 将构建的建模样本集进行归一化处理，获得归一化样本集； 根据所述归一化样本集构建前馈神经网络； 采用UKF算法对所述前馈神经网络进行训练建模，获取神经网络参数； 利用PSO粒子群算法对运用UKF算法所构建的模型进行优化，获取所构建模型的最值， 并根据所构建模型的最值确定最优控制参数； 根据所述最优控制参数获取最优控制参数成本值； 根据所述最优控制参数成本值与所述建模样本集中的最小成本值的对比结果，确定转 炉炼钢工艺的最小成本值。2. 如权利要求1所述的基于UKF神经网络的转炉炼钢工艺成本控制方法，其中， 所述控制参数包括铁水量、废钢量、造渣材料加入量、入炉铁水的温度、出钢温度、白云 石加入量、石灰石加入量、铁矿石加入量、氧气消耗量、氧枪位置。3. 如权利要求1所述的基于UKF神经网络的转炉炼钢工艺成本控制方法，其中， 在根据所述归一化样本集构建前馈神经网络的过程中， 所述前馈神经网络包括输入层、隐含层和输出层， 所述输入层与所述隐含层由权值连接，所述隐含层与所述输出层由权值连接;所述前 馈神经网络函数表达式如下：其中，函数F(X)为S型函数，为隐含层激发函数;wlik，w2kj，blk，b2j分别代表输入层与隐 含层的连接权值，隐含层和输出层的连接权值，隐含层阈值，输出层阈值;[?, K]表不归一化 的样本。4. 如权利要求3所述的基于UKF神经网络的转炉炼钢工艺成本控制方法，其中， 在采用UKF算法对所述前馈神经网络进行训练建模，获取神经网络参数的过程中， 第一步:设前馈神经网络，记M为输入层神经元数，S1为隐层神经元数，1为输出层神经元 数;输入层至隐含层神经元的连接权值<(| = 〇，1，一，4 = 1.2^4)，阈值为以〃 =1二···，^， 隐含层至输出层的连接权值= ,y = U，L，/)，阈值为&)(./ = l，2,L ,/)，则UKF神经网 络中所有权值和阈值组成的状态变量I为：其中，Xk为K时刻的神经网络输入样本;令ω k=ο，Vk = ο，Yk为神经网络输出样本； 第二步:设定UKF计算过程中控制采样点的分布状态参数a、待选参数κ，以及非负权系 数β; 第三步:计算2η+1个〇点以及σ点的相应权重，其中，η为状态矩阵的I维度，A = a2(n+K)-n，其中，κ为待选参数，a为分布状态参数;第四步:计算σ点的一步状态预测及状态变 量协方差Pk+l|k; 第五步:计算输出的一步提前预测以及协方差; 第六步:进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵； 第七步:对获取的新样本数据重新进行第二步至第六步步，直至所有样本对状态矩阵、 协方差矩阵、增益矩阵进行了更新； 第八步:对最后一组样本得到状态矩阵X作为前馈神经网络训练得到的权值和阈值； 第九步:根据获得网络参数各层权值、阈值，利用UKF神经网络构建的函数模型如下：其中，是反归一函数，}A')分别代表网络函数和归一化函数。5. 如权利要求1所述的基于UKF神经网络的转炉炼钢工艺成本控制方法，其中， 在利用PSO粒子群算法对运用UKF算法所构建的模型进行优化，获取所构建模型的最 值，并根据所构建模型的最值确定最优控制参数的过程中， 第一步：初始化一群微粒，其中，所述一群微粒的群体规模为m，包括随机的位置和速 度，随机赋值； 第二步:评价每个微粒的适应度； 第三步:对每个微粒，将它的适应值和它经历过的最好位置Pbest的作比较，如果较好，则 将其作为当前的最好位置Pht; 第四步:对每个微粒，将它的适应值和全局所经历最好位置gbest的作比较，如果较好，则 重新设置gbest的索引号； 第五步：在每一次迭代的过程中粒子通过个体极值和全局极值更新自身的速度和位 置，更新公式如下：其中，w为惯性权重;d = l, 2,…，D;i = l, 2,…，n; k为当前迭代次数;Vid为粒子的速度； Cl和C2为非负常数，为加速因子;ri和r2为分布在[0,1]之间的随机数； 第六步:重复第三步到第五步的操作，直到达到一个预设最大代数Gmax，将最后一次迭 代所得个体极值和全局极值作为优化所得最佳控制参数组合。6. -种基于UKF神经网络的转炉炼钢工艺成本控制系统，包括： 控制参数选择单元，用于根据转炉炼钢的工艺选择影响成本的控制参数； 建模样本集构建单元，用于利用模拟转炉炼钢平台采集不同控制参数的成本，构建建 模样本集； 归一化样本集获取单元，用于将构建的建模样本集进行归一化处理，获得归一化样本 集； 前馈神经网络构建单元，用于根据所述归一化样本集构建前馈神经网络； 神经网络参数获取单元，用于采用UKF算法对所述前馈神经网络进行训练建模，获取神 经网络参数； 最优控制参数获取单元，用于利用PSO粒子群算法对运用UKF算法所构建的模型进行优 化，获取所构建模型的最值，并根据所构建模型的最值确定最优控制参数； 最优控制参数成本值获取单元，用于根据所述最优控制参数获取最优控制参数成本 值； 最小成本值获取单元，用于根据所述最优控制参数成本值与所述建模样本集中的最小 成本值的对比结果，确定转炉炼钢工艺的最小成本值。7. 如权利要求6所述的基于UKF神经网络的转炉炼钢工艺成本控制系统，其中， 所述控制参数选择单元的所述控制参数包括铁水量、废钢量、造渣材料加入量、入炉铁 水的温度、出钢温度、白云石加入量、石灰石加入量、铁矿石加入量、氧气消耗量、氧枪位置。8. 如权利要求6所述的基于UKF神经网络的转炉炼钢工艺成本控制系统，其中， 所述前馈神经网络构建单元在根据所述归一化样本集构建前馈神经网络的过程中， 所述前馈神经网络包括输入层、隐含层和输出层， 所述输入层与所述隐含层由权值连接，所述隐含层与所述输出层由权值连接;所述前 馈神经网络函数表达式如下：其中，函数F(X)为S型函数，为隐含层激发函数;wlik，w2kj，blk，b2j分别代表输入层与隐 含层的连接权值，隐含层和输出层的连接权值，隐含层阈值，输出层阈值;[I； f]表示归一化 的样本。9. 如权利要求8所述的基于UKF神经网络的转炉炼钢工艺成本控制系统，其中， 神经网络参数获取单元在采用UKF算法对所述前馈神经网络进行训练建模，获取神经 网络参数的过程中， 第一步:设前馈神经网络，记M为输入层神经元数，S1为隐层神经元数，1为输出层神经元 数;输入层至隐含层神经元的连接权值4卜〇山一，财;"1，2，.^)，阈值为片(1 = 1，2，，.4)， 隐含层至输出层的连接权值<0 = 〇，以Λ;·/ = 1，2，1·，0,阈值为MC/ = 1，2，L，0,则UKF神经网 络中所有权倌和阖倌组成的状杰夺量I为： 设I中的个数为η个值；设定非线性方程：其中，Xk为K时刻的神经网络输入样本;令ω k=ο，Vk = ο，Yk为神经网络输出样本； 第二步:设定UKF计算过程中控制采样点的分布状态参数a、待选参数κ，以及非负权系 数β; 第三步:计算2η+1个〇点以及σ点的相应权重，其中，η为状态矩阵的I维度，A = a2(n+K)-n，其中，κ为待选参数，a为分布状态参数； 第四步:计算σ点的一步状态预；及状态变量协方差Pk+11 k; 第五步:计算输出的一步提前预测以及协方I; 第六步:进行滤波更新获取新的状态矩阵、协方差矩阵、增益矩阵； 第七步:对获取的新样本数据重新进行第二步至第六步步，直至所有样本对状态矩阵、 协方差矩阵、增益矩阵进行了更新； 第八步:对最后一组样本得到状态矩阵X作为前馈神经网络训练得到的权值和阈值； 第九步:根据获得网络参数各层权值、阈值，利用UKF神经网络构建的函数模型如下：其中彳1是反归一函数;f，/(X)分别代表网络函数和归一化函数。 I 〇.如权利要求6所述的基于UKF神经网络的转炉炼钢工艺成本控制系统，其中， 所述最优控制参数获取单元在利用PSO粒子群算法对运用UKF算法所构建的模型进行 优化，获取所构建模型的最值，并根据所构建模型的最值确定最优控制参数的过程中， 第一步：初始化一群微粒，其中，所述一群微粒的群体规模为m，包括随机的位置和速 度，随机赋值； 第二步:评价每个微粒的适应度； 第三步:对每个微粒，将它的适应值和它经历过的最好位置Pbest的作比较，如果较好，则 将其作为当前的最好位置Pbest; 第四步:对每个微粒，将它的适应值和全局所经历最好位置gbest的作比较，如果较好，则 重新设置gbest的索引号； 第五步：在每一次迭代的过程中粒子通过个体极值和全局极值更新自身的速度和位 置，更新公式如下：其中，w为惯性权重;d = l, 2,…，D;i = l, 2,…，n; k为当前迭代次数;Vid为粒子的速度； Cl和C2为非负常数，为加速因子;ri和r2为分布在[0,1]之间的随机数； 第六步:重复第三步到第五步的操作，直到达到一个预设最大代数Gmax，将最后一次迭 代所得个体极值和全局极值作为优化所得最佳控制参数组合。
【文档编号】G05B13/04GK106019940SQ201610452088
【公开日】2016年10月12日
【申请日】2016年6月21日
【发明人】唐海红, 辜小花, 李太福, 耿讯, 张倩影, 王坎
【申请人】重庆科技学院