Revaluating coal permeability-gas pressure relation under various gas pressure differential conditions

Identifying changes in coal permeability with gas pressure and accurately codifying mean effective stresses in laboratory samples are crucial in predicting gas-flow behavior in coal reservoirs. Traditionally, coal permeability to gas is assessed using the steady-state method, where the equivalent gas pressure in the coal is indexed to the average of upstream and downstream pressures of the coal, while ignoring the nonlinear gas pressure gradient along the gas flow path. For the flow of a compressible gas, the traditional method consistently underestimates the length/volume-averaged pressure and overestimates mean effective stress. The higher the pressure differential within the sample, the greater the error between the true mean pressure for a compressible fluid and that assumed as the average between upstream and downstream pressures under typical reservoir conditions. A correction coefficient for the compressible fluid pressure asymptotes to approximately 1.3%, representing that the error in mean pressure and effective stress can be on the order of approximately 30%, particularly for highly pressure-sensitive permeabilities and compressibilities, further amplifying errors in evaluated reservoir properties. We utilized this volume-averaged pressure and effective stress to correct permeability and compressibility data reported in the literature. Both the corrected initial permeability and the corrected pore compressibility were found to be smaller than the uncorrected values, due to the underestimation of the true mean fluid pressure, resulting in an overestimation of reservoir permeability if not corrected. The correction coefficient for the initial permeability ranges from 0.6 to 0.1 (reservoir values are only approximately 40% to 90% of laboratory values), while the correction coefficient for pore compressibility remains at approximately 0.75 (reservoir values are only approximately 25% of laboratory value). Errors between the uncorrected and corrected parameters are quantified under various factors, such as confining pressure, gas sorption, and temperature. By analyzing the evolutions of the initial permeability and pore compressibility, the coupling mechanisms of mechanical compression, adsorption swelling, and thermal expansion on the pore structure of the coal can be interpreted. These findings can provide insights that are useful for assessing the sensitivity of coal permeability to gas pressure as truly representative of reservoir conditions.


Introduction
Assessing controls on coal permeability by effective stresses and gas pressures are important in predicting gasflow behavior in coal reservoirs (Liu et al. 2020(Liu et al. , 2016a(Liu et al. , 2015)).Current laboratory measurements of rock materials rely on either steady-state or unsteady-state method.The steady-state method measures the rate of fluid flow rate in a rock core under an applied fluid pressure differential, with the permeability evaluated from Darcy's law (Darcy 1856; http.//www.coret est.com/ autom ated-perme ametedporos imeter.html; Liu et al. 2016b;Li et al. 2009).The fluid flow rate may be measured from either the influent rate or effluent rate (Gensterblum et al. 2014a;Cui et al. 2009;Ghanizadeh et al. 2014;Pei et al. 2019), and under steady state, these should be equivalent.Typically, the outlet flow rate is measured.In lowly permeable materials (< 10 -18 m 2 ) flow rates are difficult to measure using flowmeters and steady conditions must take an extended period to stabilize Thus, pulse decay metods (Brace et al. 1968) utilize the unsteady response where a transient pressure differential is applied at the two ends of a pre-saturated rock core.The resulting permeability is evaluated from the pressure equilibration rate between the upstream and downstream reservoirs or the time to reach equilibrium.
Different from single porosity media, fractured coals contain complex cleat networks and micro-pores in matrix blocks with significantly different porosities and permeabilities (Laubach et al. 1998;Chen et al. 2015;Connell et al. 2016;Gensterblum et al. 2015Gensterblum et al. , 2014b)).The maximum permeability is typically in the high porosity fractures and cleats and is strongly influenced by changes in cleat aperture (Levine 1996;Palmer and Mansoori 1998;Wang et al. 2022;Wang et al. 2021a;Shi and Durucan 2004a;Liu and Rutqvist 2010;Pan and Connell 2011;Liu et al. 2011).Since permeability is strongly effective-stressor pressure-dependent, two end-member deformability boundary conditions are typically applied (Harpalani and Chen 1997;Chen et al. 2011;Meng and Li 2017;Seomoon et al. 2015).One boundary condition is to control bulk deformation as either uniaxial deformation (Geertsma 1966) or fully constrained deformation.The other boundary condition is to control the tri-axial stress state as either a constant confining-pressure or constant effective-stress.Numerous experimental results have demonstrated that if the flow regime of the percolating fluid remains viscous, then coal permeability is positively related to gas pressure (Pini et al. 2009;Harpalani and Schraufnagel 1990;Wang et al. 2011;Kumar et al. 2015).If the gas flow regime shifts from slippage flow to viscous flow due to increasing gas pressure, the resultant coal permeability to first decrease and then partially rebound (Wang et al. 2019).
When the effective stress remains unchanged, the coal permeability varies negatively with increasing pore pressure.
Numerous theoretical models have been developed to interpret the impacts of mechanical compression and adsorptioninduced swelling on coal permeability (Gray 1987;Seidle and Huitt 1995;Palmer and Mansoori 1996;Palmer et al. 2007;Shi andDurucan 2004b, 2005;Cui and Bustin 2005).The pore deformability of coal are significantly impacted by coal compressibility.Since coal reservoirs exist in particular regimes of temperature, moisture, fluid properties, and geo-stress, the coal compressibility coefficient must be quantitatively assessed as influenced by these parameters (Robertson and Christiansen 2005;Shi andDurucan 2004c, 2010;McKee et al. 1988;Palmer 2009;Pan et al. 2010;Harpalani 1999).For example, the water-based or helium-based compressibility coefficient of coal is a positive constant, if the effective stress varies only over a small range.In contrast, the pore compressibility coefficient varies with changes in the pore pressure when the effective stress varies over a wider range.More importantly, gas-adsorption-induced swelling effects further complicate the change in porosity of the coal, whereby the compressibility coefficient may be negative (Liu and Harpalani 2014;Harpalani and Mitra 2010).Therefore, accurate characterization and measurement of the coal permeability-gas pressure relation is a prerequisite for investigating the poromechanical response of coal reservoirs.
In the most frequently used steady-state method, the fluidpressure-gradient term of Darcy's formula may be modified by the impact of fluid compressibility (Chen 1994).For steady flow with an incompressible fluid (i.e.water), coal permeability k is calculated as: where Q is gas flow rate, is gas viscosity, L is flow length within the medium, x is seepage distance and P down and P up denote the downstream and upstream pressures, respectively.The pressure distribution of the incompressible fluid along the flow distance x is For a compressible fluid, the coal permeability in the steady state is determined from: and the nonlinear pressure distribution along the flow distance x (Dana and Skoczylas 1999): (1) Figure 1 compares the two pressure distributions of incompressible and compressible fluids for identical upstream and downstream pressures, identifying linear and nonlinear distributions.Importantly, for a compressible fluid both the central and length/volume-weighted mean pressures are greater than the average of the upstream and downstream pressures.Thus, care must be taken in referencing the measured permeability magnitudes to a pressure or effective stress representative of the experiment -that is not always the average of upstream and downstream bounding magnitudes.
In the following we quantify the representative length/ volume-weighted mean pressure of the fluid under a prescribed pressure differential across a porous medium.The corrected mean pressure is always greater than the uncorrected mean pressure.We use this reevaluation to correct published magnitudes for coal permeability as a function of gas pressure and define the magnitude of the overestimation, if the mean of upstream and downstream pressures is merely used as the reference pressure.Pore compressibility coefficients from the literature are similarly corrected in the present work, and their variations with confining pressure, gas sorption, and temperature are discussed.The use of a corrected pressure method is shown capable of defining a more precise relationship between coal permeability and gas pressure, especially where significant pressure differentials are applied in experiments.

Mean pressure of gas along one-dimensional flow path
The mean pressure of gas along one-dimensional flow path may be determined from the integral of the pressure distributions shown in Fig. 1.For an incompressible fluid (Eq.( 2)), the integral domain (ABH) is the equivalent of a uniform distribution given by the rectangle (ABCD).The mean pressure p mean is the height of this rectangle (0D), as: Clearly, this result demonstrates that the mean pressure of an incompressible fluid is the average of the upstream and downstream pressures.Similarly, the mean pressure of a compressible fluid p eq is the height of the rectangle, as: Figure 2a compares the mean pressure values of compressible fluid and incompressible fluid under various combinations of upstream and downstream pressures.It is clearly seen that when the downstream pressure is 0.1 MPa (atmospheric pressure), the larger the upstream pressure is, the higher the mean pressure of compressible fluid is than the incompressible fluid pressure.If the downstream pressure is 3.0 MPa, the mean pressure of a compressible fluid is approximately equal to that of an incompressible fluid.
In order to correct for this error in evaluating mean fluidpressure for the effect of fluid compressibility, a correction coefficient p eq = p mean is defined as the ratio of the com- pressible fluid pressure p eq = p mean to the incompressible fluid pressure p eq = p mean , as expressed by: Figure 2b displays this relationship between the correction coefficient of gas pressure k eq ∞ = k k ∞ and the upstream pressure for a certain downstream pressure.The higher the downstream pressure, the closer k eq ∞ = k k ∞ is to unity.It should be noted that for the steady-state method of permeability measurement, the downstream pressure of coal sample is typically set to atmospheric, which would be 0.1 MPa (absolute) in this case.Thus,

Permeability observations corrected for gas pressure
It is apparent from Fig. 2a and b that a lower pressure differential can approximately linearize distribution curve of gas pressure.In this paper, therefore, the presented method (Eq.( 6)) is used to evaluate coal permeability for larger pressure differential case, in which case gas slip occurring in lower gas pressure is ignored.We use the McKee et al.
(5) where a measure of unjacketed permeability that is subjected to hydrostatic pressure-equivalent to the pore pressure p being equal to the confining pressure the pore compressibility coefficient.Two cor- rection coefficients for the initial permeability and the pore compressibility are defined as the ratio of the corrected value and the uncorrected value, respectively: k eq ∞ = k k ∞ for the initial permeability (Eq.( 9)) and k eq ∞ = k k ∞ for the pore compressibility (Eq.( 10)):

Impacts of confining pressure
We use three groups of published data correcting the permeability-pore pressure relation to evaluate the impacts of confining pressure.These are for non-sorbing Ar permeability at confining pressures of 10-30 MPa (Han et al. 2010) and for sorbing CH 4 permeability at 4-7 MPa (Dai 2020) and CH 4 permeability at 6-15 MPa (Teng et al. 2021), as listed in Table 1.These are corrected using Eq.(4).From ( 9)  Figs.3a, 4a, and 5a, it is clear that all of the permeabilitycorrected pore-pressure curves move to elevated pore-pressures indicating that coal-permeabilities are always overestimated if the uncorrected pressures are used.Using Eq. ( 8), both the initial permeability and the pore compressibility of the coals are fitted by the corrected and uncorrected gas-pressure relationships, respectively.The fitting parameters are plotted in Figs.3b and c.Both the initial permeability and the pore compressibility of the coal decrease as the confining pressure rising, identifying that the coal matrix skeleton contracts under the pore pressure, in which case the pore pressure is equal to the external stress (or confining pressure).Nevertheless, the two parameters that are fitted by the corrected gas pressure are less than those that are fitted by the uncorrected gas pressure.Consequently, the correction coefficient of the initial permeability increases from 0.4 to 0.7, while the correction coefficient of the pore compressibility remains at approximately 0.75.

Impacts of gas sorption
Gas-sorption effects can significantly change the evolution of coal permeability with gas pressure.CO 2 exhibits a stronger affinity for coal than does CH 4, while the adsorption of N 2 is the weakest.The larger the gas-sorption capacity is, the greater the decreasing amplitude of coal permeability becomes (Meng et al. 2015;Wang et al. 2017Wang et al. , 2021b;;Feng 2021).The pressure-corrected relationships between permeability and gas pressure from this dataset are shown in Figs.6a, b, and c.Since the mean pressure of the gas will be underestimated for a compressible gas, the corresponding coal permeability will be overestimated.Figure 7a compares the corrected initial permeability with the uncorrected permeability.For the cases using He, N 2 , CO 2 , and CH 4 , the corrected initial permeability is found to be generally lower than the uncorrected one.The relationship between the magnitude of the initial permeability and the Fig. 3 Relationship between coal permeability and Ar pressure for different confining pressures (Han et al. 2010).a Coal permeability corrected for Ar pressure; b Fitted initial permeabilities corrected for Ar pressure; c Pore compressibilityies corrected for Ar pressures gas absorbability is not clear, which is likely due both to the heterogeneity of different coals and to the applied stress states.Nevertheless, the correction coefficient for the initial permeability remains in the range between 0.6 and 0.8, as shown in Fig. 7b.Similarly, the corrected pore compressibility is also lower than the uncorrected compressibility, as shown in Fig. 8a.Although the pore compressibility of different coals is related to the gas species, the relevant correction coefficients are approximately 0.75 (Fig. 8b).

Impacts of thermally-induced expansion
The heating of coal typically results in an inhomogeneous expansion and the cracking between minerals, organic matter, and inorganic matter (Heuze 1983;Wong and Brace 1979) as well as to gas desorption (Charrière et al. 2010;Deishad et al. 2009) with these multi-physical processes potentially affecting permeability.From the literature (Gao et al. 2021;Teng et al. 2021), coal permeability can clearly be seen to decrease with increasing ambient temperature.By correcting the experimental dataset, Figs.9a and 10a indicate that all of the coal-permeability-corrected gas-pressure curves translate to higher equivalent pore-pressures.Thus, if the pressures remain uncorrected then the corresponding coal permeability will be overestimated.As shown in Figs.9b and 10b, both the corrected initial permeability and the corrected pore compressibility are lower than the uncorrected magnitudes.The correction coefficients of the initial permeability and the pore compressibility are approximately 0.5 for Fig. 9c and 0.75 for Fig. 10c, respectively.

Error analysis of the permeability-to-mean-gas-pressure relation
The aforementioned data of Figs. 3, 4, and 5 demonstrate that coal permeability will rise exponentially with increasing gas pressure, if the confining pressure remain unchanged at a certain isothermal condition.As shown in Fig. 11a, because the corrected permeability curves translate rightward to higher pressures, the absolute error between the corrected permeability and the uncorrected permeability increases correspondingly.
Based on the literature data (Dai 2020), Fig. 11b compares the relative error between the corrected permeability and the uncorrected permeability under four confining-pressure conditions ranging from 4 to 7 MPa.For any of four confining pressure condition, the rising pore pressure can cause the relative error to first increases and then decrease.At the case of a confining pressure of 4 MPa, when the pore pressure ranges below 1.5 MPa, the error is within 10%, but when the pore pressure increases to 2.5 MPa, the error peaks at 200%.
Under various isothermal conditions from Figs. 9a and 10a, it is reasonably assumed that coal permeability can be defined as a logarithmic function of gas pressure, whereby the incremental permeability decreases with the gas pressure rising.As a result, Fig. 12a demonstrates that the absolute error between the corrected permeability and the uncorrected permeability also decreases with the gas pressure rising.Changes in the corresponding relative error under four different temperatures are shown in Fig. 12b.As can be seen, the relative error first increases and then peaks from between 10% at 0.5 MPa under 20 °C and 15% at 0.7 MPa under 40 °C.It is found that the higher the isothermal condition is, the greater the relative error becomes.

Implications for pore deformation due to coalgas interaction
It is well known that the effective stress is defined by the portion of the external stress which is supported by the solid skeleton of the porous medium, while the remaining external stress carried by the pore fluid.Any change in effective stress reduces the pore radius and hence the permeability of a single pore in an impermeable medium.Unlike the singlepore medium, the matrix blocks in coal are permeable and can usually be simplified as an assemblage of discrete grains that contain connected and closed micropores.When gas  flows into the coal, the gas pressure resists not only a part of the external stress, but also the internal stresses among the grains in the matrix.Fjaer et.al. (2008) argue that as the pore pressure rises, a greater proportion of the external stress is counteracted by the pore pressure.The remaining external stress that acts on the matrix accordingly decreases, which can cause the grains in the matrix to expand.This argument can be demonstrated by the observed decrease in coal permeability with increasing pore pressure under the constant effective stress (Seomoon et al. 2015).In fact, the stress state of the presented initial permeability is an extreme condition in the constant-effective-stress case-namely zero effective stress.In this case, Han's (2010) work demonstrated that the initial permeability of coal to Ar decreases with increasing confining pressure, which is consistent with laboratory observations under non-zero effective-stress cases.This finding indicates that coal-matrix swelling not only enlarges the overall size of the coal, but also consequently reduces the fracture (or void) volume.This interplay between coal matrix and facture is shown schematically in Fig. 13.In addition to the mechanical compression shown here, gas adsorption may also cause the coal matrix to swell, whereby a part of the matrix-swelling strain can also reduce the void volume and narrow fractures, as was demonstrated in Wang's observations (Wang et al. 2021b).
Given that molecular motion in gases is sensitive to temperature, the thermal impacts of gas-bearing coal on permeability are significantly controlled via competition between the gas-desorption-induced shrinkage of the matrix and the thermal expansion of the matrix.Figure 10b reveals that the initial permeability of the coal first increases and then decreases if temperatures are raised from 20 to 50 °C.We postulate that the permeability enhancement at lower temperature is dominated by the impacts of gas-desorptioninduced shrinkage of the matrix, while the subsequent reduction might result from the thermal expansion of the matrix at higher temperature.In contrast, Teng's (Teng et al. 2021) observations suggest that the initial permeability of State I denotes an initial stage, in which case the confining pressure σ c is higher than the pore pressure p.State II denotes an equilibrium state, in which the pore pressure p is increased to the confining pressure σ c .As a result, the pore volume is reduced due to expansion of the coal volume coal declines with temperature rising from 45 to 85 °C.It seems probable that the thermal expansion of the matrix (linear with temperature and unbounded) could exceed the desorption-induced shrinkage of the matrices (nonlinear and bounded by a finite mass of desorbing gas), in which case, a portion of the expanding matrix could decrease the fracture volume.

