3炼油厂优化问题-GoogleCode

1、数学模型课程设计学校：华东交通大学课题：炼油厂的优化专业： 信息与计算科学(1) 班姓名：刘文清学号：010121组员：刘文清、刘国涛老师：朱旭生目录一、 内容摘要3二、 提出问题5三、 模型假设6四、 变量说明6五、 问题分析及模型建立7六、 模型求解9七、 参考文献12八、 个人体会12一、内容摘要内容提要：本模型针对炼油厂的生产和销售状况进行优化，在限制条件以内进行优化，争取实现炼油厂效益的最大化。运用线性规划建立了炼油厂生产的优化模型，再利用 lindo和 matlab 对本规划问题进行求解，得出最优化方案，从而实现队炼油厂生产状况的最优化方案的确定。依据题目有，某炼油厂购买两种原油（

2、C1 和 C2），经过蒸馏、重整、裂解和混合等 4 种加工过程，生成汽油和燃料油出售。 其间可以通过重整， 裂解生产中间产品，然后可以混合和生产润滑油出售。整个初级产品的生产过程如下图：原油蒸馏石油精油残油重整裂解润滑油重整汽油裂化油裂化汽油由初级产品生产最终产品如下图：汽油混合石油精重整汽油裂化汽油飞机燃料油裂化油残油燃料油轻油裂化油重油残渣二、提出问题某炼油厂购买两种原油 (C1 和C2)，经过蒸馏、重整、裂解和混合等4 种加工过程，生成汽油和燃料油出售。蒸馏 该过程将每种原油按沸点不同分离成下列分馏物：轻石油精、中石油精、重石油精、轻油、重油和残留。轻、中和重石油精的辛烷值分别为 90、

3、80和 70。一桶原油生成分馏物的桶数如下表所示：分馏物石油精油残渣原油轻中重轻重C10.10.20.20.120.20.13C20.150.250.180.080.190.12注：蒸馏过程中有少量损耗。重整 石油精可直接用于混合产生各种品级的汽油， 也可以再经过一种称作重整的加工过程，生成一种称为重整汽油的产品，其辛烷值为 115。不同的石油精每桶产生重整汽油的桶数如下表 :石油精轻中重重整汽油0.60.520.45裂解轻油和重油可以直接混合生产飞机燃料和燃料油，也可以经过催化裂解，生成裂化油和裂化汽油，后者的辛烷值为 105。一桶轻油经裂解生成裂化油 0.68 桶和裂化汽油 0.28 桶；

4、一桶重油经裂解生成裂化油 0.75 桶和裂化汽油 0.2 桶。裂化油用于混合生成飞机燃料和燃料油；裂化汽油用于混合生成汽油。残渣用于生产润滑油或混合入飞机燃料和燃料油中， 一桶残渣产生 0.5 桶润滑油。混合汽油（发动机燃料）汽油有两种： 普通汽油和特级汽油， 通过将石油精、 重整汽油和裂化汽油混合得到。普通汽油必须有不小于 84 的辛烷值；特级汽油必须有不小于 94 的辛烷值。假设辛烷值是按体积线性地混合的， 即混合物的辛烷值是各组分的辛烷值以组分的体积为权的加权平均。飞机燃料飞机燃料的蒸发压不能超过1kg/cm2。轻油、重油、裂化油和残渣的蒸发压依次为 1.0、 0.6、1.5 和 0.0

5、5kg/cm2。同样假设蒸发压是按体积线性地混合的。燃料油轻油、裂化油、重油和残渣按10:4:3:1 的比例混合生成燃料油。关于原料供应和加工能力的限制条件有:(1)C1 日供应 20000 桶；(2) C2 日供应 30000 桶；(3)日蒸馏原油最多 45000 桶；(4)日重整石油精最多10000 桶；(5)日裂解油最多 8000 桶；(6)日产润滑油必须在500 桶至 1000桶之间；(7) 特级汽油的产量必须是普通汽油的 40%。每售出一桶成品，利润数如下表所示：产品特级汽油普通汽油飞机燃料燃料油润滑油利润600 元500 元350 元300 元120 元炼油厂当如何安排生产，总利润

6、为最大？三、模型假设（ 1）假设油类产品销路很好，只要在限制条件范围内生产就不会存在积压的状况。（ 2）假设由原料加工成产品时，产品的体积是原料体积之和。（ 3）假设生产普通汽油和特级汽油时，只要辛烷值达到要求即可， 5 种原材料不需要都用到。（ 4） 假设生产飞机燃料时，只要蒸发压达到要求即可， 4 种原材料不需要都用到。（ 5）假设炼油厂购买的石油全部用于蒸馏，即买入量与蒸馏量相等。四、变量说明 ( 单位：桶 )定义以下变量：原料： 设原油 C1 产量 X1，原油 C2产量为 X2；石油精：轻 X3重整 Y3未处理 W3 生产特级汽油 W31普通汽油 W32中 X4重整 Y4未处理 W4

7、生产特级汽油 W41普通汽油 W42重 X5重整 Y5未处理 W5 生产特级汽油 W51普通汽油 W52油：轻 X6裂解 Y6未处理 W6 生产飞机燃油 W61燃料油 W62重 X7裂解 Y7未处理 W7 生产飞机燃油 W71燃料油 W72残渣： 总 X8制润滑油 Y8未处理 W8 生产飞机燃油 W81 燃料油 W82重整汽油：总 C1生产特级汽油 C11普通汽油 C12裂化油：总 C2生产飞机燃油 C21燃料油 C22裂化汽油：总 C3生产特级汽油 C31普通汽油 C32最后特级汽油总量为 Z1，普通汽油总量为 Z2，飞机燃料为 Z3，燃料油为 Z4，润滑油为 Z5;五、问题分析及建立模型该

8、炼油厂的最终生产产品是普通汽油、特级汽油、飞机燃料、燃料油、润滑油这 5 种成品油，用的最初原料是原油和原油，中间产品有重整汽油、裂化油、裂化汽油。由于各种成品油利润不同， 所以这 5 种产品各生产数量， 相互间的生产比例都会影响到总利润。 现在，本题的目标是求出这5 种成品油各生产多少时获得总利润最大。为了求出 5 种成品油各生产多少时获得总利润最大，总的求解思路如下：首先，求出总利润与5 种成品油各生产数量的函数关系；第二，找出生产这 5 种成品油所需原材料， 以及各原材料在生产时的使用比例；第三，由于炼油厂的生产实际能力有限， 以及市场需求的原因， 该炼油厂的原料油、中间产品和最终成品油

9、的生产数量是受到约束的， 所以要找出这些约束条件；第四，根据目标方程和约束条件，编程实现求解，得到最大总利润。目标函数：max=600*Z1+500*Z2+350*Z3+300*Z4+120*Z5依据题目中的约束条件有：原材料约束：X1=20000X2=30000X1+X2=45000X3=0.1*X1+0.15*X2X4=0.2*X1+0.25*X2X5=0.2*X1+0.18*X2X6=0.12*X1+0.08*X2X7=0.2*X1+0.19*X2X8=0.13*X1+0.12*X2生产能力约束：Y3+Y4+Y5=10000Y6+Y7=500Y8=94* (C11+C21+W31+W41

10、+W51)普通汽油：Z2=W32+W42+W52+C12+C32115*C12+105*C32+90*W32+80*W42+70*W52=84* (C12+C22+W32+W42+W52)飞机燃油：Z3=W61+W71+W81+C21X61+0.6*X71+1.5*C21+0.05*W81=X61+X71+C21+W81燃料油：Z4=10*W62+4*C22+3*W72+W82润滑油：Z5=Y8Z5=0.5*X8数量关系：X3=Y3+W3W3=W31+W32X4=Y4+W4W4=W41+W42X5=Y5+W5W5=W51+W52X6=Y6+W6W6=W61+W62X7=Y7+W7W7=W71+

11、W72X8=Y8+W8W8=W81+W82六、模型求解在 lingo软件中输入如下程序：model:max=600*Z1+500*Z2+350*Z3+300*Z4+120*Z5;X1=20000;X2=30000;X1+X2=45000;X3=0.1*X1+0.15*X2;X4=0.2*X1+0.25*X2;X5=0.2*X1+0.18*X2;X6=0.12*X1+0.08*X2;X7=0.2*X1+0.19*X2;X8=0.13*X1+0.12*X2;Y3+Y4+Y5=10000;Y6+Y7=500;Y8=94* (C11+C21+W31+W41+W51);115*C12+105*C32+9

12、0*W32+80*W42+70*W52=84* (C12+C22+W32+W42+W52);Z2=W32+W42+W52+C12+C32;Z3=W61+W71+W81+C21;X61+0.6*X71+1.5*C21+0.05*W81=X61+X71+C21+W81; Z4=10*W62+4*C22+3*W72+W82;X3=Y3+W3;W3=W31+W32;X4=Y4+W4;W4=W41+W42;X5=Y5+W5;W5=W51+W52;X6=Y6+W6;W6=W61+W62;X7=Y7+W7;W7=W71+W72;X8=Y8+W8;W8=W81+W82;end利用 lingo 软件编程得到结果

13、如下：VariableValueReduced CostZ16785.2010.000000Z216963.000.000000Z32900.0000.000000Z427692.310.000000Z51000.0000.000000X115384.620.000000X20.000000-160.4249X31538.4620.000000X43076.9230.000000X53076.9230.000000X61846.1540.000000X73076.9230.000000X82000.0000.000000Y30.0000001920.476Y40.0000002096.667Y

14、50.0000000.000000Y60.0000003000.000Y70.000000900.0000Y81000.0000.000000C10.0000000.000000C20.0000000.000000C30.0000000.000000W311538.4620.000000W413076.9230.000000W510.000000-2272.857C110.000000-1480.000C312169.8170.000000C211900.0000.000000C120.0000000.000000C3213886.080.000000W320.0000001920.476W4

15、20.0000002096.667W523076.9230.000000C220.000000-1200.000W610.0000002650.000W710.000000550.0000W811000.0000.000000X610.0000000.000000X710.000000-1604.952W621846.1540.000000W723076.9230.000000W820.0000003861.762W31538.4620.000000W43076.9230.000000W53076.9230.000000W61846.1540.000000W73076.9230.000000W

16、81000.0000.000000最大利润元，在获得最大利润时，炼油厂的生产计划如下表：项目数量（桶）特级汽油6785.2普通汽油16963.0飞机燃油2900.0燃料油27692.0润滑油1000.0七、参考文献1、数学模型（第三版）姜启源谢金星叶俊编高等教育出版社。八、个人体会学习数学模型这门课已经一个学期了，知道了数学建模的大概过程，其过程如下：1、模型准备 : 了解问题的实际背景，明确其实际意义，掌握对象的各种信息，用数学语言来描述问题。2、模型假设 ; 根据实际对象的特征和建模的目的， 对问题进行必要的简化，并用精确的语言提出一些恰当的假设。3、模型建立： 在假设的基础上， 利用适当的数学工具来规划各变量之间的数学关系，建立相应的数学结构。4、模型求解：利用获取的数学资料对模型的所有参数做出计算。5、模型分析：对所得的结果进行数学上的分析