A new geometrical model for non-linear fluid flow through rough fractures
TL;DR: In this paper, a polynomial expression, like Forchheimer law, was used to describe the dependence of pressure drop on flow rate for non-linear fluid flow through rough fractures.
About: This article is published in Journal of Hydrology.The article was published on 2010-07-28. It has received 118 citations till now. The article focuses on the topics: Laminar flow & Hele-Shaw flow.
Citations
More filters
11 Jun 2010
Abstract: The validity of the cubic law for laminar flow of fluids through open fractures consisting of parallel planar plates has been established by others over a wide range of conditions with apertures ranging down to a minimum of 0.2 µm. The law may be given in simplified form by Q/Δh = C(2b)3, where Q is the flow rate, Δh is the difference in hydraulic head, C is a constant that depends on the flow geometry and fluid properties, and 2b is the fracture aperture. The validity of this law for flow in a closed fracture where the surfaces are in contact and the aperture is being decreased under stress has been investigated at room temperature by using homogeneous samples of granite, basalt, and marble. Tension fractures were artificially induced, and the laboratory setup used radial as well as straight flow geometries. Apertures ranged from 250 down to 4µm, which was the minimum size that could be attained under a normal stress of 20 MPa. The cubic law was found to be valid whether the fracture surfaces were held open or were being closed under stress, and the results are not dependent on rock type. Permeability was uniquely defined by fracture aperture and was independent of the stress history used in these investigations. The effects of deviations from the ideal parallel plate concept only cause an apparent reduction in flow and may be incorporated into the cubic law by replacing C by C/ƒ. The factor ƒ varied from 1.04 to 1.65 in these investigations. The model of a fracture that is being closed under normal stress is visualized as being controlled by the strength of the asperities that are in contact. These contact areas are able to withstand significant stresses while maintaining space for fluids to continue to flow as the fracture aperture decreases. The controlling factor is the magnitude of the aperture, and since flow depends on (2b)3, a slight change in aperture evidently can easily dominate any other change in the geometry of the flow field. Thus one does not see any noticeable shift in the correlations of our experimental results in passing from a condition where the fracture surfaces were held open to one where the surfaces were being closed under stress.
1,557 citations
TL;DR: In this article, the Forchheimer equation coefficients for non-Darcy flow in deformable rough-walled fractures were evaluated using a triaxial cell under confining stresses varying from 1.0 MPa to 30 MPa.
214 citations
TL;DR: In this article, a quantitative criterion was developed to quantify the onset of nonlinear flow by comprehensive combination of Forchheimer's law and Reynolds number, and several high-precision water flow tests were carried out with different hydraulic gradients then the critical Reynolds number was determined based on the developed criterion.
Abstract: This paper experimentally investigates the role of shear processes on the variation of critical Reynolds number and nonlinear flow through rough-walled rock fractures. A quantitative criterion was developed to quantify the onset of nonlinear flow by comprehensive combination of Forchheimer's law and Reynolds number. At each shear displacement, several high-precision water flow tests were carried out with different hydraulic gradients then the critical Reynolds number was determined based on the developed criterion. The results show that (i) the Forchheimer's law was fitted very well to experimental results of nonlinear fluid flow through rough-walled fractures, (ii) the coefficients of viscous and inertial pressure drops experience 4 and 7 orders of magnitude reduction during shear displacement, respectively, and (iii) the critical Reynolds number varies from 0.001 to 25 and experiences 4 orders of magnitude enlargement by increasing shear displacement from 0 to 20 mm. These findings may prove useful in proper understanding of fluid flow through rock fractures, or inclusions in computational studies of large-scale nonlinear flow in fractured rocks.
199 citations
Cites background or methods from "A new geometrical model for non-lin..."
...The general description of Newtonian fluid motion in a rough-walled fracture is given by the Navier-Stokes (NS) equations composed of a set of coupled nonlinear partial derivatives of varying orders [Brush and Thomson, 2003; Javadi et al., 2010; Zimmerman and Bodvarsson, 1996)]....
[...]
...For instance, by increasing the fracture roughness, the critical Reynolds number will rapidly decrease [Javadi et al., 2010; Konzu and Kueper, 2004; Parrish, 1963; Ranjith and Darlington, 2007; Zimmerman et al., 2004]....
[...]
...…Ji et al., 2008; Konzu and Kueper, 2004; Nowamooz et al., 2009; Parrish, 1963; Qian et al., 2005, 2011; Quinn et al., 2011; Ranjith and Darlington, 2007; Ranjith and Viete, 2011; Zimmerman et al., 2004], and numerically [Bu es et al., 2004; Javadi et al., 2010; Kolditz, 2001; Skjetne et al., 1999]....
[...]
TL;DR: In this article, the authors investigated the nonlinear flow characteristics at low Reynolds number through rough-walled fractures subjected to a wide range of confining pressures (1.0-30.0 MPa).
166 citations
TL;DR: In this paper, a new data analysis method based on polynomial fitting was introduced to investigate the relationship between flow velocity and hydraulic gradient, and the results showed that the flow velocity versus hydraulic gradient data gradient shows a nonlinear relationship at very low hydraulic, possibly due to strong solid-water interaction, but becomes approximately linear after the hydraulic gradient is high enough.
155 citations
References
More filters
TL;DR: The validity of the cubic law for laminar flow of fluids through open fractures consisting of parallel planar plates has been established by others over a wide range of conditions with apertures ranging down to a minimum of 0.2 µm.
Abstract: The validity of the cubic law for laminar flow of fluids through open fractures consisting of parallel planar plates has been established by others over a wide range of conditions with apertures ranging down to a minimum of 0.2 µm. The law may be given in simplified form by Q/Δh = C(2b)3, where Q is the flow rate, Δh is the difference in hydraulic head, C is a constant that depends on the flow geometry and fluid properties, and 2b is the fracture aperture. The validity of this law for flow in a closed fracture where the surfaces are in contact and the aperture is being decreased under stress has been investigated at room temperature by using homogeneous samples of granite, basalt, and marble. Tension fractures were artificially induced, and the laboratory setup used radial as well as straight flow geometries. Apertures ranged from 250 down to 4µm, which was the minimum size that could be attained under a normal stress of 20 MPa. The cubic law was found to be valid whether the fracture surfaces were held open or were being closed under stress, and the results are not dependent on rock type. Permeability was uniquely defined by fracture aperture and was independent of the stress history used in these investigations. The effects of deviations from the ideal parallel plate concept only cause an apparent reduction in flow and may be incorporated into the cubic law by replacing C by C/ƒ. The factor ƒ varied from 1.04 to 1.65 in these investigations. The model of a fracture that is being closed under normal stress is visualized as being controlled by the strength of the asperities that are in contact. These contact areas are able to withstand significant stresses while maintaining space for fluids to continue to flow as the fracture aperture decreases. The controlling factor is the magnitude of the aperture, and since flow depends on (2b)3, a slight change in aperture evidently can easily dominate any other change in the geometry of the flow field. Thus one does not see any noticeable shift in the correlations of our experimental results in passing from a condition where the fracture surfaces were held open to one where the surfaces were being closed under stress.
1,729 citations
11 Jun 2010
Abstract: The validity of the cubic law for laminar flow of fluids through open fractures consisting of parallel planar plates has been established by others over a wide range of conditions with apertures ranging down to a minimum of 0.2 µm. The law may be given in simplified form by Q/Δh = C(2b)3, where Q is the flow rate, Δh is the difference in hydraulic head, C is a constant that depends on the flow geometry and fluid properties, and 2b is the fracture aperture. The validity of this law for flow in a closed fracture where the surfaces are in contact and the aperture is being decreased under stress has been investigated at room temperature by using homogeneous samples of granite, basalt, and marble. Tension fractures were artificially induced, and the laboratory setup used radial as well as straight flow geometries. Apertures ranged from 250 down to 4µm, which was the minimum size that could be attained under a normal stress of 20 MPa. The cubic law was found to be valid whether the fracture surfaces were held open or were being closed under stress, and the results are not dependent on rock type. Permeability was uniquely defined by fracture aperture and was independent of the stress history used in these investigations. The effects of deviations from the ideal parallel plate concept only cause an apparent reduction in flow and may be incorporated into the cubic law by replacing C by C/ƒ. The factor ƒ varied from 1.04 to 1.65 in these investigations. The model of a fracture that is being closed under normal stress is visualized as being controlled by the strength of the asperities that are in contact. These contact areas are able to withstand significant stresses while maintaining space for fluids to continue to flow as the fracture aperture decreases. The controlling factor is the magnitude of the aperture, and since flow depends on (2b)3, a slight change in aperture evidently can easily dominate any other change in the geometry of the flow field. Thus one does not see any noticeable shift in the correlations of our experimental results in passing from a condition where the fracture surfaces were held open to one where the surfaces were being closed under stress.
1,557 citations
TL;DR: In this article, the authors derived the cubic law of the Navier-Stokes equations for flow between smooth, parallel plates and showed that the effective hydraulic aperture is less than the mean aperture, by a factor that depends on the ratio of the mean value of the aperture to its standard deviation.
Abstract: The flow of a single-phase fluid through a rough-walled rock fracture is discussed within the context of fluid mechanics. The derivation of the ‘cubic law’ is given as the solution to the Navier-Stokes equations for flow between smooth, parallel plates - the only fracture geometry that is amenable to exact treatment. The various geometric and kinematic conditions that are necessary in order for the Navier-Stokes equations to be replaced by the more tractable lubrication or Hele-Shaw equations are studied and quantified. In general, this requires a sufficiently low flow rate, and some restrictions on the spatial rate of change of the aperture profile. Various analytical and numerical results are reviewed pertaining to the problem of relating the effective hydraulic aperture to the statistics of the aperture distribution. These studies all lead to the conclusion that the effective hydraulic aperture is less than the mean aperture, by a factor that depends on the ratio of the mean value of the aperture to its standard deviation. The tortuosity effect caused by regions where the rock walls are in contact with each other is studied using the Hele-Shaw equations, leading to a simple correction factor that depends on the area fraction occupied by the contact regions. Finally, the predicted hydraulic apertures are compared to measured values for eight data sets from the literature for which aperture and conductivity data were available on the same fracture. It is found that reasonably accurate predictions of hydraulic conductivity can be made based solely on the first two moments of the aperture distribution function, and the proportion of contact area.
1,003 citations
TL;DR: In this article, a simulation of flow between rough surfaces was done using a fractal model of surface topography and the hydraulic aperture was compared to the mean separation of the surfaces.
Abstract: Fluid flow through rock joints is commonly described by the parallel plate model where the volume flow rate varies as the cube of the joint aperture. However, deviations from this model are expected because real joint surfaces are rough and contact each other at discrete points. To examine this problem further, a computer simulation of flow between rough surfaces was done. Realistic rough surfaces were generated numerically using a fractal model of surface topography. Pairs of these surfaces were placed together to form a “joint” with a random aperture distribution. Reynolds equation, which describes laminar flow between slightly nonplanar and nonparallel surfaces, was solved on the two-dimensional aperture mesh by the finite-difference method. The solution is the local volume flow rate through the joint. This solution was used directly in the cubic law to get the so-called “hydraulic aperture.” For various surface roughnesses (fractal dimensions) the hydraulic aperture was compared to the mean separation of the surfaces. At large separations the surface topography has little effect. At small separations the flow is tortuous, tending to be channeled through high-aperture regions. The parameter most affecting fluid flow through rough joints is the ratio of the mean separation between the surfaces to the root-mean-square surface height. This parameter describes the distance the surface asperities protrude into the fluid and accounts for most of the disagreement with the parallel plate model. Variations in the fractal dimension produce only a second-order effect on the fluid flow. For the range of joint closures expected during elastic deformation these results show that the actual flow rate between rough surfaces is about 70–90% of that predicted by the parallel plate model.
826 citations
TL;DR: In this article, the effective pressure p e for measurements of fluid permeability is shown to be ( p c − sp p ) where p c is confining pressure, p p is pore pressure, and s depends on the topography of the fracture surfaces and rock type.
668 citations