# Fractal Analysis of Permeability of Unsaturated Fractured Rocks

Authors: Guoping Jiang; Wei Shi; Lili Huang
Journal: The Scientific World Journal
 (2013)

---

## Abstract

A physical conceptual model for water retention in fractured rocks is derived while taking into account the effect of pore size distribution and tortuosity of capillaries. The formula of calculating relative hydraulic conductivity of fractured rock is given based on fractal theory. It is an issue to choose an appropriate capillary pressure-saturation curve in the research of unsaturated fractured mass. The geometric pattern of the fracture bulk is described based on the fractal distribution of tortuosity. The resulting water content expression is then used to estimate the unsaturated hydraulic conductivity of the fractured medium based on the well-known model of Burdine. It is found that for large enough ranges of fracture apertures the new constitutive model converges to the empirical Brooks-Corey model.


---

## Body

## 1. Introduction

Modeling water  flow in unsaturated fractured rocks has received considerable attention in the last two decades. One of the main reasons for focussing on the study of liquid flow in this type of media is that spontaneous capillary imbibition is an important fundamental phenomenon existing extensively in a variety of processes such as oil recovery, polymer composite manufacturing, soil science, and hydrology. The other reason is that deep disposal in crystalline rocks is considered to be an effective mean of isolating radioactive wastes from the biosphere.  The study of basic transport processes has long been recognized because of  which one  can heighten the comprehensive understanding of physical phenomena, such as permeability [1, 2], heat transfer [3–5], and sorption [6, 7]. Many parameters such as the porosity, size of pore, and tortuosity of capillaries are very important for fluid flow in hard rocks. These parameters, however, are closely related to the geometric architecture of hard rocks.  Cai et al.  have analyzed the natural fractured trace maps representing a wide variety of scales, geological settings, and lithologies [8].  Cai and Yu  have reported the density of different sized fault segments within the San Andreas fault zone which is fractal [9]. The distributions of fracture aperture and fracture spacing are self-similar over a well-defined range of apertures in the Cajon Pass scientific drill hole which was found by Barton and Zoback [10]. Based on the assumption that the fracture pattern is self-similar, the Sierpinski carpet was often employed to simulate porous media [11]. The conceptual constitutive model proposed by  Wu and Yu  [12] had used its fractal dimension to the parameters of the Brooks-Corey constitutive model [13] through the Sierpinski carpet. A Sierpinski space was also adopted to characterize the spatial distribution of a drainage network in the Gardon basin, France [14].From the above review it is shown that a mechanistic model has not yet been established. In this paper, we attempt to develop a mechanistic model for unsaturated flow in fractured hard rocks based on the method using the specific fractal to describe fractured rock. The expressions of the proposed constitutive model are closed form and easy to evaluate. Another  important feature,  the tortuosity fractal dimensionDT,  which  affected  other model parameters and should not be neglected  while it was  often neglected in the past investigation, is considered.Now fractal theory has a wide variety of applications in sciences and engineering fields such as thermal science [15–19], fluid science [1, 2, 6], and industrial construction engineering [20]. For example,  Moussa  [2] has systematically investigated the transport of porous media based on fractal theory. The work [2] is open. Xiao et al. [3–5, 15, 16, 18, 19] have done much outstanding work on heat transfer of fluids by using fractal technique. In our work, we derive the analytical expressions for the relative hydraulic conductivity of fractured rock while taking into account the effect of pore size distribution based on the fractal geometry theory.


## 2. Construction of the Fractal Model

The model is presented in Figure1. It has been shown that the cumulative size distribution of contact spots on engineering surfaces is similar to islands on earth and pores in porous rock which follows the fractal scaling law [2]. 						(1)N(L≥Dn)=(Dn,maxDn)df,
where df is the fractal dimension for pores, D is diameter, L is the length scale, and N is the total number of pores whose sizes equal to and greater than Dn. The number of pores whose sizes range from Dn to D+dDn is
(2)-dN=dfDn,maxdfDn-(df+1)dDn
when water flow through the pores of porous rock, the capillaries may be tortuous. These tortuous capillaries could be expressed by fractal equation [21].  Consider
(3)La(Dn)=L0DTDn1-DT,
where DT is the tortuosity fractal dimension and lies in the range 1<DT<2, which represents the extent of convoluted ness of capillary pathways for fluid flow through a medium. Note that for a straight capillary path DT=1, and a higher value of  DT  corresponds to a highly tortuous capillary. Let the diameter of a capillary in the medium be Dn and let its tortuous length along the flow direction be La(Dn). L0 is representative length of channels. With a straight capillary, La(Dn) = L0. The total volume of  pores from Dn,min to Dn,max can be obtained from (1) as
(4)vtol=∫Dn,minDn,maxπ4Dn2La(-dN)=∫Dn,minDn,maxπ4Dn2LadfDn·maxdfDn-(df+1)dDn=π4(3-df-DT)L0DTdfDn·max3-DT×(1-(Dn,minDn,max)(3-df-DT)),
where vtol is the whole volume of  the pores, (-dN) is given by (2).Figure 1
The model presented.Similarity to the state above, the volumev(Dn,minx<Dn)of  pores/particles from Dn to Dn,max can be obtained from (1) as
(5)v(Dn,min<Dn)=π4L0DTdfDn·maxdf∫Dn,minxDn,‍Dn2-df-DTdDn=π4(3-df-DT)L0DTdfDn·maxdf(Dn3-df-DT)|Dn,minDn,.
The effective saturation S of volume v(Dn,minx<Dn) can be obtained. Consider
(6)S=v(Dn,min<Dn)vtol=|Dn,minDn,max(Dn3-df-DT)|Dn,minDn,(Dn3-df-DT)=Dn3-df-DT-Dn,min3-df-DTDn,max3-df-DT-Dn,min3-df-DT.Assuming the immiscible fluid flow in reservoir rocks, three relatively important forces are considered: the capillary pressure can be expressed. Consider(7)p=2σcos(β)ρgDn,
where σ is surface tension between the wetting and non-wetting fluids, ρ is the water density, g is the gravity acceleration, Dn is diameter of a pore, and β is contact angle between the extraneous water and solid.The relative permeabilities are usually expressed in terms of water saturationSw. The saturation curve for the proposed model of fractured rock was derived. Consider
(8)Sw=pc3-df-DT-pcmin3-df-DTpc,min3-df-DT-pc,max3-df-DT,
where
(9)pcmax=2σcos(β)ρgDn,min,pcminx=2σcos(β)ρgDn,max.


