Unconventional resources such as tight, fractured and hybrid shale gas and oil plays as well as oil or kerogen shale systems, are considered exploitable self-contained source and reservoir rocks. A better understanding of the thermal cracking of sedimentary organic matter, hydrocarbons generation, expulsion, storage and retention mechanisms constitutes a key point, estimating the oil and gas in-place, free or adsorbed, for their exploration and exploitation. Herein, we introduce a new “ready to use” method of analysis and interpretation for the Rock-Eval 6 device for better assessment of free or sorbed hydrocarbons in unconventional shale plays.
Maria-Fernanda Romero-Sarmiento1*, Daniel Pillot1, Géremie Letort1, Violaine Lamoureux-Var1, Valérie Beaumont1, Alain-Yves Huc2 and Bruno Garcia1
1IFP Energies nouvelles (IFPEN), Direction Géosciences, 1-4 avenue de Bois-Préau, 92852 Rueil-Malmaison Cedex - France. 2UPMC – Université Pierre et Marie Curie, Institut des Sciences de la Terre Paris, UMR 7193 UPMC-CNRS, 4 place Jussieu, 75005 Paris - France.
Manuscript submitted in January 2015, Manuscript accepted in March 2015, Published online in August 2015.
This method was developed at IFP Energies nouvelles (France) and was tested on 15 actual or potential unconventional shale samples from Silurian Shale (Algeria), Mississippian Barnett Shale (USA), Early Jurassic Shale (France), Late Jurassic Bazhenov Shale (Russia) and Eocene Green River Shale at different thermal maturity stages. Results indicate a better quantiﬁcation of free and/or sorbed hydrocarbons (Sh0 and Sh1 peaks) as well as a more accurate determination of the Rock-Eval Tmax maturity parameter.
Over the last 30 years, the Rock-Eval pyrolysis has been widely used to identify organic matter occurrence, type and thermal maturity level. This technique has been also used to quantify the Total Organic Carbon (TOC) content as well as the mineral carbon content of both reservoir and source rocks (Espitalié et al., 1986; Lafargue et al., 1998; Behar et al., 2001). The Rock-Eval 6 device is equipped with two ovens for pyrolysis and combustion processes, respectively. The hydrocarbons generated during a Rock-Eval analysis are monitored by a Flame Ionization Detector (FID) whereas the non-hydrocarbons compounds like CO2 and CO released during pyrolysis and oxidation stages are monitored by an infra-red detector (Behar et al., 2001). This pioneering device, thanks to the built-in “Basic/Bulk-Rock” and “Pure Organic Matter” methods, provides relevant parameters to respectively characterize any source rocks and isolated kerogens, in a petroleum system perspective. Applied to exploration purposes, the S2 peak, which is assigned to the hydrocarbon residual potential, is the main focus of these methods. It is widely accepted that the quantity of hydrocarbons released beyond 300 C is associated to the thermal cracking of the organic matter (kerogen) and is then refereed as the pyrolysis stage from Rock-Eval analyses. During this pyrolysis stage, hydrocarbons detected between 300º and 650º C or 300º and 800ºC, according to user choice, correspond to the S2 peak surface from typical “Basic/Bulk-Rock” and “Pure Organic Matter” analyses. In addition, the classical Rock-Eval thermal maturity parameter Tmax is calculated from the temperature at which the S2 peak reaches its maximum.
More details about the corresponding pyrolysis and oxidation conditions for these classical Rock-Eval methods are described in Lafargue et al. (1998) and Behar et al. (2001). When the S1 peak is the main concern, e.g. for reservoir units studies, an alternative built-in method, the “Reservoir” method, is available. This method is usually performed on cutting or core samples and is based on a pyrolysis temperature program allowing an improved recovery of “free” hydrocarbons, as a result of a lower initial temperature, and a better deconvolution of the released components.
Moreover, the “Reservoir” method was shown to provide an estimation of API gravity of hydrocarbon ﬂuids from cuttings (Trabelsi et al., 1994). The “Reservoir” method pyrolysis temperature program is also designed to estimate the occurrence and nature of heavy oils and tarsmats intervals during drilling operations (Carpentier et al., 1995, 1998). For the “Reservoir” method, the initial starting pyrolysis temperature is reduced to 180ºC, in order to monitor the low-molecular weight hydrocarbons, which are not considered with the Basic/Bulk-Rock method as it starts at 300ºC. The lower initial starting temperature is combined with a longer temperature plateau of 10 minutes for a more complete capture of the low-molecular weight hydrocarbons, resulting in a more realistic quantiﬁcation of the hydrocarbons associated to the S1 peak, comparing to the Basic/Bulk-Rock method. The “Reservoir” method therefore provides two peaks designated as S1r and S2a, which in turn can be calculated to produce a “S1 peak” more equivalent to the conventional S1 peak from the Basic/ Bulk-Rock method for reservoir oil and tar samples (Trabelsi et al., 1994). S1r represents the very lightest fractions detected during the temperature plateau at 180ºC whereas the S2a peak detects hydrocarbon compounds between 180º and 325ºC (Trabelsi et al., 1994). More details about the corresponding pyrolysis and oxidation conditions for the “Reservoir” method are described in Trabelsi et al. (1994).
Meanwhile, the increasing interest in unconventional gas and oil shale plays requires the development of a new method providing relevant parameters devoted to exploration and production of these new plays. The proposed method is designed to consider the speciﬁcity of the dual attributes of unconventional gas and oil shale plays: being a source rock and a reservoir. Herein, we introduce a speciﬁc IFPEN Shale Play method (Patent 14/55.009, Pillot et al., 2014) for any Rock-Eval 6 device, including a new pyrolysis program and associated parameters for evaluation of unconventional shale resources such as tight, fractured and hybrid shale gas and oil plays as well as oil or kerogen shale systems. This new method was elaborated and tested from 15 worldwide unconventional shale samples, i.e. Silurian Shale (Algeria), Mississippian Barnett Shale (USA), Early Jurassic Shale (France), Late Jurassic Bazhenov Shale (Russia) and Eocene Green River Shale (USA). This speciﬁc Rock-Eval method attempts to optimally quantify hydrocarbons still present within unconventional dual source/reservoir rocks.
1 Material and methods
1.1 Description: The IFPEN Shale Play Method
The IFPEN Shale Play method is characterized by a speciﬁc temperature program, which allows a more exhaustive recovery and separation of the thermovaporizable hydrocarbons (which tentatively can be assigned to both free and sorbed hydrocarbons in shale plays). Figure 1 displays the thermovaporization and pyrolysis conditions deﬁned for this new IFPEN Shale Play method in comparison to the existing Reservoir and Basic/Bulk-Rock methods:
Figure 1 Rock-Eval temperature program for hydrocarbons characterization for the new IFPEN Shale Play method compared to the classical “Basic/Bulk-Rock” and “Reservoir” methods.
The speciﬁc analytical conditions for the IFPEN Shale Play method are described as follow:
- the pyrolysis step starts at T1 = 100ºC. This temperature is chosen as the most appropriate to initiate thermovaporization and capture eventually the quantity of the lightest hydrocarbons, still present in the sample;
- from initial T1, the temperature is programmed to increase at a rate of 25ºC/min, up to T2 = 200ºC. T2 is maintained as a plateau for 3 minutes. During this step, the more easily thermovaporizable hydrocarbons are released and monitored;
- from T2 plateau, the temperature is then increased at a rate of 25ºC/min, up to T3 = 350ºC. This temperature is maintained for a plateau of 3 minutes. During this step heavier molecular weight compounds are thermovaporized. T3 is assumed to correspond to the limit for onset of kerogen thermal cracking;
- from T3, the temperature is raised again according to a linear temperature program deﬁned by a rate of 25ºC/min,
up to a ﬁnal temperature T4 = 650ºC, designed to release the pyrolyzable part of the remaining organic matter.
Depending of the scientiﬁc purpose and the type of samples to be analyzed, it is mandatory for the user to select the appropriate Rock-Eval method to ensure an useful meaning of the Rock-Eval results. Table 1 summarizes different types of geological samples that can be run using the currently available Rock-Eval methods.
1.2 New Acquired Parameters from the IFPEN Shale Play Method
The IFPEN Shale Play method provides 3 key parameters: Sh0, Sh1 and Sh2. These parameters are derived from the surface areas of the 3 recorded peaks (Fig. 2). They correspond to the quantities of HC compounds monitored by a FID. These acquired parameters are expressed in milligrams of HC compounds per gram of initial rock, as provided by all Rock-Eval methods. Figure 2 also displays the Rock-Eval pyrograms comparing the 3 now available Rock-Eval methods: Basic/Bulk-Rock, Reservoir and Shale Play.
Figure 2 Rock-Eval pyrograms comparing the 3 Rock-Eval methods: Basic/Bulk-Rock, Reservoir and Shale Play.
The Sh0 peak is assigned to the lightest thermovaporized hydrocarbons released between T1 (100ºC) and T2 (200ºC). The Sh1 peak is assigned to heavier thermovaporized hydrocarbons released between T2 (200ºC) and T3 (350ºC). Finally, the Sh2 peak is assigned to HC issued from both the pyrolysis of sedimentary organic matter, but also to the ultimate thermovaporization of the high-molecular weight hydrocarbons detected between T3 (350ºC) and T4 (650ºC). The acquired parameters for the new Rock-Eval method are shown in Table 2.
1.3 New Calculated Parameters from the IFPEN Shale Play Method
The total quantity of thermovaporized hydrocarbons detected in Sh0 and Sh1 peaks, is assumed to provide an estimation of total available free and sorbed hydrocarbons occurring in unconventional shale samples, from which is deﬁned the HC Content Index (HCcont):
HCcont = Sh0 + Sh1 [Unit: mg of HC/g of initial rock]
The relative quantity of lighter molecular weight HC provides information about the quality of reservoired ﬂuids. The HC Quality Index (HQI) is calculated from Sh0 and Sh1 as follow:
A Production Index of Shale Plays (PIShale) provides a guiding information likely to deﬁne sections of interest within a shale play unit:
The calculated parameters for the new IFPEN Shale Play method are given in Table 3.
The investigated samples are derived from 5 source rocks of different origins and ages (Tab. 4), potentially being unconventional shale resource systems at different thermal maturity stages: Silurian Shale (Algeria), Mississippian Barnett Shale (Fort Worth Basin, Texas, USA), Early Jurassic Shale (Paris Basin, France), Late Jurassic Bazhenov Shale (Russia) and Eocene Green River Shale (Uinta Basin, USA). The Silurian Shale samples come from sections in southern Algeria, North Africa Platform (Ghadamis Basin).
The Silurian silty-shale succession of the Ghadamis Basin consists of the organic-rich black shales and marls considered as a candidate for shale play production.
The geological setting of this basin is described in detail by Underdown and Redfern (2008). The Barnett Shale samples are from two boreholes (Mesquite 1, Blakely 1) and an outcrop (San Saba County) (Romero-Sarmiento et al., 2014) in the Fort Worth Basin (Texas, USA). The geological setting of the Mississippian Barnett Shale play has been the object of several publications and will not be reviewed here (Montgomery et al., 2005; Ewing, 2006; Pollastro et al., 2007; Romero-Sarmiento et al., 2013, 2014). The Torcian Shale samples come from Early Jurassic intervals of the Paris basin (France).
The current petroleum potential of this basin is mainly associated to the Liassic interval where the prospective “Banc de Roc” limestone (prone for hydraulic fracking) is stratigraphically located between the Lower Domerian Shale and the Toarcian Shale (Espitalié et al., 1987; The Bazhenov Shale samples were collected from three boreholes located in the West Siberian Basin (Russia).
The Late Jurassic Bazhenov Shale is the main source rock of the largest petroleum basin in the world. More details about the depositional setting of the Bazhenov Shale are found in Peters et al. (1993) and Littke et al. (1999). Finally, a sample from the Green River Shale (Uinta basin, USA) was also investigated. This immature sample forms part of the Parachute Creek Member of the Green River Formation (Ruble et al., 2001; Behar et al., 2010).
2 Results and Discussion
Figure 3 illustrates the Rock-Eval pyrogram generated by the IFPEN Shale Play method . For instance, the Late Jurassic Bazhenov Shale sample clearly shows the deﬁnition of Sh0, Sh1 and Sh2 peaks. We suggest that the thermovaporized hydrocarbons (Sh0 and Sh1; Fig. 3) that can be detected between T1 (100ºC) and T3 (350ºC), are assigned to the part of free or sorbed hydrocarbons still present within shale sample.
Figure 3 New Rock-Eval pyrogram generated by IFPEN Shale Play method for selected sample from Late Jurassic Bazhenov Shale (Russia).
Here, the investigated samples were run using both Basic/ Bulk-Rock and Shale Play methods in order to compare the generated Rock-Eval pyrograms corresponding to each method (Fig. 4a and 4b, respectively).
Figure 4a shows that a part of the free or sorbed low-molecular weight hydrocarbons is lost because the thermovaporization stage of the Basic/Bulk-Rock method starts at a too high temperature (300ºC). In contrast, with the new IFPEN “Shale Play” method (Fig. 4b), we can observe that the selected low temperature (T1 = 100ºC) is more appropriate to optimize the recovery of most of the thermovaporizable hydrocarbons.
The “S1 peak” equivalent (S1 = Sh0 + Sh1) acquired from IFPEN Shale Play method can be compared, even if not directly, to the concept of S1 peak obtained by the standard Basic/Bulk-Rock method (Fig. 5). Most investigated samples show that “S1 peaks” equivalents from Shale Play method are substantially higher than S1 peaks obtained by the classical Basic/Bulk-Rock method (Fig. 5). Consequently results indicate an improved assessment of extractable hydrocarbon compounds in unconventional shale systems.
Figure 4 Rock-Eval pyrograms showing a) S1 and S2 peaks from Basic/ Bulk-Rock method; b) Sh0, Sh1 and Sh2 peaks from IFPEN Shale Play method: Late Jurassic Bazhenov Shale sample.
In order to investigate the type of hydrocarbons detected by each new acquired Rock-Eval parameter, samples from Bazhenov and Silurian Shale were also selected due to their signiﬁcant Sh0 and Sh1 values. Following the extraction procedure described in Behar et al. (1989), six samples were ﬁrstly extracted with n-pentane (n-C5) in order to recover, on one hand, the n-C5 extracts and, on the other hand, the n-C5 extracted solid residues. These solid residues were then extracted with DiChloroMethane (DCM) to ﬁnally recover DCM extracts and the so-called more extensively “extracted samples” corresponding, by the end, to the solid residue of two successive extractions (n-C5 and DCM).
Figure 5 Comparison of S1 peak from Basic/Bulk-Rock method and Sh0 + Sh1 peaks from IFPEN Shale Play method for investigated samples.
The n-C5 and DCM extracts as well as the extracted sample were then analyzed by the new IFPEN Rock-Eval Shale Play method . It should be noted that aliquots of n-C5 and DCM extracts were weighed (i.e. 8 to 12 mg) in Rock-Eval crucibles, previously ﬁlled by silica following procedures also described in Trabelsi et al. (1994). Table 5 provides n-C5 and DCM extract weights from investigated samples.
Today, it is widely accepted that the n-C5 extract contains saturates, aromatics and resins that are soluble in n-pentane whereas the following DCM extract contains more polar compounds, which include the total asphaltenes and remaining resins that are not soluble in n-pentane (Behar et al., 1989, 2008). The proportions of the new Rock-Eval parameters for both the n-C5 and DCM extracts are illustrated in Figures 6a and 6b, respectively.
In this study, we assume that the n-C5 extract mostly con- tains the thermovaporized hydrocarbons of low-molecular weight released between 100 and 200 C (Sh0 peak; Fig. 6a). The n-C5 extract also contains a non-negligible part of the thermovaporized hydrocarbons of high-molecular weight (Sh1 peak; Fig. 6a). Sh0 and Sh1 peaks probably correspond to low and high-molecular weight thermovaporized hydrocarbons, respectively (i.e. tentatively saturates and aromatics soluble in n-pentane).
The low Sh2 peak from the n-C5 extract probably corresponds to resins that are soluble in n-pentane (Fig. 6a). In contrast, the DCM extract mainly contains the high-molecular weight hydrocarbons detected between 200º and 650ºC. Here, the Sh2 peak from the DCM extract probably corresponds mainly to resins and asphaltenes (Fig. 6b). Finally, the extracted sample only contains the “insoluble” organic matter as shown in Figure 7, which is cracked by pyrolysis. The Sh2 peak from the extracted sample after n-C5 and DCM extractions therefore corresponds to the kerogen fraction (Fig. 7).
Figure 6: Rock-Eval pyrograms from the new IFPEN Shale Play method: a) n-C5 and b) DCM extracts: Late Jurassic Bazhenov Shale.
Furthermore, in order to compare the new Rock-Eval Tmax (TpSh2 - DTmax) calculated by the proposed new pyrolysis method with the widely accepted Tmax from Basic/ Bulk-Rock method, we plotted the average Tmax values from “extracted samples” obtained by the classical Basic/ Bulk-Rock method versus the average Tmax values from “initial samples” obtained by the new Shale Play method (Fig. 8).
Figure 7 Rock-Eval pyrograms from the new IFPEN Shale Play method comparing initial versus extracted sample after successive n-C5 and DCM extractions: Late Jurassic Bazhenov Shale.
It should be noted that the Rock-Eval Tmax from Basic/Bulk-Rock method measured on unconventional shale samples could be affected by the possible occurrence of hydrocarbons still present in the rock after thermovaporization step. Consequently, in order to avoid any Tmax bias, solvent extractions of initial samples are performed. Results indicate that Tmax parameter of “extracted samples” obtained from Basic/Bulk-Rock method is consistent with the Tmax parameter of “initial samples” obtained by the Shale Play method.
Figure 8 shows that the Tmax parameter, depending on thermal cracking of organic matter, is not affected by the temperature program dedicated to the thermovaporization of free and sorbed hydrocarbons into Sh0 and Sh1 peaks.
Figure 8 The average Rock-Eval Tmax values from “extracted samples” obtained by the classical Basic/Bulk-Rock method versus the average Tmax values from “initial samples” obtained using the new Shale Play method.
Moreover, the optimization of thermovaporization enables a Rock-Eval Tmax parameter determination more representative of original organic matter in source rock (previous to oil formation). A noticeable implication for operational purposes is that extraction by organic solvents is probably not necessary to obtain a representative and accurate Tmax parameter directly on “impregnated” samples.
The proposed method can be easily programed within any Rock-Eval 6 instrument. It provides an efﬁcient, fast and ready to use tool informing on the quantities of “free and sorbed” hydrocarbons in shale plays. Moreover, key parameters are proposed to deﬁne sections of interest and this pyrolysis program is likely to provide a meaning-ful Tmax value by avoiding “contamination” by heavy free hydrocarbons which might bias this crucial parameter in shale plays. Even if some future works are still necessary to evaluate with precision both the amount and the type of hydrocarbons detected by each new acquired Rock-Eval parameters (Sh0, Sh1 and Sh2), the present contribution shows that the new IFPEN Shale Play method may provide original and pertinent indicators of the hydrocarbon retention within unconventional source rocks such as tight, fractured and hybrid shale plays and improves the assessment of thermal maturity (Tmax) directly from non-extracted rocks.
The authors acknowledge IFP Energies nouvelles for providing approval to publish the cited patent 14/55.009, Pillot et al., 2014. We thank GASH consortium for providing Barnett Shale samples. The editor-in-chief of OGST and two anonymous reviewers are acknowledged for useful and constructive comments on the manuscript.
- Behar F., Leblond C., Saint-Paul C. (1989) Analyse quantitative des efﬂuents de pyrolyse en milieu ouvert et fermé, Revue de l’Institut Français du Pétrole 44, 387-397.
- Behar F., Beaumont V., De B., Penteado H.L. (2001) Rock-Eval 6 Technology: Performances and Developments, Oil & Gas Science and Technology 56, 111-134.
- Behar F., Lorant F., Lewan M. (2008) Role of NSO compounds dur- ing primary cracking of a Type II kerogen and a Type III lignite, Organic Geochemistry 44, 387-397.
- Behar F., Roy S., Jarvie D. (2010) Artiﬁcial maturation of a Type I kerogen in closed system: Mass balance and kinetic modelling, Organic Geochemistry 41, 1235-1247.
- Brigaud B., Durlet C., Deconinck J.-F., Vincent B., Pucéat E., Thierry J., Trouiller A. (2009) Facies and climate/environmental changes recorded on a carbonate ramp: a sedimentological and geo- chemical approach on Middle Jurassic carbonates (Paris Basin, France), Sediment. Geol. 222, 181-206.
- Carpentier B., Huc A.Y., Hamou P., Wilhems A. (1995) Detection, distribution and origin of thin tar mats in the Miller ﬁeld (North Sea, UK), 17th International Meeting on Organic Geochemistry, San Sebastian, Spain, pp. 388-390.
- Carpentier B., Huc A.-Y., Marquis F., Badr A.E.R., Al Aidarous A.A., Al-Baker S. (1998) Distribution and origin of a Tar Mat in the S. Field (Abu Dhabi, A.E.U.), SPE Paper 49472, 1-10.
- Espitalié J., Deroo G., Marquis F. (1986) La pyrolyse Rock-Eval et ses applications, Revue de l’Institut Français du Pétrole 41, 73-89.
- Espitalié J., Marquis F., Sage L. (1987) Organic geochemistry of the Paris Basin, Brooks J., Glennie K. (eds), Petroleum Geology of North West Europe, Graham and Totman, London, pp. 71-86.
- Ewing T.E. (2006) Mississippian Barnett Shale, Fort Worth Basin: North-central Texas: gas-shale play with multi-trillon cubic foot potential. Discussion, American Association of Petroleum Geolo- gists Bulletin 90, 963-966.
- Lafargue E., Marquis F., Pillot D. (1998) Rock-Eval 6 applications in hydrocarbon exploration, production, and soil contamination studies, Oil & Gas Science and Technology 53, 421-437.
- Littke R., Cramer B., Gerling P., Lopatin N.V., Poelchau H.S., Schaefer R.G., Welte D.H. (1999) Gas generation and accumulation in the West Siberian basin, American Association of Petroleum Geologists Bulletin 83, 1642-1665.
- Montgomery S.L., Jarvie D.M., Bowker K.A., Pollastro R.M. (2005) Mississippian Barnett Shale, Fort Worth Basin, north-cen- tral Texas: gas-shale play with multitrillion cubic foot potential, American Association of Petroleum Geologists Bulletin 89, 155-175.
- Perrodon A., Zabek J. (1990) Paris Basin, in Interior Cratonic Basins, American Association of Petroleum Geologists Bulletin 51, 633-679.
- Peters K.E., Kontorovich A.E., Moldowan J.M., Andrusevich V.E., Huizinga B.J., Demaison G.J., Stasova O.F. (1993) Geochemistry of selected oils and rocks from the central portion of the West Siberian basin, Russia, American Association of Petroleum Geologists Bulletin 77, 863-887.
- Pillot D., Letort G., Romero-Sarmiento M.F., Lamoureux-Var V., Beaumont V., Garcia B. (2014) Procédé pour l’évaluation d’au moins une caractéristique pétrolière d’un échantillon de roche, Patent 14/55.009.
- Pollastro R.M., Jarvie D.M., Hill R.J., Adams C. (2007) Geologic framework of the Mississippian Barnett Shale, Barnett-Paleozoic total petroleum system, Bend Arch-Fort Worth Basin, Texas, American Association of Petroleum Geologists Bulletin 91, 405-436.
- Romero-Sarmiento M.-F., Ducros M., Carpentier B., Lorant F., Cacas M.-C., Pegaz-Fiornet S., Wolf S., Rohais S., Moretti I. (2013) Quantitative evaluation of TOC, organic porosity and gas retention distribution in a gas shale play using petroleum system modeling: Application to the Mississippian Barnett Shale, Marine and Petro- leum Geology 45, 315-330.
- Romero-Sarmiento M.-F., Rouzaud J.-N., Bernard S., Deldicque D., Thomas M., Littke R. (2014) Evolution of Barnett Shale organic carbon structure and nanostructure with increasing maturation, Organic Geochemistry 71, 7-16.
- Ruble T.E., Lewan M.D., Philp R.P. (2001) New insights on the Green River petroleum system in the Uinta Basin from hydrous pyrolysis experiments, American Association of Petroleum Geolo- gists Bulletin 85, 1333-1371.
- Trabelsi K., Espitalié J., Huc A.-Y. (1994) Characterization of Extra Heavy Oils and Tar Deposits by modiﬁed Pyrolysis Methods, Proceedings of the “Heavy Oil Technologies in a Wider Europe” Thermie EC Symposium, Berlin, pp. 30-40.
- Underdown R., Redfern J. (2008) Petroleum generation and migration in the Ghadamis Basin, North Africa: a two-dimensional basin- modeling study, American Association of Petroleum Geologists Bulletin 92, 53-76.
Cite this article as: M.-F. Romero-Sarmiento, D. Pillot, G. Letort, V. Lamoureux-Var, V. Beaumont, A.-Y. Huc and B. Garcia (2016). New Rock-Eval Method for Characterization of Unconventional Shale Resource Systems, Oil Gas Sci. Technol 71, 37.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.