Natural Gas Accumulation Characteristics in the Linxing Area, Ordos Basin, NW China: Revealed from the Integrated Study of Fluid Inclusions and Basin Modeling

Key Laboratory of Coalbed Methane Resources and Reservoir Formation Process, Ministry of Education, China University of Mining and Technology, Xuzhou 221116, China School of Resources and Geosciences, China University of Mining and Technology, Xuzhou 221116, China Shandong Provincial Key Laboratory of Depositional Mineralization & Sedimentary Minerals, College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, China Shandong Provincial Lunan Geo-engineering Exploration Instituts, Jining 272100, China


Introduction
The study of hydrocarbon accumulation characteristics is an important part of oil and gas geology research, which can indicate the formation and enrichment of oil and gas, and is of great significance for the exploration and development of oil and gas resources [1][2][3][4]. Unconventional oil and gas resources (such as tight gas) account for a growing proportion of global energy consumption, and its accumulation model is more complicated than conventional oil and gas [5][6][7]. In 2017, China's tight gas production was 343 × 10 8 m 3 , mainly from the Sichuan Basin and the Ordos Basin [8]. The tight gas resources in the Ordos Basin are 14:4 × 10 12 m 3 , accounting for 59% of the total tight gas resources in China [9]. Several large-scale tight gas fields have been found in the Ordos Basin, such as the Sulige gas field, the Mizhi gas field, and the Shenmu gas field, which are important tight gas industrial bases in China [9,10]. The Linxing area is a favorable area for the exploration and development of unconventional natural gas from coal measures in the Ordos Basin. The development of coalbed methane in the Linxing area has formed a certain industrial scale, and tight gas has become the focus of next exploration and development [11,12]. In particular, the maturities of the source rocks in the Linxing area are higher than that in the adjacent areas (the vitrinite reflectance of the source rocks is generally greater than 1.0%), and the fractures in the coal and tight sandstones are developed [12][13][14], so the prospect for exploration and development of tight gas is generally promising.
Previous studies have been carried out on the process of natural gas accumulation in the Linxing area, such as hydrocarbon generation-expulsion simulation to determine the accumulation period [13,15], but it is only a general accumulation period, and cannot explain whether natural gas is 1 multistage accumulation. Fluid inclusions are paleo-fluids trapped in the mineral lattice defects, which represent the original physical and chemical characteristics of the paleofluids in its formation period, and act as effective means to determine the hydrocarbon accumulation characteristics [1,[16][17][18][19][20][21]. The predecessors used fluid inclusion analysis to conclude that the gas accumulation in the Upper Paleozoic in the Ordos Basin was mostly continuous accumulation of 2-6 periods [22][23][24][25]. However, the previous research scopes are often large (the whole gas field), resulting in low research accuracy, and it is possible to neglect the different accumulation characteristics of different locations in the same gas field. Under the dual influence of the Yanshanian tectonic-thermal event and the thermaltectonic effect caused by Zijinshan magmatic pluton, the process of gas accumulation in the Linxing area is more complex. Previous studies are limited to the description of current phenomena, such as the vitrinite reflectance (Ro) of source rocks around Zijinshan magmatic pluton decreases annularly [12][13][14][15]26], and have not clarified the specific mechanism controlling gas accumulation.
In this paper, fluid inclusion analysis and basin modeling are used as the main research methods. Compared with previous studies, this paper is not limited to the overall general understanding of the gas accumulation characteristics in the study area but focuses on the difference of gas accumulation characteristics in different locations of the study area under the influence of Zijinshan magmatic pluton. Specifically, first of all, according to the microthermometry characteristics of fluid inclusions, combined with burial-thermal evolution and hydrocarbon generation-expulsion simulation of multiple wells, reliable gas accumulation periods and their regional differences are obtained. Secondly, through the trapping pressure simulation of fluid inclusion, the evolution of reservoir pressure and paleo-fluid potential is restored, and then, the gas accumulation process and characteristics under the influence of Zijinshan magmatic pluton are summarized in detail. Finally, the gas accumulation model under the influence of Zijinshan magmatic pluton is established to provide reference for the natural gas development in the study area.

Geological Setting
The Linxing area is located in the north of the Jinxi fold belt in the Ordos Basin ( Figure 1). The sedimentary environment of the Upper Paleozoic from the bottom to the top gradually evolved from the marine facies to the marine-terrigenous facies and continental facies. According to drilling and outcrop observation, the Upper Paleozoic from the bottom to the top is as follows: Benxi Formation (C 2 b) in Upper Carboniferous, Taiyuan Formation (P 1 t) and Shanxi Formation (P 1 s) in Lower Permian, Shihezi Formation (P 2 h) in Middle Permian, and Shiqianfeng Formation (P 3 s) in Upper Permian (Figures 1 and 2). The Upper Paleozoic is unconformity with the underlying Ordovician Majiagou Formation and the overlying Triassic strata.
Located in the south-central part of the study area, the Zijinshan magmatic pluton mainly intruded into the Ordovician, with an outcrop area of approximately 23.3 km 2 ( Figure 1) [14]. Through SHRIMP dating and isotope analysis, scholars believe that the formation period of Zijinshan magmatic pluton is mainly in the Jurassic-Cretaceous (160 Ma-120 Ma), and some samples are aged from 195 Ma to 160 Ma [27][28][29][30][31][32][33][34]. The formation of Zijinshan magmatic pluton is closely related to strong tectonic thermal events and frequent volcanic activities in the eastern part of North China Plate during the Yanshanian period [30,32,34].
The source rocks are mainly the coal measure source rocks of marine-terrigenous facies in C 2 b, P 1 t, and P 1 s, including coal and dark mudstones. The cumulative thickness of the coal seam is 15 m-24 m, with high maturity and stable distribution [12]. Multiple sets of sandstone reservoirs and mudstone caprocks were developed in the Upper Paleozoic, among which the fluvial sandstone reservoirs in P 1 s and P 2 x are stable thickness and widely gas bearing [13]. The overlying thick mudstone of P 2 s and P 3 s serve as regional caprocks, providing good preservation conditions for gas reservoirs.

Samples and Methods
3.1. Samples and Instruments. In this study, fluid inclusions samples were collected from 65 core samples of tight reservoirs from 11 wells. According to the methods for rock thin section preparation [35] and measurement of homogenization temperature and salinity of inclusions in sedimentary rock [36], the samples were prepared as 1-mm-thick double-polished thin section for microthermometric measurement.
The Axio Scope A1 microscope and Axioskop 40 Pol microscope are used to observe and analyze fluorescence and petrographic characteristics of fluid inclusions. The Linkam THMSG 600 heating/cooling stage (temperature range: -196°C to 600°C) was used to bring the fluid inclusions to a uniform state, and the measurement error of the homogenization temperatures of fluid inclusions was less than 0.1°C. The gas-liquid ratios of hydrocarbon inclusions were measured by Confocal Laser Scanning Microscopy. The gasliquid ratios of aqueous inclusions and hydrocarbon-bearing aqueous inclusions were calculated by the method proposed by Zhou et al. [37], using petrographic microscopy and fluorescent microscopy, obtain a series of two-dimensional images of single fluid inclusion. Details of the approaches, methods, and software used are given by Zhou et al. [37]. 2 Geofluids

Methods
3.2.1. Petrography. The fluid filling time during diagenesis could be established through the analysis of the textural relationships between host minerals and fluid inclusions, and fluorescence and phase characters of fluid inclusions. The hydrocarbon inclusions, which were irradiated by ultraviolet light, produced fluorescence colors in the visible light range (wavelength of 400 nm-700 nm), while the hydrocarbonbearing aqueous inclusions and aqueous inclusions exhibited no fluorescence [38]. If the fluid inclusions were trapped from a uniform state, after the temperature decreased, associated with the fluid composition and density, the single liquid phase, single vapor phase, and vapor-liquid phase may be formed.

Microthermometry.
Microthermometry includes the measurement of homogenization temperature and freezing temperature of fluid inclusions. The homogenization temperature of nonsingle-phase fluid inclusions is the temperature when the inclusion phases transform into a single uniform phase. The freezing temperature is the temperature when the frozen inclusion is heated and the last ice disappears, which is related to the salinity of fluid. The salinity of aqueous inclusions was calculated using a freezing temperature [39]. It should be noted that fluid inclusions of the same type and phase usually have similar homogenization temperature, salinity, gas-liquid ratio, and other characters. The homogenization temperature of hydrocarbon inclusions is slightly lower than its coeval aqueous inclusions, and aqueous  Figure 1: Comprehensive geological map of the Linxing area (the seismic section is cited from [14]).

Simulation of Burial-Thermal Evolution and
Hydrocarbon Generation-Expulsion. The PetroMod 2012 basin modeling software from Schlumberger ltd (http://www .software.slb.com/products/petromod) was used to simulate the burial-thermal evolution and hydrocarbon generationexpulsion.
The required data for reconstruction of burial-thermal includes stratigraphic age, stratigraphic lithology, denudation thickness, paleo-heat flow, and paleo-water depth. The stratigraphic age and stratigraphic lithology reference to Stratigraphic Chronology of Ordos Basin (Changqing Oilfield). Through the logging data analysis, the denudation thickness was restored by using the mudstone interval transit time and depth model proposed by Liu et al. [40]. The initial paleo-heat flow refers to the study of the paleo-heat flow in the northeastern Ordos Basin [41][42][43] (Figure S1). Based on previous studies of sedimentary facies and paleo-water depth of the Ordos Basin [44][45][46], the paleo-water depth evolution is set as shown in Table S1.
Based on the reconstructed burial-thermal evolution, combined with the total organic carbon (TOC) and kerogen pyrolysis analysis data of coal samples (Table S2), the Tissot in Waples [47] TIII Crack kinetic model in PetroMod 1D software was used to simulate the hydrocarbon generationexpulsion of source rock.

Determination of Gas Accumulation Periods.
According to the burial-thermal evolution, the age information of aqueous inclusions can be ascertained based on its homogenization temperatures. Then, the formation time of the coeval hydrocarbon inclusions can be ascertained. Finally, the charging time of the natural gas can be determined [48,49].
3.2.5. Trapping Pressure Simulation. Using PVTSIM 20 fluid simulation software (http://www.calsep.com/services/ multiphase-meter-customized-pvt-software.html), input fluid inclusions' homogenization temperature, gas-liquid ratio, the assumed components of hydrocarbon-bearing aqueous inclusions based on the current gas reservoir components, and the assumed components of hydrocarbon inclusions, then the pressure-temperature (P-T) phase diagrams and isochores of the hydrocarbon inclusions can be obtained. Details of the approaches, methods, and software used are given by Aplin et al. [16,17], Mi et al. [50], and Zhang et al. [51].
3.2.6. Paleo-Fluid Potential Calculation. As functions of fluid pressure, density, elevation, temperature, etc., the fluid potential and the fluid potential gradient are effective means for portraying the trend of fluid migration in petroleum geology. On the basis of pressure recovery, combined with the burial-thermal evolution, the fluid potential distribution during the main accumulation period was recovered using the Hubbert fluid potential equation [52]. Equation (1) is the Hubbert fluid potential equation: where Φ is the fluid potential, m 2 /s 2 ; g is the gravitational acceleration, m/s 2 ; z is the elevation of the calculated position relative to the datum plane (usually the ground), m; p is the fluid pressure, which is approximately equal to the trapping pressure of the fluid inclusions, Pa; ρ is the fluid density, which is approximately equal to the density of the fluids inclusions [53], kg/m 3 ; and v is the fluid velocity, m/s. Considering that the deep underground fluid migration velocity is very low, the value of v is approximately 0. In addition, water, oil, and other liquid molecules are very small and almost incompressible, which means the density can be approximated as a fixed value. Therefore, equation (1) can be simplified to: The simulation of fluid pressure and calculation of density are based on hydrocarbon-bearing aqueous inclusions and aqueous inclusions, and the fluid potential mainly indicates the migration of formation water. However, the migration of oil and gas fluids is highly consistent with the migration of formation water [54,55], and the results have implications for identifying hydrocarbon migration and accumulation. In the study area, 5 types of inclusions, such as aqueous inclusions, hydrocarbon-bearing aqueous inclusions, hydrocarbon inclusions, crystal-bearing aqueous inclusions, and aqueous-carbonic inclusions, were found ( Figure 3). Fluid inclusions are usually two-phase (vapor + liquid) at room temperature (25°C), except for the crystal-bearing aqueous inclusions, and colorless, yellow, brown, gray, and black under transmission light.

Results
The aqueous inclusions have the highest abundances in the study area and are small in size (3 μm-12 μm). Usually, the volume of the liquid phase is larger than that of the vapor phase, and the gas-liquid ratio ranges from 3% to 10%. They are colorless and light brown under transmission light (Figure 3(a)).
The hydrocarbon-bearing aqueous inclusions are generally 3 μm-12μm in diameter. The volume of the liquid phase is usually larger than that of the vapor phase, and the gasliquid ratio (5%-20%) is higher than that of aqueous inclusions. Generally, it is light brown and brown under transmission light, and a few inclusions show weak blue-white fluorescence under ultraviolet light (Figures 3(b) and 3(c)).

Geofluids
The hydrocarbon inclusions are relatively large (5 μm-18 μm) with the gas-liquid ratio of 10% to 60%. They are mainly light brown and brown under transmission light. Hydrocarbon inclusions have the property that the fluorescence under ultraviolet light changes from brown to yellow to white to blue as the increasing maturity of organic matter [56][57][58]). Most hydrocarbon inclusions in the study area show bright yellow or blue-white fluorescence, indicating that their organic matter has a certain maturity (Figures 3(d)-3(f)).
The crystal-bearing aqueous inclusions, ranging in size from 3 μm to 12 μm, contain the crystals formed by cooling in aqueous containing a solute that is supersaturated. They are three-phases (vapor + liquid + solid) at room temperature with the gas-liquid ratio of 5%-10%. The liquid phase is colorless and light brown under transmission light, and the crystals are colorless and transparent cubes (Figure 3(g)).
The aqueous-carbonic inclusions are composed of gaseous CO 2 , liquid CO 2 , and aqueous at room temperature, and the boundaries are obvious. Its diameter is generally between 10 μm and 15 μm, the gas-liquid ratio ranges from 10% to 20%, and it is light brown and brown under transmitted light (Figure 3(h)).
Under the influence of tectonic stress, differential compaction and abnormal high pressure, the Upper Paleozoic reservoir in Ordos Basin experienced multistage of microfracturing [2,59,60]. Meanwhile, organic acids generated in the hydrocarbon generation process of coal measure source rocks can dissolve unstable minerals such as pyroclasts and feldspars, and generate large amounts of Si 4+ , providing material sources for the formation of quartz overgrowth [61]. Therefore, it is a common way to divide the phases of fluid inclusions in Ordos Basin based on whether the secondary fractures with fluid inclusions cut through quartz particles and its overgrowth [62][63][64][65]. According to the results of cathode luminescence scan ( Figure S2), the development degree of quartz overgrowth is generally low, most of which are quartz overgrowth with thin local distribution, and only a few of quartz particles are completely surrounded by overgrowth. Based on the above

Microthermometry
Characteristics. The homogenization temperature of aqueous inclusions associated with hydrocarbon inclusions is closest to the real stratum temperature at the time of hydrocarbon charging. Therefore, it is preferable to select such inclusions to study the history of gas charging. According to the homogenization temperature statistics of samples (Table 1), it is considered that the distribution characteristics of homogenization temperature among different wells are quite different and can be divided into 3 groups. Group 1, including LX-12 well, LX-17 well, LX-18 well, and LX-22 well, has homogenization temperature distribution of 70°C-160°C. Group 2, including LX-1 well, LX-2 well, and LX-4 well, has homogenization temperature distribution of 80°C-200°C. Group 3, including LX-21 well, LX-26 well, LX-28 well, and LX-101 well, has homogenization temperature distribution of 60°C-150°C. The isoline maps of average homogenization temperature and salinity were drawn, and the obvious regularity distribution was found. As shown in Figure 5(a), the homogenization temperatures of LX-1, LX-2, and LX-4 wells are higher than that of other wells, and decreases to both sides. The salinity distribution shows a "northwest to southeast decline" distribution pattern ( Figure 5(b)). The above research shows that the fluid in the reservoir has obvious heterogeneity and regularity difference during the period of natural gas charging.
Based on the above analysis results, the homogenization temperatures and salinity of fluid inclusions are further analyzed. The study area is divided into 3 zones from the Zijinshan magmatic pluton outward: inner zone, middle zone, and outer zone, which correspond to 3 types of homogenization temperature and salinity distribution characteristics, respectively ( Figure 6).
In the inner zone (Group 1), the homogenization temperatures are generally low (70°C-160°C), and there are early and later fluid inclusions, but the number of later fluid inclusions is small and their distribution is not obvious. The salinity ranges from 2% to 22% and increases with the increase of homogenization temperature.
In the middle zone (Group 2), the homogenization temperatures are generally higher (80°C-200°C), and there   7 Geofluids are early and later fluid inclusions. The salinity ranges from 4% to 23% and increases with the increase of homogenization temperature.
In the outer zone (Group 3), the homogenization temperatures are generally low (60°C-150°C), and there are early and later fluid inclusions. The salinity ranges from 2% to 22%. The data distribution is discrete and the regularity is poor. that the basin has experienced at least 4 major uplifting and denudations [66,67]. In the study area, for the first 3 uplift denuded thickness, the unconformity surfaces have been denuded due to the large uplift since the Late Cretaceous, which makes it difficult to recover. According to previous studies, the total denuded thickness of the first 3 denudations in the northeastern Ordos Basin is only about 300 m, which has little effect on the recovery of sedimentary burial evolution [68]. Therefore, the first 3 denudations thickness recovered by Chen et al. [68] is used as the fixed value, which is 80 m, 60 m, and 200 m in turn. The denuded thickness is caused by uplift, since the Late Cretaceous is mainly determined by the interval transit time through the study of mudstone compaction trend [40]. The results show that the denuded thickness since the Late Cretaceous ranged from   9 Geofluids 700 m to 1600 m (Figure 7), higher than that in the southwestern part of the basin, which is consistent with the strong tectonic uplift in the northeastern part of the basin since the Late Cretaceous. The denuded thickness decreases annularly from Zijinshan magmatic pluton to its surroundings, which can be attributed to the uplift of strata around Zijinshan caused by magmatic emplacement and volcanic eruption, making it more vulnerable to be denuded.

Burial-Thermal Evolution and Hydrocarbon Generation-
Based on data such as denuded thickness and drilling lithology, the Sweeney and Burnham [69] EASY% Ro model in PetroMod 1D software was used to simulate the thermal evolution of multiple wells. The specific method is to take the Ro and the homogenization temperature of fluid inclusions as the standard values, and to continuously adjust the paleo-heat flow until the simulated values of Ro and stratum temperature match the standard values, a reliable thermal evolution history can be obtained. Taking LX-4 well as an example, the Ro in Benxi Formation of LX-4 well is 2.32%, and the maximum homogenization temperature of aqueous inclusions in the Taiyuan Formation is about 195°C. With the above data as standard values, the related data in LX-4 well are simulated as shown in Figure 8, which matches the standard values, and the sedimentary burial-thermal evolution can be obtained (Figure 9). Different burial-thermal evolution often represents different hydrocarbon accumulation characteristics. Based on the above studies, the study area can be divided into 3 areas: the inner-magma baking area, the middle-anomal thermal area, and the outer-normal thermal area ( Figure 10). The inner-magma baking area is the area where was baked by the Zijinshan magmatic pluton. The highest geothermal temperature experienced is generally greater than 200°C   10 Geofluids and is very close to the pluton, which was strongly influenced by its uplift ( Figure S4). The middle-anomal thermal area refers to the area where the highest temperature experienced is higher (180°C-200°C) and was affected by a certain degree of uplift caused by the Zijinshan magmatic pluton. The outer-normal thermal area refers to the area where the stratum temperature was weakly affected by the Zijinshan magmatic pluton (generally not exceeding 180°C) and was less affected by the pluton uplifting effect.

Hydrocarbon Generation-Expulsion. PetroMod 1D
software was used to simulate hydrocarbon generationexpulsion of coal source rocks in Taiyuan Formation, including LX-17 well in the inner-magma baking area, LX-4 well in the middle-anomal thermal area, and LX-101 well in the outer-normal thermal area ( Figure 11). The simulation results show that the burial depth of coal source rocks in the Early Jurassic generally reached about 2100 m, and the equivalent vitrinite reflectance (R eq ) of coal in LX-4 and LX-17 wells near Zijinshan magmatic pluton was 0.57%-0.68%, which had reached the hydrocarbon generation threshold. At this stage, gas generation was low, and some liquid hydrocarbons may be generated from the H-rich exinite in the coal source rocks from C 2 b to P 1 s [70,71]. During the Middle Jurassic, the R eq of coal in LX-4 and LX-17 wells increased to 0.66%-0.80% (Figure 11). The thermal evolution of source rocks gradually entered a mature stage, and the intensity of hydrocarbon generation-expulsion increased relatively and maintained until the Late Jurassic, which is the first peak period of hydrocarbon generationexpulsion in the relevant areas. During the Late Jurassic, the R eq of coal in LX-101 well in the outer-normal thermal area reached 0.61%-0.63%, which was in the stage of lowmaturity and reached the hydrocarbon generation threshold.
During the Early Cretaceous, source rocks reached the second peak period of hydrocarbon generation-expulsion, which is also the strongest period of hydrocarbon generation-expulsion. Under the dual influence of rapid sedimentary burial and Yanshanian thermal events in the North China Plate, the maximum burial depth could reach about 3000 m, and the R eq of coal in the LX-101 well could reach 1.21%, reaching the degree of high maturity ( Figure 11). Influenced by baking induced by Zijinshan magmatic pluton, R eq of coal in LX-4 and LX-17 wells could reach over 2.4%, reaching the stage of over-mature, with the highest hydrocarbon generation-expulsion intensity.

Geofluids
analyze the gas accumulation periods, improper conclusions may be drawn. Therefore, this paper analyzed the gas accumulation periods of several wells separately.
The homogenization temperature distribution of fluid inclusions in LX-17 and LX-22 wells in the inner-magma baking area is shown in Figure 12(a). There is a homogenization temperature peak range in early fluid inclusions, while the abundance of later fluid inclusions is low and the peak value is not obvious. The temperature peak range of early fluid inclusions is 110°C-120°C. Combining with burial-  12 Geofluids thermal evolution, it is considered that there was a major natural gas accumulation event that began in the Middle Jurassic ( Figure 12(a)). The homogenization temperature distribution of fluid inclusions in LX-1 and LX-2 wells in the middle-anomal thermal area is shown in Figure 12(b). Both early and later fluid inclusions in the two wells have homogenization temperature peak range, the temperature peak range of early fluid inclusions is 120°C-130°C, the temperature peak range of later fluid inclusions of LX-1 well is 150°C-160°C, and the temperature peak range of later fluid inclusions of LX-2 well is 140°C-150°C. Combining with burial-thermal evolution, it is believed that the first period of natural gas accumulation was the Middle-Late Jurassic, and the second period of natural gas accumulation was the Early Cretaceous ( Figure 12(b)).
The homogenization temperature distribution of fluid inclusions in LX-28 and LX-101 wells in the outer-normal thermal area is shown in Figure 12(c). Both early and later fluid inclusions in the two wells have homogenization temperature peak range, the temperature peak range of early fluid inclusions is 90°C -100°C, and the temperature peak range of later fluid inclusions of is 120°C-130°C. Combining with burial-thermal evolution, it is believed that the first period of natural gas accumulation was the Late Jurassic, and the second period of natural gas accumulation was the Early Cretaceous (Figure 12(c)).

Gas Accumulation Characteristics. The vertical accumulation of natural gas in the Upper Paleozoic of the Ordos
Basin is generally consistent with the characteristics of "gas accumulation in the lower reservoirs is earlier than the upper reservoirs." Natural gas fills the overlying reservoirs from the underlying coal-bearing source rocks, and the filling will inhibit the diagenesis of the filling reservoirs [72][73][74], resulting in a large number of low-temperature inclusions and a small number of high-temperature inclusions in the lower reservoirs. However, through the analysis of the homogenization temperature distribution characteristics, it is considered that the Linxing area has its unique vertical gas accumulation characteristics in the Ordos Basin.
There are unique gas accumulation characteristics in the inner-magma baking area and middle-anomal thermal area. Taking LX-22 well in inner-magma baking area and LX-4 well in middle-anomal thermal area as examples (Figure 13(a)), the homogenization temperature frequency distributions in the two wells in low-temperature range of relative lower formations (such as P 1 t in LX-4 well and P 2 s in LX-22 well) are not higher than that in relative upper formations (such as P 2 x in LX-4 well and P 3 s in LX-22 well). The data characteristics of fluid inclusions in the upper and lower formations are similar, and even the frequency of fluid inclusions in the low-temperature range of the upper formations is higher than that in the lower formations.
The outer-normal thermal area is consistent with the gas accumulation characteristics of the Ordos Basin. Taking LX-21 and LX-101 wells in the outer-normal thermal area as examples, the gas accumulation consistent with the characteristics of "gas accumulation in the lower reservoirs is earlier than the upper reservoirs." The frequency distributions of homogenization temperature in two wells (Figure 13(b)) are obviously different from that in Figure 13(a). The frequency in the low-temperature range of the relative lower formations (such as P 2 x in LX-21 well and P 1 t in LX-101 well) of the two wells is higher than that in the relative upper formations (such as P 3 s in LX-21 well and P 2 x in LX-101 well).
In case of single well contingency, taking C 2 b, P 1 t, and P 1 s as the lower formations, P 2 x, P 2 s and P 3 s as the upper formations, the homogenization temperature distributions of fluid inclusions in the inner-magma baking area, the middle-anomal thermal area, and the outer-normal thermal area were counted, respectively (Figure 13(c)). The frequency in the low-temperature range of lower formations in the inner-magma baking area and the middle-anomal thermal area is not higher than that in upper formations. The homogenization temperature distributions of fluid inclusions in upper and lower formations are similar, and their gas accumulation periods should be similar. In the outer-normal thermal area, the frequency in the low-temperature range of the lower formations is higher than that in the upper formations, and its frequency is obviously higher in the temperature range from 80°C to 100°C, which indicates that certain amounts of natural gas had charged into the reservoir when the lower formations were relatively low in temperature.

Trapping
Pressure. The homogenization temperatures of fluid inclusions are generally lower than its trapping temperatures, which is more obvious for hydrocarbon inclusions. This phenomenon can be attributed to two main factors [57,72]: (1) the homogenization temperatures are obtained under atmospheric pressure, while the trapping pressures are usually higher; (2) affected by the composition of inclusions, such as gaseous hydrocarbons. The trapping pressures of hydrocarbon inclusions and its coeval hydrocarbon-bearing aqueous inclusions in LX-17 and LX-18 wells were simulated by using PVTsim 20 software ( Figure 14; Table 2). The results show that, on average, for hydrocarbon inclusions, the trapping temperature is 17.6°C higher than the homogenization temperature, and the trapping pressure is 6.2 MPa higher than the minimum trapping pressure; for hydrocarbon-bearing aqueous inclusions, the trapping temperature is 2.0°C higher than the homogenization temperature, and the trapping pressure is 5.8 MPa higher than the minimum trapping pressure. Therefore, for samples without hydrocarbon inclusions, the trapping temperature can be approximately obtained by homogenization temperature of hydrocarbonbearing aqueous inclusions plus 2°C, and the trapping pressure can be obtained by substituting the isochores equation.
Due to the lack of later fluid inclusions in four wells in the inner-magma baking area, there are not enough samples to simulate reservoir pressure evolution. To solve this problem, by fitting the trapping temperatures and pressures of known fluid inclusions in the study area, it is considered that there is a good linear relationship between them (R 2 = 0:993) ( Figure 15). Combining with the study of thermal evolution recovery in Section 4.3.1, the strata temperatures are substituted into the fitting relationship between trapping temperatures and pressures, and the trapping pressures can be approximately obtained. On the basis of the above research, the evolution law of reservoir pressure during the gas accumulation periods was further obtained (Figure 16).
The reservoirs in the three areas were overpressured during the natural gas accumulation periods, but the development processes are different. During the Early Jurassic, the Upper Paleozoic reservoirs in the study area were in the state of hydrostatic pressure or weak overpressure, and then, the reservoir pressure state began to differentiate.
For the inner-magma baking area (Figure 16(a)), since the Middle Jurassic, the pressure has risen from 27 MPa to 37 MPa with natural gas charging. Because of the alternate occurrence of multiple small-scale strata uplift and subsidence during this time, the burial depth fluctuated but did not change much, and the pressure coefficient increased obviously, from 1.2-1.3 to 1.6-1.8. During the Early Cretaceous, with the increase of hydrocarbon generation intensity of source rocks, the reservoir pressure increased sharply, reaching a maximum of about 68 MPa. Because of the rapid subsidence, the pressure coefficient increased slightly. When For the middle-anomal thermal area (Figure 16(b)), the reservoir pressure generally increased from 23 MPa-26 MPa to 30 MPa-35 MPa during the Middle-Late Jurassic, and the pressure coefficient was about 1.6. The development of overpressure was consistent with the first gas charging during this time. During the early Cretaceous, when the second period of natural gas accumulation occurred, the reservoir pressure increased sharply from 33 MPa to 57 MPa. At the maximum burial depth, the pressure coefficient reached about 1.9-2.0.
For the outer-normal thermal area (Figure 16(c)), the reservoir pressure in the Middle Jurassic generally reached more than 20 MPa, and the pressure coefficient generally breaks through 1.0, which is in a weak overpressure state. During the Late Jurassic, the first period of natural gas charging occurred, and the reservoir pressure increased to    Figure 17). The distribution characteristics of the 4 periods are similar, and the fluid potential shows the characteristics of "annular decline from the periphery of the Zijinshan magmatic pluton." The potential difference between high and low potential areas keeps increasing. Until the Early Cretaceous, due to the high thermal degree of source rocks in the area near Zijinshan magmatic pluton, the generated natural gas charged into the reservoir rapidly and massively, and the fluid potential energy rises rapidly, and the difference of fluid potential in the study area reached the maximum.
The above results show that the fluid in the reservoir tended to migrated from the area around Zijinshan magmatic pluton to the low potential area.

Influence of Zijinshan Pluton.
Zijinshan magmatic pluton has a significant impact on the gas accumulation in the study area, which is manifested by regional differences in fluid

Geofluids
properties, hydrocarbon generation-expulsion, gas accumulation periods, vertical gas accumulation characteristics, and reservoir pressure.
The heterogeneity of fluid in gas accumulation period is mainly manifested in the difference of homogenization temperature and salinity distribution characteristics of fluid inclusions, which is supposed to be closely related to deep fluid activities triggered by Zijinshan magmatic pluton. Deep hydrothermal fluids are characterized by high temperature, high salinity, and rich in CO 2 , and often move along faults and fractures [75,76]. Because the inner and middle zones are close to the Zijinshan magmatic pluton, and there are many faults and fractures [77], the salinity of stratum water was in a relatively stable state during the interaction between surface freshwater and the upwelling of deep high salinity fluids. The salinity of fluid inclusions was mainly controlled by the increase of salinity solubility caused by the rise of strata temperatures. The interaction between strata water and surface freshwater was weak in the outer zone due to the lack of faults and fractures. The salinity of stratum water  T t : trapping temperature; P t : trapping pressure; T h : homogenization temperature; P min : minimum trapping pressure; G/L: gas-liquid ratio. 18 Geofluids was controlled by the frequent activities of deep fluids and has a poor correlation with the temperature [78][79][80].
The whole gas accumulation period in the Linxing area is from the Middle Jurassic to the Early Cretaceous, but there are differences among the 3 areas. In the inner-magma baking area, there was a major natural gas accumulation event that began in the Middle Jurassic, earlier than the other two areas. There were 2 periods of gas accumulation in the middle-anomal thermal area, but the earliest time of gas accumulation is the Middle-Late Jurassic, which is later than the inner-magma baking area. There were 2 periods of gas accumulation in the outer-normal thermal area, the first occurred in the Late Jurassic.
The fundamental factor that causes the difference in the accumulation periods of the study area is the thermaltectonic effect induced by the Zijinshan magmatic pluton. The inner-magma baking area was most strongly influenced by thermal baking and uplifting, and the short uplifting events in Yanshanian period cannot offset the thermal baking. Source rocks rapidly reached the degrees of highmaturity or even over-maturity, with high gas generation intensity and short interval of episodic hydrocarbon expulsion, and long-lasting gas charging and accumulation. The microfractures in coal seams and reservoirs caused by uplift also promoted the gas charging from source rocks to sandstone reservoirs. The lack of later high-temperature inclusions in the inner-magma baking area is also due to the rapid increase of gas saturation in reservoirs, the strata water was rapidly displaced, and the formation of later inclusions was inhibited [72][73][74], which also reflects the high gas generation in inner-magma baking area. The middleanomal thermal area suffered a certain degree of thermal baking and uplifting, which resulted in the thermal evolution of source rocks was higher than that of the basin, and the accumulation period is later than that of the inner-magma baking area. Under the influence of Zijinshan magmatic pluton, gas accumulation in the inner-magma baking area and the middle-anomal thermal area was not a "step-by-step" process from bottom to top. A large amounts of natural gas rapidly generated by the lower source rocks quickly reached the upper reservoirs through faults and fractures associated with the pluton. The interval time between the upper and lower formations is short and the accumulation time is early.
The thermal baking and uplifting effects in the outernormal thermal area were weak, and the source rock matured slowly, and the reservoirs were filled until the Late Jurassic. Affected by "episodic hydrocarbon expulsion effect" and short-term uplift events in Yanshanian period, the first period of hydrocarbon expulsion tended to stagnate, and the maturity of source rocks increased further until the Early Cretaceous rapid subsidence stage, resulting in the second period of hydrocarbon generation-expulsion. Because of the small thermal baking and uplifting, and the natural gas accumulation in this area consistent with the characteristics of "gas accumulation in the lower reservoirs is earlier than the upper reservoirs." The overpressure degree in the inner-magma baking area, the middle-anomal thermal area, and the outernormal thermal area decrease in turn, reflecting that thermal baking of Zijinshan magmatic pluton controlled the thermal evolution and hydrocarbon generation process of source rocks, and further affects the distribution of reservoir pressure. Besides, the pressurization caused by deep highpressure fluid and the gases released from the volatiles in the magma cannot be ignored [81,82].

Natural Gas Accumulation Events.
In the Linxing area, favorable source-reservoir-cap assemblages have been formed in the long process of geological and sedimentary evolution. Natural gas accumulation occurred during the Jurassic-Cretaceous, when the upper overburden rock was thick and dense, traps had been formed or were being formed in large quantities, which is a beneficial space-time matching relationship with the migration and accumulation of natural gas ( Figure 18).
According to the different influence degree of the thermal-tectonic effect caused by Zijinshan magmatic pluton, natural gas reservoirs in the 3 areas were gradually formed. The hydrocarbon generation of source rocks in the innermagma baking area started the earliest, began in the Early Jurassic, and reached the first peak of hydrocarbon generation-expulsion in the Middle Jurassic. Natural gas flooded into the reservoir in large quantities. Among them, free gas migrated along the structural updip, and dissolved gas migrated to the low potential area with strata water, which accumulated and formed reservoirs in favorable traps.

Geofluids
The source rocks in the middle-anomal thermal area also began to generate hydrocarbons in the Early Jurassic, but later than the inner-magma baking area, and the natural gas charged into the reservoirs in large quantities in the two-period accumulation of the Middle-Late Jurassic and Early Cretaceous. Source rocks in the outer-normal thermal area began to generate hydrocarbons until the Late Jurassic, and the first period of natural gas accumulation occurred. The second period of gas accumulation occurred during the Early Cretaceous.
From the Late Cretaceous to the Neogene, the Linxing area entered the period of overall uplift and gas reservoir adjustment-loss. Zijinshan magmatic pluton also rose rapidly in this period and triggered a large number of deep and large faults. Li et al. [83] calculated the gas content of sandstone reservoirs in the Linxing area according to Indonesian Formula, which shows that there exists a low gas-bearing or nongas-bearing zone around Zijinshan magmatic pluton. According to the gas reservoir profile constructed by the logging data interpretation (Figure 19), it is found that at  20 Geofluids present, the reservoirs in the wells (such as LX-3 well, LX-17 well, and LX-18 well) near Zijinshan magmatic pluton are mostly water layers. It is generally believed that a large number of gas reservoirs in the inner-magma baking area were destroyed by deep and large faults, and a large number of natural gas were lost. Gas reservoirs in thermal anomaly area and normal thermal evolution area have been lost in the process of uplift, but there are fewer faults and gas reservoirs have been preserved up to now. Although the gas reservoirs in the middle-anomal thermal area and the outer-normal thermal area have a certain amount of loss during the uplift process, there are fewer faults, and the gas reservoir has been preserved until now.

Natural Gas Accumulation
Model. According to the results of hydrocarbon generation-expulsion simulation, gas 21 Geofluids accumulation periods, reservoir pressure evolution, and fluid potential analysis, and combined with the characteristics of present gas reservoirs, it is considered that there is an accumulation model of the Upper Paleozoic tight gas reservoir in the Linxing area, which the accumulation process is affected by thermal-tectonic effect caused by the unique magmatic intrusion-the natural gas accumulation model under the control of magmatic thermal-tectonic effect.
The natural gas accumulation process in the Linxing area can be divided into 3 stages. In the first stage, the Early Jurassic, some of the source rocks near the Zijinshan magmatic pluton began to generate small amounts of natural gas and charged into the adjacent reservoirs. At this stage, the natural gas migration driving force was insufficient, and it was impossible to move laterally or vertically over long distances.
The second stage, from the Middle Jurassic to the Early Cretaceous, is the key period of gas generation, migration, and accumulation. The gas accumulation control mechanism of Zijinshan magmatic pluton is mainly reflected in its thermal baking and reservoir-transport system transformation. Thermal baking first affected the maturity evolution of source rocks in the inner-magma baking area and the  Figure 18: Petroleum system event in the Linxing area. 22 Geofluids middle-anomal thermal area, which led to the large-scale hydrocarbon generation-expulsion stage in the Middle Jurassic earlier than the outer-normal thermal area. Moreover, the hydrothermal pressurization and hydrocarbon generation pressurization caused by thermal baking not only make the gaseous stratum water and natural gas migrate in a mixed state, greatly reducing the migration resistance, but also make the inner-magma baking area in a long-term high fluid potential center, which provides the driving force for natural gas migration. In the aspect of reservoir-transport system transformation, the faults and fractures caused by the intrusion of Zijinshan magmatic pluton provide a good migration channel and reservoir space for natural gas, so that natural gas can be filled to the shallow Shangshihezi Formation and Shiqianfeng Formation reservoirs ( Figure 20). In the third stage, from the end of Early Cretaceous to Neogene, the whole study area continued to uplift, and the gas accumulation control mechanism of Zijinshan magmatic pluton is mainly reflected in the gas reservoir being destroyed by faults. After severe uplift and denudation, the gas reservoir loses its overlying thick overburden rocks. Because of a large number of deep normal faults in the inner-magma baking area, a large number of natural gas is lost, and the gas reservoir was adjusted and destroyed ( Figure 20). There are no deep faults in the middle-anomal thermal area and the outer-normal thermal area, and the gas reservoirs are well preserved, which are the favorable areas for natural gas exploration and development.

Conclusions
(1) The fluid heterogeneity in the Linxing area was strong during gas accumulation periods. According to the statistical analysis of homogenization tempera-ture and salinity of fluid inclusions, combined with the reconstructed burial-thermal evolution, the study area was divided into 3 areas: the inner-magma baking area, the middle-anomal thermal area, and the outer-normal thermal area (2) Under the influence of Zijinshan magmatic intrusive tectonic-thermal event, the gas accumulation periods are different in different areas. In the inner-magma baking area, there was a major natural gas accumulation event that began in the Middle Jurassic, earlier than the other two areas. There were 2 periods of natural gas accumulation in the middle-anomal thermal area. The first period occurred from the Middle-Late Jurassic, later than the magma baking area, and the second period occurs in the Early Cretaceous. There were also 2 periods of gas accumulation in the outer-normal thermal, the first period of gas accumulation was in the Late Jurassic, and the second period occurred in the Early Cretaceous (3) The vertical gas accumulation characteristics of natural gas show obvious regional differences. Due to the thermal-tectonic effect caused by Zijinshan magmatic pluton, natural gas accumulation in the inner-magma baking area and the middle-anomal thermal area was not a slow and gradual process from bottom to top. A large amounts of natural gas rapidly generated by the lower source rock rapidly reached the upper formation through faults and fractures caused by the intrusion of the Zijinshan magmatic pluton. The time interval between the upper and lower formations was shorter. The outer-normal thermal area was less affected by the thermal-

Gas layer Bad gas layer
Dry layer Coal seam Water layer Figure 19: N-S gas reservoir profile in the Linxing area. 23 Geofluids tectonic effect of Zijinshan magmatic pluton, and consistent with the characteristics of "gas accumulation in the lower reservoirs is earlier than the upper reservoirs." (4) The trapping pressure of fluid inclusions was simulated by PVTsim 20 fluid simulation software; it is shown that the Upper Paleozoic reservoirs had obvious overpressure characteristics in the gas accumulation  Figure 20: Schematic diagram of gas accumulation process and model of tight gas reservoir in the Linxing area. 24 Geofluids periods. The period from the Middle Jurassic to the Early Cretaceous is the key period for rapid pressure accumulation, which is consistent with the period of natural gas accumulation (5) The study on the distribution and evolution of paleo-fluid potential shows that the area near the Zijinshan magmatic pluton was the high fluid potential area during the gas accumulation periods, which indicates that natural gas and other fluids migrated from Zijinshan magmatic pluton to the surrounding area (6) The natural gas accumulation model under the control of magmatic thermal-tectonic effect in the Linxing area was proposed. The accumulation process can be divided into 3 stages. In the first stage, during the Early Jurassic, source rocks close to the Zijinshan magmatic pluton began to generate a small amount of natural gas and charged into the adjacent reservoirs.
In the second stage, the Middle Jurassic to the Early Cretaceous was the key period for gas migration and accumulation. The control mechanism of Zijinshan magmatic pluton was mainly reflected in its thermal baking and reservoir-transport system transformation. In the third stage, from the end of Early Cretaceous to Neogene, the Linxing area continued to uplift, and the formed gas reservoirs were further adjusted and destroyed due to a large number of deep faults triggered by Zijinshan magmatic pluton

Data Availability
The data used to support the findings of this study are included within the article. Table S1: sedimentary facies and paleo-water depths in the northeastern Ordos Basin. Table S2: coal sample parameters for basin modeling in the Linxing area. Figure S1: evolution of paleo-heat flow in the northeastern Ordos Basin. Figure  S2: photographs of cathode luminescence scan in the Linxing area. [(a) and (b) are 2 cathode luminescence photographs, the clastic composition is mainly quartz which glows blue and dark-blue light, and a few glows reddish-brown light; the quartz overgrowth does not glow with weak development, only locally developed in a few particles; calcite glows orange light. (c) and (d) are plane polarized and cathode luminescence photographs of the same sample, respectively, the quartz overgrowth is rarely developed. (e) and (f) are plane polarized and cathode luminescence photographs of the same sample, respectively, the quartz overgrowth is locally developed. (g) and (h) are plane polarized and cathode luminescence photographs of the same sample, respectively, the quartz overgrowth is relatively developed, and the broken and healed phenomenon of quartz particle can be seen. (i) and (j) are plane polarized and cathode luminescence photographs of the same sample, respectively, the quartz overgrowth is locally developed, and the broken and healed phenomenon of quartz particle can be seen. (k) and (l) are plane polarized and cathode luminescence photographs of the same sample, respectively, the quartz overgrowth is locally developed, and plagioclase glows yellow-green light.]. Figure S3: schematic diagram of the petrography of the early and later inclusions. Figure S4: 1200 ms seismic isochronic section in the Linxing area (cited from [14]). (Supplementary Materials)