Abstract
Knowledge of the gas sorption characteristics of a coal not only helps to explain the mechanism of enhanced coalbed methane recovery but also provides an important basis for simultaneous coal and gas extraction. In consequence, the pure and binary gas excess sorption capacity of methane, carbon dioxide, and nitrogen of bituminous coal samples derived from the Xutuan Coal Mine in Huaibei coalfield, in Anhui Province in China, was measured using the volumetric method. The fitting analysis of the pure gas Langmuir adsorption model was carried out. The binary gas excess sorption measurement showed that the final sorption capacity of bituminous samples was the same no matter what the gas adsorption order of competitive adsorption and displacement adsorption. Hence, coal gas adsorption is physical adsorption, i.e. the different adsorption and desorption process of gas molecules does not affect the final adsorption amount of coal to each component of gas. Using the fitting parameters obtained by the Langmuir equation, the extended Langmuir equation was used to predict the adsorption capacity for each component of the binary gas. The comparison between predicted adsorption capacity and measured adsorption capacity showed that the extended Langmuir equation can better describe the trend of the adsorption isotherm curves of a binary gas under different pressures. The separation coefficient and displacement coefficient were defined from Langmuir adsorption theory. The separation coefficient involves the proportion of each component in the free phase and the proportion of each component in the adsorption phase. The displacement coefficient involves the displacement ability of gas molecules at adsorption sites by free gas molecules.
Introduction
Coal is a highly porous material, and coalbed methane (CBM) is mainly adsorbed in the coal structure, hence the coal acts as both the source and as the reservoir. Natural porous media usually have an extremely complex pore structure with pore sizes that may extend over several orders of magnitude (Cai et al., 2012). As coal has many micropores, cleats, and fissures, the specific surface area of coal is generally 6–15 m2/g (Chen et al., 2013; Qi et al., 2013). As a consequence, coal exhibits a strong adsorption capacity for several gases (Yang et al., 2015) and adsorption gas accounts for more than 90% of the gas stored within the coal (Mastalerz et al., 2004). The goal of gas drainage is to release as much of the adsorbed gas as possible. In the primary coal seam, the gas is in a dynamic equilibrium state of adsorption, desorption, and migration. Reduction of the gas pressure in the coal seam leads to a reduction of its coal gas adsorption capacity.
It has long been recognized that the complex flow processes in porous media depend on the complexity of the microstructure of porous media (Cai et al., 2010). The adsorbed gas detaches from the micropore surfaces and is transmitted to the free gas phase. The free gas is transported first in the micropores of the coal matrix and then into the fissures between the matrix and then through the cleats finally into the drilled holes. Coal gas sorption studies involve an investigation of the gas storage and transport mechanism. Coal gas adsorption capacity is a key characteristic of CBM reserves. Enhanced coalbed methane (ECBM) recovery can be accomplished either by inert gas stripping or by displacement desorption methods (Ruthven, 1984; Yang, 2013) using nitrogen (N2), carbon dioxide (CO2), flue gas, compressor gas, and other types of industrial off-gas, which are used as a recovery agent (White et al., 2005). Hence, the process of ECBM recovery involves the adsorption and desorption characteristics of mixed gases. The effectiveness of the operation needs to be confirmed by multicomponent adsorption measurements, as ECBM depends on the relative affinity of the mixture of CO2 and CH4 in coal (Kurniawan et al., 2006; Pini et al., 2009). The binary gas sorption study undertaken during the present study provided a research basis for the study of ECBM migration behavior.
According to the trend of increased sorption capacity with the increasing pressure, the solid gas isotherm sorption curves generally are divided into six categories (Brunauer et al., 1940; Sing, 1998). Coal gas adsorption is physical adsorption, and its isothermal sorption curves are mainly Category I and Category II (Karacan and Okandan, 2000; Yang and Saunders, 1985). With regard to the adsorption model for gas in solids, Langmuir (1918) first proposed the single-layer molecular adsorption theory for adsorption sites of gas molecules and solid surfaces and established Langmuir adsorption isotherm equation. The simple form of the Langmuir equation and the clear meaning of the parameters allow it to be widely used in the assessment of CBM reserves. Some studies have shown that the Langmuir equation can be used with satisfactory precision to describe the coal gas sorption isotherm (George and Barakat, 2001; Lin et al., 2016; Wang et al., 2016; Weishauptová and Medek, 1998; Zhang et al., 2014).
When considering the multicomponent gas adsorption of a coal reservoir, the mixed gas adsorption capacity measurement of wet coal undertaken by Fitzgerald et al. (2005) has shown that moisture in the coal sample could cause large errors in the adsorption capacity of each component as predicted by the extended Langmuir (E-L) isotherm equation. Busch et al. (2006a) also found that coal moisture had a great impact on the adsorption capacity of each gas component. The gas adsorption study by Khosrokhavar et al. (2014) on organic carboniferous mudstone showed that the amount of excess carbon dioxide sorption was approximately seven times higher than the methane sorption capacity. Many scholars have found that the coal carbon dioxide sorption capacity was larger than that of methane and nitrogen (Krooss et al., 2002; Ottiger et al., 2008; Pan and Connell, 2007; Zhang et al., 2017, 2016). Clarkson and Bustin (2000) compared the CO2–CH4 binary gas sorption capacity on coals with different moisture contents and it was found that moisture had a larger impact on the carbon dioxide sorption capacity than that of methane. The adsorption experiments for CO2 and CH4 gas mixtures on coal showed that both CO2 and CH4 adsorption capacities performed with gas mixtures were lower than were those determined with pure gases. Thus, these results confirmed that both gases compete to be adsorbed by the coal. The comparative study of flue gas and pure carbon dioxide sorption properties of coal indicated that the flue gas adsorption capacity was much less than that of the carbon dioxide (Mazumder et al., 2006).
The bituminous coal of Xutuan Coal Mine in Huaibei coalfield in China was taken as research object for the present study. On the basis of previous work, a series of pure and binary gas adsorption tests were carried out using a bespoke Gas Adsorption and Strain Testing Apparatus (GASTA). All experiments were performed at 303 K on dry pulverized coal samples. The fitting analysis for a pure gas Langmuir adsorption model was carried out and the fitting parameters were obtained. Then, the binary gas competitive adsorption and displacement adsorption tests were carried out. Finally, based on the fitting parameters of the Langmuir equation for the pure gas, each component gas adsorption capacity was predicted by the E-L equation.
Measurement method
Research background
The Xutuan Coal Mine is located in Mengcheng County, southwest of Suzhou City, in the Anhui province of China. It is in Huaibei coalfield belonging to the Permian coal deposit. The relative gas emission of the mine is reported to be 22 m3/t, and the mine belongs to the high gas coal mine rank. The main coal seams of the mine are coal seam 32 and coal seam 82. The coal sample SA3 was derived from coal seam 32, and the coal sample SA8 was derived from coal seam 82. Samples with a particle size range of 0.25–0.18 mm were prepared in accordance with GB/T 474. Finally, the coal samples were dried completely under anaerobic conditions. The coal petrology parameters of the main coal seams are listed in Tables 1 and 2.
Maceral statistics of sampling coal seam.
The proximate analyses of sampling coal seam.
Note: “Mad” is the moisture content; “Aad” is the ash content, air-drying basis; “S.d” is the sulfur content; “P.d” is the phosphorus content; “Vdaf” is the volatile matter, dry ash-free basis.
Experiment equipment
The binary gas adsorption measurement was finished by the GASTA. A volumetric method is used to calculate the adsorptive capacity of coal samples. Inert gas is used to calibrate the free space volume of the sample chamber. The measured pressure in the chambers is used to calculate the adsorption capacity through the real gas state equation. All experiments were performed at 303 K (30°C) on dry pulverized coal samples.
The pressure range of the reference chamber and the sample chamber was 0–20 MPa. The internal diameter of the chambers was 60 mm and the inner height was 80 mm. The reference chamber, the sample chamber, and their auxiliary pressure transducers are core components of the test apparatus. The components of the apparatus are shown in Figure 1.

Schematic diagram of the GASTA.
Measurement and calculation method of adsorption capacity
In this investigation, a volumetric method was used to calculate the adsorption capacity. The pressure of the reference chamber and the sample chamber was measured by the pressure transducers. The adsorption capacity was calculated using the real gas state equation. Hence, the adsorption capacity refers to the excess (or Gibbs) adsorption capacity.
The ideal gas state equation is a gas state equation that describes the relationship between the pressure, volume, mass, and temperature. The real gas state is biased from the ideal gas state under different conditions, thus different real gas state equations are established. In this study, the gas state Soave–Redlich–Kwong (SRK) equation was used to calculate the amount of the gas in the fixed volume under a certain pressure. The SRK equation is shown as follows
The parameters of
For the sorption experiments of pure gas, the inert gas helium, which is not absorbed by coal, is injected into the reference chamber. The injected helium amount is calculated from the calibration volume of the reference chamber and measured pressure through the SRK equation. The valve between the reference chamber and sample chamber is opened after the gas pressure of the reference chamber and the sample chamber reaches equilibrium. The free volume of the sample chamber can be calculated from the injected helium amount and the equilibrium pressure. The same steps then were carried out for methane, carbon dioxide, and nitrogen and the adsorption capacity calculated using the SRK equation, with the fixed sample chamber free volume and measured pressures.
Methodology for binary gas isothermal adsorption measurements
For the binary gas adsorption experiment, the binary gas with a volume ratio of 1:1 first was prepared and introduced into the reference chamber. The next steps were the same with those for the pure gas. In order to obtain the competitive adsorption capacity of each gas into the coal samples, the composition of the exhaust gas after each competitive adsorption experiment was measured. This study used a gas sample bag to collect the exhaust gas and a gas composition analyzer was then used to measure the proportion of CO2 and CH4 in the exhaust gas.
The difference between the total amount of injected gas and the amount of free space gas was calculated using the actual gas SRK state equation. The SRK equation of state is used to calculate the molar volume of a pure gas at a certain pressure and temperature. When the gas mixture is calculated using the SRK equation, the equivalent transformation of Figure 2 can be used. A mixed gas with a fixed volume and pressure P is equivalent to having a nonresistant sliding partition around the inside of the container, which to the left side of the partition is CO2, to the right is CH4, and the pressure is P. The molar volume of CO2 and CH4 at this pressure can be calculated using the SRK equation.

