Abstract
A detailed steady-state catalytic-reforming unit (CRU) reactor process model is simulated in this work, and for the first time, different compressibility Z factor correlations have been applied using gPROMS software. The CRU has been modeled and simulated with the assumption that the gas phase behaves like an ideal gas. This is assumed for the four reactors in series and for different conditions of hydrogen–hydrocarbon ratio (HHR), operating temperature, and pressure. The results show that the Z factor varies at every point along the height of the reactors depending on reaction operating pressure, temperature, and HHR ratio. It also shows that the magnitude of deviation from ideal gas behaviour can be measured over the reactor height. The Z factor correlation of Mahmoud (J Energy Resour Technol Trans ASME 136:012903, 2014) is found to be suitable for predicting the Z factor distribution in the reactors.
Similar content being viewed by others
Introduction
The catalytic-reforming unit is an important integral part of the refinery operations. It is used to improve the quality of low-to-high-octane naphtha in a series of reactors (three or four, depending on the design of the refinery). The system can either be a semi regenerative (fixed bed) where the reaction proceeds for a particular period and the catalyst is regenerated, or a continuous one where the catalyst is continuously regenerated in a regenerator. The continuous process has the advantage of using a lower pressure and lower hydrogen usage. However, there is the need for more technical monitoring for the mechanisms of engagement and disengagement of the catalyst from the last reactor through the regenerator and back to the first reactor in continuous flow. To improve the quality and the yield of the reformate, one important aspect is the hydrodynamics study of the reaction which helps to improve the efficiency of the system by maintaining an effective pressure gradient across the reactors. The low octane naphtha meets recycled hydrogen from the recycle gas compressor and enters the first heater that precedes the reactor, a process that continues for the other four reactors and heaters in series. The hydrodynamics is important, since the variations of temperature and pressure of the reaction affect the rate of reactions, which in turn affect the entire system responses. Increase in reaction pressure decreases the reformate yield, the hydrogen yield, and the research octane number, RON [10]. The Ergun equation, as a function of density, viscosity, bed void, particle diameter, and mass flux, is used to model the behaviour of pressure in these reactors. Since these reactions’ parameters such as enthalpy, density, viscosity, and pressure change with reactor height, there is need to study the effect of the compressibility factor which affects these parameters [13]. In the past, all the CRU was modeled with compressibility factor Z as unity, implying an ideal gas system, which may not always be the case. Real systems tend to vary from ideal system. The Z factor of a fluid catalytic cracking (FCC) riser has been modeled as unity until John et al. [13] simulated the riser with different Z factor correlations and found that the riser can be modeled with Heidaryan et al. [11] correlation. Hence, the need to find the actual Z factor for the CRU and consequently investigate the response of the CRU under real system conditions.
In characterizing fluid flow behaviour in the oil and gas, it is very important to study the effect of Z factor, both upstream and downstream [12]. The fluid can be said to be compressible or incompressible depending on the type of process which it undergoes [11]. Variation in density of the gas during reaction, as in the case with catalytic-reforming, may bring about change in compressibility factor. Therefore, assuming a changing gas density system as ideal may not be always true especially when there is velocity variation as the gas density varies and the fluid could be compressible [5]. Transport and physical properties like density, viscosity, and fraction of void of the treated naphtha during reaction could change when reaction conditions [temperature, pressure, feed rate, and hydrogen–hydrocarbon ratio (HHR)] are altered. Since there is significant variation of these properties with operating conditions, there is a need to study the compressibility factor effect. An important variable in the process conditions is the HHR which preserves the catalyst activity by sweeping off amorphous carbon deposit on the catalyst, but has little effect on the aromatics and reformate yield. This variable has an effect on the reactor pressure, since it increases the partial pressure of the vapor by increasing the number of moles of hydrogen; hence, it will have effect on the reaction hydrodynamics. To have optimum and precise condition of the catalyst, the reactor bed and other equipment should adhere to professional design of the catalytic-reforming plant. In this paper, the effect of compressibility factor on reaction pressure, which is a major hydrodynamic parameter, is studied. This will help in determining an appropriate gas compressibility factor to be applied in plant design and when any modifications of the design are required such as scale up or scale down, the assumption of ideal case may not always hold. The precise compressibility factor can assist to predict accurately the hydrodynamic behaviour of operational variables like pressure drop across the reactors to ensure effective process plant design.
In this work, the effect of compressibility factor on four commercial catalytic naphtha reactors in series is investigated for the first time, by applying various correlations to model the behaviour of the reaction using gPROMS software. Hence, the Z factor variation across the reactors heights will be determined for different Z factors and a suitable correlation model for the CRU will be determined.
Compressibility factor Z
The compressibility factor of gases (Z factor) from principle of corresponding of state is defined based on the pseudo-reduced temperature (Tpr) and pseudo-reduced pressure (Ppr), which are important thermodynamic variables when determining the behaviour of gases and liquids both in upstream and downstream computations [11]. These are described in Eqs. (2) and (3):
The three equations are applied to both real and ideal gases where the Z factor is unity for ideal cases which in reality is non-existent. It is thus of paramount importance to predict the compressibility effect in gaseous phase reactions when dealing with changes of physical and transport properties. The simple definition is the ratio of actual gas volume to the gas volume of ideal gas implying a measure and extent of deviation from ideal behaviour [12].
According to Fayazi et al. [9], this parameter can be determined experimentally or from equations of state or using semi-empirical correlations. The use of experimental data seems more expensive and takes a lot of time and energy and that could be so tedious considering the number of gases in petroleum to account for Ahmed [1], whilst the use of semi-empirical correlations has proved accurate and simpler than even the use of equations of state EoS [8]. Whenever the pseudo-reduced pressure and pseudo-reduced temperature of the gas is known, the Z factor of the vapor of the hydrocarbon could be predicted and estimated [9].
In this work, Tpr and Ppr are computed using Eqs. (2) and (3). The Tpr and Ppr are shown in Fig. 1 with variation along the reactor heights. The Ppr obtained is in the range \(1.218066 \le P_{\text{pr}} \le 1.023427\) and its Tpr is within the range \(0.528144 \le T_{\text{pr}} \le 0.348992\). The values of the Tpr and Ppr may change based on the conditions of operations of the reaction. This implies that, as the variables of the process which affect the temperature and pressure of the CRU vary during plant run, the Tpr and Ppr will also vary. Conversely, the compressibility, a function of Tpr and Ppr will as well not remain constant but change.
There are some common empirical correlations [6, 14] that are not suitable when Tpr ≤ 0.92. Some of the correlations applied in this research accept Tpr above 0.92 [11, 17]. In the quest to determine the most accurate and precise Z factor for the vapor state of the reactions, a number of different empirical correlations are used. Each of the Z factors determined is compared with both the plant data and that of the literature and ascertain which of the compressibility factors predicts closely. The computed pseudo-reduced temperature here is outside the boundary of some of the empirical correlations used, but the Ppr is within the range \(0.364 \le P_{\text{pr}} \le 0.375\) which lies within the boundaries of the Ppr in the used empirical correlations in this work which are given below.
- 1.
Azizi et al. [3] Z factor:
Azizi et al. [3] established their Z factor empirical correlation by applying standing Katz chart with about 3038 points with a range of Ppr of \(0.2 \le P_{\text{pr}} \le 11\) and Tpr range of \(1.1 \le T_{\text{pr}} \le 2\). This is presented in Eq. (4):
The variables in Eq. (4) are presented in Eqs. (5–9):
The tuned coefficients for Eqs. (5–9) are presented in Appendix Table 14.
- 2.
Bahadori et al. [4] compressibility factor:
Compressibility factor of Bahadori et al. [4] is given in Eq. (10) and its coefficients are presented in Eqs. (10–14) [4]. The range is \(0.2 \le P_{\text{pr}} \le 16\) and \(1.05 \le T_{\text{pr}} \le 2.4\):
The tuned coefficients for Eqs. (10–14) are presented in Appendix Table 15.
- 3.
Compressibilty factor Heidaryan et al. [11]:
The compressibility factor for Heidaryan et al. [11] is defined as in Eq. (15) and the fine-tuned coefficients are given in Appendix Table 16. The Ppr is within \(0.2 \le P_{\text{pr}} \le 3\) and this is within the range of that of Heidaryan et al. [12]:
- 4.
Heidaryan et al. [12] compressibility factor:
The compressibility factor of Heidaryan et al. [11, 12] is defined in Eq. (16) with the tuned coefficients reported in Appendix Table 17. The boundaries of the Ppr and Tpr are \(0.20 \le P_{\text{pr}} \le 15.0\) and \(1.20 \le T_{\text{pr}} \le 3.0\) (Heidaryan et al. 2010c). The boundary of the Ppr here in this research is concurrent as that of Heidaryan et al. [11]:
- 5.
Mahmoud [16] compressibility Z factor:
The compressibility Z factor for Mahmoud [16] is defined by Eq. (17). The correlation was derived from measurements taken of 300 Z factors [16]:
- 6.
Papay [17] compressibility factor:
The compressibility factor Z of Papay [17] is defined by Eq. (18) [15]:
- 7.
Z factor correlation of Sanjari and Lay [20]:
Sanjari and Lay [17] got their compressibility factor correlation from 5844 experimental data points of different compressibility factors within the boundary of 0.010 ≤ Ppr ≤ 15.0 and 1.0 ≤ Tpr ≤ 3.0. It is reported in Eq. (19) with the fine-tuned coefficients given from Appendix Table 18:
- 8.
Shokir et al. [21] Z factor:
The Shokir et al. [21] compressibility factor is defined by Eq. (20), with its parameters given in Eqs. (21–25) [21]:
Kinetic models
The assumptions made, model equations used, and the method and procedures are presented in this section. The reactors are in series, as shown in Fig. 2. The first reactor is 5.632 m in height, 5.83 m for the second, 6.51 m for the third, and 7.26 m for the fourth reactors. The total height of the reactors in series is 25.232 m.
Some of the assumptions made are as follows:
The reactors are modeled as adiabatic processes, i.e., no heat escapes or enters the reactors due to sufficient lagging.
The reactors are modeled as plug flow, because naphtha flows as gas through the reactor at the pressure and temperature of reaction. Due to the reactor length, dispersion of matter is negligible as supported by a criterion reported elsewhere [10].
The reaction is considered as first order, because the reforming reaction uses HHR higher than the hydrocarbon concentration. Hence, its concentration is grouped with reaction coefficient.
Pseudo-homogeneous reactor model is assumed, since the vapor phase reaction is too complex to solve heterogeneously and the diffusion along the radius is negligible with the height of the reactor much greater than the radius.
Various models have been developed for both steady-state and dynamic processes of the catalytic-reforming reaction with different lumps to investigate the behaviour of the reaction and product distribution. The reaction steps and equations of the model are shown in Table 1 (a zero value is registered when the data are not available).
The kinetic rate equations for the components are given in Eqs. (26–50) where SV is space velocity and \(\frac{1}{\text{SV}}\) is residence time:
where Ki is the kinetic constants for the reactions, T0 and P0 are the reference temperatures and pressures, EAi the activation energy, and w is an exponential effect of pressure, as given in Tables 2 and 3, respectively. ki, EAi, w, T0, and P0 are values obtained from Elizalde and Ancheyta [7].
Mathematical model
Modeling the behaviour of the reactions is done by solving the model equations describing the system. The equations describing the mass balance and heat balance are solved simultaneously on gPROMS as represented in Eqs. (52–57):
Equations (26–50) are the rate equations for all the components, while Eq. (51) is the kinetic rate constant applied to all the rate equations.
Equation (52) is the material balance equation where \(\rho b\) is the bulk density and \(\varepsilon\) the bed voidage. Equation (53) is the heat balance of the reaction to determine the temperature behaviour where \(- \Delta H_{Rj}\) is the heat of the reaction of the jth component and Cpcat is the heat capacity of the catalyst. Equation (54) gives the pressure profile of the reaction where Dp is particle diameter, G the mass flux, \(\varepsilon\) the bed voidage, \(\rho\) the gas density, and A the reactor area. \(H_{ri}^{0}\), the standard heat of formation of the components and the constants of A1, B1, C1, and D1 for the Cp, heat capacity, in Eqs. (56, 57), are taken from Riazi [18]. These equations are solved simultaneously using gPROMS to determine the behaviour of the system. gPROMS is a robust mathematical software that could solve all the four reactors in series dynamically given the behaviour of the paraffins, naphthenes and aromatics as well as the dynamics of the temperature in the reactors. All the components’ behaviour was determined.
RON model
The estimation of the RON of hydrocarbons for the feed and products can be done using different methods of prediction. RON can be calculated for each pure component using a polynomial equation that correlates to the normal boiling point [18] as shown by Eq. (58):
where \(T = {\text{TBP}} \times 0.01\), TBP represents the normal boiling point (°C), and Aa, Bb, Cc, Dd, and Ee, are coefficients. The RON of a hydrocarbon mixture is calculated by assuming that the mixture consists of paraffins, naphthenic hydrocarbons, and aromatics. The equation is expressed as the sum of the RON for each pure component multiplied by the volume fractions of the components.
Model validation
The model validation was carried out by modeling and simulating the commercial data of Ancheyta et al. [2] and Rodríguez and Ancheyta [19] on gPROMS to ascertain the capability and ruggedness of the mathematical software in complex modeling. The modeling of Ancheyta et al. [2] was performed using MATLAB software and the result is compared with that obtained with gPROMS. The configuration of the commercial reformer and feed properties are given in Tables 4 and 5, respectively.
Model validation results
The commercial reformer Ancheyta et al. [2] were simulated with the properties of Tables 4 and 5, respectively, with four semi regenerative reactors in series. The reformer’s throughput is 30 MBPD at inlet temperature of 495 °C and pressure of 10.5 kg/cm2 using MATLAB mathematical tool on ODE45. The simulation was performed on gPROMS software to validate its capability in modeling this problem. The result obtained from gPROMS is compared with that obtained with MATLAB from Ancheyta et al. [2], and presented in Table 6.
From the modeling and simulation results, the Ancheyta et al.’s [2] result and that of gPROMS were compared and analyzed. The result showed a good and comparable result between the actual and that simulated with gPRMOs. This is an indication of the capability of the mathematical tool in solving complex model equations
Results and discussion
Simulation analysis
The model was validated with commercial data from Ancheyta et al. [2] and Rodríguez and Ancheyta [19] and Kaduna refinery and petrochemical company (KRPC), and shows good agreement with data from both sources. The commercial plant modeling and simulation results are in good agreement, as shown in Fig. 3 and Table 7. Tables 8 and 9 show the configuration and feed properties of KRPC commercial catalytic reformer.
Figure 4 shows the concentration profile of paraffins along the reactor height. P5, P6, P7, P8, P9, and P10 are paraffins with hydrocarbon numbers 5, 6, 7, 8, 9, and 10, respectively. Normal paraffin components basically undergo three major reactions during naphtha reforming. Firstly, they undergo hydrocracking to lighter paraffins, i.e., methane, ethane, propane, and butane, as shown in reactions 8–33 in Table 1. Secondly, the isomerization reaction to form isoparaffins. This is a slow reaction with slow reaction rate, hence, it is not considered in this work. Thirdly, the dehydrocyclization to naphthenes, a very slow reaction, leads to a decrease in the paraffins. These reactions are reactions 1–7 from Table 1. Dehydrocyclization reaction becomes easier as the molecular weight of the paraffins increases as in P8, P9, and P10, while P7 shows little increase. P5 and P6 increase due to hydrocracking. The main effects of hydrocracking are decrease of paraffins (C5+) in the reformate, decrease in hydrogen production, and increase in LPG production and hydrogenolysis. The isomerization reactions are fast, slightly exothermic and do not affect the number of carbon atoms. The thermodynamic equilibrium of isoparaffins to paraffins depends mainly on the temperature and pressure which has no effect. The paraffins isomerization results in a slight increase of the octane number. These reactions are promoted by the acidic function of the catalyst support. The paraffin dehydrocyclization step becomes easier as the molecular weight of the paraffin increases. From Table 1, the rates increase from 0.00, 0.58, 1.33, 1.81, and 2.54 for P6–P11. However, the tendency of paraffins to hydrocrack increases concurrently. Kinetically, the rate of dehydrocyclization increases with low pressure and high temperature. To sum up, the dehydrocyclization of P6 paraffins to benzene is more difficult than that of C7 paraffin to toluene, which itself is more difficult than that of C8 paraffin to xylenes. Accordingly, the most suitable fraction to feed a reforming process is the C7–C10 fraction.
Figure 5 shows the behaviour of the naphthenes along the reactor height. N6, N7, N8, N9, and N10 are naphthenes with hydrocarbon number 6, 7, 8, 9, and 10, respectively. The fastest reaction is aromatization, i.e., dehydrogenation of naphthenes, resulting in the large temperature drop due to its highly endothermic nature. The drastic reduction is due to the aromatization to aromatics. Thermodynamically, the reaction is highly endothermic and is favored by high temperature and low pressure. In addition, the higher the number of carbon atoms, the higher the aromatics production at equilibrium from N8, N9, and N10. This can be seen from Table 1 in Eqs. (34–39) where the reaction rates increase from 4.02, 9.03, 21.5, 24.5, and 24.5, respectively, for N6, N7, N8, N9, and N10. From a kinetic point of view, the rate of reaction increases with temperature. The naphthenes also undergo hydrocracking to lighter hydrocarbons leading to their decrease as shown in reactions (49–57) in Table 1.
Figure 6 shows the behaviour of aromatics along the reactor height. A6, A7, A8, A9, and A10 are paraffins with hydrocarbon numbers 5, 6, 7, 8, 9, and 10, respectively. The dehydrogenation of naphthenes, as shown in Fig. 5, increases the aromatics, as shown in Fig. 6. The sharp increase is because of drastic aromatization of the naphthenes favored by the metallic sites of the catalyst. The rate of benzene formation is lower due to lower carbon number of N6, while A7, A8, and A9 increase rapidly due to higher carbon number of N7, N8, and N9.
Figure 7 shows the temperature profile of the reactions along the reactor height and the sharp drop in first and second reactors is due to the more endothermic dehydrogenation reaction of naphthenes to aromatics. Typically, dehydrogenation and isomerization reactions take place in the first reactor, dehydrogenation, isomerization, dehydrogenation, and cracking in the second followed by dehydrogenation and cracking in the third and fourth reactors.
Figures 8 and 9 show how the increase in hydrogen yield and RON varies, respectively, along the height of the reactor. This is due to the increase in the aromatics along the reactor height. The RON increases from reactor one to four due to increase in aromatization reaction, which gives a higher research octane number, a parameter for antiknock in the gasoline engine. Treated naphtha has low RON and cannot be used as gasoline, hence the reforming reaction of the components to give a gasoline with higher research octane number.
Analysis of compressibility factor
In this study, initially, the compressibility factors of different empirical correlations were added to catalytic-reforming reactor model using HHR of 6.4. The pressure profile of the four reactors in series is given in Fig. 10. The pressure decreases continually from the first reactor and decreases from 18.65 to 12.6 kPa at the exit of the last reactor. In Fig. 11, the profile of the gas densities behaves in similar way due to the drop and decrease in the reactors pressure, but the last reactor exhibits a slightly different behaviour due to the larger pressure drop in the last reactor.
Figure 12 shows the compressibility factor profiles along the reactor heights of the different correlations used. The correlation of Shokir et al. [21] gave a non-zero Z factor along the reactor heights, because the range of Ppr and Tpr is wider than others. The Heidaryan et al.’s [12] compressibility factor was negative and thus would not give meaningful data and so it is not reported. The Z factor changes along the reactor heights due to the variation of operational conditions and transport properties aforementioned. The Z factor at each height of reactor is not the same, and it is obvious that if these variables change with the Z factors as in Fig. 12, the Z factor cannot be a constant value of unity as always assumed. The Mahmoud [16] Z factor is the closest to the ideal compared to others, although this does not signify that it is the true representation of the Z factors, but further analysis will be carried out on how it correlates or deviates from other process variables in the reactors. Factors such as the yield of reformate, hydrogen, and aromatics yield for each compressibility correlation, and the temperature and pressure profiles along the reactors’ heights will need to be considered as well.
ZB represents Bahadori et al. [4], ZH1 represents Heidaryan et al. [11], ZM represents Mahmoud [16], ZP represents Papay [17], ZS represents Sanjari and Lay [20], and ZSH represents Shokir et al. [21].
Figure 13 shows the profiles of the density of the gas phase along the reactor heights. Gas density is a function of pressure and temperature as well as feed flow rate. The pressure model equation is also a function of gas density as in Eq. (54). Therefore, the density is significant in determining the Z compressibility factors and the pressure profiles. The density profile decreases along the reactor height due to the drop in temperature across the various reactors. The behaviour is true for all the correlations and Fig. 12 shows that the Z factors have different profiles of gas densities. However, in Fig. 13, the profile of gas density for Papay [17] is the closest to that of the density for ideal vapor unlike Fig. 12 where the compressibility factor correlation of Mahmoud [16] is the closest. This implies that further analysis needs to be carried out.
The temperature variation with Z factor, as shown in Fig. 14, has effect on the product quality and yield generally as a result of the effect of the kinetic equations which are dependent on reaction temperature. Therefore, enthalpy of various compressibility factors correlations could also vary simultaneously, but the effect is obvious in the last reactor temperature decrease, as shown in Fig. 14. Figure 15 shows reformate yield along the reactors’ heights of various Z factor correlations. As with the case of temperature, where the enthalpy of reactions has a pronounced effect, the pressure, HHR, the reformate yield, and hydrogen yield also change due to the variation of the different Z factors as shown in Figs. 15 and 16. For the reformate and hydrogen yields, Table 10 gives the different exit values of the RMT and hydrogen yields with the values from the Mahmoud [16] being the closest.
For an accurate prediction of a suitable compressibility factor correlation for the CRU modeling, a study of the hydrogen-to-hydrocarbon ratio, which is a major dynamic variable that influences the reaction hydrodynamics, is performed. The variation of reaction pressure is investigated also for the compressibility empirical correlations and a comparison with the pressure values of the model results and that of the ideal case is depicted in Fig. 15
The various pressures vary along the reactor height with some above the ideal gas correlations and some below are shown in Fig. 17. The closest to ideal behaviour is the Mahmoud Z factor correlation [16].
Since the Mahmoud’s [16] compressibility factor correlation is the closest to that of the ideal case, a further analysis is performed between the two by varying the pressure of the reaction and the HHR. This will give the extent of closeness of the product distribution. Tables 7 and 8 show the pressure drops with variation in HHR, while Fig. 16 shows the effects of variation of the pressure using the Mahmoud correlation with Z equal to one. For HHR of 6.4, it has a 0.709% deviation at the exit of the first reactor, 1.734% at the exit of the second reactor, 3.118% at the exit of the third reactor, and 4.343% at the exit of the fourth reactor.
For HHR of 7.0, it has a 0.796% deviation at the exit of the first reactor, 2.00% at the exit of the second reactor, 3.720% at the exit of the third reactor, and 5.319% at the exit of the fourth reactor. The decrease in pressure for the various empirical correlations for different HHR ratios of 6.4 and 7.0 is given in Tables 12 and 13.
The increase in pressure shows a match between the Mahmoud’s [16] correlation and the ideal, as shown in Fig. 18. At the different exit of the reactor, the temperatures are different due to fall in temperature. The Z compressibility factors are different which are functions of the gas densities which are functions of pressure. The Mahmoud [16] correlation becomes more suitable and applicable at a higher pressure and lower HHR.
Conclusion
A detailed steady-state model of catalytic-reforming unit (CRU) of four reactors in series is modeled in this work and a simulation with various compressibility factors is performed using gPROMS, a mathematical and modeling software. The conclusions are:
- 1.
The simulated results from this study are compared with data from the KRPC industrial plant. The parameters of both kinetic and thermodynamic are obtained from the open literature [7]. Good agreement between the simulated and plant data shows the robust strength and capability of the gPROMS model builder in modeling and simulating four reactors in series. Therefore, gPROMS can be recommended for modeling complex processes such as simulation of the CRU unit.
- 2.
Mahmoud’s [16] compressibility Z factor is found to be a more suitable correlation in predicting the Z factor across the four reactors.
- 3.
Different yields of profiles for process variables such as density, reaction pressure, and enthalpy of reaction are obtained with corresponding different compressibility factors in the simulation of the CRU using the same process conditions. These profiles are as a result of varying temperature profiles and varying RMT yield as well as hydrogen yield.
- 4.
The pressure in each reactor is not the same for different HHR ratios. The reaction pressure is also not the same in each of the reactors when the Mahmoud [16] compressibility Z factor correlation is applied rather than assuming and treating the vapor as an ideal gas.
- 5.
The pressure drops across the four reactors are similar and comparable when Mahmoud’s [16] correlation is applied and gives almost the same result at an inlet pressure of 2265 kPa. Hence, Mahmoud’s [16] correlation can be used in modeling the pressures in CRU.
Abbreviations
- \(a,b,c\) :
-
Parameters from hydrogen reaction rate equation
- P :
-
Paraffins
- N :
-
Naphthenes
- \(A\) :
-
Aromatics
- \(A,B,C,D\) :
-
Constants for calculating heat capacities
- \(Aa,Bb,Cc,Dd,Ee\) :
-
Constants for calculating research octane number (RON)
- \(A_{10}\) :
-
Aromatics having ten atoms of carbon
- \(N_{10}\) :
-
Naphthenes having ten atoms of carbon
- \(P_{10}\) :
-
Paraffins having ten atoms of carbon
- \(Cp\) :
-
Heat capacity (kJ/kmol k)
- \(d_{\text{p}}\) :
-
Particle diameter (m)
- \(E_{\text{A}}\) :
-
Activation energy (kJ/kmol)
- \(F\) :
-
Molar flow (kmol/h)
- \(g_{\text{c}}\) :
-
Force to mass conversion factor, 9.8066 (\({\text{kg}}_{\text{m}}\) m)/(\({\text{kg}}_{\text{f}} \,{\text{s}}^{2}\))
- \(G\) :
-
Mass velocity (kg/m2 h)
- \(\Delta G^{^\circ }\) :
-
Reaction standard Gibbs energy (kJ/kmol k)
- \(\Delta H\) :
-
Heat of reaction (kJ/kmol k)
- \(k_{i}\) :
-
Kinetic constant at T (kmol/h)
- \(k_{i}^{^\circ }\) :
-
Kinetic constant at T (kmol/h)
- \(k_{\text{e}}\) :
-
Equilibrium constant
- \({\text{WHSV}}\) :
-
Weight hourly space velocity (h−1)
- \(MW\) :
-
Molecular weight (g/gmol)
- \(n\) :
-
Reaction order
- \(N\) :
-
Number of reactions naphthenes
- \({\text{NC}}\) :
-
Number of components
- \(P_{i}\) :
-
Partial pressure of component i (Pa)
- \(P\) :
-
Pressure (Pa)
- \(P_{^\circ }\) :
-
Standard base pressure (Pa)
- \(r_{i}\) :
-
Rate of reaction of component i (kmol/h)
- \(R_{\text{g}}\) :
-
Universal constant of gases (kJ/mol k)
- \(A\) :
-
Cross sectional area (m2)
- \({\text{SV}}\) :
-
Space velocity (h−1)
- \(T\) :
-
Reaction temperature (K)
- \(T_{^\circ }\) :
-
Base reaction temperature (K)
- \(y_{i}\) :
-
Molar composition of component i (mol%)
- \(z\) :
-
Reactor height (m)
- \(\varepsilon\) :
-
Void fraction of catalyst bed
- \(\rho\) :
-
Density of gas mixture
- \(\rho_{\text{c}}\) :
-
Density of catalyst
- \(\mu\) :
-
Viscosity of gas mixture
References
Ahmed T (2001) Reservoir engineering handbook, 2nd edn. Elsevier, Amsterdam
Ancheyta J, Villafuerte-Macias E, Diaz-Garcia L, Gonzalez-Arredondo E (2001) Modeling and simulation of four catalytic reactors in series for naphtha reforming. Energy Fuels 15(4):887–893
Azizi N, Behbahani R, Isazadeh MA (2010) An efficient correlation for calculating compressibility factor of natural gases. J Nat Gas Chem 19(6):642–645
Bahadori A, Mokhatab S, Towler BF (2007) Rapidly estimating natural gas compressibility factor. J Nat Gas Chem 16(4):349–353
Bar-Meir G (2013) Fundamentals of compressible fluid_mechanics Version 0.3.3.0. Chicago, p557
Beggs DH, Brill JP (1973) A study of two-phase flow in inclined pipes. J Pet Technol 25(5):607–617
Elizalde I, Ancheyta J (2015) Dynamic modeling and simulation of a naphtha catalytic reforming reactor. Appl Math Model 39(2):764–775
Elsharkawy AM (2004) Efficient methods for calculations of compressibility, density and viscosity of natural gases. Fluid Phase Equilibria 218:1–13
Fayazi A, Arabloo M, Mohammadi AH (2014) Efficient estimation of natural gas compressibility factor using a rigorous method. J Nat Gas Sci Eng 16:8–17
George JA, Abdullah MA (2004) Catalytic naphtha reforming, 2nd edn. Marcel Dekker Inc., New York
Heidaryan E, Moghadasi J, Rahimi M (2010) New correlations to predict natural gas viscosity and compressibility factor. J Pet Sci Eng 73(1–2):67–72
Heidaryan E, Salarabadi A, Moghadasi J (2010) A novel correlation approach for prediction of natural gas compressibility factor. J Nat Gas Chem 19(2):189–192
John YM, Patel R, Mujtaba IM (2018) Effects of compressibility factor on fluid catalytic cracking unit riser hydrodynamics. Fuel 223:230–251
Kumar N (2004) Compressibility factors for natural and sour reservoir gases by correlations and cubic equations of state. Texas Tech University, Lubbock
Li C, Peng Y, Dong J (2014) Prediction of compressibility factor for gas condensate under a wide range of pressure conditions based on a three-parameter cubic equation of state. J Nat Gas Sci Eng 20:380–395
Mahmoud M (2014) Development of a new correlation of gas compressibility factor (Z-factor) for high pressure gas reservoirs. J Energy Resour. Technol. Trans ASME 136(1):012903–0129011
Papay JA (1968) A termelestechnologiai Parameterek valtozasa a gastelepek muve-lesesoran. OGIL Musz. Kuzl, Budapest: Tud
Riazi M (2005) Characterization and properties of petroleum fractions, 1st edn. International Standards (ASTM), Philadelphia
Rodríguez MA, Ancheyta J (2011) Detailed description of kinetic and reactor modeling for naphtha catalytic reforming. Fuel 90(12):3492–3508
Sanjari E, Lay EN (2012) Estimation of natural gas compressibility factors using artificial neural network approach. J Nat Gas Sci Eng 9:220–226
Shokir EME-M, El-Awad MN, Al-Quraishi AA, Al-Mahdy OA (2012) Compressibility factor model of sweet, sour, and condensate gases using genetic programming. Chem Eng Res Des 90(6):785–792
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Yusuf, A.Z., John, Y.M., Aderemi, B.O. et al. Effect of compressibility factor on the hydrodynamics of naphtha catalytic-reforming reactors. Appl Petrochem Res 9, 147–168 (2019). https://doi.org/10.1007/s13203-019-0233-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13203-019-0233-1