High-Pressure Adsorption of Pure Gases on Coals and Activated

Dec 22, 2011 - ash %, 6.6, 6.3, 14.3, 13.0, 4.7, 6.5 ... (22) used different adsorption models to study the volumetric effects of .... (40) Thus, the ...
0 downloads 0 Views 937KB Size
Article pubs.acs.org/EF

High-Pressure Adsorption of Pure Gases on Coals and Activated Carbon: Measurements and Modeling Sayeed A. Mohammad, Arunkumar Arumugam, Robert L. Robinson, Jr., and Khaled A. M. Gasem* School of Chemical Engineering, Oklahoma State University, Stillwater, Oklahoma 74078, United States S Supporting Information *

ABSTRACT: High-pressure gas adsorption of pure methane, nitrogen, and CO2 were measured on five coals of varying rank and an activated carbon. All measurements were conducted at 328.2 K and pressures to 13.8 MPa using a volumetric method. The adsorption isotherm data are for dry Beulah Zap, Illinois no. 6, Pocahontas, Upper Freeport, and Wyodak coals and Filtrasorb (F-400) activated carbon. At about 7 MPa pressure, the ratio of excess adsorptions for N2/CH4/CO2 varied from about 1/1.9/2.8 for Pocahontas coal to 1/2.1/5.9 for Beulah Zap coal. A distinct maximum in excess adsorption was observed for CO2 on each coal between 6 and 8 MPa. For activated carbon at 7 MPa, the ratio of excess adsorptions for N2/CH4/CO2 was 1/ 1.5/2.3, which is consistent with our earlier data on this activated carbon at a different temperature (318.2 K). The isotherm measurements have expected experimental uncertainties that vary from 1 to 7% for the adsorption of the three gases. The simplified local-density (SLD) model was fit to describe the pure-gas adsorption data. Overall, the weighted average absolute deviation (WAAD) was about 0.7 for both the coals and activated carbon. Generalized expressions were then developed to facilitate the prediction of nitrogen and CO2 adsorption isotherms, solely on the basis of the methane adsorption measurements and coal characterization information. In general, the newly developed generalized expressions were capable of predicting the nitrogen and CO2 adsorption within two times the experimental uncertainties for coals. The generalized expressions were further tested on an external data set comprised of three dry coals from the literature. For each coal, the CO2 adsorption isotherms were predicted on the basis of the methane isotherms. The generalized predictions provided average absolute percentage deviation (% AAD) of about 10−12% for CO2 adsorption on the three coals, which were about two times the deviations observed through direct regressions of the SLD model for these systems.

1. INTRODUCTION Geologic sequestration of CO2 is being studied internationally as a means to reduce the carbon footprint and limit the increase of anthropogenic greenhouse gases such as CO2 into the atmosphere. Although several options for CO2 sequestration are being considered, a potentially attractive avenue is the storage of CO2 in deep, unmineable coalbeds.1 Such coalbeds frequently contain large amounts of recoverable methane (or natural gas), and the recovery of this natural gas can be enhanced by injecting CO2 into the coalbeds. The injection of CO2 can serve dual purposesto provide an increased supply of natural gas (an important and less greenhouse-intensive fossil fuel) and to simultaneously sequester large amounts of CO2. The design of optimal recovery of natural gas and CO2 sequestration in coalbeds relies greatly on the availability of high-pressure, supercritical adsorption data for natural gas components as well as reliable adsorption models that are capable of accurate predictions of adsorption phenomena. Several investigators, including the authors, 2−8 have presented high-pressure, pure-gas adsorption isotherms on coals for coalbed gases (methane, nitrogen, and CO2) under conditions encountered in coalbed-related work. Mixture adsorption isotherms on coals have also been presented in some studies.5−7,9 In the current work, we present adsorption of pure methane, nitrogen, and CO2 on five coals covering a wide range of rank and an activated carbon, all under dry conditions. The coals are from the premium coal sample program at the Argonne National Laboratory,10 and the © 2011 American Chemical Society

activated carbon is a commercial carbon (F-400) from Calgon Carbon Company. The main objectives of this work are to (a) present new raw data for the adsorption of pure gases on five well-characterized coals and activated carbon and (b) develop a new generalized model that is capable of predicting the adsorption of other gases (such as nitrogen and CO2) based on adsorption measurements for methane on a given adsorbent. The puregas adsorption isotherms on coals presented here have appeared earlier in graphical format in the context of adsorption modeling studies8,11 or in an interlaboratory study on CO2 isotherms;12 however, the “raw” experimental data for these isotherms have not been reported previously. The current work also presents two generalized models that were developed on the basis of the adsorption data presented here. Such generalized models, if successful, can reduce significantly the number of isotherms that need to be measured experimentally. In fact, several studies present adsorption isotherms for either one or two gases, and therefore, a generalized model can be helpful in predicting isotherms for components that have not been measured. Further, a generalized adsorption model would be especially useful in coalbed reservoir simulations that require prediction of adsorption amounts under varied conditions. Received: September 15, 2011 Revised: November 28, 2011 Published: December 22, 2011 536

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Figure 1. Schematic diagram of the experimental apparatus.

The CO2 adsorption data presented here was measured as part of a Department of Energy/National Energy Technology Laboratory (NETL) interlaboratory study.12 That “round robin” study was organized to (a) assess the reproducibility of CO2/coal adsorption isotherm measurements among several laboratories and (b) develop a standard method for measuring coal adsorption isotherms. However, significant discrepancies were observed among the different laboratories in the reported adsorption amounts for CO2 in the NETL interlaboratory study. Goodman et al.12 attributed these discrepancies to the different procedures used to remove moisture in obtaining the dried coal samples. Some of the adsorption systems reported here have also been measured independently by other investigators at similar conditions. In particular, Sakurovs et al.13 measured adsorption of methane, nitrogen, and CO2 on three Argonne coals. In addition, Busch et al.5 reported adsorption isotherms for methane and CO2 on five Argonne coals; however, the temperature of their study (22 °C) was quite different than the one used in the present work (55 °C). As mentioned earlier, the main aim of the present work is to develop generalized correlations to facilitate predictions of gas adsorption on coals. The remainder of this paper is organized in the following manner: section 2 discusses the experimental method used in this work; section 3 presents the simplified local-density model; section 4 contains the experimental results; and section 5 presents details on the model generalization efforts undertaken as part of this work.

ments.6−8,14 The essential details of our measurement method are described in the following paragraph. The entire apparatus is maintained in a constant temperature air bath. The equilibrium cell (Figure 1) is filled with the adsorbent under study, and the cell is initially placed under vacuum. The void volume in the cell, Vvoid, is determined by injecting a known quantity of helium from a calibrated injection pump (Ruska). Because helium is not significantly adsorbed, the void volume can be determined from measured values of the temperature, pressure, and amount of helium injected into the cell. The material balance equation is

2. EXPERIMENTAL METHODS AND PROCEDURES

The Gibbs adsorption (also known as the excess adsorption) can be calculated directly from experimentally measured quantities. For puregas adsorption measurements, a known quantity, ninj, of gas (e.g., methane) is injected from the pump section into the cell section. Some Gibbs , will exist of the injected gas will be adsorbed and the remainder, nunads in the equilibrium bulk gas phase in the cell. A material balance is used

Vvoid =

( PZTΔV )pump ⎛ P2 P ⎞ ⎜ − 1⎟ Z1T ⎠ ⎝ Z 2T cell

(1)

where ΔV is the volume of the gas injected from the pump, Z is the compressibility factor of helium, T is the temperature, P is the pressure, subscripts “cell” and “pump” refer to conditions in the cell and pump sections of the apparatus, respectively, and 1 and 2 refer to conditions in the cell before and after injection of gas from the pump, respectively. The helium void volume was measured at the temperature of the gas adsorption isotherms (328.2 K) and over a range of pressures from atmospheric to about 13.8 MPa (2000 psia) in intervals of 1.4 MPa (200 psia). Several injections made into the cell at different pressures showed consistency in the calculated void volume. Generally, the void volume calculated from sequential injections varied less than 0.3 cm3 from the average value of approximately 80 cm3. The helium void volume includes all the volume of the cell section exclusive of the adsorbent volume that is impenetrable to helium gas. The constancy of the calculated void volume from the incremental injections over a range of pressures confirmed the validity of our assumption that adsorption of helium is negligible at the conditions of the measurements. The typical variation of void volume as a function of pressure is shown graphically in the Supporting Information.

