A natural analogue study in CO2-rich coalfield (Haishiwan, China) provides a strong support for safe, reliable, and long-term storage by analyzing the mechanism of CO2 migration, entrapment, and storage in coal reservoir. Thus, effects of geological tectonism on reservoir properties were investigated. Simultaneously, coal and oil shale samples before and after supercritical CO2 (SCCO2) treatment via geochemical reactor were collected to analyze changes in pore structure, functional group distributions, and SCCO2 extraction. Observations from in situ properties of coal seam indicate that there is a positive relationship with CH4 contents and F19 fault whereas CO2 and carbonate contents decrease as the distance from F19 increases. Analysis of pore properties reveals that SCCO2 enlarges the development of coal pore and facilitates the diffusion and seepage channel of coal reservoir, while no changes in larger pores are found in oil shale, which may restrain fluids from passing through. Then, oxygen-containing functional groups are mobilized by SCCO2 from oil shale, associated with a decrease in sorption sites. The sealing capacity of cap rock (oil shale) and geological tectonism (F19 fault), as the major contributors to CO2 enrichment and accumulation, provides insights into the suitable selection of CCGS site for long geological time.

1. Introduction

Atmospheric CO2 originates mainly from fossil fuel combustion, which is responsible for about 70% of the greenhouse effect [1, 2]. China’s CO2 emissions in 2016 continue to be number one in the world with 27.3% of total [3]. To reduce anthropogenic CO2 emissions, CCGS is put forward as a promising approach by CO2 capture, transport, and injection into a storage site, such as depleted oil and gas reservoirs, saline aquifers, and unmineable coal seams [46]. As one of the most common potential disposal sites, CO2 geological storage in coal mines has a significant meaning related to CO2-enhanced coalbed methane (ECBM) recovery [7]. However, storage in coal seam faces tremendous challenges: (1) the potential geological risks related to seismicity and ground movement in unstable coal reservoirs; (2) the health, safety, and environmental risks causing CO2 leakage to cap rock, tectonic zone, and fractures [8, 9]. Therefore, special attention should be given on the effective and safe storage of CO2 in suitable sites, with the purpose of CO2 leakage control, which serves as a primary risk factor.

In fact, only several pilot CCGS projects have been conducted and it is apparent that CO2 storage in coal is not widely accepted in China [10]. A knowledge on CO2 behavior of long-term storage is limited for the field experiments in coalfields [11]. Previous studies have been demonstrated that CO2 in Yaojie coalfield can be naturally stored over long geological period under suitable conditions [12]. To some extent, the presence of CO2 in Yaojie coalfield can be analogous to long-term storage of CO2, which may yield insights into the mechanisms of CO2 migration, entrapment, and storage in coal reservoirs. Therefore, the suitability of CO2 injection sites and other activities in this area need to be fully considered from the macroscopic and microscopic views.

From the macroscopic view, a group of geological factors including occurrence of coal seam, subsurface condition, cap rock types, and regional tectonic condition should be analyzed and evaluated [1315]. The microscopic view refers to the interactions between CO2 and rocks (coal or cap rock) and the effects of CO2 on the physical (structure of pore or fracture) and chemical (molecules) properties [16]. For coal, the interactions with CO2 are generally divided into three categories: (1) swelling effect, caused by CO2 absorbing into coal matrix, leading to a decrease in the fracture aperture and permeability of coal reservoir [17, 18]; (2) the transformation and rearrangement of coal structure after SCCO2 treatment [19, 20]; (3) SCCO2 extraction of organic matter and inorganic material in coal, accompanied by the existence of water [21]. Although the above findings have been published in the literature over the past decade in the area of CCGS, a knowledge gap exists on the effects of SCCO2 on cap rocks, that is, oil shale in this study. In Haishiwan coalfield, the roof of coal seam is enclosed by a mass of oil shale, which serves as a better sealant for SCCO2 pressure [12]. Actually, cap rock, one of the most important factors, dominates the sealing capacity of coal reservoirs [22]. Hence, an effort should be made to study the changes of oil shale structure after SCCO2 treatment.

