Fluid flow in synthetic rough fractures and application to the Hachimantai geothermal hot dry rock test site

Fracture profiles from the Hachimantai geothermal hot dry rock (HDR) test site in northern Japan have been measured and analyzed to characterise their geometrical properties. These properties have been used to create a population of numerical synthetic fractures that were tuned to imitate all the geometric and statistical properties of the natural fracture (described in a companion paper). Such fractures have been used as input boundary conditions in three types of modeling. (1) Simple elastic normal closure relating the aperture of the fracture to applied normal load. (2) Hagen-Poiseuille calculations of fluid transmissivity in the fracture as a function of normal load, fracture fluid pressure, and temperature. (3) Two-dimensional flow modeling within the rough walled fracture using Reynolds equation. The modeled closure of these fractures provided a realistic relationship between normal and fracture fluid pressure and aperture. When this pressure/aperture relationship was combined with a Hagen-Poiseuille approach to calculating fluid transmissivity in the fracture, we obtained results which, when compared with data from field transmissivity tests at the Hachimantai site, showed that the functional dependence of fluid transmissivity with fracture fluid pressure was well modeled, but overestimated by a factor of about 2. Reynolds equation flow modeling was carried out in the synthetic fracture to ascertain the extent to which the Hagen-Poiseuille law was overestimating the transmissivity due to the rough surfaces affecting the fluid flow in the fracture. When the calculations were corrected for this effect, the overestimation in the fluid transmissivity was reduced considerably.

[1]  P. Kitanidis,et al.  Prediction of single phase transport parameters in a variable aperture fracture , 1995 .

[2]  J. R. Macmillan,et al.  Stochastic analysis of spatial variability in subsurface flows: 2. Evaluation and application , 1978 .

[3]  D. Saupe Algorithms for random fractals , 1988 .

[4]  Y. Guéguen,et al.  Pressure dependence of permeability: a model for cracked rocks , 1989 .

[5]  Stephen R. Brown,et al.  Fluid flow through rock joints: The effect of surface roughness , 1987 .

[6]  Stephen R. Brown,et al.  Applicability of the Reynolds Equation for modeling fluid flow between rough surfaces , 1995 .

[7]  Stephen R. Brown,et al.  Closure of random elastic surfaces in contact , 1985 .

[8]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

[9]  David D. Nolte,et al.  Fluid percolation through single fractures , 1988 .

[10]  J. Archard Elastic deformation and the laws of friction , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[11]  Stephen R. Brown,et al.  Transport of fluid and electric current through a single fracture , 1989 .

[12]  J. Greenwood,et al.  Contact of nominally flat surfaces , 1966, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[13]  Stephen R. Brown Simple mathematical model of a rough fracture , 1995 .

[14]  Vincent C. Tidwell,et al.  Laboratory investigation of matrix imbibition from a flowing fracture , 1995 .

[15]  T. Illangasekare,et al.  A mathematical model for flow and solute transport in non-homogeneous rock fractures , 1994 .

[16]  K. Matsuki,et al.  Fluid flow in fractally rough synthetic fractures , 1997 .

[17]  Paul W. J. Glover,et al.  Synthetic rough fractures in rocks , 1998 .