CN110619931A - 一种基于多时间尺度优化的烧结过程碳效优化方法 - Google Patents
一种基于多时间尺度优化的烧结过程碳效优化方法 Download PDFInfo
- Publication number
- CN110619931A CN110619931A CN201910888222.4A CN201910888222A CN110619931A CN 110619931 A CN110619931 A CN 110619931A CN 201910888222 A CN201910888222 A CN 201910888222A CN 110619931 A CN110619931 A CN 110619931A
- Authority
- CN
- China
- Prior art keywords
- optimization
- carbon efficiency
- parameters
- time scale
- sintering process
- 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.)
- Granted
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 187
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 151
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 122
- 238000005245 sintering Methods 0.000 title claims abstract description 76
- 230000008569 process Effects 0.000 title claims abstract description 72
- 239000002994 raw material Substances 0.000 claims abstract description 48
- 239000000571 coke Substances 0.000 claims abstract description 36
- 239000000463 material Substances 0.000 claims abstract description 13
- 238000004519 manufacturing process Methods 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000007246 mechanism Effects 0.000 claims abstract description 7
- 238000012352 Spearman correlation analysis Methods 0.000 claims abstract description 5
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 34
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 18
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 16
- 230000000694 effects Effects 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 230000000630 rising effect Effects 0.000 claims description 13
- 238000013528 artificial neural network Methods 0.000 claims description 10
- 229910052681 coesite Inorganic materials 0.000 claims description 10
- 229910052906 cristobalite Inorganic materials 0.000 claims description 10
- 230000002068 genetic effect Effects 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 10
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 10
- 229910052682 stishovite Inorganic materials 0.000 claims description 10
- 229910052905 tridymite Inorganic materials 0.000 claims description 10
- 229910052742 iron Inorganic materials 0.000 claims description 9
- 238000002485 combustion reaction Methods 0.000 claims description 8
- 239000000377 silicon dioxide Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 4
- 239000001569 carbon dioxide Substances 0.000 claims description 4
- 229910002091 carbon monoxide Inorganic materials 0.000 claims description 4
- ZOKXUAHZSKEQSS-UHFFFAOYSA-N tribufos Chemical compound CCCCSP(=O)(SCCCC)SCCCC ZOKXUAHZSKEQSS-UHFFFAOYSA-N 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000000843 powder Substances 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 abstract description 6
- 230000009467 reduction Effects 0.000 abstract description 6
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 238000004134 energy conservation Methods 0.000 abstract description 4
- 238000013178 mathematical model Methods 0.000 abstract 2
- 229910000831 Steel Inorganic materials 0.000 description 9
- 239000010959 steel Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 4
- 238000003723 Smelting Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
- C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
- C22B1/00—Preliminary treatment of ores or scrap
- C22B1/14—Agglomerating; Briquetting; Binding; Granulating
- C22B1/16—Sintering; Agglomerating
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Crystallography & Structural Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Manufacturing & Machinery (AREA)
- Health & Medical Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Analytical Chemistry (AREA)
- Materials Engineering (AREA)
- Mechanical Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Manufacture And Refinement Of Metals (AREA)
Abstract
本发明提供了一种基于多时间尺度优化的烧结过程碳效优化方法,包括:首先确定CO/CO2和综合焦比作为衡量碳效的指标,通过机理分析和Spearman关联性分析方法对影响碳效的状态参数进行分析,确定碳效指标预测的状态参数,并建立碳效指标与关键过程参数的数学模型;然后基于上述数学模型,考虑原料参数与操作参数的不同时间尺度问题,建立不同时间尺度上的优化策略;最后针对不同时间尺度优化的需求,采取不同的多目标优化方法,提出不同时间尺度优化策略的转换策略,以协调两种策略的实施。采用本发明所述的多时间尺度优化策略,能够在不同时间尺度下合理优化原料参数和操作参数,达到最优的碳效,为烧结过程节能降耗和绿色制造提供行之有效的解决办法。
Description
技术领域
本发明涉及钢铁烧结过程生产节能降耗领域,尤其涉及一种基于多时间尺度优化的烧结过程碳效优化方法。
背景技术
虽然我国钢铁行业在节能环保方面已经取得了非常大的进步,但与国际先进企业还存在一定的差距。由此可见,我国钢铁工业节能减排仍然存在巨大的改进空间,实现我国钢铁工业的绿色制造任重而道远。烧结是人工富矿的主要途径,产生的烧结矿用于高炉冶炼,烧结矿质量及产量对高炉生产有直接影响,以至于烧结生产对整个炼铁生产至关重要。钢铁冶炼中烧结工序需要巨大的能耗,通常都能够占到总能耗的10%-20%,仅次于高炉炼铁工序。由此可见,实现钢铁生产过程碳效在线优化、有效利用能源和减少排放的关键在于揭示能耗、排放、质量和产量与过程状态参数和操作参数之间的关系,从而运用智能优化算法优化与碳效相关的过程参数,对于降低我国吨钢综合能耗,提高钢铁工业能源利用率,节约成本有重大意义。
发明内容
为了解决上述问题,本发明提供了一种基于多时间尺度优化的烧结过程碳效优化方法,主要包括以下步骤:
S101:设定综合焦比和CO/CO2作为衡量烧结过程的碳效指标,并通过烧结过程的机理分析和数据驱动的分析方法得到影响所述碳效指标的关键过程参数,进而根据所述关键过程参数,采用RBF神经网络建立碳效指标预测模型;所述关键过程参数包括原料参数、操作参数和状态参数;
S102:结合所述碳效指标预测模型,采用一种基于多时间尺度的碳效优化方法对烧结过程进行碳效优化,即在不同的时间尺度上,采用不同的优化策略对烧结过程进行碳效优化,得到最优碳效指标、最佳操作参数和最佳原料参数,进而根据所述最佳操作参数和最佳原料参数对烧结过程进行控制,以达到实际碳效优化的目的。
进一步地,步骤S101中,具体包括如下步骤:
S201:定义CO/CO2和综合焦比为烧结过程的碳效指标;其中,CO/CO2为一氧化碳与二氧化碳的比值,其大小与烧结过程中碳的燃烧程度密切相关;综合焦比的含义是生产一吨成品烧结矿所消耗的碳量,其计算公式如下:
上式中,η为综合焦比,PC表示混合料中焦粉的配比,γb表示烧损率,Y表示成品率,QL、QS和QR分别表示大成矿、小成矿产量以及返矿,QG表示铺底料;
S202:选取碱度R、转鼓强度J和铁品位TFe作为衡量烧结矿的质量指标,并将该质量指标作为碳效优化的约束条件;
S203:通过烧结过程的机理分析,确定影响状态参数和碳效指标的原料参数和操作参数,并通过Spearman相关性分析方法计算状态参数与碳效指标之间的相关性,进而得到影响碳效指标中的综合焦比和CO/CO2的具体状态参数;
S204:在确定了状态参数与碳效指标之间的相关性的基础上,采用RBF神经网络预测模型作为数据驱动模型,以原料参数和操作参数为输入,以状态参数中的BRP和TBRP为输出,建立状态参数预测模型,以得到原料参数、操作参数与BRP和TBRP之间的对应关系,并在该对应关系的基础上,采用RBF神经网络预测模型,以状态参数为输入,碳效指标为输出,建立碳效指标预测模型。
进一步地,步骤S101中,所述原料参数包括:返矿PRS,焦粉配比PC、SiO2含量CSiO2、CaO含量CCaO、TFe含量TFe和MgO含量CMgO;所述操作参数包括:台车速度VT、料层厚度H、水分配比M和点火温度T;所述状态参数包括:风箱负压CN、垂直燃烧速度Vavg、温度上升点位置BRP、温度上升点温度TBRP、烧结终点位置BTP和烧结终点温度TBTP。
进一步地,步骤S102中,所述不同的时间尺度包括短时间尺度和长时间尺度;所述的不同的优化策略包括短时间尺度优化策略和长时间尺度优化策略;所述短时间尺度为x小时,所述长时间尺度为y天,且所述x和y均大于0,根据实际情况自行预先设定,所述操作参数每间隔一个所述短时间尺度更新一次,所述原料参数每间隔一个所述长时间尺度更新一次;每隔一个所述短时间尺度,采用所述短时间尺度优化策略对烧结过程进行一次碳效优化,每隔一个所述长时间尺度,采用所述长时间尺度优化策略对烧结过程进行依次碳效优化;根据不同时间尺度的周期不同,交替进行不同时间尺度的优化策略,达到多时间尺度优化烧结过程碳效的目的。
进一步地,所述短时间尺度优化策略,具体为:在短时间尺度上,根据所述碳效指标预测模型,以操作参数为决策变量,质量指标为约束条件,进行短时间尺度的优化;具体步骤如下:
S301:以综合焦比和CO/CO2最小为优化目标,建立如下式所示的优化目标:
min{Fi 1}
上式中,Fi 1表示第i个目标的优化函数;i=1,2;i=1时,代表CCR;i=2时,代表CO/CO2;
S302:采用多目标转单目标的线性加权和的方法,将步骤S301中的所述优化目标进行转化,并结合所述碳效指标预测模型,得到如下的优化策略:
上式中,BRP=YBRP(·)和分别表示状态参数预测模型中原料参数、操作参数与温度上升点位置BRP和温度上升点温度TBRP之间的对应关系;Fi 1=Yi(·)表示碳效指标预测模型中状态参数与碳效指标之间的对应关系;Rmin、Rmax分别为预设的碱度最小值和最大值,Jmin和Jmax分别为预设的转鼓强度最小值和最大值,TFemin和TFemax分别为预设的铁品位最小值和最大值,Opj min和Opj max分别为预设的第j个操作参数的操作范围;j=1,2,3,4;j=1时,代表H,j=2时,代表VT;j=3时,代表M;j=4时,代表T;
S303:采用遗传优化方法对步骤S302中的所述优化策略进行求解,得到碳效指标最优时的最佳操作参数。
进一步地,步骤S303中,采用遗传优化方法对步骤S302中的所述优化策略进行求解时,以优化目标线性加权和的值为优化目标,进行遗传优化算法的迭代更新,当迭代次数到了预设的数目或者输出趋于稳定时,结束迭代,获得最佳操作参数和最优碳效指标。
进一步地,所述长时间尺度优化策略,具体为:在长时间尺度上,根据所述碳效指标预测模型,以原料参数和操作参数为决策变量,质量指标为约束条件,进行长时间尺度优化具体步骤如下:
S401:以状态参数作为决策变量,以综合焦比和CO/CO2最小为优化目标,质量指标作为约束条件,并结合所述碳效指标预测模型,得到以下第一多目标优化表达式:
上式中,F1 2和分别表示综合焦比和CO/CO2;Sn min和Sn max分别表示预设的第n个状态参数的最小值和最大值,n=1,2,…,6,且n=1,2,…,6分别代表BTP、TBTP、BRP、TBRP、VAVG、CN;求解该第一多目标优化表达式,得到最佳状态参数;
S402:将α1和α2最小为优化目标,以原料参数和操作参数为决策变量,得到如下第二多目标优化表达式:
min{α1,α2}
上式中,α1和α2为中间变量;OBRP和OTBRP表示BRP和TBRP在步骤S401中所述最佳状态参数中所对应的优化值,Rmk min和Rmk max表示预设的第k个原料参数的最小值和最大值,k=1,2,…,5,且k=1,2,…,5分别代表PC、PRS、CCaO、CMgO、;Opj min和Opj max分别表示预设的第j个操作参数;j=1,2,3,4;j=1时,代表H,j=2时,代表VT;j=3时,代表M;j=4时,代表T;
求解该第二多目标优化表达式,得到碳效指标最优时的最佳操作参数和最佳原料参数。
进一步地,步骤S401中,采用多目标粒子群算法求解所述第一多目标优化表达式:
以综合焦比和CO/CO2最小为优化目标,进行多目标粒子群算法的迭代更新,当到达迭代设定次数或者输出趋于稳定时停止迭代更新,得到Pareto最优解集,采用如下公式进行解集筛选:
S=|b1i-b2i|
上式中,b1i和b2i分别表示优化的解和当前的状态参数的值,i=1,2,…,6;当S最小时的优化解即为筛选后的最佳状态参数。
进一步地,步骤S402中,采用多目标粒子群算法求解所述第二多目标优化表达式:
以α1和α2最小为优化目标,进行多目标粒子群算法的迭代更新,当到达迭代预设次数或者输出趋于稳定时停止更新,得到Pareto最优解集。
本发明提供的技术方案带来的有益效果是:本发明所提出的技术方案利用影响烧结碳效的原料参数与操作参数作为决策变量,考虑原料参数与操作参数的不同变化周期,在不同时间尺度上实现碳效优化,该优化技术能够为烧结生产降低能耗,提高碳的利用率提供有效解决方案,进而为实现钢铁工业绿色制造及节能降耗奠定基础。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例中一种基于多时间尺度优化的烧结过程碳效优化方法的流程图;
图2是本发明实施例中碳效指标预测模型原理框图;
图3是本发明实施例中碳效指标预测结果与实际值的对比图;
图4是本发明实施例中短时间尺度优化效果图;
图5是本发明实施例中长时间尺度优化效果图;
图6是本发明实施例中不同时间尺度优化策略转换示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
本发明的实施例提供了一种基于多时间尺度优化的烧结过程碳效优化方法。
请参考图1,图1是本发明实施例中一种基于多时间尺度优化的烧结过程碳效优化方法的流程图,具体包括如下步骤:
S101:设定综合焦比和CO/CO2作为衡量烧结过程的碳效指标,并通过烧结过程的机理分析和数据驱动的分析方法得到影响所述碳效指标的关键过程参数,进而根据所述关键过程参数,采用RBF神经网络建立碳效指标预测模型;所述关键过程参数包括原料参数、操作参数和状态参数;
S102:结合所述碳效指标预测模型,采用一种基于多时间尺度的碳效优化方法对烧结过程进行碳效优化,即在不同的时间尺度上,采用不同的优化策略对烧结过程进行碳效优化,得到最优碳效指标、最佳操作参数和最佳原料参数,进而根据所述最佳操作参数和最佳原料参数对烧结过程进行控制,以达到实际碳效优化的目的。
步骤S101中,具体包括如下步骤:
S201:定义CO/CO2和综合焦比为烧结过程的碳效指标;其中,CO/CO2为一氧化碳与二氧化碳的比值,其大小与烧结过程中碳的燃烧程度密切相关;综合焦比的含义是生产一吨成品烧结矿所消耗的碳量,其计算公式如下:
上式中,η为综合焦比(kg/t),PC表示混合料中焦粉的配比(%),γb表示烧损率(%),Y表示成品率(%),QL、QS和QR分别表示大成矿(kg/h)、小成矿产量(kg/h)以及返矿(kg/h),QG表示铺底料(kg/h);
S202:选取碱度R、转鼓强度J和铁品位TFe作为衡量烧结矿的质量指标,并将该质量指标作为碳效优化的约束条件;
S203:通过烧结过程的机理分析,确定影响状态参数和碳效指标的原料参数和操作参数,并通过Spearman相关性分析方法计算状态参数与碳效指标之间的相关性,进而得到影响碳效指标中的综合焦比和CO/CO2的具体状态参数;
S204:在确定了状态参数与碳效指标之间的相关性的基础上,采用RBF神经网络预测模型作为数据驱动模型,以原料参数和操作参数为输入,以状态参数中的BRP和TBRP为输出,建立状态参数预测模型,以得到原料参数、操作参数与BRP和TBRP之间的对应关系,并在该对应关系的基础上,采用RBF神经网络预测模型,以状态参数为输入,碳效指标为输出,建立碳效指标预测模型。
步骤S101中,所述原料参数包括:返矿PRS,焦粉配比PC、SiO2含量CSiO2、CaO含量CCaO、TFe含量TFe和MgO含量CMgO;所述操作参数包括:台车速度VT、料层厚度H、水分配比M和点火温度T;所述状态参数包括:风箱负压CN、垂直燃烧速度Vavg、温度上升点位置BRP、温度上升点温度TBRP、烧结终点位置BTP和烧结终点温度TBTP;其中,所述原料参数和所述操作参数均从烧结过程的历史数据中获取,所述状态参数中除温度上升点位置BRP和温度上升点温度TBRP以外的参数均从烧结过程的历史数据中获取。
步骤S102中,所述不同的时间尺度包括短时间尺度和长时间尺度;所述的不同的优化策略包括短时间尺度优化策略和长时间尺度优化策略;所述短时间尺度为x小时,所述长时间尺度为y天,且所述x和y均大于0,根据实际情况自行预先设定,所述操作参数每间隔一个所述短时间尺度更新一次,所述原料参数每间隔一个所述长时间尺度更新一次;每隔一个所述短时间尺度,采用所述短时间尺度优化策略对烧结过程进行一次碳效优化,每隔一个所述长时间尺度,采用所述长时间尺度优化策略对烧结过程进行依次碳效优化;根据不同时间尺度的周期不同,交替进行不同时间尺度的优化策略,达到在不同时间尺度上利用不同的参数,达到多时间尺度优化烧结过程碳效的目的;在本发明实施例中,x的值为2,y的值为15。
所述短时间尺度优化策略,具体为:
在短时间尺度上,根据所述碳效指标预测模型,以操作参数为决策变量,质量指标为约束条件,进行短时间尺度的优化;
所述长时间尺度优化策略,具体为:
在长时间尺度上,根据所述碳效指标预测模型,以原料参数和操作参数为决策变量,质量指标为约束条件,进行长时间尺度优化;且在长时间尺度优化上,考虑解的多样性,使其最佳参数可以根据不同的实际需求条件进行筛选,因而采用多目标进化的方法求解,得到Pareto解集,进而根据实际生产需要再所述Pareto解集中获取最佳原料参数和最佳操作参数;不同时间尺度的优化策略根据原料参数和操作参数的更新周期不同,交替进行,达到在不同时间尺度上利用不同的参数,进行优化碳效的目的。
在短时间尺度上,根据所述碳效指标预测模型,以操作参数为决策变量,质量指标为约束条件,进行短时间尺度的优化;具体步骤如下:
S301:以综合焦比和CO/CO2最小为优化目标,建立如下式所示的优化目标:
min{Fi 1}
上式中,Fi 1表示第i个目标的优化函数;i=1,2;i=1时,代表CCR;i=2时,代表CO/CO2;
S302:采用多目标转单目标的线性加权和的方法,将步骤S301中的所述优化目标进行转化,并结合所述碳效指标预测模型,得到如下的优化策略:
上式中,BRP=YBRP(·)和分别表示状态参数预测模型中原料参数、操作参数与温度上升点位置BRP和温度上升点温度TBRP之间的对应关系;Fi 1=Yi(·)表示碳效指标预测模型中状态参数与碳效指标之间的对应关系;Rmin、Rmax分别为预设的碱度最小值和最大值,Jmin和Jmax分别为预设的转鼓强度最小值和最大值,TFemin和TFemax分别为预设的铁品位最小值和最大值,Opj min和Opj max分别为预设的第j个操作参数的操作范围j=1,2,3,4;j=1时,代表H,j=2时,代表VT;j=3时,代表M;j=4时,代表T;
S303:采用遗传优化方法对步骤S302中的所述优化策略进行求解,得到碳效指标最优时的最佳操作参数。
步骤S303中,采用遗传优化方法对步骤S302中的所述优化策略进行求解时,以优化目标线性加权和的值为优化目标,进行遗传优化算法的迭代更新,当迭代次数到了预设的数目或者输出趋于稳定时,结束迭代,获得最佳操作参数和最优碳效指标。
在长时间尺度上,根据所述碳效指标预测模型,以原料参数和操作参数为决策变量,质量指标为约束条件,进行长时间尺度优化具体步骤如下:
S401:以状态参数作为决策变量,以综合焦比和CO/CO2最小为优化目标,质量指标作为约束条件,并结合所述碳效指标预测模型,得到以下第一多目标优化表达式:
上式中,F1 2和分别表示综合焦比和CO/CO2;Sn min和Sn max分别表示预设的第n个状态参数的最小值和最大值,n=1,2,…,6,且n=1,2,…,6分别代表BTP、TBTP、BRP、TBRP、VAVG、CN;求解该第一多目标优化表达式,得到最佳状态参数;
S402:将α1和α2最小为优化目标,以原料参数和操作参数为决策变量,得到如下第二多目标优化表达式:
min{α1,α2}
上式中,α1和α2为中间变量;OBRP和OTBRP表示BRP和TBRP在步骤S401中所述最佳状态参数中所对应的优化值,Rmk min和Rmk max表示预设的第k个原料参数的最小值和最大值,k=1,2,…,5,且k=1,2,…,5分别代表PC、PRS、CCaO、CMgO、;Opj min和Opj max分别表示预设的第j个操作参数;j=1,2,3,4;j=1时,代表H,j=2时,代表VT;j=3时,代表M;j=4时,代表T;
求解该第二多目标优化表达式,得到碳效指标最优时的最佳操作参数和最佳原料参数。
步骤S401中,采用多目标粒子群算法(MOPSO)求解所述第一多目标优化表达式:
以综合焦比和CO/CO2最小为优化目标,进行MOPSO算法的迭代更新,当到达迭代设定次数或者输出趋于稳定时(相邻两次输出的差值的绝对值小于预设的阈值)停止迭代更新,得到Pareto最优解集,采用如下公式进行解集筛选:
S=|b1i-b2i|
上式中,b1i和b2i分别表示优化的解和当前的状态参数的值,,i=1,2,…,6;当S最小时的优化解即为筛选后的最佳状态参数。
步骤S402中,采用多目标粒子群算法(MOPSO)求解所述第二多目标优化表达式:
以α1和α2最小为优化目标,进行MOPSO算法的迭代更新,当到达迭代预设次数或者输出趋于稳定(相邻两次输出的差值的绝对值小于预设的阈值)时停止更新,得到Pareto最优解集。
在本发明实施例中,进行了如下仿真实验,共分为以下六个步骤:
(1)确定烧结碳效指标
CO/CO2被定义为一氧化碳与二氧化碳的比值,其大小与碳的燃烧程度密切相关;综合焦比表示生产吨成品烧结矿所消耗的碳量。本发明采用综合焦比和CO/CO来作为衡量烧结过程碳效的指标;考虑烧结矿的质量,选取碱度(R)、转鼓强度(J)和铁品位(TFe)作为衡量烧结矿的质量指标。
(2)确定影响综合焦比的关键烧结参数
烧结过程具有多种物理化学变化,原料参数与操作参数影响烧结过程的状态,状态参数是反映烧结燃烧情况的一类参数,而状态参数能够反映碳效的情况。采用Spearman相关性分析方法,分析状态参数与碳效指标的关系,得到碳效指标预测模型的输入参数,如表1所示:
表1碳效指标预测模型输入参数
(3)烧结过程碳效指标预测模型
在确定影响碳效指标的输入参数以后,将烧结实际生产现场收集到的数据进行时序配准以及预处理工作:设置采样周期为45min,经过处理后的数据总共有450组,其中,350组数据用来对模型进行训练,剩下的100组数据中,50组数据采集于同一个配料单中,另外50组数据采集于多个配料单。这些数据归一化到[0,1]的范围内;
以涵盖多个配料单的数据,共350组数据作为训练数据,以原料参数与操作参数作为输入,状态参数(BRP、TBRP)作为输出,采用RBF神经网络的建模方法,建立状态参数预测模型;再以状态参数(CN、Vavg、BTP、TBTP、BRP、TBRP)作为碳效指标预测模型的输入参数,以碳效指标(CCR、CO/CO2、J、R、TFe)作为输出,碳效指标预测模型原理框图见附图2。
所建碳效指标预测模型的预测效果如图3所示,由图3经过计算得出,碳效指标综合焦比、CO/CO2以及碱度的预测范围在[-4%,4%],铁品位和转鼓强度预测精度在[-3%,3%]。因此,本模型能够较好地预测碳效指标,能够满足实际生产的需求,为后续的碳效多时间尺度优化提供了模型基础。
(4)短时间尺度优化策略
基于上述碳效指标预测模型,在短时间尺度上,以CO/CO2和综合焦比最小为决优化目标,以操作参数为决策变量,质量参数为约束条件,采用线性加权和的多目标转单目标方法,采用遗传算法实现上述优化问题,获得最佳的综合焦比和CO/CO2,并得到最佳操作参数。短时间尺度优化的结果如图4所示。前六组优化的操作参数如表2所示。
表2短时间尺度优化结果
由图4经过计算后可知,平均综合焦比实际值为52.95kg/t,优化后的综合焦比平均值为48.91kg,综合焦比平均降低了4.04kg/t。CO/CO2的实际值平均为0.15,平均优化值为0.132,平均降低了0.018。这表明了短时间尺度优化策略的有效性。同时,最佳操作参数用来指导操作人员进行基础设定。
(5)长时间尺度优化策略
基于碳效指标预测模型,在长时间尺度上,以CO/CO2和综合焦比最小为决优化目标,以操作参数和原料参数为决策变量,质量参数为约束条件,采用NSGA-II的多目标优化方法,实现长时间尺度优化策略的优化,获得最佳的综合焦比和CO/CO2,并得到最佳的操作参数和原料配比。长时间尺度优化的结果如图5所示。第一组获得的优化参数如表3所示。
表3长时间尺度优化结果
由图5所示,用了三组数据对长时间尺度优化策略进行了验证,综合焦比优化平均值为48.34kg/t,实际值为51.13kg/t,平均降低了2.79kg/t。而CO/CO2,优化结果表明,其实际平均值为0.129,实际值为0.146,平均降低了0.17。这表明了长时间尺度优化策略能有效降低碳效。
(6)不同时间尺度优化策略的转换
由于存在两种不同的优化策略,需要制定一个协调实施策略,使不同的优化策略能够共同实施,达到最优碳效指标,转换策略如图6所示。假设tl0和tl1是长时间尺度优化策略实施的时间点,ts0和tsn是短时间尺度优化策略实施的时间点。tl0和tl1之间的时间间隔是15天,ts0和ts1之间的时间间隔是2小时。在tl0时实施一次长时间尺度优化策略,得到最佳的原料参数和操作参数,作为原料配比的设定值和基础控制回路的设定值。2小时以后,在ts0时刻,实施一次短时间优化策略来得到最佳的操作参数。在tl1时刻到来之前,每两个小时实施一次短时间尺度优化策略。在tl1进行第二次长时间尺度优化策略,以此类推协调进行两种优化策略,以达到最佳的碳效优化。
本发明的有益效果是:本发明所提出的技术方案利用影响烧结碳效的原料参数与操作参数作为决策变量,考虑原料参数与操作参数的不同变化周期,在不同时间尺度上实现碳效优化,该优化技术能够为烧结生产降低能耗,提高碳的利用率提供有效解决方案,进而为实现钢铁工业绿色制造及节能降耗奠定基础。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:包括以下步骤:
S101:设定综合焦比和CO/CO2作为衡量烧结过程的碳效指标,并通过烧结过程的机理分析和数据驱动的分析方法得到影响所述碳效指标的关键过程参数,进而根据所述关键过程参数,采用RBF神经网络建立碳效指标预测模型;所述关键过程参数包括原料参数、操作参数和状态参数;
S102:结合所述碳效指标预测模型,采用一种基于多时间尺度的碳效优化方法对烧结过程进行碳效优化,即在不同的时间尺度上,采用不同的优化策略对烧结过程进行碳效优化,得到最优碳效指标、最佳操作参数和最佳原料参数,进而根据所述最佳操作参数和最佳原料参数对烧结过程进行控制,以达到实际碳效优化的目的。
2.如权利要求1所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:步骤S101中,具体包括如下步骤:
S201:定义CO/CO2和综合焦比为烧结过程的碳效指标;其中,CO/CO2为一氧化碳与二氧化碳的比值,其大小与烧结过程中碳的燃烧程度密切相关;综合焦比的含义是生产一吨成品烧结矿所消耗的碳量,其计算公式如下:
上式中,η为综合焦比,PC表示混合料中焦粉的配比,γb表示烧损率,Y表示成品率,QL、QS和QR分别表示大成矿、小成矿产量以及返矿,QG表示铺底料;
S202:选取碱度R、转鼓强度J和铁品位TFe作为衡量烧结矿的质量指标,并将该质量指标作为碳效优化的约束条件;
S203:通过烧结过程的机理分析,确定影响状态参数和碳效指标的原料参数和操作参数,并通过Spearman相关性分析方法计算状态参数与碳效指标之间的相关性,进而得到影响碳效指标中的综合焦比和CO/CO2的具体状态参数;
S204:在确定了状态参数与碳效指标之间的相关性的基础上,采用RBF神经网络预测模型作为数据驱动模型,以原料参数和操作参数为输入,以状态参数中的BRP和TBRP为输出,建立状态参数预测模型,以得到原料参数、操作参数与BRP和TBRP之间的对应关系,并在该对应关系的基础上,采用RBF神经网络预测模型,以状态参数为输入,碳效指标为输出,建立碳效指标预测模型。
3.如权利要求1所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:步骤S101中,所述原料参数包括:返矿PRS,焦粉配比PC、SiO2含量CSiO2、CaO含量CCaO、TFe含量TFe和MgO含量CMgO;所述操作参数包括:台车速度VT、料层厚度H、水分配比M和点火温度T;所述状态参数包括:风箱负压CN、垂直燃烧速度Vavg、温度上升点位置BRP、温度上升点温度TBRP、烧结终点位置BTP和烧结终点温度TBTP。
4.如权利要求1所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:步骤S102中,所述不同的时间尺度包括短时间尺度和长时间尺度;所述的不同的优化策略包括短时间尺度优化策略和长时间尺度优化策略;所述短时间尺度为x小时,所述长时间尺度为y天,且所述x和y均大于0,根据实际情况自行预先设定,所述操作参数每间隔一个所述短时间尺度更新一次,所述原料参数每间隔一个所述长时间尺度更新一次;每隔一个所述短时间尺度,采用所述短时间尺度优化策略对烧结过程进行一次碳效优化,每隔一个所述长时间尺度,采用所述长时间尺度优化策略对烧结过程进行依次碳效优化;根据不同时间尺度的周期不同,交替进行不同时间尺度的优化策略,达到多时间尺度优化烧结过程碳效的目的。
5.如权利要求4所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:所述短时间尺度优化策略,具体为:在短时间尺度上,根据所述碳效指标预测模型,以操作参数为决策变量,质量指标为约束条件,进行短时间尺度的优化;具体步骤如下:
S301:以综合焦比和CO/CO2最小为优化目标,建立如下式所示的优化目标:
min{Fi 1}
上式中,Fi 1表示第i个目标的优化函数;i=1,2;i=1时,代表CCR;i=2时,代表CO/CO2;
S302:采用多目标转单目标的线性加权和的方法,将步骤S301中的所述优化目标进行转化,并结合所述碳效指标预测模型,得到如下的优化策略:
上式中,BRP=YBRP(·)和分别表示状态参数预测模型中原料参数、操作参数与温度上升点位置BRP和温度上升点温度TBRP之间的对应关系;Fi 1=Yi(·)表示碳效指标预测模型中状态参数与碳效指标之间的对应关系;Rmin、Rmax分别为预设的碱度最小值和最大值,Jmin和Jmax分别为预设的转鼓强度最小值和最大值,TFemin和TFemax分别为预设的铁品位最小值和最大值,Opjmin和Opjmax分别为预设的第j个操作参数的操作范围;j=1,2,3,4;j=1时,代表H,j=2时,代表VT;j=3时,代表M;j=4时,代表T;
S303:采用遗传优化方法对步骤S302中的所述优化策略进行求解,得到碳效指标最优时的最佳操作参数。
6.如权利要求5所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:步骤S303中,采用遗传优化方法对步骤S302中的所述优化策略进行求解时,以优化目标线性加权和的值为优化目标,进行遗传优化算法的迭代更新,当迭代次数到了预设的数目或者输出趋于稳定时,结束迭代,获得最佳操作参数和最优碳效指标。
7.如权利要求4所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:所述长时间尺度优化策略,具体为:在长时间尺度上,根据所述碳效指标预测模型,以原料参数和操作参数为决策变量,质量指标为约束条件,进行长时间尺度优化具体步骤如下:
S401:以状态参数作为决策变量,以综合焦比和CO/CO2最小为优化目标,质量指标作为约束条件,并结合所述碳效指标预测模型,得到以下第一多目标优化表达式:
上式中,F1 2和分别表示综合焦比和CO/CO2;Snmin和Snmax分别表示预设的第n个状态参数的最小值和最大值,n=1,2,…,6,且n=1,2,…,6分别代表BTP、TBTP、BRP、TBRP、VAVG、CN;求解该第一多目标优化表达式,得到最佳状态参数;
S402:将α1和α2最小为优化目标,以原料参数和操作参数为决策变量,得到如下第二多目标优化表达式:
min{α1,α2}
上式中,α1和α2为中间变量;OBRP和OTBRP表示BRP和TBRP在步骤S401中所述最佳状态参数中所对应的优化值,Rmkmin和Rmkmax表示预设的第k个原料参数的最小值和最大值,k=1,2,…,5,且k=1,2,…,5分别代表PC、PRS、CCaO、CMgO、CSiO2;Opjmin和Opjmax分别表示预设的第j个操作参数;j=1,2,3,4;j=1时,代表H,j=2时,代表VT;j=3时,代表M;j=4时,代表T;
求解该第二多目标优化表达式,得到碳效指标最优时的最佳操作参数和最佳原料参数。
8.如权利要求7所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:步骤S401中,采用多目标粒子群算法求解所述第一多目标优化表达式:
以综合焦比和CO/CO2最小为优化目标,进行多目标粒子群算法的迭代更新,当到达迭代设定次数或者输出趋于稳定时停止迭代更新,得到Pareto最优解集,采用如下公式进行解集筛选:
S=|b1i-b2i|
上式中,b1i和b2i分别表示优化的解和当前的状态参数的值,i=1,2,…,6;当S最小时的优化解即为筛选后的最佳状态参数。
9.如权利要求7所述的一种基于多时间尺度优化的烧结过程碳效优化方法,其特征在于:步骤S402中,采用多目标粒子群算法求解所述第二多目标优化表达式:
以α1和α2最小为优化目标,进行多目标粒子群算法的迭代更新,当到达迭代预设次数或者输出趋于稳定时停止更新,得到Pareto最优解集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910888222.4A CN110619931B (zh) | 2019-09-19 | 2019-09-19 | 一种基于多时间尺度优化的烧结过程碳效优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910888222.4A CN110619931B (zh) | 2019-09-19 | 2019-09-19 | 一种基于多时间尺度优化的烧结过程碳效优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110619931A true CN110619931A (zh) | 2019-12-27 |
CN110619931B CN110619931B (zh) | 2021-11-30 |
Family
ID=68923606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910888222.4A Expired - Fee Related CN110619931B (zh) | 2019-09-19 | 2019-09-19 | 一种基于多时间尺度优化的烧结过程碳效优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110619931B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112506057A (zh) * | 2020-12-02 | 2021-03-16 | 郑州轻工业大学 | 不确定奇异摄动***在线多时间尺度快速自适应控制方法 |
CN115200370A (zh) * | 2022-07-14 | 2022-10-18 | 中冶华天工程技术有限公司 | 一种烧结过程多特征点分段控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106834662A (zh) * | 2017-01-19 | 2017-06-13 | 中国地质大学(武汉) | 一种基于多工况烧结过程的co/co2比值预测方法 |
CN109376500A (zh) * | 2018-12-19 | 2019-02-22 | 中国地质大学(武汉) | 基于预测模型的烧结过程综合焦比在线优化方法及*** |
CN109918702A (zh) * | 2019-01-03 | 2019-06-21 | 上海交通大学 | 一种高炉配料与操作的协同多目标优化方法 |
-
2019
- 2019-09-19 CN CN201910888222.4A patent/CN110619931B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106834662A (zh) * | 2017-01-19 | 2017-06-13 | 中国地质大学(武汉) | 一种基于多工况烧结过程的co/co2比值预测方法 |
CN109376500A (zh) * | 2018-12-19 | 2019-02-22 | 中国地质大学(武汉) | 基于预测模型的烧结过程综合焦比在线优化方法及*** |
CN109918702A (zh) * | 2019-01-03 | 2019-06-21 | 上海交通大学 | 一种高炉配料与操作的协同多目标优化方法 |
Non-Patent Citations (1)
Title |
---|
王海东等: "铁矿烧结过程能耗数学模型的研究进展", 《钢铁》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112506057A (zh) * | 2020-12-02 | 2021-03-16 | 郑州轻工业大学 | 不确定奇异摄动***在线多时间尺度快速自适应控制方法 |
CN112506057B (zh) * | 2020-12-02 | 2021-09-17 | 郑州轻工业大学 | 不确定奇异摄动***在线多时间尺度快速自适应控制方法 |
CN115200370A (zh) * | 2022-07-14 | 2022-10-18 | 中冶华天工程技术有限公司 | 一种烧结过程多特征点分段控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110619931B (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109583118B (zh) | 一种烧结配比计算及烧结矿成本优化方法 | |
CN110619931B (zh) | 一种基于多时间尺度优化的烧结过程碳效优化方法 | |
CN107368125B (zh) | 一种基于cbr与rbr并行混合推理的高炉炉温控制***及方法 | |
CN109376500A (zh) | 基于预测模型的烧结过程综合焦比在线优化方法及*** | |
CN110610255B (zh) | 基于智能算法的炼铁过程配料优化方法 | |
CN102031319B (zh) | 一种高炉铁水含硅量的预报方法 | |
CN109447346A (zh) | 基于灰色预测与神经网络组合模型转炉氧耗量预测方法 | |
CN103194553A (zh) | 一种基于最小二乘支持向量机的钢铁冶炼高炉氧气使用量控制方法 | |
CN103544273A (zh) | 利用模式识别技术对炉况整体状态进行评估的方法 | |
CN106777684B (zh) | 一种建立综合焦比预测模型和预测综合焦比的方法 | |
CN110570019A (zh) | 一种基于宽度学习的烧结过程综合焦比时序预测方法 | |
CN111308894B (zh) | 一种烧结过程耦合配矿及控制过程的参数调节方法 | |
CN115522012B (zh) | 一种大型转炉控制tsc阶段吹氧量的方法 | |
CN111639820A (zh) | 一种水泥炉窑生产的能耗参数优化方法及*** | |
CN111047103B (zh) | 基于多工况模式辨识的烧结过程碳效在线预测方法及设备 | |
CN108062583A (zh) | 一种面向节能减排的加热炉工艺参数寻优方法 | |
CN108004368A (zh) | 智能化自动炼钢方法及装置 | |
CN111798023B (zh) | 一种炼钢烧结生产中综合焦比预测的方法 | |
CN115081319B (zh) | 铝电解生产过程智能决策、多目标优化方法 | |
CN110361972B (zh) | 一种基于水泥生料燃烧的分解炉炉温快速优化控制方法 | |
CN110617716B (zh) | 一种混钢加热的生产节奏及炉温在线控制方法 | |
CN1594601A (zh) | 基于平衡点的稳定高炉操作的方法 | |
CN116127345B (zh) | 基于深度聚类生成对抗网络的转炉炼钢工艺模式设计方法 | |
CN112410486B (zh) | 一种高炉炉缸护炉效果的评判方法 | |
Chen et al. | Prediction model of the sinter comprehensive performance based on neural network |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211130 |