3D geostatistical modeling of fracture system in a granitic massif to characterize hydraulic properties and fracture distribution

Abstract This study integrates 3D models of rock fractures from different sources and hydraulic properties aimed at identifying relationships between fractures and permeability. The Tono area in central Japan, chiefly overlain by Cretaceous granite, was examined because of the availability of a unique dataset from deep borehole data at 26 sites. A geostatistical method (GEOFRAC) that can incorporate orientations of sampled data was applied to 50,900 borehole fractures for spatial modeling of fractures over a 12 km by 8 km area, to a depth of 1.5 km. GEOFRAC produced a plausible 3D fracture model, in that the orientations of simulated fractures correspond to those of the sample data and the continuous fractures appeared near a known fault. Small-scale fracture distributions with dominant orientations were also characterized around the two shafts using fracture data from the shaft walls. By integrating the 3D model of hydraulic conductivity using sequential Gaussian simulation with the GEOFRAC fractures from the borehole data, the fracture sizes and directions that strongly affect permeable features were identified. Four fracture-related elements: lineaments from a shaded 10-m DEM, GEOFRAC fractures using the borehole and shaft data, and microcracks from SEM images, were used for correlating fracture attributes at different scales. The consistency of the semivariogram models of distribution densities was identified. Using an experimental relationship between hydraulic conductivity and fracture length, the fractures that typically affect the hydraulic properties at the drift scale were surmised to be in the range 100–200 m. These results are useful for a comprehensive understanding of rock fracture systems and their hydraulic characteristics at multiple scales in a target area.

[1]  Philipp Blum,et al.  Evaluation of the representative elementary volume (REV) of a fractured geothermal sandstone reservoir , 2010 .

[2]  T. Mitchell,et al.  Towards quantifying the matrix permeability of fault damage zones in low porosity rocks , 2012 .

[3]  三枝 博光,et al.  Final Report on the Surface-based Investigation Phase (Phase I) at the Mizunami Underground Research Laboratory Project , 2011 .

[4]  M. Riley,et al.  Estimating the hydraulic properties of the fracture network in a sandstone aquifer. , 2007, Journal of contaminant hydrology.

[5]  Chaoshui Xu,et al.  A new computer code for discrete fracture network modelling , 2010, Comput. Geosci..

[6]  Allen M. Shapiro,et al.  A Stochastic Model of a Fractured Rock Conditioned by Measured Information , 1984 .

[7]  R. D. Terzaghi.  Sources of Error in Joint Surveys , 1965 .

[8]  S. P. Neuman.  Multiscale relationships between fracture length, aperture, density and permeability , 2008 .

[9]  Chrystel Dezayes,et al.  3D model of fracture zones at Soultz-sous-Forêts based on geological data, image logs, induced microseismicity and vertical seismic profiles , 2010 .

[10]  James P. Evans,et al.  Fault zone architecture and permeability structure , 1996 .

[11]  K. Koike,et al.  Auto-detection and integration of tectonically significant lineaments from SRTM DEM and remotely-sensed geophysical data , 2011 .

[12]  Neil A. Chapman,et al.  Geological Disposal of Radioactive Wastes - Concept, Status and Trends Almacenamiento Geológico de Residuos Radiactivos - Concepto, Situación y , 2006 .

[13]  Maria José Jurado,et al.  Two-dimensional geostatistical modeling and prediction of the fracture system in the Albala Granitic Pluton, SW Iberian Massif, Spain , 2001 .

[14]  Herbert H. Einstein,et al.  Characterizing rock joint geometry with joint system models , 1988 .

[15]  M. Vinches,et al.  Application of geostatistical characteristics of rock mass fracture systems in 3D model generation , 2008 .

[16]  Tayfun Babadagli,et al.  Estimation of equivalent fracture network permeability using fractal and statistical network properties , 2012 .

[17]  Yuichi Ichikawa,et al.  Spatial correlation structures of fracture systems for deriving a scaling law and modeling fracture distributions , 2006, Comput. Geosci..

[18]  R. Holdsworth,et al.  The spatial heterogeneity of structures in high porosity sandstones: Variations and granularity effects in orientation data , 2009 .

