Next Article in Journal
A Review of Catalyst Integration in Hydrothermal Gasification
Previous Article in Journal
Comparative Analysis of Aeroshell 500 Oil Effects on Jet A and Diesel-Powered Aviation Microturbines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculation of the Rate Constants of Vacuum Residue Hydrogenation Reactions in the Presence of a Chrysotile/NiTi Nanocatalyst

by
Nazerke Balpanova
1,*,
Murzabek Baikenov
1,*,
Assanali Ainabayev
1,
Aikorkem Kyzkenova
2,
Gulzhan Baikenova
3 and
Almas Tusipkhan
1
1
Department of Chemical Technologies and Petrochemicals, Faculty of Chemistry, Karaganda Buketov University, Karaganda 100028, Kazakhstan
2
School of Pharmacy, Karaganda Medical University, Karaganda 100008, Kazakhstan
3
Department of Ecology and Assessment, Karaganda Economic University Kazpotrebsoyuz, Karaganda 1000000, Kazakhstan
*
Authors to whom correspondence should be addressed.
Fuels 2024, 5(3), 364-374; https://doi.org/10.3390/fuels5030021
Submission received: 2 July 2024 / Revised: 1 August 2024 / Accepted: 20 August 2024 / Published: 23 August 2024

Abstract

:
The paper presents the results of an investigation into the kinetics of catalytic hydrogenation of vacuum residue at temperatures of 380, 400 and 420 °C and different durations, ranging from 30 to 70 min, using a nanocatalyst containing the active metals nickel and titanium supported on chrysotile. It was found that the yield of oils from 30 to 50 wt.% and tars from 12 to 18 wt.% increased with increasing temperatures and reaction times. A slight increase in the proportion of solids in the range of 2.0 to 6.0 wt.% is explained by the activity of the nanocatalyst used. In the study of the kinetics of vacuum residue hydrogenation, using the nanocatalyst developed by the authors, we were able to achieve a low yield of solids with a short contact time as well as a high yield of low-molecular-weight compounds such as oils and tars. To determine the kinetic parameters (rate constants and activation energies), Simpson’s integral method and a random search engine optimization method were used. High values of rate constants are characteristic of reactions in the formation of oils k1, tars k2 and asphaltenes k3 in the temperature range of 380–420 °C. The high values of the rate constants k1, k2 and k3 in the catalytic hydrogenation of the vacuum residue indicate the high reaction rate and activity of the nanocatalyst used. With an increase in temperature from 380 to 420 °C, the rate constant of the formation of gas products from vacuum residue and the conversion of asphaltenes into oils significantly increase, which indicates the accumulation of low-molecular-weight compounds in oils. The activation energy for reactions leading to the formation of oils, tars, asphaltenes, gas and solid products was 75.7, 124.8, 40.7, 205.4 and 57.2 kJ/mol, respectively. These data indicate that the processes of vacuum residue hydrogenation with the formation of oils and asphaltenes require the lowest energy inputs. Reducing the process temperature to increase the selectivity of the vacuum residue hydrogenation process when using the prepared nanocatalyst is recommended. The formation of oils at the initial stage plays a key role in the technology of the heavy hydrocarbon feedstock (HHF) hydrogenation process. Perhaps the resulting oils can serve as an additional solvent for high-molecular-weight products such as asphaltenes, as evidenced by the low activation energy of the process.

1. Introduction

The vacuum residue is a heavy fraction that is formed after vacuum distillation of crude oil at temperatures of 500–520 °C. The share of vacuum residues can reach up to 30% of the total volume of oil refining, and this share is expected to grow [1]. Various methods are used to process heavy hydrocarbon feedstocks, including hydrocracking, thermal and catalytic cracking, and other technologies. Despite the many existing industrial methods of vacuum residue processing, the conversion depth of these processes remains insufficient. The main problem is the high concentration of high-molecular-weight components such as asphaltenes [2]. In oil and residual oil fractions, asphaltenes are in a solid state in the form of a dispersed phase, the precipitation of which is complicated by the presence of solvate shells of tars and polar heteroatomic compounds contained in the feedstock. In addition, as the viscosity of the produced oil increases, the content of asphaltenes and tars increases. These components contribute to the formation of a significant amount of solid coke-like products, deactivation of catalysts used in conversion processes, and reduction in distillate fraction yields [3]. Heteroatoms and metals mainly remain in asphaltenes and tars, which increases refining costs and leads to contamination and loss of activity of expensive catalysts [4]. In this regard, it is necessary to conduct an in-depth study of the chemical properties of asphaltenes, tars and oils in order to understand their behavior in the processes of processing heavy-oil residues. This will help to find a solution to eliminate the above-mentioned problem. There is a considerable amount of research devoted to analyzing the composition, properties and molecular structure of these components [5,6,7,8].
A large number of studies are devoted to the kinetics of the catalytic hydrogenation of vacuum residue, but the chemical aspects of these transformations remain insufficiently studied. This is due to the complexity of vacuum residue itself, the difficulty of comparing data obtained using different heterogeneous catalysts [9,10,11,12,13,14,15,16,17,18] and the need for a comprehensive analysis of the kinetics of hydrogenation of heavy hydrocarbon feedstock using the nanocatalysts. In the literature, there are few studies devoted to the kinetics of the process of vacuum residue hydrogenation with the use of nanocatalysts. In [19], the authors demonstrated the effectiveness of the nanocatalyst in the hydrocracking of vacuum residue. However, during hydrocracking, nanoparticles tend to agglomerate, forming spherical and agglomerated structures, which reduces the activity of the nanocatalyst. To reduce the agglomeration of nanoparticles in the works [20,21], chrysotile, which consists of nanotubes, was used as a catalyst carrier.
The study of the kinetics of vacuum residue hydrogenation is complicated by the following two parallel processes: the formation of low-molecular compounds (oils) due to degradation and condensation of asphaltenes and tars into solid products, which leads to the formation of coke. The solution of this problem will make it possible to acquire new knowledge about carrying out the process of catalytic hydrogenation of vacuum residue in industrial conditions. Therefore, the problem of obtaining new experimental data on the thermal destruction of asphaltenes, tars and oils in the presence of the nanocatalysts remains unsolved. Calculation of kinetic parameters (rate constants and activation energy) of catalytic hydrogenation of vacuum residue will significantly help to control the processes of formation of low-molecular compounds and reduce the formation of solid products (coke).
Differential equations are used to model the chemical kinetics of the hydrogenation of coal and coal tar as well as heavy-oil residues. In previous studies [22,23,24,25,26], when studying the kinetics of hydrogenation of heavy and solid hydrocarbon raw materials, rigid types of differential equations were used to calculate kinetic parameters. It is known that the chemical kinetics of complex reactions is characterized by the presence of rapidly and slowly changing variables. Due to the presence of reaction stages with different rates, the solution of direct kinetic problems is complicated by the rigidity of the differential equation systems describing the mechanism of these reactions [27,28,29,30]. In recent years, Runge–Kutta methods, previously considered reliable, began to cede their dominant position among algorithms for solving ordinary differential equations in favor of modified methods capable of coping with rigid problems [31,32]. To solve the problems of chemical kinetics, the method of random search optimization [33] and the Simpson’s integral method [34] were first used. Previously, these methods were used separately to solve systems of differential equations and were not used together to solve direct kinetic problems.
The purpose of this work is to establish the kinetic patterns of the vacuum residue hydrogenation in the presence of a chrysotile-based nanocatalyst modified with metals such as nickel and titanium.

2. Materials and Methods

The study was conducted on vacuum residue produced by the Pavlodar Petrochemical Plant (PPCHP). The main physical and chemical characteristics of this vacuum residue are presented in Table 1. The kinematic viscosity of the initial material was determined using an automatic viscometer SYD-265B-1, and the density was measured in accordance with GOST 1481-81.
Preparation of the nanocatalyst. As a catalyst and carrier, mountain flax was used (also known as chrysotile), which is produced by Kostanay Minerals JSC, a mining company specializing in the extraction of chrysotile asbestos and the production of chrysotile fiber. The chemical formula of chrysotile is Mg3Si2O5(OH)4. Chrysotile fibers are formed from twisted layers of MgO and SiO2 [35].
In this work, a chrysotile-based nanocatalyst with additives of active metals, such as nickel and titanium (Chrysotile/NiTi), was used. The nanocatalyst (Chrysotile/NiTi) was prepared by a wet-mixing method using drain-free technology, as follows: chrysotile was previously treated with a 20% hydrochloric acid solution to remove magnesium, potassium and other alkaline earth elements; then, a solution of nickel nitrate and titanyl sulfate was mixed with chrysotile at 100 ± 5 °C for 60 min; then, the precipitate was separated, dried at 100 ± 5 °C and calcined at 550 °C.
Physical and chemical characteristics of the prepared nanocatalyst.
X-ray spectral fluorescence, gravimetric analysis of the initial chrysotile. Using X-ray spectral fluorescence analysis on spectrometer CPM-25 (X-ray MULTICHANNEL SPECTROMETER, Laborant Industrial Group, Russia, ), the chemical composition of the initial chrysotile was determined (Table 2).
From Table 2, it follows that the main components of chrysotile include silicon, magnesium and iron oxides.
Transmission electron microscopy. The particle size of the deposited nickel and titanium on the surface of chrysotile nanotubes was determined using a Joel Jem-1400Plus transmission electron microscope (JEOL, Freising, Germany) at an accelerating voltage of 120 kV. Figure 1 shows that the average diameter of the nanotubes in the nanocatalyst sample is 40 nm.
In the prepared nanocatalyst, TEM revealed the presence of finely dispersed particles of nickel and titanium oxide (Figure 2a) and, based on statistical processing, a curve of the distribution of nickel particles by size was plotted (Figure 2b). Figure 2b shows that nickel particles of 6–9 nm and 29–35 nm are evenly distributed on the surface of chrysotile nanotubes, while the distribution peak falls at 12 nm, which is 28% of the total amount of adsorbed nickel and titanium metals.
Atomic emission spectral analysis of the nanocatalyst. Using an inductively coupled plasma Profile Plus atomic emission spectrometer (Teledyne Instruments Leeman Labs Inc., Hudson, NH, USA) under standard conditions (temperature 21 °C and atmospheric pressure 727 mmHg), the amount of supported metals on chrysotile was measured.
According to atomic emission spectral analysis, the content of nickel and titanium in the nanocatalyst was 5%.
Investigation of the textural characteristics of the nanocatalyst. The porous properties of the initial chrysotile and the nanocatalyst based on it, containing nickel and titanium nanolayers, were studied. For this purpose, a low-temperature nitrogen adsorption method at 77 K on an automated Digisorb-2600 Micromeritics unit was used. Prior to measurements, the sample was vacuum-treated at 350 °C for 5 h. The classical method of analysing the isotherms of polymolecular adsorption (BET method) using the instrument software was also applied. The BET analysis gave a surface area for the initial chrysotile and nanocatalyst of 39.4 and 54 m2/g, respectively.
Thermo-programmable desorption. To confirm the acidic properties of metal-supported chrysotile, a thermo-programmable ammonia desorption assay was performed. The ammonia desorption rate curve shows a pronounced peak at 200 °C, which indicates the presence of acid centers at a concentration of 267 μmol/g.
Vacuum residue hydrogenation. The vacuum residue hydrogenation experiments were carried out in a 0.3 L high-pressure reactor with an internal stirrer (HXCHEM, Shanghaiha, China). The starting materials (vacuum residue and nanocatalyst) were premixed and placed in the reactor. The mass fraction of the nanocatalyst was 1.0% of the initial fraction. The reactor was closed and filled with hydrogen at an overpressure of 4.0 MPa. The temperature in the reactor was then raised at a rate of 10 °C/min to a predetermined value and held there for the specified time. The reactor was then cooled to room temperature. After 24 h, the reactor was opened, the solid product was separated from the liquid and the solid hydrogenation product was thoroughly washed with benzene.
Determination of component composition. To isolate the asphaltenes, the sample was diluted with a 40-fold volume of hexane and allowed to stand for a day, and then the precipitate was filtered off. The resulting precipitate was placed in a paper cartridge and washed with hexane in a Soxhlet apparatus to remove oils and tars, after which the asphaltenes were removed from the cartridge with chloroform. Deasphalted samples were applied to activated coarse silica gel, hydrocarbon components (oils) were sequentially extracted using hexane in a Soxhlet apparatus and the tars were isolated with a 1:1 mixture of benzene and ethanol [36].
Method for determining kinetic parameters of vacuum residue hydrogenation. To determine the kinetic parameters (rate constant and activation energy), the Simpson’s integral method and the random search method were used [33,34]. The solution of the differential equations system is carried out using the Simpson’s method, which allows us to calculate rate constants with a calculation error of up to 12%, and further calculation using the random search optimization method allows us to reduce the discrepancy between the obtained calculated data and the experimental values of the concentrations of the obtained products with an accuracy of 5%. As part of the study, a formalized mechanism of mutual transformations of asphaltenes, tars, oils, gas, and vacuum residue product in the hydrogenation process was developed (Figure 3).
When developing the scheme, it was assumed that, during the hydrogenation of vacuum residue, direct and parallel reactions of decay and condensation of its main components occur. Based on this mechanism, the following kinetic model of the vacuum residue hydrogenation process was developed:
d C 1 d τ = ( k 1 + k 1 + k 3 + k 4 + k 5 ) C 1 ,
d C 2 d τ = k 1 C 1 + k 6 C 3 + k 7 C 4 ,
d C 3 d τ = k 2 C 1 k 6 C 3 k 8 C 3 ,
d C 4 d τ = k 3 C 1 k 7 C 4 ,
d C 5 d τ = k 4 C 1 + k 8 C 3 ,
d C 6 d τ = k 5 C 1 .
where, C1, C2, C3, C4, C5 and C6—mass fractions of vacuum residue, oil, tar, asphaltene, gas and solid product (coke); τ—process duration, min; k1, k2, k3, k4 and k5-transformation of vacuum residue into oils, tars, asphaltenes, gas and solid products; k6 and k7—reactions from vacuum residue and asphaltene transformations into oils; k8—reaction of vacuum residue transformation into gas products.

3. Results and Discussion

