Thermal maturity evaluation using Raman spectroscopy for oil shale samples of USA: comparisons with vitrinite reflectance and pyrolysis methods

Thermal maturity is commonly assessed by various geochemical screening methods (e.g., pyrolysis and organic petrology). In this contribution, we attempt to establish an alternative approach to estimating thermal maturity with Raman spectroscopy, using 24 North American oil shale samples with thermal maturity data generated by vitrinite reflectance (VRo%) and pyrolysis (Tmax)-based maturity calculation (VRe%). The representative shale samples are from the Haynesville (East Texas), Woodford (West Texas), Eagle Ford and Pearsall (South Texas) Formations, as well as Gothic, Mancos, and Niobrara Formation shales (all from Colorado). The Raman spectra of disordered carbonaceous matter (D1 and G bands separation) of these samples were directly obtained from the rock chips without prior sample preparation. Using the Gaussian and Lorentzian distribution approach, thermal maturities from VR were correlated with carbon G and D1. We found that the Raman band separation (RBS) displayed a better correlation for equivalent VRe% than vitrinite reflectance VRo%. The RBS (D1–G) distance versus total organic carbon, free hydrocarbons from thermal extraction (S1), and the remaining hydrocarbon generating potential (S2) indicate that the RBS (D1–G) distance is also related to kerogen type. Data presented here from three methods of maturity determination of shale demonstrate that Raman spectroscopy is a quick and valid approach to thermal maturity assessment.


Introduction
Geochemical evaluation is a critical aspect in evaluating source rock reservoir plays, specifically assessing the quantity, quality, and thermal maturity of sedimentary organic matter, and the implication for hydrocarbon generation, retention, and expulsion (Carvajal-Ortiz and Gentzis 2015). The thermal maturity assessment of source rock-reservoir rocks is routinely performed using pyrolysis (T max ) and vitrinite reflectance (VR o %) data, though there are some other less common methods available, such as spore coloration index, solid bitumen reflectance, and spectral microfluorescence (e.g., Tissot 1984;Taylor et al. 1998).
Vitrinite is the remains of woody materials, commonly found in source rock-reservoir rocks, such as shales. Vitrinite reflectance (VR o %) values increase with maturity (generally a function of depth) due to progressive changes in aromatization, condensation, and rearrangement of the vitrinite molecular structure (McCartney and Teichmüller 1972;Levine and Davis 1984;Carr 1999). Although this method is effective and widely used, it is expensive, labor-intensive, and dependent on an organic petrographer's judgment and experience that can significantly vary from person to person. Hackley and Cardott (2016) showed that the identification of vitrinite macerals versus certain other macerals, such as various types of solid bitumens, inertinite, and semi-fusinite, is challenging. Moreover, large variations (20%-30%) in vitrinite reflectance of the same sample determined by different laboratories are a common problem.

3
Another widely used maturity assessment method, pyrolysis, involves thermal decomposition of organic matter in an inert atmosphere (Peters 1986;Lafargue et al. 1998), yielding a temperature at the maximum rate of hydrocarbon generation (T max ) from pyrolysis of the S2 peak, along with other parameters, such as S1, S2, and S3. Additionally, a broad S2 peak in a pyrogram often compromises the accuracy and precision of the T max value. Peters (1986) showed that when S 2 ≤ 1 mg HC/g rock or the formation has a low TOC, the measured T max value could be unreliable.
Raman spectroscopy is potentially a better technique for nondestructive, rapid, and precise maturity evaluation of hydrocarbon source rocks (Kelemen and Fang 2001;Marshall et al. 2012;Wilkins et al. 2014). In recent years, Raman spectroscopy has emerged as a indispensable useful method that provides information about the structural and chemical properties of carbon materials (e.g., Ferrari and Robertson 2000;Beyssac et al. 2002;Xie et al. 2008;Wilkins et al. 2014Wilkins et al. , 2015Zhou et al. 2014;Schito et al. 2017;Sauerer et al. 2017;Schmidt et al. 2017;Khatibi et al. 2018;An et al. 2018;Hackley and Lünsdorf 2018). This method can also be used to study small, dispersed maceral grains that are difficult to study by other methods. This technique relies on peaks of carbon being resolved into Gaussian bands (Wilkins et al. 2014). The Raman spectrum of the organic carbon materials shows two primary bands: G band (graphitic) near ~ 1580 cm −1 and D band (disordered) around ~ 1350 cm −1 (Tuinstra and Koenig 1970;Beny-Bassez and Rouzaud 1985;Gruber et al. 1994;Spötl et al. 1998;Kelemen and Fang 2001;Marshall et al. 2012;Schito et al. 2017). However, with an increase in thermal maturity, the D and G peaks shift toward lower and higher wavelengths, respectively. The distance between G and D peaks is called Raman band separation (RBS), and it has been used as a reliable maturity indicator (Kelemen and Fang 2001;Marshall et al. 2012;Wilkins et al. 2014;Sauerer et al. 2017;Schito et al. 2017).
Another important issue in recent works is that only the vitrinite reflectance is considered as the primary benchmark for correlation with Raman spectral data. The main problem is that the vitrinite reflectance itself shows a 10%-30% error, as the measurement is based on personal judgment. Recent studies such as Sauerer et al. (2017) conducted on 11 samples that are all type II organic matter, while Kelemen and Fang (2001) worked on shale samples of type II and Spötl et al. (1998) on type III organic matter (OM), respectively. The studies of Wilkins et al. (2015) and Schito et al. (2016) are based on core samples at different depths from a single well and were correlated with Raman spectroscopy, while Khatibi et al. (2018) concentrated only on Williston Basin, North Dakota. Our work examined Raman spectroscopybased maturity assessment on samples from several different basins across the USA (Table 1). Raman spectral data were compared with both vitrinite reflectance-based maturity and pyrolysis-derived equivalent reflectance maturity of the same samples. In contrast with other studies, we examine here the correlation of different Raman spectral parameters data with maximum, minimum, and average vitrinite reflectance thermal maturity and equivalent vitrinite reflectance (derived from pyrolysis experiments T max value). Another unique aspect of our present study is that it focuses on the various types of kerogen materials-type II-III, type III, and type IV.

Samples
Sixteen samples with known VR o % and 21 samples with VR e % were selected from multiple organic-rich mudstones, shale and chalk formations in active USA source rock-reservoir plays (Table 1). The age of these samples ranges from Late Devonian to Cretaceous. Samples were obtained as 10-20 g splits from mostly drill cores (sidewall plugs), other than Eagle Ford samples which were from outcrops. Cores were chipped into approximately 10-15 mm size cube pieces for the Raman spectroscopy analyses. The powdered shale sample fractions (< 75 μm) were used to determine TOC using both Shimadzu TOC-VWS wet chemical total organic carbon analyzer at the University of Texas at Arlington, USA, and a Leco CS-200 carbon/sulfur analyzer in Nanjing University, China. The pyrolysis (< 75 μm) and vitrinite reflectance data (core chips) were obtained from external laboratories, as described below.

Pyrolysis and organic richness
Pyrolysis was performed on an OCE-II oil-gas evaluation workstation following Chinese Technical Standard GB/T 18602-2012(SAC 2012 at the China University of Geosciences in Wuhan. Previous work has explained in detail about different pyrolysis apparatuses-their operations, measured parameters, potential pitfalls, and applications to petroleum exploration (Behar and Pelet 1985;Espitalié 1986;Peters and Cassa 1994;Lafargue et al. 1998;Cheshire et al. 2017). The OCE-II oil-gas evaluation workstation performs temperature-programmed heating to (1) thermally distil "free" or adsorbed hydrocarbons (S1 peak) from the sample, (2) pyrolyze the kerogen to produce hydrocarbons (S2 peak) in a second oven, and (3) obtain residual organic carbon content in mg/g (S4 unit mass) of hydrocarbon source rock after pyrolysis. This instrument is not equipped with the S3 measurement unit which can detect unit mass of released carbon dioxide (mg/g) from thermal breakdown of kerogen.
The T max is estimated based on S2 peak temperature that coincides with the maximum generation of hydrocarbons. The total hydrocarbon potential of the sample is defined by (S1 + S2), the hydrogen index (HI) = S2/TOC, and production index (PI) = S1/(S1 + S2). Other parameters obtainable from pyrolysis are the bitumen index (BI), defined as S1/TOC (Killops et al. 1998), and the quality index (QI) ~ (S1 + S2)/TOC (Pepper and Corvi 1995). The equivalent vitrinite reflectance is deduced from T max (°C) using the following equation (Jarvie et al. 2001;Jarvie, 2012): VR e % = [0.018 × T max ] − 7.16. A calculated estimation of the uncertainty in T max based upon duplicate analyses of identical sample splits is 3-5 °C for samples with S2 values > 10 mg HC/g rock, equivalent to an internal precision of ± 0.05-0.10 for the estimated VR e %.

Vitrinite reflectance measurements
Vitrinite reflectance measurements were conducted at two places. The analyses performed at Nanjing University followed the procedure of SY/T 5124-2012 (2012)

Raman spectroscopy measurements
Raman spectroscopic measurements were performed on the flat surface of rock chips with minimum sample preparation (such as level the surface to fit into the stage on Raman spectroscopy), with the aim to establish a fast and convenient thermal maturation determination method with minimal laboratory preparation and analysis. All data in this study were obtained with a Thermo Scientific DXRxi Raman imaging spectroscopy. Spectra were generated mainly using 10 × and sometimes using 50 × objectives, by monochromatic excitation with 532 nm wavelength laser dispersed by a fixed 900 lines/mm holographic grating. The laser spot size for analyses was 10 μm using 10 × objective and 2 μm while using 50 × objectives. The operating conditions were 1-5 s for spectral acquisition, 1-10 accumulations, and 0.7-5 mW for laser power. Though the DXRxi Raman imaging microscope is calibrated through auto-alignment, Si (520.5 cm −1 peaks) is used to check the calibration. The Raman spectra were obtained and processed using the OMNICxi Raman Imaging software. Sometimes, carbonaceous matter does not produce useful Raman spectra because of the "fluorescence" effect, from the presence of certain organic matter types. The fluorescence effect depends upon the interaction between the laser and the organic materials being examined and is sensitive to acquisition time and laser power. To minimize this effect, statistically large amounts of data were generated from each sample under different operating conditions (e.g., scan time, laser power, and magnification) so that all the RBS data given in this manuscript are an average of 20 spectra collected from different positions of the same rock chip. Here, we did not differentiate maceral type for collecting different spectra.

Curve fitting of Raman spectroscopy
Raman spectra of highly disordered carbonaceous materials are affected by interference due to fluorescence caused by the presence of long-chain alkenes and aromatic compounds. The first step in the analysis of the spectra is the removal of the high fluorescence background by a baseline subtraction procedure. The baseline subtraction was performed using a third-order polynomial curve that best represents the real trend of the fluorescence in the most immature samples with high background. To avoid errors induced by differences in the spectral range considered by the subtraction, we fixed baseline points at 1100 and 1700 cm −1 for all spectra. Once the process was optimized, it must be then reproduced for the entire set of samples. Deconvolution and fitting of the Raman spectra were made using five peaks with mixed Gaussian/Lorentzian profiles and linear baseline correction (Origin 2017 software). No constraints were put on the peak parameters in the curvefitting procedure to allow an optimal fit to the measured spectra. We observed that the maximum difference resulting from baseline corrections was detected for the D band intensity, whose values tend to be higher by 10%-15%. After removal of the background, the spectra were deconvoluted using Origin software. As the Raman spectra of organic matter (OM) in diagenetic conditions show broad overlapping bands, and from the application of multiple curve-fitting solutions (Lünsdorf and Lünsdorf 2016;Lupoi et al. 2017Lupoi et al. , 2018Jubb et al. 2018;Hackley and Lünsdorf 2018), we found that the best solution using the minimum number of components could be attained using a five-band deconvolution with a mixed Gaussian-Lorentzian band profile (Fig. 4). All spectra were split using the five-peak deconvolution procedure. The band positions, amplitudes (intensities), full width at half maximum (FWHM), and integrated areas of all peaks were extracted from the resulting curve fits. Although the 20 spectra for each sample showed significant ranges in the peak widths, intensities, and areas, the mean G and D1 band positions and the error estimations for G and D1 positions were determined with a 1σ error ranges between 0 and 5 cm −1 . Table 2 shows TOC, S1, S2, kerogen type, T max values, and calculated values, such as HI, BI, QI, and PI, from pyrolysis of 21 samples. The TOC value ranges between 0.74 and 7.8 wt%, S1 values are between 0.18 and 4.24 mg HC/g rock, and S2 values range around 0.18-10 mg HC/g rock (Table 2). Hydrogen index (HI) of these samples is between 12 and 891 mg HC/g TOC. The HI values are lower than 100 mg HC/g TOC, suggesting that most of the samples are thermally mature to postmature.

Pyrolysis and equivalent reflectance
Note that five Niobrara samples yield high S1, S2, HI and (S1/TOC) × 100, and respective T max values range between 419 and 440 °C. Vitrinite reflectance equivalence from T max is estimated as 0.40 to 0.78VR e % (Table 1). Calculated PI values are within the range of 0.23-0.69. Based on the pyrolysis data, one sample had type I kerogen, three samples with type II-III kerogen, nine samples with type III kerogen and eight samples with type IV kerogen. The T max versus HI cross-plot of pyrolysis data showed that half the samples are II-III and III type and thermally mature to the postmature window which represents an oilto-gas generation condition. Three samples are immature, and the rest of the samples show an inert state (Fig. 1). A graph of S2 versus TOC indicates a direct correlation as the S2 value increases with the TOC value (Fig. 2a). The T max versus bitumen index (BI) and quality index (QI) shows that with an increase in T max value, the value of BI and QI decreases (Fig. 2b, c). A cross-plot between TOC and QI, where QI is combining S1 and S2, shows slightly sifted no correlation between the amount of then TOC and S2QI (Fig. 2d).

Raman spectroscopy and thermal maturity
The described fitting procedure was carried out on all sample spectra. An example of fitting spectra (Sample 5: TX P-3987) is shown in Fig. 4 about the position and band names of carbonaceous material on Raman spectra. Best-fit results were achieved by considering G and Dl as purely Gaussian and D2, D3, and D4 as mixed Gaussian-Lorentzian functions. It is worth mentioning here that RBS is found as a more robust parameter than G and D1 band area ratio while comparing with equivalent reflectance (%VR e ). The efficacy of RBS with thermal maturity was also observed by Sauerer et al. (2017). The RBS versus vitrinite reflectance (%VR o ) shows a poor correlation suggesting a higher error in petrographic measurements.

Correlation of Raman band separation (RBS) with maturity
The D1 and G Raman bands represent processes occurring at different carbonization ranks. The D1 band shift toward lower wavelength can be assigned to the growth of larger aromatic clusters passing from disordered to ordered organic material (Ferrari and Robertson 2000  organic shales in different basins shows an increase in RBS with sample depth (R 2 = 0.38) (Fig. 5a), and this increment in RBS suggests the higher degree of thermal maturity. The RBS distance increases with depth from about 225 to 270 cm −1 . When correlating all the Raman data with both vitrinite reflectance VR o % and equivalent vitrinite reflectance VR e %, the higher shift of G position and lower shift of D1 position with vitrinite reflectance VR o % shows QI, S1+S2 TOC Fig. 2 Correlation between the measured parameters by pyrolysis: a S2 versus TOC; b T max versus BI (S1/TOC); c T max versus QI (S1 + S2); d TOC versus QI (S1 + S2)  Fig. 4 Raman spectral characteristics of common five-band deconvolution outcome for an organic matter spectrum; D = disordered bands (1, 2, 3, 4) and G = graphitic band a weak relationship (Fig. 5c, e). The D1 band position shifts toward lower values with the increase in VR e % from about 1370 to 1330 cm −1 (Fig. 5d). Similarly, the G band position shows a shift from 1592 to 1608 cm −1 with an increase in thermal maturity (Fig. 5f). The data show a modest correlation for an increase in the RBS distance with equivalent vitrinite reflectance determined by the shift of D1 position toward lower wavenumber (Fig. 5d, R 2 = 0.70) and G position toward higher wavenumber position (Fig. 5f, R 2 = 0.83). This observation is similar for RBS distance with maturity in other work (Ferrari and Robertson 2000;Kelemen and Fang 2001;Guedes et al. 2010;Liu et al. 2012;Lünsdorf and Lünsdorf 2016). Our data clearly show that vitrinite reflectance-based maturity does not show a good correlation with RBS, as emphasized by previous researchers (e.g., Khatibi et al. 2018). Tuinstra and Koenig (1970) reported that a vertical intensity ratio between G and D1 bands (ID1/IG) is related to the in-plane crystallized size in disordered graphite and gave a weak proportional correlation against T max in a range between 420 and 500 °C. We observe that the intensity ID1/ IG ratio also has a weak correlation with the T max value. Some researchers showed a significant correlation between G and D1 band intensity versus thermal maturity (Wilkins et al. 2014(Wilkins et al. , 2015Schito et al. 2017;Sauerer et al. 2017;Schmidt et al. 2017). However, this study does not show any significant correlation (Fig. 6). The Raman band intensity largely depends on the analytical condition, whereas the wavenumber position is independent of the analytical condition and better represents the structural ordered/disordered mode of matter (Matthews et al. 1999).

Correlation of Raman band area and band separation against thermal maturity
In this section, we use the distance and full width at half maximum (FWHM) of the D and G bands area to correlate with thermal maturity. Here, we use mainly two area ratios based on Raman spectra: 1. area ratio 1 = D1 area/G area 2. area ratio 2 = D1 area/(D1 + D2 + G) area These parameters are plotted against the maturity profile versus VR o % and VR e % (Fig. 6). VR o % versus area ratio 1 (D1 area/G area) shows a moderate correlation (R 2 = 0.49). In Fig. 6a, an increase in area1 ratio indicates an increase in thermal maturity. VR e % versus area1 also shows a moderate correlation (R 2 = 0.48) (Fig. 6b). VR o % versus area ratio 2 [D1 area/(D1 + D2 + G) area] also shows a moderate correlation (R 2 = 0.46) representing an increase in area ratio 2 proportionally related to the increase in VR o % (Fig. 6c). VR e % versus area ratio 2 shows a very poor correlation (Fig. 6d).
From Fig. 5c-f, it is observed that D1 wavelengths at peak maxima move to lower numbers and G peak maxima move to higher wavelength maxima with progressively higher thermal maturity. The correlation between RBS distance and VR o % and VR e % demonstrates that the VR e % data are more comparable with the Raman maturity method. The plot of VR o % versus RBS distance shows a weak correlation (Fig. 7a). This suggests that the RBS distance and VR o % are only poorly proportional because the increase in RBS distance corresponds only very generally to higher VR o %. In comparison, VR e % versus RBS distance plot displays a good nonlinear correlation (R 2 = 0.76) (Fig. 7b). The graph shows that the equivalent VR e % maturity increases with higher RBS distances. This study indicates that based on the relationship between equivalent reflectance VR e % and RBS distance, thermal maturity can be estimated by direct Raman spectroscopic analyses. Pyrolysis-based T max value versus RBS distance represents a fair linear correlation (R 2 = 0.73), with an increase between D1-G distance and T max .

Correlation of RBS, pyrolysis, and kerogen type
The relationship between RBS distance and pyrolysis data shows that RBS distance versus S1 (free remaining hydrocarbon) plot does not exhibit any direct relation. However, in general, it shows that small S1 values have large RBS distances. RBS distance versus S2 (remaining generating hydrocarbon potential) graph shows a similar trend, with a smaller S2 value indicating higher RBS distance. Figure 8 shows the relationship of thermal maturity and RBS distance of different kerogen types. Vitrinite reflectance VR o % versus kerogen-type RBS distance graph represents that there is a weak correlation between kerogen type and vitrinite reflectance (Fig. 8a). A plot of equivalent vitrinite reflectance VR e % versus kerogen-type RBS distance represents an excellent linear correlation for some kerogen type. Kerogen type II-III signifies an excellent relationship with RBS distance (R 2 = 0.98) (Fig. 8b). Kerogen type III and IV also show a good relationship with RBS distance. Equivalent vitrinite reflectance VR e % maturity increases with higher RBS distance (R 2 = 0.84) (Fig. 8b).  VRo%   I  II-III  III  IV  I  II-III  III

Pyrolysis and vitrinite reflectance
The samples used for the present study show a marked difference in thermal maturity between vitrinite reflectance and equivalent reflectance from T max . A correlation of thermal maturity with Raman band separation indicates that an equivalent VR e % is more appropriate for correlation purposes. It is essential to consider both T max and vitrinite reflectance VR o % values within a geological context for a proper understanding of maturity due to complex relations between paleoenvironment of deposition, organic matter types, and diagenetic development of the sedimentary formation (Snowdon 1995;Sykes and Snowdon 2002;Dembicki 2009;Carvajal-Ortiz and Gentzis 2015). However, Fig. 9 (vitrinite reflectance VR o % vs. HI relationship) and Fig. 1  Several researchers have shown that application of vitrinite reflectance to the determination of the thermal maturity of marine and lacustrine shales introduces frequent sources of error (Dennis 1974;Carvajal-Ortiz and Gentzis 2015). One of the critical issues is that the vitrinite reflectance scale is based on vitrinite, which is abundant and well preserved in coal. Due to the sedimentary reworking process of shale, the properties and volume of carbonaceous matter are different from coal. For this, the thermal maturity assessment  1 3 method of shale should be different with respect to coal. Barker (1996) suggested a minimum of 20 maceral measurements for better statistical representation of constituent macerals. However, only 20 measurements can limit the validity of the calculated average value. For example, our sample 19 (Fig. 3) has 20 maceral observations. Meanwhile, there are no adequate data to outline a normalized distribution. Hackley and Cardott (2016) described that the most efficient method of the vitrinite reflectance estimation is to look not only at the mean value and the number of measurements, but also to the shape of the reflectance histogram, the spread of values, and photomicrographs of identified vitrinite.

Raman band parameters, area ratios, and their application to dispersed organic matter
With the aid of Raman spectroscopy, D1 and G band positions are controlled by the development of different types of carbonaceous material within the shale samples. The longchain aromatic clusters are transformed from disordered to ordered clusters when D1 bands shift toward lower wavenumbers (Fig. 5b;Ferrari and Robertson 2000). The change in Raman shift of D1 position toward a lower wavenumber (Fig. 5b) and G position toward higher wavenumber values controls the RBS distance (Fig. 5d), which clearly shows a direct relation with thermal maturity (Fig. 7b). This observation supports previous work on Raman band shifts (Ferrari and Robertson 2000; Kelemen and Fang 2001;Guedes et al. 2010;Liu et al. 2012;Schito et al. 2017). In this study, Raman spectral data fit well with VR e % (i.e., T max ), but not VR o %. Figure 7b demonstrates that the RBS distance is a good parameter to use for determining thermal maturity. In this work, the G band is fitted around 1603 cm −1 (Fig. 4) surrounded by a Lorentzian peak at 1610 cm −1 and a Gaussian D2 peak around 1500 cm −1 . Heise et al. (2009) observed that the D2 band corresponds to the same vibration as the G band in small-size aromatic hydrocarbons. The D band area is decoupled to a central Lorentzian D band at 1340.45 cm −1 , and two outer area Gaussian bands on D4, to the left of D, and D2 to the right of D, at 1250 and 1500 cm −1 , respectively, that were assigned to identical vibrations of the D band but in small-size aromatic domains.
In this study, we have not observed a strong correlation between Raman intensity and maturity, but rather a weak correlation between D1/G area ratio and VR e %. Moreover, several spot analyses at different analytical conditions from a single sample show a wide variation in Raman spectral intensity. The low correlation coefficient in band intensity versus T max suggests a lack of definite relationship, while Kelemen and Fang (2001) and Schito et al. (2017) found the most significant maturity relation with the G and D1 band intensity. Previous workers also noticed that the D2 and D3 band areas decrease with the increase in thermal maturity in shale (Li 2007). The area ratio of D1 and G bands and parameters related to this ratio increase with the increase in thermal maturity (Fig. 6). Our results and earlier research show that an increase in structural ordering increases the Raman shift value of D1 band, and the Raman shift wavelength of G, D2, and D3 band decreases (Ferrari and Robertson 2000;Li 2007;Rebelo et al. 2016). Rebelo et al. (2016) showed that D4 and D bands have a similar pattern of vibrational modes and structural units for D2 and D3 bands for nanotube. Schito et al. (2017) observed similar coupled behavior for these bands at increasing temperatures. However, our study does not indicate any maturation trend for this complex band (D3 and D4). This result could be due to the different saturated and unsaturated hydrocarbon components in OM of the shale samples, or a variation of aromatic OM enrichment which results in different responses of the excitation bands with increasing thermal maturity. At experimental lower temperature (about 200-320 °C), carbonaceous material displays band D4 from sp 3 to sp 2 shells around 1250 cm −1 (Lahfid et al. 2010), which is also observed in this study (Fig. 4). Lis et al. (2005) suggested that the heterogeneities in C-C and C-H bonds in aliphatic compounds at low maturity stages (VR o % < 1.5) are complex. In general, the changes in the Raman spectra are possibly due to a reduction in branched aliphatic chains, shortening of aliphatic chains, and formation of cyclic/aromatic compounds. Beyssac et al. (2002) observed the presence of two additional bands in the Raman spectrum of organic matter. One of the bands appears as a right end shoulder on the G band at approximately 1620 cm −1 if the sample is over-mature. In the present study, we do not have any postmature samples; the D2 band is a wideband existing at around 1500 cm −1 .

Thermal maturity and kerogen type correlation from Raman spectroscopy
Estimated vitrinite reflectance equivalences of the samples were correlated with the Raman band separation (G band shift minus D1 band shift), as described in the earlier sections, and were found to correlate with thermal maturity (e.g., Kelemen and Fang 2001;Schito et al. 2016;Sauerer et al. 2017;Cheshire et al. 2017). In this study, thermal maturity is determined using the equivalent vitrinite reflectance (VR e %) versus RBS correlation from Raman measurements (Fig. 7b). Based on Fig. 7b correlation graph unknown samples RBS can convert to thermal maturity in the range of 0.5%-1.8%VR e . This correlation is based on type II-III, type III, and type IV kerogen samples. It is worth mentioning here that the type IV is not important for hydrocarbon generation. Types II-III represent a nearly perfect correlation for thermal maturity (Fig. 8b) with RBS (R 2 = 0.98), while type III showed a good correlation (R 2 = 0.84; Fig. 8b). This correlation can also be helpful to categorize kerogen types based on their RBS value trends for unknown samples. The most significant result of this study is that equivalent VR e % has a much stronger correlation than vitrinite reflectance with thermal maturity, accompanying RBS increase with the increase in thermal maturity of the samples. The higher the maturity, the larger the RBS distance. Another important outcome of our study is that the Raman spectral intensity is not a useful parameter to estimate the thermal maturity.

Conclusion
Raman spectroscopy can be used as a method for the assessment of the thermal maturity of shale samples. The present study demonstrates that with an increase in thermal maturity, D1 and G peaks shift toward lower and higher wavelengths, respectively, which increases the Raman band separation (RBS). A total of 24 shale samples from different petroleum basins in the USA with a wide range of maturity (0.6-2 VR o %) were used to document that RBS has a good correlation (R 2 = 0.76) with equivalent thermal maturity (VR e %). Therefore, we can potentially gauge the thermal maturity of shales by this method of curve fitting for unknown samples using laser Raman spectroscopy. All these correlations show that with the aid of Raman spectroscopy alone, the maturation windows for immature kerogen, oil, wet gas, and dry gas can be evaluated for unknown samples. The present study suggests that the Raman-based thermal maturity can be better correlated with equivalent VR e % instead of vitrinite reflectance VR o %. The conventional vitrinite reflectance (VR o %) determination by petrography might produce a huge error. For this, we claim that the Raman spectroscopy-based maturity estimation is a better alternative. We also claim here that the RBS method can be helpful to determine the T max . Additionally, Raman spectra can be generated directly on drill cores and cutting chips without any sample preparation. Further data mining of Raman spectra collection from different basin samples will be needed to refine and establish a fast, less strenuous, and quantitative method for measuring thermal maturity for unconventional petroleum prospects.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.