Conclusions
We identify the need to define effectives stresses in samples via a volume-averaged pore pressure representative of the flow path.This is necessary to correctly index the effective-stress-dependent and pressure-dependent nature of permeability, pore compressibility and sorption to the correct reference stresses and allow laboratory observations to be matched to reservoir data.We develop a straightforward method to calculate the mean (volume-averaged) pressure of compressible fluids transiting a porous medium.We apply this method to existing data representative of laboratory observations of permeability, pore compressibility and sorption conducted at equivalent reservoir conditions and show that reference pressures may be in error by 0% to approximately 30%.Since permeability is a sensitive and nonlinear function of effective stress, permeabilities may be overestimated by 40%-90% and pore compressibilities by approximately 25%-thus the effect is significant.We draw the following conclusions: (1) The higher the pressure differential between the two ends of a porous medium, the greater the mean pressure of the compressible fluid relative to the average for an incompressible fluid.A correction coefficient for the mean pressure for the compressible fluid pressure is defined relative to the assumed mean pressure of the incompressible fluid (i.e. the average of upstream and downstream pressures).Where the effluent/downstream pressure is atmospheric the correction coefficient asymptotes to approximately 1.3 for reasonable reservoir conditions and gases (i.e.upstream pressure > 3 MPa)-defining a 30% error.(2) Laboratory measured permeabilities and pore compressibilities present in the literature data were corrected with the method.If the mean pressure of the gas under a certain pressure differential is not corrected for compressibility, then the laboratory-interpreted permeability will overestimate the real permeability in the reservoir.If the confining pressure remain unchanged, the absolute error between the corrected permeability and the uncorrected permeability begin to accumulate, as the pore pressure increases; the relative error of the permeability under a certain confining pressure first increases and then decreases.As the isothermal level is elevated, the incremental permeability decreases with the gas pressure rising.The absolute error between the corrected permeability and the uncorrected permeability also decreases with increasing gas pressure.
The relative error first increases with increasing gas pressure and then peaks from between 10% at 0.5 MPa under 20 °C and 15% at 0.7 MPa under 40 °C.The higher the coal temperature is, the greater the relative error becomes.(3) Poroelastic parameters of the McKee (McKee et al. 1988) model are fitted from both the corrected and uncorrected datasets.Both the corrected initial permeability and the corrected pore compressibility were found to be smaller than the uncorrected magnitudes, due to the underestimation of true mean fluid pressure.
The correction coefficient for the initial permeability ranges widely but in the range approximately 0.6 to 0.1 (approximately 40% to 90%), while the correction coefficient for pore compressibility remains at approximately 0.75 (approximately 25%).( 4) According to the fitted the initial permeability and the pore compressibility of coal under different conditions, such as, confining pressure, gas adsorption and temperature, it is reasonably demonstrated that the swelling or expansion of the coal matrix block can not only increase the overall size of the coal, but also reduce the coal fracture aperture to a certain extent.This finding can provide a critical insight into evaluating the poroelastic behavior of a coal.
Page 13 of 14 55 need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.

Fig. 1
Fig. 1 Schematic of fluid-pressure-distribution along a flow path for compressible and incompressible fluids

Fig. 2
Fig. 2 Relationship between mean incompressible and compressible fluid pressures with varying upstream and downstream pressures.a Mean pressure of incompressible and compressible fluids; b Change in correction-coefficient

Fig. 4
Fig. 4 Coal permeability-CH 4 pressure curves under different confining pressures (Dai 2020).a Coal permeability corrected for CH 4 pressure; b Initial permeability corrected for CH 4 pressure; c Pore compressibility corrected for CH 4 pressure

Fig.Fig. 7
Fig. CH 4 pressure-dependent coal permeability under different confining pressures (Teng et al. 2021).a Coal permeability corrected for CH 4 pressure; b Initial permeability corrected for CH 4 pressure; c Pore compressibility corrected for CH 4 pressure

Fig. 8 Fig. 9
Fig. 8 Comparison between uncorrected pore compressibility and corrected compressibility for different gases.a Corrected pore compressibility; b Correction coefficient for pore compressibility

Fig. 10
Fig. 10 Correlations permeability with CH 4 pressure for variable temperatures (Gao et al. 2021).a Coal permeability corrected for CH 4 pressure; b Fitted initial permeability for corrected CH 4 pressure; c Fitted pore compressibility for corrected CH 4 pressure

Fig. 12 Fig. 13
Fig. 12Error analysis for the relationship between coal permeability and corrected gas pressure under different isothermal conditions(Gao et al. 2021).a Relative error calculation for permeability; b Relative error of permeability with gas pressure and temperature

Table 1
Compilation of published experimental data