Pore Structure Characteristics and Permeability Stress Sensitivity of Jurassic Continental Shale of Dongyuemiao Member of Ziliujing Formation, Fuxing Area, Eastern Sichuan Basin

: A shale condensate gas reservoir with a high clay content and a high formation pressure was found in the Jurassic shale of the Dongyuemiao Member in the Fuxing area of the eastern Sichuan Basin. Reservoir characteristics and formation pressure have a signiﬁcant inﬂuence on optimal development. The present study investigated the continental shale of the Dongyuemiao Member in Well F. The petrological properties, physical properties, and pore structure of the Dongyuemiao Member were studied using X-ray diffraction (XRD), ﬁeld emission scanning electron microscopy (FE-SEM), N 2 adsorption, and mercury intrusion porosimetry (MIP). The permeability stress sensitivity characteristics of the shale reservoirs are discussed based on the change in shale porosity and permeability under overburden pressure. The tested shale samples yielded total organic carbon (TOC) and S 1 + S 2 values ranging mainly from 1.0 wt.% to 1.5 wt.% and from 0.39 to 2.28 mg/g, respectively, which was in the high maturity stage of the thermal evolution of organic matter (OM). The shales of the Dongyuemiao Member were found to contain high average clay mineral contents (more than 50%) of calcite and quartz, as well as albite, pyrite, dolomite, and halite. The main developments were identiﬁed as silica-rich argillaceous shale lithofacies, argillaceous shale lithofacies, and mixed argillaceous shale lithofacies. The pores were found to mainly be plate-like and ﬂake-like interlayer pores of clay minerals and OM pores with various shapes. The pore size was mainly concentrated below 110 nm, and the pore volume increment increased in ﬂakes with pore diameter. The average porosity and permeability of shale were found to be 4.827% and 0.243 mD, respectively. Clay minerals and quartz are beneﬁcial for improving the porosity and permeability of reservoirs, while carbonate minerals have the opposite effect. The permeability of the shale showed a negative exponential change with increasing effective stress under overburden pressure. When the effective conﬁning pressure was greater than 20 MPa, the decline rate of the shale permeability decreased with increases in the effective stress. The higher the clay mineral and TOC content, the stronger the stress sensitivity of shale permeability. The higher the carbonate mineral content, the weaker the stress sensitivity of shale permeability. The porosity sensitivity exponent indicates that matrix pores and micro-fractures are both developed in the Dongyuemiao Member, and the development of internal fractures is the main factor in the strong stress sensitivity of the shale permeability in the study area.


Introduction
With the successful revolution of marine shale hydrocarbon in the United States, the world energy landscape has transformed [1][2][3][4][5]. Micro-nano scale pores and fractures in shale affect oil hydrocarbon occurrence form and reservoir capacity [6,7]. Domestic and international scholars have analyzed pore structure from many aspects, including pore type, pore size and shape, pore volume and specific surface area, spatial distribution, and connectivity [8][9][10][11]. Large-scale organic-rich shale has developed in terrestrial strata such The Fuxing area is located in the Eastern Sichuan fold belt. It is geographically locat at the junction of Liangping, Dianjiang, Shizhu, Fengdu, and Zhongxian in Chongqi and northeast of Fuling. The main body of the study area is located in the Wanxian sy clinorium between the Datianchi and Fangdoushan anticline zones (Figure 1a). The Fu ing area in the eastern Sichuan Basin experienced multiple lake transgressions during t Jurassic. Three sets of organic-rich shales in Dongyuemiao, Da'anzhai, and Lianggaosh Member 2 were formed from bottom to top [31][32][33]. During the period of lake transgr sion in the Dongyuemiao section, sedimentary water bodies generally changed from sh low to deep to shallow in Fuxing. A set of littoral-shallow, lacustrine-semi-deep, and custrine gray-black shale facies with thin shell limestone was deposited. The limesto stratum was shell-rich, the clay mineral content of the shale was generally high, and t thickness of the sedimentary stratum was 60-70 m (Figure 1b). Combined with core o servations, logging responses, and whole-rock analysis, the Dongyuemiao Formation the study area can be further subdivided into three sub-members and seven sub-lay from bottom to top. The Dong1 sub-member is a gray-black shale with thin shell lim stone, including thin layers 1, 2, 3, and 4. The Dong2 sub-member is gray-black argil ceous limestone corresponding to thin layer 5; the Dong3 sub-member is dark gray a gray-black massive mudstone corresponding to thin layers 6 and 7 (Figure 1b) [34].

Samples
The research samples were derived from Dongyuemiao Member core samples c lected in Well F. A total of 38 samples were collected from 2752 and 2813 m (Figure and Table 1). The lithology was mainly gray-black shale with a thin shell-rich limeston Table 1. Stratigraphic organization and sampling information of the Dongyuemiao Member in Fuxing area of the Sichuan Basin.

Samples
The research samples were derived from Dongyuemiao Member core samples collected in Well F. A total of 38 samples were collected from 2752 and 2813 m ( Figure 1d and Table 1). The lithology was mainly gray-black shale with a thin shell-rich limestone.

Analytical Methods
XRD experiments were performed at the State Key Laboratory of Geological Process and Mineral Resources, China University of Geosciences (Wuhan, China) on an X-'Pert PRODY2198 X-ray diffractometer (PANalytical, Amsterdam, The Netherlands) with 40 kV of CuKα radiation and a 30 mA Ni filter.
The TOC content test sample was a 0.1-0.5 g powder sample with a particle size of 200 mesh. After the inorganic carbon content was removed with a 60 • C-80 • C acidic solution (1:7 hydrochloric acid (HCl): distilled water), the TOC was determined using a LECOCS-200 analyzer (LECO, St. Joseph, MI, USA). The experiment was conducted at the Key Laboratory of Tectonics and Petroleum Resources Ministry of Education, China University of Geosciences, Wuhan, China. A Rock-Eval instrument was used to conduct rock pyrolysis in order to determine the Tmax, HI, S1, and S2 values. The experiment was conducted at the Petroleum Exploration and Development Research Institute of Jianghan Oilfield Company (Wuhan, China).
A Zeiss-Merlin FE-SEM (Carl Zeiss AG, Battenwürburg Canton, Germany) was used for the experiment, with a resolution of 0.8 nm, a maximum magnification of 2 million times, and a test acceleration voltage of 0.02-30 kV. The samples were polished with argon prior to the experiment. The experiment was conducted at the Institute of Geology and Geophysics, Chinese Academy of Sciences.
The low-temperature N 2 adsorption-desorption test was performed using an ASAP2020 automatic specific surface area and microporous/mesoporous analyzer produced by the American Mac Company (Anton, Graz, Austria). Shale samples were powdered to 0.28-0.18 mm (60-80 mesh), and experiments were conducted at −195.8 • C, a relative pressure (P/P 0 ) range of 0.001-0.998, and a pore size range of 0.35-500.00 nm. The sample pore volume was calculated with the BJH method [35] by measuring the amount of nitrogen adsorption and desorption under equilibrium vapor pressure. The experiment was completed in the Key Laboratory of Tectonics and Petroleum Resources Ministry of Education, China University of Geosciences, Wuhan, China.
MIP was performed using an Auto-pore IV 9520 high-pressure mercury injection instrument produced by United States Instruments (Micromeritics, GA, USA), with a maximum pressure of 413 MPa and a pore throat test range of 800 µm-3 nm. All samples were cut into 1 cm 3 cubes, polished, and baked in an oven at 60 • C for 48 h to remove moisture and volatile substances. The experiment was conducted at the Key Laboratory of Tectonics and Petroleum Resources Ministry of Education, China University of Geosciences, Wuhan, China.
The overburden porosity and permeability experiment employed a PoroPDP-200 overburden porosity and permeability measuring instrument produced by Core Laboratories (Houston, TX, USA). The helium porosity measurement range was 0.01%-40%, the permeability measurement range was 0.00001-10 mD, and the overburden pressure was up to 70 MPa. Porosity measurements were based on Boyle's law, and the shale permeability was obtained with the pulse attenuation method. After the core was baked in an oven at 60 • C for 48 h, the porosity and permeability under different effective stress conditions were measured. The experiment was conducted at the Key Laboratory of Tectonics and Petroleum Resources Ministry of Education, China University of Geosciences, Wuhan, China.

Geochemical Characteristics
The Dongyuemiao shale in Eastern Sichuan has a high TOC content and hydrocarbon generation potential, slightly better organic matter type, and higher maturity, making it a better exploration horizon for shale oil than that in western Sichuan [30]. The experimental data showed that the TOC content of the Dongyuemiao shale is closely related to depth and stratigraphic horizon. It was found that the TOC content was mainly distributed between 1.0 wt.% and 1.5 wt.% (Figure 2a)

Geochemical Characteristics
The Dongyuemiao shale in Eastern Sichuan has a high TOC content and hydrocarb generation potential, slightly better organic matter type, and higher maturity, making better exploration horizon for shale oil than that in western Sichuan [30]. The exp mental data showed that the TOC content of the Dongyuemiao shale is closely related depth and stratigraphic horizon. It was found that the TOC content was mainly dist uted between 1.0 wt.% and 1.5 wt.% (Figure 2a)  The results of rock pyrolysis in the Dongyuemiao Member of the study area show that the shale hydrocarbon generation potential (S1 + S2) ranged from 0.39 to 2.28 mg and the hydrocarbon generation potential (S1 + S2) and TOC were positively correla ( Table 2). Most were high-quality source rocks, and some were good source rocks ( Fig  3). Previous studies have suggested that the vitrinite reflectance Ro of the shale in Dongyuemiao Member is between 1.52 and 1.58%, which generally indicates a high th mal maturity in the condensate oil wet gas-wet gas generation stage [34]. Most OM w of type II1-II2, predominantly layered and dispersed [33]. ble 2. Rock pyrolysis parameters of the Dongyuemiao shale in the Fuxing Area, Eastern Sichuan.

ormation
Depth (m) S1 (mg/g) S2 (mg/g) S1 + S2 (mg/g) TOC (wt%) Tmax (°C ) HI The results of rock pyrolysis in the Dongyuemiao Member of the study area showed that the shale hydrocarbon generation potential (S 1 + S 2 ) ranged from 0.39 to 2.28 mg/g, and the hydrocarbon generation potential (S 1 + S 2 ) and TOC were positively correlated ( Table 2). Most were high-quality source rocks, and some were good source rocks ( Figure 3). Previous studies have suggested that the vitrinite reflectance R o of the shale in the Dongyuemiao Member is between 1.52 and 1.58%, which generally indicates a high thermal maturity in the condensate oil wet gas-wet gas generation stage [34]. Most OM was of type II 1 -II 2 , predominantly layered and dispersed [33].  . Relationship between S1 + S2 and TOC of the Dongyuemiao shale in the Fuxing area, Eastern Sichuan. Area (a) represents poor hydrocarbon source rocks, area (b) represents medium-quality hydrocarbon source rocks, area (c) represents good hydrocarbon source rocks, and area (d) represents high-quality hydrocarbon source rocks.
With clay minerals, carbonate minerals, and felsic minerals as three-terminal ele- Sichuan. Area (a) represents poor hydrocarbon source rocks, area (b) represents medium-quality hydrocarbon source rocks, area (c) represents good hydrocarbon source rocks, and area (d) represents high-quality hydrocarbon source rocks.

Porosity and Permeability Distribution
According to the statistical analysis of the physical properties of shale reservoirs in Fuxing, the helium porosity of the Dongyuemiao Member is generally lower than 10%. Vertically, the porosity of thin layers 3 and 4 of Well F was found to be advantageous ( Figure 5). The porosity of shale samples was found to be concentrated at 2%-5%, accounting for 57.1% of the samples. The porosity of 23.8% of shale samples was 1%-2%, and the porosity of 14.3% of shale samples was 5%-10%. The overall permeability of the Dongyuemiao shale was found to be very low, with more than 95% of the shale samples having a permeability of less than 1 mD. In comparison, the permeability of thin layers 2-4 of the Dong1 sub-member was found to be higher; more than 80% of the shale samples showed a permeability of less than 0.5 mD, and the permeability concentration distribution range was 0.1-0.5 mD.  With clay minerals, carbonate minerals, and felsic minerals as three-terminal elements, the shale of the Dongyuemiao Member can be divided into eight lithofacies: silica-rich argillaceous shale, argillaceous shale, mixed argillaceous shale, mixed shale, clay-rich calcareous shale, siliceous/argillaceous mixed shale, clay-rich siliceous shale, and mixed calcareous shale. The upper part of the Dong1 sub-member mainly developed argillaceous shale lithofacies, while the lower part mainly developed mixed argillaceous shale lithofacies. The Dong2 sub-member mainly developed mixed shale lithofacies, and the Dong3 submember mainly developed silica-rich argillaceous shale lithofacies.

Porosity and Permeability Distribution
According to the statistical analysis of the physical properties of shale reservoirs in Fuxing, the helium porosity of the Dongyuemiao Member is generally lower than 10%. Vertically, the porosity of thin layers 3 and 4 of Well F was found to be advantageous ( Figure 5). The porosity of shale samples was found to be concentrated at 2%-5%, accounting for 57.1% of the samples. The porosity of 23.8% of shale samples was 1%-2%, and the porosity of 14.3% of shale samples was 5%-10%. The overall permeability of the Dongyuemiao shale was found to be very low, with more than 95% of the shale samples having a permeability of less than 1 mD. In comparison, the permeability of thin layers 2-4 of the Dong1 sub-member was found to be higher; more than 80% of the shale samples showed a permeability of less than 0.5 mD, and the permeability concentration distribution range was 0.1-0.5 mD.

Porosity and Permeability Distribution
According to the statistical analysis of the physical properties of shale reservoirs in Fuxing, the helium porosity of the Dongyuemiao Member is generally lower than 10% Vertically, the porosity of thin layers 3 and 4 of Well F was found to be advantageous ( Figure 5). The porosity of shale samples was found to be concentrated at 2%-5%, accounting for 57.1% of the samples. The porosity of 23.8% of shale samples was 1%-2%, and the porosity of 14.3% of shale samples was 5%-10%. The overall permeability of the Dongyuemiao shale was found to be very low, with more than 95% of the shale samples having a permeability of less than 1 mD. In comparison, the permeability of thin layers 2-4 of the Dong1 sub-member was found to be higher; more than 80% of the shale samples showed a permeability of less than 0.5 mD, and the permeability concentration distribution range was 0.1-0.5 mD.

Scanning Electron Microscopy
FE-SEM images showed that the main pore types of the shale in the Dongyuemiao Member include interlayer pores of clay minerals, intragranular pores of calcite, intergranular pores of pyrite, OM pores, and microfractures. Clay minerals were shown to be the most developed in the Dongyuemiao Member, and the interlayer pores of clay minerals were shown to be the most developed pore types. The pores were mostly plate-like and flake-like, with some connectivity (Figure 6a,b). Intragranular pores mainly developed inside the calcite particles as isolated pores (Figure 6c). Pyrite developed intergranular pores, and, due to the strong compressive resistance of rigid particles, intergranular pores could be seen in the periphery of pyrite. Some intergranular pores were found to be filled with OM, and OM pores developed in them. Pyrite protected the interlayer pores of clay minerals from being completely compacted to a certain extent, mostly as curved sheets (Figure 6d,e). The OM pores were irregular, slit shaped, and oval ( Figure 6f). The formation of microfractures in the Dongyuemiao shale is related to clay minerals and OM. Layered clay minerals are more likely to form microfractures with a greater degree of bending under the action of in situ stress (Figure 6g,h). Most of the micro-fractures related to OM developed at the contact between OM and other minerals, which is related to the hydrocarbon generation evolution of OM ( Figure 6i). ImageJ software was used to process the pore area and pore diameter curves corresponding to the FE-SEM images. The pore diameters were mainly found to be 0.01-10 µm, the pore area of the sample began to increase after 1 um, and the pore area in the picture with microfractures and clay mineral interlayer holes more obviously changed with the jagged change in pore size.

Scanning Electron Microscopy
FE-SEM images showed that the main pore types of the shale in the Dongyuemiao Member include interlayer pores of clay minerals, intragranular pores of calcite, intergranular pores of pyrite, OM pores, and microfractures. Clay minerals were shown to be the most developed in the Dongyuemiao Member, and the interlayer pores of clay minerals were shown to be the most developed pore types. The pores were mostly plate-like and flake-like, with some connectivity (Figure 6a,b). Intragranular pores mainly developed inside the calcite particles as isolated pores (Figure 6c). Pyrite developed intergranular pores, and, due to the strong compressive resistance of rigid particles, intergranular pores could be seen in the periphery of pyrite. Some intergranular pores were found to be filled with OM, and OM pores developed in them. Pyrite protected the interlayer pores of clay minerals from being completely compacted to a certain extent, mostly as curved sheets (Figure 6d,e). The OM pores were irregular, slit shaped, and oval ( Figure 6f). The formation of microfractures in the Dongyuemiao shale is related to clay minerals and OM. Layered clay minerals are more likely to form microfractures with a greater degree of bending under the action of in situ stress (Figure 6g,h). Most of the micro-fractures related to OM developed at the contact between OM and other minerals, which is related to the hydrocarbon generation evolution of OM ( Figure 6i). ImageJ software was used to process the pore area and pore diameter curves corresponding to the FE-SEM images. The pore diameters were mainly found to be 0.01-10 μm, the pore area of the sample began to increase after 1 um, and the pore area in the picture with microfractures and clay mineral interlayer holes more obviously changed with the jagged change in pore size.

Pore Structure Characterization
Samples were selected in each thin layer for N2 adsorption and MIP to reflect the pore structure of the Dongyuemiao shale (Table 3). According to the IUPAC classification [36], the adsorption isotherms of N 2 was of type IV, the maximum adsorption capacity was between 5 and 12 cm 3 /g, and the hysteresis loop showed H 4 type characteristics, indicating that parallel plate-like slit-type pores were developed in the shale. This type of pore is conducive to gas flow and easy to compress. The adsorption isotherm of the Jurassic Dongyuemiao shale hydrocarbon reservoir in the Fuxing area of Eastern Sichuan showed an anti-'S' shape. In the low-pressure stage (P/P 0 < 0.05), the gas adsorption capacity rapidly increased with increases in the relative pressure, which is the single molecule adsorption of micropores. In the medium-pressure stage (0.95 < P/P 0 < 0.05), adsorption mainly occurred through multi-molecular layer adsorption. In the final high-pressure stage, the adsorption capacity rapidly increased, and there was no adsorption saturation phenomenon, which is a characteristic of macropores and microfractures. In the N 2 desorption process, the relative pressure decreased from the maximum value. When 0.45 < P/P 0 < 1, the adsorption amount corresponding to the same pressure was higher than that in the adsorption process, forming a hysteresis loop that indicated the existence of nanoscale slit-like pores [37]. When P/P 0 < 0.45, the desorption and adsorption curves almost coincided (Figure 7a). Based on the BJH theoretical model, the pore size distribution curves of shale samples in different pore size ranges were obtained by processing the experimental data of N2 desorption. The curve shows that the pore size of the shale was mainly distributed within 10 nm (Figure 7b).
The MIP experiments could better reflect the characteristics of pore throat distribution and pore throat connectivity. Many nanoscale pores were developed in the Dongyuemiao shale, including some mesopores and nanoscale/micron macropores. When the mercury inlet pressure was less than 17,000 Psia (pore size range greater than 10 nm), the mercury inlet curve showed a large difference. The sudden increases in mercury saturation at 25 Psia and 100 Psia may have been caused by microfractures. When the mercury injection pressure was greater than 17,000 Psia (pore size range greater than 10 nm), more mercury flowed into the sample, and the different sample curves were not significantly different. The hysteresis loop showed that there was nearly 60% mercury retention, which may have been because some funnel-shaped pores in the sample limited the mercury removal process (Figure 7c). The pore-size distribution of the Dongyuemiao Member calculated from the MIP data had multi-peak characteristics, and the peaks were distributed from nanometers to micrometers. The maximum pore volume increment of shale in Well F was approximately 0.0035 cm 3 /g, the pore volume displaying 3-30 nm pore diameters was approximately 70%-80%, and the pore volume displaying 50-200 nm pore diameters was approximately 10%-15% ( Figure 7d). and the silica-rich argillaceous shale lithofacies F-1-2 showed a lower pore volume increment over the whole pore size range, indicating that the contribution of carbonate minerals to pore volume was not high compared with clay minerals. The macropore volumes in argillaceous shale lithofacies and mixed argillaceous shale lithofacies were larger, indicating that clay minerals were the main contributors to macropore volume ( Figure 8).   The combination of N2 adsorption and MIP experiments to characterize the pore size distribution of shale reservoirs has always attracted attention. In this study, considering the advantages of N 2 adsorption and MIP experiments in the characterization of mesopores and macropores, the boundary between mesopores and macropores (pore size of 50 nm) was selected for splicing. Splicing showed that the pore volume increased with increases in pore diameter. These jagged peaks indicated that the pore sizes in the shale were not evenly distributed, the volume increment of the mesoporous pores was high, and the pore volume increment showed a multi-modal state. The full pore size characterization showed that the Dongyuemiao shale had the most developed mesopores. The pore volume of the samples generally accounted for more than 50%, followed by the macro pore volume accounting for 18.4%-39.89% and the micropore volume accounting for the lowest fraction ( Table 3). The mixed shale lithofacies F-5-2 with a higher calcite content and the silica-rich argillaceous shale lithofacies F-1-2 showed a lower pore volume increment over the whole pore size range, indicating that the contribution of carbonate minerals to pore volume was not high compared with clay minerals. The macropore volumes in argillaceous shale lithofacies and mixed argillaceous shale lithofacies were larger, indicating that clay minerals were the main contributors to macropore volume (Figure 8).

Factors Influencing Reservoirs Properties
Shale reservoir pore types are diverse and the pore size ranges are large. Different pore types and sizes affect the occurrence of shale gas [38]. The physical and chemical properties of minerals and the development of OM ultimately determine the ability of shale to form pores. Here, clay minerals were found to maintain a good positive correlation with micropores, macropores, and the total shale pore volume (Figure 9), which was conducive to pore development. The correlation between quartz and micropores, mesopores, and total pore volume in the continental shale of the Dongyuemiao Member was not obvious, and it showed a significant negative correlation with macropores in the shale, which may be related to genesis and quartz content. The high quartz content (>35%) in biogenic marine siliceous shale is conducive to the development of mesopores and micropores [21]. The quartz content in the Dongyuemiao Member was generally shown to be lower than 35%, and the quartz minerals have mainly terrigenous inputs and support clay transformations [39]. Quartz is poorly correlated with the development of micropores and mesopores, and it is important in forming a rigid skeleton for the transformation of widely developed interlayer pore clay minerals, which is not conducive to the development of macropores. Calcite showed a good negative correlation with micropores, mesopores, macropores, and the total pore volume of the shale, indicating that calcite contributed little to the shale pores. The development of calcite according to FE-SEM was mainly in intragranular dissolved pores, and the overall surface porosity of the mineral surface was not high (Figure 6c).
Organic nanopores and clay mineral interlayer pores are the most widely developed in shale [6,40], and OM is accompanied by the formation of pores during hydrocarbon generation (Figure 6f,i). In particular, shale with a high TOC content and a high thermal evolution degree has a high secondary porosity during hydrocarbon generation and transformation [41]. The OM pores formed in shale are mostly micropores and mesopores [40,42,43]. The TOC content in the Dongyuemiao Member showed a positive correlation with the micropore volume. The positive correlation between organic carbon content and macropore volume may be related to the overpressured formation environment. The formation pressure coefficient obtained after the fracturing of Well F was approximately 1.93, and it is speculated that the original formation pressure coefficient was approximately 1.70, which indicates an ultra-high pressure system [34]. The thermal evolution of OM in the formation creates high pore fluid pressure, which is conducive to the interlayer pores of clay minerals not being compacted or closed. The mesoporous development of the Dongyuemiao continental shale has been affected by many factors, and the correlation with clay minerals, quartz, and OM is not obvious. The use of gas to measure the porosity and permeability of shale under macroscopic conditions can well-reflect real reservoir characteristics. The TOC content of the Dong1 and Dong2 sub-members was found to be high, which is of great significance to actual production. The effects of the mineral composition and TOC on the porosity and permeability of the Dong1 and Dong2 sub-members were experimentally studied ( Figure 10). The relationship among clay mineral content, porosity, and permeability was as follows: when the content was less than 20%, the porosity of most samples was less than 2% and the permeability was less than 0.2 mD. When the content of clay minerals was more than 20%, the porosity of most samples was more than 4% and the permeability of some samples was more than 0.2 mD, indicating that the interlayer pores developed by clay minerals were conducive to shale porosity development. However, owing to the influence of compaction, some interlayer pores in clay minerals were compressed, resulting in a decrease in permeability. Quartz is beneficial to the improvement of porosity and permeability: when its content was less than 15%, the porosity of most samples was less than 2% and the permeability was less than 0.2 mD. When the quartz mineral content was more than 15%, the porosity of most samples was more than 4% and the permeability was more than 0.2 mD. This shows that the presence of quartz particles is beneficial for the development of shale porosity, mainly because quartz particles can resist the overlying pressure and protect the pore space [44][45][46]. Overall, carbonate minerals are not conducive to reservoir porosity and permeability development. High-carbonate mineral samples were often found to correspond to thin-shell limestone interlayers developed in the Dongyuemiao Member, indicating that such interlayers have lower porosity and permeability, which is not conducive to hydrocarbon migration. Because it is affected by the thermal evolution of shale OM and its anti-compaction ability, the porosity of shale often first increases and then decreases with the TOC [47][48][49][50], and the TOC content of the Dongyuemiao Member was generally found to be low, which was slightly positively correlated with porosity and permeability. In contrast, the minerals corresponding to the shale with The use of gas to measure the porosity and permeability of shale under macroscopic conditions can well-reflect real reservoir characteristics. The TOC content of the Dong1 and Dong2 sub-members was found to be high, which is of great significance to actual production. The effects of the mineral composition and TOC on the porosity and permeability of the Dong1 and Dong2 sub-members were experimentally studied ( Figure 10). The relationship among clay mineral content, porosity, and permeability was as follows: when the content was less than 20%, the porosity of most samples was less than 2% and the permeability was less than 0.2 mD. When the content of clay minerals was more than 20%, the porosity of most samples was more than 4% and the permeability of some samples was more than 0.2 mD, indicating that the interlayer pores developed by clay minerals were conducive to shale porosity development. However, owing to the influence of compaction, some interlayer pores in clay minerals were compressed, resulting in a decrease in permeability. Quartz is beneficial to the improvement of porosity and permeability: when its content was less than 15%, the porosity of most samples was less than 2% and the permeability was less than 0.2 mD. When the quartz mineral content was more than 15%, the porosity of most samples was more than 4% and the permeability was more than 0.2 mD. This shows that the presence of quartz particles is beneficial for the development of shale porosity, mainly because quartz particles can resist the overlying pressure and protect the pore space [44][45][46]. Overall, carbonate minerals are not conducive to reservoir porosity and permeability development. High-carbonate mineral samples were often found to correspond to thin-shell limestone interlayers developed in the Dongyuemiao Member, indicating that such interlayers have lower porosity and permeability, which is not conducive to hydrocarbon migration. Because it is affected by the thermal evolution of shale OM and its anti-compaction ability, the porosity of shale often first increases and then decreases with the TOC [47][48][49][50], and the TOC content of the Dongyuemiao Member was generally found to be low, which was slightly positively correlated with porosity and permeability.
In contrast, the minerals corresponding to the shale with high porosity and permeability had a certain combination relationship: the content of clay minerals was 20%-70%, the content of quartz was 15%-30%, and the content of carbonate minerals was less than 40%.
Minerals 2022, 12, x 13 of 22 high porosity and permeability had a certain combination relationship: the content of clay minerals was 20%-70%, the content of quartz was 15%-30%, and the content of carbonate minerals was less than 40%.

Porosity-Permeability Characteristics of Continental Shale under Overburden Pressure
According to the functional relationship, the core stress sensitivity of a logarithmic model is the strongest, that of a power model is second, and that of an exponential or binomial model is the weakest [51]. In addition, the function model of stress sensitivity also reflects the pore geometry model of rocks, which reflects the stress sensitivity model of logarithmic and power relationships. Elliptical tube and round tube pores reflect the stress sensitivity model of a binomial relationship, and star tube pores reflect the stress sensitivity model of an exponential relationship [52]. The samples were collected from the shale core of the F well of the Dongyuemiao Member and were machined to approximately 25 mm in diameter and 57 mm in height. To better compare the stress sensitivity characteristics of the shale, a total of six samples were selected from the different lithologies of the Dong1, Dong2, and Dong3 sub-members to ensure that there was a large difference between the samples (Table 4). This experimental analysis was based on the Method for Determination of Porosity and Permeability of Rocks under Overburden Pressure (SY/T 6385-2016), and the effective stress was calculated using Formula (1).
where Peff is the effective stress, Pc is the confining pressure, Pp is the pore pressure, and nk is the permeability effective stress coefficient, which reflects the different sensitivities of rock permeability to external pressure and pore pressure. Bernabe [53] first proposed the concept of the effective stress coefficient of permeability when studying the effective stress of granite. Later, research on the effective stress coefficient of permeability showed that lithology is an important factor affecting the effective stress of permeability. Kwon et al. [54] obtained an effective stress coefficient of approximately 1 for illite-rich shale in the Wilcox Formation. The permeability effective stress coefficient of high-clay sandstone with a low stiffness ratio is around 1 [55]. Combining the characteristics of the sample and the existing research suggested that the effective stress coefficient was equal to 1.

Porosity-Permeability Characteristics of Continental Shale under Overburden Pressure
According to the functional relationship, the core stress sensitivity of a logarithmic model is the strongest, that of a power model is second, and that of an exponential or binomial model is the weakest [51]. In addition, the function model of stress sensitivity also reflects the pore geometry model of rocks, which reflects the stress sensitivity model of logarithmic and power relationships. Elliptical tube and round tube pores reflect the stress sensitivity model of a binomial relationship, and star tube pores reflect the stress sensitivity model of an exponential relationship [52]. The samples were collected from the shale core of the F well of the Dongyuemiao Member and were machined to approximately 25 mm in diameter and 57 mm in height. To better compare the stress sensitivity characteristics of the shale, a total of six samples were selected from the different lithologies of the Dong1, Dong2, and Dong3 sub-members to ensure that there was a large difference between the samples (Table 4). This experimental analysis was based on the Method for Determination of Porosity and Permeability of Rocks under Overburden Pressure (SY/T 6385-2016), and the effective stress was calculated using Formula (1).
where P eff is the effective stress, P c is the confining pressure, P p is the pore pressure, and n k is the permeability effective stress coefficient, which reflects the different sensitivities of rock permeability to external pressure and pore pressure. Bernabe [53] first proposed the concept of the effective stress coefficient of permeability when studying the effective stress of granite. Later, research on the effective stress coefficient of permeability showed that lithology is an important factor affecting the effective stress of permeability. Kwon et al. [54] obtained an effective stress coefficient of approximately 1 for illite-rich shale in the Wilcox Formation. The permeability effective stress coefficient of high-clay sandstone with a low stiffness ratio is around 1 [55]. Combining the characteristics of the sample and the existing research suggested that the effective stress coefficient was equal to 1. The burial depth of the shale in the Dongyuemiao Member of Well F was found to be between 2750 and 2815 m, the original formation pressure coefficient was approximately 1.70, the original formation pressure was approximately 48.48 MPa, and the maximum overburden pressure of this experimental design was 30 MPa. The porosity and permeability data of continental shale under overburden pressure and the variation characteristics with pressure are presented in Table 5 and Figure 11. There was a good correlation between the porosity and permeability of the sample and the effective stress. The correlation coefficient between porosity and effective stress greatly varied, ranging from 0.7746 to 0.9161, while the correlation coefficient between permeability and effective stress was above 0.92. The burial depth of the shale in the Dongyuemiao Member of Well F was found to be between 2750 and 2815 m, the original formation pressure coefficient was approximately 1.70, the original formation pressure was approximately 48.48 MPa, and the maximum overburden pressure of this experimental design was 30 MPa. The porosity and permeability data of continental shale under overburden pressure and the variation characteristics with pressure are presented in Table 5 and Figure 11. There was a good correlation between the porosity and permeability of the sample and the effective stress. The correlation coefficient between porosity and effective stress greatly varied, ranging from 0.7746 to 0.9161, while the correlation coefficient between permeability and effective stress was above 0.92.
With an increase in effective stress, the porosity and permeability of the sample exponentially decreased. The relationship between the shale porosity and effective stress satisfied the following equation: where φi is the porosity under the given stress condition (%), φ0 is the porosity when the initial stress is 0 (%), c is the sample pore compressibility (MPa −1 ), and p is the effective stress (MPa). According to the regression analysis equation between the porosity of the sample and the applied effective stress, the initial porosity φ0 and the compression coefficient c of the sample were obtained. The pore compressibility coefficient of the Dongyuemiao shale was 0.021-0.035 MPa −1 , with an average of 0.028 MPa −1 (Figure 11a). The relationship between the shale permeability and effective stress satisfied the following equation: where Ki is the permeability under a given stress condition (mD), K0 is the permeability when the initial stress is 0 (mD), γ is the permeability stress sensitivity coefficient (MPa −1 ), and p is the effective stress (MPa). The permeability stress sensitivity coefficient of Dongyuemiao shale was found to be 0.082-0.11 MPa −1 , with an average of 0.092 MPa −1 . In addition, when the effective stress exceeded 20 MPa, the porosity and permeability of the Dongyuemiao continental shale decreased with increasing pressure (Figure 11b). With an increase in effective stress, the porosity and permeability of the sample exponentially decreased. The relationship between the shale porosity and effective stress satisfied the following equation: where ϕ i is the porosity under the given stress condition (%), ϕ 0 is the porosity when the initial stress is 0 (%), c is the sample pore compressibility (MPa −1 ), and p is the effective stress (MPa). According to the regression analysis equation between the porosity of the sample and the applied effective stress, the initial porosity ϕ 0 and the compression coefficient c of the sample were obtained. The pore compressibility coefficient of the Dongyuemiao shale was 0.021-0.035 MPa −1 , with an average of 0.028 MPa −1 (Figure 11a). The relationship between the shale permeability and effective stress satisfied the following equation: where K i is the permeability under a given stress condition (mD), K 0 is the permeability when the initial stress is 0 (mD), γ is the permeability stress sensitivity coefficient (MPa −1 ), and p is the effective stress (MPa). The permeability stress sensitivity coefficient of Dongyuemiao shale was found to be 0.082-0.11 MPa −1 , with an average of 0.092 MPa −1 .
In addition, when the effective stress exceeded 20 MPa, the porosity and permeability of the Dongyuemiao continental shale decreased with increasing pressure (Figure 11b).

Stress Sensitivity Parameters of Permeability of Continental Shale Reservoir
In this study, we referenced the Method for Determination of Porosity and Permeability of Rocks under Overburden Pressure (SY/T 6385-2016) and the Method for Experimental Evaluation of Reservoir Sensitivity Flow (SY/T 5358-2010) for the stress sensitivity analysis of continental shale permeability. The stress sensitivity of shale was evaluated with the permeability damage rate and stress sensitivity coefficient ( Table 6). The permeability damage rate reflects the percentage of reservoir permeability damage under effective stress, defined as: where D k is the permeability stress damage rate corresponding to the effective stress (%), K 0 represents the permeability under initial effective stress (mD), and K i is the permeability under certain effective stress (mD). The permeability stress sensitivity coefficient reflects the degree of change of the formation permeability with the pressure. The larger the value of the permeability stress sensitivity coefficient, the more sensitive the permeability of the shale formation is to the pressure change, which is defined as: where γ is the permeability stress sensitivity coefficient (MPa −1 ), dK represents the change in permeability (mD), and dP represents the change in pressure (MPa).

Effect of Petrological Characteristics on Permeability Stress Sensitivity
Different mineral compositions have different mechanical properties. The mechanical properties of rocks determine the degree of deformation of the minerals that constitute the core under pressure, which affects the permeability stress sensitivity of shale. Wang et al. [56] observed that clay minerals are an important factor affecting stress sensitivity in a study of the stress sensitivity of tight sandstone reservoirs. Clay minerals not only make pore throats smaller but also make them more irregular in shape. Li et al. [57] found that the reason for the high permeability damage rate of mudstone samples with a high clay content is that stressing rocks with a high clay content makes the pore throats smaller. Furthermore, after being subject to stress, the clay easily spalls off and migrates to block pore throats. The shale in the Dongyuemiao Member is a typical continental shale with a high clay mineral content, and the pore types are mainly interlayer pores within clay minerals. Flat interlayer pores are beneficial for improving the porosity and permeability of the reservoir, but they are more easily compressed and deformed under compression. Accordingly, the permeability damage rate and permeability stress sensitivity coefficient were found to be positively correlated with the clay mineral content (Figure 12), indicating that the higher the clay mineral content, the stronger the permeability stress sensitivity of shale, which is similar to the response observed in marine shale [58]. There is high formation pressure in actual extraction, and there are more holes and fractures in shale after hydraulic fracturing. Combined with the characteristics of its condensate gas reservoir, to ensure the sustainability of shale gas exploitation, it is recommended to adopt pressure-holding techniques to reduce the overlying formation pressure borne by the gas during extraction. Pressure is applied to a shale reservoir so that many pores and fractures in the reservoir will be closed, resulting in a decrease in extraction efficiency. Quartz and other siliceous minerals are conducive to the formation of pores in marine shale [59], but the correlation between their content and permeability stress sensitivity parameters in continental shale is poor (Figure 12), indicating that the low content of quartz in continental shale has a limited contribution to pores and that pores are more derived from other factors. The porosity and permeability of limestone formations are generally low and without fracture development. The SEM images in this study show that carbonate minerals mainly developed isolated dissolution pores, which were less likely to be compressed and closed than the flat interlayer pores developed in clay minerals. Therefore, there was a negative correlation between carbonate mineral content and permeability stress sensitivity parameters ( Figure 12). The high-carbonate minerals in the Dongyuemiao Member improve the brittleness and permeability stress sensitivity of the reservoir and provide a strong fracturing and production capacity for the Dongyuemiao Member shale.

Effect of TOC Content on Permeability Stress Sensitivity
Under FE-SEM, these OM pores mostly appeared as round pores, stretched pores, and shrinkage cracks. The development of OM pores is controlled by many factors, such as the type and thermal maturity of the OM. The OM of the Dongyuemiao Member was found to mainly be of type II 1 and type II 2 , and the R O was found to be between 1.52% and 1.58%. Overall, it was shown to be in the high-maturity stage of OM thermal evolution and the formation stage of condensate oil-wet gas [34]. OM pores are formed after hydrocarbon discharge. Compared with other pore types, it is easier for OM pores to form interconnected pore networks and enhance reservoir porosity and permeability. Wang et al. [27] noted that the presence of OM in shale can lead to the significant deformation of matrix pores, resulting in a more pronounced change in permeability with pressure, and some OM pores could be seen as flat under the scanning electron microscope in this study (Figure 6f). A high OM content reduces the mechanical properties of shale, and OM pores are more easily destroyed during compaction [60]. Therefore, the permeability of shale with a high TOC content more obviously changes with effective stress. A good positive correlation between TOC content and permeability stress sensitivity parameters was also found in this study, indicating that the presence of OM in shale enhances the permeability stress sensitivity in this reservoir ( Figure 12).

Figure 12.
Relationship among TOC content, mineral composition, and permeability stress sensitivity parameters.

Effect of TOC Content on Permeability Stress Sensitivity
Under FE-SEM, these OM pores mostly appeared as round pores, stretched pores, and shrinkage cracks. The development of OM pores is controlled by many factors, such as the type and thermal maturity of the OM. The OM of the Dongyuemiao Member was found to mainly be of type II1 and type II2, and the RO was found to be between 1.52% and 1.58%. Overall, it was shown to be in the high-maturity stage of OM thermal evolution and the formation stage of condensate oil-wet gas [34]. OM pores are formed after hydrocarbon discharge. Compared with other pore types, it is easier for OM pores to form interconnected pore networks and enhance reservoir porosity and permeability. Wang et al. [27] noted that the presence of OM in shale can lead to the significant deformation of matrix pores, resulting in a more pronounced change in permeability with pressure, and some OM pores could be seen as flat under the scanning electron microscope in this study (Figure 6f). A high OM content reduces the mechanical properties of shale, and OM pores are more easily destroyed during compaction [60]. Therefore, the permeability of shale with a high TOC content more obviously changes with effective stress. A good positive correlation between TOC content and permeability stress sensitivity parameters was also found in this study, indicating that the presence of OM in shale enhances the permeability stress sensitivity in this reservoir ( Figure 12).

Effect of Pore Type on Permeability Stress Sensitivity
Fluid seepage in shale reservoirs is dependent on not only the matrix pores but also the fracture pores developed in shale. In an experiment on shale reservoirs in the Ordos Basin, Gao et al. [61] found that some artificial fractures were easier to close when the effective stress increased. When the effective pressure reached 24 MPa, the permeability loss rate of the fractured core was close to 100%, which was significantly larger than that of the matrix core. When the shale developed micrometer fractures, Li et al. [22] noted that micrometer fractures could not reopen after being closed under pressure, resulting in unrecoverable permeability. Ju et al. [26] observed that the permeability loss rate of unfilled natural microfractures is relatively high and that the solubility of fracture fillings has a Figure 12. Relationship among TOC content, mineral composition, and permeability stress sensitivity parameters.

Effect of Pore Type on Permeability Stress Sensitivity
Fluid seepage in shale reservoirs is dependent on not only the matrix pores but also the fracture pores developed in shale. In an experiment on shale reservoirs in the Ordos Basin, Gao et al. [61] found that some artificial fractures were easier to close when the effective stress increased. When the effective pressure reached 24 MPa, the permeability loss rate of the fractured core was close to 100%, which was significantly larger than that of the matrix core. When the shale developed micrometer fractures, Li et al. [22] noted that micrometer fractures could not reopen after being closed under pressure, resulting in unrecoverable permeability. Ju et al. [26] observed that the permeability loss rate of unfilled natural microfractures is relatively high and that the solubility of fracture fillings has a significant influence on the permeability stress sensitivity of fractured reservoirs. Previous studies have shown that pore type affects the stress sensitivity of rock permeability. To reflect the pore types developed in shale, Zhang et al. [62,63] simplified the various types of pores in shale cores. Using the circular pore, elliptical pore, fracture, and dual-porosity models, a method for judging pore types using the porosity sensitivity exponent was established. Using Equation (6), the ratio of permeability K i to initial permeability K 0 under certain pressure conditions is logarithmically dimensionless, the ratio of porosity ϕ i to initial porosity ϕ 0 is logarithmically dimensionless, and the porosity permeability power exponent α can be obtained.
The influence of the morphological characteristics of porous media on stress sensitivity can be reflected by α. In this study's experiment, the porosity sensitivity exponent was found to be 2.1567-4.5295 ( Figure 13). The pore structure developed in shale can be determined using the porosity sensitivity exponent. It can be seen that the type of shale pore development was mainly of the dual-porosity type, and microfractures and matrix pores were developed in the continental shale ( Figure 14). Implications from a pulse decay experiment on a fractured core revealed that a direct reduction in gas viscosity due to nanopore confinement may not have a straightforward impact on flow behavior in fractured shale reservoirs and that fracture permeability is much larger than matrix permeability [64]. Interlayer pores in clay minerals can be considered special microfractures ( Figure 6). The width and number of microfractures are important factors affecting shale permeability [65]. Here, the permeability damage rate of the experimental samples was much higher than the porosity damage rate. Most of the pores lost in this experiment were fracture pores or flat clay mineral interlayer pores. Theoretically, in a fracture model [62,63], the relationship between permeability, fracture porosity, and fracture width in the direction of fracture parallel to bedding is as follows: where b represents the average fracture aperture of shale formations (m), ϕ f represents the average porosity of shale formation fractures (%), and K f represents the average porosity of shale formation fractures (m 2 ).
φi to initial porosity φ0 is logarithmically dimensionless, and the porosity permeability power exponent α can be obtained.
The influence of the morphological characteristics of porous media on stress sensitivity can be reflected by α. In this study's experiment, the porosity sensitivity exponent was found to be 2.1567-4.5295 ( Figure 13). The pore structure developed in shale can be determined using the porosity sensitivity exponent. It can be seen that the type of shale pore development was mainly of the dual-porosity type, and microfractures and matrix pores were developed in the continental shale ( Figure 14). Implications from a pulse decay experiment on a fractured core revealed that a direct reduction in gas viscosity due to nanopore confinement may not have a straightforward impact on flow behavior in fractured shale reservoirs and that fracture permeability is much larger than matrix permeability [64]. Interlayer pores in clay minerals can be considered special microfractures (Figure 6). The width and number of microfractures are important factors affecting shale permeability [65]. Here, the permeability damage rate of the experimental samples was much higher than the porosity damage rate. Most of the pores lost in this experiment were fracture pores or flat clay mineral interlayer pores. Theoretically, in a fracture model [62,63], the relationship between permeability, fracture porosity, and fracture width in the direction of fracture parallel to bedding is as follows: where b represents the average fracture aperture of shale formations (m), φf represents the average porosity of shale formation fractures (%), and Kf represents the average porosity of shale formation fractures (m 2 ).  Equation (7) shows that the larger the fracture aperture, the larger the fracture porosity and the greater the permeability of the shale formation. Additionally, the influence of the fracture aperture was found to be greater than that of the fracture porosity on permeability. The matrix pores in shale were mainly distributed within 110 nm, and the microfracture aperture was shown to be the key to the stress sensitivity of the permeability. The permeability stress sensitivity coefficient of dual-porosity type II was found to be larger than that of dual-porosity type I. The change in permeability under the condition of overburden pressure also illustrates this point. When the effective stress was less than 5 MPa, the measured permeability exceeded the fitting curve. The permeability stress sensitivity in the low-pressure zone was stronger than in the high-pressure zone and showed the control of the microfractures in the shale on the permeability, whereas the permeability stress sensitivity in the high-pressure zone was lower than in the low-pressure zone ( Table 6). This showed that the influence of the matrix pores in the shale on permeability was more obvious. Therefore, at present, the fracturing technology of high-density fracture volume reconstruction in Well F in the Fuxing area [34] is suitable for the reconstruction of shale reservoirs in the Dongyuemiao Member. Subsequently, reasonable and efficient development should be carried out considering the gas reservoir characteristics and permeability stress sensitivity characteristics of the reservoirs.
than that of dual-porosity type I. The change in permeability under the condition of overburden pressure also illustrates this point. When the effective stress was less than 5 MPa, the measured permeability exceeded the fitting curve. The permeability stress sensitivity in the low-pressure zone was stronger than in the high-pressure zone and showed the control of the microfractures in the shale on the permeability, whereas the permeability stress sensitivity in the high-pressure zone was lower than in the low-pressure zone ( Table  6). This showed that the influence of the matrix pores in the shale on permeability was more obvious. Therefore, at present, the fracturing technology of high-density fracture volume reconstruction in Well F in the Fuxing area [34] is suitable for the reconstruction of shale reservoirs in the Dongyuemiao Member. Subsequently, reasonable and efficient development should be carried out considering the gas reservoir characteristics and permeability stress sensitivity characteristics of the reservoirs. Figure 14. Pore type discrimination diagram (The yellow area represents the matrix pore and the white part represents the crack) [26].

Conclusions
The TOC content, mineral composition, and physical properties of the Dongyuemiao Member showed obvious heterogeneity in the vertical direction. The TOC content was mainly distributed between 1.0 wt.% and 1.5 wt.%, and the TOC of thin layer 4 was the highest. The mineral composition was found to mainly comprise clay minerals, followed by quartz and calcite, and the content of clay minerals was the highest in the Dong1 submember (average = 72.79%). The calcite content was the highest in the Dong2 sub-member (average = 50.55%). The porosity was concentrated in 2~5%, and the permeability was concentrated in 0.1~0.5 mD. The porosity and permeability of thin layers 3 and 4 were highest.
The inorganic pores in the Dongyuemiao Member mainly include clay mineral interlayer pores, calcite intragranular pores, pyrite intercrystalline pores, and microfractures. Clay minerals were found to be the main contributors to pore volume. The pore diameter of the Dongyuemiao Member is mainly concentrated below 110 nm, and the mesopores are the most developed, followed by the macropore volume and the micropore volume. A large number of parallel plate-like clay minerals interlayer pores provide space for petroleum enrichment. Figure 14. Pore type discrimination diagram (The yellow area represents the matrix pore and the white part represents the crack) [26].

Conclusions
The TOC content, mineral composition, and physical properties of the Dongyuemiao Member showed obvious heterogeneity in the vertical direction. The TOC content was mainly distributed between 1.0 wt.% and 1.5 wt.%, and the TOC of thin layer 4 was the highest. The mineral composition was found to mainly comprise clay minerals, followed by quartz and calcite, and the content of clay minerals was the highest in the Dong1 sub-member (average = 72.79%). The calcite content was the highest in the Dong2 submember (average = 50.55%). The porosity was concentrated in 2~5%, and the permeability was concentrated in 0.1~0.5 mD. The porosity and permeability of thin layers 3 and 4 were highest.
The inorganic pores in the Dongyuemiao Member mainly include clay mineral interlayer pores, calcite intragranular pores, pyrite intercrystalline pores, and microfractures. Clay minerals were found to be the main contributors to pore volume. The pore diameter of the Dongyuemiao Member is mainly concentrated below 110 nm, and the mesopores are the most developed, followed by the macropore volume and the micropore volume. A large number of parallel plate-like clay minerals interlayer pores provide space for petroleum enrichment.
Clay minerals and OM can enhance the permeability stress sensitivity of continental shale, and carbonate minerals can reduce the permeability stress sensitivity of continental shale. The average content of brittle minerals such as quartz, calcite, feldspar, pyrite, and dolomite in the Dongyuemiao Member was found to be about 40%, which provides a certain fracturing and exploitation capacity for continental shale. Both the permeability variation and the porosity sensitivity exponent under overburden pressure indicate that microfractures and matrix pores are also developed in continental shale and have a high permeability stress sensitivity. Combined with the characteristics of the gas reservoir and the permeability stress sensitivity of the reservoir, it is recommended that Well F should be subjected to pressure-holding mining under the existing fracturing treatment featuring a "dense-fracture volume stimulation".