Abstract

Tectonism is one of the major controlling factors of shale gas accumulation and enrichment in China. To explore the relationship between tectonism and composition and pore characteristics of shale reservoirs, this research carried out mineralogy tests, organic geochemistry tests, field emission scanning electron microscopy (FE-SEM) experiments, and low-pressure gas adsorption (LPGA, N2 and CO2) experiments on the shale samples of various deformation intensities from Southwestern China. Based on the FE-SEM image analyses, it can be found that there are large differences in pore characteristics in shale samples with different deformation intensities. The samples with strong deformation have more organic pores, mainly related to the clay-organic aggregates and rigid grains. Tectonism can cause organic matter (OM) and clay minerals to be mixed or OM to fill in the clay layers, resulting in the retention of some organic pores. It is the presence of pressure shadows around the rigid grains that can resist tectonic extrusion and protect some organic pores. LPGA experiment results also show that micropore-specific surface areas and pore volumes of the samples with strong deformation are larger than those with weak deformation. The shale samples with strong deformation also have more microchannels and microfractures. Tectonism can also cause some micropores to become macropores; for example, tectonism can cause the rigid grains to slide and rotate, enlarging the dissolution pores at the edges of rigid grains. Shale samples with strong deformation have a smaller mesopore volume; but due to the presence of organic-clay aggregates, a larger mesopore-specific surface area embarks on these samples. According to fractal dimension calculations, it is found that in strong deformed shale, more multiple dimensions of the pore system tend to represent rougher pore surfaces and more irregular shapes. Besides, rougher pore surfaces are eager to provide more adsorption sites and enhance the adsorption capacity of the deformed shale. This study investigates the relationship between tectonism and composition and pore characteristics of shale reservoirs and may promote understanding of the accumulation of shale gas in highly deformed areas.

1. Introduction

Shale gas, taken as a very important unconventional natural gas, is essential to oil and gas production. Shale gas is mainly stored in pores and fractures in the state of free gas or as gas adsorbed onto kerogen and clay particle surface [14]. The pore structure strikes a remarkable impact on the gas storage and flow capacity of shale [513]. Unlike the gas shale field in the United States, the marine shale in Southern China has experienced multiple tectonism, as tectonism has been regarded as one of the most important factors affecting the shale pore structure and gas occurrence state [4, 1417]. Nowadays, more and more Chinese researchers have approved such viewpoints that tectonics is the key factor for shale gas accumulation in shale, and they set about conducting extensive research on tectonics in gas shale area [1831]. However, the relationship between tectonism and composition and pore characteristics of shale reservoirs is still unclear.

Two main ways can be observed concerning the influence of tectonism on the reservoir, namely, the vertical burial and uplift of the stratum and the extrusion and tension on the plane [4, 18, 23, 26]. And tectonism can also cause changes in minerals, especially clay minerals and OM [3234]. Therefore, tectonism can change the pores and fractures of shale reservoirs by stress and heat stress, which then further alters the occurrence of shale gas [18, 23, 34, 35]. Multistage tectonism together with dramatic changes in the pore structure took place in some shale-deposited zones in Sichuan Basin, Southwestern China [21, 2327]. Therefore, it is important to figure out the relationship between tectonism and pore structure. It is worth noting that there are still many controversies about the relationship between tectonic stress and the pore volume and specific surface area of shale. Some scholars believe that tectonic stress will cause the reduction of pores; that is, as the tectonism increases, the degree of pore reduction will increase, especially for mesopores and macropores [4, 23]. However, it is argued among experts that the tectonism can also produce the cleavage domains and increase the pore volume and specific surface area [32, 36]. Ma et al. found that the Lujiaping shale developed three-dimensional connected pore systems composed of nanometer-sized intergranular pore spaces, aggregate pore spaces in clay flakes, and a pore network in the cleavage domains, which provides the main specific surface area and adsorption point [32]. Although researchers pay more and more attention to the role of tectonism in shale gas accumulation, there exist many blurry doubts and unsolved issues such as the mechanism of the influence of tectonism on reservoir pore structure, which needs further research.

Recently, some researchers have got down to notice natural tectonic deformed samples and conducted experiments to study the relationship between tectonism and shale composition and pore structure features through experiments [32, 3640]. Ju et al. made inner comparison among the material composition and pore characteristics of brittle deformation, ductile deformation, and brittle-ductile shale samples from various areas and found that brittle shear results in microfractures and large pores and thus has an impact on the desorption and percolation capability of shale gas; the ductile deformation increased the specific surface area and enhanced the adsorption capacity of deformed shale [36]. Liang et al. found that the tectonism has little effect on shale mineral composition, maturity, porosity, and total pore volume by comparing the deformed sample with the undeformed sample but has a great influence on specific surface area, adsorption capacity, and pore size distribution [38]. Zhu et al. pointed out that the pores and microfracture systems will be combined during the deformation process so that these open pores become good reservoir space and migration channel [37, 39]. This study is of far-reaching importance of guiding the future exploration and development of structural regions, especially the weaker one.

2. Materials and Methods

2.1. Sample Selection

All samples were collected from the Lower Cambrian Lujiaping Formation in the Changzhang region. Located in the Dabashan thrust belt and undergone multistage tectonic movements, the Changzhang region is the cradle of anticline structures [41, 42]. The samples were all taken from a symmetrical anticline, and the deformation intensity is weakened from the core to the two wings (Figure 1). Based on the macroscopic characteristics, including hardness, scratch, rubbing mirrors, and deformation intensity, the samples were divided into strong deformed shale and weak deformed shale. Strong deformed samples include CZ-4, CZ-5, and CZ6, and weak deformed samples include CZ-1 to CZ-3 and CZ-7 to CZ-8.

2.2. Petrophysical Test

The samples were crushed to below 60 mesh. TOC content test and mineral composition test were measured by a LECO CS230 carbon and sulfur analyzer and Rigaku D/Max-3B diffractometer, respectively. The experiments were conducted in accordance with Chinese national standards GB/T19145-2003 and SY/T5163-2010. The whole rock analysis and clay minerals were estimated by semiquantitative mineral spectrum. Both experiments were carried out at Experimental Research Center of East China Branch, Sinopec.

The equivalent vitrinite reflectance test () was measured by 3Y-Leica DMR XP microscopy equipped with a microphotometer, and the results were analyzed in the Key Laboratory of Coalbed Methane Resources and Reservoir Formation Process of the Ministry of Education (CUMT). Owing to the deficiency in vitrinite in the Lower Cambrian Lujiaping Formation shale, the is converted into vitrinite-like material reflectance () by calculation [43].

2.3. FE-SEM Experiment

The samples were processed by argon ion beam and then sprayed with 1 mm thick gold. The FE-SEM experiment was performed at the Jiangsu Design Institute of Geology for Minerals Resources, taking use of FE-SEM (Quanta 200F) equipped with an energy-dispersive spectrometer (EDS) to identify the mineral types and help identify pore types. These analyses were presented at the temperature of 24°C with the humidity level of 35%.

2.4. Low-Pressure Gas Adsorption (LPGA)

LPGA experiments were strictly performed at the Physical and Chemical Center of Beijing, according to the Chinese National Standard GB/T 19587-2004. All experimental samples were smashed to 40-60 mesh, and then, were automatically degassed under a vacuum of ~110° for more than 16 h. The instrument used was called the Micromeritics ASAP 2020 HD88 analyzer. N2 adsorption experiment was implemented in the bath temperature of -196°C (77 K) under the relative pressure () ranging from 0.009 to 0.995. And the bath temperature in CO2 adsorption experiment is 0°C (273 K), and the relative pressure () of CO2 adsorption range is between 0.0001 and 0.032. The experiment data of LPGA experiment were interpreted by the field density functional theory (DFT) model for pore size distribution of the shale.

2.5. Method of Fractal Calculation

The fractal dimensions were calculated based on the data obtained from LP-N2-GA experiment. The equation is as follows: where is the adsorbed gas volume at the equilibrium pressure (), is the monolayer coverage volume, is a characteristic constant, is the universal gas constant, is the absolute temperature, is the saturated vapor pressure, and the parameter of depends on fractal dimensions and the specific mechanism of gas adsorption.

Equation (1) can be written in log-log terms:

According to the fractal FHH theory, the regression straight slope () of can be used to calculate fractal dimensions :

3. Results

3.1. Geochemical and Mineralogical Characterization of Different Shale Samples

The geochemical and mineralogical data for all samples are shown in Table 1. The TOC contents of the Lujiaping shale scope from 1.01% to 1.79% with average values of 1.35%. The values between 3.93% and 4.05% are calculated with the mean of 3.98%, and all the samples are in the overmature stage.

The XRD data have proven that quartz and clay are examined to be the dominant minerals as their overall contents occupy more than 55%. Specifically, the content of quartz is usually relatively high, between 30.1% and 48.8% with average values of 40.1% (Table 1). All the Lujiaping shales boast a certain amount of carbonate, feldspar, and pyrite with greatly varied contents seen from different samples. It is worth noting that the clay mineral content ranged from 19.1% to 27.7% with a mean value of 24.2%. The illite and chlorite represent the primary clay minerals, and the former contents of which scope from 50% to 72% of total clay content with a mean value of 57.4%. The total amount of quartz and clay of the CZ-4 and CZ-6 with strong deformation can reach about 70%; in particular, the quartz content is much higher than other samples.

3.2. Pore Structure Analyses with the LPGA

Figure 2 shows the N2 adsorption-desorption isotherm of the selected deformed shale samples. Combining with previous studies, obviously, to summarize the pore morphology characteristics in shale can be inferred based on the type of N2 adsorption-desorption isotherm and the shape of the hysteresis loop [8, 11, 44]. In accordance with the IUPAC classification of the hysteresis loops, the N2 adsorption-desorption isotherm of the shale samples belongs to Types H2 and H3 [11, 4548]. It can be seen that the N2 adsorption-desorption isotherm of strong deformed samples has large hysteresis loops (Type H2), indicating that these samples contain a large number of ink bottle-shaped pores. The N2 adsorption-desorption isotherm of the weak deformed samples (CZ-2 and CZ-8) is similar to Type H3, indicating that these samples contain slit-like pores. These pores may be related to the plate structure of clay particles, which are usually flat in the SEM images. What should be paid attention is that the hysteresis loop of N2 adsorption is a superposition of multiple standard loops and is also a manifestation of the complex pore structure characteristics of the samples.

Figure 3 is the CO2 adsorption isotherms of all shale samples, typically as the I adsorption isotherms. Shale samples with various deformation intensities have a significant difference in carbon dioxide adsorption capacities. In CZ-6, CO2 adsorption is seen as the largest amount while the smallest amount in the CZ-1. In general, the shale samples with strong deformation have a larger adsorption amount. The relationship between the pressure and the CO2 adsorption amount of these samples can be converted into the differences in specific surface area and micropore volume.

The pore volume and specific surface area of these samples were calculated by the DFT model using the N2 and CO2 adsorption data (Table 2). Figure 4 shows that micropores and mesopores contribute to the major pore volume, which can reach almost 80%. The total pore volume is mainly contributed by pores smaller than 2 nm and 10-50 nm, and pores of other pore sizes also contribute a part of the pore volume. And it can be found that the pore volume contributed by the micropores in CZ-4, CZ-5, and CZ-6 with strong tectonic deformation is larger than that of the mesopores. Conversely, in other samples with weak deformation, the mesopores contribute more to the pore volume. Surprisingly, pores of 2-10 nm are not well developed in these samples and may be related to the tectonism.

Figure 5 shows that in all samples, the specific surface area gradually decreased as the pore size increased. The specific surface area is mainly contributed by the micropores, which provide up to 70% of the total specific surface area.

Although the total specific surface area of all samples is created by micropores, the micropores in the strong deformed shale samples provide more specific surface area.

3.3. Pore Morphology Analyses Based on FE-SEM Images

Shale has a variety of pore types, and each pore type has a complex morphology and origin [4, 813, 4951]. In recent years, some researchers have used FE-SEM experiments to observe the appearance and pore types of shale pores but rarely observed the pore characteristics in the deformed shale samples [9, 1113]. FE-SEM observation of the pore characteristics of the Lujiaping Formation shale reveals that there are significant variations in the pore characteristics of the samples among various deformation intensities. The tectonism has a significant influence on OM pores, interP pores, and microchannels and microfractures, and the differences among these pores are emphasized.

There are relatively few OM pores in most samples (Figures 6(a)–6(d)), which is consistent with previous research results, probably due to the overmature stage of the Lujiaping Formation and related to strong compaction of tectonic deformation [32, 37, 50]. The weak deformed shale is possessed with some scattered organic pores and microfractures at the edges of OM and clay minerals (Figure 6(a)). However, comparatively more OM pores in the strong deformed shale samples are mainly associated with the OM-clay aggregate and the supporting effect of rigid grains (Figures 6(c) and 6(d)). Multiple types of interP pores can be seen in the SEM images, which have different shapes and sizes and vary greatly (Figures 6(e)–6(h)). It shows that there are usually interP pores between OM and clay minerals (Figure 6(a)). These pores are generally large in length and can effectively improve the pore connectivity. InterP pores are always closely related to rigid grains (quartz, carbonate, feldspar, and mica), especially around rigid grains (Figures 6(e)–6(h)). There are always some triangular regions at the junction of rigid grains and other minerals. These interP pores are either irregular in shape or nearly rectangular outlines due to the structure of the rigid grains (Figures 6(f) and 6(h)).

SEM images also reveal that there are more porous OM pores between the pyrite framboids, although these organic pores are small, but the number is relatively large (Figure 6(c)). For one aspect, it may be due to the tectonic compaction that many pyrite grains are loose from the surrounding clay, resulting in some interP pores (Figure 6(j)). Microfractures and microchannels are typically linear, with lengths between nanometers and micrometers. For another aspect, microfractures are usually cut along the edges of grains, whereas microchannels are usually perpendicular to the plane (Figures 6(i)–6(l)). The microfractures generally feature better connectivity and larger size than those of the microchannels. The multistage tectonic movements developed more microchannels and microfractures. The strong deformed shale samples have more microchannels and microfractures. In the meantime, these microchannels and microfractures play a key role in shale gas migration and storage. They can increase porosity and permeability and contribute to later fracturing production [1, 51].

3.4. Fractal Characteristics of Pores in Deformed Shale

The fractal curves of the shale samples are similar, but there are some differences in the slopes of the liner fitting lines (Figure 7). The fractal dimension of all samples obtained from the FHH model using nitrogen adsorption is mentioned in Table 2. Each sample has two fractal dimensions, including at 0-0.5 and 0.5-1. The slopes of the correlation curves of the two parts are different, and the correlation coefficient is greater than 0.99. This phenomenon may be caused by different mechanisms of nitrogen adsorption at diverse stages. The linear regression values and and the correlation coefficients and of the correlation curve are included in Table 3.

represents the fractal dimension of the relatively small nanopore in the low-pressure stage (), while represents the fractal dimension of the relatively large nanopore in the high-pressure stage (). And is between 2.5216 and 2.6836, with an average of 2.6357. is between 2.6141 and 2.7443, and the average is 2.7122. According to Table 2 and Figure 7, it can be found that is usually larger than . Strong deformed samples generally have relatively larger and , suggesting that tectonism makes the pores in the shale more complex, especially for relatively large nanopores.

4. Discussion

4.1. Effect of Tectonism on Shale Composition

The primary compositions of organic-rich shale embrace OM, clay minerals, and brittle minerals, such as quartz, feldspar, carbonate, and pyrite. There are many factors which can affect the composition of shale, including material sources, sedimentary environment, diagenesis, and tectonism [2, 4, 7, 13, 36]. According to Table 1, all samples collected from the Lujiaping Formation shale in the Changzhang region generally contain relatively high TOC and brittle mineral contents. Previous studies have shown that brittle deformed shale (BDS) samples always contain relatively high TOC and brittle mineral content, while ductile deformed shale (DDS) samples generally have higher clay mineral content [3639]. However, this conclusion is not very consistent with the XRD results of the deformed samples in the Changzhang region. Shale samples with strong deformation generally have high quartz content and relatively low clay mineral content than samples with weak deformed shale samples (Figure 8).

In fact, tectonism can change the diagenetic stage by stress and heat stress and thus change the composition of shale. On the one hand, as the buried depth increases and the diagenesis strengthens, some unstable minerals will be transformed into other minerals [45]. For example, feldspar will be transformed into quartz and kaolinite under the action of organic acid [52]. On the other hand, extrusion and tension on the plane can change the stress structure of the reservoir by forming different structural styles, and the microstructure and chemical composition of minerals of the shale were changed to a large extent [23, 26, 31, 53]. Ramsay first analyzed the anticline fold and its internal strain characteristics, finding that the stratum above the neutral plane exhibits tensile stress and the stratum below the neutral plane exhibits compressive stress (Figure 9). During the process of stress extrusion to form an anticline, the part above the neutral surface of the Lujiaping shale will be subjected to strong tensile stress, which will change the composition of the shale [3335, 5254]. As shown in Figure 6, under tensile stress, more ductile clay and OM are compact in the space around the brittle minerals and increase the heterogeneity of the shale. Due to the tectonism, the hydrothermal fluid is filled into the microfracture formed by the tectonic deformation, which also may cause the shale with strong deformation to have high quartz content [37].

4.2. Effect of Tectonism on Pore Characteristics

Tectonic stress may have destroyed the mesopores to macropores within the OM in the deformed area, while OM-hosted micropores are retained due to their greater mechanical resilience [3335, 54]. And tectonism also can cause the shale to develop a “mylonite structure” in the cleavage domains, which has a pore system composed of nanoscale intergranular gap, lamella gap aggregate, and pore network, thus greatly increasing the specific surface area and adsorption capacity of the shale [32]. FE-SEM images show that strong deformed shale samples have more OM-hosted micropores (Figures 6(a)–6(d)). From the results of LPGA experiment, it can be found that strong deformed shale samples generally have a larger micropore specific surface area and a larger micropore volume (Figures 10(b) and 10(d)). Ju et al. also found that strong ductile deformation enables the pore size shift its distribution in shale, and a strong tectonic shear can improve the coal-bearing shale pore specific surface area [36]. These changes in pore characteristics are originated from the rheological behavior of the ductile minerals (clay minerals and OM). Tectonism can stretch the OM and fill the clay mineral layers to form more micropores (Figure 11, I). There is the presence of pressure shadows around the rigid grains (quartz, feldspar, dolomite, calcite, and pyrite) [33], which can resist the tectonic compression and also protect the surrounding OM pores (Figures 6(j) and 11, II).

During the compaction process, a large number of the interP pores and intraP pores can be destroyed; this is why the pore volume is reduced by up to 88% [49]. The late tectonic extrusion is essential to the interP pores and the intraP pores [23, 36, 37]. Although the intraP pores are relatively developed, considering the independent development and light connectivity of the intraP pores, they contribute less to the overall pore system [7, 13]. The most common interP pores are usually elongated strips along the edges of rigid grains (Figures 6(f) and 6(h)). Due to the dissolution, some small interP pores are formed at the edge of the rigid grains, ranging from a few tens to several hundreds of nanometers (Figure 6(e)–6(h)). The tectonism can cause the rigid grains to rotate and slide and further enlarge the dissolution pores at the rigid grains’ edge (Figure 11, III). Tectonism can also cause some mesopores to become macropores; for example, tectonism can connect the interP pores at the rigid grains’ edge to form microfractures (Figure 6(h)). Therefore, the sample with strong deformation has a larger mesopore specific surface area and a smaller mesopore volume (Figures 10(a) and 10(c)).

Shale deformation is classified in three main forms of brittle deformation, brittle-ductile deformation, and ductile deformation [3640]. Whether it is brittle deformation or ductile deformation, as the deformation intensity increases, both microchannels and microfractures increase [37, 38, 40]. Compared with the development characteristics of microchannels and microfractures of shale samples with various deformation intensities, some points may be concluded that samples with strong deformation have more microchannels and microfractures (Figures 6(i)–6(l)). Through linking with microfractures, tectonism is able to form the microfracture system (Figure 11, IV). Although some microchannels and microfractures will close under tectonic stress, the principal stress of different directions in multiple phases will cause the microfractures to unfold again [23]. Most noteworthy is that these microfractures and microchannels also change the mechanical properties and migration capabilities of shale, which is more conducive to future development [51].

4.3. Effect of Tectonism on Pore Fractal

To characterize the complexity of the pore network in shale reservoirs, some researchers reach this goal in a way to calculate fractal dimensions [8, 11, 26, 45]. N2 has different adsorption mechanisms under different pressure conditions, so the fractal dimension of the low pressure stage is usually used to indicate the surface roughness of the pore structure, while the fractal dimension of the high-pressure stage represents the irregularity of the pore structure [11, 47, 48]. The tectonism can change the relationship between the pores in the shale reservoir, making the pore system of the shale more complex [32, 3640]. The samples with strong deformation generally have a relatively large and than samples with weak deformation (Table 2 and Figure 12). The samples with strong deformation that always have a relatively large may be related to the mixing behavior of clay and OM. The OM-clay aggregate has many organic pores, and the organic pores have neither a smooth surface nor a regular shape, which also leads the shale with strong deformation having a large . The rough surface can also form some additional pores that increase the adsorption capacity of the shale. The irregularity of pore structure has an important influence on the spatial geometrical dimension, pore heterogeneity, and interconnectivity [11, 13, 49]. And multiple tectonism has an important influence on the pore network, which makes the nm-μm scale interP gap particularly developed [37, 40]. This also causes the shale with strong deformation to have a high .

5. Conclusion

This paper mainly studies the effects of tectonics on the composition and pore characteristics of shale reservoirs. A series of experiments were carried out on 8 deformed shale samples from the Lujiaping Formation in the northeastern of Chongqing, including mineralogy, organic geochemistry, low-pressure N2/CO2 adsorption, and FE-SEM. The major conclusions of this paper are as follows: (1)Tectonism has an important influence on the pores in the shale reservoir, especially OM pores, interP pores, microchannels, and microfractures. The strong deformed shale samples have more organic pores, larger interP pores, and more microfractures and microchannels(2)The organic-clay aggregate produced by the tectonism provides more specific surface area and increases the adsorption capacity of the deformed shale. However, due to the tectonism, some micropores become macropores, which also leads the strong deformed shale to have a smaller mesopore volume(3)Different structural parts have different stress environment, and different minerals have different mechanical properties, so different pore characteristics exist in different structural parts. Tectonism can cause the rheology of the clay mineral and OM around the rigid grains so that the rigid grains can protect the surrounding micropores. The tectonism can cause the rigid grains to slide and rotate and also enlarge the dissolution pores at the rigid grains’ edge(4)Tectonism makes the pore network in shale reservoirs more complex. The strong deformed samples generally have larger and than weak deformed samples. As the shale deformation intensity increases, the surface roughness and irregularity of the pore structure are enhanced

Data Availability

The original data used to support the findings of this study are available from the corresponding author ([email protected]) upon request.

Additional Points

Highlights. (1) The effects of tectonism on the composition of shale reservoirs are investigated to specify the mode of influence. (2) Analysis of the effect of tectonism on the composition is used for clarifying the influence of tectonism on the pore characteristics. (3) The tectonism will complicate the pore network in the shale reservoir as well as cast a crucial impact and have an important effect on the pore fractal dimension.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research is sponsored by the National Natural Science Foundation of China (Project Nos. 41872132 and 41802183), National Science and Technology Major Project (Project No. 2017ZX05035004-002), and Scientific Research Foundation of Key Laboratory of Coalbed Methane Resources and Reservoir Formation Process Ministry of Education (China University of Mining and Technology) (No. 2019-008). We really appreciate the help from other experimenters in the Jiangsu Design Institute of Geology for Minerals Resources labs.