The main purpose of this study is to investigate the feasibility of effective and safe geological sequestration of CO2 in coal seams underling cap rock (oil shale), which takes Haishiwan coal seam in China with high CO2 content as a natural analogue. A knowledge of the mechanisms of CO2 migration, entrapment, and storage in Haishiwan coalfield contributes to a better understanding of long-term storage of CO2 in the field of CCGS. Therefore, the structure in this paper is as follows. Firstly, geological factors related to tectonism and its influence on in situ properties of coal seam were investigated and analyzed. Subsequently, samples of coal and oil shale were treated by SCCO2 geochemical reactor to perform comparative tests, including pore properties, functional group distributions, and Gas Chromatography-Mass Spectrum (GC-MS) analysis. Then, the mechanisms of CO2 migration, entrapment, and storage were discussed to demonstrate how to effectively and safely store CO2 in coalfield.

2. Regional Geology

2.1. Geological Setting of Yaojie Coalfield

As shown in Figure 1, Yaojie coalfield, as a secondary tectonic unit in Minhe basin, is located in the west margin of Gansu and Qinghai provinces. There are many well-developed faults and constructions in Minhe, which is a Mesozoic–Cenozoic mountain depression basin expanded by the Qilian orogenic belt. The tectonic evolution of Minhe basin is composed of Indosinian movement, Yanshan movement, and Himalayan movement, resulting in the basin formation and transformation. Yaojie coalfield belongs to middle Jurassic coal-bearing series in Mesozoic. The strike, dip, and area are 2.6 km, 11 km, and 28.6 km2, respectively. The coal bearing basin extended in a NNW direction with the F19 fault zone located to the east of it. Yaojie coalfield, belonging to the Yaojie Coal-Electricity Corporation, consists of Zhangergou coal mine, Yaojie NO. III coal mine, Jinhe coal mine, and Haishiwan coal mine. Haishiwan coalfield, as a concealed coalfield, is located in the south of Yaojie coal fields in Minhe basin. The main coal seam containing CO2 is No. 2 coal seam; the average thickness of which is 20 m with thicker in the west and thin in the east. The roof of No. 2 coal seam is carbon mudstone, siltstone, and oil shale. The floor of No. 2 coal seam is carbon mudstone.

2.2. F19 Faults

There is no direct evidence exactly to show the formation time of CO2 fluids were affected by F19 faults, but the isotopic age chlorite thermal alteration and the fission track age of the calcite veins of ultrabasic have proved that CO2 hydrothermal activity was concentrated in this period [12]. In Neogene period with Himalayan movement, F19 fault was a tensional activity period, the permeability and fluid potential of which increased with the possibility of the upward migration of CO2 fluids [23]. Therefore, CO2 fluids were mainly formed in Neogene period. Whereafter, fault tensional activity quitted or transformed into pressure shear state, making the vertical migration of CO2 fluids slow down to stagnation. Up to the present, the vertical closure and principal compressive stress of NE direction suggest F19 fault has been still in a closed state. For this area, ductile shear fault is classified as damage zone and fault nucleation zone, which is presented in Figure 2. The fault nucleation was filled with breccia, fault gouge, and cataclasite, and the old faults have been blocked by quartz, calcite, and zeolite, which could result in a poor permeability in F19 fault. However, there are amount of minor faults and fractures in damage zone, indicating that good permeability in this area may be the main lateral migration channel of CO2 fluids.

2.3. Sources of CO2

It has been proven that CO2 of Earth’s crust derived from stratal carbonate pyrolysis generating inorganic CO2, which are influenced by contact metamorphism of carbonates and dynamic metamorphism of fault activities [24]. According to geological survey in Geological Institute (Gansu, China), magmatic intrusion has never happened in the region since the Paleozoic [25]. Therefore, the possibility of contact metamorphism can be excluded during inorganic CO2 generating process. However, due to the existence of F19 fault, the fault zone tends to suffer intense extrusion and shear, leading to the hydrothermal activities with structural and tectonic characteristics [26]. Given that carbonate rocks in active fault zone may result in dynamometamorphism under the effect of ductile-brittle shearing action, element differentiation may occur and CO2 will be generated after migration [27]. This is the main formation mechanism of inorganic CO2 in this region.

3. Experimental Section

3.1. Sample Preparation

Coal samples were sampled from the air-return way of coalface in Haishiwan coal mine, while oil shale samples were derived from the roof of No. 2 coal seam. The samples were immediately sealed and sent to the laboratory as soon as possible to prevent oxidation. Subsequently, the samples were crushed and screened to the appropriate quantity and sizes for each test. Before tests, samples were dried at 80°C for 48 h using turbo molecular pump vacuum.

Following China National Standard GB/T 212-2008, proximate analyses of moisture, ash, and volatile matter were performed using a 5E-MAG6600 proximate analyzer (Changsha Kaiyuan Instruments, China). The mean maximum reflectance of vitrinite (, %) was determined using a Zeiss microscope-photometer (German) following China National Standard GB/T 6948-2008. The basic properties of the coal and oil shale samples are listed in Table 1.

3.2. SCCO2 High-Pressure Geochemical Reactor

To better understand the changes in physical and chemical structure, the samples (coal and oil shale) before and after ScCO2 treatment were conducted by SCCO2 high-pressure geochemical reactor, which can simulate CO2 injection in the reservoir (Figure 3). Firstly, samples were degassed in the vacuum pump system (60°C and 4 Pa for 48 h). Secondly, CO2 was injected to pipeline while air compressor started to provide power. 10 min later, gas booster pump was opened and then high-pressure CO2 was filled in the reactor. Eventually, the pressure and temperature of reaction were set at 10 MPa and 50°C, respectively, and the reaction time is 360 h. Actually, This condition could not lead to chemical degradation of the coal matrix [28]. To some degree, the pressure-temperature conditions in this study were used to simulate CO2 sequestration scenarios in coal beds (1 km depth), though it could not be representative of CO2 sequestration in deep coal beds based on field [28].

3.3. Experimental Methods

Pore structure analysis of macropore was performed using A PoreMaster 33 automated mercury intrusion porosimeter (Quantachrome Instruments, United States), which could effectively measure pore diameters of 6 to 100 μm over a pressure range from 0.14 to 215 MPa. The mesopore and micropore analysis were obtained using an AUTOSORB-1 (Quantachrome Instruments, United States). With these two methods, specific surface area (SSA), total pore volume (TPV), pore size distributions (PSDs), and porosity of micropores, mesopores, and macropores can be precisely analysed.

GC-MS analysis was conducted by a Hewlett–Packard 6890/5973 GC/MS associated with a capillary column applied with HP-5MS (cross-link 5% PH MEsiloxane,  i.d., 0.25 mm film thickness), a quadrupole analyzer, and operated in electron-impact (70 eV) mode. The data can be performed automatically by Chemstation software. The compounds were characterized by comparing mass spectra according to NIST05 library data.

FTIR spectra of the samples were measured using a device that combined a VERTEX 80 Fourier transform infrared spectrometer with a HYPERION 2000 infrared microscope (Bruker, Germany) at a resolution of 0.06 cm−1 and spectral region of 8000–350 cm−1.

4. Results

4.1. Relationships with F19 Fault, CO2, and Carbonate Distributions

Figure 1 exhibits the locations of 28 surface boreholes for the No. 2 coal seam during geological exploration period in Haishiwan coalfield. Gas (CO2, CH4) and carbonate mineral contents were obtained from surface borehole coal samples, which were shown in Figure 4.

The on-site CO2 contents are between 0.3 m3/g and 10 m3/g, whereas CH4 contents are between 0.03 m3/g and 7.23 m3/g. Then, a negative correlation was found between CO2 contents and the distance from F19 faults, while it seems to be a positive correlation between CH4 contents and F19 faults. Because CO2 molecule has stronger capacity to diffuse and adsorb than CH4 molecule, the highest content of CO2 may be found near F19 faults. That is, it can be regarded as the input of CO2. It is summarized that CH4 in the coal seam will be recovered after CO2 sequestration through CO2 injection.

The distribution of carbonate mineral ranges from 0.2% to 9.7%. The data in Figure 4 suggests that carbonate contents decrease with increasing distance from F19 fault. According to the literature, carbonate, as a tracing mineral of inorganic CO2 migration, is related to the existence of water, pH value, pressure, and former mineral [29]. The distribution of carbonate mineral indicates that under a certain condition CO2 is likely to transform into carbonate minerals near the F19 fault zone.

4.2. Pore Characteristics

Studies on the pore structure of oil shales and coals indicate the development and connectivity of smaller pores (micropores) and larger pores (macropores) directly relate to the desorption, diffusion, and seepage, respectively [30, 31]. The International Union of Pure and Applied Chemistry (IUPAC) classified pores into micropores (<2 nm), mesopores (2−50 nm), and macropores (>50 nm), which is widely recognized by the scholars [32].

For mesopores, N2 adsorption at 77 K can be regarded as the standard method to determine TPV and SSA using the Barrett−Joyner−Halenda (BJH) and Brunauer−Emmett−Teller (BET) models, respectively [3234]. For micropores, CO2 adsorption at 273 K may overcome the disadvantages that N2 adsorption is inaccessible to ultramicropores (<0.7 nm), which can be interpreted by the density functional theory (DFT) analyses for TPV and SSA [35, 36]. In this study, physisorption method (N2/CO2 adsorption) can be automatically analyzed by the ASiQwin computer software from Quantachrome (United States). However, physisorption method is not suitable for larger macropores (>300 nm) analysis [37]. Thus, mercury intrusion porosimetry (MIP) is more effective for TPV and SSA of macropores through PoreMaster from Quantachrome (United States).

As shown in Table 2, in general, micropores volume and surface area of oil shale and coal samples occupy a large proportion in total pore volume. However, there is less development of macropores, mesopores, and micropores in oil shale, compared to coal sample. Especially for MIP, the porosity of coal is 5.82%, which seems to be more than twice as large as oil shale of 2.31%. It may indicate that coal has a better migration channel than oil shale, due to a more complex and developed pore structure. Also, the SSA of micropores dominate in coal, compared with oil shales, suggesting that the most fluid molecules (CH4 or CO2) can be inclined to adsorb and store in micropores.

Generally, the PSD of micropore and lower pore size range of mesopore can be characterized by N2 adsorption using DFT method [38, 39]. MIP with a function of Incremental pore volume and pore diameter can reflect the PSD of larger pore size range mesopore and macropore. These methods show a good representation from micropore to macropore, which almost cover the whole range of PSDs.

Figure 5 depicts the PSDs of untreated and SCCO2-treated oil shale. No obvious changes were observed in the PSD of MIP, especially in the PSD of macropores. This indicates that SCCO2 has little influence on the larger pores (mesopore and macropore) of oil shale. The inset in Figure 5 highlights that smaller pores (micropore and mesopore) decrease obviously after SCCO2 treatment, suggesting that SCCO2 may alter the inherent structure of oil shale by reducing the TPV [40]. The decrease in micropores and mesopores may be attributed to the adsorption-induced deformation in pores after SCCO2 treatment [41]. Figure 6 shows the PSDs of coal before and after SCCO2 treatment. Unlike the PSD of oil shale, the PSDs of coal in Figure 6 exhibit a bipolar distribution, which contains more micropores and macropores, and mesopores are not developed in coal. This may indicate that coal has a strong ability for gas desorption, diffusion, and seepage, compared with oil shale. Furthermore, the PSD of MIP after SCCO2 treatment shows more macropores and mesopores than untreated coal sample, except for the pore diameter range around 100 nm. For DFT method, two well-defined peaks and some of disorderly peaks are found in micropores and mesopores, suggesting the well-developed micropores in coal structure and better adsorbed ability. With SCCO2 treatment, micropores and mesopores show an increase in the PSDs of DFT method. Eventually, SCCO2 promotes the development of micropores, mesopores, and macropores in coal, which may facilitate the movement of fluids from the coal matrix via desorption and diffusion and transport in the fracture via seepage.

The cumulative pore volume curves of the untreated and SCCO2-treated oil shale analyzed by MIP are displayed in Figure 7. In general, the concave shape is shown in the ejection curves, which may result from a mass of open pores and semiclosed pores. This may contribute to the connectivity of pores. Then, the injection volume of untreated and SCCO2 treated oil shale shows similar trends with the same increases, which means the same amount of mercury was injected into pores before and after SCCO2 treatment. Also, the changes in the ejection volume for oil shale show a decrease of 6.28% (0.00175 to 0.00164 cc/g), indicating a decrease in effective pores after SCCO2 treatment. Therefore, the hysteresis loop becomes narrow with SCCO2 treatment, which may result from the decrease in semiclosed pores. The reason may be attributed to adsorption-induced deformation and physical constriction caused by the SCCO2 dissolving and mobilizing hydrocarbons [42, 43]. Taken in conjunction with the little changes in macropores and a decrease in micropores and mesopores of oil shale as mentioned above, the results can be summarized as the decrease in the characteristics and connectivity of the pore structure related to desorption, diffusion, and seepage.

