Transport Behavior of Methane Conﬁned in Nanoscale Porous Media: Impact of Pore Evolution Characteristics

: As a key technical aspect contributing to shale gas development, nanoconﬁned methane ﬂow behavior has received tremendous research interest, which remains challenging to understand clearly. The majority of previous contributions put emphasis on the mechanism model for methane conﬁned in a single nanopore; at the same time, the other part focusing on an upscaling approach fails to capture the spatial pore-network characteristics as well as the way to assign pressure conditions to methane ﬂow behavior. In light of the current knowledge gap, pore-network modeling is performed, in which a pore coordination number, indicating the maximum pores a speciﬁed pore can connect, gas ﬂow regimes classiﬁed by Knudsen numbers, as well as different assigned pressure conditions, are incorporated. Notably, the pore-network modeling is completely self-coded, which is more ﬂexible in adjusting the spatial features of a constructed pore network than a traditional one. In this paper, the nanoconﬁned methane ﬂow behavior is elaborated ﬁrst, then the pore network modeling method based on the mass conservation principle is introduced for upscaling, and in-depth analysis is implemented after that. Results show that (a) as for porous media with pore sizes ranging from 5~80 nm, dramatic advancement on apparent gas permeability takes place while pressure is less than 1 MPa; (b) apparent gas permeability evaluated at a speciﬁed pressure shall be underestimated by as much as 31.1% on average compared with that under the pressure-difference condition; (c) both a large pore size and a high coordination number are beneﬁcial for strong gas ﬂow capacity through nanoscale porous media, and the rising ratio can reach about 6 times by altering the coordination number from 3 to 7, which is quantiﬁed and presented for the ﬁrst time.


Introduction
The development of unconventional gas reservoirs, such as shale gas, coal seam gas, and tight gas, is a feasible approach to reduce coal combustion as well as CO 2 emissions [1,2]. As reported by the EIA at the end of 2020, the geological reserve of shale gas reaches 317.76 trillion cubic feet, accounting for 70.68% of total proved natural gas reserves in the United States. As a result, both the petroleum industry and the ecological environment will benefit from the successful exploitation of shale gas [3,4]. With the combination of hydraulic fracturing technology and multi-lateral horizontal wells [5,6], gas production performance from ultra-tight shale formations with rich nanoscale pores has been surprisingly enhanced in the past few decades. At the same time, theoretical breakthroughs in shale gas production, particularly for methane confined in nanoscale pore space, are still challenging. Although profound knowledge of multi-scale gas flow physics in shale is the best guarantee for shale Impact of spatial structure of pore network is not elaborated.
microscopic oil flow mechanism and volumetric toc.
The way pressure condition is assigned on flow has not been revealed.
Wang and Sheng (2019) [27] Gas-water flow in shale porous media. nanopore tortuosity, pore size distribution of organic matter and inorganic matter.
Impact of spatial structure of pore network is not elaborated.
The way pressure condition is assigned on flow has not been revealed. Yi et al. (2021) [29] Gas transport in nanoscale porous media.
microscopic gas flow mechanism, spatial structure of pore network.
The way pressure condition is assigned on flow has not been revealed.
Based on the above analysis, although current research on PNM has achieved fairly accurate prediction of fluid flow in nanoscale porous media, several key shortcomings remain to be overcome. First, the effect of the spatial structure of the pore network, such as the pore size distribution (PSD) and coordination number, on the gas flow behavior has not been investigated completely. Furthermore, the most common method for determining fluid flow capacity is to assign a typical pressure with a negligible pressure difference, while the presence of a pressure difference during the depressurization development process is obvious. As a result, the frequently utilized way to assign pressure condition is different from that in field condition, while the associated impact resulting from the difference is unclear. In light of current research updates, research is performed to try to fill the above-mentioned knowledge gap. It should be underlined that the pore network modeling in this work is totally self-coded, with great flexibility to alter the spatial structure of the pore network, which facilitates identifying the influences of each parameter. Both the PSD and the coordination number of gas transport in nanoscale porous media are investigated. Moreover, the way of assigning pressure conditions is discussed, and a detailed analysis is provided. Therefore, the novel technical part of this research basically resides in the incorporation of spatial pore-network structure and sensible pressure conditions for realistic field exploitation.
The whole paper context is arranged in the following way. The physical model making clear the essential advantage of pore network modeling over the traditional mechanism model or capillary bundle model is presented first. Then, reliable formulas reproducing the nanoconfined gas flow behavior are introduced, and the primary concepts establishing the pore network modeling are briefly provided, including the part to generate concrete spatial structure and that for flow capacity calculation. After that, gas transport through nanoscale porous media is seriously discussed. Finally, several conclusions are drawn.

Physical Model
In order to reach a good understanding of gas flow in the shale matrix, which is ultratight and rich in nanopores with pore sizes ranging from several nanometers to hundreds of nanometers, typical procedures to achieve the upscaling objective are presented in Figure 1. Generally, research efforts start with the primary mechanism model, representing the gas flow behavior in a single nanopore. As for this technical part, molecular dynamics (MD), theoretical modeling, as well as advanced experimental measurement, are used to investigate the nanoconfined gas transport behavior. With the intention to bridge nanoconfined flow and reservoir simulation, it entails an effective tool to transition the mechanism model to that suitable for shale core scale, which can be applied as reservoir permeability. At the very beginning, the classic capillary bundle model is utilized, consisting of a certain number of nanopores with various pore radii, pore lengths, and tortuosities. Although the capillary bundle model is a fairly direct upscaling method, there is an evident difference between the capillary bundle model and the realistic shale matrix, limiting its accuracy on gas flow in porous media. In the capillary bundle model, gas flux in a specific nanopore does not affect that in neighboring nanopores, which seriously conflicts with the actual case. In light of this, pore network modeling is developed, which generally includes two parts. The first part is the construction of the pore network, which is a spatial structure composed of massive nanopores with assigned pore network information, including pore size distribution, coordination number, and so on. Then, the second part is the simulation of gas flow in porous media, which is inherently dependent on the mechanism model and the mass conservation principle. From the perspective of flow physics, pore network modeling is evidently superior to the capillary bundle model, which is able to capture gas flux exchange between nanopores and has a more complicated spatial structure.
Processes 2022, 10, x FOR PEER REVIEW 4 of 14 through nanoscale porous media is seriously discussed. Finally, several conclusions are drawn.

Physical Model
In order to reach a good understanding of gas flow in the shale matrix, which is ultratight and rich in nanopores with pore sizes ranging from several nanometers to hundreds of nanometers, typical procedures to achieve the upscaling objective are presented in Figure 1. Generally, research efforts start with the primary mechanism model, representing the gas flow behavior in a single nanopore. As for this technical part, molecular dynamics (MD), theoretical modeling, as well as advanced experimental measurement, are used to investigate the nanoconfined gas transport behavior. With the intention to bridge nanoconfined flow and reservoir simulation, it entails an effective tool to transition the mechanism model to that suitable for shale core scale, which can be applied as reservoir permeability. At the very beginning, the classic capillary bundle model is utilized, consisting of a certain number of nanopores with various pore radii, pore lengths, and tortuosities. Although the capillary bundle model is a fairly direct upscaling method, there is an evident difference between the capillary bundle model and the realistic shale matrix, limiting its accuracy on gas flow in porous media. In the capillary bundle model, gas flux in a specific nanopore does not affect that in neighboring nanopores, which seriously conflicts with the actual case. In light of this, pore network modeling is developed, which generally includes two parts. The first part is the construction of the pore network, which is a spatial structure composed of massive nanopores with assigned pore network information, including pore size distribution, coordination number, and so on. Then, the second part is the simulation of gas flow in porous media, which is inherently dependent on the mechanism model and the mass conservation principle. From the perspective of flow physics, pore network modeling is evidently superior to the capillary bundle model, which is able to capture gas flux exchange between nanopores and has a more complicated spatial structure. In this work, self-coded pore network modeling is utilized so that pore size distribution, coordination number, and skeleton porosity can be accurately and conveniently controlled. Moreover, microscopic flow physics is fully coupled in this pore network modeling, and gas mass flux between nanopores is considered to ensure the mass conservation of the whole system by an efficient algorithm. Notably, the gas flow capacity through porous media is evaluated in a steady state, which is suitable to describe gas permeability in shale matrix during the long-term depressurization development process. In the modeling of this research, the calculation node, representing pores, is over 200, and the number of connecting throats is over 300, which mitigate the uncertainty caused by a discontinuous medium and ensure the reliability of the obtained results. In this work, self-coded pore network modeling is utilized so that pore size distribution, coordination number, and skeleton porosity can be accurately and conveniently controlled. Moreover, microscopic flow physics is fully coupled in this pore network modeling, and gas mass flux between nanopores is considered to ensure the mass conservation of the whole system by an efficient algorithm. Notably, the gas flow capacity through porous media is evaluated in a steady state, which is suitable to describe gas permeability in shale matrix during the long-term depressurization development process. In the modeling of this research, the calculation node, representing pores, is over 200, and the number of connecting throats is over 300, which mitigate the uncertainty caused by a discontinuous medium and ensure the reliability of the obtained results.

