scholarly journals Influence of deformation-band fault damage zone on reservoir performance

2017 ◽  
Vol 5 (4) ◽  
pp. SP41-SP56 ◽  
Author(s):  
Dongfang Qu ◽  
Jan Tveranger ◽  
Muhammad Fachri

Access to 3D descriptions of fault zone architectures and recent development of modeling techniques allowing explicit rendering of these features in reservoir models, provide a new tool for detailed implementation of fault zone properties. Our aim is to assess how explicit rendering of fault zone architecture and properties affects performance of fluid flow simulation models. The test models use a fault with a maximum 100 m displacement and a fault damage zone with petrophysical heterogeneity caused by the presence of deformation bands. The distribution pattern of deformation bands in fault damage zones is well-documented, which allows generation of realistic models. A multiscale modeling workflow is applied to incorporate these features into reservoir models. Model input parameters were modulated to provide a range of property distributions, and the interplay between the modeling parameters and reservoir performance was analyzed. The influence of deformation-band damage zone on reservoir performance in the presence of different fault core transmissibility-multipliers was investigated. Two configurations are considered: one in which the fault terminates inside the model domain, representing a case in which the fluid can flow around the fault, and one in which the fault dissects the entire model domain, representing a case in which the fluid is forced to cross the fault. We observed that the impact of deformation-band fault damage zone on reservoir performance changes when the fault core transmissibility multiplier is changed. Reservoir performance is insensitive to changing damage zone heterogeneity in a configuration in which flow can move around the fault. Where flow cannot bypass the fault, the influence of fault damage zone heterogeneity on reservoir performance is significant even when the fault core transmissibility multiplier is low.

2019 ◽  
Vol 132 (5-6) ◽  
pp. 1183-1200 ◽  
Author(s):  
Mattia Pizzati ◽  
Fabrizio Balsamo ◽  
Fabrizio Storti ◽  
Paola Iacumin

Abstract In this work, we report the results of a multidisciplinary study describing the structural architecture and diagenetic evolution of the Rocca di Neto extensional fault zone developed in poorly lithified sandstones of the Crotone Basin, Southern Italy. The studied fault zone has an estimated displacement of ∼90 m and consists of: (1) a low-deformation zone with subsidiary faults and widely spaced deformation bands; (2) an ∼10-m-wide damage zone, characterized by a dense network of conjugate deformation bands; (3) an ∼3-m-wide mixed zone produced by tectonic mixing of sediments with different grain size; (4) an ∼1-m-wide fault core with bedding transposed into foliation and ultra-comminute black gouge layers. Microstructural investigations indicate that particulate flow was the dominant early-stage deformation mechanism, while cataclasis became predominant after porosity loss, shallow burial, and selective calcite cementation. The combination of tectonic compaction and preferential cementation led to a strain-hardening behavior inducing the formation of “inclined conjugate deformation band sets” inside the damage zone, caused by the kinematic stress field associated with fault activity. Conversely, conjugate deformation band sets with a vertical bisector formed outside the damage zone in response to the regional extensional stress field. Stable isotope analysis helped in constraining the diagenetic environment of deformation, which is characterized by mixed marine-meteoric signature for cements hosted inside the damage zone, while it progressively becomes more meteoric moving outside the fault zone. This evidence supports the outward propagation of fault-related deformation structures in the footwall damage zone.


2017 ◽  
Vol 5 (4) ◽  
pp. SP9-SP26 ◽  
Author(s):  
Charlotte Botter ◽  
Nestor Cardozo ◽  
Dongfang Qu ◽  
Jan Tveranger ◽  
Dmitriy Kolyukhin

Faults play a key role in reservoirs by enhancing or restricting fluid flow. A fault zone can be divided into a fault core that accommodates most of the displacement and a surrounding damage zone. Interpretation of seismic data is a key method for studying subsurface features, but the internal structure and properties of fault zones are often at the limit of seismic resolution. We have investigated the seismic response of a vertical fault zone model in sandstone, populated with fault facies based on deformation band distributions. Deformation bands reduce the porosity of the sandstone, and they condition its elastic properties. We generate synthetic seismic cubes of the fault facies model for several wave frequencies and under realistic conditions of reservoir burial and seismic acquisition. Seismic image quality and fault zone definition are highly dependent on wave frequency. At a low wave frequency (e.g., 10 Hz), the fault zone is broader and no information about its fault facies distribution can be extracted. At higher wave frequencies (e.g., 30 and 60 Hz), seismic attributes, such as tensor and envelope, can be used to characterize the fault volume and its internal structure. Based on these attributes, we can subdivide the fault zone into several seismic facies from the core to the damage zone. Statistical analyses indicate a correlation between the seismic attributes and the fault internal structure, although seismic facies, due to their coarser resolution, cannot be matched to individual fault facies. The seismic facies can be used as input for reservoir models as spatial conditioning parameters for fault facies distributions inside the fault zone. However, relying only on the information provided by seismic analyses might not be enough to create high-resolution fault reservoir models.


2021 ◽  
Author(s):  
Irène Aubert ◽  
Juliette Lamarche ◽  
Philippe Leonide

<p>Understanding the impact of fault zones on reservoir trap properties is a major challenge for a variety of geological ressources applications. Fault zones in cohesive rocks are complex structures, composed of 3 components: rock matrix, damage zone fractures and fault core rock. Despite the diversity of existing methods to estimate fault zone permeability/drain properties, up to date none of them integrate simultaneously the 3 components of fracture, fault core and matrix permeability, neither their evolution with time. We present a ternary plot that characterizes the fault zones permeability as well as their drainage properties. The ternary plot aims at (i) characterizing the fault zone permeability between the three vertices of matrix, fractures and fault core permeability ; and at (ii) defining the drain properties among 4 possible hydraulic system: (I) good horizontal and vertical, fault-perpendicular and -parallel; (II) moderate parallel fluid pathway; (III) good parallel fault-core and (IV) good parallel fractures. The ternary plot method is valid for 3 and 2 components fault zones. The application to the Castellas Fault case study show the simplicity and efficiency of the plot for studying underground and/or fossil, simple or polyphase faults in reservoirs with complete or limited permeability data.</p>


2009 ◽  
Vol 14 (2) ◽  
pp. 231-248 ◽  
Author(s):  
Dmitry Kolyukhin ◽  
Sylvie Schueller ◽  
Magne S. Espedal ◽  
Haakon Fossen

2020 ◽  
Author(s):  
Alba Zappone ◽  
Antonio Pio Rinaldi ◽  
Melchior Grab ◽  
Quinn Wenning ◽  
Clément Roques ◽  
...  

Abstract. The success of geological carbon storage depends on the assurance of a permanent confinement of the injected CO2 in the storage formation at depth. One of the critical elements of the safekeeping of CO2 is the sealing capacity of the caprock overlying the storage formation, despite faults and/or fractures, which may occur in it. In this work, we present an ongoing injection experiment performed in a fault hosted in clay at the Mont Terri underground rock laboratory (NW Switzerland). The experiment aims at improving our understanding on the main physical and chemical mechanisms controlling i) the migration of CO2 through a fault damage zone, ii) the interaction of the CO2 with the neighbouring intact rock, and iii) the impact of the injection on the transmissivity in the fault. To this end, we inject a CO2-saturated saline water in the top of a 3 m think fault in the Opalinus Clay, a clay formation that is a good analogue of common caprock for CO2 storage at depth. The mobility of the CO2 within the fault is studied at decameter scale, by using a comprehensive monitoring system. Our experiment aims to the closing of the knowledge gap between laboratory and reservoir scales. Therefore, an important aspect of the experiment is the decameter scale and the prolonged duration of observations over many months. We collect observations and data from a wide range of monitoring systems, such as a seismic network, pressure temperature and electrical conductivity sensors, fiber optics, extensometers, and an in situ mass spectrometer for dissolved gas monitoring. The observations are complemented by laboratory data on collected fluids and rock samples. Here we show the details of the experimental concept and installed instrumentation, as well as the first results of the preliminary characterization. Analysis of borehole logging allow identifying potential hydraulic transmissive structures within the fault zone. A preliminary analysis of the injection tests helped estimating the transmissivity of such structures within the fault zone, as well as the pressure required to mechanically open such features. The preliminary tests did not record any induced microseismic events. Active seismic tomography enabled a sharp imaging the fault zone.


2021 ◽  
Author(s):  
Anita Torabi ◽  
Behzad Alaei ◽  
Audun Libak

<p>Understanding fault geometry and processes of faulting are important research areas for many applications such as petroleum exploration and production; geothermal energy managements; hydrogeology; waste disposal and CO2 storage underground; earthquake seismology and geological hazard studies. Faults can be described as comprising a core and an enveloping damage zone (e.g. Caine et al. 1996).  The fault core accommodates most of the displacement along multiple slip surfaces and may include fault rocks such as fault gouge, cataclasites, breccia, clay smear, fractures, diagenetic features, and lenses of deformed and undeformed rocks trapped between slip surfaces. Whereas, the deformation is less intense in the damage zone and may include fractures and/or deformation bands depending on the initial porosity of the host rock, minor faults, and folds (Torabi et al., 2020). Fault geometric attributes include fault shape, fault displacement, length, damage zone width and fault core thickness (Caine et al., 1996; Torabi and Berg, 2011). Currently, there are uncertainties in defining and understanding of fault 3D geometry. These uncertainties are to some extent related to the accessibility of the fault geometric attributes and the methodological constraints, utilizing biased data. Details of fault damage zone and fault core structures can be mapped at outcrop, however, their descriptions and statistical handling are usually constrained by their accessibility in the field and their definitions by individual researchers.</p><p>Reflection seismic data is used to study faults in the subsurface, although the interpretation of faults could be affected by the seismic resolution and the accuracy of interpretation (Marchal et al., 2003; Lohr et al., 2008; Iacopini et al., 2016; Torabi et al., 2016). Utilizing seismic attributes, we are able to directly images faults from seismic without a need for interpretation. Using this method, we extracted fault geometric attributes directly from fault images in the fault attribute volumes and studied the 3D shape and displacement distribution of faults (Torabi et al., 2019). By integrating spectral decomposition with seismic attribute workflows, we created enhanced fault attribute volumes with a high resolution, allowing us to detect, and map fault damaged zone (fault damage zone plus fault core in outcrop scale) in seismic data (Alaei and Torabi, 2017). Finally, we integrated the data from outcrop and seismic study in the scaling relations between the faults geometric attributes in order to predict the fault geometry in the subsurface.</p><p> </p><p> </p>


Solid Earth ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 319-343
Author(s):  
Alba Zappone ◽  
Antonio Pio Rinaldi ◽  
Melchior Grab ◽  
Quinn C. Wenning ◽  
Clément Roques ◽  
...  

Abstract. The success of geological carbon storage depends on the assurance of permanent containment for injected carbon dioxide (CO2) in the storage formation at depth. One of the critical elements of the safekeeping of CO2 is the sealing capacity of the caprock overlying the storage formation despite faults and/or fractures, which may occur in it. In this work, we present an ongoing injection experiment performed in a fault hosted in clay at the Mont Terri underground rock laboratory (NW Switzerland). The experiment aims to improve our understanding of the main physical and chemical mechanisms controlling (i) the migration of CO2 through a fault damage zone, (ii) the interaction of the CO2 with the neighboring intact rock, and (iii) the impact of the injection on the transmissivity in the fault. To this end, we inject CO2-saturated saline water in the top of a 3 m thick fault in the Opalinus Clay, a clay formation that is a good analog of common caprock for CO2 storage at depth. The mobility of the CO2 within the fault is studied at the decameter scale by using a comprehensive monitoring system. Our experiment aims to close the knowledge gap between laboratory and reservoir scales. Therefore, an important aspect of the experiment is the decameter scale and the prolonged duration of observations over many months. We collect observations and data from a wide range of monitoring systems, such as a seismic network, pressure temperature and electrical conductivity sensors, fiber optics, extensometers, and an in situ mass spectrometer for dissolved gas monitoring. The observations are complemented by laboratory data on collected fluids and rock samples. Here we show the details of the experimental concept and installed instrumentation, as well as the first results of the preliminary characterization. An analysis of borehole logging allows for identifying potential hydraulic transmissive structures within the fault zone. A preliminary analysis of the injection tests helped estimate the transmissivity of such structures within the fault zone and the pressure required to mechanically open such features. The preliminary tests did not record any induced microseismic events. Active seismic tomography enabled sharp imaging the fault zone.


SPE Journal ◽  
2007 ◽  
Vol 12 (02) ◽  
pp. 167-178 ◽  
Author(s):  
Mohsen Masihi ◽  
Peter Robert King ◽  
Peyman Reza Nurafza

