• 《工程索引》(EI)刊源期刊
  • 综合性科学技术类中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于改进差分进化算法的加热炉调度方法

闫祺 李文甲 王稼晨 马凌 赵军

闫祺, 李文甲, 王稼晨, 马凌, 赵军. 基于改进差分进化算法的加热炉调度方法[J]. 工程科学学报. doi: 10.13374/j.issn2095-9389.2020.02.19.004
引用本文: 闫祺, 李文甲, 王稼晨, 马凌, 赵军. 基于改进差分进化算法的加热炉调度方法[J]. 工程科学学报. doi: 10.13374/j.issn2095-9389.2020.02.19.004
YAN Qi, LI Wen-jia, WANG Jia-chen, MA Ling, ZHAO Jun. Reheat furnace production scheduling based on the improved differential evolution algorithm[J]. Chinese Journal of Engineering. doi: 10.13374/j.issn2095-9389.2020.02.19.004
Citation: YAN Qi, LI Wen-jia, WANG Jia-chen, MA Ling, ZHAO Jun. Reheat furnace production scheduling based on the improved differential evolution algorithm[J]. Chinese Journal of Engineering. doi: 10.13374/j.issn2095-9389.2020.02.19.004

基于改进差分进化算法的加热炉调度方法

doi: 10.13374/j.issn2095-9389.2020.02.19.004
基金项目: 国家重点研发计划资助项目(2018YFB0605901)
详细信息

Reheat furnace production scheduling based on the improved differential evolution algorithm

More Information
  • 摘要: 提出一种以燃料消耗量最小为优化目标的加热炉生产调度新方法。首先基于热力学第一定律分析了流入及流出加热炉的各项能量,并对燃料消耗量的计算式进行了理论推导。进而根据加热炉区实际生产调度特点归纳各约束条件,以多台加热炉总燃料消耗量最小为优化目标,构建调度优化数学模型。采用自适应差分进化算法搭配禁忌搜索算法进行综合求解,并通过9组实际钢坯生产案例模拟验证了该算法的可行性和有效性。同时,为了探究加热炉燃料消耗量的影响因素,提出了分别衡量加热炉区缓冲等待、炉内加热两部分时间同理想生产时间匹配程度的评价参数μ1μ2,并分析了燃料消耗量对二者的敏感性,结果表明:当连铸坯到达加热炉节奏与热轧工序出坯节奏之比由0.5增至2时,燃料消耗量对两评价参数的敏感性逐渐减弱。
  • 图 1  加热炉生产流程

    Figure 1.  Reheat furnace production process

    图 2  加热炉区钢坯调度甘特图[3]

    Figure 2.  Billet scheduling Gantt chart in reheat furnace area[3]

    图 3  加热炉能量流模型

    Figure 3.  Energy flow model of reheat furnace

    图 4  差分进化算法框架

    Figure 4.  Differential evolutionary algorithm framework

    图 5  ζ=0.5时不同惩罚因子组合下燃料消耗量变化。(a)三维图;(b)C1C2平面投影图

    Figure 5.  Fuel consumption under the different punishment factor combinations when ζ = 0.5: (a) three-dimensional figure; (b) C1 and C2 plane projection

    图 6  ζ=1时不同惩罚因子组合下燃料消耗量变化。(a)三维图;(b)C1C2平面投影图

    Figure 6.  Fuel consumption under the different punishment factor combinations when ζ = 1: (a) three-dimensional figure; (b) C1 and C2 plane projection

    图 7  ζ=2时不同惩罚因子组合下燃料消耗量变化。(a)三维图;(b)C1C2平面投影图

    Figure 7.  Fuel consumption under the different punishment factor combinations when ζ = 2: (a) three-dimensional figure; (b) C1 and C2 plane projection

    图 8  燃料消耗量随μ1 (a)和μ2 (b)变化关系

    Figure 8.  Relationship between fuel consumption and μ1 (a) and μ2 (b)

    表 1  本文算法与其他算法燃耗对比结果

    Table 1.  Fuel consumption comparison between the proposed algorithm and other algorithms

    TestBillet numberDE/rand/1 fuel consumption/(104 m3)DE/best/1 fuel consumption/(104 m3)DE/current-to-best/1 fuel consumption/(104 m3)Proposed algorithm fuel consumption/(104 m3)
    Test 1823.45653.45993.48723.2736
    Test 2904.18494.21244.16673.8481
    Test 3985.03694.9184.94264.7187
    Test 41075.86065.78855.81025.1345
    Test 51186.12295.91446.01775.8095
    Test 6835.64685.65645.69755.0639
    Test 7947.10357.04697.03916.9521
    Test 81068.31898.17088.23647.8882
    Test 912010.243810.380710.469610.1119
    下载: 导出CSV

    表 2  本文目标与其他调度目标燃耗对比结果

    Table 2.  Fuel consumption comparison between the proposed objective and other objectives

    TestBillet numberFuel consumption/(104 m3)
    Proposed objectiveObjective 1Objective 2Objective 3
    Test 1823.23745.06644.39644.1741
    Test 2903.80584.95675.364.3039
    Test 3984.62165.37265.43194.9666
    Test 41075.07626.23966.01145.2608
    Test 51185.74527.52457.04765.9771
    Test 6835.00787.27397.21765.6939
    Test 7946.87588.67388.79917.8711
    Test 81067.75639.09668.09228.3118
    Test 91209.972610.84310.87410.53
    Average1005.78877.22757.02566.3425
    下载: 导出CSV
  • [1] Ma S H, Wen Z G, Chen J N, et al. Mode of circular economy in China’s iron and steel industry: a case study in Wu’an city. J Clean Prod, 2014, 64: 505 doi:  10.1016/j.jclepro.2013.10.008
    [2] Lu B, Chen G, Chen D M, et al. An energy intensity optimization model for production system in iron and steel industry. Appl Therm Eng, 2016, 100: 285 doi:  10.1016/j.applthermaleng.2016.01.064
    [3] Tang L X, Ren H Z, Yang Y. Reheat furnace scheduling with energy consideration. Int J Prod Res, 2015, 53(6): 1642 doi:  10.1080/00207543.2014.919418
    [4] 杨业建, 姜泽毅, 张欣欣. 钢坯热轧加热炉区生产调度模型与算法. 北京科技大学学报, 2012, 34(7):841

    Yang Y J, Jiang Z Y, Zhang X X. Model and algorithm of furnace area production scheduling in slab hot rolling. J Univ Sci Technol Beijing, 2012, 34(7): 841
    [5] 屠乃威, 罗小川, 柴天佑. 基于蚁群优化算法的步进式加热炉调度. 东北大学学报: 自然科学版, 2011, 32(1):1

    Tu N W, Luo X C, Chai T Y. Scheduling of walking beam reheating furnaces based on ant colony optimization algorithm. J Northeast Univ Nat Sci, 2011, 32(1): 1
    [6] 李铁克, 王柏琳, 赵艳艳. 求解并行加热炉群调度问题的三阶段算法. 系统工程学报, 2011, 26(1):105

    Li T K, Wang B L, Zhao Y Y. Three-stage algorithm for the scheduling problem of parallel reheating furnaces. J Syst Eng, 2011, 26(1): 105
    [7] Ilmer Q, Haeussler S, Missbauer H. Optimal synchronization of the hot rolling stage in steel production. IFAC-PapersOnLine, 2019, 52(13): 1615 doi:  10.1016/j.ifacol.2019.11.431
    [8] Chen D M, Lu B, Zhang X H, et al. Fluctuation characteristic of billet region gas consumption in reheating furnace based on energy apportionment model. Appl Therm Eng, 2018, 136: 152 doi:  10.1016/j.applthermaleng.2018.03.007
    [9] 刘青, 刘倩, 杨建平, 等. 炼钢‒连铸生产调度的研究进展. 工程科学学报, 2020, 42(2):144

    Liu Q, Liu Q, Yang J P, et al. Progress of research on steelmaking‒continuous casting production scheduling. Chin J Eng, 2020, 42(2): 144
    [10] 朱光俊, 曾红, 阮开军, 等. 冶金热工基础. 北京: 冶金工业出版社, 2007

    Zhu G J, Zeng H, Ruan K J, et al. Basis of Metallurgical Thermal Engineering. Beijing: Metallurgical Industry Press, 2007
    [11] Chen D M, Lu B, Dai F Q, et al. Bottleneck of slab thermal efficiency in reheating furnace based on energy apportionment model. Energy, 2018, 150: 1058 doi:  10.1016/j.energy.2018.02.149
    [12] 余志祥. 连铸坯热送热装技术. 北京: 冶金工业出版社, 2002

    Yu Z X. Hot Transfer and Hot Charging Technology of Continuous Casting Billet. Beijing: Metallurgical Industry Press, 2002
    [13] 柴艳红. 钢铁连轧加热工序的节能优化方案设计与应用[学位论文]. 哈尔滨: 哈尔滨工业大学, 2016

    Cai Y H. Design and Application on Optimization of Energy-Saving in Steel Rolling Heating Process[Dissertation]. Harbin: Harbin Institute of Technology, 2016
    [14] 殷瑞钰. 冶金流程工程学. 2版. 北京: 冶金工业出版社, 2009

    Yin R Y. Metallurgical Process Engineering. 2nd Ed. Beijing: Metallurgical Industry Press, 2009
    [15] 夏琼. 基于差分进化算法的加热炉生产过程操作调度[学位论文]. 沈阳: 东北大学, 2010

    Xia Q. Operation and Scheduling of Reheating Furnace Production Process based on Differential Evolution Algorithm[Dissertation]. Shenyang: Northeastern University, 2010
    [16] Bilal, Pant M, Zaheer H, et al. Differential evolution: A review of more than two decades of research. Eng Appl Artif Intell, 2020, 90: 103479 doi:  10.1016/j.engappai.2020.103479
    [17] Marcic T, Stumberger B, Stumberger G. Differential-evolution-based parameter identification of a line-start IPM synchronous motor. IEEE Trans Ind Electron, 2014, 61(11): 5921 doi:  10.1109/TIE.2014.2308160
    [18] Kadhar K M A, Baskar S, Amali S M J. Diversity controlled self adaptive differential evolution based design of non-fragile multivariable PI controller. Eng Appl Artif Intell, 2015, 46: 209 doi:  10.1016/j.engappai.2015.09.015
    [19] Sakr W S, El-Sehiemy R A, Azmy A M. Adaptive differential evolution algorithm for efficient reactive power management. Appl Soft Comput, 2017, 53: 336 doi:  10.1016/j.asoc.2017.01.004
    [20] Wu G H, Mallipeddi R, Suganthan P N, et al. Differential evolution with multi-population based ensemble of mutation strategies. Inform Sci, 2016, 329: 329 doi:  10.1016/j.ins.2015.09.009
    [21] Zhou X G, Zhang G J, Hao X H, et al. Differential evolution with multi-stage strategies for global optimization//2016 IEEE Congress on Evolutionary Computation. Vancouver, 2016: 2550
    [22] Guo Q X, Tang L X. Modelling and discrete differential evolution algorithm for order rescheduling problem in steel industry. Comput Ind Eng, 2019, 130: 586 doi:  10.1016/j.cie.2019.03.011
    [23] Dong Y, Zhao R. Solve train stowage planning problem of steel coil using a pointer-based discrete differential evolution. Int J Prod Res, 2018, 56(22): 6937 doi:  10.1080/00207543.2017.1413260
    [24] Xia Q, Wang X P, Tang L X. Furnace operation optimization with hybrid model based on mechanism and data analytics. Soft Comput, 2019, 23: 9551 doi:  10.1007/s00500-018-3519-9
    [25] Ghosh A, Das S, Chowdhury A, et al. An improved differential evolution algorithm with fitness-based adaptation of the control parameters. Inform Sci, 2011, 181(18): 3749 doi:  10.1016/j.ins.2011.03.010
    [26] Al-Dabbagh R D, Neri F, Idris N, et al. Algorithmic design issues in adaptive differential evolution schemes: review and taxonomy. Swarm Evol Comput, 2018, 43: 284 doi:  10.1016/j.swevo.2018.03.008
    [27] Das S, Mullick S S, Suganthan P N. Recent advances in differential evolution–an updated survey. Swarm Evol Comput, 2016, 27: 1 doi:  10.1016/j.swevo.2016.01.004
    [28] Tang L X, Zhao Y, Liu J Y. An improved differential evolution algorithm for practical dynamic scheduling in steelmaking-continuous casting production. IEEE Trans Evol Comput, 2014, 18(2): 209 doi:  10.1109/TEVC.2013.2250977
    [29] 宁树实, 王伟, 刘全利. 钢铁生产中的加热炉优化调度算法研究. 控制与决策, 2006, 21(10):1138 doi:  10.3321/j.issn:1001-0920.2006.10.012

    Ning S S, Wang W, Liu Q L. An optimal scheduling algorithm for reheating furnace in steel production. Control Decis, 2006, 21(10): 1138 doi:  10.3321/j.issn:1001-0920.2006.10.012
    [30] 谭园园, 宋健海, 刘士新. 加热炉优化调度模型及算法研究. 控制理论与应用, 2011, 28(11):1549

    Tan Y Y, Song J H, Liu S X. Model and algorithm for scheduling reheating furnace. Control Theory Appl, 2011, 28(11): 1549
    [31] 刘健, 刘青, 王彬, 等. 加热炉工序的配置与调度优化研究. 中国冶金, 2014, 24(11):53

    Liu J, Liu Q, Wang B, et al. Research on configuration and scheduling optimization for heating furnace process. China Metall, 2014, 24(11): 53
  • [1] 孔福林, 童莉葛, 魏鹏程, 张培昆, 王立, 吴冰, 陈恩军.  基于粒子群算法的转炉用氧节能优化调度 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2020.04.02.002
    [2] 刘青, 刘倩, 杨建平, 张江山, 高山, 李强笃, 王宝, 王柏琳, 李铁克.  炼钢‒连铸生产调度的研究进展 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2019.04.30.002
    [3] 阙志刚, 吴胜利, 艾仙斌.  基于优化粗粒级固体燃料赋存形态的铁矿烧结过程NOx减排 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2019.02.21.001
    [4] 梁舒, 彭开香.  分布式一致性最优化的梯度算法与收敛分析 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2019.09.05.005
    [5] 刘倩, 杨建平, 王柏琳, 刘青, 高山, 李宏辉.  基于“炉−机对应”的炼钢−连铸生产调度问题遗传优化模型 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2019.08.02.004
    [6] 罗振敏,  刘利涛,  王涛,  张江,  程方明.  温度、其他可燃气体对甲烷爆炸特性及动力学特性研究 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2020.10.22.002
    [7] 刘小娟, 王联国.  一种基于差分进化的正弦余弦算法 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2020.07.26.002
    [8] 苏晓峰, 吴博威, 刘建华, 刘建, 庄昌凌.  微波加热金属液体的实验研究 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2018.12.10.004
    [9] 林歆悠, 夏玉田, 魏申申.  基于增强学习算法的插电式燃料电池电动汽车能量管理控制策略 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2018.10.15.001
    [10] 汪清, 滕立平, 崔鑫, 王阔, 王娴, 李文健.  基于自适应双冗余的TTE调度方法及性能分析 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2019.03.013
    [11] 孟宇, 肖小凤, 赵坤.  基于UWB的地下定位算法和拓扑优化 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2018.06.013
    [12] 李国清, 侯杰, 胡乃联.  地下矿山生产接续与设备调度集成优化模型 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2018.09.005
    [13] 梁江涛, 赵征志, 尹鸿祥, 路洪洲, 陈伟健, 唐荻.  超高强热成形钢的应变速率敏感性 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2018.09.009
    [14] 安富强, 周伟男, 李平.  锂离子电芯用电极对温度与SOC的敏感性 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2018.06.011
    [15] 赵抢抢, 侯保林.  函数型数据分析与优化极限学习机结合的弹药传输机械臂参数辨识 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2017.04.017
    [16] 臧勇, 王远, 秦勤, 管奔.  基于渐进变分法的蛋盒型结构等效刚度分析及优化 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2017.09.012
    [17] 张水平, 王碧, 陈阳.  基于逐层演化的群体智能算法优化 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2017.03.020
    [18] 焦克新, 张建良, 刘征建, 刘彦祥, 梁利生, 贾国栋.  高炉炉缸凝铁层物相分析 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2017.06.004
    [19] 关鸿鹏, 林振娴, 李瑜仙, 刘青, 邢云颖, 王晶, 王修云.  X70管线钢及焊缝在模拟煤制气含氢环境下的氢脆敏感性 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2017.04.008
    [20] 余稳, 向嵩, 石维, 马国强, 于明飞.  20CrMnTi齿轮钢的点蚀敏感性及裂纹萌生风险 . 工程科学学报, doi: 10.13374/j.issn2095-9389.2017.05.011
  • 加载中
图(8) / 表 (2)
计量
  • 文章访问数:  818
  • HTML全文浏览量:  265
  • PDF下载量:  21
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-19
  • 网络出版日期:  2020-12-29

基于改进差分进化算法的加热炉调度方法

doi: 10.13374/j.issn2095-9389.2020.02.19.004
    基金项目:  国家重点研发计划资助项目(2018YFB0605901)
    通讯作者: E-mail:zhaojun@tju.edu.cn
  • 中图分类号: TF087

摘要: 提出一种以燃料消耗量最小为优化目标的加热炉生产调度新方法。首先基于热力学第一定律分析了流入及流出加热炉的各项能量,并对燃料消耗量的计算式进行了理论推导。进而根据加热炉区实际生产调度特点归纳各约束条件,以多台加热炉总燃料消耗量最小为优化目标,构建调度优化数学模型。采用自适应差分进化算法搭配禁忌搜索算法进行综合求解,并通过9组实际钢坯生产案例模拟验证了该算法的可行性和有效性。同时,为了探究加热炉燃料消耗量的影响因素,提出了分别衡量加热炉区缓冲等待、炉内加热两部分时间同理想生产时间匹配程度的评价参数μ1μ2,并分析了燃料消耗量对二者的敏感性,结果表明:当连铸坯到达加热炉节奏与热轧工序出坯节奏之比由0.5增至2时,燃料消耗量对两评价参数的敏感性逐渐减弱。

English Abstract

闫祺, 李文甲, 王稼晨, 马凌, 赵军. 基于改进差分进化算法的加热炉调度方法[J]. 工程科学学报. doi: 10.13374/j.issn2095-9389.2020.02.19.004
引用本文: 闫祺, 李文甲, 王稼晨, 马凌, 赵军. 基于改进差分进化算法的加热炉调度方法[J]. 工程科学学报. doi: 10.13374/j.issn2095-9389.2020.02.19.004
YAN Qi, LI Wen-jia, WANG Jia-chen, MA Ling, ZHAO Jun. Reheat furnace production scheduling based on the improved differential evolution algorithm[J]. Chinese Journal of Engineering. doi: 10.13374/j.issn2095-9389.2020.02.19.004
Citation: YAN Qi, LI Wen-jia, WANG Jia-chen, MA Ling, ZHAO Jun. Reheat furnace production scheduling based on the improved differential evolution algorithm[J]. Chinese Journal of Engineering. doi: 10.13374/j.issn2095-9389.2020.02.19.004
  • 轧钢加热炉是钢铁热轧工序中的重要生产设备,也是钢铁行业中的耗能大户,其能耗约占钢铁生产总能耗的10%~20%[1-2]。降低加热炉能耗对于钢铁企业具有重要意义。加热炉是钢坯热轧前的加热设备,为连铸‒热轧环节起到缓冲作用。由于连铸机、加热炉、热轧机之间生产能力,生产周期存在不协调不匹配的情况,且一个轧制计划中的钢坯规格、质量、入炉温度等参数往往存在差异,如果调度计划制定的不合理,必然会影响上下游工序的匹配衔接。因此优化加热炉生产调度计划对提高生产效率,降低成本具有重要作用。

    近年来,针对加热炉优化调度问题,一些学者提出了行之有效的方法。Tang等[3]以最小化实际在炉时间同理想在炉时间的偏差为目标,采用分散搜索算法求解。杨业建等[4]以极小化板坯温降时间和同炉离散度为目标,采用改进型遗传算法进行求解。屠乃威等[5]以最小化加热炉内冷热板坯混装程度,板坯在炉时间为目标,采用改进蚁群算法求解。李铁克等[6]以最大完工时间和炉内停留时间为优化目标,采用3阶段启发式算法求解。Ilmer等[7]提出一种考虑不同类型加热炉协同生产的调度方法,以最小化钢坯出炉时间和轧机等待时间为目标求解。通过文献分析可以看出,目前对加热炉调度的优化目标大多为缩短某段时间差。实际上,加热炉的生产能耗与入炉钢坯的温度、在炉时间和加热炉产量等诸多因素有关[8],且各因素彼此存在非线性耦合的特征。因此,等待时间总和最小或总完工时间最短这类目标难以完整反映加热炉生产能耗特征,以此目标优化出的调度方案也不一定是能耗最低的生产方案。若推导出对加热炉生产调度方案进行能耗评价的函数表达式,并以此作为优化目标,则其优化出的生产方案在能耗方面定会优于现有方法。

    鉴于以上分析,本文通过热力学第一定律推导出加热炉调度方案的能耗评价函数,并以此作为优化目标汇同生产约束建立调度数学模型,通过多个算例的计算分别测试算法和能耗目标的效果。同时本文考虑到,钢坯入炉前的缓冲等待时间和炉内的保温等待时间是造成能耗高的主要原因,因此提出两个衡量这两类时间相对大小的评价参数,并讨论二者对能耗的影响机制。

    • 轧钢加热炉的生产流程如图1所示,从连铸机至热轧机之间,有4种不同的连接方式:连铸坯直接轧制工艺(HDR),连铸坯直接热装工艺(DHCR),连铸坯热装工艺(HCR)和传统的冷装工艺(CCR)。其中HDR对生产设备和管理水平要求较高,目前国内应用程度较低,其他3种都需要经过加热炉的加热最后送至热轧机。由于上游钢坯来源不同,导致钢坯之间的温度差异性较大,给编排计划的制定带来很大的困难。

      图  1  加热炉生产流程

      Figure 1.  Reheat furnace production process

      加热炉的调度通常是在轧制单元计划,连铸出坯计划已定的条件下,确定待加热的板坯在多台加热炉上的生产顺序,其决策的内容包括计划内所有板坯的入炉号,入炉时间,出炉时间和入炉顺序。在实际生产中,加热炉区的调度计划编排应符合多种约束:钢坯的入炉时间不得早于到达时间;钢坯的出炉顺序遵循轧制计划的顺序;同一加热炉内先轧制的板坯先入炉;每台加热炉有炉容量的限制,任何时间可同时加热的最大钢坯数量是一定的。

      钢铁企业通常根据生产数据的分析总结出不同种类、不同出炉温度、入炉温度的钢坯所需的加热时间,记为标准加热时间。每块钢坯在炉内的加热时间应不小于标准加热时间;同时,钢坯的在炉时间也不应过大,以避免迟滞生产节奏和加剧钢坯表面的氧化烧损,因此钢坯总的加热时间应低于最大住炉时间;钢坯在加热完成后在等待轧制的过程中需待在加热炉内保温;完整的加热炉区钢坯调度计划应如图2所示。

      图  2  加热炉区钢坯调度甘特图[3]

      Figure 2.  Billet scheduling Gantt chart in reheat furnace area[3]

      加热炉生产调度问题属于组合优化难题,由于决策变量维数高,目标函数和约束具有非线性特征,传统的运筹学方法难以快速解决此类问题[9],因此本文选择智能优化算法求解此问题。

    • 加热炉内部结构及各能量交换如图3所示,根据热力学第一定律,建立能量守恒方程式(1)[10-11]。其中,Q1Q5为各项热收入,含义依次为燃料化学热、燃料物理热、空气物理热、钢坯氧化反应热及钢坯带入物理热,Q1'到Q8'为各项热支出,含义依次为钢坯带出物理热、烟气带出物理热、化学不完全燃烧热损失、机械不完全燃烧热损失、炉壁散热、炉门开启辐射热、炉门开启逸气热及冷却系统带出热。

      图  3  加热炉能量流模型

      Figure 3.  Energy flow model of reheat furnace

      $$ \sum\limits_{i=1}^{5} Q_{i}=\sum\limits_{j=1}^{8} Q_{j}^{\prime} $$ (1)

      图3中流入加热炉各项能量Q1Q5由式(2)~(6)计算

      $${Q_1} = B{Q_{{\rm{ne}}}}$$ (2)
      $${Q_2} = B{c_{{\rm{fu}}}}{t_{{\rm{fu}}}}$$ (3)
      $${Q_3} = B{n_{\rm{0}}}{L_0}{c_{{\rm{ai}}}}{t_{{\rm{ai}}}}$$ (4)
      $${Q_4} = 5652Ga$$ (5)
      $${Q_5} = G{c_{{\rm{me}}}}{t_{{\rm{me}}}}$$ (6)

      式中:B为单位燃料消耗量,m3·h−1Qne为燃料的低位发热量,kJ·m−3cfu为燃料平均体积热容,kJ·m−3·℃−1tfu为燃料的温度,℃;n0为空气消耗系数;L0为单位燃料理论所需空气量,m3cai为空气平均体积热容,kJ·m−3·℃−1tai为空气温度,℃;G为加热炉产量,kg·h−1a为烧损率;cme为入炉钢坯的比热容,kJ·kg−1·℃−1tme表示入炉钢坯的平均温度,℃。

      图3中流出加热炉各的项能量Q1'到Q8'由式(7)~(14)计算

      $$Q_1' = Gc_{{\rm{me}}}'t_{{\rm{me}}}'$$ (7)
      $$Q_2' = B{V_{\rm{n}}}{c_{{\rm{wa}}}}{t_{{\rm{wa}}}}$$ (8)
      $$Q_3' = 184.06B{V_{\rm{n}}}\varphi ({\rm{CO}})$$ (9)
      $$Q_4' = KB{Q_{{\rm{ne}}}}$$ (10)
      $$Q_5' = \frac{{({t_1} - {t_2})A}}{{\dfrac{{{\delta _1}}}{{{\lambda _1}}} + \dfrac{{{\delta _2}}}{{{\lambda _2}}} + \cdots + \dfrac{{{\delta _n}}}{{{\lambda _n}}} + {\textit{z}}}}$$ (11)
      $$Q_6' = {\rm{1}}{\rm{.575}}{\left( {\frac{{{t_1}}}{{100}}} \right)^{\rm{4}}}\varphi {A'}\psi $$ (12)
      $$Q_7' = {V_0}{c_{{\rm{lu}}}}{t_{{\rm{lu}}}}\psi $$ (13)
      $$Q_8' = {G_{\rm{s}}}\left( {c_{\rm{s}}^{{\rm{''}}}t_{\rm{s}}^{''} - c_{\rm{s}}^{\rm{'}}t_{\rm{s}}'} \right) + {G_{\rm{g}}}\left( {{H_{\rm{q}}} - c_{\rm{s}}^{\rm{'}}t_{\rm{s}}' - {\rm{0}}{\rm{.01}}{\gamma '}{\omega '}} \right)$$ (14)

      式中:cme'为钢坯的比热容,kJ·kg−3·℃−1tme'为出炉钢坯的平均温度,℃;Vn为单位燃料排出废气量,m3cwa为出炉废气的体积热容,kJ·m−3·℃−1twa为出炉烟气的温度,℃;$ \varphi \left(\mathrm{C}\mathrm{O}\right) $为烟气中CO的体积分数,%;K为机械不完全燃烧系数;t1t2分别为炉内温度和周围环境温度,℃;A为炉壁的面积,m2δ1δ2,...δn为炉壁各层材料厚度,m;λ1λ2,…λn为各层材料的导热系数,kJ·m−1·h−1·K−1z为炉壁外表面与空间的热阻,kJ· m−1·h−1·K−1$ \varphi $为综合角度系数;A'为炉门的面积,m2ψ为炉门开启时间占总生产时间的比例;V0为炉门开启时的逸气量,m3·h−1clu为炉内气体的体积热容,kJ·m−3·K−1tlu为炉内气体的温度,℃;Gs为冷却水流量,t·h−1; Gg为蒸汽流量,t·h−1cs',cs''分别为冷却水进,出口比热容,kJ·kg−1·℃−1ts',ts''分别为进出口水温,℃;Hq为蒸汽的焓值,kJ·kg−1γ'为水的汽化潜热,kJ·kg−1ω'为蒸汽的湿度,%。

      由于钢坯在等候装炉时不断向环境散热,其最终入炉温度通常小于其初始温度,因此通过式(15)[12]计算钢坯的最终入炉温度。

      $${T_{\rm{s}}} = \left( {{T_{{\rm{s0}}}} - {T_{{\rm{a}}0}}} \right){\operatorname{e} ^{ - 0.000028\theta {\tau _{\rm{0}}}/V}} + {T_{{\rm{a}}0}}$$ (15)

      式中:Ts0为初始时刻钢坯的平均温度,℃;θ为钢坯的表面积,m2V为钢坯的体积,m3τ0为钢坯在空气中的停留时间,s;Ta0为空气温度,℃。

    • 选择总燃料消耗量Ctf(m3)作为衡量加热炉能耗高低的指标。为便于叙述,下文简称为“燃耗”。其计算式为(16):

      $${C_{{\rm{tf}}}} = B\Delta \tau $$ (16)

      式中:Δτ为加热炉运行总时间,h。综合(1~16)各式得出加热炉燃耗的计算式(17)

      $${C_{{\rm{tf}}}} = \frac{{\Delta \tau G\left( {c_{{\rm{me}}}'t_{{\rm{me}}}^{\rm{'}} - {c_{{\rm{me}}}}{t_{{\rm{me}}}} - 5652a} \right) + \Delta \tau \left( {1.575{{\left( {\dfrac{{{t_1}}}{{100}}} \right)}^4}\varphi {A'} + {V_0}{c_{{\rm{lu}}}}{t_{{\rm{lu}}}}} \right)\psi + \Delta \tau \left( {Q_5' + Q_8'} \right)}}{{\left( {1 - k} \right){Q_{{\rm{ne}}}} + {c_{{\rm{fu}}}}{t_{{\rm{fu}}}} + {n_{\rm{0}}}{L_0}{c_{{\rm{ai}}}}{t_{{\rm{ai}}}} - {V_n}{c_{{\rm{wa}}}}{t_{{\rm{wa}}}} - 184.06{V_{\rm{n}}}\phi \left( {{\rm{CO}}} \right)}}$$ (17)

      为便于计算和分析,考虑加热炉生产调度计划制定的特点和决策对象的特征,化简(17)式可得:

      $${C_{{\rm{tf}}}} = \Delta \tau {A_1}\left[ {G\left( {{A_2} - c_{{\rm{me}}}^{}t_{{\rm{me}}}^{}} \right) + {A_3}\psi + {A_4}} \right]$$ (18)

      式中,A1A2A3A4为化简后的常数,其表达式为(19)~(22)。

      $$\begin{split} & {A_1} = \left[ \left( {{\rm{1 - }}k} \right){Q_{{\rm{ne}}}} + {c_{{\rm{fu}}}}{t_{{\rm{fu}}}} + {n_{\rm{0}}}{L_0}{c_{{\rm{ai}}}}{t_{{\rm{ai}}}} - {V_{\rm{n}}}{c_{{\rm{wa}}}}{t_{{\rm{wa}}}} -\right. \\ & \left. 184.06{V_{\rm{n}}}\phi \left( {{\rm{CO}}} \right) \right]^{{\rm{ - 1}}} \end{split}$$ (19)
      $${A_2} = c_{{\rm{me}}}'t_{{\rm{me}}}' - 5652a$$ (20)
      $${A_3} = 1.575{\left( {\frac{{{t_1}}}{{100}}} \right)^4}\varphi {A'} + {V_0}c{}_{_{{\rm{lu}}}}{t_{{\rm{lu}}}}$$ (21)
      $${A_4} = Q_5' + Q_8'$$ (22)
    • 假设在轧制计划已知的条件下,选择DHCR和CCR两种加热炉生产调度模式进行研究。考虑到加热炉上游工序繁多,为便于计算,在保证模型的通用性前提下做出适当简化,说明如下:

      (1)调度计划中各钢坯到达加热炉前缓冲区的时刻,温度已知。

      (2)每块钢坯的标准加热时间已知,该值由企业根据以往生产数据总结而成[13],钢种,目标出炉温度等已知。

      (3)忽略加热炉故障导致的停产等突发情况,认为炉体始终保持稳定生产状态。

      (4)加热炉装钢,出钢操作时间为定值。

    • 为便于描述,定义以下符号和参数

      (1)索引与集合。S为加热炉生产调度计划中的钢坯集合,S={1, 2,…, i,…, n},其中n表示钢坯数。L表示加热炉集合,L={1, 2,…, k,…, l},其中l表示加热炉数量。

      (2)问题参数。mi表示钢坯i的质量。phi表示钢坯i在加热炉中标准加热时间。Zi表示钢坯i到达加热炉前缓冲区时刻。pdi表示钢坯i在加热炉中的最大住炉时间。pri表示钢坯i在热轧机上的轧制时间。T1i表示钢坯i从等候区域输送至加热炉的传输时间。T2i表示钢坯i从加热炉输送至热轧机的传输时间。pk表示第一块进入加热炉k的钢坯序号。qk表示最后一块离开加热炉k的钢坯序号。tbi表示钢坯i入炉时温度。Ii表示在同一加热炉中,排在钢坯i后第C块(C为加热炉炉容)入炉加热的钢坯号。ε表示热轧机最大空转等待时间。δmin表示同一加热炉加热的相邻钢坯的最小入炉时间间隔。Y表示一个足够大正数。Δt0表示加热炉装钢,出钢的操作时间。

      (3)决策变量。γik为二元变量,如果钢坯i在加热炉k中加热该值记为1,其他情形为0。bi表示钢坯i的入炉时间。di表示钢坯i的出炉时间。Oijk为三元变量,如果钢坯j紧随钢坯i装入加热炉k中加热该值为1,其他情形为0。

      基于上述参数,变量的定义,构建加热炉调度的整体数学模型:

      $$\begin{split} & {\rm{Min }}{A}_{1}\left[ \sum\limits _{i\in S}{m}_{i}{A}_{2}- \sum \limits_{i\in S}{m}_{i}{c}_{\rm{me},}{}_{{\rm{tb}}_{i}}\cdot {\rm{tb}}_{i}+2{A}_{3}n\Delta {t}_{0}+\right. \\ & \left. {A}_{4}{ \sum\limits _{k\in L}{ \sum \limits_{i\in S}{\gamma }_{ik}({d}_{{q}_{k}}-{b}_{{p}_{k}})}}\right] \end{split}$$ (23)
      $$ { \sum \limits_{k\in L}{\gamma }_{ik}=1,\forall i\in S}$$ (24)
      $$ {\rm{ph}}_{i}\leqslant {d}_{i}-{b}_{i}\leqslant {\rm{pd}}_{i},\forall i\in S$$ (25)
      $$ 0\leqslant {d}_{i+1}-{d}_{i}-{\rm{pr}}_{i}\leqslant \varepsilon ,\forall i\in S$$ (26)
      $$ {d}_{j}-{d}_{i}-{T}_{2}+Y\left(3-{\gamma }_{ik}-{\gamma }_{jk}-{O}_{ijk}\right)\geqslant 0,\forall i,j\in S,\forall k\in L$$ (27)
      $$ {b}_{j}-{b}_{i}-{T}_{1}+Y\left(\!3-{\gamma }_{ik}-{\gamma }_{jk}-{O}_{ijk}\!\right)\!\geqslant\! {\delta }_{\mathrm{min}},\forall i,j\!\in\! S,\forall k \!\in \! L$$ (28)
      $$ {b}_{{I}_{i}}\geqslant {d}_{i},\forall i\in S$$ (29)
      $$ {\gamma }_{ik}\in \left\{0,1\right\},\forall i\in S,\forall k\in L$$ (30)
      $$ {O}_{ijk}\in \left\{0,1\right\},\forall i,j\in S,\forall k\in L$$ (31)

      目标函数(23)表示多台加热炉总燃料消耗量最小。约束(24)表示每块钢坯只能分配给一台加热炉加热。约束(25)表示每块钢坯在加热炉内的加热时间应不低于标准加热时间且不高于最大住炉时间。约束(26)表示板坯在热轧机组空闲前需一直在炉内停留且热轧机组停机等待板坯的时间不能过长。约束(27)表示在同一加热炉加热的钢坯按照轧制顺序出炉。约束(28)表示在同一加热炉加热的钢坯按轧制顺序入炉。约束(29)表示每台加热炉同时加热的钢坯数不得超过其炉容量。约束(30)和约束(31)表示决策变量γikOijk的取值范围。

    • 分析每块钢坯在加热炉区停留的时间可知,由于加热炉的出钢顺序受轧制计划的约束,因此每块钢坯从出炉至输送到热轧机上的时间间隔较固定;但另外两部分时间(入炉前等待时间,炉内加热时间)由于受到上游连铸出坯时间、加热炉炉容量、轧制计划等制约,往往波动性较大。因此有必要探究这两部分时间的变化同加热炉燃料消耗量的关系。殷瑞钰[14]提出用“连续化程度”这一物流指标来衡量各工序匹配程度,本文考虑加热炉调度实际,引入加热匹配系数μ1和等待匹配系数μ2两个评价参数,分别衡量一个调度计划内所有钢坯炉内加热时间和炉前等待时间同理想生产时间的差距,并分析二者对燃料消耗量的影响规律,μ1μ2计算公式分别为式(32)和(32)。

      $${\mu _{\rm{1}}} = \frac{{\displaystyle\sum\limits_{i \in S} {\left( {{d_i} - {b_i}} \right)} }}{{\displaystyle\sum\limits_{i \in S} {{\rm{p}}{{\rm{h}}_i}} }}$$ (32)
      $${\mu _{\rm{2}}} = \frac{{\displaystyle\sum\limits_{i \in S} {\left( {{b_i} - {Z_i}} \right)} }}{{\displaystyle\sum\limits_{i \in S} {{\rm{p}}{{\rm{h}}_i}} }}$$ (33)
    • 差分进化算法(Differential evolution algorithm)是一种随机的启发式搜索算法,凭借其编程简便,收敛速度快,寻优能力强[15]等优势,已连续多年在进化计算大会(CEC)竞赛中跻身前3名[16],并且在模式识别[17]、工业控制[18-19]、函数优化[20-21]等诸多领域获得广泛应用。特别是在钢铁生产领域,差分进化算法已被成功用于解决企业订单重调度[22]、钢卷配载规划[23]、加热炉操作优化[24]等实际问题。

      为避免差分进化算法的种群在迭代过程中易出现的早熟现象,减轻算法参数调试的工作量,Ghosh等[25]在差分进化算法(DE)的基础上开发了FiADE(Fitness-adaptive DE)算法,它可根据每代种群个体的适应度函数值变化,对DE算法的缩放因子、交叉概率动态调整,平衡算法的开发和探索能力[26]。本文以FiADE算法作为主体求解算法,同时为提高算法寻优效率,对种群中表现较优的个体进一步执行禁忌搜索。算法整体流程如图4所示。其主要步骤如下:

      图  4  差分进化算法框架

      Figure 4.  Differential evolutionary algorithm framework

      (1)初始化种群。

      种群采用实数编码,随机产生NPD维解向量Xi={xi1$,\cdots $ xiD},i = 1$,\cdots $,NPD为决策变量个数,NP为种群数。各决策变量分布在其上边界Xmin={xmin1$, \cdots $xminD}和下边界Xmax={xmax1$,\cdots $ xmaxD}构成的解空间中,计算公式如(34)所示:

      $$x_{i,0}^j = x_{\min }^j + {\rm{rand}}\left( {0,1} \right) \cdot \left( {x_{\max }^j - x_{\min }^j} \right)$$ (34)

      rand(0, 1)指[0, 1]之间服从均匀分布的随机数,xi,0j中的0指此时种群处于初始代。

      (2)变异操作。

      对种群中的每个目标向量Xi,EE表示算法运行至第E代)执行差分变异操作(35),生成的变异向量记为Vi,E=[vi,G1$, \cdots $vi,GD]

      $${V_{i,E}} = {X_{r_1^i,E}} + {F_i} \cdot \left( {{X_{r_2^i,E}} - {X_{r_3^i,E}}} \right) + {F_i} \cdot \left( {{X_{r_4^i,E}} - {X_{r_5^i,E}}} \right)$$ (35)

      式中,下标r1ir2ir3ir4ir5i为随机抽取的位于[1, NP]之间的不同于i的个体序号,Fi为个体i所对应的缩放因子,其计算流程如下

      Step1:计算F1,对于每个个体Xi

      $${F_1} = {F_{\max }} \cdot \left( {\frac{{\Delta {f_i}}}{{{\lambda _{{\rm{con}}}} + \Delta {f_i}}}} \right)$$ (36)

      式中,∆fi=|f(Xi)‒f(Xbest)|,Xbest为当前种群中适应度最低的个体,λcon = 0.1·∆fi + 10−14Fmax为给定的参数。

      Step2:计算F2,对于每个个体Xi

      $${F_2} = {F_{\max }} \cdot \left( {1 - {{\rm{e}}^{ - \Delta {f_i}}}} \right)$$ (37)

      Step3:计算Fi

      $${F_i} = \max \left( {{F_1},{F_2}} \right)$$ (38)

      (3)交叉操作。

      将目标向量Xi,E和其对应的变异向量Vi,E执行交叉策略,新生成的个体记为Ui,E={ui,E1ui,E2$, \cdots $,ui,ED},个体i中第j个点ui,Ej计算公式为

      $$u_{i,G}^j = \left\{ \begin{array}{l} v_{i,E}^j\begin{array}{*{20}{c}} {}&{{\rm{if}}\mathop {}\limits_{} {\rm{rand}}\left( {0,1} \right) \leqslant {\rm{C}}{{\rm{R}}_i}} \end{array} \\ x_{i,E}^j\begin{array}{*{20}{c}} {}&{{\rm{otherwise}}} \end{array} \\ \end{array} \right.$$ (39)

      式中,CRi为个体i所对应的交叉概率,其计算流程如下:

      比较变异向量Vi和原始向量Xbest的适应度函数值,如果f(Vi) < f(Xbest),执行(40)式,否则执行(41)式。

      $${\rm{C}}{{\rm{R}}_i} = {\rm{C}}{{\rm{R}}_{{\rm{const}}}}$$ (40)
      $${\rm{C}}{{\rm{R}}_i} = {\rm{C}}{{\rm{R}}_{\min }} + \frac{{{\rm{C}}{{\rm{R}}_{\max }} - {\rm{C}}{{\rm{R}}_{\min }}}}{{1 + \left| {f\left( {{{{V}}_i}} \right) - f\left( {{{{X}}_{{\rm{best}}}}} \right)} \right|}}$$ (41)

      式中,CRconst,CRmin,CRmax为设定的参数。

      (4)选择操作。

      比较原始个体Xi,E和经变异,交叉操作后所得新个体Ui,E的适应度函数值,以(42)式进行选择,优者进入下一轮循环,弱者淘汰出局。

      $${{{X}}_{i,E + 1}} = \left\{ \begin{array}{l} {{{U}}_{i,E}}\begin{array}{*{20}{c}} {}&{{\rm{if}}\begin{array}{*{20}{c}} {f\left( {{{{U}}_{i,E}}} \right) \leqslant f\left( {{{{X}}_{i,E}}} \right)} \end{array}} \end{array} \\ {{{X}}_{i,E}}\begin{array}{*{20}{c}} {}&{{\rm{otherwise}}} \end{array} \\ \end{array} \right.$$ (42)

      (5)禁忌搜索。

      选择操作结束后,对新生成的子代个体,按照适应度函数值由小到大进行排序,并选择前10%的优良个体执行禁忌搜索。具体步骤如下:

      步骤1:设定算法运行的各输入参数,即候选解个数、禁忌表长度、迭代次数、初始解等,初始禁忌表为空,同时将初始解设为历史最优解。

      步骤2:对当前解进行局部搜索生成候选解,并计算每个候选解的适应度函数值。

      步骤3:判断每个候选解是否属于禁忌且选出适应度函数值最优的候选解。

      步骤4:比较最优候选解与历史最优解的适应值大小,若优于历史最优解则执行步骤5,否则执行步骤6。

      步骤5:忽略最优候选解的禁忌状态,将其作为新的历史最优解和当前解,加入禁忌表中,并更新禁忌表的状态,执行步骤7。

      步骤6:在候选解中选择非禁忌的最优适应度函数值的解,将其作为当前解并加入禁忌表中,更新禁忌表的状态,执行步骤7。

      步骤7:判断是否达到最大迭代次数,若达到执行步骤8,否则迭代次数加一并返回步骤2。

      步骤8:输出历史最优解。

    • 为验证本文所提方法在降低加热炉生产能耗方面的效果,以国内某热轧厂生产线实际数据建立数学模型,该模型包含3台并行布置且容量相同的加热炉,下游配置1台热轧机,上游钢坯来源为连铸机和板坯库。选取9组轧制作业计划中的钢坯数据作为计算案例,其中前5个算例采用DHCR调度方式,后四个算例为DHCR与CCR耦合的调度方式。

      选取国内某加热炉热平衡测试时的数据,其燃料低位发热量为33812 kJ·m−3,理论空气需要量8.39 m3,空气消耗系数为1.5,单位燃料废气量为11.03 m3,机械不完全燃烧损失系数为0.01,烟气中CO体积分数为1%,烧损率为0.03,冷却水流量250 t·h−1。为便于计算,各处温度取加热一批钢坯时间内的平均值,燃料温度307 ℃,空气温度355 ℃,出炉烟气温度833 ℃,炉门开启逸出炉气温度为787 ℃。

      算法参数设置方面,设本文差分进化算法的种群个体数为100,迭代次数为200,Fmax为0.8,CRconst为0.9,CRmin为0.1,CRmax为0.8;禁忌搜索算法的候选解长度50,禁忌列表长度20,迭代次数20。

    • 选择LENOVO510S(Intel Core i5-7200U CPU@2.50GHZ,4 GB RAM)作为测试设备。为测试本文改进DE算法的效果,选择DE/rand/1,DE/best/1,DE/current-to-best/1 3种典型差分算法[27]作为对比算法,3种算法的种群数,迭代次数与本文算法一致,缩放因子F和交叉概率CR通过4水平全因子设计实验[28]获得,它们的参数分别设置为:DE/rand/1:F=0.4、CR=0.7;DE/best/1:F=0.7、CR=0.4;DE/current-to-best/1:F=0.6、CR=0.9。当算法运行停止返回出最终的调度方案后,依据式(16)计算各方案对应的加热炉燃耗。每种算法均运行10次取平均值,计算结果见表1。由表1可以看出,在所有测试算例中,本文所提改进差分进化算法均获得最低燃耗的生产方案,这验证了本算法的有效性。

      表 1  本文算法与其他算法燃耗对比结果

      Table 1.  Fuel consumption comparison between the proposed algorithm and other algorithms

      TestBillet numberDE/rand/1 fuel consumption/(104 m3)DE/best/1 fuel consumption/(104 m3)DE/current-to-best/1 fuel consumption/(104 m3)Proposed algorithm fuel consumption/(104 m3)
      Test 1823.45653.45993.48723.2736
      Test 2904.18494.21244.16673.8481
      Test 3985.03694.9184.94264.7187
      Test 41075.86065.78855.81025.1345
      Test 51186.12295.91446.01775.8095
      Test 6835.64685.65645.69755.0639
      Test 7947.10357.04697.03916.9521
      Test 81068.31898.17088.23647.8882
      Test 912010.243810.380710.469610.1119

      同时,考虑到优化目标对于调度模型的重要性,为验证文中所提燃耗目标的有效性,以总住炉等待时间[29],总完工时间[30],总热轧机[3]等待时间作为对比目标,3种目标的数学表达式为式(43)~(45),为便于表述,三种目标依次记为目标1,目标2,目标3。采用本文提出的改进差分进化算法作为求解算法,对每组算例用燃耗目标和3种对比目标分别进行求解,每种调度目标计算时均按上文所设参数连续运行10次,取10次燃耗结果的最小值,计算结果见表2

      表 2  本文目标与其他调度目标燃耗对比结果

      Table 2.  Fuel consumption comparison between the proposed objective and other objectives

      TestBillet numberFuel consumption/(104 m3)
      Proposed objectiveObjective 1Objective 2Objective 3
      Test 1823.23745.06644.39644.1741
      Test 2903.80584.95675.364.3039
      Test 3984.62165.37265.43194.9666
      Test 41075.07626.23966.01145.2608
      Test 51185.74527.52457.04765.9771
      Test 6835.00787.27397.21765.6939
      Test 7946.87588.67388.79917.8711
      Test 81067.75639.09668.09228.3118
      Test 91209.972610.84310.87410.53
      Average1005.78877.22757.02566.3425
      $${\rm{Min}}\sum\limits_{i \in S} {\left( {{d_i} - {b_i} - {\rm{p}}{{\rm{h}}_i}} \right)} $$ (43)
      $${\rm{Min}}\sum\limits_{i \in S} {\left( {{d_i} - {b_i}} \right)} $$ (44)
      $${\rm{Min}}\sum\limits_{i \in S} {\left( {{d_{i + 1}} - {d_i} - {\rm{p}}{{\rm{r}}_i}} \right)} $$ (45)

      表2可以看出,在所有测试算例中,依据本文所提目标求得的调度方案,其加热炉燃料消耗量最低。根据所有算例测试结果的平均值表明,若依本文所提优化目标进行生产调度,则理论上加热炉能耗将比调度目标1降低约20%,比目标2降低约17.6%,比目标3降低约8.7%,可见本文提出的调度目标在降低加热炉能耗方面具有一定优势。

    • 现通过算例1中的钢坯数据,探究加热炉燃耗对工序匹配评价参数μ1μ2的敏感性,由于μ1μ2并非完全独立,故同时考虑二者对能耗的影响。设计新的优化目标(46),其中C1C2分别为μ1μ2的惩罚因子,C1表示优化目标中μ1的重要程度,C2表示优化目标中μ2的重要程度。C1C2取值区间均设为[0,1],间隔为0.1,共考虑100种不同的惩罚因子组合。设计优化目标(46)的目的是通过算法寻优,使得调度计划对应的μ1μ2不断下降,并观察对应的燃耗变化情况以分析两参数同燃耗的敏感性;当C1大于C2时,最终两参数的优化程度通常为μ1>μ2,因此通过不同的C1C2组合,获得μ1μ2下降幅度不同的调度方案,可进一步观察二者对燃耗的影响程度。算法参数保持上文中的设置不变。

      $$ {\rm{Min}}\left(C_{1} \mu_{1}+C_{2} \mu_{2}\right) $$ (46)

      由于加热炉区的生产调度会受到连铸工序的钢坯到达节奏和下游热轧机生产节奏的影响,因此本文考虑划分连铸坯到达节奏大于,等于和小于热轧机出坯节奏3种情形讨论,将一块连铸坯经辊道运输至加热炉的平均时间定义为连铸坯到达节奏β1 min,设定热轧机出坯节奏记为β2 min[31],定义连铸坯到达节奏与轧制节奏之比ζ=β1/β2。将轧机节奏设为1 min,改变β1以构建3种不同ζ的情形,讨论当ζ为0.5、1、2时加热炉燃耗对μ1μ2的敏感性。其计算结果分别如图57所示。对比图5(a)图6(a)图7(a)可知,当ζ从0.5逐渐增至2时,加热炉生产燃耗的波动区间呈逐渐减小的趋势,燃耗对μ1μ2的敏感性逐渐减弱。

      图  5  ζ=0.5时不同惩罚因子组合下燃料消耗量变化。(a)三维图;(b)C1C2平面投影图

      Figure 5.  Fuel consumption under the different punishment factor combinations when ζ = 0.5: (a) three-dimensional figure; (b) C1 and C2 plane projection

      图  6  ζ=1时不同惩罚因子组合下燃料消耗量变化。(a)三维图;(b)C1C2平面投影图

      Figure 6.  Fuel consumption under the different punishment factor combinations when ζ = 1: (a) three-dimensional figure; (b) C1 and C2 plane projection

      图  7  ζ=2时不同惩罚因子组合下燃料消耗量变化。(a)三维图;(b)C1C2平面投影图

      Figure 7.  Fuel consumption under the different punishment factor combinations when ζ = 2: (a) three-dimensional figure; (b) C1 and C2 plane projection

      图5(b)图6(b)可知,当上游连铸坯到达加热炉节奏快于或等于轧机的生产节奏时,往往当C2大于C1时,对应燃耗更低。相反,当C2小于C1时,燃耗往往偏高。这表明,若仅考虑获得能耗最低的生产调度方案,降低μ2比降低μ1更重要,即降低钢坯在入炉前的等待时间,提高入炉温度,相对于减少钢坯在炉内的加工时间“收益”更大。

      图7(b)可知,当ζ=2时,无论怎么改变C1C2的搭配,能耗变化非常小。这主要是由于当ζ>1时,连铸坯到达节奏慢于轧机生产节奏,使得连铸工序物质流通量低于下游热轧工序物质流通量。此工况下钢坯在加热炉区的停留时间接近其标准加热时间,加热炉充当“缓冲活套”的作用被削弱。从生产实际角度分析,当上游物质流通量低于下游物质流通量时,可能会造成轧机空转等待,影响热轧机的正常运行。为保证连铸至热轧区间的稳定生产,应避免出现这种工况。

      通过敏感性分析可知,当ζ ≤1时,能耗对μ1μ2的变化敏感。然而两参数对能耗的具体影响规律需进一步讨论,因此本文以算例1的钢坯信息为例,并设定ζ=1。选择本文目标优化出的调度方案作为初始点(initial point),该方案对应燃料消耗量32374 m3μ1为1.51,μ2为0.036。依次以最大化μ1(记为优化方向1)和最小化μ1(记为优化方向2)为优化目标进行寻优计算,并保存算法每代最优解,算法迭代次数设为2000。之后计算每代最优解对应的μ1μ2,燃耗;汇总并绘制出燃耗随μ1μ2变化的曲线如图8所示。

      图  8  燃料消耗量随μ1 (a)和μ2 (b)变化关系

      Figure 8.  Relationship between fuel consumption and μ1 (a) and μ2 (b)

      综合分析图8可以看出,当调度方案对应μ1<1.51或μ2>0.036时,加热炉燃耗随着μ1的减小而增大,随着μ2的减小而减小,这表明在此情形下,由工序间节奏不匹配引起的钢坯入炉前的等待时间,成为了加热炉能耗高的主要矛盾,应设法降低所有钢坯的炉前等待时间提升钢坯的平均入炉温度。反之,当μ1>1.51或μ2<0.036时,加热炉燃耗随着μ1的增大先增大后趋于稳定,随着μ2的增大而减小,这表明在此情形下,钢坯在炉内的停留时间过长成了主要矛盾,应设法加快钢坯在加热炉区间的生产节奏。

    • (1)以加热炉燃料消耗量最小为优化目标,建立了加热炉区生产调度数学模型,选用自适应差分进化算法搭配禁忌搜索算法对问题进行求解。为了探究加热炉燃料消耗量的影响因素,提出了两个衡量钢坯炉前等待和炉内生产两部分时间同理想生产时间匹配程度的评价参数μ1μ2

      (1)九组算例的模拟结果表明:本文提出的改进自适应差分进化算法的性能优于经典差分进化算法;以燃料消耗量最小为目标进行优化调度理论上可比其他目标降低能耗8.7%至20%。

      (2)μ1μ2同能耗的敏感性分析表明,当连铸坯到达节奏与热轧出坯节奏之比ζ由0.5增至2时,燃料消耗量对两参数的敏感性由强变弱。且当ζ≤1时,降低μ2比降低μ1更能有效降低加热炉能耗。

      (3)通过对本文算例1的模拟计算发现,当ζ=1时,其最优调度方案所需燃料消耗量为32374 m3,此时μ1为1.51,μ2为0.036。对于其他调度方案,当其μ1<1.51或μ2>0.036时,燃料消耗量随μ1的减小而增大,随μ2的减小而减小。反之,燃料消耗量随着μ1的增大先增大后趋于稳定,随着μ2的增大而减小。

参考文献 (31)

目录

    /

    返回文章
    返回