Types and Quantitative Characterization of Microfractures in the Continental Shale of the Da’anzhai Member of the Ziliujing Formation in Northeast Sichuan, China

: A number of wells in the Sichuan Basin of China have tested industrial gas ﬂow pressure arising from the shale of the Da’anzhai section of the Ziliujing Formation, revealing good exploration potential. Microfractures in shales affect the enrichment and preservation of shale gas and are important storage spaces and seepage channels for gas. In order to increase productivity and to reduce the risks associated with shale gas exploration, the types, connectivity, and proportion of microfractures in the Da’anzhai Member have been studied in this work by core and thin section observations, micro-CT, scanning electron microscopy, nitrogen adsorption, and high-pressure mercury intrusion. The results show that four types of fractures have developed in the shale of the Da’anzhai section: microfractures caused by tectonic stress, diagenetic shrinkage fractures of clay minerals, marginal shrinkage fractures of organic matter, and microfractures inside mineral particles. Among these, structural fractures and organic matter contraction fractures are the main types and are signiﬁcant for shale reservoirs and seepage. The structural microfractures are mainly opened and are well-developed in the shale, with a straight shape, mainly between bedding, with the fracture surface being curved, fully opened, and mainly tensile. Organic matter fractures often develop on the edge of the contact between organic matter and minerals, presenting a slit-like appearance. The fractures related to bedding in the shale are particularly developed, with larger openings, wider extensions, intersecting and expanding, and forming a three-dimensional interconnected pore-fracture system. Based on image recognition, generally speaking, microfractures account for about 20% of the total pore volume. However, the degree of the microfractures’ development varies greatly, depending upon the structural environment, with the proportion of microfractures in fault-wrinkle belts and high-steep zones reaching 40% to 90% of the total pore space. On the other hand, micro-fractures in areas with underdeveloped structures account for about 10% of the total pore space.


Introduction
Fractures in shale affect the enrichment and preservation of shale gas [1][2][3][4][5], where they form the seepage channels for gas. Fracture studies in the shale gas exploration stage can therefore provide a reference for the optimal areas for shale gas exploration [6][7][8][9]. In the development stage of shale gas reservoirs, the effectiveness and permeability anisotropy of fractures directly affects the pattern of the deployment of the extraction wells and the reservoir fracturing reformation [10][11][12][13], which is an important reference for horizontal well deployment and geosteering [14]. In addition, the existence of shale cracks provides storage space for shale gas [15][16][17]. Fractures can be classified according to their cause and scale, and there are various classification schemes [18]. Due to their various scales and development characteristics, these fracture types have different effects on shale gas accumulation and enrichment [19][20][21]. For example, shrinkage fractures and abnormally high-pressure fractures account for a small proportion and are mostly filled by later mineralization [22][23][24], so they have little impact on the enrichment and preservation of shale gas. However, interlayer shear fractures, bedding detachment fractures, interlayer expansion fractures, and interlayer fractures have a greater impact on gas enrichment and preservation [25][26][27][28].
Microfractures are mainly micro-scale cracks smaller than 100 µm [29], which are difficult to identify with the naked eye, and their evaluation requires observation and statistical analysis with the help of microscopy instrumentation. In association with macroscopic fractures, microfractures also have an important influence on the physical and gas-bearing properties of shales [30,31]. In the process of shale gas development, natural gas is usually produced by gradual release, and the free gas in microfracture pores is discharged first [32]. Then, the adsorbed gas attached to the microfracture surface is released [33,34]. Finally, adsorbed gas escapes from the shale matrix. The desorption process of adsorbed gas is very slow, but when microfractures develop near the bottom of a well, the free gas in the microfractures quickly discharges, forming a low-pressure zone near the microfractures, in turn promoting the desorption of adsorbed gas in the area near the fractures and improving the gas production efficiency [35,36]. Microfractures in the shale can be used as penetration channels. Pores in a shale matrix are extremely undeveloped, being mostly microcapillary pores, and the permeability is much lower than that of tight sandstone; however, microfractures can significantly increase the permeability of this material [37][38][39][40]. Therefore, the development of effective microfractures is of great significance for increasing the permeability of shale. Previous descriptions of fractures have mainly focused on the core to outcrop scale, with few studies describing micro-scale fractures or the influence of microfracture development on shale reservoir space, with no quantitative studies to the authors' knowledge of the influence of microfractures on the physical properties of shales.
In this work, we take the Yuanba 21, Yuanba 102, Yuanlu 4, and another 10 wells in Northeast Sichuan and carry out a detailed description, classification, and quantitative characterization of microfractures based on a large number of observations of core and thin section samples, micro-CT, scanning electron microscopy, nitrogen adsorption, and high-pressure mercury injection tests. These results potentially play an important role in clarifying the effective reservoir space and gas-bearing properties of continental shale gas in the Sichuan Basin and are of great significance in clarifying the enrichment mechanism of continental shale gas and in supporting the exploration and deployment of continental shale gas resources.

