CN102653801B - 基于遗传算法建立的无料钟高炉炉顶布料控制方法 - Google Patents
基于遗传算法建立的无料钟高炉炉顶布料控制方法 Download PDFInfo
- Publication number
- CN102653801B CN102653801B CN201210125508.5A CN201210125508A CN102653801B CN 102653801 B CN102653801 B CN 102653801B CN 201210125508 A CN201210125508 A CN 201210125508A CN 102653801 B CN102653801 B CN 102653801B
- Authority
- CN
- China
- Prior art keywords
- cloth
- population
- ore
- individuality
- coke ratio
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于遗传算法建立的无料钟高炉炉顶布料控制方法,该方法应用与于无料钟高炉炉顶布料控制。本发明公开了一种基于遗传算法建立的无料钟高炉炉顶布料控制方法,根据期望的矿焦比,利用遗传算法计算能够实现期望矿焦比的最优布料矩阵,从而实现高炉顺行,提高煤气利用率,达到节能减排的目标。其特点在于:定义了矿焦比数学表达式;给出了期望矿焦比的数学模型;结合当前炉矿信息,找到合适的期望矿焦比,利用自适应遗传算法找到最优布料矩阵,完成布料控制,实现了高炉生产的在线控制。
Description
技术领域
本发明涉及一种基于自适应遗传算法的控制方法,该方法应用与于无料钟高炉炉顶布料控制。
背景技术
节能减排是当今工业发展的方向。炼铁行业作为冶金行业中比重较大的领域,更是实现节能减排不可或缺的部分。实现炼铁领域节能减排的关键因素是保障高炉炉内煤气流的合理分布。矿焦比分布是煤气流分布的内在表现。炉顶布料又是调节矿焦比分布的重要手段。因此,通过炉顶布料实现期望的矿焦比,从而实现高炉炉内煤气流的合理分布,该方法在很大程度上可以实现节能减排的目标。无料钟高炉的炉顶布料方式调节手段丰富,更能充分发挥无钟炉顶布料的优越性,从而更好地实现期望的矿焦比和煤气流分布。
由于高炉内部高温、高压、复杂多变等恶劣条件的限制,高炉布料控制的优化问题成为了一种非线性、多峰值的复杂优化问题,因此很多经典的控制方法都不能有效的指导高炉炉顶的布料控制。目前的高炉生产中,在布料方式选择上,通常是根据炉长经验结合布料模型分析,进行布料矩阵选择的。也有部分企业是根据现场数据和炉长的经验建立了布料专家***,该方法在一定程度上对高炉的生产起到了积极作用,提高了矿焦比,提高了煤气利用率。但是这些方法具有一定的滞后性,也就是只有当现场数据反映出问题(煤气流分布恶化、矿焦比低等)后才能进行布料控制调节,炉内状况不能保证始终处于合理阶段,而且这些方法在如何选择多环布料中的布料矩阵(环数、溜槽倾角和布料圈数等)的问题上没有严谨的理论依据。
因此,有必要设计一种无料钟高炉炉顶布料控制方法,根据当前矿焦比分布实时调节布料矩阵,使高炉矿焦比、煤气流分布始终处于合理,进一步提高煤气利用率和降低焦比,最终实现节能减排的目的。
发明内容
本发明的目的是提供一种基于遗传算法的无钟高炉布料闭环控制方法,以找到最优布料布阵,以保障矿焦比分布始终处于当前最佳的状态,实现炉况长期稳定顺行、提高煤气利用率、降低焦炭比以及节能减排的目的。
所述遗传算法(genetic algorithms,GA)是在达尔文的生物进化论和孟德尔的遗传学说的基础上建立起来的一种模拟自然选择和遗传机制的寻优方法。该算法模拟了生物的遗传、进化过程,并结合了随机统计原理,采用启发式的全局搜索方法,具有很强的全局搜索和优化能力,特别适应于处理高炉布料控制这样多维多峰的非线性优化问题。因此,遗传算法成为解决无钟高炉炉顶布料控制方法的不二选择。
遗传算法包括以下几个基本步骤:初始化,选择,交叉,变异等,其中选择、交叉和变异操作都是基于个体的适应度来操作的。
本发明的核心思想是利用自适应遗传算法,根据当前所需的矿焦比分布,找到最优布料矩阵,实现期望矿焦比分布。
本发明的技术方案如下:一种基于遗传算法建立的无料钟高炉炉顶布料控制方法,特征在于:采用自适应遗传算法,结合期望的矿焦比曲线,求得所需的布料矩阵,通过多环布料方式实现布料控制,使得所得的矿焦比与期望矿焦比尽可能吻合,进而达到提高高炉煤气利用率的目的,所述方法的具体步骤如下:
(1)定义矿焦比
假设布料次序为JK,即先布焦后布矿,我们定义矿焦比为:
式中,fb(x)表示已知初始料面形状,g(x)表示实施布料操作后的焦炭料面形状,f(x)表示实施布料操作后的矿石料面形状,x表示炉喉的径向位置;g(x)和f(x)都是在布料模型的基础上,根据当前的布料矩阵计算得到的;
(2)期望矿焦比曲线数学模型
其数学模型可描述为:
式中[0,r1)表示炉心敞开区,[r2,R]表示边缘抑制区,表示大部分截面上最大程度利用煤气时的矿焦比,δ表示敞开中心煤气和抑制边缘煤气的程度;根据专家经验和数据挖掘方法,构建参数r1,r2,和δ与当前炉况信息包括温度、煤气成分、矿焦比和顶压等的关联规则,并建立小型的专家***;然后,根据当前的炉况信息包括温度、煤气成分、顶压和矿焦比等,自动设置下次应达到的期望矿焦比曲线;
(3)利用遗传算法计算布料矩阵
根据布料公式以及高炉的具体数据,可以得知布料角度α的取值范围,从而对个体的可行性进行判断;α的取值范围如下所示:
α∈[αmin,αmax]
其中,α表示布料角度,αmin表示布料角度的最小值,αmax表示布料角度的最大值;
本专利中,假定布料所需的矿石体积VO和焦炭体积VC都是确定的,则自适应遗传算法的具体流程及实施方法如下:
(3-1)算法的初始化:设定种群大小NP,最大遗传代数GM,选择概率PS,初始交叉概率PC0,最小交叉概率PCmin(0<PCmin<PC0<1),初始变异概率PM0,最小变异概率PMmin(0<PMmin<PM0<1),初始交叉系数α0,最小交叉系数αmin(α0>αmin>0),解的误差限δ等,并令当前遗传代数G=1,k=1;
(3-2)以当前生产所使用的布料矩阵M0为基础,在给定的变化范围内随机产生NP个个体,判断个体是否为可行解并对不可行解进行处理;
(3-4)确定适应度函数:根据适应度函数计算当前种群中各个个体的适应度值,其中适应度函数的表达式如下:
f=∫|O/C*(x)-O/C(x)|2dx,x∈(0,R)
(3-5)选择:以个体的适应度值为依据,选择定量的个体进行进化操作;
其中而每个个体的变异概率与个体的适应度值有关,即
其中,favg表示当代种群pk的平均适应度值,fmax表示当代种群Pk中最大的适应度值,f′表示进行交叉操作的两个个体中的较大适应度值。
在实施交叉操作时运用非均匀算术交叉,交叉系数αk是由当前遗传代数G、最大遗传代数GM以及初始交叉系数α0确定的:
αk=α0-(α0-αmin)*G/GM
判断新产生的个体是否为可行解并对不可行解进行处理;
判断新产生的个体是否为可行解并对不可行解进行处理;
(3-9)记录当前适应度值最小的个体Mk,令G=G+1,k=k+1;
(3-10)进行遗传算法操作,直到满足算法结束条件,找到最优解M。
附图说明
图1:无料钟高炉炉顶布料控制结构框图。
图2:理想矿焦比分布示意图
具体实施方式
1)定义矿焦比
定义新的矿焦比曲线数学模型,从而使得矿焦比的定义更加科学,从而更好的反映煤气流分布。假设布料次序为JK,即先布焦后布矿,我们定义矿焦比为:
式中,fb(x)表示已知初始料面形状,g(x)表示实施布料操作后的焦炭料面形状,f(x)表示实施布料操作后的矿石料面形状,x表示炉喉的径向位置。g(x)和f(x)都是在布料模型的基础上,根据当前的布料矩阵计算得到的。因此,矿焦比可以很好的反映当前布料矩阵的优劣情况。
2)期望矿焦比曲线数学模型
经过众多冶金专家大量的观察和研究,认为理想的矿焦比分布的大致形状应如图2所示。其数学模型可描述为:
式中[0,r1)表示炉心敞开区,[r2,R]表示边缘抑制区,表示大部分截面上最大程度利用煤气时的矿焦比,δ表示敞开中心煤气和抑制边缘煤气的程度。根据专家经验和数据挖掘方法,构建参数r1,r2,和δ与当前炉况信息(温度、煤气成分、矿焦比和顶压等)的关联规则,并建立小型的专家***。然后,便可根据当前的炉况信息(温度、煤气成分、顶压和矿焦比等),自动设置下次应达到的期望矿焦比曲线。
3)利用遗传算法计算布料矩阵
利用遗传算法进行迭代计算,首先要确定采用那种遗传算法。本专利中,我们采用了自适应遗传算法。这是由于高炉布料矩阵的求解问题搜索空间较大,搜索结果有可能局部收敛或者收敛速度太慢,因此,需要所用的算法能够在初始阶段有较大的交叉概率以加大搜索空间,在中期阶段有较大的变异概率以便突破局部收敛的限制,而在后期有较小的交叉概率和变异概率以加快算法收敛。自适应遗传算法的交叉概率、变异概率随着遗传代数的变化而变化的趋势能够很好得吻合我们的要求。而且我们在算法中加入了优秀基因保留策略,从而保证了算法的收敛性。另外,在遗传算法的过程中,不可避免会产生不可行解,其在无钟高炉布料时具体表现为布料角度过大时炉料在下料过程中撞击炉墙。因此,进行算法操作时,需要对不可行解进行筛选和删除工作。根据布料公式以及高炉的具体数据,可以得知布料角度α的取值范围,从而对个体的可行性进行判断。α的取值范围如下所示:
α∈[αmin,αmax]
其中,α表示布料角度,αmin表示布料角度的最小值,αmax表示布料角度的最大值。
本专利中,假定布料所需的矿石体积VO和焦炭体积VC都是确定的,则自适应遗传算法的具体流程及实施方法如下:
a)算法的初始化。编码是遗传算法的基础工作,本专利中,我们采用实值编码方式,主要是由于该方法适用于求解像高炉布料矩阵这样搜索空间大、精度要求高的应用问题。求解的变量设定为布料矩阵M(包括m个焦炭布料角度和布料圈数以及m个矿石布料角度和布料圈数),即个体的基因数为NV=4m。例如,当前生产所使用的布料矩阵M0的具体的表现形式如下:
其中,上标“0”表示还没有进行遗传算法迭代,下标“c”表示与焦炭有关的量,下标“o”表示与矿石有关的量,表示第i个焦炭布料角度,表示其对应的布料圈数,表示第i个矿石布料角度,表示其对应的布料圈数,i=1,2,…m。
设定适应度函数f为期望矿焦比O/C*与计算的矿焦比O/C之间的误差的平方和,即
f=∫|O/C*(x)-O/C(x)|2dx,x∈(0,R)。
其中,f值越小,说明计算的矿焦比与期望矿焦比越接近,解(当前布料矩阵)的质量也就越高。设定种群大小NP,最大遗传代数GM,选择概率PS,初始交叉概率PC0,最小交叉概率PCmin(0<PCmin<PC0<1),初始变异概率PM0,最小变异概率PMmin(0<PMmin<PM0<1),初始交叉系数α0,最小交叉系数αmin(α0>αmin>0),解的误差限δ,当前遗传代数G=1;
其中下标“I”表示为当前种群中的第I个个体(I=1,2,…,NP),表示第1代种群中第I个个体的第i个焦炭布料角度,表示其对应的布料圈数,表示第1代种群中第I个个体的第i个矿石布料角度,表示其对应的布料圈数,i=1,2,…m,且有:
这NP个个体组成初始种群P1。判断每个个体所包括的布料角度是否在可行解范围内,删除不可行解,用上面所示的方法产生新的个体,直至所有个体都符合要求;计算每个个体的适应度值并记录适应度值最小的个体为M1;或者可以先设定一个较小的遗传代数g,令其运行过程中最好的一代个体为整个算法的初始种群P1;
c)遗传算法算子的实现:在得到初始种群以及所需的参数后,就可以进入遗传算法算子的实现即遗传算法的迭代过程了,其具体步骤如下:
i.判断G=GM是否成立,成立则退出循环,MG就是最优布料矩阵;
否则进入下一步:
iii.交叉:对种群包含的个体(临近的个体,没有重复)进行交叉操作,本专利中,当前种群的最大交叉概率与遗传代数G有关,即
本专利中运用非均匀算术交叉法进行交叉操作,因此需要确定交叉系数αG,αG是由当前遗传代数G、最大遗传代数GM以及初始交叉系数α0确定的:
αG=α0-(α0-αmin)*G/GM
然后利用轮盘赌法确定个体是否进行交叉操作,对要实施交叉操作的两个个体,随机产生交叉位置LC(1≤LC≤4m),对两个交叉个体LC后的基因实施非均匀算数交叉操作,产生新的个体,组成子代种群具体操作如下:
判断新产生的个体是否为可行解,如果为不可行解,则取消该次交叉操作,重新进行交叉操作,直至产生的个体符合要求,进入下一步;
在进行变异操作时,假设要实施变异操作的基因为X,首先确定X的可取值范围[X1,X2],然后用符合此范围内的均匀分布的随机数XM来替换X;利用轮盘赌法确定是否对该个体进行变异操作,对要进行变异操作的个体,首先设定一个误差限δ以及变异次数TM;在变异时,首先在选中个体的布矿圈数对应的基因位置即中随机选择一个位置进行变异,生成新的种群pM1,令T1=1,计算其个体的适应度函数,记录最优个体(即适应度值最小的个体)MM以及其适应度值fM,如果fM≤δ,则变异操作结束;如果fM≥δ,则如上继续变异,并令T1=T1+1,直到fM≤δ或者T1=TM,此时如果fM≤δ,则算法结束,MM即为最优解;否则,依次在种群的基础上,对选中的个体的布矿角度布焦圈数布焦角度对应的基因位置中随机选择一个位置进行变异,得到种群pM2、pM3、pM4直到fM≤δ或者T=TM,记录最优种群PM、最优个体MM以及其适应度值fM,如果fM≤δ,则算法结束,MM即为最优解(所求解),否则进入下一步;判断新产生的个体是否为可行解,如果为不可行解,则取消该次交叉操作,重新进行交叉操作,直至产生的个体符合要求;
vi.令G=G+1,记录当前适应度值最小的个体为MG,返回到(i);
d)MG即为我们所求的布料矩阵。
Claims (1)
1.一种基于遗传算法建立的无料钟高炉炉顶布料控制方法,特征在于:采用自适应遗传算法,结合期望的矿焦比曲线,求得所需的布料矩阵,通过多环布料方式实现布料控制,使得所得的矿焦比与期望矿焦比尽可能吻合,进而达到提高高炉煤气利用率的目的,所述方法的具体步骤如下:
(1)定义矿焦比
假设布料次序为JK,即先布焦后布矿,定义矿焦比为:
(2)期望矿焦比曲线数学模型
其数学模型可描述为:
式中表示炉心敞开区,表示边缘抑制区,表示大部分截面上最大程度利用煤气时的矿焦比,表示敞开中心煤气和抑制边缘煤气的程度;根据专家经验和数据挖掘方法,构建参数和与当前炉况信息包括温度、煤气成分、矿焦比和顶压的关联规则,并建立小型的专家***;然后,根据当前的炉况信息包括温度、煤气成分、顶压和矿焦比,自动设置下次应达到的期望矿焦比曲线;
(3)利用遗传算法计算布料矩阵
(3-4)确定适应度函数:根据适应度函数计算当前种群中各个个体的适应度值,其中适应度函数的表达式如下:
(3-5)选择:以个体的适应度值为依据,选择定量的个体进行进化操作;
判断新产生的个体是否为可行解并对不可行解进行处理;
判断新产生的个体是否为可行解并对不可行解进行处理;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210125508.5A CN102653801B (zh) | 2012-04-25 | 2012-04-25 | 基于遗传算法建立的无料钟高炉炉顶布料控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210125508.5A CN102653801B (zh) | 2012-04-25 | 2012-04-25 | 基于遗传算法建立的无料钟高炉炉顶布料控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102653801A CN102653801A (zh) | 2012-09-05 |
CN102653801B true CN102653801B (zh) | 2014-01-01 |
Family
ID=46729530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210125508.5A Expired - Fee Related CN102653801B (zh) | 2012-04-25 | 2012-04-25 | 基于遗传算法建立的无料钟高炉炉顶布料控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102653801B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103173584A (zh) * | 2013-04-01 | 2013-06-26 | 鞍山市海汇自动化有限公司 | 具有自学习控制功能的高炉布料控制*** |
CN103455003B (zh) * | 2013-09-03 | 2016-01-13 | 邵光震 | 用于石化行业生产的监控方法及*** |
CN103950930B (zh) * | 2013-12-03 | 2015-09-23 | 国家电网公司 | 一种用于电石生产配料的控制方法 |
CN104133945B (zh) * | 2014-07-15 | 2017-03-08 | 东北大学 | 一种高炉布料过程径向矿焦比的控制方法 |
CN105243190B (zh) * | 2015-09-14 | 2018-08-07 | 内蒙古科技大学 | 一种高炉布料过程料面输出形状的建模方法 |
CN107119159B (zh) * | 2017-06-14 | 2019-02-22 | 内蒙古科技大学 | 一种存在整数规划问题的高炉布料过程布料矩阵的优化方法 |
CN109146233A (zh) * | 2018-06-29 | 2019-01-04 | 首钢京唐钢铁联合有限责任公司 | 一种无钟炉顶矿石自学习布料的方法及装置 |
CN109284942B (zh) * | 2018-10-31 | 2021-08-03 | 山西太钢不锈钢股份有限公司 | 用于确定高炉原燃料参数合理范围的方法 |
CN109977814A (zh) * | 2019-03-13 | 2019-07-05 | 武汉大学 | 一种基于一致化LBP的AdaBoost行人检测方法 |
CN114139438A (zh) * | 2021-10-26 | 2022-03-04 | 中冶南方工程技术有限公司 | 一种高炉炉料轨迹模型构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102181592A (zh) * | 2011-05-16 | 2011-09-14 | 北京科技大学 | 基于多点雷达数据的无料钟高炉炉顶布料闭环控制方法 |
-
2012
- 2012-04-25 CN CN201210125508.5A patent/CN102653801B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102653801A (zh) | 2012-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102653801B (zh) | 基于遗传算法建立的无料钟高炉炉顶布料控制方法 | |
Li et al. | Composite particle swarm optimizer with historical memory for function optimization | |
Gaing | Constrained dynamic economic dispatch solution using particle swarm optimization | |
CN106599519A (zh) | 面向生产订单组合优化的中厚板母板与板坯协同设计方法及*** | |
CN105786759A (zh) | 一种改进标准混合蛙跳算法的方法 | |
CN108038306A (zh) | 一种面向海量高维数据的电站锅炉燃烧建模及优化的方法 | |
CN106119458A (zh) | 基于bp神经网络的转炉炼钢工艺成本控制方法及*** | |
CN110400009A (zh) | 基于自适应遗传算法的高炉炼铁多目标智能优化方法 | |
CN106019935B (zh) | 基于约束模糊关联规则的多目标锅炉燃烧优化 | |
Wang et al. | A hybrid discrete imperialist competition algorithm for fuzzy job-shop scheduling problems | |
Xu et al. | An efficient hybrid multi-objective particle swarm optimization with a multi-objective dichotomy line search | |
CN103279793A (zh) | 一种确定环境下的无人飞行器编队任务分配方法 | |
CN105447567B (zh) | 基于bp神经网络与mpso算法的铝电解节能减排控制方法 | |
CN110109358B (zh) | 一种基于反馈的混合多智能体协同控制方法 | |
Che et al. | An enhanced seagull optimization algorithm for solving engineering optimization problems | |
CN103927580A (zh) | 一种基于改进人工蜂群算法的工程约束参数优化方法 | |
CN109685264A (zh) | 火电机组运行优化方法、装置及计算机设备 | |
CN102181592A (zh) | 基于多点雷达数据的无料钟高炉炉顶布料闭环控制方法 | |
CN104268712A (zh) | 基于改进的混合多种群进化算法的能源平衡与调度方法 | |
Li et al. | A two-archive algorithm with decomposition and fitness allocation for multi-modal multi-objective optimization | |
CN109740724A (zh) | 基于参考方向的代理差分演化算法 | |
CN109886448A (zh) | 采用变学习率bp神经网络以及nsga-ii算法的热泵多目标优化控制方法 | |
CN111046559A (zh) | 一种基于双层精英协同进化的多目标优化方法 | |
CN108062583A (zh) | 一种面向节能减排的加热炉工艺参数寻优方法 | |
Zhai et al. | A two-layer algorithm based on PSO for solving unit commitment problem |
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 | ||
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: 20140101 Termination date: 20180425 |