Effects of local rock heterogeneities on the hydromechanics of fractured rocks using a digital-image-based technique

A digital-image-based (DIB) finite element approach is developed based on the numerical code rock failure process analysis (RFPA) to characterize micro-scale rock heterogeneity, and to understand the impact of micro-scale rock heterogeneity on the macro-scale hydromechanical response of rocks. The DIB technique incorporates small-scale spatial variability of initial deformation modulus, strength and permeability directly into a coupled hydromechanical model. Variability in Young's modulus, strength, and permeability is applied by a property map defined from the pixel-scale of a digital image. In the RFPA, mechanical deformation is followed, including the accumulation of damage applied in individual elements, which modifies modulus, strength, and permeability with the intensity of damage. The RFPA simulates progressive failure in fractured rocks, representing both the growth of existing fractures and the formation of new fractures, without having to identify crack tips and their interaction explicitly. In this DIB simulation approach, image voxels are used to give equivalent mechanical and flow properties. These property maps are ported to the model capable of solving directly for the evolving deformation, and fluid flow fields. The model is validated through comparisons of the simulated results with phenomenological observations documented in previous studies. The validated model is then applied to investigate the hydromechanical response of fractured rock characterized by digital image. The model is able to reproduce the spatial evolution of damage in the sample, the coalescence of existing cracks, and the formation of new cracks.

[1]  Leslie George Tham,et al.  Digital image-based numerical modeling method for prediction of inhomogeneous rock failure , 2004 .

[2]  Leslie George Tham,et al.  Numerical studies of the influence of microstructure on rock failure in uniaxial compression — Part I: effect of heterogeneity , 2000 .

[3]  N. Kikuchi,et al.  Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue , 1994, Biotechnology and bioengineering.

[4]  L. Tham,et al.  Influence of Heterogeneity of Mechanical Properties on Hydraulic Fracturing in Permeable Rocks , 2004 .

[5]  Z. T. Bieniawski,et al.  Brittle fracture propagation in rock under compression , 1965 .

[6]  R. Zimmerman,et al.  EFFECTIVE STRESS LAW FOR THE PERMEABILITY OF CLAY-RICH SANDSTONES , 2004 .

[7]  C. A. Tang,et al.  Coupled analysis of flow, stress and damage (FSD) in rock failure , 2002 .

[8]  Chun’an Tang,et al.  Micromechanical Model for Simulating the Fracture Process of Rock , 2004 .

[9]  E. Z. Lajtai Brittle fracture in compression , 1974 .

[10]  T. R. Reid,et al.  A semi-automated methodology for discontinuity trace detection in digital images of rock mass exposures , 2000 .

[11]  Zhong Qi Yue,et al.  Digital image processing for aggregate orientation in asphalt concrete mixtures , 1997 .

[12]  C. Tang,et al.  Numerical simulation of progressive rock failure and associated seismicity , 1997 .

[13]  W. F. Brace,et al.  A note on brittle crack growth in compression , 1963 .

[14]  Teng-fong Wong,et al.  The transition from brittle faulting to cataclastic flow: Permeability evolution , 1997 .

[15]  J. Chermant,et al.  Some fields of applications of automatic image analysis in civil engineering , 2001 .

[16]  Leslie George Tham,et al.  Finite element modeling of geomaterials using digital image processing , 2003 .

[17]  W. F. Brace A note on permeability changes in geologic material due to stress , 1978 .

[18]  Chin-Fu Tsang,et al.  Linking Thermal, Hydrological and Mechanical Processes in Fractured Rocks , 1999 .

[19]  Leslie George Tham,et al.  A replica method for observing microcracks on rock surfaces , 2003 .

[20]  S. C. Blair,et al.  Analysis of compressive fracture in rock using statistical techniques: Part I. A non-linear rule-based model , 1998 .

[21]  M. Zoback,et al.  The effect of microcrack dilatancy on the permeability of westerly granite , 1975 .

[22]  P. Philippi,et al.  3D reconstitution of porous media from image processing data using a multiscale percolation system , 2004 .

[23]  Jonny Rutqvist,et al.  The role of hydromechanical coupling in fractured rock engineering , 2003 .

[24]  J. P. Harrison,et al.  Development of a hydro-mechanical local degradation approach and its application to modelling fluid flow during progressive fracturing of heterogeneous rocks , 2005 .

[25]  M. Zako,et al.  Microstructure-based stress analysis and evaluation for porous ceramics by homogenization method with digital image-based modeling , 2003 .

[26]  Jianwen Liang,et al.  Digital image processing for non-linear system identification , 2004 .

[27]  J. Mier Fracture Processes of Concrete , 1997 .

[28]  Xavier Maldague,et al.  An Evaluation of Image Analysis Algorithms for Constructing Discontinuity Trace Maps , 2003 .

[29]  S. P. Neuman,et al.  Trends, prospects and challenges in quantifying flow and transport through fractured rocks , 2005 .

[30]  Teng-fong Wong,et al.  The transition from brittle faulting to cataclastic flow in porous sandstones: Mechanical deformation , 1997 .

[31]  R. Kranz,et al.  Crack-crack and crack-pore interactions in stressed granite , 1979 .

[32]  G. Zyvoloski,et al.  A numerical model for thermo-hydro-mechanical coupling in fractured rock , 1997 .

[33]  Jean-Louis Chermant,et al.  Why automatic image analysis? An introduction to this issue , 2001 .

[34]  Paulo Cesar Philippi,et al.  Using three-dimensional reconstructed microstructures for predicting intrinsic permeability of reservoir rocks based on a Boolean lattice gas method , 2002 .