Next Article in Journal
Development of Gold Nanoparticle Micropatterns for the Electrical Detection of Proteins
Next Article in Special Issue
Dynamics of Stress-Driven Two-Phase Elastic Beams
Previous Article in Journal
Photocatalytic Activity of Magnetic Nano-β-FeOOH/Fe3O4/Biochar Composites for the Enhanced Degradation of Methyl Orange Under Visible Light
Previous Article in Special Issue
Dynamic Behavior of Magnetically Affected Rod-Like Nanostructures with Multiple Defects via Nonlocal-Integral/Differential-Based Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Full-Scale Pore Structure Characteristics and the Main Controlling Factors of Mesoproterozoic Xiamaling Shale in Zhangjiakou, Hebei, China

1
School of Earth and Space Sciences, Peking University, Beijing 100871, China
2
Hebei Key Laboratory of Strategic Critical Mineral Resources, Hebei GEO University, Shijiazhuang 050031, China
3
State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Changping, Beijing 102249, China
4
Experimental Practice Teaching Centre of Hebei GEO University, Shijiazhuang 050031, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(2), 527; https://doi.org/10.3390/nano11020527
Submission received: 19 January 2021 / Revised: 14 February 2021 / Accepted: 15 February 2021 / Published: 18 February 2021
(This article belongs to the Special Issue Mechanics of Nanostructures and Nanomaterials)

Abstract

:
Nanoscale pore structure characteristics and their main controlling factors are key elements affecting the gas storage capacity, permeability, and the accumulation mechanism of shale. A multidisciplinary analytical program was applied to quantify the pore structure of all sizes of Xiamaling shale from Zhangjiakou, Hebei. The result implies that Mercury injection porosimetry (MIP) and low-pressure N2 curves of the samples can be divided into three and four types, respectively, reflecting different connectivity performances. The maximum CO2 adsorbing capacity increases with increasing total organic carbon (TOC) content, pore volume (PV), and surface area (SA) of the micropores are distributed in a three-peak type. The full-scale pore structure distribution characteristics reveal the coexistence of multiple peaks with multiple dominant scales and bi-peak forms with mesopores and micropores. The porosity positively correlates with the TOC and quartz content, but negatively correlates with clay mineral content. Organic matter (OM) is the main contributor to micropore and mesopore development. Smectite and illite/smectite (I/S) assist the development of the PV and SA of pores with different size. Illite promotes the development of the nanoscale PV, but is detrimental to the development of the SA. Thermal maturity controls the evolution of pores with different size, and the evolution model for the TOC-normalized PVs of different diameter scales is established. Residual hydrocarbon is mainly accumulated in micropores sized 0.3 to 1.0 nm and mesopores sized 40 nm, 2 nm and less than 10 nm. Since the samples were extracted, the pore space occupied by residual hydrocarbon was released, resulting in a remarkable increase in PV and SA.

1. Introduction

Shale gas is mainly adsorbed on the surface of mineral particles or organic matter (OM), or is stored in a free state in nanoscale pores or a dissolved state in kerogen and bitumen [1,2,3,4]. Nanoscale pores not only provide storage space for shale gas but also provide a geological basis for later reservoir exploitation. Therefore, nanoscale pore structure is an important element in the evaluation of the gas storage capacity and seepage performance of shale gas and a key indicator of whether shale gas has commercial exploitation value [1,5,6].
However, owing to the complexity and variability in shale pores, conventional reservoir porosity testing methods are far from sufficient for the requirements of characterization. To better characterize the pore structure of shale, a variety of technical means should be combined. Accordingly, a range of qualitative and quantitative methods have previously been applied by scholars. The qualitative methods include computerized tomography (CT) scanning, transmission electron microscopy (TEM), field emission scanning electron microscopy (FE-SEM), broad-ion beam SEM (BIB-SEM), and focused ion beam polishing in combination with SEM (FIB-SEM) [7,8,9,10,11]. These qualitative image observation technologies have previously been widely applied to observe pore sizes and distributions and the geometry and connectivity of pore networks. Image-based visible nanoscale pore information can be extracted via point counting in Image J software (V1.8.0.112, 22 July 2019) [12,13]. Quantitative methods, such as small angle and ultra-small angle neutron scattering, gas adsorption (CH4, CO2, N2), and mercury intrusion porosimetry (MIP), have been used to measure pore structure parameters at nm-scale resolution. To quantify pore structures of all sizes, combinations of MIP, CO2, and N2 adsorption have been successfully used by numerous scholars [2,6,14,15,16].
Under actual geological conditions, the nanoscale pore structure of shale is constantly changing due to the combined effect of diagenesis and the thermal transformation of OM and kerogen [17,18]. OM and kerogen are the material bases for the hydrocarbon generation of shale, and play a decisive role in the generation and development of organic pores [19]. According to Jiarvie et al. [20], the porosity of shale increases by 4.9% with organic carbon consumption of 35%. According to the simulation results of Modica and Lapierre [21], each gram of organic carbon generally contributes 0.8–0.85% to the pores of shale at the end of the gas generation of OM. Therefore, convertible OM results in pyrolysis and leads to the generation and discharge of hydrogen [22] at a certain thermal maturity. As large numbers of OM pores are generated, the residual hydrogen in the shale inevitably influences the pore system. The inorganic minerals in shale mainly include brittle minerals such as quartz and feldspar, whereas clay minerals are dominated by smectite, illite/smectite (I/S), and illite. These minerals determine the formation and evolution of a framework of mineral pores in shale during the process of diagenesis [1,22]. However, current research on factors affecting the full-scale pore structure primarily involves postmature or extremely mature marine gas generation shales, with a thermal maturity greater than 2.0% [6,23,24,25], whereas research on the main factors controlling the full-scale pore structure of low- and middle-mature marine shale with favourable aerogenesis is scarce. High- and over-mature shale is the product of low- and middle-mature shale during diagenesis and hydrocarbon generation, and immature shale with different pore structure characteristics inevitably generates high- and over-mature shale with different accumulation effects. In particular, low- and middle-mature shale formations with strong hydrocarbon generation ability contain large amounts of residual hydrocarbon, which leads to nonnegligible influences on the formation and development of the pore structure and subsequently changes the gas storage and filtration ability of shale, influencing exploitation and gas generation.
Therefore, based on the aforementioned status, a combination of MIP, low-pressure N2 and CO2 adsorption was applied (1) to characterize the complex full-scale pore size distribution of Xiamaling shale originating from Zhangjiakou, Hebei, China; (2) to reveal the relationship between the pore structure and affecting factors; and (3) to provide evidence for evaluating the reservoir performance of Xiamaling shale and the shale of different date and formation of the world. The MIP and N2 low pressure adsorption curve classification of this research can be used to evaluate the development degree, morphological characteristics, and connectivity of pores with different pore sizes in shale, and to show the contribution degree of pores with different pore sizes to PV and SA. This study provides important evidence for understanding and evaluating the desorption, diffusion, seepage, migration, and later exploitation of shale gas in different occurrence states.

2. Geological Setting

The research area is located in Zhangjiakou, Hebei (Figure 1). Xiamaling shale is widely distributed in the villages of Xiahuayuan, Zhaojiashan, Jingjitun, and Xiajiagou. The Zhaojiashan and Xiajiagou areas are the deposition centre of Xiamaling Formation, with a maximum cumulative thickness of 587 m and thinning from west to east, and the section exposure is almost complete, in which a standard section has developed. The Xiamaling Formation unconformably contacts the overlying Changlongshan Formation and underlying Tieling Formation and can be divided into four sections [26,27]. The thick of the third section is up to 350 m and contains black carbonaceous, silty and grey-black siliceous shale, as well as oil shale accumulated by red algae and brown algae, indicating that the Xiamaling shale in this region has superior hydrocarbon generation potential, with favourable shale gas exploration prospects. Additionally, according to the latest isotopic timing results, the Ximaling shale formation was deposited 1368 Ma years ago, and formed at Jixian system in Precambrian, indicating that the Ximaling shale is the oldest marine shale found in China, and it is one of the oldest ancient marine shales in the world [28,29].

3. Samples and Experiments

In total, 18 fresh outcrop shale samples from the six villages of Zhaojiashan, Jizhentun, Xiajiagou, Huangtugang, Jimingshan, and Xiahuayuan were collected for the precursory studies. The six sections all include organic-rich Xiamaling marine shale, and the locations are marked in Figure 1. The experimental samples were mainly black and carbonaceous shale. During the sample collection process, fresh samples were collected 2–3 m below the surface. As the experimental samples were collected from the outcrops of the Xiamaling shale in Zhangjiakou, any weathered surface material on the collected samples was removed, and the samples were thoroughly cleaned using ethanol to remove possible surface contamination. Each of the selected samples was split into several aliquots for analysis of the mineralogical composition, pore morphology, and pore size distributions, for organic petrographic analysis, and the extraction of residual hydrocarbon.

3.1. Basic Geochemistry and Mineralogy

The TOC content was measured by an LEC CS-230 carbon-sulphur analyser after treatment with HCl to remove carbonates according to the Chinese National Standard GB/T18602-2001 and GB/T19145-2003. Because no “true” vitrinite was found in the shale, vitrinite-like maceral reflectance was used to characterize the thermal maturity. The measurements were conducted on polished blocks under reflected light using a 3Y microphotometric system equipped with an oil-immersion objective lens and photometer. The vitrinite-like maceral reflectance was then converted to the equivalent vitrinite reflectance (EqRo) using the equation proposed by Xiao et al. [30].
X-ray diffraction (XRD) was performed using a Bruker D8 Advance X-ray diffractometer with a scan range of 3° to 85° (2θ) and a scan speed of 4°/min. The analysis was conducted on powdered bulk shales and clay fractions to determine the mineralogical composition and identify the clay minerals, respectively. The area under the curve for the major peaks of each mineral was used to estimate the relative mineral percentages of the shale with correction for Lorentz polarization [31].

3.2. Mercury Injection Porosimetry

Porosity and macropore characterization was performed on cubic samples sized approximately 1 × 1 × 1 cm3 using an AutoPore IV 9510 Micrometerics instrument (Micromeritics Instrument, Norcross, GA, USA) at 19 °C and a humidity of 18%. The bulk and skeletal densities used for determining the porosity of the sample could be obtained from the mercury intrusion measurements performed at 0 and 413 MPa, respectively. The measured spans of the pore diameter were approximately 3 nm–1000 μm, the mercury contact angle was 130°, and the surface tension was 485 dyne/cm. The accuracy of the mercury injection volume was 0.1 mL, and the macropore size distribution was calculated by means of the Washburn formula [32]. Each sample prepared for the MIP test was oven dried for 48 h at 110 °C to remove adsorbed moisture and volatile matter and then cooled to room temperature before measurements were obtained.

3.3. Residual Hydrocarbon Extraction

Prior to the Soxhlet extraction experiments, free water was removed from the shale samples for 24 h at 48 °C. Then, the samples were crushed and sieved by 80–120 meshes. Approximately 40 g of the measured sample grains was used for Soxhlet extraction for 72 h with a 25:2 vol/vol mixture of dicholoromethane and methanol. After the extracted liquid mixture was concentrated using a rotary evaporator, the mixture was transferred and weighed in a 4 cm3 glass vial, and the weight of the concentrated liquid mixture represented the content of residual hydrocarbon. The original samples and their corresponding extracted solid residues were placed in a vacuum oven and held for 24 h at 90 °C. Then, low-pressure N2 and CO2 adsorption measurements were obtained.

3.4. Low-Pressure N2 and CO2 Adsorption

Low-pressure N2 and CO2 adsorption were performed using a Quantachrome Autosorb-1 apparatus SA analyser at the State Key Laboratory of Heavy Oil Processing in China University of Petroleum, Beijing. The 2–3 g shale sample was crushed into grains of approximately 60–100 mesh (250–150 μm) to obtain nanoscale pore structure parameters. Prior to the low-pressure N2 and CO2 adsorption measurements, the selected samples were automatically degassed under vacuum for approximately 12 h at 110 °C to remove air, free water, and other gases. The conditions of N2 and CO2 adsorption measurements were maintained at 77.35 K at 101.3 kPa and 273.15 K, respectively. The equilibrium times of the measurements, over which the pressure had to be stable, were short. Equilibrium times of 30 s and 45 s were set for N2 and CO2, respectively. The N2 and CO2 adsorption experiment relative pressures ranged from 0.009 to 0.995 and 0.0001 to 0.032, respectively. A density functional theory (DFT) molecular model was used to determine the pore size distributions (PSDs) of the micropores and mesopores since this technique can provide a more accurate approach for pore size analysis. The N2 adsorption data were interpreted by applying Langmuir analyses and the multi-point Barrett-Emmett-Teller (BET) model for the SA and the Barrett-Joyner-Halenda (BJH) model for the PV of the mesopores. The SA and PV characteristics of micropores measured by CO2 adsorption were analysed using the BET, Langmuir, and DFT models, and detailed descriptions of these theories were previously reported in various studies [9,33,34].

3.5. Determination of Total Porosity

The porosity of the shale samples was calculated as per the differences between bulk density (ρbulk) and skeletal density (ρskeletal); detailed descriptions of these experimental procedures have already been reported in various literature works [9,15,23,24,33]. ρbulk of the cylindrical plug samples were determined by caliper measurements. The drilled plug samples with diameter 2.54 cm and height 4 cm were dried at 110 °C overnight to remove free water. The samples were then weighed in air before and after coated by paraffin of known density, then the paraffin-coated samples were weighed both in air and water of known density to obtain the samples’ bulk volume, finally the ρbulk was calculated by weight in air and bulk volume. After the analysis, the paraffin was crushed between 20 and 40 mesh sizes (380–830 μm) and dried at 110 °C in a vacuum for 24 h. ρskeletal was measured with a helium pycnometer at pressure less than 17.4 psi according to Boyle’s Law. Then the total porosity (φ) of the Xiamaling shales were calculated by the differences between the ρskeletal and ρbulk, according to the following equation:
φ   =   ( 1     ρ bulk ρ skeletal )   ×   100 %

4. Results

4.1. Basic Organic Geochemistry and Mineralogy

The TOC content of the Xiamaling shale varies from 0.55% to 7.60%, with an average value of 3.10%, and the equivalent vitrinite reflectance ranges between 0.52% and 1.88%, indicating that the shale is a “good” source rock and is in the low- and high-maturation stages. Six samples were extracted by dicholoromethane and methanol, and the extracted OM content ranges between 0.79% and 1.93%, with an average value of 1.16% (Table 1).
The Xiamaling shale samples are dominated by clay minerals (21.5–53.8%, average 40.57%) and quartz (21.8–48.4%, average 34.26%), followed by calcite (5.2–32.4%, average 16.59%) and feldspar (1.7–11.4%, average 5.76%). A low content of pyrite (0–7.2%, average 2.82%) is also present. The clay minerals mainly consist of I/S (8.5–40.0%, average 19.86%), followed by illite (3.8–27.4%, average 9.75%) and smectite (3.6–17.9%, average 8.84%), and the content of chlorite is relatively low, with average values less than 3.0% (Table 1).

4.2. Pore Characterization Based on MIP

The macro-PV of the Xiamaling shale samples ranges from 0.0016 cm3/g, as noted in sample JMS-4, to 0.0136, as noted sample XHY-4. The macropore SA ranges from 0.21 m2/g, as noted in sample JZT-6, to 2.67 m2/g, as noted in sample JZT-1. The porosity of the Xiamaling shale samples is in the range of 0.57–3.22% (Table 2).
The different hysteresis loops of MIP curves reflect different pore characteristics and serve as the basis for revealing the pore characteristics in the samples [35,36]. The MIP curves and reflected pores of the Xiamaling shale samples can be divided into three types according to the morphology of the hysteresis loops (Figure 2).
The first type is represented by XHY-3, JMS-4, JMS-5, and XHY-7 and is predominantly composed of pores with widths less than 100 nm, which implies that the mesopores and macropores are more developed. The mercury-in curve increases sharply at approximately 5 MPa, and the mercury-out curve exhibits a convex or a horizontal shape. The hysteresis loop of the MIP curve is wide, with a large difference between the mercury-in volume and the mercury-out volume. This finding implies that most of the pores are open, with favourable pore connectivity. There are large numbers of connected pores with small pore widths; however, macropores do not develop, which is beneficial to the desorption, diffusion, and transfusion of shale gas. The second type is represented by XJG-9, ZJS-6, XJG-3, and ZJS-1 and is characterized by full-scale pores. The hysteresis loop of the MIP curve is wide, with a certain amount of mercury-out as well as a large difference between the mercury-in volume and mercury-out volume. This characteristic indicates that most of the pores are open, with favourable connectivity. Both the macropores and mesopores are developed in the shale, which is also beneficial to the desorption, diffusion, and transfusion of shale gas. The third type includes HTG-3, JZT-1, HTG-2, and HTG-5 and is characterized by pores with diameters greater than 10,000 nm and less than 10 nm. During the initial stage, the mercury-in pressure is low, and the mercury-in speed is relative high. There is a flat stage in the range of 0.1–10 MPa, and the mercury-in amount and the mercury-out speed increase after 10 MPa. The hysteresis loop of the MIP curve is narrow, with a certain mercury-out volume and a small difference between the mercury-in volume and mercury-out volume. This characteristic indicates that there are relatively small numbers of open pores and large numbers of disconnected pores, which are caused by undeveloped pores with diameters of 10–10,000 nm. This type of pore is detrimental to the desorption, diffusion, and transfusion of shale gas.
The main pore structure parameters measured by the low-pressure gas adsorption of the Xiamaling shale are summarized in Table 3. The PV of micropores and mesopores are in the range of 0.0021–0.0245 cm3/g and 0.0046–0.0222 cm3/g, respectively. The SAs of micropores and mesopores are in the ranges of 7.3622–41.2702 m2/g and 4.2128–19.3278 m2/g, respectively.
In this study, the morphological characteristics from the low-pressure N2 adsorption curves of shale samples from the Xiamaling Formation were statically analysed. According to the IUPAC classification, the samples were generally classified into four types (Figure 3). The first type is represented by HTG-2, HTG-3, and HTG-5. The curves exhibit a certain amount of N2 adsorption and a wide hysteresis loop, which corresponds to type H2 isotherms as defined by IUPAC [37]. The second type is represented by JZT-1, ZJS-6, and ZJS-1. The curve shape is similar to that of the first type of curve, but the curves exhibit a large adsorption capacity at P/P0 of approximately 0.995, demonstrating an unsaturated adsorption state and a large hysteresis loop, which is similar to type H2 but has the characteristics of type H3 as defined by IUPAC [37]. The results indicate that these samples have larger pores than the first type, including wedge or V-type pores with both ends open. The third type is represented by XHY-3, JMS-5, and XJG-9. These curves exhibit a certain adsorption content at the low-pressure and high-pressure stage and are steeper at saturated vapour pressure. The hysteresis is narrower than those of the aforementioned two types and close to that of type H3, as defined by IUPAC [37]. The results show that the pores are mainly flat slit structures, cracks, and wedge structures, and the hysteresis loop is mainly related to clay and other flaky granular materials. The fourth type is represented by XJG-3, XHY-7, and JMS-4. These curves are less to the third type of curve, except the increment at the saturated vapor pressure is smaller than that of the third type. These curves are similar to the H3-type curve defined by IUPAC but exhibit the characteristics of the H4-type curve [37] The results show that the pores are slit-shaped pores and wedge-shaped semi-closed pores with parallel walls, and a certain amount of ink bottle-shaped pores are noted.
The PSDs of the Xiamaling shales calculated by the BJH model are presented in Figure 4. The plot of dV/d(logD) versus D clearly reveals that the PV is mainly contributed by the mesopores and micropores and can be divided into three types. The first type is represented by JMS-5, XHY-7, JMS-4, and XHY-3. The TOC of the samples is less than 1.0%, and the clay mineral content is high. The PV distribution is post-peak type, indicating that fewer micropores were generated and that the mesopores were more developed. This finding is attributed to the fact that the samples with low TOC content did not contain the numerous micropores generated by hydrocarbon generation of OM. Clay minerals provide a small number of micropores and a large number of mesopores. The second type is represented by ZJS-1, ZJS-6, XJG-3, and XJG-9. The TOC content varies from 1.0% to 3.0%, and there is a low content of clay minerals. The PV exhibits a bimodal distribution. The main maximum value is distributed in the range of 1–10 nm diameter, and the other is distributed at 2–10 nm diameter. The results showed that the micropores and narrow mesopores in these samples were more developed and contributed the most to the PV. The third type is represented by JZT-1, HTG-2, HTG-3, and HTG-5. The TOC content is greater than 3.0%, and the clay mineral content is much lower than that of the second type. The PV exhibits a pre-peak distribution, indicating that the micropores of these samples are more developed and constitute most of the PV.
The log differential SAs of all the studied samples exhibit broad porosity ranges and approximately similar variation tendencies of bimodal distributions with modes of approximately 1.5 nm and 2–3 nm (Figure 4d–f). The peak values of the log differential SA obviously increase with the TOC content. The pores with diameters <10 nm contribute relatively significantly to the total SA.
The CO2 adsorption isotherms of the Xiamaling shale exhibit similar trends and variations (Figure 5), which manifest as type I according to the IUPAC and are characterized by microporous solids [37]. Obviously as TOC content increases, the maximum volume of the CO2 adsorption content increased from 0.83 to 2.96 cm3/g, indicating that the number of micropores increased with increasing TOC content.
The micropore size distribution of the Xiamaling shale obtained by the DFT model exhibits a high log differential PV and SA and similar distribution characteristics (Figure 6). There are three stable peaks at approximately 0.3~0.4 nm, 0.45–0.6 nm and 0.7~1.0 nm, which exhibit the same increasing and decreasing tendencies. The peak values of PV and SA generally increase with TOC content, which indicates a high abundance of micropore distributions of the corresponding pore size, and the abundance of micropores increases with increasing TOC content. Furthermore, micropores sized 0.45–0.7 nm are the most prevalent in the distribution, and pores with diameters of approximately 0.3~0.45 nm and 0.7~1.0 nm also significantly contribute to the PV and SA, respectively.

5. Discussion

5.1. Relationship between TOC and Quartz and Clay Minerals

The TOC content of 18 shale samples exhibits a weak correlation with the quartz content (R2 = 0.32, Figure 7a), which is also observed in Longmaxi shale from the Sichuan Basin in South China and Devonian shale samples from the Horn River Basin [23,38]. The positive relationship indicates that the quartz is of biogenic origin. Additionally, the TOC content exhibits a negative correlation with clay minerals (R2 = 0.66, Figure 7b). Pan et al. [39] reported that a negative correlation is observed for samples with TOC content < 12.0%, whereas a positive correlation exists for the samples with TOC content > 12.0%. Accordingly, more samples with higher TOC contents are needed to completely estimate the relationship between the TOC content and clay minerals in the Xiamaling shales. In addition, the clay minerals of the Xiamaling shale exhibited marked negative relationships with calcite content (R2 = 0.51, Figure 7c) and quartz content (R2 = 0.42, Figure 7d).

5.2. Full-Scale Pore Structure Characteristics

As reported in previous studies, pore throats are the same as pore bodies, and the PSD curves of MIP and low-pressure N2/CO2 adsorption exhibit good correlations [6,14,24,40]. Therefore, we connect the incremental PV curves and neglect the intrinsic characteristics of the MIP and low-pressure N2/CO2 adsorption techniques. According to the combined representation results, the pore diameters of the shale mainly include the following three types:
Type I: multimodal with multiple scales of pore types, which are represented by XHY-3, XHY-7, JMS-4, and JMS-5 (Figure 8a). This type includes pores with a multi-peak diameter distribution, in addition to micropores and mesopores, and macropores contribute a significant proportion. Micropores, mesopores, and macropores are developed in these types of pores at equivalent scales, which is beneficial to shale gas reservation, migration, and exploitation. Type II: bimodal with primarily mesopores, which are represented by HTG-2, HTG-3, HTG-5, and JMS-1 (Figure 8b). This type includes pores with bimodal distributions of micropores and mesopores, which can provide large amounts of reservation space for shale gas in adsorbed and free states. Type III: bimodal with primarily micropores, which are represented by ZJS-1, ZJS-6, XJG-3, and XJG-9 (Figure 8c). This type includes pores with bimodal distributions of mesopores and micropores, which can provide large amounts of SA and are extremely beneficial for reserving shale gas in the adsorbed state.
In addition, it is worth noting that the pore classification reflected by the MIP curve and the low-pressure gas adsorption and desorption curve in this study is also applicable to shale reservoirs of different ages and strata in other regions, and the full-scale pore classification is also generally and comprehensively representative and applicable. Compared with the previous classification, these classification and pore diameter distribution features can be used to evaluate the development degree, morphological characteristics, and connectivity of pores with different pore sizes in shale more comprehensively, intuitively, and deeply, and reveal their influence on shale gas desorption, diffusion, and seepage more clearly [14,24,35,40]. Furthermore, these classifications can clearly demonstrate the contribution of pores with different sizes to PV and SA, and distinguish their influences on the occurrence, migration, and later exploitation of shale gas with different occurrence states.

5.3. Control Effects on Pore Structure

5.3.1. Porosity and Controlling Factors

Figure 9a demonstrates that the total porosity is weakly positively correlated with the TOC content, indicating that OM plays a significant positive role in the total porosity but is not the most essential controlling factor. The positive relationship between the total porosity and TOC content was previously mentioned by many scholars [23,24,41], and some studies have reported that the total porosity is positively correlated with the TOC content but is negatively correlated with increasing TOC content for samples with TOC content higher than a certain value owing to diagenesis and compaction [19,39].
As shown in Figure 9b, the porosity is positively correlated with quartz. The biogenic quartz in marine shale may contain a large number of intraparticle pores [1,2], but this is unlikely to be the main reason for the elevated porosity of the Xiamaling shale since the TOC content increases with increasing quartz content (Figure 7a). The relationship between the total porosity and quartz content actually reflects the correlation of the porosity with the TOC content. The total porosity is negatively correlated with clay content according to the illustration of Figure 9c, indicating that the contribution of clay minerals to the total porosity of the Ximaling shale could be negligible, even though clay minerals contribute to the porosity of the shale [2,31]. The interpretation of these relationships needs to be confirmed with more detailed research, given that this is a preliminary result.

5.3.2. Effect of TOC Content on Pore Structure

This research investigated the relationship between the TOC content and nanoscale pore structure parameters in Xiamaling shale (Figure 10). A positive linear correlation is observed between the TOC content and PV of the micropores, mesopores, and macropores (R2 = 0.41, 0.40 and 0.62, respectively) (Figure 10a). Strong positive relationships are also found between the TOC contents and micropore and mesopore SAs (R2 = 0.85, 0.56, respectively). However, the macropore SA remains stable with variation of the TOC content. The aforementioned appearance suggests that OM is the most important contributor to micro- and mesopores and that OM can be used as an important parameter to evaluate the micropores and mesopores in Xiamaling shale. The shale with higher TOC content tends to have more micropores and mesopores with larger PV and SA. The thermal maturation and hydrocarbon expulsion of OM results in the generation and development of abundant OM pores, which are mainly micro- and mesopores.

5.3.3. Effect of Mineral Composition on Pore Structure

As shown in Figure 11a,b, a relatively weak negative correlation is observed between the clay minerals and PV of the micropores, mesopores, and macropores (R2 = 0.14, 0.12 and 0.13, respectively). The SA of the micropores and mesopores is negatively correlated with the clay minerals (R2 = 0.55, 0.27), whereas the macro-PV exhibits no correlation with the clay minerals. The negative correlations between the clay mineral content and the PV and SA may be caused by the decrease in pore space that occurs because various pores are filled as the clay mineral content increases.
In Figure 11c,d, we observe that almost all of the relationships between the quartz and the PV of the micropores, mesopores, and macropores are relatively weakly positive (R2 = 0.07, 0.08 and 0.05, respectively). The relationships between the quartz and the SA of micropores and mesopores are weakly positive (R2 = 0.18, 0.17, respectively). The SA of the macropores remains stable as the quartz content changes. The positive relationships between the quartz and the SA and PV of the nanoscale pores are similar to those of the Wufeng–Longmaxi shale from the Jiashiba area, Sichuan Basin of China and the Devonian gas shales in the Horn River Basin of Canada [15,24,38]. However, Yang et al. [42] reported that the Lower Permian Shanxi shale in the Ordos Basin exhibits the opposite tendency. The different relationships between quartz and the pore structure parameters of the PV and SA in different study areas may be due to the diagenetic or sedimentary environment and the origin of quartz [43]. According to the aforementioned discussions, we know that the TOC content is the main factor controlling nanoscale pore structures. To thoroughly analyse the influence of clay minerals and quartz on the nanoscale pore structure, the PV and SA are normalized by the TOC content to remove the obscure TOC content influencing the relationship between the pore structure and clay minerals, calcites, and quartz content. The analysis results show that the clay mineral content exhibits strong positive relationships with the TOC-normalized PVs of the micropores, mesopores, and macropores (R2 = 0.63, 0.72 and 0.8, respectively) (Figure 12a). The TOC-normalized PVs of the micropores, mesopores, and macropores exhibit weak positive relationships with the smectite (R2 = 0.12, 0.04 and 0.08, respectively) and I/S contents (R2 are 0.27, 0.54 and 0.54, respectively) and a relatively strong positive correlation with illite (R2 = 0.5, 0.48 and 0.3, respectively) (Figure 13a–c).
The clay minerals are relatively weakly positively correlated with the TOC-normalized SA of micro- and mesopores (R2 = 0.28 and 0.28, respectively), whereas they exhibit no obvious relationship with the TOC-normalized SA of macropores (Figure 12b). Further analysis reveals that the TOC-normalized SA of the micropores, mesopores, and macropores exhibits an obvious positive correlation with smectite (R2 = 0.03, 0.03 and 0.06, respectively) and a good correlation with I/S (R2 = 0.28, 0.34 and 0.30, respectively) but a weak negative correlation with illite (R2 = 0.39, 0.29 and 0.41, respectively) (Figure 13d–f). This relationship indicates that I/S and illite are the major contributors to the PV and SA, whereas smectite makes a negligible contribution to the pore structure. The microporosity of clays in shale is related to the clay crystal size and irregular pores between clay frameworks; these results were confirmed by previous studies [1,15,24,44]. In addition, a positive correlation exists between the clay mineral content and the PV and SA, as shown in Figure 12a,b and Figure 13 but not the negative correlation shown in Figure 11a,b. This finding is attributed to the fact that the TOC content obscures the influences of the clay minerals on the PV and SA of the micropores, mesopores, and macropores. The clay minerals containing smectite and I/S in clay-rich shale with low TOC content make a significant contribution to the pore structure in organic-lean shales.
Furthermore, the TOC-normalized PV and SA of the micropores, mesopores, and macropores are negatively correlated with calcite (Figure 12c,d). However, Yang et al. [45] reported that theoretically, the composition cannot have a negative effect on the PV or SA. In reality, the passively negative correlations between calcite and pores of different diameter scales originated from calcite’s negative relationship with clay (Figure 7c).
Additionally, an extremely weak negative correlation is observed between the quartz content and the TOC-normalized PV and SA of the micropores, mesopores, and macropores (Figure 12e,f). Loucks et al. [7] found that the intragranular pores associated with quartz are less abundant in the siliceous mudstones of the Mississippian Barnett shale and that the mechanical compaction and pore blocking of secondary OM and clay floccus results in the destruction of these pores. Therefore, theoretically, quartz makes a minimal contribution to the pore structures of different diameter scales in shale.
It should be noted that although many scholars have analysed the influence of mineral composition on pore structure, most of the research results have not removed the concealment of TOC content on pore structure, and the correlation between each clay mineral composition and pore structure of different scales is rarely discussed under the removal of TOC concealment [23,29,46]. After removing the concealing effect of TOC content on the micro-nano pore structure, this study comprehensively and deeply reveals the influence of brittle minerals such as quartz and calcite, clay minerals such as smectite, I/S and illite, on the formation and development of the PV and SA of micropores, mesopores, and macropore. The influence of the PV and SA of macropores on the formation and development of shale pores is carefully discussed. The control mechanism of mineral composition on pore structure formation and development of shale is discussed in detail. These conclusions are of great significance in revealing the formation and development mechanism of micro-nano pores in Marine shale.