Summary Investigating the impact of geological uncertainty (i.e., spatial distribution of fractures) on reservoir performance may aid management decisions. The conventional approach to address this is to build a number of possible reservoir models, upscale them, and then run flow simulations. The problem with this approach is that it is computationally very expensive. In this study, we use another approach based on the permeability contrasts that control the flow, called percolation approach. This assumes that the permeability disorder of a rock can be simplified to either permeable or impermeable. The advantage is that by using some universal laws from percolation theory, the effect of the complex geometry which influences the global properties (e.g., connectivity or conductivity) can be easily estimated in a fraction of a second on a spread sheet. The aim of this contribution is to establish the percolation framework to examine the connectivity of fracture systems at a given finite observation scale in 2D and 3D. In particular, we use numerical simulation to show how the scaling laws of the connectivity derived originally for constant-length isotropic systems can be expanded to cover more realistic cases including fracture systems with anisotropy and fracture-length distribution. Finally, the outcrop data of mineralized fractures exposed on the southern margin of the Bristol Channel Basin was used to show that the predictions from the percolation approach are in agreement with the results calculated from field data but can be obtained very quickly. As a result, this may be used for practical engineering purposes for decision making. Introduction Fractured reservoirs are very complex, containing geological heterogeneities (i.e., fractures) on various length scales from microns to kilometers. These heterogeneities have significant impact on the flow behavior and have to be modeled to make reliable prediction of reservoir performance. However, we have very few direct measurements of the flow properties (e.g., core and image-log data) that are 1D and represent a very small volume of a typical reservoir. Other type of data are more widespread (e.g., well-test or seismic data) but generally are related indirectly to fracture distribution. The consequence is a great deal of uncertainty about the spatial distribution of the fractures that influence the flow and affect the reservoir performance. A major factor in analysis of flow and transport in these reservoirs is the appropriate representation of the heterogeneities that control flow (Bear et al. 1993). The conventional approach to investigate the impact of geological uncertainty on reservoir recovery is to build a detailed reservoir model using geophysical and geological data, upscale it, and then perform flow simulation. This is typically done by assuming either equivalent continuum models (i.e., dual porosity or dual permeability), discrete network models, or an integration of both (Warren and Root 1963; Dershowitz et al. 2000). The fractures can be assumed to be infinite (Snow 1969), which means that they are perfectly connected, or finite in length (Sagar and Runchal 1982; Long and Witherspoon 1985). If fractures are poorly interconnected and the matrix rock is relatively impermeable, the network formed by the fractures may control the flow. On the other hand, if the matrix is relatively permeable and the fractures are regular and highly interconnected, fractures and matrix could be treated as two separate continuums occupying the entire domain (Warren and Root 1963). In order to have a reliable estimation of reservoir performance parameters, it is necessary to construct a number of possible reservoir models (with associated probabilities) and then run flow simulations many times. The problem with this approach is that it is computationally very expensive. Therefore, there is a great incentive to produce much simpler physically-based models to predict uncertainty in performance very quickly, especially for engineering purposes.


2020 ◽  
Author(s):  
Ulrich Kelka ◽  
Thomas Poulet ◽  
Luk Peeters

<p>Fault and fracture networks can govern fluid flow patterns in the subsurface and predicting fluid flow on a regional scale is of interest in a variety of fields like groundwater management, mining engineering, energy, and mineral resources. Especially the pore fluid pressure can have a strong impact on the strength of fault zones and might be one of the drivers for fault reactivation. Reliable simulations of the transient changes in fluid pressure need to account for the generic architecture of fault zones that comprises strong permeability contrast between the fault core and damage zone.</p><p>Particularly, the distribution and connectivity of large-scale fault zones can have a strong impact on the flow field. Yet, modelling numerically such features in their full complexity remains challenging. Often faults zones are conceptualized as forming exclusively either barriers or conduits to fluid flow. However, a generic architecture of fault zones often comprises a discrete fault core surrounded by a diffuse damage zone and conceptualizing large scale discontinuities simply as a barrier or conduit is unlikely to capture the regional scale fluid flow dynamics. It is known that if the fault zone is hosted in low-permeability strata, such as clays or crystalline rocks, a transversal flow barrier can form along the fault core whereas the fracture-rich fault damage zone represents a longitudinal conduit. In more permeable host-rocks (i.e. sandstones or carbonates) the reverse situation can occur, and the permeability distributions in the damage zones can be governed by the abundance of low-permeability deformation features. A reliable numerical model needs to account for the difference and strong contrasts in fluid flow properties of the core and the damage zone, both transversally and longitudinally, in order to make prediction about the regional fluid flow pattern.</p><p>Here, we present a numerical method that accounts for the generic fault zone architecture as lower dimensional interfaces in conforming meshes during fluid flow simulations in fault networks. With this method we aim to decipher the impact of fault zone architecture on subsurface flow pattern and fluid pressure evolution in fractures and faulted porous media. The method is implemented in a finite element framework for Multiphysics simulations. We demonstrate the impact of considering the more generic geological structure of individual faults on the flow field by conceptualizing discontinuities either as barriers, conduits or as a conduit-barrier system and show were these conceptualizations are applicable in natural systems. We further show that a reliable regional scale fluid flow simulation in faulted porous media needs to account for the generic fault zone architecture. The approach is finally used to evaluate the fluid flow response of statistically parameterised faulted media, in order to investigate the impact and sensitivity of each variable parameter.</p>


Sign in / Sign up

Export Citation Format

Share Document