4.3. Functional Group Distribution

Infrared spectrum analysis has been widely used to study the macromolecular structure of the reservoir rocks, which is derived from the important information regarding the functional groups of organic compounds. Due to the specific molecular vibration spectrum of functional groups, the molecular structure of oil shale can be confirmed by the position and intensity of absorption peaks in the infrared spectrogram. To better understand the sealing ability of cap rocks, changes in the chemical structure of the original and SCCO2 untreated oil shale above the coal reservoir were analyzed using FT-IR.

Generally, the infrared spectrum is split into the functional group region (spectral peak range: 4000–1500 cm−1) and the fingerprint region (spectral peak range: 1500–400 cm−1), which are shown in Figure 8. After SCCO2 treatment, the absorption bands of oil shale change significantly in the region of –CH2, –CH3 symmetrical or antisymmetric stretching vibrations (2935–2918, 2858–2847, and 2882–2862 cm−1), and –CH3 antisymmetric deformation vibrations (1460–1435 cm−1) in naphthenic or aliphatic bonds. Also, the spectra are characterized by a broad hydroxyl region (3697–3610 cm−1), carbonyl in aldehyde, ketone, and ester (1736–1722 cm−1), C–O in epoxy compounds or ethers (1330–1060 cm−1), and bands in the carboxyl region (1715–1690 cm−1). Moreover, the spectra contain C=C stretching vibrations of aromatic or fused rings (1635–1595 cm−1), Si–O stretching vibrations (1060–1020 cm−1), –SH with a peak at 475 cm−1, and –S–S– with a peak at 540 cm−1.

As mentioned above, the spectra, corresponding to functional groups, are grouped into three categories: (1) methyl, methylene in naphthenic, or aliphatic bonds; (2) oxygen-containing functional groups, such as hydroxyl, carbonyl, and carboxyl; (3) other functional groups. Generally, oil shale, similar to coal, has a dual pore structure that includes both matrix pores and fractures and can adsorb gas or fluids in the surface area. Thus, the first category may be thought as aliphatic and aromatic hydrocarbons in the pore structure or matrix of oil shale. The second category refers to oxygen-containing functional groups on the external surface of oil shale. The third category may be related to the silicate minerals (Si–O), organic sulfur compounds (–SH), and so on. Therefore, SCCO2 may extract some hydrocarbons or other containing functional groups compounds in the matrix pores or fractures, resulting in the changes on the surface site and body structure of oil shale.

4.4. GC-MS Analysis

An investigation of the potential of SCCO2 extraction was carried out for coal sample. The organic compounds, such as hydrocarbons or other substances, were presented clearly by the GC-MS detector. The partial results are exhibited in Figure 9.

According to the analysis of Chemstation software, SCCO2 is capable of mobilizing small hydrocarbon molecules from the coal matrix. The hydrocarbons in the coal sample extract consisted largely of the small molecular weight n-alkanes, aromatic hydrocarbons (C6H6), and, most notably, nC5, nC6, nC7, nC8, and nC9. These results are consistent with the research of Kolak and Burruss [21]. That is, SCCO2 mobilized comparable amounts of hydrocarbons from coal sample. According to the previous studies, the SCCO2 extraction of hydrocarbons is classified as aliphatic hydrocarbons and polycyclic hydrocarbons.

It has been demonstrated that hydrocarbons, existing in the macropores, mesopores, and micropores of coal structure, may generate high activation energy barriers within the coal matrix, which may hinder gas or fluids from transporting through physical constrictions. Thus, gas molecules (such as CO2) have insufficient capacity to pass through coal beds. However, once it interacts with SCCO2, hydrocarbons and mineral matter may be extracted by SCCO2, contributing to reopen pores and fractures. Given the results of pore structure and GC-MS analysis, important information need to be put forward that SCCO2 may facilitate coal permeability and fluid transport, resulting in the development of pore structure and migration channel of coal.

5. Discussion

5.1. Contributions of Coalbed to CO2 Migration for Geological Sequestration

The measured CO2 pressures in No. 2 coal seam are more than the critical pressure (7.38 MPa), and reservoir temperature of No. 2 coal seam is over critical temperature (34°C) according to the calculations, which have been demonstrated by previous studies [44]. The investigation has confirmed that SCCO2 can occur in the coal seam, even in the roof of coal seam (under cap rock) [45]. Because the source of inorganic CO2 is derived from crust and flows along F19 faults, which will be discussed in detail in Section 5.2, CO2 may migrate and accumulate to SCCO2 in the coal seam. In this process, coal fractures may initially narrow and a decrease of permeability may occur in coal reservoir. However, the results of coal pore properties and GC-MS analysis from Figures 59 indicate that SCCO2 enlarges the fracture width and promotes pore development by extracting hydrocarbons and mineral matter. The reasons are as follows: the migration channel, that is, fracture aperture, is a complex factor, which is related to external stress, CO2 injection pressure, adsorption swelling effect, and SC-CO2 extraction [46]. After a long time of evolution, an irreversible change may occur in the channel for gas diffusion and seepage, contributing to the permeability of coal reservoir [16]. Eventually, the influence of SCCO2 on physical and chemical of coal may occur companied with the changes of pore structure, CH4 displacement, and carbonate evolution after a period of CO2 injection.

5.2. Effect of Oil Shale on Seal Capacity of CO2 Sequestration

The findings of FTIR analysis in this study confirm that the amount of oxygen-containing functional groups were extracted from oil shale during SCCO2 injection. Generally, oxygen-containing functional groups dominate the adsorption capacity through the sorption sites [47]. The decrease of oxygen-containing functional groups in oil shale may weaken the capacity of CO2 adsorption, causing the decrease in the CO2 sorption sites. The direct effect may cause less free CO2 changing into adsorbed state. With migrating towards oil shale, more free CO2 may gather beneath the cap rock. Section 5.1 has confirmed that SCCO2 may facilitate the diffusion and seepage channel of coal reservoir. In this case, plenty of hydrocarbons after SCCO2 extraction can move in the channel of coal bed, accompanied by free CO2 flow. After a long injecting time, hydrocarbons will gather below oil shale. Studies indicate that hydrocarbons may restrict physical access (fractures and larger pores) to oil shale matrix [48]. A large amount of hydrocarbons and other mineral matters, existing in the oil shale structure, may have a negative effect on the permeability, which in turn can hinder CO2 escaping from cap rock.

5.3. Implications of Tectonic Deformation on CO2 Entrapment

Original coal seam is known to be a typical dual-porosity system that contains amount of coal matrix blocks divided by face cleat and butt cleat, as shown in Figure 10(a). Coal matrix, as a part of the blocks, consists of hierarchical pores with different morphology. However, because of the long-term geological formation movement, some areas may be subjected to structural changes. Taking Haishiwan coalfield as an example, undeformed coal and ductile deformed coal can be distinguished clearly in the field, which is shown in Figure 10(b). With fault tectonism, the migration of tectonic deformation may occur locally in coal seam, leading to the existence of coal reservoir and surrounding rocks. Moreover, a great quantity of winding and rub can not only produce more mylonite or other broken coals but also seal original cleats, micro-fractures, and seepage pores for CO2 transport [49]. Therefore, poor connectivity and low permeability in the area of ductile deformed coal seam may block CO2 migration along the bedding, which is liable to preserve more CO2 fluids. This phenomenon serves as a role for CO2 entrapment although geological sequestration in coal depends on multiple factors.

5.4. How to Effectively and Safely Store CO2 in Haishiwan Coalfield

From a sequestration viewpoint, effective and safe storing in coal reservoir refers to the stable injection, migration, and storage of CO2 in the appropriate spots for avoiding CO2 leakage to cause Environmental Safety and Health (ES&H) issues. In Haishiwan coalfield, the favourable conditions of gas sources, migration and conduction of reservoir, sealing ability of cap rock, traps, and storage can bring about a steady environment of CO2 gas reservoir.

According to the source and occurrence regularity of CO2 in Yaojie coalfield as mentioned above, high concentration of CO2 in No. 2 coal seam is derived from the indirect accumulation mode with gas production, gas guidance, gas migration, and gas storage through the brittle-ductile shear dynamic metamorphism of the active F19 fault zone [50]. Based on the analysis of F19 fault tectonic movement and CO2 distribution in No. 2 coal seam, F19 fault played a negative role on gas source and gas storage and leads to gas entrapment during the process of CO2 accumulation.

Large amount of inorganic CO2, produced by the dynamic metamorphism of the active F19 fault zone, could move upwards along F19 fault and accumulate in the border of No. 2 coal seam and F19 fault. Driven by pressure gradient, CO2 flow through the bedding orientation of No. 2 coal seam [51, 52]. With the continuous injection of CO2, the pressure of coal reservoir increased significantly, generating SCCO2 condition. The interaction between coal and SCCO2 may contribute to the development of pore structure and fracture system of coal seam, which may provide conditions for the further migration. Simultaneously, CO2 migration may provide a vector by which Polycyclic Aromatic Hydrocarbons (PAHs) and other hydrocarbons are extracted from coal matrix and transported to the direction of cap rock (oil shale) [18, 21]. These may cause the decrease in the permeability of oil shale by these substances obstructing seepage channel, which contributes to the increase in sealing capacity of oil shale [45]. Besides, oil shale and tight sandstone in the roof and floor of No. 2 coal seam have a low permeability. Also, the coal seam near cap rock may has a low permeability due to the swelling effect of CO2 adsorption in coal matrix. The above factors may result in the CO2 accumulation below the cap rock (the roof of No. 2 coal seam), which is exhibited in Figure 11. When CO2 fluids transport from injection point to coal reservoir, competitive adsorption between CO2 and CH4 molecules occurs, leading to CH4 desorption and migration. However, the sealing ability of cap rock (oil shale) and geological tectonism (F19 fault) contributes to gas enrichment and accumulation during a long geological time. Therefore, on-site conditions, that is, the existence of cap rock and geological tectonism, need to be considered, which have a positive implication for safe and effective CO2 geological sequestration. The stable regional tectonic conditions and cap rock with low permeability are key factors for Haishiwan coalfield, providing an insight into the suitable selection of CCGS site.

6. Conclusions

Natural analogue studies of Haishiwan coalfield demonstrated that CO2 can be naturally stored in deep coal seams for long period on condition that cap rock (oil shale) with good sealing capacity and stable regional tectonism, which contribute to gas entrapment and accumulation during a long geological period. Simultaneously, No. 2 coal seam, influenced by consecutive CO2 injection and existence of SCCO2, may generate amount of fractures and will facilitate fluids migration. To study the suitable selection of CCGS site, in situ conditions and experiments were investigated and analyzed through geological tectonism on reservoir properties, pore structure, functional group distributions, and SCCO2 extraction. Major findings are summarized as follows.

(1) CH4 contents increase and CO2 contents decrease with increasing distance from F19 faults, indicating that CH4 can be recovered after CO2 injection. Also, higher carbonate content is closer to F19 faults, which may be derived from the transformation of CO2 near tectonism.

(2) Changes in smaller pores and seepage-flow pores (mesopores and macropores) are facilitated by SCCO2, but macropores in oil shale show no significant change and a slight decrease in micropores and mesopores was found after SCCO2 treatment. Thus, transport in coal matrix and seepage ability of coal seam will be promoted. Meanwhile, changes in pore structure properties of oil shale can lead to a low permeability and connectivity, which hinder the movement of fluid desorption, diffusion, and seepage.

(3) For FTIR analysis, amount of matters containing hydrocarbons or functional groups compounds in oil shale matrix are mobilized by SCCO2, causing the changes on the surface site and body structure. Then, aliphatic hydrocarbons and polycyclic hydrocarbons obtained from SCCO2 extraction enhance the development of pore structure of coal, generating more migration channel.

Conflicts of Interest

The authors declare no competing financial interests.


The authors are grateful to the Fundamental Research Funds for the Central Universities (Grant no. 2015XKMS006) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).