To determine the kinetic parameters of the catalytic hydrogenation of vacuum residue, experiments were carried out at temperatures of 380, 400 and 420 °C with durations from 30 to 70 min (Table 3). During hydrogenation, the destruction of tars and asphaltenes is observed, which increases the yield of oils, gases and solids.
With increasing temperatures and durations, an increase in the yield of oils from 30 to 50 wt.% and tars from 12 to 18 wt.%, respectively, was observed. The slight increase in solid products under the given conditions from 2.0 to 6.0 wt.% seems to be due to the activity of the prepared nanocatalyst (Chrysotile/NiTi).
Comparison of experimental data on the hydrocracking kinetics of vacuum residues without a catalyst and in the presence of a catalyst obtained in [2], where vacuum residues from different refineries were used without a catalyst, the yield of solid product varied from 19 to 38 wt.%, and, when the nanocatalyst was used, the yield of solid product decreased by 3–3.5 times, which shows the high activity of the nanocatalyst.
However, analysis of the literature has shown that the low yield of solids is associated with the duration of the process (72 h) [9]. By investigating the kinetics of vacuum residue hydrogenation in the presence of the nanocatalyst we prepared, a low solid yield was achieved with a short contact time, as well as a high yield of low-molecular-weight compounds such as oils and tars (Table 3).
Similarly, an increase in gas yield from 9 to 14% and asphaltenes from 3 to 8.4% is observed at the given temperatures and contact time duration. Increasing the duration and temperature of catalytic hydrogenation of vacuum residue contributes to a slight increase in the yield of new asphalt-containing substances, which is associated with the condensation reaction of asphaltenes and leads to the formation of solid products. The increase in the yield of gaseous products is due to catalytic destruction of vacuum residue.
Based on the above experimental data, a scheme of the vacuum residue hydrogenation in the presence of the nanocatalyst is proposed (Figure 2). Based on the kinetic scheme, a kinetic model of vacuum residue hydrogenation was compiled and rate constants were calculated. Table 4 shows the calculated values of the rate constants and the activation energies of the vacuum residue hydrogenation reaction.
High values of rate constants (Table 4) are observed for the reaction of formation of oils (k1 = 0.53·10−2), resins (k2 = 0.38·10−2) and asphaltenes (k3 = 0.48·10−2) at a temperature of 380 °C. Apparently, the high values of the rate constants k1, k2, k3 during the catalytic hydrogenation of vacuum residue indicate a high rate of the hydrogenation reaction and the activity of the prepared nanocatalyst. With an increase in temperature to 420 °C, the values of the rates constants during the formation of gas products from vacuum residue (k4 = 0.60·10−2) and the reaction of converting asphaltenes to oils (k7 = 1.10·10−2) sharply increase, which indicates the accumulation of low-molecular-weight compounds in oils, which is consistent with the data presented in [36]. As the temperature increases, there is a slight increase in the rate constants from 0.10·10−2 to 33·10−2 min−1 for the conversion of tars to oils (k6) and a sharp increase in the rate constants for the conversion of tars to gas products from 0.03·10−2 to 0.95·10−2 min−1 (k8), which is apparently due to the formation of a low-molecular-weight part of the products. Low values of rate constants from 0.19·10−2 to 0.36·10−2 min−1 are characteristic of the reaction of solid formation (k5), which once again emphasizes the high activity and selectivity of the nanocatalyst. Furthermore, it should be noted that the nanocatalyst used has an advantage over known nanocatalyzers which undergo agglomeration. Figure 1 and Figure 2 show that the prepared nanocatalyst does not undergo agglomeration, since metal nanoparticles are adsorbed on the surface and inside nanotubes with dimensions of 5–40 nm, which complicates the process of agglomeration of nickel and titanium nanoparticles.
To determine the activation energy of the main chemical transformations in the process of the vacuum residue hydrogenation, the dependences of reaction constants in the Arrhenius coordinates were plotted (Figure 3).
Figure 4 shows a graphical display of the dependence of the logarithm of the constant rate of conversion of the vacuum residue to oils with a reverse in temperature. Using the method of least squares in the angle of inclination on the line in the Arrhenius coordinates, the values of the activation energy of the process of converting the vacuum residue into components were calculated. The activation energies for the oil, tar, asphaltene, gas and solid formation reactions are presented in Table 4 and are 75.7, 124.8, 40.7, 205.4 and 57.2 kJ/mol, respectively. From the data obtained, it follows that the processes of the vacuum residue hydrogenation with the formation of oils and asphaltenes require the lowest activation energy. High values of activation energy are characteristic of reactions during the formation of tar (124.8 kJ/mol) and gas products from vacuum residue (205.4 kJ/mol), as well as the formation of oils from tars and asphaltenes (110.5 and 254.5 kJ/mol) and the conversion of tars into gas products (310.8 kJ/mol). It follows from these data that the processes of destruction of the organic mass of the vacuum residue with the formation of oils and asphaltenes (75.7 and 40.7 kJ/mol) are the least difficult energetically. Therefore, to increase the selectivity of the catalytic vacuum residue hydrogenation process, the lower process temperature should be limited. It is believed that the preferential formation of oils at the initial stage is important for the technology of the HHF hydrogenation process, since the formation of oils is an additional solvent for high-molecular-weight products such as asphaltenes, which confirms the low activation energy.

