Next Article in Journal
Thermal Transmittance Measurements of the Historical Masonries: Some Case Studies
Next Article in Special Issue
A Numerical Study on the Diversion Mechanisms of Fracture Networks in Tight Reservoirs with Frictional Natural Fractures
Previous Article in Journal
Decomposition and Forecasting of CO2 Emissions in China’s Power Sector Based on STIRPAT Model with Selected PLS Model and a Novel Hybrid PLS-Grey-Markov Model
Previous Article in Special Issue
Evolution of Coal Permeability during Gas Injection—From Initial to Ultimate Equilibrium
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Non-Linear Flow Model for Porous Media Based on Conformable Derivative Approach

1
College Petroleum of Engineering & Geoscience, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia
2
College of Petroleum Engineering, China University of Petroleum, Beijing 102249, China
3
School of Mechanics and Civil Engineering, China University of Mining and Technology, Beijing 100083, China
4
College of Engineering, Peking University, Beijing 100871, China
*
Author to whom correspondence should be addressed.
Submission received: 11 September 2018 / Revised: 24 October 2018 / Accepted: 29 October 2018 / Published: 1 November 2018
(This article belongs to the Special Issue Flow and Transport Properties of Unconventional Reservoirs 2018)

Abstract

:
Prediction of the non-linear flow in porous media is still a major scientific and engineering challenge, despite major technological advances in both theoretical and computational thermodynamics in the past two decades. Specifically, essential controls on non-linear flow in porous media are not yet definitive. The principal aim of this paper is to develop a meaningful and reasonable quantitative model that manifests the most important fundamental controls on low velocity non-linear flow. By coupling a new derivative with fractional order, referred to conformable derivative, Swartzendruber equation and modified Hertzian contact theory as well as fractal geometry theory, a flow velocity model for porous media is proposed to improve the modeling of Non-linear flow in porous media. Predictions using the proposed model agree well with available experimental data. Salient results presented here include (1) the flow velocity decreases as effective stress increases; (2) rock types of “softer” mechanical properties may exhibit lower flow velocity; (3) flow velocity increases with the rougher pore surfaces and rock elastic modulus. In general, the proposed model illustrates mechanisms that affect non-linear flow behavior in porous media.

1. Introduction

Ever since Henry Darcy (1865) developed his famous linear flow model (the classical Darcy’s law), based on a series of sand pack experiments, the linear flow through porous media has drawn tremendous attention in various scientific and engineering field [1,2]. However, it’s a common phenomenon that experiments on low velocity flow in tight porous media, deviate from the Darcy’s law and the flow velocity is lower than that predicted from Darcy’s law. As stated in the literature, the existence of low velocity non-Darcy flow (or low velocity non-linear flow) in tight porous media (e.g., shale gas/oil reservoirs, coalbed, or tight gas/oil reservoirs) is due to the interaction forces between the fluid and tight pores [3,4]. Many scholars have documented that, there existed threshold Reynolds number or pressure gradient, which could be used to well describe low velocity non-linear flow [5,6,7]. And they concluded that there is no flow in tight porous media when the pressure gradient is beyond the certain value (i.e., threshold pressure gradient). However, Li provided contradictory evidence for the threshold pressure gradient [8]. He suggested that the threshold pressure gradient measured in labs can be probably ascribed to the difficulty in measuring lower flow velocity, and the false phenomenon of the existence of threshold pressure gradient is strengthened by the skin effect. Until now, appropriate non-linear model for fluid flow through tight porous media remains unclear, though some more formulas have been established to describe the non-linear flow, such as power function model [9], exponential function model [10], incomplete Gamma function model [11,12] and fractional derivative approach [13,14,15]. As stated in the literature, in some extent, these models above are suitable for the description of non-linear flow in porous media with lower permeability [9,10,11,12,13,14]. Most recently, Yang [14] and Zhou [15] also suggested that the conformable derivative approach is suitable for the describing non-linear flow in low-permeability porous media. However, these models above never took effective stress into account. It is reported that the porous media will be compressed as the effective stress increases, causing fluid flow behavior in porous media to be strongly stress-dependent [16,17,18,19,20,21,22,23,24].
As implied by this brief literature review, we suggest that the characteristic behavior of non-linear flow is still not definitively determined. Therefore, a major goal of this research was to develop an analytical model in a closed form for the description of low velocity non-linear flow. The specific objectives of this work were: (1) to establish a reasonable quantitative model to quantify the essential controls on non-linear flow; (2) to verify the model with available experimental data. Compared with the previous models, our model takes into account more factors, including the influence of the effective stress and the microstructural parameters of the pore space. The proposed models can reveal more mechanisms that affect the low velocity non-linear flow in porous media.

2. Mathematical Model

In this section, the analytical low velocity non-linear flow model for porous media is detailed. The conformable derivative is used to develop the Swartzendruber model for description of low velocity non-linear flow in pores, the fractal geometry theory and modified Hertzian contact theory are used to describe the complex pore structure of porous media under stress condition.

2.1. Model Assumptions

The following assumptions are made to simplify the flow system:
  • The porous media is composed by a bundle of capillary bundles and a single capillary with the equivalent radius r is made up of a packing of equivalent spherical grains.
  • The interspaces in porous media have fractal characteristics.
  • The single phase flow is under isothermal and stress condition, which is fully developed and at steady state.
  • The deformation of porous media obeys Hertzian contact theory.
  • During the flow, the fluid has constant viscosity and density.

2.2. Conformable Derivative Approach to Swartzendruber Equation

As suggested by decades of literature, the Swartzendruber equation can well describe the non-Darcian flow in tight porous media with low permeability [13,14]. Based on the Swartzendruber equation, the following equation can be written as:
d u d i = K ρ g μ ( 1 e i I ) ,
where u is the flow velocity in the cross section; K is the permeability of porous media; μ is fluid viscosity; ρ is fluid density, g is the gravitational acceleration, i and I represent hydraulic gradient and threshold hydraulic gradient, respectively. According to Equation (1), the flow velocity in a single capillary with radius r can be written as:
d u d i = r 2 ρ g 8 μ ( 1 e i I ) ,
where r is the radius of the capillary.
As the linear operator does not inherit all the operational behaviors from the typical first derivative, the Swartzendruber equation fails to capture the full range of non-Darcy flow behavior in porous media [13,14,15]. Fortunately, as a well-behaved and efficient method, conformable derivative approach with real order can be applied to address this problem. The conformable derivative of the flow velocity u(i): [0, ∞) → R for all i > 0 with order α ∈ (0, 1] can be defined by [10]:
T α u ( i ) = lim ε 0 u ( i + ε i 1 α ) u ( i ) ε ,
and the conformable derivative at 0 is given by ( T α u ) ( 0 ) = lim i 0 ( T α u ) ( i ) .
As stated in the literature [10,14,25], the relationship between the conformable derivative and the first derivative can be written as:
T α u ( i ) = i 1 α d u ( i ) d i .
Equation (3b) shows that the conformable derivative coincides with the classical first derivative with a given differential order α = 1, which means the conformable derivative is a modification of classical derivative in direction and magnitude [25,26,27]. As stated in the literature [25,26,27], the physical interpretation of the conformable derivative is a modification of classical derivative indirection and magnitude. Replacing the first order derivative in Equation (2) with conformable derivative, the Swartzendruber equation can be rewritten as:
T α u ( i ) = r 2 ρ g 8 μ ( 1 e i I ) .
By solving Equation (4) with Laplace transform and inverse Laplace transform, the flow velocity in the cross section can be determined as:
u ( i ) = r 2 ρ g 8 μ i α α [ 1 F 1 1 ( α ; α + 1 ; i I ) ] ,
where F 1 1 ( α ; α + 1 ; i I ) = j = 0 α ( j ) ( α + 1 ) ( j ) ( i / I ) j j ! is the Kummer confluent hypergeometric function [23].
Based on Equation (5), the flow rate q in the cross section can be written as:
q = π r 4 ρ g 8 μ i α α [ 1 F 1 1 ( α ; α + 1 ; i I ) ] ,
where q is the flow rate in the capillary with the radius r.

2.3. Non-Linear Flow Model

The fractal theory is used to develop the non-linear flow model for tight porous media. Based on fractal theory and modified Hertzian contact theory described in detail in [28,29], the total volumetric flow rate Q in the cross section under stress condition can be calculated by integrating the flow rate over the radius ranging from the minimum radius rmin to the maximum porous radius rmax [28,29,30,31,32,33]:
Q = N r min r max q f d r ,
with:
{ N = ( r max / r min ) D f ; f ( r ) = D f r min D f r ( D f + 1 ) , r = r 0 { 1 4 [ 3 π ( 1 v 2 ) p eff 4 E ] β } , r min = r min 0 { 1 4 [ 3 π ( 1 v 2 ) p eff 4 E ] β } , r max = r max 0 { 1 4 [ 3 π ( 1 v 2 ) p eff 4 E ] β } .
In Equation (8), N is the number of pores; r0 and r are the initial pore radius and pore radius under stress condition, respectively. f is the probability density function for pore size distribution; β is the power law index which is related to the structure of pore surface. rmax0 and rmin0 are the maximum and minimum pore radius of porous media at zero stress, respectively. rmax and rmin are the stress-dependent maximum and the stress-dependent minimum pore radius of porous media, respectively. peff is the effective stress, E is rock elastic modulus and v is rock Poisson’s ratio. Df is the fractal dimension for pore size distribution which can be determined as [32,33]:
{ D f = 2 ( 2 D f 0 ) r max 0 ( 3 D f 0 ) r max ( 2 D f 0 ) r max 0 , D f 0 = 2 ln φ 0 ln ( r min 0 / r max 0 ) ,
where φ0 is the initial porosity of the porous media, and Df0 is the fractal dimension for pore size distribution at zero stress.
Substituting Equations (6), (8) and (9) into Equation (7), the flow rate can be rewritten as:
Q = π ρ g D f r max D f 8 μ i α α [ 1 F 1 1 ( α ; α + 1 ; i I ) ] r max 4 D f r min 4 D f 4 D f .
As stated in the literature [32,33,34,35], the cross sectional area of a unit cell A can be expressed as:
A = π D f r max 2 [ 1 ( r min / r max ) 2 D f ] ( 2 D f ) ( r min / r max ) 2 D f .
Then, based on Equations (10) and (11), the average flow velocity uav can be written as:
u av = Q A = ρ g φ i α ( 2 D f ) ( r max 2 φ r min 2 ) 8 α μ ( 1 φ ) ( 4 D f ) [ 1 F 1 1 ( α ; α + 1 ; i I ) ]
It is evident that the flow rate (or average flow velocity) is the function of pore structural parameters, power law index, rock elastic modulus, effective stress, and differential order α as well as hydraulic gradient and threshold hydraulic gradient.

3. Results and Discussion

This section aims at studying the novel analytical models in detail. In the following, we first compare our results with those from experimental data. Then, in order to analyze essential controls on non-linear flow in tight porous media, the effects of relevant parameters on average flow velocity are studied in detail.
The availability of the proposed model Equation (12) depends on its ability to adequately fit experimental data [15]. To verify our quantitative model, the measured average flow velocity versus hydraulic gradient relationship in [36] and that predicted by our proposed model are compared (Figure 1). In the experiment of Prakash K. et al. [36], the flow velocity tests were conducted on soils at an effective consolidation stress of 6.25 kPa during loading process. In our proposed model, the initial porosity of rock is 15%, the rock elastic modulus is 45 GPa, the rock Poisson’s ratio is 0.23 and power law index is 3/4. Furthermore, to ensure the effective consolidation stress is 6.25 kPa, the effective stress assigned is 6.25 kPa. The values of other parameters (e.g., rmax0, rmin0, differential order α and threshold hydraulic gradient I) are listed in the Figure 1. Results displayed in Figure 1 suggest that the proposed model calculated relationship between average flow velocity and hydraulic gradient is in good agreement with that determined by experimental data [36]. Results (Figure 1) also suggest a definitive positive correlation between the average flow velocity and hydraulic gradient.
Figure 2 provides a comparison of the relationship between average flow velocity and hydraulic gradient predicted by the proposed model with experimental data of [37]. In the experiment of Zhang et al. [37], the permeability tests were conducted on gap-graded sands (e.g., sand A, sand B and sand C) to determine the critical hydraulic gradient of piping in sands. Sand A, sand B and sand C were prepared by mixing a coarse sand with particle size of 3–5 mm and a fine sand with particle size of 0.25–0.50 mm at various ratios of 4:1, 5:1 and 6:1 respectively.
For sand A, the hydraulic conductivity and void ratio are 2.6 × 10−3 m/s and 0.56 respectively. The hydraulic conductivity and void ratio of sand B are 3.6 × 10−3 m/s and 0.60 respectively. In addition, for sand C, the hydraulic conductivity and void ratio are 4.0 × 10−3 m/s and 0.63 respectively. To ensure the porous media simulated in the proposed model exhibits the same physical properties as the porous media tested in the experiments of [37], the initial void ratio e0 and initial porosity φ0 (i.e., φ0 = e0/(1 + e0)) applied in the proposed model are the same with that of the sands in the experiments and the effective stress assigned is 0 MPa. The values of other parameters (e.g., rmax0, rmin0, differential order α and threshold hydraulic gradient I) are listed in the Figure 2. It can be seen from Figure 2 that our predicted values agree well with the corresponding experimental data [37].
As the average flow velocity is related to the threshold hydraulic gradient, differential order α, effective stress, rock elastic modulus, power law index, pore structural parameters, and the hydraulic gradient. We will then study the effects of these relevant parameters (e.g., threshold hydraulic gradient, differential order α, effective stress, rock elastic modulus, power law index and initial porosity) on average flow velocity in detail to analyze essential controls on non-linear flow in tight porous media. We plotted average flow velocity vs. hydraulic gradient with different threshold hydraulic gradient I and different differential order α (Figure 3a,b).
The parameters assigned in the proposed model were listed in the Figure 3. As suggested by the results (Figure 3a), average flow velocity decreases as threshold hydraulic gradient increases. In addition, Figure 3b shows that larger different differential order leads to smaller average flow velocity with smaller hydraulic gradient, however, on the contrary, when the hydraulic gradient increases to a certain value, average flow velocity increases as differential order increases.
Figure 4 illustrates the average flow velocity vs. hydraulic gradient with different effective stress. The parameters assigned in the proposed model were listed in the Figure 4. As suggested by the results (Figure 4), average flow velocity decreases as effective stress increases. This may be attributed to the decrease of pore radius which is resulted from the pore compaction. Therefore, fluid flow behavior in tight porous media is strongly stress-dependent.
Plotting average flow velocity vs. hydraulic gradient with different rock elastic modulus, power law index and initial porosity (Figure 5) was also useful. For the necessary calculations, the parameters assigned in the proposed model were listed in the Figure 5. As suggested by the results (Figure 5a), the average flow velocity increases as rock elastic modulus increases. We interpret this result to indicate that the larger rock elastic modulus decreases the contact surface radius of a given particle, leading to reduced pore volume compressibility and larger hydraulic conductivity. Therefore, rock types of “soft” lithology can yield lower flow velocity. As suggested by the results (Figure 5b), the average flow velocity increases as power law index increases. Correspondingly, we suggest that the larger power law index β implies rougher pore surfaces, leading to only a limited number of pores being compressed and the larger hydraulic conductivity. Figure 5c shows the average flow velocity decreases as rock initial porosity decreases, which is expected. Correspondingly, we suggest that the smaller initial porosity implies narrower pore radius, leading to the smaller hydraulic conductivity.

4. Conclusions

In this study, we developed a novel non-linear flow model for tight porous media. The model is based on Swartzendruber equation and conformable derivative approach and as well as the modified Hertzian contact theory and fractal geometry, and allowed us to analyze essential controls on non-linear flow in tight porous media. An advantage of this model is that it lacks empirical constants, and, more importantly, every parameter in the model has specific physical significance. Predictions from the proposed analytical model exhibit similar variation trends as experimental data, suggesting validity of the model to predict the average flow velocity. Moreover, we analyzed resulting model predictions in detail, to confirm the model’s robustness in this context. Results of the new model show the following salient conclusions:
  • The proposed models indicate that average flow velocity in tight porous media is a function of microstructural parameters of the pore space, rock lithology and differential order α as well as hydraulic gradient and threshold hydraulic gradient.
  • The parametric study reveals that average flow velocity increases with the rougher pore surfaces and rock elastic modulus, and decreases with increasing effective stress. “Softer” rock lithology may yield lower average flow velocity.
  • This non-linear model presented here considers microstructural parameters of pore space and rock lithology; we have shown that its forecasted values are robust, at least compared to experimental data, and thus may be useful for performance predictions of non-linear flow behavior in tight porous media. Results also reveal more information about the details of specific parameters (and therefore mechanisms) that affect non-linear flow behavior in porous media. The new model presented in this work can be used to depict the non-linear flow in tight porous media, and may provide meaningful applications for design and development of tight reservoirs. In addition, as the model takes effective stress into account, it is also useful for performance predictions of the coupled flow deformation behavior (stress sensitivity) in tight porous media.

Author Contributions

Conceptualization, G.L.; Funding acquisition, H.W.; Investigation, N.C. and D.L.; Methodology, G.L.; Project administration, H.W.; Visualization, N.C.; Writing—original draft, G.L.; Writing—review & editing, H.W.

Funding

This research was funded by the State Major Science and Technology Special Project of China during the 13th Five-Year Plan (grant number: 2016ZX05025-003-007 and 2016ZX05014004-006). Appreciation and thanks are given to the organizations for their finical funding herein.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

Latin symbols
ACross sectional area of a unit cell, µm2
Df0Initial pore area fractal dimension at zero stress, dimensionless
DfPore area fractal dimension, dimensionless
e0Initial void ratio of porous media, dimensionless
ERock elastic modulus of porous media, GPa
fProbability density function for pore size distribution, dimensionless
1F1Kummer confluent hypergeometric function
gGravitational acceleration, N/kg
KAbsolute permeability of porous media, µm2
iHydraulic gradient, dimensionless
IThreshold hydraulic gradient, dimensionless
NNumber of pores of a unit cell, dimensionless
peffEffective stress, MPa
qFlow rate in the cross section, m3/s
QTotal volumetric flow rate in the cross section under stress condition, m3/s
r0Initial equivalent pore radius of capillary at zero stress, µm
rEquivalent pore radius of capillary under effective stress, µm
uFlow velocity in the cross section, m/s
uavAverage flow velocity, m/s
Greek symbols
αDifferential order
βPower law index, dimensionless
μFluid viscosity, mPa·s
ρFluid density, kg/m3
φ0Initial porosity of porous media, dimensionless
φPorosity under effective stress, dimensionless
νPoisson’s ratio, dimensionless
Subscript
avAverage
effEffective
maxmaximum values
max0Initial maximum values at zero stress
minminimum values
min0Initial minimum values at zero stress

References

  1. Freeze, R.A.; Cherry, J.A. Groundwater; Prentice-Hall: Upper Saddle River, NJ, USA, 1979. [Google Scholar]
  2. Zhang, J.; Xu, Q.; Chen, Z. Seepage analysis based on the unified unsaturated soil theory. Mech. Res. Commun. 2001, 28, 107–112. [Google Scholar] [CrossRef]
  3. Xu, C.; Kang, Y.; Chen, F.; You, Z. Analytical model of plugging zone strength for drill-in fluid loss control and formation damage prevention in fractured tight reservoir. J. Pet. Sci. Eng. 2017, 149, 686–700. [Google Scholar] [CrossRef]
  4. Cai, J.; Lin, D.; Singh, H.; Wei, W.; Zhou, S. Shale gas transport model in 3D fractal porous media with variable pore sizes. Mar. Geol. 2018, 98, 437–447. [Google Scholar] [CrossRef]
  5. Thomas, L.K.; Katz, D.L.; Tek, M.R. Threshold pressure phenomena in porous media. Soc. Petrol. Eng. J. 1968, 8, 174–184. [Google Scholar] [CrossRef]
  6. Tian, W.; Li, A.; Ren, X.; Josephine, Y. The threshold pressure gradient effect in the tight sandstone gas reservoirs with high water saturation. Fuel 2018, 226, 221–229. [Google Scholar] [CrossRef]
  7. Morozov, P.; Abdullin, A.; Khairullin, M. An analytical model of SAGD process considering the effect of threshold pressure gradient. IOP Conf. Ser. Earth Environ. Sci. 2018, 155, 012001. [Google Scholar] [CrossRef]
  8. Li, C. Is a starting pressure gradient necessary for flow in porous media? Acta Petrol. Sin. 2010, 31, 867–870. [Google Scholar]
  9. Hansbo, S. Consolidation equation valid for both Darcian and non-Darcian flow. Geotech 2001, 51, 51–54. [Google Scholar] [CrossRef]
  10. Swartzendruber, D. Modification of Darcy’s law for the flow of water in soils. Soil Sci. 1962, 93, 22–29. [Google Scholar] [CrossRef]
  11. Liu, H.H.; Birkholzer, J. On the relationship between water flux and hydraulic gradient for unsaturated and saturated clay. J. Hydrol. 2012, 475, 242–247. [Google Scholar] [CrossRef] [Green Version]
  12. Liu, H.H. Non-Darcian flow in low-permeability media: Key issues related to geological disposal of high-level nuclear waste in shale formations. Hydrogeol. J. 2014, 22, 1525–1534. [Google Scholar] [CrossRef]
  13. Wang, R.; Zhou, H.W.; Zhong, J.C.; Yi, H.; Chen, C.; Zhao, Y. The study on non-Darcy seepage equation of low velocity flow. Sci. Sin. Phys. 2017, 47, 064702. (In Chinese) [Google Scholar] [CrossRef]
  14. Yang, S.; Wang, L.; Zhang, S. Conformable derivative: Application to non-Darcian flow in low-permeability porous media. Appl. Math. Lett. 2018, 79, 105–110. [Google Scholar] [CrossRef]
  15. Zhou, H.W.; Yang, S.; Zhang, S.Q. Conformable derivative approach to anomalous diffusion. Physica A 2018, 491, 1001–1013. [Google Scholar] [CrossRef]
  16. Zhang, H.J.; Jeng, D.S.; Barry, D.A.; Seymour, B.R.; Li, L. Solute transport in nearly saturated porous media under landfill clay liners: A finite deformation approach. J. Hydrol. 2013, 479, 189–199. [Google Scholar] [CrossRef] [Green Version]
  17. Srinivasacharya, D.; Srinivasacharyulu, N.; Odelu, O. Flow and heat transfer of couple stress fluid in a porous channel with expanding and contracting walls. Int. Commun. Heat Mass Transf. 2009, 36, 180–185. [Google Scholar] [CrossRef]
  18. Neto, L.B.; Kotousov, A.; Bedrikovetsky, P. Elastic properties of porous media in the vicinity of the percolation limit. J. Pet. Sci. Eng. 2011, 78, 328–333. [Google Scholar] [CrossRef]
  19. Mokni, N.; Olivella, S.; Li, X.; Smets, S.; Valcke, E. Deformation induced by dissolution of salts in porous media. Phys. Chem. Earth Parts A/B/C 2008, 33, S436–S443. [Google Scholar] [CrossRef]
  20. Gangi, A.F. Variation of whole and fractured porous rock permeability with confining pressure. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1978, 15, 249–257. [Google Scholar] [CrossRef]
  21. Archer, R.A. Impact of stress sensitive permeability on production data analysis. In Proceedings of the SPE Unconventional Reservoirs Conference, Keystone, CO, USA, 10–12 January 2008. SPE-114166-MS. [Google Scholar]
  22. Schreyer-Bennethum, L. Theory of flow and deformation of swelling porous materials at the macroscale. Comput. Geotech. 2007, 34, 267–278. [Google Scholar] [CrossRef]
  23. Jennings, J.B.; Carroll, H.B.; Raible, C.J. The relationship of permeability to confining pressure in low permeability rock. In Proceedings of the SPE/DOE Low Permeability Gas Reservoirs Symposium, Denver, CO, USA, 27–29 May 1981; p. SPE-9870-MS. [Google Scholar]
  24. Tan, X.H.; Li, X.P.; Liu, J.Y.; Zhang, L.H.; Fan, Z. Study of the effects of stress sensitivity on the permeability and porosity of fractal porous media. Phys. Lett. A 2015, 379, 2458–2465. [Google Scholar] [CrossRef]
  25. Khalil, R.; Al Horani, M.; Yousef, A.; Sababheh, M. A new definition of fractional derivative. J. Comput. Appl. Math. 2014, 264, 65–70. [Google Scholar] [CrossRef]
  26. Abdeljawad, T. On conformable fractional calculus. J. Comput. Appl. Math. 2015, 279, 57–66. [Google Scholar] [CrossRef]
  27. Zhao, D.; Luo, M. General conformable fractional derivative and its physical interpretation. Calcolo 2017, 54, 903–917. [Google Scholar] [CrossRef]
  28. Lei, G.; Dong, Z.; Li, W.; Wen, Q.; Wang, C. Theoretical study on stress sensitivity of fractal porous media with irreducible water. Fractals 2018, 26, 1850004. [Google Scholar] [CrossRef]
  29. Lei, G.; Mo, S.; Dong, Z.; Wang, C.; Li, W. Theoretical and experimental study on stress-dependency of oil-water relative permeability in fractal porous media. Fractals 2018, 26, 1840010. [Google Scholar] [CrossRef]
  30. Yu, B.; Li, J. Some fractal characters of porous media. Fractals 2001, 9, 365–372. [Google Scholar] [CrossRef]
  31. Cai, J.; Yu, B.; Zou, M.; Luo, L. Fractal characterization of spontaneous co-current imbibition in porous media. Energy Fuels 2010, 24, 1860–1867. [Google Scholar] [CrossRef]
  32. Lei, G.; Dong, P.; Wu, Z.; Mo, S.; Gai, S.; Zhao, C.; Liu, Z.K. A fractal model for the stress-dependent permeability and relative permeability in tight sandstones. J. Can. Pet. Technol. 2015, 54, 36–48. [Google Scholar] [CrossRef]
  33. Cai, J.; Wei, W.; Hu, X.; Wood, D.A. Electrical conductivity models in saturated porous media: A review. Earth Sci. Rev. 2017, 171, 419–433. [Google Scholar] [CrossRef]
  34. Sheng, M.; Li, G.; Tian, S.; Huang, Z.; Chen, L. A fractal permeability model for shale matrix with multi-scale porous structure. Fractals 2016, 24, 1650002. [Google Scholar] [CrossRef]
  35. Lu, T.; Duan, Y.; Fang, Q.; Dai, X.; Wu, J. Analysis of fractional flow for transient two-phase flow in fractal porous medium. Fractals 2016, 24, 1650013. [Google Scholar] [CrossRef]
  36. Prakash, K.; Sridharan, A.; Prasanna, H.S. Dominant parameters controlling the permeability of compacted fine-grained soils. Indian Geotech. J. 2016, 46, 408–414. [Google Scholar] [CrossRef]
  37. Zhang, J.; Jiang, S.; Wang, Q.; Hou, Y.; Chen, Z. Critical hydraulic gradient of piping in sand. In Proceedings of the Twentieth International Offshore and Polar Engineering Conference, Beijing, China, 20–25 June 2010. [Google Scholar]
Figure 1. A comparison between the experimental data [36], and results of the proposed model.
Figure 1. A comparison between the experimental data [36], and results of the proposed model.
Energies 11 02986 g001
Figure 2. A comparison between the experimental data [37], and results of the proposed model.
Figure 2. A comparison between the experimental data [37], and results of the proposed model.
Energies 11 02986 g002
Figure 3. The average flow velocity curves: (a) for different threshold hydraulic gradient I; (b) for different differential order α.
Figure 3. The average flow velocity curves: (a) for different threshold hydraulic gradient I; (b) for different differential order α.
Energies 11 02986 g003
Figure 4. The average flow velocity curves versus hydraulic gradient with different effective stress.
Figure 4. The average flow velocity curves versus hydraulic gradient with different effective stress.
Energies 11 02986 g004
Figure 5. The average flow velocity curves: (a) for different rock elastic modulus; (b) for different power law index; (c) for different initial porosity.
Figure 5. The average flow velocity curves: (a) for different rock elastic modulus; (b) for different power law index; (c) for different initial porosity.
Energies 11 02986 g005aEnergies 11 02986 g005b

Share and Cite

MDPI and ACS Style

Lei, G.; Cao, N.; Liu, D.; Wang, H. A Non-Linear Flow Model for Porous Media Based on Conformable Derivative Approach. Energies 2018, 11, 2986. https://0-doi-org.brum.beds.ac.uk/10.3390/en11112986

AMA Style

Lei G, Cao N, Liu D, Wang H. A Non-Linear Flow Model for Porous Media Based on Conformable Derivative Approach. Energies. 2018; 11(11):2986. https://0-doi-org.brum.beds.ac.uk/10.3390/en11112986

Chicago/Turabian Style

Lei, Gang, Nai Cao, Di Liu, and Huijie Wang. 2018. "A Non-Linear Flow Model for Porous Media Based on Conformable Derivative Approach" Energies 11, no. 11: 2986. https://0-doi-org.brum.beds.ac.uk/10.3390/en11112986

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop