Investigation of Competitive Adsorption Properties of CO/CO2/O2 onto the Kailuan Coals by Molecular Simulation

To reveal the CO, CO2, and O2 adsorption properties of two bituminous coals at different pressures and temperatures, the molecular unit-cell structures of two types of bituminous coal are constructed (C1180H960O120N20 and C1160H860O80N20) by Fourier transform infrared (FTIR) spectroscopy. The bituminous coal molecular FTIR spectroscopic curve is calculated by quantum chemistry, and the results are consistent with the experimental curve. The isothermal adsorption curves of the single-component gases CO, CO2, and O2 conform to the Langmuir equation from 20 to 60 °C. The adsorption simulations are mainly performed using grand canonical Monte Carlo (GCMC) methods. The amount of adsorption decreases with increasing temperature at the same pressure, and CO2 can be the first to reach adsorption saturation at the same temperature. The CO2/CO adsorption selectivity for binary gas mixtures has apparent advantages in low-pressure or shallow buried coal seams. The adsorption selectivity of O2/CO varying under different pressures is not obvious. The high amount of CO inhibits the adsorption capacity of CO2 and O2. In other words, the effect of injecting CO2 to control fire extinguishing in bituminous coal seams with high abnormal CO concentrations is not significant.


INTRODUCTION
Bituminous coal is the most widely distributed coal in nature. 1 It is mainly distributed in northern provinces (autonomous regions) in China, of which the bituminous coal reserves in North China account for more than 60% of the national reserves. Coal seams are prone to oxidative spontaneous combustion during mining and operations, which increases the risk of fire accidents. It is necessary to study the microstructure of coal to explore its adsorption mechanism and ignition mechanism. 2−4 The differences in the coal structure affect the permeation and adsorption of coalbed methane (CBM) in coal seams. 5 Chen et al., 6 Chen et al., 7 and Solomon et al. 8 used the analytical technique of Fourier transform infrared (FTIR) spectroscopy to determine the molecular differences in coal. The ability of coal with different molecular structures to adsorb gases differs. 9,10 O 2 is the dominant factor in the violent reaction stage of spontaneous coal combustion. 11−13 However, O 2 is physically adsorbed on the surface of coal functional groups before participating in the coal oxygen reaction. To prevent O 2 from reaching saturation adsorption and participating in chemical reactions, Fang et al. 14 designed an experiment to displace O 2 with inert gas. Li 15 and Chen et al. 16 used experiments to determine the physical adsorption of O 2 by lignite. They concluded that the oxygen uptake of coals with a low metamorphic and large specific surface area has higher oxygen absorption. Lu et al. 17 and Cheng et al. 18 found that O 2 was mainly adsorbed by van der Waals forces and compared the adsorption states of CO 2 , O 2 , and other small molecules in lignite. CO 2 is the greenhouse gas generated by the coal oxidation reaction and a common fire extinguishing component. 19 Wu et al. 20 used the grand canonical Monte Carlo (GCMC) method to simulate and compare CO 2 and O 2 . CO 2 is adsorbed by van der Waals forces and electrostatic action, so the competitiveness is usually CO 2 > O 2 . Most scholars generally compare the competition and adsorption laws between CO 2 and CH 4 . Zhou et al., 21,22 Gao et al., 23 Wang et al., 24 Zhang et al., 25 Ding, 26 and Sui et al. 27,28 discussed the competitive adsorption behavior of CO 2 and CH 4 in lignite or organic matter. However, for some mines with abnormal CO sources, the change in the CO index gas should be given more attention. Therefore, Zhu et al. 29 and Zhang et al. 30 used experimental instruments to analyze the performance of coal adsorbing CO. Deng et al. 31,32 simulated the explosion process and the explosion concentration limit of CH 4 and CO mixtures based on experiments. Zhou 33 used quantum chemistry to calculate the adsorption characteristics of CO and CO 2 on the coal surface, clarified the competitive adsorption process of CO mixed with other gases, and found that the affinity sequence of adsorption is CO 2 > CO.
The source of underground CO is oxidation, spontaneous combustion, and the original CO in coal, 34−36 which scholars have verified. The original CO in coal formation will also make CO exceed the standard of security. However, coal mine safety regulations stipulate that the maximum allowable value of the underground CO concentration is 0.0024%. 37 It is also common to use CO as an index gas to determine coal seam spontaneous combustion in the actual production process. 38,39 Therefore, it is necessary to study the adsorption of CO in coal. Although Zhang et al. 40 analyzed the competition between CO and other small molecules in lignite by molecular simulation, there is a lack of comparison of adsorption selectivity between CO and other small molecules.
Therefore, we aimed to clarify the competitive characteristics and adsorption capacity between O 2 , CO 2 , and CO of bituminous coal, aiming at the occurrence and diffusion behavior of CO in coal seams. The molecular structure parameters were calculated by Fourier infrared spectroscopy experiments, and a simplified molecular model and supercell structure of two types of bituminous coal were constructed. The adsorption capacity of the single-component gases CO, CO 2 , and O 2 in the adsorption system under different temperatures and burial depths was compared. The difference in gas adsorbed by different bituminous coals was analyzed. The binary adsorption competition relationship between CO, CO 2 , and O 2 was simulated to obtain the adsorption selectivity of CO, CO 2 , and O 2 at different concentrations and compare the adsorption selectivity at different pressures. This study provides a theoretical basis for mine CO anomalies and coal spontaneous combustion fire prevention.  41 The coal sample particle size was controlled below 200 mesh by crushing, screening, and grinding. Samples were placed into sealed bags and labeled. A Vario El III organic element tester and an automatic sulfur tester were used to measure the proportion of major elements in the coal samples. 42 A large gap resulted when the 5E-MAG6700 automatic industrial analyzer was used in the industrial analysis because of the influence of experimental conditions. Therefore, based on the industrial analysis method of coal (GB/T 212-2008), naturally dried coal with no additional moisture removal was taken as the reference coal, and the particle size was controlled to below 0.2 mm.