2.1. Adsorption Isotherm Measurements. The experimental technique used in this study is based on the volumetric method of measuring adsorption. The experimental apparatus, shown schematically in Figure 1, has been used successfully in previous measure537

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Table 1. Analyses of Adsorbents Used in This Study analysesa

Beulah Zap

Wyodak

carbon % hydrogen % oxygen % sulfur %

72.9 4.83 20.3 0.80

75.0 5.35 18.0 0.63

moisture % vol. matter % fixed carbon % ash %

32.2 30.5 30.7 6.6

28.1 32.2 33.0 6.3

Illinois no. 6

Upper Freeport

Ultimate (Dry, Ash-free basis) 77.7 85.5 5.00 4.70 13.5 7.5 4.83 2.32 Proximate (As-Received Basis) 8.0 1.1 36.9 27. 1 40.9 58.7 14.3 13.0

Pocahontas

activated carbon (F-400)

91.1 4.44 2.5 0.66

94.8 0.79 3.2 0.78

0.7 18.5 76.1 4.7

1.0b 3.7 89.9 6.5

a Coal analyses provided by Argonne National Laboratory; activated carbon analyses provided by Hoffman Laboratories, Inc. bEquilibrium moisture content for activated carbon determined at Oklahoma State University to be 27%.

Gibbs to calculate the amount adsorbed, nads , as

Gibbs Gibbs nads = ninj − n unads

because the coal samples are not exposed to the atmosphere during drying. 2.4. Error Analysis. Instrument calibration tests were performed routinely during the course of the experiments. Generally, the calibrations were performed and/or repeated before the adsorption experiments on a new adsorbent sample. The thermocouples and RTDs were calibrated against a platinum reference RTD. Sensotec Super TJE pressure transducers (range: 0−13.8 MPa) were calibrated using helium as the working fluid against a Ruska dead weight tester with a calibration traceable to the National Institute of Standards and Technology. Detailed information on calibration procedures is available elsewhere.19 The uncertainties in the experimentally measured quantities after calibrations were estimated as follows: temperature, 0.1 K; pressure, 6.9 kPa; and injected gas volume, 0.02 cm3. A detailed error analysis was performed to estimate the uncertainty associated with each experimental data point by propagating the errors from the primary measurements of pressure, temperature and volume. The error analysis indicated that the average uncertainties for the gas adsorption measurements vary between 1 and 7% for activated carbon and coals. Additional information on the error analysis used in our measurements is available elsewhere.20 2.5. Coal Swelling. The adsorption of certain gases (such as CO2) on coals may cause swelling of the coal matrix. Some investigators believe that adsorption of CO2 can alter the porous coal structure and these changes, if left unaccounted for, can result in large errors in the calculation of CO2 adsorption on coals. Several researchers have attempted to account for potential swelling of coal by incorporating volumetric corrections to the adsorption isotherm equations. Ozdemir et al.21 and Dutta et al.22 used different adsorption models to study the volumetric effects of CO2 adsorption on coals. Romanov et al.23 have also attempted to interpret the volumetric changes in coals under CO2 pressure. Pan and Connell24 developed a theoretical model to describe adsorption-induced coal swelling by balancing the change in surface energy due to adsorption to the change in elastic energy of the coal matrix. Day et al.25 measured swelling of coals caused by CO2 and corrected their adsorption measurements to account for volumetric changes to the sample. These corrections involved adjusting the void volume to account for an increased volume of coal sample. In our study, the helium void volume was measured before and after each adsorption isotherm experiment. The constancy of the calculated void volume within its experimental uncertainty of 0.3% indicated that there was no irreversible change to the volume of the sample. This result is also supported by the findings of Day et al.,25 who found the coal swelling to be entirely reversible on release of CO2 gas pressure. Although Day et al.25 then applied a correction to the isotherm for coal swelling, we have used a constant void volume in our data reduction procedures. Thus, the CO2 adsorption data reported in this study are under the assumption that there is no appreciable swelling of the coal. The adsorption of methane and nitrogen is generally not expected to cause significant swelling of the coal matrix.

(2)

The amount injected can be determined from pressure, temperature, and volume measurements of the pump section:

ninj =

⎛ P ΔV ⎞ ⎜ ⎟ ⎝ ZRT ⎠pump

(3)

Similarly, the amount of unadsorbed gas is calculated from conditions at equilibrium in the cell Gibbs = n unads

⎛ PVvoid ⎞ ⎜ ⎟ ⎝ ZRT ⎠cell

(4)

where the pressure P is measured after equilibrium is reached in the cell (usually within 6 to 12 h, depending on the adsorption capacity of the adsorbent), which occurs when no further change in pressure is observed. In eqs 3 and 4, Z is the compressibility factor of the gas at the applicable conditions of temperature and pressure.

The above steps are repeated at sequentially higher pressures to yield a complete adsorption isotherm. In this study, all measurements of excess adsorption are reported in terms of mmol/g adsorbent on a dry basis. Further, all measurements reported in this work were made after complete drying of the adsorbent under vacuum in an equilibrium cell at 353 K for 36 h. 2.2. Gas Compressibility Factors. As is evident from the above equations, accurate compressibility factors are required for pure gases for proper adsorption data analysis. The compressibility factors for pure methane, nitrogen, and CO2 were calculated from highly accurate equations of state.15−17 Further, for void volume determination, the helium compressibility factor was calculated with an expression based on experimental data from the National Bureau of Standards Technical Note 631 for helium.18 2.3. Materials. The pure gases used in this work were obtained from Airgas (PA) with reported purity of 99.99% and were used as received. The Argonne coal samples were obtained from the Argonne National Laboratory (Argonne, IL) in ampules containing 5 g of 100mesh material of each coal. The activated carbon (F-400) was obtained from the Calgon Carbon Company. The compositional analyses of coals and activated carbon used in this study are presented in Table 1. The Illinois no. 6 coal is a high-volatile bituminous coal from the Illinois no. 6 or Herrin seam, whereas Wyodak coal is a subbituminous coal from the Wyodak−Anderson seam. The Upper Freeport coal is a medium-volatile bituminous coal, Pocahontas coal is a low-volatile bituminous coal, and Beulah Zap coal is lignite and is the lowest rank coal in this study.10 All isotherm measurements were conducted after complete drying of the coal. The coal samples were dried under vacuum at 80 °C (353 K) for 36 h to a constant weight. Further, the drying was conducted in the same equilibrium cell, which was used to measure adsorption isotherms. These steps prevent oxidation/degradation of the samples 538

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

where NA is Avogadro’s number, Ψ(z) and Ψ(L − z) are the fluid−solid interactions for the two surfaces of a slit of length L. Substituting eqs 6, 7, and 8 into eq 5 provides the equilibrium relationship for adsorption within the slit:

3. SIMPLIFIED LOCAL DENSITY (SLD) ADSORPTION MODEL In this study, we utilized the simplified local-density (SLD) model to describe the adsorption behavior of pure gases on coals and activated carbon. The SLD model accounts for the fluid−fluid and fluid−solid interactions in a slit-shaped pore. The model was first developed by Rangarajan et al.26 who used the van der Waals equation of state (EOS) to provide the fluid−fluid interaction information. Several researchers have used different equations of state, including the Peng−Robinson, Bender, and Elliot−Suresh−Donohue EOS to provide fluid− fluid interaction information.27−30 Following our earlier work with the SLD model,11,31,32 the Peng−Robinson EOS is used in this study. The essential details of the SLD model are outlined below. The SLD model envisions the adsorbent to be composed of a rectangular-shaped slit and the adsorbate molecules reside within this two-surface slit, as illustrated in Figure 2. The

⎛ Ψfs(z) + Ψfs(L − z) ⎞ ⎟⎟ fff (z) = fbulk exp⎜⎜ − kT ⎠ ⎝

