1 Introduction

Heavy oil and oil sand bitumen dominate the world’s oil inventory (Marcano et al. 2013), even exceed the normal oil in quantity (Wang et al. 2016; Li and Huang 2020). Biodegradation may be responsible for these unconventional resources, which primarily formed from the microbial degradation of the conventional oils over geological timescales (Larter et al. 2012). With the proceeding of biodegradation, the oil will exhibit decreases in mass and net volume, increases in oil density, oil viscosity, oil acidity and sulfur content, enrichments in nitrogen, sulfur and oxygen-containing organic compounds, and trace mental elements, which heavily affect the oil’s chemical compositions and physical properties (Larter et al. 2006; Zhou et al. 2008; Forsythe et al. 2019). Even within a single petroleum reservoir, biodegradation may create orders of magnitude variation in oil viscosity and fluid property laterally and vertically (Adams et al. 2012). Allowing for their great impact on the economic value and predictability, knowledge of the petroleum biodegradation is critical in selecting exploitation strategies (Marcano et al. 2013; Yin et al. 2013), evaluating the altering extent (Bautista et al. 2015) and forecasting the oil compositional characteristics and viscosity (Larter et al. 2012; Zhang et al. 2014a, b).

Experimental examinations and case studies indicate that relative sensitivity to microbial attack exist among different compound class, even among the individual compound in each class. Therefore, biodegradation is a quasi-stepwise process with consume of biomarkers in a preferential order (Peters et al. 2005). Based on the alteration/removal of saturated biomarkers, i.e., normal alkanes, acyclic isoprenoids, terpanes and steranes, and selected aromatic compounds (mainly aromatic steroids), Peters and Moldowan (1993) developed a classical scale to assess the extent of oil biodegradation, which is typically abbreviated as PM scale. Five terms, i.e., light (PM1–3), moderate (PM4–5), heavy (PM6–7), very heavy (PM8–9) and severe (PM10) are assigned to ranges of the PM scale. Another scheme was established by Wenger et al. (2002) fundamentally based on the alteration extent with a compound class and the presence or depletion of single key compound class, only with the terms of heavy and severe to describe the PM 4–10 range.

The normal alkanes are generally more susceptible to biodegradation than the branched alkanes and isoprenoids (Chen et al. 2017), thus ratios of i-C5/n-C5, pristane/n-C17 and phytane/n-C18 are widely used to assess the alteration content of light to moderate biodegraded oils (Peters and Moldowan 1993). However, some case studies showed no evident correlations between Ph/n-C18 and oil API, and then the parameter “mean degradative loss” was proposed to assess the biodegradation extent by means of quantifying the depletion in volumetrically important individual oil constituents, especially for the light to moderate level of alteration (Elias et al. 2007). Additionally, the production and biodegradation of oxygen-containing compounds also follow a preferential order, thereafter the acyclic (DBE 1)/cyclic (DBE 2-4) can reveal the biodegradation level, which well correlated with the PM scale when it less than PM < 6 (Angolini et al. 2015).

Many case investigations suggest that the sequential microbial degradation of oil constitutes does not occur in a true stepwise fashion and strictly follow the schemes previously developed (Bennett and Larter 2008; Wang et al. 2013; Chang et al. 2018). All these scales are unsuitable for heavy oil or super heavy oil, and for the mixtures of oils (Larter et al. 2012; Zhang et al. 2014a, b). Selecting 8 compound classes (alkyl toluenes, C0-1 naphthalenes, C2 naphthalenes, C3 naphthalenes, methyl dibenzothiophenes, C4 naphthalenes, C0-2 phenanthrenes and steranes) to reflect the increasing resistance to biodegradation, assigning 0–4 scores to describe the removal extent, Larter et al. (2012) proposed a Manco scale to assess the biodegradation level (PM4–8), with its MN2 value positively correlates with oil viscosity, and effectively used (López 2014).

Till now, published studies mainly focused on the geochemical behaviors of molecular biomarkers at the biodegradation level of PM8 or less, for those biodegradation level higher than PM8, rare studies were reported. With the reduction in conventional resources, extremely biodegraded oils continuously discovered in the petroleum exploration attracted much more attention. Using semi-quantitative parameter-stripping method and refined Manco scale, this paper aims to deeply understand and quantitatively determine the biodegradation extent of heavy oils in the Chepaizi Uplift, Junggar Basin, NW China.

