Gas Flow Behavior of Nanoscale Pores in Shale Gas Reservoirs

The gas transport in shale nanopores is always one of the major concerns in terms of the development of shale gas reservoirs. In this study, the gas flow regimes in shale nanopores were classified and analyzed according to Knudsen number. Then the gas flow model considering Darcy flow, slip flow, transition flow, molecular free flow and adsorption effect was proposed to evaluate the gas flow behavior in shale nanopores. The result shows that the contributions of Darcy flow, slip flow and transition flow in shale nanopores are reciprocal, and are mainly dominated by pore radius and pressure. The adsorption effect greatly influences the total mass flux. The total mass flux will increase as Langmuir pressure and temperature increase while it will decrease with reservoir pressure and the adsorption thickness. These results can provide insights for a better understanding of gas flow in the shale nanopores so as to optimize the production performance of shale gas reservoirs.


Introduction
In the last few years, shale gas is playing an important role among energy sources and has attracted wide attention [1,2].According to a survey from U.S. Energy Information Administration, 34% of gas production in 2011 in the U.S.A is from shale and the proportion will reach 45% by 2035.Unlike conventional gas reservoirs, shale gas reservoirs are characterized by extremely low porosity, ultra-low permeability and high clay content [3,4].Moreover, the pore structures in shale gas reservoirs are varied and heterogeneous, including organic matter, nonorganic matrix, natural fractures and pore space induced by hydraulic fractures [5][6][7].The giant variation of pores scales makes gas flow in shale gas reservoirs very complex.Thus an understanding of gas flow and transport in shale pores is great significance for gas productivity and for optimizing the hydraulic-fracturing design in shale gas reservoirs.
The gas in shale gas reservoirs includes gas adsorbed on organic matter and inorganic mineral surfaces and free gas in fractures and pores [8,9].Using Scanning Electron Microscopy (SEM) methods Javadpour et al. [10], Locks et al. [11] found that the main pore size of shale matrix is in the range from 1 to 200 nm.Due to the existence of nanopores, the widely used Darcy's law cannot accurately describe gas flow in shales.Based on the work from Beskok and Karniadakis [12], Swami et al. [13] divided gas flow in shales into four regimes, including viscous flow (K n ≤ 0.001), slip flow (0.001 ≤ K n ≤ 0.1), transition flow (0.1 ≤ K n ≤ 10) and Knudsen flow (K n ≥ 10).The gas transport in shale nanopores is a complex combination of all these flow mechanisms, which cannot be described by Darcy's law [14].
In addition, the adsorption and desorption process also exist in the pores with adsorbed gas on organic matter and inorganic mineral surfaces.Some transport models have been developed to understand and quantify gas flow in shale gas reservoirs over the past years.Javadpour [15] and Ozkan et al. [16] proposed a gas transport model in shales including viscous flow and Knudsen flow.Moridis et al. [17] and Freeman et al. [18] presented a very comprehensive review of flow models in shale gas reservoirs considering the four regimes.In order to analyze production dynamics in shale gas reservoirs, an extended dynamic-slippage model was put forward by Clarkson et al. [19,20].Civan et al. [21], Sakhaee and Bryant [22] extended Knudsen diffusion into slippage factor to describe gas transport in shales.Anderson et al. [23] proposed a model to describe the transport mechanisms, which is related to the experiment empirical coefficients.Although there are lots of models proposed to describe gas transport in shale nanopores, a detailed comparison and analysis of different flow patterns is lacking.Moreover, the effects of these parameters on gas flow patterns, such as pressure, temperatures, and adsorption parameters, were also not studied systematically.Therefore, it is extremely necessary to understand gas flow behavior in shale nanopores and the effects on different flow patterns so as to predict gas production and optimize the fracturing treatment for shale gas reservoirs.
In this work, the gas flow regimes were described according to Knudsen number in shale gas reservoirs.The gas flow mathematical models under different flow patterns were analyzed and proposed to understand the flow behavior.Then the effects of pressure and temperature on flow pattern versus pore diameters were analyzed.Besides, a series of sensitivity studies to quantify the effects of these parameters such as adsorption, adsorption thickness, Langmuir pressure, temperature, pressure on the gas flow behavior were performed.This work should provide insights for a better understanding of gas flow behavior in shale nanopores.

