Next Article in Journal
Comprehensive Review on Smart Techniques for Estimation of State of Health for Battery Management System Application
Previous Article in Journal
Health Assessment and Remaining Useful Life Prediction of Wind Turbine High-Speed Shaft Bearings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Hydrodearomatisation Reaction Kinetics on the Modelling of Sulphur and Aromatics Removal from Diesel Fuel in an Industrial Hydrotreating Process

by
Sandra B. Glišić
* and
Aleksandar M. Orlović
Faculty of Technology and Metallurgy, University of Belgrade, Karnegijeva 4, 11120 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Energies 2021, 14(15), 4616; https://doi.org/10.3390/en14154616
Submission received: 9 June 2021 / Revised: 12 July 2021 / Accepted: 19 July 2021 / Published: 30 July 2021
(This article belongs to the Special Issue Emissions from Internal Combustion Engines)

Abstract

:
Over the years, the hydrotreating process has been considerably improved in order to facilitate the production of environmentally friendly diesel fuels by reducing sulphur and aromatics concentrations, as mandated by contemporary emissions regulations. In this study, different kinetic models for the hydrodearomatisation reaction and the influence of reaction rate on performance of the industrial trickle bed reactor for hydrotreating of gas oil and light cycle oil fractions were analysed. The impact on reactor temperature, catalyst wetting efficiency, and conversion of sulphur and aromatics were determined. The results of simulations were compared with experimental data from an industrial test run and the best model for the observed process is proposed. Reactor performance and overall efficiency of the process is strongly dependent on the kinetics of hydrodearomatisaton with respect to aromatics conversion but even more so with respect to the temperature increase in the reactor, which affects all key catalytic reaction parameters, catalyst wetting efficiency, and thus the sulphur conversion. Based on the obtained simulation results, it could be concluded that reactor performance is strongly dependent on the hydrodearomatisation reaction. The best predictions of outlet temperature as well as sulphur and aromatic conversion (deviation from the experimental value 0.87 K, 0.01% and 2.57%, respectively) are achieved with the Langmuir–Hinshelwood kinetic model proposed by Owusu-Boakye.

Graphical Abstract

1. Introduction

Impurities like sulphur and aromatics contained in refined petroleum products lower the quality of fuel and have an important impact on soot emissions with serious effects on human health and potential carcinogenic effects [1,2,3,4]. Several studies have shown that undesired soot emissions from diesel fuel are proportional to the aromatics content of the fuel [2,3,4,5].
The current sulphur concentration limits are very low, and consequently the process of desulphurization is highly optimised to reduce sulphur content, but elimination of aromatic compounds is still difficult, especially for diesel feedstocks with high aromatics content. Nowadays, more attention is paid to the content of aromatics in the product, and aromatics specifications for diesel fuel vary widely (depending on the location and feedstock composition). According to the regulatory standard for fuel specifications, the aromatic hydrocarbon content should not exceed 10 vol%, while at the same time, for small refineries, the reference fuel specification limit is 20 vol% (aromatic hydrocarbon content ≤21.0 wt% and polycyclic aromatic hydrocarbon content ≤3.5 wt% by the D5186-96 test method). The saturation and content of aromatic compounds directly impact the cetane number and the fuel emissions characteristics [3].
Stability and quality of the fuels is improved in the process of hydrotreating as it is the most important process in refineries for sulphur removal and aromatics reduction through the reactions of hydrodesulphurisation (HDS) and hydrodearomatisaton (HDA). In the literature, a broad range of different research papers analysing the hydrotreating process can be found [6,7,8,9,10,11,12] in which hydrodesulphurisation and hydrodearomatisation reactions in trickle-bed reactors are the key step [6,7]. Study and further improvement of hydrotreating processes could be analysed through adequate investigations of kinetics and reactor modelling. The previously developed mathematical model used in this study was developed with input data for complex mixtures of hydrocarbons being processed in industrial units and not from studies of bench-scale and pilot-scale processes, which are often found to be unreliable for industrial-scale applications. A reliable and well-established mathematical model can be crucial for process design, determination of optimum operating parameters, for the understanding of the entire process, and for a realistic representation of the process [6,7,8]. Complexity of the mathematical model implies: detailed analysis of the inlet and outlet mixtures, well-established material and energy balances, the adequate kinetic equations and parameters, the phase equilibrium and thermodynamic calculations for multicomponent mixtures, the rate equations of reactions that are taking place in very complex reaction mixtures in liquid and vapour phases, the evaporation rate as well as the effectiveness factors, and catalyst wetting efficiency. In such a complex reaction network, the reactivity of components is affected by many factors [6,7,9].
Comparisons of the reactivity of different aromatic compounds have been reported in a review study [10]. Most of the quantitative reactivity data characterise reactions of individual aromatic hydrocarbons with hydrogen and few studies have been reported with mixtures of hydrocarbons chosen so that inhibition effects were negligible. However, inhibition effects due to competitive adsorption can affect reactivity significantly. Several review studies about HDA reaction networks and kinetics can be found in the literature [11,12,13].
Hydrodearomatisation reactions are reactions in which aromatic rings are saturated with hydrogen. For polycondensed aromatic hydrocarbons, the hydrogenation of the first ring is in general the fastest step, and the rate of hydrogenation for subsequent rings tend to be lower-rate steps, with the last ring being the least reactive [12,13]. The rate of hydrogenation of the last ring is significantly lower than that of the first one. Overall, aromatic conversion is therefore often limited by thermodynamic equilibrium, since inside the liquid-filled pores of the catalyst particles diffusion of hydrogen from the gas phase through a liquid film and into the particle has to occur before a reaction can take place. Phenanthrene is a good model compound to represent polyaromatics. An important aspect of the hydrodearomatisation reactions is that in typical hydrotreating conditions the conversion can be limited by the thermodynamic equilibrium and by mass transport of the reacting species depending on the process conditions.
There have been several attempts toward a thorough description of the HDA reaction, which can be found in the literature [14,15,16,17,18,19,20,21,22,23,24,25,26]. Most of the HDA kinetic models account for one sum reaction and very few are developed to accommodate three classes of reacting aromatic compounds, mono-, di-, and triaromatics. Kinetic models for the HDA reaction in a gas/oil system mainly assume that hydrogenation and dehydrogenation reactions occur according to the Langmuir–Hinshelwood mechanisms and the HDA reaction is represented as a first-order reversible reaction. Cheng et al. [27] investigated the performance of a fixed-bed reactor in concurrent and counter-current flows to remove sulphur and aromatics in diesel fuel. The model presented by this group is a one-dimensional heterogeneous model that accounts for the HDS and HDA reactions to simulate the concentration profiles of the reactants and products in the gas, liquid, and solid phases.
A complex and reliable mathematical model for hydrotreating of straight run gas oil (SRGO) blended with fluid catalytic cracking naphtha and light cycle oil (FCC-N-LCO) in an industrial trickle bed reactor, where several types of sulphur and aromatic compounds coexist and several classes of reactions occur simultaneously, was previously developed and verified [6,7]. In this study, a previously developed model [6] was used for the investigation of different kinetic models for the HDA reaction. Different kinetic models for the HDA reaction were incorporated in the model [6] and the results of the simulations were analysed in order to obtain the most accurate prediction of reactor performance. Concentrations of sulphur and aromatic compounds in products along with temperature difference in the reactor were verified by comparison with industrial test run data.