2 Geological setting

The Chepaizi Uplift, covering an area of 10,500 km2, is a petroliferous target of the northwestern Junggar Basin (Zhao et al. 2019), which is surrounded by two hydrocarbon generative source kitchens, i.e., the Changji Sag to the east and Sikeshu sag to the south (Fig. 1). The early stage of late Hercynian tectonic movement witnessed the formation of prototype of the Chepaizi Uplift, which subsequently experienced the intense uplift during the Indosinian and Yanshan movements, slow subsidence during the Himalayan tectonic movement. As an inherited paleo-uplift, Chepaizi Uplift was characterized by the seriously lack of deposition succession in the structurally high parts, with the Cretaceous, Paleogene, Neogene and Quaternary sedimentary rocks overly Carboniferous volcanics (Song et al. 2007).

Fig. 1
figure 1

Map showing structural elements of the Junggar Basin (a), with location of the Chepaizi Uplift (b), and a field-scale map showing the sampling wells in the Chepaizi Uplift (c). This figure was modified after (Chang et al. 2018)

Due to the multi-staged transgression and uplifting, three sets of favorable reservoir rocks were developed in the Chepaizi Uplift, i.e., medium-fine sandstone and glutenite occurred in the Neogene Shawan Formation, fine sandstones and siltstone occurred in the Cretaceous interval, and basalt, andesite, tuff and volcanic breccia in the Carboniferous interval.

The Hongche Fault, an active fault separating the Chepaizi Uplift from the Changji Sag (Dong et al. 2017; Ni et al. 2019), cut through sandbodies developed in the Cretaceous, Neogene Shawan Formation and volcanics in the Carboniferous, providing a favorable vertical conduit for the migration of hydrocarbon generated in the Changji Sag to the Chepaizi Uplift and yielding three vertical oil-bearing intervals (Miao et al. 2015; Meng et al. 2016).

3 Samples and experimental

3.1 Sample preparation

Eleven DST oil samples from the Carboniferous volcanic interval were collected from wells in the Chepaizi Uplift for geochemical analysis. Asphaltenes were removed from the oil samples and source rocks using n-hexane precipitation, and the deasphaltened oil was divided into two aliquots. The first of these underwent column chromatography using a routine silica gel and alumina column from which aliphatic and aromatic fractions were obtained using n-hexane and dichloromethane (DCM).

3.2 Gas chromatography–mass spectrometry

Gas chromatography–mass spectrometry (GC–MS) analysis of the aliphatic and aromatic fractions was performed with a Finnigan Model SSQ-710 quadrupole analytical system coupled to a DB-5 fused silica column (30 m × 0.32 mm i.d.) and linked to an IAIS data processing system. GC temperature operating conditions for the aliphatic fraction were as follows: 100 °C (1 min) to 220 °C at 4 °C/min and, then to 300 °C (held 5 min) at 2 °C/min; for the aromatic fraction: 80 °C (1 min) to 300 °C (held 15 min) at 3 °C/min. MS conditions were as follows: electron impact (EI) ionization mode; 70-eV electron energy; 300-mA emission current; and 50–550 amu/s scan range.

4 Results and discussion

4.1 Oil bulk compositions and oil families

The Carboniferous oils from the eastern Chepaizi Uplift are characterized by higher oil density (0.9285–0.9590 g/cm3) and viscosity (154–8968 mPa·s) than those from the western Chepaizi Uplift, which can be classified as heavy crude oils (Table 1). The Carboniferous oils are predominantly aliphatic, as indicated by their saturate/aromatic (ST/AR) ratio (1.92–3.56) and saturate fraction abundance (44.45%–64.38%). The oil density showed roughly positive correlation with the viscosity and the NSO (resin + asphaltene) fraction content, and negative correlation with the burial depth. Progressive biodegradation of crude oils may be responsible, as it decreases the content of saturate and aromatic hydrocarbons and enriches the resins and asphaltenes, resulting in an increase in oil density (López et al. 2015; Wenger et al. 2002).

Table 1 Bulk compositions of crude oils from Chepaizi Uplift. ST, saturate hydrocarbon; AR, aromatic hydrocarbon; NSO, resin + asphaltene; G/C30H, gammacerane/C30 hopane; *, data cited from Chang et al. (2018) and Xu et al. (2018)

As for the study of oil origins, oil family is a widely used term to distinguish oils with different genetic affinity. An oil family, a group of oils that derived from a same source rock, possibly experienced similar reservoir-forming history, belonged to a same oil system, and possessed same or similar chemical compositions. Combing the ratios of C27 diasterane 20S/20R (C27DS 20S/20R), C27 diasterane/C27 regular sterane (C27DS/C27RS), C26/C28triaromatic steroid (20S) (C26/C28TAS(20S)), and C27/C28TAS(20R) with the stable carbon isotope distribution, previous studies concluded that the Carboniferous oils in the eastern Chepaizi Uplift were mainly derived from the mudstone of the Middle Permian Wuerhe Formation (P2w) in the Changji Sag, whereas the ones in the western Chepaizi Uplift were essentially originated from the Jurassic mudstone in the Sikeshu Sag (Zhang et al. 2012; Xu et al. 2018; Mao et al. 2020). Although two oil charging episodes were defined in the Carboniferous reservoirs based on the oil geochemistry, fluid inclusions and basin modeling (Chang et al. 2019; Shi et al. 2020), the later charge actually was the remigration of early reservoired oils due to the tectonic adjustment (Cao et al. 2010; Song et al. 2016; Chang et al. 2019). Although the Carboniferous oils all exhibited the characteristics of lacustrine source facies, source-diagnostic and redox potential of depositional environment-related molecular biomarkers showed marked distinction between the eastern and western parts (Xu et al. 2018), implying at least two oil groups.

Cluster analysis, a statistical method, is an effective tool to classify studied samples into different groups by their similarity distances which were calculated from the different variables investigated. Allowing for the severe oil alteration, eight biomarker parameters strongly resistant to biodegradation were selected to act as the variables in the clustering analysis, i.e., C26/C28TAS (20S), C27/C28TAS (20R), C24 tetracyclic terpane/C30 hopane (C24Tet/C30H), gammacerane/C30 hopane (G/C30H), C35H (22S)/C34H (22S), C29 18α(H)-30- norneohopane/C29 Hopane (C29Ts/C29H), C24 tricyclic terpane (TT)/C23TT, C22TT/C21TT. In the clustering tree graph (Fig. 2), the Carboniferous oils can be clearly divided into two groups, that is, Group I for the western Chepaizi Uplift, Group II for the eastern Chepaizi Uplift. In addition, within the Groups II, three oils families can be further subdivided (II1, II2 and II3) according to their similarity distances.

Fig. 2
figure 2

Clustering tree graphs showing the different oil groups (a) and oil families (b)

4.2 Biodegradation level by PM scale

4.2.1 Qualitative evaluation by molecular compositions


  1. 1.

    Oil Group I


Group I (wells P70 and 702), Carboniferous oils from the western Chepaizi Uplift, was least biodegraded among the investigated oils. These oils showed faintly “UCM” and relatively intact n-alkanes on the TIC chromatogram (Fig. 3a). Hopanes were slightly higher than the tricyclic terpanes (TTs) in abundance with C30 hopane as the peak compound and reversed “L” -shaped distribution of C20TT-C21TT-C23 TT (Fig. 3b). Gammacerane (G) was low in content as evidenced by the quietly low G/C30H values (0.08–0.09). Fully developed 25-norhopanes were detected indicating heavy biodegradation. Pregnane, homopregnane, and diasteranes (DS) were lower than the regular steranes (RS) in abundance with “V”-shaped distribution of ααα20R C27-C28-C29 steranes (Fig. 3c). Stable carbon isotope varied from −27.5 ‰ ~ −27.8  ‰ (Table 1). Naphthalenes and phenanthrenes were nearly intact and triaromatic steroids (TAS) were essentially unchanged (Fig. 3d). Thereby, Group I can be assigned to biodegradation level of PM6.

Fig. 3
figure 3

