Seismic imaging and statistical analysis of fault facies models

2017 ◽  
Vol 5 (4) ◽  
pp. SP71-SP82 ◽  
Author(s):  
Dmitriy R. Kolyukhin ◽  
Vadim V. Lisitsa ◽  
Maxim I. Protasov ◽  
Dongfang Qu ◽  
Galina V. Reshetova ◽  
...  

Interpretation of seismic responses from subsurface fault zones is hampered by the fact that the geologic structure and property distributions of fault zones can generally not be directly observed. This shortcoming curtails the use of seismic data for characterizing internal structure and properties of fault zones, and it has instead promoted the use of interpretation techniques that tend to simplify actual structural complexity by rendering faults as lines and planes rather than volumes of deformed rock. Facilitating the correlation of rock properties and seismic images of fault zones would enable active use of these images for interpreting fault zones, which in turn would improve our ability to assess the impact of fault zones on subsurface fluid flow. We use a combination of 3D fault zone models, based on empirical data and 2D forward seismic modeling to investigate the link between fault zone properties and seismic response. A comparison of spatial statistics from the geologic models and the seismic images was carried out to study how well seismic images render the modeled geologic features. Our results indicate the feasibility of extracting information about fault zone structure from seismic data by the methods used.

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>


2020 ◽  
Author(s):  
Doris Rupprecht ◽  
Sven Fuchs ◽  
Andrea Förster ◽  
Mariella Penz-Wolfmayr

<p>The GeoTief EXPLORE project aims to explore the geothermal potential and quantify the geothermal resources of the Vienna Basin (Austria) and the underlying Northern Calcareous Alpine basement. The main target of geothermal interest is the massive and tectonically remolded Hauptdolomite facies that has been identified as potential geothermal reservoir in previous studies. Now, this formation is studied using outcrop analogues for the investigation of their petrophysical characterization and specific thermal properties (thermal conductivity and thermal diffusivity).</p><p> </p><p>Here, we report new measurements on a total of 60 samples from 6 outcrops in and around the area of Vienna applying different methods for the laboratory measurement of thermal and hydraulic rock properties. The petrophysical analysis considers the impact of deformation along and across fault zones, which introduces heterogeneity of storage properties and consequently in the thermophysical properties. Using the standard fault core and damage zone model, outcrop samples were grouped into unfractured and fractured protoliths, as well as in fault rocks, like breccias and cataclasites. Rock samples are then classified by their fracture density (m² fracture surface per m³ rock) and by their matrix content and differences in grain sizes, respectively.</p><p> </p><p>The measured thermal rock properties vary significantly between the selected rock groups. The total range [90 % of values] is between 3.2 and 5.0 W/(mK) for thermal conductivity and between 1.3 and 2.7 mm²/s for thermal diffusivity. The results generally met the expected trend for fractured rocks as conductivity and diffusivity decreases with increasing porosity under unsaturated and saturated conditions. The total porosities are less than 5%. The variability of thermal conductivity under saturated conditions shows complex trends depending on the different rock classifications where fault rocks and highly fractured rocks of the damage zone show lower increase in thermal conductivities.</p><p> </p><p>The new petrophysical characterization will be the base for further numerical investigations of the hydraulic and thermal regime as well as for the analysis of the geothermal resources of the Hauptdolomite.</p><p> </p><p> </p><p> </p><p> </p>


Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. U139-U149
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali ◽  
Yi Luo

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.


2020 ◽  
Author(s):  
Laurent Guillou-Frottier ◽  
Hugo Duwiquet ◽  
Gaëtan Launay ◽  
Audrey Taillefer ◽  
Vincent Roche ◽  
...  

Abstract. In the first kilometres of the subsurface, temperature anomalies due to heat conduction processes rarely exceed 20–30 °C. When fault zones are sufficiently permeable, fluid flow may lead to thermal anomalies much higher, as evidenced by the emergence of thermal springs or by fault-related geothermal reservoirs. Hydrothermal convection triggered by buoyancy effects creates thermal anomalies whose morphology and amplitude are not well known, especially when depth- and time-dependent permeability are considered. Exploitation of shallow thermal anomalies for heat and power production partly depends on the volume and on the temperature of the hydrothermal reservoir. This study presents a non-exhaustive numerical investigation of fluid flow models within and around simplified fault zones, where realistic fluid and rock properties are accounted for, as well as appropriate boundary conditions. 2D simplified models point out relevant physical mechanisms for geological problems, such as thermal inheritance or splitting plumes showing a pulsating behaviour. When permeability is increased, the classic finger-like upwellings evolve towards a bulb-like geometry, resulting in a large volume of hot fluid at shallow depth. In the simplified 3D models, where fault zone dip angle and fault zone thickness are varied, the anomalously hot reservoir exhibits a kilometre-sized hot air balloon morphology, or, when permeability is depth-dependent, a funnel-shape geometry. For thick faults, the number of thermal anomalies increases but not the amplitude. The largest amplitude (up to 80–90 °C) is obtained for vertical fault zones. At the top of a vertical, 100 m wide, fault zone, temperature anomalies greater than 30 °C may extend laterally over more than 1 km from the fault boundary. These preliminary results should motivate further geothermal investigations of more elaborated models where topography and fault intersections would be accounted for.


2021 ◽  
Vol 9 ◽  
Author(s):  
S.C. Fabbri ◽  
C. Affentranger ◽  
S. Krastel ◽  
K. Lindhorst ◽  
M. Wessels ◽  
...  

Probabilistic seismic hazard assessments are primarily based on instrumentally recorded and historically documented earthquakes. For the northern part of the European Alpine Arc, slow crustal deformation results in low earthquake recurrence rates and brings up the necessity to extend our perspective beyond the existing earthquake catalog. The overdeepened basin of Lake Constance (Austria, Germany, and Switzerland), located within the North-Alpine Molasse Basin, is investigated as an ideal (neo-) tectonic archive. The lake is surrounded by major tectonic structures and constrained via the North Alpine Front in the South, the Jura fold-and-thrust belt in the West, and the Hegau-Lake Constance Graben System in the North. Several fault zones reach Lake Constance such as the St. Gallen Fault Zone, a reactivated basement-rooted normal fault, active during several phases from the Permo-Carboniferous to the Mesozoic. To extend the catalog of potentially active fault zones, we compiled an extensive 445 km of multi-channel reflection seismic data in 2017, complementing a moderate-size GI-airgun survey from 2016. The two datasets reveal the complete overdeepened Quaternary trough and its sedimentary infill and the upper part of the Miocene Molasse bedrock. They additionally complement existing seismic vintages that investigated the mass-transport deposit chronology and Mesozoic fault structures. The compilation of 2D seismic data allowed investigating the seismic stratigraphy of the Quaternary infill and its underlying bedrock of Lake Constance, shaped by multiple glaciations. The 2D seismic sections revealed 154 fault indications in the Obersee Basin and 39 fault indications in the Untersee Basin. Their interpretative linkage results in 23 and five major fault planes, respectively. One of the major fault planes, traceable to Cenozoic bedrock, is associated with a prominent offset of the lake bottom on the multibeam bathymetric map. Across this area, high-resolution single channel data was acquired and a transect of five short cores was retrieved displaying significant sediment thickness changes across the seismically mapped fault trace with a surface-rupture related turbidite, all indicating repeated activity of a likely seismogenic strike-slip fault with a normal faulting component. We interpret this fault as northward continuation of the St. Gallen Fault Zone, previously described onshore on 3D seismic data.


Author(s):  
Jennifer Cunningham ◽  
Wiktor Weibull ◽  
Nestor Cardozo ◽  
David Iacopini

<p>PS seismic data from the Snøhvit field are compared with forward seismic modelling to understand the effect of azimuthal separation and incidence angle on the imaging of faults. Two faults, one oriented oblique to the survey and one approximately parallel to the survey were chosen. Azimuthally separated W (source is W relative to receivers) and E (source E relative to receivers) data demonstrate that fault imaging is more affected by azimuth when the faults are oblique to the survey orientation, and W data image the faults better. Partial stack data show that with increasing incidence angle there is a systematic improvement in the quality of fault imaging in both the E and W data. In addition, the frequency content of seismic waves back-scattered from within and around fault zones is analysed in the Snøhvit data. Low-medium frequencies are dominant within fault zones, compared with higher frequencies in adjacent areas and haloes of medium frequencies surrounding the faults. Two synthetic experiments support the azimuth, incidence angle and frequency observations. In the first experiment, the fault is modelled as a planar discontinuity and the data were processed in the same way as the Snøhvit data (into separate azimuths and incidence angle stacks). The first experiment confirms a strengthening in the seismic signal from faults in the W data. This is due to the interaction of specular waves and diffractions which are more abundant in the W data. The second experiment had three parts modelling the fault zone with different layering complexity. It proved that frequencies in the fault and adjacent areas increase with fault zone complexity, and that the internal architecture of faults can impact the frequencies in the data adjacent to faults. </p>


2020 ◽  
Author(s):  
Hagen Deckert ◽  
Steffen Abe ◽  
Wolfgang Bauer

<p>In the course of hydrocarbon or geothermal exploration the characterisation of fault zone architectures is of interest for fluid flow modelling and geomechanical studies. Seismic data normally offer the best information for the identification of fault zone geometries in sedimentary basins. However, the internal structure or the damage zone of a fault can be hardly resolved with seismic data as displacements along single fault strands or fractures are by far too small. Thus, it is not possible to directly map small scale faults with seismic methods, though these structures might significantly influence fluid flow. We try to examine the architecture of extensional fault zones in carbonate rocks at subseismic scales by using discrete element method (DEM) techniques to numerically simulate the evolution of fault zones including their associated damage zones.</p><p>As a case study we have analysed the geometry, displacement and fault width of normal faults in fine grained jurassic limestones in a quarry in Franconia, Germany. The quarry shows a rather simple set of conjugated 60deg dipping normal faults. Displacement is rather small and varies between c. 5cm up to c. 2m, some faults show almost no offset. The fault thickness varies between 2cm and c. 1m. A closer investigation of the fault geometries reveals, next to planar parts, sometimes complex fault zone structures including restraining and releasing bends, multiple fault strands as well as lenses and associated riedel shears. Analysis of high resolution photogrammetric data revealed a high number of small scale fractures between neighbouring discrete fault surfaces which are interpreted as highly fractured damage zones. Some faults with rather small displacement suggest that the overall inclination of the fault is a result of small subvertical sections which are connected in a staircase like appearance. </p><p>The DEM models simulate normal faulting in a layered marl-limestone sequence driven by the displacement of an underlying basement fault. Different layer geometries and effective vertical stresses in the range of 15-45 MPa, equivalent to an overburden thickness of c. 1000-3000m, have been used in the models. The stress range covers the maximum burial depth of the carbonates, which is assumed to be c. 1500m. Material properties used in the DEM were calibrated based on laboratory data, i.e. results of triaxial deformation tests on the studied limestones.</p><p>Results of the models show fault geometries which resemble those observed in the studied outcrop. In particularly under low stress, small offsets and with strongly decoupled layers we observe steeply dipping faults (>70deg) which also show staircase structures composed of sub-vertical fractures within each of the layers and horizontal offsets along the layer interfaces. We also observe the development of multiple fault strands and associated damage zones. </p><p>Our study shows that the DEM models are capable to reproduce observed fault geometries and damage zones. The results help to understand fault zone architectures and depict highly fractured areas in a sub-seismic scale.</p>


2018 ◽  
Vol 9 (3) ◽  
pp. 629-652
Author(s):  
Yu. L. Rebetsky

The problem of forecasting seismic hazards is discussed. The stress state data characterizing various aspects of brittle failure are reviewed in detail. It is shown that the most convenient tool for analyzing such data is the Mohr stress diagram and the Coulomb criterion. Noted is the role of a fluid in not only reducing the normal stresses responsible for brittle failure, but also predetermining the major processes in fault zones. In each fault body, a node can be distinguished as a fault part wherein the main structural and material transformations take place. The node contains narrow elongated zones of modification of mylonites, from protomylonites to ultramylonites and blastomylonites, that are related to the localization of continuous and discontinuous shear deformation. Due to the metamorphic processes, fault zones are less strong than the surrounding consolidated blocks of the crust. A theoretical analysis of the mechanism of displacements along the discontinuities of different scale ranks shows differences in their manifestation. Tectonic and seismic displacements along the rupture occupy the entire area at once, while displacements along the fault zone occur in stages along its extent and follow the ‘rolling-carpet’ principle that is also typical of intra-crystal dislocations. The stress state in the vicinity of ruptures and faults has different characteristic features. Based on the seismological and tectonophysical data on earthquake focal parameters and discontinuities, it is possible to identify two or three ranks of stresses, which differ in the laws predetermining their mutual relationships. Actually, this conclusion contradicts the hypothesis of self-similarity of discontinuities in their continuous range, from a dislocation to a fault zone, which length amounts to tens of kilometers. Besides, it imposes a restriction on the use of statistical analysis of seismic data. The seismic data show that in the source of a large earthquake, displacement develops as a running band (‘rolling-carpet’ principle). In the source of a weak earthquake, it occupies the entire earthquake focal area at once. The differences in the types of shearing in the sources of weak and strong earthquakes are related to the relationships between three dynamic parameters of the medium: velocity of seismic wave propagation, rate of rupture propagation, and displacement rate of the sides of the fracture. Using tectonophysical methods, the stress state was reconstructed for the seismically active regions of the planet and the sources of the mega-earthquakes of the 21st century. Based on the reconstructions, the mean strength and stress values were calculated, and the specific features of the stress fields were revealed. It is established that the strongest regional earthquakes ‘avoided’ the areas with increased effective isotropic pressure. The sizes of the sources of the strongest earthquakes were controlled by the size of the region with decreased effective pressure. The sites, wherefrom the earthquake were initiated, were often located in the zones of the highest stress gradients. These regularities support the term “metastability of the state of fault zone” (introduced to seismology from the physics of the states of matter) and justify it by a specific distribution pattern of stress values prior to the mega-earthquake. Based on the tectonophysical definition of the metastable state of faults, the important role is outlined for a stress gradient zone that represents a location wherein a trigger earthquake occurs. The ‘maturity’ of the zone with increased stress gradient values is, in essence, a characteristic of the time interval of metastability of the fault zone.


2019 ◽  
Vol 496 (1) ◽  
pp. 39-73 ◽  
Author(s):  
M. Nogueira Kiewiet ◽  
C. Lima ◽  
A. Giwelli ◽  
C. Delle Piane ◽  
V. Lemiale ◽  
...  

AbstractTo contribute to the understanding of the impacts of fault reactivation induced by reservoir exploitation, we describe the final series of laboratory experiments, numerical simulations and microstructural analysis conducted during the ‘Fault Reactivation in Carbonates’ research project. In the project, the structure and hydromechanical properties of carbonate-hosted fault zones were investigated. For the analyses here reported, faults were artificially generated by direct shearing composite blocks made of layers of reservoir analogue rocks (outcrop travertine or synthetic grainstone) intercalated with one layer of a sealing analogue rock (synthetic carbonate mudstone). Post-direct shearing, cylindrical plugs containing the fault zone and parts of intact rock were cored out from the blocks and tested in a triaxial test rig, simulating fault reactivation. Varied stress paths and pore-pressure conditions representative of fluid depletion and injection were considered. In parallel, two-dimensional mechanical models representative of the direct shear experiments were developed using smoothed particle hydrodynamics (SPH). We observed a continuous reduction in fault transmissibility during direct shearing, followed by a permeability reduction of 50–80% with increasing mean effective stress in the subsequent fault reactivation tests. Experimental fault zone geometries produced during direct shear were broadly reproduced by the two-dimensional modelling approach. We also detected that the inclusion of the carbonate mud sealing rock into the fault zone caused greater compaction of the fault materials when compared to experiments conducted without carbonate mud layers. We conclude that with fault displacement, increasing incorporation of carbonate mud sealing material into the fault zone and the concomitant development of gouge results in the continuous reduction of fault transmissibility/permeability. This occurs in the two very different limestone host-rock types and for all the stress configurations investigated. Discussions on these results and also on the outcomes of the research project as a whole are presented in the paper.


Solid Earth ◽  
2020 ◽  
Vol 11 (4) ◽  
pp. 1571-1595 ◽  
Author(s):  
Laurent Guillou-Frottier ◽  
Hugo Duwiquet ◽  
Gaëtan Launay ◽  
Audrey Taillefer ◽  
Vincent Roche ◽  
...  

Abstract. In the first kilometers of the subsurface, temperature anomalies due to heat conduction processes rarely exceed 20–30 ∘C. When fault zones are sufficiently permeable, fluid flow may lead to much larger thermal anomalies, as evidenced by the emergence of thermal springs or by fault-related geothermal reservoirs. Hydrothermal convection triggered by buoyancy effects creates thermal anomalies whose morphology and amplitude are not well known, especially when depth- and time-dependent permeability is considered. Exploitation of shallow thermal anomalies for heat and power production partly depends on the volume and temperature of the hydrothermal reservoir. This study presents a non-exhaustive numerical investigation of fluid flow models within and around simplified fault zones, wherein realistic fluid and rock properties are accounted for, as are appropriate boundary conditions. 2D simplified models point out relevant physical mechanisms for geological problems, such as “thermal inheritance” or pulsating plumes. When permeability is increased, the classic “finger-like” upwellings evolve towards a “bulb-like” geometry, resulting in a large volume of hot fluid at shallow depth. In simplified 3D models wherein the fault zone dip angle and fault zone thickness are varied, the anomalously hot reservoir exhibits a kilometer-sized “hot air balloon” morphology or, when permeability is depth-dependent, a “funnel-shaped” geometry. For thick faults, the number of thermal anomalies increases but not the amplitude. The largest amplitude (up to 80–90 ∘C) is obtained for vertical fault zones. At the top of a vertical, 100 m wide fault zone, temperature anomalies greater than 30 ∘C may extend laterally over more than 1 km from the fault boundary. These preliminary results should motivate further geothermal investigations of more elaborated models wherein topography and fault intersections would be accounted for.


Sign in / Sign up

Export Citation Format

Share Document