Fractal analysis of permeability of unsaturated fractured rocks.

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.

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 dimension D , 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 Figure 1. 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]. where d 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 D . The number of pores whose sizes range from D to D + dD is 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 where D is the tortuosity fractal dimension and lies in the range 1 < D < 2, which represents the extent of convoluted ness of capillary pathways for fluid flow through a medium. Note that for a straight capillary path D = 1, and a higher value of  D   corresponds to a highly tortuous capillary. Let the diameter of a capillary in the medium be D and let its tortuous length along the flow direction be L (D ). L 0 is representative length of channels. With a straight capillary, L (D ) = L 0. The total volume of pores from D to D can be obtained from (1) as where v tol is the whole volume of the pores, (−dN) is given by (2).

Figure 1

The model presented.

Similarity to the state above, the volume v (of pores/particles from D to D can be obtained from (1) as The effective saturation S of volume v ( can be obtained. Consider

Assuming the immiscible fluid flow in reservoir rocks, three relatively important forces are considered: the capillary pressure can be expressed. Consider where σ is surface tension between the wetting and non-wetting fluids, ρ is the water density, g is the gravity acceleration, D 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 saturation S . The saturation curve for the proposed model of fractured rock was derived. Consider where

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 where K (s) is the relative hydraulic conductivity. Inserting (8) to (10) we obtain the following form for K (s): 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 D , D and residual water content β.

The novel constitutive model has some similarities with the well-known Brooks-Corey model, which is S 0 = (p /p )− and K (s) = S 3+, where p is the reciprocal of air entry pressure and λ is a model parameter related to pore size distribution.

When D = 1, λ = 2 − d , and p ≪ p , p 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 d = 1.8, D = 10−1 cm, and D = 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.

Figure 2

Comparison between the Brooks-Corey model and the new relations for different ranges of fracture apertures when D = 1.

The 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.25 dy/pa. In order to analyze the influence of the fracture density we consider different values of D . With (8) and (11), the relationship is obtained (Figure 3). With the physical constants D = 0.1 cm,  D = 1 cm, and D = 0.01 cm, the relationship between d and S 0 can also be determined (Figure 4).

Figure 3

The new relations for K (S) and D .

Figure 4

The new relations for S 0 and d .

In order to test the proposed invasion depth model, it is crucial to correctly determine the fractal dimensions D and d .

According to the following formula proved recently by Yu and Li [22, 23], where D is the Euclidean dimension and D = 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, where K 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 where From (12) and (13), the relationship among ϕ, D , and d is determined (Figure 5).

Figure 5

The relationship among ϕ,D , and d .

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 dimension D T which affected other model parameters and should not be neglected while it was often neglected in the past investigation is considered in this paper.