Numerical Simulation of the Non-Darcy Flow Based on Random Fractal Micronetwork Model for Low Permeability Sandstone Gas Reservoirs

Li, Juhua; Chen, Chen

doi:https://doi.org/10.1155/2020/8884885

Geofluids

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Modeling and Simulation of Micro/Nano-Scale Oil and Gas Migration in Tight Rocks and Shales

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 8884885 | https://doi.org/10.1155/2020/8884885

Numerical Simulation of the Non-Darcy Flow Based on Random Fractal Micronetwork Model for Low Permeability Sandstone Gas Reservoirs

Juhua Li¹and Chen Chen¹

Academic Editor: Jinze Xu

Received17 Jul 2020

Revised24 Aug 2020

Accepted21 Sept 2020

Published14 Oct 2020

Abstract

Darcy’s law is not suit for describing high velocity flow in the near wellbore region of gas reservoirs. The non-Darcy coefficient of the Forchheimer’s equation is a main parameter for the evaluation of seepage capacity in gas reservoirs. The paper presented a new method to calculate by performing gas and con-water flow simulations with random 3D micropore network model. Firstly, a network model is established by random fractal method. Secondly, based on the network simulation method of non-Darcy flow in the literature of Thauvin and Mohanty, a modified model is developed to describe gas non-Darcy flow with irreducible water in the porous medium. The model was verified by our experimental measurements. Then, we investigated the influence of different factors on the non-Darcy coefficient, including micropore structure (pore radius and fractal dimension), irreducible water saturation (), tortuosity, and other reservoir characteristics. The simulation results showed that the value of the non-Darcy coefficient decreases with the increase in all: the average pore radius, fractal dimension, irreducible water saturation, and tortuosity. The non-Darcy coefficients obtained by the fractal method of microparameters are estimated more precisely than the conventional methods. The method provides theoretical support for the productivity prediction of non-Darcy flow in gas reservoirs.

1. Introduction

The Darcy law describes the flow of the subsurface fluid and shows the linear relationship between the pressure gradient and the volume flow (Darcy’s velocity). However, Darcy’s law only applies to laminar and viscous flows. At high velocities, such as high-velocity flows near the wellbore, fluid flow behaves as an inertial flow, and Darcy’s law fails to follow this high-velocity seepage law. Forchheimer proposed a classical Forchheimer equation with a non-Darcy coefficient to estimate non-Darcy effects in Equation (1) and accounted for the impact of both viscous and inertial effects.

Since then, the non-Darcy coefficient has attracted a lot of research. Coefficient can be obtained from analytical methods, physical experiments, and numerical simulations. Geertsma [1] conducted an experiment to establish the quantitative relationship between non-Darcy coefficient and porosity and permeability as well as influence of the immobile fluids on non-Darcy coefficients. Frederick and Graves [2] also studied experimentally the seepage process of gases in porous media considering the influence of both mobile and immobile water on nonlinear seepage of gases. Reid et al. [3] researched high-speed non-Darcy flow process of gas-water fluids in porous media and determined that the presence of irreducible and mobile water affects the gas flow in porous media, both of which reduce the gas permeability and increase the non-Darcy coefficient. However, at the same time, the effect of movable water on gas permeability and non-Darcy coefficient was found to be more significant than that of irreducible water.

Thauvin and Mohanty [4] reported the nonlinear seepage phenomena of gas by establishing the pore structure model considering the influence of pore characteristics on the non-Darcy coefficient. Firoozabadi [5], Jones [6], and others obtained the non-Darcy coefficient for different rocks by physical methods. A 2D random pore network developed by Neeman et al. [7] was found that the porous medium structure has a strong effect on the flow properties. Later, Chukwudozie et al. [8] used LBM to estimate the permeability, tortuosity, and factor of porous rock.

Dou [9] used analytic methods and determined that non-Darcy coefficient is a function of pore size, permeability, friction coefficient, and drag coefficient. Li [10] expanded his analytical research and determined that the is also dependent on the fluid density and fluid viscosity. Wang et al. [11] obtained the relationship between , permeability, and porosity in high-pressure gas reservoirs by statistical analysis of data from multiple sets of rock samples under different overburden tests. Ren et al. [12] established the nonlinear seepage law affected by irreducible water in low permeability reservoirs based on the laboratory seepage experiment results.

Thus, just from the short overview of the previous literature mentioned above, one can see that a lot of theoretical and experimental work was conducted to obtain different forms of non-Darcy coefficient expression, of which the main form is the quantitative relationship between permeability, porosity, and the non-Darcy coefficient. However, due to the different methods applied as well as various porous media and fluids, the non-Darcy coefficient expressions are quite different [13–15]. There are two orders of magnitude difference in the results calculated by different expressions with the same parameters. The fundamental cause of this difference is the intrinsic complexity of micropore structure.

The goal of our work is to obtain an inner relationship of the non-Darcy flow in porous media by a random fractal pore-scale network model. We carried out some flow experiments at different pressure drops and showed the validity of the Forchheimer’s equation. We vary certain pore-scale parameters, study the way it affects various flow parameters, and identify the correlation between different parameters. The mathematical characterization was carried out to confirm the essence of the nonlinear flow. This study is limited to single-phase flow in isotropic porous media.

2. Methodology

The section mainly presents the dynamic simulation method of non-Darcy flow using random fractal network model.

2.1. Construction of Network Model by Random Fractal Theory

Since the first pore network by Fatt [16], a pore numerical model is a versatile studying platform for a variety of subjects [17]. Piri and Blunt [18] developed a quasistatic random network model for researching numerous fluid configurations for two- and three-phase flow. Even since then, the complex geometry of the pore space may be reconstructed [19, 20]. Fang et al. [21] established a fundamental network model, and Mohammadi et al. [22] modified the model to calculate Scri. Consequently, the construction of pore networks is extremely effective for the research of various type of flowing process in porous medium [23–25].

The pore sizes and the throats in the pore network model are also the key to be determined for researchers [26, 27]. It is difficult to describe the irregular pores in traditional geometry. Various studies have shown that the intrinsic properties of the pore meet the statistical significance of fractal distribution [28–31]. We have recently established a constructing digital core method based on fractal theory [32]. The fractal characteristics of pore size are calculated by capillary pressure curve measured from the mercury injection experiment [33]. We then developed a procedure for constructing a 3D pore-throat structure.

Our last paper [27, 32] gives the detailed derivation procedure for the fractal relationship between the pore radius (), fractal dimension (), minimum pore radius (), and random number ().

Porous media are considered to be consisted of large amounts of pore throats which interconnect large and small pores and shrinkage intervals. These pores and pore throats can be described as “points” and “edges” in a 3D network model. The pore nodes are randomly distributed, and the pore throats are connected with a pore as the pore’s coordination number which is generally between 1 and 6. Based on the node position coordinates, the pore’s coordination number, and the principle that the closest pores should be priorly connected, one can first establish the main structural framework of a 3D PNM by connecting the “points” and “edges” and then obtain the geometrical parameters of pores and throats including its radius and length, shape factor, and volume to form a PNM with the topology structure and geometric characteristics of actual core pores.

2.2. Numerical Simulation of Gas Phase Non-Darcy Flow

The dynamic model is adopted using the Matlab program. The pore-throat structure shown in Figure 1 was used as a unit to analyze the pressure drop caused by the fluid flowing into the pores [34]. The fluid flow corresponds to the Hagen-Poiseuille equation during the laminar flow in a porous medium. Its pressure drop is a primary mechanism to overcome the viscous drag. In our model, in addition to considering the viscous drag, the tortuosity and the pressure drop of inertial resistance caused by the gradual pore-throat radius gradient were taken into account.

Based on the pore level law presented by Thauvin and Mohanty [4], the viscous pressure drop between the centers of two adjacent pores can be described by the following equation in single-phase Darcy (slow or Stokes) flow,

In the near wellbore region, the velocity is higher, and Reynolds number increases to the order of Re ~10 in typical gas wells. In porous media, the inertial terms become so important that it affect the additional pressure loss. It was also taken into account:

The pressure drop resulting from the narrowing of the tunnel can be expressed as:

The pressure drop caused by widening of the tunnel is expressed by the following equation:

Based on the above equations, the total pressure drop across the structural unit can be obtained as:

The pressure and velocity relationship in pore level is combined in the network model to calculate macroscopic flow parameters.

Equation (7) is shown below: wherein and are shown, respectively:

Mass balance is applied in one of the principal directions of the pore. The relation between velocity, in a throat, and the pressure drop between the two adjacent body centers is given by Equation (8). To solve this system of equations, we used the Newton-Raphson method. The calculations are repeated at several different pressure drops to verify the Forchheimer equation. The calculations are also repeated for pressure drops applied in each of the pore directions.

In this work, numerical simulations in pore scale are run at specified pressure gradients; meanwhile, the boundary conditions in the network consist of fixed inlet and outlet pressures.

3. Model Validation

We designed six groups of non-Darcy experiments for low permeability sandstone samples. To ensure that the experiments are carried out in porous media, CT scanning was used to scan the samples to avoid the existence of microfractures in the stage of preparation and the end. Finally, two real cores with #1 and #2 were selected.

3.1. Establishment of Pore Network Model

The original data and calculated fractal feature parameters of the rock samples in the certain gas reservoir are given in Table 1 by the capillary pressure curve Figure 2. The first step gets the information of the pore and throat size distribution using Equation (2). Through the relationship between porosity and the pore number, the porosity of rock samples #1 and #2 is 15.2% and 17.2%, respectively, so the corresponding pore numbers are 1053 and 1268. The second step is to generate a pore network model that best fits the original shape of the rock sample, the best model that accurately represents the true capillary pressure characteristics of the core sample.

The size of the random fractal micropore network established in this paper is , as shown in Figure 3, in which the pore numbers are 1053 and 1268; the throat numbers are 2476 and 2982, respectively; the average coordination numbers are 4.362 and 4.614; and the tortuosities are 1.46 and 1.62. In this model, we assume that the shape of pore and throat is triangle, and the value ranges from 0 to 0.0418. The shape factor of pore and throat was given by the Weibull distribution.

(a)

(b)

3.2. Simulation of the Non-Darcy Experiment

Nitrogen and distilled water corresponded to the gas and liquid phases, respectively. The experiments were performed at room temperature. The core was placed in a holder and fixed with a pressure ring. The outlet pressure was fixed at a certain value. Inlet pressure was regulated to change the pressure gradient as needed. Gas flow was monitored and adjusted by flow meters. The experimental results are shown in Table 2.

We proved the validity of the new method by comparing it with the experiments result. The porous network in our model was filled with water and gas as well as the experiments. Water in the pores was in the film form; thus, it did not participate in the flow. The irreducible water in the network model was established according to the principle of small channel occupancy. The gas state was chosen to be constant because of the very small pressure changes under the simulated conditions. The network model had four closed sides and two open ones (the entrance and exit sides). There was no fluid flow through the four closed sides. The flow process was in quasisteady state, and the flow parameters only changed with space but not over time [35]. The model ambient conditions were room temperature. The nitrogen viscosity is 0.011 mPa·s, and the density is 150 kg/m³.

Based on the two sets of network models established earlier in the text, we simulated the single-phase nonlinear gas flow and obtained the average gas flow rate under different pressure gradients. The relationship between gas flow rate and pressure gradient was obtained. The experimental results were compared with the simulated ones (see Figure 4).

(a)

(b)

As can be seen from Figure 4, the flow rate gradually increases with the pressure gradient increase, showing a significant nonlinear relationship. The experimental results are relatively close to the simulated ones.

According to Darcy’s law, the corresponding gas permeability at each flow rate can be obtained using the following equation:

As the gas velocity increases, the gas permeability calculated from the Darcy’s formula decreases mainly because the gas shows nonlinear seepage flow type, and the flow resistance increases together with the flow rate. The calculated permeability was less than the intrinsic absolute permeability of the model.

According to the Forchheimer equation, the seepage velocity and pressure gradient of a fluid can be expressed as a binomial equation, in which the coefficient of the binomial first power phase is the ratio of the fluid viscosity to its permeability . The quadratic phase coefficient is the product of non-Darcy coefficient and fluid density. The quantitative relationship between the pressure gradient and the flow rate can be determined using the above equations. The pressure gradient and the flow rate are obtained from the flow simulation. where is the gas permeability for different gas velocities, is the absolute permeability of the network model, is a non-Darcy coefficient, is the gas density, is the gas viscosity, and is the gas flow rate. The slope of the linear relationship between and in Cartesian coordinates is the non-Darcy coefficient . It was compared with the calculated result of Geertsma formula, which is expressed as:

In the formula, represents the irreducible water saturation, and represents the porosity.

The non-Darcy coefficients calculated by the network model and the experiments are shown in Table 3. However, the results of the two samples calculated by the Geertsma formula are close due to the similar values. The difference between two samples in the internal microscopic structure causes the different coefficients. The non-Darcy coefficients obtained by the fractal method of microparameters are estimated more precisely than the conventional methods. It shows that the non-Darcy flow network model is effective.

4. Factors Affecting the Non-Darcy Coefficient by the Fractal Method

The non-Darcy flow model in random fractal pore network can effectively evaluate the non-Darcy flow coefficient. The internal factors affecting the non-Darcy coefficient are focused on in this section.

4.1. The Impact of Irreducible Water Saturation

To study the influence of irreducible water saturation, we simulated the nonlinear seepage flow under five irreducible water saturation levels (, 0.1, 0.2, 0.35, and 0.5) with different fractal dimensions and obtained the non-Darcy coefficients under different conditions. Several data points from the literature [36] were added to the graph as well (see Figure 5).

Non-Darcy coefficient increased with irreducible water saturation [37–39] was mainly because the irreducible water filled the corners of the pores, only partially leaving the rest of the space for the gas. The indirectly reduced throat radius decreases the flow passage and aggravates the nonlinear flow of gas. Gas flow radius has a quantitative relationship with the gas saturation . Since there are only gas and water phases in the pores, . We used the power function to fit the relationship between non-Darcy coefficient and irreducible water saturation: where is the non-Darcy coefficient in m^-1; and are parameters ( is ~-4.5 according to the fitting).

4.2. The Impact of Fractal Dimension

To analyze the influence of fractal dimension on non-Darcy coefficients, the value of fractal dimension was changed during the reconstruction of the rock sample data by the random fractal method. All other parameters were unchanged. Five groups of three-dimensional network models with , , , , and were established. The simulation of nonlinear gas flow in the presence of irreducible water was carried out, and the non-Darcy coefficients under different fractal dimensions are shown in Figure 6.

The fractal dimension reflects the complexity of pore radius distribution. The larger the fractal dimension, the greater the heterogeneity of pore radii and the larger the non-Darcy coefficients are. The relationship between the non-Darcy coefficient and the fractal dimension can be fitted using exponential function: where is the non-Darcy coefficient in m^-1; and are the parameters (the value of obtained from the fitting was ~4.5).

4.3. The Impact of Average Pore Size

To analyze the effect of pore size on non-Darcy coefficients, the value of was changed during the reconstruction of the rock sample data by using the random fractal method. Five three-dimensional network models with average pore radii of , 3, 5, 7, and 9 μm were established. The simulation of nonlinear seepage of gas under different irreducible water conditions was carried. The corresponding non-Darcy coefficients are shown in Figure 7.

The relationships of non-Darcy coefficient and the average pore radius are consistent under different irreducible water saturations (see Figure 7): the non-Darcy coefficient decreases with the increase of the average pore radius, which agrees with the literature data [36]. This relationship can be fitted using the following power function: where is non-Darcy coefficient in m^-1; and are parameters (the value of obtained from the fitting was ~-4.5).

4.4. The Impact of Tortuosity

To analyze the effect of tortuosity on non-Darcy coefficients, the value of tortuosity was changed, while the other parameters remained unchanged. Five sets of three-dimensional network models with tortuosity of , , , , and were established and obtained non-Darcy coefficients are shown in Figure 8.

The relationship between the non-Darcy coefficient and the tortuosity behaves similarly, and the non-Darcy coefficient increased with the increase of the tortuosity. We fitted this relationship using the following power function: where is non-Darcy coefficient in m^-1; and are parameters (the value of obtained from the fitting was ~1.52).

4.5. Fractal Characterization of Non-Darcy Coefficients

According to the analysis in the previous section, the non-Darcy coefficient can be characterized as a function of the average pore radius, fractal dimension, tortuosity, and irreducible water saturation. The mathematical characterization form is shown below. where is the proportionality constant, and its value depends on the fluid properties (viscosity, density) as well as other parameters. Non-Darcy coefficients and their related parameters were inserted into Equation (17) to obtain for different conditions. The value of was calculated to be equal to through trial and error calculations. Therefore, the mathematical expression of non-Darcy coefficient of this model is:

The non-Darcy coefficient calculated by Equation (18) and the corresponding experimental results are shown in Table 3. We can see that the order of magnitude of non-Darcy coefficient calculated by the experimental data was similar to Equation (18). Equation (18) we conducted by the fractal method of microparameters reflects the influence of rock internal structure on seepage flow. It is more precise than the conventional methods that are obtained by the physical parameters of porosity and permeability. The mathematical characterization was carried out to confirm the essence of the nonlinear flow. This study is limited to single-phase flow in isotropic porous media.

5. Conclusion

(1)Fractal theory was used to analyze the microscopic pore size distribution obtained by two sets of real core mercury injection methods to determine the fractal dimensions and pore sizes. The constructed pore network models are in good agreement with real cores(2)The gas single-phase nonlinear seepage model was established. Besides the viscous resistance, both the tortuosity and the gradual change of the pore-throat radius can produce extraresistance. Based on the three-dimensional network model, the relation curve of flow rate and pressure gradient was obtained by the simulation of the single-phase seepage with irreducible water, and these results were compared with the non-Darcy tests. The validity of the model was confirmed by the good agreement between the obtained non-Darcy coefficient with the calculated results of non-Darcy experiments(3)The sensitivity analysis demonstrated that the value of the non-Darcy coefficient decreases with the increase in all parameters: the average pore radius, fractal dimension, irreducible water saturation, and tortuosity. A quantitative relationship between the non-Darcy coefficients of gas nonlinear seepage and these factors was established. We demonstrated non-Darcy flow characteristics by the fractal method

Nomenclature

:	Fraction of the fluid from a throat that bends at the next body
:	Fraction of the pore throats closed
:	Permeability (Darcy)
:	The gas permeability
:	Permeability tensor (Darcy)
:	Path length from one body center to another (cm)
:	Pore-throat length (cm)
:	Pressure drop (atm)
:	Pressure drop due to bending in the path (atm)
:	Pressure drop due to contraction (atm)
:	Pressure drop due to expansion (atm)
:	Pressure drop for Stokes’ flow (atm)
:	Pressure (atm)
:	Flow rate (cm³/s)
:	Body radius (cm), is the radius of pore 1, is the radius of pore 2
:	Throat radius (cm)
:	Minimum value of the radius (cm)
:	Maximum value of the radius (cm)
:	The pore radius (cm)
:	Fractal dimension
:	Random number
:	Interstitial velocity (cm/s)
:	Superficial velocity (cm/s)
:	Non-Darcy coefficient (cm^-1)
:	Viscosity (g·cm^-1·s^-1)
:	Porosity (%)
:	Density (g·cm^-3)
:	Tortuousity.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is financially supported by the National Natural Science Foundation of China (No. 51504039) and Major National Science and Technology Projects of China (No. 2016ZX05060-019).

References

J. Geertsma, “Estimating the coefficient of inertial resistance in fluid flow through porous media,” Society of Petroleum Engineers Journal, vol. 14, no. 5, pp. 445–450, 2013.
View at: Publisher Site | Google Scholar
D. C. Frederick and R. M. Graves, “New Correlations To Predict Non-Darcy Flow Coefficients at Immobile and Mobile Water Saturation,” in SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 1994.
View at: Publisher Site | Google Scholar
R. B. Grigg and M. K. Hwang, “High velocity gas flow effects in porous gas-water system,” in SPE Gas Technology Symposium, Calgary, Alberta, Canada, 1998.
View at: Publisher Site | Google Scholar
F. Thauvin and K. K. Mohanty, “Network modeling of non-Darcy flow through porous media,” Transport in Porous Media, vol. 31, no. 1, pp. 19–37, 1998.
View at: Publisher Site | Google Scholar
A. Firoozabadi and D. L. Katz, “An analysis of high-velocity gas flow through porous media,” Journal of Petroleum Technology, vol. 31, no. 2, pp. 211–216, 2013.
View at: Publisher Site | Google Scholar
F. O. Jones and W. W. Owens, “A laboratory study of low-permeability gas sands,” Journal of Petroleum Technology, vol. 32, no. 9, pp. 1631–1640, 2013.
View at: Publisher Site | Google Scholar
H.-W. Lao, H. J. Neeman, and D. V. Papavassiliou, “A pore network model for the calculation of non-Darcy flow coefficients in fluid flow through porous media,” Chemical Engineering Communications, vol. 191, no. 10, pp. 1285–1322, 2004.
View at: Publisher Site | Google Scholar
C. P. Chukwudozie, M. Tyagi, S. O. Sears, and C. D. White, “Prediction of non-Darcy coefficients for inertial flows through the Castlegate sandstone using image-based modeling,” Transport in Porous Media, vol. 95, no. 3, pp. 563–580, 2012.
View at: Publisher Site | Google Scholar
H. Dou, “Further understanding on turbulence coefficient of high-speed non-Darcy flow of gas reservoir,” Fault Blcok Oil and Gas Field, vol. 20, no. 4, pp. 465–469, 2013.
View at: Google Scholar
C. Li, “High rate non-Darcy flow is easier to occur in low permeability reservoirs,” Lithologic Reservoirs, vol. 11, no. 6, pp. 111–119, 2011.
View at: Google Scholar
X. Wang, D. Zhang, and M. Liao, “Turbulence factor of high pressure gas reservoir in Keyi tectonic belt,” Drilling & Production Technology, vol. 25, no. 1, pp. 44–49, 2002.
View at: Google Scholar
X. Ren, Q. Yan, and Q. He, “Experimental study of low permeability gas reservoir gas seepage characteristics,” Journal of Xi' an Shiyou University (Natural Science Edition), vol. 12, no. 3, pp. 22–25, 1997.
View at: Google Scholar
D. Li and T. W. Engler, “Literature review on correlations of the non-Darcy coefficient,” in SPE Permian Basin Oil and Gas Recovery Conference, Midland, Texas, 2001.
View at: Publisher Site | Google Scholar
K. Xu, X. Lei, and Q. Meng, “Study of inertial coefficient of non-Darcy seepage flow,” Chinese Journal of Rock Mechanics and Engineering, vol. 31, no. 1, pp. 164–170, 2012.
View at: Google Scholar
Y. Chen and N. Dong, “A method of determining high-velocity coefficient β and its correlation,” Fault Blcok Oil and Gas Field, vol. 5, no. 6, pp. 20–26, 1998.
View at: Google Scholar
I. Fatt, “The network model of porous media I. capillary pressure characteristics,” Transactions of AIME, vol. 207, no. 1, pp. 144–181, 2013.
View at: Publisher Site | Google Scholar
H. Okabe and M. J. Blunt, “Pore space reconstruction using multiple-point statistics,” Journal of Petroleum Science and Engineering, vol. 46, no. 1-2, pp. 121–137, 2005.
View at: Publisher Site | Google Scholar
M. Piri and M. J. Blunt, “Three-dimensional mixed-wet random pore-scale network modeling of two- and three-phase flow in porous media. I. Model description,” Physical Review E, vol. 71, no. 2, article 026301, 2005.
View at: Publisher Site | Google Scholar
P. A. Byrnes, R. M. Cluff, and J. C. Webb, “Analysis of critical permeability, capillary and electrical properties for Mesaverde tight gas sandstones from Western US basins: final scientific,” Tech. Rep., Technical report submitted to DOE and NETL, 2009.
View at: Google Scholar
N. A. Idowu, “Pore-scale modeling: stochastic network generation and modeling of rate effects in water flooding,” Tech. Rep., PhD thesis, Imperial College London, UK, 2009.
View at: Google Scholar
F. Fang, A. Firoozabadi, M. Abbaszadeh, and C. Radke, “A phenomenological modeling of critical condensate saturation,” in SPE Annual Technical Conference and Exhibition, Denver, Colorado, 1996.
View at: Publisher Site | Google Scholar
S. Mohammadi, K. S. Sorbie, and A. Danesh, “Pore-level modelling of gas-condensate flow through horizontal porous media,” in SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 2009.
View at: Publisher Site | Google Scholar
X. Zhao, Z. Yang, W. Lin et al., “Study on pore structures of tight sandstone reservoirs based on nitrogen adsorption, high-pressure mercury intrusion, and rate-controlled mercury intrusion,” Journal of Energy Resources Technology, vol. 141, no. 11, article 112903, 2019.
View at: Publisher Site | Google Scholar
S. Zhang, X. Xian, J. Zhou et al., “Experimental study of the pore structure characterization in shale with different particle size,” Journal of Energy Resources Technology, vol. 140, no. 5, article 054502, 2018.
View at: Publisher Site | Google Scholar
C. Lu, W. Zhao, Y. Liu, and X. Dong, “Pore-scale transport mechanisms and macroscopic displacement effects of in-situ oil-in-water emulsions in porous media,” Journal of Energy Resources Technology, vol. 140, no. 10, article 102904, 2018.
View at: Publisher Site | Google Scholar
F. Alreshedan and A. Kantzas, “Investigation of permeability, formation factor, and porosity relationships for Mesaverde tight gas sandstones using random network models,” Journal of Petroleum Exploration and Production Technology, vol. 6, no. 3, pp. 545–554, 2016.
View at: Publisher Site | Google Scholar
Y. Yi, J. Li, and L. Ji, “Numerical Determinationof critical condensate saturation in gas condensate reservoirs,” Journal of Energy Resources Technology, vol. 139, no. 6, article 062801, 2017.
View at: Publisher Site | Google Scholar
B. Yu and P. Cheng, “A fractal permeability model for bi-dispersed porous media,” International Journal of Heat and Mass Transfer, vol. 45, no. 14, pp. 2983–2993, 2002.
View at: Publisher Site | Google Scholar
J. Wu and B. Yu, “A fractal resistance model for flow through porous media,” Heat and Mass Transfer, vol. 50, no. 19-20, pp. 3925–3932, 2007.
View at: Publisher Site | Google Scholar
F. Wang, Z. Liu, L. Jiao, C. Wang, and H. Guo, “A fractal permeability model coupling boundary-layer effect for tight oil reservoirs,” Fractals, vol. 25, no. 5, article 1750042, 2017.
View at: Publisher Site | Google Scholar
F. Wang, Z. Liu, J. Cai, and J. Gao, “A fractal model for low-velocity non-Darcy flow in tight oil reservoirs considering boundary-layer effect,” Fractals, vol. 26, no. 5, article 1850077, 2018.
View at: Publisher Site | Google Scholar
J. Li and B. Zheng, “Digital core and pore network model reconstruction based on random fractal theory,” International Journal of Energy and Statistics, vol. 3, no. 1, article 1550001, 2015.
View at: Publisher Site | Google Scholar
F. Wang, L. Jiao, Z. Liu, X. Tan, C. Wang, and J. Gao, “Fractal analysis of pore structures in low permeability sandstones using mercury intrusion porosimetry,” Journal of Porous Media, vol. 21, no. 11, pp. 1097–1119, 2018.
View at: Publisher Site | Google Scholar
M. T. Balhoff and M. F. Wheeler, “A predictive pore-scale model for non-Darcy flow in porous media,” SPE Journal, vol. 14, no. 4, pp. 579–587, 2009.
View at: Publisher Site | Google Scholar
G. Wang, X. Yang, and X. Zhang, “Numerical simulation on non-Darcy seepage of CBM by means of 3D reconstruction based on computed tomography,” Journal of China Coal Society, vol. 41, no. 4, pp. 931–940, 2016.
View at: Google Scholar
J. M. Lombard, D. Longeron, and F. Kalaydjian, “Influence of connate water and condensate saturation on inertial effects in gas-condensate fields,” in SPE Annual Technical Conference and Exhibition, Houston, Texas, 1999.
View at: Publisher Site | Google Scholar
W. Li, X. Ma, and Y. Wang, “Effect of irreducible water saturation on deliverability of hydraulic fractured gas well taking into account non-Darcy flow,” Oil Drilling & Production Technology, vol. 33, no. 3, pp. 35–37, 2011.
View at: Google Scholar
N. Li, X. Tang, and Q. Zhang, “Experimental study on gas low-velocity non-Darcy seepage in low permeability gas reservoir,” Natural Gas Exploration & Development, vol. 26, no. 2, pp. 49–55, 2003.
View at: Google Scholar
R. D. Evans, C. S. Hudson, and J. E. Greenlee, “The effect of an immobile liquid saturation on the non-Darcy flow coefficient in porous media,” SPE Production Engineering, vol. 2, no. 4, pp. 331–338, 1987.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Juhua Li and Chen Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies