Tectonic interpretation of the connectivity of a multiscale fracture system in limestone

This paper studies the statistics and tectonism of a multiscale natural fracture system in limestone. The fracture network exhibits a self‐similar characteristic with a correlation between its power law length exponent a and fractal dimension D, i.e., a ≈ D + 1. Contradicting the scale‐invariant connectivity of idealized self‐similar systems, the percolation state of trace patterns mapped at different scales and localities of the study area varies significantly, from well to poorly connected. A tectonic interpretation based on a polyphase fracture network evolution history is proposed to explain this discrepancy. We present data to suggest that the driving force for fracture formation may be dissipated at the end of a tectonic event when the system becomes connected. However, the “effective” connectivity can successively be reduced by cementation of early fractures and reestablished by subsequent cracking, rendering a variable “apparent” connectivity that can be significantly above the percolation threshold.ISSN:0094-8276ISSN:1944-800

A new approach to upscaling fracture network models while preserving geostatistical and geomechanical characteristics

A new approach to upscaling two-dimensional fracture network models is proposed for preserving geostatistical and geomechanical characteristics of a smaller-scale “source” fracture pattern. First, the scaling properties of an outcrop system are examined in terms of spatial organization, lengths, connectivity, and normal/shear displacements using fractal geometry and power law relations. The fracture pattern is observed to be nonfractal with the fractal dimension D ≈ 2, while its length distribution tends to follow a power law with the exponent 2 < a < 3. To introduce a realistic distribution of fracture aperture and shear displacement, a geomechanical model using the combined finite-discrete element method captures the response of a fractured rock sample with a domain size L = 2 m under in situ stresses. Next, a novel scheme accommodating discrete-time random walks in recursive self-referencing lattices is developed to nucleate and propagate fractures together with their stress- and scale-dependent attributes into larger domains of up to 54 m × 54 m. The advantages of this approach include preserving the nonplanarity of natural cracks, capturing the existence of long fractures, retaining the realism of variable apertures, and respecting the stress dependency of displacement-length correlations. Hydraulic behavior of multiscale growth realizations is modeled by single-phase flow simulation, where distinct permeability scaling trends are observed for different geomechanical scenarios. A transition zone is identified where flow structure shifts from extremely channeled to distributed as the network scale increases. The results of this paper have implications for upscaling network characteristics for reservoir simulation