Type I kerogen-rich oil shale from the Democratic Republic of the Congo: mineralogical description and pyrolysis kinetics

The Democratic Republic of the Congo holds important reserves of oil shale which is still under geological status. Herein, the characterization and pyrolysis kinetics of type I kerogen-rich oil shale of the western Central Kongo (CK) were investigated. X-ray diffraction, Fourier-transform infrared spectroscopy and thermal analysis (TG/DTA) showed that CK oil shale exhibits a siliceous mineral matrix with a consistent organic matter rich in aliphatic chains. The pyrolysis behavior of kerogen revealed the presence of a single mass loss between 300 and 550 °C, estimated at 12.5% and attributed to the oil production stage. Non-isothermal kinetics was performed by determining the activation energy using the iterative isoconversional model-free methods and exhibits a constant value with E = 211.5 ± 4.7 kJ mol−1. The most probable kinetic model describing the kerogen pyrolysis mechanism was obtained using the Coats–Redfern and Arrhenius plot methods. The results showed a unique kinetic triplet confirming the nature of kerogen, predominantly type I and reinforcing the previously reported geochemical characteristics of the CK oil shale. Besides, the calculation of thermodynamic parameters (ΔH*, ΔS* and ΔG*) corresponding to the pyrolysis of type I kerogen revealed that the process is non-spontaneous, in agreement with DTA experiments.


Introduction
In order to satisfy the global energy demand and overcome the regression of conventional oil and natural gas reserves, oil shale constitutes a promising alternative resource from which fossil fuels can be produced (Ngo and Natowitz 2016). Oil shale is a sedimentary rock comprising a large portion of solid organic matter (OM) called kerogen, heterogeneously distributed in a major mineral matrix composed mainly of silicates, carbonates and clays (Yen and Chilingarian 1976). The kerogen part which is a three-dimensional macromolecular complex and insoluble in ordinary solvents can be converted by thermal treatment up to certain critical temperature into oil, gas and residual carbon (Demirbas 2016). Oil shale deposits are found all over the world, the largest being in USA, China, Brazil and Estonia (Dyni 2003), and several industrial plants were taken place to generate electrical energy by direct combustion or high-value fuels by pyrolysis (Dyni 2003;Yen and Chilingarian 1976). The Democratic Rebublic of the Congo (DRC) country also has important oil shale reserves, ranked first in the African continent but it is still under the geological exploration status (Delvaux and Fernandez-Alonso 2015).
The fundamental study of the pyrolysis of oil shale kerogen in laboratory conditions is crucial as it provides relevant information on hydrocarbon potential and predicts thermal behavior in order to design an efficient operating protocol for large-scale applications. Therefore, several pyrolysis processes were the subject of extensive investigation in the past decades (Bai et al. 2017;Maaten et al. 2018;Saif et al. 2017). The composition and yield of Edited by Xiu-Qiu Peng pyrolysis products depend on many factors, such as final temperature, heating rate and particle size. (Geng et al. 2017;Jiang et al. 2015;Maaten et al. 2018). Therefore, these factors should be optimized through the understanding of reaction mechanisms involved during pyrolysis and the determination of kinetic parameters describing the process as the activation energy, pre-exponential factor and reaction model (Li and Yue 2004;Wang et al. 2013). Several tools have been used to optimize the energy potential of oil shale pyrolysis, including thermogravimetric analysis (TGA), which allows to study in detail the thermal kinetics of different processes (Khraisha and Shabib 2002;Tiwari and Deo 2012a;Wang et al. 2014). For the most part, the main mass losses monitored during non-isothermal oil shale pyrolysis consist of moisture removal at low temperatures, kerogen degradation at around 200 °C-600 °C and clays and carbonates minerals decomposition at above 700 °C (Maaten et al. 2016;Strizhakova and Usova 2008).
In the literature, the kerogen pyrolysis has been reported as a complex reaction process depending on the nature of oil shale studied and its preservation environment (Maaten et al. 2016). Generally, the pyrolysis mechanism is described by two or more steps with different kinetic parameters (Bai et al. 2015a; Moine et al. 2016;Ren et al. 2015;Wang et al. 2015). Using model fitting methods, several reports admitted one or multiple first-order reaction kinetic models to describe the oil shale pyrolysis processes (Maaten et al. 2016;Wang et al. 2015). In this context, Maaten et al. (2016) found two steps in the kerogen pyrolysis process of American, Chinese and Estonian oil shales and adopted the first-order reaction for each one to estimate the activation energies and preexponential factors. On the other hand, using model-free isoconversional methods, many authors have found a variable activation energy tendency with the reaction advancement due to the multistep character of oil shale pyrolysis (Bai et al. 2015b;Janković 2013). They reported nth-order reaction model for oil shale pyrolysis mechanism, with n equal or different from unity (Moine et al. 2016;Wang et al. 2015). According to Janković (2013), the pyrolysis of Brazilian oil shale was found to follow nth-order reaction model with n = 1.23. While, Moine et al. (2016) used the isoconversional Model Free methods followed by master plots Fitting method to examine the three steps mechanism of pyrolysis of Moroccan Rif oil shale. The authors concluded that all processes follow the nth-order reaction model, with n = 1 for the first and third processes, corresponding to the pyrolysis of bitumen and pyrite, respectively, and n = 1.6 for the second process, which corresponds to the decomposition of kerogen. However, the pyrolysis of oil shale containing type I kerogen remains very limited as this type of material is very rare, except those of American and Turkish oil shales for which the thermal behavior generally leads to a single process with first-order reaction model and with an activation energies ranging between 155.2 and 363.2 kJ mol −1 (Karabakan and Yürüm 1998;Tiwari and Deo 2012b).
In the present work, the pyrolysis kinetics of oil shale originating from the western Central Kongo Formation of the DRC was investigated to illustrate the degradation mechanism governing the thermal process of OM and to complete the mineralogical and type I kerogen characterization, recently described in our previous report ). The pyrolysis kinetic study was monitored using TG instrument under non-isothermal conditions and the iterative forms of isoconversional model-free methods to estimate the activation energy with conversion degree, followed by Coats-Redfern and Arrhenius plot methods to get the exact reaction model f(α) and the pre-exponential factor A. Finally, the thermodynamic parameters (ΔH*, ΔG* and ΔS*) related to the thermal degradation of kerogen were determined using the activated complex theory. To our knowledge, the pyrolysis kinetic of western Central Kongo oil shale has never been the subject of a previous study.

Oil shale samples
The oil shale (CK) samples used in the present study were collected from the western part of the Central Kongo as shown in Fig. 1  , by our partners in the Center of Geological and Mining Research (CGMR) of Kinshasa (DRC). The sampling region is part of the Neoproterozoic PanAfrican-Brasiliano Orogenic Belts that contain two N-S-trending orogens (Cailteux et al. 2015). It takes the core of the Central African shield and can be determined as an assembly of Archean fragments, including the Kasai, Mbomou and the Chaillu-Gabon cratons, which were once overtaken with the Sao Francisco (Brazil) craton, to form the amalgamated Central African landmass at the time of Gondwana assembly (Daly et al. 1992;Roberts et al. 2015). Samples were taken from unexposed layers to air in order to avoid weathering impacts in an area with great similarity of rock deposits. Before the experiments, the sample was ground and sieved to the size of less than 180 μm (80 mesh), according to the International Confederation for Thermal Analysis and Calorimetry (ICTAC) committee recommendations and American ASTM standard (E11) (Vyazovkin et al. 2014). The Rock-Eval parameters and mineralogical analysis of the CK samples obtained in our previous geochemical study are summarized in Table 1  . As can be seen from the results in Table 1, the OM of CK oil shale presented a very good organic richness (TOC = 10.77%) with a very good petroleum potential (S2 = 80.4 mgHC/g rock). In addition, the elemental analysis showed a high value of the H/C molar ratio (1.71), which corresponded to a sapropelic OM with type I kerogen. These findings were corroborated by the high value of hydrogen index (HI = 746 mgHC/g TOC) and low value of oxygen index (OI = 13 mgCO 2 /g TOC). The S1 parameter which was equal to 2.7 mgHC/g rock showed a low content of organic volatiles like bitumen. From these results, we can conclude that the shale rock studied presents a very promising petroleum potential and could be of significant commercial interest by comparing with other reported oil shales (Tissot and Welte 1984;Muhammad et al. 2011).

Characterization
X-ray diffraction (XRD) analysis was performed at room temperature using a Shimadzu 6100 diffractometer equipped with a copper anticathode (λ CuKα = 1.541838 Å) with a step scan of 0.02 °s −1 , from 5 to 70 (°2θ) and at (40 kV and 30 mA). The identification of minerals in the samples was performed by comparing the obtained XRD patterns with those of the Joint Committee on Powder Diffraction Standards (JCPDS-ICDD) database. Fourier-transform infrared (FTIR) spectroscopy was carried out on a Jasco FT/IR 4600 spectrometer equipped with ATR (Pro One) module. FTIR spectra of the sample were obtained at 4 cm −1 resolution and collected in the range of 4000-400 cm −1 . Simultaneous thermogravimetry/differential thermal analysis (TG/ DTA) was carried out using Labsys™ Evo (1F) SETARAM equipment. ICTAC kinetics committee recommendations for collecting experimental data and performing kinetics computations were used to evaluate the kinetic parameters (Vyazovkin et al. 2011(Vyazovkin et al. , 2014. The thermal analysis system was standardized for temperature reading with reference  1 3 metals of 99.99% purity to calibrate the microbalance by minimizing the buoyancy effects for better estimation of mass changes. The pyrolysis experiments were performed with a mass of 10 mg (± 0.1 mg) of oil shale from room temperature to 1000 °C at different heating rates of 5, 10, 15 and 20 °C min −1 and under a nitrogen flow of 45 mL min −1 . The mass of sample was placed in alumina pan with well dispersion and low depth to reduce heat and mass transfers. The thermal experiments were performed in triplicate to ensure repeatability and reproducibility of the obtained data.

XRD analysis
The oil shale is a complex mixture of organic and mineral substances forming a very hard and nonporous material. The XRD pattern of the CK oil shale and its calcined form at 900 °C (CK 900) are displayed in Fig. 2. The CK oil shale is found to be predominantly containing quartz with main observed diffraction peaks corresponding to JCPDS Card file No. 01-079-1910. Small amounts of clay minerals, such as montmorillonite, illite and muscovite, were also detected which are indicative of a great hydrocarbon potential of the source rock (Wei et al. 2006). On the other hand, no traces of carbonate minerals such as calcite and dolomite were detected. This result is rarely found in known oil shales, but similar to those of Puertollano Spanish oil shale and Arbaa Ayacha Moroccan oil shale (Moine et al. 2018;Torrente and Galán 2001). Additionally, a broad continuum is observed in the XRD pattern, which is certainly due to the presence of an amorphous phase, corresponding to the complex OM in oil shale (Altun et al. 2009). The appearance of pyrite at 26.18 and 33.32° is frequently observed in various carbon-rich sedimentary rocks (Gai et al. 2014) and identified by SEM microscopy in our previous report . The presence of pyrite in CK oil shale indicates the anoxic environment of the source rock, in agreement with marine or lacustrine deposition (Dyni 2003). The calcined CK oil shale at 900 °C (CK 900) under air atmosphere is mostly composed of quartz with the subsistence of pyrite. At high temperature, clay minerals are destroyed into oxides and OM volatilizes.

FTIR-ATR spectroscopy
FTIR spectra of CK and CK 900 samples are shown in Fig. 3. As can be seen, the CK sample exhibits a complex spectrum formed by the superposition of bands due to the OM and to the mineral matrix. Therefore, the OM was identified by the main unsaturated and saturated aliphatic C-H absorption bands around 2930, 2800 and 1455 cm −1 , revealing the dominance of sapropelic kerogen (Solomon and Miknis 1980). The existence of aromatic compounds is excluded because of the absence of their characteristic narrow band at about 3000 cm −1 (Alstadt et al. 2012). On the other side, the occurrence of small shoulder at 1715 cm −1 indicates the probable presence of carbonyl kinds ). The mineral matter was revealed by the absorption bands at 3703 and 3265 cm −1 which are assigned to OH groups linked to Si and Al of clays. Other absorption bands which appear around 3700 and 1630 cm −1 were attributed to internal water molecules of mineral clays (Palayangoda and Nguyen 2012). The multitude of bands observed in the region 1125-420 cm −1 were assigned to different vibration Si-O modes of quartz (Palayangoda and Nguyen 2012 (type I kerogen) and corroborates our previous Rock-Eval analysis ( Table 1).
As for the calcined sample of CK oil shales (CK 900), the FTIR spectrum shows only the absorption bands due to Si-O and Al-O. However, a slight shift of bands toward high frequencies is observed. This is due to the structural arrangement of silica and alumina after the thermal degradation of organic matter and complete dehydroxylation of clay minerals ). Figure 4 shows the thermal behavior of CK sample from room temperature up to 1000 °C, at heating rate of 10 °C min −1 under nitrogen gas flow (45 mL min −1 ) and with initial mass of 10.0 ± 0.1 mg. The presented TG, derivative TG (DTG) and DTA profiles show three mass losses. At temperature less than 200 °C, the observed mass loss of about 4% and its corresponding endothermic peak are attributed to the removal of moisture as well as water molecules associated with the mineral matrix. The second and main mass loss occurring between 300 and 550 °C, evaluated at 12.5% and accompanied by an endothermic effect in DTA, corresponds to the pyrolysis of kerogen present in oil shale to form oil, gas and residual carbon (Yan et al. 2013). This second stage is considered as the important pyrolysis step in processing oil shale. The temperature corresponding to the maximum pyrolysis rate and evaluated at 469 °C is similar to Kentucky and Green River oil shales (Maaten et al. 2016), but lower than many world oil shales (Janković 2013). According to the literature, low maximum decomposition temperature informs about the nature of OM and probably indicates a high kerogen reactivity in agreement with type I kerogen (Maaten et al. 2016). In the temperature range between 550 and 650 °C, the low mass loss observed is assigned to the decomposition of pyrite contained in CK oil shale (Chen et al. 2000;Gai et al. 2014). The absence of distinct mass loss beyond 600 °C justifies the absence of carbonates and confirms the XRD and FTIR characterization results.

Theoretical basis
The general theoretical basis of thermal decomposition kinetics in solid state is expressed by the Arrhenius equation as follows (Khawam and Flanagan 2006;Vyazovkin et al. 2011): where α is the reacted fraction, expressed as = m 0 − m t ∕ m 0 − m ∞ , k is the rate constant, E a (kJ mol −1 ) is the apparent activation energy, and A (min −1 ) is the pre-exponential factor, R is the gas constant (R = 8.314 J mol −1 K −1 ) and f(α) is the reaction model function that gives crucial information about the mechanism. The most known functions f(α) used in solid-state kinetics as well as their corresponding integral forms g(α) are recapitulated in Table 2.
The combination A, E a and f(α) is called the kinetic triplet. In case of non-isothermal modes, the temperature varies at a constant rate β, given by = dT∕dt . Equation (1) can then be transformed into: The function g(α) can be obtained by integrating Eq. (2) as follows: The function p(x) has no analytical solutions but is solved using different mathematical approximation methods (Flynn 1997).
3.4.1.1 Determination of apparent activation energy E a According to ICTAC recommendations, the most accurate procedure to determine the activation energy in solid-state kinetics is the isoconversional method, which enables the calculation of activation energy without prior knowledge about the reaction model (Farjas and Roura 2011;Vyazovkin et al. 2011). This method is based on the isoconversional principle which states the linear depend- . 4 Pyrolysis profiles of CK oil shale at heating rate of 10 °C min −1 ence between the reaction rate and temperature at constant conversion degree. Furthermore, it has the advantage to detect the presence of eventual complicated processes (Moine et al. 2016).
In this study, the iterative form of the integral methods of Ozawa-Flynn-Wall (OFW) and Kissinger-Akahira-Sunose (KAS) was used to calculate the apparent activation energy as function of the conversion degree α.
The iterative integral methods are used to allow more precise approximation of the temperature integral [Eq. (3)]. These methods assume that in a slight variation in temperature, the corresponding activation energy is almost constant (Deng et al. 2009). Among the most used, we find: The Ozawa-Flynn-Wall method (OFW) is based on Doyle's approximation to solve the temperature integral (Cai and Chen 2012;Flynn and Wall 1966). Rearrangement of Eq. (3) yields the following expression: Its iterative form IT-OFW is written as follows: The Kissinger-Akahira-Sunose method (KAS) is based on Murray and White's approximation to approach the temperature integral (Cai and Chen 2012;Liu et al. 2009). After rearrangement of Eq. (3), the KAS method can be written as follows: Its iterative form IT-KAS can be expressed as: where h(x) = x 4 +18x 3 +86x 2 +96x x 4 +20x 3 +120x 2 +240x+120 Plotting the left side of Eq. (5), ln(β/H(x)) versus (1/T), and of Eq. (7), ln(β/h(x)T 2 ) versus (1/T), at each conversion Reaction models Symbol f(α) g(α) Chemical process or mechanism non-invoking equations First-order F1 ( degree α, leads to an estimate activation energy E α which is deduced from the slope of the obtained straight lines. (iii) Iteration procedure The iteration procedure is realized by assuming that h(x)=1 or H(x)=1 to estimate the initial value of activation energy E α,1 . The energy calculation using the isoconversional methods gives then the value of E α,1 corresponding to the conventional energy of the integral methods. This energy is then used to calculate new values of h(x) and H(x) and thus E α,2 . This operation is repeated continuously by replac- value of E ,i is then the exact value of activation energy of the pyrolysis process.

Reaction model determination
The conversion function f(α) is the mathematical and theoretical description of the decomposition reaction mechanism (Table 2). Different methods were developed to evaluate the kinetic model describing the experimental data (Vyazovkin et al. 2011), upon which is the Coats-Redfern method that was developed from Eq. (3) using the temperature integral approximation of Coats and Redfern (1964) and different reaction models f(α). It is expressed as: The plot ln[g(α)/T 2 ] against 1/T should give a straight line if the proposed reaction model corresponds to the pyrolysis mechanism involved.
In order to refine the kinetic parameters obtained by Coats-Redfern method, the direct Arrhenius plot method was used. Considering the general kinetic equation and after rearrangement, Eq. (2) becomes: The linear plot of Eq. (9),ln 1 T , is obtained after supposing the reaction model f(α). The slope and intercept will provide the activation energy E and the pre-exponential factor A, respectively. The most probable reaction model f(α), which describes the pyrolysis process, will have an activation energy near the one obtained from the iterative isoconversional methods.

Pre-exponential factor determination
The preexponential factor A value was estimated using the equation proposed by Málek (1989) [Eq. (10)] and compared to the values obtained by Arrhenius plot method at each heating rate: where T max is the temperature and α max is the conversion degree at the maximum reaction rate.

Thermodynamic parameters
The thermodynamic parameters of CK oil shale pyrolysis were estimated using the equation developed by the transition state theory (Font and García 1995). They are expressed using the following equations: where e = 2.7183 is the Neper Number, χ is the transmission factor, which is equal to unity for monomolecular reactions, k B is the Boltzmann constant (1.381 × 10 −23 J K −1 ), h is the Plank constant (6.626 × 10 −34 J s), T p is the temperature at the maximum reaction rate in DTA curve and E a is the average activation energy obtained from IT-KAS method.

Thermal kinetics of CK oil shale
Kinetic study of CK oil shale pyrolysis was conducted for the main thermal step corresponding to the OM decomposition as illustrated in Fig. 4. Thermal measurements were realized under non-isothermal conditions using four sets of heating rates (5, 10, 15 and 20 °C min −1 ), in accordance with ASTM standard (E11) and ICTAC kinetics recommendations. The obtained curves of conversion degree α and reaction rate dα/dt as function of temperature are shown in Fig. 5. As can be seen, all thermal curves present a single thermal transformation corresponding to the decomposition process of OM to hydrocarbons. Similar trend of singlestep pyrolysis was also observed in other oil shales, such as Green River oil shale and Kentucky oil shale (Maaten et al. 2016(Maaten et al. , 2017Tiwari and Deo 2012b). The single-step mechanism in pyrolysis of oil shales is uncommon and related to the type of kerogen (Williams and Ahmad 2000). In fact, the type I kerogen which is predominantly composed of aliphatic hydrocarbons can lead to the same thermal behavior, while the type II and type III kerogens exhibit complex behaviors because of the additional presence of bitumen and aromatics in oil shales (Bai et al. 2015a;Moine et al. 2018). The displacement of the thermal curves toward higher temperatures when increasing the heating rates is a common result of several previous studied world's oil shales (Bai et al. 2015a;Li and Yue 2004). This is due to the temperature gradient inside and outside the particles of OM, which is lower at low heating rates and higher at longer exposure of sample to a particulate temperature (Sánchez-Rodríguez et al. 2014), (Table 3).

Activation energy determination
The iterative isoconversional methods IT-OFW and IT-KAS [Eqs. (5) and (7)] were used to evaluate the distribution of apparent activation energy during the pyrolysis process. The (E α − α) curves illustrated in Fig. 6 show comparable behavior for both iterative methods. In the conversion range between 0.2 and 0.8, the apparent activation energy is practically constant and its variation does not exceed 5%. This result reflects the presence of a unique OM degradation process, characterized by a single kinetic triplet. A similar finding was also obtained on Huadian oil shale (Qing et al. 2009). The average activation energy corresponding to the OM pyrolysis process of CK oil shale is estimated to be 211.5 ± 4.7 kJ mol −1 . This value is consistent with the energies obtained during the type I kerogen pyrolysis of the Green River, Colorado and Huadian oil shales which were evaluated to be 201, 219.4 and 231 kJ mol −1 , respectively (Campbell et al. 1978;Han et al. 2015;Maaten et al. 2017).

Determination of reaction mechanism
The reaction mechanism is a description of the pathway by which the decomposition reaction occurs. It is generally expressed by a theoretical and mathematical equation called the reaction model. Hence, for CK oil shale, the reaction model was determined using the Coats-Redfern procedure [Eq. (6)].
The linear regression of ln[g(α)/T 2 ] versus (1/T) for different heating rates gives the results summarized in Table 4. Many reaction models show good correlation with the experimental data (R 2 = 0.999), whereas the calculated activation energies are different between the supposed models. When comparing these values with the average activation energy obtained from the isoconversional IT-KAS method, evaluated to be 211.5 kJ mol −1 , it enables identification of the suitable kinetic model which has the closest value to the isoconversional value. Therefore, since the IT-KAS activation energy is found between the mean energy values of F1 and F2 models (164.1 kJ mol −1 and 247.9 kJ mol −1 , respectively) compared to other models (Table 4), the mechanism that describes the pyrolysis kinetics of CK oil shale is more probably following the Fn model, with n value ranging from 1 to 2. This Fn model has also been reported to describe the mechanism of many oil shales' pyrolysis or combustion (Kök and Pamir 2000;Kuang et al. 2018;Ma and Li 2018).
To obtain the exact value of n, the Arrhenius linear plot method (Eq. (9)) was used, with the plot of ln 1 f ( ) d dT versus 1000/T (Fig. 6). According to the results obtained by Coats-Redfern method, the reaction model f(α) used was (1 − α) n with n estimated in the range of 1 and 2. The reliable and precise reaction order n was assumed to generate the activation energy equal to that of IT-KAS isoconversional method. Table 5 displays the calculated values of E and lnA of CK oil shale pyrolysis as function of n using the experimental data. Figure 7 shows the Arrhenius linear curves obtained with the exact value of the reaction order n at which the slope gives an activation energy equal to that estimated by the IT-KAS method. The obtained results at different  The pre-exponential factor values obtained by the Arrhenius linear plot method at different heating rates were compared to the calculated values using Eq. (10). The obtained values of lnA are close to each other for all heating rates (Table 6) and confirm the proposed reaction order model Fn (n = 1.275). The mean value of A was evaluated at 3.23 10 14 , which is a high value that indicates the presence of "loose" complex process during pyrolysis (Cordes 1968). The same magnitude was observed for most world oil shales degradations (Janković 2013;Maaten et al. 2016;Moine et al. 2016). The high pre-exponential values were also found for the complex decomposition of polymeric organic macromolecules such as poly(tetrafluoroethene) and chitin (Genieva et al. 2010;Georgieva et al. 2012).
The comparison of kinetic parameters of CK oil shale with other type I kerogen oil shale is illustrated in Table 7. Although the pyrolysis kinetic approach of CK oil shale is different compared to other oil shale reported in the literature, the experimental kinetic parameters are in good agreement. From these results, it can be concluded that thermogravimetric pyrolysis is a useful complementary technique for studying the nature and composition of kerogen in oil shales.

Modeling of experimental data
The above-obtained kinetic parameters at each heating rate were substituted in Eq. (1) to simulate the experimental reaction rate curves of CK oil shale pyrolysis. Figure 8 displays the plots of the experimental and the obtained theoretical data at different  and volatile bitumen present in CK oil shale, in agreement with the Rock-Eval analysis. The kinetic triplet found by the methods used in this work leads to a satisfactory modeling of the experimental data of CK oil shale pyrolysis.

Thermodynamic parameters
The application of transition state theory to CK oil shale pyrolysis was performed through the estimation of thermodynamic parameters using Eqs. (11)-(13). The activation energy used was the mean value of IT-KAS method. The pre-exponential factor A and reaction order n were taken from the Arrhenius plot method. The obtained average results are summarized in Table 8. The positive values of activated entropies ΔS* indicate the presence of malleable activated complexes and large number of freedom degrees for rotation and vibration (Shannon 1964). The positive values of activated enthalpies ΔH* are consistent with the endothermic effect observed in CK oil shale pyrolysis. However, the positive values of free Gibbs energies ΔG* suggest thermodynamically stable nature of CK oil shale under atmosphere pressure and its decomposition occurs upon heating.

Conclusion
The physicochemical characterization of Central Kongo oil shale (CK) was carried out using XRD, FTIR and TG/DTA techniques and showed that it mainly consists of quartz and organic matter (OM) with sapropelic chains (Type I kerogen). The non-isothermal pyrolysis of CK sample was investigated using thermogravimetric analysis, and the kinetic parameters were estimated using the iterative isoconversional IT-KAS and IT-OFW methods followed by the Coats-Redfern and Arrhenius plot methods. The main pyrolysis step was observed in the temperature range of 300 °C-550 °C and characterized by a single mass loss of 12.5% which corresponds to degradation of type I kerogen.
Using isoconversional methods, the apparent activation energy of the pyrolysis, evaluated by IT-KAS method, was found constant in the main range of conversion [0.2-0.8] and equal to 211.5 ± 4.7 kJ mol −1 , in compatibility with the type I kerogen. The pyrolysis mechanism of kerogen was estimated using Coats-Redfern followed by Arrhenius plot methods, and the results give a single kinetic triplet with Eα = 210.0 kJ mol −1 , lnA = 34.26 and f(α) = (1 − α) 1.275 . These values are consistent with type I kerogen which is formed mainly by sapropelic chains. On the other hand, the thermodynamic parameters evaluated were fully correlated with the endothermic effect and the non-spontaneous nature of the CK oil shale during pyrolysis.