where k is the Boltzmann’s constant. The SLD model is a simplification of the more general localdensity theory. The term “local” refers to the thermodynamic properties of a fluid at any local point z, where an average single density value is calculated, ρ(z).26 In addition, the SLD model uses mean-field theory in calculating the chemical potential. The mean-field theory replaces all interactions with an effective or average interaction so that no fluctuations are considered within the slit. Applying the SLD model, the excess adsorption (nEx) is given as

A Right Side of Slit (ρ(z) − ρ bulk ) dz (10) 2 Left Side of Slit Ex Here, n is the excess adsorption of adsorbate in number of moles per unit mass of adsorbent, and A is the surface area of the adsorbate on a particular solid. The lower limit in eq 10 is 3 /8σff, which is 3/8 of the diameter of an adsorbed molecule touching the left plane surface. The upper limit is L − 3/8σff, the location of an adsorbed molecule touching the right plane surface. The local density is assumed to be zero for the distances less than 3/8σff away from the wall. The value 3/8σff is chosen to account for most of the adsorbed gas; details are given elsewhere by Fitzgerald.33 The left and right sides of the slit each comprise half of the total surface area, A/2. Following previous studies at Oklahoma State University,11,32 the Peng−Robinson equation of state34 was used to provide the bulk fluid density and fugacity. The EOS, expressed in terms of density, is given as P 1 = ρRT (1 − ρb) a (T )ρ − RT[1 + (1 − 2 )ρb][1 + (1 + 2 )ρb] (11) n Ex =

Figure 2. SLD model slit geometry.

distance between the slit surfaces is L, and the position of a molecule within the slit is z. The position, z, is orthogonal to the solid surface formed by carbon atoms. A molecule within the slit has interactions with both walls of the adsorbent slit. At equilibrium, the chemical potential of the fluid, μ, is expressed as the sum of the fluid−fluid and fluid−solid potentials at a position, z, as follows:

μ(z) = μff (z) + μfs(z) = μ bulk (5) where the subscript “bulk” refers to the bulk fluid and “ff” and “fs” refer to the fluid−fluid and fluid−solid interactions, respectively. The chemical potential of the bulk fluid can be expressed in terms of fugacity as

⎛f ⎞ μ bulk = μ0(T ) + RT ln⎜⎜ bulk ⎟⎟ ⎝ f0 ⎠

a (T ) = (6)

b=

0.457535α(T )R2TC2 PC

(12)

0.077796RTC PC

(13)

The term, α(T), in eq 12 is calculated using the following expression developed at Oklahoma State University in an earlier work:35

(7)

2

α(T ) = exp((A + BTr)(1 − TrC + Dω+ E ω )) (14) where A, B, C, D, and E are correlation parameters and their values are 2.0, 0.8145, 0.134, 0.508, and −0.0467, respectively. The values used were based on an accurate description of saturation pressures for the pure gases under conditions encountered in coalbed operations. Other fluid properties used in this study are listed in Table 2.

where f ff(z) is fluid fugacity at a position z and f 0 refers to the same arbitrary reference state as in eq 6. The fluid−solid interactions in the model are accounted for through a potential energy function. In particular, the fluid− solid potential is given as

μfs(z) = NA[Ψfs(z) + Ψfs(L − z)]



where

where subscript “0” designates an arbitrary reference state and f refers to fugacity. Similarly, the chemical potential from fluid− fluid interactions is given as

⎛ f (z ) ⎞ ⎟⎟ μff (z) = μ0(T ) + RT ln⎜⎜ ff ⎝ f0 ⎠

(9)

(8) 539

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

fluid at high pressures. The covolume is corrected by an adjustable parameter, Λb

Table 2. Physical Properties of Fluids Used in SLD Model TC (K) PC (MPa) σff (nm) εff/k (K)

methane

nitrogen

CO2

190.56 4.599 0.3758 148.6

126.19 3.396 0.3798 71.4

304.13 7.377 0.3941 195.2

bads = b(1 + Λb) Equation 16 then becomes

f (z ) badsρ(z) ln ff = P 1 − badsρ(z)

The fugacity of the bulk fluid using PR EOS is

ln

fbulk P

=

bρ a (T )ρ − 1 − bρ RT (1 + 2bρ − b2ρ2) ⎡ P a (T ) Pb ⎤ − ln⎢ − ⎥− ⎣ RT ρ RT ⎦ 2 2 bRT ⎡ 1 + (1 + ln⎢ ⎣ 1 + (1 −

2 )ρb ⎤ ⎥ 2 )ρb ⎦

2 RT[1 + 2badsρ(z) − bads ρ(z)2 ] ⎡ 1 − badsρ(z) ⎤ aads(z) − ln⎢ ⎥− ⎣ RT ρ(z) ⎦ 2 2 badsRT ⎡ 1 + (1 + 2 )ρ(z)b ⎤ ads ln⎢ ⎥ ⎣ 1 + (1 − 2 )ρ(z)bads ⎦

(15)

(18)

In this work, we used a fixed value of Λb = 0.1 for coals and 0.3 for activated carbon. The fluid−solid interaction, Ψfs(z), was represented by Lee’s partially integrated 10−4 potential,36 which is a truncated form of Steele’s 10−4−3 potential:37

f (z ) bρ(z) ln ff = P 1 − bρ(z)

⎛ σ10 Ψfs(z) = 4πρatomsεfsσ2fs⎜⎜ fs 10 ⎝ 5(z′)

aads(z)ρ(z)

RT (1 + 2bρ(z) − b2ρ2(z)) ⎡ P aads(z) Pb ⎤ − ln⎢ − ⎥− ⎣ RT ρ(z) RT ⎦ 2 2 bRT ⎡ 1 + (1 + 2 )ρ(z)b ⎤ ln⎢ ⎥ ⎣ 1 + (1 − 2 )ρ(z)b ⎦

aads(z)ρ(z)



For the adsorbing fluid, the fugacity for fluid−fluid interactions is



(17)



4 ⎞ σ4fs 1 ⎟ ∑ 2 i = 1 (z′ + (i − 1)σss)4 ⎟⎠

(19)

where εfs is the fluid−solid interaction energy parameter and ρatoms = 0.382 atoms/ Å2 is the solid atom density. The parameters σff and σss signify, respectively, the molecular diameter of the adsorbate and the carbon interplanar distances. The carbon interplanar distance was taken to be the value for graphite, 0.335 nm38 and values of σff and εff were taken from Reid et al.39 The fluid−solid molecular diameter, σfs and dummy coordinate z′ used in numerical integration of eq 10 are defined as

(16)

The parameter aads(z) in eq 16 varies with the position within the slit. Chen et al.27 provided equations for aads(z), which depends on the ratio of slit length L to the molecular diameter σff. In our earlier work with the SLD model, the covolume b in the PR EOS was adjusted to improve the predictive capability for adsorption of pure gases on activated carbon and coals. The covolume has a significant effect on the local density of the adsorbed fluid near the surface. A simple empirical correction was used to account for the repulsive interactions of adsorbed

σ + σss σfs = ff 2

(20)

Table 3. Adsorption of Pure Gases on Dry Pocahontas Coal at 328.2 K methane

a

nitrogen

CO2

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorptiona (mmol/g)

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.71 1.46 2.85 4.22 5.61 7.01 8.38 9.73 11.12 12.48 13.80

0.311 0.451 0.601 0.679 0.729 0.756 0.779 0.796 0.802 0.808 0.810

0.022 0.021 0.021 0.021 0.021 0.021 0.022 0.023 0.023 0.025 0.026

0.76 1.46 2.84 4.23 5.61 6.99 8.36 9.72 11.13 12.51 13.79

0.102 0.167 0.255 0.318 0.360 0.397 0.422 0.441 0.460 0.475 0.483

0.018 0.018 0.018 0.018 0.018 0.018 0.019 0.019 0.020 0.021 0.022

1.07 1.66 2.85 4.25 5.55 6.97 8.33 9.84 11.10 12.38 13.84

0.715 0.831 0.979 1.067 1.110 1.124 1.089 1.001 0.860 0.758 0.695

0.037 0.036 0.036 0.035 0.035 0.035 0.037 0.043 0.066 0.077 0.105

σ is the estimated experimental uncertainty. 540

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Table 4. Adsorption of Pure Gases on Dry Upper Freeport Coal at 328.2 K methane

nitrogen

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.73 1.46 2.84 4.23 5.62 6.98 8.37 9.74 11.12 12.47 13.76

0.241 0.340 0.448 0.512 0.549 0.576 0.593 0.603 0.612 0.617 0.617

0.019 0.019 0.019 0.018 0.019 0.019 0.019 0.020 0.021 0.022 0.023

pressure (MPa) 0.77 1.44 2.83 4.25 5.61 7.03 8.41 9.74 11.11 12.50 13.77

CO2

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.079 0.124 0.193 0.237 0.273 0.300 0.319 0.335 0.349 0.360 0.371

0.016 0.016 0.016 0.016 0.016 0.016 0.017 0.017 0.018 0.019 0.020

0.71 1.52 2.88 4.26 5.60 6.96 8.33 9.67 10.98 12.44 13.82

0.485 0.646 0.789 0.866 0.906 0.919 0.904 0.851 0.751 0.655 0.608

0.030 0.029 0.029 0.028 0.028 0.028 0.030 0.033 0.049 0.058 0.077

Table 5. Adsorption of Pure Gases on Dry Illinois No. 6 Coal at 328.2 K methane

nitrogen

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.67 2.45 3.93 5.54 6.94 8.30 9.69 11.03 12.44 13.80

0.265 0.459 0.565 0.649 0.689 0.736 0.769 0.788 0.802 0.818

0.028 0.027 0.027 0.027 0.027 0.028 0.029 0.030 0.032 0.039

CO2

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.69 1.39 2.78 4.19 5.52 6.95 8.32 9.70 11.08 12.46 13.87

0.084 0.136 0.211 0.267 0.310 0.344 0.373 0.397 0.416 0.431 0.446

0.024 0.024 0.023 0.023 0.023 0.024 0.024 0.025 0.026 0.027 0.029

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.70 1.41 2.78 4.16 5.57 6.94 8.25 9.64 10.98 12.29 13.81

0.592 0.822 1.096 1.282 1.403 1.485 1.502 1.432 1.232 1.033 0.913

0.042 0.041 0.040 0.040 0.039 0.040 0.043 0.053 0.080 0.098 0.134

Table 6. Adsorption of Pure Gases on Dry Wyodak Coal at 328.2 K methane

nitrogen

CO2

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.69 1.42 2.80 4.19 5.57 6.95 8.34 9.70 11.10 12.46 13.71

0.242 0.355 0.481 0.565 0.624 0.671 0.701 0.721 0.740 0.770 0.782

0.027 0.027 0.026 0.026 0.026 0.027 0.027 0.028 0.029 0.040 0.047

0.73 1.45 2.83 4.21 5.58 6.97 8.33 9.70 11.09 12.47 13.72

0.091 0.147 0.221 0.272 0.313 0.344 0.373 0.394 0.416 0.431 0.445

0.024 0.023 0.023 0.023 0.023 0.023 0.024 0.025 0.025 0.027 0.028

0.69 1.45 2.78 4.18 5.60 6.98 8.34 9.65 10.99 12.33 13.82

0.791 1.078 1.388 1.602 1.749 1.820 1.827 1.735 1.587 1.430 1.324

0.046 0.045 0.044 0.043 0.043 0.043 0.046 0.052 0.082 0.098 0.135

σ z′ = z + ss 2

area, A; fluid−solid interaction energy, εfs/k; and the slit length, L.

(21)

4. EXPERIMENTAL RESULTS

For the adsorbed phase, the slit is divided into two halves and each half is subdivided into 50 intervals. The local density is then calculated by solving eqs 15 and 18 simultaneously for each interval. Once the local density is determined across the slit, the excess adsorption is calculated by integrating eq 10 numerically using Simpson’s rule.40 Thus, the SLD model contains the following three regressible parameters: surface

The raw experimental data for pure-gas adsorption on dry Pocahontas, Upper Freeport, Illinois no. 6, Wyodak, and Beulah Zap coals are presented in Tables 3−7, respectively. The adsorption data for dry activated carbon are presented in Table 8. The tables include the pressure, excess adsorption values, and the estimated uncertainties in excess adsorption for each data point for the three gases (methane, nitrogen, and CO2). For each isotherm, two replicate runs were 541

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Table 7. Adsorption of Pure Gases on Dry Beulah Zap Coal at 328.2 K methane

nitrogen

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.72 1.44 2.83 4.22 5.60 6.97 8.33 9.71 11.09 12.48 13.78

0.257 0.361 0.476 0.544 0.591 0.621 0.651 0.664 0.681 0.695 0.708

0.033 0.032 0.032 0.031 0.032 0.032 0.033 0.034 0.035 0.037 0.042

pressure (MPa) 0.77 1.46 2.83 4.22 5.60 6.98 8.34 9.73 11.09 12.47 13.80

CO2

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.089 0.137 0.201 0.244 0.277 0.297 0.316 0.331 0.339 0.346 0.358

0.028 0.028 0.027 0.027 0.027 0.028 0.028 0.029 0.030 0.032 0.033

0.61 1.42 2.78 4.17 5.54 6.96 8.30 9.59 11.01 12.22 13.81

0.819 1.108 1.397 1.591 1.706 1.760 1.752 1.661 1.470 1.344 1.248

0.048 0.047 0.046 0.045 0.045 0.045 0.047 0.054 0.086 0.101 0.141

Table 8. Adsorption of Pure Gases on Dry Activated Carbon (F-400) at 328.2 K methane

nitrogen

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

1.50 2.78 4.11 5.59 7.07 8.38 9.18 9.77 11.11 12.43 13.74

2.845 3.555 3.936 4.167 4.277 4.310 4.306 4.306 4.273 4.221 4.145

0.047 0.046 0.045 0.045 0.045 0.045 0.045 0.045 0.046 0.046 0.047

CO2

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

0.81 1.46 2.93 4.19 5.53 6.98 8.36 9.69 11.08 12.54 13.70

1.015 1.473 2.075 2.407 2.651 2.834 2.945 3.018 3.068 3.100 3.108

0.041 0.040 0.039 0.039 0.039 0.039 0.039 0.040 0.039 0.040 0.040

conducted to test for the reproducibility of the data. The two runs were well within the experimental uncertainties of the adsorption data. A more general and detailed analysis of experimental errors in adsorption measurements using a volumetric method have been reported elsewhere.20 As mentioned earlier, the CO2 adsorption data on dry coals were measured to complement an interlaboratory study conducted by NETL.12 The main objective of that study was to investigate the reproducibility of CO2/coal adsorption isotherm measurements among different laboratories. In contrast, the objective of our study is to develop gas/matrix calibrated model to facilitate generalized predictions of nitrogen and CO2 adsorption isotherms based solely on methane isotherm measurements and readily available coal characterization information (such as the ultimate and proximate analyses of coals). Such a capability, if successful, would reduce significantly the experimental burden of measuring gas adsorption isotherms on each coal. In other words, only methane isotherms would need to be measured on a given coal, and then, the generalized model may be used to predict, a priori, the adsorption of other two gases (nitrogen and CO2). Figures 3−7 illustrate the excess adsorption of pure gases on dry Pocahontas, Upper Freeport, Illinois no. 6, Wyodak, and Beulah Zap coals, respectively. These figures also depict the model representations obtained from the simplified local-density (SLD) model, as described in section 5. A distinct maximum is observed for CO2 adsorption isotherms on coals between 6 and 8 MPa. At about 7 MPa, the ratio of excess adsorptions for N2/CH4/CO2 ranged from about 1/1.9/2.8 for Pocahontas to 1/2.1/5.9 for Beulah Zap coal. The N2/CH4 ratio did not show large variations among the five coals; however, the CO2 excess adsorption maxima varied considerably among the five coals. The coals used in this work cover a wide range of rank. Specifically, the

pressure (MPa)

excess adsorption (mmol/g)

σ in excess adsorption (mmol/g)

1.49 2.85 4.22 5.62 7.07 8.31 9.62 11.11 12.49

4.887 5.885 6.321 6.462 6.396 6.134 5.616 4.524 3.522

0.113 0.110 0.108 0.107 0.106 0.105 0.106 0.170 0.174

Figure 3. Adsorption of pure gases on dry Pocahontas coal at 328.2 K: experimental data and SLD model representations. percent carbon of the coals ranged from about 72% to 91% on a dry, ash-free basis. In addition, the oxygen content of coals varied from about 2.5% to 20%. The oxygen content of coals is an indication of the presence of polar functional groups on the coal surface, and this also appears to affect the CO2 adsorption capacity on these coals to some extent. This aspect is discussed in more detail in section 5 on the SLD modeling results on coals. The adsorption of methane, nitrogen, and CO2 was also measured on a dry activated carbon (F-400). Figure 8 presents the pure-gas adsorption data on activated carbon. At 7 MPa, the ratio of excess 542

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Figure 4. Adsorption of pure gases on dry Upper Freeport coal at 328.2 K: experimental data and SLD model representations.

Figure 7. Adsorption of pure gases on dry Beulah Zap coal at 328.2 K: experimental data and SLD model representations.

Figure 5. Adsorption of pure gases on dry Illinois no. 6 coal at 328.2 K: experimental data and SLD model representations.

Figure 8. Adsorption of pure gases on dry activated carbon (F-400) at 328.2 K: experimental data and SLD model representations. and methane adsorption than for the adsorption of CO2. These trends are quite typical of the adsorbent systems considered in this work.

5. SLD MODELING The aim of this work was to develop a predictive, generalized model capable of predicting the adsorption of other gases (such as nitrogen and CO2) from the knowledge of methane adsorption isotherms for a given adsorbent. Four specific case studies were designed to guide this effort, and the generalized predictions were compared with direct regressions. The case studies 1−4 are described below: 1 Five parameters were regressed. They were surface area A, fluid−solid interaction energy of each gas εfs/k (total of three), and slit length L. The surface area and slit length were taken as common for all gases, and the fluid−solid interaction energy was the only gas-specific parameter regressed. 2 Three parameters were determined from the methane adsorption isotherm for each coal. The parameters were the surface area, slit length and fluid−solid interaction energy of methane. Then, retaining the surface area and slit length as fixed, only the fluid−solid interaction energies of nitrogen and CO2 were determined from their respective adsorption data.

Figure 6. Adsorption of pure gases on dry Wyodak Coal at 328.2 K: experimental data and SLD model representations. adsorptions for N2/CH4/CO2 was 1/1.5/2.3, which was consistent with our earlier data14 on this activated carbon at a different temperature (318.2 K). On average, the expected experimental uncertainties ranged from about 1 to 7% for the adsorption of three gases on activated carbons and coals. Lower uncertainties were obtained on activated carbon when compared to coals. Further, the uncertainties are generally smaller (on a percentage basis) for nitrogen 543

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Table 9. SLD Model Parameters for Pure-Gas Adsorption on Dry Coals and Activated Carbon case 1

case 2

εfs/k (K)

εfs/k (K)

coal

A (m2/g)

CH4

N2

CO2

Pocahontas Upper Freeport Illinois no. 6 Wyodak Beulah Zap

78.1 58.6 83.7 74.1 74.7

81.4 81.1 64.0 67.8 62.7

38.7 38.4 30.2 33.5 29.5

417.8

117.0

59.2 case 3

114.0 121.9 118.5 213.7 219.1 Activated 186.5

coal

A (m2/g)

CH4

N2

CO2

L (nm)

A (m2/g)

CH4

Pocahontas Upper Freeport Illinois no. 6 Wyodak Beulah Zap

85.3 63.5 79.2 74.5 76.7

73.0 73.5 66.8 68.0 61.5

35.6 35.9 32.1 32.8 29.1

73.0 73.5 66.8 68.0 61.5

105.2

53.8

1.46 1.56 2.87 2.64 2.17 Carbon 1.68

85.3 63.5 79.2 74.5 76.7

448.9

99.8 101.4 128.0 205.7 221.6 Activated 201.5

448.9

105.2

53.8

F-400

L (nm) 1.53 1.65 2.18 3.10 2.46 Carbon 1.87

A (m2/g)

CH4

N2

CO2

L (nm)

85.3 63.5 79.2 74.5 76.7

73.0 73.5 66.8 68.0 61.5

35.8 35.8 29.4 33.6 29.1

99.8 108.2 123.8 218.7 221.6

1.46 1.56 2.87 2.64 2.17

448.9

105.2

55.1 case 4

169.3

1.68

N2

CO2

L (nm)

35.6 35.9 32.1 32.8 29.1

99.7 108.2 128.8 190.7 221.6

1.46 1.56 2.87 2.64 2.17

192.9

1.68

εfs/k (K)

F-F00

εfs/k (K)

energy parameter for CO2 was much larger for two coals (Beulah Zap and Wyodak coals) that contained high levels of moisture and oxygen content. These aspects are investigated further in the two generalization cases presented in cases 3 and 4. Case 2. For case 2, the methane adsorption isotherm data alone was used to determine the common surface area, slit length and fluid−solid interaction energy of methane for each coal. Then, fluid−solid interaction energies of nitrogen and CO2 were determined from their respective adsorption data. Case 2 served as a precursor to the development of generalized models that are presented in cases 3 and 4. Tables 9 and 10 present the model parameters and statistics obtained for each coal and activated carbon for case 2. As shown in the tables, the statistics for case 2 were very similar to those for case 1, and they yielded an overall WAAD of 0.7 and 0.8 for the five coals and activated carbon. Generalization: Gas/Matrix Calibrated Models. The model parameters from case 2 were used to develop a predictive, generalized model to facilitate prediction of nitrogen and CO2 adsorption isotherms based on the methane adsorption isotherms. In many studies, the adsorption isotherms of all three gases (methane, nitrogen, and CO2) are not measured, and the adsorption data for only one (or two) gases are available. In such cases, it would be beneficial if the measured isotherms for one gas can be used to predict, a priori, the adsorption isotherms for other gases. Because the adsorption of methane appears to be the most frequently measured isotherms in coalbed-related work, we developed our generalization based on methane. Two alternative (and potentially equivalent) generalized models were developed, and they are discussed below in cases 3 and 4. The generalized models were developed only on the basis of adsorption data for the five coals measured in this work. Then, the adsorption data on activated carbon and an external data set comprised of three coals from the literature3 were used to test/validate the performance of the newly developed generalized correlations. In the following, we present the gas/matrix calibrated,

3 The fluid−solid interaction energies of nitrogen and CO2 were generalized in terms of fluid−solid interaction energy of methane and moisture content of the asreceived samples. The surface area, slit length, and the fluid−solid interaction energy for methane were used from case 2. 4 The fluid−solid interaction energy of CO2 was also generalized in terms of fluid−solid interaction energy for methane and the oxygen content of the coal. The generalized expression for nitrogen was the same as in case 3. Case 1 provides the ultimate correlative capability of the SLD model. Case 2 served as the precursor to the two generalized models that are presented in cases 3 and 4. The regressions used to determine the model coefficients in cases 1 and 2 minimized the weighted root-mean-square (WRMS) of the predicted amounts adsorbed, where the weights were the expected experimental uncertainties. Case 1. The parameters for each coal included the common surface area A, slit length L, and fluid−solid interaction energy parameter εfs/k for each of the three gases. Tables 9 and 10 present the model parameters and statistics for the cases, obtained for the five coals and activated carbon. Overall, the SLD model represented the adsorption data with a WAAD of 0.6 for five coals and 0.7 for activated carbon. Thus, the SLD model was capable of representing the adsorption data within the experimental uncertainties. Figures 3−8 illustrate the modeling results obtained for case 1 for pure-gas adsorption on coals and activated carbon. The εfs/k values were the highest for CO2 and lowest for nitrogen, indicating that the fluid−solid interaction energy is related to the affinity of the adsorbate for the coal surface. The εfs/k values also appear to be indicative of the polarity of the adsorbate, because CO2 is the most polar adsorbate of the three gases considered in this work. Further, the regressed εfs/k for nitrogen exhibited a linear correlation with εfs/k for methane (R2 of 0.98). The εfs/k for CO2, however, did not exhibit a significant linear correlation with εfs/k for either methane or nitrogen for the five coals. Further, the fluid−solid interaction 544

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Table 10. SLD Model Statistics for Pure-Gas Adsorption on Dry Coals and Activated Carbon case 1