Nanoconfined Flow Behavior
The mathematic part of this research starts with the characterization of nanoconfined gas flow behavior, which is apparently different from regular fluid flow described by the classic Navier-Stokes equation [30]. The origin of the difference is caused by the Processes 2022, 10, 2746 5 of 13 nanoconfinement effect induced by the greatly reduced pore size. From the microscopic viewpoint, the molecular flow velocity near the wall is close to 0 as in the case of bulk scale, and its movability is considered negligible. In contrast, when collisions between gas molecules and the wall surface become dominant, molecular slip takes place at the boundary, manifesting that the movability of molecules near the wall cannot be neglected. It is experimentally reported that nanoconfined flow capacity can reach as great as several orders of magnitude as evaluated by the no-slip equation [31,32].
A number of theoretical models have been proposed to elucidate the discrepancy compared with regular equations for gas flow, and four flow regimes are summarized, consisting of continuum flow, slip flow, transition flow, and free molecular diffusion. The flow regime is heavily dependent on the surrounding pressure and temperature conditions, which are identified by the Knudsen number, with the following formula: where Kn is the Knudsen number, dimensionless; λ is the mean free path, representing the distance a gas molecule can move between two adjacent collision events, its calculation formula is Equation (2), m; r is the pore radius, m; µ is gas viscosity, cp; p is the pressure, MPa; R is the gas universal constant, which is equal to 8.314 J/mol/K; Z is the gas compressibility factor, dimensionless; T is the temperature, K; M is the gas molar mass, 16 g/mol. After carefully checking existing models developed for nanoconfined gas flow, the primary constitutive equation for the majority of current models describing bulk-gas flow in nanopores is Beskok and Karniadakis (1999), covering the entire Knudsen number range. Thus, the reliable model is adopted in this work, notably the stress dependence and surface diffusion of the adsorption phase, which are not considered currently but can be incorporated and discussed in the future. Concrete formulas for gas transport in nanopores considering gas slip are provided below.
where α is the rarefaction coefficient, dimensionless; α o is the rarefaction coefficient with infinite Kn, dimensionless; α 1 and β are constants, which are 1.19 and 0.4, respectively. On the basis of Equations (1)-(4), the apparent permeability for gas transport through a nanopore can be calculated. Furthermore, in order to implement the conservation principle of gas mass rate, it is crucial to adopt a suitable approach to obtaining gas physical properties, such as viscosity and density. The classic PVT equation [33,34] is provided below as a beginning point.
where v is the molar volume, m 3 /mol. Notably, a and b are basic PVT parameters, which are described as a function of critical pressure and critical temperature in the bulk state previously. The authors have noticed variations in gas critical properties under the nanoconfinement effect, also its effect on gas density can result in errors as high as 23% [35]. Thus, the shift of critical properties shall be incorporated in the calculation of gas density. Formulas describing a and b are below.
Processes 2022, 10, 2746 where P c , T c are the critical nanoconfined gas pressure and temperature, concrete formulas suggesting the shift compared with the original critical properties may refer to our previous research. T r is the reduced temperature, dimensionless, and ω is the methane acentric factor, dimensionless. Then, the gas density can be sensibly obtained with the given pressure and temperature conditions. After calculating sensible gas density under the nanoconfinement effect, the formulas below are available for gas viscosity [36,37], which is described as a function of gas density.
where ρ is gas density, kg/m 3 ; µ is gas viscosity, cp. Then, for a single nanopore with gas transport, the mass flow rate has the following expression.
where ∆p is the pressure difference between the two ends of the nanopore, Pa; l is the length of nanopore, m; q is the mass flow rate, kg/s.

