Thermal Behavior of Heavy Oil Catalytic Pyrolysis and Aquathermolysis

: There is still considerable controversy surrounding the mechanisms, thermodynamics, and kinetics of heavy oil aquathermolysis and pyrolysis processes. The present paper aims to widen our knowledge about the effect of iron tallates on pyrolysis and aquathermolysis of Cuban heavy oil. The obtained SARA (S: saturates, A: aromatics, R: resins, A: asphaltenes) analysis has shown a signiﬁcant increase in light hydrocarbon content during aquathermolysis. Moreover, the elemental analysis has indicated an increase in C and H content by almost 4% and 6%, respectively, with a signiﬁcant decrease in S and O content by up to 23% in the presence of iron tallates. These results have been further conﬁrmed by infrared spectrometry. The obtained IR data indicated that asphaltene and resin compounds transform into light hydrocarbons after aquathermolysis. On another hand, the activation energy of heavy oil pyrolysis decreased in the presence of the utilized catalyst; meanwhile, the reaction rate increased, especially in the temperature range of 200–480 ◦ C, which may validate a signiﬁcant effect of the used catalyst in real conditions. Moreover, the obtained thermodynamic data showed a decrease in the enthalpy and entropy of activation of oil pyrolysis in the presence of iron tallates. Our results are encouraging in terms of energy consumption, optimization, and process control and should be validated by a larger sample size. pyrolysis process and the fit DSC curves for this process based on a model corresponding to n-order autocatalysis reaction with m-power by-product (Cnm) in the presence and the absence of iron tallates.


Introduction
Unconventional hydrocarbons are widely considered to be a potential alternative resource for satisfying energy demand associated with population growth and the shortage of conventional resources of energy such as light oil [1][2][3][4]. However, these hydrocarbons have been neglected due to the lack of efficient methods [5]. Heavy and extra-heavy oils are believed to be the most available source of energy on the planet [6][7][8]. These oils are characterized by a high viscosity which complicates the process of their extraction and processing. Nowadays, different types of methods, known as enhanced oil recovery methods, have been developed for the extraction of heavy and extra-heavy oil [9][10][11][12][13][14]. Among them, thermal methods still generate tremendous interest from experts due to their economic and environmental prospects [15].
It is common knowledge that thermal enhanced oil recovery methods are widely used to decrease the viscosity of heavy oil in situ. Moreover, these methods are able to increase significantly oil mobility toward producing wells. Perhaps, the main known methods widely used for heavy oil extraction nowadays are the processes of cyclic steam stimulation (CSS) [16] and steam-assisted gravity drainage (SAGD) [17,18]. Nevertheless, one of the major drawbacks of these processes is the low recovery factor and the high viscosity of produced oil which requires further processing and complicates transportation [19].
Within the next few years, catalyst application is likely to become an important component in the processes of steam injection [20] because catalytic agents can decrease the activation energies of associated reactions and increase the rate at which they occur, as shown by several previous works [21,22]. To improve the process of steam injection, many experts [23][24][25] now argue that a deeper understanding of pyrolysis, aquathermolysis, and cracking reactions may be the main key, especially when it comes to the application of catalytic agents for this process.
Recently, a growing body of literature has investigated the effect of different catalysts on aquathermolysis and pyrolysis reactions [20,26,27]. In their works, Vakhin et al. [28][29][30] have studied the effect of different oil-soluble metal-based catalysts and found that Ni, Co, Fe, and Cu increase the saturate and aromatic contents in the oil after the steam effect.
It has been suggested that oil-soluble catalyst precursors transform during the steam injection processes and turn into nanoparticles able to break down the organosulfur compounds in asphaltenes by cleaving C-S bonds [19]. Unfortunately, most previous studies on oil aquathermolysis and pyrolysis reaction mechanisms neglected the kinetic and thermodynamic aspects of the process and did not focus on the energy's variation in the presence and absence of catalysts. The kinetic and thermodynamic parameters such as activation energy, enthalpy, and entropy are the main parameters that allow estimating the efficiency of the adopted catalyst.
This paper aims to validate the results obtained at larger laboratory scales regarding catalytic aquathermolysis and pyrolysis reactions in the presence of iron tallates from an analytical, non-isothermal kinetic, and thermodynamic point of view. In this context, we investigated the effect of iron tallates on oil composition during aquathermolysis and pyrolysis reactions. In addition, we calculated activation energies, enthalpy, and entropy of activation of the process of heavy oil pyrolysis by using a set of thermal analyses based on non-isothermal kinetics.

Materials
In order to investigate the samples' composition and the kinetic and thermodynamic behavior of aquathermolysis and pyrolysis reactions, we have opted for a Cuban extraheavy oil provided by JSC Zarubezhneft, Moscow, Russia from Boca de Jaruco's extra-heavy oil field. The physical properties of the used oil are described in Table 1. Organic solvents purer than 99.5% were obtained from Component Reactiv and utilized without additional purification. Inorganic salts were obtained from Sigma-Aldrich, Moscow, Russia. The physical characteristics of the used oil are presented in Table 1.

Catalyst Synthesis
The steps regarding iron tallate synthesis proceed very much in the same way as indicated in [31] with some modifications. Distilled tall oil (DTO) was chosen as a ligand- The obtained iron tallates (61.8 g) have been used further for catalytic pyrolysis experiments.

Thermal Processing of the Extra-Heavy Oil and Its SARA and Elemental Analysis
To model the steam injection process of the extra-heavy oil thermal processing, we used a batch reactor equipped with a heated cylindrical thick-walled vessel made of steel, a thermocouple, an electric drive, and a manometer. The reactor temperature control was provided by a heating jacket and a cooling coil. The composition content of the initial oil has been obtained by SARA analysis. The obtained results are presented in Tables 2 and 3. It should be noted that the obtained oil from the thermal processing in the batch reactor is referred to as the control experiment throughout the text.

Infrared Spectrometry Essays
So that we could study the structural group composition of the investigated products, we calculated the spectral coefficients. These coefficients define the ratio of the optical density values at the maxima of the corresponding absorption bands where C 1 = D 1600 /D 720 (aromaticity), C 2 = D 1710 /D 1465 (oxidation), C 3 = D 1380 /D 1465 (branching), C 4 = (D 720 + D 1380 )/D 1600 (aliphaticity), and C 5 = D 1030 /D 1465 (saturation) using a baseline in the 2000-400 cm −1 spectrum. The first coefficient represents aromaticity which results from stretching vibrations of C=C arom bonds of aromatic rings. Next, the second coefficient (C 2 ) represents the oxidation which is associated with the presence of carbonyl groups (CO). Moreover, the third coefficient (C 3 ) is related to the structure of paraffin chains which can be estimated by their branching from CH 3 /CH 2 . The fourth factor is aliphaticity (CH 3 + CH 2 )/C = S arom , and it shows the proportion of methyl and methylene groups to aromatic groups; meanwhile, the fifth coefficient C 5 shows the proportion of sulfoxide groups.

Sample Preparation
So that we would be able to study the effect of iron tallates on extra-heavy oil pyrolysis reactions, we prepared a sample of pure extra-heavy oil and 2 wt.% of iron tallates in initial oil for thermal analysis.

Thermal Analysis
At a temperature range of 30-600 • C, differential scanning calorimetry (DSC) analysis was performed by using an STA 449 F1 Jupiter (Netzsch, Selb, Germany) thermoanalyzer at linear heating rates of 5, 10, 15, and 20 • C·min −1 , under 50 mL·min −1 flow of argon gas. Moreover, the mass of the studied samples was~10 mg in the Al crucible. Data processing was performed using Proteus Analysis 5.2.1 and NETZSCH Kinetics Neo 2.1.2.2 program package.
Broadly speaking, steam-based processes are believed to include different reactions types such as pyrolysis [32,33] (thermolysis, thermal cracking), aquathermolysis, and gasification [34]. Moreover, these reactions are not widely understood in detail due to their complex nature and the heterogeneous mediums in which they occur. Therefore, determining the kinetic behavior of such reactions requires modern and sophisticated physical and chemical methods which may allow a detailed description of the processes associated with steam-based enhanced oil recovery methods. One of the most promising methods in the field of kinetics is the non-isothermal kinetic approach, which is consistent with the characteristics of such processes as steam injection. It is based on studying the aquathermolysis and pyrolysis reactions at different heating rates, unlike the conventional approach, which is based on studying different reactions at a constant temperature (isothermal principle), due to significant heat variations associated with these processes which lead to uncertain results.

Kinetic Theory
The process of extra-heavy oil pyrolysis is considered quite complicated for kinetic description, and it is often described as a function of conversion degree α, the preexponential factor A, and the apparent activation energy E α in addition to the reaction model f (α) as follows:

