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.
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 . Cai and Yu have reported the density of different sized fault segments within the San Andreas fault zone which is fractal . 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 . Based on the assumption that the fracture pattern is self-similar, the Sierpinski carpet was often employed to simulate porous media . The conceptual constitutive model proposed by Wu and Yu  had used its fractal dimension to the parameters of the Brooks-Corey constitutive model  through the Sierpinski carpet. A Sierpinski space was also adopted to characterize the spatial distribution of a drainage network in the Gardon basin, France .
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 , 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 . For example, Moussa  has systematically investigated the transport of porous media based on fractal theory. The work  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 . where is the fractal dimension for pores, is diameter, is the length scale, and is the total number of pores whose sizes equal to and greater than . The number of pores whose sizes range from to is when water flow through the pores of porous rock, the capillaries may be tortuous. These tortuous capillaries could be expressed by fractal equation . Consider where is the tortuosity fractal dimension and lies in the range , which represents the extent of convoluted ness of capillary pathways for fluid flow through a medium. Note that for a straight capillary path , and a higher value of corresponds to a highly tortuous capillary. Let the diameter of a capillary in the medium be and let its tortuous length along the flow direction be . is representative length of channels. With a straight capillary, = . The total volume of pores from to can be obtained from (1) as where is the whole volume of the pores, is given by (2).
Similarity to the state above, the volume of pores/particles from to can be obtained from (1) as The effective saturation of volume 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, is the gravity acceleration, 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 . 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 is the relative hydraulic conductivity. Inserting (8) to (10) we obtain the following form for 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 and residual water content .
The novel constitutive model has some similarities with the well-known Brooks-Corey model, which is and , where is the reciprocal of air entry pressure and is a model parameter related to pore size distribution.
When , , and , 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 , cm, and cm, cm, and 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.
(a) Effective saturation
(b) Relative hydraulic conductivity
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 . In order to analyze the influence of the fracture density we consider different values of . With (8) and (11), the relationship is obtained (Figure 3). With the physical constants cm, cm, and cm, the relationship between and can also be determined (Figure 4).
In order to test the proposed invasion depth model, it is crucial to correctly determine the fractal dimensions and .
According to the following formula proved recently by Yu and Li [22, 23], where is the Euclidean dimension and 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 is permeability, is Kozeny constant, which considers the tortuosity of capillaries and pore non-uniformity. And the parameters and are all proved . Consider where From (12) and (13), the relationship among , , and is determined (Figure 5).
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 which affected other model parameters and should not be neglected while it was often neglected in the past investigation is considered in this paper.
B. Berkowitz and A. Hadad, “Fractal and multifractal measures of natural and synthetic fracture networks,” Journal of Geophysical Research B, vol. 102, no. 6, pp. 12205–12218, 1997.View at: Google Scholar
R. Moussa, “Is the drainage network a fractal Sierpinskispace,” Water Resources Research, vol. 33, pp. 2399–2408, 1997.View at: Google Scholar
B. Q. Xiao, J. T. Fan, and F. Ding, “Prediction of relative permeability of unsaturated porous media based on fractal theory and monte carlo simulation,” Energy & Fuels, vol. 26, pp. 6971–6978, 2012.View at: Google Scholar
P. G. Obuko and K. Aki, “Fractal geometry in the San Andreas fault system,” Journal of Geophysical Research, vol. 92, pp. 345–355, 1987.View at: Google Scholar
C. A. Barton and M. D. Zoback, “Self-similar distribution and properties of macroscopic fractures at depth in crystalline rock in the Cajon Pass scientific drill hole,” Journal of Geophysical Research, vol. 97, no. 4, pp. 5181–5200, 1992.View at: Google Scholar
J. S. Wu and B. M. Yu, “A fractal resistance model for flow through porous media,” International Journal of Heat and Mass Transfer, pp. 3925–3932, 2007.View at: Google Scholar
R. H. Brooks and A. T. Corey, Hydraulic Properties of Porous Media. Hydrology Paper no. 3, lorado State University, Fort Collins, CO, 1964.
N. T. Burdine, “Relative permeability from pore size distribution data,” Transactions. American Institute of Mining, Metallurgical and Petroleum Engineers, vol. 198, pp. 71–78, 1953.View at: Google Scholar
T. M. Nordberg and S. T. Thorolfsson, “Low impact development and bioretention areas in cold climates,” in Proceedings of the World Water and Environmental Resources Congress: Critical Transitions in Water and Environmental Resources Management, pp. 3409–3418, Salt Lake City, Utah, USA, July 2004.View at: Google Scholar
A. J. Katz and A. H. Thompson, “Fractal sandstone pores: implications for conductivity and formation,” Physical Review Letters, vol. 54, pp. 1325–1328, 1985.View at: Google Scholar
W. F. Hunt and A. R. Jarrett, “Evaluating bioretention areas from two field sites in North Carolina,” in Proceedings of the World Water and Environmental Resources Congress: Critical Transitions in Water and Environmental Resources Management, pp. 797–806, Salt Lake City, Utah, USA, July 2004.View at: Google Scholar