1 Introduction

Tight sandstone gas (TSG) as an unconventional gas source has an extensive distribution and large resource potential worldwide (Jia et al. 2012; Dai et al. 2012a, b; Zou et al. 2012a, b; Arthur and Cole 2014). The term “tight sandstone gas” reservoirs usually refers to reservoirs with low permeability (less than 0.6 md), low porosity (less than 10%) and a complicated pore structure (Guo et al. 2012; Lyu et al. 2017). It is estimated in 2012 that the current proven reserve of tight sandstone gas in the world is 3921×1012 m3 (Zou et al. 2012a, b). To date sixteen, large tight sandstone gas fields (TSGFs) have been discovered in China in areas such as Ordos Basin, Sichuan Basin and Tarim Basin and the potential for the discovery of very large resources has led to increasing exploration emphasis on these plays (Dai et al. 2012a, b; Qiu et al. 2013). TSG reservoirs are frequently distributed continuously or semi-continuously on the slopes and in the centers of the basins and commonly have no distinct trap boundaries (Li et al. 2013).

TSG exploration in Sichuan Basin has been dominated by Xujiahe Formation of Upper Triassic in West Sichuan Depression since the 1940s. At present, hydrocarbon resources are basically proven, and explorative extent remains stagnant (Xie et al. 2014; Wang et al. 2016). In 2019, with the discovery of new strata (Jurassic Shaximiao Formation) of high-yield industrial airflow in QL-16 and QL-17 wells, several Jurassic fluvial facies sand body in Jinhua, Qiulin and Yanting areas broke through successively (4.2–8.8 million m3/day in single well; data come from Petro China Southwest Oil and Gasfield Company) and becomes the newest and largest producing TSG field in Sichuan basin, indicating the exploration and exploitation potential of Jurassic TSG in central Sichuan basin. Structure, fracture, and source rock are associated with the TSG distribution (Li et al. 2018; Zhao et al. 2019). However, the mechanism of tight gas enrichment in the central Sichuan Basin is not clear, multistage fluid activity is complex. The geochemical index of accurate restoration of paleo-fluid evolution can be used to study the mechanism of oil and gas accumulation (Haszeldine et al. 1984; Kareem et al. 2019). In particular, the indication of paleo-pressure evolution for the long-term effective preservation conditions as the key evidence (Guo et al. 2016) provides theoretical basis for further exploration and steady production of TSG in Sichuan Basin. The main control factors of natural gas accumulation can be systematically discussed by fluid quantification (Huang et al. 2020), in this regard, the pressure, temperature and compositional data on the geological fluids (diagenetic fluids and hydrocarbon fluids) can be determined by analysis of host minerals and fluid inclusions (Roedder and Bodnar 1980).

Laser Raman micro-spectroscopy (LRM) is a rapid, high precision, and nondestructive micro-analytical technique capable of in situ analysis (Zozulya et al. 1998). Because Raman spectroscopy is capable of quantitative analyses over a wide range of pressure-temperature conditions, it can directly determine the compositions and total densities of inclusions at both low and high temperatures (Seitz et al. 1996; Wang et al. 2019). Microthermometric measurement and laser Raman spectroscopy analyses performed on fluid inclusions can provide quantitative data on fluid inclusions (e.g., compositions, entrapped pressure, and bulk density), which when combined with thermodynamic modeling, can reconstruct P-T entrapment circumstance (Dubessy et al. 1989, 2001; Azbej et al. 2007). Laser confocal Raman micro-spectroscopy is used to measure CH4 and CO2 solubility, salinity, pressure, and density in individual gaseous and aqueous inclusions by means of peak shift, peak separation, peak area ratio, peak height ratio, and other parameters (Dubessy et al. 2001; Burke 2001; Baumgartner and Bakker 2009; Huang et al. 2020). Significantly, the paleo-pressures of methane inclusions can be accurately calculated with the approach of laser Raman spectroscopy, which provides direct evidence of the phase transformations and fluid evolution even in ultra-deep reservoirs.

The objectives of this study are to quantitatively characterize the changes in the reservoir fluid compositions, pressures, properties, and phases in the reservoirs and to define the fluid accumulation processes in tight sandstone reservoirs in the central Sichuan Basin. The principle research methods were mainly microthermometry, petrographic and optical observations (e.g., transmission light; cathodoluminescence), and laser Raman micro-spectrometry (LRM). These techniques were used to characterize the fluid inclusions according to their compositions, distributions, and host minerals. Micro-thermometry and quantitative Raman spectroscopy were used to determine the compositions, pressures, and thermal properties of the CH4 and CO2 gas fluid inclusions and these data were integrated with thermodynamic modeling to reconstruct the pore fluid pressure and overpressure history of the Shaximiao Formation. These results can then be applied to determine the timing of gas charging events, the reservoir pressure history, and analysis of the gas retention and gas loss in tight gas reservoirs in the Jurassic Shaximiao Formation.

2 Geological setting

The Sichuan Basin is a large oil and gas bearing intracratonic basin, which located in the Upper Yangtze stable craton, western China (Fig. 1a). The Jurassic Shaximiao Formation reserves reach to 2744.95×108 m3 and the daily gas production of individual intervals vary from 2 to 83 million cubic meters (data reported from Petro China Southwest Oil and Gasfield Company). The Qiulin and Jinhua tight gas fields located in the gentle structural belt of the Central Sichuan Basin (Fig. 1b). The Shaximiao Formation can be divided into the Upper Shaximiao Formation (J2s2) and Lower Shaximiao Formation (J2s1; Fig. 2). The Shaximiao Formation is composed primarily of purplish red and light-gray mudstone, gray-green fine sandstone and siltstone. Three groups of channel sand bodies (named No. 6, 7, and 8: Petro China Southwest Oil and Gasfield Company) are the main producing intervals of tight sandstone gas in the Shaximiao Formation. The different sand bodies are strongly heterogeneous longitudinally (Fig. 1c) and are in direct contact with Jurassic and Triassic source rocks through a series of faults (Figs. 2 and 3).

Fig. 1
figure 1

a: Location of Sichuan Basin in China. b: Location of study area in central Sichuan Basin. c: Distribution of sand bodies of Jurassic gas fields, sampled wells, and the location of the geological cross section Fig. 2) in the central Sichuan Basin

Fig. 2
figure 2

Geological cross section of the wells in Jurassic Qiulin and Jinhua gas field

Fig. 3
figure 3

Jurassic stratigraphic column of central Sichuan Basin (adopted from Li et al. 2011; Xiao et al. 2019)

The tectonics and structure of the Upper Triassic-Jurassic sediments in the Sichuan Basin were influenced by three tectonic episodes: the Indosinian Event (late Triassic), the Yanshan Event (Jurassic), and the Himalayan (early Cenozoic). The formation of tectonic depressions occurred mainly in the Indosinian period. During the Yanshanian period, the strata underwent uplift and erosion, and the degree of denudation increased (Jiang et al. 2018). The Himalayan movement was the most intense tectonic event in basin and led to the occurrence of local tectonic deformation in the Upper Triassic and Jurassic strata. This local tectonism had significant influence on the migration, occurrences, and fluid properties of the natural gas in the basin (Guo et al. 2012; Wang et al. 2015).

Exploration results demonstrate that the reservoir conditions in the channel sands are favorable for shallow oil and gas accumulations (Guo et al. 2018). The porosity of the reservoirs in the Shaximiao reservoirs ranges from 0.4 to 9.6 %, with an average of 4.5 %, and the average permeability is approximately 0.23 mD, both of which are characteristic of classical ultra-low porosity and ultra-low permeability tight reservoirs. The primary types of storage space include residual primary intergranular pores, intergranular pores, and micro-fractures (Liu et al. 2010). Gas/source correlations show that the reservoir gases were primarily generated in the lower Triassic Xujiahe source rocks (Xiao et al. 2019). The reservoired gas is rich in methane (84.89–97.34%), and C2+ concentrations range from 2.42-13.53%. The non-hydrocarbon gases are primarily N2 and CO2, with a total content ranging from 0.00-1.49% (Xiao et al. 2019).

3 Samples and analytical methods

Approximately 15 TSG core samples were collected from the QL-17, QL-18, and JH-9 borehole of Shaximiao Formation (J2s) in Central Sichuan Basin (Fig. 1). Standard logging interpretation was used to select specific gas producing intervals for sampling. The reservoirs samples are mainly multiperiod channel sand bodies. Sample depths from the QL-17 well range from 2166.8 to 2168.2 m, sample depths from the QL-18 well vary from 2087.1 to 2093.1 m and sample depths from the JH-9 range from 2911.82 to 2912.38 m (Table 1).

Table 1 Pure CH4 inclusions Raman quantitative analysis summary

Doubly polished thin sections (~200 µm) were prepared for fluid inclusion studies using Canada balsam. Residual organic matter and migrated hydrocarbons in the samples can produce fluorescence effects which can obscure the weaker Raman signals from the compositional material (Frezzotti et al. 2012). To mitigate organics interference, core sheets were immersed in methylene chloride and methanol for 3–5 hours until the sheets are detached from the glass.

A petrographic examination was carried out under transmitted visible light with up to 10X, 20X, 50X, and 100X magnification using an Olympus BX51 microscope. Thermodynamics measurements of the fluid inclusions were carried out using TH-600 heating–cooling stage (Linkam). The homogenization temperatures (Th) of fluid inclusions were measured by thermodynamic cycle and standard deviation of instrument keep in 1°C (Shepherd et al. 1985). Temperature controlled in 0.5 °C/min for unit until the phase differentiation disappeared or freezing point appeared. The heat–cooling stage stood stable for 2 min and obtained the Th, repetitive increase or decrease temperature for verification. The phase transitions during thermal measurement and molar volumes of binary CO2-CH4 system are described by Thiery et al. (1994) and van den Kerkhof (2001).

Cathode luminescence (CL) analysis can be used to distinguish quartz from calcite and to identify different generations of cement (Mason 2006). A selection of polished thin sections was analyzed using a CL8200 MK5 cathode luminescence system coupled to a Leica DM2500 camera. The beam conditions were: beam voltage: 17 kV, beam current: 500 μA. All CL images were obtained with the same exposure time (70 ms) and gain (1.005 s).

Laser Raman micro-spectroscopy (LRM) was used to analyze the compositions of the gases trapped within fluid inclusions in fracture-healed quartz or cements. The Raman system was a JY/Horiba LabRam HR800 (frequency multiplication neodymium-doped yttrium aluminum garnet laser, 532.06 nm, output laser power of 400–500 mW) at Wuhan University of Geosciences, China and objective with a 0.5 numerical aperture. Gaseous fluid inclusions were investigated using a 200 µm confocal aperture, and Raman peak calibration was accomplished with the 520.7 cm−1 band of a polished silicon wafer. The compositions of the gases in the fluid inclusions were analyzed using the 300 gr·mm-1 grating and 2650 cm-1 grating center. Spectrum was collected in the 1200–4000 cm−1 range. The collection time ranged from 0 to 50 seconds and 200–500 scans were acquired per sample in order to acquire the signal-to-noise ratio. For single-phase inclusions, high-resolution techniques (an 1800 gr·mm-1 grating and narrow slit of 50 μm) were used to enhance the resolution of the peak positions. A Ne lamp was laid under the objective to collect the Ne emission line spectra for wavelength calibration purposes (Kawakami et al. 2003; Dubessy et al. 2012; Zhang et al. 2016).

The regular shift of methane stretching vibration peak (v1) with density fluctuation can be used to determine the density of pure methane gas inclusions (Seitz et al. 1993, 1996; Lin 2005; Lu et al. 2007) at 25 °C using homogenized samples. After the acquisition of the Raman spectra, PVT-x of pure methane gas inclusions can be determined by following the methodology described by Zhang et al. (2016). The densities of single-phase CO2 inclusions were determined from the differences in wavenumber between the two peaks of the Fermi diad of CO2, with the application of the neon peak shift correction (Rosso and Bodnar 1995). The spectra of the CO2 inclusions were collected at 160 °C to ensure complete homogenization and the densities of CO2 calculated using the multiple regression methods described by Huang et al. (2018).

The burial and thermal histories were modeled using the BasinMod-1D (Version 7.06) software, which integrates the stratigraphic thickness, lithology, absolute age, erosional thickness, and measured borehole temperature data in order to calculate the heat flow and vitrinite reflectance histories through time. Measured vitrinite reflectance values and borehole temperature data were obtained from well completion reports supplied by the Petro China Southwest Oil and Gasfield Company. The modeling process was primarily calibrated using measured vitrinite reflectance values and corrected borehole temperatures.

4 Results

4.1 Fluid inclusion characteristics

4.1.1 Petrology of fluid inclusions

Both CH4 and CO2 bearing fluid inclusions can be trapped together in micro-fractures and overgrowths on quartz crystals (Fig. 4). Primary fluid inclusions are those that are trapped during the initial diagenesis and crystallization of the rock, while secondary origin inclusions trapped during overgrowth and precipitation of fracture- and vein-filling cements. Cathode luminescence photomicrographs show that the color of quartz overgrowths is different from that of quartz grains, which assists in identifying fluid inclusions, and the entrapment zone of secondary fluids can be further defined by combining cathode luminescence photomicrographs with those from single polarized light (Fig. 5).

Fig. 4
figure 4

a Pure CH4 gas inclusions exist as linear distribution in healed micro-fracture coeval with aqueous inclusion; b CH4 gaseous inclusions exist in quartz overgrowth and CO2 gaseous inclusions exist isolated in grain; c Pure CO2 gas inclusions exist in overgrowths of quartz; d Abundant pure CO2 inclusions distribute as clusters in quartz grain; e Hybrid CH4-CO2 gaseous inclusions exist in quartz overgrowths; f and g Hybrid CH4-CO2 gaseous inclusions exist in healed micro-fractures; h Macro perspective shows that secondary CH4 gas inclusions cross over the quartz grain and overgrowths. CO2 gas inclusions only exist in overgrowths; i: Two types of hybrid inclusions with different molar ratios of CH4 and CO2 are existed in two intersecting fractures. All components were identified by micro-Raman spectroscopy

Fig. 5
figure 5

Cathodoluminescence and plane-polarized light photomicrographs of reservoir sandstones in Shaximiao Formation; Qz: Quartz, Cc: Carbonate cements, Fel: Feldspar

Fluid inclusions assemblages (FIA) can be classified into four types based upon petrographic observations at room temperature and Raman analysis.

Inclusion type I: This FIA type is single-phase (VCH4) gas inclusions found along healed microfractures and coeval with aqueous inclusions (Fig. 4a, b) implying secondary origin. Diameters of gas inclusions range from 5 to 11 μm.

Inclusion type II: These vapor-rich inclusions are two-phase inclusions (LCO2+VCO2) and coeval with aqueous inclusions. They primarily occur as clusters in quartz grains or linear arrays in quartz overgrowths (Fig. 4b, c, d). The diameters vary from 8–15 μm. These inclusions represent both primary and secondary origins, but with different timing of formation. They generally have elliptical to circular shapes. In Fig. 4h, secondary pure CH4 gas inclusions cut across the quartz grain and overgrowths and isolated inclusions of pure CO2 occur in overgrowths, which indicates CH4 fluid was trapped later than CO2 fluid.

Inclusion type III: These inclusions contain CH4 vapor (V) or CO2 vapor. These inclusions are commonly coeval with aqueous inclusions (Fig. 4e, f, g) with sizes ranging from 7–13 μm and generally occur along microfractures and in quartz overgrowths. In Fig. 4i, a mixed gas inclusion exists in the position of the intersection of two conjugated healed fractures (horizontal micro-fracture and the vertical micro-fracture). On account of long axis direction of fluid inclusions in accordance to the fracture distribution of FIA2, which indicates that FIA2 was grew later than FIA1 from petrographic examination (Goldstein 2003).

Inclusion type IV: Rich nitrogen and bearing small amounts of CO2 inclusions (Fig. 8d) distributes along microfractures. The diameters range in size from 5 to 9 μm.

4.1.2 Microthermometry

Microthermometric measurements were performed on aqueous inclusions coexisting with the four types of gas inclusions described above. The salinities (wt.% NaCl equivalent) of these fluid inclusions were calculated based on their final ice-melting temperatures (Tm, ice; Hall et al. 1988). The Raman and microthermometric data and the calculated density and pressure parameters for these fluid inclusions are summarized in Tables 1, 2 and illustrated in Fig. 6. The homogenization temperatures of aqueous inclusions coeval with type I inclusions range from 115.0 to 160.0 °C, with a well-defined maximum distribution between 120 and 130 °C. The homogenization temperatures of aqueous inclusions coeval with type II inclusions range from 175.0 to 230.0 °C, with a well-defined maximum distribution between 185 and 190 °C and 210.0–220.0 °C. Th of aqueous inclusions coeval with type IV inclusions range from 200.0 to 210.0 °C (Fig. 6). Tm of aqueous inclusions coeval with pure CH4 inclusions vary between − 10.1 and − 13.8 °C, with salinities ranging from 14.0 to 17.6 wt.% NaCl equiv. The ice-melting temperature of aqueous inclusions coeval with pure CO2 inclusions varies between − 4.9 and − 1.7 °C, with calculated salinities ranging from 0.92 to 2.30 wt.% NaCl equiv.

Table 2 Pure CO2 inclusions Raman quantitative analysis summary
Fig. 6
figure 6

Homogenization temperature statistics of aqueous inclusions in Shaximiao Formation

Inclusions with mixed CH4-CO2 gaseous components were also found to occur. These samples were initially distinguished by their C–H symmetric stretching band (~2917 cm-1) and Fermi diads (~1280 cm-1 and ~1380 cm-1) at 100.0 °C (keep inclusions at homogeneous state). Subsequently, decreased temperature to − 120.0 °C, complete freezing inclusions into solid (S) phase. During this thermal measurement, the carbonic inclusions homogenized to a liquid (L) phase, on the contrary, CH4-rich inclusion homogenized to a gas phase. In the terminal-origin hybrid gas inclusions (e.g., FIA2), the melting of CO2 (Tm, carb) (L+V→SV→LV) occurred between − 67.9 and − 69.4 °C. Since Tm, carb is lower than that of single-phase CO2 (− 62.3 to − 63 °C), the reduction of the freezing point indicates that CH4 must also be present in the inclusion and take responsibility of the observed Tm reduction. Homogenization temperatures (Th, carb) of hybrid gas inclusions (yellow dots in Fig. 7) vary from − 12.3 to − 19.8 °C and − 63.1 to − 69.2 °C. The concentration of CH4 can be determined from the Tm, carb and Th, carb as described by Thiery et al. (1994). Thermodynamics measured results show molar proportion of CH4 not exceeding 30% and the molar proportion of CO2 vary from 75% to 77%. Contrast with their petrographic characteristics, significantly higher concentrations of CH4 (75%–77%) inclusions were entrapped later than the low concentrations of CH4 (23%–26%) inclusions (Fig. 7).

Fig. 7
figure 7

adopted from Huang et al. 2020

Photomicrographs of phase transition processes during determination of melting temperature (Tm, carb) and homogenization temperatures (Th, carb) for hybrid gas inclusions in the Shaximiao Formation. Inclusions of different petrographic generations shown different gas concentrations. Black lines define the melting temperature of the solid phase, dashed lines represent isothermal conditions of homogenization; V(L) = stability field of the vapor phase; L = V = liquid and vapor phase are in equilibrium; L(V) = stability field of the liquid phase; L(SV) = tie-line joining the field of liquid phase in coexistence with a vapor and a solid phase; for a detailed explanation refer to Thiery et al. 1994, the original chart is

4.1.3 Gaseous inclusion quantitative analysis by Raman spectroscopy

Pressure and density of methane fluid inclusions in quartz minerals can be measured using quantitative Raman analysis (Huang et al. 2020; Zhang et al. 2016). Figure 8a-a shows a typical Raman spectrum of single-phase methane inclusions acquired at homogeneous state. Table 1 lists the data for 13 type I inclusions, including the data of densities, determined from the v1 symmetric stretching band of CH4. Overall, the density of secondary pure CH4 fluid inclusions in quartz minerals range from 0.195 to 0.237 g/cm3 (average value: 0.222 g/cm3).

Fig. 8
figure 8

a: Spectrum of pure CH4 collected using 300 gr.mm−1 grating in quartz; (a’): wave number calibration using the Ne emission spectrum; (a’’): CH4 symmetric stretching band (v1) peak was collected using 1800 gr.mm−1 grating with Ne lamp; b: Spectrum of pure CO2 acquired using a 300 gr.mm−1 grating in a homogeneous state; (b’) Ne emission spectrum for wave number calibration; (b”): CO2 Fermi double peaks were collected using 1800 gr.mm−1 grating with Ne lamp; c: Spectrum of a hybrid gas inclusion collected using a 300 gr. mm−1 grating; d: Spectrum of a N2-rich inclusion collected using 300 gr.mm−1 grating in quartz

Table 2 shows the data for nine pure CO2 inclusions, which were analyzed by the Fermi diads (double peaks) of CO2 at homogenization condition (Fig. 8b). The pure CO2 inclusions can be divided into two classes primary CO2 inclusions with densities ranging from 0.874 to 1.020 g/cm3 (average value: 0.941 g/cm3) and secondary inclusions with densities varying between 0.514 and 0.715 g/cm3 (average value: 0.651 g/cm3). The homogenization pressures were calculated using the Peng-Robinson equation of state (EOS) model (1960), which uses the densities and homogenization temperatures of the coexisting aqueous inclusions.

The occurrence of fluid inclusions with highly variable gas-liquid ratios in a single fluid inclusion assemblage (FIA) indicates that the pure CH4 and CO2 inclusions and the aqueous saline inclusions were trapped under immiscible entrapment conditions. Consequently, homogenization P-T is correspond to the entrapment conditions (Diamond 2001; Goldstein 2003). As a result, the thermodynamic systems integrate with the Raman spectroscopy calibration can determine the compositions and densities of the multiphase and multicomponent fluids (Peng and Robinson 1960; Duan and Mao 2006; Becker et al. 2008; Lecumberri-Sanchez et al. 2012; Mao et al. 2013).

Pressure coefficient (Pc) here refers to the actual formation fluid pressure (Pf) divided by normal hydrostatic pressure (Phyd) in the same position, and it is one of the primary indicators of abnormal reservoir pressure (Law and Spencer 1998). Trapped pressures of pure CH4 inclusions vary between approximately 44.0 and 58.5 MPa (Table 1). Estimated paleo-depth against the trapped pressure of fluid inclusions plotted the paleo Pc (Fig. 9, hydrostatic gradient of 10.00 MPa/km and the lithostatic gradient of 24.99 MPa/ km for comparison). Take into accounts of density log data for the well QL-17 calculated with a depth-averaged rock density and a water column with an average seawater density, the lithostatic gradient was 2.55 g/cm3, and the hydrostatic gradient was 1.01 g/cm3. As can be seen in Fig. 9, most of the samples are weakly overpressured.

Fig. 9
figure 9

Calculated entrapment pressures plotted against inferred paleo-burial depth (Pc: Pressure coefficient)

4.2 Basin model

The burial and thermal histories were simulated to determine the present-day porosity–depth relationship in order to correct the present-day burial depths for decompaction effects (Falvey 1981). The initial porosities and lithologies (shale/mudstone, sandstone, siltstone, and limestone) were adopted as the default values in the BasinMod software. The depth and absolute age data for each stratigraphic unit were obtained from Petro China Southwest Oil and Gasfield Company. The eroded thickness, tectonic subsidence, hydrocarbon generation, paleo-burial depth and paleo-heat flow data were obtained from Guo et al. (2018) and Lin et al. (2017). The thermal conductivities of the rock units can be used to calculate the present-day heat flow based on the transient heat flow model of the BasinMod 1D software. The good correlation between the modeled vitrinite reflectance and temperature and the measured data implies that the final calculated thermal history model is suitable for the study area (Fig. 10).

Fig. 10
figure 10

Burial history and thermal history of the Shaximiao Formation in well QL-17

The reconstructed burial and thermal histories for QL-17 show that the oil generation (Early mature) in the J2s Formation began at ~165 Ma and a temperature of 85 °C. The maturity of the J2s-T3x Formation reached 0.7% Ro at ~145 Ma. The thermal maturity of the J2s-T3x Formation reached 1.3% Ro at 85 Ma (later Cretaceous), which is consistent with the end of the oil generation stage. The uplift and erosion occurring during the Himalayan tectonic event resulted in the cooling of the rocks in the J2s-T3x Formations and the cessation of the thermally driven processes, including petroleum generation. The Himalayan tectonic movements caused erosional thicknesses exceeding 1700 m, beginning in the later Cretaceous (Wang et al. 2016), which resulted in a reduction in overpressure. In the modeled pressure history of J2s1 (Fig. 11), the yellow-filled points are the results for the CH4 inclusions, and the red-filled points are pressure data from wells. In the model, the pressure began to rise at the end of the middle (160 Ma) and increased rapidly in late Jurassic to Cretaceous reaching the maximum of ~ 53 MPa at ~ 70 Ma before the onset of the Himalayan tectonic uplift. With the rapid uplift beginning in the Late Cretaceous, the pressure in the J2s1 declined rapidly, reaching normal pressure conditions in the Neogene. In addition to the reduction in overburden stress, factors such as reduced pore fluid pressure, temperature reduction, trap space enlargement, and the pore volume expansion all also contributed to the pressure decline (Liu et al. 2016).

Fig. 11
figure 11

The modeled pressure evolution of the J2s1 Formation (dashed lines). The yellow-filled dots represent the trapping pressure of CH4 inclusions, and the red-filled dot are measured well data (primarily from drill stem tests)

5 Discussion

5.1 Fluid evolution in tight sandstone reservoirs

The combination of petrographic observations and in situ laser Raman micro-spectroscopy indicates that high concentrations of primary and secondary CO2 and CH4 bearing fluid inclusions can occur in the quartz grains in this study. The primary CO2 inclusions, with high densities (0.874 to 1.020 g/cm3) and high homogenization temperatures (>210 °C), were trapped during early diagenesis and mineral crystallization. In contrast, the secondary CO2 inclusions, with low densities (0.514–0.715 g/cm3) and low homogenization temperatures (~180–200 °C), were trapped during secondary fluid migration, after the primary diagenetic processes were complete. The Th throwing-dot method can be used to estimate the gas and oil accumulation times (Haszeldine et al. 1984b; Horsfield and McLimans 1984). The relatively narrow range of homogenization temperatures for the aqueous inclusions coexisting with methane inclusions indicates that there was a single episode of charging of the hydrocarbon gas into the tight sandstone was accumulated at ~ 75–65 Ma (Fig. 10). The timing of accumulation of the hydrocarbon gas is later than the fluid charging of the CO2. The homogenization temperatures of aqueous inclusions coeval with CO2 gas inclusions are significantly higher than the modeled geothermal gradients, which is an indication of possible hydrothermal activity. The CO2-rich hydrothermal fluid is thought to be the result of the upward migration of inorganic mantle CO2 gas along the faults, since Permian igneous rocks are found in this area.

The petrography occurrence of mixed CH4 and CO2 inclusions with variable relative concentrations (Figs. 4i and 7) suggests that the CH4 was charged by a process of displacement, which resulted in the occurrence of primarily CH4 in the present-day Jurassic TSG in the Shaximiao Formation, with only small concentrations of non-hydrocarbon gases (CO2 and N2).

5.2 Pore fluid pressure evolution

Aqueous inclusions and methane-rich gaseous inclusions were observed to occur together in some of the same veins. These pore fluids are immiscible under the inclusion trapping conditions at saturated methane concentrations. These immiscible homogenization temperatures and homogenization pressures are equal to the temperatures and pressures at which the inclusions were trapped in saturated paleoenvironments (Goldstein 1986; Diamond 2001; Hurai 2010). Pore fluid pressure is a key parameter in the assessment of the reservoir conditions (porosity and permeability), and the formation fluid pressures and pressure coefficients are critical parameters in pore pressure calculations.

The trapping temperatures and pressures of fluid inclusion represent a record of the pressure-temperature conditions of the Jurassic TSG. By combining the depth of the formation of the quartz minerals and the trapping pressures of the fluid inclusions, the overpressure conditions of the Shaximiao Formation at depths near to its maximum burial can be determined. The range of trapping pressures of the methane inclusions is 44.0–58.5 MPa (Table 1) and the calculated pressure coefficients range from ~1.05 to 1.43 (average of 1.29; Table 1). The trapping pressures imply that the reservoir was in a low-to-medium overpressure state at maximum burial depths during the early stages of Himalayan event. Due to the very low porosities and permeabilities of tight sandstones, undercompaction can be eliminated as a cause of the overpressure, and hydrocarbon generation is most probable cause. The persistence of overpressure to the present-day also suggests that the preservation of favorable reservoir conditions is the main reason for the occurrence of the TSG gas reservoirs.

5.3 Accumulation processes for charging of TSG

The characteristics of the mixed gas inclusions that were trapped along two groups of healed micro-fractures with different timing demonstrate that the CH4 gradually displaced the early emplaced CO2 to form the present-day methane-filled reservoirs. The source rocks of the Xujiahe Formation entered the gas generation stage at ~ 90–70 Ma (Fig. 10), and the inclusion data and modeling results indicate there was only one charging event. The gas reservoirs reached their maximum fill prior to the Himalayan tectonic uplift and overpressure ensued (Pc: 1.05–1.43, avg. 1.29).

By integrating the analytical data with the modeling results, a model of hydrocarbon accumulation in Shaximiao Formation of Jurassic can be proposed (Fig. 12). The maximum Th interval (120.0–130.0 °C) of the aqueous inclusions indicates late Cretaceous (~ 75–65 Ma) charging event, which is close to the maximum burial depth. The sand body at the bottom of J2s1 is in direct contact with the mudstone of J1l (Fig. 3), and the oil and gas can migrate into the trap by means of face-loading. However, most of the gases in the J2s1 reservoirs are coal-type gases (Xiao et al. 2019) derived from the Xujiahe Formation and which migrated into the J2s1 Formation reservoirs through faults. The tectonic movements also led to the development of a number of stratigraphic traps in the Jurassic Shaximiao, which were well-sealed horizontally and vertically (Wang et al. 2016). It can be seen from the pressure history (Fig. 11) that the tectonic movements during and after the Himalayan period led to the degradation of the gas reservoir conditions and sealing capacity in the TGS reservoirs. While the development of the fault networks enhanced the migration and charging of shallow gas reservoirs, it also led to increased fault leakage leads to the partial loss of shallow gas.

Fig. 12
figure 12

Accumulation mechanism for Jurassic tight sandstone gas in central Sichuan Basin

6 Conclusions

In this study, multiple paleo-fluids were observed in fluid inclusions in quartz grains, which provided data to determine the pressure–temperature–time–composition (P-T-t-x) of the inclusions and subsequently the trapping conditions. This made it possible to determine the pressure, composition, and fill histories of the reservoirs. Based on our results, the following conclusions can be reached.

  1. 1.

    In J2s tight gas reservoir of central Sichuan Basin, hydrothermal CO2 activities can be divided into two stages before Himalayan movement: high density fluid (0.874–1.020 g/cm3) and the low density fluid (0.514–0.715 g/cm3). Pure CH4 FIAs are secondary in origin with densities of 0.195–0.237 g/cm3. The integration of the Th throwing-dot method with the burial history indicates a late Cretaceous accumulation event (methane fluid, ~75–65 Ma) close to the maximum burial depth. Methane displaced the CO2 and became the primary reservoired gas. Entrapment of different mixtures of CH4 and CO2 in inclusions in related petrographic sequences indicates the dynamic nature of the process of methane emplacement and the gradual displacement of CO2.

  2. 2.

    Quantitative Raman analysis can be used to recover paleo-pressures. Pore fluid pressure ranged from 44.0 to 58.5 MPa for CH4-bearing fluid inclusions and the paleo-pressure coefficients varied from 1.05 to 1.43. The low-medium overpressure is noticed to start at the end of the middle Jurassic (~ 160 Ma) and increased rapidly in late Jurassic to Cretaceous caused by thermocatalytic gas generation. With the tectonic uplift in late Cretaceous, the pore pressure released rapidly to normal pressure by the formation temperature reduction and the loss of tight gas. The continued occurrence of overpressure to the present-day also suggests that the favorable preservation condition take responsibility for the enrichment of Jurassic tight gas in central Sichuan Basin.

  3. 3.

    The coal-type gas derived from the T3x source rocks enters the J2s1 by migration through faults and fractures and CH4 gradually displaced the early emplaced CO2 in reservoirs. The initial stage of gas reservoir is characterized by good storage conditions and high oil and gas filling degree. Subsequently, Himalayan movement affected the continual decline in reservoir pressure to the current mild-to-moderately overpressured state and the reservoir sealing properties. Fracture development not only provides fluid transport channels, but also leads to loss of pore pressure.