Olivier
Said-Aizpuru
*ab,
Florent
Allain
a,
Aurélie
Dandeu
*a,
Fabrice
Diehl
a,
David
Farrusseng
b and
Jean-François
Joly
a
aIFP Énergies nouvelles, Rond-point de l'échangeur de Solaize – BP 3, 69360 Solaize, France. E-mail: olivier.saidaizpuru@ugent.be
bIRCELYON, CNRS, 2 Avenue Albert Einstein, 69100 Villeurbanne, France
First published on 19th April 2021
Bridging the gap between kinetic model conception and catalyst design is targeted in catalytic naphtha reforming process development. New catalysts are continuously optimised in order to achieve higher selectivity in C5+ products. An adequate description of catalytic transformations by the kinetic models would therefore provide clues for catalyst design and accelerate the time to market implementation of process simulators. This study investigates the influence of site density and location on n-heptane reforming selectivity. It identifies the nature of the limiting steps for the different reforming pathways on a broad range of catalyst formulations. A common lumped model using power law kinetics is developed to describe already published experimental observations on the set of selected catalysts. Linear free energy relationships are used in order to handle a reduced number of statistically relevant adjustable parameters. The dependence between reference rate constants and active phase formulation is then unravelled. Unexpected results indicate that chlorine content and repartition at the crystallite scale affects the hydrogenolysis activity. Within the range of tested formulations, this study suggests that hydroisomerisation reactions are limited by acid sites transformations whereas the aromatisation pathways seem to proceed through a metal/acid bi-functional scheme. The further elaboration of a kinetic model that is able to predict the effect of an industrial catalyst active phase formulation change in full naphtha cut reforming lies beyond the scope of this article.
Heat and mass transport within catalyst pores, chemical transformation of reagents as well as active sites deactivation are examples of physico-chemical phenomena that are related to catalyst structure and that rule process performance. Given the complexity and the entanglement of these phenomena, the conception of kinetic models is often long and usually requires advanced characterisations of the catalyst. There is therefore an industrial need for robust kinetic modelling methodologies that would: (i) fasten the time to market development of simulators (ii) be able to predict the performances obtained with the use of a fictive solid in order to achieve retro-design of new catalysts.
Naphtha catalytic reforming would highly benefit from such a methodology.5 This process is aimed at producing high-octane gasoline bases, aromatics and hydrogen.6–8 Naphtha reforming chemistry requires metal/acid bi-functional catalysis.9 It consists in converting the linear paraffins and the naphthenes of a C6–C12 hydrocarbon cut into isoparaffins and aromatics. These targeted transformations require bi-functional catalytic reactions. Hydrocracking, hydrogenolysis and coking are parasite reactions. Industrial catalysts are made of platinum nano-particles (possibly associated with other metals such as tin or rhenium) that are supported on chlorinated gamma alumina. Metallic sites catalyse hydro/dehydrogenation reactions whereas chlorinated gamma alumina exhibits mild Brønsted acid sites that can be involved in H+ transfer.
Since the pioneering work of Smith in 1959,10 naphtha reforming kinetic models underwent a significant sophistication including a complexification of their reaction scheme and the introduction of deactivation kinetics. Detailed kinetic models based on power-law rate expressions are widespread approaches used in naphtha reforming. These models are used to obtain a good understanding of naphtha reforming reactions11 and are successfully applied to the conception of industrial reactors.12,13 A current challenge consists in finding a compromise between the simplicity and the accuracy of reforming detailed kinetic model.14
New reforming catalysts active phase formulations are continuously developed in order to reach higher C5+ compounds yields. Active phase optimisation is typically achieved by catalyst screening techniques followed by extrapolation to industrial conditions. This leads to a frequent renewal of the commercial catalysts offer and therefore requires a fast development of associated kinetic models. Reciprocally, current models are hardly able to predict the effect of active phase formulation changes on reforming kinetics and are still not used for the retro-design of new catalysts formulations.
A common industrial practice also consists in adjusting chlorine catalyst concentration (in the 0.7–1%wt range) during operation when a new feedstock is treated or when new product specifications are targeted. This practise still lacks dedicated models.
The objective of this study is to explore the impact of Pt/γ-Al2O3–Cl active sites density and repartition at the crystallite scale on n-heptane reforming intrinsic kinetics. A preliminary experimental study provided suitable data obtained in kinetic regime. Results have already been published and describe the effect of catalyst formulation over a broad range of mild n-heptane reforming conditions.15 The present article presents a kinetic model using power law kinetics and gathering reforming reactions into families through linear free energy relationships (LFER). This model describes the impact of an active phase formulation change on a limited number of apparent rate constants.
This study brings valuable information on the influence of site density and location over reforming selectivity. It also discusses the nature of limiting steps for the different reforming pathways on a broad range of catalyst formulations. This work focuses on intrinsic kinetic studies performed on a model reaction in order to simplify the reaction scheme. Operating conditions were chosen to lower the impact of intra-granular transport as well as deactivation phenomena. Consequently, the model that has been developed cannot be directly used for full naphtha cut reforming catalyst optimisation. Even if some experimental observation remain unexplained, this modelling study deepened the interpretation of experimental results by providing descriptors for the aromatisation and the hydroisomerisation pathways.
Platinum concentrations vary from 0.3%wt (close to the industrial catalyst content) to 1%wt. Platinum dispersion is high (>85%) and kept for all samples. A partition of the metal between sub-nanometre particles and single atoms was observed on fresh reduced catalyst.16 However, given the conservation of the dispersion level, metallic sites are considered as unique. Therefore, platinum content variation between catalysts is interpreted as a variation in the density of metallic sites.
Chlorine content ranges from 0.1%wt to 1.4%wt among the selected samples. This includes and largely extends the industrial reforming catalyst chlorine concentrations. The nature and the number of mild Brønsted acid sites is affected by chlorine concentration at the surface of the alumina crystallites. In this study, it has been chosen to consider a common mean Brønsted acid site for all the catalysts and to assume that chlorine content can somehow be correlated to the global acid sites concentration.
The 19 tested catalysts are split in two families regarding the nature of the γ-Al2O3 supports. 9 solids are prepared from Sasol® TH100™ boehmite gels. They present parallelepiped shaped crystallites that are referred to as T-flat. The 10 remaining catalysts are synthesised from Sasol® PuralSB3™ boehmite gels and exhibit more roundish alumina crystallites that are called P-egg. Compared to T-flat systems, P-egg crystallites present a higher ratio between edge length and facet surface for a given volume unit. Multi-technique structural investigations revealed that both chlorine and platinum particles are preferentially located on edges.16,17 Therefore, in this study, the repartition of platinum and chlorine at equal concentration is affected by the nature of the support.
Mass balances for each species on a control volume within the catalytic bed leads to a set of ordinary differential equations. Stationary regime is reached. This problem is then discretised and solved thanks to the LSODE program.19,20
Reactions are listed in Fig. 1 and on Table 1. They are grouped in 8 families. Hydrogenolysis and Hydrocracking reactions are assumed to be direct. Cyclisations, hydroisomerisations, naphthene interconversion and naphthene dehydrogenations are considered as reversible. Cyclisations are split in three families, C6 cyclisation, cyclisation into 6N7 and cyclisation into 5N7. It is to be noted that the proportion in 6N6 and A6 is rarely higher than 1% at the reactor outlet, which makes the model prediction insensitive to the presence of the benzene aromatisation network. Contrary to what is observed experimentally, the transformation between 6N7 and 5N7 is set at equilibrium in the kinetic model and the value of the corresponding equilibrium constant was fitted. We found it impossible to obtain both sensitive naphthene interconversion kinetic parameters and satisfactory fit of experimental results. This might be explained by a different reactivity of 5N7 isomers, that are not accurately quantified in these experiments.
Number | Reaction | Type | Rate expression |
---|---|---|---|
Family 1: hydrogenolysis | |||
1 | 6P6 + H2 → 2·P3 | Direct | r 1 = k1(P6P6PH2) |
2 | 6P6 + H2 → 4P4 + P2 | ′′ | r 2 = k2(P6P6PH2) |
3 (ref) | 7P7 + H2 → 6P6 + P1 | ′′ | r 3 = k 3 (P 7P7 P H2 ) |
4 | 7P7 + H2 → 5P5 + P2 | ′′ | r 4 = k4(P7P7PH2) |
5 | 7P7 + H2 → P3 + P4 | ′′ | r 5 = k5(P7P7PH2) |
6 | 2,5-P6 + H2 → 2·P3 | ′′ | r 6 = k6(P2,5-P6PH2) |
7 | 2-6P7 + H2 → 3P4 + P3 | ′′ | r 7 = k7(P2-6P7PH2) |
8 | 2-6P7 + H2 → 4P4 + P3 | ′′ | r 8 = k8(P2-6P7PH2) |
9 | 2,3-5P7 + H2 → 4P4 + P3 | ′′ | r 9 = k9(P2,3-5P7PH2) |
10 | 2,4-5P7 + H2 → 3P4 + P3 | ′′ | r 10 = k10(P2,4-5P7PH2) |
11 | 2-5P6 + H2 → 5P5 + P1 | ′′ | r 11 = k11(P2-5P6PH2) |
12 | 2-5P6 + H2 → 4P5 + P1 | ′′ | r 12 = k12(P2-5P6PH2) |
13 | 2-5P6 + H2 → 3P4 + P2 | ′′ | r 13 = k13(P2-5PCPH2) |
14 | 3-5P6 + H2 → 4P5 + P1 | ′′ | r 14 = k14(P3-5P6PH2) |
15 | 3-5P6 + H2 → 3P4 + P2 | ′′ | r 15 = k15(P3-5P6PH2) |
16 | 2-6P7 + H2 → 2-5P6 + P1 | ′′ | r 16 = k16(P2-6P7PH2) |
17 | 2-6P7 + H2 → 6P6 + P1 | ′′ | r 17 = k17(P2-6P7PH2) |
18 | 2-6P7 + H2 → 4P5 + P2 | ′′ | r 18 = k18(P2-6P7PH2) |
19 | 3-6P7 + H2 → 6P6 + P1 | ′′ | r 19 = k19(P3-6P7PH2) |
20 | 3-6P7 + H2 → 3-5P6 + P1 | ′′ | r 20 = k20(P3-6P7PH2) |
21 | 3-6P7 + H2 → 2-5P6 + P1 | ′′ | r 21 = k21(P3-6P7PH2) |
22 | 3-6P7 + H2 → 4P5 + P2 | ′′ | r 22 = k22(P3-6P7PH2) |
23 | 3-6P7 + H2 → 4P4 + 3P3 | ′′ | r 23 = k23(P3-6P7PH2) |
24 | 2,3-5P7 + H2 → 2-5P6 + P1 | ′′ | r 24 = k24(P2,3-5P7PH2) |
25 | 2,3-5P7 + H2 → 3-5P6 + P1 | ′′ | r 25 = k25(P2,3-5P7PH2) |
26 | 2,3-5P7 + H2 → 4P5 + P2 | ′′ | r 26 = k26(P2,3-5P7PH2) |
27 | 2,4-5P7 + H2 → 2-5P6 + P1 | ′′ | r 27 = k27(P2,4-5P7PH2) |
28 | Ethyl-5P7 + H2 → 3-5P6 + P1 | ′′ | r 28 = k28(Pethyl-5P7PH2) |
29 | Ethyl-5P7 + H2 → 5P5 + P2 | ′′ | r 29 = k29(Pethyl-5P7PH2) |
Family 2: hydrocracking | |||
---|---|---|---|
30 | 6P6 + H2 → 2·P3 | Direct | r 30 = k30(P6P6PH2) |
31 (ref) | 7P7 + H2 → 4P4+P3 | ′′ | r 31 = k 31 (P 7P7 P H2 ) |
32 | 2-5P6 + H2 → 2·P3 | ′′ | r 32 = k32(P2-5P7PH2) |
33 | 2-6P7 + H2 → 3P4 + P3 | ′′ | r 33 = k33(P2-6P7PH2) |
34 | 2-6P7 + H2 → 4P4 + P3 | ′′ | r 34 = k34(P2-6P7PH2) |
35 | 3-6P7 + H2 → 4P4 + P3 | ′′ | r 35 = k35(P3-6P7PH2) |
36 | 2,3-5P7 + H2 → 4P4 + P3 | ′′ | r 36 = k36(P2,3-5P7PH2) |
37 | 2,4-5P7 + H2 → 3P4 + P3 | ′′ | r 37 = k37(P2,4-5P7PH2) |
Family 3: 6N6 cyclisation | |||
---|---|---|---|
38 | 6N6 + H2 ⇄ 6P6 | Reversible | r 38 = k38(P6P6PH2 − P6P6/Keq (38)) |
39 (ref) | 5N6 + H2 ⇄ 2-5P6 | ′′ | r 39 = k 39 (P 5N6 P H2 − P 2-5P6 /K eq (39)) |
40 | 5N6 + H2 ⇄ 3-5P6 | ′′ | r 40 = k40(P5N6PH2 − P3-5P6/Keq (40)) |
41 | 5N6 + H2 ⇄ 6P6 | ′′ | r 41 = k41(P5N6PH2 − P6P6/Keq (41)) |
Family 4: 6N7 cyclisation | |||
---|---|---|---|
42 (ref) | 6N7 + H2 ⇄ 7P7 | Reversible | r 42 = k 42 (P 6N7 P H2 − P 7P7 /K eq (42)) |
43 | 6N7 + H2 ⇄ 2-6P7 | ′′ | r 43 = k43(P6N7PH2 − P2-6P7/Keq (43)) |
44 | 6N7 + H2 ⇄ 3-6P7 | ′′ | r 44 = k44(P6N7PH2 − P3-6P7/Keq (44)) |
Family 5: 5N7 cyclisation | |||
---|---|---|---|
45 (ref) | 5N7 + H2 ⇄ 2-6P7 | Reversible | r 45 = k 45 (P 5N7 P H2 − P 2-6P7 /K eq (45)) |
46 | 5N7 + H2 ⇄ 3-6P7 | ′′ | r 46 = k46(P5N7PH2 − P3-6P7/Keq (46)) |
47 | 5N7 + H2 ⇄ 2,3-5P7 | ′′ | r 47 = k47(P5N7PH2 − P2,3-5P7/Keq (47)) |
48 | 5N7 + H2 ⇄ 2,4-5P7 | ′′ | r 48 = k48(P5N7PH2 − P2,4-5P7/Keq (48)) |
Family 6: hydroisomerisation | |||
---|---|---|---|
49 | 5P5 ⇄ 4P5 | Reversible | r 49 = k49(P5P5 − P4P5/Keq (49)) |
50 | 6P6 ⇄ 2-5P6 | ′′ | r 50 = k50(P6P6 − P2-5P6/Keq (50)) |
51 | 2-5P6 ⇄ 3-5P6 | ′′ | r 51 = k51(P2-5P6 − P3-5P6/Keq (51)) |
52 | 6P6 ⇄ 3-5P6 | ′′ | r 52 = k52(P6P6 − P3-5P6/Keq (52)) |
53 | 7P7 ⇄ 2-6P6 | ′′ | r 53 = k53(P7P7 − P2-6P7/Keq (53)) |
54 (ref) | 7P7 ⇄ 3-6P7 | ′′ | r 54 = k 54 (P 7P7 − P3-6P7/Keq(54)) |
55 | 7P7 ⇄ ethyl-5P7 | ′′ | r 55 = k55(P7P7 − Pethyl-5P7/Keq (55)) |
56 | 2-6P7 ⇄ 3-6P7 | ′′ | r 56 = k56(P2-6P7 − P2,3-5P7/Keq (56)) |
57 | 2-6P7 ⇄ 2,3-5P7 | ′′ | r 57 = k57(P2-6P7 − P2,3-5P7/Keq (57)) |
58 | 2-6P7 ⇄ 2,4-5P7 | ′′ | r 58 = k58(P2-6P7 − P2,4-5P7/Keq (58)) |
Power law kinetics are used and partial orders are equal to the stoichiometric coefficient of the reagents. The expressions handle partial pressures in gas phase and equilibrium constants. Rate expressions are provided on Table 1. Excepted the reactions of methylcyclohexane dehydrogenation and naphthenes interconversion, the equilibrium constant of reaction j(Keq(j)) is given by .
Methylcyclohexane dehydrogenation experiments were carried out on the same set-up in order to determine the value of the related equilibrium constant. Results indicated strong deviations between the experimental dehydrogenation equilibrium and the expected value considering the ideal gas tabulated free enthalpy. This observation has been reported by other authors.21,22 The correlation proposed by Schildhauer et al. was well verified on methylcyclohexane dehydrogenation test and was therefore implemented in the present n-heptane reforming kinetic model. This gives, with T expressed in Kelvins, . A dihydrogen partial pressure order was introduced and fitted in the expression of methylcyclohexane dehydrogenation rate constant. The value of this order was found to be close to 0.3. This suggests that hydrogen adsorption on the catalytic surface cannot be neglected and that a power law expression is hardly valid for this equation.
The teams of Yoneda and Klein made a similar use of LFER in order to address catalyst formulation change issues.23–29 They observed that, when ln(ki) is plotted over Δ(ΔrH) the strait lines corresponding to an identical family of reactions for different catalysts are parallel (see Fig. 2). This suggests that the proportionality coefficient αj is kept from a catalyst to another. This is an approximation that was made all along this study.
Reference activation energies, reference rate constants and αj parameters were first fitted on P-egg tests. Then, activation energies and αj coefficients values were extrapolated to T-flat catalysts before the adjustment of reference rate constants on this second family of catalysts. N6 cyclisation (family 3) parameters were not statistically significant. N6 compounds are trace species and family 3 could have been removed from the model.
The optimisation routine is written as a Fortran code and is based on a modified Levenberg–Marquardt algorithm.30,31 The optimiser identifies a set of parameters values that minimises an objective function measuring the deviations between model-predicted and experimental observables. In this case, the outlet molar flows of the 25 selected hydrocarbons are the observables. As presented on Table 2, the objective function is defined as the sum of square deviations for each observable weighted by the carbon atom number of each species. The experiments selected for the optimisation procedure were randomly chosen at a total pressure of 10 barg. The number of selected points is equal to 75% of the validation basis (see Table 3). The set of common fitted parameters estimates is displayed on Table 4.
Fit on P-egg catalysts | Fit on T-flat catalysts | |
---|---|---|
Number of fitted parameters | 82 | 63 |
Validation basis size/pts | 405 | 355 |
Fitting basis size/pts | 622 | 551 |
Objective function | 2.8 × 10−4 | 3.5 × 10−4 |
Fig. 5 Model prediction of temperature and contact time variations for a given catalyst (P-egg_0.3%wtPt_1.4%wtCl). |
The individual significance of each parameter was assessed by performing a student's t-test. In this study, the parameters were considered statistically different from zero, when the estimates dispersion exceeds 95% of the reference tabulated value at the same degree of freedom. Parameters for which the absolute t-value is lower than ten do not satisfy the t-test. 6N cyclisation kinetic parameters, hydrocracking and 6N7 cyclisation LFER α coefficients are therefore considered as statistically irrelevant whereas all the other parameters satisfy the t-test.
Hydrogenolysis and naphthene dehydrogenation pathways were expected to be catalysed by the metallic phase only. As illustrated on Fig. 8, at a given Cl content, hydrogenolysis activity is roughly proportional to Pt concentration. On the contrary, an increase in Pt is not correlated to a significant rise in naphthene dehydrogenation reference rate constant. This could be explained by the rapidity of naphthene dehydrogenation transformations. Even at 0.3%wt Pt, cyclisation appears as the limiting paraffin aromatisation step. Similarly, the apparent dependence of dehydrogenation reference constant on Cl content might be an artefact inherited from the preceding cyclisation step. The proper distinction between naphthene cyclisation and dehydrogenation kinetics is made difficult by the low concentrations in cycloalkane intermediates. However, a clear dependence of the overall paraffin aromatisation pathway on chlorine content is observed. This supports the bi-functional description of the aromatisation pathway in naphtha reforming.32–35
Hydroisomerisation activity is controlled by chlorine catalyst content and there is no significant dependence of this parameter on Pt content. Moreover, hydroisomerisation activity tends to disappear when chlorine concentration tends to zero. Therefore, it seems that in the range of formulations tested in this study, hydroisomerisation is limited by transformation on acid sites.
Direct comparison between T-flat and P-egg catalysts is provided on Fig. 7. The striking difference between the catalytic behaviour of the two families of catalysts lies in the hydrogenolysis activity, hydrogenolysis being enhanced on T-flat catalysts. The impact of the support nature on the kinetic parameters of the other reaction families is small or negligible.
The selectivity in C1 + C2 products singularly drops at intermediate chlorine content (around 0.5%wt). According to this study, this phenomenon is specifically due to an inhibition of the hydrogenolysis rate constants. This observation was reproducible and systematically made on different supports and at different Pt contents (see Fig. 8). The origin of the non-monotonous trends linking catalysts performance and their active phase composition is still unclear. A more accurate description of chlorine impact on active phase structure and a more sophisticated description of the active sites could bring some explanations. Chlorine was found to stabilise Pt single atoms that can possibly act as secondary metallic sites. As in the case of the support effect, the influence of chlorine on metallic site structure and strength may be involved in the non-monotonous relation between hydrogenolysis activity and chlorine content. Another explanation is proposed: assuming the formation of a coke deposit on acid sites, the preferential location of chlorine could impact the preferential repartition of coke and could lead to an inhibition of hydrogenolysis activity at intermediate chlorine levels.
This study digs two debated issues in naphtha reforming with Pt/γ-Al2O3–Cl catalysts: the nature of the catalytic limiting steps for different reaction pathways as well as the physico-chemical properties of the catalyst impacting reaction selectivity. As observed experimentally, this study confirms that the main impact of alumina support change is made on the hydrogenolysis activity. The rate constants related to the other reaction families are poorly affected by the nature of the support. Hydroisomerisation activity is very preferentially correlated to chlorine content. Over the catalyst panel, it is therefore likely that hydroisomerisation reactions proceed through acid step limited mechanisms.
Aromatisation of paraffins through dehydrocyclisation into naphthenes followed by dehydrogenation was found to be sensitive on both platinum and chlorine amounts. This result excludes an aromatisation pathway catalysed exclusively on the metallic sites.
As reported by the experimental study, a spectacular inhibition of the hydrogenolysis activity is systematically observed at intermediate chlorine contents (around 0.5%wt Cl). The kinetic model does not provide explanations about the origin of this phenomenon. Complementary testing and investigations of post-mortem catalysts could bring some hints.
One can suggest to apply the kinetic modelling methodology proposed in n-heptane reforming to a full naphtha cut. It requires the extension and the validation of LFER reaction families to heavier compounds (up to C12). The addition of a deactivation model is also needed in these conditions. Then, a descriptor based model could be extrapolated at the process scale. Some conclusions can still already be applied to the optimisation of new catalysts. The dependence of the hydrogenolysis activity on chlorine content and gamma alumina support morphology opens new catalysts development perspectives in order to meet the demand for higher C5+ reformate products yields.
α | LFER proportionality coefficient |
ΔrG | Reaction free enthalpy |
ΔrH | Reaction enthalpy |
ω | Parameter weight in the regression procedure |
E a | Activation energy |
exp | Experimental value |
F n | Molar flow of species n |
h | Catalyst |
i | Family of reaction |
j | Reaction |
K eq | Equilibrium constant |
k | Rate constant |
n | Reaction partial order |
obs | Observable |
P | Partial pressure |
o | Reference reaction within a family |
R | Ideal gas constant |
S | Selectivity |
sim | Simulated value |
X | Conversion |
Y | Yield |
v | Speed of reaction |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1re00073j |
This journal is © The Royal Society of Chemistry 2021 |