Mauro Cacace

Modelling of fractured carbonate reservoirs: outline of a novel technique via a case study from the Molasse Basin, southern Bavaria, Germany

Fluid flow in low-permeable carbonate rocks depends on the density of fractures, their interconnectivity and on the formation of fault damage zones. The present-day stress field influences the aperture hence the transmissivity of fractures whereas paleostress fields are responsible for the formation of faults and fractures. In low-permeable reservoir rocks, fault zones belong to the major targets. Before drilling, an estimate for reservoir productivity of wells drilled into the damage zone of faults is therefore required. Due to limitations in available data, a characterization of such reservoirs usually relies on the use of numerical techniques. The requirements of these mathematical models encompass a full integration of the actual fault geometry, comprising the dimension of the fault damage zone and of the fault core, and the individual population with properties of fault zones in the hanging and foot wall and the host rock. The paper presents both the technical approach to develop such a model and the property definition of heterogeneous fault zones and host rock with respect to the current stress field. The case study describes a deep geothermal reservoir in the western central Molasse Basin in southern Bavaria, Germany. Results from numerical simulations indicate that the well productivity can be enhanced along compressional fault zones if the interconnectivity of fractures is lateral caused by crossing synthetic and antithetic fractures. The model allows a deeper understanding of production tests and reservoir properties of faulted rocks.

Sensitivity of 3D thermal models to the choice of boundary conditions and thermal properties: a case study for the area of Brandenburg (NE German Basin)

Based on newly available data of both, the structural setting and thermal properties, we compare 3D thermal models for the area of Brandenburg, located in the Northeast German Basin, to assess the sensitivity of our model results. The structural complexity of the basin fill is given by the configuration of the Zechstein salt with salt diapirs and salt pillows. This special configuration is very relevant for the thermal calculations because salt has a distinctly higher thermal conductivity than other sediments. We calculate the temperature using a FEMethod to solve the steady state heat conduction equation in 3D. Based on this approach, we evaluate the sensitivity of the steady-state conductive thermal field with respect to different lithospheric configurations and to the assigned thermal properties. We compare three different thermal models: (a) a crustal-scale model including a homogeneous crust, (b) a new lithosphere-scale model including a differentiated crust and (c) a crustal-scale model with a stepwise variation of measured thermal properties. The comparison with measured temperatures from different structural locations of the basin shows a good fit to the temperature predictions for the first two models, whereas the third model is distinctly colder. This indicates that effective thermal conductivities may be different from values determined by measurements on rock samples. The results suggest that conduction is the main heat transport mechanism in the Brandenburg area.

Influence of fluid flow on the regional thermal field: results from 3D numerical modelling for the area of Brandenburg (North German Basin)

We analyse the effect of fluid flow on the recent thermal field for the Brandenburg region (North German Basin) which is strongly affected by salt structures. The basin fill is modified by a thick layer of mobilized salt (Zechstein, Upper Permian) that decouples the overburden from deeper parts of the lithosphere and is responsible for thermal anomalies since salt has a distinctly higher thermal conductivity than the surrounding sediments and is impermeable to fluid flow. Numerical simulations of coupled fluid flow and heat transfer are carried out to investigate the influence of fluid flow on the shallow temperature field above the Zechstein salt, based on the finite element method. A comparison of results from conductive and coupled modelling reveals that the temperature field down to the low-permeable Triassic Muschelkalk is influenced by fluids, where the shallow low-permeable Tertiary Rupelian-clay is absent. Overall cooling is induced by forced convective forces, the depth range of which is controlled by the communication pathways between the different aquifers. Moreover, buoyancy-induced effects are found in response to temperature-dependent differences in the fluid density where forced convective forces are weak. The range of influence is controlled by the thickness and the permeability of the permeable strata above the Triassic Muschelkalk. With increasing depth, thermal conduction mainly controls the short-wavelength pattern of the temperature distribution, whereas the long-wavelength pattern results from interaction between the highly conductive crust and low-conductive sediments. Our results provide generic implications for basins affected by salt tectonics.

Characterization of main heat transport processes in the Northeast German Basin: Constraints from 3-D numerical models

To investigate and quantify main physical heat driving processes affecting the present-day subsurface thermal field, we study a complex geological setting, the Northeast German Basin (NEGB). The internal geological structure of the NEGB is characterized by the presence of a relatively thick layer of Permian Zechstein salt (up to 5000 m), which forms many salt diapirs and pillows locally reaching nearly the surface. By means of three-dimensional numerical simulations we explore the role of heat conduction, pressure, and density driven groundwater flow as well as fluid viscosity related effects. Our results suggest that the regional temperature distribution within the basin results from interactions between regional pressure forces as driven by topographic gradients and thermal diffusion locally enhanced by thermal conductivity contrasts between the different sedimentary rocks with the highly conductive salt playing a prominent role. In contrast, buoyancy forces triggered by temperature-dependent fluid density variations are demonstrated to affect only locally the internal thermal configuration. Locations, geometry, and wavelengths of convective thermal anomalies are mainly controlled by the permeability field and thickness values of the respective geological layers.

Impact of single inclined faults on the fluid flow and heat transport: results from 3-D finite element simulations

The impact of inclined faults on the hydrothermal field is assessed by adding simplified structural settings to synthetic models. This study is innovative in carrying out numerical simulations because it integrates the real 3-D nature of flow influenced by a fault in a porous medium, thereby providing a useful tool for complex geothermal modelling. The 3-D simulations for the coupled fluid flow and heat transport processes are based on the finite element method. In the model, one geological layer is dissected by a dipping fault. Sensitivity analyses are conducted to quantify the effects of the fault’s transmissivity on the fluid flow and thermal field. Different fault models are compared with a model where no fault is present to evaluate the effect of varying fault transmissivity. The results show that faults have a significant impact on the hydrothermal field. Varying either the fault zone width or the fault permeability will result in relevant differences in the pressure, velocity and temperature field. A linear relationship between fault zone width and fluid velocity is found, indicating that velocities increase with decreasing widths. The faults act as preferential pathways for advective heat transport in case of highly transmissive faults, whereas almost no fluid may be transported through poorly transmissive faults.

MeshIt—a software for three dimensional volumetric meshing of complex faulted reservoirs

Flow, mass and energy transport processes in natural reservoirs are controlled to a large degree by the presence of geological heterogeneities including structures such as fractures and fault zones embedded in a spatially varying three-dimensional (3D) porous matrix of the reservoir. Despite recent advances, currently, state-of-the-art models rely on a number of simplifications partly related to our inability to represent heterogeneities as observed in the field into dynamic model realizations. In this respect, an adequate geometric representation of the discrete system is a basic requirement. In this study, we show how fundamental concept from computational geometry can be assembled and used to bridge the gap between geological and dynamic forward models. The result is an automated, open source software solution (MeshIt) to generate quality 3D meshes suitable for the study of flow and transport processes in faulted and fractured reservoirs. The software enables us to integrate into a 3D volumetric representation dipping structures, comprising fault zones and fractures as well as inclined well paths. This permits us to correctly simulate interactions between discrete flow paths along these interacting components and the 3D flow within the reservoir matrix. The crucial factor that makes the approach applicable to real case reservoirs is that all algorithms are local and scalable parallel and have computing times increasing approximately linearly with data volumes. We test the performance and the robustness of the software against three different scenarios of increasing complexity and further discuss current limitations and range of applicability of the software. Although all examples describe geothermal applications, it is worth mentioning that the approach is equally valid for other applications in geoscience from oil and gas industry to carbon capture and sequestration issues.

Controls on the deep thermal field: implications from 3-D numerical simulations for the geothermal research site Groß Schönebeck

The deep thermal field in sedimentary basins can be affected by convection, conduction or both resulting from the structural inventory, physical properties of geological layers and physical processes taking place therein. For geothermal energy extraction, the controlling factors of the deep thermal field need to be understood to delineate favorable drill sites and exploitation compartments. We use geologically based 3-D finite element simulations to figure out the geologic controls on the thermal field of the geothermal research site Groß Schönebeck located in the E part of the North German Basin. Its target reservoir consists of Permian Rotliegend clastics that compose the lower part of a succession of Late Carboniferous to Cenozoic sediments, subdivided into several aquifers and aquicludes. The sedimentary succession includes a layer of mobilized Upper Permian Zechstein salt which plays a special role for the thermal field due to its high thermal conductivity. Furthermore, the salt is impermeable and due to its rheology decouples the fault systems in the suprasalt units from subsalt layers. Conductive and coupled fluid and heat transport simulations are carried out to assess the relative impact of different heat transfer mechanisms on the temperature distribution. The measured temperatures in 7 wells are used for model validation and show a better fit with models considering fluid and heat transport than with a purely conductive model. Our results suggest that advective and convective heat transport are important heat transfer processes in the suprasalt sediments. In contrast, thermal conduction mainly controls the subsalt layers. With a third simulation, we investigate the influence of a major permeable and of three impermeable faults dissecting the subsalt target reservoir and compare the results to the coupled model where no faults are integrated. The permeable fault may have a local, strong impact on the thermal, pressure and velocity fields whereas the impermeable faults only cause deviations of the pressure field.

