Permeability prediction in deep coal seam: a case study on the no. 3 coal seam of the Southern Qinshui Basin in China.
Coal bed methane (CBM) is a natural product in the coalification process [1, 2]. CBM is a serious threat to safety in underground coal mining and can cause disasters, such as coal and gas outbursts and gas explosions [1, 3]. However, CBM is also an unconventional natural gas resource that has been exploited worldwide in countries such as USA, Australia, and China [4, 5].
Coal permeability is an important parameter in mine methane control and CBM exploitation because it determines the practicability of methane extraction. The permeability determined in the laboratory is not the real value of coal seam tested in the field. The coal seam permeability is completed based on the well logs in the CBM exploitation and is calculated based on the gas flow theory of the hole in or through the coal seam during coal mining [2,6]. In the region where the exploration or mining does not conduct, especially in deep coal seam, it requires a large number of engineering works to test the permeability, which is very expensive. The permeability of deep coal seam predicted by a permeability model is fast, reliable, and economical for preassessment of CBM exploitation and methane extraction.
The permeability of coal depends on the fracture characteristics, including the size, spacing, connectivity, width, mineral fill, and distribution , which are affected by the strata stress, gas pressure, and temperature of the coal seam.
Several models have been proposed to explain the variability of coal permeability [8-19]. The coal permeability models can be divided into two important classes : those under uniaxial strain conditions [8-13] and those under triaxial stress conditions [14-19].
However, uniaxial strain conditions are a simplified homogenisation of the stress-strain states of coal during mining and exploitation and may be valid at the scale of a relatively large basin; the mechanical conditions at the local scale are expected to be much more complex in coal seams . And the tri-axial strata stress conditions change with depth. Therefore, it is appropriate for predicting the permeability of deep coal seam using a permeability model under conditions of tri-axial stress-strain.
The effects of effective stress and matrix sorption deformation are considered by most coal permeability models [8, 9,11-19] except for Seide and Huitt's model who assumed that cleat deformation was caused entirely by desorption shrinkage . The coal matrix is assumed to be incompressible by assuming that the bulk modulus of the coal matrix is much larger than the coal bulk modulus and then Biot's coefficient a is assumed to be 1 [8-13,16,17]. However, the compression of the coal matrix by the pore pressure could not be ignored [20, 21]. Therefore, the Biot's coefficient for coal is less than 1 [17, 22, 23].
Most models consider the matrix deformation to be equal to the fracture deformation. However, only part of the matrix deformation contributes to the fracture deformation . Liu and Rutqvist  and Connell et al.  established permeability models in which the sorption deformation applied to the fracture.
The strata stress, gas pressure, and temperature of the coal seam change with the depth [25-29]. As a result, the permeability of the coal seam changes with depth. In this paper, we discussed the effect of stress, gas pressure, and temperature on the porosity of the coal (the coal matrix and the fracture) under tri-axial stress-strain conditions. Then we established a permeability model for deep coal seam. At last, we predicted the permeability of the No. 3 coal seam of South Qinshui Basin (SQB) in China.
2. Establishment of the Permeability Model
The coal has a natural dual porosity structure that consists of the coal matrix and the fracture in which there are numerous inorganic minerals, mainly kaolinite, pyrite, and illite, as shown in Figure 1. More than 95% of the gas occurs as adsorbed gas in the sorption space of the abundant micropores . The gas migrates by diffusion in the micropore system and follows Fick's Law. The closely spaced natural fractures surrounding the coal matrix, which form the cleat system, determine the mechanical properties of the coal and the flow paths for the methane; this flow follows Darcy's Law. Therefore, the coal fracture permeability is closely related to the characteristics of the fractures, which are controlled by the coal rank, geologic structure, mining, strata stress, gas pressure, formation temperature, and other factors. The coal fracture porosity is controlled by the strata stress, gas pressure, and temperature of the coal seam which change with the depth. Therefore, we analyzed the effect of the strata stress, gas pressure, and temperature of the coal seam on the fracture in the later section.
Before analysis, we make the following assumptions.
(1) Coal is considered to be a dual continuous isotropic elastic medium even though the coal consists of coal matrix and fracture.
(2) The strain is elastic and infinitesimal, so the second and higher order terms can be ignored. Therefore, the strains induced by different factors can be added.
(3) The effect of methane on the coal is divided into normal effective stress effect of the gas like nonadsorptive gas and the effect of sorption deformation, which are added directly.
2.1. Effective Stress. According to the effective stress principle , the bulk volumetric strain increment and the pore volume strain increment can be expressed as [19, 31]
d[[epsilon].sub.e] = d[V.sub.e]/V = -1/K([bar.d[sigma]] - [alpha]dp), (1)
d[[epsilon].sub.ep] = d[V.sub.ep]/V = -1/[K.sub.p]([bar.d[sigma]] - [beta]dp), (2)
where [[epsilon].sub.e] and [V.sub.e] are the coal bulk strain and volume caused by the effective stress, respectively; [[epsilon].sub.ep] and [V.sub.ep] are the pore strain and volume caused by the effective stress, respectively; <r = (1/3)([[sigma].sub.11] + [[sigma].sub.22] + [[sigma].sub.33]) is the mean stress, MPa; K = E/3(1-2v) is the coal bulk modulus, MPa; [K.sub.p] is the coal pore modulus, MPa; [alpha] = 1 -K/[K.sub.m] is Biot's coefficient; [beta] = 1 - [K.sub.p]/[K.sub.m] is the effective coefficient for the pore system; Km is the coal matrix modulus, MPa; E is the elastic modulus of the coal, MPa; and v is Poisson's ratio.
Without the gas sorption effect, the volumetric change of the porous medium satisfies the Betti-Maxwell reciprocal theorem [([partial derivative]V/[partial derivative]p).sub.[bar.[sigma]]] = [([partial derivative][V.sub.p]/[partial derivative] [bar.[sigma]]).sub.p], and we obtain
[K.sub.p] = [phi]/[alpha] K, (3)
where [phi] is the porosity.
2.2. The Coal Matrix Deformation. The coal matrix swells due to adsorb gas and temperature increase underground.
(1) Thermal Deformation. The deformation of the coal matrix due to temperature change is expressed as
d[[epsilon].sub.mt] = [eta]dT, (4)
where [[epsilon].sub.mt] is the thermal deformation of the coal matrix; [eta] is the coefficient of the thermal deformation, [K.sup.-1]; T is the temperature of the strata, K.
(2) Sorption Deformation. The gas sorption capacity of coal increases with pressure and closely follows Langmuir type isotherm. Langmuir isotherm describes the amount of adsorbed gas on the coal as a function of pressure and has the following form :
[V.sub.ad] = abp/1 + bp, (5)
where [V.sub.ad] is the adsorbed gas volume, [m.sup.3]/t; a is the Langmuir volume, [m.sup.3]/t; b is the adsorption coefficient, [MPa.sup.-1].
Adsorption coefficient indicates coal's ability to adsorb gases and is exponentially increased with the energy of adsorption and decreased with the temperature of adsorption. Adsorption coefficient in a Langmuir-type isotherm can be estimated using the Arrhenius rate equation at equilibrium condition, in which affinity is a function of temperature and is given by the following equation :
b = [b.sub.0]/[square root of (T)] exp ([DELTA]H/RT), [b.sub.0] = S[alpha]/[k.sub.d[infinity]] [square root of (2[pi]MR)], (6)
where [DELTA]H is the heat of sorption, J/mol; R is the gas universal constant, J/mol x K.
This is to say that the adsorption capacity is independent of temperature, and as a result the heat of adsorption is a constant, independent of loading. Increase in the temperature will decrease the adsorbed amount at a given pressure. This is due to the greater energy acquired by the adsorbed molecule to evaporate.
It has been tested and verified that the coal swells it when adsorbing gas and the adsorption deformation is proportional to the adsorbed gas volume [22, 24, 34]:
[[epsilon].sub.s] = [chi][V.sub.ad] = [chi] abp/1 + bp, (7)
where [[epsilon].sub.s] is the sorption deformation of the coal matrix; [chi] is the coefficient of sorption deformation, t/[m.sup.3].
However, the coal matrix is also compressed by the sorptive gas in addition to swelling due to adsorption [20, 21, 35]. When testing the sorption deformation, we need the following equation to calibrate the experimental data:
[[epsilon].sub.s] = [[epsilon].sub.exp] + P/[K.sub.s], (8)
where [[epsilon].sub.exp] is the deformation measured in the experiment.
Therefore, the coal matrix deformation is the sum of adsorption deformation and thermal deformation:
[[epsilon].sub.m] = [[epsilon].sub.s] + [[epsilon].sub.T]. (9)
(3) Fracture Deformation Caused by the Coal Matrix. Deformation of the coal matrix can affect the deformation of both the bulk coal and the fractures in the coal [10, 11, 14, 36]. The coal matrix deformation is assumed to contribute entirely to the fracture deformation [10,11,15]. However, the contribution of coal matrix deformation to the fracture has been significantly overestimated [16, 17, 24]. For example, Roberston and Christiansen  demonstrated that the most commonly used models [11,12] significantly overestimate the effects of matrix swelling on the permeability changes which was observed in laboratory experiments.
The coal is divided into the matrixes by the fracture, but there are bridges between the matrixes, as shown in Figure 1. When the coal matrix swells, the bridge limits the coal matrix deformation to the fracture . Karacan  and Dawson et al.  found that numerous inorganic minerals, mainly kaolinite, pyrite, and illite, are present in coal fractures. These minerals prevent the coal matrix from completely closing the fracture. Therefore, only part of the matrix deformation contributes to the fracture deformation.
The effective coal matrix deformation factor, [f.sub.m], is introduced to measure the degree of influence of the coal matrix deformation on fracture deformation. The factor [f.sub.m] is a parameter of the coal structure and depends on the distribution of fractures, the filling characteristics of the fracture, and other factors. For a particular coal, [f.sub.m] is a constant between 0 and 1. If there is no fracture in the coal, the parameter [f.sub.m] is equal to 0. The parameter [f.sub.m] would be equal to 1 when two surfaces of the fracture are smooth and parallel.
Thus, the fracture deformation due to the coal matrix deformation is expressed as
dV[m.sub.f] = = [f.sub.m]d[V.sub.m] = [f.sub.m][V.sub.m] (d[[epsilon].sub.s] + d[[epsilon].sub.T]), (10)
where V[m.sub.f] is the fracture volume deformation due to deformation of coal matrix; [V.sub.m] is the volume of coal matrix.
2.3. The Permeability Model under Triaxial Stress Conditions. As a porous medium, the coal bulkvolume V is composed of the matrix volume [V.sub.m] and the pore volume [V.sub.p]:
V = [V.sub.p] + [V.sub.m]. (11)
Based on the definition of porosity, [phi] = [V.sub.p]/V, we obtain
d[phi] = d([V.sub.p]/V) = [phi](d[V.sub.p]/[V.sub.p] - dV/V). (12)
The bulk volume deformation of coal is equal to the sum of the deformation due to the effective stress and the coal matrix deformation due to adsorption and temperature change:
dV = d[V.sub.e] + d[V.sub.mv] = -V/K ([bar.d[sigma]] - [alpha]dp) + (1 - [f.sub.m]) [V.sub.m] d[[epsilon].sub.m], (13)
where [V.sub.mv] is the bulk volume deformation due to the matrix deformation.
Dividing both sides of (13) by the coal bulk volume, we obtain
dV/V = -1/K ([bar.d[sigma]] - [alpha]dp) + (1 - [f.sub.m]) (1 - [phi]) d[[epsilon].sub.m]. (14)
Similarly, from (2) and 10), we obtain
d[V.sub.p]/[V.sub.p] = -1/[K.sub.p]([bar.d[sigma]] - [beta]dp) - 1/- [phi]/[phi] [f.sub.m] d[[epsilon].sub.m]. (15)
By substituting (14) and 15) into 12), we obtain
d[phi]/[phi] = -1/[K.sub.p] ([bar.[sigma]] - [beta]p) + 1/K ([bar.[sigma]] - [alpha]p) - (1 - [phi]/[phi] [f.sub.m] + (1 - [f.sub.m]) (1 - [phi])) d[[epsilon].sub.m]. (16)
Then, substituting [K.sub.p] = ([phi]/[alpha])K and [beta] = 1 - [K.sub.p]/[K.sub.m] into (16) and considering that [phi] [much less than] 1 (0 < 10%), we can rearrange and simplify the equation to obtain
d[phi] = -[alpha]/K ([bar.d[sigma]] - dp) -[f.sub.m]d[[epsilon].sub.m]. (17)
Integrating (17) gives
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)
The widely used cubic relationship between permeability and porosity [8-17] is given
k/[k.sub.0] = [([phi]/[[phi].sub.0]).sup.3], (19)
where k is the coal permeability.
Substituting (18) into (19), the coal permeability model that considers the effect of the effective stress and coal matrix deformation (ESMD model) is given
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)
It is clear that the model contains an effective stress term and a coal matrix deformation term. The factor [f.sub.m] measures the degree of influence of the coal matrix deformation on the fracture deformation.
3. Model Validation and Evaluation
3.1. Experimental Data. Numerous laboratory experiments have been conducted on coal permeability [24, 39, 40]. Pini et al.  conducted experiments that tested the mechanical parameters, porosity, adsorption swelling parameters, and coal permeability of a coal sample (Sulcis coal sample) from the Monte Sinni coal mine in the Sulcis Coal Province (Sardinia, Italy). We use the experimental data to validate and evaluate the ESMD model because of the comprehensive set of parameters available for the coal sample and the detailed experimental data. The coal permeability experiments were conducted under hydrostatic conditions at a constant confining pressure (10 MPa) and various gas pressures between 0 MPa and 8 MPa at 45[degrees]C using [N.sub.2] and C[O.sub.2]. The adsorption swelling parameters of the Sulcis coal sample for [N.sub.2] and C[O.sub.2] were corrected using (8). The parameters of the Sulcis coal sample are shown in Table 1.
3.2. Validation. The experimental data were matched by the ESMD model, the Palmer-Mansoori (P-M) model , the Shi-Durucan (S-D) model , and the RobertsonChristiansen (R-C) model  using the parameters in Table 1. The results are shown in Figure 3. The ESMD model can predict the experimental data well.
Only part of the matrix deformation contributes to the fracture deformation. The factor [f.sub.m], which ranges from 0 to 1, is introduced to measure the degree of influence of the coal matrix deformation on the fracture deformation in the ESMD model. The factor fm is a parameter of the coal structure and does not vary with the type of gas. The factor [f.sub.m] of the Sulcis coal sample is 0.1723 for both [N.sub.2] and C[O.sub.2]. Biot's coefficient a of coal is less than 1, and a = 0.925 for the Sulcis coal sample, which has a bulk modulus of 778 MPa and a matrix modulus of 10,340 MPa.
The three models poorly match the experimental data for two reasons: in all three models, Biot's coefficient is assumed to be 1 in the P-M model and the S-D model by assuming that the coal matrix is incompressible. Deformation of the coal matrix contributes to the fracture deformation entirely in the three models, which is an overestimation. As shown in Figure 2, the R-C model matches the experimental data well for [N.sub.2] but poorly for C[O.sub.2]. The decline rate of fracture compressibility with increasing effective stress [theta] for the R-C model is 2.65 x [10.sup.-14] [MPa.sup.-1], which implies that the fracture compressibility does not vary with the effective stress. However, the decline rate [theta] varies between 2.45 x [10.sup.-2] [MPa.sup.-1] and 2.61 x [10.sup.-1] [MPa.sup.-1] [14,24, 41].
4. The Parameters of the No. 3 Coal Seam in the SQB
The SQB has become the China's first commercial CBM reservoir. CBM reservoir is an unconventional gas reservoir and located at shallow depths, compared to the conventional gas reservoir. The SQB refers to a region, including Changzhi, Gaoping, Jincheng, Yangcheng, Qinshui, and Anze in the southeast of Shanxi Province. It is the most important production base for high quality anthracite in China. Coal seams, generated in Carboniferous and Permian periods, contain abundant methane. Permeability in the coalbed reservoir is relatively high compared to other CBM reservoirs in China. The exploration and production tests in this field have been conducted since 1990s. The results show that the Qinshui Basin is a very promising coalbed methane reservoir with the most exploration wells, the best development prospect, and a higher commercialized production in China's CBM reservoirs. The No. 3 coal seam is the major coal seam mined and CBM reservoir in the SQB. The thickness of No. 3 coal seam in the SQB changes from 0.5 m to 7.0 m with an average value of 5.5 m. The maximum depth is about 2000 m. The coal is anthracite mainly. In this paper, we predicted the permeability of the No. 3 coal seam.
4.1. Gas Pressure. There are two different directions of gas migration in the coal seam with an outcrop. For example, the generated methane migrates upward and the surface air migrates downward in the coal seam. Thus, there are two vertical zones of the CBM occurrence in the coal seam, which are the gas weathered zone and the methane zone . Generally, they are zoned by the methane concentration of the coal seam with the value of 80%. The methane concentration change of the No. 3 coal seam in the SQB with the depth is present in Figure 4 . The depth of the gas weathered zone of the No. 3 coal seam is about 160 m.
The gas pressure of the No. 3 coal seam in the SQB is present in Figure 5. It has been found that the relationship between gas pressure and depth is linear by analyzing the numerous measured gas pressure values [1, 2, 26-28]. There are many factors that influence the gas pressure which could lead to a deviation from actual values, and the measured data do not possess basic conditions for regression methods [27, 28]. A pressure prediction method, the safety line method, was used to analyse the variation of the gas pressure with depth . Two true symbol points were selected to make the safetyline. All the other points except the abnormal points are below the line, as shown in Figure 5. The gas pressure gradient is 0.0136 MPa/m. The relationship between the gas pressure and the depth is expressed as
p = 0.0136 x H - 2.1003, (21)
where H is the depth of coal seam, m.
The gas pressure calculated by (21) at the depth of 161 m is 0.1 MPa. The result is consistent with the depth of the gas weathered zone because the gas pressure at the depth of the gas weathered zone is from 0.1 to 0.15 MPa statistically in China [1,27].
4.2. The Strata Stress. It has been practically verified that the strata stress increases linearly with depth [25,43]. Meng et al.  investigated the stratum stresses in the SQB and they are expressed as
[[sigma].sub.V] = 0.027 x H, [[sigma].sub.H] = 0.0304 x H -3.1975, [[sigma].sub.h] = 0.0235 x H -3.5127, (22)
where [[sigma].sub.V] is the vertical stress, MPa; [[sigma].sub.H] is the maximum horizontal stress, MPa; [[sigma].sub.h] is the minimum horizontal stress, MPa.
4.3. The Strata Temperature. The temperature of the strata also increases linearly with the depth. Sun et al.  investigated the strata temperature in SQB and found that the relationship between the strata temperature and the depth is
T = G x (H - 20) + 282.15. (23)
4.4. The Parameters of No. 3 Coal Seam. The parameters of the No. 3 coal seam in the SQB have been measured by many researchers [26, 42, 44-47]. The parameters on average are presented in Table 2. Bangham and Franklin  and Kelemen and Kwiatek  measured the thermal deformation coefficient of different coal samples at the temperature less than 100[degrees]C. The thermal deformation coefficient was 38.13 x [10.sup.-6] [K.sup.-1], which was used in the prediction.
The prediction of the permeability of the No. 3 coal seam in the SQB was conducted by using (6) and (20) to (23) based on the parameters in Table 2. The results were shown in Figure 6.
(1) The statistics of the permeability of the No. 3 coal seam tested in the field was investigated by Meng et al.  and was presented in Figure 7. The permeability model could predict the deep coal permeability of the No. 3 coal seam in the SQB efficiently.
(2) The average strata stress, the gas pressure, and the strata temperature of the No. 3 coal seam in the SQB increase linearly with the depth. But the gradient of the average strata stress is larger than that of the gas pressure (Figure 8). So the effective stress also increases with the depth, which compresses and closes the fracture narrower. The fracture is also closed more narrowly by the adsorption swelling and thermal expansion increasing with the depth increase. As a result, the permeability of the No. 3 coal seam in the SQB decreases with the increase of the depth.
(3) There is a cubic relationship between the coal permeability and the coal porosity. The different effects of the effective stress, the adsorption deformation, and the thermal deformation of the coal seam on the permeability could be represented by different effects of them on the porosity which are calculated easily and simply. The different effects of the effective stress, the adsorption deformation, and the thermal deformation were calculated using (18), shown in Figure 9. The negative means that the effect of the facts closes the fracture.
The effect of the adsorption swelling on the porosity increases rapidly at first slowly with the increase of depth. The effect of the thermal expansion and the effective stress compression increases linearly with the depth increase. The effect of the effective stress compression is larger than the others.
As a result, the most effective way to improve the permeability in deep and low permeability coal bed for exploiting CBM or extracting methane is to reduce the effective stress. For example, the method of protective layer mining, which is widely used in China, is very effective and it reduces the stress of the overlying and underlying strata .
(4) The temperature affects the sorption capacity of the coal. The sorption capacity decreases with the increase of the temperature. The adsorption volume increases with the increase of the depth to a maximum value at the depth of 1600 m, and then it decreases. However, if the effect of the temperature on the sorption capacity is neglected, the adsorption volume increases always with the increase of the depth. And the adsorption volume is overestimated (Figure 10). As a result, the effect of the adsorption swelling is also overestimated (Figure 9).
In the ESMD model proposed in this paper, only part of the matrix deformation contributes to the fracture deformation. If the matrix deformation contributes to the fracture deformation, the effect of the matrix deformation due to adsorption swelling and thermal expansion is overestimated significantly (Figure 9).
The Biot's coefficient is less than 1 in the ESMD model, which is true for the coal. But in mangy coal permeability models, it is equal to 1 based on the assumption that the bulk modulus of the coal matrix is much larger than that of the coal, which overestimates the effect of the effective stress (Figure 9).
The coal permeability is an important parameter in mine methane control and CBM exploitation, which determines the practicability of methane extraction. Permeability prediction in deep coal seam plays a significant role in evaluating the practicability of CBM exploitation.
The coal permeability depends on the coal fractures controlled by strata stress, gas pressure, and strata temperature which change with depth. The effect of the strata stress, gas pressure, and strata temperature on the coal (the coal matrix and fracture) under tri-axial stress and strain conditions was studied. Then we got the change of coal porosity with strata stress, gas pressure, and strata temperature and established a coal permeability model (ESMD model) under tri-axial stress and strain conditions. The coal permeability is controlled by the effective stress, sorption deformation, and thermal expansion. The sorption capacity decreases with the increase of the temperature and thus the sorption deformation is affected by temperature. And only part of the matrix deformation contributes to the fracture deformation. The Biot's coefficient is less than 1 in the ESMD model, which is true for the coal.
The permeability of the No. 3 coal seam of the SQB in China was predicted, which is consistent with that tested in the field.
The effect of the sorption deformation on porosity (permeability) firstly increases rapidly and then slowly with the increase of depth. However, the effect of thermal expansion and effective stress compression on porosity (permeability) increases linearly with the increase of depth. When exploiting CBM or extracting methane, the most effective way to improve the permeability in CBM exploitation or mine methane control is to reduce the effective stress. For example, the method of protective layer mining, which is widely used in China, is very effective, which reduces the stress of the overlying and underlying strata.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are grateful for the support of the National Foundation of China (nos. 51004106, 41202118, 51204173, 51374204, and 51304188) and the National Basic Research Program of China (973 Program, no. 2011CB201204).
 Q. X. Yu, Mine Gas Prevention and Control, China University of Mining and Technology Press, XuZhou, China, 1992.
 S. N. Zhou and B. Q. Lin, The Theory of Gas Flow and Storage in Coal Seams, China Coal Industry Publishing House, Beijing, China, 1997.
 C. O. Karacan, F. A. Ruiz, M. Cote, and S. Phipps, "Coal mine methane: a review of capture and utilization practices with benefits to mining safety and to greenhouse gas reduction," International Journal of Coal Geology, vol. 86, no. 2-3, pp. 121-156, 2011.
 J. Liu, Z. Chen, D. Elsworth, H. Qu, and D. Chen, "Interactions of multiple processes during CBM extraction: a critical review," International Journal of Coal Geology, vol. 87, no. 3-4, pp. 175-189, 2011.
 T A. Moore, "Coalbed methane: a review," International Journal of Coal Geology, vol. 101, pp. 36-81, 2012.
 M.-J. Liu and X.-Q. He, "Optimistic calculation method of gas seepage coefficients," Journal of the China Coal Society, vol. 29, no. 1, pp. 74-77, 2004.
 S. E. Laubach, R. A. Marrett, I. E. Olson, and A. R. Scott, "Characteristics and origins of coal cleat: a review," International Journal of Coal Geology, vol. 35, no. 1-4, pp. 175-207, 1998.
 I. Gray, "Reservoir engineering in coal seams. Part 1: the physical process of gas storage and movement in coal seams," SPE Reservoir Engineering, vol. 2, no. 1, pp. 28-34, 1987
 W. K. Sawyer, G. W. Paul, and R. A. Schraufnagel, "In development and application of A 3-D coalbed simulator," in Proceedings of the Annual Technical Meeting, Calgary, Alberta, 1990.
 J. P. Seidle and L. G. Huitt, "Experimental measurement of coal matrix shrinkage due to gas desorption and implications for cleat permeability increases," in Proceedings of the International Meeting on Petroleum Engineering, pp. 575-582, Beijing, China, November 1995.
 I. Palmer and J. Mansoori, "How permeability depends on stress and pore pressure in coalbeds: a new model," SPE Reservoir Engineering, vol. 1, no. 6, pp. 539-543, 1998.
 J. Q. Shi and S. Durucan, "Drawdown induced changes in permeability of coalbeds: a new interpretation of the reservoir response to primary recovery," Transport in Porous Media, vol. 56, no. 1, pp. 1-16, 2004.
 X. Cui and R. M. Bustin, "Volumetric strain associated with methane desorption and its impact on coalbed gas production from deep coal seams," AAPG Bulletin, vol. 89, no. 9, pp. 1181-1202, 2005.
 E. P. Robertson and R. L. Christiansen, "A permeability model for coal and other fractured, sorptive-elastic media," in Proceedings of the SPE Eastern Regional Meeting, pp. 89-100, Canton, Ohio, USA, October 2006.
 H. Zhang, J. Liu, and D. Elsworth, "How sorption-induced matrix deformation affects gas flow in coal seams: a new FE model," International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 8, pp. 1226-1236, 2008.
 H.-H. Liu and J. Rutqvist, "A new coal-permeability model: internal swelling stress and fracture-matrix interaction," Transport in Porous Media, vol. 82, no. 1, pp. 157-171, 2010.
 L. D. Connell, M. Lu, and Z. Pan, "An analytical coal permeability model for tri-axial strain and stress conditions," International Journal of Coal Geology, vol. 84, no. 2, pp. 103-114, 2010.
 J. Liu, Z. Chen, D. Elsworth, X. Miao, and X. Mao, "Linking gassorption induced changes in coal permeability to directional strains through a modulus reduction ratio," International Journal ofCoal Geology, vol. 83, no. 1, pp. 21-30, 2010.
 J. Liu, Z. Chen, D. Elsworth, X. Miao, and X. Mao, "Evolution of coal permeability from stress-controlled to displacementcontrolled swelling conditions," Fuel, vol. 90, no. 10, pp. 2987-2997, 2011.
 Z. Pan and L. D. Connell, "A theoretical model for gas adsorption-induced coal swelling," International Journal of Coal Geology, vol. 69, no. 4, pp. 243-252, 2007
 S. Hol and C. J. Spiers, "Competition between adsorption-induced swelling and elastic compression of coal at C[O.sub.2] pressures up to 100 MPa," Journal of the Mechanics and Physics of Solids, vol. 60, no. 11, pp. 1862-1882, 2012.
 S. Durucan, M. Ahsanb, and J. Q. Shia, "Matrix shrinkage and swelling characteristics of European coals," Energy Procedia, vol. 1, no. 1, pp. 3055-3062, 2009.
 S. George J.D. and M. A. Barakat, "The change in effective stress associated with shrinkage from gas deportation in coal," International Journal of Coal Geology, vol. 45, no. 2-3, pp. 105-113, 2001.
 E. P. Robertson and R. L. Christiansen, "Modeling permeability in coal using sorption-induced strain data," in Proceedings of the SPE Annual Technical Conference and Exhibition (ATCE 05),pp. 4133-4142, Dallas, Tex, USA, October 2005.
 H. Kang, X. Zhang, L. Si, Y. Wu, and F. Gao, "In-situ stress measurements and stress distribution characteristics in underground coal mines in China," Engineering Geology, vol. 116, no. 3-4, pp. 333-345, 2010.
 Z. Meng, J. Zhang, and R. Wang, "In-situ stress, pore pressure and stress-dependent permeability in the Southern Qinshui Basin," International Journal of Rock Mechanics and Mining Sciences, vol. 48, no. 1, pp. 122-131, 2011.
 Y. P. Cheng, H. F. Wang, L. Wang, and H. X. Zhou, Theories and Engineering Applications on Coal Mine Gas Control, China University of Mining and Technology Press, Xuzhou, China, 2010.
 L. Wang, Y.-P. Cheng, L. Wang, P.-K. Guo, and W. Li, "Safety line method for the prediction of deep coal-seam gas pressure and its application in coal mines," Safety Science, vol. 50, no. 3, pp. 523-529, 2012.
 Z. X. Sun, W. Zhang, B. Q. Hu, W. J. Li, and T. Y. Pan, "The relationship between the geothermal field characteristics and the distribution of the CBM in the Quishui basin," Chinese Science Bulletin, vol. 50, pp. 93-98, 2005 (Chinese).
 M. A. Biot, "General theory of three-dimensional consolidation," Journal of Applied Physics, vol. 12, no. 2, pp. 155-164, 1941.
 E. Detournay and A. H. D. Cheng, "Fundamentals of poroelasticity," in Comprehensive Rock Engineering: Principles, Practice and Projects, Vol. II, Analysis and Design Method, C. Fairhurst, Ed., pp. 113-171, Pergamon, Oxford, UK, 1993.
 A. Salmachi and M. Haghighi, "Feasibility study of thermally enhanced gas recovery of coal seam gas reservoirs using geothermal resources," Energy and Fuels, vol. 26, no. 8, pp. 5048-5059, 2012.
 D. D. Do, Adsorption Analysis: Equilibria and Kinetics, Imperial College Press , London, UK, 1998.
 S. Harpalani and A. Mitra, "Impact of CO2 injection on flow behavior of coalbed methane reservoirs," Transport in Porous Media, vol. 82, no. 1, pp. 141-156, 2010.
 S. Liu and S. Harpalani, "A new theoretical approach to model sorption induced coal shrinkage or swelling," AAPG Bulletin, vol. 97, no. 7, pp. 1033-1049, 2013.
 X. Cui, R. M. Bustin, and L. Chikatamarla, "Adsorption-induced coal swelling and stress: implications for methane production and acid gas sequestration into coal seams," Journal of Geophysical Research B, vol. 112, no. 10, Article ID B10202, 2007.
 C. O. Karacan, "Swelling-induced volumetric strains internal to a stressed coal associated with CO2 sorption," International Journal of Coal Geology, vol. 72, no. 3-4, pp. 209-220, 2007
 G. K. W. Dawson, S. D. Golding, J. S. Esterle, and P. Massarotto, "Occurrence of minerals within fractures and matrix of selected Bowen and Ruhr Basin coals," International Journal of Coal Geology, vol. 94, pp. 150-166, 2012.
 Z. Chen, Z. Pan, J. Liu, L. D. Connell, and D. Elsworth, "Effect of the effective stress coefficient and sorption-induced strain on the evolution of coal permeability: experimental observations," International Journal of Greenhouse Gas Control, vol. 5, no. 5, pp. 1284-1293, 2011.
 R. Pini, S. Ottiger, L. Burlini, G. Storti, and M. Mazzotti, "Role of adsorption and swelling on the dynamics of gas injection in coal," Journal of Geophysical Research B, vol. 114, no. 4, Article ID B04203, 2009.
 C. R. McKee, A. C. Bumb, and R. A. Koenig, "Stress-dependent permeability and porosity of coal and other geologic formations," SPE Formation Evaluation, vol. 3, no. 1, pp. 81-91, 1988.
 Y. Song, X. M. Zhang, and S. B. Liu, The Basic Theory of Coalbed Methane Geology and Exploration in China, Science Press, Beijing, China, 2012.
 E. T. Brown and E. Hoek, "Trends in relationships between measured in-situ stresses and depth," International Journal of Rock Mechanics and Mining Sciences and, vol. 15, no. 4, pp. 211-215, 1978.
 S. Tao, Y. Wang, D. Tang et al., "Dynamic variation effects of coal permeability during the coalbed methane development process in the Qinshui Basin, China," International Journal of Coal Geology, vol. 93, pp. 16-22, 2012.
 X. Fu, Y. Qin, B. Jiang, W. Wang, and G. Li, "Physical and numerical simulations of permeability of coal reservoirs in central and southern part of the Qinshui Basin, Shanxi," Scientia Geologica Sinica, vol. 38, no. 2, pp. 229-239, 2003 (Chinese).
 G.-Z. Yin, C.-B. Jiang, J. Xu, S.-J. Peng, and W.-P. Li, "Experimental study of thermo-fluid-solid coupling seepage of coal containing gas," Journal of the China Coal Society, vol. 36, no. 9, pp. 1495-1500, 2011.
 S. H. Zhang, S. H. Tang, Y. Wan, Z. C. Li, and J. P. Zhang, "The Migration of CH4 and CO2 in jincheng anthracite coal," Journal of China University of Mining and Technology, vol. 41, no. 5, pp. 770-775, 2012.
 D. H. Bangham and R. E. Franklin, "Thermal expansion of coals and carbonised coals," Transactions of the Faraday Society, vol. 42, pp. B289-B294, 1946.
 S. R. Kelemen and L. M. Kwiatek, "Physical properties of selected block Argonne Premium bituminous coal related to C[O.sub.2], C[H.sub.4], and [N.sub.2] adsorption," International Journal of Coal Geology, vol. 77, no. 1-2, pp. 2-9, 2009.
Pinkun Guo (1,2) and Yuanping Cheng (1,2)
(1) School of Safety Engineering, China University of Mining & Technology, Xuzhou 221116, China
(2) National Engineering Research Center for Coal & Gas Control, China University of Mining & Technology, Xuzhou 221116, China
Correspondence should be addressed to Yuanping Cheng; firstname.lastname@example.org
Received 1 September 2013; Accepted 18 November 2013
Academic Editors: Z. Dai and A. Gudmundsson
TABLE 1: Parameter magnitudes for matching the experimental data. Parameter Value Elastic modulus, E (GPa) 1.12 Poisson's ratio, v 0.26 Bulk modulus, K (MPa) 778 Matrix modulus, [K.sub.m] (MPa) 10340 Constrained axial modulus, M (MPa) 1369 Biot's coefficient, [alpha] 0.925 Initial porosity (for [N.sub.2]), 0.585 [[phi].sub.0] (%) Initial porosity (for C[O.sub.2]), 0.42 [[phi].sub.0] (%) Maximum sorption volume(for [N.sub.2]), 30.05 [alpha] ([m.sup.3]/t) Langmuir parameter (for [N.sub.2]), b 0.07 ([MPa.sup.-1]) Sorption swelling coefficient (for [N.sub.2]), 4.99 x [10.sup.-4] [chi] (t/[m.sup.3]) Langmuir volume (for C[O.sub.2]), [alpha] 55.78 ([m.sup.3]/t) Sorption parameter (for C[O.sub.2]), b 0.48 ([MPa.sup.-1])] Sorption swelling coefficient (for C[O.sub.2]), 9.30 x [10.sup.-4] [chi] (t/[m.sup.3]) Effective coal matrix deformation factor, 0.1723 [f.sub.m] Empirical parameter for P-M model, f 0.1 Matrix compressibility, [gamma] ([MPa.sup.-1]) 9.67 x [10.sup.-5] Fracture compressibility, [C.sub.f] 0.013 ([MPa.sup.-1]) Initial fracture compressibility, [C.sub.0] 0.3422 ([MPa.sup.-1]) Decline rate of fracture compressibility with 2.65 x [10.sup.-14] increasing effective stress, [theta] ([MPa.sup.-1]) TABLE 2: The parameters for predicting permeability in the SQB. Parameters Values Elastic modulus, E (GPa) 3.36 Poisson's ratio, v 0.30 Biot's coefficient, [alpha] 0.80 Coefficient of the thermal deformation, 38.13 x [10.sup.-6] [eta] ([K.sup.-1]) Langmuir volume, a ([m.sup.3]/t) 31.43 Adsorption heat, [DELTA]H (J) 1.70 x [10.sup.4] Adsorption parameter coefficient, [b.sup.0] 8.95 x [10.sup.-3] ([Mpa.sup.-1] x [K.sup.-0.5]) Sorption swelling coefficient, [chi] 2.23 x [10.sup.-4] (t/[m.sup.3]) Initial permeability, [k.sup.0] (mD) 3.00 Initial porosity, [[phi].sup.0] (%) 0.90 Effective coal matrix deformation factor, 0.25 [f.sub.m]
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Guo, Pinkun; Cheng, Yuanping|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2013|
|Previous Article:||Uniformity masks design method based on the shadow matrix for coating materials with different condensation characteristics.|
|Next Article:||Recognition of emotions in Mexican Spanish speech: an approach based on acoustic modelling of emotion-specific vowels.|