5.3.4. Effect of Maturation on Pore Evolution

To obtain information about the effect of thermal maturity on the nanoscale pore structure, this research estimated the relationship between the TOC-normalized PV and SA of different nanoscale pores of Ximaling shale (Figure 14). The results show that the TOC-normalized PV of micropores, mesopores, and macropores are negatively correlated with EqRo (R2 = 0.21, 0.28 and 0.32, respectively). The TOC-normalized micropores and mesopore SAs have a relatively weak negative relationship with EqRo (R2 = 0.16 and 0.24, respectively), whereas EqRo has a weaker negative effect on the TOC-normalized macropore SA (R2 = 0.05). The thermal maturity reflects the dual effects of compaction and time experienced by shale. According to Liu et al. [29], the Xiamaling shale in Zhangjiakou, Hebei has been in a slow subsidence and deep burying state for a long time since its deposition to the end of the Carboniferous; then, this shale started lifting at the beginning of the Permian and subsided rapidly in the middle of the Cretaceous. The negative relationship between EqRo and the TOC-normalized PV and SA of different nanoscale pores of the Ximaling shale is predominantly due to the deep progressive burial and compaction of the shale.
Various studies have reported on nanoscale pore evolution with progressing thermal maturity in outcrop shale samples [19,34,47,48,49]. Topór et al. [47] proposed a pore evolution model based on the Silurian Baltic Basin shales with maturation sequences varying from 0.4% to 3.5%. However, the evolution model is ambiguous for Ro > 2.0% given the lack of research dates. Pommer and Milliken [48] established a pore evolution pattern for the Eagle Ford shales, with Ro values ranging between 0.5 and 1.3% but did not interpret the evolution characteristics for shales with maturity greater than 2.0%. To interpret the pore evolution trend with thermal maturity, the TOC-normalized PVs of different diameter scales with progressing EqRo are established based on the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei (Figure 15).The increases in the TOC-normalized PV of different diameter scales from the immature to early stage (0.4% < EqRo < 0.7%) mainly result from the progressing burial and compaction, which results in the adjustment of the pore size proportions, especially the conversion of macropores to mesopores and micropores [4,23,50]. Furthermore, the liquid hydrocarbon expulsion of OM and kerogen at a maturity of approximately 0.4–0.7% could also result in an increase in the PV owing to the appearance of organic bubble pores [22]. Additionally, the earlier dissolution and transformation of unstable minerals can contribute to the PV increases at different diameter scales [50,51].
The clear decrease in the TOC-normalized PV of different diameter scales in the middle mature stage with EqRo in the range of 0.7∼1.2% resulted from the oil and bitumen infill in pore spaces, which is also confirmed by the significant increase in the meso- and macro-PVs in the “oil window” after extraction reported by Guo et al. [52]. The obvious decreasing tendency of the TOC-normalized PV of different diameter scales may also result from the adjustment of clay aggregates and continuous compaction.
During the late mature stage (1.2% < EqRo < 2.0%), the significant increase in the TOC-normalized PV of different diameter scales is mainly caused by the secondary cracking of oil and bitumen, accompanied by wet gas generation, which increases the pore space. Hu et al. [53] observed that the porosity increased after Ro 0.9% and reached the maximum value at 1.5∼2.0%. Kuila et al. [54] proposed a significant increase in the PV<5nm for Baltic basin shale with maturity in the range of 1.2∼2.0%. Various similar conclusions can also be found in previous studies [24,34,47,55].
This study found a negative correlation between maturity and shale pore structure, and further analysis found that the TOC-normalized pore volume of micropores, mesopores, and macropores of the shale does not exhibit monotonous increase or decrease. These pores increased first, and then increased again after decrease as the maturity increased to 2.0%. Thus, there is a larger value and a minimum value at 0.7% and 1.2%, respectively. The research conclusion of the influence factors of maturity on pore structure in this study complements the previous research results of shale in high and over matured area, and especially has general applicability to other immature to highly matured marine shale formations of different ages.

5.4. Effect of Residual Hydrocarbon on the Pore Structure

5.4.1. Adsorption Desorption Curve Characteristics

In this research, low-pressure CO2 and N2 adsorption were conducted on six original Xiamaling shales and the samples extracted using a 25:2 vol/vol mixture of dicholoromethane and methanol (Figure 16 and Figure 17). The detected samples are of similar OM maturity, type, and mineral composition and content but with different TOC. The result show that CO2 and N2 adsorption contents of the extracted shales are greater than those of the original shale, indicating that the residual hydrocarbon fills in the pores or blocks the pore throats. When the residual hydrocarbon in the samples is extracted, the pore space is released, and the gas adsorption capacity of the sample increases. Valenza II et al. [56] researched the pore characteristics of shale in the oil generation stage and found that the SA and PV of the extracted samples increased significantly. Jarvia et al. [20] also found that residual bitumen blocked the pores and throats of shale and restricted the adsorption and migration of shale gas. With increasing maturity, the residual bitumen gradually became solid, and its influence on the nanoscale pore structure decreased. Various scholars have found that residual hydrocarbon can block the pores of shale, resulting in a decrease in the SA and PV, supporting our research results [4,34,57].
As shown in Figure 18a,b we found that the TOC content is positively correlated with CO2/N2 adsorption content. The higher the content of TOC, the more micropores and mesopores are enriched in OM at the low maturity stage. In addition, a positive correlation exists between the TOC content and residual hydrocarbon content (R2 = 0.84, Figure 18c), and the difference between the adsorption content of the extracted sample and original shale is positively correlated with the TOC content (R2 = 0.54 of CO2 and 0.17 of N2, Figure 18d). These results imply that under the same conditions for OM type, maturity, and mineral components, a higher TOC content indicates a larger amount of residual hydrocarbon generated by shale, with a greater occupied pore space and greater adsorbing capacity after extraction. Furthermore, the negative correlation between the residual hydrocarbon and TOC-normalized adsorbing capacity of CO2 and N2 (R2 = 0.42, 0.07, Figure 18d,e) indicate that a higher content of residual hydrocarbon means a larger extent of blocking of the nanoscale pores of shale. These research results have already been presented in various studies [4,57].
As listed in Table 4, gap values of the micro- and meso-PVs of the extracted shale and original sample are in the ranges of 0.001–0.011 cm3/g and 0.001–0.006 cm3/g, respectively. Gap values of the micropore and mesopore SAs are in the ranges of 2.309–4.676 m2/g and 1.579–4.992 m2/g, respectively. The results imply that the pore structures of micropores and mesopores in the extracted samples are significantly larger than those of the original samples. These results indicate that the presence of residual hydrocarbon would occupy the pore space of micropores and mesopores and that the extraction could increase the PV and SA of micropore and mesopores, dredge the micropores and mesopores in shale, and improve the storage capacity and transfusion of shale gas.

5.4.2. Variation Characteristics of the Micropore Structures

According to the variation characteristics of the micro-PVs in Table 2 and Figure 19, the micro-PVs of the various pore ranges of the extracted samples increase, and the three-peak distribution characteristics of the PV are constant. In addition, the pores with widths of 0.3~0.4 nm, 0.45–0.6 nm, and 0.7~1.0 nm are the main contributors to the micro-PV, and the rates of change in the PVs of the 3 main peaks are relatively high. These results indicate that the pores with widths of 0.3–1.0 nm exhibit a relatively large PV and provide most of the PV of the micropores, whereas the rate of change in the PV of the pores with widths greater than 1.0 nm is 0, indicating that the micropores with diameters larger than 1 nm make almost no contribution to the PV.
The SA of the micropores of the six original shales and the extracted samples are distributed in a three-peak manner, and the main peak values are situated at approximately 0.3–0.4 nm, 0.45–0.6 nm, and 0.7–1.0 nm (Figure 20). As the pore diameter increases gradually from 0.3 nm to approximately 1.0 nm, the micropore SA experiences a relatively high rate of change, indicating that the pores of the three peaks have a relatively large SA, contributing most of the SA of the micropores. The rate of change in the SA of the pores with widths greater than 1 nm in various samples is near 0, indicating that the micropores with diameters larger than 1 nm contribute minimally to the SA.
The micropore SA and PV of the 6 extracted shales are greater than those of the original samples. The six samples have different TOC contents and residual hydrocarbon contents, leading to large changes in the micropore SA and PV (Figure 21), indicating that the contents of TOC and residual hydrocarbon greatly influence the micropore structure of shale. The micropore SA and PV of the six samples have similar variation characteristics to that of the pore diameter distribution. The micropore SA and PV of the original shale and the corresponding extracted samples have the same distribution characteristics. The micropore SA and PV are larger at 0.3–1.0 nm, and there is a maximum peak value in the range of 0.45–0.6 nm. However, these parameters are relatively smaller for pore diameters greater than 1 nm and less than 0.3 nm. This finding indicates that the residual hydrocarbon mainly fills micropores with widths of 0.3–1.0 nm, resulting in a significant difference between the micropore SAs and PVs for the original shale and the extracted samples with widths of 0.3–1.0 nm.
As shown in Figure 22, the extractable OM (EOM) is significantly positively correlated with the gap value of the micro-PV (R2 = 0.48) and SA (R2 = 0.69). This finding indicates that residual hydrocarbon is one of the main controlling factors for the formation and development of the micropore PV and SA of Xiamaling shale; a higher residual hydrocarbon content occupies a larger micropore SA and PV and results in a larger increase in the micropore SA and PV in the extracted samples.