Evaluation of three exploitation concepts for a deep geothermal system in the North German Basin

In this paper, we compare three exploitation strategies for the deep geothermal system of Gro?s Schonebeck in the North German Basin. Investigating optimum reservoir designs is one of the key issues for efficient and sustainable utilization of geothermal resource. With this objective we simulate the hydraulic-thermal coupled subsurface processes related to the provision of geothermal energy. The presented application, including visualization, mesh generation and numerical simulation, is based on open source software. The numerical investigations of the three exploitation concepts take into account all geological layers, major natural fault zones, hydraulic fractures, geothermal wells and related hydraulic-thermal coupled processes. In the current exploitation concept, the fluid flows through the rock matrix between the injection and the production well (matrix dominated). The related numerical model is compared and calibrated to available field data. Then, the model is used to investigate two alternative stimulation concepts. All three concepts were evaluated taking into account the evolution of the production temperature as well as the hydraulic conductivity between production and injection well. As an alternative to the current situation, a fracture dominated system is investigated where the fluid flows through hydraulically induced fractures between injection and production well. Compared to the reference model, a twofold increase in productivity could be observed together with a significantly reduced time before the onset of a thermal breakthrough. The second alternative is a hybrid concept combining both matrix and fracture-dominated flow paths between the production and the injection well. We show that this hybrid approach could significantly increase the reservoir productivity and prolongs the time before the onset of thermal breakthrough. HighlightsBuilding numerical forward simulations based on geological models.Simulation of hydraulic-thermal coupled processes of faulted and fractured reservoirs.Evaluation of exploration concepts for deep geothermal systems.

Coupled thermo-mechanical 3D subsidence analysis along the SW African passive continental margin

Sedimentary basins along the SW African margin deliver information about processes, which occurred in the past and serve as a good starting point for reconstructing the margin’s evolution. By integrating detailed information on the present-day configuration of the SW African margin into a 3D thermo-mechanical forward modeling framework, we attempt a margin-wide reconstruction of its subsidence history since its onset during breakup. The 3D forward modeling approach as applied on the SW African margin area makes use of the coupling between thermal relaxation of the lithosphere and flexural isostatic balance in response to sediment deposition and build-up of thermal stresses during lithosphere cooling. On the one hand, our results provide useful information about the behavior of the lithosphere during breakup and the subsequent post-rift phase. On the other hand, restored paleobathymetries for specific time intervals during the post-rift evolution give evidence for possible uplift events superposed on the long-term subsidence history of the margin. Restored paleobathymetries are in agreement with conclusions derived from former studies and provided strong indications for the presence of a rather heterogeneous crustal configuration of the margin marked by a mechanical decoupling of the upper and lower crustal domains during most of the rifting phase.

Evaluating Micro-Seismic Events Triggered by Reservoir Operations at the Geothermal Site of Groß Schönebeck (Germany)

This study aims at evaluating the spatial and temporal distribution of 26 micro-seismic events which were triggered by hydraulic stimulation at the geothermal site of Groß Schönebeck (Germany). For this purpose, the alteration of the in-situ stress state and the related change of slip tendency for existing fault zones due to stimulation treatments and reservoir operations is numerical simulated. Changes in slip tendency can potentially lead to reactivation of fault zones, the related movement can lead to the occurrence of seismic events. In the current numerical study, results obtained based on the thermal–hydraulic–mechanical coupled simulation are compared to field observations. In particular, the study focuses on describing the fault reactivation potential: (1) under in-situ stress conditions; (2) during a waterfrac stimulation treatment; and (3) during a projected 30 years production and injection period at the in-situ geothermal test-site Groß Schönebeck. The in-situ stress state indicates no potential for fault reactivation. During a waterfrac stimulation treatment, micro-seismic events were recorded. Our current evaluation shows an increase of slip tendency during the treatment above the failure level in the direct vicinity of the micro-seismic events. During the projected production and injection period, despite increased thermal stress, the values for slip tendency are below the threshold for fault reactivation. Based on these results, and to prove the applied method to evaluate the observed micro-seismic events, a final discussion is opened. This includes the in-situ stress state, the role of pre-existing fault zones, the adopted criterion for fault reactivation, and a 3D rock failure criterion based on true triaxial measurements.