Representative fragmentograms of oil group I showing molecular compositions (well P70). a Total ion chromatogram (TIC); Pr, pristane; Ph, phytane. b m/z 191, C19-C26, tricyclic terpanes with different carbon number; Tet, tetracyclic terpane; Ts, 18α(H)- trisnorneohopane; Tm, 17α(H)- trisnorhopane; 25-NH, C29 25-norhopane; C29H, C29 hopane; C30H, C30 hopane; G, gammacerane. c m/z 217, C21P, C21 pregnane; C22P, C22 homopregnane; and d m/z 231, C20–C21, C20–C21 triaromatic steroids; C26–C28 20S and 20R,20S and 20R isomers for triaromatic steroids with different carbon numbers

  1. 2.

    Oil group II- family II1


Family II1 (wells P661, P666 and P66), featured substantially removed n-alkanes and iso-alkanes and prominent “UCM” on the TIC chromatograms (Fig. 4a). C31–35 homohopanes were heavily depleted. Tricyclic terpanes were far higher than hopanes in abundance with reversed “V” distribution of C20TT–C21TT–C23TT (Fig. 4b). Gammacerane was relatively high in content as showed by the high G/C30H values (0.49–0.54). Pregnane and homopregnane were enriched and nearly equal to the regular steranes in abundance with reversed “L” distribution of ααα20R C27–C28–C29 sterane (Fig. 4c). Stable carbon isotope varied from −29.9‰ to −30.0  ‰ (Table 1). Naphthalenes and phenanthrenes were substantially depleted (Fig. 4d). The biodegradation can be classified into PM6 + to PM7.

Fig. 4
figure 4

Representative fragmentograms of oil family II1 showing molecular compositions (well P661). a Total ion chromatogram (TIC), RS, regular sterane, b m/z 191, c m/z 217, and d m/z 231. Abbreviations are the same to Fig. 3

  1. 3.

    Oil group II- family II2


Family II2 (wells P663, P668, P61 and P60), showed heavily removed hopane series with C29 hopane as the peak compound, equal tricyclic terpanes to hopanes in abundance with sequentially increased distribution of C20TT–C21TT–C23 TT (Fig. 5b), and high G/C30H ratios (0.41–0.66). Pregnane and homopregnane far exceeded the regular steranes with “L” distribution of ααα20R C27–C28–C29 sterane (Fig. 5c). The diasteranes were prominently altered. Stable carbon isotope (−29.7‰ to −30.7‰) was like that of the family II1 (Table 1). Triaromatic steroids were still unchanged (Fig. 5d). The biodegradation level reached PM8 to PM8+ .

Fig. 5
figure 5

Representative fragmentograms of oil family II2 showing molecular compositions (well P663). a Total ion chromatogram (TIC), b m/z 191, c m/z 217, and d m/z 23. All abbreviations are the same to Fig. 3

  1. 4.

    Oil group II- family II3

Family II3 (wells P665 and P685), by contrast, exhibited the strongest biodegradation characteristics. Hopanes were essentially removed with C2925-norhopane as the peak compound, and tricyclic terpanes were prominently reduced (Fig. 6b). Gammacerane was abnormally high with G/C30H values ranging from 3.93–4.27. Pregnane and homopregnane were far beyond the majorly depleted regular steranes (Fig. 6c). However, no changes of the stable carbon isotope (−29.7‰ to −30.7‰) can be seen (Table 1). The essentially unaltered triaromatic steroids (Fig. 6d) confirmed the biodegradation level of PM9+ .

Fig. 6
figure 6

Representative fragmentograms of oil family II3 showing molecular compositions (well P665). a Total ion chromatogram (TIC), b m/z 191, c m/z 217, and d m/z 23. All abbreviations are the same to Fig. 3

  1. 5.

    Sequential biodegradation revealed by biomarkers

Although many studies reported that the microbial degradation of oil did not appear to be strictly consistent with the stepwise fashion in established schemes, relative variations of biomarkers in oil families II1, II2 and II3 can be prominently distinguished (Fig. 7). The values of C22T/C21TT, C24TT/C23TT, C24C24Tet/C26TT and (C20 + C21)/C26TT ratios increased with the biodegradation extent increased from PM7 to PM8 and finally to PM9+ (Fig. 7a–d), suggesting the alterations of tricyclic terpanes, preferential removal of lower molecular weight homologues and more bioresistance of tetracyclic terpane (Huang and Li 2017). Hopane/sterane ratio varied with biodegradation in a zigzag-shaped fashion (Fig. 7e), i.e., essentially unchanged from the level of PM6 + to PM7, sharply increased from the level of PM7 to PM8, and progressively decreased from the level of PM8 to PM9+, confirming the different susceptibility to biodegradation for hopanes and steranes at different biodegradation level (Reed 1977;). Ratios of C29/C30H, G/30H and C30 diahopane/C30 hopane (C30D/C30H) kept relatively unchanged from level PM7 to PM8, then sharply increased to PM9 + (Fig. 7f–h), indicating approximately the same bioresistance of C30 hopane to C29 hopane, gammacerane and C30 diahopane before PM8 level, but more susceptibility at level PM9 + . Gradually increased 18α(H)- trisnorneohopane/17α(H)- trisnorhopane (Ts/Tm) ratio with biodegradation indicated the relatively more susceptible of Tm than Ts to microbial alteration (Fig. 7i). Increasing C29 25-norhopane/gammacerane (C29NH/G) value with biodegradation confirmed the formation of 25-norhopane (Fig. 7j) at level above PM6, which was consistent with the increasing C28NH/C29H and C29NH/C30H (Fig. 7k, 7l). However, its decrease from PM8 to PM9 + revealed the degradation of 25-norhopanes occurred under extreme biodegradation level (Huang and Li 2017; Chang et al. 2018; Killops et al. 2019), as validated by the increasing C29NH/C28NH (Fig. 7m). 25-norhopane ratio (∑C30–C34 25-norhopanes/(∑C30–C34 25-norhopanes + ∑C31–C35 homohopanes) generally exhibited a tendency to be increased with biodegradation, which was consistent with the formation of 25-norhopanes and reduction in homohopanes (Fig. 7n). Increasing C27DS/C27RS with biodegradation displayed the faster degradation of regular steranes than the diasteranes after PM7 level (Fig. 7o). Unchanged C26/C28TAS(20S) and C27/C28TAS(20R) well documented the strongest bioresistance of triaromatic steroids (Fig. 7q, 7r), yet the decreasing (C20 + C21)/C26–28TAS demonstrated that the short-chained counterparts were slightly destroyed above PM7 (Fig. 7p).

Fig. 7
figure 7

Correlations of biomarker parameters among three oil families of Group II showing the alterations of biomarkers with different susceptibility to biodegradation. TT, tricyclic terpane; Tet, tetracyclic terpane; H, hopane; G, gammacerane; Ts, 18α(H)-trisnorneohopane; Tm, 17α(H)-trisnorhopane; NH, 25-norhopane; C30D, C30 diahopane; ∑25-NHs, sum of C30-C34 25-norhopanes; ∑HHs, sum of C31-C35 homohopanes; DS, diasterane; RS, regular sterane; TAS, triaromatic steroid. All the ratio data were cited after (Chang et al. 2018)

4.2.2 Semi-quantitative evaluation by parameter-stripping method


Ideally, biomarkers were sequentially removed by biodegradation in a stepwise fashion and can be described by the ten-point scales (Peters and Moldowan 1993). To further proof the extent of biodegradation four biodegradation stages were defined following the general alteration tendency in PM scheme (Fig. 8). Stage 1 covers approximately from PM1 to PM5, stage 2 from PM5 to PM7, stage 3 from PM7 to PM8, and stage 4 involves PM8 and higher. Biodegradation levels of different oil families were semi-quantitatively defined by sequentially stripping the samples in the cross-plots of biodegradation-resistant parameters.

Fig. 8
figure 8

Schematic plot of the four-stages of the PM scale

  1. 1.

    Stage 4

Diasteranes show particular resistance to biodegradation and remain where steranes and hopanes are totally removed in case of no 25-norhopanes are present (PM9) (Peters et al. 2005). Pregnane (P) and homopregnane have high resistance to biodegradation, comparable to the diasterane (PM9). Tricyclic terpanes can be altered at about the same level as the diasteranes (Lin et al. 1989). Non-hopanoid triterpanes, such as gammacerane, diahopanes, oleanane, are highly resistant to biodegradation (Peters et al. 2005), and even beyond the point (PM ≥ 9) where the tricyclic terpanes have been removed (Wenger and Isaksen 2002). Although aromatized steroids keep unaltered in all but the most severely biodegraded oils (PM10) and can be effectively used to determine the origin and thermal maturity for extremely biodegraded oils (Peters and Moldowan 1993), low molecular weight triaromatic steroids (C20-C21) are among the first aromatic steroids to be depleted during biodegradation (Wardroper et al., 1984). From the cross-plots of gammacerane/C30 diahopane (G/C30D) vs. C27DS(20S)/C27DS(20R), C26/C28TAS(20S) vs. C27/C28TAS(20R), Carboniferous oils in the western and eastern Chepaizi Uplift were clearly distinguished exactly corresponding to the Group I and II, respectively (Fig. 9a, 9b).

Fig. 9
figure 9

Scattering plots of the Carboniferous oils by the parameter-striping method. DS, diasterane; G, gammacerane;/DH, diahopane; TAS, triaromatic steroid; TT, tricyclic terpane; P, pregnane; RS, regular sterane; 25-NH, 25-norhopane; H, hopane; C29Ts, C29 18α-30-norneohopane; HHI, C35 homohopane/sum of C31-C35 homohopanes; ETR, (C28 tricyclic terpane + C29 tricyclic terpane)/(C28 tricyclic terpane + C29 tricyclic terpane +18α(H)-trisnorneohopane); TeBR, 1,3,6,7-/1,3,5,7 tetramethylnaphthalene; DNR1, (2,6- +2,7-)/1,5-dimethylnaphthalene; *, data cited from Chang et al. (2018)


Using triaromatic steroids ratios, previous studies concluded that the Carboniferous oils from the eastern Chepaizi Uplift had a common origin and similar maturities (Xi et al. 2014; Xiao et al. 2014; Xu et al. 2018; Chang et al. 2019; Mao et al. 2020); therefore, the biodegradation extent may be responsible for the observed differences of oil compositions. The tricyclic terpanes feature similarly high bioresistance to the diasteranes (PM > 8; Lin et al. 1989), slightly lower than the non-hopanoid terpanes. In the extreme biodegradation cases, tricyclic terpanes with lower molecular weight would be preferentially degraded (Huang and Li 2017; Chang et al. 2018). Oil family II3 exhibited prominently higher C24/C23TT and C22/C21TT ratios than other oil families, implying a biodegradation level at least up to PM 9–PM9+ (Fig. 9c). Obviously, by contrast, the biodegradation rank of family II1 and II2 were lower than PM9.

  1. 2.

    Stage 3

Hopanes are removed before or after steranes, 25-norhopanes occur in oils where the hopanes are preferentially removed (Reed 1977; Peters et al. 2005), as supported by this investigation (Fig. 7e). C31 and C32 homohopanes were more susceptible to biodegradation than C30 hopane in the asphalts from Madagascar (Rullkotter and Wendisch 1982), while C28-C30 17α-hopanes are typically biodegraded in the same manner and at approximately the same rate as the C31-C35 extended hopanes (Williams et al. 1986). The 25-norhopane ratio, that is ratio of the total C30–C34 25-norhopanes to the sum of these compounds plus the C31–C35 homohopanes, could be used to evaluate the alternation extent among severely biodegraded oils of PM ~ 6–9 (Peters et al. 1996). Oil family II2, featured medium ratios of C21 P/∑C27–29 RS and C22 P/∑C27–29 RS (Fig. 9d), which clearly distinguished from those of the family II3 (quite high) and family II1 (rather low), and further verified by the varying range of C2925-NH/C30H and C3025-NH/C31H ratios (Fig. 9e). Hence, the biodegradation level of the family II2 can be evaluated at PM8–PM8+ , while the family II1 less than PM8.

  1. 3.

    Stage 2

Normal alkanes are preferentially removed at the biodegradation level of PM1–2; however, selective biodegradation of the isoprenoid over more bioresistant steranes or hopanes is used to determine the level of PM3–4 (Peters et al. 2005). Long-chained alkylated cyclopentanes and cyclohexanes are about as susceptible as branched alkanes and monocyclic alkanes will be depleted or in trace quantities at PM3–4 level (Peters and Moldowan 1993). C14-C16 bicyclic terpanes are less susceptible to biodegradation than isoprenoid and will be completely eliminated before the start of sterane and hopane biodegradation (Peters et al. 2005). Regular steranes are removed after the complete removal of C15-C20 isoprenoids and before or after the hopanes (Peters et al. 2005). Comparatively, slight alteration of methyl and dimethylnapthalenes occurs during the removal of n-alkanes, trimethylnaphthalenes are altered during the removal of the isoprenoids, and tetramethylnaphthalenes persist until steranes are largely depleted (Fisher et al. 1998). Phenanthrenes are generally more resistant to biodegradation than alkylnaphthalenes, and methylbiphenyls, dimethylbiphenyls, and methyldiphenylmethanes lacks in the biodegraded oils at PM7 (Peters and Moldowan 1993). Higher-hopane homologues, particularly the C35 pentahomohopanes, are preferentially bioresistant, and the alteration of C35 homhopanes are possibly at the rank of PM ≥ 7 (Seifert et al. 1984; Moldowan et al. 1995). The G/C30H and C29 hopane/C30 hopane (C29H/C30H) varied from roughly unchanged to rapidly increased with the increasing biodegradation level, confirming the faster degradation of C30H than C29H, especially at level higher than PM8 (Bennett et al. 2006; Chang et al. 2018); however, the substantial constancy of C29 18α-30-norneohopane/C29 hopane (C29Ts/C29H) implies the similar susceptibility of C29Ts and C29H to biodegradation (Chang et al. 2018). Oil family II3 was biodegraded to PM7 by its lower C35/C34 hopane (22S) and HHI values (Fig. 9f, g). Furthermore, the tetramethylnaphthalene ratio (TeBR = 1,3,6,7-TeMN/1,3,5,7TeMN, TeMN = tetramethylnaphthalene; Peters et al. 2005) begins to be altered significantly at levels PM5–6 (Fisher et al. 1998). Therefore, oil family II1 distinctly distinguished from other families and showed biodegradation level of PM7 as evidenced by the higher TeBR values (Fig. 9h).

  1. 4.

    Stage 1

No biodegradation occurred.

4.3 Quantitative evaluation by refined Manco scale

4.3.1 Refined Manco scale

The present biodegradation scales are essentially established according to the presence or absence of single compound classes and the alteration extent within a compound class, which encountered issues for super heavy oils (Larter et al. 2012). Manco scale emerged as the requirement, integrating the extent of degradation of various sets of compound classes. This quantitative Manco scale provided a higher resolution, however, only effective for the biodegradation of PM4–8.

More biodegradation refractory compound class, i.e., 25-norhopane, tricyclic terpanes and triaromatic steroids, were added to refine the Manco scale to meet the situation of very heavy to severe biodegraded oils from Chepaizi Uplift. Followed the academic principles of Manco scale, eight compound classes that covered the whole range of biodegradation rank (PM0–10) and possessed increasing resistance to biodegradation were selected as vector elements. The eight compound classes and their GC–MS detection m/z values are listed in Table 2.

Table 2 Refined compound classes used to represent increasing bioresistance

Five levels of refined Manco score (RMS) from 0–4 were assigned to the 8 compound classes, referencing after the Manco scale (Table 3). The refined Manco number (RMN1 and RMN2) can be calculated by the following formulas.

$$RMN_{1} = \sum m_{i} 5^{i}$$
$$RMN_{2} = [n + (\log_{5} \left( {RMN_{1} } \right)\cdot\left( {S_{\hbox{max} } - 1} \right)]/n$$

where m refers to the refined Manco score for each out of the eight vector elements (0–4), i means the class number (0–7), and n is the number of compound classes. Smax, maximum for RMN2, is designate to be 1000 to avoid confusion with the currently existing scales, and to ensure enough resolution at different levels of biodegradation when using integer values.

Table 3 Refined Manco Scores used to qualitatively distinguish the biodegradation level of compound classes

When applied the refined Manco scale to evaluate the Carboniferous oils in the Chepaizi Uplift (Table 4), the RMN2 generally well correlated with the viscosity and NSO content (Fig. 10a, b), implying the loss of oils, especially those light ends or susceptible compounds with increasing biodegradation extent. Summed pentacyclic terpanes (∑PTs), summed tricyclic terpanes (∑TTs) and diasteranes concentrations decreased with the increasing RMN2 (Fig. 10c, d, f), suggesting the slight alterations of pentacyclic terpanes, tricyclic terpanes and diasteranes from the PM7 to PM8 level, and sharp depletion from PM8 to PM9 + . However, the variation of diahopane concentrations (Fig. 10e) possibly indicated its high bioresistance or other mechanism needed to further investigate.

Table 4 Refined Manco numbers for Carboniferous oils of Chepaizi Uplift calculated by the refined Manco scores. P, pregnane; DH, diahopane
Fig. 10
figure 10

Correlations of the refined Manco number (MN2) with oil viscosity (a), NSO content (b), summed pentacyclic terpane (c) and tricyclic terpane concentrations (d), and absolute concentrations of diahopane (e) and diasteranes (f). All the concentration data were cited after (Chang et al. 2018)

Generally, RMN2 ranged from 546.04 to 909.63, which approximately distinguished the oils at the level of PM6 to PM9 + (Table 4). Exceptionally, by contrast, the oils from well P60 exhibited higher oil density but lower biodegradation extent than that from well P685. Factually, there is not a simple relationship between the oil viscosity and the Manco Number (Larter et al. 2012). Processes other than biodegradation, i.e., secondary oil charge, water washing, mixing of multiple maturity oil charges, and loss of light ends from heavy oils could produce variations in oil viscosity and API gravity (López 2014). The two episodes of oil charging, early biodegraded oils mixed with the later remigration of preexisting oils due to the structural adjustment, yet the same oil origin in the Chepaizi Uplift maybe responsible for this case (Chang et al. 2019; Shi et al., 2020).

Notably, although the RMN2 calculated by the refined Manco scale showed roughly positive correlation with oil density, some samples, especially those three biodegraded at PM8, exhibited much the same RMN2 but with different oil densities. This indicated that biodegradation refractory compound class added in this refined scale could not perfectly differentiate the appearances among oils biodegraded at about PM8. More assemblages of compound classes need to be established to refine the Manco scales in further trials.

4.4 Geochemical implications

Although the Carboniferous oils in the western and eastern Chepaizi Uplift were derived from different source kitchens, oils within the individual groups had common origin and similar maturities (Zhang et al. 2014a, b; Xi et al. 2014; Xiao et al. 2014; Xu et al. 2018; Mao et al. 2020). The differences observed in PM and refined Manco scales were mainly attributed to the biodegradation extent. By contrast, the refined Manco scale showed higher resolution than the PM scale, which can differentiate the biodegradation extent of super heavy oils with same PM values but different oil viscosities (Table 4).

For the eastern Chepaizi Uplift, Carboniferous oils were mostly biodegraded above PM8; however, their oil viscosity displayed orders of magnitude variations. Known to all, the physical properties, particularly the oil viscosity, were critical to the choice of exploration strategies. Subtle differences revealed by the refined Manco scale could be helpful. Besides the practical application, refined Manco scale also allow geochemists conducting more basic study to understand the relative differences in the extent of the biodegradation process among related samples.

5 Conclusions

The heavy to severe biodegradation was responsible for the Carboniferous heavy oils in the Chepaizi Uplift, with PM6 in the western and PM7–PM9+ in the eastern part, respectively. According the oil-source correlation, the oils in the western Chepaizi Uplift were derived from the Jurassic source rock in the Sikeshu sag and those in the eastern Chepaizi Uplift were originated from the Permian source rocks in the Changji sag. Therefore, two oil groups (I and II) were distinguished. Three oil families (II1, II2 and II3) were subdivided with the biodegradation level of PM7, PM8–8+ , PM9+ , respectively, based on molecular compositions and parameter-stripping method of strongly bioresistant parameters. Three biodegradation refractory compound class (25-norhopane, tricyclic terpanes and triaromatic steroids) were added to establish a refined Manco scale to quantify the biodegradation extent of Carboniferous oils from Chepaizi Uplift. The evaluation results clearly differentiate the biodegradation extent of heavy oils with same PM values but different oil viscosities, indicating a high resolution and potential prospect.