CN106971243B - 一种降低生产用电成本的热轧优化调度方法 - Google Patents
一种降低生产用电成本的热轧优化调度方法 Download PDFInfo
- Publication number
- CN106971243B CN106971243B CN201710198430.2A CN201710198430A CN106971243B CN 106971243 B CN106971243 B CN 106971243B CN 201710198430 A CN201710198430 A CN 201710198430A CN 106971243 B CN106971243 B CN 106971243B
- Authority
- CN
- China
- Prior art keywords
- rolling
- electricity
- production
- time
- rolling unit
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 99
- 238000005098 hot rolling Methods 0.000 title claims abstract description 31
- 238000005457 optimization Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005096 rolling process Methods 0.000 claims abstract description 185
- 230000005611 electricity Effects 0.000 claims abstract description 106
- 210000000349 chromosome Anatomy 0.000 claims abstract description 20
- 230000035772 mutation Effects 0.000 claims abstract description 16
- 230000002068 genetic effect Effects 0.000 claims abstract description 10
- 101100348084 Drosophila melanogaster CDase gene Proteins 0.000 claims description 33
- 238000004364 calculation method Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000006798 recombination Effects 0.000 claims description 9
- 238000005215 recombination Methods 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 6
- 238000007667 floating Methods 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 2
- 230000003044 adaptive effect Effects 0.000 abstract description 3
- 229910000831 Steel Inorganic materials 0.000 description 6
- 239000010959 steel Substances 0.000 description 6
- 238000005265 energy consumption Methods 0.000 description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000004134 energy conservation Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 230000033764 rhythmic process Effects 0.000 description 1
Images
Classifications
-
- 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"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Genetics & Genomics (AREA)
- Manufacturing & Machinery (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Primary Health Care (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- Educational Administration (AREA)
- Control Of Metal Rolling (AREA)
Abstract
本发明公开了一种降低生产用电成本的热轧优化调度方法,包括以下步骤:加载分时电价费率数据、轧制速度设定数据和待排产的板坯数据;确定最小化生产用电成本的热轧优化调度的目标函数和约束条件;批量计划编制、轧制单元待产时间和轧制速度信息的混合染色体编码;确定混合染色体编码的选择、交叉、变异算子及其操作;执行非支配排序多目标遗传算法优化热轧调度目标函数;解集的多目标决策及推荐方案输出。本发明在批量计划编制时同步构造与电价时段相适应的生产负荷单元,通过合理安排轧制单元待产时间以及调整板坯的轧制速度进行生产负荷的时间转移,以错峰和避峰生产方式组织生产,优化相邻板坯间惩罚值和生产用电成本,促进了用电负荷削峰。
Description
技术领域
本发明涉及冶金领域,特别涉及一种降低生产用电成本的热轧优化调度方法。
背景技术
钢铁行业是我国十大重点能耗行业之首,由于电力成本居高不下,钢铁生产节能降耗需求迫切。作为钢铁生产的关键工序,热轧过程能耗较大且可调度性强,是节能降耗重点关注环节,特别在考虑峰值负荷管理和分时电价下的分时段电力负荷管理时,热轧生产的优化调度对于实现负荷优化分配和降低电力成本尤为重要。
热轧生产调度主要通过批量计划进行组织,批量计划编制结果直接影响产品质量、生产效率、生产成本等。考虑分时电价环境,通过批量计划编制优化降低生产用电成本在已有的文献和技术实现中均有涉及。以热轧生产过程为对象实施经济负荷调度,在批量计划编制时,同步构造与电价时段相适应的生产负荷单元,在轧制生产能力裕量下安排轧制生产的待产时间进行自身生产负荷的时间转移,以错峰和避峰生产方式组织生产,可以节约生产用电成本并促进用电负荷削峰。
值得注意的是,轧制速度对轧机负荷影响很大,在考虑分时电价的情况下,轧机负荷对用电成本的影响同样很大。在保证生产考虑轧制生产节奏的前提下,对轧制速度在一定的裕量范围内进行调节,可以进一步节约生产用电成本并促进用电负荷削锋。分时电价环境下,如何考虑热轧批量计划编制、轧制待产时间、轧制速度调节之间的作用关系,并通过一体化的热轧优化调度降低生产用电成本还缺乏有效的方法,有待研究和实现。
发明内容
为了解决上述技术问题,本发明提供一种通过一体化的批量计划编制、轧制单元待产时间分配和轧制速度调节,以错峰和避峰生产方式组织生产的降低生产用电成本的热轧优化调度方法。
本发明解决上述问题的技术方案是:一种降低生产用电成本的热轧优化调度方法,包括以下步骤:
步骤一:加载分时电价费率数据、轧制速度设定数据和待排产的板坯数据;
步骤二:确定最小化生产用电成本的热轧优化调度的目标函数和约束条件;
步骤三:批量计划编制、轧制单元待产时间和轧制速度信息的混合染色体编码;
步骤四:确定混合染色体编码的选择、交叉、变异算子及其操作;
步骤五:执行非支配排序多目标遗传算法优化热轧调度目标函数;
步骤六:解集的多目标决策及推荐方案输出。
上述降低生产用电成本的热轧优化调度方法,所述步骤一中,分时电价费率数据表示为πj,即电价时段j的电价,其中j∈T,T为电价时段集合,T={1,2,…,t},t为一个电价费率周期划分的电价时段个数;轧制速度设定数据表示为V'={v1',v'2,…,vt'},v'j为电价时段j设置的轧制速度。
上述降低生产用电成本的热轧优化调度方法,所述步骤二中,所述目标函数为:
①最小化相邻板坯间宽度、厚度、硬度跳变惩罚值的目标函数为:
②最小化生产用电成本的目标函数为:
所述约束条件为:
v'min≤vi≤v'max,i∈M;
所述约束条件和目标函数中,N为板坯库中的板坯序号集合,N={1,2,...,n},n为轧制板坯数;M为轧制单元序号集合,M={1,2,...,m},m为轧制单元数;Cij为相邻板坯因宽度、厚度、硬度跳变引起的惩罚值,其中和分别为相邻板坯i和j的宽度、厚度、硬度跳变惩罚值;vi为轧制单元i的平均轧制速度;li为第i个板坯的轧制长度;v'min和v'max分别为V’中的最小值和最大值;L为轧制单元内板坯连续轧制长度下限值;U为轧制单元内板坯连续轧制长度上限值;R为轧制单元内相同宽度板坯连续轧制长度上限值;bi表示轧制单元i的待产时间;Ts为可安排生产时间;pi为轧制单元i生产所需最大时间,即平均轧制速度为v'min时所需的生产时间;Wi为板坯i轧制生产所需电能,i∈N;
多目标优化的决策变量为:
上述降低生产用电成本的热轧优化调度方法,所述步骤三中,混合染色体编码表示为G由自然数序列C、浮点数序列B和整数序列V构成:G=(C,B,V)=(c1,c2,...,cm×n,b1,b2,...,bm,v1,v2,...,vm);其中,C中元素ci(i=1,2,…,m×n)为[1,m×n]范围内不重复的自然数,通过编码映射由序列C确定批量计划编制方案,m为轧制单元数,n为轧制板坯数;B中元素bk(1≤k≤m)表示第k个轧制单元生产前的待产时间;V中元素vk(1≤k≤m)表示第k个轧制单元的平均轧制速度。
上述降低生产用电成本的热轧优化调度方法,所述步骤三中,编码映射的具体步骤为:
1)令板坯i的利用标志变量fi=0(i=1,2,...,n);对于轧制单元k,令板坯数numk=0(k=1,2,...,m),板坯轧制总长度dk=0(k=1,2,...,m),相同宽度板坯轧制总长度qk=0(k=1,2,...,m);设置循环变量j=1;
2)确定自然数cj对应的板坯s与轧制单元k,计算方法为:
其中[]表示向下取整;
3)判断fs=0?
i.若fs=0成立,表明板坯s为自由板坯,判断ws≠w'k是否成立,ws≠w'k成立则设置qk=0,其中ws表示板坯s的宽度,w'k表示轧制单元k中最近加入板坯的宽度,ws≠w'k不成立则进一步判断dk+ls≤U且qk+ls≤R是否成立,不成立则转步骤4),成立则将板坯s加入轧制计划k,令numk=numk+1,dk=dk+ls,qk=qk+ls,fs=1,轧制单元k中板坯序列Rk=Rk∪{s};
ii.如果fs=0不成立,直接转步骤4);
4)令j=j+1,转步骤2)重复执行上述步骤,直到j=m×n+1;
5)检查是否所有fi=1(i=1,2,...,n)且dk≥L(k=1,2,...,m)成立;
i.若成立,表明在满足约束条件前提下所有板坯均分配到轧制单元中,集合R={R1,R2,…,Rm}即为批量计划编制的一个可行解;
ii.若不成立,表明该染色体编码表示的批量计划不满足约束,为不可行解,对该编码赋予一个适应度函数值,使其不被选入新种群。
上述降低生产用电成本的热轧优化调度方法,所述步骤四中,选择算子的操作为:采用基于非支配排序和拥挤距离计算的锦标赛排序方法选择种群内个体;交叉算子的操作为:采用部分交叉匹配策略对父代染色体编码的序列C进行交叉重组,交叉重组后在序列C的基础上规划序列B和V;变异算子的操作为:采用子串重组变异策略对父代染色体编码的序列C进行变异操作,变异操作后在序列C的基础上规划序列B和V。
上述降低生产用电成本的热轧优化调度方法,所述步骤四中,在序列C的基础上规划序列B和V的具体步骤为:
3)将t个电价时段按电价降序排序,得到排序后的电价时段集合T'=(t'1,t'2,...,t't),t'k对应的电价为π'k(k=1,2,...,t),与该电价时段相对应的轧制速度为令循环变量z=1;
4)初始化每个轧制单元的平均轧制速度vi=v'min,i∈M;
5)调整轧制单元的平均轧制速度以及待产时间,对从电价时段t'k开始生产的轧制单元i,若满足且bi+1>0,i∈M,且t'k+1≤t,则令 bi=bi+1+bi+Δti,bi+1=0,其中Δti表示轧制单元i提高平均轧制速度后缩短的轧制时间,i∈M;对从电价时段t'k结束生产的轧制单元i,若满足且bi>0,且t'k-1≥1,则令bi+1=bi+1+bi+Δti,bi=0;
上述降低生产用电成本的热轧优化调度方法,所述步骤五中,非支配排序多目标遗传算法为NSGA-II算法。
上述降低生产用电成本的热轧优化调度方法,所述步骤六中,解集的多目标决策为采用TOPSIS双基点法对优化解集内的有限个可行解进行排序,以选取与理想解接近程度最大的解作为最优可行解,选取最优可行解作为热轧优化调度的推荐方案输出。
上述降低生产用电成本的热轧优化调度方法,所述步骤六中,对优化解集内的有限个可行解进行排序的具体步骤如下:
步骤1:初始决策矩阵X可以表示为:X是一个大小为ROW×COL的二维矩阵,其中ROW为多目标算法的可行解个数,COL为优化目标个数,这里COL=2,X中的元素xrow,col为第row个可行解第col个目标的值,用向量规范法求得规范决策矩阵Z,Z中元素zrow,col的计算方法为:
步骤2:构成加权规范阵Q,Q中元素qrow,col=wcol·zrow,col,其中wcol为第col个目标的权重;
步骤3:确定理想解和负理想解;
步骤4:计算各个可行解到正理想解与负理想解的距离;
本发明的有益效果在于:本发明在批量计划编制时同步构造与电价时段相适应的生产负荷单元,在轧制生产能力裕量下通过合理安排轧制单元待产时间以及调整板坯的轧制速度进行生产负荷的时间转移,以错峰和避峰生产方式组织生产,基于多目标遗传算法优化相邻板坯间惩罚值和生产用电成本,可在保证生产前提下节约生产用电成本并促进用电负荷削峰。
附图说明
图1为本发明的流程图。
图2为本发明非支配排序的多目标遗传算法流程图。
图3为对从电价时段t'k开始生产的轧制单元i,若开始生产时的电价大于结束生产时的电价时,其平均轧制速度以及待产时间的调整示意图。
图4为对从电价时段t'k结束生产的轧制单元i,若开始生产时的电价小于结束生产时的电价时,其平均轧制速度以及待产时间的调整示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示,本发明一种降低生产用电成本的热轧优化调度方法的具体实现步骤如下:
步骤S101:加载分时电价费率数据、轧制速度设定数据和待排产的板坯数据。
加载表1所示的某钢铁厂执行的分时电价费率,该费率表由供电企业根据地区电网负荷将每日划分为尖、高、平、谷4类共计8个电价时段。
表1分时电价费率表
加载每个电价段轧制速度的设定数据,轧制速度设定数据表示为V'={v1',v'2,…,vt'},v'j为电价时段j设置的轧制速度。
加载从某钢铁厂生产数据中选取的一组板坯数据。
步骤S102、确定考虑电价和轧制速度调节的热轧优化调度的目标函数和约束条件。
所述目标函数为:
①最小化相邻板坯间宽度、厚度、硬度跳变惩罚值的目标函数为:
②最小化生产用电成本的目标函数为:
所述约束条件为:
v'min≤vi≤v'max,i∈M
约束条件和目标函数中,N为板坯库中的板坯序号集合,N={1,2,...,n},n为板坯数,M为轧制单元序号集合,M={1,2,...,m},m为轧制单元数;Cij为相邻板坯因宽度、厚度、硬度跳变引起的惩罚值,其中和分别为相邻板坯i和j的宽度、厚度、硬度跳变惩罚值,相邻板坯间宽度、厚度和硬度跳变惩罚值按表2~表4所示惩罚系数表确定:
表2宽度跳变惩罚系数表
表3厚度跳变惩罚系数表
表4硬度跳变惩罚系数表
vi为轧制单元i的平均轧制速度;v'min和v'max分别为V’中的最小值和最大值;li为第i个板坯的轧制长度,i∈N;L为轧制单元内板坯连续轧制长度下限值,取值为5km;U为轧制单元内板坯连续轧制长度上限值,取值为10km;R为轧制单元内相同宽度板坯连续轧制长度上限值,取值为1km;bi表示轧制单元i的待产时间;Ts为可安排生产时间;pi为轧制单元i生产所需最大时间,即平均轧制速度为v'min时所需的生产时间;T为电价时段集合,T={1,2,…,t},t为一个电价费率周期划分的电价时段个数,一个费率周期取值24h;πj为电价时段j的电价,j∈T;Wi为板坯i轧制生产所需电能,i∈N;单个板坯的轧制长度和轧制生产所需电能、所需时间通过钢铁厂热轧过程控制系统预测取得;
多目标优化的决策变量为:
步骤S103、批量计划编制、轧制单元待产时间和轧制速度信息的混合染色体编码。
混合染色体编码表示为G由自然数序列C、浮点数序列B和整数序列V构成:G=(C,B,V)=(c1,c2,...,cm×n,b1,b2,...,bm,,v1,v2,...,vm);
其中,C中元素ci(i=1,2,…,m×n)为[1,m×n]范围内不重复的自然数,通过编码映射由序列C确定批量计划编制方案,m为轧制单元数,n为轧制板坯数;B中元素bk(1≤k≤m)表示第k个轧制单元生产前的待产时间;V中元素vk(1≤k≤m)表示第k个轧制单元的平均轧制速度。
编码映射的具体步骤为:
1)令板坯i的利用标志变量fi=0(i=1,2,…,n);对于轧制单元k,令板坯数numk=0(k=1,2,…,m),板坯轧制总长度dk=0(k=1,2,…,m),相同宽度板坯轧制总长度qk=0(k=1,2,…,m);设置循环变量j=1;
2)确定自然数cj对应的板坯s与轧制单元k,计算方法为:
其中[]表示向下取整;
3)判断fs=0?
i.若fs=0成立,表明板坯s为自由板坯,判断ws≠w'k是否成立,ws≠w'k成立则设置qk=0,其中ws表示板坯s的宽度,w'k表示轧制单元k中最近加入板坯的宽度,ws≠w'k不成立则进一步判断dk+ls≤U且qk+ls≤R是否成立,不成立则转步骤4),成立则将板坯s加入轧制计划k,令numk=numk+1,dk=dk+ls,qk=qk+ls,fs=1,轧制单元k中板坯序列Rk=Rk∪{s};
ii.如果fs=0不成立,直接转步骤4);
4)令j=j+1,转步骤2)重复执行上述步骤,直到j=m×n+1;
5)检查是否所有fi=1(i=1,2,...,n)且dk≥L(k=1,2,...,m)成立;
i.若成立,表明在满足约束条件前提下所有板坯均分配到轧制单元中,集合R={R1,R2,…,Rm}即为批量计划编制的一个可行解;
ii.若不成立,表明该染色体编码表示的批量计划不满足约束,为不可行解,对该编码赋予一个较大适应度函数值,使其不被选入新种群。
步骤S104、确定混合染色体编码的选择、交叉、变异算子及其操作。
所述选择算子的操作为:采用基于非支配排序和拥挤距离计算的锦标赛排序方法选择种群内个体;
所述交叉算子的操作为:采用部分交叉匹配策略对父代染色体编码的序列C进行交叉重组,交叉重组后在序列C的基础上规划序列B和V;
所述变异算子的操作为:采用子串重组变异策略对父代染色体编码的序列C进行变异操作,变异操作后在序列C的基础上规划序列B和V。
所述在序列C的基础上规划序列B和V的具体步骤为:
3)将t个电价时段按电价降序排序,得到排序后的电价时段集合T'=(t'1,t'2,...,t't),t'k对应的电价为π'k(k=1,2,...,t),与该电价时段相对应的轧制速度为令循环变量z=1;
4)初始化每个轧制单元的平均轧制速度vi=v'min,i∈M;
5)调整轧制单元的平均轧制速度以及待产时间,图3所示为对从电价时段t'k开始生产的轧制单元i,若开始生产时的电价大于结束生产时的电价时,其平均轧制速度以及待产时间的调整过程,S301表示构成轧制单元的板坯,方框中序号表示其在轧制单元内的生产顺序,对从电价时段t'k开始生产的轧制单元i,若满足且bi+1>0,i∈M,且t'k+1≤t,则令bi=bi+1+bi+Δti,bi+1=0,Δti表示轧制单元i提高平均轧制速度后缩短的轧制时间,i∈M;
图4所示为对从电价时段t'k结束生产的轧制单元i,若开始生产时的电价小于结束生产时的电价时,其平均轧制速度以及待产时间的调整过程。S401表示构成轧制单元的板坯,方框中序号表示其在轧制单元内的生产顺序,对从电价时段t'k结束生产的轧制单元i,若满足且bi>0,且t'k-1≥1,则令 bi+1=bi+1+bi+Δti,bi=0;
步骤S105、执行非支配排序多目标遗传算法优化热轧调度目标函数,该非支配排序多目标遗传算法为NSGA-II算法。
步骤S106、解集的多目标决策及推荐方案输出:解集的多目标决策为采用TOPSIS双基点法对优化解集内的有限个可行解进行排序,以选取与理想解接近程度最大的解作为最优可行解,选取最优可行解作为热轧优化调度的推荐方案输出。
对优化解集内的有限个可行解进行排序的具体步骤如下具体步骤如下:
步骤1:初始决策矩阵X可以表示为:X是一个大小为ROW×COL的二维矩阵,其中ROW为多目标算法的可行解个数,COL为优化目标个数,这里COL=2,X中的元素xrow,col为第row个可行解第col个目标的值,用向量规范法求得规范决策矩阵Z,Z中元素zrow,col的计算方法为:
步骤2:构成加权规范阵Q,Q中元素qrow,col=wcol·zrow,col,其中wcol为第col个目标的权重;
步骤3:确定理想解和负理想解;
步骤4:计算各个可行解到正理想解与负理想解的距离;
选取根据TOPSIS方法得到的最优可行解作为热轧优化调度的推荐方案输出。
在分时电价环境下,该方法通过轧制批量计划编制时同步调节各个轧制单元的待产时间和轧制速度,实现了错峰生产,进而降低用电成本。除节约生产用电费用外,该方法的意义还在于通过响应电价信号构造了与电价时段相适应的生产负荷单元,并通过生产负荷的时间转移,参与了电网“移峰填谷”的运行调峰,可促进电网稳定。
Claims (4)
1.一种降低生产用电成本的热轧优化调度方法,包括以下步骤:
步骤一:加载分时电价费率数据、轧制速度设定数据和待排产的板坯数据;
步骤一中,分时电价费率数据表示为πj,即电价时段j的电价,其中j∈T,T为电价时段集合,T={1,2,...,t},t为一个电价费率周期划分的电价时段个数;轧制速度设定数据表示为V'={v′1,v′2,...,v′t},v′j为电价时段j设置的轧制速度;
步骤二:确定最小化生产用电成本的热轧调度目标函数和约束条件;所述目标函数为:
①最小化相邻板坯间宽度、厚度、硬度跳变惩罚值的目标函数为:
②最小化生产用电成本的目标函数为:
所述约束条件为:
v'min≤vi≤v'max,i∈M;
所述约束条件和目标函数中,N为板坯库中的板坯序号集合,N={1,2,...,n},n为轧制板坯数;M为轧制单元序号集合,M={1,2,...,m},m为轧制单元数;Cij为相邻板坯因宽度、厚度、硬度跳变引起的惩罚值,其中和分别为相邻板坯i和j的宽度、厚度、硬度跳变惩罚值;vi为轧制单元i的平均轧制速度;li为第i个板坯的轧制长度;v'min和v'max分别为V’中的最小值和最大值;L为轧制单元内板坯连续轧制长度下限值;U为轧制单元内板坯连续轧制长度上限值;R为轧制单元内相同宽度板坯连续轧制长度上限值;bi表示轧制单元i的待产时间;Ts为可安排生产时间;pi为轧制单元i生产所需最大时间,即平均轧制速度为v'min时所需的生产时间;Wi为板坯i轧制生产所需电能,i∈N;
多目标优化的决策变量为:
步骤三:批量计划编制、轧制单元待产时间和轧制速度信息的混合染色体编码;
所述步骤三中,混合染色体编码表示为G由自然数序列C、浮点数序列B和整数序列V构成:G=(C,B,V)=(c1,c2,...,cm×n,b1,b2,...,bm,v1,v2,...,vm);其中,C中元素ci为[1,m×n]范围内不重复的自然数,i=1,2,…,m×n,通过编码映射由序列C确定批量计划编制方案,m为轧制单元数,n为轧制板坯数;B中元素bk表示第k个轧制单元生产前的待产时间,1≤k≤m;V中元素vk表示第k个轧制单元的平均轧制速度,1≤k≤m;
编码映射的具体步骤为:
1)令板坯i的利用标志变量fi=0,i=1,2,...,n;对于轧制单元k,令板坯数numk=0,k=1,2,...,m,板坯轧制总长度dk=0,k=1,2,...,m,相同宽度板坯轧制总长度qk=0,k=1,2,...,m;设置循环变量j=1;
2)确定自然数cj对应的板坯s与轧制单元k,计算方法为:
其中[]表示向下取整;
3)判断fs=0是否成立;
i.若fs=0成立,表明板坯s为自由板坯,判断ws≠w'k是否成立,ws≠w'k成立则设置qk=0,其中ws表示板坯s的宽度,w'k表示轧制单元k中最近加入板坯的宽度,ws≠w'k不成立则进一步判断dk+ls≤U且qk+ls≤R是否成立,不成立则转步骤4),成立则将板坯s加入轧制单元k,令numk=numk+1,dk=dk+ls,qk=qk+ls,fs=1,轧制单元k中板坯序列Rk=Rk∪{s};
ii.如果fs=0不成立,直接转步骤4);
4)令j=j+1,转步骤2)重复执行上述步骤,直到j=m×n+1;
5)检查是否所有fi=1,i=1,2,...,n且dk≥L,k=1,2,...,m成立;
i.若成立,表明在满足约束条件前提下所有板坯均分配到轧制单元中,集合R={R1,R2,...,Rm}即为批量计划编制的一个可行解;
ii.若不成立,表明该染色体编码表示的批量计划不满足约束,为不可行解,对该编码赋予一个适应度函数值,使其不被选入新种群;
步骤四:确定混合染色体编码的选择、交叉、变异算子及其操作;
选择算子的操作为:采用基于非支配排序和拥挤距离计算的锦标赛排序方法选择种群内个体;交叉算子的操作为:采用部分交叉匹配策略对父代的混合染色体编码的序列C进行交叉重组,交叉重组后在序列C的基础上规划序列B和V;变异算子的操作为:采用子串重组变异策略对父代的混合染色体编码的序列C进行变异操作,变异操作后在序列C的基础上规划序列B和V;
在序列C的基础上规划序列B和V的具体步骤为:
4)初始化每个轧制单元的平均轧制速度vi=v'min,i∈M;
5)调整轧制单元的平均轧制速度以及待产时间,对从电价时段t'k开始生产的轧制单元i,若满足且bi+1>0,i∈M,且t'k+1≤t,则令 bi=bi+1+bi+Δti,bi+1=0,其中Δti表示轧制单元i提高平均轧制速度后缩短的轧制时间,i∈M;对从电价时段t'k结束生产的轧制单元i,若满足且bi>0,且t'k-1≥1,则令bi+1=bi+1+bi+Δti,bi=0;
步骤五:执行非支配排序多目标遗传算法优化热轧调度目标函数;
步骤六:优化解集的多目标决策及推荐方案输出。
2.根据权利要求1所述的降低生产用电成本的热轧优化调度方法,其特征在于:所述步骤五中非支配排序多目标遗传算法为NSGA-II算法。
3.根据权利要求2所述的降低生产用电成本的热轧优化调度方法,其特征在于:所述步骤六中,优化解集的多目标决策为采用TOPSIS双基点法对优化解集内的有限个可行解进行排序,以选取与理想解接近程度最大的解作为最优可行解,选取最优可行解作为热轧优化调度的推荐方案输出。
4.根据权利要求3所述的降低生产用电成本的热轧优化调度方法,其特征在于:所述步骤六中,对优化解集内的有限个可行解进行排序的具体步骤如下:
步骤1:初始决策矩阵X可以表示为:X是一个大小为ROW×COL的二维矩阵,其中ROW为非支配排序多目标遗传算法的可行解个数,COL为优化目标个数,这里COL=2,X中的元素xrow,col为第row个可行解第col个目标的值,用向量规范法求得规范决策矩阵Z,Z中元素zrow,col的计算方法为:
步骤2:构成加权规范阵Q,Q中元素qrow,col=wcol·zrow,col,其中wcol为第col个目标的权重;
步骤3:确定正理想解和负理想解;
步骤4:计算各个可行解到正理想解与负理想解的距离;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710198430.2A CN106971243B (zh) | 2017-03-29 | 2017-03-29 | 一种降低生产用电成本的热轧优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710198430.2A CN106971243B (zh) | 2017-03-29 | 2017-03-29 | 一种降低生产用电成本的热轧优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106971243A CN106971243A (zh) | 2017-07-21 |
CN106971243B true CN106971243B (zh) | 2020-06-02 |
Family
ID=59335515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710198430.2A Active CN106971243B (zh) | 2017-03-29 | 2017-03-29 | 一种降低生产用电成本的热轧优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106971243B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110648050B (zh) * | 2019-08-21 | 2023-05-02 | 大连理工大学 | 一种传统流水线装配向单元式装配方式转换的重构方法 |
CN112462699A (zh) * | 2020-11-02 | 2021-03-09 | 湖南华菱涟源钢铁有限公司 | 热轧生产电耗管理系统及方法 |
CN115178590B (zh) * | 2022-09-07 | 2022-12-06 | 承德建龙特殊钢有限公司 | 一种无缝钢管生产方法 |
CN116882593B (zh) * | 2023-09-05 | 2024-02-23 | 美云智数科技有限公司 | 工单排产方法、装置、电子设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593719A (zh) * | 2013-11-29 | 2014-02-19 | 湘潭大学 | 一种基于板坯与合同优化匹配的轧制节能方法 |
CN104209354A (zh) * | 2014-07-17 | 2014-12-17 | 镇江金昶坤能源设备有限公司 | 轧钢除鳞系统的喷水泵节能控制方法 |
CN104463414A (zh) * | 2014-10-23 | 2015-03-25 | 湘潭大学 | 一种分时电价下热轧生产经济负荷调度方法 |
CN104599055A (zh) * | 2015-01-04 | 2015-05-06 | 湘潭大学 | 一种钢铁企业的电力负荷需量管理方法 |
-
2017
- 2017-03-29 CN CN201710198430.2A patent/CN106971243B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103593719A (zh) * | 2013-11-29 | 2014-02-19 | 湘潭大学 | 一种基于板坯与合同优化匹配的轧制节能方法 |
CN104209354A (zh) * | 2014-07-17 | 2014-12-17 | 镇江金昶坤能源设备有限公司 | 轧钢除鳞系统的喷水泵节能控制方法 |
CN104463414A (zh) * | 2014-10-23 | 2015-03-25 | 湘潭大学 | 一种分时电价下热轧生产经济负荷调度方法 |
CN104599055A (zh) * | 2015-01-04 | 2015-05-06 | 湘潭大学 | 一种钢铁企业的电力负荷需量管理方法 |
Non-Patent Citations (1)
Title |
---|
钢铁企业电力合理生产与优化研究;刘向斌;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20140715;第1-79页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106971243A (zh) | 2017-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106971243B (zh) | 一种降低生产用电成本的热轧优化调度方法 | |
CN110112767B (zh) | 广域多形态需求侧负荷参与系统调峰的荷源优化控制方法 | |
CN105956689B (zh) | 一种基于改进粒子群优化的运输和生产协同调度方法 | |
CN104808636B (zh) | 柔性流水车间能耗优化调度方法 | |
CN104463414B (zh) | 一种分时电价下热轧生产经济负荷调度方法 | |
CN101604409B (zh) | 应用于多重入复杂制造系统的调度方法 | |
CN110598941A (zh) | 一种基于仿生策略的粒子群优化制造系统双目标排产方法 | |
CN105652791A (zh) | 订单驱动的离散制造过程能耗优化方法 | |
CN108710970B (zh) | 一种巨型梯级水电系统的多目标调度并行降维方法 | |
CN110414863A (zh) | 一种智能制造车间资源调度方法 | |
CN110247428A (zh) | 一种基于源网荷协同优化的配电网光储联合配置方法 | |
CN107766980A (zh) | 基于用户行为习惯的电动汽车分时充电电价优化方法 | |
CN115147245B (zh) | 一种工业负荷参与调峰辅助服务的虚拟电厂优化调度方法 | |
CN104182804A (zh) | 一种预测出力不确定小水电与大中型水电站协调的日前发电方法 | |
CN103617470A (zh) | 一种设备能力变化条件下的选矿综合生产指标优化方法 | |
CN112952847A (zh) | 考虑用电需求弹性的多区域主动配电系统调峰优化方法 | |
CN101989743B (zh) | 基于直流潮流的节能发电调度优化方法 | |
CN111783254A (zh) | 一种基于多目标混合整数规划的钢材切割控制方法及装置 | |
CN104616071B (zh) | 一种风光储互补发电系统优化配置方法 | |
CN113537568A (zh) | 一种钢铁热轧生产调度计划确定方法及系统 | |
CN111626492A (zh) | 一种电动汽车充电网络模糊多目标机会约束规划方法 | |
CN114597907A (zh) | 一种基于水泥企业不连续生产过程及仓储能力的电力需求响应方法 | |
CN102206741A (zh) | 一种罩式退火炉钢卷堆垛组合优化方法 | |
CN112148446B (zh) | 一种用于多技能资源受限项目调度的进化策略方法 | |
CN111160384B (zh) | 一种移动储能车调度优化方法、系统及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |