Abstract
Chemometrics has been widely used to cope with the problems of oil-oil and oil-source correlations because of its unique advantages in the comprehensive consideration of multiple parameters and the classification of samples or variables. In this paper, three chemometric methods, especially multidimensional scaling, were used to revisit the genetic oil family and the relationship between the crude oil and the source rock, because the oil source in the Wushi Sag, a significant petroliferous sag in the Beibu Gulf Basin of South China Sea, is still controversial. Two genetic families of crude oils, namely group A and group B, have been identified based on chemometric results. Group A oils are characterized by relatively higher Pr/Ph ratios and a high abundance of C27 ααα 20R steranes and C30-methylsteranes than those of group B oils, suggesting that this group of oils was deposited under a more oxic condition with more contribution of algae organic matter. Group A oils have been interpreted to be a mixture derived from the member 2 and member 3 of the Liushagang Formation (LS-2 and LS-3), whereas group B oils can be ascribed to the LS-2 member. The contribution of LS-3 mudstone member to the Wushi oils in previous studies may have been underestimated to some extent, which was inferred from the chemometric oil-source correlation results. The results of oil-source rock correlation may be used to guide future petroleum exploration activities with the incorporation of geological evidence. The spatial distribution of oil and gas reservoirs varies with burial depth. Taking into consideration other geological evidence, we may infer that the reservoir in eastern Wushi Sag was mainly distributed in the deep of Liushagang Formation, whereas the reservoir in southeast Wushi Sag was likely located in the shallow Liushagang Formation.
Introduction
Chemometrics is a useful tool for recognizing patterns and extracting useful information from measure data (; ). The method is use of multivariate statistics to identify and remove noise from the data, and show affinities among samples or variables (). Chemometric methods have been used to evaluate data in many disciplines (; ; ), including multivariate geochemical data (Zhan et al., 2019; Wang et al., 2020a). The application of chemometric method in geochemical correlation has been paid more and more attention due to the complexity of petroleum exploration and geological background. Therefore, many chemometric methods have been used to cope with the problem of oil-gas sources in petroleum systems according to practical demands (; Wang et al., 2019; Wang et al., 2020b). Among them, hierarchical clustering analysis (HCA) and principal component analysis (PCA) are the two commonly used methods (; ). Additionally, the method for oil-oil and oil-source rock correlations also includes factor analysis (Zumberge, 1987; ), star diagram (; ), K-nearest neighbor (; ), multidimensional scaling (; Wang et al., 2020a), discriminant analysis (Zhang et al., 2019; ), t-distributed stochastic neighbor embedding (). Meanwhile, geochemists try to introduce new methods (), indicating that chemometrics has great potential in coping with the problem of geochemical correlation. Multidimensional scaling (MDS) is a newly introduced chemometric method (), which is a nonlinear dimensionality reduction technique (), especially suitable for studying oil-gas sources in complex areas (; Wang et al., 2020b). Therefore, this paper intends to use the chemometrics method (e.g., MDS) to revisit the oil source in the Wushi Sag of the Beibu Gulf Basin of the South China Sea.
The Wushi sags of the Beibu Gulf Basin were confirmed to be hydrocarbon-rich sag by drilling data (; ). As of 2009, six commercial oil wells have been drilled in the Wushi Sag (Yang et al., 2011). Multiple studies have confirmed that three sets of source rocks exist in the Wushi Sag, namely, the LS-1 member (the first member of the Liushagang Formation), LS-2 member (the second member of the Liushagang Formation), and LS-3 member (the third member of the Liushagang Formation) (Yang et al., 2016; ; ). Controversies still remain about crude oil sources and hydrocarbon potential for each member of the Liushagang Formation in the Wushi Sag, although previous research has confirmed that the discovered oils in Beibu Gulf Basin may have been mainly derived from the second member of the Liushagang Formation (LS-2) (Zhao et al., 2013; ). Some investigators suggested that the Wushi oils are divided into three groups, and were primarily derived from the LS-2 member (Yang et al., 2016), whereas others believed that there are two categories of crude oils (I and II), in which group I oils were generated from the LS-2 member and group II oils representing mixtures generated from the second and third members of the Liushagang Formation (). Therefore, this paper intends to revisit the oil source of the Wushi oils using hierarchical clustering analysis (HCA), principal component analysis (PCA), and multidimensional scaling (MDS) of multivariate biomarkers. On the one hand, oil source rock correlation between oils and source rocks usually become more reliable when a composite analysis of multiple biomarker parameters are compared (; ). On the other hand, the result will help better understand the petroleum system in the Wushi Sag and reduce the exploration risk in this area. HCA calculates the distance matrices of the studied data objects, and organize objects with similar features into clusters (). Both PCA and MDS transform multiple correlated variables (i.e., biomarker ratios) into a small number of new uncorrelated variables (principal component). PCA is based on linear projection, whereas MDS is via non-linear dimensionality transformation ().
Geological Background
The Wushi Sag, as an important petroliferous sag, is located in the south of the Beibu Gulf Basin (Figure 1A). The sag is bounded by the Weixinan Sag in the north, the Haitoubei Sag in the west, the Maichen Sag in the south, and the Leidong Sag in the east (Figure 1A), and the wells distribution in the study area are shown in Figure 1B. It is approximately 2,680 km2 in area. The sag has undergone three stages of tectonic evolution, including the initial fault-depression stage, intense fault-depression stage, and late depression stage (; Yang et al., 2009; Yuan et al., 2020).
FIGURE 1
The bedrock in the Wushi Sag is a pre-Paleogene metamorphic rock, with a sedimentary cover of Cenozoic sediment. The Paleogene Changliu Formation, Eocene Liushagang Formation, Oligocene Weizhou Formation, Miocene Xiayang, Jiaowei, and Dengloujiao Formations, and Pliocene Wangliugang Formation were subsequently developed in the Cenozoic deposit from bottom to top (Figure 2). The Liushagang Formation was developed during the maximum expansion of the lake basin, in a lake system of alternating expansion and contraction and forming predominantly lacustrine and deltaic deposits (
FIGURE 2

General stratigraphic column of the Wushi Sag in the Beibu Gulf Basin (modified after (
The LS-3 member was deposited during lake expansion. The lithology of this member is interbedded sandstones and shales (
Samples and Methods
Samples
Five core samples of mudstone have been collected from potential Eocene source rocks, including the LS-1, LS-2, and LS-3 members from the WZ22-3-1, WZ10-7-1, and WZ10-3-30 wells in the study area (Table 1). In addition, four oil samples collected from the WS22-9-3D, WS16-1-HD, and WS16-1W-7 wells have been analyzed using gas chromatography-mass spectrometry. Additionally, twenty-two published samples, including source rocks and oils, have been used for chemometric analysis (
TABLE 1
| Well | Depth (m) | Formation | Litholoy | TOC (%) | S1 | S2 | S1+S2 | HI | OI | Tmax |
|---|---|---|---|---|---|---|---|---|---|---|
| WZ22-3-1 | 2,454 | LS-1 | Mudstone | 0.72 | 0.01 | 0.51 | 0.52 | 91 | 73 | 433 |
| WZ23-3-1 | 2,472.5–2,477.5 | LS-1 | Mudstone | 0.66 | 0.02 | 0.55 | 0.57 | 120 | 131 | 434 |
| WZ22-3-1 | 2,670–2,680 | LS-2 | Mudstone | 1.91 | 0.07 | 4.08 | 4.14 | 250 | 44 | 433 |
| WZ10-7-1 | 2,078–2,080 | LS-2 | Mudstone | 2.32 | 0.13 | 6.92 | 7.05 | 347 | 55 | 430 |
| WZ10-3-30 | 2,186–2,188 | LS-3 | Mudstone | 11.06 | 1.54 | 53.89 | 55.43 | 586 | 19 | 429 |
Rock-Eval data of core samples from the Liushagang Formation in the study area.
Note: The units of the Rock-Eval pyrolysis parameters and indices: TOC: wt%; Tmax: °C; S1: mg HC/g rock; S2: mg HC/g rock; S3: mg CO2/g rock; HI: mg HC/g TOC; OI: mg CO2/g TOC. LS-1, LS-2, and LS-3, represent the first, second, and third members of the Liushagang Formation.
TABLE 2
| Sample No. | Formation | Depth (m [ft]) | Lithology | No. | Group | Pr/Ph# | Pr/n-C17 | Ph/n-C18 | Ts/Tm# | Ol/C30H# | Ga/C30H # | C35/C34 | C27(%)# | C28(%)# | C29(%)# | C30-4MSI# | C29 20S/(20S + 20R) | C32 22S/(22S + 22R) | References |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| WZ111- 1R1 | LS-1 | 2,884–2,888 (9,462–9,476) | MD | 1 | — | 2.34 | — | — | 3.54 | 0.10 | 0.07 | — | 52 | 11 | 37 | 1.48 | 0.29 | — | [29] |
| WZ111- 1R2 | LS-1 | 2,928 (9,607) | MD | 2 | — | 2.19 | — | — | 4.25 | 0.12 | 0.05 | — | 56 | 9 | 35 | 0.93 | 0.58 | — | — |
| WZ122N- 1R1 | LS-1 | 2,835–2,840 (9,302–9,318) | MD | 3 | — | 2.88 | — | — | 1.57 | 0.06 | 0.07 | — | 48 | 8 | 44 | 1.68 | 0.28 | — | — |
| WZ22-3-1 | LS-1 | 2,454 | MD | 4 | — | 2.17 | 3.03 | 0.76 | 1.94 | 0.03 | 0.08 | 0.43 | 11 | 14 | 75 | 0.53 | 0.27 | 0.58 | This study |
| WZ23-3-1 | LS-1 | 2,472.5-2,477.5 | MD | 5 | — | 0.44 | 0.88 | 0.67 | 2.04 | 0.02 | 0.10 | 0.55 | 28 | 14 | 59 | 0.40 | 0.17 | 0.47 | This study |
| WZ122- 2R1 | LS-2 | 2,480 (8,137) | MD | 6 | — | 3.89 | — | — | 1.3 | 0.06 | 0.06 | — | 39 | 15 | 47 | 1.07 | 0.37 | — | [29] |
| WZ122- 2R2 | LS-2 | 2,798–2,800 (9,180–9,187) | MD | 7 | — | 1.64 | — | — | 1.64 | 0.07 | 0.05 | — | 34 | 14 | 51 | 3.75 | 0.51 | — | — |
| WZ122-2R3 | LS-2 | 2,788–2,790 (9,147–9,154) | MD | 8 | — | 1.56 | — | — | 1.43 | 0.06 | 0.05 | — | 35 | 13 | 52 | 3.88 | 0.53 | — | — |
| WZ122-2R4 | LS-2 | 2,764–2,766 (9,069–9,075) | MD | 9 | — | 1.83 | — | — | 1.28 | 0.05 | 0.05 | — | 37 | 15 | 48 | 3.12 | 0.39 | — | — |
| WZ1211-7R1 | LS-2 | 2,538–2,540 (8,327–8,334) | MD | 10 | — | 3.52 | — | — | 1.21 | 0.04 | 0.02 | — | 53 | 13 | 34 | 5.19 | 0.68 | — | — |
| WZ1211-7R2 | LS-2 | 2,558–2,560 (8,393–8,399) | MD | 11 | — | 2.29 | — | — | 0.96 | 0.06 | 0.03 | — | 43 | 15 | 42 | 7.64 | 0.66 | — | — |
| WS172-2R1 | LS-2 | 2,140–2,145 (7,021–7,038) | MD | 12 | — | 1.66 | — | — | 0.41 | 0.02 | 0.02 | — | 42 | 11 | 47 | 1.48 | 0.14 | — | — |
| WS172-2R2 | LS-2 | 2,280–2,285 (7,481–7,497) | MD | 13 | — | 3.06 | — | — | 0.46 | 0.02 | 0.03 | — | 37 | 14 | 49 | 4.32 | 0.25 | — | — |
| WS172-2R3 | LS-2 | 2,250 (7,382) | MD | 14 | — | 2.2 | — | — | 0.52 | 0.02 | 0.03 | — | 42 | 14 | 44 | 3.15 | 0.25 | — | — |
| WZ22-3-1 | LS-2 | 2,670-2,680 | MD | 15 | — | 2.68 | 2.78 | 0.76 | 2.39 | 0.08 | 0.07 | 0.53 | 29 | 17 | 54 | 0.66 | 0.20 | 0.44 | This study |
| WZ10-7-1 | LS-2 | 2,078-2,080 | MD | 16 | — | 0.94 | 0.99 | 0.85 | 2.48 | 0.03 | 0.05 | 0.84 | 35 | 16 | 49 | 0.94 | 0.23 | 0.51 | This study |
| WS172-2R4 | LS-3 | 2,518–2,520 (8,262–8,268) | MD | 17 | — | 2.41 | — | — | 1.09 | 0.06 | 0.06 | — | 38 | 14 | 48 | 1.04 | 0.21 | — | [29] |
| WZ122-2R5 | LS-3 | 2,832–2,834 (9,292–9,298) | MD | 18 | — | 1.77 | — | — | 1.37 | 0.08 | 0.07 | — | 34 | 14 | 51 | 2.61 | 0.4 | — | — |
| WZ118-1R1 | LS-3 | 3,392–3,394 (11,129–11,136) | MD | 19 | — | 1.76 | — | — | 2.19 | 0.09 | 0.06 | — | 36 | 19 | 45 | 1.18 | 0.52 | — | — |
| WZ118-1R2 | LS-3 | 3,312 (10,867) | MD | 20 | — | 2.04 | — | — | 1.96 | 0.10 | 0.08 | — | 37 | 22 | 41 | 1.21 | 0.56 | — | — |
| WZ10-3-30 | LS-3 | 2,186-2,188 | MD | 21 | — | 1.39 | 1.08 | 0.94 | 5.19 | 0.14 | 0.04 | 0.67 | 32 | 17 | 50 | 3.03 | 0.30 | 0.48 | This study |
| WS15-1M | XY | 1,230 (4,036) | CO | 22 | A | 2.12 | — | — | 1.2 | 0.06 | 0.03 | — | 39 | 15 | 47 | 2.33 | 0.47 | — | [29] |
| WS172- 7M | LS-2 | 1,615 (5,299) | CO | 23 | A | 2.39 | — | — | 1.33 | 0.05 | 0.03 | — | 31 | 12 | 56 | 2.87 | 0.41 | — | — |
| WS162-1DM | LS-3 | 1,308 (4,292) | CO | 24 | A | 2.54 | — | — | 1.41 | 0.14 | 0.10 | — | 40 | 11 | 49 | 2.86 | 0.44 | — | — |
| WS172-1DM | LS-3 | 2,627.2 (8,620) | CO | 25 | A | 2.19 | — | — | 1.59 | 0.07 | 0.04 | — | 41 | 14 | 45 | 2.67 | 0.41 | — | — |
| WS172- 1D2 | LS-3 | 2,619–2,646 (8,593–8,682) | CO | 26 | A | 2.27 | — | — | 2.05 | 0.09 | 0.05 | — | 37 | 16 | 47 | 2.89 | 0.45 | — | — |
| WS172-1D1 | LS-3 | 2,634–2,646 (8,642–8,682) | CO | 27 | A | 2.15 | — | — | 2.18 | 0.09 | 0.05 | — | 33 | 20 | 48 | 2.9 | 0.49 | — | — |
| WS229-3D | LS-1 | 2,990 | CO | 28 | B | 1.46 | 0.50 | 0.33 | 1.37 | 0.09 | 0.13 | 0.50 | 37 | 16 | 46 | 0.79 | 0.64 | 0.59 | This study |
| WS229-3D | LS-1 | 3,463.8 | CO | 29 | B | 1.37 | 0.41 | 0.30 | 1.25 | 0.13 | 0.11 | 0.45 | 32 | 15 | 53 | 0.72 | 0.59 | 0.59 | — |
| WS161-HD | LS-1 | 1,954 | CO | 30 | B | 1.18 | 0.49 | 0.42 | 1.68 | 0.11 | 0.14 | 0.45 | 31 | 13 | 56 | 0.50 | 0.53 | 0.60 | — |
| WS161W-7 | LS-2 | 2,658.5 | CO | 31 | B | 1.17 | 0.45 | 0.39 | 1.58 | 0.12 | 0.14 | 0.53 | 35 | 16 | 49 | 0.50 | 0.57 | 0.62 | — |
Selected geochemical parameters for the source rocks and crude oils in the Wushi Sag.
Note: LS-1, LS-2, LS-3, represent the first, second, and third member of the Liushagang Formation. “#” biomarker ratios for chemometric analysis; MD, Mudstone; CO, Crude oil; R1 = Pr/Ph; R2 = Pr/n-C17; R3 = Ph/n-C18; R4 = Ts/Tm; R5 = Ol/C30H; R6 = Ga/C30H; R7 = C35/C34 homohopanes; R8 = C27(%); R9 = C28(%); R10 = C29(%); R11 = C30-4MSI/C29 regular steranes; R12 = C29 steranes 20S/(20S + 20R); R13 = C32 22S/(22S + 22R) homohopanes. Pr/Ph = Pristane/Phytane; C27 (%), C28 (%), C29 (%) = relative percentage of C27, C28, and C29 ααα 20R steranes within the C27–C29 steranes (20R); 4MSI, C30-4-methylsterane index (4-methylsteranes/C29 ααα steranes); Ol/C30H = oleanane/C30 hopane; Ga/H = gammacerane/C30 hopane; Ts/Tm = 18a(H),21b(H)-22,29,30-trinorhopane/17a(H),21b(H)-22,29,30-trinorhopane; C35/C34 = C35/C34 homohopanes; C29 20S/(20S + 20R) = C29 steranes 20S/(20S + 20R); C32 22S/(22S + 22R) = C32 22S/(22S + 22R) homohopanes. These parameters can be referred to
Rock-Eval Pyrolysis and Bitumen Extraction
The core samples were cleaned first with distilled water and then crushed into powder. The powdered samples were subjected to IFP Rock-Eval 6 analyzer to obtain geochemical parameters such as total organic matter (TOC), hydrogen index (HI) and hydrocarbon generation potential (S1 + S2), following the same protocols described previously (
Gas Chromatography-Mass Spectrometry
Gas chromatography-mass spectrometry (GC-MS) analysis of oils and source rock extracts were analyzed using a Thermo Fisher Trace 1,300 gas chromatography coupled to an ISQ 7000 mass spectrometer equipped with a DB-5 capillary column (30 m × 0.25 mm × 0.25 μm). Pure helium gas (99.999%) was used as the carrier gas. The GC oven temperature was initially held at 50°C for 2 min, programmed to 120°C at a rate of 20°C/min, then to 250°C at 4°C/min and to 310°C at 3°C/min, and held at 310°C for 30 min. The mass spectrometer ionization energy was 70 eV and the ion source temperature was 230°C. Full-scan and selected ion recording were used for analysis, with the full-scan detection range of 50–550 (m/z). The peak area was integrated by Xcalibur software (Thermo Scientific). The biomarker compounds were identified in the studied rock and oil samples in accordance with the relative retention time, mass spectral feature, peak sequence, and identification standards from prior studies (
Chemometric Analysis
In this paper, three chemometric methods (e.g., HCA, PCA and MDS) were applied to yield a detailed oil-oil and oil-source rock correlations in the Wushi Sag using eight biomarker ratios (labelled by “#” in Table 2), similar to previous studies (
Results and Discussion
Source Rock Characteristics
Three sets of source rocks are present in the Wushi Sag: the first, second, and third members of the Liushagang Formation, i.e., LS-1, LS-2, and LS-3 (
FIGURE 3

Crossplots of (A) Rock-Eval S1 + S2 vs. total organic carbon (TOC) and (B) hydrogen index vs. Tmax for the studied potential source rocks in the Liushagang Formation. LS-1 = first member of the Liushagang Formation; LS-2 = second member of the Liushagang Formation; LS-3 = third member of the Liushagang Formation. The criteria distinguish poor, fair, good, very good and excellent source rock are from (
The thermal maturity of the studied rock samples was evaluated based on Rock-Eval Tmax and biomarker parameters. Tmax values for the Wushi samples range from 429 to 434°C, which indicates immature to early mature organic matter (
The distribution characteristics of n-alkanes and isoprenoid hydrocarbons in the Liushagang Formation can be used to interpret sources and depositional conditions of organic matter. The distribution of n-alkanes for the studied mudstone samples typically ranges from n-C15 to n-C31 (Figure 4). All three samples have bimodal distribution, suggesting a mixture of one mature source and one immature source in different proportions. The pristane/phytane (Pr/Ph) ratio is commonly used as an index of the depositional environment (
FIGURE 4

Representative n-alkane (m/z 85), hopane (m/z 191), and sterane (m/z 217) mass chromatograms for the source rocks from the first (LS-1), second (LS-2) and third (LS-3) member of the Liushagang Formation. C27, C28 and C29 = C27, C28 and C29 ααα 20R steranes; C30-4MSI = C30 4-methylsteranes. The symbolic meanings of other compounds can be found in Table 2.
The representative sterane chromatograms (m/z 217) for studied rock samples are illustrated in Figure 4. The relative content of C27, C28, and C29 ααα 20R steranes of the studied samples are in the range of 11–35%, 14–17%, and 49–75%, respectively (Table 2). This high C29 steranes contents may imply a strong contribution of terrigenous organic matter input (
Terpane mass chromatograms (m/z 191) of the representative source rocks are characterized by a dominant C30 hopane with a higher abundance of pentacyclic terpanes relative to tricyclic terpanes (Figure 4). The C35/C34 ratios of the studied source rocks range from 0.43-0.84, probably suggesting anoxic to suboxic depositional conditions (
Crude Oil Geochemistry
A complete set of acyclic isoprenoids (e.g., Pr and Ph) and low molecular weight n-alkanes are present in the analyzed oils, probably indicating no distinct signs of biodegradation (
FIGURE 5

Mass chromatograms of m/z 85 (A), 191 (B), and 217 (C) for the representative Wushi oil. C27, C28 and C29 = C27, C28 and C29 ααα 20R steranes; C30-4MSI = C30 4-methylsteranes; LS-1 = first member of the Liushagang Formation; LS-2 = second member of the Liushagang Formation; LS-3 = third member of the Liushagang Formation. The symbolic meanings of other compounds can be found in Table 2.
The m/z 217 mass chromatogram for the representative oil is shown in Figure 5B. The relative contents of C27, C28, and C29 ααα 20R steranes for the oils range from 31–37%, 13–16%, and 46–56%, respectively (Table 2). These data may indicate a substantial contribution of terrigenous organic matter input (
The m/z 191 mass chromatogram for the representative oil is shown in Figure 5C. These studied oils have C35/C34 ratios ranging from 0.45 to 0.53, indicating that oxidation is mainly occurring on the source (
Oil-Oil and Oil-Source Rock Correlations
Hierarchical cluster analysis (HCA) and principal component analysis (PCA) are valuable tools to identify genetic relationships among crude oils and source rocks based on age- and source-related biomarker ratios (
FIGURE 6

Chemometric analysis of oil genetic families for the Wushi Sag of the Beibu Gulf Basin based on (A) hierarchical cluster analysis (HCA) and (B) principal component analysis (PCA). The first three components of PCA, including PC1, PC2, and PC3, account for 92% of the total variance in the original dataset.
Group A oils are primarily produced from the LS-3 reservoirs (Table 2). This group of oils have relatively high Pr/Ph ratios, ranging from 2.12-2.54, and low Ga/C30 hopane ratios of 0.03–0.09, indicating the source rocks of Group B oils were probably deposited in more oxic condition (Table 2). The relative abundance of C29 ααα (20R) steranes and C30-4MSI ratios are in the range of 45–56% and 2.33–2.90%, respectively (Table 2). These data may indicate source rocks of Group B oils may have received more contribution of marine organic matter when compared to the source rocks of Group B oils. The ratios of C29 20S/(20S + 20R) range from 0.41 to 0.49 (Table 2), indicating that Group A oils are mature and have relatively lower maturity than that of Group B oils.
Group B oils principally consist of oils from the LS-1 member. The Pr/Ph ratios are relatively low, ranging from 1.16 to 1.46 (Table 2) and suggesting that the source rocks of Group A oils were deposited under more anoxic depositional environments. Group A oils have a relatively low abundance of C27 ααα (20R) steranes and a high abundance of C29 ααα (20R) steranes. Relatively low C30-4MSI ratios range from 0.50–0.79 (Table 2), probably indicating that the source rocks may have received more contributions from terrigenous organic matter. The maturity parameter of C29 20S/(20S + 20R) has relatively high values, ranging from 0.53–0.64 (Table 2), which is indicative of relatively high maturity.
Similarly, we have also identified genetic affinities between oils and source rocks using chemometric methods. Figures 7, 8 show the PCA and MDS scenario of multi-parameter oil-source rock correlation for the Wushi Sag. The biomarker parameters used for chemometric methods are the same as those in the HCA and PCA of Figure 6, marked by stars in Table 2. The three-dimensional view of PCA indicates two genetic oil families (Figure 7), inconsistent with the results discussed above. Group A oils have a strong affinity to the LS-2 and LS-3 members. The LS-2 mudstone member is related to group B oils. This conclusion is also confirmed by the MDS plot (Figure 8).
FIGURE 7

Chemometric results of oil-source rock correlation from the Wushi Sag of the Beibuwan Basin based on principal component analysis. The first three components of PCA, including PC1, PC2, and PC3, account for 92% of the total variance in the original dataset.
FIGURE 8

Chemometric analysis of oil-source rock correlation for the Wushi Sag of the Beibu Gulf Basin based on multidimensional scaling (MDS). The first two components of MDS, including MDS1 and MDS2, account for 96.5% of the total variance in the original dataset.
Implications for Petroleum Exploration
Oil and gas exploration activities have been primarily in oils derived from the second member of the Liushagang Formation (LS-2) and structure traps in the Wushi Sag over the last three decades (
As described previously, two categories of oils at different maturity levels have been identified using chemometric methods. Group A oils represent mixtures generated from the LS-2 and LS-3 members, whereas group B oils were derived from the LS-3 member. Group A oils with relatively lower maturation were distributed in eastern Wushi Sag (Figure 9) and primarily produced from the LS-3 reservoir of the WS162- and WS172- wells (Table 1). The corresponding maturity of Group A oils is consistent with the maturity of the nearby source rocks in the study area (Figure 9). The group of oil is close to the WS172 N Fault. The fault was in a shallower depth and had relatively weak activities during the Paleogene period (Yuan et al., 2020). Furthermore, the combination of oil-source rock correlation results and low permeability reservoirs (
FIGURE 9

The distribution of mudstone thickness and maturity of the third of the Liushagang Formation (LS-3 member) in the Wushi Sag (modified after (
Group B oils with relatively high maturation are distributed in southeast Wushi Sag and are mainly produced from the shallower LS-1 reservoir of the WS22-9 and WS16-1 wells (Table 1). The corresponding maturity of Group B oils is essentially equivalent to the maturity of the nearby source rocks from the LS-2 and LS-3 members (Figures 9, 10). Similarly, the distribution of oil is in the vicinity of the WS16-1 Fault, which holds across the entire source rock members of the study area (Figure 11) and is a good pathway connecting oil and source (Yuan et al., 2020). This evidence implies that group B oils were probably migrated vertically through faults, contributing to the shallower distribution of crude oil.
FIGURE 10

The distribution of mudstone thickness and maturity of the second of the Liushagang Formation (LS-2 member) in the Wushi Sag (modified after (
FIGURE 11

Geological profile of the Wushi Sag in south-north trending (modified after (Yuan et al., 2020)). Fm. = Formation, WZ = Weizhou Formation, LS = Liushagang Formation.
In summary, the results in the present study provide some inspiration for future petroleum exploration. For example, the oil-source rock correlation results suggest that the contribution of the LS-3 member to the Wushi oils has been underestimated in the past, and thus greater emphasis should be deserved in future exploration. In addition, exploration of the eastern Wushi Sag, where close to the distribution of Group A oils, should probably be focused on the relatively deep members of the Liushagang Formation, considering a combination of fault features and reservoir characteristics. In contrast, the exploration of southeast Wushi Sag, wherein the proximity of Group B oils should be paid more attention to the shallow layer of the Liushagang Formation.
Conclusions
The studied rocks from the three candidate source rock members, LS-1, LS-2, and LS-3, have a variable abundance of organic matter content and predominantly Type II kerogen. The LS-1 mudstone member of the studied sample is a fair source rock, whereas the LS-2 and LS-3 mudstone members of the studied sample have good to excellent hydrocarbon generation potential. Two oil families with different thermal maturity have been identified based on HCA and PCA. Group A oils with relatively low maturity are characterized by relatively high Pr/Ph ratios, low C27 ααα 20R steranes contents, and low C30-4MSI ratios. Conversely, the opposite is true for Group B oils. The oil-source rock correlation results show that group A oils were mainly derived from a mixture of LS-2 and LS-3 members, and group B oils were likely originated from the LS-2 mudstone member. These results also provide some inspiration for future petroleum exploration when combined with geological evidence. For example, the eastern Wushi Sag has a good potential for petroleum exploration in the deep Liushagang Formation (e.g., LS-3), whereas more attention should be given to the shallow Liushagang Formation in the southeast Wushi Sag.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
Y-PW conceived and wrote the manuscript; YG, XiZ, and ZS supervised and revise the manuscript; SW and JX provided technical support and conceptual advice; XuZ performed all the geochemical analyses and synthesized all the data. All authors discussed the results and implications and commented critically on important intellectual content on the manuscript at all stages.
Funding
This work was supported by the Doctoral Research Initiation Project of Guangdong Ocean University (Grant Nos R20030 and R17001), National Science Foundation of China (Grant No. 41602139), and the Special Financial Aid for Talents of Guangdong Ocean University (Grant No. 002026002004).
Acknowledgments
We acknowledge Yan-Rong Zou for his kind help with this study.
Conflict of interest
XiZ was employed by the company CNOOC International Limited.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AsemaniM.RabbaniA. R. (2021). A Novel and Efficient Chemometric Approach to Identifying Oil Families by Saturate Biomarker Data and FTIR Spectroscopy of Asphaltene Subfractions. Mar. Pet. Geology.124, 104838. 10.1016/j.marpetgeo.2020.104838
2
BrassellS. C.EglintonG.ShengG.FuJ. (1988). Biological Markers in Lacustrine Chinese Oil Shales. Geol. Soc. Lond. Spec. Publications40 (1), 299–308. 10.1144/gsl.sp.1988.040.01.24
3
ChakhmakhchevA.SuzukiN.SuzukiM.TakayamaK. (1996). Biomarker Distributions in Oils from the Akita and Niigata Basins, Japan. Chem. Geology.133 (1-4), 1–14. 10.1016/s0009-2541(96)00061-7
4
CheongK.-L.WuD.-T.DengY.LeongF.ZhaoJ.ZhangW.-J.et al (2016). Qualitation and Quantification of Specific Polysaccharides from Panax Species Using GC-MS, Saccharide Mapping and HPSEC-RID-MALLS. Carbohydr. Polym.153, 47–54. 10.1016/j.carbpol.2016.07.077
5
DidykB. M.SimoneitB. R. T.BrassellS. C.EglintonG. (1978). Organic Geochemical Indicators of Palaeoenvironmental Conditions of Sedimentation. Nature272 (5650), 216–222. 10.1038/272216a0
6
DongM.-Y.WuH.-L.LongW.-J.WangT.YuR.-Q. (2021). Simultaneous and Rapid Screening and Determination of Twelve Azo Dyes Illegally Added into Food Products by Using Chemometrics-Assisted HPLC-DAD Strategy. Microchemical J.171, 106775. 10.1016/j.microc.2021.106775
7
FuJ.ShengG. Y.PengP. G.BrassellS. C.EglintonG.JiangJ. G. (1986). Peculiarities of Salt lake Sediments as Potential Source Rocks in China. Org. Geochem.10 (1-3), 119–126. 10.1016/0146-6380(86)90015-x
8
FuN.LiuJ. (2018). Hydrocarbon Generation and Accumulation Characteristics of Three Type Source Rock of Liu II Segment in Beibuwan Basin. Nat. Gas Geosci.29 (07), 932–941. (in Chinese with English abstract). 10.11764/j.issn.1672-1926.2018.06.002
9
GanH.WangH.ShiY.MaQ.PanZ. (2019). Geochemical Characteristics and Genetic Origin of Crude Oil in the Fushan Sag, Beibuwan Basin, South China Sea. Mar. Pet. Geology.112, 104114. 10.1016/j.marpetgeo.2019.104114
10
GanJ.HuC.LiX.JinQ. (2017). Genetic Types and Differential Distribution Characteristics of Crude Oil in Wushi Sag. China Offshore Oil and Gas29 (06), 23–31. (in Chinese with English abstract). 10.11935/j.issn.1673-1506.2017.06.003
11
GaoZ.YangX.HuC.WeiL.JiangZ.YangS.et al (2019). Characterizing the Pore Structure of Low Permeability Eocene Liushagang Formation Reservoir Rocks from Beibuwan Basin in Northern South China Sea. Mar. Pet. Geology.99, 107–121. 10.1016/j.marpetgeo.2018.10.005
12
GongZ.LiS. (1998). Basin Analysis and Hydrocarbon Accumulation in the Northern continental Margin of the South China Sea.
13
GoodwinN. S.MannA. L.PatienceR. L. (1988). Structure and Significance of C30 4-methyl Steranes in Lacustrine Shales and Oils. Org. Geochem.12 (5), 495–506. 10.1016/0146-6380(88)90159-3
14
HarrisN. B.FreemanK. H.PancostR. D.WhiteT. S.MitchellG. D. (2004). The Character and Origin of Lacustrine Source Rocks in the Lower Cretaceous Synrift Section, Congo Basin, West Africa. Bulletin88 (8), 1163–1184. 10.1306/02260403069
15
HeJ.WuW.ZhuY.ChenS.CuiS.GongX. (2010). Genetic Types, Migration and Accumulation of Oil & Gas and Exploration Proposal in North Marginal Basins, South China Sea. Nat. Gas Geosci.21 (01), 7–17. (in Chinese with English abstract). 10.11764/j.issn.1672-1926.2010.01.7
16
HuangB.TianH.WilkinsR. W. T.XiaoX.LiL. (2013). Geochemical Characteristics, Palaeoenvironment and Formation Model of Eocene Organic-Rich Shales in the Beibuwan Basin, South China Sea. Mar. Pet. Geology.48, 77–89. 10.1016/j.marpetgeo.2013.07.012
17
HuangB.ZhuW.TianH.JinQ.XiaoX.HuC. (2017). Characterization of Eocene Lacustrine Source Rocks and Their Oils in the Beibuwan Basin, Offshore South China Sea. Bulletin101 (9), 1395–1423. 10.1306/10171615161
18
HuangW.-Y.MeinscheinW. G. (1979). Sterols as Ecological Indicators. Geochimica et Cosmochimica Acta43(5),739–745. 10.1016/0016-7037(79)90257-6
19
HuntJ. M. (1995). Petroleum Geochemistry and Geology. USA: W.H. Freeman and Company.
20
JustwanH.DahlB.IsaksenG. H. (2006). Geochemical Characterisation and Genetic Origin of Oils and Condensates in the South Viking Graben, Norway. Mar. Pet. Geology.23 (2), 213–239. 10.1016/j.marpetgeo.2005.07.003
21
KumarN.BansalA.SarmaG. S.RawalR. K. (2014). Chemometrics Tools Used in Analytical Chemistry: An Overview. Talanta123, 186–199. 10.1016/j.talanta.2014.02.003
22
LiF.LuL.MaW.DiaoH.XuZ.BaiY.et al (2021). Geochemical Characteristics of Source Rocks in Liushagang Formation of Wushi Sag. J. China Univ. Petroleum( Edition Nat. Science)45 (01), 31–40. (in Chinese with English abstract). 10.3969/j.issn.1673-5005.2021.01.004
23
LiM. (1994). “Paleoclimate, and Paleoenvironment in the Petroliferous Region of the north continental Shelf Region of South China Sea,” in Tertiary In Petroliferous Region of China (VIII). Editors JiangZ.ZenL.LiM. (Beijing: Petroleum Industry Press), 80–90. (in Chinese).
24
LiuM. (2004). The Assumed Gas Accumulation Systems in Beibuwan basin. China Offshore Oil & Gas16, 93–97. (in Chinese with English abstract). 10.3969/j.issn.1673-1506.2004.02.004
25
LiuY.RenY.YangX.LiuH.XuX. (2018). Geochemical Signatures of the Source Rocks from the Liushagang Formation in the Wushi Depression, Beibuwan Basin. Sediment. Geology. Tethyan Geology.38 (01), 103–112. (in Chinese with English abstract). 10.3969/j.issn.1009-3850.2018.01.012
26
LiuY.ZhaoZ.LiJ.XuJ. (2008). Hydrocarbon Play Division and Major Control Analysis in Wushi Sag, Beibuwan Basin. China Offshore Oil and Gas20 (02), 82–86. (in Chinese with English abstract). 10.3969/j.issn.1673-1506.2008.02.003
27
MashhadiZ. S.RabbaniA. R. (2015). Organic Geochemistry of Crude Oils and Cretaceous Source Rocks in the Iranian Sector of the Persian Gulf: An Oil-Oil and Oil-Source Rock Correlation Study. Int. J. Coal Geology.146, 118–144. 10.1016/j.coal.2015.05.003
28
OliveiraL. G.AraújoK. C.BarretoM. C.BastosM. E. P. A.LemosS. G.FragosoW. D. (2021). Applications of Chemometrics in Oil Spill Studies. Microchemical J.166, 106216. 10.1016/j.microc.2021.106216
29
PanseriyaH. Z.GosaiH. B.ValaA. K.GavaliD. J.DaveB. P. (2021). Assessment of Surface Water of Gulf of Kachchh, West Coast of India: A Chemometric Approach. Mar. Pollut. Bull.170, 112589. 10.1016/j.marpolbul.2021.112589
30
PetersK. E.CassaM. R. (1994). “Applied Source Rock Geochemistry,” in The Petroleum System: From Source to Trap, 60. Editors MagoonL. B.DowW. G. (American Association of Petroleum Geologists Memoir), 93–120.
31
PetersK. E.MoldowanJ. M. (1991). Effects of Source, thermal Maturity, and Biodegradation on the Distribution and Isomerization of Homohopanes in Petroleum. Org. Geochem.17 (1), 47–61. 10.1016/0146-6380(91)90039-m
32
PetersK. E.MoldowanJ. M. (1993). The Biomarker Guide: Interpreting Molecular Fossils in Petroleum and Ancient Sediments. Englewood Cliffs, NJ (United States): Prentice-Hall.
33
PetersK. E.WaltersC. C.MoldowanJ. M. (2005). The Biomarker Guide: Biomarkers and Isotopes in Petroleum and Earth History. Cambridge, United State: Cambridge University Press.
34
PetersK. E.Scott RamosL.ZumbergeJ. E.ValinZ. C.ScoteseC. R.GautierD. L. (2007). Circum-Arctic Petroleum Systems Identified Using Decision-Tree Chemometrics. Bulletin91 (6), 877–913. 10.1306/12290606097
35
PetersK. E.HostettlerF. D.LorensonT. D.RosenbauerR. J. (2008). Families of Miocene Monterey Crude Oil, Seep, and Tarball Samples, Coastal California. Bulletin92 (9), 1131–1152. 10.1306/04180807113
36
PetersK. E.CoutrotD.NouvelleX.RamosL. S.RohrbackB. G.MagoonL. B.et al (2013). Chemometric Differentiation of Crude Oil Families in the San Joaquin Basin, California. Bulletin97 (1), 103–143. 10.1306/05231212018
37
SeifertW. K.Michael MoldowanJ. (1978). Applications of Steranes, Terpanes and Monoaromatics to the Maturation, Migration and Source of Crude Oils. Geochimica et Cosmochimica Acta42 (1), 77–95. 10.1016/0016-7037(78)90219-3
38
SeifertW. K.MoldowanJ. M. (1980). The Effect of thermal Stress on Source-Rock Quality as Measured by Hopane Stereochemistry. Phys. Chem. Earth12, 229–237. 10.1016/0079-1946(79)90107-1
39
SeifertW. K.MoldowanJ. M. (1986). Use of Biological Markers in Petroleum Exploration. Methods Geochem. Geophys.24, 261–290.
40
ShiB.ChangX.XuY.WangY.MaoL.WangY. (2020). Origin and Migration Pathway of Biodegraded Oils Pooled in Multiple-Reservoirs of the Chepaizi Uplift, Junggar Basin, NW China: Insights from Geochemical Characterization and Chemometrics Methods. Mar. Pet. Geology.122, 104655. 10.1016/j.marpetgeo.2020.104655
41
Sinninghe DamstéJ. S.KenigF.KoopmansM. P.KösterJ.SchoutenS.HayesJ.et al (1995). Evidence for Gammacerane as an Indicator of Water Column Stratification. Geochimica et Cosmochimica Acta59 (9), 1895–1900. 10.1016/0016-7037(95)00073-9
42
SumithraV. S.SubuS. (2015). A Review of Various Linear and Non Linear Dimensionality Reduction Techniques. Int. J. Comp. Sci. Inf. Tech.6 (3), 2354–2360.
43
TaoK.CaoJ.WangY.MiJ.MaW.ShiC. (2020). Chemometric Classification of Crude Oils in Complex Petroleum Systems Using T-Distributed Stochastic Neighbor Embedding Machine Learning Algorithm. Energy Fuels34 (5), 5884–5899. 10.1021/acs.energyfuels.0c01333
44
TissotB. P.WelteD. H. (1984). Petroleum Formation and Occurrence. New York: Springer-Verlag.
45
VolkmanJ. K.AlexanderR.KagiR. I.NobleR. A.WoodhouseC. W. (1983). A Geochemical Reconstruction of Oil Generation in the Barrow Sub-basin of Western Australia. Geochimica et Cosmochimica Acta47 (12), 2091–2105. 10.1016/0016-7037(83)90034-0
46
WangY.-P.ZhangF.ZouY.-R.ZhanZ.-W.PengP. a. (2016). Chemometrics Reveals Oil Sources in the Fangzheng Fault Depression, NE China. Org. Geochem.102, 1–13. 10.1016/j.orggeochem.2016.09.008
47
WangY.-P.ZhangF.ZouY.-R.SunJ.-N.LinX.-H.LiangT. (2018a). Oil Source and Charge in the Wuerxun Depression, Hailar Basin, Northeast China: A Chemometric Study. Mar. Pet. Geology.89 (3), 665–686. 10.1016/j.marpetgeo.2017.10.032
48
WangY.-P.ZhangF.ZouY.-R.ZhanZ.-W.CaiY. (2018b). Origin and Genetic Family of Huhehu Oil in the Hailar Basin, Northeast China. Acta Geochim37 (6), 820–841. 10.1007/s11631-018-0288-1
49
WangY.-P.ZhanX.ZouY.-R.RodriguezC. G. V.LiaoW. (2019). Chemometric Methods as a Tool to Reveal Genetic Types of Natural Gases - A Case Study from the Turpan-Hami Basin, Northwestern China. Pet. Sci. Tech.37 (3), 310–316. 10.1080/10916466.2018.1542444
50
WangY.-P.ZhanX.LuoT.GaoY.XiaJ.WangS.et al (2020a). Oil Chemometrics and Geochemical Correlation in the Weixinan Sag, Beibuwan Basin, South China Sea. Energy Exploration & Exploitation38, 2695–2710. 10.1177/0144598720950467
51
WangY.-P.ZhanX.GaoY.WangS.XiaJ.ZouY.-R. (2020b). Chemometric Differentiation of Natural Gas Types in the Northwestern Junggar Basin, NW China. Energy Exploration & Exploitation38 (6), 2128–2142. Unsp 0144598720913972. 10.1177/0144598720913972
52
YangH.ZhaoZ.LiJ.LiuY.XuJ.XiaoL. (2009). An Anlysis of Petrleum Geology and Potential Exploration Targets in Wushi Sag. China Offshore Oil and Gas21 (04), 227–231. (in Chinese with English abstract). 10.3969/j.issn.1673-1506.2009.04.003
53
YangH.LianJ.HuW. (2011). Structural Features and Impacts on Hydrocarbon Accumulation in Wushi Sag. J. Southwest Pet. Univ. (Science Tech. Edition)33 (03), 41–46+191. 192. 10.3863/j.issn.1674-5086.2011.03.006
54
YangX.HuL.JinQ. (2016). Analysis of Hydrocarbon Origin and its Accumulation Time in the Middle and East Zone of Wushi Sag. J. Northeast Pet. Univ.40 (05), 9–17. (in Chinese with English abstract). 10.3969/j.issn.2095-4107.2016.05.002
55
YuanL.KongL.ZengX.ShangJ.RaoZ.LiuD. (2020). The Major Faults and its Control on Hydrocarbon Distribution in the Wushi Sag. J. Shengli Coll. China Univ. Pet.34 (04), 5–10. (in Chinese). 10.3969/j.issn.1673-5935.2020.04.002
56
ZhanZ.-W.LinX.-H.ZouY.-R.LiZ.WangD.LiuC.et al (2019). Chemometric Differentiation of Crude Oil Families in the Southern Dongying Depression, Bohai Bay Basin, China. Org. Geochem.127, 37–49. 10.1016/j.orggeochem.2018.11.004
57
ZhangL.BaiG.ZhaoX.ZhouL.ZhouS.JiangW.et al (2019). Oil-source Correlation in the Slope of the Qikou Depression in the Bohai Bay Basin with Discriminant Analysis. Mar. Pet. Geology.109, 641–657. 10.1016/j.marpetgeo.2019.06.055
58
ZhaoZ.WuJ.LiC. (2013). Sub-sag Sorting and Petroleum Distribution in Beibuwan Basin. Pet. Geology. Exp.35 (03), 285–290. (in Chinese with English abstract). 10.11781/sysydz201303285
59
ZumbergeJ. E. (1987). Prediction of source rock characteristics based on terpane biomarkers in crude oils: A multivariate statistical approach. Geochimica et Cosmochimica Acta51 (6), 1625–1637. 10.1016/0016-7037(87)90343-7
Summary
Keywords
Oil-Oil Correlation, Oil-Source Rock Correlation, Chemometrics, Wushi Sag, Beibu Gulf Basin
Citation
Wang Y-P, Zhan X, Zhou X, Gao Y, Wang S, Xia J and Song Z (2022) The Geochemical Study of Oil-Oil and Oil-Source Rock Correlations in the Wushi Sag of the Beibu Gulf Basin, South China Sea. Front. Earth Sci. 10:850922. doi: 10.3389/feart.2022.850922
Received
08 January 2022
Accepted
10 February 2022
Published
04 March 2022
Volume
10 - 2022
Edited by
Shuang Yu, Guangzhou Institute of Geochemistry (CAS), China
Reviewed by
Taotao Cao, Hunan University of Science and Technology, China
Lu Hong, Guangzhou Institute of Geochemistry (CAS), China
Qingtao Wang, Guangzhou Institute of Energy Testing, China
Updates

Check for updates
Copyright
© 2022 Wang, Zhan, Zhou, Gao, Wang, Xia and Song.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yao-Ping Wang, wangyp@gdou.edu.cn
This article was submitted to Geochemistry, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.