5.4.3. Variation Characteristics of the Mesopore and Macropore Structures

As shown in Figure 23, the meso- and macro-PV of the extracted samples are significantly greater than those of the original shales. ZJS-3, JMS-9, XHY-4 and their extracted samples are of the pre-peak type (Figure 23a,e,f). The pores with smaller diameters contribute relatively more to the PV, and the pores with diameters less than 10 nm contribute the most to the PV. The micropores and mesopores with small diameters constitute the majority of the PV. The PV distributions of JZT-6, XJG-6, and HTG-4 are bi-peak with 2 maximum values at 2 nm and 40 nm and are dominated by the latter (Figure 23b–d).
The SA of the six extracted samples is greater than that of the original shale. The SA of the mesopores and macropores exhibits a unimodal distribution, and the pores with smaller pore diameters contribute more to the SA. Pores with widths less than 10 nm contribute the most to the SA (Figure 24). The SA variation characteristics of the mesopores and macropores are similar to that of the PVs, with some differences, dominated by the following three aspects.
First, the PV distribution of the mesopores and macropores exhibits pre-peak and bi-peak characteristics, whereas the SA distribution shows only pre-peak characteristics. The PV of the pre-peak pores is mainly contributed by pores with widths less than 10 nm, whereas that of the bi-peak pores is mainly contributed by pores with diameters of 2 nm and 40 nm. The SA of the mesopores is mainly contributed by pores with widths less than 10 nm. Second, the mesopore structure variation characteristics of the original shale are similar to those of the extracted samples, and the SAs of the extracted samples are larger than those of the original shale. This finding is attributed to the fact that the pore space is released and increased due to the extraction of the solid residual hydrocarbon in shale, such as solid asphalt. Third, residual hydrocarbon may gather in pores with widths less than 10 nm, as well as those sized 2 nm and 40 nm, so there are significant changes between the original shale and extracted samples regarding the PV and SA of pores with diameters of 2 nm and 40 nm.
As shown in Figure 25, the mesopore SA and PV of the six extracted samples are greater than that of the original shale when they were extracted with Soxhlet extraction. There are great differences between the changes in mesopore SA and PV given the different contents of TOC and residual hydrocarbon in the six samples, indicating that the mesopore structure of the shale is also controlled by the TOC and residual hydrocarbon contents. The mesopore SA and PV distribution characteristics of the six original shale samples are the same as those of the extracted samples, and pores with diameters less than 10 nm, 2 nm, and 40 nm may be where residual hydrocarbon is stored. These characteristics result in significant difference between the SA and PV of the pre-peak type mesopores with widths less than 10 nm of the original shale and those of the extracted samples, as well as that of the PV of the bi-peak type mesopores with widths of 2 nm and 40 nm.
The EOM and gap value of meso-PV and SA exhibits an obvious positive correlation with coefficients of 0.63 and 0.63, respectively (Figure 26a,b). This finding implied that residual hydrocarbon have a significant effect on the formation and development of the mesopore SA and PV of Xiamaling shale. The higher the content of residual hydrocarbon, the larger the mesopore SA and PV it occupies, and the larger the SA and PV of the mesopore of the sample after extraction, which is consistent with the correlation between residual hydrocarbon and micropore structure.
It is believed that residual hydrocarbons can block pores and throats of shale, and have a certain restrictive effect on gas adsorption and migration. However, few scholars have reported that residual hydrocarbons specifically affect the pore of that pore size, and reports about the correlation between residual hydrocarbon content and SA and PV of different scales are scarce [4,20,34,56,58,59]. This study clearly shows that residual hydrocarbons mainly affect the changes in the pore structure parameters, reveals the degree of restriction on the SA and PV of pores with different diameter, clarifies positive correlation between residual hydrocarbon content and SA and PV of pores with different scales, and reveals the main storage space of residual hydrocarbon in the reservoir and its influence on reservoir porosity characterization. Residual hydrocarbons have an important impact on the gas-bearing properties of matured shale. These studies provide important evidence for understanding and evaluating the adsorption, desorption, diffusion, and migration of shale gas.

6. Conclusions

(1) MIP curve of the Xiamaling shale can be divided into three types, which represent three different sets of pore characteristics and reflect different pore connectivity performances. The N2 isotherm curve can be roughly classified into four types, denoting four different types of mesopores. The N2 PV distribution exhibits a pre-peaks, bi-peaks, and post-peaks due to the influences of the TOC content and mineral components, whereas the micropore SA and PV exhibits a three-peak type distribution. The full-scale distribution characteristics of the samples are multimodal with multiple scales, bi-peaks with primarily mesopores or micropores.
(2) OM contributes mostly to the development of micro- and mesopores. Clay minerals, mainly smectite and I/S, are the main contributors to the TOC-normalized PV and SA of the micropores, mesopores, and macropores. Illite is positively correlated with the TOC-normalized PV, but is negatively correlated with the TOC-normalized SA of pores with different diameter scales.
(3) The thermal maturity controls the generation and development of pores at different diameter scales. Additionally, this research establishes a general TOC-normalized PV evolution mode.
(4) The micropore and mesopore parameters are significant affected by TOC and residual hydrocarbon content. A increased residual hydrocarbon content indicates a more occupied pore space. The residual hydrocarbon mainly fills micropores with diameters of 0.3–1.0 nm and mesopores with widths less than 10 nm and those of 2 nm and 40 nm. The SA and PV of the extracted sample increase significantly due to the release of the occupied pore space.

Author Contributions

Conceptualization, L.X., K.Y. and L.L.; methodology, L.X., K.Y., H.W. and L.L.; software, X.L., T.X. and X.W.; validation, K.Y. and H.W.; data curation, L.C.; writing—original draft preparation, L.X., K.Y. and L.L.; writing—review and editing, L.X., K.Y. and H.W.; supervision, L.L.; project administration, K.Y. and L.L.; funding acquisition, K.Y. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