case 2

WAADa

WAAD

coal

CH4

N2

CO2

overall

CH4

N2

CO2

overall

Pocahontas Upper Freeport Illinois no. 6 Wyodak Beulah Zap overall

0.58 0.46 0.48 0.21 0.64

0.69 0.58 0.67 0.64 1.06 0.6

0.68 0.48 0.70 1.03 0.45

0.65 0.51 0.62 0.63 0.72

0.38 0.30 0.43 0.17 0.64

0.89 0.73 0.91 0.66 1.07 0.7

0.76 0.51 0.99 1.30 0.85

0.67 0.51 0.78 0.71 0.86

F-400

0.96

0.33

0.75

0.70

0.85

0.94

0.8

Activated Carbon 0.7

case 3

case 4

WAAD

WAAD

coal

CH4

N2

CO2

overall

CH4

N2

CO2

overall

Pocahontas Upper Freeport Illinois no. 6 Wyodak Beulah Zap overall

0.38 0.30 0.43 0.17 0.64

0.87 0.74 0.80 0.62 1.07 0.7

0.76 0.78 1.09 1.53 0.85

0.67 0.60 0.78 0.77 0.86

0.38 0.30 0.43 0.17 0.64

0.87 0.74 0.80 0.62 1.07 0.8

0.76 0.51 1.12 2.33 0.85

0.67 0.51 0.78 1.04 0.86

F-400

0.70

1.38

2.64

0.70

1.38

2.04

1.4

Activated Carbon 1.6

case 1

case 2

% AADb

% AAD

coal

CH4

N2

CO2

overall

CH4

N2

CO2

overall

Pocahontas Upper Freeport Illinois no. 6 Wyodak Beulah Zap overall

2.36 2.04 2.68 1.04 4.66

5.37 5.55 9.83 7.78 16.46 5.0

3.27 2.67 4.25 4.30 2.40

3.66 3.42 5.58 4.37 7.84

1.28 1.11 2.26 0.86 4.75

6.85 6.89 12.59 7.90 16.59 5.5

3.63 2.86 5.39 5.52 4.66

3.92 3.62 6.75 4.76 8.67

F-400

1.18

0.53

2.03

0.88

1.72

2.13

1.6

Activated Carbon 1.3

case 3

case 4

% ADD

a

% AAD

coal

CH4

N2

CO2

overall

CH4

N2

CO2

overall

Pocahontas Upper Freeport Illinois no. 6 Wyodak Beulah Zap overall

1.28 1.11 2.26 0.86 4.75

6.94 6.85 9.66 8.27 16.59 5.6

3.63 4.72 6.24 6.30 4.66

3.95 4.23 6.06 5.14 8.67

1.28 1.11 2.26 0.86 4.75

6.94 6.85 9.66 8.27 16.59 5.7

3.63 2.86 6.43 8.55 4.66

3.95 3.61 6.12 5.89 8.67

F-400

0.88

2.98

Activated Carbon 6.00 3.3

0.88

2.98

4.69

2.9

b

Weighted average absolute deviation. Average absolute percentage deviation.

had large errors for two of the five coals (Beulah Zap and Wyodak coals). These two coals contain high-levels of oxygen and the as-received coals contained large amounts of equilibrium moisture.10 In fact, as shown in Table 9 for cases 1 and 2, the fluid−solid interaction energy of CO2 on Beulah Zap and Wyodak coals were much larger than for other coals. This may indicate a potential effect of the polar functional groups present in these coals on the (enhanced) molecular interactions of CO2 with the coal surface. Our analysis further showed that εfs/k for CO2 can be predicted within about 3% average absolute error as a linear function of equilibrium moisture of the as-received coal samples. Thus, for case 3, the

generalized models that were developed in this work (cases 3 and 4). Case 3. For the generalized model in case 3, the parameters obtained in case 2 for methane adsorption were retained. These included the common surface area A and slit length L. Then, the fluid−solid interaction energies εfs/k of nitrogen and CO2 on the five coals were generalized. The εfs/k for nitrogen was generalized as a linear function of εfs/k of methane and is given by eq 22. This expression was found to predict εfs/k for nitrogen with an average absolute error of 2.6% for the five coals used to develop the generalized expressions. The εfs/k for CO2 predicted from εfs/k of methane alone (similar to eq 22) 545

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Table 11. SLD Model Parameters for Pure-Gas Adsorption on Dry Coals (Data from Day et al.3) case 1

case 2

εfs/k (K)

εfs/k (K)

coal

A (m2/g)

CH4

CO2

L (nm)

A (m2/g)

CH4

CO2

L (nm)

Hunter Valley Illawarra Tashan

38.4 26.6 21.0

71.5 76.5 58.3

193.5 191.6 203.2

1.74 1.64 1.78

34.6 23.6 19.1

81.1 88.4 64.7

204.6 177.3 209.1

1.98 2.00 2.05

case 3

case 4

εfs/k (K)

εfs/k (K)

coal

A (m2/g)

CH4

CO2

L (nm)

A (m2/g)

CH4

CO2

L (nm)

Hunter Valley Illawarra Tashan

34.6 23.6 19.1

81.1 88.4 64.7

123.8 105.6 108.3

1.98 2.00 2.05

34.6 23.6 19.1

81.1 88.4 64.7

133.5 141.6 105.0

1.98 2.00 2.05

Table 12. SLD Model Statistics for Pure-Gas Adsorption on Dry Coalsa

a

case 1

case 2

case 3

case 4

% AAD

% AAD

% AAD

% AAD

Coal

CH4

CO2

overall

CH4

CO2

overall

CH4

CO2

overall

CH4

CO2

overall

Hunter Valley Illawarra Tashan

1.52 0.68 0.80

5.81 6.08 5.46

3.7 3.4 3.1

0.89 0.77 1.33

6.65 9.94 7.10

3.8 5.4 4.2

0.89 0.77 1.33

10.95 13.75 12.13

5.9 7.3 6.7

0.89 0.77 1.33

9.80 10.82 12.62

5.4 5.8 7.0

Data from ref 3.

εfs/k for CO2 was generalized based on the equilibrium moisture of coals, as shown in eq 23.

⎛ εfs ⎞nitrogen ⎛ ε ⎞methane ⎜ ⎟ = C1⎜ fs ⎟ + C2 ⎝k ⎠ ⎝k ⎠

(22)

⎛ εfs ⎞CO2 ⎜ ⎟ = C3(EM) + C4 ⎝k ⎠

(23)

for nitrogen and CO2 adsorption on coals were quite similar to direct regressions obtained in cases 1 and 2 and there appeared to be very little loss in accuracy from the use of generalized expressions. The generalization in case 3 was evaluated further by predicting the nitrogen and CO2 adsorption on activated carbon (F-400). The results (presented in Tables 9 and 10) provided predictions with a WAAD of 1.38 and 2.64 for nitrogen and CO 2 adsorption, respectively. Thus, the generalized expressions developed from coals alone were capable of predicting adsorption on activated carbon within three times the experimental uncertainties. This corresponded to an average absolute percentage deviation (% AAD) of 3% and 6% for nitrogen and CO2 adsorption, respectively. Day et al.3 presented data for the adsorption of pure methane and CO2 on three bituminous coals under dry and moist conditions. Day et al.3 dried their coal samples under vacuum at 60 °C (333 K) for 2 days. In our measurements reported in this work, the coal samples were dried under vacuum at 80 °C (353 K) for 36 h. Thus, the drying conditions were comparable, which aids in the analysis of these coal samples originating from different sources. The adsorption data on dry coals reported by Day et al.3 were used to further test/validate the generalized expressions developed in this work. Specifically, the methane adsorption data was used to predict the adsorption of CO2 on the three coals. The generalized predictions for CO2 adsorption using eq 23 are presented in Tables 11 and 12. As shown in Table 12, the generalized expression for case 3 yielded a % AAD of 11−14% for the adsorption of CO2 on three coals. This compares with direct regressions of CO2 data that yielded % AADs of about 5−6% (case 1) and about 7−10% (case 2), as shown in Table 12. Thus, the generalized expression for CO2 based on equilibrium moisture alone (eq 23) provides predictions with about two times the deviations observed through direct regressions for these coals.

where C1 to C4 were 0.566, −5.694, 3.866, and 97.10, respectively, and EM is percent equilibrium moisture of the sample. The generalization of fluid−solid interaction energy of CO2 on dry coals by utilizing equilibrium moisture of coals may, at first glance, seem counterintuitive to some readers. However, the equilibrium moisture level of a coal contains vital information about the presence of polar functional groups on the coal surface. The adsorption of CO2 and its molecular interactions with the coal surface are sensitive to the presence of polar (or oxygenated) functional groups on the coal surface. In fact, this has also been observed in the literature for CO2 adsorption on coals. Nishino et al.41 found the adsorption of CO2 on coals to be correlated with the carboxylic functional groups on the coal surface. Thus, the correlation of fluid−solid interaction energy of CO2 based on equilibrium moisture content of as-received coal samples appeared reasonable, because equilibrium moisture may be viewed as a coal characterization parameter. This generalization was tested further with an independent data set and will be discussed below. Tables 9 and 10 present the modeling/generalization results on coals for case 3. The results for methane adsorption on coals were identical to case 2, because the same parameters were used. The generalized expressions given in eqs 22 and 23 were used with the common surface area and slit length (based on methane only) to predict the adsorption of nitrogen and CO2 on five coals. As shown in Table 10, the generalized predictions 546

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

Case 4. A second approach was also investigated to predict the fluid−solid interaction energy of CO2 on coals and activated carbon. Because the polar functional groups on the coal surface appeared to have an effect on the molecular interactions of CO2 with the coal surface, a new expression was developed that utilizes the oxygen content of the coal. Specifically, the εfs/k for CO2 was generalized as a function of εfs/k for methane and an additional term to account for the oxygen content of coals. The generalized expression, given in eq 24, was found to predict εfs/k for CO2 with an average absolute error of 3.4% for the five coals (when compared with case 2). Further, the εfs/k for nitrogen was predicted with eq 22. In summary, case 4 comprised of predicting the εfs/k for nitrogen and CO2 from eqs 22 and 24, respectively.

6. CONCLUSIONS Adsorption isotherms for pure methane, nitrogen, and CO2 were measured on five coals of varying rank and an activated carbon at 328.2 K and pressures to 13.8 MPa using a volumetric method. The isotherms yielded expected experimental uncertainties between 1 and 7% depending on the gas adsorbed and the adsorbent considered (activated carbon or coals). The ratios of excess adsorption for the three gases showed larger variations among the coals covering wide ranges of rank and oxygen content. The SLD model was used to describe the adsorption data on coals and activated carbon. The overall WAAD for pure-gas adsorption on coals and activated carbon was 0.62 and 0.68 for the coals and activated carbon, respectively. Two generalized models based solely on the methane adsorption isotherm measurements and coal characterization information were presented that enable prediction of adsorption isotherms for nitrogen and CO2. The generalized models were found capable of predicting the nitrogen and CO2 adsorption within two times the experimental uncertainties for coals, on average. The generalized models were further tested on an external data set comprised of three dry coals from the literature. For each coal, the CO2 adsorption isotherms were predicted based on the methane isotherms. The generalized predictions yielded average absolute percentage deviation (% AAD) of about 10−12% for CO2 adsorption on the three coals, which was about two times the errors observed through direct model regressions for these coals.

⎛ εfs methane ⎞2 ⎛ εfs ⎞CO2 ⎜ ⎟ = C5⎜ ⎟ + C6(% Oxygen)3 ⎝k ⎠ ⎝k ⎠ + C7

(24)

where C5 to C7 were 0.0162, 0.0176, and 13.151, respectively. Tables 9 and 10 present the generalization results on coals for case 4. The results for methane and nitrogen adsorption on coals were identical to case 3, because the same parameters were used. In other words, the only difference between case 3 and 4 is in the prediction of CO2 adsorption. For case 4, the generalized expression given by eq 24 predicted the CO2 adsorption on the five coals with similar accuracy as in case 3, as evident from Table 10. The generalized expression (eq 24) was also used to predict the adsorption of CO2 on activated carbon. As shown in Table 10, the WAAD and % AAD for CO2 adsorption on activated carbon was 2.0 and 4.7%, respectively, for case 4. This compares with a WAAD and % AAD of 2.6 and 6%, respectively, for case 3. Thus, the predictions from case 4 were slightly better than those in case 3 for activated carbon. Table 11 and 12 also present the generalized predictions for case 4 obtained for the CO2 adsorption on three coals from an independent data set presented by Day et al.3 The model deviations were marginally better than in case 3 and yielded a % AAD of 10−13% for CO2 adsorption. This compares with a % AAD of about 5−6% (case 1) and 7−10% (case 2), which refers to direct regressions of CO2 adsorption data, as shown in Table 12. Note that the generalized expressions in cases 3 and 4 were developed only on the basis of the adsorption data for the five coals measured in this work, and the CO2 adsorption on activated carbon and the three coals from the literature was then predicted directly. Both cases 3 and 4 provided prediction results with similar accuracy, as evident from Tables 10 and 12. Thus, there appears to be two alternative and apparently equivalent methods to predict the adsorption of other gases from the knowledge of methane adsorption isotherms and the coal characterization information of the adsorbent. Notwithstanding the promising results obtained to date, additional testing of these expressions with larger data sets is needed to establish the viability of this modeling approach. Further, the expressions in this work have been developed for pure-gas adsorption on dry adsorbents alone. Because most coalbeds contain water/moisture, an extension of this approach to wet adsorbents would be beneficial in predicting supercritical gas adsorption for systems of interest in coalbed methane work.



ASSOCIATED CONTENT

S Supporting Information *

Helium void volume as a function of pressure at 328.2 K. This material is available free of charge via the Internet at http:// pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*Telephone: (405) 744-5280. Fax: (405) 744-6338. E-mail: [email protected].



ACKNOWLEDGMENTS The financial support of U.S. Department of Energy and Advanced Resources International, Inc. is gratefully acknowledged.



REFERENCES

(1) White, C. M.; Smith, D. H.; Jones, K. L.; Goodman, A. L.; Jikich, S. A.; LaCount, R. B.; DuBose, S. B.; Ozdemir, E.; Morsi, B. I.; Schroeder, K. T. Sequestration of Carbon Dioxide in Coal with Enhanced Coalbed Methane RecoveryA Review. Energy Fuels 2005, 19 (3), 659−724. (2) Day, S.; Duffy, G.; Sakurovs, R.; Weir, S. Effect of Coal Properties on CO2 Sorption Capacity under Supercritical Conditions. Int. J. Greenhouse Gas Control 2008, 2 (3), 342−352. (3) Day, S.; Sakurovs, R.; Weir, S. Supercritical Gas Sorption on Moist Coals. Int. J. Coal Geol. 2008, 74 (3−4), 203−214. (4) Busch, A.; Gensterblum, Y.; Krooss, B. M.; Siemons, N. Investigation of High-Pressure Selective Adsorption/Desorption Behavior of CO2 and CH4 on Coals: An Experimental Study. Int. J. Coal Geol. 2006, 66 (1−2), 53−68. (5) Busch, A.; Gensterblum, Y.; Krooss, B. M. Methane and CO2 Sorption and Desorption Measurements on Dry Argonne Premium Coals: Pure Components and Mixtures. Int. J. Coal Geol. 2003, 55 (2− 4), 205−224. 547

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548

Energy & Fuels

Article

