CN110427715B - 基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 - Google Patents
基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 Download PDFInfo
- Publication number
- CN110427715B CN110427715B CN201910728160.0A CN201910728160A CN110427715B CN 110427715 B CN110427715 B CN 110427715B CN 201910728160 A CN201910728160 A CN 201910728160A CN 110427715 B CN110427715 B CN 110427715B
- Authority
- CN
- China
- Prior art keywords
- furnace
- blast furnace
- hearth
- cooling wall
- thermal state
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Game Theory and Decision Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于时间序列和高炉多维度的炉缸热状态趋势预测方法,利用物质流运动规律、动态时间序列相关性分析法依次给出高炉冷却壁立体水温差对炉缸温度响应的滞后时间,计算泊松相关系数确定有效的历史炉温时间序列及与炉温预测相关的其它过程参数;再利用立体高炉各指标参数及其相关历史信息,建立基于时间序列和高炉多维度的PSO‑LSSVM炉缸热状态预测模型,验证高炉冷却壁立体水温差的多维性变化对炉缸热状态变化趋势预测的准确性。本发明充分利用了高炉冷却壁段水温差的立体分布特征及其依次对炉缸温度的滞后性影响,可协助高炉操作者,准确预判高炉炉缸热状态,为高炉的高效、节能运行奠定基础。
Description
技术领域
本发明属于复杂高炉炼铁过程中实时不可测关键参数预测领域,尤其是涉及一种基于时间序列和高炉多维度的炉缸热状态趋势预测的方法,高炉冷却壁立体水温差的多维性变化对炉缸热状态变化趋势的超前预知,可准确预判高炉炉缸热状态,协助高炉操作者,为高炉的高效、节能运行奠定基础。
背景技术
高炉生产运行目标是炉况平稳下的高效、节能运行。稳定的炉温是炉况平稳的先决条件,其高低直接影响炼铁过程中高炉稳定顺行、低能耗、高质量等性能指标,因而现有生产操作制度的基础就是确保炉温的稳定。高炉炼铁进行节能降耗的最有效措施是高煤比降焦,而制约其有效性的则是炉缸热状态。不可实时测得的炉温及炉缸热状态所反映的高炉运行行为热状态,是直接影响喷煤量的关键因素,故实时准确地预测炉温及炉缸热状态,进而掌握其变化趋势,对于具有大滞后、大惯性的高炉实施节能降耗优化操作至关重要。
传统的高炉炉缸热状态通常用铁水[Si]、铁水温度(MIT)或炉缸冷却壁热电偶温度表征,然而,单一参数或单一参数建立的炉温预测模型皆难以准确描述高炉炉缸热状态。据此,本发明利用高炉冷却壁段的空间立体分布特性,定义冷却壁立体水温差概念。高炉炉体冷却壁从上到下分为多段,圆周上呈立体式分布,各冷却壁段水温差的变化可反映对应炉内温度的变化。本文基于纵向物质流运动规律,通过分析高炉冷却壁立体水温差对炉缸温度响应的滞后性,提出高炉冷却壁立体水温差的多维性变化可超前反映炉缸热状态的变化趋势,从而提出了一种基于时间序列和高炉多维度的炉缸热状态趋势预测方法。
发明内容
为了解决上述技术问题,本发明的目的是提供一种基于时间序列和高炉多维度的炉缸热状态趋势预测的方法。
为解决上述技术问题,本发明采用如下技术方案:
一种基于时间序列和高炉多维度的炉缸热状态趋势预测的方法,包括如下步骤:利用纵向物质流运动规律和动态时间序列相关性分析法依次给出高炉冷却壁立体水温差对炉缸温度响应的滞后时间,然后计算泊松相关系数提取有效的历史炉温时间序列及与炉温预测相关的其它过程参数;利用立体高炉各指标参数及其相关历史信息,建立基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型,验证高炉冷却壁立体水温差的多维性变化对炉缸热状态变化趋势预测的准确性。
进一步的,具体步骤如下:
S1、采集某高炉实时运行过程数据
包括历史炉温数据、各冷却壁段历史水温差及与炉温预测相关的其它过程数据;
S2、过程数据预处理
用拉依达准则剔除所选原始样本集中的异常数据,对异常数据赋空值并用牛顿插值法修补,最终使选取的过程数据分布在(μ-3σ,μ+3σ)区间内;
S3、炉缸热状态预测相关特征参数提取
计算所选过程参数与炉温的泊松相关系数,依据相关系数的阈值及顶峰值选取炉缸热状态预测相关特征参数;
S4、建立模型
将选取的相关特征参数数据作为模型的输入,建立基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型;
S5、对比分析无水温差时和加入高炉冷却壁立体水温差时所建模型的相对误差和命中率。
进一步的,所述步骤(S1)中,高炉炉体的各冷却壁段具体为:从上到下分为多段,圆周上呈立体式分布,每段冷却壁在横向圆周上设置若干个热电偶测温点,实现多点测温;选取炉身的上、中、下、炉腰、炉腹、和炉缸六大段冷却壁段,利用高炉冷却壁纵向分段水温差的立体多维性,依次分析立体水温差变化对炉缸温度响应的滞后性。
进一步的,所述步骤S2中,剔除过程及实际各冷却壁段水温差的计算如下:
μ和σ分别为样本数据xk的平均值和标准差;
对修补后的冷却壁段所有测温点数据,计算高炉每段冷却壁横向圆周上所有测温点的平均水温差作为冷却壁段的实际水温差。
进一步的,所述步骤S3中,泊松相关系数计算、相关系数顶峰值及阈值设定过程如下:
泊松相关系数计算过程描述为
选取相关系数最大的尖峰点处所对应的各冷却壁段水温差的历史时间序列作为高炉冷却壁立体水温差对炉缸温度响应的滞后时间;选择相关系数|r|≥0.1的历史炉温时间序列、与炉温预测相关的其它过程参数作为炉缸热状态预测相关特征参数。
进一步的,所述步骤S4中,建模过程如下:
S4-1、将预处理后的指标数据分为训练集和测试集,将训练数据输送到LSSVM预测模型中训练,同时采用粒子群算法对LSSVM的惩罚因子γ和核参数宽度σ'进行优化,最终得到用于炉缸热状态预测的PSO-LSSVM模型;
S4-2、利用训练好的PSO-LSSVM模型对测试集进行预测。
进一步的,所述步骤S4-1中,采用PSO算法对LSSVM中的惩罚因子γ和核参数宽度σ'进行优化的具体过程为:
S4-1-1、随机产生一个初始种群,种群中每个粒子都有自己的速度和位置,初始的速度和位置随机生成;
S4-1-3、粒子通过跟踪两个“极值”:个体极值Pibest和全局极值PGibest,更新粒子的当前速度和位置;
S4-1-4、反复迭代,当适应度满足预定精度或达到最大迭代次数时停止迭代;
S4-1-5、记录此时对应的γ和σ'即为优化结果;
其中所述步骤S4-1-3对粒子的速度、位置的更新过程如下:
速度调整:
Vi(t+1)=wVi(t)+c1rand(t)(Pibest-Xi(t))+c2rand(t)(PGibest-Xi(t))
位置调整:
Xi(t+1)=Xi(t)+Vi(t+1)
式中,Vi和Xi分别为粒子的当前速度和位置;c1、c2为学习因子,且c1、c2∈[0,2];rand(t)为(0,1)之间的随机函数;w为惯性权重;Pibest、PGibest分别为粒子和粒子群的最优位置;搜索中,通过将Xi对应的Pibest、PGibest的适应度值比较和位置调整方程不断调整最优的Pibest和PGibest;
其中所述步骤S4-1-3中更新Pibest和PGibest的过程为:
①将各粒子的当前适应度值F(Xi)与该粒子自身的最优适应度值F(Pibest)比较,若F(Xi)<F(Pibest),调整F(Pibest)=F(Xi),并将粒子的当前位置作为该粒子的最优位置;
②将各粒子的F(Pibest)与粒子群的最优位置适应度值F(Pibest)比较,若F(Pibest)<F(PGibest),调整F(Pibest)=F(Xi),将此粒子的最优位置作为所有粒子的最优位置;
③由粒子群优化算法的进化方程对粒子的速度、位置进化调整,得到新的粒子位置;判断所有粒子最优位置的适应度值或迭代次数满足要求,结束计算并保存当前粒子群整体的最优位置值,若不满足则返回流程②继续计算;
④经参数优化调整,将最终粒子群整体最优位置所对应的参数作为最小二乘支持向量机建模所需的参数,建立PSO-LSSVM模型对测试样本进行预测;
其中所述步骤S4-1-3中,惯性权重设置如下:
较大的惯性权重w有较好的全局搜索能力,反之有较强的局部搜索能力,随迭代次数的增加,惯性权重应不断减少,粒子群算法在初期应具有较强的全局收敛能力,而晚期有较强的局部收敛能力;取w随迭代次数T增加线性递减的效果最好,w与T的关系表达式为
式中,wmax为初始权重,wmin为最终权重,Tmax为最大迭代次数。
进一步的,所述步骤S5中,模型评判标准如下:
相对误差(MSE)定义为
命中率定义为
式中,l为样本总数,ε为模型输出的允许误差。
与现有技术相比,本发明的有益技术效果:
本发明利用高炉冷却壁段水温差的立体分布特征,对横向圆周上的每段冷却壁求所有测温点的平均水温差,以此作为冷却壁段实际水温差,纵向利用物质流运动规律和动态时间序列相关性分析法,依次分析高炉冷却壁立体水温差的多维性变化对炉缸热状态响应的滞后性,结合历史炉温信息及与炉温预测相关的其它过程参数对炉缸热状态的影响,建立了基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型,并以某大型高炉为研究对象,给出了基于实际运行数据的预测与分析,验证高炉冷却壁立体水温差的多维性变化对炉缸热状态变化趋势预测的准确性。
附图说明
下面结合附图说明对本发明作进一步说明。
图1为本发明中基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型建立流程示意图。
图2为本发明中粒子寻优算法的流程示意图。
图3为本发明实施例中高炉结构示意图。
图4为本发明实施例中高炉冷却壁立体水温差概率密度统计分布直方图。
图5为本发明实施例中高炉冷却壁立体水温差参数与炉温的相关系数矩阵可视化图。
图6为本发明实施例中炉缸热状态预测相关特征参数提取三维可视化图。
图7为本发明实施例中基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型跟踪效果图。
图8为本发明实施例中加入立体水温差的PSO-LSSVM模型粒子群收敛过程图。
具体实施方式
下面结合以上附图及实施例对本发明进一步详述,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
建模流程如图1所示,具体步骤如下:
S1、采集某高炉实时运行过程数据
包括历史炉温数据、各冷却壁段历史水温差及与炉温预测相关的其它过程数据;
S2、过程数据预处理
用拉依达准则(3σ准则)剔除所选原始样本集中的异常数据(规定突跳点数据与均值的偏差大于3σ),对异常数据赋空值并用牛顿插值法修补,最终使选取的过程数据分布在(μ-3σ,μ+3σ)区间内;
S3、炉缸热状态预测相关特征参数提取
计算所选过程参数与炉温的泊松相关系数,依据相关系数的顶峰值及阈值选取炉缸热状态预测相关特征参数(包括高炉冷却壁立体水温差时间序列、历史炉温序列及与炉温预测相关的其它过程参数的选取);
S4、建立模型
将选取的相关特征参数数据作为模型的输入,建立基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型;
S5、对比分析无水温差时和加入高炉冷却壁立体水温差时所建模型的相对误差(MSE)和命中率(H)。
其中所述步骤S2中,剔除过程及实际各冷却壁段水温差的计算如下:
μ和σ分别为样本数据xk的平均值和标准差;
对修补后的冷却壁段所有测温点数据,计算高炉每段冷却壁横向圆周上所有测温点的平均水温差作为冷却壁段的实际水温差。
其中所述步骤(S3)中,泊松相关系数计算、相关系数顶峰值及阈值设定过程如下:
泊松相关系数计算过程描述为
选取相关系数最大的尖峰点处所对应的各冷却壁段水温差的历史时间序列作为高炉冷却壁立体水温差对炉缸温度响应的滞后时间。选择相关系数|r|≥0.1的历史炉温时间序列、与炉温预测相关的其它过程参数作为炉缸热状态预测相关特征参数。
其中所述步骤(S4)中,建模过程如下:
S4-1、将预处理后的指标数据分为训练集和测试集,将训练数据输送到LSSVM预测模型中训练,同时采用粒子群算法对LSSVM的惩罚因子γ和核参数宽度σ'进行优化,最终得到用于炉缸热状态预测的PSO-LSSVM模型;
S4-2、利用训练好的PSO-LSSVM模型对测试集进行预测。
其中所述步骤S4-1中,如图2所示为粒子寻优算法的流程图,采用PSO算法对LSSVM中的惩罚因子γ和核参数宽度σ'进行优化的具体过程为:
S4-1-1、随机产生一个初始种群,种群中每个粒子都有自己的速度和位置,初始的速度和位置随机生成;
S4-1-3、粒子通过跟踪两个“极值”(个体极值Pibest和全局极值PGibest)更新粒子的当前速度和位置;
S4-1-4、反复迭代,当适应度满足预定精度或达到最大迭代次数时停止迭代;
S4-1-5、记录此时对应的γ和σ'即为优化结果。
其中所述步骤S4-1-2对粒子的速度、位置的更新过程如下:
速度调整:
Vi(t+1)=wVi(t)+c1rand(t)(Pibest-Xi(t))+c2rand(t)(PGibest-Xi(t))
位置调整:
Xi(t+1)=Xi(t)+Vi(t+1)
式中,Vi和Xi分别为粒子的当前速度和位置;c1、c2为学习因子,且c1、c2∈[0,2];rand(t)为(0,1)之间的随机函数;w为惯性权重;Pibest、PGibest分别为粒子和粒子群的最优位置。搜索中,通过将Xi对应的Pibest、PGibest的适应度值比较和位置调整方程不断调整最优的Pibest和PGibest。
其中所述步骤S4-1-3中更新Pibest和PGibest的过程为:
①将各粒子的当前适应度值F(Xi)与该粒子自身的最优适应度值F(Pibest)比较,若F(Xi)<F(Pibest),调整F(Pibest)=F(Xi),并将粒子的当前位置作为该粒子的最优位置;
②将各粒子的F(Pibest)与粒子群的最优位置适应度值F(Pibest)比较,若F(Pibest)<F(PGibest),调整F(Pibest)=F(Xi),将此粒子的最优位置作为所有粒子的最优位置;
③由粒子群优化算法的进化方程对粒子的速度、位置进化调整,得到新的粒子位置(同时得到新的最小二乘支持向量机参数)。判断所有粒子最优位置的适应度值或迭代次数满足要求,结束计算并保存当前粒子群整体的最优位置值,若不满足则返回流程②继续计算;
④经参数优化调整,将最终粒子群整体最优位置所对应的参数作为最小二乘支持向量机建模所需的参数,建立PSO-LSSVM模型对测试样本进行预测。
其中所述步骤S4-1-3中,惯性权重设置如下:
较大的惯性权重w有较好的全局搜索能力,反之有较强的局部搜索能力。随迭代次数的增加,惯性权重应不断减少,粒子群算法在初期应具有较强的全局收敛能力,而晚期有较强的局部收敛能力。本发明经多次尝试,取w随迭代次数T增加线性递减的效果最好,w与T的关系表达式为
式中,wmax为初始权重,wmin为最终权重,Tmax为最大迭代次数。
其中所述步骤S5中,模型评判标准如下:
相对误差(MSE)定义为
命中率定义为
式中,l为样本总数,ε为模型输出的允许误差。
实施例:
高炉冷却壁空间立体分布描述:本发明以某高炉(容积为4150m3)为例,图3为高炉内部结构示意图。高炉自上而下纵向间隔设置18段冷却壁段,每段冷却壁呈横向圆周设置若干个测温点(热电偶测温,每段温度检测段包含44~58个测温点)实现多点测温,图3右侧为炉身中X4(第12段)冷却壁水温差监测点横向安置剖面简图。本发明选取炉身X4(H6:16~18段(上)、H5:12~15段(中)、H4:9~11段(下))、炉腰X3(H3:8段)、炉腹X2(H2:7段)、炉缸X1(H1:1~5段)六大段冷却壁分析,由于风口区的冷却壁(第6段)水温差对炉长的参考价值较小,本发明将忽略该部分的影响。
本发明选取某高炉不同工况下的生产数据386组,其中选取训练样本240组,剩余的146组数作为测试样本。指标参数有风量(BV)、压差(*P)、风温(BT)、富氧率(XO2)、燃料比(FR)、透气性指数(K)、顶温(TT)和铁水[P](%)、铁水[Mn](%)和MIT(℃)、铁水[Si](%)、冷却壁段H1、H2、H3、H4、H5、H6水温差的k时刻及其历史时刻数据。
本发明利用3σ准则剔除所选原始样本集中的异常数据,对修补后的高炉冷却壁段所有测温点数据,计算高炉每段冷却壁横向圆周上所有测温点的平均水温差作为冷却壁段的实际水温差。图4为高炉冷却壁立体水温差概率密度统计分布直方图,由于高炉长期运行于稳定状态,各参数近似服从以预期工作点为中心的正态分布。高炉运行过程中,炉身中部的冷却壁段水温差最高,炉身下到炉缸部分逐段下降,且从炉腹部位开始下降幅度增大。
本发明通过计算各指标参数与炉温的泊松相关系数,选取相关系数最大的尖峰点处所对应的各冷却壁段水温差的历史时间序列,以此作为冷却壁立体水温差对炉缸温度响应的滞后时间,选择相关系数|r|≥0.1的历史炉温时间序列及与炉温预测相关的其它过程参数作为炉缸热状态预测相关特征参数。
图5为冷却壁立体水温差参数与炉温的相关系数矩阵可视化图,纵轴代表高炉冷却壁段H1、H2、H3、H4、H5、H6水温差,图中横坐标分别为高炉铁水[Si]和MIT的历史时间序列,不同颜色的方块为相关系数的可视化,右侧的条形颜色图代表相关系数对应的颜色范围,相关系数为正且越大时,方块儿颜色趋于红色系,反之,则趋于蓝色系。相关系数的绝对值越大,方块面积越大,反之越小。
水温差数据采样周期为1小时,如图5所示,H6、H5、H4、H3、H2冷却壁段水温差对铁水[Si]和MIT的滞后时间依次为6、5、5、3、3小时和8、6、6、5、3小时,H1与高炉铁水[Si]和MIT当前时刻值均呈显性相关。综上分析,高炉冷却壁立体水温差参数对炉缸温度响应的滞后时间在7~8小时内,刚好与该高炉冶炼周期相同,表明在此期间,高炉下降料速平稳,高炉处于规则、稳定的运行状态。此时选取高炉铁水[Si]和MIT预测模型的高炉冷却壁立体水温差输入变量分别为28个和34个。
图6为炉缸热状态预测相关特征参数提取三维可视化图,数据采样间隔为1炉次(1~2小时),本发明中设定相关系数阈值为0.1。图6(a)验证了前四炉炉温对当前炉次炉温存在影响,包含了6~8小时的历史炉温信息即一个冶炼周期,可涵盖历史炉温对实际炉温的影响。图6(b)中,铁水[Si]和MIT对应横坐标段(k-3~k)分别为MIT历史时刻温度和铁水[Si]历史时刻值。综上所述,由于阈值限定,选取高炉铁水[Si]和MIT的过程输入变量分别为15个和16个。
综上分析,无水温差时,高炉铁水[Si]和MIT预测模型输入变量分别为15和16个,加入高炉冷却壁立体水温差时,高炉铁水[Si]和MIT预测模型输入变量分别为43和50个。设置粒群最大迭代次数Tmax=200,c1=1.5,c2=1.7,初始权重wmax=0.9,最终权重wmin=0.45,取铁水[Si]和MIT预测允许误差ε分别为(±0.1%)和(±17℃)。
图7(a~b)为无水温差的PSO-LSSVM炉缸热状态预测模型(NW PSO-LSSVM)和加入高炉冷却壁立体水温差的PSO-LSSVM炉缸热状态预测模型(WPSO-LSSVM)的仿真效果图。NWPSO-LSSVM模型和WPSO-LSSVM模型中,铁水[Si]和MIT预测模型的命中率分别为84.93、86.33%和97.58%、97.24%。图7(c~d)为WPSO-LSSVM预测模型的误差曲线图,预测误差基本在允许误差内。显然,WPSO-LSSVM炉缸热状态预测模型较传统的NW PSO-LSSVM模型预测效果更为理想。图8为WPSO-LSSVM模型参数的优化过程,此时铁水[Si]和MIT的种群数分别为6和10,适应度值随迭代次数的增加逐渐收敛到0.0012和0.8142,最终搜索得到的优化参数分别为:
铁水[Si]:k=0.27,σ'=0.1,γ=93.0
MIT:k=0.27,σ'=0.1,γ=87.5
本发明突破传统的单一参数或单一参数建立炉温预测模型表征炉缸热状态的局限性,利用高炉冷却壁立体水温差的多维性变化对炉缸热状态变化趋势预知的准确性,建立了基于时间序列和高炉多维度的PSO-LSSVM高炉炉缸热状态预测模型。模型的仿真结果表明,加入高炉冷却壁立体水温差的预测模型更能准确地预测炉温。高炉冷却壁立体水温差的多维性变化可超前预知炉缸热状态变化趋势,可协助高炉操作者,更好地使炼铁过程参数平稳运行于优化目标的小范围内,对高炉实现节能降耗具有重大意义。
以上所述的实施例仅是对本发明的优选方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (7)
1.一种基于时间序列和高炉多维度的炉缸热状态趋势预测的方法,其特征在于,具体步骤如下:
S1、采集某高炉实时运行过程数据
包括历史炉温数据、各冷却壁段历史水温差数据;
S2、过程数据预处理
用拉依达准则剔除所选原始样本集中的异常数据,对异常数据赋空值并用牛顿插值法修补,最终使选取的过程数据分布在(μ-3σ,μ+3σ)区间内,其中μ和σ分别为样本数据的平均值和标准差;
S3、炉缸热状态预测相关特征参数提取
计算所选过程参数与炉温的泊松相关系数,依据相关系数的阈值和顶峰值选取炉缸热状态预测相关特征参数;
S4、建立模型
将选取的相关特征参数数据作为模型的输入,建立基于时间序列和高炉多维度的PSO-LSSVM炉缸热状态预测模型;
S5、对比分析无水温差时和加入高炉冷却壁立体水温差时所建模型的相对误差和命中率。
2.根据权利要求1所述的基于时间序列和高炉多维度的炉缸热状态趋势预测的方法,其特征在于,所述步骤S1中,高炉炉体的各冷却壁段具体为:从上到下分为多段,圆周上呈立体式分布,每段冷却壁在横向圆周上设置若干个热电偶测温点,实现多点测温;选取炉身的上、中、下、炉腰、炉腹、和炉缸六大段冷却壁段,利用高炉冷却壁纵向分段水温差的立体多维性,依次分析立体水温差变化对炉缸温度响应的滞后性。
5.根据权利要求1所述的基于时间序列和高炉多维度的炉缸热状态趋势预测的方法,其特征在于,所述步骤S4中,建模过程如下:
S4-1、将预处理后的指标数据分为训练集和测试集,将训练数据输送到LSSVM预测模型中训练,同时采用粒子群算法对LSSVM的惩罚因子γ和核参数宽度σ'进行优化,最终得到用于炉缸热状态预测的PSO-LSSVM模型;
S4-2、利用训练好的PSO-LSSVM模型对测试集进行预测。
6.根据权利要求5所述的基于时间序列和高炉多维度的炉缸热状态趋势预测的方法,其特征在于,所述步骤S4-1中,采用PSO算法对LSSVM中的惩罚因子γ和核参数宽度σ'进行优化的具体过程为:
S4-1-1、随机产生一个初始种群,种群中每个粒子都有自己的速度和位置,初始的速度和位置随机生成;
S4-1-3、粒子通过跟踪两个“极值”:个体极值Pibest和全局极值PGibest,更新粒子的当前速度和位置;
S4-1-4、反复迭代,当适应度满足预定精度或达到最大迭代次数时停止迭代;
S4-1-5、记录此时对应的γ和σ'即为优化结果;
所述步骤S4-1-3对粒子的速度、位置的更新过程如下:
速度调整:
Vi(t+1)=wVi(t)+c1rand(t)(Pibest-Xi(t))+c2rand(t)(PGibest-Xi(t))
位置调整:
Xi(t+1)=Xi(t)+Vi(t+1)
式中,Vi和Xi分别为粒子的当前速度和位置;c1、c2为学习因子,且c1、c2∈[0,2];rand(t)为(0,1)之间的随机函数;w为惯性权重;Pibest、PGibest分别为粒子和粒子群的最优位置;搜索中,通过将Xi对应的Pibest、PGibest的适应度值比较和位置调整方程不断调整最优的Pibest和PGibest;
所述步骤S4-1-3中更新Pibest和PGibest的过程为:
①将各粒子的当前适应度值F(Xi)与该粒子自身的最优适应度值F(Pibest)比较,若F(Xi)<F(Pibest),调整F(Pibest)=F(Xi),并将粒子的当前位置作为该粒子的最优位置;
②将各粒子的F(Pibest)与粒子群的最优位置适应度值F(Pibest)比较,若F(Pibest)<F(PGibest),调整F(Pibest)=F(Xi),将此粒子的最优位置作为所有粒子的最优位置;
③由粒子群优化算法的进化方程对粒子的速度、位置进化调整,得到新的粒子位置;判断所有粒子最优位置的适应度值或迭代次数满足要求,结束计算并保存当前粒子群整体的最优位置值,若不满足则返回流程②继续计算;
④经参数优化调整,将最终粒子群整体最优位置所对应的参数作为最小二乘支持向量机建模所需的参数,建立PSO-LSSVM模型对测试样本进行预测;
所述步骤S4-1-3中,惯性权重设置如下:
较大的惯性权重w有较好的全局搜索能力,反之有较强的局部搜索能力,随迭代次数的增加,惯性权重应不断减少,粒子群算法在初期应具有较强的全局收敛能力,而晚期有较强的局部收敛能力;取w随迭代次数T增加线性递减的效果最好,w与T的关系表达式为
式中,wmax为初始权重,wmin为最终权重,Tmax为最大迭代次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910728160.0A CN110427715B (zh) | 2019-08-08 | 2019-08-08 | 基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910728160.0A CN110427715B (zh) | 2019-08-08 | 2019-08-08 | 基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110427715A CN110427715A (zh) | 2019-11-08 |
CN110427715B true CN110427715B (zh) | 2022-10-14 |
Family
ID=68414900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910728160.0A Active CN110427715B (zh) | 2019-08-08 | 2019-08-08 | 基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110427715B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111695735B (zh) * | 2020-06-12 | 2023-06-16 | 成都国铁电气设备有限公司 | 一种基于流计算的铁路弓网实时预警方法、系统及装置 |
CN113849333B (zh) * | 2021-09-26 | 2023-07-14 | 中国地质大学(武汉) | 一种基于WN-Spline基线域算法的数据中心自驱排障方法及系统 |
CN115034370B (zh) * | 2022-06-24 | 2024-02-27 | 山东钢铁股份有限公司 | 基于bp网络模型预测高炉炉缸活性的方法 |
CN115584371A (zh) * | 2022-10-26 | 2023-01-10 | 河钢股份有限公司唐山分公司 | 高炉炉温闭环控制系统的设计方法及应用 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663412A (zh) * | 2012-02-27 | 2012-09-12 | 浙江大学 | 基于最小二乘支持向量机的电力设备载流故障趋势预测方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4948304B2 (ja) * | 2006-11-13 | 2012-06-06 | 株式会社神戸製鋼所 | 高炉の溶銑温度予測方法 |
KR101246513B1 (ko) * | 2011-08-30 | 2013-03-25 | 현대제철 주식회사 | 고로의 노열 변동 예측방법 |
CN103472865B (zh) * | 2013-09-22 | 2015-09-30 | 浙江大学 | 智能最小二乘的农药废液焚烧炉炉温最佳化系统及方法 |
CN104267610B (zh) * | 2014-08-29 | 2017-05-17 | 内蒙古科技大学 | 高精度的高炉冶炼过程异常数据检测及修补方法 |
CN105605610A (zh) * | 2016-01-13 | 2016-05-25 | 上海交通大学 | 一种基于粒子群算法的锅炉燃烧优化方法 |
CN107133460A (zh) * | 2017-04-26 | 2017-09-05 | 中国能源建设集团广东省电力设计研究院有限公司 | 一种锅炉飞灰含碳量在线动态预测方法 |
CN107119159B (zh) * | 2017-06-14 | 2019-02-22 | 内蒙古科技大学 | 一种存在整数规划问题的高炉布料过程布料矩阵的优化方法 |
CN109146141B (zh) * | 2018-07-25 | 2021-11-09 | 上海交通大学 | 一种燃煤电站回转式空气预热器漏风率预测方法 |
-
2019
- 2019-08-08 CN CN201910728160.0A patent/CN110427715B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663412A (zh) * | 2012-02-27 | 2012-09-12 | 浙江大学 | 基于最小二乘支持向量机的电力设备载流故障趋势预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110427715A (zh) | 2019-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110427715B (zh) | 基于时间序列和高炉多维度的炉缸热状态趋势预测的方法 | |
CN109063358B (zh) | 一种新型高炉顺行评价的方法 | |
CN103544273A (zh) | 利用模式识别技术对炉况整体状态进行评估的方法 | |
CN111915080B (zh) | 一种基于铁水质量约束的原燃料成本最优配比方法 | |
CN106521059B (zh) | 用相控阵雷达测量高炉料面矿焦比来控制高炉气流分布的方法 | |
CN101109950A (zh) | 高炉生产过程控制信息智能化系统 | |
CN103631140B (zh) | 基于性能评估的焦炉加热燃烧过程火道温度自动调节方法 | |
CN110400009A (zh) | 基于自适应遗传算法的高炉炼铁多目标智能优化方法 | |
CN112066355B (zh) | 基于数据驱动的余热锅炉阀门自适应调节方法 | |
Li et al. | Dynamic time features expanding and extracting method for prediction model of sintering process quality index | |
CN104615856B (zh) | 基于热风炉组的煤气消耗量预测模型建立方法及装置 | |
CN107622279A (zh) | 高炉内部状态的分类方法 | |
CN116502769A (zh) | 高炉炉况的评分方法与系统 | |
CN116127760A (zh) | 基于混合建模与改进沙丘猫群算法的加热炉炉温优化方法 | |
CN112359159A (zh) | 一种基于深度强化学习的热风炉自动烧炉方法及系统 | |
CN110373508B (zh) | 软熔带识别方法、装置及电子设备 | |
CN112329269B (zh) | 一种基于工况识别的烧结点火温度建模预测方法 | |
WO2023070694A1 (zh) | 一种判定高炉最优气流分布的方法、电子设备和存储介质 | |
CN115456264B (zh) | 一种中小型转炉的终点碳含量和终点温度预测方法 | |
CN103160629B (zh) | 一种预报高炉趋热的方法 | |
CN115186900A (zh) | 适用于多工况类别的动态高炉煤气产生量预测方法及系统 | |
CN113408186A (zh) | 一种冶金焦炭质量预测影响因素确定方法 | |
CN103160626B (zh) | 一种判断高炉炉缸过凉的方法 | |
CN116127345B (zh) | 基于深度聚类生成对抗网络的转炉炼钢工艺模式设计方法 | |
CN111679584B (zh) | 一种高炉冶炼的调控方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |