Constructing a discrete fracture network constrained by seismic inversion data

Rock fractures are of great practical importance to petroleum reservoir engineering because they provide pathways for fluid flow, especially in reservoirs with low matrix permeability, where they constitute the primary flow conduits. Understanding the spatial distribution of natural fracture networks is thus key to optimising production. The impact of fracture systems on fluid flow patterns can be predicted using discrete fracture network models, which allow not only the 6 independent components of the second‐rank permeability tensor to be estimated, but also the 21 independent components of the fully anisotropic fourth‐rank elastic stiffness tensor, from which the elastic and seismic properties of the fractured rock medium can be predicted. As they are stochastically generated, discrete fracture network realisations are inherently non‐unique. It is thus important to constrain their construction, so as to reduce their range of variability and, hence, the uncertainty of fractured rock properties derived from them. This paper presents the underlying theory and implementation of a method for constructing a geologically realistic discrete fracture network, constrained by seismic amplitude variation with offset and azimuth data. Several different formulations are described, depending on the type of seismic data and prior geologic information available, and the relative strengths and weaknesses of each approach are compared. Potential applications of the method are numerous, including the prediction of fluid flow, elastic and seismic properties of fractured reservoirs, model‐based inversion of seismic amplitude variation with offset and azimuth data, and the optimal placement and orientation of infill wells to maximise production.     

FracKfinder: A MATLAB Toolbox for Computing Three-Dimensional Hydraulic Conductivity Tensors for Fractured Porous Media

Fractures in porous media have been documented extensively. However, they are often omitted from groundwater flow and mass transport models due to a lack of data on fracture hydraulic properties and the computational burden of simulating fractures explicitly in large model domains. We present a MATLAB toolbox, FracKfinder, that automates HydroGeoSphere (HGS), a variably saturated, control volume finite-element model, to simulate an ensemble of discrete fracture network (DFN) flow experiments on a single cubic model mesh containing a stochastically generated fracture network. Because DFN simulations in HGS can simulate flow in both a porous media and a fracture domain, this toolbox computes tensors for both the matrix and fractures of a porous medium. Each model in the ensemble represents a different orientation of the hydraulic gradient, thus minimizing the likelihood that a single hydraulic gradient orientation will dominate the tensor computation. Linear regression on matrices containing the computed three-dimensional hydraulic conductivity (K) values from each rotation of the hydraulic gradient is used to compute the K tensors. This approach shows that the hydraulic behavior of fracture networks can be simulated where fracture hydraulic data are limited. Simulation of a bromide tracer experiment using K tensors computed with FracKfinder in HGS demonstrates good agreement with a previous large-column, laboratory study. The toolbox provides a potential pathway to upscale groundwater flow and mass transport processes in fractured media to larger scales.     

Impact of fracture network geometry on free convective flow patterns

The effect of fracture network geometry on free convection in fractured rock is studied using numerical simulations. We examine the structural properties of fracture networks that control the onset and strength of free convection and the patterns of density-dependent flow. Applicability of the equivalent porous medium approach (EPM) is also tested, and recommendations are given, for which situations the EPM approach is valid. To date, the structural properties of fracture networks that determine free convective flow are examined only in few, predominantly simplified regular fracture networks. We consider fracture networks containing continuous, discontinuous, orthogonal and/or inclined discrete fractures embedded in a low-permeability rock matrix. The results indicate that bulk permeability is not adequate to infer the occurrence and magnitude of free convection in fractured rock. Fracture networks can inhibit or promote convection depending on the fracture network geometry. Continuous fracture circuits are the crucial geometrical feature of fracture networks, because large continuous fracture circuits with a large vertical extent promote convection. The likelihood of continuous fracture circuits and thus of free convection increases with increasing fracture density and fracture length, but individual fracture locations may result in great deviances in strength of convection between statistically equivalent fracture networks such that prediction remains subject to large uncertainty.     

Effect of seismic waves on the hydro-mechanical properties of fractured rock masses

The transmission of seismic waves in a particular region may influence the hydraulic properties of a rock mass, including permeability, which is one of the most important. To determine the effect of a seismic wave on the hydraulic behavior of a fractured rock mass, systematic numerical modeling was conducted. A number of discrete fracture network(DFN) models with a size of 20 m × 20 m were used as geometrical bases, and a discrete element method(DEM) was employed as a numerical simulation tool. Three different boundary conditions without(Type Ⅰ) and with static(Type Ⅱ) and dynamic(Type Ⅲ) loading were performed on the models, and then their permeability was calculated. The results showed that permeability in Type Ⅲ models was respectively 62.7% and 44.2% higher than in Type I and Type Ⅱ models. This study indicates that seismic waves can affect deep earth, and, according to the results, seismic waves increase the permeability and change the flow rate patterns in a fractured rock mass.     

Effect of NAPL exposure on the wettability and two‐phase flow in a single rock fracture

Surface‐wetting properties are an important cause of changing the groundwater and two‐phase fluid flows. Various factors affecting the surface wettability were investigated in a parallel‐walled glass fracture with non‐aqueous phase liquid (NAPL) (gasoline, diesel, trichloroethylene, and creosote) wetted surfaces. First, the effect of the duration of NAPL exposure on wettability change was considered at pre‐wet fracture surfaces using the various NAPL species, and the result showed that the surface became hydrophobic after the exposure time of NAPL exceeded 2000 min. Second, the initial wetting state of the surface affected the timing when the wettability change begins as well as the extent of the wettability change in an NAPL‐wetted rock fractures. Under the dry condition, the wettability change was completed within a very short time of exposure to NAPL (~5 min), and then it finally reached the intermediate and weakly NAPL wetting (contact angle of 118°). Under the pre‐wet condition, a relatively long time of exposure (~5000 min) was needed to observe the obvious change of the surface wettability, which was changed up to strongly NAPL wetting (contact angle of 142°). Third, the wettability changed by NAPL exposure was stable and maintained for a long time, regardless of water flushing rate and temperature. Finally, the wettability change by the exposure of NAPL on parallel fracture surfaces was evaluated at various groundwater flow velocities. Result showed that groundwater flow velocity has an important impact upon measured contact angle. Although fracture surfaces were exposed to NAPL at the low groundwater flow velocity, the wettability was not changed from hydrophilic to hydrophobic when the contact time between NAPL and mineral surfaces was not sufficient owing to the pulse‐type movement of NAPL. This implies that the variation of exposure pattern due to groundwater flow on the wettability change can be an important factor affecting the wettability change of fracture surface and migration behaviour at natural fractured rock aquifers in case of NAPL spill. Copyright © 2015 John Wiley & Sons, Ltd.     

Network properties of a two-dimensional natural fracture pattern

The influence that fractures exert on the permeability of a fractured rock is, to a large extent, controlled by the nature of the network formed by the fracture system. Here, the network properties of a two-dimensional natural pattern, mapped from the surface of a sandstone layer, are investigated and compared to those of realizations of spatially randomly distributed line segments with similar orientation and length distributions and line segment density (line length per unit area) to the natural pattern. These patterns are composed of clusters of varying size and shape, made up of interconnected fracture traces or line segments. Comparing the natural pattern with the realizations, the natural pattern was found to contain roughly half the number of clusters while the mass (total line length) of the largest cluster is approximately double that of the realizations. The size of the largest cluster controls the connectivity of the patterns, as can be seen by comparing the largest cluster of the natural pattern, which connects all four sides of the region, with those of the realizations, which are unconnected or connect only two sides. Cluster scaling characteristics were found to be similar in the natural pattern and the realizations and show a crossover from a dimension of one (their topological dimension) to two (the dimension of the embedding medium) at a point that corresponds to the fracture spacing. An investigation of the self-similarity dimension, using the box-counting method, showed similar characteristics with a broad transition zone between one- and two-dimensional behaviour at smaller box sizes. The patterns are therefore found to be non-fractal. The effect of the spatial distribution shown by the natural pattern is thus to modify the manner in which fractures are distributed among clusters, increasing connectivity (and permeability in the case of open fractures), but does not affect the cluster scaling characteristics or the self-similarity dimension of the fracture patterns.     

Development of a new semi-analytical model for cross-borehole flow experiments in fractured media

Analysis of borehole flow logs is a valuable technique for identifying the presence of fractures in the subsurface and estimating properties such as fracture connectivity, transmissivity and storativity. However, such estimation requires the development of analytical and/or numerical modeling tools that are well adapted to the complexity of the problem. In this paper, we present a new semi-analytical formulation for cross-borehole flow in fractured media that links transient vertical-flow velocities measured in one or a series of observation wells during hydraulic forcing to the transmissivity and storativity of the fractures intersected by these wells. In comparison with existing models, our approach presents major improvements in terms of computational expense and potential adaptation to a variety of fracture and experimental configurations. After derivation of the formulation, we demonstrate its application in the context of sensitivity analysis for a relatively simple two-fracture synthetic problem, as well as for field-data analysis to investigate fracture connectivity and estimate fracture hydraulic properties. These applications provide important insights regarding (i) the strong sensitivity of fracture property estimates to the overall connectivity of the system; and (ii) the non-uniqueness of the corresponding inverse problem for realistic fracture configurations.     

Induced Temperature Gradients to Examine Groundwater Flowpaths in Open Boreholes

Techniques for characterizing the hydraulic properties and groundwater flow processes of aquifers are essential to design hydrogeologic conceptual models. In this study, rapid time series temperature profiles within open‐groundwater wells in fractured rock were measured using fiber optic distributed temperature sensing (FO‐DTS). To identify zones of active groundwater flow, two continuous electrical heating cables were installed alongside a FO‐DTS cable to heat the column of water within the well and to create a temperature difference between the ambient temperature of the groundwater in the aquifer and that within the well. Additional tests were performed to examine the effects of pumping on hydraulic fracture interconnectivity around the well and to identify zones of increased groundwater flow. High‐ and low‐resolution FO‐DTS cable configurations were examined to test the sensitivities of the technique and compared with downhole video footage and geophysical logging to confirm the zones of active groundwater flow. Two examples are presented to demonstrate the usefulness of this new technique for rapid characterization of fracture zones in open boreholes. The combination of the FO‐DTS and heating cable has excellent scope as a rapid appraisal tool for borehole construction design and improving hydrogeologic conceptual models.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

A constitutive model for water flow in unsaturated fractured rocks

A conceptual model for describing effective saturation in fractured hard rock is presented. The fracture network and the rock matrix are considered as an equivalent continuum medium where each fracture is conceptualized as a porous medium of granular structure and the rock matrix is assumed to be impermeable. The proposed model is based on the representation of a rough‐walled fracture by an equivalent porous medium, which is described using classical constitutive models. A simple closed‐form equation for the effective saturation is obtained when the van Genuchten model is used to describe saturation inside fractures and fractal laws are assumed for both aperture and number of fractures. The relative hydraulic conductivity for the fractured rock is predicted from a simple relation derived by Liu and Bodvarsson. The proposed constitutive model contains three independent parameters, which may be obtained by fitting the proposed effective saturation curve to experimental data. Two of the model parameters have physical meaning and can be identified with the reciprocal of the air entry pressure values in the fractures of minimum and maximum apertures. Effective saturation and relative hydraulic conductivity curves match fairly well the simulated constitutive relations obtained by Liu and Bodvarsson. Copyright © 2008 John Wiley & Sons, Ltd.     

Are open fractures necessarily aligned with maximum horizontal stress?

Fluid flow in fractured rock is an increasingly central issue in recovering water and hydrocarbon supplies and geothermal energy, in predicting flow of pollutants underground, in engineering structures, and in understanding large-scale crustal behaviour. Conventional wisdom assumes that fluids prefer to flow along fractures oriented parallel or nearly parallel to modern-day maximum horizontal compressive stress, or S_Hmax. The reasoning is that these fractures have the lowest normal stresses across them and therefore provide the least resistance to flow. For example, this view governs how geophysicists design and interpret seismic experiments to probe fracture fluid pathways in the deep subsurface. Contrary to these widely held views, here we use core, stress measurement, and fluid flow data to show that S_Hmax does not necessarily coincide with the direction of open natural fractures in the subsurface (>3 km depth). Consequently, in situ stress direction cannot be considered to predict or control the direction of maximum permeability in rock. Where effective stress is compressive and fractures are expected to be closed, chemical alteration dictates location of open conduits, either preserving or destroying fracture flow pathways no matter their orientation.     

Characterization of a hydraulically induced crystalline‐rock fracture

Hydraulic fracturing has become an important technique for enhancing the permeability of hydrocarbon source rocks and increasing aquifer transmissivity in many hard rock environments where natural fractures are limited, yet little is known about the nature or behaviour of these hydraulically induced fractures as conduits to flow and transport. We propose that these fractures tend to be smooth based on observed hydraulic and transport behaviour. In this investigation a multi‐faceted approach was used to quantify the properties and characteristics of an isolated hydraulically induced fracture in crystalline rocks. Packers were used to isolate the fracture that is penetrated by two separate observation wells located approximately 33 m apart. A series of aquifer tests and an induced gradient tracer test were performed to better understand the nature of this fracture. Aquifer test results indicate that full recovery is slow because of the overall low permeability of the crystalline rocks. Drawdown tests indicate that the fracture has a transmissivity of 1–2 m²/day and a specific storage on the order of 2–9 × 10^?7/m. Analysis of a potassium–bromide tracer test break through curve shows classic Fickian behaviour with minimal tailing analogous to parallel plate flow. Virtually all of the tracer was recovered, and the breakthrough curve dilution indicates that the swept area is only about 11% of a radial flow field and the estimated aperture is ≤0.5 mm, which implies a narrow linear flow region. These outcomes suggest that transport within these hydraulically induced ‘smooth’ fractures in crystalline rocks is rapid with minimal mixing, small local velocity fluctuations and no apparent diffusion into the host rock or secondary fractures. Copyright © 2016 John Wiley & Sons, Ltd.     

Dependency of flow and transport properties on aperture distributions and compression states

Fluid conductivity and elastic properties in fractures depend on the aperture geometry – in particular, the roughness of fracture surfaces. In this study, we have characterized the surface roughness with a log-normal distribution and investigated the transport and flow behaviour of the fractures with varying roughness characteristics. Numerical flow and transport simulations have been performed on a single two-dimensional fracture surface, whose aperture geometry changes with different variances and correlation lengths in each realization. We have found that conventional measurement of hydraulic conductivity alone is insufficient to determine these two parameters. Transient transport measurements, such as the particle breakthrough time, provide additional constraints to the aperture distribution. Nonetheless, a unique solution to the fracture aperture distribution is still under-determined with both hydraulic conductivity and transport measurements. From numerical simulations at different compression states, we have found that the flow and transport measurements exhibit different rates of changes with respect to changes in compression. Therefore, the fracture aperture distribution could be further constrained by considering the flow and transport properties under various compression states.     

Influence of source region properties on scaling relations forM<0 events

Excavation induced seismic events with moment magnitudesM<0 are examined in an attempt to determine the role geology, excavation geometry, and stress have on scaling relations. Correlations are established based on accurate measurements of excavation geometry and methodology, stress regime, rock mass structure, local tectonics, and seismic locations. Scaling relations incorporated seismic moments and source radii obtained by spectral analysis, accounting for source, propagation, and site effects, and using Madariaga's dynamic circular fault model. Observations suggest that the interaction of stresses with pre-existing fractures, fracture complexity and depth of events are the main factors influencing source characteristics and scaling behaviour. Self-similar relationships were found for events at similar depths or for weakly structured rock masses with reduced clamping stresses, whereas a non-similar behaviour was found for events with increasing depth or for heavily fractured zones under stress confinement. Additionally, the scaling behaviour for combined data sets tended to mask the non-similar trends. Overall, depth and fracture complexity, initially thought as second order effects, appear to significantly influence source characteristics of seismic events withM<0 and consequently favour a non-similar earthquake generation process.     

A new upscaling method for fractured porous media

We present a method to determine equivalent permeability of fractured porous media. Inspired by the previous flow-based upscaling methods, we use a multi-boundary integration approach to compute flow rates within fractures. We apply a recently developed multi-point flux approximation Finite Volume method for discrete fracture model simulation. The method is verified by upscaling an arbitrarily oriented fracture which is crossing a Cartesian grid. We demonstrate the method by applying it to a long fracture, a fracture network and the fracture network with different matrix permeabilities. The equivalent permeability tensors of a long fracture crossing Cartesian grids are symmetric, and have identical values. The application to the fracture network case with increasing matrix permeabilities shows that the matrix permeability influences more the diagonal terms of the equivalent permeability tensor than the off-diagonal terms, but the off-diagonal terms remain important to correctly assess the flow field.     

Fault zone fabrics and geofluid properties as indicators of rock deformation modes

In this paper we present the results of a geostructural study on active faults in central Italy, where seismogenic fault zones occur as part of a Quaternary network dissecting and/or inverting earlier tectonic features of the central Apennines fold and thrust belt. In our work we focus on the possibility of using structurally-oriented quantitative analysis of fault fabrics and fluid inclusion studies for assessing the hydraulic properties and scaling relations of fault zones in order to evaluate the role and effects of the interaction between rock and fluids in the brittle deformation of strained crustal rock volumes. The results of our study show that this approach is appropriate for (i) assessing the structural permeability of faulted and fractured rock volumes, (ii) defining the conduit/barrier behaviour of fault zones to fluid flow, (iii) mapping spatial variations of the fluid pressure across different fault segments, (iv) evaluating the maturity of a structural network and the degree of interaction of linked structural discontinuities, (v) assessing fluid composition and the conditions of deformation by means of microstructural and fluid inclusion data.     

Numerical evaluation of the effect of fracture network connectivity in naturally fractured shale based on FSD model

In this paper, the effect of pre-existing discrete fracture network (DFN) connectivity on hydraulic fracturing is numerically investigated in a rock mass subjected to in-situ stress. The simulation results show that DFN connectivity has a significant influence on the hydraulic fracture (HF) & DFN interaction and hydraulic fracturing effectiveness, which can be characterized by the total interaction area, stimulated DFN length, stimulated HF length, leak-off ratio, and stimulated total length. In addition, even at the same fluid injection rate, simulation models exhibit different responses that are strongly affected by the DFN connectivity. At a low injection rate, total interaction area decreases with increasing DFN connectivity; at a high injection rate, total interaction area increases with the increase of DFN connectivity. However, for any injection rate, the stimulated DFN length increases and stimulated HF length decreases with the increase of connectivity. Generally, this work shows that the DFN connectivity plays a crucial role in the interaction between hydraulic fractures, the pre-existing natural fractures and hydraulic fracturing effectiveness; in return, these three factors affect treating pressure, created microseismicity and corresponding stimulated volume. This work strongly relates to the production technology and the evaluation of hydraulic fracturing effectiveness. It is helpful for the optimization of hydraulic fracturing simulations in naturally fractured formations.     

Numerical evaluation of the shear stimulation effect in naturally fractured formations

In shale gas fracking, the stimulated natural fracture system is often critical to the gas production. In this paper, we present the results of state-of-the-art modeling of a detailed parametric evolution of the shear stimulation effect in discrete fracture network (DFN) formations. Two-dimensional computational modeling studies have been used in an attempt towards understanding how naturally fractured reservoirs response in hydraulic fracturing. Simulations were conducted as a function of: (1) the in-situ stress ratio; (2) internal friction angle of DFN; (3) DFN orientation with the stress field; and (4) operational variables such as injection rate. A sensitivity study reveals a number of interesting observations resulting from these parameters on the shear stimulation in natural fracture system. This work strongly links the production technology, geomechanical evaluation and aids in the understanding and optimization of hydraulic fracturing simulations in naturally fractured reservoirs.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Review of Modeling Approaches to Groundwater Flow in Deformed Carbonate Aquifers

We discuss techniques to represent groundwater flow in carbonate aquifers using the three existing modeling approaches: equivalent porous medium, conduit network, and discrete fracture network. Fractures in faulted stratigraphic successions are characterized by dominant sets of sub-vertical joints. Grid rotation is recommended using the equivalent porous medium to match higher hydraulic conductivity with the dominant orientation of the joints. Modeling carbonate faults with throws greater than approximately 100 m is more challenging. Such faults are characterized by combined conduit-barrier behavior. The barrier behavior can be modeled using the Horizontal Flow Barrier Package with a low-permeability vertical barrier inserted to represent the impediment of horizontal flow in faults characterized by sharp drops of the piezometric surface. Cavities can occur parallel to the strike of normal faults generating channels for the groundwater. In this case, flow models need to account for turbulence using a conduit network approach. Channels need to be embedded in an equivalent porous medium due to cavities a few centimeters large, which are present in carbonate aquifers even in areas characterized by low hydraulic gradients. Discrete fracture network modeling enables representation of individual rock discontinuities in three dimensions. This approach is used in non-heavily karstified aquifers at industrial sites and was recently combined with the equivalent porous medium to simulate diffusivity in the matrix. Following this review, we recommend that the future research combines three practiced modeling approaches: equivalent porous medium, discrete fracture network, and conduit network, in order to capture structural and flow aspects in the modeling of groundwater in carbonate rocks.