Pore Network Modeling
First, it is necessary to elaborate pore network construction, in which numerous nanopores with various pore sizes and lengths are connected; the number of nanopores exceeds 200 to mitigate the uncertainty caused by discontinuous medium. Notably, pore size distribution, porosity, and coordination of the ultimate pore network are controlled by a powerful algorithm, ensuring that the generated spatial parameters can have excellent matches against the assigned targets. At the beginning of construction, an assigned number of nodes, suggesting the pores in a tight matrix, are placed randomly in a specific space. Then, the interpore connections take place, suggesting that the nodes find their neighboring nodes and try to connect them at a certain distance. At the same time, the formation of throats is subject to a coordination number, which indicates the average number of throats a node should connect. After that, the basic pore network is established, and a fine algorithm shall complete the assignments of pore size to all throats and pores.
The inherent physics driving proper modeling is the mass conservation of gas flow rates. As for the case of i-th node in the pore network, it may have several connections, like n, indicating that the inflow and outflow of the i-th node must happen in the n connections. Furthermore, for a stable state, the mass flow rates from inflow and outflow should be equal after the whole balance is reached. On the basis of Equation (13), which describes the mass flow rate in a single connection, then the summation of the mass flow rates from n connections has the following expression.
Once the pressure of two connected nodes is given, the mass flow rate between them can be calculated, in which gas density and viscosity are evaluated at the arithmetic mean pressure of the two connected nodes. Notably, for the whole pore network, where the total number of nodes is m, as the pressure at each node is unknown from the beginning of the calculation, there exist m unknown pressure values. Moreover, for each node, there exists Equation (14) as a function of its unknown pressure; therefore, m equations exist for m unknown values, indicating the solutions of all nodes can be obtained properly.
After all the unknown pressures are calculated, the nodes connecting to the boundary nodes at the exit side of the whole pore network are identified, and the mass flow rates from these nodes are summarized and regarded as the total gas mass rate. On the basis of the total rate, the apparent permeability of the entire porous medium is obtained using the equation below [38][39][40], in which gas density and viscosity are evaluated at the mean pressures of inlet pressure and outlet pressure.
where q t is the total gas mass rate through the pore network, kg/s; K app is apparent gas permeability in the constructed pore network, nD.

Results and Discussion
In this research, the slip feature and spatial characteristics of nanoscale porous media on apparent gas permeability at the core scale are focused. According to a typical tight gas reservoir located in Xinjiang, China, the pore size in this work falls in the range from 5~80 nm, and the pressure ranges of 0.1~50 MPa [41], covering the entire pressure variation during actual production.

