Please wait a minute...
Frontiers of Engineering Management

ISSN 2095-7513

ISSN 2096-0255(Online)

CN 10-1205/N

邮发代号 80-905

Frontiers of Engineering Management  2018, Vol. 5 Issue (2): 214-226   https://doi.org/10.15302/J-FEM-2017022
  本期目录
Online gasoline blending with EPA Complex Model for predicting emissions
Stefan JANAQI1(), Mériam CHÈBRE2, Guillaume PITOLLAT3
1. Ecole des Mines d’Alès, LGI2P, Parc G. Besse, 30035 Nîmes, France
2. TOTAL SA, 24, crs Michelet, 92069 PARIS LA DEFENSE Cedex, France
3. 3Manufacturing Engineering, BECKMAN COULTER, Georgia Institute of Technology, Marseille, France
 全文: PDF(169 KB)   HTML
Abstract

The empirical Complex Model developed by the US Environmental Protection Agency (EPA) is used by refiners to predict the toxic emissions of reformulated gasoline with respect to gasoline properties. The difficulty in implementing this model in the blending process stems from the implicit definition of Complex Model through a series of disjunctions assembled by the EPA in the form of spreadsheets. A major breakthrough in the refinery-based Complex Model implementation occurred in 2008 and 2010 through the use of generalized disjunctive and mixed-integer nonlinear programming (MINLP). Nevertheless, the execution time of these MINLP models remains prohibitively long to control emissions with our online gasoline blender. The first objective of this study is to present a new model that decreases the execution time of our online controller. The second objective is to consider toxic thresholds as hard constraints to be verified and search for blends that verify them. Our approach introduces a new way to write the Complex Model without any binary or integer variables. Sigmoid functions are used herein to approximate step functions until the measurement precision for each blend property is reached. By knowing this level of precision, we are able to propose an extremely good and differentiable approximation of the Complex Model. Next, a differentiable objective function is introduced to penalize emission values higher than the threshold emissions. Our optimization module has been implemented and tested with real data. The execution time never exceeded 1 s, which allows the online regulation of emissions the same way as other traditional properties of blended gasoline.

Key wordsemissions    reformulated gasoline    online control    global optimization
收稿日期: 2017-04-05      出版日期: 2018-06-28
Corresponding Author(s): Stefan JANAQI   
 引用本文:   
. [J]. Frontiers of Engineering Management, 2018, 5(2): 214-226.
Stefan JANAQI, Mériam CHÈBRE, Guillaume PITOLLAT. Online gasoline blending with EPA Complex Model for predicting emissions. Front. Eng, 2018, 5(2): 214-226.
 链接本文:  
https://academic.hep.com.cn/fem/CN/10.15302/J-FEM-2017022
https://academic.hep.com.cn/fem/CN/Y2018/V5/I2/214
VOCNOxTOX
xb1466.30134086.34
C11200130090
C21700140095
C31600130095
Tab.1  
SeasonRegionGasoline typeNOx penaltyVOC penaltyTOX penalty
Summer1RFG234
Tab.2  
OXYSULRVPE200E300AROOLEBENMTBETBETHTAM
0.00339.008.7041.0083.0032.009.201.530.000.000.000.00
Tab.3  
PropB1B2B3B4B5B6
MON52.380.479.4109.886.792.882.190
ROAD51.984.483.4114.988.1597.987.1100
RVP2.9913.8515.460.64153.19112.2115.37
V/L182.07134.4136.4216.24-24.7234.3120150
T90289348.4365.428534372100371
EP317444.4420.428535466100427
DI1,253.51,027.41,128.41,557- 561,8871001,230
E20016.664.448.4010000100
E30094.780.484.60.910000100
OLE3.121.430.40.101.72030
ARO8.922.423.499.9091.63040
BEN0.70.40.44000.403.7
SUL049.432.4084080
OXY00000000.1
01020002
214560100620
Tab.4  
Optimal values
F1F2F3F4
0.39222.36091.14870.1766
Optimal recipes
Basesu1u2u3u4
B10.00%0.00%3.32%2.80%
B20.00%0.00%54.65%65.05%
B387.25%79.04%26.94%19.67%
B411.57%6.61%13.02%12.48%
B51.18%1.83%2.08%0.00%
B60.00%12.52%0.00%0.00%
Emissions
BaselineE1E2E3E4
NOx: 1340.01392.621362.221314.791313.76
VOC: 1466.31479.451610.301480.561218.20
TOX: 86.3470.0673.0065.6965.01
Execution time [s]0.040.061450.62
Optimal gasolines
Propertiesx1x2x3x4
MON82.1083.0083.2283.16
ROAD87.1087.3887.1087.10
RVP15.3715.1915.1012.21
V/L142.56149.15141.79146.34
T90352.19354.85336.22342.17
EP400.19410.11404.45416.22
DI1,164.011,230.001108.541119.69
E20043.4140.0950.8651.88
E30075.1068.7672.0771.70
OLEF30.0026.5424.2520.00
AROM31.9736.5631.8431.89
BENZ0.380.400.360.37
SULF28.3626.2635.8938.51
OXY0.000.000.000.00
Tab.5  
xbxbRFGCFG
SummerWinter
OXY0.00.00.04.00.04.0
SUL3393380.0500.00.01,000.0
RVP8.711.56.410.06.411.0
E20041.050.030.070.030.070.0
E30083.083.070.0100.070.0100.0
ARO32.026.40.050.00.055.0
OLE9.211.90.025.00.030.0
BEN1.531.640.02.00.04.9
  
VOC+ TOXNOx
Normal emitters (wN)0.4440.738
High emitters (wH)0.5560.262
  
SummerWinter
Region 1Region 2Region 1Region 2
NOx1,340.01,340.01,540.01,540.0
VOC1,466.31,399.11,341.01,341.0
TOX86.3485.61120.55120.55
  
Property of xMinMax
SUL10.0450.0
OLE3.7719.0
ARO18.036.8
  
IF List
If Season= Winter, then y(RVP)=8.7.
If x(E300)>95 then y(E300)=95.
If x(SUL)<10 then y(SUL)=10, z(SUL)=x(SUL)–10.
If x(SUL)>450, then y(SUL)=450, z(SUL)=x(SUL)–450.
If x(OLE)<3.77 then y(OLE)=3.77, z(OLE)=x(OLE)–3.77.
If x(OLE)>19 then y(OLE)=19, z(OLE)=x(OLE)–19.
If x(ARO)<10 then y(ARO)=10, z(ARO)=–8.
If x(ARO)<18 then y(ARO)=18, z(ARO)= x(ARO)–18.
If x(ARO)>36.8 then y(ARO)=36.8, z(ARO)= x(ARO)–36.8.
  
1 Balas E (2010). Disjunctive programming. In: 50 Years of Integer Programming 1958–2008: From the Early Years to the State-of-the-Art, 283–340
2 Cason W W (1997). The illusory flexibility of the complex model: a graphical analysis of the implications of the complex model for refinery blending flexibility. In: NPRA Annual Meeting. San Antonio, Texas: NPRA Annual Meeting, 16–18
3 Chèbre M, Creff Y, Petit N (2010). Feedback control and optimization for the production of commercial fuels by blending. Journal of Process Control, 20(4): 441–451
https://doi.org/10.1016/j.jprocont.2010.01.008
4 EPA (2015). Complex emission model.
5 EPA (2016). Template spreadsheet for calculating emissions from combustion plants.
6 Furman K C, Androulakis I P (2008). A novel MINLP-based representation of the original complex model for predicting gasoline emissions. Computers & Chemical Engineering, 32(12): 2857–2876
https://doi.org/10.1016/j.compchemeng.2008.02.002
7 Hirshfield D S, Kolb J A (1994). Minimize the cost of producing reformulated gasoline by integrating EPA’s Complex Model into refinery linear programming (LP) models. NPRA Annual Meeting, 64–67
8 Hirshfield D S, Kolb J A (1997). The economics of gasoline reformulation: refining economics, emissions standards, and the Complex Model. In: NPRA Annual Meeting, San Antonio, USA
9 Jackson J, Vittachi K (1995). CITGO Corpus Christi refinery gasoline blender upgrade. In: NPRA Annual Meeting, San Francisco, USA
10 Janaqi S, Chebre M, Pitollat G (2015). Method and device for monitoring induced properties of a mixture of components, in particular emission properties, 69.
11 Korotney D J (1995). Reformulated gasoline effects on exhaust emissions: Phase III; investigation on the effects of sulfur, olefins, volatility, and aromatics and the interactions between olefins and volatility or sulfur. SAE Special Publications, 179–186
12 Misener R, Floudas C A (2009). Advances for the pooling problem: modeling, global optimization, and computational studies Survey. Applied and Computational Mathematics, 8(1): 3–22
13 Misener R, Gounaris C E, Floudas C A (2010). Mathematical modeling and global optimization of large-scale extended pooling problems with the (EPA) complex emissions constraints. Computers & Chemical Engineering, 34(9): 1432–1456
https://doi.org/10.1016/j.compchemeng.2010.02.014
14 Naman B T (1999). Linear models help refiners develop RFG [(reformulated gasoline)] recipes. Oil & Gas Journal, 97(8): 64–66
15 Nocca J L, Forestiere A, Cosyns J (1994). Diversify process strategies for reformulated gasoline. Fuel Reformulation, 4(5): 18–22
16 Raman R, Grossmann I E (1994). Modelling and computational techniques for logic based integer programming. Computers & Chemical Engineering, 18(7): 563–578
https://doi.org/10.1016/0098-1354(93)E0010-7
17 Treiber S (1998). RFG [(reformulated gasoline)]: the challenge to conventional blending technology. In: Gulf Publishing, eds. Hydrocarbon Processing 2nd International Conference on Process Optimization. Houston: Hydrocarbon Processing 2nd International Conference on Process Optimization, 101–103
18 Trierwiler L D (1995). Representing the simple and complex models in linear programming with respect to reformulated gasoline. In: American Energy Week ’95 “Pipelines, Terminals Storage, and Reformulated Fuels” Conference’ Houston, USA
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed