Fractured and karstified carbonate units are key exploration targets for the hydrocarbon industry as they represent important reservoirs. Furthermore, large water reserves and geothermal systems are hosted in carbonate aquifers. This paper documents the relationships between stratigraphy, structural patterns, silicification, and the spatial-morphological organization of a 3D multistorey cave system developed in a Neoproterozoic mixed carbonate-siliciclastic sequence. We found that the combination of lithology, silicification, fracture patterns (controlled by lithostratigraphic variability), and petrophysical properties control the formation of high or low permeability zones; their distribution was fundamental for the spatial organization of dissolution and the compartmentalization of the resulting conduit system in different speleogenetic storeys. We propose a deep-seated hydrothermal origin for the fluids involved in the main phases of karst formation. Warm and alkaline hydrothermal fluids caused silica dissolution, followed by chalcedony and quartz reprecipitation in pore space and fractures. Rising fluids concentrated along through-going vertical fracture zones in the lower storey, whereas sub-horizontal bedding-parallel fluid flow was focused on sedimentary packages containing highly silicified dolostones (SiO₂>80 wt%) characterized by high permeability. The Calixto Cave is an enlightening example for the complex speleogenetic history affecting a mixed carbonate-siliciclastic succession where the combined effect of silicification and hydrothermal karst dissolution can potentially generate high-quality reservoirs.

@en

The present study used a multitool approach to characterize fractures of several orders of magnitude in large fracture corridors, caves, and canyons to investigate their impact on fluid flow in carbonate units. The study area is the Brejões carbonate karst system that is located in the Neoproterozoic Salitre Formation in the Irecê Basin, São Francisco Craton, Brazil. The approach included satellite imagery, used for interpreting the regional structural context, Unmanned Aerial Vehicle (UAV) and ground-based Light Detection And Ranging (LiDAR) imagery, used for detailed structural interpretation. Regional interpretation revealed that fracture corridors, caves and canyons occur along a N–S-oriented anticline hinge. An advanced stage of karstification caused fracture enlargement and intrabed dissolution, and the formation of caves and canyons. A river captured by the highly fractured zone along the anticline hinge played an important role as an erosive agent. Detailed characterization of fracture corridors comprised structural analysis, topological studies, persistence estimations, power-law fitting of fracture trace length distributions, and identification of network backbones. Our results indicate that fracture corridors comprise four subvertical fracture sets: N–S and E-W and a conjugate pair, NNE-SSW and NW-SE. Fractures observed in the caves show the same dominant directions. Fracture directions are consistent with a common origin associated with the anticline folding. Fracture traces range from 1.0 m to 300 m, comprising both subseismic (<50 m) and seismic scale fractures (>50 m). Networks have dominance of node terminations Y and X (notably Y), C_B values higher than 1.8, high P₂₀ and P₂₁ persistence values, and highly interconnected backbones. Fracture network connectivity is associated with power-law exponents greater than 2.5 for the fracture trace distributions, indicating large influence of subseismic-scale fractures on fluid flow. As the final result of folding and karstification, large volumes of secondary macroporosity were created, particularly in the zone of maximum fracture intensity around the hinge zone of the anticline. This scenario can be used to understand better oil reservoirs formed in similar structural controls in near-surface conditions.

@en

Caves developed in carbonate units have a significant role in fluid flow, but most of these subsurface voids are below seismic resolution. We concentrated our study on four caves to determine the roles of fractures and folds in the development of karst conduits that may form flow pathways in carbonate reservoirs. We performed structural field investigations, petrographic analyses, and geometric characterization using Light Detection and Ranging (LIDAR) for caves in Neoproterozoic carbonates of the Salitre Formation, central part of the São Francisco Craton, Brazil. We found that the conduit shape, usually with an ellipsoidal cross-section, reflects the tectonic features and textural variations. Carbonate layers containing pyrite and low detritic mineral contents are generally karstified and appear to act as favorable flow pathways. Our results indicate that the development of the karst system is related to fracture corridors formed along parallel and orthogonal sets of fold hinges, which provide preferential pathways for fluid flow and contribute to the development of super-K zones. This study provides insights into the prediction of subseismic-scale voids in carbonate reservoirs, with direct application for the hydrocarbon and hydrogeology flow and storage.

@en

This study combines multiscale analyses of geological, fault, fracture, and stable isotope data to investigate strike-slip deformation and channeling of hydrothermal fluids along the Cafarnaum fault and calcite veins at different distances from the fault, which is a structure in the São Francisco Craton, northeastern Brazil. Meteoric fluids with δD values near −45‰ and δ¹⁸O values near −6.5‰ and temperatures at 40–70 °C precipitated as calcite veins in the host carbonate units. The Cafarnaum fault, a N-S-striking vertical, ~170 km long fault zone, juxtaposes Neoproterozoic carbonate rocks in the western block and Mesoproterozoic siliciclastic rocks in the eastern block. A zone of restraining bends occurs at the central part of the fault, whereas termination zones of horsetail geometry occur at both ends of the Cafarnaum fault. These zones are marked by NW-SE-striking extensional faults that are oblique to the main N-S-striking fault zone, where hydrothermal deposits occur. The zone of influence of the Cafarnaum fault is ~ 20 km wide around the main fault. The fault formed during the Brasiliano orogeny (740–560 Ma) after Neoproterozoic carbonate platform deposition. In contrast with the host units, fluids along the fault zone originated in deeper levels of the crust and show much lower δ¹⁸O values, indicating higher crystallization temperatures. These fluids caused brecciation in the Neoproterozoic carbonate host rocks, whereas a subsequent decrease in fluid pressure and cooling near the surface resulted in the precipitation of a hydrothermal paragenesis in veins, also affecting the host rock.

@en

Rifting and related normal stress regime in the equatorial continental margin of Brazil ceased during the Late Cretaceous, when the stress regimes in eastern South America and West Africa changed to induce strike-slip or reverse motion. In this study, we explore the postrift tectonic, geomorphic, magmatic, and sedimentary responses to stress changes using the Potiguar Basin, the easternmost basin in the equatorial margin of Brazil, as a case study. We use field and topographic data, 2D seismic reflection lines, vertical electric soundings, and geochronological and borehole data to constrain the stress evolution of the Potiguar Basin from the Late Cretaceous to the Quaternary, discussing the role of basin inversion on sedimentation and landforms. Our results indicate the presence of two strike-slip stress regimes after rifting. The first stress field (SF1) occurred from Late Cretaceous to the middle Miocene and consisted of a N-S-oriented maximum subhorizontal compression and an E-W-oriented extension. The second stress field (SF2) took place from the middle Miocene to the present day and included subhorizontal E-W to NW-SE compression combined with N-S and NE-SW subhorizontal extension. Emplacement of volcanic rocks occurred along transtensional faults, with a principal peak during SF1 at 20–30 Ma and a subordinate peak during SF2 at 5–10 Ma. In response to shortening during SF2, a 70-km-long and 50-km-wide dome formed, where marine Miocene strata were uplifted to ~250 m asl. This uplift induced the displacement of alluvial channels away from the dome. Anticlines formed by transpression along the main NE-SW-striking faults during both SF1 and ST2 acted as traps in the petroleum system. Similar shifts and stress field inversions documented in other areas of the Brazilian continental margin are consistent with the Neogene rise of the Andes and may have implications for reconstructing the tectonic history of the Equatorial Atlantic margin of South America.

@en

The hypogenic caves developed in carbonate units have a significant structural control but most of their features are not detect by conventional methods due to their size below seismic resolution. This contribution focuses on the structural, petrographic and geometric characterization of karst conduits in Neoproterozoic carbonates of the Salitre Formation, central part of São Francisco Craton, Brazil. We address the influence of fractures and folds on the development of karst conduits through field and laboratory analysis and the application of Light Detection Ranging to characterize cave/conduit geometry. The preliminary results indicate that the process of karstification are intensified in fractures corridors developed along fold hinges, which create fluid flow corridors in carbonate units and may change petrophysical reservoir properties.

@en

Multiscale karst systems, typically below seismic resolution, are important contributors to permeability but also provide significant challenges in drilling operations and in flow simulations. A substantial portion of subsurface karsts is of hypogenic origin, that is, is associated with upward movement of chemically aggressive fluids. The Morro Vermelho Karst System and Cave developed within Proterozioc shallow water carbonates and is an inspiring analog for buried hypogenic karsts. MVKS formed in association with a regional strike-slip fault. In an initial stage, deformation in the “softer” carbonates overlying more rigid quartzites and basement was accommodated by a 100s of m wide semi-ductile shear zone. Wholesale dolomitization took place and numerous dolomite-filled veins formed. With progressing strike-slip, the fault propagated upward creating an anticline parallel to the fault zone, causing the development of m-long extensional fractures, the arrival of SiO2 rich fluids in the karts and, eventually, the formation of the cave. Vertical fracture corridors and/or faults capturing deep, generally layer parallel flow are present in different settings, from foredeeps to rifted margins and are primary candidates for widespread hypogenic karst development.@en

The carbonate rocks exposed in the Irecê (Brazil) are pervasively affected by hydrothermal silicification and dolomitization. These mineralizing events drastically changed the original petrophysical properties of the host rock. In order to understand the role played by deformation mechanisms and tectonic evolution of these carbonates on the ingression of the mineralizing fluids, we focused our study on the mineralogical, petrographic, stratigraphic, structural and geochemical characterization of the faulted Neoproterozoic carbonates of the Salitre Formation, central Brazil. In this contribution, we present the first results of a combined outcrop- to micro-scale investigation of the structural evolution of the fracture network responsible for hydrothermal silicification and dolomitization processes.

@en

The Morro Vermelho Cave (MVC) (Brazil) developed within the Morro Vermelho karst system, which affected Neoproterozoic limestones (Salitre Formation). The MVC experienced little interactions with meteoric processes and is an example of a hypogenic cave formed during strike-slip deformation. The Salitre carbonates in the MVC experienced distributed deformation along an elongated domain overlying a buried strike-slip fault. Gently dipping, semiductile shear zones formed with decimeter-scale (3.9 in.) dolomitic veins. In our model, Mg-rich fluids flowing along the Salitre aquifer caused at the same time extensive dolomitization of the body of rock (100-m [328-ft] scale) experiencing distributed deformation. With progressive displacement, the deep strike-slip fault propagated upward causing the development of an anticline pop-up, steepening sedimentary layers, and steep 1-10-m-long (3.3-33.8-ft) fractures, which served as pathways for upward fluid flow. These steep extensional fractures made it possible for fluids flowing in lower, quartzitic aquifers to enter the carbonate aquifer causing silica deposition in rock cavities and in fractures and fault planes. Following the main stage of speleogenesis, silica deposition took over again depositing on the cave walls a continuous silica crust, rarely observed in other settings worldwide. The interplay between regional bedding-parallel flow and focused circulation of fluids along steep faults and dipping layers, and the associated rock-fluid interactions are not unique to the contractional settings presented but can also occur in association with similar faults in rifted continental margins.

@en

Geological reservoirs can be extensively fractured but the well-test signatures observed in the wells may not show a pressure transient response that is representative of naturally fractured reservoirs (NFRs): for example, one that indicates two distinct pore systems (i.e. the mobile fractures and immobile matrix). Yet, the production behaviour may still be influenced by these fractures. To improve the exploitation of hydrocarbons from NFRs, we therefore need to improve our understanding of fluid-flow behaviour in fractures.

Multiple techniques are used to detect the presence and extent of fractures in a reservoir. Of particular interest to this work is the analysis of well-test data in order to interpret the flow behaviour in an NFR. An important concept for interpreting well-test data from an NFR is the theory of dual-porosity model. However, several studies pointed out that the dual-porosity model may not be appropriate for interpreting well tests from all fractured reservoirs.

This paper therefore uses geological well-testing insights to explore the limitations of the characteristic flow behaviour inherent to the dual-porosity model in interpreting well-test data from Type II and III NFRs of Nelson's classification. To achieve this, we apply a geoengineering workflow with discrete fracture matrix (DFM) modelling techniques and unstructured-grid reservoir simulations to generate synthetic pressure transient data in both idealized fracture geometries and real fracture networks mapped in an outcrop of the Jandaira Formation. We also present key reservoir features that account for the classic V-shape pressure derivative response in NFRs. These include effects of fracture skin, a very tight matrix permeability and wells intersecting a minor, unconnected fracture close to a large fracture or fracture network. Our findings apply to both connected and disconnected fracture networks.@en

In carbonate rocks, channelized fluid flow through fracture conduits can result in the development of large and connected karst networks. These cavity systems have been found in multiple hydrocarbon and geothermal reservoirs, and are often associated with high-permeability zones, but also pose significant challenges in drilling and reservoir management. Here, we expand on the observed interplay between fractures, fluid flow and large cave systems, using outcrop analysis, drone imagery and fluid-flow modelling. The studied carbonate rocks are heavily fractured and are part of the Salitre Formation (750–650 Ma), located in central Bahia (NE Brazil). Firstly, the fracture and cave network data show a similar geometry, and both systems depict three main orientations, namely; NNE–SSW, NW–SE and ESE–WNW. Moreover, the two datasets are dominated by the longer NNE–SSW features. These observed similarities suggest that the fractures and caves are related. The presented numerical results further acknowledge this observed correlation. These results show that open fractures act as the main fluid-flow conduits, with the aperture model defining the fracture-controlled flow contribution. Furthermore, the performed modelling highlights that geometrical features such as length, orientation and connectivity play an important role in the preferred flow orientations.

@en

Regional natural fracture networks often show variations on a scale below that captured by seismic reflection data. This variability is not considered in most reservoir models, but likely impacts uncertainties in permeability. We quantify this uncertainty using a database of 13,000 fractures in nine outcrops digitised in the carbonate Jandaíra Formation (Potiguar basin, Brazil). Distance between outcrops is on average 11 km, with a minimum of 300 m, which is comparable to the distance between wells in naturally fractured reservoirs. In between outcrops, significant variations exist in orientation, intensity, length and topology. Using discrete fracture-matrix flow models, we model the permeability of each deterministic pattern and find that small changes in geometry and topology result in permeability variations that are not captured by connectivity-based analyses such as percolation probabilities, particularly when the matrix is permeable. The permeability variations associated with subseismic-scale fracture variability are not captured in conventional stochastic models, but can be captured using deterministic outcrop models with flow through discrete fractures. The deterministic models provide a permeability range associated with subseismic fracture variability, that can be assigned to grid cells of fractured reservoir flow models, as an alternative to assuming constant permeability in the absence of subseismic-scale deformation.@en

The shallow-marine carbonate rocks of the Jandaíra Formation have been subject to significant permeability variations through time due to various events of fracturing and calcite cementation. As a consequence, the Jandaíra Formation accommodated fluid flow only during specific moments in time. We reconstructed these episodes of fluid flow based on isotope characterizations and microscope characteristics of calcite veins and host rock cements. The Jandaíra Formation, which belongs to the post-rift sequence of the Potiguar Basin in northeast Brazil, was deposited from the Turonian onward until a marine regression exposed it in the Campanian. Due to the subaerial exposure, meteoric waters flushed out marine connate waters, leading to an event of early diagenesis and full cementation of the Jandaíra Formation. Fluid flow through the resulting impermeable carbonate formation appears to be closely related to fracturing. Fracturing in the Late Cretaceous induced a drastic increase in permeability, giving rise to extensive fluid circulation. Host rock dissolution associated to the circulating fluids led to calcite vein cementation within the fracture network, causing it to regain an impermeable and sealing character. In the research area, fluid flow occurred during early burial of the Jandaíra Formation at estimated depths of 400–900 m. This study documents the first application of fluid inclusion isotope analysis on vein precipitates, which allowed full isotopic characterization of the paleo-fluids responsible for calcite vein cementation. The fluid inclusion isotope data indicate that upwelling of groundwater from the underlying Açu sandstones provided the fluids to the fracture network. In Miocene times, renewed tectonic compression of a lower intensity created a secondary fracture network in the Jandaíra Formation. The density of this fracture network, however, was too low to induce a new episode of fluid circulation. As a result, this tectonic event is associated with the development of barren extensional fractures.

@en

We investigated the Cretaceous Potiguar Basin in the Equatorial margin of Brazil to understand how the geometry of major faults evolved to form the basin internal architecture. Previous studies pointed out that the rift is an asymmetrical half-graben elongated along the NE-SW direction. We used 2D seismic, well logs and 3D gravity modeling to analyze faults that constitute the rift boundary and determine their maximum displacement (D_max) and length (L) ratio in the Potiguar Rift. We constrained the 3D gravity modeling with well data and the interpretation of seismic sections. The difference of the fault displacement and depth of the basement obtained in the gravity model is in the order of 10% compared to seismic and well data. The fault-growth curves allowed us to divide the faulted rift border into four main fault systems, which provide roughly similar D_max/L ratios. Fault-growth curves suggest that a regional uniform tectonic mechanism influenced growth of these faults. These fault systems are composed of minor faults that we define as segments. The variation of the displacements along the fault segments indicates that the fault systems were formed independently during rift initiation and were linked by hard and soft linkages. The latter formed relay ramps. In the interconnection zones the D_max/L ratios are highest due to interference of fault segment motions. We divided the evolution of the Potiguar Rift into five stages based on these ratios and correlated them with the major tectonic stages of the breakup between South America and Africa in the Early Cretaceous.

@en