Rarefaction Phenomenon
The unique feature of nanoconfined flow is the molecular slip at the boundary; other than that, gas flow in nanoscale porous media is affected by network structure. In this section, we put attention on the slip at the boundary, also called the rarefaction phenomenon, the main mechanism contributing to the enhancement of gas flow capacity. The first step in performing pore network modeling is constructing a fancy network with desirable key parameters, like pore size distribution and coordination number. At this point, in order to remove the influence induced by network characteristics, relevant parameters are assigned in advance.
As presented in Figure 2a, from the comparison between the assigned pore size distribution and the actual one generated by the self-coded algorithm provided, it can be demonstrated the algorithm is able to reach a good agreement with the assigned goal. From Figure 2b, it can be observed that the pore size in the network ranges from 5 to 80 nm, also different pores connect in such a way that the coordination number is 3. The gas flow through nanoscale porous media is simulated under the condition that the pressure difference between the inlet and outlet is set at 0.001 MPa. In this case, the gas flow capacity at a specified pressure condition can be mimicked, which makes it fairly clear to explore the rarefaction phenomenon with regards to pressure variation. In Figure  3a, the red dotted line represents the apparent permeability for gas flow in porous media The gas flow through nanoscale porous media is simulated under the condition that the pressure difference between the inlet and outlet is set at 0.001 MPa. In this case, the gas flow capacity at a specified pressure condition can be mimicked, which makes it fairly clear to explore the rarefaction phenomenon with regards to pressure variation. In Figure 3a, the red dotted line represents the apparent permeability for gas flow in porous media without considering the slip boundary, which remains a constant value and is insensitive to pressure change as indicated. The phenomenon that the permeability remains unchanged is in line with the classic Darcy's law, suggesting that the permeability is an inherent property of reservoirs and shall not change with pressure variation. After considering boundary slip taking place under the nanoconfinement effect, the traditional principle about unchanged permeability breaks down. As illustrated in Figure 3a, the apparent permeability (blue line) shows an evident negative trend with the pressure rise; the apparent permeability can reach 8103 nD at 0.1 MPa while declining to 626 nD at 10 MPa, close to the value without considering boundary slip. Firstly, at a low-pressure range, the rarefaction effect plays as a greatly positive role in enhancing the apparent permeability. The gas velocity in nanopores is greatly advanced due to gas slippage at boundary, in which the molecules close to the boundary gain a certain degree of mobility. Meanwhile, the sharp decrease in apparent permeability with increasing pressure indicates that high pressure will inhibit the rarefaction effect greatly. As can be observed in Figure 3a, the apparent permeability at 50 MPa is exactly equal to the value from the classic theory, suggesting the gas slip vanishes at high pressure. The gas flow through nanoscale porous media is simulated under the condition that the pressure difference between the inlet and outlet is set at 0.001 MPa. In this case, the gas flow capacity at a specified pressure condition can be mimicked, which makes it fairly clear to explore the rarefaction phenomenon with regards to pressure variation. In Figure  3a, the red dotted line represents the apparent permeability for gas flow in porous media without considering the slip boundary, which remains a constant value and is insensitive to pressure change as indicated. The phenomenon that the permeability remains unchanged is in line with the classic Darcy's law, suggesting that the permeability is an inherent property of reservoirs and shall not change with pressure variation. After considering boundary slip taking place under the nanoconfinement effect, the traditional principle about unchanged permeability breaks down. As illustrated in Figure 3a, the apparent permeability (blue line) shows an evident negative trend with the pressure rise; the apparent permeability can reach 8103 nD at 0.1 MPa while declining to 626 nD at 10 MPa, close to the value without considering boundary slip. Firstly, at a low-pressure range, the rarefaction effect plays as a greatly positive role in enhancing the apparent permeability. The gas velocity in nanopores is greatly advanced due to gas slippage at boundary, in which the molecules close to the boundary gain a certain degree of mobility. Meanwhile, the sharp decrease in apparent permeability with increasing pressure indicates that high pressure will inhibit the rarefaction effect greatly. As can be observed in Figure 3a, the apparent permeability at 50 MPa is exactly equal to the value from the classic theory, suggesting the gas slip vanishes at high pressure. Notably, the assigned pressure condition of negligible difference between the inlet and outlet pressure failed to fit the realistic pressure condition during the depressurization development process. The pressure difference may make a difference on gas flow capacity, as nanoconfined gas flow is sensitive to pressure variation as indicated in Figure 3a. The spatial pressure distribution with the given inlet pressure (10 MPa) and outlet pressure (9.999 MPa) is illustrated in Figure 3b. In order to address this issue, in the following investigation, outlet pressure is assigned as a constant of 0.001 MPa, at the same time, the inlet pressure ranges from 0.1~50 MPa. As presented in Figure 4a, similar to that under the specified pressure condition, apparent permeability rises with the pressure reduction, and the rise amplitude increases in the low-pressure condition. Compared with the specified pressure condition that the gas flow regime in the whole network is almost the same, the large pressure difference in this case will result in several flow regimes in the network. For example, the continuum flow dominates the region close to the inlet side of the network, while slip flow or transition flow is the main flow mechanism for the region close to the outlet side. As a result, the apparent gas permeability is controlled by multiple flow regimes in the network in this pressure condition rather than a simple single flow regime with a negligible pressure difference. time, the inlet pressure ranges from 0.1~50 MPa. As presented in Figure 4a, similar to that under the specified pressure condition, apparent permeability rises with the pressure reduction, and the rise amplitude increases in the low-pressure condition. Compared with the specified pressure condition that the gas flow regime in the whole network is almost the same, the large pressure difference in this case will result in several flow regimes in the network. For example, the continuum flow dominates the region close to the inlet side of the network, while slip flow or transition flow is the main flow mechanism for the region close to the outlet side. As a result, the apparent gas permeability is controlled by multiple flow regimes in the network in this pressure condition rather than a simple single flow regime with a negligible pressure difference. To highlight the difference in gas flow capacity under both pressure conditions, a comparison is performed. The discrepancy, the vertical axis value in Figure 5, is calculated by the following equation, quantitatively displaying the difference in apparent permeability.
_ _ _ app l app n app n K K Dis K − = (16) where Kapp_l is the apparent permeability through the pore network with a large pressure difference, nd; Kapp_n is the apparent permeability with a negligible pressure difference, nD; Dis is the discrepancy. To highlight the difference in gas flow capacity under both pressure conditions, a comparison is performed. The discrepancy, the vertical axis value in Figure 5, is calculated by the following equation, quantitatively displaying the difference in apparent permeability.
where K app_l is the apparent permeability through the pore network with a large pressure difference, nd; K app_n is the apparent permeability with a negligible pressure difference, nD; Dis is the discrepancy.
Processes 2022, 10, x FOR PEER REVIEW 10 of 14 Figure 5. Discrepancy caused by different ways to assign pressure condition.
As indicated in Figure 5, there exists a negative relationship between pressure augment and discrepancy; the discrepancy can be as great as 31.1% at 1 MPa, which ranges from 2.8% to 31.1%. It means that the apparent permeability evaluated with the neglected pressure difference will underestimate the essential gas flow capacity through nanoscale porous media. Moreover, it can be observed from Figure 5 that the decline rate is greatly mitigated after pressure exceeds 10 MPa. Furthermore, the discrepancy can be overlooked when the pressure is over 30 MPa, in which case, the slip phenomenon is drastically weakened.

Pore-Network Characteristics
Various approaches have been applied to characterize microscopic pore structures, As indicated in Figure 5, there exists a negative relationship between pressure augment and discrepancy; the discrepancy can be as great as 31.1% at 1 MPa, which ranges from 2.8% to 31.1%. It means that the apparent permeability evaluated with the neglected pressure difference will underestimate the essential gas flow capacity through nanoscale porous media. Moreover, it can be observed from Figure 5 that the decline rate is greatly mitigated after pressure exceeds 10 MPa. Furthermore, the discrepancy can be overlooked when the pressure is over 30 MPa, in which case, the slip phenomenon is drastically weakened.

Pore-Network Characteristics
Various approaches have been applied to characterize microscopic pore structures, particularly for shale gas reservoirs or coalbed methane reservoirs, which is the inherent factor resulting in the nanoconfined gas flow. In this research, pore-network characteristics, including pore size distribution and coordination number, are studied. Specifically, as displayed in Figure 6, three sets of pore size distribution are provided, the case_1 symbol represents porous media dominated by small nanopores, the case_2 symbols represent the medium nanopores, and case_3 symbols represent the majority of large nanopores. from 2.8% to 31.1%. It means that the apparent permeability evaluated with the neglected pressure difference will underestimate the essential gas flow capacity through nanoscale porous media. Moreover, it can be observed from Figure 5 that the decline rate is greatly mitigated after pressure exceeds 10 MPa. Furthermore, the discrepancy can be overlooked when the pressure is over 30 MPa, in which case, the slip phenomenon is drastically weakened.

Pore-Network Characteristics
Various approaches have been applied to characterize microscopic pore structures, particularly for shale gas reservoirs or coalbed methane reservoirs, which is the inherent factor resulting in the nanoconfined gas flow. In this research, pore-network characteristics, including pore size distribution and coordination number, are studied. Specifically, as displayed in Figure 6, three sets of pore size distribution are provided, the case_1 symbol represents porous media dominated by small nanopores, the case_2 symbols represent the medium nanopores, and case_3 symbols represent the majority of large nanopores. Then, the impact of pore size distribution on apparent gas permeability is illustrated in Figure 7. Similarly, the general variation trend remains that the apparent gas permeability will decline with increasing pressure, regardless of the alteration of pore size distribution. In addition, at any pressure, apparent gas permeability with large pore size is always stronger than that with small pore size. The discrepancy caused by different pore size distributions has a negative correlation with increasing pressure. The containment effect of boundary slip mechanisms caused by increasing pressure may take responsibility for the above phenomenon. At low pressure, the slip phenomenon intensifies gas permeability more in small nanopores than in large nanopores. As a result, the difference between case_1 and case_2 increases due to the low-pressure conditions. Then, the impact of pore size distribution on apparent gas permeability is illustrated in Figure 7. Similarly, the general variation trend remains that the apparent gas permeability will decline with increasing pressure, regardless of the alteration of pore size distribution. In addition, at any pressure, apparent gas permeability with large pore size is always stronger than that with small pore size. The discrepancy caused by different pore size distributions has a negative correlation with increasing pressure. The containment effect of boundary slip mechanisms caused by increasing pressure may take responsibility for the above phenomenon. At low pressure, the slip phenomenon intensifies gas permeability more in small nanopores than in large nanopores. As a result, the difference between case_1 and case_2 increases due to the low-pressure conditions. The coordination number denotes the maximum number of pores a single pore can connect to in the pore network, which is generally believed to be in the range of 3 to 7. In this research, three sets of coordination numbers are assigned, and the pore size distribution remains unchanged on the basis of case_2, as presented in Figure 7. It can be observed from Figure 8 that the apparent permeability in the case of a coordination number of 7 can The coordination number denotes the maximum number of pores a single pore can connect to in the pore network, which is generally believed to be in the range of 3 to 7. In this research, three sets of coordination numbers are assigned, and the pore size distribution remains unchanged on the basis of case_2, as presented in Figure 7. It can be observed from Figure 8 that the apparent permeability in the case of a coordination number of 7 can reach as high as 6164 nD. There is an evident tendency that apparent permeability rises dramatically with increasing coordination number. The increasing coordination number undoubtedly enhances the connectivity of the whole pore network; as a result, the possibility for the presence of isolated pores is greatly reduced, which can be attributed as a positive aspect for permeability enhancement. In addition, for a molecule flowing in the network with a large coordination number, there will be more passageways available, and the optimal one would be selected based on the minimum resistance principle.
The coordination number denotes the maximum number o connect to in the pore network, which is generally believed to be this research, three sets of coordination numbers are assigned, an tion remains unchanged on the basis of case_2, as presented in Fig  from Figure 8 that the apparent permeability in the case of a coord reach as high as 6164 nD. There is an evident tendency that app dramatically with increasing coordination number. The increasi undoubtedly enhances the connectivity of the whole pore netwo bility for the presence of isolated pores is greatly reduced, whic positive aspect for permeability enhancement. In addition, for a network with a large coordination number, there will be more pas the optimal one would be selected based on the minimum resista

Conclusions
(1) A self-coded pore network modeling platform for gas flow in nanoscale porous media is established, in which nanoconfined flow mechanisms covering the entire Knudsen number range and pore evolution features, such as pore size distribution and coordination number, are well incorporated. The platform is more flexible on adjusting key parameters compared with the existing tools for PNM. (2) The rarefaction phenomenon acts as a positive factor in advancing gas flow capacity through porous media, particularly in the case of low pressure, where the apparent permeability shall sharply decline from 8103 nD to 598 nD after taking no account of the boundary slip. Moreover, the pressure difference across the inlet and outlet imposes a great influence on gas flow; the apparent permeability would be underestimated by as much as 31.1% with a neglected pressure difference. (3) The variation trend for apparent gas permeability versus pressure remains regardless of the assigned different pore size distribution. Moreover, apparent gas permeability with a large pore size is always stronger than that with a small pore size, discrepancy caused by alteration of pore size distribution has a negative correlation with increasing pressure. Apparent gas permeability rises dramatically with the increasing coordination number, which can be attributed to the reduction of isolated pores as well as the optimized gas flow passageways triggered by the minimum resistance principle.