10 Blending Optimization
10 Blending Optimization
10 Blending Optimization
1
Technology Assessment
Conceptual
Process
2
Modeling Hierarchy
unit
operations
single single
plant Process
model
process Simulation
multi-plant refinery
model model
multiple
processes LP
multi-refinery
Simulation
model
3
What is “Linear Programming”?
z = a01
• Subject to the primary x1 + a02 x2 + ⋯ + a0N xN
constraints
x1 ≥ 0, x2 ≥ 0,
• and to M additional constraints (all
…b, xNpositive)
≥ 0,
n
ai1 x1 + ai 2 x2 + ⋯ + aiN xN ≤ bi
aj1 x1 + aj 2 x2 + ⋯ + ajN xN ≥ bj
ak1 x1 + ak 2 x2 + ⋯ + akN xN = bk
4
Linear Programming
5
Linear Programming Terminology
• Objective Function
» Function z to be maximized
• Feasible Vector
» A set of values x1, x2, …, xN that satisfies all constraints
• Optimal Feasible Vector
» Feasible vector that maximizes the objective function
6
Linear Programming Problem Solutions
7
Sample 2D LP Problem
• Objective function
Maximize z = x1 − 3x2
• Constraints
−2 x1 + x2 ≤ 4
0.5x1 + x2 ≤ 20
1.5x1 − x2 ≤ 10
3x1 + x2 ≤ 30
8
Sample 2D LP Problem
Maximize z = x1 − 3x2
−2 x1 + x2 ≤ 4
0.5x1 + x2 ≤ 20
1.5x1 − x2 ≤ 10
3x1 + x2 ≤ 30
9
Sample 2D LP Problem
Maximize z = x1 − 3x2
−2x1 + x2 ≤ 4
0.5x1 + x2 ≤ 20
1.5x1 − x2 ≤ 10
3x1 + x2 ≤ 30
10
Sample 2D LP Problem
Maximize z = x1 − 3x2
−2 x1 + x2 ≤ 4
0.5x1 + x2 ≤ 20
1.5x1 − x2 ≤ 10
3x1 + x2 ≤ 30
11
Sample 2D LP Problem
Maximize z = x1 − 3x2
−2 x1 + x2 ≤ 4
0.5x1 + x2 ≤ 20
1.5x1 − x2 ≤ 10
3x1 + x2 ≤ 30
12
Sample 2D LP Problem
Maximize z = x1 − 3x2
−2 x1 + x2 ≤ 4
0.5x1 + x2 ≤ 20
1.5x1 − x2 ≤ 10
x2
3x1 + x2 ≤ 30
x1
13
Sample 2D LP Problem
• Solution
» Z=6.67 @ x1=6.67 & x2=0.00
• Values at “corners”
x1 X2 z
5.2 14.4 -38
8.889 3.333 -1.111
0 4 -12
0 0 0
6.667 0 6.667
14
Sample LP Problem
• Objective function
Maximize z = x1 + x2 + 3x3 − 0.5x4
• Constraints
x1 + 2 x3 ≤ 740
2 x2 − 7 x4 ≤ 0
x2 − x3 + 2 x4 ≥ 0.5
x1 + x2 + x3 + x4 = 9
• Solution
» z= @ x1=0, x2=3.325, x3=4.725, & x4=0.95
15
Simplex Method
16
Restricted Normal Form
17
Restricted Normal Form
18
Restricted Normal Form
z ′ = − z1 − z2 − z3 − z4
= − (749.5 − 2 x1 − 4 x2 − 2 x3 + 4 x4 − y1 − y2 + y3 )
19
Restricted Normal Form
The Tableau
x1 x2 x3 x4 y1 y2 y3
z 0 1 1 3 -0.5 0 0 0
z1 740 -1 0 -2 0 -1 0 0
z2 0 0 -2 0 7 0 -1 0
z3 0.5 0 -1 1 -2 0 0 1
Z4 9 -1 -1 -1 -1 0 0 0
Z’ -749.5 2 4 2 -4 1 1 -1
20
Tableau Rows & Columns
• Rows are constraints
» Represent limits of operation
» Only enter optimal solution if actually limiting
» Material balances, utility balances, unit capacities, etc.
• Columns are variables
» An LP has more columns than rows
» Column Activity — value of variable
» Columns generally represent flows in the process
» Solution vector — column activities which maximize or minimize
objective function
21
Simplex Solution Method
22
Solved Form — Final Updated Tableau
x1 y2 y3
23
Interpretation of Solved Tableau
24
LP Software
25
Comparison of LP Models
Columns Rows
(variables) (equations)
26
Problems with LPs
• Nonlinear blending
» Linearize with blending factors
1.25 1.25
(RVP ) = ∑ v (RVP )
i i
γo = ∑ v i γ o, i = ∑V γ
i o, i
⇒ γ o ∑ Vi = ∑V γ
i o, i
∑V i
0= ∑V (γ
i o, i − γo )
27
Product Blending Example
Revenue $0
Cost $0
Profit $0.0
Profit per gal sold #DIV/0!
Raw Materials
Properties for Blending Calculations Cost & Availability Usage
Octane Cost Minimum Maximum
Number RVP VPBI ($/gal) Required Available Regular Premium Total
Butane 120 60 167.0 0.20 0 25,000 0 0 0
SR Gasoline/Naphtha 62 10.8 0.70 0 50,000 0 0 0
Reformate 100 RON 105 2 1.80 0 35,000 0 0 0
FCC Gasoline 87 5.1 0.70 0 50,000 0 0 0
Alkylate 92 4.6 1.05 0 12,000 0 0 0
Products
Lower & Upper Limits on Properties Price & Production Requirements
Price Minimum Maximum
Octane Number RVP VPBI ($/gal) Required Allowed
Regular 87 9.0 1.50 75,000
Premium 92 9.0 1.70 40,000
28