Calculation of molar volume of mixed gas.
This investigation used the first adsorption equilibrium pressure
According to the SRK equation, the molar volume of carbon dioxide
From the equivalent transformation of Figure 2, the amount of CH4 injected into the cylinder system can be calculated as
The amount of free CO2 in the cylinder system was
The amount of free CH4 in the cylinder system was
And the amount of carbon dioxide and methane adsorbed per unit mass of coal was
In Formulae (2) to (7), the molar volume of the gas under standard conditions was 22.4 × 103cm3/mol, and the adsorption amount of CO2 and CH4 per unit mass of coal sample was
The method for calculating the adsorptive capacity of the displacement adsorption gas was approximately the same as that for competitive adsorption and therefore is not repeated here. As the displacement adsorption procedure first performs the adsorption test for pure gas, and then the adsorption test for the mixed gas is performed, the cylinder system must be evacuated after the adsorption equilibrium of the mixed gas.
The fitting analysis of the pure gas Langmuir adsorption model
Based on the hypothesis of surface adsorption sites, Langmuir (1918) derived a most popular solid/gas two-phase adsorption model. The Langmuir adsorption isotherm equation is usually expressed as the following formula
The pure gas experimental data for coal samples SA3 and SA8 were fitted by the Langmuir equation.
The fitting parameters of coal sample SA3 by Langmuir equation are shown in Table 3.
The Langmuir fitting parameters of coal sample SA3.
The fitting parameters of coal sample SA8 by Langmuir equation are shown in Table 4.
The Langmuir fitting parameters of coal sample SA8.
From the comparison between the experimental data above, it was found that the Langmuir equation can describe properly the isotherm sorption of CO2, CH4, and N2. From the goodness of fit for the data in Tables 3 and 4, the correlation coefficient of the Langmuir equation to the three kinds of gases was greater than 0.985, which meant that the sorption capacity prediction of the Langmuir equation achieved satisfactory results and the Langmuir equation gives significant guidance as to the potential gas reserves in the coal sample(s).
The Langmuir adsorption model assumes that the adsorption sites on the surface of the adsorbent are one to one corresponding to the adsorbate molecules. When the adsorption sites of all adsorbents are occupied by gas molecules, the maximum adsorption capacity of the adsorbent is reached, i.e. the Langmuir volume. According to the fitting data in Tables 3 and 4, coal has different maximum adsorption capacities for different gases. It was determined that the order of VL of the three gases was CO2 > CH4 > N2.
Adsorption analysis of binary gas
Coal matrix molecules show different adsorption affinity with different gas molecules, which means that coal will show different adsorption characteristics for different gases. The competitive adsorption and displacement adsorption of a binary gas are the main research content in this instance.
Competitive adsorption refers to the process by which the components of the binary gas are adsorbed at the same time. Displacement adsorption refers to the process whereby one gas first is adsorbed and then the other gas is injected to be absorbed by the coal sample. Although the essence of both experiments is binary gas adsorption, the former is an “on the same starting line” competition, and the latter is a “catch up from behind” scenario. In the present study, competitive adsorption and displacement adsorption tests with a 1:1 volume ratio binary gas were carried out on coal sample SA3.
Binary gas adsorption capacity
Comparative analysis of displacement adsorption and competitive adsorption
When the binary gas adsorption determinations had been completed, the binary gas adsorption isotherm curves were plotted (Figure 3).

The adsorption isotherms of binary CO2–CH4 gas.
As shown in Figure 3, the adsorption isotherm curves of competitive adsorption and displacement adsorption were basically overlapped, i.e. the binary gas adsorption capacity on the coal sample was the same. Gas adsorption for coal is a physical reaction of pore surfaces to the gas molecules. The sorption process will not change the pore surface characteristics of the coal. On the other hand, the equilibrium state of the coal gas adsorption is a dynamic equilibrium where the adsorption rate of gas molecules is equal to the desorption rate. In the process of displacement adsorption, methane was first injected into the closed chamber to achieve adsorption equilibrium. After injection of carbon dioxide into the chamber system, the adsorption equilibrium state of the methane gas molecules was broken. At this time, when the methane molecules were desorbed from the microporous surface of the coal, the carbon dioxide molecules, which have more affinity to porous surfaces, will occupy the adsorption sites. Eventually, the displacement adsorption condition will reach the same equilibrium state of competitive adsorption. In fact, the adsorption capacity of the coal sample to each component of the binary gas is related only to its partial pressure. Therefore, the competitive adsorption and displacement adsorption characteristics are only different for the sorption process of the gas molecules and do not affect the final adsorption capacity of coal to the mixture. Other scholars (Yang et al., 2015) have corroborated this conclusion.
Discussion contrasting pure and binary gas adsorption capacity
In order to analyze the contrasting characteristics of binary gas adsorption capacity, the pure gas and binary gas adsorption capacity are plotted in Figure 4.

Comparison of the characteristics of pure gas and binary gas adsorption. (a) Competitive adsorption and pure gas adsorption and (b) displacement adsorption and pure gas adsorption.
As shown in Figure 4, the shapes of the isothermal adsorption curves for pure gas and binary gas were basically the same, and both of them belong to the type I adsorption curve. Under the same adsorption equilibrium pressure, the total adsorption capacity of the binary gas was higher than that of pure methane and lower than that of pure carbon dioxide. For pure gas adsorption, the ratio of the adsorption capacity of the coal to the two gases was CO2:CH4 = 2.6:1. For the volume ratio 1:1 competitive adsorption, the ratio of the adsorption capacity of the coal to the two gases was CO2:CH4 = 3.3:1. For the volume ratio 1:1 displacement adsorption, the ratio of adsorption capacity of coal to two gases was CO2:CH4 = 3.31:1. The CO2 and CH4 adsorption capacity ratio difference of pure gas and binary gas adsorption indicates that carbon dioxide has a stronger adsorption advantage than methane in the binary gas. The adsorption and desorption rates of a gas molecule are involved only in the adsorption equilibrium of the single component gas. For binary gas adsorption, when the adsorbed gas molecules are separated from the adsorption site to form an empty adsorption site, the two different gas molecules compete for the empty adsorption sites. As carbon dioxide shows a higher adsorption affinity for coal, carbon dioxide has an advantage in competitive adsorption and when the binary gas eventually reaches the equilibrium state, carbon dioxide molecules occupy a higher proportion of the adsorption sites.
E-L equation and its prediction analysis
The Langmuir equation, which incorporates also the E-L equation, is the most widely used adsorption prediction model in the multicomponent adsorption field of CBM. Based on the fitting parameters for the pure gas adsorption isotherm curves, combined with the proportion of the binary gas in the state of adsorption equilibrium, the predicted adsorption capacity of coal samples to each component of the binary gas can be obtained from the E-L equation. In this study, the predicted binary gas adsorption capacity was obtained using the E-L equation and there was little difference between the predicted and the measured adsorption capacity.
E-L equation, separation coefficient, and displacement coefficient
Markham and Benton (1931) proposed a multicomponent adsorption model, i.e. the E-L equation based on the Langmuir adsorption theory of single component gas. The E-L equation preserves all the assumptions of the Langmuir adsorption theory for single component gases. It is considered that the adsorption of gas molecules on the surface of an adsorbent on mixed gases is competitive adsorption.
Combined with the fitted parameters of the Langmuir volume and the Langmuir pressure for the pure gas adsorption isotherm, the adsorption capacity of each component at the same temperature can be calculated using equation (7), and the total adsorption capacity is equal to the sum of the adsorption capacity of each component
The partial pressure of the mixed gas can be determined from the adsorption equilibrium pressure and the component volume fraction
The separation coefficient characterizes the relative adsorption ability of the two components in the mixture
The separation coefficient characterizes the relative adsorption ability of the two components of the binary gas. The separation coefficient is close to 1, indicating that the adsorption ability of the two components is similar. If the separation coefficient is greater than 1, it indicates that the adsorption ability of the component i to the adsorption site is stronger than that of the component j. If the reverse relationship holds good, then the adsorption ability of component j is stronger than that of component i.
From the E-L equation, the expression of the separation coefficient can be transformed into the following equation
As shown in equation (10), the separation coefficient that is derived from the E-L equation is a constant. The separation coefficient characterizes the proportion of the adsorption phase and the free phase in the different components, but it does not show the degree of competition between the two components at the adsorption sites. In order to characterize the degree of competition of the two components at the adsorption sites, that is the capability of one gas’s disturbance to the other gas in the binary gas sorption process, the displacement coefficient may be defined based on the single component adsorption
By comparing the equilibrium pressure with the behavior of a single component equal to the adsorption amount under the component pressure, the displacement coefficient defines the replacement capacity of the component in the mixture. When the displacement coefficient is greater than 1, compared with the pure gas adsorption, the component
With the Langmuir equation and E-L equation, the expression of the displacement coefficient can be transformed into the following equation
As shown in equation (12), the displacement coefficient is related to the Langmuir pressure and the partial pressure of the component and is unrelated to the Langmuir volume.
Prediction and analysis of binary gas adsorption by E-L equation
The E-L prediction of the adsorption capacity and the measured binary gas adsorption capacity are presented in Figure 5.

The measured values of binary gas adsorption and the predicted values of the E-L equation. (a) Competitive adsorption and E-L prediction and (b) displacement adsorption and E-L prediction.
The relative deviation of predicted adsorption capacity and measured adsorption capacity is shown in Tables 5 and 6.
The predicted and measured values for competitive adsorption of CO2–CH4.
The predicted and measured values for displacement adsorption of CO2–CH4.
As shown in Figure 5, the trend of the predicted and measured isothermal sorption curves is basically the same and the adsorption capacity of E-L equation prediction generally is a little greater than the measured value. From the data of Tables 5 and 6, it was found that although the relative error of some test points was more than 15%, that of the majority of the test points was less than 10%. The relative error for competitive adsorption was lower than the relative error for displacement adsorption. The predicted adsorption capacity of CO2 is always greater than the measured adsorption capacity. At lower equilibrium pressure, the predicted adsorption capacity of CH4 is less than the measured adsorption capacity but at higher equilibrium pressure, the predicted adsorption capacity of CH4 is greater than the measured adsorption capacity.
Generally speaking, the relative error of the E-L equation in the prediction adsorption capacity is within the control range. Hence, the E-L equation has a direct relationship to the adsorption isotherm curve plot for a mixed gas.
From the data of the Langmuir pressure and Langmuir volume in Table 3, the separation coefficient for CO2 to CH4 in coal sample SA3 can be calculated to be 3.43 using equation (10) and the measured separation coefficient is listed in Tables 7 and 8.
Measured separation coefficient under different equilibrium pressures for the competitive adsorption of CO2–CH4.
Measured separation coefficient under different equilibrium pressures for the displacement adsorption of CO2–CH4.
As shown in Tables 7 and 8, the measured separation coefficient under different equilibrium pressures was not constant. These data vary at around 3.43, which means that when the binary CO2–CH4 free gas ratio is 1:1 the adsorbed phase ratio is 3.43:1.
The displacement coefficient
Measured displacement coefficient under different equilibrium pressures for the competitive adsorption of CO2–CH4.
Measured displacement coefficient under different equilibrium pressures for the displacement adsorption of CO2–CH4.
Conclusions
The experimental study of gas sorption characteristics helps not only to explain the mechanism of ECBM but also provides an important basis for simultaneous coal and gas extraction. The bituminous coal of Xutuan Coal Mine in Huaibei coalfield in China was taken as the research object in the present study. A series of pure and binary gas adsorption tests were carried out, using the bespoke GASTA. The main conclusions were as follows:
By comparing the test results, it was confirmed that the Langmuir equation has a good fitting behavior for the adsorption processes of CO2, CH4, and N2, and the fitting degree was greater than 0.985. This confirmed the rationality and relevance of the Langmuir equation to the gas adsorption/desorption process. The competitive and displacement adsorption mechanisms of a binary gas involve the differences in their adsorption sequence. Their overall adsorption capacity generally remains unchanged. This phenomenon indicates that coal gas sorption is a physical adsorption phenomenon. Based on the fitting parameters of the pure gas Langmuir equation, binary gas adsorption capacity could be predicted using the E-L equation, which was consistent with the measured test results. The separation coefficient of the test points and for the Langmuir theory was calculated, respectively. The component ratio of the free phase and the adsorbed phase was considered. Furthermore, it was confirmed that the displacement coefficient was defined by the ability of a component molecule to displace another component molecule from the adsorption sites. The displacement coefficient of CO2 to CH4 was greater than 1, and hence, a portion of the adsorbed CH4 can be displaced by CO2.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is supported by National Natural Science Foundation of China (51704274, 51502339, 51704278); the Young Talent Lifting Project from CAST, Beijing Natural Science Foundation (8184082); the Open Projects of State Key Laboratory of Coal Resources and Safe Mining CUMT (SKLCRSM16KF10); the State Key Laboratory Cultivation Base for Gas Geology and Gas Control (Henan Polytechnic University) (WS2017B02); Key Laboratory of Coal-based CO2 Capture and Geological Storage, Jiangsu Province CUMT (no: 2016B09); and the Priority Academic Programme Development of Higher Education Institutions in Jiangsu Province.
