ORIGINAL RESEARCH article

Front. Earth Sci., 19 September 2025

Sec. Geochemistry

Volume 13 - 2025 | https://doi.org/10.3389/feart.2025.1650751

Occurrence and recoverability of shale oil in Paleogene formation, Subei Basin

  • Wuxi Petroleum Geology Research Institute, SINOPEC Exploration & Production Research Institute, Wuxi, Jiangsu, China

Abstract

In recent years, significant breakthroughs have been made in the exploration and development of lacustrine shale oil, but the stable and guaranteed production of shale oil still faces huge challenges. Quantitative characterization of the occurrence state of shale oil and accurate evaluation of the movable oil content are the key scientific issues that need to be solved for the efficient exploration and development of lacustrine shale oil. The second and fourth sections of the Funing Formation in the northern Jiangsu Basin are preferred layers for the exploration and development of continental shale oil in eastern China. However, there are few reports on the occurrence status and availability of shale oil in these areas. This study systematically investigates the oil-bearing characteristics and mobility mechanisms of shale and interlayer systems in the second and fourth membersof the Funing Formation (E1f2, E1f4), Subei Basin, through an integrated approach combining multi-temperature pyrolysis, nuclear magnetic resonance (NMR) analysis, and limited production data from 23 wells. Key findings reveal:1) Shale matrices exhibit low total retained oil content (<5.0 mg/g), dominated by adsorption-miscible oil (40%–95% of total hydrocarbons), with its proportion inversely correlated to burial depth (R2 = 0.78). Light free oil content remains exceptionally low (<0.1 mg/g), yielding a mobility ratio <3%. In contrast, interlayers and adjacent sandstones demonstrate significantly higher oil saturation (>5 mg/g), predominantly as free-phase oil (60%–82%), with light free oil exceeding 0.25 mg/g and mobility ratios of 4–7%—comparable to pre-stimulation Bakken tight sandstones (5%); 2) NMR-derived T2 cutoff values (6 m) and throat radius thresholds (>0.18 μm) indicate movable shale oil primarily resides in fracture networks (contributing 68%–85% of permeability) and secondary micropores. Phase-state analysis reveals free oil > adsorbed oil within these conductive pathways, corroborated by production data from fractured carbonate reservoirs.3) Fractured lacustrine carbonates exhibit high initial yields (e.g., Well XuX38: 15.8 t/d) but rapid decline rates (68% within 90 days). Conversely, intra-source interlayer systems demonstrate sustained production (Well TianX96: <15% decline over 18 months), attributed to stable pressure maintenance and effective matrix-fracture connectivity. Economic analysis suggests 2.1× higher EUR (Estimated Ultimate Recovery) in interlayer systems compared to fractured carbonates.

1 Introduction

Shale oil, as an important component of unconventional oil and gas resources, has become a new focus in unconventional exploration in recent years (; ; ; Zou et al., 2015; ; ; ; ). In 2022, the United States emerged as the resource-producing country with the highest growth in global crude oil production, with tight oil—primarily shale oil—contributing the most, increasing by 275.0786 million tons year-on-year (). In China, continental shale oil has achieved industrial breakthroughs in multiple basins, including the Junggar, Ordos, Bohai Bay, Songliao, Sichuan, and Subei basins (; ; ; ; ; ; ; ; ; ).

The second and fourth members of the Funing Formation in the Subei Basin are the principal source rocks for conventional oil and gas, and they also constitute the primary target intervals for the formation of unconventional shale oil reservoirs (; ). According to several scholars and experts, the shale oil resources within the Funing Formation could be the most promising replacement for conventional resources in the Subei Basin (Zheng and Peng, 2013; ). With ongoing research and exploration into shale oil, earlier studies have made significant advances in understanding the formation conditions, reservoir characteristics, fracture prediction, and enrichment patterns of these shale systems (; ; ; ). However, there is still a relative lack of research regarding the occurrence states and mobility of shale oil within the Funing Formation. Accurately evaluating how shale oil exists within reservoirs—its forms, states, and recoverability—is currently a key bottleneck in shale oil exploration and one of the urgent issues that needs to be resolved (; ). This paper is based on a systematic analysis of shale, interlayers, and adjacent carbonate and sandstone units from the second and fourth members of the Funing Formation. By employing multi-step pyrolysis and nuclear magnetic resonance (NMR) techniques, the study aims to investigate and evaluate the oil-bearing properties, occurrence states, and mobility of shale oil in both intra-source interlayers and adjacent rocks (including the first and third members of the formation), with the goal of providing insights into the exploration and development of shale oil in the Subei Basin.

2 Regional geological overview

