运筹学案例1.1和1.2.doc

第一章

案例1.1炼油厂生产计划安排

某炼油厂的工艺流程图如下图所示。

炼油厂输入两种原油（原油1和原油2）.原油先进入蒸馏装置，每桶原油经蒸馏后的产品及份额见下表，其中轻、中、重石脑油的辛烷值分别是90、80和70.

?

轻石脑油

中石脑油

重石脑油

轻??油

重油

残油

原油1

0.10

0.20

0.20

0.12

0.20

0.13

原油2

0.15

0.25

0.18

0.08

0.19

0.12

问题:应如何制定炼油厂的生产计划，以得到最大利润。

问题分析

该炼油厂的最终生产产品是普通汽油、特级汽油、飞机燃料、燃料油、润滑油这5种成品油，用的最初原料是原油和原油，中间产品有重整汽油、裂化油、裂化汽油。由于各种成品油利润不同，所以这5种产品各生产数量，相互间的生产比例都会影响到总利润。现在，本题的目标是求出这5种成品油各生产多少时获得总利润最大。

为了求出5种成品油各生产多少时获得总利润最大，总的求解思路如下：

首先，求出总利润与5种成品油各生产数量的函数关系；

第二，找出生产这5种成品油所需原材料，以及各原材料在生产时的使用比例；

第三，由于炼油厂的生产实际能力有限，以及市场需求的原因，该炼油厂的原料油、中间产品和最终成品油的生产数量是受到约束的，所以要找出这些约束条件；

第四，根据目标方程和约束条件，编程实现求解，得到最大总利润。

二、模型建立与求解

模型的建立：

变量设置

用z1，z2，代表原油1和原油2的输入量，x1，x2，x3分别为由蒸馏后得到的轻、中、重石脑油的数量，xi1，xi2为i种石脑油用于高档和普通发动机油混合的数量，xi3为i种石脑油输入重整装置的量（i=1,2,3）。x4为由重整装置得到的重整汽油量，x41，x42为分别用于高档和普通发动机油混合的数量。轻油（x6）和重油（x7）有一部分输入裂解装置，一部分直接输入煤油、燃料油的混合装置。由蒸馏装置得到的渣油（x9）一部分输入渣油处理装置。提炼得润滑油；一部分输入煤油、燃料油混合装置。该炼油最终产品为y1，y2，……，y5。

目标函数

Maxz=700y1+600y2+400y3+350y4+150y5

(3)约束条件

=1\*GB3①各装置投入产出关系

如蒸馏装置有x1=0.7z1+0.15z2，x2=0.2z1+0.25z2，x3=0.2z1+0.18z2，

x7=0.12z1+0.08z2，x8=0.2z1+0.19z2，x9=0.13z1+0.12z2

又x1=x11+x12+x13，x2=x21+x22+x23，x3=x31+x32+x33，

x7=x71+x72+x73，x8=x81+x82+x83，x9=x91+x92+x93

对重整装置有x4=0.6x13+0.52x23+0.45x33，又x4=x41+x42+x43

对裂解装置有x5=0.28x73+0.2x83，又x5=x51+x52

x6=0.68x73+0.75x83，又x6=x61+x62

对渣油处理装置有y5=0.5x93

对发动机油混合装置有y1=x11+x21+x31+x41+x51

y2=x12+x22+x32+x42+x52

对煤油、燃料油混合装置有y3=x61+x71+x81+x91

y4=x62+x72+x82+x92

=2\*GB3②各装置能力限制

z1+z245000，x13+x23+x3310000

x73+x838000，

=3\*GB3③发动机油辛烷值限制

90x11+80x21+70x31+115x41+105x5194(x11+x21+x31+x41+x51)

90x12+80x22+70x32+115x42+105x5284(x12+x22+x32+x42+x52)

=4\*GB3④煤油气压的限制

1.5x61+1.0x71+0.6x81+0.05x91x61+x71+x81+x91

=5\*GB3⑤燃料油比例的限制

=6\*GB3⑥原油供应限制

=7\*GB3⑦最终产品数量限制

=8\*GB3⑧变量非负限制，所有变量均

模型的求解

采用li