Knudsen Number
Compared with conventional gas reservoirs, shale gas reservoirs are characterized by low porosity and ultra-low permeability.With decreasing shale pore radius, the gas transport in shales will change from Darcy flow to other flow patterns.Thus the Darcy flow equation does not correctly describe gas transport in shale gas reservoirs.When the shale pore radius is less than the gas molecular mean free path, the collisions between gas molecules and hole-walls are more frequent than those between gas molecules, and the Knudsen diffusion will occur.In order to divide gas flow patterns, Hadjiconstanttinou [24] used Knudsen number (K n ) to distinguish gas flow regions in the process of micro-nano pipe flow.It can be expressed as follows: where K n is Knudsen number; λ is mean free path of gas molecular; and r is pore radius.Knudsen number is a characteristic parameter of gas flow in different scale channels [25].It is very important to correctly determine the mean free path of gas molecules.The mean free path of gas molecules can be written as: where k B is the Boltzmann constant, k B = 1.3805 × 10 −23 J/K; T is temperature; P is pressure; and δ is collision diameter.

Gas Flow Regimes
According to Knudsen number, the gas flow in shales can be divided four regimes, including continuous flow (known as Darcy flow), slip flow, transition flow and molecular free flow, as shown in Figure 1.At a low Knudsen number (K n < 0.001), the gas transport in shales is Darcy flow.Compared with the pore throat size, the mean free path can be ignored.The molecules' collisions with pore-walls are negligible but the intermolecular collisions are important, so the gas flow is controlled by gas viscosity while Knudsen diffusion can be ignored.The widely used Darcy flow equation and Navier-Stokes equation without slippage both describe gas flow in porous media.With increasing Knudsen number (0.001 < K n < 0.1), the regime is slip flow.For this Knudsen number range, the flow velocity at the pore boundary cannot be ignored and the slippage effect need to be considered in viscous flow.For the mean free path of gas molecules becomes important, the molecules' collisions with pore-walls need to be considered.The Navier-Stokes equation considering slippage effect can describe gas transport in porous media.
in Figure 1.At a low Knudsen number (Kn < 0.001), the gas transport in shales is Darcy flow.Compared with the pore throat size, the mean free path can be ignored.The molecules' collisions with pore-walls are negligible but the intermolecular collisions are important, so the gas flow is controlled by gas viscosity while Knudsen diffusion can be ignored.The widely used Darcy flow equation and Navier-Stokes equation without slippage both describe gas flow in porous media.With increasing Knudsen number (0.001 < Kn < 0.1), the regime is slip flow.For this Knudsen number range, the flow velocity at the pore boundary cannot be ignored and the slippage effect need to be considered in viscous flow.For the mean free path of gas molecules becomes important, the molecules' collisions with pore-walls need to be considered.The Navier-Stokes equation considering slippage effect can describe gas transport in porous media.When Kn ranges from 0.1 to 10, this is the transition flow.The mean free path and pore throat size are the same in scale.The molecules' collisions with pore-walls and the intermolecular collisions are both important.In this condition of flow regime, the continuity assumption and Navier-Stokes equation cannot describe gas flow, which is the most difficult regime.Slip flow and Knudsen diffusion must be considered.When Kn is over 10, this is the molecular free flow.For this range, the mean free path of gas molecules becomes predominant, the molecules' collisions with pore-walls are important and the intermolecular collisions can be ignored.The gas flow is controlled by Knudsen diffusion while gas viscosity are negligible.Thus the Knudsen diffusion equation can be used to describe gas flow in this regime.

Darcy Flow
When Kn is less than 0.001, the gas flow in shale nanopores can be described by the Darcy equation, which is expressed as: where νg is gas velocity; kg is gas permeability; μg is gas viscosity; and ∇ P is pressure gradient.
According to the Hagen-Poiseuille equation, the Darcy permeability can be expressed as: When K n ranges from 0.1 to 10, this is the transition flow.The mean free path and pore throat size are the same in scale.The molecules' collisions with pore-walls and the intermolecular collisions are both important.In this condition of flow regime, the continuity assumption and Navier-Stokes equation cannot describe gas flow, which is the most difficult regime.Slip flow and Knudsen diffusion must be considered.When K n is over 10, this is the molecular free flow.For this range, the mean free path of gas molecules becomes predominant, the molecules' collisions with pore-walls are important and the intermolecular collisions can be ignored.The gas flow is controlled by Knudsen diffusion while gas viscosity are negligible.Thus the Knudsen diffusion equation can be used to describe gas flow in this regime.