4. Conclusions

  • It was found that finely dispersed particles of nickel oxide and titanium were present in the prepared nanocatalyst, and a particle size distribution curve of these metals was also plotted. The photomicrograph shows that metal particles of 6–9 nm and 29–35 nm are uniformly distributed on the surface of chrysotile nanotubes, with the maximum distribution being observed at 12 nm, which is 28% of the total amount of adsorbed metals. The difference of this nanocatalyst from known ones lies in the fact that the deposited nickel and titanium oxides do not agglomerate, which is explained by the adsorption of nanoparticles both on the surface and inside the nanotubes.
  • A study of the kinetics of vacuum residue hydrogenation showed that the composition of high-molecular-weight compounds such as tars and asphaltenes significantly affects the material balance of the process. It was found that, with increasing temperatures and reaction times, an increase in the yield of oils from 30 to 50 wt.% and tars from 12 to 18% by weight was observed, respectively. The slight increase in the proportion of solids from 2.0 to 6.0 wt.% under these conditions is probably due to the activity of the nanocatalyst used (Chrysotile/NiTi). In the early stages of vacuum residue hydrogenation, the preferential formation of oils is of key importance for heavy-oil feedstock processing technology, since these oils serve as additional hydrogen donors for the formation of low-molecular-weight compounds and prevent precipitation of high-molecular-weight compounds such as asphaltenes and tars on the reactor walls.
  • Using the Simpson’s integral method and the random search optimization method, the rate constants of the catalytic hydrogenation of vacuum residue were calculated. High values of constants of speed are observed for the reaction in the formation of oils (k1 = 0.53·10−2), tars (k2 = 0.38·10−2) and asphaltenes (k3 = 0.48·10−2) at a temperature of 380 °C. These high values of constants indicate the rapid progress of the hydrogenation reaction and the high activity of the nanocatalyst used. When the temperature rises to 420 °C, there is an increase in the rate constants for the formation of gas products from vacuum residue (k4 = 0.60·10−2) and for the reaction of converting asphaltenes to oils (k7 = 1.10·10−2), which indicates the accumulation of low-molecular-weight compounds in oils. The low values of speed constants for the reaction in the formation of solid products (k5), varying from 0.19·10−2 up to 0.36·10−2 min−1, demonstrate the high activity and selectivity of the nanocatalyst. Calculations have shown that the value of the rate constants for the formation of oils, tars and asphaltenes depends not only on the initial concentration of components of heavy petroleum residues but also on the activity and selectivity of the nanocatalyst used.
  • Based on the obtained data on the kinetics of vacuum residue hydrogenation, activation energy values were calculated for the conversion of the vacuum residue into various components. The activation energy for the formation reactions of oils, tars, asphaltenes, gases and solids are 75.7, 124.8, 40.7, 205.4 and 57.2 kJ/mol, respectively. The highest values of activation energy are observed for the reactions during the formation of tars (124.8 kJ/mol) and gas products from vacuum residue (205.4 kJ/mol), as well as for the transformation of tars into gas products (310.8 kJ/mol) and the formation of oils from resins and asphaltenes (110.5 and 254.5 kJ/mol). The processes of destruction of the organic mass of vacuum residue with the formation of oils and asphaltenes (75.7 and 40.7 kJ/mol) are the least energy-consuming. Therefore, lower heating temperatures are recommended to increase the selectivity of the catalytic hydrogenation of the vacuum residue.
  • The obtained kinetic data can be used in the design of the reactor and the selection of active and selective catalysts for the destructive hydrogenation process of HHF; in addition, it will enable researchers to study the mechanism of conversion of HHF to low-molecular-weight compounds.

Author Contributions

Conceptualization, N.B. and M.B.; methodology, N.B. and M.B.; software, A.A. and A.K.; validation, G.B., A.A. and A.K.; formal analysis, N.B., M.B. and A.T.; investigation, N.B., M.B. and A.T.; resources, N.B., M.B., G.B., A.A. and A.K.; data curation, N.B. and M.B.; writing—original draft preparation, N.B. and M.B.; writing—review and editing, N.B. and M.B.; visualization, G.B., A.A. and A.K.; supervision, A.T., G.B., A.A. and A.K.; project administration, N.B. and M.B.; funding acquisition, N.B. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the financial support of the Ministry of Science and Higher Education of the Republic of Kazakhstan, project No. AP13268918 “Nanocatalytic System for the Hydrotreatment of Heavy Hydrocarbons”, agreement No. 143/ZhG-1-22-24 of 21 June 2022.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kim, S.H.; Kim, K.D.; Lee, Y.K. Effects of dispersed MoS2 catalysts and reaction conditions on slurry phase hydrocracking of vacuum residue. J. Catal. 2017, 347, 127–137. [Google Scholar] [CrossRef]
  2. Goncharov, A.V.; Krivtsov, E.B. Calculation of the rate constants of thermal cracking and condensation reactions of high-sulfur tar asphaltenes. Solid Fuel Chem. 2022, 56, 108–115. [Google Scholar] [CrossRef]
  3. Zhao, Y.; Xu, C.; Zhao, S.; Shi, Q. Pattern Recognition Technology Application in Intelligent Processing of Heavy Oil. Energy Fuels 2012, 26, 7251–7256. [Google Scholar] [CrossRef]
  4. Strausz, O.P.; Safarik, I.; Lown, E.M.; Morales-Izquierdo, A. A Critique of Asphaltene Fluorescence Decay and Depolarization-Based Claims about Molecular Weight and Molecular Architecture. Energy Fuels 2008, 22, 1156–1166. [Google Scholar] [CrossRef]
  5. Antipenko, V.R.; Grinko, A.A. Grinko Parameters of macrostructure of insoluble products obtained by thermolysis of resins and asphaltenes of the Usinskaya oil. Bull. Tomsk. Polytech. Univ. Geo Assets Eng. 2021, 332, 123–131. [Google Scholar] [CrossRef]
  6. Korneev, D.S.; Pevneva, G.S.; Golovko, A.K. Thermal transformations of asphaltenes at a temperature of 120 °C. J. Sib. Fed. Univ. Chem. 2019, 12, 101–117. [Google Scholar] [CrossRef]
  7. Krivtsov, E.B.; Goncharov, A.V. Effect of Styrene Additives on the Cracking Kinetics of Components of High-Sulfur Vacuum Residue. Pet. Chem. 2020, 60, 358–364. [Google Scholar] [CrossRef]
  8. Goncharov, A.V.; Krivtsov, E.B. Calculation of the kinetic parameters of the formation and destruction of the thiophene derivatives in high-sulphur vacuum residue cracking. Khim. Interes. Ustoich. Razvit. 2021, 29, 699–717. [Google Scholar] [CrossRef]
  9. Coronel-García, M.A.; de la Torre, A.R.; Domínguez-Esquivel, J.M.; Melo-Banda, J.A.; Martínez-Salazar, A.L. Heavy oil hydrocracking kinetics with nano-nickel dispersed in PEG300 as slurry phase catalyst using batch reactor. Fuel 2021, 283, 118930. [Google Scholar] [CrossRef]
  10. Alonso-Ramírez, G.; Cuevas-García, R.; Sánchez-Minero, F.; Ramirez, J.; Moreno-Montiel, M.; Silva-Oliver, G.; Ancheyta, J.; Carbajal-Vielman, R. Catalytic hydrocracking of a Mexican heavy oil on a MoS2/Al2O3 catalyst: II. Study of the transformation of isolated aromatics fraction obtained from SARA analysis. Fuel 2021, 288, 119541. [Google Scholar] [CrossRef]
  11. Félix, G.; Ancheyta, J. Using separate kinetic models to predict liquid, gas, and coke yields in heavy oil hydrocracking. Ind. Eng. Chem. Res. 2019, 58, 7973–7979. [Google Scholar] [CrossRef]
  12. Mishra, P.; Yadav, A. Modelling of a vacuum residue hydrocracking in an industrial slurry phase reactor. Can. J. Chem. Eng. 2023, 101, 7275–7292. [Google Scholar] [CrossRef]
  13. de Oliveira, L.P.; Verstraete, J.J.; Kolb, M. Simulating vacuum residue hydroconversion by means of Monte-Carlo techniques. Catal. Today 2014, 220, 208–220. [Google Scholar] [CrossRef]
  14. Kohli, K.; Prajapati, R.; Maity, S.K.; Sau, M.; Sharma, B.K. Accelerated pre-coking of NiMo/γ-Al2O3 catalyst: Effect on the hydroprocessing activity of vacuum residue. Fuel 2019, 235, 437–447. [Google Scholar] [CrossRef]
  15. Manek, E.; Haydary, J. Hydrocracking of vacuum residue with solid and dispersed phase catalyst: Modeling of sediment formation and hydrodesulfurization. Fuel Process. Technol. 2017, 159, 320–327. [Google Scholar] [CrossRef]
  16. Puron, H.; Arcelus-Arrillaga, P.; Chin, K.; Pinilla, J.; Fidalgo, B.; Millan, M. Kinetic analysis of vacuum residue hydrocracking in early reaction stages. Fuel 2014, 117, 408–414. [Google Scholar] [CrossRef]
  17. Nguyen, N.T.; Kang, K.H.; Pham, H.H.; Go, K.S.; Van Pham, D.; Seo, P.W.; Nho, N.S.; Lee, C.W.; Park, S. Catalytic hydrocracking of vacuum residue in a semi-batch reactor: Effect of catalyst concentration on asphaltene conversion and product distribution. J. Ind. Eng. Chem. 2021, 102, 112–121. [Google Scholar] [CrossRef]
  18. Armas, J.D.; Garcia-Vila, A.; Ortega, L.C.; Scott, C.E.; Maini, B.; Pereira-Almao, P. In-situ upgrading in a dolomite porous medium: A kinetic model comparison and nanocatalyst deposition study. J. Petrol. Sci. Eng. 2021, 205, 108799. [Google Scholar] [CrossRef]
  19. Zhang, H.; Duan, T.; Zhu, W.; Yao, W.T. Chrysotile-Based Nanowires Decorated with Monodispersed Ag Nanoparticles as a Highly Active and Reusable Hydrogenation Catalyst. J. Phys. Chem. 2015, 119, 21465–21472. [Google Scholar] [CrossRef]
  20. Swift, H. Catalytic Properties of Synthetic Layered Silicates. In Advanced Materials in Catalysis; Academic Press: Cambridge, MA, USA, 1977; pp. 209–233. [Google Scholar] [CrossRef]
  21. Balpanova, N.Z.; Gyul’maliev, A.M.; Pankin, Y.N.; Aitbekova, D.E.; Ma, F.; Su, K.; Kaikenov, D.A.; Baikenova, G.G.; Borsynbaev, A.S.; Baikenov, M.I. Kinetics of Hydrogenation of Heavy and Solid Hydrocarbon Raw Materials. Solid Fuel Chem. 2019, 53, 319–323. [Google Scholar] [CrossRef]
  22. Gulmaliev, A.M.; Golovin, G.S.; Gladun, T.G. Theoretical Foundations of Coal Chemistry; Publishing House of Moscow State Mining University: Moscow, Russia, 2003; p. 556. [Google Scholar]
  23. Dai, F.; Gong, M.; Li, C.; Li, Z.; Zhang, S. New kinetic model of coal tar hydrogenation process via carbon number component approach. Appl. Energy 2015, 137, 265–272. [Google Scholar] [CrossRef]
  24. Luo, X.; Guo, Q.; Zhang, D.; Zhou, H.; Yang, Q. Simulation, exergy analysis and optimization of a shale oil hydrogenation process for clean fuels production. Appl. Therm. Eng. 2018, 140, 102–111. [Google Scholar] [CrossRef]
  25. Dai, F.; Zhang, Y.; Xia, E.; Zhang, Z.; Zhang, Z.; Li, C. Modeling of fixed bed reactor for coal tar hydrogenation via the kinetic lumping approach. Carbon Resour. Convers. 2018, 1, 279–283. [Google Scholar] [CrossRef]
  26. Yunusov, A.A.; Gubaidullin, I.M.; Bainazarova, N.M. Optimization of numerical methods for solving rigid problems of chemical kinetics. In Proceedings of the Scientific Service on the Internet: Search for New Solutions: Third International Supercomputer Conference, Moscow, Russia, 17–22 September 2012; p. 278. [Google Scholar]
  27. Niemeyer, K.E.; Curtis, N.J.; Sung, C.J. pyJac: Analytical Jacobian generator for chemical kinetics. Comput. Phys. Commun. 2017, 215, 188–203. [Google Scholar] [CrossRef]
  28. Sun, J.; Li, D.; Yao, R.; Sun, Z.; Li, X.; Li, W. Modeling the hydrotreatment of full range medium temperature coal tar by using a lumping kinetic approach. React. Kinet. Mech. Catal. 2015, 114, 451–471. [Google Scholar] [CrossRef]
  29. Niemeyer, K.E.; Sung, C.-J. Accelerating moderately stiff chemical kinetics in reactive-flow simulations using GPUs. J. Comput. Phys. 2014, 256, 854–871. [Google Scholar] [CrossRef]
  30. Shi, Y.; Green, W.H.; Wong, H.-W.; Oluwole, O.O. Redesigning combustion modeling algorithms for the Graphics Processing Unit (GPU): Chemical kinetic rate evaluation and ordinary differential equation integration. Combust. Flame 2011, 158, 836–847. [Google Scholar] [CrossRef]
  31. Le, H.P.; Cambier, J.-L.; Cole, L.K. GPU-based flow simulation with detailed chemical kinetics. Comput. Phys. Commun. 2013, 184, 596–606. [Google Scholar] [CrossRef]
  32. Reizlin, V.I. Numerical Optimization Methods; Publishing House Tomsk Polytechnic University: Tomsk, Russia, 2011; p. 105. [Google Scholar]
  33. Verzhbitsky, V.M. Fundamentals of Numerical Methods; Publishing House High School: Moscow, Russia, 2002; p. 840. [Google Scholar]
  34. Balpanova, N.; Baikenov, M. Thermal Degradation Kinetics of Vacuum Residues in the Presence of Chrysotile Supported Ni-Ti Catalyst. Catalysts 2023, 13, 1361. [Google Scholar] [CrossRef]
  35. Dmitriev, D.E.; Golovko, A.K. Transformations of resins and asphaltenes during the thermal treatment of heavy oils. Pet. Chem. 2010, 50, 106–113. [Google Scholar] [CrossRef]
  36. Sviridenko, N.N. Cracking of heavy oil in supercritical water in the presence of iron oxide nanopowder: Asphaltene transformations and process kinetics. Pet. Chem. 2023, 63, 585–593. [Google Scholar] [CrossRef]