Isoconversional Kinetic Analysis
Thermal analysis in this study aims to determine the kinetic triplets of the processes occurring during steam-based methods of enhancing oil recovery. Therefore, we chose Kissinger-Akahira-Sunose (KAS) [35] and Friedman [36] methods as model-free analyses to obtain the energy of activation E α and preexponential factor A of these processes.
The Kissinger-Akahir-Sunose (KAS) method is presented by the following equation: The Friedman method is described as follows: where the index i identifies an individual heating rate and T α,i is the temperature at which the extent of conversion α is reached under ith heating rate. Then for any given α, the value of E α is estimated from the slope of a plot of ln dα which depends mainly on the obtained values of activation energy, pre-exponential factor, reaction order, and rate. These models were optimized utilizing KineticsNeo software (Netzsch). Kinetic models have been performed with the lowest number of steps to obtain an acceptable fit and consistent approach with the pyrolysis process of Cuban extra-heavy oil. Table 4 presents the classical solid reaction kinetic description. Table 4. Models' methods for calculating kinetic parameters.

Model Equation
Reaction of nth order (Fn) The reaction of nth order with m-power autocatalysis by-product (Cnm)

Chemical Composition and Elemental Analysis
The results obtained from SARA analysis and elemental analysis of the initial heavy oil and oil after the batch reactor processing at higher temperature (250 • C) and pressure (90 atm) in the absence (control experiment) and the presence of iron tallates are presented in Figure 1 and Table 5 Next, for assuming the kinetic model and the reaction type, we proceeded to apply the model-based analysis which depends mainly on the obtained values of activation energy, pre-exponential factor, reaction order, and rate. These models were optimized utilizing KineticsNeo software (Netzsch). Kinetic models have been performed with the lowest number of steps to obtain an acceptable fit and consistent approach with the pyrolysis process of Cuban extra-heavy oil. Table 4 presents the classical solid reaction kinetic description.

Chemical Composition and Elemental Analysis
The results obtained from SARA analysis and elemental analysis of the initial heavy oil and oil after the batch reactor processing at higher temperature (250 °C) and pressure (90 atm) in the absence (control experiment) and the presence of iron tallates are presented in Figure 1 and Table 5, respectively.   As seen in Figure 1, the quantity of saturated hydrocarbons in the presence of iron tallates is relatively higher than the amount found in the control experiment and even in the initial sample. This is probably due to the destruction of aromatic compounds which are in lower quantities in the catalytic process compared to the control and initial samples thermal treating. Moreover, Figure 1 shows a significant opposite relationship between asphaltene and resin content changes. In other words, the amount of asphaltenes in the catalytic experiments is less than the amount of asphaltenes in control and initial samples. Therefore, the increase in the amount of resins could be explained by the conversion of a part of asphaltene compounds into resin fractions.
On another hand, the elemental analysis of catalytic and non-catalytic aquathermolysis products (Table 5) indicates higher carbon and hydrogen amounts (73.60% and 9.95%) in the presence of iron tallates, unlike the non-catalytic process which showed low values for these elements (69.34% and 9.30% for carbon and hydrogen, respectively). The good effect of the used catalyst is expressed by the amounts of sulfur and oxygen as well, which decreased from 8.40% and 12.60%, respectively, in its presence to 6.43% and 9.65% in the non-catalytic process.
The same pattern has been established for the SARA elemental analysis which has indicated an increase in C amount in the saturates and aromatics in the presence of the catalyst compared to its absence as an indicator of C-C destruction reactions from aromatics into saturates. This could be confirmed by the amount of H in the saturates and aromatics, which decreased from 11.21% in the absence of the catalyst to 10.79% in its presence, and this is an indicator of hydrogen consumption during aromatic ring destruction in the presence of iron tallates.
The decrease in sulfur and oxygen amount in the presence of iron tallates is a result of catalyst precursor transformation into iron sulfide and iron oxygen molecules via thermal degradation at higher temperatures. The obtained iron sulfides and oxides probably present nanoparticle form and, in turn, play the role of the catalytic agent of the occurring reactions. There was a significant correlation between the elemental analysis of oil samples and their SARA fraction elemental data. Table 5 shows the significant increase in sulfur and oxygen content of the resin and asphaltene samples obtained in the presence of iron tallates, confirming our aforementioned hypotheses about the conversion of asphaltenes into resins which are found to enrich the amounts of the saturate and aromatic fractions; as a result, the heteroatoms would be present in greater quantities compared to carbon and hydrogen atoms. So, to comprehend the composition change behavior and to support our findings, we have applied infrared spectrometry. The obtained results of the infrared spectrometry are presented in Table 6. According to the obtained results, we can observe the rise in aromaticity C 1 from 0.93 in the absence of the catalyst to 1.13 in its presence, which evidences the process of asphaltene destruction into aromatic-ring-containing resins. Moreover, it should be noted that in both thermal experiments, the aromaticity factor C 1 decreases, which also evidences the conversion of aromatic compounds into saturated hydrocarbons, as has been previously demonstrated by SARA and elemental analysis. Moreover, the increase in the oxidation factor C 2 indicates the occurrence of slight intermolecular oxidation processes leading to some oxygenated compounds contributing to improving the quality of the produced oil by forming some active catalytic agents such as iron oxide nanoparticles. On another hand, the obtained data show a significant increase in the branched aliphatic compounds as indicated by the branching factor C 3 values, which increase from 0.56 in the absence of catalyst for both the initial and the control experiment to 0.68 after the catalytic effect of iron tallates. This is in good agreement with elemental and SARA analysis, which suggests that the increase in saturated hydrocarbons is the result of the thermal cracking of the aromatic compounds within the oil composition. The aliphatic factor C 4 did not show a significant difference between the control experiment and the experiment conducted in the presence of iron tallates. On contrary, it decreased relative to the initial oil sample, which is expected because of the evaporation of light gases and hydrocarbons during the thermal processing. Finally, the saturation factor value is 3 times higher for the heavy oil catalytic process compared to those factors related to the initial and control experiment. This is a good indicator of the formation of sulfoxide groups in the presence of iron tallate catalysts, which enhance the catalytic processes of heavy oil pyrolysis, aquathermolysis, and cracking via the formation of iron sulfide nanoparticles at higher temperatures.

