Next Article in Journal
Stages of Development and Solvents Determine the Anticancer Potential of Mountain Arnica (Arnica montana L.) Inflorescence Extracts
Next Article in Special Issue
Experimental Study of the Multiple Fractalisation of Coal and Rock Failure Subjected to the Coupled Effects of Water, Temperature and Dynamic Loads
Previous Article in Journal
Optimising Amber Processing Using 3D Scanning: New Perspectives in Cultural Heritage
Previous Article in Special Issue
Study on Active Support Parameters for Surrounding Rock with Ultra-Large Span Open-Off Cut in Thick Coal Seam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Simulation Study on the Effect of Co-Associated Minerals on Methane Adsorption and Mechanical Properties of Coal

1
Chinese Institute of Coal Science, Beijing 100013, China
2
Key Laboratory of In-Situ Property Improving Mining of Ministry of Education, Taiyuan University of Technology, Taiyuan 030024, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(24), 12975; https://doi.org/10.3390/app132412975
Submission received: 13 November 2023 / Revised: 30 November 2023 / Accepted: 2 December 2023 / Published: 5 December 2023
(This article belongs to the Special Issue Advanced Methodology and Analysis in Coal Mine Gas Control)

Abstract

:
When rockbursts and coal and gas outbursts simultaneously occur in a coal mine, changes in gas adsorption (concentration of ambient methane) and displacement of coal and rock must occur. The co-associated minerals in coal reservoirs can affect the mechanical properties and methane adsorption capacity, which are commonly disregarded. It is important to construct compound molecular structure models of coal and rock and conduct molecular dynamic simulations to gain a microscopic understanding of underground disasters. In this work, the molecular structure models of anthracite and coking coal–rock compound models containing different contents of calcite and kaolinite were constructed, and the methane adsorption amount and mechanical properties considering temperature, pressure, and mineral contents were simulated and analysed. The results showed that the methane adsorption amount of the compound models increased rapidly, then increased moderately, and stabilized eventually with increasing adsorption pressure, and the Langmuir fitting findings were good. The saturation adsorption amount of methane in the coal models linearly decreased with increasing temperature, while the methane adsorption heat increased. The presence of minerals adsorbed a certain amount of methane, and the methane adsorption amount increased with increasing mineral contents. The mechanical properties of coal molecules changed when mineral molecules such as calcite and kaolinite were present, which had opposite contribution effects. The addition of kaolinite minerals to the coal molecular model always increased the bulk modulus and shear modulus, while the addition of calcite decreased the bulk modulus of the anthracite, causing an increase in the brittleness of the models. The results of the study further explain the adsorption behaviour and mechanical properties of methane in coal and minerals.

1. Introduction

Since the 21st century, coal mining technology has been greatly upgraded, integrated mining and caving mining methods have been promoted, and the number of rockburst coal mines in China increased to 121 by 2008. With the increasing coal mining intensity and depth in China, the number of mines with rockbursts rapidly increased from 142 in 2012 to 177 in 2017 [1]. The accident investigations revealed that two disasters, a rock burst and a coal and gas outburst, simultaneously occurred. For example, the average mining depth of the Pingdingshan No. 12 Coal Mine reached 1100 m, and the rockburst and coal and gas outbursts occurred on 29 June 2005 and 19 March 2006 [2]. The occurrence mechanism of compound disasters is complex, and the prediction index for a single disaster may be invalid [3]. Research on the mechanism of compound disasters is one of the most fundamental research contents in the study of compound disasters and provides a theoretical basis for the evaluation, prediction, and prevention of compound disasters [4]. Just like the human gene sequence, the molecular structures determine the thermodynamic, optical, electromagnetic, and surface properties of coal. The mechanism of methane adsorption and mechanical properties can be reflected through molecular dynamic simulation, so varied coal molecular structure models under different geological conditions need to be built urgently, and the effects of pressure, temperature and mineral contents on methane adsorption need to be clarified. To clarify the mechanism of rockbursts, coal and gas outbursts, and the relationship between the two disasters, this study constructs coal–rock compound structure models, which are intended to provide a reference for the dynamic mechanism of compound disasters.
Analysis and interpretation of the occurrence mechanism of rockburst include strength theory [5], stiffness theory [6], energy theory [7], impact propensity theory [8], “three criteria” theory [9], deformation system instability theory [10], and “three factors” theory [11]. There are three important factors in the occurrence of a rockburst hazard, the propensity to impact is the intrinsic factor, the high-stress concentration and dynamic disturbance is the source factor, and the presence of soft layers in the coal rock is the structural factor. Qi et al. pointed out that a coal rock structure with impact tendency deforms under the action of high stress, forming a high-stress concentration and locally gathering energy. Under the disturbance of mining stress, stick–slip and the release of a large amount of energy along the weak surface or contact surface of the coal rock structure will cause a rockburst disaster [12]. Analysis and interpretation of the occurrence mechanism of coal and gas outbursts include the comprehensive action hypothesis [13], rheological hypothesis [14], spherical shell instability hypothesis [15], and solid–fluid coupling instability theory [16]. Wang et al. reported that coal and gas outbursts are the result of the comprehensive action of stress, the gas contained in the coal seam, and the physical and mechanical properties of the coal seam itself [17]. Field monitoring data, laboratory tests, and molecular dynamics simulation results confirmed that a large amount of gases such as CH4 and CO burst out when the shear toughness zone of the coal seam was damaged [18,19,20,21]. Therefore, the construction of the atomic representation of coal–rock compound structure models and the adsorption characteristics and mechanical properties of methane in the model need further development.
Molecular dynamics simulation was first proposed by Alder and Wainwright in 1957 [22]. To date, molecular simulation technology has been successfully applied in the chemical [23], materials [24], and energy [25,26] industries. In 1942, the University of Pennsylvania built the first coal molecular structure model [27]. Hundreds of coal molecular structure models have since been successfully constructed worldwide. Recently, molecular dynamics simulation in coal and rock models has become an important tool for analysing the adsorption and desorption behaviour of CH4, CO2, and other gases in extreme environments [28,29,30,31,32]. In addition, Wang et al. used molecular dynamics to simulate the change in coal molecular structures and the type of gas generated under shear stress and pointed out that shear gas production is the main source of excessive gas emissions when coal and gas outburst disasters occur [20,21]. To explore the mineralization mechanism of coal-bearing graphite, Ma et al. constructed a high-rank coal molecular structure model of Fengxian in Shaanxi [33]. The CH4 adsorption capacity of coal is influenced by different dynamic and static factors. The dynamic factors generated by the redistribution of the initial rock stress by coal mining or methane extraction, static factors such as coal ranks, moisture contents, temperatures, and minerals of methane adsorption in coal also have attracted extensive attention. Molecular dynamics simulation of coal or rock reflects the influence mechanism of temperature, pressure, water content, and other parameters on the molecular structures and the adsorbed gas therein through volume expansion, interaction energy, and radial distribution function. Zhang et al. constructed an anthracite molecular structure model of the Qinshui Basin and simulated the methane adsorption capacity under different adsorption pressure conditions and the changes in the mechanical properties of anthracite after saturated adsorption [30].
In this work, atomistic representations of a coal–rock compound structure of anthracite in the Qinshui Basin and coking coal in the Ordos Basin were constructed, and molecular dynamic simulations were performed to characterize methane adsorption and mechanical properties considering temperature, pressure, and mineral content. The response characteristics of different associated mineral types and contents to the adsorption characteristics and mechanical properties of methane in coal from the molecular level were explained, and a micro mechanism for the occurrence mechanism of rockburst and coal and gas outburst combined disasters was provided.

2. Molecular Structure Models

Combined with elemental analysis, Fourier transform infrared (FTIR) spectrometry, 13C NMR spectroscopy, X-ray photoelectron spectroscopy (XPS), and transmission electron microscopy of coal samples, the molecular structures of anthracite and coking coal samples, are constructed [34], so the detailed construction process has been overlapped. According to the test results, the 3D models of the molecular structure of the anthracite and coking coal are constructed in Materials Studio software 2020 (as shown in Figure 1). The chemical formulas of the anthracite and coking coal are C220H82N2O13S2 and C213H114O16N2S, respectively, and the dimensions of their 3D models are 1.53 × 1.53 × 1.53 nm and 1.56 × 1.56 × 1.56 nm, respectively. The minerals on the coal surface are mainly kaolinite and calcite [35], thus, calcite and kaolinite molecular are added to the molecular structure models of the anthracite and coking coal to build the coal–rock compound molecular structure model with different mineral contents in Materials Studio software 2020 [36]. It is noted that the methane adsorption on the edges of kaolinite and calcite should be studied before the adsorption simulations in the compound models; the omission of this section is due to the methane isothermal adsorption on a single pure mineral surface (0 0 1) and coal surfaces having been presented and validated in our previous work [37].
In Figure 1, carbon, hydrogen, nitrogen, oxygen, sulphur, calcium, silicon, and aluminium atoms are represented by grey, white, blue, red, yellow, green, dark yellow, and magenta, respectively. Eight coal–rock compound molecular structure models with different kaolinite and calcite mineral contents were constructed. When adding calcite with mass fractions of 25.08w% and 40.11 wt% to the anthracite coal model, the chemical formulas of the model are C230H82N2O43S2Ca10 and C240H82N2O73S2Ca20, and the model sizes are 1.68 × 1.68 × 1.68 nm and 1.81 × 1.81 × 1.81 nm. When kaolinite with mass fractions of 24.18w% and 38.94 wt% is added to the anthracite coal model, the chemical formulas of the model are C220H90N2O45S2Al8Si8 and C220H98N2O77S2Al16Si16, with model sizes of 1.68 × 1.68 × 1.68 nm and 1.80 × 1.80 × 1.80 nm. When adding calcite with mass fractions of 25.08w% and 40.11 wt% to the coking coal model, the chemical formulas of the model are C230H82N2O43S2Ca10 and C240H82N2O73S2Ca20, and the model sizes are 1.72 × 1.72 × 1.72 nm and 1.85 × 1.85 × 1.85 nm. When kaolinite with mass fractions of 24.40w% and 39.17 wt% is added to the coking coal model, the chemical formulas of the model are C213H122O48N2SAl8Si8 and C213H130O80N2SAl16Si16, with model sizes of 1.72 × 1.72 × 1.72 nm and 1.84 × 1.84 × 1.84 nm. Note that the size of the coal models created in the manuscript was indeed small (<2 nm), and readers should critically refer to refer to the simulation results.

3. Computational Methods

3.1. Simulation Scheme

As mentioned above, the minerals on the coal surface are mainly kaolinite and calcite; thus, calcite and kaolinite molecules were added to the molecular structure models of anthracite and coking coal to build the coal–rock compound molecular structure model with different mineral contents. The molecular dynamics simulation schemes include two coal molecular structure models and eight coal–rock compound models. For each coal molecular model, five adsorption temperatures (273.15, 283.15, 293.15, 303.15, and 313.15 K) and thirteen adsorption pressures (0.1, 0.3, 0.5, 1, 2, 3, 4, 5, ..., 10 MPa) were selected and calculated, Refs. [30,32], as showed in Table 1. Under each constant temperature, methane adsorption simulations with different adsorption pressures were carried out in each coal molecular model. Each coal–rock compound model was used to simulate different methane adsorption pressures at 273.15 K.

3.2. Implementation of Molecular Simulations

The simulation details of molecular dynamics have been elaborated in detail in our previous work [30,37]. The force field utilized in the model construction and methane adsorption simulations was the Dreiding force field, and the simulations were completed by the Amorphous cell module, Sorption module, and Force module in the Materials Studio software 2020. The mechanical property simulation was completed in Mechanical Properties in the Force module, also using the Dreiding force field [38]. Note that a suitable force field selection is the key to the accuracy of molecular simulation results. The force field form selected in the simulation was only one of many force fields [38,39,40,41,42,43,44,45], and the simulation parameters and results were expected to be critically referenced by researchers. The cutoff distance for methane adsorption is 12.5 Å, and the amount of absolute adsorption was collected. To calculate the saturation sorption amount of methane, VL, on modified coal, the Langmuir isothermal sorption equation was fitted to the test results.
V a b = V L P P L + P
where P is the methane gas pressure, MPa; Vab is the methane adsorption volume of coal at pressure P, mL/g; VL is the methane saturation adsorption volume of coal at a certain temperature, mL/g, and PL is the methane gas pressure at half of the saturation adsorption volume, MPa. The simulation results of mechanical properties were mainly presented in the form of a stiffness matrix, and then the mechanical parameters were obtained.

4. Results and Discussion

4.1. Methane Adsorption in Coal Models

To verify the correctness and applicability of the models and the applied force field, the adsorption of methane on the coal was fitted by the Langmuir equation. Figure 2 shows the methane isotherm adsorption curves and Langmuir fitting parameter, a, of the anthracite and coking coal at 273.15–313.15 K. When the adsorption temperature is constant, the adsorption capacity of methane molecules in the coal molecular structure model rapidly increases and then stabilises with increasing adsorption pressure, which is consistent with the Langmuir adsorption equation. When the adsorption pressure is constant, the amount of methane adsorption decreases with increasing temperature. In the molecular structure model of anthracite, the saturated adsorption amount of methane is 5.80 N/uc at 273.15 K; it decreases to 4.85 N/uc at 313.15 K. In the molecular structure model of coking coal, the methane adsorption amount is 5.23 N/uc at 273.15 K; it decreases to 4.78 N/uc at 313.15 K. Regarding the methane adsorption amount in the molecular structure model of coking coal under different temperature and pressure conditions, Zhu et al. built a coking coal model of the Chiyu coal mine and also found that the methane adsorption amount decreased with the increase in temperature [46]. A probe with a diameter equal to the methane molecular diameter (0.38 nm) was used to detect and count the micropores of the two coal molecular structure models. The results show that the free volumes of methane adsorbed by the anthracite and coking coal are 16.25 × 10−3 and 8.98 × 10−3 nm3/uc, respectively, and that the accessible surface areas of methane are 0.45 and 0.56 nm2/uc, respectively. The simulation results are consistent with Ref [47], and Cheng et al. reported that the micropores in coal are the main places where methane occurs, that is, the larger the micropore volume, the greater the methane adsorption capacity.
The adsorption heat of methane in coal is an important thermodynamic parameter for characterizing the adsorption behaviour of methane in coal; it serves as a quantitative indicator to evaluate the adsorption affinity of coal to methane [48]. The isothermal adsorption heat of methane in coal is generally calculated by the Clausius–Clapeyron equation [46,48]. Based on the methane adsorption as an exothermic process, Tang et al. used negative values to characterize the methane adsorption heat [48]. The simulation results show that the adsorption heat is expressed as a positive value. Figure 3 shows the isothermal adsorption heat simulation results of methane in the anthracite and coking coal at 273.15–313.15 K. The methane adsorption heat of anthracite is 28.13–28.57 Kj/mol, and that of coking coal is 26.23–26.48 Kj/mol. As the adsorption temperature increases from 273.15 K to 313.15 K, the methane adsorption heat of both coals shows an increasing trend, which is consistent with the results of Tang et al. and Zhang et al. [48,49].

4.2. Effect of Minerals on Methane Adsorption of Coal-Rock Models

The adsorption of methane in the coal–rock compound model can also be fitted by the Langmuir equation [49]. To express the adsorption results clearly, a temperature of 273.15 K was selected to simulate the methane adsorption of compound coal–mineral models. Figure 4 shows the methane isotherm adsorption curves and Langmuir fitting parameter, a, in anthracite–rock and coking coal–rock compound models at 273 K. When the adsorption temperature is constant, the adsorption capacity of methane molecules in the coal–rock compound molecular structure model rapidly increases and then stabilizes with increasing adsorption pressure, which still conforms to the Langmuir adsorption model. When the adsorption pressure is constant, the methane adsorption capacity increases with increasing mineral content. In the anthracite–rock molecular structure model, the saturated adsorption capacity of methane reaches a maximum of 39.82 wt% calcite, which is 24.82 N/uc. In the coking coal–rock molecular structure model, the methane adsorption reaches a maximum of 40.11 wt% calcite, which is 15.18 N/uc. A probe with a diameter equal to the methane molecular diameter (0.38 nm) was selected to detect and count the micropores of the two coal–rock molecular structure models. In the molecular structure model of anthracite–rock, when the mineral content is 39.82 wt% calcite, the micropore volume reaches 0.70 nm3/uc, which shows the large amount of methane adsorption when minerals are present. The results indicate that the highly developed micropores in clay minerals in coal are the main places for methane adsorption in coal, which is also reported by Feng et al. [50]. As can be seen in our previous paper [37,49] and the simulation results from Wang et al. [29], the surface of clay minerals like kaolinite also has the capacity to adsorb some methane. Besides, Feng et al. reported that, from the Scanning Electron Microscopy (SEM) micrographs of the same coal sample with different densities determined by the X-ray CT scan, the mesostructures of cell cavity pores with non-compact packing of the clay minerals appear to be the primary sites of methane adsorption in coal. Therefore, it is more generally accepted by researchers that the microporosity and surface structure parameters in clay minerals jointly affect the adsorption capacity.
The isothermal adsorption heat of methane in the coal–rock compound model can also be calculated by the Clausius–Clapeyron equation [49,51]. Zhang et al. used the methane adsorption heat to characterize the energy release information of methane in kaolinite and determined that the reduction in the adsorption heat under the condition of high water content means that the interaction energy between methane molecules and kaolinite molecules was weakened [49]. Figure 5 shows the isothermal adsorption thermal simulation results of methane in the anthracite–rock model and coking coal–rock model at 273.15 K. When calcite and kaolinite minerals are added, the methane adsorption heat of anthracite is 28.13–23.09 Kj/mol, and that of coking coal is 26.23–23.28 Kj/mol. With an increase in the mineral content, the methane adsorption heat of the two coal rock compounds decreases, and the contribution of calcite is greater than that of kaolinite.

4.3. Mechanical Properties of Coal-Rock Models

The bulk modulus and shear modulus in the coal–rock compound models were simulated in the Mechanical Properties of the Force module. Figure 6 shows the simulation results of the bulk modulus and shear modulus of the anthracite–rock model and coking coal–rock model under different calcite and kaolinite mineral contents. The bulk modulus of the anthracite–rock model decreases with the addition of calcite and tends to decrease with increasing calcite mineral content. The bulk modulus increases with the addition of the kaolinite mineral and tends to increase with increasing mineral content. The bulk modulus of the coking coal–rock model increases with the addition of calcite and kaolinite and tends to increase with increasing mineral content. Note that under the same mass percentage, the influence of kaolinite on the change in the bulk modulus and shear modulus of the two coal rock models is greater than that of calcite on the change in the two mechanical properties. The ratio of bulk modulus to shear modulus (E/G) is often used to estimate the brittle or ductile behaviour of materials. A high E/G represents plasticity, and a low E/G represents brittleness [52]. When minerals exist, the brittleness of the molecular model increases, and the brittleness index increases with increasing mineral content [53]. In addition, the atomic representation of coal–rock compound models is presented on a nanometre scale, which cannot illustrate the macroscopic mechanism of rockbursts and coal and gas outbursts. Larger coal–rock compound models will be constructed in our future work.

5. Conclusions

(1)
The atomistic representations of the coking coal model and anthracite coal model considering the influence of minerals were constructed, which is an approach to characterize the methane adsorption and mechanical characteristics in coal.
(2)
The amount of methane adsorption follows the order anthracite coal > coking coal, and the presence of minerals increases the methane adsorption capacity; the increased amount of methane adsorption follows the order calcite > kaolinite.
(3)
The presence of calcite and kaolinite greatly increased the shear modulus of compound coal and mineral models. Notably, after the addition of calcite and kaolinite minerals, the brittleness of the model increases, and the brittleness index increases with increasing mineral content.
(4)
Only calcite and kaolinite were considered typical co-associated minerals in coal, which is limited to explaining the methane adsorption and mechanical properties in all the coal–rock compound models. In addition, there are limited coal molecular structure models for molecular dynamic simulation work, and the construction of diversified coal molecular structure models is difficult but imperative. Thus, based on Graph Representation Learning and Graph Neural Networks, our team is committed to building diverse molecular structure models for almost all coal and rock, providing a theoretical basis for coal and rock dynamic disasters.

Author Contributions

Methodology, writing—original draft, X.Z.; writing—review and editing, Y.L.; investigation, T.Z.; conceptualization, writing—review and editing, Y.C.; formal analysis, B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the youth project of National Natural Science Foundation of China (No. 52204220, No. 52304088, No. 42302202), the Science and Technology Innovation and Entrepreneurship project of TDTEC (No. 2022-QN001, No. 2022-2-MS003, No. 2023-TD-MS010), the Natural Science Foundation of Liaoning Province, China (No. 2022-KF-23-09), the Independent Research fund of Joint National-Local Engineering Research Centre for Safe and Precise Coal Mining(Anhui University of Science and Technology) (No. EC2022021).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

Authors Xiaoyu Zhang, Yingjie Liu, Tianbai Zhou, Yongbo Cai, were employed by the company Chinese Institute of Coal Science. The remaining author Bin Zhang declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Qi, Q.; Li, Y.; Zhao, S.; Zhang, N.; Zheng, W.; Li, H.; Li, H. Seventy Years Development of Coal Mine Rockburst in China: Establishment and Consideration of Theory and Technology System. Coal Sci. Technol. 2019, 47, 1–40. [Google Scholar]
  2. Pan, Y. Rockburst in Coal Mine; Science Press: Beijing, China, 2018. [Google Scholar]
  3. Sun, J.; Cheng, J. Study on the Perception and Alarm Method of Coal Mine Rock Burst and Coal and Gas Outburst. Ind. Mine Autom. 2022, 48, 1–6. [Google Scholar]
  4. Pan, Y. Integrated Study on Compound Dynamic Disaster of Coal-gas Outburst and Rockburst. J. China Coal Soc. 2016, 41, 105–112. [Google Scholar]
  5. Yu, M. Linear and Nonlinear Unified Strength Theory. Chin. J. Rock Mech. Eng. 2007, 26, 662–669. [Google Scholar]
  6. Petukhov, M.; Linkov, M. The Theory of Post-Failure Deformations and the Problem of Stability in Rock Mechanics. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1979, 16, 57–76. [Google Scholar] [CrossRef]
  7. Calder, P.N.; Madsen, D. High Frequency Precursor Analysis Prior to a Rockburst. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1989, 26, 3–4. [Google Scholar]
  8. Singh, S.P. Burst Energy Release Index. Rock Mech. Rock Eng. 1988, 21, 149–155. [Google Scholar] [CrossRef]
  9. Li, Y. Rockburst Mechanism and its Preliminary Application. J. China Univ. Min. Technol. 1985, 14, 37–43. [Google Scholar]
  10. Zhang, M. Instability Theory and Mathematical Model for Coal/Rock Bursts. Chin. J. Rock Mech. Eng. 1987, 6, 197–204. [Google Scholar]
  11. Li, H.; Qi, Q.; Zhao, S.; Li, H.; Shu, H.; Chen, L. Discussion on Generalized “Three Factors” Mechanism of Coal Mine Dynamic Disaster. Coal Sci. Technol. 2021, 49, 42–52. [Google Scholar]
  12. Qi, Q.; Pan, Y.; Li, H.; Jiang, D.; Shu, L.; Zhao, S.; Zhang, Y.; Pan, J.; Li, H.; Pan, P. Theoretical Basis and Key Technology of Prevention and Control of Coal-Rock Dynamic Disasters in Deep Coal Mining. J. China Coal Soc. 2020, 45, 1567–1584. [Google Scholar]
  13. Li, X.; Lin, B. Status of Research and Analysis on Coal and Gas Outburst Mechanism. Coal Geol. Explor. 2010, 38, 7–13. [Google Scholar]
  14. Ma, Y.; Nie, B.; He, X.; Li, X.; Meng, J.; Song, D. Mechanism Investigation on Coal and Gas Outburst: An Overview. Int. J. Miner. Metall. Mater. 2020, 27, 872–887. [Google Scholar] [CrossRef]
  15. Jiang, C.; Yu, Q. Spherical Shell Losing Stability” Hypothesis about Outburst Mechanism and Prevention Technology; China University of Mining and Technology Press: Xuzhou, China, 1998. [Google Scholar]
  16. Liang, B.; Zhang, M.; Pan, Y.; Wang, L. Theory of Instability of Flow Fixation Coupling for Coal and Gas Outburst. J. China Coal Soc. 1995, 20, 492–496. [Google Scholar]
  17. Wang, E.; Zhang, G.; Zhang, C.; Li, Z. Research Progress and Prospect on Theory and Technology for Coal and Gas Outburst Control and Protection in China. J. China Coal Soc. 2022, 47, 297–322. [Google Scholar]
  18. Hou, Q.; Li, H.; Fan, J.; Ju, Y.; Wang, T.; Li, X.; Wu, Y. Structure and Coalbed Methane Occurrence in Tectonically Deformed Coals. Sci. China Earth Sci. 2012, 55, 1755–1763. [Google Scholar] [CrossRef]
  19. Hou, Q.; Luo, Y.; Song, C.; Han, Y.; Du, J.; Xu, R. Gas Generation during Middle-Rank Coal Deformation and the Preliminary Discussion of the Mechanism. J. China Coal Soc. 2014, 39, 1675–1682. [Google Scholar]
  20. Wang, J.; Guo, G.J.; Han, Y.; Hou, Q.; Geng, M.; Zhang, Z. Mechanolysis Mechanisms of the Fused Aromatic Rings of Anthracite Coal under Shear Stress. Fuel 2019, 253, 1247–1255. [Google Scholar] [CrossRef]
  21. Wang, J.; Hou, Q.; Zeng, F.; Guo, G.-J. Stress Sensitivity for the Occurrence of Coalbed Gas Outbursts: A Reactive Force Field Molecular Dynamics Study. Energy Fuels 2021, 35, 5801–5807. [Google Scholar] [CrossRef]
  22. Alder, B.J.; Wainwright, T.E. Phase Transition for a Hard Sphere System. J. Chem. Phys. 1957, 27, 1208–1209. [Google Scholar] [CrossRef]
  23. Dai, X.; Bai, J.; Yuan, P.; Du, S.; Li, D.; Wen, X.; Li, W. The Application of Molecular Simulation in Ash Chemistry of Coal. Chin. J. Chem. Eng. 2020, 28, 2723–2732. [Google Scholar] [CrossRef]
  24. Jiang, J.; Babarao, R.; Hu, Z. Molecular Simulations for Energy, Environmental and Pharmaceutical Applications of Nanoporous Materials: From Zeolites, Metal-Organic Frameworks to Protein Crystals. Chem. Soc. Rev. 2011, 40, 3599–3612. [Google Scholar] [CrossRef] [PubMed]
  25. Yang, Y.; Liu, J.; Yao, J.; Kou, J.; Li, Z.; Wu, T.; Zhang, K.; Zhang, L.; Sun, H. Adsorption Behaviors of Shale Oil in Kerogen Slit by Molecular Simulation. Chem. Eng. J. 2020, 387, 124054. [Google Scholar] [CrossRef]
  26. Ungerer, P.; Nieto-Draghi, C.; Lachet, V.; Wender, A.; di Lella, A.; Boutin, A.; Rousseau, B.; Fuchs, A.H. Molecular Simulation Applied to Fluid Properties in the Oil and Gas Industry. Mol. Simul. 2007, 33, 287–304. [Google Scholar] [CrossRef]
  27. Fuchs, W.; Sandhoff, A.G. Theory of Coal Pyrolysis. Ind. Eng. Chem. 1942, 34, 567–571. [Google Scholar] [CrossRef]
  28. Xiang, J.; Zeng, F.; Liang, H.; Li, B.; Song, X. Molecular Simulation of the CH4/CO2/H2O Adsorption onto the Molecular Structure of Coal. Sci. China Earth Sci. 2014, 57, 1749–1759. [Google Scholar] [CrossRef]
  29. Wang, K.; Zhang, B.; Kang, T. Effect of Geological Depths on CH4 Adsorption, Diffusion, and Swelling in Kaolinite by Molecular Simulations. Energy Fuels 2020, 34, 1620–1626. [Google Scholar] [CrossRef]
  30. Zhang, B.; Zeng, F.; Wang, D.; Kang, G.; Zhang, X.; Kang, T. Macromolecular Structure Model of Anthracite in Southern Qinshui Basin and its Methane Bearing Mechanical Properties. J. China Coal 2021, 46, 534–543. [Google Scholar]
  31. Taheri, Z.; Pour, A.N. Studying of the Adsorption and Diffusion Behaviors of Methane on Graphene Oxide by Molecular Dynamics Simulation. J. Mol. Model. 2021, 27, 1–8. [Google Scholar] [CrossRef]
  32. Zhang, X.; Kang, T.; Zhang, B.; Kang, G.; Kang, J.; Zhang, R.; Zhang, L. Combined Molecular Simulations and Experimental Study of Methane Adsorption on Anthracite. Energy Fuels 2020, 34, 12118–12125. [Google Scholar] [CrossRef]
  33. Ma, R.; Zhang, S.; Hou, D.; Liu, W.; Yuan, L.; Liu, X. Model Construction and Optimization of Molecule Structure of High-Rank Coal in Feng County, Shaanxi Province. J. China Coal Soc. 2019, 44, 1827–1835. [Google Scholar]
  34. Zhang, L.; Kang, T.; Kang, J.; Zhang, X.; Zhang, B.; Guo, J.; Chai, Z. Response of Molecular Structures and Methane Adsorption Behaviors in Coals Subjected to Cyclical Microwave Exposure. ACS Omega 2021, 6, 31566–31577. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, X.; Zhang, R.; Kang, T.; Hu, Y.; Li, C. Experimental and Mechanistic Research on Methane Adsorption in Anthracite Modified by Electrochemical Treatment using Selected Electrode Materials. Sci. Rep. 2019, 9, 1–12. [Google Scholar] [CrossRef] [PubMed]
  36. Materials Studio. 2020. Available online: https://cgzb.tyut.edu.cn/info/1078/4523.htm (accessed on 16 December 2020).
  37. Zhang, B.; Kang, J.; Kang, T. Effect of Water on Methane Adsorption on the Kaolinite (0 0 1) Surface based on Molecular Simulations. Appl. Surf. Sci. 2018, 439, 792–800. [Google Scholar] [CrossRef]
  38. Skipper, N.T.; Refson, K.; McConnell, J.D.C. Computer Calculation of Water-clay Interactions using Atomic Pair Potentials. Clay Miner. 1989, 24, 411–425. [Google Scholar] [CrossRef]
  39. Cygan, R.T.; Liang, J.J.; Kalinichev, A.G. Molecular Models of Hydroxide, Oxyhydroxide, and Clay Phases and the Development of a General Force Field. J. Phys. Chem. B 2004, 108, 1255–1266. [Google Scholar] [CrossRef]
  40. Heinz, H.; Koerner, H.; Anderson, K.L.; Vaia, R.A.; Farmer, B.L. Force Field for Mica-Type Silicates and Dynamics of Octadecylammonium Chains Grafted to Montmorillonite. Chem. Mater. 2005, 17, 5658–5669. [Google Scholar] [CrossRef]
  41. Xiao, S.; Edwards, S.A.; Gräter, F. A New Transferable Forcefield for Simulating the Mechanics of CaCO3 Crystals. J. Phys. Chem. C 2011, 115, 20067–20075. [Google Scholar] [CrossRef]
  42. Abramov, A.; Iglauer, S. Application of the CLAYFF and the DREIDING Force Fields for Modeling of Alkylated Quartz Surfaces. Langmuir 2019, 35, 5746–5752. [Google Scholar] [CrossRef]
  43. Lammers, L.N.; Bourg, I.C.; Okumura, M.; Kolluri, K.; Sposito, G.; Machida, M. Molecular Dynamics Simulations of Cesium Adsorption on Illite Nanoparticles. J. Colloid Interface Sci. 2017, 490, 608–620. [Google Scholar] [CrossRef]
  44. Pouvreau, M.; Greathouse, J.A.; Cygan, R.T.; Kalinichev, A.G. Structure of Hydrated Kaolinite Edge Surfaces: Dft Results and Further Development of the Clayff Classical Force Field with Metal–O–H Angle Bending Terms. J. Phys. Chem. C 2019, 123, 11628–11638. [Google Scholar] [CrossRef]
  45. Meng, J.; Niu, J.; Meng, H.; Xia, J.; Zhong, R. Insight on Adsorption Mechanism of Coal Molecules at Different Ranks. Fuel 2020, 267, 117234. [Google Scholar] [CrossRef]
  46. Zhu, H.; Zhang, Y.; Qu, B.; Liao, Q.; Wang, H.; Gao, R. Thermodynamic Characteristics of Methane Adsorption about Coking Coal Molecular with Different Sulfur Components: Considering the Influence of Moisture Contents. J. Nat. Gas Sci. Eng. 2021, 94, 104053. [Google Scholar] [CrossRef]
  47. Cheng, Y.; Hu, B. Main Occurrence Form of Methane in Coal: Micropore filling. J. China Coal Soc. 2021, 46, 2933–2948. [Google Scholar]
  48. Tang, X.; Wang, Z.; Ripepi, N.; Kang, B.; Yue, G. Adsorption Affinity of Different Types of Coal: Mean Isosteric Heat of Adsorption. Energy Fuels 2015, 29, 3609–3615. [Google Scholar] [CrossRef]
  49. Zhang, B.; Kang, J.; Kang, T. Molecular Simulation of Methane Adsorption and its Effect on Kaolinite Swelling as Functions of Pressure and Temperature. Mol. Simul. 2018, 44, 789–796. [Google Scholar] [CrossRef]
  50. Feng, Z.; Cai, T.; Zhou, D.; Zhao, Y.-S.; Wang, C. Temperature and Deformation Changes in Anthracite Coal after Methane Adsorption. Fuel 2017, 192, 27–34. [Google Scholar] [CrossRef]
  51. Chen, J.; Min, F.; Liu, L.; Liu, C. Mechanism research on surface hydration of kaolinite, insights from DFT and MD simulations. Appl. Surf. Sci. 2019, 476, 6–15. [Google Scholar] [CrossRef]
  52. Zhou, X.; Zhang, X.; Xu, S.; Wu, S.; Liu, Q.; Fan, Z. Evaluation of Thermo-Mechanical Properties of Graphene/Carbon-Nanotubes Modified Asphalt with Molecular Simulation. Mol. Simul. 2017, 43, 312–319. [Google Scholar] [CrossRef]
  53. Guo, H.; Yu, Y.; Wang, K.; Yang, Z.; Wang, L.; Xu, C. Kinetic characteristics of desorption and diffusion in raw coal and tectonic coal and their influence on coal and gas outburst. Fuel 2023, 343, 127883. [Google Scholar] [CrossRef]
Figure 1. Schematic of coal–rock compound structure models.
Figure 1. Schematic of coal–rock compound structure models.
Applsci 13 12975 g001
Figure 2. Methane adsorption isotherm of anthracite (a) and coking coal (b) at temperatures in the range of 273.15–313.15 k, and Langmuir fitting constant, a, of anthracite (c) and coking coal (d).
Figure 2. Methane adsorption isotherm of anthracite (a) and coking coal (b) at temperatures in the range of 273.15–313.15 k, and Langmuir fitting constant, a, of anthracite (c) and coking coal (d).
Applsci 13 12975 g002
Figure 3. Adsorption heat of methane at different temperatures in anthracite and coking coal models. (a) Adsorption heat of methane in anthracite. (b) Adsorption heat of methane in coking coal.
Figure 3. Adsorption heat of methane at different temperatures in anthracite and coking coal models. (a) Adsorption heat of methane in anthracite. (b) Adsorption heat of methane in coking coal.
Applsci 13 12975 g003
Figure 4. Methane adsorption isotherm of the anthracite–mineral model (a) and coking coal–mineral model (b) at a temperature of 273.15 k, and Langmuir fitting constant, a, of the anthracite–mineral model (c) and coking coal–mineral model (d).
Figure 4. Methane adsorption isotherm of the anthracite–mineral model (a) and coking coal–mineral model (b) at a temperature of 273.15 k, and Langmuir fitting constant, a, of the anthracite–mineral model (c) and coking coal–mineral model (d).
Applsci 13 12975 g004
Figure 5. Adsorption heat of methane at 273.15 K in anthracite–mineral model (a) and coking coal–mineral model (b).
Figure 5. Adsorption heat of methane at 273.15 K in anthracite–mineral model (a) and coking coal–mineral model (b).
Applsci 13 12975 g005
Figure 6. Change in mechanical properties after the addition of calcite and kaolinite.
Figure 6. Change in mechanical properties after the addition of calcite and kaolinite.
Applsci 13 12975 g006
Table 1. Simulation Schemes.
Table 1. Simulation Schemes.
SchemeMolecular ModelsAdsorption Pressure (MPa)Adsorption Temperature (K)
1Anthracite coal0.1, 0.3, 0.5, 1, 2, 3...10273.15, 283.15, 293.15, 303.15, 313.15
2Anthracite + 24.86 wt% Calcite 0.1, 0.3, 0.5, 1, 2, 3...10273.15
3Anthracite + 39.82 wt% Calcite 0.1, 0.3, 0.5, 1, 2, 3...10273.15
4Anthracite + 24.18 wt% Kaolinite0.1, 0.3, 0.5, 1, 2, 3...10273.15
5Anthracite + 38.94 wt% Kaolinite0.1, 0.3, 0.5, 1, 2, 3...10273.15
6Coking coal0.1, 0.3, 0.5, 1, 2, 3...10273.15, 283.15, 293.15, 303.15, 313.15
7Coking coal + 25.08 wt% Calcite 0.1, 0.3, 0.5, 1, 2, 3...10273.15
8Coking coal + 40.11 wt% Calcite 0.1, 0.3, 0.5, 1, 2, 3...10273.15
9Coking coal + 24.40 wt% Kaolinite0.1, 0.3, 0.5, 1, 2, 3...10273.15
10Coking coal + 39.17 wt% Kaolinite0.1, 0.3, 0.5, 1, 2, 3...10273.15
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, X.; Liu, Y.; Zhou, T.; Cai, Y.; Zhang, B. Molecular Simulation Study on the Effect of Co-Associated Minerals on Methane Adsorption and Mechanical Properties of Coal. Appl. Sci. 2023, 13, 12975. https://doi.org/10.3390/app132412975

AMA Style

Zhang X, Liu Y, Zhou T, Cai Y, Zhang B. Molecular Simulation Study on the Effect of Co-Associated Minerals on Methane Adsorption and Mechanical Properties of Coal. Applied Sciences. 2023; 13(24):12975. https://doi.org/10.3390/app132412975

Chicago/Turabian Style

Zhang, Xiaoyu, Yingjie Liu, Tianbai Zhou, Yongbo Cai, and Bin Zhang. 2023. "Molecular Simulation Study on the Effect of Co-Associated Minerals on Methane Adsorption and Mechanical Properties of Coal" Applied Sciences 13, no. 24: 12975. https://doi.org/10.3390/app132412975

APA Style

Zhang, X., Liu, Y., Zhou, T., Cai, Y., & Zhang, B. (2023). Molecular Simulation Study on the Effect of Co-Associated Minerals on Methane Adsorption and Mechanical Properties of Coal. Applied Sciences, 13(24), 12975. https://doi.org/10.3390/app132412975

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