A permeability model for naturally fractured carbonate reservoirs

Abstract Based on a detailed structural study performed on a reservoir surface analogue, a fracture permeability model for carbonate reservoirs is proposed. This involves four hierarchical systems, each one assumed to convey fluids exclusively to that of immediately higher order. New basic equations are then provided, simulating hydrodynamic reservoir response. The large scale, first-order structures consist of faults, to which high-permeability structures (damage zones) are associated. At the meter scale, stratabound joint systems, together with bedding joints bounding the mechanical layers, form a connected network transporting fluids to the fault system. At a lower scale, down to crystal size, non-stratabound joints constitute a pervasive and capillary fracture network which conveys fluids within the rock mass. The non-fractured host rock constitutes the lower permeability system (at the crystal scale). As no published works include an effective integration of detailed structural and numerical models, the present study aims at covering a significant gap in the literature, also providing appropriate theoretical basis for subsequent studies on fracture network modeling and reservoir simulation. Besides its application in the field of reservoir development and management, the illustrated model may also be useful in environmental studies involving ground fluids such as e.g. in the field of special/toxic/nuclear waste management.

[1]  Muhammad Sahimi,et al.  Fractal-wavelet neural-network approach to characterization and upscaling of fractured reservoirs , 2000 .

[2]  K. B. Krystinik,et al.  Selected Characteristics of Limestone and Dolomite Reservoirs in the United States , 1985 .

[3]  Rodolfo Camacho Velazquez,et al.  Decline-Curve Analysis of Fractured Reservoirs With Fractal Geometry , 2008 .

[4]  A. Soares Direct Sequential Simulation and Cosimulation , 2001 .

[5]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[6]  M. R. Gross The origin and spacing of cross joints: examples from the Monterey Formation, Santa Barbara Coastline, California , 1993 .

[7]  S. Mazzoli,et al.  Tectonic evolution of Pliocene-Pleistocene wedge-top basins of the southern Apennines: new constraints from magnetic fabric analysis , 2012 .

[8]  E. Tondi,et al.  The crustal fault structure responsible for the 1703 earthquake sequence of central Italy , 1998 .

[9]  Herlander Mata-Lima Reservoir characterization with iterative direct sequential co-simulation: Integrating fluid dynamic data into stochastic model , 2008 .

[10]  F. Dati,et al.  Modes and timing of fracture network development in poly-deformed carbonate reservoir analogues, Mt. Chianello, southern Italy , 2012 .

[11]  F. Delay,et al.  Numerical simulations of geological reservoirs: improving their conditioning through the use of entropy , 2000 .

[12]  Agust Gudmundsson Geometry, formation and development of tectonic fractures on the Reykjanes Peninsula, southwest Iceland , 1987 .

[13]  Abbas Seifi,et al.  An Improved Model for Geostatistical Simulation of Fracture Parameters and their Effect on Static and Dynamic Models , 2008 .

[14]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[15]  Michel Dekking,et al.  A Modern Introduction to Probability and Statistics: Understanding Why and How , 2007 .

[16]  M. Zoback,et al.  Self‐similar distribution and properties of macroscopic fractures at depth in crystalline rock in the Cajon Pass Scientific Drill Hole , 1992 .

[17]  V. Di Donato,et al.  The Pliocene‐Quaternary wedge‐top basins of southern Italy: an expression of propagating lateral slab tear beneath the Apennines , 2012 .

[18]  S. Mazzoli,et al.  Thin-skinned versus thick-skinned structural models for Apulian carbonate reservoirs: constraints from the Val d'Agri Fields, S Apennines, Italy , 2004 .

[19]  E. Tondi,et al.  Fault zone characteristics and scaling properties of the Val d’Agri Fault System (Southern Apennines, Italy) , 2000 .

[20]  Frederik Michel Dekking,et al.  A Modern Introduction to Probability and Statistics , 2005 .

[21]  S. Laubach Practical approaches to identifying sealed and open fractures , 2003 .

[22]  W. Narr Estimating Average Fracture Spacing in Subsurface Rock , 1996 .

[23]  D. Scrocca,et al.  Applying Thick-skinned Tectonic Models to the Apennine Thrust Belt of Italy—Limitations and Implications , 2004 .

[24]  Benoit B. Mandelbrot,et al.  Fractal Geometry of Nature , 1984 .

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

[26]  M. R. Gross,et al.  Strain accommodated by brittle failure in adjacent units of the Monterey Formation, U.S.A.: scale effects and evidence for uniform displacement boundary conditions , 1995 .

[27]  C. C. Barton Fractal Analysis of Scaling and Spatial Clustering of Fractures , 1995 .

[28]  M. Cacas,et al.  Distribution of joints: probabilistic modelling and case study near Cardiff (Wales, U.K.) , 1997 .

[29]  W. Narr,et al.  Joint spacing in sedimentary rocks , 1991 .

[30]  S. Santini,et al.  Applying the Multiple Inverse Method to the analysis of earthquake focal mechanism data: New insights into the active stress field of Italy and surrounding regions , 2012 .

[31]  W. Belfield Multifractal characteristics of natural fracture apertures , 1994 .

[32]  N. P. Christensen,et al.  Variations in fracture system geometry and their implications for fluid flow in fractures hydrocarbon reservoirs , 1999, Petroleum Geoscience.

[33]  Ridha Gharbi,et al.  Use of reservoir simulation for optimizing recovery performance , 2004 .

[34]  R. Schlische,et al.  The geometric and statistical evolution of normal fault systems: an experimental study of the effects of mechanical layer thickness on scaling laws , 2001 .

[35]  Zhixi Chen,et al.  Integrated conditional global optimisation for discrete fracture network modelling , 2006, Comput. Geosci..

[36]  J. J. Walsh,et al.  Influence of layering on vein systematics in line samples , 1999, Geological Society, London, Special Publications.

[37]  Stephen E. Laubach,et al.  A scale-independent approach to fracture intensity and average spacing measurement , 2006 .

[38]  G. Massonnat,et al.  Hierarchical Organization of Flow Network in Fractured Carbonate Reservoirs: Identification and Characterization of Key Parameters , 2002 .

[39]  M. Giorgioni,et al.  Quantifying uncertainties in multi-scale studies of fractured reservoir analogues: Implemented statistical analysis of scan line data from carbonate rocks , 2010 .

[40]  M. D'Errico,et al.  Tectonic burial and 'young' (< 10 Ma) exhumation in the southern Apennines fold and thrust belt (Italy) , 2008 .

[41]  C. Margottini,et al.  ‘Diffuse faulting’ in the Machu Picchu granitoid pluton, Eastern Cordillera, Peru , 2009 .

[42]  D. Pollard,et al.  Fracture spacing in layered rocks: a new explanation based on the stress transition , 2000 .

[44]  Ronald A. Nelson,et al.  Geologic Analysis of Naturally Fractured Reservoirs , 1986 .

[45]  Nibir Mandal,et al.  Evidence for a non-linear relationship between fracture spacing and layer thickness , 1994 .

[46]  Improved statistical multi-scale analysis of fractured reservoir analogues , 2011 .

[47]  J. Angelier,et al.  Fracture spacing and its relation to bed thickness , 1989, Geological Magazine.

[48]  A double-porosity model for a fractured aquifer with non-Darcian flow in fractures , 2008 .

[49]  S. Mazzoli,et al.  Late Dolomitization in Basin Limestones of the Southern Apennines Fold and Thrust Belt (Italy) , 2011 .

[50]  Lanru Jing,et al.  A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering , 2003 .

[51]  Stephen E. Laubach,et al.  Aperture-size scaling variations in a low-strain opening-mode fracture set, Cozzette Sandstone, Colorado , 2009 .

[52]  E. Tondi,et al.  Reconstruction of continental margin architecture deformed by the contraction of the Lagonegro Basin, southern Apennines, Italy , 2001, Journal of the Geological Society.

[53]  Tom Manzocchi,et al.  Scaling relationships of joint and vein arrays from The Burren, Co. Clare, Ireland , 2001 .

[54]  Vincenzo Guerriero,et al.  Power Law Distribution: Method of Multi-scale Inferential Statistics , 2012 .

[55]  D. Pollard,et al.  An experimental study of the relationship between joint spacing and layer thickness , 1995 .

[56]  C. Invernizzi,et al.  Fault properties and fluid flow patterns from Quaternary faults in the Apennines, Italy , 2001 .