The Subei (Northern Jiangsu) Late Cretaceous-Cenozoic Basin overlies the Mesozoic basins of the Lower Yangtze region, with its internal structural units controlled by the Indosinian-Yanshanian tectonic framework (; ; ). Based on the relationship between the Subei Basin and surrounding regional structural units, its interior can be divided into four second-order structural units: the Dongtai Depression, Jianhu Uplift, Yanfu Depression, and Binhai Uplift. Each depression can be further subdivided into several sub-level sags and low uplifts (Figure 1a). Among them, the four major sags within the Dongtai Depression—Jinhu, Gaoyou, Qintong, and Hai’an—are the key exploration areas for oil and gas discoveries.

FIGURE 1

Although the stratigraphy of the Funing Formation in the Subei Basin exhibits notable lithological and electrical property variations across different depressions, a general vertical rhythmic pattern of sandstone–mudstone–sandstone–mudstone is observed from the first to the fourth members. Member E1f1 and E1f3 are characterized by the development of large deltaic depositional systems, featuring widespread sandstone reservoirs, which predominantly consist of single sets of sandstone beds. Member E1f2 and E1f4, in contrast, are dominated by extensive lacustrine depositional systems, hosting widely distributed source rocks primarily composed of dark gray mudstone and shale. These members also contain thin carbonate and tight sandstone interlayers embedded within the shale source rocks. Overall, these lithological features provide the material foundation necessary for the formation of shale oil. Current research indicates that the Funing Formation in the Subei Basin comprises:Four regionally continuous and mature source rock intervals (E1f21, E1f22, E1f23, E1f41), Three tight sandstone reservoirs closely associated with source rocks (E1f11, E1f23 and E1f33), Three intra-source carbonate or sandstone interlayers (E1f11, E1f23 and E1f33), As well as two sets of regional interlayers (E1f21 caprock and E1f42), which serve as effective seals.

Considering that the source rocks of E1f42 and E1f23 (in the Gaoyou and Jinhu Sag) have low organic richness and poor hydrocarbon generation potential, and due to their considerable thickness, they act as vertical barriers to the overlying oil-prone layers and underlying tight reservoirs. Taking into account the lateral migration capacity of faults and sand bodies, five favorable shale oil accumulation models can be identified (Table 1; Figure 1b), including:1. Source-below model: where source rocks (E1f22) underlie tight sandstone reservoirs (E1f11), 2. Intra-source model: formed by source rocks (E1f21, E1f22, E1f23) and their internal carbonate or sandstone interbeds, 3. Distal-source fault–sandstone-guided model: involving source rocks in E1f41 and E1f2 and tight sandstone reservoirs in E1f2, E1f3, and E1f1 connected through faults and lateral migration pathways.

TABLE 1

StratumReservoir layerOil generation layerOil and gas geological propertiesDevelop layers, regions and lithology
MemberSubmemberThickness/mΦ
%
K mDThickness/mTOC
%
Sub-sourceIn-sourceFar-source
E1f4E1f41>1001.2⑤ Far-source fault-sand conduit hydrocarbon accumulationDense sandstone of Gaoyou Sag
E1f42>1000.74
E1f3E1f31<308–120.1–1
E1f32
E1f33<20
Laminar carbonate rocksBiological limestone and sandstone interlayers in the Jinhu and Gaoyou Sag
E1f2E1f215–103.240–802.14② In-source
E1f225–108–262.2–79.860–801.29① Sub-source③ In-source
E1f235–100.98④ In-source
<408–120.1–1
E1f1E1f11<808–120.1–1Dense sandstone of Gaoyou Sag

Distribution of source - reservoir superimposed relationship and accumulation association in E1f, Subei Basin.

3 Samples and experimental methods

3.1 Sample collection

To investigate the occurrence states of shale oil in the Funing Formation under different lithological combinations, systematic sampling was conducted on shale, carbonate rock (tight sandstone) interbeds within source rocks, and adjacent tight sandstone layers in key sags (Gaoyou and Jinhu) of the study area. A total of 165 core samples and 213 cuttings samples were collected from 17 wells. Subsequent experimental analyses included multi-temperature step pyrolysis on 60 samples, organic carbon testing on 170 samples, and nuclear magnetic resonance (NMR) analysis on 174 samples.2.2 Experimental Test Methods.

3.1.1 Organic carbon

The total organic carbon content of the core samples was determined using the LECO CS-230 carbon analyzer of China University of Petroleum (Beijing). To remove inorganic carbon, the core samples were crushed and sieved through an 80-mesh sieve, and then reacted with a 10% (volume fraction) HCL solution for more than 2 h until the reaction was complete. The measurement accuracy is estimated to be 0.45% of the measured value.

3.1.2 Multi-temperature pyrolysis

This study obtained the occurrence state of shale oil in the samples through multi-temperature pyrolysis experiments. The experimental analysis was conducted using the improved Rock-Eval 6 pyrolysis instrument of Wuxi Petroleum Geology Research Institute of Sinopec Petroleum Exploration and Development Research Institute. S1-1 (light oil component) was measured at a constant temperature of 200 °C for 1 min. Then, it was heated to 350 °C at a rate of 25 °C/min and kept at a constant temperature for 1 min to measure S1-2 (light - medium oil component). Then, the temperature was raised to 450 °C at a rate of 25 °C/min and maintained at a constant temperature for 1 min to measure S2-1 (heavy hydrocarbon and gelatinous asphaltene components). Finally, the temperature was raised to 600 °C at the same rate to measure S2-1 (kerogen pyrolytic hydrocarbon).

3.1.3 NMR experiments

The two-dimensional nuclear magnetic resonance (NMR) experiments were conducted using an American Core Company MRCore-040V 23 MHz NMR instrument to quantitatively measure the longitudinal relaxation time (T1) and transverse relaxation time (T2) of 1H-containing compounds in shale samples under conventional processing conditions, with a total of 5 samples analyzed.

For the organic matter extraction experiments, dichloromethane (DCM) solvent was employed in a Soxhlet extraction apparatus to isolate solvent-soluble components from the samples. The extraction process was carried out at 70 °C for 48 h to quantify the retained oil within the shale samples. Following the extraction experiments, the two-dimensional NMR measurements were repeated.

4 Occurrence state of shale oil

4.1 Composition characteristics of shale oil

During the room-temperature storage of shale, a significant portion of the light hydrocarbon components (C6–C13) in the samples is lost. These lost light hydrocarbons constitute an important part of the free oil in shale oil. Directly using pyrolysis parameters to represent residual hydrocarbons often underestimates the actual free oil content. Therefore, light hydrocarbon correction is necessary to restore the lost light hydrocarbons before pyrolysis experiments, thereby accurately characterizing the occurrence state of shale oil.

A comparative analysis was conducted between the pyrolysis parameter S1 of freeze-preserved samples and the pyrolysis parameter S1’ of samples stored at room temperature to calculate the light hydrocarbon correction factor (K = S1/S1’) (). Due to the excessive number of data points, it was difficult to fit a curve. Instead, the samples were sorted by depth, and averages were taken every 10 m to establish a chart showing the variation of the light hydrocarbon correction factor for shale oil with depth (see Figure 2). The light hydrocarbon correction factor mainly ranges from 1.1 to 4.2, with an average of 2.6. As shown in Figure 2, the correction factor increases with depth, reflecting a trend of higher organic matter maturity and a significant increase in the proportion of light hydrocarbons as depth increases.

FIGURE 2

The total quantity of retained oil in various occurrence states within shale or sandstone serves as a critical indicator of oil-bearing potential. Two primary methods—solvent extraction and thermal desorption—are commonly used for quantitatively analyzing retained oils in free and adsorbed states (; ). In this study, systematic analyses were conducted using multi-step pyrolysis on samples from the Funing Formation, including shale, intra-source carbonate (tight sandstone) interbeds, and adjacent tight sandstones (Table 2). The results were used to evaluate the oil-bearing characteristics and oil mobility in shale, interbeds, and neighboring tight sandstones. Among them, the pyrolytic hydrocarbon S1-1 is primarily composed of light oil components existing in a free state, representing the most readily mobilizable oil currently available. S1-2 mainly consists of light-medium oil components and also occurs in a free state. S2-1 is dominated by heavy oil, resin-asphaltene components, coexisting with organic matter in an adsorbed-dissolved state. S2-2 primarily originates from regenerated hydrocarbons through kerogen pyrolysis in shale. Therefore, the total free hydrocarbon content is the sum of S1-1 and S1-2, while the total retained oil volume comprises both free hydrocarbons and adsorbed-dissolved oil, serving as critical parameters for characterizing oil-bearing properties.

TABLE 2

LithologyWell numberStratumSample positionFree light oil S′1–1 (mg/g)Free medium - heavy oil S′1–2 (mg/g)Adsorb miscible oils S′2–1 (mg/g)Total free oil content (mg/g)Total retained oil volume (mg/g)Free hydrocarbons S′1-2/Total retained oil (%)Adsorption - miscible oil/Total retained oil (%)Free hydrocarbons S′1-1/Total retained oil (%)Pyrolysis regenerated oil S′2-2 (mg/g)Sag
Light grey siltstoneCheng2E1f2interlayer0.072.461.022.533.5569.3028.731.970.50Jinhu
Grey silty mudstoneCheng2E1f2interlayer0.000.010.020.010.0333.3366.670.000.05
Light grey calcareous mudstoneCui2E1f3interlayer0.000.151.840.151.997.5492.460.002.73
Grey siltstoneDai1E1f2interlayer0.011.187.591.198.7813.4486.450.117.22
Dark grey marlHeX4E1f4Adjacent layer0.000.122.030.122.155.5894.420.004.32
Grey siltstoneHeX4E1f2interlayer0.416.113.256.529.7762.5433.274.201.55
Dark grey mudstoneHeX4E1f2Shale0.000.010.020.010.0333.3366.670.000.02
Grayish-black mudstoneHeCan1E1f2Shale0.020.981.291.002.2942.7956.330.872.52
Grayish-black mudstoneHeCan 1E1f2Shale0.061.333.071.394.4629.8268.831.354.96
Grayish-black mudstoneHeCan 1E1f2Shale0.000.010.040.010.0520.0080.000.000.08
Grey siltstoneTang5E1f2interlayer0.618.315.108.9214.0259.2736.384.353.07
Dark grey sandy mudstoneTian89E1f2interlayer0.000.190.350.190.5435.1964.810.000.39
Grey sandstoneTianX79E1f4interlayer0.516.292.976.809.7764.3830.405.222.51
Dark grey mudstoneTianX79E1f4Shale0.000.030.060.030.0933.3366.670.000.10
Grey mudstoneDaX3E1f2Mud shale0.000.655.520.656.1710.5389.470.007.99Gaoyou
Grey mudstoneFang4E1f3Adjacent layer0.213.567.583.7711.3531.3766.781.8513.90
Grayish-black mudstoneFuShen1E1f1Adjacent layer0.020.340.330.360.6949.2847.832.900.29
Grayish-black mudstoneFuShen 1E1f2Mud shale0.070.921.100.992.0944.0252.633.351.96
Grayish-black mudstoneFuShen 1E1f4Mud shale0.020.773.290.794.0818.8780.640.499.21
Dark grey mudstoneHuaShen 1E1f30.020.140.120.160.2850.0042.867.140.59
Grey siltstoneLian5E1f2interlayer0.000.020.010.020.0366.6733.330.000.01
Dark grey marlLian 5E1f3Adjacent layer0.000.050.230.050.2817.8682.140.000.56
Black mudstoneLian 5E1f4Mud shale0.000.020.090.020.1118.1881.820.000.16
Grayish-black mudstoneLinX3E1f2Mud shale0.000.040.180.040.2218.1881.820.000.24
Grayish-white sandstoneSha20E1f1Adjacent layer0.8614.495.5715.3520.9269.2626.634.112.73
Dark grey sandstoneSha20E1f1Adjacent layer0.7212.906.0213.6219.6465.6830.653.672.82
Dark grey siltstoneWei5E1f1Adjacent layer0.599.437.6910.0217.7153.2543.423.334.50

Evaluation of oily and mobility in shale, interlayer and adjacent sandstone in E1f2 and E1f4 in key sag, Subei Basin.

Based on the analysis of shale samples from the E2 and E4 Members in the Gaoyou and Jinhu Sags (the main exploration areas of the Subei Basin), along with their interlayers and adjacent tight sandstone layers (Table 2; Figures 3a, b), the total retained oil volume in the E2 and E4 shale of the Jinhu Sag is generally less than 5.0 mg/g, with some samples below 0.5 mg/g (Figure 3a). Adsorbed-dissolved oil dominates the retained oil, accounting for 40%–95% of the total. The proportion of adsorbed-dissolved oil exhibits a decreasing trend with increasing burial depth. Light free oil content remains extremely low (<0.1 mg/g), and the current mobilizable oil ratio is below 2%. In contrast, the total retained oil volume in interbedded shales and adjacent oil-bearing tight sandstones of the E2 and E4 Members typically exceeds 5 mg/g (maximum 14.02 mg/g), with free oil as the dominant phase. Light free oil content generally surpasses 0.25 mg/g, yielding a current mobilizable oil ratio of 4%–7%. This aligns closely with the ∼5% recovery rate observed in unstimulated Middle Bakken tight sandstone intervals in North America ().

FIGURE 3

In the E2 and E4 shale of the Gaoyou Sag, the total retained oil volume generally remains below 5.0 mg/g, with some samples registering less than 0.5 mg/g. Adsorbed-dissolved oil predominates, accounting for 50%–92% of the total retained oil. This proportion exhibitis a decreasing trend with increasing burial depth. Light free oil content is consistently low (<0.2 mg/g), resulting in a current mobilizable oil ratio below 3%. Conversely, interbedded shales and adjacent oil-bearing tight sandstones within the E2 and E4 Members typically contain retained oil volumes exceeding 15 mg/g (Figure 3b), dominated by free oil. Light free oil content in these intervals generally surpasses 0.50 mg/g, yielding a current mobilizable oil ratio of 3%–4%.

As shown in Figure 4 illustrating the relationships between free oil (S1-1, S1-2), adsorbed-dissolved oil (S2-1), and TOC, free oil primarily resides in inorganic pore-fracture systems (Figure 4a,b), whereas adsorbed-dissolved oil predominantly occurs within Figure 4c organic matter, demonstrating a moderate positive correlation with TOC (Figure 4c).

FIGURE 4

4.2 Shale oil saturation characteristics

Experience from North American shale oil exploration shows that productive intervals—whether mudstone or sandstone—typically exhibit an oil saturation index (S1 × 100/TOC) greater than 100 (). Statistical results from pyrolysis analyses of E1f2 and E1f4 samples from the main shale oil sags in Gaoyou and Jinhu indicate (Figure 5): The E1f2 section is generally moderately to highly oil-bearing, with some samples showing an oil saturation index >100, signifying shale oil potential (Figure 5a). The E1f4 section is mostly low in oil content, with only a few samples indicating shale oil potential (Figure 5b). The data suggest that the E1f2 mudstone–shale system has good exploration potential, and that the oil saturation index tends to increase with thermal maturity (Tmax).

FIGURE 5

5 Shale oil mobility

5.1 NMR characteristics of shale oil mobility

Mobility and movable fluid content are critical parameters for shale oil evaluation. Zhou et al. (2014) applied nuclear magnetic resonance (NMR) technology to analyze movable fluids in shale gas reservoirs, tight sandstone, and tight limestone reservoirs. Zhang et al. (2014) combined experimental analysis with well-logging data to investigate shale oil mobility through systematic studies of shale porosity, compressibility, rock mechanical properties, oil saturation, gas-oil ratio, and crude oil saturation pressure. This study primarily utilizes NMR data from 174 E1f4 shale samples of Well Huang 158 to analyze the mobility and occurrence state of shale oil in the E1f lacustrine carbonate-rich shales of the Subei Basin. The NMR results reveal that most shale samples exhibit unimodal T2 spectra (Figure 6), indicating no detectable movable fluids. Only 62 samples show bimodal or trimodal T2 spectra, suggesting the presence of limited movable fluids. The movable fluid content in these samples is generally <20.68% (average 6.11%; Table 3), with a T2 cutoff time of 6 m (corresponding to a pore-throat radius of 0.18 μm). These findings imply that:Pore-throat radii in E1f4 shales are predominantly <0.18 μm.Movable fluids (interpreted as crude oil based on fracturing and production testing) are only retained in micropores with radii >0.18 μm within a limited subset of shales.

FIGURE 6

TABLE 3

Crack developmentMovable fluid,%Φ,%
All samplesThere are movable fluid samplesAll samples
Developmental crack0.17 (17)0.14∼2.52/1.47 (2)0.24∼6.58/2.10 (17)
Multiple development cracks3.34 (89)0.01∼20.68/6.75 (44)0.26∼13.2/4.20 (89)
Thin-layered, 2 page-shaped0.27 (32)1.06∼3.80/2.15 (4)1.5∼10.18/6.00 (32)
Thin-layered, developed cracks0.81 (36)0.14∼11.58/4.86 (6)1.39∼10.36/6.20 (36)

Movable fluid content in shales in E1f4, well Huang158.

However, the linear relationship between the content of movable fluid in the E1f4 shale of Well Huang 158 and its matrix porosity (Figure 7a) and organic matter abundance (Figure 7b) is not obvious, reflecting that movable shale oil should mainly occur in fractures, that is, structural fractures are the main occurrence space of movable shale oil, and some micropores (with a roar radius greater than 0.18 μm) are secondary occurrence Spaces. This determines that the occurrence phase state of shale oil should mainly be the free state, followed by the adsorbed state.

FIGURE 7

5.2 Production test characteristics of shale oil mobility

This study takes well locations Xu X38 and Tian 96 as examples to elaborate on the characteristics of shale oil mobility test production in two different types of shale reservoirs: intra-source fractured type and intra-source interlayer type. The burial depth of the Xu X38 E1f2 shale oil reservoir ranges from 2580 m to 2840 m, and the oil-bearing layers are E1f2 page 2 (3,019.3m–3021.3 m) and E1f2 page 3 (3,076.3m–3077.6 m). According to the drilling data of Well Xu X38, the lithology of the upper oil-bearing layer is mainly dark gray massive grey mudstone, with underdeveloped margins but rich in biological debris. The lower oil-bearing layer is a stratified gray-bearing and cloud-bearing shale with well-developed margins and a thickness of 41 m. The average TOC of the upper oil-bearing layer is 2.80%. The average TOC of the lower oil-bearing layer is 1.39%, all of which have excellent oil generation capacity. Logging interpretation: The porosity of the reservoir ranged from 5.6% to 7.9%, with an average of 6.9%. Permeability: (0.2–6.6) ×10−3 μm2, with an average of 1.1 × 10-3 μm2; The oil saturation of the reservoir is 65.8%. In addition, the sensing resistivity of the oil-bearing layer is lower than that of the upper and lower surrounding rocks. The resistivity is in the shape of a “V” or a sharp knife, and the cracks are well developed. Well logging interpretation: There are multiple fracture development sections in Well X38 E1f2. During the fracture development sections, the oil and gas are active, with the highest peak of total hydrocarbons reaching 65.8% and 99.9% respectively. It belongs to the intra-source fracture type shale oil reservoir.

The oil-bearing strata of the Tian96 E1f2 shale oil reservoir include the sandstone interlayer within the E1f21 source and the sand body of the Xiatan dam in the E1f22 source. In the Tongcheng area where the oil reservoir is located, E1f21 is a subfacies deposition of shallow lakeside. Within the thick dark mudstone layer, thin-layer beach dam sand bodies are developed. The percentage content of sandstone is about 20%. The maximum single-layer thickness is 6.1 m and the minimum is 2 m. The rock types are feldspar siltstone and fine sandstone. The porosity ranges from 4.61% to 13.25%, with an average of 9.30%. The permeability was (0.10–7.37) ×10−3 μm2, with an average of 3.51 × 10-3 μm2. The oil saturation ranged from 14.28% to 32.19%, with an average of 23.24%. The lithology of the source rocks is mainly grayish-black striated marl, with some areas being striated to thin marl. Fossils of ostracods can be found in the striated marl, which is 19 m thick. The average TOC of shale was 2.51%, the average “A” of chloroform asphalt was 0.26%, and the average S1 of pyrolyzed hydrocarbons was 0.44 mg/g. It belongs to the intra-source interlayer type shale oil reservoir.

From the perspective of the later oil test and production results of shale oil, although the initial output of fractured lacustrine carbonate shale oil is relatively high, the output decreases rapidly (Figure 8a), the production cycle is short, and the economic recoverability is limited. In contrast, although the initial output of the intra-source interlayer shale oil is not high, it decreases slowly (Figure 8b), has a long production cycle, and has better economic recoverability.

FIGURE 8

6 Conclusion

  • 1. The occurrence state and mobility of shale oil in the Funing Formation of the Subei Basin are jointly controlled by organic matter abundance, type, and porosity-permeability characteristics. Research indicates a significant positive correlation between organic matter abundance (TOC) and adsorbed-miscible oil, suggesting that high-TOC intervals are more conducive to oil adsorption and enrichment. In contrast, free oil primarily occurs in inorganic pore-fracture systems and shows a weaker relationship with TOC. Increasing organic matter maturity promotes the generation of light free oil, but the overall movable oil ratio remains relatively low. Nuclear magnetic resonance (NMR) analysis reveals that micropores with throat radii >0.18 μm serve as secondary storage spaces for movable fluids, while fracture systems act as the primary flow channels.

  • 2. In the E1f2 and E1f4 mudstone–shale members, retained oil is dominated by adsorbed–miscible oil with low mobility, suggesting limited shale oil exploration potential. However, interbeds (carbonate or sandstone) and adjacent tight sandstones are dominated by free oil, have relatively higher movable oil content, and thus possess greater exploration potential.

  • 3. In contrast, the initial production of fractured shale oil reservoirs within the source is high, but the production decreases rapidly, the production cycle is short, and the economic recoverability is limited. The initial production of tight sandstone oil reservoirs in interlayers or adjacent layers within the source is relatively low, but the decline is slow and the production cycle is long. They have good economic recoverability, providing a basis for the selection of favorable exploration areas for shale oil in the northern Jiangsu Basin.

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

PJ-N: Project administration, Methodology, Writing – original draft. DC: Writing – review and editing, Validation, Visualization. LL: Investigation, Conceptualization, Writing – review and editing. DM: Data curation, Investigation, Writing – review and editing

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by the Integrated Project of the National Natural Science Foundation of China's Enterprise Innovation and Development Joint Fund, “Enrichment Mechanism and Three-dimensional Development Method of Shale Oil in Continental Fault Depression Lake Basins” (U24B6002), and the project of Sinopec Group, “Research on the Development Characteristics and Evolution of the Paleogene Funing Basin in Northern Jiangsu” (P14073).

Conflict of interest

The 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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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

    CuiJ. W.ZhangZ. Y.LiuG. L.ZhangY.QiY. (2022). Breakthrough pressure anisotropy and intra-source migration model of crude oil in shale. Mar. Petroleum Geol.135, 105433. 10.1016/j.marpetgeo.2021.105433

  • 2

    FengD. J. (2022). Geological characteristics and exploration direction of Continental shale gas in Jurassic daanzhai member, sichuan basin. Petroleum Geol. and Exp.44 (2), 219–230. 10.11781/sysydz202202219

  • 3

    GongH. J.JiangZ. X.ZhuF. (2024). Quantitative characterization and control factors of shale oil occurrence state in the shale of member 2 of funing formation in the gaoyou Sag,Subei basin. J. Northeast Petroleum Univ.48 (02), 142–143. 10.3969/j.issn.2095-4107.2024.02.005

  • 4

    GuoX. S.LiM. W.ZhaoM. Y. (2023). Shale oil development and utilization and its role in energy industry. Bull. Chin. Acad. Sci.38 (1). 10.16418/j.issn.1000-3045.20221027001

  • 5

    HuT.JiangF. J.PangX. Q.LiuY.WuG.ZhouK.et al (2024). Identification and evaluation of shale oil micro-migration and its petroleum geological significance. Petroleum Explor. Dev.51 (01), 127–140. 10.1016/s1876-3804(24)60010-8

  • 6

    JiaC. Z.ZouC. N.LiJ. Z. (2012). Assesment criteria,main types,basic features and resource prospects of the shale oil in China. Acta Pet. Sin.33 (3), 333–350.

  • 7

    JiangQ. G.LiM. W.QianM. H. (2016). Quantitative characterization of shale oil in different occurrence states and its application. Petroleum Geol. and Exp.38 (6), 842–599.

  • 8

    JiangF. J.HuoL. N.ChenD.CaoL.ZhaoR.LiY.et al (2023a). The controlling factors and prediction model of pore structure in global shale sediments based on random forest machine learning. Earth-Science Rev.241, 104442. 10.1016/j.earscirev.2023.104442

  • 9

    JiangF. J.JiaC. Z.PangX. Q.JiangL.ZhangC.MaX.et al (2023b). Upper Paleozoic total petroleum system and geological model of natural gas enrichment in ordos basin, NW China. Petroleum Explor. Dev. Online50 (2), 281–292. 10.1016/s1876-3804(23)60387-8

  • 10

    JiangF. J.HuM. L.HuT.LyuJ.HuangL.LiuC.et al (2023c). Controlling factors and models of shale oil enrichment in Lower Permian fengcheng formation, mahu sag, junggar basin, NW China. Petroleum Explor. Dev.50 (04), 812–825. 10.1016/s1876-3804(23)60430-6

  • 11

    LangD. S.GuoS. S.MaD. H. (1996). The correlation method of pyrolysis parameters to evaluate hydrocarbon-beating samples from reservoirs. Mar. Orig. Pet. Geol.1 (1), 53–55.

  • 12

    LiuS. L.DuanH. L.ZhangY. (2014). Analysis of oil and gas exploration potential in F2 M ember Continental shale of subei basin. Offshore oil.34 (3), 27–33. 10.3969/j.issn.1008-2336.2014.03.027

  • 13

    LiuX.AnF.ChenQ. H. (2016). Technical analysis of enhancing crude oil recovery in tight oil reservoirs: a case study of tight oil in the bakken formation. Petroleum Geol. Dev. Daqing35 (06), 164–169. 10.3969/j.issn.1000-3754.2016.06.031

  • 14

    QiuX. M.LiuY. R.FuQ. (2005). Sequence stratigraphy and sedimentary evolution of the Upper Cretaceous and third systems in subei basin. Beijing: Geological Publishing House, 1–59.

  • 15

    SonnenbergS. A.PramuditoA. (2009). Petroleum geology of the giant elm coulee field, williston basin. AAPG Bull.93 (9), 1127–1153. 10.1306/05280909006

  • 16

    U.S. Energy Information Administration (2016). Lower 48 states shale plays. Washington D C: U.S. Energy Information Administration.

  • 17

    WangZ. L. (2013). Research progress, existing problem and development trend of tight rock oil. Petroleum Geol. and Exp.35 (6), 587–596. 10.11781/sysydz201306587

  • 18

    WangH. F. (2016). Recognition of effective fractures within the oil shale in the fourth member of funing formation in northern Jiangsu basin. J. Southwest Petroleum Univ. Sci. and Technol. Ed.38 (3), 21–27. 10.11885/j.issn.1674-5086.2014.08.27.02

  • 19

    WangH. W.DuanH. L. (2016). Formation condition and enrichment rule of shale oil in the second member of funing formation in yancheng sag. Complex Hydrocarb. Reserboirs9 (3), 14–19. 10.16181/j.cnki.fzyqc.2016.03.004

  • 20

    WangS. Q.SunX. M.DuJ. Y. (2012). Analysis of structural styles in northern segment of tancheng—lujiang fault zone. Geol. Rev.58 (3), 414–425. 10.16509/j.georeview.2012.03.007

  • 21

    WangX.ZhangJ. F.JiangW. L. (2016). Transverse structures features of different depths derived from bouguer gravity anomalies in the southern segment of tan-lu fault zone. Seismol. Geol.38 (2), 370–385.

  • 22

    WangR.HeW. J.ZhaoX. M. (2022). Geological section analysis of shale oil in lucaogou formation of well-Ji-174, junggar basin. Reserv. Eval. Ddvelopment12 (1), 192–203. 10.13809/j.cnki.cn32-1825/te.2022.01.017

  • 23

    WangZ. Q.FanZ.ChenX. (2023). Global oil and gas development in 2022: situation, trend and enlightenment. Petroleum Explor. Dev.50 (5), 1016–1031. 10.1016/S1876-3804(23)60456-2

  • 24

    WuS. T.ZhuR. K.LuoZ. (2022). Laminar structure of typical Continental shales and reservoir quality evaluation in central-western basins in China. China Pet. Explor.27 (5), 62–72. 10.3969/j.issn.1672-7703.2022.05.006

  • 25

    WuY. Q.JiangF. J.HuT.XuY.GuoJ.XuT.et al (2024). Shale oil content evaluation and sweet spot prediction based on convolutional neural network. Mar. Petroleum Geol.167, 106997. 10.1016/j.marpetgeo.2024.106997

  • 26

    WuY. Q.HuT.JiangF. J.GuoJ.WangF.QiZ.et al (2025a). Lacustrine records of Paleocene-Eocene thermal maximum (PETM) triggered by volcanic activity. Org. Geochem.200, 104899. 10.1016/j.orggeochem.2024.104899

  • 27

    WuY. Q.JiangF. J.XuY. L.GuoJ.XuT.HuT.et al (2025b). Middle Eocene climatic optimum drove palaeoenvironmental fluctuations and organic matter enrichment in lacustrine facies of the Bohai Bay basin, China. Palaeogeogr. Palaeoclimatol. Palaeoecol.659, 112665. 10.1016/j.palaeo.2024.112665

  • 28

    XiongS. C.ChuS. S.PiS. H. (2017). Micro-pore characteristics and recoverability of shale oil reservoirs. Earth Sci.42 (8), 1379–1385. 10.3799/dqkx.2017.550

  • 29

    YanL.RanQ. Q.GaoY. (2017). Shale oil occurrence form and recoverability evaluation of shale oil reservoir in lucaogou formation of Xinjiang. Reserv. Eval. Dev.7 (6), 20–26. 10.13809/j.cnki.cn32-1825/te.2017.06.004

  • 30

    YangH.LiS. X.LiuX. Y. (2013). Characteristics and resource prospects of shale oil and shale oil in ordos basin. Acta Pet. Sin.34 (1), 1–11. 10.7623/syxb201301001

  • 31

    YunL.HeX. P.HuaC. X. (2023). Accumulation characteristics and resource potential of Paleogene Continental shale oil in qintong sag of subei basin. Acta Pet. Sin.44 (1), 174–187. 10.7623/syxb202301011

  • 32

    ZanL.LuoW. F.MaX. D. (2016). Hydrocarbon generation potential and genetic environments of second member of funing formation in qintong Sag,Subei basin. Unconv. Oil and Gas3 (3), 1–8.

  • 33

    ZhangX. L.ZhuX. M.ZhougD. K. (2004). The character of sequence framework of tertiary and Upper Cretaceous in gaoyou sag,subei basin. Acta Sedimentol. Sin.22 (3), 393–399. 10.14027/j.cnki.cjxb.2004.03.004

  • 34

    ZhangJ. D.YangC. C.LiuC. Z. (2010). The deep structures of strike-slip and extension faults and their composite relationship in the southern of tanlu fault zone. Chin. J. Geophys.53 (4), 864–873. 10.3969/j.issn.0001-5733.2010.04.011

  • 35

    ZhangL. Y.BaoY. S.LiJ. Y.LiZ.ZhuR.ZhangJ. (2014). Movability of lacustrine shale oil:a case study of dongying Sag,Jiyang Depression,Bohai Bay basin. Petroleum Explor. Dev.41 (6), 703–711. 10.1016/s1876-3804(14)60084-7

  • 36

    ZhengK. F.PengX. L. (2013). Hydrocarbon accumulation and favorable zone of shale oil and gas in Upper Cretaceous Neogene of north Jiangsu basin. Jiangsu Geol.37 (1), 147–154. 10.14027/j.cnki.cjxb.2004.03.004

  • 37

    ZhouS. W.XueH. Q.GuoW. (2014). Experimenta1 of movable fluid of shale oil reservoir in Jurassic shale oil reservoirs in central sichuan basin. J. Liaoning Tech. Univ. Nat. Sci.33 (6), 768–772. 10.3969/j.issn.1008-0562.2014.06.010

  • 38

    ZouC. N.ZhuR. K.BaiB. (2015). Significance,geologic characteristics.resource potential and future challenges of shale oil and shale oil. Bull. Mineralogy,Petrology Geochem.34 (1), 3–17. 10.3969/j.issn.1007-2802.2015.01.001

Summary

Keywords

funing formation, Subei Basin, shale oil, occurrence state, recoverability

Citation

Jin-Ning P, Chongjiao D, Longlong L and Mingzhe D (2025) Occurrence and recoverability of shale oil in Paleogene formation, Subei Basin. Front. Earth Sci. 13:1650751. doi: 10.3389/feart.2025.1650751

Received

20 June 2025

Accepted

18 August 2025

Published

19 September 2025

Volume

13 - 2025

Edited by

Tao Hu, China University of Petroleum, Beijing, China

Reviewed by

Shuai Yin, Xi’an Shiyou University, China

Yuqi Wu, China University of Petroleum, Beijing, China

Kong Deng, China University of Geosciences, China

Updates

Copyright

*Correspondence: Peng Jin-ning,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics