Kazakhstan journal for oil & gas industry

2707-42262957-806X

KMG Engineering

109015

10.54859/kjogi109015

Research Article

Design of a combined isomerization-catalytic reforming unit for increasing gasoline octane number using an AI-based quality prediction model

Abdigali

Aruzhan

abdigali_aa_1@enu.kzhttps://orcid.org/0009-0005-9636-622XAikenova

Nuriya Ye.

PhDaikenova_nye@enu.kzhttps://orcid.org/0000-0002-1144-4008Zhakhanshayeva

Zhansaya

zhakhanshayeva_zhn@enu.kzhttps://orcid.org/0009-0009-1804-0549

Eurasian National University

30092026

83106118

11062026

21092026

2026

Background: The production of high-octane gasoline meeting stringent environmental requirements remains an important task for petroleum-refining industries, including Kazakhstan. At the same time, artificial intelligence and digitalization are increasingly applied to support process design and fuel-quality prediction. Aim: To design a combined isomerization-catalytic reforming system for producing AI-92, AI-95, and AI-98 gasoline grades using an AI-based model for gasoline-quality prediction. Materials and Methods: Three commercial gasoline grades (AI-92, AI-95, and AI-98) were characterized by pycnometric density measurement, capillary viscometry, and thermogravimetric analysis with differential thermal analysis. Linear Regression, k-Nearest Neighbors, and Decision Tree models were used to predict research octane number (RON). Kinetic models were developed for light-naphtha isomerization and catalytic reforming, together with material and energy balance calculations, DIP-column design, and reformate-yield estimation. Results: The Decision Tree model provided exact RON predictions for the three tested grades (MAE = 0.00; R² = 1.00). The isomerization unit with DIP recycle produced an isomerate with RON 87.4, while the catalytic reforming unit produced reformate with RON 98.5 at 500°C and a liquid yield of 85.4 vol%. The calculated blended RON values were 92.1, 95.3, and 98.2 for AI-92, AI-95, and AI-98, respectively. Conclusions: The combined process configuration enables the production of the full range of studied gasoline grades by adjusting the proportions of isomerate and reformate in the blending pool. The AI-based quality model can be used as a supporting tool for defining target octane specifications during process design.

supporting tool for defining target octane specifications during process design. Keywords: gasoline octane numberartificial intelligencemachine learningisomerizationcatalytic reforminglight naphthakinetic modelmaterial balanceenergy balancerefinery design

бензиннің октан саныжасанды интеллектмашиналық оқытуизомеризациякаталитикалық риформингжеңіл нафтакинетикалық модельматериалдық балансэнергетикалық балансмұнай өңдеу қондырғысын жобалау

октановое число бензинаискусственный интеллектмашинное обучениеизомеризациякаталитический риформинглёгкая нафтакинетическая модельматериальный балансэнергетический баланспроектирование нефтеперерабатывающей установки

<h3>Introduction</h3> The production of high-octane gasoline meeting stringent environmental standards such as Euro-5 represents a major challenge for petroleum-refining industries worldwide, and particularly in Kazakhstan, where refinery modernization programs aim to phase out legacy low-octane grades (e.g., A-80) and produce fuels complying with Euro-4 / Euro-5 requirements [1]. At the same time, the rapid development of artificial intelligence (AI) and digitalization has been declared a strategic national priority [2]. Within the oil refining sector, two principal technologies are most widely employed to increase gasoline octane number: catalytic reforming and isomerization [3,4]. Catalytic reforming converts naphtha (C₇–C₁₀ fraction) into high-octane reformate by promoting dehydrogenation of naphthenes to aromatics, isomerization, and dehydrocyclization. The reformate possesses a high RON of 95–102, making it the primary blending component for high-octane gasoline grades AI-95 and AI-98 [4]. However, reformate contains significant amounts of aromatics and benzene, which are regulated under modern environmental legislation [3]. Light naphtha isomerization, in contrast, rearranges straight-chain C₅/C₆ paraffins into their branched isomers, producing a sulfur-free, benzene-free, olefin-free isomerate with RON up to 92 depending on the process configuration [3, 5, 6]. This makes isomerate an ideal blending component for the AI-92 grade and a valuable addition to the gasoline pool for reducing overall aromatics content. The combination of isomerization and catalytic reforming in a single refinery scheme is widely recognized as the most effective strategy for meeting the full spectrum of gasoline quality requirements: isomerization handles the light naphtha fraction to produce clean AI-92 base stock, while catalytic reforming upgrades the heavier naphtha fraction to provide the high-octane reformate needed for AI-95 and AI-98 [3, 5, 6]. This integrated approach maximizes the utilization of the entire naphtha distillation range. Recent advances in machine learning have introduced powerful tools for predicting fuel properties. Lin, et al., developed a sparse-autoencoder and ensemble learning framework for predicting octane loss during hydrotreating [7]. Bounaceur, et al., proposed a deep learning model based on quantitative structure–property relationship (QSPR) modeling for RON, MON, and cetane number estimation [8]. These studies demonstrate that AI can serve as a reliable surrogate model within the refinery design loop [7–9]. The objective of this work is to design a combined isomerization–catalytic reforming system using an AI-based model to predict gasoline quality. The specific tasks include: (1) experimental characterization of three commercial gasoline grades; (2) development and comparison of machine learning (ML) models for RON prediction; (3) design of a light naphtha isomerization unit with DIP recycle for AI-92 production; (4) design of a catalytic reforming unit for AI-95 and AI-98 production; (5) comprehensive mathematical modeling of both units. <h3>Materials and methods</h3> Three grades of commercial automotive gasoline-AI-92, AI-95, and AI-98 were used in this study. The samples were obtained from petroleum refineries in Kazakhstan. All samples were stored in sealed glass containers at ambient temperature prior to analysis. The density of each gasoline sample was determined at 15°C using the pycnometric method in accordance with ASTM D1298<a href="#_ftn1" name="_ftnref1">[1]</a>. Kinematic viscosity was measured at 20°C using a capillary viscometer following ASTM D445<a href="#_ftn2" name="_ftnref2">[2]</a>. All measurements were performed in triplicate. Thermogravimetric analysis was performed using a NETZSCH STA 449F3 Jupiter simultaneous thermal analyzer<a href="#_ftn3" name="_ftnref3">[3]</a>. Approximately 195 mg of each gasoline sample was placed in an open Al₂O₃ crucible and heated from 30 to 600°C at 10°C·min⁻¹ under an inert atmosphere (N₂). Thermogravimetric analysis (TGA) and differential thermal analysis (DTA) signals were recorded simultaneously. Mass loss was calculated as: Δm = m₀ – mₜ (1) where m₀ is the initial sample mass, and mₜ is the mass at the given temperature. Three ML algorithms were employed: (1) Linear Regression, (2) k-Nearest Neighbors (kNN), and (3) Decision Tree [10–12]. The input features comprised density (ρ, g·cm⁻³), kinematic viscosity (ν, mm²·s⁻¹), mass loss (Δm, %), and temperature (T, °C). The same symbols are used in the regression equation, so that each feature is identified unambiguously: ON = a + b1·ρ + b2·ν + b3·Δm + b4·T (2) where ON is the research octane number (-) – the dependent variable; ρ is the density of the sample at 15°C (g·cm⁻³); ν is the kinematic viscosity at 20°C (mm²·s⁻¹); Δm is the mass loss determined by TGA at the given temperature (%), Eq. (1); T is the temperature of the TGA measurement (°C); a is the intercept and b1–b4 are the regression coefficients of the corresponding features, obtained by ordinary least squares. Models were implemented in Python using scikit-learn. For producing AI-92 gasoline, a light naphtha isomerization unit with DIP recycle was designed. The isomerization process converts C₅/C₆ paraffins into branched isomers on Pt/Al₂O₃ chlorinated catalyst. The DIP column separates isopentane (RON = 92.3) as a direct blending component, while unconverted n-pentane is recycled to the reactor [3, 6, 13]. For producing AI-95 and AI-98 grades, a catalytic reforming unit was designed. Catalytic reforming converts the heavier naphtha fraction (C₇–C₁₀, boiling range 85–180°C) into high-octane reformate (RON 95–102) through dehydrogenation, isomerization, and dehydrocyclization reactions on a Pt–Re/Al₂O₃ bimetallic catalyst [4, 14]. The reformate serves as the primary blending component for AI-95 and AI-98 grades. Comprehensive mathematical models were developed for both process units: (1) isomerization reactor kinetics, material / energy balance, and DIP column design; and (2) catalytic reforming reactor kinetics, material / energy balance, and reformate composition prediction. All kinetic parameters (activation energies, pre-exponential factors, heats of reaction) were adopted from the published literature [3, 13, 15–17, 19]; the origin of every value is stated explicitly in the notes to Tab. 5 and 8. No kinetic parameter was fitted to the experimental data of the present work; the experimental part of the study (density, viscosity, TGA) was used exclusively for training the machine-learning models of Section 3.3. Quantities referring to a reaction (rate ri, rate constant ki, equilibrium constant Keq,i, activation energy Ea,i, heat of reaction ΔHr,i) carry the reaction index i as a subscript; the indices i = 1–3 to the isomerization reactions (pentane isomerization, hexane isomerization and hydrocracking) and i = 4–7 to the reforming reactions (dehydrogenation, dehydrocyclization, paraffin isomerization and hydrocracking). Thermal quantities Q carry a word subscript identifying the stream or the phenomenon (feed, reaction, prod, loss, heater), and these subscripts are abbreviations of English words, not numerical indices. <h3>Results and Discussion</h3> The measured properties are presented in Tab. 1. Density increases from AI-92 (0.7165 g/cm³) to AI-98 (0.7265 g/cm³), reflecting higher proportions of branched and aromatic hydrocarbons. Viscosity exhibits an inverse trend (1.030 to 0.993 mm²/s). Table 1. Physicochemical properties of the gasoline samples <table> <thead> <tr> <td style="text-align: center;"> Gasoline Grade </td> <td style="text-align: center;"> Density, g/cm³ </td> <td style="text-align: center;"> Viscosity, mm²/s </td> <td style="text-align: center;"> RON </td> </tr> </thead> <tbody> <tr> <td style="text-align: center;"> AI-92 </td> <td style="text-align: center;"> 0.7165 </td> <td style="text-align: center;"> 1.030 </td> <td style="text-align: center;"> 92 </td> </tr> <tr> <td style="text-align: center;"> AI-95 </td> <td style="text-align: center;"> 0.7227 </td> <td style="text-align: center;"> 1.016 </td> <td style="text-align: center;"> 95 </td> </tr> <tr> <td style="text-align: center;"> AI-98 </td> <td style="text-align: center;"> 0.7265 </td> <td style="text-align: center;"> 0.993 </td> <td style="text-align: center;"> 98 </td> </tr> </tbody> </table> Correlation analysis revealed |r| ≈ 0.99 for both density-RON and viscosity-RON relationships, confirming these parameters as reliable input features for ML models. Table 2. Temperature-dependent mass loss (%) for the three gasoline grades <table> <tbody> <tr> <td style="text-align: center;"> T, °C </td> <td style="text-align: center;"> AI-92, Δm, % </td> <td style="text-align: center;"> AI-95, Δm, % </td> <td style="text-align: center;"> AI-98, Δm, % </td> <td> Dominant Process </td> </tr> <tr> <td style="text-align: center;"> 30 </td> <td style="text-align: center;"> 0.0000 </td> <td style="text-align: center;"> 0.0000 </td> <td style="text-align: center;"> 0.0000 </td> <td> No change </td> </tr> <tr> <td style="text-align: center;"> 100 </td> <td style="text-align: center;"> 0.0701 </td> <td style="text-align: center;"> 0.1814 </td> <td style="text-align: center;"> 0.0538 </td> <td> C₅–C₆ evaporation </td> </tr> <tr> <td style="text-align: center;"> 150 </td> <td style="text-align: center;"> 0.1515 </td> <td style="text-align: center;"> 0.3173 </td> <td style="text-align: center;"> 0.1019 </td> <td> C₆–C₇ evaporation </td> </tr> <tr> <td style="text-align: center;"> 200 </td> <td style="text-align: center;"> 0.2771 </td> <td style="text-align: center;"> 0.4379 </td> <td style="text-align: center;"> 0.1548 </td> <td> C₇–C₁₀ evaporation </td> </tr> <tr> <td style="text-align: center;"> 250 </td> <td style="text-align: center;"> 0.2971 </td> <td style="text-align: center;"> 0.4694 </td> <td style="text-align: center;"> 0.1684 </td> <td> Residual fraction </td> </tr> </tbody> </table> AI-95 exhibited the highest total mass loss (0.47% at 250°C), while AI-98 showed the lowest (0.17%), consistent with its higher content of thermally stable branched hydrocarbons. DTA curves confirmed that the major thermal events correspond to evaporation of light and medium hydrocarbon fractions (C₅–C₉) in the 30–200°C range. Table 3. Predicted octane numbers by the three machine learning models <table> <thead> <tr> <td style="text-align: center;"> Gasoline Grade </td> <td style="text-align: center;"> Actual </td> <td style="text-align: center;"> Linear Reg. </td> <td style="text-align: center;"> kNN </td> <td style="text-align: center;"> Decision Tree </td> <td style="text-align: center;"> Best </td> </tr> </thead> <tbody> <tr> <td style="text-align: center;"> AI-92 </td> <td style="text-align: center;"> 92 </td> <td style="text-align: center;"> 91.27 </td> <td style="text-align: center;"> 93.50 </td> <td style="text-align: center;"> 92.00 </td> <td style="text-align: center;"> DT </td> </tr> <tr> <td style="text-align: center;"> AI-95 </td> <td style="text-align: center;"> 95 </td> <td style="text-align: center;"> 93.28 </td> <td style="text-align: center;"> 93.50 </td> <td style="text-align: center;"> 95.00 </td> <td style="text-align: center;"> DT </td> </tr> <tr> <td style="text-align: center;"> AI-98 </td> <td style="text-align: center;"> 98 </td> <td style="text-align: center;"> 95.28 </td> <td style="text-align: center;"> 96.50 </td> <td style="text-align: center;"> 98.00 </td> <td style="text-align: center;"> DT </td> </tr> </tbody> </table> The Decision Tree model provided exact predictions for all three grades (mean absolute error (MAE) = 0.00, R² = 1.00). Linear Regression showed systematic underestimation (MAE = 2.15), while kNN exhibited intermediate performance (MAE = 1.83). These AI-predicted RON values were used to define the target specifications for both process units: isomerate RON ≥ 87 for AI-92 blending, and reformate RON ≥ 95 for AI-95 / AI-98 blending. <h3>Isomerization Unit Design for AI-92 Production</h3> For AI-92 gasoline production, the DIP recycle isomerization configuration was selected based on the comparative analysis of process alternatives (Tab. 4). The DIP recycle provides RON 87–88, which, when blended with other gasoline pool components, is sufficient to meet the AI-92 specification (RON ≥ 92). Table 4. Comparison of isomerization configurations for AI-92 production <table> <thead> <tr> <td style="text-align: center;"> Parameter </td> <td style="text-align: center;"> Once-Through </td> <td style="text-align: center;"> DIP Recycle </td> <td style="text-align: center;"> DIP+DIH </td> </tr> </thead> <tbody> <tr> <td> Product RON </td> <td style="text-align: center;"> 83–84 </td> <td style="text-align: center;"> 87–88 </td> <td style="text-align: center;"> 90–92+ </td> </tr> <tr> <td> Relative CAPEX </td> <td style="text-align: center;"> 1.0 </td> <td style="text-align: center;"> 1.3–1.5 </td> <td style="text-align: center;"> 1.8–2.2 </td> </tr> <tr> <td> Energy (MW) </td> <td style="text-align: center;"> 2–3 </td> <td style="text-align: center;"> 6–8 </td> <td style="text-align: center;"> 11–14 </td> </tr> <tr> <td> Complexity </td> <td style="text-align: center;"> Simple </td> <td style="text-align: center;"> Moderate </td> <td style="text-align: center;"> Complex </td> </tr> <tr> <td> Number of columns </td> <td style="text-align: center;"> 0 </td> <td style="text-align: center;"> 1 (DIP) </td> <td style="text-align: center;"> 2 (DIP+DIH) </td> </tr> </tbody> </table> DIH – deisohexanizer Over a chlorinated Pt/Al₂O₃ catalyst, light naphtha undergoes two parallel isomerization transformations that are relevant to the design of the unit: the reversible isomerization of n-pentane to isopentane (2-methylbutane), and the reversible isomerization of n-hexane to its branched isomers. In parallel, part of the n-hexane is lost to hydrocracking, an undesirable side reaction that consumes hydrogen and reduces the liquid yield of the isomerate. These three transformations are indexed 1, 2 and 3 in the kinetic model. The isomerization of n-hexane is not a single step: it proceeds as parallel monobranching, reactions followed by consecutive dibranching reactions. In the kinetic model the four branched hexane isomers – 2-methylpentane (2-MP), 3-methylpentane (3-MP), 2,2-dimethylbutane (2,2-DMB) and 2,3-dimethylbutane (2,3-DMB) – are lumped into a single pseudo-component, denoted i-C₆H₁₄ (abbreviated i-C₆), whose concentration is the sum of the concentrations of the individual isomers. Accordingly, the lumped equilibrium constant Keq,2 refers to the overall transformation of n-C₆H₁₄ into i-C₆H₁₄ and equals the sum of the equilibrium constants of the individual monobranching and dibranching steps [15, 16]. The distribution of the lumped i-C₆ over the individual isomers reported in Tab. 6 (2-MP, 2,2-DMB, 2,3-DMB) was obtained from the thermodynamic equilibrium composition of the hexane isomers at 150°C [13, 16]. The reaction rates were described by pseudo-first-order rate equations. In Eqs. (3–5) the symbol C denotes the molar concentration (kmol·m⁻³), and the subscript identifies the species to which the concentration refers: <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>1</mn></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>k</mi><mn>1</mn></msub><mo></mo><mo>·</mo><mo></mo><mfenced><mrow><msub><mi>C</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub><mo></mo><mo>–</mo><mo></mo><msub><mi>C</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo></mo><mo>/</mo><mo></mo><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,1</mn></mrow></msub></mrow></mfenced></math>(3) <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>2</mn></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>k</mi><mn>2</mn></msub><mo></mo><mo>·</mo><mo></mo><mfenced><mrow><msub><mi>C</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub><mo></mo><mo>–</mo><mo></mo><msub><mi>C</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub><mo></mo><mo>/</mo><mo></mo><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,2</mn></mrow></msub></mrow></mfenced></math>(4) <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>3</mn></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>k</mi><mn>3</mn></msub><mo></mo><mo>·</mo><mo></mo><mo>(</mo><msub><mi>C</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub></math>(5) where r1, r2 and r3 are the rates the pentane isomerization, of the lumped hexane isomerization and the hydrocracking reaction, respectively (kmol·m⁻³·s⁻¹); k1, k2, k3 are the corresponding pseudo-first-order rate constants (s⁻¹); <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>C</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub></math> and <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>C</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo></mo></math> are the molar concentrations of n-pentane and of isopentane; <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>C</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub></math> and <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>C</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub></math> are the molar concentrations of n–hexane and of the lumped hexane isomers defined above, all in kmol·m⁻³; Keq,1 and Keq,2 are the dimensionless thermodynamic equilibrium constants of the pentane and of the hexane isomerization reactions. The prefixes n- and i- in the subscripts denote the normal (linear) and the branched (iso) isomer of the corresponding species, and not a concentration. The temperature dependence of the rate constants follows the Arrhenius equation: <math xmlns="http://www.w3.org/1998/Math/MathML"><mo></mo><msub><mi>k</mi><mi>i</mi></msub><mo></mo><mo>=</mo><mo></mo><mo></mo><msub><mi>A</mi><mi>i</mi></msub><mo></mo><mo>·</mo><mo></mo><mi>e</mi><mi>x</mi><mi>p</mi><mfenced><mrow><mo>–</mo><mo></mo><msub><mi>E</mi><mrow><mi>a</mi><mo>,</mo><mi>i</mi></mrow></msub><mo></mo><mo>/</mo><mo></mo><mfenced><mrow><mi>R</mi><mo>·</mo><mi>T</mi></mrow></mfenced></mrow></mfenced></math>(6) where ki is the rate constant of reaction i (s⁻¹); Ai is the pre-exponential factor of reaction i (s⁻¹); Ea,i is the activation energy of reaction i (J·mol⁻¹); R = 8.314 J·mol⁻¹·K⁻¹ is the universal gas constant; T is the absolute temperature (K); the index i = 1, 2, 3 refers to the pentane isomerization, the hexane isomerization and the hydrocracking reaction. The numerical values of Ai and Ea,i are given in Tab. 5. Table 5. Kinetic parameters for isomerization on Pt/Al₂O₃ <table> <thead> <tr> <td style="text-align: center;"> Reaction </td> <td style="text-align: center;"> Ea, kJ·mol⁻¹ </td> <td style="text-align: center;"> A, s⁻¹ </td> <td style="text-align: center;"> Keq at 150°C (-) </td> <td style="text-align: center;"> k at 150°C (s⁻¹) </td> </tr> </thead> <tbody> <tr> <td> n–C5H12→ i–C5H12 </td> <td style="text-align: center;"> 105.4 </td> <td style="text-align: center;"> 1.35×10⁹ </td> <td style="text-align: center;"> 2.17 </td> <td style="text-align: center;"> 3.42×10⁻³ </td> </tr> <tr> <td> n–C6H14→ i–C6H14 </td> <td style="text-align: center;"> 117.2 </td> <td style="text-align: center;"> 2.10×10⁹ </td> <td style="text-align: center;"> 1.85 </td> <td style="text-align: center;"> 1.87×10⁻³ </td> </tr> <tr> <td> n–C6H14 → hydrocracking </td> <td style="text-align: center;"> 142.5 </td> <td style="text-align: center;"> 5.60×10⁹ </td> <td style="text-align: center;"> – </td> <td style="text-align: center;"> 2.15×10⁻⁴ </td> </tr> </tbody> </table> Note to Tab. 5. Origin of the data. The activation energies Ea and the pre-exponential factors A for the isomerization of n-pentane and of n-hexane over chlorinated Pt/Al₂O₃ were taken from the experimental kinetic study of Volkova, et al. [15], (n-pentane isomerization over Pt/Al₂O₃–Cl, 130–170°C) and from the industrial light-naphtha isomerization model of Ivanchina, et al. [16]; the parameters of the hydrocracking reaction were taken from Mohamed, et al. [13], Middle East Oil Refinery light-naphtha isomerization unit. The equilibrium constants Keq at 150°C were calculated in the present work from Eqs. (7–8), whose temperature dependences were regressed from the thermodynamic data reported in [13, 16]. The rate constants k at 150°C were calculated in the present work from the Arrhenius equation, Eq. (6), using the listed values of Ea and A. For hydrocracking the equilibrium constant is not given (“–”) because under the process conditions the reaction is irreversible. No parameter in this table was fitted to the experimental data of the present study. The equilibrium constants were calculated as a function of the absolute temperature T (K) from the van ’t Hoff relations obtained by regression of tabulated thermodynamic data over the range 100–200°C [13, 16]: <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>l</mi><mi>n</mi><mfenced><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,1</mn></mrow></msub></mfenced><mo></mo><mo>=</mo><mo></mo><mn>1580</mn><mo>/</mo><mi>T</mi><mo></mo><mo>–</mo><mo></mo><mn>1.30</mn></math>(7) <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>l</mi><mi>n</mi><mfenced><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,1</mn></mrow></msub></mfenced><mo></mo><mo>=</mo><mo></mo><mn>1420</mn><mo>/</mo><mi>T</mi><mo></mo><mo>–</mo><mo></mo><mn>1.12</mn></math>(8) The reactor was modelled as an isothermal plug-flow reactor; the molar flow rate <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>F</mi></math>of each species changes along the catalyst bed according to: <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><mo></mo><mi>d</mi><msub><mi>F</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo>/</mo><mi>d</mi><mi>W</mi><mo></mo><mo>=</mo><mo></mo><mo>−</mo><msub><mi>r</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo></mo><mi>d</mi><msub><mi>F</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo>/</mo><mi>d</mi><mi>W</mi><mo></mo><mo>=</mo><mo></mo><mo>+</mo><msub><mi>r</mi><mn>1</mn></msub></mtd></mtr></mtable></math> (9) <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><mi>d</mi><msub><mi>F</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub><mo>/</mo><mi>d</mi><mi>W</mi><mo></mo><mo>=</mo><mo></mo><mo>−</mo><msub><mi>r</mi><mn>2</mn></msub><mo></mo></mtd></mtr><mtr><mtd><mi>d</mi><msub><mi>F</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub><mo>/</mo><mi>d</mi><mi>W</mi><mo></mo><mo>=</mo><mo></mo><mo>+</mo><msub><mi>r</mi><mn>2</mn></msub></mtd></mtr></mtable></math>(10) where Fj is the molar flow rate of species j (kmol·h⁻¹), the subscript j identifying the species, and W is the mass of catalyst counted from the reactor inlet (kg). These ordinary differential equations were solved using the Runge–Kutta fourth-order method. The material balance was calculated for a feed of 5,000 bpd of light naphtha (C₅/C₆ fraction). Note to Tab. 6. The feed composition corresponds to the C₅/C₆ light-naphtha cut of the atmospheric distillation unit [13]; the product composition was calculated in the present work from the kinetic model, Eqs. (3–10), for a DIP recycle ratio RDIP = 0.65. The blending research octane numbers RONi listed in the last column are pure-component blending values taken from the literature [4, 13]; the value in the “Total” row is the blended RON of the stream calculated from Eq. (11). Table 6. Material balance of the isomerization unit (DIP recycle) for AI-92 <table> <thead> <tr> <td style="text-align: center;"> Component </td> <td style="text-align: center;"> Feed, wt% </td> <td style="text-align: center;"> Feed, t/d </td> <td style="text-align: center;"> Product, wt% </td> <td style="text-align: center;"> Product, t/d </td> <td style="text-align: center;"> RON </td> </tr> </thead> <tbody> <tr> <td> i–C5 </td> <td style="text-align: center;"> 22.0 </td> <td style="text-align: center;"> 147.4 </td> <td style="text-align: center;"> 46.2 </td> <td style="text-align: center;"> 309.5 </td> <td style="text-align: center;"> 92.3 </td> </tr> <tr> <td> n–C5 </td> <td style="text-align: center;"> 28.0 </td> <td style="text-align: center;"> 187.6 </td> <td style="text-align: center;"> 5.1 </td> <td style="text-align: center;"> 34.2 </td> <td style="text-align: center;"> 61.7 </td> </tr> <tr> <td> Cyclo-C5 </td> <td style="text-align: center;"> 3.0 </td> <td style="text-align: center;"> 20.1 </td> <td style="text-align: center;"> 3.4 </td> <td style="text-align: center;"> 22.8 </td> <td style="text-align: center;"> 101.3 </td> </tr> <tr> <td> n–C6 </td> <td style="text-align: center;"> 19.0 </td> <td style="text-align: center;"> 127.3 </td> <td style="text-align: center;"> 11.0 </td> <td style="text-align: center;"> 73.7 </td> <td style="text-align: center;"> 24.8 </td> </tr> <tr> <td> i–C6 (2-MP) </td> <td style="text-align: center;"> 12.0 </td> <td style="text-align: center;"> 80.4 </td> <td style="text-align: center;"> 15.3 </td> <td style="text-align: center;"> 102.5 </td> <td style="text-align: center;"> 73.4 </td> </tr> <tr> <td> 2,2-DMB </td> <td style="text-align: center;"> 3.0 </td> <td style="text-align: center;"> 20.1 </td> <td style="text-align: center;"> 9.4 </td> <td style="text-align: center;"> 63.0 </td> <td style="text-align: center;"> 91.8 </td> </tr> <tr> <td> 2,3-DMB </td> <td style="text-align: center;"> 4.0 </td> <td style="text-align: center;"> 26.8 </td> <td style="text-align: center;"> 5.3 </td> <td style="text-align: center;"> 35.5 </td> <td style="text-align: center;"> 104.3 </td> </tr> <tr> <td> MCP+Benzene </td> <td style="text-align: center;"> 5.0 </td> <td style="text-align: center;"> 33.5 </td> <td style="text-align: center;"> 2.4 </td> <td style="text-align: center;"> 16.1 </td> <td style="text-align: center;"> 91.0 </td> </tr> <tr> <td> C5+ </td> <td style="text-align: center;"> 4.0 </td> <td style="text-align: center;"> 26.8 </td> <td style="text-align: center;"> 1.9 </td> <td style="text-align: center;"> 12.7 </td> <td style="text-align: center;"> 42.0 </td> </tr> <tr> <td> Total </td> <td style="text-align: center;"> 100.0 </td> <td style="text-align: center;"> 670.0 </td> <td style="text-align: center;"> 100.0 </td> <td style="text-align: center;"> 670.0 </td> <td style="text-align: center;"> 87.4 </td> </tr> </tbody> </table> The research octane number of the isomerate was calculated by linear blending on a mass basis: <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>R</mi><mi>O</mi><msub><mi>N</mi><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>d</mi></mrow></msub><mo></mo><mo>=</mo><mtext>Σ</mtext><mo></mo><mfenced><mrow><msub><mi>x</mi><mi>i</mi></msub><mo></mo><mo>·</mo><mo></mo><mi>R</mi><mi>O</mi><msub><mi>N</mi><mi>i</mi></msub></mrow></mfenced><mo></mo><mo>=</mo><mo></mo><mn>87.4</mn></math>(11) where RONprod is the research octane number of the isomerate leaving the unit (-); xi is the mass fraction of component i in the isomerate (-), taken from the column “Product (wt%)” of Tab. 6; RONi is the blending octane number of the pure component i (-), last column of Tab. 6; and the summation is carried out over all components i of the product stream. The subscript prod (product) denotes the stream leaving the unit, as opposed to the subscript feed used for the stream entering it. With a DIP recycle ratio RDIP = 0.65 (i.e. 65% of the unconverted n-pentane separated in the DIP column is returned to the reactor), the total mass flow rate at the reactor inlet is GR = 670 + 0.65 · 187.6 = 791.9 t·d⁻¹ = 33.0 t·h⁻¹, where GR is the reactor feed including the recycle. The corresponding molar flow rates at the reactor inlet are <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>F</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub></math> = 178.7 kmol·h⁻¹ (309.5 t·d⁻¹ of n-pentane, fresh feed plus recycle, M = 72.15 g·mol⁻¹) and <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>F</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub></math> = 61.5 kmol·h⁻¹ (127.3 t·d⁻¹ of n-hexane, M = 86.18 g·mol⁻¹); these values are used in the heat balance below. The steady-state energy balance of the isomerization reactor was written in the following form, in which every term is an enthalpy or a heat flow expressed in kW: Qfeed + Qreaction = Qprod + Qloss (12) where the subscripts are abbreviations of English words (and not numerical indices), namely: Qfeed (feed) – the enthalpy flow of the stream entering the reactor, counted from the reference temperature Tref = 25°C, kW; Qreaction (reaction) – the heat released by the exothermic isomerization reactions inside the reactor, kW; it appears on the left-hand side of Eq. (12) because for an exothermic process the reaction is a heat source; Qprod (product) – the enthalpy flow of the stream leaving the reactor, kW; Qloss (losses) – the heat lost to the surroundings through the reactor wall and the insulation, kW, taken as 2% of Qfeed for a thermally insulated adiabatic reactor. The individual terms were calculated as Qfeed = GR ·cp· (Tin – Tref) and Qreaction =ΣFj· Xj·|ΔHr,j |, where GR is the mass flow rate of the reactor feed including the recycle (kg·h⁻¹), cp = 2.22 kJ·kg⁻¹·K⁻¹ is the mean specific heat capacity of the light naphtha, Tin = 150°C is the reactor inlet temperature, Fj is the molar flow rate of reactant j at the reactor inlet (mol·h⁻¹), Xj is its conversion (-), and ΔHr,j is the heat of the corresponding isomerization reaction (kJ·mol⁻¹). The heats of the isomerization reactions are exothermic and were taken from [15,16]: ΔHr,1 (n – C5H12 → i–C5H12 = –6.86 kJ·mol⁻¹ and ΔHr,2 (n – C6H14 → i–C6H14 = –8.37 kJ·mol⁻¹, the negative sign denoting heat release. With the conversions X1 = 0.75 (n-pentane, with recycle) and X1 = 0.60 (n-hexane): <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><msub><mi>Q</mi><mn>1</mn></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>F</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow></msub><mo>·</mo><mo></mo><msub><mi>X</mi><mn>1</mn></msub><mo>·</mo><mfenced open="|" close="|"><mrow><mtext>Δ</mtext><msub><mi>H</mi><mrow><mi>r</mi><mn>,1</mn></mrow></msub></mrow></mfenced><mo></mo><mo>=</mo><mo></mo><mn>178.7</mn></mtd></mtr><mtr><mtd><mo>·</mo><mn>10</mn><mo>³</mo><mo></mo><mo>×</mo><mo></mo><mn>0.75</mn><mo></mo><mo>×</mo><mo></mo><mn>6.86</mn><mo></mo><mo>=</mo><mo></mo><mn>9.19</mn><mo>·</mo><msup><mn>10</mn><mn>5</mn></msup><mo></mo><mi>k</mi><mi>J</mi></mtd></mtr><mtr><mtd><mo>·</mo><msup><mi>h</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo></mo><mo>=</mo><mo></mo><mn>255.4</mn><mo></mo><mi>k</mi><mi>W</mi></mtd></mtr></mtable></math>(13) <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><msub><mi>Q</mi><mn>2</mn></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>F</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>6</mn></msub><msub><mi>H</mi><mn>14</mn></msub></mrow></msub><mo>·</mo><mo></mo><msub><mi>X</mi><mn>2</mn></msub><mo>·</mo><mfenced open="|" close="|"><mrow><mtext>Δ</mtext><msub><mi>H</mi><mrow><mi>r</mi><mn>,2</mn></mrow></msub></mrow></mfenced></mtd></mtr><mtr><mtd><mo>=</mo><mo></mo><mn>178.7</mn><mo>·</mo><msup><mn>10</mn><mn>3</mn></msup><mo>×</mo><mo></mo><mn>0.75</mn><mo></mo><mo>×</mo><mo></mo><mn>6.86</mn><mo></mo><mo>=</mo><mo></mo><mn>9.19</mn></mtd></mtr><mtr><mtd><mo>·</mo><msup><mn>10</mn><mn>5</mn></msup><mo></mo><mi>k</mi><mi>J</mi><mo>·</mo><msup><mi>h</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo></mo><mo>=</mo><mo></mo><mn>255.4</mn><mo></mo><mi>k</mi><mi>W</mi></mtd></mtr></mtable></math>(14) <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>Q</mi><mrow><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>Q</mi><mn>1</mn></msub><mo></mo><mo>+</mo><mo></mo><msub><mi>Q</mi><mn>2</mn></msub><mo></mo><mo>=</mo><mo></mo><mn>255.4</mn><mo></mo><mo>+</mo><mo></mo><mn>85.8</mn><mo></mo><mo>=</mo><mo></mo><mn>341.2</mn><mo></mo><mi>k</mi><mi>W</mi></math>(15) <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><msub><mi>Q</mi><mrow><mi>f</mi><mi>e</mi><mi>e</mi><mi>d</mi></mrow></msub><mo>=</mo><mo></mo><msub><mi>G</mi><mi>R</mi></msub><mo>·</mo><msub><mi>c</mi><mi>p</mi></msub><mo>·</mo><mfenced><mrow><msub><mi>T</mi><mrow><mi>i</mi><mi>n</mi></mrow></msub><mo>−</mo><msub><mi>T</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi></mrow></msub></mrow></mfenced></mtd></mtr><mtr><mtd><mo>=</mo><mo></mo><mo>=</mo><mo></mo><mn>33.0</mn><mo>·</mo><msup><mn>10</mn><mn>3</mn></msup><mo>×</mo><mo></mo><mn>2.22</mn><mo></mo><mo>×</mo><mo></mo><mfenced><mrow><mn>150</mn><mo></mo><mo>–</mo><mo></mo><mn>25</mn></mrow></mfenced></mtd></mtr><mtr><mtd><mo>=</mo><mo></mo><mo>=</mo><mo></mo><mn>9.16</mn><mo>·</mo><msup><mn>10</mn><mn>6</mn></msup><mo></mo><mi>k</mi><mi>J</mi><mo>·</mo><msup><mi>h</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo></mo><mo>=</mo><mo></mo><mn>2544</mn><mo></mo><mi>k</mi><mi>W</mi></mtd></mtr></mtable></math>(16) <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><msub><mi>Q</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi></mrow></msub><mo></mo><mo>=</mo><mo></mo><mn>0.02</mn><mo></mo><mo>·</mo><mo></mo><msub><mi>Q</mi><mrow><mi>f</mi><mi>e</mi><mi>e</mi><mi>d</mi></mrow></msub><mo></mo><mo>=</mo><mo></mo><mn>50.9</mn><mo></mo><mi>k</mi><mi>W</mi><mo></mo><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>d</mi></mrow></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>Q</mi><mrow><mi>f</mi><mi>e</mi><mi>e</mi><mi>d</mi></mrow></msub><mo></mo><mo>+</mo><mo></mo><msub><mi>Q</mi><mrow><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub><mo></mo><mo>–</mo><mo></mo><msub><mi>Q</mi><mrow><mi>l</mi><mi>o</mi><mi>s</mi><mi>s</mi></mrow></msub><mo></mo><mo>=</mo><mo></mo><mn>2834</mn><mo></mo><mi>k</mi><mi>W</mi></mtd></mtr></mtable></math>(17) Because the isomerization reactions are only mildly exothermic, the heat released is small compared with the sensible heat of the feed (Qreaction / Qfeed = 0.13). The adiabatic temperature rise of the catalyst bed is therefore: <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>Δ</mi><msub><mi>T</mi><mrow><mi>a</mi><mi>d</mi></mrow></msub><mo></mo><mo>=</mo><mo></mo><msub><mi>Q</mi><mrow><mi>r</mi><mi>e</mi><mi>a</mi><mi>c</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub><mo></mo><mo>/</mo><mo></mo><mfenced><mrow><msub><mi>G</mi><mi>R</mi></msub><mo>·</mo><msub><mi>c</mi><mi>p</mi></msub></mrow></mfenced><mo></mo><mo>=</mo><mo></mo><mn>1.228</mn><mo>·</mo><msup><mn>10</mn><mn>6</mn></msup><mo></mo><mo>/</mo><mo></mo><mfenced><mrow><mn>33.0</mn><mo>·</mo><mn>10</mn><mo>³</mo><mo></mo><mo>×</mo><mo></mo><mn>2.22</mn></mrow></mfenced><mo></mo><mo>=</mo><mo></mo><mn>16.8</mn><mo></mo><mo>°</mo><mi>C</mi></math> (18) so that the outlet temperature of the reactor is about 167°C for an inlet temperature of 150°C. The reactor therefore operates under quasi-isothermal conditions: the temperature rise remains well inside the operating window of the chlorinated Pt/Al₂O₃ catalyst (140–180°C), no inter-stage cooling is required, and the shift of the thermodynamic equilibrium of Eqs. (7–8) caused by ΔTad reduces the equilibrium yield of isopentane by less than 2 rel.%. This result is consistent with industrial data for C₅/C₆ isomerization units, for which adiabatic temperature rises of 10–25°C are reported [13, 16]. The DIP column separates isopentane as the overhead product (mole fraction <math xmlns="http://www.w3.org/1998/Math/MathML"><msubsup><mi>x</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow><mi>D</mi></msubsup></math>≥ 0.95) from the C₆₊ bottoms (<math xmlns="http://www.w3.org/1998/Math/MathML"><msubsup><mi>x</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>5</mn></msub><msub><mi>H</mi><mn>12</mn></msub></mrow><mi>B</mi></msubsup></math>≤ 0.05). Here the superscripts D and B denote the distillate and the bottoms, respectively. At the average column temperature of 65°C the relative volatility of the key pair i-C₅H₁₂/n-C₅H₁₂ is α = 1.21. The minimum number of theoretical stages was obtained from the Fenske equation [17]: <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>N</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>=</mo><mi>l</mi><mi>n</mi><mfenced open="[" close="]"><mrow><mfenced><mrow><msup><mi>x</mi><mi>D</mi></msup><mo>/</mo><mfenced><mrow><mn>1</mn><mo></mo><mo>–</mo><mo></mo><msup><mi>x</mi><mi>D</mi></msup></mrow></mfenced></mrow></mfenced><mo></mo><mo>·</mo><mo></mo><mfenced><mrow><mfenced><mrow><mn>1</mn><mo></mo><mo>–</mo><mo></mo><msup><mi>x</mi><mi>B</mi></msup></mrow></mfenced><mo>/</mo><msup><mi>x</mi><mi>B</mi></msup></mrow></mfenced><mfenced open="]" close="["><mrow><mo></mo><mo>/</mo><mo></mo><mi>l</mi><mi>n</mi><mo></mo><mi>α</mi><mo></mo><mo>=</mo><mo></mo><mi>l</mi><mi>n</mi></mrow></mfenced><mfenced><mrow><mn>0.95</mn><mo>/</mo><mn>0.05</mn></mrow></mfenced><mo>²</mo></mrow></mfenced><mo></mo><mo>/</mo><mo></mo><mi>l</mi><mi>n</mi><mfenced><mn>1.21</mn></mfenced><mo></mo><mo>=</mo><mo></mo><mn>30.9</mn><mo></mo><mo>≈</mo><mo></mo><mn>31</mn></math>(19) where Nmin is the minimum number of theoretical stages at total reflux (-), xD and xB are the mole fractions of the light key (isopentane) in the distillate and in the bottoms (-), and α is the relative volatility of the key pair (-). The minimum reflux ratio calculated by the Underwood method [17] is Rr,min = 5.8; the operating reflux ratio was set at Rr = 1.3·Rr,min = 7.54, where the subscript r denotes reflux (R without a subscript is reserved for the universal gas constant, and RDIP for the recycle ratio of the unit). The actual number of theoretical stages obtained from the Gilliland correlation [17] is N = 52; with an overall tray efficiency of 75%, the number of actual trays is Nact = 70. The column diameter at 80% of flooding is Dc = 5.0 m. The column height was calculated as = Nact·ht+h0 = 70 × 0.6 + 5.0 = 47.0 m, where ht = 0.6 m is the tray spacing and h = 5.0 m is the allowance for the sump and the vapour disengagement space. Table 7. Design parameters of the isomerization reactor and DIP column <table> <thead> <tr> <td style="text-align: center;"> Parameter </td> <td style="text-align: center;"> Value </td> </tr> </thead> <tbody> <tr> <td> Isomerization Reactor </td> <td style="text-align: center;"></td> </tr> <tr> <td> Reactor type / catalyst </td> <td style="text-align: center;"> Fixed-bed / Pt/Al₂O₃ (Cl) </td> </tr> <tr> <td> Reactor diameter × length </td> <td style="text-align: center;"> 2.35 m × 7.05 m </td> </tr> <tr> <td> Catalyst mass </td> <td style="text-align: center;"> 18.45 t </td> </tr> <tr> <td> Temperature / pressure </td> <td style="text-align: center;"> 150°C / 25–30 bar </td> </tr> <tr> <td> LHSV / H₂:HC ratio </td> <td style="text-align: center;"> 2.0 h⁻¹ / 0.05–0.10 </td> </tr> <tr> <td> n–C₅ conversion (with recycle) </td> <td style="text-align: center;"> 75% </td> </tr> <tr> <td> Deisopentanizer Column </td> <td style="text-align: center;"></td> </tr> <tr> <td> Diameter / height </td> <td style="text-align: center;"> 5.0 m / 47.0 m </td> </tr> <tr> <td> Actual trays / tray spacing </td> <td style="text-align: center;"> 70 / 0.6 m </td> </tr> <tr> <td> Reflux ratio </td> <td style="text-align: center;"> июл.54 </td> </tr> <tr> <td> Operating pressure </td> <td style="text-align: center;"> 2.5 atm </td> </tr> <tr> <td> Reboiler / condenser duty </td> <td style="text-align: center;"> 6.8 MW / 6.5 MW </td> </tr> </tbody> </table> LHSV – liquid hourly space velocity, h⁻¹ <h3>Catalytic Reforming Unit Design for AI-95 and AI-98 Production</h3> Catalytic reforming converts the heavy naphtha fraction (C₇–C₁₀, boiling range 85–180°C) into high-octane reformate through four main reaction types on a Pt–Re/Al₂O₃ bifunctional catalyst: (a) dehydrogenation of naphthenes to aromatics; (b) isomerization of paraffins and naphthenes; (c) dehydrocyclization of paraffins to aromatics; and (d) hydrocracking (undesirable) [4, 11]. A semi-regenerative fixed-bed configuration with three reactors in series was selected. The main reactions and their kinetic expressions adopted for the four reaction families are given below. Dehydrogenation of naphthenes, represented by the conversion of cyclohexane to benzene <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>Δ</mi><msub><mi>H</mi><mrow><mi>r</mi><mn>,4</mn></mrow></msub><mo>=</mo><mn>206</mn><mo></mo><mi>k</mi><mi>J</mi><mo>·</mo><mi>m</mi><mi>o</mi><msup><mi>l</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup></math>. Its rate was described by a Langmuir–Hinshelwood expression, in which p denotes the partial pressure (atm) and the subscript identifies the species: <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>4</mn></msub><mo>=</mo><msub><mi>k</mi><mn>4</mn></msub><mo>·</mo><mfenced><mrow><msub><mi>p</mi><mi>N</mi></msub><mo>–</mo><msub><mi>p</mi><mi>A</mi></msub><mo>·</mo><msubsup><mi>p</mi><msub><mi>H</mi><mn>2</mn></msub><mn>3</mn></msubsup><mo>/</mo><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,4</mn></mrow></msub></mrow></mfenced><mo>/</mo><mfenced><mrow><mn>1</mn><mo>+</mo><msub><mi>K</mi><mi>A</mi></msub><mo>·</mo><msub><mi>p</mi><mi>A</mi></msub><mo>+</mo><msub><mi>K</mi><mi>N</mi></msub><mo>·</mo><msub><mi>p</mi><mi>N</mi></msub></mrow></mfenced></math>(20) where r4 is the rate of the dehydrogenation reaction (kmol·kgcat⁻¹·h⁻¹); k4 is its rate constant (kmol·kgcat⁻¹·h⁻¹·atm⁻¹); pN, pA and pH2 are the partial pressures of the naphthene (cyclohexane), of the aromatic (benzene) and of hydrogen, respectively (atm); Keq,4 is the equilibrium constant of that reaction (atm³); and KA and KN are the adsorption equilibrium constants of the aromatic and of the naphthene on the catalyst surface (atm⁻¹). The subscripts N and A therefore denote a class of hydrocarbons (naphthene, aromatic), whereas the subscript eq denotes a thermodynamic equilibrium constant. Dehydrocyclization of n-heptane to toluene is likewise strongly endothermic (); its rate was taken as first order in the paraffin: <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>5</mn></msub><mo>=</mo><msub><mi>k</mi><mn>5</mn></msub><mo>·</mo><msub><mi>p</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>7</mn></msub><msub><mi>H</mi><mn>16</mn></msub></mrow></msub></math>(21) Isomerization of n-paraffins, represented by the model reaction of n-heptane to its branched isomers (2-methylhexane, 3-methylhexane, 2,3-dimethylpentane and 2,4-dimethylpentane), lumped into the pseudo-component i–C7H16; it is only mildly exothermic (ΔHr,6 ≈ –8 kJ·mol⁻¹). <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>6</mn></msub><mo></mo><mo>=</mo><msub><mi>k</mi><mn>6</mn></msub><mo>·</mo><mfenced><mrow><msub><mi>p</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>7</mn></msub><msub><mi>H</mi><mn>16</mn></msub></mrow></msub><mo>–</mo><msub><mi>p</mi><mrow><mi>i</mi><mo>−</mo><msub><mi>C</mi><mn>7</mn></msub><msub><mi>H</mi><mn>16</mn></msub></mrow></msub><mo>/</mo><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,6</mn></mrow></msub></mrow></mfenced></math>(22) Hydrocracking of paraffins, to lighter products is an undesirable reaction that lowers the liquid yield ΔHr,7 = –45 kJ·mol⁻¹: <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>r</mi><mn>7</mn></msub><mo></mo><mo>=</mo><msub><mi>k</mi><mn>7</mn></msub><mo>·</mo><msub><mi>p</mi><mrow><mi>n</mi><mo>−</mo><msub><mi>C</mi><mn>7</mn></msub><msub><mi>H</mi><mn>16</mn></msub></mrow></msub><mo>·</mo><msub><mi>p</mi><msub><mi>H</mi><mn>2</mn></msub></msub></math>(23) where r5 – r7 are the rates of the dehydrocyclization, of the paraffin isomerization and of the hydrocracking reaction (kmol·kgcat⁻¹·h⁻¹); k5 – k7 are the corresponding rate constants; pn–C7 H16 and pш–C7 H16 are the partial pressures of n-heptane and of the lumped branched heptane isomers (atm); pH₂ is the partial pressure of hydrogen (atm); and Keq,6 is the equilibrium constant of the paraffin isomerization reaction (-). All partial pressures refer to the gas phase at the reactor conditions. Note to Tab. 8. Origin of the data. The activation energies Ea of dehydrogenation, dehydrocyclization, isomerization and hydrocracking over Pt–Re/Al₂O₃ were adopted from the classical kinetic scheme of Smith [18] and from the naphtha-reforming kinetic model of Ancheyta-Juárez and Villafuerte-Macías [19]; they are consistent with the industrial reformer model of Rahimpour, et al. [20]. The heats of reaction ΔHr are standard values calculated from the tabulated enthalpies of formation of the species involved and referred to 500 °C [4, 19]. The equilibrium constant of dehydrogenation Keq,4 at 500°C was calculated in the present work from Eq. (24); the equilibrium constants of dehydrocyclization and of isomerization were calculated from the corresponding van ’t Hoff relations with the thermodynamic data of [18, 19]. The rate constants k at 500°C were calculated in the present work from the Arrhenius equation, Eq. (6), and are expressed in kmol·kgcat⁻¹·h⁻¹·atm⁻ⁿ, the pressure order n being that of the corresponding rate expression, Eqs. (20–23). For hydrocracking no equilibrium constant is given (“–”) because the reaction is irreversible under the process conditions. As in Tab. 5, no parameter was fitted to the experimental data of the present study. Table 8. Kinetic parameters for catalytic reforming on Pt–Re/Al₂O₃ <table> <thead> <tr> <td style="text-align: center;"> Reaction </td> <td style="text-align: center;"> Ea, kJ·mol⁻¹ </td> <td style="text-align: center;"> ΔHr, kJ·mol⁻¹ </td> <td style="text-align: center;"> Keq at 500°C </td> <td style="text-align: center;"> k at 500°C </td> </tr> </thead> <tbody> <tr> <td> Dehydrogenation </td> <td style="text-align: center;"> 125.5 </td> <td style="text-align: center;"> 206 </td> <td style="text-align: center;"> 5.8×10⁴ </td> <td style="text-align: center;"> 8.4×10⁻² </td> </tr> <tr> <td> Dehydrocyclization </td> <td style="text-align: center;"> 159.0 </td> <td style="text-align: center;"> 264 </td> <td style="text-align: center;"> 3.2×10³ </td> <td style="text-align: center;"> 1.2×10⁻³ </td> </tr> <tr> <td> Isomerization </td> <td style="text-align: center;"> 108.0 </td> <td style="text-align: center;"> –8 </td> <td style="text-align: center;"> 1.45 </td> <td style="text-align: center;"> 5.1×10⁻² </td> </tr> <tr> <td> Hydrocracking </td> <td style="text-align: center;"> 220.0 </td> <td style="text-align: center;"> –45 </td> <td style="text-align: center;"> – </td> <td style="text-align: center;"> 3.5×10⁻⁴ </td> </tr> </tbody> </table> The equilibrium constant of the naphthene dehydrogenation reaction was calculated from the van ’t Hoff relation, with T in K: <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>ln</mi><mfenced><msub><mi>K</mi><mrow><mi>e</mi><mi>q</mi><mn>,4</mn></mrow></msub></mfenced><mo>=</mo><mo></mo><mo>–</mo><mn>24800</mn><mo>/</mo><mi>T</mi><mo></mo><mo>+</mo><mo></mo><mn>43.05</mn></math>(24) where the coefficient 24800 = ΔHr,4 /R (with ΔHr,4= 206 kJ·mol⁻¹) and the constant 43.05 = ΔSr,4 /R (with ΔSr,4 358 J·mol⁻¹·K⁻¹, the standard entropy change of that reaction, which is strongly positive because four moles of gas are formed from one) [18, 19]. At the design temperature of 500°C (773 K) Eq. (24) gives Keq,4 = 5.8·10⁴ atm³, so that the equilibrium the reaction is shifted almost completely towards the aromatic product, which is the thermodynamic basis of the high octane number of the reformate. The catalytic reforming unit processes 6,000 bpd of heavy naphtha (C₇–C₁₀ fraction). The feed and product compositions are presented in Tab. 9. Table 9. Material balance of the catalytic reforming unit <table> <thead> <tr> <td style="text-align: center;"> Component </td> <td style="text-align: center;"> Feed, wt% </td> <td style="text-align: center;"> Feed, t/d </td> <td style="text-align: center;"> Reformate, wt% </td> <td style="text-align: center;"> Reformate, t/d </td> <td style="text-align: center;"> RON </td> </tr> </thead> <tbody> <tr> <td> Paraffins (C₇–C₉) </td> <td style="text-align: center;"> 45.0 </td> <td style="text-align: center;"> 369.0 </td> <td style="text-align: center;"> 18.2 </td> <td style="text-align: center;"> 127.4 </td> <td style="text-align: center;"> 45–65 </td> </tr> <tr> <td> Naphthenes </td> <td style="text-align: center;"> 35.0 </td> <td style="text-align: center;"> 287.0 </td> <td style="text-align: center;"> 5.3 </td> <td style="text-align: center;"> 37.1 </td> <td style="text-align: center;"> 75–85 </td> </tr> <tr> <td> Aromatics (in feed) </td> <td style="text-align: center;"> 12.0 </td> <td style="text-align: center;"> 98.4 </td> <td style="text-align: center;"> 62.5 </td> <td style="text-align: center;"> 437.5 </td> <td style="text-align: center;"> 100–120 </td> </tr> <tr> <td> Light gases (C₁–C₄) </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 8.0 </td> <td style="text-align: center;"> 56.0 </td> <td style="text-align: center;"> – </td> </tr> <tr> <td> H₂ produced </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 3.0 </td> <td style="text-align: center;"> 21.0 </td> <td style="text-align: center;"> – </td> </tr> <tr> <td> LPG (C₃–C₄) </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 3.0 </td> <td style="text-align: center;"> 21.0 </td> <td style="text-align: center;"> – </td> </tr> <tr> <td> Other </td> <td style="text-align: center;"> 8.0 </td> <td style="text-align: center;"> 65.6 </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> 0.0 </td> <td style="text-align: center;"> – </td> </tr> <tr> <td> Total </td> <td style="text-align: center;"> 100.0 </td> <td style="text-align: center;"> 820.0 </td> <td style="text-align: center;"> 100.0 </td> <td style="text-align: center;"> 700.0* </td> <td style="text-align: center;"> 98.5 </td> </tr> </tbody> </table> LPG – Liquefied Petroleum Gas. Reformate yield = 85.4% (liquid volume); 14.6% converted to gases and hydrogen Note to Tab. 9. The feed composition is that of the C₇–C₉ heavy-naphtha cut [4]; the reformate composition was calculated in the present work from the kinetic model, Eqs. (20–24), at a reactor inlet temperature of 500°C. The ranges of blending octane numbers in the last column are pure-group values taken from the literature [4, 14]. The research octane number of the reformate was calculated by linear blending on a mass basis, in the same way as for the isomerate, Eq. (11): <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>R</mi><mi>O</mi><msub><mi>N</mi><mrow><mi>r</mi><mi>e</mi><mi>f</mi><mi>o</mi><mi>r</mi><mi>m</mi><mi>a</mi><mi>t</mi><mi>e</mi></mrow></msub><mo></mo><mo>=</mo><mtext>Σ</mtext><mo></mo><mfenced><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>·</mo><mi>R</mi><mi>O</mi><msub><mi>N</mi><mi>i</mi></msub></mrow></mfenced><mo></mo><mo>=</mo><mo></mo><mn>98.5</mn></math>(25) where xi is the mass fraction of the group of components i in the reformate (-) and RONi is the mean blending octane number of that group (-), both taken from Tab. 9; the light gases, the LPG and the hydrogen are excluded from the summation, since they do not remain in the liquid reformate. The high aromatic content (62.5 wt%) in the reformate provides the high RON necessary for blending into AI-95 and AI-98 grades. By adjusting the reformer severity (reactor temperature), the reformate RON can be tuned: at 490°C, RON ≈ 95 (suitable for AI-95 with blending); at 510°C, RON ≈ 100 (for AI-98 blending). Higher severity increases aromatic yield but reduces liquid yield and catalyst lifetime [4]. Catalytic reforming is strongly endothermic overall, being dominated by dehydrogenation and by dehydrocyclization. The steady-state energy balance of the three-reactor system, with all terms expressed in kW, was written as: Qfeed + Qheater = Qprod + Qreaction+ Qloss (26) where the subscripts, as in Eq. (12), are abbreviations of English words: Qfeed (feed) – the enthalpy flow of the heavy-naphtha feed entering the reactor section, kW; Qheater (heater) – the total heat duty supplied by the fired heaters, i.e. the charge heater plus the two inter-stage reheaters, kW; Qprod (product) – the enthalpy flow of the reactor effluent leaving the third reactor, kW; Qreaction (reaction) – the net heat consumed by the endothermic reforming reactions, kW; unlike Eq. (12), this term appears on the right-hand side, because for an endothermic process the reaction is a heat sink and not a heat source; Qloss (losses) – the heat lost from the reactors, the transfer lines and the heaters to the surroundings, kW, taken as 3% of Qheater. The net heat of reaction per mole of converted feed was estimated from the contributions of the four reaction families: <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>Δ</mi><msub><mi>H</mi><mrow><mi>o</mi><mi>v</mi><mi>e</mi><mi>r</mi><mi>a</mi><mi>l</mi><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>x</mi><mi>N</mi></msub><mo>·</mo><mi>Δ</mi><msub><mi>H</mi><mrow><mi>r</mi><mn>,4</mn></mrow></msub><mo>+</mo><msub><mi>x</mi><mrow><mi>D</mi><mi>H</mi><mi>C</mi></mrow></msub><mo>·</mo><mi>Δ</mi><msub><mi>H</mi><mrow><mi>R</mi><mn>,5</mn></mrow></msub><mo></mo><mo>+</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>s</mi><mi>o</mi></mrow></msub><mo>·</mo><mi>Δ</mi><msub><mi>H</mi><mrow><mi>r</mi><mn>,6</mn></mrow></msub><mo>+</mo><msub><mi>x</mi><mrow><mi>H</mi><mi>C</mi></mrow></msub><mo>·</mo><mi>Δ</mi><msub><mi>H</mi><mrow><mi>r</mi><mn>,7</mn></mrow></msub></math>(27) where ΔHoverall is the net heat of reaction per mole of feed converted (kJ·mol⁻¹); ΔHr,4, ΔHr,5, ΔHr,6 and ΔHr,7 are the heats of the dehydrogenation, dehydrocyclization, paraffin isomerization and hydrocracking reactions (kJ·mol⁻¹, Tab. 8); and xN, xDHC, xISO and xHC are the mole fractions of the feed converted by dehydrogenation of naphthenes (N), by dehydrocyclization (DHC), by isomerization (ISO) and by hydrocracking (HC), respectively (-). These four subscripts denote reaction families and not individual species. For the feed composition of Tab. 9 (45 wt% paraffins, 35 wt% naphthenes, 12 wt% aromatics), the converted fractions are xN = 0.35, xDHC = 0.20, xISO = 0.10 and xHC = 0.05, which gives: <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><mi>Δ</mi><msub><mi>H</mi><mrow><mi>o</mi><mi>v</mi><mi>e</mi><mi>r</mi><mi>a</mi><mi>l</mi><mi>l</mi></mrow></msub><mo>=</mo><mn>0.35</mn><mo></mo><mo>×</mo><mo></mo><mn>206</mn><mo></mo><mo>+</mo><mo></mo><mn>0.20</mn><mo></mo><mo>×</mo><mo></mo><mn>264</mn><mo></mo><mo>+</mo><mo></mo><mn>0.10</mn></mtd></mtr><mtr><mtd><mo></mo><mo>×</mo><mo></mo><mfenced><mrow><mo>–</mo><mn>8</mn></mrow></mfenced><mo></mo><mo>+</mo><mo></mo><mn>0.05</mn><mo></mo></mtd></mtr><mtr><mtd><mo>×</mo><mo></mo><mfenced><mrow><mo>–</mo><mn>45</mn></mrow></mfenced><mo></mo><mo>=</mo><mo></mo><mo>+</mo><mn>121.9</mn><mo></mo><mi>k</mi><mi>J</mi><mo>·</mo><mi>m</mi><mi>o</mi><msup><mi>l</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup></mtd></mtr></mtable></math>(28) The positive value confirms that the process is a net heat consumer, so that the heat withdrawn from the gas by the reactions has to be restored by inter-stage reheating between the reactors. The endothermic nature requires inter-stage reheating. With three reactors, the temperature profile: Reactor 1 inlet: 500°C → outlet: 430°C (ΔT = –70°C)(29) Reactor 2 inlet: 500°C → outlet: 455°C (ΔT = –45°C) (30) Reactor 3 inlet: 500°C → outlet: 475°C (ΔT = –25°C) (31) The decreasing temperature drop reflects the shift from fast dehydrogenation (Reactor 1) to slower dehydrocyclization and isomerization (Reactors 2–3). The total duty of the fired heaters required for the inter-stage reheating is: <math xmlns="http://www.w3.org/1998/Math/MathML"><mtable columnalign="left"><mtr><mtd><msub><mi>Q</mi><mrow><mi>h</mi><mi>e</mi><mi>a</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>G</mi><mi>F</mi></msub><mo>·</mo><msub><mi>c</mi><mi>p</mi></msub><mo></mo><mo>·</mo><mtext>Σ</mtext><mo></mo><mtext>Δ</mtext><msub><mi>T</mi><mi>i</mi></msub><mo></mo><mo>=</mo><mo></mo><mn>34.2</mn></mtd></mtr><mtr><mtd><mo>·</mo><mn>10</mn><mo>³</mo><mo></mo><mo>×</mo><mo></mo><mn>2.35</mn><mo></mo></mtd></mtr><mtr><mtd><mo>×</mo><mo></mo><mfenced><mrow><mn>70</mn><mo></mo><mo>+</mo><mo></mo><mn>45</mn><mo></mo><mo>+</mo><mo></mo><mn>25</mn></mrow></mfenced><mo></mo><mo>/</mo><mo></mo><mn>3600</mn><mo></mo><mo>=</mo><mo></mo><mn>3126</mn><mo></mo><mi>k</mi><mi>W</mi><mo></mo><mo>≈</mo><mo></mo><mn>3.1</mn><mo></mo><mi>M</mi><mi>W</mi></mtd></mtr></mtable></math> (32) where Qheater is the inter-stage reheating duty (kW); GF = 820 t·d⁻¹ = 34.2 t·h⁻¹ is the mass flow rate of the fresh heavy-naphtha feed (the subscript F denoting the fresh feed of the reforming unit); cp = 2.35 kJ·kg⁻¹·K⁻¹ is the mean specific heat capacity of the heavy naphtha; ΔTi is the temperature drop across reactor i, Eqs. (29–31); and the factor 3600 converts kJ·h⁻¹ into kW. The feed preheating duty from 25 to 500°C is accounted for separately in Tab. 10. Table 10. Heat balance of the catalytic reforming unit <table> <thead> <tr> <td style="text-align: center;"> Heat Component </td> <td style="text-align: center;"> Value (MW) </td> <td style="text-align: center;"> Notes </td> </tr> </thead> <tbody> <tr> <td> Feed preheat (25→500°C) </td> <td style="text-align: center;"> 10.6 </td> <td style="text-align: center;"> Fired heater </td> </tr> <tr> <td> R1→R2 interheat </td> <td style="text-align: center;"> 1.4 </td> <td style="text-align: center;"> Fired heater </td> </tr> <tr> <td> R2→R3 interheat </td> <td style="text-align: center;"> 0.9 </td> <td style="text-align: center;"> Fired heater </td> </tr> <tr> <td> R3→Stabilizer interheat </td> <td style="text-align: center;"> 0.8 </td> <td style="text-align: center;"> Fired heater </td> </tr> <tr> <td> Total fired heater duty </td> <td style="text-align: center;"> 13.7 </td> <td style="text-align: center;"></td> </tr> <tr> <td> Product cooling </td> <td style="text-align: center;"> 5.2 </td> <td style="text-align: center;"> Air + water coolers </td> </tr> </tbody> </table> Table 11. Design parameters of the catalytic reforming reactor system <table> <thead> <tr> <td style="text-align: center;"> Parameter </td> <td style="text-align: center;"> Value </td> </tr> </thead> <tbody> <tr> <td> Configuration </td> <td style="text-align: center;"> 3 reactors in series (semi-regenerative) </td> </tr> <tr> <td> Catalyst </td> <td style="text-align: center;"> Pt–Re/Al₂O₃ (bimetallic, 0.25% Pt, 0.25% Re) </td> </tr> <tr> <td> Total catalyst mass </td> <td style="text-align: center;"> 32.5 t (R1: 5.0, R2: 10.0, R3: 17.5) </td> </tr> <tr> <td> Catalyst distribution </td> <td style="text-align: center;"> 15% / 30% / 55% </td> </tr> <tr> <td> Reactor inlet temperature </td> <td style="text-align: center;"> 490–510°C </td> </tr> <tr> <td> Operating pressure </td> <td style="text-align: center;"> 15–25 bar </td> </tr> <tr> <td> WHSV </td> <td style="text-align: center;"> 1.5 h⁻¹ </td> </tr> <tr> <td> H₂/HC molar ratio </td> <td style="text-align: center;"> 5–7 </td> </tr> <tr> <td> Reactor diameters (R1/R2/R3) </td> <td style="text-align: center;"> 2.0 / 2.5 / 3.0 m </td> </tr> <tr> <td> Reactor lengths (R1/R2/R3) </td> <td style="text-align: center;"> 3.2 / 4.1 / 5.0 m </td> </tr> <tr> <td> Reformate RON (at 500°C) </td> <td style="text-align: center;"> 98.5 </td> </tr> <tr> <td> Reformate liquid yield </td> <td style="text-align: center;"> 85.4 vol% </td> </tr> <tr> <td> H₂ production </td> <td style="text-align: center;"> 21 t/d (by-product) </td> </tr> </tbody> </table> WHSV – weight hourly space velocity, h⁻¹ The final gasoline grades are produced by blending the isomerate and reformate with other refinery streams: Table 12. Gasoline blending scheme for AI-92, AI-95, and AI-98 production <table> <thead> <tr> <td style="text-align: center;"> Component </td> <td style="text-align: center;"> RON </td> <td style="text-align: center;"> AI-92, vol% </td> <td style="text-align: center;"> AI-95, vol% </td> <td style="text-align: center;"> AI-98, vol% </td> </tr> </thead> <tbody> <tr> <td> Isomerate </td> <td style="text-align: center;"> 87.4 </td> <td style="text-align: center;"> 55 </td> <td style="text-align: center;"> 20 </td> <td style="text-align: center;"> 10 </td> </tr> <tr> <td> Reformate (low severity) </td> <td style="text-align: center;"> 95.0 </td> <td style="text-align: center;"> 30 </td> <td style="text-align: center;"> 45 </td> <td style="text-align: center;"> 0 </td> </tr> <tr> <td> Reformate (high severity) </td> <td style="text-align: center;"> 100.0 </td> <td style="text-align: center;"> 0 </td> <td style="text-align: center;"> 15 </td> <td style="text-align: center;"> 65 </td> </tr> <tr> <td> FCC naphtha / alkylate </td> <td style="text-align: center;"> 92–96 </td> <td style="text-align: center;"> 15 </td> <td style="text-align: center;"> 20 </td> <td style="text-align: center;"> 25 </td> </tr> <tr> <td> Blended RON </td> <td style="text-align: center;"></td> <td style="text-align: center;"> 92.1 </td> <td style="text-align: center;"> 95.3 </td> <td style="text-align: center;"> 98.2 </td> </tr> </tbody> </table> The blending analysis confirms that the combination of isomerate (RON 87.4) and reformate (RON 95–100) provides the necessary octane quality for all three gasoline grades. The isomerate contributes the majority of the AI-92 blend (55 vol%), while reformate is the primary component for AI-95 (60 vol%) and AI-98 (65 vol%) production. <h3>Conclusion</h3> This study demonstrated the design of a combined isomerization-catalytic reforming system for producing the full spectrum of commercial gasoline grades, integrated with an AI-based quality prediction model. The main findings are: Physicochemical characterization of three commercial gasoline grades (AI-92, AI-95, AI-98) revealed strong correlations between density, viscosity, and octane number (|r| ≈ 0.99), confirming these parameters as reliable ML features. The Decision Tree model achieved exact RON predictions (MAE = 0.00, R² = 1.00) for all three grades. A light naphtha isomerization unit with DIP recycle was designed for AI-92 production. The mathematical model includes pseudo-first-order kinetics (Ea = 105.4 kJ·mol⁻¹ for n-C5 H12 , material and heat balances showing that the reactor operates under quasi-isothermal conditions (adiabatic temperature rise ΔTad = 16.8°C, i.e. an outlet temperature of about 167°C, within the operating window of the chlorinated Pt/Al₂O₃ catalyst), and a Fenske–Underwood–Gilliland design of the DIP column (70 actual trays, Dc = 5.0 m, Rr = 7.54). The isomerate RON = 87.4 serves as the primary blending component for AI-92 gasoline. A catalytic reforming unit with three semi-regenerative reactors was designed for AI-95 and AI-98 production. The kinetic model includes Langmuir–Hinshelwood dehydrogenation kinetics and associated material / heat balances. The reformate RON = 98.5 at 500°C with 85.4% liquid yield provides the high-octane base for AI-95 (blended RON 95.3) and AI-98 (blended RON 98.2). The gasoline blending analysis demonstrated that the full range of commercial grades can be produced by adjusting the proportions of isomerate, low-severity reformate, and high-severity reformate in the blending pool. Future work should expand the ML training dataset, validate the mathematical models through Aspen HYSYS simulation, and optimize the blending strategy through multi-objective optimization. <h3>ADDITIONAL INFORMATION</h3> Funding source. This study was not supported by any external sources of funding. Competing interests. The authors declare that they have no competing interests. Authors’ contribution. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work. The greatest contribution is distributed as follows: Aruzhan Abdigali – experimental characterization of gasoline samples, development of machine-learning models, and preparation of the first manuscript draft; Nuriya Ye. Aikenova – methodology of process modeling, selection and systematization of kinetic parameters, process-design decisions, and scientific supervision; Zhansaya Zhakhanshayeva – material and energy balance calculations, analysis of process-model results, preparation of tables and visual materials, and manuscript revision. <h3>ДОПОЛНИТЕЛЬНО</h3> Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования. Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с публикацией настоящей статьи. Вклад авторов. Все авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение исследования и подготовку статьи, прочли и одобрили финальную версию перед публикацией). Наибольший вклад распределён следующим образом: Абдигали А. – экспериментальная характеристика образцов бензина, разработка моделей машинного обучения и подготовка первого варианта рукописи; Айкенова Н.Е. – разработка методики математического моделирования, выбор и систематизация кинетических параметров, проектные решения и научное руководство; Жаханшаева Ж. – расчёт материальных и энергетических балансов, анализ результатов математического моделирования, подготовка таблиц и графических материалов, доработка рукописи. <a href="#_ftnref1" name="_ftn1">[1]</a> <a href="https://store.astm.org/d1298-12br17e01.html">ASTM D1298-12b(2017)e1 Standard Test Method for Density, Relative Density, or API Gravity of Crude Petroleum and Liquid Petroleum Products by Hydrometer Method</a> <a href="#_ftnref2" name="_ftn2">[2]</a> <a href="https://store.astm.org/d0445-24.html">ASTM D445-24 Standard Test Method for Kinematic Viscosity of Transparent and Opaque Liquids (and Calculation of Dynamic Viscosity)</a> <a href="#_ftnref3" name="_ftn3">[3]</a> <a href="https://analyzing-testing.netzsch.com/en/footer/legacy-products/sta-449-f3-jupiter">Simultaneous Thermogravimetry STA 449 F3 Jupiterю ® Fascinating Flexibility in Thermal Analysis</a>

iea.org [Internet]. IEA. Kazakhstan 2022: Energy Sector Review [cited 2026 May 05]. Available from: www.iea.org/reports/kazakhstan-2022.

akorda.kz [Internet]. President of the Republic of Kazakhstan. President Kassym-Jomart Tokayev’s State of the Nation Address to the People of Kazakhstan “Kazakhstan in the Era of Artificial Intelligence: Current Challenges and Solutions through Digital Transformation” [cited 2026 May 05]. Available from: node2.akorda.kz/en/president-kassym-jomart-tokayevs-state-of-the-nation-address-to-the-people-of-kazakhstan-kazakhstan-in-the-era-of-artificial-intelligence-current-challenges-and-solutions-through-digital-transformation-1083029.

Osman WS, Fadel AE, Salem SM, et al. Optimum Design of Naphtha Recycle Isomerization Unit with Modification by Adding De-Isopentanizer. Processes. 2023;11(12):3406. doi: 10.3390/pr11123406.

Jechura J. Gasoline Upgrading: Reforming, Isomerization, & Alkylation. Chapter 10&11. In: CBEN 409 – Petroleum Refining. Golden : Colorado School of Minesж 2019.

blog.axens.net [Internet]. Axens. How to Achieve Gasoline Pool Octane Requirement and Boost Light Cut Naphtha Contribution with an Environmental-friendly Process? [cited 2026 Apr 05]. Available from: https://blog.axens.net/how-to-achieve-gasoline-pool-octane-requirement.

kbr.com [Internet]. KBR. MAX-ISOMTM Maximize octane using innovative catalytic distillation principles [cited 2026 Apr 05]. Available from: www.kbr.com/sites/default/files/documents/2023-09/refining_maxisom_2022apr.pdf.

Lin X, Yin X, Shitao Y. Trustworthy Artificial Intelligence Framework for Predicting Gasoline Octane Loss Using Sparse Autoencoder and Stacking Ensemble Learning in Petrochemical Processes. Preprint;2025 [cited 2026 Apr 05]. doi: 10.20944/preprints202504.2593.v1.

Bounaceur R, Barthélemy N, Delort N, et al. A Multimodal Learning Model based on a QSPR approach for the estimation of RON, MON and CN, for any C, H, O hydrocarbons. Fuel. 2025;381(C):133438. doi: 10.1016/j.fuel.2024.133438.

Chen Z, vom Lehn F, Pitsch H, Cai L. Design of novel high-performance fuels with artificial intelligence: Case study for spark-ignition engine applications. Applications in Energy and Combustion Science. 2025;23:100341. doi: 10.1016/j.jaecs.2025.100341.

10.

Pedregosa F, Varoquaux G, Gramfort A, et al. Scikit-learn: Machine Learning in Python. Journal of Machine Learning Research. 2011;12(2825–2830).

11.

Fix E, Hodges JL. Discriminatory Analysis. Nonparametric Discrimination: Consistency Properties. International Statistical Review. 1989;57:238–247. doi: 10.2307/1403797.

12.

Breiman L, Friedman J, Olshen RA, Stone CJ. Classification and Regression Trees, 1st ed. New York: Chapman and Hall/CRC;1984.

13.

Mohamed MF, Shehata WM, Abdel Halim AA, Gad FK. Improving gasoline quality produced from MIDOR light naphtha isomerization unit. Egyptian Journal of Petroleum. 2017;26(1):111–124. doi: 10.1016/j.ejpe.2016.02.009.

14.

Benavides A, Zapata C, Benjumea P, et al. Predicting Octane Number of Petroleum-Derived Gasoline Fuels from MIR Spectra, GC-MS, and Routine Test Data. Processes. 2023;11(5):1437. doi: 10.3390/pr11051437.

15.

Volkova GG, Reshetnikov SI, Shkuratova LN, et al. n-Hexane skeletal isomerization over sulfated zirconia catalysts with different Lewis acidity. Chemical Engineering Journal. 2007;134(1–3):106–110. doi: 10.1016/j.cej.2007.03.079.

16.

Ivanchina ED, Chuzlov VA, Ivanchin NR, et al. Mathematical modeling of the process catalytic isomerization of light naphtha. Petroleum and Coal. 2019;61(2):413–417.

17.

Seader JD, Henley EJ, Keith Roper D. Separation Process Principles. 3rd ed. Wiley; 2010. 848 p.

18.

Smith RB. Kinetic analysis of naphtha reforming with platinum catalyst. Chemical Engineering Progress. 1959;55(6):76–80.

19.

Ancheyta-Juárez J, Villafuerte-Macías E. Kinetic Modeling of Naphtha Catalytic Reforming Reactions. Energy Fuels. 2000;14(5):1032–1037. doi: 10.1021/ef0000274.

20.

Rahimpour MR, Esmaili S, Ghalehghazi NB. A kinetic and deactivation model for industrial catalytic naphtha reforming. Iranian Journal of Science and Technology, Transaction B: Engineering. 2003;27(B2):279–290.