[19]  A. Thoraval,et al.  Use of an Integrated Discrete Fracture Network Code for Stochastic Stability Analyses of Fractured Rock Masses , 2012, Rock Mechanics and Rock Engineering.

[20]  D. McIlroy,et al.  Bioturbation influence on reservoir quality: A case study from the Cretaceous Ben Nevis Formation, Jeanne d'Arc Basin, offshore Newfoundland, Canada , 2010 .

[21]  Ki-Bok Min,et al.  Determining the equivalent permeability tensor for fractured rock masses using a stochastic REV approach: Method and application to the field data from Sellafield, UK , 2004 .

[22]  Katsuaki Koike,et al.  Morphotectonics inferred from the analysis of topographic lineaments auto-detected from DEMs: Application and validation for the Sinai Peninsula, Egypt , 2011 .

[23]  A. G. Journel,et al.  Indicator principal component kriging , 1991 .

[24]  P. Smith.  Santa Fe, New Mexico , 1969 .

[25]  D. Sornette,et al.  Hierarchical geometry of faulting , 1996 .

[26]  Olivier Bour,et al.  A note on the angular correction applied to fracture intensity profiles along drill core , 2006 .

[27]  H. Koide,et al.  Geological Disposal of Carbon Dioxide and Radioactive Waste in the Geotectonically Active Country of Japan , 2011 .

[28]  Y. Tsang,et al.  Some anomalous features of flow and solute transport arising from fracture aperture variability , 1990 .

[29]  Noelle E. Odling,et al.  Scaling and connectivity of joint systems in sandstones from western Norway , 1997 .

[30]  S. Ingebritsen,et al.  Permeability of the continental crust: Implications of geothermal data and metamorphic systems , 1999 .

[31]  María Asunción Soriano,et al.  Characterizing lineaments from satellite images and field studies in the central Ebro basin (NE Spain) , 1998 .

[32]  Yasuhiro Yamada,et al.  Geological Investigations for Geological Model of Deep Underground Geoenvironment at the Mizunami Underground Research Laboratory (MIU) , 2013 .

[33]  Sevket Durucan,et al.  Carbon dioxide storage risk assessment: Analysis of caprock fracture network connectivity , 2011 .

[34]  L. Jing,et al.  Shear-Induced Flow Channels in a Single Rock Fracture and Their Effect on Solute Transport , 2011 .

[35]  Katsuaki Koike,et al.  Incorporation of fracture directions into 3D geostatistical methods for a rock fracture system , 2011, Environmental Earth Sciences.

[36]  Paul A. Witherspoon,et al.  A Model for Steady Fluid Flow in Random Three‐Dimensional Networks of Disc‐Shaped Fractures , 1985 .

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

[38]  Analysis of the Spatial Variation In Rock Mass Properties Through Geostatistics , 1980 .

[39]  Zoe K. Shipton,et al.  A review of recent developments concerning the structure, mechanics and fluid flow properties of fault zones , 2010 .

[40]  J. Chilès.  Fractal and geostatistical methods for modeling of a fracture network , 1988 .

[41]  Roger Thunvik,et al.  Predicting mass transport in discrete fracture networks with the aid of geometrical field data , 1986 .

[42]  K. Koike,et al.  Lineament analysis of satellite images using a segment tracing algorithm (STA) , 1995 .

[43]  Peter A. Dowd,et al.  A Comparison of Methods for the Stochastic Simulation of Rock Fractures , 2007 .

[44]  Philippe Renard,et al.  Three-dimensional geometric modeling of a faulted domain: the Soultz Horst example (Alsace, France) , 1994 .

[45]  Spatial Distribution, Scaling and Self-similar Behavior of Fracture Arrays in the Los Planes Fault, Baja California Sur, Mexico , 2005 .

[46]  T. Tsuruta,et al.  Three-dimensional fracture distribution in relation to local cooling in a granitic body: An example from the Toki granitic pluton, Central Japan , 2012 .

[47]  M. Oda,et al.  Numerical experiments on permeability tensor and its application to jointed granite at Stripa Mine, Sweden , 1987 .

[48]  A. Torabi,et al.  Scaling of fault attributes: A review , 2011 .

[49]  Richard Gloaguen,et al.  Derivation of groundwater flow-paths based on semi-automatic extraction of lineaments from remote sensing data , 2011 .

[50]  D. Pollard,et al.  Imaging 3-D fracture networks around boreholes , 2002 .

[51]  D. Marcotte,et al.  GEOSTATISTICAL ANALYSIS OF FRACTURES IN SHEAR ZONES IN THE CHIBOUGAMAU AREA : APPLICATIONS TO STRUCTURAL GEOLOGY , 1997 .

[52]  W. Dershowitz,et al.  Interpretation of fracture spacing and intensity , 1992 .

[53]  M. Zoback,et al.  Stress perturbations associated with active faults penetrated by boreholes: Possible evidence for near-complete stress drop and a new technique for stress magnitude measurement , 1994 .

[54]  Kunio Ota,et al.  Development of the Quality Management System for Borehole Investigations: Part 1—Quality Assurance and Quality Control Methodology for Hydraulic Packer Testing , 2010 .

[55]  Maria José Jurado,et al.  Architecture of fault zones determined from outcrop, cores, 3-D seismic tomography and geostatistical modeling: example from the Albalá Granitic Pluton, SW Iberian Variscan Massif , 2003 .

[56]  G. M. Laslett,et al.  Censoring and edge effects in areal and line transect sampling of rock joint traces , 1982 .

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

[58]  Günter Zimmermann,et al.  Scale Dependence of Hydraulic and Structural Parameters in the Crystalline Rock of the KTB , 2003, Pure and Applied Geophysics.

[59]  J. Chilès,et al.  Taking into account the complexity of natural fracture systems in reservoir single-phase flow modelling , 2002 .

[60]  Diana M. Allen,et al.  Estimating regional-scale fractured bedrock hydraulic conductivity using discrete fracture network (DFN) modeling , 2012, Hydrogeology Journal.

[61]  J. Almeida,et al.  Construction of high-resolution stochastic geological models and optimal upscaling to a simplified layer-type hydrogeological model , 2012 .

[62]  James P. Evans,et al.  Permeability of fault-related rocks, and implications for hydraulic structure of fault zones , 1997 .

[63]  J. Yamamoto.  Correcting the Smoothing Effect of Ordinary Kriging Estimates , 2005 .

[64]  D. Kurtzman,et al.  Correlating fracture trends and hydraulic head using semivariogram cloud analysis , 2005, Ground water.

[65]  Olivier Bour,et al.  A statistical scaling model for fracture network geometry, with validation on a multiscale mapping of a joint network (Hornelen Basin, Norway) , 2002 .

[66]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[67]  C. Renshaw,et al.  Effect of mechanical interactions on the scaling of fracture length and aperture , 1997, Nature.

[68]  P. White,et al.  Three-dimensional geological modelling and multivariate statistical analysis of water chemistry data to analyse and visualise aquifer structure and groundwater composition in the Wairau Plain, Marlborough District, New Zealand , 2012 .

[69]  Davide Elmo,et al.  An Integrated Numerical Modelling–Discrete Fracture Network Approach Applied to the Characterisation of Rock Mass Strength of Naturally Fractured Pillars , 2010 .

[70]  Jan M. Nordbotten,et al.  Application of simplified models to CO2 migration and immobilization in large-scale geological systems , 2012 .

[71]  Katsuaki Koike,et al.  Tectonic architecture through Landsat-7 ETM+/SRTM DEM-derived lineaments and relationship to the hydrogeologic setting in Siwa region, NW Egypt , 2006 .

[72]  Kevin Hestir,et al.  Three-dimensional statistical modelling of a fractured rock mass—an example from the Fanay-Augères mine , 1989 .

[73]  Emmanuel Gringarten,et al.  3-D geometric description of fractured reservoirs , 1996 .

[74]  K. Koike,et al.  Construction and analysis of interpreted fracture planes through combination of satellite-image derived lineaments and digital elevation model data , 1998 .

[75]  N. Odling,et al.  Scaling of fracture systems in geological media , 2001 .