TEST AND
According to the national standard for coal classification in China (GB5751-86), no. 1 and 2 coal samples are gas coal and coking coal, both of which are bituminous coals. The moisture, ash, and various element compositions of the two coal samples are shown in Table 1. Under high-temperature conditions, small molecular side chains and active oxygen-containing functional groups in the coal molecular structure produce substances, such as H 2 O and CO 2 . Therefore, when the degree of coalification increases, polycondensation of the molecular structure of coal reduces the number of decomposition products that are formed during the thermal reaction. The volatile content of no. 1 coal is 41.43% and that of no. 2 coal is 28.39%. No. 1 coal has a low coalification degree and a high volatile content, and the percentage of fixed carbon and elemental carbon is less than that of no. 2 coal. Because of the high rank, the coal density increases, the pores are poorly developed, the specific surface area decreases, the carbon condensation level increases, the number of functional groups of high metamorphic coal relative to the low metamorphic coal decreases, and the waterabsorption capacity is inferior to that of the low metamorphic coal.
2.1.2. Fourier Transform Infrared Spectroscopy Experiment. Coal samples with a particle diameter below 200 mesh were mixed completely with KBr. The halogenating agent tablet pressing method was selected for the infrared spectrum test on a Shimadzu FTIR-8400 FTIR spectrometer made in Japan. The main parameters were set as follows: the image resolution was 4.0 cm −1 , the wavenumber was varied between 400 and 4000 cm −1 , each coal sample was scanned 30 times, the two coals were analyzed by FTIR, and the absorbance curve that corresponds to each wavenumber was obtained. 43 The functional group region in the infrared spectrum is 1300−4000 cm −1 , and the wavenumber in the fingerprint region is 650−1300 cm −1 . The infrared fingerprint region is like a human fingerprint and represents the characteristic peaks of some functional groups to distinguish subtle differences in the material structure. Figure 1 shows the trend chart of the peak Nomenclature: M ad represents moisture on an air dry basis; A d represents ash on a dry basis; V daf represents volatile on a dry ash free basis; FC d represents fixed carbon; S td stands for total sulfur on a dry basis; and O daf , C daf , H daf , and N daf represent the percentages of oxygen, carbon, hydrogen, and nitrogen elements, respectively. spectra of the coal samples after baseline calibration. The metamorphic degree of the coal sample affects the shape and position of the absorption peak. In the fingerprint area, the absorption peak is 700−900 cm −1 , which represents the change in the aromatic ring. The absorbance tends to increase with an increase in the coalification degree. The region 1000−1300 cm −1 represents ether bond destruction and recombination. In the functional group area, the peak pattern that corresponds to the same wavenumber in the curve in Figure 1 shows that the aliphatic side chain, the oxygen-containing functional group, and the hydrogen bond changed, and differences existed in the structure of the two bituminous coals.
The segmented FTIR spectra of the coal are compared and calculated quantitatively by peak fitting so that the peak group with a wavenumber of 400−4000 cm −1 is divided into four regions. Figure 2a shows aromatic-ring-substituted hydrogen in the range of 700−900 cm −1 . The characteristic functional groups of 1000−1800 cm −1 are shown in Figure 2b. Figure 2c shows the aliphatic hydrocarbon absorption zone of 2800−3000 cm −1 . Figure 2d shows hydroxyl or hydrogen bond absorption peak groups of 3000−3700 cm −1 . 44 As shown in Figure 2a, the peak number of two bituminous coals at 700−900 cm −1 is 3, and the absorbance of no. 2 coal increased significantly at ∼750 cm −1 . In Figure 2b, the absorbance of no. 2 coal at ∼1025 cm −1 is greater than that of no. 1 coal. Continuous peaks occurred at ∼1185 and 1250 cm −1 , and the peak intensity was weak. Shoulder peaks occurred at ∼1370 cm −1 , and two peaks existed near 1440 and 1600 cm −1 . After a weak peak appears near 1740 cm −1 , the absorbance of the two coals gradually approached 0. In Figure 2c, the absorbance curves of the two coal samples are roughly parallel without intersection. At the same wavenumber, the absorbance of the higher metamorphic coal is less than that of the lower metamorphic coal. Two broad peaks formed at ∼2851 and 2918 cm −1 . At 2851 cm −1 , the peak of no. 2 coal is lower than that of no. 1 coal, which indicates that fewer methylene groups are present. At 2918−3000 cm −1 , the slope of the absorbance curve of no. 2 coal is less than that of no. 1 coal, which indicates that a large methyl content exists in no. 2 coal. In Figure 2d, at ∼3034 cm −1 , the peak shape of no. 2 coal is sharper than that of no. 1 coal. At ∼3435 cm −1 , the peak neck of no. 1 coal is longer than that of no. 2 coal and the peak shoulder is wider than that of no. 2 coal. Accounting for the high metamorphic degree of no. 2 coal, the number of fused rings increases because of carbon condensation during coalification, which reduces the spatial distance of the hydroxyl groups and continuously formed selfassociated hydrogen bonds. The weak peak at ∼3518 cm −1 belongs to the hydrogen bond that is composed of a hydroxyl and a π bond.
In Figures 3−5, (a) represents no. 1 coal and (b) represents no. 2 coal. The baseline of the peak fitting is consistent, the fitting degree R 2 > 99.6%, and all peaks are Gaussian.  Four aromatic ring-substitution types existed in the experimental coal samples, namely mono-, di-, tri-, and pentasubstitution, and the corresponding wavenumbers are shown in Table 2. The change in the coal aromatic structure is shown in Figure 3. The single substitution contents of no. 1 and 2 coal are 1.3 and 5.8%, respectively. The main mode of no. 1 coal is metadi-substitution, with a content of 81.6%. No. 2 coal is mainly mono-, di-, and trisubstituted, with a content of 43.1%. No peak sample existed in the wavenumber range of 710−750 cm −1 for no. 1 coal. The algebraic quantity of mono-, tri-, and pentasubstituted no. 1 coal is 12.6% more than that of no. 2 coal, and the pentasubstituted content of no. 2 coal is 31.1% more than that of no. 1 coal. Therefore, a high rank of bituminous coal yields more substitution sites with a more stable structure and a greater ring-forming rate. The substitution site content changes from double to tri-and penta-substitution because of the disconnection of C−H bonds of the aromatic ring and substitution by some atoms or groups, dehydrogenation of naphthenic aliphatic hydrocarbons into rings, dehydrogenation of aromatic ring branches, and dehydroxylation and decarboxylation of the benzene structure. 45 The characteristic atoms or atomic groups in coal include ether bonds, carbonyls, carboxyls, hydroxyls, esters, and anhydrides. The wavenumber range of the main characteristic functional groups is 1000−1800 cm −1 , 46 including symmetric and antisymmetric stretching vibrations, such as C−O, CO, and CC bonds. 47 Figure 4 shows the number of oxygen- The peak spectra of small molecular aliphatic hydrocarbon groups of bituminous coals with two metamorphic degrees are shown in Figure 5, which are divided into eight peaks in the 2800−3000 cm −1 band for quantitative analysis. 48 The wavenumber is methyl antisymmetric stretching vibration near 2940−2975 cm −1 and methyl symmetric stretching vibration near 2865 ± 5 cm −1 . The wavenumber of 2880−2940 cm −1 represents the methylene antisymmetric stretching vibration, and the methylene symmetric stretching vibration occurs at ∼2845 ± 10 cm −1 . No. 1 coal showed the methyl group vibration at peaks of ∼2943, 2962, and 2864 cm −1 , and no. 2 coal also    The molecular composition of coal is centered on an aromatic nucleus, and many structural units with similar but different structures are connected by bridge bonds. There are also small molecular compounds. 49 The marginal atomic groups of the basic structural units of coal include carboxyl, phenolic hydroxyl, carbonyl, and methoxy groups with oxygen atoms, and alkyl branched side chains. The length of the molecular branch side chain is closely related to the rank of coal metamorphism. A deeper metamorphism yields a shorter side chain length and a smaller proportion of aliphatic group carbons and total carbons. When the carbon content is ∼70%, the alkyl carbon accounts for 8% of the total carbon, and the number of carbon atoms in the side chain is ∼2−3. When the carbon content is ∼80%, alkyl carbon accounts for 6% of the total carbon with ∼2.2 carbon atoms. For an elemental carbon content of ∼84%, the number of alkyl carbon atoms is ∼1.8. The main body of the coal structural unit is often expressed by parameters such as the number of condensed aromatic rings, the ratio of hydrogen to carbon atoms, the aromatic carbon rate, and the aromatic hydrogen rate. The bituminous coal with an 80% carbon content has two aromatic rings. A higher carbon content results in more rings. There are three aromatic rings in 85% coal. When the carbon content exceeds 90%, the coal structure may tend to graphitization or the number of closed rings exceeds 40. In addition to aromatic rings, nitrogen or sulfur atomic heterocycles may appear in the main structure. Bridge bonds occur at the connection between the molecular unit structures. The types of bridge bonds include −CH 2 −methylene, −O−, −S− iso-ether bonds, sulfide bonds, −CH 2 −O−, −CH 2 −S− iso-methylene ether bonds, and methylene−sulfide bonds. The distribution quantity and position of bridge bonds in molecules differ, and the content is uneven. Low-rank bituminous coal is dominated by −CH 2 − or −CH 2 −O− groups, whereas medium metamorphic bituminous coal is dominated by −O− or −CH 2 − in small quantities. As a vulnerable group of molecules, bridge bonds are prone to thermal or oxidative fracture, so bridge bonds reflect the stability of the coal molecular structure.
The parameters of the unit body are constructed through experimental data. , where C al represents the amount of aliphatic carbon, C represents the total carbon, H al is the amount of aliphatic hydrogen, H is the total hydrogen, and H al /C al is the number ratio of aliphatic hydrogen to carbon atoms, which is generally 1.8. 50,51 The f ar of no. 1 coal is 0.68 and that of no. 2 coal is 0.78.
(4) According to the percentage content of elements, the number of atoms that constitutes the coal structure is obtained, as shown in Table 3. To facilitate calculation and simulation, the total number of carbon atoms of coal molecules is set to 60, and the basic structures of two bituminous coals are constructed. Therefore, theoretically, the molecular formulae of no. 1 and 2 coal are C 60 H 47 O 7 N and C 60 H 40 O 4 N.
The coal molecular difference shown in Table 4 is the absolute value of the difference between the simulated and experimental values, which shows the error in building the model. If sulfur atoms appear in the unit model, the error between the element content of coal molecules and the experimental value is large. To enhance the authenticity of the simulated structure, through infrared analysis and structural parameters, after continuous adjustment, optimization, and modification on the basis of previous work, the structural configuration plan of the no. 1 coal molecule is obtained, as shown in Figure 7a, and no. 2 coal is shown in Figure 7b  2.3. Simulation Verification of the Molecular Structure. The COMPASSII force field is a molecular force field to unify the force field of organic molecular systems and that of inorganic molecular systems, which can be used to simulate  organic and inorganic small molecules, polymers, some metal ions, etc. The force field parameters come from the empirical parameters of quantum mechanics calculated from ab initio. The COMPASSII force field can analyze and calculate the molecular structure, vibration frequency, conformational energy, crystal structure, and binding energy density of the system and analyze and predict isolated and condensed molecules' structural energy characteristics. 52 The structure at the lowest point of molecular structure energy has research significance. First, the structure of the coal molecular simplified model is optimized, the calculation is run using Materials Studio (MS) software, and the structure at the local energy minimum is optimized using force tools. The parameter settings are as follows: the force field is uniformly set to COMPASSII, and the maximum number of iterations is 5000. The first molecular dynamics simulation was carried out. The parameters are set as follows: the ensemble is NVT, the temperature is 600.00 K, the control method is Nose, and the number of steps is 50000. To obtain a stable configuration, annealing dynamics simulation is required. The parameters are set as follows: annual, number of annual cycles is 5, initial temperature is 300.00 K, and mid-cycle temperature is 600.00 K. Finally, the second dynamic calculation is carried out, and the parameters are consistent with the first one. The calculation formula of molecular potential energy is given as follows 53 where E V is the valence electron energy, kcal/mol; E N is the nonbonding energy, kcal/mol; E B is the key stretching energy, kcal/mol; E A is the bond angle energy, kcal/mol; E T is the  torsional energy, kcal/mol; E I is the inversion energy, kcal/mol; E VAN is the van der Waals energy, kcal/mol; E E is the Coulombic energy; and E H is the hydrogen bond energy, kcal/mol. After calculation, the molecular spatial structure of no. 1 and no. 2 coal is shown in Figure 6. The gray, white, red, and blue balls are carbon, hydrogen, oxygen, and nitrogen. The spatial structure of the coal molecules deformed and twisted at the bridge bond position. In no. 1 coal, the pyrrole structure is perpendicular to the aromatic ring in space, and a large angle torsion existed in the oxygen-containing functional group. An approximate vertical state exists between the aromatic ring of the no. 2 coal and the adjacent aromatic ring, the rotation range at the aliphatic group is large, and the spatial voids become more obvious. According to formulae (1), (2), and (3), the total molecular energy decreases significantly from 845.8 kcal/mol in the first dynamic simulation to 766.2 kcal/mol in the second simulation, indicating that the molecular structure tends to be stable, and E B decreases from 93.4 to 78.8 kcal/mol. The stretching range of the bond becomes narrower, E A decreases from 140.7 to 106.4 kcal/mol, E VAN decreases from 42.7 to 37.0 kcal/mol, the bond angle is relatively contracted, and the molecular spacing becomes smaller.
The structure with the lowest energy is imported into the VAMP calculator in MS, the AM1 semiempirical Hamiltonian function in neglect of a diatomic differential overlap (NDDO) is used for geometric optimization, the calculation attribute is set to frequency, and the simulated infrared spectrum is obtained. In Figure 7, A is the infrared spectrum curve of the simulated molecules and B is the experimental infrared spectrum. The peak shape of the simulation results is consistent with the experimental results with a slight difference in the peak position. Some simulated peak positions are close to the experimental peak positions. For example, the peak group at the aromatic ringsubstitution position of no. 1 coal is caused by the bending vibration outside the −CH− bond plane. As another example, the peak group at 1000−1800 cm −1 of the no. 2 coal spectrum that is caused by stretching or bending vibration of oxygencontaining functional groups simulated that the peak position deviation is small. The position of absorption peak A in no. 1 coal is slightly higher than that of B, with a quantity of ∼170 units. The peak position after translation is consistent with the experimental spectrum. B has a wide peak at ∼3400 cm −1 , and the peak shape of the calculated spectrum is sharp. The characteristic frequency of the group in the infrared spectrum shows a self-associated hydroxyl hydrogen bond region, which exhibits mainly an intermolecular force. A shows the infrared spectrum of only one molecule, and the hydrogen bond in the molecule is weak. Therefore, the peak is shown in the infrared spectrum, and the peak width accounts for a relatively narrow proportion of the total spectrum width. The situation of no. 1 and 2 coal at ∼3400 cm −1 is basically the same, and a single molecule results in the simulation with a peak at the same position. The simulated peaks of no. 2 coal at 1050, 1400, and 1600 cm −1 are basically consistent with the experimental peaks. Excluding the influence of error and intermolecular force, the simulation results are satisfactory, which indicates that the structural unit is reasonable. 54

Adsorption Simulation Design.
After the simplified molecular model is optimized, the amorphous cell module in MS adds periodic boundary conditions to the simplified model. 21 The structure is optimized until the total energy converges. The task is construction, the quality is fine, the density is 1.400 g/ cm 3 , force field is COMPASSII, and the electrostatics is Ewald. After a series of geometric optimization and dynamic calculations, the system with the lowest energy is selected for subsequent calculation, the unit-cell size of the no.  20 ), as shown in Figure 6. It should be noted that the construction of periodic coal molecules does not represent the actual real molecules but only the statistical results of the distribution of main functional groups in coal. The simulated infrared spectrum shows that the chemical bond and composition in the model are consistent with the actual situation. Therefore, coal molecules with periodic boundary conditions can intuitively show the adsorption state of coal, which makes the simulation study easier to understand. 55,56 Different temperatures and pressures corresponding to the geological burial depth are selected. The surface temperature is set at 20°C; the temperature gradients are set at 30,40,50 25 The pressure is calculated at different burial depths under five temperature gradients of CO, CO 2 , and O 2 single gases. Then, the absolute adsorption capacity of mixed gas with CO ratios of 0.1, 0.3, 0.5, 0.7, and 0.9 is calculated. Isothermal adsorption curves are formed at different temperatures, and the adsorption differences of bituminous coal with different metamorphic degrees are compared. The competitive behavior of two bituminous coals under different pressures is judged by the adsorption competitive selectivity of CO and other gases.
The gas isothermal adsorption curve is calculated by the grand canonical Monte Carlo simulation method and realized by Materials Studio software sorption tools. 57 Task: fixed pressure, quality: ultrafine, and force field: COMPASSII. The conversion between fugacity and pressure is converted by the Peng Robinson formula, 23 in which fugacity is directly proportional to pressure.
Through the calculation of competitive adsorption energy by the Sorption module in MS software, the adsorption position and state of the adsorbate gas molecules in the model can also be observed, and the adsorption characteristics can be further clarified. Specific simulation parameter settings are as follows: task: locate, Monte Carlo method: Metropolis, force field: COMPASSII, and charges: force field assigned. The temperature was set to 293.15 K and the pressure to 8 MPa.

RESULTS AND DISCUSSION
3.1. Single-Component Gas Adsorption of CO, O 2 , and CO 2 . The single-component gas adsorption curves of the no. 1 coal and no. 2 coal are shown in Figure 8. The adsorption curves of CO, CO 2 , and O 2 of no. 1 coal are expressed by (a 1 ), (b 1 ), and (c 1 ), respectively. The adsorption curves of no. 2 coal are (a 2 ), (b 2 ), and (c 2 ). The fitting results of the adsorption curves in Figure 8 comply with the Langmuir equation, 58 and the Langmuir fitting formula is given as follows: where a represents the saturated adsorption capacity at infinite pressure, mmol/g; b is the adsorption constant, MPa −1 ; and R 2 is the fitting degree. The closer the value of R 2 is to 1, the closer the fitting degree is to the real value. The parameter values are shown in Table 5.
Under the same pressure, the adsorption capacity of both coal samples decreases with increasing temperature. The no. 1 coal generally has more significant adsorption capacity for CO and O 2 than the no. 2 coal under identical temperature and pressure. Because the rank of no. 1 coal metamorphism is low, the molecular structure is loose, many functional diagrams easily participate in the oxidation reaction, and the occurrence probability of CO and O 2 near the short-chain alkane structure becomes greater. CO is the product of the coal oxygen reaction and the oxide involved in the oxidation reaction. Compared with coal, which is not easy to oxidize, CO more easily exists in the structure of easily oxidized coal through molecular gaps. The adsorption between CO 2 and coal molecules occurs mainly through dispersion force, and an increase in the number of condensation rings leads to an increase in adsorption potential. Therefore, the adsorption capacity of CO 2 is stronger in coal molecules with a higher coalification degree. 59 Figure 9 shows the variation curve of the single-component gas adsorption capacity of two coal samples with pressure at the same temperature. The adsorption capacity is related to the physical properties of individual molecules and other factors. The no. 2 coal sample has the highest CO 2 adsorption capacity, which is greater than the CO and O 2 adsorption capacity. Higher critical temperature and pressure correspond to greater adsorption capacity. 60 The no. 1 coal has the largest CO 2 adsorption capacity in the range of 0−6.5 MPa. The adsorption law is that a larger molecular dynamics diameter corresponds to a smaller adsorption capacity. 24 The saturated adsorption pressure of CO 2 is approximately 8 MPa, which is significantly lower than O 2 and CO. The CO 2 and O 2 adsorption capacity cross point moves backward with increasing temperature. After the crossing point, the adsorption capacity of O 2 significantly increases, which may be caused by the smallest molecular dynamics diameter of CO 2 . O 2 continues to occupy the adsorption site after it first tends to saturate. A deeper coal seam experiences a greater pressure, which is not conducive to gas adsorption. When the burial depth is less than approximately 650 m, CO 2 injection should be used to prevent the coal seam oxygen concentration from being too high and prevent the physical adsorption of oxygen from further transforming into a coal oxygen chemical reaction. 61 Figure 9 shows that the adsorption capacity of CO 2 and O 2 is much greater than that of CO.
3.2. Binary Competitive Adsorption of CO with O 2 and CO 2 . In Figure 10, the solid line represents the adsorption curves of CO, O 2 , and CO 2 of no. 1 coal at 293.15 K, and the molar ratios are 0.1:0.9, 0.3:0.7, 0.5:0.5, 0.7:0.3, and 0.9:0.1. Figure 10 represents the adsorption curve of no. 2 coal. Figure 10 shows an apparent correlation between the number of small molecules adsorbed by the crystal cell and the molar ratio, and there is a general law in direct proportion. For no. 1 coal, the CO 2 adsorption capacity rapidly increases in the stage of 0−4 MPa, and the adsorption capacity tends to be flat after 8 MPa. The rapid growth stage of CO is at 0−8 MPa. Obviously, in the competitive adsorption of CO 2 and CO, the low-pressure stage can be close to the saturation state. When the molar ratio of CO is much greater than that of CO 2 , the adsorption capacity of CO and CO 2 is approximately 6 MPa. When the molar ratio of CO to CO 2 in the coal seam is less than 0.7:0.3, CO 2 injection can be used for fire prevention. 52 However, when the molar ratio of CO is approximately 9 times that of CO 2 , the effect of the CO 2 injection fire prevention technology is not significant. With increasing buried depth pressure, the effect is less obvious. In terms of the competitive adsorption capacity of different molar ratios of CO and O 2 , when the molar ratio of CO to O 2 is small or 1, the adsorption capacity of O 2 is higher than that of CO. The adsorption capacity is close at 0.7:0.3 and 0.9:0.1, and the adsorption capacity of CO is better than that of O 2 . As a polar molecule, CO will inhibit the adsorption of O 2 when CO occupies a high molar ratio. When the adsorption capacity of O 2 is higher than that of CO, air leakage should be controlled to prevent CO in coal seams containing primary CO from desorbing to improve the accuracy of predicting spontaneous coal combustion when CO is used as the index gas. 40 For no. 2 coal, compared with the CO 2 competitive adsorption capacity, the stronger adsorption capacity is shown by CO 2 . With increasing CO molar ratio, the CO adsorption capacity also increases. The decrease in CO 2 adsorption is more significant than CO adsorption. Changing the molar ratio significantly affects the adsorption capacity of CO 2 , so a high proportion of CO 2 will affect the adsorption capacity of CO. The CO 2 adsorption capacity of no. 2 coal at any molar ratio is stronger than that of no. 1 coal under identical conditions, which is consistent with the conclusion of single-gas adsorption. The change law of the competitive adsorption of CO and O 2 is consistent with that of the no. 1 coal. At 0.7:0.3, the isothermal adsorption curves of CO and O 2 at 0−4 MPa are close to coincidence, which indicates that in this case, the ability of CO and O 2 to occupy the adsorption site is equivalent. With increasing burial depth pressure, the CO adsorption capacity is slightly less than the O 2 adsorption capacity. Even if the oxygen concentration is small, the oxidation reaction of coal remains. 62 Therefore, when the spontaneous combustion of the coal  oxygen reaction occurs, oxygen is consumed for a short time, which rapidly decreases the oxygen concentration, and the CO concentration, which is the product of the coal oxygen reaction, increases. The CO adsorption index actively changes with the decrease in oxygen, so it can respond to the degree of the coal oxidation reaction.
3.3. Adsorption Selectivity of CO. Adsorption selectivity refers to the ability of adsorbents to preferentially adsorb some substances due to their different components and structures. The adsorption selectivity of binary mixtures of other gases G and CO is defined as where x G is the mole fraction of gas G in the adsorption component, x CO is the mole fraction of CO in the adsorption component, and y G and y CO are the mole fractions of G and CO in the free state. When S G/CO is greater than 1, in the binary mixed gas, G is preferentially adsorbed by the adsorbent compared with CO, and the adsorption capacity of the adsorbent for G is stronger. G is easily enriched in the adsorbent. 63 Figure 11a,c shows the difference in adsorption selectivity between CO, CO 2 , and O 2 at 293.15 K. The adsorption selectivity of CO 2 /CO decreases with increasing pressure, and the competitiveness of CO relative to CO 2 increases. When the CO concentration is 0.1, the changing trend of the ordinate corresponding to 0∼17 MPa is the largest, and the decrease range of adsorption selectivity at 0.5 CO is less than that of 0.1 CO. At 0.9 CO, the downward trend is gradually gentle, both coal samples have an evident gentle trend with the increase in the amount and concentration of CO, and the competitive selectivity of CO increases. However, S CO 2 /CO remains greater than 1, and CO 2 remains significantly more competitive than CO, which is consistent with the conclusion of their binary adsorption curves. When the pressure is less than 8 MPa, the  adsorption selectivity values of the two coal samples fluctuate in a certain range. After 8 MPa, the adsorption selectivity of CO 2 / CO is significantly inhibited, and it is evident at 0.1 CO. When the CO concentration is low, the competitiveness after 8 MPa is considerably more substantial than that before 8 MPa.
Coal with a high coal mineralization degree showed a greater adsorption selectivity for CO 2 /CO, and the adsorption selectivity increased when the amount of CO species increased. The selectivity for O 2 /CO shows a slightly increasing trend. All competitive adsorption capacities of O 2 are larger than those of CO at different species concentrations. 64 Adsorption selectivity for O 2 /CO: S O 2 /CO > S CO 2 /CO at pressures greater than 8 MPa in 0.1 CO.
3.4. Competitive Adsorption Energy of the Adsorbate. There are 10 adsorbate gas molecules, and the number of CO molecules is set to 1, 3, 5, 7, and 9. Correspondingly, the number of CO 2 is set to 9, 7, 5, 3, and 1. The number of O 2 is equal to that of CO 2 . Figure 12A shows the total energy change of CO and CO 2 competitive adsorption in the no. 1 coal. The adsorption process is an exothermic process. In the competition between CO and CO 2 , the CO amount is positively correlated with the total energy. The nonbonding energy accounts for the major part of the total energy, while the van der Waals energy accounts for the major part of the nonbonding energy. The van der Waals energy is consistent with the changing trend of the total energy and is positively correlated with the number of CO atoms. At 1CO + 9CO 2 , the nonbond energy is −71.34 kcal/mol, and the van der Waals energy is −59.78 kcal/mol. At 5CO + 5CO 2 , the van der Waals energy is −51.99 kcal/mol. At 9CO + 1CO 2 , the van der Waals energy is −47.77 kcal/mol. Therefore, most adsorption occurs due to van der Waals energy. The adsorption energies of CO and O 2 are also mainly dependent on the van der Waals energy. Although the curve in Figure 12B shows an overall downward trend, the difference between the lowest energy and the highest energy is only approximately 0.20 kcal/mol. The increase in the CO amount did not greatly change the total energy. Figure 13A shows the competitive adsorption energy of CO and CO 2 in the no. 2 coal, and the changing trend is basically identical to that of the no. 1 coal. The total energy significantly changes with the change in the amount of CO and CO 2 . In Figure 13B, the total energy of CO and O 2 competitive adsorption slightly fluctuates. Figures 12 and 13 clearly show that the angle, aggregation, and dispersion of the adsorbate change and the quantity and the concentration of the adsorbate also affect the variation in adsorption sites. The total energy can reflect the stability of the system's equilibrium state. All systems tend to decrease in energy. A greater absolute value of the total energy corresponds to a greater decrease in energy, i.e., adsorption will more likely occur in a more stable system. 20 Therefore, in the competition between CO and CO 2 , more CO 2 makes the system more stable. In other words, the increase in molarity of CO will increase the instability of the system. The total energy of the mixed system containing CO and O 2 has no significant difference, and the energy difference has no obvious rule with the change in the adsorbent molar concentration. However, the no. 2 coal has a larger total energy than the no. 1 coal, so the competitive adsorption of CO and O 2 may be greatly affected by the coal structure itself. Under an identical buried depth pressure, the increase in temperature will inhibit the adsorption capacity of a single-component gas, which shows a negative correlation between temperature and adsorption capacity. When the temperature is 20°C and the pressure is less than 6.5 MPa, the adsorption capacity is CO 2 > O 2 > CO. The pressure when CO 2 reaches saturated adsorption is smaller than that of CO and O 2 . (3) In binary competitive adsorption, the adsorption curves of 90% CO and 10% CO 2 are close, and the adsorption capacity of CO 2 is significantly inhibited with increasing CO concentration. A higher coalification degree of bituminous coal corresponds to a greater adsorption selectivity of CO 2 /CO. The adsorption capacity between 70% CO and 30% O 2 is close. When CO is higher than 70%, CO has greater adsorption competitiveness than O 2 . Therefore, when the CO content of the original coal seam is high or CO anomalies occur underground, it is not suitable to reduce the CO adsorption by injecting gas. For coal seam fire prevention, grouting or other methods should be changed to control the air leakage of the working face, prevent excessive O 2 from occupying the adsorption site, coal from further spontaneous combustion, and CO adsorbed by coal from diffusing and seeping into the working face.