## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Brazilian Journal of Chemical Engineering

##
*Print version* ISSN 0104-6632*On-line version* ISSN 1678-4383

### Braz. J. Chem. Eng. vol.21 no.2 São Paulo Apr./June 2004

#### http://dx.doi.org/10.1590/S0104-66322004000200022

**PROCESS OPTIMIZATION**

**Flowsheet optimization of a lubricant base oil hydrotreatment process**

**J. L. Medeiros ^{I}; L. C. Barbosa^{I}; F. M. Vargas^{I}; O. Q. F. Araújo^{I}; R. M. F. Silva^{II} **

^{I}Escola de Química, Departamento de Engenharia Química, Universidade Federal do Rio de Janeiro, Phone (21) 2562-7419, Fax (21) 2562-7637, Ilha do Fundão, C.T. Bloco E, CEP 21949-909, Rio de Janeiro – RJ, Brasil. E-mail: jlm@h2oeq.ufrj.br

^{II}Centro de Pesquisas Leopoldo A. Miguez de Mello, Petróleo Brasileiro S.A, Phone (21) 3865-6824, E-mail: raissa@cenpes.petrobras.ufrj.br

**ABSTRACT**

Unsaturated, nitrogenated and sulfured compounds may reach undesirable levels in lubricant base oils, requiring hydrotreatment (HDT) at high temperatures and pressures. HDT processes are well known for their high capital and operational costs due to the use of hydrogen, compressors and multistage heterogeneous reactors. Process costs are thus highly dependent on the applied conditions. An oversized process entails unnecessary costs and capital investment. On the other hand, mild reaction conditions lead to unspecified products. This work analyzes the process using a compositional modeling and a kinetic framework developed previously (Barbosa et al., 2002) for the HDT of lubricant base oils. The industrial reactor is rigorously modeled on an adiabatic and multistage configuration. For the remaining equipment in the flowsheet, we adopt shortcut models for compressors and exchangers. The process is then optimized in terms of its variables subject to product specification constraints.

**Keywords:** lubricant oil, hydrotreatment, compositional modeling, optimization.

**INTRODUCTION **

Lubricant base oils are a complex mixture of paraffinic, aromatic and naphtenic hydrocarbons with molecular weights ranging from medium to high values, which produce oils with desirable viscosities, densities and distillation curves. However, these oils may have olefinic, nitrogenated and sulfured compound contents that jeopardize the quality of the final product. In order to meet composition specifications, these oils are submitted to the hydrotreatment (HDT) process before addition of other compounds for commercial production. This process – herein referred to as HDT-LUB – is characterized by the use of operational conditions of low or medium severity so that aromatic distribution is not seriously affected and sulfur compounds are not completely removed, as these compounds affect positively resistance to oxidation under service.

In addition, process costs have nontrivial characteristics associated with: (i) the high value of product and raw material; (ii) mandatory compliance with the product oil composition requirements; (iii) the use of an expensive reactant (hydrogen); (iv) operation at high temperatures and pressures; (v) need for investing heavily in compressors and specific heterogeneous catalysts.

This discussion is important in establishing that observation of operational conditions near an optimal point is a critical factor in HDT-LUB profitability. Problems in obtaining an optimum process include the description of: (i) composition of base oil which, like other crude oil fractions, has a complex structure; (ii) reactive process involving heterogeneous catalytic reaction networks; (iii) adiabatical operation of the industrial reactor in various hydrogenation stages, under medium and high pressure conditions, multiple hydrogen feeds and with liquid and vapor phases in constant equilibrium; and (iv) the overall process including, in addition to the reactor, peripheral structures for feed preheating, heat recovery, compression and hydrogen recycling.

This paper analyzes HDT-LUB, aiming to optimize the process and operational aspects under the constraints of product specifications. Previously developed compositional modeling (Barbosa et al., 2002), based on feed and product characterization data in a pilot plant (Valle et al., 1998), is used for the base oil and reaction network.