## 3. The Permeability of Unsaturated Fractured Rocks

The Burdine and Mualem models are the two most widely used models predicting relative hydraulic conductivity. For the particular state of flow in fractured hard rock, the Burdine model seems to be more consistent and to be often adopted. The expression of the Burdine model is(10)KD(s)=s2∫0s(ds/pc2)∫01(ds/pc2),
where KD(s) is the relative hydraulic conductivity. Inserting (8) to (10) we obtain the following form for KD(s):(11)KD(s)=(pcdf+DT-3-pcmaxdf+DT-3pc,minxdf+DT-3-pcmaxdf+DT-3)2×pcdf+DT-5-pcmaxdf+DT-5pcminxdf+DT-5-pcmaxdf+DT-5.
The expressions of (6) and (10) represent the proposed constitutive model for fractured hard rocks. Note that all model parameters are determined by geometric parameters of Dn,max,,Dn,minx and residual water content β.The novel constitutive model has some similarities with the well-known Brooks-Corey model, which isS0=(pc/pd)-λ and KD(s)=S3+λ/2,  where pd is the reciprocal of air entry pressure and λ is a model parameter related to pore size distribution.WhenDT=1, λ=2-df, and pc,max≪pc,minx, pc,max is neglected here, the model derived here is simplified to the Brooks-Corey model. Comparisons of the proposed and the Brooks-Corey models for three different ranges of fracture apertures are depicted in Figure 2. The assumed parameters are df=1.8, Dn,max=10-1 cm, and Dn,minx=10-3 cm, 10-5 cm, and 10-7 cm. According to Figure 2, the Brooks-Corey model seems to be adequate to describe the hydraulic properties of fractured rocks for large ranges of fracture apertures and low values of  pressure head.Comparison between the Brooks-Corey model and the new relations for different ranges of fracture apertures whenDT=1.
(a)
Effective saturation
(b)
Relative hydraulic conductivityThe relationship between the fractal dimension and effective saturation is very important in the study of unsaturated flow in fractured hard rocks. The geometric parameters and physical constants used for the analysis areβ=0,σ=72.25dy/pa. In order to analyze the influence of the fracture density we consider different values of DT. With (8) and (11), the relationship is obtained (Figure 3). With the physical constants Dn=0.1cm,  Dn·max=1cm, and Dn·minx=0.01cm, the relationship between df and S0 can also be determined (Figure 4).The new relations forKD(S) and DT.
(a)
D
n
,
min
=
10
-
3 cm
(b)
D
n
,
min
=
10
-
7 cmFigure 4
The new relations forS0  and df.In order to test the proposed invasion depth model, it is crucial to correctly determine the fractal dimensionsDT and df.According to the following formula proved recently by Yu and Li [22, 23],
(12)ϕ=(Dn,minDn,max)DE-df,
where DE is the Euclidean dimension and DE=2 is used in this work. The maximum pore diameter can be calculated based on the model of square arrangement of particles.According to the following formula proved recently by Jiangchao Cai, whereK is permeability, k is Kozeny constant, which considers the tortuosity of capillaries and pore non-uniformity. And the parameters K and k are all proved [23]. Consider
(13)DT=1+lntavlnL0/tav,Dn,max=2(1-ϕ)ϕKkϕ(ϕ1-ϕ+π4(1-ϕ)-1),
where
(14)tav=12[1+121-ϕ+(1-ϕ-1)2+((1-ϕ)/4)1-1-ϕ],L0tav=df-121-ϕϕπdf(2-df)Dn,minDn,max.
From (12) and (13), the relationship among ϕ, DT, and df is determined (Figure 5).Figure 5
The relationship amongϕ,DT,  and df.


## 4. Conclusions

In this paper, with the consideration of pore size distribution and tortuosity of capillaries, a new fractal model for relative hydraulic conductivity of fractured rock is developed. The derived constitutive model is an effort to understand and characterize unsaturated flow in fractured rocks. The expressions of water content and relative hydraulic conductivity curves have analytical closed forms. The parameters can be completely determined by the geometry of the fractal model. Every parameter of the proposed formulas of calculating relative hydraulic conductivity of fractured rock has a clear physical meaning. The fractal model can reveal the mechanisms of hydraulic conductivity for fluids flow through unsaturated rocks. The tortuosity fractal dimensionDT which affected other model parameters and should not be neglected while it was often neglected in the past investigation is considered in this paper.



---

*Source: 490320-2013-03-27.xml*