Darcy Flow
When K n is less than 0.001, the gas flow in shale nanopores can be described by the Darcy equation, which is expressed as: where ν g is gas velocity; k g is gas permeability; µ g is gas viscosity; and ∇ P is pressure gradient.
According to the Hagen-Poiseuille equation, the Darcy permeability can be expressed as: Substituting Equation (4) into Equation ( 3), the mass flux in Darcy flow can be expressed as: where ρ g is gas density; and r is radius.

Slip Flow
When K n ranges from 0.001 to 0.1, the molecules' collisions with pore-walls are not negligible and there exists slip flow.The gas velocity along pore-walls is not zero, so the slippage effect needs to be considered.An equation proposed by Klinkenberg [26] can be expressed as: where k l is equivalent liquid permeability; c is constant; P is pore pressure; and b k is Klinkenberg coefficient.
For the slip flow regime, Roy et al. [25] presented a similar equation as follows: where σ v is tangential momentum coefficient, σ v ≈ 0.9; and D is pore diameter.Then the permeability in slip flow can be expressed as: Thus the mass flux in the slip flow regime can be written as: where F is slip flow coefficient, F = 1 + 5K n .

Transition Flow
When K n is in the range between 0.1 and 10, the mean free path of gas molecules is the same order as the pore size of porous media.The regime is the transition flow.It not only considers the slippage effect, but also includes the Knudsen diffusion, which is a complex flow process.For the transition flow regime, Ho and Webb [27] used the dusty gas model to describe gas flow, which is expressed as: 0.5 (11) where D k is the Knudsen diffusion coefficient.
Based on the model from Beskok and Karniadakis [28], Civan [29] revised the rarefaction coefficient to describe the gas flow in this regime as follows: α = 1.358 1 + 0.17 K n 0.4348 (13) where α is the rarefaction coefficient; and b is the slip coefficient.
Shi et al. [30] used a ratio coefficient to describe the slippage and diffusion effect in the transient gas flow, and the mass flux can be expressed as: where f is the ratio coefficient; J s and J Kn are the mass flux in the slip and Knudsen diffusion, separately.

Molecular Free Flow
When K n exceeds 10, the gas transport is the molecular free flow.This regime is controlled by Knudsen diffusion.For this regime, the mass flux can be expressed as: where D Kn is the Knudsen diffusivity; and ∇ n i is the gas density gradient of gas component i.
Based on a circular tube, Roy et al. [25] proposed the Knudsen diffusivity which can be written as follows: Thus the mass flux in the molecular free flow can be expressed as:

Gas Mathematical Model
According to the above analytical result, the gas flow in shale gas reservoirs is a complex flow process, including the Darcy flow, slip flow, transition flow and molecular free flow.For the gas transport in shale nanopores, these flow patterns may exist at the same time, and it is not comprehensive to consider one or two flow process.Thus the total mass flux with multi-flow regimes may be expressed as: Substituting Equations ( 5), ( 9), ( 14), ( 18) and ( 19), we obtain: When considering the effect of gas adsorption, the adsorption thickness based on the Langmuir equation can be expressed as: where d is the adsorption thickness; d 0 is the max adsorption thickness; and P L is Langmuir pressure.
And the effective radius in the shale nanopores and Knudsen number considering the adsorption effect can be written as: Thus the total mass flux with multi-flow regimes and adsorption effect can be written as:

Case Study
According to the Equations ( 20) and ( 24), it can be seen that the total mass flux in shale gas flow is influenced by lots of properties and parameters, such as Knudsen number, gas density, pore size, gas viscosity temperature, gas molar mass and pressure gradient.In order to understand the gas flow behavior in shale nanopores, here we will analyze different flow patterns under different conditions.Table 1 presents the mole fraction, mole mass and collision diameter for each component in a typical shale gas reservoir.Table 2 shows the properties of gas flow in shale nanopores.According to the calculations, the ratio mass flux in different flow patterns versus pore size is illustrated in Figure 2. From the result of Figure 2, with increasing pore size, the mass flux ratio of Darcy flow rises while that of slip flow slows down.The transition flow increases firstly, and then decreases after the pore size reaches 10 nm.When the pore size is less than 10 nm, the slip flow is predominant.In the range from 10 nm and 1000 nm, the slip and transient flow dominate the gas transport, and the Darcy flow plays a role in the process.It is noteworthy although the molecular free flow has an effect on gas flow, its mass flux is very small.

Effect of Temperature on Flow Regimes
The relationship between different flow patterns and pore radius for different temperatures from 320 K to 400 K is shown in Figure 4 where it is seen that as temperatures increase, the ratio of the Darcy, slip, transition and molecular free flow changes little at given pore radius.When temperatures are given, with the increasing of pore radius, the Darcy and molecular free flow increases while the slip and transition flow decreases.

Effect of Temperature on Flow Regimes
The relationship between different flow patterns and pore radius for different temperatures from 320 K to 400 K is shown in Figure 4 where it is seen that as temperatures increase, the ratio of the Darcy, slip, transition and molecular free flow changes little at given pore radius.When temperatures are given, with the of pore radius, the Darcy and molecular free flow increases while the slip and transition flow decreases.

Effect of Langmuir Pressure
Four cases of the Langmuir pressure values from 10 MPa to 25 MPa are selected to study the gas flow in shale nanopores, as shown in Figure 5.The result shows the Jad/Jno-ad (Jad and Jnon-ad are the mass flux with adsorption and without adsorption, respectively) ratio will increase with Langmuir pressure when the pore radius is constant.That means that the total mass flux rises with increasing Langmuir pressure.For a given Langmuir pressure, with the pore radius increasing, the Jad/Jno-ad ratio will rise.It says that the total mass flux increases with the pore radius.

Effect of Langmuir Pressure
Four cases of the Langmuir pressure values from 10 MPa to 25 MPa are selected to study the gas flow in shale nanopores, as shown in Figure 5.The result shows the J ad /J no-ad (J ad and J non-ad are the mass flux with adsorption and without adsorption, respectively) ratio will increase with Langmuir pressure when the pore radius is constant.That means that the total mass flux rises with increasing Langmuir pressure.For a given Langmuir pressure, with the pore radius increasing, the J ad /J no-ad ratio will rise.It says that the total mass flux increases with the pore radius.

Effect of Adsorption Thickness
The adsorption thickness values from 0.3 nm to 1 nm are chosen to study the effect of adsorption thickness on gas flow in shale nanopores.The Jad/Jno-ad ratio versus pore radius is illustrated in Figure 6.From the results of Figure 6, it is seen that the Jad/Jno-ad ratio will decrease with adsorption thickness for a given pore size.In other words, with increasing adsorption thickness, the total mass flux will decrease.When pressure is constant, the Jad/Jno-ad ratio will increase with the rising pore radius.It means that the total mass flux will increase as the pore radius increases.

Effect of Pressure
Figure 7 shows the effects in different pressures where four cases from 5 MPa to 40 MPa are considered to study gas flow in shale nanopores.From the result of Figure 7, with increasing pressure in shale nanopores, the Jad/Jno-ad ratio will decrease when the pore radius is constant.The total mass flux decreases as pressure increases.However, for a given pressure the Jad/Jno-ad ratio will increase with the increasing pore radius.It means that the total mass flux will increase as pore size increases when pressure is given.

Effect of Adsorption Thickness
The adsorption thickness values from 0.3 nm to 1 nm are chosen to study the effect of adsorption thickness on gas flow in shale nanopores.The J ad /J no-ad ratio versus pore radius is illustrated in Figure 6.From the results of Figure 6, it is seen that the J ad /J no-ad ratio will decrease with adsorption thickness for a given pore size.In other words, with increasing adsorption thickness, the total mass flux will decrease.When pressure is constant, the J ad /J no-ad ratio will increase with the rising pore radius.It means that the total mass flux will increase as the pore radius increases.

Effect of Adsorption Thickness
The adsorption thickness values from 0.3 nm to 1 nm are chosen to study the effect of adsorption thickness on gas flow in shale nanopores.The Jad/Jno-ad ratio versus pore radius is illustrated in Figure 6.From the results of Figure 6, it is seen that the Jad/Jno-ad ratio will decrease with adsorption thickness for a given pore size.In other words, with increasing adsorption thickness, the total mass flux will decrease.When pressure is constant, the Jad/Jno-ad ratio will increase with the rising pore radius.It means that the total mass flux will increase as the pore radius increases.

Effect of Pressure
Figure 7 shows the effects in different pressures where four cases from 5 MPa to 40 MPa are considered to study gas flow in shale nanopores.From the result of Figure 7, with increasing pressure in shale nanopores, the Jad/Jno-ad ratio will decrease when the pore radius is constant.The total mass flux decreases as pressure increases.However, for a given pressure the Jad/Jno-ad ratio will increase with the increasing pore radius.It means that the total mass flux will increase as pore size increases when pressure is given.

Effect of Pressure
Figure 7 shows the effects in different pressures where four cases from 5 MPa to 40 MPa are considered to study gas flow in shale nanopores.From the result of Figure 7, with increasing pressure in shale nanopores, the J ad /J no-ad ratio will decrease when the pore radius is constant.The total mass flux decreases as pressure increases.However, for a given pressure the J ad /J no-ad ratio will increase with the increasing pore radius.It means that the total mass flux will increase as pore size increases when pressure is given.

Effect of Temperature
The values of different temperatures between 320 K and 400 K are conducted to investigate the effect of adsorption thickness on the gas flow in shale pores.The variation of the Jad/Jno-ad ratio versus pore radius is shown in Figure 8.As illustrated in Figure 8, it is seen that as temperatures increase, the Jad/Jno-ad ratio increases for a given pore radius.It means that the elevated temperature favors total mass flux.When pressure is constant, the Jad/Jno-ad ratio rises with increasing pore radius.That is, the corresponding total mass flux increases as the pore radius increases.

Conclusions
Based on Knudsen number, the gas flow regimes in shale nanopores were classified.Then the gas flow model considering Darcy flow, slip flow, transition flow, molecular free flow and adsorption effect was proposed to describe the gas flow behavior in shale nanopores.The proportion of flow patterns and the effects on flow patterns in shale pores were analyzed.According to the above analytical results, the following conclusions can be drawn: (1) As the pore radius increases, the Darcy flow rises steadily in shale nanopores while the slip flow decreases.The transition flow increases, and then decreases after the pore radius reaches 10 nm.When the pore radius is less than 10 nm, the slip flow is predominant.In the range from 10 nm and 1000 nm, the slip and transient flow dominates the gas transport, and the Darcy flow plays a role in the process.It is noteworthy although the molecular free flow has an effect on gas flow, its mass flux is very small; (2) For a given pore radius,

Effect of Temperature
The values of different temperatures between 320 K and 400 K are conducted to investigate the effect adsorption thickness on the gas flow in shale pores.The variation of the J ad /J no-ad ratio versus pore radius is shown in Figure 8.As illustrated in Figure 8, it is seen that as temperatures increase, the J ad /J no-ad ratio increases for a given pore radius.It means that the elevated temperature favors total mass flux.When pressure is constant, the J ad /J no-ad ratio rises with increasing pore radius.That is, the corresponding total mass flux increases as the pore radius increases.

Effect of Temperature
The values of different temperatures between 320 K and 400 K are conducted to investigate the effect of adsorption thickness on the gas flow in shale pores.The variation of the Jad/Jno-ad ratio versus pore radius is shown in Figure 8.As illustrated in Figure 8, it is seen that as temperatures increase, the Jad/Jno-ad ratio increases for a given pore radius.It means that the elevated temperature favors total mass flux.When pressure is constant, the Jad/Jno-ad ratio rises with increasing pore radius.That is, the corresponding total mass flux increases as the pore radius increases.

Conclusions
Based on Knudsen number, the gas flow regimes in shale nanopores were classified.Then the gas flow model considering Darcy flow, slip flow, transition flow, molecular free flow and adsorption effect was proposed to describe the gas flow behavior in shale nanopores.The proportion of flow patterns and the effects on flow patterns in shale pores were analyzed.According to the above analytical results, the following conclusions can be drawn: (1) As the pore radius increases, the Darcy flow rises steadily in shale nanopores while the slip flow decreases.The transition flow increases, and then decreases after the pore radius reaches 10 nm.When the pore radius is less than 10 nm, the slip flow is predominant.In the range from 10 nm and 1000 nm, the slip and transient flow dominates the gas transport, and the Darcy flow plays a role in the process.It is noteworthy although the molecular free flow has an effect on gas flow, its mass flux is very small; (2) For a given pore radius, with increasing pressure, the Darcy and mole free flow will increase while the slip and transition flow will decrease.The temperature has little effect on the four flows in shale nanopores; (3) The

Conclusions
Based on Knudsen number, the gas flow regimes in shale nanopores were classified.Then the gas flow model considering Darcy flow, slip flow, transition flow, molecular free flow and adsorption effect was proposed to describe the gas flow behavior in shale nanopores.The proportion of flow patterns and the effects on flow patterns in shale pores were analyzed.According to the above analytical results, the following conclusions can be drawn: (1) As the pore radius increases, the Darcy flow rises steadily in shale nanopores while the slip flow decreases.The transition flow increases, and then decreases after the pore radius reaches 10 nm.When the pore radius is less than 10 nm, the slip flow is predominant.In the range from 10 nm and 1000 nm, the slip and transient flow dominates the gas transport, and the Darcy flow plays a role in the process.It is noteworthy although the molecular free flow has an effect on gas flow, its mass flux is very small; (2) For a given pore radius, with increasing pressure, the Darcy and mole free flow will increase while the slip and transition flow will decrease.The temperature has little effect on the four flows in shale nanopores; (3) The adsorption process has a strong influence on the total mass flux in shale pores.When the pore radius is constant, the total mass flux will increase with increasing Langmuir pressure and temperature.However, it will decrease as reservoir pressure and the adsorption thickness increase.

Figure 1 .
Figure 1.Gas flow regimes in shale reservoirs according to Knudsen number.

Figure 1 .
Figure 1.Gas flow regimes in shale reservoirs according to Knudsen number.

Figure 3
Figure 3 presents the effects in different pressures where four cases were calculated by varying pressure from 5 MPa to 40 MPa to study different flow patterns in shale nanopores.From the results of Figure 3, it is known that the flow fraction of Darcy and molecular free flow increases with increasing pressure, but the slip and transition flow decreases.When pressures are constant, with increasing pore radius, the slip and transition flow slows down while the Darcy and molecular free flow gradually rise.

Figure 3 .
Figure 3. Gas flow fraction versus pore radius for various pressures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.

Figure 3
Figure 3 presents the effects in different pressures where four cases were calculated by varying pressure from 5 MPa to 40 MPa to study different flow patterns in shale nanopores.From the results of Figure 3, it is known that the flow fraction of Darcy and molecular free flow increases with increasing pressure, but the slip and transition flow decreases.When pressures are constant, with increasing pore radius, the slip and transition flow slows down while the Darcy and molecular free flow gradually rise.

Figure 3
Figure 3 presents the effects in different pressures where four cases were calculated by varying pressure from 5 MPa to 40 MPa to study different flow patterns in shale nanopores.From the results of Figure 3, it is known that the flow fraction of Darcy and molecular free flow increases with increasing pressure, but the slip and transition flow decreases.When pressures are constant, with increasing pore radius, the slip and transition flow slows down while the Darcy and molecular free flow gradually rise.

Figure 3 .
Figure 3. Gas flow fraction versus pore radius for various pressures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.Figure 3. Gas flow fraction versus pore radius for various pressures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.

Figure 3 .
Figure 3. Gas flow fraction versus pore radius for various pressures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.Figure 3. Gas flow fraction versus pore radius for various pressures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.

Figure 4 .
Figure 4. Gas flow fraction versus pore radius for various temperatures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.

Figure 4 .
Figure 4. Gas flow fraction versus pore radius for various temperatures.(a) The fraction of Darcy flow; (b) the fraction of Slip flow; (c) the fraction of Transition flow; and (d) the fraction of Diffusion flow.

Figure 5 .
Figure 5. J ad /J non-ad ratio versus pore radius for various Langmuir pressures.

Figure 6 .
Figure 6.J ad /J non-ad ratio pore radius for various adsorption thickness.

Figure 7 .
Figure 7. J ad /J non-ad ratio versus pore radius for various pressures.

Figure 8 .
Figure 8. J ad /J non-ad ratio versus pore radius for various temperatures.

Table 1 .
Mole fraction, mole mass and collision diameter for each component in typical shale gas reservoir.

Table 2 .
Properties of gas flow in shale nanopores.