Kinetic Calculations
We opted for a small differential scanning calorimetry (DSC) sample size in the present study because this type of analysis is well known for better estimation of kinetic parameters [37,38] of extra-heavy oil pyrolysis process in the presence and absence of different catalysts.
So that we could compare the effect generated by the obtained catalyst on the process of extra-heavy oil pyrolysis, we utilized DSC. The DSC curves of the non-catalytic and catalytic pyrolysis are presented in Figure 2.
The obtained curves (endothermic peaks from 350 to 500 • C) indicate the pyrolysis zone. They show an adequate behavior for both catalytic and non-catalytic pyrolysis processes, which allowed us to implement kinetic calculations.  The obtained curves (endothermic peaks from 350 to 500 °C) indicate the pyrolysis zone. They show an adequate behavior for both catalytic and non-catalytic pyrolysis processes, which allowed us to implement kinetic calculations. Table 7 contains the different values of activation energy and pre-exponential factor obtained during the pyrolysis processes in the presence and in the absence of iron tallates.  Table 7 contains the different values of activation energy and pre-exponential factor obtained during the pyrolysis processes in the presence and in the absence of iron tallates. In fact, the activation energies and the preexponential factors associated with Friedman, KAS, and ASTM E2890 values have been selected for a conversion degree of 0.5. These values are illustrated in Figure 3. Table 7. Kinetic parameters of Cuban extra-heavy oil pyrolysis process obtained from the isoconversional methods for a conversion degree of 0.5.

ASTM E2890
ASTM E2890 ASTM E2890 In fact, the activation energies and the preexponential factors associated with Friedman, KAS, and ASTM E2890 values have been selected for a conversion degree of 0.5. These values are illustrated in Figure 3. Table 7. Kinetic parameters of Cuban extra-heavy oil pyrolysis process obtained from the isoconversional methods for a conversion degree of 0.5.  From the obtained data shown in Table 7 and Figure 3, we observe a decrease in activation energy in the presence of the catalyst, which reflects the efficiency of iron tallates as a pyrolysis catalyst. From another side, the curves of the non-catalytic pyrolysis demonstrate an increase in activation energies with conversion degree evolution. From the obtained data shown in Table 7 and Figure 3, we observe a decrease in activation energy in the presence of the catalyst, which reflects the efficiency of iron tallates as a pyrolysis catalyst. From another side, the curves of the non-catalytic pyrolysis demonstrate an increase in activation energies with conversion degree evolution. However, the catalytic pyrolysis demonstrates almost a constant value with conversion degree evolution, which is explained by the adopted autocatalytic mechanism of the process caused by the presence of iron tallates. The values of activation energy and preexponential factors in Table 8 confirm our hypothesis.  The obtained values of activation energy and preexponential factors obtained from model-based methods demonstrate a similar approach in terms of increasing the energy of activation of the oil pyrolysis process in the presence of iron tallates for all models except the KS model. Broadly speaking, the parameters obtained using the isoconversional approach should coincide with the parameters obtained using the model method as recommended by International Confederation for Thermal Analysis and Calorimetry (ICTAC). Based on this recommendation and taking into account the correlation coefficients, it has been found that the corresponding reaction models and isoconversional methods that adequately describe the process of Cuban extra-heavy oil pyrolysis are the Bna and Cnm models, as well as the Friedman isoconversional method in addition to ASTM E2890. The next figure (Figure 4) demonstrates the correlation between experimental DSC curves obtained for the extra-heavy oil pyrolysis process and the fit DSC curves for this process based on a model corresponding to n-order autocatalysis reaction with m-power by-product (Cnm) in the presence and the absence of iron tallates.