Please add: This research was funded by Chinese National Science and Technology Major Project (No. 2017ZX05035-002); Key Research and Development Project of Hebei Province (No. 19224205D); Natural Science Foundation of Hebei Province (No. D2020403019).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ross, D.J.K.; Bustin, R.M. The importance of shale composition and pore structure upon gas storage potential of shale gas reservoirs. Mar. Pet. Geol. 2009, 26, 916–927. [Google Scholar] [CrossRef]
  2. Chalmers, G.R.L.; Bustin, R.M.; Power, I.M. Characterization of gas shale pore systems by porosimetry, pycnometry, surface area, and field emission scanning electron microscopy/transmission electron microscopy image analyses: Examples from the Barnett, Woodford, Haynesville, Marcellus, and Doig units. AAPG Bull. 2012, 96, 1099–1119. [Google Scholar]
  3. Li, J.; Zhou, S.X.; Li, Y.J.; Yang, Y.N.; Li, C.C. Effect of organic matter on pore structure of mature lacustrine organic-rich shale: A case study of the Triassic Yanchang shale, Ordos Basin, China. Fuel 2016, 185, 421–431. [Google Scholar] [CrossRef]
  4. Guo, H.J.; Jia, W.L.; Peng, P.A.; Lei, Y.H.; Luo, X.R.; Cheng, M.; Wang, X.Z.; Zhang, L.X.; Jiang, C.F. The composition and its impact on the methane sorption of lacustrine shales from the Upper Triassic Yanchang Formation, Ordos Basin, China. Mar. Pet. Geol. 2014, 57, 509–520. [Google Scholar] [CrossRef]
  5. Curtis, J.B. Fractured shale-gas systems. AAPG Bull. 2002, 11, 1921–1938. [Google Scholar]
  6. Clarkson, C.R.; Solano, N.; Bustin, R.M.; Chalmers, G.R.L.; Melnichenko, Y.B.; Radlinski, A.P. Pore structure characterization of North American shale gas reservoirs using USANS/SANS, gas adsorption, and mercury intrusion. Fuel 2013, 103, 606–616. [Google Scholar] [CrossRef]
  7. Loucks, R.G.; Reed, R.M.; Ruppel, S.C.; Jarvie, D.M. Morphology, genesis, and distribution of nanometer-scale pores in siliceous mudstones of the Mississippian Barnett Shale. J. Sediment. Res. 2009, 79, 848–861. [Google Scholar] [CrossRef] [Green Version]
  8. Bernard, S.; Horsfield, B.; Schulz, H.M.; Wirth, R.; Schreiber, A.; Sherwood, N. Geochemical evolution of organic-rich shales with increasing maturity: A STXM and TEM study of the Posidonia Shale (Lower Toarcian, northern Germany). Mar. Pet. Geol. 2012, 31, 70–89. [Google Scholar] [CrossRef]
  9. Mastalerz, M.; He, L.; Melnichenko, Y.B.; Rupp, J.A. Porosity of coal and shale: Insights from gas adsorption and SANS/USANS techniques. Energy Fuels 2012, 8, 5109–5120. [Google Scholar] [CrossRef]
  10. Klaver, J.; Desbois, G.; Littke, R.; Urai, J.L. BIB-SEM characterization of pore space morphology and distribution in postmature to overmature samples from the Haynesville and Bossier Shales. Mar. Pet. Geol. 2015, 59, 451–466. [Google Scholar] [CrossRef]
  11. Gou, Q.Y.; Xu, S.; Hao, F.; Yang, F.; Zhang, B.Q.; Shu, Z.G.; Zhang, A.H.; Wang, Y.X.; Lu, Y.B.; Cheng, X.; et al. Full-scale pores and micro-fractures characterization using FE-SEM, gas adsorption, nano-CT and micro-CT: A case study of the Silurian Longmaxi Formation shale in the Fuling area, Sichuan Basin, China. Fuel 2019, 253, 167–179. [Google Scholar] [CrossRef]
  12. Wang, P.F.; Jiang, Z.X.; Ji, W.M.; Zhang, C.; Yuan, Y.; Chen, L.; Yin, L.S. Heterogeneity of intergranular, intraparticle and organic pores in Longmaxi shale in Sichuan Basin, South China: Evidence from SEM digital images and fractal and multifractal geometries. Mar. Pet. Geol. 2016, 72, 122–138. [Google Scholar] [CrossRef]
  13. Fu, C.Q.; Zhu, Y.M.; Chen, S.B.; Xiao, H.X. Structure and fractal characteristics of nano-micro pores in organic-rich Qiongzhusi formation shales in Eastern Yunnan Province. J. Nanosci. Nanotechnol. 2017, 17, 5996–6013. [Google Scholar] [CrossRef]
  14. Wang, Y.; Zhu, Y.M.; Chen, S.B.; Li, W. Characteristics of the nanoscale pore structure in Northwestern Hunan Shale gas reservoirs using field emission scanning electron microscopy, high-pressure mercury intrusion, and gas adsorption. Energy Fuels 2014, 28, 945–955. [Google Scholar] [CrossRef]
  15. Yang, F.; Ning, Z.F.; Qiang, W.; Zhang, R.; Krooss, B.M. Pore structure characteristics of lower Silurian shales in the southern Sichuan Basin, China: Insights to pore development and gas storage mechanism. Int. J. Coal Geol. 2016, 156, 12–24. [Google Scholar] [CrossRef]
  16. Ma, Y.; Ardakani, O.H.; Zhong, N.N.; Liu, H.L.; Huang, H.P.; Larter, S.; Zhang, C. Possible pore structure deformation effects on the shale gas enrichment: An example from the Lower Cambrian shales of the Eastern Upper Yangtze Platform, South China. Int. J. Coal Geol. 2020, 217, 1–22. [Google Scholar] [CrossRef]
  17. Hao, F.; Zou, H.Y.; Lu, Y.C. Mechanisms of shale gas storage: Implications for shale gas exploration in China. AAPG Bull. 2013, 97, 1325–1346. [Google Scholar] [CrossRef]
  18. Liu, H.M.; Zhang, S.; Song, G.Q.; Wang, X.J.; Teng, J.B.; Wang, M.; Bao, Y.S.; Yao, S.P.; Wang, W.Q.; Zhang, S.P.; et al. Effect of shale diagenesis on pores and storage capacity in the Paleogene Shahejie Formation, Dongying Depression, Bohai Bay Basin, east China. Mar. Pet. Geol. 2019, 103, 738–752. [Google Scholar] [CrossRef]
  19. Milliken, K.L.; Rudnicki, M.; Awwiller, D.N.; Zhang, T.W. Organic matter-hosted pore system, Marcellus Formation (Devonian), Pennsylvania. AAPG Bull. 2013, 97, 177–200. [Google Scholar] [CrossRef]
  20. Jarvie, D.M.; Hill, R.J.; Ruble, T.E.; Pollastro, R.M. Unconventional shale-gas systems: The Mississippian Barnett Shale of north-central Texas as one model for thermogenic shale-gas assessment. AAPG Bull. 2007, 91, 475–499. [Google Scholar] [CrossRef]
  21. Modica, C.J.; Lapierre, S.G. Estimation of kerogen porosity in source rocks as a function of thermal transformation: Example from the Mowry Shale in the Powder River Basin of Wyoming. AAPG Bull. 2012, 96, 87–108. [Google Scholar] [CrossRef]
  22. Loucks, R.G.; Reed, R.M.; Ruppel, S.C. Spectrum of poretypes and networks in mudrocks and a descriptive classification for matrix-related mudrock pores. AAPG Bull. 2012, 96, 1071–1098. [Google Scholar] [CrossRef] [Green Version]
  23. Chen, L.; Jiang, Z.X.; Liu, K.Y.; Tan, J.Q.; Gao, F.L.; Wang, P.F. Pore structure characterization for organic-rich Lower Silurian shale in the Upper Yangtze Platform, South China: A possible mechanism for pore development. J. Nat. Gas Sci. Eng. 2017, 46, 1–15. [Google Scholar] [CrossRef]
  24. Wang, Y.; Liu, L.F.; Zheng, S.S.; Luo, Z.H.; Sheng, Y.; Wang, X.M. Full-scale pore structure and its controlling factors of the Wufeng-Longmaxi shale, southern Sichuan Basin, China: Implications for pore evolution of highly overmature marine shale. J. Nat. Gas Sci. Eng. 2019, 67, 134–146. [Google Scholar] [CrossRef]
  25. Wang, P.F.; Jiang, Z.X.; Chen, L.; Yin, L.S.; Li, Z.; Zhang, C.; Tang, X.L.; Wang, G.Z. Pore structure characterization for the Longmaxi and Niutitang shales in the Upper Yangtze Platform, South China: Evidence from focused ion beam–He ion microscopy, nano-computerized tomography and gas adsorption analysis. Mar. Pet. Geol. 2016, 72, 1323–1337. [Google Scholar] [CrossRef]
  26. Bian, L.Z.; Zhang, S.C.; Zhang, B.M.; Wang, D.R. Red algal fossils discovered from the Neoproterozoic Xiamaling oil shales, Xiahuayuan town of Hebei province. Acta Micropalaeontol. Sin. 2005, 22, 209–216. [Google Scholar]
  27. Zhang, S.C.; Zhang, B.M.; Bian, L.Z.; Jin, Z.J.; Wang, D.R.; Chen, J.F. The Xiamaling oil shale generated through rhodophyal over 800 Ma ago. Sci. China Ser. D Earth Sci. 2007, 37, 636–643. [Google Scholar] [CrossRef]
  28. Gao, L.Z.; Zhang, C.H.; Shi, X.Y.; Zhou, H.R.; Wang, Z.Q. Zircon SHRIMP U-Pb dating of the tuffbed in the Ximaling formation of the Qinbaikou system in North China. Geol. Bull. China 2007, 26, 249–255. [Google Scholar]
  29. Liu, A.Q.; Tang, D.J.; Shi, X.Y.; Zhou, L.M.; Zhou, X.Q.; Shang, M.H.; Li, Y.; Song, H.Y. Growth mechanisms and environmental implications of carbonate concretions from the ~1.4 Ga Xiamaling Formation, North China. J. Palaeogeogr. 2019, 8, 1–16. [Google Scholar] [CrossRef] [Green Version]
  30. Xiao, X.M.; Wilkins, R.W.T.; Liu, D.H.; Liu, Z.F.; Fu, J.M. Investigation of thermal maturity of lower Palaeozoic hydrocarbon source rocks by means of vitrinite-like maceral reflectance—A Tarim Basin case study. Org. Geochem. 2000, 31, 1041–1052. [Google Scholar]
  31. Chalmers, G.R.L.; Bustin, R.M. Lower Cretaceous gas shales in northeastern British Columbia, Part II: Evaluation of regional potential gas resources. Bull. Can. Pet. Geol. 2008, 56, 22–61. [Google Scholar] [CrossRef]
  32. Washburn, E.W. The dynamics of capillary flow. Phys. Rev. 1921, 17, 273–283. [Google Scholar] [CrossRef]
  33. Tian, H.; Pan, L.; Xiao, X.M.; Wilkins, R.W.T.; Meng, Z.P.; Huang, B.J. A preliminary study on the pore characterization of Lower Silurian black shales in the Chuandong Thrust Fold Belt, southwestern China using low pressure N2 adsorption and FE-SEM methods. Mar. Pet. Geol. 2013, 48, 8–19. [Google Scholar] [CrossRef]
  34. Mastalerz, M.; Schimmelmann, A.; Drobniak, A.; Chen, Y.Y. Porosity of Devonian and Mississippian New Albany Shale across a maturation gradient: Insights from organic petrology, gas adsorption, and mercury intrusion. AAPG Bull. 2013, 97, 1621–1643. [Google Scholar] [CrossRef]
  35. Chen, S.B.; Zhu, Y.M.; Li, W.; Wang, H. Influence of magma intrusion on gas outburst in a low rank coal mine. Int. J. Min. Sci. Technol. 2012, 22, 259–266. [Google Scholar] [CrossRef]
  36. Schmitt, M.; Fernandes, C.P.; Da Cunha Neto, J.A.B.; Wolf, F.G.; Santos, V.D. Characterization of pore systems in seal rocks using Nitrogen Gas Adsorption combined with Mercury Injection Capillary Pressure techniques. Mar. Pet. Geol. 2013, 39, 138–149. [Google Scholar] [CrossRef]
  37. Sing, K.; Everett, D.; Haul, R.; Moscou, L.; Pierotti, R.A.; Rouquerol, J.; Siemieniewsha, T. Reporting physisorption data for gas/solid systems with special reference to the determination of surface area and porosity. Pure Appl. Chem. 1985, 11, 603–619. [Google Scholar] [CrossRef]
  38. Chalmers, G.R.L.; Ross, D.J.K.; Bustin, R.M. Geological controls on matrix permeability of Devonian Gas Shales in the Horn River and Liard basins, northeastern British Columbia, Canada. Int. J. Coal Geol. 2012, 103, 120–131. [Google Scholar] [CrossRef]
  39. Pan, L.; Xiao, X.M.; Tian, H.; Zhou, Q.; Chen, J.; Li, T.F.; Wei, Q. A preliminary study on the characterization and controlling factors of porosity and pore structure of the Permian shales in Lower Yangtze region, Eastern China. Int. J. Coal Geol. 2015, 146, 68–78. [Google Scholar] [CrossRef]
  40. Wang, Y.; Zhu, Y.M.; Wang, H.Y.; Feng, G.J. Nanoscale pore morphology and distribution of lacustrine shale reservoirs: Examples from the Upper Triassic Yanchang Formation, Ordos Basin. J. Energy Chem. 2015, 24, 512–519. [Google Scholar] [CrossRef]
  41. Tian, H.; Pan, L.; Zhang, T.W.; Xiao, X.M.; Meng, Z.P.; Huang, B.J. Pore characterization of organic-rich Lower Cambrian shales in Qiannan Depression of Guizhou Province, Southwestern China. Mar. Pet. Geol. 2015, 62, 28–43. [Google Scholar] [CrossRef]
  42. Yang, C.; Zhang, J.C.; Tang, X.; Ding, J.H.; Zhao, Q.R.; Dang, W.; Chen, H.Y.; Su, Y.; Li, B.W.; Lu, D.F. Nanoscale pore structure and fractal characteristics of marine-continental transitional shale: A case study from the Lower Permian Shanxi Shale in the southeastern Ordos Basin, China. Mar. Pet. Geol. 2017, 88, 54–68. [Google Scholar] [CrossRef]
  43. Zeng, W.T.; Zhang, J.C.; Ding, W.L.; Wang, X.Z.; Zhu, D.W.; Liu, Z.J. The gas content of continental Yanchang shale and its main controlling factores: A case study of Liuping-171 well in Ordos Basin. Nat. Gas Geosci. 2014, 25, 291–301. [Google Scholar]
  44. Ji, L.M.; Zhang, T.W.; Milliken, K.L.; Qu, J.L.; Zhang, X.L. Experimental investigation of main controls to methane adsorption in clay-rich rocks. Appl. Geochem. 2012, 27, 2533–2545. [Google Scholar] [CrossRef]
  45. Yang, C.; Zhang, J.C.; Tang, X.; Ding, J.H.; Zhao, Q.R.; Dang, W.; Chen, H.Y.; Su, Y.; Li, B.W.; Lu, D.F. Comparative study on micro-pore structure of marine, terrestrial, and transitional shales in key areas, China. Int. J. Coal Geol. 2017, 171, 76–92. [Google Scholar] [CrossRef]
  46. Shao, X.H.; Pang, X.Q.; Li, Q.W.; Wang, P.W.; Chen, D.; Shen, W.B.; Zhao, Z.F. Pore structure and fractal characteristics of organic-rich shales: A case study of the lower Silurian Longmaxi shales in the Sichuan Basin, SW China. Mar. Pet. Geol. 2017, 80, 192–202. [Google Scholar] [CrossRef]
  47. Topór, T.; Derkowski, A.; Ziemiański, P.; Xzczurowski, J.; McCarty, D.K. The effect of organic matter maturation and porosity evolution on methane storage potential in the Baltic Basin (Poland) shale-gas reservoir. Int. J. Coal Geol. 2017, 180, 46–56. [Google Scholar] [CrossRef]
  48. Pommer, M.; Milliken, K. Pore types and pore-size distributions across thermal maturity, Eagle Ford Formation, southern Texas. AAPG Bull. 2015, 99, 1713–1744. [Google Scholar] [CrossRef]
  49. Liu, W.P.; Liu, J.; Cai, M.L.; Luo, C.; Shi, X.W.; Zhang, J. Pore evolution characteristic of shale in the Longmaxi Formation, Sichuan Basin. Pet. Res. 2017, 2, 291–300. [Google Scholar] [CrossRef]
  50. Xu, L.W.; Wang, Y.; Liu, L.F.; Chen, L.; Chen, J. Evolution characteristics and model of nanopore structure and adsorption capacity in organic-rich shale during artificial thermal maturation: A pyrolysis study of the Mesoproterozoic Xiamaling marine shale with type II kerogen from Zhangjiakou, Hebei, China. Energy Explor. Exploit. 2019, 37, 493–518. [Google Scholar]
  51. Guo, S.B.; Mao, W.J. Division of diagenesis and pore evolution of a Permian Shanxi shale in the Ordos Basin, China. J. Pet. Sci. Eng. 2019, 182, 1–9. [Google Scholar] [CrossRef]
  52. Guo, H.J.; Jia, W.L.; Peng, P.A.; Lei, Y.H.; Luo, X.R.; Cheng, M.; Wang, X.Z.; Zhang, L.X.; Jiang, C.F. Evolution of organic matter and nanometer-scale pores in an artificially matured shale undergoing two distinct types of pyrolysis: A study of the Yanchang Shale with Type II kerogen. Org. Geochem. 2017, 105, 56–66. [Google Scholar] [CrossRef]
  53. Hu, H.Y.; Hao, F.; Liu, J.F.; Lu, Y.C.; Ma, Y.Q.; Li, Q. Organic matter-hosted pore system in the Wufeng-Longmaxi (O3W-S1L) shale, Jiaoshiba, Eastern Sichuan Basin, China. Int. J. Coal Geol. 2017, 173, 40–50. [Google Scholar] [CrossRef]
  54. Kuila, U.; McCarty, D.K.; Derkowski, A.; Fischer, T.B.; Topór, T. Nano-scale texture and porosity of organic matter and clay minerals in organic-rich mudrocks. Fuel 2014, 135, 359–373. [Google Scholar] [CrossRef] [Green Version]
  55. Curtis, M.E.; Cardott, B.J.; Sondergeld, C.H.; Chandra, S.R. Development of organic porosity in the Woodford Shale with increasing thermal maturity. Int. J. Coal Geol. 2012, 103, 26–31. [Google Scholar] [CrossRef]
  56. Valenza, J.J., II; Drenzek, N.; Marques, J.M.; Pagels, M.; Mastalerz, M. Geochemical controls on shale microstructure. Geology 2013, 41, 611–614. [Google Scholar] [CrossRef]
  57. Xu, L.W.; Liu, L.F.; Jiang, Z.X.; Chen, L.; Wang, Y. Methane adsorption in the low-middle-matured Neoproterozoic Xiamaling marine shale in Zhangjiakou, Hebei. Aust. J. Earth Sci. 2018, 65, 691–710. [Google Scholar] [CrossRef]
  58. Löhr, S.C.; Baruch, E.T.; Hall, P.A.; Kennedy, M.J. Is organic pore development in gas shales influenced by the primary porosity and structure of thermally immature organic matter? Org. Geochem. 2015, 87, 119–132. [Google Scholar] [CrossRef] [Green Version]
  59. Luo, P.; Zhong, N.N. The role of residual bitumen on the pore structure of organic-rich shales from low to over mature Insight from shale and coal samples after the hydrous pyrolysis. Int. J. Coal Geol. 2020, 226, 92–116. [Google Scholar] [CrossRef]
Figure 1. Structural outline of the research area and the sampling point locations.
Figure 1. Structural outline of the research area and the sampling point locations.
Nanomaterials 11 00527 g001
Figure 2. Mercury intrusion curve characteristics and pore type classification of Mesoproterozoic Xiamaling shale samples in Zhangjiakou, Hebei.4.3. Mesopore and micropore characterization based on gas adsorption. (a,b): the first type of MIP curve and its corresponding pores; (c,d): the second type of MIP curve and its corresponding pores; (e,f): the third type of MIP curve and its corresponding pores).
Figure 2. Mercury intrusion curve characteristics and pore type classification of Mesoproterozoic Xiamaling shale samples in Zhangjiakou, Hebei.4.3. Mesopore and micropore characterization based on gas adsorption. (a,b): the first type of MIP curve and its corresponding pores; (c,d): the second type of MIP curve and its corresponding pores; (e,f): the third type of MIP curve and its corresponding pores).
Nanomaterials 11 00527 g002
Figure 3. N2 adsorption and desorption isotherms of representative samples of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei ((a), the first type of N2 adsorption curve; (b), the second type of N2 adsorption curve; (c), the third type of N2 adsorption curve; (d), the fourth type of N2 adsorption curve).
Figure 3. N2 adsorption and desorption isotherms of representative samples of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei ((a), the first type of N2 adsorption curve; (b), the second type of N2 adsorption curve; (c), the third type of N2 adsorption curve; (d), the fourth type of N2 adsorption curve).
Nanomaterials 11 00527 g003
Figure 4. PV (ac) and SA (df) distribution with pore size derived from the N2 adsorption branch for the isotherms of representative samples of Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Figure 4. PV (ac) and SA (df) distribution with pore size derived from the N2 adsorption branch for the isotherms of representative samples of Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Nanomaterials 11 00527 g004
Figure 5. CO2 adsorption isotherms of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei ((a), TOC < 1.0%; (b), 1.0% < TOC < 3.0%; (c), TOC > 3.0%).
Figure 5. CO2 adsorption isotherms of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei ((a), TOC < 1.0%; (b), 1.0% < TOC < 3.0%; (c), TOC > 3.0%).
Nanomaterials 11 00527 g005
Figure 6. Micropore PV (ac) and SA (df) distribution derived from the CO2 adsorption isotherms of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Figure 6. Micropore PV (ac) and SA (df) distribution derived from the CO2 adsorption isotherms of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Nanomaterials 11 00527 g006
Figure 7. Cross plots of (a) quartz vs. TOC; (b) clays vs. TOC; (c) clay vs. calcite; (d) clay vs. quartz.
Figure 7. Cross plots of (a) quartz vs. TOC; (b) clays vs. TOC; (c) clay vs. calcite; (d) clay vs. quartz.
Nanomaterials 11 00527 g007
Figure 8. Full-scale PV distribution obtained by MIP, low-pressure N2 and CO2 adsorption of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei ((a), multimodal with multiple scale pores; (b), bimodal with primarily mesopores; (c), bimodal with primarily micropores).
Figure 8. Full-scale PV distribution obtained by MIP, low-pressure N2 and CO2 adsorption of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei ((a), multimodal with multiple scale pores; (b), bimodal with primarily mesopores; (c), bimodal with primarily micropores).
Nanomaterials 11 00527 g008
Figure 9. Relationships between the porosity and TOC content (a), quartz (b) and clay (c).
Figure 9. Relationships between the porosity and TOC content (a), quartz (b) and clay (c).
Nanomaterials 11 00527 g009
Figure 10. Relationship between the TOC content and pore structure in Ximaling shale ((a): pore volume; (b) surface area).
Figure 10. Relationship between the TOC content and pore structure in Ximaling shale ((a): pore volume; (b) surface area).
Nanomaterials 11 00527 g010
Figure 11. PVs and SAs of different scales versus the clay mineral content (a,b) and quartz content (c,d).
Figure 11. PVs and SAs of different scales versus the clay mineral content (a,b) and quartz content (c,d).
Nanomaterials 11 00527 g011
Figure 12. TOC-normalized PVs and SAs of different scales versus the (a,b) clay mineral content, (c,d) calcite content, and (e,f) quartz content.
Figure 12. TOC-normalized PVs and SAs of different scales versus the (a,b) clay mineral content, (c,d) calcite content, and (e,f) quartz content.
Nanomaterials 11 00527 g012
Figure 13. Different clay mineral compositions versus the TOC-normalized PV of the (a) micropores, (b) mesopores, (c) macropores and the TOC-normalized SA of the (d) micropores, (e) mesopores, (f) macropores.
Figure 13. Different clay mineral compositions versus the TOC-normalized PV of the (a) micropores, (b) mesopores, (c) macropores and the TOC-normalized SA of the (d) micropores, (e) mesopores, (f) macropores.
Nanomaterials 11 00527 g013
Figure 14. Thermal maturity (EqRo) versus the (a) TOC-normalized PV and (b) SA of the micropores, mesopores, and macropores.
Figure 14. Thermal maturity (EqRo) versus the (a) TOC-normalized PV and (b) SA of the micropores, mesopores, and macropores.
Nanomaterials 11 00527 g014
Figure 15. Evolution of TOC-normalized (a) micropore, (b) mesopore, and (c) macropore volume with the equivalent vitrinite reflectance.
Figure 15. Evolution of TOC-normalized (a) micropore, (b) mesopore, and (c) macropore volume with the equivalent vitrinite reflectance.
Nanomaterials 11 00527 g015
Figure 16. CO2 adsorption isotherms of the original shale and extracted shale ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Figure 16. CO2 adsorption isotherms of the original shale and extracted shale ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Nanomaterials 11 00527 g016
Figure 17. N2 adsorption and desorption isotherms of the original shale and extracted shale ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Figure 17. N2 adsorption and desorption isotherms of the original shale and extracted shale ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Nanomaterials 11 00527 g017
Figure 18. Correlations between gas adsorption content and (a,b,d) TOC content and (c,e) EOM.
Figure 18. Correlations between gas adsorption content and (a,b,d) TOC content and (c,e) EOM.
Nanomaterials 11 00527 g018
Figure 19. Micro-PV of the original shale and its corresponding extracted sample varies with pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Figure 19. Micro-PV of the original shale and its corresponding extracted sample varies with pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Nanomaterials 11 00527 g019
Figure 20. Micropore SA of the original shale and its corresponding extracted sample varies with pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Figure 20. Micropore SA of the original shale and its corresponding extracted sample varies with pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Nanomaterials 11 00527 g020
Figure 21. Micropore SA and PV variation characteristics of the original shale and its corresponding extracted samples.
Figure 21. Micropore SA and PV variation characteristics of the original shale and its corresponding extracted samples.
Nanomaterials 11 00527 g021
Figure 22. Relationships between the EOM and gap value of the (a) micropore volume and (b) surface area.
Figure 22. Relationships between the EOM and gap value of the (a) micropore volume and (b) surface area.
Nanomaterials 11 00527 g022
Figure 23. Meso- and macro-PV of the original shale and its corresponding extracted sample varies with the pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Figure 23. Meso- and macro-PV of the original shale and its corresponding extracted sample varies with the pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Nanomaterials 11 00527 g023
Figure 24. Mesopore and macropore SA of the original shale and its corresponding extracted sample varies with the pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Figure 24. Mesopore and macropore SA of the original shale and its corresponding extracted sample varies with the pore size ((a), ZJS-3; (b), JZT-6; (c), XJG-6; (d), HTG-4; (e), JMS-9; (f), XHY-4).
Nanomaterials 11 00527 g024
Figure 25. Mesopore SA and PV variation characteristics of the original shale and its corresponding extracted samples.
Figure 25. Mesopore SA and PV variation characteristics of the original shale and its corresponding extracted samples.
Nanomaterials 11 00527 g025
Figure 26. Relationships between the EOM and gap value of the (a) mesopore volume and (b) surface area.
Figure 26. Relationships between the EOM and gap value of the (a) mesopore volume and (b) surface area.
Nanomaterials 11 00527 g026
Table 1. Geochemical information and mineralogical composition for the investigated Xiamaling shales.
Table 1. Geochemical information and mineralogical composition for the investigated Xiamaling shales.
Sample IDFormationTOC
(%)
EqRo
(%)
Mineralogical Composition Relative Percent (%)EOM
(%)
QuartzCalciteFeldsparPyriteClaysSmectiteI/SIlliteChlorite
ZJS-1Qnx2.230.8931.510.23.12.452.86.715.227.43.5
ZJS-3Qnx4.480.7727.819.62.32.248.18.215.820.14.00.86
ZJS-6Qnx2.580.8135.28.14.73.648.44.040.04.40
JZT-1Qnx3.410.8337.819.14.05.134.04.920.18.10.9
JZT-6Qnx4.690.8546.115.88.24.925.05.112.36.01.61.02
XJG-3Qnx1.310.9332.210.14.91.950.912.321.89.96.9
XJG-6Qnx2.450.9836.15.23.51.453.87.327.617.41.50.79
XJG-9Qnx1.471.0132.615.76.74.140.914.620.85.50
HTG-2Qnx4.121.4627.228.63.71.339.28.414.416.40
HTG-3Qnx3.671.1827.632.41.71.736.67.719.87.51.6
HTG-4Qnx4.661.4348.424.31.93.122.37.88.55.10.90.94
HTG-5Qnx3.281.8828.529.68.61.132.23.713.113.32.1
JMS-4Qnx0.550.5228.35.89.47.249.312.725.66.34.7
JMS-5Qnx0.970.6221.89.611.34.452.917.926.78.30
JMS-9Qnx6.801.8242.322.511.4023.85.612.13.82.31.42
XHY-3Qnx0.820.5726.712.310.73.147.213.225.15.73.2
XHY-4Qnx7.601.6245.824.35.13.321.53.612.84.11.01.93
XHY-7Qnx0.741.7440.75.42.5051.415.425.86.24.0
Note: “EqRo” indicates equivalent vitrinite reflectance; “–” indicates no data; “I/S” indicates illite-smectite mixed layer; “EOM” indicates extracted OM content.
Table 2. MIP measurement results of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Table 2. MIP measurement results of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Sample IDMacropore PV (cm3/g)Macropore SA (m2/g)Bulk Density
(g/cm3)
Skeletal Density
(g/cm3)
Porosity (%)
ZJS-10.00910.24142.22012.26141.84
ZJS-30.01270.29092.19612.23451.72
ZJS-60.01031.06172.31392.34561.35
JZT-10.00622.67032.24212.26440.99
JZT-60.00910.21222.28022.33322.27
XJG-30.00430.31282.21462.23450.89
XJG-60.00810.21392.30622.35111.91
XJG-90.00421.53412.17912.21631.68
HTG-20.01120.36292.35772.37120.57
HTG-30.00810.58282.25232.26540.58
HTG-40.00430.22272.19872.26252.82
HTG-50.00561.08392.35132.37430.97
JMS-40.00160.56422.33892.36831.24
JMS-50.00411.26182.35812.38251.03
JMS-90.01110.36282.30792.38473.22
XHY-30.00320.55092.63412.66411.13
XHY-40.01360.26422.27552.34542.98
XHY-70.00261.31332.39032.43141.69
Table 3. Main pore structure parameters measured by low-pressure gas adsorption of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Table 3. Main pore structure parameters measured by low-pressure gas adsorption of the Mesoproterozoic Xiamaling shale in Zhangjiakou, Hebei.
Sample IDN2 AdsorptionCO2 Adsorption
VBJH (cm3/g)SBET (m2/g)APD (nm)VDFT (cm3/g)SDFT (m2/g)APD (nm)
ZJS-10.0170111.722.690.0090714.40.885
ZJS-30.0169910.712.340.0188216.40.845
ZJS-3-E0.0188114.616.870.0244918.71.236
ZJS-60.0166716.94.110.0102619.70.554
JZT-10.0171511.416.230.0094219.70.886
JZT-60.0162612.811.630.0193119.80.995
JZT-6-E0.0178416.414.320.0214822.51.689
XJG-30.010754.810.480.0060114.90.596
XJG-60.0193510.74.450.0113714.30.972
XJG-6-E0.0215412.26.890.0163916.61.554
XJG-90.012348.14.220.005848.21.511
HTG-20.0166414.17.840.0166623.80.885
HTG-30.0129313.59.360.0083323.01.287
HTG-40.0128611.48.040.0106318.70.844
HTG-4-E0.0145614.212.890.0194522.31.248
HTG-50.012718.89.660.0080320.41.567
JMS-40.004595.43.980.002117.41.264
JMS-50.009254.25.540.004288.71.563
JMS-90.0208115.48.530.0100229.31.267
JMS-9-E0.0222119.313.050.0208634.31.994
XHY-30.005128.33.660.0042210.10.875
XHY-40.0132414.34.090.0126236.60.866
XHY-4-E0.0196319.39.290.0209341.31.632
XHY-70.005895.13.540.003598.60.894
Note: SDFT = SA using density functional theory; VDFT = PV using density functional theory; SBET = SA using the Brunauer-Emmett-Teller method; VBJH = PV using the Barrett-Joyner-Halenda method; APD = average pore diameter.
Table 4. Gap values between the micropore and mesopore PVs and SAs in the extracted samples and original shales.
Table 4. Gap values between the micropore and mesopore PVs and SAs in the extracted samples and original shales.
Sample IDCO 2 AdsorptionN2 Adsorption
Gap of VDFT (cm3/g)Gap of SDFT (m2/g)Gap of VBJH (cm3/g)Gap of SBET (m2/g)
ZJS-30.0062.3090.0023.890
JZT-60.0102.6200.0023.551
XJG-60.0052.3110.0021.579
HTG-40.0093.6250.0022.807
JMS-90.0115.0600.0013.967
XHY-40.0014.6760.0064.992
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, L.; Yang, K.; Wei, H.; Liu, L.; Li, X.; Chen, L.; Xu, T.; Wang, X. Full-Scale Pore Structure Characteristics and the Main Controlling Factors of Mesoproterozoic Xiamaling Shale in Zhangjiakou, Hebei, China. Nanomaterials 2021, 11, 527. https://doi.org/10.3390/nano11020527

AMA Style

Xu L, Yang K, Wei H, Liu L, Li X, Chen L, Xu T, Wang X. Full-Scale Pore Structure Characteristics and the Main Controlling Factors of Mesoproterozoic Xiamaling Shale in Zhangjiakou, Hebei, China. Nanomaterials. 2021; 11(2):527. https://doi.org/10.3390/nano11020527

Chicago/Turabian Style

Xu, Liangwei, Keji Yang, Hao Wei, Luofu Liu, Xiao Li, Lei Chen, Tong Xu, and Ximeng Wang. 2021. "Full-Scale Pore Structure Characteristics and the Main Controlling Factors of Mesoproterozoic Xiamaling Shale in Zhangjiakou, Hebei, China" Nanomaterials 11, no. 2: 527. https://doi.org/10.3390/nano11020527

APA Style

Xu, L., Yang, K., Wei, H., Liu, L., Li, X., Chen, L., Xu, T., & Wang, X. (2021). Full-Scale Pore Structure Characteristics and the Main Controlling Factors of Mesoproterozoic Xiamaling Shale in Zhangjiakou, Hebei, China. Nanomaterials, 11(2), 527. https://doi.org/10.3390/nano11020527

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop