Abstract
Recent exploration work in the Tazhong district has gradually transferred to the exploitation of high and over mature oils in deep and ultra-deep layers. This has proved problematic, however, as the distribution of crude oils in the Tazhong is complex. This means that the geochemical characterization of high and over mature oils, especially for light crude oils, have become increasingly important. The stability of concerted ring structure of aromatics makes them having stronger thermal stability and resistance to biodegradation. This means that there are abundant aromatic compounds in high and over mature oils. This study presents a series of geochemical analyses of the maturity parameters of 89 crude oils from the Tazhong area, including stable carbon and hydrogen isotope analyses of compounds from 43 light crude oils. These analyses are then compared with other data from the Tazhong Number I fault zone, as well as the Tazhong Number 10 and Tazhong Uplift structural zones. Results show that the geochemical parameters of oils from Tazhong Number I fault zone generally encompass a wider range than those from the Tazhong Number 10 structural zone, which indicates that the Tazhong Number I slope belt is more active than its counterpart structural belt and generates oils with more complex geochemical characteristics. The positive correlation between the toluene/methyl cyclohexane ratio and the dibenzothiophene/phenanthrene ratio, as well as with the naphthalene/phenanthrene ratio indicates that aromatization parameters can be used to evaluate the maturity of light crude oils, and there may be inherited relationships between toluene and methyl cyclohexane in crude oils.
Introduction
The Tarim Basin has experienced a number of structural movements, which has led to a series of complex and diverse structural styles, and consequently large differences in oil characteristics between different zones. The fact that large volumes of light crude oil have been extracted from the Tazhong (TZ) area within the Tarim Basin highlights the great potential for massive petroleum exploration from deep layers in this region (He et al., 2016; Wang et al., 2014). Zhu et al. (2015) reported the newly discovered Rewapu oilfield in Tarim Basin, which was well preserved due to low geothermal stress and quick deep burial. Wang et al. (2015) analyzed sulfur-containing compounds in Ha9 and ZS1C crude oils, identifying DBT (dibenzothiophene) series compounds in Ha9 crude oil and thiadiamondoids with 1–3 cages in ZS1C crude oil. Zhang et al. (2015a) conducted a series of geochemical analysis on oil reservoir fluids, suggesting that thermal cracking and TSR have contributed a lot to the deep formation oil. However, because of the high maturity of light crude oils, a number of normal geochemical evaluation parameters cannot be applied. In addition, several structural movements have changed, or even destroyed, crude oils within the Tarim Basin to different degrees; the geothermal changes induced by these movements have affected the thermal evolution of hydrocarbon source rocks and caused large differences in the maturity of crude oils at different phases.
Two sets of source rocks from the CambrianLower Ordovician and Middleto-Upper Ordovician mainly contributed to the oils from Tarim Basin (Zhang and Huang, 2005; Hu et al., 2016; Huang et al., 2017), but the problem which one is the main source have been debated for a few of years (Huang et al., 2016). While some workers believe the main source of hydrocarbons is Middle-to-Upper Ordovician rocks (Hanson et al., 2000; Zhang et al., 2000; Zhao et al., 2005), others have argued that the Cambrian-Lower Ordovician is the main source (Liang et al., 2000; Sun et al., 2003). As this debate remains unresolved, the main aim of this work is to provide a geochemical characterization and discussion of the usage of aromatic compounds in the TZ area rather than to discuss the main source of hydrocarbons within the Tarim Basin.
Compared with other molecules, aromatic hydrocarbons with higher thermal stability are an important research object as their highly delocalized electronic conjugated system makes their energy system low (Xing et al., 2005). This is important because deep layer source rocks in the Tarim Basin are in high mature; source rocks from Cambrian and lower Ordovician strata are considered over mature, and most of these rocks from the middle-to-upper Ordovician fall within the peak oil generation and later oil windows. One good example is the TC 1 well; the equivalent vitrinite reflectance of this well between 5713m and 7124m depth within Ordovician layers is between 1.29% and 2.25% (Zhang et al., 2004), which means that most normal geochemical parameters are invalid because of high maturity (Liang and Chen, 2005). Therefore, this work also discusses the rules governing changes in aromatic hydrocarbons within TZ crude oils.
Light hydrocarbon is another series of compounds for petroleum research. Liu et al. (2015) reported the light composition, hydrogen and carbon isotope of the upper Paleozoic gases in the Ordos Basin and suggested that they were all sourced from Carboniferous–Permian coal-bed source rocks, distinguishing two type of natural gas. Han et al. (2017) used MCH (methyl cyclohexane) index to identify Upper Paleozoic gas source rocks in Shenmu gas field and think they are mainly of humic genesis of swamp facies.
Samples and methods
Samples
A total of 89 oils samples were selected for analysis in this work from the TZ Number I fault zone, the TZ Number 10 structural zone, and the central uplift zone. These samples comprise 49 black oils and 40 light crude oils, collected from different districts of TZ area. The black oils are mainly from depths between 3000 m and 5000 m, while light crude oils mainly from depths between 5000 m and 6500 m. The distribution of oil samples is shown in Figure 1.
Tectonic units that comprise the Tazhong Uplift and locations of the oil samples discussed in this work.
Methods
The samples were treated by following steps.
About 3 mg of the TZ black oils were taken and separated into saturates and aromatics by successively using hexane and a mixture solvents (hexane: dichloromethane = 1:3) in a 9 cm column filled with 3:1 silica gel and aluminum oxide. All separated components were analyzed with a Thermo TRACE GC Ultra combined to a DSQ II mass instrument which was fitted with a HP-5 mass spectrometry (MS) capillary column (30 m × 0.32 mm × 0.25 µm). Helium was used as the carrier gas at a current speed of 1.2 mL/min, and the injection temperature was set at 295℃. The heating procedure for saturates is set as 80℃ for 4 min, followed by heating of 4℃/min up to 295℃, and then held for 20 min. The heating procedure for aromatics as 80℃ for 4 min, followed by heating of 3℃/min up to 295℃, then held for 20 min. The light crude oils were directly used for the analyses without separation, and then the heating procedure is set as 35℃ held steady for 4 min, followed by heating at a rate of 4℃/min up to 295℃, and then held for 30 min. TZ light crude oils were analyzed using a GV isoprime GC–MS to determine the specific stable carbon (C) and hydrogen (H) isotope values of their toluene and methyl cyclohexane compounds. Thus, compounds specific C isotope values were analyzed using a GV isoprime GC–MS fitted with a capillary column HP-1MS (30 m × 0.25 mm × 0.25 µm). Finally, compound specific C isotope values were analyzed using a ThermoFisher gas chromatography GC-IRMS fitted with a HP-1MS capillary column (30 m × 0.25 mm × 0.25 µm). The heating procedure and carrier gas are the same as previously described for GC–MS.
Results and discussion
Geochemical parameters of Tazhong black oils.
G/C30 H: gammacerane/C30 hopane; C23/C21 TT: C23 tricyclic terpane/C21 tricyclic terpane; C29/C30 H: C29 hopane/C30 hopane; C27 (%), C28 (%), and C29 (%) refer to C27, C28, and C29 steranes; MDBT: methyldibenzothiophene; DBT: dibenzothiophene; MDR = 4-/1-MDBT; MPI1: methylphenanthrene indicator 1; MPI2: methylphenanthrene indicator 2.
Geochemical parameters of Tazhong light crude oils.
Abbreviations (in addition to those already listed in the text and below Table 1): Tol: toluene; MCH: methyl cyclohexane; MN: methylnaphthalene; DMN: dimethylnaphthalene; TMN: trimethylnaphthalene; N: naphthalene; P: phenanthrene.
Distribution characteristics of saturates and aromatics in TZ black oils
A total ion chromatogram (TIC) and the compounds identified in selected TZ black oil samples are presented in Figure 2; within these samples, TZ26, TZ44 and TZ241 are located within the TZ Number I fault zone, while TZ4, TZ104 and TZ404 are located within TZ Number 10 structural zone, and TZ12, TZ35 and TZ162 are located within the central uplift zone.
Selected TICs of some Tazhong black oil samples. Wells TZ26, TZ44 and TZ241 are located within the Tazhong Number I fault zone, while TZ4, TZ104 and TZ404 are located within the Tazhong Number 10 structural zone, and TZ12, TZ35 and TZ162 are located within the central uplift zone. The numbers above peaks denote n-alkane carbon numbers.
The distribution of n-alkanes showed that TZ black oils mainly ranges between C13 and C30, with the main peak between C15 and C16. The pristane (Pr)/phytane (Ph) ratio falls between 0.83 and 1.32; this ratio for samples from TZ Number 10 structural zone mainly ranges between 0.80 and 1.00, while mainly ranges between 1.00 and 1.20 for samples from TZ Number I fault zone. The UCM (Unidentified Complex Mixture) peaks in some samples suggested their suffering from biodegrading, which means that these TZ oil reservoirs may have experienced more than once charges. And it seems that the early charged oils were biodegraded before being mixed with the later charged oils.
Table 1 lists the biomarker parameters of tricyclic terpanes, hopanes, steranes, and gammaceranes. These data show that the gammacerane/C30 αβ hopane ratio for oils from the TZ Number 10 structural zone as well as the central uplift zone mostly range between 0.05 and 0.31, while samples from the TZ Number I fault zone mainly range between 0.08 and 0.21. Data also reveal a low abundance of tricyclic terpanes as well as a relatively high abundance of pentacyclic terpanes; the C23/C21 tricyclic terpane ratio ranges between 1.05 and 2.42, while that of C29/C30 hopanes ranges between 0.42 and 1.39. The distribution of steranes in these samples is C27 regular steranes between 19.00% and 50.14%, while C28 steranes occur between 16.66% and 37.63% abundance, and C29 steranes occur between 33.13% and 55.58%, which is generally characterized by an overall ratio relationship C27 > C28 < C29. At the same time, however, a few sterane parameters of these samples display the ratio C27 < C28 < C29. These results indicate hydrocarbon precursors derived from bacteria and seaweed, while differences in sterane distribution might be the result of several charging phases.
Results (Table 1) show little difference between the oil samples from different structural zones concerning the distribution of regular biomarker parameters because these parameters have been homogenized by the facts that the oil generating rocks within the Tarim Basin comprise several sources (Jin, 2005; Yun et al., 2014), which have experienced several charging phases (Chang et al., 2013; Zhao et al., 2012) and characterized by the adjustment and destruction of generated hydrocarbons (Zhuo and Wang, 1999). Thus, three tectonic phases can be identified corresponding with the end of the Early Paleozoic, the Late Paleozoic-to-Early Mesozoic, and the Late Cenozoic; these match with the oil charging phases of the Tarim Basin, specifically named as the late Caledonian-to-Hercynian, the late Hercynian and Himalayan periods (Zhang et al., 2011). Oils from the Cambrian source rocks (late Caledonian-to-Hercynian) have generally been destroyed, which is accordant to that the Late Hercynian is the main oil reserved period of the two source rocks. The Himalayan period was the main time for gas generation, as gas invading condensate gas reservoirs formed during this period (Zhang et al., 2012). In contrast, samples containing crude oils from these three generating phases also have particular characteristics as they formed in distinct structural zones, all of which were highly saline marine facies derived from bacterial and seaweed precursors.
The aromatic hydrocarbon components of samples mainly comprise naphthalene and alkylnaphthalene, as well as phenanthrene and alkylphenanthrene, fluorene and alkylfluorene, and dibenzothiophenes, particularly among which the naphthalene and phenanthrene homologues occur at higher abundances. Earlier work concerning methylphenanthrene parameters have been related to maturity (Radke et al., 1982). Indeed, both the methylphenanthrene indicator 1 (MPI1) and methylphenanthrene indicator 2 (MPI2) values of source rocks initially increase and then decrease as depth increases; calculated MPI1 values [MPI = 1.5 × (2-methylphenanthrene + 3-methylphenanthrene)/(phenanthrene + 1-methylphenanthrene + 9-methylphenanthrene)] for TZ black oils range between 0.42 (TZ452) and 1.22 (TZ35) (although mainly concentrated between 0.55 and 0.85), while MPI2 values [MPI2 = 3 × (2-methylphenanthrene)/(phenanthrene + 1-methylphenanthrene + 9-methylphenanthrene)] of these samples range between 0.47 (TZ452) and 1.36 (TZ35) (mainly between 0.69 and 0.91). Equivalent vitrinite reflectances (Ro%) calculated based on MPI (Ro% = 0.60 × MPI1 + 0.40) (i.e. 0.65 ≤ Ro% < 1.35) range between 0.73 and 1.13.TZ oils derived from marine source rock while MPI mainly works for type III kerogen, thus calculated Ro cannot well correlate to source rock and then less reliable.
In terms of sulfur-containing aromatic compounds, it is informative of different locations of methyl groups within methyldibenzothiophenes as they are indicative of thermal stability (4-MDBT > 1-MDBT). Research in this area was performed by Dzou et al. (1995) who reported that the relationship between MDR (4-MDBT/1-MDBT) and equivalent vitrinite reflectance is Ro% = 0.2663 × Ln (MDR) + 0.9034. Compounds of this type have been studied in details (Chakhmakhchev et al., 1997; Luo et al., 2001; Santamaría-Orozco et al., 1998). Zhou et al. (2008), for example, calculated the maturity of source rocks at different depths within the TC1 well using this formula. The MDR of oil samples ranged between 1.93 and 24.66, while MDR4 (4-MDBT/DBT) ranged between 0.62 and 3.66, and Ro% ranged between 1.08 and 1.76, which is more reliable than the Ro calculated by MPI. At the same time, it is also noteworthy that, in addition to maturity, aromatic sulfur-containing compound parameters may also be influenced by other factors including biodegradation.
Geochemical characterization of Tazhong light crude oils
Overall geochemical properties of the light crude oils
Some TIC results of TZ light crude oils are illustrated in Figure 3, among which ZG10, ZG70 and ZG12 are located within the TZ Number I fault zone, while ZG46, ZG47 and ZG48 are located within the TZ Number 10 structural zone, and the ZS1 reservoir is located in the central uplift zone deposited in Cambrian strata. As these samples mainly belong to light crude oils, the n-alkanes in samples are light hydrocarbons with low carbon number. As an example, the sample from well ZG12 has a peak n-alkane at C8 which demonstrates high maturity; this also means that most normally detected biomarkers are not present and so the regular geochemical parameters usually applied to sample characterization are also invalid. Rich aromatic hydrocarbons are detectable in all these samples including alkylbenzene, alkylnaphthalene and alkylphenanthrene as well as biphenyl, alkyldibenzothiophene.
Partial TICs of selected light crude oil samples. Numbers above peaks represent n-alkane carbon numbers.
For ZG17 light crude oil, the mass chromatograms exemplified by ions of m/z 91, 98, 128, 142, 156, 170, 178, 184, 192 and 198 are presented in Figure 4. Toluene and methylcyclohexane in the samples may share some common geochemical origin as they may have been derived from the same precursor molecules and even the toluene may be derived from methylcyclohexane in the oil reservoirs, and then the ratio of toluene to methyl cyclohexane will be further discussed concerning their geochemical characteristics including maturity.
Partial mass chromatograms of light crude oils from ZG17.
Geochemical characteristics of the aromatics of the light crude oils
Some geochemical parameters of aromatic compounds in TZ light crude oils are shown in Figure 5. In terms of the first group, the ratio between toluene and methyl cyclohexane ranges between 0.63 and 3.03 (mainly between 0.98 and 2.01). The ratio between these compounds in samples from the TZ Number I fault zone ranges more widely than those for the TZ Number 10 structural zone, with ratios for the latter concentrated between 0.79 and 1.30, however with the exception for ZG46 sample (also see the following discussion).
Cross-plots for Tazhong light crude oils: (A) 4-MDBT/1-MDBT ratios versus 4-MDBT/DBT ratios; (B) toluene/MCH ratios versus DBT/phenanthrene ratios; (C) naphthalene/phenanthrene ratios versus toluene/MCH ratios; and (D) DMN/TMN ratios versus MN/DMN ratios.
Some previous research has suggested that a high abundance of naphthalene series compounds is indicative of highly mature source rocks (Price, 1993). In TZ45 condensate oils and light crude oils, the abundance of naphthalene series is relatively high (Zhou et al., 2008). Similarly, in terms of methylnaphthalene (MN) molecular characteristics, methyl in location β tend to be more stable than those in location α, which means that 2-MN tend to exhibit higher thermal stability than their 1-MN counterparts. Thus, the ratio between 2-MN and 1-MN compounds can be applied for maturity comparisons of light crude oils.
Our results show that the MN indicators (2-MN/1-MN) in TZ light crude oils range between 1.21 and 2.81. Comparing these results with research on MN compounds from Li and Wang (2005) suggests that the Ro% of samples should be more than 1.0%; however, this previous result was based on low maturity samples and so a difference should be expected when compared to the high maturity oils discussed in this work. Indeed, in the case of multiple-substituted naphthalene parameters, MN/C2 substituted naphthalenes (DMN) ratios of our samples range between 0.49 and 1.30; specifically, this ratio for samples from the TZ Number I fault mainly between 0.51 and 0.88 while those from the TZ Number 10 structural zone between 0.59 and 0.75. At the same time, DMN/C3 substituted naphthalenes (TMN) ratios range between 0.83 and 3.25; this ratio in samples from the TZ Number I fault zone mainly between 0.90 and 1.78, while for samples from the TZ Number 10 structural zone range between 0.83 and 1.44.
Values of MPI were chosen to illustrate the parameters of phenanthrene series compounds. Results show that MPI1 values for light crude oils range between 0.61 and 1.18 while MPI2 values range between 0.56 and 1.35. Thus, equivalent Ro% values calculated using MPI (i.e. Ro%a = 0.60 × MPI1 + 0.40; 0.65 ≤ Ro% < 1.35; Ro%b = −0.60 × MPI1 + 2.30; 1.35 ≤ Ro% < 2.00) range between 0.72% and 1.11% and between 1.59% and 2.00%, respectively. At the same time, the ratios between naphthalene and phenanthrene range between 2.48 and 30.25; results show that N/phenanthrene ratios for samples from the TZ Number I fault zone range between 3.18 and 30.25, a much wider variation than those from the TZ Number 10 structural zone between 2.48 and 8.13 (with the exception for ZG46).
The MDR of oil samples ranged between 2.06 and 10.84, while MDR4 (4-MDBT/DBT) ranged between 0.60 and 4.20, and Ro% ranged between 1.10 and 1.51.
Our results show that the DBT/phenanthrene ratios of TZ oil samples range between 0.24 and 9.6; these ratios in samples from the TZ Number I fault zone are between 0.64 and 9.64, a wider range than samples from the TZ Number 10 structural zone of 0.80 and 3.91 (with the exception for ZG46). Reported by Zhang et al. (2015b), the extremely high DBT/phenanthrene ratios may be influenced by TSR (Thermal Sulfate Reduction). The results of MN/DMN, DMN/TMN, toluene/MCH, DBT/phenanthrene, and naphthalene/phenanthrene ratios of light crude oils from the TZ Number I fault zone range widely, while those from the TZ Number 10 fault zone encompass a relatively narrow range (Figure 5). This result indicates that oil samples from the TZ Number I fault zone have undergone more oil charging than their counterparts; indeed, even oils generated from the same hydrocarbon source rock with the same original material exhibit differences in maturity because they were generated in different hydrocarbon expulsion phases. In other words, the oil samples from the TZ Number I fault zone are more complicated in composition than their counterparts from the TZ Number 10 structural zone, as illustrated by the geochemical parameters discussed in this study.
Qiu et al. (1997, 2010, 2012) and Liu et al. (2017) performed a lot of work in investigating hydrocarbon charge. In earlier work, Qiu et al. (1997, 2010, 2012) analyzed a series of integrated thermal indicators of apatite and zircon (U–Th)/He ages, as well as apatite fission tracks, and equivalent vitrinite reflectance data from TZ source rocks. Their results were used to hypothesize that Indosinian and Yanshan movements lifted strata in this region and made multistage accumulation possible as source beds moved into a detained stage and stopped generating oil and gas following rapid source rock maturation. In a similar study, Yang et al. (2011) performed research on fluid inclusions within a TZ Ordovician carbonate reservoir and suggested three distinct phases of oil charging, namely Caledonian-to-early Hercynian, late Hercynian, and Himalayan stages. The last phase of the Himalayan, was mainly charged with gas. Thus, as the characteristics of TZ light crude oils were influenced by oil charging during Himalayan gas invading stage, they share similar features with condense oils as the late Hercynian. In contrast, as oil charging during the Caledonian was destroyed this leads to degradation characteristics if original crude oils are considered. In earlier work, Zhang et al. (2009) also considered different structural zones during a variety of evolutionary phases can influence oil accumulation. These researchers showed that the TZ Number 10 structural zone shared the same accumulation phase as the TZ Number I fault zone, but that the tectonic activity within the structural zone was less than that within the fault zone. This resulted in an intact clastic rock layer which hindered oil migration above the limestone. The advantage of this migration within the TZ Number I fault zone is that oils generated there tend to be more complex in their composition than their counterparts formed within the TZ Number 10 structural zone. These differences also mean that the geochemical characteristics of oils formed within the TZ Number I fault zone tend to vary more widely.
A large number of strike slip faults formed within the TZ area caused marked improvements in the growth of reservoir carbonate voids and stems linking different reservoirs (Wu et al., 2012; Yang et al., 2007; Zhou et al., 2013). At the same time, these strike slip faults controlled the migration and accumulation of oil into reservoirs. Most notably, the end of Silurian-to-early Devonian and late Permian periods were two most active phases of strike slip faults within the Tarim Basin (Jia, 2004). Zhang et al. (2008) reported that middle-to-upper Ordovician gas reservoirs formed by the pyrolysis of paleo-oil accumulations underneath Cambrian salt deposits as the growth of main and strike slip faults provided a number of possibilities for hydrocarbons to migrate upwards. As shown on a geological map of this area (Figure 1), light crude oils within the TZ Number 10 structural zone are mostly located within the boundaries of the main fault while there are more strike slip faults within the TZ Number I fault zone. These smaller faults provided a large number of effective channels connecting different reservoirs and resulted in several phases of oil charging (Li et al., 2006); this range of crude oils with markedly different geochemical characteristics is the cause of the parameter variation seen within this region. Notably, although well ZG46 is located within the TZ Number 10 structural zone, a number of strike slip faults connect it to the TZ Number I fault zone nearby; this might be the cause of several distinct episodes of oil charging. In addition, structures within the Tarim Basin tend to become less active from east-to-west and the ZG46 well is located within the western TZ Number 10 fault, a possible advantage that prevented the degradation of high maturity oils. Indeed, a highly mature oil reservoir in the eastern TZ Number 10 structural zone was destroyed as it lacked good cover even though strike slip faults were present in the vicinity. These structural observations explain why the toluene/MCH ratio of ZG46 oils sample is 3.03, the naphthalene/phenanthrene ratio is 14.71, and the DBT/phenanthrene ratio is 3.91, markedly higher than the values recorded for other oil samples within the TZ Number 10 structural zone.
Geochemical features of toluene and methylcyclohexane in the light crude oils
Regarding the geochemical characterization of toluene and methyl cyclohexane in TZ light crude oils, data demonstrate that saturated ring hydrocarbons can transform into aromatic hydrocarbons by dehydrogenating as maturity gets higher. This is at least one important way in which aromatic hydrocarbons are generated. Several recent studies have shown that alkyl cyclohexane can transform into alkyl benzene series compounds by dehydrogenating (Darwin et al., 2016a, 2016b). Based on these results, it is supposed that a close evolutionary origin may exist between toluene and methyl cyclohexane in light crude oils, maybe some toluene can be directly derived from methylcyclohexane, and basically in crude oils toluene and methylcyclohexane should have been derived from the same or some likewise structural units of precursor molecules. Thus, during geological evolution process, changes in relative concentration and stable isotope values of toluene and methylcyclohexane can be used to reveal the geochemical characteristics of high mature oils. Therefore, some more discussion has been made to the TZ light crude oils concerning their toluene and methylcyclohexane.
There are more methylcyclohexane than n-heptanes in samples from the TZ Number 10 fault zone (Figure 6), which is accordant to the characterization of marine crude oil. But this relationship is different in samples from the TZ Number I fault zone, among which some samples have more methylcyclohexanes while others have more n-heptane (Figure 6). These results are accordant to the fact that some light crude oils from the TZ Number I fault zone had contributed to some oil-cracking gas, and then depleted in methylcyclohexane (Hu et al., 2005; Yu et al., 2013), and therefore enriched in toluene showing a higher maturity as shown in Figure 6. Thus, these results reveal that the TZ Number I fault zone had undergone more complicated geological evolution. As discussed above, ZG46 was obviously enriched in toluene because it is at a very special tectonic location, so the ratio between toluene and methylcyclohexane can be used to distinguish depositional environments and differentiate the maturity of light crude oils.
Plot of concentration relationship between n-heptane, methylcyclohexane and toluene in Tazhong light crude oil.
Stable carbon/hydrogen isotope values of toluene and methylcyclohexane of the light crude oils.

Graphs to show: (A) δ13C of toluene versus methylcyclohexane; (B) δD of toluene versus methylcyclohexane; (C) δ13C versus δD of toluene; and (D) δ13C versus δD of methylcyclohexane in Tazhong light crude oils.
The data distribution of δ13C and δD in toluene and methyl cyclohexane (Figure 7A,B) fall around a line with a slope of 1, while the δD of toluene and methyl cyclohexane in Figure 7(B) fall below this line, especially the samples from TZ Number 10 structural zone. This result suggested that the δD of toluene is generally heavier than those of methyl cyclohexane. As discussed above, one important way that aromatic hydrocarbons are generated via the transmutation of saturated ring hydrocarbons in crude oils via dehydrogenation as maturity increases. This will cause the breakage of C–H bonds in the molecule and then enriched in deuterium for toluene.
No distinct correlation between δ13C and δD in either toluene or methylcyclohexane is evident from the TZ light crude oils, as the former is sensitive to both precursors and depositional environments but less so to thermal evolution, whereas the latter is more sensitive to thermal stress. After oil charging, crude oils containing different kinds of compounds from different charging phases will exhibit different stable isotope features, which will be commingled into the geochemical outcomes of the light crude oils. So in Figure 8 these characteristics are reflected on plots of toluene and methylcyclohexane ratios versus δ13C and δD in TZ light crude oils. No matter how these ratios vary, no obvious regular pattern is found in our data between toluene and methylcyclohexane and their δ13C or δD, so long as crude oils contain compounds charged during different geological phases when they were derived from different precursors and depositional environments. However, in the geochemical study on light crude oils, toluene and methylcyclohexane should be investigated in details as they are abundant in the light crude oils and may be with genetic relationships, so further work will be helpful to the study on light crude oils from the composition, stable carbon/hydrogen isotope analyses of toluene and methylcyclohexane.
Plots of toluene/methylcyclohexane ratios versus: (A) δ13C values of toluene; (B) δD values of toluene; (C) δ13C values of methylcyclohexane; and (D) δD values of methyl cyclohexane in Tazhong light crude oils.
Conclusions
A total of 89 crude oils were analyzed in this study with emphasis on TZ light crude oil aromatic hydrocarbons from different structural zones. The aim of this study was to evaluate the geochemical parameters of these crude oils and to investigate their stable carbon/hydrogen isotope ratios.
Results reveal that the crude oils from TZ area are highly mature, especially true for the light crude oils. At high mature stage the normal biomarker parameters are invalid and then the aromatic hydrocarbon compounds are the main study targets in this work. The distribution of aromatic hydrocarbons is different between crude oils from different structural zone. Thus, the MN/DMN, DMN/TMN, MCH/Tol, DBT/phenanthrene, and naphthalene/phenanthrene ratios of crude oils from the TZ Number I fault zone range widely while those from the TZ Number 10 structural zone encompass a relatively narrow range of variation. This difference highlights the fact that the TZ Number I fault zone experienced more intense structural movements than the TZ Number 10 structural zone. The former also formed a more effective migration pathway as a result; thus, the crude oils that formed within the TZ Number I fault zone have undergone several charging phases and are more complicated in composition.
The distribution of stable carbon/hydrogen isotope values revealed by this study indicates a relationship between the origin of toluene and methylcyclohexane. Toluene can be derived from the dehydrogenation of methylcyclohexane in crude oils. Nevertheless, the ratios of toluene and methylcyclohexane are complicated within the TZ area and no regular pattern is evident due to multiple charging events and the complex geological evolution of the TZ area. The results of this study nevertheless reveal that the geochemical characteristics and significance of toluene and methylcyclohexane in light crude oils are worthy of further study.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship and/or publication of this article: This work was financially supported by the National Natural Science Foundation of China (Grant No. 41472109), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA14010103) and the National Oil and Gas Project of China (2017ZX05008002).
