Next Article in Journal
Radially Symmetric Solution for Fractional Laplacian Systems with Different Negative Powers
Previous Article in Journal
Asymptotic and Finite-Time Synchronization of Fractional-Order Memristor-Based Inertial Neural Networks with Time-Varying Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Permeability Prediction of Saturated Geomaterials with Revised Pore–Solid Fractal Model and Critical Path Analysis

1
School of Water Conservancy and Civil Engineering, Zhengzhou University, Zhengzhou 450002, China
2
School of Mathematics (Zhuhai), Sun Yat-sen University, Zhuhai 519082, China
*
Author to whom correspondence should be addressed.
Submission received: 21 April 2022 / Revised: 7 June 2022 / Accepted: 16 June 2022 / Published: 23 June 2022
(This article belongs to the Topic Mathematical Modeling in Physical Sciences)

Abstract

:
The revised pore–solid fractal (PSF) model is presented by using the largest inscribed circle-based geometries of squares or cubes to replace the original pore or solid subregions as the new pore or solid phase in porous media. The revised PSF model changes the discrete lacunar pore and solid phases in the original PSF model to integrated. Permeability is an intrinsic property of geomaterials and has broad applications in exploring fluid flow and species transport. Based on the revised PSF model and critical path analysis, a fractal model for predicting the permeability of saturated geomaterials is proposed. The permeability prediction model is verified by comparison with the existing predicted model and the laboratory testing. The results show that the predicted permeabilities match the measured values very well. This work provides a theoretical framework for the revised PSF model and its application in predicting the permeability of geomaterials.

1. Introduction

Accurate prediction of the permeability for porous media such as geomaterials is an interesting topic in hydrology, geophysics, petroleum engineering, geotechnical engineering, etc. Permeability is the reflection and measurement of the complex relationship between the geometry and the topology of the pore phase in porous media, and it has been commonly used in exploring fluid flow and species transportation in porous media. The literature on predicting permeability is extensive. In terms of experimental testing, mercury intrusion porosimetry [1,2] and nuclear magnetic resonance [3,4] have been widely used in quantifying the pore properties of a porous medium and have become the effective verification methods for the theoretical models of estimating permeability. In terms of calculation methods, empirical models [5,6,7,8] and theoretical models considering porosity and pore or solid particle size distribution have been proposed; examples include percolation theory [9,10,11,12], effective medium approximation [11,12,13], critical path analysis [14,15] and fractal geometry [16]. Critical path analysis (CPA) can explore both throat connectivity and tortuosity of pore channels in porous media, which has been used to estimate the permeability of heterogeneous geomaterials [17,18,19,20].
The micro geometric structures of porous media have self-similar properties that are better interpreted and described with fractal geometry. A porous medium can be regarded as a fractal body, and its micro geometric structures can be analyzed by fractal models such as a solid or pore fractal model (Sierpinski carpet or Menger sponge) [21], pore–solid fractal model [22], bunched capillaries fractal model [16,23,24], multifractal model [25] and intermingled fractal units model [26]. The bunched capillaries fractal model [16,23,24] is mostly applied to estimate the permeability of a porous medium. This model considers the porous medium to consist of a bunch of cylindrical pore capillaries of different diameters; the cumulative pore size distribution and the length of capillaries follow the fractal scaling law. However, the permeability prediction models based on bunched capillaries fractal model lack a characterization of the pore interconnections in real existence in a porous medium, which is an idealization of a natural porous medium.
The pore or solid fractal model based on the standard two-dimensional Sierpinski carpet or three-dimensional Menger sponge has been mainly applied to focus on the cumulative pore or solid particle size distribution, the interface between the pore and solid, the mass (volume) and the water retention curve. However, based on the mono-fractal structure, these simplified fractal models cannot reflect the real porous media, especially geomaterials.
A geometrical, multiscale model of soil structure, which is named a pore–solid fractal (PSF) model [22], characterizes the structure of soil exhibits self-similarity to a degree; that is, where local structure appears, it is regarded as similar to the whole structure. It is achieved through the incorporation of a pore phase, solid phase and fractal phase that can be decomposed into smaller pore, solid and fractal subregions when iteration increases. The pore phase and the solid phase both have volume fractions greater than 0 and less than 1 even for infinite iterations, and they coexist at each level in the hierarchy, which is a main difference from the classical pore or solid fractal model with only the pore and solid phases based on the standard Sierpinski carpet or Menger sponge. In addition, when the fractal phase is considered as a solid phase or pore phase, the PSF model is simplified to a solid fractal model or a pore fractal model, respectively. At present, the PSF model is mostly used to estimate the soil water retention function [27,28] and determine the unsaturated hydraulic conductivity [29].
In this article, a fractal model for predicting the permeability of saturated geomaterials is proposed. Section 2 establishes a PSF model with circular-based pores to reflect the pore phase of different geometric pores and a PSF model with circular-based solid particles to describe the solid phase of different geometric solid particles. The permeability for the fragmented pore or solid phase is derived using the Kozeny–Carman equation, and the permeability for continuous geomaterials is predicted by combining the PSF model with circular-based pores and critical path analysis in Section 3. The revised PSF models degenerate to the general solid or pore fractal model, and the predicted model is verified by comparison with the existing predicted model and the laboratory testing in Section 4. The conclusion is expressed in Section 5.

2. Revised PSF Model

The original PSF model uses a square (or cube) as a basic unit that is the same as the pore or solid fractal model, which makes the contact between the pore and the solid complicated and makes theoretical analysis and numerical simulation difficult. In order to improve the applicability of the PSF model, it should be revised.

2.1. PSF Model with Circular-Based Pores

The model construction is consistent with the original PSF model. Firstly, divide a region of size L into N equal subregions in the Euclidean space of dimension d as the initiator. N is calculated by
N = n d .
n = L / r .
where r is the linear length of each subregion and n is the number of divisions in each direction.
Define Nz subregions as the self-similar replicated phase, Ns subregions as the solid phase of porous media and Np subregions as the pore phase.
p + s + z = 1 .
where p , s and z are the proportions of the pore, solid and self-similar replicated phases in the initiator, respectively.
Take the largest inscribed circles (cylinders) in the Np subregions as the new pore phase and consider the remaining parts of the Np subregions as the solid phase. The above operations are performed in 1st generation as a generator. In the next ith (i 2 ) generation, a recursive process replaces ( N z ) i 1 subregions with the reduced replicate of the generator with the ratio ( 1 / n ) i 1 . An example of the construction process of the PSF model with circular-based pores is illustrated in Figure 1.
The cumulative pore size distribution P ( r i p ) of the PSF model with circular-based pores is still calculated using the corresponding equation of the original PSF model [22]:
P ( r i p ) = P i 1 + ( N z ) i 1 N p = N p ( N z ) i 1 N z 1 N p N z 1 L D ( r i p ) D .
where D is the fractal dimension of the original PSF model; P ( r i p ) is the cumulative number of pores after ith generation; r i p is the pore diameter after ith generation, r i p = L n i .
An inscribed circle (cylinder) of maximum diameter in a square (cube) subregion of size r i instead of the square (cube) pore subregion leads to a decrease in porosity as
ϑ = r i d α r i d r i d = 1 α .
where α is geometric parameter, circle (cylinder) α = π / 4 , semicircle (semicylinder) α = π / 8 and square (cube) α = 1 .
At the end of the ith generation, the decreased porosity is
d ϕ i = d ϕ i 1 + p ϑ z i 1 = p ϑ j = 0 i 1 z j = p ϑ p + s ( 1 z i ) .
Therefore, the porosity ϕ i of the PSF model with circular-based pores is shown as
ϕ i = p p + s ( 1 z i ) p ϑ p + s ( 1 z i ) = p α p + s [ 1 ( r i p L ) d D ] .
When i , r i p 0 , and then ϕ = p α p + s .
The surface area of circular-based pores S p ( i ) is obtained by
S p ( i ) = S p ( i 1 ) + π N p ( N z ) i 1 ( L n i ) d 1 = π N p ( L n ) d 1 j = 1 i ( N z n d 1 ) j 1 .
Noting N = n d , Equation (6) changes to
S p ( i ) = π p n L d 1 j = 1 i ( n z ) j 1 .
Since the construction process of the PSF model with circular-based pores is consistent with the original PSF, refer to the original PSF model (Perrier et al., 1999 [22]); the pore–solid interface area S s p ( i ) is achieved by
S s p ( i ) = ( 1 ϕ ) S p ( i ) = ϕ S s ( i ) .
So, the surface area of the solid phase S s ( i ) is written as
S s ( i ) = 1 ϕ ϕ S p ( i ) = p ϑ + s p α S p ( i ) .
Assume M s ( r i p ) is the cumulative mass of solid phase when the pore diameter is less than or equal to r i p , which is obtained by multiplying the solid volume of pore diameter smaller than or equal to r i p and solid density ρ s . After T 1 iterations (Tth generation) of the PSF model with circular-based pores, M s ( r i p ) is derived as follows:
M s ( r i p ) = ρ s L d [ s + ( 1 α ) p ] j = i T z j 1 = ρ s L d s + ( 1 α ) p s + p ( z i 1 z T ) = ρ s L d s + ( 1 α ) p s + p [ ( r i p r max p ) d D ( r T p L ) d D ] .
When T , r T p 0 , and then
M s ( r i p ) = s + ( 1 α ) p s + p ρ s L d ( r i p r max p ) d D .

