Multi-Phases Fluid Activity Characteristics of Longmaxi Formation and Its Impact on Resistivity in the Changning Area, Southern Sichuan Basin, Southwest China

Wells with low gas content and low resistivity in the Changning area, southern Sichuan Basin were selected for this study. The burial-thermal history was reconstructed and the characteristics of multi-phase fluid activity were clarified using microscopic observation and testing of fluid inclusions in the Longmaxi shale fracture veins. Compared with wells with a high gas content and high resistivity, the influence of fluid activity on resistivity was analyzed. The results showed that the thermal evolution of the bituminous inclusions trapped in the veins has reached the stage of carbonaceous-metamorphic bitumen, and the organic matter is fully cracked for gas generation, with some organic matter exhibiting the phenomenon of “graphitization.” The synchronous fluid with bitumen was existed due to shallow burial with a middle-low maturity stage of about 280 and 292 Ma. Two phases of fluids existed in the deep burial stage, thus maturing for about 103 Ma, and the uplift stage at about 28 and 32 Ma, with high homogenization temperatures (Th) (varying from 185 to 195°C and 165–180°C). The corresponding pressure coefficients varied between 1.67 and 2.09, 1.56 and 1.92 in a moderate-strong high-pressure state. The last two phases of fluid formation in the late uplift stage for about 4 to 19 Ma and 6 to 10 Ma were characterized by low salinity at medium-low Th (varied from 140 to 155°C and 120–135°C), with pressures of 57.47–74.50 MPa and 51.44–59.41 MPa (pressure coefficients of 1.09–1.41 and 1.18–1.37), in an atmospheric-weak overpressure state. In the initial uplift stage after deep burial, the fluid closure in the Changning area was good. In contrast, the wells are filled with low gas content because of the strong tectonic forces causing the shale gas to be released and the multi-phase fluid activity that happens during the late uplift stage. New evidence indicates that the emergence of low resistance in the localized Changning area is not only related to the high degree of evolution of organic matter but is also affected by the multi-phase fluid modification in the late uplift stage.

Wells with low gas content and low resistivity in the Changning area, southern Sichuan Basin were selected for this study. The burial-thermal history was reconstructed and the characteristics of multi-phase fluid activity were clarified using microscopic observation and testing of fluid inclusions in the Longmaxi shale fracture veins. Compared with wells with a high gas content and high resistivity, the influence of fluid activity on resistivity was analyzed. The results showed that the thermal evolution of the bituminous inclusions trapped in the veins has reached the stage of carbonaceous-metamorphic bitumen, and the organic matter is fully cracked for gas generation, with some organic matter exhibiting the phenomenon of "graphitization." The synchronous fluid with bitumen was existed due to shallow burial with a middle-low maturity stage of about 280 and 292 Ma. Two phases of fluids existed in the deep burial stage, thus maturing for about 103 Ma,and the uplift stage at about 28 and 32 Ma, with high homogenization temperatures (Th) (varying from 185 to 195°C and 165-180°C). The corresponding pressure coefficients varied between 1.67 and 2.09, 1.56 and 1.92 in a moderate-strong high-pressure state. The last two phases of fluid formation in the late uplift stage for about 4 to 19 Ma and 6 to 10 Ma were characterized by low salinity at medium-low Th (varied from 140 to 155°C and 120-135°C), with pressures of 57.  MPa (pressure coefficients of 1.09-1.41 and 1.18-1.37), in an atmospheric-weak overpressure state. In the initial uplift stage after deep burial, the fluid closure in the Changning area was good. In contrast, the wells are filled with low gas content because of the strong tectonic forces causing the shale gas to be released and the multi-phase fluid activity that happens during the late uplift stage. New evidence indicates that the emergence of low resistance in the localized Changning area is not only related to the high degree of evolution of organic matter but is also affected by the multi-phase fluid modification in the late uplift stage.

INTRODUCTION
The southern Sichuan Basin in China is a typical demonstration area with remarkable achievements in marine shale gas exploration and development. Shale gas is mainly enriched in the Upper Ordovician (Wufeng formation)-Lower Silurian (Longmaxi formation) shales (Ma and Xie, 2018;Yang et al., 2019;Gou et al., 2021;Han et al., 2022). Nevertheless, with relatively similar initial sedimentary conditions, exploration wells with low or even no gas content have emerged in some high gas content areas (Bachrach, 2011;Wang et al., 2014;Yang et al., 2015;Xu et al., 2020a;Xu et al., 2020b;Jia et al., 2021). The organic shale sections in the failed exploration wells often exhibit extremely low resistivity logging, with logging resistivities below 10 Ωm, and some even below 1 Ωm. In most studies related to shale gas, attention has been given to this special phenomenon (Hou et al., 2021).
Since the time of deposition, marine shales in the Changning area, southern Sichuan Basin have generally undergone the modifications of deep burial, diagenesis, and hydrocarbon generation. However, multi-stage and complex tectonic processes have led to different degrees of fluid activities and tectonic deformation among different tectonic units, resulting in significant differences in shale gas preservation conditions and huge differences in shale gas shows or production of exploration and development wells Nie et al., 2020), with lowproducing wells often accompanied by low resistivity characteristics. Several studies have been conducted about high clay content, very mature organic matter, high pyrite content, high water salinity, poor physical properties in the reservoir and excellent geological conditions which cause abnormally low resistance in the reservoir (Sun et al., 2018;Hou et al., 2021). However, the influencing factors vary in different regions and reservoirs. Most researchers believe that the "carbonization" of organic matter is the main cause of the phenomenon of low resistivity in the Changning region (Wang et al., 2014;Wang et al., 2018;Hou et al., 2021), but some researchers still hold different views and believe that the mosaic of rock particles in close contact with each other is the main cause (Gao et al., 2016). The above studies are mostly based on the properties of the reservoir rock itself. Low resistivity as a parameter for characterizing shale gas content needs to be further researched from the perspective of the relationship between the changes in shale gas content caused by fluid activities in the reservoir as well as low resistivity. As the carrier of original information regarding geological fluids, fluid inclusion (FI) has become an important means of studying paleo-fluid activities (Chen et al., 2002;Liu et al., 2010;Xu et al., 2020c;Wang et al., 2020).
To clarify the influence of fluid activity on gas bearing property and resistivity, wells with extremely low gas content and resistivity in the Changning area were selected for this study. By using microscopic observation, homogenization temperature (Th) measurement of fluid inclusion assemblages (FIAs), and Raman spectroscopy testing, this study focuses on analyzing the tectonic evolution characteristics of low-producing wells, defines the multiphases fluid properties and activities, and compares and analyzes the fluid activity characteristics of previously studied wells with high gas content and resistivity.

GEOLOGICAL SETTING
Sichuan Basin is located in Southwest China, east of the Qinghai Tibet Plateau. It belongs to the first-class tectonic-sedimentary basin in the northwest of the Upper Yangtze platform. The Longmenshan fold orogenic belt is situated in its northwest direction; in the north, there is Micangshan-Dabashan fold orogenic belt; the Qiyueshan fold belt is in the southeast and the Daloushan-Daliangshan fold belt is in the south, creating an overall rhombic NE trend. The Changning area is located on the low and steep deformation belt of the southern margin of Sichuan Basin and the Daloushan-Daliangshan fold belt ( Figure 1A). Bounded by the Changning anticline, the northern part belongs to the low and gentle fold belt of southern Sichuan Basin, and the Daloushan-Daliangshan fold belt is in the south . The study area has mainly existed east-west, northwest and northeast structures, including Xuyong, Jianwu, Luochang, and Shuanglongchang synclines and the Changning anticline. The folded shape is relatively wide and gentle ( Figure 1B).
The study area was subjected to tectonic transformation in many stages and different directions during its geological history Liu et al., 2021). The early Caledonian movement led to the extension of the Craton Basin in southern China, the development of large tensile faults in the basement, and the extensional tectonic stress in the Changning area. The strong Dongwu uplift caused the Qianzhong paleouplift to transport from south to north, and the whole Changning area was in the uplift stage. After entering the Indosinian tectonic period, the Changning area received stable sedimentation. The paleouplift ceased with the north-south tectonic form being laid down. During the Yanshanian and Himalayan periods, the Changning area was compressed by the Dabashan in the east to north-south, and the structure was deformed into an east-west structure in the south of the Changning area, finally evolving into the present-day fault-fold belt (Feng et al., 2022). The multi-stage structural superposition causes the reverse faults that existed in the Changning area, and their scale is mainly small-medium faults, most of which disappear within the Silurian strata ( Figure 1C).
The surface of the study area is covered by sedimentary rocks, except for the Permian Emeishan basalt, which is outcropped sporadically in the southwest. The overall stratigraphic sequence is normal, lacking the Carboniferous and Devonian layers, and the exposed strata on the surface have existed from the Cambrian to Paleogene periods. The Silurian is denuded in the core of the Changning anticline, and the Cambrian strata are exposed ( Figure 1B).
The buried depth of the Longmaxi formation in the Changning area is mainly between 1,500 and 4,000 m. The burial depth is generally less than 3,500 m in the high part of the Changning structure and the south flank, and gradually increases to the north and west. The Longmaxi formation received marine finegrained sedimentation and existed members 1 and 2 from bottom to top. Among them, sub-member 1 of member 1 mainly existed in deep-water shelf facies organic matter rich shale, with a high organic carbon content (>3%). The degree of thermal evolution is in the stage of high and over maturity, experiencing complex thermal evolution processes such as oil generation, oil cracking gas, and kerogen cracking gas (Ma and Xie, 2018).

SAMPLES AND METHODS
When considering the differences in resistivity and gas production of different exploration wells in the Changning area, well N228 located in the Luochang syncline with extremely low resistivity and almost no gas production was selected for the research. This selection was based on the collection and summary of the fracture veins of the shales studied by previous researchers (Gong, 2019;Liang et al., 2021). Five core samples from five veins of the Longmaxi shales at depths ranging from 3,456.3 to 3,479.4 m were tested ( Figure 2). The Luochang syncline where well N228 occurs is in the Changning deep depressional area, and the buried depth of the Longmaxi shales in this area is more than 3,500 m. Wells N209, N215, and N216 which were studied by previous researchers are in the Jianwu syncline, in the shallow depression of the Changning region. The buried depth of Longmaxi shales is between 2,000 and 3,500 m. The bottom depth of the Longmaxi formation of well N228 is 3,480 m whereas the bottom depths of the Longmaxi formation of wells N209, N215, and N216 in Jianwu syncline are 3,170, 2,950, and 2,600 m, respectively.
Fluid inclusion analysis was performed after the thin sections of the shale fracture vein sample were doubly polished. The characteristics, types, and occurrence minerals with different generations of fluid inclusions were observed by transmission polarized light microscopy, cathode luminescence microscopy, fluorescence microscopy, and Raman spectral. All experiments were conducted in the Experimental Research Center, School of Geosciences, Yangtze University, Wuhan, China.
The experimental equipment used for transmission polarization and fluorescence is a Nikon Eclipse 80i dual-channel fluorescencetransmission light microscope with a UV excitation wavelength of 330-380 nm. The equipment used for cathode luminescence is the CL8200 MK5 cathode luminometer. Under the conditions of 14-15 kV and 300-350 mA, it emits light onto the thin-section to determine the formation period of the veins. Hydrocarbon inclusions usually emit different colors of fluorescence under UV excitation. The organic molecular structure in the inclusions controls the color and intensity of fluorescence, reflecting the characteristics of the mixture (Barres et al., 1987;Pradier et al., 1990;Permanyer et al., 2002;Li et al., 2004). For the gas inclusions in these samples, a LabRAM HR800 laser Raman spectrometer with a YAG solid-state laser at a wavelength of 532 nm was used for non-destructive characterization of individual inclusions larger than 1 μm to verify the gas-phase components. A different period of aqueous inclusions was selected to determine the Th, and the measured minimum Th of FIA was used to approximate the formation temperature at the time when oil and gas inclusions were trapped. Based on the homogeneous method and the freezing method, the Th and the ice melting temperatures (Tm) were determined by using the Linkam THMSG600 microscope cold-heat table, and the error of the temperature measurement results was controlled to about ± 0.1°C. The Th of methane inclusions can be determined by rapidly cooling the inclusions to the phase change point using liquid nitrogen, continuing cooling the inclusions to observe bubble changes, and then increasing the temperature to reach a homogeneous phase of the methane inclusions and recording the Th.

Vein and Mineral Characteristics
The Longmaxi shales are mostly buried at depths greater than 3,500 m in the Luochang syncline, located in the deep depressional area of the Changning area, whereas the Jianwu syncline is in the shallow depression of the Changning area where the Longmaxi shales are buried at depths between 2,000 and 3,500 m. The Longmaxi shales of well N228 are black and grayishblack with foliation formation. Many natural fractures can be observed in the core, some of which are filled with calcite, mainly in a high-angle attitude and some in a horizontal attitude. The horizontal vein is wider than the high angle vein, with the width of the high angle vein being approximately 0.2 cm, and that of the horizontal vein being approximately 1 cm. High-angle calcite veins mostly appear in groups, and some can develop in conjugate X-type formations ( Figure 2).
For macroscopic observation, the shale veins were sliced and microscopic petrographic observation was conducted to clarify the mineral composition, symbiotic relationship, formation phase, and other information. It is evident from the results that the calcite in the direct contact part between the edge of the vein and the surrounding rock is irregular and the cleavage is not obvious. According to the crystal size, it is considered to be fine crystal calcite ( Figure 3C). The calcite crystals become larger toward the formation direction inside the vein, which is coarse crystal calcite, and the surrounding rock fragments are a mixture of calcites ( Figure 3A). The cleavage inside the crystal is obvious and more than 1-2 groups have existed. The edges of some coarse crystal calcites are rough and existed in a zigzag shape. In addition, calcite veinlets connected with the main vein have existed in some surrounding rocks ( Figure 3C). For ascertaining the different phases of calcite, cathodoluminescence detection of different veins was conducted since it would be affected by Fe 2+ and Mn 2+ in the fluid (Pierson, 1981;Ten-Have and Heijnen, 1985). Some calcite veins emit overall dark cathodoluminescence, and the cathodoluminescence of sample N228-3463.5 m is enhanced from the surrounding rock to the interior of the vein ( Figures  3A,B). Comparing the cathodoluminescence of different calcite crystals from different samples, it was found that some coarse crystal calcites with brightly colored cathodoluminescence may have high Fe 2+ and Mn 2+ content (Huang, 1992;Huang et al., 2008). Different cathodoluminescence indicates that the fluid composition changed during the crystal growth. Additionally, combined with the observation results of cathodoluminescence and transmitted light, it was found that coke bitumen existed simultaneously at the edge of the crystal, while fracture formation was visible inside the vein ( Figures 3C,D).

Petrographic Characteristics of Fluid Inclusion
Fluid inclusions can be divided into primary, secondary, and pseudo-secondary inclusions according to their genesis. They can also be divided into gas, liquid, solid and multiphase inclusions according to the phase state, and can be divided into aqueous, hydrocarbon, CO 2 inclusions, etc., (Liu, 2005). In this study, the fluid inclusions are existed in euhedral and semi-euhedral calcite crystals. Hydrocarbon inclusions are mainly methane inclusions in a single gas phase, methane-bearing bituminous inclusions in solid-gas phases, and bituminous inclusions in a single solid phase ( Figure 4). The symbiotic relationship between them and aqueous inclusions can be classified into the following six types.
The first type is gas-liquid two-phase aqueous inclusions, which are linear, clustered, densely distributed, or existed separately, with a quadrilateral or elliptical morphology and a long axis of 3-20 μm. The gas-liquid ratio is of two kinds, one is larger, between 15% and −30%, and the other is smaller, between 5% and −10%. The second is linear methane inclusions clustered, densely distributed, or existing separately, with a quadrilateral, elliptical or irregular shape and a long axis between 3 and −30 μm. The methane inclusions with highdensity characteristics are uniformly translucent-transparent under the transmitted light, bright in the middle and surrounded by light gray, while other methane inclusions are less transparent and gray in transmitted light (Liu et al., 2009). Not all methane inclusions are fluorescent ( Figures 4A,B). The third is the symbiosis of gas-liquid two-phase aqueous inclusions and methane inclusions. This combination is referred to as a group of FIA in micro-temperature measurement, in which the methane inclusions are densely existed in lines or clusters, with a regular shape (quadrilateral or elliptical), and the long axis being between 3 and −15 μm. The fourth is bituminous inclusions, whose size and distribution characteristics are like those of the second type, which are grayblack or black under transmitted light, opaque and nonfluorescent ( Figures 4C,D). The fifth is the symbiosis of gas-liquid two-phase aqueous inclusions and bituminous inclusions, which is called a group of FIA in micro-temperature measurement. The size and distribution characteristics of these are like those of the third type (Figures 4E,G,H). The sixth is methane-bearing bituminous inclusions. The bitumen is attached to the sidewalls of the inclusions with methane in the center. Under transmitted light, the inclusions have one or more black lines on the sidewalls and are brighter in the middle, but remain gray-black and non-fluorescent, with similar size and distribution characteristics to the second type of methane inclusions.

Methane Inclusion and Methane-Bearing Bituminous Inclusion
Raman spectroscopy allows qualitative analysis of the composition of inclusions as well as the host mineral composition ( Figure 5). The Raman scattering peak shifts of methane inclusions measured in this experiment have a wide range of displacement, between 2,909.07 cm −1 and 2,915.49 cm −1 ( Figure 5B). Since the Raman scattering peak of pure methane inclusion is closely related to the internal pressure, the methane Raman scattering peak deviates to 2,910.00 cm −1 which indicates supercritical high-density inclusion, and that deviates to 2,918.00 cm −1 which usually indicates lowpressure gas-phase methane (Lu et al., 2007), with generally less than 2,912.00 cm −1 being supercritical state high-density methane inclusions (Liu et al., 2009).
Meanwhile, the methane Raman scattering peak of methanebearing bituminous inclusions are like those of pure methane inclusions, ranging from 2,907.05 cm −1 to 2,917.14 cm −1 . It also shows a pair of bitumen Raman characteristic peaks, with Raman shifts near 1,355 cm −1 (D-peak) and 1,580 cm −1 (G-peak) ( Figure 5A).

Over Mature Bituminous Inclusion
The Raman spectrum of the bituminous inclusions trapped in well N228 shows a pair of strong carbonaceous bitumen Raman characteristic peaks, with Raman shifts near 1,355 cm −1 (D-peak) and 1,580 cm −1 (G-peak). According to the variation pattern of D-peak and G-peak heights in the thermal evolution of bitumen obtained by Liu et al. (2013), the bituminous inclusions in the study area can be divided into three types. The D peak height (Dh) of the first type of inclusions is significantly lower than the G-peak height (Gh), and the shape of the G-peak is relatively obvious, which is the stage of carbonaceous bitumen ( Figure 5C); The Dh of the second type of inclusion is like the Gh, which is the stage of carbonaceous-metamorphic bitumen ( Figure 5D); the Dh of the third type of inclusion is significantly higher than the Gh, and the shape of the D-peak is relatively obvious, and this is the stage of metamorphic bitumen ( Figure 5E).
Meanwhile, the Raman spectrum can reflect the atomic and molecular vibration information in the structure of aromatic carbon rings of organic matter, that is, the spectral peak parameters of the Raman spectrum of organic matter have a good response relationship with the changes in the chemical structure of organic matter during thermal evolution, and thus characterize the maturity (Kelemen and Fang, 2001;Zhou et al., 2014). Vitrinite reflectance of kerogen is the most intuitive parameter to characterize the maturity of organic matter (Hunt, 1996), which is not contained in the highly mature hydrocarbon source rocks of the Lower Paleozoic, and therefore cannot be used to evaluate the maturity of its hydrocarbon source rocks. However, solid bitumen is the residue left after oil cracking, and the peak spacing between G-and D-peaks has been proven to be the most reliable maturity index (Kelemen and Fang, 2001;Liu et al., 2013;Zhou et al., 2014). Therefore, the bitumen reflectance can be converted into vitrinite reflectance to evaluate the maturity of Lower Paleozoic hydrocarbon source rocks. Use the following formula to calculate the same (Liu et al., 2013): where Rmc Ro% is the reflectivity calculated by Raman analysis parameters, %; d (G-D) is the peak spacing of Raman shift. According to Formula 1, we calculated the equivalent vitrinite reflectance of shale bitumen as well as bituminous inclusions in the Longmaxi Formation in the Changning area, which ranges from 2.73% to 3.91%. In contrast, the peak spacing is no longer considered to change at the over-mature stage observed in previous studies, which is only suitable for the calculation of high-over mature bituminous inclusions. In contrast, the peak height is suitable for the calculation of the maturity of solid organic matter from the over-mature to the pre-graphite stage. Therefore, Liu et al. (2013) proposed a new formula: Rmc Ro% 1.1659h(Dh Gh) + 2.7588 (2) where Rmc Ro% is the reflectivity calculated by Raman analysis parameters, %; h (Dh/Gh) is the peak height ratio (Dh is the D-peak height and Gh is the G-peak height). According to Formula 2, the maturity ranges from 3.54% to 4.05%. The results calculated by the two methods are similar, indicating that the maturity of organic matter in the study area is extremely high and reaches the stage of carbonaceousmetamorphic bitumen.

Th and Salinity of Symbiotic Aqueous Inclusions
Single-phase methane inclusions or bituminous inclusions are commonly found around the gas-liquid two-phase aqueous inclusions in the vein, indicating that these types of inclusions were formed in an immiscible two-phase system with saturated hydrocarbons. The Th and salinity of the gas-liquid two-phase aqueous inclusions represent the capture temperature and salinity of hydrocarbon inclusions. Therefore, the gas-liquid two-phase aqueous inclusions symbiotic with hydrocarbon inclusions FIA were selected for micro-temperature measurement in this study. The measured Th and salinity are shown in Figures 6, 7 and Since some of the inclusions may be metastable, not all aqueous inclusions can be converted to salinity by measuring the Tm (salinity calculation based on Bodnar, 1993). Therefore, the salinity of the aqueous inclusions in samples N228-3456.3 m, N228-3463.5 m, N228-3469.1 m, N228-3470.1 m, and N228-3479.4 m ranges from 1.4 to 3.7 wt% NaCl equivalent, 2.9 to 7.4 wt% NaCl equivalent, 1.4 to 7.3 wt% NaCl equivalent, 3.1 to 18.9 wt% NaCl equivalent, and 5.6 to 21.5 wt% NaCl equivalent, respectively. A salinity greater than 5 wt% NaCl equivalent is considered a highly saline fluid. The salinity of some aqueous inclusions measured in the experiment is high, which indicates that the shale is a good sealant. In contrast, the low salinity may be affected by the external fluid, that is, atmospheric fresh water.

Density Measurement of Methane Inclusions by Low Temperature Method
Methane inclusion is an isovolumetric system and its thermodynamic change occurs under the isovolumetric condition; hence, the phase transition point (Th) of methane inclusion can be used to determine its density (Liu and Shen, 1999;Liu, 2005).
In the process of micro-temperature measurement of methane inclusions, according to the micro-temperature measurement method of methane-rich inclusions proposed by previous studies (Kerkhof, 1990;Kerkhof and Thiéry, 2001), the slightly larger and easily observable methane inclusion was selected under the microscope and rapidly frozen using liquid nitrogen. The lowest temperature that can be achieved in the test is −196°C, and the three-phase point temperature of methane is −182.5°C. Theoretically, the appearance of the methane solid phase should be observed during the measurement of the pure methane inclusions, but this phenomenon is difficult to observe owing to the metastable state of the system (Liu, 2005).
Therefore, during the micro-temperature measurement of methane inclusions, the low-temperature phase transition process of methane inclusions was observed. Liquid nitrogen was used to rapidly freeze the methane inclusion, and when frozen to approximately −100°C, a small bubble pops out from the single-phase methane inclusion. There is no other phenomenon in the methane inclusion except that the bubble volume becomes larger when cooling is continued to −190°C, which is still a gas-liquid two-phase methane inclusion. The bubbles gradually decrease until they disappear and finally homogenize into the liquid phase when the temperature is slowly increased (Figure 8).
Therefore, the Th measurement results of the methane inclusions in this study varied from −93.9 to −92.1°C (Figure 8). From the low-temperature phase transformation process of methane inclusions, most methane inclusions in calcite veins in the Changning area appear to be pure methane inclusions, and the Th of this kind of inclusion is close to that of a pure methane system (−82.6°C) (Andersen and Burke, 1996).
The density of methane inclusions can be determined by methane inclusions Th (Liu and Shen, 1999 where ρ is the density of methane inclusions, g/cm 3 ; Th is the homogenization temperature of methane inclusions,°C. The results calculated by Formula 3 show that the methane inclusions density in Longmaxi shales fracture calcite veins in well N228 varied from 0.275 to 0.285 g/cm 3 , with an average of 0.281 g/cm 3 ( Figure 9; Table 2).

Density Measurement of Methane Inclusions by Raman Spectroscopy
As discussed above, the reliable acquisition of the traditional calculation of methane inclusion density data requires the selection of methane inclusions with regular morphology and long axis lengths greater than 6 μm for Th testing. Moreover, there are too many test steps that take a long time, and do not apply to all methane inclusions. Therefore, previous studies used the Raman scattering peak v 1 of methane inclusions to calculate the density (Seitz et al., 1996;Lin et al., 2007;Lu et al., 2007;Hansen and Berg, 2009). Lu et al. (2007) fitted a good linear relationship between the v 1 shift of methane Raman scattering peak and methane density, which is suitable for the density calculation of methane inclusions with methane content of 90%-100%, and the correlation coefficient is 0.9987: where ρ is the density of methane inclusions, g/cm 3 ; D = v 1 − v 0 , v 1 is the measured methane Raman scattering peak of methane inclusions, and v 0 is the methane Raman scattering peak of methane inclusions when the pressure is close to 0. The value of v 0 varies in different laboratories because of different calibration methods. In this study, the value of v 0 is calibrated by the Raman Laboratory in the Key Laboratory of the School of Geosciences, Yangtze University. Based on Formula 3, the density of methane inclusions varies greatly and can be divided into low-pressure gas-phase methane inclusions and supercritical high-density methane inclusions, taking the density of 0.162 g/cm 3 as the boundary. The density of lowpressure gas-phase methane inclusions ranges from 0.050 to 0.162 g/cm 3 , and that of supercritical high-density methane inclusions ranges from 0.210 to 0.293 g/cm 3 , with an average of 0.245 g/cm 3 ; whereas the density of methane-bearing bituminous inclusions ranges from 0.013 to 0.201 g/cm 3 ( Figure 9; Table 2). The density of methane inclusions trapped in sample N228-3479.4 m calculated by this method is similar to that obtained by the micro-temperature measurement of methane inclusions.

Trapped Pressure of Methane Inclusions
The trapped temperature and pressure of fluid inclusions record the paleo-temperature and paleo-pressure of the formation (Chen et al., 2002). The calculation of the trapped pressure is based on the equation of state applicable for the supercritical methane system established by Duan et al. (1992): Pc ; P is the pressure, bar; T is the temperature, K; R is the gas constant, 0.08314467 bar dm 3 K −1 mol −1 ; V is the molar volume, which can be determined by the density of methane inclusions and molar mass, dm 3 /mol; Z is the compression factor; P r and T r are the comparative pressure and temperature, respectively, and their dimensions are 1; P c and T c are critical pressure (46 bar) and critical temperature (190.4 K), respectively, and the unit is the same as P and T; a 1 = 0.0872553928; a 2 = −0.752599476; a 3 = 0.375419887; a 4 = 0.0107291342; a 5 = 0.0054962636; a 6 = −0.0184772802; a 7 = 0.000318993183; a 8 = 0.000211079375; a 9 = 0.0000201682801; a 10 = −0.0000165606189; a 11 = 0.000119614546; a 12 = −0.000108087289; α = 0.0448262295; β = 0.75397; γ = 0.077167.
The calculation of this formula requires the supercritical highdensity methane inclusions. Therefore, the trapped pressure of well N228 calculated using Formula 5 ranges from 51.44 to 129.53 MPa ( Table 2). The Th of FIAs trapped in well N209 studied by previous researchers is divided into two stages, ranging from 127.7 to 143.5°C and 159.2-171.2°C, with the corresponding pressures being    (Zhang et al., 2014;Gao et al., 2017), it is roughly calculated that the pressure coefficient at the maximum buried depth of Longmaxi shale in the Changning area before or at the initial stage of the Yanshanian uplift was between 1.51 and 1.88, and the shale was in a medium-strong overpressure state.

Characteristics of Tectonic Subsidence in the Changning Area
BasinMod basin simulation software was employed to restore the burial-thermal evolution history of Longmaxi shales in the deep depressional area of the Changning area, and to study the fluid evolution history of key drilling wells in combination with the results of fluid inclusions tests. For the burial history simulation of well N228, the data of drilling stratification and depth, the content of sandstone, mudstone and limestone of each layer, geothermal gradient, rock thermal conductivity, the measured TOC value of Wufeng-Longmaxi formation, kerogen type and hydrocarbon generation potential coefficient were mainly based on the actual data of the oilfield. The paleo-denudation amounts of Caledonian, Hercynian, and Indosinian were set as 200, 150 and 200 m, respectively (Zhu et al., 2009;Yuan et al., 2013;Huang et al., 2019). The amount of uplift and denudation in the late Mesozoic and Cenozoic was based on the apatite fission-track results of previous studies (Zhu et al., 2010;Liang et al., 2021), and the collected vitrinite reflectance was taken as a constraint. The burial history of the Changning shale gas field represented by well N228 reveals that the sedimentary process of the Wufeng-Longmaxi formation in the region has the sedimentary evolution characteristics of a typical multicycle craton basin Qiu et al., 2020) (Figure 10). After the initial deposition of the Longmaxi formation, it entered the rapid burial stage, with a burial depth of 1,800 m. Subsequent Garidonian tectonic movement in the middle-late Silurian led to the overall uplift and denudation of the strata, and the sedimentation and subsidence rate decreased rapidly. The strata re-entered the sedimentary subsidence stage after the Permian, while the Dongwu movement in the middle-Late Permian of Hercynian caused a brief uplift of the strata, and then it re-entered the sedimentary subsidence stage, which had the maximum stratum subsidence rate. Affected by the Indosinian movement in the middle-Late Triassic, the strata were briefly uplifted again. Later, the strata of the middle Triassic were denuded to different degrees at this stage. From the Late Triassic to Late Cretaceous, the Longmaxi formation continued to subside, with a maximum buried depth of about 6,600 m. At the end of the Late Cretaceous, the strata deposited from the Middle Jurassic to the Late Cretaceous were completely denuded by the influence of Himalayan tectonics. Some strata of the early Jurassic are currently exposed, and the bottom of the Longmaxi formation of well N228 is buried at a depth of 3480 m. The formation and evolution of the Emeishan mantle plume in the Late Permian led to a high heat flow (i.e., 60-80 mW m −2 ) in local areas of the Sichuan Basin. However, in general, the basin in the Paleozoic Cenozoic had a relatively stable geothermal flow (i.e., 60-65 mW m −2 ) and geothermal field characteristics (Zhu et al., 2010). Therefore, the simulation of the hydrocarbon generation history of the Longmaxi shales in the Changning area reveals that it began to enter a low maturity stage with Ro ranging from 0.5% to 0.7% at the end of the Early Silurian at about 420 Ma. However, controlled by Caledonian tectonic movement, the thermal evolution rate in the Early Paleozoic was slow. In the early Permian, at about 290 Ma, the Longmaxi shales entered a medium maturity stage with Ro ranges from 0.7% to 1.3%, during which a large amount of oil was generated from the source rock. This was proven by a large amount of bitumen observed in the veins and symbiotic aqueous inclusions, with bituminous inclusions measured in the veins (Figure 4). This oil generation stage lasted until the Middle Triassic for about 220 Ma. Subsequently, the organic matter entered a high maturity stage with Ro ranges from 1.3% to 2.0%, and the former oil was cracked to generate natural gas. Until the deep burial period in the early Cretaceous for about 135 Ma, the organic matter entered an over maturity stage with a Ro over 2.0%, and many chemical bonds of heavy hydrocarbons being broken to generate light hydrocarbons. Methane existed as a separate immiscible in the formation owing to supersaturation, a process evidenced by the observed methane inclusions. After the end of the Late Cretaceous, controlled by uplift and denudation, the earlier generated fracture veins opened, and the shale gas entered a stage of dynamic adjustment. The tectonic compression during the Himalayan period formed a series of tectonic fractures, making the migration, adjustment, and dissipation of shale gas in this area more intense. A series of low-temperature aqueous inclusions symbiotic with methane inclusions recorded the fluid activity information in this period ( Figure 4; Table 1).

Multi-phases Fluid Activity Process
Horizontal and high angle calcite veins are common in the Longmaxi black shales, and aqueous inclusions symbiotic with hydrocarbon inclusions have existed in the veins, all of which trapped the original fluids during the hydrocarbon generation and enrichment of the shale. This reflects the paleotemperaturepressure and salinity characteristics of the formation in different tectonic periods. The Th of FIAs trapped in well N228 has five peak characteristics of 95-115°C, 120-135°C, 140-155°C, 165-180°C, and 185-195°C, respectively. Combined with burial-thermal history, the trapped depth was defined and the time axis projection was made. The projection results were shown in Figure 11. The green arrow represents the trapped time of bituminous inclusion, that is, the trapped time of oil, and the red arrow represents the trapped time of the methane inclusion.
The gas-liquid two-phase aqueous inclusions symbiotic with bituminous inclusions of the first phase were enriched in the high angle fractured calcite veins, and the Th ranges from 95 to 115°C, indicating that they were trapped in the low-middle maturity stage of early shallow burial. The corresponding salinity has the characteristics of low salinity (1.4-3.9 wt% NaCl equivalent) and high salinity (5.1-18.9 wt% NaCl equivalent) (Figure 12), indicating that it may be affected by the contamination of high and low salinity fluids, that is, it can be significantly different from high Th-high salinity fluids and medium Th-low salinity fluids due to the mixing of freshwater fluids (Gao et al., 2017;Nie et al., 2020). This phase of fluid filling occurred during 280-292 Ma ( Figure 11). Therefore, it is ascribed to the Caledonian tectonic deformation and uplift denudation, resulting in the exposure of some faults to the surface, which is partially connected with shallow atmospheric freshwater.
The gas-liquid two-phase aqueous inclusions symbiotic with high-density methane inclusions of the second phase were enriched in horizontal fractured calcite veins, and the Th ranges from 185 to 195°C, the corresponding pressure ranges from 103.41 to 129.53 MPa, with the pressure coefficient ranges from 1.67 to 2.09, illustrating that they were trapped in the over maturity stage of deep burial. They were then trapped for about 103 Ma before the uplift, with the corresponding salinity being FIGURE 11 | Trapping time of FI of wells N228, N216, N209, and N215 in different synclines of the Changning area, southern Sichuan Basin. The data of wells N209, N215, and N216 collected by Gong (2019) and Liang et al. (2021).
FIGURE 12 | Scatter diagram of Th-salinity of FI trapped in wells N228, N216, N209, and N215, and data of wells N209, N215, and N216 collected by Gong (2019) and Liang et al. (2021). characterized by high salinity of 5.7 wt% NaCl equivalent. This demonstrates a better closure of the Longmaxi shales during the deep burial and a relatively closed fluid system. The formation was in a medium-strong overpressure state, that is, the starting point of pressure evolution was a medium-strong overpressure state.
The gas-liquid two-phase aqueous inclusions symbiotic with high-density methane inclusions of the third phase were enriched in the same horizontal fractured calcite veins, and the Th ranges from 165 to 180°C, and the corresponding pressure ranges from 88.67 to 109.22 MPa, with the pressure coefficient ranges from 1.56 to 1.92, manifesting that they were trapped in the uplift stage after deep burial. They were trapped for about 28 to 32 Ma, with the corresponding salinity being high at 11.1-21.5 wt% NaCl equivalent. This shows that in the uplift stage after deep burial, although the gas escape adjustment is caused by the opening fractures, it was still in a relatively closed fluid system and the formation is still in a medium-strong overpressure state.
The gas-liquid two-phase aqueous inclusions symbiotic with high-density methane inclusions of the fourth and fifth phases were enriched in the high angle fractured calcite veins, and the Th ranges from 140 to 155°C and 120 to 135°C. The corresponding pressure ranges from 57.47 to 74.50 MPa and 51.44 to 59.41 MPa, with the pressure coefficient ranging from 1.09 to 1.41 and 1.18 to 1.37, respectively. They were trapped for about 14 to 19 Ma and 6 to 10 Ma at the end of the uplift stage. The salinity corresponding to these two phases was low, indicating that in the later stage of uplift, owing to the openness of the fluid system and the uncertainty of preservation conditions, the underground fluid was affected by atmospheric freshwater to a certain extent, and the formation was in the state of the atmospheric-weak overpressure.
The trapping time and Th-salinity of wells N209, N215, and N216 in the shallow depressional area studied by previous researchers are projected in Figures 11, 12. As can be identified from Figures 11, 12, like well N228, these three wells have the same characteristics of low Th with high and low salinity in the early oil generation stage. It shows that during the Caledonian and Indosinian periods, due to the strata uplift, the fault may be connected to the surface, so it was affected by the contamination of atmospheric freshwater to a certain extent, while in the gas generation period after deep burial, the strata had good closure. Unlike well N228 deep in the depressional area, the wells in the shallow depressional area only experienced fluid filling and migration in the earliest stage of uplift after deep burial, with salinity ranging from 4.51 wt% NaCl equivalent to 5.41 wt% NaCl equivalent with a pressure coefficient ranging from 1.48 to 2.01, indicating that the formation closure was still good and the formation was in a medium-strong overpressure state in the early stage of Yanshanian uplift.

Impact of Fluid Activity on Resistivity
In the Changning area of the southern Sichuan Basin, it is found that the gas content is low or almost non-gaseous in the deeper buried Longmaxi shales of Luochang syncline, and high in the shallower buried Longmaxi shales of Jianwu syncline. That is, the gas content of well N228 located in the Luochang syncline is only 0.17 m 3 /t, whereas the gas content of wells N209, N215, and N216 located in the Jianwu syncline is high, being 2.99, 2.05, and 2.03 m 3 /t, respectively. In addition, low resistivity is found in wells with low or no gas content, compared to that of wells with high gas content is higher. The resistivity of wells N228, N209, N215, and N216 is 0.49 Ωm, 43.5 Ωm, 144.75 Ωm, and 61.05 Ωm, respectively ( Figure 13). In view of this phenomenon, two main controlling factors were proposed. One is that the low resistance is caused by the high conductivity mineral graphite produced by high evolution (Wang et al., 2018;Zhao et al., 2018), and the other is that the low resistance is mainly caused by inlaid close contact between particles due to diagenesis, resulting in lower gas saturation, increased water saturation, and increased mineralization within the shale to enhance conductivity (Gao et al., 2016;Liang et al., 2021).
It is evident from this investigation that there is an abnormal phenomenon of Dh > Gh with Rmc Ro > 4% in some bituminous inclusions, based on the variation of Dh/Gh in the thermal evolution of bitumen and the equivalent vitrinite reflectance calculated by the thermal maturity formula, which confirms the existence of high thermal evolution, i.e., "graphitization" of organic matter in the Luochang syncline of the Changning area, southern Sichuan Basin. Meanwhile, bituminous inclusions were also found in the Jianwu syncline of the Changning area, but these kinds of inclusions did not develop organic matter "graphitization." In addition, as discussed previously, it is considered that well N228 located in Luochang syncline existed more phases of fluid activities in the late uplift period than wells N209, N215, and N216 located in Jianwu syncline. Moreover, the fluid activity in this period corresponds to lower salinity, indicating that the late uplift and denudation have a greater impact on the region. The mixing of external fluids may destroy the preservation of shale gas, and the migration and escape of natural gas, water, and other fluids in the shale reservoir to varying degrees. The pore pressure decreases owing to gas dissipation, and the pressure of the overlying formation compacts the shale skeleton and reduces the pores and microfractures of the shale, resulting in a significant decrease in gas saturation and an increase in water saturation. The porosity of wells N209, N215, N216, and N228 is 4.50%, 4.56%, 4.25%, and 3.44%, respectively, and the gas saturation of wells N209, N215, N216, and N228 is 60.67%, 85.08%, 81.31%, and 33.80%, respectively ( Figure 13). The frequent fluid activity in the late uplift stage, the increase of water saturation in the formation, and the local area experiencing strong water washing (Liang et al., 2021), all led to the failure of the fluid self-closure system and the destruction of preservation conditions; thus, eventually leading to enhanced electrical conductivity and low resistance characteristics. This indicates that the emergence of low resistance in localized deep depressional areas, such as the Luochang syncline of the Changning area, is not only related to the high evolution degree of organic matter, but is also affected by the later multi-phases fluid modification.
The early subsidence rate of the Changning area was stable, with two oil generation peaks of 200-250 Ma and 290 Ma and a gas generation peak of 160-185 Ma ( Figure 11). It reached the maximum burial depth in the Late Cretaceous, when the pressure coefficient could reach two, with superior preservation conditions, good fluid storage and high enrichment of shale gas. In contrast, the multi-phase uplift and denudation experienced after the Late Cretaceous resulted in the reduction of the temperature and pressure of the Longmaxi shales. Meanwhile, the degrees of pore-fracture formation in the early stage and the stress differences of rocks affected by uplift in different tectonic sites led to differences in the degree of pore and fracture formation in the late stage. Compared with other wells in the shallow depressional area, well N228 in the deep depressional area was affected by tectonic uplift to a greater extent, with low porosity and strong fluid activity, resulting in strong adjustment and modification of shale gas, low gas saturation, and no effective shale gas production.

CONCLUSION
The FIAs trapped in the calcite veins in the Longmaxi shales of the Changning area were observed and measured to provide compelling evidence of the influence of multi-phase fluid activities on resistivity. The conclusions are summarized as follows: (1) The calcite vein orientation of the Longmaxi shales in the Changning area is considered mainly in a highly angular and horizontal fashion, with a huge quantity of bituminous inclusions and methane inclusions captured. The thermal evolution degree of bituminous inclusions reached the stage of carbon bitumen to metamorphic bitumen, the organic matter was fully cracked and gas was generated, and some of them exhibited the "graphitization" phenomenon. (3) The low resistance is not only caused by the "graphitization" of organic matter, but also related to the water washing of shale strata. The early burial situations in the Changning area were similar. The preservation conditions were superior, with better fluid storage and high shale gas enrichment. The Luochang syncline, i.e., the deep depressional area is more affected by tectonic uplift, with low porosity, strong fluid activities caused by fracture opening, and strong adjustment and transformation of shale gas, which makes it difficult to form an effective shale gas enrichment area compared with the Jianwu syncline, i.e., shallow depressional area.

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 authors. Frontiers in Earth Science | www.frontiersin.org June 2022 | Volume 10 | Article 931012