2. Experimental Part

2.1. Industrial Test Run Data

The industrial test run was carried out in a catalytic hydrotreater reactor, which is normally a steady-state operating tubular flow reactor in which an inlet stream of a straight run gas oil stream was blended with a stream of light cycle oil (LCO) and fluid catalytic cracking naphtha (FCC N) upstream of the reactor. The blended inlet stream was hydrotreated in the presence of a commercial Co-Mo/γ-Al2O3 catalyst. The inlet stream was analysed using GC-MS and different classes of sulphur compounds were found to be present in the stream: substituted benzothiophenes (C1–C4-BTs), dibenzothiophene and naphthothiophene (DBT and NT), and substituted DBTs and NTs (C1-C3-DBTs and -NTs). Properties of the feedstock and output stream used in the test run relevant for the model development and validation are given in the Table 1. All chemical characterisations and measurements were obtained using methods specified by the relevant technical standards, Table 1 in [7].
The composition of the initial feed into the catalytic reactor for this experiment was 87 vol% SRGO and 13 vol% FCC N-LCO. Pressure in the reactor was 40 bar, inlet temperature 606 K, and amount of the catalyst 18,000 kg. A more detailed description of the reactor, catalyst, and the test run can be found in the literature [6,7]. Experimental points used for the model validation in this study were: overall conversion of sulphur of 99.49 wt%, overall conversion of aromatics of 70 wt%, and the reactor temperature difference of 11.5 K.

2.2. The Reactor Model

The developed mathematical model used for the simulation of the industrial hydrotreater is deterministic and one-dimensional. The model takes into consideration the volatility of the reaction mixture and therefore the HDS and HDA reactions in the liquid and in the vapour phases. The composition of the reaction mixture was approximated using pseudocomponents according to the distillation curve of the reaction mixture. Pseudocomponents for thermodynamic representation in the vapour liquid equilibrium of different classes of sulphur compounds (BT, DBT1, DBT2, DBT3) were: n-dodecane, n-hexadecane, and n-octadecane, respectively. Aromatics were represented by toluene, tetraline, and phenanthrene as pseudo components. Two additional components in this model were methane and hydrogen as the main constituents of the gaseous phase. Toluene represents the monoaromatics and tetraline stands for diaromatics typically present in diesel fractions. Phenantrene was included to provide information on the successive hydrogenation steps of polyaromatic compounds.
The complex reaction network is even more complicated when vapour–liquid equilibrium (VLE) is included. VLE is crucial in the hydrotreater simulation due to the volatile nature of the oil and high solubility of hydrogen in the oil [6,9]. The presence of reactants in both phases creates a complex reaction system. In this study, vapour–liquid equilibria through the catalyst bed were simulated using the Peng–Robinson equation of state (PR-EOS) with thermodynamic properties and components distribution calculated by UniSim and AspenPlus software. The model assumes trickle flow in the adiabatic catalytic reactor with negligible axial dispersion, while catalyst deactivation during the industrial test run was neglected. Material balance equations in the trickle-bed reactor (TBR) for all reacting components in vapour and in liquid phases are described by the set of ordinary differential equations (ODEs). The material and overall energy balance equations, calculation of catalyst wetting efficiency, and overall catalyst effectiveness factors for vapour and liquid phases are incorporated in the model and described in a previously published paper [6]. Below are the summarised main equations used in the model (Table 2). All equations and systems of equations were solved in a program developed in MATLAB® & Simulink® Release 2010b.

2.3. Hydrodearomatisation Reactions

Kinetic equations for different classes of sulphur compounds involved in the HDS reaction network are represented by the Langmuir–Hinshelwood rate equations proposed in the literature [19,20,21,22]. These equations were developed on a commercial CoMo/Al2O3 catalyst under operating conditions significant to industrial applications. The kinetic equations for hydrogenation of aromatics used in this study are summarized within Table 3.
Hydrogenation of aromatic rings is reversible, highly exothermic, and a very important reaction in hydroprocessing, and therefore development of HDA reaction rate equations requires great attention. Several researchers have proposed models based on a simple first-order reversible reaction. Typical industrial feeds contain mixtures of tri-, di-, and monoaromatics, and hydrogenation proceeds via consecutive reversible reactions; however. such an approach has certain limitations. The HDA reaction is equilibrium-limited at high temperatures, under which the reverse reaction of naphthalene dehydrogenation occurs. Since HDA reactions are reversible, therefore, an equilibrium conversion is the upper limit for conversion. Cheng and co-authors [27] found that 380 °C was the upper temperature limit that can be practically employed.
Kinetic expressions used in this study are proposed in Table 3 and briefly described below.
Chowdhury [14] experimentally investigated desulphurisation and hydrogenation of aromatics in diesel oil in an isothermal lab-scale TBR where three groups of aromatics were considered: mono-, di-, and polyaromatics. All polyaromatics (tri-, tetra-, and penta-aromatics) behave like triaromatics. In this model, HDA reactions are first-order reversible with the assumption of a completely wetted catalyst bed. The heat of HDA reactions is set at 67 kJ/mol. The temperature dependencies of all the HDA reaction rate constants have been described by the Arrhenius law, where activation energies and the frequency factors have been determined via regression analysis.
Avraam [19] developed a model for trickle-bed reactors for the hydroprocessing of light oil feedstocks containing volatile compounds where reactions of desulphurisation, denitrogenation, saturation of olefins, and hydrogenation of aromatics were considered. Hydrogenation of aromatics was treated as a reversible reaction where A1 is hydrogenated into a totally or partially saturated compound, A2, with a substantial heat of reaction of 180 kJ/mol. This model was validated using pilot-scale data.
In the work of Owusu-Boakye [22], experiments in a lab-scale TBR for single- and two- stage hydrogenation of aromatic compounds in light gas oil were performed. Using the single-site mechanism form of the Langmuir–Hinshelwood rate of reaction, a kinetic model for the HDA reaction was developed. The apparent kinetic parameters of the single- and two-stage processes were obtained using the nonlinear least squares approach. Apparent activation energy and heats of adsorption for calculation of kinetic constants were determined directly from the slopes of the Arrhenius and Van’t Hoff plots. Heat of reaction used in this model is 85 kJ/mol.
Yui [23] developed a simple power law kinetic model for HDA reactions based on first-order reversible reactions. Kinetic parameters are calculated using published catalytic aromatics hydrogenation data from a variety of sources. Both forward- and reverse-rate constants are calculated, and heat of reaction used is 255 kJ/mol.
De Oliveira [24] brought a novel approach for kinetic modelling of hydrotreating reactions. In this study, hydrotreating of LCO gas oil was simulated using a stochastic simulation method. In view of the fact that the HDA is a monomolecular reaction, the stochastic rate constants of forward and reverse HDA reactions were calculated using quantitative structure/reactivity correlations (QS/RCs). Heat of HDA used in the de Oliveira model is 51 kJ/mol.
A dynamic mathematical model for the commercial hydrotreating process, including hydrodesulphurisation, hydrodenitrogenation, and hydrodearomatisation, was developed by Liu and co-workers [25]. The hydrogenation of each aromatic group (mono-, di-, polyaromatics) is presented as power-low reaction kinetics. The reaction rate constants are calculated according to the Arrhenius law and heat of reaction used is 125 kJ/mol.
Values of calculated kinetic parameters for the inlet temperature in a reactor are given in Table 4.
Simulations using different kinetic models were carried out in a user-friendly program developed in MATLAB and Simulink Release 2010b. The system of ordinary differential equations was solved using a fourth-order Runge–Kutta routine with variable step size. Results of simulations were compared mutually and with experimental data.

3. Results and Discussion

As the result of simulations, molar flows of different classes of sulphur and aromatic components in liquid and in vapor phases were calculated, as well as temperature changes through the reactor, catalyst wetting efficiency, and effectiveness factors. All variables were calculated through the reactor as a function of catalyst weight. Key results to be observed were overall conversions of aromatics and sulphur compounds, since the final goal of the hydrotreating process is the reduction of aromatics and removal of sulphur in the final product. Figure 1 shows the conversions of aromatics calculated using different HDA kinetic models.
Experimental data of total aromatics conversion from the industrial test run was compared to all proposed kinetic models. The general conclusion is that there is a strong difference in the results of simulations when Langmuir–Hinshelwood kinetic models are used when compared to the power law kinetic models.
In Figure 1, conversions of total aromatics can be compared. The reaction rate of HDA is higher as the number of aromatic rings increases (comparing kinetic constants in Models 1 and 6, Table 4). Therefore, the HDA reaction of monoaromatic molecules was the slowest, the rate of hydrogenation of diaromatic molecules was faster, and tri- and polyaromatic molecules were the most reactive compounds. In addition, many studies of HDA reactions have shown that no partially hydrogenated ring compounds were observed in the product analyses, which proves that HDA occurs in a ring-by-ring manner [8,28].
The best prediction of total aromatics conversion in this study was achieved using the Langmuir–Hinshelwood kinetic model for HDA reactions proposed by Owusu Boakye (Figure 1). The kinetic model proposed by Chowdhury showed the highest deviation from the experimental point due to the lowest energy activation and consequently the low temperature increase in the reactor (Table 4 and Figure 1 and Figure 2). The power law kinetic model proposed by Yui, as shown in Figure 1, resulted in considerably higher deviation from experimental point.
Since HDS and HDA reactions occur simultaneously and in parallel, in any hydrotreater reactor it is very important to monitor both sulphur and aromatics conversion, as well as temperature increase in the reactor (impacted by both exothermic reactions). Overall sulphur conversion along the catalyst bed is shown in Figure 2. Almost all of the applied HDA kinetic models predicted correctly the reduction of sulphur content in the reaction mixture along the catalyst bed. It can be concluded that all HDA kinetic models except those proposed by Liu and Avraam match very well with the experimental data for sulphur conversion (Figure 2a).
The temperature increase within the catalyst bed is shown in Figure 2b. Since the reactions occurring in the reactor (HDA and HDS) are exothermic, a temperature increase is to be expected. The degree of the temperature increase is very dependent on the applied kinetic model. The outlet temperature ranged from 611.7 K to 627.1 K and the best agreement with the experimental data was achieved using the Owusu-Boakye model.
The temperature increase and deviation from the experimental value was the highest for the Avraam and Yui kinetic models, which correspond to the highest proposed value of reaction heat for HDA reaction heat within the kinetic model. The calculated temperature increase was lowest for the de Oliveira kinetic model. Temperature change in the reactor had crucial impact on conversion and also on vapour–liquid equilibrium since equilibrium constants and liquid-to-vapour flow ratios are polynomial functions dependent on the temperature. At lower temperatures, more liquid is present in the reactor; as the temperature increases through the catalyst bed, liquid evaporates until catalyst particles are dried out.
The difference in catalyst wetting efficiency between the predictions of the models, using different kinetic equations for the HDA reaction, is very significant (Figure 3). Wetting of the catalyst bed is calculated depending on the reactor design and also the catalyst particle shape and size, but mostly on the liquid flow rate and velocity. The degree of liquid evaporation is dependent on the degree of temperature increase, while temperature increase is dependent on reaction rate and heat released by chemical reactions. The lowest wetting of the catalyst bed was predicted by simulations with the Yui and Avraam HDA models, where liquid evaporates almost instantly. For these models the temperature increase was the highest. The de Oliveira model predicted the lowest temperature increase (Figure 2), therefore wetting efficiency was very high with a substantial quantity of liquid through the entire reactor length.
Experimental data for the wetting efficiency in the simulated industrial reactor during the industrial test run were not available. It was formerly reported that very poor wetting efficiency is to be expected in the industrial-scale reactors that operate at low LHSV [29]. A different study showed that the degree of the wetting efficiency of the catalyst bed affects the reaction rates greatly [30] and decreasing the amount of the liquid phase increases reaction rates and an intensification of the whole process can be achieved [31,32,33,34,35].
The thermodynamic model predicted catalyst bed wetting efficiency in the reactor, taking into account the flows of vapour and liquid phase, which are calculated through VLE calculations based on the PR-EOS. At the beginning of the process, the catalyst bed was partially wetted, while along the catalyst bed the evaporation of the liquid phase occurred. As this model is based on the inlet mixture that contains considerable amounts of volatile compounds, this phenomenon is likely and expected to occur (high methane content in the inlet hydrogen stream).
Overall effectiveness factors in the reactor for all models in vapour and in liquid phases are shown in Figure 4. By calculation of effectiveness factors, the impact of mass transfer limitations of the process is included in the model.
Overall effectiveness factors decreased along the reactor due to the temperature increase, and therefore reaction rates increased, resulting in a growing influence of the mass transfer resistance. In liquid phase, values of effectiveness factors were approaching zero due to the diminishing amounts of liquid in the catalyst bed, while in the gas phase effectiveness factors were higher. The values of calculated effectiveness factors are in agreement with the literature data [28,29].
Experimental data obtained in the industrial reactor were: total sulphur conversion of 99.49 wt%, conversion of aromatic compounds of 70 wt%, and temperature increase in the reactor of 11.5 K. Shown in Table 5 are the results of the total sulphur conversion, temperature increase in the reactor, and conversion of aromatics obtained by the six models and the comparison with experimental results.
Values of overall sulphur conversions in Table 5 are slightly higher or lower than the experimental values. Model 6 calculated the lowest and Model 4 the highest values of sulphur conversion. Apart from the reaction kinetics, overall sulphur conversions were following the values of effectiveness factors, which were applied in material balance equations of each model, for calculation of the amounts of all specific sulphur and aromatic compounds through the reactor (Figure 3 and Figure 4).
All kinetic models except the Owusu-Boakye (Model 3) and de Oliveira (Model 5) models predicted poorly the conversion of aromatics in the simulations of the industrial-scale reactor. The enormous deviation of outlet temperature in Model 4 occurred due to the very high value of the suggested HDA reaction heat, which was applied in the energy balance equation of the model. In addition, a large deviation from the temperature experimental data was observed for the Avraam and Liu models, which again proposed very high values for heat of the reaction.
Keeping in mind the experimental data, Model 3 could describe this catalytic reactor with the highest accuracy, or lowest deviation from the experimental data, among the investigated models. Model 3 described the system quite well, with very good predictions of sulphur and aromatics conversion, as well as the reactor temperature increase. Model 5 predicted well the conversion of aromatics, but it underestimated the outlet temperature and overestimated the sulphur conversion.
Novel research focused on blending different pyrolysis oils with light cycle oil [36,37] could be interesting to apply to this model for evaluating such technology in future research.

4. Conclusions

In this study, six different kinetic models for hydrodearomatisation reactions were compared through simulations of the industrial catalytic reactor for hydrotreating of straight run gas oil-containing FCC naphtha and light cycle oil fractions. The influence of hydrodearomatisation reaction kinetics on the conversion of sulphur, temperature increase, and wetting efficiency in the reactor were investigated. Experimental data from the industrial test run was used for validation of different models. Based on the obtained simulation results, it could be concluded that reactor performance is strongly dependent on the hydrodearomatisation reaction with respect to aromatics conversion but even more so with respect to the temperature increase in the reactor, which affects all key catalytic reaction parameters and thus also the sulphur conversion. Calculated total conversions of aromatic compounds differed considerably between models used, depending on the type of hydrodearomatisation kinetics total aromatic conversion ranges from to 13% to 63%. Different hydrodearomatisation kinetics do not affect conversion of the sulphur significantly when expressed as the percentage deviation; maximal absolute deviation from the experimental point was 1.2% conversion, but this relatively low deviation is critical for the quality of prediction since real conversions are expected within the narrow range of around 99.5 +/− 0.5%. Depending on the used kinetic model, the outlet temperature ranged from 607.7 K to 625 K, and the best prediction of experimental temperature was achieved using the Owusu-Boakye model, while the Avraam and Yui models overestimated outlet temperature due to very high values of the reaction heat proposed by these models. Catalyst wetting efficiency is an important parameter in the reactor and affects the reaction rates and overall efficiency of the process, but it cannot be evidenced directly in the industrial reactor. Good predictions of total sulphur and aromatics conversions, along with accurate prediction of temperature increases in the reactor, confirmed the relevance of overall efficiency and phase equilibrium calculations within the model. The best predictor of outlet temperature, sulphur, and aromatics conversion was achieved with the Langmuir–Hinshelwood kinetic model proposed by Owusu-Boakye.

Author Contributions

Conceptualization, methodology, validation, investigation, resources, writing—original draft preparation, writing—review and editing, visualization, supervision, project administration, S.B.G. and A.M.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

This study was supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia, project III45019 (Contract No. 451-03-9/2021-14/200135).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Ttemperature, K
Ppressure, bar
Kivaur–liquid equilibrium constant for component i
FL/FVliquid-to-vapour molar flow ratio
Cpheat capacity, kJ/mol/K
fcatalyst wetting efficiency
ReReynolds number
ΔPpressure drop, Pa/m
ggravitational acceleration, m/s2
GaGalileo number
dpequivalent particle diameter, m
Wcatalyst weight, kg
Cimolar concentration of component i in liquid phase, kmol/m3
rireaction rate for component i, kmol/kgcat/h
ki, σkinetic parameter for component i, kmol/kgcat/h
Ki,σkinetic adsorption parameter for component i, m3/kmol
DEN σoverall adsorption parameter
Henthalpy, kJ/kmol
D e f f , i . effective diffusivity for component i, m2/s
DAB,imolecular diffusivity for component i, m2/s
DK,iKnudsen diffusivity for component i, m2/s
kiAPPIpseudo first-order kinetic constant for component i
kgigas phase external mass transfer coefficient for component i, m/s
kliliquid phase external mass transfer coefficient for component i, m/s
XLockhart–Martinelli number
WeWeber number
ScSchmidt number
Subscripts:
icomponent i
BTbenzothiophene
DBT1methyldibenzothiophene
DBT2dimethyldibenzothiophene
DBT3trimethyldibenzothiophene
H2hydrogen
H2Shydrogen sulphide
Superscripts:
Lliquid
Vvapour
Greek symbols:
ρdensity, kmol/m3
Ωioverall effectiveness factor
ηiinternal effectiveness factor
ΦiThiele modulus
εporosity of the catalyst particle
εbporosity of the catalyst bed
ρpdensity of the catalyst pellet, kg/m3

References

  1. Hayashida, K.; Haji, K. Effects of Fuel Properties on Diffusion Combustion and Deposit Accumulation, Fossil Fuel and the Environment; Khan, S., Ed.; IntechOpen: London, UK, 2012; pp. 1–18. ISBN 978-953-51-0277-9. Available online: https://www.intechopen.com/books/fossil-fuel-and-the-environment/effects-of-fuel-properties-on-diffusion-combustion-and-deposit-accumulation (accessed on 20 June 2021).
  2. Ghosh, P.; Jaffe, S.B. Detailed composition-based model for predicting the cetane number of diesel fuels. Ind. Eng. Chem. Res. 2006, 45, 346–351. [Google Scholar] [CrossRef]
  3. Kidoguchi, Y.; Yang, C.; Miwa, K. Effects of Fuel Properties on Combustion and Emission Characteristics of a Direct-Injection Diesel Engine. J. Fuels Lubr. 2000, 109, 1149–1157. Available online: www.jstor.org/stable/44745922 (accessed on 8 June 2021).
  4. Danaher, W.J.; Palmer, L.D. Chemical changes and ignition quality improvement resulting from hydrotreating light cycle oil. Fuel 1988, 67, 1441–1445. [Google Scholar] [CrossRef]
  5. Quality of Petrol and Diesel Fuel Used for Road Transport in the European Union, Report from the Commission to The European Parliament and the Council; 23.3.2012 COM(2012) 127 Final; European Commission: Belgium, Brussels, 2012; Available online: https://data.consilium.europa.eu/doc/document/ST-8168-2012-INIT/en/pdf (accessed on 20 June 2021).
  6. Mijatovic, I.M.; Glisic, S.B.; Orlovic, A.M. Modelling a catalytic reactor for hydrotreating of straight-run gas oil blended with FCC naphtha and light cycle oil: The influence of vapour-liquid equilibrium. Ind. Eng. Chem. Res. 2014, 53, 19104–19116. [Google Scholar] [CrossRef]
  7. Dukanović, Z.; Glišić, S.B.; Čobanin, V.J.; Nićiforović, M.; Georgiou, C.A.; Orlović, A.M. Hydrotreating of straight-run gas oil blended with FCC naphtha and light cycle oil. Fuel Process. Technol. 2013, 106, 160–165. [Google Scholar] [CrossRef]
  8. Anchyeta, J. Modeling and Simulation of Catalytic Reactors for Petroleum Refining, 1st ed.; John Wiley & Sons: New York, NY, USA, 2011; ISBN 978-0-470-18530-8. [Google Scholar]
  9. Chen, J.; Mulgundmath, V.; Wang, N. Accounting for Vapor-Liquid Equilibrium in the Modeling and Simulation of a Commercial Hydrotreating Reactor. Ind. Eng. Chem. Res. 2011, 50, 1571–1579. [Google Scholar] [CrossRef]
  10. Girgis, M.J.; Gates, B.C. Reactivities, reaction networks, and kinetics in high pressure catalytic hydroprocessing. Ind. Eng. Chem. Res. 1991, 30, 2021–2058. [Google Scholar] [CrossRef]
  11. Stanislaus, A.; Cooper, B.H. Aromatic hydrogenation catalysis: A review. Catal Rev. 1994, 36, 75–123. [Google Scholar] [CrossRef]
  12. Korre, S.C.; Klein, M.T.; Quann, R.J. Polynuclear aromatic hydrocarbons hydrogenation. 1. Experimental reaction pathways and kinetics. Ind. Eng. Chem. Res. 1995, 34, 101–117. [Google Scholar] [CrossRef]
  13. Korre, S.C.; Neurock, M.; Klein, M.T.; Quann, R.J. Hydrogenation of polynuclear aromatic hydrocarbons. 2. Quantitative structure/reactivity correlations. Chem. Eng. Sci. 1994, 49, 4191–4210. [Google Scholar] [CrossRef]
  14. Chowdhury, R.; Pedernera, E.; Reimert, R. Trickle-Bed Reactor Model for Desulfurization and Dearomatization of Diesel. AIChE J. 2002, 48, 126–135. [Google Scholar] [CrossRef]
  15. Wilson, M.F.; Kriz, J.F. Upgrading of Middle Distillate Fractions of a Syncrude from Athabasca Oil Sands. Fuel 1984, 63, 190–196. [Google Scholar] [CrossRef]
  16. Bhaskar, M.; Valavarasu, G.; Sairam, B.; Balaraman, K.S.; Balu, K. Three-Phase Reactor Model to Simulate the Performance of Pilot-Plant and Industrial Trickle-Bed Reactors Sustaining Hydrotreating Reactions. Ind. Eng. Chem. Res. 2004, 43, 6654–6669. [Google Scholar] [CrossRef]
  17. Gunjal, P.R.; Ranade, V.V. Modeling of laboratory and commercial scale hydro-processing reactors using CFD. Chem. Eng. Sci. 2007, 62, 5512–5526. [Google Scholar] [CrossRef]
  18. Mederos, F.S.; Ancheyta, J.; Elizalde, I. Dynamic modeling and simulation of hydrotreating of gas oil obtained from heavy crude oil. Appl. Catal. A Gen. 2012, 425-426, 13–27. [Google Scholar] [CrossRef]
  19. Avraam, D.G.; Vasalos, I.A. HdPro: A mathematical model of trickle-bed reactors for the catalytic hydroprocessing of oil feedstocks. Catal. Today 2003, 79–80, 275–283. [Google Scholar] [CrossRef]
  20. Jiménez, F.; Kafarov, V.; Nuňez, M. Modeling of industrial reactor for hydrotreating of vacuum gas oils—Simultaneous hydrodesulfurization, hydrodenitrogenation and hydrodearomatization reactions. Chem. Eng. J. 2007, 134, 200–208. [Google Scholar] [CrossRef]
  21. Jiménez, F.; Nuňez, M.; Kafarov, V. Study and modeling of simultaneous hydrodesulfurization, hydrodenitrogenation and hydrodearomatization on vacuum gas oil hydrotreatment. Comput. Aided Chem. Eng. 2005, 20, 619–624. [Google Scholar]
  22. Owusu-Boakye, A.; Dalai, A.K.; Ferdous, D.; Adjaye, J. Experimental and Kinetics Studies of Aromatic Hydrogenation in a Two-Stage Hydrotreating Process using NiMo/Al2O3 and NiW/Al2O3 Catalysts. Can. J. Chem Eng. 2006, 84, 572–580. [Google Scholar] [CrossRef]
  23. Yui, S.M.; Sanford, E.C. Kinetics of aromatics hydrogenation and prediction of cetane number of synthetic distillates, Proceedings—Refining Department. In Proceedings of the 50th Midyear Meeting, Kansas City, MO, USA, 13–16 May 1985; American Petroleum Institute: Washington, DC, USA, 1985; pp. 290–297. [Google Scholar]
  24. De Oliveira, L.P.; Verstraete, J.J.; Kolb, M. A Monte Carlo modeling methodology for the simulation of hydrotreating processes. Chem Eng. J. 2012, 207–208, 94–102. [Google Scholar] [CrossRef]
  25. Liu, Z.; Zheng, Y.; Wang, W.; Zhang, Q.; Jia, L. Simulation of hydrotreating of light cycle oil with a system dynamics model. Appl. Catal. A Gen. 2008, 339, 209–220. [Google Scholar] [CrossRef]
  26. Magnabosco, L.M. A mathematical model for cataytic hydrogenation of aromatics in petroleum refining feedstocks. Stud. Surf. Sci. Catal. 1989, 53, 481–496, ISBN 9780080887005. [Google Scholar]
  27. Cheng, Z.; Fang, X.; Zeng, R.; Han, B.; Huang, L.; Yuan, W. Deep removal of sulfur and aromatics from diesel through two-stage concurrently and counter currently operated fixed-bed reactors. Chem. Eng. Sci. 2004, 59, 5465–5472. [Google Scholar] [CrossRef]
  28. Papadias, D.; Edsberg, L.; Bjornbom, P. Simplifed method for effectiveness factor calculations in irregular geometries of washcoats. Chem. Eng. Sci. 2000, 55, 1447–1459. [Google Scholar] [CrossRef]
  29. Sie, S.T.; Krishna, R. Process development and scale up: III. Scale-up and scale-down of trickle bed processes. Rev. Chem. Eng. 1998, 14, 203–252. [Google Scholar] [CrossRef]
  30. Cheng, Z.-M.; Anter, A.M.; Yuan, W.-K. Intensification of phase transition on multiphase reactions. AIChE J. 2001, 47, 1185–1192. [Google Scholar] [CrossRef]
  31. Cheng, Z.-M.; Zhou, Z.-M.; Yuan, W.-K. Determination of catalyst wetting fraction on the molecular level. AIChE J. 2007, 53, 741–745. [Google Scholar] [CrossRef]
  32. Cheng, Z.M.; Anter, A.M.; Fang, X.C.; Xiao, Q.; Yuan, W.K.; Bhatia, S.K. Dryout phenomena in a three-phase fixed-bed reactor. AIChE J. 2003, 49, 225–231. [Google Scholar] [CrossRef]
  33. Naranov, E.R.; Maximova, A.L. Selective conversion of aromatics into cis-isomers of naphthenes using Ru catalysts based on the supports of different nature. Catal. Today 2019, 329, 94–101. [Google Scholar] [CrossRef]
  34. Ramírez-Castelán, C.E.; Hidalgo-Vivas, A.; Brix, J.; Jensen, A.D.; Huusom, J.K. Mathematical modelling and simulation of a trickle-bed reactor for hydrotreating of petroleum feedstock. Int. J. Chem. React. Eng. 2019, 17, 20180176. [Google Scholar] [CrossRef] [Green Version]
  35. Neto, A.T.P.; Fernandes, T.C.R.L.; Silva, H.B., Jr.; de Araújo, A.C.B.; Alves, J.J.N. Three-phase trickle-bed reactor model for industrial hydrotreating processes: CFD and experimental verification. Fuel Process. Technol. 2020, 208, 106496. [Google Scholar] [CrossRef]
  36. Palos, R.; Kekäläinen, T.; Duodu, F.; Gutiérrez, A.; Arandes, J.M.; Jänis, J.; Castañoa, P. Screening hydrotreating catalysts for the valorization of a light cycle oil/scrap tires oil blend based on a detailed product analysis. Appl. Catal. B Environ. 2019, 256, 117863. [Google Scholar] [CrossRef]
  37. Palos, R.; Kekäläinen, T.; Duodu, F.; Gutiérrez, A.; Arandes, J.M.; Jänis, J.; Castaño, P. Detailed nature of tire pyrolysis oil blended with light cycle oil and its hydroprocessed products using a NiW/HY catalyst. Waste Manag. 2021, 128, 36–44. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Overall conversion of aromatics through the reactor calculated by different model equations proposed by: Chowdhury [14], Avraam [19], Owusu-Boakye [22], Yui [23], de Oliveira [24], Liu [25], and experimental data at the reactor outlet for total aromatics conversion.
Figure 1. Overall conversion of aromatics through the reactor calculated by different model equations proposed by: Chowdhury [14], Avraam [19], Owusu-Boakye [22], Yui [23], de Oliveira [24], Liu [25], and experimental data at the reactor outlet for total aromatics conversion.
Energies 14 04616 g001
Figure 2. Overall sulphur conversion simulation (a) and temperature increase along the reactor bed (b) with different kinetic models for aromatics conversion.
Figure 2. Overall sulphur conversion simulation (a) and temperature increase along the reactor bed (b) with different kinetic models for aromatics conversion.
Energies 14 04616 g002
Figure 3. Catalyst wetting efficiency through the reactor.
Figure 3. Catalyst wetting efficiency through the reactor.
Energies 14 04616 g003
Figure 4. Overall effectiveness factor for: (a) liquid phase, (b) vapour phase calculated using six different kinetics models for HDA reactions.
Figure 4. Overall effectiveness factor for: (a) liquid phase, (b) vapour phase calculated using six different kinetics models for HDA reactions.
Energies 14 04616 g004
Table 1. Properties of input and output streams used in the industrial scale test run.
Table 1. Properties of input and output streams used in the industrial scale test run.
ParameterInlet StreamOutlet Stream
Straight Run Gas OilFluid Catalytic Cracking Naphtha and Light Cycle OilHydrotreated Gas Oil
Distillation range, °C 33–335
Sulphur content, % wt.0.73051.1400
Thiol/mercaptan sulphur, % wt.0.01620.02710
Paraffin and naphthenic content, % vol.68.732.173.0
Olefins content, % vol.5.216.90
Aromatics content, % vol.26.151.027.0
Monoaromatic content, % wt.19.410.927.4
Diaromatic content, % wt.11.55.14.6
Triaromatic content, % wt.1.4116.00.70
Table 2. Equations used in the model.
Table 2. Equations used in the model.
Parameter:Equation:
Equilibrium constants K i = A + B · T + C · T 2 + D · T 3 ; i = pseudocomponent
Liquid-to-vapour flow ratios F L F V = a + b · T + c · T 2 + d · T 3
Densities of the vapour and liquid phases ρ i = A i + B i · T ; i = V , L
Heat capacities C p = α + β · T + γ · T 2 + δ · T 3
Catalyst wetting efficiency f = 1.104 · Re L 1 / 3 · [ 1 + [ ( Δ p / Δ z ) / ρ L g ] G a L ] 1 / 9
Overall catalyst effectiveness; the internal effectiveness factor;
generalized Thiele modulus
Ω f = X f i Φ i f 2 η i f X f i Φ f i 2 ; η i f = tanh ( Φ i f ) Φ i f ; Φ i f = V p S p k i , A P P I , f C i f ρ D e f f , i f
where i = BT, DBT1, DBT2, DBT3, and A1, f = phase (liquid or vapour)
Material balance equations for reacting components in the process of hydrodesulphurisation: BT, DBT1, DBT2, and DBT3 d F i V d W = ( 1 f ) · Ω i V · r i , H D S ; d F i L d W = f · Ω i L · r i , H D S
Material balance equations for reacting components in the process of hydrodearomatisation d F A 1 V d W = ( 1 f ) · Ω A 1 V · r A 1 , A R M ρ p · ( 1 ε p ) ; d F A 1 L d W = f · Ω A 1 L · r A 1 , A R M ρ p · ( 1 ε p )
The overall energy balance Δ E H D S = 1 i F i V · ( 1 f ) · Ω i V · r i , H D S V · Δ H H D S + 1 i F i L · f · Ω i L · r i , H D S L · Δ H H D S Δ E H D A = 1 i F j V · ( 1 f ) · Ω j V · r j , H D A V · Δ H H D A + 1 i F j L · f · Ω j L · r j , H D A L · Δ H H D A d T d W = Δ E H D S + Δ E H D A i F i V · C p i , V + i F i L · C p i , L
Table 3. Selected kinetic models for HDA reactions used in this study.
Table 3. Selected kinetic models for HDA reactions used in this study.
ModelRef.CatalystPressure, Temperature, and Gas/Diesel Fraction RatioKinetic ExpressionsStoichiometric Equations
Model 1
Chowdhury
[14,15]NiMo/Al2O32–8 MPa 300–380 °C
100–500 m3(NTP)/m3
r p o l y = k p o l y · c p o l y · P H 2 n 3 + k p o l y · c d i
r d i = k d i · c d i · P H 2 n 2 + k d i · c m o n o
r m o n o = k m o n o · c m o n o · P H 2 n 1 + k m o n o · c N a p h
P o l y a r o m a t i c + H 2 k P o l y k P o l y D i a r o m a t i c D i a r o m a t i c + 2 H 2 k D i k D i y M o n o a r o m a t i c M o n o a r o m a t i c + 3 H 2 k M o n o y k M o n o N a p h t h e n e
[16]NiMo/Al2O34 MPa 320–360 °C
200 L/L
[17]-20–80 MPa
320–380 °C
100–500 m3/m3
[18]CoMo/Al2O35.3 MPa
340–380 °C
356 std m3/m3
[26]W-Mo-Ni-Co/Al2O3
Ni-Mo/Al2O3
6.5 MPa
350–360 °C
Model 2
Avraam
[19,20]CoMo/&-Al2O33–7 MPa
350–390 °C
0.1–3 Nm3/m
r H D A = K H 2 · C A 1 1 + K 1 , H 2 S · C H 2 S + K D · C A 2 1 + K 1 , H 2 S · C H 2 S A 1 + α H 2 A 2
[21]NiMo/Al2O36–10 MPa
330–390 °C
Ratio = 4.5–6.25
Model 3
Owusu-Boakye
[22]NiMo/Al2O3
NiW/Al2O3
11.0 MPa
350–390 °C
550 mL/ml
r A i = d C A i d t = k A i · K A i · K H i · P H 2 · C A i 1 + K A i · C A i + K H 2 S i · P H 2 S i A r o m a t i c s + n H 2 S a t u r a t e s
Model 4
Yui
[23]NiMo/Al2O35.8 Mpa
380 °C
2000 (scf/Bbl)
r H D A = k f · P H 2 · c A k r · ( 1 c A ) A k f k r B
Model 5
de Oliveira
[24]CoMo/Al2O3NiMo/Al2O37.0 MPa (H2)
320 °C
r H y d r o = C H y d r o · P H 2 · x i = k s r · K a d s · P H 2 · x i
r D e h y d r o = C D e y d r o · x i = C H y d r o K e q · x i = k s r · K a d s · K e q x i
Hydrogenation of an aromatic ring
D i a r o m a t i c + 2 H 2 M o n o a r o m a t i c
Dehydrogenation of a saturated ring
M o n o a r o m a t i c D i a r o m a t i c + 2 H 2
Model 6
Liu
[25]-7.1 MPa
350 °C
1000 NL/L
r S a t u r a t e s = k a 1 · C M o n o n a 1 · η M o n o
r M o n o = k a 1 · C M o n o n a 1 · η M o n o k a 2 · C D i n a 2 · η D i
r D i = k a 2 · C D i n 2 · η D i k a 3 · C T r i n a 3 · η T r i
r T r i = k a 3 · C T r i n a 3 · η T r i k a 4 · C T e t r a n a 4 · η T e t r a
r T e t r a = k a 4 · C T e t r a n a 4 · η T e t r a
T e t r a k a 4 T r a k a 3 D i k a 2 M o n o k a 1 S a t u r a t e s
Table 4. Calculated kinetic parameters and rate constants for proposed models.
Table 4. Calculated kinetic parameters and rate constants for proposed models.
ModelKinetic ConstantEa, kJ/molValue
Model 1
Chowdhury
k*poly, s−1
k*di, s−1
k*mono, s−1
645.9 × 10−5
2.8 × 10−5
2.8 × 10−5
Model 2
Avraam
kH, s−1
kD, s−1
1801.2 × 10−1
2.1 × 10−2
Model 3
Owusu-Boakye
kA, s−1
KA
KH2
KH2S
858.7 × 10−3
1.3 × 10−5
2.5 × 10−2
1.3 × 10−5
Model 4
Yui
kf
kr
2553.9 × 10−6
7.6 × 10−7
Model 5
de Oliveira
ksr, kmol kg−1 s−1
Kads, m3 kmol−3
keq, m3 kmol−3
512.2 × 10−5
1.5 × 10−5
5.9 × 10−8
Model 6
Liu
k1, s−1
k2, s−1
k3, s−1
k4, s−1
1251.5 × 10−4
1.7 × 10−4
1.5 × 10−4
1.9 × 10−4
Table 5. Comparison of model results and experimental data.
Table 5. Comparison of model results and experimental data.
ModelΔT in Reactor, KDeviation from the Experimental Value, %Overall Sulphur Conversion,
mass%
Deviation from the Experimental Value, %Aromatics Conversion, mass%Deviation from the Experimental Value, %
Model 1
Chowdhury
10.112.0899.610.1238.6031.40
Model 2
Avraam
16.846.0898.60.8948.6253.84
Model 3
Owusu-Boakye
11.40.8799.410.0171.802.57
Model 4
Yui
20.981.7499.800.3146.3033.86
Model 5
de Oliveira
8.823.599.740.2568.342.37
Model 6
Liu
15.938.2697.701.8052.0725.62
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Glišić, S.B.; Orlović, A.M. The Influence of Hydrodearomatisation Reaction Kinetics on the Modelling of Sulphur and Aromatics Removal from Diesel Fuel in an Industrial Hydrotreating Process. Energies 2021, 14, 4616. https://doi.org/10.3390/en14154616

AMA Style

Glišić SB, Orlović AM. The Influence of Hydrodearomatisation Reaction Kinetics on the Modelling of Sulphur and Aromatics Removal from Diesel Fuel in an Industrial Hydrotreating Process. Energies. 2021; 14(15):4616. https://doi.org/10.3390/en14154616

Chicago/Turabian Style

Glišić, Sandra B., and Aleksandar M. Orlović. 2021. "The Influence of Hydrodearomatisation Reaction Kinetics on the Modelling of Sulphur and Aromatics Removal from Diesel Fuel in an Industrial Hydrotreating Process" Energies 14, no. 15: 4616. https://doi.org/10.3390/en14154616

APA Style

Glišić, S. B., & Orlović, A. M. (2021). The Influence of Hydrodearomatisation Reaction Kinetics on the Modelling of Sulphur and Aromatics Removal from Diesel Fuel in an Industrial Hydrotreating Process. Energies, 14(15), 4616. https://doi.org/10.3390/en14154616

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop