CN112669912A - 考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法 - Google Patents
考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法 Download PDFInfo
- Publication number
- CN112669912A CN112669912A CN202011340916.3A CN202011340916A CN112669912A CN 112669912 A CN112669912 A CN 112669912A CN 202011340916 A CN202011340916 A CN 202011340916A CN 112669912 A CN112669912 A CN 112669912A
- Authority
- CN
- China
- Prior art keywords
- cracking
- batch
- ethylene
- time
- raw material
- 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.)
- Pending
Links
- 238000005336 cracking Methods 0.000 title claims abstract description 123
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 title claims abstract description 70
- 239000005977 Ethylene Substances 0.000 title claims abstract description 69
- 239000002994 raw material Substances 0.000 title claims abstract description 63
- 238000004939 coking Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004519 manufacturing process Methods 0.000 claims abstract description 26
- 238000005235 decoking Methods 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 230000008859 change Effects 0.000 claims abstract description 10
- 239000000571 coke Substances 0.000 claims description 17
- 239000000463 material Substances 0.000 claims description 15
- 238000010923 batch production Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 4
- 230000007613 environmental effect Effects 0.000 abstract description 8
- 230000008901 benefit Effects 0.000 abstract description 6
- 230000009467 reduction Effects 0.000 abstract description 5
- 238000004134 energy conservation Methods 0.000 abstract description 3
- 239000003344 environmental pollutant Substances 0.000 abstract description 3
- 231100000719 pollutant Toxicity 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 18
- 239000003915 liquefied petroleum gas Substances 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000004140 cleaning Methods 0.000 description 3
- 238000003912 environmental pollution Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000010485 coping Effects 0.000 description 2
- WSLDOOZREJYCGB-UHFFFAOYSA-N 1,2-Dichloroethane Chemical compound ClCCCl WSLDOOZREJYCGB-UHFFFAOYSA-N 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- IAYPIBMASNFSPL-UHFFFAOYSA-N Ethylene oxide Chemical compound C1CO1 IAYPIBMASNFSPL-UHFFFAOYSA-N 0.000 description 1
- XTXRWKRVRITETP-UHFFFAOYSA-N Vinyl acetate Chemical compound CC(=O)OC=C XTXRWKRVRITETP-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 229920001903 high density polyethylene Polymers 0.000 description 1
- 239000004700 high-density polyethylene Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 229920001684 low density polyethylene Polymers 0.000 description 1
- 239000004702 low-density polyethylene Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 239000003348 petrochemical agent Substances 0.000 description 1
- 238000004227 thermal cracking Methods 0.000 description 1
Images
Landscapes
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明公开了一种考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法,包括:获取不同裂解原料投放进不同炉型裂解炉中的平均结焦量与乙烯产品量,获得调度模型的目标函数,建立约束条件,构造MINLP模型,使用DICOPT求解器对MINLP模型进行优化求解。本发明综合考虑乙烯生产中产率随时间衰减的特性,兼顾裂解原料负荷的变化,为每个裂解炉规划调度时间范围内的批次、批处理时间、清焦顺序以及批次原料进料量的最佳安排。此外,本发明可以在清焦阶段减少污染物的排放,在牺牲少量利润的代价下,获取更加可观的环境效益,为乙烯工厂的节能减排、优化排产提供了理论依据。
Description
技术领域
本发明涉及乙烯生产技术领域,尤其涉及一种考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法。
背景技术
乙烯是石化/化工行业最重要的有机化合物,其产量大大超过其他石化产品。广泛的乙烯衍生物在人们的日常生活中起着极其重要的作用,如环氧乙烷,乙酸乙烯酯,二氯化乙烯和高/低密度聚乙烯等。蒸汽热裂解是乙烯装置的第一个操作部分,在很大程度上决定了下游产品的收率和整个装置的能耗。如今,随着利润率越来越紧缩和原料市场波动,鼓励乙烯工厂多样化原料组合,并在面对原材料不确定性的情况下提高其运营弹性。
随着目前我国对于乙烯工业生产中节能减排、低碳环保日益严格的背景下,如何有效减少乙烯工厂的污染物排放是需要特别重视,而减少裂解炉结焦量可以从根源上达到环保要求。因此,对于裂解炉炉群***的改进与优化已经提上日程。由于先前的研究大多集中在单台裂解炉的操作优化上,实际工业中是由多台裂解炉并行生产乙烯产品。因此,需要对裂解炉炉群生产进行原料调度,从而达到相应的优化目标。
发明内容
为解决现有技术存在的局限和缺陷,本发明提供一种考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法,包括:
获取不同裂解原料投放进不同炉型裂解炉中的平均结焦量与乙烯产品量;
根据所述平均结焦量与所述乙烯产品量,获得调度模型的目标函数,所述目标函数用于最小化单位乙烯产量的平均结焦量,所述目标函数的表达式如下:
建立约束条件,所述约束条件包括原料平衡约束、时间约束、整数约束、非同步清焦约束、变量上下限约束和附加约束;
考虑裂解原料负荷变化,通过模拟获得不同负荷水平下乙烯收率模型,所述裂解原料负荷仅在操作条件发生变化时在进料范围内进行调整优化,在单个调度周期内保持不变;
根据所述目标函数、所述约束条件和所述乙烯收率模型,形成炉群调度的MINLP模型;
使用DICOPT求解器对所述MINLP模型进行优化求解,使用子求解器CPLEX和子求解器CONOPT分别处理MILP问题和NLP问题。
可选的,所述乙烯收率模型是关于时间与负荷的函数关系,所述乙烯收率模型的表达式如下:
ai,j(Di,j)=pmai,j×Fi,j+pnai,j (2)
bi,j(Di,j)=pmbi,j×Fi,j+pnbi,j (3)
cij(Dij)=pmcij×Fij+pncij (4)
可选的,所述原料平衡约束的表达式如下:
Gi≤(Dupi-Dloi)H (9)
Floi,jgyi,j,k≤Fi,j,k≤Fupi,jgyi,j,k (10)
所有裂解炉裂解的总进料量在原料供应的范围内,调度时间范围H内的进料i的消耗量,Gi代表所述进料i超出下限的过量,所述裂解原料的进料量在预设范围内进行调整,yi,j,k为二元逻辑变量,若所述进料i未分配给第j炉的第k个批次生产运行,所述裂解原料的进料量Fi,j,k为0。
可选的,所述时间约束的表达式如下:
所述进料i的连续处理时间设置在预设的最小值和最大值之内,若所述进料i未分配给第j炉的第k个批次生产运行,所述处理时间为0;
Tsj,k为所有裂解炉的裂解批次的开始时间,Tej,k为所有裂解炉的裂解批次的结束时间,当前批次的裂解开始时间点等于上一批次的裂解结束时间点加上两批次之间的停炉清焦时间,当前批次的结束时间点等于当前批次的开始时间点加上批处理时间,所有裂解炉第一批次的开始时间均为0,若批处理启动时间点大于调度时间范围H,开始时间Tsj,k等于0;
引入二进制变量yi,j,k表示所述进料i是否分配给第j炉的第k个批次生产运行,若所述裂解原料已经分配,二进制变量yi,j,k为1,若所述裂解原料尚未分配,二进制变量yi,j,k为0,若第k个批次没有生产运行,所述第k个批次的开始和结束时间相同。
可选的,所述整数约束的表达式如下:
所有裂解炉的第一个批次始终用来裂解进料,每个批次最多只能裂解一种进料,在整个调度时间范围内,所有类型的原料至少处理一次,整个调度过程依次先后使用前一批次和后一批次,若裂解炉的某个批次未使用,所述裂解炉的后续批次都不会使用。
可选的,所述非同步清焦约束的表达式如下:
不同裂解炉清焦的时间间隔禁止具有重叠部分,使用两个起点之间的时间差和两个终点之间的时间差对所述非同步清焦约束进行描述,所述两个起点的时间差和所述两个终点的时间差之间的乘积小于或等于0。
可选的,所述变量上下限约束的表达式如下:
Tpi,j,k≥0,Gi≥0,yi,j,k∈{0,1},H>0 (22)
所述开始时间Tsj,k的上限为0,所述结束时间Tej,k的下限H,所述调度时间范围H大于0,所述二元逻辑变量yi,j,k为0或1,其他变量不小于0。
可选的,所述附加约束的表达式如下:
所述附加约束为每种调度方案的日平均利润。
本发明具有下述有益效果:
本发明考虑裂解炉实际生产过程,摒弃了原有模型中某种原料进入某裂解炉所有批次时间相等的假设,引入二元变量来表示原料是否进入裂解炉的第k个批次进行裂解,更加符合实际生产情况。同时,考虑裂解炉负荷变化的情况能够更好地优化每种原料对于各台裂解炉的进料负荷,以及每台裂解炉所投放原料的组合。本发明还充分考虑乙烯生产过程中非同步清焦的生产特性,从而保证调度的正常稳健进行。本发明可以有效地模拟乙烯化工炉群调度生产,有效提升环境效益,从而有效应对目前愈加严重的环境污染,可以对实际生产起到参考与指导作用。
附图说明
图1为本发明实施例一提供的裂解炉炉群调度模型的调度示意图。
图2a为本发明实施例一提供的乙烯裂解产率随时间变化的一种示意图。
图2b为本发明实施例一提供的乙烯裂解产率随时间变化的另一种示意图。
图3a为本发明实施例一提供的结焦量随运行时间的变化的一种趋势图。
图3b为本发明实施例一提供的结焦量随运行时间的变化的另一种趋势图。
图3c为本发明实施例一提供的结焦量随运行时间的变化的又一种趋势图。
图4a为本发明实施例一提供的基础模型调度方案的一种示意图。
图4b为本发明实施例一提供的基础模型调度方案的另一种示意图。
图4c为本发明实施例一提供的基础模型调度方案的又一种示意图。
图5a为本发明实施例一提供的改进模型调度方案的一种示意图。
图5b为本发明实施例一提供的改进模型调度方案的另一种示意图。
图5c为本发明实施例一提供的改进模型调度方案的又一种示意图。
具体实施方式
为使本领域的技术人员更好地理解本发明的技术方案,下面结合附图对本发明提供的考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法进行详细描述。
实施例一
本实施例的目的在于:以单位乙烯产品的平均结焦量为目标函数,兼顾考虑原料负荷变化,构建改进的乙烯裂解炉炉群调度模型,求解后给出可靠、高效的炉群调度解决方案,为乙烯企业提高生产效率、降低环境污染提供技术支撑与技术参考。本实施例确立以单位乙烯产品的平均结焦量最小化为目标函数,通过乙烯实际生产过程建立相应的约束条件、考虑原料负荷变化、对炉群调度问题进行建模、对所建模型进行优化求解。
本实施例提供的目标函数需要引入平均结焦量的概念。本实施例定义平均结焦量avecoke为裂解炉每生产1吨乙烯,炉管内壁所生成的焦炭质量。这里使用变量output来表示在调度时间范围内裂解炉***中乙烯产品的总产量,变量cokemass表示在调度时间范围内裂解炉***中生成的总焦炭质量。通过COILSIM1D所仿真得到的裂解炉炉管结焦量随时间变化的拟合函数关系,其中P1、P2、P3、P4都是拟合参数。本实施例对裂解炉炉群***中所有裂解炉,所有批次的乙烯产量与炉管结焦量求和,从而计算炉群***总的平均结焦量,进一步得到平均结焦量。
本实施例提供的约束条件包括原料平衡约束、时间约束、整数约束、非同步清焦约束、变量上下限约束以及附加约束。本实施例考虑原料负荷变化,裂解原料的进料量(原料负荷)可根据实际工况在一定范围内进行调整。通过模拟得到不同负荷水平下乙烯收率模型,原料负荷仅在操作条件变化下才在进料范围内进行调整优化,在一个调度周期内保持不变。
本实施例建立调度优化模型,综合上述目标函数和约束变量,建立基本的裂解炉炉群调度的MINLP数学模型。本实施例提供的模型求解使用GAMS求解,求解MINLP问题的求解器则选用DICOPT,并采用子求解器CPLEX和CONOPT来分别处理其中的MILP和NLP问题。本实施例提供一种考虑平均结焦量与负荷变化的乙烯裂解炉炉群调度建模方法,用于优化炉群调度过程,通过改进模型的目标函数与约束条件,可以提高模型的精度,并得到更好的优化结果。本实施例可以实现最优目标值下调度方案的获取,从而起到优化的作用并给出指导方案。
本实施例可以有效应用于并行多进料、多产品的乙烯裂解炉炉群调度***中。图1为本发明实施例一提供的裂解炉炉群调度模型的调度示意图。如图1所示。通过***的最优调度,可以在牺牲少量利润的代价下有效地减少单位产品的平均结焦量。本实施例提出一种新的MINLP模型,以考虑产品产率随时间呈指数衰减的限制下,获取最大裂解炉***环境效益的调度策略。图2a为本发明实施例一提供的乙烯裂解产率随时间变化的一种示意图。图2b为本发明实施例一提供的乙烯裂解产率随时间变化的另一种示意图。乙烯产率与时间的关系如图2a和图2b所示。
本实施例给定以下信息:原料负荷上下限;不同裂解炉裂解不同原料的产品收率模型;各裂解炉裂解不同进料时的清焦时间;批处理时间的上下限;清焦与配置费用;产品价格指数;给定的调度时间范围;每个裂解炉炉管内壁结焦速率随时间变化的拟合函数。可通过最佳调度确定的信息包括:为每个裂解炉分配的批次数量;在每批操作中处理的进料类型;每批次裂解操作的开始时间点Tsj,k和结束时间点Tej,k;整个裂解炉***的具体除焦操作顺序。
本实施例建立目标函数,调度模型的目标函数是最小化每单位乙烯产品的平均结焦量,表达式如下:
本实施例兼顾考虑原料负荷变化,乙烯收率关于时间和负荷的函数关系如下所示:
ai,j(Di,j)=pmai,j×Fi,j+pnai,j (2)
bi,j(Di,j)=pmbi,j×Fi,j+pnbi,j (3)
cij(Dij)=pmcij×Fij+pncij (4)
本实施例建立约束条件,所述约束条件包括原料平衡约束、时间约束、整数约束、非同步清焦约束、变量上下限约束以及附加约束。
对于原料平衡约束,所有裂解炉裂解的总进料量应该在原料供应的范围内。调度时间范围H内的进料i的消耗量,其中Gi代表进料i超出下限的过量。本实施例中,裂解原料的进料量(原料负荷)可根据实际工况在一定范围内进行调整,yi,j,k为二元逻辑变量,如果进料i未分配给第j炉的第k个批次生产运行,则裂解原料的进料量Fi,j,k为0,表达式如下:
Gi≤(Dupi-Dloi)H (9)
Floi,jgyi,j,k≤Fi,j,k≤Fupi,jgyi,j,k (10)
对于时间约束,从管理和操作的角度看,进料i的连续处理时间在实际上允许的最小值和最大值之内。如果进料i未分配给第j炉的第k个批次生产运行,则该批次的处理时间为0。定义了所有裂解炉的裂解批次的开始时间Tsj,k和结束时间Tej,k。当前批次的裂解开始时间点等于上一批次的裂解结束时间点加上两批次之间的停炉清焦时间。当前批次的结束时间点等于其开始时间点加上批处理时间。
假设在时间范围开始时所有裂解装置都是干净的,且所有裂解炉第一批次的开始时间均为0。如果批处理启动时间点大于调度时间范围H,则说明该批处理实际未使用,即Tpj,k等于0。引入二进制变量yi,j,k来表述进料i是否分配给第j炉的第k个批次生产运行。如果原料分配,那么二进制变量yi,j,k为1,反之,则为0。需要注意,一个批次的开始和结束时间是这样安排的:如果第k个批次实际上没有被利用,那么它的开始和结束时间点将是相同的,表达式如下:
对于整数约束,在实际生产中,所有裂解炉的第一个批次必须始终使用来裂解进料。对于所有的裂解炉,每个批次最多只能裂解一种进料。在整个调度时间范围内,所有类型的原料必须至少处理一次。整个调度过程中,使用后一批次之前必须先使用前一批次。若裂解炉的某个批次未使用,则该裂解炉的后续批次都不会使用,表达式如下:
对于非同步清焦约束,不同裂解炉清焦的时间间隔不应该有重叠部分,如图1所示。可以使用两个起点之间的时间差和两个终点之间的时间差来描述该约束,这两类时间差之间的乘积应该小于或等于0。注意,此约束不能应用于每个炉子的最后一次清理,即所有装置的最后生产运行中同时清焦是可行的,表达式如下:
对于变量上下限约束,Tsj,k与Tej,k的上、下限分别为0和H;调度时域H应大于0;二元变量yi,j,k只能取0或1;其他变量则需要不小于0,表达式如下:
Tpi,j,k≥0,Gi≥0,yi,j,k∈{0,1},H>0 (22)
对于附加约束,将每种调度方案下的日平均利润作为该调度模型的附加约束,表达式如下:
通过上述步骤,可以建立以公式(1)为目标函数,公式(2)-(23)为约束条件的新的MINLP优化模型。为了突出本实施例的效果,将未考虑原料负荷变化的旧的基础模型与新模型进行对比分析。
本实施例提供的案例来源于国内某乙烯工厂,研究了3种炉型的裂解炉,分别为GK-VI型、GK-V型以及GK-III型,分别用1、2、3表示;加工三种裂解原料,轻石脑油(NAP),轻石脑油(LNAP)以及液化石油气(LPG),分别用A、B、C表示。调度时域H为200天。调度模型中所需要的相关参数如表1所示:
表1裂解炉炉群***的相关参数值
图3a为本发明实施例一提供的结焦量随运行时间的变化的一种趋势图。图3b为本发明实施例一提供的结焦量随运行时间的变化的另一种趋势图。图3c为本发明实施例一提供的结焦量随运行时间的变化的又一种趋势图。其中,NAP为石脑油,LNAP为轻石脑油,LPG为液化石油气,GK-VI代表GK-VI型号裂解炉,GK-V代表GK-V型号裂解炉,CBL-III代表CBL-III型号裂解炉。图4a为本发明实施例一提供的基础模型调度方案的一种示意图。图4b为本发明实施例一提供的基础模型调度方案的另一种示意图。图4c为本发明实施例一提供的基础模型调度方案的又一种示意图。原料批处理次序为:1号裂解炉:A-B-C;2号裂解炉:A-B-C;3号裂解炉:A-B-C。图5a为本发明实施例一提供的改进模型调度方案的一种示意图。图5b为本发明实施例一提供的改进模型调度方案的另一种示意图。图5c为本发明实施例一提供的改进模型调度方案的又一种示意图。原料批处理次序为:1号裂解炉:C-A-C-B;2号裂解炉:B-C-A;3号裂解炉:A-B-A-C。
如图4a、图4b、图4c与图5a、图5b、图5c所示,基础模型与改进模型均采用GAMS进行求解,求解器采用DICOPT。MINLP问题分解为MILP和NLP子问题,分别通过CPLEX求解器和CONOPT求解器完成求解。总循环时间为200天。基础模型平均结焦量为0.173kg/t乙烯,日平均利润为735070元/天。新的改进模型平均结焦量为0.164kg/t乙烯,日平均利润为764890元/天。日平均利润与平均结焦量均可得改进。
两个模型的对比在图3a、图3b、图3c与图4a、图4b、图4c最优调度优化方案可以发现,改进后的模型原料处理顺序和改进前有很大差异。实际生产原料并不会按照规定顺序一批一批加入,且同种原料不同批次的处理时间也不会完全相同。因此,对于原料批次的安排更加合理,其平均利润对比原模型也有一定的提高。
由表2的计算结果可知,考虑平均结焦量的裂解炉炉群调度建模优化方法可以在牺牲少量利润前提下,取得更高的结焦量减少,从而达到更好的环境效益。本发明充分考虑乙烯生产过程中非同步清焦的生产特性,从而保证调度的正常稳健进行。本发明可以有效地模拟乙烯化工炉群调度生产,有效提升环境效益,从而有效应对目前愈加严重的环境污染,可以对实际生产起到参考与指导作用。
表2旧模型与新模型的计算结果
本实施例提供一种考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法,包括:获取不同裂解原料投放进不同炉型裂解炉中的平均结焦量与乙烯产品量,获得调度模型的目标函数,建立约束条件,构造MINLP模型,使用DICOPT求解器对MINLP模型进行优化求解。本实施例综合考虑乙烯生产中产率随时间衰减的特性,兼顾裂解原料负荷的变化,为每个裂解炉规划调度时间范围内的批次、批处理时间、清焦顺序以及批次原料进料量的最佳安排。此外,本实施例可以在清焦阶段减少污染物的排放,在牺牲少量利润的代价下,获取更加可观的环境效益,为乙烯工厂的节能减排、优化排产提供了理论依据。
可以理解的是,以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式,然而本发明并不局限于此。对于本领域内的普通技术人员而言,在不脱离本发明的精神和实质的情况下,可以做出各种变型和改进,这些变型和改进也视为本发明的保护范围。
Claims (8)
1.一种考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法,其特征在于,包括:
获取不同裂解原料投放进不同炉型裂解炉中的平均结焦量与乙烯产品量;
根据所述平均结焦量与所述乙烯产品量,获得调度模型的目标函数,所述目标函数用于最小化单位乙烯产量的平均结焦量,所述目标函数的表达式如下:
建立约束条件,所述约束条件包括原料平衡约束、时间约束、整数约束、非同步清焦约束、变量上下限约束和附加约束;
考虑裂解原料负荷变化,通过模拟获得不同负荷水平下乙烯收率模型,所述裂解原料负荷仅在操作条件发生变化时在进料范围内进行调整优化,在单个调度周期内保持不变;
根据所述目标函数、所述约束条件和所述乙烯收率模型,形成炉群调度的MINLP模型;
使用DICOPT求解器对所述MINLP模型进行优化求解,使用子求解器CPLEX和子求解器CONOPT分别处理MILP问题和NLP问题。
4.根据权利要求3所述的考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法,其特征在于,所述时间约束的表达式如下:
所述进料i的连续处理时间设置在预设的最小值和最大值之内,若所述进料i未分配给第j炉的第k个批次生产运行,所述处理时间为0;
Tsj,k为所有裂解炉的裂解批次的开始时间,Tej,k为所有裂解炉的裂解批次的结束时间,当前批次的裂解开始时间点等于上一批次的裂解结束时间点加上两批次之间的停炉清焦时间,当前批次的结束时间点等于当前批次的开始时间点加上批处理时间,所有裂解炉第一批次的开始时间均为0,若批处理启动时间点大于调度时间范围H,开始时间Tsj,k等于0;
引入二进制变量yi,j,k表示所述进料i是否分配给第j炉的第k个批次生产运行,若所述裂解原料已经分配,二进制变量yi,j,k为1,若所述裂解原料尚未分配,二进制变量yi,j,k为0,若第k个批次没有生产运行,所述第k个批次的开始和结束时间相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011340916.3A CN112669912A (zh) | 2020-11-25 | 2020-11-25 | 考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011340916.3A CN112669912A (zh) | 2020-11-25 | 2020-11-25 | 考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112669912A true CN112669912A (zh) | 2021-04-16 |
Family
ID=75403690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011340916.3A Pending CN112669912A (zh) | 2020-11-25 | 2020-11-25 | 考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112669912A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113380339A (zh) * | 2021-06-17 | 2021-09-10 | 中国环境科学研究院 | 区域空间污染特征的检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103289725A (zh) * | 2012-03-05 | 2013-09-11 | 中国石油天然气股份有限公司 | 乙烯裂解炉的节能优化控制方法 |
CN103524284A (zh) * | 2013-10-14 | 2014-01-22 | 中国石油化工股份有限公司 | 一种乙烯裂解原料配置的预测和优化方法 |
CN103605821A (zh) * | 2013-09-16 | 2014-02-26 | 华东理工大学 | 乙烯裂解炉群负荷分配优化方法 |
CN104537430A (zh) * | 2014-12-15 | 2015-04-22 | 浙江大学 | 一种基于启发式算法的乙烯厂的生产控制方法 |
CN104611000A (zh) * | 2014-12-25 | 2015-05-13 | 东北大学 | 一种提高大型乙烯裂解炉作业效率的生产批调度控制方法 |
CN111639793A (zh) * | 2020-05-13 | 2020-09-08 | 新奥数能科技有限公司 | 一种锅炉群组调度优化方法及装置 |
CN111931975A (zh) * | 2020-06-19 | 2020-11-13 | 北京化工大学 | 一种考虑下游扰乱约束下的裂解炉炉群调度建模与方法 |
-
2020
- 2020-11-25 CN CN202011340916.3A patent/CN112669912A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103289725A (zh) * | 2012-03-05 | 2013-09-11 | 中国石油天然气股份有限公司 | 乙烯裂解炉的节能优化控制方法 |
CN103605821A (zh) * | 2013-09-16 | 2014-02-26 | 华东理工大学 | 乙烯裂解炉群负荷分配优化方法 |
CN103524284A (zh) * | 2013-10-14 | 2014-01-22 | 中国石油化工股份有限公司 | 一种乙烯裂解原料配置的预测和优化方法 |
CN104537430A (zh) * | 2014-12-15 | 2015-04-22 | 浙江大学 | 一种基于启发式算法的乙烯厂的生产控制方法 |
CN104611000A (zh) * | 2014-12-25 | 2015-05-13 | 东北大学 | 一种提高大型乙烯裂解炉作业效率的生产批调度控制方法 |
CN111639793A (zh) * | 2020-05-13 | 2020-09-08 | 新奥数能科技有限公司 | 一种锅炉群组调度优化方法及装置 |
CN111931975A (zh) * | 2020-06-19 | 2020-11-13 | 北京化工大学 | 一种考虑下游扰乱约束下的裂解炉炉群调度建模与方法 |
Non-Patent Citations (3)
Title |
---|
CHAOWEI LIU, ET AL.: "Cyclic scheduling for best profitability of industrial cracking furnace system", 《COMPUTERS & CHEMICAL ENGINEERING》, 22 October 2009 (2009-10-22), pages 544 - 553 * |
JIANG P, ET AL.: "Multi-objective modeling and optimization for scheduling of cracking furnace systems", 《CHINESE JOURNAL OF CHEMICAL ENGINEERING》, 8 July 2017 (2017-07-08), pages 992 - 998 * |
简海东 等: "兼顾负荷变化的裂解炉炉群调度建模与优化", 《控制工程》, vol. 21, no. 6, 30 November 2014 (2014-11-30), pages 944 - 948 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113380339A (zh) * | 2021-06-17 | 2021-09-10 | 中国环境科学研究院 | 区域空间污染特征的检测方法 |
CN113380339B (zh) * | 2021-06-17 | 2022-02-08 | 中国环境科学研究院 | 区域空间污染特征的检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Leung et al. | A review on the development and commercialization of biomass gasification technologies in China | |
Jiang et al. | Conserving water by optimizing production schedules in the dyeing industry | |
CN103605821B (zh) | 乙烯裂解炉群负荷分配优化方法 | |
US20090240378A1 (en) | Concentrating control method of generating sets for grid energy saving and reduction of so2 discharge | |
CN112183811B (zh) | 一种用于炼厂的多周期计划调度协同排产优化方法及*** | |
Zhao et al. | Production planning optimization of an ethylene plant considering process operation and energy utilization | |
CN102520705A (zh) | 一种炼化生产过程优化分析方法及*** | |
CN102213949A (zh) | 一种乙烯装置价值优化方法 | |
CN104537430B (zh) | 一种基于启发式算法的乙烯厂的生产控制方法 | |
CN104611000B (zh) | 一种提高大型乙烯裂解炉作业效率的生产批调度控制方法 | |
Su et al. | Scheduling of cracking production process with feedstocks and energy constraints | |
CN112669912A (zh) | 考虑平均结焦量与原料负荷的乙烯裂解炉炉群调度方法 | |
Liu et al. | A mixed-integer programming approach to strategic planning of chemical centres: A case study in the UK | |
CN111461406A (zh) | 一种基于混合建模的乙烯生产原料优化方法 | |
Cheung et al. | Short-term site-wide maintenance scheduling | |
CN117540972A (zh) | 一种针对化工行业电能供需波动的调节及电能能耗比提升的方法 | |
CN111931975A (zh) | 一种考虑下游扰乱约束下的裂解炉炉群调度建模与方法 | |
Foo et al. | A heuristic-based algebraic targeting technique for aggregate planning in supply chains | |
Lin et al. | Data-driven scheduling optimization of ethylene cracking furnace system | |
CN113780655A (zh) | 一种基于智慧供应链的钢材多品种需求预测方法 | |
Gong et al. | Interval-parameter bi-level programming for energy system management under uncertainty: Towards a deep-decarbonized and sustainable future in China | |
CN101441466A (zh) | 一种基于预测控制的聚丙烯多牌号生产计划的控制方法 | |
CN115455694A (zh) | 一种乙烯裂解炉炉群节能调度的改进人工蜂群优化方法 | |
Kychkin et al. | Automated Cross Channel Temperature Predictions for the PFR Lime Kiln Operating Support | |
Lin et al. | Cycle Scheduling of Ethylene Cracking Furnace System with Inventory Constraints |
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 |