Geological Setting
The study area is located in the northeastern part of the Sichuan Basin, namely around the Longmen Mountain, Micang Mountain, and Daba Mountain front [41][42][43]. The structural division is located in the gentle fold belt of northern Sichuan and the highland fold belt of eastern Sichuan. North Sichuan is the main area for the exploration and development of shale gas reservoirs for the entire Sichuan Basin [44][45][46][47]. The study area is located on the northern margin of the Yangtze Plate and is adjacent to the Qinling fold belt. The Lower Jurassic Da'anzhai Member of the North Sichuan Basin is mainly distributed across the Micangshan and Dabashan arc structures [48][49][50]. The belt overlaps with the arcshaped skirt belt in eastern Sichuan. During the Jurassic period, due to the disappearance of the Songpan Ganzi Sea and the Western Sea, all of the surrounding areas were surrounded and evolved into inland lakes [51][52][53]. The Early to Middle Jurassic basins were deep lakes, semi-deep lakes, and shallow lakes, where rivers, lakes, and deltas formed [54,55]. The facies are mainly composed of interbedded sand and mudstone. Among them, there are dark clastic rocks and carbonate strata from the Early to early Middle Jurassic and continental red strata from the late Middle Jurassic to the Paleogene [56,57]. The Lower Jurassic Da'anzhai Member is composed of lakeside, shallow lake, and semi-deep lake deposits, with a thickness ranging between 100 to 160 m in northeastern Sichuan ( Figure 1). These deposits are mainly composed of shell limestone and dark shale interbedded with each other with unequal thicknesses, and siltstone can be seen locally [58]. These deposits are rich in freshwater bivalves and ostracod fossils, and are common in scale deposits. Among these deposits, the second member of the Da'anzhai Member is a semi-deep to deep lake deposit, and is the main mudstone development interval in this section [59,60]. disappearance of the Songpan Ganzi Sea and the Western Sea, all of the surrounding areas were surrounded and evolved into inland lakes [51][52][53]. The Early to Middle Jurassic basins were deep lakes, semi-deep lakes, and shallow lakes, where rivers, lakes, and deltas formed [54,55]. The facies are mainly composed of interbedded sand and mudstone. Among them, there are dark clastic rocks and carbonate strata from the Early to early Middle Jurassic and continental red strata from the late Middle Jurassic to the Paleogene [56,57]. The Lower Jurassic Da'anzhai Member is composed of lakeside, shallow lake, and semi-deep lake deposits, with a thickness ranging between 100 to 160 m in northeastern Sichuan ( Figure 1). These deposits are mainly composed of shell limestone and dark shale interbedded with each other with unequal thicknesses, and siltstone can be seen locally [58]. These deposits are rich in freshwater bivalves and ostracod fossils, and are common in scale deposits. Among these deposits, the second member of the Da'anzhai Member is a semi-deep to deep lake deposit, and is the main mudstone development interval in this section [59,60].

Sampling and Laboratory Methods
In this study, samples were collected from 15 wells in the Yuanba area of the Sichuan Basin (YB 21, YB102, YB 104, YB 122, YB 221, YB 273, YL 4, YL 30. YL 17, YL 171, YL 175, YL 176, FY 1, FY 4, FY 5), with cores from the Da'anzhai section of the Ziliujing Formation. Some of the well positions are shown in Figure 1. The clay minerals in the Da'anzhai lake facies mud shale are very developed, making up from 21% to 56%, with most being about 40% (illite and Imonite mixed layer account for 50% and 25% of the clay minerals, respectively). The second most common are brittle minerals such as quartz and feldspar with a content of 30% to 45%, and 5% to 25% of carbonate minerals. The organic content is about 0.8%, and the mineral composition of the shale constitutes the material basis for the development of microscopic pores and fractures. The organic matter content of the lacustrine mudstone in the Da'anzhai section of the Yuanba area varies significantly, with the TOC total organic carbon (TOC) values of the analyzed samples ranging from 0.04 to 2.30%, of which more than half of the samples were higher than 0.5%, and 27% of the samples were higher than 1.0%, with some belonging to type III kerogen, and some to type II kerogen. The Ro of the samples shows little variation, ranging from 1.29% to 1.97%, with an average of 1.63%. Readers should note that the analyses undertaken for this work were carried out at the State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum (Beijing), except where stated otherwise.

Sampling and Laboratory Methods
In this study, samples were collected from 15 Figure 1. The clay minerals in the Da'anzhai lake facies mud shale are very developed, making up from 21% to 56%, with most being about 40% (illite and Imonite mixed layer account for 50% and 25% of the clay minerals, respectively). The second most common are brittle minerals such as quartz and feldspar with a content of 30% to 45%, and 5% to 25% of carbonate minerals. The organic content is about 0.8%, and the mineral composition of the shale constitutes the material basis for the development of microscopic pores and fractures. The organic matter content of the lacustrine mudstone in the Da'anzhai section of the Yuanba area varies significantly, with the TOC total organic carbon (TOC) values of the analyzed samples ranging from 0.04 to 2.30%, of which more than half of the samples were higher than 0.5%, and 27% of the samples were higher than 1.0%, with some belonging to type III kerogen, and some to type II kerogen. The Ro of the samples shows little variation, ranging from 1.29% to 1.97%, with an average of 1.63%. Readers should note that the analyses undertaken for this work were carried out at the State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum (Beijing), except where stated otherwise.

Rock Slices and Identification
Rock sections from the samples were polished and identified. To prevent man-made cracks from occurring, the samples were cut into sections of 3 to 5 cm in length and width with a wire cutting device. Then, sample which was processed was then glued onto the ground glass with epoxy resin. After gluing, the sample was placed on a fixed platform and baked at a low temperature of 50 • C for about 6 to 8 h to consolidate and harden to avoid rock fracturing caused by later observation. After the glue hardened, the specimen was cut on a petro-thin machine and ground to a thickness of between 100 and 150 µm, then ground on a Lapping plate to a standard thickness of 30 µm. A Leica DM4P polarized light microscope was used to observe the prepared rock slices, to determine the mineral composition of the sample, to study its structure, to analyze the sequence of mineral generation, and to determine the rock type and its microcrack characteristics.

Micro-CT Scanning
The micro-CT experiment was carried out at the Petroleum Geology Experiment Center, PetroChina Research Institute of Petroleum Exploration and Development. The instrument was a XRadia Ultra-XRM L200 stereoscopic microscope, with the voltage set to 8 kV, and the instrument's limit resolution being 0.7 µm. The sample was made by vertical shale line cutting, and the samples' shapes were cylinders about 1 cm in length and 4 mm in width. A sample would be fixed in the device vertically and X-ray scanning was used to obtain the three-dimensional data volume of the core. In this paper, the software named Avizo 9.0 was used for the routine analysis of the three-dimensional data. According to the different X-ray absorption capacities of materials with different densities and thicknesses, the materials in the shale could be divided into three types: iron-bearing minerals, mineral matrix, and pores and microcracks.