2.2. PSF Model with Circular-Based Solid Particles

The construction process is similar to the PSF model of circular-based pores except that the largest inscribed circles (spheres) in the N s ( N z ) i 1 subregions of size r i are used to represent the new solid phase and the remaining parts of the N s ( N z ) i 1 subregions of size r i are regarded as the pore phase.
The cumulative solid particle size distribution S ( r i s ) of the PSF model with circular-based solid particles is still calculated as
S ( r i s ) = S i 1 + ( N z ) i 1 N s = N s ( N z ) i 1 N z 1 N s N z 1 L D ( r i s ) D .
where S ( r i s ) is the cumulative number of solid particles after ith generation; r i s is the solid particle diameter after ith generation, r i s = L n i .
An inscribed circle (sphere) of maximum diameter in a square (cube) subregion of size r i replaces the square (cube) solid subregion, which causes an increase in porosity as
υ = r i d β r i d r i d = 1 β .
where β is a geometric parameter. When the subregion is a circle, β = π / 4 ; sphere, β = π / 6 ; square (cube), β = 1 .
At the end of the ith generation, the increased porosity is
d ϕ i = d ϕ i 1 + s υ z i 1 = s υ j = 0 i 1 z j = s υ p + s ( 1 z i ) .
Therefore, the porosity ϕ i of the PSF model with circular-based solid particles is obtained as
ϕ i = p p + s ( 1 z i ) + s υ p + s ( 1 z i ) = p + s υ p + s [ 1 ( r i s L ) d D ] .
when i , r i s 0 , and then ϕ = p + s υ p + s .
According to the PSF model with circular-based pores, the surface area of circular-based solid particles S s ( i ) is written as
S s ( i ) = π s n L d 1 j = 1 i ( n z ) j 1 .
The surface area of the pore phase S p ( i ) is obtained as
S p ( i ) = ϕ 1 ϕ S s ( i ) = p + s υ s β S s ( i )
The cumulative mass of the solid phase M s ( r i s ) when the solid particle diameter is less than or equal to r i s is derived as
M s ( r i s ) = ρ s s β L d j = i r z i 1 = s s + p ρ s β L d [ ( r i s r max s ) d D ( r T s L ) d D ]
When T , then
M s ( r i s ) = s β s + p ρ s L d ( r i s r max s ) d D .
When the subregions of solid phase are squares (cubes), β = 1 . Equation (20) coincides with the cumulative mass of the solid phase given by Bird et al. [30].

3. Permeability for Revised PSF Model

3.1. Fragmented Pore or Solid Phase

The Kozeny–Carman equation is widely used to estimate the permeability of porous media [5,6], which is expressed as
k = ϕ 3 C ( 1 ϕ ) 2 S 2 .
where k is the permeability, ϕ is the porosity, C is the empirical parameter and S is the specific surface area per unit volume of the solid phase.
From the revised PSF model construction, a porous medium can be fragmented into circular-based pores of different sizes embedded in an aggregated solid phase or circular-based solid particles surrounded by pore space.
For circular-based pores, according to Equation (11), the specific surface area per unit volume of solid phase in ith iteration S i is
S i = S s ( i ) V = π ( 1 ϕ ) p n ϕ L j = 1 i ( n z ) j 1 .
If n z = 1 ,
S i = π ( 1 ϕ ) p n i ϕ L   k i = ϕ 5 ( r max p ) 2 C ( 1 ϕ ) 4 ( π p i ) 2 .
If n z 1 ,
S i = π p n L ( n z ) i 1 n z 1   k i = ϕ 5 ( r max p ) 2 ( n z 1 ) 2 C ( 1 ϕ ) 4 ( π p ) 2 [ ( n z ) i 1 ] 2 .
For circular-based solid particles, based on Equation (18), S i is obtained by
S i = S s ( i ) V = π s n L j = 1 i ( n z ) j 1 .
If n z = 1 ,
S i = π s n i L   k i = ϕ 3 ( r max s ) 2 C ( 1 ϕ ) 2 ( π s i ) 2 .
If n z 1 ,
S i = π s n L ( n z ) i 1 n z 1   k i = ϕ 3 ( r max s ) 2 ( n z 1 ) 2 C ( 1 ϕ ) 2 ( π s ) 2 [ ( n z ) i 1 ] 2 .

3.2. Continuous Geomaterials

Equations (4) and (13) show that the cumulative discrete pore size distribution and solid particle size distribution of the revised PSF model obey P r D and S r D , respectively. Hunt argues that pore size is commonly used to characterize continuous porous media and can be translated into the case of using solid particles [31,32]. Learning from this, the continuous probability density function of pore size f ( r ) is obtained as
f ( r ) = α p ( d D ) ( p + s ) r max d D r 1 D .
where r max is the maximum pore diameter.
When α = 1 , Equation (29) degenerates to the formula proposed by Ghanbarian et al. [11,29].
The cumulative probability distribution function for pore diameter smaller than r is achieved by integrating Equation (29) as
V ( A ) = α p ( d D ) ( p + s ) r max d D r min r r d 1 D d r = α p p + s [ ( r r max ) d D ( r min r max ) d D ] .
If r = r max ,
V ( A ) = α p p + s [ 1 ( r min r max ) d D ] .
It is obvious that Equation (31) is the equivalent of Equation (7).
The water retention function θ is obtained with the process of the original PSF model [30] as
θ = ϕ α p p + s [ 1 ( h h min ) d D p ] h min h h max .
where h min and h max are the water tension heads at which the maximum and minimum pore diameters occur in the porous media, respectively.
According to critical path analysis [17,19,20], the percolation threshold p c of pore volume for saturated condition ( θ = ϕ ) is defined as
p c = r c r max f ( r ) d r = α p ( p + s ) [ 1 ( r c r max ) 3 D ] .
where r c is the critical pore diameter for saturated porous media.
Transform Equation (33) to obtain
r c = r max [ 1 p c ( p + s ) α p ] 1 3 D .
The percolation threshold p c and the critical pore diameter r c for an unsaturated condition ( θ ϕ ) are
p c = r c r f ( r ) d r = α p ( p + s ) [ ( r r max ) 3 D ( r c r max ) 3 D ] .
and
r c = [ r 3 D p c ( p + s ) α p r max 3 D ] 1 3 D .
The generalized equation for permeability k derived by CPA is expressed as [17,18,20,29,33]
k = γ c 2 4 C F .
where F is the formation factor and C is a constant according to the specific porous materials.
The formation factor F of saturated porous media is approximated from the water retention curve for the saturated condition ( θ = ϕ ) and is written as [20,33]
1 F = ϕ 3 [ 1 ( r c 3 r max ) 3 D ] .
Equation (38) has been successfully used for predicting the permeability of geomaterials [20,33].
Substitute the critical pore diameter Equation (34) and the formation factor Equation (38) into Equation (37) to yield the permeability as
k = r max 2 ϕ 12 C [ 1 ( 1 3 ) 3 D ( 1 p c ( p + s ) α p ) ] [ 1 p c ( p + s ) α p ] 2 3 D .

4. Discussion and Results

4.1. General Pore or Solid Fractal Model

When s = 0 , p + z = 1 , and the PSF model with circular-based pores degenerates to the general solid fractal model. P and z are the proportions of pore and solid. An example of a general solid fractal model with circular-based pores is illustrated in Figure 2.
The porosity of the solid fractal model with circular-based pores is
ϕ i = α ( 1 z i ) = α [ 1 ( r i L ) d D ] .
where D is the fractal dimension of the general solid fractal model.
When i , ϕ = α .
For the solid fractal model based on a standard two-dimensional Sierpinski carpet or three-dimensional Menger sponge, similarly, take the largest inscribed circles (cylinders) in the square (cube) pore subregions of size r i as the new pore phase; the porosity is
ϕ i = α ( 1 ( 8 9 ) i ) = α ( 1 z i ) .
It can clearly be seen that the solid fractal model with circular-based pores based on standard Sierpinski carpet or Menger sponge is a special case of the general solid fractal model changed from PSF.
When p = 0 , s + z = 1 , and the PSF model with circular-based solid particles changes to the general pore fractal model. s and z are the proportions of solid and pore. The porosity of the pore fractal model with circular-based solid particles is
ϕ i = υ ( 1 z i ) + z i = 1 β + β z i = 1 β + β ( r i L ) d D .
where D is the fractal dimension of the general pore fractal model.
When i , ϕ = υ = 1 β .
For the pore fractal model based on a Sierpinski carpet or Menger sponge, the largest inscribed circles (spheres) in the square (cube) solid subregions of size r i are used to represent the new solid phase; the porosity is
ϕ i = υ [ 1 ( 8 9 ) i ] + ( 8 9 ) i = υ ( 1 z i ) + z i .
Equation (43) coincides with the Khabbazi formula [34].
Obviously, the pore fractal model with circular-based solid particles based on a Sierpinski carpet or Menger sponge is also a special case of the general pore fractal model transformed from PSF.

4.2. Permeability

Equation (39) is applied to predict the permeabilities of the geomaterials samples used in Daigle’s prediction [17]. The comparison of predicted permeabilities and the measured values is shown in Figure 3.
It can be seen from Figure 3 that the predicted permeabilities with Equation (36) are almost consistent with the measured permeabilities and k m e a k p r e = ± 0.4 × 10 11 is the upper and lower boundary. Meanwhile, the ratio k p r e k m e a is between 0.9 and 1.1 when log ( r max r min ) 2.7 ; for log ( r max r min ) < 2.7 , the error between the predicted permeabilities and the measured values is 20% to 10%.
The difference between the predicted values in this article and Daigle’s predictions is mainly due to the difference in the calculation method of formation factor F. Therefore, the accurate expression of formation factor F determines the accuracy of permeability prediction.
Figure 4 shows the predicted permeability k p r e is highly correlated with the critical pore diameter r c . The permeability is not only influenced by micro geometric characteristics such as pore diameter, pore size distribution and interface area of pores and solid particles, but also influenced by topological properties such as connectivity between pore and solid phases in the porous medium, which were also obtained by Ghanbarian et al. [20].

5. Conclusions

The revised PSF models use the largest inscribed circle-based geometries in pore or solid subregions instead of original squares or cubes to represent pores or solid particles in porous media, which increases the integrated contact between pores and the solid phase.
A fractal model for predicting the permeability of geomaterials based on the PSF model with circular-based pores and critical path analysis is proposed; the permeability is highly correlated with the critical pore diameter.

Author Contributions

Funding acquisition, J.W.; Conceptualization, J.W.; Methodology, L.K.; Supervision, J.W.; validation, L.K.; Writing—original draft, L.K.; Writing—review and editing, L.K. and W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Fundamental Research Funds for the Central Universities of Sun Yat-sen University (Grant no. 22qntd290) and the National Natural Science Foundation of China (Grant no. 51708512).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Katz, A.J.; Thompson, A.H. Quantitative prediction of permeability in porous rock. Phys. Rev. B 1986, 34, 8179–8181. [Google Scholar] [CrossRef]
  2. Ghanbarian, B.; Torres-Verdín, C.; Lake, L.W.; Marder, M. Gas permeability in unconventional tight sandstones: Scaling up from pore to core. J. Pet. Sci. Eng. 2018, 173, 1163–1172. [Google Scholar] [CrossRef] [Green Version]
  3. Coates, G.; Xiao, L.; Prammer, M. NMR Logging: Principles and Applications; Haliburton Energy Services: Houston, TX, USA, 1999; Volume 344. [Google Scholar]
  4. Glover, P.W.; Zadjali, I.I.; Frew, K.A. Permeability prediction from MICP and NMR data using an electrokinetic approach. Geophysics 2006, 71, F49–F60. [Google Scholar] [CrossRef]
  5. Kozeny, J. Ueber kapillare leitung des wassers im boden. Stizungsber. Akad. Wien Akad. Wiss 1927, 136, 271–306. [Google Scholar]
  6. Carman, P.C. Permeability of saturated sands, soils and clays. J. Agric. Sci. 1939, 29, 262–273. [Google Scholar] [CrossRef]
  7. Brutsaert, W. Some Methods of Calculating Unsaturated Permeability. Trans. ASAE 1967, 10, 400–404. [Google Scholar] [CrossRef]
  8. Mishra, S.; Parker, J.C. On the Relation Between Saturated Conductivity and Capillary Retention Characteristics. Ground Water 1990, 28, 775–777. [Google Scholar] [CrossRef]
  9. Sahimi, M. Applications of Percolation Theory; Taylor and Francis: London, UK, 1994. [Google Scholar]
  10. Stauffer, D.; Aharony, A. Introduction to Percolation Theory; Taylor and Francis: London, UK, 1994. [Google Scholar]
  11. Ghanbarian, B.; Hunt, A.; Skinner, T.; Ewing, R. Saturation dependence of transport in porous media predicted by percolation and effective medium theories. Fractals 2015, 23, 1540004. [Google Scholar] [CrossRef]
  12. Richesson, S.; Sahimi, M. Hertz-Mindlin Theory of Contacting Grains and the Effective-Medium Approximation for the Permeability of Deforming Porous Media. Geophys. Res. Lett. 2019, 46, 8039–8045. [Google Scholar] [CrossRef]
  13. Sahimi, M. Heterogeneous Materials I: Linear Transport and Optical Properties; Springer: New York, NY, USA, 2003. [Google Scholar]
  14. Hunt, A. Applications of percolation theory to porous media with distributed local conductances. Adv. Water Resour. 2001, 24, 279–307. [Google Scholar] [CrossRef]
  15. Hunt, A. Comments on “Fractal Fragmentation, Soil Porosity, and Soil Water Properties: I. Theory”. Soil Sci. Soc. Am. J. 2007, 71, 1418–1419. [Google Scholar] [CrossRef]
  16. Cai, J.; Luo, L.; Ye, R.; Zeng, X.; Zhang, L.; Hu, Z. Recent advances on fractal modeling of permeability for fibrous porous media. Fractals 2015, 23, 1540006. [Google Scholar] [CrossRef]
  17. Daigle, H. Application of critical path analysis for permeability prediction in natural porous media. Adv. Water Resour. 2016, 96, 43–54. [Google Scholar] [CrossRef]
  18. Skaggs, T. Assessment of critical path analyses of the relationship between permeability and electrical conductivity of pore networks. Adv. Water Resour. 2011, 34, 1335–1342. [Google Scholar] [CrossRef]
  19. Hunt, A.; Ghanbarian, B.; Saville, K.C. Unsaturated hydraulic conductivity modeling for porous media with two fractal regimes. Geoderma 2013, 207–208, 268–278. [Google Scholar] [CrossRef]
  20. Ghanbarian, B.; Hunt, A.; Skaggs, T.H.; Jarvis, N. Upscaling soil saturated hydraulic conductivity from pore throat characteristics. Adv. Water Resour. 2017, 104, 105–113. [Google Scholar] [CrossRef] [Green Version]
  21. Sergeyev, Y. Evaluating the exact infinitesimal values of area of Sierpinski’s carpet and volume of Menger’s sponge. Chaos Solitons Fractals 2009, 42, 3042–3046. [Google Scholar] [CrossRef] [Green Version]
  22. Perrier, E.; Bird, N.; Rieu, M. Generalizing the fractal model of soil structure: The pore–solid fractal approach. Geoderma 1999, 88, 137–164. [Google Scholar] [CrossRef]
  23. Yu, B.; Cheng, P. A fractal permeability model for bi-dispersed porous media. Int. J. Heat Mass Transf. 2002, 45, 2983–2993. [Google Scholar] [CrossRef]
  24. Kou, L.; Miao, R.H.; Miao, F.Y. Fractal analysis of non-Newton fluid grouting through soil composed of arbitrary cross-sectional capillaries. Fractals 2021, 29, 2150139. [Google Scholar] [CrossRef]
  25. Perfect, E.; Gentry, R.; Sukop, M.; Lawson, J. Multifractal Sierpinski carpets: Theory and application to upscaling effective saturated hydraulic conductivity. Geoderma 2006, 134, 240–252. [Google Scholar] [CrossRef]
  26. Pia, G.; Sanna, U. An intermingled fractal units model and method to predict permeability in porous rock. Int. J. Eng. Sci. 2014, 75, 31–39. [Google Scholar] [CrossRef]
  27. Huang, G.; Zhang, R. Evaluation of soil water retention curve with the pore–solid fractal model. Geoderma 2005, 127, 52–61. [Google Scholar] [CrossRef]
  28. Hwang, S.I.; Yun, E.Y.; Ro, H.M. Estimation of soil water retention function based on asymmetry between particle- and pore-size distributions. Eur. J. Soil Sci. 2011, 62, 195–205. [Google Scholar] [CrossRef]
  29. Ghanbarian-Alavijeh, B.; Hunt, A. Unsaturated hydraulic conductivity in porous media: Percolation theory. Geoderma 2012, 187–188, 77–84. [Google Scholar] [CrossRef]
  30. Bird, N.; Perrier, E.; Rieu, M. The water retention function for a model of soil structure with pore and solid fractal distributions. Eur. J. Soil Sci. 2000, 51, 55–63. [Google Scholar] [CrossRef]
  31. Hunt, A. Continuum percolation theory for transport properties in porous media. Philos. Mag. 2005, 85, 3409–3434. [Google Scholar] [CrossRef]
  32. Hunt, A.; Sahimi, M. Flow, Transport, and Reaction in Porous Media: Percolation Scaling, Critical-Path Analysis, and Effective Medium Approximation. Rev. Geophys. 2017, 55, 993–1078. [Google Scholar] [CrossRef]
  33. Ghanbarian, B. Predicting Single-Phase Permeability of Porous Media Using Critical-Path Analysis. In Complex Media and Percolation Theory; Sahimi, M., Hunt, A.G., Eds.; Encyclopedia of Complexity and Systems Science Series; Springer: New York, NY, USA, 2021. [Google Scholar]
  34. Khabbazi, A.E.; Hinebaugh, J.; Bazylak, A. Analytical tortuosity–porosity correlations for Sierpinski carpet fractal geometries. Chaos Solitons Fractals 2015, 78, 124–133. [Google Scholar] [CrossRef]
Figure 1. Construction process of PSF model with circular-based pores: (a) 0th generation; (b) 1st generation; (c) 2nd generation.
Figure 1. Construction process of PSF model with circular-based pores: (a) 0th generation; (b) 1st generation; (c) 2nd generation.
Fractalfract 06 00351 g001
Figure 2. Construction process of solid fractal model with circular-based pores: (a) 0th generation; (b) 1st generation; (c) 2nd generation.
Figure 2. Construction process of solid fractal model with circular-based pores: (a) 0th generation; (b) 1st generation; (c) 2nd generation.
Fractalfract 06 00351 g002
Figure 3. (a) Comparison of predicted permeabilities, measured values and existing predicted values. (b) Distribution of the ratio of predicted to measured permeabilities.
Figure 3. (a) Comparison of predicted permeabilities, measured values and existing predicted values. (b) Distribution of the ratio of predicted to measured permeabilities.
Fractalfract 06 00351 g003
Figure 4. Relationship of predicted permeability and critical pore diameter.
Figure 4. Relationship of predicted permeability and critical pore diameter.
Fractalfract 06 00351 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kou, L.; Li, W.; Wu, J. Permeability Prediction of Saturated Geomaterials with Revised Pore–Solid Fractal Model and Critical Path Analysis. Fractal Fract. 2022, 6, 351. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070351

AMA Style

Kou L, Li W, Wu J. Permeability Prediction of Saturated Geomaterials with Revised Pore–Solid Fractal Model and Critical Path Analysis. Fractal and Fractional. 2022; 6(7):351. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070351

Chicago/Turabian Style

Kou, Lei, Wuxue Li, and Jujie Wu. 2022. "Permeability Prediction of Saturated Geomaterials with Revised Pore–Solid Fractal Model and Critical Path Analysis" Fractal and Fractional 6, no. 7: 351. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6070351

Article Metrics

Back to TopTop