成都信息工程学院运筹学课程设计2013/11/23炼油厂生产计划安排2010062037 张靖2010062036 张磊2010062035 刘宏1 问题的提出炼油厂输入两种（原油1和原油2）.原油先进入蒸馏装置, 每桶原油经蒸馏后的产品份额见表1-29, 其中轻、中、重的分别为90、80和70.石脑油部分直接用于发动机油混合, 部分输入重整装置, 得辛烷值为115的重整汽油.1桶轻、中、重石脑油经重整后得到的重整汽油分别为0.6、0.52、0.45桶.蒸馏得到的和 , 一部分直接用于燃料油和燃料油的混合, 一部分经裂解装置得到和裂解油.裂解汽油的辛烷值为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：4：1.已知每天可供原油1为20000桶, 原油2为30000桶.蒸馏装置能力每天最大为45000桶, 重整装置每天最多重整10000桶石脑油, 裂解装置能力每天最大为8000桶.润滑油每天产量应在500至1000桶之间, 高档发动机油产量应不低于普通发动机油的40%.又知最终产品的利润（元/桶）分别为：高档发动机油700, 普通发动机油600, 燃料油400, 燃料油350, 润滑油150.试为该炼油厂制订一个使总盈利为最大的计划.2 建立模型1)定义变量原油1 X1中石脑油X2重油X3用于生产高档发动机油的轻石脑油 X4用于生产高档发动机油的中石脑油 X5用于生产高档发动机油的重石脑油 X6用于裂化生产的轻油X7用于裂化生产的重油X8用于生产润滑油的残油X9重整汽油X10用于生产普通发动机油的重整汽油 X11用于生产高档发动机油的裂化汽油 X12普通发动机油X13煤油X14原油2 X15重石脑油X16轻油X17用于生产普通发动机油的轻石脑油 X18用于生产普通发动机油的中石脑油 X19用于生产普通发动机油的重石脑油 X20用于生产煤油的轻油X21用于生产煤油的重油X22用于生产煤油的残油X23裂化油X24用于生产高档发动机油的重整汽油 X25用于生产燃料油的裂化油X26高档发动机油X27燃料油X28轻石脑油X29残油X30用于生产重整汽油的轻石脑油X31用于生产重整汽油的中石脑油X32用于生产重整汽油的重石脑油X33用于生产燃料油的轻油X34用于生产燃料油的重油X35用于生产燃料油的残油X36裂化汽油X37用于生产普通发动机油的裂化汽油 X38用于生产煤油的裂化油X39润滑油X402) 目标函数所求目标为利润最大化：max=0.7*X27+0.6*X13+0.4*X14+0.35*X28+0.15*X403) 约束条件X1<=20000;X15<=30000;X1+X15<=45000;X29-0.10*X1-0.15*X15<=0;X2-0.20*X1-0.25*X15<=0;X16-0.20*X1-0.18*X15<=0;X17-0.12*X1-0.08*X15<=0;X3-0.20*X1-0.19*X15<=0;X30-0.13*X1-0.12*X15<=0;X31+X32+X33<=10000;X10-0.6*X31-0.52*X32-0.45*X33<=0;X7+X8<=8000;X37-0.28*X7-0.20*X8<=0;X24-0.68*X7-0.75*X8<=0;-X40<=-500;X40<=1000;X27-X4-X5-X6-X25-X12<=0;X13-X18-X19-X20-X11-X38<=0;X14-X21-X22-X23-X39<=0;X28-X34-X35-X36-X26<=0;X34/X26=10/4;X26/X35=4/3;X35/X36=3/1;X40-0.5*X9<=0;X27-0.40*X13 = 0;84*X13-90*X18-80*X19-70*X20-115*X10-105*X37<=0; 94*X27-90*X4-80*X5-70*X6-115*X10-105*X37<=0;-X14+1.0*X21+0.6*X22+1.5*X39+0.05*X23<=0;X34/X26=10/4;X35/X36=3/1;X34-10*X36 = 0;X26-4*X36 = 0;X35-3*X36 = 0;X4 +X18+X31-X29<=0;X5+X19+X32-X2<=0;X6+X20+X33-X16<=0;X7+X21+X34-X17<=0;X8+X22+X35-X3<=0;X9+X23+X36-X30<=0;X11+X25-X10<=0;X38+X12-X37<=0;X26+X39-X24<=0;3 模型求解1)利用LINGO软件编程求解Max = 0.7 * X27 + 0.6 * X13 + 0.4 * X14 + 0.35 * X28 + 0.15 * X40;X1 <= 20000;X15 <= 30000;X1 + X15 <= 45000;X29 - 0.10 * X1 - 0.15 * X15 <= 0;X2 - 0.20 * X1 - 0.25 * X15 <= 0;X16 - 0.20 * X1 - 0.18 * X15 <= 0;X17 - 0.12 * X1 - 0.08 * X15 <= 0;X3 - 0.20 * X1 - 0.19 * X15 <= 0;X30 - 0.13 * X1 - 0.12 * X15<=0;X31 + X32 + X33 <= 10000;X10 - 0.6 * X31 - 0.52 * X32 - 0.45 * X33 <= 0;X7 + X8 <= 8000;X37 - 0.28 * X7 - 0.20 * X8 <= 0;X24 - 0.68 * X7 - 0.75 * X8 <= 0;-X40 <= -500;X40 <= 1000;X27 - X4 - X5 - X6 - X25 - X12 <= 0;X13 - X18 - X19 - X20 - X11 - X38 <= 0;X14 - X21 - X22 - X23 - X39 <= 0;X28 - X34 - X35 - X36 - X26 <= 0;X34 / X26 = 10 / 4;X26 / X35 = 4 / 3;X35 / X36 = 3 / 1;X40 - 0.5 * X9 <= 0;X27 - 0.40 * X13 = 0;84 * X13 – 90 * X18 – 80 * X19 – 70 * X20 – 115 * X10 – 105 * X37 <= 0;94 * X27 – 90 * X4 – 80 * X5 – 70 * X6 – 115 * X10 – 105 * X37 <= 0;-X14 + 1.0 * X21 + 0.6 * X22 + 1.5 * X39 + 0.05 * X23 <= 0;X34 / X26 = 10 / 4;X35 / X36 = 3 / 1;X34 – 10 * X36 = 0;X26 – 4 * X36 = 0;X35 – 3 * X36 = 0;X4 + X18 + X31 - X29 <= 0;X5 + X19 + X32 - X2 <= 0;X6 + X20 + X33 - X16 <= 0;X7 + X21 + X34 - X17 <= 0;X8 + X22 + X35 - X3 <= 0;X9 + X23 + X36 - X30 <= 0;X11 + X25 - X10 <= 0;X38 + X12 - X37 <= 0;X26 + X39 - X24 <= 0;2)数据求解结果Variable Value Reduced CostX27 7667.429 0.000000X13 19168.57 0.000000X14 15156.00 0.000000X28 0.000000 0.6238095E-01X40 500.0000 0.000000X1 15000.00 0.00000015 30000.00 0.000000X29 6000.000 0.000000X2 10500.00 0.000000X16 8400.000 0.000000X17 4200.000 0.000000X3 8700.000 0.000000X30 5550.000 0.000000X31 0.000000 0.2514286X32 0.000000 0.3017143X33 0.000000 0.3457143X10 0.000000 0.000000X7 4200.000 0.000000X8 3800.000 0.000000X37 1936.000 0.000000X24 5706.000 0.000000X4 1739.872 0.000000X5 2548.312 0.000000X6 2523.973 0.000000X25 0.000000 0.000000X12 855.2720 0.000000X18 4260.128 0.000000X19 7951.688 0.000000X20 5876.027 0.000000X11 0.000000 0.000000X38 1080.728 0.000000X21 0.000000 0.2228571E-01X22 4900.000 0.000000X23 4550.000 0.000000X39 5706.000 0.000000X34 0.000000 0.000000X35 0.000000 0.000000X36 0.000000 0.000000X26 0.000000 0.000000X9 1000.000 0.000000Row Slack or Surplus Dual Price1 23005.74 1.0000002 5000.000 0.0000003 0.000000 0.2539429E-014 0.000000 0.49696005 0.000000 0.62857146 0.000000 0.62857147 0.000000 0.62857148 0.000000 0.42228579 0.000000 0.400000010 0.000000 0.400000011 10000.00 0.00000012 0.000000 0.628571413 0.000000 0.2571429E-0114 0.000000 0.628571415 0.000000 0.400000016 0.000000 0.650000017 500.0000 0.00000018 0.000000 0.628571419 0.000000 0.628571420 0.000000 0.400000021 0.000000 0.412381022 0.000000 0.00000023 0.000000 0.00000024 0.000000 0.00000025 0.000000 0.800000026 0.000000 0.7142857E-0127 23988.49 0.00000028 19673.22 0.00000029 3429.500 0.00000030 0.000000 0.00000031 0.000000 0.00000032 0.000000 0.00000033 0.000000 0.00000034 0.000000 0.00000035 0.000000 0.628571436 0.000000 0.628571437 0.000000 0.628571438 0.000000 0.422285739 0.000000 0.400000040 0.000000 0.400000041 0.000000 0.628571442 0.000000 0.628571443 0.000000 0.40000003)最优生产方案炼油厂为了获得最大利润应计划原油1的使用量为15000桶、原油2的使用量为30000桶。