Figure 1. Micrograph of the nanocatalyst sample.
Figure 1. Micrograph of the nanocatalyst sample.
Fuels 05 00021 g001
Figure 2. Photomicrograph of nickel particles on the surface of the nanocatalyst obtained through TEM (a) and the particle size distribution curve of nickel particles on the surface (b).
Figure 2. Photomicrograph of nickel particles on the surface of the nanocatalyst obtained through TEM (a) and the particle size distribution curve of nickel particles on the surface (b).
Fuels 05 00021 g002
Figure 3. Scheme of generalized kinetic model of vacuum residue hydrogenation.
Figure 3. Scheme of generalized kinetic model of vacuum residue hydrogenation.
Fuels 05 00021 g003
Figure 4. Dependence of the logarithm of rate constant of vacuum residue to oil conversion on the inverse temperature.
Figure 4. Dependence of the logarithm of rate constant of vacuum residue to oil conversion on the inverse temperature.
Fuels 05 00021 g004
Table 1. Physicochemical characteristics of vacuum residue (PPCHP).
Table 1. Physicochemical characteristics of vacuum residue (PPCHP).
IndicatorContent
Viscosity at 70 °C, cSt262
Density, kg/m31.000
Atomic ratio H/C1.95
Ash content, %1.64
Hydrogen, %11.8
Carbon, %72.5
Nitrogen, %0.9
Sulphur, %2.6
Oxygen, %2.3
Asphaltenes, %9.8
Tars, %12.5
Oils, %78.2
Table 2. X-ray spectral fluorescence analysis of the initial chrysotile.
Table 2. X-ray spectral fluorescence analysis of the initial chrysotile.
Defined ComponentsContent, %
SiO266.595
TiO20.027
Al2O3<0.95
Fe2O33.641
CaO0.270
MgO25.555
MnO0.033
P2O50.021
K2O<0.1
Na2O0.557
LDC3.40
LDC—losses during calcination: H2O, CO2, carbon, hydrogen, partly SO3, trace elements.
Table 3. Component composition of catalytic vacuum residue hydrogenation products.
Table 3. Component composition of catalytic vacuum residue hydrogenation products.
T, Kτ, minContent, wt.%
VR,
C1
Oil,
C2
Tar,
C3
Asphaltene,
C4
Gas,
C5
Coke,
C6
6533044.730.111.72.78.82.0
4037.534.512.33.39.03.4
5030.040.312.53.59.24.5
6026.041.614.73.79.44.6
7023.243.015.54.09.54.8
6733037.837.113.54.25.22.2
4028.641.114.24.67.44.1
5022.245.114.94.88.54.5
6017.746.016.25.29.85.1
7012.948.217.15.510.85.5
6933024.142.115.15.78.64.4
4010.247.116.07.913.35.5
506.249.517.18.013.55.7
604.650.217.48.313.75.8
703.650.517.68.414.05.9
Table 4. Kinetic parameters of the vacuum residue hydrogenation process.
Table 4. Kinetic parameters of the vacuum residue hydrogenation process.
T, °CRate Constants, × 10−2 min−1
k1k2k3k4k5k6k7k8
3800.530.380.480.060.190.100.070.03
4000.840.670.570.290.270.220.190.13
4201.201.470.740.600.360.331.100.95
E, kJ/mol75.7124.840.7205.457.2110.5254.5310.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Balpanova, N.; Baikenov, M.; Ainabayev, A.; Kyzkenova, A.; Baikenova, G.; Tusipkhan, A. Calculation of the Rate Constants of Vacuum Residue Hydrogenation Reactions in the Presence of a Chrysotile/NiTi Nanocatalyst. Fuels 2024, 5, 364-374. https://doi.org/10.3390/fuels5030021

AMA Style

Balpanova N, Baikenov M, Ainabayev A, Kyzkenova A, Baikenova G, Tusipkhan A. Calculation of the Rate Constants of Vacuum Residue Hydrogenation Reactions in the Presence of a Chrysotile/NiTi Nanocatalyst. Fuels. 2024; 5(3):364-374. https://doi.org/10.3390/fuels5030021

Chicago/Turabian Style

Balpanova, Nazerke, Murzabek Baikenov, Assanali Ainabayev, Aikorkem Kyzkenova, Gulzhan Baikenova, and Almas Tusipkhan. 2024. "Calculation of the Rate Constants of Vacuum Residue Hydrogenation Reactions in the Presence of a Chrysotile/NiTi Nanocatalyst" Fuels 5, no. 3: 364-374. https://doi.org/10.3390/fuels5030021

APA Style

Balpanova, N., Baikenov, M., Ainabayev, A., Kyzkenova, A., Baikenova, G., & Tusipkhan, A. (2024). Calculation of the Rate Constants of Vacuum Residue Hydrogenation Reactions in the Presence of a Chrysotile/NiTi Nanocatalyst. Fuels, 5(3), 364-374. https://doi.org/10.3390/fuels5030021

Article Metrics

Back to TopTop