Reservoir Characteristics of the Lower Jurassic Lacustrine Shale in the Eastern Sichuan Basin and Its E ﬀ ect on Gas Properties: An Integrated Approach

: The exploration of shale gas in Fuling area achieved great success, but the reservoir characteristics and gas content of the lower Jurassic lacustrine in the northern Fuling areas remain unknown. We conducted organic geochemical analyses, Field Emission Scanning Electron Microscope (FE-SEM), X-ray di ﬀ raction (XRD) analysis, high-pressure mercury intrusion (MIP) and CH 4 adsorption experimental methods, as well as NMR logging, to study mineral composition, geochemical, pore structure characteristics of organic-rich shales and their e ﬀ ects on the methane adsorption capacity. The Da’anzhai shale member is generally a set of relatively thick (avg. 75 m) and high carbonate-content (avg. 56.89%) lacustrine sediments with moderate total organic carbon (TOC) (avg. 1.12%) and thermal maturation (Vitrinite reﬂectance (VR): avg. 1.19%). Five types of lithofacies can be classiﬁed: marl (ML), calcareous shale (CS), argillaceous shale (AS), muddy siltstone (MS), and silty shale (SS). CS has good reservoir quality with a high porosity (avg. 4.72%). The small pores with the transverse relaxation time of 0.6–1 ms and 1–3 ms comprised the major part of the porosity in the most lithofacies from Nuclear magnetic resonance (NMR) data, while the large pore ( > 300 ms) accounts for a small porosity proportion in the CS. The pores mainly constitute of mesopores (avg. 23.2 nm). The clay minerals with a large number of interparticle pores in the SEM contributes most to surface area in the shale lithofacies with a moderate TOC. The adsorption potential of shale samples is huge with an average adsorption capacity of 4.38 mL / g and also has high gas content (avg. 1.04 m 3 / t). The adsorption capacity of shale samples increases when TOC increases and temperature decreases. Considered reservoir properties and gas properties, CS with the laminated structures in the medium-upper section of the Da’anzhai member is the most advantage lithofacies for shale gas exploitation.


Introduction
A great success of shale gas exploration has been made in a Fuling gas field. Recently, many more oil and gas discoveries in lacustrine shale of the Lower Jurassic succession (e.g., Da'anzhai member shale) in the northern Fuling area gives a rise to the possibility of nonmarine shale gas accumulation in neighboring areas [1,2]. Recent studies focus on the pore structure and geochemical characteristics of the Da'anzhai member lacustrine shale [3][4][5]. There is little research concerning the influence of lacustrine shale reservoir characteristics on the adsorption capacity and total gas content. However, gas content is not only a key factor of favorable shale gas area evaluation, but also controls whether shale reservoir has a commercial exploration value or not [2,6].
In this study, we present a comprehensive approach of mineral composition and geochemical tests, high-resolution imaging investigation, high-pressure mercury intrusion (MIP) and CH 4 adsorption experiments, as well as NMR well logging, to (1) determine the mineral composition and geochemical, pore structure and petrophysical characteristics of the different shale lithofacies in the Da'anzhai members; (2) study the influence of continental reservoir characteristics on shale gas properties; (3) optimize the advantageous shale lithofacies to provide guidance for selecting the target shale gas layers.

Geological Setting
The study area is located in the eastern part of the Sichuan Basin, which is one of the major petroliferous basins in China. This area is dominated by the huge thrust-fold belts with an NNE or NE striking. It includes the narrow Datianchi, Huangnitang, Dachigan, Fangdoushan anticlines and gentle Liangping, Bashansi synclines between these structures from east to west [7]. The North Fuling gas field is located in the Bashansi Synclines ( Figure 1A,B). The Ziliujing formation in the Lower Jurassic is divided into Da'anzhai, Maanzhan, and Dongyuemiao members. The target layer of theDa'anzhai member consists of three lithological members: the shell limestone in the lower layer (J 1 zDa 1 ); the interbed in the middle layer with black shale, silty shale, and shell limestone (J 1 zDa 2 ); and the shell limestone and marl in the top (J 1 zDa 3 ) ( Figure 1C). The middle part of the Da'anzhai member (J 1 zDa 2 ) is the most important source rock with a thickness of 30-80 m. This layer was widely deposited when the lacustrine system was rapidly expanded and the water within the lake was deepened [8], which produced a positive impact on the accumulation and preservation of organic matter.
Energies 2020, 13, x FOR PEER REVIEW 2 of 17 accumulation in neighboring areas [1,2]. Recent studies focus on the pore structure and geochemical characteristics of the Da'anzhai member lacustrine shale [3][4][5]. There is little research concerning the influence of lacustrine shale reservoir characteristics on the adsorption capacity and total gas content. However, gas content is not only a key factor of favorable shale gas area evaluation, but also controls whether shale reservoir has a commercial exploration value or not [2,6].
In this study, we present a comprehensive approach of mineral composition and geochemical tests, high-resolution imaging investigation, high-pressure mercury intrusion (MIP) and CH4adsorption experiments, as well as NMR well logging, to (1) determine the mineral composition and geochemical, pore structure and petrophysical characteristics of the different shale lithofacies in the Da'anzhai members; (2) study the influence of continental reservoir characteristics on shale gas properties; (3) optimize the advantageous shale lithofacies to provide guidance for selecting the target shale gas layers.

Geological Setting
The study area is located in the eastern part of the Sichuan Basin, which is one of the major petroliferous basins in China. This area is dominated by the huge thrust-fold belts with an NNE or NE striking. It includes the narrow Datianchi, Huangnitang, Dachigan, Fangdoushan anticlines and gentle Liangping, Bashansi synclines between these structures from east to west [7]. The North Fuling gas field is located in the Bashansi Synclines ( Figure 1A,B). The Ziliujing formation in the Lower Jurassic is divided into Da'anzhai, Maanzhan, and Dongyuemiao members. The target layer of theDa'anzhai member consists of three lithological members: the shell limestone in the lower layer (J1zDa 1 ); the interbed in the middle layer with black shale, silty shale, and shell limestone (J1zDa 2 ); and the shell limestone and marl in the top (J1zDa 3 ) ( Figure 1C). The middle part of the Da'anzhai member (J1zDa 2 ) is the most important source rock with a thickness of 30-80 m. This layer was widely deposited when the lacustrine system was rapidly expanded and the water within the lake was deepened [8], which produced a positive impact on the accumulation and preservation of organic matter.

Samples and Experiments
In this paper, all the samples in the Da'anzhai member were selected from fresh cores in five wells, including FY1, FY4, FY5, FY3-2, and XL101 wells. Most of the samples were analyzed in the following methods: Total organic carbon determination, Rock-Eval pyrolysis analyses, kerogen vitrinite reflectance analysis, kerogen maceral composition microscopy, X-ray diffraction (XRD) testing, and high-pressure mercury intrusion (MIP), that were performed by the Marine Geology Lab, PetroChina Jianghan Oilfield Company. High-pressure intrusion experiments were conducted by AutoPore IV 9520 instruments, which can obtain a pore size range of 0.003-1000 µm. The pore micromorphology of 15 samples were determined at the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation of Chengdu University of Technology. These samples were investigated using an "FEI Quanta model 200F from FEI corporation, Netherlands" field-emission scanning electron microscopy (FE-SEM) under a working condition of 20 kv and distance of 8-9 mm. The methane adsorption isotherm experiments of five samples were conducted by the Shale Gas Lab, Sichuan Coal Geology Bureau.
Nuclear magnetic resonance (NMR) well logging is an approach that confirms the nuclear spin states of hydrogen in water and oil within a thin annulus situated deeply into the formation with several inches [9]. Under the pulsed and static magnetic field, the transverse relaxation time (T 2 ) of protons generally shows a positive relationship with the amount of hydrogen protons in the pore fluid [10]. The NMR well logging curves generally shows 2 to 3 main peaks and the peaks becomes larger with higher relaxation time: the first peak in the shale reservoir normally is the predominant one indicating clay-bound water and capillary-bound water, while free-fluid volume and micro-fractures are only revealed in some small peaks [11].
Field desorption is a relatively accurate method to directly obtain the data of gas content, including desorption gas content, losing gas content and residual gas content. Desorption gas content was measured by the ISO-300 system under a temperature of 80 • C ( Figure 2) when the sealed shale samples were placed in the desorption tank. Losing gas content was calculated by the United States Bureau of Mine (USBM) [12]. The samples (100-200 g) were placed in the residual gas tank, and then crushed for 30 min. Finally, residual gas content can be measured by the ISO-300 system ( Figure 2).

Samples and Experiments
In this paper, all the samples in the Da'anzhai member were selected from fresh cores in five wells, including FY1, FY4, FY5, FY3-2, and XL101 wells. Most of the samples were analyzed in the following methods: Total organic carbon determination, Rock-Eval pyrolysis analyses, kerogen vitrinite reflectance analysis, kerogen maceral composition microscopy, X-ray diffraction (XRD) testing, and high-pressure mercury intrusion (MIP), that were performed by the Marine Geology Lab, PetroChina Jianghan Oilfield Company. High-pressure intrusion experiments were conducted by AutoPore IV 9520 instruments, which can obtain a pore size range of 0.003-1000 μm. The pore micromorphology of 15 samples were determined at the State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation of Chengdu University of Technology. These samples were investigated using an "FEI Quanta model 200F from FEI corporation, Netherlands" field-emission scanning electron microscopy (FE-SEM) under a working condition of 20 kv and distance of 8-9 mm. The methane adsorption isotherm experiments of five samples were conducted by the Shale Gas Lab, Sichuan Coal Geology Bureau.
Nuclear magnetic resonance (NMR) well logging is an approach that confirms the nuclear spin states of hydrogen in water and oil within a thin annulus situated deeply into the formation with several inches [9]. Under the pulsed and static magnetic field, the transverse relaxation time (T2) of protons generally shows a positive relationship with the amount of hydrogen protons in the pore fluid [10]. The NMR well logging curves generally shows 2 to 3 main peaks and the peaks becomes larger with higher relaxation time: the first peak in the shale reservoir normally is the predominant one indicating clay-bound water and capillary-bound water, while free-fluid volume and micro-fractures are only revealed in some small peaks [11].
Field desorption is a relatively accurate method to directly obtain the data of gas content, including desorption gas content, losing gas content and residual gas content. Desorption gas content was measured by the ISO-300 system under a temperature of 80 °C ( Figure 2) when the sealed shale samples were placed in the desorption tank. Losing gas content was calculated by the United States Bureau of Mine (USBM) [12]. The samples (100-200 g) were placed in the residual gas tank, and then crushed for 30 min. Finally, residual gas content can be measured by the ISO-300 system ( Figure 2). The ISO-300 filed testing system of desorption gas content in the Da'anzhai shale member.1-thermostatic apparatus; 2-desorption tank; 3-dehumidication apparatus; 4-testing instruments of gas flow; 5-flexible pipe; 6-device for collecting gas; 7-pressure gauge; 8-thermometer; 9-data-processing system.

Organic Geochemical Characteristics
According to the TOC determination and Rock-Eval data of 143 samples in wells FY3-2 and XL101, The TOC of the Da'anzhai member widely ranged from 0.05 wt.% to 3.08 wt.% (avg. 1.12%) Energies 2020, 13, 4495 4 of 16 ( Figure 3 and Table 1). Some 65% of samples have a TOC value of more than (>0.5 wt.%). The TOC level over 1.0 wt.% and 2.0 wt.% account for 45.1% and 8.9% of all samples, respectively. The variation of TOC value in vertical is related to lithology ( Figure 3). The laminated calcareous shale and massive argillaceous shale in the upper part of the J 1 zDa 2 have high TOC value of avg. 1.58 wt.% [13], while the massive or bedded marl and silty mudstone in the J 1 zDa 1 and J 1 zDa 3 have a low TOC of <0.50 wt.%. The potential yield (P), defined as the sum of Rock-Eval S1 and S2 values, is generally used to evaluate the genetic potential of a source rock [14]. P (S1+S2) and TOC are the most important indicators of organic matter abundance [15], and they show a good positive correlation (R 2 = 0.82) for all samples ( Figure 3). P has a range of 1.03-6.84 mg/g (avg. 2.34 mg/g). The chloroform bitumen"A" content of 20 samples vary between 0.01% and 0.19% with an average value of 0.11% ( Figure 3). Based on the thresholds of nonmarine source rocks from [15], the Da'anzhai member shale is a fair-good source rock and provides a good quality for hydrocarbon generation.
The kerogen composition analysis of 42 samples indicate all samples are dominated by sapropelinite (65-74%), along with moderate vitrinite (25-51.7%) ( Figure 4). The sapropel organic matter is mostly correlated with initial production of algae or bacterial phytoplanktonic sources at the surface of lake [16]. The plot of HI vs. Tmax ( Figure 4) indicates organic matter (OM) in the Da'anzhai member shale is dominated by type II, mainly typeII 2 [17]. The OM type in vertical transitioned from mainly typeII 2 in J 1 zDa 3 to typeII 1 in the J 1 zDa 2 and J 1 zDa 1 . This suggests the change of the OM origins in vertical result in the diversity of kerogen compositions and OM types.
Energies 2020, 13, x FOR PEER REVIEW 4 of 17 ( Figure 3 and Table 1). Some 65% of samples have a TOC value of more than (>0.5 wt.%). The TOC level over 1.0 wt.% and 2.0 wt.% account for 45.1% and 8.9% of all samples, respectively. The variation of TOC value in vertical is related to lithology ( Figure 3). The laminated calcareous shale and massive argillaceous shale in the upper part of the J1zDa 2 have high TOC value of avg. 1.58 wt.% [13], while the massive or bedded marl and silty mudstone in the J1zDa 1 and J1zDa 3 have a low TOC of <0.50 wt.%. The potential yield (P), defined as the sum of Rock-Eval S1 and S2 values, is generally used to evaluate the genetic potential of a source rock [14]. P (S1+S2) and TOC are the most important indicators of organic matter abundance [15], and they show a good positive correlation (R 2 = 0.82) for all samples ( Figure 3). P has a range of 1.03-6.84 mg/g (avg. 2.34 mg/g). The chloroform bitumen"A" content of 20 samples vary between 0.01% and 0.19% with an average value of 0.11% ( Figure 3). Based on the thresholds of nonmarine source rocks from [15], the Da'anzhai member shale is a fair-good source rock and provides a good quality for hydrocarbon generation.   sapropelinite (65-74%), along with moderate vitrinite (25-51.7%) (Figure 4). The sapropel organic matter is mostly correlated with initial production of algae or bacterial phytoplanktonic sources at the surface of lake [16]. The plot of HI vs. Tmax ( Figure 4) indicates organic matter (OM) in the Da'anzhai member shale is dominated by type II, mainly typeII2 [17]. The OM type in vertical transitioned from mainly typeII2 in J1zDa 3 to typeII1 in the J1zDa 2 and J1zDa 1 . This suggests the change of the OM origins in vertical result in the diversity of kerogen compositions and OM types. The vitrinite reflectance (Ro) values of 17 samples ranged from 1.0% to 1.46% and averaged 1.19% (Figure 3). The temperature of maximum development of the S2 peak (Tmax) mainly ranges from 439 °C to 510 °C (avg. 452 °C). This indicates most of the samples reach a moderate-high maturity with a high potential for shale gas generation [18].

Mineral Composition and Reservoir Properties
X-ray diffraction analysis of core samples from the two wells showed that the Da'anzhai Member shale is dominated by carbonate minerals (avg. 56.8%), followed by clay minerals (avg. 29.8%), quartz and feldspars (avg. 19.6%), and pyrite (avg. 3.05%) ( Figure 3 and Table 1). These kinds of mineral compositions confirmed that the Da'anzhai member shale is not a typical shale of mainly terrigenous origin (generally defined by clay contents greater than 75%). The brittle mineral content (avg. 68.5%) of the Da'anzhai member shale has reached at the minimum standards of fracturing feasibility (40%) [19]. Illite and mixed-layerillite/smectite (avg. 65.8%) are the predominant clay minerals from the corresponding clay mineral data (Table 1) with comparatively small amounts of chlorite (avg. 15.8%) and kaolinite (avg. 18.2%). According to the mineral compositions, sedimentary genesis, and structure, the lithofacies of the Da'anzhai member can be classified as five groups: marl (ML), calcareous shale (CS), argillaceous shale (AS), muddy siltstone (MS), and silty shale (SS).
Porosity and permeability analyses from MIP in well FY3-2 indicates that the differences of reservoir properties between several lithofacies are dramatic. There is a weakly positive relationship between Porosity and permeability (R 2 = 0.31) ( Figure 5). The permeability of different lithofacies influenced by microfractures shows a wide range of 0.07-28.76 mD (Table 1). A permeability of 0.1-1 mD account for 66.7% in all samples. The permeability of CS with the laminated structures is relatively high. Moreover, its porosity (avg. 4.72%) is also higher than AS (avg. 2.72%), followed by The vitrinite reflectance (R o ) values of 17 samples ranged from 1.0% to 1.46% and averaged 1.19% (Figure 3). The temperature of maximum development of the S2 peak (T max ) mainly ranges from 439 • C to 510 • C (avg. 452 • C). This indicates most of the samples reach a moderate-high maturity with a high potential for shale gas generation [18].

Mineral Composition and Reservoir Properties
X-ray diffraction analysis of core samples from the two wells showed that the Da'anzhai Member shale is dominated by carbonate minerals (avg. 56.8%), followed by clay minerals (avg. 29.8%), quartz and feldspars (avg. 19.6%), and pyrite (avg. 3.05%) ( Figure 3 and Table 1). These kinds of mineral compositions confirmed that the Da'anzhai member shale is not a typical shale of mainly terrigenous origin (generally defined by clay contents greater than 75%). The brittle mineral content (avg. 68.5%) of the Da'anzhai member shale has reached at the minimum standards of fracturing feasibility (40%) [19]. Illite and mixed-layerillite/smectite (avg. 65.8%) are the predominant clay minerals from the corresponding clay mineral data (Table 1) with comparatively small amounts of chlorite (avg. 15.8%) and kaolinite (avg. 18.2%). According to the mineral compositions, sedimentary genesis, and structure, the lithofacies of the Da'anzhai member can be classified as five groups: marl (ML), calcareous shale (CS), argillaceous shale (AS), muddy siltstone (MS), and silty shale (SS).
Porosity and permeability analyses from MIP in well FY3-2 indicates that the differences of reservoir properties between several lithofacies are dramatic. There is a weakly positive relationship between Porosity and permeability (R 2 = 0.31) ( Figure 5). The permeability of different lithofacies influenced by microfractures shows a wide range of 0.07-28.76 mD (Table 1). A permeability of 0.1-1 mD account for 66.7% in all samples. The permeability of CS with the laminated structures is relatively high. Moreover, its porosity (avg. 4.72%) is also higher than AS (avg. 2.72%), followed by MS (avg. 1.89%) and SS (avg. 1.75), whereas porosity of ML is relatively poor (avg. 1.02%). The information of pore throat size, the sorting and connectivity within pores also can be obtained from the MIP data and curves ( Figure 6). The distribution of pore throat size indicates that AS, CS, and SS contains two types of pores: small pores (8-100 nm) and large pores (1-15 µm) with a good connectivity (Figure 7). But the tight ML and MS composed of a large number of small pores (<20 nm).
Energies 2020, 13, 4495 7 of 16 MS (avg. 1.89%) and SS (avg. 1.75), whereas porosity of ML is relatively poor (avg. 1.02%). The information of pore throat size, the sorting and connectivity within pores also can be obtained from the MIP data and curves ( Figure 6). The distribution of pore throat size indicates that AS, CS, and SS contains two types of pores: small pores (8-100 nm) and large pores (1-15 μm) with a good connectivity (Figure 7). But the tight ML and MS composed of a large number of small pores (<20 nm).   information of pore throat size, the sorting and connectivity within pores also can be obtained from the MIP data and curves ( Figure 6). The distribution of pore throat size indicates that AS, CS, and SS contains two types of pores: small pores (8-100 nm) and large pores (1-15 μm) with a good connectivity (Figure 7). But the tight ML and MS composed of a large number of small pores (<20 nm).   NMR is also an effective tool to determine the probe reservoir characteristics, including bound fluid, free fluid porosity, effective porosity, total porosity, and hydrocarbon existence [20]. The reservoir evaluation results of the Da'anzhai member from NMR well logging were shown in Figure 8 and Table 2. The heterogeneity of shale lithofacies strongly influenced the T 2 spectrum and porosity characteristics from NMR logging data. CS generally constitute of 2-3 peaks in the T 2 spectrum: The first peak (around 3 ms) is the largest peak, followed by a second peak (around 60 ms), while the third (1000 ms) peak is generally missing or it is not significant. This indicates that the pore types mainly consist of micropores and meso/macro pore. Sometimes, the fractures can be developed in these kinds of lithofacies. This is consistent with the distribution of pore throat size from the MIP results. The porosity of bound water is larger than that of free fluid. The T 2 spectrum of AS, SS, and MS shows an amplitude peak and a minor peak. However, the value of T 2 in the AS is generally higher than that of SS and MS. This indicates the AS has better reservoir quality. The T 2 spectrum of ML differs significantly from other lithofacies. It only has one peak around 1 ms, indicating that micropores prevail. Interestingly, Energies 2020, 13, 4495 8 of 16 free water volume does not exist in the ML because of high carbonate cementation. This is further confirmed by the carbonate content from XRD results. NMR is also an effective tool to determine the probe reservoir characteristics, including bound fluid, free fluid porosity, effective porosity, total porosity, and hydrocarbon existence [20]. The reservoir evaluation results of the Da'anzhai member from NMR well logging were shown in Figure  8 and Table 2. The heterogeneity of shale lithofacies strongly influenced the T2 spectrum and porosity characteristics from NMR logging data. CS generally constitute of 2-3 peaks in the T2 spectrum: The first peak (around 3 ms) is the largest peak, followed by a second peak (around 60 ms), while the third (1000 ms) peak is generally missing or it is not significant. This indicates that the pore types mainly consist of micropores and meso/macro pore. Sometimes, the fractures can be developed in these kinds of lithofacies. This is consistent with the distribution of pore throat size from the MIP results. The porosity of bound water is larger than that of free fluid. The T2 spectrum of AS, SS, and MS shows an amplitude peak and a minor peak. However, the value of T2 in the AS is generally higher than that of SS and MS. This indicates the AS has better reservoir quality. The T2 spectrum of ML differs significantly from other lithofacies. It only has one peak around 1 ms, indicating that micropores prevail. Interestingly, free water volume does not exist in the ML because of high carbonate cementation. This is further confirmed by the carbonate content from XRD results.    FE-SEM images are the techniques that are generally used to qualitatively determine the pore morphology characteristics of shales (e.g., the shapes, sizes, and distributions) [21][22][23]. The FE-SEM images show that a large quantity of micro-nano pores and microfracture exits may occur in the lacustrine shale samples. Based on the relationship between pores and minerals, pores can be classified into the three groups: organic matter-hosted pores (OM pores), the pores in the framework minerals (FM pores), and clay minerals-associated pores (CM pores). OM pores usually form during the thermal evolution and hydrocarbon generation process of kerogen and other OM. They mainly have aspherical, ellipsoidal-shape with a pore size of 2-100 nm ( Figure 9A,B). Compared with the marine shale in the Fuling areas [5], OM pores are scarce in the Da'anzhai Lacustrine shale because of a moderate maturity (R o : avg. 1.19%) and low TOC content. A large number of clay mineral pores was investigated in the SEM. The intrapartical and interlayer pores were developed between the clay mineral. The pore size of the observed wedge-shaped pore is 50 nm to 700 nm, and the length of the slits is up to 6 µm ( Figure 9C,E,F). These pores provide a large surface for adsorbed gas. FM pores are mostly a certain amount of nano-scale pyrite framboids or calcite intercrystalline pores ( Figure 9G,H). These pores are generally irregular with poor connectivity. The FM pores also contain the nano-scale pores generated by the dissolution of minerals (e.g., quartz and calcite) ( Figure 9I). Moreover, a large number of microcracks could be investigated in the CS and AS. These shrinkage cracks were generally caused by dehydration of clay minerals and thermal pressurization of hydrocarbons [24,25], forming a crack-pore network system ( Figure 9A,D). The shrinkage crack and the edge of the mineral particle form a bend-ridge-like distribution, with crack spacing of 20-30 nm and length of 8-40 µm. These microcracks provide important channels for gas migration.
Pore structure parameters measured by MIP are shown in Table 1. The average pore size and medium pore size show a main range of 17.8-57.9 nm and 2.7-18 nm, respectively. This indicates the pore size in the Da'anzhai member shale is dominated by mesopores, which is in consistent with the full pore-size distribution results of previous studies in Da'anzhai member shale of the central or northeastern Sichuan Basin [3][4][5]. The specific surface area varies from 1.42 to 13.92 m 2 /g (avg. 5.81 m 2 /g). The pore volume ranges from 0.35 to 7.01 cm 3 /100g (avg. 2.59 cm 3 /100g). Compared with the marine Longmaxi Formation shale in the same area (i.e., Fuling area; [5]), the surface area and pore volume is much lower because a lower TOC of lacustrine shale in Da'anzhai member contributes less to the surface and pore volume than marine shale. In Figure 10, TOC values are positively correlated with surface area and pore volume with a coefficient of 0.46 and 0.23, respectively. However, this correlation is much poorer because of a low-medium TOC value (avg. 1.12 wt.%) and moderate maturity. The clay mineral and surface area exit a strong positive correlation with a coefficient of 0.62. This indicates clay mineral play a more important role on the surface area [26] and pore volume of lacustrine shale with a low TOC (Figure 10). Many interparticle pores among clay minerals are also investigated in the SEM ( Figure 9C,E,F). On thecontrary, the carbonate content is negatively correlate with surface area and pore volume, because high carbonate minerals (e.g., calcite) with a strong cementation can greatly influence reservoir quality and usually have adverse impacts. The quartz minerals show a weakly positive relationship with pore volume (R 2 = 0.49). The total pore volume Energies 2020, 13, 4495 10 of 16 increases with quartz and feldspar content in the lacustrine shale due to the generation of pores by dissolution and weakening of compaction [27].
slits is up to 6 μm ( Figure 9C,E,F). These pores provide a large surface for adsorbed gas. FM pores are mostly a certain amount of nano-scale pyrite framboids or calcite intercrystalline pores ( Figure  9G,H). These pores are generally irregular with poor connectivity. The FM pores also contain the nano-scale pores generated by the dissolution of minerals (e.g., quartz and calcite) ( Figure 9I). Moreover, a large number of microcracks could be investigated in the CS and AS. These shrinkage cracks were generally caused by dehydration of clay minerals and thermal pressurization of hydrocarbons [24,25], forming a crack-pore network system ( Figure 9A,D). The shrinkage crack and the edge of the mineral particle form a bend-ridge-like distribution, with crack spacing of 20-30 nm and length of 8-40 μm. These microcracks provide important channels for gas migration. Pore structure parameters measured by MIP are shown in Table 1. The average pore size and medium pore size show a main range of 17.8-57.9 nm and 2.7-18 nm, respectively. This indicates the

Methane Sorption and Adsorption Potential
The methane adsorption isotherms of five shale samples measured at different pressure and temperatures are shown in Figure 11. We set a temperature group of 30 °C, 80 °C, 100 °C. The pressure values were measured at the 13 balance pressure points. The balancing time at each spot was set by 12.0 h. The amount of methane adsorbed into the shale at early stage of <5 MPa rapidly increased under different temperature and reached 85% of the maximum adsorption capacity.

Methane Sorption and Adsorption Potential
The methane adsorption isotherms of five shale samples measured at different pressure and temperatures are shown in Figure 11. We set a temperature group of 30 • C, 80 • C, 100 • C. The pressure values were measured at the 13 balance pressure points. The balancing time at each spot was set by 12.0 h. The amount of methane adsorbed into the shale at early stage of <5 MPa rapidly increased under different temperature and reached 85% of the maximum adsorption capacity. However, the shale with a high TOC content at same temperature (80 • C) has much faster methane adsorption rate with a large gradient (Figure 11A), because a high TOC of organic-rich shale provide a large amount of surface area for methane adsorption. The increased amount of methane adsorption became less from 5 MPa up to 12 MPa, because the methane adsorption reached equilibrium state. This process can be described by the Langmuir isotherm equation [28]: where V is the absorbed methane volume, cm 3 /g; V L is the Langmuir volume, cm 3 /g; b is the reciprocal of Langmuir pressure, MPa −1 ; p is the pressure, MPa. Langmuir volume of Da'anzhai member shale ranges from 1.34 mL/g to 2.28 mL/G (avg. 1.70 mL/g), indicating medium-high adsorbed capacity (Table 3). However, shale with different TOC content has quite different V L value at the same temperature (80 • C). When TOC content increased from 0.41 wt.% to 1.48 wt.%, the Langmuir volume of shale changed from 1.34 mL/g to 1.77 mL/g ( Figure 11). Even with a similar TOC level, the Langmuir volume decreased rapidly at a high temperature. When the experimental temperature jumped from 30 • C to 100 • C, the Langmuir volume of shale decreased from 2.28 mL/g to 1.44 mL/g. This can be interpreted by adsorption potential theory. However, the shale with a high TOC content at same temperature (80 °C) has much faster methane adsorption rate with a large gradient ( Figure 11A), because a high TOC of organic-rich shale provide a large amount of surface area for methane adsorption. The increased amount of methane adsorption became less from 5 MPa up to 12 MPa, because the methane adsorption reached equilibrium state. This process can be described by the Langmuir isotherm equation [28]: where V is the absorbed methane volume, cm 3 /g; VL is the Langmuir volume, cm 3 /g; b is the reciprocal of Langmuir pressure, MPa −1 ; p is the pressure, MPa. Langmuir volume of Da'anzhai member shale ranges from 1.34 mL/g to 2.28 mL/G (avg. 1.70 mL/g), indicating medium-high adsorbed capacity (Table 3). However, shale with different TOC content has quite different VL value at the same temperature (80 °C). When TOC content increased from 0.41 wt.% to 1.48 wt.%, the Langmuir volume of shale changed from 1.34 mL/g to 1.77 mL/g ( Figure 11). Even with a similar TOC level, the Langmuir volume decreased rapidly at a high temperature. When the experimental temperature jumped from 30 °C to 100 °C, the Langmuir volume of shale decreased from 2.28 mL/g to 1.44 mL/g. This can be interpreted by adsorption potential theory.   Based on Polanyi's theory of adsorption potential [29], the adsorption potential correlated with pressure can be defined as: where ε is the adsorption potential, J/mol; P i is the equilibrium pressure of methane at temperature T, MPa; P 0 is the saturated vapor pressure of methane at T, MPa; R is a constant value of 8.3114, J/(mol·K); T is the thermodynamic temperature, K; P is the pressure, MPa; under a certain temperature, the saturated vapor pressure of methane keeps constant. For example, the saturated vapor pressure of methane is generally 15.8 MPa at 80 • C. So according to Equation (2), the adsorption potential of shale samples can be easily determined from isothermal adsorption experiment ( Figure 11B). When the temperature increased under the same pressure, the adsorption potential obviously increased. It indicates the result of an increase of temperature and gas molecular kinetic energy. The increase of adsorption capacity cannot keep more gas molecules absorbed by the adsorbent, because the increasing rate of energy change is much lower than the kinetic energy. Therefore, the adsorption capacity increases with the decrease of temperature.

Gas Content of Field Desorption
Field desorption analysis of 18 samples indicates total gas content of the Da'anzhai member shale range from 0.14 to 2.29 m 3 /t (avg. 1.04 m 3 /t). It is much greater than the minimum standard of 0.5 m 3 /t for a favorable lacustrine shale gas [2]. Moreover, the CES in the medium-upper section of the J 1 zDa 2 has high gas content with a range of 1.15-2.29 m 3 /t (avg. 1.57 m 3 /t). Desorption gas content accounts for 10.56-60.44% (avg. 26.75%) of total gas content. This is much lower than the ratio of desorption gas content to total gas content due to a moderate TOC and maturity. TOC is positively related to the total gas content with a correlation coefficient of 0.71 ( Figure 11C). It shows a strongly positive relationship between TOC and desorption gas content (R 2 = 0.83). This suggests OM contributes most to desorption gas content. There is also a weakly positive correlation between porosity and total gas content ( Figure 11D), because shale with a high porosity can provide more sites for free gas.

Conclusions
(1) The lacustrine Da'anzhai member shale in the north Fuling area has a medium-high OM content (avg. TOC 1.12 wt.%). The laminated calcareous shale in the upper part of the J 1 zDa 2 has high TOC value of avg. 1.58 wt.%. The average potential yield P (S1 + S2) and chloroform bitumen content "A" are 2.34 mg hydrocarbons per g TOC and 0.11% with a medium hydrogen generation potential. The kerogen type of mainly typeII 2 and a moderate OM maturity (VR = avg. 1.19) indicate a huge potential of gas generation. The mineral content is dominated by calcite, clay and quartz. The clay minerals mainly consist of illite and mixed-layerillite/smectite. The lithofacies of the Da'anzhai member includes marl (ML), calcareous shale (CS), argillaceous shale (AS), muddy siltstone (MS), and silty shale (SS). (2) Heterogeneity of shale lithofacies strongly influences the reservoir quality. The porosity and permeability of shale is overall low, but the laminated calcareous shale has good reservoir properties with a high porosity (avg. 4.72%). The pore size is dominated by mesopores, and the medium pore diameter is 23.2 nm. Better NMR data show that the small pores with the transverse relaxation time of 0.6-1 ms and 1-3 ms comprised most of the porosity of the Da'anzhai shale member, while the larger pores with a T 2 value of 300-1000 ms and 1000-3000 ms only accounts for a small porosity proportion. However, it is very important for gas migration and generally showed a third peak in the NMR spectrum of CS, and this peak generally missed after centrifuge processing. (3) The clay minerals in Da'anzhai member shale contributes most to surface area. The average Langmuir volume and total gas content is 1.7 mL/g and 1.04 m 3 /t. The adsorption capacity of methane decreases with decreasing TOC and increasing temperature. The total gas content of lacustrine shale is controlled by TOC and porosity. CS have the highest gas content (avg. 1.57 m 3 /t) than other shale lithofacies. The reservoir properties, mineral composition and gas content data suggest the laminated calcareous shale in the medium-upper section of the J 1 zDa 2 are the most advantage lithofacies for shale gas potential production.