CN103439889B - 基于排放烟气分析的烧结终点预测方法 - Google Patents
基于排放烟气分析的烧结终点预测方法 Download PDFInfo
- Publication number
- CN103439889B CN103439889B CN201310381670.8A CN201310381670A CN103439889B CN 103439889 B CN103439889 B CN 103439889B CN 201310381670 A CN201310381670 A CN 201310381670A CN 103439889 B CN103439889 B CN 103439889B
- Authority
- CN
- China
- Prior art keywords
- end point
- sintering
- flue gas
- sintering end
- data
- 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
Landscapes
- Manufacture And Refinement Of Metals (AREA)
Abstract
本发明提供了基于排放烟气分析的烧结终点预测方法。该方法采用在线烟气分析仪实时分析烧结过程中CO和CO2的浓度,根据在线烟气分析仪监测获得的CO和CO2的浓度,同时监测烧结终点的位置数据,经过构建数据库、数据预处理及建立烧结终点预测模型最终实现烧结终点的预测。本发明可以准确预测烧结终点,提高烧结生产过程的控制精度,极大地减弱烧结生产过程控制的“大滞后性”,从而实现生产过程的及时、准确调整,进而达到提高***效率和过程质量的作用。
Description
技术领域
本发明涉及一种烧结终点预测方法,尤其是一种分析排放烟气成分实现的烧结终点方法。
背景技术
烧结是将各种粉状含铁原料,配入适量的燃料和熔剂,加入适量的水,经混合和造球后在烧结设备上使物料发生一系列物理化学变化,将矿粉颗粒黏结成块的过程。在烧结过程中,烧结终点是烧结结束时的位置,是判断烧结过程的重要参数之一。烧结终点预测是一个多变量、非线性、大时滞问题,传统的判断方法为风箱烟气温度判断法、负压法及理论计算法,但都存在不同程度的偏差和不足。
烧结排放烟气是烧结过程中产生的气体和烧结烟气的残留,主要成分为O2、氮氧化物(NOX)、CO、CO2、硫氧化物(SOX)以及水汽等。排放烟气的组分浓度可以反映烧结过程情况,根据烧结烟气的浓度变化情况来与烧结终点关联,则有望基于烧结烟气的监测实现烧结终点的预测。
发明内容
本发明所要解决的技术问题是:提供一种基于排放烟气分析的烧结终点预测方法,对烧结排放烟气进行实时检测,实现烧结终点的准确预测。
本发明通过实时检测排放烟气中COx的浓度,建立烧结过程中烧结终点的预测模型,根据预测结果优化烧结工艺,提高烧结产品质量,降低烧结能耗,控制烟气排放。
为了解决上述技术问题,本发明提供的技术方案是,基于排放烟气分析的烧结终点预测方法,该方法采用在线烟气分析仪实时分析烧结过程中CO和CO2的浓度,同时监测烧结终点的位置数据,并将上述监测值传送给计算机,烧结终点预测模型安装在计算机中,根据在线烟气分析仪发送的CO和CO2的浓度实现烧结终点的预测,基于排放烟气分析的烧结终点预测方法包括如下步骤:
(1)构建烧结烟气浓度-烧结终点数据库
实时监测、采集烧结烟气CO和CO2的浓度和烧结终点位置,构建烧结烟气CO和CO2浓度-烧结终点数据库;
(2)数据预处理
对数据库中的数据进行异常数据的剔除;
(3)建立烧结终点预测模型
根据烧结烟气中CO与CO2浓度值,采用核偏最小二乘法建立烧结终点预测模型。
前述的基于排放烟气分析的烧结终点预测方法,当烧结过程完成后,将监测、采集的排放烟气CO和CO2的浓度和烧结终点位置纳入数据库中。
前述的基于排放烟气分析的烧结终点预测方法,剔除异常数据采用两种方法:1)对数据库中机速为<0.9m/min的数据直接剔除;和2)将数据不完整、缺失的数据进行直接剔除。
前述的基于排放烟气分析的烧结终点预测方法,所述核偏最小二乘法建立烧结终点预测模型是将数据库中数据映射到高维特征空间,然后在该高维特征空间中执行偏最小二乘。
前述的基于排放烟气分析的烧结终点预测方法,所述核偏最小二乘法建立烧结终点预测模型的步骤包括:
步骤1 求取X的潜变量,其中X为CO+CO2的浓度值
KYYTt=λt (1)
其中,K表示由所有输入数据内积构成的n×n的核矩阵,Kij=K(xi,xj),其中K为核函数;Y为烧结终点,YT为Y的转置,t为X的潜变量,λ为单位矩阵。
步骤2 估计Y的潜变量
u=YYTt (2)
其中,u为Y的潜变量,Y为烧结终点,YT为Y的转置,t为X的潜变量;
步骤3 分解核矩阵K和输出矩阵Y
K=(In-ttT)K(In-ttT) (3)
Y=(In-ttT)Y (4)
其中In为n维的单位矩阵;
步骤4 在特征空间中对核变换数据进行中心化处理
其中1n代表各元素为1的n维向量;
步骤5循环执行式(1)到式(5)各步,直到最佳数量的潜变量t被完全提取;由最佳数量潜变量构成的潜变量矩阵为T,核偏最小二乘法的回归系数矩阵B为
B=ΦTU(TTKU)-1TTY (6)
其中,ΦT为高维特征空间中的输入矩阵,T和U分别为从输入气体浓度矩阵X和输出烧结终点矩阵Y中提取得到的潜变量矩阵,K为核矩阵;
步骤6对检验集样本数据回归
其中,Φt为高维空间中的检验集样本数据矩阵,Kt为检验集样本数据构造的核矩阵,Tt为从检验集样本数据中提取出的潜变量;
其中,高斯核函数的宽度L和核偏最小二乘法提取的最佳潜变量个数t采用交叉证实法确定。该交叉证实法为现有常规方法。
采用本发明的技术方案,可以准确预测烧结终点,提高烧结生产过程的控制精度,极大地减弱烧结生产过程控制的“大滞后性”,从而实现生产过程的及时、准确调整,进而达到提高***效率和过程质量的作用。
附图说明
图1为本发明基于排放烟气分析的烧结终点预测方法流程图;
图2为本发明烧结终点预测值与实际值对比;
图3为本发明烧结终点预测的绝对误差。
具体实施方式
为充分了解本发明之目的、特征及功效,借由下述具体的实施方式,对本发明做详细说明。
烧结矿在烧结机上进心烧结的过程中排放烟气的主要成分是O2、氮氧化物(NOX)、CO、CO2、硫氧化物(SOX)以及水汽等。这些成分的变化能够反映烧结过程的变化情况及其变化特点。本发明通过实时检测烧结过程排放烟气中COx的浓度,建立烧结过程中烧结终点的预测模型,根据预测结果优化烧结工艺,提高烧结产品质量,降低烧结能耗,控制烟气排放。
如图1所示的流程图,一种基于排放烟气分析的烧结终点预测方法,该方法包括如下步骤:
(1)构建烧结烟气浓度-烧结终点数据库
实时监测、采集烧结烟气组分CO、CO2浓度和烧结终点位置,构建烧结烟气CO和CO2浓度-烧结终点数据库。
烧结烟气CO和CO2浓度的监测和采集可以采用现有常规技术完成,例如由烟气分析仪完成实时监测,由网络通讯方式实现采用光纤进行数据传输。
在获得烧结烟气CO和CO2浓度后,基于各组分浓度和与其对应的各烧结终点位置构建实时数据库,该数据库的构建可以采用Oracle10g(企业版)进行。
(2)数据预处理
在参数采集过程中,将异常数据剔除。剔除异常数据采用两种方法:1)对数据库中机速为<0.9m/min的数据直接剔除;和2)将数据不完整、缺失的数据进行直接剔除。
(3)建立烧结终点预测模型
根据烧结烟气中CO与CO2浓度值,采用核偏最小二乘法建立烧结终点预测模型,当核偏最小二乘法模型出现最小值时,对应的最佳数量的潜变量t和高斯核函数的宽度L为最终优化结果。
具体的,所述核偏最小二乘法建立烧结终点预测模型是将数据库中数据映射到高维特征空间,然后在该高维特征空间中执行偏最小二乘。
所述核偏最小二乘法建立烧结终点预测模型的步骤包括:
步骤1求取X的潜变量,其中X为CO+CO2的浓度值
KYYTt=λt (1)
其中,K表示由所有输入数据内积构成的n×n的核矩阵,Kij=K(xi,xj),其中K为核函数;Y为烧结终点,YT为Y的转置,t为X的潜变量,λ为单位矩阵;
步骤2估计Y的潜变量
u=YYTt (2)
其中,u为Y的潜变量,Y为烧结终点,YT为Y的转置,t为X的潜变量;
步骤3分解核矩阵K和输出矩阵Y
K=(In-ttT)K(In-ttT) (3)
Y=(In-ttT)Y (4)
其中In为n维的单位矩阵;
步骤4在特征空间中对核变换数据进行中心化处理
其中1n代表各元素为1的n维向量;
步骤5循环执行式(1)到式(5)各步,直到最佳数量的潜变量t被完全提取;由最佳数量潜变量构成的潜变量矩阵为T,核偏最小二乘法(KPLS)的回归系数矩阵B为
B=ΦTU(TTKU)-1TTY (6)
其中,ΦT为高维特征空间中的输入矩阵,T和U分别为从输入气体浓度矩阵X和输出烧结终点矩阵Y中提取得到的潜变量矩阵,K为核矩阵;
步骤6对检验集样本数据回归
其中,Φt为高维空间中的检验集样本数据矩阵,Kt为检验集样本数据构造的核矩阵,Tt为从检验集样本数据中提取出的潜变量;
其中,高斯核函数的宽度L和核偏最小二乘法提取的最佳潜变量个数t采用现有交叉证实法确定。
下面通过具体的实施例来阐述本发明,本领域技术人员应当理解的是,这不应被理解为对本发明权利要求范围的限制。
实施例
烧结机:450m2烧结机,总吨位14700吨;
烧结原料:自产精矿:58%,巴西粉:26%,澳洲矿粉:11%,杂料:5%;
自产精矿配比:磁铁矿:20-30%,赤铁矿:20-30%,铁酸钙:30-35%,硅酸二钙:3-5%,钙铁橄榄石:5-10%,玻璃质:3-5%。
采用西克麦哈克公司产的烟气分析仪采集烧结烟气组分浓度,从烧结控制***中采集其对应的烧结终点值,采用Oracle10g(企业版)进行构建数据库。按照如下方法剔除异常数据:1)对数据库中机速为<0.9m/min的数据直接剔除;和2)将数据不完整、缺失的数据进行直接剔除。
采用核偏最小二乘法建立烧结终点预测模型,
步骤1 求取X的潜变量
KYYTt=λt (1)
其中,K表示由所有输入数据内积构成的n×n的核矩阵,Kij=K(xi,xj),其中,核函数K为高斯核函数;
步骤2 估计Y的潜变量
u=YYTt (2)
步骤3 分解核矩阵K和输出矩阵Y
K=(In-ttT)K(In-ttT) (3)
Y=(In-ttT)Y (4)
其中In为n维的单位矩阵;
步骤4在特征空间中对核变换数据进行中心化处理
其中1n代表各元素为1的n维向量;
步骤5 循环执行式(1)到式(5)各步,直到最佳数量的潜变量t被完全提取;由最佳数量潜变量构成的潜变量矩阵为T,核偏最小二乘法的回归系数矩阵B为
B=ΦTU(TTKU)-1TTY (6)
其中,ΦT为高维特征空间中的输入矩阵,T和U分别为从输入气体浓度矩阵X和输出烧结终点矩阵Y中提取得到的潜变量矩阵,K为核矩阵;
步骤6 对检验集样本数据回归
其中,Φt为高维空间中的检验集样本数据矩阵,Kt为检验集样本数据构造的核矩阵,Tt为从检验集样本数据中提取出的潜变量;
其中,高斯核函数的宽度L和核偏最小二乘法提取的最佳潜变量个数t采用分段交叉证实法确定。
图2所示为烧结终点预测与实际值的对比图,其绝对误差如图3所示。根据本发明所述方法预测所得烧结终点与实际烧结终点的平均绝对误差小于0.13m。
Claims (4)
1.基于排放烟气分析的烧结终点预测方法,该方法采用在线烟气分析仪实时分析烧结过程中CO和CO2的浓度,同时监测烧结终点的位置数据,并将上述监测值传送给计算机,烧结终点预测模型安装在计算机中,根据在线烟气分析仪发送的CO和CO2的浓度实现烧结终点的预测,其特征在于:基于排放烟气分析的烧结终点预测方法包括如下步骤:
(1)构建烧结烟气浓度-烧结终点数据库
实时监测、采集烧结烟气CO和CO2的浓度和烧结终点位置,构建烧结烟气CO和CO2浓度-烧结终点数据库;
(2)数据预处理
对数据库中的数据进行异常数据的剔除;
(3)建立烧结终点预测模型
根据烧结烟气中CO与CO2浓度值,采用核偏最小二乘法建立烧结终点预测模型。
2.根据权利要求1所述的基于排放烟气分析的烧结终点预测方法,其特征在于:当烧结过程完成后,将监测、采集的排放烟气CO和CO2的浓度和烧结终点位置纳入数据库中。
3.根据权利要求1或2所述的基于排放烟气分析的烧结终点预测方法,其特征在于:剔除异常数据采用两种方法:1)对数据库中机速为<0.9m/min的数据直接剔除;和2)将数据不完整、缺失的数据进行直接剔除。
4.根据权利要求1所述的基于排放烟气分析的烧结终点预测方法,其特征在于:所述核偏最小二乘法建立烧结终点预测模型是将数据库中数据映射到高维特征空间,然后在该高维特征空间中执行偏最小二乘。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310381670.8A CN103439889B (zh) | 2013-08-28 | 2013-08-28 | 基于排放烟气分析的烧结终点预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310381670.8A CN103439889B (zh) | 2013-08-28 | 2013-08-28 | 基于排放烟气分析的烧结终点预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103439889A CN103439889A (zh) | 2013-12-11 |
CN103439889B true CN103439889B (zh) | 2015-10-21 |
Family
ID=49693584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310381670.8A Active CN103439889B (zh) | 2013-08-28 | 2013-08-28 | 基于排放烟气分析的烧结终点预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103439889B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103697699B (zh) * | 2013-12-26 | 2015-03-11 | 中冶长天国际工程有限责任公司 | 烧结终点控制方法及*** |
CN104833234B (zh) * | 2014-02-10 | 2016-08-24 | 山西太钢不锈钢股份有限公司 | 一种烧结机烧结过程物料仓位动态平衡控制方法 |
CN105628866A (zh) * | 2014-12-01 | 2016-06-01 | 鞍钢股份有限公司 | 一种烧结均质度检测装置及方法 |
CN108267008B (zh) * | 2016-12-30 | 2019-12-27 | 中冶长天国际工程有限责任公司 | 一种氮氧化物排放控制方法及装置 |
CN106834662B (zh) * | 2017-01-19 | 2018-09-18 | 中国地质大学(武汉) | 一种基于多工况烧结过程的co/co2比值预测方法 |
CN109211429B (zh) * | 2017-07-03 | 2020-12-22 | 佛山市顺德区美的电热电器制造有限公司 | 温度采样方法、装置及烹饪电器 |
CN112380779B (zh) * | 2020-11-27 | 2023-11-24 | 东北大学 | 一种烧结终点的鲁棒软测量方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6445963B1 (en) * | 1999-10-04 | 2002-09-03 | Fisher Rosemount Systems, Inc. | Integrated advanced control blocks in process control systems |
CN102540889A (zh) * | 2011-12-30 | 2012-07-04 | 中冶长天国际工程有限责任公司 | 一种烧结终点预测方法及*** |
CN103017530A (zh) * | 2012-12-27 | 2013-04-03 | 中冶长天国际工程有限责任公司 | 烧结终点预测方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6915172B2 (en) * | 2001-11-21 | 2005-07-05 | General Electric | Method, system and storage medium for enhancing process control |
-
2013
- 2013-08-28 CN CN201310381670.8A patent/CN103439889B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6445963B1 (en) * | 1999-10-04 | 2002-09-03 | Fisher Rosemount Systems, Inc. | Integrated advanced control blocks in process control systems |
CN102540889A (zh) * | 2011-12-30 | 2012-07-04 | 中冶长天国际工程有限责任公司 | 一种烧结终点预测方法及*** |
CN103017530A (zh) * | 2012-12-27 | 2013-04-03 | 中冶长天国际工程有限责任公司 | 烧结终点预测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN103439889A (zh) | 2013-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103439889B (zh) | 基于排放烟气分析的烧结终点预测方法 | |
Hu et al. | A multilevel prediction model of carbon efficiency based on the differential evolution algorithm for the iron ore sintering process | |
Li et al. | Dynamic time features expanding and extracting method for prediction model of sintering process quality index | |
CN104777830A (zh) | 一种基于kpca混合模型的多工况过程监控方法 | |
CN105004542A (zh) | 基于主元分析的卷烟制丝掺配加香段在线监测和故障诊断方法 | |
Yang et al. | Response surface optimization of process parameters for reduction roasting of low-grade pyrolusite by bagasse | |
CN104134120A (zh) | 一种选矿生产指标监控***及方法 | |
Li et al. | Improvements in the decision making for Cleaner Production by data mining: Case study of vanadium extraction industry using weak acid leaching process | |
CN110070217A (zh) | 一种基于过程参数的烧结矿质量预测方法 | |
CN104808648A (zh) | 一种基于k近邻的批次过程在线实时监测方法 | |
CN104793606A (zh) | 基于改进的kpca和隐马尔科夫模型的工业故障诊断方法 | |
CN110782067B (zh) | 一种基于模糊信息粒化的烧结终点波动范围预测方法 | |
CN104199389B (zh) | 金湿法冶金置换过程的优化方法 | |
CN106777684B (zh) | 一种建立综合焦比预测模型和预测综合焦比的方法 | |
CN110533082B (zh) | 一种基于双模型协同预测的烧结混合加水控制方法 | |
CN105242660A (zh) | 基于相对变化分析的多模态卷烟制叶丝过程在线监测与故障诊断方法 | |
CN104846192A (zh) | 一种计算铁矿石烧结最佳配加比例的方法 | |
CN109885012A (zh) | 一种金湿法冶金全流程实时优化补偿方法 | |
CN101423348A (zh) | 水泥回转窑烧成工况综合识别方法 | |
CN104536396A (zh) | 一种水泥生料在分解炉内分解过程的软测量建模方法 | |
CN104913639A (zh) | 基于数据融合的烧结终点控制***及控制方法 | |
CN111612210A (zh) | 一种磨煤机煤粉细度的在线优化方法 | |
CN107942873A (zh) | 一种家具制造生产线的运营成本智能化核算及监控方法 | |
CN104263960B (zh) | 金湿法冶金置换过程中置换率和金泥品位的在线预测方法 | |
Assis et al. | Artificial neural network-based committee machine for predicting fuel rate and sulfur contents of a coke 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 | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 030003, Shanxi, Taiyuan, pointed lawn area, lawn street, No. 2 Patentee after: Shanxi Taigang Stainless Steel Co.,Ltd. Address before: 030003 No. 2, sharp lawn, Taiyuan City, Shanxi Province Patentee before: Shanxi Taigang Stainless Steel Co.,Ltd. |