炼油厂案例
炼油厂生产计划安排摘要
问题:本案例主要论述了炼油厂生产计划的优化问题,通过对炼油厂生产流程的分析,将原问题转
化为一个线性规划问题,进而建立了一个模型。
方法:解决本案例的方法是利用线性规划求解最优解,以及Lindo软件辅助求解。
结论:通过求得最优解,使该厂利润达到最大化,实现最优生产,并且通过敏感分析知道当各个终
产品的利润在什么变化幅度内可以保证生产计划不变,避免厂家因利润波动而草率的更改生产方案而导致费时费力。
关键词 炼油 线性规划 敏感性分析
背景分析
全球金融危机,使得人们对价格不断飙升的油产品需求急剧减少,虽然现在危机的云雾刚刚淡去,但炼油厂面临着的形势依旧是十分严峻的,要渡过这不景气的阶段,就需要炼油厂进行生产上的优化。在面对倒闭亦或是发展这两种结果的时候,该厂做出了自己的决策:在已有生产能力和市场条件的限制下,通过分析所有的生产流程,得到了总利润与成品油产量的函数关系,以及成品油产量所需要满足的约束条件,从而将问题转化为线性规划问题,然后转换建立了一个模型,靠这一模型得出的成产计划,能否帮助此厂脱离困境甚至收益颇丰呢,我们拭目以待……
理论综述
(名词解释及理论概述)
1.各组员分工:梁晓雲、黄桂凤,录入文本及编排文档;陈剑锋、管诗广,设计和制作幻灯片;梁莹,进行了背景分析,并根据整个案例总结摘要,提炼关键词;蔡婷婷,查阅文献,制作综述部分,为后面的解题方法做了理论上的铺垫;林穗阳,建模求解,归纳结论;廖冬燕,俯瞰通篇,审核,批改,修正案例。 2.原油:原油即石油,也称黑色金子,是一种粘稠的、深褐色(有时有点绿色的)液体。是石油刚开采出来未经提炼或加工的物质。地壳上层部分地区有石油储存。它由不同的碳氢化合物混合组成,其主要组成成分是烷烃,此外石油中还含硫、氧、氮、磷、钒等元素。不过不同油田的石油成分和外貌可以有很大差别。 3.石脑油:石脑油(naphtha):一部分石油轻馏分的泛称。因用途不同有各种不同的馏程。我国规定馏程自初镏点至220℃左右。主要用作重整和化工原料。作为生产芳烃的重整原料,采用70~145℃馏分,称轻石脑油;当以生产高辛烷值汽油为目的时,采用70~180℃馏分,称重石脑油。用作溶剂时,则称溶剂石脑油,来自煤焦油的芳香族溶剂也称重石脑油或溶剂石脑油。
4.辛烷值:汽油辛烷值是汽油在稀混合气情况下抗爆性的表示单位,在数值上等于在规定条件下与试样抗爆性相同时的标准燃料中所含异辛烷的体积百分数。
5. 轻油:轻油,又称石脑油、白电油,英文名:Naphtha。是石油提炼后的一种油质的产物。它由不同的碳氢化合物混合组成。它的主要成分是含5到11个碳原子的链烷、环烷或芳烃。石脑油可用作提炼煤气之用。
油富含芳烃,经过加氢精制后可作为高辛烷值汽油组分或用于萃取苯、甲苯、乙苯、二甲苯等化工原料。 8.线性规划:线性规划是运筹学中研究较早、发展较快、应用广泛、方法较成熟的一个重要分支,它是辅助人们进行科学管理的一种数学方法.研究线性约束条件下线性目标函数的极值问题的数学理论和方法,英文缩写LP。它是运筹学的一个重要分支,广泛应用于军事作战、经济分析、经营管理和工程技术等方面。为合理地利用有限的人力、物力、财力等资源作出的最优决策,提供科学的依据。
9.LINDO:LINDO(Linear, INteractive, and Discrete Optimizer)是一个解决二次线性整数规划问题的方便而强大的工具。
10. 敏感性分析:敏感性分析是指从定量分析的角度研究有关因素发生某种变化对某一个或一组关键指标影响程度的一种不确定分析技术。其实质是通过逐一改变相关变量数值的方法来解释关键指标受这些因素变动影响大小的规律。
(前人之研究及相应文献)
11.石油化工过程人工介入控制系统的应用研究:
1).对国内外石油化工工业控制现状及发展趋势进行了介绍。
2).从石油化工生产过程的特点出发,分析了石油化工生产过程人机系统中人的因素,并对人的可靠性进行了
定量分析。 3).在分析石油化工过程控制的基础上,提出人工介入控制系统的实质是监督控制,其重要功能在于优化操作。
对人工介入控制系统的基本理论进行了介绍,提出了石油化工过程的人工介入控制系统模型和结构,对操作行为的三个阶段及过程实时状态进行了描述,提出了石化企业的组织分析方法和石化企业的组织再设计的原则。
4).介绍了人工介入生产过程控制系统在中石化长岭分公司的下属企业---湖南长盛石油化工公司的实践。电
子巡检系统与设备状态网络诊断系统的成功应用表明:在石油化工过程中引入人工介入控制系统,对提高企业经济效益、提升企业管理水平均具有积极意义,可以达到了安全、经济、科学、实效管理的目的。
1 张志檩;严格在线模拟技术[J];石油化工;1996年05期 2 3 5
;2001年《石油化工动态》将改为月刊并更名为《当代石油化工》[J];炼油设计;2000年10期
徐用懋,杨尔辅;石油化工流程模拟、先进控制与过程优化技术的现状与展望[J];工业控制计算机;2001年09期
张盈珍 ,魏昭彬;第六届全国石油化工催化会议在浙江镇海召开[J];催化学报;1992年02期
4 ;《石油化工工学丛书》陆续出版[J];石油化工设计;1995年03期
6 高文宇,张力;数据挖掘技术在人因分析中的应用[J];人类工效学;2001年04期 7 ;著名石油化工专家王松汉教授[J];科学中国人;2004年07期
8 白尔铮;第三届石油化工催化会议在昆明举行[J];催化学报;1986年01期 9
裴瑞凌,荣冈;炼油过程的智能工厂流程模拟仿真平台[J];化工自动化及仪表;2005年02期
10 谢月洁;计算机在化工设计中的应用[J];内蒙古石油化工;2004年03期
12.炼油工业生产计划系统的研究
为了改进炼油厂制定生产计划的决策过程,提高炼厂的科学管理水平,研究并开发了基于炼油过程仿真的优化生产计划系统,它使计划编制和炼油工艺过程及市场对油品的需求状况紧密结合起来,从而提高了生产计划系统的实用性,使制定的生产计划更准确地指导炼厂的实际生产,提高了企业的经济效益。本文给出了基于炼油过程仿真优化生产计划系统的框架结构、优化原理、应用情况和其在FR-CIMS中的位置。
1
李雷,江燕斌,钱宇,罗金生;多应用集成的原油混炼操作支持系统[J];炼油技术与工程;2003年01期
2 杨家本,陈武,康飚,兰鸿森;计算机集成制造系统环境下流程工业企业管理模式重构的研究[J];石油炼制与化工;2000年03期 3 4
丁锋,杨家本,兰鸿森,黄河清,胡鑫;基于炼油过程仿真的优化生产计划系统的开发与应用[J];石油炼制与化工;1998年11期 丁锋,杨家本,兰鸿森,黄河清,毛福忠,胡鑫,肖纲,苏晓俊;连续过程生产调度5日滚动作业计划的研究[J];系统工程理论与实
践;1998年11期
5 6 7
宋洁蔚;油品生产与储运调度问题研究[D];浙江大学;2003年
赵向海;石化企业综合自动化信息集成平台及应用研究[D];浙江大学;2004年
赵小强;炼厂生产调度问题研究[D];浙江大学;2005
建模
某炼油厂的工艺流程图如图1所示。
图1
炼油厂输入两种原油 2(原油1和原油2)。原油先进入蒸馏装置,每桶原油经蒸馏后的产品份额见表1,其中轻、中、重石脑油 3的辛烷值 4分别为90、80和70。
cji
原油1 原油2
轻石脑油 0.1 0.15
中石脑油 0.2 0.25
重石脑油 0.2 0.18
表1
轻油 0.12 0.08
重油 0.2 0.19
渣油 0.13 0.12
石脑油部分直接用于发动机油混合,部分输入重整装置,得辛烷值为115的重整汽油。1桶轻、中、重石脑油经重整后得到的重整汽油分别为0.6、0.52、0.45桶。
蒸馏得到的轻油 5和重油 6,一部分直接用于煤油和燃料油的混合,一部分经裂解装置得到裂解汽油 7和裂解油。裂解汽油的辛烷值为105。1桶轻油经裂解后得到0.28桶裂解汽油和0.68桶裂解油;1桶重油裂解后得到0.2桶裂解汽油和0.75桶裂解油。其中裂解汽油用于发动机油混合,裂解油用于煤油和燃料油的混合。
渣油可直接用于煤油和燃料油的混合,或用于生产润滑油。1桶渣油经处理后可得0.5桶润滑油。 混合成的发动机油高档的辛烷值应不低于94,普通的辛烷值不低于84。混合物的辛烷值按混合前各油料辛烷值和所占比例线性加权计算。
规定煤油的气压不准超过1kg/cm2,而轻油、重油、裂解油和渣油的气压分别为1.0、0.6、1.5和0.05kg/cm2。而气压的计算按各混合成分的气压和比例线性加权计算。
燃料油中,轻油、重油、渣油和裂解油的比例应为10:3:1:4。
已知每天可供原油1为20000桶,原油2为30000桶。蒸馏装置能力每天最大为45000桶,重整装置每天最多重整10000桶石脑油,裂解装置能力每天最大为8000桶。润滑油每天产量应在500~1000桶之间,高档发动机油产量应不低于普通发动机油的40%。
又知最终产品的利润(元/桶)分别为:高档发动机油700,普通发动机油600,煤油400,燃料油350,润滑油150.
试为该炼油厂制订一个使总盈利为最大的计划,并求当各种终产品的利润分别在什么范围时,生产计划保持不变。
解:
这是一个线性规划 8的问题。
设使用原油桶数为
zi(i=1,2);分馏后的产物轻、中、重石脑油、轻油、重油桶数分别为xi(i=1,2,…,6);
石脑油输入到重整装置后得到的重整汽油桶数为x7,轻油和重油经裂解装置得到裂解汽油和裂解油分别为
x8、x9;石脑油和重整汽油、裂解汽油直接用于发动机油混合来生产高档、普通发动机油的桶数为xi,j
(i=1,2,3,7,8;j=1,2,其中1表示高档,2表示普通);轻油、重油、渣油、裂解油用于煤油和燃料油的混合生产从而得到煤油或者燃料油的桶数为
xm,n (m=4,5,6,9;n=1,2,其中1表示煤油,2表示燃料油);石脑油
用于输入重整装置的桶数为
xi3(i=1,2,3),轻、重油用于输入裂解装置的桶数分别为x43、x53,渣油用
于生产润滑油的桶数为
x63;最终产物高档、普通发动机油、煤油、燃料油、润滑油桶数为yj (j=1,…,5);
cji ( j =1,2; i =1,2,…,6),1桶轻、中、重石脑油经重整后
j
又设1桶原油经过蒸馏后得到的产物桶数为
得到的重整汽油桶数分别为i(i=1,2,3), 1桶轻油、重油经裂解后得到的裂解汽油为m(m=4,5),1桶轻油、重油经裂解后得到的裂解油桶数为m (m=4,5),轻、中、重石脑油、重整汽油、裂解汽油的辛烷值为
r
k
wi(i=1,2,3,7,8), 轻油、重油、裂解油和渣油的气压为pm(m=4,5,6,9),最终产物的利润分别为lj
(j=1,2,…,5),则:
Max L=
ly
jj1
2j1
5
j
St:蒸馏装置(来源): xizjcji i1,2…,,6
又(去向):
xixi t i1,2…,,6
t1
3
3
重整装置(来源):
x7rixi3
i1
又(去向):
x7x7j
j1
5
2
裂解装置(来源):
x8jmxm3 , x9kmxm3
m4
m4
5
又(去向):
x8x8j , x9x9n
j1
n1
2
2
对发动机油混合装置、煤油、燃料油混合装置以及润滑油装置有:
yj[(j3)(j4)(j5)](j1)(j2)xi,j[(j1)(j2)(j5)](j3)(j4)xm,n[(j1)(j2)(j3)(j4)]j50.5x63
(j=1,2,…,5;i=1,2,3,7,8;m=4,5,6,9;n=1,2)
xi1
wi94 (i=1,2,3,7,8); 普通——辛烷值: 高档——iy1
煤油气压:
xi2
wi84 (i=1,2,3,7,8) iy2
xm1
pm1 (m=4,5,6,9) my3
燃料油比例:
x42x52x62x10314
3
92
原油供应: 装置能力限制:
最终产品产量:
z120000,z230000
z1z245000 ,xi310000 ,x43x538000
i1
500y51000,y10.4y2
另外所有变量均非负。
模型分析
软件求解
将上述模型中的变量代入具体的数值,列齐所有方程,并用Lindo 9求解并进行敏感性分析 10。Lindo的录入方程如下:
max 700y1+600y2+400y3+350y4+150y5st
蒸馏装置:x1-0.1z1-0.15z2=0 x2-0.2z1-0.25z2=0 x3-0.2z1-0.18z2=0 x4-0.12z1-0.08z2=0 x5-0.2z1-0.19z2=0 x6-0.13z1-0.12z2=0 x1-x11-x12-x13=0 x2-x21-x22-x23=0 x3-x31-x32-x33=0 x4-x41-x42-x43=0 x5-x51-x52-x53=0 x6-x61-x62-x63=0重整装置及裂解装置:x7-0.6x13-0.52x23-0.45x33=0 x7-x71-x72=0 x8-0.28x43-0.2x53=0 x8-x81-x82=0 x9-0.68x43-0.75x53=0 x9-x91-x92=0另外的三个装置:y1-x11-x21-x31-x71-x81=0 y2-x12-x22-x32-x72-x82=0 y3-x41-x51-x61-x91=0 y4-x42-x52-x62-x92=0 y5-0.5x63=0辛烷值: 高档—— 90x11+80x21+70x31+115x71+105x81-94y1>0 普通—— 90x12+80x22+70x32+115x72+105x82-84y2>0气压: 1.5x91+x41+0.6x51+0.05x61-y3
比例: 3x42-10x52=0 4x52-3x92=0 x92-4x62=0 原油供应:z1
装置能力限制:z1+z2500 y1-0.4y2>0end
结果如下:
LP OPTIMUM FOUND AT STEP 26 OBJECTIVE FUNCTION VALUE 1) 0.2113651E+08
VARIABLE VALUE REDUCED COST Y1 6817.778809 0.000000 Y2 17044.447266 0.000000 Y3 15156.000000 0.000000 Y4 0.000000 0.000000 Y5 500.000000 0.000000
ROW SLACK OR SURPLUS DUAL PRICES 2) 0.000000 665.376221 3) 0.000000 548.270020 4) 0.000000 431.163818 5) 0.000000 439.282837 6) 0.000000 400.000000 X1 6000.000000 0.000000 Z1 15000.000000 0.000000 Z2 30000.000000 0.000000 X2 10500.000000 0.000000 X3 8400.000000 0.000000 X4 4200.000000 0.000000 X5 8700.000000 0.000000 X6 5550.000000 0.000000 X11 0.000000 0.000000 X12 6000.000000 0.000000 X13 0.000000 90.491142 X21 3537.524414 0.000000 X22 6962.475586 0.000000 X23 0.000000 50.036312 X31 0.000000 0.000000 X32 2993.138184 0.000000 X33 5406.861816 0.000000 X41 0.000000 39.282845 X42 0.000000 0.000000 X43 4200.000000 0.000000 X51 4900.000000 0.000000 X52 0.000000 0.000000 X53 3800.000000 0.000000 X61 4550.000000 0.000000 X62 0.000000 0.000000 X63 1000.000000 0.000000 X7 2433.087891 0.000000 X71 1344.254395 0.000000 X72 1088.833374 0.000000 X8 1936.000000 0.000000 X81 1936.000000 0.000000 X82 0.000000 0.000000 X9 5706.000000 0.000000 X91 5706.000000 0.000000 X92 0.000000 323.207123
7) 0.000000 400.000000 8) 0.000000 -665.376221 9) 0.000000 -548.270020 10) 0.000000 -431.163818 11) 0.000000 -439.282837 12) 0.000000 -400.000000 13) 0.000000 -400.000000 14) 0.000000 958.141785 15) 0.000000 -958.141785 16) 0.000000 841.035583 17) 0.000000 -841.035583 18) 0.000000 400.000000 19) 0.000000 -400.000000 20) 0.000000 -388.579712 21) 0.000000 -388.579712 22) 0.000000 400.000000 23) 0.000000 350.000000 24) 0.000000 800.000000 25) 0.000000 -11.710622 26) 0.000000 -11.710622 27) 3429.500000 0.000000 28) 0.000000 -29.760948 29) 0.000000 -86.902374 30) 0.000000 12.500000 31) 5000.000000 0.000000 32) 0.000000 26.487722 33) 0.000000 447.138336 34) 4593.138184 0.000000 35) 0.000000 68.207115 36) 500.000000 0.000000 37) 0.000000 -650.000000 38) 0.000000 -12.218727
NO. ITERATIONS= 26
RANGES IN WHICH THE BASIS IS UNCHANGED: OBJ COEFFICIENT RANGES VARIABLE CURRENT ALLOWABLE ALLOWABLE COEF INCREASE DECREASE Y1 700.000000 18.067978 892.088684 Y2 600.000000 INFINITY 15.097159 Y3 400.000000 272.828461 74.729973 Y4 350.000000 71.823807 INFINITY Y5 150.000000 650.000000 INFINITY X1 0.000000 INFINITY 529.754456 Z1 0.000000 26.487722 447.138336 Z2 0.000000 INFINITY 26.487722 X2 0.000000 INFINITY 529.754456 X3 0.000000 1324.386108 2235.691650 X4 0.000000 662.193054 3726.152832 X5 0.000000 2648.772217 2235.691650 X6 0.000000 2648.772217 3439.525879 X11 0.000000 0.000000 INFINITY X12 0.000000 INFINITY 0.000000 X13 0.000000 90.491127 INFINITY X21 0.000000 0.000010 0.000000 X22 0.000000 0.000000 0.000010 X23 0.000000 50.036308 INFINITY X31 0.000000 0.000002 INFINITY X32 0.000000 68.731743 0.000002 X33 0.000000 50.500000 68.731743 X41 0.000000 39.282837 INFINITY X42 0.000000 129.282852 INFINITY X43 0.000000 662.193054 39.282837 X51 0.000000 68.207115 39.282837 X52 0.000000 430.942841 INFINITY X53 0.000000 39.282837 68.207115 X61 0.000000 2648.772217 325.000000 X62 0.000000 1292.828491 INFINITY X63 0.000000 325.000000 INFINITY X7 0.000000 112.222221 6894.907715 X71 0.000000 0.000004 0.000000 X72 0.000000 0.000000 0.000004 X8 0.000000 8277.413086 341.035583 X81 0.000000 8277.413086 0.000003 X82 0.000000 0.000003 INFINITY X9 0.000000 561.183411 90.942818 X91 0.000000 561.183411 90.942818 X92 0.000000 323.207123 INFINITY
单纯形表最终格式不再列出,篇幅较大,可在附加的txet
结果分析
RIGHTHAND SIDE RANGES
ROW CURRENT ALLOWABLE ALLOWABLE RHS INCREASE DECREASE 2 0.000000 6447.203613 6000.000000 3 0.000000 12886.083008 8186.603516 4 0.000000 8043.830078 6097.761719 5 0.000000 3800.000000 4200.000000 6 0.000000 INFINITY 4900.000000 7 0.000000 INFINITY 3610.000000 8 0.000000 6000.000000 6447.203613 9 0.000000 8186.603516 12886.083008 10 0.000000 6097.761719 8043.830078 11 0.000000 4200.000000 3800.000000 12 0.000000 4900.000000 INFINITY 13 0.000000 3610.000000 INFINITY 14 0.000000 5671.715820 2743.992920 15 0.000000 2743.992920 5671.715820 16 0.000000 2394.595459 1936.000000 17 0.000000 1936.000000 2394.595459 18 0.000000 6859.000000 5706.000000 19 0.000000 5706.000000 6859.000000 20 0.000000 1017.322388 606.869507 21 0.000000 1017.322388 780.733582 22 0.000000 INFINITY 3429.500000 23 0.000000 INFINITY 0.000000 24 0.000000 500.000000 1805.000000 25 0.000000 88361.710938 50837.417969 26 0.000000 88361.710938 62676.164062 27 0.000000 INFINITY 3429.500000 28 0.000000 12600.000000 0.000000 29 0.000000 4523.076660 0.000000 30 0.000000 0.000000 14440.000000 31 20000.000000 INFINITY 5000.000000 32 30000.000000 15000.000000 5000.000000 33 45000.000000 5000.000000 13412.202148 34 10000.000000 INFINITY 4593.138184 35 8000.000000 4425.161621 3800.000000 36 1000.000000 INFINITY 500.000000 37 500.000000 500.000000 500.000000 38 0.000000 6736.749023 4969.412598
文件中打开查阅,文件名为”lindo分析结果——炼油厂案例”
1.根据第一部分(Objective function value 目标函数值)结果的第一列可知,最优生产计划是:每日使用的原油1和原油2的桶数分别为15000、30000,可分别得到轻、中、重石脑油,轻、重油,渣油,重整汽油,裂解汽油,裂解油6000、10500、8400、4200、8700、5550、2433.087891、1936、5706桶,以及最终产物高档、普通发动机油,煤油,燃料油,润滑油分别6817.778809、17044.447266、15156、0、500桶,可以注意到该厂将不生产燃料油,在这种生产计划下,公司获得的最大利润是21136510元/天。(第二列是松弛变量及其对应的对偶价格(影子价格),很明显由于有37个方程,故有37个松弛变量,而本题所设的变量有40个,故不能确定各个对偶价格分别对应哪个变量,本人能力尚浅,不会分析,请老师指导)
2.第二部分(RANGES IN WHICH THE BASIS IS UNCHANGED 最优基不变的范围)敏感性分析结果的第一列是目标函数中的变量系数范围,目标函数中各变量的系数在可行增量和可行减量之间变动可保证最优解不变,以Y1为例,其当前系数(current coef)为700,可行增量是18.067978,可行减少量是892.088684,所以在最优解不变的前提下Y1的系数可以变化的范围是(-192.089,718.067978),其他分析情况同上,注意,在临界值处可能会出现奇异性。第二部分结果的第二列是限制条件中右边系数的变化范围,右边系数在可行增量和可行减量之间变动可保证最优基不变,以第30个限制条件(即该列的第31行)为例,其当前系数(current rhs)是20000,可行增量是+∞,可行减少量是5000,所以在最优基不变的前提下该限制条件可以变化的范围是(15000, +∞),其他分析情况同上,注意,在临界值处可能会出现奇异性。
结论
最优生产计划是:每日使用的原油1和原油2的桶数分别为15000、30000,可得到最大利润21136510元/天;高档发动机油,普通发动机油,煤油,燃料油,润滑油的利润分别再(-192.089,718.067978),( 584.9028,+∞),( 325.27,672.828461),(-∞,421.823807),(-∞, 800)范围内时生产计划不变。