Non-Catalytic Pyrolysis
Interestingly, the obtained results have shed light on the nature of the process of Cuban extra-heavy oil pyrolysis in the presence and the absence of iron tallates. It is worth noting that the obtained Cnm model provides the same results regarding activation energy decreasing in the presence of iron tallates. Moreover, according to the Cnm model, the process of Cuban extra-heavy oil pyrolysis could be described by a set of autocatalytic reactions. In addition, the parameters related to these autocatalytic reactions are well improved by the presence of iron tallates, where the autocatalytic reaction order is greater (4.3) in the presence of iron tallates compared to the non-catalytic pyrolysis (0.01). Besides, autocatalytic reaction parameter m-power increases to 0.15 in the presence of iron tallates compared to when pyrolysis occurs without any catalyst (0.01). Interestingly, the obtained results have shed light on the nature of the process of Cuban extra-heavy oil pyrolysis in the presence and the absence of iron tallates. It is worth noting that the obtained Cnm model provides the same results regarding activation energy decreasing in the presence of iron tallates. Moreover, according to the Cnm model, the process of Cuban extra-heavy oil pyrolysis could be described by a set of autocatalytic reactions. In addition, the parameters related to these autocatalytic reactions are well improved by the presence of iron tallates, where the autocatalytic reaction order is greater (4.3) in the presence of iron tallates compared to the non-catalytic pyrolysis (0.01). Besides, autocatalytic reaction parameter m-power increases to 0.15 in the presence of iron tallates compared to when pyrolysis occurs without any catalyst (0.01).
Regardless of the good effect highlighted by the values of activation energy, autocatalytic order, and power, we have noted that the preexponential factors, however, demonstrate a decrease in values in the presence of iron tallates for all isoconversional and model-based methods. Therefore, to confirm our hypothesis about the efficiency of iron tallates in the process of Cuban extra-heavy oil pyrolysis we calculated the half-life times of the pyrolysis processes in the presence and absence of the catalyst to predict the associated reactions rates, as shown in Figure 5, at a temperature range of 200 °C to 400 °C, which is believed the range at which oil pyrolysis occurs. Regardless of the good effect highlighted by the values of activation energy, autocatalytic order, and power, we have noted that the preexponential factors, however, demonstrate a decrease in values in the presence of iron tallates for all isoconversional and model-based methods. Therefore, to confirm our hypothesis about the efficiency of iron tallates in the process of Cuban extra-heavy oil pyrolysis we calculated the half-life times of the pyrolysis processes in the presence and absence of the catalyst to predict the associated reactions rates, as shown in Figure 5, at a temperature range of 200 • C to 400 • C, which is believed the range at which oil pyrolysis occurs. Interestingly, the obtained results have shed light on the nature of the process of Cuban extra-heavy oil pyrolysis in the presence and the absence of iron tallates. It is worth noting that the obtained Cnm model provides the same results regarding activation energy decreasing in the presence of iron tallates. Moreover, according to the Cnm model, the process of Cuban extra-heavy oil pyrolysis could be described by a set of autocatalytic reactions. In addition, the parameters related to these autocatalytic reactions are well improved by the presence of iron tallates, where the autocatalytic reaction order is greater (4.3) in the presence of iron tallates compared to the non-catalytic pyrolysis (0.01). Besides, autocatalytic reaction parameter m-power increases to 0.15 in the presence of iron tallates compared to when pyrolysis occurs without any catalyst (0.01).
Regardless of the good effect highlighted by the values of activation energy, autocatalytic order, and power, we have noted that the preexponential factors, however, demonstrate a decrease in values in the presence of iron tallates for all isoconversional and model-based methods. Therefore, to confirm our hypothesis about the efficiency of iron tallates in the process of Cuban extra-heavy oil pyrolysis we calculated the half-life times of the pyrolysis processes in the presence and absence of the catalyst to predict the associated reactions rates, as shown in Figure 5, at a temperature range of 200 °C to 400 °C, which is believed the range at which oil pyrolysis occurs.   The obtained half-life times demonstrate higher pyrolysis reaction rates in the presence of iron tallates, especially at lower temperatures (less than 300 • C), which are believed to be the most intense zone of oil pyrolysis processes. Our results are consistent with previous works on heavy oil pyrolysis at larger scales [28][29][30] and our conclusions suggest the same hypothesis about C-S bonds breaking at a temperature less than 300 • C, transforming higher oil components into light fractions and therefore enhancing oil quality and decreasing oil viscosity.

Thermodynamic Functions of Activated Complex Formation
To calculate thermodynamic parameters of activated complex formation, we have used the Eyring equation: where K B and h are the Boltzmann and Planck constants and ∆ ‡ G 0 is the standard Gibbs energy of activation complex formation. The sophisticated Eyring equation is useful for giving insight into the nature of any transition state. In fact, the following equations could be useful for relating the thermodynamic and kinetic parameters obtained by the Arrhenius equation if we consider pyrolysis as a one-step process at each conversion degree: where ∆ ‡ H 0 and ∆ ‡ S 0 are standard enthalpy and entropy of activation complex formation, and T st is a standard temperature which was defined as the peak temperature at the minimum heating rate (716 K). According to Equations (5) and (6), the activation energy is related to the energy difference between reactants and the transition state, and the pre-exponential factor is related to the change in order degree.
To calculate the thermodynamic functions, the Arrhenius parameters derived from KAS and Friedman methods were used. The calculated results are presented in Table 9. Table 9. Thermodynamic parameters for catalytic and non-catalytic pyrolysis processes based on KAS and Friedman methods.

Pure Oil
Oil with Catalyst α (%) From Table 9, we can see that thermodynamic functions calculated based on Arrhenius parameters derived from the KAS and Friedman methods converge well with each other within uncertainties. Comparison of the activation enthalpies for the process in the presence and absence of the catalyst indicates that the use of iron tallates significantly reduces the energy consumption for the formation of an activated complex and also makes it possible to avoid a significant increase in the activation enthalpy as pyrolysis proceeds. However, at the same time, when using the catalyst, the activation entropy is significantly reduced, which indicates a greater degree of arrangement of the transition state. Therefore, at high temperatures, the beneficial decrease in enthalpy will be largely suppressed by the entropy contribution (∆ ‡ G 0 is about the same for both processes). At the same time, at lower temperatures, the entropy contribution will be much smaller, so the process in the presence of iron tallates will proceed much faster (see Figure 5).

Conclusions
Taken together, the results of the present study have shed light on the kinetic behavior of Cuban oil pyrolysis in the presence and absence of iron tallates. During this study, iron tallates were obtained and used as a material of investigation alongside oil samples. We applied thermal analysis to study the kinetic behavior of oil under different heating rates. The obtained SARA analysis has shown a significant increase in the content of light hydrocarbons (saturates and aromatics) during the oil catalytic aquathermolysis reactions. Moreover, the elemental analysis has indicated an increase in C and H content by almost 4% and 6%, respectively, with a significant decrease in S and O content by up to 23% in the presence of iron tallates. These results have been further confirmed by infrared spectrometry, which evidenced the conversion of asphaltene and resin fractions into light hydrocarbons. Moreover, our results have indicated a significant decrease in activation energy in the presence of iron tallates. Moreover, our hypothesis about catalyst efficiency has been confirmed further by predicting the pyrolysis reaction rate in the presence and absence of iron tallates via calculating pyrolysis half-life times at 200 • C to 400 • C and finding that iron tallates increase the pyrolysis process rate in a significant manner, especially at lower temperatures, which are supposed to be the temperatures of pyrolysis initiation and evolution. In addition, the evidence from this study suggests that the process of Cuban extra-heavy oil pyrolysis includes reactions of n th order with m-power autocatalysis by-product. The evidence from the thermodynamic study suggests that the adopted catalyst for the process of Cuban extra-heavy oil pyrolysis reduces the enthalpy and entropy of activation, which may be useful for energy consumption optimization and process control. We believe that these findings add to a growing body of literature on extra-heavy oil upgrading technologies via thermal enhanced oil recovery methods. The obtained results are encouraging and may improve the application of steam injection processes, especially in the presence of oil-soluble catalysts.