Estimation of REV size and three-dimensional hydraulic conductivity tensor for a fractured rock mass through a single well packer test and discrete fracture fluid flow modeling

Abstract A new methodology is presented to determine the representative elementary volume (REV) size and three-dimensional (3-D) hydraulic conductivity tensor for a fractured rock mass. First, a 3-D stochastic fracture network model was built and validated for a gneissic rock mass based on the fracture data mapped from scanline surveys at the site. This validated fracture network model was combined with the fracture data observed on a borehole to generate a stochastic-deterministic fracture network system in a cubic block around each packer test conducted at a different depth region in the same borehole. Each packer test was simulated numerically applying a developed discrete fracture fluid flow model to estimate the influenced region or effective range for the packer test. A cubic block of size 18 m, with the packer test interval of length about 6.5 m located at the centre of this block, was found to be suitable to represent the influenced region. Using this block size, the average flow rate per unit hydraulic gradient (defined as the transmissivity multiplied by mean width of flow paths) field for fractures was calibrated at different depth regions around the borehole by numerically simulating the packer tests conducted at different depth regions. The average flow rate per unit hydraulic gradient of the fractures that intersect the borehole was considered to be quite different to the average flow rate per unit hydraulic gradient of the fractures that do not intersect the borehole. A relation was developed to quantify the ratio between these two parameters. By studying the directional hydraulic conductivity behaviour of different cubic block sizes having the validated stochastic fracture network and calibrated hydraulic parameters, a REV for the hydraulic behaviour of the rock mass was estimated to be a block size of 15 m. The hydraulic conductivity tensor in 3-D computed through regression analysis using the calculated directional hydraulic conductivity values in many directions was found to be significantly anisotropic. The principal directions of the hydraulic conductivity tensor were found to be agreeable with the existing fracture system of the site. Further, the geometric hydraulic conductivity calculated was found to be comparable to the hydraulic conductivity estimated through the radial flow assumption in continuum porous media.

[1]  Pinnaduwa Kulatilake,et al.  Relations between Fracture Tensor Parameters and Jointed Rock Hydraulics , 1999 .

[2]  Björn Dverstorp,et al.  Conditional simulations of fluid flow in three-dimensional networks of discrete fractures , 1987 .

[3]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[4]  Pinnaduwa Kulatilake,et al.  RELATION BETWEEN DISCONTINUITY SIZE AND TRACE LENGTH. , 1986 .

[5]  Pinnaduwa H.S.W. Kulatilake,et al.  Analysis of structural homogeneity of rock masses , 1990 .

[6]  K. Karasaki,et al.  Hydraulic well testing inversion for modeling fluid flow in fractured rocks using simulated annealing: a case study at Raymond field site, California , 2000 .

[7]  Pinnaduwa H.S.W. Kulatilake,et al.  EFFECT OF BLOCK SIZE AND JOINT GEOMETRY ON JOINTED ROCK HYDRAULICS AND REV , 2000 .

[8]  P. Kulatilake,et al.  Effect of joint geometry and transmissivity on jointed rock hydraulics , 1999 .

[9]  Franklin W. Schwartz,et al.  A Stochastic Analysis of Macroscopic Dispersion in Fractured Media , 1983 .

[10]  A. Rouleau,et al.  Stochastic discrete fracture simulation of groundwater flow into an underground excavation in granite , 1987 .

[11]  Paul A. Witherspoon,et al.  The relationship of the degree of interconnection to permeability in fracture networks , 1985 .

[12]  Pinnaduwa Kulatilake,et al.  Estimation of mean trace length of discontinuities , 1984 .

[13]  G. Marsily,et al.  Modeling fracture flow with a stochastic discrete fracture network: calibration and validation: 1. The flow model , 1990 .

[14]  K. Karasaki,et al.  A multidisciplinary fractured rock characterization study at Raymond field site, Raymond, CA. , 2000 .

[15]  Chin-Fu Tsang,et al.  Channel model of flow through fractured media , 1987 .

[16]  E. A. Sudicky,et al.  The Laplace Transform Galerkin Technique for large-scale simulation of mass transport in discretely fractured porous formations , 1992 .

[17]  C. Fidelibus,et al.  Derivation of equivalent pipe network analogues for three‐dimensional discrete fracture networks by the boundary element method , 1999 .

[18]  S. P. Neuman,et al.  Field Determination of the Three-Dimensional Hydraulic Conductivity Tensor of Anisotropic Media: 2. Methodology and Application to Fractured Rocks , 1985 .

[19]  Pinnaduwa Kulatilake,et al.  Discontinuity geometry characterization in a tunnel close to the proposed permanent shiplock area of the three gorges dam site in China , 1996 .

[20]  Jim Gallanes,et al.  On Target: The Arrowhead East and West Tunnels , 1996 .

[21]  Peter Clive. Robinson,et al.  Connectivity, flow and transport in network models of fractured media , 1984 .

[22]  G. W. Wathugala,et al.  A general procedure to correct sampling bias on joint orientation using a vector approach , 1990 .

[23]  Ove Stephansson,et al.  Joint network modelling with a validation exercise in Stripa mine, Sweden , 1993 .

[24]  Ove Stephansson,et al.  Effect of finite size joints on the deformability of jointed rock in three dimensions , 1993 .

[25]  Pinnaduwa Kulatilake,et al.  Groundwater resources evaluation case study via discrete fracture flow modeling , 2001 .

[26]  P. Kulatilake,et al.  Box fractal dimension as a measure of statistical homogeneity of jointed rock masses , 1997 .

[27]  S. P. Neuman,et al.  Field Determination of the Three‐Dimensional Hydraulic Conductivity Tensor of Anisotropic Media: 1. Theory , 1985 .