HDT-LUB is approached on three levels:

(i) development of a steady state simulator to predict the final state of the product and hydrogen and utilities consumption from process variables; (ii) definition of the economic objective, taking into account realistic aspects of process costs and the need to achieve product specification; and (iii) optimization of design and operation.

The process was modeled rigorously on an industrial reactor, formed of several adiabatic stages with liquid and vapor phases in constant equilibrium. Adiabatic separation vessels and the recycling structure were also modeled rigorously. The rest of the process was modeled using simple compressor, heat exchanger, scrubber and furnace models (Douglas, 1988). Investment was modeled according to Guthrie's proposition (Douglas, 1988).

**THE HDT-LUB PROCESS**

The HDT process is shown in Figure 1. A fixed base oil feed of 25000 kg/h, pumped and mixed with hydrogen under pressure, is assumed. The feed is preheated by a reactor outlet stream and then is heated to the reactor inlet temperature in the furnace. HDT takes place in a fixed bed with multiple gas recycle injections to supply hydrogen and control reaction temperature. After leaving the reactor and losing heat in the next water-refrigerated heat exchanger, the mixture is depressurized and separated from hydrogen in a high-pressure separator. The gas generated is sent to a water scrubber for hydrocarbon, H_{2}S and NH_{3} retention. The recycle compressor pumps the recycle gas obtained at the top of the scrubber. The oil at the bottom of the high-pressure separator is expanded in the low-pressure separator, producing residual gas and treated oil. Key variables, which have decisive roles in the process performance, are the following:

1) Temperature at the furnace exit (T_{FURNACE}). Its increase favors high rates of HDT reaction, although one has to consider the impact on furnace costs and the highly exothermic characteristics of HDT reactions, which increase the possibility of generating elevated thermal levels at some point in the reactor. Also, some exothermic equilibrium reactions, such as hydrogenation of tri- and tetra-aromatics, are shifted unfavorably with an increase in reactor temperatures.

2) Feed pressure at reaction stage 1 (P_{COMPRESSOR}). Hydrogen fugacity directly influences the kinetics and severity of the process. Its increase favors all HDT reactions, although it creates investment problems due to feed and recycle compressor design.

3) Reaction stage spatial velocities (WHSV_{1}, WHSV_{2} and WHSV_{3}). An increase in WHSV (weight hourly spatial velocity) translates into reduction of catalyst mass, reducing the severity of the process and investments. WHSV reduction can be used to maintain the level of severity as feed pressure is reduced.

4) Hydrogen/feed ratios RH2HC_{1}, RH2HC_{2} and RH2HC_{3} in reaction stages. The increase in hydrogen recycling has advantages, such as the quenching effect at injection points, and favors heavy hydrocarbon vaporization, thus helping to maintain thermal control over the reactor. Also, a higher concentration of hydrogen inside the reactor results in easier maintenace of fugacity level, guaranteeing constant reaction rates. Problems that arise are the need for larger compressors and equipment in the recycle loop.

** COMPOSITIONAL MODEL AND REACTION NETWORK **

Barbosa et al. (2002) characterized lube base oil by data reconciliation, using a compositional model (CM), based on experimental data obtained from Valle et al. (1998). Adjustment of CM parameters is based on making the values of estimated properties equal to experimental data, such as density and distillation curve. These calculations use correlations and thermodynamic models, as well as the Joback group contribution method (Reid et al., 1987). The CM developed by the authors has 23 key components with five different paraffinic side chains. The authors also proposed a reaction network for the HDT-LUB process, with 14 (ncin=14) kinetic laws, 39 (n=39) chemical reactions and 51 (nc=51) chemical species. Kinetic and adsorption parameters for the network were estimated based on pilot plant data obtained by Valle et al. (1998). All results generated by Barbosa et al. (2002) were used in the present paper: (i) CM feed parameters; (ii) the reaction network matrix (: stoichiometric matrix (nc by nr); : component selection matrix for reaction rate definition (nr by nc); : matrix for association of kinetic laws with reactions (nr by ncin); (iii) vector of Arrhenius constants (__K__(T)) ((gmol/h)/(kg catalyst . bar)) for kinetic laws (ncin by 1); and (iv) Langmuir constants for adsorption on the catalyst ().

**HDT-LUB MODELING**

The flowsheet shown in Figure 1 was modeled as follows. After solving the equipment equations, their installed costs (ISBL, US$) were determined by the Guthrie models (Douglas, 1988) for year 1988 using the corresponding Marshall-Swift (M&S) index. Annual operational costs (UTIL, US$/year) were defined based on respective consumption of utilities. Table 1 contains the premises adopted:

1) Thermodynamic properties:

2) Fugacities, densities and enthalpies are estimated using cubic equations of state (Soave/Peng-Robinson) and the Joback method (critical constants and ideal gas enthalpies).

3) Liquid-vapor separation via Flash Rigorous Algorithms (P,H) and Flash (P,T): applied to separation vessels, reactor and energy balance in two-phase streams.

4) Heat exchangers: once temperature approximation intervals (DT_{AP}) are defined, thermal duty for each stream is calculated (via Flash (P,T)). Limiting stream, heat exchanger duty and final temperatures are determined. Area and ISBL are calculated for constant thermal coefficients (U). Operational costs (UTIL) are calculated only for water-refrigerated exchangers.

5) Reactor: The three adiabatic stages model assumes: (i) liquid-vapor equilibrium in the axial direction; (ii) plug-flow; (iii) stages calculated individually by integration of mass balance equations (Equation 1) associated with isenthalpic liquid-vapor equilibrium for determining temperature profiles.

where

F_{0} is the oil feed flow rate (kg/h), t is the spatial time (kg catalyst/(kg/h)) and N is the species rate vector (gmol/h). Equation 1 is integrated into stage k from 0 to 1/WHSVk. Reactor ISBL is determined by the catalyst inventory (catalyst price in Table 1) and the pressure vessel ISBL.

Compressor: Power (BHP) in N_{S} stages, efficiency h, feed (T_{E} (K)) and interstage (T_{I }(K)) temperatures, compression ratio Q, flow rate Q (gmol/s) and polytropic exponent (g = C_{P}/C_{V}); Equation 2 follows (R=8.314 J/gmol.K):

ISBL is determined with BHP, stipulating a turbine-operated centrifuge compressor. Operational cost corresponds to high-pressure vapor in the turbine.

The process design is generated from eight-variable process vector __V__: (i) P_{COMPRESSOR} (bar); (ii) T_{FURNACE} (ºC); (iii) WHSV1, WHSV2 and WHSV3 (kg/h/kg catalyst); and (iv) RH2HC_{1}, RH2HC_{2} and RH2HC_{3} (Nl/kg). Search intervals are defined by __V___{min} and __V___{max} so that 100 __<__ RH2HC1 __<__ 400; 50 __<__ RH2HC_{2}, RH2HC_{3 } __<__ 400; 1 __<__ WHSV_{1}, WHSV_{2}, WHSV_{3} __<__ 6; 300 __<__ T_{FURNACE < } 350; and 40 __<__ P_{COMPRESSOR} __<__ 120).

**ECONOMIC OBJECTIVE **

It can be shown that a ten-year process with 12% annual rate of return, negligible labor costs and "instantaneous" implementation, has an annualized profit function (AP, US$/year) given approximately by Equation 3 (Douglas, 1998):

where ISBL_{k} and UTIL_{k} are installed capital (US$) and annual utilities costs (US$/year) of unit k and REV and RAW represent annual income and raw material costs (US$/year). In HDT-LUB these costs are derived from H_{2} and fixed-price oil consumption (P_{H2} and P_{RAW}). Income results from product rates and its market price (P_{PROD}), which has a maximum (P_{LUB}) if specifications are achieved. These constants are shown in Table 1. Equation 4 gives the product price model, where C_{k} is the k^{th} valid assay for the commercial fitness decision, while [L_{k}, U_{k}] represent the acceptable interval for C_{k}, and b is a large positive constant.

This function has the role of eliminating specification constraints from the objective function of the process, so that the values of products that do not agree with * L_{k} < C_{k} < U_{k}* for any k depreciate. Specifications used for characterizing commercial oil are listed in Table 2. Equation 5 gives the optimization formulation. Search variables in the unconstrained space q are defined in Equation 6. The HDT-LUB process was optimized using the SIMPLEX Nelder-Mead method (Reiklaitis et al., 1983), with Equations 5 and 6, in MATLAB R12 environment. Final values obtained were P

_{COMPRESSOR }= 40.8 bar, T

_{FURNACE }= 328.6 ºC,

__WHSV__= [1.66, 3.2, 3.7] h

^{–1}, and

__RH2HC__=[177.8, 213.9, 130] (Nl/kg) .

**RESULTS**

Starting from a viable loosely defined initial condition, trajectory of function AP (Figure 2a) evolves along the feasible region contours, exploring viable and nonviable designs (Figures 2a and 2c). Each nonviable condition found creates a sudden reduction in AP, causing the method to increase severity and thus return to a viable design. The objective value was markedly increased from the initial condition, tightening process operational conditions and obtaining a more economic combination of variables that influence severity. When the search was interrupted near 14 MMUS$/year, it was sloping up, with 2.6 MMU$ of total ISBL and 3.8 MMU$/year for H_{2} costs. The search is directed towards lowering operating pressure and reactor size (Figures 2b and 2d), increasing furnace temperature to maintain the necessary severity for meeting the specification. It can be seen that 40ºC viscosity fitness was the main factor that defined the applied severity level (Figure 2f). Results for the optimum design simulation showing some property profiles for the reactor are given in Figure 3. Some reductions such as the following can be highlighted: (i) density; (ii) sulfur and nitrogen contents; (iii) simulated distillation curve temperatures and (iv) effects of the three hydrogen injections on temperature profiles, enthalpy fluxes and fraction vaporized inside the reactor.

** CONCLUSIONS**

A methodology for the lube base oil HDT modeling and optimization process, based on the following: (i) compositional feed modeling; (ii) a preliminary kinetic survey of the extensive reaction network; and (iii) a set of quality aims that must be met. It was shown that an efficient compositional model, together with a rigorous reactor description and simple models for peripheral equipment, is a practical tool for project analysis near optimum conditions.

**ACKNOWLEDGMENTS **

The authors are grateful to the National Petroleum Agency - ANP - for the scholarship received and to FINEP/CTPETRO and CNPq for the financial support.

**REFERENCES **

Barbosa, L.C.; Vargas, F.M., Medeiros, J.L., Araújo, O.Q.F. and Silva, R.M.F., Modelagem de Hidrotratamento de Óleos Bases de Lubrificantes, presented at the XIV Congresso Brasileiro de Engenharia Química, 2002. [ Links ]

Douglas, J.M., Conceptual Design of Chemical Processes, McGraw-Hill Book Company, 1988. [ Links ]

Reid, R.C., Prausnitz, J.M., Poling, B.E., The Properties of Gases and Liquids, McGraw-Hill, N.Y., 1987. [ Links ]

Reiklaitis, G.V., Ravindran, A. and Ragsdell, K.M., Engineering Optimization, J. Wiley & Sons, 1983. [ Links ]

Valle, A.A.F., Rego, P.S.S., Carvalho, R.L., Produção de óleo mineral isolante por hidrotratamento catalítico na PETROBRAS. Proceedings of the Rio Oil & Gas Conference, 1998, Rio de Janeiro, Brazil, IBP03098. [ Links ]

Received: February 13, 2003;

Accepted: October 4, 2003