Field Emission Scanning Electron Microscope (FE-SEM) Analysis
FE-SEM analysis was performed using a Quanta200F cold-field scanning electron microscope. The Quanta 200 F is a high-performance multi-purpose high-resolution FE-SEM with a high comprehensive field emission electron microscope. Resolution and environmental scanning electron microscopes (ESEM) are suitable as they have the advantages of sample diversity and can perform static and dynamic observations and analyses of various samples in high vacuum mode, low vacuum mode, and environmental scanning mode. Each sample was cut along the direction perpendicular to the bedding and polished with an argon ion beam to produce a smooth and flat surface with little surface change. In a high vacuum with 30 kV, the highest resolution can reach 2 nm, while the low vacuum with 3 kV can reach 3 nm.

Nitrogen Adsorption
Nitrogen adsorption was performed using a Kangta AbsorbIQ gas adsorption analyzer. The experiment required the use of high-purity nitrogen, with the sample needing to be degassed before the test, as the extent of the degassing directly affects the test results. Generally, a vacuum degassing system is adopted. High-resolution physical adsorption, chemical adsorption, and vapor adsorption isotherms can be measured to obtain the most accurate estimates of the pore sizes, the surface areas, and the specific gas/solid interactions. The sample was grinded to 60 to 80 mesh particles and degassing at 150 • C for more than 12 h.

High-Pressure Mercury Injection
High-pressure mercury intrusion was performed using the Mike 9505 automatic mercury intrusion instrument. As mercury does not wet general solids, to make mercury enter the pores, an external pressure must be applied. The greater the external pressure, the smaller the radius of the pore that the mercury can enter. By measuring the amount of mercury entering the pores under different external pressures, the pore volume of the corresponding pore size can be found. The maximum pressure of the mercury porosimeter is about 228 MPa, and the diameters of the measurable holes range between 50 to 360 µm. These measurements are used to draw capillary pressure curves and can also be used to describe the characteristics of multiple reservoirs, especially the size distribution of pore bellows in porous media.

Feature of Fracture Morphology under Rock Slice
Mud shale is the main component of unconventional shale reservoirs in the second sub-member of the Sophomore stage. Through a microscope, we can see the clay structure, silt-containing clay structure, and silty-clay structure, among which shale is often seen. A small amount of scale fragments (two-petal species, mainly ostracods), and some silt particles can be seen. When the orientation of those detrital particles is more obvious, the direction of bedding can be identified, and some authigenic minerals are common ( Figure 2). For example, in pyrite particles, a large amount of flaky organic matter orientated in the direction of the bedding can often be seen. The development forms of microcracks are mainly hairline, dendritic, and fractured networks. They are mainly distributed in the bedding and develop in the layers. In the well-defined shale, some shell-bearing shale and silt-bearing shale also have a few microcracks. Among them, silt-bearing shale has the least proportion of microcracks. When the silt content is high, it is difficult to find microcracks under the microscope.
the smaller the radius of the pore that the mercury can enter. By measuring the amount of mercury entering the pores under different external pressures, the pore volume of the corresponding pore size can be found. The maximum pressure of the mercury porosimeter is about 228 MPa, and the diameters of the measurable holes range between 50 to 360 μm. These measurements are used to draw capillary pressure curves and can also be used to describe the characteristics of multiple reservoirs, especially the size distribution of pore bellows in porous media.

Feature of Fracture Morphology under Rock Slice
Mud shale is the main component of unconventional shale reservoirs in the second sub-member of the Sophomore stage. Through a microscope, we can see the clay structure, silt-containing clay structure, and silty-clay structure, among which shale is often seen. A small amount of scale fragments (two-petal species, mainly ostracods), and some silt particles can be seen. When the orientation of those detrital particles is more obvious, the direction of bedding can be identified, and some authigenic minerals are common ( Figure  2). For example, in pyrite particles, a large amount of flaky organic matter orientated in the direction of the bedding can often be seen. The development forms of microcracks are mainly hairline, dendritic, and fractured networks. They are mainly distributed in the bedding and develop in the layers. In the well-defined shale, some shell-bearing shale and silt-bearing shale also have a few microcracks. Among them, silt-bearing shale has the least proportion of microcracks. When the silt content is high, it is difficult to find microcracks under the microscope.

Type, Occurrence, and Scale of Micro-CT Scan Cracks
At the micrometer scale, the bedding phenomenon of the shale samples is not significant, but we can see from the three-dimensional display images that the microcracks are almost parallel and spread out in the layers. Excluding other influencing factors, we conclude that the microcracks we observe are bedding-related microcracks, which either extend along the bedding or the main body extends along the shale bedding, but intersects with it at a low angle. This is fully illustrated in the sample FY1-13 ( Figure 3). mudstone, networked micro-fractures in mudstone.

Type, Occurrence, and Scale of Micro-CT Scan Cracks
At the micrometer scale, the bedding phenomenon of the shale samples is not significant, but we can see from the three-dimensional display images that the microcracks are almost parallel and spread out in the layers. Excluding other influencing factors, we conclude that the microcracks we observe are bedding-related microcracks, which either extend along the bedding or the main body extends along the shale bedding, but intersects with it at a low angle. This is fully illustrated in the sample FY1-13 ( Figure  3).

Types of Microcracks Observed under a Scanning Electron Microscope
According to our observations, there are mainly the following types of microfractures in the shale of the Da'anzhai section: structural microfractures, diagenetic shrinkage fractures of clay minerals, marginal shrinkage fractures of organic matter, and microfractures inside mineral particles. Among them, the structural microcracks that are of greatest significance are the extension of microcracks at the thin slice scale to the nanoscale. We have analyzed and counted them, and other types of microfractures are either too underdeveloped, or not typical enough, or even unable to be called fractures. While they have little impact on the reservoir, we nonetheless present here a brief introduction to them ( Figure 4).

Types of Microcracks Observed under a Scanning Electron Microscope
According to our observations, there are mainly the following types of microfractures in the shale of the Da'anzhai section: structural microfractures, diagenetic shrinkage fractures of clay minerals, marginal shrinkage fractures of organic matter, and microfractures inside mineral particles. Among them, the structural microcracks that are of greatest significance are the extension of microcracks at the thin slice scale to the nanoscale. We have analyzed and counted them, and other types of microfractures are either too underdeveloped, or not typical enough, or even unable to be called fractures. While they have little impact on the reservoir, we nonetheless present here a brief introduction to them ( Figure 4). mudstone, networked micro-fractures in mudstone.

Type, Occurrence, and Scale of Micro-CT Scan Cracks
At the micrometer scale, the bedding phenomenon of the shale samples is not significant, but we can see from the three-dimensional display images that the microcracks are almost parallel and spread out in the layers. Excluding other influencing factors, we conclude that the microcracks we observe are bedding-related microcracks, which either extend along the bedding or the main body extends along the shale bedding, but intersects with it at a low angle. This is fully illustrated in the sample FY1-13 ( Figure  3).

Types of Microcracks Observed under a Scanning Electron Microscope
According to our observations, there are mainly the following types of microfractures in the shale of the Da'anzhai section: structural microfractures, diagenetic shrinkage fractures of clay minerals, marginal shrinkage fractures of organic matter, and microfractures inside mineral particles. Among them, the structural microcracks that are of greatest significance are the extension of microcracks at the thin slice scale to the nanoscale. We have analyzed and counted them, and other types of microfractures are either too underdeveloped, or not typical enough, or even unable to be called fractures. While they have little impact on the reservoir, we nonetheless present here a brief introduction to them ( Figure 4).

Joint Evaluation of Pore Structure by High-Pressure Mercury Intrusion and Nitrogen Adsorption
According to the International Union of Pure and Applied Chemistry (IUPAC) classification, the measured nitrogen adsorption curve belongs to a type IV loop, with the hysteresis loop similar to H2 and H4. It has the characteristics of flaky pores and ink bottlelike throat, with small pore packing, point B, and multi-layer adsorption. The capillary agglomerates multiple processes, reflecting the existence of micropores and mesopores, and finally, no platform appears, indicating the existence of macropores, which may be caused by microcracks ( Figure 5A). In N 2 adsorption, the NLDFT model pays more attention to the characterization of micro-mesopores, and the BJH model pays more attention to the characterization of meso-macropores. In this paper, N 2 adsorption and high-pressure mercury injection were used to characterize the pore size, and the range of pore size selected by high-pressure mercury injection was greater than 50 nm. Therefore, NLDFT model, which focuses more on micro-mesopore characterization, was selected for characterization. We used the NLDFT model to explain its pore size distribution, and found that the pore size is mainly distributed between 1 to 2 nm ( Figure 5B). The combination of mercury intrusion and nitrogen adsorption can measure the distribution of full-aperture pores and cracks from 0.8 to 1500 nm, with the high-pressure mercury intrusion being used to evaluate the development of microcracks. The flaky clay is plastic, and high pressure will produce cracks. There are errors in interpretation. Cracks and microcracks generally have a width of more than 100 nm, and a macro-hole radius of more than 50 nm is just in the macro-hole interval ( Figure 6A,B). Therefore, the macro-hole is used to evaluate the proportion of microcracks, and the development degree of microcracks should not exceed this value.

Joint Evaluation of Pore Structure by High-Pressure Mercury Intrusion and Nitrogen Adsorption
According to the International Union of Pure and Applied Chemistry (IUPAC) classification, the measured nitrogen adsorption curve belongs to a type IV loop, with the hysteresis loop similar to H2 and H4. It has the characteristics of flaky pores and ink bottlelike throat, with small pore packing, point B, and multi-layer adsorption. The capillary agglomerates multiple processes, reflecting the existence of micropores and mesopores, and finally, no platform appears, indicating the existence of macropores, which may be caused by microcracks ( Figure 5A). In N2 adsorption, the NLDFT model pays more attention to the characterization of micro-mesopores, and the BJH model pays more attention to the characterization of meso-macropores. In this paper, N2 adsorption and high-pressure mercury injection were used to characterize the pore size, and the range of pore size selected by high-pressure mercury injection was greater than 50 nm. Therefore, NLDFT model, which focuses more on micro-mesopore characterization, was selected for characterization. We used the NLDFT model to explain its pore size distribution, and found that the pore size is mainly distributed between 1 to 2 nm ( Figure 5B). The combination of mercury intrusion and nitrogen adsorption can measure the distribution of full-aperture pores and cracks from 0.8 to 1500 nm, with the high-pressure mercury intrusion being used to evaluate the development of microcracks. The flaky clay is plastic, and high pressure will produce cracks. There are errors in interpretation. Cracks and microcracks generally have a width of more than 100 nm, and a macro-hole radius of more than 50 nm is just in the macro-hole interval ( Figure 6A,B). Therefore, the macro-hole is used to evaluate the proportion of microcracks, and the development degree of microcracks should not exceed this value.

Comprehensive Qualitative Evaluation of Multi-Scale Microcracks
Through multi-scale analyses of the shale from the Da'anzhai Member, we have observed numerous micro-scale fractures. Generally speaking, the microfractures identified at the slice scale are similar to those observed at the core scale and are mainly microfractures caused by tectonic stress, followed by smaller clay mineral diagenetic shrinkage joints, organic matter edge shrinkage joints, and microcracks inside mineral particles.
Through the combined observation of single-biased and orthogonal light, we found that the structural microfractures are mainly open, developed in the bedding-developed

Comprehensive Qualitative Evaluation of Multi-Scale Microcracks
Through multi-scale analyses of the shale from the Da'anzhai Member, we have observed numerous micro-scale fractures. Generally speaking, the microfractures identified at the slice scale are similar to those observed at the core scale and are mainly microfractures caused by tectonic stress, followed by smaller clay mineral diagenetic shrinkage joints, organic matter edge shrinkage joints, and microcracks inside mineral particles.
Through the combined observation of single-biased and orthogonal light, we found that the structural microfractures are mainly open, developed in the bedding-developed shale, with straight shapes and mainly bedding. In some areas, interlayer fractures have developed and are connected into dendrites and nets in fine-grained mud shale deposits, since mud shale is composed of a large number of lamellar clay minerals which have cleavage. The development of most microcracks is bedding-compliant, and we believe this may be related to the sheet-like arrangement of clay minerals ( Figure 2B-I). However, in other lithologies without obvious bedding, such as siltstone and fine sandstone, the microcracks are not developed. In addition, in the calcareous cement layers distributed along the beds in some mudstones and the calcareous cement near some organic fragments, microcracks (laminas) are particularly developed, which may have formed by stress. Under an electron microscope, the widths of the structural microfractures usually range between 0.1 and 3 µm, the fracture surfaces are curved, fully opened, mainly tensile, and often extend along the edges of the clay mineral flakes and brittle mineral particles ( Figure 4A,B). They occasionally cut through particles, which can be simply divided into marginal fractures and interlaminar fractures according to their development positions. They are the nanoscale extension of microfractures observed at the slice scale, which play a major role in the quality of shale reservoirs. The communication effect is conducive to the improvement of its porosity and permeability, which may be caused by the stress generated by the source rock overpressure and the regional tectonic stress. Based on micro-CT imaging, we can clearly see that the seams related to the bedding are particularly developed, with larger openings, wider extensions, and some interlayer seams that cross and extend, constituting a threedimensional interconnected hole-slit system ( Figure 3C,D). We believe this is sufficient to show that the cracks we have identified are mainly bedding-related.
The other three fracture types have different degrees of development in different samples. The clay mineral composition of the Da'anzhai section contains more than 50% of illite, and more than 25% of the Imonite mixed layer. In the process of converting illite to montmorillonite, the syneresis of clay minerals will occur. During this process, a large number of diagenetic shrinkage microfractures will form, mainly in the clay mineral flakes. The edges without structural microfractures are rough, their extension length is short, and the fracture width is narrow. These are mainly less than 100 nm, concentrated in the range of 30 to 80 nm ( Figure 4E). These appear less frequently than structural microcracks, and only in large numbers in a few samples. They are also identified as clay mineral interlayer pores in other pore and fracture classifications, and their pore size distribution interval is much smaller than that of conventional microfractures.
Organic matter edge shrinkage joints are formed by the volume shrinkage of organic matter during the thermal evolution associated with the burial and conversion of organic matter to hydrocarbons. They generally occur at the edges where organic matter and minerals are in contact, forming a slit-like shape. In the past, some researchers believed that these may be an important storage space for shale gas. In our observation of organic matter, their occurrence frequency is very high, and the fracture width is generally about 30 to 100 nm ( Figure 4C,D).
The common minerals in the shale of the Da'anzhai Member are terrigenous quartz, feldspar debris, and some shell particles. In addition to some structural microcracks that are prone to occur on the edges of minerals, some microcracks also occur inside mineral particles. The widths of the fractures are very small and they are difficult to find, being more common in shells, quartz, and pyrite particles. Since shell limestone is also an important potential reservoir in the Da'anzhai Member, the pore structure inside shell particles will have a certain reference value, hence some examples are presented ( Figure 4F).

Quantitative Characteristics of Micro-CT Microcracks
Micro-CT images were segmented and reconstructed using Avizo software. Firstly, the units of interest are selected and the image is segmented by using the threshold tool in Avizo software. The brightness of the high-resolution image is proportional to the atomic number of the sample components, which can clearly show the difference in the gray value of shale pores and cracks, organic matter and inorganic mineral matrix. According to the above principle, manual threshold method is adopted to extract the pores and cracks, organic matter and minerals respectively, and display their three-dimensional spatial structure. Then the Avizo Pore Network Model (PNM) module software was used to analyze the pore structure. The data type stored by the PNM module represents a grid composed of multiple linear lines in three-dimensional space. The branches or endpoints of the grid represent pores, and the straight lines connecting the pores are called throats. For each pore and throat, parameters such as radius, throat length, and coordination number can be calculated. Using the model that comes with the software, we extracted the relevant parameters of the pores. We calculated the average porosity of the sample to be 1.49%. It is almost the same as the hole and seam distribution observed in our images and the actual test data. The volume and area percentages of the holes and seams were calculated, and the statistical results are as follows. The distribution first shows a two-peak situation, indicating that the sample is a double-porosity medium with small micropores and large cracks. We believe that the large pores identified by the software are the cracks that can be seen, but the pores with an equivalent diameter of 50 to 100 nm only account for 16.9% of the entire pore volume (Figure 7). This seems to be very low compared to the other peak, but the ratio is still high. Due to the problem of threshold segmentation, some of the fractures are divided into many disconnected points (as can be seen from the model), with the actual contribution of the fractures to the porosity being much greater than this value. We can see from the three-dimensional display of the pore model that the pore distribution in this sample is mainly related to pores of a size between 1 to 2 µm and 2 to 5 µm (Figure 8). They contribute more than 71% of the pore volume and 81% of the pore surface area, and show that the peak of the pore size distribution of 1 to 2 µm is higher and that the software cannot identify pores less than 1 µm in size. This means that although we have identified a large number of pores ranging between 1 to 2 µm, it may mean that there are more pores less than 1 µm in size that have not been identified. This is also one of the limitations of this method; however, there appears to be no issue when analyzing pores 1 µm or more in size. ranging between 1 to 2 μm, it may mean that there are more pores less than 1 μm in size that have not been identified. This is also one of the limitations of this method; however, there appears to be no issue when analyzing pores 1 μm or more in size.

Quantitative Identification of Microcracks by Scanning Electron Microscopy
The microfractures observed under an electron microscope are all open, which is of great significance to the improvement of the porosity and permeability of shale reservoirs. Among them, a small number of fractures become channels for the migration and accumulation of hydrocarbons, which are contaminated by the migration of organic matter. Therefore, parameters such as the width and length of the fractures are critical to their connectivity. We have conducted statistical analyses on the distribution of the widths of the microcracks and found that they are mainly distributed between 200 to 400 nm, followed by between 400 to 1000 nm, and the number of microcracks can be ignored. To measure the degree of development of the microcracks, their face ratios in the electron microscope field of view is a parameter worth assessing ( Figure 9A). We used a magnification of 500 to 2000 times to observe and photograph microcracks. After calculating their face ratios, we found that they are mainly distributed below 0.8%. Although this value is not high, it is indeed an objective microscopic channel for oil and gas migration and accumulation and is of great potential significance to improving shale storage and permeability ( Figure 9B).

Quantitative Identification of Microcracks by Scanning Electron Microscopy
The microfractures observed under an electron microscope are all open, which is of great significance to the improvement of the porosity and permeability of shale reservoirs. Among them, a small number of fractures become channels for the migration and accumulation of hydrocarbons, which are contaminated by the migration of organic matter. Therefore, parameters such as the width and length of the fractures are critical to their connectivity. We have conducted statistical analyses on the distribution of the widths of the microcracks and found that they are mainly distributed between 200 to 400 nm, followed by between 400 to 1000 nm, and the number of microcracks can be ignored. To measure the degree of development of the microcracks, their face ratios in the electron microscope field of view is a parameter worth assessing ( Figure 9A). We used a magnification of 500 to 2000 times to observe and photograph microcracks. After calculating their face ratios, we found that they are mainly distributed below 0.8%. Although this value is not high, it is indeed an objective microscopic channel for oil and gas migration and accumulation and is of great potential significance to improving shale storage and permeability ( Figure 9B).

Comprehensive Quantitative Characterization of Multi-Scale Microfractures
Although we have obtained hole-fracture distribution information, we face new problems given that we intend to describe and characterize the microfractures in detail. A method such as micro-CT gives a too-low resolution in terms of its accuracy in recognizing microcracks, with the connected three-dimensionally distributed microcracks being identified as many isolated small holes. This leads to a decrease in the contribution of cracks to porosity and an increase in the contribution of porosity to porosity. In addition, the scale of various experiments is different, and it is difficult to determine the proportion of microfractures in the pore-fracture system.
On this basis, a comprehensive quantitative characterization of multi-scale microfractures was carried out. First, the combined quantitative characterization of highpressure mercury intrusion and nitrogen adsorption can measure the distribution of pores

Comprehensive Quantitative Characterization of Multi-Scale Microfractures
Although we have obtained hole-fracture distribution information, we face new problems given that we intend to describe and characterize the microfractures in detail. A method such as micro-CT gives a too-low resolution in terms of its accuracy in recognizing microcracks, with the connected three-dimensionally distributed microcracks being identified as many isolated small holes. This leads to a decrease in the contribution of cracks to porosity and an increase in the contribution of porosity to porosity. In addition, the scale of various experiments is different, and it is difficult to determine the proportion of microfractures in the pore-fracture system.
On this basis, a comprehensive quantitative characterization of multi-scale microfractures was carried out. First, the combined quantitative characterization of high-pressure mercury intrusion and nitrogen adsorption can measure the distribution of pores over a range of 0.8 to 1500 nm. The pore size data interpreted by nitrogen adsorption give values below 50 nm, and the data interpreted by high-pressure mercury injection provide values above 50 nm. These were sampled and statistics at uniform intervals were undertaken to obtain the total pore size distribution of the shale ( Figure 10A). It is found that the shale pore size has two peaks, both occupying the mesopore interval (2 to 50 nm), meaning the pore volume is mainly provided by mesopores. These peaks were at 1 to 10 nm and 10 to 100 nm, indicating that there are multiple types of pores. There were also multiple peaks above 100 nm, between 100 and 1000 nm, and between 1000 to 10,000 nm, but the ratios were small. This shows that the shale pore structure is complex, and multi-scale pores and fractures develop ( Figure 10B). Observed by light and electron microscopes, microfractures are basically macro-pores with a diameter greater than 50 µm. From the point of view of pore size distribution, the error associated with using macropores to estimate the proportion of microfractures will not be too large. However, due to the low overall diagenesis of the continental shale in the study area, there are still some primary pores and clay mineral pores in the shale that will be counted among them.

Comprehensive Quantitative Characterization of Multi-Scale Microfractures
Although we have obtained hole-fracture distribution information, we face new problems given that we intend to describe and characterize the microfractures in detail. A method such as micro-CT gives a too-low resolution in terms of its accuracy in recognizing microcracks, with the connected three-dimensionally distributed microcracks being identified as many isolated small holes. This leads to a decrease in the contribution of cracks to porosity and an increase in the contribution of porosity to porosity. In addition, the scale of various experiments is different, and it is difficult to determine the proportion of microfractures in the pore-fracture system.
On this basis, a comprehensive quantitative characterization of multi-scale microfractures was carried out. First, the combined quantitative characterization of highpressure mercury intrusion and nitrogen adsorption can measure the distribution of pores over a range of 0.8 to 1500 nm. The pore size data interpreted by nitrogen adsorption give values below 50 nm, and the data interpreted by high-pressure mercury injection provide values above 50 nm. These were sampled and statistics at uniform intervals were undertaken to obtain the total pore size distribution of the shale ( Figure 10A). It is found that the shale pore size has two peaks, both occupying the mesopore interval (2 to 50 nm), meaning the pore volume is mainly provided by mesopores. These peaks were at 1 to 10 nm and 10 to 100 nm, indicating that there are multiple types of pores. There were also multiple peaks above 100 nm, between 100 and 1000 nm, and between 1000 to 10,000 nm, but the ratios were small. This shows that the shale pore structure is complex, and multiscale pores and fractures develop ( Figure 10B). Observed by light and electron microscopes, microfractures are basically macro-pores with a diameter greater than 50 μm. From the point of view of pore size distribution, the error associated with using macropores to estimate the proportion of microfractures will not be too large. However, due to the low overall diagenesis of the continental shale in the study area, there are still some primary pores and clay mineral pores in the shale that will be counted among them. Thin slices can allow macroscopic microcracks above 5 µm to be identified, and electron microscopy can identify microcracks below this size, with some overlap between the observable scales of these methods. According to the statistics of fracture width, the widths of microcracks are not arbitrary. The slit width of the sheet is mainly 10 to 100 µm, and the aperture ratio is about 1%. The slit widths under the electron microscope are mainly 100 to 1000 µm, and the aperture ratio is about 0.2%. The distribution of slit widths does not overlap, indicating that the observed objects are of different scales, so they can complement each other. Therefore, a full-scale quantitative evaluation of microcracks can be performed by combining thin slices and electron microscope crack identification methods. Due to different magnifications used under microscopes, differing minimum crack widths can be observed, i.e., the higher the magnification, the wider the scale of the microcrack that can be observed, and the greater the number of microcracks observed, the higher the calculated face rate. Therefore, the large image stitching function of the NIS-Elements software was used to continuously stitch multiple microscope images at high magnifications into a single large image. The image obtained in this way not only shows more microscopic microfractures, but also avoids the influence of heterogeneity of microfracture development in shale.
For example, using a 5× objective lens and a 10× eyepiece, stitching slices of the full field of view for statistical analysis, the average aperture ratio of microfractures is 1.08%, and the aperture ratio of most shale samples is mainly distributed below 0.7%. However, there are still some samples with exceptionally developed microcracks, with face ratios between 1.5% and 3%, and the standard deviation reaches 1.163%, which is higher than the average value, indicating that the data are very discrete. It also shows that the degree of differentiation of microfracture development is very high. We evaluated the development of microcracks in the microstructure by randomly taking more than ten photos of each sample under the electron microscope at magnifications of 400 to 2000 times. The average aperture ratio of the nano-scale microcracks under the electron microscope is 0.53%, and the difference is not large, with all of them concentrated below 0.8%, and mostly below 0.2%. We believe that the development of these nanoscale microfractures will play a positive role in shale reservoir properties. The microfracture contribution rate (the ratio of fracture porosity to total porosity), is calculated according to the porosity of different samples, including the total porosity of shale in the Yuanba area of northern Sichuan and the Fuling area of eastern Sichuan, which are 5% and 7%, respectively. Then, the microfractures' contribution rates for different regions and wells are calculated. The data are highly dispersed, and in general, microfractures account for about 20% of the total pore volume ( Figure 11). microcrack that can be observed, and the greater the number of microcracks observed, the higher the calculated face rate. Therefore, the large image stitching function of the NIS-Elements software was used to continuously stitch multiple microscope images at high magnifications into a single large image. The image obtained in this way not only shows more microscopic microfractures, but also avoids the influence of heterogeneity of microfracture development in shale.
For example, using a 5× objective lens and a 10× eyepiece, stitching slices of the full field of view for statistical analysis, the average aperture ratio of microfractures is 1.08%, and the aperture ratio of most shale samples is mainly distributed below 0.7%. However, there are still some samples with exceptionally developed microcracks, with face ratios between 1.5% and 3%, and the standard deviation reaches 1.163%, which is higher than the average value, indicating that the data are very discrete. It also shows that the degree of differentiation of microfracture development is very high. We evaluated the development of microcracks in the microstructure by randomly taking more than ten photos of each sample under the electron microscope at magnifications of 400 to 2000 times. The average aperture ratio of the nano-scale microcracks under the electron microscope is 0.53%, and the difference is not large, with all of them concentrated below 0.8%, and mostly below 0.2%. We believe that the development of these nanoscale microfractures will play a positive role in shale reservoir properties. The microfracture contribution rate (the ratio of fracture porosity to total porosity), is calculated according to the porosity of different samples, including the total porosity of shale in the Yuanba area of northern Sichuan and the Fuling area of eastern Sichuan, which are 5% and 7%, respectively. Then, the microfractures' contribution rates for different regions and wells are calculated. The data are highly dispersed, and in general, microfractures account for about 20% of the total pore volume ( Figure 11).

Controlling Factors of Microfractures Development
Through our analysis of the previously calculated microfracture contribution rate, we can see that the degree of the microfractures' influence on the development of porosity is still very large. Generally speaking, the influencing factors of the development of microfractures in shale include the brittle mineral content, lithological combination, and regional structure. Therefore, we will discuss these three factors to determine which controls the development of microfractures.
We use the brittleness index (BI) to evaluate the evaluation index of brittle minerals. There are many calculation methods, with generally the method of rock mineralogy being used. Quartz, carbonate minerals account for percentage of quartz, carbonate minerals, and clay minerals. The percentage of these is the brittleness index. During our analysis, we found that in addition to these minerals, pyrite is also a component of our samples. However, pyrite is also a brittle mineral and is not accounted for. Therefore, we used the content of other minerals, except for clay minerals, to describe the brittleness index of shale, and calculated the brittleness index of each sample separately, and hence determined the brittleness index, porosity, fracture porosity, and fracture contribution rate (Table 1). These parameters were then, respectively, analyzed for correlations. We analyzed the correlation between the three factors of total porosity, fracture porosity, and the contribution rate of fractures to porosity presented in Table 1 and the brittleness index. The analysis results are shown in Figure 12. From the results of the correlation analysis, these three parameters are seen to be very similar, and the resulting graphs are close to each other. However, from these results, we cannot prove that the total porosity, fracture porosity, and fracture-to-porosity of the rock have a significant correlation with the brittleness index of the rock. and clay minerals. The percentage of these is the brittleness index. During our analysis, we found that in addition to these minerals, pyrite is also a component of our samples. However, pyrite is also a brittle mineral and is not accounted for. Therefore, we used the content of other minerals, except for clay minerals, to describe the brittleness index of shale, and calculated the brittleness index of each sample separately, and hence determined the brittleness index, porosity, fracture porosity, and fracture contribution rate (Table 1). These parameters were then, respectively, analyzed for correlations. We analyzed the correlation between the three factors of total porosity, fracture porosity, and the contribution rate of fractures to porosity presented in Table 1 and the brittleness index. The analysis results are shown in Figure 12. From the results of the correlation analysis, these three parameters are seen to be very similar, and the resulting graphs are close to each other. However, from these results, we cannot prove that the total porosity, fracture porosity, and fracture-to-porosity of the rock have a significant correlation with the brittleness index of the rock. In terms of the lithological combination, for the samples examined, there is no mutation from mudstone to limestone or sandstone in terms of their lithology. In a In terms of the lithological combination, for the samples examined, there is no mutation from mudstone to limestone or sandstone in terms of their lithology. In a previous article, it was found that a small amount of biological shells and sandstone bands exist in some samples of the fracture morphology. In shell-bearing shale and silt-bearing shale, a small amount of microcracking can be seen, but their proportions are less than that of shale and cannot play a significant role in the development of microfractures.
In terms of regional structural distribution, wells YB21 and YB102 are located in structurally stable areas, and large regional faults are not developed in their vicinity. The porosity in this area is well developed, but the contribution rate of microcracks to porosity is not high, with the microcracks accounting for about 10% of the total pores. Wells YL4, YL171, and YL176 are located in structurally developed areas, with fault-folded belts and high-steep belts. Their porosity is similar, but their microfractures lead to very high porosity, with the microfractures accounting for up to 40% to 90% (Figure 13). Through the above discussion of the contribution rate of microfractures, we gain the following basic understanding that the contribution rate of microfractures is less correlated with the content of brittle minerals, and the contribution of microfractures is higher in areas with developed structures than in areas with underdeveloped structures (Figure 14). The development of shale microfractures is mostly shale-related microfractures, and their degree of development is related to the development of the shale, and the regional tectonic environment. the above discussion of the contribution rate of microfractures, we gain the following basic understanding that the contribution rate of microfractures is less correlated with the content of brittle minerals, and the contribution of microfractures is higher in areas with developed structures than in areas with underdeveloped structures (Figure 14). The development of shale microfractures is mostly shale-related microfractures, and their degree of development is related to the development of the shale, and the regional tectonic environment. Figure 13. Relationship between the contribution of microfracture porosity and the degree of regional tectonic development.  porosity, with the microfractures accounting for up to 40% to 90% (Figure 13). Through the above discussion of the contribution rate of microfractures, we gain the following basic understanding that the contribution rate of microfractures is less correlated with the content of brittle minerals, and the contribution of microfractures is higher in areas with developed structures than in areas with underdeveloped structures (Figure 14). The development of shale microfractures is mostly shale-related microfractures, and their degree of development is related to the development of the shale, and the regional tectonic environment. Figure 13. Relationship between the contribution of microfracture porosity and the degree of regional tectonic development.

Conclusions
Through thin slice section, micro-CT scanning, scanning electron microscopy analysis, and nitrogen adsorption and high-pressure mercury intrusion measurements, we analyzed the characteristics of microfractures in the Da'anzhai shale from Lower Jurassic artesian wells in northeastern Sichuan. The conclusions are as follows: (1) Four types of fractures appear to have developed in the shale of the Da'anzhai section: mainly microfractures caused by tectonic stress, followed by smaller clay mineral diagenetic shrinkage fractures, organic marginal shrinkage fractures, and microfractures inside mineral particles. Among them, structural fractures and organic matter contraction fractures are the main development types and are of significance for shale reservoir development and seepage.
(2) The structural microfractures are mainly open and have developed in the beddingdeveloped shale, with straight shapes and bedding orientation. In some areas, interlayer seams are developed and formed into dendritic and net-like shapes, with curved seam surfaces, fully opened, and mainly tensile. Organic matter cracks often develop on the edge of the contacts between organic matter and minerals, showing a slit-like shape. Through micro-CT imaging, we clearly see that the seams related to the bedding are particularly developed, with larger openings, wider extensions, intersecting and expanding, forming a three-dimensional interconnected hole-slit system.
(3) Multi-scale stitching allows the calculation of the contribution rate of microfractures for different regions and wells. Generally, microfractures account for about 20% of the total pore volume. However, the degree of development of microfractures in different structural regions varies greatly. Microfractures in fault-wrinkle belts and high-steep belts can account for up to 40% to 90%, and microfractures in areas with underdeveloped structures account for about 10% of the total pore volume.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to some basic research involves confidentiality.