(26) Rangarajan, B.; Lira, C. T.; Subramanian, R. Simplified Local Density Model for Adsorption over Large Pressure Ranges. AIChE J. 1995, 41 (4), 838−845. (27) Chen, J. H.; Wong, D. S. H.; Tan, C. S.; Subramanian, R.; Lira, C. T.; Orth, M. Adsorption and Desorption of Carbon Dioxide onto and from Activated Carbon at High Pressures. Ind. Eng. Chem. Res. 1997, 36 (7), 2808−2815. (28) Puziy, A. M.; Herbst, A.; Poddubnaya, O. I.; Germanus, J.; Harting, P. Modeling of High-Pressure Adsorption Using the Bender Equation of State. Langmuir 2003, 19 (2), 314−320. (29) Soule, A. D.; Smith, C. A.; Yang, X.; Lira, C. T. Adsorption Modeling with the ESD Equation of State. Langmuir 2001, 17 (10), 2950−2957. (30) Yang, X.; Lira, C. T. Theoretical Study of Adsorption on Activated Carbon from a Supercritical Fluid by the SLD-ESD Approach. J. Supercrit. Fluids 2006, 37 (2), 191−200. (31) Fitzgerald, J. E.; Sudibandriyo, M.; Pan, Z.; Robinson, R. L.; Gasem, K. A. M. Modeling the Adsorption of Pure Gases on Coals with the SLD Model. Carbon 2003, 41 (12), 2203−2216. (32) Fitzgerald, J. E.; Robinson, R. L.; Gasem, K. A. M. Modeling High-Pressure Adsorption of Gas Mixtures on Activated Carbon and Coal Using a Simplified Local-Density Model. Langmuir 2006, 22 (23), 9610−9618. (33) Fitzgerald, J. E. Adsorption of Pure and Multi-Component Gases of Importance to Enhanced Coalbed Methane Recovery: Measurements and Simplified Local-Density Modeling. Ph.D. Dissertation, Oklahoma State University, Stillwater, OK 2005. (34) Peng, D.-Y.; Robinson, D. B.; New Two-Constant, A Equation of State. Ind. Eng. Chem. Fund. 1976, 15 (1), 59−64. (35) Gasem, K. A. M.; Gao, W.; Pan, Z.; Robinson, R. L. Jr. A Modified Temperature Dependence for the Peng−Robinson Equation of State. Fluid Phase Equilib. 2001, 181 (1−2), 113−125. (36) Lee, L. L. Molecular Thermodynamics of Nonideal Fluids; Butterworth: Stoneham, MA, 1988. (37) Steele, W. A. The Interaction of Gases with Solid Surfaces; Pergamon Press: Oxford, 1974. (38) Subramanian, R.; Pyada, H.; Lira, C. T. Engineering Model for Adsorption of Gases onto Flat Surfaces and Clustering in Supercritical Fluids. Ind. Eng. Chem. Res. 1995, 34 (11), 3830. (39) Reid, R. C.; Prausnitz, J. M.; Poling, B. E. The Properties of Gases and Liquids; McGraw-Hill: New York, 1987. (40) Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; Vetterling, W. T. Numerical Recipes in Fortran: The Art of Scientific Computing; Cambridge University Press: Cambridge, 1992. (41) Nishino, J. Adsorption of Water Vapor and Carbon Dioxide at Carboxylic Functional Groups on the Surface of Coal. Fuel 2001, 80 (5), 757−764.

(6) Hall, F.; Zhou, C.; Gasem, K. A. M.; Robinson, R. L., Jr. Adsorption of Pure Methane, Nitrogen, and Carbon Dioxide and Their Binary Mixtures on Wet Fruitland Coal, SPE Paper 29194. SPE Eastern Regional Conference & Exhibition, Charleston, SC, 1994. (7) Fitzgerald, J. E.; Pan, Z.; Sudibandriyo, M.; Robinson, R. L. Jr.; Gasem, K. A. M.; Reeves, S. Adsorption of Methane, Nitrogen, Carbon Dioxide and Their Mixtures on Wet Tiffany Coal. Fuel 2005, 84 (18), 2351−2363. (8) Mohammad, S. A.; Chen, J. S.; Fitzgerald, J. E.; Robinson, R. L.; Gasem, K. A. M. Adsorption of Pure Carbon Dioxide on Wet Argonne Coals at 328.2 K and Pressures up to 13.8 MPa. Energy Fuels 2009, 23 (2), 1107−1117. (9) Clarkson, C. R.; Bustin, R. M. Binary Gas Adsorption/Desorption Isotherms: Effect of Moisture and Coal Composition upon Carbon Dioxide Selectivity over Methane. Int. J. Coal Geol. 2000, 42 (4), 241− 271. (10) Vorres, K. S. The Argonne Premium Coal Sample Program. Energy Fuels 1990, 4 (5), 420−426. (11) Mohammad, S. A.; Chen, J. S.; Robinson, R. L.; Gasem, K. A. M. Generalized Simplified Local-Density/Peng-Robinson Model for Adsorption of Pure and Mixed Gases on Coals. Energy Fuels 2009, 23 (12), 6259−6271. (12) Goodman, A. L.; Busch, A.; Duffy, G. J.; Fitzgerald, J. E.; Gasem, K. A. M.; Gensterblum, Y.; Krooss, B. M.; Levy, J.; Ozdemir, E.; Pan, Z.; Robinson, R. L.; Schroeder, K.; Sudibandriyo, M.; White, C. M. An Interlaboratory Comparison of CO2 Isotherms Measured on Argonne Premium Coal Samples. Energy Fuels 2004, 18 (4), 1175−1182. (13) Sakurovs, R.; Day, S.; Weir, S.; Duffy, G. Temperature Dependence of Sorption of Gases by Coals and Charcoals. Int. J. Coal Geol. 2008, 73 (3−4), 250−258. (14) Sudibandriyo, M.; Pan, Z.; Fitzgerald, J. E.; Robinson, R. L. Jr.; Gasem, K. A. M. Adsorption of Methane, Nitrogen, Carbon Dioxide, and Their Binary Mixtures on Dry Activated Carbon at 318.2 K and Pressures up to 13.6 MPa. Langmuir 2003, 19 (13), 5323−5331. (15) Angus, S.; Armstrong, B.; de Reuck, K. M. International Thermodynamic Tables of the Fluid State-5: Methane; Pergamon Press: New York, 1978. (16) Angus, S.; de Reuck, K. M.; Armstrong, B. International Thermodynamic Tables of the Fluid State-6: Nitrogen.; Pergamon Press: New York, 1979. (17) Span, R.; Wagner, W.; New, A Equation of State for Carbon Dioxide Covering the Fluid Region from the Triple Point Temperature to 1100 K at Pressures up to 800 MPa. J. Phys. Chem. Ref. Data 1996, 25, 1509−1590. (18) McCarty, R. D. Thermophysical Properties of Helium-4 from 2 to 1500 K with Pressures to 1000 atm, NBS Technical Note 631; U.S. Department of Commerce: 1972. (19) Hall, F. E. Adsorption of Pure and Multicomponent Gases on Wet Fruitland Coal, M.S. Thesis, Oklahoma State University, Stillwater, 1993. (20) Mohammad, S. A.; Fitzgerald, J. E.; Robinson, R. L. Jr.; Gasem, K. A. M. Experimental Uncertainties in Volumetric Methods for Measuring Equilibrium Adsorption. Energy Fuels 2009, 23 (5), 2810− 2820. (21) Ozdemir, E.; Morsi, B. I.; Schroeder, K. Importance of Volume Effects to Adsorption Isotherms of Carbon Dioxide on Coals. Langmuir 2003, 19 (23), 9764−9773. (22) Dutta, P.; Harpalani, S.; Prusty, B. Modeling of CO2 Sorption on Coal. Fuel 2008, 87 (10−11), 2023−2036. (23) Romanov, V.; Soong, Y.; Schroeder, K. Volumetric Effects in Coal Sorption Capacity Measurements. Chem. Eng. Technol. 2006, 29 (3), 368−374. (24) Pan, Z.; Connell, L. D. A Theoretical Model for Gas Adsorption-Induced Coal Swelling. Int. J. Coal Geol. 2007, 69 (4), 243−252. (25) Day, S.; Fry, R.; Sakurovs, R. Swelling of Australian Coals in Supercritical CO2. Int. J. Coal Geol. 2008, 74 (1), 41−52. 548

dx.doi.org/10.1021/ef201393p | Energy Fuels 2012, 26, 536−548